交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析

韓 冬 杜啟振, 張明強(qiáng) 郭成鋒 公緒飛

1) 中國(guó)山東青島266580中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院2) 中國(guó)天津300451中海油田服務(wù)股份有限公司

?

交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析

1) 中國(guó)山東青島266580中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院2) 中國(guó)天津300451中海油田服務(wù)股份有限公司

利用交錯(cuò)網(wǎng)格有限差分和旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分進(jìn)行各向異性介質(zhì)彈性波場(chǎng)數(shù)值模擬時(shí)， 質(zhì)點(diǎn)振動(dòng)速度分量與應(yīng)力張量的網(wǎng)格節(jié)點(diǎn)定義方式均不相同， 從而對(duì)各向異性波場(chǎng)分離效果產(chǎn)生不同的影響． 針對(duì)這一問(wèn)題， 本文以具有垂直對(duì)稱軸的橫向各向同性(VTI)介質(zhì)的波場(chǎng)分離為例， 首先分析了兩種網(wǎng)格的參數(shù)定義方式以及VTI介質(zhì)波場(chǎng)的分離過(guò)程； 其次， 詳細(xì)研究和分析了這兩種網(wǎng)格的參數(shù)定義方式對(duì)各向異性介質(zhì)波場(chǎng)分離的影響， 并依據(jù)波前面連續(xù)性以及波場(chǎng)分離效果等方面， 通過(guò)數(shù)值模擬實(shí)驗(yàn)對(duì)該影響進(jìn)行分析驗(yàn)證. 結(jié)果表明， 旋轉(zhuǎn)交錯(cuò)網(wǎng)格的參數(shù)定義方式更有利于進(jìn)行各向異性介質(zhì)波場(chǎng)數(shù)值模擬和波場(chǎng)分離．

波場(chǎng)分離 交錯(cuò)網(wǎng)格 旋轉(zhuǎn)交錯(cuò)網(wǎng)格 有限差分 VTI介質(zhì)

引言

隨著勘探技術(shù)的發(fā)展， 多分量地震技術(shù)逐漸受到業(yè)界的重視(孫歧峰， 杜啟振， 2011)． 為了充分利用多分量資料的縱橫波信息， 處理流程通常都包括波場(chǎng)分離環(huán)節(jié)(董敏煜， 2002)． 依據(jù)不同的處理方法， 可將波場(chǎng)分離技術(shù)分為兩大類(李慶忠， 王建花， 2000). 第一類是基于標(biāo)量波場(chǎng)理論的波場(chǎng)分離技術(shù)， 即對(duì)地震記錄實(shí)現(xiàn)縱波與橫波的波場(chǎng)分離． 工業(yè)界通常將多分量地震記錄的垂直分量視為P波， 并利用坐標(biāo)旋轉(zhuǎn)得到徑向(R)和切向(T)分量， 然后采用相對(duì)成熟的聲波處理流程對(duì)上述分量分別進(jìn)行處理． 除此之外， 還可利用極化濾波(Etgen， 1988； 李志遠(yuǎn)等， 2013)、f-k濾波(Amundsenetal， 1998)、τ-P變換(許世勇等， 1999)、 拉東(Randon)變換(馮晅等， 2011)、 平面波分解(Devaney， Oristaglio， 1986)等方法實(shí)現(xiàn)對(duì)地震記錄縱橫波的波場(chǎng)分離． 第二類是基于矢量波場(chǎng)理論的分離技術(shù)， 即在波場(chǎng)傳播過(guò)程中實(shí)現(xiàn)對(duì)波場(chǎng)快照的波場(chǎng)分離． 在各向同性介質(zhì)方面， Sun和McMechan(2001)基于亥姆霍茲(Helmholtz)分解理論提出利用散度和旋度算子實(shí)現(xiàn)P波與S波的波場(chǎng)分離， 但會(huì)引入90°相移． 針對(duì)這一問(wèn)題， Sun等(2001)利用希爾伯特(Hilbert)變換實(shí)現(xiàn)了相位校正． 相對(duì)于各向同性介質(zhì)， 各向異性介質(zhì)波場(chǎng)分離的難點(diǎn)在于質(zhì)點(diǎn)振動(dòng)方向與波場(chǎng)傳播方向并不完全平行或正交， 因此不能直接利用散度和旋度算子進(jìn)行波場(chǎng)分離(Yan， Sava， 2009a， b)． Dellinger和Etgen(1990)依據(jù)克里斯托弗爾(Christoffel)方程， 在波數(shù)域中將彈性波場(chǎng)向偏振方向投影， 從而得到qP波和qS波． 在此基礎(chǔ)上， Yan和Sava(2009a， b)進(jìn)一步計(jì)算得到空間域不同尺寸空間窗的波場(chǎng)分離算子， 并利用該算子與波場(chǎng)值卷積得到qP波和qS波． 此外， Zhang和McMechan(2010)還提出了一種在波數(shù)域獲得每個(gè)笛卡爾坐標(biāo)分量上保幅qP波和qS波的波場(chǎng)分解方法． 與第一類波場(chǎng)分離技術(shù)相比， 第二類波場(chǎng)分離技術(shù)更能夠描述波場(chǎng)在地下介質(zhì)中傳播的矢量特性． 因此， 針對(duì)第二類波場(chǎng)分離技術(shù)展開(kāi)研究具有重要意義．

