CN104179499A - 考虑油藏参数时变的数值模拟方法 - Google Patents

考虑油藏参数时变的数值模拟方法 Download PDF

Info

Publication number
CN104179499A
CN104179499A CN201310200461.9A CN201310200461A CN104179499A CN 104179499 A CN104179499 A CN 104179499A CN 201310200461 A CN201310200461 A CN 201310200461A CN 104179499 A CN104179499 A CN 104179499A
Authority
CN
China
Prior art keywords
numerical simulation
permeability
reservoir
model
oil
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.)
Pending
Application number
CN201310200461.9A
Other languages
English (en)
Inventor
曹绪龙
王延忠
谭保国
陈燕虎
许强
郭兰磊
唐从见
戴涛
张世明
董亚娟
宋勇
胡慧芳
侯玉培
王成锋
柳士成
黄迎松
郭士博
单联涛
王红艳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Shengli Geological Scientific Reserch Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Shengli Geological Scientific Reserch Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201310200461.9A priority Critical patent/CN104179499A/zh
Publication of CN104179499A publication Critical patent/CN104179499A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种考虑油藏参数时变的数值模拟方法,该方法包括:步骤1,建立油藏数值模拟模型,在油藏数值模拟的每一步迭代计算过程中,根据该油藏数值模拟模型每个网格块的属性参数,计算得到一组新的注水冲刷倍数的数据场;步骤2,建立储层渗透率随注水冲刷倍数的变化关系模型,在迭代计算过程中时时校正油藏模型静态参数的动态变化范围;步骤3,计算基于物性变化的相对渗透率曲线的时变规律,描述迭代计算中每个网格的相对渗透率曲线;以及步骤4,在获得该相对渗透率曲线后,用以每一迭代时间步内计算油水运动规律。该考虑油藏参数时变的数值模拟方法具有原理简单、描述更加准确、可操作性强等特点,因而具有很好的推广应用价值。

Description

