CO2井筒相態(tài)分布3學習資料

1、Good is good, but better carries it.精益求精，善益求善。CO2井筒相態(tài)分布3-CO2井筒壓力溫度的分布張勇唐人選（中石化華東分公司采油廠，江蘇泰洲，225300）摘要：本文根據(jù)垂直管流的能量平衡方程，結合Ramy井筒溫度分布計算方法，推導出CO2在不同的注入速度、注入溫度、注入時間等因素下的井筒壓力、溫度分布。文章以蘇北工區(qū)草8井試注實測數(shù)據(jù)為例加以驗證，其中壓力誤差3.8%，溫度誤差2.5%。在此基礎上進一步討論了影響CO2井底溫度的因素，并預測了3口注入井井筒壓力溫度分布。關鍵題：二氧化碳井筒壓力井筒溫度預測前言蘇北工區(qū)目前正在進行中石化草舍油田泰州組油

2、藏CO2驅(qū)提高采收率先導試驗項目，其中CO2注入井井底的壓力和溫度是大家關心主要問題之一。因為其大小直接決定了油田是否混相，如注入井井底附近均不能混相，那么其它地方更不能混相，因為注入井井底壓力要高于油藏的壓力。另外，實測壓力溫度也很費工費時，文章基于地層熱力學計算及CO2物性數(shù)據(jù)，編制了相應的預測軟件，預測井底壓力和溫度，實例表明該方法是可行的。1、注入井壓力計算方法根據(jù)能量平衡方程，注入井單位管長的壓降損失可寫為（見文獻1）：（1）其中：為流體的重力；為加速度損失；為摩擦損失?？紤]加速度損失主要取決于速度的改變，對于可壓縮流體而言主要取決于流體密度的改變，從實際數(shù)據(jù)分析，CO2在井筒中的密

3、度變化不是很大，因此該項損失可忽略不計。為井筒軸線與水平方向的夾角，實際計算可用水平位移與相應的垂深來確定。在上式中有關參數(shù)的確定：1.1、注入速度的計算（2）1.2、密度的計算為井筒中CO2的密度，與溫度壓力值有關。根據(jù)真實氣體壓縮方程，考慮氣體偏差系數(shù)，則在不同壓力溫度下的氣體密度為（3）1.3、氣體偏差系數(shù)的計算氣體的偏差系數(shù)可采用Redlich-Kwong方程求得，Redlich-Kwong方程如下：（4）給定實際壓力和溫度，求解方程（4）即得不同壓力溫度下的偏差系數(shù)（圖1）。將偏差系數(shù)代入（3），即可求得不同壓力溫度下的井筒密度（圖2）。圖1不同壓力溫度下CO2的壓縮系數(shù)圖2不同壓力

4、溫度下CO2的密度4、CO2在井筒中溫度的計算在CO2注入過程中，由于CO2溫度與周圍井筒溫度不同而存在溫差，因此注入的CO2必然與周圍存在熱交換，導致注入的CO2的溫度發(fā)生改變，Ramy在文獻2中推導出流體溫度沿井筒的變化關系式，該關系式考慮了注入時間、注入速度等因素。（5）其中：（6）（7）由于油管和套管的熱阻很小，從注入流體到地層總的傳熱系數(shù)Ut可表示為（8）Ut的求法見文獻3，是經(jīng)過反復迭代求得。另外，在求Ut的過程中需要用到若干參數(shù)，如環(huán)空是充滿防腐液的水，需求不同溫度下的水的密度、水的粘度、水的熱膨脹系數(shù)、水的熱容及水的導熱系數(shù)等參數(shù)，這些參數(shù)均是壓力溫度的函數(shù)，其參數(shù)大小參見文獻

5、3。求出Ut后，再求出A，代入公式（5）即可求得某一深度下的溫度T。在采用Ramy法計算井筒中流體的溫度時，需要對以下參數(shù)進行確定：CO2的熱容C文獻4列出了不同壓力溫度下CO2的熱容值，經(jīng)分析CO2的熱容隨溫度的升高而降低，隨壓力的升高先升后降。實際應用時，可將其簡化，在壓力P10MPa，溫度0t100條件下，CO2的熱容變化范圍在0.8181.38kJ/(kg.K)，取其平均值1.09kJ/(kg.K)；在P10MPa時，溫度0t100條件下，CO2的熱容變化范圍在1.4272.66kJ/(kg.K)，取其平均值2.04kJ/(kg.K)。地層的導熱系數(shù)與地層熱擴散系數(shù)文獻3列出了美國熱采

