裂縫誘導(dǎo)雙相HTI介質(zhì)地震波場錯格偽譜法模擬與波場特征分析

劉財， 遲喚昭， 高煒， 鹿琪*， 蘭慧田

1 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院， 長春　130026 2 吉林大學(xué)地球科學(xué)學(xué)院， 長春　130061 3 大慶油田有限責(zé)任公司勘探開發(fā)研究院， 黑龍江大慶　163712

裂縫誘導(dǎo)雙相HTI介質(zhì)地震波場錯格偽譜法模擬與波場特征分析

劉財1,2， 遲喚昭2， 高煒1， 鹿琪1*， 蘭慧田3

1 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院， 長春130026 2 吉林大學(xué)地球科學(xué)學(xué)院， 長春130061 3 大慶油田有限責(zé)任公司勘探開發(fā)研究院， 黑龍江大慶163712

摘要裂縫誘導(dǎo)的雙相具有水平對稱軸的橫向各向同性(HTI)介質(zhì)模型是由一組平行排列的垂直裂縫嵌入到統(tǒng)計各向同性的流體飽和多孔隙巖石中而組成的，它綜合考慮了裂縫型儲層巖石的各向異性和孔隙性.高精度的地震波場數(shù)值模擬技術(shù)是研究該介質(zhì)中地震波傳播規(guī)律的主要方法.本文結(jié)合錯格偽譜法和時間分裂法，求解描述該介質(zhì)中地震波傳播的一階速度-應(yīng)力方程.模擬了單層和雙層模型中的地震波場，并對其進(jìn)行了特征分析.研究結(jié)果表明：錯格偽譜法能有效消除標(biāo)準(zhǔn)網(wǎng)格偽譜法波場模擬結(jié)果中出現(xiàn)的數(shù)值偽影現(xiàn)象，與時間分裂法結(jié)合能夠獲得穩(wěn)定的、高精度的模擬結(jié)果；裂縫誘導(dǎo)雙相HTI介質(zhì)中的地震波場兼具裂縫各向異性介質(zhì)和雙相介質(zhì)中傳播的地震波的波場特征.

關(guān)鍵詞裂縫誘導(dǎo)； 雙相HTI介質(zhì)； 錯格偽譜法； 時間分裂法； 波場特征

1引言

天然裂縫型油氣藏作為一類重要的油氣藏類型已經(jīng)引起了越來越多的關(guān)注.基于地震各向異性理論發(fā)展的相關(guān)技術(shù)是裂縫性儲層識別和預(yù)測的一類重要手段，而地震各向異性與裂縫型儲層參數(shù)的聯(lián)系是通過等效介質(zhì)理論模型來建立的，因而，等效介質(zhì)模型間接地對裂縫型油氣藏的勘探起到重要作用.理論模型與實(shí)際儲層介質(zhì)的接近程度，在一定程度上，決定了以此為理論基礎(chǔ)所發(fā)展起來的儲層預(yù)測和流體識別技術(shù)的有效性.傳統(tǒng)的裂縫介質(zhì)等效介質(zhì)模型基于單相介質(zhì)理論建立，主要包括Hudson模型(Hudson, 1981)和線性滑動模型(Schoenberg and Sayers, 1995)，它們都沒有考慮儲層巖石的孔隙性和滲透性.而實(shí)際上，裂縫型儲層巖石通常是既包含裂縫又包含孔隙的裂縫型孔隙介質(zhì)，而且，隨著研究的深入，人們逐漸認(rèn)識到裂縫型油氣藏中裂縫定向排列所引起的各向異性以及孔隙流體的存在所引起的雙相或多相性對地震波的傳播都有很大影響(Chapman, 2003; Gurevich, 2003; Gurevich et al., 2009; Sil et al., 2011).所以，建立能夠綜合考慮裂縫型儲層巖石的各向異性和孔隙型的等效介質(zhì)理論模型對于裂縫型儲層的精細(xì)勘探尤為重要.

近年來，研究者們提出了一些綜合考慮裂縫各向異性、孔隙性以及裂縫與孔隙之間流體交換作用的介質(zhì)模型，并對彈性波的傳播特征進(jìn)行了研究.Thomsen(1995)、Gurevich(2003)假設(shè)裂縫與球形等徑孔隙之間流體壓力平衡，提出了考慮孔隙和裂縫液壓連通的裂縫介質(zhì)模型；Sil等(2011)、唐杰等(2015)將Gurevich模型應(yīng)用于多孔隙HTI介質(zhì)的流體替換分析中.Hudson等(1996, 2001)考慮裂縫與孔隙之間的流體壓力不完全平衡，推導(dǎo)了與頻率相關(guān)的動態(tài)等效剛度張量.Chapman(2003)、Jakobsen等(2003)和Gurevich等(2009)提出了含有中觀尺度(遠(yuǎn)大于孔隙尺寸，同時又小于彈性波波長)裂縫的裂縫介質(zhì)模型，并通過波誘導(dǎo)的流體中觀流動機(jī)制解釋了一些裂縫型儲層中地震頻段內(nèi)所觀測到的頻散和衰減現(xiàn)象.Parra(2000)結(jié)合BISQ模型(Dvorkin and Nur, 1993)和Thomsen模型，提出了一種考慮微裂縫與等徑孔隙之間流體局部噴射流動的等效介質(zhì)模型，并研究了滲透率各向異性對地震波頻散和衰減的影響.結(jié)合BISQ模型和Hudson模型，軒義華等(2006)、張顯文等(2010)分別建立了雙相HTI介質(zhì)和雙相裂隙正交各向異性介質(zhì)模型.劉財?shù)?2013)結(jié)合改進(jìn)的BISQ模型和Gurevich模型提出了又一種考慮微裂縫與等徑孔隙之間流體局部噴射流動的等效介質(zhì)模型，并采用偽譜法進(jìn)行了波場數(shù)值模擬.杜啟振等(2009)、Du等(2012)、孔麗云等(2012)結(jié)合孔隙彈性Biot模型和線性滑動模型建立了裂縫誘導(dǎo)HTI雙孔隙介質(zhì)模型，并進(jìn)行了波場數(shù)值模擬與裂縫參數(shù)分析.高煒等(2014)結(jié)合Biot模型和Gurevich模型建立了一種裂縫誘導(dǎo)雙相HTI介質(zhì)模型，給出了彈性波的傳播方程.不同模型在地震頻段的適用性主要取決于巖石中流體的流動性(Batzle et al., 2006).

