陳可洋，范興才，吳清嶺，陳樹(shù)民，李來(lái)林，劉振寬

（1.中國(guó)石油天然氣集團(tuán)公司大慶油田有限責(zé)任公司勘探開(kāi)發(fā)研究院，黑龍江大慶163712；2.中國(guó)石油天然氣集團(tuán)公司大慶油田有限責(zé)任公司勘探事業(yè)部，黑龍江大慶163453）

地震波逆時(shí)偏移成像技術(shù)是目前復(fù)雜構(gòu)造高精度成像中最精確、理論較為成熟的成像方法之一［1］。由于逆時(shí)偏移方法對(duì)地震波動(dòng)方程的近似較少，適用于任意陡傾角和速度在縱、橫向變化較為劇烈情況的疊前偏移成像問(wèn)題，可以實(shí)現(xiàn)多次波、回轉(zhuǎn)反射波等通常認(rèn)為是干擾波類型波場(chǎng)的準(zhǔn)確成像問(wèn)題，這正是當(dāng)前地震波疊前逆時(shí)成像技術(shù)受到地球物理學(xué)界的充分重視并得到廣泛應(yīng)用的原因［2－6］。目前，隨著計(jì)算機(jī)技術(shù)的快速發(fā)展，特別是基于CPU／GPU高性能集群并行快速計(jì)算技術(shù)和大容量磁盤(pán)的快速存儲(chǔ)技術(shù)的出現(xiàn)，較大程度地改善了逆時(shí)偏移技術(shù)的工業(yè)化應(yīng)用現(xiàn)狀，同時(shí)也為基于逆時(shí)偏移處理的地震波全波形反演方法的工業(yè)化應(yīng)用奠定了基礎(chǔ)［7－13］。

如何提高地震波疊前逆時(shí)成像精度的研究是當(dāng)前地球物理學(xué)界關(guān)注的核心內(nèi)容之一，是疊前逆時(shí)成像技術(shù)廣泛應(yīng)用于實(shí)際地震資料處理的理論基礎(chǔ)。目前，這方面的工作可以歸納為如下幾個(gè)方面［14－31］：①地震資料的預(yù)處理；②逆時(shí)成像條件的改進(jìn)；③成像結(jié)果的后續(xù)處理。地震資料預(yù)處理主要包括噪聲的保幅、高保真壓制和速度高精度建模。Pavloske等［14］對(duì)此進(jìn)行了研究。逆時(shí)偏移要求輸入的預(yù)處理炮集為未經(jīng)球面擴(kuò)散補(bǔ)償?shù)臄?shù)據(jù)，這明顯不同于克?；舴蚱品椒?，是炮域波動(dòng)方程偏移固有的要求。目前工業(yè)界的逆時(shí)成像條件主要采用炮檢點(diǎn)波場(chǎng)相關(guān)處理，但這會(huì)引入較強(qiáng)的低頻背景噪聲，降低了計(jì)算結(jié)果的信噪比。目前針對(duì)逆時(shí)成像條件也已有較多的研究工作，且取得了較好的應(yīng)用效果。如Yoon等［16］提出了在零延遲互相關(guān)條件中加入波印廷矢量來(lái)消除成像噪聲；Liu等［17－18］把全波場(chǎng)分解成單程波分量，并運(yùn)用逆時(shí)成像條件結(jié)合這些單程波波場(chǎng)分量達(dá)到消除成像噪聲的目的。逆時(shí)成像結(jié)果的后續(xù)處理是提高逆時(shí)成像結(jié)果精度的有效辦法［19］，如 Guitton等［20］采用最小二乘預(yù)測(cè)誤差濾波器消除成像噪聲；Zhang等［21］指出拉普拉斯算子濾波相當(dāng)于成像波場(chǎng)角度域衰減；陳可洋［22］研究了低頻噪聲的來(lái)源，并在疊加域、成像點(diǎn)域、共炮點(diǎn)域進(jìn)行了拉普拉斯算子去噪分析。