彈性波場(chǎng)數(shù)值模擬是實(shí)現(xiàn)第二類波場(chǎng)分離的前提． 有限差分法由于其計(jì)算效率高、 占用內(nèi)存小以及便于實(shí)現(xiàn)等優(yōu)點(diǎn)而成為彈性波場(chǎng)數(shù)值模擬的重要工具． 相對(duì)于中心規(guī)則網(wǎng)格有限差分法， 交錯(cuò)網(wǎng)格有限差分法(Madariaga， 1976)的精度更高且頻散壓制效果更好， 因此被廣泛應(yīng)用于各種介質(zhì)的波場(chǎng)數(shù)值模擬中. 例如： Virieux(1984， 1986)利用交錯(cuò)網(wǎng)格有限差分法模擬了各向同性介質(zhì)中SH波、 P波及SV波的傳播； Carcione和Helle(1999)利用交錯(cuò)網(wǎng)格有限差分法實(shí)現(xiàn)了孔隙介質(zhì)中彈性波場(chǎng)的數(shù)值模擬； 裴正林(2004)利用交錯(cuò)網(wǎng)格高階有限差分法實(shí)現(xiàn)了三維各向異性介質(zhì)彈性波數(shù)值模擬． 但交錯(cuò)網(wǎng)格有限差分法的不足之處在于質(zhì)點(diǎn)振動(dòng)的速度分量定義在不同的網(wǎng)格點(diǎn)上， 利用這種差分格式進(jìn)行數(shù)值模擬得到的波場(chǎng)值不能夠直接應(yīng)用到Y(jié)an和Sava(2009a， b)所提出的各向異性介質(zhì)波場(chǎng)分離方法中． 而旋轉(zhuǎn)交錯(cuò)網(wǎng)格(Saengeretal， 2000)將所有的應(yīng)力分量定義在同一網(wǎng)格點(diǎn)上， 將所有的質(zhì)點(diǎn)振動(dòng)速度分量定義在另一網(wǎng)格點(diǎn)上， 且由旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分得到的波場(chǎng)值能夠直接應(yīng)用到各向異性介質(zhì)波場(chǎng)分離過(guò)程中． 同時(shí)， 采用旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分進(jìn)行數(shù)值模擬時(shí)， 彈性模量無(wú)需進(jìn)行插值處理(嚴(yán)紅勇， 劉洋， 2012)． 綜上可見(jiàn)， 旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分更適于各向異性介質(zhì)波場(chǎng)數(shù)值模擬及其波場(chǎng)分離過(guò)程．