波動方程數(shù)值模擬是認(rèn)識裂縫型儲層中地震波傳播規(guī)律的重要手段.在諸多波動方程數(shù)值解法中，偽譜法以其精度高、壓制網(wǎng)格數(shù)值頻散效果好的特點(diǎn)，被廣泛應(yīng)用于復(fù)雜介質(zhì)的波場數(shù)值模擬中(劉洋和李承楚, 2000; 郭智奇等, 2007; 劉財?shù)? 2007; 單啟銅和樂友喜， 2007；巴晶等, 2008; Ba et al., 2008；張軍舵等, 2008； 李紅星和陶春輝, 2009; Carcione and Gurevich, 2011)，尤其適用于黏彈性介質(zhì)、孔隙介質(zhì)等耗散介質(zhì)的波場模擬問題，因為在耗散介質(zhì)中傳播的彈性波會因能量衰減而伴隨產(chǎn)生速度頻散，波場數(shù)值模擬中采用偽譜法則可最大程度地避免任何頻段內(nèi)這種與介質(zhì)性質(zhì)相關(guān)的物理頻散與數(shù)值計算帶來的頻散的混淆.然而，研究發(fā)現(xiàn)，采用偽譜法求解一階速度-應(yīng)力方程時，存在Nyquist誤差問題，導(dǎo)致波場模擬結(jié)果中出現(xiàn)數(shù)值偽影現(xiàn)象.Witte和Richards(1987)提出了一種交錯網(wǎng)格偽譜法(簡稱錯格偽譜法)來克服這一問題.?zdenvar和McMechan(1996, 1997)詳細(xì)分析了Nyquist誤差問題，從數(shù)學(xué)上論證了錯格偽譜法較標(biāo)準(zhǔn)網(wǎng)格偽譜法的改進(jìn)，并將其應(yīng)用到標(biāo)量方程、聲波方程、彈性波方程和孔隙彈性方程的數(shù)值求解中.之后，該方法被其他一些學(xué)者應(yīng)用于波場模擬中，例如：Chen(1996)采用錯格偽譜法進(jìn)行了黏聲介質(zhì)中波場模擬；Carcione等(1999)對黏彈各向異性介質(zhì)和孔隙黏彈性介質(zhì)中的波場進(jìn)行了數(shù)值模擬；劉炯等(2008;Liu et al., 2010)模擬了縫洞型儲層模型、隨機(jī)介質(zhì)中的地震波場；巴晶等(2010)模擬了礁、灘相儲層的地震波場；杜增利等(2010)、吳寶年等(2012)模擬了聲學(xué)介質(zhì)和彈性介質(zhì)中的地震波場.

與前人的模型相比(杜啟振等, 2009; Du et al., 2012; 孔麗云等, 2012)，裂縫誘導(dǎo)雙相HTI介質(zhì)模型具有物理意義明確、彈性參數(shù)計算更為簡潔的優(yōu)點(diǎn).本文對裂縫誘導(dǎo)雙相HTI介質(zhì)中的地震波場進(jìn)行研究，為了解決該類介質(zhì)中彈性波傳播一階速度-應(yīng)力方程的剛性問題并獲得高精度的波場模擬結(jié)果，將錯格偽譜法與時間分裂法相結(jié)合對波傳播方程進(jìn)行數(shù)值求解，并對模擬波場進(jìn)行特征分析，為進(jìn)一步深入認(rèn)識實(shí)際裂縫型儲層介質(zhì)的地震波傳播規(guī)律奠定理論基礎(chǔ).

2裂縫誘導(dǎo)雙相HTI介質(zhì)模型及彈性波傳播方程

實(shí)際裂縫型儲層巖石通常既包含裂縫又包含孔隙，前者孔隙度很小(通常<0.1%)主要影響儲層的滲透性，而后者主要存儲儲層中的流體.在波作用下，裂縫與孔隙之間產(chǎn)生流體壓力梯度，流體發(fā)生流動以使壓力達(dá)到新的平衡，這一過程對波的傳播有重要影響.對于某些高滲透性裂隙巖石，在地震頻段內(nèi)，流體壓力有充足時間在裂縫與孔隙之間達(dá)到平衡.基于這一考慮，杜啟振等(2009, 2012)將一組垂直定向排列的裂縫嵌入到流體飽和各向同性孔隙背景介質(zhì)中，建立了裂縫誘導(dǎo)HTI雙孔隙介質(zhì)模型.采用類似的方式，高煒等(2014)建立了裂縫誘導(dǎo)雙相HTI介質(zhì)模型(圖1).

圖1　裂縫誘導(dǎo)雙相HTI介質(zhì)模型示意Fig.1　Diagram of fracture-induced two-phase HTI medium

該模型是基于以下假設(shè)建立的(高煒, 2014)： ① 干燥的裂縫型孔隙巖石由一組平行排列的垂直裂縫嵌入到滿足統(tǒng)計各向同性的孔隙背景巖石中構(gòu)成，孔隙之間以及孔隙與裂縫之間相互連通，裂縫和孔隙的最大尺寸遠(yuǎn)小于彈性波波長； ② 孔隙和裂縫由一種各向同性、具有黏滯性和可壓縮性的流體所飽和； ③ 在波傳播方向上，固體骨架和流體之間存在相對位移，流體相對于固體的流動屬于Poiseuille型流動，流體滲流遵守達(dá)西定律；孔隙和裂縫之間的流體壓力在波的半個周期內(nèi)有足夠時間達(dá)到完全平衡； ④ 忽略熱彈性效應(yīng)，孔隙流體與巖石基質(zhì)也不發(fā)生化學(xué)作用.

基于以上假設(shè)，二維裂縫誘導(dǎo)雙相HTI介質(zhì)一階速度-應(yīng)力彈性波傳播方程為(高煒,2014)：

(1a)

(1b)

(1c)

(1d)

(1e)

(1f)

(1g)

(1h)

式中，

3波傳播方程的數(shù)值解法

3.1剛性問題

由于雙相介質(zhì)中存在傳播速度與單相介質(zhì)相似的快縱波和傳播速度非常慢的慢縱波，傳播矩陣具有兩種大小相差很大的不同特征值，表明一階速度-應(yīng)力微分方程組是剛性的.在數(shù)值模擬時需要采用非常小的時間步長才能滿足穩(wěn)定性條件，這將大大增加計算量.為了解決剛性問題，以便可以采用大步長的顯式時間積分法對其進(jìn)行數(shù)值求解，Carcione和Quiroga-Goode(1995)引入了一種時間分裂解法.在時間分裂法中，剛性一階偏微分方程組被分成兩部分，一部分為剛性，另一部分為非剛性，交替求解兩個方程組，其中一個方程組的解作為另一個方程組解的初始值.

