李　東， 吳玉國， 潘　振， 段騰龍

(遼寧石油化工大學(xué) 石油天然氣工程學(xué)院，遼寧 撫順 113001)

?

氣舉法消除嚴(yán)重段塞流的模擬分析研究

李東， 吳玉國， 潘振， 段騰龍

(遼寧石油化工大學(xué) 石油天然氣工程學(xué)院，遼寧 撫順 113001)

在海洋油氣資源開采中，特殊的管路系統(tǒng)、重力、流量等因素會(huì)形成嚴(yán)重段塞流。嚴(yán)重段塞流會(huì)引發(fā)次生危害，威脅安全生產(chǎn)。因而需要采取方法消除嚴(yán)重段塞流。采用OLGA軟件和PVTsim軟件，對氣舉法消除嚴(yán)重段塞流進(jìn)行模擬研究。研究結(jié)果表明，氣舉法在一定程度上可以消除嚴(yán)重段塞流,注氣位置不同，效果也有一定的差別,注氣位置在立管底部的上游位置更能有效地消除嚴(yán)重段塞流。

氣舉法；嚴(yán)重段塞流；OLGA軟件；模擬

海底的油氣采出后，從井口傳輸?shù)胶Ｑ笃脚_(tái)，由于地形等原因，管道總會(huì)經(jīng)過水平—下傾—立管形成L形立管。此特殊的L形立管在油田生產(chǎn)的早期和后期時(shí)，油氣流量較低，在重力等的作用下，使得下傾管內(nèi)形成分層流，之后，液體會(huì)堆積在立管底部堵塞氣體的流出，從而形成液塞，并且隨著時(shí)間的推移液塞長度會(huì)不斷增加，達(dá)到一個(gè)甚至幾個(gè)立管高度，形成嚴(yán)重段塞流。嚴(yán)重段塞流會(huì)使管內(nèi)的壓力和流量周期性的劇烈波動(dòng)，進(jìn)而導(dǎo)致管線的接頭等構(gòu)件受到應(yīng)力腐蝕。化學(xué)腐蝕和應(yīng)力腐蝕之間的協(xié)同作用，會(huì)加劇管道腐蝕。同時(shí)嚴(yán)重段塞流引發(fā)管路的背壓升高，使得油氣產(chǎn)量降低。在壓力劇烈降低時(shí)，會(huì)有蠟晶析出，堵塞管道；在液塞噴發(fā)階段，會(huì)出現(xiàn)溢流現(xiàn)象；在氣液交替階段，會(huì)出現(xiàn)斷流、汽蝕現(xiàn)象，損害設(shè)備，甚至導(dǎo)致下游設(shè)備停產(chǎn)。因此對消除嚴(yán)重段塞流的方法進(jìn)行深入研究具有重要意義。

在消除嚴(yán)重段塞流的方法研究，早在1973年B.T.Yocum[1]提出了增加回壓法來消除嚴(yán)重段塞流。1977年，Z.Schmidt[2]提出了節(jié)流法來消除嚴(yán)重段塞流。之后一些學(xué)者利用反饋法，即在立管底部和管道出口處加裝一差壓變送器來控制節(jié)流閥的開度來消除嚴(yán)重段塞流。1985年，B.M.Pots[3]提出了氣舉法。1990年，Kaasa[4]提出了多加一根立管的氣液分離法。1994年，M.W.J.Wyllie等[5]提出在立管內(nèi)加裝套管的方法來消除嚴(yán)重段塞流。1997年，K.S.Johal[6]提出了加裝多相泵的方法來消除嚴(yán)重段塞流，并提出了將其他管線的氣體引入管路來實(shí)現(xiàn)氣舉的方法。1999年，A.R.Almeida等[7]提出了在上升管底部使用文丘里管的方法來消除嚴(yán)重段塞流。2000年，S.Song等[8]提出了海底預(yù)分離法。2000年，C.Sarica等[9]提出了自供氣舉法。2003年K.Kovalev等[10]報(bào)道了殼牌公司在段塞流捕集器前加裝小型預(yù)分離器的S3(Slug Suppresion System)法。2006年，羅曉明等[11]提出在下傾管底部加裝節(jié)流孔板的方法來消除嚴(yán)重段塞流。2011年，呂宇玲等[12]提出在立管頂端加裝潛油泵的方法來消除嚴(yán)重段塞流。2015年，邱偉偉等[13]提出了分流法。增加回壓法、節(jié)流法、氣液分離法都會(huì)導(dǎo)致產(chǎn)量降低。加裝多相泵方法中，泵的選型、具體安裝位置為此方法應(yīng)用帶來局限。泡沫法中，選用合適的泡沫劑和怎么形成泡沫是難題。海底預(yù)分離法中，海底情況復(fù)雜多變，預(yù)分離器的安裝，運(yùn)行維護(hù)是個(gè)難題。氣舉法與之相比具有很大優(yōu)勢，其原理是，在立管上游管段或者立管中部的合適位置注入壓縮氣體，來增大立管內(nèi)氣體的壓力和液體的流速，增大立管內(nèi)的真實(shí)含氣率，從而達(dá)到減小壓力波動(dòng)的幅度、系統(tǒng)穩(wěn)定的目的。