本文擬以具有垂直對(duì)稱軸的橫向各向同性(transversely isotropic with a vertical axis， 簡(jiǎn)寫為VTI)介質(zhì)的波場(chǎng)分離為例， 從波場(chǎng)分離的實(shí)現(xiàn)過(guò)程出發(fā)， 對(duì)比分析交錯(cuò)網(wǎng)格和旋轉(zhuǎn)交錯(cuò)網(wǎng)格的參數(shù)定義方式對(duì)VTI介質(zhì)波場(chǎng)分離(Yan， Sava， 2009a)造成的影響， 并通過(guò)數(shù)值試驗(yàn)對(duì)其進(jìn)行分析驗(yàn)證， 以進(jìn)一步說(shuō)明旋轉(zhuǎn)交錯(cuò)網(wǎng)格在各向異性波場(chǎng)數(shù)值模擬及其波場(chǎng)分離方面的優(yōu)勢(shì).

1 VTI介質(zhì)一階速度-應(yīng)力方程及兩種網(wǎng)格定義方式

二維情況下， VTI介質(zhì)一階速度-應(yīng)力方程(Virieux， 1986)為

(1)

式中：vx和vz分別為x方向和z方向的質(zhì)點(diǎn)速度分量；τxx和τzz為正應(yīng)力，τxz為切應(yīng)力；C11,C13,C33和C55為剛度系數(shù)；ρ為密度． 分別采用交錯(cuò)網(wǎng)格有限差分(Madariaga， 1976)和旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分(Saengeretal， 2000)對(duì)式(1)進(jìn)行數(shù)值模擬．

圖1給出了交錯(cuò)網(wǎng)格(Virieux， 1986)和旋轉(zhuǎn)交錯(cuò)網(wǎng)格(Saengeretal， 2000)應(yīng)力場(chǎng)和速度場(chǎng)的定義方式示意圖， 其中： 實(shí)線表示定義在時(shí)間上的主網(wǎng)格， 虛線表示在時(shí)間t上交錯(cuò)了Δt/2的半網(wǎng)格；O點(diǎn)為參考點(diǎn)，P1點(diǎn)和P2點(diǎn)分別代表某一點(diǎn)x方向和z方向上的質(zhì)點(diǎn)振動(dòng)速度分量． 比較圖1a與圖1b可知， 質(zhì)點(diǎn)振動(dòng)的速度分量在交錯(cuò)網(wǎng)格中沒(méi)有定義在同一個(gè)網(wǎng)格節(jié)點(diǎn)上， 而在旋轉(zhuǎn)交錯(cuò)網(wǎng)格中則定義在同一個(gè)網(wǎng)格節(jié)點(diǎn)上．

圖1 交錯(cuò)網(wǎng)格(a)和旋轉(zhuǎn)交錯(cuò)網(wǎng)格(b)示意圖

2 VTI介質(zhì)波場(chǎng)分離過(guò)程

VTI介質(zhì)中， qP波與qS波的偏振方向相互垂直， 利用偏振矢量將波場(chǎng)投影至偏振方向可得到qP波和qS波(Dellinger， Etgen， 1990)． 在波數(shù)域中， 通過(guò)求解克里斯托弗爾方程， 可得qP波和qS波的偏振矢量. 二維情況下求解

(2)

可得qP波的偏振矢量UqP為

(3)

qSV波的偏振矢量UqSV為

(4)

圖2 空間域波場(chǎng)分離算子. (a) x方向; (b) z方向

3 兩種交錯(cuò)網(wǎng)格對(duì)波場(chǎng)分離的影響

這里以qP波的波場(chǎng)分離過(guò)程為例， 詳細(xì)闡述兩種交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離所產(chǎn)生的影響．

波數(shù)域和空間域中的qP波(Yan，Sava， 2009a)分別表示為

(5)

(6)

波數(shù)域和空間域中各速度分量的傅里葉變換對(duì)為

(7)

(8)

如圖1a所示， 在交錯(cuò)網(wǎng)格中， O點(diǎn)處的質(zhì)點(diǎn)振動(dòng)速度分量可利用P1和P2點(diǎn)處的質(zhì)點(diǎn)振動(dòng)速度分量表示為

(9)