根據(jù)時間分裂法，可以將方程(1)分裂為剛性方程組和非剛性方程組，其中剛性方程組為

(2a)

(2b)

(2c)

(2d)

(3a)

(3b)

(3c)

(3d)

即為方程(1)式的數(shù)值解.

(4a)

(4b)

(4c)

(4d)

(4e)

(4f)

(4g)

(4h)

3.2錯格偽譜法離散遞推格式

對波傳播平面的空間區(qū)域和傳播時間進(jìn)行離散化，取x=iΔx或x=(i±1/2)Δx，z=jΔz或z=(j±1/2)Δz，t=nΔt或t=(n±1/2)Δt，其中Δx、Δz和Δt分別表示空間離散采樣步長和時間采樣步長，i,j為整數(shù)，分別表示x和z方向的空間離散采樣點(diǎn)號，n表示離散時間采樣點(diǎn)號.場分量(速度或應(yīng)力分量)和介質(zhì)物性參數(shù)以圖2的方式分配在整網(wǎng)格點(diǎn)或半網(wǎng)格點(diǎn)上.半網(wǎng)格點(diǎn)處的介質(zhì)物性參數(shù)通過整網(wǎng)格點(diǎn)處的平均值來計算，則(i+1/2,j)和(i,j+1/2)處的ρ,ρf,T,φ和η/k值分別為

這里a代表上述介質(zhì)物性參數(shù).而在(i+1/2,j+1/2)處的c55的值為

場量的一階空間微分由錯格傅里葉偽譜微分算子計算，若以φ表示場量，則其沿x1方向的一階空間微分為

(5)

因此，非剛性方程組(4)式的錯格偽譜法離散格式為

圖2　二維裂縫誘導(dǎo)雙相HTI介質(zhì)速度-應(yīng)力方程交錯網(wǎng)格示意Fig.2　Diagram of velocity-stress equation staggered grid in 2D fracture-induced two-phase HTI medium

(6a)

(6b)

(6c)

(6d)

(6e)

(6f)

(6g)

4數(shù)值算例與分析

4.1單層模型

設(shè)計一個單層裂縫誘導(dǎo)雙相HTI介質(zhì)模型，模型網(wǎng)格數(shù)為255×255，網(wǎng)格大小為10 m×10 m.震源采用流體體積源加垂直分量固相源，即將震源能量加載在流體壓力和固相垂直分量上.震源中心點(diǎn)位置為(1270 m, 1270 m)，震源時間函數(shù)采用Ricker子波，震源函數(shù)形式如下：

×e-λ[(x-x0)2+(z-z0)2],

(7)

式中,f0為震源中心頻率(震源主頻)，t0=1/f0為震源子波延遲時間，(x0,z0)表示震源中心位置，λ表示震源力空間作用的集中系數(shù).震源子波主頻為40 Hz，延遲時間為0.025 s.時間采樣間隔為1 ms.

為消除人工邊界反射，采用簡便實(shí)用的Cerjan衰減邊界(Cerjan et al.， 1985)，即沿人工邊界向外擴(kuò)充N個網(wǎng)格點(diǎn)，進(jìn)入擴(kuò)充區(qū)域內(nèi)的波場通過乘以下面的因子逐漸衰減為0，其中a為衰減系數(shù)，可通過試驗確定最佳值.

(8)

由于所建立的裂縫誘導(dǎo)雙相HTI介質(zhì)模型考慮了流體的全局流動，根據(jù)Biot理論，其波場中可能會包含有慢縱波，但雙相介質(zhì)中流-固界面的黏滯摩擦阻礙流-固相對運(yùn)動，將導(dǎo)致慢縱波快速衰減.為了研究慢縱波的傳播特征，本文除進(jìn)行含標(biāo)準(zhǔn)黏滯流體的黏滯相界情況下的彈性波傳播數(shù)值模擬外，也開展流體為黏滯系數(shù)非常小的近似理想相界情況下的波場模擬.兩種情況下的數(shù)值模擬所采用的介質(zhì)物性參數(shù)除黏滯系數(shù)η不同以外，其余物性參數(shù)相同，見表1，其中Ks、μs和Kf分別表示固體基質(zhì)的體積模量、剪切模量和流體體積模量，ΔN和ΔT表示裂縫法向和切向弱度參數(shù)，φp和φc表示背景巖石孔隙度和裂縫孔隙度，其余參數(shù)物理意義同前文.在模擬近似理想相界時，取η=10-10Pa·s，而在模擬黏滯相界時，取水的黏滯系數(shù)η=0.001 Pa·s.圖3給出了近似理想相界情況下，固體和流體(相對于固體)質(zhì)點(diǎn)速度垂直分量和水平分量在350 ms時的波場傳播快照.相應(yīng)地，圖4給出了黏滯相界情況下的波場傳播快照.圖5是近似理想相界情況下，位于(108，108)網(wǎng)格點(diǎn)處的質(zhì)點(diǎn)速度垂直分量時間記錄圖，圖6是黏滯相界情況下，該質(zhì)點(diǎn)速度垂直分量時間記錄圖.

表1　裂縫誘導(dǎo)雙相HTI介質(zhì)物性參數(shù)

從圖3和圖4中各分量的波場快照可以看出，波場中各波前面都非常清晰，沒有可見網(wǎng)格數(shù)值頻散的產(chǎn)生，表明了錯格偽譜法在裂縫誘導(dǎo)雙相HTI介質(zhì)波場數(shù)值模擬中的有效性和精確性.

分析以上波場數(shù)值模擬結(jié)果可知：(1)震源在裂縫誘導(dǎo)雙相HTI介質(zhì)中激發(fā)了三種類型的波，從外向內(nèi)分別為快準(zhǔn)P波、準(zhǔn)SV波和慢P波；(2)快準(zhǔn)P波和準(zhǔn)SV波類似于單相各向異性介質(zhì)中的準(zhǔn)縱波和準(zhǔn)SV波，波前面呈橢圓形，準(zhǔn)SV波的波前面出現(xiàn)波面尖角現(xiàn)象，這體現(xiàn)了裂縫誘導(dǎo)雙相HTI介質(zhì)中裂縫各向異性所引起的波場特征；(3)慢準(zhǔn)P波在不同相界情況下具有不同的特征，當(dāng)流體為近似理想相界時(圖3)，慢縱波呈傳播模式，波前面也呈橢圓向外傳播，其傳播速度明顯低于快準(zhǔn)P波和準(zhǔn)SV波，而當(dāng)流體為黏滯性時(圖4)，慢縱波表現(xiàn)為擴(kuò)散性，以靜態(tài)模式出現(xiàn)在震源位置，如圖4中箭頭所標(biāo)示的慢P波；(4)結(jié)合圖3和圖4的波場快照以及圖5和圖6的質(zhì)點(diǎn)速度垂直分量的時間記錄可以看出：快準(zhǔn)P波和準(zhǔn)SV波引起固相質(zhì)點(diǎn)和流相質(zhì)點(diǎn)的同相位振動，慢準(zhǔn)P波引起固相質(zhì)點(diǎn)和流相質(zhì)點(diǎn)反相位振動；近似理想相界情況下，快準(zhǔn)P波和準(zhǔn)SV波在固相中的振幅明顯大于在流相中的振幅，慢準(zhǔn)P波在固相中的振幅明顯小于在流相中的振幅，即在流相中更容易觀測到慢縱波；在黏滯相界情況下，快準(zhǔn)P波和準(zhǔn)SV波在固相和流相中振幅差異非常小，慢準(zhǔn)P波在固相和流相中均呈擴(kuò)散性，在震源附近快速衰減.

