南海北部似海底反射層速度結(jié)構(gòu)全波形反演

霍元媛,張 明

(1.中石化華北分公司勘探開發(fā)研究院,鄭州 450006;2.廣州海洋地質(zhì)調(diào)查局,廣州 510760)

0 前言

天然氣水合物是一種新型能源,并作為一種潛在的石油天然氣替代品得到了廣泛的研究[1～3]。其中,似海底反射(BSR)是識(shí)別天然氣水合物最直觀的證據(jù),它是一個(gè)從較高速度降至較低速度的界面反射標(biāo)志,一般與現(xiàn)代海底近于平行,并且多與海底沉積層反射斜交。相對(duì)于海底反射,具有較強(qiáng)的反射振幅和極性反轉(zhuǎn)特征[4～8]。

為了正確識(shí)別BSR,利用BSR上方、下方的彈性參數(shù)結(jié)構(gòu),來識(shí)別并分析天然氣水合物及其下方游離氣的分布特征和形成機(jī)制,許多學(xué)者對(duì)BSR層上方、下方的地震速度結(jié)構(gòu)進(jìn)行了研究[9～13]。通常層速度是基于常規(guī)疊加速度分析,由Dix公式計(jì)算均方根速度得到的,但是利用地震疊加速度進(jìn)行水合物沉積層的層速度預(yù)測(cè),存在較大的誤差,其主要原因是在現(xiàn)有的地震勘探頻率下,地震疊加速度的變化不能由地震疊加振幅敏感地反映出來,使得地震疊加速度的選取精度難以保證,造成由此求取的層速度分辨率低,誤差較大[14]。Singh[15]在1993年首次提出了全波形反演的方法,并將這一方法應(yīng)用于溫哥華岸外的天然氣水合物反射層速度結(jié)構(gòu)的研究中,此后,全波形反演的方法得到了廣泛的應(yīng)用和發(fā)展[16～20]。在我國天然氣水合物的勘探中,全波形反演方法還未發(fā)揮過作用。因此,作者在前人研究的基礎(chǔ)上,提出了一種新的全波形反演方法,并將其應(yīng)用到對(duì)我國南海北部海域天然氣水合物BSR速度結(jié)構(gòu)的研究中,為我國天然氣水合物首鉆目標(biāo)優(yōu)選提供了技術(shù)支持。

1 全波形反演

波形反演是一個(gè)強(qiáng)非線性問題,目標(biāo)函數(shù)存在眾多極值,傳統(tǒng)的迭代方法使用的局部最優(yōu)化法,易于收斂于局部極值,往往得不到全局最優(yōu)解。全局優(yōu)化算法對(duì)整個(gè)模型空間進(jìn)行搜索,通常能夠?qū)ふ业饺肿顑?yōu)值。而采用全局搜索的方法尋找最優(yōu)解時(shí),當(dāng)模型參數(shù)很多時(shí),運(yùn)算量太大。因此作者在本文采用了一種折中的辦法,即先用全局搜索方法來解決非線性程度較高的旅行時(shí)反演,再用局部搜索方法來解決非線性程度較低的波形反演問題。

目前,比較理想的用于全局優(yōu)化問題求解的算法主要是:Monte Carlo算法、模擬退火算法以及遺傳算法。Sambridge[21]通過一個(gè)多元二次函數(shù)的最優(yōu)化問題,闡述了當(dāng)未知數(shù)個(gè)數(shù)增加時(shí),遺傳算法的效率明顯高于Monte Carlo算法。因此作者利用遺傳算法來進(jìn)行旅行時(shí)反演,獲得背景速度模型,以進(jìn)一步提高傳統(tǒng)波形反演方法的計(jì)算效率。作者在本文建立的全波形反演方法主要分二個(gè)步驟:①用全局搜索的遺傳算法構(gòu)建基于旅行時(shí)反演的低頻背景速度模型;②用步驟①得到的模型作為初始模型,利用波形資料應(yīng)用共軛梯度算法,來求取高頻速度結(jié)構(gòu)。這樣最后得到的模型中就既有低頻(長(zhǎng)波長(zhǎng))信息,又有高頻(短波長(zhǎng))信息。