為了進(jìn)一步探究嚴(yán)重段塞流消除方法的效果及其影響因素，本文采用OLGA軟件對其中的氣舉法進(jìn)行模擬研究，并分析注氣位置對氣舉法消除嚴(yán)重段塞流效果的影響。

1　OLGA軟件介紹

OLGA軟件是由IFE和SINTEF聯(lián)合在幾個(gè)石油公司的合作下開發(fā)的，其最早版本成型于1984年，之后經(jīng)過一系列的修改，形成雙流體模型。它包含的模塊有水三相模塊、段塞流模塊、組分跟蹤模塊、三維熱傳導(dǎo)模塊、水合物抑制劑追蹤模塊等，以及MATLAB工具箱。OLGA軟件[14]計(jì)算中的閉合定律需要以流型為基礎(chǔ)。模型基于最小滑移概念將流型劃分為分離流和分散流。分離流又可劃分為環(huán)狀-彌散流和分層流，分散流可劃分為段塞流和氣泡流。流型之間靠最小滑移理論和最小準(zhǔn)則來確定其是否轉(zhuǎn)換。分層流與環(huán)狀流的轉(zhuǎn)換是以濕周為參數(shù)，當(dāng)液膜與內(nèi)圓周相同時(shí)，流型轉(zhuǎn)換為分層流。當(dāng)液塞空隙率比段塞平均空隙率大時(shí)，液塞會(huì)產(chǎn)生連續(xù)的氣泡流。當(dāng)氣液兩相中，波浪高度到達(dá)管內(nèi)頂部，層狀流轉(zhuǎn)變?yōu)榄h(huán)狀流。

在模型求解過程中，OLGA模型求解采用有限容積法，在此過程中，管線離散成許多控制體單元，壓力、溫度等變量儲(chǔ)存在控制體的中心，流量、流速等存儲(chǔ)在控制體邊界。對于瞬態(tài)多相流問題，模型中采用拉格朗日網(wǎng)格，可由液塞前段位置修正控制體平衡方程，減小偏離真實(shí)解的誤差。之后，開發(fā)者依據(jù)SINTEF兩相流環(huán)道實(shí)驗(yàn)的數(shù)據(jù)對OLGA軟件進(jìn)行了改進(jìn)。K.H.Bendiksen等[15]將OLGA軟件的預(yù)測結(jié)果與SINTEF數(shù)據(jù)進(jìn)行對比，結(jié)果非常好。H.Dhulesia等[16]利用Boussens和SINTEF環(huán)道實(shí)驗(yàn)數(shù)據(jù)，將OLGA軟件和PEPITE、TACITE等軟件在壓降和持液率等方面做了全面對比，發(fā)現(xiàn)OLGA軟件的預(yù)測效果比較好。

2　數(shù)值模擬研究

2.1模型描述

