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被遍历完，完成所有钻孔虚拟地层的推断，更新地层信息，保存含有虚拟地层的钻孔数据；

其中，步骤(5)具体包括：

(5‑1)从推断钻孔集合ND中获取任一推断钻孔dk，判断推断钻孔dk是否存在虚拟地层；

(5‑2)若存在虚拟地层，则根据钻孔集合UD内每个钻孔的地层厚度，采用反距离加权插值法计算虚拟地层的厚度；

(5‑3)循环执行步骤(5‑1)‑(5‑2)，直至完成推断钻孔集合ND中所有推断钻孔的遍历，得到所有推断钻孔的虚拟地层推断；

步骤(5‑1)具体包括：

(5‑1‑1)从推断钻孔集合ND中获取任一推断钻孔dk，根据钻孔dk与三角网M内三角形的位置推断是否存在虚拟地层，具体方法如下：(a)读取三角网M内任一三角形mt，根据三角形mt三个顶点与钻孔dk点坐标采用同向法判断钻孔dk是否位于三角形mt内；(b)若钻孔dk位于三角形mt内，则获取mt三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为

0，则表示推断钻孔dk存在虚拟地层，执行步骤(5‑2)，反之表示推断钻孔dk不存在虚拟地层，则设置推断钻孔dk对应的虚拟地层厚度为0；若钻孔dk不位于三角形mt内，则返回执行步骤(a)；(c)若三角网M被遍历完，钻孔dk不位于所有三角形mt内，则执行步骤(5‑1‑2)；

(5‑1‑2)根据推断钻孔dk类型推断虚拟地层，具体方法如下：(a)若dk位于非边界钻孔集合FBD内，则表示dk不存在虚拟地层，则设置dk对应的虚拟地层厚度为0；(b)若dk位于边界钻孔集合BD内，则计算三角网M内每一三角形的中心与dk点的距离，根据距离大小获取离dk最近的三角形mp；(c)以钻孔dk点坐标为中心，dk与三角形mp中心的距离为半径构建几何圆；若几何圆内存在非边界钻孔，则表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0；若几何圆内不存在非边界钻孔，则获取mp三个顶点对应的钻孔dr、dr+1与dr+2，若dr、dr+1与dr+2的地层厚度都不为0，则表示dk存在虚拟地层，执行步骤(5‑2)，反之表示dk不存在虚拟地层，设置dk对应的虚拟地层厚度为0。

2.根据权利要求1所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(1)具体包括：(1‑1)读取钻孔数据到钻孔集合D＝{di|i＝1,2,…,DN}；其中，i表示钻孔序号，DN表示钻孔数量；

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

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

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

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＝{mt|t＝1,2,…,TM}，其中，mt表示三角网内第t个三角形，TM表示三角网内三角形数量。

6.根据权利要求5所述的面向地质钻孔的虚拟地层推断方法，其特征在于：步骤(4‑1)具体包括：(4‑1‑1)根据非推断钻孔集合ID中每一钻孔坐标信息，构建Delaunary三角网T0，并获取T0内的所有三角边，存入三角边集合E；

(4‑1‑2)从三角边集合E中读取任一三角边e，获取位于三角边e两端的钻孔dl、dl+1，在两钻孔之间的1/2处创建一尖灭钻孔do，若钻孔dl、dl+1的地层厚度slc.h、s(l+1)c.h均不为0，则将尖灭钻孔do的地层厚度soc.h赋值为(slc.h+s(l+1)c.h)/2，否则soc.h＝0，(4‑1‑3)重复执行步骤(4‑1‑2)，直至实现集合E中所有三角边的遍历，得到所有尖灭钻孔，并将所有尖灭钻孔存入尖灭钻孔集合AD中。

7.一种面向地质钻孔的虚拟地层推断装置，包括处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于：所述处理器执行所述程序时实现权利要求1‑6中任意一项所述的方法。