CN112160734B - 注采井相关性分析方法、装置、存储介质及计算机设备 - Google Patents

注采井相关性分析方法、装置、存储介质及计算机设备 Download PDF

Info

Publication number
CN112160734B
CN112160734B CN202011102067.8A CN202011102067A CN112160734B CN 112160734 B CN112160734 B CN 112160734B CN 202011102067 A CN202011102067 A CN 202011102067A CN 112160734 B CN112160734 B CN 112160734B
Authority
CN
China
Prior art keywords
well
production
injection
water injection
correlation
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
Application number
CN202011102067.8A
Other languages
English (en)
Other versions
CN112160734A (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 Beijing
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Original Assignee
China University of Petroleum Beijing
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
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 Beijing, China National Offshore Oil Corp CNOOC, CNOOC China Ltd Tianjin Branch filed Critical China University of Petroleum Beijing
Priority to CN202011102067.8A priority Critical patent/CN112160734B/zh
Publication of CN112160734A publication Critical patent/CN112160734A/zh
Application granted granted Critical
Publication of CN112160734B publication Critical patent/CN112160734B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/20—Displacing by water
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH DRILLING; MINING
    • E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00—Survey of boreholes or wells
    • E21B47/06—Measuring temperature or pressure
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Mining
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Business, Economics & Management (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Marketing (AREA)
  • Agronomy & Crop Science (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Animal Husbandry (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

本发明涉及一种注采井相关性分析方法、装置、存储介质及计算机设备，其中，注采井相关性分析方法包括以下步骤：1)获取不同时间注入井的注水量、采出井的产液量和压力数据，其中，所述压力数据包括：注入井压力、采出井压力和地层压力；2)并根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；3)计算调整后的注水量和产液量；4)建立调整后的注水量和调整后的产液量之间的关系；5)根据调整后的注水量和调整后的产液量之间的关系，分析注采井相关性。本发明能够准确估计注采井相关性，提高油田生产效率，减少无效作业的成本。

Description

