国产日韩欧美一区二区三区三州_亚洲少妇熟女av_久久久久亚洲av国产精品_波多野结衣网站一区二区_亚洲欧美色片在线91_国产亚洲精品精品国产优播av_日本一区二区三区波多野结衣 _久久国产av不卡

?

海上CO2埋存井井筒溫度壓力影響因素研究

2023-05-30 20:10馬維臻管志川李成許玉強(qiáng)韓超萬飛鵬
石油機(jī)械 2023年2期
關(guān)鍵詞:井筒

馬維臻 管志川 李成 許玉強(qiáng) 韓超 萬飛鵬

摘要:準(zhǔn)確預(yù)測(cè)CO2埋存過程中井筒溫度壓力場(chǎng)以及CO2的物性參數(shù)變化對(duì)安全埋存至關(guān)重要。為此,建立了埋存井井筒溫度、壓力與CO2物性參數(shù)的耦合計(jì)算模型,計(jì)算得到了實(shí)例井井筒溫度壓力分布以及CO2物性參數(shù)隨井深的變化規(guī)律,并對(duì)注入?yún)?shù)對(duì)于井筒溫度壓力分布影響規(guī)律進(jìn)行分析。研究結(jié)果表明:井筒內(nèi)CO2流體的流速、努塞爾特?cái)?shù)和對(duì)流換熱系數(shù)隨井深的增加而增大,密度、黏度、摩阻系數(shù)、導(dǎo)熱系數(shù)和普朗特?cái)?shù)隨井深增加而減小,定壓比熱容在溫度壓力綜合作用下有一定波動(dòng);注入溫度對(duì)井筒壓力和井底溫度壓力影響很小;注入速率增大會(huì)使相同井深處溫度降低、壓力升高,調(diào)節(jié)注入速率可以在對(duì)壓力影響較小的同時(shí)有效調(diào)節(jié)井筒溫度分布;注入壓力的變化對(duì)壓力梯度幾乎無影響,在壓力較大時(shí)對(duì)溫度梯度影響較大,通過調(diào)節(jié)注入壓力可以有效調(diào)節(jié)井筒壓力分布。研究結(jié)果可為海上CO2埋存井井筒完整性的準(zhǔn)確評(píng)價(jià)提供理論基礎(chǔ)。

關(guān)鍵詞:海上CO2埋存;井筒;溫度壓力場(chǎng);物性參數(shù);耦合計(jì)算

0 引 言

“碳捕獲、利用與封存”(Carbon Capture Utilization and Storage,CCUS)是目前世界公認(rèn)最好的溫室氣體減排方案之一。CCUS技術(shù)把捕捉到的高純度CO2通過加壓至超臨界態(tài)注入到油氣儲(chǔ)層等地質(zhì)結(jié)構(gòu)中進(jìn)行封存。在海上CO2埋存過程中,一旦CO2發(fā)生泄露,會(huì)對(duì)海洋生態(tài)環(huán)境造成極其嚴(yán)重的影響。埋存井作為CO2泄露的重要途徑,其密封完整性對(duì)于安全埋存至關(guān)重要;而精確地計(jì)算CO2注入井井筒的溫度、壓力是進(jìn)行井筒密封完整性評(píng)價(jià)的重要基礎(chǔ)[1-2]。

雖然井筒溫度壓力場(chǎng)計(jì)算模型的研究已經(jīng)相對(duì)成熟[3-6],但由于CO2具有可壓縮性,其密度、黏度、比熱容等參數(shù)受溫度、壓力影響較大,所以在CO2注入地層的過程中其各項(xiàng)物性參數(shù)都在不斷變化[7-9]。而目前大多數(shù)研究中物性參數(shù)采用固定值或通過Span-Wagner模型等方法計(jì)算,這樣便存在計(jì)算誤差較大和計(jì)算量繁重等問題[10-12]。針對(duì)這些問題,筆者基于能量守恒方程、傳熱學(xué)等理論建立了適用于CO2注入井的井筒溫度壓力計(jì)算模型,并通過Matlab軟件直接調(diào)用Refprop9.0數(shù)據(jù)庫(kù)中各溫度、壓力下對(duì)應(yīng)的CO2物性參數(shù),實(shí)現(xiàn)了耦合計(jì)算,分析注入過程中井筒溫度壓力和CO2主要物性參數(shù)的變化規(guī)律以及注入?yún)?shù)對(duì)于井筒溫度壓力的影響。研究結(jié)果可為后續(xù)的井筒完整性評(píng)價(jià)等研究提供理論支撐。

