天然氣水合物生成條件預(yù)測模型及應(yīng)用畢業(yè)設(shè)計

1、天然氣水合物生成條件預(yù)測模型及應(yīng)用1引言1.1研究的意義和目的水合物是目前科學(xué)領(lǐng)域中的熱門課題，不僅與石油天然氣開采、儲存和運輸密切相關(guān)，而且與環(huán)境保護(hù)、氣候變遷，特別是人類未來賴以生存的能源有關(guān)。一方面，由于水合物特殊的物理化學(xué)性質(zhì)潛藏著巨大的工業(yè)應(yīng)用潛能，另一方面，在天然氣開采及儲運過程中，其壓力和溫度會發(fā)生變化，尤其當(dāng)天然氣通過突縮孔眼時，由于摩擦耗能會產(chǎn)生較大的節(jié)流壓降和溫降。同時在天然氣管道中，天然氣的飽和含水量隨天然氣壓力增加或溫度降低而降低，在天然氣經(jīng)過變徑管、閥門或過濾器、油嘴等節(jié)流時，雖然壓力下降，但溫度也下降，管道中天然氣的飽和含水量隨天然氣壓力增加或溫度降低而降低，另外由

2、于節(jié)流產(chǎn)生渦流，就有可能達(dá)到水合物生成的條件。水合物會導(dǎo)致閥門堵塞、氣井停產(chǎn)、管道停輸?shù)葒?yán)重事故，所以研究預(yù)測水合物的生成溫度和壓力條件以及如何防止水合物的生成是重要以及必要的。通常，在輸送天然氣過程中清除水合物的方法是用熱水或熱蒸汽對管道進(jìn)行加熱，在水合物和金屬接觸點上，將溫度提高到3040，使水合物很快分解。據(jù)統(tǒng)計防止水合物生成的費用約占生產(chǎn)總成本的58%。在工程上對抑制劑用量不能準(zhǔn)確計算，抑制劑的用量往往大于實際需求量，這樣一方面不利于節(jié)約成本，另一方面導(dǎo)致不必要的環(huán)境污染針對上述問題，需要用科學(xué)的實驗方法，準(zhǔn)確測定天然氣水合物的生成條件，并篩選和評價抑制劑的抑制效果，從而為天然氣集輸管

3、道水合物防治工作提供科學(xué)依據(jù)。同時，針對旅大10-1天然氣壓力高、流速高以及組分復(fù)雜的特殊情況，需要建立正確的天然氣水合物生成預(yù)測模型及加入抑制劑后水合物的生成預(yù)測模型，并開發(fā)出相應(yīng)的實用程序。1.2國內(nèi)外發(fā)展概況及分析國內(nèi)外預(yù)測水合物形成壓力和溫度的方法大致可分為經(jīng)驗公式法、圖解法、相平衡計算法和統(tǒng)計熱力學(xué)模型法四大類。經(jīng)驗公式法是通過室外取樣室內(nèi)實驗或現(xiàn)場測試水合物形成壓力和溫度,在此基礎(chǔ)上,擬合得到的半經(jīng)驗或經(jīng)驗?zāi)Ｐ?。?jīng)驗、半經(jīng)驗?zāi)Ｐ途哂泻唵?、計算方便等?yōu)點,但使用范圍狹窄,計算精度不高。經(jīng)驗圖解法是根據(jù)甲烷及不同相對密度天然氣形成水合物的平衡曲線，大致確定出天然氣形成水合物的溫度和壓力

4、。但對含有硫化氫的天然氣誤差較大，若相對密度在兩條曲線之間,需采用內(nèi)插法進(jìn)行近似計算。為了便于計算機(jī)運算，圖解法被回歸成幾個公式，若已知天然氣的相對密度和溫度,可選擇其中合適的公式計算水合物形成壓力。若已知相對密度和壓力可選擇其中合適的公式進(jìn)行迭代求得水合物形成溫度。但是這種方法也存在適用范圍狹窄，計算精度不高的缺點。相平衡計算法是1940年katz根據(jù)氣固平衡常數(shù),提出了一種估算天然氣水合物生成條件的方法,可用于計算含有典型烷烴組成的無硫天然氣,而對非烴含量多的氣體及壓力高于6.9mpa時誤差較大。水合物的相平衡理論模型大多是在van der waals-platteeuw模型的基礎(chǔ)上發(fā)展起

5、來的,具有很強的理論基礎(chǔ),計算精度高,但計算較復(fù)雜,難掌握。目前預(yù)測水合物生成條件的熱力學(xué)模型幾乎都是以經(jīng)典統(tǒng)計熱力學(xué)為基礎(chǔ)的。vander waals和platteeuw 在1959年提出的 langmuir氣體吸附模型，沒有考慮到客體分子間的相互作用，而假定洞穴為球形對稱等等。使用此模型，1964年，saito等人提出了一種預(yù)測水合物生成條件的方法。1972年，parrish和pransnitz考慮到langmuir參數(shù)僅為溫度的函數(shù)，提出了簡單的經(jīng)驗關(guān)系式計算langmuir常數(shù)，大大簡化了vander waals模型，并首次將vander waals模型推廣到多元體系。1977年ng&

