CN109918769B - 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 - Google Patents
利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 Download PDFInfo
- Publication number
- CN109918769B CN109918769B CN201910161366.XA CN201910161366A CN109918769B CN 109918769 B CN109918769 B CN 109918769B CN 201910161366 A CN201910161366 A CN 201910161366A CN 109918769 B CN109918769 B CN 109918769B
- Authority
- CN
- China
- Prior art keywords
- water
- oil
- formation
- time
- 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.)
- Expired - Fee Related
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 113
- 230000009545 invasion Effects 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000001052 transient effect Effects 0.000 title description 7
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 41
- 230000006835 compression Effects 0.000 claims abstract description 14
- 238000007906 compression Methods 0.000 claims abstract description 14
- 238000005065 mining Methods 0.000 claims abstract description 5
- 239000007788 liquid Substances 0.000 claims abstract description 4
- 238000004519 manufacturing process Methods 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 abstract description 11
- 230000008569 process Effects 0.000 abstract description 8
- 238000000605 extraction Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 4
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 239000003129 oil well Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000037237 body shape Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
Images
Landscapes
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，先依据油气田日常计量数据获得相关参数；然后利用公式求出采液速度；再利用公式，根据初始地层压力、开采任意时间的地层压力计算出水侵系数；最后利用公式，根据环形水层函数、无因次时间、油层外半径、孔隙度、综合压缩系数、水的粘度，计算出水侵量。本发明提供了一种建立在瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，该方法可以解决缝洞型油藏在非稳态水侵过程中定量水侵量的问题，计算方法繁琐，计算精度误高。
Description
技术领域
本发明涉及采矿技术领域，具体来说是涉及油藏开发时定量非稳态水侵水侵量的方法。
背景技术
在油藏开发过程中，影响开发效果的主要原因就是水体的侵入，油井一旦见水，产油量会不断降低，产水量不断增加导致油井水淹，严重影响油藏的开发效果及最终采收率。
缝洞型储层种类多，大致可分为溶洞型，裂缝型及二者复合型。而每种类型由于其储层规模不同，接触方式不同导致其非均质性很强，没有一定规律性。如溶洞型储集层的水侵过程多为稳态水侵和拟稳态水侵，而裂缝型及二者复合型则非稳态水侵居多。因此在油藏开采过程中，压力供给方式较为复杂，水侵量难于系统计算。
水侵量分析是油藏动态分析、注水采油方案确定及控水措施实施的前期基础工作。非稳态水侵阶段计算水侵量时要考虑时间的影响，油水界面的压力是不断变化的，但和时间没有一个确定的连续函数，压力与每一时刻水侵量的关系呈现为散乱的点。
目前，常用的非稳态水侵量计算模型主要van Everdingen-Hurst非稳态模型，该模型假设油藏的厚度、渗透率、孔隙度、岩石和水的压缩性不变，且系统的驱动能量来源于水的膨胀作用和岩石的压缩性。需要对水体形态和大小进行一定的猜测，然后进行试算，计算方法繁琐，存在一定的误差。
发明内容
本发明要解决的技术问题在于，针对现有技术中常用的非稳态水侵量计算模型-van Everdingen-Hurst非稳态模型，需要对水体形态和大小进行一定的猜测，然后进行试算，计算方法繁琐，存在一定的误差的技术缺陷，提供了一种利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法。
本发明解决其技术问题所采用利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法包含如下步骤：
S1、获得油藏地质储量N、地层压力、地层压力为pi时的地层油体积系数Boi、累积产油量NP、地层压力为p时的地层油体积系数Bo、地层束缚水压缩系数Cw、地层束缚水的饱和度Swi、地层原始含油饱和度Soi、累积产水量Wp、水的体积系数Bw、环形水层函数QD、无因次时间tD、油层外半径r0、孔隙度
综合压缩系数ct以及水的粘度μw；其中，pi为初始地层压力，p为开采任意时间地层压力；
S2、利用下述公式求出采液速度dWe/dt；
其中，We为水侵量，t为时间，压差△P＝pi-p；
S3、利用下述公式，根据初始地层压力、开采任意时间的地层压力计算出水侵系数k′：
S4、利用下述公式，根据环形水层函数、无因次时间、油层外半径、孔隙度、综合压缩系数、水的粘度，计算出水侵量：
其中：QD为环形水层函数；tD为无因次时间；r0为油层外半径；
为孔隙度；ct为综合压缩系数；μw为水的粘度，pj是表示时间tj时的地层压力，压力增量
pi＝p-1，时间j与时间j+1之间的时间段为无限小时间段，ΔtD表示在j时间步长下的无因次时间,n表示无限小时间段的总段数。
进一步地，在本发明的利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，其特征在于，获得油藏地质储量N、地层压力、地层压力为pi时的地层油体积系数Boi、累积产油量NP、地层压力为p时的地层油体积系数Bo、地层束缚水压缩系数Cw、地层束缚水的饱和度Swi、地层原始含油饱和度Soi、累积产水量Wp、水的体积系数Bw、油层外半径r0、孔隙度
综合压缩系数ct以及水的粘度μw依据油气田日常计量数据获得。
实施本发明的利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，具有以下有益效果：本发明提供了一种建立在瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，该方法可以解决缝洞型油藏在非稳态水侵过程中定量水侵量的问题，计算方法繁琐，计算精度误高。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明的利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法一实施例的流程图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
参考图1，其为本发明的利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法一实施例的流程图。本发明的利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法包含如下步骤：
S1、获得油藏地质储量N、地层压力、地层压力为pi时的地层油体积系数Boi、累积产油量NP、地层压力为p时的地层油体积系数Bo、地层束缚水压缩系数Cw、地层束缚水的饱和度Swi、地层原始含油饱和度Soi、累积产水量Wp、水的体积系数Bw、环形水层函数QD、无因次时间tD、油层外半径r0、孔隙度
综合压缩系数ct以及水的粘度μw；其中，pi为初始地层压力，p为开采任意时间地层压力；其中，获得油藏地质储量N、地层压力、地层压力为pi时的地层油体积系数Boi、累积产油量NP、地层压力为p时的地层油体积系数Bo、地层束缚水压缩系数Cw、地层束缚水的饱和度Swi、地层原始含油饱和度Soi、累积产水量Wp、水的体积系数Bw、油层外半径r0、孔隙度
综合压缩系数ct以及水的粘度μw依据油气田日常计量数据获得。
S2、根据物质平衡方程的基本思想，在不考虑地层条件下溶解气体积的情况下可列出如下表达式：
整理后将各式对时间t进行求导有：
利用上述公式即可求出采液速度dWe/dt；其中，We为水侵量，t为时间，压差△P＝pi-p；
S3、由Schilthuis水侵模型可知，水侵速度dWe/dt与压差(pi-p)成正比。利用下述公式，根据初始地层压力、开采任意时间的地层压力计算出水侵系数k′：
S4、基于van Everdingen和Hurst所提出的理论，即在内边界压力恒定，无岩体弹性能释放，无水弹性能释放，无流固耦合，无井间干扰的情况下，水侵量可以写成下述公式。利用下述公式，根据环形水层函数、无因次时间、油层外半径、孔隙度、综合压缩系数、水的粘度，计算出水侵量：
其中：QD为环形水层函数；tD为无因次时间；r0为油层外半径；
为孔隙度；ct为综合压缩系数；μw为水的粘度，pj是表示时间tj时的地层压力，压力增量
pi＝p-1，时间j与时间j+1之间的时间段为无限小时间段，ΔtD表示在j时间步长下的无因次时间,n表示无限小时间段的总段数。
将步骤S3已求出的水侵系数k′代入本步骤中的上式即可求出每一个无限短时间内的水侵量。由于瞬时方程计算的是无限小段时间内的水侵量，所以依次叠加求出总水侵量即可。本领域人员可以理解的是，无限小时间段只存在于理论中，实际在进行本步骤时，无限小时间段只要取得足够小就行了，具体取多小本领域人员是可以判断得出的。
在油藏开发过程中，注采比是表征油田注水开发过程中注采平衡状况,反映产液量，注水量与地层压力之间联系的综合指标，且对油田开发措施及提高采收率有重要意义。边底水油藏注采比的确定需要有准确的水侵量为依托，从而推导边水油藏地层压力恢复速度与注采比的关系式,实现了对注采比的计算和合理的预测。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。
Claims (2)
1.一种利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，其特征在于，包含如下步骤：
S1、获得油藏地质储量N、地层压力、地层压力为pi时的地层油体积系数Boi、累积产油量NP、地层压力为p时的地层油体积系数Bo、地层束缚水压缩系数Cw、地层束缚水的饱和度Swi、地层原始含油饱和度Soi、累积产水量Wp、水的体积系数Bw、环形水层函数QD、无因次时间tD、油层外半径r0、孔隙度
综合压缩系数ct以及水的粘度μw；其中，pi为初始地层压力，p为开采任意时间地层压力；
S2、利用下述公式求出采液速度dWe/dt；
其中，We为水侵量，t为时间，压差△P＝pi-p；
S3、利用下述公式，根据初始地层压力、开采任意时间的地层压力计算出水侵系数k′：
S4、利用下述公式，根据环形水层函数、无因次时间、油层外半径、孔隙度、综合压缩系数、水的粘度，计算出水侵量：
其中：QD为环形水层函数；tD为无因次时间；r0为油层外半径；
为孔隙度；ct为综合压缩系数；μw为水的粘度，pj是表示时间tj时的地层压力，压力增量
pi＝p-1，时间j与时间j+1之间的时间段为无限小时间段，ΔtD表示在j时间步长下的无因次时间,n表示无限小时间段的总段数。
2.根据权利要求1所述的利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法，其特征在于，获得油藏地质储量N、地层压力、地层压力为pi时的地层油体积系数Boi、累积产油量NP、地层压力为p时的地层油体积系数Bo、地层束缚水压缩系数Cw、地层束缚水的饱和度Swi、地层原始含油饱和度Soi、累积产水量Wp、水的体积系数Bw、油层外半径r0、孔隙度
综合压缩系数ct以及水的粘度μw依据油气田日常计量数据获得。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910161366.XA CN109918769B (zh) | 2019-03-04 | 2019-03-04 | 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910161366.XA CN109918769B (zh) | 2019-03-04 | 2019-03-04 | 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109918769A CN109918769A (zh) | 2019-06-21 |
| CN109918769B true CN109918769B (zh) | 2020-11-27 |
Family
ID=66963127
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910161366.XA Expired - Fee Related CN109918769B (zh) | 2019-03-04 | 2019-03-04 | 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109918769B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112487594B (zh) * | 2019-08-23 | 2022-11-01 | 中国石油天然气股份有限公司 | 油藏水体倍数计算方法及装置 |
| CN110984974B (zh) * | 2019-11-29 | 2022-01-25 | 中国石油集团川庆钻探工程有限公司 | 基于水侵速度确定有水气藏合理采气速度的方法 |
| CN112329243B (zh) * | 2020-11-09 | 2023-08-25 | 中国海洋石油集团有限公司 | 非均质储层高含水水平井水侵过程反演方法和装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6842700B2 (en) * | 2002-05-31 | 2005-01-11 | Schlumberger Technology Corporation | Method and apparatus for effective well and reservoir evaluation without the need for well pressure history |
| CN104975827A (zh) * | 2014-04-03 | 2015-10-14 | 中国石油化工股份有限公司 | 预测二氧化碳驱油藏指标的物质平衡方法 |
| CN105156102A (zh) * | 2015-09-28 | 2015-12-16 | 中国石油大学(北京) | 底水油藏水体能量三维物理模拟装置及方法 |
| CN105464652A (zh) * | 2015-12-03 | 2016-04-06 | 恒泰艾普石油天然气技术服务股份有限公司 | 缝洞型碳酸盐岩油藏缝洞单元动态储量计算方法和系统 |
| CN105956712A (zh) * | 2016-05-18 | 2016-09-21 | 恒泰艾普（北京）能源科技研究院有限公司 | 缝洞型碳酸盐岩油藏见水预警方法 |
| CN108133080A (zh) * | 2017-12-04 | 2018-06-08 | 陕西延长石油(集团)有限责任公司研究院 | 一种考虑非稳态吸附的非均质裂缝页岩气藏数值模拟方法 |
| CN109339775A (zh) * | 2018-10-25 | 2019-02-15 | 西南石油大学 | 一种确定水驱气藏水体大小的方法 |
-
2019
- 2019-03-04 CN CN201910161366.XA patent/CN109918769B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6842700B2 (en) * | 2002-05-31 | 2005-01-11 | Schlumberger Technology Corporation | Method and apparatus for effective well and reservoir evaluation without the need for well pressure history |
| CN104975827A (zh) * | 2014-04-03 | 2015-10-14 | 中国石油化工股份有限公司 | 预测二氧化碳驱油藏指标的物质平衡方法 |
| CN105156102A (zh) * | 2015-09-28 | 2015-12-16 | 中国石油大学(北京) | 底水油藏水体能量三维物理模拟装置及方法 |
| CN105464652A (zh) * | 2015-12-03 | 2016-04-06 | 恒泰艾普石油天然气技术服务股份有限公司 | 缝洞型碳酸盐岩油藏缝洞单元动态储量计算方法和系统 |
| CN105956712A (zh) * | 2016-05-18 | 2016-09-21 | 恒泰艾普（北京）能源科技研究院有限公司 | 缝洞型碳酸盐岩油藏见水预警方法 |
| CN108133080A (zh) * | 2017-12-04 | 2018-06-08 | 陕西延长石油(集团)有限责任公司研究院 | 一种考虑非稳态吸附的非均质裂缝页岩气藏数值模拟方法 |
| CN109339775A (zh) * | 2018-10-25 | 2019-02-15 | 西南石油大学 | 一种确定水驱气藏水体大小的方法 |
Non-Patent Citations (2)
| Title |
|---|
| Jiong Wang;Weijia Ma;Jungcheng Liu;Yu Wang.The Prediction Methods of the Water Influx Intensity of the Non-.《2012 Fourth International Conference on Computational and Information Sciences》.2012,第1384-1388页. * |
| 计算水驱气藏动态储量和水侵量的简易新方法;胡俊坤，李晓平，张健涛，张修明;《天然气地球科学》;20121231;第3卷(第6期);第1175-1178页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109918769A (zh) | 2019-06-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109918769B (zh) | 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法 | |
| CN110905472B (zh) | 确定基于复合暂堵体系的实时转向压裂参数的方法 | |
| CN108868748A (zh) | 一种页岩气水平井重复压裂裂缝开启压力的计算方法 | |
| CN108518218B (zh) | 一种非常规油气藏多段压裂水平井单井动态储量确定方法 | |
| CN107066769B (zh) | 适用于超深层裂缝型碳酸盐岩储层的高效酸化设计方法 | |
| CN105626006A (zh) | 低渗透油藏co2驱技术极限井距确定方法 | |
| CN111794740B (zh) | 一种适用在缝洞型碳酸盐岩油藏动态储量计算的方法 | |
| CN110096718B (zh) | 一种获取碳酸盐岩油藏中溶洞的体积的方法 | |
| CN106677750A (zh) | 碳酸盐岩油藏溶洞‑裂缝储集体注水指示曲线解释模型 | |
| CN106545336A (zh) | 考虑致密气藏渗流机理的产能计算方法 | |
| CN103510940B (zh) | 机械采油井工况综合诊断分析方法及装置 | |
| CN106499370A (zh) | 井组同步分注分采的分层采油井各层段产液量计算方法及装置 | |
| CN104405374A (zh) | 一种致密气藏储层应力敏感性的测量方法 | |
| CN111950112B (zh) | 一种适用于底部封闭的碳酸盐岩储层动态分析方法 | |
| CN113486537B (zh) | 一种基于气水交替驱油的井距设计方法、装置及设备 | |
| CN102392636A (zh) | 水驱开发砂岩油藏地下流体通道参数的描述方法 | |
| CN111582532A (zh) | 应力敏感性油藏水平井产液能力预测方法及装置 | |
| CN109522579B (zh) | 水平井压裂施工破裂压力预测方法 | |
| CN112253102A (zh) | 油井套管放气压力的确定方法和装置 | |
| CN112377182A (zh) | 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置 | |
| CN111950111B (zh) | 一种适用于底部开放的碳酸盐岩储层动态分析方法 | |
| CN112035993B (zh) | 一种底部定压的碳酸盐岩储层测试评价方法 | |
| CN111963161A (zh) | 确定隐性不正常油井的方法及装置 | |
| Mott | Calculating well deliverability in gas condensate reservoirs | |
| CN109915122B (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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20221102 Address after: 100000 22 Chaoyangmen North Street, Chaoyang District, Beijing. Patentee after: CHINA PETROLEUM & CHEMICAL Corp. Patentee after: CHINA PETROLEUM & CHEMICAL CORPORATION EXPLORATION & PRODUCTION Research Institute Address before: 430000 Lu Mill Road, Hongshan District, Wuhan, Hubei Province, No. 388 Patentee before: CHINA University OF GEOSCIENCES (WUHAN CITY) |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201127 |