馬維臻,等:海上CO2埋存井井筒溫度壓力影響因素研究

1 注入井溫度壓力計(jì)算模型

1.1 模型基本條件

溫度壓力模型的建立基于以下假設(shè):①井筒任意截面上各層結(jié)構(gòu)均無偏心現(xiàn)象,管柱密封條件良好,無泄漏;②CO2在井筒為一維流動(dòng),只考慮徑向傳熱,不考慮井深方向的軸向傳熱,同一截面上各點(diǎn)的溫度、壓力相等且流體物性不發(fā)生變化;③井筒內(nèi)的傳熱均為一維穩(wěn)態(tài)傳熱,水泥環(huán)外緣至地層傳熱為一維非穩(wěn)態(tài)傳熱;④不考慮套管接箍處對(duì)于傳熱的影響;⑤地層中的物理參數(shù)為常數(shù),不隨著溫度和深度的變化而變化。

1.2 環(huán)境溫度分布

以我國(guó)南海為例建立海水段環(huán)境溫度計(jì)算模型,我國(guó)南海位于北緯34°~11°55′,可通過南海表面月平均溫度統(tǒng)計(jì)數(shù)據(jù)選取海面溫度。該區(qū)域海水按垂直方向上的溫度變化可以分為混合層、躍溫層、恒溫層。水深小于200 m的混合層和躍溫層的溫度受海平面溫度影響較大,通過下式計(jì)算[13]:

3 CO2注入井井筒溫度壓力分布及流體參數(shù)變化規(guī)律

3.1 井筒溫度壓力分布

本文計(jì)算所用實(shí)例為南海某油田的淺層大位移CO2埋存井。該井為三開結(jié)構(gòu),井深為3 527 m,主要由油管、套管、水泥環(huán)、隔水管以及尾管構(gòu)成。套管底端通過封隔器坐封,油套環(huán)空中為完井液,CO2注入溫度為-10 ℃,注入壓力為20 MPa,注入速率為20 t/d。經(jīng)過計(jì)算可以得到該井井筒內(nèi)CO2流體溫度、壓力剖面圖,如圖1所示。實(shí)例井中CO2流體在1 255 m井深處由液態(tài)轉(zhuǎn)變?yōu)槌R界態(tài),井筒溫度和壓力均隨井深增大而升高,井筒內(nèi)外溫差減小使得溫度梯度有所減小,而井斜角的增大對(duì)溫度幾乎無影響,但是對(duì)壓力影響較大。這主要是因?yàn)樵搮^(qū)塊地溫梯度較小,使得井斜角發(fā)生改變時(shí)井筒內(nèi)外溫差的變化微乎其微,而井斜角的增大使得重力因素對(duì)壓力梯度的影響減小,摩擦因素影響增大[24],最終使得壓力梯度隨著井斜角的增大而降低,直到水平段重力的影響徹底消失,井筒壓力由于流體與管壁的摩擦有所下降。

3.2 CO2流體物性參數(shù)變化規(guī)律

3.2.1 密度和流速

井筒內(nèi)CO2流體密度和流速的變化曲線如圖2所示。由于液態(tài)和超臨界態(tài)下CO2主要因溫度升高而膨脹,所以可以看到隨著井深增加,密度逐漸減小,流速逐漸增大;同時(shí)在進(jìn)入尾管后由于管徑突然的增大使得流速迅速變緩。這里CO2流體密度以及流速的整體變化規(guī)律與孫寶江等[25]采用Span-Wagner模型計(jì)算的到的結(jié)果相符。

3.2.2 黏度和摩阻系數(shù)

圖3為CO2流體黏度和摩阻系數(shù)變化曲線。通過圖3可以看到,CO2黏度隨著井深的增加逐漸降低,這是由于在高壓下,流體黏度與溫度呈負(fù)相關(guān)所造成的[26]。而通過王志遠(yuǎn)所提出的計(jì)算公式(式(15))可知,CO2的摩阻系數(shù)受到密度、流速和黏度的共同影響,且密度、流速變化相對(duì)較小,所以摩阻系數(shù)在進(jìn)入尾管段發(fā)生突變前后都呈現(xiàn)與黏度相似的變化規(guī)律。

3.2.3 熱物性參數(shù)