如圖1b所示， 在旋轉(zhuǎn)交錯(cuò)網(wǎng)格中， O點(diǎn)處的質(zhì)點(diǎn)振動(dòng)速度分量可利用P1和P2點(diǎn)處的質(zhì)點(diǎn)振動(dòng)速度分量表示為

(10)

將式(9)和式(10)帶入式(6)， 可得

(11)

(12)

此即為空間域中的qP波.

Yan和Sava(2009a)提出的VTI介質(zhì)波場(chǎng)分離方法是將Vx(kx, kz)和Vz(kx, kz)以及vx(x, z)和vz(x, z)視為同一個(gè)網(wǎng)格節(jié)點(diǎn)上的質(zhì)點(diǎn)振動(dòng)速度分量進(jìn)而實(shí)現(xiàn)波場(chǎng)分離． 對(duì)比式(11)與式(12)可得， 由于交錯(cuò)網(wǎng)格質(zhì)點(diǎn)振動(dòng)速度分量未定義在同一個(gè)網(wǎng)格節(jié)點(diǎn)上， 導(dǎo)致用于VTI介質(zhì)波場(chǎng)分離的x方向質(zhì)點(diǎn)振動(dòng)速度分量比z方向多出平移算子exp(-i2πkxΔx/2), z方向質(zhì)點(diǎn)振動(dòng)速度分量比x方向多出平移算子exp(-i2πkzΔz/2)， 從而造成波場(chǎng)分離不徹底. 針對(duì)這一問(wèn)題，Zhang和McMechan(2010)利用兩點(diǎn)平滑估算其中間值， 從而得到同一個(gè)網(wǎng)格節(jié)點(diǎn)上的質(zhì)點(diǎn)振動(dòng)速度分量， 但該方法只是對(duì)波場(chǎng)中間值的近似， 并不能完全表示兩點(diǎn)中間的波場(chǎng)值．Du等(2014)提出在波數(shù)域中對(duì)波場(chǎng)值進(jìn)行平移校正， 從而實(shí)現(xiàn)兩點(diǎn)中間值的準(zhǔn)確求取； 但鑒于本文采用的VTI介質(zhì)波場(chǎng)分離方法(Yan，Sava， 2009a)的處理流程是在空間域中進(jìn)行的， 這種波數(shù)域中的插值方法將會(huì)造成額外的計(jì)算量， 故不適用于本文．

如式(12)所述， 由旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分得到的用于波場(chǎng)分離的x方向與z方向的質(zhì)點(diǎn)振動(dòng)速度分量在網(wǎng)格節(jié)點(diǎn)上的對(duì)應(yīng)位置保持一致， 不需要進(jìn)行額外的插值處理， 同時(shí)可避免因波場(chǎng)插值造成的波場(chǎng)值不準(zhǔn)確及效率降低等問(wèn)題． 另外， 在使用旋轉(zhuǎn)交錯(cuò)網(wǎng)格進(jìn)行波場(chǎng)數(shù)值模擬的過(guò)程中也不需要對(duì)彈性模量進(jìn)行插值處理.

4 數(shù)值算例

本節(jié)將采用均勻VTI介質(zhì)模型進(jìn)行彈性波場(chǎng)數(shù)值模擬和波場(chǎng)分離， 通過(guò)模擬測(cè)試結(jié)果對(duì)比分析交錯(cuò)網(wǎng)格和旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離(Yan，Sava， 2009a)造成的影響. 本文采用的模型網(wǎng)格大小為600×300， 水平方向和垂直方向的空間采樣間隔均為10m， 時(shí)間采樣間隔為0.001s， 模型參數(shù)如表1所示． 交錯(cuò)網(wǎng)格有限差分和旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分均采用空間8階、 時(shí)間2階進(jìn)行彈性波場(chǎng)數(shù)值模擬， 利用非分裂完全匹配層(non-splittingperfectmatchedlayer， 簡(jiǎn)寫為NPML)對(duì)邊界反射進(jìn)行吸收衰減． 震源子波選取主頻為35Hz的雷克子波， 其解析式為f(t)=[1-2(πfpt)2]exp[-(πfpt)2]， 式中fp為子波主頻； 震源位置坐標(biāo)為(2800m, 500m)．

