CN105239976B - 考虑水驱油藏渗流参数动态变化的数值模拟方法 - Google Patents
考虑水驱油藏渗流参数动态变化的数值模拟方法 Download PDFInfo
- Publication number
- CN105239976B CN105239976B CN201510750095.3A CN201510750095A CN105239976B CN 105239976 B CN105239976 B CN 105239976B CN 201510750095 A CN201510750095 A CN 201510750095A CN 105239976 B CN105239976 B CN 105239976B
- Authority
- CN
- China
- Prior art keywords
- permeability
- flux
- reservoir
- face flux
- numerical simulation
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种考虑水驱油藏渗流参数动态变化的数值模拟方法，包括如下步骤：S1、建立油藏数值模拟模型，在迭代计算过程中，首先计算求解得到每个网格块的压力和饱和度数据场；S2、计算每个网格块的水相流量，计算每个网格块的方向面通量和总面通量，得到方向面通量和总面通量数据场；S3，建立储层渗透率随方向面通量的变化关系模型，对渗透率数据场进行更新；S4，建立储层相对渗透率曲线随总面通量的变化关系模型，对相对渗透率曲线的端点值进行重新标定；S5，在获得修改后的储层渗透率和相对渗透率曲线，进行下一个时间步的循环计算。本发明可为准确预测剩余油分布、定量计算可动剩余油、制定油田开发方案、指导油田挖潜潜力评价提供理论支持。
Description
技术领域
本发明涉及油藏数值模拟应用及油气田开发提高油藏采收率领域，其针对水驱油藏渗流参数动态变化进行定量评价，形成了一套新的数值模拟方法。
背景技术
长期水驱条件下储层宏观、微观结构及流体成分特征不断变化，影响渗流的关键参数——储层渗透率、原油粘度、岩石润湿性都将发生较大变化，这些物性参数的变化无疑对剩余油的分布及油田最终开发效果具有重要影响。油藏数值模拟技术，由于对实际储层非均质性有比较好的刻画，因而是目前描述剩余油分布预测的最好方法。但是，目前国内外开发的油藏数值模拟软件的数值模型中，一般都认为注水过程中油藏的储层参数不发生变化，考虑储层参数变化的数值模型中，普遍存在储层参数变化不能连续表征、不能方向性表征、受网格划分大小影响等问题。因此，该技术对储层参数时变特别是连续时变的过程难以客观的表达，预测潜力区剩余油宏观分布存在较大误差，因而油藏物理特征参数时变带来的一系列问题面临的挑战，
上述论述内容目的在于向读者介绍可能与下面将被描述和/或主张的本发明的各个方面相关的技术的各个方面，相信该论述内容有助于为读者提供背景信息，以有利于更好地理解本发明的各个方面，因此，应了解是以这个角度来阅读这些论述，而不是承认现有技术。
发明内容
本发明的目的在于避免现有技术中的不足而提供一种描述合理、可靠性强、操作简单的考虑水驱油藏渗流参数动态变化的数值模拟方法。
本发明的目的通过以下技术方案实现：
提供一种考虑水驱油藏渗流参数动态变化的数值模拟方法，包括如下步骤：
S1、建立油藏数值模拟模型，在油藏数值模拟的每一步迭代计算过程中，首先计算求解得到每个网格块的压力和饱和度数据场；
S2、计算每个网格块的各方向水相流量，计算每个网格块的方向面通量和总面通量，得到方向面通量数据场、总面通量数据场；
S3，建立储层渗透率随方向面通量的变化关系模型，对渗透率数据场进行更新；
S4，建立储层相对渗透率曲线随总面通量的变化关系模型，对相对渗透率曲线的端点值进行重新标定；
S5，在获得修改后的储层渗透率和相对渗透率曲线后，进行下一个时间步的循环计算。
其中，在步骤S2中，面通量指累计通过单位面积的水相体积，方向面通量指累计通过X、Y、Z各方向上的单位面积的水相体积，总面通量指累计通过X、Y、Z方向上的单位面积的总水相体积。
其中，在步骤S3中，根据矿场取心井统计和室内试验研究，建立该储层渗透率随方向面通量的变化关系模型，定量表征不同方向面通量冲刷下的渗透率变化范围，从而在迭代计算过程中更新不同方向上的渗透率数据场。
其中，在步骤S4中，根据矿场取心井统计和室内试验研究，建立该储层相对渗透率曲线随总面通量的变化关系模型，对不同面通量冲刷下的相对渗透率曲线重新进行端点标定，描述迭代计算中每个网格的相对渗透率曲线。
本发明的有益效果：本发明所述考虑水驱油藏渗流参数动态变化的数值模拟方法，选取面通量表征水驱冲刷强度，建立储层渗流参数随面通量的变化规律，对每一个迭代步的储层参数动态变化进行更新，可以实现对油藏参数动态变化的稳定、准确、连续性表征，保证了油藏数值模拟结果的稳定性和可靠性，为准确预测剩余油分布，定量计算可动剩余油，制定油田开发方案，指导油田挖潜潜力评价提供理论支持。
附图说明
利用附图对本发明作进一步说明，但附图中的实施例不构成对本发明的任何限制，对于本领域的普通技术人员，在不付出创造性劳动的前提下，还可以根据以下附图获得其它的附图。
图1为本发明的流程图。
图2为本发明建立的渗透率变化倍数与方向面通量模型的示意图。
图3为本发明建立的相渗端点值与总面通量模型的示意图。
图4为本发明的一具体实例中不考虑水驱油藏渗流参数动态变化的剩余油饱和度分布图。
图5为本发明的一具体实例中考虑水驱油藏渗流参数动态变化的剩余油饱和度分布图。
图6为本发明的一具体实例中不考虑和考虑水驱油藏渗流参数动态变化的采收率对比图。
具体实施方式
为了使本领域的技术人员更好地理解本发明的技术方案，下面结合附图和具体实施例对本发明作进一步详细的描述，需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
如图1所示，本发明所述的考虑水驱油藏渗流参数动态变化的数值模拟方法，包括以下步骤：S1、建立油藏数值模拟模型，在油藏数值模拟的每一步迭代计算过程中，首先计算求解得到每个网格块的压力和饱和度数据场。S2、计算每个网格块的各方向水相流量，计算每个网格块的方向面通量和总面通量，得到方向面通量数据场、总面通量数据场；在步骤S2中，面通量指累计通过单位面积的水相体积，方向面通量指累计通过X、Y、Z各方向上的单位面积的水相体积，总面通量指累计通过X、Y、Z方向上的单位面积的总水相体积。S3，建立储层渗透率随方向面通量的变化关系模型，对渗透率数据场进行更新；在步骤S3中，根据矿场取心井统计和室内试验研究，建立该储层渗透率随方向面通量的变化关系模型，定量表征不同方向面通量冲刷下的渗透率变化范围，从而在迭代计算过程中更新不同方向上的渗透率数据场，即传导率数据场。S4，建立储层相对渗透率曲线(相渗曲线)随总面通量的变化关系模型，对相对渗透率曲线的端点值进行重新标定；在步骤S4中，根据矿场取心井统计和室内试验研究，建立该储层相对渗透率曲线随总面通量的变化关系模型，对不同面通量冲刷下的相对渗透率曲线重新进行端点标定，描述迭代计算中每个网格的相对渗透率曲线。S5，在获得修改后的储层渗透率和相对渗透率曲线后，进行下一个时间步的循环计算。
为使本发明的操作流程、发明意义及必要性更加通俗易懂，现举出一实施概念模型为例，并配合附图，作详细说明如下。
建立油藏数值模拟模型，模型X方向25个网格，Y方向25个网格，Z方向5个网格。平面网格步长10米，纵向网格步长2米。模型孔隙度取0.25，平面方向渗透率均为100mD，垂直方向渗透率10mD，初始含油饱和度0.7，初始含水饱和度0.3，注采单元为一注一采，定压差生产10年。
建立的渗透率变化倍数与方向面通量模型见图2；建立的相渗端点值与总面通量模型见图3，在计算过程中考虑高面通量的冲刷作用下储层渗透率和相对渗透率曲线发生变化，图4为不考虑油藏参数时变计算的剩余油饱和度分布图，图5为考虑油藏参数时变计算的剩余油饱和度分布图，图6为不考虑和考虑水驱油藏渗流参数动态变化的采收率对比图。
计算结果表明，当不考虑油藏参数时变时，如图4所示，主流线附近的含油饱和度比较低，大部分剩余油接近残余油饱和度，开发效果最好，而非主流线区域冲刷比较少，含油饱和度比较高，开发效果差。当考虑水驱油藏渗流参数动态变化时，随着面通量的增大，储层渗透率逐渐变大，相对渗透率端点残余油饱和度降低，对应的水相渗透率增加，如图5所示，导致主流线附近的含油饱和度不断降低，非主流线区域的开发效果也变好，剩余油富集部位大大减小，整体采收率增加，如图6所示。研究结果表明，考虑水驱油藏渗流参数动态变化对驱油效率、剩余油分布、开发效果的影响非常大，必须在实际油藏数值模拟研究中加以重视，以准确预测剩余油分布，加深对水驱油藏开发规律的认识，指导油田开发方案的调整与优化。
本发明所述的考虑水驱油藏渗流参数动态变化的数值模拟方法，选取面通量表征水驱冲刷强度，建立储层渗流参数随面通量的变化规律，对每一个迭代步的储层参数动态变化进行更新，可以实现对油藏参数动态变化的稳定、准确、连续性表征，比较客观反映长期水驱油藏渗流参数的变化和油水运动规律的影响，保证了油藏数值模拟结果的稳定性和可靠性，为准确预测剩余油分布，定量计算可动剩余油，制定油田开发方案，指导油田挖潜潜力评价提供理论支持。
上面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，不能理解为对本发明保护范围的限制。
总之，本发明虽然例举了上述优选实施方式，但是应该说明，虽然本领域的技术人员可以进行各种变化和改型，除非这样的变化和改型偏离了本发明的范围，否则都应该包括在本发明的保护范围内。
Claims (4)
1.一种考虑水驱油藏渗流参数动态变化的数值模拟方法，其特征在于，包括如下步骤：
S1、建立油藏数值模拟模型，在油藏数值模拟的每一步迭代计算过程中，首先计算求解得到每个网格块的压力和饱和度数据场；
S2、计算每个网格块的各方向水相流量，计算每个网格块的方向面通量和总面通量，得到方向面通量数据场、总面通量数据场；
S3，建立储层渗透率随方向面通量的变化关系模型，对渗透率数据场进行更新；
S4，建立储层相对渗透率曲线随总面通量的变化关系模型，对相对渗透率曲线的端点值进行重新标定；
S5，在获得修改后的储层渗透率和相对渗透率曲线后，进行下一个时间步的循环计算。
2.根据权利要求1所述的考虑水驱油藏渗流参数动态变化的数值模拟方法，其特征在于：在步骤S2中，面通量指累计通过单位面积的水相体积，方向面通量指累计通过X、Y、Z各方向上的单位面积的水相体积，总面通量指累计通过X、Y、Z方向上的单位面积的总水相体积。
3.根据权利要求1所述的考虑水驱油藏渗流参数动态变化的数值模拟方法，其特征在于：在步骤S3中，根据矿场取心井统计和室内试验研究，建立该储层渗透率随方向面通量的变化关系模型，定量表征不同方向面通量冲刷下的渗透率变化范围，从而在迭代计算过程中更新不同方向上的渗透率数据场。
4.根据权利要求1所述的考虑水驱油藏渗流参数动态变化的数值模拟方法，其特征在于：在步骤S4中，根据矿场取心井统计和室内试验研究，建立该储层相对渗透率曲线随总面通量的变化关系模型，对不同面通量冲刷下的相对渗透率曲线重新进行端点标定，描述迭代计算中每个网格的相对渗透率曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510750095.3A CN105239976B (zh) | 2015-11-04 | 2015-11-04 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510750095.3A CN105239976B (zh) | 2015-11-04 | 2015-11-04 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105239976A CN105239976A (zh) | 2016-01-13 |
| CN105239976B true CN105239976B (zh) | 2018-01-23 |
Family
ID=55037770
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510750095.3A Active CN105239976B (zh) | 2015-11-04 | 2015-11-04 | 考虑水驱油藏渗流参数动态变化的数值模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105239976B (zh) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107762495A (zh) * | 2016-08-18 | 2018-03-06 | 中国石油化工股份有限公司 | 特高含水期油藏模型纵向网格尺寸的优化方法 |
| CN107102968A (zh) * | 2017-03-28 | 2017-08-29 | 中国石油大学(华东) | 一种非均质致密油藏渗流的时间尺度分析方法 |
| CN109577942B (zh) * | 2017-09-27 | 2022-07-19 | 中国石油化工股份有限公司 | 一种优势渗流通道发育油藏剩余油的挖潜方法 |
| CN110020453B (zh) * | 2018-01-09 | 2023-03-10 | 姚君波 | 地层渗透率三维分解模拟器 |
| CN108487904B (zh) * | 2018-02-06 | 2020-05-19 | 中国石油大学(华东) | 一种基于图版的消除末端效应的相渗曲线校正方法 |
| CN108825217B (zh) * | 2018-04-19 | 2021-08-20 | 中国石油化工股份有限公司 | 适用于油藏数值模拟的综合井指数计算方法 |
| CN108665534B (zh) * | 2018-05-10 | 2022-05-13 | 昆仑数智科技有限责任公司 | 一种注水油藏流线数值模拟的在线三维可视化方法 |
| CN109268005A (zh) * | 2018-10-30 | 2019-01-25 | 中国石油大学(华东) | 一种基于储层时变性的剩余油预测方法及工业化流程 |
| CN109356564B (zh) * | 2018-11-08 | 2021-11-23 | 中国海洋石油集团有限公司 | 一种基于三元动态参数的水驱油相渗时变规律表征方法 |
| CN109858177B (zh) * | 2019-02-23 | 2022-06-07 | 中国海洋石油集团有限公司 | 一种基于拟流线法的底水油藏水平井水驱数值模拟方法 |
| CN112069701B (zh) * | 2019-05-22 | 2022-10-04 | 中国石油天然气股份有限公司 | 油藏开发方法及装置 |
| CN112230278A (zh) * | 2019-07-15 | 2021-01-15 | 中国石油天然气集团有限公司 | 渗流场特征参数确定方法及装置 |
| CN110424944B (zh) * | 2019-07-22 | 2022-08-23 | 中国海洋石油集团有限公司 | 巨厚碳酸盐岩油藏油水拟相渗曲线的建立方法 |
| CN110593849B (zh) * | 2019-08-22 | 2022-07-12 | 魏峰 | 一种考虑相渗曲线动态变化的油藏数值模拟方法 |
| CN112990631B (zh) * | 2019-12-17 | 2024-03-29 | 中国石油天然气股份有限公司 | 渗流场的评价方法及装置、计算机可读存储介质 |
| CN111428375A (zh) * | 2020-03-31 | 2020-07-17 | 中国石油大学(华东) | 一种高含水油田层系重组划分新方法 |
| CN112560357A (zh) * | 2020-10-21 | 2021-03-26 | 中海石油(中国)有限公司 | 一种考虑聚合物驱过程相渗时变的油藏数值模拟方法 |
| CN112560310B (zh) * | 2020-12-11 | 2022-08-30 | 中国海洋石油集团有限公司 | 一种多层油藏渗透率动态变化定量表征的方法 |
| CN112685970B (zh) * | 2020-12-25 | 2022-08-30 | 中国石油大学(华东) | 一种水驱油藏流动单元渗流界面定量表征方法及系统 |
| CN112945829B (zh) * | 2021-02-07 | 2023-05-26 | 西安石油大学 | 一种致密砂岩储层水驱剩余油分析方法及系统 |
| CN112796749B (zh) * | 2021-03-10 | 2023-05-09 | 青岛石大恒瑞石油新技术开发有限公司 | 一种基于毛管数的油藏水驱数值模拟方法 |
| CN113255247B (zh) * | 2021-06-04 | 2022-03-25 | 中国石油大学(华东) | 一种高含水期油藏多尺度数值模拟方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2259473C2 (ru) * | 2003-09-30 | 2005-08-27 | Открытое акционерное общество "Татнефть" им. В.Д. Шашина | Способ разработки обводненной нефтяной залежи (варианты) |
| CN102562048A (zh) * | 2010-12-30 | 2012-07-11 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN103573234A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种确定完整油水相对渗透率曲线的方法 |
| CN103821485A (zh) * | 2014-03-17 | 2014-05-28 | 中国石油大港油田勘探开发研究院 | 预测水驱油田含水上升率变化的方法 |
| CN103924966A (zh) * | 2013-01-10 | 2014-07-16 | 中国石油化工股份有限公司 | 基于储层物性时变模型的层系井网再建方法 |
| CN104179498A (zh) * | 2013-05-22 | 2014-12-03 | 中国石油化工股份有限公司 | 分带模拟饱和度方法 |
| CN104834807A (zh) * | 2015-03-18 | 2015-08-12 | 成都北方石油勘探开发技术有限公司 | 一种基于分形理论的应力敏感储层相对渗透率计算方法 |
-
2015
- 2015-11-04 CN CN201510750095.3A patent/CN105239976B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2259473C2 (ru) * | 2003-09-30 | 2005-08-27 | Открытое акционерное общество "Татнефть" им. В.Д. Шашина | Способ разработки обводненной нефтяной залежи (варианты) |
| CN102562048A (zh) * | 2010-12-30 | 2012-07-11 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN103573234A (zh) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | 一种确定完整油水相对渗透率曲线的方法 |
| CN103924966A (zh) * | 2013-01-10 | 2014-07-16 | 中国石油化工股份有限公司 | 基于储层物性时变模型的层系井网再建方法 |
| CN104179498A (zh) * | 2013-05-22 | 2014-12-03 | 中国石油化工股份有限公司 | 分带模拟饱和度方法 |
| CN103821485A (zh) * | 2014-03-17 | 2014-05-28 | 中国石油大港油田勘探开发研究院 | 预测水驱油田含水上升率变化的方法 |
| CN104834807A (zh) * | 2015-03-18 | 2015-08-12 | 成都北方石油勘探开发技术有限公司 | 一种基于分形理论的应力敏感储层相对渗透率计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105239976A (zh) | 2016-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105239976B (zh) | 考虑水驱油藏渗流参数动态变化的数值模拟方法 | |
| CN102507412B (zh) | 碳酸盐岩油藏等效模型的裂缝-基质渗透率级差判别方法 | |
| CN103590818B (zh) | 一种井筒内温度分布半解析确定方法及装置 | |
| CN104318052B (zh) | 一种低渗透油藏水驱波及系数的评价方法 | |
| CN105740563B (zh) | 一种成熟油田二次开发之优势通道识别方法 | |
| CN103149600B (zh) | 一种基于优化控制点确定水文地质参数的自动配线方法 | |
| CN103577886A (zh) | 一种低渗气藏水平井分段压裂产能预测方法 | |
| CN103399970B (zh) | 用离散裂缝线处理进行数模测算油藏流动状况的方法 | |
| CN105401927A (zh) | 油藏井网井距确定方法及装置 | |
| CN103867175A (zh) | 一种蒸汽驱井网结构及其蒸汽驱开发方法 | |
| CN104346511A (zh) | 一种全新的油藏动态监测方法及装置 | |
| CN110593849B (zh) | 一种考虑相渗曲线动态变化的油藏数值模拟方法 | |
| CN105095642A (zh) | 油藏数值模拟中历史拟合调整参数确定方法及装置 | |
| CN104933268A (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
| CN105631138A (zh) | 油气藏开发单元及井组数值模拟动态分析方法 | |
| CN110439502A (zh) | 一种油藏水驱带分析方法 | |
| CN110083852B (zh) | 井底压力的确定方法、装置及存储介质 | |
| CN108595782B (zh) | 一种离散裂缝中基质与裂缝间的传质计算方法 | |
| CN110990924B (zh) | 一种注浆工程加固路基浆液扩散范围的确定方法 | |
| CN106401549A (zh) | 裂缝与井轨迹夹角对水平井生产效果影响的图版建立方法 | |
| Li et al. | A semi‐elliptical surface compound diffusion model for synchronous grouting filling stage in specially shaped shield tunnelling | |
| CN107288595A (zh) | 一种注水利用率的评价方法 | |
| CN105302979A (zh) | 两相流体网络模型中阀门组的建模方法和系统 | |
| CN104747154A (zh) | 一种利用驱油效率比提高蒸汽驱剩余油研究精度的方法 | |
| CN106600075A (zh) | 一种考虑地层各向异性的水平井产能预测新方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Applicant after: China Offshore Oil Group Co., Ltd. Applicant after: CNOOC (China) Limited Zhanjiang Branch Applicant after: QINGDAO SHIDA HENGRUI PETROLEUM NEW TECHNOLOGY DEVELOPMENT CO., LTD. Address before: China CNOOC Building No. 25 Beijing 100000 Chaoyang District Chaoyangmen North Street Applicant before: China National Offshore Oil Corporation Applicant before: CNOOC (China) Limited Zhanjiang Branch Applicant before: QINGDAO SHIDA HENGRUI PETROLEUM NEW TECHNOLOGY DEVELOPMENT CO., LTD. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |