一种复杂边水油藏的水体分布智能识别方法

摘要

本发明公开了一种复杂边水油藏的水体分布智能识别方法，属于油气田开发技术领域；它解决了现今复杂边水油藏因对边水分布情况认识不清楚而导致无法提出针对性的治水对策等问题；其技术方案是：基于单井地质数据和生产资料，建立可模拟复杂边水水侵流入动态的水侵单元数值模拟模型，结合遗传算法，修正水体单元的特征参数，包含水体单元体积和水体单元的水侵量，以实际的生产动态数据自动拟合模型计算的生产动态数据，得到最优水体单元的特征参数解，并赋值给水体单元，反演确定水体分布情况。本发明水体分布智能识别方法步骤简单，水体分布反演结果与数值模拟器的结果对比后，证明了该方法的准确性。

权利要求

1.一种复杂边水油藏的水体分布智能识别方法，其特征在于，该方法包括下列步骤：S100、水侵单元数值模型建立步骤，以水体单元‑产井单元‑加密单元为基础，采用物质平衡方程和水驱前缘推进方程计算各单元的压力、饱和度和含水率，建立水侵单元数值模拟模型；

S200、水体特征参数自动拟合步骤，在水侵单元数值模型的基础上，以实际的生产动态数据为参考，通过遗传算法自动修正水体单元体积和水体单元水侵量两个特征参数使模型预测得到的数据与实际的生产动态数据相匹配，修正后的水体单元体积和水体单元水侵量作为模型自动拟合后的水体特征参数；

S300、水体分布反演步骤，在得到水体特征参数最优解后，将水体单元体积和水体单元水侵量两个特征参数赋值到被离散的水体单元中，构建水体单元柱状图，通过每个单元的柱状高度，确定水体分布情况。

2.根据权利要求1所述的一种复杂边水油藏的水体分布智能识别方法，其特征在于：所述步骤S100具体包括以下步骤：S101、将边水油藏进行简化表征，将边部水体离散成与产井单元相连通的水体单元，将地质认识清晰的位置离散成加密单元，构建以水体单元‑产井单元‑加密单元为基础的简化边水油藏单元系统，并用水体单元体积和水体单元水侵量两个参数表征水体单元；

S102、以水体单元‑产井单元‑加密单元为基础，采用物质平衡方程计算各单元的压力；

S103、根据计算得到的各单元压力，采用前缘推进方程计算各单元的饱和度和含水率。

3.根据权利要求1所述的一种复杂边水油藏的水体分布智能识别方法，其特征在于：所述步骤S200具体包括以下步骤：S201、根据建立的水侵单元数值模型建立水体单元的特征参数向量和目标函数；

S202、设置水体单元体积的约束条件，认为水体单元的体积等于整个边水水体的大小；

S203、采用遗传算法自动拟合水体单元体积和水体单元水侵量两个特征参数，直到模型预测数据与实际生产动态数据误差值最小，得到修正后的水体单元体积和水体单元水侵量两个水体特征参数。

4.根据权利要求1或2所述的一种复杂边水油藏的水体分布智能识别方法，其特征在于：所述是S100中的物质平衡方程为 ，3

其中， 为水体单元i与产井单元j之间在t时刻的平均传导率，m /（d·MPa）；nwell为产井单元及加密单元数，无因次； 为水体单元i在t时刻的水体压力，MPa； 为产井单元j3

在t时刻的井底平均压力，MPa； 为水体单元i在t时刻的水侵量，m /d；ct为地层综合压‑1 3

缩系数，MPa ； 为水体单元i在t时刻的水体体积，m。

5.根据权利要求1或2所述的一种复杂边水油藏的水体分布智能识别方法，其特征在于：所述S100中的水驱前缘推进方程为 ，其中，x为水侵通道的3

某一位置，m；Q为流体的总流量，m ；Sw为位置x处的含水饱和度，无因次； 为水相的分2

流量，无因次；A为渗流截面积，m； 为孔隙度，无因次。

6.根据权利要求1所述的一种复杂边水油藏的水体分布智能识别方法，其特征在于：所述的生产动态数据为单井的井底压力和单井的含水率。