6、工程常用的油藏巖石熱物性參數(shù)，典型的孔隙性巖石的導熱系數(shù)K=1.72.4(W/(m.K)，取其平均值2.09(W/(m.K)，地層熱擴散系數(shù)約為0.0037m2/h。水泥漿導熱系數(shù)和水泥環(huán)導熱系數(shù)，文獻3列出了泥漿的導熱系數(shù)0.60W/(m.K)，水泥環(huán)的導熱系數(shù)0.35W/(m.K)。以上這些參數(shù)已為美國某些油公司廣泛采用，并列入某些熱采工程手冊中。1.5摩阻系數(shù)的確定摩阻系數(shù)的大小一般根據(jù)雷諾數(shù)由經(jīng)驗公式確定，不同粗糙管的摩阻系數(shù)的經(jīng)驗公式見文獻5。在計算井筒中CO2雷諾數(shù)的大小時，需要用到CO2的粘度，其大小同樣與壓力溫度有關，文獻4列出了不同壓力溫度下CO2的粘度（圖3）。圖3不同壓力

7、溫度下CO2的粘度在上述參數(shù)均確定后，即可求得在不同注入速度、注入溫度、注入壓力及注入時間下井筒中壓力溫度的分布。2、實際應用分析2.1壓力、溫度計算值與實測值的對比實例：草舍油田草8井2005年8月24日進行CO2試注試驗。用MPS97電子壓力計于24日9為22下至3100m，所測靜壓40.18MPa/3100m（井筒水壓力），靜溫108/3100m（井底溫度）。試驗具體過程及數(shù)據(jù)見下表（表1）。由于總體注入時間較短，四個測試點又在同一深度測試，計算結果可能重復，下面以第3個測試點為例。表1草8井注入?yún)?shù)及注入過程對表1數(shù)據(jù)整理，第4個測試點累計注入13個小時（t=0.54d），累積注入量為

8、30t，折算平均每天注入量為Gt=55.4t/d。井口注入壓力24.5MPa，注入溫度-20。該井的水泥返高2074m。井眼結構參數(shù)如下：井眼直徑rh=215.9mm，套管外徑rco=137mm，套管內(nèi)徑rci=124.37mm。將以上注入?yún)?shù)、地層熱力學參數(shù)及井筒結構參數(shù)代入壓力溫度程序進行計算，其中深度迭代步長取50m（不能太大否則誤差較大），求得井深3100m處的壓力50.06MPa(實測壓力52.02MPa)，CO2的溫度101.25(實測溫度101)。壓力、溫度誤差分別為3.76%，2.47%。井筒壓力溫度相圖見圖4。井筒壓力溫度數(shù)據(jù)預測見表2。草8井試注CO2壓力溫度計算誤差均在工

9、程誤差允許范圍之內(nèi)，說明本文收集的地層熱力學數(shù)據(jù)、CO2流體物性數(shù)據(jù)及壓力溫度計算方法是可行的。采用此方法可以討論不同注入溫度、注入速度及注入時間對井底溫度的影響，并可預測其它井在注入狀態(tài)和關井狀態(tài)下井筒壓力溫度數(shù)據(jù)分布。2.2、影響井底溫度因素的討論2.2.1注入溫度、注入時間對井底溫度的影響圖4列出了草8井不同注入時間下的井底溫度的變化，從圖形趨勢分析，隨著注入時間的增加，井底溫度有下降的趨勢，起初溫度變化幅度較大，隨著時間的增加，溫度變化幅度較小，但溫度總體下降幅度較小。從注入時間分析，當注入時間大于50d后，溫度下降幅度不大。從注入溫度的對比分析，注入溫度的大小對井底溫度的影響甚微。說

10、明注入溫度的大小和注入時間對井底溫度的影響均很小。圖4井底溫度隨時間的變化關系曲線圖5草8井試注CO2井筒相圖分布注入速度對井底溫度的影響仍以草8井為例，從圖6計算結果可以看出，注入速度的變化對井底溫度的影響較大。注入速度20t/d與70t/d比較，兩者井底溫度相差5.86。從上述井口注入溫度、注入時間及注入速度對井底溫度影響的分析比較可知，注入速度影響較大，其次是注入時間的影響，注入溫度影響最小。圖6不同注入速度對井底溫度的影響(數(shù)據(jù)取自草8井)表2草8井試注CO2井筒壓力溫度預測數(shù)據(jù)表序號井深(m)溫度()壓力(MPa)密度(Kg/m3)相態(tài)10-20.0024.501006.86液相21

