李曉波,董良國(guó)

(同濟(jì)大學(xué)海洋地質(zhì)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海200092)

斑狀飽和是沉積巖中重要的流體賦存形式,在近地表的巖石中以氣、水斑狀分布為主,在深部?jī)?chǔ)集層中還存在油或其它烴類的分布形式[1]。Cadoret等通過CT技術(shù)呈現(xiàn)了石灰?guī)r中厘米尺度的氣、水斑塊狀分布特征[2];Quintal等也通過CT技術(shù)證明了砂巖中存在氣、油斑狀分布的客觀事實(shí)[3]。地震波在斑狀飽和介質(zhì)中傳播時(shí)會(huì)發(fā)生較強(qiáng)能量衰減。White(1975)首次考慮介觀尺度(mesoscale,遠(yuǎn)大于孔隙尺寸且遠(yuǎn)小于地震波波長(zhǎng))斑狀飽和介質(zhì),成功預(yù)測(cè)了地震勘探頻帶(1～1000Hz)內(nèi)的地震波強(qiáng)衰減和速度頻散特征[4]。理論分析[5-8]和實(shí)驗(yàn)室測(cè)量[9]均表明,介觀尺度波致流體流動(dòng)(wave-induced fluid flow)是導(dǎo)致斑狀飽和介質(zhì)中地震波本征衰減的主要原因。斑狀飽和介質(zhì)中的地震波傳播理論考慮了地震波衰減的物理機(jī)制,且衰減與儲(chǔ)層參數(shù)之間關(guān)系密切[10-11],因此研究斑狀飽和介質(zhì)中的地震波傳播特征對(duì)于儲(chǔ)層預(yù)測(cè)非常重要。利用孔隙彈性理論[12-14]可以研究斑狀飽和介質(zhì)中的地震波傳播特征[15-16],基于該理論的數(shù)值模擬技術(shù)常被用于模擬斑狀飽和介質(zhì)中的地震波傳播[17-18],并準(zhǔn)確求出地震勘探頻帶內(nèi)的地震波強(qiáng)衰減和速度頻散特征[1,19-21]。然而,描述斑狀飽和介質(zhì)所需的儲(chǔ)層參數(shù)比較多,并且需要采用厘米甚至毫米級(jí)的空間網(wǎng)格才能描述介觀尺度的非均質(zhì)特征[19,22],導(dǎo)致地震波數(shù)值模擬中的計(jì)算量超過常規(guī)彈性波正演的千萬倍,使得基于斑狀飽和介質(zhì)的地震波傳播數(shù)值模擬技術(shù)無法真正應(yīng)用于實(shí)際。為此,一些學(xué)者考慮用等效粘彈模型定量表征斑狀飽和介質(zhì)中的地震波衰減和速度頻散特征[10,23-25]。

斑狀飽和介質(zhì)中的地震波衰減和速度頻散可由低頻極限(弛豫狀態(tài))模量M0以及高頻極限(非弛豫狀態(tài))模量M∞和特征頻率fc(衰減峰值發(fā)生的頻率)3個(gè)重要參數(shù)描述[26],標(biāo)準(zhǔn)線性體(SLS)[26-27]則是結(jié)合這3個(gè)參數(shù)的一個(gè)理想粘彈模型。該模型雖然未特別強(qiáng)調(diào)地震波衰減發(fā)生的物理機(jī)制[26],但可以預(yù)測(cè)出斑狀飽和介質(zhì)中的地震波衰減和速度頻散等粘彈性質(zhì),而且使描述孔隙介質(zhì)的參數(shù)大幅度減少,便于我們通過數(shù)值模擬來研究斑狀飽和介質(zhì)中的地震波傳播規(guī)律,極大地降低了數(shù)值計(jì)算量,因而具有更強(qiáng)的實(shí)用性。Sidler等[28]首次對(duì)比了基于斑狀飽和與粘滯彈性理論模擬得到的地震波記錄。我們首先將復(fù)雜的儲(chǔ)層參數(shù)簡(jiǎn)化為高、低頻極限彈性模量及特征頻率3個(gè)變量,采用SLS模型定量表征斑狀飽和介質(zhì)中的地震波衰減及速度頻散特征,進(jìn)而通過數(shù)值模擬分析孔隙度和含氣飽和度的變化對(duì)斑狀飽和介質(zhì)中地震波振幅和相位特征的影響。