6、robinson模型(簡稱nr)，考慮不同客體間的相互作用，修正van der waals模型。1983年holderjohn模型（簡稱hj）考慮到洞穴并非球形，水分子與洞穴中心距離不等的實際情況，采用三層球模型來描述水分子和客體分子間的相互作用，并引人擾動因子來矯正球形分子的langmuir常數(shù)；1988年杜亞和和郭天民（簡稱du&guo88模型）利用十點gauss數(shù)值積分法計算出不同溫度下langmuir常數(shù)；1996年，chen- guo基于水合物生成動力學(xué)，提出了一個完全不同于vdw- p模型的新水合物預(yù)測模型。他們認(rèn)為，在生成水合物時，體系存在準(zhǔn)化學(xué)平衡和物理吸附平衡。左有祥、郭天民

7、基于水合物兩步生成動力學(xué)機(jī)理提出了新的水合物熱力學(xué)模型，與傳統(tǒng)vdw-p模型相比，表達(dá)式得到簡化且更接近普通的溶液熱力學(xué)模型，該預(yù)測模型都己成功應(yīng)用于純水體系、含醇或含鹽體系。1.3本文的主要工作本文研究的主要內(nèi)容是天然氣水合物生成條件預(yù)測模型及應(yīng)用。綜合考慮了天然氣管流及嘴流的壓力和溫度等水合物生成的環(huán)境因素、抑制劑及其濃度對水合物生成條件的影響，分別介紹了基于熱力學(xué)相平衡理論的水合物生成條件預(yù)測模型和基于統(tǒng)計力學(xué)的水合物生成條件預(yù)測模型并建立了適合含醇類抑制劑體系的天然氣水合物預(yù)測綜合數(shù)學(xué)模型。本文還建立了抑制劑用量計算模型并討論了管道伴熱的一般辦法，為論證和優(yōu)選合理的水合物生成防止工藝措

8、施提供了重要的技術(shù)依據(jù)。本文主要完成了以下工作：1) 基于質(zhì)量、動量、能量守恒原理，導(dǎo)出了天然氣管流壓力、溫度預(yù)測模型；2) 在對天然氣節(jié)流機(jī)理研究的基礎(chǔ)上，考慮焦耳-湯姆遜系數(shù)對天然氣節(jié)流壓降、溫降的影響，綜合討論確定了焦耳-湯姆遜系數(shù)的bwrs方程計算方法；3) 對天然氣水合物的生成條件及機(jī)理進(jìn)行了深入研究，總結(jié)了目前常用的天然氣水合物生成條件預(yù)測方法，著重比較了應(yīng)用統(tǒng)計熱力學(xué)的理論基礎(chǔ)的幾種水合物生成條件預(yù)測模型，為本文應(yīng)用以chen&guo96模型作為基礎(chǔ)的改進(jìn)模型提供理論依據(jù)。4) 應(yīng)用基于統(tǒng)計力學(xué)的理論建立天然氣水合物生成條件預(yù)測模型，并結(jié)合天然氣管流、節(jié)流壓力溫度預(yù)測模型建立管道

9、某處水合物形成預(yù)測模型。5) 考慮了醇類抑制劑加入對體系相平衡的影響，改進(jìn)了pr狀態(tài)方程并推導(dǎo)了水活度系數(shù)的計算方程，使其適用于含醇類抑制劑體系水合物生成條件的預(yù)側(cè)；對于電解質(zhì)抑制劑的體系，考慮用左有祥等改進(jìn)的pt狀態(tài)方程計算電解質(zhì)溶液中水的分逸度，推導(dǎo)出電解質(zhì)體系天然氣水合物生成條件預(yù)測的模型。6）對天然氣水合物的防止措施進(jìn)行了總結(jié)和討論，并且針對不同的現(xiàn)場情況提出了不同的水合物防止措施以及解堵措施。討論了天然氣加熱伴熱的幾種方式，并且進(jìn)行了天然氣最大容許水含量計算，為加熱、脫水或注抑制劑防止水合物形成提供了技術(shù)支持；7) 基于hammerschmidt方程，結(jié)合ng和chen提出的用溶解度

10、系數(shù)描述醇類損失量的方法，導(dǎo)出了可較準(zhǔn)確計算醇類抑制劑用量的算法，同時對電解質(zhì)抑制劑加入量也進(jìn)行了計算；提出了基于水合物兩步生成機(jī)理的醇類抑制劑和電解質(zhì)抑制劑注入量的嚴(yán)格算法；8）利用vi.suai basic語言將上述理淪方法軟件化，編制了有關(guān)天然氣管流及節(jié)流的壓力、溫度和天然氣水合物生成條件預(yù)測的系統(tǒng)軟件。2天然氣管流和節(jié)流壓力溫度預(yù)測 正確的預(yù)測天然氣在管道中流動的壓力溫度分布，不僅有利于掌握天然氣流動過程的動態(tài)信息，也為天然氣的安全生產(chǎn)和輸送提供了技術(shù)支持。根據(jù)流體力學(xué)的知識，建立氣體流動相應(yīng)的連續(xù)性方程、運動方程和能量方程，并且在此基礎(chǔ)上導(dǎo)出并求解了描述天然氣管流壓力、溫度、流速及密

