劉春花 郝忠獻(xiàn) 劉新福 何鴻銘 魏松波 周超

(1. 中國(guó)石油大學(xué)(華東) 機(jī)電工程學(xué)院 2.中國(guó)石油勘探開發(fā)研究院3. 青島理工大學(xué)機(jī)械與汽車工程學(xué)院 4. 中海油研究總院有限責(zé)任公司)

0 引 言

天然氣水合物開發(fā)、運(yùn)輸及加工與相平衡息息相關(guān)[1-3]。天然氣在井筒和管道流動(dòng)過程中處于高壓、低溫的環(huán)境，極易形成水合物并不斷聚集，輕則降低井液的流量，重則堵塞井筒，導(dǎo)致關(guān)井[4-6]。因此，準(zhǔn)確預(yù)測(cè)水合物的生成條件對(duì)防治致密天然氣水合物生成和保障井筒流動(dòng)安全具有重要意義[7-8]。

天然氣水合物生成的預(yù)測(cè)方法主要有經(jīng)驗(yàn)圖解法、相平衡計(jì)算法、統(tǒng)計(jì)熱力學(xué)計(jì)算法和智能算法。除智能算法外，其他模型大都是在經(jīng)典Vander Waals-Platteeuw(VDW-P)理論上發(fā)展和補(bǔ)充的。VDW-P模型主要預(yù)測(cè)天然氣水合物穩(wěn)定存在時(shí)的溫度和壓力條件，并計(jì)算氣體組分逸度和Langmuir常數(shù)等參數(shù)。該模型對(duì)含氮?dú)夂退嵝詺怏w的天然氣水合物平衡條件預(yù)測(cè)誤差較大，且模型計(jì)算較為復(fù)雜。各國(guó)學(xué)者提出了經(jīng)驗(yàn)公式來預(yù)測(cè)水合物生成條件，比如：J.JAVANMARDI等[9]考慮純液態(tài)水和空氣水合物晶格體積變化對(duì)水合物生成的影響，探討了高壓時(shí)VDW-P模型的適用性；H.DELAVAR等[10]考慮氣體在液相中的溶解度并建立相平衡模型，將Huron-Vidar混合規(guī)則應(yīng)用于SRK方程并修正了UNIQUAC模型參數(shù)；J.B. KLAUDA等[11]應(yīng)用PRSV狀態(tài)方程計(jì)算了水與其他組分在氣相的逸度，并修正了UNIFAC模型的參數(shù)；孫寶江等[12]采用UNIQUAC模型確定醇類體系的γw項(xiàng)，建立了含醇類天然氣的水合物生成預(yù)測(cè)模型。DU Y.H.等[13]改進(jìn)了Parrish-Prausnitz模型，并預(yù)測(cè)了含有丙烷或丙烷以上烴組分及酸性氣體的天然氣水合物生成條件。

目前，采用Van Laar混合規(guī)則的PRSV2狀態(tài)方程預(yù)測(cè)并分析天然氣水合物生成條件的研究較少，為此，本文提出基于Van Laar混合規(guī)則PRSV2狀態(tài)方程的液相和氣相逸度計(jì)算方法，建立了致密天然氣水合物相平衡動(dòng)態(tài)模型，并依據(jù)數(shù)值模擬和試驗(yàn)測(cè)試分析預(yù)測(cè)了致密天然氣水合物生成條件，揭示了不同組分及混合組分氣體對(duì)致密氣水合物生成的影響。

1 致密天然氣水合物相平衡模型

水合物相平衡模型可以模擬水合物-液相-氣相以及水合物-冰相-氣相的三相平衡方程，如式(1)所示。其中下標(biāo)w代表水相，上標(biāo)H代表水合物相，上標(biāo)π代表液相或冰相，上標(biāo)V代表氣相，判斷三相共存的依據(jù)為平衡階段各組分逸度相等。

(1)

1.1 液相與氣相的熱力學(xué)模型