1 斑狀飽和與粘滯彈性理論簡(jiǎn)介

1.1 斑狀飽和理論

一般情況下,孔隙儲(chǔ)層的儲(chǔ)集空間中總有多種流體(如氣、水、油等)同時(shí)充填。斑狀飽和介質(zhì)考慮兩種流體同時(shí)充填,采用巖石彈性模量、孔隙度、滲透率、流體粘滯系數(shù)、流體飽和度、密度等儲(chǔ)層參數(shù)來定量刻畫儲(chǔ)層單元的特征?？紤]一斑狀飽和的孔隙儲(chǔ)層巖石單元,其骨架的體積模量為Kdry,剪切模量為μdry,孔隙度為φ。組成巖石骨架的固體顆粒的體積模量為Ks,孔隙中充填氣、水兩種流體,其體積模量分別為Kg和Kw,含氣飽和度為Sg。在斑狀飽和介質(zhì)中,流體分布是非均質(zhì)的,存在一個(gè)臨界弛豫尺度Lc,與孔隙流體的體積模量Kf和粘滯系數(shù)ηf以及巖石骨架的孔隙度φ和滲透率κ共同決定了地震波衰減發(fā)生的特征頻率fc[29]:

(1)

當(dāng)頻率足夠低時(shí),弛豫尺度無限大,由波動(dòng)誘發(fā)的孔隙壓力在整個(gè)孔隙空間中可以有足夠的時(shí)間達(dá)到平衡,此時(shí)巖石的體積模量可由Gassmann理論[30]準(zhǔn)確計(jì)算得到,即

(2)

式中:α=1-Kdry/Ks;M=Ks/(α-φ+φKs/Kf)。此時(shí)Kf為極限低頻條件下的流體體積模量,可由有效流體模型預(yù)測(cè)得到[26]

(3)

當(dāng)頻率無限高時(shí),弛豫尺度無限小,由波動(dòng)引起的孔隙壓力在含不同流體的區(qū)域無法達(dá)到平衡,含水和含氣的巖石具有不同的體積模量K(Kw)和K(Kg)(可用Gassmann理論作局部性描述),但是整個(gè)巖石的剪切模量不受流體影響,即μ=μdry。因此基于Hill平均[31]可準(zhǔn)確求取高頻條件下斑狀飽和巖石的有效體積模量:

(4)

根據(jù)(5)式和(6)式即可求得斑狀飽和模型中P波模量的高、低頻極限:

1.2 粘彈介質(zhì)模型

粘彈介質(zhì)模型常被用于描述復(fù)雜介質(zhì)中的地震波傳播。SLS模型(圖1)是一個(gè)較理想的線性粘彈性體,由一個(gè)彈性單元與一個(gè)粘性單元并聯(lián)后,再與一個(gè)彈性單元串聯(lián)組成,其應(yīng)力-應(yīng)變關(guān)系為[26]

(7)

圖1 三元件標(biāo)準(zhǔn)線性粘彈性體模型[26]

當(dāng)頻率趨近于0和趨近于無窮大時(shí),根據(jù)(8)式可以得到SLS模型復(fù)模量的極限：

M∞=E2ω→∞

(10)

M0和M∞都是與頻率無關(guān)的實(shí)數(shù),因此，在無限高頻和無限低頻條件下SLS模型呈現(xiàn)為彈性性質(zhì)。復(fù)模量也可以簡(jiǎn)單地用兩個(gè)極限和特征頻率來表示[26]:

(11)

其中，ωc為特征圓頻率，即

(12)

1.3 粘彈性表征

逆品質(zhì)因子Q-1常被用于定量表征介質(zhì)的粘彈性質(zhì),具有多種表達(dá)形式,其中復(fù)模量的虛部與實(shí)部之比是一種精確的表達(dá)方式[26]:

(13)

