[发明专利]三维地质覆盖层界面建模方法、装置、设备及存储介质

权利要求书

1.一种三维地质覆盖层界面建模方法，其特征在于，包括：

步骤1：根据覆盖层的边界线、该边界线内的至少一个已知勘探点创建一平面三角网格，所述平面三角网格的角点包括边界线节点、已知勘探点以及多个新增节点；

步骤2：设定计算覆盖层深度的拟合函数；

步骤3：当已知勘探点为一个时，根据已知勘探点的覆盖层深度、已知勘探点与距离该已知勘探点最近的边界线节点之间的水平距离，求出所述拟合函数的常量；当已知勘探点为两个或两个以上时，根据已知勘探点的覆盖层深度、已知勘探点与距离该已知勘探点最近的边界线节点之间的水平距离、以及已知勘探点与待求新增节点之间的水平距离，求出所述拟合函数的常量或动态参数；

步骤4：根据所述步骤3中的拟合函数计算出该待求新增节点的覆盖层深度；

步骤5：重复步骤3和4，计算出每个新增节点的覆盖层深度；

步骤6：将所述平面三角网格的所有角点投影到地表，再分别下降各角点对应的覆盖层深度，即将边界线节点、已知勘探点以及多个新增节点三角化，生成覆盖层界面的三角网模型。

2.如权利要求1所述的一种三维地质覆盖层界面建模方法，其特征在于，所述拟合函数为幂函数，所述幂函数的数学表达式为：

f(d)＝φ·d^a＝h

其中，φ为动态参数，0＜a＜1，a为常量，d为已知勘探点与距离该已知勘探点最近的边界线节点之间的水平距离，或者d为新增节点与距离该新增节点最近的边界线节点之间的水平距离，h为已知勘探点的覆盖层深度或新增节点的覆盖层深度，且：

当已知勘探点为一个或两个时，动态参数φ等于1，常量a根据所述步骤3求出；当已知勘探点为三个及三个以上时，常量a根据经验设置，动态参数φ根据所述步骤3求出。

3.如权利要求2所述的一种三维地质覆盖层界面建模方法，其特征在于，当已知勘探点为三个及三个以上时，所述常量a为0.6。

4.如权利要求2所述的一种三维地质覆盖层界面建模方法，其特征在于，所述步骤3中，当已知勘探点为两个时，常量a采用如下公式来求得：

其中，A₁为根据一个已知勘探点K₁计算得到的常量，A₂为根据另一个已知勘探点K₂计算得到的常量，s₁为待求新增节点到一个已知勘探点K₁的水平距离，s₂为待求新增节点到另一个已知勘探点K₂的水平距离。

5.如权利要求2或3所述的一种三维地质覆盖层界面建模方法，其特征在于，所述步骤3中，当已知勘探点为三个或三个以上时，动态参数φ的计算步骤为：

步骤3.1：根据每个已知勘探点计算出一组动态参数值φ₁,φ₂,…,φ_i,…,φ_n，φ_i为根据第i个已知勘探点得到的动态参数，n为已知勘探点的个数，n≥3；

步骤3.2：找出待求新增节点与n个已知勘探点之间水平距离的最大值，计算水平距离的最大值与每个水平距离的比值，所述每个水平距离是指每个已知勘探点与待求新增节点之间的水平距离，即：

s_max/s_i＝t_i

其中，s_max为待求新增节点与n个已知勘探点之间水平距离的最大值，s_i为第i个已知勘探点与待求新增节点之间的水平距离，t_i为第i个已知勘探点对应的比值；

步骤3.3：计算每个已知勘探点对待求新增节点的影响程度，即权重，具体计算公式如下：

其中，ε_i为第i个已知勘探点对待求新增节点的权重；

步骤3.4：动态参数φ的计算公式为：

φ＝φ₁ε₁+φ₂ε₂+…+φ_iε_i+…+φ_nε_n。