區(qū)別于立方狀態(tài)方程描述氣相逸度和活度系數(shù)描述液相，為了計(jì)算平衡階段非水合物相中各組分的逸度，建立相應(yīng)的熱力學(xué)模型，并采用Stryjek和Vera修正的Peng-Robinso狀態(tài)方程(PRSV)(見式(2))對(duì)液相和氣相同時(shí)進(jìn)行描述。

(2)

其中：

(3)

(4)

(5)

系數(shù)k與溫度有關(guān)，且系數(shù)k的計(jì)算式為：

(6)

k0=0.379+1.490ω-0.171ω2+0.019 7ω3

(7)

式中：a、b分別為吸引力參數(shù)和體積參數(shù)；R為氣體常數(shù)，J/(K·mol)；T為熱力學(xué)溫度，K；p為壓力，MPa；V為摩爾體積，L/mol；pc為臨界壓力，MPa；Tc為臨界溫度，K；ω為偏心因子；α為計(jì)算中間變量；k0為偏心因子函數(shù)；k1、k2和k3均為純組分系數(shù)；k為與溫度相關(guān)的系數(shù)。

區(qū)別于PR方程，PRSV狀態(tài)方程對(duì)于各純化合物都有額外的k1、k2和k3參數(shù)可調(diào)[14-15]，且混合組分氣體的PRSV方程需結(jié)合Van Laar型混合規(guī)則使用。Van Laar型混合規(guī)則的PRSV狀態(tài)方程分別為：

(8)

(9)

(10)

式中：am為混合組分的吸引力參數(shù)，aij和aji為二元相互作用項(xiàng)，bm為混合組分的體積參數(shù)，i和j分別為組分i和組分j，kij與kji為二元交互參數(shù)[16]。

PRSV狀態(tài)方程中，組分i逸度的計(jì)算式為：

fi(T，p，x)=pφixi

(11)

其中：

(12)

式中：φi為組分i的逸度系數(shù)；Z為壓縮因子[17]；Am與Bm均為無量綱系數(shù)，Am=ap/(R2T2)，Bm=bp/(RT)。

1.2 水合物相的熱力學(xué)模型

(13)

其中：

(14)

(15)

(16)

(17)

Γ(r)為勢(shì)能模型并采用Kihara勢(shì)能，即：

(18)

式中：a為分子核半徑，m；r為兩分子核間距，m；σ為勢(shì)能等于0時(shí)兩分子間距，m；ε為最大引力勢(shì)，J/kg。

客體分子與孔穴壁z個(gè)水分子間總勢(shì)能為：

(19)

(20)

式中：N為指數(shù)，D為孔穴直徑，z為孔穴配位數(shù)[18]，δ為計(jì)算中間變量。

Kihara勢(shì)能的相關(guān)計(jì)算參數(shù)見表1。

表1 Kihara勢(shì)能的計(jì)算參數(shù)Table 1 Calculation parameters of Kihara potential energy

1.3 相平衡動(dòng)態(tài)模型求解方法

生成Ⅰ型致密天然氣水合物時(shí)，輸入溫度T與混合體系的摩爾組成xi，設(shè)定初始?jí)毫0，并采用等溫閃蒸方法計(jì)算氣液平衡[19]。致密氣水合物相平衡動(dòng)態(tài)模型求解方法為：

(1)依據(jù)式(17)、式(19)和式(20)計(jì)算溫度T時(shí)每種組分的Langmuir常數(shù)Cml，并通過式(11)和式(12)計(jì)算流體相中水的逸度；

(2)由式(16)計(jì)算m型孔穴組分l占有率θml；

(3)由式(13)～式(15)計(jì)算空氣水合物逸度；

(5)生成Ⅱ型致密天然氣水合物時(shí)，則按照步驟(1)～(4)依次進(jìn)行計(jì)算，直至體系達(dá)到平衡；