速度頻散常伴隨地震波衰減發(fā)生,根據(jù)復(fù)彈性模量和介質(zhì)密度ρ即可得到不同頻率條件下的P波傳播速度:

(14)

為了分析SLS模型對(duì)斑狀飽和儲(chǔ)層單元粘滯性的表征能力,考慮一斑狀飽和巖石,其儲(chǔ)層參數(shù)如表1所示(根據(jù)實(shí)驗(yàn)室測(cè)量結(jié)果[32]計(jì)算所得),孔隙度為20%,含氣飽和度為10%。其中每個(gè)含氣斑塊的儲(chǔ)層巖石單元為一邊長(zhǎng)1m的立方體,此時(shí)含氣斑塊的直徑約為0.58m(與介觀尺度的定義一致),特征頻率約為32Hz。

表1 斑狀飽和儲(chǔ)層單元的物性參數(shù)

基于斑狀飽和介質(zhì)模型和SLS模型,從理論上預(yù)測(cè)了表1中孔隙介質(zhì)的地震波衰減和速度頻散特征,見圖2。圖2中藍(lán)色曲線表示SLS模型的預(yù)測(cè)結(jié)果,紅色曲線表示斑狀飽和理論[33]的預(yù)測(cè)結(jié)果?？梢钥闯?地震波衰減在低頻段隨著頻率增加而增加,在特征頻率達(dá)到最大;在高頻段又隨著頻率的增加逐漸降低。P波速度的低頻極限約為3250m/s。隨著頻率的增加，P波速度逐漸增大,在特征頻率附近P波速度增加最快,在高頻段又逐漸趨近于3470m/s?；赟LS模型預(yù)測(cè)的地震波最大衰減略高于斑狀飽和理論的預(yù)測(cè)結(jié)果(圖2a),P波速度更快地達(dá)到高頻極限(圖2b),但二者的總體特征表現(xiàn)出了很強(qiáng)的一致性,說明至少?gòu)目v波的衰減和頻散角度看,利用SLS粘彈模型可以較好地表征比較復(fù)雜的斑狀飽和介質(zhì)模型。

圖2 斑狀飽和介質(zhì)(表1所示參數(shù))中縱波衰減(a)和速度頻散(b)理論曲線

2 基于SLS模型的斑狀飽和介質(zhì)中地震波傳播模擬

從上述分析可知,SLS模型可以有效地表征斑狀飽和介質(zhì)的粘彈(地震波衰減和速度頻散)特征。與斑狀飽和介質(zhì)模型相比,SLS模型描述介質(zhì)的粘彈性質(zhì)所需要的參數(shù)大幅度減少?；赟LS模型的波動(dòng)方程易于數(shù)值求解,極大地降低了地震波傳播數(shù)值模擬的計(jì)算量,更便于研究孔隙介質(zhì)的粘彈特性。為此,本文在建立斑狀飽和介質(zhì)模型與SLS模型在表征孔隙介質(zhì)的粘彈特性關(guān)系的基礎(chǔ)上,采用基于SLS模型的地震波傳播數(shù)值模擬方法,研究斑狀飽和介質(zhì)中的地震波傳播特征。

對(duì)于表1中的均勻斑狀飽和介質(zhì)模型,為了得到地震勘探頻帶內(nèi)所有頻率成份的地震波信息,數(shù)值模擬試驗(yàn)中選取寬頻Ricker子波[34]作為子波函

數(shù),即

(15)

兩個(gè)主頻分別為f0=50Hz，f1=500Hz,子波延遲時(shí)t0=20ms。時(shí)間、空間采樣間隔分別采用0.05ms和1m;震源點(diǎn)(圖3中紅色五角星標(biāo)注)坐標(biāo)為(0,0),模型深度為200m。

采用VSP觀測(cè)系統(tǒng),記錄201個(gè)不同深度位置(圖3中黃色五邊形標(biāo)注)接收的地震記錄。由圖3可見,隨著傳播距離的增加,地震波振幅快速減小,并且發(fā)生了嚴(yán)重的物理頻散現(xiàn)象,地震波波形變化明顯。從對(duì)應(yīng)的振幅譜(圖4)上可以看出,隨著傳播距離的增大,地震波能量發(fā)生了較強(qiáng)的衰減,且主要能量的頻帶快速向低頻移動(dòng)。