注采井相关性分析方法、装置、存储介质及计算机设备
技术领域
本发明涉及油气开采的技术领域，具体是关于一种注采井相关性分析方法、装置、存储介质及计算机设备。
背景技术
水驱是目前世界上应用最广、技术经济效益最好的油田开发方式，约有60％～70％的油田主要是依靠水驱开发的，水驱年采油量占世界原油产量的60％以上，水驱管理是水驱开发油田稳产或减缓产量递减的重要途径之一，水驱管理的关键除了层系划分外，主要是井网部署、射孔优化及注采参数优化调整。在开发过程中，需要根据注采井相关性进行注入井的注水量和采出井的产液量进行调配，注入井进行调剖，采油井堵水、补孔、打加密井来完善注入井和采出井之间的相关性(简称注采相关性)。因此，注采井相关性研究是水驱管理的核心内容，注采井相关性认识是否清楚直接影响水驱管理效率和效果。
小层对比、沉积相等油田地质研究及注入井和采出井动态分析是目前注采关系的主要研究方法，经常采用油田地质研究与开发动态分析相结合的方法。一般是由专家综合地质信息(例如地层、井位、断层等信息)和产量实现的，这种方法具有如下不足：
1)复杂的油藏环境和油田海量数据使得专家进行人工分析方法非常困难且效率低；
2)采出井附近注入井注采关系分析的局限性常常无法估计出因裂缝等强非均质性造成的远距离注采关系。
根据注水井的注入量与采油井的产液量来确定注水井与采油井之间的相关性是推断注采井相关性的一个重要方法。有许多学者基于注入井的注水量和采出井的产液量来估计注采相关性，把油藏视为一个由连续脉冲响应所表征的系统，将油田生产近似看成由输入信号(即注水量)产生输出信号(即产液量)的过程。可将油藏模型建成一个由权值代表的电阻模型、由表示井间连通性和耗散性参数表征的电容模型或电容-电阻模型，井间连通性可由模型参数多线性回归量化估计。但是，考虑到现场井底压力变化、修井和天然或人工所导致的地质影响，实际注水量和产液量的变化较大，含有的噪声较多，导致注采井相关性估计结果浮动较大，准确度较低。
发明内容
针对上述问题，本发明的目的是提供一种注采井相关性分析方法、装置、存储介质及计算机设备，能够准确估计注采井相关性，提高油田生产效率，减少无效作业的成本。
为实现上述目的，本发明采取以下技术方案：
本发明所述的注采井相关性分析方法，包括以下步骤：
1)获取不同时间注入井的注水量、采出井的产液量和压力数据，其中，所述压力数据包括：注入井压力、采出井压力和地层压力；
2)并根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；
3)计算调整后的注水量和产液量；
4)建立调整后的注水量和调整后的产液量之间的关系；
5)根据调整后的注水量和调整后的产液量之间的关系，分析注采井相关性。
所述的注采井相关性分析方法，优选地，所述步骤3)包括如下子步骤：
3.1)注采井组的选取
以注水井W为中心，选定两个井距范围内的采油井O，形成注采井组合井对(W,O)；如果井距大于两个井距或者注采井之间存在有效封闭断层，则注采井相关性IPR(W,O)的值为零；
3.2)动态资料选取
针对所选注采井组中IPR(W,O)值不为零的注采井，选取射孔数据、生产动态数据、注水井管柱设计参数、生产井压力测试数据和油藏平均压力数据；并应用注水井管柱设计参数按照式(1-6)计算井筒压力损失；
Pf＝Pf1+Pf2+Pf3 (1)
式中，Pf为井筒压力损失；Pf1为垂直管损；Pf2为水嘴损失；Pf3为炮眼摩阻；
其中，垂直管损是从井口到井底炮眼处的沿程压力损失，根据水力学公式(2-4)计算垂直管损：
Figure BDA0002725710670000021
Figure BDA0002725710670000022
Figure BDA0002725710670000023
式中，λ为沿程阻力系数，无因次；L为井口到井底炮眼处的距离；D为油管直径；r0为井筒半径；v为油管中注入水平均流速；g为重力加速度；μ为流体粘度；ρ为流体密度；Re为雷诺数，用来表征流体流动情况的无量纲数；
其中，水嘴损失利用实验得到的经验公式(5)求得或参照嘴损图表查得：
Figure BDA0002725710670000031
式中，Q为单井日注水量；Dw为水嘴直径；
其中，炮眼摩阻通过公式(6)求得，具体为：
Figure BDA0002725710670000032
式中，ρ为流体密度；d为炮眼直径；N为炮眼个数；q为注水量；
3.3)根据式(7)计算调整后的注水量：
Figure BDA0002725710670000033
3.4)根据式(8)计算调整后的产液量：
Figure BDA0002725710670000034
所述的注采井相关性分析方法，优选地，所述步骤4)包括如下子步骤：
4.1)预估脉冲响应参数的初始值；
4.2)根据所述调整后的注水量和产液量，以及脉冲响应参数的初始值，建立状态方程，并根据状态方程推出辅助函数Ki；
4.3)利用递归算法获得所述状态方程中各脉冲响应参数的最终值，并计算所选时间段内所有注采井相关性估值和注采井相关性估值平均值。
所述的注采井相关性分析方法，优选地，所述步骤4.1)包括如下子步骤：
4.1.1)根据条件a1∈[1，10]和β1∈[0，1]下选取参数α1和脉冲响应参数β1；
4.1.2)应用参数α1根据下式计算脉冲响应参数α1：
Figure BDA0002725710670000035
其中α1为脉冲响应参数，a1为参数，T为时间间隔；
4.1.3)根据条件γ1∈(0，1)选取脉冲响应参数γ1。
所述的注采井相关性分析方法，优选地，所述步骤4.2)包括如下子步骤：
4.2.1)应用各脉冲响应参数的αi，βi和γi，以及调整后的注水量xi和产液量yi建立状态方程Fi为：
Figure BDA0002725710670000041
其中，F为状态方程，αi、βi和γi分别为脉冲响应参数，xi为调整后的注水量，yi为调整后的产液量，i＝1,2,3...；ni为噪声参数；
4.2.2)由状态方程Fi推出辅助函数Ki为：
Figure BDA0002725710670000042
式中，Ki为辅助函数，Fi为状态方程，yi为调整后的产液量，Var(yi)为yi值的方差。
所述的注采井相关性分析方法，优选地，所述步骤4.3)包括如下子步骤：
4.3.1)应用递归算法计算调整后注水量变化引起的调整后产液量的变化量；
4.3.2)获取各脉冲响应参数的最终值α0，β0和γ0，其中，所述脉冲响应参数的最终值为调整后产液量响应变化量与实际产液量变化量相同时的脉冲响应参数；
4.3.3)应用各脉冲响应参数的最终值α0，β0和γ0通过式(12)计算注采井间相关程度估值IPR(W，O)，
Figure BDA0002725710670000043
式中，IPR(W,O)为注采井相关性估值；α0，β0和γ0为脉冲响应参数的最终值；
4.3.4)求取所选时间段的所有IPR(W，O)，计算注采井相关性估值平均值
Figure BDA0002725710670000044
所述的注采井相关性分析方法，优选地，所述步骤5)包括如下子步骤：
5.1)利用归一化公式以注水井为中心的n个相关性估值平均值
Figure BDA0002725710670000045
值进行100以内的归一化处理，统一量化注水井与周围各采油井的相关性，所述归一化公式为：
Figure BDA0002725710670000046
5.2)将归一化后的评分结果排序，通过排序结果显示出注水井与周围采油井的相关性强弱，用于划分注水井与周边采油井的相关性级别，判定与注水井相关性较大的采油井；
5.3)结合油藏地质资料对评分结果进行验证。
本发明所述的注采井相关性分析装置，包括：
数据获取模块，用于获取注入井的注水量，采出井的产液量和压力数据，所述压力数据包括：注入井压力，采出井压力和地层压力；
数据调整模块，用于根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；
分析模块，用于建立调整后的注水量和产液量之间的关系，根据调整后的注水量和产液量之间的关系，分析注采井相关性。
本发明所述的计算机可读存储介质，其中，所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明所述的计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其中，所述处理器执行所述计算机程序时实现上述方法。
本发明由于采取以上技术方案，其具有以下优点：
1)应用调整后的注入井注水量和采出井产液量来估计油藏的注采相关性，而不是简单分析注水量与产液量之间的动态响应程度。综合考虑油水井射孔情况、井底压力分布情况等实际因素，与油田实际情况匹配程度更高。
2)循环使用递归算法时引入噪声函数，减小了因井距和渗透率等因素造成的估计误差，很大程度提高了估计准确性。
3)可以将一个井网单元中的多口采出井单独考虑来估计注采关系，也可将多口采出井作为一个采出井组来估计注采关系。
4)不受注采井数量限制，尤其适用于大规模井网部署下的多口注入井和多口采出井相关性估计。
5)利用计算机存储和调用数据、计算调整后注水量变化产生的调整后产液量变化响应量、分析注采井相关性、绘制注采井相关性示意图，与费时费力的传统地质研究方法相比，大幅度提高了效率。
附图说明
图1为本发明所述的注采井相关性分析方法示意图；
图2为本发明所述的注采井相关性分析方法的流程图；
图3为本发明主力开发层系井位图；
图4为本发明注水井调整后的注水量和产液量曲线图；
图5为本发明注采井相关性估值曲线图；
图6为本发明注采井相关性示意图。
具体实施方式
以下将结合附图对本发明的较佳实施例进行详细说明，以便更清楚理解本发明的目的、特点和优点。应理解的是，附图所示的实施例并不是对本发明范围的限制，而只是为了说明本发明技术方案的实质精神。
如图1和图2所示，本发明提供的一种注采井相关性分析方法，包括以下步骤：
步骤1：获取不同时间注入井的注水量、采出井的产液量和压力数据，其中压力数据包括：注入井压力、采出井压力和地层压力；
步骤2：根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；
步骤3：计算调整后的注水量和产液量，具体包括如下内容：
3.1：注采井组的选取：
以注水井W为中心，选定两个井距范围内的采油井O，形成注采井组合井对(W,O)；如果井距大于两个井距或者注采井之间存在有效封闭断层，则注采井相关性IPR(W,O)的值为零；
3.2：动态资料选取：
针对所选注采井组中IPR(W,O)值不为零的注采井，选取射孔数据、生产动态数据、注水井管柱设计参数、生产井压力测试数据和油藏平均压力等数据，并应用注水井管柱设计参数按照式(1)计算井筒压力损失：
Pf＝Pf1+Pf2+Pf3 (1)
式中，Pf为井筒压力损失；Pf1为垂直管损；Pf2为水嘴损失；Pf3为炮眼摩阻；
其中，垂直管损是从井口到井底炮眼处的沿程压力损失，根据水力学公式(2-4)计算垂直管损：
Figure BDA0002725710670000061
Figure BDA0002725710670000062
Figure BDA0002725710670000063
式中，λ为沿程阻力系数，无因次；L为井口到井底炮眼处的距离；D为油管直径；r0为井筒半径；v为油管中注入水平均流速；g为重力加速度；μ为流体粘度；ρ为流体密度；Re为雷诺数，用来表征流体流动情况的无量纲数；
水嘴损失可用实验得到的经验公式(5)求得或参照嘴损图表1查得：
Figure BDA0002725710670000071
式中，Q为单井日注水量；Dw为水嘴直径。
表1不同排量下的水嘴压力损失
排量/m<sup>3</sup>/d 压力损失/MPa
≤200 0.1
200～300 0.3
300～400 0.5
500～600 1.0
炮眼摩阻可通过公式(6)求得，具体为：
Figure BDA0002725710670000072
式中，ρ为流体密度；d为炮眼直径；N为炮眼个数；q为注水量。
3.3：根据式(7)计算调整后的注水量：
Figure BDA0002725710670000073
3.4：根据式(8)计算调整后的产液量：
Figure BDA0002725710670000074
步骤4：建立调整后的注水量和产液量之间的关系，具体包括如下内容：
4.1：预估各脉冲响应参数的初始值：
4.1.1：根据条件a1∈[1，10]和β1∈[0，1]下选取参数a1和脉冲响应参数β1，若地质条件好，则a1和β1的数值可选取大些，若地质条件较差，则a1和β1的数值选取小些，若地质条件未知，则a1和β1的数值可选取适中。
4.1.2：应用参数a1根据下式计算脉冲响应参数α1：
Figure BDA0002725710670000075
其中α1为脉冲响应参数，a1为参数，T为时间间隔；
4.1.3：根据条件γ1∈(0，1)选取脉冲响应参数γ1。
4.2：根据所述调整后的注水量和产液量，以及脉冲响应参数的初始值，建立状态方程，并根据状态方程推出辅助函数Ki；
4.2.1：应用各脉冲响应参数的αi，βi和γi，以及调整后的注水量xi和产液量yi建立状态方程Fi为：
Figure BDA0002725710670000081
式中，αi，βi和γi分别为脉冲响应参数，xi为调整后的注水量，yi为调整后的产液量，i＝1,2,3...；ni为噪声参数。
4.2.2：由状态方程Fi推出辅助函数Ki为：
Figure BDA0002725710670000082
式中，Ki为辅助函数，Fi为状态方程，yi为调整后的产液量，Var(yi)为yi值的方差。
4.3：利用递归算法获得所述状态方程中各脉冲响应参数的最终值，并计算所选时间段内所有注采井相关性估值和注采井相关性估值平均值：
4.3.1：应用递归算法计算调整后注水量变化引起的调整后产液量的变化量；
4.3.2：获取各脉冲响应参数的最终值α0，β0和γ0，其中脉冲响应参数的最终值为调整后产液量响应变化量与实际产液量变化量相同时的脉冲响应参数；
4.3.3：应用各脉冲响应参数的最终值α0，β0和γ0通过式(12)计算注采井间相关程度估值IPR(W，O)：
Figure BDA0002725710670000083
式中，IPR(W,O)为注采井相关性估值；α0，β0和γ0为脉冲响应参数的最终值；
4.3.4：求取所选时间段的所有IPR(W，O)，计算注采井相关性估值平均值
Figure BDA0002725710670000084
需要说明的是：循环使用递归算法时引入辅助函数，减小了因井距和渗透率等因素造成的估计误差，很大程度提高了估计准确性。
步骤5：根据调整后的注水量和产液量之间的关系，分析注采井相关性，具体包括如下内容：
5.1：利用归一化公式以注水井为中心的n个相关性估值平均值
Figure BDA0002725710670000091
值进行100以内的归一化处理，统一量化注水井与周围各采油井的相关性，所述归一化公式为：
Figure BDA0002725710670000092
5.2：将归一化后的评分结果排序，通过排序结果显示出注水井与周围采油井的相关性强弱，用于划分注水井与周边采油井的相关性级别，判定与注水井相关性较大的采油井；
5.3：结合油藏地质资料对评分结果进行验证。
本发明还提供了一种注采井相关性分析装置，该装置包括：
数据获取模块，用于获取注入井的注水量，采出井的产液量和压力数据，所述压力数据包括：注入井压力，采出井压力和地层压力；
数据调整模块，用于根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；
分析模块，用于建立调整后的注水量和产液量之间的关系，根据调整后的注水量和产液量之间的关系，分析注采井相关性。
本发明还提供了一种计算机可读存储介质，其中，计算机可读存储介质存储有执行上述方法的计算机程序。
本发明还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其中，处理器执行计算机程序时实现上述方法。
实施例1
下面给出一个具体实施例，说明本发明实施例中注采井相关性分析方法的具体应用。某一实际油田为例，选取某一区域进行研究，处理流程如图2所示。该区域主力开发层系井位图如图3所示。设以注水井Inj1为中心注水井，选定两个井距范围内的采油井，形成注采井组合对，共选取9个注采井组合对，分别是(Inj1，P2)，(Inj1，P3)，(Inj1，P4)，(Inj1，P5)，(Inj1，P6)，(Inj1，P7)，(Inj1，P8)，(Inj1，P9)，(Inj1，P10)，各注采井组合对的注采井相关性估值IPR(W,O)待定。井距大于两个井距的注采对，包括(Inj1，P1)和(Inj1，P11)的IPR(W,O)＝0。选取注采井组之后，针对所选注采井组的注采井，获取注入井的注水量，采出井的产液量和压力数据，所述压力数据包括：注入井压力，采出井压力和地层压力。根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量。其中，调整后的注水量和产液量结果如图4所示。由于该油藏地质条件未知，为确保计算稳定性，选取脉冲响应参数初值α1＝3，β1＝0.2和γ1＝0.2，脉冲响应参数在循环运算过程中不断收敛，最终进入最佳范围后替代原始初值重新计算，将初值误差降至最低。
应用脉冲响应参数初值、调整后的注水量和产液量建立状态方程F与辅助函数K。状态方程F与辅助函数K相结合，应用递归算法计算调整后的注水量变化引起的调整后的产液量响应变化量。将调整后的产液量的响应变化量与实际产液量变化量相同时的脉冲响应参数作为脉冲响应参数的最终值，表示为α0，β0和γ0。应用脉冲响应参数计算注采井组合(Inj1，P4)的注采井相关性估值IPR(W,O)，结果如图5所示。以注水井Inj1井为中心选取的9个注采井组合，分别应用上述方法计算得到的
Figure BDA0002725710670000101
值。以注水井为中心的所有
Figure BDA0002725710670000102
值进行100以内的归一化处理，将归一化后的评分结果排序，计算结果如表2所示。由上述分析结果可绘制注采井相关性示意图，如图6所示。
表2归一化后的评分结果排序
Figure BDA0002725710670000103
结合注采井相关性估计结果可知，与注水井Inj1相关性强的采油井为P8、P9、P10。与油田现场已有的注采井相关性认知一致，具有较高的准确性。
综上所述，本发明实施例通过获取注入井的注水量，采出井的产液量和压力数据，所述压力数据包括：注入井压力，采出井压力和地层压力，根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量，建立调整后的注水量和产液量之间的关系，根据调整后的注水量和产液量之间的关系，分析注采井相关性。本发明实施例考虑了压力变化对注入井的注水量和采出井的产液量的影响，利用压力数据对注水量和产液量进行调整，从而抵消了由压力变化引起的误差，然后根据调整后的注水量和产液量之间的关系，分析注采井相关性，获得稳定、准确度高的相关性分析结果，为提高注入井调剖、采出井堵水及酸化压裂等油层改造措施效果，提高油田生产效率，减少无效作业成本奠定基础。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (9)

