[發明專利]一種基于三維傾角屬性體校正的古地貌分析方法

有效

• 申請號: 202010087886.3
• 申請日: 2020-02-12
• 公開（公告）號: CN111239818B
• 公開（公告）日: 2022-02-25
• 發明（設計）人: 熊曉軍;陳容;袁野
• 申請（專利權）人: 成都理工大學
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 暫無信息
• 代理人: 暫無信息

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)進行平滑處理，并將處理結果進行平面成圖，表征古地貌的形態。