11、度分布的數(shù)值模型。并考慮焦耳-湯姆遜系數(shù)對節(jié)流溫降的影響，建立的節(jié)流溫降壓降的求解方法。 2.1天然氣管流壓力溫度計算2.1.1基本方程天然氣在直管道中流動可視為一元穩(wěn)定流動，在管路中取截面積為a、長為，流向為的管柱，如圖2. 1所示。質(zhì)量、動量和能量守恒方程如下：質(zhì)量守恒方程若不考慮管道截面的變化，由質(zhì)量守恒定律，氣體流動連續(xù)性方程為 式中 氣體的密度，kg/m3； 氣體的流速，m/s； 時間變量，s； 沿管長變量，m。對于穩(wěn)定流動，其流動參數(shù)不隨時間而變化，式中左邊第一項為0。因此，連續(xù)性方程為 即質(zhì)量守恒方程 動量守恒方程根據(jù)牛頓第二定律，由流體力學(xué)所建立的運動方程形式可寫為 式中 重力

12、加速度，m/s2；管道與水平面間的傾角，rad；水力摩阻系數(shù)；管道內(nèi)徑，m；管道中氣體壓力，pa。對于穩(wěn)定流動，方程中左邊第一項為0，并考慮，得動量守恒方程 能量方程根據(jù)能量守恒定律，由流體力學(xué)建立的能量方程為 式中 單位質(zhì)量氣體向外界放出的熱量，j/kg； 氣體內(nèi)能，j/kg； 氣體的焓，j/kg； 管道位置高度，m。對于穩(wěn)定流動，方程右邊第一項為0，同時考慮方程，則可化為 表示單位質(zhì)量氣體在管長為上的熱量損失，由傳熱學(xué)關(guān)系有 式中 管道的總傳熱系數(shù)，w/(m2k)； 管道埋深處地溫，k； 氣體質(zhì)量流量，kg/s。則可寫為 為了方便使,式可寫為 2.1.2比焓梯度由熱力學(xué)一般關(guān)系，焓是溫度和

13、壓力的狀態(tài)函數(shù)，其微分形式如下式 則比焓梯度可表示為 由于氣體在直管道中流動，其管徑變化一般不大，故焦湯系數(shù)甚小可忽略。所以，可認(rèn)為式中 氣體的比定壓熱容，j/(kgk) 。 比焓梯度可寫為 2.1.3水力摩阻系數(shù)目前水力摩阻系數(shù)的計算公式有無數(shù)種，為了較全面考慮摩阻因素，對于輸氣管道而言，流動一般處于紊流區(qū)，包括光滑管紊流(水力光滑區(qū))、部分粗糙管紊流(混合摩擦區(qū))和完全粗糙管紊流(阻力平方區(qū))三種情況。對于距離較短，沿線工藝簡單的輸氣管道，摩阻系數(shù)采用適用范圍為5000107和0.000040.05的條件下jain公式計算： 式中 雷諾數(shù)，無因次； 輸氣管道的內(nèi)表面粗糙度，mm； 輸氣管道

14、內(nèi)徑，式中單位為mm。2.1.4氣體總傳熱系數(shù)總傳熱系數(shù)表示氣體至周圍介質(zhì)的散熱強弱，由下式表示 式中 氣體至管內(nèi)壁的放熱系數(shù)，w/(m2k)； 管道外壁至周圍介質(zhì)的放熱系數(shù)，w/(m2k)； 第層（管壁、防護(hù)層、絕熱層）導(dǎo)熱系數(shù)，w/(m2k)； 管道內(nèi)徑，m； 官道上第層（管壁、防護(hù)層、絕熱層）的外徑，m； 確定總傳熱系數(shù)的計算管徑，m。當(dāng)時，取外徑；當(dāng)時，取平均值，即內(nèi)外直徑之和的一半；若，取內(nèi)徑。2.1.5數(shù)學(xué)模型及其求解將比焓梯度方程代入能量方程，則上述質(zhì)量、動量和能量守恒方程式、和中含氣體的壓力、溫度、流速及密度四個未知數(shù)，補充實際氣體的狀態(tài)方程 。上述基本方程可進(jìn)一步表示為壓力、

15、溫度、流速和密度的梯度方程組，根據(jù)管道入口處壓力和溫度，由氣體狀態(tài)方程可求出相應(yīng)的氣體密度和流速，并以此作為方程組的邊界條件式。 邊界條件式 式中 氣體偏差系數(shù)，無因次量； 氣體分子量，kg/mol； 分別將欲求的四個未知量、和記為（=1，2，3，4），相應(yīng)的梯度方程的右函數(shù)記為，則上述方程組式可表示為 （=1，2，3，4） 對于這類常微分方程組的初值問題可應(yīng)用計算精度較高的四階龍格庫塔法進(jìn)行數(shù)值求解。 將起點位置的函數(shù)值記為，取步長為，節(jié)點處的解可用四階龍格庫塔法表示為 （=1，2，3，4） 式中若未達(dá)到預(yù)計距離，再將節(jié)點的計算值作為下步計算的起點值，重復(fù)上述步驟，如此連續(xù)向前推算直達(dá)到預(yù)計