波形反演的目標(biāo)函數(shù),是使實(shí)際的觀測(cè)記錄波形與計(jì)算合成記錄的波形差異最小。反演所使用的地震數(shù)據(jù),經(jīng)過嚴(yán)格的拉冬變換,轉(zhuǎn)化為截距時(shí)間～慢度域(τ-P域)數(shù)據(jù)。在τ-P域中,反射以橢圓形軌跡為特征。全波形反演所使用的二種算法如下所述。

1.1 遺傳算法

遺傳算法是一種應(yīng)用廣泛的隨機(jī)搜索方法,它是模仿生物進(jìn)化的過程,按照適者生存的法則建立的一種最優(yōu)化方法。遺傳算法與傳統(tǒng)反演方法相比,具有以下一些特點(diǎn)[22]:

(1)遺傳算法對(duì)模型參數(shù)的編碼進(jìn)行操作,而不對(duì)參數(shù)進(jìn)行直接操作。

(2)遺傳算法同時(shí)對(duì)多個(gè)模型進(jìn)行操作,而不是對(duì)單獨(dú)的一個(gè)模型進(jìn)行操作。

(3)遺傳算法使用隨機(jī)的優(yōu)化策略,而不是確定的優(yōu)化規(guī)則。

作者在本文主要通過以下幾個(gè)步驟來實(shí)現(xiàn)旅行時(shí)反演:

(1)確定旅行時(shí)反演參數(shù),并指定模型參數(shù)的搜索范圍。在旅行時(shí),反演的參數(shù)為層速度,搜索范圍為1 400 m/s到2 600 m/s,允許變化范圍為±10 m/s。

(2)隨機(jī)產(chǎn)生初始模型種群。通過隨機(jī)選擇產(chǎn)生的初始模型種群。群體的維數(shù)越大,其代表性越廣泛,最終進(jìn)化到最優(yōu)解的可能性越大。但維數(shù)大的群體,勢(shì)必會(huì)造成計(jì)算時(shí)間的增加。通過對(duì)不同模型規(guī)模的試驗(yàn)分析,當(dāng)模型個(gè)數(shù)為30時(shí),基本能得到滿意的反演結(jié)果,因此可將模型規(guī)模設(shè)置為30。

(3)確定目標(biāo)函數(shù)。將目標(biāo)函數(shù)定義為實(shí)測(cè)的旅行時(shí)數(shù)據(jù)與理論計(jì)算旅行時(shí)數(shù)據(jù)的誤差平方和:

其中 tmi和tni分別為第i個(gè)接收點(diǎn)上實(shí)測(cè)的和理論計(jì)算的旅行時(shí);n為數(shù)據(jù)個(gè)數(shù)。

(4)模型參數(shù)編碼。編碼是應(yīng)用遺傳算法時(shí)要解決的首要問題,采用格雷碼對(duì)參數(shù)進(jìn)行編碼,這樣有利于提高遺傳算法的局部搜索能力,并使得交叉、變異等遺傳操作便于實(shí)現(xiàn)。

(5)目標(biāo)函數(shù)值轉(zhuǎn)換為適應(yīng)度值其中 σ是目標(biāo)函數(shù)的標(biāo)準(zhǔn)偏差。

(6)算子描述。首先計(jì)算個(gè)體的相對(duì)適應(yīng)度pi=fi/∑fi,然后根據(jù)概率分布,從初始種群中隨機(jī)選一些染色體,構(gòu)成一個(gè)新的種群。為了避免少數(shù)個(gè)體霸占整個(gè)群體,這里采用了從子代和父代中挑選出的最好的n個(gè)可以直接進(jìn)入子代的策略,這種選擇方式收斂速度也比較快。