表1 均勻VTI介質(zhì)參數(shù)表

在彈性波場(chǎng)數(shù)值模擬過(guò)程中， 提取第400個(gè)時(shí)間步長(zhǎng)(即第400ms)的質(zhì)點(diǎn)振動(dòng)速度x分量與z分量的波場(chǎng)快照進(jìn)行波場(chǎng)分離， 以闡明這兩種交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離所造成的影響．

圖3和圖6分別為交錯(cuò)網(wǎng)格和旋轉(zhuǎn)交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬得到的波場(chǎng)快照. 對(duì)比圖3與圖6可知， 采用兩種交錯(cuò)網(wǎng)格有限差分進(jìn)行波場(chǎng)數(shù)值模擬均能得到清晰的波場(chǎng)快照． 圖4給出了直接利用交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬得到的波場(chǎng)快照進(jìn)行波場(chǎng)分離得到的qP波和qS波； 圖5給出了利用交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬并進(jìn)行兩點(diǎn)平均插值處理后的波場(chǎng)快照進(jìn)行波場(chǎng)分離得到的qP波和qS波(Zhang，McMechan， 2010)． 可以看出， 經(jīng)過(guò)兩點(diǎn)插值處理后進(jìn)行波場(chǎng)分離得到的結(jié)果要優(yōu)于未采用波場(chǎng)插值處理直接進(jìn)行波場(chǎng)分離的結(jié)果； 但如圖4和圖5中箭頭所指位置所示，qP波中仍然殘留有qS波，qS波中同樣殘留有qP波， 且qS波中的殘留現(xiàn)象更為嚴(yán)重． 圖7為旋轉(zhuǎn)交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬得到的波場(chǎng)快照進(jìn)行波場(chǎng)分離的結(jié)果． 比較圖4、 圖5及圖7中箭頭和方框所示可知， 使用旋轉(zhuǎn)交錯(cuò)網(wǎng)格數(shù)值模擬得到的波場(chǎng)值進(jìn)行波場(chǎng)分離得到的qP波和qS波較為干凈徹底， 且qP波的連續(xù)性較好．

圖3 交錯(cuò)網(wǎng)格數(shù)值模擬得到的波場(chǎng)快照(400 ms). (a) x分量； (b) z分量

圖4 交錯(cuò)網(wǎng)格數(shù)值模擬得到的波場(chǎng)快照直接進(jìn)行波場(chǎng)分離的結(jié)果. (a) qP波; (b) qS波

圖5 交錯(cuò)網(wǎng)格數(shù)值模擬得到波場(chǎng)快照并進(jìn)行兩點(diǎn)平均插值處理后進(jìn)行波場(chǎng)分離得到的結(jié)果. (a) qP波; (b) qS波

圖6 旋轉(zhuǎn)交錯(cuò)網(wǎng)格數(shù)值模擬得到的波場(chǎng)快照(400 ms). (a) x分量； (b) z分量

圖7 旋轉(zhuǎn)交錯(cuò)網(wǎng)格數(shù)值模擬得到的波場(chǎng)快照直接進(jìn)行波場(chǎng)分離得到的結(jié)果. (a) qP波; (b) qS波

圖8 VTI 介質(zhì)Hess模型. (a) δ; (b) ε; (c) ρ; (d) vP0; (e) vS0

圖9 Hess模型交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬結(jié)果及波場(chǎng)分離結(jié)果(a) x分量波場(chǎng)快照； (b) z分量波場(chǎng)快照； (c) qP波； (d) qS波

圖10Hess模型旋轉(zhuǎn)交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬結(jié)果及波場(chǎng)分離結(jié)果(a) x分量波場(chǎng)快照； (b) z分量波場(chǎng)快照； (c)qP波； (d)qS波

