CN111648759B - 一种页岩气水平井积液位置判断方法 - Google Patents
一种页岩气水平井积液位置判断方法 Download PDFInfo
- Publication number
- CN111648759B CN111648759B CN202010449572.3A CN202010449572A CN111648759B CN 111648759 B CN111648759 B CN 111648759B CN 202010449572 A CN202010449572 A CN 202010449572A CN 111648759 B CN111648759 B CN 111648759B
- Authority
- CN
- China
- Prior art keywords
- gas
- critical
- pressure
- horizontal well
- calculating
- 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
- 239000007788 liquid Substances 0.000 title claims abstract description 63
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000009825 accumulation Methods 0.000 title description 4
- 239000007789 gas Substances 0.000 claims abstract description 105
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 34
- 239000003345 natural gas Substances 0.000 claims abstract description 21
- 239000000203 mixture Substances 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000001133 acceleration Effects 0.000 claims description 15
- 239000012530 fluid Substances 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 238000004519 manufacturing process Methods 0.000 claims description 8
- 239000012071 phase Substances 0.000 claims description 6
- 239000007791 liquid phase Substances 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/06—Measuring temperature or pressure
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/10—Geothermal energy
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本申请提供一种页岩气水平井积液位置判断方法,涉及页岩气领域;其包括以下步骤:S1.计算页岩气水平井沿井筒分布的压力温度剖面,判断沿所述井筒方向的流型分布,计算所述井筒不同位置处的气体流速;S2.计算天然气偏差系数随压力温度的变化,已知天然气组分组成,计算出天然气拟临界压力、拟临界温度;S3.计算页岩气水平井的水平段、倾斜段以及垂直段三处的临界携液流量;S4.计算临界携液流量并对比不同位置处的气体真实流量,根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图,根据所述对比变化图判断出积液的位置。该方法建立了一套系统的判断水平井中的积液位置的方法,提高了判断积液位置的准确度。
Description
技术领域
本申请涉及页岩气领域,尤其涉及一种页岩气水平井积液位置判断方法。
背景技术
页岩气是一种储量大、开采难的非常规油气资源,近年来,国内长水平井分段压裂理论、技术和装备的快速发展、成熟和完善,页岩气实现了工业化、商业化的开发。而页岩气一般具有初期产量高、中后期产量递减快、返排率低的特点。压裂液长期同气体从地层到井底再到井口,需要消耗更多的地层能量,后期地层压力不足时,将会有液体停滞在井底或水平段,对页岩气的排采稳产工作造成技术难题,尤其是积液位置的不确定性,对排采工艺的选择无指导性。
国内外页岩气开发时间较短,后期排采理论与技术一般均是借鉴常规天然气手段,无法同页岩气长水平段的特殊性建立联系,目前尚未见到页岩气水平井积液位置判断方法的公开报道。
发明内容
本申请提供一种页岩气水平井积液位置判断方法,解决现有技术中页岩气水平井中无法准确且系统的判断出积液的位置问题。
本申请的技术方案是:
一种页岩气水平井积液位置判断方法,包括以下步骤:
S1.计算页岩气水平井沿井筒分布的压力温度剖面,判断沿所述井筒方向的流型分布,计算所述井筒不同位置处的气体流速;
S2.计算天然气偏差系数随压力温度的变化,已知天然气组分组成,计算出天然气拟临界压力、拟临界温度;
S3.计算页岩气水平井的水平段、倾斜段以及垂直段三处的临界携液流量;
S4.计算临界携液流量并对比不同位置处的气体真实流量,根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图,根据所述对比变化图判断出积液的位置;若是对比变化图中的临界携液气量大于实际产气量,则所述页岩气水平井段中有积液;若是对比变化图中的临界携液气量小于实际产气量,则所述页岩气水平井段中没有积液。
作为本申请的一种技术方案,在所述步骤S1中,采用Beggs-Brill方法计算页岩气水平井沿井筒分布的压力温度剖面,所述Beggs-Brill方法包括以下公式:
式中:p-压力;z-沿井筒方向的长度;-压力梯度;ρl-液体密度;Hl-持液率;ρg-气体密度;g-重力加速度;θ-管柱与水平方向的夹角;λ-流动阻力系数;G-混合物的质量流量;v-混合物平均流速;D-油管内径;A-管的流通截面积;Vsg-气相表观流速。
作为本申请的一种技术方案,在所述步骤S1中,包括以下步骤:
S11.采用Beggs-Brill方法的公式,根据基础数据,确定计算的起始点压力p、计算段数n、段深度△h;
S13.计算所述计算段内的压力差Δp′与末端压力ph′1,若满足ph′1-ph1<0.0001,则将计算的末端压力作为下一段的起点压力;若不满足ph′1-ph1<0.0001,则重新将ph′1替代ph1,继续计算,直到满足误差要求为止;
S14.重复步骤S11至S13来继续计算下一段压力分布,直到所有计算分割段计算完成为止。
作为本申请的一种技术方案,步骤S11中的所述基础数据包括液体密度、气体密度、液相体积流量、气相体积流量、液体表面张力、重力加速度、液体黏度、气体黏度、井筒倾斜角、油管内径、井口温度以及温度梯度。
作为本申请的一种技术方案,在步骤S2中,采用Dranchuk-Abu-Kassem方法计算天然气偏差系数随压力温度的变化,所述Dranchuk-Abu-Kassem方法包括以下公式:
采用牛顿迭代法对Z进行迭代计算,式中:Z-,A1=0.3265,A2=-1.0700,Tpr-拟临界温度,A3=-0.5339,A4=0.01569,A5=-0.05165,ρpr-无因次对比密度,A6=0.5457,A7=-0.7361,A8=0.1844,A9=0.1056,A10=0.6134,A11=0.7210,Ppr-拟临界压力。
作为本申请的一种技术方案,在步骤S3中,所述页岩气水平井的临界产气量为:
作为本申请的一种技术方案,在步骤S3中,所述页岩气水平井的垂直段临界速度为:
式中:g-重力加速度;σ-气水界面张力,N/cm;ρl-液体密度,kg/m3;ρg-气体密度,kg/m3。
作为本申请的一种技术方案,在步骤S3中,所述页岩气水平井的倾斜段临界速度为:
式中:g-重力加速度;σ-气水界面张力,N/cm;ρl-液体密度,kg/m3;ρg-气体密度,kg/m3;α-井斜角。
作为本申请的一种技术方案,在步骤S3中,所述页岩气水平井的水平井段临界速度为:
式中:g-重力加速度;σ-气水界面张力,N/cm;ρl-液体密度,kg/m3;ρg-气体密度,kg/m3。
本申请的有益效果:
本申请通过利用Dranchuk-Abu-Kassem方法计算天然气偏差系数随压力温度的变化,采用Beggs-Brill方法计算沿井筒分布的压力温度剖面,判断沿井筒方向的流型分布,计算不同位置处的气体流速,通过计算临界携液流量并对比不同位置出的气体真实流量,根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图,根据该对比变化图来系统有效的判断出页岩气水平井段中的积液的位置;该方法科学系统,能够有效地判断出页岩气水平井段的积液的位置,方法科学,且使用效率高。
具体实施方式
实施例:
本申请实施例提供了一种页岩气水平井积液位置判断方法,其主要包括以下步骤:
S1.采用Beggs-Brill方法计算页岩气水平井沿井筒分布的压力温度剖面,判断沿井筒方向的流型分布,计算井筒不同位置处的气体流速;
S2.采用Dranchuk-Abu-Kassem方法计算天然气偏差系数随压力温度的变化,已知天然气组分组成,计算出天然气拟临界压力、拟临界温度;
S3.计算页岩气水平井的水平段、倾斜段以及垂直段三处的临界携液流量;
S4.计算临界携液流量并对比不同位置处的气体真实流量,根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图,根据对比变化图判断出积液的位置;若是对比变化图中的临界携液气量大于实际产气量,则页岩气水平井段中有积液;若是对比变化图中的临界携液气量小于实际产气量,则页岩气水平井段中没有积液。
需要说明的是,在本实施例中,Beggs-Brill方法包括以下公式:
式中:p-压力;z-沿井筒方向的长度;-压力梯度;ρl-液体密度;Hl-持液率;ρg-气体密度;g-重力加速度;θ-管柱与水平方向的夹角;λ-流动阻力系数;G-混合物的质量流量;v-混合物平均流速;D-油管内径;A-管的流通截面积;Vsg-气相表观流速。
需要说明的是,在本实施例中,Dranchuk-Abu-Kassem方法包括以下公式:
采用牛顿迭代法对Z进行迭代计算,式中:Z-,A1=0.3265,A2=-1.0700,Tpr-拟临界温度,A3=-0.5339,A4=0.01569,A5=-0.05165,ρpr-无因次对比密度,A6=0.5457,A7=-0.7361,A8=0.1844,A9=0.1056,A10=0.6134,A11=0.7210,Ppr-拟临界压力。需要说明的是,该方法适用于1.0≤Tpr≤3.0,或0.7≤Tpr≤1.0,ppr<1.0的情况。
进一步地,在步骤S1中,其主要包括以下步骤:
S11.采用Beggs-Brill方法的公式,根据基础数据,确定计算的起始点压力p、计算段数n、段深度△h;
S13.计算计算段内的压力差Δp′与末端压力p′h1,若满足|p′h1-ph1|<0.0001,则将计算的末端压力作为下一段的起点压力,若不满足|p′h1-ph1|<0.0001,则重新将p′h1替代ph1,继续计算直到满足误差要求为止;
S14.重复步骤S11至S13来继续计算下一段压力分布,直到所有计算分割段计算完成。
需要说明的是,在本实施例中,步骤S11中的基础数据包括液体密度、气体密度、液相体积流量、气相体积流量、液体表面张力、重力加速度、液体黏度、气体黏度、井筒倾斜角、油管内径、井口温度以及温度梯度。
需要说明的是,在本实施例中,在步骤S3中,页岩气水平井的临界产气量为:
需要说明的是,在本实施例中,在步骤S3中,页岩气水平井的垂直段临界速度为:
式中:g-重力加速度;σ-气水界面张力,N/cm;ρl-液体密度,kg/m3;ρg-气体密度,kg/m3。
需要说明的是,在本实施例中,在步骤S3中,页岩气水平井的倾斜段临界速度为:
式中:g-重力加速度;σ-气水界面张力,N/cm;ρl-液体密度,kg/m3;ρg-气体密度,kg/m3;α-井斜角。
需要说明的是,在本实施例中,在步骤S3中,页岩气水平井的水平井段临界速度为:
式中:g-重力加速度;σ-气水界面张力,N/cm;ρl-液体密度,kg/m3;ρg-气体密度,kg/m3。
需要说明的是,在本实施例中,本申请通过利用Dranchuk-Abu-Kassem方法计算天然气偏差系数随压力温度的变化,采用Beggs-Brill方法计算沿井筒分布的压力温度剖面,判断沿井筒方向的流型分布,计算不同位置处的气体流速,通过计算临界携液流量并对比不同位置出的气体真实流量,根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图,根据该对比变化图来系统有效的判断出页岩气水平井段中的积液的位置;该方法科学系统,能够有效地判断出页岩气水平井段的积液的位置,方法科学,且使用效率高。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,应包含在本申请的保护范围之内。
Claims (9)
1.一种页岩气水平井积液位置判断方法,其特征在于,包括以下步骤:
S1.计算页岩气水平井沿井筒分布的压力温度剖面,判断沿所述井筒方向的流型分布,计算所述井筒不同位置处的气体流速;
S2.计算天然气偏差系数随压力温度的变化,已知天然气组分组成,计算出天然气拟临界压力、拟临界温度;
S3.计算页岩气水平井的水平段、倾斜段以及垂直段三处的临界携液流量;
S4.通过计算的临界携液流量并对比页岩气水平井的不同位置处的气体真实流量,根据井眼轨迹分布绘制不同参数下产量与临界携液流量沿井眼轨迹对比变化图,根据所述对比变化图判断出积液的位置;若是对比变化图中的临界携液流量大于实际产气量,则所述页岩气水平井段中有积液;若是对比变化图中的临界携液流量小于实际产气量,则所述页岩气水平井段中没有积液。
3.根据权利要求2所述的页岩气水平井积液位置判断方法,其特征在于,在所述步骤S1中,包括以下步骤:
S11.采用Beggs-Brill方法的公式,根据基础数据,确定计算的起始点压力p、计算段数n、段深度△h;
S13.计算所述计算段内的压力差Δp′与末端压力p′h1,若满足|p′h1-ph1|<0.0001,则将计算的末端压力作为下一段的起点压力;若不满足|p′h1-ph1|<0.0001,则重新将p′h1替代ph1,继续计算,直到满足误差要求为止;
S14.重复步骤S11至S13来继续计算下一段压力分布,直到所有计算分割段计算完成为止。
4.根据权利要求3所述的页岩气水平井积液位置判断方法,其特征在于,步骤S11中的所述基础数据包括液体密度、气体密度、液相体积流量、气相体积流量、液体表面张力、重力加速度、液体黏度、气体黏度、井筒倾斜角、油管内径、井口温度以及温度梯度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010449572.3A CN111648759B (zh) | 2020-05-25 | 2020-05-25 | 一种页岩气水平井积液位置判断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010449572.3A CN111648759B (zh) | 2020-05-25 | 2020-05-25 | 一种页岩气水平井积液位置判断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111648759A CN111648759A (zh) | 2020-09-11 |
CN111648759B true CN111648759B (zh) | 2023-06-27 |
Family
ID=72352740
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010449572.3A Active CN111648759B (zh) | 2020-05-25 | 2020-05-25 | 一种页岩气水平井积液位置判断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111648759B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113468826B (zh) * | 2021-06-17 | 2024-05-28 | 西南石油大学 | 基于真实液膜分布的页岩气水平井临界携液流量预测方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9297247B2 (en) * | 2011-06-20 | 2016-03-29 | James F. Lea, Jr. | Plunger lift slug controller |
US9777562B2 (en) * | 2013-09-05 | 2017-10-03 | Saudi Arabian Oil Company | Method of using concentrated solar power (CSP) for thermal gas well deliquification |
CN106401570B (zh) * | 2015-07-30 | 2019-05-07 | 中国石油化工股份有限公司 | 页岩气井产水的确定方法、积液的确定方法及排液方法 |
CN107446566B (zh) * | 2016-05-31 | 2020-09-18 | 中国石油化工股份有限公司 | 用于深层页岩气水平井压裂的前置酸液及酸预处理方法 |
CN206860155U (zh) * | 2017-04-24 | 2018-01-09 | 贵州省地质矿产勘查开发局一一五地质大队 | 一种煤系气u型井钻进及压裂结构 |
-
2020
- 2020-05-25 CN CN202010449572.3A patent/CN111648759B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111648759A (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Takacs | Sucker-rod pumping handbook: production engineering fundamentals and long-stroke rod pumping | |
CN104895560B (zh) | 一种深水测试井筒压力、温度场模拟及水合物预测方法 | |
CN110866321B (zh) | 一种适用于油气水同产水平井积液判别方法 | |
Brito et al. | Effect of well trajectory on liquid removal in horizontal gas wells | |
CN111648759B (zh) | 一种页岩气水平井积液位置判断方法 | |
CN110878688A (zh) | 一种高气液比水平井临界携液流量的计算方法 | |
Nguyen | Artificial lift methods: design, practices, and applications | |
CN112347601A (zh) | 一种优化气井泡沫排水采气在线实时起泡剂加量的方法 | |
US11099584B2 (en) | Method and apparatus for stabilizing gas/liquid flow in a vertical conduit | |
CN110580656B (zh) | 一种水平井筒地下极限携液流量预测方法 | |
CN114396255A (zh) | 一种地下储气库气井的储层温度预测方法 | |
CN108756830B (zh) | 一种水平井砾石充填模拟计算方法 | |
CN102108851B (zh) | 水平井注汽方法及系统 | |
Dinata et al. | A methodology of end-of-tubing location optimization for horizontal shale gas wells with and without deliquification | |
CN111206919B (zh) | 长井段高产气井储层段井筒压力计算方法 | |
McCoy et al. | A laboratory study with field data of downhole gas separators | |
CN115204504A (zh) | 一种衰竭式开采自喷油井停喷预测方法 | |
Beauquin et al. | World's first metal PCP SAGD field test shows promising artificial-lift technology for heavy-oil hot production: Joslyn field case | |
US11859469B2 (en) | Utilizing natural gas flaring byproducts for liquid unloading in gas wells | |
Nair et al. | Existence and prediction of severe slugging in toe-down horizontal wells | |
CN109784003A (zh) | 一种基于水平井修井管柱的优化组合方法 | |
CN114462658A (zh) | 一种页岩气井积液判断方法 | |
Carpenter | Study Models Onset of Liquid Loading in Large-Diameter Deviated Gas Wells | |
CN110863812B (zh) | 一种基于油套分离比的油井产气量计量方法 | |
Dye et al. | Flow Regulation in Horizontal Wells: Evaluating a Tailpipe System Designed to Optimize Artificial Lift Performance in Horizontal Wells |
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 |