超臨界二氧化碳射流計(jì)算模型的建立與射流特性分析*

孫 雪 倪紅堅(jiān) 王恵文 沈忠厚

(1.中國(guó)石油大學(xué)(華東)非常規(guī)油氣與新能源研究院 山東青島 266580； 2.中國(guó)石化中原石油工程有限公司鉆井技術(shù)研究院 河南濮陽(yáng) 457001)

天然水合物開(kāi)采技術(shù)已成為目前世界科技前沿創(chuàng)新技術(shù)研究熱點(diǎn)[1-4]，其中二氧化碳置換法開(kāi)采海底天然氣水合物是近年來(lái)海洋石油工程領(lǐng)域的研究熱點(diǎn)和前沿技術(shù)，而超臨界二氧化碳射流作為一種輔助技術(shù)在開(kāi)采天然氣水合物方面有許多優(yōu)勢(shì)。超臨界二氧化碳是指溫度壓力高于臨界點(diǎn)(臨界溫度304.13 K，臨界壓力7.38 MPa)的二氧化碳相態(tài)，兼具氣態(tài)和液態(tài)的特性，具有高密度、低黏度、滲透性好的特點(diǎn)[5]。利用超臨界二氧化碳射流開(kāi)采海底天然氣水合物，可以增大二氧化碳和天然氣水合物的接觸面積，提高置換效率[6]；同時(shí)，由于射流的節(jié)流降溫效應(yīng)，降低了二氧化碳的溫度，更有利于二氧化碳水合物的形成，有助于開(kāi)采天然氣水合物及二氧化碳深海封存技術(shù)的實(shí)現(xiàn)。

超臨界二氧化碳射流技術(shù)首次在石油工程中提出是在頁(yè)巖氣領(lǐng)域。2000年，Kolle[7]提出了連續(xù)油管超臨界二氧化碳射流鉆井技術(shù)，并通過(guò)實(shí)驗(yàn)對(duì)比，發(fā)現(xiàn)超臨界二氧化碳射流可以在更低的門(mén)限壓力下破碎曼柯斯頁(yè)巖，且破碎效率是水射流的3倍。后來(lái)，在國(guó)內(nèi)杜玉昆 等[8]、黃飛 等[9]分別針對(duì)人工巖心和四川龍馬溪頁(yè)巖開(kāi)展超臨界二氧化碳射流破巖實(shí)驗(yàn)，證明了超臨界二氧化碳射流輔助破巖技術(shù)的可行性；程宇雄 等[10-11]提出了超臨界二氧化碳連續(xù)油管噴射壓裂技術(shù)，并通過(guò)數(shù)值模擬論證了該技術(shù)的可行性，認(rèn)為超臨界二氧化碳噴射壓裂的增壓效果好，有助于降低起裂壓力，且不會(huì)引起黏土膨脹、水鎖等儲(chǔ)層傷害，是在頁(yè)巖氣壓裂改造中應(yīng)用前景良好的一種無(wú)水壓裂技術(shù)；杜玉昆 等[12]提出了一種利用超臨界二氧化碳射流開(kāi)采海底天然氣水合物的裝置和方法，利用二氧化碳射流鉆穿水合物上部的凍土層，并將二氧化碳水平噴射進(jìn)水合物層，置換出甲烷氣后，二氧化碳和甲烷從環(huán)空上返至平臺(tái)并分離,可以實(shí)現(xiàn)天然氣水合物鉆井和開(kāi)采的一體化，且易于實(shí)現(xiàn)欠平衡開(kāi)發(fā)，降低地層失穩(wěn)的風(fēng)險(xiǎn)。

