CN114139386B - 一种基于多层合采的连通性评估方法及装置 - Google Patents

一种基于多层合采的连通性评估方法及装置 Download PDF

Info

Publication number
CN114139386B
CN114139386B CN202111465675.XA CN202111465675A CN114139386B CN 114139386 B CN114139386 B CN 114139386B CN 202111465675 A CN202111465675 A CN 202111465675A CN 114139386 B CN114139386 B CN 114139386B
Authority
CN
China
Prior art keywords
production
connectivity
layer
index
coefficient
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
CN202111465675.XA
Other languages
English (en)
Other versions
CN114139386A (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
Original Assignee
China University of Petroleum Beijing
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 filed Critical China University of Petroleum Beijing
Priority to CN202111465675.XA priority Critical patent/CN114139386B/zh
Publication of CN114139386A publication Critical patent/CN114139386A/zh
Application granted granted Critical
Publication of CN114139386B publication Critical patent/CN114139386B/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
    • 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)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本说明书实施例提供一种基于多层合采的连通性评估方法及装置。所述方法包括：获取对应于目标工区至少两个开发地层的地质参数和生产参数；确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；结合所述单层参数影响关系和合采参数影响关系计算连通性系数；基于所述连通性系数评估所述目标工区的连通性。上述方法保证了连通性系数的准确性以及与实际生产工况之间的关联性，提高了生产方案的开采效果，保证了生产开发的有效进行。

Description

