(中國石油化工股份有限公司華北油氣分公司,河南 鄭州 450006)
紅河油田注CO2井筒及油層溫度場分布研究*
劉學(xué)全
(中國石油化工股份有限公司華北油氣分公司,河南 鄭州 450006)
在注CO2提高原油采收率的過程中,注氣井井筒及油層溫度場分布是影響注CO2效果的主要因素之一。結(jié)合注CO2井筒及油層溫度場分布模型對紅河油田注CO2井筒及油層溫度分布進(jìn)行了研究,對影響井筒及油層溫度場分布的因素進(jìn)行了分析。結(jié)果表明:注CO2過程中,井口注入溫度和注入壓力對井底溫度和油層溫度場分布幾乎沒有影響;井口日注入量越大,注入時間越長,井筒溫度和油層溫度越低,CO2波及到的油層范圍越大。為保證紅河油田注CO2的效果,現(xiàn)場試驗中總注入量應(yīng)該有一個合適范圍,井口日注入量應(yīng)控制在30 t以內(nèi),采用3Cr防腐蝕油管可有效防止注CO2對管材腐蝕。
CO2井筒油層溫度場腐蝕
紅河油田長8油藏位于甘肅省鎮(zhèn)原-涇川縣,東北方向與西峰油田相鄰,是中國石油化工股份有限公司華北分公司在鄂爾多斯盆地的主力區(qū)塊。長8油藏儲層溫度65~72 ℃,地層水型為CaCl2水型,礦化度75 406 mg/L,Ca2+質(zhì)量濃度5 486 mg/L;平均孔隙度為10.8%,平均滲透率為0.4×10-3μm2,是典型的低孔超低滲儲層。長8油藏目前主要采用天然能量開發(fā),產(chǎn)量能力遞減較快;儲層滲透率低,注水滲流阻力大,難以驅(qū)動基質(zhì)中的原油并且注水難以見效。
國內(nèi)外研究及現(xiàn)場試驗均表明,注CO2是提高低滲透儲層開發(fā)效果的一種很具潛力的增油措施。CO2可通過降低油水界面張力、降低原油黏度、原油體積膨脹、提高地層滲透率以及溶氣驅(qū)等機理,大幅度提高原油采收率。但是注入CO2后會對儲層溫度造成一定影響。注CO2后井底溫度降低,影響注入井周圍油層溫度場分布,進(jìn)而影響原油物性,決定儲層是否發(fā)生冷傷害。因此注氣過程中注入井井底及井筒周圍油層溫度場分布是影響注CO2效果的主要因素之一[1-4]。
注CO2氣井壓力較高,而且存在腐蝕,容易損壞儀器,對井底溫度進(jìn)行測試費工費時,同時監(jiān)測井筒周圍油層溫度場分布也存在困難。為了有效預(yù)測注CO2過程中注氣井井底及周圍油層溫度場變化規(guī)律,利用井筒溫度梯度模型及Lauwerier油層傳熱模型,結(jié)合CO2熱力學(xué)數(shù)據(jù),對紅河油田注CO2井筒及油層溫度場分布進(jìn)行了研究。
1.1井筒溫度場分布模型
自20世紀(jì)50年代以來,國內(nèi)外學(xué)者對注CO2井井筒溫度場分布進(jìn)行了研究。在已開展注CO2井筒溫度預(yù)測研究中,在狀態(tài)方程選擇時,不同條件下CO2物性參數(shù)(密度、黏度、比熱容和導(dǎo)熱系數(shù)等)的計算模型誤差較大;部分學(xué)者將傳熱對井筒溫度影響簡單處理,沒有將壓力與流體物性進(jìn)行耦合計算,誤差較大[5-6]。選擇陳林[7]等人的井筒溫度、壓力場分布耦合模型,同時對不同溫度、壓力條件下的CO2物性參數(shù)(密度、黏度、比熱容和導(dǎo)熱系數(shù)等)進(jìn)行數(shù)值擬合[8-10],建立不同溫度、壓力下的CO2物性參數(shù)方程,計算注CO2井井筒溫度場分布。
(1)
(2)
(3)
式中:Tf為井筒溫度,℃;Tsur為地表溫度,℃;GDC為地溫梯度,℃/m;Z為井筒深度,m;rto為油管外半徑,m;λe為地層導(dǎo)熱系數(shù),W/(m·K);W為注入流量,kg/s;Cp為注入流體定壓比熱容,kJ/(kg·℃);Uto為以油管外表面為基準(zhǔn)面積的總傳熱系數(shù),W/(m2·K),可依據(jù)文獻(xiàn)[11-12],將井身分成多段,在每一段上根據(jù)不同的溫度、壓力采用迭代法求解;f(τD)為傳熱函數(shù),無因次,可依據(jù)Hasan[13]公式計算;P為井筒壓力,Pa;ρ為油管內(nèi)流體密度,kg/m3;g為重力加速度,9.807 m/s2;θ為注氣井筒傾角,(°);v為油管內(nèi)流體流速,m/s;f為摩阻系數(shù),無因次,可依據(jù)文獻(xiàn)[14]計算。
1.2油層溫度場分布模型
沿井筒注入的CO2進(jìn)入油層多孔介質(zhì)后,熱量傳遞方式既有直接的熱傳導(dǎo),又有隨著注入流體滲流而伴隨的熱量傳遞。注入的CO2在油層中溫度分布可采用Lauwerier[15]方法計算,距注入井任意半徑處的溫度為
(4)
(5)
(6)
(7)
(8)
(9)
式中:xD,tD和λ為無因次參數(shù);λm為頂?shù)讓訉?dǎo)熱系數(shù),W/(m·℃);t為注入時間,s;R為徑向流模型中微元體內(nèi)側(cè)面距井眼的徑向距離,m;h為油層厚度,m;ρinj,ρf和ρm分別為井底流體、油層和頂?shù)讓用芏?,kg/m3;Cinj,Cf和Cm分別為井底流體、油層和頂?shù)讓颖葻崛?,J/(kg·℃);Qinj為井底流體流量,m3/s;Tri為油層初始溫度,℃;Tinj為井底溫度,℃;erfc為余誤差函數(shù)。
2013年6月紅河油田紅河156井組開始注CO2先導(dǎo)試驗,設(shè)計注入量5 000 t。基于紅河156井井筒及地層基本參數(shù),利用井筒及油層溫度場分布模型,改變注入溫度、注入壓力和日注入量等注入工藝參數(shù),對紅河156井注CO2井筒及油層溫度場分布進(jìn)行研究。紅河156井基本計算參數(shù)見表1。
表1 紅河156井基本計算參數(shù)
2.1注入溫度
根據(jù)紅河156井現(xiàn)場注入數(shù)據(jù),計算注入壓力15 MPa,日注入量15 t,總注入量5 000 t時不同注入溫度下的井筒溫度及油層溫度場分布(見圖1和圖2),研究其注入溫度對井筒溫度及油層溫度場影響規(guī)律。
圖1 井筒溫度場分布
由圖1可知,注入溫度在井口附近顯著影響井筒溫度剖面,但隨著井深增大,影響程度減小。當(dāng)井深達(dá)到1 200 m時,注入溫度對井筒溫度幾乎沒有影響;由圖2可知,注入溫度改變對油層溫度場分布幾乎沒有影響,這是因為,在不同的注入溫度下,井底溫度幾乎相等,井底CO2的物性參數(shù)也相差不大,在相同日注入量下CO2與油層之間的熱量傳遞也相差不大,因此油層溫度場分布與注入溫度關(guān)系不大。以井筒為中心,油層溫度隨距井筒徑向距離的增加而升高,當(dāng)增加至一定距離時,溫度恢復(fù)至原始地層溫度。
圖2 油層溫度場分布
2.2注入壓力
根據(jù)紅河156井現(xiàn)場注入數(shù)據(jù),計算注入溫度10 ℃,日注入量15 t,總注入量5 000 t時不同注入壓力下的井筒溫度及油層溫度場分布(見圖3和圖4),研究其對井筒溫度和油層溫度場影響規(guī)律。
圖3 井筒溫度場分布
圖4 油層溫度場分布
從圖3和圖4可以看出,井筒溫度隨注入壓力的增加略有增加,油層溫度也隨注入壓力的增加略有增加,注入的CO2所能波及到的徑向距離油層溫度隨注入壓力的增加略有減小,但總體變化不大。這是由于CO2的比熱容隨著壓力增大而略有減小,CO2在井筒內(nèi)流動過程中地層傳遞給CO2同樣的熱量,CO2溫度升高幅度略有增大;同樣,CO2在油層內(nèi)流動過程中油層傳遞給CO2同樣熱量,CO2溫度升高幅度略有增大。
2.3日注入量
根據(jù)紅河156井現(xiàn)場注入數(shù)據(jù),計算注入溫度10 ℃,注入壓力15 MPa,總注入量5 000 t時不同日注入量條件下井筒溫度及油層溫度場分布(見圖5和圖6),研究其對井筒溫度和油層溫度場的影響規(guī)律。
圖5 井筒溫度場分布
圖6 油層溫度場分布
由圖5可知,隨井口日注入量增大,井筒溫度降低,并且降低幅度明顯。日注入量增大,相同時間內(nèi)注入CO2量增大,從地層吸收同樣熱量引起溫度升高的幅度隨之降低;從圖6可知,日注入量對近井筒地帶的油層溫度場分布影響是比較明顯的。日注入量越大,井底溫度越低,在相同徑向距離下油層的溫度越低。日注入量增大對井筒壓力分布的影響并不明顯。
根據(jù)楊勝來等[16]的研究,在注CO2過程中,油層中可能會出現(xiàn)4種油層溫度與原油黏度帶模型。根據(jù)上述計算結(jié)果,結(jié)合長8油藏原油黏溫曲線及CO2溶于原油后的物性變化,可以得出每日注入不同量的CO2條件下油層溫度與原油黏度帶,見圖7。
圖7 油層溫度與原油黏度帶
由圖7可知,在日注入量10 t時,CO2降低的油層最低溫度略低于但接近原始油層溫度,雖然降溫使原油增黏,但CO2降黏效果比降溫增黏作用更明顯,因此原油黏度降低,CO2增加原油采收率效果最明顯,是最理想的情況;當(dāng)日注入量達(dá)到30 t時,CO2降低的油層最低溫度低于油層原始溫度,但高于長8油藏原油凝固點(21 ℃)。低溫導(dǎo)致原油黏度增加的儲層傷害顯著于其降黏作用,此時,對儲層形成了一定程度的傷害,降低了驅(qū)油效率,井筒周圍出現(xiàn)模型一定半徑范圍的低溫增黏帶,日注入量越大,低溫增黏帶的半徑越大,對增產(chǎn)越不利;當(dāng)日注入量增加到500 t時,CO2降低的油層最低溫度低于長8油藏原油凝固點,油層中出現(xiàn)一個原油結(jié)蠟帶,原油結(jié)蠟會堵塞地層,造成油井產(chǎn)液量降低,甚至停產(chǎn)。因此為提高注CO2的開發(fā)效果,現(xiàn)場日注入量應(yīng)控制在30 t之內(nèi)。
2.4注入時間
根據(jù)紅河156井現(xiàn)場的注入數(shù)據(jù),計算注入溫度10 ℃,注入壓力15 MPa,日注入量15 t,總注入量5 000 t時不同注入時間下井筒溫度和油層溫度場分布的影響規(guī)律(見圖8和圖9)。
從圖8和圖9可以看出,起初隨井口注入時間增加,井底溫度下降幅度較大;隨著注入時間的繼續(xù)增加,井底溫度降幅逐漸變小。當(dāng)注入時間大于2 000 h時,注入時間對井底溫度影響很??;初始階段近井油層溫度下降較快,隨著時間增加(即注入量增大),低溫帶波及到的體積逐漸增大,且油層溫度逐漸降低,但推進(jìn)速度逐漸變慢。這主要是因為注入初期,注入流體波及到的油層半徑小,在相同的注入流量及其攜帶的熱量作用下,低溫帶的推進(jìn)速度較快;隨著時間的增加,注入流體波及到的油層半徑增加,使得降低相同半徑油層的溫度所需要的低溫流體也越多,因而推進(jìn)速度變慢。此外,隨著半徑增加,與相鄰油層、井底交換的熱量也就越多,對注入流體前沿的影響也就越小。因此現(xiàn)場試驗中的總注入量應(yīng)該有一合適的范圍,并不是越大越好。
圖8 井筒溫度場分布
圖9 油層溫度場分布
根據(jù)前期調(diào)研試驗掛片評價結(jié)果,3Cr材質(zhì)油管抗腐蝕性較強,因此注入井采用3Cr防腐蝕油管和耐腐蝕高壓井口以及環(huán)空隔離保護(hù)液技術(shù),能有效防止注CO2過程中對管材的腐蝕。
(1)注CO2過程中,井口注入溫度在井口附近顯著影響井筒溫度剖面,但對井底溫度及油層溫度場分布影響不大,注入壓力對井筒和油層溫度場分布幾乎沒有影響;
(2)注CO2過程中,井口日注入量越大,井底溫度和油層溫度越低,為保證紅河油田注CO2開發(fā)效果,井口日注入量應(yīng)控制在30 t以內(nèi);
(3)注CO2過程中,注入時間越長,注入量越大,井底溫度越低,CO2波及到的油層范圍越大,且油層溫度越低,現(xiàn)場試驗中的總注入量應(yīng)該有一合適的范圍,并不是越大越好;
(4)采用3Cr防腐蝕油管可有效防止注CO2對管材腐蝕。
[1] 吳曉東,王慶,何巖峰.考慮相態(tài)變化的注CO2井井筒溫度壓力場耦合計算模型[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2009,33(1):73-77.
[2] 陳林,孫雷,彭彩珍,等.注CO2井筒溫度場分布規(guī)律模擬研究[J].斷塊油氣田,2009,16(6):82-84.
[3] 王慶,吳曉東,劉長宇,等.CO2井筒相變流動溫度壓力計算模型研究[J].鉆采工藝,2009,32(4):39-41.
[4] 劉洪波,程林松,宋立新,等.注超臨界氣體井筒溫度壓力場計算方法[J].石油大學(xué)學(xué)報(自然科學(xué)版),2004,28(1):52-54.
[5] 唐人選.CO2注入井井筒壓力溫度分布研究[J].試采技術(shù),2007,28(1):53-57.
[6] 張勇,唐人選.CO2井筒壓力溫度的分布[J].海洋石油,2007,27(2):59-64.
[7] 陳林.注CO2井筒及油層溫度場分布規(guī)律模擬研究[D].成都:西南石油大學(xué),2008.
[8] Span R,Wagner W.A new equation of state for carbon dioxide covering the fluid region from the triple-point temperature to 1 100 K at pressures up to 800 MPa[J].Journal of Physical and Chemical Reference Data,1996,25(6):1509-1596.
[9] Fenghour A,Wakeham W A.The viscosity of carbon dioxide[J].Journal of Physical and Chemical Reference Data,1998,27(1):31-44.
[10] 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.
[11] Ramey,JR H J.Wellbore heat transmission[J].Journal of Petroleum Technology,1962,14(4):427-435.
[12] Dropkin D,Somerscales E.Heat transfer by natural convection in liquids confined by two parallel plates which are inclined at various angles with respect to tie horizontal[J].Journal of Heat Transfer,1965,87(1):77-82.
[13] Hasan A R,Kabir C S.Aspects of wellbore heat transfer during two-phase flow[J].SPE Production & Facilities,1994,9(3):211-216.
[14] Jain A K.Accurate explicit equation for friction factor[J].Journal of the Hydraulics Division,1976,102(5):674-677.
[15] Lauwerier H A.The transport of heat in an oil layer caused by the injection of hot fluid[J].Applied Scientific Research,1955(5):145-150.
[16] 楊勝來,何建軍,榮光迪.CO2吞吐期間井筒及油層溫度場及其對CO2吞吐效果的影響[J].石油鉆探技術(shù),2004,32(2):47-49.
[17] 郭炳仁.CO2對石油及氣體管線的腐蝕[J].石油化工環(huán)境保護(hù),1994(2):41-44.
[18] 姚曉.氣田開發(fā)中CO2對井內(nèi)管材的腐蝕與預(yù)防[J].鉆采工藝,1996,19(6):31-37.
(編輯 寇岱清)
ResearchonTemperatureDistributionofCO2InjectionWellboreandOilReservoirinHongheOilField
LiuXuequan
(SINOPECNorthChinaCompany,Zhengzhou450006,China)
Temperature distribution of CO2injection wellbore and oil reservoir is a key factor in the oil recovery process enhanced by CO2injection. Related influencing factors were analyzed through the simulation of temperature distribution. The results showed that temperature and pressure in the inlet of wellbore had little effect on the bottom-hole temperature and temperature distribution of oil reservoir; temperature of wellbore and oil reservoir decreased with the increasing of injection rate and injection time, which was opposite with the tendency of reservoir range swept by CO2. In order to ensure the effect of CO2injection, the injection volume should be kept in a reasonable range of 30 t/d, and 3Cr anticorrosion pipe should be used to prevent the CO2corrosion effectively.
CO2, wellbore, temperature distribution of oil reservoir, corrosion
2017-03-09;修改稿收到日期:2017-07-26。
劉學(xué)全(1980—),本科,工程師,2005年畢業(yè)于中國石油大學(xué)(北京)石油工程學(xué)院油氣儲運工程專業(yè),現(xiàn)在該公司工程技術(shù)研究院從事機械采油和提高采收率工作。E-mail:lxq811223@163.com
國家科技重大專項“低豐度致密低滲油氣藏開發(fā)關(guān)鍵技術(shù)”(2016ZX05048)。