地震道插值是提高地震資料處理和偏移成像精度的一種有效方法，在地震資料處理流程中得到廣泛使用。由于野外采集常受到施工條件的綜合影響（如地形、障礙物等），致使采集數(shù)據(jù)中常常在某些位置出現(xiàn)地震記錄的道間距過(guò)大、異常道、數(shù)據(jù)不完整或數(shù)據(jù)缺失等問(wèn)題，這意味著一些地層的地震響應(yīng)信息的丟失，并可能在各種處理過(guò)程中產(chǎn)生許多不必要的噪聲。從實(shí)際應(yīng)用狀況來(lái)看，地震道插值方法可以在一定程度上恢復(fù)這些缺失的地震波場(chǎng)信息，從而提高偏移成像的精度，并切實(shí)有效地壓制噪聲和抑制空間假頻。目前，地震道插值已成為疊前偏移資料預(yù)處理中的一個(gè)關(guān)鍵步驟，但其在逆時(shí)偏移資料預(yù)處理中的應(yīng)用尚未見(jiàn)相關(guān)的研究報(bào)道。

我們?cè)谇叭搜芯康幕A(chǔ)上，開(kāi)展了針對(duì)逆時(shí)偏移的地震資料預(yù)處理方法研究——疊前地震道插值處理?；诶碚撃Ｐ蛿?shù)據(jù)和實(shí)際三維地震資料，在滿足逆時(shí)偏移方法數(shù)值頻散關(guān)系的條件下，通過(guò)加密疊前地震數(shù)據(jù)的道間距和線間距，來(lái)論證地震道插值處理對(duì)逆時(shí)偏移成像結(jié)果的影響。

1 理論模型分析

以Marmousi速度模型（圖1）為例，模型總大小為3 400m×1 400m；縱向和橫向空間網(wǎng)格步長(zhǎng)均為5m；模型的最小速度為1 028m／s，最大速度為4 670m／s；密度均為1g／cm3。采用主頻為40Hz的Ricker子波作為震源在地表激發(fā)，時(shí)間步長(zhǎng)為0.2ms，滿足計(jì)算所需的穩(wěn)定性條件，數(shù)值模擬的有限差分精度為時(shí)間2階、空間12階，可以保證地震波正演模擬結(jié)果具有較高的數(shù)值計(jì)算精度（無(wú)明顯的數(shù)值頻散現(xiàn)象）［23－24］。在人工截?cái)噙吔缣幉捎猛耆ヅ鋵樱≒ML）吸收邊界條件［25－26］，吸收邊界網(wǎng)格節(jié)點(diǎn)數(shù)為20個(gè)，并采用正弦和余弦型阻尼衰減因子，能夠保證計(jì)算結(jié)果最大程度地壓制邊界反射波，提高有效模擬區(qū)域計(jì)算結(jié)果的信噪比。每炮接收時(shí)間均為2s，在整個(gè)地表均布置檢波點(diǎn)（圖1），排列長(zhǎng)度固定為3 400m，道間距為5m，第1炮置于理論模型地表最左側(cè)0處，炮間距為20m，按照該觀測(cè)系統(tǒng)采集得到171個(gè)共炮點(diǎn)道集。采用常規(guī)相關(guān)法疊前逆時(shí)成像條件對(duì)產(chǎn)生的171個(gè)共炮點(diǎn)道集進(jìn)行抽稀疊前逆時(shí)深度偏移處理，采用的逆時(shí)偏移頻率、網(wǎng)格步長(zhǎng)等參數(shù)與正演模擬參數(shù)相同，且滿足數(shù)值頻散關(guān)系［27－28］，能夠保證逆時(shí)偏移過(guò)程中不存在明顯的數(shù)值頻散問(wèn)題。