16、終點。2.2天然氣節(jié)流壓降溫降計算天然氣通過突縮孔眼時，由于摩擦耗能會產(chǎn)生較大的節(jié)流壓降和溫降，可能生成水合物。正確地預(yù)測天然氣節(jié)流后的溫度，可為預(yù)測水合物生成條件提供技術(shù)依據(jù),避免盲目性。2.2.1天然氣管道節(jié)流壓降計算天然氣流過節(jié)流裝置（如閥門）時，因過流截面突縮，其流速會迅速增大，造成局部阻力增大，使其壓力顯著下降。節(jié)流前后的壓力與其流量的關(guān)系為： 式中 氣體標(biāo)況下體積流量，m3/d； ，天然氣節(jié)流前、后的壓力，pa；節(jié)流孔眼直徑，mm；節(jié)流前的溫度，k；節(jié)流前的氣體偏差系數(shù)；天然氣相對密度；氣體絕熱指數(shù)。在已知天然氣的流量、節(jié)流前的溫度和壓力時,由上式可確定節(jié)流后的壓力。2.2.2天然

17、氣管道節(jié)流溫降計算天然氣通過地面節(jié)流孔眼時流速很高，在孔眼附近的氣流與外界的熱交換甚小，一般可忽略不計。若節(jié)流前后兩截面取在距離孔眼稍遠(yuǎn)的位置，則節(jié)流動能較其焓值甚小，也可忽略。故天然氣節(jié)流過程的穩(wěn)定流動能量方程可簡化為： 絕熱節(jié)流過程中的溫度效應(yīng)與氣體性質(zhì)有關(guān)，也與氣體所處狀態(tài)有關(guān)。絕熱節(jié)流效應(yīng)（又稱焦耳-湯姆遜效應(yīng)）常用溫度下降的數(shù)值與壓力下降數(shù)值的比表示，此時要引入焦耳-湯姆遜系數(shù)，即： 式中 焦耳-湯姆遜系數(shù)，k/pa； 氣體的熱力學(xué)溫度，k； 氣體的壓力，pa； 氣體比焓，j/kmol； 氣體的比容，m3/kmol； 氣體定壓摩爾比熱容，j/(kmolk)。節(jié)流前后氣體的溫度降的計算

18、式為： 式中 節(jié)流前后氣體的溫度降，k。焦耳-湯姆遜系數(shù)計算方法很多文獻(xiàn)中推薦所用的的計算方法通常是將的求值裝換成計算（為氣體的壓縮因子）或和（為氣體密度，單位是kg/m3），這樣使得計算量比較大，計算方法也較復(fù)雜，特別是應(yīng)用來計算，如果不引用的關(guān)聯(lián)式進(jìn)行修正，結(jié)果精度相對會較低。對于實際氣體，通常采用范德瓦爾方程（vdw）對其進(jìn)行描述，即： 其中 式中 摩爾氣體常數(shù)，8.314j/(molk)； 氣體的摩爾體積，m3/mol； 臨界溫度，k； 臨界壓力，pa。將代入得 范德瓦爾方程雖是描述實際氣體的狀態(tài)方程，但是大量的實踐證明其在定量上還不夠準(zhǔn)確，故需要根據(jù)已知的、求得，再代入計算。為了提高

19、計算精度，的求解也可由vdw方程的修正式推出。vdw方程的修正式主要有rk(redlich-kwong)方程、spk(soave-redlich-kwong)方程、pr(peng-robinson)方程、bwrs(benedict-webb-rubin-starling)方程等。bwrs狀態(tài)方程是公認(rèn)最精確的經(jīng)驗型狀態(tài)方程，適用性廣泛，計算精度高。為此使用bwrs方程建立了焦湯系數(shù)計算模型?？紤]到密度與比體積互成倒數(shù)的關(guān)系，式可進(jìn)一步寫為： 從式可以看出，只要能求出，就能夠求出焦耳-湯姆遜系數(shù)。這里采用bwrs方程計算，方程形式如下： 式中 系統(tǒng)壓力，kpa； 系統(tǒng)溫度，k； 汽相或液相的密度

20、，kmol/m3； 氣體常數(shù)，=8.3143kj/(kmolk)。、狀態(tài)方程的11個參數(shù)。式兩端在保持不變的條件下，對求導(dǎo)數(shù)，可得到： 式中，、和為符號參數(shù)，由下式計算: 將代入中，可以得到： 天然氣的定壓摩爾熱容與其組成、壓力、溫度有關(guān)，可近似按下式計算： 這個公式中定壓摩爾熱容的單位是kj/（mol）；對于混合氣體，計算定壓摩爾熱容時要用到當(dāng)量摩爾質(zhì)量： 式中 當(dāng)量摩爾質(zhì)量，kg/kmol； 混合氣體的組分?jǐn)?shù)； 組分的摩爾分?jǐn)?shù)，無因此小數(shù)； 組分當(dāng)量摩爾質(zhì)量，kg/kmol。已知混合氣體各組分的摩爾分?jǐn)?shù)以及節(jié)流前所處的壓力和溫度，可以從一些文獻(xiàn)和相關(guān)手冊中查出所給氣體的臨界參數(shù)（組分氣體臨