交叉運(yùn)算決定了遺傳算法的全局搜索能力。為了加速收斂,作者在本文中選擇了多點(diǎn)交叉,隨機(jī)選擇交叉位置,交叉算子設(shè)為0.9。

變異運(yùn)算是指將個(gè)體染色體編碼串中的某些基因,用其它等位基因來替換,從而形成一個(gè)新的個(gè)體。經(jīng)過試驗(yàn),可將變異算子設(shè)定為0.05。

(7)終止條件。設(shè)置最大遺傳代數(shù)K,若已完成K次計(jì)算,則停止。在本文中K=50。

1.2 共軛梯度算法

采用共軛梯度法[23]得到短波長(zhǎng)分量,目標(biāo)函數(shù)為

其中 dcal為由模型M計(jì)算出的波場(chǎng);dobs為觀測(cè)到的波場(chǎng);m0為初始模型;模為用數(shù)據(jù)協(xié)方差矩陣cD加權(quán)的L2模;為用模型協(xié)方差矩陣cM加權(quán)的L2模。

作者在本文測(cè)試了幾種不同的cD和cM結(jié)合方式,最后設(shè)定為1.0 I和2.0 I,該模型參數(shù)為縱波速度。此外在制作合成地震記錄時(shí),反射體處出現(xiàn)的不連續(xù),會(huì)引入強(qiáng)振幅和高頻能量,使得反演不穩(wěn)定。因此由遺傳算法得到的背景速度模型必須進(jìn)行平滑,在本文中采用三角平均濾波器對(duì)模型進(jìn)行平滑。

1.3 正演算法

進(jìn)行波形反演,首先要選擇一種合適的正演模擬算法。正演模擬有多種方法,其中射線追蹤法計(jì)算成本比較低,但精確度低并且難以處理復(fù)雜介質(zhì)中不同類型波之間的耦合[24];差分法應(yīng)用范圍比較廣泛,能夠較精確地模擬任意非均勻介質(zhì)中的地震波場(chǎng),但計(jì)算成本較高,且需要采用吸收邊界條件。作者在本文中采用的是慢度法[25]進(jìn)行正演,這種方法能夠提供全波場(chǎng)解,包括地震波在地層內(nèi)傳播時(shí)形成的一次波、多次波和轉(zhuǎn)換波。它的實(shí)現(xiàn)步驟為:

(1)通過傅立葉變換和漢克爾變換,將本構(gòu)方程和彈性動(dòng)力學(xué)方程變換到頻率～慢度域。

(2)在頻率～慢度域得到反射,透射系數(shù)矩陣遞推公式。

(3)分別對(duì)頻率和慢度進(jìn)行積分變換,得到時(shí)～空域的地震道集。

圖1 某線偏移剖面(BSR位于1.65 s處,具有強(qiáng)振幅、反極性并與周圍地層斜交的特征)Fig.1 The migrated section(BSR is shown clearly on the migrated section at 1.65s,which has characteristics of strong amplitude,reverse polarity and crossing obliquely with surrounding sediments)

2 南海北部海域天然氣水合物BSR速度結(jié)構(gòu)反演

作者在本文采用了以上方法,對(duì)南海海域的地震數(shù)據(jù)進(jìn)行了反演。預(yù)處理包括觀測(cè)系統(tǒng)定義,速度分析,動(dòng)校正與疊前時(shí)間偏移,得到了偏移剖面(見圖1)。從偏移剖面上可以看到,箭頭所指處的BSR具有與海底近似平行,并與海底沉積層反射斜交,以及強(qiáng)振幅和負(fù)極性的特征。作者抽取四個(gè)CMP道集形成一個(gè)超道集,以增強(qiáng)信噪比。并且此處的BSR與海底幾乎都是水平的,滿足了在反演方法中,假設(shè)速度僅隨深度變化的條件。