CO2熱物性參數(shù)的變化規(guī)律如圖4所示。由圖4可以看到,導(dǎo)熱系數(shù)逐漸減小,由于定壓比熱與溫度呈正相關(guān)、與壓力呈負(fù)相關(guān)且主導(dǎo)因素不斷變化[26],所以變化規(guī)律相對(duì)復(fù)雜,但單一井段下比熱的變化規(guī)律同樣與孫寶江等[25]的計(jì)算結(jié)果相符。CO2流體的對(duì)流換熱系數(shù)是井筒溫度分布的重要影響因素,在油管尺寸不變的情況下,流體對(duì)流換熱系數(shù)為努塞爾特?cái)?shù)和導(dǎo)熱系數(shù)的函數(shù),且努塞爾特?cái)?shù)起主導(dǎo)作用,所以對(duì)流換熱系數(shù)與努塞爾特?cái)?shù)一樣在管徑變化的時(shí)候發(fā)生突變。

4 CO2注入?yún)?shù)對(duì)井筒溫度壓力分布的影響

4.1 不同CO2流體注入溫度下井筒溫壓分布

設(shè)置CO2流體注入壓力為20 MPa、注入速率為50 t/d,分別計(jì)算注入溫度為-35、-20、-10、0和10 ℃下的井筒溫度和壓力剖面,結(jié)果如圖5所示。通過圖5可以看到,CO2流體注入溫度的升高使得前期的溫度梯度和壓力梯度略微降低,但是對(duì)井底溫度、壓力的影響較小。上述注入溫度下的井底溫度分別為83.067、83.211、83.308、83.408和83.510 ℃,井底壓力分別為34.997、34.659、34.424、34.179和33.922 MPa??梢娬{(diào)節(jié)注入溫度可以在對(duì)井底溫度壓力影響較小的同時(shí)調(diào)節(jié)前期的井筒溫度。

4.2 不同CO2流體注入速率下井筒溫壓分布

設(shè)置CO2流體注入溫度為-20 ℃、注入壓力為20 MPa,分別計(jì)算注入速率分別為10、30、50、50和100 t/d時(shí)井筒的溫度和壓力剖面,結(jié)果如圖6所示。通過圖6可以看到,注入速率的增大使得溫度梯度明顯降低,而壓力梯度在斜井段略微升高。這主要是隨著注入速率的增加,流速相較于重力和摩擦阻力成為壓力梯度的主導(dǎo)因素。同時(shí),注入速率對(duì)于井底溫度的影響較為明顯。上述注入速率下的井底的溫度分別為88.933、87.466、83.211、77.402和68.163 ℃,而井底壓力分別為33.397、34.172、34.659、34.994和35.333 MPa。說明調(diào)節(jié)注入速率可以在對(duì)壓力影響較小的同時(shí)有效調(diào)節(jié)全井段的溫度分布以及井底溫度。

4.3 不同CO2流體注入壓力下井筒溫壓分布

設(shè)置CO2流體注入溫度為-20 ℃、注入速率為50 t/d,計(jì)算注入壓力分別5、10、20、30和40 MPa時(shí)井筒的溫度和壓力剖面,結(jié)果如圖7所示。通過圖7可以看到,注入壓力的增大使得溫度梯度在前期有一定增大。這是由于在CO2熱物性參數(shù)在低壓環(huán)境下受壓力影響較大[26],而壓力梯度幾乎不受注入壓力影響。注入壓力對(duì)于井底溫度的影響較為有限,對(duì)井底壓力影響較大。在上述注入壓力下,井底的溫度分別為77.789、80.867、83.211、84.152和84.658 ℃,井底壓力為18.499、23.967、34.659、45.178和55.599 MPa??梢娬{(diào)節(jié)注入壓力可在對(duì)井筒溫度不產(chǎn)生較大影響的同時(shí)有效調(diào)節(jié)井筒的壓力分布。

5.結(jié) 論

(1)綜合考慮溫度和壓力對(duì)于CO2流體物性參數(shù)的影響,基于傳熱學(xué)等原理并通過MATLAB調(diào)用Refprop9.0數(shù)據(jù)庫(kù)中CO2物性參數(shù),建立海上大位移CO2埋存井井筒溫度、壓力以及CO2物性參數(shù)耦合計(jì)算模型。以南海某埋存井為例,分析了井筒溫度、壓力分布和CO2流體物性參數(shù)隨井深的變化規(guī)律。

(2)在管徑相同的情況下,CO2流體的流速、努塞爾特?cái)?shù)和對(duì)流換熱系數(shù)隨井深的增加而增大,密度、黏度、摩阻系數(shù)、導(dǎo)熱系數(shù)和普朗特?cái)?shù)隨井深增加而減小,定壓比熱容在溫度壓力綜合作用下有一定波動(dòng)。