21、界溫度、組分氣體臨界壓力、組分氣體臨界密度和組分氣體當(dāng)量摩爾質(zhì)量）以及組分相互作用系數(shù)，可計算 bwrs方程的、這11個參數(shù)，方法如下：對于純組分的這11個參數(shù)和臨界參數(shù)，及偏心因子關(guān)聯(lián)如下 式中 、通用常數(shù)（=1，2，11）。starling通過對正構(gòu)烴烷的多性質(zhì)分析（-、焓和蒸氣壓）得到了和值，見表2.1。根據(jù)旅大10-1注氣系統(tǒng)天然氣組分，關(guān)聯(lián)時對正構(gòu)烷烴（甲烷至己烷）所用的臨界溫度、臨界密度及偏心因子見表2.2。表2.1 bwrs方程中的通用常數(shù)和值10.4436900.11544921.284380-0.92073130.3563061.7087140.544979-0.270896

22、50.5286290.34926160.4840110.75413070.0705233-0.04444880.5040871.3224590.03074520.179433100.07328280.463492110.006450-0.022143表2.2 使用bwrs方程時推薦使用的臨界性質(zhì)等數(shù)值名稱分子式分子量臨界溫度()臨界密度(kmol/ m3)偏心因子甲烷ch416.042-82.46110.0500.013乙烷c2h630.06832.2306.7560.1018丙烷c3h844.09496.7394.9990.157異丁烷i-c4h1058.120134.993.6000.18

23、4正丁烷n-c4h1058.120152.033.9210.197異戊烷i-c5h1272.146187.243.2470.2286正戊烷n-c5h1272.146196.343.2150.252己烷c6h1486.172234.132.7170.302氮氣n228.013-47.1801.2800.0372二氧化碳co244.01031.107.2900.2667水h2o18.015373.8521.8300.3434將bwrs方程應(yīng)用于混合物時各參數(shù)的混合規(guī)則如下 式中 汽相或液相混合物組分的摩爾分?jǐn)?shù)； 組分間的交互作用參數(shù)（=，=0），表2.3列出了starling所給出的18種組分間的

24、數(shù)據(jù)。表2.3 bwrs方程中使用的二元交互作用參數(shù)值（=）物質(zhì)甲烷乙烷丙烷異丁烷正丁烷異戊烷正戊烷己烷氮氣二氧化碳甲烷0.00.010.0230.02750.0310.0360.0410.050.0250.05乙烷0.010.00.00310.0040.00450.0050.0060.0070.070.048丙烷0.0230.00310.00.0030.00350.0040.00450.0050.100.045異丁烷0.02750.0040.0030.00.00.0080.0010.00150.110.05正丁烷0.0310.00450.00350.00.00.0080.0010.00150

25、.120.05異戊烷0.0360.0050.0040.0080.0080.00.00.00.1340.05正戊烷0.0410.0060.00450.0010.0010.00.00.00.1480.05己烷0.050.0070.0050.00150.00150.00.00.00.2000.05氮氣0.0250.070.100.110.120.1340.1480.2000.00.0二氧化碳0.050.0480.0450.050.050.050.050.050.00.0將計算出的bwrs方程的114個狀態(tài)參數(shù)以及已知節(jié)流前的氣體壓力和溫度代入中計算出氣體密度。為了方便求解，可將式改為的形式，并采用式

26、的迭代格式進(jìn)行迭代求解： （=1，2， ） 取迭代初值為=0.0，,迭代到為止。這里取=，一般情況迭代46次即可收斂。將定壓摩爾熱容和已知節(jié)流前的壓力、溫度以及求出的密度代入，即可求出所給氣體在所處壓力和溫度條件下的焦耳湯姆遜系數(shù)。節(jié)流溫降確定可寫為以下這個式子： 將求得的焦耳湯姆遜系數(shù)代入即可求得天然氣節(jié)流后的溫度。3天然氣水合物生成條件及機(jī)理研究水合物是目前國際上的熱門研究領(lǐng)域，由于其不僅與石油天然氣生產(chǎn)、儲藏和運輸密切相關(guān)，而且與環(huán)境保護(hù)、氣候變遷，特別是人類未來賴以生存的能源有關(guān)系。美國地質(zhì)調(diào)查局經(jīng)過多年調(diào)查研究認(rèn)為，天然氣水合物將可能成為本世紀(jì)的一種主要能源，是能源的“一個新領(lǐng)域”。

