CN115600509A - 一种气窜通道体积计算方法及系统 - Google Patents

一种气窜通道体积计算方法及系统 Download PDF

Info

Publication number
CN115600509A
CN115600509A CN202211495725.3A CN202211495725A CN115600509A CN 115600509 A CN115600509 A CN 115600509A CN 202211495725 A CN202211495725 A CN 202211495725A CN 115600509 A CN115600509 A CN 115600509A
Authority
CN
China
Prior art keywords
production
production well
well
gas
tracer
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.)
Granted
Application number
CN202211495725.3A
Other languages
English (en)
Other versions
CN115600509B (zh
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 University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN202211495725.3A priority Critical patent/CN115600509B/zh
Publication of CN115600509A publication Critical patent/CN115600509A/zh
Application granted granted Critical
Publication of CN115600509B publication Critical patent/CN115600509B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Software Systems (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Medical Informatics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及油气田开发技术领域，目的是提供一种气窜通道体积计算方法及系统，首先根据目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料时的示踪剂前缘推进速度，以使目标区块中所有生产井均具有示踪剂前缘推进速度，然后将目标区块的地质资料和各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到目标区块中各生产井方向的日注入气体积；根据各生产井示踪剂前缘推进速度和日注入气体积，计算出目标区块中各生产井的气窜通道体积。本发明全面考虑了注气过程中各种开发因素和复杂的地质因素带来的影响，可精确计算目标区块中各生产井的气窜通道体积。

Description

一种气窜通道体积计算方法及系统
技术领域
本发明涉及油气田开发技术领域，特别是涉及一种气窜通道体积计算方法与系统。
背景技术
注气提高采收率以其驱油效率高等优点被广泛应用于油田开发过程中。然而，由于油气间的流度差异，注入气易于在原油中发生黏性指进从而过早地到达生产井形成气窜，造成油田开发效果变差。因此，合理调控气窜通道、抑制注入气快速窜进对提高注入气波及程度、改善气驱开发效果至关重要。气窜通道体积是其封堵方案设计中的重要参数。准确计算气窜通道体积对气窜调控方案的设计具有重要意义。
现有技术中，专利CN113863920A公开了一种气窜通道体积监测方法，将油藏气驱开采时产生的气窜通道近似简化为圆形毛细管束模型进而对气窜通道体积进行计算。但是该方法仅能对一注一采条件下的气窜通道体积进行计算，无法应用于实际矿场中多注多采的情况。专利CN113653475A公开了一种泡沫驱气窜通道的治理与利用方法，其中涉及使用采出井的示踪剂浓度及井筒参数计算气窜通道体积的方法。然而由于井间示踪剂测试成本较高，现场实际的示踪剂资料较为有限，仅有部分生产井会进行示踪剂检测，因此该方法无法应用于现场所有井的气窜通道体积计算。专利CN111852417A公开了一种二氧化碳驱油藏气窜通道快速反演方法和装置，但该方法使用简单的算术平均将注气井的注气量劈分到不同生产井的方向，没有考虑实际油藏中地质因素和开发因素对气体流动过程的影响，因此计算结果存在较大偏差。
发明内容
本发明的目的是提供一种准确性高的气窜通道体积计算方法及系统。
为实现上述目的，一种气窜通道体积计算方法，包括以下步骤：
获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度；所述第一生产井为所述目标区块中具有示踪剂资料的生产井；
根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度；所述第二生产井为所述目标区块中不具有示踪剂资料的生产井；
将所述目标区块的地质资料和所述目标区块中各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到所述目标区块中各生产井方向的日注入气体积；所述日注入气体积预测模型为以所述区块的地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下油藏数值模拟得到的各生产井方向的日注入气体积作为标签训练得到的模型；
根据所述目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到所述目标区块中各生产井的气窜通道体积。
可选的，根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度，具体包括：
根据所述第一生产井的示踪剂前缘推进速度和各所述第一生产井与注气井之间的井间连通性系数，确定示踪剂前缘推进速度与对应井间连通性系数之间的拟合关系式；
将所述第二生产井的井间连通性系数代入所述拟合关系式，得到所述第二生产井对应的示踪剂前缘推进速度。
可选的，所述方法还包括对所述日注入气体积预测模型进行训练，具体为：
获取区块数据资料，所述区块数据资料包括所述目标区块的地质资料以及历史生产动态资料；
构建注入气体积样本库，所述注入气体积样本库的构建包括：
根据所述地质资料构建地质模型；
根据所述历史生产动态资料设置注采参数进行油藏数值模拟，得到所述地质资料和所述历史生产动态资料对应条件下的各生产井方向的日注入气体积；
根据所述地质资料、所述历史生产动态资料和对应条件下所述各生产井方向的日注入气体积，构成注入气体积样本库中的一组数据；
以所述注入气体积样本数据库中所述地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下各生产井方向的日注入气体积作为标签，对日注入气体积预测模型进行训练，得到训练好的日注入气体积预测模型。
可选的，在所述获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度之前，还包括：
获取所述目标区块的地质资料和生产动态资料；
根据所述目标区块的地质资料和生产动态资料，计算所述目标区块中各生产井与注气井之间的连通性系数。
可选的，所述生产动态资料包括目标区块中各口注气井的注气量、各口生产井的产量和井底流压。
可选的，所述地质资料包括目标区块的孔隙度、渗透率、油层厚度、地层压力、油气最小混相压力、油气粘度比和相对渗透率曲线。
可选的，所述目标区块中各生产井与注气井的连通性系数的数学模型如下：
Figure 791958DEST_PATH_IMAGE001
Figure 54312DEST_PATH_IMAGE002
Figure 93681DEST_PATH_IMAGE003
其中，
Figure 520115DEST_PATH_IMAGE004
为数学模型求解的目标函数，λ j 为第j口生产井与注气井之间的连通性系数；N为生产井的个数；q j (t)为公式计算得到的第j口生产井在t时刻的日产液量；
Figure 306062DEST_PATH_IMAGE005
为实际生产动态资料中第j口生产井在t时刻的日产液量；
Figure 422922DEST_PATH_IMAGE006
为第j口生产井的时间常数；t 0为初始生产时刻；t为当前生产时刻；i(ξ)为实际生产动态资料中注气井在ξ时刻对应的地下日注气量；v j 为第j口生产井对应的采油指数；
Figure 462553DEST_PATH_IMAGE007
为第j口生产井在t时刻的井底流压。
可选的，根据所述各生产井的示踪剂前缘推进速度和日注入气体积，计算得到各生产井的气窜通道体积，具体包括：
根据A j =Q j /v j 计算各生产井方向的气窜通道截面积；
其中，A j 为第j口生产井方向的气窜通道截面积；Q j 为第j口生产井方向的日注入气体积；v j 为第j口生产井方向的示踪剂前缘推进速度；
根据V j =A j L j 计算气窜通道体积；
其中，V j 为第j口生产井的气窜通道体积；L j 为第j口生产井对应的注采井距。
本发明还提供了一种气窜通道体积计算系统，包括：
获取模块，用于获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度；所述第一生产井为所述目标区块中具有示踪剂资料的生产井；
示踪剂前缘推进速度计算模块，用于根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度；所述第二生产井为所述目标区块中不具有示踪剂资料的生产井；
日注入气体积计算模块，用于将所述目标区块的地质资料和所述目标区块中各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到所述目标区块中各生产井方向的日注入气体积；所述日注入气体积预测模型为以所述区块的地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下油藏数值模拟得到的各生产井方向的日注入气体积作为标签训练得到的模型；
气窜通道体积计算模块，用于根据所述目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到所述目标区块中各生产井的气窜通道体积。
可选的，所述示踪剂前缘推进速度计算模块，具体包括：
拟合关系式确定单元，用于根据所述第一生产井的示踪剂前缘推进速度和各所述第一生产井与注气井之间的井间连通性系数，确定示踪剂前缘推进速度与对应井间连通性系数之间的拟合关系式；
示踪剂前缘推进速度计算单元，用于将所述第二生产井的井间连通性系数代入所述拟合关系式，得到所述第二生产井对应的示踪剂前缘推进速度。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明提供了一种气窜通道体积计算方法及系统，其中获取目标区块中各生产井与注气井之间的连通性系数，并根据现场已有示踪剂资料的生产井的示踪剂前缘推进速度以及对应目标区块中各生产井与注气井之间的连通性系数，预测出不具有示踪剂资料的生产井的示踪剂前缘推进速度，通过将目标区块的地质资料和生产动态资料输入至训练好的日注入气体积预测模型，得到目标区块中各生产井方向的日注入气体积，根据目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到目标区块中各生产井的气窜通道体积。由于本发明对不具有示踪剂资料的生产井进行示踪剂前缘推进速度的预测，且目标区块中各生产井方向的日注入气体积是通过训练好的日注入气体积预测模型计算得到。因此，本发明能够对现场所有生产井的气窜通道体积进行计算，同时还可以全面考虑注气过程中各种地质因素和开发因素的影响，在一注一采的条件和多注多采的条件下，均能够对目标区块的气窜通道体积进行精确计算，使本发明提供的气窜通道体积计算方法具有准确度高、适用性强的特性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明提供的气窜通道体积计算方法流程图；
图2为本发明提供的示踪剂前缘推进速度和井间连通性系数关系图；
图3为本发明提供的不同生产井气窜通道体积计算时对应计算区域的示意图；
图4为本发明提供的XGBoost模型预测得到的不同生产井方向的日注入气体积和实际日注入气体积关系图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种准确性高的气窜通道体积计算方法及系统。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
本申请提供的气窜通道体积计算方法包括以下步骤：
步骤1：获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度；所述第一生产井为所述目标区块中具有示踪剂资料的生产井。
步骤2：根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度；所述第二生产井为所述目标区块中不具有示踪剂资料的生产井。
步骤3：将所述目标区块的地质资料和所述目标区块中各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到所述目标区块中各生产井方向的日注入气体积；所述日注入气体积预测模型为以所述区块的地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下油藏数值模拟得到的各生产井方向的日注入气体积作为标签训练得到的模型。
步骤4：根据所述目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到所述目标区块中各生产井的气窜通道体积。
在一些实施例中，在步骤1之前，还需要对目标区块中各生产井与注气井的连通性系数进行计算。
其中，生产动态资料可以包括注气井的注气量、生产井的产气量和井底流压、示踪剂资料等，地质资料可以包括孔隙度、渗透率、油层厚度、地层压力、油气最小混相压力、油气粘度比和相对渗透率曲线等。相对渗透率曲线可以直接得到残余油饱和度、束缚水饱和度。
利用生产动态资料中的注气井注气量、生产井产气量和井底流压进行各生产井与注气井之间动态连通性的分析，具体数学模型为：
Figure 15763DEST_PATH_IMAGE008
Figure 884362DEST_PATH_IMAGE009
Figure 340882DEST_PATH_IMAGE010
其中，
Figure 738366DEST_PATH_IMAGE011
为数学模型求解的目标函数，λ j 为第j口生产井与注气井之间的连通性系数；N为生产井的个数；q j (t)为公式计算得到的第j口生产井在t时刻的日产液量；
Figure 781802DEST_PATH_IMAGE012
为实际生产动态资料中第j口生产井在t时刻的日产液量；
Figure 204824DEST_PATH_IMAGE006
为第j口生产井的时间常数；t 0为初始生产时刻；t为当前生产时刻；i(ξ)为实际生产动态资料中注气井在ξ时刻对应的地下日注气量；v j 为第j口生产井对应的采油指数；
Figure 499539DEST_PATH_IMAGE013
为第j口生产井在t时刻的井底流压。
具体的，可以采用生物地理学优化算法（BBO）对上述连通性系数进行优化求解，本申请采用的优化算法中，以初始种群规模设为5，最大迭代次数设为12为例，求解得到不同生产井与注气井之间的连通性系数如表1所示。
表1 不同生产井与注气井间的连通性系数
Figure 582770DEST_PATH_IMAGE014
表2 第一生产井示踪剂前缘推进速度
Figure 392463DEST_PATH_IMAGE015
表2展示了一个具体示例中第一生产井的示踪剂前缘推进速度数据。通过获取第一生产井的示踪剂前缘推进速度以及第一生产井与注气井间的连通性系数，绘制第一生产井的示踪剂前缘推进速度和对应的井间连通性系数的散点图，如图2所示，并使用线性关系对其进行拟合，得到目标区块示踪剂前缘推进速度和井间连通性系数的关系式，在本示例中，该关系式为v j =27.8326λ j +0.7454，然后将第二生产井中的连通性系数代入目标区块示踪剂前缘推进速度和井间连通性系数的关系式，从而得到第二生产井的示踪剂前缘推进速度，如表3所示。
表3 第二生产井示踪剂前缘推进速度
Figure 291280DEST_PATH_IMAGE016
在一些实施例中，步骤3中日注入气体积预测模型的训练步骤如下：
获取区块数据资料，所述区块数据资料是根据目标区块的地质资料和历史生产动态资料中数据的分布范围，通过采样获得的具有目标区块特征的数据样本。
构建注入气体积样本库，所述的注入气体积样本库构建包括根据上述区块数据样本库中的地质资料构建地质模型，再根据上述区块数据样本库中的历史生产动态资料设置注采参数进行油藏数值模拟，进而得到该地质条件和生产动态条件下的各生产井方向的日注入气体积，地质条件、生产动态条件和各生产井方向的日注入气体积一同构成注入气体积样本库中的一组数据。
以所述注入气体积样本数据库中的地质资料、生产动态资料为输入，以所述注入气体积样本库中各生产井方向的日注入气体积作为标签，对日注入气体积预测模型进行训练，得到训练好的日注入气体积预测模型。
其中，目标区块的地质资料包括孔隙度、渗透率、油层厚度、地层压力、油气最小混相压力以及油气粘度比，目标区块中各生产井与注气井的历史生产动态资料包括注气井的注气量、生产井的产液量和井底流压等。
区块数据样本库构建时的采样方式可以选取拉丁超立方采样；根据区块数据样本库中的地质资料构建地质模型，其中构建地质模型的软件可以是斯伦贝谢公司的Petrel软件、贝克休斯公司的Jewel Suite软件或其他；得到油藏数值模拟结果后，根据
Figure 706081DEST_PATH_IMAGE017
计算不同生产井方向的日注入气体积V j ，其中，I，J，K分别对应数值模拟网格中x，y，z方向的位置；S g I,J,K 为（I，J，K）网格位置处的含气饱和度；Φ I,J,K 网格位置处的孔隙度；Δx，Δy，Δz分别是数值模拟中x，y，z方向的网格长度；Δt为数值模拟的天数；S j 为第j口生产井对应的计算区域，详见图3。
在一些实施例中，日注入气体积预测模型可以采用XGBoost模型。
模型采用Python中的XGBoost包进行构建，调节模型参数使得预测误差最小，最终构建模型所用的参数具体如下：迭代次数500；学习率0.08；损失函数类型RMSE；最大树深度6；最大叶子数量36；随机样本比例0.8。预测得到的不同生产井方向的日注入气体积和实际日注入气体积如图4所示。从图4中可以看出预测结果与实际结果吻合较好（分布在45°线附近），证实了构建的XGBoost模型的准确性。
本实施例中，步骤3得到的各生产井方向的日注入气体积计算结果如表4所示。
表4 实际区块中各生产井方向的日注入气体积
Figure 217004DEST_PATH_IMAGE018
本实施例中，步骤4采用下式对气窜通道截面积进行计算，计算结果如表5所示。
A j =Q j /v j， 其中，Aj为第j口生产井方向的气窜通道截面积；Qj为第j口生产井方向的日注入气体积；vj为第j口生产井方向的示踪剂前缘推进速度。
表5 目标区块中各生产井方向气窜通道截面积
Figure 920518DEST_PATH_IMAGE019
不同生产井与注入井之间的距离如表6所示，不同生产井的气窜通道体积采用下式计算：
V j =A j L j
其中，V j 表示第j口生产井的气窜通道体积；L j 表示第j口生产井对应的注采井距。
表6 目标区块中各生产井与注气井的距离
Figure 669031DEST_PATH_IMAGE020
将注采井间距离代入上述公式计算得到目标区块中各生产井的气窜通道体积，如表7所示。
表7 目标区块中各生产井对应的气窜通道体积
Figure 157912DEST_PATH_IMAGE021
本发明还提供了一种气窜通道体积计算系统，包括：
获取模块，用于获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度；所述第一生产井为所述目标区块中具有示踪剂资料的生产井。
示踪剂前缘推进速度计算模块，用于根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度；所述第二生产井为所述目标区块中不具有示踪剂资料的生产井。
日注入气体积计算模块，用于将所述目标区块的地质资料和所述目标区块中各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到所述目标区块中各生产井方向的日注入气体积；所述日注入气体积预测模型为以所述区块的地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下油藏数值模拟得到的各生产井方向的日注入气体积作为标签训练得到的模型。
气窜通道体积计算模块，用于根据所述目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到所述目标区块中各生产井的气窜通道体积。
所述示踪剂前缘推进速度计算模块，具体包括：
拟合关系式确定单元，用于根据所述第一生产井的示踪剂前缘推进速度和各所述第一生产井与注气井之间的井间连通性系数，确定示踪剂前缘推进速度与对应井间连通性系数之间的拟合关系式。
示踪剂前缘推进速度计算单元，用于将所述第二生产井的井间连通性系数代入所述拟合关系式，得到所述第二生产井对应的示踪剂前缘推进速度。
综上所述，本发明具有以下优点：利用已有示踪剂资料的生产井的示踪剂前缘推进速度以及井间动态连通性反演结果预测未进行示踪剂检测的生产井的示踪剂前缘推进速度，解决了由于示踪剂检测成本高、流程复杂而导致的仅有部分生产井有示踪剂资料的问题，进而可对目标区块中的各生产井进行气窜通道体积的计算。并且通过将目标区块的地质资料和生产动态资料输入至训练好的日注入气体积预测模型，可以全面地考虑注气过程中各种地质因素和开发因素的影响，在此基础上，利用XGBoost得到的日注入气体积准确性高，而且预测模型构建完成后可以在目标区块反复使用，适用性强。将二者结合极大程度地提高了气窜通道体积计算的精度，且方便简单，便于现场应用。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述得比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种气窜通道体积计算方法，其特征在于，包括以下步骤：
获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度；所述第一生产井为所述目标区块中具有示踪剂资料的生产井；
根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度；所述第二生产井为所述目标区块中不具有示踪剂资料的生产井；
将所述目标区块的地质资料和所述目标区块中各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到所述目标区块中各生产井方向的日注入气体积；所述日注入气体积预测模型为以所述区块的地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下油藏数值模拟得到的各生产井方向的日注入气体积作为标签训练得到的模型；
根据所述目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到所述目标区块中各生产井的气窜通道体积。
2.根据权利要求1所述的气窜通道体积计算方法，其特征在于，根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度，具体包括：
根据所述第一生产井的示踪剂前缘推进速度和各所述第一生产井与注气井之间的井间连通性系数，确定示踪剂前缘推进速度与对应井间连通性系数之间的拟合关系式；
将所述第二生产井的井间连通性系数代入所述拟合关系式，得到所述第二生产井对应的示踪剂前缘推进速度。
3.根据权利要求1所述的气窜通道体积计算方法，其特征在于，所述方法还包括对所述日注入气体积预测模型进行训练，具体为：
获取区块数据资料，所述区块数据资料包括所述目标区块的地质资料以及历史生产动态资料；
构建注入气体积样本库，所述注入气体积样本库的构建包括：
根据所述地质资料构建地质模型；
根据所述历史生产动态资料设置注采参数进行油藏数值模拟，得到所述地质资料和所述历史生产动态资料对应条件下的各生产井方向的日注入气体积；
根据所述地质资料、所述历史生产动态资料和对应条件下所述各生产井方向的日注入气体积，构成注入气体积样本库中的一组数据；
以所述注入气体积样本数据库中所述地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下各生产井方向的日注入气体积作为标签，对日注入气体积预测模型进行训练，得到训练好的日注入气体积预测模型。
4.根据权利要求1所述的气窜通道体积计算方法，其特征在于，在所述获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度之前，还包括：
获取所述目标区块的地质资料和生产动态资料；
根据所述目标区块的地质资料和生产动态资料，计算所述目标区块中各生产井与注气井之间的连通性系数。
5.根据权利要求1所述的气窜通道体积计算方法，其特征在于，所述生产动态资料包括目标区块中各口注气井的注气量、各口生产井的产量和井底流压。
6.根据权利要求1所述的气窜通道体积计算方法，其特征在于，所述地质资料包括目标区块的孔隙度、渗透率、油层厚度、地层压力、油气最小混相压力、油气粘度比和相对渗透率曲线。
7.根据权利要求1所述的气窜通道体积计算方法，其特征在于，所述目标区块中各生产井与注气井的连通性系数的数学模型如下：
Figure 39257DEST_PATH_IMAGE001
Figure 267982DEST_PATH_IMAGE002
Figure 375615DEST_PATH_IMAGE003
其中，
Figure 787136DEST_PATH_IMAGE004
为数学模型求解的目标函数，λ j 为第j口生产井与注气井之间的连通性系数；N为生产井的个数；q j (t)为公式计算得到的第j口生产井在t时刻的日产液量；
Figure 132667DEST_PATH_IMAGE005
为实际生产动态资料中第j口生产井在t时刻的日产液量；
Figure 17446DEST_PATH_IMAGE006
为第j口生产井的时间常数；t 0为初始生产时刻；t为当前生产时刻；i(ξ)为实际生产动态资料中注气井在ξ时刻对应的地下日注气量；v j 为第j口生产井对应的采油指数；
Figure 688792DEST_PATH_IMAGE007
为第j口生产井在t时刻的井底流压。
8.根据权利要求1所述的气窜通道体积计算方法，其特征在于，根据所述各生产井的示踪剂前缘推进速度和日注入气体积，计算得到各生产井的气窜通道体积，具体包括：
根据A j =Q j /v j 计算各生产井方向的气窜通道截面积；
其中，A j 为第j口生产井方向的气窜通道截面积；Q j 为第j口生产井方向的日注入气体积；v j 为第j口生产井方向的示踪剂前缘推进速度；
根据V j =A j L j 计算气窜通道体积；
其中，V j 为第j口生产井的气窜通道体积；L j 为第j口生产井对应的注采井距。
9.一种气窜通道体积计算系统，其特征在于，包括：
获取模块，用于获取目标区块中各生产井与注气井之间的连通性系数以及第一生产井的示踪剂前缘推进速度；所述第一生产井为所述目标区块中具有示踪剂资料的生产井；
示踪剂前缘推进速度计算模块，用于根据所述目标区块中各生产井与注气井之间的连通性系数以及所述第一生产井的示踪剂前缘推进速度，确定第二生产井在假设具有示踪剂资料的情况下的示踪剂前缘推进速度；所述第二生产井为所述目标区块中不具有示踪剂资料的生产井；
日注入气体积计算模块，用于将所述目标区块的地质资料和所述目标区块中各生产井与注气井的生产动态资料输入至训练好的日注入气体积预测模型，得到所述目标区块中各生产井方向的日注入气体积；所述日注入气体积预测模型为以所述区块的地质资料和所述历史生产动态资料为输入，以所述地质资料和所述历史生产动态资料对应条件下油藏数值模拟得到的各生产井方向的日注入气体积作为标签训练得到的模型；
气窜通道体积计算模块，用于根据所述目标区块中各生产井的示踪剂前缘推进速度和日注入气体积，计算得到所述目标区块中各生产井的气窜通道体积。
10.根据权利要求9所述的气窜通道体积计算系统，其特征在于，所述示踪剂前缘推进速度计算模块，具体包括：
拟合关系式确定单元，用于根据所述第一生产井的示踪剂前缘推进速度和各所述第一生产井与注气井之间的井间连通性系数，确定示踪剂前缘推进速度与对应井间连通性系数之间的拟合关系式；
示踪剂前缘推进速度计算单元，用于将所述第二生产井的井间连通性系数代入所述拟合关系式，得到所述第二生产井对应的示踪剂前缘推进速度。
CN202211495725.3A 2022-11-28 2022-11-28 一种气窜通道体积计算方法及系统 Active CN115600509B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211495725.3A CN115600509B (zh) 2022-11-28 2022-11-28 一种气窜通道体积计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211495725.3A CN115600509B (zh) 2022-11-28 2022-11-28 一种气窜通道体积计算方法及系统

Publications (2)

Publication Number Publication Date
CN115600509A true CN115600509A (zh) 2023-01-13
CN115600509B CN115600509B (zh) 2023-03-10

Family

ID=84852255

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211495725.3A Active CN115600509B (zh) 2022-11-28 2022-11-28 一种气窜通道体积计算方法及系统

Country Status (1)

Country Link
CN (1) CN115600509B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110863806A (zh) * 2019-11-28 2020-03-06 西安石油大学 一种二氧化碳驱气体前缘动态变化预测方法
CN111814347A (zh) * 2020-07-20 2020-10-23 中国石油大学(华东) 预测油藏内气窜通道的方法及系统
CN111852417A (zh) * 2019-04-30 2020-10-30 中国石油化工股份有限公司 二氧化碳驱油藏气窜通道快速反演方法和装置
CN112343587A (zh) * 2020-09-03 2021-02-09 中国石油天然气股份有限公司 一种特低渗透油藏优势渗流通道识别表征方法
CN113536653A (zh) * 2021-08-06 2021-10-22 中国石油大学(华东) 一种基于生产动态资料的气窜通道识别方法与系统
CN113743023A (zh) * 2021-10-12 2021-12-03 中国石油化工股份有限公司 一种二氧化碳驱油气窜通道分级表征方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111852417A (zh) * 2019-04-30 2020-10-30 中国石油化工股份有限公司 二氧化碳驱油藏气窜通道快速反演方法和装置
CN110863806A (zh) * 2019-11-28 2020-03-06 西安石油大学 一种二氧化碳驱气体前缘动态变化预测方法
CN111814347A (zh) * 2020-07-20 2020-10-23 中国石油大学(华东) 预测油藏内气窜通道的方法及系统
CN112343587A (zh) * 2020-09-03 2021-02-09 中国石油天然气股份有限公司 一种特低渗透油藏优势渗流通道识别表征方法
CN113536653A (zh) * 2021-08-06 2021-10-22 中国石油大学(华东) 一种基于生产动态资料的气窜通道识别方法与系统
CN113743023A (zh) * 2021-10-12 2021-12-03 中国石油化工股份有限公司 一种二氧化碳驱油气窜通道分级表征方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FENG LIU等: "Experimental Study of Oil Displacement and Gas Channeling during CO2 Flooding in Ultra—Low Permeability Oil Reservoir", 《ENERGIES》 *
张世明: "低渗透油藏CO2驱气窜通道识别方法", 《油气地质与采收率》 *
秦润森等: "秦皇岛32-6油田南区明下段Ⅰ油组3小层河道砂体叠置类型及其动态响应特征", 《油气地质与采收率》 *
陈德坡等: "利用井间动态连通性模型定量描述优势通道", 《大庆石油地质与开发》 *

Also Published As

Publication number Publication date
CN115600509B (zh) 2023-03-10

Similar Documents

Publication Publication Date Title
CN106837297B (zh) 一种识别井间连通性及油水动态预测的方法
CN104879103B (zh) 一种分层注水效果分析方法
US4273187A (en) Petroleum recovery chemical retention prediction technique
CN111581854A (zh) 一种考虑非平衡各向异性相对渗透率的油藏状态预测方法
Li et al. Modeling tracer flowback in tight oil reservoirs with complex fracture networks
Behmanesh Rate-transient analysis of tight gas condensate and black oil wells exhibiting two-phase flow
Da Silva et al. A review of coupled dynamic well-reservoir simulation
Yousefi et al. Interwell connectivity identification in immiscible gas-oil systems using statistical method and modified capacitance-resistance model: a comparative study
Nasriani et al. Investigation into the effect of capillary number on productivity of a lean gas condensate reservoir
Ganesh et al. Reduced physics modeling of CO2 injectivity
Dietrich Relative permeability during cyclic steam stimulation of heavy-oil reservoirs
CN106014365A (zh) 一种预测水驱开发油田产量递减率的方法
CN115600509B (zh) 一种气窜通道体积计算方法及系统
Begum et al. A case study of reservoir parameter estimation in Norne oil field, Norway by using Ensemble Kalman Filter (EnKF)
Kim Development of linear capacitance-resistance models for characterizing waterflooded reservoirs
Zakirov et al. Well test for in-situ determination of oil and water relative permeabilities
Du et al. New Gas Tracer Convection–Diffusion Model between Wells in Heavy Oil Reservoirs
Chen et al. The application of Stefan problem in calculating the lateral movement of steam chamber in SAGD
Zakirov et al. In-situ determination of displacement efficiency and oil and water relative permeability curves through integrated well test study at exploration-to-pilot stage of the oilfield development project
CN111364955A (zh) 注采井间流场演化模拟方法
Dong et al. A semi-analytical method for optimizing the gas and water bidirectional displacement in the tilted fault block reservoir
Mahmud et al. Investigating relative permeability measurements using unsteady-state core flooding method
CN110778313A (zh) 一种稀土元素示踪剂井间监测技术
CN113605883B (zh) 一种剩余油饱和度解析的方法
Yu et al. Rate transient analysis methods for water-producing gas wells in tight reservoirs with mobile water

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