Fig.10SnapshotsforHessmodelbyusingrotatedstaggered-gridnumericalsimulationandtheseparationresults. (a) x-componentsnapshot； (b) z-componentsnapshot； (c)qP-wave； (d)qS-wave

為了進(jìn)一步說(shuō)明使用旋轉(zhuǎn)交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬得到的波場(chǎng)值進(jìn)行波場(chǎng)分離的優(yōu)勢(shì)， 本文采用如圖8所示的VTI介質(zhì)Hess模型進(jìn)行波場(chǎng)數(shù)值模擬和波場(chǎng)分離． 模型網(wǎng)格大小為600×350， 水平方向和垂直方向的采樣間隔為10m， 時(shí)間采樣間隔為0.001s； 采用主頻為35Hz的雷克子波， 震源位置坐標(biāo)為(3200m, 700m)； 交錯(cuò)網(wǎng)格有限差分和旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分均采用空間8階、 時(shí)間2階進(jìn)行彈性波場(chǎng)數(shù)值模擬； 采用非分裂完全匹配層對(duì)邊界反射進(jìn)行吸收衰減， 同時(shí)選取600ms的波場(chǎng)快照進(jìn)行波場(chǎng)分離． 圖9和圖10分別給出了利用交錯(cuò)網(wǎng)格和旋轉(zhuǎn)交錯(cuò)網(wǎng)格進(jìn)行數(shù)值模擬得到的波場(chǎng)快照以及對(duì)波場(chǎng)快照進(jìn)行波場(chǎng)分離的結(jié)果， 對(duì)比其中紅色箭頭標(biāo)示處可驗(yàn)證， 使用旋轉(zhuǎn)交錯(cuò)網(wǎng)格波場(chǎng)數(shù)值模擬得到的波場(chǎng)值進(jìn)行波場(chǎng)分離得到的qP波與qS波分離效果良好．

5 討論與結(jié)論

本文針對(duì)交錯(cuò)網(wǎng)格和旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離造成的影響進(jìn)行了詳細(xì)分析， 并選取簡(jiǎn)單模型和復(fù)雜模型進(jìn)行了相應(yīng)的數(shù)值實(shí)驗(yàn)， 得到如下結(jié)論：

1) 直接采用交錯(cuò)網(wǎng)格有限差分所得波場(chǎng)快照進(jìn)行波場(chǎng)分離， 無(wú)法徹底分離qP波與qS波； 盡管通過(guò)對(duì)模擬結(jié)果進(jìn)行簡(jiǎn)單的算術(shù)平均插值處理后再進(jìn)行波場(chǎng)分離， 能夠在一定程度上提高波場(chǎng)分離的效果， 但仍不能徹底解決波場(chǎng)無(wú)法完全分離的問(wèn)題．

2) 采用旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分所得波場(chǎng)快照進(jìn)行波場(chǎng)分離， 能夠得到分離較為徹底的qP波和qS波， 且無(wú)需進(jìn)行任何插值處理．

因此， 當(dāng)對(duì)各向異性介質(zhì)進(jìn)行波場(chǎng)數(shù)值模擬并進(jìn)行相應(yīng)的波場(chǎng)分離時(shí)， 旋轉(zhuǎn)交錯(cuò)網(wǎng)格有限差分是一種較為理想的選擇．

評(píng)審專家和編輯部對(duì)本文提出了寶貴意見(jiàn)和建議， 中國(guó)石油大學(xué)(華東)EWEP課題組提供了幫助， 作者在此表示衷心的感謝.

董敏煜. 2002. 多波多分量地震勘探[M]. 北京: 石油工業(yè)出版社: 76--83.

DongMY. 2002. Multi-Wave and Multi-Component Seismic Exploration[M].Beijing:PetroleumIndustryPress: 76--83 (inChinese).

馮晅， 張先武， 劉財(cái)， 王典， 楊慶節(jié). 2011. 帶有多道相關(guān)的拋物線Radon變換法分離P-P、P-SV波[J]. 地球物理學(xué)報(bào)， 54(2): 304--309.