27、按保守估計，“全世界的天然氣水合物形式存在的碳的總量是地球上已知化石燃料（包括煤）中碳含量的2倍”。所以一方面由于水合物特殊的物理化學(xué)性質(zhì)，使其潛藏著開發(fā)工業(yè)應(yīng)用技術(shù)的巨大潛能，因此它不僅得到眾多學(xué)術(shù)研究機(jī)構(gòu)的重視，而且得到許多政府部門的重視；但是另一方面因為天然氣水合物會給天然氣的開采儲運造成很大的影響甚至危害，需要采取措施防止其生成，因此研究天然氣水合物生成條件及其機(jī)理是非常有必要的。3.1天然氣水合物的物理性質(zhì)及結(jié)構(gòu)天然氣水合物是由甲烷、乙烷、丙烷、異丁烷、正丁烷及n2、co2、h2s等分子在一定溫度和壓力條件下，與游離水結(jié)合而形成的結(jié)晶籠狀固體。其中水分子（主體分子）借助氫鍵形成主體結(jié)

28、晶網(wǎng)絡(luò)，晶格中的孔穴充滿輕烴或非輕烴分子（客體分子），主體分子和客體分子之間依靠范德華力形成穩(wěn)定性不同的水合物，氣體分子和水分子之間沒有化學(xué)計量關(guān)系，其穩(wěn)定性和結(jié)構(gòu)與客體分子的大小、種類及外界等因素有關(guān)。水合物在外觀上是白色或其它顏色的結(jié)晶物體。盡管天然氣水合物是由無色的碳?xì)錃怏w和水分子組成，但它并非全都呈現(xiàn)為白色，還有多種其他的色彩。一些從墨西哥灣海底獲取的天然氣水合物就具有黃色、橙色以及紅色等多種很鮮艷的顏色，而從大西洋海底取得的夭然氣水合物則呈現(xiàn)為灰色或藍(lán)色。由天然氣與水組成的結(jié)晶狀具有籠形構(gòu)造的天然氣水合物，是在特定的低溫高壓物理化學(xué)條件下形成的產(chǎn)物，因此其嚴(yán)格受溫壓參數(shù)所控制，目前所

29、公認(rèn)的形成天然氣水合物的溫壓條件是：溫度2.5250，但如果中重?zé)N和非烴類雜質(zhì)含量增加，則溫度范圍尚可拓寬，壓力一般為540mpa。據(jù)估算，水合物水分子和氣體分子的摩爾數(shù)分別占85%和15%，因此其密度和冰大致相等，但硬度和剪切模量小于冰，熱導(dǎo)率和電阻率遠(yuǎn)小于冰且具有多孔性。其化學(xué)成分不穩(wěn)定，通式為mnh2o，m為水合物中的氣體分子，n為水分子數(shù)。水合物的物理性質(zhì)見表3.1。表3.1 甲烷天然氣水合物和冰的性質(zhì)性質(zhì)甲烷天然氣水合物泥沙沉積物中海底甲烷天然氣水合物冰硬度(mohs)2-474剪切強度(mpa)12.27剪切模量2.43.9密度(g/cm3)0.9110.917聲學(xué)速率(m/s)3

30、30038003500熱容(kj/cm3)-273k2.322.3熱傳導(dǎo)率(w/mk)0.50.52.23電阻率(km)51005003.2天然氣水合物的分類天然氣水合物可以有兩種分類方法：一種是按產(chǎn)出環(huán)境或溫度壓力機(jī)制分類，一種是按其結(jié)構(gòu)類型分類。按產(chǎn)出環(huán)境，天然氣水合物可以分為海底天然氣水合物和極地天然氣水合物兩種類型。這兩種產(chǎn)出環(huán)境也代表著兩種截然不同的壓力溫度機(jī)制。通常把在海洋過渡帶、邊緣海和內(nèi)陸海等海洋底蘊藏的天然氣水合物都稱為海底天然氣水合物。在海洋中，甲烷天然氣水合物是儲量最豐富的一種類型，常出現(xiàn)在深海中或極地大陸上。1m3完全飽和的甲烷天然氣水合物（在其中100%的客體位置均被

31、占據(jù)）包含164m3的甲烷（在標(biāo)準(zhǔn)溫度、壓力下）和0.87m3的水。極地天然氣水合物是在較低的壓力和溫度下形成的，蘊藏的深度相對比較淺。極地天然氣水合物可作為水冰混合物出現(xiàn)在陸地的永久凍土帶或大陸架的永久凍土帶，在永久凍土帶之下的油氣田中也可能出現(xiàn)。由于極地天然氣水合物的分布在很大程度上受地域的限制，因而其總量少于海底天然氣水合物。按結(jié)構(gòu)類型，天然氣水合物可分為型、型和h型三種結(jié)構(gòu)。天然氣水合物結(jié)構(gòu)與性能是一個研究比較成熟的領(lǐng)域，它可被看作是一類主客體材料。結(jié)構(gòu)型氣水合物為立方晶體結(jié)構(gòu)，其在自然界分布最為廣泛，僅能容納甲烷（c1）、乙烷這兩種小分子的烴以及n2、co2、h2s等非烴分子，這種水