考虑油藏参数时变的数值模拟方法
技术领域
本发明涉及油藏数值模拟应用及油气田开发提高油藏采收率领域,特别是涉及到一种考虑油藏参数时变的数值模拟方法。
背景技术
对于陆相沉积的多层砂岩整装油藏,储层非均质严重,纵向和平面矛盾突出,剩余油分布复杂,特别是对于经过长期注水冲刷,综合含水超过98%的多层系高采出程度油藏,油藏储层目前的非均质性及其流体的渗流规律都随注水冲刷而发生了很大的变化。常规油藏数值模拟研究油藏模拟已经作为油藏工程研究和管理的标准技术被业内认同,它可以辅助油藏工程师开展划分层系井网结构、确定注采方式及驱替机理、定量研究剩余油分布规律、研究剩余油挖潜措施方案、开展开发方案预测与优化研究及提高采收率研究等。但常规数值模拟方法没有反映高注入倍数特征条件下油藏储层参数及渗流规律的变化,从而不能精确定量表征强注水冲刷油藏剩余油分布,直接影响了油田大幅度提高油藏最终采收率的准确性。为此我们发明了一种新的考虑油藏参数时变的数值模拟方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种原理简单、描述更加准确、可操作性强的考虑油藏参数时变的数值模拟方法。   
本发明的目的可通过如下技术措施来实现:考虑油藏参数时变的数值模拟方法,该考虑油藏参数时变的数值模拟方法包括:步骤1,建立油藏数值模拟模型,在油藏数值模拟的每一步迭代计算过程中,根据该油藏数值模拟模型每个网格块的属性参数,计算得到一组新的注水冲刷倍数的数据场;步骤2,建立储层渗透率随注水冲刷倍数的变化关系模型,在迭代计算过程中时时校正油藏模型静态参数的动态变化范围;步骤3,计算基于物性变化的相对渗透率曲线的时变规律,描述迭代计算中每个网格的相对渗透率曲线;以及步骤4,在获得该相对渗透率曲线后,用以每一迭代时间步内计算油水
运动规律。
本发明的目的还可通过如下技术措施来实现:
在步骤1中,在时变油藏数值模拟研究过程中,定义该注水冲刷倍数为流过单位网格体积的累积注水量与该网格体积的比值。
该比值由每一步迭代以计算得到该组新的注水冲刷倍数的数据场。
在步骤2中,根据地质研究及动态分析,建立该储层渗透率随注水冲刷倍数的变化关系模型,以精确计算在高注水冲刷倍数下储层渗透率场的变化大小,从而在迭代计算过程中时时校正油藏模型静态参数的动态变化范围。
在步骤3中,根据实验及动态分析研究,建立束缚水饱和度变化值与渗透率倍数的关系模型,残余油饱和度变化值与渗透率倍数的关系模型,油相渗透率与绝对渗透率的关系模型,水相渗透率与油相渗透率的关系模型,以计算该基于物性变化的相对渗透率曲线的时变规律,描述迭代计算中每个网格的该相对渗透率曲线。
本发明中的考虑油藏参数时变的数值模拟方法,针对传统的油藏数值模拟方法没有考虑油藏物性参数及相对渗透率曲线的时变性,不能反映油藏特高含水期后期物性分布及流体渗流规律,提出了在油藏数值模拟计算过程中考虑注水冲刷倍数与渗透率及相对渗透率的变化关系,从而建立动态变化模型的处理方法。采用此改进的数值模拟方法,能够有效反映油藏储层参数随注水冲刷倍数变化而变化的现象,从而保证了油藏数值模拟结果的可靠性,进而为油田开发方案的调整与优化提供合理的参考依据。其在计算过程中修改储层参数随注水冲刷倍数的变化,精确描述储层参数及流体在特高含水后期高注入倍数情况下的变化规律,定量表征储层经过长期注水冲刷后剩余油分布,探索该类油藏特高含水后期大幅度提高采收率技术,为同类油藏进一步有效挖潜提供借鉴和指导。
附图说明
图1为本发明的考虑油藏参数时变的数值模拟方法的一具体实施例的流程图;
图2为本发明的一具体实施例中建立的渗透率变化倍数与注水冲刷孔隙体积倍数模型的示意图;
图3为本发明的一具体实施例中建立的束缚水饱和度与渗透率变化倍数模型的示意图;
图4为本发明的一具体实施例中建立的残余油饱和度与渗透率变化倍数模型的示意图;
图5为本发明的一具体实施例中建立的油相渗透率与渗透率变化倍数模型的示意图;
图6为本发明的一具体实施例中建立的水相渗透率与油相渗透率模型的示意图;
图7为本发明的—具体实施例中不考虑油藏参数时变计算的剩余油饱和度分布图;
图8为本发明的—具体实施例中考虑油藏参数时变计算的剩余油饱和度分布图。
 
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合所附图式,作详细说明如下。
如图1所示,图1为本发明的考虑油藏参数时变的数值模拟方法的流程图。
在步骤101,建立数值模拟模型,在油藏数值模拟的每一步迭代计算过程中,根据该数值模拟模型每个网格块的属性参数,计算得到一组新的注水冲刷倍数的数据场。也就是说,在时变油藏数值模拟研究过程中,引入注水冲刷倍数的概念。定义注水冲刷倍数为流过单位网格体积的累积注水量与该网格体积的比值。该值由每一步迭代计算得到一组新的注水冲刷倍数的数据场。流程进入到步骤102。
在步骤102,根据地质研究及动态分析,建立储层渗透率随注水冲刷倍数的变化关系模型,以精确计算在高注水冲刷倍数下储层渗透率场的变化大小,从而在迭代计算过程中时时校正油藏模型静态参数的动态变化范围。流程进入到步骤103。
在步骤103,根据实验及动态分析研究,建立束缚水饱和度变化值与渗透率倍数的关系模型,残余油饱和度变化值与渗透率倍数的关系模型,油相
渗透率与绝对渗透率的关系模型,水相渗透率与油相渗透率的关系模型,计算基于物性变化的相对渗透率曲线的时变规律,描述迭代计算中每个网格的相对渗透率曲线。流程进入到步骤104。
在步骤104,在已有油藏数值模拟软件的基础上,得到每一个迭代时间步的注水冲刷倍数数据场后,按照油藏开发的含水阶段进行分段,分阶段利用步骤102建立的渗透率随注水冲刷倍数的变化关系模型计算每个网格渗透率变化倍数,在模型中时时修改模型渗透率大小,然后根据步骤103建立的相对渗透率时变模型计算每个网格的相对渗透率曲线的变化,获得新的相对渗透率曲线后,用以每一迭代时间步内计算油水运动规律。流程结束。
本发明的考虑油藏参数时变的数值模拟方法,通过修改渗透率及相对渗透率曲线,数值模拟模型计算充分反映了储层参数及流体流动规律的变化,有效模拟了特高含水后期近极限含水条件下的水驱油运动规律,保证了特高含水后期油藏数值模拟的准确性,为开展剩余油分布研究与挖潜、开发调整方案的制定提供了有利的技术保障。
在应用本发明的考虑油藏参数时变的数值模拟方法的一具体实施例中,建立数值模拟概念模型,模型X方向45个网格,Y方向1个网格,Z方向10个网格。平面网格步长10米,纵向网格步长1米。模型孔隙度取0.28,均质模型渗透率2000×10-3μm2。初始含油饱和度0.67,初始含水饱和度0.33,地层压力取21.3MPa。注采单元为一注一采,注采比控制为1。
建立的渗透率变化倍数与注水冲刷孔隙体积倍数模型见图2,如图3到图6,图3为本发明的一具体实施例中建立的束缚水饱和度与渗透率变化倍数模型的示意图;图4为本发明的一具体实施例中建立的残余油饱和度与渗透率变化倍数模型的示意图;图5为本发明的一具体实施例中建立的油相渗透率与渗透率变化倍数模型的示意图;图6为本发明的一具体实施例中建立的水相渗透率与油相渗透率模型的示意图。
注采井分别位于模型两端。生产井以每天0.392立方米的液量进行定液量生产,最小井底流压为11MPa,注水井注水量按油藏注采比1进行软件自动计算,极限含水取0.98。在计算过程中考虑高注水冲刷倍数下渗透率和相对渗透率发生一定的变化。图7为本发明的—具体实施例中不考虑油藏参数时变计算的剩余油饱和度分布图,图8为本发明的—具体实施例中考虑油藏参
数时变计算的剩余油饱和度分布图。计算结果表明,当不考虑油藏参数时变时,纵向水驱油效率相对接近,纵向下部大部分剩余油接近残余油饱和度,上部只有局部存在剩余油富集。当考虑油藏参数时变时,储层底部渗透率随注水冲刷倍数增大而变大,相对渗透率端点残余油饱和度降低,注入水主要沿储层底部窜流到生产井,纵向上驱油效率差别相对较大,顶部剩余油富集范围相对较大。研究结果表明,考虑油藏参数时变对剩余油分布影响非常大,需要在实际油藏数值模拟研究中加以重视,以获得更加准确的研究结论,指导油田开发方案的调整与优化。
本发明中的考虑油藏参数时变的数值模拟方法,能够很好地反映随着开发的进行储层参数的变化及相对渗透率的端点变化对油水运动规律的影响,从而保证了特高含水后期油藏数值模拟的准确性,进而为油田开发方案的调整与优化提供合理的参考依据。