说明书

一种复杂边水油藏的水体分布智能识别方法

技术领域

[0001] 本发明涉及一种复杂边水油藏的水体分布智能识别方法，属于油气田开发技术领域。

背景技术

[0002] 对于存在边水的油藏而言，会在油藏开发阶段由于水侵的影响造成油藏水淹。生产井在见水后会影响整体的生产制度，从而影响油藏的最终采收率。因此能够准确地认识
油藏水体分布特征对于油藏选择合理的治水措施、开发方案调整具有重要的指导意义。在
水体特征研究方面，许多学者通过水体体积和水侵量等特征参数来表征水体特征，但是这
些研究对于水体的分布情况无法表征，因而无法满足对发生水侵的油藏进行针对性治理与
调整的需求。

[0003] 目前，部分学者能够在构建物理模型方法和数值模拟方法中表征出水体的分布位置。但两种方法仍然具有一定的局限性，对于构建物理模型方法来讲，油藏水侵过程属于大
范围的流体运动，物理模型实现与原型的相似存在着难以克服的问题；随着计算机技术的
发展，数值模拟方法适用范围越来越广，计算效率也越来越高，且成本较低，许多学者会采
用数值模拟方法模拟水侵动态过程，目前数值模拟方法模拟水侵动态过程多是基于地质资
料构建水体，通过历史拟合验证模型的准确性。在这过程中存在两方面问题：（1）对于地质
资料缺乏或者存在测量误差的数值模型，建立的水体模型结果说服力不足；（2）模型历史拟
合过程中，可调整的参数过多，无法充分说明建立的水体模型的准确性。

[0004] 总的来说，目前关于复杂边水油藏水体分布刻画的方法均存在一定的的局限性。因此，迫切地需要一种适用于复杂边水油藏的水体分布智能识别方法。

发明内容

[0005] 本发明目的是：为了解决复杂边水油藏因水体分布认识不清楚而造成油藏水侵过程中无法开展针对性治水对策，从而影响油藏最终采收率的问题，基于单井地质数据和生
产资料，在数值模拟方法和智能优化算法的基础上反演出水体特征，实现对水体分布的再
识别。

[0006] 为实现上述目的，本发明提供了一种复杂边水油藏的水体分布智能识别方法，该方法包括下列步骤：

[0007] 第一，水侵单元数值模型建立步骤，以水体单元‑产井单元‑加密单元为基础，采用物质平衡方程和水驱前缘推进方程计算各单元的压力、饱和度和含水率，建立水侵单元数
值模拟模型；

[0008] 第二，水体特征参数自动拟合步骤，在水侵单元数值模型的基础上，以实际的生产动态数据为参考，通过遗传算法自动修正水体单元体积和水体单元水侵量两个特征参数使
模型预测得到的数据与实际的生产动态数据相匹配，修正后的水体单元体积和水体单元水
侵量作为模型自动拟合后的水体特征参数；

[0009] 第三，水体分布反演步骤，在得到水体特征参数最优解后，将水体单元体积和水体单元水侵量两个特征参数赋值到被离散的水体单元中，构建水体单元柱状图，通过每个单
元的柱状高度，确定水体分布情况。

[0010] 上述一种复杂边水油藏的水体分布智能识别方法中，水侵单元数值模型的建立具体包括下列几个步骤：

[0011] 第一，将边水油藏进行简化表征，将边部水体离散成与产井单元相连通的水体单元，将地质认识清晰的位置离散成加密单元，构建以水体单元‑产井单元‑加密单元为基础
的简化边水油藏单元系统，并用水体单元体积和水体单元水侵量两个参数表征水体单元；

[0012] 第二，以水体单元‑产井单元‑加密单元为基础，采用物质平衡方程计算各单元的压力，物质平衡方程表示为：

[0013]3

[0014] 其中， 为水体单元i与产井单元j之间在t时刻的平均传导率，m /（d·MPa）；nwell为产井单元及加密单元数，无因次； 为水体单元i在t时刻的水体压力，MPa； 为
3
产井单元j在t时刻的井底平均压力，MPa； 为水体单元i在t时刻的水侵量，m/d；ct为地
‑1 3
层综合压缩系数，MPa ； 为水体单元i在t时刻的水体体积，m ；

