[发明专利]地层曲面构建方法和装置、设备及储存介质

权利要求书

1.一种地层曲面构建方法，其特征在于，包括：

获取待构建的地层曲面的钻孔数据，并基于获取到的所述钻孔数据构建不规则三角网；其中，所述钻孔数据的个数为多个；

基于所述不规则三角网，对多个所述钻孔数据进行综合预处理，得到预处理结果；其中，所述预处理结果包括划分为n级的钻孔数据；

对所述预处理结果中的钻孔数据进行抽稀处理，得到抽稀处理后的钻孔数据；其中，所述抽稀处理后的钻孔数据参与构建所述地层曲面；

利用曲面样条函数对所述抽稀处理后的钻孔数据进行拟合，构建得到所述地层曲面。

2.根据权利要求1所述的方法，其特征在于，所述基于获取到的所述钻孔数据构建不规则三角网，包括：

将获取到的所述钻孔数据的坐标作为离散点，在多个所述离散点中选取任意一点作为第一点，并由多个所述离散点中查找距所述第一点最近的离散点作为第二点，将所述第一点与所述第二点之间的连线作为初始基线；

在所述初始基线右侧运用Delaunay法则，由多个所述离散点中选取第三点，并依次连接所述第一点、所述第二点和所述第三点，生成相应的三角形；

由生成的所述三角形中，选取除所述初始基线之外的另外两边分别作为重新选定的初始基线，生成相应的三角形，直至所重新选定的初始基线的右侧无所述离散点为止；

其中，所生成的各三角形构建成所述不规则三角网。

3.根据权利要求1所述的方法，其特征在于，基于所述不规则三角网，对多个所述钻孔数据进行综合预处理，包括：

获取所述不规则三角网中各三角形的三个边长；

在当前获取的三角形的三个边长均小于预设阈值时，将当前获取的三角形的三个顶点综合为一个点，并将综合后得到的点作为综合后的钻孔数据；

其中，所述预设阈值表征待构建的所述地层曲面在每一级别下的分辨率；

所述地层曲面的级别通过谷歌瓦片地图算法进行切分得到，且所述地层曲面的级别为n级，n的取值范围为：0＜n≤22。

4.根据权利要求3所述的方法，其特征在于，对所述预处理结果中的钻孔数据进行抽稀处理，包括：

利用谷歌地图的应用程序编程接口将划分为n级后的所述钻孔数据的经纬度坐标换算为谷歌像素索引；

利用所述谷歌像素索引，计算所述第m级下钻孔数据的左上角坐标(x1，y1)和右下角坐标(x2，y2)；

根据所述左上角坐标(x1，y1)和所述右下角坐标(x2，y2)生成第n级下钻孔数据的p行、q列的抽稀格网；

遍历每个所述抽稀格网内的钻孔数据，在当前获取的抽稀格网内存在多个所述钻孔数据时，保留距当前获取的抽稀格网中心距离最近的钻孔数据，得到抽稀处理后的钻孔数据；

其中，m的取值范围为：0≤m≤n-4。

5.根据权利要求1至4任一项所述的方法，其特征在于，利用曲面样条函数对所述抽稀处理后的钻孔数据进行拟合，构建得到所述地层曲面，包括：

读取根据平面弹性理论推导出的位移公式：

其中，A_i表征所述地层曲面i点上设置的垂直位移；Z_A(r)表征在所述地层曲面i点上设置垂直位移A_i，并限制距i点半径为R以外位移为0时，由A_i引起距i点为r处的位移；i为参与构建所述地层曲面的钻孔数据中的任一钻孔数据；

其中，r＝ri，x，y为自变量；

利用所述位移公式，对所述抽稀处理后的钻孔数据进行叠加，得到所述曲面样条函数：

其中，N为所述抽稀处理后的钻孔数据的个数；

读取所建立的实测高程方程：

其中，h_i为实测高程，

基于抽稀处理后的钻孔数据，对所述实测高程方程进行求解，得到所述曲面样条函数中的系数；

利用确定系数后的曲线样条函数，对抽稀处理后的钻孔数据进行拟合，构建得到所述地层曲面。