顯然，超臨界二氧化碳射流技術(shù)在石油工業(yè)中應(yīng)用前景廣闊，而超臨界二氧化碳射流計(jì)算和特性研究是各項(xiàng)技術(shù)成功研制并應(yīng)用的基礎(chǔ),但目前超臨界二氧化碳射流特性的研究方法局限于實(shí)驗(yàn)測(cè)試和數(shù)值模擬。例如，Wang等[13]、Tian等[14]分別通過(guò)自主設(shè)計(jì)的實(shí)驗(yàn)裝置測(cè)試了超臨界二氧化碳射流下巖石的壓力和溫度分布，但實(shí)驗(yàn)裝置無(wú)法測(cè)量超臨界二氧化碳在噴嘴內(nèi)的流場(chǎng)變化；Wang等[15]通過(guò)數(shù)值模擬方法求解了超臨界二氧化碳沖擊射流流場(chǎng)，并分析了射流參數(shù)對(duì)出口流速的影響規(guī)律，但其模擬案例僅限于射流圍壓30 MPa以上的情況，而當(dāng)射流圍壓接近臨界壓力時(shí)二氧化碳可壓縮性強(qiáng)，數(shù)值模擬耗時(shí)長(zhǎng)，收斂精度差。由此可見(jiàn)，目前實(shí)驗(yàn)測(cè)試和數(shù)值模擬方法都不適合用于現(xiàn)場(chǎng)工藝和方案的設(shè)計(jì)，迫切需要一種精確而簡(jiǎn)便的射流特性求解方法。筆者將超臨界二氧化碳考慮為真實(shí)氣體，采用Span-Wagner方程描述二氧化碳的狀態(tài)方程和熱力學(xué)性質(zhì)，建立超臨界二氧化碳在噴嘴中的一維流動(dòng)模型，進(jìn)而得到了一種簡(jiǎn)單實(shí)用的超臨界二氧化碳射流特性的計(jì)算方法，可為超臨界二氧化碳射流技術(shù)的實(shí)驗(yàn)方案設(shè)計(jì)和工藝參數(shù)優(yōu)化提供理論支持。

1 超臨界二氧化碳射流計(jì)算模型建立與驗(yàn)證

1.1 模型假設(shè)

鉆井中常用的噴嘴形狀為漸縮型噴管接等直徑的圓管，因此將超臨界二氧化碳在噴嘴的流動(dòng)過(guò)程分為收縮段和圓柱段2個(gè)部分，并做以下假設(shè)：

1) 由于噴嘴中軸向流速很大，徑向流動(dòng)和切向流動(dòng)可忽略，因此，噴嘴的流動(dòng)過(guò)程簡(jiǎn)化為一維流動(dòng)。

2) 忽略流體與外界環(huán)境傳熱的影響，整個(gè)流動(dòng)過(guò)程可視為絕熱過(guò)程。

3) 由于收縮段流速和壓力變化大，主要是由于動(dòng)量改變引起的，因此收縮段可近似為等熵膨脹過(guò)程，而圓柱段可視為絕熱摩擦管流[16]。

1.2 模型建立

沿噴嘴軸線，任取一微元，則微元內(nèi)滿足一般流體的質(zhì)量守恒方程，即

d(ρuA)=0

(1)

式(1)中：ρ為流體的密度，kg/m3；u為流體的流速，m/s；A為噴嘴的截面積，m2。

由于噴嘴很短，忽略流體與外界的熱量交換，因此流體微元滿足一維定常絕熱流的能量守恒方程[17],即

(2)

式(2)中：H為流體的焓，J。

收縮段流體微元遵循等熵膨脹規(guī)律,即

dS=0

(3)

式(3)中：S為流體的熵，J/(mol·K)。

圓柱段滿足考慮摩擦的動(dòng)量守恒方程,即

(4)

式(4)中：p為流體壓力，MPa；D為噴嘴直徑，mm；f為二氧化碳流動(dòng)的摩擦阻力系數(shù)。

采用Wang等[18]通過(guò)實(shí)驗(yàn)擬合的超臨界二氧化碳摩阻系數(shù)的計(jì)算公式，即

(5)

由于水平段流速變化緩慢，任一微元內(nèi)du≈0，則式(2)、(4)可簡(jiǎn)化為

dH=0

(6)

(7)

采用Span和Wagner給出的二氧化碳物性方程組[20]求解二氧化碳的物性，即

(8)

(9)

-φ0-φr]

(10)

(11)

1.3 模型求解

在噴射鉆井施工過(guò)程中，噴嘴入口溫度、噴嘴壓降和圍壓(即射流出口壓力)是由地層深度和施工參數(shù)決定的。為求解噴射參數(shù)，可先假設(shè)1個(gè)入口馬赫數(shù)，即已知入口壓力、溫度、流速情況下，結(jié)合上述公式求解收縮段和圓柱段流動(dòng)參數(shù)。整個(gè)求解過(guò)程通過(guò)迭代入口馬赫數(shù)，以噴嘴出口壓力計(jì)算值與圍壓相等為收斂條件，求解步驟如下：