為了進(jìn)一步驗(yàn)證通過SLS模型模擬斑狀飽和介質(zhì)中地震波傳播思路的正確性,我們利用基于SLS模型模擬的地震波記錄求取地震波衰減和速度頻散曲線,并將其與斑狀飽和介質(zhì)理論的預(yù)測(cè)結(jié)果進(jìn)行對(duì)比(圖5)。由圖5可以看出,除少數(shù)頻率成分的計(jì)算結(jié)果存在較小的偏差外,根據(jù)SLS粘彈模型模擬的地震波記錄可以較好地預(yù)測(cè)出斑狀飽和介質(zhì)中整個(gè)地震勘探頻帶的P波逆品質(zhì)因子和P波傳播速度隨頻率的變化特征。因此,根據(jù)SLS模型模擬的地震記錄可較好地計(jì)算斑狀飽和介質(zhì)中的地震波衰減和速度頻散特征,為研究斑狀飽和介質(zhì)中地震波場(chǎng)特征與儲(chǔ)層參數(shù)的定量關(guān)系提供了一種簡(jiǎn)單、快速的數(shù)值分析工具。

圖3 一維均勻斑狀飽和介質(zhì)中VSP觀測(cè)方式及模擬記錄

圖4 零偏VSP模擬記錄(壓力場(chǎng))的頻譜

圖5 斑狀飽和介質(zhì)(表1所示參數(shù))中縱波衰減(a)和速度頻散(b)理論曲線(藍(lán)色散點(diǎn)線為基于SLS模型模擬的地震記錄的數(shù)值計(jì)算結(jié)果,紅色曲線為斑狀飽和理論的預(yù)測(cè)結(jié)果)

3 儲(chǔ)層參數(shù)對(duì)斑狀飽和介質(zhì)中地震波場(chǎng)特征的影響分析

前人對(duì)斑狀飽和介質(zhì)中的地震波衰減和速度頻散特征研究結(jié)果表明,含氣飽和度和孔隙度是影響斑狀飽和介質(zhì)中地震波衰減和速度頻散的兩個(gè)重要儲(chǔ)層參數(shù)[10-11],也是重要的儲(chǔ)層表征參數(shù),因此,研究含氣飽和度和孔隙度對(duì)地震記錄的影響顯得尤為重要。但是,采用斑狀飽和介質(zhì)模型進(jìn)行地震波模擬所需參數(shù)多,計(jì)算量異常龐大,為此,我們?cè)诙勘碚靼郀铒柡徒橘|(zhì)模型粘彈效應(yīng)的基礎(chǔ)上,采用基于SLS模型的地震波數(shù)值模擬方法,研究斑狀飽和介質(zhì)中含氣飽和度和孔隙度對(duì)地震波傳播特征的影響。

3.1 含氣飽和度的影響

根據(jù)表1所示儲(chǔ)層參數(shù)建立均勻斑狀飽和模型,固定孔隙度(20%),改變含氣飽和度。由于單元體尺度固定,所以當(dāng)含氣飽和度增大時(shí),含氣斑塊的尺度也相應(yīng)增大。采用SLS模型定量表征該組斑狀飽和模型,圖6顯示了不同含氣飽和度條件下(不同顏色)的地震波衰減和速度頻散特征。對(duì)比可知,含氣飽和度為20%時(shí)(圖6中綠色曲線),地震波衰減的峰值最大,相應(yīng)的速度頻散效應(yīng)也最強(qiáng)；均勻含水(圖6中藍(lán)色曲線)和均勻含氣(圖6中黃色曲線)情況下,幾乎不發(fā)生衰減和速度頻散。

