• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    VSP時空域高角度單程波方程偏移及其應用研究

    2016-11-16 00:55:52范廷恩黃旭日馬淑芳
    地球物理學報 2016年9期

    范廷恩, 黃旭日, 馬淑芳

    1 西南石油大學, 成都 610500 2 中海油研究總院, 北京 100028

    ?

    VSP時空域高角度單程波方程偏移及其應用研究

    范廷恩1,2, 黃旭日1, 馬淑芳2*

    1 西南石油大學, 成都 610500 2 中海油研究總院, 北京 100028

    偏移成像是VSP數(shù)據(jù)處理中的一個重要環(huán)節(jié),常規(guī)的VSP成像方法通常利用VSP-CDP轉換或Kirchhoff偏移,均存在保幅性差及成像精度低等問題,而波動方程疊前深度偏移被認為是對地下復雜構造進行成像的精確偏移方法.任意廣角波動方程作為一種高精度的空間域單程波波動方程,同時由于只含有二階偏導數(shù)項,易于數(shù)值實現(xiàn),與其他單程波波動方程相比,具有更大的成像傾角,因此是偏移成像的有力工具之一.本文將AWWE推廣應用到VSP數(shù)據(jù)成像中,實現(xiàn)了VSP時空域高角度單程波方程偏移.首先從三維標量任意廣角波動方程出發(fā),推導了完全匹配層吸收邊界條件,在基本不增加計算量的前提下有效地壓制了邊界反射成像噪音,同時利用非線性反演算法優(yōu)選參考速度來提高平方根算子的近似程度,從而提高高角度地層的成像精度.模型數(shù)值模擬實驗驗證了該方法的有效性,同時表明該方法在陡傾角構造情況下能取得很好的成像效果.最后對某地區(qū)實際觀測的VSP資料進行了偏移成像,并與地面地震偏移結果進行了對比,顯示出VSP波動方程偏移在成像分辨率上的優(yōu)勢.關鍵詞 任意廣角波動方程; VSP偏移; PML

    1 引言

    VSP 技術在井旁構造精細成像、井旁斷層識別、井旁地層巖性描述和地震波衰減等方面都取得了許多實際的應用效果.VSP成像是VSP 數(shù)據(jù)處理中的一個重要環(huán)節(jié),VSP成像方法主要包括常規(guī)的VSP-CDP轉換、基于射線理論的Kirchhoff積分法偏移和基于波動方程的偏移三種,前兩種方法都有其局限性,VSP-CDP轉換僅對水平層狀介質或橫向變速緩慢的介質才能取得比較好的處理結果,Kirchhoff 積分方法因它固有的高頻近似假設和射線理論在復雜介質中的缺陷(如焦散和多路徑等)使其在復雜構造區(qū)難以取得理想的結果.而波動方程偏移成像方法相對于射線追蹤技術,能更準確地描述波在復雜介質中的傳播以及波場能量的變化情況,這使得基于波動方程的偏移方法更受青睞.波動方程偏移方法又可以進一步分為基于單程波方程的偏移方法和基于雙程波方程的偏移方法.基于雙程波的逆時偏移不受地下界面傾角的限制,并且在保幅方面具有優(yōu)勢.孫文博和孫贊東(2010)開展了VSP資料偽譜法逆時偏移研究,李文杰等(2012)研究了二維各向同性介質中兩分量的非零偏VSP疊前逆時深度偏移方法,蔡曉慧等(2015)實現(xiàn)了基于自適應優(yōu)化有限差分方法的全波VSP逆時偏移,但該方法主要存在的應用瓶頸是:計算量大、存儲量大和噪聲問題.在波場的正向和逆時外推過程中需要存儲各個時刻的波場,非常耗時.同時將針對單程波方程提出的互相關成像條件用于逆時偏移,會產(chǎn)生低頻噪音干擾.因此三維逆時偏移的計算量和對速度模型的高要求使得它在目前工業(yè)界的應用效果并不是十分理想.

    Claerbout(1976)首次提出了時空域單程波方程,然而該方程具有明顯的傾角限制,無法對高陡構造準確成像,一些學者(馬在田,1983)采用不同的高階近似方法來提高成像傾角,但由于方程中往往包含高階偏導數(shù)項,而高階差分計算量明顯增加,使得數(shù)值計算難以實現(xiàn).Guddati(2006)提出了任意廣角波動方程(AWWE),它是一種在空間域表述的高精度單程波方程.AWWE的優(yōu)點是精度高,要提高AWWE的成像傾角,只需增加參考速度的個數(shù)即可,AWWE的形式簡單,只包含二階偏導數(shù)項,易于數(shù)值實現(xiàn),因此AWWE被廣泛應用于偏移成像,何兵壽等(2008a)用AWWE對地面地震資料做逆時偏移成像研究,同時將其與基于雙程聲波方程的逆時偏移方法做比較,指出基于AWWE的偏移方法不僅可以對陡傾角的地層準確成像,而且還可以很好地避免雙程波方程逆時偏移的低頻成像噪音;孫歧峰和杜啟振(2011)推導了頻率-空間域的AWWE,并利用有限差分高階分裂法求解方程,實現(xiàn)地震偏移成像,指出該方法能夠適應強橫向變速地區(qū)和陡傾角地層的成像.在吸收邊界條件方面,Heidari和Guddati(2006)導出廣角方程正向傳播的吸收邊界條件;何兵壽等(2008b)將加權校正及插值預測吸收邊界條件引入任意廣角方程的逆時偏移領域,這種邊界條件的本質是對Clayton 吸收邊界條件的改進,它難以有效地吸收地震波逆時延拓過程中投射到邊界的外行波;何兵壽等(2010)給出了適用于二維任意廣角單程聲波方程逆時延拓的吸收邊界條件,給出的算法只需在少量增加計算量的前提下有效地壓制了截斷邊界的偽反射,消除了邊界處的偏移噪聲;Chen等(2013)將完全匹配層(PerfectlyMatched Layer,PML)用于AWWE,有效地解決了AWWE數(shù)值計算中的吸收邊界問題.盡管AWWE具有高精度的特點,但是其計算過程包含重復的矩陣求逆和矩陣相乘,使得其計算代價比其他單程波方程高出許多,特別是增加參考速度個數(shù)時,其計算量顯著增加.周輝等(2014)改進了現(xiàn)有的有限差分計算方案,大幅度提高了計算效率,同時采用了一種有效的方法在f-k域比較準確地區(qū)分出倏逝波和非倏逝波,從而獲得一個合理的視速度作為f-k濾波器的門檻來壓制倏逝波干擾,提高成像質量.本文將AWWE推廣應用至VSP領域,實現(xiàn)了VSP資料的疊前深度偏移,采用非線性算法提高方程的成像精度,并應用了PML吸收邊界條件,很好地壓制了邊界反射成像噪音.

    2 基本原理

    2.1 高角度單程波VSP疊前深度偏移算法

    三維標量AWWE下行波方程的一般形式為(Guddati and Heidari,2005)

    (1)

    其中c為背景速度,dT=(1,0,…,0)n×1,u=(u,u1,u2,…,un-1)T,u1,u2…,un-1為輔助變量,T表示轉置,x、y、z為空間坐標,t表示時間變量,

    (2)

    (3)

    式中,c1,c2,…,cn為根據(jù)實際介質速度選取的參考速度.

    2.2 高角度單程波時域PML匹配方程

    在進行高角度單程波方程有限差分偏移時會遇到邊界反射問題,因此針對下行波式(1)推導其PML匹配方程,考慮沿x方向和y方向衰減,將其分裂成三項,即

    (4)

    (5)

    (6)

    返回時間域得到

    (7)

    將式(7)與式(4)的第三式寫在一起,即組成了三維任意廣角單程波下行波的時域PML匹配方程,即

    (8)

    由于是單程波方程不需要考慮Z方向的邊界反射問題,在計算時只需考慮前后左右四個邊界和四周四個棱邊,X方向兩個邊界只需要考慮X方向衰減,Y方向兩個邊界只需要考慮Y方向衰減,而四個棱邊需要同時考慮X和Y方向的衰減.

    2.3 非線性反演算法提高方程成像精度

    在推導高角度單程波方程時,不可避免地要對平方根算子做近似,對平方根算子的近似程度直接影響偏移成像的精度.從平方根算子表達的頻散關系出發(fā),對參考速度的選擇進行優(yōu)化.優(yōu)化的準則是在最小平方約束下使高角度方程對平方根算子的近似與準確的單程波平方根算子達到最佳的匹配.下面選定兩個參考速度來詳細說明高角度聲波方程參考速度的優(yōu)化過程.

    下行單程波平方根算子的準確表達式可以寫為

    (9)

    (10)

    其中θ表示波的傳播方向與界面法線之間的夾角.

    二維AWWE下行波方程:

    (11)

    為了評價AWWE的精度,對上式中的x、t做傅里葉變換,就可以得到頻率波數(shù)域中的高角度方程

    (12)

    進而可以整理出高角度方程對平方根算子的近似函數(shù)Fn(ξ,c1,…,cn),當選取兩個參考速度,并對速度進行歸一化后,得到F2(ξ,c1,c2)為(詳細推導過程參見附錄A):

    (13)

    (8)與式(10)的相對誤差為

    (14)

    在最小平方約束條件下,使得它們的相對誤差Er最小.相對誤差越小,平方根算子的近似程度越高,越有利于大角度地層準確成像.因此研究方程的精度對實際地震資料處理有一定的價值.可以通過非線性反演算法(牛頓法)求得Er最小時的參考速度(參見附錄B).

    在兩個參考速度的情況下,當固定第一個參考速度c1=1時,得到的優(yōu)化參考速度是c1=1,c2=5.24;當不固定參考速度c1的數(shù)值時,得到的優(yōu)化參考速度是c1=1.23,c2=6.46.在固定c1=1時的三個參考速度的情況下得到的優(yōu)化參考速度是c1=1,c2=2.82,c3=21.15.圖1是四種參考速度下AWWE頻散曲線,圖中OWWE代表的是下行單程波算子準確表達式的頻散曲線.圖1a是參考速度c1=1,c2=4的AWWE頻散曲線,AWWE的精度可以達到75°左右.圖1b是參考速度為c1=1,c2=5.24的AWWE頻散曲線,在80°以下它與準確的頻散曲線基本上是符合的,與圖1a相比,圖1b的最大準確成像角度由75°提高到80°.圖1c是參考速度為c1=1.23,c2=6.46的AWWE頻散曲線,與圖1a和圖1b相比,其精度可以達到82°左右,圖1d是參考速度為c1=1,c2=2.82,c3=21.15的AWWE頻散曲線,最大準確成像角度可接近90°,同時可以看到AWWE的頻散曲線與準確的頻散曲線整體上都可以精確地吻合,成像精度大大得到了提高.

    從上面的分析中可以看出:優(yōu)化的三個參考速度AWWE頻散曲線的精度最高;優(yōu)化的兩個參考速度AWWE頻散曲線的最大準確成像角度都得到了提高,但是會損失小角度處的精度;與Guddati給出的c1=1,c2=4參考速度相比,優(yōu)化方法選取的參考速度可使大角度處AWWE的精度得到提高.采用優(yōu)化參考速度參數(shù)的AWWE做疊前深度偏移處理,而非采用增加參考速度個數(shù)的方法,就可以在避免增加運算量的情況下進一步解決大角度成像的問題.

    圖1 不同參考速度對應的AWWE頻散曲線及其誤差(a) c1=1,c2=4; (b) c1=1,c2=5.24; (c) c1=1.23,c2=6.46; (d) c1=1,c2=2.82,c3=21.15,圖中OWWE代表的是下行單程波算子準確表達式的頻散曲線.Fig.1 The dispersion relations and errors with different reference velocities

    3 數(shù)值模擬實驗

    3.1 層狀介質模型

    用一個四層介質模型(如圖2a所示)的疊前深度偏移實驗分析高角度單程波VSP數(shù)據(jù)的偏移效果.模型大小為1800 m×3000 m,界面分別在900 m、1800 m和2700 m深度處,四個水平層速度分別為2000 m·s-1、2200 m·s-1、2400 m·s-1和2600 m·s-1.空間網(wǎng)格大小為3 m,時間采樣間隔為0.5 ms,合成地震記錄長度為2 s.炮點深度為3 m,水平方向600個網(wǎng)格點上均勻放置600炮激發(fā).井水平坐標903 m,從深度為300 m到深度1200 m,以15 m為間隔均勻布置60個檢波點.

    我們利用二維聲波方程正演得到原始的共炮集數(shù)據(jù)之后,由于直達波(下行波)是先于反射波(上行波)到達檢波器的,所以我們可以在共炮道集上對每個道沿著時間增大的方向尋值,根據(jù)直達波信號的大小,當尋到第一個與直達波信號量級相當?shù)臄?shù)值時,沿著時間增大方向將一個子波長度內(nèi)的值變?yōu)榱悖@樣就可以去除直達波.去除了直達波之后抽取共檢波點道集,再利用共檢波點道集數(shù)據(jù)通過高角度單程波疊前深度偏移得到60個共檢波點道集偏移結果,通過疊加得到最后的結果如圖2b所示,可以看出三個水平界面得到準確成像,偏移噪音很小,但是成像范圍并不是很大,這是由VSP本身的觀測范圍所限制的.

    3.2 復雜模型

    我們用復雜模型(如圖3a所示)的疊前深度偏移實驗分析高角度單程波VSP數(shù)據(jù)的偏移效果.模型大小為3000 m×3000 m,水平方向網(wǎng)格大小為5 m,深度方向網(wǎng)格大小為5 m,時間采樣間隔為0.5 ms,合成地震記錄長度為2 s.從上到下各地層中聲波速度分別為4000 m·s-1、4100 m·s-1、4400 m·s-1、4100 m·s-1、4000 m·s-1、4200 m·s-1和4300 m·s-1.炮點深度為5 m,水平方向600個網(wǎng)格點上均勻放置600炮激發(fā).井水平坐標2000 m,從深度為250 m到深度850 m,以20 m為間隔均勻布置30個檢波點.

    圖3b是多炮偏移疊加成像的結果,可以看出,對于離井較近的高陡地層得到了很好地成像,偏移的成像噪音也得到了很好的壓制.

    圖2 簡單模型偏移結果(a) 水平層狀模型; (b) 偏移結果.Fig.2 The migration result of simple model(a) Horizontal layered model; (b) Its migration result.

    圖3 復雜模型偏移結果(a) 復雜模型; (b) 偏移結果.Fig.3 The migration result of complex model(a) Complex model; (b) Its migration result.

    4 實際資料測試

    選取Q油田Walk-away VSP測線進行了疊前深度偏移成像,本次Walk-away VSP 觀測系統(tǒng)震源布設范圍為-1900~1900 m,炮間距25 m,檢波器布設范圍為600~800 m,采用間距為10 m的20級檢波器接收,每個“排列”上提5 m重復測量,最終達到5 m的接收間距,采集時一共布置了八條測線,這里對其中一條測線進行處理.圖4為偏移速度模型,由于只是對其中一條測線的處理,所以這里給出的只是截取了的一部分速度模型,并將坐標值均轉化為正值,以便于描述.抽取共成像點道集進行偏移成像處理,同時為了更好地進行同相疊加,抽取了40個共檢波點道集不同偏移距的共成像點道集,利用互相關方法對得到的共成像點道集進行整道校正,從而提高成像質量.圖5為水平偏移距400 m校正前后的共成像點道集,對比可以看出,校正后的共成像點道集同相軸連續(xù)性得到增強,圖6為共成像點校正后的偏移結果,可以看出同相軸連續(xù)性良好,分辨率也很高,這主要是得益于VSP數(shù)據(jù)具有較高的主頻和信噪比.

    圖4 偏移速度模型Fig.4 Velocity model for migration

    圖5 共成像點道集校正(偏移距400 m處)(a) 校正前; (b) 校正后.Fig.5 The correction of CIP gathers (offset=400 m) (a) Before; (b) After.

    為了進一步分析成像結果,將地面地震偏移成像結果與此次VSP偏移結果進行對比,如圖7,可以看出VSP偏移結果與地面地震剖面產(chǎn)狀吻合較好,并具有更高的分辨率和保幅性,主頻達到100 Hz,頻帶寬度達到8~180 Hz,為后續(xù)波組特征的識別和斷層的精細刻畫提供了高品質資料基礎.

    圖6 共成像點校正后的VSP偏移結果

    圖7 VSP與地面地震偏移結果的對比(a) VSP剖面與地面地震剖面疊合圖; (b) 地面地震剖面頻譜; (c) VSP剖面頻譜.Fig.7 The comparison of VSP migration result and surface seismic section(a) The overlapping graph of VSP migration result and surface seismic section;(b) Spectrum of surface seismic section; (c) Spectrum of VSP migration result.

    5 結論

    與其他單程波波動方程相比,AWWE是一種高精度的單程波方程,可以實現(xiàn)高陡構造的準確成像.本文將這種時空域高角度單程波方程偏移推廣應用到VSP數(shù)據(jù)中,從三維標量任意廣角波動方程出發(fā),推導了PML公式,很好地壓制了邊界反射成像噪音,同時應用牛頓法在最小平方意義下求解方程組來選取參考速度,從而提高了地層成像傾角,優(yōu)化后成像傾角接近90°.模型試算結果表明該方法在陡傾角情況下能取得較好的成像效果,實際資料的試處理也表明VSP時空域高角度單程波方程偏移成像結果相比地面地震偏移剖面具有更高的分辨率.

    致謝 感謝中海油研究總院地球物理總師李緒宣、開發(fā)研究院院長胡光義的指導;感謝中海油天津分公司、中海油服提供的VSP資料;感謝重大專項項目組及聯(lián)合單位中國石油大學(北京)提供的幫助.

    附錄A 提高精度的理論分析

    二維時空域下行波AWWE為

    (A1)

    AWWE作為一種單程波方程(OWWE),不可避免地要對平方根算子做近似,其近似程度直接影響偏移成像的精度和所能準確成像的界面傾角范圍.用水平波數(shù)kx,角頻率ω表示的準確的下行單程波方程表達形式如下:

    (A2)

    令ckx/ω=sinθ,θ表示波的傳播方向與界面法線之間的夾角,θ的范圍是0到90°,所以(A2)可以簡化為

    (A3)

    所以我們對方程(A1)中的x,t做傅里葉變換,就可以得到

    (A5)

    移項可得

    (A6)

    (A7)

    (A8)

    附錄B 非線性反演算法選取參考速度

    相對誤差Er最小時,平方根算子的近似函數(shù)Fn(ξ,c1,…,cn)精度最高.Er最小時的參考速度即為我們要求取的優(yōu)化參考速度.

    (B1)

    定義目標函數(shù)

    (B2)

    其中θ1和θ2是所選擇的角度范圍,通常選取θ1=0°,θ2=90°.目標函數(shù)Ir表示的是在一定角度范圍內(nèi)對誤差Er的最小平方約束.目標函數(shù)越小,則F2與R越接近,此時AWWE的成像精度越高.

    通過求取目標函數(shù)Ir的最小值,得到參考速度的值.同樣以兩個參考速度的情況為例,在最小平方意義下由方程(B2)可得

    (B3)

    (B4)

    (B5)

    (B6)

    (B7)

    (B8)

    (B9)

    (B10)

    把(B1)式代入到(B6)式至(B10)式中,再根據(jù)(B5)式就可以得到每一次的參考速度更正值.

    當參考速度的更正值小于給定的值(本文選取10-6)時,則終止迭代,并輸出此時的參考速度.輸出的參考速度可以使AWWE的精度在所選擇的角度范圍內(nèi)達到全局最優(yōu).也就是我們要求的優(yōu)化參考速度.

    Cai X H, Liu Y, Wang J M, et al. 2015. Full-wavefield VSP reverse-time migration based on the adaptive optimal finite-difference scheme.ChineseJ.Geophys.(in Chinese), 58(9): 3317-3334, doi: 10.6038/cjg20150925.

    Chen H M, Zhou H, Lin H, et al. 2013. Application of perfectly matched layer for scalar arbitrarily wide-angle wave equations.Geophysics, 78(1): T29-T39.

    Claerbout J F. 1976. Fundamentals of Geophysical Data Processing. Oxford: Scientific Publications.

    Guddati M N, Heidari A H. 2005. Migration with arbitrarily wide-angle wave equations.Geophysics, 70(3): S61-S70.

    Guddati M N. 2006. Arbitrarily wide-angle wave equations for complex media.ComputerMethodsinAppliedMechanicsandEngineering, 195(1-3): 65-93.

    He B S, Zhang H X, Zhang J. 2008a.Prestack reverse-time depth migration of arbitrarily wide-angle wave equations.ActaSeismologicaSinica(in Chinese), 30(5): 491-499. He B S, Zhang H X, Han L H, et al. 2008b.Prestack reverse-time depth migration of two acoustic wave equation and comparison between both.OilGeophysicalProspecting(in Chinese), 43(6): 661-684.

    He B S, Zhang H X, Zhang J J. 2010.Absorbing boundary conditions of arbitrarily wide-angle acoustic equations reverse-time migration.JournalofChinaCoalSociety(in Chinese), 35(1): 106-109.

    Heidari A H, Guddati M N. 2006.Highly accurate absorbing boundary conditions for wide-angle wave equations.Geophysics, 71(3): S85-S97.

    Li W J, Qu S L, Wei X C, et al. 2012. Discussion about reverse-time depth migration technology of prestack elastic wavefield for offset VSP data.ChineseJ.Geophys.(in Chinese), 55(1): 238-251, doi: 10.6038/j.issn.0001-5733.2012.01.023.

    Ma Z T. 1983. A splitting-up method for solution of higher-order migration equation by finite-difference scheme.ActaGeophysicaSinica(in Chinese), 26(4): 377-388.

    Sun Q F, Du Q Z. 2011.The frequency-space domain prestack depth migration with arbitrarily wide-angle wave equation.OilGeophysicalProspecting(in Chinese), 46(6): 890-896. Sun W B, Sun Z D. 2010. VSP reverse time migration based on the pseudo-spectral method and its applications.ChineseJ.Geophys.(in Chinese), 53(9): 2196-2203, doi: 10.3969/j.issn.0001-5733.2010.09.020. Zhou H, Chen H M, Li X X, et al. 2014. Methods for improving computational efficiency of arbitrarily wide-angle wave equation and suppressing evanescent wave noises.ChineseJ.Geophys. (in Chinese), 57(4): 1188-1198, doi: 10.6038/cjg20140416.

    附中文參考文獻

    蔡曉慧, 劉洋, 王建民等. 2015. 基于自適應優(yōu)化有限差分方法的全波VSP逆時偏移. 地球物理學報, 58(9): 3317-3334, doi: 10.6038/cjg20150925.

    何兵壽, 張會星, 張晶. 2008a. 任意廣角波動方程疊前逆時深度偏移.地震學報, 30(5): 491-499.

    何兵壽, 張會星, 韓令賀等. 2008b. 兩種聲波方程的疊前逆時深度偏移及比較.石油地球物理勘探, 43(6): 661-684.

    何兵壽, 張會星, 張建軍. 2010. 任意廣角聲波方程逆時偏移的吸收邊界條件. 煤炭學報, 35(1): 106-109.

    李文杰, 曲壽利, 魏修成等. 2012. 非零偏VSP彈性波疊前逆時深度偏移技術探討. 地球物理學報, 55(1): 238-251, doi: 10.6038/j.issn.0001-5733.2012.01.023.

    馬在田. 1983. 高階方程偏移的分裂算法. 地球物理學報, 26(4): 377-388.

    孫歧峰, 杜啟振. 2011. 任意廣角波動方程頻率—空間域疊前深度偏移成像. 石油地球物理勘探, 46(6): 890-896.

    孫文博, 孫贊東. 2010. 基于偽譜法的VSP逆時偏移及其應用研究. 地球物理學報, 53(9): 2196-2203, doi: 10.3969/j.issn.0001-5733.2010.09.020.

    周輝, 陳漢明, 李緒宣等. 2014. 提高任意廣角波動方程計算效率及壓制倏逝波干擾方法. 地球物理學報, 57(4): 1188-1198, doi: 10.6038/cjg20140416.

    (本文編輯 胡素芳)

    VSP migration based on high angle one-way wave equation in time-space domain and its applications

    FAN Ting-En1,2, HUANG Xu-Ri1, MA Shu-Fang2*

    1SouthwestPetroleumUniversity,Chengdu610500,China2CNOOCResearchInstitute,Beijing100028,China

    Migration is a key step in VSP data processing, VSP imaging commonly uses VSP-CDP conversion or Kirchhoff migration which has issues of poor amplitude conservation and low imaging precision, while wave equation prestack depth migration is an effective method to image complex structures.Arbitrarily wide-angle wave equation (AWWE) is a space domain one-way wave equation (OWWE) with a high accuracy and easy numerical implementation because of only including two-order partial derivative term.Compared with other OWWEs, AWWE is capable of imaging steeper dips, which enables it to be one of powerful migration tools. In this paper, we extended AWWE to VSP imaging and realized VSP migration based on high angle one-way wave equation in time-space domain. And the perfectly matched layer (PML) absorbing boundary condition was derived to provide a good approach to suppress artifacts arising at truncation boundaries at a price of limited additional computation cost. Meanwhile, we optimized reference velocities using nonlinear algorithm to enhance the approximation degree of square root operator and improved imaging precision of high-dip layers. Model tests proved the effectivity of the method and showed that it can obtain good imaging result in the case of steep-dip structures. The method was also applied to field VSP data, the comparison of VSP migration result and surface seismic section showed that VSP migration provides a higher resolution image.

    Arbitrarily wide-angle wave equation; VSP migration; Perfectly matched layer

    10.6038/cjg20160928.

    國家科技重大專項(2011ZX05024)資助.

    范廷恩,男,1972年生,高級工程師,主要從事地震反演、地震屬性分析及油藏精細描述等方面的研究. E-mail:fante@cnooc.com.cn*通訊作者 馬淑芳,女,1984年生,博士,主要從事井中地震及井地聯(lián)合高分辨率處理等方面的研究. E-mail:mashf2@cnooc.com.cn

    10.6038/cjg20160928

    P631

    2016-04-12,2016-07-06收修定稿

    范廷恩, 黃旭日, 馬淑芳. 2016. VSP時空域高角度單程波方程偏移及其應用研究. 地球物理學報,59(9):3459-3469,

    Fan T E, Huang X R, Ma S F. 2016. VSP migration based on high angle one-way wave equation in time-space domain and its applications.ChineseJ.Geophys. (in Chinese),59(9):3459-3469,doi:10.6038/cjg20160928.

    一二三四社区在线视频社区8| 欧美日韩一级在线毛片| 免费久久久久久久精品成人欧美视频| 精品高清国产在线一区| 亚洲中文字幕日韩| 国产91精品成人一区二区三区| 长腿黑丝高跟| 午夜免费激情av| 女人精品久久久久毛片| 丰满饥渴人妻一区二区三| 无遮挡黄片免费观看| 日韩欧美一区二区三区在线观看| 一级片免费观看大全| 99久久精品国产亚洲精品| 我的亚洲天堂| 不卡一级毛片| 亚洲国产精品合色在线| 午夜日韩欧美国产| 日本黄色日本黄色录像| 少妇被粗大的猛进出69影院| 亚洲精品粉嫩美女一区| 99香蕉大伊视频| 国产片内射在线| 久久精品亚洲精品国产色婷小说| 天天躁夜夜躁狠狠躁躁| 国产一区二区激情短视频| www.www免费av| 亚洲免费av在线视频| 脱女人内裤的视频| 久久婷婷成人综合色麻豆| 黄色女人牲交| 这个男人来自地球电影免费观看| 亚洲第一欧美日韩一区二区三区| 美国免费a级毛片| 免费在线观看影片大全网站| 99国产精品免费福利视频| 男女之事视频高清在线观看| 亚洲成国产人片在线观看| 国产成人精品久久二区二区免费| 日本免费一区二区三区高清不卡 | videosex国产| 99久久国产精品久久久| 看免费av毛片| 亚洲成人久久性| 亚洲激情在线av| 精品福利永久在线观看| 老司机福利观看| 午夜日韩欧美国产| 女人高潮潮喷娇喘18禁视频| 日韩视频一区二区在线观看| 黄色片一级片一级黄色片| 看免费av毛片| 国产男靠女视频免费网站| 国产黄a三级三级三级人| 午夜精品在线福利| 一区二区三区国产精品乱码| 无遮挡黄片免费观看| 欧美日韩乱码在线| 成人18禁在线播放| 国产精品亚洲一级av第二区| 国产高清激情床上av| 18禁美女被吸乳视频| 两人在一起打扑克的视频| 日日摸夜夜添夜夜添小说| 久久久久久免费高清国产稀缺| 男女下面进入的视频免费午夜 | 91成人精品电影| 久久久久久大精品| 日本撒尿小便嘘嘘汇集6| 两个人看的免费小视频| 国产精品免费视频内射| 免费av毛片视频| 免费人成视频x8x8入口观看| www.熟女人妻精品国产| 国产精品亚洲一级av第二区| 长腿黑丝高跟| 操美女的视频在线观看| 精品久久久精品久久久| 国产av在哪里看| 天天添夜夜摸| 一级片'在线观看视频| 美女国产高潮福利片在线看| 人人妻人人爽人人添夜夜欢视频| 久99久视频精品免费| 亚洲精品av麻豆狂野| 亚洲一区二区三区色噜噜 | √禁漫天堂资源中文www| www日本在线高清视频| 天天影视国产精品| 日韩人妻精品一区2区三区| 可以免费在线观看a视频的电影网站| 一本综合久久免费| 国产精品香港三级国产av潘金莲| 亚洲av熟女| 狂野欧美激情性xxxx| 亚洲一区二区三区不卡视频| 久久狼人影院| 亚洲av熟女| av网站在线播放免费| 亚洲午夜精品一区,二区,三区| 久久狼人影院| 十八禁人妻一区二区| xxx96com| 日本五十路高清| 性色av乱码一区二区三区2| 可以免费在线观看a视频的电影网站| 色婷婷久久久亚洲欧美| 日韩欧美一区二区三区在线观看| 国产欧美日韩一区二区精品| 免费在线观看日本一区| 亚洲精品在线观看二区| 久久久久久久精品吃奶| 人妻丰满熟妇av一区二区三区| 黄色 视频免费看| 国产三级黄色录像| 中文字幕色久视频| 亚洲熟女毛片儿| 长腿黑丝高跟| 亚洲人成伊人成综合网2020| 国产免费男女视频| 丰满的人妻完整版| av福利片在线| 热99国产精品久久久久久7| 男女下面进入的视频免费午夜 | 午夜免费观看网址| 69精品国产乱码久久久| 精品国产国语对白av| 国产男靠女视频免费网站| 免费观看人在逋| 中出人妻视频一区二区| 侵犯人妻中文字幕一二三四区| 国产av一区二区精品久久| 亚洲国产欧美一区二区综合| 午夜免费鲁丝| 亚洲av成人不卡在线观看播放网| 成在线人永久免费视频| 亚洲精品一区av在线观看| 国产精品自产拍在线观看55亚洲| 老司机午夜福利在线观看视频| 又黄又爽又免费观看的视频| 99国产精品免费福利视频| 极品人妻少妇av视频| 欧美国产精品va在线观看不卡| 久久精品国产亚洲av香蕉五月| 国产精品久久久久久人妻精品电影| 久久久精品国产亚洲av高清涩受| 精品久久蜜臀av无| 国产日韩一区二区三区精品不卡| 欧美日韩瑟瑟在线播放| 国产一区二区三区视频了| 日韩 欧美 亚洲 中文字幕| 天天添夜夜摸| 丝袜人妻中文字幕| 首页视频小说图片口味搜索| 一级毛片精品| 日韩中文字幕欧美一区二区| 在线视频色国产色| 久久久久久大精品| 成人手机av| 色综合欧美亚洲国产小说| 窝窝影院91人妻| 变态另类成人亚洲欧美熟女 | 精品一品国产午夜福利视频| 精品免费久久久久久久清纯| 91老司机精品| 天天影视国产精品| 欧美日韩一级在线毛片| 超碰成人久久| 在线天堂中文资源库| 免费在线观看亚洲国产| 亚洲专区国产一区二区| 日韩视频一区二区在线观看| 精品久久久久久久久久免费视频 | 久9热在线精品视频| 精品卡一卡二卡四卡免费| 亚洲国产看品久久| bbb黄色大片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲视频免费观看视频| 国产国语露脸激情在线看| 超碰成人久久| 欧美黑人欧美精品刺激| aaaaa片日本免费| 女性被躁到高潮视频| 1024香蕉在线观看| 国产精品亚洲av一区麻豆| 777久久人妻少妇嫩草av网站| 亚洲男人天堂网一区| 精品无人区乱码1区二区| 欧美日韩瑟瑟在线播放| 宅男免费午夜| 免费在线观看亚洲国产| 国产无遮挡羞羞视频在线观看| 亚洲狠狠婷婷综合久久图片| 欧美大码av| 老熟妇仑乱视频hdxx| 久久中文字幕人妻熟女| 久久久久久久久免费视频了| 新久久久久国产一级毛片| 久久亚洲真实| av电影中文网址| 国产精品久久久人人做人人爽| 亚洲一卡2卡3卡4卡5卡精品中文| 九色亚洲精品在线播放| 少妇裸体淫交视频免费看高清 | 婷婷六月久久综合丁香| 精品一区二区三区四区五区乱码| 麻豆成人av在线观看| 国产一区二区在线av高清观看| 女人精品久久久久毛片| 亚洲精品在线美女| 国产精品自产拍在线观看55亚洲| 我的亚洲天堂| 美女午夜性视频免费| 如日韩欧美国产精品一区二区三区| 高清欧美精品videossex| 在线永久观看黄色视频| 女人精品久久久久毛片| 12—13女人毛片做爰片一| 国内毛片毛片毛片毛片毛片| av国产精品久久久久影院| 国产精品一区二区在线不卡| 精品国产乱子伦一区二区三区| 久久久国产欧美日韩av| 嫁个100分男人电影在线观看| 国产成人精品久久二区二区免费| 在线永久观看黄色视频| 精品日产1卡2卡| 日本免费一区二区三区高清不卡 | 性少妇av在线| 亚洲 欧美 日韩 在线 免费| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 1024香蕉在线观看| 中文字幕色久视频| 人人妻,人人澡人人爽秒播| 中国美女看黄片| 亚洲欧美一区二区三区久久| 免费女性裸体啪啪无遮挡网站| 日本三级黄在线观看| 91字幕亚洲| 国产精品久久视频播放| 欧美国产精品va在线观看不卡| 国产午夜精品久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品999在线| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲三区欧美一区| 人成视频在线观看免费观看| 精品久久久精品久久久| 91成年电影在线观看| 久久影院123| 欧美在线黄色| 桃色一区二区三区在线观看| 国产一区二区三区视频了| 交换朋友夫妻互换小说| 久久精品国产综合久久久| 精品乱码久久久久久99久播| 精品午夜福利视频在线观看一区| 亚洲精品在线观看二区| 国产国语露脸激情在线看| 国产精品亚洲一级av第二区| 精品福利永久在线观看| 亚洲一区二区三区不卡视频| 午夜福利在线免费观看网站| 国产精品电影一区二区三区| 深夜精品福利| 久久国产精品人妻蜜桃| 男女做爰动态图高潮gif福利片 | 最新美女视频免费是黄的| 在线视频色国产色| 人人妻人人澡人人看| 亚洲熟女毛片儿| 高清黄色对白视频在线免费看| 看黄色毛片网站| 日本精品一区二区三区蜜桃| 久久热在线av| 国产一卡二卡三卡精品| av网站在线播放免费| 国产av在哪里看| 国产熟女午夜一区二区三区| 精品一区二区三区视频在线观看免费 | 欧美成狂野欧美在线观看| 国产精品乱码一区二三区的特点 | 女警被强在线播放| 久久婷婷成人综合色麻豆| 夜夜爽天天搞| 又黄又爽又免费观看的视频| 日日干狠狠操夜夜爽| netflix在线观看网站| 国产精品二区激情视频| 久久影院123| 国产真人三级小视频在线观看| 老司机午夜十八禁免费视频| 久久婷婷成人综合色麻豆| 日韩人妻精品一区2区三区| 又黄又粗又硬又大视频| 五月开心婷婷网| 国产日韩一区二区三区精品不卡| 国产一卡二卡三卡精品| 国产精品永久免费网站| av视频免费观看在线观看| 免费在线观看亚洲国产| 亚洲伊人色综图| 成人黄色视频免费在线看| 一边摸一边抽搐一进一出视频| 天天添夜夜摸| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久 | 亚洲五月婷婷丁香| 午夜两性在线视频| 午夜福利在线观看吧| 99国产精品一区二区三区| 成人免费观看视频高清| 日日爽夜夜爽网站| 99在线人妻在线中文字幕| 日本三级黄在线观看| 亚洲第一av免费看| 欧美成人免费av一区二区三区| 成人av一区二区三区在线看| 久久午夜亚洲精品久久| 国产视频一区二区在线看| av网站免费在线观看视频| 日韩免费高清中文字幕av| 欧美日本中文国产一区发布| 亚洲精品久久午夜乱码| 美女高潮喷水抽搐中文字幕| 国产成人av激情在线播放| 午夜免费观看网址| 亚洲精品国产色婷婷电影| 男人舔女人下体高潮全视频| 精品日产1卡2卡| 最好的美女福利视频网| 在线av久久热| 精品久久久久久,| 免费在线观看黄色视频的| 午夜两性在线视频| 国产色视频综合| 午夜两性在线视频| 国产精华一区二区三区| 激情视频va一区二区三区| www.www免费av| 亚洲一区中文字幕在线| 欧美日韩乱码在线| 香蕉久久夜色| 成在线人永久免费视频| 成人永久免费在线观看视频| 露出奶头的视频| 欧美久久黑人一区二区| 国产一卡二卡三卡精品| 最新美女视频免费是黄的| 日韩欧美一区二区三区在线观看| 男女高潮啪啪啪动态图| 欧美老熟妇乱子伦牲交| 成人18禁在线播放| 久久精品国产亚洲av香蕉五月| 精品欧美一区二区三区在线| 国产精品久久久av美女十八| 婷婷丁香在线五月| 欧美黑人精品巨大| 国产精品永久免费网站| 亚洲,欧美精品.| 成人18禁在线播放| 久久亚洲精品不卡| 高清av免费在线| 欧美在线黄色| 黑人猛操日本美女一级片| av在线播放免费不卡| 大码成人一级视频| 99re在线观看精品视频| 老司机福利观看| 男人舔女人下体高潮全视频| 欧洲精品卡2卡3卡4卡5卡区| 在线永久观看黄色视频| 久久人人爽av亚洲精品天堂| www.999成人在线观看| 免费在线观看日本一区| 琪琪午夜伦伦电影理论片6080| 国产激情欧美一区二区| 免费高清在线观看日韩| 久久精品国产综合久久久| 成在线人永久免费视频| 在线十欧美十亚洲十日本专区| 亚洲国产精品一区二区三区在线| 亚洲国产精品合色在线| 18美女黄网站色大片免费观看| 欧美不卡视频在线免费观看 | 大型黄色视频在线免费观看| 午夜老司机福利片| 亚洲精品一区av在线观看| 亚洲一区高清亚洲精品| 成人三级黄色视频| 高清av免费在线| 巨乳人妻的诱惑在线观看| 成人三级做爰电影| 久久精品国产亚洲av香蕉五月| 国产高清videossex| 他把我摸到了高潮在线观看| 高潮久久久久久久久久久不卡| 国产精品二区激情视频| 日日夜夜操网爽| 国产精品日韩av在线免费观看 | 麻豆久久精品国产亚洲av | 最好的美女福利视频网| aaaaa片日本免费| 色在线成人网| 亚洲成国产人片在线观看| 亚洲国产精品999在线| 成年女人毛片免费观看观看9| 精品一区二区三卡| 搡老乐熟女国产| 亚洲av美国av| 日日爽夜夜爽网站| 亚洲狠狠婷婷综合久久图片| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 成人特级黄色片久久久久久久| tocl精华| 97碰自拍视频| 久99久视频精品免费| 一二三四在线观看免费中文在| 99精国产麻豆久久婷婷| 人妻丰满熟妇av一区二区三区| 精品人妻在线不人妻| 99久久国产精品久久久| 动漫黄色视频在线观看| 国产一区二区三区在线臀色熟女 | 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 免费在线观看影片大全网站| 久久久久九九精品影院| 美女国产高潮福利片在线看| 国产av精品麻豆| 高潮久久久久久久久久久不卡| 久久精品国产亚洲av香蕉五月| 亚洲免费av在线视频| 黑人巨大精品欧美一区二区mp4| 久久久精品欧美日韩精品| svipshipincom国产片| 99久久人妻综合| 久久精品亚洲av国产电影网| www日本在线高清视频| 国产精品爽爽va在线观看网站 | 国产高清激情床上av| 十分钟在线观看高清视频www| 精品国产亚洲在线| 日日干狠狠操夜夜爽| 男男h啪啪无遮挡| 91麻豆精品激情在线观看国产 | 久久午夜综合久久蜜桃| 一级a爱视频在线免费观看| 国产精品av久久久久免费| 亚洲在线自拍视频| av欧美777| 国产单亲对白刺激| 久久人人97超碰香蕉20202| 国产成人精品无人区| 中文字幕高清在线视频| 亚洲熟妇中文字幕五十中出 | 亚洲成人精品中文字幕电影 | 黄色毛片三级朝国网站| 国产精品免费一区二区三区在线| av网站在线播放免费| 亚洲少妇的诱惑av| 女警被强在线播放| 久久精品91蜜桃| 在线观看午夜福利视频| 天堂俺去俺来也www色官网| tocl精华| 男女午夜视频在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产91精品成人一区二区三区| 丁香六月欧美| 久久久精品欧美日韩精品| 欧美国产精品va在线观看不卡| www.www免费av| 12—13女人毛片做爰片一| 日本vs欧美在线观看视频| 村上凉子中文字幕在线| 国产激情久久老熟女| 88av欧美| 免费一级毛片在线播放高清视频 | 亚洲一区二区三区欧美精品| 国产精品亚洲一级av第二区| 操出白浆在线播放| 美女大奶头视频| 最新美女视频免费是黄的| 国产免费av片在线观看野外av| 女人精品久久久久毛片| 欧美亚洲日本最大视频资源| 亚洲精品中文字幕一二三四区| 精品一区二区三区av网在线观看| 国产成人系列免费观看| 国产亚洲欧美98| 欧美日韩亚洲高清精品| 日韩大码丰满熟妇| 热99国产精品久久久久久7| 精品国产乱码久久久久久男人| 国产成人av激情在线播放| av在线播放免费不卡| 校园春色视频在线观看| 他把我摸到了高潮在线观看| 国产成人免费无遮挡视频| 日日干狠狠操夜夜爽| 亚洲五月色婷婷综合| 又大又爽又粗| 韩国精品一区二区三区| 成人18禁在线播放| 国产三级黄色录像| 精品福利观看| 两性夫妻黄色片| 亚洲少妇的诱惑av| 国内久久婷婷六月综合欲色啪| 欧美日本中文国产一区发布| 亚洲一区二区三区不卡视频| 最近最新免费中文字幕在线| 老司机靠b影院| 亚洲第一欧美日韩一区二区三区| 亚洲成人国产一区在线观看| 真人做人爱边吃奶动态| cao死你这个sao货| 黄色 视频免费看| 成人国语在线视频| 免费av毛片视频| 麻豆一二三区av精品| 免费搜索国产男女视频| 老汉色∧v一级毛片| 欧美丝袜亚洲另类 | 天堂中文最新版在线下载| 亚洲成人免费电影在线观看| 日韩一卡2卡3卡4卡2021年| 丰满的人妻完整版| 中文字幕高清在线视频| 中国美女看黄片| 午夜免费成人在线视频| 国产aⅴ精品一区二区三区波| 日韩视频一区二区在线观看| 一二三四在线观看免费中文在| 无人区码免费观看不卡| 日本 av在线| 一级毛片高清免费大全| 老鸭窝网址在线观看| 嫩草影视91久久| 美女扒开内裤让男人捅视频| 99国产精品免费福利视频| 国产精品乱码一区二三区的特点 | 淫妇啪啪啪对白视频| 中文字幕人妻丝袜制服| 麻豆久久精品国产亚洲av | 岛国视频午夜一区免费看| 色播在线永久视频| 精品国产美女av久久久久小说| 脱女人内裤的视频| 最近最新中文字幕大全免费视频| 少妇粗大呻吟视频| av网站在线播放免费| 久热爱精品视频在线9| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区三| 99国产精品99久久久久| 久久天躁狠狠躁夜夜2o2o| 午夜免费成人在线视频| 国产aⅴ精品一区二区三区波| 91av网站免费观看| 侵犯人妻中文字幕一二三四区| 久久久精品国产亚洲av高清涩受| av视频免费观看在线观看| 午夜91福利影院| 电影成人av| 国产深夜福利视频在线观看| 制服人妻中文乱码| 欧美日韩乱码在线| 岛国在线观看网站| 最近最新免费中文字幕在线| 亚洲av五月六月丁香网| 亚洲中文av在线| 日本wwww免费看| 久久欧美精品欧美久久欧美| a级毛片在线看网站| 热re99久久国产66热| av在线天堂中文字幕 | 精品国产乱码久久久久久男人| 一区二区三区激情视频| 少妇的丰满在线观看| 欧美激情极品国产一区二区三区| 97人妻天天添夜夜摸| 亚洲成a人片在线一区二区| 国产精品二区激情视频| 9热在线视频观看99| 国产成人一区二区三区免费视频网站| 日韩欧美一区二区三区在线观看| 国产伦一二天堂av在线观看| 久9热在线精品视频| 久久精品国产亚洲av香蕉五月| 女人爽到高潮嗷嗷叫在线视频| 男人操女人黄网站| 婷婷六月久久综合丁香| 欧美日韩国产mv在线观看视频| 国产成人影院久久av| 精品久久久久久,| 丁香六月欧美| 亚洲一区二区三区欧美精品| 亚洲国产精品999在线| 国产1区2区3区精品| 视频在线观看一区二区三区| 91成人精品电影| 日日夜夜操网爽| www.www免费av| 啦啦啦免费观看视频1| 一级毛片女人18水好多| 亚洲一区二区三区欧美精品| 在线观看一区二区三区激情| 日日爽夜夜爽网站| 久久久久九九精品影院| 免费少妇av软件|