[发明专利]一种基于各向异性插值的地层分层方法

权利要求书

1.一种基于各向异性插值的地层分层方法，其特征在于：

S1：在工程区进行地质勘探，采集各勘探钻孔的岩性界面的信息数据，所述岩性界面的信息数据包含地层岩性名称及地层岩性深度范围；

S2：导入各所述钻孔的岩性界面的信息数据，在整个勘探空间构建大立方网，该大立方网包含多个小立方网格，参考小立方网格大小把钻孔沿深度方向离散得到多个离散数据点；

S3:选择一地层岩性名称，在各所述钻孔内，将落在该地层岩性深度范围内的离散数据点的属性值赋值为1，其余数据点的属性值赋值为0，最后将这些初始化的离散数据点作为样本点；

S4：根据所述样本点的属性值，采用各向异性Kriging插值方法，得到每个小立方网格位置中心点的属性值，每个属性值的计算包含与样本数据的关联程度，所述关联程度包括水平方向和垂直方向的权重关系，周围样本位置分布；遍历所有的小立方网格位置中心点的属性值，将属性值大于预设中间值的小立方网格位置中心点的属性值改为1，小于这一中间值的小立方网格位置中心点的属性值改为0，并把所述样本点的属性值赋给距离样本点最近的小立方网格的位置中心点；若各小立方网格位置中心点的属性值为1，则认为该小立方网格位置中心点属于该岩性界面对应的地层，若各小立方网格位置中心点的属性值为0，则认为该小立方网格位置中心点不属于该岩性界面对应的地层；

S5：在大立方网中，把所有属性值为1的小立方网格的位置中心点连接起来生成等值面；

S6：重复步骤S3至S5，依次创建所述钻孔内各地层岩性的等值面；

S7：依据各地层岩性的等值面，完成所述勘探区的地层分层。

2.根据权利要求1所述的基于各向异性插值的地层分层方法，其特征在于：所述预设中间值为0.5。

3.根据权利要求1所述的基于各向异性插值的地层分层方法，其特征在于，所述步骤S4还包括：以地层岩性界面的产状面为参考面插值计算得到各所述小立方网格位置中心点的属性值；当岩性界面没有产状面时，以地表趋势面为参考面插值计算得到各所述小立方网格位置中心点的属性值。

4.根据权利要求1所述的基于各向异性插值的地层分层方法，其特征在于：所述地层岩性名称包括素填土、淤泥、粉质粘土、粉土及粉砂、砂岩、泥岩。

5.根据权利要求1所述的基于各向异性插值的地层分层方法，其特征在于：将小立方网格位置中心点连成生成等值面，所述等值面内的空间形成一个等值层。

6.根据权利要求1所述的基于各向异性插值的地层分层方法，其特征在于，所述步骤S5之后还包括：利用光滑的离散插值方法将所述等值面进行优化处理得到光滑的等值面。

7.根据权利要求1所述的基于各向异性插值的地层分层方法，其特征在于，所述地层分层方法还包括：将岩性界面的厚度为零的钻孔位置加入至等值面中，来确定地层的尖灭位置。