分別對(duì)6種含氣飽和度條件下的斑狀飽和介質(zhì)進(jìn)行地震波傳播模擬實(shí)驗(yàn),子波函數(shù)為主頻50Hz的Ricker子波,時(shí)間、空間采樣間隔分別為0.05ms和1m,記錄距震源200m處的質(zhì)點(diǎn)震動(dòng)速度(如圖7所示)。由圖7可見，當(dāng)介質(zhì)均勻含氣或含水時(shí),孔隙介質(zhì)中的流體性質(zhì)完全相同,地震波傳播時(shí)無法產(chǎn)生較大的孔隙壓力梯度,因而流體流動(dòng)及引起的地震波衰減較弱,地震記錄的振幅和相位幾乎沒有發(fā)生變化(圖7中含氣飽和度為0和1時(shí))。當(dāng)介質(zhì)中同時(shí)充填氣、水兩種流體時(shí),由于氣、水的可壓縮性相差較大,地震波傳播過程中在氣、水界面會(huì)產(chǎn)生較大的壓力梯度,從而導(dǎo)致較強(qiáng)的流體流動(dòng)和地震波衰減,地震記錄的波形也會(huì)發(fā)生相應(yīng)的變化。分析圖6可知,含氣飽和度為20%時(shí)的地震波衰減和速度頻散最明顯,從而導(dǎo)致該條件下地震記錄的振幅和相位變化也最劇烈(圖7);含氣飽和度為40%時(shí),地震記錄也發(fā)生了較明顯的變化,其它條件下地震波波形變化較小。

圖6 孔隙度為20%時(shí)不同含氣飽和度條件下地震波衰減(a)和速度頻散(b)特征

圖7 不同含氣飽和度條件下直達(dá)波波形對(duì)比(孔隙度為20%)

3.2 孔隙度的影響

孔隙度變化可以顯著改變孔隙介質(zhì)儲(chǔ)層巖石的彈性模量和滲透率,并且影響孔隙介質(zhì)的本征衰減特征。為了研究孔隙度變化對(duì)地震波傳播的影響,改變孔隙度大小,不同孔隙度對(duì)應(yīng)的巖石物性參數(shù)如表2所示(根據(jù)實(shí)驗(yàn)室測(cè)量結(jié)果[34]計(jì)算所得),儲(chǔ)層中充填氣、水兩種流體,固定含氣飽和度(10%),建立不同的斑狀飽和儲(chǔ)層模型。圖8顯示了相同含氣飽和度(10%)、不同孔隙度條件下的地震波衰減和速度頻散特征。由圖8可見,低孔隙度條件下速度頻散效應(yīng)與地震波衰減顯著減小,孔隙度為5%時(shí)(圖8中藍(lán)色曲線)的衰減峰值最小。孔隙度的增大,顯著增強(qiáng)了介質(zhì)對(duì)地震波的耗散能力,地震波衰減和速度頻散效應(yīng)也隨著孔隙度的增大逐漸增強(qiáng)。

分別對(duì)6種孔隙度條件下的斑狀飽和介質(zhì)進(jìn)行數(shù)值模擬實(shí)驗(yàn),子波函數(shù)為主頻50Hz的Ricker子波,延遲時(shí)為20ms,時(shí)間采樣間隔為0.05ms,空間采樣間隔為1m。圖9顯示了不同孔隙度條件下質(zhì)點(diǎn)震動(dòng)速度的波形記錄,每道中的兩個(gè)震相分別對(duì)應(yīng)傳播距離為0和200m時(shí)的直達(dá)波。對(duì)比分析可知,孔隙度為30%對(duì)應(yīng)的地震波振幅衰減量明顯高于低孔隙度條件下的地震記錄,且相位也發(fā)生了明顯的變化,而低孔隙度條件下的地震波能量和相位變化均較小?？紫抖仍龃笫沟酶蓭r石的體積模量和剪切模量顯著降低,巖石骨架的可壓縮性增強(qiáng),相同的地震波傳播經(jīng)過高孔隙度巖石時(shí)可以產(chǎn)生較大的孔隙壓力梯度,引起較強(qiáng)的流體流動(dòng)和地震波衰減,從而導(dǎo)致地震記錄的波形發(fā)生較大的變化。此外,高孔隙度條件下地震波傳播速度較低,傳播200m后到達(dá)的波形有較大的延遲量。

表2 不同孔隙度條件下干巖石的物性參數(shù)

圖8 含氣飽和度為10%時(shí)不同孔隙度條件下地震波衰減(a)和速度頻散(b)特征

