介觀孔隙中天然氣水合物生成過程模擬

摘 要

摘 要:目前,對沉積物中天然氣水合物(以下簡稱水合物)的生成過程研究大部分集中在實驗研究以及水合物生成的本征動力學(xué)方面,在孔隙尺度方面缺乏對水合物生成過程的全方位考慮。

摘 要:目前,對沉積物中天然氣水合物(以下簡稱水合物)的生成過程研究大部分集中在實驗研究以及水合物生成的本征動力學(xué)方面,在孔隙尺度方面缺乏對水合物生成過程的全方位考慮。為此,從介觀角度出發(fā),結(jié)合實驗?zāi)M結(jié)果,研究了沉積物孔隙中甲烷及二氧化碳水合物的生成過程。首先對甲烷水合物和二氧化碳水合物在沉積物中的生成過程進(jìn)行了實驗?zāi)M(沉積物樣品選用南海淺表層沉積物,粒徑為60100),結(jié)果表明:水合物在沉積物孔隙中的生成過程比較平緩,體系溫度基本沒有大的變化;在初始階段水合物生成量比較大,隨著反應(yīng)的進(jìn)行,水合物生成速率逐漸減小。在實驗的基礎(chǔ)上,以孔隙水中溶解氣體的濃度為變量,從介觀角度數(shù)值模擬了甲烷水合物和二氧化碳水合物在沉積物中的生成過程,并以單個沉積物孔隙空間的水合物生成表征沉積物體系內(nèi)水合物的生成特性。結(jié)論認(rèn)為:沉積物顆粒堆積孔隙內(nèi)部的水合物先在沉積物壁面處生成,然后水合物層逐漸加厚,向孔隙中心生長,水合物呈層狀生成,最后填滿整個沉積物顆??紫?;伴隨著水合物的生成,沉積物體系孔隙率降低。通過模擬計算得到的水合物轉(zhuǎn)化率與實驗結(jié)果進(jìn)行對比,其誤差范圍介于3%~7%,表明該模型具有較強(qiáng)的可靠性。

關(guān)鍵詞:天然氣水合物  孔隙  生成過程  模擬  介觀角度  甲烷水合物  二氧化碳水合物  轉(zhuǎn)化率  實驗室試驗

A simulation study of natural gas hydrate generation in mesoscopic-scale pores

AbstractThe studies of natural gas hydrate formation process in deposits are still at the experimental stage and mostly focus on that in terms of intrinsic kinetics,but rarely consider the scale of pores in the depositsIn view of this,we made an investigaffon,based on the simulation results,into the formation process of methane and carbon dioxide hydrate in mesoscopiC seale pores of depositsFirstwe conducted such an experimental simulation study,in which the shallow deposits in the South China Sea were chosen as the samples with the pore size of 60-100 meshesThe following results were achievedaThe process of natural gas hydrate formation in the pores of the deposits is relatively stable and slow and the temperature of the system basically keeps stable without significant fluctuation  bThe rate of gas hydrate formation is relatively great at the initial stage and then decreases progressively along with the progress of reactionOn this basis,with the concentration of the dissolved gas in pore water as the variablea mesoscopic pore space model was established to characterize the gas hydrate formation process in the pores of depositsThus,the other findings were obrainedaGas hydrate first is generated at the interface of a sediment particle,then becomes thickening and grows to the ccnter of pore space,and finally forms in layers until it fills up the whole pore spacebRelative porosity of sediments gradually decreases with the formation of gas hydrateIn comparison with the experimental datathe calculated resuhs from the simulation modeling produce an error of 37%,indicating that the simulation model is of rather high reliability

Key wordsnatural gas hydrate,pore,formation process,simulation,mesoscopic scale,methane hydrate,CO2 hydrate,conversion rate

近年來的資源調(diào)查、評價成果表明,中國南海北部天然氣水合物(以下簡稱水合物)賦存天然氣資源量約為185×108t油當(dāng)量[1-4],若能加以開采利用則能有效緩解中國能源需求與供給之間的矛盾。

自然界水合物發(fā)育區(qū)一般都位于深海海底以及大陸永久凍土帶中,水合物在沉積物多孔介質(zhì)中生成,其沉積物堆積體系孔隙尺寸比較大(通常大于50nm),介予大孔的范疇。水合物在沉積物孔隙里的生成是由許多不同的參數(shù)共同作用的,除了受海底溫度、壓力等熱力學(xué)條件的影響外,氣體組分以及孔隙水組成也控制著水合物的生成過程,其結(jié)晶和生長還取決于沉積物顆粒的大小、形狀和組成[5]ClarkeBishnoi[6]提出了一些水合物形成和分解的動力學(xué)模型,將水合過程視為鏈?zhǔn)椒磻?yīng)來研究其生長動力學(xué)。BuffettZatsepina[7-8]較真實地模擬了海洋環(huán)境中天然氣水合物的生成過程,實驗發(fā)現(xiàn)在沒有游離氣體的多孔介質(zhì)中,氣體以溶解態(tài)擴(kuò)散作用進(jìn)入孔隙流體生成水合物,并由此建立了溶解氣體生成水合物模型。Maddena[9]研究了在不同直徑沉積物顆粒體系中的水合物生成過程,實驗結(jié)果說明沉積物粒徑對水合物生成有一定的影響。Makogon[10]較早模擬了不同粒級砂巖中水合物形成的過程,發(fā)現(xiàn)若溫度不變,天然氣水合物在孔徑為3.55.7mm的孔隙中,其形成壓力隨孔徑的減小而上升,而在更大的孔隙中其形成和分解條件完全與塊層狀水合物相同。Prasadu婦利用粒徑為50mm150mm200mm的沉積物開展實驗,結(jié)果表明粒徑對沉積物中水合物分解并沒有特別的影響,同時,較之于純水體系來說,沉積物更能促進(jìn)水合物結(jié)晶和生長。

但是,沉積物中水合物的生成過程仍舊是一個未解的難題?,F(xiàn)有的研究大部分都集中在實驗研究以及水合物生成的本征動力學(xué)方面,在孔隙尺度方面缺乏對沉積物中水合物生成過程的全方位考慮。同時,在中國南海北部水合物發(fā)育區(qū)內(nèi)發(fā)現(xiàn)不僅有甲烷水合物的存在,也有大量非烴氣體水合物如二氧化碳水合物的賦存。因此,筆者從介觀角度出發(fā),結(jié)合實驗?zāi)M結(jié)果,研究了沉積物孔隙中甲烷及二氧化碳水合物的生成過程,以期為水合物勘探開發(fā)提供理論基礎(chǔ)和支撐。

1 實驗?zāi)M

在進(jìn)行數(shù)值模擬之前,對甲烷和二氧化碳?xì)怏w水合物的生成過程先進(jìn)行了實驗室模擬,實驗?zāi)M裝置主要由管路、氣瓶、增壓裝置、真空泵、恒溫水浴和反應(yīng)釜組成,具體見本文參考文獻(xiàn)[12]。沉積物樣品選用南海淺表層沉積物,人工進(jìn)行清洗并干燥,后進(jìn)行篩分,按照一定的粒徑進(jìn)行分類。本次實驗粒徑分布介于60100目的沉積物,測得其孔隙度為36.7%。甲烷水合物和二氧化碳水合物生成過程分別在不同的溫度和壓力條件下進(jìn)行,觀測水合物生成過程中體系內(nèi)溫度和壓力的變化。實驗結(jié)果如圖1所示。

 

從甲烷水合物和二氧化碳水合物生成過程的壓力變化規(guī)律可以看出(1),水合物生成過程中,體系溫度基本沒有大的變化。初始階段反應(yīng)釜內(nèi)壓力下降比較明顯,這表明水合物在短時間內(nèi)生成,消耗了大量的孔隙水溶解氣。同時由于水合物容易在氣相和液相的交界面處形成,界而處的氣體也被大量消耗,這就導(dǎo)致氣相與沉積物體系之間存在比較明顯的氣體濃度差,氣體擴(kuò)散進(jìn)入沉積物體系中的驅(qū)動力增大,氣體消耗量也就比較大。因此,氣相壓力下降比較迅速。當(dāng)反應(yīng)進(jìn)行到一定階段后,壓力下降速度變緩,此時,水合物生成量比較大,氣體擴(kuò)散阻力增加,氣相中的氣體若想進(jìn)入水合物籠中需要克服較大的擴(kuò)散阻力,反應(yīng)速度下降,壓力變化也就隨之變小。

2 模型處理

本次的研究對象為沉積物堆積體系內(nèi)的單個孔隙空間,屬于多孔介質(zhì)內(nèi)的孔隙結(jié)構(gòu),幾何邊界比較復(fù)雜,計算區(qū)域比較小,而格子玻爾茲曼方法(LBM)在處理復(fù)雜幾何邊界的問題中具有優(yōu)勢[13-14]。因此采用LBM方法建立二維D2Q9模型,表征孔隙空間內(nèi)水合物生成過程中各參數(shù)的變化。

由于沉積物孔隙內(nèi)水合物實際生成過程比較復(fù)雜,因而模型對物理過程進(jìn)行了簡化,假定所有沉積物顆粒都為正網(wǎng)形且粒徑分布均一,沉積物體系內(nèi)各向同性。自然界中,沉積物層位于海底,由于海水和沉積層的壓力,沉積物顆粒被壓實,孔隙水充填于沉積物孔隙內(nèi)部,氣體溶解于填充的孔隙水中。沉積物顆粒間彼此錯落排列,形成沉積物間孔隙的是顆粒的圓弧面,計算區(qū)域的邊界為沉積物顆粒的弦長,如圖2所示。邊界條件采用反彈邊界。

 

根據(jù)Kang[15-20]的模型,結(jié)合實際沉積物顆粒的條件以及水合物反應(yīng)的特殊性質(zhì),選用D2Q9模型,分別對沉積物孔隙水和孔隙水中的溶解氣的濃度參數(shù)建立不同的控制方程,耦合求解。

3 模擬結(jié)果

31 沉積物孔隙內(nèi)濃度分布規(guī)律

實際計算區(qū)域僅為沉積物顆粒之間的縫隙,即4顆沉積物顆粒中間形成的區(qū)域,如圖3、4中所示為該區(qū)域的二維切片,顯示的是區(qū)域中間溶解氣體的濃度場。

 

 

 

沉積物孔隙面相對比較平整,水合物先在沉積物壁面生成,而后沿著壁面法線方向增長,水合物在生成過程中呈層狀生成,最后在孔隙中間會合,從而充滿整個沉積物孔隙,如圖3所示。甲烷水合物層的表面(即沉積物壁面)是光滑的,無液體滲漏進(jìn)沉積物表面,當(dāng)水合物形成后,沉積物顆粒表面即被水合物層所覆蓋。在水合物層表面,存在一個濃度梯度,氣體不斷通過水合物表面擴(kuò)散進(jìn)入水合物,存在一個穩(wěn)定的質(zhì)量傳遞過程,直至反應(yīng)結(jié)束。

當(dāng)采用CO2氣體作為水合物生成氣時,水合物生成反應(yīng)時間相對甲烷來說比較短,但生成過程與甲烷水合物生成過程比較相似,如圖4所示。水合物先在沉積物壁面和壁面的交界處生成,而后沉積物在初始生成處聚集,并呈層狀向外擴(kuò)展,往沉積物孔隙中間匯集。因此,沉積物孔隙中水合物的生成過程呈層狀生成,最后在孔隙中間會合,從而充滿整個沉積物孔隙。

32 水合物轉(zhuǎn)化率的變化

在模型假設(shè)中,假定沉積物顆??讖绞蔷坏模覟闃?biāo)準(zhǔn)的正圓形,孔隙水均勻地分布于沉積物孔隙中,整個沉積物和孔隙水體系為各向均勻體系。整個體系內(nèi)可以看作是由同樣的沉積物顆粒和顆粒組成的孔隙所組成,而水合物的生成區(qū)域也即孔隙水占據(jù)的空間就是由一個一個的孔隙組成,孔隙之間具有相似性。因此,可以用一個孔隙內(nèi)的情況代表整個體系內(nèi)的情況。單個孔隙內(nèi)水合物的轉(zhuǎn)化率代表整個體系內(nèi)所有孔隙內(nèi)水合物的轉(zhuǎn)化率。將不同時刻消耗的濃度進(jìn)行積分,得到整個反應(yīng)區(qū)間內(nèi)所消耗的氣體量,也就可以計算得出沉積物孔隙內(nèi)孔隙水轉(zhuǎn)化為水合物的比率。

當(dāng)水合物反應(yīng)氣體分別為甲烷和二氧化碳時,水合物反應(yīng)的水合數(shù)為5.75,表達(dá)式為

 

因此水合物轉(zhuǎn)化率計算公式為

 

式中hh為水合物轉(zhuǎn)化率;ng為氣體消耗摩爾量;nw為沉積物孔隙水的摩爾量。

為了驗證模擬結(jié)果的可靠性,將模擬結(jié)果與實驗結(jié)果進(jìn)行對比,如圖5、6所示。實驗得到的甲烷水合物和二氧化碳水合物的轉(zhuǎn)化率均大于模擬結(jié)果。這主要由以下幾個因素造成。

 

 

首先,實驗采用的沉積物并非是孔徑為150mm的均一顆粒,而是介于l50250mm的區(qū)間,因此,形成的顆粒間的孔隙尺寸也就不統(tǒng)一,孔隙尺寸分布呈一定區(qū)間,對孔隙空間內(nèi)水合物生成過程的影響也就不盡相同,而模型模擬的是沉積物粒徑為l50mm的孔隙空間,孔隙尺寸是恒定的。其次,實際實驗過程中,氣體與沉積物體系的交界面處最容易生成水合物,同時,界面處的水合物容易通過毛細(xì)作用將沉積物體系內(nèi)孔隙水運移到水合物周圍,而水合物生成區(qū)域的擴(kuò)散傳遞作用是最強(qiáng)烈的。氣體的擴(kuò)散時刻都在發(fā)生,可能直接擴(kuò)散進(jìn)入沉積物孔隙水里,也可能直接擴(kuò)散進(jìn)入水合物中,生成水合物,而模型的假設(shè)是氣體必須首先擴(kuò)散進(jìn)入沉積物孔隙水中,再從孔隙水中通過水合物膜表面擴(kuò)散生成水合物。因此,綜合考慮,模擬結(jié)果與實驗結(jié)果存在一定的偏差,但誤差比較小,誤差范圍為3%~7%,模型結(jié)果基本可靠。

4 結(jié)論

1)針對甲烷水合物和_二氧化碳水合物在沉積物中的生成過程,分別進(jìn)行了實驗?zāi)M和介艦數(shù)值模擬,并利用實驗?zāi)M結(jié)果驗證了數(shù)值模擬的可靠性。實驗?zāi)M結(jié)果表明,水合物在沉積物孔隙中的生成過程比較平緩,體系溫度基本沒有大的變化且在初始階段水合物生成量比較大,隨著反應(yīng)的進(jìn)行,水合物生成速率逐漸減小。

2)介觀尺度模型主要研究的是沉積物顆粒組成的孔隙,以孔隙水中溶解氣體的濃度為變量,模擬了由沉積物顆粒組成的單個孔隙內(nèi)水合物的生成過程。結(jié)果表明,在沉積物顆粒體系內(nèi)部,孔隙面相對比較平整,水合物先在沉積物壁面生成,而后沿著壁面法線方向增長,水合物在生成過程中呈層狀生成,最后在孔隙中間會合,從而充滿整個沉積物孔隙。

3)單個沉積物孔隙中水合物的轉(zhuǎn)化率結(jié)果計算值與實驗值的對比說明模擬結(jié)果與實驗結(jié)果存在一定的偏差,誤差范圍為3%~7%,模型具有一定的可靠性。

 

參考文獻(xiàn)

[1]WU Nengyou,ZHANG Haiqi,SU Kinet alHigh concentrations of hydrate in disseminated forms found in very fine grained sediments of Shenhu area,South China Sea[J]Terra Nostra2007,(12)236-237

[2]趙宗彬,仇性啟,許俊良.深海天然氣水合物鉆探取心鉆柱振動模態(tài)分析[J].天然氣工業(yè),2011,31(1)73-76

ZHAO Zongbin,QIU XingqiXU JunliangFinite element analysis of modal vibrations of drill strings in deep-sea natural gas hydrate coringdrilling[J]Natural Gas Industry,2011,31(1)73-76