Claims (5)

1.考虑油藏参数时变的数值模拟方法,其特征在于,该考虑油藏参数时变的数值模拟方法包括:
步骤1,建立油藏数值模拟模型,在油藏数值模拟的每一步迭代计算过程中,根据该油藏数值模拟模型每个网格块的属性参数,计算得到一组新的注水冲刷倍数的数据场;
步骤2,建立储层渗透率随注水冲刷倍数的变化关系模型,在迭代计算过程中时时校正油藏模型静态参数的动态变化范围;
步骤3,计算基于物性变化的相对渗透率曲线的时变规律,描述迭代计算中每个网格的相对渗透率曲线;以及
步骤4,在获得该相对渗透率曲线后,用以每一迭代时间步内计算油水运动规律。
2.根据权利要求1所述的考虑油藏参数时变的数值模拟方法,其特征在于,在步骤1中,在时变油藏数值模拟研究过程中,定义该注水冲刷倍数为流过单位网格体积的累积注水量与该网格体积的比值。
3.根据权利要求2所述的考虑油藏参数时变的数值模拟方法,其特征在于,该比值由每一步迭代以计算得到该组新的注水冲刷倍数的数据场。
4.根据权利要求1所述的考虑油藏参数时变的数值模拟方法,其特征在于,在步骤2中,根据地质研究及动态分析,建立该储层渗透率随注水冲刷倍数的变化关系模型,以精确计算在高注水冲刷倍数下储层渗透率场的变化大小,从而在迭代计算过程中时时校正油藏模型静态参数的动态变化范围。
5.根据权利要求1所述的考虑油藏参数时变的数值模拟方法,其特征在于,在步骤3中,根据实验及动态分析研究,建立束缚水饱和度变化值与渗透率倍数的关系模型,残余油饱和度变化值与渗透率倍数的关系模型,油相渗透率与绝对渗透率的关系模型,水相渗透率与油相渗透率的关系模型,以计算该基于物性变化的相对渗透率曲线的时变规律,描述迭代计算中每个网格的该相对渗透率曲线。
CN201310200461.9A 2013-05-27 2013-05-27 考虑油藏参数时变的数值模拟方法 Pending CN104179499A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310200461.9A CN104179499A (zh) 2013-05-27 2013-05-27 考虑油藏参数时变的数值模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310200461.9A CN104179499A (zh) 2013-05-27 2013-05-27 考虑油藏参数时变的数值模拟方法