在進(jìn)行正演計(jì)算時(shí),震源子波由距震源一定距離點(diǎn)上直接觀測(cè)得到,再通過慢度法得到合成記錄?？紤]到可能有游離氣的存在,使得縱波速度顯著下降,而橫波和密度受到的影響較小,為了使反演結(jié)果更加穩(wěn)定,可以將橫波速度和密度參數(shù)固定不變,而只將縱波速度作為自由變量。在反演中,根據(jù)Hamilton[26]的經(jīng)驗(yàn)公式,計(jì)算出橫波速度和密度?？v波速度搜索范圍為1 400 m/s到2 600 m/s。

經(jīng)過大量試驗(yàn)結(jié)果對(duì)比,作者認(rèn)為將遺傳算法的參數(shù)中,交叉概率設(shè)為0.9;變異概率設(shè)為0.05比較合適。在實(shí)際應(yīng)用時(shí),可在此基礎(chǔ)上,根據(jù)具體情況進(jìn)行調(diào)整。模型群體的大小,可視問題的復(fù)雜程度,計(jì)算機(jī)容量,以及計(jì)算速度來進(jìn)行綜合考慮。太大占用內(nèi)存多且計(jì)算速度慢,太小則搜索不夠徹底,應(yīng)根據(jù)具體情況權(quán)衡選擇,在本文中將初始種群個(gè)數(shù)設(shè)為30。

圖2是遺傳算法迭代過程中的平均目標(biāo)函數(shù)的統(tǒng)計(jì)圖。由圖2中可以看到,在迭代的初始階段,群體具有較快的收斂速度,在迭代十八次后,收斂速度趨于平緩?？梢哉J(rèn)為,此時(shí)種群中的個(gè)體已經(jīng)達(dá)到真實(shí)數(shù)值附近。

圖2 目標(biāo)函數(shù)平均值統(tǒng)計(jì)圖Fig.2 Average value of objective function changes with iteration

在反演時(shí),非常容易受到噪聲的影響。當(dāng)數(shù)據(jù)含噪聲水平過高時(shí),遺傳算法會(huì)出現(xiàn)種群中差別很大的模型的目標(biāo)函數(shù)值都相等,由此造成自動(dòng)停止迭代;或者收斂速度非常緩慢,造成已達(dá)到最大迭代次數(shù)卻未得到滿意解的情況。因此,通過對(duì)不同信噪比的數(shù)據(jù)進(jìn)行試驗(yàn)發(fā)現(xiàn),信噪比大于10的數(shù)據(jù),其反演得到的結(jié)果較為可靠。

圖3是波形反演的最終結(jié)果。由圖3可知,一個(gè)速度明顯高于周圍速度的高速帶(箭頭1所指處)和一個(gè)低速帶(箭頭2所指處)。在300 m處高速向低速轉(zhuǎn)換,由2 000 m/s驟降至1 450 m/s,這個(gè)速度轉(zhuǎn)換邊界對(duì)應(yīng)著地震剖面上的BSR,其上方大約為15 m厚的高速帶,這表明此處存在天然氣水合物,并暗示水合物層的厚度可能在15 m左右。在高速帶的下方,可以觀測(cè)到一個(gè)1 450 m/s的速度極小值,遠(yuǎn)低于沉積背景速度,這有力地證明了游離氣的存在。

圖3 反演結(jié)果圖(箭頭1指示可能的水合物層,箭頭2指示可能的游離氣層)Fig.3 Inversion result(Arrow 1 implies gas hydrate sediments,while Arrow 2 implies free gas sediments)