11、000.8725.45956.95液相365030.8530.26871.24液相470032.3930.69868.73氣相5200071.5341.38815.82氣相6300098.2549.28794.13氣相7305099.7549.67792.87氣相83100101.2550.06791.62氣相2.2.3、對草8井正常注入CO2時井筒壓力及溫度分布預測井口注入溫度為20，井口注入壓力30MPa，平均每天注CO221.17t/d，連續(xù)注入時間26.58d，則井筒壓力溫度分布數(shù)據(jù)見表3。從表3數(shù)據(jù)可以看出，注入26.58d之后，井底溫度為103.14，比正常地溫108略有下降。表3

12、草8井注CO2井筒壓力溫度分布（2006.08數(shù)據(jù)）序號井深(m)溫度()壓力(MPa)密度(Kg/m3)相態(tài)1020.0030.00903.71液相210018.2230.89912.64液相365030.8435.80901.60液相470032.3436.24899.45氣相5100041.3438.86886.82氣相6195069.8746.93850.83氣相7200071.3747.35849.23氣相83000100.1455.54822.60氣相93100103.1456.35820.29氣相3、結論(1)、介紹了井底壓力溫度分布的計算方法，并以實際數(shù)據(jù)進行了驗證，誤差均小于

13、5%。(2)、列出了地層熱力學參數(shù)，CO2部分物性參數(shù)，氣體偏差系數(shù)的計算，經(jīng)實例驗證，這些數(shù)據(jù)是可行的。(3)、討論了不同注入時間、注入溫度及注入速度對井底溫度的影響。其中注入速度的影響最大，注入溫度的影響最小。(4)、預測了草8井正常注入時的井筒壓力溫度分布。符號注釋a為地溫悌度，取0.03/m；A為注入管橫截面積，m2；b為地面平均溫度，；C為CO2的熱容，KJ/(Kg.K)；D為注入管內(nèi)徑，m；hc為傳熱系數(shù)，W/(m.K)；hr為輻射熱系數(shù)，W/(m2.K)；f(t)為隨時間變化的導熱傳熱函數(shù)；H為地層深度,m；Gt為注入速度，t/d；K為地層導熱系數(shù)，W/(m.K)；P為壓力，MP

14、a；Pc為氣體臨界壓力,MPa；Q為流量，m3/d；To為注入流體的溫度，；T為氣體絕對溫度,K；Tc為氣體臨界溫度，K；t為注入時間,d；R為氣體常數(shù)，0.008314kJ/kmol.Krto為油管外徑，m；rco為套管外徑，m；rh為井眼半徑，m；V為流體在管中的流速，m/s；為氣體粘度,Pa.s；Ut為總傳熱系數(shù)，W/(m2.K)；為流體的密度，Kg/m3；為熱擴散系數(shù)，m2/h；為摩阻系數(shù)，無因次。Z為氣體的壓縮系數(shù)，無因次。主要參考文獻1、王鴻勛，張琪.采油工藝原理(修訂本）M北京：石油工業(yè)出版社，1994.P282、Ramy,H.J.Jr為“WellboreHeatTransmis

15、sion”，Trans，AIME，225（1962）。JPT(April1962)P427-4353、萬仁溥、羅英俊.采油技術手冊(第8分冊)（修訂本）M北京：石油工業(yè)出版社，2001。P131；P136-1504、石油化學工業(yè)部化工設計院.氮肥工藝設計手冊M北京：石油化學工業(yè)出版社1977.P50，P865、袁恩熙.工程流體力學M.北京石油工業(yè)出版社，1986.P116STUDYOFDISTRIBUTIONOFWELLBOREPRESSUREANDTEMPERATUREOFCO2INJECTIONWELLWTPT2007,28(1)ZhangYongTangRenxuan(OilProduc

16、tionPlantofHuadongBranch,Sinopec)InthelightoftheenergybalanceequationofverticalconduitflowandcombinedwiththecomputationalmethodfortheRamywellboretemperaturedistribution,thepapergivesthedistributionsofwellborepressureandtemperaturewhenCO2isinjectedatdifferentinjectionrates,injectiontemperaturesandinj

17、ectiontimeetc.andverifiesitbyactualdataofinjectiontestatCao-8WellinSubeiOilfield.Theresultisthatthepercentageerrorofpressureis3.8%andthepercentageerroroftemperatureis2.5%.Onthebasisofthis,thepaperfurtherdiscussesthefactorsinfluencingwellboretemperatureofCO2injectionwellandbelievesthattheinjectionrat

18、eofCO2isamaininfluencefactor.KeywordsCO2,wellborepressure,wellboretemperature,forecastFig-1:TheCO2compressmodulusesvsthedifferentpressuresandtempreturesFig-2:TheCO2densitiesvsthethedifferentpressuresandtempreturesFig-3:TheCO2viscostiesvsthethedifferentpressuresandtempreturesFig-4:ThebottomtemperaturevstheinjectingtimesFig-5:TheCO2phasestatechartofCAO-8wellunderinjec