FengX，ZhangXW，LiuC，WangD，YangQJ. 2011.SeparatingP-PandP-SVwavebyparabolicRadontransformwithmultiplecoherence[J]. Chinese Journal of Geophysics， 54(2): 304--309 (inChinese).

李慶忠， 王建花. 2000. 多波地震勘探的難點(diǎn)與展望[M]. 青島: 中國(guó)海洋大學(xué)出版社: 35--36.

LiQZ，WangJH. 2000. Problems of Multi-Component Seismic Exploration and Its Solutions[M].Qingdao:OceanUniversityofChinaPress: 35--36 (inChinese).

李志遠(yuǎn)， 梁光河， 谷丙洛. 2013. 基于散度和旋度縱橫波分離方法的改進(jìn)[J]. 地球物理學(xué)報(bào)， 56(6): 2012--2022.

LiZY，LiangGH，GuBL. 2013.ImprovedmethodofseparatingP-andS-wavesusingdivergenceandcurl[J]. Chinese Journal of Geophysics， 56(6): 2012--2022 (inChinese).

裴正林. 2004. 三維各向異性介質(zhì)中彈性波方程交錯(cuò)網(wǎng)格高階有限差分法數(shù)值模擬[J]. 石油大學(xué)學(xué)報(bào)： 自然科學(xué)版， 28(5): 23--29.

PeiZL. 2004.Three-dimensionalnumericalsimulationofelasticwavepropagationin3-Danisotropicmediawithstaggered-gridhigh-orderdifferencemethod[J]. Journal of China University of Petroleum： Edition of Natural Sciences， 28(5): 23--29 (inChinese).

孫歧峰， 杜啟振. 2011. 多分量地震數(shù)據(jù)處理技術(shù)研究現(xiàn)狀[J]. 石油勘探與開(kāi)發(fā)， 38(1): 67--73.

SunQF，DuQZ. 2011.Areviewofthemulti-componentseismicdataprocessing[J]. Petroleum Exploration and Development， 38(1): 67--73 (inChinese).

許世勇， 李彥鵬， 馬在田. 1999. τ-q變換法波場(chǎng)分離[J]. 中國(guó)海上油氣: 地質(zhì)， 13(5): 334--337.

XuSY，LiYP，MaZT. 1999.SeparationofP-andS-wavefieldsviatheτ-qtransform[J]. China Offshore Oil and Gas: Geology， 13(5): 334--337 (inChinese).

嚴(yán)紅勇， 劉洋. 2012. 黏彈TTI介質(zhì)中旋轉(zhuǎn)交錯(cuò)網(wǎng)格高階有限差分?jǐn)?shù)值模擬[J]. 地球物理學(xué)報(bào)， 55(4): 1354--1365.

YanHY，LiuY. 2012.Rotatedstaggeredgridhigh-orderfinite-differencenumericalmodelingforwavepropagationinviscoelasticTTImedia[J]. Chinese Journal of Geophysics， 55(4): 1354--1365 (inChinese).

AmundsenL，IkelleLT，MartinJ. 1998.MultipleattenuationandP/SsplittingofOBCdata:Aheterogeneousseafloor[C]∥Proceedings of the 68th Annual International Meeting， SEG， Expanded Abstracts.NewOrleans:SocietyofExplorationGeophysicists: 722--725.

CarcioneJM，HelleHB. 1999.Numericalsolutionoftheporoviscoelasticwaveequationonastaggeredmesh[J]. J Comput Phys， 154(2): 520--527.

DellingerJ，EtgenJ. 1990.Wave-fieldseparationintwo-dimensionalanisotropicmedia[J]. Geophysics， 55(7): 914--919.

DevaneyAJ，OristaglioML. 1986.Aplane-wavedecompositionforelasticwavefieldsappliedtotheseparationofP-wavesandS-wavesinvectorseismicdata[J]. Geophysics， 51(2): 419--423.

DuQZ，ZhangMQ，ChenXR，GongXF，GuoCF. 2014.True-amplitudewavefieldseparationusingstaggered-gridinterpolationinthewavenumberdomain[J]. Appl Geophys， 11(4): 437--446.