1) 將噴嘴內(nèi)部流動(dòng)考慮為一維，沿射流軸向劃分為n段，每段長(zhǎng)l。

2) 已知噴嘴入口壓力p1、溫度T1和圍壓pa，假設(shè)一個(gè)入口馬赫數(shù)M1，由式(8)、(11)可計(jì)算出入口流體的聲速C1，則入口流速u1=C1M1。

(12)

(13)

Si=Si-1

(14)

4) 假設(shè)i節(jié)點(diǎn)的壓力pi，焓Hi，依據(jù)噴嘴圓柱段守恒條件，由式(4)、(6)、(7)可得i節(jié)點(diǎn)的流動(dòng)參數(shù)，即

(15)

Hi=Hi-1

(16)

(17)

1.4 模型驗(yàn)證

為了驗(yàn)證模型的正確性，采用Wang等[15]提出的數(shù)值模擬方法模擬超臨界二氧化碳射流在噴嘴內(nèi)的流場(chǎng)，與本文模型的計(jì)算結(jié)果相對(duì)比。噴嘴幾何模型和計(jì)算參數(shù)均采用文獻(xiàn)[15]。噴嘴為漸縮型，該噴嘴幾何結(jié)構(gòu)如圖1所示。噴嘴的結(jié)構(gòu)參數(shù)為：入口直徑0.016 m、出口直徑0.006 m、收縮段長(zhǎng)度0.020 m、圓柱段長(zhǎng)度0.012 m；計(jì)算參數(shù)為：入口壓力60 MPa、入口溫度360 K、出口壓力30 MPa。

圖1 噴嘴幾何模型Fig.1 Geometric model of nozzle

超臨界二氧化碳射流在噴嘴內(nèi)的流速及溫度模擬計(jì)算結(jié)果如圖2所示。從圖2可以看出，本文建立的計(jì)算模型計(jì)算結(jié)果與參考文獻(xiàn)數(shù)值模擬結(jié)果吻合較好，噴嘴內(nèi)速度場(chǎng)誤差小于5.23%，溫度場(chǎng)誤差小于0.51%，說(shuō)明本文的計(jì)算模型可靠度高，可適用于現(xiàn)場(chǎng)工藝的參數(shù)優(yōu)化和方案設(shè)計(jì)中。

圖2 本文建立的計(jì)算模型與文獻(xiàn)[15]數(shù)值模擬結(jié)果對(duì)比Fig.2 Comparison of calculation results of established model in this paper and numerical simulation results from reference[15]

2 超臨界二氧化碳射流特性分析

與高壓水射流不同，射流速度并不是評(píng)價(jià)超臨界二氧化碳射流特性的唯一指標(biāo)。超臨界二氧化碳通過(guò)噴嘴后，不僅射流速度急劇上升，流體的物性也發(fā)生了很大的變化。與噴嘴入口流體相比，出口流體的溫度和密度下降，射流馬赫數(shù)增加。射流馬赫數(shù)是射流速度和流體聲速的比值，反映了流體的壓縮性。射流速度、射流馬赫數(shù)、流體溫度和密度共同反映了超臨界二氧化碳射流的性能，并決定了二氧化碳射流的沖蝕、混合和攜帶等能力。在超臨界二氧化碳射流的實(shí)際應(yīng)用中，射流特性受到噴嘴壓降、入口溫度和圍壓的共同影響，因此研究這些工況參數(shù)對(duì)射流特性分析是非常必要的。

2.1 噴嘴壓降對(duì)射流特性的影響

