[发明专利]一种基于三维倾角属性体校正的古地貌分析方法

权利要求书

1.一种基于三维倾角属性体校正的古地貌分析方法，其特征在于采用以下步骤：

步骤1：地层残余视厚度计算，具体如下，

步骤1-1：根据研究区地层沉积与层序研究成果，选取目的层下方在地震剖面上表现为强反射特征的、全区稳定分布的、可连续追踪的且具有地质意义的基准面，得地层残余视厚度；

步骤1-2：基于根据研究区井震标定结果得到的时深关系进行速度场建立，将步骤1-1得到的时间域的目的层位及基准面转换到深度域；

步骤1-3：将步骤1-2得到的深度域的目的层位与基准面相减，并剔除结果中异常值，得到目的层位与基准面之间地层的残余视厚度数据H(x,y)；

步骤2：基于梯度结构张量的地层倾角数据提取，具体如下，

步骤2-1：采用公式(1)构建梯度结构张量矩阵，记为GST(x,y,t)，并根据矩阵特征分解，进行矩阵分解：

式中：u(x,y,t)为地震数据；g₁，g₂，g₃分别是三维地震数据在x，y，t方向上的梯度向量体；Δx，Δy，Δt三维地震数据在x，y，t方向上的间隔；λ₁，λ₂，λ₃分别是梯度结构张量矩阵的三个特征值；v₁，v₂，v₃分别是特征值λ₁，λ₂，λ₃对应的特征向量；

步骤2-2：利用步骤2-1得到的梯度结构张量的最大特征值对应的特征向量求取局部数据的视倾角，采用公式(2)进行计算：

式中：p(x,y,t)，q(x,y,t)分别为三维地震数据x和y方向的视倾角；v_1x，v_1y，v_1t分别为梯度结构张量矩阵的最大特征值λ₁对应的特征向量v₁的三个分量；

步骤3：构建地层厚度校正系数，计算地层残余真实厚度，提高古地貌分析精度，具体如下，

步骤3-1：计算平面上每一点目的层与其下方基准面之间地层的x和y方向的平均视倾角p_mean(x,y)，q_mean(x,y)，采用公式(3)计算：

式中：t₀为目的层地震时间层位数据；t_e为基准面地震时间层位数据；dt为地震数据采样率；

步骤3-2：采用平面上每一点x、y方向的平均视倾角p_mean(x,y)，q_mean(x,y)数据，采用公式(4)构建平面上每一点地层厚度的校正系数，计算地层的真实厚度：

式中：K(x,y)为平面上每一点地层厚度的校正系数；h(x,y)利用地震倾角属性数据校正得到的地层真实厚度；

步骤3-3：对利用地层倾角数据校正后的地层厚度h(x,y)进行平滑处理，并将处理结果进行平面成图，表征古地貌的形态。