32、合物中甲烷普遍存在的形式是構(gòu)成ch4nh2o的幾何格架；結(jié)構(gòu)型水合物屬于金剛石晶體立方結(jié)構(gòu)，除包容c1、c2等小分子外，較大的“籠子”（水合物晶體中水分子間的空穴）還可以容納丙烷（c3）及異丁烷（i-c4）等烴類；結(jié)構(gòu)h型，屬于六面體結(jié)構(gòu)（取英文字母h代表），其大的“籠子”甚至可以容納直徑超過異丁烷（i-c4）的分子，如i-c5和其它直徑在7.58.6之間的分子，也可由揮發(fā)油和汽油等大分子形成。結(jié)構(gòu)h型氣水合物早期僅存在于實驗室，1993年才在墨西哥灣大陸斜坡發(fā)現(xiàn)其天然產(chǎn)物。實驗表明型和h型水合物比型水合物更穩(wěn)定。正五邊形的十二面體（d-512）在三種結(jié)構(gòu)中均存在，而十四面體（t-51262）

33、、十六面體（h-51264）和二十面體（e-51264）分別是、和h型結(jié)構(gòu)中的大籠子。另外，在h型結(jié)構(gòu)水合物中尚有另外一種十二面體（d-435663）籠子。實際三類結(jié)構(gòu)中決定晶格的籠的排列次序是非常不同的：型結(jié)構(gòu)水合物由特定的t籠頂點確定；而和h型結(jié)構(gòu)是基于被氫鍵連接的d籠不同的堆迭次序。h型水合物的晶胞由3個d籠，2個d籠和1個e籠堆砌而成，共含34個分子，晶胞分子式用s3s2l34h2o來表示。圖3-1是三種水合物晶體結(jié)構(gòu)和相應(yīng)的籠，五種單元。表3.2是水合物的晶格常數(shù)。表3.2 水合物的晶格常數(shù)參 數(shù)結(jié)構(gòu)結(jié)構(gòu)結(jié)構(gòu)h單位晶胞中水分子數(shù)4613634單位晶胞中校孔穴數(shù)2163單位晶胞中中等孔

34、穴數(shù)/2單位晶胞中大孔穴數(shù)681小孔穴常數(shù)1/232/17/大孔穴常數(shù)3/231/17/晶格孔穴直徑（）小7.957.827.82中/8.12大8.609.4011.423.2.1天然氣工業(yè)中可能形成的水合物結(jié)構(gòu)天然氣水合物型、型的研究理論已經(jīng)比較成熟，1987年ripmeester等人發(fā)現(xiàn)h型天然氣水合物，以來h型水合物的研究不斷受到重視，而我國目前對h型水合物的研究還不多。一般認(rèn)為含有58個碳原子的烴類物質(zhì)(分子直徑小于9)在小分子(如甲烷、乙烷)烴類氣體的幫助下可形成h型水合物。在油氣中常含有一些能形成h型水合物的烴類物質(zhì) (如甲基環(huán)戊烷、異戊烷等)，因此討論三種類型型水合物是否在一般的

35、油氣輸送管道以及設(shè)備中穩(wěn)定存在是有必要的。對于自然產(chǎn)出的氣體，氣體分子只要能進(jìn)入三種水合物晶體結(jié)構(gòu)的五種氣體占據(jù)空間，并且使天然氣水合物結(jié)構(gòu)達(dá)到穩(wěn)定，那么這種氣體就可以形成天然氣水合物。圖3.1表示了自然產(chǎn)出的氣體分子大小范圍、氣體所占據(jù)的空間類型及形成的水合物類型。從圖可見，直徑小于0.35nm的分子對于任何氣體占據(jù)空間都是太小了，因而不能形成天然氣水合物，如氫、氦和氖。直徑在0.350.42nm之間的氣體，如氮、氧、氬和氪氣體不僅能充填512空間，而且也能充填51264空間，形成二種空間均有氣體充填的型純氣體水合物。直徑為0.420.6nm氣體硫化氫、甲烷同時能充填512和51262二種空

36、間，形成型結(jié)構(gòu)天然氣水合物；乙烷和二氧化碳僅能占據(jù)51262空間，但不能穩(wěn)定512空間，形成512空間空的型天然氣水合物。氣體直徑為0.60.7nm的氣體的丙烷和異丁烷能占據(jù)型結(jié)構(gòu)51262空間，512小空間沒有氣體占據(jù)，形成512小空間空缺的型純氣體水合物。直徑為0.70.9nm的大分子烴對于和型水合物的氣體占據(jù)空間是太大了，不能形成和型水合物，在有小分子氣體加入時才能形成h型水合物，如正丁烷、2-甲基丁烷、甲基環(huán)戊烷、甲基環(huán)已烷、環(huán)庚烷、環(huán)辛烷、金剛烷等充填51268大空間，512和435663二種空間需要其它氣體小分子如氮、硫化氫和甲烷來充填才能穩(wěn)定。因此，大分子烴類不能形成純水合物。0