由于縱波速度與水飽和度之間的高度非線性關(guān)系,直接根據(jù)縱波速度預(yù)測(cè)游離氣的含量是不可靠的。而常用于估算沉積層中游離氣含量的Domenico[27]的方法,對(duì)于我國南海海域的深部沉積層也不一定適用。因此,我們將縱波速度作為孔隙空間中水飽和度的函數(shù),采用Biot-Geerts ma[28,29]和Gass mann[30]方程計(jì)算其值。在計(jì)算時(shí),我們采用的模型參數(shù),包括了全部含水飽和時(shí)的縱波速度(2 000 m/s)、全部含水飽和時(shí)的橫波速度(600 m/s)、孔隙度(0.55)、水的體積模量(2.1×109Pa)、氣體的體積模量(0.096×109Pa)、固體物質(zhì)的密度(2.6 g/cm3)、水的密度(1.035 g/cm3)以及氣體的密度(0.3 g/cm3)。沉積物種類包括方解石33.06%、石英和斜長(zhǎng)石31.75%、云母類34.53%,沉積物骨架的有效體積模量為53×109Pa。其中,氣體的屬性是根據(jù)BSR處的靜巖壓力計(jì)算得到,并考慮了與理想氣體的偏差。模型數(shù)值參考了前人的研究成果[27、31～34]。我們給定一個(gè)縱波速度(1 450 m/s),計(jì)算出孔隙空間中的水飽和度大約低于99%,這說明氣體含量至少為1%孔隙體積。

至于甲烷的來源,主要有三種可能的來源:①地球深部非生物成因(無機(jī)成因)的甲烷;②海底沉積巖中有機(jī)質(zhì)生成的甲烷;③綜合以上二種成因的混合來源。我們的研究區(qū)位于南海北部陸緣,珠江口盆地珠二坳陷白云凹陷南側(cè),處在陸架到深海的陸坡位置,其南與南部神狐一統(tǒng)暗沙隆起區(qū)相接,海底地形總體呈東北高、西南低的斜坡形態(tài)。根據(jù)古氣候條件,古地理環(huán)境及地震相特征綜合分析認(rèn)為:白云凹陷在始新世～早漸新世(斷陷期),具有形成大型中深湖相沉積的條件——潮濕的氣候環(huán)境,全封閉的深洼陷,高的沉積速率[36]。在該時(shí)期形成了巨厚的文昌～恩平組烴源巖,目前已處于產(chǎn)生裂解氣階段。近年來,在白云凹陷北坡～番禺低隆起的一系列天然氣發(fā)現(xiàn)及鉆探結(jié)果表明,白云北坡氣藏區(qū)氣藏充滿度高,氣源來自深部的恩平～文昌組烴源巖。通過對(duì)區(qū)內(nèi)地質(zhì)調(diào)查站位資料的分析,在區(qū)內(nèi)淺表層沉積物中,普遍存在游離氣。甲烷碳同位素δ13C測(cè)試結(jié)果顯示,研究區(qū)海域δ13C(‰)值在-46.2‰～-74.3‰(PDB)之間,平均為-60.9‰(PDB),除二個(gè)樣品的δ13C值為-46.2‰和-51‰(PDB)外,大多數(shù)樣品的δ13C值小于-57‰(PDB),證實(shí)了這一區(qū)域淺表層沉積物天然氣為生物氣。許多站位頂空氣甲烷含量在垂向上保持相對(duì)較高的豐度,特別是白云凹陷內(nèi)的甲烷,含量接近120μl/kg和200μl/kg,暗示其深部可能有持續(xù)穩(wěn)定的游離甲烷供應(yīng),來源于深部的熱解氣。由此推測(cè),研究區(qū)淺部地層中的天然氣可能兼有生物氣和熱解氣。

3 結(jié)論

作者在本文中將基于遺傳算法的全波形反演方法,應(yīng)用于我國南海北部海域的地震資料,反演得到了高分辨率的BSR速度結(jié)構(gòu),清晰的識(shí)別出了BSR的速度反轉(zhuǎn)特征。波形反演對(duì)于天然氣水合物的識(shí)別具有非常重要的作用,如果可以在鉆井目標(biāo)區(qū)進(jìn)行有針對(duì)性的波形反演,從而得到該區(qū)域精細(xì)的速度結(jié)構(gòu),就可以有效地指導(dǎo)井位的確定。在利用遺傳算法進(jìn)行反演時(shí),正演過程非常耗時(shí),而種群優(yōu)化的過程又需要反復(fù)進(jìn)行正演計(jì)算。遺傳算法非常適合于并行計(jì)算,為了加快計(jì)算速度,我們正在進(jìn)行對(duì)遺傳并行算法的研究。