在國內(nèi)外眾多學(xué)者研究的基礎(chǔ)之上，對典型的L型立管進(jìn)行模擬。管道的流徑只有一條，在出現(xiàn)嚴(yán)重段塞流時(shí)，流徑中L型管道的水平段和下傾段、下傾段底部、垂直管段和垂直管段下游的水平管段的各項(xiàng)參數(shù)具有很大的差異，據(jù)此將流徑劃分為5條管道，依次為PIPE-1、PIPE-2、PIPE-3、PIPE-4、 PIPE-5。OLGA的空間積分是由專門的網(wǎng)格劃分組件劃分網(wǎng)格后實(shí)現(xiàn)的。為了確保準(zhǔn)確性，每一管線至少劃分兩節(jié)以上，并且相鄰管線比值介于0.5到2，此舉對于管段模擬非常重要。鑒于此以及管道的長度，管道劃分為如下幾段。管道1和2全部劃分為10段，管道3劃分為2段，管道4劃分為5段，管道5劃分為3段，立管長度的具體數(shù)據(jù)如圖1所示，效果圖如圖2所示。

圖1　管道地形編輯器界面圖

Fig.1Pipeline geometry editor interface figure

圖2　管道立體效果圖

Fig.2Pipeline perspective renderings

為了方便模擬，對于管道中流體組分所采用的是軟件自帶模型中的水動(dòng)力段塞流的組分，如表1所示。注氣時(shí)所用的氣源為OLGA軟件中的一個(gè)物性源替代，其物性源的文件可由組分文件來取代，其組分為空氣。模擬研究過程中，假設(shè)管壁是絕熱的。節(jié)點(diǎn)NODE2的壓力PRESSURE 為65×105Pa，節(jié)點(diǎn)溫度TEMPERATURE為20 ℃，是質(zhì)量型。管道中的流體SOURCE_1為質(zhì)量型，MASSFLOW為2 kg/s。對于注氣所用的氣源SOURCE_2的類型SOURCETYPE也為質(zhì)量型。注氣法的兩個(gè)關(guān)鍵注氣位置點(diǎn)，依次是PIPE-3的第2部分和PIPE-4的第1部分。閃蒸模型為WATER、流體相態(tài)為三相，并且相間無滑移。

表1　油氣組分

2.2模擬結(jié)果及分析

注氣的位置在管道的PIPE-3的第2部分(428,-19.947)處，注氣的流量為0.5、5、50、500 m3/h時(shí),分別運(yùn)行模擬。管道的PIPE-4第3段處的持液率隨時(shí)間的變化曲線如圖3所示。注氣位置在PIPE-4的第1部分(430,-18.122)處時(shí)，注氣的流量為0.5、5、50、500 m3/h時(shí)，管道的PIPE-4第3段處的持液率隨時(shí)間的變化曲線如圖4所示。

圖3　注氣位置在3管2段時(shí)觀測點(diǎn)持液率隨時(shí)間變化

Fig.3Versus time curve of the liquid holdup gas injection position in PIPE-3 second section at the observation point

圖4　注氣位置在4管1段時(shí)觀測點(diǎn)持液率隨時(shí)間變化曲線

Fig.4Versus time curve of the liquid holdup gas injection position in PIPE-4 first section at the observation point

注氣位置在管道PIPE-3第2段的時(shí)候，注氣所用流量變化時(shí)，對圖3所示持液率的數(shù)據(jù)變化曲線分析可知，管道PIPE-4第3段處，持液率依然呈周期性變化。對注氣量增大到500 m3/h后，嚴(yán)重段塞流的現(xiàn)象依然存在，并且對數(shù)據(jù)具體分析后由表2可以看到，段塞流的一個(gè)周期時(shí)間變短，時(shí)間變短表明可以在一定程度上消除嚴(yán)重段塞流，但是仍然不能完全消除段塞流。

表2　模擬結(jié)果關(guān)鍵參數(shù)值對比

當(dāng)注入的氣體流量相同時(shí)，對于注氣位置在管道PIPE-3第2段和注氣位置在管道PIPE-4第1段時(shí)，對于同一氣源流量，在相同的情況下，注氣位置選在管道PIPE-3第2段的時(shí)候，持液率的周期時(shí)間變的更短，在總過流斷面面積中液相的過流斷面面積占總過流斷面面積的周期時(shí)間變短。在液塞形成階段，持液率從0到1的時(shí)間變短，液塞由長液塞變?yōu)槎桃喝布匆喝拈L度變短，嚴(yán)重段塞流的流型改變[17]。由此表明，注氣的位置選在管道PIPE-3第2段比注氣位置在管道PIPE-4第1段能更好的消除嚴(yán)重段塞流。

