CN111648759A - 一种页岩气水平井积液位置判断方法 - Google Patents
一种页岩气水平井积液位置判断方法 Download PDFInfo
- Publication number
- CN111648759A CN111648759A CN202010449572.3A CN202010449572A CN111648759A CN 111648759 A CN111648759 A CN 111648759A CN 202010449572 A CN202010449572 A CN 202010449572A CN 111648759 A CN111648759 A CN 111648759A
- Authority
- CN
- China
- Prior art keywords
- gas
- liquid
- horizontal well
- critical
- pressure
- 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.)
- Granted
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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 DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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-气相表观流速。
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来继续计算下一段压力分布,直到所有计算分割段计算完成为止。
作为本申请的一种技术方案,步骤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-拟临界压力。
6.根据权利要求1所述的页岩气水平井积液位置判断方法,其特征在于,在步骤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 true CN111648759A (zh) | 2020-09-11 |
CN111648759B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113468826A (zh) * | 2021-06-17 | 2021-10-01 | 西南石油大学 | 基于真实液膜分布的页岩气水平井临界携液流量预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120318524A1 (en) * | 2011-06-20 | 2012-12-20 | Lea Jr James F | Plunger lift slug controller |
US20150060073A1 (en) * | 2013-09-05 | 2015-03-05 | Saudi Arabian Oil Company | Method of using concentrated solar power (csp) for thermal gas well deliquification |
CN106401570A (zh) * | 2015-07-30 | 2017-02-15 | 中国石油化工股份有限公司 | 页岩气井产水的确定方法、积液的确定方法及排液方法 |
CN107446566A (zh) * | 2016-05-31 | 2017-12-08 | 中国石油化工股份有限公司 | 用于深层页岩气水平井压裂的前置酸液及酸预处理方法 |
CN206860155U (zh) * | 2017-04-24 | 2018-01-09 | 贵州省地质矿产勘查开发局一一五地质大队 | 一种煤系气u型井钻进及压裂结构 |
-
2020
- 2020-05-25 CN CN202010449572.3A patent/CN111648759B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120318524A1 (en) * | 2011-06-20 | 2012-12-20 | Lea Jr James F | Plunger lift slug controller |
US20150060073A1 (en) * | 2013-09-05 | 2015-03-05 | Saudi Arabian Oil Company | Method of using concentrated solar power (csp) for thermal gas well deliquification |
CN106401570A (zh) * | 2015-07-30 | 2017-02-15 | 中国石油化工股份有限公司 | 页岩气井产水的确定方法、积液的确定方法及排液方法 |
CN107446566A (zh) * | 2016-05-31 | 2017-12-08 | 中国石油化工股份有限公司 | 用于深层页岩气水平井压裂的前置酸液及酸预处理方法 |
CN206860155U (zh) * | 2017-04-24 | 2018-01-09 | 贵州省地质矿产勘查开发局一一五地质大队 | 一种煤系气u型井钻进及压裂结构 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113468826A (zh) * | 2021-06-17 | 2021-10-01 | 西南石油大学 | 基于真实液膜分布的页岩气水平井临界携液流量预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111648759B (zh) | 2023-06-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110866321B (zh) | 一种适用于油气水同产水平井积液判别方法 | |
CN110878688B (zh) | 一种高气液比水平井临界携液流量的计算方法 | |
CN111648759A (zh) | 一种页岩气水平井积液位置判断方法 | |
CN111400978A (zh) | 一种考虑液滴形变和多参数影响的临界携液流量计算方法 | |
CN110580656B (zh) | 一种水平井筒地下极限携液流量预测方法 | |
US11099584B2 (en) | Method and apparatus for stabilizing gas/liquid flow in a vertical conduit | |
CN111779468B (zh) | 一种水平井斜井段涡流排水采气装置及排水采气方法 | |
CN111506864B (zh) | 一种分簇射孔中控制工具串与电缆作业安全的方法 | |
CN116341423A (zh) | 油水两相流滑移速度模型的计算方法 | |
CN115204504A (zh) | 一种衰竭式开采自喷油井停喷预测方法 | |
CN115936258A (zh) | 一种井筒结垢动态沉积堵塞预测模型构建方法 | |
Croce et al. | Intermittent Gas Lift for Liquid Loaded Horizontal Wells in Tight Gas Shale Reservoirs | |
US11649704B2 (en) | Processes and systems for injection of a liquid and gas mixture into a well | |
CN114462658A (zh) | 一种页岩气井积液判断方法 | |
Singh et al. | Unconventional cyclone gas lift completion for offshore wells of Cambay Basin: A smart completion to optimize production and well intervention | |
CN117332723B (zh) | 一种页岩气水平井井筒压降计算方法 | |
CN112031741B (zh) | 东海地区大斜度气井连续临界携液流速确定方法 | |
US20120211238A1 (en) | Gas production using a pump and dip tube | |
CN113688533B (zh) | 一种预测固井注水泥环空当量循环密度的计算方法 | |
CN115708102A (zh) | 一种新型泵口流入动态分析方法 | |
CN110863812A (zh) | 一种基于油套分离比的油井产气量计量方法 | |
Dye et al. | Flow Regulation in Horizontal Wells: Evaluating a Tailpipe System Designed to Optimize Artificial Lift Performance in Horizontal Wells | |
CN113516326A (zh) | 一种页岩气水平井临界排液能力评价方法 | |
CN118008212A (zh) | 一种页岩气水平井生产管柱下入优化方法 | |
Zadnorouzi et al. | Effective Parameters on Oil Production Rate using Artificial Gas Lift: Case Study in the One of South Iranian Oil Fields |
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 |