[發(fā)明專利]一種疊后地震數(shù)據(jù)繞射波屬性提取方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201810896115.1
• 申請(qǐng)日: 2018-08-08
• 公開（公告）號(hào): CN108594303B
• 公開（公告）日: 2019-12-13
• 發(fā)明（設(shè)計(jì)）人: 舒夢(mèng)珵;康洪全;程濤;賈懷存
• 申請(qǐng)（專利權(quán)）人: 中國海洋石油集團(tuán)有限公司;中海油研究總院有限責(zé)任公司
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28
• 代理公司: 11245 北京紀(jì)凱知識(shí)產(chǎn)權(quán)代理有限公司
• 代理人: 徐寧;孫楠

1.一種疊后地震數(shù)據(jù)繞射波屬性提取方法，其特征在于，包括以下步驟：

1)地震數(shù)據(jù)預(yù)處理；

2)根據(jù)實(shí)際數(shù)據(jù)建立正演模型，通過計(jì)算，確定主成分分析技術(shù)提取繞射波時(shí)，最適參數(shù)及主成分的維度；

3)將本征矢量計(jì)算公式改進(jìn)為適用于三維地震數(shù)據(jù)計(jì)算式，并應(yīng)用步驟2)中的最適參數(shù)及主成分維度，提取實(shí)際地震數(shù)據(jù)中的繞射波場(chǎng)；

根據(jù)三維地震數(shù)據(jù)特點(diǎn)，通過下式提取空間繞射波場(chǎng)數(shù)據(jù)：

式中,是計(jì)算矩陣中位置為i,j,k的采樣點(diǎn)上第r個(gè)正交分量；s_i+n,j+m,k+l是計(jì)算矩陣中位置為i,j,k的地震數(shù)據(jù)采樣點(diǎn)值；n,m,l是i,j,k方向上的延遲距；N,M,L為最大延遲距；s_avr為計(jì)算矩陣中樣點(diǎn)振幅的平均值；s_max,s_min分別為計(jì)算矩陣中樣點(diǎn)振幅的最大值和最小值；分別為矩陣特征值和特征向量。

2.如權(quán)利要求1所述方法，其特征在于：所述步驟1)中，預(yù)處理包括以下步驟：

1.1)定義如下變量：A＝{A₁,A₂,...,A_n}是一個(gè)m×n的矩陣，n為輸入的地震道數(shù)，每個(gè)地震道為m維向量，代表m個(gè)采樣點(diǎn)；A_n表示第n個(gè)地震道的地震數(shù)據(jù)采集序列；進(jìn)行地震數(shù)據(jù)體的標(biāo)準(zhǔn)化預(yù)處理：

式中，x_pq為標(biāo)準(zhǔn)化后的地震數(shù)據(jù)采集序列；為樣本均值：A_pq為原始采集的地震數(shù)據(jù)采集序列；S_p為樣本標(biāo)準(zhǔn)差：p＝1,…,n；q＝1,…,m；

1.2)處理后得到標(biāo)準(zhǔn)化地震數(shù)據(jù)體為X＝{x₁,x₂,...,x_n}；其中，X是一個(gè)m×n的矩陣，n為輸入的地震道數(shù)，每個(gè)地震道為m維向量，代表m個(gè)采樣點(diǎn)，x_n表示第n個(gè)地震道標(biāo)準(zhǔn)化后的地震數(shù)據(jù)采集序列。

3.如權(quán)利要求1所述方法，其特征在于：所述步驟2)中，最適參數(shù)及主成分的維度確定方法包括以下步驟：

2.1)通過下式將輸入數(shù)據(jù)矢量x_g變換成新的矢量s_g：

s_g＝U^Tx_g

其中，U是一個(gè)n×n正交矩陣，它的第g列U_g是樣本協(xié)方差矩陣；g＝1，…，n；

2.2)計(jì)算第g個(gè)本征矢量：

λ_gu_g＝Cu_g，g＝1,...,n (1)

其中，C為一個(gè)n階方陣；λ_g是C的一個(gè)本征值，u_g是對(duì)應(yīng)的本征矢量；

2.3)通過上述步驟對(duì)地震模型中繞射波場(chǎng)進(jìn)行提取，當(dāng)達(dá)到壓制反射波同相軸同時(shí)充分保留小尺度異常體繞射波信息的效果時(shí)，確定最適參數(shù)及主成分維度。