專利名稱：一種提高垂直地震剖面成像質(zhì)量的方法
技術(shù)領(lǐng)域：
本發(fā)明涉及石油物探技術(shù)，是垂直地震剖面(VSP)數(shù)據(jù)處理方
法，具體為提高垂直地震剖面成像質(zhì)量的方法。
背景技術(shù)：
常規(guī)VSP資料處理通常包括初至拾取、速度反演、波場(chǎng)分離、成 像等。VSP資料的成像處理一般基于水平速度模型。然而，在反射界 面存在傾角的情況下，基于水平層狀速度模型進(jìn)行VSP成像會(huì)使歸
位不準(zhǔn)確。對(duì)這一問(wèn)題通常的考慮是通過(guò)傾角時(shí)差校正或迭代成像來(lái)
解決。VSP數(shù)據(jù)的傾角時(shí)差校正方法是基于地下介質(zhì)為單一均勻傾斜 層的假設(shè)，并且需要先由傾角測(cè)井或VSP資料獨(dú)立地估算出反射界 面的傾角來(lái)。VSP數(shù)據(jù)的迭代成像實(shí)際上隱含著對(duì)反射界面的傾角估
計(jì)(通過(guò)不斷修改速度模型來(lái)降低反射界面傾角對(duì)成像結(jié)果的影響)。 如果能在成像前獲得反射界面的傾角信息，則可以形成一個(gè)合乎邏輯
的VSP數(shù)據(jù)處理思路，即基于速度模型(反射層的幾何形態(tài)和速度 屬性)的VSP數(shù)據(jù)處理。
VSP的波場(chǎng)結(jié)構(gòu)與反射界面的傾角緊密相關(guān)，因此，可通過(guò)VSP 數(shù)據(jù)來(lái)估算反射界面的傾角，用以修正速度模型，來(lái)提高VSP成像質(zhì)量。
Kennett、 Wyatt等曾提出利用VSP資料基于解析公式估算反射 界面傾角的方法。這一方法實(shí)際上是通過(guò)測(cè)量?jī)蓚€(gè)接收深度位置之間的傾角時(shí)差來(lái)計(jì)算反射界面傾角。隨后，Malloy、 N叩onen等將以一 道或兩道傾斜時(shí)差進(jìn)行反射界面傾角估計(jì)變?yōu)橐远嗟纼A斜時(shí)差進(jìn)行 反射界面的傾角估計(jì)，提高了傾角估算的準(zhǔn)確性。但這類傾斜時(shí)差法 的理論基礎(chǔ)實(shí)際上是基于單層傾斜平反射界面，這意味著實(shí)際估計(jì)反 射界面傾角時(shí)要使用平均速度，沒(méi)有考慮到射線偏折的影響，使得估 算的反射界面傾角存在較大誤差。
為了避免使用平均速度對(duì)反射界面傾角估算結(jié)果的影響，Lins 等提出了基于層速度模型進(jìn)行反射界面傾角估算的模型比較方法。以 零傾角作為層速度模型的初始估計(jì)值，計(jì)算射線從激發(fā)點(diǎn)傳播到接收 點(diǎn)的旅行時(shí)，再比較根據(jù)模型計(jì)算的旅行時(shí)與實(shí)際觀測(cè)的旅行時(shí)之 差，不斷調(diào)整模型的傾角估計(jì)值，直到旅行時(shí)差值達(dá)到給出的精度范 圍內(nèi)。
為了減少反演計(jì)算的變量個(gè)數(shù)和反演過(guò)程中的多解性，Lins采 用所謂"剝層"法，從淺層到深層逐步進(jìn)行反射界面的傾角估算。但 用"剝層"法反演深層反射界面的傾角時(shí)，淺層誤差會(huì)累積進(jìn)去，使 得對(duì)于深層反射界面傾角的估算結(jié)果并不理想。另一方面，相比于初 至波旅行時(shí)， 一次反射波旅行時(shí)對(duì)于界面傾角的變化敏感，使用一次 反射波旅行時(shí)進(jìn)行反射界面傾角反演結(jié)果會(huì)更精確。但在Lins的模 型比較方法中，由于進(jìn)行比較的是射線旅行時(shí)間，這樣，多次波的干 擾及在界面附近下行直達(dá)波的干涉就會(huì)影響到一次反射波旅行時(shí)的 準(zhǔn)確拾取，從而影響到反射界面傾角的反演成像結(jié)果。利用VSP數(shù)據(jù)估算反射界面傾角還有一類偏振測(cè)量法(Spencer 等)，通過(guò)對(duì)三分量VSP數(shù)據(jù)進(jìn)行偏振分析估算出射線的傳播方向， 再由射線傳播方向和反射界面傾角的關(guān)系獲得反射界面的傾角。這一 方法估算反射層傾角的準(zhǔn)確性受偏振測(cè)量誤差的影響很大。但由于波 場(chǎng)干涉的影響，要想準(zhǔn)確確定波場(chǎng)的偏振方向并非易事。 發(fā)明目的
本發(fā)明目的是提供一種用修正速度模型，提高利用垂直地震剖面 成像質(zhì)量的方法。

