[發(fā)明專(zhuān)利]波動(dòng)方程偏移直接產(chǎn)生角道集方法

有效

• 申請(qǐng)?zhí)枺?/td>: 200910238212.2
• 申請(qǐng)日: 2009-11-17
• 公開(kāi)（公告）號(hào): CN101923175A
• 公開(kāi)（公告）日: 2010-12-22
• 發(fā)明（設(shè)計(jì)）人: 張劍鋒;張輝;劉禮農(nóng)
• 申請(qǐng)（專(zhuān)利權(quán)）人: 中國(guó)科學(xué)院地質(zhì)與地球物理研究所
• 主分類(lèi)號(hào): G01V1/28
• 分類(lèi)號(hào): G01V1/28;G01V1/36
• 代理公司: 北京市中實(shí)友知識(shí)產(chǎn)權(quán)代理有限責(zé)任公司 11013
• 代理人: 李玉明

1.一種波動(dòng)方程深度偏移直接產(chǎn)生角道集方法，其特征在于采用以下步驟：A)用單邊或雙邊采集方式，由檢波器獲取人工震源激發(fā)的反射地震信號(hào)，記錄到磁帶上；B)從磁帶上讀取反射地震信號(hào)，按疊前炮集排序存放到計(jì)算機(jī)中；C)利用波動(dòng)方程偏移的單程波算法分別將炮點(diǎn)的震源波場(chǎng)和檢波點(diǎn)的接收波場(chǎng)沿深度軸進(jìn)行波場(chǎng)深度延拓，對(duì)延拓后的炮點(diǎn)和檢波點(diǎn)波場(chǎng)應(yīng)用反褶積成像條件進(jìn)行成像；D)利用單程波方法計(jì)算炮點(diǎn)的出射波對(duì)地下各成像點(diǎn)的入射角度；E)依據(jù)步驟D得到的各炮點(diǎn)出射波的入射角度，將由步驟C得到的各個(gè)單炮偏移結(jié)果按角度累加并補(bǔ)償非均勻覆蓋對(duì)幅值的影響，形成角道集；F)對(duì)角道集做剩余動(dòng)校，將出現(xiàn)明顯拉伸部分的對(duì)應(yīng)數(shù)值置為零，用奇異值分解方法壓制隨機(jī)噪音，即得到反應(yīng)介質(zhì)界面物理參數(shù)變化的疊前偏移的角道集；G)由角道集中讀取不同空間位置處成像幅值隨角度的變化，可由反演方法獲得地層的巖石物理參數(shù)，巖石物理參數(shù)圖像即能指示地下構(gòu)造的含油氣情況。

2.根據(jù)權(quán)力要求1所述的一種波動(dòng)方程深度偏移直接產(chǎn)生角道集方法，其特征在于：所述的對(duì)延拓后的炮點(diǎn)和檢波點(diǎn)波場(chǎng)應(yīng)用反褶積成像條件進(jìn)行成像是這樣實(shí)現(xiàn)的：令T₀是地震記錄的時(shí)長(zhǎng)，則頻率采樣Δω＝2π/T₀(弧度/秒)；設(shè)Δz是深度方向的空間采樣，n是深度上采樣點(diǎn)的順序數(shù)，x和y是水平位置坐標(biāo)，m是頻率采樣點(diǎn)的順序數(shù)；若炮點(diǎn)和檢波點(diǎn)波場(chǎng)經(jīng)深度延拓后在nΔz深度上分別為P_D(x，y，nΔz，mΔω)和P_U(x，y，nΔz，mΔω)，計(jì)算

A(x，y，nΔz，mΔω)＝P_D(x，y，nΔz，mΔω){P_D(x，y，nΔz，mΔω)}^*

式中上標(biāo)*是復(fù)共額，炮點(diǎn)和檢波點(diǎn)波場(chǎng)P_D和P_U的單位依據(jù)檢波器類(lèi)型可以是速度或壓力，頻率的單位是弧度/秒，深度和水平坐標(biāo)的單位為米。令A(yù)₀(nΔz)是A(x，y，nΔz，mΔω)在深度nΔz上的平均值，當(dāng)A(x，y，nΔz，mΔω)≥1.2A₀(nΔz)，計(jì)算

B(x,y,nΔz,mΔω)=PU(x,y,nΔz,mΔω){PD(x,y,nΔz,mΔω)}*A(x,y,nΔz,mΔω)]]>

當(dāng)A(x，y，nΔz，mΔω)＜1.2A₀(nΔz)，計(jì)算

B(x,y,nΔz,mΔω)=PU(x,y,nΔz,mΔω){PD(x,y,nΔz,mΔω)}*A0(nΔz)(2-a)(1+(1-a)2),]]>其中

a=A(x,y,nΔz,mΔω)A0(nΔz)]]>

令地震資料有效頻帶的下、上限分別是l₁Δω和l₂Δω，這里l₁和l₂是正整數(shù)，則空間各位置處反褶積成像條件的成像幅值φ(x，y，nΔz)為

φ(x,y,nΔz)=Σm=l1l2Re{βB(x,y,nΔz,mΔω)}]]>

若j是單位虛數(shù)，上式中做三維偏移計(jì)算時(shí)β＝j(luò)mΔω，二維偏移計(jì)算時(shí)β=jmΔω;]]>符號(hào)Re代表求復(fù)波場(chǎng)值的實(shí)部。