[0015] 第三，采用有限差分方法求解物质平衡方程从旧时间步t到新时间步t+1推进的压力解，构造的压力矩阵为：

[0016]

[0017] 其中， 为单元n为研究对象时的其他单元，包括水体单元i和产井单元j，无量纲； 代表单元n在t时刻的流量，由于单元n包含水体单元i和产井单元j，因此 值为正
3
数时，等于水侵量，为负数时，等于生产井产量，m/d； 表示时间步长，d；

[0018] 第四，根据计算得到的各单元压力，采用水驱前缘推进方程计算各单元的饱和度和含水率，水驱前缘推进方程为：

[0019]

[0020] 其中，x为水侵通道的某一位置，m；Q为流体的总流量，m3；Sw为位置 处的含水饱和2
度，无因次； 为水相的分流量，无因次；A为渗流截面积，m； 为孔隙度，无因次；

[0021] 上述一种复杂边水油藏的水体分布智能识别方法中，水体特征参数自动拟合具体包括下列几个步骤：

[0022] 第一，根据建立的水侵单元数值模型建立水体单元的特征参数向量和目标函数，建立的特征参数向量包含水体单元体积参数和水体单元水侵量参数，表达式为：
，建立的目标函数表
达式为：

[0023] 其中，m为水体单元的特征参数向量；Y(m)为水体单元的目标函数；ydata为实测动态数据；F(m)为模型修正水体单元的特征参数计算得到的动态数据；

[0024] 第二，设置水体单元体积的约束条件，认为水体单元的体积等于整个边水水体的大小；

[0025] 第三，采用遗传算法自动拟合水体单元体积和水体单元水侵量两个特征参数，直到模型预测数据与实际生产动态数据误差值最小，得到修正后的水体单元体积和水体单元
水侵量两个水体特征参数。

附图说明

[0026] 图1是本方法的技术路线图。

[0027] 图2是水体分布反演示意图。

[0028] 图3是建立的数值模拟相渗曲线图。

[0029] 图4是凹型水体分布示意图。

[0030] 图5a是P1井的含水率拟合曲线图。

[0031] 图5b是P2井的含水率拟合曲线图。

[0032] 图5c是P3井的含水率拟合曲线图。

[0033] 图5d是P4井的含水率拟合曲线图。

[0034] 图5e是P5井的含水率拟合曲线图。

[0035] 图5f是P6井的含水率拟合曲线图。

[0036] 图6是模型自动拟合后得到的水体分布反演结果图。

[0037] 图7是模型计算得到的含水饱和度分布图。

具体实施方式

[0038] 下面结合实施方式和附图对本发明做进一步说明。

[0039] 本发明提供了一种复杂边水油藏的水体分布智能识别方法，图1为本发明的技术路线图，该方法包括下列步骤：

[0040] 第一，水侵单元数值模型建立步骤，以水体单元‑产井单元‑加密单元为基础，采用物质平衡方程和水驱前缘推进方程计算各单元的压力、饱和度和含水率，建立水侵单元数
值模拟模型；

[0041] 第二，水体特征参数自动拟合步骤，在水侵单元数值模型的基础上，以实际的生产动态数据为参考，通过遗传算法自动修正水体单元体积和水体单元水侵量两个特征参数使
模型预测得到的数据与实际的生产动态数据相匹配，修正后的水体单元体积和水体单元水
侵量作为模型自动拟合后的水体特征参数；

[0042] 第三，水体分布反演步骤，在得到水体特征参数最优解后，将水体单元体积和水体单元水侵量两个特征参数赋值到被离散的水体单元中，构建水体单元柱状图，通过每个单
元的柱状高度，确定水体分布情况，最终的水体分布反演结果如图2所示。

[0043] 进一步地，水侵单元数值模型的建立，其具体包括下列几个步骤：

[0044] 第一，将边水油藏进行简化表征，将边部水体离散成与产井单元相连通的水体单元，将地质认识清晰的位置离散成加密单元，构建以水体单元‑产井单元‑加密单元为基础
的简化边水油藏单元系统，并用水体单元体积和水体单元水侵量两个参数表征水体单元；

