[發明專利]一種巖性識別方法及識別系統

有效

• 申請號: 201711319523.2
• 申請日: 2017-12-12
• 公開（公告）號: CN108680955B
• 公開（公告）日: 2019-09-17
• 發明（設計）人: 李岳桐;盧宗盛
• 申請（專利權）人: 中國地質大學（武漢）
• 主分類號: G01V1/30
• 分類號: G01V1/30;G01V1/40
• 代理公司: 武漢知產時代知識產權代理有限公司 42238
• 代理人: 龔春來

1.一種巖性識別方法，其特征在于，包括如下步驟：

S1、獲取測井資料；

S2、獲得巖性曲線，以及根據所述測井資料計算獲得波阻抗曲線A；

S3、根據所述巖性曲線獲得至少一個巖性的巖性敏感曲線B；

S4、對所述波阻抗曲線A進行低通濾波處理，通過公式(1)計算得到低頻趨勢波阻抗A_lpf，以及對所述敏感曲線B進行高通濾波處理，通過公式(2)計算得到高頻敏感曲線B_hpf；所述公式(1)為：A_lpf＝F(A)，所述公式(2)為：B_hpf＝F(B)；其中，F(A)、F(B)均為傅里葉變換函數；

以及S5、根據所述低頻趨勢波阻抗A_lpf與高頻敏感曲線B_hpf，首先，將所述高頻敏感曲線Bhpf轉化到高頻波阻抗Ahpf值域范圍內，得到高頻敏感波阻抗Chpf；然后，在頻率域內根據公式(5)對所述高頻敏感波阻抗Chpf與高頻波阻抗Ahpf進行加權計算，計算得到巖性識別敏感因子；最后，通過所述巖性識別敏感因子進行地震反演，以對巖性進行識別；所述公式(5)為：

S＝f[(F(Ahpf)+F(Chpf)]；

其中，F(x)為傅里葉變換函數，f(x)為傅里葉逆變換函數。

2.根據權利要求1所述的巖性識別方法，其特征在于，步驟S1中，所述測井資料包括：聲波時差、密度、自然伽馬、自然電位中的一種或幾種。

3.根據權利要求1所述的巖性識別方法，其特征在于，所述步驟S5包括：

S51、統計高頻波阻抗A_hpf值域變化范圍，得到其最小值A_min和最大值A_max，以及統計高頻敏感曲線B_hpf值域變化范圍，得到其最小值B_min和最大值B_max；

S52、將所述高頻敏感曲線B_hpf與所述高頻波阻抗A_hpf進行交匯分析，根據公式(3)擬合求取相關函數Fx(B_hpf)，所述公式(3)為：Fx(B_hpf)＝a*(B_hpf)²+b*(B_hpf)+c；其中，所述a、b、c均為相關常數；

S53、根據公式(4)將所述高頻敏感曲線B_hpf轉化到所述高頻波阻抗A_hpf值域范圍內，得到高頻敏感波阻抗C_hpf；所述公式(4)為：

以及S54、在頻率域內根據公式(5)對所述高頻敏感波阻抗C_hpf與高頻波阻抗A_hpf進行加權計算，由此得到巖性敏感因子S。

4.根據權利要求3所述的巖性識別方法，其特征在于，所述步驟S52中，所述高頻敏感曲線B_hpf與高頻波阻抗A_hpf的相關系數大于或等于0.7。

5.一種實現權利要求1-4任一項所述巖性識別方法的巖性識別系統，其特征在于，包括：

第一曲線獲取單元，其用于獲取巖性曲線，以及根據所述測井資料計算獲得波阻抗曲線A；

第二曲線獲取單元，其用于根據所述巖性曲線獲得至少一個巖性的巖性敏感曲線B；

第一曲線處理單元，其用于對所述波阻抗曲線A進行低通濾波處理，計算得到低頻趨勢波阻抗A_lpf；

第二曲線處理單元，其用于對所述敏感曲線B進行高通濾波處理，計算得到高頻敏感曲線B_hpf；

以及敏感因子獲取單元，其用于根據所述低頻趨勢波阻抗A_lpf與高頻敏感曲線B_hpf計算得到巖性識別敏感因子。

6.根據權利要求5所述的巖性識別系統，其特征在于，還包括：數據輸出單元，其用于根據所述巖性識別敏感因子輸出巖性統計數據。