CN111222261B - 动静参数结合的产量劈分方法 - Google Patents
动静参数结合的产量劈分方法 Download PDFInfo
- Publication number
- CN111222261B CN111222261B CN202010322869.3A CN202010322869A CN111222261B CN 111222261 B CN111222261 B CN 111222261B CN 202010322869 A CN202010322869 A CN 202010322869A CN 111222261 B CN111222261 B CN 111222261B
- Authority
- CN
- China
- Prior art keywords
- layer
- gas production
- well
- gas
- coefficient
- 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
- 238000000034 method Methods 0.000 title claims abstract description 44
- 230000003068 static effect Effects 0.000 title claims abstract description 29
- 230000000694 effects Effects 0.000 claims abstract description 28
- 238000012795 verification Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 22
- 230000015572 biosynthetic process Effects 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 9
- 230000035699 permeability Effects 0.000 claims description 6
- 239000003245 coal Substances 0.000 claims description 3
- 210000002615 epidermis Anatomy 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 210000003491 skin Anatomy 0.000 claims 1
- 238000011156 evaluation Methods 0.000 abstract description 2
- 238000011161 development Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mining & Mineral Resources (AREA)
- Human Resources & Organizations (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- Development Economics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及动静参数结合的产量劈分方法,属于油气藏产能评价领域;它解决现今产量劈分误差较大,劈分效果不理想等问题;其技术方案是:多次测量产气剖面立标准,首先调整各层产能系数及单井控制储量初值使迭代计算的各层产量贡献率与实测多次产气剖面数据进行拟合,拟合成功后计算产能系数以及单井控制储量初值与实际值相差倍数以作标准;第二,一次测量产气剖面以验证,利用倍数,结合地层参数计算实际各层单井控制储量及产能系数,迭代计算出各层产量贡献率与该次产气剖面数据进行验证;第三,没有产剖的井作预测,结合所得的倍数,迭代预测各层产量贡献率。本发明结合产气剖面数据和地层参数实现产量劈分,拟合效果较好,可推广性强。
Description
技术领域
本发明涉及动静参数结合的产量劈分方法,属于油气藏产能评价领域。
背景技术
在油气田开发过程中,储层动用情况的研究是调整开发方案、综合挖潜和提高油藏采收率的基础和根据。但如何进一步提高储层动用状况的研究水平并能有效应用到具体的油田开发过程中始终是一个需要不断研究、不断完善的重大问题。找到一种准确、现实可行的产量劈分方法是油气藏储量动用状况研究的重点。
油气藏储量的动用状况大体上受到油气藏地质条件和开发方案的影响。客观存在的油气藏地质条件影响着油气水运动规律,人为的开发方案影响着储层平面和空间的动用程度大小。因此储层的产量劈分要综合考虑各方面的因素。
现场常用测量产气剖面的方法确定各各层产量贡献率,但是产气剖面测量成本较高,并且只能测出某时间点的瞬态产量贡献率,无法反映整个生产阶段的产量贡献率变化。
常规的产量劈分方法主要有三大类:第一类为参数法,主要包括有效厚度法、地层系数法等,该类方法考虑因素单一,仅考虑到地层参数,误差较大;第二类为数值模拟法,该类方法计算模型复杂,需要经过建模、历史拟合等过程,多解性强,计算速度慢;第三类为突变法,尽管考虑因素较多,但其本质仍是考虑静态参数实现产量劈分,精度有所提高,但仍存在较大误差。
发明内容
本发明目的是:为了解决现今产量劈分误差较大,劈分效果不理想等问题,本发明结合动静参数实现产量劈分,可以结合动态产气剖面数据和静态地层参数实现产量劈分,拟合效果较好,可推广性强。
为实现上述目的,本发明提供了动静参数结合的产量劈分方法,该方法包括下列步骤:
第一,多次测量产气剖面确立标准,首先利用多次测量产气剖面的井数据,结合气井产能方程以及定容气藏物质平衡方程,进行迭代,计算出各层产量;以地层参数所计算得到的产能系数以及单井控制储量为初值,不断调整各层产能系数及单井控制储量使迭代计算的各层产量贡献率与多次产气剖面数据进行拟合,拟合成功后,得到实际的各层产能系数及单井控制储量;再分别计算产能系数以及单井控制储量的初值与实际值的相差倍数,将该倍数作为周边井区计算实际产能系数以及单井控制储量的标准;
第二,一次测量产气剖面以验证,导入一次测量产气剖面的井数据,根据地层参数计算产能系数以及单井控制储量的初值,除以第一步中计算出的倍数,即可得到实际的各层产能系数以及单井控制储量;再进行迭代计算,计算出各个时刻下各层产量贡献率,利用一次测量的产气剖面进行验证;
第三,没有测量产气剖面的井作预测,根据地层参数计算产能系数以及单井控制储量的初值,除以第一步中计算出的倍数,即可得到实际的各层产能系数以及单井控制储量,将实际的各层产能系数以及单井控制储量导入进行迭代,拟合单井日产气量,预测各层产量贡献率。
上述动静参数结合的产量劈分方法中,所述迭代思路为,
第一步,利用总产气量Q sc 及套压P t 反算井底流压P wf ;
再利用第四步所得各层地层压力(),从第二步开始继续进行迭代;利用一次迭代所得的各层产量()计算各层产量贡献率,第n层产量贡献率计算公式为,生产数据中的一天为一个迭代周期,不断迭代将得到整个生产阶段的气井各层产量贡献率,迭代完成后最终的各层可采储量即为各层单井控制储量。
上述动静参数结合的产量劈分方法中,计算各层产能系数公式为,其中K为各层渗透率,h为各层厚度,T为各层温度,为各层气体粘度,为各层气体偏差因子,S为各层表皮系数;计算各层单井控制储量初值的步骤为,首先利用容积法计算该井周边半井距范围内的原始地质储量,再利用储量动静比与原始地质储量相乘即可得到各层单井控制储量的初值。
上述动静参数结合的产量劈分方法中,针对有产气剖面数据的井,所述的拟合效果包括两方面,一是各层产量贡献率与实测产气剖面数据的拟合效果,二是迭代计算得到的单井日产气量与实测日产气量的拟合效果;针对没有产气剖面数据的井,所述的拟合效果是指迭代计算得到的单井日产气量与实测日产气量的拟合效果,拟合效果较好是指拟合成功率80%以上。
上述动静参数结合的产量劈分方法中,所述的动静参数中动参数是指产气剖面数据,静参数是指地层参数如渗透率、地层厚度。
与现有技术相比,本发明具有以下有益效果:(1)结合动态产气剖面数据和静态地层参数实现产量劈分,拟合效果较好;(2)编程实现劈分,劈分过程省时省力;(3)可推广性强。
附图说明
在附图中:
图1是本方法技术路线图。
图2是某标准井日产气量拟合图。
图3是某一次测量产气剖面井验证图。
图4是某井各层产量贡献率预测图。
具体实施方式
下面结合实施方式和附图对本发明做进一步说明。
本发明提供了动静参数结合的产量劈分方法,图1为本方法的技术路线图,该方法包括下列步骤:
第一,多次测量产气剖面确立标准,首先利用多次测量产气剖面的井数据,结合气井产能方程以及定容气藏物质平衡方程,进行迭代,计算出各层产量;以地层参数所计算得到的产能系数以及单井控制储量为初值,不断调整各层产能系数及单井控制储量使迭代计算的各层产量贡献率与多次产气剖面数据进行拟合,拟合成功后,得到实际的各层产能系数及单井控制储量;再分别计算产能系数以及单井控制储量的初值与实际值的相差倍数,将该倍数作为周边井区计算实际产能系数以及单井控制储量的标准,如图2所示,该图为某标准井日产气量拟合图,整体拟合效果较好。
第二,一次测量产气剖面以验证,导入一次测量产气剖面的井数据,根据地层参数计算产能系数以及单井控制储量的初值,除以第一步中计算出的倍数,即可得到实际的各层产能系数以及单井控制储量;再进行迭代计算,计算出各个时刻下各层产量贡献率,利用一次测量的产气剖面进行验证;如果迭代计算得到的产量贡献率数据与该次产气剖面数据拟合效果较好,并且迭代计算后的单井日产气量与实际日产气量拟合效果较好,则该方法在该区块得到了验证,方法可靠;如图3所示,该图为标准井周边一口一次测量产气剖面的井,迭代计算得到的产量贡献率数据与一次测量的产气剖面数据拟合效果较好。
第三,没有测量产气剖面的井作预测,根据地层参数计算产能系数以及单井控制储量的初值,除以第一步中计算出的倍数,即可得到实际的各层产能系数以及单井控制储量,将实际的各层产能系数以及单井控制储量导入进行迭代,拟合单井日产气量,预测各层产量贡献率;若迭代计算得到的单井日产气量与实测日产气量的拟合效果,则认为所计算得到的各层产量贡献率及各层控制储量可靠,如图4所示,该井未进行产气剖面测试,应用该方法得到各层产量贡献率预测图。
进一步的,所述迭代思路为,
第一步,利用总产气量Q sc 及套压P t 反算井底流压P wf ;
再利用第四步所得各层地层压力(),从第二步开始继续进行迭代;利用一次迭代所得的各层产量()计算各层产量贡献率,第n层产量贡献率计算公式为,生产数据中的一天为一个迭代周期,不断迭代将得到整个生产阶段的气井各层产量贡献率,迭代完成后最终的各层可采储量即为各层单井控制储量。
进一步的,所述计算各层产能系数公式为,其中K为各层渗透率,h为各层厚度,T为各层温度,为各层气体粘度,为各层气体偏差因子,S为各层表皮系数;计算各层单井控制储量初值的步骤为,首先利用容积法计算该井周边半井距范围内的原始地质储量,再利用储量动静比与原始地质储量相乘即可得到各层单井控制储量的初值。
进一步的,针对有产气剖面数据的井,所述的拟合效果包括两方面,一是各层产量贡献率与实测产气剖面数据的拟合效果,二是迭代计算得到的单井日产气量与实测日产气量的拟合效果;针对没有产气剖面数据的井,所述的拟合效果是指迭代计算得到的单井日产气量与实测日产气量的拟合效果,拟合效果较好是指拟合成功率80%以上。
进一步的,所述的动静参数中动参数是指产气剖面数据,静参数是指地层参数如渗透率、地层厚度。
与现有技术相比,本发明具有以下有益效果:(1)结合动态产气剖面数据和静态地层参数实现产量劈分,拟合效果较好;(2)编程实现劈分,劈分过程省时省力;(3)可推广性强。
最后所应说明的是:以上实施例仅用以说明而非限制本发明的技术方案,尽管参照上述实施例对本发明进行了详细说明,本领域的普通技术人员应该理解:依然可以对本发明进行修改或者等同替换,而不脱离本发明的精神和范围的任何修改或局部替换,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.动静参数结合的产量劈分方法,其特征在于,该方法包括下列步骤:
第一,多次测量产气剖面确立标准,首先利用多次测量产气剖面的井数据,结合气井产能方程以及定容气藏物质平衡方程,进行迭代,计算出各层产量;以地层参数所计算得到的产能系数以及单井控制储量为初值,不断调整各层产能系数及单井控制储量使迭代计算的各层产量贡献率与多次产气剖面数据进行拟合,拟合成功后,得到实际的各层产能系数及单井控制储量;再分别计算产能系数以及单井控制储量的初值与实际值的相差倍数,将该倍数作为周边井区计算实际产能系数以及单井控制储量的标准;
第二,一次测量产气剖面以验证,导入一次测量产气剖面的井数据,根据地层参数计算产能系数以及单井控制储量的初值,除以第一步中计算出的倍数,即可得到实际的各层产能系数以及单井控制储量;再进行迭代计算,计算出各个时刻下各层产量贡献率,利用一次测量的产气剖面进行验证;
第三,没有测量产气剖面的井作预测,根据地层参数计算产能系数以及单井控制储量的初值,除以第一步中计算出的倍数,即可得到实际的各层产能系数以及单井控制储量,将实际的各层产能系数以及单井控制储量导入进行迭代,拟合单井日产气量,预测各层产量贡献率。
4.根据权利要求1所述的动静参数结合的产量劈分方法,其特征在于:所述迭代思路为,
第一步,利用总产气量Q sc 及套压P t 反算井底流压P wf ;
6.根据权利要求1所述的动静参数结合的产量劈分方法,其特征在于:针对有产气剖面数据的井,所述的拟合效果包括两方面,一是各层产量贡献率与实测产气剖面数据的拟合效果,二是迭代计算得到的单井日产气量与实测日产气量的拟合效果;针对没有产气剖面数据的井,所述的拟合效果是指迭代计算得到的单井日产气量与实测日产气量的拟合效果。
7.根据权利要求1所述的动静参数结合的产量劈分方法,其特征在于:所述的动静参数中动参数是指产气剖面数据,静参数是指渗透率、地层厚度、原始地层压力、地层温度、气体粘度、气体偏差因子以及表皮系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010322869.3A CN111222261B (zh) | 2020-04-22 | 2020-04-22 | 动静参数结合的产量劈分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010322869.3A CN111222261B (zh) | 2020-04-22 | 2020-04-22 | 动静参数结合的产量劈分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111222261A CN111222261A (zh) | 2020-06-02 |
CN111222261B true CN111222261B (zh) | 2020-07-14 |
Family
ID=70827360
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010322869.3A Active CN111222261B (zh) | 2020-04-22 | 2020-04-22 | 动静参数结合的产量劈分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111222261B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111625757B (zh) * | 2020-06-09 | 2024-03-01 | 承德石油高等专科学校 | 一种数字化油田预测单井产量的算法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106127599A (zh) * | 2016-07-04 | 2016-11-16 | 中国石油大学(华东) | 用于致密气藏在集气站生产模式下的气井产量劈分方法 |
CN110617062A (zh) * | 2019-10-14 | 2019-12-27 | 西南石油大学 | 一种储层剩余气分布特征的分析方法 |
-
2020
- 2020-04-22 CN CN202010322869.3A patent/CN111222261B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106127599A (zh) * | 2016-07-04 | 2016-11-16 | 中国石油大学(华东) | 用于致密气藏在集气站生产模式下的气井产量劈分方法 |
CN110617062A (zh) * | 2019-10-14 | 2019-12-27 | 西南石油大学 | 一种储层剩余气分布特征的分析方法 |
Non-Patent Citations (4)
Title |
---|
A calculation model of critical liquid-carrying velocity of gas wells considering the influence of droplet shapes;Pan Jie等;《ScienceDirect》;20181231;第337-343页 * |
Flow model of a multi-stage hydraulic fractured horizontal well based on tree-shaped fractal fracture networks;Xiao-Hua Tan等;《Journal of Petroleum Science and Engineering》;20181231;494-503页 * |
三重介质复合气藏水平井不稳定产量递减动态分析;黄雨等;《天然气地球科学》;20180831;第29卷(第8期);第1190-1197页 * |
多层合采气藏产量劈分数值模拟研究;史进等;《海洋石油》;20150630;第35卷(第2期);第56-60、81页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111222261A (zh) | 2020-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102041995B (zh) | 复杂油藏水淹状况监测系统 | |
CN111927411B (zh) | 一种智能化有水气藏水侵跟踪和预警方法 | |
US8473268B2 (en) | Method for comparing and back allocating production | |
CA2793825C (en) | Automated field development planning | |
CN110334431A (zh) | 一种低渗透致密气藏单井控制储量计算及剩余气分析方法 | |
CN108984886B (zh) | 一种反演多层油藏井间动态连通性的方法 | |
CN103912248A (zh) | 水驱油田预测含水率方法 | |
CN111706317B (zh) | 一种确定加密调整区低渗储层剩余油分布状况的方法 | |
CN114427432B (zh) | 一种气藏剩余气开发潜力确定方法 | |
CN111222261B (zh) | 动静参数结合的产量劈分方法 | |
CN113297740A (zh) | 一种水平井凝胶调剖参数优化方法 | |
CA3116482C (en) | Reservoir fluid property modeling using machine learning | |
CN111706318B (zh) | 一种确定低渗储层剩余油分布状况的方法 | |
CN111506865B (zh) | 一种页岩气井生产拟合与预测方法 | |
CN104948177A (zh) | 一种基于米产液指数统计特征规律的产能预测方法 | |
CN110486008A (zh) | 一种径向复合油藏的参数解释方法及系统 | |
CN110486009A (zh) | 一种无限大地层的参数自动反求方法及系统 | |
Ali et al. | A semi-analytical method for history matching and improving geological models of layered reservoirs: CGM analytical method | |
CN113496301B (zh) | 油气田资产评价方法及装置 | |
Chen et al. | Practical shale gas decline analysis in Changning shale gas play-a variable pressure drop study | |
Xie et al. | Research of Dynamic Reserves Estimation Methods for Complex Gas Wells in Low Permeability and Tight Gas Reservoir | |
CN116184493A (zh) | 一种面向水平井开发的井震联合高精度速度建场方法 | |
Jreou | Using Ann for Well Type Identifying and Increasing Production from Sadi Formation of Halfaya Oil Field South of Iraq–Case Study | |
CN115271182A (zh) | 海上油田水驱采收率预测方法 | |
Danilović et al. | GAS-LIFT WELLS OPTIMIZATION AT THE OIL FIELD “K” |
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 |