圖3　裂縫誘導(dǎo)雙相HTI介質(zhì)中固相(A)和流相(相對于固相)(B)質(zhì)點(diǎn)速度水平分量(a1,b1)和垂直分量(a2,b2)在350 ms時的波場傳播快照(η=10-10 Pa·s)Fig.3　Wavefield snapshots of solid phase (A) and fluid phase (relative to solid phase) (B) particle velocity horizontal component (a1,b1) and vertical component (a2,b2) in fracture-induced two-phase HTI medium at 350 ms (η=10-10 Pa·s)

圖4　裂縫誘導(dǎo)雙相HTI介質(zhì)中固相(A)和流相(相對于固相)(B)質(zhì)點(diǎn)速度水平分量(a1,b1)和垂直分量(a2,b2)在350 ms時的波場傳播快照(η=0.001 Pa·s)Fig.4　Wavefield snapshots of solid phase (A) and fluid phase (relative to solid phase) (B) particle velocity horizontal component (a1,b1) and vertical component (a2,b2) in fracture-induced two-phase HTI medium at 350 ms (η=0.001 Pa·s)

三種波的波前面?zhèn)鞑ヌ卣鳌⒐滔嗪土飨噘|(zhì)點(diǎn)的振動特性以及慢準(zhǔn)P波的存在共同表明：裂縫誘導(dǎo)雙相HTI介質(zhì)具有裂縫各向異性和雙相性，模型中背景孔隙介質(zhì)主導(dǎo)了介質(zhì)的雙相性，而裂縫系統(tǒng)則誘導(dǎo)了介質(zhì)的各向異性.

為了分析不同流體類型對裂縫誘導(dǎo)雙相HTI介質(zhì)中波傳播特征的影響，將上述模型中的流體由水替換為氣，介質(zhì)模型的物性參數(shù)分別為：流體體積模量Kf=0.025 GPa、流體密度ρf=65 kg·m-3、流體黏滯系數(shù)η=2×10-5Pa·s，其余參數(shù)同表1.圖7給出了裂縫誘導(dǎo)雙相HTI介質(zhì)含水和含氣兩種情況下的固相質(zhì)點(diǎn)速度垂直分量的波場快照，圖8給出了固相質(zhì)點(diǎn)(108,108)速度垂直分量的時間記錄.從圖中可以看出，快準(zhǔn)P波在含氣飽和裂縫誘導(dǎo)雙相HTI介質(zhì)中的傳播速度比含水時的傳播速度小，準(zhǔn)SV波則與之相反.快準(zhǔn)P波和準(zhǔn)SV波在含不同流體的裂縫誘導(dǎo)雙相HTI介質(zhì)中傳播時的速度大小關(guān)系，符合Gassmann方程，對于準(zhǔn)P波，速度的大小由含流體多孔隙巖石的彈性模量和密度共同決定，因而，流體體積模量和密度共同影響地震波的傳播速度，對于低頻范圍，含不同流體的巖石，

圖5　近似理想相界情況固相(實(shí)線)和流相(虛線)質(zhì)點(diǎn)(108，108)速度垂直分量時間記錄(附圖為固相質(zhì)點(diǎn)速度分量時間記錄中慢準(zhǔn)P波的局部放大圖)Fig.5　Vertical component time record of solid phase(solid line) and fluid phase(dashed line) particle(108,108) in the ideal phase boundary case (top right corner is the enlarged figure of slow quasi-P wave in solid phase particle velocity component time record)

圖6　黏滯相界情況固相(實(shí)線)和流相(虛線)質(zhì)點(diǎn)(108，108)速度垂直分量時間記錄Fig.6　Vertical component time record of solid phase (solid line) and fluid phase(dashed line)particle (108,108) in the viscous phase boundary case

圖7　裂縫誘導(dǎo)雙相HTI介質(zhì)含水(a)和含氣(b)時固相質(zhì)點(diǎn)速度垂直分量350 ms波場快照Fig.7　Wavefield snapshots of solid phase particle velocity vertical component in fracture-induced two-phase HTI medium at 350 ms when the pore fluid is water (a) and gas (b)

其體積模量的差異比密度的差異對速度的影響更大，所以含水時速度大于含氣時；而對于準(zhǔn)SV波，不同流體密度的差異對速度的影響更大，所以，含氣時速度大于含水時.另外，由于兩種流體均是黏滯性流體，慢準(zhǔn)P波表現(xiàn)出彌散性，快速衰減，波場中無法觀測到.

4.2雙層模型

通過一個兩層模型來分析彈性波在雙層裂縫誘導(dǎo)雙相HTI介質(zhì)分界面處的反射和透射特征.模型網(wǎng)格數(shù)為255×255，網(wǎng)格大小為10 m×10 m.界面位于1680 m處.上、下層均為裂縫誘導(dǎo)雙相HTI介質(zhì)，物性參數(shù)見表2.震源采用流體體積注入源加垂直分量固相源，震源中心點(diǎn)位置為(1270 m,1480 m)，震源函數(shù)和震源參數(shù)同上.時間采樣間隔為0.5 ms.接收排列在深度1180 m處水平分布.為了研究慢縱波在介質(zhì)分界面處的反射和透射特征，假設(shè)上層介質(zhì)中含有非黏滯性流體，即上層介質(zhì)為理想相界情況.對于下層介質(zhì)，考慮含非黏滯性和黏滯性流體兩種情況，即除考慮表2中下層介質(zhì)黏滯系數(shù)η=0外，還考慮黏滯系數(shù)η=0.001的情況.

