[發(fā)明專利]一種基于地震屬性的裂縫定量化預(yù)測(cè)的方法

有效

• 申請(qǐng)?zhí)枺?/td>: 201611177648.1
• 申請(qǐng)日: 2016-12-19
• 公開（公告）號(hào): CN106842299B
• 公開（公告）日: 2019-07-09
• 發(fā)明（設(shè)計(jì)）人: 馬曉宇;張研;徐光成;李文科;于豪
• 申請(qǐng)（專利權(quán)）人: 中國(guó)石油天然氣股份有限公司
• 主分類號(hào): G01V1/28
• 分類號(hào): G01V1/28;G01V1/30
• 代理公司: 北京三友知識(shí)產(chǎn)權(quán)代理有限公司 11127
• 代理人: 姚亮;沈金輝

1.一種基于地震屬性的裂縫定量化預(yù)測(cè)的方法，其包括：

獲取目的層的三維地震數(shù)據(jù)和測(cè)井?dāng)?shù)據(jù)，基于這兩項(xiàng)數(shù)據(jù)獲得目的層的構(gòu)造解釋；

對(duì)三維地震數(shù)據(jù)進(jìn)行去噪處理，得到去噪處理后的三維地震數(shù)據(jù)體；

對(duì)去噪處理后的三維地震數(shù)據(jù)體進(jìn)行相干屬性計(jì)算，得到相干屬性體；

對(duì)去噪處理后的三維地震數(shù)據(jù)體進(jìn)行曲率屬性計(jì)算，得到曲率屬性體；

對(duì)相干屬性體進(jìn)行值域校正，得到校正后的相干屬性體；

基于曲率屬性體以及校正后的相干屬性體，計(jì)算得到裂縫預(yù)測(cè)屬性體；

根據(jù)裂縫預(yù)測(cè)屬性體和目的層的構(gòu)造解釋對(duì)地下裂縫進(jìn)行預(yù)測(cè)；

對(duì)相干屬性體進(jìn)行值域校正包括：

設(shè)定門檻值，相干屬性體中采樣點(diǎn)的值小于所述門檻值時(shí)，則確定該采樣點(diǎn)進(jìn)行值域校正后的值為1；相干屬性體中采樣點(diǎn)的值大于或等于所述門檻值時(shí)，則確定該采樣點(diǎn)進(jìn)行值域校正后的值為0；

所述門檻值的取值為0≤所述門檻值≤1；

所述裂縫預(yù)測(cè)屬性體的計(jì)算公式如式13所示

Z_i＝E_i′×Q_i 式13

在式13中，i＝1，2，…，M；Z_i為裂縫預(yù)測(cè)屬性體中的采樣點(diǎn)的值；E_i′為校正后的相干屬性體中對(duì)應(yīng)位置的采樣點(diǎn)的值；Q_i為曲率屬性體中對(duì)應(yīng)位置的采樣點(diǎn)的值。

2.根據(jù)權(quán)利要求1所述的方法，其中，對(duì)三維地震數(shù)據(jù)進(jìn)行去噪處理時(shí)，所采用的去噪方法包括中值濾波法或均值濾波法。

3.根據(jù)權(quán)利要求1所述的方法，其中，對(duì)去噪處理后的三維地震數(shù)據(jù)體進(jìn)行相干屬性計(jì)算，得到相干屬性體包括以下過(guò)程：

過(guò)程1，將去噪處理后的三維地震數(shù)據(jù)體作為一個(gè)分析窗口，該分析窗口中有N道地震數(shù)據(jù)，K個(gè)采樣點(diǎn)，則所述去噪處理后的三維地震數(shù)據(jù)體表示為式3所示的矩陣D：

在式3中，d_KN為第N道的第K個(gè)樣點(diǎn)值；

過(guò)程2，按式4所示的公式計(jì)算上述去噪處理后的三維地震數(shù)據(jù)體的協(xié)方差矩陣C

在式4中，為矩陣D中的第k行向量

過(guò)程3，確定相干屬性的計(jì)算公式，由此得到相干屬性體。

4.根據(jù)權(quán)利要求3所述的方法，其中，在過(guò)程3中，所述相干屬性的計(jì)算公式如式5所示

在式5中，λ_k是協(xié)方差矩陣C的第k個(gè)特征值，λ_m是其中最大的特征值。

5.根據(jù)權(quán)利要求1所述的方法，其中，對(duì)去噪處理后的三維地震數(shù)據(jù)體進(jìn)行曲率屬性計(jì)算，得到曲率屬性體包括：

將去噪處理后的三維地震數(shù)據(jù)體轉(zhuǎn)化為傾角數(shù)據(jù)體；

根據(jù)傾角數(shù)據(jù)體計(jì)算得到數(shù)據(jù)體內(nèi)任意反射點(diǎn)的曲率。

6.根據(jù)權(quán)利要求5所述的方法，其中，將去噪處理后的三維地震數(shù)據(jù)體轉(zhuǎn)化為傾角數(shù)據(jù)體包括以下過(guò)程：

過(guò)程1，計(jì)算三維地震數(shù)據(jù)體的瞬時(shí)頻率；

過(guò)程2，計(jì)算三維地震數(shù)據(jù)體的瞬時(shí)波數(shù)；

過(guò)程3，根據(jù)得到的三維數(shù)據(jù)體的瞬時(shí)頻率和瞬時(shí)波數(shù)，計(jì)算三維地震數(shù)據(jù)體的傾角。

7.根據(jù)權(quán)利要求6所述的方法，其中，

所述過(guò)程1包括：三維地震數(shù)據(jù)體在空間任意點(diǎn)是時(shí)間標(biāo)量u(t,x,y)，根據(jù)式6所述的計(jì)算公式，計(jì)算得到瞬時(shí)頻率ω

式6中，u^H(t,x,y)為所述時(shí)間標(biāo)量u關(guān)于時(shí)間t的希爾伯特變換；

所述過(guò)程2包括：按式7所示的公式，分別計(jì)算三維地震數(shù)據(jù)體沿x方向和y方向的瞬時(shí)波數(shù)

式7中，k_x為x方向的瞬時(shí)波數(shù)，k_y為y方向的瞬時(shí)波數(shù)；

所述過(guò)程3包括：根據(jù)得到的三維數(shù)據(jù)體的瞬時(shí)頻率和瞬時(shí)波數(shù)，利用式8所示的公式計(jì)算三維地震數(shù)據(jù)體的傾角：

式8中，p_x、q_y分別為x方向和y方向的傾角分量。