(6)體系平衡最低壓力為水合物生成壓力。

2 致密氣水合物生成試驗(yàn)測(cè)試分析

2.1 致密天然氣水合物生成試驗(yàn)裝置及方法

為了量化純組分、混合組分、溫度及壓力等試驗(yàn)參數(shù)對(duì)水合物生成的影響程度，同時(shí)為防治水合物生成和保障井筒流動(dòng)安全提供依據(jù)，開展了致密氣水合物生成規(guī)律試驗(yàn)。該試驗(yàn)裝置主要由高壓反應(yīng)釜、氣體控制單元、溫度控制單元、溫度和壓力數(shù)據(jù)采集系統(tǒng)及其輔助模塊等組成，如圖1所示。

圖1 致密天然氣水合物生成規(guī)律試驗(yàn)流程圖Fig.1 Flow chart for test of tight gas hydrate generation

試驗(yàn)采用DTA技術(shù)測(cè)定水合物生成，高壓反應(yīng)釜(容積50 mL)、進(jìn)口閥及出口閥均置入恒溫試驗(yàn)控制箱，試驗(yàn)設(shè)備容積提供等溫環(huán)境，最大程度減小水合物分解試驗(yàn)過程中外部條件的影響。試驗(yàn)前將蒸餾水純化至電阻率18 MΩ·cm，并在真空中脫氣12 h，所有測(cè)試氣體混合并在試驗(yàn)前靜置1周時(shí)間。

致密天然氣水合物生成規(guī)律試驗(yàn)不同混合組分氣體的質(zhì)量分?jǐn)?shù)見表2。

表2 致密天然氣水合物生成試驗(yàn)參數(shù)Table 2 Test parameters of tight gas hydrate generation

試驗(yàn)流程為：高壓反應(yīng)釜抽真空并注入20 mL蒸餾水，系統(tǒng)溫度穩(wěn)定時(shí)，打開截止閥并提供純組分/混合組分氣體，隨后系統(tǒng)加壓直至水合物生成；高壓反應(yīng)釜內(nèi)流體以250 r/min的速度混合，稍微開啟出口閥并使氣體緩慢流過高壓反應(yīng)釜且連續(xù)流動(dòng)6 h，排除氣體濃度對(duì)水合物生成條件的影響；氣體飽和后，轉(zhuǎn)速增至400 r/min，測(cè)定水合物生成點(diǎn)后，加熱系統(tǒng)并解離水合物；系統(tǒng)升溫并高于試驗(yàn)生成點(diǎn)0.5 K，隨后每3 h降溫0.1 K，直至系統(tǒng)達(dá)到水合物生成溫度，重復(fù)該過程3次并取平均值作為水合物在該壓力時(shí)的生成溫度。

該試驗(yàn)采用控制變量法，在設(shè)定好各試驗(yàn)參數(shù)的情況下，改變其中某個(gè)影響因素，氣體通過壓縮機(jī)注入，氣體流量由氣體控制單元自動(dòng)調(diào)節(jié)并保持恒定，試驗(yàn)溫度由溫度控制單元自動(dòng)控制。試驗(yàn)過程中，對(duì)模型的壓力、壓差、溫度及流量進(jìn)行實(shí)時(shí)監(jiān)測(cè)，并由數(shù)據(jù)采集系統(tǒng)記錄和存儲(chǔ)。

2.2 純組分水合物生成試驗(yàn)結(jié)果與分析

