[發(fā)明專利]全波形反演噪音濾波算子的提取及使用其噪音濾波的方法

在審

• 申請?zhí)枺?/td>: 201310522885.7
• 申請日: 2013-10-29
• 公開（公告）號: CN104570090A
• 公開（公告）日: 2015-04-29
• 發(fā)明（設(shè)計）人: 孟憲軍;王玉梅;慎國強;何兵紅;王希萍;苗永康;王振濤;鈕學(xué)民;張洪;張達
• 申請（專利權(quán)）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/36
• 代理公司: 濟南日新專利代理事務(wù)所 37224
• 代理人: 謝省法

1.全波形反演噪音濾波算子的提取方法，其特征在于包括如下步驟：

1）在全波形反演中獲取單炮的全波形反演梯度之后，通過二維傅里葉變換得到多尺度二維波數(shù)下全波形反演參數(shù)的波數(shù)域分布特征，從而將數(shù)據(jù)從空間域轉(zhuǎn)化到波數(shù)域；

2）針對全波形反演速度噪音在波數(shù)域中的可存在區(qū)間，以及有效全波形反演速度在波數(shù)域中的分布，確立波數(shù)域中垂直波數(shù)與水平波數(shù)的關(guān)系；

3）將地質(zhì)構(gòu)造的角度引入了虛假全波形反演信息，當(dāng)?shù)刭|(zhì)構(gòu)造傾角為θ時，在波數(shù)域中滿足：tanθ=k_x/k_z，因此當(dāng)角度大于地質(zhì)構(gòu)造傾角θ時，存在關(guān)系式k_x＞k_ztanθ，在該區(qū)間內(nèi)為噪音，將該區(qū)間定義為噪音區(qū)；在k_x＜k_ztanθ時，有效全波形反演梯度占據(jù)主導(dǎo)地位，將該區(qū)間定義為有效區(qū)；在有效區(qū)與噪音區(qū)之間建立過渡區(qū)，區(qū)角度范圍為

其中，θ為地層傾角，k_x為水平波數(shù)，k_z為垂直波數(shù)；為微小擾動角度；

4）建立全波形反演噪音濾波算子

在理想情況下期望大于地質(zhì)地層角度的噪音完全消除，同時小于地層傾角的有效信息完全保留，該種情況下波數(shù)域濾波算子為：

h(kx,kz)=1kx≤kztanθ0kx>kztanθ---(1)]]>

式中：h(k_x,k_z)為波數(shù)域濾波算子，θ為地層傾角，k_x為水平波數(shù)，k_z為垂直波數(shù)；

h(k_x,k_z)一階導(dǎo)數(shù)為非連續(xù)函數(shù)，傅里葉變換存在嚴重的吉普斯現(xiàn)象，引入從噪音區(qū)到有效區(qū)的過渡區(qū)域，過渡區(qū)角度范圍為在該角度區(qū)間的斷點處波數(shù)域濾波算子為:

式中：θ為地層傾角，為微小擾動角度，k_x為水平波數(shù)，k_z為垂直波數(shù)；

根據(jù)線性擬合函數(shù)：

式中：G為濾波算子值隨水平波數(shù)變化的斜率，P為水平波數(shù)為零時濾波算子值。

則全波形反演噪音濾波算子可表示為:

其中

引入擾動量ε，綜上，可得穩(wěn)定的全波形反演噪音濾波算子：

其中