一种基于多层合采的连通性评估方法及装置
技术领域
本说明书实施例涉及地质勘探开发技术领域，特别涉及一种基于多层合采的连通性评估方法及装置。
背景技术
在油气田开采过程中，可能会出现工区下方存在多个储层的情况。为了降低开采成本，提高储层动用程度，一般会采用合层开采的方式进行开发，即同时针对多个储层进行开采，以改善钻井效果。
在油藏开发的后期，单纯依靠天然能量进行开发已经不能满足实际的开发需求量，因此，往往通过注水、注气以及生物采油等方式向地层中补充能量。而在实际应用中，由于各层段储层性质存在差异，流通性不同，流体渗流状态也不相同，导致向储层中注入液体时，不同储层的流体推进速度、见水时间都存在差异，流体渗流状态复杂。
因此，多层非均质油藏在开发生产过程中层间矛盾突出，层间干扰现象也比较严重。而确定不同储层之间的油藏连通性对于查找开发后期油井低产的原因具有重要影响，并能够相应地制定合理的开发方案，进而能够有效降低层间干扰以提高油气采收率。但是，目前在确定油藏连通性时，难以同时保证计算结果的准确性和计算量大小，使得目前的计算方法无法在实际生产中的得到有效应用。因此，目前亟需一种能够快速准确地对工区中各个储层之间的连通性进行评估以制定开发方案的方法。
发明内容
本说明书实施例的目的是提供一种基于多层合采的连通性评估方法及装置，以解决如何快速准确地对工区中各个储层之间的连通性进行评估以制定开发方案的问题。
为了解决上述技术问题，本说明书实施例提出一种基于多层合采的连通性评估方法，所述方法包括：获取对应于目标工区至少两个开发地层的地质参数和生产参数；确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；基于所述连通性系数评估所述目标工区的连通性。
本说明书实施例还提出一种基于多层合采的连通性评估装置，所述装置包括：参数获取模块，用于获取对应于目标工区至少两个开发地层的地质参数和生产参数；单层参数影响关系确定模块，用于确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；合采参数影响关系确定模块，用于确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；连通性系数计算模块，用于结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；连通性评估模块，用于基于所述连通性系数评估所述目标工区的连通性。
本说明书实施例还提出一种基于多层合采的连通性评估设备，包括存储器和处理器；所述存储器，用于存储计算机程序指令；所述处理器，用于执行所述计算机程序指令以实现以下步骤：获取对应于目标工区至少两个开发地层的地质参数和生产参数；确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；基于所述连通性系数评估所述目标工区的连通性。
由以上本说明书实施例提供的技术方案可见，本说明书实施例在获取到开发地层的地质参数和生产参数后，首先通过实验模拟结果确定地质参数与生产参数之间的单层参数影响关系和合采参数影响关系，并进一步结合所述单层参数影响关系和合采参数影响关系计算连通性系数，从而有效对不同开发地层之间的连通性进行评估。通过上述方法，在考虑单采、合采以及实验模拟、实际开发的不同情况下，生产效果与地质的对应关系后，确定不同开发地层之间的连通性系数，保证了连通性系数的准确性以及与实际生产工况之间的关联性，提高了生产方案的开采效果，保证了生产开发的有效进行。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本说明书实施例一种基于多层合采的连通性评估方法的流程图；
图2A为本说明书实施例一种单层参数影响关系变化曲线的示意图；
图2B为本说明书实施例一种单层参数影响关系变化曲线的示意图；
图3A为本说明书实施例一种合采参数影响关系变化曲线的示意图；
图3B为本说明书实施例一种合采参数影响关系变化曲线的示意图；
图4为本说明书实施例一种连通性系数变化规律的示意图；
图5是本说明书实施例一种基于多层合采的连通性评估装置的模块图；
图6是本说明书实施例一种基于多层合采的连通性评估设备的结构图。
具体实施方式
下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书保护的范围。
为了解决上述技术问题，本说明书实施例提出了一种基于多层合采的连通性评估方法。所述基于多层合采的连通性评估方法的执行主体为基于多层合采的连通性评估设备，所述基于多层合采的连通性评估设备包括但不限于服务器、工控机、PC机等。如图1所示，所述基于多层合采的连通性评估方法具体可以包括以下步骤。
S110：获取对应于目标工区至少两个开发地层的地质参数和生产参数。
目标工区是需要进行生产作业的工区，例如可以是油气田所对应的区域。所述目标工区可以是还未进行开发的工区，需要制定新的开发方案，也可以是处于开发中的工区，需要根据工区的实际情况调整当前的施工方案，对此不做限制。
在本实施例中，目标工区对应有至少两个开发地层，开发地层即为目标工区内的储层，即在所述目标工区的不同深度存在有多个开发地层。为了提高开发效率，可以针对目标工区内的这些开发地层进行多层合采。
地质参数可以是针对目标工区进行探测后，得到的相应地质条件所对应的参数，例如针对不同的开发地层，地质参数可以是对应的渗透率、含水率及地层厚度等。实际应用中根据计算需要可以获取其他参数作为所述地质参数，对此不做限制。
生产参数是对应于目标工区生产过程中各个开发地层的生产效果所体现出的参数值。所述生产参数可以是通过实验的方式模拟目标工区的开采过程，所采集得到的参数，也可以是目标工区在开发的过程中通过直接测量而获得的数据。在一些实施方式中，所述生产参数具体的可以包括产油量、产水量、产液量、生产压差等。实际应用中可以根据需要设置其他类型的生产参数，对此不做限制。
在一些实施方式中，在获取到地质参数和生产参数后，还可以确定各个开发地层的渗透率变异系数。由于各个开发地层的渗透率一般均会有所差异，因此可以通过渗透率变异系数来衡量不同开发地层的渗透率之间的差异性。
具体的，可以向获取各个开发地层的渗透率，并将其中最低渗透层的渗透率值作为基准渗透率，以Kmin表示。在计算渗透率变异系数时，可以先计算平均渗透率，再根据平均渗透率确定各层的渗透率变异系数。
在实际应用中，可以利用公式
Figure BDA0003388694520000041
计算渗透率变异系数，式中，VK为渗透率变异系数，n为开发地层的层数，Ki为开发地层的渗透率，
Figure BDA0003388694520000042
为平均渗透率，
Figure BDA0003388694520000043
hi为开发地层厚度。
在后续过程中对连通性进行分析时，可以结合所述渗透率变异系数来对地层的生产效果进行分析，进而帮助制定开发方案。
S120：确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数。
生产参数可以包括单层生产指数和合采生产指数。单层生产指数可以表示针对单一开发地层进行开发时所对应的生产参数，相应的，各个开发地层都可以对应有各自的单层生产指数。合采生产指数可以表示针对多个开发地层进行合采时所对应的生产参数，根据合采的开发地层的不同，合采生产指数也可以存在一定的差异，优选的，所述合采生产指数为针对当前所有开发地层同时进行开采时对应的生产指数，实际应用中也可以根据生产实际对所针对的开发地层进行调整，对此不做限制。
在获取到单层生产指数之后，可以确定单层参数影响关系。所述单层参数影响关系用于表示当前的地质参数和单层生产指数之间的对应关系，具体的可以通过拟合的方式来获取。
由于本说明书实施例主要针对的是合采情况，因此一般不会单独针对各个开发地层分别进行开发以获取所述单层生产指数，即所获取的单层生产指数一般情况下指代的是通过实验模拟所获取到的单层实验模拟数据。
具体的，所述单层实验模拟数据可以包括单层采油指数模拟数据、单层采液指数模拟数据中的至少一种，分别用于描述采油量的大小和采液量的大小。相应的，可以根据地质参数和单层采油指数模拟数据拟合得到单层采油指数影响关系，以及根据地质参数和单层采液指数模拟数据拟合得到单层采液指数影响关系。
利用一个具体的示例进行说明，根据单层采油指数模拟数据和地质参数构建的关系式可以为
Figure BDA0003388694520000051
式中，Jdo为单层采油指数模拟数据，ml/min×MPa-1，K为渗透率，10-3μm2，fw为含水率，％，am+3n为待拟合的系数。
此外，根据单层采液指数模拟数据和地质参数构建的关系式可以为
Figure BDA0003388694520000052
Figure BDA0003388694520000053
式中，Jdl为单层采液指数模拟数据，ml/min×MPa-1，K为渗透率，10-3μm2，fw为含水率，％，bm+4n为待拟合的系数。
在一个具体的实际应用中，可以得到上述关系式中待拟合的系数的具体数值如下表1所示。实际应用中拟合得到的系数值也可以与下表中的数值不同，对此不做限制。
表1
Figure BDA0003388694520000054
如图2A和图2B所示，为在单采情况下，确定出的含水率与采油指数、采液指数之间的对应关系，从而直观地表示参数之间的对应情况。
在拟合得到相应的参数后，即完成了关系式的构建，可以通过所获取到的地质参数来推算对应的采油指数或采液指数。
S130：确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数。
基于步骤S120中的说明，生产参数还可以包括合采生产指数。优选的，所述合采生产指数也可以是通过实验模拟而确定的生产参数。
所述合采实验模拟数据也可以包括合采采油指数模拟数据、合采采液指数模拟数据中的至少一种，在确定合采参数影响关系时，可以根据地质参数和合采采油指数模拟数据拟合得到合采采油指数影响关系，以及根据地质参数和合采采液指数模拟数据拟合得到合采采液指数影响关系。
在一些实施方式中，合采采油指数影响关系可以通过公式
Figure BDA0003388694520000055
Figure BDA0003388694520000056
进行计算，式中，Jho为合采采油指数模拟数据，ml/min×MPa-1，Kmax为最大渗透率，10-3μm2，Kmin为基准渗透率，10-3μm2，VK为渗透率变异系数，fw为含水率，％，ci(0≤i≤7)为待拟合的系数。
相应的，合采采液指数影响关系可以通过公式
Figure BDA0003388694520000061
Figure BDA0003388694520000062
进行计算，式中，Jhl为合采采液指数模拟数据，ml/min×MPa-1，Kmax为最大渗透率，10-3μm2，Kmin为基准渗透率，10-3μm2，VK为渗透率变异系数，fw为含水率，％，di(0≤i≤11)为待拟合的系数。
在一个具体的实际应用中，可以得到上述关系式中待拟合的系数的具体数值如下表2所示。实际应用中拟合得到的系数值也可以与下表中的数值不同，对此不做限制。
表2
Figure BDA0003388694520000063
如图3A和图3B所示，为在单采情况下，确定出的含水率与采油指数、采液指数之间的对应关系，从而直观地表示参数之间的对应情况。
在拟合得到相应的参数后，即完成了关系式的构建，可以通过所获取到的地质参数来推算合采过程所对应的采油指数或采液指数。
S140：结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度。
在获取到单层参数影响关系和合采参数影响关系，可以计算对应的连通性系数。
在一些实施方式中，可以基于地质参数，利用单层参数影响关系和合采参数影响关系分别计算实际单采指数和实际合采指数后，再拟合所述实际单采指数与地质参数之间的实际单采影响关系，并拟合所述实际合采指数与地质参数之间的实际合采影响关系，最终综合所述实际单采影响关系、实际合采影响关系得到连通性系数。
在根据上述过程获取到单层参数影响关系和合采参数影响关系后，可以根据上述关系对不同情况下的地质参数所对应的生产情况进行模拟，进而确定出相应的实际单采指数和实际合采指数。
所述实际单采影响关系包括实际单采米采油影响关系和实际单采米采液影响关系；所述实际合采影响关系包括实际合采米采油影响关系和实际合采米采液影响关系。相应的，在拟合关系时，也可以针对这些参数分别拟合对应的关系。
在一个具体的示例中，针对实际单采米采油影响关系的拟合公式可以为J′do＝-0.00147Ln(K)×fw+0.14585Ln(K)，式中，J′do为实际单采米采油指数，m3/d·MPa-1·m-1，K为渗透率，10-3μm2，fw为含水率，％。
针对实际单采米采液影响关系的拟合公式可以为
Figure BDA0003388694520000071
Figure BDA0003388694520000072
式中，J′dl为实际单采米采液指数，m3/d·MPa-1·m-1，K为渗透率，10-3μm2，fw为含水率，％。
针对合采状态下的合采参数影响关系，可以拟合得到公式如下
Figure BDA0003388694520000073
式中，Jhl为合采采液指数模拟数据，ml/min×MPa-1，Kmax为最大渗透率，10-3μm2，Kmin为基准渗透率，10-3μm2，VK为渗透率变异系数，fw为含水率，％，a1、b1、b2、c1、c2、c3、d1、d2、d3、d4、e1、e2、e3、e4、e5为待拟合的系数。上述待拟合的系数针对采油和采液可以具有不同的参数，以适应不同的情况。
具体的对于上述参数的拟合情况可以如下表3所示。实际应用中拟合得到的系数值也可以与下表中的数值不同，对此不做限制。
表3
Figure BDA0003388694520000074
在计算得到实际单采影响关系、实际合采影响关系，即可进行具体的计算得到连通性系数，具体的，所述连通性系数包括米采油连通性系数和米采液连通性系数，可以针对不同的情况分别进行计算。
在一些实施方式中，可以利用公式
Figure BDA0003388694520000075
Figure BDA0003388694520000081
计算米采油连通性系数，式中，α′0为米采油连通性系数，n为开发地层的层数，Jdoi′[fw(t)]为含水率为fw(t)时第i层的实际单采米采油指数，J′ho[fw(t+Δt)]为含水率为fw(t+Δt)时第i层的实际合采米采油指数。
相应的，可以利用公式
Figure BDA0003388694520000082
计算米采油连通性系数，式中，α′0为米采油连通性系数，n为开发地层的层数，Jdoi′[fw(t)]为含水率为fw(t)时第i层的实际单采米采油指数，J′ho[fw(t+Δt)]为含水率为fw(t+Δt)时第i层的实际合采米采油指数。
S150：基于所述连通性系数评估所述目标工区的连通性。
在获取到连通性系数后，可以利用所述连通性系数对目标工区的连通性进行评估。例如，可以预先设置相应的评估标准，将不同连通性系数的数值大小对应至不同的连通性评级，相应的连通性评级用于描述开发底层之间不同程度的连通性，具体的可以对应至不同的开发方案，从而在有效对工区中储层之间连通性进行评估的同时，指导了开发方案的有效制定。
在一些实施方式中，可以基于所述连通性系数和地质参数确定目标工区的连通性变化规律，以更为直观准确地对目标工区中各个开发底层之间的连通性进行评估。连通性变化规律用于描述连通性随其他参数的综合变化情况，例如，可以是描述连通性随地质参数和/或生产参数的变化情况。
具体的，可以如图4所示，利用连通性系数表征与基准渗透率、渗透率极差以及渗透率变异系数之间的关系，其中，基准渗透率选取范围为20mD～150mD，渗透率级差选取范围为2～20，渗透率变异系数变化范围为0.27～1.22。图4的结果表明，随着渗透率变异系数的增加，连通性系数随之增大；且基准渗透率越低，渗透率级差越大，连通性系数受渗透率变异系数的影响越剧烈；当基准渗透率处于低水平(20mD左右)且渗透率变异系数也处于低水平时，随着渗透率级差的增大，层段连通性系数变小，其原因在于低渗层段由于受到层段连通性影响而处于抑制状态，但同时低渗产能贡献比例也越小，因此，层段连通性系数变小；当基准渗透率处于中高水平时，随渗透率级差的增大，连通性系数变大；渗透率级差一定时，随着基准渗透率增大，连通性系数变大，且渗透率级差越大，连通性系数增长幅度越大。
结合上述变化规律，可以对实际生产开发的过程进行分析，进而拟定相应的开发方案。对于开发方案的具体制定过程可以基于实际应用的需求进行设置，在此不再赘述。
根据上述实施例的介绍，可以看出，所述方法在获取到开发地层的地质参数和生产参数后，首先通过实验模拟结果确定地质参数与生产参数之间的单层参数影响关系和合采参数影响关系，并进一步结合所述单层参数影响关系和合采参数影响关系计算连通性系数，从而有效对不同开发地层之间的连通性进行评估。通过上述方法，在考虑单采、合采以及实验模拟、实际开发的不同情况下，生产效果与地质的对应关系后，确定不同开发地层之间的连通性系数，保证了连通性系数的准确性以及与实际生产工况之间的关联性，提高了生产方案的开采效果，保证了生产开发的有效进行。
基于上述基于多层合采的连通性评估方法，本说明书还提出一种基于多层合采的连通性评估装置的实施例。所述基于多层合采的连通性评估装置可以设置于所述基于多层合采的连通性评估设备。如图5所示，所述基于多层合采的连通性评估装置具体包括以下模块。
参数获取模块510，用于获取对应于目标工区至少两个开发地层的地质参数和生产参数。
单层参数影响关系确定模块520，用于确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数。
合采参数影响关系确定模块530，用于确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数。
连通性系数计算模块540，用于结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度。
连通性评估模块550，用于基于所述连通性系数评估所述目标工区的连通性。
基于上述基于多层合采的连通性评估方法，本说明书实施例还提出一种基于多层合采的连通性评估设备。如图6所示，所述基于多层合采的连通性评估设备包括存储器和处理器。
在本实施例中，所述存储器可以按任何适当的方式实现。例如，所述存储器可以为只读存储器、机械硬盘、固态硬盘、或U盘等。所述存储器可以用于存储计算机程序指令。
在本实施例中，所述处理器可以按任何适当的方式实现。例如，处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application SpecificIntegrated Circuit，ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。所述处理器可以执行所述计算机程序指令实现以下步骤：获取对应于目标工区至少两个开发地层的地质参数和生产参数；确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；基于所述连通性系数评估所述目标工区的连通性。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
虽然上文描述的过程流程包括以特定顺序出现的多个操作，但是，应当清楚了解，这些过程可以包括更多或更少的操作，这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
本申请是参照根据本说明书实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
本领域技术人员应明白，本说明书的实施例可提供为方法、系统或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本说明书实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施例可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施例，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书实施例的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。

