鄧 霖，高 陽，李志敏，張 沂，陶俊文，馬海軍

隨機(jī)介質(zhì)理論天然氣水合物羽狀流正演模擬

鄧 霖，高 陽，李志敏，張 沂，陶俊文，馬海軍

（中國石油長慶油田分公司第一采氣廠，陜西 西安 710018）

羽狀流發(fā)育地區(qū)通常會發(fā)現(xiàn)富含天然氣水合物的儲層，但其非均質(zhì)性很強(qiáng)，刻畫難度大，目前的常規(guī)建模方法已經(jīng)無法完整地描述天然氣羽狀流，需要進(jìn)一步發(fā)展天然氣水合物羽狀流建模理論。介紹了隨機(jī)介質(zhì)建模的基本理論。其次，在建立天然氣水合物羽狀流模型時(shí)，引入了隨機(jī)介質(zhì)建模理論。最后，通過正演模擬結(jié)果發(fā)現(xiàn)，該方法可以靈活、準(zhǔn)確地描述天然氣水合物羽狀流的地震波形特征，從而為尋找天然氣水合物的提供了一定的參考依據(jù)。

天然氣水合物； 羽狀流； 正演； 隨機(jī)介質(zhì)

天然氣水合物也可稱為“可然冰”，具有使用方便，燃燒值高，清潔無污染等特點(diǎn)，被科學(xué)家譽(yù)為21世紀(jì)能源。早在200多年前人類就發(fā)現(xiàn)了天然氣水合物，2007年，我國科學(xué)家證明中國海域存在豐富的天然氣水合物資源。經(jīng)研究發(fā)現(xiàn)，天然氣在地殼動力作用下，會滲入到海水中，形成羽狀流，而羽狀流發(fā)育地區(qū)通常會發(fā)現(xiàn)富含天然氣水合物的儲層[1]。由于羽狀流的非均質(zhì)性很強(qiáng)烈，主要表現(xiàn)氣泡半徑和含量與所處的海水深度有關(guān)，而常規(guī)的數(shù)值模擬方法是以均勻介質(zhì)或?qū)訝罱橘|(zhì)理論為基礎(chǔ)，沒有考慮到羽狀流的非均質(zhì)性，無法得到與這些非均質(zhì)性相聯(lián)系的地震波場特征[2-3]。因此，應(yīng)該需要尋求一種能靈活、方便、完整地描述羽狀流的非均勻性的模型。本文主要利用隨機(jī)介質(zhì)理論來模擬天然氣水合物羽狀流，并進(jìn)行正演模擬分析，為該天然氣水合物羽狀流的識別提供了一定的參考依據(jù)。

1 基本原理

1.1 隨機(jī)介質(zhì)

非均勻介質(zhì)可由大尺度和小尺度的非均勻的兩部分所組成，大尺度部分描述介質(zhì)的平均特性，即傳統(tǒng)意義上的地質(zhì)模型；小尺度部分是加在地質(zhì)模型上的隨機(jī)擾動 ，如圖1所示。 是均值為零且具有一定自相關(guān)函數(shù)及方差的空間平穩(wěn)隨機(jī)過程，通常采用各種自相關(guān)函數(shù)來求取，其中比較典型的有指數(shù)型、高斯型、Von Karman型[4]。根據(jù)Brich原理，假設(shè)介質(zhì)的縱波、橫波速度的相對擾動量是相同的，而密度的相對擾動量與其有線性關(guān)系，因此，可用縱波速度的相對擾動量來描述隨機(jī)介質(zhì)在小尺度上的非均勻性，即：

由（1）式得：

1.2 天然氣水合物羽狀流

羽狀流現(xiàn)象是由于海水中存在滲漏甲烷氣體產(chǎn)生的氣泡導(dǎo)致的，這將改變海水內(nèi)壓力分布，從而使海水的聲學(xué)特性發(fā)生改變。溶入氣泡的海水聲速與氣泡半徑和氣泡含量有關(guān)，含氣泡介質(zhì)內(nèi)聲速表達(dá)式[1]：

式中：c——?dú)庖夯旌象w的聲速，m/s；

——液體體積模量，N/m2；

K——?dú)怏w體積模量，N/m2；

——液體密度，kg/m3；

——?dú)怏w密度，kg/m3；

——頻率，Hz；

——?dú)馀莅霃剑琺；

——液體表面張力，N/m2；

——假定含氣泡兩相混合區(qū)為球形時(shí)的半徑，m；

——?dú)馀莺炕驓馀蒹w積分?jǐn)?shù)，%，即半徑為的球形區(qū)域內(nèi)

氣泡所占據(jù)的體積分?jǐn)?shù)。

其中，、K、、、、為固定參量，、、和為給定的可變參量。