注氣量為0.5、5、50、500 m3/h的情況下，運(yùn)行模擬后，得到如圖5、6所示的曲線。圖5為在管道PIPE-3第2段(428,-19.947)處注氣時(shí)，管道PIPE-4第3段處的壓力隨時(shí)間變化曲線；圖6為注氣位置在管道PIPE-4第1段(430,-18.122)處時(shí)，管道PIPE-4第3段處的壓力隨時(shí)間變化曲線圖。

從圖5和圖6中可以看到，壓力的最大峰值隨著注氣量的增大呈現(xiàn)出變小的總體趨勢，但是當(dāng)注氣量為5 m3/h時(shí)，壓力的最大峰值不僅未變小，反而變大，這會(huì)導(dǎo)致壓力波動(dòng)更大，此注氣量不利于消除嚴(yán)重段塞流。在注氣量為50 m3/h和注氣量為500 m3/h穩(wěn)定狀態(tài)時(shí)，平均一個(gè)液塞周期內(nèi)，注氣位置在PIPE-3管第2段相對于注氣位置在PIPE-4管第1段相同狀況下，在液塞形成階段所用時(shí)間短；在液塞發(fā)展階段所用時(shí)間長；在液塞流出階段所形成的壓力最高峰值更低，表明在PIPE-3管第2段注氣時(shí)管道內(nèi)的壓力變化更平穩(wěn)，此結(jié)論在表2中也可看出。在一個(gè)段塞周期內(nèi)壓力波動(dòng)變小，周期時(shí)間變短，但是仍然不能消除嚴(yán)重段塞流。從圖5和圖6中還可以得到，在同一注氣流量下，圖5中的一個(gè)段塞的周期時(shí)間會(huì)變更短，這表明在一定程度上減弱了嚴(yán)重段塞流。注氣位置在管道PIPE-3第2段的位置比管道PIPE-4第1段的位置注氣效果更好。此結(jié)論與上文對持液率的分析一致，表明結(jié)果可靠。

圖5　注氣位置在PIPE-3管第2段3管2段時(shí)觀測點(diǎn)壓力隨時(shí)間變化曲線

Fig. 5Versus time curve of the pressure gas injection position in PIPE-3 second section at the observation point

圖6　注氣位置在PIPE-4管第1段時(shí)觀測點(diǎn)壓力隨時(shí)間變化曲線

Fig.6Versus time curve of the pressure gas injection position in PIPE-4 first section at the observation point

3　結(jié)論

利用OLGA軟件對氣舉法消除嚴(yán)重段塞流的效果進(jìn)行了模擬研究，并對結(jié)果進(jìn)行了分析。模擬研究結(jié)果表明，注氣的位置在立管底部的上游比立管底部的下游消除嚴(yán)重段塞流的效果更好；注氣量增大時(shí)，在一定程度上對嚴(yán)重段塞流有一定的消除作用，但并不能完全消除嚴(yán)重段塞流,在注氣量增大到使得立管內(nèi)流型由嚴(yán)重段塞流轉(zhuǎn)化為其他流型時(shí)，可消除嚴(yán)重段塞流，此注氣量需要通過流型圖判斷和計(jì)算得到。此結(jié)論表明，在實(shí)際生產(chǎn)中，如果周邊有合適的氣源，可以采用氣舉法來消除嚴(yán)重段塞流，可以明顯增加油氣輸送量，保障正常安全生產(chǎn)。除此之外，本文的研究并沒有在折算氣速和折算液速的角度分析，在之后的研究中，需要進(jìn)一步分析。

[1]Yocum B T．Offshore riser slug flow avoidance[R].SPE 4312, 1973.

[2]Schmidt Z．Experimental study of two-phase slug flow in a pipeline—riser pipe system[D]．Tulsa: University of Tulsa,1977．

[3]Pots B M,Bromilow I G,Konijn M J .Severe slug flow in offshore flow-line/riser systems[J].SPE J., 1987,2(4):319-324．

[4]Kaasa ?.A subsea slug catcher to prevent severe slugging[C]//6th Underwater Tech .Int .Conf. Bergen,Norway:[s.n.],1990．

[5]Wyllie M W J,Brackenridge A A. Retrofit solution to reduce slugging effects in multiphase subsea pipelines.The internal riser insert system(IRIS)[C].[S.L.]：Subsea International Conference，1994．

[6]Johal K S. An alternative economic method to riser-base gas lift for deep water subsea oil/gas field developments [R]. SPE 38541,1997.