1.一种注采井相关性分析方法，其特征在于，包括以下步骤：
1)获取不同时间注入井的注水量、采出井的产液量和压力数据，其中，所述压力数据包括：注入井压力、采出井压力和地层压力；
2)并根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；
3)计算调整后的注水量和产液量；
4)建立调整后的注水量和调整后的产液量之间的关系；
5)根据调整后的注水量和调整后的产液量之间的关系，分析注采井相关性；
所述步骤3)包括如下子步骤：
3.1)注采井组的选取
以注水井W为中心，选定两个井距范围内的采油井O，形成注采井组合井对(W,O)；如果井距大于两个井距或者注采井之间存在有效封闭断层，则注采井相关性IPR(W,O)的值为零；
3.2)动态资料选取
针对所选注采井组中IPR(W,O)值不为零的注采井，选取射孔数据、生产动态数据、注水井管柱设计参数、生产井压力测试数据和油藏平均压力数据；并应用注水井管柱设计参数按照式(1)-(6)计算井筒压力损失；
Pf＝Pf1+Pf2+Pf3 (1)
式中，Pf为井筒压力损失；Pf1为垂直管损；Pf2为水嘴损失；Pf3为炮眼摩阻；
其中，垂直管损是从井口到井底炮眼处的沿程压力损失，根据水力学公式(2)-(4)计算垂直管损：
Figure FDA0003902455230000011
Figure FDA0003902455230000012
Figure FDA0003902455230000013
式中，λ为沿程阻力系数，无因次；L为井口到井底炮眼处的距离；D为油管直径；v为油管中注入水平均流速；g为重力加速度；μ为流体粘度；ρ为流体密度；Re为雷诺数，用来表征流体流动情况的无量纲数；
其中，水嘴损失利用实验得到的经验公式(5)求得或参照嘴损图表查得：
Figure FDA0003902455230000021
式中，Q为单井日注水量；Dw为水嘴直径；
其中，炮眼摩阻通过公式(6)求得，具体为：
Figure FDA0003902455230000022
式中，ρ为流体密度；d为炮眼直径；N为炮眼个数；q为注水量；
3.3)根据式(7)计算调整后的注水量：
Figure FDA0003902455230000023
3.4)根据式(8)计算调整后的产液量：
Figure FDA0003902455230000024
2.根据权利要求1所述的注采井相关性分析方法，其特征在于，所述步骤4)包括如下子步骤：
4.1)预估脉冲响应参数的初始值；
4.2)根据所述调整后的注水量和产液量，以及脉冲响应参数的初始值，建立状态方程，并根据状态方程推出辅助函数Ki；
4.3)利用递归算法获得所述状态方程中各脉冲响应参数的最终值，并计算所选时间段内所有注采井相关性估值和注采井相关性估值平均值。
3.根据权利要求2所述的注采井相关性分析方法，其特征在于，所述步骤4.1)包括如下子步骤：
4.1.1)根据条件a1∈[1，10]和β1∈[0，1]下选取参数a1和脉冲响应参数初始值β1；
4.1.2)应用参数a1根据下式计算脉冲响应参数初始值a1：
Figure FDA0003902455230000025
其中α1为脉冲响应参数初始值，a1为参数，T为时间间隔；
4.1.3)根据条件γ1∈(0，1)选取脉冲响应参数初始值γ1。
4.根据权利要求3所述的注采井相关性分析方法，其特征在于，所述步骤4.2)包括如下子步骤：
4.2.1)应用各脉冲响应参数的αi，βi和γi，以及调整后的注水量xi和产液量yi建立状态方程Fi为：
Figure FDA0003902455230000031
其中，F为状态方程，αi、βi和γi分别为脉冲响应参数，xi为调整后的注水量，yi为调整后的产液量，i＝1,2,3...；ni为噪声参数；
4.2.2)由状态方程Fi推出辅助函数Ki为：
Figure FDA0003902455230000032
式中，Ki为辅助函数，Fi为状态方程，yi为调整后的产液量，Var(yi)为yi值的方差。
5.根据权利要求4所述的注采井相关性分析方法，其特征在于，所述步骤4.3)包括如下子步骤：
4.3.1)应用递归算法计算调整后注水量变化引起的调整后产液量的变化量；
4.3.2)获取各脉冲响应参数的最终值α0，β0和γ0，其中，所述脉冲响应参数的最终值为调整后产液量响应变化量与实际产液量变化量相同时的脉冲响应参数；
4.3.3)应用各脉冲响应参数的最终值α0，β0和γ0通过式(12)计算注采井相关性IPR(W，O)，
Figure FDA0003902455230000033
式中，IPR(W,O)为注采井相关性估值；α0，β0和γ0为脉冲响应参数的最终值；
4.3.4)求取所选时间段的所有IPR(W，O)，计算注采井相关性平均值
Figure FDA0003902455230000034
6.根据权利要求5所述的注采井相关性分析方法，其特征在于，所述步骤5)包括如下子步骤：
5.1)利用归一化公式对以注水井为中心的n个相关性估值平均值
Figure FDA0003902455230000035
值进行100以内的归一化处理，统一量化注水井与周围各采油井的相关性，所述归一化公式为：
Figure FDA0003902455230000041
5.2)将归一化后的评分结果排序，通过排序结果显示出注水井与周围采油井的相关性强弱，用于划分注水井与周边采油井的相关性级别，判定与注水井相关性较大的采油井；
5.3)结合油藏地质资料对评分结果进行验证。
7.一种基于权利要求1-6任一项所述的注采井相关性分析方法的注采井相关性分析装置，其特征在于，包括：
数据获取模块，用于获取注入井的注水量，采出井的产液量和压力数据，所述压力数据包括：注入井压力，采出井压力和地层压力；
数据调整模块，用于根据所述注入井的注水量，注入井压力和地层压力，调整注水量；根据所述采出井的产液量，采出井压力和地层压力，调整产液量；
分析模块，用于建立调整后的注水量和产液量之间的关系，根据调整后的注水量和产液量之间的关系，分析注采井相关性。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至6任一所述方法的计算机程序。
9.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
CN202011102067.8A 2020-10-15 2020-10-15 注采井相关性分析方法、装置、存储介质及计算机设备 Active CN112160734B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011102067.8A CN112160734B (zh) 2020-10-15 2020-10-15 注采井相关性分析方法、装置、存储介质及计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011102067.8A CN112160734B (zh) 2020-10-15 2020-10-15 注采井相关性分析方法、装置、存储介质及计算机设备