Claims (8)

1.一种基于多层合采的连通性评估方法，其特征在于，所述方法包括：
获取对应于目标工区至少两个开发地层的地质参数和生产参数；地质参数包括渗透率和含水率中的至少一种；
确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；
确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；
结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；
基于所述连通性系数评估所述目标工区的连通性；其中，所述确定地质参数与单层生产指数之间的单层参数影响关系之前，还包括：
确定各个开发地层的渗透率变异系数；所述渗透率变异系数用于表示开发地层的渗透率相对于平均渗透率的大小；
其中，所述确定各个开发地层的渗透率变异系数，包括：
利用公式
Figure FDA0003805864280000011
计算渗透率变异系数，式中，VK为渗透率变异系数，n为开发地层的层数，Ki为开发地层的渗透率，
Figure FDA0003805864280000012
为平均渗透率，
Figure FDA0003805864280000013
hi为开发地层厚度。
2.如权利要求1所述的方法，其特征在于，所述确定地质参数与单层生产指数之间的单层参数影响关系，包括：
获取针对单一开发地层进行开采时的单层实验模拟数据；所述单层实验模拟数据包括单层采油指数模拟数据、单层采液指数模拟数据中的至少一种；
根据地质参数和单层采油指数模拟数据拟合得到单层采油指数影响关系；
根据地质参数和单层采液指数模拟数据拟合得到单层采液指数影响关系。
3.如权利要求1所述的方法，其特征在于，所述确定地质参数与合采生产指数之间的合采参数影响关系，包括：
获取针对至少两个开发地层进行开采时的合采实验模拟数据；所述合采实验模拟数据包括合采采油指数模拟数据、合采采液指数模拟数据中的至少一种；
根据地质参数和合采采油指数模拟数据拟合得到合采采油指数影响关系；
根据地质参数和合采采液指数模拟数据拟合得到合采采液指数影响关系。
4.如权利要求1所述的方法，其特征在于，结合所述单层参数影响关系和合采参数影响关系计算连通性系数，包括：
基于地质参数，利用单层参数影响关系和合采参数影响关系分别计算实际单采指数和实际合采指数；
拟合所述实际单采指数与地质参数之间的实际单采影响关系；
拟合所述实际合采指数与地质参数之间的实际合采影响关系；
综合所述实际单采影响关系、实际合采影响关系得到连通性系数。
5.如权利要求4所述的方法，其特征在于，所述实际单采影响关系包括实际单采米采油影响关系和实际单采米采液影响关系；所述实际合采影响关系包括实际合采米采油影响关系和实际合采米采液影响关系；所述连通性系数包括米采油连通性系数和米采液连通性系数；
所述综合所述实际单采影响关系、实际合采影响关系得到连通性系数，包括：
利用公式
Figure FDA0003805864280000021
计算米采油连通性系数，式中，α′0为米采油连通性系数，n为开发地层的层数，Jdoi′[fw(t)]为含水率为fw(t)时第i层的实际单采米采油指数，J′ho[fw(t+Δt)]为含水率为fw(t+Δt)时第i层的实际合采米采油指数；
利用公式
Figure FDA0003805864280000022
计算米采油连通性系数，式中，α′0为米采油连通性系数，n为开发地层的层数，Jdoi′[fw(t)]为含水率为fw(t)时第i层的实际单采米采油指数，J′ho[fw(t+Δt)]为含水率为fw(t+Δt)时第i层的实际合采米采油指数。
6.如权利要求1所述的方法，其特征在于，所述基于所述连通性系数评估所述目标工区的连通性，包括：
基于所述连通性系数和地质参数确定目标工区的连通性变化规律；所述连通性变化规律用于描述连通性随地质参数和/或生产参数的变化情况；
利用所述连通性变化规律评估目标工区的连通性。
7.一种基于多层合采的连通性评估装置，其特征在于，所述装置包括：
参数获取模块，用于获取对应于目标工区至少两个开发地层的地质参数和生产参数；
单层参数影响关系确定模块，用于确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；
合采参数影响关系确定模块，用于确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；
连通性系数计算模块，用于结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；
连通性评估模块，用于基于所述连通性系数评估所述目标工区的连通性，
其中，所述确定地质参数与单层生产指数之间的单层参数影响关系之前，所述装置还执行以下方法：
确定各个开发地层的渗透率变异系数；所述渗透率变异系数用于表示开发地层的渗透率相对于平均渗透率的大小；
其中，所述确定各个开发地层的渗透率变异系数，包括：
利用公式
Figure FDA0003805864280000031
计算渗透率变异系数，式中，VK为渗透率变异系数，n为开发地层的层数，Ki为开发地层的渗透率，
Figure FDA0003805864280000032
为平均渗透率，
Figure FDA0003805864280000033
hi为开发地层厚度。
8.一种基于多层合采的连通性评估设备，包括存储器和处理器；
所述存储器，用于存储计算机程序指令；
所述处理器，用于执行所述计算机程序指令以实现以下步骤：获取对应于目标工区至少两个开发地层的地质参数和生产参数；确定地质参数与单层生产指数之间的单层参数影响关系；所述单层生产指数包括针对单一开发地层进行开发时的生产参数；确定地质参数与合采生产指数之间的合采参数影响关系；所述合采生产指数包括针对至少两个开发地层进行合采时的生产参数；结合所述单层参数影响关系和合采参数影响关系计算连通性系数；所述连通性系数用于表示不同开发地层之间的连通程度；基于所述连通性系数评估所述目标工区的连通性；
其中，所述确定地质参数与单层生产指数之间的单层参数影响关系之前，还包括：
确定各个开发地层的渗透率变异系数；所述渗透率变异系数用于表示开发地层的渗透率相对于平均渗透率的大小；
其中，所述确定各个开发地层的渗透率变异系数，包括：
利用公式
Figure FDA0003805864280000041
计算渗透率变异系数，式中，VK为渗透率变异系数，n为开发地层的层数，Ki为开发地层的渗透率，
Figure FDA0003805864280000042
为平均渗透率，
Figure FDA0003805864280000043
hi为开发地层厚度。
CN202111465675.XA 2021-12-02 2021-12-02 一种基于多层合采的连通性评估方法及装置 Active CN114139386B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111465675.XA CN114139386B (zh) 2021-12-02 2021-12-02 一种基于多层合采的连通性评估方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111465675.XA CN114139386B (zh) 2021-12-02 2021-12-02 一种基于多层合采的连通性评估方法及装置