圖9 不同孔隙度、相同含氣飽和度(10%)條件下質(zhì)點(diǎn)振動(dòng)速度對(duì)比(每道中的兩個(gè)震相分別對(duì)應(yīng)傳播距離為0和200m時(shí)的直達(dá)波)

4 結(jié)束語

本文將斑狀飽和介質(zhì)的儲(chǔ)層參數(shù)簡(jiǎn)化為高頻極限模量、低頻極限模量和特征頻率3個(gè)參數(shù),采用SLS模型定量表征了介觀尺度非均質(zhì)儲(chǔ)層單元的衰減和速度頻散特征,同時(shí),利用地震波傳播數(shù)值模擬技術(shù),分析了斑狀飽和對(duì)地震波傳播的影響。數(shù)值模擬研究結(jié)果表明:

1) SLS模型可有效地定量表征斑狀飽和介質(zhì)中地震波衰減和速度頻散特征。

2) 利用較簡(jiǎn)單的SLS粘彈模型,可以方便快速地模擬斑狀飽和介質(zhì)中的地震波傳播。

3) 斑狀飽和介質(zhì)中的強(qiáng)衰減和頻散特征顯著影響地震波記錄的振幅和相位,根據(jù)SLS模型模擬的地震波記錄可以比較準(zhǔn)確地提取出斑狀飽和介質(zhì)中的地震波衰減和頻散特征。

4) 在含氣飽和度20%附近和高孔隙度的斑狀飽和介質(zhì)中,地震波的振幅衰減和相位變化比較大；含氣飽和度和孔隙度漸變時(shí),地震波的振幅和相位呈現(xiàn)相應(yīng)的漸變趨勢(shì)。這些波場(chǎng)變化特征可以為儲(chǔ)層參數(shù)的預(yù)測(cè)研究提供分析依據(jù)。

參 考 文 獻(xiàn)

[1] Masson Y J,Pride S R.Seismic attenuation due to patchy saturation[J].Journal of Geophysical Research,2011,116:B03206

[2] Cadoret T,Mavko G,Zinszner B.Fluid distribution effect on sonic attenuation in partially saturated limestones[J].Geophysics,1998,63(1):154-160

[3] Quintal B,Frehner M,Madonna C,et al.Integrated numerical and laboratory rock physics applied to seismic characterization of reservoir rocks[J].The Leading Edge,2011,30(12):1360-1367

[4] White J E.Computed seismic speeds and attenuation in rocks with partial gas saturation[J].Geophysics,1975,40(2):224-232

[5] Pride S R,Berryman J G,Harris J M.Seismic attenuation due to wave-induced flow[J].Journal of Geophysical Research,2004,109:B01201

[6] Toms J,Müller T M,Ciz R,et al.Comparative review of theoretical models for elastic wave attenuation and dispersion in partially saturated rocks[J].Soil Dynamics and Earthquake Engineering,2006,26(6-7):548-565

[7] Müller T M,Gurevich B,Lebedev M.Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks:a review[J].Geophysics,2010,75(5):A147-A164

[8] Johnson D L.Theory of frequency dependent acoustics in patchy-saturated porous media[J].The Journal of the Acoustical Society of America,2001,110(2):682

[9] Tisato N,Quintal B.Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone:evidence of fluid flow on the mesoscopic scale[J].Geophysical Journal International,2013,195:342-351

[10] Carcione J M,Picotti S.P-wave seismic attenuation by slow-wave diffusion:effects of inhomogeneous rock properties[J].Geophysics,2006,71(3):O1-O8

[11] Ciz R,Gurevich B,Markov M.Seismic attenuation due to wave-induced fluid flow in a porous rock with spherical heterogeneities[J].Geophysical Journal International,2006,165(3):957-968

[12] Biot M A.Theory of propagation of elastic waves in a fluid-saturated porous solid I:low frequency range[J].The Journal of the Acoustical Society of America,1956,28(2):168-178

[13] Biot M A.Theory of propagation of elastic waves in a fluid-saturated porous solid II:higher frequency range[J].The Journal of the Acoustical Society of America,1956,28(2):179-191