計(jì)算了入口溫度320 K、圍壓10 MPa條件下射流特性和射流物性隨噴嘴壓降的變化規(guī)律，如圖3、4所示。從圖3可以看出，隨著噴嘴壓降的增大，射流速度顯著增大，但噴嘴壓降對(duì)射流馬赫數(shù)的影響甚微。這是由于噴嘴壓降越大，意味著更多的壓能轉(zhuǎn)化為射流動(dòng)能，因此射流速度顯著增大。另一方面，噴嘴壓降越大，導(dǎo)致焦耳-湯姆遜節(jié)流效應(yīng)更強(qiáng)，射流出口流體溫度更低，流體密度增大(圖4)，噴射出的流體性質(zhì)更趨向于高密度的不可壓縮流體，噴射流體聲速的增加，因此射流馬赫數(shù)變化不大。從圖3、4可以看出，當(dāng)噴嘴壓降從20 MPa增加至50 MPa，射流速度增長(zhǎng)了49.1%，射流馬赫數(shù)增長(zhǎng)了19.19%，出口溫度降低了3.9%，出口密度增加了10.62%?？梢?jiàn)，噴嘴壓降對(duì)射流速度的影響顯著，而對(duì)出口流體的物性影響較小。

圖3 入口溫度320 K、圍壓10 MPa條件下射流特性隨噴嘴壓降的變化曲線Fig.3 Relation curve of jet speed and nozzle pressure drop (inject temperature 320 K and ambient pressure 10 MPa)

圖4 入口溫度320 K、圍壓10 MPa條件下射流物性隨噴嘴壓降的變化曲線Fig.4 Relation curve of jet properties and nozzle pressure drop (inject temperature 320 K and ambient pressure 10 MPa)

2.2 入口溫度對(duì)射流特性的影響

圖5、6分別為噴嘴壓降40 MPa、圍壓10 MPa條件下射流特性和射流物性隨入口溫度的變化曲線。從圖5、6可以看出，入口溫度的增加提高了入口流體的總能量，因此射流速度和射流馬赫數(shù)均有所增長(zhǎng)，并且入口溫度的變化對(duì)射流出口的流體物性影響較大。入口溫度越大，出口流體的密度越小，流體性質(zhì)更趨向于低密度的可壓縮流體，流體聲速減小，因此,隨著入口溫度的增加，射流馬赫數(shù)增長(zhǎng)幅度比射流速度高。在本文分析中，入口溫度從280 K增至340 K，射流速度增加了10.9%，射流馬赫數(shù)增加了103.25%，出口溫度增長(zhǎng)了15.28%，出口密度降低了23.8%?？梢?jiàn)，入口溫度對(duì)射流速度的影響較小，而對(duì)射流馬赫數(shù)的影響顯著，并且入口溫度的改變對(duì)出口流體的物性影響較大。

圖5 噴嘴壓降40 MPa、圍壓10 MPa條件下射流特性隨入口溫度的變化曲線Fig.5 Relation of jet speed and inject temperature (nozzle pressure drop 40 MPa and ambient pressure 10 MPa)

圖6 噴嘴壓降40 MPa、圍壓10 MPa條件下射流物性隨入口溫度的變化曲線Fig.6 Relation of jet properties and inject temperature (nozzle pressure drop 40 MPa and ambient pressure 10 MPa)

2.3 圍壓對(duì)射流特性的影響

分析了入口溫度320 K、噴嘴壓降40 MPa條件下射流特性和射流物性隨圍壓的變化規(guī)律，如圖7、8所示。

從圖7、8可以看出，隨著圍壓的增大，流體不可壓縮性增強(qiáng)，焦耳-湯姆遜節(jié)流效應(yīng)減弱，因此噴嘴溫降減小，導(dǎo)致出口溫度增加，射流密度略有增加，而射流速度稍有降低。在計(jì)算范圍內(nèi)，隨著圍壓從7.5 MPa增加至20.0 MPa，射流速度降低了2.5%，出口溫度和密度分別增加了2.40%和8.15%，而射流馬赫數(shù)減少了26.4%?？梢?jiàn)，與噴嘴壓降和入口溫度相比，圍壓對(duì)射流速度、出口溫度和密度的影響較小；而圍壓對(duì)射流馬赫數(shù)的影響顯著，隨著圍壓增加，射流馬赫數(shù)急劇減小，體現(xiàn)出射流特性向不可壓縮射流轉(zhuǎn)化的趨勢(shì)。

圖7 入口溫度320 K、噴嘴壓降40 MPa條件下射流特性隨圍壓的變化曲線Fig.7 Relation of jet speed and ambient pressure (inject temperature 320 K and nozzle pressure drop 40 MPa)