EtgenJT. 1988.PrestackmigrationofP-andSV-waves[C]∥Proceedings of the 58th Annual International Meeting， SEG， Expanded Abstract.Tulsa:SocietyofExplorationGeophysicists: 972--975.

MadariagaR. 1976.Dynamicsofanexpandingcircularfault[J]. Bull Seismol Soc Am， 66(3): 639--666.

SaengerEH，GoldN，ShapiroSA. 2000.Modelingthepropagationofelasticwavesusingamodifiedfinite-differencegrid[J]. Wave Motion， 31(1): 77--92.

SunR，ChowJD，ChenKJ. 2001.PhasecorrectioninseparatingP-andS-wavesinelasticdata[J]. Geophysics， 66(5): 1515--1518.

SunR，McMechanGA. 2001.Scalarreverse-timedepthmigrationofprestackelasticseismicdata[J]. Geophysics， 66(5): 1519--1527.

VirieuxJ. 1984.SH-wavepropagationinheterogeneousmedia:Velocity-stressfinite-differencemethod[J]. Geophysics， 49(11): 1933--1942.

VirieuxJ. 1986.P-SVwavepropagationinheterogeneousmedia:Velocity-stressfinite-differencemethod[J]. Geophy-sics， 51(4): 889--901.

YanJ，SavaP. 2009a.Elasticwave-modeseparationforVTImedia[J]. Geophysics， 74(5):WB19--WB32.

YanJ，SavaP. 2009b.Elasticwave-modeseparationforTTImedia[C]∥Proceedings of the 79th Annual International Meeting， SEG， Expanded Abstracts.Houston:SocietyofExplorationGeophysicists: 1197--1021.

ZhangQS，McMechanGA. 2010. 2Dand3DelasticwavefieldvectordecompositioninthewavenumberdomainforVTImedia[J]. Geophysics， 75(3):D13--D26.

Effects of staggered-grid and rotated staggered-grid on the wavefield separation in VTI media

1)SchoolofGeoscience，ChinaUniversityofPetroleum，ShandongQingdao266580，China2)ChinaOilfieldServicesLimited，Tianjin300451，China

Wavefield separation is an important part of multi-component seismic data processing， and the wavefield numerical simulation method is one of the necessary tools to this process. When staggered-grid finite difference scheme and rotated staggered-grid finite difference scheme are utilized to elastic wavefield simulation in anisotropic media， the definition of the particle velocity components and that of the stress tensor are different among these two grids. Furthermore， these two types of definitions will affect the wavefield separation. So we take the case of transversely isotropic media with a vertical axis (VTI media) as an example to analyze. Firstly， we analyze the parameter definition of these two grids and details about the wavefield separation in VTI media. Then， we focus on the analysis on the separation effects caused by the different parameter definition of these two grids. Considering the aspects of wavefront continuation and the purity of wavefield separation， numerical experiments are utilized to verify the results. Finally we draw the conclusion that the parameter definition of rotated staggered-grid is more beneficial to the anisotropic numerical simulation and the wavefield separation.

wavefield separation； staggered-grid； rotated staggered-grid； finite difference； VTI media

國(guó)家自然科學(xué)基金(41174100)、 國(guó)家科技重大專項(xiàng)(2011ZX05019-008-08)和中國(guó)石油天然氣集團(tuán)公司(2014A-3609)聯(lián)合資助.

2015-07-16收到初稿， 2015-10-10決定采用修改稿.

e-mail: multi-wave@163.com

10.11939/jass.2016.01.011

P315.3+1

A

韓冬， 杜啟振， 張明強(qiáng)， 郭成鋒， 公緒飛. 2016. 交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析. 地震學(xué)報(bào), 38(1): 111--121. doi:10.11939/jass.2016.01.011.

Han D, Du Q Z, Zhang M Q, Guo C F, Gong X F. 2016. Effects of staggered-grid and rotated staggered-grid on the wavefield separation in VTI media.ActaSeismologicaSinica, 38(1): 111--121. doi:10.11939/jass.2016.01.011.