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

    地震波逆時偏移中不同域共成像點道集偏移噪聲分析

    2014-02-10 03:04:16陳可洋吳清嶺范興才陳樹民李來林劉振寬
    巖性油氣藏 2014年2期
    關(guān)鍵詞:波場角度界面

    陳可洋,吳清嶺,范興才,陳樹民,李來林,劉振寬

    地震波逆時偏移中不同域共成像點道集偏移噪聲分析

    陳可洋1,吳清嶺1,范興才1,陳樹民1,李來林1,劉振寬2

    (1.中國石油大慶油田有限責(zé)任公司勘探開發(fā)研究院,黑龍江大慶163712;2.中國石油大慶油田有限責(zé)任公司勘探事業(yè)部,黑龍江大慶163453)

    地震波逆時成像方法通常輸出2種共成像點道集:一種是共成像點偏移距道集,另一種是共成像點角度道集。開展了針對地震波逆時偏移的2種共成像點道集波場特征及其差異的數(shù)值實驗。以傾斜界面模型和復(fù)雜的Marmousi模型為例,合成了2種共成像點道集,其中角度道集是在逆時延拓過程中采用波印廷矢量制作的。計算結(jié)果表明,在速度模型準(zhǔn)確的情況下,2種共成像點道集均可拉平,其中偏移噪聲在偏移距道集上的分布規(guī)律較差,不利于后續(xù)的精細處理,而在角度道集上的偏移噪聲主要集中在90°附近的高角度區(qū)域,對小角度范圍內(nèi)的地震道進行疊加時,逆時成像剖面的信噪比和地層的刻畫能力均可得到有效提高,同時低頻噪聲能量與界面的反射系數(shù)存在正相關(guān)性。

    地震波;疊前逆時偏移;共成像點偏移距道集;共成像點角度道集;噪聲分析

    0 引言

    地震波逆時成像技術(shù)是目前復(fù)雜構(gòu)造高精度成像中理論較為成熟的地震成像方法之一,其主要優(yōu)點在于能夠?qū)崿F(xiàn)多次波、回轉(zhuǎn)反射波等通常認(rèn)為是干擾波類型波場的成像,因此,更適合于解決任意陡傾角、速度變化較為劇烈情況時的地震波成像問題,并能夠解決克?;舴蚍e分法和單程波偏移方法的一些缺陷[1]?;谏鲜鲋T多優(yōu)點,目前把地震成像方法定位為逆時偏移,并得到了地球物理學(xué)界的充分重視和廣泛應(yīng)用[2-6]。這項技術(shù)隨著計算機技術(shù)的快速發(fā)展,較大程度地改善了逆時偏移技術(shù)工業(yè)化應(yīng)用的現(xiàn)狀,也為基于逆時偏移地震波全波形反演方法的工業(yè)化應(yīng)用奠定了基礎(chǔ)[7-8]。

    地震波逆時偏移是逐炮進行偏移的,無法實現(xiàn)對任意目標(biāo)線進行偏移,這一點類似于單程波算法。研究還表明,影響地震波逆時偏移計算效率的主要因素是由偏移網(wǎng)格、偏移頻率、偏移孔徑所確定的網(wǎng)格點個數(shù),三維情況通??蛇_到1 000萬個以上,這對目前計算機硬件資源及其計算能力提出了更高的要求。這正是20世紀(jì)80年代逆時偏移方法提出以來,近20多年沒有得到廣泛重視和發(fā)展的根本原因[9-14]。地震波逆時偏移的炮成像體按照共成像點進行抽取,并按照與炮點的距離進行排序,可以形成共成像點偏移距道集,計算結(jié)果為海量數(shù)據(jù),且所需磁盤空間和數(shù)據(jù)提取的計算量均較大,但制作過程相對比較簡單,因此在實際地震資料處理中得到了廣泛應(yīng)用。而地震波逆時偏移的角度道集制作過程較為復(fù)雜,計算量較大,隨著計算機硬件技術(shù)的發(fā)展,使得逆時偏移角度道集的制作和應(yīng)用逐漸成為現(xiàn)實。Save等[15-16]利用單程波偏移首先提取了地下局部偏移距域的共成像點偏移距道集,然后再根據(jù)速度轉(zhuǎn)換為角度道集,后來又提出了一種基于時延的共成像點偏移距道集轉(zhuǎn)換為角度道集的新方法,但所需的數(shù)據(jù)存儲量巨大;Fomel[17]將這種方法推廣到三維;Xu等[18-19]提出一種直接從頻率波數(shù)域提取角度道集的方法,這種方法提取的角度道集質(zhì)量較好,但需保存所有時刻的炮點和檢波點波場,計算量和存儲量巨大,因此不適合進行實際地震工區(qū)的逆時偏移等處理;Yoon等[20]采用坡印廷矢量來修改逆時成像條件在成像過程中實現(xiàn)逆時偏移低頻噪聲的有效壓制;王保利等[21]根據(jù)坡印廷矢量來確定波場傳播方向的問題,從而實現(xiàn)逆時偏移角度道集的提取。綜上可知,基于時延的共成像點角度道集提取方法的存儲量和計算量巨大,同時從共成像點偏移距道集轉(zhuǎn)為角度道集方式存在較大的誤差,因此,在實際地震資料逆時偏移角度道集提取中受到限制。而基于坡印廷矢量的角度道集提取方法具有計算過程簡單、計算量較小的優(yōu)點,在實際逆時偏移角度道集提取中具有潛在的應(yīng)用價值,同時,目前國內(nèi)外未見針對地震波逆時偏移方法的2種共成像點道集進行對比分析的報道,為此,筆者在王保利等[21]研究工作的基礎(chǔ)上,開展了2個理論模型的數(shù)值實驗研究,詳細探討地震波逆時偏移的共成像點偏移距道集和角度道集的波場特征及其差異,同時對這2種共成像點道集上的偏移噪聲特征及其成因進行對比分析,以期為地震資料的高精度逆時偏移處理提供指導(dǎo)。

    1 共成像點道集提取原理

    根據(jù)應(yīng)力和應(yīng)變關(guān)系,地震波波動方程通??梢员硎緸槿缦碌囊浑A雙曲型方程組形式[9]:

    式中:u為質(zhì)點的振動速度,m/s;m和n為便于應(yīng)用高階交錯網(wǎng)格有限差分法和PML吸收邊界條件而引入的2個中間變量;v為地球介質(zhì)速度,m/s;f(xs,zs,t)為震源函數(shù);xs和zs分別為點震源在x和z方向的空間坐標(biāo);t為時間,s。式(1)中的具體數(shù)值離散方法、吸收邊界條件及其穩(wěn)定性條件可參見文獻[1]。

    Yoon等[20]采用坡印廷矢量來確定波前的法線方向,即地震波場的傳播方向,計算公式為其中P為質(zhì)點的應(yīng)力場。以Ss和Sr分別表示炮點和檢波點波場的坡印廷矢量[21],根據(jù)余弦定理,可求得炮點檢波點這2個矢量之間的夾角(θ),計算公式為而實際地層界面的反射角為β=θ/2。這種角度道集的制作方法需要在每一步波場延拓過程中進行角度計算和角度道集的提取。當(dāng)角度道集內(nèi)的道數(shù)較多時,對應(yīng)的計算量和存儲量均較大。

    2 數(shù)值模擬

    2.1含傾斜層模型

    為了研究逆時成像技術(shù)在復(fù)雜構(gòu)造條件下的地震波逆時成像效果及其共成像點道集提取等問題,設(shè)計了含多個角度的傾斜界面速度模型[圖1(a)],用以模擬和指導(dǎo)實際地震資料逆時成像處理。該模型的橫向距離為2.5 km,共501個網(wǎng)格點,橫向網(wǎng)格間距為5 m;縱向深度為1 km,共201個網(wǎng)格點,縱向網(wǎng)格間距為5 m。模型頂層①和底層⑦厚度分別為150 m和250 m,對應(yīng)的介質(zhì)速度分別為1 800 m/s和3 500 m/s。模型中間部分包含一組不同傾角的傾斜界面,傾斜界面自右向左的傾角依次為30°,45°,60°和90°,劃分了5層介質(zhì),對應(yīng)的介質(zhì)序號及其速度分別是:層②為2 000 m/s,層③為2 300 m/s,層④為2 600 m/s,層⑤為3 000 m/s,層⑥為2 700 m/s。震源采用最大頻率為60 Hz的雷克子波,從模型最左側(cè)5 m處開始往右側(cè)激發(fā),震源深度均為10 m,炮間距為50 m,共激發(fā)50炮,檢波器布置于整個地表,道間距為5 m,共501個檢波器。數(shù)值模擬記錄的道長為2 s,時間步長為0.5 ms,滿足計算所需的穩(wěn)定性條件[9]。采用交錯網(wǎng)格有限差分法,其差分精度分別為時間2階和空間10階[1],在邊界處采用20個網(wǎng)格點的PML吸收邊界條件來壓制邊界反射波,可以保證截斷處具有較好的邊界吸收效果[22-23]。此外,逆時偏移角度道集的角度間隔為2°,因此,每個角度道集共含有45個地震道,且角度道集從左到右按照成像角度從小到大進行排序。而逆時偏移共成像點偏移距道集中含有的道數(shù)與炮數(shù)相關(guān),因此每個共成像點偏移距道集的總道數(shù)為50道,從左到右也按照炮點到成像點的絕對值距離由小到大進行排序。

    圖1 傾斜界面模型、逆時偏移剖面及其共成像點道集和角道集Fig.1 Tilted interface model and reverse-time migration section and its common imaging gathers and angle gathers

    分析圖1(b)可知,介質(zhì)⑥和介質(zhì)⑦的頂邊界均存在同相軸粗大和近似直流分量的波場現(xiàn)象(低頻逆時噪聲),同時還存在一些偏移劃弧噪聲,在淺層尤為嚴(yán)重。圖1(c)和圖1(d)分別為0°~45°和45°~90°的逆時偏移角道集疊加結(jié)果;圖1(e)和圖1(f)分別為500 m位置處的逆時偏移共成像點道集和角度道集;圖1(g)和圖1(h)分別為1 250 m位置處的逆時偏移共成像點道集和角度道集。分析逆時偏移共成像點偏移距道集[參見圖1(e)和圖1(g)]可知,偏移噪聲出現(xiàn)的位置是隨著分析點位置的變化而變化。例如,500 m位置處的偏移噪聲主要集中于共成像點偏移距道集的較大偏移距處,而1 250 m位置處的偏移噪聲主要集中于共成像點偏移距道集的較小偏移距處。同時,逆時偏移共成像點偏移距道集的同相軸還存在抖動現(xiàn)象,這主要是因為成像點位置兩側(cè)波場的逆時延拓路徑和成像波場類型不同造成的,在水平界面處差異較?。蹍⒁妶D1(e)],但在傾斜界面處差異較大[參見圖1(g)]。差異較大位置左側(cè)成像的波場主要來自于回轉(zhuǎn)反射波(經(jīng)介質(zhì)⑥與⑦的分界面向上層介質(zhì)回轉(zhuǎn)傳播)和多次波等波場,其成像能量相對較弱;而該位置右側(cè)成像的波場主要來自于一次反射波或多次波等,其成像能量相對較強。當(dāng)這2個方向的地震波場組合在一起進行排序就形成了共成像點偏移距道集同相軸的抖動現(xiàn)象。分析可知,在共成像點偏移距道集上無法有效地開展精細切除等后續(xù)處理,這主要是因為偏移噪聲在該道集上的分布范圍是不固定的。但在逆時偏移共成像點角度道集上[參見圖1(f)和圖1(h)],偏移噪聲(包括低頻噪聲)主要集中于高角度區(qū)域,因此,有利于在角度域進行切除等精細處理,以達到提高逆時成像精度的目的。為此,我們在逆時偏移角度道集上進行了分角度范圍疊加分析,在0°~45°角度道集的疊加結(jié)果上[參見圖1(c)],地層界面清晰,同相軸的一致性更好,并且同相軸粗大和偏移噪聲等問題得到了有效壓制。而偏移噪聲主要集中于45°~90°角度道集的疊加結(jié)果上,由此可見,在逆時偏移共成像點角度道集上開展精細處理工作有利于成像精度的提高。偏移噪聲與地層界面存在一定的相關(guān)性[參見圖1(d)],這主要是由炮點和檢波點波場在逆時延拓過程中沿著成像路徑對所有不相干的波場進行成像而引入的,在阻抗差異界面處既是不同類型波場(反射和透射等)形成的位置,也是不相干波場參與常規(guī)相關(guān)型逆時偏移成像條件計算的位置,因此,在這些阻抗差異界面附近形成了偏移噪聲,特別是低頻噪聲。此外,該低頻噪聲能量與分界面的反射系數(shù)存在正相關(guān)性,例如:介質(zhì)①與②的分界面、介質(zhì)⑥與⑦的分界面的反射系數(shù)分別為0.148 9和0.129 0,而介質(zhì)②與③的分界面、介質(zhì)⑤與⑦的分界面的反射系數(shù)分別為0.069 8和0.769 0,結(jié)合圖1(b)分析可知,反射系數(shù)越大,對應(yīng)的低頻能量越強,反之亦反。目前,這種低頻噪聲能量可通過應(yīng)用上下行波波場分離的逆時成像條件實現(xiàn)有效壓制[3]。

    2.2復(fù)雜的Marmousi模型

    為了進一步驗證共成像點角度道集在實際地震資料逆時偏移中的應(yīng)用效果,以國際標(biāo)準(zhǔn)的、更為復(fù)雜的Marmousi模型為例[圖2(a)、圖3(a)]來開展研究和分析。模型總大小為3.4 km×1.4 km,縱、橫向空間網(wǎng)格大小均為5 m,最小速度為1 028 m/s,最大速度為4 670 m/s,震源采用最大頻率為60 Hz的雷克子波,在地表距離模型最左側(cè)5 m處開始往右側(cè)激發(fā),震源深度為10 m,炮間距為50 m,共激發(fā)65炮。合成記錄的道長為2 s,時間步長為0.2 ms,滿足計算所需的穩(wěn)定性條件[9]。因此,每個角度道集共含有45道記錄,且角度道集從左到右按照角度從小到大進行排序。文中逆時偏移共成像點偏移距道集的道數(shù)為65道,也從左到右按照炮點到成像點的絕對值距離由小到大進行排序。這里采用2種地震觀測方式:一種是常規(guī)地面觀測系統(tǒng)采集(圖2),檢波器布置于整個地表,道間距為5 m,共680個檢波器;另一種是VSP觀測系統(tǒng)采集(圖3),檢波器布置于距離模型左側(cè)1.7 km處,道間距為5 m,共280個檢波器。

    圖2 Marmousi模型、逆時偏移剖面及其共成像點道集和角道集Fig.2 Marmousi model and reverse-time migration section and its common imaging gathers and angle gathers

    分析圖2(b)和圖3(b)可知,逆時偏移剖面中的偏移噪聲(包括低頻噪聲)均較為嚴(yán)重,模糊了地層的細節(jié)特征,因此信噪比較低。圖2(c)、圖3(c)和圖2(d)、圖3(d)分別為0°~45°和45°~90°的逆時偏移角道集疊加結(jié)果;圖2(e)和圖2(f)分別為1 750 m位置處的逆時偏移共成像點道集和角度道集;圖3(e)和圖3(f)分別為1 775 m位置處的逆時偏移共成像點道集和角度道集;圖2(g)和圖2(h)分別為2 750 m位置處的逆時偏移共成像點道集和角度道集;圖3(g)和圖3(h)分別為1 850 m位置處的逆時偏移共成像點道集和角度道集。分析逆時偏移共成像點偏移距道集[參見圖2(e)、圖2(g)、圖3(e)和圖3(g)]可知,偏移噪聲的分布在該道集上較為雜亂,無規(guī)律,若用統(tǒng)一的切除函數(shù)進行切除將較大程度地丟失有效信號的能量,從而影響最終的處理結(jié)果。與此同時,逆時偏移共成像點偏移距道集上也存在類似同相軸抖動現(xiàn)象,其成因機理是一致的,均是由不同空間方位成像的波場類型及其成像能量引起的。而在對應(yīng)的逆時偏移共成像點角度道集上[參見圖2(f)、圖2(h)、圖3(f)和圖3(h)],偏移噪聲基本集中于高角度位置,這有利于開展對偏移噪聲的精細切除等后續(xù)處理工作。筆者也對逆時偏移角度道集進行了分角度范圍疊加分析。在0°~45°角度道集的疊加結(jié)果上[參見圖2(c)和圖3(c)],斷面、小斷層、不整合界面以及地層層序界面均得到清晰地刻畫,偏移噪聲問題得到了有效壓制,且偏移噪聲主要集中于45°~90°的疊加結(jié)果上[參見圖2(d)和圖3(d)]。分析圖2(d)和圖3(d)可知,剖面中的低頻噪聲能量也與分界面的反射系數(shù)存在正相關(guān)性。由此可見,在逆時偏移共成像點角度道集上開展精細處理工作有利于提高逆時成像精度。

    圖3 Marmousi模型、VSP逆時偏移剖面及其共成像點道集和角道集Fig.3 Marmousi model and VSP reverse-time migration section and its common imaging gathers and angle gathers

    3 結(jié)束語

    筆者基于傾斜界面模型和復(fù)雜的Marmousi模型開展了針對地震波逆時偏移方法的2種共成像點道集(偏移距道集和角度道集)的數(shù)值實驗及其偏移噪聲分析研究。

    (1)在逆時偏移共成像點偏移距道集上的偏移噪聲分布不均勻、位置不固定,因此,較難在該道集上開展精細切除等處理工作。而在逆時偏移共成像點角度道集上,由于偏移噪聲主要集中于高角度區(qū)域,因此,能夠比較方便地進行精細切除和優(yōu)化處理,從而可提高逆時偏移的疊加質(zhì)量。

    (2)在去除的偏移噪聲中,低頻噪聲能量強弱與速度模型中的阻抗差異界面的反射系數(shù)存在正相關(guān)性,即反射系數(shù)越大,低頻能量越強,反之亦反。

    (3)無論是用地面觀測方式,還是用VSP觀測方式,地震波逆時偏移共成像點角度道集的優(yōu)越性是顯而易見的。同時逆時偏移角度道集還可為AVA分析和速度更新提供重要的數(shù)據(jù)基礎(chǔ),雖然目前其計算量仍較大,但隨著計算機技術(shù)的進步,必將在石油勘探開發(fā)中得到推廣應(yīng)用。

    [1]陳可洋.高階彈性波波動方程正演模擬及逆時偏移成像研究[D].大慶:大慶石油學(xué)院,2009.

    [2]Huang Y,Lin D,Bai B,et al.Challenges in presalt depth imaging of the deepwater Santos Basin,Brazil[J].The Leading Edge,2010,29(8):820-825.

    [3]陳可洋,吳清嶺,范興才,等.地震波疊前逆時偏移脈沖響應(yīng)研究與應(yīng)用[J].石油物探,2013,52(2):163-170.

    [4]Liu F Q,Zhang G,Morton S A,et al.Reverse-time migration using one way wavefield imaging condition[G].Expanded Abstracts of 77th Annual International SEG Meeting,2007:2170-2174.

    [5]Liu F Q,Zhang G,Morton S A,et al.An effective imaging condition for reverse-time migration using wave-field decomposition[J].Geophysics,2011,76(1):S29-S39.

    [6]陳可洋,范興才,吳清嶺,等.提高逆時偏移成像精度的疊前插值處理研究與應(yīng)用[J].石油物探,2013,52(4):409-416.

    [7]陳可洋.時限時移相關(guān)法疊前逆時成像條件及其應(yīng)用[J].石油物探,2011,50(1):22-26.

    [8]陳可洋.地震波逆時偏移方法研究綜述[J].勘探地球物理進展,2010,33(3):153-159.

    [9]陳可洋.基于高階有限差分的波動方程疊前逆時偏移方法[J].石油物探,2009,48(5):475-478.

    [10]Pavloske K,Beyoglu J.Pre-processing Improvements for Enhanced PSTM and TTI PSDM Images in the Santos Basin:A case history[G].Expanded Abstracts of 81st Annual International SEG Meeting,2011:3258-3262.

    [11]Yoon K,Marfurt K J.Reverse-time migration using the Poynting vector[J].Exploration Geophysics,2006,37(1):102-107.

    [12]Fletcher R P,F(xiàn)owler P J,Kitchenside P.Suppressing artifacts in prestack reverse time migration[G].Expanded Abstracts of 75th Annual International SEG Meeting,2005:2049-2051.

    [13]Guitton A,Kaelin B,Biondi B.Least-squares attenuation of reversetime migration artifacts[J].Geophysics,2007,72(1):S19-S23.

    [14]凌云.地震數(shù)據(jù)采集·處理·解釋·一體化實踐與探索[M].北京:石油工業(yè)出版社,2007.

    [15]Sava P C,F(xiàn)ormel S.Angle-domain common-image gathers by wavefield continuation methods[J].Geophysics,2003,68(3):1065-1074.

    [16]Sava P,F(xiàn)ormel S.Time-shift imaging condition in seismic migration[J].Geophysics,2006,71(6):S209-S217.

    [17]Fomel S.Theory of 3-D angle gathers in wave-equation imaging[G]. Expanded Abstracs of 74th Annual International SEG Meeting,2004:1053-1056.

    [18]Xu S,Zhang Y,Tang B.3D common image gathers from Reverse time migration[G].Expanded Abstracts of 80th Annual International SEG Meeting,2010:3257-3260.

    [19]Xu S,Zhang Y,Tang B.3D angle gathers from Reverse time migration[J].Geophysics,2011,76(2):S77-S92.

    [20]Yoon K,Marfurt Kurt J,Starr W.Challenges in reverse-time migration[G].Expanded Abstracts of 74th Annual International SEG Meeting,2003:1057-1060.

    [21]王保利,高靜懷,陳文超,等.逆時偏移中用Poynting矢量高效地提取角道集[J].地球物理學(xué)報,2013,56(1):262-268.

    [22]陳可洋.邊界吸收中鑲邊法的評價[J].中國科學(xué)院研究生院學(xué)報,2010,27(2):170-175.

    [23]陳可洋.完全匹配層吸收邊界條件研究[J].石油物探,2010,49(5):472-477.

    (本文編輯:楊琦)

    《巖性油氣藏》期刊被收錄信息

    《中國科學(xué)引文數(shù)據(jù)庫》(CSCD)來源期刊

    《中國核心期刊(遴選)數(shù)據(jù)庫》

    《中國期刊全文數(shù)據(jù)庫》

    《中國優(yōu)秀博碩士學(xué)位論文全文數(shù)據(jù)庫》

    《中國引文數(shù)據(jù)庫》

    《中文科技期刊數(shù)據(jù)庫》

    《中國石油文摘》

    Seismic wave reverse-time migration noise analysis within different common imaging point gathers

    CHEN Keyang1,WU Qingling1,F(xiàn)AN Xingcai1,CHEN Shumin1,LI Lailin1,LIU Zhenkuan2
    (1.Research Institute of Exploration and Development,PetroChina Daqing Oilfield Company Ltd.,Daqing 163712,Heilongjiang,China;2.Company of Exploration,PetroChina Daqing Oilfield Company Ltd.,Daqing 163453,Heilongjiang,China)

    There are two kinds of common imaging point gathers of the seismic imaging methods,one is the common imaging point offset gather,and the other is common imaging point angle gather.We carried out numerical experiments on the wave field characteristics and difference of two kinds of common imaging point gathers with relation of seismic wave reverse-time migration method.Taking tilted interface model and complex Marmousi model as examples,we composed two kinds of common imaging point gathers,of which the angle gather was made up by using Poynting vector during the reverse-time migration process.The computational results show that both the two kinds of common imaging point gathers are flat with correct velocity model.The noise distribution in the offset gather is not fixed,therefore,so it is not suitable for refining processing afterward,while the migration noise of the angle gather is almost focused on higher angle gather area,and with the branch angle stack processing,the low angle stacking section has the better S/N and the formation information is improved greatly.What’s more,the low frequency noise energy is positive proportional to the reflection coefficient of the interface.

    seismicwave;pre-stackreverse-timemigration;commonimagingpointoffsetgather;commonimagingpoint angle gather;noise analysis

    P631

    A

    1673-8926(2014)02-0118-07

    2013-11-14;

    2013-12-16

    國家重點基礎(chǔ)研究發(fā)展計劃(973)項目“火山巖油氣藏的形成機制與分布規(guī)律”(編號:2009CB219307)資助

    陳可洋(1983-),男,碩士,工程師,主要從事高精度地震波傳播模擬與逆時成像、多線程并行計算與模塊開發(fā)、地震資料數(shù)字處理方法與應(yīng)用等研究工作。地址:(163712)黑龍江省大慶市讓胡路區(qū)大慶石油勘探開發(fā)研究院地震處理二室。電話:(0459)5508524。E-mail:keyangchen@163.com。

    猜你喜歡
    波場角度界面
    神奇的角度
    國企黨委前置研究的“四個界面”
    一個涉及角度和的幾何不等式鏈的改進
    彈性波波場分離方法對比及其在逆時偏移成像中的應(yīng)用
    基于FANUC PICTURE的虛擬軸坐標(biāo)顯示界面開發(fā)方法研究
    角度不同
    37°女人(2017年8期)2017-08-12 11:20:48
    人啊
    滇池(2017年7期)2017-07-18 19:32:42
    人機交互界面發(fā)展趨勢研究
    交錯網(wǎng)格與旋轉(zhuǎn)交錯網(wǎng)格對VTI介質(zhì)波場分離的影響分析
    基于Hilbert變換的全波場分離逆時偏移成像
    啦啦啦中文免费视频观看日本| 国产视频内射| 97在线人人人人妻| 美女脱内裤让男人舔精品视频| 91aial.com中文字幕在线观看| 免费av中文字幕在线| 国产亚洲精品久久久com| 国产一级毛片在线| 欧美+日韩+精品| 91午夜精品亚洲一区二区三区| 色吧在线观看| 春色校园在线视频观看| 少妇精品久久久久久久| 久久久久久久久久人人人人人人| 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦中文免费视频观看日本| 韩国高清视频一区二区三区| 亚洲,欧美,日韩| 亚洲精品乱码久久久久久按摩| 777米奇影视久久| 国产成人免费无遮挡视频| 欧美3d第一页| 国产欧美亚洲国产| 国产精品秋霞免费鲁丝片| 日韩成人av中文字幕在线观看| 秋霞在线观看毛片| 午夜免费观看性视频| 免费观看在线日韩| 国产精品免费大片| 国产成人aa在线观看| 黄色欧美视频在线观看| 欧美精品一区二区大全| 自拍欧美九色日韩亚洲蝌蚪91 | 赤兔流量卡办理| h视频一区二区三区| 久久这里有精品视频免费| 国产成人午夜福利电影在线观看| 国产欧美亚洲国产| 欧美精品一区二区大全| 一级,二级,三级黄色视频| 国产免费福利视频在线观看| 亚洲无线观看免费| 国产精品人妻久久久久久| 大香蕉久久网| 肉色欧美久久久久久久蜜桃| 2022亚洲国产成人精品| 久久久久久久久久久免费av| 女性生殖器流出的白浆| 久久人妻熟女aⅴ| 两个人的视频大全免费| 欧美日韩在线观看h| 少妇丰满av| 色吧在线观看| 一本—道久久a久久精品蜜桃钙片| av又黄又爽大尺度在线免费看| 国产精品免费大片| 日本av手机在线免费观看| 久久国产亚洲av麻豆专区| 久久97久久精品| 精品少妇久久久久久888优播| 好男人视频免费观看在线| av不卡在线播放| 国产日韩欧美视频二区| 99久久人妻综合| 国产成人午夜福利电影在线观看| 国产熟女欧美一区二区| 亚洲熟女精品中文字幕| 国产欧美另类精品又又久久亚洲欧美| 在线观看美女被高潮喷水网站| 国产成人免费无遮挡视频| av女优亚洲男人天堂| 婷婷色综合www| 最近手机中文字幕大全| 国产淫片久久久久久久久| 亚洲精品中文字幕在线视频 | 亚洲精品国产av蜜桃| 免费看日本二区| av线在线观看网站| 久久精品熟女亚洲av麻豆精品| 国产成人一区二区在线| 亚洲欧美精品自产自拍| 日韩av在线免费看完整版不卡| 深夜a级毛片| 婷婷色综合大香蕉| 丰满人妻一区二区三区视频av| 一级a做视频免费观看| 一本色道久久久久久精品综合| 国产欧美日韩综合在线一区二区 | av天堂久久9| 精品一区在线观看国产| 一边亲一边摸免费视频| 亚洲精品久久午夜乱码| 91精品伊人久久大香线蕉| 大香蕉97超碰在线| a 毛片基地| 欧美成人午夜免费资源| 免费人成在线观看视频色| 国产在线视频一区二区| 亚洲精品乱久久久久久| 国产精品嫩草影院av在线观看| 日韩电影二区| 成人黄色视频免费在线看| 嫩草影院新地址| 免费黄网站久久成人精品| h日本视频在线播放| 亚洲国产av新网站| av女优亚洲男人天堂| 国产精品国产av在线观看| 免费高清在线观看视频在线观看| 在线观看免费高清a一片| 蜜桃久久精品国产亚洲av| 成人亚洲欧美一区二区av| 精品久久久久久电影网| 国产精品熟女久久久久浪| 大又大粗又爽又黄少妇毛片口| 天天操日日干夜夜撸| 在线观看国产h片| 免费看av在线观看网站| 国产精品久久久久久久电影| 精品一区在线观看国产| 国产无遮挡羞羞视频在线观看| 亚洲欧美日韩东京热| 免费播放大片免费观看视频在线观看| 亚洲高清免费不卡视频| 日韩欧美一区视频在线观看 | 亚洲精品国产av蜜桃| 国产av精品麻豆| 国产欧美另类精品又又久久亚洲欧美| 国产精品一区www在线观看| 另类亚洲欧美激情| 在线观看免费日韩欧美大片 | 亚洲欧美成人综合另类久久久| 国产毛片在线视频| 国产精品久久久久久久电影| av免费在线看不卡| 国产精品熟女久久久久浪| 久久国产精品男人的天堂亚洲 | 免费看av在线观看网站| kizo精华| 精品人妻偷拍中文字幕| 亚洲天堂av无毛| 国产精品三级大全| 最近2019中文字幕mv第一页| 一级毛片aaaaaa免费看小| 波野结衣二区三区在线| 国产伦精品一区二区三区视频9| 偷拍熟女少妇极品色| 少妇熟女欧美另类| 国产精品国产av在线观看| 纵有疾风起免费观看全集完整版| 十分钟在线观看高清视频www | 亚洲av免费高清在线观看| 在线免费观看不下载黄p国产| h视频一区二区三区| 夜夜看夜夜爽夜夜摸| 一区二区三区精品91| 欧美bdsm另类| 爱豆传媒免费全集在线观看| 丰满迷人的少妇在线观看| 在线观看免费高清a一片| 国产中年淑女户外野战色| 中文字幕亚洲精品专区| 在线观看一区二区三区激情| 久久av网站| 全区人妻精品视频| 久久99一区二区三区| 桃花免费在线播放| 国产精品国产三级国产专区5o| 看免费成人av毛片| 有码 亚洲区| 91精品国产九色| 免费观看在线日韩| 一级二级三级毛片免费看| 狠狠精品人妻久久久久久综合| 亚洲欧美一区二区三区黑人 | 国产爽快片一区二区三区| 亚洲国产日韩一区二区| 欧美bdsm另类| 97精品久久久久久久久久精品| 高清视频免费观看一区二区| 丰满人妻一区二区三区视频av| 亚洲精品第二区| 精品酒店卫生间| 国产男人的电影天堂91| 午夜精品国产一区二区电影| 国产永久视频网站| 性色avwww在线观看| 久久久欧美国产精品| 久久精品国产自在天天线| 五月天丁香电影| 国产毛片在线视频| 国产男女内射视频| 亚洲精品一区蜜桃| 成人黄色视频免费在线看| 国语对白做爰xxxⅹ性视频网站| 亚洲人成网站在线观看播放| 97超视频在线观看视频| 成年av动漫网址| 91久久精品电影网| 在线观看免费高清a一片| 成年女人在线观看亚洲视频| 久久久久久伊人网av| 久久人人爽人人片av| 少妇高潮的动态图| 我要看日韩黄色一级片| 青春草视频在线免费观看| 精品久久久久久久久亚洲| 国产精品嫩草影院av在线观看| 一个人看视频在线观看www免费| 女人精品久久久久毛片| 另类精品久久| 亚洲av日韩在线播放| 免费观看性生交大片5| 国产一区有黄有色的免费视频| 久久99蜜桃精品久久| 亚洲中文av在线| 日韩制服骚丝袜av| 人妻少妇偷人精品九色| 亚洲成人av在线免费| 久久久精品94久久精品| 精品人妻一区二区三区麻豆| 国产午夜精品久久久久久一区二区三区| 一本大道久久a久久精品| 欧美人与善性xxx| 自拍欧美九色日韩亚洲蝌蚪91 | 一本—道久久a久久精品蜜桃钙片| 日本黄色片子视频| 亚洲欧美精品自产自拍| videossex国产| 熟女av电影| 日韩大片免费观看网站| 成年av动漫网址| 国产黄色免费在线视频| 国产精品久久久久久精品电影小说| 大香蕉久久网| 色视频在线一区二区三区| 99热这里只有是精品50| 国产成人91sexporn| 日韩一区二区三区影片| 国产成人精品无人区| 久久精品久久精品一区二区三区| 少妇的逼水好多| 成人国产麻豆网| 人妻系列 视频| 熟女电影av网| 18禁在线无遮挡免费观看视频| 国产又色又爽无遮挡免| 日韩av免费高清视频| 99久久人妻综合| 亚洲精品,欧美精品| 色婷婷久久久亚洲欧美| 全区人妻精品视频| 久久97久久精品| av视频免费观看在线观看| 菩萨蛮人人尽说江南好唐韦庄| a级毛色黄片| 三级国产精品欧美在线观看| 丰满少妇做爰视频| 久久久国产精品麻豆| 丰满乱子伦码专区| 亚洲第一av免费看| 久久久久久久久大av| 91精品国产九色| 男女啪啪激烈高潮av片| 国产精品久久久久久精品电影小说| 免费av不卡在线播放| 天堂8中文在线网| 精品一品国产午夜福利视频| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久精品古装| 美女大奶头黄色视频| av国产久精品久网站免费入址| 亚洲国产精品999| 免费久久久久久久精品成人欧美视频 | a级毛片在线看网站| 能在线免费看毛片的网站| 欧美人与善性xxx| 国产精品99久久99久久久不卡 | 极品人妻少妇av视频| 久久99热6这里只有精品| 亚洲欧美一区二区三区黑人 | 成人毛片a级毛片在线播放| 热re99久久国产66热| 久久久久久久久久成人| av天堂久久9| 精品一区二区三卡| 777米奇影视久久| 乱人伦中国视频| 日本av免费视频播放| 欧美日韩av久久| 日本免费在线观看一区| 少妇猛男粗大的猛烈进出视频| 亚洲无线观看免费| 丰满迷人的少妇在线观看| 久久婷婷青草| 大码成人一级视频| 国产日韩欧美在线精品| 少妇人妻久久综合中文| 乱人伦中国视频| 婷婷色综合大香蕉| 女的被弄到高潮叫床怎么办| 69精品国产乱码久久久| 另类精品久久| 国产真实伦视频高清在线观看| 老司机亚洲免费影院| 人人妻人人澡人人爽人人夜夜| 天美传媒精品一区二区| 伦精品一区二区三区| av一本久久久久| 国产日韩欧美在线精品| 亚洲av福利一区| 丁香六月天网| 18禁裸乳无遮挡动漫免费视频| 国产精品蜜桃在线观看| 欧美bdsm另类| 内地一区二区视频在线| 国产精品秋霞免费鲁丝片| 国产精品蜜桃在线观看| 国产一区二区在线观看日韩| 亚洲精品国产av蜜桃| 欧美日韩视频精品一区| 一级av片app| 日日啪夜夜撸| 亚洲欧美日韩另类电影网站| 天天躁夜夜躁狠狠久久av| 免费观看的影片在线观看| 99视频精品全部免费 在线| 啦啦啦视频在线资源免费观看| 国产伦理片在线播放av一区| 人人澡人人妻人| 国产又色又爽无遮挡免| 99九九线精品视频在线观看视频| 亚洲人成网站在线观看播放| 国产高清国产精品国产三级| 中文资源天堂在线| 一级毛片久久久久久久久女| 亚洲无线观看免费| 黄色怎么调成土黄色| 欧美老熟妇乱子伦牲交| 在线观看美女被高潮喷水网站| 久久亚洲国产成人精品v| 免费久久久久久久精品成人欧美视频 | 久久精品国产亚洲av涩爱| 精品久久久久久久久av| 久久久久久久亚洲中文字幕| 免费观看a级毛片全部| 国产黄色视频一区二区在线观看| 欧美xxⅹ黑人| 精品国产国语对白av| 中文资源天堂在线| 久久久久人妻精品一区果冻| 少妇的逼水好多| 欧美精品国产亚洲| 免费黄频网站在线观看国产| 最新中文字幕久久久久| 嘟嘟电影网在线观看| 日韩人妻高清精品专区| 色婷婷久久久亚洲欧美| 我的老师免费观看完整版| 久久午夜福利片| 午夜福利,免费看| 一级毛片电影观看| 亚洲丝袜综合中文字幕| 国产欧美另类精品又又久久亚洲欧美| 日本91视频免费播放| 搡女人真爽免费视频火全软件| 妹子高潮喷水视频| 最新的欧美精品一区二区| 国产精品久久久久久精品电影小说| 黄色日韩在线| 好男人视频免费观看在线| 一个人看视频在线观看www免费| 久久ye,这里只有精品| 大又大粗又爽又黄少妇毛片口| 极品人妻少妇av视频| 精品99又大又爽又粗少妇毛片| 深夜a级毛片| 韩国高清视频一区二区三区| 一区二区三区精品91| 51国产日韩欧美| 三级国产精品欧美在线观看| 最近中文字幕高清免费大全6| 80岁老熟妇乱子伦牲交| 一级二级三级毛片免费看| 午夜视频国产福利| 夫妻午夜视频| 一级毛片aaaaaa免费看小| av又黄又爽大尺度在线免费看| 极品少妇高潮喷水抽搐| 亚洲av成人精品一二三区| 18禁动态无遮挡网站| 少妇人妻一区二区三区视频| 热re99久久精品国产66热6| a级毛片在线看网站| 五月天丁香电影| 夫妻性生交免费视频一级片| 99久久综合免费| 精品酒店卫生间| 久久97久久精品| 亚洲国产精品999| 久久午夜综合久久蜜桃| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 久久毛片免费看一区二区三区| 久久久久视频综合| 又黄又爽又刺激的免费视频.| 亚洲怡红院男人天堂| 精品99又大又爽又粗少妇毛片| 91午夜精品亚洲一区二区三区| 熟女人妻精品中文字幕| 日本黄大片高清| 少妇猛男粗大的猛烈进出视频| 中文字幕亚洲精品专区| 国产极品粉嫩免费观看在线 | 亚洲美女黄色视频免费看| 熟女人妻精品中文字幕| 亚洲精品乱码久久久久久按摩| 丝袜脚勾引网站| 边亲边吃奶的免费视频| 国产片特级美女逼逼视频| 欧美三级亚洲精品| 亚洲精品日本国产第一区| 久久女婷五月综合色啪小说| 国语对白做爰xxxⅹ性视频网站| 国产高清不卡午夜福利| 最近手机中文字幕大全| 三级国产精品欧美在线观看| 我要看黄色一级片免费的| 99精国产麻豆久久婷婷| 99久久综合免费| 亚洲经典国产精华液单| 一本久久精品| 最近中文字幕2019免费版| 六月丁香七月| 欧美日韩一区二区视频在线观看视频在线| 国产亚洲精品久久久com| 制服丝袜香蕉在线| 乱系列少妇在线播放| 在线看a的网站| 午夜免费鲁丝| 大又大粗又爽又黄少妇毛片口| 最新中文字幕久久久久| .国产精品久久| 国产精品成人在线| 色婷婷久久久亚洲欧美| 人人妻人人爽人人添夜夜欢视频 | 五月伊人婷婷丁香| 99九九在线精品视频 | 最新中文字幕久久久久| 成年人午夜在线观看视频| 国产精品无大码| 亚洲va在线va天堂va国产| 精品少妇黑人巨大在线播放| 视频中文字幕在线观看| 久久久国产精品麻豆| 国内精品宾馆在线| 亚洲精品乱码久久久久久按摩| 婷婷色av中文字幕| 99久久精品热视频| av专区在线播放| 天堂8中文在线网| 黄色日韩在线| 丝袜脚勾引网站| 91aial.com中文字幕在线观看| 亚洲美女搞黄在线观看| 久久久国产一区二区| 亚洲欧美成人精品一区二区| 在线观看免费高清a一片| 免费观看在线日韩| 三级国产精品片| 视频区图区小说| 丝袜喷水一区| 久久影院123| 黄色一级大片看看| av不卡在线播放| 亚洲精品久久久久久婷婷小说| 一级av片app| 国产成人免费观看mmmm| 久久久精品免费免费高清| 久久综合国产亚洲精品| 精品国产一区二区三区久久久樱花| 曰老女人黄片| 婷婷色麻豆天堂久久| 久久久久久久国产电影| 一级片'在线观看视频| 欧美老熟妇乱子伦牲交| 91精品伊人久久大香线蕉| 亚洲国产精品一区二区三区在线| 国产视频内射| 国产成人免费观看mmmm| 老女人水多毛片| 国产成人免费观看mmmm| av在线app专区| 久久精品夜色国产| 欧美xxⅹ黑人| 热99国产精品久久久久久7| 久久国产精品大桥未久av | videos熟女内射| 亚洲怡红院男人天堂| 蜜桃久久精品国产亚洲av| 亚洲欧洲国产日韩| 免费黄网站久久成人精品| 久久精品久久久久久噜噜老黄| 国产精品久久久久久久久免| 欧美日韩视频精品一区| 国产精品久久久久久久久免| 一区在线观看完整版| 欧美日韩亚洲高清精品| 欧美日韩视频精品一区| 亚洲国产欧美日韩在线播放 | 亚洲精品一二三| 中文精品一卡2卡3卡4更新| 久久久a久久爽久久v久久| 国产欧美另类精品又又久久亚洲欧美| 又黄又爽又刺激的免费视频.| av天堂久久9| 久久久久久久久久人人人人人人| 熟女人妻精品中文字幕| 少妇人妻 视频| 成人综合一区亚洲| 久久青草综合色| 亚洲av日韩在线播放| 男人狂女人下面高潮的视频| 五月玫瑰六月丁香| 亚洲欧美日韩东京热| 国产精品嫩草影院av在线观看| 深夜a级毛片| 一级爰片在线观看| 国产黄片视频在线免费观看| 国产精品一区二区性色av| 嘟嘟电影网在线观看| 五月玫瑰六月丁香| 国产欧美亚洲国产| 黄色视频在线播放观看不卡| 亚洲欧美成人综合另类久久久| 午夜免费观看性视频| 在线观看免费高清a一片| 在线观看av片永久免费下载| 九九爱精品视频在线观看| 国产精品免费大片| 国产黄片美女视频| 黄色视频在线播放观看不卡| 亚洲美女视频黄频| 成人午夜精彩视频在线观看| 精品卡一卡二卡四卡免费| 日韩一本色道免费dvd| 中文字幕免费在线视频6| 有码 亚洲区| h日本视频在线播放| 99热6这里只有精品| 日日撸夜夜添| 亚洲,欧美,日韩| 五月天丁香电影| av福利片在线| 欧美变态另类bdsm刘玥| 成人无遮挡网站| 国产极品天堂在线| 亚洲图色成人| 国产精品一区二区在线不卡| 夫妻性生交免费视频一级片| 三级国产精品片| 爱豆传媒免费全集在线观看| 国产一区二区三区综合在线观看 | 天堂俺去俺来也www色官网| 少妇人妻久久综合中文| 人妻一区二区av| 日韩亚洲欧美综合| 少妇丰满av| 国产又色又爽无遮挡免| 国产在视频线精品| 中文欧美无线码| 亚洲三级黄色毛片| 美女福利国产在线| 男人舔奶头视频| av在线app专区| 人人妻人人添人人爽欧美一区卜| 亚洲欧美精品自产自拍| 精品99又大又爽又粗少妇毛片| 久久久精品94久久精品| 熟女电影av网| 欧美日韩亚洲高清精品| 国国产精品蜜臀av免费| 久久久久久久久久久免费av| 青春草国产在线视频| 亚洲av福利一区| 久久久久久久久久久免费av| 国产精品国产三级国产专区5o| 国产欧美另类精品又又久久亚洲欧美| 成年人免费黄色播放视频 | 国产成人精品久久久久久| 国产欧美日韩一区二区三区在线 | 在线观看www视频免费| 午夜激情久久久久久久| 黄色日韩在线| 免费看av在线观看网站| 国产精品三级大全| 99视频精品全部免费 在线| 91精品国产九色| 97在线视频观看| 熟女电影av网| 亚洲国产成人一精品久久久| 伦理电影免费视频| 欧美xxⅹ黑人| 国产深夜福利视频在线观看| 偷拍熟女少妇极品色| 成人毛片a级毛片在线播放| av黄色大香蕉| 少妇裸体淫交视频免费看高清| 国产精品国产av在线观看| 久久久a久久爽久久v久久| 中国国产av一级| 永久网站在线| 国产深夜福利视频在线观看| 人人澡人人妻人| 天天躁夜夜躁狠狠久久av|