通過對(duì)反演結(jié)果的分析發(fā)現(xiàn),研究區(qū)處的BSR,是由其上的高速層和下方的低速層共同作用所造成的。這個(gè)高速層可以解釋為含天然氣水合物的沉積層,而低速層可以解釋為含至少1%游離氣的薄層。并且,通過與此處AVO屬性剖面,瞬時(shí)屬性剖面,以及波阻抗反演剖面的對(duì)比,可以得出相同的結(jié)論。

另外,通過對(duì)研究區(qū)域附近鉆探結(jié)果,及區(qū)內(nèi)地質(zhì)調(diào)查站位資料的分析,我們認(rèn)為研究區(qū)淺部地層中的天然氣,可能兼有生物氣和熱解氣。致 謝:感謝雷新華教授對(duì)我的指導(dǎo)和幫助;

感謝廣州海洋地質(zhì)調(diào)查局各位老師的大力支持。

[1] LEE J H,BAEK Y S,RYU B J,et al.A seismic survey to detect natural gas hydrate in the East Sea of Korea[J].Mar Geophys Res,2005,3,26(1):51.

[2] TREHU A M,TORRESM E,LONG P E,et al.Three-dimensional distribution of gas hydrate beneath southern Hydrate Ridge:constraints from ODP Leg 204[J].Earth Planet SciLet,2004,222:845.

[3] TATSUO SAEKI,MASAO HAYASH I,TAKAO INAMOR I,et al.Velocity structure of the Kumano basin in the Nankai trough[M].Proceedings of the Fifth International Conference on Gas Hydrates,2005.

[4] SHIPLEY T H,HOUSTON M H,BULLER R T,et al.Seismic evidence for wide-spread possible gas hydrate horizons on continental slopes and margins[J].Am Assoc Pet GeolBull,1979,63:2204.

[5] MACKAYM E,JARRARD R D,WESTBROOK G K,et al.Origin of bottom-simulating reflectors;geophysical evidence from the Cascadian accretinary prism[J].Geology,1994,22:459.

[6] HOLBROOKW S,HOSK INS H,WOOD W T,et al.Methane hydrate and free gas on the Blake Ridge from vertical seismic profiling[J].Science,1996,273:1840.

[7] RE MPEL A W,BUFFETT B A.Formation and accumulation of gas hydrate in porous media[J].Geophys Res,1997,102:10151.

[8] HYNDMAN R,SPENCE G A.Seismic study of methane hydrates marine bottom simulating reflectors[J].Geophys Res,1992,97:6638.

[9] CHENGW B,LEE C S,L IU C S.Velocity structure in marine sediments with gas hydrate reflectors in offshore S W Taiwan,from OBS data tomography[J].Ter Atm Oce,2006,17(4):739.

[10]ZILLMER M,FLUEH E R,PETERSEN J.seismic investigation of a bottom simulating reflector and quantification of gas hydrate in the Black Sea[J].Geophys J Int,2005,161(3):662.

[11]JAISWAL P,ZELT C A,PECHER IA.seismic characterization of a gas hydrate system in the Gulf ofMexico using wide-aperture data[J].Geophy J Int,2006,165(1):108.

[12]NETZEBAND GL,HUBSCHER C P,GAJEWSKID J,et al.Seismic velocities from the Yaquina forearc basinoff Peru:evidence for free gas within the gas hydrate stability zone[J].Int J Earth Sci,2005,94(3):420.

[13]ECKER C,DVORK IN J,NUR A M.Est imating the amount of hydrate and free gas from marine seismic data[J].Geophysics,2000,65(2):565.

[14]陳建文,吳志強(qiáng),龔建明.天然氣水合物的地球物理識(shí)別技術(shù)[J].海洋地質(zhì)動(dòng)態(tài),2004,20(6):1.