[14] Biot M A.Mechanics of deformation and acoustic propagation in porous media[J].Journal of Applied Physics,1962,33(4):1482-1498

[15] Dutta N C,Odé H.Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation (white model)—part Ⅰ:Biot theory[J].Geophysics,1979,44(11):1777-1788

[16] Dutta N C,Odé H.Attenuation and dispersion of compressional waves in fluid-filled porous rocks with partial gas saturation (white model)—part Ⅱ:results[J].Geophysics,1979,44(11):1789-1805

[17] Carcione J M,Helle H B,Pham N H.White’s model for wave propagation in partially saturated rocks:comparison with poroelastic numerical experiments[J].Geophysics,2003,68(4):1389-1398

[18] Rubino J G,Santos J E,Picotti S,et al.Simulation of upscaling effects due to wave-induced fluid flow in Biot media using the finite-element method[J].Journal of Applied Geophysics,2007,62(3):193-203

[19] Masson Y J,Pride S R.Poroelastic finite difference modeling of seismic attenuation and dispersion due to mesoscopic-scale heterogeneity[J].Journal of Geophysical Research,2007,112:B03204

[20] Wenzlau F,Müller T M.Finite-difference modeling of wave propagation and diffusion in poroelastic media[J].Geophysics,2009,74(4):T55-T66

[21] Rubino J G,Ravazzoli C L,Santos J E.Equivalent viscoelastic solids for heterogeneous fluid-saturated porous rocks[J].Geophysics,2009,74(1):N1-N13

[22] 王東,張海瀾,王秀明.部分飽和孔隙巖石中聲波傳播數(shù)值研究[J].地球物理學(xué)報(bào),2006,49(2):524-532

Wang D,Zhang H L,Wang X M.A numerical study

of acoustic wave propagation in partially saturated poroelastic rock[J].Chinese Journal of Geophysics,2006,49(2):524-532

[23] Chapman M,Liu E,Li X.The influence of fluid-sensitive dispersion and attenuation on AVO analysis[J].Geophysical Journal International,2006,167(1):89-105

[24] Rubino J G,Ravazzoli C L,Santos J E.Equivalent viscoelastic solids for heterogeneous fluid-saturated porous rocks[J].Geophysics,2009,74(1):N1-N13

[25] Picotti S,Carcione J M,Rubino J,et al.A viscoelastic representation of wave attenuation in porous media[J].Computers & Geosciences,2010,36(1):44-53

[26] Mavko G,Mukerji T,Dvorkin J.The rock physics handbook:tools for seismic analysis of porous media[M].2nd ed.Cambridge:Cambridge University Press,2009:121-127

[27] Carcione J M.Wave fields in real media:wave propagation in anisotropic,anelastic porous & electromagnetic media (Handbook of geophysical exploration,seismic exploration)[M].Amsterdam,Netherlands:Elsevier,2007:74-75

[28] Sidler R,Rubino J G,Holliger K.Quantitative comparison between simulations of seismic wave propagation in heterogeneous poro-elastic media and equivalent visco-elastic solids for marine-type environments[J].Geophysical Journal International,2013,193(1):463-474

[29] Batzle M L,Han D,Hofmann R.Fluid mobility and frequency-dependent seismic velocity—direct measurements[J].Geophysics,2006,71(1):N1-N9

[30] Gassmann F.Elastic waves through a packing of spheres[J].Geophysics,1951,16(4):673-685

[31] Hill R.Elastic properties of reinforced solids:some theoretical principles[J].Journal of the Mechanics and Physics of Solids,1963,11(5):357-372

[32] Gomez C T,Dvorkin J,Vanorio T.Laboratory measurements of porosity,permeability,resistivity,and velocity on Fontainebleau sandstones[J].Geophysics,2010,75(6):E191-E204

[33] Dutta N C,Seriff A J.White’s model of attenuation in rocks with partial gas saturation[J].Geophysics,1979,44(11):1806-1812

[34] 俞壽朋.寬帶Ricker子波[J].石油地球物理勘探,1996,31(5):605-615

Yu S P.Wide-band Ricker wavelet[J].Oil Geophysical Prospecting,1996,31(5):605-615