電加熱稠油熱采井筒溫度場數值計算方法

朱廣海 劉章聰 熊旭東 宋洵成 王軍恒 翁博

引用本文:

電加熱稠油熱采井筒溫度場數值計算方法

    作者簡介: 朱廣海(1985—),男,河北衡水人,2008年畢業于中國石油大學(北京)石油工程專業,工程師,主要從事采油采氣研究工作。Email:[email protected]
    通訊作者: 宋洵成, [email protected]
  • 基金項目:

    國家科技重大專項課題“超重油油藏冷采穩產與改善開發效果技術”(編號:2016ZX05031-001)資助

  • 中圖分類號: TE357.44

Numerical Calculation Method of the Wellbore Temperature Field for Electric Heating Heavy Oil Thermal Recovery

    Corresponding author: SONG Xuncheng, [email protected] ;
  • CLC number: TE357.44

  • 摘要: 電加熱稠油熱采井筒溫度場是熱采作業參數設計的核心依據,基于傳熱學理論和氣液兩相流井筒溫壓場計算方法,考慮溫度對稠油熱物性影響,建立了連續電加熱和電磁短節加熱工藝井筒溫度場的數值計算方法,并以大港油田X井為例,計算了不同加熱功率下連續電加熱和電磁短節加熱工藝的井筒溫度場。計算結果表明:井口溫度的模型計算結果與實測值相對誤差僅為3.10%,滿足工程設計精度要求,也驗證了計算方法的有效性和準確性;連續電加熱工藝的井筒溫度剖面平滑連續,而電磁短節加熱工藝的井筒溫度剖面呈鋸齒形,且溫度波動更大;連續電加熱工藝的井口溫度高于電磁短節加熱,而連續電加熱工藝的平均溫度則低于電磁短節加熱工藝。該研究結果可為電加熱稠油熱采工藝選擇、作業參數設計提供指導和借鑒。
  • 圖 1  油管內產液與井筒/地層界面間的換熱剖面

    Figure 1.  Heat transfer profile between the produced fluids in the tubing and the wellbore/formation interface

    圖 2  連續電加熱與電磁短節加熱井筒網格劃分

    Figure 2.  Wellbore meshing of continuous tube electric heating and electromagnetic nipple heating

    圖 3  連續電加熱工藝的井筒溫度場剖面

    Figure 3.  Wellbore temperature field profile of continuous electric heating

    圖 4  電磁短節加熱工藝的井筒溫度場剖面

    Figure 4.  Wellbore temperature field profile of electromagnetic nipple heating

    圖 5  連續電加熱和電磁短節加熱工藝的井口溫度與平均溫度對比

    Figure 5.  Comparison on the wellhead temperatures and average temperatures of continuous electric heating and electromagnetic nipple heating processes

    表 1  大港油田X井實鉆井身結構

    Table 1.  Actual casing program of Well X in Dagang Oilfield

    套管外徑/mm井眼直徑/mm套管下深/m水泥返高/m
    表層套管244.50311.1290.00地面
    生產套管139.70215.91 388.00 865.00
    下載: 導出CSV

    表 2  連續電加熱和電磁短節加熱工藝的井筒溫度場剖面均方差比較

    Table 2.  Comparison of the mean square errors of wellbore temperature field profiles formed by continuous electric heating and electromagnetic nipple heating processes

    加熱方法不同加熱功率的井筒溫度場剖面均方差
    20 kW40 kW60 kW80 kW100 kW
    連續電加熱11.638.405.312.973.82
    電磁短節加熱14.4010.96 7.404.404.01
    下載: 導出CSV
    南国七星彩高手论坛
  • [1] 劉文章. 普通稠油油藏二次熱采開發模式綜述[J]. 特種油氣藏, 1998, 5(2): 1–7.LIU Wenzhang. Summary of secondary thermal recovery for conventional heavy oil reservoir[J]. Special Oil & Gas Reservoirs, 1998, 5(2): 1–7.
    [2] 梁丹,馮國智,曾祥林,等. 海上稠油兩種熱采方式開發效果評價[J]. 石油鉆探技術, 2014, 42(1): 95–99. doi: 10.3969/j.issn.1001-0890.2014.01.019LIANG Dan, FENG Guozhi, ZENG Xianglin, et al. Evaluation of two thermal methods in offshore heavy oilfields development[J]. Petroleum Drilling Techniques, 2014, 42(1): 95–99. doi: 10.3969/j.issn.1001-0890.2014.01.019
    [3] 唐述凱,李明忠,劉陳偉,等. 閉式熱流體循環井筒溫度分布規律研究[J]. 石油鉆探技術, 2014, 42(5): 104–108.TANG Shukai, LI Mingzhong, LIU Chenwei, et al. Temperature distribution in closed wellbore with hot fluid circulation[J]. Petroleum Drilling Techniques, 2014, 42(5): 104–108.
    [4] 任瑛. 加熱開采稠油工藝的探討:井筒中的熱流體循環[J]. 華東石油學院學報(自然科學版), 1982, 6(4): 53–66.REN Ying. Thermal production of viscous oil with recycling of hot fluid in well bore[J]. Journal of East China Petroleum Institute (Edition of Natural Science), 1982, 6(4): 53–66.
    [5] 闞慶山,江永全,徐凌云,等. 大港油田火成巖油藏稠油熱采工藝[J]. 石油鉆采工藝, 2009, 31(2): 115–117.KAN Qingshan, JIANG Yongquan, XU Lingyun, et al. Technology of thermal recovery of volcanic high-viscosity reservoirs in Dagang Oilfield[J]. Oil Drilling & Production Technology, 2009, 31(2): 115–117.
    [6] 李偉超,劉平,于繼飛,等. 渤海稠油油田井筒電加熱技術可行性分析[J]. 斷塊油氣田, 2012, 19(4): 513–516.LI Weichao, LIU Ping, YU Jifei, et al. Feasibility analysis on electric heating of wellbore in heavy oilfield of Bohai[J]. Fault-Block Oil & Gas Field, 2012, 19(4): 513–516.
    [7] 余五星, 潘建華, 林軍, 等. 超稠油中頻電加熱采油技術[J]. 特種油氣藏, 2000, 7(增刊1): 31–33.YU Wuxing, PAN Jianhua, LIN Jun, et al. Thermal producing of super heavy oil with frequency current[J]. Special Oil & Gas Reservoirs, 2000, 7(Supplement1): 31–33.
    [8] 宋時權,張公社,李晶晶,等. 電加熱桿抽油井溫度分布計算[J]. 斷塊油氣田, 2008, 15(3): 121–123.SONG Shiquan, ZHANG Gongshe, LI Jingjing, et al. Calculation of temperature distribution of electrical heating rod in pumping well[J]. Fault-Block Oil & Gas Field, 2008, 15(3): 121–123.
    [9] 姚傳進,雷光倫,吳川,等. 高凝原油井筒溫度場影響因素研究[J]. 石油鉆探技術, 2011, 39(5): 74–78. doi: 10.3969/j.issn.1001-0890.2011.05.016YAO Chuanjin, LEI Guanglun, WU Chuan, et al. Study of the factors impacting on wellbore temperature in high pour point oil production[J]. Petroleum Drilling Techniques, 2011, 39(5): 74–78. doi: 10.3969/j.issn.1001-0890.2011.05.016
    [10] 盧智慧,何雪芹,劉志恒,等. 考慮環境因素的井筒流動溫度場計算[J]. 斷塊油氣田, 2016, 23(5): 652–654.LU Zhihui, HE Xueqin, LIU Zhiheng, et al. Calculation of wellbore flowing temperature distribution considering environmental factors[J]. Fault-Block Oil & Gas Field, 2016, 23(5): 652–654.
    [11] 李夢博,柳貢慧,李軍,等. 考慮非牛頓流體螺旋流動的鉆井井筒溫度場研究[J]. 石油鉆探技術, 2014, 42(5): 74–79.LI Mengbo, LIU Gonghui, LI Jun, et al. Research on wellbore temperature field with helical flow of non-Newtonian fluids in drilling operation[J]. Petroleum Drilling Techniques, 2014, 42(5): 74–79.
    [12] 王江帥,李軍,柳貢慧,等. 循環鉆進過程中井筒溫度場新模型[J]. 斷塊油氣田, 2018, 25(2): 240–243.WANG Jiangshuai, LI Jun, LIU Gonghui, et al. New model of wellbore temperature field during drilling process[J]. Fault-Block Oil and Gas Field, 2018, 25(2): 240–243.
    [13] RAMEY H J Jr. Wellbore heat transmission[J]. Journal of Petroleum Technology, 2013, 14(4): 427–435.
    [14] 宋洵成,韋龍貴,何連,等. 氣液兩相流循環溫度和壓力預測耦合模型[J]. 石油鉆采工藝, 2012, 34(6): 5–9. doi: 10.3969/j.issn.1000-7393.2012.06.004SONG Xuncheng, WEI Longgui, HE Lian, et al. Acoupled model for wellbore temperature and pressure prediction of gas-liquid drilling fluid[J]. Oil Drilling & Production Technology, 2012, 34(6): 5–9. doi: 10.3969/j.issn.1000-7393.2012.06.004
    [15] HASAN A R, KABIR C S, SAYARPOUR M. Simplified two-phase flow modeling in wellbores[J]. Journal of Petroleum Science and Engineering, 2010, 72(1/2): 42–49.
    [16] BEGGS H D, BRILL J P. A study of two-phase flow in inclined pipes[J]. Journal of Petroleum Technology, 1973, 25(5): 607–617. doi: 10.2118/4007-PA
    [17] ANSARI A M, SYLVESTER N D, SHOHAM O, et al. A comprehensive mechanistic model for upward two-phase flow in wellbores[R]. SPE 20630, 1990.
  • [1] 孫小輝孫寶江王志遠王金堂 . 超臨界CO2鉆井井筒水合物形成區域預測. 石油鉆探技術, 2015, 43(6): 13-19. doi: 10.11911/syztjs.201506003
    [2] 羅艷艷程林松高山王宏遠 . 水平井蒸汽驅溫度場與黏度場平面物理模擬. 石油鉆探技術, 2012, 40(1): 74-77. doi: 10.3969/j.issn.1001-0890.2012.01.015
    [3] 李夢博柳貢慧李軍魏曉強高海軍 . 考慮非牛頓流體螺旋流動的鉆井井筒溫度場研究. 石油鉆探技術, 2014, 42(5): 74-79. doi: 10.11911/syztjs.201405013
    [4] 姚傳進雷光倫吳川高達蔣寶云劉海慶 . 高凝原油井筒溫度場影響因素研究. 石油鉆探技術, 2011, 39(5): 74-78. doi: 10.3969/j.issn.1001-0890.2011.05.016
    [5] 于法浩蔣召平白健華劉義剛孟祥海 . 渤海油田稠油水平井防砂篩管耐溫能力的確定. 石油鉆探技術, 2018, 46(6): 65-70. doi: 10.11911/syztjs.2018010
    [6] 萬立夫李根生田守嶒宋先知王海柱 . 鉆井過程中硫化氫侵入井筒后的賦存狀態研究. 石油鉆探技術, 2013, 41(6): 29-33. doi: 10.3969/j.issn.1001-0890.2013.06.006
    [7] 鄭春峰魏琛張海濤李昂孟紅霞 . 海上油井井筒結蠟剖面預測新模型. 石油鉆探技術, 2017, 45(4): 103-109. doi: 10.11911/syztjs.201704018
    [8] 唐述凱李明忠劉陳偉張國棟 . 閉式熱流體循環井筒溫度分布規律研究. 石油鉆探技術, 2014, 42(5): 104-108. doi: 10.11911/syztjs.201405019
    [9] 吳雪婷鄒韻陸彥穎趙增義周城漢 . 漏失循環條件下井筒溫度預測與漏層位置判別. 石油鉆探技術, 2019, 47(6): 1-6. doi: 10.11911/syztjs.2019119
    [10] 王通孫永濤白健華馬增華蔣召平李海濤 . 海上耐高溫井下安全控制管柱系統的研制. 石油鉆探技術, 2017, 45(6): 49-54. doi: 10.11911/syztjs.201706009
    [11] 楊順輝豆寧輝趙向陽柯珂王志遠 . 多層合采智能井井筒溫度場預測模型及應用. 石油鉆探技術, 2019, 47(4): 83-91. doi: 10.11911/syztjs.2019049
    [12] 馬帥張風波洪楚僑劉雙琪鐘家峻王世朝 . 多層合采階梯井產能計算模型的建立與求解. 石油鉆探技術, 2015, 43(5): 94-99. doi: 10.11911/syztjs.201505016
    [13] 何淼柳貢慧李軍李夢博查春青李根 . 多相流全瞬態溫度壓力場耦合模型求解及分析. 石油鉆探技術, 2015, 43(2): 25-32. doi: 10.11911/syztjs.201502005
    [14] 李根生竇亮彬田守嶒宋先知 . 酸性氣體侵入井筒瞬態流動規律研究. 石油鉆探技術, 2013, 41(4): 8-14. doi: 10.3969/j.issn.1001-0890.2013.04.003
    [15] 何漢平 . 油氣井井筒完整性系統風險評估方法. 石油鉆探技術, 2017, 45(3): 72-76. doi: 10.11911/syztjs.201703013
    [16] 董長銀張清華崔明月王鵬高彥才李效波 . 復雜條件下疏松砂巖油藏動態出砂預測研究. 石油鉆探技術, 2015, 43(6): 81-86. doi: 10.11911/syztjs.201506015
    [17] 陳威 . 井下外磁轉子式渦輪發電機設計與試驗研究. 石油鉆探技術, 2015, 43(6): 114-119. doi: 10.11911/syztjs.201506021
    [18] 李猛賀會群辛永安婁爾標張云飛 . 連續管鉆井電液定向裝置工具面調整方法. 石油鉆探技術, 2016, 44(6): 48-54. doi: 10.11911/syztjs.201606008
    [19] 王正旭高德利 . 高頻電磁加熱稠油儲層溫度分布及其影響因素分析. 石油鉆探技術, 2019, 47(6): 1-8. doi: 10.11911/syztjs.2019128
    [20] 何保生劉華亮劉剛楊全枝耿占力 . 基于球面波傳播理論的鉆頭趨近鄰井防碰預測模型. 石油鉆探技術, 2013, 41(3): 62-66. doi: 10.3969/j.issn.1001-0890.2013.03.012
  • 加載中