37、.9nm0.7nm0.6nm0.35nm0.42nmn- c4h10iso-c4h10c-c3h6(ch2)3oc3h8c2h6co2ch4xe，h2sarkro2n20.50.4分子直徑/nm0.30.60.70.8無水合物區(qū)型水合物512+51264型水合物512+51262型水合物51264h型水合物512+435663+51262圖3.1 分子大小和氣體占據(jù)空間類型及形成的水合物類型通過上述討論，我們能夠發(fā)現(xiàn)石油天然氣工業(yè)中的天然氣水合物結(jié)構(gòu)一般為型和型。對于管輸天然氣組成來說，生成結(jié)構(gòu)型的氣體主要是ch4、c2h6等碳分子數(shù)小于3的烴類和n2、o2的單組分或它們的混合物；一般當(dāng)氣體分

38、子直徑比c3h8大而比c4h8小或氣體混合物中含這種氣體組分時，生成的氣體水合物通常為型結(jié)構(gòu)；而在輸氣管道工程中形成h型水合物的幾率是很低的。3.3天然氣水合物生成條件一般水合物的生成條件有：天然氣中含有足夠的水分，以形成空穴結(jié)構(gòu)；具有一定的溫度與壓力；氣體處于脈動紊流等激烈擾動中，并有結(jié)晶中心存在。對輸氣管道中水合物的形成條件，從以下兩方面即熱力學(xué)方面和動力學(xué)方面進(jìn)行分析。3.3.1天然氣水合物形成的熱力學(xué)條件天然氣水合物的生成除與天然氣的組分、組成和游離水含量有關(guān)外，還需要一定的熱力學(xué)條件，即溫度和壓力?？捎孟旅娣匠瘫硎境鏊衔锷蓷l件：因而，生成水合物的第一個條件是：也就是說，只要當(dāng)系統(tǒng)

39、中氣體壓力（）大于它的水合物分解壓力（）時，才能有被水蒸氣飽和的氣體（）自發(fā)的生成水合物。嚴(yán)格地講，上式應(yīng)用逸度表示如下：水合物生成的第二個條件是：由第二個條件可以看出：從熱力學(xué)觀點來看，水合物的自發(fā)形成不是必須使氣體（）被飽和，只要系統(tǒng)中水的蒸汽壓大于水合物晶格表面的水蒸氣壓力就夠了?？傊?，水合物的生成需要三個條件，兩個主要條件是：天然氣中含有足夠的水分，以形成空穴結(jié)構(gòu)；具有一定的溫度與壓力條件，如高壓和低溫；一個輔助條件：氣體處于脈動紊流等激烈擾動中、酸性氣體的存在、微小水合晶核的誘導(dǎo)以及晶核停留的特定物理位置如彎頭、孔板、閥門、粗糙的管壁等。從理論上講，在形成水合物時，不一定需要游離水，

40、只要存在氣相或冷凝碳?xì)浠衔镏杏行纬伤衔锏慕M分共同存在，壓力和溫度條件滿足（即高壓和低溫），水和一些組分就會形成固體水合物，遠(yuǎn)高于水的冰點或游離水。水合物形成的臨界溫度，可能是水合物存在的最高溫度，高于此溫度，不管壓力多大，也不會形成水合物。天然氣各組分生成水合物的臨界溫度見表3.3。表3.3 天然氣各組分生成水合物的臨界形成溫度組分名稱ch4c2h6c3h8ic4h10nc4h10co2n2臨界溫度（）21.514.55.52.51.010.029.0曾以為甲烷水合物的臨界溫度是21.5，但經(jīng)研究表明當(dāng)壓力在3376mpa條件下，溫度為28.8時甲烷水合物仍然存在，而在390mpa高壓下甲

41、烷水合物形成溫度可提升到47。從熱力學(xué)方面講，天然氣水合物的形成主要與以下條件有關(guān)：先決條件是管道中含有形成水合物的主、客體物質(zhì)組成水分，小分子烴類，h2o，co2 等；輸送過程的工況條件即溫度壓力條件符合形成穩(wěn)定水合物的要求；氣流中存在形成水合物的結(jié)晶中心。形成水合物的客體物質(zhì)(烴類及雜質(zhì))大量研究表明，烴類物質(zhì)并不是全部都可以形成水合物的，直鏈烷烴中只有 ch4，c2h6 和 c3h8 能支持水合物的形成，支鏈烷烴中只有異丁烷能支持。此外，天然氣中的雜質(zhì)組分如 h2s，co2，n2 和 o2 等也支持水合物的形成。不同組成的天然氣其相對密度也不同，相對密度愈高的天然氣愈容易形成水合物，而且研究表明：重質(zhì)烴的很少含量就會使開始水合形成需要的壓力大幅度降低。形成水合物的主體物質(zhì)(水分)天然氣在管道中流動隨著壓力溫度的變化，其飽和含水量也發(fā)生著變化，如圖3.2所示曲線由 、為天然氣管道壓力溫度分布曲線，曲線為對應(yīng)輸氣管道的壓力溫度下的飽和含水量曲線。在輸氣管的前半段，壓力下降不大，而溫度急劇下降，飽和含水量也隨之下降，如段；在輸氣管后半段，溫度下降平緩，接近于周圍介質(zhì)溫度，而壓力則急劇下降，對應(yīng)的飽和含水量逐步上升，如段。其中點的飽和含水量最小為。 氣體的壓力溫度和含水量/mpa/距離/m