圖8 入口溫度320 K、噴嘴壓降40 MPa條件下射流物性隨圍壓的變化曲線Fig.8 Relation of jet properties and ambient pressure (inject temperature 320 K and nozzle pressure drop 40 MPa)

3 結(jié)論

1) 基于一維可壓縮流動(dòng)理論，考慮二氧化碳?xì)怏w真實(shí)性質(zhì)，建立了超臨界二氧化碳射流計(jì)算模型。該模型易于求解，計(jì)算結(jié)果與數(shù)值模擬結(jié)果吻合較好，精確度高，可適用于現(xiàn)場(chǎng)工藝的參數(shù)優(yōu)化和方案設(shè)計(jì)。

2) 超臨界二氧化碳射流特性分析結(jié)果表明，隨著噴嘴壓降從20 MPa增至50 MPa，射流速度和射流密度分別增加了49.1%和10.62%；當(dāng)入口溫度由280 K增至340 K，射流速度增加了10.9%，而射流密度降低了23.8%；隨著射流圍壓由7.5 MPa增至20.0 MPa，射流速度降低了2.5%，而射流密度增加了8.15%?？梢?jiàn)，噴嘴壓降對(duì)射流速度的影響顯著，而入口溫度對(duì)射流物性影響更大，射流圍壓對(duì)射流速度和物性影響均不明顯。因此，提高噴嘴壓降和圍壓都使射流的性質(zhì)向高密度不可壓縮流體轉(zhuǎn)化，而入口溫度的增加使射流性質(zhì)更趨向于低密度可壓縮流體。

[1] 徐純剛，李小森，蔡晶，等．二氧化碳置換法模擬開(kāi)采天然氣水合物的研究進(jìn)展[J]．化工學(xué)報(bào)，2013,64(7)：2309-2315．

XU Chungang,LI Xiaosen,CAI Jing,et al.Advance on simulation exploitation of natural gas hydrate by replacement with CO2[J].CIESC Journal,2013,64(7):2309-2315.

[2] 周守為，陳偉，李清平.深水淺層天然氣水合物固態(tài)流化綠色開(kāi)采技術(shù)[J].中國(guó)海上油氣，2014，26(5)：1-7.

ZHOU Shouwei,CHEN Wei,LI Qingping.The green solid fluiduzation development principle of natural gas hydrate stored in shallow layers of deep water[J].China Offshore Oil and Gas,2014,26(5):1-7.

[3] 周守為，李清平，陳偉，等.天然氣水合物開(kāi)采三維實(shí)驗(yàn)?zāi)M技術(shù)研究[J].中國(guó)海上油氣，2016，28(2)：1-9.DOI:10.11935/j.issn.1673-1506.2016.02.001.

ZHOU Shouwei,LI Qingping,CHEN Wei,et al.Research on 3D experiment technology of natural gas hydrate exploitation[J].China Offshore Oil and Gas,2016,28(2):1-9.DOI:10.11935/j.issn.1673-1506.2016.02.001.

[4] 周守為，陳偉，李清平，等.深水淺層非成巖天然氣水合物固態(tài)流化試采技術(shù)研究及進(jìn)展[J].中國(guó)海上油氣，2017，29(4)：1-8.DOI:10.11935/j.issn.1673-1506.2017.04.001.

ZHOU Shouwei,CHEN Wei,LI Qingping,et al.Research on the solid fluidization well testing and production for shallow non-diagenetic natural gas hydrate in deep water area [J].China Offshore Oil and Gas,2017,29(4):1-8.DOI:10.11935/j.issn.1673-1506.2017.04.001.

[5] 李根生，王海柱，沈忠厚，等．超臨界CO2射流在石油工程中應(yīng)用研究與前景展望[J]．中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2013,37(5)：76-80,87．

LI Gensheng,WANG Haizhu,SHEN Zhonghou,et al.Application investigations and prospects of supercritical Carbon dioxide jet in petroleum engineering[J].Journal of China University of Petroleum(Edition of Natural Science),2013,37(5):76-80,87.

[6] 李菊，劉妮，軒小波，等．二氧化碳水合物生成過(guò)程的強(qiáng)化方法研究進(jìn)展[J]．環(huán)境保護(hù)科學(xué)，2010,36(3)：14-18．

