[发明专利]面向地质钻孔的虚拟地层推断方法及装置

权利要求书

1.一种面向地质钻孔的虚拟地层推断方法，其特征在于该方法包括：

(1)读取钻孔数据、钻孔地层信息和约束边界矢量数据，生成钻孔集合D、地层集合S、约束边界集合Boundary和地层厚度集合H；

(2)基于约束边界集合Boundary，将钻孔集合D划分为边界钻孔集合BD和非边界钻孔集合FBD，并根据地层集合S确定需推断地层的地层编号区间C；

(3)从地层编号区间C中获取任一地层编号c，构建推断钻孔集合ND与非推断钻孔集合ID；

(4)基于钻孔集合ID，构建尖灭钻孔集合AD，合并集合ID和集合AD为钻孔集合UD，并根据集合UD创建三角网M；

(5)基于三角网M、边界钻孔集合BD和非边界钻孔集合FBD，获取推断钻孔集合ND中所有钻孔缺失的虚拟底层，并采用反距离加权插值计算虚拟地层的厚度；

(6)循环执行步骤(3)-(5)，直至地层编号区间C被遍历完，完成所有钻孔虚拟地层的推断，更新地层信息，保存含有虚拟地层的钻孔数据。

2.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(1)具体包括：

(1-1)读取钻孔数据到钻孔集合D＝{d_i|i＝1，2，…,DN}；其中，i表示钻孔序号，DN表示钻孔数量；

(1-2)读取所有钻孔的地层信息，存入地层集合S＝{s_ij|i＝1,2,…,DN,j＝1,2,…,SN_i}，其中，s_ij表示钻孔d_i的第j个地层，SN_i表示钻孔d_i的地层数量；

(1-3)读取约束边界矢量数据到约束边界集合Boundary；

(1-4)计算每一地层的厚度，存入地层厚度集合H＝{s_ij.h|i＝1,2,…,DN,j＝1,2,…,SN_i},其中，s_ij.h＝s_ij.et-s_ij.eb表示s_ij的地层厚度，s_ij.et为s_ij的层顶深度，s_ij.eb为s_ij的层底深度。

3.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(2)具体包括：

(2-1)获取钻孔集合D内每一钻孔的坐标，判断钻孔是否位于约束边界集合Boundary中的约束边界上，若是则标记钻孔为边界钻孔，并存入边界钻孔集合BD中，否则标记为非边界钻孔，并存入非边界钻孔集合FBD中；

(2-2)从地层集合S中获取钻孔集合D内所有钻孔的末尾地层的地层编号，得到地层编号的最小值α与最大值β，构建需推断地层的地层编号区间C＝[α+1，β],区间内取正整数。

4.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(3)具体包括：

(3-1)从地层编号区间C中获取任一地层编号c；

(3-2)从地层集合S中获取钻孔集合D内每一钻孔的末尾地层的地层编号，若该地层编号小于c，则代表该地层编号对应的钻孔为需推断钻孔，将其保存到需推断钻孔集合ND，否则代表该地层编号对应的钻孔为非推断钻孔，将其保存到非推断钻孔集合ID中；

(3-3)若非推断钻孔集合ID中的非推断钻孔数量LN≥3，则返回执行步骤(3-1)，否则执行步骤(4)。

5.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(4)具体包括：

(4-1)根据非推断钻孔集合ID，采用1/2尖灭规则构建尖灭钻孔集合AD；

(4-2)合并非推断钻孔集合ID与尖灭钻孔集合AD为钻孔集合UD；

(4-3)根据钻孔集合UD内的钻孔坐标信息，创建Delaunary三角网M＝{m_t|t＝1,2,…,TM}，其中，m_t表示三角网内第t个三角形，TM表示三角网内三角形数量。