由于羽狀流分布在海水中，氣泡半徑和氣泡含量與所處海水深度有關(guān)，呈現(xiàn)出很強(qiáng)的非均質(zhì)性，所以可以根隨機(jī)介質(zhì)理論求出氣泡半徑和氣泡含量的隨機(jī)擾動，然后結(jié)合（2）式建立羽狀流速度模型，如圖2所示。從圖中可以看出，羽狀流速度是隨機(jī)變化的，因?yàn)闅馀莅霃胶秃渴请S機(jī)變化的，所以得到的速度也是隨機(jī)變化的。

圖2 羽狀流速度模型

2 模型正演

2.1 天然氣水合物羽狀流建模

羽狀流處于海水中，把圖2中的羽狀流速度模型放入海水速度模型中得到羽狀流速度模型，模型的規(guī)模為：縱向深1 350 m，橫向?qū)? 000 m，如圖3所示。圖中黑色部分為速度均一的海水，速度為1 500 m/s，中間部分為羽狀流，羽狀流位于海底深350~1 350 m處，橫向位于模型中間，即900~1 100 m。與實(shí)際探測到的羽狀流對比，圖3模型基本能夠較好地反映出實(shí)際羽狀流的特征。

2.2 地震正演

對上述天然氣水合物羽狀流速度模型，進(jìn)行模型正演，得到正演模擬結(jié)果，如圖4所示。從圖中可以看出，羽狀流模型部分反射波雜亂，剖面的中部能量強(qiáng)，而上下兩端較弱，當(dāng)遇到該種反射特征，天然氣水合物羽狀流的可能性比較大。該成果為進(jìn)一步研究羽狀流產(chǎn)生的地震響應(yīng)特征奠定了基礎(chǔ)，也為需找天然氣水合物提供了一條新的思路。

圖3 天然氣水合物羽狀流速度模型

圖4 羽狀流模型正演剖面

3 結(jié)束語

本文在前人的研究基礎(chǔ)上，針對目前天然氣水合物羽狀流的非均質(zhì)性的特點(diǎn)，基于隨機(jī)介質(zhì)理論，通過正演模擬技術(shù)，對天然氣水合物羽狀流進(jìn)行了應(yīng)用研究，得出了以下結(jié)論：

（1）從上述模擬結(jié)果可以看出，利用隨機(jī)介質(zhì)理論模擬的天然氣水合物羽狀流，具有適應(yīng)性強(qiáng)、使用方便靈活的優(yōu)點(diǎn)，能有效地描述天然氣水合物羽狀流的非均質(zhì)性。

（2）通過正演模擬發(fā)現(xiàn)，羽狀流模型部分反射波雜亂，剖面的中部能量強(qiáng)，而上下兩端較弱，當(dāng)遇到該種反射特征，天然氣水合物羽狀流的可能性比較大。

上述結(jié)論為進(jìn)一步研究羽狀流產(chǎn)生的地震響應(yīng)特征奠定了基礎(chǔ)，可以為尋找天然氣水合物的提供了一定的參考依據(jù)。

［1］李燦蘋，劉學(xué)偉，勾麗敏，等．冷泉活動區(qū)天然氣水合物上覆水體中氣泡羽狀流的數(shù)值模擬[J]．巖石力學(xué)與工程學(xué)報(bào)，2013，43(2)：391-399．

［2］楊志力，呂福亮，吳時(shí)國，等．西沙海域天然氣水合物的地庚響應(yīng)特征及分布[J]．地球物理學(xué)進(jìn)展，2013，28(2)：3307-3312．

［3］姚姚，奚先．區(qū)域多尺度隨機(jī)介質(zhì)模型及其波場分析[J]．石油物探，2004，43(1)：1-7．

［4］馬靈偉，顧漢明，趙迎月，等．應(yīng)用隨機(jī)介質(zhì)正演模擬刻畫深水區(qū)臺緣礁碳酸鹽巖儲層[J]．石油地球物理勘探，2013，48(4)：583-590+676+502．

Forwarding Modeling of Natural Gas Hydrate Plumes Based on Random Medium Theory

(PetroChina Changqing Oilfield Company No.1 Gas Production Plant, Shaanxi Xi'an 710018, China)

There are always rich nature gas hydrate reservoirs in plume development areas, but the plume has a strong heterogeneity, so it is difficult to describe it. Current conventional model-building is difficult to completely describe the nature gas hydrate plume, so we need develop the modeling-building theory of nature gas hydrate plume. In this paper, the basic theory of random medium model-building was firstly described. And then the theory of random medium model-building wasintroducedinto model-building of nature gas hydrate plume. The result of forwarding modeling shows that this method can describe the seismic wave characteristics of nature gas hydrate plume flexibly and accurately, which can provide some reference for searching nature gas hydrate.

nature gas hydrate; plume; forward; random medium

TE19

A

1004-0935（2017）03-0251-03

2017-01-16

鄧霖（1986-），男，助理工程師，陜西西安人，研究方向：從事天然氣勘探與開發(fā)工作。