CN113468476A - 一种聚合物驱增油量预测方法及装置 - Google Patents

一种聚合物驱增油量预测方法及装置 Download PDF

Info

Publication number
CN113468476A
CN113468476A CN202010244033.6A CN202010244033A CN113468476A CN 113468476 A CN113468476 A CN 113468476A CN 202010244033 A CN202010244033 A CN 202010244033A CN 113468476 A CN113468476 A CN 113468476A
Authority
CN
China
Prior art keywords
stage
polymer
polymer flooding
oil
water content
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
CN202010244033.6A
Other languages
English (en)
Other versions
CN113468476B (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010244033.6A priority Critical patent/CN113468476B/zh
Publication of CN113468476A publication Critical patent/CN113468476A/zh
Application granted granted Critical
Publication of CN113468476B publication Critical patent/CN113468476B/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
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK 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
    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; 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)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Data Mining & Analysis (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • General Business, Economics & Management (AREA)
  • Pure & Applied Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Mathematical Physics (AREA)
  • Marketing (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Operations Research (AREA)
  • Mathematical Optimization (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Evolutionary Biology (AREA)
  • Development Economics (AREA)
  • Software Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Databases & Information Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Quality & Reliability (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)

Abstract

本发明提供一种聚合物驱增油量预测方法及装置，所述方法包括：根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。所述装置用于执行上述方法。本发明实施例提供的聚合物驱增油量预测方法及装置，能够提高聚合物驱增油量预测的准确性。

Description

一种聚合物驱增油量预测方法及装置
技术领域
本发明涉及石油开采技术领域，具体涉及一种聚合物驱增油量预测方法及装置。
背景技术
水驱油田进入开发中后期，由于储层非均质性、流体流度差异以及注采差异等原因，地层中逐渐形成水驱优势通道，低效及无效注入水循环严重。由此油藏含水逐渐上升，而石油产量逐年下降。稳油控水和提高水驱效率成为中高含水非均质油藏亟需解决的难题。聚合物驱作为老油田稳产及增产的重要接替技术，在国内外油田都得到广泛应用。
聚合物驱是一种驱替流体流度控制技术。该技术主要是在注入水中掺入一定数量的聚丙烯酰胺或者黄原胶，以达到提高驱替相粘度和降低驱替相相对渗透率的目的。聚合物驱可以降低水油流度比，同时聚合物在岩石的吸附作用可增加水相的渗流阻力，较大幅度的扩大注入水的波及体积，改善粘性指进和舍进现象，以提高油藏的最终采收率。现有技术中，不同油藏的沉积特征、储层物性及流体性质不尽不同，不同油藏聚合物驱的开发特征也不相同，关于聚合物驱增油量预测主要以油藏数值模拟方法为主，其研究周期较长，且数值模拟历史拟合的好坏直接影响预测精度。
发明内容
针对现有技术中的问题，本发明实施例提供一种聚合物驱增油量预测方法及装置，能够至少部分地解决现有技术中存在的问题。
一方面，本发明提出一种聚合物驱增油量预测方法，包括：
根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；
根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；
根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
另一方面，本发明提供一种聚合物驱增油量预测装置，包括：
第一获得单元，用于根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；
第二获得单元，用于根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；
预测单元，用于根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述聚合物驱增油量预测方法的步骤。
又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述聚合物驱增油量预测方法的步骤。
本发明实施例提供的聚合物驱增油量预测方法及装置，能够根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量，根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量，根据聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量，由于通过聚合物特征曲线计算聚合物驱的累积产油量，能够提高聚合物驱增油量预测的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明一实施例提供的聚合物驱增油量预测方法的流程示意图。
图2是本发明另一实施例提供的聚合物驱增油量预测方法的流程示意图。
图3是本发明一实施例提供的注聚区块的聚合物驱特征曲线的示意图。
图4是本发明一实施例提供的目标区块的聚合物驱特征曲线的示意图。
图5是本发明一实施例提供的26-106井区聚合物驱和水驱的采出程度和含水率曲线。
图6是本发明一实施例提供的26-106井区的两种预测方法的采出程度的对比示意图。
图7是本发明一实施例提供的聚合物驱增油量预测装置的结构示意图。
图8是本发明另一实施例提供的聚合物驱增油量预测装置的结构示意图。
图9是本发明一实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
图1是本发明一实施例提供的聚合物驱增油量预测方法的流程示意图，如图1所示，本发明实施例提供的聚合物驱增油量预测方法，包括：
S101、根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；
具体地，对于目标区块，基于所述目标区块的聚合物的注入孔隙体积倍数和注聚阶段的采出程度，可以建立目标区块的聚合物特征曲线，聚合物特征曲线表明聚合物的注入孔隙体积倍数和阶段采出程度在半对数坐标系中呈现线性关系。然后根据目标区块的聚合物驱相关数据以及聚合物特征曲线，可以获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量。其中，所述目标区块是指采用聚合物驱进行石油开采的区块。所述聚合物驱相关数据包括但不限于目标区块的预测起点的聚合物累积注入量、目标区块的预测起点的聚合物日注入量、各个阶段的聚合物累积注入量，各个阶段的日注聚量、各个阶段的采出程度、各个阶段的时间步、油藏的总孔隙体积、目标区块的地质储量等，根据实际需要进行设置，本发明实施例不做限定。本发明实施例提供的聚合物驱增油量预测方法的执行主体包括但不限于计算机。
聚合物驱特征曲线可以表示为：
R＝alnQiPV+b (1)
其中，R为阶段采出程度，可以以百分数表示；QiPV为聚合物的注入孔隙体积倍数，a和b为回归系数。
S102、根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；
具体地，在获得每个阶段的聚合物驱日产油量之后，将每个阶段的聚合物驱日产油量与每个阶段的时间步相乘，获得每个阶段的累积产油量，再对各个阶段的累积产油量求和，将求和的结果作为所述目标区块的聚合物驱累积产油量。其中，每个阶段的时间步以天为基本单位，时间步可以设置为30天，根据实际需要进行设置，本发明实施例不做限定。
S103、根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
具体地，在获得所述目标区块的聚合物驱的累积产油量之后，可以计算所述目标区块的水驱累积产油量，然后用所述聚合物驱的累积产油量减去所述水驱累积产油量，获得所述目标区块的聚合物驱增油量。其中，与聚合物驱的累积产油量计算不同的是，水驱日产油量采用产量递减法计算的。水驱日产液量及含水率的计算过程与聚合物驱日产液量日产液量及含水率的计算过程相同。
例如，根据公式q'o＝qoie-Dt，计算获得所述目标区块的水驱日产油量，其中，q'o为水驱日产油量，单位可以采用吨/天，qoi为预测起始点的水驱日产油量，单位可以采用吨/天，t为预测时间，单位可以采用天，D为递减率，单位为1/天，递减率可以由拟合水驱历史生产数据得到。
本发明实施例提供的聚合物驱增油量预测方法，能够根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量，根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量，根据聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量，由于通过聚合物特征曲线计算聚合物驱的累积产油量，能够提高聚合物驱增油量预测的准确性。此外，相对于油藏数值模拟方法，提高了聚合物驱增油量预测的便捷性。准确地预测聚合物驱增油量还为聚合物驱开发规划及经济效益分析奠定基础。
图2是本发明另一实施例提供的聚合物驱增油量预测方法的流程示意图，如图2所示，在上述各实施例的基础上，进一步地，所述根据目标区块的聚合物驱相关数据，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量包括：
S1011、根据上一阶段的聚合物累积注入量、当前阶段的日注聚量、油藏的总孔隙体积以及聚合物驱特征曲线，计算获得当前阶段的采出程度；
具体地，通过上一阶段的聚合物累积注入量、当前阶段的日注聚量和油藏的总孔隙体积可以计算获得当前阶段的注入孔隙体积倍数，然后将当前阶段的注入孔隙体积倍数代入到聚合物驱特征曲线的表达式中，可以计算获得当前阶段的采出程度。其中，所述聚合物驱特征曲线是预先获得的。所述目标区块的预测起点的聚合物累积注入量和所述目标区块的预测起点的聚合物日注入量是预先获得的。
例如，采出程度计算公式为：
Figure BDA0002433480860000051
其中，Rn为第n个阶段的采出程度，可以以百分数表示，
Figure BDA0002433480860000052
为第n-1个阶段的聚合物累积注聚量，单位可以采用m3，
Figure BDA0002433480860000053
为第n个阶段的日注聚量，单位可以采用立方米/每天，Δt为第n个阶段的时间步，Vp为油藏的总孔隙体积，单位可以采用m3，a和b为所述目标区块的聚合物驱特征曲线的回归系数。
例如，在计算第1阶段的采出程度时，根据公式(2)获得：
Figure BDA0002433480860000054
其中，
Figure BDA0002433480860000055
为第0个阶段的聚合物累积注聚量，即所述目标区块的预测起点的聚合物累积注入量，
Figure BDA0002433480860000056
为第0个阶段的日注聚量，即所述目标区块的预测起点的聚合物日注入量。
S1012、根据所述当前阶段的采出程度、上一阶段的采出程度、所述目标区块的地质储量以及所述当前阶段的时间步，获得所述当前阶段的聚合物驱日产油量；
具体地，在获得所述当前阶段的采出程度之后，可以根据所述当前阶段的采出程度、上一阶段的采出程度、所述目标区块的地质储量以及所述当前阶段的时间步，获得所述当前阶段的聚合物驱日产油量。
例如，所述当前阶段的聚合物驱日产油量可以根据如下公式计算获得：
Figure BDA0002433480860000061
其中，
Figure BDA0002433480860000062
为第n个阶段的聚合物驱日产油量，单位可以采用吨/天，Rn为第n个阶段的采出程度，可以以百分数表示，Rn-1为第n-1个阶段的采出程度，可以以百分数表示，N为所述目标区块的地质储量，单位可以采用吨，Δt为第n个阶段的时间步。
S1013、根据所述当前阶段的聚合物驱日产油量、所述当前阶段的日产液量以及含水率计算公式，计算获得所述当前阶段的含水率；
具体地，在获得当前阶段的聚合物驱日产油量之后，可以根据所述当前阶段的聚合物驱日产油量、所述当前阶段的日产液量以及含水率计算公式，计算获得所述当前阶段的含水率。其中，所述含水率计算公式是预设的。
例如，所述含水率计算公式为：
Figure BDA0002433480860000063
其中，
Figure BDA0002433480860000064
为第n个阶段的含水率，
Figure BDA0002433480860000065
为第n个阶段的聚合物驱日产油量，
Figure BDA0002433480860000068
为第n个阶段的日产液量，单位可以采用立方米/每天。
日产液量可以根据如下计算公式获得：
Figure BDA0002433480860000066
其中，
Figure BDA0002433480860000067
为第n个阶段的日产液量，单位可以采用立方米/每天，zn为第n个阶段的注采比，无量纲，Bo为原油体积系数，单位可以采用m3/m3，ρo为原油密度，单位可以采用g/cm3。
S1014、若判断获知所述当前阶段的含水率小于所述极限含水率，则计算下一阶段的聚合物驱日产油量和下一阶段的含水率，直到所述当前阶段的含水率大于等于极限含水率。
具体地，在获得所述当前阶段的含水率之后，将当前阶段的含水率与所述极限含水率件比较，如果所述当前阶段的含水率小于所述极限含水率，说明目标区块还可以继续开采，计算下一阶段的聚合物驱日产油量和下一阶段的含水率，所述下一阶段的聚合物驱日产油量与和下一阶段的含水率与所述当前阶段的聚合物驱日产油量以及当前阶段的含水率的计算方法类似，此处不进行赘述。如果所述当前阶段的含水率大于等于所述极限含水率，说明目标区块无法继续采用聚合物驱进行开采，无需计算下一阶段的聚合物驱日产油量以及下一阶段的含水率。
例如，从第1阶段开始，采用公式(2)和(3)可以计算出所述目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量，采用公式(4)和(5)可以计算出所述目标区块的含水率达到极限含水率之前各个阶段的含水率。
在上述各实施例的基础上，进一步地，所述聚合物驱特征曲线是基于所述目标区块的聚合物驱阶段采出程和聚合物注入孔隙体积倍数通过曲线拟合获得的。
具体地，在所述目标区块采用聚合物驱进行采油时，获得聚合物驱采油过程中聚合物的注入孔隙体积倍数和阶段采出程度，在半对数坐标系中以聚合物的注入孔隙体积倍数为横坐标，以阶段采出程度为纵坐标，绘制出坐标点，然后进行曲线拟合，可以获得所述目标区块的聚合物驱特征曲线。
例如，统计分析国内某油田已经结束的40多个注聚合物区块(简称注聚区块)，结果发现在聚合物驱过程中，聚合物的注入孔隙体积倍数和阶段采出程度在半对数坐标系中呈现线性关系，该线性关系称之为聚合物驱特征曲线。图3是本发明一实施例提供的注聚区块的聚合物驱特征曲线的示意图，如图3所示，在半对数坐标系上，当注聚区块的注入孔隙体积倍数大于0.22后，聚合物的注入孔隙体积倍数和阶段采出程度呈现出线性关系，且线性相关性系数高达0.9995。
下面以本发明实施例提供的聚合物驱增油量预测方法在Nuraly油田的应用为例，来说明所述聚合物驱增油量预测方法的具体实现过程。
Nuraly油田，是一个多层状砂岩油藏。主力层M-II-3层以粉砂岩和细砂岩为主，含少量砾石，属中孔中渗储层，平均孔隙度为17.1％，渗透率为426mD，为扇三角洲前缘沉积。由于M-II-3层河道砂纵向非均质性强，局部发育优势通道，水驱纵向动用程度仅为50％，即纵向上将近一半的油层还处于未动用状态。为此，M-II-3层于2014年在26-106井区开展聚合物驱，预测起点的注聚合物日注入量354m3/d，累积注入聚合物溶液0.12PV。自26-106井区实施聚合物驱以来，注入井吸水剖面得到明显改善，高渗层吸水量减少，低渗层吸水量增加，吸水剖面非均质性减弱。自开展聚合物驱以来，26-106井区含水稳定，日产油、日产液均大幅提高，产量递减分析表明井组累增油8.9万吨，平均单井日增油7.8t/d，聚合物驱开发效果较好。
以26-106井区为目标区块，将26-106井区的聚合物的注入孔隙体积倍数与阶段采出程度绘制在半对数坐标系内，如图4所示，获得26-106井区在注聚过程中存在的一条聚合物驱特征曲线，其表达式为：
R＝3.112lnQiPV+13.051 (6)
其中，R为阶段采出程度，QiPV为聚合物的注入孔隙体积倍数，3.112和13.051为回归系数。
基于本发明实施例提供的聚合物驱增油量预测方法，便可以计算出26-106井区按注采比0.6和注聚速度354m3/d条件下聚合物驱和水驱的采出程度及含水率的变化曲线，如图5所示。极限含水98％条件下，26-106井区聚合物驱阶段采出程度为10.2％，可以计算出26-106井区累积产油量为19.8万吨；而普通水驱的阶段采出程度仅为1.5％，累产油2.9万吨。聚合物驱比普通水驱提高采出程度8.7％，可以预测出26-106井区的聚合物驱增油量为16.9万吨。截止到数据统计日，26-106井区已注入0.12PV的聚合物，聚合物驱至含水率达到98％时，未来将还会有7.6万吨的原油被采出。另外，当聚合物的注入孔隙体积倍数大于0.3后，聚合物驱的阶段采出程度增加缓慢，因此建议聚合物驱的总注入量为0.3PV。
对于26-106井区，采用本发明实施例提供的聚合物驱增油量预测方法预测各阶段的采出程度，并采用油藏数值模拟的方法计算各阶段的采出程度，以上述两种方法获得的各阶段的采出程度相对于水驱下各阶段的采出程度的提高量为纵坐标，以聚合物的注入孔隙体积倍数为横坐标，绘制出图6所示的各阶段的采出程度的对比示意图。从图6中可以看出，本发明实施例提供的方法预测至含水98％时，累积注入孔隙体积倍数为0.382PV，聚合物驱阶段的采出程度为10.2％；而数值模拟预测至含水98％时，累积注入孔隙体积倍数为0.390PV，聚合物驱阶段的采出程度为10.05％，略低于10.2％。
图7是本发明一实施例提供的聚合物驱增油量预测装置的结构示意图，如图7所示，本发明实施例提供的聚合物驱增油量预测装置包括第一获得单元701、第二获得单元702和预测单元703，其中：
第一获得单元701用于根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；第二获得单元702用于根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；预测单元703用于根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
具体地，对于目标区块，基于所述目标区块的聚合物的注入孔隙体积倍数和注聚阶段的采出程度，可以建立目标区块的聚合物特征曲线，聚合物特征曲线表明聚合物的注入孔隙体积倍数和阶段采出程度在半对数坐标系中呈现线性关系。第一获得单元701根据目标区块的聚合物驱相关数据以及聚合物特征曲线，可以获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量。其中，所述目标区块是指采用聚合物驱进行石油开采的区块。所述聚合物驱相关数据包括但不限于目标区块的预测起点的聚合物累积注入量、目标区块的预测起点的聚合物日注入量、各个阶段的聚合物累积注入量，各个阶段的日注聚量、各个阶段的采出程度、各个阶段的时间步、油藏的总孔隙体积、目标区块的地质储量等，根据实际需要进行设置，本发明实施例不做限定。
在获得每个阶段的聚合物驱日产油量之后，第二获得单元702将每个阶段的聚合物驱日产油量与每个阶段的时间步相乘，获得每个阶段的产油量，再对各个阶段的产油量求和，将求和的结果作为所述目标区块的聚合物驱累积产油量。其中，每个阶段的时间步以天为基本单位，时间步可以设置为7天，根据实际需要进行设置，本发明实施例不做限定。
在获得所述目标区块的聚合物驱的累积产油量之后，预测单元703可以计算所述目标区块的水驱累积产油量，然后用所述聚合物驱的累积产油量减去所述水驱累积产油量，获得所述目标区块的聚合物驱增油量。其中，与聚合物驱的累积产油量计算不同的是，水驱日产油量采用产量递减法计算的。水驱日产液量及含水率的计算过程与聚合物驱日产液量日产液量及含水率的计算过程相同。
本发明实施例提供的聚合物驱增油量预测装置，能够根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量，根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量，根据聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量，由于通过聚合物特征曲线计算聚合物驱的累积产油量，能够提高聚合物驱增油量预测的准确性。此外，相对于油藏数值模拟方法，提高了聚合物驱增油量预测的便捷性。准确地预测聚合物驱增油量还为聚合物驱开发规划及经济效益分析奠定基础。
图8是本发明另一实施例提供的聚合物驱增油量预测装置的结构示意图，如图8所示，在上述各实施例的基础上，进一步地，第一获得单元701包括第一计算子单元7011、获得子单元7012、第二计算子单元7013和判断子单元7014，其中：
第一计算子单元7011用于根据上一阶段的聚合物累积注入量、当前阶段的日注聚量、油藏的总孔隙体积以及聚合物驱特征曲线，计算获得当前阶段的采出程度；获得子单元7012用于根据所述当前阶段的采出程度、上一阶段的采出程度、所述目标区块的地质储量以及所述当前阶段的时间步，获得所述当前阶段的聚合物驱日产油量；第二计算子单元7013用于根据所述当前阶段的聚合物驱日产油量、所述当前阶段的日产液量以及含水率计算公式，计算获得所述当前阶段的含水率；判断子单元7014用于在判断获知所述当前阶段的含水率小于所述极限含水率之后，计算下一阶段的聚合物驱日产油量和下一阶段的含水率，直到所述当前阶段的含水率大于等于所述极限含水率。
具体地，第一计算子单元7011通过上一阶段的聚合物累积注入量、当前阶段的日注聚量和油藏的总孔隙体积可以计算获得当前阶段的注入孔隙体积倍数，然后将当前阶段的注入孔隙体积倍数代入到聚合物驱特征曲线的表达式中，可以计算获得当前阶段的采出程度。其中，所述聚合物驱特征曲线是预先获得的。所述目标区块的预测起点的聚合物累积注入量和所述目标区块的预测起点的聚合物日注入量是预先获得的。
在获得所述当前阶段的采出程度之后，获得子单元7012可以根据所述当前阶段的采出程度、上一阶段的采出程度、所述目标区块的地质储量以及所述当前阶段的时间步，获得所述当前阶段的聚合物驱日产油量。
在获得当前阶段的聚合物驱日产油量之后，第二计算子单元7013可以根据所述当前阶段的聚合物驱日产油量、所述当前阶段的日产液量以及含水率计算公式，计算获得所述当前阶段的含水率。其中，所述含水率计算公式是预设的。
在获得所述当前阶段的含水率之后，判断子单元7014将当前阶段的含水率与所述极限含水率件比较，如果所述当前阶段的含水率小于所述极限含水率，说明目标区块还可以继续开采，计算下一阶段的聚合物驱日产油量和下一阶段的含水率，所述下一阶段的聚合物驱日产油量与和下一阶段的含水率与所述当前阶段的聚合物驱日产油量以及当前阶段的含水率的计算方法类似，此处不进行赘述。如果所述当前阶段的含水率大于等于所述极限含水率，说明目标区块无法继续采用聚合物驱进行开采，无需计算下一阶段的聚合物驱日产油量以及下一阶段的含水率。
在上述各实施例的基础上，进一步地，所述当前阶段的日产液量根据如下公式计算获得：
Figure BDA0002433480860000111
其中，
Figure BDA0002433480860000112
为第n个阶段的日产液量，
Figure BDA0002433480860000113
为第n个阶段的日注聚量，zn为第n个阶段的注采比，Bo为原油体积系数，
Figure BDA0002433480860000114
为第n个阶段的聚合物驱日产油量，ρo为原油密度。
在上述各实施例的基础上，进一步地，所述聚合物驱特征曲线是基于所述目标区块的聚合物驱阶段采出程和聚合物的注入孔隙体积倍数通过曲线拟合获得的。
具体地，在所述目标区块采用聚合物驱进行采油时，获得聚合物驱采油过程中聚合物的注入孔隙体积倍数和阶段采出程度，在半对数坐标系中以聚合物的注入孔隙体积倍数为横坐标，以阶段采出程度为纵坐标，绘制出坐标点，然后进行曲线拟合，可以获得所述目标区块的聚合物驱特征曲线。
本发明实施例提供的装置的实施例具体可以用于执行上述各方法实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
图9是本发明一实施例提供的电子设备的实体结构示意图，如图9所示，该电子设备可以包括：处理器(processor)901、通信接口(Communications Interface)902、存储器(memory)903和通信总线904，其中，处理器901，通信接口902，存储器903通过通信总线904完成相互间的通信。处理器901可以调用存储器903中的逻辑指令，以执行如下方法：根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
此外，上述的存储器903中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种聚合物驱增油量预测方法，其特征在于，包括：
根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；
根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；
根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
2.根据权利要求1所述的方法，其特征在于，所述根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量包括：
根据上一阶段的聚合物累积注入量、当前阶段的日注聚量、油藏的总孔隙体积以及聚合物驱特征曲线，计算获得当前阶段的采出程度；
根据所述当前阶段的采出程度、上一阶段的采出程度、所述目标区块的地质储量以及所述当前阶段的时间步，获得所述当前阶段的聚合物驱日产油量；
根据所述当前阶段的聚合物驱日产油量、所述当前阶段的日产液量以及含水率计算公式，计算获得所述当前阶段的含水率；
若判断获知所述当前阶段的含水率小于所述极限含水率，则计算下一阶段的聚合物驱日产油量和下一阶段的含水率，直到所述当前阶段的含水率大于等于所述极限含水率。
3.根据权利要求2所述的方法，其特征在于，所述当前阶段的日产液量根据如下公式计算获得：
Figure FDA0002433480850000011
其中，
Figure FDA0002433480850000012
为第n个阶段的日产液量，
Figure FDA0002433480850000013
为第n个阶段的日注聚量，zn为第n个阶段的注采比，Bo为原油体积系数，
Figure FDA0002433480850000014
为第n个阶段的聚合物驱日产油量，ρo为原油密度。
4.根据权利要求1至3任一项所述的方法，其特征在于，所述聚合物驱特征曲线是基于所述目标区块的聚合物驱阶段采出程和聚合物的注入孔隙体积倍数通过曲线拟合获得的。
5.一种聚合物驱增油量预测装置，其特征在于，包括：
第一获得单元，用于根据目标区块的聚合物驱相关数据以及聚合物特征曲线，获得目标区块的含水率达到极限含水率之前各个阶段的聚合物驱日产油量；
第二获得单元，用于根据每个阶段的聚合物驱日产油量以及每个阶段的时间步，获得所述目标区块的聚合物驱累积产油量；
预测单元，用于根据所述聚合物驱的累积产油量以及水驱累积产油量，预测所述目标区块的聚合物驱增油量。
6.根据权利要求5所述的装置，其特征在于，所述第一获得单元包括：
第一计算子单元，用于根据上一阶段的聚合物累积注入量、当前阶段的日注聚量、油藏的总孔隙体积以及聚合物驱特征曲线，计算获得当前阶段的采出程度；
获得子单元，用于根据所述当前阶段的采出程度、上一阶段的采出程度、所述目标区块的地质储量以及所述当前阶段的时间步，获得所述当前阶段的聚合物驱日产油量；
第二计算子单元，用于根据所述当前阶段的聚合物驱日产油量、所述当前阶段的日产液量以及含水率计算公式，计算获得所述当前阶段的含水率；
判断子单元，用于在判断获知所述当前阶段的含水率小于所述极限含水率之后，计算下一阶段的聚合物驱日产油量和下一阶段的含水率，直到所述当前阶段的含水率大于等于所述极限含水率。
7.根据权利要求6所述的装置，其特征在于，所述当前阶段的日产液量根据如下公式计算获得：
Figure FDA0002433480850000021
其中，
Figure FDA0002433480850000022
为第n个阶段的日产液量，
Figure FDA0002433480850000023
为第n个阶段的日注聚量，zn为第n个阶段的注采比，Bo为原油体积系数，
Figure FDA0002433480850000024
为第n个阶段的聚合物驱日产油量，ρo为原油密度。
8.根据权利要求5至7任一项所述的装置，其特征在于，所述聚合物驱特征曲线是基于所述目标区块的聚合物驱阶段采出程和聚合物的注入孔隙体积倍数通过曲线拟合获得的。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至4任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至4任一项所述方法的步骤。
CN202010244033.6A 2020-03-31 2020-03-31 一种聚合物驱增油量预测方法及装置 Active CN113468476B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010244033.6A CN113468476B (zh) 2020-03-31 2020-03-31 一种聚合物驱增油量预测方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010244033.6A CN113468476B (zh) 2020-03-31 2020-03-31 一种聚合物驱增油量预测方法及装置

Publications (2)

Publication Number Publication Date
CN113468476A true CN113468476A (zh) 2021-10-01
CN113468476B CN113468476B (zh) 2023-07-28

Family

ID=77865533

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010244033.6A Active CN113468476B (zh) 2020-03-31 2020-03-31 一种聚合物驱增油量预测方法及装置

Country Status (1)

Country Link
CN (1) CN113468476B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116201538A (zh) * 2023-03-15 2023-06-02 西南石油大学 一种基于采出程度的全生命周期储层损害评价方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100300682A1 (en) * 2009-05-27 2010-12-02 Ganesh Thakur Computer-implemented systems and methods for screening and predicting the performance of enhanced oil recovery and improved oil recovery methods
CN105089589A (zh) * 2015-07-27 2015-11-25 中国石油大学(华东) 一种实现聚合物驱合层与分层注入的可视化实验装置及其工作方法
CN105225019A (zh) * 2015-11-11 2016-01-06 四川建筑职业技术学院 特高含水期累积产油量及采收率预测方法
CN106014359A (zh) * 2016-06-08 2016-10-12 西南石油大学 一种海上油藏早期注聚最早转注时机判断方法
CN110610288A (zh) * 2019-08-01 2019-12-24 裴雪皓 一种油气井生产数据智能系统分析方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100300682A1 (en) * 2009-05-27 2010-12-02 Ganesh Thakur Computer-implemented systems and methods for screening and predicting the performance of enhanced oil recovery and improved oil recovery methods
CN105089589A (zh) * 2015-07-27 2015-11-25 中国石油大学(华东) 一种实现聚合物驱合层与分层注入的可视化实验装置及其工作方法
CN105225019A (zh) * 2015-11-11 2016-01-06 四川建筑职业技术学院 特高含水期累积产油量及采收率预测方法
CN106014359A (zh) * 2016-06-08 2016-10-12 西南石油大学 一种海上油藏早期注聚最早转注时机判断方法
CN110610288A (zh) * 2019-08-01 2019-12-24 裴雪皓 一种油气井生产数据智能系统分析方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116201538A (zh) * 2023-03-15 2023-06-02 西南石油大学 一种基于采出程度的全生命周期储层损害评价方法
CN116201538B (zh) * 2023-03-15 2023-09-12 西南石油大学 一种基于采出程度的全生命周期储层损害评价方法

Also Published As

Publication number Publication date
CN113468476B (zh) 2023-07-28

Similar Documents

Publication Publication Date Title
Yuncong et al. Performance and gas breakthrough during CO2 immiscible flooding in ultra-low permeability reservoirs
CN105587298B (zh) 多流管模式的水驱油井含水率反演方法
CN102287172B (zh) 一种稠油油藏采油方法
CN104915512B (zh) 一种预测油田采出程度与含水率的方法
CN105626006A (zh) 低渗透油藏co2驱技术极限井距确定方法
CN103939066B (zh) 一种一注多采井组定注水量确定油井产液量的方法
CN109958413B (zh) 一种特高含水期油藏动态流动单元划分方法
CN109858177A (zh) 一种基于拟流线法的底水油藏水平井水驱数值模拟方法
CN104392091B (zh) 一种气驱开发下的油藏工程动态分析方法
CN103967458A (zh) 一种防砂段水驱方法
CN111810101A (zh) 一种水驱油藏动态分析方法及装置
Pan et al. Twelve years field applications of offshore heavy oil polymer flooding from continuous injection to alternate injection of polymer-water
CN110984974B (zh) 基于水侵速度确定有水气藏合理采气速度的方法
CN113468476B (zh) 一种聚合物驱增油量预测方法及装置
CN109918769A (zh) 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法
CN114528780A (zh) 一种考虑稠油非牛顿特性的数值模拟方法
CN109296363A (zh) 特低渗透油藏co2驱初期产能预测方法
CN104675372B (zh) 一种用于聚合物驱注入采出量分配的方法
CN109915123A (zh) 利用瞬时方程计算缝洞型油藏拟稳态水侵水侵量的方法
CN113326465B (zh) 油藏开发动态分析方法及装置
CN112324413B (zh) 一种提高注入井注入量的化学施工方法
CN112528573B (zh) 一种特高含水期油井的堵剂参数优化方法及系统
CN111322054B (zh) 一种砂岩油藏在化学驱阶段三三结合挖潜优化开采方法
CN116025331A (zh) 一种无效水识别及量化表征方法
CN111456695A (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