亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        大洋鉆探天然氣水合物儲層測井評價研究進(jìn)展

        2020-09-08 02:41:40鐘廣法趙巒嘯
        天然氣工業(yè) 2020年8期
        關(guān)鍵詞:航次水合物站位

        鐘廣法 張 迪 趙巒嘯

        海洋地質(zhì)國家重點(diǎn)實驗室·同濟(jì)大學(xué)

        1 研究背景

        國際科學(xué)大洋鉆探計劃始于1968年,經(jīng)歷了深海 鉆 探(Deep Sea Drilling Project, 縮 寫 為 DSDP,1968—1983 年)、大洋鉆探(Ocean Drilling Program,縮寫為ODP,1985—2003年)、綜合大洋鉆探(Integrated Ocean Drilling Program,縮寫為 IODP,2003—2013年)及國際大洋發(fā)現(xiàn)(International Ocean Discovery Program,縮寫為 IODP,2013—2023 年)等 4個階段,歷時50余年。這是地球科學(xué)領(lǐng)域歷時最長、影響最為深遠(yuǎn)的國際科學(xué)研究合作計劃。

        該計劃已在全球各大洋總共執(zhí)行了400個鉆探航次,其中涉及天然氣水合物(以下簡稱水合物)的航次為23個,已在大西洋、太平洋和印度洋累計53個站位鉆遇了水合物(圖1、2)。這些站位分布于美國東部大西洋被動大陸邊緣(布萊克海脊)[1-4];墨西哥灣北部被動大陸邊緣[5];太平洋東部活動大陸邊緣,包括智利三聯(lián)點(diǎn)[6]、秘魯海溝[7-8]、中美洲海溝(危地馬拉[9-10]、墨西哥[11]和哥斯達(dá)黎加邊緣[12-13])及俄勒岡(水合物脊)[14-15]和溫哥華岸外的Cascadia邊緣[16];太平洋西部活動大陸邊緣,包括日本海東部[17]、日本南海(Nankai)海槽[18-19]及新西蘭東部Hikurangi邊緣[20];印度東部大陸邊緣(孟加拉灣)[21]及地中海[22]等地(圖1、2)。大洋鉆探水合物鉆探航次的主要科學(xué)目標(biāo)包括:獲取水合物實物樣品;了解水合物及其宿主沉積物的原位性質(zhì);估算水合物飽和度;研究擬海底反射(BSR)的成因及水合物和游離氣的分布;研究水合物的形成條件和形成機(jī)制等。

        大洋鉆探采取一站多孔和全取心策略。盡管如此,由于水合物在地表條件下不穩(wěn)定、易分解,取心難度大,而且在取心過程中水合物的分解會導(dǎo)致巖心發(fā)生擾動,呈奶油狀/湯狀,破壞了沉積物的原生沉積組構(gòu)。因此,單憑巖心分析很難全面了解水合物與宿主沉積物之間的共生關(guān)系。大洋鉆探測井技術(shù)有效地彌補(bǔ)了巖心分析之不足,它們能提供井筒附近含水合物地層連續(xù)的高分辨率地球物理觀測記錄,為理解海洋水合物及其宿主沉積物的原位特性、估算水合物飽和度、預(yù)測水合物分布提供了重要的依據(jù)。為此,筆者對大洋鉆探水合物相關(guān)航次報告和文獻(xiàn)資料進(jìn)行了較系統(tǒng)的調(diào)研與分析,在此基礎(chǔ)上,對測井資料在大洋鉆探水合物儲層評價研究方面所取得的進(jìn)展進(jìn)行綜述,并對存在的問題及未來可能的發(fā)展方向進(jìn)行了探討。

        圖1 國際科學(xué)大洋鉆探計劃水合物鉆探歷程示意圖

        圖2 國際科學(xué)大洋鉆探計劃水合物站位分布圖

        我國海洋水合物勘查起步較晚,但發(fā)展迅速,已在南海北部圈定出了神狐、東沙及瓊東南等多個富含水合物的有利區(qū)塊,初步查明了南海北部水合物資源量。南海北部水合物試采也取得了重大突破性進(jìn)展。水合物資源勘查與商業(yè)試采已被列為我國能源發(fā)展的重大戰(zhàn)略目標(biāo)。盡管如此,較之于發(fā)達(dá)國家,我國在海洋水合物勘查技術(shù)方面仍然存在著一定的差距?!八街梢怨ビ瘛?,筆者希望本文對大洋鉆探水合物測井解釋進(jìn)展的回顧以及發(fā)展方向的探討,能為我國開展相關(guān)研究、發(fā)展具有自主知識產(chǎn)權(quán)的海洋水合物儲層測井評價技術(shù)提供參考與借鑒。

        2 大洋鉆探水合物測井概況

        深海測井始于1968年DSDP執(zhí)行的第1個航次[23]。1979年4月,DSDP66航次進(jìn)行了首次深海水合物測井試驗,但因鉆頭釋放故障未獲成功[11]。同年6月,DSDP67航次在497站位裸眼測井獲得成功,獲取了最早的包括自然伽馬、地層密度、井徑等在內(nèi)的深海水合物測井資料[9]。1982年,DSDP84航次在危地馬拉活動陸緣獲得了首個經(jīng)取心證實的、有明確結(jié)論的深海水合物測井記錄[10]。1989年,高分辨率FMS地層微電阻率掃描成像測井技術(shù)首次在大洋鉆探水合物站位獲得應(yīng)用(ODP127航次[17,24])。1996年,隨鉆測井技術(shù)在ODP170航次首次應(yīng)用于深海水合物鉆探[12](圖1)。迄今為止,在53個鉆遇水合物的大洋鉆探站位中,共有30個站位采集了不同類型的測井?dāng)?shù)據(jù),其中有15個站位還采集了隨鉆測井?dāng)?shù)據(jù)(圖2)。

        大洋鉆探所采用的測井儀器主要有標(biāo)準(zhǔn)測井和隨鉆測井兩大系列[25]。標(biāo)準(zhǔn)測井儀器包括三組合測井儀(Triple Combo)和FMS(地層微電阻率掃描成像測井)—聲波組合測井儀(FMS-Sonic)兩個類別。其中,三組合測井儀由惡劣環(huán)境自然伽馬測井儀(HNGS)、高分辨率陣列側(cè)向測井儀(HRLA)和相量雙感應(yīng)—球形聚焦電阻率測井儀(DIT)、中子孔隙度測井儀(APS)、惡劣環(huán)境巖性密度測井儀(HLDS)及磁化率測井儀(MSS)中的任意3種儀器組合而成,能提供包括自然伽馬和U、Th、K含量,淺、中、深電阻率,中子孔隙度,地層密度和光電吸收截面指數(shù),以及磁化率等常規(guī)測井參數(shù)。FMS—聲波測井儀由地層微電阻率掃描成像測井儀(FMS)、偶極聲波成像儀(DSI-2)、通用目的井斜測井儀(GPIT)及惡劣環(huán)境自然伽馬測井儀(HNGS)等組成,主要提供井壁地層高分辨率微電阻率圖像,縱、橫波速度,井斜方位及自然伽馬能譜等參數(shù)。大洋鉆探隨鉆測井儀器主要由隨鉆復(fù)合電阻率測井儀(geoVISION)、隨鉆補(bǔ)償陣列電阻率測井儀(arcVISION)、隨鉆補(bǔ)償方位密度和中子測井儀(adnVISION)、隨鉆多極子聲波測井儀(sonicVISION)、隨鉆核磁共振測井儀(proVISION)及多功能隨鉆測井儀(EcoScope或NeoScope)等組成,能提供井壁地層的自然伽馬能譜,多種測量方式與探測深度的電阻率,地層密度和光電吸收截面指數(shù),縱、橫波時差和聲波全波列,中子孔隙度及核磁共振等測井參數(shù)(圖3)。

        圖3 國際科學(xué)大洋鉆探計劃所使用的主要測井儀器和測井參數(shù)的垂向分辨率匯總圖

        標(biāo)準(zhǔn)測井與隨鉆測井的主要區(qū)別在于:標(biāo)準(zhǔn)測井在完鉆后進(jìn)行,隨鉆測井則將儀器組合在鉆具上,與鉆井同步進(jìn)行。前者的數(shù)據(jù)質(zhì)量容易受鉆完井液侵入、井壁垮塌和井孔不規(guī)則變形等因素的影響;在海底淺層數(shù)十米至100多米以內(nèi)的套管井段無法測量;局部井孔變形嚴(yán)重時還會導(dǎo)致下伏井段測井工作無法進(jìn)行,所以標(biāo)準(zhǔn)測井?dāng)?shù)據(jù)很難提供覆蓋整個水合物穩(wěn)定帶的連續(xù)觀測記錄。后者則不受上述因素的制約,由于儀器直接綁定在鉆具上,隨鉆測井實現(xiàn)了真正意義上的全井段測井,而且由于測井在地層鉆開后馬上進(jìn)行,鉆完井液污染和水合物分解等因素對測井?dāng)?shù)據(jù)質(zhì)量的影響程度可以降至最低。因此,隨鉆測井比標(biāo)準(zhǔn)測井能更加真實地反映水合物及其宿主地層的原位特征。此外,隨鉆測井中很多傳感器隨著鉆頭的旋轉(zhuǎn)而旋轉(zhuǎn),可以實現(xiàn)360°掃描測量。這些數(shù)據(jù)經(jīng)過處理后可以轉(zhuǎn)化為覆蓋整個井壁的各種成像測井圖像,主要包括自然伽馬,淺、中、深方位電阻率,地層密度,光電吸收截面指數(shù),視中子孔隙度及井徑等。這些隨鉆成像測井圖像的分辨率很高,多介于3~15 cm(圖3),為研究井壁地層的巖性、結(jié)構(gòu)、層理和沉積構(gòu)造、斷層和裂縫、水合物的產(chǎn)狀和分布等特征提供了重要的材料。相比之下,除FMS外,標(biāo)準(zhǔn)測井方法僅能提供分辨率較低(多數(shù)大于30 cm,圖3)的一維測井曲線。

        3 含水合物沉積物的測井響應(yīng)

        水合物具有一系列獨(dú)特的物理和化學(xué)性質(zhì),包括不導(dǎo)電、聲波速度較高、含氫量大等。這就為根據(jù)測井資料識別水合物、預(yù)測其豐度提供了可能[26-28]。

        電阻率測井為識別水合物提供了一個重要的依據(jù)。模擬實驗表明,水合物的存在會導(dǎo)致地層電阻率增加1~3個數(shù)量級[26]。因此,含水合物沉積物的電阻率通常要明顯高于其鄰近不含水合物的沉積物(圖4)。在鉆井過程中,由于水合物分解可能會導(dǎo)致井徑擴(kuò)大和鉆完井液侵入,在標(biāo)準(zhǔn)測井曲線上,含水合物沉積物的深、淺電阻率可能會出現(xiàn)分離。而隨鉆測井則不同,由于水合物來不及分解,隨鉆深、淺電阻率曲線通常重疊在一起[28]。在FMS及各種隨鉆電阻率成像測井圖像上,含水合物沉積物表現(xiàn)為白色高阻特征[15-16,29](圖4)。

        聲波測井為識別水合物提供了另一項重要的依據(jù)。水合物的存在會導(dǎo)致沉積物聲波速度明顯增大(圖4)。模擬實驗表明,含水合物沉積物的聲波速度可能高出同類含水層80%以上[26]。飽和水合物的粗粒沉積物,其縱波速度要明顯高于飽和水或氣體的同類沉積物。含水合物粗粒沉積物的縱波速度要明顯高于細(xì)粒沉積物[30]。此外,水合物的產(chǎn)狀也會對聲波速度產(chǎn)生影響,尤其是以膠結(jié)物形式產(chǎn)出的水合物對縱波速度影響較大,而呈“漂浮”狀散布在孔隙中的水合物對速度的影響則較小。這就為根據(jù)聲波速度預(yù)測水合物的生長習(xí)性提供了可能[30]。

        含水合物沉積物在密度、中子孔隙度及核磁共振(NMR)等其他測井資料上也有著比較明顯的響應(yīng)(圖4)。由于水合物的密度(在0.9 g/cm3左右)比水低,較之于相同條件下的含水層,含水合物層的密度較低。中子孔隙度反映地層的含氫量,含氫量越高,視中子孔隙度越大。理論計算結(jié)果表明,單位體積水、甲烷、Ⅰ型水合物和Ⅱ型水合物的含氫量分別為[31]:6.70×1022、0.01×1022、7.18×1022、7.55×1022氫原子/cm3。顯然,水合物的含氫量為最高,其次為水,而甲烷則最低。所以,在孔隙度一定的情況下,含水合物沉積物的視中子孔隙度要高于水層和游離氣層。與中子孔隙度測井類似,NMR測井也是對地層含氫量的響應(yīng)。但不同于中子孔隙度測井,NMR測井僅對液體中的氫反應(yīng)靈敏[28,32]。因此含水合物沉積物的NMR孔隙度要低于中子孔隙度(圖4-b)。

        盡管如此,由于水合物的密度和含氫指數(shù)與水接近,而電阻率則與油氣類似,所以單純根據(jù)密度或中子測井不能很好地區(qū)分水合物和水層,而單純根據(jù)電阻率測井又很難將水合物與油氣層分開。因此,要準(zhǔn)確地識別水合物,需要根據(jù)多種測井資料來綜合判定。

        4 測井儲層參數(shù)計算

        大洋鉆探水合物儲層測井評價基本上沿用了油氣勘探儲層測井評價的方法。

        含水合物沉積物的孔隙度主要根據(jù)密度、聲波或中子孔隙度測井方法獲得,但在具體計算時需要針對水合物的存在并考慮含水合物沉積物的未固結(jié)或弱固結(jié)特性進(jìn)行適當(dāng)?shù)男U4送?,由于水合物多發(fā)育于海底淺部未固結(jié)或弱固結(jié)沉積物中,井孔條件較差時(如大井徑或井徑劇烈變化),各種孔隙度測井資料的質(zhì)量會變差,此時可以考慮根據(jù)電阻率測井曲線利用阿爾奇公式估算沉積物的孔隙度[33]。

        水合物飽和度主要根據(jù)電阻率測井資料利用阿爾奇公式加以計算。20世紀(jì)90年代中后期以來,密度—NMR測井聯(lián)合方法、各種形式的三相聲波方程及建立在各種巖石物理模型基礎(chǔ)上的速度正演方法等也逐漸應(yīng)用于大洋鉆探水合物飽和度的計算(表1、2)。

        圖4 含天然氣水合物沉積物的測井響應(yīng)圖

        4.1 阿爾奇公式法

        Archie[34]建立了水飽和的未固結(jié)砂或砂巖地層的電阻率、孔隙度與孔隙水電阻率之間的經(jīng)驗關(guān)系(表1)。當(dāng)沉積物中含有水合物時,其電阻率值會升高。假定電阻率的增加全部由水合物的存在而引起,則可以根據(jù)阿爾奇公式計算水合物飽和度[3,36,65](表1)。標(biāo)準(zhǔn)阿爾奇公式涉及孔隙度、孔隙水電阻率及經(jīng)驗參數(shù)a、m和n的求取(表1)。其中,孔隙度可以根據(jù)密度、聲波或中子孔隙度測井求得,也可以根據(jù)巖心測試結(jié)果擬合得到的孔隙度隨深度變化趨勢進(jìn)行估算。原位孔隙水的電阻率可以根據(jù)孔隙水鹽度和地溫梯度估算[66]。a和m分別被稱為“彎曲度系數(shù)”和“膠結(jié)指數(shù)”,是反映沉積物孔隙結(jié)構(gòu)特征的兩個經(jīng)驗參數(shù),一般根據(jù)完全水飽和沉積物的地層因素(水飽和地層的電阻率與孔隙水電阻率之比)與孔隙度在雙對數(shù)交會圖上通過線性擬合確定,其中最優(yōu)擬合直線的斜率代表m,該直線在地層因素軸上的截距(孔隙度=1)為a[36,38,67]。Archie[34]給出的固結(jié)砂巖的m值介于1.8~2.0,未固結(jié)砂的m值在1.3左右,部分固結(jié)砂巖的m值介于1.3~2.0。n為飽和度指數(shù),可以用不同巖性的平均值代替[32,39,43];Malinverno等[36]采用巖心飽和度標(biāo)定方法確定n的值,Hyndman等[65]則假定n=m;n的取值還可用缺省值2代替[67]。此外,Archie[34]認(rèn)為,純凈的未固結(jié)砂和固結(jié)砂巖的n值一般在2左右,并由此提出了簡化的阿爾奇公式,亦稱快速評價阿爾奇公式(表1)。當(dāng)孔隙水電阻率等參數(shù)未知時,可以用該公式來估算水合物飽和度。

        表1 根據(jù)測井資料計算水合物飽和度的方法統(tǒng)計表

        表2 建立在各種巖石物理模型模量計算基礎(chǔ)上的估算水合物飽和度并預(yù)測水合物生長習(xí)性的速度正演方法統(tǒng)計表

        續(xù)表2-1

        續(xù)表2-2

        阿爾奇公式已被廣泛應(yīng)用于大洋鉆探水合物飽和度計算。Collett等[35]根據(jù)標(biāo)準(zhǔn)和簡化阿爾奇公式對布萊克海脊994、995、997三個站位(ODP164航次)的水合物飽和度進(jìn)行了估算,其結(jié)果介于2%~11%。Guerin等[33]根據(jù)簡化阿爾奇公式計算得到的995站位水合物飽和度介于5%~10%,與Collett等[35]的結(jié)果相近。不過,Guerin等[33]認(rèn)為,由于鉆井過程中井壁附近水合物分解會使孔隙流體鹽度降低,從而導(dǎo)致常規(guī)電阻率測井值增大。因此根據(jù)常規(guī)電阻率測井資料計算得到的水合物飽和度代表了水合物飽和度的最大估計值。Tréhu等[15]利用阿爾奇公式對水合物脊鉆探的1245、1248、1249、1251 四個站位(ODP204航次)的水合物飽和度進(jìn)行了計算,并將計算結(jié)果與巖心熱紅外掃描成像、孔隙水氯離子濃度及沉積學(xué)觀察等資料所揭示的水合物分布進(jìn)行了對比,結(jié)果具有良好的一致性,其中1245站位計算的水合物飽和度介于10%~30%;1248站位計算的水合物飽和度介于0~58%,但除最上部20 mbsf塊狀水合物發(fā)育段外, 20 mbsf 以下水合物飽和度一般僅為2%~5%;1251站位計算的水合物飽和度為18%;1249站位計算的水合物飽和度介于10%~92%,是該航次水合物豐度最高的站位。Riedel等[16]根據(jù)隨鉆測井資料認(rèn)為,Cascadia 邊緣北部 U1325、U1326、U1327、U1328四個站位(IODP311航次)中,高濃度的水合物主要分布于高電阻率、高聲波速度和低密度薄砂層內(nèi),并根據(jù)阿爾奇公式計算得出這些薄砂層內(nèi)水合物飽和度高達(dá)80%~95%。

        阿爾奇公式法存在的主要問題如下:①該公式建立在不含泥質(zhì)的純凈砂和凈砂巖基礎(chǔ)之上,此時地層的導(dǎo)電性主要受孔隙流體的導(dǎo)電性控制。當(dāng)?shù)貙又泻写罅筐ね恋V物時,需要對黏土礦物表面電荷引起的附加導(dǎo)電效應(yīng)進(jìn)行修正[33,67-68]。不過,也有觀點(diǎn)認(rèn)為,對于孔隙度很高的海洋泥質(zhì)沉積物而言,這種附加導(dǎo)電效應(yīng)并不明顯,阿爾奇公式可以無須修正而直接使用[69-70]。②阿爾奇公式僅對水濕巖石有效,適用于水合物“漂浮”在孔隙中間,即水合物的存在不影響顆粒表面水濕特性的情形[28]。當(dāng)水合物以粒間膠結(jié)物或結(jié)核狀、脈狀、塊狀等形式產(chǎn)出時,可能會出現(xiàn)較大的偏差[67]。③阿爾奇公式所需的經(jīng)驗參數(shù)a、m、n并不總是很容易求得的,而且這些參數(shù)的可靠性很少經(jīng)過嚴(yán)格的標(biāo)定和論證[28,36,38,43,67]。以布萊克海脊889/890站位為例,Hyndman等[65]利用巖心樣品測試結(jié)果擬合確定的a和m值分別為:a=1.41、m=1.76;而Collett[32]根據(jù)中子測井和電阻率測井?dāng)?shù)據(jù)擬合得到的參數(shù)值分別為:a=0.967、m=2.81。顯然這兩組參數(shù)之間存在著較大的偏差,而這些經(jīng)驗參數(shù)取值上的差異有可能導(dǎo)致水合物飽和度估算結(jié)果相差很大。

        4.2 密度—NMR測井聯(lián)合法

        密度—NMR測井聯(lián)合法或DMR法[28,39,40-41](表1)利用NMR測井反映液體含氫量這一特性,將視NMR測井孔隙度(水合物是固體骨架的一部分)與視密度測井孔隙度(水合物是孔隙空間的一部分)聯(lián)合,求解沉積物的總孔隙度和水合物飽和度。該方法假定,完全水飽和沉積物中測得的NMR回波幅度反映總孔隙度。當(dāng)沉積物中含有很多小孔隙(如細(xì)粒富含黏土的沉積物)時,T2弛豫時間會短于NMR儀器的探測極限,此時NMR孔隙度將低估沉積物的孔隙度。這種情況容易出現(xiàn)在細(xì)粒富含黏土沉積物的長死時間隨鉆NMR測井中。常規(guī)NMR測井儀的死時間要短得多,可以用來測定泥質(zhì)沉積物的總孔隙度。DMR方法的誤差來源主要有兩個:①骨架密度不準(zhǔn)引起的密度孔隙度計算誤差;②當(dāng)沉積物含大量黏土礦物束縛水時,T2弛豫時間會加快,視NMR孔隙度將被低估,從而影響水合物飽和度計算結(jié)果的準(zhǔn)確度[28]。

        Collett等[41]假定NMR孔隙度相對于密度孔隙度的下降全部由水合物的存在而引起,提出了一個根據(jù)密度和NMR測井孔隙度計算水合物飽和度的簡單公式(表1)。他們利用該公式,根據(jù)ODP204航次在水合物脊采集的隨鉆NMR和隨鉆密度測井孔隙度資料,計算出了該航次全部9個站位水合物穩(wěn)定帶內(nèi)水合物的飽和度,其結(jié)果介于-10%~50%,多數(shù)站位介于20%~30%,總體上低于根據(jù)阿爾奇公式計算得出的水合物飽和度。此外,在非水合物層段內(nèi)計算的水合物飽和度介于-10%~10%,可能代表了該方法估算水合物飽和度時的誤差范圍[41]。

        4.3 各種形式的三相聲波方程法

        根據(jù)聲波測井計算含水合物沉積物的孔隙度和飽和度時,通常將水合物作為獨(dú)立組分,采用各種形式的三相(礦物骨架+水合物+地層水)聲波方程,包括三相Wyllie時間平均方程[42-43]、三相Wood方程[45-46]及三相時間平均—Wood加權(quán)方程[45,47](表1)。Wyllie時間平均方程是根據(jù)固結(jié)砂巖總結(jié)出來的經(jīng)驗方程,而Wood方程描述的則是懸浮顆粒的速度與孔隙度的關(guān)系。應(yīng)用上述兩個方程估算深海未固結(jié)或弱固結(jié)且含大量黏土及有機(jī)質(zhì)的海洋沉積物的速度時,均存在著較大的誤差,其中時間平均方程的估計值偏高,而Wood方程的估計值則偏低。為此,Nobes等[71]提出,可以將Wyllie時間平均方程和Wood方程結(jié)合起來,采用加權(quán)平均方法計算海洋沉積物的縱波速度,進(jìn)而提出了所謂的兩相加權(quán)方程。Lee等[47]將該方程拓展至三相,提出了三相加權(quán)方程(表1)。但這類加權(quán)平均方程本身并沒有明確的物理意義,而且Wyllie時間平均方程和Wood方程的權(quán)系數(shù)如何設(shè)定也沒有明確的依據(jù)。此外,考慮到深海含水合物沉積物中黏土含量較高,Lee等[47]還建議,進(jìn)一步將黏土從骨架中分離出來作為一個獨(dú)立組分,提出了所謂的四相聲波方程,包括四相Wyllie時間平均方程、四相Wood方程和四相加權(quán)方程,但由于所涉及的未知變量過多,四相聲波方程在實際中很少得到應(yīng)用。Lee[72]用三相加權(quán)方程根據(jù)聲波測井資料對布萊克海脊994、995、997三個站位(ODP164航次)的水合物飽和度進(jìn)行了估算。計算時分別采用了三種不同的孔隙度資料,包括密度測井孔隙度、巖心孔隙度及巖心孔隙度線性擬合趨勢。計算得到的平均水合物飽和度介于3.8%(巖心孔隙度)至12.1%(使用體積密度孔隙度)之間。

        4.4 基于巖石物理模型的彈性波速度模擬估算水合物飽和度

        沉積物的縱、橫波速度取決于體積模量、剪切模量和密度等彈性參數(shù)(表2)。這些彈性參數(shù)與沉積物的組成、結(jié)構(gòu)、孔隙度、水合物的產(chǎn)狀和飽和度等具有密切的關(guān)系。沉積物的密度可以由密度測井獲得;沉積物的組成可以根據(jù)薄片鑒定及X射線衍射礦物成分分析等方法確定;孔隙度可以根據(jù)測井資料或巖心樣品測試獲得。給定含水合物沉積物的基本物理性質(zhì)和相關(guān)參數(shù),可以通過巖石物理模型計算得到不同飽和度取值時沉積物的彈性模量和速度;然后將模擬計算得到的速度與實際聲波測井速度進(jìn)行對比,二者最吻合時,計算彈性模量所采用的飽和度值便代表了水合物飽和度的預(yù)測結(jié)果。因此,選擇合適的能夠代表實際地質(zhì)特征的巖石物理模型,是通過速度建模估算水合物飽和度的關(guān)鍵。

        圖5 根據(jù)縱波速度和水合物飽和度預(yù)測水合物生長習(xí)性圖[72]

        迄今已提出的一系列用于計算含水合物沉積物體積模量和剪切模量的巖石物理模型,大致可以劃分為3類(表2):①接觸與膠結(jié)理論模型[49-52];②有效介質(zhì)理論模型[47,51,53-54,55-56];③考慮礦物顆粒、水合物和流體相互作用的理論模型[57-58,59-64]。不同的巖石物理模型在計算含水合物沉積物的模量時,通常將水合物在沉積物中的賦存狀態(tài)概括為以下4種概念模型[73](圖5):①水合物是孔隙流體的一部分,它們對沉積物骨架的彈性響應(yīng)幾乎不產(chǎn)生實質(zhì)性影響;②水合物以顆粒形式產(chǎn)出,與沉積物顆粒一樣是沉積物骨架的一部分,共同起到應(yīng)力支撐作用;③水合物以膠結(jié)物形式在顆粒接觸部位生長(接觸式膠結(jié)物);④水合物以膠結(jié)物形式包繞顆粒表面生長(顆粒襯邊式膠結(jié)物)。不同的巖石物理模型對水合物產(chǎn)狀的處理存在著一定的差異,比如接觸與膠結(jié)模型將水合物視為膠結(jié)物,以接觸式或襯邊式產(chǎn)出;而有效介質(zhì)理論則將水合物視為孔隙充填物或巖石骨架的一部分,作為流體的一部分或以水合物顆粒形式產(chǎn)出。

        接觸與膠結(jié)模型一般將沉積物顆粒近似視為理想的彈性球的集合體,顆粒介質(zhì)的彈性模量主要是由顆粒之間的接觸方式、顆粒接觸面積、接觸點(diǎn)數(shù)、球體尺寸、介質(zhì)彈性參數(shù)以及所施加的壓力來決定的。在含水合物沉積物中,水合物被視作接觸膠結(jié)物時可以被看成是很薄的彈性介質(zhì)[49],其膠結(jié)面積小于礦物顆粒的面積,利用彈性力學(xué)方法分析膠結(jié)物和顆粒之間的相互作用可以得出其相對應(yīng)的法向剛度和剪切剛度,并結(jié)合孔隙度與配位數(shù)計算水合物為接觸式膠結(jié)時的沉積物骨架的體積模量和剪切模量(表2)[50]。而當(dāng)水合物被視為襯邊式膠結(jié)物時,Dvorkin等[51]通過改變膠結(jié)層半徑與礦物顆粒半徑的比值,利用相同的方法可以計算出含水合物沉積物骨架的彈性模量(表2)。礦物顆粒的體積模量和剪切模量通常根據(jù)Viogt-Ruess-Hill平均計算,而得到沉積物骨架的彈性模量后,可以通過Gassmann方程來計算流體對整體含水合物沉積物彈性性質(zhì)的影響。Guerin等[33]根據(jù)膠結(jié)理論對ODP164航次955站位的水合物飽和度進(jìn)行估算,當(dāng)水合物為接觸式膠結(jié)時,估算的水合物飽和度范圍為0~5%;而當(dāng)水合物為襯邊式膠結(jié)時,估算的水合物飽和度范圍為0~10%,后者與保壓巖心樣品測量結(jié)果相吻合。接觸與膠結(jié)理論雖然能夠較好的刻畫水合物作為膠結(jié)物形式對彈性特征的影響,但其缺陷在于不能有效地解釋速度隨圍壓的變化[51]。

        表征含水合物沉積物彈性特征的有效介質(zhì)理論模型主要基于界限模型(如Hashin-Shtrikman邊界等)和包體模型(如Kuster-Toks?z模型、自洽模型、微分等效介質(zhì)模型等)(表2)。這類模型主要刻畫水合物作為顆粒骨架或孔隙充填物時對沉積物整體彈性性質(zhì)的貢獻(xiàn)。根據(jù)含水合物沉積物在超過臨界孔隙度和低于臨界孔隙度時的力學(xué)結(jié)構(gòu)特征,Dvorkin等[53]利用改進(jìn)的Hashin-Shtrikman上邊界與下邊界[51]建立了分段函數(shù),分別計算這兩個地質(zhì)狀態(tài)下的沉積物骨架彈性模量(表2)[54]。在該模型中,當(dāng)水合物作為流體的一部分時,它們不影響礦物顆粒的體積模量和剪切模量,只影響沉積物中流體的體積模量;當(dāng)水合物以顆粒形式產(chǎn)出時,它們是沉積物骨架的一部分,影響礦物顆粒的體積模量和剪切模量,并改變干沉積物的等效彈性性質(zhì),進(jìn)而最終改變含水合物沉積物的速度特征(表2)。Chand等[74]假設(shè)水合物以顆粒形式產(chǎn)出,根據(jù)有效介質(zhì)理論估算ODP164航次995和997站位的水合物飽和度,其結(jié)果介于10%~20%,與根據(jù)電阻率測井估算的飽和度結(jié)果相當(dāng)。該理論模型在使用過程中需要注意的是,其所采用的臨界孔隙度(36%~40%)與砂巖相當(dāng)[75],當(dāng)沉積物中含有不同數(shù)量黏土?xí)r會導(dǎo)致配位數(shù)出現(xiàn)很大的變化[76]。

        Kuster-Toks?z方程(以下簡稱K-T方程)根據(jù)基質(zhì)相和包含相的彈性模量以及包含物的濃度,利用波在兩相介質(zhì)中的傳播理論計算沉積物的體積模量和剪切模量(表2)[55]。Zimmerman等[56]將其應(yīng)用于加拿大北極的巖心樣本并預(yù)測縱、橫波速度;Lee等[57]將其擴(kuò)展到含水合物沉積物中預(yù)測水合物飽和度。利用K-T方程計算含水合物沉積物的體積模量與剪切模量可以分兩步進(jìn)行:首先,將孔隙水視為基質(zhì)相中的球體包含物,根據(jù)K-T方程計算得到由水合物和孔隙水合并組成的新基質(zhì)相的體積模量和剪切模量;然后,將沉積物顆粒視為由水合物和孔隙水合并組成的新的基質(zhì)相中的球體包含物,計算得到含水合物沉積物的體積模量和剪切模量。Guerin等[33]根據(jù)K-T方程估算ODP164航次955站位的水合物飽和度時,假定含水合物沉積物由水合物、水和黏土組成,計算得到該站位水合物飽和度介于10%~20%,遠(yuǎn)高于保壓巖心樣品測量值。與大多數(shù)基于理想橢球體的包體理論一樣,K-T方程假設(shè)包體相內(nèi)的應(yīng)力分布是獨(dú)立的,并且沒有考慮到不同包體之間的相互作用。這與實際沉積物中水合物與顆粒之間復(fù)雜的力學(xué)接觸狀態(tài)往往相違。與水合物以膠結(jié)物形式產(chǎn)出相比,要獲得相同的模量增加量,K-T方程所需要的水合物的體積會更大[33],故根據(jù)K-T方程估算得到的水合物飽和度往往偏高[47]。此外,K-T方程將沉積物中的包含物如石英、黏土等視為懸浮體,這一假設(shè)前提僅在沉積物孔隙度很高(超過80%)時才成立[47]。

        由于水合物特殊的賦存狀態(tài)和彈性特征,上述第三類理論模型主要考慮的是礦物顆粒、水合物和流體之間的相互作用對含水合物沉積物整體彈性性質(zhì)的影響,這方面代表性的工作包括改進(jìn)的Biot-Gassmann(BGTL)和三相Biot方程(TPBE,表2)。Lee[58]假定沉積物速度比(橫波速度比縱波速度)取決于骨架的速度比和孔隙度,給出了根據(jù)骨架剪切模量、Biot系數(shù)、速度比關(guān)系和孔隙度計算沉積物剪切模量的方法(表2)。改進(jìn)的Biot-Gassmann模型利用Biot[77]和Gassmann[78]方法計算沉積物的體積模量,利用Lee[58]的方法計算沉積物的剪切模量。Biot系數(shù)取決于孔隙度,在孔隙度為零時,如果已知骨架的性質(zhì)、沉積物的剪切模量或橫波速度與孔隙度的關(guān)系,可以通過實際橫波速度與計算得到的橫波速度擬合獲得Biot系數(shù),其中沉積物的剪切模量可以根據(jù)加權(quán)方程或有效介質(zhì)理論計算。由于速度比與孔隙度、差異壓力、固結(jié)程度和黏土含量等因素有關(guān),Lee[58]引入指數(shù)n減少差異壓力對速度比的影響,Lee[57]根據(jù)Prasad[79]匯編的實驗室數(shù)據(jù)(頻率為100 kHz~1 MHz)導(dǎo)出了指數(shù)n與差異壓力的關(guān)系式,并引入系數(shù)G消除黏土對縱、橫波速度比的影響。Han等[80]通過實驗室數(shù)據(jù)導(dǎo)出了G與黏土含量的關(guān)系式。Lee[57,81]、Lee等[82]利用該方法,根據(jù)標(biāo)準(zhǔn)測井橫波速度(縱波速度易受游離氣存在的影響)對ODP204航次1244( 76~127 mbsf )、1245( 73~129 mbsf)和 1247(74~129 mbsf)3個站位的水合物飽和度進(jìn)行了估算,平均值分別為10.2%、10.4%、6.1%。根據(jù)標(biāo)準(zhǔn)測井橫波速度計算得到的飽和度在變化趨勢上與根據(jù)隨鉆電阻率測井估算結(jié)果相當(dāng),但橫波速度估算結(jié)果總體要高于電阻率估算結(jié)果。BGTL方法最適合于超聲波頻率下聲波測井測得的無明顯頻散的速度[57]。該方法存在的問題在于利用加權(quán)方程或有效介質(zhì)理論擬合Biot系數(shù)時,其對應(yīng)的問題也存在于改進(jìn)的Biot模型中。

        三相Biot方程(表2)基于滲流理論將Biot理論推廣到由沉積物顆粒、孔隙水和水合物組成的均勻的、相互交織的骨架中,假設(shè)每種骨架的體積模量和剪切模量均對沉積物的縱、橫波速度產(chǎn)生影響,根據(jù)廣義的Biot彈性系數(shù)推導(dǎo)出沉積物的體積模量和剪切模量[59-60]。其中,沉積物骨架和水合物骨架的體積模量和剪切模量由Pride等[61]和Lee[62]提出的干沉積物骨架公式計算(表2)[63-64]。固結(jié)參數(shù)α可以根據(jù)不含水合物沉積物的速度與深度曲線擬合獲得。參數(shù)ε是衡量水合物形成對礦物骨架支撐程度的經(jīng)驗參數(shù)。Lee等[64]用三相Biot理論模擬了IODP311航次1326D孔含水合物沉積物的聲波速度并估算水合物飽和度,其結(jié)果與根據(jù)核磁共振和電阻率測井估算的飽和度值基本吻合。Miyakawa等[83]根據(jù)聲波測井?dāng)?shù)據(jù)利用三相Biot理論計算得到C0002站位(IODP314/315航次)水合物的飽和度介于0~30%,與根據(jù)阿爾奇公式計算得到的水合物飽和度(0~35%)相當(dāng)??傮w上,TPBE方程更適合于含水合物的砂層,而不是黏土含量較高的泥質(zhì)沉積,因為Biot理論的假設(shè)前提是不含束縛水[64]。

        此外,有些學(xué)者根據(jù)含水合物沉積物的彈性和電阻率巖石物理模型,開展了基于聲波和電阻率數(shù)據(jù)聯(lián)合反演水合物飽和度的相關(guān)研究。Sava等[84]提出了在貝葉斯理論框架下,根據(jù)阿爾奇方程和水合物作為承載分量的巖石物理彈性模型來反演水合物飽和度的方法。此方法在墨西哥灣綠色峽谷含水合物井的應(yīng)用結(jié)果表明,與單獨(dú)使用電阻率或速度數(shù)據(jù)相比,使用聯(lián)合反演的方法提高了計算水合物飽和度的準(zhǔn)確度。Lee[85]在最小二乘框架下提出了聯(lián)合聲波慢度(縱、橫波慢度)和電阻率測井?dāng)?shù)據(jù)反演水合物飽和度的方法,并將其應(yīng)用于加拿大麥肯齊三角洲(Mackenzie Delta)Mallik 2L-38 研究井的水合物飽和度估算。Pan等[86]提出利用縱波阻抗、泊松比和對數(shù)電阻率數(shù)據(jù)并建立相應(yīng)的三維巖石物理模板來反演水合物飽和度的方法,并應(yīng)用于ODP204航次1247和1250站位,其反演結(jié)果和巖心數(shù)據(jù)測得的水合物飽和度結(jié)果吻合較好。相對而言,僅用電阻率估算的結(jié)果偏低,而用單一彈性屬性估算的結(jié)果則偏高。

        5 水合物分布與宿主沉積物巖性之間的相關(guān)性

        海底水合物穩(wěn)定帶內(nèi),水合物的分布受多種因素的制約,包括低溫、高壓環(huán)境和充足的形成水合物的氣體來源等。一般認(rèn)為,在氣源供應(yīng)充足的前提下,宿主沉積物巖性對水合物的形成和分布具有重要的控制作用[87]。巖心觀察結(jié)果表明,海底水合物穩(wěn)定帶內(nèi),水合物的分布具有明顯的非均質(zhì)性,主要表現(xiàn)在水合物分布對宿主沉積物巖性的選擇性上,也體現(xiàn)在水合物在宿主沉積物內(nèi)部對成核和生長部位的選擇性上。前者可以稱之為“層間非均質(zhì)性”,后者可以稱之為“層內(nèi)非均質(zhì)性”。

        一般認(rèn)為,水合物傾向于分布在粒度較粗的宿主沉積物中。Ginsburg等[88]對布萊克海脊533站位(DSDP76航次)及994、995和997站位(ODP164航次)沉積物的粒度分析結(jié)果表明,盡管該海脊水合物穩(wěn)定帶內(nèi)的沉積物基本上由均一的半遠(yuǎn)洋細(xì)粒沉積物(粉砂質(zhì)泥和泥質(zhì)粉砂)組成,但含水合物的沉積物(對應(yīng)于孔隙水氯離子負(fù)異常帶)中,粗組分(粉砂)的含量要略高于不含水合物的沉積物,或粒度小于 0.05 mm、0.01 mm、0.005 mm、0.001 mm 的細(xì)組分含量略低于不含水合物的沉積物。Weinberger等[89]將ODP204航次根據(jù)巖心熱紅外掃描成像記錄所確定的水合物分布位置與沉積物巖性剖面進(jìn)行對比后發(fā)現(xiàn),該航次所有站位的水合物分布均與厚度不足0.5 cm的薄砂層有關(guān)。Pi?ero等[90]對ODP204航次在1245、1248、1249三個站位采集的湯狀和奶油狀含水合物沉積物樣品做了進(jìn)一步的粒度、礦物組分和磁化率分析。結(jié)果表明,含水合物沉積物的平均粒徑較大、磁化率高,對應(yīng)于粗粒濁流沉積。溫哥華岸外 Cascadia 大陸邊緣 U1325、U1326、U1327等站位(IODP311航次)水合物主要分布于粒度較粗、單層厚度一般小于5 cm(最大值可達(dá)23 cm)的濁積砂中[87,91]。這些薄砂層中水合物飽和度可以高達(dá)50%~80%[87,91-92],但在分隔這些薄砂層、平均厚度高達(dá)2.5 m的細(xì)粒泥質(zhì)沉積物(黏土和粉砂質(zhì)黏土)中卻很少或不含水合物[91]。水合物分布于較粗粒沉積物中的實例還有:①墨西哥灣北部陸緣618站位(DSDP96航次)20~40 mbsf深度段內(nèi),巖心觀察到的1~2 mm直徑的球狀水合物主要分布于砂質(zhì)層段[5];②哥斯達(dá)黎加邊緣1040和1041站位(ODP170航次)中的水合物在較粗粒的沉積物和富含火山灰的層段中濃度較高[12];③日本南海海槽活動邊緣增生楔C0002站位(IODP314/315航次)的水合物以膠結(jié)物形式優(yōu)先分布于BSR之上的砂質(zhì)濁積層中[19];④新西蘭東部Hikurangi活動大陸邊緣的U1517站位(IODP372航次)的水合物也優(yōu)先富集于粒度較粗的濁積砂中[20]。關(guān)于水合物優(yōu)先分布于粗粒沉積物中的成因,存在著不同的解釋。Ginsburg等[93-94]認(rèn)為,當(dāng)水溶液從細(xì)粒沉積物進(jìn)入粗粒沉積物時,由于相平衡和化學(xué)動力學(xué)條件發(fā)生改變,水合物容易從水溶液中析出。Torres等[92]認(rèn)為,粗粒沉積層內(nèi)孔隙毛細(xì)管壓力較小,有利于甲烷氣體的運(yùn)移和水合物的形成。Malinverno[91]認(rèn)為,水合物優(yōu)先在粗粒層中發(fā)育的原因還可能與水合物在細(xì)粒沉積物中的形成過程容易被抑制有關(guān)。在僅考慮原地微生物成因甲烷氣源的前提下,Malinverno[91]應(yīng)用物質(zhì)平衡模型模擬了U1325站位(IODP311航次)粗粒、細(xì)粒沉積物交替序列中水合物的形成。結(jié)果表明,該站位細(xì)粒層內(nèi)甲烷的濃度從未達(dá)到水合物形成所必需的過飽和條件。

        另一方面,大洋鉆探證實,水合物在很多細(xì)粒沉積物中也大量的存在[68,95]。最著名例子便是布萊克海脊。該海脊的水合物穩(wěn)定帶內(nèi),沉積物巖性非常均一。雖然前人的研究[88]認(rèn)為,該海脊含水合物的沉積物中粗粒組分的含量稍多一些或細(xì)粒組分含量略低一些,但在這種粒度如此均一的細(xì)粒沉積序列中水合物的分布所表現(xiàn)出來的強(qiáng)烈的非均質(zhì)性,僅用如此細(xì)微的粒度組分含量上的變化來加以解釋,多少還是難以令人信服的。除了布萊克海脊之外,ODP112航次在秘魯活動陸緣鉆探的4個水合物站位中,水合物穩(wěn)定帶內(nèi)的沉積物基本上由細(xì)粒的硅質(zhì)生物軟泥(硅藻泥)和泥組成,水合物在這些泥質(zhì)沉積物中主要呈分散狀分布;在該航次688站位141 mbsf處采到的一枚數(shù)厘米大的水合物樣品也由水合物與泥的混合物組成[7]。此外,ODP204航次在位于水合物脊南高點(diǎn)附近1248站位最上部20 mbsf獲得的塊狀水合物,也主要分布于細(xì)粒粉砂質(zhì)黏土和含硅藻的粉砂質(zhì)黏土中[15]。很顯然,水合物在粗粒和細(xì)粒沉積物中均可以形成。研究表明,粒度粗細(xì)對水合物形成和分布的影響,似乎主要體現(xiàn)在水合物的產(chǎn)狀和飽和度上[68,96-98]。在粗粒沉積物(粗粉砂和砂,如日本南海海槽)中,水合物多呈孔隙充填式,一般具有較高的飽和度(>50%~80%);細(xì)粒沉積物(如布萊克海脊)中的水合物則多呈分散狀且飽和度低(通常不足10%)。此外,細(xì)粒沉積物中的水合物還常呈結(jié)核或透鏡體分布于局部沉積物粒度稍微增大、滲透率相對較高的部位,其飽和度值變化大[88,92,95,99-100]。

        在利用大洋鉆探測井資料研究水合物分布與沉積物巖性相關(guān)性方面已經(jīng)開展了一定量的工作,但多數(shù)僅局限于水合物飽和度的計算。從上文介紹的測井水合物飽和度計算實例來看,海洋水合物穩(wěn)定帶內(nèi),水合物的豐度無論在縱向上(站位內(nèi))還是橫向上(站位間)均表現(xiàn)出了明顯的非均質(zhì)性。與水合物飽和度計算相比,測井巖性預(yù)測方面的工作做得還很少,現(xiàn)有的測井巖性評價研究基本上還停留在簡單地用自然伽馬測井歸一化方法估計沉積物泥質(zhì)含量方面(如 Guerin 等[33]、Lee 等[82])。

        6 天然氣水合物生長習(xí)性預(yù)測

        巖心觀察結(jié)果表明,水合物的產(chǎn)狀主要有塊狀、(紋)層狀、結(jié)核狀或透鏡狀、脈狀及分散狀等類型。這表明水合物在沉積物中并不是均勻地從所有孔隙中析出,而是優(yōu)先從沉積物中的某些部位開始生長。巖心分析結(jié)果進(jìn)一步證明,水合物的成核與生長更容易發(fā)生在沉積物顆粒之間開闊的孔隙空間或孔洞中,而不是更為封閉和狹小的區(qū)域[93,101]。Lorenson[102]利用透射和掃描電鏡對布萊克海脊(ODP164航次)和哥斯達(dá)黎加陸緣(ODP170航次)含水合物沉積物的觀察發(fā)現(xiàn),分散狀水合物主要形成于沉積物中原生或次生孔隙較發(fā)育的部位,包括含有豐富的微體化石(有孔蟲、硅藻等)體腔孔、被松散砂或砂級微體化石碎屑充填的潛穴、未固結(jié)的草莓狀黃鐵礦集合體及疏松多孔狀的火山灰中,而不是缺少這些大孔隙空間的沉積物中。Lorenson[102]認(rèn)為,狹窄孔隙空間中的毛細(xì)管力對水合物形成的溫度—壓力穩(wěn)定場具有抑制作用。Kraemer等[103]根據(jù)布萊克海脊994站位(ODP164航次)沉積物樣品分析發(fā)現(xiàn),該站位185~260 mbsf間的水合物發(fā)育帶內(nèi),沉積物的碳酸鹽含量較低,但硅質(zhì)微體化石含量較高。他們認(rèn)為,硅質(zhì)微體化石豐度的增加使得孔隙的大小和圓度增加;大而圓的孔隙毛細(xì)管力小,有利于水合物的成核。

        聲波測井在預(yù)測海底沉積物中水合物生長習(xí)性方面已經(jīng)有了可喜的開端。上文提到的建立在各種巖石物理模型模量計算基礎(chǔ)之上的速度正演方法,通過模擬速度與聲波測井速度的對比,除了可以估算水合物飽和度外,還可以用于預(yù)測水合物的生長習(xí)性,評估水合物究竟是作為孔隙流體的一部分“懸浮”于孔隙中生長,還是以不同形式的膠結(jié)物生長在顆粒接觸部位或包繞顆粒表面生長(圖5)。水合物生長習(xí)性不同,對沉積物的強(qiáng)度、彈性模量和聲波速度具有不同的影響。實驗證明,當(dāng)水合物以粒間膠結(jié)物形式產(chǎn)出時,即使含量很少,但它們作為沉積物承載系統(tǒng)的一部分也會極大地增加沉積物的硬度;而當(dāng)水合物作為孔隙流體的一部分產(chǎn)出時,水合物的存在對沉積物骨架的硬度沒有影響[104]。速度模擬方法在實際應(yīng)用中也遇到了一些問題,以布萊克海脊995站位(ODP164航次)為例,不同的人采用不同的巖石物理模型得到了不同的水合物生長習(xí)性預(yù)測結(jié)果。如Helgerud等[54]基于有效介質(zhì)理論,分別假定水合物以漂浮式(作為流體的一部分)和承載式(作為骨架的一部分)兩種情形進(jìn)行模擬。結(jié)果表明,漂浮式模擬得到的水合物飽和度明顯偏高,而承載式模擬得到的水合物飽和度則與根據(jù)電阻率測井及孔隙水地化分析得出的水合物飽和度非常吻合,說明該站位水合物以承載式產(chǎn)出,水合物與固體礦物骨架一道構(gòu)成了含水合物沉積物的承載體系。Guerin等[33]考慮了接觸式和襯邊式兩種可能的水合物膠結(jié)形式,采用膠結(jié)理論進(jìn)行模擬,結(jié)果表明,水合物以接觸式膠結(jié)物產(chǎn)出時,計算得到的水合物飽和度過低;但以襯邊式膠結(jié)物形式產(chǎn)出時,計算結(jié)果與保壓巖心測試的水合物飽和度相吻合。Jakobsen等[105]以富含黏土的沉積物為例,綜合應(yīng)用了各向異性自洽理論,和微分等效介質(zhì)理論探討了水合物在懸浮式和膠結(jié)式兩種生長模式下對彈性特征的影響,并將理論模擬結(jié)果與布萊克海脊995站位的實際數(shù)據(jù)進(jìn)行對比,認(rèn)為該站位水合物更傾向于懸浮式生長模式。

        7 存在的問題與展望

        綜上所述,大洋鉆探測井資料在海底含水合物沉積物的識別和水合物儲層評價等方面已經(jīng)取得了長足的進(jìn)展,基本建立起了水合物儲層測井識別依據(jù)和水合物儲層參數(shù)(包括孔隙度和水合物飽和度)計算方法體系,在水合物儲層巖性和非均質(zhì)性評價特別是水合物生長習(xí)性預(yù)測方面獲得了初步應(yīng)用。相關(guān)研究目前存在的主要問題與展望分述如下。

        1)測井資料在表征水合物儲層的巖性和非均質(zhì)性,特別是水合物分布與宿主沉積物之間相互關(guān)系方面的研究程度還較低。海底水合物穩(wěn)定帶內(nèi)水合物的垂向分布規(guī)律還不十分清楚,巖性控制水合物形成和分布的機(jī)理也不夠明確。現(xiàn)有的水合物與宿主沉積物之間關(guān)系的認(rèn)識主要建立在巖心觀察及樣品分析(主要是粒度分析),并與巖心水合物標(biāo)志(包括湯狀/奶油狀組構(gòu)、巖心熱紅外掃描冷斑異常和孔隙水氯離子濃度低值異常等)及巖心—測井估算的水合物飽和度進(jìn)行對比分析的基礎(chǔ)之上。由于巖心和標(biāo)準(zhǔn)測井資料固有的缺陷(巖心不完整,易擾動變形;標(biāo)準(zhǔn)測井受井孔條件制約很難真實反映原位特征且分辨率不高),根據(jù)這些資料很難對穩(wěn)定帶內(nèi)水合物的非均質(zhì)分布特征及其巖性控制機(jī)理得以全面的理解。相比較而言,連續(xù)的、覆蓋整個水合物穩(wěn)定帶且?guī)缀醪皇芫鬃冃魏湍酀{污染影響的高分辨率隨鉆測井資料,在反演含水合物沉積物的巖性、預(yù)測水合物非均質(zhì)分布方面的應(yīng)用還很少。

        2)現(xiàn)有的應(yīng)用于水合物儲層測井評價的地層模型還過于簡單。以孔隙度和水合物飽和度計算為例,現(xiàn)有模型基本上僅考慮了兩種(骨架+孔隙水,水合物被視為骨架或孔隙流體的一部分)或三種地層組分(骨架+水合物+孔隙流體)。雖然Lee等[47]提出了所謂的四組分模型(骨架+黏土+水合物+孔隙流體),但因涉及參數(shù)過多,實際中很少獲得應(yīng)用。況且,即便是四組分模型也難以刻畫深海含水合物沉積物復(fù)雜的巖性和流體組分的全貌。導(dǎo)致上述局面的一個重要原因就是所使用的測井資料相對單一。比如阿爾奇公式計算水合物飽和度時基本上只用到電阻率和孔隙度兩種測井資料;各種形式的三相聲波方程也只用到了聲波速度和孔隙度測井資料。將水合物與其宿主沉積物視為一個整體,考慮更為復(fù)雜的、能夠反映含水合物沉積物真實地質(zhì)特點(diǎn)的地層模型,綜合利用多種高分辨率隨鉆測井?dāng)?shù)據(jù),聯(lián)合反演地層的巖性組分、孔隙度和水合物飽和度,或許是未來水合物儲層測井評價技術(shù)發(fā)展的一個重要方向。

        3)盡管目前大洋鉆探已積累了來自不同大陸邊緣共30個站位的測井資料,但對于不同大陸邊緣或同一大陸邊緣不同站位之間水合物穩(wěn)定帶內(nèi)的水合物分布特征還缺乏必要的橫向?qū)Ρ妊芯?。不同大陸邊緣之間或相同大陸邊緣不同站位之間水合物縱向分布存在著哪些共性和差異,目前并不清楚。海底水合物穩(wěn)定帶內(nèi)水合物的垂向非均質(zhì)分布的模式尚未建立起來,由此限制了我們對水合物分布和形成機(jī)理的深入了解。

        4)大洋鉆探水合物研究中,巖心和測井資料的綜合研究還很欠缺。地質(zhì)學(xué)家基于巖心的研究大多偏重于地質(zhì)樣品的分析、測試和鑒定;地球物理學(xué)家基于測井資料的分析往往側(cè)重于孔隙度和飽和度的定量參數(shù)的計算。水合物與宿主沉積物是一個整體,單純的巖心或測井研究均只能“管中窺豹略見一斑”,唯有測井資料特別是高分辨率隨鉆測井資料與巖心資料相結(jié)合的綜合研究,才有可能真正了解水合物—宿主沉積物系統(tǒng)的全貌。從這個角度來講,巖心—測井綜合研究亟待加強(qiáng)。

        劉堂宴教授閱讀本文初稿,王秀娟研究員、韓同城教授審閱本文,并提出了寶貴的修改意見。在此一并致謝!

        猜你喜歡
        航次水合物站位
        提高政治站位 對標(biāo)國內(nèi)一流
        建黨百年說“站位”
        水上消防(2021年3期)2021-08-21 03:12:00
        氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
        提升站位講政治 創(chuàng)新?lián)?dāng)爭出彩
        熱水吞吐開采水合物藏數(shù)值模擬研究
        天然氣水合物保壓轉(zhuǎn)移的壓力特性
        我國海域天然氣水合物試采成功
        我國集裝箱航運(yùn)企業(yè)實施作業(yè)成本管理法面臨的困難及解決方案
        集裝箱化(2014年10期)2014-10-31 18:26:46
        電梯站位暗示什么
        中文亚洲成a人片在线观看| 久久www免费人成精品| 97久人人做人人妻人人玩精品| 亚洲日韩精品一区二区三区| 国产午夜精品一区二区三区嫩草| 无码国产69精品久久久孕妇| 久久国产亚洲AV无码麻豆| 亚洲成a人v欧美综合天堂麻豆| 欧美成人免费高清视频| 国产欧美另类精品久久久| 免费观看视频在线播放| 午夜一区二区三区在线观看| 午夜少妇高潮在线观看视频| 极品美女一区二区三区免费| 色老板美国在线观看| 精品国内自产拍在线观看| 欧美成aⅴ人高清免费| 国产高潮流白浆免费观看不卡| 亚洲一区二区三区在线观看| 美女一区二区三区在线视频| 男女主共患难日久生情的古言| 日韩视频在线观看| 日韩国产成人无码av毛片蜜柚| 国产裸体歌舞一区二区| 久久青青草原亚洲AV无码麻豆| 中文人妻av大区中文不卡| 国产精品日本一区二区三区在线| 亚洲天堂精品一区入口| 丰满人妻被两个按摩师| 国产成人亚洲精品无码h在线| 福利一区二区三区视频午夜观看| 激,情四虎欧美视频图片| 国产丝袜美腿在线播放| 亚洲色偷偷综合亚洲avyp| 无码人妻久久一区二区三区app| vr成人片在线播放网站| 国产亚洲曝欧美不卡精品| 色和尚色视频在线看网站| 日韩网红少妇无码视频香港| 精品人妻系列无码一区二区三区| 亚洲av美女在线播放啊|