[0045] 第二，以水体单元‑产井单元‑加密单元为基础，采用物质平衡方程计算各单元的压力，物质平衡方程表示为：

[0046]

[0047] 其中， 为水体单元i与产井单元j之间在t时刻的平均传导率，m3/（d·MPa）；nwell为产井单元及加密单元数，无因次； 为水体单元i在t时刻的水体压力，MPa；
3
为产井单元j在t时刻的井底平均压力，MPa； 为水体单元i在t时刻的水侵量，m /d；ct
‑1 3
为地层综合压缩系数，MPa ； 为水体单元i在t时刻的水体体积，m ；

[0048] 第三，采用有限差分方法求解物质平衡方程从旧时间步t到新时间步t+1推进的压力解，构造的压力矩阵为：

[0049]

[0050] 其中， 为单元n为研究对象时的其他单元，包括水体单元i和产井单元j，无量纲； 代表单元n在t时刻的流量，由于单元n包含水体单元i和产井单元j，因此 值为
3
正数时，等于水侵量，为负数时，等于生产井产量，m/d； 表示时间步长，d；

[0051] 第四，根据计算得到的各单元压力，采用水驱前缘推进方程计算各单元的饱和度和含水率，水驱前缘推进方程为：

[0052]3

[0053] 其中，x为水侵通道的某一位置，m；Q为流体的总流量，m ；Sw为位置 处的含水饱和2
度，无因次； 为水相的分流量，无因次；A为渗流截面积，m ； 为孔隙度，无因次；

[0054] 上述一种复杂边水油藏的水体分布智能识别方法中，水体特征参数自动拟合具体包括下列几个步骤：

[0055] 第一，根据建立的水侵单元数值模型建立水体单元的特征参数向量和目标函数，建立的特征参数向量包含水体单元体积参数和水体单元水侵量参数，表达式为：
，建立的目标函数表
达式为：

[0056] 其中，m为水体单元的特征参数向量；Y(m)为水体单元的目标函数；ydata为实测动态数据；F(m)为模型修正水体单元的特征参数计算得到的动态数据；

[0057] 第二，设置水体单元体积的约束条件，认为水体单元的体积等于整个边水水体的大小；

[0058] 第三，采用遗传算法自动拟合水体单元体积和水体单元水侵量两个特征参数，直到模型预测数据与实际生产动态数据误差值最小，得到修正后的水体单元体积和水体单元
水侵量两个水体特征参数。

[0059] 将本发明所述方法所建模型应用到具体实例中如下：

[0060] 应用本文所建模型，笔者将水侵动态特征的计算结果与数值模拟器得到的结果进行对比验证。验证数值模拟模型采用黑油模型，模型网格设置为50×50×1，每个网格块横
3
向步长设置为10m，纵向步长设置为20m，模拟时长为2000天，原油密度为776kg/m ，原油粘
‑1
度为20cp，原油压缩系数为0.005MPa ，设置油藏的相渗曲线如图3所示，数值模拟模型设置
7 3
数值边水水体大小为1.8×10m ，水体分布类型为凹型水体分布，如图4所示，设置产井定液
3 3 3 3 3 3
量生产，分别为P1：100m /d，P2：50m/d，P3：100m/d，P4：200m /d，P5：100m/d，P6：200m/d，
模拟生产2000d。

[0061] 提取2000d的生产数据，将数值模拟器预测得到的生产井P1‑P6的含水率作为模型的目标函数，进行水体单元体积及水体单元水侵速度两项特征参数的历史拟合，可以看出
在凹型水体模型验证中，含水率曲线拟合效果较好，生产井的见水时间和最终含水率一致，
如图5a‑图5f所示。

[0062] 根据拟合结果，得到了凹型水体单元体积和水体单元水侵速度的特征值，组合这些特征参数，可以直观的看出凹型边水水体分布情况，如图6所示。在拟合得到的水体单元
体积和水体单元水侵速度的特征参数的基础上，模型反演得到了水侵动态特征，如图7所
示。