圖8　裂縫誘導(dǎo)雙相HTI介質(zhì)含水(虛線)和含氣(實(shí)線)時固相質(zhì)點(diǎn)(108,108)速度垂直分量時間記錄Fig.8　Vertical component time record of solid phase particle(108,108) in fracture-induced two-phase HTI medium when the pore fluid is water(dotted line) and gas(solid line)

為了驗證錯格偽譜法在強(qiáng)波阻抗界面存在情況下波場模擬的優(yōu)勢，所設(shè)計的模型界面兩側(cè)波阻抗差較大.分別采用標(biāo)準(zhǔn)網(wǎng)格偽譜法和錯格偽譜法進(jìn)行波場模擬，并將兩種方法的模擬結(jié)果進(jìn)行對比分析.

表2　裂縫誘導(dǎo)雙相HTI介質(zhì)雙層模型的物性參數(shù)

圖9和圖10分別是采用標(biāo)準(zhǔn)網(wǎng)格偽譜法和錯格偽譜法進(jìn)行波場模擬所得到的裂縫誘導(dǎo)雙相HTI介質(zhì)固相和流相質(zhì)點(diǎn)速度分量的波場快照.從圖9可以看出，在標(biāo)準(zhǔn)網(wǎng)格偽譜法模擬得到的波場快照中存在明顯的振鈴?fù)衔铂F(xiàn)象，主要集中在垂直方向上，而在圖10所給出的錯格偽譜法波場模擬結(jié)果中則不存在著這樣的現(xiàn)象.這表明，采用錯格偽譜法能夠適應(yīng)地下存在強(qiáng)波阻抗反射界面的波場數(shù)值模擬問題，與標(biāo)準(zhǔn)網(wǎng)格偽譜法相比，能夠獲得更高的精度.

圖11和圖12分別給出了當(dāng)下層介質(zhì)為理想相界時，雙層裂縫誘導(dǎo)雙相HTI介質(zhì)波場模擬的波場快照和合成地震記錄，而圖13則給出了下層介質(zhì)為黏滯相界時的波場快照.

從圖11可知：當(dāng)兩層裂縫誘導(dǎo)雙相HTI介質(zhì)中的流體均為非粘滯性時，快準(zhǔn)P波、準(zhǔn)SV波和慢準(zhǔn)P波在介質(zhì)分界面處均發(fā)生了反射和透射現(xiàn)象，并發(fā)生了波型轉(zhuǎn)換，波場快照中可同時觀測到多種波，這些波大致可以分為以下4類：(1)直達(dá)波類：直達(dá)快準(zhǔn)P波、直達(dá)準(zhǔn)SV波和直達(dá)慢準(zhǔn)P波；(2)反射波類：反射快準(zhǔn)P波、反射準(zhǔn)SV波、反射慢準(zhǔn)P波；(3)透射波類：透射快準(zhǔn)P波、透射準(zhǔn)SV波、透射慢準(zhǔn)P波；(4)轉(zhuǎn)換波類：由快準(zhǔn)P波在界面反射之后形成的反射轉(zhuǎn)換準(zhǔn)SV波和慢準(zhǔn)P波，由準(zhǔn)SV波在界面折射之后形成的折射轉(zhuǎn)換快準(zhǔn)P波，由準(zhǔn)SV波在界面反射之后形成的反射轉(zhuǎn)換慢準(zhǔn)P波，由慢準(zhǔn)P波在界面反射之后形成的反射轉(zhuǎn)換快準(zhǔn)P波和準(zhǔn)SV波，由快準(zhǔn)P波在界面透射之后形成的透射轉(zhuǎn)換準(zhǔn)SV波和慢準(zhǔn)P波，由準(zhǔn)SV波在界面透射之后形成的透射轉(zhuǎn)換快準(zhǔn)P波和慢準(zhǔn)P波，由慢準(zhǔn)P波在界面透射之后形成的透射轉(zhuǎn)換快準(zhǔn)P波和準(zhǔn)SV波.在合成地震記錄中可以觀測到反射域的各種彈性波，如圖12所示.從圖13中可以看出，當(dāng)下層介質(zhì)中的流體為粘滯性時，透射的慢準(zhǔn)P波以及由快準(zhǔn)P波和準(zhǔn)SV波形成的透射轉(zhuǎn)換慢準(zhǔn)P波由于強(qiáng)衰減已觀測不到.

5結(jié)論與認(rèn)識

本文采用錯格偽譜法和時間分裂法求解裂縫誘導(dǎo)雙相HTI介質(zhì)的一階速度-應(yīng)力方程進(jìn)行地震波場數(shù)值模擬，以此研究了裂縫誘導(dǎo)雙相HTI介質(zhì)單層模型和雙層介質(zhì)模型中的地震波傳播規(guī)律，取得了以下結(jié)論和認(rèn)識：

(1) 錯格偽譜法能夠有效消除標(biāo)準(zhǔn)網(wǎng)格偽譜法波場模擬結(jié)果中存在的數(shù)值偽影現(xiàn)象，與時間分裂法相結(jié)合，能夠獲得穩(wěn)定的、高精度的波場模擬結(jié)果，是一種進(jìn)行裂縫誘導(dǎo)雙相HTI介質(zhì)中波場數(shù)值模擬的高精度算法.

圖9　標(biāo)準(zhǔn)網(wǎng)格偽譜法波場模擬所得的固相(A)和流相(B)質(zhì)點(diǎn)速度水平分量(a1,b1)和垂直分量(a2,b2)在500 ms時的波場傳播快照Fig.9　Wavefield snapshots of solid phase (A) and fluid phase (B) particle velocity horizontal component (a1,b1) and vertical component (a2,b2) in standard grid pseudo-spectral method at 500 ms

圖10　錯格偽譜法波場模擬所得的固相(A)和流相(B)質(zhì)點(diǎn)速度水平分量(a1,b1)和垂直分量(a2,b2)在500 ms時的波場傳播快照Fig.10　Wavefield snapshots of solid phase (A) and fluid phase (B) particle velocity horizontal component (a1,b1) and vertical component (a2,b2) in staggered-grid pseudo-spectral method at 500 ms

