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

有效

• 申請(qǐng)?zhí)枺?/td>: 202010087886.3
• 申請(qǐng)日: 2020-02-12
• 公開（公告）號(hào): CN111239818B
• 公開（公告）日: 2022-02-25
• 發(fā)明（設(shè)計(jì)）人: 熊曉軍;陳容;袁野
• 申請(qǐng)（專利權(quán)）人: 成都理工大學(xué)
• 主分類號(hào): G01V1/30
• 分類號(hào): G01V1/30
• 代理公司: 暫無(wú)信息
• 代理人: 暫無(wú)信息

1.一種基于三維傾角屬性體校正的古地貌分析方法，其特征在于采用以下步驟：

步驟1：地層殘余視厚度計(jì)算，具體如下，

步驟1-1：根據(jù)研究區(qū)地層沉積與層序研究成果，選取目的層下方在地震剖面上表現(xiàn)為強(qiáng)反射特征的、全區(qū)穩(wěn)定分布的、可連續(xù)追蹤的且具有地質(zhì)意義的基準(zhǔn)面，得地層殘余視厚度；

步驟1-2：基于根據(jù)研究區(qū)井震標(biāo)定結(jié)果得到的時(shí)深關(guān)系進(jìn)行速度場(chǎng)建立，將步驟1-1得到的時(shí)間域的目的層位及基準(zhǔn)面轉(zhuǎn)換到深度域；

步驟1-3：將步驟1-2得到的深度域的目的層位與基準(zhǔn)面相減，并剔除結(jié)果中異常值，得到目的層位與基準(zhǔn)面之間地層的殘余視厚度數(shù)據(jù)H(x,y)；

步驟2：基于梯度結(jié)構(gòu)張量的地層傾角數(shù)據(jù)提取，具體如下，

步驟2-1：采用公式(1)構(gòu)建梯度結(jié)構(gòu)張量矩陣，記為GST(x,y,t)，并根據(jù)矩陣特征分解，進(jìn)行矩陣分解：

式中：u(x,y,t)為地震數(shù)據(jù)；g₁，g₂，g₃分別是三維地震數(shù)據(jù)在x，y，t方向上的梯度向量體；Δx，Δy，Δt三維地震數(shù)據(jù)在x，y，t方向上的間隔；λ₁，λ₂，λ₃分別是梯度結(jié)構(gòu)張量矩陣的三個(gè)特征值；v₁，v₂，v₃分別是特征值λ₁，λ₂，λ₃對(duì)應(yīng)的特征向量；

步驟2-2：利用步驟2-1得到的梯度結(jié)構(gòu)張量的最大特征值對(duì)應(yīng)的特征向量求取局部數(shù)據(jù)的視傾角，采用公式(2)進(jìn)行計(jì)算：

式中：p(x,y,t)，q(x,y,t)分別為三維地震數(shù)據(jù)x和y方向的視傾角；v_1x，v_1y，v_1t分別為梯度結(jié)構(gòu)張量矩陣的最大特征值λ₁對(duì)應(yīng)的特征向量v₁的三個(gè)分量；

步驟3：構(gòu)建地層厚度校正系數(shù)，計(jì)算地層殘余真實(shí)厚度，提高古地貌分析精度，具體如下，

步驟3-1：計(jì)算平面上每一點(diǎn)目的層與其下方基準(zhǔn)面之間地層的x和y方向的平均視傾角p_mean(x,y)，q_mean(x,y)，采用公式(3)計(jì)算：

式中：t₀為目的層地震時(shí)間層位數(shù)據(jù)；t_e為基準(zhǔn)面地震時(shí)間層位數(shù)據(jù)；dt為地震數(shù)據(jù)采樣率；

步驟3-2：采用平面上每一點(diǎn)x、y方向的平均視傾角p_mean(x,y)，q_mean(x,y)數(shù)據(jù)，采用公式(4)構(gòu)建平面上每一點(diǎn)地層厚度的校正系數(shù)，計(jì)算地層的真實(shí)厚度：

式中：K(x,y)為平面上每一點(diǎn)地層厚度的校正系數(shù)；h(x,y)利用地震傾角屬性數(shù)據(jù)校正得到的地層真實(shí)厚度；

步驟3-3：對(duì)利用地層傾角數(shù)據(jù)校正后的地層厚度h(x,y)進(jìn)行平滑處理，并將處理結(jié)果進(jìn)行平面成圖，表征古地貌的形態(tài)。