一種繞射波與反射波分離成像方法
【專利摘要】本發(fā)明提供了一種繞射波與反射波分離成像方法，屬于地震資料成像處理領(lǐng)域。本發(fā)明方法包括以下步驟：(1)對(duì)每一條線內(nèi)的每一個(gè)道集進(jìn)行三維τ-p正變換得到平面波數(shù)據(jù)d(py，px，y，x，τ)；(2)對(duì)由步驟(1)輸出的平面波數(shù)據(jù)d(py，px，y，x，τ)采用平面波域解構(gòu)濾波進(jìn)行處理，得到主要含繞射波能量的平面波數(shù)據(jù)體dr(py，px，y，x，τ)；(3)從步驟(2)得到的主要含繞射波能量的平面波數(shù)據(jù)體dr(py，px，y，x，τ)中抽取每一個(gè)道集對(duì)應(yīng)的三維τ-p域數(shù)據(jù)，并從原道集中獲取每個(gè)道集的炮檢點(diǎn)坐標(biāo)，然后對(duì)所述三維τ-p域數(shù)據(jù)做三維τ-p反變換得到波場(chǎng)分離后的疊前道集；(4)對(duì)由步驟(3)輸出的波場(chǎng)分離后的疊前道集進(jìn)行炮域疊前深度偏移獲得繞射波分離成像結(jié)果。本發(fā)明實(shí)現(xiàn)了繞射波信號(hào)提取和小地質(zhì)異常體或斷層的清晰成像。
【專利說明】一種繞射波與反射波分罔成像方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明屬于地震資料成像處理領(lǐng)域，具體涉及一種繞射波與反射波分離成像方法，用于繞射波信號(hào)提取和小地質(zhì)異常體或斷層清晰成像。
【背景技術(shù)】
[0002]由于地震記錄中來自小地質(zhì)異常體的繞射波信號(hào)能量較弱，要對(duì)其精確成像通常需要在成像前進(jìn)行波場(chǎng)分離或繞射波提取。目前已有的平面波域解構(gòu)濾波繞射波與反射波分離成像方法是一種有效的方法，其包含τ -ρ變換、共P平面波數(shù)據(jù)分選、平面波域解構(gòu)濾波、平面波疊前偏移等步驟。然而要使平面波疊前偏移與傳統(tǒng)炮域疊前深度偏移的成像效果相當(dāng)，需要在τ-ρ變換時(shí)分解出大量的平面波數(shù)據(jù)，這往往又為后面的偏移帶來巨大的計(jì)算量，甚至使得平面波疊前偏移的效率反而比傳統(tǒng)炮域偏移的效率還低。為此，可以將波場(chǎng)分離后的平面波數(shù)據(jù)通過τ -ρ反變換重新獲得疊前道集，再采用傳統(tǒng)炮域疊前偏移實(shí)現(xiàn)地震成像。常規(guī)二維τ-ρ正反變換所用信息僅僅是偏移距的值，故當(dāng)兩道的偏移距值相近而炮檢位置差異很大時(shí)，兩道地震數(shù)據(jù)在τ-ρ正反變換中將混疊在一起，出現(xiàn)τ-ρ反變換數(shù)據(jù)信號(hào)失真的問題，影響繞射波與反射波分離成像的效果。

【發(fā)明內(nèi)容】