圖2為根據(jù)觀測(cè)系統(tǒng)采集的炮集數(shù)值模擬記錄在道間距抽稀前（道間距為5m）和抽稀后（道間距為20m）的逆時(shí)偏移結(jié)果。其中，圖2a和圖2b分別為第1炮數(shù)據(jù)道間距抽稀前和抽稀后的逆時(shí)偏移結(jié)果；圖2c和圖2d分別為所有炮集數(shù)據(jù)道間距抽稀前和抽稀后的逆時(shí)偏移疊加結(jié)果；圖2e和圖2f分別為圖2c和圖2d虛線框內(nèi)的局部放大圖。對(duì)比圖2a和圖2b可知，單炮數(shù)據(jù)道間距抽稀后的逆時(shí)偏移結(jié)果存在一些空間假頻，淺層較強(qiáng)（圖2b中紅色圓圈），深層相對(duì)較弱（圖2b中藍(lán)色圓圈），但線性干擾特征也較為明顯，而疊前數(shù)據(jù)道間距抽稀前的逆時(shí)偏移結(jié)果（圖2a）不存在該現(xiàn)象，精度較高。對(duì)比圖2c和圖2d可知，從宏觀上看，兩張剖面的構(gòu)造特征基本相似，僅在偏移結(jié)果的能量分布上存在差異，但通過(guò)對(duì)比局部放大圖（圖2e和圖2f）可見(jiàn)，疊前數(shù)據(jù)道間距抽稀后的逆時(shí)偏移疊加結(jié)果存在明顯的線性干擾（圖2f中藍(lán)色圓圈），類似于空間假頻，影響了同相軸的連續(xù)性；而抽稀前的逆時(shí)偏移疊加剖面同相軸的連續(xù)性較好，無(wú)空間假頻存在。同時(shí)注意到，該理論模型的逆時(shí)偏移參數(shù)是完全依據(jù)數(shù)值頻散關(guān)系，不存在由數(shù)值頻散問(wèn)題引入的干擾波存在。由此表明，疊前道集數(shù)據(jù)的道間距稀疏程度對(duì)逆時(shí)偏移結(jié)果的精度（特別是淺層成像）存在較大的影響，有必要對(duì)疊前道集進(jìn)行道間距插值來(lái)達(dá)到提高逆時(shí)成像精度的目的。

圖1 Marmousi速度模型及其觀測(cè)方式

圖2 模擬記錄道間距抽稀前、后的逆時(shí)偏移結(jié)果

2 實(shí)際資料處理

研究對(duì)象為松遼盆地某區(qū)塊的實(shí)際三維地震資料，該地區(qū)近地表?xiàng)l件較為簡(jiǎn)單，工區(qū)原始地震資料的信噪比較高，中淺層成像是該區(qū)塊的逆時(shí)成像難點(diǎn)。地震資料的原始采集面元為20m×40m，檢波點(diǎn)方向道間距為40m，檢波站方向線間距為160m。采用精細(xì)的地震速度優(yōu)化建模技術(shù)和高保真地震資料預(yù)處理技術(shù)，為該工區(qū)地震資料的逆時(shí)偏移處理提供了高保真的炮集數(shù)據(jù)和最終的深度域速度模型。應(yīng)用GPU高性能并行計(jì)算技術(shù)實(shí)現(xiàn)加速逆時(shí)偏移處理（偏移網(wǎng)格為10m×10m，深度步長(zhǎng)為5m），同時(shí)保證該工區(qū)采用的逆時(shí)偏移參數(shù)均滿足數(shù)值頻散關(guān)系，由此開(kāi)展疊前道間插值對(duì)逆時(shí)偏移結(jié)果的影響研究。

我們采用三維傅里葉頻率－空間域疊前插值方法，共包括4個(gè)具體計(jì)算步驟［29］：①對(duì)原始地震數(shù)據(jù)沿時(shí)間方向作傅里葉變換，得到頻率－空間域的地震數(shù)據(jù)；②估計(jì)原始地震數(shù)據(jù)中線性同相軸的數(shù)量，用最小平方算法求取預(yù)測(cè)算子；③依據(jù)內(nèi)插算子與預(yù)測(cè)算子的對(duì)偶關(guān)系，計(jì)算內(nèi)插算子；④對(duì)內(nèi)插算子應(yīng)用最小平方原理計(jì)算內(nèi)插出的地震數(shù)據(jù)。