[3]王智鋒,管志川,許俊良.天然氣水合物深水深孔鉆探取心系統(tǒng)研制[J].天然氣工業(yè),2012,32(5)46-48

WANG ZhifengGUAN Zhichuan,XU JunliangResearch and development of a gas hydrate coring system in the deep sea water[J]Natural Gas Industry2012,32(5)46-48

[4]傅寧,林青,劉英麗,等.白云凹陷天然氣水合物氣源條件再認(rèn)識[J].中國石油勘探,201116(4)48-53

FU NingLIN Qing,LIU Yingli,et alA Rethinking upon gas source condition of gas hydrate in Baiyun Sag[J]China Petroleum Exploration,2011,16(4)48-53

[5]CLENNELL M B,HOVLAND M,BOOTH J Set alFormation of natural gas hydrates in marine sediments1Conceptual model of gas hydrate growth conditioned by host sediment properties[J]Journal of Geophysical Research,1999104(B10)22985-23003

[6]CLARKE M A,BISHNOI P RDetermination of the intrinsic kinetics of CO2 gas hydrate formation using in situ particle size analysis[J]Chemical Engineering Science,2005,60(3)695-709

[7]ZATSEPINA O YBUFFETT B APhase equilibriunl of gas hydrateImplications for the formation of hydrate in the deep sea floor[J]Geophysical Research Letters,l997,24(13)1567-1570

[8]ZATSEPINA O YBUFFETT B ANucleation of gas hydrate in marine environments[J]Geophysical Research Letters,2003,30(9)1451

[9]MADDENA M E,ULRICHB S,SZYMCEKB Pet alExperimental formation of massive hydrate deposits from accumulation of CH4 gas bubbles within synthetic and natural sediments[J]Marine and Petroleum Geology,2009,26(3)369-378

[10]MAKOGON Y FHydrate of natural gas[M]TulsaPenn Well1981

[11]PRASAD P S R,CHARI V D,SHARMA D Vet alEffect of silica particles on the stability of methane hydrates[J]Fluid Phase Equilibria,2012318110-114

[12]ZANG Xiaoya,IAANG Deqing,WU NengyouGas hydrate formation in fine sand[J]Science ChinaEarth Seiences2013,56(4)549-556

[13]郭照立,鄭楚光.格子Boltzmann方法的原理及應(yīng)用[M].北京:科學(xué)出版社,2008

GUO Zhaoli,ZHENG ChuguangTheroy and applications of lattice Boltzmann method[M]BeijingScience Press,2008

[14]何雅玲,王勇,李慶.格子Boltzmann方法的理論及應(yīng)用[M].北京:科學(xué)出版社,2008

HE Yaling,WANG YongLI QingLattice Boltzmann methodTheory and applications[M]BeijingScience Press,2008

[15]JEONG NCHOI D H,LIN C LEstimation of thermal and mass diffusivity in a porous medium of complex structure using a lattice Boltzmann method[J]International Journal of Heat and Mass Transfer2008,51(1516)3913-3923

[16]KANG Qinjun,ZHANG Dongxiao,LICHTNER P Cet alLattice Boltzmann model for crystal growth from supersaturated solution[J]Geophysical Research Iretters,200431(21)604

[17]KANG Qinjun,TSIMPANOGIANNIS I NZHANG Dongxiao,et alNumerical modeling of pore-scale phenomena during CO2 sequestration in oceanic sediments[J]Fuel Processing Technology2005,86(1415)1647-1665

[18]KANG Qinjun,ZHANG Dongxiao,CHEN Shiyi,et alLattice Boltzmann simulation of chemical dissolution in porous media[J]Physical Review E2002,65(3)6318-6325

[19]KANG Qinjun,ZHANG DongxiaoCHEN ShiyiSimulation of dissolution and precipitation in porous media[J]Journal of Geophysical Research,2003,108(B10)2505

[20]LU G PDEPAOLO D J,KANG Qinjun,et alLattice Bohzmann simulation of snow crystal growth in clouds[J]Journal of Geophysical ResearchAtmospheres,2009,114(D7)

 

本文作者:梁德青  臧小亞  吳能友

作者單位:中國科學(xué)院廣州能源研究所

  中國科學(xué)院可再生能源與天然氣水合物重點實驗室