Publications (2)

Publication Number Publication Date
CN114139386A CN114139386A (zh) 2022-03-04
CN114139386B true CN114139386B (zh) 2022-12-09

Family

ID=80387492

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111465675.XA Active CN114139386B (zh) 2021-12-02 2021-12-02 一种基于多层合采的连通性评估方法及装置

Country Status (1)

Country Link
CN (1) CN114139386B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112746836A (zh) * 2021-01-13 2021-05-04 重庆科技学院 基于层间干扰的油井各层产量计算方法
CN112749835A (zh) * 2020-12-22 2021-05-04 中国石油大学(北京) 一种储层产能预测方法、装置及设备
CN113486537A (zh) * 2021-08-04 2021-10-08 中国石油大学(北京) 一种基于气水交替驱油的井距设计方法、装置及设备

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112749835A (zh) * 2020-12-22 2021-05-04 中国石油大学(北京) 一种储层产能预测方法、装置及设备
CN112746836A (zh) * 2021-01-13 2021-05-04 重庆科技学院 基于层间干扰的油井各层产量计算方法
CN113486537A (zh) * 2021-08-04 2021-10-08 中国石油大学(北京) 一种基于气水交替驱油的井距设计方法、装置及设备

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
层间干扰问题的评价及解决方法研究;李禹等;《山东化工》;20180423(第08期);全文 *