由于沿檢波站方向的道間距為40m，道數(shù)為80道，而沿檢波線方向的道間距為160m，且線數(shù)僅為12條，因此，優(yōu)選道數(shù)較多的方向（檢波站方向）先進(jìn)行插值，再對(duì)道數(shù)稀疏方向（檢波線方向）進(jìn)行插值，從而保證每次插值結(jié)果具有較高的可靠性。具體插值處理分為3步：第1步為沿檢波站方向的道間插值（內(nèi)插1道），使道間距由原來(lái)的40m變?yōu)?0m，該方向的道數(shù)則由原80道變?yōu)?59道；第2步是在第1步道間插值的基礎(chǔ)上再沿檢波線方向的道間內(nèi)插1道（一次線內(nèi)插），該方向的道間距由原160m變?yōu)?0m，道數(shù)則由原來(lái)的12道變?yōu)?3道；第3步是在第2步一次線內(nèi)插的基礎(chǔ)上仍沿檢波線方向的道間再內(nèi)插1道（二次線內(nèi)插），該方向的道間距再由80m變?yōu)?0m，道數(shù)則由23道變?yōu)?5道。

圖3a為第1步道間插值前、后的地震道位置變化圖；圖4a和圖4b分別為第1步道間插值前、后的單炮數(shù)據(jù)。對(duì)比分析圖4a和圖4b可知，疊前三維地震道插值后單炮記錄的同相軸連續(xù)性更好，插值道的波形過(guò)渡自然，無(wú)畸變現(xiàn)象。圖3b為第2步一次線內(nèi)插前、后的地震道位置變化圖；圖4c和圖4d分別為第2步一次線內(nèi)插前、后的單炮數(shù)據(jù)。圖3c為第3步二次線內(nèi)插前、后的地震道位置變化圖；圖4e為第3步二次線內(nèi)插后的單炮數(shù)據(jù)。對(duì)比圖4c，圖4d和圖4e可見(jiàn)，疊前單炮數(shù)據(jù)由稀疏變?yōu)槊芗?，單炮記錄上的同相軸連續(xù)性得到更為明顯的提高，同時(shí)插值道的波形過(guò)渡均較為自然，無(wú)任何畸變現(xiàn)象。

圖5為圖4對(duì)應(yīng)的單炮記錄經(jīng)逆時(shí)偏移處理后的結(jié)果。對(duì)比圖5a和圖5b可知，在第1步道間插值前的單炮記錄逆時(shí)偏移結(jié)果上，淺層存在較強(qiáng)能量的線性偏移噪聲（紅色橢圓內(nèi)），類似于第1章理論分析部分疊前數(shù)據(jù)抽稀后的逆時(shí)偏移結(jié)果中的干擾波場(chǎng)特征；而在第1步道間插值后的單炮記錄逆時(shí)偏移結(jié)果上，該線性噪聲得到較好的壓制，同相軸的連續(xù)性更好，同時(shí)逆時(shí)偏移結(jié)果（藍(lán)色矩形部分）的同相軸連續(xù)性也得到較為明顯的提高。對(duì)比圖5c，圖5d和圖5e可知，在第2步一次線內(nèi)插前的單炮記錄逆時(shí)偏移結(jié)果（圖5c）中，淺層存在較強(qiáng)能量的線性偏移噪聲（藍(lán)色矩形和紅色橢圓內(nèi)），而經(jīng)第2步和第3步兩次線內(nèi)插后，線間距變?yōu)?0（圖5d）和40m（圖5e），該線性偏移噪聲得到了有效壓制，且同相軸的連續(xù)性和偏移結(jié)果的信噪比均得到有效提高，同時(shí)地層的細(xì)節(jié)特征也更加清晰。

圖5 實(shí)際單炮記錄插值前、后沿檢波站和檢波線方向逆時(shí)偏移結(jié)果

圖6a和圖6b分別給出了疊前三維插值處理前沿檢波站方向和沿檢波線方向的逆時(shí)偏移疊加剖面；圖6c和圖6d分別給出了疊前三維插值處理后（檢波站方向道間距插值為20m，檢波線方向線間距插值為40m）沿檢波站方向和沿檢波線方向的逆時(shí)偏移疊加剖面，最終偏移網(wǎng)格為10m×10m。分析圖6可知，疊前插值處理前，沿檢波站方向和沿檢波線方向的逆時(shí)偏移疊加結(jié)果均存在一定能量的線性干擾（藍(lán)色矩形和紅色橢圓內(nèi)），影響了同相軸的連續(xù)性，不利于準(zhǔn)確刻畫(huà)中淺層的薄互層特征；而經(jīng)插值處理后，該線性干擾現(xiàn)象得到了明顯的壓制，同時(shí)淺層的信噪比得到提高，同相軸的連續(xù)性更好，這與第1章中對(duì)疊前數(shù)據(jù)抽稀前的逆時(shí)偏移結(jié)果分析得到的結(jié)論相一致。由此可見(jiàn)，疊前插值處理在一定程度上能夠提高逆時(shí)偏移結(jié)果的精度和橫向連續(xù)性。

