[发明专利]地层厚度的确定方法

权利要求书

1.一种地层厚度的确定方法，其特征在于，所述确定方法包括以下步骤：

获取地层的原始数据点，将所述原始数据点进行平面拟合得到地层拟合平面；

通过插值拟合将地层拟合平面进行网格化得到具有网格的地层网格拟合平面；

通过所述原始数据点与所述地层网格拟合平面确定所述原始数据点对应的第一插值在地层深度方向上的绝对差值，对于所述绝对差值处于置信区间以外的原始数据点进行处理并重新进行插值计算；

通过双三次拟合差分方法，确定所述地层网格拟合平面上的并位于网格线交点处的点的切面的法向量；

根据所述法向量获得地层的顶界面与底界面之间的时间域值；

根据所述时间域值与弹性波在地层内传播的速度确定地层的真实厚度。

2.根据权利要求1所述的确定方法，其特征在于，所述地层网格拟合平面通过Voronoi三角网插值获得。

3.根据权利要求1所述的确定方法，其特征在于，所述地层网格拟合平面上的所述网格线交点的插值通过以下公式获得：

其中，

a＝y₂-y₁；

b＝z₂(y₁-y₀)-z₁(y₂-y₀)；

c＝x₁-x₂；

d＝z₁(x₂-x₁)-z₂(x₂-x₀)；

e＝(x₁-x₀)(y₂-y₀)-(x₂-x₀)(y₁-y₀)；

其中，x₁、y₁、z₁为原始数据点P₁的坐标，x₂、y₂、z₂为原始数据点P₂的坐标，x₀、y₀、z₀为待插值点N的坐标，z₀为待插值的未知高程值。

4.根据权利要求1所述的确定方法，其特征在于，所述确定方法包括：

采用距离加权平均算法对不在所述网格线交点处的所述原始数据点进行第二插值计算。

5.根据权利要求4所述的确定方法，其特征在于，所述第二插值通过以下公式获得：

其中，

所述V为不在所述网格线交点处的所述原始数据点至在所述网格线交点处的所述原始数据点的距离；

所述V_i为所述网格线交点处的所述原始数据点的第i个点；

d_i为不在所述网格线交点处的所述原始数据点至在所述网格线交点处的第i个所述原始数据点的第i个距离。

6.根据权利要求1所述的确定方法，其特征在于，所述置信区间为σ，其中，σ＝95％。

7.根据权利要求1所述的确定方法，其特征在于，在确定所述真实厚度之前还包括，将位于网格线交点处的所述原始数据点进行坐标变换，对经坐标变换后的所述原始数据点做法线，所述法线与地层的底界面相交，根据所述法线与地层的底界面的交点确定地层的顶界面与底界面之间的所述真实厚度，其中，所述法线的法向量的确定方法包括：

获取所述地层网格拟合平面的地层界面三维空间数字模型后，通过双三次拟合差分方法的离散法线沿x轴的自左而右的方向，以及沿y轴的自上而下的方向求取三维空间上的顶界面上的所述原始数据点的所述法向量。