LI Ju,LIU Ni,XUAN Xiaobo,et al.Advance on methods for promoting Carbon dioxide hydrate formation[J].Environmental Protection Science,2010,36(3):14-18.

[7] KOLLE J J.Coiled-tubing drilling with supercritical carbon dioxide[R].SPE 65534,2000.

[8] 杜玉昆，王瑞和，倪紅堅(jiān)，等．超臨界二氧化碳射流破巖試驗(yàn)[J]．中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2012,36(4)：93-96．

DU Yukun,WANG Ruihe,NI Hongjian,et al.Rock-breaking experiment with supercritical Carbon dioxide jet[J].Journal of China University of Petroleum(Edition of Natural Science),2012,36(4):93-96.

[9] 黃飛，盧義玉，湯積仁，等．超臨界二氧化碳射流沖蝕頁(yè)巖試驗(yàn)研究[J]．巖石力學(xué)與工程學(xué)報(bào)，2015,34(4)：787-794．

HUANG Fei,LU Yiyu,TANG Jiren,et al.Research on erosion of shale impacted by supercritical Carbon dioxide jet[J].Chinese Journal of Rock Mechanics and Engineering,2015,34(4):787-794.

[10] 程宇雄，李根生，王海柱，等．超臨界CO2連續(xù)油管噴射壓裂可行性分析[J]．石油鉆采工藝，2013,35(6)：73-77．

CHENG Yuxiong,LI Gensheng,WANG Haizhu,et al.Feasibility analysis on coiled-tubing jet fracturing with supercritical CO2[J].Oil Drilling & Production Technology,2013,35(6):73-77.

[11] 程宇雄，李根生，王海柱，等．超臨界二氧化碳噴射壓裂孔內(nèi)流場(chǎng)特性[J]．中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014,38(4)：81-86．

CHENG Yuxiong,LI Gensheng,WANG Haizhu,et al.Flow field character in cavity during supercritical Carbon dioxide jet fracturing[J].Journal of China University of Petroleum(Edition of Natural Science),2014,38(4):81-86.

[12] 杜玉昆，王瑞和，陳曉紅，等.一種超臨界二氧化碳射流噴射鉆井開(kāi)采水合物裝置及方法:201510044567.3[P].2015-01-29.

[13] WANG Ruihe,HUO Hongjun,HUANG Zhiyuan,et al.Experimental and numerical simulations of Bottom hole temperature and pressure distributions of supercritical CO2jet for well-drilling[J].Journal of Hydrodynamics,2014,26(2):226-233.

[14] TIAN S C,HE Z G,LI G S,et al.Influences of ambient pressure and nozzle-to-target distance on SC-CO2jet impingement and perforation[J].Journal of Natural Gas Science and Engineering,2016,29:232-242.

[15] WANG Haizhu,LI Gensheng,TIAN Shouceng,et al.Flow field simulation of supercritical carbon dioxide jet:comparison and sensitivity analysis[J].Journal of Hydrodynamics,2015,27(2):210-215.

[16] 陳鴻雁，蔡建國(guó)，鄧修，等．超臨界溶液快速膨脹過(guò)程中噴嘴的流動(dòng)模型[J]．高?；瘜W(xué)工程學(xué)報(bào)，2000,14(3)：218-223．

CHEN Hongyan,CAI Jianguo,DENG Xiu,et al.A nozzle flow model of rapid expansion of supercritical solutions process[J].Journal of Chemical Engineering of Chinese Universities,2000,14(3):218-223.

[17] 童秉綱，孔祥言，鄧國(guó)華.氣體動(dòng)力學(xué)[M].北京:高等教育出版社，2012:58.

[18] WANG Zhiyuan,SUN Baojiang,WANG Jintang,et al.Experimental study on the friction coefficient of supercritical carbon dioxide in pipes[J].International Journal of Greenhouse Gas Control,2014,25:151-161.

[19] VESOVIC A,WAKEHAM W A.The transport properties of carbon dioxide[J].J.Phys.Chem.Ref.Data,1990,19(3):763-808.

[20] SPAN R,WAGNER W.A new equation of state for CO2covering the fluid region from the triple-point temperature to 1100 K at pressure up to 800 MPa[J].J.Phys.Chem.Ref.Data,1996,25(6):1509-1596.