值得指出的是，常規(guī)克?；舴蚍e分法深度偏移可采用反假頻算子較好地克服空間假頻的問(wèn)題，單程波深度偏移也可在頻率－波數(shù)域引入反假頻算子來(lái)實(shí)現(xiàn)這一目的，但目前逆時(shí)偏移方法無(wú)法引入反假頻算子；因此，要提高逆時(shí)偏移的頻率（可提高橫、縱向的分辨率），滿足精細(xì)地質(zhì)勘探的要求，必須減小偏移網(wǎng)格（如面元和深度步長(zhǎng)），以滿足逆時(shí)偏移算法本身所要求的數(shù)值頻散關(guān)系（與速度、偏移網(wǎng)格、差分精度等相關(guān)）。我們建議在信噪比較高的實(shí)際地震資料逆時(shí)偏移處理中，可以適當(dāng)減小偏移網(wǎng)格，如減小為原始采集面元的1／2或1／4，并對(duì)原始的疊前地震數(shù)據(jù)進(jìn)行1道或3道插值（如圖3所示），能夠獲得較好的逆時(shí)偏移處理效果。此外，在逆時(shí)偏移參數(shù)（如偏移網(wǎng)格、孔徑、深度及其步長(zhǎng)、頻率等）給定的情況下，疊前地震數(shù)據(jù)插值結(jié)果對(duì)逆時(shí)偏移計(jì)算效率的影響較小，僅在疊前數(shù)據(jù)讀取時(shí)消耗較少的I／O時(shí)間。

圖6 實(shí)際三維地震資料插值前、后疊前逆時(shí)偏移疊加剖面

3 結(jié)束語(yǔ)

我們從理論模型和高信噪比的實(shí)際地震資料出發(fā)，詳細(xì)分析研究了疊前插值對(duì)提高逆時(shí)偏移成像精度的影響，取得了以下認(rèn)識(shí)與成果。

1）對(duì)比基于Marmousi理論模型的合成炮集數(shù)據(jù)抽稀前、后的逆時(shí)偏移結(jié)果可知，在保證逆時(shí)偏移參數(shù)滿足數(shù)值頻散關(guān)系條件下，數(shù)據(jù)抽稀直接降低了逆時(shí)偏移結(jié)果的精度，引入了線性干擾波場(chǎng)，影響了橫向分辨率，并形成了錯(cuò)斷特征的虛假小斷層，這對(duì)中淺層如薄互層等地質(zhì)特征的逆時(shí)成像是不利的。

2）經(jīng)疊前頻率－空間域插值處理后，炮域數(shù)據(jù)插值道的同相軸過(guò)渡自然無(wú)畸變，對(duì)應(yīng)逆時(shí)偏移結(jié)果的線性干擾得到有效壓制，同相軸的連續(xù)性和橫向分辨率也得到有效提高。由此表明，對(duì)于資料信噪比較高的地震工區(qū)，可以采用疊前插值方法來(lái)提高逆時(shí)成像精度，從而彌補(bǔ)檢波線線間距和檢波站道間距稀疏采樣對(duì)逆時(shí)成像結(jié)果的影響。

3）必須指出，在逆時(shí)偏移參數(shù)給定的情況下，疊前插值道數(shù)的多少對(duì)逆時(shí)偏移效率的影響較?。坏逯档罃?shù)過(guò)多會(huì)引入虛假的成像信息，這是因?yàn)椴逯档啦⒎钦鎸?shí)的地震道信息，僅是地震響應(yīng)波場(chǎng)的某種近似。此外，對(duì)于資料信噪比較低的三維地震工區(qū)，不建議采用疊前插值處理，此時(shí)容易使得插值道波形發(fā)生畸變，并引入較多的干擾波場(chǎng)，從而降低逆時(shí)偏移成像質(zhì)量。