Publications (2)

Publication Number Publication Date
CN112160734A CN112160734A (zh) 2021-01-01
CN112160734B true CN112160734B (zh) 2023-02-07

Family

ID=73867036

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011102067.8A Active CN112160734B (zh) 2020-10-15 2020-10-15 注采井相关性分析方法、装置、存储介质及计算机设备

Country Status (1)

Country Link
CN (1) CN112160734B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113807677A (zh) * 2021-09-01 2021-12-17 中控智网(北京)能源技术有限公司 油田能耗指标确定方法、装置、设备及存储介质
CN114183106A (zh) * 2021-11-29 2022-03-15 华鼎鸿基采油技术服务（北京）有限公司 一种化学驱替液的注入控制方法和系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2873722C (en) * 2012-05-14 2017-03-21 Landmark Graphics Corporation Method and system of predicting future hydrocarbon production
CN103939066B (zh) * 2014-04-04 2016-04-27 中国石油大学(华东) 一种一注多采井组定注水量确定油井产液量的方法
CN104747147B (zh) * 2015-03-24 2017-05-10 中国石油天然气股份有限公司 一种水驱开发油藏井网注采关系确定方法

Also Published As

Publication number Publication date
CN112160734A (zh) 2021-01-01

Similar Documents

Publication Publication Date Title
CN104879103B (zh) 一种分层注水效果分析方法
CN107044277B (zh) 低渗透非均质油藏水平井重复压裂增产潜力评价方法
CN112160734B (zh) 注采井相关性分析方法、装置、存储介质及计算机设备
CN104504604B (zh) 一种定性气井井筒积液的方法
WO2012051196A2 (en) Lift-gas optimization with choke control
CN112343587A (zh) 一种特低渗透油藏优势渗流通道识别表征方法
WO2019242296A1 (zh) 注水与采油井关系模型获得方法及产量和注水量确定方法
CN103912248A (zh) 水驱油田预测含水率方法
CN104453834A (zh) 一种井组注采关系优化调整方法
CN106150477A (zh) 一种确定单井控制储量的方法
CN108518219B (zh) 一种基于天然能量开发油田规模提液选井的评价方法
CN113792479B (zh) 一种基于物理约束的煤层气藏压裂效果评价方法
CN111810101B (zh) 一种水驱油藏动态分析方法及装置
CN105386751A (zh) 一种基于油藏渗流模型的水平井测井产能预测方法
CN106250984A (zh) 油井的油水关系模式的判断方法及装置
CN107120111A (zh) 一种基于多重分形的油藏井间连通程度的评估方法及系统
CN103075142A (zh) 一种注水开发油田堵水油井选井方法
CN108150160A (zh) 一种地层中欠压实作用超压的求取方法
CN109063403B (zh) 一种滑溜水压裂优化设计方法
CN106875286A (zh) 一种聚合物驱油田全过程注聚参数分阶段优化决策方法
CN107895092A (zh) 一种基于复杂非线性注采建模的井间连通定量评价方法
CN107035348B (zh) 一种油田调剖用多因素分析选井方法及装置
CN109441415A (zh) 基于邻井干扰的聚合物驱油藏测试井的试井解释方法
CN110390154A (zh) 一种提高复杂断块油气田油藏数值模拟效率的方法
CN106481315A (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