• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看

      ?

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

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

      馬維臻 管志川 李成 許玉強 韓超 萬飛鵬

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

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

      0 引 言

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

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

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

      1 注入井溫度壓力計算模型

      1.1 模型基本條件

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

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

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

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

      3.1 井筒溫度壓力分布

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

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

      3.2.1 密度和流速

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

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

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

      3.2.3 熱物性參數(shù)

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

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

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

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

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

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

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

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

      5.結(jié) 論

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

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

      (3)注入溫度的升高使前期的溫度梯度和壓力梯度略微減小;注入速率的增大使溫度梯度明顯降低,而壓力梯度在斜井段略微升高;注入壓力的增大使溫度梯度在前期有一定增大,而壓力梯度幾乎不受注入壓力影響。

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

      參考文獻:

      [1] 胡永樂,郝明強,陳國利,等.中國CO2驅(qū)油與埋存技術(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泄漏途徑及風險分析[J].石油學報,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].石油機械,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)空壓力預測及其調(diào)控機制[D].青島:中國石油大學(華東),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].石油學報,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] 王軻,劉彪,張俊,等.高溫高壓氣井井筒溫度場計算與分析[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].大慶:東北石油大學,2019.

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

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

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

      [10] 王瑞和,倪紅堅.二氧化碳連續(xù)管井筒流動傳熱規(guī)律研究[J].中國石油大學學報(自然科學版),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井筒溫度壓力預測模型及影響因素研究[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] 蘭建平,龔群,徐治國.CO2壓裂參數(shù)對井內(nèi)溫度和壓力的影響[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] 劉萍.深水氣井測試的溫度壓力分布研究[D].成都:西南石油大學,2016.

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

      [14] 高永海.深水油氣鉆探井筒多相流動與井控的研究[D].東營:中國石油大學,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].石油大學學報(自然科學版),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)場應用[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] 孫寶江,孫小輝,王志遠,等.超臨界CO2鉆井井筒內(nèi)流動參數(shù)變化規(guī)律[J].中國石油大學學報(自然科學版),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] 呂欣潤.頁巖氣水平井CO2壓裂井筒溫度場研究[D].北京:中國石油大學(北京),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 方法的排碴井筒堵塞概率分析
      再談滿井放礦
      雙機抬吊法吊裝某井筒的過程靜力學分析
      金川礦區(qū)龍首礦西二副井井筒加固
      金屬礦山(2021年6期)2021-07-10 06:08:38
      礦井井筒煤柱開采技術(shù)措施
      黑龍江科學(2016年1期)2016-03-15 08:22:30
      滑溜水壓裂支撐劑在水平井筒內(nèi)沉降規(guī)律研究
      煤峪口礦西三井筒提升中心的測定
      復雜地段副斜井井筒施工方法的選擇
      人間(2015年21期)2015-03-11 15:24:48
      煤礦井筒施工技術(shù)的探討
      河南科技(2014年18期)2014-02-27 14:14:46
      河南省新豐煤礦主井井筒中心確定方法探討
      河南科技(2014年16期)2014-02-27 14:13:07
      淳化县| 治县。| 宁强县| 巴南区| 教育| 曲阜市| 赤城县| 防城港市| 红桥区| 张家港市| 镇坪县| 拜城县| 铜梁县| 手游| 都匀市| 新巴尔虎左旗| 九台市| 高青县| 清水县| 甘德县| 贡觉县| 安新县| 龙江县| 长治市| 四子王旗| 林州市| 桂东县| 广水市| 通城县| 呼图壁县| 永清县| 石嘴山市| 德昌县| 天门市| 南昌县| 陇川县| 闽清县| 山阳县| 松滋市| 仙游县| 荆州市|