［1］Yoon K，Marfurt K，Starr E W.Challenges in reverse－time migration［J］.Expanded Abstracts of 74thAnnual Internat SEG Mtg，2004，1057－1060

［2］陳可洋.高階彈性波波動(dòng)方程正演模擬及逆時(shí)偏移成像研究［D］.大慶：大慶石油學(xué)院，2009 Chen K Y.High－order elastic wave equation forward modeling and reverse－time migration［D］.Daqing：Daqing Petroleum Institute，2009

［3］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

［4］劉紅偉，劉洪，鄒振.地震疊前逆時(shí)偏移中的去噪與存儲(chǔ)［J］.地球物理學(xué)報(bào)，2010，53（9）：2171－2180 Liu H W，Liu H，Zou Z，et al.The problems of denoise and storage in seismic reverse time migration［J］.Chinese Journal of Geophysics，2010，53（9）：2171－2180

［5］陳可洋.兩種類型共炮點(diǎn)域相關(guān)型疊前深度成像方法［J］.內(nèi)陸地震，2012，26（1）：17－27 Chen K Y.Two kinds of correlation type pre－stack depth migration methods in common shot domain［J］.Inland Earthquake，2012，26（1）：17－27

［6］陳可洋.地震波逆時(shí)偏移方法研究綜述［J］.勘探地球物理進(jìn)展，2010，33（3）：153－159 Chen K Y.Reviews on seismic wave reverse－time migration methods［J］.Progress in Exploration Geophysics，2010，33（3）：153－159

［7］楊仁虎，常旭，劉伊克.疊前逆時(shí)偏移影響因素分析［J］.地球物理學(xué)報(bào)，2010，53（8）：1902－1913 Yang R H，Chang X，Liu Y K.The influence factors analysis of imaging precision in pre－stack reverse time migration［J］.Chinese Journal of Geophysics，2010，53（8）：1902－1913

［8］陳可洋.基于高階有限差分的波動(dòng)方程疊前逆時(shí)偏移方法［J］.石油物探，2009，48（5）：475－478 Chen K Y.Wave equation pre－stack reverse－time migration scheme based on high－order finite－difference［J］.Geophysical Prospecting for Petroleum，2009，48（5）：475－478

［9］陳可洋.正演子波響應(yīng)特征及逆時(shí)成像分析［J］.石油物探，2011，50（5）：455－462 Chen K Y.Seismic wavelet response characteristic and reverse－time migration analysis［J］.Geophysical Prospecting for Petroleum，2011，50（5）：455－462

［10］陳可洋.地震波旅行時(shí)計(jì)算方法及其模型試驗(yàn)分析［J］.石油物探，2010，49（2）：153－157 Chen K Y.Seismic wave travel－time computational method and its model experiment analysis［J］.Geophysical Prospecting for Petroleum，2010，49（2）：153－157

［11］陳可洋.地震波數(shù)值模擬中差分階數(shù)的邊界效應(yīng)［J］.高原地震，2011，23（1）：20－23 Chen K Y.Difference order’s boundary effect of the seismic wave numerical simulation［J］.Plateau Earthquake Research，2011，23（1）：20－23

［12］陳可洋.地震波數(shù)值模擬中優(yōu)化的通量校正傳輸方法［J］.內(nèi)陸地震，2012，26（3）：169－179 Chen K Y.An optimized flux－corrected transport scheme in the seismic wave numerical simulation［J］.Inland Earthquake，2012，26（2）：169－179

［13］陳可洋，吳清嶺，范興才，等.地震波疊前逆時(shí)偏移脈沖響應(yīng)研究與應(yīng)用［J］.石油物探，2013，52（2）：163－170 Chen K Y，Wu Q L，F(xiàn)an X C，et al.Study on seismic wave pre－stack reverse－time migration’s impulse response and its application［J］.Geophysical Prospecting for Petroleum，2013，52（2）：163－170