Publications (1)

Publication Number Publication Date
CN104179499A true CN104179499A (zh) 2014-12-03

Family

ID=51960812

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310200461.9A Pending CN104179499A (zh) 2013-05-27 2013-05-27 考虑油藏参数时变的数值模拟方法

Country Status (1)

Country Link
CN (1) CN104179499A (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089582A (zh) * 2015-05-28 2015-11-25 中国石油天然气股份有限公司 基于井下流量控制设备的油藏数值模拟方法及装置
CN107762495A (zh) * 2016-08-18 2018-03-06 中国石油化工股份有限公司 特高含水期油藏模型纵向网格尺寸的优化方法
CN108197366A (zh) * 2017-12-26 2018-06-22 中国石油大学(华东) 考虑注入水水质对储层伤害的注采参数优化方法
CN108487904A (zh) * 2018-02-06 2018-09-04 中国石油大学(华东) 一种基于图版的消除末端效应的相渗曲线校正方法
CN109268005A (zh) * 2018-10-30 2019-01-25 中国石油大学(华东) 一种基于储层时变性的剩余油预测方法及工业化流程
CN109522583A (zh) * 2018-08-13 2019-03-26 唐山冀油瑞丰化工有限公司 一种用于多层油藏开发层系组合的非均质性界限确定方法
CN110424944A (zh) * 2019-07-22 2019-11-08 中国海洋石油集团有限公司 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法
CN110984979A (zh) * 2020-01-06 2020-04-10 中国石油化工股份有限公司 针对角度不整合油藏建模数模的网格设计方法
CN111191183A (zh) * 2020-01-09 2020-05-22 中国石油大学(北京) 一种储层参数计算方法、装置及设备
CN111220509A (zh) * 2020-01-20 2020-06-02 中国石油天然气股份有限公司 一种考虑渗透率时变性的油水相对渗透率曲线修正方法
CN112560310A (zh) * 2020-12-11 2021-03-26 中国海洋石油集团有限公司 一种多层油藏渗透率动态变化定量表征的方法
CN112796749A (zh) * 2021-03-10 2021-05-14 青岛石大恒瑞石油新技术开发有限公司 一种基于毛管数的油藏水驱数值模拟方法
CN113674100A (zh) * 2020-05-14 2021-11-19 中国石油化工股份有限公司 油藏注采优化方法、装置、存储介质及电子设备
CN114673494A (zh) * 2022-05-06 2022-06-28 重庆科技学院 蒸汽吞吐后储层孔渗预测方法
CN115526067A (zh) * 2022-11-28 2022-12-27 中国石油大学(华东) 一种水驱油藏储层参数时变规律预测方法及***
CN116859021A (zh) * 2023-07-06 2023-10-10 西南石油大学 一种通过产能损伤率评价储层敏感性的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060047432A1 (en) * 2004-09-02 2006-03-02 Patrick Egermann Method of determining multiphase flow parameters of a porous medium Taking account of the local heterogeneity
CN101661514A (zh) * 2008-05-21 2010-03-03 中国石化股份胜利油田分公司地质科学研究院 一种油藏黑油模型数值模拟***
CN101942994A (zh) * 2010-09-16 2011-01-12 中国石油天然气股份有限公司 水淹层产水率定量预测方法及其***

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060047432A1 (en) * 2004-09-02 2006-03-02 Patrick Egermann Method of determining multiphase flow parameters of a porous medium Taking account of the local heterogeneity
CN101661514A (zh) * 2008-05-21 2010-03-03 中国石化股份胜利油田分公司地质科学研究院 一种油藏黑油模型数值模拟***
CN101942994A (zh) * 2010-09-16 2011-01-12 中国石油天然气股份有限公司 水淹层产水率定量预测方法及其***

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刘显太: "中高渗透砂岩油藏储层物性时变数值模拟技术", 《油气地质与采收率》 *
崔传智等: "水驱油藏高含水期渗透率的动态分布计算模型及应用", 《中国石油大学学报( 自然科学版)》 *
邴绍献: "水驱油藏单井可采储量影响因素权重定量研究", 《油气地质与采收率》 *

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105089582A (zh) * 2015-05-28 2015-11-25 中国石油天然气股份有限公司 基于井下流量控制设备的油藏数值模拟方法及装置
CN105089582B (zh) * 2015-05-28 2017-12-05 中国石油天然气股份有限公司 基于井下流量控制设备的油藏数值模拟方法及装置
CN107762495A (zh) * 2016-08-18 2018-03-06 中国石油化工股份有限公司 特高含水期油藏模型纵向网格尺寸的优化方法
CN108197366A (zh) * 2017-12-26 2018-06-22 中国石油大学(华东) 考虑注入水水质对储层伤害的注采参数优化方法
CN108487904A (zh) * 2018-02-06 2018-09-04 中国石油大学(华东) 一种基于图版的消除末端效应的相渗曲线校正方法
CN109522583A (zh) * 2018-08-13 2019-03-26 唐山冀油瑞丰化工有限公司 一种用于多层油藏开发层系组合的非均质性界限确定方法
CN109522583B (zh) * 2018-08-13 2023-10-31 唐山冀油瑞丰化工有限公司 一种用于多层油藏开发层系组合的非均质性界限确定方法
CN109268005A (zh) * 2018-10-30 2019-01-25 中国石油大学(华东) 一种基于储层时变性的剩余油预测方法及工业化流程
CN110424944B (zh) * 2019-07-22 2022-08-23 中国海洋石油集团有限公司 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法
CN110424944A (zh) * 2019-07-22 2019-11-08 中国海洋石油集团有限公司 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法
CN110984979A (zh) * 2020-01-06 2020-04-10 中国石油化工股份有限公司 针对角度不整合油藏建模数模的网格设计方法
CN111191183A (zh) * 2020-01-09 2020-05-22 中国石油大学(北京) 一种储层参数计算方法、装置及设备
CN111220509A (zh) * 2020-01-20 2020-06-02 中国石油天然气股份有限公司 一种考虑渗透率时变性的油水相对渗透率曲线修正方法
CN113674100A (zh) * 2020-05-14 2021-11-19 中国石油化工股份有限公司 油藏注采优化方法、装置、存储介质及电子设备
CN113674100B (zh) * 2020-05-14 2024-04-30 中国石油化工股份有限公司 油藏注采优化方法、装置、存储介质及电子设备
CN112560310A (zh) * 2020-12-11 2021-03-26 中国海洋石油集团有限公司 一种多层油藏渗透率动态变化定量表征的方法
CN112560310B (zh) * 2020-12-11 2022-08-30 中国海洋石油集团有限公司 一种多层油藏渗透率动态变化定量表征的方法
CN112796749A (zh) * 2021-03-10 2021-05-14 青岛石大恒瑞石油新技术开发有限公司 一种基于毛管数的油藏水驱数值模拟方法
CN114673494A (zh) * 2022-05-06 2022-06-28 重庆科技学院 蒸汽吞吐后储层孔渗预测方法
CN115526067A (zh) * 2022-11-28 2022-12-27 中国石油大学(华东) 一种水驱油藏储层参数时变规律预测方法及***
CN116859021A (zh) * 2023-07-06 2023-10-10 西南石油大学 一种通过产能损伤率评价储层敏感性的方法

Similar Documents

Publication Publication Date Title
CN104179499A (zh) 考虑油藏参数时变的数值模拟方法
Lie et al. Successful application of multiscale methods in a real reservoir simulator environment
CN103939066B (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN104750896B (zh) 一种缝洞型碳酸盐岩油藏数值模拟方法
CN105740563B (zh) 一种成熟油田二次开发之优势通道识别方法
CN104318052B (zh) 一种低渗透油藏水驱波及系数的评价方法
CN105089582B (zh) 基于井下流量控制设备的油藏数值模拟方法及装置
Durlofsky et al. Uncertainty quantification for subsurface flow problems using coarse-scale models
CN103400020B (zh) 一种测算多条相交离散裂缝流动状况的油藏数值模拟方法
Kruger Determining areal permeability distribution by calculations
CN106484933A (zh) 一种用于确定页岩气井井控动态储量的方法及***
CN105401927B (zh) 油藏井网井距确定方法及装置
CN103577886A (zh) 一种低渗气藏水平井分段压裂产能预测方法
CN103279647A (zh) 一种基于孔隙特征参数计算致密砂岩储层渗透率的方法
CN104453834A (zh) 一种井组注采关系优化调整方法
Yao et al. Fractured vuggy carbonate reservoir simulation
CN105089566A (zh) 一种气井***配产方法
CN104915530B (zh) 油藏井间连通关系的建立方法
CN103912248A (zh) 水驱油田预测含水率方法
Baker Streamline technology: reservoir history matching and forecasting= its success, limitations, and future
CN105631078A (zh) 天然裂缝性油藏自适应介质的数值模拟方法
Izgec et al. Maximizing volumetric sweep efficiency in waterfloods with hydrocarbon F–Φ curves
CN106014365B (zh) 一种预测水驱开发油田产量递减率的方法
CN110863818B (zh) 一种剩余油/气分布的描述方法及装置
Lian et al. Evaluation and applicability study on prediction methods of water inflow in mines

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20141203