[15]S INGH S C,M INSHULL TA,SPENCE G D.Velocity structure of gas hydrate reflector[J].Science,1993,260:204.

[16]M INSHULL T A,SINGH S C,WESTBROOK G K.seismic velocity structure at a gas hydrate reflector,offshore western Columbia,from full waveform inversion[J].Geophys Res,1994,99:4715.

[17]PECHER I A,MINSHULL T A,SINGH S C,et al.Velocity structure of a bottom s imulating reflector offshore Peru:Results from full waveform and travel time inversion of wide-angle seismic data[J].Geophys,1997,102:15345.

[18]PECHER I A,RANERO C R,HUENE R,et al.The nature and distribution of bottom simulating reflectors at the Costa Rican convergent margin[J].Geophys J Int,1998,133:219.

[19]YUAN T,SPENCE GD,HYNDMAN R D,et al.Seismic velocity studies of a gas hydrate simulating reflector on the northern Cascadia continent margin:Amplitude modeling and full wave for minversion[J].Geophys Res.,1999,104:1179.

[20]宋海斌,OSAMU M,KURAMOTO S.天然氣水合物似海底反射層的全波形反演[J].地球物理學(xué)報(bào),2003,46(1):42.

[21]SAMBRIDGE M,DRIJKONINGEN G.Genetic algorithms in seismic waveform inversion[J].Geophys,JInt.,1992,109:323.

[22]石琳珂,孫銘心,王廣國,等.地球物理遺傳反演方法[M].北京:地震出版社,2000.

[23]KORMEND I F,D IETR ICH M.Non-linear waveform inversion of plane wave seismograms in stratified elastic media[J].Geophysics,1991,56:664.

[24]KENNETT B L N,KERRY N J.seismic waves in a stratified half-space[J].Geophys J R Astron Soc,1979,57:557.

[25]張繁昌,印興耀,趙劍.慢度法全波場(chǎng)模擬及其水平慢度積分的實(shí)現(xiàn)[J].石油地球物理勘探,2003,38(6):597.

[26]HAM ILTON E L.Vp/Vs and poisson’s ratios in marine sediments and rocks[J].Acoust Soc Am,1979,66(4):1093.

[27]DOMEN ICO S N.Elastic properties of unconsolidated porous sand reservoirs[J].Geophysics,1977,42:1339.

[28]B IOTM A.Theory of propagation of elastic waves in a fluid-saturated porous solid[J].Acoustic Soc Am,1956,28:168.

[29]GEERTS MA J.Velocity-log interpretation:The effect of rock bulk compressibility[J].Soc Pet Eng J,1961,1:235.

[30]GASS MANN F.Uber die Elastizitat poroser Medien[J].Veirteljahrsschr Natur forsche Gesellschaft Ges Zurich,1951,96:1.

[31]ANGENHEISTER G,ED.,LANNDOLT-BOERNSTE IN.Numerical data and functional relationships in science and techno ledge[M].Physical properties of rocks.Springer,Berlin,1982.

[32]ATK INS P W.Physical chemistry[M].Oxford Univ Press,Oxford,1978.

[33]WY MAN R E.Petrogeo physics,the interrelationships of petrophysics,geology and geophysics,AAPG continuing education-Stratigraphic interpretation of seismic data[M].Am Assoc Pet Geol,1982.

[34]房殿勇,王汝建,邵磊,等.南海ODP184站深海相漸新統(tǒng)硅質(zhì)成巖作用[J].海洋地質(zhì)與第四紀(jì)地質(zhì),2002,22(2):75.

[35]CASTAGNA J P,BATZLE M L,EAST WOOD R L.Relationships be tween compress ional-wave and shear-wave velocities in elastic silicate rocks[J].Geophysics,1985,50:571.

[36]張光學(xué),黃永樣,陳邦彥.海域天然氣水合物地震學(xué)[M].北京:海洋出版社,2003.