圖(5)表(2)
計量
  • 文章訪問數:  219
  • HTML全文瀏覽量:  126
  • PDF下載量:  16
  • 被引次數: 0
出版歷程
  • 收稿日期:  2019-06-25
  • 錄用日期:  2019-08-25
  • 網絡出版日期:  2019-09-02
  • 刊出日期:  2019-09-01

電加熱稠油熱采井筒溫度場數值計算方法

    通訊作者: 宋洵成, [email protected]
    作者簡介: 朱廣海(1985—),男,河北衡水人,2008年畢業于中國石油大學(北京)石油工程專業,工程師,主要從事采油采氣研究工作。Email:[email protected]
  • 1. 中國石油大港油田分公司石油工程研究院,天津 300280
  • 2. 中國石油勘探開發研究院,北京 100083
  • 3. 中國石油新疆油田分公司工程技術處,新疆克拉瑪依 834000
  • 4. 中國石油大學(華東)石油工程學院,山東青島 266580
基金項目:  國家科技重大專項課題“超重油油藏冷采穩產與改善開發效果技術”(編號:2016ZX05031-001)資助

摘要: 電加熱稠油熱采井筒溫度場是熱采作業參數設計的核心依據,基于傳熱學理論和氣液兩相流井筒溫壓場計算方法,考慮溫度對稠油熱物性影響,建立了連續電加熱和電磁短節加熱工藝井筒溫度場的數值計算方法,并以大港油田X井為例,計算了不同加熱功率下連續電加熱和電磁短節加熱工藝的井筒溫度場。計算結果表明:井口溫度的模型計算結果與實測值相對誤差僅為3.10%,滿足工程設計精度要求,也驗證了計算方法的有效性和準確性;連續電加熱工藝的井筒溫度剖面平滑連續,而電磁短節加熱工藝的井筒溫度剖面呈鋸齒形,且溫度波動更大;連續電加熱工藝的井口溫度高于電磁短節加熱,而連續電加熱工藝的平均溫度則低于電磁短節加熱工藝。該研究結果可為電加熱稠油熱采工藝選擇、作業參數設計提供指導和借鑒。