［14］Pavloske K，Beyoglu J.Pre－processing improvements for enhanced PSTM and TTI PSDM images in the Santos Basin：a case history［J］.Expanded Abstracts of 81stAnnual Internat SEG Mtg，2011，3258－3262

［15］陳可洋.時(shí)限時(shí)移相關(guān)法疊前逆時(shí)成像條件及其應(yīng)用［J］.石油物探，2011，50（1）：22－26 Chen K Y.Time－limited time－shift correlation prestack reverse－time image condition and its application［J］.Geophysical Prospecting for Petroleum，2010，50（1）：22－26

［16］Yoon K，Marfurt K J.Reverse－time migration using the Poynting vector［J］.Exploration Geophysics，2006，37（1）：102－107

［17］Liu F Q，Zhang G，Morton S A，et al.Reverse－time migration using one way wavefield imaging condition［J］.Expanded Abstracts of 77thAnnual Internat SEG Mtg，2007，2170－2174

［18］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

［19］Fletcher R P，F(xiàn)owler P J，Kitchenside P.Suppressing artifacts in prestack reverse time migration［J］.Expanded Abstracts of 75thAnnual Internat SEG Mtg，2005，2049－2051

［20］Guitton A，Kaelin B，Biondi B.Least－squares attenuation of reverse－time migration artifacts［J］.Geophysics，2007，72（1）：S19－S23

［21］Zhang Y，James S.Practical issues of reverse time migration：true amplitude gathers，noise removal and harmonic source encoding［C］∥CPS／SEG Beijing 2009International Geophysical Cofference Technical Program Expanded Abstracts.Beijing：CPS＆SEG，2009：5

［22］陳可洋.基于拉普拉斯算子的疊前逆時(shí)噪聲壓制方法［J］.巖性油氣藏，2011，23（5）：87－95 Chen K Y.Pre－stack reverse－time noise suppressing method based on Laplacian operator［J］.Lithologic Reservoirs，2011，23（5）：87－95

［23］陳可洋.標(biāo)量聲波波動(dòng)方程高階交錯(cuò)網(wǎng)格有限差分法［J］.中國(guó)海上油氣，2009，21（4）：232－236 Chen K Y.High－order staggered－grid finite difference scheme for scalar acoustic wave equation［J］.China Offshore Oil and Gas，2009，21（4）：232－236

［24］陳可洋.地震波數(shù)值模擬中差分近似的各向異性分析［J］.石油物探，2010，49（1）：19－22 Chen K Y.Anisotropic analysis of difference approximation in seismic wave numerical modeling［J］.Geophysical Prospecting for Petroleum，2010，49（1）：19－22

［25］陳可洋.完全匹配層吸收邊界條件研究［J］.石油物探，2010，49（5）：472－477 Chen K Y.Study on perfectly matched layer absorbing boundary condition［J］.Geophysical Prospecting for Petroleum，2010，49（5）：472－477

［26］陳可洋.邊界吸收中鑲邊法的評(píng)價(jià)［J］.中國(guó)科學(xué)院研究生院學(xué)報(bào)，2010，27（2）：170－175 Chen K Y.Evaluation on the bordering method of the absorbing boundary condition［J］.Journal of the Graduate School of the Chinese Academy of Sciences，2010，27（2）：170－175

［27］陳可洋.網(wǎng)格剖分及其精度和計(jì)算量分析［J］.內(nèi)陸地震，2011，25（1）：12－20 Chen K Y.Mesh generation and its accuracy and computational amount analysis［J］.Inland Earthquake，2010，25（1）：12－20

［28］陳可洋.數(shù)值模擬的尺度調(diào)節(jié)機(jī)理［J］.內(nèi)陸地震，2011，25（4）：321－328 Chen K Y.Scale adjusting mechanism on numerical simulation［J］.Inland Earthquake，2011，25（4）：321－328

［29］張文穎.保幅性三維傅里葉變換疊前道內(nèi)插技術(shù)研究［D］.北京：中國(guó)石油大學(xué)，2011 Zhang W Y.Amplitude－preserved prestack interpolation of seismic data with 3－D Fourier Transform［D］.Beijing：China University of Petroleum，2011