Also Published As

Publication number Publication date
CN114139386A (zh) 2022-03-04

Similar Documents

Publication Publication Date Title
Peters et al. Results of the Brugge benchmark study for flooding optimization and history matching
US8862450B2 (en) Method of developing a petroleum reservoir from history matching of production data and seismic data
RU2590265C2 (ru) Системы и способы для оценки моментов прорыва флюида в местонахождениях добывающих скважин
FR2823877A1 (fr) Methode pour contraindre par des donnees dynamiques de production un modele fin representatif de la repartition dans le gisement d'une grandeur physique caracteristique de la structure du sous-sol
CN107451311B (zh) 一种油水相对渗透率曲线计算方法及装置
CN111665560A (zh) 油气层识别方法、装置、计算机设备及可读存储介质
CN106875286A (zh) 一种聚合物驱油田全过程注聚参数分阶段优化决策方法
CN110007344B (zh) 一种断溶体储层连通性的地震识别方法及装置
Al-Khazraji et al. Development of heterogeneous immature Brownfield with Waterdrive using dynamic opportunity index: a case study from Iraqi oilfields
Nakashima et al. Near-well upscaling for three-phase flows
CN106960264B (zh) 低渗透砂砾岩地层压力预测方法及装置
CN114139386B (zh) 一种基于多层合采的连通性评估方法及装置
US20150226061A1 (en) System and method for estimating flow capacity of a reservoir
CN111472763B (zh) 地层厚度预测方法及装置
Alghazal Development and testing of artificial neural network based models for water flooding and polymer gel flooding in naturally fractured reservoirs
CN113622908B (zh) 一种水侵气藏废弃地层压力的确定方法
Hasbollah et al. A preliminary basin scale evaluation framework of potential sedimentary basins in Malaysia for carbon dioxide sequestration
CN111582532A (zh) 应力敏感性油藏水平井产液能力预测方法及装置
Davudov et al. Machine learning assisted well placement optimization
CN113486537A (zh) 一种基于气水交替驱油的井距设计方法、装置及设备
Yadavalle et al. Use of pressure transient data to obtain geostatistical parameters for reservoir characterization
CN109033508A (zh) 一种油藏含水率及含水上升率的确定方法和装置
CN109063228A (zh) 一种水驱油藏含水上升率变化的确定方法及装置
CN111856567B (zh) 基于地震叠加速度的压实系数确定方法及装置
CN108920781A (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