English Abstract

  • 稠油的黏度大、流動性差,且其黏度對溫度特別敏感,溫度每升高8~9 ℃,黏度可降低50%,因而提高溫度是改善稠油流動性的有效措施[1],因此,普遍采用熱采工藝開發稠油。目前,稠油熱采工藝主要有電加熱、熱流體循環、蒸汽吞吐和蒸汽驅等[26],其中,電加熱工藝是稠油開發的主要選擇[7]。現場常用電加熱工藝主要有空心油桿電加熱、伴熱帶電加熱和電磁短節加熱等,按照加熱介質和加熱功率是否連續,可將電加熱工藝分為連續電加熱和電磁短節加熱兩類。由于稠油黏溫特性、油井井身結構的不同,電加熱工藝選擇、加熱功率和加熱時長等作業參數設計,均需要精確計算井筒溫度場。因此,深入研究電加熱工藝的井筒與儲層間的換熱機理,建立換熱模型和溫度場計算方法,進而獲取不同加熱工藝和作業參數下井筒溫度場的分布特征,對電加熱稠油熱采工藝選擇、作業參數設計和提高稠油開采效果具有重要意義。

    國內外對電加熱稠油熱采換熱問題的研究主要集中于連續電加熱工藝,而對電磁短節加熱工藝井筒溫度場的研究較少。此外,溫度不僅對稠油黏度影響較大,還對其比熱容和熱導率2個熱物性參數有較大影響,而現有模型未考慮溫度對稠油熱物性的影響[89]。為此,筆者考慮溫度場工程計算精度需求和數值計算方法的可靠性,耦合半瞬態換熱分析方法[1014]和基于流型的氣液兩相流機理模型[1517],建立了考慮溫度對稠油熱物性影響的電加熱稠油熱采流動與換熱控制方程,形成了連續電加熱和電磁短節加熱井筒溫度場的數值計算方法,并用計算實例分析了2種電加熱工藝的井筒溫度場剖面特征、加熱功率對2種工藝井口溫度及平均溫度的影響。

    • 基于電加熱稠油熱采工藝,做以下基本假設:1)油管內流體為一維穩態流動和傳熱,流速、壓力、溫度只隨軸向位置變化而變化;2)地層內僅發生徑向換熱,相同深度地層為均質地層;3)產出液可壓縮,熱物性隨溫度變化而變化;4)忽略生產期間的軸向熱傳導換熱。

    • 由于產出液可壓縮,其物性參數受溫度場和壓力場共同影響,故需耦合求解質量守恒方程、動量守恒方程和能量守恒方程。以井筒中心線為Z軸,井口指向井底的方向為正,建立一維坐標系,Z軸原點為井深參考點,井斜角為θ

    • 油管內流體流動方向與Z軸正方向相反,兩相流質量流量不隨軸向位置變化,質量守恒方程為:

      $\frac{{{{\partial}} \left( {{\rho _{\rm{m}}}{v_{\rm{m}}}A} \right)}}{{{{\partial}} Z}} = 0$

      式中:$\,{\rho _{\rm{m}}}$油管內為產出液平均密度,kg/m3${v_{\rm{m}}}$為產出液平均流速,m/s;$A$為油管橫截面積,m2

    • 油管內流體動量守恒方程為:

      $ - \frac{{{{\partial}} ({\rho _{\rm{m}}}v_{\rm{m}}^2)}}{{{{\partial}} Z}} = \frac{{{{\partial}} p}}{{{{\partial}} Z}} - {\rho _{\rm{m}}}g\cos\;\theta - \frac{{2{\rho _{\rm{m}}}{f_{\rm{m}}}v_{\rm{m}}^2}}{{{d_{{\rm{ti}}}}}}$

      式中:p為產出液壓力,Pa;θ為井斜角,(°);fm為范寧摩阻系數;dti為油管內徑,m。

    • 微元控制體發生的能量傳遞過程有:Z方向對流換熱,即單位時間內流入和流出控制體的能量(包括動能、勢能和焓);產出液與地層的換熱量以及電加熱的生熱量。依據能量守恒原理,得電加熱能量守恒方程為:

      $G\frac{\partial }{{\partial {{Z}}}}\left( { - gZ\cos\;\theta + \frac{{v_{\rm{m}}^2}}{2} + {H_{\rm{m}}}} \right) + \frac{{{T_{\rm{f}}} - {T_{\rm{m}}}}}{{{R_{\rm{f}}}}} + q = 0$

      式中:G為油管內產出液質量流量,kg/s;Hm為油管內產出液比焓,J/kg;Tf為地層溫度,K;Tm為油管內產出液溫度,K;Rf為產出液到地層總熱阻,K/(W·m);q為油管內單位長度的生熱量,W/m。

    • 1)稠油黏溫關系方程。由于在一定溫度下稠油密度變化較小,近似地認為稠油的動力黏度與溫度的關系在ASTM坐標圖上也呈直線關系,其精度能滿足熱采工程計算的要求:

      $ \lg \lg{\mu _{\rm{oD}}} = A - B\lg \left( {{T_{\rm{m}}} - 273.15} \right) $

      式中:μoD為稠油脫氣黏度,Pa·s;AB為常數。

      2)稠油比熱容方程。稠油比熱容主要受溫度和密度的影響,采用Gambill關系式計算稠油比熱容:

      $ {c_{\rm{o}}} = \frac{1}{{\sqrt {{\rho _{15}/1000}} }}\left[ {1.684\;8 + 0.003\;39\left( {T - 273.15} \right)} \right] $

      式中:co為稠油比熱容,kJ/(kg·K);ρ15為15 ℃時稠油密度,kg/m3T為稠油溫度,K。

      3)稠油熱導率方程。稠油的熱導率隨溫度升高而減小,且受稠油密度影響,計算公式為:

      $ {k_{\rm{o}}} = \left( {134.257\;5 - 0.063\;18T} \right)/{\rho _{15}} $

      式中:ko為稠油熱導率,W/(m·K)。

    • 油管內產出液與地層間的換熱剖面如圖1所示,其換熱主要為油管壁面處強迫對流換熱,油管壁、水泥環、套管等的傳導換熱,油套環空內自然對流換熱和輻射換熱。穩態換熱工程意義為軸向位置相同介質(產液、油管壁面、油套環空、套管、水泥環及地層)中發生的徑向換熱熱流量相等,即產出液與井壁的熱流量等于井壁與地層的換熱量,據此可得到產出液與地層的總換熱熱阻。

      圖  1  油管內產液與井筒/地層界面間的換熱剖面

      Figure 1.  Heat transfer profile between the produced fluids in the tubing and the wellbore/formation interface

      產出液到井壁界面(rw)的熱阻Rw為:

      $\begin{array}{l} \quad\! {R_{\rm{w}}} = \dfrac{1}{{\Delta {{Z}}}}\left[ {\dfrac{1}{{2{\text{π}} {r_{{\rm{ti}}}}{h_{\rm{m}}}}} + \dfrac{{\ln ({r_{{\rm{to}}}}/{r_{{\rm{ti}}}})}}{{2{\text{π}} {k_{\rm{t}}}}} + } \right.\\ \quad\quad\quad \left. {\dfrac{{\ln ({r_{{\rm{ci}}}}/{r_{{\rm{to}}}})}}{{2{\text{π}} {k_{\rm{a}}}}} + \dfrac{{\ln ({r_{{\rm{co}}}}/{r_{{\rm{ci}}}})}}{{2{\text{π}} {k_{\rm{c}}}}} + \dfrac{{\ln ({r_{\rm{w}}}/{r_{{\rm{co}}}})}}{{2{\text{π}} {k_{\rm{s}}}}}} \right] \end{array}$

      式中:Rw為產出液到井壁的熱阻,K/(W·m);ktkakcks分別為油管、油套環空、套管、套管與地層環空的熱導率,W/(m·K);hm為產出液與管壁間的對流換熱系數,W/(m2·K);rti是油管內半徑,m;rto是油管外半徑,m;rci是套管內半徑,m; rco是套管外半徑,m。

      由于井筒與地層間的換熱量等于井筒內產出液與井壁間的換熱量,則得:

      $\frac{{{T_{\rm{w}}} - {T_{\rm{m}}}}}{{{R_{\rm{w}}}}} = \frac{{2{\text{π}} {k_{\rm{f}}} ({T_{\rm{f}}} - {T_{\rm{w}}})\Delta {{Z}}}}{{f(t)}}{\rm{ }}$

      式中:Tw為井壁溫度,K;f(t)為Ramey時間函數;kf為地層熱導率,W/(m·K)。

      消去井壁溫度,得到用地層溫度和產出液溫度表示的熱流量計算公式:

      $q = \frac{{{T_{\rm{f}}} - {T_{\rm{m}}}}}{{{R_{\rm{w}}} + \dfrac{{f(t)}}{{2{\text{π}} {k_{\rm{f}}}\Delta {{Z}}}}}}$

      因此,油管內產出液與地層間的總換熱熱阻計算公式為:

      ${R_{\rm{f}}} = {R_{\rm{w}}} + \frac{{f(t)}}{{2{\text{π}} {k_{\rm{f}}}\Delta {{Z}}}}$

    • 稠油自井底向井口流動過程中,溫度、壓力隨井深變化,且溫度和壓力與稠油的密度、黏度、比熱容、熱導率相互影響,因此,求解電加熱稠油熱采控制方程時,需要將全井段按照井深進行網格劃分,應用離散格式控制方程,耦合求解每個網格的溫度、壓力、熱物性、截面含氣率等參數。

    • 連續電加熱和電磁短節加熱2種稠油熱采工藝的熱源布置方式不同,連續電加熱工藝是由空心油桿內加熱電纜或油桿旁邊的加熱電纜沿井筒軸線方向連續提供熱量,故泵掛深度以下井段和油桿段均可按照自定義軸向間距ΔZ均勻劃分網格,如圖2(a)所示。電磁短節加熱工藝是由電纜供電,主要熱量由電磁加熱短節分散提供,同時,電纜也會發熱而成為連續低功率熱源,網格劃分時除按照軸向步長均勻劃分網格外,還需在電磁加熱短節處增加相應長度的網格,如圖2(b)所示。

      圖  2  連續電加熱與電磁短節加熱井筒網格劃分

      Figure 2.  Wellbore meshing of continuous tube electric heating and electromagnetic nipple heating

      此外,為避免壓力計算溢出,采用交錯網格方式將溫度節點布置在網格控制體中心,壓力和速度節點布置在網格控制體上下2個邊界處。

    • 由于能量守恒方程中既有溫度項又有焓項,不利于溫度場求解,因此應用焓的溫壓依賴關系對能量守恒方程進行變換,得到用溫度表示的能量守恒方程。

      真實氣體的焓熱力學微分關系式為:

      $\frac{{{\rm{d}}H}}{{{\rm{d}}Z}} = {c_{\rm{p}}}\frac{{{\rm{d}}T}}{{{\rm{d}}Z}} + {H_{\rm{p}}}\frac{{{\rm{d}}p}}{{{\rm{d}}Z}}$

      式中:cp為定壓比熱容,cp =(dH/dT)pHp為焓變,Hp= (dH/dp)TcpHp可應用真實氣體狀態方程計算。

      將式(11)代入式(3),化簡可得用溫度表示的能量守恒方程:

      $\begin{array}{l} - G{c_{\rm{p}}}\dfrac{{{\rm{d}}{T_{\rm{m}}}}}{{{\rm{d}}Z}} + \dfrac{{{T_{\rm{f}}} - {T_{\rm{m}}}}}{{{R_{\rm{f}}}{\rm{d}}Z}} - Gg\cos\;\theta +\\ \quad G\dfrac{{\rm{d}}}{{{\rm{d}}Z}}\left( {\dfrac{{v_{\rm{m}}^2}}{2}} \right) + G{H_{\rm{p}}}\dfrac{{{\rm{d}}p}}{{{\rm{d}}Z}} + q = 0 \end{array}$

      對于節點i而言,已知控制體入口溫度Ti+1和對應垂深處地層原始溫度Tf,i,則節點i的溫度Ti計算公式為:

      $\begin{array}{l} \quad\qquad{T_i} = \dfrac{{G{c_{{\rm{pm}}}}{T_{i + 1}} + {T_{{\rm{f}},i}}/{R_{{\rm{f}},i}} - Gg\cos \theta \Delta {z_i}}}{{G{c_{{\rm{pm}}}} + 1/{R_{{\rm{f}},i}}}} + \\ \dfrac{{G\left( {\frac{{v_{{\rm{m}},i + 1/2}^2}}{2}\! -\! \dfrac{{v_{{\rm{m}},i - 1/2}^2}}{2}} \right) \!+\! G{H_{\rm{p}}}({p_{{\rm{m}},i + 1/2}} - {p_{{\rm{m}},i - 1/2}}) + {{\dot q}_{{\rm{m}},i}}\Delta {z_i}}}{{G{c_{{\rm{pm}}}} \!+\! 1/{R_{{\rm{f}},i}}}} \end{array}$

    • 電加熱溫壓場耦合數值計算基本步驟為:

      1)依據井身結構、油管串及電加熱工藝進行網格劃分;

      2)確定網格中心對應的地層原始溫度剖面;

      3)初始化網格節點溫度、壓力和產出液密度與熱物性;

      4)依據地層水靜液柱壓力設置泵掛處的壓力;

      5)自下而上采用試算法迭代求解質量守恒方程,判別氣液兩相流流型,求解動量守恒方程和能量守恒方程,計算每個網格內中心節點溫度和下游邊界處的壓力,直至井口;

      6)比較計算的井口壓力與設定井口回壓,若滿足收斂條件,計算結束,輸出計算結果;否則,調整泵掛處的壓力,重復步驟5)和6)。

    • 大港油田X井為生產井,井身結構見表1,泵掛深度為1 300.00 m,油管外徑73.0 mm,日產油量5.37 m3,日產氣量53.00 m3,日產水量11.40 m3,溫度50 ℃時原油黏度911 mPa·s。該井采用連續電加熱工藝,加熱功率40 kW,加熱深度1 300 m,連續電加熱7 d后的平均井口溫度為59.60 ℃,應用上述井筒溫度場計算模型計算的井口溫度為61.45 ℃,較實測值略高,但相對誤差為3.10%,滿足工程精度要求,表明建立的井筒溫度場計算模型具有較好的可靠性。

      套管外徑/mm井眼直徑/mm套管下深/m水泥返高/m
      表層套管244.50311.1290.00地面
      生產套管139.70215.91 388.00 865.00

      表 1  大港油田X井實鉆井身結構

      Table 1.  Actual casing program of Well X in Dagang Oilfield

    • 應用建立的電加熱井筒溫度場數值計算方法,計算了加熱功率分別為20,40,60,80和100 kW時的連續電加熱和電磁短節加熱(電磁短節長度為5 m,分別布置在井深400.00,600.00,800.00,1 000.00和1 200.00 m處)2種加熱工藝的溫度場剖面特征,分別見圖3圖4

      圖  3  連續電加熱工藝的井筒溫度場剖面

      Figure 3.  Wellbore temperature field profile of continuous electric heating

      圖  4  電磁短節加熱工藝的井筒溫度場剖面

      Figure 4.  Wellbore temperature field profile of electromagnetic nipple heating

      圖3圖4可以看出,2種電加熱工藝的井筒溫度場剖面具有以下共同特征:1)隨著加熱功率增加,溫度場剖面逐漸向右偏移,產出液溫度升高;2)下部井段溫度梯度較高,上部井段溫度梯度較低;3)隨加熱功率增加,產出液溫度最高點對應的井深逐漸上移。2種加熱工藝的井筒溫度場剖面特征的主要差異為:1)連續電加熱工藝的井筒溫度場剖面平滑連續,而電磁短節加熱工藝的井筒溫度場剖面呈鋸齒形,電磁短節部位產出液溫度明顯升高,其上部產出液的溫度則快速下降;2)連續電加熱工藝的井筒溫度場剖面均方差小于電磁短節加熱工藝(見表2),表明連續電加熱工藝的溫度場剖面更均勻。

      加熱方法不同加熱功率的井筒溫度場剖面均方差
      20 kW40 kW60 kW80 kW100 kW
      連續電加熱11.638.405.312.973.82
      電磁短節加熱14.4010.96 7.404.404.01

      表 2  連續電加熱和電磁短節加熱工藝的井筒溫度場剖面均方差比較

      Table 2.  Comparison of the mean square errors of wellbore temperature field profiles formed by continuous electric heating and electromagnetic nipple heating processes

    • 連續電加熱和電磁短節加熱工藝的井口溫度和平均溫度計算結果如圖5所示。

      圖  5  連續電加熱和電磁短節加熱工藝的井口溫度與平均溫度對比

      Figure 5.  Comparison on the wellhead temperatures and average temperatures of continuous electric heating and electromagnetic nipple heating processes

      圖5可以看出,連續電加熱和電磁短節加熱的井口溫度和平均溫度均隨著加熱功率增大呈線性升高,連續電加熱工藝的井口溫度略高于電磁短節加熱工藝,而電磁短節加熱工藝的平均溫度略高于連續電加熱工藝。

      計算結果表明,電磁短節加熱功率為100 kW時,X井多個井深處產出液的溫度超過100.00 ℃,最高溫度為111.83 ℃,而采用連續電加熱工藝時產出液的最高溫度為96.68 ℃。可見,連續電加熱和電磁短節加熱稠油熱采過程中,井下產出液的溫度會有較大的波動,從而會對油管、井下工具和儀器的安全使用造成不利影響。

    • 1)考慮溫度對稠油熱物性影響,建立了連續電加熱和電磁短節加熱工藝的井筒溫度場數值計算方法,計算實例表明,模型計算結果與實測值相對誤差為3.10%,滿足工程設計精度要求。

      2)連續電加熱和電磁短節加熱工藝的井筒溫度場剖面均表現出下部井段溫度梯度較高、而上部井段溫度梯度較低的特征,但連續電加熱的溫度場剖面平滑連續,電磁短節加熱工藝的溫度場剖面呈鋸齒形,且溫度波動較大。

      3)加熱功率相同條件下,連續電加熱工藝的井口溫度略高于電磁短節加熱工藝,而連續電加熱工藝的平均溫度低于電磁短節加熱工藝。

      4)稠油黏溫關系對稠油熱采井筒溫度場預測與作業參數確定影響較大,建議在進行稠油熱采溫度場分析前進行5個溫度點以上的黏度測試,以提高分析精度。

參考文獻 (17)

目錄

    /

    返回文章
    返回