圖11　固相(a1,b1)和流相(相對于固相)(a2,b2)質(zhì)點(diǎn)速度水平分量(B)和垂直分量(A)在500 ms時的波場傳播快照(下層介質(zhì)為理想相界)1 直達(dá)快準(zhǔn)P波； 2 直達(dá)準(zhǔn)SV波； 3 直達(dá)慢準(zhǔn)P波； 4 反射快準(zhǔn)P波； 5 反射準(zhǔn)SV波； 6 反射慢準(zhǔn)P波； 7 透射快準(zhǔn)P波； 8 透射準(zhǔn)SV波； 9 透射慢準(zhǔn)P波； 10 快準(zhǔn)P波產(chǎn)生的反射轉(zhuǎn)換準(zhǔn)SV波； 11 快準(zhǔn)P波產(chǎn)生的反射轉(zhuǎn)換慢準(zhǔn)P波12 準(zhǔn)SV波產(chǎn)生的折射轉(zhuǎn)換快準(zhǔn)P波； 13 準(zhǔn)SV波產(chǎn)生的反射轉(zhuǎn)換慢準(zhǔn)P波； 14 慢準(zhǔn)P波產(chǎn)生的反射轉(zhuǎn)換快準(zhǔn)P波； 15 慢準(zhǔn)P波產(chǎn)生的反射轉(zhuǎn)換準(zhǔn)SV波； 16 快準(zhǔn)P波產(chǎn)生的透射轉(zhuǎn)換準(zhǔn)SV； 17 快準(zhǔn)P波產(chǎn)生的透射轉(zhuǎn)換慢準(zhǔn)P波； 18 準(zhǔn)SV波產(chǎn)生的透射轉(zhuǎn)換快準(zhǔn)P波； 19 準(zhǔn)SV波產(chǎn)生的透射轉(zhuǎn)換慢準(zhǔn)P波； 20 慢準(zhǔn)P波產(chǎn)生的透射轉(zhuǎn)換快準(zhǔn)P波； 21 慢準(zhǔn)P波產(chǎn)生的透射轉(zhuǎn)換準(zhǔn)SV波.Fig.11　Wavefield snapshots of solid phase (left) and fluid phase (right) particle velocity horizontal component (down) and vertical component (up) at 500 ms (the lower layer medium is the ideal phase boundary)1 Direct fast quasi-P wave (qP1 wave)； 2 Direct quasi-SV wave (qSV wave)； 3 Direct slow quasi-P wave (qP2 wave)； 4 Reflected qP1 wave； 5 Reflected qSV wave； 6 Reflected qP2 wave； 7 Transmitted qP1 wave； 8 Transmitted qSV wave； 9 Transmitted qP2 wave； 10 Reflected converted quasi-SV wave by fast quasi-P wave (qP1-qSV)； 11 Refracted qSV-qP1； 12 Reflected qSV-qP1； 13 Reflected qSV-qP2； 14 Reflected qP2-qP1； 15 Reflected qP2-qSV； 16 Transmitted qP1-qSV； 17 Transmitted qP1-qP2； 18 Transmitted qSV-qP1； 19 Transmitted qSV-qP2； 20 Transmitted qP2-qP1, 21 Transmitted qP2-qSV.

圖12　固相(a1,b1)和流相(相對于固相)(a2,b2)速度水平分量(B)和垂直分量(A)合成地震記錄Fig.12　Synthetic seismogram of solid phase (relative to solid phase) (a1,b1) and fluid phase (a2,b2)velocity horizontal component (B) and vertical component (A)

圖13　固相(a1,b1)和流相(相對于固相)(b2,b2)質(zhì)點(diǎn)速度水平分量(B)和垂直分量(A)在500 ms時的波場傳播快照(下層介質(zhì)為黏滯相界)Fig.13　Wavefield snapshots of solid phase (a1,b1) and fluid phase (a2,b2) particle velocity horizontal component (B) and vertical component (A) at 500 ms (the lower layer medium is the viscous phase boundary)

(2) 快準(zhǔn)P波、準(zhǔn)SV波和慢準(zhǔn)P波在近似理想相界和黏滯相界情況下的波前面?zhèn)鞑ヌ卣骱唾|(zhì)點(diǎn)振動特性表明，裂縫誘導(dǎo)雙相HTI介質(zhì)中的地震波場兼具裂縫各向異性介質(zhì)和雙相介質(zhì)中傳播的地震波的波場特征.

(3) 各種波在介質(zhì)分界面處均發(fā)生反射和透射現(xiàn)象，并發(fā)生波型轉(zhuǎn)換，波場快照中可同時觀測到多種波，這些波可以分為直達(dá)波、反射波、透射波和轉(zhuǎn)換波4類，快準(zhǔn)P波和準(zhǔn)SV波與慢P波之間可發(fā)生相互轉(zhuǎn)換.對裂縫誘導(dǎo)雙相HTI介質(zhì)分界面處地震波反射和透射現(xiàn)象的研究，為進(jìn)一步加深對裂縫介質(zhì)儲層地震響應(yīng)的認(rèn)識奠定基礎(chǔ).

References

Ba J, Lu M H, Hu B, et al. 2008. The skeleton-relaxed model for poroviscoelastic media.ChineseJ.Geophys. (in Chinese), 51(5): 1527-1537, doi: 10.3321/j.issn:0001-5733.2008.05.027. Ba J, Nie J X, Cao H, et al. 2008. Mesoscopic fluid flow simulation in double-porosity rocks.GeophysicalResearchLetters, 35(4): L04303.

Ba J, Cao H, Yao F C, et al. 2010. Application of staggerred pseudospectrum method in simulating seismic wave field of reef beach facies reservoirs.OGP(in Chinese), 45(2): 177-184.

Batzle M L, Han D H, Hofmann R. 2006. Fluid mobility and frequency-dependent seismic velocity-Direct measurements.Geophysics, 71(1): N1-N9.

Carcione J M, Quiroga-Goode G. 1995. Some aspects of the physics and numerical modeling of Biot compressional waves.JournalofComputationalAcoustics, 3(34): 261-280. Carcione J M, Helle H B. 1999. Numerical solution of the poroviscoelastic wave equation on a staggered mesh.J.Comput.Phys., 154(2): 520-527. Carcione J M. 1999. Staggered mesh for the anisotropic and viscoelastic wave equation.Geophysics, 64(6): 1863-1866.

Carcione, J M., Gurevich, B., 2011. Differential form and numerical implementation of Biot′s poroelasticity equations with squirt dissipation.Geophysics, 76(6),: N55-N64.

Cerjan C, Kosloff D, Kosloff R, et al. 1985. A nonreflecting boundary condition for discrete acoustic and elastic wave equation.Geophysics, 50(4): 705-708.

Chapman M. 2003. Frequency-dependent anisotropy due to meso-scale fractures in the presence of equant porosity.GeophysicalProspecting, 51(5): 369-379.