[0003]本發(fā)明針對(duì)疊前偏移成像中繞射波信號(hào)弱、成像不清晰的問題，提供一種繞射波與反射波分離成像方法，根據(jù)地震數(shù)據(jù)中的炮檢坐標(biāo)信息進(jìn)行三維時(shí)間域τ -ρ正反變換，將平面波解構(gòu)濾波后的平面波域繞射波數(shù)據(jù)完整地轉(zhuǎn)換為疊前道集，進(jìn)而對(duì)繞射波信息進(jìn)行炮域疊前偏移成像，從而避免二維τ -ρ正反變換信號(hào)失真問題。
[0004]本發(fā)明是通過以下技術(shù)方案實(shí)現(xiàn)的:
[0005]一種繞射波與反射波分離成像方法，包括以下步驟:
[0006](I)對(duì)每一條線內(nèi)的每一個(gè)道集進(jìn)行三維1-1)正變換得到平面波數(shù)據(jù)(1(&，？!￡，7，
χ, τ );
[0007](2)對(duì)由步驟⑴輸出的平面波數(shù)據(jù)d(py，px, y，X，τ )采用平面波域解構(gòu)濾波進(jìn)行處理，得到主要含繞射波能量的平面波數(shù)據(jù)體c^(py，px，y，x，τ)；
[0008](3)從步驟⑵得到的主要含繞射波能量的平面波數(shù)據(jù)體dr(py，px，y，X，τ)中抽取每一個(gè)道集對(duì)應(yīng)的三維τ -ρ域數(shù)據(jù)，并從原道集中獲取每個(gè)道集的炮檢點(diǎn)坐標(biāo)，然后對(duì)所述三維τ-ρ域數(shù)據(jù)做三維τ-ρ反變換得到波場(chǎng)分離后的疊前道集；
[0009](4)對(duì)由步驟(3)輸出的波場(chǎng)分離后的疊前道集進(jìn)行炮域疊前深度偏移獲得繞射波分離成像結(jié)果(即分離掉了反射波而只包括繞射波的成像結(jié)果)。
[0010]所述步驟⑴具體如下:
[0011](11)采用下面的三維τ -P正變換公式對(duì)每一條線內(nèi)的每一個(gè)道集進(jìn)行三維τ -ρ正變換得到τ -ρ域數(shù)據(jù)ψ:`
[0012]WiPr^ = YM: ” A +PyK)[0013]其中，C^i表示第i道數(shù)據(jù)，hix、hiy是第i道數(shù)據(jù)的兩個(gè)方向的偏移距分量，此兩個(gè)偏移距分量是由第i道數(shù)據(jù)的炮檢點(diǎn)坐標(biāo)計(jì)算得到的，px、Py表示兩個(gè)方向的射線參數(shù)；
[0014](12)把步驟(11)得到的τ-ρ域數(shù)據(jù)V存入五維數(shù)據(jù)體內(nèi)形成平面波數(shù)據(jù)d (py，Px，y，X，τ)(一個(gè)三維工區(qū)包括N條線，每條線包括M個(gè)道集，Ψ是一個(gè)道集的變換結(jié)果，所有V就組成了五維數(shù)據(jù)體)，其中，y代表線號(hào)，X代表線內(nèi)道集號(hào)。
[0015]所述步驟(3)中所述從步驟(2)得到的平面波數(shù)據(jù)體dr(py，px, y，X，τ )中抽取每一個(gè)道集對(duì)應(yīng)的三維τ -ρ域數(shù)據(jù)是這樣實(shí)現(xiàn)的:按線號(hào)y和道集號(hào)χ從平面波數(shù)據(jù)體dr(Py，Px，y，X，τ )取值得到每一個(gè)道集對(duì)應(yīng)的三維τ-ρ域數(shù)據(jù)。
[0016]所述步驟(3)中對(duì)所述三維τ -ρ域數(shù)據(jù)做三維τ -ρ反變換得到波場(chǎng)分離后的疊前道集是這樣實(shí)現(xiàn)的:
[0017]采用下面的三維τ -ρ反變換公式將所述三維τ -ρ域數(shù)據(jù)變換為波場(chǎng)分離后的疊前道集$i(t):
【權(quán)利要求】
1.一種繞射波與反射波分離成像方法，其特征在于:所述方法包括以下步驟: (1)對(duì)每一條線內(nèi)的每一個(gè)道集進(jìn)行三維τ-P正變換得到平面波數(shù)據(jù)d(py，px, y, x,τ ); (2)對(duì)由步驟⑴輸出的平面波數(shù)據(jù)d(py，px，y，x，τ)采用平面波域解構(gòu)濾波進(jìn)行處理，得到主要含繞射波能量的平面波數(shù)據(jù)體cUpy，px，y，X，τ)； (3)從步驟(2)得到的主要含繞射波能量的平面波數(shù)據(jù)體cUpy，px，y，x，τ)中抽取每一個(gè)道集對(duì)應(yīng)的三維τ -P域數(shù)據(jù)，并從原道集中獲取每個(gè)道集的炮檢點(diǎn)坐標(biāo)，然后對(duì)所述三維τ-P域數(shù)據(jù)做三維τ-P反變換得到波場(chǎng)分離后的疊前道集； (4)對(duì)由步驟(3)輸出的波場(chǎng)分離后的疊前道集進(jìn)行炮域疊前深度偏移獲得繞射波分離成像結(jié)果。
2.根據(jù)權(quán)利要求1所述的繞射波與反射波分離成像方法，其特征在于:所述步驟(1)具體如下: (11)采用下面的三維τ-P正變換公式對(duì)每一條線內(nèi)的每一個(gè)道集進(jìn)行三維τ-ρ正變換得到τ -P域數(shù)據(jù)ψ:
3.根據(jù)權(quán)利要求2所述的繞射波與反射波分離成像方法，其特征在于:所述步驟(3)中所述從步驟(2)得到的平面波數(shù)據(jù)體cUpy，px, y, χ, τ)中抽取每一個(gè)道集對(duì)應(yīng)的三維τ-P域數(shù)據(jù)是這樣實(shí)現(xiàn)的淑線號(hào)丫和道集號(hào)乂從平面波數(shù)據(jù)體火匕也小^ τ)取值得到每一個(gè)道集對(duì)應(yīng)的三維τ-ρ域數(shù)據(jù)。
4.根據(jù)權(quán)利要求3所述的繞射波與反射波分離成像方法，其特征在于:所述步驟(3)中對(duì)所述三維τ-ρ域數(shù)據(jù)做三維τ-ρ反變換得到波場(chǎng)分離后的疊前道集是這樣實(shí)現(xiàn)的: 采用下面的三維τ-P反變換公式將所述三維τ-ρ域數(shù)據(jù)變換為波場(chǎng)分離后的疊前道集 Φ“?):
【文檔編號(hào)】G01V1/28GK103675896SQ201210316656
【公開日】2014年3月26日 申請(qǐng)日期:2012年8月30日 優(yōu)先權(quán)日:2012年8月30日
【發(fā)明者】段心標(biāo), 方伍寶, 呂秋玲, 王鵬燕, 宋林 申請(qǐng)人:中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術(shù)研究院