(3)注入溫度的升高使前期的溫度梯度和壓力梯度略微減?。蛔⑷胨俾实脑龃笫箿囟忍荻让黠@降低,而壓力梯度在斜井段略微升高;注入壓力的增大使溫度梯度在前期有一定增大,而壓力梯度幾乎不受注入壓力影響。

(4)調(diào)節(jié)注入溫度可以在對(duì)井底溫度壓力影響較小的同時(shí)調(diào)節(jié)前期的井筒溫度;調(diào)節(jié)注入速率可以在對(duì)壓力影響較小的同時(shí)有效調(diào)節(jié)全井段的溫度分布以及井底溫度;調(diào)節(jié)注入壓力可在對(duì)井筒溫度不產(chǎn)生較大影響的同時(shí)有效調(diào)節(jié)井筒的壓力分布。

參考文獻(xiàn):

[1] 胡永樂,郝明強(qiáng),陳國(guó)利,等.中國(guó)CO2驅(qū)油與埋存技術(shù)及實(shí)踐[J].石油勘探與開發(fā),2019,46(4):716-727.

HU Y L,HAO M Q,CHEN G L,et al.Technologies and practice of CO2 flooding and sequestration in China[J].Petroleum Exploration and Development,2019,46(4):716-727.

[2] 任韶然,李德祥,張亮,等.地質(zhì)封存過程中CO2泄漏途徑及風(fēng)險(xiǎn)分析[J].石油學(xué)報(bào),2014,35(3):591-601.

REN S R,LI D X,ZHANG L,et al.Leakage pathways and risk analysis of carbon dioxide in geological storage[J].Acta Petrolei Sinica,2014,35(3):591-601.

[3] 李成,管志川,張波,等.高溫高壓井油套環(huán)空早期圈閉壓力研究[J].石油機(jī)械,2020,48(9):43-50.

LI C,GUAN Z C,ZHANG B,et al.Study on trapped pressure in tubing-casing annulus at early production stage in HTHP wells[J].China Petroleum Machinery,2020,48(9):43-50.

[4] 張波.深水油氣井套管環(huán)空壓力預(yù)測(cè)及其調(diào)控機(jī)制[D].青島:中國(guó)石油大學(xué)(華東),2018.

ZHANG B.Prediction and mitigation mechanism of casing annular pressure in deepwater wells[D].Qingdao:China University of Petroleum (East China),2018.

[5] 宋洵成,管志川.深水鉆井井筒全瞬態(tài)傳熱特征[J].石油學(xué)報(bào),2011,32(4):704-708.

SONG X C,GUAN Z C.Full transient analysis of heat transfer during drilling fluid circulation in deep-water wells[J].Acta Petrolei Sinica,2011,32(4):704-708.

[6] 王軻,劉彪,張俊,等.高溫高壓氣井井筒溫度場(chǎng)計(jì)算與分析[J].石油機(jī)械,2019,47(1):8-13.

WANG K,LIU B,ZHANG J,et al.Calculation and analysis of wellbore temperature field in HTHP gas wells[J].China Petroleum Machinery,2019,47(1):8-13.

[7] RAYWARD-SMITH W J,WOODS A W.Some implications of cold CO2 injection into deep saline aquifers[J].Geophysical Research Letters,2011,38(6):1-6.

[8] 李圓圓.二氧化碳注入井筒內(nèi)相態(tài)分析[D].大慶:東北石油大學(xué),2019.

LI Y Y.Phase analysis in CO2 injection well[D].Daqing:Northeast Petroleum University,2019.

[9] 韓布興.超臨界流體科學(xué)與技術(shù)[M].北京:中國(guó)石化出版社,2005.

HAN B X.Supercritical fluid science & technology[M].Beijing:China Petrochemical Press,2005.

