CN109403961B - 多层合采油藏无因次采液指数预测方法 - Google Patents
多层合采油藏无因次采液指数预测方法 Download PDFInfo
- Publication number
- CN109403961B CN109403961B CN201811153031.5A CN201811153031A CN109403961B CN 109403961 B CN109403961 B CN 109403961B CN 201811153031 A CN201811153031 A CN 201811153031A CN 109403961 B CN109403961 B CN 109403961B
- Authority
- CN
- China
- Prior art keywords
- dimensionless
- index
- layer
- water
- production
- 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 116
- 239000007788 liquid Substances 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 33
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 100
- 239000012530 fluid Substances 0.000 claims abstract description 88
- 230000035699 permeability Effects 0.000 claims abstract description 47
- 238000000605 extraction Methods 0.000 claims abstract description 28
- 238000011161 development Methods 0.000 claims abstract description 16
- 238000005065 mining Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000013461 design Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000000704 physical effect Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 abstract description 3
- 238000011084 recovery Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000630 rising effect Effects 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
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)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种多层合采油藏无因次采液指数预测方法，该多层合采油藏无因次采液指数预测方法包括：步骤1，利用相渗数据计算各个分层的无因次采液指数；步骤2，获得分层的产液剖面及含水剖面；步骤3，插值得到分层不同含水对应的无因次采液指数；步骤4，液量加权平均求得合采无因次采液指数；步骤5，拟合外推得到全过程无因次采液指数。该多层合采油藏无因次采液指数预测方法充分反映特高含水开发阶段油藏的动、静态叠加造成高、低渗层段渗流能力的差异，更为准确的描述特高含水后期的储层供液能力，有效指导中高渗水驱油藏特高含水开发后期的开发生产实践。
Description
技术领域
本发明涉及油藏开发技术领域，特别是涉及到一种多层合采油藏无因次采液指数预测方法。
背景技术
无因次采液指数为某一含水率下的采液指数与含水率为零时的采液指数(即采油指数)之比,是评价不同含水率时采液能力的指标，与储集层类型和油藏流体性质有关。准确认识油田无因次采液指数变化规律,为油田产能预测、提液时机优选以及油田的开发方案设计、开发调整及开发动态分析提供重要依据。
目前通用计算算法：石油大学出版社于2002年出版的《油气藏工程方法与应用》中，给出了由实验室得出的相对渗透率资料，利用多条相渗曲线直接平均法得到一条相渗曲线，绘出无因次采液指数曲线。
俞启泰，罗洪，陈素珍在《断块油气田》杂志1998年06期发表文章——河流相储层油藏无因次采液指数计算的典型通用公式。该文章根据建立的河流相储层渗透率平面非均质抽象模型以及河流相储层的典型参数，考虑油水粘度比、油水井数比、井距不均匀系数、井网密度、平面渗透率对数正态分布变异系数、概率50％处渗透率、油藏大小、断层边界、含水率共9个因素，进行27个方案的数值模拟。经处理数值模拟结果，得到无因次采液指数与含水率关系的典型通用公式。
上述方法在中高含水阶段较好地指导了油田开发生产实践。然而，随着含水的升高，理论计算或经验公式的适用性逐渐降低。分析认为，中高含水开发阶段，多层合采油藏中的高、低渗层段渗流能力的差距较小，渗透率级差起着较强的控制作用。随着含水上升，多层合采油藏中的高、低渗层段渗流能力的差距逐渐加大。在特高含水后期，剩余油饱和度的微小变化，会引起水相渗透率的急剧变化，相对渗透率变化引起的动态非均质和绝对渗透率引起的静态非均质迭加后，多层合采油藏高、低渗层段的渗流能力差距急剧增加，增加到上百倍。这种渗流能力的显著差异在特高含水阶段是不可忽略的，而传统的无因次采液指数计算方法均无法体现该阶段性特征。为此我们发明了一种新的多层合采油藏无因次采液指数预测方法，解决了以上技术问题。
发明内容
本发明的目的是提供一种更接近多层合采油藏实际的特高含水期无因次采液指数预测，在特高含水后期能准确预测多层合采油藏无因次采液指数的方法。
本发明的目的可通过如下技术措施来实现：
多层合采油藏无因次采液指数预测方法，该多层合采油藏无因次采液指数预测方法包括：步骤1，利用相渗数据计算各个分层的无因次采液指数；步骤2，获得分层的产液剖面及含水剖面；步骤3，插值得到分层不同含水对应的无因次采液指数；步骤4，液量加权平均求得合采无因次采液指数；步骤5，拟合外推得到全过程无因次采液指数。
本发明的目的还可通过如下技术措施来实现：
在步骤1，根据计算油藏的参数，选取多层合采油藏各层的典型油水两相相对渗透率数据，计算得到各层的含水-无因次采液指数数据。
在步骤1，根据实验室中做出的相对渗透率曲线资料，绘制出无因次采油指数αo、无因次采液指数αl随含水fw变化曲线；
计算无因次采油指数αo的公式如下：
式中：
Kro(sw)：不同含水饱和度sw下的油相相对渗透率；
Kromax：束缚水Swi下的油相相对渗透率；
K：fw＝0时的油层绝对渗透率；
Kw：含水为fw时的油层绝对渗透率。
如果不考虑注水开发过程中绝对渗透率的变化，令K＝Kw，则上式变为：
此式即为无因次采油指数的计算公式；
无因次采液指数αl的计算公式如下：
在步骤2，通过矿场产液剖面测试或根据多层合采油藏物性、流体特征开展多管并联物理模拟试验，确定多层合采油藏合采时不同含水时刻各小层的产液量及含水点分布数据。
在步骤3，根据步骤2得到的各小层含水剖面的点分布数据，在相邻含水点之间各个分层的无因次采液指数近似线性变化，依据相应小层在步骤1中得到的含水-无因次采液指数变化规律，采用线性差值的方法到各小层不同含水对应的无因次采液指数。
在步骤3，某一小层在含水f1对应的求得相渗数据计算无因次采液指数为αl1，在含水f3对应的求得相渗数据计算无因次采液指数为αl3，要求取在含水f1与含水f3之间的含水f2对应的相渗数据计算无因次采液指数为αl2，应用如下公式计算：
在步骤4，根据不同合采时刻各小层产液量比例，对各小层相应含水对应的无因次采液指数进行加权平均，计算得到多层合采油藏的无因次采液指数。
在步骤4，对分层无因次采液指数进行液量加权平均的计算公式为：
式中：Jos权为多层合采油藏的无因次采液指数，为小数；fw为多层合采油藏综合含水，小数；fw1、fw2分别对分层的综合含水，小数；Jos1、Jos2分别对分层的无因次采液指数，小数，QB为两个层的产液倍数，为层位二/层位一，小数。
在步骤5，将合采的含水-无因次采液指数的计算结果进行多项式拟合，外推至含水100％，建立驱替全过程的含水-无因次采液指数关系。
本发明中的多层合采油藏无因次采液指数预测方法，涉及到中高渗透水驱开发油藏在特高含水后期无因次采液指数的预测方法，充分反映特高含水开发阶段油藏的动、静态叠加造成高、低渗层段渗流能力的差异，更为准确的描述特高含水后期的储层供液能力，有效指导中高渗水驱油藏特高含水开发后期的开发生产实践。
附图说明
图1为本发明的多层合采油藏无因次采液指数预测方法的一具体实施例的流程图；
图2为本发明的一具体实施例中A油藏特高含水期无因次采液指数预测结果示意图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特举出较佳实施例，并配合附图所示，作详细说明如下。
如图1所示，图1为本发明的多层合采油藏无因次采液指数预测方法的流程图。
步骤101，计算分层的无因次采液指数。选取多层合采油藏各层的典型油水两相相对渗透率数据，计算得到各层的含水-无因次采液指数数据；
根据实验室中做出的相对渗透率曲线资料，绘制出无因次采油指数αo、无因次采液指数αl随含水fw变化曲线。
计算无因次采油指数(αo)的公式如下：
式中：
Kro(sw)：不同含水饱和度sw下的油相相对渗透率；
Kromax：束缚水Swi下的油相相对渗透率；
K：fw＝0时的油层绝对渗透率；
Kw：含水为fw时的油层绝对渗透率。
如果不考虑注水开发过程中绝对渗透率的变化，令K＝Kw，则上式变为：
此式即为无因次采油指数的计算公式。
无因次采液指数αl的计算公式如下：
步骤102，获得分层的产液剖面及含水剖面。通过矿场产液剖面测试或根据多层合采油藏物性、流体特征开展多管并联物理模拟试验，确定多层合采油藏合采时不同含水时刻各小层的产液量及含水分布数据；
步骤103，插值得到分层不同含水对应的无因次采液指数。根据步骤102得到的各小层含水剖面，依据相应小层在步骤101中得到的含水-无因次采液指数变化规律，插值得到各小层不同含水对应的无因次采液指数；
例如某一小层在含水f1对应的求得相渗数据计算无因次采液指数为αl1，在含水f3对应的求得相渗数据计算无因次采液指数为αl3，要求取在含水f1与含水f3之间的含水f2对应的相渗数据计算无因次采液指数为αl2，应用如下公式计算：
步骤104，液量加权平均求得合采无因次采液指数。根据不同合采时刻各小层产液量比例，对各小层相应含水对应的无因次采液指数进行加权平均，计算得到多层合采油藏的无因次采液指数；
对分层无因次采液指数进行液量加权平均的计算公式为：
式中：Jos权为多层合采油藏的无因次采液指数，为小数；fw为多层合采油藏综合含水，小数；fw1、fw2分别对分层的综合含水，小数；Jos1、Jos2分别对分层的无因次采液指数，小数，QB为两个层的产液倍数，为层位二/层位一，小数。
步骤105，拟合外推得到全过程无因次采液指数。将合采的含水-无因次采液指数的计算结果进行多项式拟合，外推至含水100％，建立驱替全过程的含水-无因次采液指数关系。
在应用本发明的一具体实施例中，包括了以下步骤：
在步骤1中，A油藏为两层合采，层位一渗透率800毫达西，层位二渗透率2100毫达西，选取具备代表性的室内测得的相对渗透率曲线得到各层的无因次采液指数与含水的变化关系(见表1)。
表1.A油藏不同层位的无因次采液指数数据
在步骤2中，根据层位一、层位二的物性特征(层位一渗透率800毫达西，层位二渗透率2100毫达西)，开展双管并联室内物理模拟实验，确定合采不同含水条件下，各层的产液量和含水分布状况(见表2)。
表2.A油藏合采及各层产液剖面、含水剖面状况表
在步骤3中，根据两个层合采时各个层的含水(见表2)对应与步骤1提供的无因次采液指数数据，插值得到表3各个层对应的含水-无因次采液指数关系数据，具体数据见表3。
表3.A油藏各层含水-无因次采液指数数据表
在步骤4中，按如下公式对分层无因次采液指数进行液量加权平均，计算合采油藏特高含水期无因次采液指数，计算结果见表4。
式中：Jos权为多层合采油藏的无因次采液指数，为小数；fw为多层合采油藏综合含水，小数；fw1、fw2分别对分层的综合含水，小数；Jos1、Jos2分别对分层的无因次采液指数，小数，QB为两个层的产液倍数，为层位二/层位一，小数。
表4.A油藏合采情况下含水-无因次采液指数数据表
在步骤5中，按步骤4中计算得到的结果与合采综合含水进行拟合，外推至含水100％，得到全部含水过程的无因次采液指数(见图2)。
同过对比A油藏无因次采液指数的实际值与计算值(见表5)，计算误差小于10％，满足矿场上对产液量决策的计算需求。
表5.A油藏无因次采液指数实际值与计算值对比情况表
Claims (4)
1.多层合采油藏无因次采液指数预测方法，其特征在于，该多层合采油藏无因次采液指数预测方法包括：
步骤1，利用相渗数据计算各个分层的无因次采液指数；
步骤2，获得分层的产液剖面及含水剖面；
步骤3，插值得到分层不同含水对应的无因次采液指数；
步骤4，液量加权平均求得合采无因次采液指数；
步骤5，拟合外推得到全过程无因次采液指数；
在步骤3，根据步骤2得到的各小层含水剖面的点分布数据，在相邻含水点之间各个分层的无因次采液指数近似线性变化，依据相应小层在步骤1中得到的含水-无因次采液指数变化规律，采用线性差值的方法到各小层不同含水对应的无因次采液指数；某一小层在含水f1对应的求得相渗数据计算无因次采液指数为αl1，在含水f3对应的求得相渗数据计算无因次采液指数为αl3，要求取在含水f1与含水f3之间的含水f2对应的相渗数据计算无因次采液指数为αl2，应用如下公式计算：
在步骤4，根据不同合采时刻各小层产液量比例，对各小层相应含水对应的无因次采液指数进行加权平均，计算得到多层合采油藏的无因次采液指数；对分层无因次采液指数进行液量加权平均的计算公式为：
式中：Jos权为多层合采油藏的无因次采液指数，为小数；fw为多层合采油藏综合含水，小数；fw1、fw2分别对分层的综合含水，小数；Jos1、Jos2分别对分层的无因次采液指数，小数，QB为两个层的产液倍数，为层位二/层位一，小数；
在步骤5，将合采的含水-无因次采液指数的计算结果进行多项式拟合，外推至含水100％，建立驱替全过程的含水-无因次采液指数关系。
2.根据权利要求1所述的多层合采油藏无因次采液指数预测方法，其特征在于，在步骤1，根据计算油藏的参数，选取多层合采油藏各层的典型油水两相相对渗透率数据，计算得到各层的含水-无因次采液指数数据。
3.根据权利要求2所述的多层合采油藏无因次采液指数预测方法，其特征在于，在步骤1，根据实验室中做出的相对渗透率曲线资料，绘制出无因次采油指数αo、无因次采液指数αl随含水fw变化曲线；
计算无因次采油指数αo的公式如下：
式中：
Kro(sw)：不同含水饱和度sw下的油相相对渗透率；
Kromax：束缚水Swi下的油相相对渗透率；
K：fw＝0时的油层绝对渗透率；
Kw：含水为fw时的油层绝对渗透率。
如果不考虑注水开发过程中绝对渗透率的变化，令K＝Kw，则上式变为：
此式即为无因次采油指数的计算公式；
无因次采液指数αl的计算公式如下：
4.根据权利要求1所述的多层合采油藏无因次采液指数预测方法，其特征在于，在步骤2，通过矿场产液剖面测试或根据多层合采油藏物性、流体特征开展多管并联物理模拟试验，确定多层合采油藏合采时不同含水时刻各小层的产液量及含水点分布数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811153031.5A CN109403961B (zh) | 2018-09-29 | 2018-09-29 | 多层合采油藏无因次采液指数预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811153031.5A CN109403961B (zh) | 2018-09-29 | 2018-09-29 | 多层合采油藏无因次采液指数预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109403961A CN109403961A (zh) | 2019-03-01 |
| CN109403961B true CN109403961B (zh) | 2021-09-28 |
Family
ID=65466618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811153031.5A Active CN109403961B (zh) | 2018-09-29 | 2018-09-29 | 多层合采油藏无因次采液指数预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109403961B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110863814B (zh) * | 2019-11-04 | 2022-10-04 | 中国石油天然气股份有限公司 | 巨厚型生屑灰岩油藏单井分段比采液指数确定方法及装置 |
| CN111749664B (zh) * | 2020-07-01 | 2022-06-07 | 中国海洋石油集团有限公司 | 一种油田全程化学驱组合接替时机优化方法和系统 |
| US11499424B2 (en) | 2021-02-18 | 2022-11-15 | Saudi Arabian Oil Company | Systems and methods to determine the productivity index of individual laterals under commingled flow |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101942994A (zh) * | 2010-09-16 | 2011-01-12 | 中国石油天然气股份有限公司 | 水淹层产水率定量预测方法及其系统 |
| CN105625992A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 针对屋脊断块油藏的合理提液方法 |
| CN106014359A (zh) * | 2016-06-08 | 2016-10-12 | 西南石油大学 | 一种海上油藏早期注聚最早转注时机判断方法 |
| CN107366529A (zh) * | 2017-06-28 | 2017-11-21 | 中国石油化工股份有限公司 | 非均质水驱油藏无因次采液指数曲线表征方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7893692B2 (en) * | 2008-11-03 | 2011-02-22 | Schlumberger Technology Corporation | Method for estimating the formation productivity from nuclear magnetic resonance measurements |
| CN105626036B (zh) * | 2014-11-07 | 2019-05-31 | 中国石油化工股份有限公司 | 一种确定油藏合理产液量油藏工程计算方法 |
| US20160154907A1 (en) * | 2014-12-01 | 2016-06-02 | Schlumberger Technology Corporation | Integrated network asset modeling |
-
2018
- 2018-09-29 CN CN201811153031.5A patent/CN109403961B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101942994A (zh) * | 2010-09-16 | 2011-01-12 | 中国石油天然气股份有限公司 | 水淹层产水率定量预测方法及其系统 |
| CN105625992A (zh) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | 针对屋脊断块油藏的合理提液方法 |
| CN106014359A (zh) * | 2016-06-08 | 2016-10-12 | 西南石油大学 | 一种海上油藏早期注聚最早转注时机判断方法 |
| CN107366529A (zh) * | 2017-06-28 | 2017-11-21 | 中国石油化工股份有限公司 | 非均质水驱油藏无因次采液指数曲线表征方法 |
Non-Patent Citations (2)
| Title |
|---|
| Index system evaluating water flooding development effect of oilfield at ultra-high water cut stage;Hua Wen等;《Journal of Petroleum Exploration and Production Technology 》;20160509;第111-123页 * |
| 多层合采复杂水淹层产能预测新方法;任晓宁等;《重庆科技学院学报（自然科学版）》;20121231;第14卷(第6期);第79-82页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109403961A (zh) | 2019-03-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106777651B (zh) | 注采平衡原理的油水井产量劈分方法 | |
| CN106894814B (zh) | 复杂断块油藏高含水后期剩余油二次富集的快速识别方法 | |
| CN109403961B (zh) | 多层合采油藏无因次采液指数预测方法 | |
| CN107038268B (zh) | 一种确定非均质储层五点井网水驱波及系数的方法 | |
| CN110259444B (zh) | 基于流场诊断的水驱油藏渗流场可视化表征与评价方法 | |
| CN103939066A (zh) | 一种一注多采井组定注水量确定油井产液量的方法 | |
| CN105626010A (zh) | 一种分段注水井中注水层段合理划分方法 | |
| CN111810101A (zh) | 一种水驱油藏动态分析方法及装置 | |
| CN113297740A (zh) | 一种水平井凝胶调剖参数优化方法 | |
| CN110782187B (zh) | 基于含水率的枯竭油气藏型储气库动态评价方法及系统 | |
| CN106481315B (zh) | 陆上砂岩油藏单井可采储量快速确定模型及建立方法 | |
| CN114462262A (zh) | 一种基于时间和空间双重维度的历史拟合预测方法 | |
| CN112746836B (zh) | 基于层间干扰的油井各层产量计算方法 | |
| CN110851982B (zh) | 一种中低渗油藏减氧空气驱可行性分析方法 | |
| CN112324429B (zh) | 确定有水气藏的合理生产压差的方法 | |
| CN115705452A (zh) | 整装砂岩油藏开发中后期的新型采收率预测方法 | |
| CN111287739B (zh) | 一种基于地层原油粘度的剩余油分布预测方法 | |
| CN105089616A (zh) | 改进的技术可采储量标定方法 | |
| CN110410044B (zh) | 气驱co2、n2开发方式下区块产油量计算方法 | |
| CN111428375A (zh) | 一种高含水油田层系重组划分新方法 | |
| CN112360410A (zh) | 一种提高低渗透油藏水驱驱油效率的方法 | |
| CN114004078B (zh) | 一种表征薄互层油藏多层合采渗流拟相渗计算方法 | |
| Bui et al. | Material balance analysis in complex mature reservoirs–experience in Samarang field, Malaysia | |
| CN113762559A (zh) | 多层合采气井智能产量劈分方法及系统 | |
| Alem et al. | Predicting the performance of tight gas reservoirs |
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 |