圖2為純甲烷與純乙烷組分時(shí)室內(nèi)試驗(yàn)和理論分析的水合物生成對(duì)比情況。從圖2可以看出，純組分天然氣的水合物生成壓力隨著溫度的升高而不斷增加，純甲烷組分天然氣的水合物生成曲線近似呈指數(shù)增長(zhǎng)，而純乙烷組分天然氣的水合物生成曲線在290 K附近出現(xiàn)轉(zhuǎn)折點(diǎn)并迅速上升，其主要原因?yàn)?90 K時(shí)乙烷氣體發(fā)生了相變，且相態(tài)由Lw-H-V變?yōu)長(zhǎng)w-H-LC2H6。溫度低于290 K時(shí)，甲烷與乙烷的水合物生成壓力大致相同，而溫度為291～320 K時(shí)，乙烷更難生成水合物，此時(shí)甲烷相態(tài)為L(zhǎng)w-H-V，而乙烷相態(tài)則變?yōu)長(zhǎng)w-H-LC2H6。由于液態(tài)烴更難以生成水合物，所以溫度高于321 K后，甲烷更難生成水合物。

圖2 純甲烷與純乙烷組分致密天然氣水合物生成曲線Fig.2 Generating curve of tight gas hydrate with pure methane and pure ethane components

圖3為純丙烷與純異丁烷組分時(shí)試驗(yàn)和理論分析的水合物生成對(duì)比情況。溫度低于275 K時(shí)，純丙烷與純異丁烷組分天然氣水合物的生成壓力幾乎一致，純異丁烷組分天然氣水合物生成曲線在275 K附近出現(xiàn)轉(zhuǎn)折點(diǎn)，隨后驟然上升，而丙烷在280 K后水合物生成壓力才開始驟然上升，其主要原因仍為純組分天然氣發(fā)生相變，且相同條件的丙烷比異丁烷更易生成水合物。

圖3 純丙烷與純異丁烷組分致密天然氣水合物生成曲線Fig.3 Generating curve of tight gas hydrate with pure propane and pure isobutane components

2.3 混合組分水合物生成試驗(yàn)結(jié)果與分析

圖4為甲烷與乙烷的混流一和甲烷與丙烷混流二的水合物生成試驗(yàn)和理論分析結(jié)果對(duì)比情況，各氣體組分質(zhì)量分?jǐn)?shù)見表2。從圖4可見，溫度為280～305 K時(shí)，甲烷與乙烷的混流一比甲烷與丙烷的混流二更難生成水合物。溫度低于285 K時(shí)，兩種混流的天然氣水合物生成壓力相差不大，而溫度高于285 K之后，兩者的水合物生成壓力差值變大。對(duì)比分析圖1～圖3可知，相對(duì)于純組分氣體，混合了乙烷或丙烷的天然氣水合物生成壓力差值更易在某個(gè)溫度拐點(diǎn)之后急劇增長(zhǎng)。

圖4 甲烷與乙烷和丙烷混流致密天然氣水合物生成曲線Fig.4 Generating curve of tight gas hydrate with mixed flow of methane, ethane and propane

圖5為甲烷與正丁烷的混流三和甲烷與異丁烷的混流四的水合物生成試驗(yàn)和理論分析結(jié)果對(duì)比情況，各氣體組分質(zhì)量分?jǐn)?shù)見表2。從圖5可以看出，甲烷與正丁烷的混流三和甲烷與異丁烷混流四的天然氣水合物生成曲線差別較大，甲烷與正丁烷的混流三在溫度升至279 K以后，水合物生成壓力開始快速上升，而甲烷與異丁烷混流四形成天然氣水合物的生成壓力上升趨勢(shì)不明顯。

圖5 甲烷與正丁烷和異丁烷不同混流水合物生成曲線Fig.5 Generating curve of hydrate with different mixed flow of methane, normal butane and isobutane

圖6為CH4與CO2的混流五、CH4和N2的混流六、CH4與H2S的混流七以及CH4與CO2、N2和H2S混流八的水合物生成試驗(yàn)和理論分析結(jié)果對(duì)比情況，各氣體組分質(zhì)量分?jǐn)?shù)見表2。從圖6可以看出，相同條件下，天然氣混入H2S后的水合物生成壓力顯著減小，說明H2S對(duì)天然氣水合物的生成有促進(jìn)作用。溫度為275～300 K時(shí)，CH4與CO2的混流五與CH4和N2混流六的水合物生成壓力均高于純CH4組分天然氣，說明CO2與N2對(duì)水合物的生成均有抑制作用。相同溫度時(shí)，CH4與H2S的混流七和CH4與CO2、N2和H2S混流八的水合物生成壓力均低于純CH4的水合物生成壓力，說明天然氣混入CO2、N2和H2S后，作為酸性氣體的H2S對(duì)天然氣水合物生成的影響最大并起主導(dǎo)作用。