Chen H W. 1996. Staggered-grid pseudospectral viscoacoustic wave field simulation in two-dimensional media.J.Acoust.Soc.Am., 100(1): 120-131. Du Q Z, Kong L Y, Han S C. 2009. Wavefield propagation characteristics in the fracture-induced anisotropic double-porosity medium.ChineseJ.Geophys. (in Chinese), 52(4): 1049-1058, doi: 10.3969/j.issn.0001-5733.2009.04.022.

Du Q Z, Wang X M, Ba J, et al. 2012. An equivalent medium model for wave simulation in fractured porous rocks.GeophysicalProspecting, 60(5): 940-956. Du Z L, Xu F, Gao H L. 2010. Pseudo spectral seismic wavefield simulation with staggered grid.GPP(in Chinese), 49(5): 430-437. Dvorkin J, Nur A. 1993. Dynamic poroelasticity: A unified model with the squirt and the Biot mechanisms.Geophysics, 58(4): 524-533.

Gao W, Liu C, Guo Z Q, et al. 2014. Fracture-induced two-phase HTI medium model and its elastic wave propagation equations.GlobalGeology(in Chinese), 33(4): 904-909.

Guo Z Q, Liu C, Yang B J, et al. 2007. Seismic wave fields modeling and feature in viscoelastic anisotropic media.ProgressinGeophysics(in Chinese), 22(3): 804-810. Gurevich B. 2003. Elastic properties of saturated porous rocks with aligned fractures.JournalofAppliedGeophysics, 54(3-4): 203-218. Gurevich B, Brajanovski R, Galvin R J, et al. 2009. P-wave dispersion and attenuation in fractured and porous reservoirs—poroelasticity approach.GeophysicalProspecting, 57(2): 225-237. Hudson J A. 1981. Wave speeds and attenuation of elastic waves in material containing cracks.GeophysicalJournaloftheRoyalAstronomicalSocietyInternational, 64(1): 133-150.

Hudson J A, Liu E, Crampin S. 1996. The mechanical properties of materials with interconnected cracks and pores.GeophysicalJournalInternational, 124(1): 105-112.

Hudson J A, Pointer T, Liu E. 2001. Effective-medium theories for fluid-saturated materials with aligned cracks.GeophysicalProspecting, 49(5): 509-522.

Jakobsen M, Johansen T A, McCann C. 2003. The acoustic signature of fluid flow in complex porous media.JournalofAppliedGeophysics, 54(3-4): 219-246.

Kong L Y, Wang Y B, Yang H Z. 2012. Fracture parameters analyses in fracture-induced HTI double-porosity medium.ChineseJ.Geophys. (in Chinese), 55(1): 189-196, doi: 10.6038/j.issn.0001-5733.2012.01.018.

Li H X, Tao C H. 2009. Features analysis of seismic wave field in two-phase anisotropic random medium with the pseudo-spectral method.ActaPhysicaSinica(in Chinese), 58(4): 2836-2842.

Liu C, Guo Z Q, Yang B J, et al. 2007. Analysis of reflection and transmission problems of waves in viscoelastic anisotropic media.ChineseJ.Geophys. (in Chinese), 50(4): 1216-1224.

Liu C, Lan H T, Guo Z Q, et al. 2013. Pseudo-spectral modeling and feature analysis of wave propagation in two-phase HTI medium based on reformulated BISQ mechanism.ChineseJ.Geophys. (in Chinese), 56(10): 3461-3473, doi: 10.6038/cjg20131021.Liu J, Ma J W, Yang H Z, et al. 2008. Staggered-grid pseudo-spectrum simulation of fracture and cave reservoir.OGP(in Chinese), 43(6): 723-727.

Liu J, Ba J, Ma J W, et al. 2010. An analysis of seismic attenuation in random porous media. Science.ChinaSerG-Phys.Mech.Astron., 53(4): 628-637.

Liu Y, Li C C. 2000. The study on elastic wave propagation with the pseudo-spectral method in two-phase anisotropic medium.ActaSeismologicaSinica(in Chinese), 22(2): 132-138.

?zdenvar T, McMechan G A. 1996. Causes and reduction of numerical arteifacts in pseudo-spectral wavefield extrapolation.Geophy.J.Int., 126(3): 819-828.

?zdenvar T, McMechan G A. 1997. Algorithms for staggered-grid computations for poroelastic, elastic, acoustic, and scalar wave equations.GeophysicalProspecting, 45(3): 403-420. Parra J O. 2000. Poroelastic model to relate seismic wave attenuation and dispersion to permeability anisotropy.Geophysics, 65(1): 202-210.Schoenberg M, Sayers C M. 1995. Seismic anisotropy of fractured rock.Geophysics, 60(1): 204-211.

Shan Q T, Yue Y X. 2007. Wavefield simulation of 2-D viscoelastic medium in Perfectly Matched Layer boundary.GPP(in Chinese), 46(2): 126-130.Sil S, Sen M K, Gurevich B. 2011. Analysis of fluid substitution in a porous and fractured medium.Geophysics, 76(3): WA157-WA166. Tang J, Fang B, Sun C Y, et al. 2015. Study of seismic wave propagation characteristics based on anisotropic fluid substitution in fractured medium.GeophysicalProspectingforPetroleum(in Chinese), 54(1): 1-8.

Thomsen L. 1995. Elastic anisotropy due to aligned cracks in porous rock.GeophysicalProspecting, 43(6): 805-829.

Witte D, Richards P G. 1987. Contribution to the pseudo-spectral method for computing synthetic seismograms.∥ Expanded Abstracts of 57th SEG Annual Int. Mtg., 517-519.

Wu B N, Wu S Q, Huang Y, et al. 2012. Application of staggered grids in elastic wavefield modeling by pseudo-spectral method.GPP(in Chinese), 51(5): 440-445.

Xuan Y H, He J D, Meng Q S, et al. 2006. Wavefield analysis of biphase EDA medium based on BISQ mechanism.OGP(in Chinese), 41(5): 550-556.

Zhang J D, Yue Y X, Wang Y X. 2008. Numerical simulation wave field by pseudo-spectrum method in isotropic two-phase media.GPP(in Chinese), 47(4): 338-345.

Zhang X W, Wang D L, Wang Z J, et al. 2010. The study on azimuth echaraecterlistiecs of attenuation and dispersion in 3D two-phase orthotropic crack medium based on BISQ mechanism.ChineseJ.Geophys. (in Chinese), 53(10): 2452-2459, doi: 10.3969/j.issn.0001-5733.2010.10.019.

附中文參考文獻(xiàn)

巴晶, 盧明輝, 胡彬等. 2008. 黏彈雙相介質(zhì)中的松弛骨架模型. 地球物理學(xué)報, 51(5): 1527-1537, doi: 10.3321/j.issn:0001-5733.2008.05.027.

