CN112487594B - 油藏水体倍数计算方法及装置 - Google Patents

油藏水体倍数计算方法及装置 Download PDF

Info

Publication number
CN112487594B
CN112487594B CN201910782014.6A CN201910782014A CN112487594B CN 112487594 B CN112487594 B CN 112487594B CN 201910782014 A CN201910782014 A CN 201910782014A CN 112487594 B CN112487594 B CN 112487594B
Authority
CN
China
Prior art keywords
reservoir
water
oil
invasion
water invasion
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
CN201910782014.6A
Other languages
English (en)
Other versions
CN112487594A (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 CN201910782014.6A priority Critical patent/CN112487594B/zh
Publication of CN112487594A publication Critical patent/CN112487594A/zh
Application granted granted Critical
Publication of CN112487594B publication Critical patent/CN112487594B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Data Mining & Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Fluid Mechanics (AREA)
  • Tourism & Hospitality (AREA)
  • Agronomy & Crop Science (AREA)
  • Strategic Management (AREA)
  • Primary Health Care (AREA)
  • Marine Sciences & Fisheries (AREA)
  • General Business, Economics & Management (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Marketing (AREA)
  • Animal Husbandry (AREA)
  • General Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Operations Research (AREA)
  • Health & Medical Sciences (AREA)
  • Algebra (AREA)
  • Environmental & Geological Engineering (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种油藏水体倍数计算方法及装置，该方法包括：获取油藏地质特征参数和油藏生产开发参数；根据油藏地质特征参数、油藏生产开发参数和物质平衡方程计算第一油藏水侵量；第一油藏水侵量包括不同时刻的油藏水侵量的值；根据油藏地质特征参数、油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量；第二油藏水侵量包括不同时刻的油藏水侵量的值；若每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值，则根据第二油藏水侵量所使用的预设水体体积的值计算油藏水体倍数。本发明计算过程简单，能够快速计算出油藏的水体倍数，可为油藏天然能量评价提供数据依据。

Description

油藏水体倍数计算方法及装置
技术领域
本发明涉及地球物理勘探技术领域，尤其是涉及一种油藏水体倍数计算方法及装置。
背景技术
能量是油气田开发的核心问题，是决定油气田最终采收率的主要因素。油气藏的天然能量一般有水驱、弹性驱、溶解气驱、气顶驱和重力驱。其中水驱又包含弹性水驱和刚性水驱，属于油气藏的外部天然能量；而弹性驱、溶解气驱、气顶驱和重力驱等属于油气藏的内部天然能量。
国内外油藏的开发实践表明，很多油藏都与外部水体相互连通，该油藏的驱动能量包括外部水体的压能和弹性能。这些水体可能是具有外部供给的无限水体，也可能是存在边界的有限水体。在油藏的开发过程中，随着地层原油及天然气的不断采出，油藏地层压力不断下降。受油藏与外部水体之间压力差的影响，外部水体的水会不断侵入油藏，在补充油藏能量的同时，也会造成油水边界附近的油井过快水淹，从而对油藏的开发效果造成显著的影响。
当油藏的外部水体较小时，油藏开发引起的压力降可以很快波及到整个水体范围，外部水体的侵入无法完全补充油藏的地层亏空，因此油藏的地层压力会不断迅速下降；而当油藏的水体较大时，油藏的压力降仅能波及到油藏附近的部分水体，仍有一部分天然水体保持原始地层压力，外部水体可弥补大部分的油藏地层亏空，因此油藏的地层压力会趋于稳定。因此，外部水体的大小对油藏地层压力的影响较大。
对于外部水体较小的油藏，为了保持地层压力稳定，可通过人工注水来补充地层能量，改善油气藏的整体开发效果。人工注水自20世纪二三十年代在美国获得工业化应用以来，在全世界范围内得到了广泛应用，已经成为目前油气田开发的主要开采方式。按照注水井在油藏中的位置和注水井与生产井之间的排列关系，人工注水可以分为边缘注水、切割注水以及面积注水。
可见，水体大小计算是油藏动态分析、开发调整方案制定及稳油控水措施实施的重要基础性工作，而现有油藏水体倍数计算方法计算过程复杂，适用性差。
发明内容
本发明提供了一种油藏水体倍数计算方法及装置，可以简单、快速定量评价油藏水体大小，为天然水驱油藏开发提供技术支撑。
第一方面，本发明实施例提供了一种油藏水体倍数计算方法，该方法包括：获取油藏地质特征参数和油藏生产开发参数；根据所述油藏地质特征参数、所述油藏生产开发参数和物质平衡方程计算第一油藏水侵量；所述第一油藏水侵量包括不同时刻的油藏水侵量的值；根据所述油藏地质特征参数、所述油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量；所述第二油藏水侵量包括不同时刻的油藏水侵量的值；根据不同的所述预设水体体积生成的第二油藏水侵量不同；若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数。
第二方面，本发明实施例还提供一种油藏水体倍数计算装置，该装置包括：获取模块，用于获取油藏地质特征参数和油藏生产开发参数；第一水侵量计算模块，用于根据所述油藏地质特征参数、所述油藏生产开发参数和物质平衡方程计算第一油藏水侵量；所述第一油藏水侵量包括不同时刻的油藏水侵量的值；第二水侵量计算模块，用于根据所述油藏地质特征参数、所述油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量；所述第二油藏水侵量包括不同时刻的油藏水侵量的值；根据不同的所述预设水体体积生成的第二油藏水侵量不同；油藏水体倍数计算模块，用于若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数。
第三方面，本发明实施例还提供一种计算机设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述油藏水体倍数计算方法。
第四方面，本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，所述程序代码使所述处理器执行上述油藏水体倍数计算方法。
本发明实施例带来了以下有益效果：本发明实施例提供了一种油藏水体倍数计算方法及装置，该方法首先获取油藏地质特征参数和油藏生产开发参数，再分别结合物质平衡方程和Fetkovitch方程计算第一油藏水侵量和第二油藏水侵量，第一油藏水侵量和第二油藏水侵量均包括不同时刻的油藏水侵量的值，基于不同时刻的油藏水侵量的值对第一油藏水侵量和第二油藏水侵量的关系进行分析，当每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值时，根据第二油藏水侵量所使用的预设水体体积计算油藏水体倍数。本发明实施例计算过程简单，能够快速计算出油藏的水体倍数，可为油藏天然能量评价提供数据依据。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的油藏水体倍数计算方法流程图；
图2为本发明实施例提供的边水油藏边缘注水开发示意图；
图3为本发明实施例提供的油藏水侵量交会图示意图；
图4为本发明实施例提供的油藏累积水侵量随开采时间的变化曲线；
图5为本发明实施例提供的一种油藏水体倍数计算装置结构框图；
图6为本发明实施例提供的另一种油藏水体倍数计算结构框图；
图7为本发明实施例提供的计算机设备结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前，关于水体大小的计算方法主要是基于非稳态水侵计算，需要通过试凑法确定无因次水侵量函数，再根据无因次半径确定外部水体的大小。该方法计算过程复杂，而且结果偏差较大，实用性较差。
本发明实施例提供的一种油藏水体倍数计算方法及装置，基于物质平衡原理和Fetkovitch水侵量方程，建立一种新的油藏水体计算模型，以为油藏水体计算提供借鉴。
为便于对本实施例进行理解，首先对本发明实施例所公开的一种油藏水体倍数计算方法进行详细介绍。
本发明实施例提供了一种油藏水体倍数计算方法，参见图1所示的一种油藏水体倍数计算方法流程图，该方法包括以下步骤：
步骤S102，获取油藏地质特征参数和油藏生产开发参数。
在本发明实施例中，油藏地质特征参数是指用于描述油藏所在位置的地层压力、体积、面积以及温度等与油藏开发相关的地质特征的数据。油藏生产开发参数包括与油藏开发相关的生产数据，例如采油井参数以及注水井参数等。
步骤S104，根据油藏地质特征参数、油藏生产开发参数和物质平衡方程计算第一油藏水侵量。
在本发明实施例中，物质平衡方程是指根据物质平衡原理确定的可以用于计算油藏水侵量的方程。在计算的过程中，将油藏地质特征参数和油藏生产开发参数作为物质平衡方程中的已知量，可以计算出第一油藏水侵量。
需要说明的是，第一油藏水侵量包括多个不同时刻的油藏水侵量的值。例如，在时刻1，第一油藏水侵量的值为a1；在时刻2，第一油藏水侵量的值为a2；在时刻3，第一油藏水侵量的值为a3。
步骤S106，根据油藏地质特征参数、油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量。
在本发明实施例中，使用不同的预设水体体积的值可以得到不同的第二油藏水侵量。
需要说明的是，第二油藏水侵量包括多个不同时刻的油藏水侵量的值。例如，在时刻1，第二油藏水侵量的值为b1；在时刻2，第二油藏水侵量的值为b2；在时刻3，第二油藏水侵量的值为b3。
步骤S108，若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数。
在本发明实施例中，预设值是可以根据实际需求或者经验进行设置或修改的，本发明实施例对具体的预设条件不作具体限定。计算出第一油藏水侵量和多个第二油藏水侵量之后，对各个时刻对应的第一油藏水侵量的值和多个第二油藏水侵量的值进行比较：例如，对于时刻1，第一油藏水侵量的值为a1，第二油藏水侵量的值为b1比较a1和b1是否相同或者二者的差值是否小于预设值，然后，对于时刻2的油藏水侵量的值也进行比较，在对各个时刻的油藏水侵量的值进行比较之后，若存在一个第二水侵量使得每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值，则将该第二油藏水侵量使用的预设水体体积的值作为计算油藏水体倍数的一个参数。
本发明实施例提供了一种油藏水体倍数计算方法，该方法首先获取油藏地质特征参数和油藏生产开发参数，再分别结合物质平衡方程和Fetkovitch方程计算第一油藏水侵量和第二油藏水侵量，第一油藏水侵量和第二油藏水侵量均包括不同时刻的油藏水侵量的值，基于不同时刻的油藏水侵量的值对第一油藏水侵量和第二油藏水侵量的关系进行分析，当每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值时，根据第二油藏水侵量所使用的预设水体体积的值计算油藏水体倍数。本发明实施例计算过程简单，能够快速计算出油藏的水体倍数，可为油藏天然能量评价提供数据依据。
为了提高确定第一油藏水侵量和第二油藏水侵量关系的效率，若每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值，则根据第二油藏水侵量所使用的预设水体体积的值计算油藏水体倍数之前，还可以执行如下步骤：
(1)根据第一油藏水侵量和每个第二油藏水侵量分别生成多个油藏水侵量交会图；不同的油藏水侵量交会图对应的预设水体体积的值不同；
在本发明实施例中，参见图3所示的油藏水侵量交会图示意图，侵量交会图用于描述第一油藏水侵量和第二油藏水侵量之间的关系趋势。图中，横轴和纵轴分别表示第一油藏水侵量和第二油藏水侵量在不同时刻的水侵量值。由于根据不同的预设水体体积的值计算得到的第二油藏水侵量不同，因此，对于每一个预设水体体积的值得到的第二油藏水侵量均可以得到与第一油藏水侵量的水侵量交会图，水侵量交会图的个数与所使用的预设水体体积的值的个数相同。
(2)获取第一油藏水侵量和每个第二油藏水侵量的关系判断结果；关系判断结果基于油藏水侵量交会图生成；关系判断结果为第一油藏水侵量第二油藏水侵量相关或第一油藏水侵量第二油藏水侵量不相关。
在本发明实施例中，参见图3所示的油藏水侵量交会图示意图，如果某时刻第一油藏水侵量和每个第二油藏水侵量的值均为950，在油藏水侵量交会图中，(950，950)这个点将落在分别与横轴、纵轴呈45°的对角线上，如果多个时刻的第一油藏水侵量和每个第二油藏水侵量的值相等或几乎相等，则图3中落在对角线上的点或落在对角线附近的点会更多，因此，关系判断结果可以用于描述第一油藏水侵量和每个第二油藏水侵量在不同时刻值的大小关系。当关系判断结果为第一油藏水侵量第二油藏水侵量相关时，表明多个时刻的第一油藏水侵量和第二油藏水侵量的值相等或几乎相等；当关系判断结果为第一油藏水侵量第二油藏水侵量不相关时，表明多个时刻的第一油藏水侵量和第二油藏水侵量的值不相等或相差较大。基于油藏水侵量交会图可以得到关系判断结果。获取第一油藏水侵量和每个第二油藏水侵量之间的关系判断结果。
为了快速定量评价油藏水体大小，若每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值，则根据第二油藏水侵量所使用的预设水体体积的值计算油藏水体倍数，可以按照如下步骤执行：
当关系判断结果为第一油藏水侵量第二油藏水侵量相关时，将相应第二油藏水侵量对应的预设水体体积的值作为目标水体体积；根据目标水体体积和目标油藏水体计算公式计算油藏水体倍数。
在本发明实施例中，当第一油藏水侵量第二油藏水侵量相关时，获取计算该第二油藏水侵量时所使用的预设水体体积的值，并将该预设水体体积的值作为目标水体体积，将目标水体体积作为目标油藏水体计算公式的一个计算参数，从而得到油藏水体倍数。目标油藏水体计算公式可以根据实际需求选择，本发明实施例对此不作具体限定。
例如，根据目标水体体积和目标油藏水体计算公式计算油藏水体倍数，可以包括：根据如下公式计算油藏水体倍数：
Figure BDA0002176900000000071
其中，k为油藏水体倍数，Vw为目标水体体积，Swc为束缚水饱和度，N为油藏地质储量，Boi为原始条件下的原油体积系数。
根据油藏地质特征参数、油藏生产开发参数和物质平衡方程计算第一油藏水侵量，包括：根据如下公式计算第一油藏水侵量：
Figure BDA0002176900000000072
其中，Np为标准状态下累计产油量，单位为m3；N为油藏地质储量，单位为m3；Bo为目前地层压力下的原油体积系数，单位为m3/m3；Bg为目前地层压力下的气体体积系数，单位为m3/m3；Boi为原始条件下的原油体积系数，单位为m3/m3；Rp为生产气油比，单位为m3/m3；Rs为目前地层压力下原油的溶解气油比，单位为m3/m3；Rsi为原始条件下原油的溶解气油比，单位为m3/m3；Cf为孔隙压缩系数，单位为MPa-1；Cw为地层水压缩系数，单位为MPa-1；Bw为目前地层压力下水的体积系数，单位为m3/m3；Swc为束缚水饱和度，为分数；Δp为地层压力降落差，单位为MPa；We为水侵量，单位为m3；Wi为累计注水量，单位为m3；Wp为累计产水量，单位为m3。
在本发明实施例中，参见图2所示的边水油藏边缘注水开发示意图，根据物质平衡原理，边水油藏的地面累计产量转换到地层条件下，应等于油藏中因地层压力下降所引起流体膨胀量、注入流体量和外部水侵量之和，即有：
Figure BDA0002176900000000073
根据该式可得到：
Figure BDA0002176900000000074
根据油藏地质特征参数及油藏生产开发参数以及上述物质平衡方程即可得到第一油藏水侵量。其中，油藏地质特征参数包括油藏地质储量、孔隙压缩系数以及地层水压缩系数等，油藏生产开发参数包括累计注水量以及生产气油比等多种参数。
根据油藏地质特征参数、油藏生产开发参数、预设水体体积和目标水侵量方程计算第二油藏水侵量，包括：根据如下公式计算第二油藏水侵量：
Figure BDA0002176900000000081
Wei＝Vw(Cw+Cf)pwi；
Figure BDA00021769000000000810
其中，Wei为水体的最大水侵潜量；Vw为水体的体积，单位为m3；Cf为孔隙压缩系数，单位为MPa-1；Cw为地层水压缩系数，单位为MPa-1；pwi为水体的原始地层压力，单位为MPa；J为水侵指数，单位为m3/(d·MPa)；
Figure BDA0002176900000000082
为油藏的圆周系数，
Figure BDA0002176900000000083
θ为油藏的水侵角，单位为度；k为储集层渗透率，单位为μm2；h为水层厚度，单位为m；a为换算系数，a＝86.4；μw为地层水的粘度，单位为mPa·s；re、ro分别为水体和油藏的半径，单位为m；t为开发时刻，单位为分；p为目前条件下油藏的平均压力，单位为MPa。
考虑到油藏开发过程中，油藏的实际平均压力是不断变化的，可以做离散化处理，因此，根据油藏地质特征参数、油藏生产开发参数、预设水体体积和目标水侵量方程计算第二油藏水侵量，包括：根据如下公式计算第二油藏水侵量：
Figure BDA0002176900000000084
Wei＝Vw(Cw+Cf)pwi；
Figure BDA0002176900000000089
Figure BDA0002176900000000085
其中，Wei为水体的最大水侵潜量，单位为m3；Vw为水体的体积，单位为m3；Cf为孔隙压缩系数，单位为MPa-1；Cw为地层水压缩系数，单位为MPa-1；pwi为水体的原始地层压力，单位为MPa；J为水侵指数，单位为m3/(d·MPa)；
Figure BDA0002176900000000086
为油藏的圆周系数，
Figure BDA0002176900000000087
θ为油藏的水侵角，单位为度；k为储集层渗透率，单位为μm2；h为水层厚度，单位为m；a为换算系数，a＝86.4；μw为地层水的粘度，单位为mPa·s；re、ro分别为水体和油藏的半径，单位为m；t为开发时刻，单位为分；pn为时间段Δtn内油藏的平均压力，单位为MPa；
Figure BDA0002176900000000088
为第n时间段Δtn内的水体平均压力，单位为MPa。
本发明实施例提供的一种油藏水体倍数计算方法及装置，可以利用半解析方法快速计算出油藏的水体倍数，可为油藏天然能量评价提供依据。
某强边水油藏为一多层状背斜型砂岩油藏，油藏埋深1550-1950m，以粗砂岩和中砂岩为主，为缓坡快速堆积型扇三角洲沉积，平均油层厚度约11m，平均孔隙度为25.8％，平均渗透率为1320md，为高孔高渗储层。地层原油粘度0.456mPa·s，地层原油密度为0.676g/cm3，为低密度低粘度原油。油藏原始地层压力19.17MPa，目前平均地层压力为14.5MPa。原油初始溶解气油比为120.7m3/m3，原油初始体积系数为1.389。地层束缚水饱和度0.14，地层水等温压缩系数4.1×10-4MPa-1，孔隙压缩系数为4.9×10-4MPa-1。
油藏于1996年投入开发，2002年开始边外注水开发。油藏目前综合含水42.5％，地质储量采出程度66.3％，采油速度3.4％。基于油藏的生产数据及流体物性，利用上述计算模型计算该油藏的水体倍数。当边水水体体积为2.8×10-9m3时，物质平衡法和Fetkovitch法得到的水侵在交会图上大致呈450的直线关系，参见附图3。由此根据公式计算，可得到油藏的水体倍数约为35。在该水体倍数下，由物质平衡法和Fetkovitch法分别得到的累积水侵量参见4所示的油藏累积水侵量随开采时间的变化曲线。从图4可以看出，在开发早期，两种方法的水侵量计算结果差别较大；而在开发中后期，两种方法的结果基本吻合。由此可见，前面所建立的油藏水体倍数计算模型适用于开发中后期。
本发明实施例还提供一种油藏水体倍数计算装置，参见图5所示的一种油藏水体倍数计算装置结构框图，该装置包括：
获取模块51，用于获取油藏地质特征参数和油藏生产开发参数；第一水侵量计算模块52，用于根据油藏地质特征参数、油藏生产开发参数和物质平衡方程计算第一油藏水侵量；第一油藏水侵量包括不同时刻的油藏水侵量的值；第二水侵量计算模块53，用于根据油藏地质特征参数、油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量；第二油藏水侵量包括不同时刻的油藏水侵量的值；根据不同的预设水体体积生成的第二油藏水侵量不同；油藏水体倍数计算模块54，用于若每个相同时刻的第一油藏水侵量的值和第二油藏水侵量的值相同或差值小于预设值，则根据第二油藏水侵量所使用的预设水体体积的值计算油藏水体倍数。
参见图6所示的另一种油藏水体倍数计算装置结构框图，该装置还包括关系判断模块55，用于：根据第一油藏水侵量和每个第二油藏水侵量分别生成多个油藏水侵量交会图；不同的油藏水侵量交会图对应的预设水体体积的值不同；获取第一油藏水侵量和每个第二油藏水侵量的关系判断结果；关系判断结果基于油藏水侵量交会图生成；关系判断结果为第一油藏水侵量第二油藏水侵量相关或第一油藏水侵量第二油藏水侵量不相关。
本发明实施例还提供一种计算机设备，参见图7所示的计算机设备结构示意框图，该计算机设备包括存储器71、处理器72，存储器中存储有可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一种方法的步骤。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的计算机设备的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述
本发明实施例还提供一种具有处理器可执行的非易失的程序代码的计算机可读介质，程序代码使处理器执行上述任一种方法的步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种油藏水体倍数计算方法，其特征在于，包括：
获取油藏地质特征参数和油藏生产开发参数；
根据所述油藏地质特征参数、所述油藏生产开发参数和物质平衡方程计算第一油藏水侵量；所述第一油藏水侵量包括不同时刻的油藏水侵量的值；
根据所述油藏地质特征参数、所述油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量；所述第二油藏水侵量包括不同时刻的油藏水侵量的值；根据不同的所述预设水体体积生成的第二油藏水侵量不同；
若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数；
其中，根据所述油藏地质特征参数、所述油藏生产开发参数和物质平衡方程计算第一油藏水侵量，包括：
根据如下公式计算第一油藏水侵量：
Figure FDA0003786298510000011
其中，Np为标准状态下累计产油量；N为油藏地质储量；Bo为目前地层压力下的原油体积系数；Bg为目前地层压力下的气体体积系数；Boi为原始条件下的原油体积系数；Rp为生产气油比；Rs为目前地层压力下原油的溶解气油比；Rsi为原始条件下原油的溶解气油比；Cf为孔隙压缩系数；Cw为地层水压缩系数；Bw为目前地层压力下水的体积系数；Swc为束缚水饱和度；Δp为地层压力降落差；We为水侵量；Wi为累计注水量；Wp为累计产水量。
2.根据权利要求1所述的方法，其特征在于，若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数之前，还包括：
根据所述第一油藏水侵量和每个所述第二油藏水侵量分别生成多个油藏水侵量交会图；不同的油藏水侵量交会图对应的所述预设水体体积的值不同；
获取所述第一油藏水侵量和每个所述第二油藏水侵量的关系判断结果；所述关系判断结果基于所述油藏水侵量交会图生成；所述关系判断结果为第一油藏水侵量第二油藏水侵量相关或第一油藏水侵量第二油藏水侵量不相关。
3.根据权利要求2所述的方法，其特征在于，若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数，包括：
当所述关系判断结果为第一油藏水侵量第二油藏水侵量相关时，将相应第二油藏水侵量对应的预设水体体积的值作为目标水体体积；
根据所述目标水体体积和目标油藏水体计算公式计算油藏水体倍数。
4.根据权利要求3所述的方法，其特征在于，根据所述目标水体体积和目标油藏水体计算公式计算油藏水体倍数，包括：
根据如下公式计算油藏水体倍数：
Figure FDA0003786298510000021
其中，k为油藏水体倍数，Vw为目标水体体积，Swc为束缚水饱和度，N为油藏地质储量，Boi为原始条件下的原油体积系数。
5.根据权利要求1-4任一项所述的方法，其特征在于，根据所述油藏地质特征参数、所述油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量，包括：
根据如下公式计算第二油藏水侵量：
Figure FDA0003786298510000022
Wei＝Vw(Cw+Cf)pwi
Figure FDA0003786298510000023
其中，Wei为水体的最大水侵潜量；Vw为水体的体积；Cf为孔隙压缩系数；Cw为地层水压缩系数；pwi为水体的原始地层压力；J为水侵指数；
Figure FDA0003786298510000024
为油藏的圆周系数，
Figure FDA0003786298510000025
θ为油藏的水侵角；k为储集层渗透率；h为水层厚度；a为换算系数，a＝86.4；μw为地层水的粘度；re、ro分别为水体和油藏的半径；t为开发时刻；p为目前条件下油藏的平均压力。
6.根据权利要求1-4任一项所述的方法，其特征在于，根据所述油藏地质特征参数、所述油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量，包括：
根据如下公式计算第二油藏水侵量：
Figure FDA0003786298510000031
Wei＝Vw(Cw+Cf)pwi
Figure FDA0003786298510000032
Figure FDA0003786298510000033
Figure FDA0003786298510000034
其中，Wei为水体的最大水侵潜量；Vw为水体的体积；Cf为孔隙压缩系数；Cw为地层水压缩系数；pwi为水体的原始地层压力；J为水侵指数；
Figure FDA0003786298510000035
为油藏的圆周系数，
Figure FDA0003786298510000036
θ为油藏的水侵角；k为储集层渗透率；h为水层厚度；a为换算系数，a＝86.4；μw为地层水的粘度；re、ro分别为水体和油藏的半径；t为开发时刻；pn为时间段Δtn内油藏的平均压力；
Figure FDA0003786298510000037
为第n时间段Δtn内的水体平均压力。
7.一种油藏水体倍数计算装置，其特征在于，包括：
获取模块，用于获取油藏地质特征参数和油藏生产开发参数；
第一水侵量计算模块，用于根据所述油藏地质特征参数、所述油藏生产开发参数和物质平衡方程计算第一油藏水侵量；所述第一油藏水侵量包括不同时刻的油藏水侵量的值；
第二水侵量计算模块，用于根据所述油藏地质特征参数、所述油藏生产开发参数、预设水体体积和Fetkovitch方程计算第二油藏水侵量；所述第二油藏水侵量包括不同时刻的油藏水侵量的值；根据不同的所述预设水体体积生成的第二油藏水侵量不同；
油藏水体倍数计算模块，用于若每个相同时刻的所述第一油藏水侵量的值和所述第二油藏水侵量的值相同或差值小于预设值，则根据所述第二油藏水侵量所使用的所述预设水体体积的值计算油藏水体倍数；
其中，根据所述油藏地质特征参数、所述油藏生产开发参数和物质平衡方程计算第一油藏水侵量，包括：
根据如下公式计算第一油藏水侵量：
Figure FDA0003786298510000041
其中，Np为标准状态下累计产油量；N为油藏地质储量；Bo为目前地层压力下的原油体积系数；Bg为目前地层压力下的气体体积系数；Boi为原始条件下的原油体积系数；Rp为生产气油比；Rs为目前地层压力下原油的溶解气油比；Rsi为原始条件下原油的溶解气油比；Cf为孔隙压缩系数；Cw为地层水压缩系数；Bw为目前地层压力下水的体积系数；Swc为束缚水饱和度；Δp为地层压力降落差；We为水侵量；Wi为累计注水量；Wp为累计产水量。
8.根据权利要求7所述的装置，其特征在于，还包括关系判断模块，用于：
根据所述第一油藏水侵量和每个所述第二油藏水侵量分别生成多个油藏水侵量交会图；不同的油藏水侵量交会图对应的所述预设水体体积的值不同；
获取所述第一油藏水侵量和每个所述第二油藏水侵量的关系判断结果；所述关系判断结果基于所述油藏水侵量交会图生成；所述关系判断结果为第一油藏水侵量第二油藏水侵量相关或第一油藏水侵量第二油藏水侵量不相关。
9.一种计算机设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至6任一项所述的方法的步骤。
10.一种具有处理器可执行的非易失的程序代码的计算机可读介质，其特征在于，所述程序代码使所述处理器执行上述权利要求1至6任一项所述的方法。
CN201910782014.6A 2019-08-23 2019-08-23 油藏水体倍数计算方法及装置 Active CN112487594B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910782014.6A CN112487594B (zh) 2019-08-23 2019-08-23 油藏水体倍数计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910782014.6A CN112487594B (zh) 2019-08-23 2019-08-23 油藏水体倍数计算方法及装置

Publications (2)

Publication Number Publication Date
CN112487594A CN112487594A (zh) 2021-03-12
CN112487594B true CN112487594B (zh) 2022-11-01

Family

ID=74920297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910782014.6A Active CN112487594B (zh) 2019-08-23 2019-08-23 油藏水体倍数计算方法及装置

Country Status (1)

Country Link
CN (1) CN112487594B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113882850B (zh) * 2021-10-18 2023-05-09 中国石油大学(北京) 气藏动态储量预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858701A (zh) * 2019-02-12 2019-06-07 中国石油大学(北京) 裂缝性边水气藏水侵量的定量识别方法及系统
CN109918769A (zh) * 2019-03-04 2019-06-21 中国地质大学（武汉） 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法
CN109915122A (zh) * 2019-03-04 2019-06-21 中国地质大学（武汉） 利用全局方程计算缝洞型油藏稳态水侵水侵量的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109858701A (zh) * 2019-02-12 2019-06-07 中国石油大学(北京) 裂缝性边水气藏水侵量的定量识别方法及系统
CN109918769A (zh) * 2019-03-04 2019-06-21 中国地质大学（武汉） 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法
CN109915122A (zh) * 2019-03-04 2019-06-21 中国地质大学（武汉） 利用全局方程计算缝洞型油藏稳态水侵水侵量的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
西58-8小断块边水油藏天然能量评价研究;周春香等;《岩性油气藏》;20091215(第04期);全文 *
边底水油藏水体及储层参数定量评价;于清艳等;《油气地质与采收率》;20170525(第03期);全文 *

Also Published As

Publication number Publication date
CN112487594A (zh) 2021-03-12

Similar Documents

Publication Publication Date Title
CN106777651B (zh) 注采平衡原理的油水井产量劈分方法
CN105041307B (zh) 一种碎屑岩油气储层优势渗流通道识别的工业化流程
CN105569646B (zh) 一种油气井技术可采储量预测方法
Ehrenberg et al. A megascale view of reservoir quality in producing sandstones from the offshore Gulf of Mexico
CN104453834A (zh) 一种井组注采关系优化调整方法
CN101806224B (zh) 井下煤层气抽采产能预测方法
CN106600443A (zh) 动态含水饱和度的油井产量劈分方法
CN106894814A (zh) 复杂断块油藏高含水后期剩余油二次富集的快速识别方法
CN103455667B (zh) 充气法治理承压含水层海水入侵的数值模拟方法
CN105370269B (zh) 一种基于压实校正等时地层界面追踪对比方法
CN104481520A (zh) 页岩气井可采储量早期评价方法
CN104122602A (zh) 一种获取沉积盆地剥蚀厚度的方法及装置
CN104712328B (zh) 快速评价复杂油藏中单个流动单元动用状况的方法
CN104533519B (zh) 立井井筒通过强含水厚岩层时涌水水害的治理方法
CN111706318B (zh) 一种确定低渗储层剩余油分布状况的方法
CN110984974B (zh) 基于水侵速度确定有水气藏合理采气速度的方法
CN112487594B (zh) 油藏水体倍数计算方法及装置
CN115034489A (zh) 一种考虑溶解作用的气藏co2埋存潜力预测方法
CN115270662A (zh) 一种欠饱和储层煤层气井生产边界定量预测的计算方法
Al-Obaidi et al. The Effect of anisotropy in formation permeability on the efficiency of cyclic water flooding
CN110318741B (zh) 一种异常高压、高含水致密低渗气藏剩余气分布描述方法
CN109918769A (zh) 利用瞬时方程计算缝洞型油藏非稳态水侵水侵量的方法
CN104459820B (zh) 多尺度压力场的建立方法
RU2535545C1 (ru) Способ разработки нефтяной залежи
CN106529199B (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