[7]Almeida A R，Goncalves M A L.Venture for severe slugging elimination[C].[S.L.]：BHR Group 1999 Multiphase,1999：149-158.

[8]Song S, Kouba G.Fluids transport optimization using seabed separation[J].Journal of Energy Resources Technology,2000,122(3):105-109.

[9]Sarica C，Tengesdal J.A new technique to eliminate severe slugging in pipeline/riser systems[R]. SPE 63185,2000.

[10]Kovalev K,Cruickshank A, Purvis J.The slug suppression system in operation[R].SPE 84947,2003.

[11]羅小明,何利民,趙越超,等.上升管中強(qiáng)烈段塞流消除方法試驗(yàn)研究[J].工程熱物理學(xué)報(bào), 2006,27(3):445-447.

Luo Xiaoming，He Limin，Zhao Yuechao，et al. Experimental study on the elimination of severe slug ging in riser．pipes[J]. Journal of Engineering Thermophysics，2006,27(3):445-447.

[12]呂宇玲，何利民，牛殿國,等.自控泵吸法消除強(qiáng)烈段塞流[J].化工學(xué)報(bào),2011,62(7): 1846-1851.

Lv Yuling，He Limin， Niu Dianguo，et al. Auto-pumping method for severe slugging elimination[J]. CIESC Journal，2011,62(7): 1846-1851.

[13]邱偉偉,宮敬.分流法消除嚴(yán)重段塞流實(shí)驗(yàn)研究[J].科學(xué)技術(shù)與工程, 2015,15(1):135-138.

Qiu Weiwei， Gong Jing. Elimination experiment of severe slugging in the riser[J]. Science Technology and Engineering，2015,15(1):135-138.

[14]劉新野.海底混輸管線嚴(yán)重段塞流的預(yù)測與控制[D].東營：中國石油大學(xué)(華東),2011.

[15]Bendiksen K H,Moe R,Nuland S. The dynamic two-fluid model OLGA: theory and application[J].SPE Prpduction Engineering, 1991,6(2):171-180.

[16]Dhulesia H,Lopez D.Critical evaluation of mechanistic two-phase flow pipeline and well simulation models[R].SPE 36611, 1996.

[17]常方圓，馬貴陽，安利姣，等. 傾斜管內(nèi)多相流流型研究[J]. 遼寧石油化工大學(xué)學(xué)報(bào),2016,36(1):34-36.

Chang Fangyuan，Ma Guiyang，An Lijiao，et al. The multiphase flow in inclined pipe[J]. Journal of Liaoning Shihua University，2016, 36(1):34-36.

(編輯王亞新)

Simulation Analysis of the Gas-Lift to Eliminate Severe Slugging

Li Dong, Wu Yuguo, Pan Zhen, Duan Tenglong

(SchoolofPetroleumandGasEngineering,LiaoningShihuaUniversity,FushunLiaoning113001,China)

In the exploitation of offshore oil and gas resources, the severe slugging will appear due to the special piping system, gravity, flow, and other factors. The severe slugging will cause secondary hazards and threaten the safety of production. Therefore, it is needed to adopt the method to eliminate the severe slugging. In this paper, OLGA software and PVTsim software were used to simulate the gas-lift to eliminate the severe slugging. The results showed that the severe slugging could be eliminate by gas-lift in a certain extent. Different gas injection position had different effect to eliminate severe slugging. When the position of gas injection was in the upstream at the bottom of the riser, severe slugging could be eliminated effectively.

The gas-lift method; Severe slugging; OLGA software; Simulation

1006-396X(2016)03-0087-05

2016-01-05

2016-05-30

遼寧省教育廳科學(xué)研究項(xiàng)目(L2015306)；遼寧省高等學(xué)校優(yōu)秀人才支持計(jì)劃資助(LJQ2014038)。

李東(1990-)，男，碩士研究生，從事油氣管道流動(dòng)保障研究；E-mail：447674340@qq.com。

吳玉國(1977-)，男，博士，副教授，從事油氣儲(chǔ)運(yùn)方面研究；E-mail：wyg0413@126.com。

TE53

Adoi:10.3969/j.issn.1006-396X.2016.03.017

投稿網(wǎng)址：http://journal.lnpu.edu.cn