巴晶, 曹宏, 姚逢昌等. 2010. 利用錯格虛譜法模擬礁、灘相儲層地震波場. 石油地球物理勘探, 45(2): 177-184.

杜啟振, 孔麗云, 韓世春. 2009. 裂縫誘導(dǎo)各向異性雙孔隙介質(zhì)波場傳播特征. 地球物理學(xué)報, 52(4): 1049-1058, doi: 10.3969/j.issn.0001-5733.2009.04.022.

杜增利, 徐峰, 高宏亮等. 2010. 虛譜法交錯網(wǎng)格地震波場數(shù)值模擬. 石油物探, 49(5): 430-437.

高煒, 劉財, 郭智奇等. 2014. 裂縫誘導(dǎo)雙相HTI介質(zhì)模型及其彈性波傳播方程. 世界地質(zhì), 33(4): 904-909.

郭智奇, 劉財, 楊寶俊等. 2007. 粘黏彈各向異性介質(zhì)中地震波場模擬與特征. 地球物理學(xué)進(jìn)展, 22(3): 804-810.

孔麗云, 王一博, 楊慧珠. 2012. 裂縫誘導(dǎo)HTI雙孔隙介質(zhì)中的裂縫參數(shù)分析. 地球物理學(xué)報, 55(1): 189-196, doi: 10.6038/j.issn.0001-5733.2012.01.018.

李紅星, 陶春輝. 2009. 雙相各向異性隨機(jī)介質(zhì)偽譜法地震波場特征分析. 物理學(xué)報, 58(4): 2836-2842.

劉財, 郭智奇, 楊寶俊等. 2007. 黏彈各向異性介質(zhì)中波的反射與透射問題分析. 地球物理學(xué)報, 50(4): 1216-1224

劉財, 蘭慧田, 郭智奇等. 2013. 基于改進(jìn)BISQ機(jī)制的雙相HTI介質(zhì)波傳播偽譜法模擬與特征分析. 地球物理學(xué)報, 56(10): 3461-3473, doi: 10.6038/cjg20131021.

劉炯, 馬堅偉, 楊慧珠等. 2008. 縫洞型儲層錯格偽譜法模擬. 石油地球物理勘探, 43(6): 723-727.

劉洋, 李承楚. 2000. 雙相各向異性介質(zhì)中彈性波傳播偽譜法數(shù)值模擬研究. 地震學(xué)報, 22(2): 132-138.

單啟銅, 樂友喜. 2007. PML邊界條件下二維粘彈性介質(zhì)波場模擬. 石油物探, 46(2): 126-130.

唐杰, 方兵, 孫成禹等. 2015. 基于各向異性流體替換的裂隙介質(zhì)波傳播特征研究. 石油物探, 54(1): 1-8.

吳寶年, 吳肅琴, 黃義等. 2012. 交錯網(wǎng)格在偽譜法彈性波場數(shù)值模擬中的應(yīng)用. 石油物探, 51(5): 430440-437445.

軒義華, 何樵登, 孟慶生等. 2006. 基于BISQ機(jī)制的雙相EDA介質(zhì)的波場分析. 石油地球物理勘探, 41(5): 550-556.

張軍舵, 樂友喜, 王艷香. 2008. 雙相各向同性介質(zhì)偽譜法地震波場數(shù)值模擬. 石油物探, 47(4): 338-345.

張顯文, 王德利, 王者江等. 2010. 基于BISQ機(jī)制三維雙相正交裂隙各向異性介質(zhì)衰減及頻散方位特性研究. 地球物理學(xué)報, 53(10): 2452-2459, doi: 10.3969/j.issn.0001-5733.2010.10.019.

(本文編輯汪海英)

Seismic wavefield simulation and feature analysis for a fracture-induced two-phase HTI medium based on the staggered-grid pseudo-spectral method

LIU Cai1,2, CHI Huan-Zhao2, GAO Wei1, LU Qi1*, LAN Hui-Tian3

1CollegeofGeo-explorationScienceandTechnology,JilinUniversity,Changchun130026,China2CollegeofEarthSciences,JilinUniversity,Changchun130061,China3ExplorationandDevelopmentResearchInstituteofDaqingOilfieldCompanyLtd.,HeilongjiangDaqing163712,China

AbstractA model for a fracture-induced horizontal transverse isotropic (HTI) double-porosity medium has been established by embedding a set of vertical fractures into the isotropic porous background rock in previous work, which considers the fractured reservoir rock anisotropy and porosity. High-precision seismic wavefield simulation is the main technique for the investigation of waves propagating in this media. In this paper, seismic wavefield is simulated using a single-layer model and two-layer model by solving the one-order velocity-stress equations describing seismic wave propagation in this medium by combining staggered-grid pseudo-spectral method and time-splitting method. Then the wavefield features are analyzed. The results of numerical simulation indicate that the artifacts in the wavefield simulation results by the standard-grid pseudo-spectral method can be removed by the staggered-grid pseudo-spectral method, and the stable and high-precision simulations results can be obtained by combining it with the time-splitting method. The seismic wavefield in the fractured-induced two-phase HTI medium has the characteristics of wavefield of seismic waves propagating in a fracture anisotropic medium and two-phase medium.

KeywordsFracture-induced; Two-phase HTI media; Staggered-grid pseudo-spectral method; Time-splitting method; Wavefield characteristics

基金項目國家重點(diǎn)基礎(chǔ)研究發(fā)展計劃(2013CB429805)和國家自然科學(xué)基金項目(41430322，41340039)共同資助.

作者簡介劉財,男,1963年生,教授,博士生導(dǎo)師,主要從事地震波場正反演理論、綜合地球物理等研究. E-mail:liucai@jlu.edu.cn *通訊作者鹿琪,女,1972年生,博士,教授,主要從事探地雷達(dá)和地震數(shù)據(jù)處理研究. E-mail:luqi@jlu.edu.cn

doi:10.6038/cjg20160521 中圖分類號P631

收稿日期2015-05-14，2016-03-02收修定稿

劉財， 遲喚昭， 高煒等. 裂縫誘導(dǎo)雙相HTI介質(zhì)地震波場錯格偽譜法模擬與波場特征分析. 地球物理學(xué)報，59(5):1776-1789,doi:10.6038/cjg20160521.

Liu C, Chi H Z, Gao W, et al. 2016. Seismic wavefield simulation and feature analysis for a fracture-induced two-phase HTI medium based on the staggered-grid pseudo-spectral method.ChineseJ.Geophys. (in Chinese),59(5):1776-1789,doi:10.6038/cjg20160521.