圖6 甲烷與酸性氣體不同混流致密氣水合物生成曲線Fig.6 Generating curve of tight gas hydrate with different mixed flow of methane and sour gas

3 致密天然氣井場(chǎng)試驗(yàn)與分析

3.1 致密天然氣井基本參數(shù)

選取鄂爾多斯盆地東緣大寧-吉縣區(qū)塊3口致密天然氣井開展天然氣水合物生成井場(chǎng)試驗(yàn)，這3口致密天然氣井依次為DJ-04氣井、DJ-12氣井和DJ-50氣井，測(cè)定不同混合組分氣體的水合物生成溫度和壓力等參數(shù)，依據(jù)川慶鉆探長(zhǎng)慶井下技術(shù)作業(yè)公司的氣樣檢測(cè)報(bào)告，3口井的井場(chǎng)測(cè)試基本參數(shù)見表3。

表3 大寧-吉縣區(qū)塊3口致密天然氣井排采參數(shù)Table 3 Parameters for drainage gas recovery in 3 tight gas wells in Daning-Jixian block

3.2 井場(chǎng)試驗(yàn)測(cè)試與實(shí)例計(jì)算分析

圖7為不同致密天然氣井混流的水合物生成井場(chǎng)試驗(yàn)和理論分析結(jié)果對(duì)比情況。從圖7可以看出，隨著溫度的上升，致密天然氣的水合物生成壓力逐漸增大，DJ-04氣井和DJ-12氣井中，異丁烷組分的質(zhì)量分?jǐn)?shù)差別較大，且異丁烷質(zhì)量分?jǐn)?shù)高的DJ-04氣井致密天然氣更易生成水合物，說明異丁烷對(duì)天然氣水合物的生成有促進(jìn)作用。DJ-04氣井和DJ-50氣井中，N2組分的質(zhì)量分?jǐn)?shù)差別較大，且N2質(zhì)量分?jǐn)?shù)低的DJ-50氣井致密天然氣更易生成水合物，說明N2對(duì)天然氣水合物的生成有抑制作用，這與生成試驗(yàn)和理論分析結(jié)果一致。

圖7 井場(chǎng)測(cè)試不同混流天然氣水合物生成曲線Fig.7 Generating curve of gas hydrate with different mixed flow obtained at wellsite test

4 結(jié) 論

(1)基于Van Laar混合規(guī)則的PRSV2狀態(tài)方程，建立了致密天然氣水合物相平衡和熱力學(xué)模型，預(yù)測(cè)了致密氣水合物的生成條件，這對(duì)防治致密天然氣水合物生成和保障井筒流動(dòng)安全具有重要意義。

(2)溫度是影響水合物生成的主要因素，且致密天然氣水合物生成壓力隨著溫度的升高而呈近似指數(shù)的方式增加，純乙烷、純丙烷和純異丁烷的水合物分別在290、280和275 K附近發(fā)生相變并出現(xiàn)轉(zhuǎn)折點(diǎn)，隨后生成壓力迅速上升。

(3)以甲烷為主的二元組分和多元組分混流與純甲烷組分天然氣的對(duì)比分析結(jié)果表明，CO2與N2對(duì)水合物生成有抑制作用，而天然氣混入H2S或其他烷烴會(huì)極大降低水合物生成壓力，且在CO2、N2和H2S這3種氣體中H2S對(duì)水合物生成的影響最大并起主導(dǎo)作用。