CN110197040B - 一种基于雷诺数的环空压力计算方法 - Google Patents
一种基于雷诺数的环空压力计算方法 Download PDFInfo
- Publication number
- CN110197040B CN110197040B CN201910493699.2A CN201910493699A CN110197040B CN 110197040 B CN110197040 B CN 110197040B CN 201910493699 A CN201910493699 A CN 201910493699A CN 110197040 B CN110197040 B CN 110197040B
- Authority
- CN
- China
- Prior art keywords
- annular
- pressure
- reynolds number
- gas
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Fluid Pressure (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及的是一种基于雷诺数的环空压力计算方法,它包括:一、通过现场测试或实验收集注采井的基本参数;二、计算注采井环空保护液下部泄露流量;三、计算窜流气体脱离水泥环表面时的初生气泡半径;四、判断气泡雷诺数,根据雷诺数的大小选用不同气泡上升速度模型;五、计算环空保护液中气体上升速度;步骤六:计算环空带压值,环空压力的计算是一个迭代过程,在得知某一时间气泡上升速度后,求解此时环空压力;在求得初始时刻i=1时的环空压力后,若P f =P c ,计算结束;若P f ≠P c ,i=i+1,返回步骤二重算,直至P f =P c 。本发明考虑了雷诺数对气泡上升速度的影响,计算出的环空压力值更能科学的反映环空压力实际上升情况。
Description
技术领域:
本发明涉及的是井筒完整性评价技术,具体涉及的是一种基于雷诺数的环空压力计算方法。
背景技术:
环空带压不仅会影响注采井的注采效果,同时还会严重影响生产安全。由于气体的性质要更为活泼,因此气井的环空带压现象更为普遍。据统计,美国墨西哥湾外大陆架有超过43%的井存在环空带压现象;中国塔里木油田有超过93%的高压气井环空带压,最高环空带压值甚至超过了50MPa。环空带压对注采井安全生产提出了严峻的挑战,因此为治理环空带压现象,有必要采用精确的环空压力计算方法,实现压力值的精确预测,从而指导环空带压现象的治理,保障注采井井筒完整性。
现有环空压力计算方法主要分为以下两类:
(1)忽略气体在环空保护液中的运移变化,仅考虑气体在水泥环中窜流及在气顶处聚集。
该类计算方法认为水泥环是质地均匀的多孔介质,利用达西公式求得从水泥环进入环空保护液的泄露流量,并将此泄露流量当作最终进入气顶的气体流量,而后通过PVT状态方程迭代计算得出环空压力值。该类方法可计算不同时间下环空压力上升情况,但窜流气体进入环空保护液后,气体在环空保护液中运移状态将发生变化,这将导致环空保护液下部泄露流量与上部进入气顶的气体流量产生差异。由于该类计算方法忽略气体在环空保护液中的运移变化,由此所得到单位时间内进入气顶的气体体积与实际情况不符,最终所计算出的环空压力值无法反映环空的真实压力,该类计算方法仅适用于固井水泥浆返至地面的情况。
(2)考虑窜流气体在环空保护液中的运移状态的变化。
该类方法能够计算水泥浆未返至地面情况下的环空压力。该类计算方法为描述气体在环空保护液中运移速度的变化规律,假定窜流气体在环空保护液中的上浮属于泡状流,通过对单一气泡受力分析建立平衡方程,基于该方程确定了最终气泡进入气顶的上升速度。相比于不考虑气体在环空保护液中运移变化的计算方法,该类计算方法更为符合井下实际特征。然而,该类计算方法仍存在不足:
首先,不同工况下窜流气体的气泡雷诺数将产生差异。当气泡雷诺数小于1时,气泡以球状上升,且在上升过程中形貌几乎不发生改变;其次,当气泡雷诺数大于1时,气泡在上升过程中极易由球型变为帽型。由于上升过程中气泡的形状直接影响运移速度,该类计算方法无法对不同雷诺数气泡的上升速度进行区分计算,因此现有的计算方法无法对环空压力上升规律进行精确预测。
注采井服役过程中环空带压不可避免,为了能够真实评价套管环空带压程度,需要提供一种与气体在环空保护内流动规律更为接近的计算方法。
发明内容:
本发明的目的是提供一种基于雷诺数的环空压力计算方法,这种基于雷诺数的环空压力计算方法用于解决现有的计算方法无法对环空压力上升规律进行精确预测问题。
本发明解决其技术问题所采用的技术方案是:这种基于雷诺数的环空压力计算方法:
步骤一:通过现场测试或实验收集注采井的基本参数,包括:地层的温度压力、水泥环孔隙度、水泥环渗透率、井口和井筒温度、环空保护液的相关物性参数和窜流气体的相关物性参数;
步骤二:计算注采井环空保护液下部泄露流量qg:
式中:K为水泥环渗透率,m2;A为环空截面积,m2;Tsc为标况温度,K;μg为气相粘度,Pa·s;Tf为地层温度,K;Lf为水泥环长度,m;Psc为标况压力,Pa;Pf为地层压力,Pa;Pc为环空保护液下部面压力,Pa。
步骤三:计算窜流气体脱离水泥环表面时的初生气泡半径rb:
式中:ρl为环空保护液密度,kg/m3;ρg为窜流气体密度,kg/m3;g为重力加速度,9.8m/s2;μl为环空保护液粘度,Pa·s;σ为环空保护液与气体的表面张力系数,N/m;φ为水泥环孔隙度,无量纲;为水泥颗粒直径,m。
步骤四:判断气泡雷诺数,根据雷诺数的大小选用不同气泡上升速度模型;气泡雷诺数公式为:
①当气泡雷诺数Re<1时,选用小雷诺数气泡上升速度计算公式:
式中:CD为阻力系数,其值为24/Re。
②当气泡雷诺数Re>1时,选用大雷诺数气泡上升速度计算公式:
步骤五:计算环空保护液中气体上升速度vt,w:
其中:
a3=-(a1+a2+a4)
步骤六:计算环空带压值Pt:
式中:Vm为环空保护液体积,m3;Vt为气顶体积,m3;Cm为环空保护液压缩系数,Pa-1;Twh为井口温度,K;A为环空截面积,m2;Twb为井筒温度,K;Lf为环空保护液高度,m。
环空压力的计算是一个迭代过程,在得知某一时间气泡上升速度后,求解此时环空压力;在求得初始时刻i=1时的环空压力Pt 1后,若Pf=Pc,则计算结束;若Pf≠Pc,则i=i+1,返回步骤二重新计算,直至地层压力Pf与环空保护液下部压力Pc相等。
本发明具有以下有益效果:
(1)本发明考虑了雷诺数对气泡上升速度的影响,以窜流的气泡雷诺数大小为评价指标,所计算出的环空压力值更能科学的反映环空压力实际上升情况;
(2)本发明可对影响环空带压的众多参数进行敏感性分析,并基于不同条件下环空压力值的变化规律,实现对工程参数的优选;
(3)本发明提出的环空压力计算方法具有普遍适用性,且能够依据环空压力动态曲线,确定合理的放压时机,从而为注采井井筒完整性提供保障。
具体实施方式
下面对本发明做进一步的说明:
这种基于雷诺数的环空压力计算方法:
步骤一:通过现场测试或实验收集注采井的基本参数,包括:地层的温度压力、水泥环孔隙度、水泥环渗透率、井口和井筒温度、环空保护液的相关物性参数和窜流气体的相关物性参数;
步骤二:计算注采井环空保护液下部泄露流量qg:
式中:K为水泥环渗透率,m2;A为环空截面积,m2;Tsc为标况温度,K;μg为气相粘度,Pa·s;Tf为地层温度,K;Lf为水泥环长度,m;Psc为标况压力,Pa;Pf为地层压力,Pa;Pc为环空保护液下部面压力,Pa。
步骤三:计算窜流气体脱离水泥环表面时的初生气泡半径rb:
式中:ρl为环空保护液密度,kg/m3;ρg为窜流气体密度,kg/m3;g为重力加速度,9.8m/s2;μl为环空保护液粘度,Pa·s;σ为环空保护液与气体的表面张力系数,N/m;φ为水泥环孔隙度,无量纲;为水泥颗粒直径,m。
步骤四:判断气泡雷诺数,根据雷诺数的大小选用不同气泡上升速度模型;气泡雷诺数公式为:
①当气泡雷诺数Re<1时,选用小雷诺数气泡上升速度计算公式:
式中:CD为阻力系数,其值为24/Re。
②当气泡雷诺数Re>1时,选用大雷诺数气泡上升速度计算公式:
步骤五:计算环空保护液中气体上升速度vt,w:
其中:
a3=-(a1+a2+a4)
式中:α为含气率,无量纲;n为流性指数,无量纲。
步骤六:计算环空带压值Pt:
式中:Vm为环空保护液体积,m3;Vt为气顶体积,m3;Cm为环空保护液压缩系数,Pa-1;Twh为井口温度,K;A为环空截面积,m2;Twb为井筒温度,K;Lf为环空保护液高度,m。
环空压力的计算是一个迭代过程,在得知某一时间气泡上升速度后,求解此时环空压力;在求得初始时刻i=1时的环空压力Pt 1后,若Pf=Pc,则计算结束;若Pf≠Pc,则i=i+1,返回步骤二重新计算,直至地层压力Pf与环空保护液下部压力Pc相等。
本发明分别计算不同雷诺数条件下的气泡运移速度,更为符合环空保护液内气体的真实运移规律,准确描述环空压力恢复期间各时间段内压力上升情况,从而指导环空带压诊断及井筒完整性补救措施的制定,延缓环空带压对安全生产造成的危害。
Claims (1)
1.一种基于雷诺数的环空压力计算方法,其特征在于:
步骤一:通过现场测试或实验收集注采井的基本参数,包括:地层的温度压力、水泥环孔隙度、水泥环渗透率、井口和井筒温度、环空保护液的相关物性参数和窜流气体的相关物性参数;
步骤二:计算注采井环空保护液下部泄露流量qg:
式中:K为水泥环渗透率,m2;A为环空截面积,m2;Tsc为标况温度,K;μg为气相粘度,Pa·s;Tf为地层温度,K;Lf为水泥环长度,m;Psc为标况压力,Pa;Pf为地层压力,Pa;Pc为环空保护液下部面压力,Pa;
步骤三:计算窜流气体脱离水泥环表面时的初生气泡半径rb:
式中:ρl为环空保护液密度,kg/m3;ρg为窜流气体密度,kg/m3;g为重力加速度,9.8m/s2;μl为环空保护液粘度,Pa·s;σ为环空保护液与气体的表面张力系数,N/m;φ为水泥环孔隙度,无量纲;为水泥颗粒直径,m;
步骤四:判断气泡雷诺数,根据雷诺数的大小选用不同气泡上升速度模型;气泡雷诺数公式为:
①当气泡雷诺数Re<1时,选用小雷诺数气泡上升速度计算公式:
式中:CD为阻力系数,其值为24/Re;
②当气泡雷诺数Re>1时,选用大雷诺数气泡上升速度计算公式:
步骤五:计算环空保护液中气体上升速度vt,w:
其中:
a3=-(a1+a2+a4)
式中:α为含气率,无量纲;n为流性指数,无量纲;
步骤六:计算环空带压值Pt:
式中:Vm为环空保护液体积,m3;Vt为气顶体积,m3;Cm为环空保护液压缩系数,Pa-1;Twh为井口温度,K;A为环空截面积,m2;Twb为井筒温度,K;Lf为环空保护液高度,m;
环空压力的计算是一个迭代过程,在得知某一时间气泡上升速度后,求解此时环空压力;在求得初始时刻i=1时的环空压力Pt 1后,若Pf=Pc,则计算结束;若Pf≠Pc,则i=i+1,返回步骤二重新计算,直至地层压力Pf与环空保护液下部压力Pc相等。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910493699.2A CN110197040B (zh) | 2019-06-06 | 2019-06-06 | 一种基于雷诺数的环空压力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910493699.2A CN110197040B (zh) | 2019-06-06 | 2019-06-06 | 一种基于雷诺数的环空压力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110197040A CN110197040A (zh) | 2019-09-03 |
CN110197040B true CN110197040B (zh) | 2023-04-07 |
Family
ID=67754149
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910493699.2A Active CN110197040B (zh) | 2019-06-06 | 2019-06-06 | 一种基于雷诺数的环空压力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110197040B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498214A (zh) * | 2009-03-06 | 2009-08-05 | 中国海洋石油总公司 | 一种偏心圆管环空幂律流体流态的测算方法 |
CN202628028U (zh) * | 2012-07-10 | 2012-12-26 | 东北石油大学 | 水平井剖面顶替效率测定装置 |
CN103195394A (zh) * | 2013-04-15 | 2013-07-10 | 中国石油大学(华东) | 增加避台风期间海洋钻井安全作业周期的方法 |
CN104453813A (zh) * | 2014-11-19 | 2015-03-25 | 东北石油大学 | 一种利用毛细管单向阀控制起泡剂加注装置及其方法 |
CN204255802U (zh) * | 2014-12-17 | 2015-04-08 | 中国科学院工程热物理研究所 | 液体参数测量系统 |
CA3003510A1 (en) * | 2015-10-29 | 2017-05-04 | Nexen Energy Ulc | Emulsion composition sensor |
CN107060731A (zh) * | 2017-07-04 | 2017-08-18 | 中国石油集团钻井工程技术研究院 | 一种基于井涌余量的深水钻井套管下深修正方法 |
CN107575214A (zh) * | 2016-07-04 | 2018-01-12 | 中国石油天然气股份有限公司 | 用于注采过程的井筒内温度与压力的预测方法 |
CN108131130A (zh) * | 2017-12-05 | 2018-06-08 | 中国石油天然气集团公司 | 对气井井口套管环空压力监测数据的分析方法及其装置 |
CN208140203U (zh) * | 2018-05-08 | 2018-11-23 | 东北石油大学 | 一种基于改进温度补偿方法的气体质量流量测量装置 |
CN109162697A (zh) * | 2018-10-13 | 2019-01-08 | 中石化石油工程技术服务有限公司 | 一种天然气井油套环空流动动态描述方法 |
CN109815513A (zh) * | 2017-11-22 | 2019-05-28 | 四川省雅安市羌江机械有限责任公司 | 一种航空发动机滑油系统设计方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7027968B2 (en) * | 2002-01-18 | 2006-04-11 | Conocophillips Company | Method for simulating subsea mudlift drilling and well control operations |
US8577626B2 (en) * | 2008-07-22 | 2013-11-05 | General Electric Company | System and method for assessing fluid dynamics |
US8620611B2 (en) * | 2009-08-13 | 2013-12-31 | Baker Hughes Incorporated | Method of measuring multi-phase fluid flow downhole |
MY149416A (en) * | 2009-11-17 | 2013-08-30 | Petroliam Nasional Berhad Petronas | Enhanced dynamic well model for reservoir pressure determination |
-
2019
- 2019-06-06 CN CN201910493699.2A patent/CN110197040B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101498214A (zh) * | 2009-03-06 | 2009-08-05 | 中国海洋石油总公司 | 一种偏心圆管环空幂律流体流态的测算方法 |
CN202628028U (zh) * | 2012-07-10 | 2012-12-26 | 东北石油大学 | 水平井剖面顶替效率测定装置 |
CN103195394A (zh) * | 2013-04-15 | 2013-07-10 | 中国石油大学(华东) | 增加避台风期间海洋钻井安全作业周期的方法 |
CN104453813A (zh) * | 2014-11-19 | 2015-03-25 | 东北石油大学 | 一种利用毛细管单向阀控制起泡剂加注装置及其方法 |
CN204255802U (zh) * | 2014-12-17 | 2015-04-08 | 中国科学院工程热物理研究所 | 液体参数测量系统 |
CA3003510A1 (en) * | 2015-10-29 | 2017-05-04 | Nexen Energy Ulc | Emulsion composition sensor |
CN107575214A (zh) * | 2016-07-04 | 2018-01-12 | 中国石油天然气股份有限公司 | 用于注采过程的井筒内温度与压力的预测方法 |
CN107060731A (zh) * | 2017-07-04 | 2017-08-18 | 中国石油集团钻井工程技术研究院 | 一种基于井涌余量的深水钻井套管下深修正方法 |
CN109815513A (zh) * | 2017-11-22 | 2019-05-28 | 四川省雅安市羌江机械有限责任公司 | 一种航空发动机滑油系统设计方法 |
CN108131130A (zh) * | 2017-12-05 | 2018-06-08 | 中国石油天然气集团公司 | 对气井井口套管环空压力监测数据的分析方法及其装置 |
CN208140203U (zh) * | 2018-05-08 | 2018-11-23 | 东北石油大学 | 一种基于改进温度补偿方法的气体质量流量测量装置 |
CN109162697A (zh) * | 2018-10-13 | 2019-01-08 | 中石化石油工程技术服务有限公司 | 一种天然气井油套环空流动动态描述方法 |
Non-Patent Citations (9)
Title |
---|
CO_2注入井封隔器早期密封失效机制研究;张德平,冯福平,李清,严茂森,艾池,丛子渊.;《润滑与密封》;第44卷(第5期);全文 * |
Influence of various aspects of low Reynolds number k-epsilon turbulence models on predicting in-tube buoyancy affected heat transfer to supercritical pressure fluids;Zhao, CR;Zhang, Z;Jiang, PX;《NUCLEAR ENGINEERING AND DESIGN》;第313卷;全文 * |
Interfacial and wall friction factors of swirling annular flow in a vertical pipe;Funahashi, H;Kirkland, KV;Hayashi, K;《NUCLEAR ENGINEERING AND DESIGN》;第330卷;全文 * |
Laminar and turbulent comparisons for channel flow and flow control;Marusic, Ivan;Joseph, D. D;Mahesh, Krishnan;《JOURNAL OF FLUID MECHANICS》;第570卷;全文 * |
Prediction of natural gas hydrate formation region in wellbore during deep-water gas well testing;Wang, ZY;Sun, BJ;Wang, XR;《JOURNAL OF HYDRODYNAMICS》;第26卷(第4期);全文 * |
低雷诺数下层流分离的等离子体控制;孟宣市,杨泽人,陈琦;《航空学报》;第37卷(第7期);全文 * |
循环流化床提升管内压力脉动特性;胡小康,刘小成,徐俊;《化工学报》;第61卷(第4期);全文 * |
气井井筒气液两相环雾流压降计算新方法;刘通,王世泽,郭新江,潘国华,陈海龙,伍洲.;《石油钻采工艺》;第39卷(第3期);全文 * |
深水气井测试临界携液条件的优化设计;何玉发,李紫晗,张滨海,高飞,李莹莹.;《天然气工业》;第37卷(第9期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110197040A (zh) | 2019-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shuhong et al. | New understanding on water-oil displacement efficiency in a high water-cut stage | |
CN107301306B (zh) | 用于致密砂岩气藏压裂水平井的动态无阻流量预测方法 | |
CN106285646B (zh) | 基于多信息融合的钻井漏失层位识别方法 | |
CN104504604B (zh) | 一种定性气井井筒积液的方法 | |
CN107462936B (zh) | 利用压力监测资料反演低渗透储层非达西渗流规律的方法 | |
CN112343587A (zh) | 一种特低渗透油藏优势渗流通道识别表征方法 | |
CN109001438A (zh) | 一种封缝堵气实验模拟装置及测试方法 | |
CN106522928B (zh) | 一种酸化压裂后停泵测井口压降不稳定试井方法 | |
CN109061766B (zh) | 一种断-砂配置油气侧向分流能力的定量评价方法 | |
CN109915128B (zh) | 地层承压能力动态测试方法及固井方法 | |
CN107448197B (zh) | 一种海相砂岩油藏窜流通道空间发育定量识别的方法 | |
CN105672997A (zh) | 钻井液地层漏失量监测方法 | |
CN105019875B (zh) | 人工隔层隔离剂评价方法 | |
CN109459795A (zh) | 一种孔隙性砂岩断裂带封闭性评价方法 | |
CN114370269B (zh) | 深层碳酸盐岩气藏有效储层物性下限综合确定方法 | |
CN110197040B (zh) | 一种基于雷诺数的环空压力计算方法 | |
CN105547960B (zh) | 一种基于透明砂土的基坑降水地下水渗流可视化模拟试验方法 | |
CN107503717B (zh) | 一种砾石充填防砂储层保护效果的评价方法 | |
CN106321076B (zh) | 一种注水井启动压力测试方法 | |
CN113109162A (zh) | 一种基于热流固耦合的岩石起裂压力计算方法 | |
CN109932296B (zh) | 一种定量表征贾敏效应动态变化的方法 | |
CN110593856B (zh) | 一种固井安全作业密度窗口测定方法 | |
CN106484930A (zh) | 用于确定非均质油藏复杂结构井产量的方法及系统 | |
CN110929447A (zh) | 一种稠化酸酸化过程中井筒温度场数值计算方法 | |
CN114169204B (zh) | 一种用于海上油气田开发生产的防砂时机确定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |