CN110984976B - 多层合采气井产量劈分方法 - Google Patents
多层合采气井产量劈分方法 Download PDFInfo
- Publication number
- CN110984976B CN110984976B CN201911325055.9A CN201911325055A CN110984976B CN 110984976 B CN110984976 B CN 110984976B CN 201911325055 A CN201911325055 A CN 201911325055A CN 110984976 B CN110984976 B CN 110984976B
- Authority
- CN
- China
- Prior art keywords
- gas
- scj
- zone
- production
- 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
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000007789 gas Substances 0.000 claims description 61
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 19
- 239000003345 natural gas Substances 0.000 claims description 15
- 230000000704 physical effect Effects 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 3
- 230000005251 gamma ray Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000011161 development Methods 0.000 abstract description 2
- XQCFHQBGMWUEMY-ZPUQHVIOSA-N Nitrovin Chemical compound C=1C=C([N+]([O-])=O)OC=1\C=C\C(=NNC(=N)N)\C=C\C1=CC=C([N+]([O-])=O)O1 XQCFHQBGMWUEMY-ZPUQHVIOSA-N 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005065 mining Methods 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling 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
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及油气田开发技术领域,具体涉及一种多层合采气井产量劈分方法。一种多层合采气井产量劈分方法,该过程如下:(1)根据气井标准状态下产气量Qsc和Kh值劈分得到产层j的产气量初始值;(2)得到产层j到产层n在标准状态下的气体体积流量之和;(3)计算气体流速初始值vj(0)下的雷诺数初始值Rej(0)和摩阻因子初始值fj(0);(4)计算得到产层j到产层n在标准状态下的气体体积流量之和;(5)计算产层j的产量,并进行误差判定。本发明根据实际井筒压力数据和气井管流方程进行产量劈分,理论依据更充分,计算精度更准确。
Description
技术领域
本发明涉及油气田开发技术领域,具体涉及一种多层合采气井产量劈分方法。
背景技术
为提高单井产量,很多气藏采用多层合采方式开采,各层的产出情况是气藏剩余气分布研究的理论依据,而不同产气量、不同时刻气井各层贡献情况也不同,因此产量劈分一直是气藏动态分析的难点。
目前产量劈分常用的方法有产出剖面测试、Kh值劈分法、数值模拟法等。产出剖面测试可以准确测量气井该时刻的产量贡献情况,但由于不同时刻各层产量贡献不同,需要定期对所有气井进行测试,操作成本高;Kh值劈分法主要研究静态参数对产量的影响,无法考虑产量贡献随时间的变化,且忽略了技术、管理等方面的影响,计算结果与实际生产差距较大;数值模拟方法可以得到不同时刻各层的贡献情况,但计算过程复杂,且具有多解性,应用效果差。
发明内容
本发明旨在针对上述问题,提出一种通过较易获得的井底流压测试数据,快速计算气井该时刻各小层的产量,从而实现对多层合采气井的产量劈分。
本发明的技术方案在于:
一种多层合采气井产量劈分方法,包括如下步骤:
(2)气井生产层位共有n层,根据气井标准状态下产气量Qsc和Kh值劈分得到产层j的产气量初始值qscj(0)
(3)计算得到产层j底部到产层j-1底部井筒段的气体流速初始值vj(0)
(4)计算气体流速初始值vj(0)下的雷诺数初始值Rej(0)和摩阻因子初始值fj(0)
式中,e为管壁绝对粗糙度,单位m;R为摩尔气体常数,8.315×10-3MPa·m3/(kmol·k);Re为雷诺数;D为生产管柱内径,单位m;γg为天然气相对密度;f为摩阻因子,无量纲;
式中,θ为气井管斜角,单位为°;sj为指数,无量纲;
(6)计算产层j的产量qscj′
(7)若|qscj′-qscj(0)|/qscj′≤ε,其中,ε=5%;则产层j的产量qscj即为qscj′;
若|qscj′-qscj(0)|/qscj>ε,则以qscj′作为qscj(1),重复步骤(3)-(6),直至计算出符合误差条件的qscj。
所述e取0.016×10-3m。
本发明的技术效果在于:
1.不测产气剖面的情况下,本发明通过井筒压力数据实现多层合采气井的产量劈分,大大节约了测试成本;
2.相对于Kh值劈分法,本发明可根据井筒压力测试数据得到气井不同阶段各层的产量贡献;
3.Kh值劈分法仅考虑了储层静态参数对产量的影响,数值模拟模型具有多解性,而本发明根据实际井筒压力数据和气井管流方程进行劈分产量,理论依据更充分,且计算精度更准确。
具体实施方式
具体实验例
以鄂尔多斯盆地东南部S265井为例,该井为山2、山1、盒8三层合采,日产气量Qsc=2.7610×104m3/d,产层处生产管柱为套管,套管内径D=0.1594m,管壁绝对粗糙度e=0.016×10-3m;三个产层底部垂深分别为H3=2675.0m、H2=2620.6m和H1=2569.8m,平均渗透率K3=0.15mD、K2=0.11mD、K1=0.13mD,有效厚度h3=15m、h2=5m、h3=10m。通过高精度压力计测试得到三个产层底部的井筒压力值分别为Pwf3=16.128776MPa、Pwf2=16.077504MPa和Pwf1=16.029593MPa,井筒温度分别为T3=94.60℃、T2=93.88℃和T1=93.14℃,天然气相对密度γg=0.5871。
利用本发明提供的多层合采气井产量劈分方法的具体实施过程如下。
一、计算产层3的产气量
(3)计算指数s3
(4)根据气井标准状态下产气量Qsc和Kh值劈分得到产层3产气量的初始值qsc3(0)
(5)计算产层3底部到产层2底部井筒段的气体流速初始值v3(0)
(6)计算气体流速初始值v3(0)下的雷诺数Re3(0)和摩阻因子f3(0)
(7)计算产层3在标准状态下的气体体积流量qsc3′为
误差判定:设定误差ε=5%;|qsc3′-qsc3(0)|/qsc3′=|1.3553-0.8754|/1.3553=35.41%>5%;不满足误差判定条件,进行第一次迭代,将qsc3′作为产层3的产量初始值qsc3(1),即qsc3(1)=qsc3′,重复步骤(5)-(7),计算得到v3(1)=0.057621m/s,Re3(1)=51690.08,f3(1)=0.0210,qsc3″=1.4214×104m3/d,再次判断|qsc3″-qsc3(1)|/qsc3″=|1.4214-1.3553|/1.4214=4.65%<5%,满足误差判定条件;因此产层3的产量qsc3=qsc3″=1.4214×104m3/d。
二、计算产层2的产气量
(3)计算指数s2
(4)根据气井标准状态下产气量Qsc和Kh值劈分得到产层2产气量的初始值qsc2(0)
(5)计算产层2底部到产层1底部井筒段的气体流速初始值v2(0)
(6)计算气体流速初始值v2(0)下的雷诺数Re2(0)和摩阻因子f2(0)
(7)计算产层2到产层3在标准状态下的气体体积流量之和为
(8)计算产层2的产量qsc2′
误差判定:设定误差ε=5%;|qsc2′-qsc2(0)|/qsc2′=|0.6831-0.3704|/0.6831=45.78%>5%;
不满足误差判定条件,进行第一次迭代:将qsc2′作为产层2的产量初始值qsc2(1),即qsc2(1)=qsc2′;作为产层2到产层3体积流量之和初始值即重复步骤(5)-(8),计算得到v2(1)=0.089571m/s,Re2(1)=80263.98,f2(1)=0.0192,
qsc2″=0.7933×104m3/d,再次判断|qsc2″-qsc2(1)|/qsc2″=|0.7933-0.6831|/0.7933=13.89%>5%,不满足误差判定条件,继续迭代。将qsc2″作为产层2的产量初始值qsc2(2),即qsc2(2)=qsc2″,作为产层2到产层3体积流量之和初始值即重复步骤(5)-(8),计算得到
v2(2)=0.094261m/s,Re2(2)=84466.93,f2(2)=0.0190,qsc2″′=0.8049×104m3/d,再次判断|qsc2″′-qsc2(2)|/qsc2″′=|0.8049-0.7933|/0.8049=1.44%<5%,因此产层2的产量qsc2=qsc2″′=0.8049×104m3/d。
三、计算产层1的产气量
四、结果分析
对比发现,本发明计算结果与产气剖面测试误差较小,计算精度远大于Kh值劈产法,见表1。
表1本发明比Kh值劈产法对比
Claims (5)
1.一种多层合采气井产量劈分方法,其特征在于:包括如下步骤:
(2)气井生产层位共有n层,根据气井标准状态下产气量Qsc和Kh值劈分得到产层j的产气量初始值qscj(0)
(3)计算得到产层j底部到产层j-1底部井筒段的气体流速初始值vj(0)
(4)计算气体流速初始值vj(0)下的雷诺数初始值Rej(0)和摩阻因子初始值fj(0)
式中,e为管壁绝对粗糙度,单位m;R为摩尔气体常数,8.315×10-3MPa·m3/(kmol·k);Re为雷诺数;D为生产管柱内径,单位m;γg为天然气相对密度;f为摩阻因子,无量纲;
式中,θ为气井管斜角,单位为°;sj为指数,无量纲;
(6)计算产层j的产量qscj′
(7)若|qscj′-qscj(0)|/qscj′≤ε,其中,ε=5%;则产层j的产量qscj即为qscj′。
2.根据权利要求1所述多层合采气井产量劈分方法,其特征在于:若|qscj′-qscj(0)|/qscj>ε,则以qscj′作为qscj(1),重复步骤(3)-(6),直至计算出符合误差条件的qscj。
5.根据权利要求4所述多层合采气井产量劈分方法,其特征在于:所述e取0.016×10- 3m。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911325055.9A CN110984976B (zh) | 2019-12-20 | 2019-12-20 | 多层合采气井产量劈分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911325055.9A CN110984976B (zh) | 2019-12-20 | 2019-12-20 | 多层合采气井产量劈分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110984976A CN110984976A (zh) | 2020-04-10 |
CN110984976B true CN110984976B (zh) | 2023-02-28 |
Family
ID=70073528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911325055.9A Active CN110984976B (zh) | 2019-12-20 | 2019-12-20 | 多层合采气井产量劈分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110984976B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111411946B (zh) * | 2020-05-12 | 2021-11-16 | 中国石油大学(北京) | 一种致密气藏气井开采方式的确定方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102747991A (zh) * | 2011-04-22 | 2012-10-24 | 长江大学 | 一种确定合采井单层产量的方法 |
CN106127599A (zh) * | 2016-07-04 | 2016-11-16 | 中国石油大学(华东) | 用于致密气藏在集气站生产模式下的气井产量劈分方法 |
CN108729912A (zh) * | 2018-04-19 | 2018-11-02 | 中国石油化工股份有限公司 | 适用于油藏数值模拟的产量劈分方法 |
CN108804819A (zh) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | 一种低渗气藏动态储量评价方法 |
CN108894777A (zh) * | 2018-07-06 | 2018-11-27 | 西南石油大学 | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 |
CN109858075A (zh) * | 2018-12-13 | 2019-06-07 | 朱亚婷 | 一种基于动静耦合的油井分层产量劈分方法 |
-
2019
- 2019-12-20 CN CN201911325055.9A patent/CN110984976B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102747991A (zh) * | 2011-04-22 | 2012-10-24 | 长江大学 | 一种确定合采井单层产量的方法 |
CN106127599A (zh) * | 2016-07-04 | 2016-11-16 | 中国石油大学(华东) | 用于致密气藏在集气站生产模式下的气井产量劈分方法 |
CN108729912A (zh) * | 2018-04-19 | 2018-11-02 | 中国石油化工股份有限公司 | 适用于油藏数值模拟的产量劈分方法 |
CN108804819A (zh) * | 2018-06-10 | 2018-11-13 | 西南石油大学 | 一种低渗气藏动态储量评价方法 |
CN108894777A (zh) * | 2018-07-06 | 2018-11-27 | 西南石油大学 | 一种分层压裂多层合采油气藏储层物性及裂缝特性参数的确定方法 |
CN109858075A (zh) * | 2018-12-13 | 2019-06-07 | 朱亚婷 | 一种基于动静耦合的油井分层产量劈分方法 |
Non-Patent Citations (3)
Title |
---|
协调点分析在多层合采井产量劈分中的应用;张连枝等;《石油化工应用》;20200131;第39卷(第1期);第47-50页 * |
基于气水相渗的合采气井产量劈分方法;姜宇玲;《大庆石油地质与开发》;20151031;第34卷(第5期);第73-76页 * |
苏里格气田苏东区块多层产量劈分新方法;林孟雄;《新疆地质》;20190930;第37卷(第3期);第419-421页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110984976A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7580797B2 (en) | Subsurface layer and reservoir parameter measurements | |
CN108804819A (zh) | 一种低渗气藏动态储量评价方法 | |
US20100241410A1 (en) | Relative and Absolute Error Models for Subterranean Wells | |
CN103912269B (zh) | 页岩气储层地层破裂压力梯度录井确定方法 | |
CN104504604B (zh) | 一种定性气井井筒积液的方法 | |
CN101906966A (zh) | 储层产能的预测方法及装置 | |
CN106194154B (zh) | 一种非常规油气藏中长期产能预测方法 | |
CN103590824A (zh) | 经过多段压裂改造后的致密气藏水平井的产能计算方法 | |
US10215002B2 (en) | Determining wellbore leak crossflow rate between formations in an injection well | |
CN106869917A (zh) | 基于井底和井口流入动态曲线共同约束的产能评价方法 | |
CN109488282A (zh) | 可动储量物性下限的确定方法 | |
CN109815543A (zh) | 计算气田动态储量的方法 | |
EP3707345B1 (en) | Determining wellbore leak crossflow rate between formations in an injection well | |
CN108729912B (zh) | 适用于油藏数值模拟的产量劈分方法 | |
CN110984976B (zh) | 多层合采气井产量劈分方法 | |
CN109958431A (zh) | 一种基于流入动态曲线计算油井日产液量的方法 | |
US20110276271A1 (en) | Method to determine current condensate saturation in a near-wellbore zone in a gas-condensate formation | |
Kumar et al. | Machine learning applications for a qualitative evaluation of the fracture network in the Wolfcamp shale using tracer and completion data | |
CN107725035B (zh) | 一种用于非均匀产液水平井判别来水方向的压力监测方法 | |
CN112836438B (zh) | 一种非均质边水气藏气井见水时间预测方法及装置 | |
CN109594965A (zh) | 一种不同生产压差下二氧化碳驱技术极限井距优化方法 | |
CN116050626B (zh) | 一种多级压裂水平井生产分级优化与效果预测方法 | |
CN113109162A (zh) | 一种基于热流固耦合的岩石起裂压力计算方法 | |
CN110630243B (zh) | 基于压裂井生产数据确定压裂裂缝参数的方法 | |
CN103590813B (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 |