發(fā)明內(nèi)容
本發(fā)明的具體實(shí)施步驟為
1、 在地表人工激發(fā)地震波，在井中放置地震波接收器串接收地 震波場(chǎng)信號(hào)，再由電纜傳到地面記錄儀器記錄VSP地震波波場(chǎng)；
2、 在零偏及非零偏VSP地震波波場(chǎng)數(shù)據(jù)上拾取不同接收點(diǎn)處的 下行縱波初至?xí)r間；
3、 對(duì)VSP地震波波場(chǎng)進(jìn)行波場(chǎng)分離，得到實(shí)際觀測(cè)的反射波場(chǎng);
4、 根據(jù)已知的聲波測(cè)井資料及零偏VSP下行縱波初至?xí)r間建立 無(wú)傾角的初始層速度模型；
步驟4)所述的建立無(wú)傾角的初始層速度模型為通過(guò)聲波測(cè)井資 料劃分層位，利用零偏VSP下行縱波初至?xí)r間反演各層速度。
5、 根據(jù)下式將非零偏VSP縱波初至?xí)r間外推至井旁點(diǎn)，獲得下
行波的時(shí)間場(chǎng)t(x，z):
/(x,= max{f(/z) - r(/z; x, z)}式中"L(h;X，Z)是用射線追蹤方法確定的下行波由深度h處接收
點(diǎn)到井旁點(diǎn)(x，z)的旅行時(shí)間，t(h)指非零偏VSP深度h處接收點(diǎn)下行
縱波初至?xí)r間，X指井旁點(diǎn)橫坐標(biāo)，Z指井旁點(diǎn)縱坐標(biāo)，h代表接收
點(diǎn)深度，t和T用來(lái)表示時(shí)間，max表示找到沿不同射線路徑所得一 系列t(h)與i;(h;x，z)時(shí)間差中最大者；
步驟5)所述射線追蹤過(guò)程采用的是對(duì)初始層速度模型用界面傾 角假設(shè)值cp修正后的層速度模型。
6、計(jì)算下行波時(shí)間場(chǎng)t(x，z)沿水平x和垂直z兩個(gè)方向上的偏導(dǎo) 數(shù)§和§，并根據(jù)下式來(lái)求取下行波射線角度a(x，z):
a(x，z) = arcsin
其中arcsin表示反正弦函數(shù)；
7、 利用下行波射線角度a(x，z)求取出上行波射線角度P(x，z):
風(fēng);c, z) = p + arcsm"j-^-f
式中Sin表示正弦函數(shù)。V。入射速度，Vp為反射速度；Cp為界 面傾角假設(shè)值，X、 Z分別為計(jì)算點(diǎn)處的水平坐標(biāo)和垂直坐標(biāo)。
步驟7)所述的入射速度，反射速度(Va、 Vp)是縱波(P)波
速度并且相等。
8、 利用下行波的時(shí)間場(chǎng)t(x，z)、上行波射線角度0(x，z)，采用射 線追蹤方法形成理論上行波場(chǎng)u"h，t);
9、 根據(jù)下式計(jì)算對(duì)應(yīng)界面傾角cp時(shí)u"h，t)與u(h，t)的相關(guān)性度
量K(h，Cp)值式中利用了理論上行波場(chǎng)U"h，t)和實(shí)際觀測(cè)的反射波場(chǎng) U(h，t)，其中h表深度，t為時(shí)間，(p為界面傾角，W表示3維實(shí)空間。
步驟9)所述的K(h，cp)是在近井范圍內(nèi)沿依賴于u"h，t)的射線 路徑對(duì)上行波場(chǎng)u(h，t)求和。
10、 改變界面傾角的值cp，重復(fù)步驟5到步驟9，以獲取一系列
的對(duì)應(yīng)不同界面傾角Cp的K(h，Cp)值；
11、 在一系列的K(h，(p)值中拾取最大，其所對(duì)應(yīng)的界面傾角為 所要求取的井旁反射界面傾角；
12、 用井旁反射界面傾角值對(duì)步驟4中所得初始層速度模型進(jìn)行 修正，利用修改后的帶有界面傾角信息的層速度模型進(jìn)行垂直地震剖 面成像。
本發(fā)明通過(guò)外推下行縱波(P)波初至?xí)r間以獲得接收點(diǎn)附近上 行P波時(shí)間場(chǎng)的方法，可有效降低上覆地層對(duì)反射界面傾角估計(jì)結(jié)果 的影響和一次反射波旅行時(shí)拾取誤差對(duì)反射界面傾角估計(jì)結(jié)果的影 響，提高了利用垂直地震剖面成像質(zhì)量。