[10] 王瑞和,倪紅堅(jiān).二氧化碳連續(xù)管井筒流動(dòng)傳熱規(guī)律研究[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,37(5):65-70.

WANG R H,NI H J.Wellbore heat transfer law of carbon dioxide coiled tubing drilling[J].Journal of China University of Petroleum (Edition of Natural Science),2013,37(5):65-70.

[11] 竇亮彬,李根生,沈忠厚,等.注CO2井筒溫度壓力預(yù)測(cè)模型及影響因素研究[J].石油鉆探技術(shù),2013,41(1):76-81.

DOU L B,LI G S,SHEN Z H,et al.Wellbore pressure and temperature prediction model and its affecting factors for CO2 injection wells[J].Petroleum Drilling Techniques,2013,41(1):76-81.

[12] 蘭建平,龔群,徐治國(guó).CO2壓裂參數(shù)對(duì)井內(nèi)溫度和壓力的影響[J].石油機(jī)械,2018,46(11):97-103.

LAN J P,GONG Q,XU Z G.Effect of CO2 fracturing parameters on well temperature and pressure[J].China Petroleum Machinery,2018,46(11):97-103.

[13] 劉萍.深水氣井測(cè)試的溫度壓力分布研究[D].成都:西南石油大學(xué),2016.

LIU P.Study on temperature-pressure distribution of deepwater gas well testing[D].Chengdu:Southwest Petroleum University,2016.

[14] 高永海.深水油氣鉆探井筒多相流動(dòng)與井控的研究[D].東營(yíng):中國(guó)石油大學(xué),2007.

GAO Y H.Study on multi-phase flow in wellbore and well control in deep water drilling[D].Dongying:China University of Petroleum (East China),2007.

[15] RAMEY H J,Jr.Wellbore heat transmission[J].Journal of Petroleum Technology,1962,14(4):427-435.

[16] PAN L H,OLDENBURG C M,WU Y S,et al.Wellbore flow model for carbon dioxide and brine[J].Energy Procedia,2009,1(1):71-78.

[17] SOAVE G.Rigorous and simplified procedures for determining the pure-component parameters in the Redlich-Kwong-soave equation of state[J].Chemical Engineering Science,1980,35(8):1725-1730.

[18] SPAN R,WAGNER W.A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1100 K at pressures up to 800 MPa[J].Journal of Physical and Chemical Reference Data,1996,25(6):1509-1596.

[19] VESOVIC V,WAKEHAM W A,OLCHOWY G A,et al.The transport properties of carbon dioxide[J].Journal of Physical and Chemical Reference Data,1990,19(3):763-808.

[20] FENGHOUR A,WAKEHAM W A,VESOVIC V.The viscosity of carbon dioxide[J].Journal of Physical and Chemical Reference Data,1998,27(1):31-44.

[21] 王彌康.注蒸汽井井筒熱傳遞的定量計(jì)算[J].石油大學(xué)學(xué)報(bào)(自然科學(xué)版),1994(4):77-82.

WANG M K.Quantitative calculation of wellbore heat transmission for steam injection wells[J].Journal ofthe University of Petroleum,China (Natural Science Edition),1994(4):77-82.

[22] HASAN A R,KABIR C S.Aspects of wellbore heat transfer during two-phase flow[J].SPE Production & Facilities,1994,9(3):211-216.

[23] WANG Z Y,SUN B J,WANG J T,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.

[24] 張永剛,羅懿,劉岳龍,等.注CO2井筒溫度和壓力分布模型研究及現(xiàn)場(chǎng)應(yīng)用[J].巖性油氣藏,2014,26(2):108-113.

ZHANG Y G,LUO Y,LIU Y L,et al.Research and application of wellbore temperature and pressure distribution models for CO2 injection well[J].Lithologic Reservoirs,2014,26(2):108-113.

[25] 孫寶江,孫小輝,王志遠(yuǎn),等.超臨界CO2鉆井井筒內(nèi)流動(dòng)參數(shù)變化規(guī)律[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,40(3):88-95.

SUN B J,SUN X H,WANG Z Y,et al.Flow behavior analysis during supercritical CO2 drilling:consideration of varying thermo-physical properties of CO2 in wellbore[J].Journal of China University of Petroleum (Edition of Natural Science),2016,40(3):88-95.

[26] 呂欣潤(rùn).頁(yè)巖氣水平井CO2壓裂井筒溫度場(chǎng)研究[D].北京:中國(guó)石油大學(xué)(北京),2019.

LYU X R.Study on the temperature field of horizontal wellbores during CO2 fracturing in shale gas reservoirs[D].Beijing:China University of Petroleum (Beijing),2019.

猜你喜歡
井筒
基于DDA 方法的排碴井筒堵塞概率分析
再談滿井放礦
雙機(jī)抬吊法吊裝某井筒的過程靜力學(xué)分析
金川礦區(qū)龍首礦西二副井井筒加固
礦井井筒煤柱開采技術(shù)措施
滑溜水壓裂支撐劑在水平井筒內(nèi)沉降規(guī)律研究
煤峪口礦西三井筒提升中心的測(cè)定
復(fù)雜地段副斜井井筒施工方法的選擇
煤礦井筒施工技術(shù)的探討
河南省新豐煤礦主井井筒中心確定方法探討