圖1為用來(lái)產(chǎn)生理論VSP數(shù)據(jù)的二維構(gòu)造模型。 圖2為根據(jù)圖1給出的模型所產(chǎn)生的理論三分量VSP上行波場(chǎng)。 圖3為傾角掃描結(jié)果。其中a為傾角掃描圖譜，b為根據(jù)振幅最 大準(zhǔn)則在圖譜中所拾取的各反射層傾角值。
具體實(shí)施例方式
以下結(jié)合附圖進(jìn)一步說(shuō)明本發(fā)明。
在地表人工激發(fā)地震波，在井中放置地震波接收器串接收地震波
場(chǎng)信號(hào)，再由電纜傳到地面記錄儀器記錄VSP地震波波場(chǎng)。激發(fā)位置
包括零偏和非零偏兩種方式。
本發(fā)明實(shí)施例中使用的數(shù)據(jù)來(lái)自于一個(gè)二維構(gòu)造模型。圖1給出
了該構(gòu)造模型，圖2是根據(jù)這一模型所產(chǎn)生的理論三分量VSP上行 波場(chǎng)數(shù)據(jù)，炮點(diǎn)到井口的距離為IOOO米。
分別在所記錄的零偏和非零偏兩種數(shù)據(jù)中拾取下行縱波初至。然 后，通過(guò)聲波測(cè)井資料劃分層位，利用零偏VSP下行縱波初至?xí)r間 反演各層速度，建立無(wú)傾角信息的初始層速度模型，進(jìn)行傾角掃描時(shí) 所采用的層速度模型數(shù)據(jù)如下表1所示；
序號(hào)X坐標(biāo) (m)Y坐標(biāo) (m)z坐標(biāo) (m)P波速 度(m/s)s波速 度(m/s)界面傾 角(° )
10.00.00.01750.00712.7300.0
20.00.0300.01800.00717.420.0
30.00.0700.02200.001044.800.0
40,00.01200.01900.00882.850.0
50.00.01800.02600.00914.150.0
60.00.02100.03000.00953.880.0
70.00.02400.03400.001024.560.0
表中的坐標(biāo)代表層界面與井軌跡交點(diǎn)處的坐標(biāo)值。
對(duì)非零偏VSP數(shù)據(jù)進(jìn)行波場(chǎng)分離，以得到實(shí)際觀測(cè)的反射波場(chǎng)， 如圖2所示，根據(jù)圖1給出的模型所產(chǎn)生的理論三分量VSP上行波 場(chǎng)中的X、 Z分量?？v坐標(biāo)為深度(米)，橫坐標(biāo)為時(shí)間(毫秒)。
給出界面傾角的假設(shè)值cp，根據(jù)下式將所拾取的非零偏VSP初至?xí)r間外推至井旁點(diǎn)(x， z)處，以獲得下行波的時(shí)間場(chǎng)t(x，z): f (x， z) = max{/(/z) - x， z)}
式中x指井旁點(diǎn)橫坐標(biāo)，z指井旁點(diǎn)縱坐標(biāo)，h代表接收點(diǎn)深度，t 和t用來(lái)表示時(shí)間。"c(h;x，z)是用射線追蹤方法確定的下行波由深度h 處接收點(diǎn)到井旁點(diǎn)(x， z)的旅行時(shí)間，射線追蹤過(guò)程采用的是對(duì)初始 層速度模型用界面傾角假設(shè)值cp修正后的層速度模型；t(h)指非零偏 VSP深度h處接收點(diǎn)下行波初至?xí)r間；max表示找到沿不同射線路徑 所得一系列t(h)與T(h;x，z)時(shí)間差中最大者。
計(jì)算下行波時(shí)間場(chǎng)t(x，z)沿水平x和垂直z兩個(gè)方向上的偏導(dǎo)數(shù) ^和^，并根據(jù)下式來(lái)求取下行波射線角度a(x，z):
or(;c,z) = arcsin
其中arcsin表示反正弦函數(shù)。
然后，依據(jù)Snell定律，利用下行波射線角度a(x，z)求取出上行波 射線角度P(x，z):
sin[p + "(x，z)]1 "(x， z) = p + arcsin〗~^-^-^
這里sin表示正弦函數(shù)。V。入射速度，Vp為反射速度；cp為界面傾角 假設(shè)值，x、 z分別為計(jì)算點(diǎn)處的水平坐標(biāo)和垂直坐標(biāo)。本實(shí)施例中 采用的是縱波波場(chǎng)，則V。、 Vp都是縱波速度并且相等。
利用上面得到的下行波的時(shí)間場(chǎng)t(x，z)、上行波射線角度P (x，z)， 采用射線追蹤方法形成理論上行波場(chǎng)u"h，t)。射線追蹤過(guò)程采用的 是對(duì)初始層速度模型用界面傾角假設(shè)值cp修正后的層速度模型。計(jì)算對(duì)應(yīng)界面傾角假設(shè)值cp時(shí)U"h，t)與實(shí)際觀測(cè)的反射波場(chǎng) U(h,t)的相關(guān)性度量K(h,Cp)值
式中h表深度，t為時(shí)間，cp為界面傾角，W表示3維實(shí)空間。這里，
K(h，cp)是在近井范圍內(nèi)沿依賴于u"h，t)的射線路徑對(duì)上行波場(chǎng)u(h，t)求和。
改變界面傾角的假設(shè)值cp,重復(fù)上面的步驟，可以獲取一系列的對(duì)應(yīng) 不同界面傾角假設(shè)值cp的K(h，cp)值，這里稱之為傾角掃描譜圖，如圖 3中a所示，圖中縱坐標(biāo)為深度(米)，橫坐標(biāo)為角度(度)。在傾角 掃描譜圖中最大振幅值位置處對(duì)應(yīng)的傾角即為井旁反射界面的傾角。
最后，用依據(jù)振幅最大準(zhǔn)則在傾角掃描譜圖中所拾取的各反射界 面傾角值(如圖3中b所示)對(duì)初始層速度模型進(jìn)行修正，利用修改 后的帶有界面傾角信息的層速度模型進(jìn)行垂直地震剖面成像。
本發(fā)明實(shí)施例中使用的數(shù)據(jù)來(lái)自于一個(gè)二維構(gòu)造模型。圖1給出 了該構(gòu)造模型，圖2是根據(jù)這一模型所產(chǎn)生的理論三分量VSP上行 波場(chǎng)數(shù)據(jù)，炮點(diǎn)到井口的距離為1000米。本發(fā)明使用圖2中的理論 上行波場(chǎng)數(shù)據(jù)及表1的層速度模型數(shù)據(jù)作為輸入，對(duì)反射界面的傾角 進(jìn)行掃描。傾角掃描結(jié)果如圖3所示，其中a為掃描圖譜，b為根據(jù) 振幅最大準(zhǔn)則在圖譜中所拾取的各反射層傾角值。將圖3中b所示的 傾角拾取結(jié)果和圖1給出的理論模型中各層傾角進(jìn)行對(duì)比，可知拾取 的傾角值是準(zhǔn)確的，表明本發(fā)明是行之有效的。
權(quán)利要求
1、一種提高垂直地震剖面成像質(zhì)量的方法，其特征在于具體實(shí)施步驟為1)在地表人工激發(fā)地震波，在井中放置地震波接收器串接收地震波場(chǎng)信號(hào)，再由電纜傳到地面記錄儀器記錄VSP地震波波場(chǎng)；2)在零偏及非零偏VSP地震波波場(chǎng)數(shù)據(jù)上拾取不同接收點(diǎn)處的下行縱波初至?xí)r間；3)對(duì)VSP地震波波場(chǎng)進(jìn)行波場(chǎng)分離，得到實(shí)際觀測(cè)的反射波場(chǎng)；4)根據(jù)已知的聲波測(cè)井資料及零偏VSP下行縱波初至?xí)r間建立無(wú)傾角的初始層速度模型；5)根據(jù)下式將非零偏VSP縱波初至?xí)r間外推至井旁點(diǎn)，獲得下行波的時(shí)間場(chǎng)t(x，z)<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><mi>t</mi><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>z</mi> <mo>)</mo></mrow><mo>=</mo><munder> <mi>max</mi> <mi>h</mi></munder><mo>{</mo><mi>t</mi><mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo></mrow><mo>-</mo><mi>&tau;</mi><mrow> <mo>(</mo> <mi>h</mi> <mo>;</mo> <mi>x</mi> <mo>,</mo> <mi>z</mi> <mo>)</mo></mrow><mo>}</mo> </mrow>]]></math></maths>式中τ(h；x，z)是用射線追蹤方法確定的下行波由深度h處接收點(diǎn)到井旁點(diǎn)(x，z)的旅行時(shí)間，t(h)指非零偏VSP深度h處接收點(diǎn)下行縱波初至?xí)r間，x指井旁點(diǎn)橫坐標(biāo)，z指井旁點(diǎn)縱坐標(biāo)，h代表接收點(diǎn)深度，t和τ用來(lái)表示時(shí)間，max表示找到沿不同射線路徑所得一系列t(h)與τ(h；x，z)時(shí)間差中最大者；6)計(jì)算下行波時(shí)間場(chǎng)t(x，z)沿水平x和垂直z兩個(gè)方向上的偏導(dǎo)數(shù) id="icf0002" file="A2008101167090002C2.tif" wi="3" he="5" top= "226" left = "29" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>和 id="icf0003" file="A2008101167090002C3.tif" wi="5" he="5" top= "226" left = "38" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>并根據(jù)下式來(lái)求取下行波射線角度α(x，z)<maths id="math0002" num="0002" ><math><![CDATA[ <mrow><mi>&alpha;</mi><mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>z</mi> <mo>)</mo></mrow><mo>=</mo><mi>arcsin</mi><mo>[</mo><mfrac> <mfrac><mrow> <mo>&PartialD;</mo> <mi>t</mi></mrow><mrow> <mo>&PartialD;</mo> <mi>x</mi></mrow> </mfrac> <msqrt><msup> <mrow><mo>(</mo><mfrac> <mrow><mo>&PartialD;</mo><mi>t</mi> </mrow> <mrow><mo>&PartialD;</mo><mi>x</mi> </mrow></mfrac><mo>)</mo> </mrow> <mn>2</mn></msup><mo>+</mo><msup> <mrow><mo>(</mo><mfrac> <mrow><mo>&PartialD;</mo><mi>t</mi> </mrow> <mrow><mo>&PartialD;</mo><mi>z</mi> </mrow></mfrac><mo>)</mo> </mrow> <mn>2</mn></msup> </msqrt></mfrac><mo>]</mo> </mrow>]]></math></maths>其中arcsin表示反正弦函數(shù)；7)利用下行波射線角度α(x，z)求取出上行波射線角度β(x，z)式中sin表示正弦函數(shù)。Vα入射速度，Vβ為反射速度； id="icf0006" file="A2008101167090003C2.tif" wi="3" he="3" top= "70" left = "161" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>為界面傾角假設(shè)值，x、z分別為計(jì)算點(diǎn)處的水平坐標(biāo)和垂直坐標(biāo)；8)利用下行波的時(shí)間場(chǎng)t(x，z)、上行波射線角度β(x，z)，采用射線追蹤方法形成理論上行波場(chǎng) id="icf0007" file="A2008101167090003C3.tif" wi="17" he="5" top= "102" left = "93" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>9)根據(jù)下式計(jì)算對(duì)應(yīng)界面傾角 id="icf0008" file="A2008101167090003C4.tif" wi="3" he="3" top= "114" left = "104" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>時(shí) id="icf0009" file="A2008101167090003C5.tif" wi="14" he="5" top= "113" left = "114" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>與u(h，t)的相關(guān)性度量 id="icf0010" file="A2008101167090003C6.tif" wi="14" he="5" top= "124" left = "33" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>值式中理論上行波場(chǎng) id="icf0012" file="A2008101167090003C8.tif" wi="14" he="5" top= "146" left = "84" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>和實(shí)際觀測(cè)的反射波場(chǎng)u(h，t)，其中h表深度，t為時(shí)間， id="icf0013" file="A2008101167090003C9.tif" wi="3" he="3" top= "159" left = "80" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>為界面傾角，R3表示3維實(shí)空間；10)改變界面傾角的值 id="icf0014" file="A2008101167090003C10.tif" wi="4" he="3" top= "170" left = "87" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>重復(fù)步驟5到步驟9，以獲取一系列的對(duì)應(yīng)不同界面傾角 id="icf0015" file="A2008101167090003C11.tif" wi="3" he="3" top= "181" left = "72" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>的 id="icf0016" file="A2008101167090003C12.tif" wi="14" he="5" top= "179" left = "81" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>值；11)在一系列的 id="icf0017" file="A2008101167090003C13.tif" wi="14" he="5" top= "190" left = "74" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/>值中拾取最大，其所對(duì)應(yīng)的界面傾角為所要求取的井旁反射界面傾角；12)用井旁反射界面傾角值對(duì)步驟4中所得初始層速度模型進(jìn)行修正，利用修改后的帶有界面傾角信息的層速度模型進(jìn)行垂直地震剖面成像。
2. 根據(jù)權(quán)利要求1所述的提高垂直地震剖面成像質(zhì)量的方法， 其特征在于步驟4)所述的建立無(wú)傾角的初始層速度模型為通過(guò)聲波 測(cè)井資料劃分層位，利用零偏VSP下行縱波初至?xí)r間反演各層速度。根據(jù)權(quán)利要求1所述的提高垂直地震剖面成像質(zhì)量的方法，其特征在于步驟5)所述射線追蹤過(guò)程采用的是對(duì)初始層速度模型用界面傾角假設(shè)值(P修正后的層速度模型。根據(jù)權(quán)利要求1所述的提高垂直地震剖面成像質(zhì)量的方法，其特征在于步驟7)所述的入射速度，反射速度(V。、 Vp)是縱波(P) 波速度并且相等。根據(jù)權(quán)利要求1所述的提高垂直地震剖面成像質(zhì)量的方法，其特 征在于步驟9)所述的K(h，cp)是在近井范圍內(nèi)沿依賴于u"h,t)的射 線路徑對(duì)上行波場(chǎng)u(h，t)求和。
全文摘要
本發(fā)明是石油物探中提高垂直地震剖面成像質(zhì)量的方法，在零偏及非零偏垂直地震剖面波場(chǎng)數(shù)據(jù)上拾取不同接收點(diǎn)處的下行縱波初至?xí)r間，建立無(wú)傾角的初始層速度模型，縱波初至?xí)r間外推至井旁點(diǎn)，獲得下行波的時(shí)間場(chǎng)，求取下、上行波射線角度利用下行波射線角度，采用射線追蹤方法形成理論上行波場(chǎng)，獲取一系列的對(duì)應(yīng)不同界面傾角值，確定最大對(duì)應(yīng)的界面傾角為所要求取的井旁反射界面傾角。本發(fā)明可有效降低上覆地層對(duì)反射界面傾角估計(jì)結(jié)果的影響和一次反射波旅行時(shí)拾取誤差對(duì)反射界面傾角估計(jì)結(jié)果的影響，提高利用垂直地震剖面成像質(zhì)量。
文檔編號(hào)G01V1/28GK101630016SQ200810116709
公開日2010年1月20日 申請(qǐng)日期2008年7月16日 優(yōu)先權(quán)日2008年7月16日
發(fā)明者侯愛(ài)源, 王寶彬 申請(qǐng)人:中國(guó)石油天然氣集團(tuán)公司;中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司