CN113657046A - 油藏复杂结构井试井解释方法、装置、电子设备及介质 - Google Patents

油藏复杂结构井试井解释方法、装置、电子设备及介质 Download PDF

Info

Publication number
CN113657046A
CN113657046A CN202110919976.9A CN202110919976A CN113657046A CN 113657046 A CN113657046 A CN 113657046A CN 202110919976 A CN202110919976 A CN 202110919976A CN 113657046 A CN113657046 A CN 113657046A
Authority
CN
China
Prior art keywords
well
oil reservoir
heavy oil
complex structure
zonal
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
CN202110919976.9A
Other languages
English (en)
Other versions
CN113657046B (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 National Offshore Oil Corp Shenzhen Branch
Original Assignee
China National Offshore Oil Corp Shenzhen 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 National Offshore Oil Corp Shenzhen Branch filed Critical China National Offshore Oil Corp Shenzhen Branch
Priority to CN202110919976.9A priority Critical patent/CN113657046B/zh
Publication of CN113657046A publication Critical patent/CN113657046A/zh
Application granted granted Critical
Publication of CN113657046B publication Critical patent/CN113657046B/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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明实施例公开了一种油藏复杂结构井试井解释方法、装置、电子设备及介质。所述方法包括：确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；依据试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。采用本申请方案，在考虑出砂引起的分区物性变化及多种复杂渗流特征的综合作用来对这类疏松砂岩稠油油藏复杂结构井进行试井解释，能够适用于任意分区形态及复杂井型。

Description

油藏复杂结构井试井解释方法、装置、电子设备及介质
技术领域
本发明实施例涉及油气开发技术领域，尤其涉及一种油藏复杂结构井试井解释方法、装置、电子设备及介质。
背景技术
试井是对油气井进行专门的测试工作，是为了确定井的生产能力和研究储层的参数及动态情况，试井能提供可靠的地层动态信息。试井解释是利用渗流理论分析测试资料，评价地层或井下参数的方法，是油气渗流理论在油气田开发中的实际应用。但是，目前针对分区试井压力曲线的研究主要是在复合分区模型中考虑单个渗流特征，针对疏松砂岩稠油油藏存在的出砂及复杂渗流特征多重因素综合影响的问题，缺少相关试井理论分析，导致常规的试井解释方法得到的解释结果与实际储层特征往往不符。
发明内容
本发明实施例中提供了一种油藏复杂结构井试井解释方法、装置、电子设备及介质，以实现对疏松砂岩稠油油藏复杂结构井进行试井解释。
第一方面，本发明实施例中提供了一种油藏复杂结构井试井解释方法，所述方法包括：
确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
第二方面，本发明实施例中还提供了一种油藏复杂结构井试井解释装置，所述装置包括：
信息确定模块，用于确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
试井压力曲线确定模块，用于依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
流动阶段识别模块，用于依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
试井解释模块，用于分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
第三方面，本发明实施例中还提供了一种电子设备，包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如本发明任意实施例中提供的油藏复杂结构井试井解释方法。
第四方面，本发明实施例中还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明任意实施例中提供的油藏复杂结构井试井解释方法。
本发明实施例中提供了一种油藏复杂结构井试井解释方法，确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；分区渗流特征信息包括应力敏感和启动压力梯度；依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；依据试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。采用本申请方案，在考虑出砂引起的分区物性变化及多种复杂渗流特征(包括应力敏感、启动压力梯度)的综合作用来对这类疏松砂岩稠油油藏复杂结构井进行试井解释，能够适用于任意分区形态及复杂井型，克服了传统早期计算精度低及计算量大的缺点。
上述发明内容仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1是本发明实施例中提供的一种油藏复杂结构井试井解释方法的流程图；
图2是本发明实施例中提供的一种疏松砂岩稠油油藏出砂示意图；
图3是本发明实施例中提供的一种疏松砂岩稠油油藏的复合分区示意图；
图4是本发明实施例中提供的一种疏松砂岩稠油油藏的水平井试井压力曲线示意图；
图5是本发明实施例中提供的一种疏松砂岩稠油油藏的直井典型流动形态示意图；
图6是本发明实施例中提供的一种疏松砂岩稠油油藏的试井解释技术流程图；
图7是本发明实施例中提供的另一种油藏复杂结构井试井解释方法的流程图；
图8是本发明实施例中提供的一种压力数据格式的局部示意图；
图9是本发明实施例中提供的一种产量数据格式的局部示意图；
图10a是本发明实施例中提供的一种预处理后的压力及压力导数的局部示意图；
图10b是本发明实施例中提供的一种实际测压数据曲线示意图；
图11是本发明实施例中提供的一种利用疏松砂岩稠油油藏试井解释流程进行理论图版曲线示意图；
图12是本发明实施例中提供的一种测压数据拟合图；
图13是本发明实施例中提供的一种油藏复杂结构井试井解释装置的结构图；
图14是本发明实施例中提供的一种电子设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
在更加详细地讨论示例性实施例之前应当提到的是，一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各项操作(或步骤)描述成顺序的处理，但是其中的许多操作可以被并行地、并发地或者同时实施。此外，各项操作的顺序可以被重新安排。当其操作完成时所述处理可以被终止，但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。
图1是本发明实施例中提供的一种油藏复杂结构井试井解释方法的流程图，本实施例可适用于对疏松砂岩稠油油藏复杂结构井进行试井解释的情况。该方法可由油藏复杂结构井试井解释装置来执行，该装置可由软件和/或硬件实现，并集成在任何具有网络通信功能的电子设备上。如图1所示，本实施例中提供的油藏复杂结构井试井解释方法，可包括以下步骤：
S110、确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；分区渗流特征信息包括应力敏感和启动压力梯度。
参见图2，在疏松砂岩稠油油藏中，由于出砂引起渗透率变化的作用，从水平井筒到油藏边界可以划分为淤塞区、疏通区及不出区三个区域。淤塞区考虑为表皮因子，出砂会引起储层物性改善(K1)，出砂区渗透率大于储层初始渗透率，不出砂区为储层原始渗透率(K2)，形成疏松砂岩稠油油藏复合分区模型。
参见图3，对于疏松砂岩稠油油藏，疏松砂岩储层由出砂引起储层物性分区变化，将井筒附近淤塞区考虑为表皮因子(S)，从内到外分为疏通区和不出砂区，其中疏通区渗透率为K1，不出砂区渗透率为K2，并且由于疏通区的出砂导致疏通区渗透率为K1大于不出砂区渗透率为K2，即可得到疏松砂岩稠油油藏出砂引起的分区物性变化信息。可选地，分区物性变化信息包括疏松砂岩稠油油藏中由出砂引起渗透率变化作用产生的储层物性分区变化。
分区渗流特征信息包括应力敏感和启动压力梯度。应力敏感性可以是指多孔介质孔隙体积(孔隙度)、渗透率随有效应力变化的改变量；启动压力梯度可以是流体在低渗透油藏中渗流时必须有一个附加的压力梯度克服岩石表面吸附膜或水化膜引起的阻力才能流动，该附加压力梯度称为启动压力梯度。
S120、依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息。
针对疏松砂岩稠油油藏复杂结构井，在开展海上疏松砂岩稠油试井解释时，不仅考虑疏松砂岩稠油油藏中存在应力敏感与启动压力梯度，同时还考虑疏松砂岩稠油油藏出砂引起的典型圆形复合分区特征，针对海上疏松砂岩油藏生产过程出砂引起的渗透率变化，将油藏划分为近井出砂区及远井不出砂区，结合出砂引起的分区物性变化信息，从多个维度来分析关键参数对试井压力曲线的影响，得到考虑应力敏感、启动压力梯度以及复合分区的疏松砂岩稠油油藏复杂结构井的试井压力曲线信息。
S130、依据试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段。
参见图4，疏松砂岩稠油油藏复杂结构井的试井压力曲线可包括以下流动阶段：井储及表皮阶段、早期径向流阶段、第一过渡流阶段、内区中期线性流阶段、内区整体拟径向流阶段、内外区过渡流阶段以及外边界响应阶段。参见图4与图5，针对各个流动阶段的描述如下：
第Ⅰ阶段：井储及表皮阶段。在流动刚开始发生时，内区流体未动用，只有井筒内部储存的流体流动，生产井的产量完全由近井地区供给，此时表皮对产量的影响相对较大，内区几乎无动用。如图4中第Ⅰ阶段。
第Ⅱ阶段：早期径向流阶段。流动开始一段时间后，近井地区流体的流动可以看成平面径向流，如图4中第Ⅱ阶段。
第Ⅲ阶段：第一过渡流阶段。流动经过早期径向流之后，内区流体开始被动用，此时早期径向流和内区向近井地区流动同时存在。如图4中第Ⅲ阶段。
第Ⅳ阶段：内区中期线性流阶段。随着流动的进行，近井地区流体大多数被采出，此时内区流体的流动状态类似于一维直线渗流，如图4中第Ⅳ阶段。
第Ⅴ阶段：内区整体拟径向流阶段。此时内区流体大多数被动用，整合内区的流动可以看成平面径向流，形成了拟径向流，如图中第Ⅴ阶段。
第Ⅵ阶段：内外区过渡流阶段。流动经过内区整体拟径向流之后，外区流体开始被动用，此时内区整体拟径向流和外区向内区流动同时存在。如图4中第Ⅵ阶段。
第Ⅶ阶段：外边界响应阶段。流动进行到最后，外区流体大多数被采出，外边界开始受到影响，此时根据边界不同可以将其分为三类。第一类：地层为无限大地层，这种情况下可以看作边界对油藏一直存在供给，因此试井压力曲线保持正常；第二类：外边界处有边水，油井迅速见水；第三类：边界为断层，此时产量迅速降低，但是油井不见水，因此其试井压力曲线上升。
S140、分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
在本实施例的一种可选方案中，参见图6，分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果，可包括以下步骤A1-A3：
步骤A1、对疏松砂岩稠油油藏复杂结构井的油田生产数据进行整化处理。
步骤A2、依据整化处理的油田生产数据分阶段进行流动分析，执行疏松砂岩稠油油藏复杂结构井的储层参数反演初值。
步骤A3、依据全测试阶段的压力测试数据，全测压拟合调整疏松砂岩稠油油藏复杂结构井的储层参数，以得到油藏复杂结构井的试井解释结果。
整理疏松砂岩稠油油藏复杂结构井的实际油田生产数据，将其进行规范化处理，接着识别其各个流动阶段，然后分阶段进行流动处理，接着进行储层参数反演初值，然后加入其约束条件，根据全测历史拟合调整储层参数，最后判断其是否满足误差，如果满足误差，则得出DST试井解释结果。如果不满足误差，则对各个流动阶段重新进行分析，重复上述步骤，一直到全测历史拟合调整储层参数满足误差，此时得出疏松砂岩稠油油藏复杂结构井的试井解释结果。
可选地，在得到考虑应力敏感、启动压力梯度以及复合分区的疏松砂岩稠油油藏复杂结构井的试井压力曲线的基础上，针对海上疏松砂岩稠油油藏复杂结构井的测试数据的“因素多、段数多、时间短”特殊性，引入反褶积技术利用整化处理及筛选后的压力测试数据，结合流动阶段分析结果及全测试阶段的压力和产量拟合结果，形成了一套适用于疏松砂岩稠油油藏试井解释。反褶积方法可包括但不限于以下内容：脉冲反褶积、预测反褶积及最小平方法反褶积。
根据本发明实施例中提供的油藏复杂结构井试井解释方法，在考虑出砂引起的分区物性变化及多种复杂渗流特征(包括应力敏感、启动压力梯度)的综合作用来对这类疏松砂岩稠油油藏复杂结构井进行试井解释，特别是能有效处理具有“多、短”特点的海上测试压力数据，形成了一套适用于疏松砂岩稠油油藏试井解释流程，能够适用于任意分区形态及复杂井型，克服了传统早期计算精度低及计算量大的缺点。
图7是本发明实施例中提供的另一种油藏复杂结构井试井解释方法的流程图，本发明实施例在上述实施例的基础上对前述实施例进行进一步优化，本发明实施例可以与上述一个或者多个实施例中各个可选方案结合。如图7所示，本实施例中提供的油藏复杂结构井试井解释方法，可包括以下步骤：
S710、确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度。
S720、依据分区物性变化信息和分区渗流特征信息，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型。
在本实施例的一种可选方案中，依据分区物性变化信息和分区渗流特征信息，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型，可包括以下步骤B1-B4：
步骤B1、将疏松砂岩稠油油藏储层的应力敏感变量代入运动方程，并通过摄动法处理得到应力敏感指标下的线性方程。
对于疏松砂岩稠油油藏的储层应力敏感，可以利用指数式应力敏感进行表征，对于应力敏感引起的非线性问题，经过变量代换之后，可以运用摄动法进行处理得到应力敏感指标下的线性方程，具体涉及以下方式与方程。
K＝Kie-γΔP (1-1)
Figure BDA0003207098870000101
其中，非线性项的处理过程具体如下：
引入中间变量：
Figure BDA0003207098870000102
消除二次方项：
Figure BDA0003207098870000103
引入摄动原理：
Figure BDA0003207098870000104
线性模型：
Figure BDA0003207098870000105
在上述公式中，K为渗透率，mD；Ki为初始渗透率，mD；γ为应力敏感系数，无量纲；ΔP为生产压差，MPa；γD为应力敏感系数，无量纲；pD为无因此压力，无量纲；tD为无因次时间，无量纲。
步骤B2、将疏松砂岩稠油油藏储层的启动压力梯度无因次化，并将无因次化的启动压力梯度方程带入运动方程，得到启动压力梯度指标下的线性方程。
基于稠油流变特征，考虑启动压力梯度(等效为多孔介质中的达西流动)，并将启动压力梯度无因次化，然后将无因次化的启动压力梯度方程带入运动方程，再将运动方程无因次化，具体过程如下：
启动压力梯度：
Figure BDA0003207098870000111
启动压力梯度无因次化：
Figure BDA0003207098870000112
处理办法是先将启动压力梯度代入运动方程：
Figure BDA0003207098870000113
将运动方程无因次化：
Figure BDA0003207098870000114
在上述公式中：v为渗流速度，cm/s；λ为启动压力梯度，MPa/m；λD为无因次启动压力梯度，无量纲；μ为原油粘度，mPa·s；rD为无因次半径，无量纲；Kv为垂向渗透率，mD；Kh为径向渗透率，mD；hD为无因次油藏厚度，无量纲；LD为无因次水平段长度，无量纲；Ct为总压缩系数，MPa-1；p为地层压力，MPa；pi为原始地层压力MPa。r为半径，m；rw为井筒半径，m。
步骤B3、确定疏松砂岩稠油油藏出砂引起的分区物性变化指标下的分区渗流方程。
参见图3，疏松砂岩稠油油藏储层由出砂引起储层物性分区变化，将井筒附近淤塞区考虑为表皮因子(S)，从内到外分为疏通区和不出砂区，其中疏通区渗透率为K1，不出砂区渗透率为K2，并且由于疏通区的出砂导致疏通区渗透率为K1大于不出砂区渗透率为K2，可以推出分区渗流方程如下：
内区：
Figure BDA0003207098870000121
外区：
Figure BDA0003207098870000122
分区界面连接条件如下：
压力连续：
Figure BDA0003207098870000123
流量连续：
Figure BDA0003207098870000124
在上述公式中：λD为无因次启动压力梯度，无量纲。μ为原油粘度，mPa·s；rD为无因次半径，无量纲；Kv为垂向渗透率，mD；Kh为径向渗透率，mD；hD为无因次油藏厚度，无量纲；LD为无因次水平段长度，无量纲；Ct为总压缩系数，MPa-1；r为半径，m；rw为井筒半径，m；M12是渗透率比，无量纲。
步骤B4、基于应力敏感指标下的线性方程、启动压力梯度指标下的线性方程以及分区物性变化指标下的分区渗流方程，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型。
通过拉普拉斯变换及正交变换之后，考虑出砂分区及应力敏感、启动压力梯度的综合渗流微分方程，即可得到疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型如下。
Figure BDA0003207098870000131
在上述模型中：M12为渗透率比，无量纲；s为拉式因子，无量纲；λ为启动压力梯度，MPa/m；λD为无因次启动压力梯度，无量纲。LD为无因次水平段长度，无量纲；zD为无因次垂向距离，无量纲；Z为z方向坐标，m；rD为无因次半径，无量纲。
S730、基于分区边界元求解方法对无因次试井解释数学模型进行求解，得到同时考虑应力敏感、启动压力梯度以及物性变化的试井压力曲线信息。
S740、依据试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段。
针对以上考虑出砂分区及应力敏感、启动压力梯度的综合渗流微分方程的求解如下：
Figure BDA0003207098870000132
考虑井储和表皮：
Figure BDA0003207098870000133
其中：
Figure BDA0003207098870000141
Figure BDA0003207098870000142
在上述公式中：M12为渗透率比，无量纲；s为拉式因子，无量纲；λ为启动压力梯度，MPa/m；λD为无因次启动压力梯度，无量纲。LD为无因次水平段长度，无量纲；rD为无因次半径，无量纲；K0(x)为零阶第二类修正贝塞尔函数；I0(x)为零阶第一类修正贝塞尔函数；K1(x)为一阶第二类修正贝塞尔函数；I1(x)为一阶第一类修正贝塞尔函数。
S750、分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
在上述实施例的基础上，以某疏松砂岩稠油油藏的现场测试压力数据为基础，结合该井的油田物性参数进行分析，具体流程及步骤如下：
步骤C1、原始数据预处理：按时间序列整理测压数据，然后将实际日期时间转换为为数据点格式，产量数据转换为持续数据格式(生产井为正，注水井为负)，压力数据转换为连续数据格式，如图8与图9所示。以及，利用sufer前处理软件将对准压力及产量的时间点，并进行降噪及预处理，得到压力及压力导数，如图10a与图10b所示。
步骤C2、理论曲线绘制：设置包括无因次井储系数、表皮系数、分区渗透率比、无因次分区半径、应力敏感系数及启动压力梯度等关键参数，利用疏松砂岩稠油油藏试井解释流程进行理论图版计算，如图11所示。
步骤C3、测压数据拟合：调整表皮系数、内区渗透率、分区半径三个值、应力敏感系数及启动压力梯度可根据渗透率参数计算得出，手动拟合实际压力及压力导数数据，表皮系数初值3-5、内区渗透率是外区渗透率的2倍、分区半径初值设置20m左右，其结果如图12所示。
根据本发明实施例中提供的油藏复杂结构井试井解释方法，在考虑出砂引起的分区物性变化及多种复杂渗流特征(包括应力敏感、启动压力梯度)的综合作用来对这类疏松砂岩稠油油藏复杂结构井进行试井解释，特别是能有效处理具有“多、短”特点的海上测试压力数据，形成了一套适用于疏松砂岩稠油油藏试井解释流程，建立了非牛顿流体下考虑应力敏感、启动压力梯度的出砂分区拉式空间渗流数学模型，形成的离散井筒-分区边界元混合高效求解方法能够适用于任意分区形态及复杂井型(不同井型)，能够适用于任意分区形态及复杂井型，克服了传统早期计算精度低及计算量大的缺点。
图13是本发明实施例中提供的一种油藏复杂结构井试井解释装置的结构图，本实施例的技术方案可适用于对疏松砂岩稠油油藏复杂结构井进行试井解释的情况。该装置可由软件和/或硬件实现，并集成在任何具有网络通信功能的电子设备上。如图13所示，本实施例中提供的油藏复杂结构井试井解释装置可包括以下：信息确定模块1310、试井压力曲线确定模块1320、流动阶段识别模块1330和试井解释模块1340。其中：
信息确定模块1310，用于确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
试井压力曲线确定模块1320，用于依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
流动阶段识别模块1330，用于依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
试井解释模块1340，用于分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
在上述实施例的基础上，可选地，所述分区物性变化信息包括疏松砂岩稠油油藏中由出砂引起渗透率变化作用产生的储层物性分区变化。
在上述实施例的基础上，可选地，依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息包括：
依据分区物性变化信息和分区渗流特征信息，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型；
基于分区边界元求解方法对无因次试井解释数学模型进行求解，得到同时考虑应力敏感、启动压力梯度以及物性变化的试井压力曲线信息。
在上述实施例的基础上，可选地，依据分区物性变化信息和分区渗流特征信息，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型，包括：
将疏松砂岩稠油油藏储层的应力敏感变量代入运动方程，并通过摄动法处理得到应力敏感指标下的线性方程；
将疏松砂岩稠油油藏储层的启动压力梯度无因次化，并将无因次化的启动压力梯度方程带入运动方程，得到启动压力梯度指标下的线性方程；
确定疏松砂岩稠油油藏出砂引起的分区物性变化指标下的分区渗流方程；
基于应力敏感指标下的线性方程、启动压力梯度指标下的线性方程以及分区物性变化指标下的分区渗流方程，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型。
在上述实施例的基础上，可选地，疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段包括：井储及表皮阶段、早期径向流阶段、第一过渡流阶段、内区中期线性流阶段、内区整体拟径向流阶段、内外区过渡流阶段以及外边界响应阶段。
在上述实施例的基础上，可选地，分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果，包括：
对疏松砂岩稠油油藏复杂结构井的油田生产数据进行整化处理
依据整化处理的油田生产数据分阶段进行流动分析，执行疏松砂岩稠油油藏复杂结构井的储层参数反演初值；
依据全测试阶段的压力测试数据，全测压拟合调整疏松砂岩稠油油藏复杂结构井的储层参数，以得到油藏复杂结构井的试井解释结果。
本发明实施例中所提供的油藏复杂结构井试井解释装置可执行上述本发明任意实施例中所提供的油藏复杂结构井试井解释方法，具备执行该油藏复杂结构井试井解释方法相应的功能和有益效果，详细过程参见前述实施例中的油藏复杂结构井试井解释方法的相关操作。
图14是本发明实施例中提供的一种电子设备的结构示意图。如图14所示结构，本发明实施例中提供的电子设备包括：一个或多个处理器1410和存储装置1420；该电子设备中的处理器1410可以是一个或多个，图14中以一个处理器1410为例；存储装置1420用于存储一个或多个程序；所述一个或多个程序被所述一个或多个处理器1410执行，使得所述一个或多个处理器1410实现如本发明实施例中任一项所述的油藏复杂结构井试井解释方法。
该电子设备还可以包括：输入装置1430和输出装置1440。
该电子设备中的处理器1410、存储装置1420、输入装置1430和输出装置1440可以通过总线或其他方式连接，图14中以通过总线连接为例。
该电子设备中的存储装置1420作为一种计算机可读存储介质，可用于存储一个或多个程序，所述程序可以是软件程序、计算机可执行程序以及模块，如本发明实施例中所提供的油藏复杂结构井试井解释方法对应的程序指令/模块。处理器1410通过运行存储在存储装置1420中的软件程序、指令以及模块，从而执行电子设备的各种功能应用以及数据处理，即实现上述方法实施例中油藏复杂结构井试井解释方法。
存储装置1420可包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据电子设备的使用所创建的数据等。此外，存储装置1420可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储装置1420可进一步包括相对于处理器1410远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
输入装置1430可用于接收输入的数字或字符信息，以及产生与电子设备的用户设置以及功能控制有关的键信号输入。输出装置1440可包括显示屏等显示设备。
并且，当上述电子设备所包括一个或者多个程序被所述一个或者多个处理器1410执行时，程序进行如下操作：
确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
当然，本领域技术人员可以理解，当上述电子设备所包括一个或者多个程序被所述一个或者多个处理器1410执行时，程序还可以进行本发明任意实施例中所提供的油藏复杂结构井试井解释方法中的相关操作。
本发明实施例中提供了一种计算机可读介质，其上存储有计算机程序，该程序被处理器执行时用于执行油藏复杂结构井试井解释方法，该方法包括：
确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
可选的，该程序被处理器执行时还可以用于执行本发明任意实施例中所提供的油藏复杂结构井试井解释定方法。
本发明实施例的计算机存储介质，可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(Random AccessMemory，RAM)、只读存储器(Read Only Memory，ROM)、可擦式可编程只读存储器(ErasableProgrammable Read Only Memory，EPROM)、闪存、光纤、便携式CD-ROM、光存储器件、磁存储器件、或者上述的任意合适的组合。计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于：电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：无线、电线、光缆、无线电频率(RadioFrequency，RF)等等，或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++，还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)——连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (10)

1.一种油藏复杂结构井试井解释方法，其特征在于，所述方法包括：
确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
2.根据权利要求1所述的方法，其特征在于，所述分区物性变化信息包括疏松砂岩稠油油藏中由出砂引起渗透率变化作用产生的储层物性分区变化。
3.根据权利要求1所述的方法，其特征在于，依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息包括：
依据分区物性变化信息和分区渗流特征信息，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型；
基于分区边界元求解方法对无因次试井解释数学模型进行求解，得到同时考虑应力敏感、启动压力梯度以及物性变化的试井压力曲线信息。
4.根据权利要求3所述的方法，其特征在于，依据分区物性变化信息和分区渗流特征信息，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型，包括：
将疏松砂岩稠油油藏储层的应力敏感变量代入运动方程，并通过摄动法处理得到应力敏感指标下的线性方程；
将疏松砂岩稠油油藏储层的启动压力梯度无因次化，并将无因次化的启动压力梯度方程带入运动方程，得到启动压力梯度指标下的线性方程；
确定疏松砂岩稠油油藏出砂引起的分区物性变化指标下的分区渗流方程；
基于应力敏感指标下的线性方程、启动压力梯度指标下的线性方程以及分区物性变化指标下的分区渗流方程，建立疏松砂岩稠油油藏复杂结构井的无因次试井解释数学模型。
5.根据权利要求1所述的方法，其特征在于，疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段包括：井储及表皮阶段、早期径向流阶段、第一过渡流阶段、内区中期线性流阶段、内区整体拟径向流阶段、内外区过渡流阶段以及外边界响应阶段。
6.根据权利要求1所述的方法，其特征在于，分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果，包括：
对疏松砂岩稠油油藏复杂结构井的油田生产数据进行整化处理
依据整化处理的油田生产数据分阶段进行流动分析，执行疏松砂岩稠油油藏复杂结构井的储层参数反演初值；
依据全测试阶段的压力测试数据，全测压拟合调整疏松砂岩稠油油藏复杂结构井的储层参数，以得到油藏复杂结构井的试井解释结果。
7.一种油藏复杂结构井试井解释装置，其特征在于，所述装置包括：
信息确定模块，用于确定疏松砂岩稠油油藏出砂引起的分区物性变化信息以及存在的分区渗流特征信息；所述分区渗流特征信息包括应力敏感和启动压力梯度；
试井压力曲线确定模块，用于依据分区物性变化信息和分区渗流特征信息，确定疏松砂岩稠油油藏复杂结构井的试井压力曲线信息；
流动阶段识别模块，用于依据所述试井压力曲线信息识别疏松砂岩稠油油藏复杂结构井试井所需的各个流动阶段；
试井解释模块，用于分阶段进行流动分析，依据分阶段的流动分析结果与全测试阶段的压力测试数据，确定油藏复杂结构井的试井解释结果。
8.根据权利要求7所述的装置，其特征在于，所述分区物性变化信息包括疏松砂岩稠油油藏中由出砂引起渗透率变化作用产生的储层物性分区变化。
9.一种电子设备，其特征在于，包括：
一个或多个处理器；
存储装置，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现权利要求1-6中任一所述的油藏复杂结构井试井解释方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1-6中任一所述的油藏复杂结构井试井解释方法。
CN202110919976.9A 2021-08-11 2021-08-11 油藏复杂结构井试井解释方法、装置、电子设备及介质 Active CN113657046B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110919976.9A CN113657046B (zh) 2021-08-11 2021-08-11 油藏复杂结构井试井解释方法、装置、电子设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110919976.9A CN113657046B (zh) 2021-08-11 2021-08-11 油藏复杂结构井试井解释方法、装置、电子设备及介质

Publications (2)

Publication Number Publication Date
CN113657046A true CN113657046A (zh) 2021-11-16
CN113657046B CN113657046B (zh) 2023-07-07

Family

ID=78491435

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110919976.9A Active CN113657046B (zh) 2021-08-11 2021-08-11 油藏复杂结构井试井解释方法、装置、电子设备及介质

Country Status (1)

Country Link
CN (1) CN113657046B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107133452A (zh) * 2017-04-18 2017-09-05 中国石油大学(北京) 油藏渗流数值模拟方法及装置
CN108019196A (zh) * 2017-11-22 2018-05-11 中国石油大学(北京) 一种体积压裂水平井试井解释处理方法及装置
AU2020101364A4 (en) * 2020-03-26 2020-08-20 University Of Science And Technology Beijing Simulation method for flow field of multi-stage fracturing on horizontal well in tight oil reservoir
CN111764890A (zh) * 2020-06-08 2020-10-13 中海石油(中国)有限公司深圳分公司 一种特低渗储层的试井分析方法
US10815778B1 (en) * 2019-12-02 2020-10-27 China University Of Petroleum (East China) Method for obtaining formation parameters of gas hydrate reservoir through well testing interpretation
CN112412411A (zh) * 2020-11-20 2021-02-26 中海石油(中国)有限公司 压裂井多井系统试井模型的建立及分析方法和装置
CN113111582A (zh) * 2021-04-12 2021-07-13 中国石油大学(北京) 基于深度学习的页岩储层试井智能解释分析方法及装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107133452A (zh) * 2017-04-18 2017-09-05 中国石油大学(北京) 油藏渗流数值模拟方法及装置
CN108019196A (zh) * 2017-11-22 2018-05-11 中国石油大学(北京) 一种体积压裂水平井试井解释处理方法及装置
US10815778B1 (en) * 2019-12-02 2020-10-27 China University Of Petroleum (East China) Method for obtaining formation parameters of gas hydrate reservoir through well testing interpretation
AU2020101364A4 (en) * 2020-03-26 2020-08-20 University Of Science And Technology Beijing Simulation method for flow field of multi-stage fracturing on horizontal well in tight oil reservoir
CN111764890A (zh) * 2020-06-08 2020-10-13 中海石油(中国)有限公司深圳分公司 一种特低渗储层的试井分析方法
CN112412411A (zh) * 2020-11-20 2021-02-26 中海石油(中国)有限公司 压裂井多井系统试井模型的建立及分析方法和装置
CN113111582A (zh) * 2021-04-12 2021-07-13 中国石油大学(北京) 基于深度学习的页岩储层试井智能解释分析方法及装置

Also Published As

Publication number Publication date
CN113657046B (zh) 2023-07-07

Similar Documents

Publication Publication Date Title
CN109838230B (zh) 油藏水淹层的定量评价方法
CN103775057B (zh) 一种致密油气藏有效储层的判识方法及装置
US9898559B2 (en) Method for predicting changes in properties of a formation in a near-wellbore area exposed to a drilling mud
US11573348B1 (en) Method and system using nuclear magnetic resonance well logging for T2 cutoff value estimation
US8301428B2 (en) Modeling the nonlinear hysteresis response of reservoir media
US10620340B2 (en) Tuning digital core analysis to laboratory results
CN105386751A (zh) 一种基于油藏渗流模型的水平井测井产能预测方法
RU2385413C1 (ru) Способ определения текущей газонасыщенности в призабойной зоне скважины в залежи летучей нефти
US11614417B2 (en) Determining saturation in low resistivity pay zones
CN110552694B (zh) 一种考虑多因素影响的泥质白云岩油藏油井产能评价方法
Watson et al. Identifiability of estimates of two-phase reservoir properties in history matching
CN106501146B (zh) 一种致密油藏物性上限确定方法及装置
CN111577264A (zh) 裂缝孔隙型油藏水平井产能预测方法及装置
CN112257349B (zh) 一种判断致密砂岩可动水气藏气井是否具有开发价值的方法
CN113657046B (zh) 油藏复杂结构井试井解释方法、装置、电子设备及介质
CN114293961B (zh) 海上压裂选层选段方法、装置、电子设备及存储介质
CN111894570A (zh) 一种开发中后期油藏老井水淹程度定量评价方法
Parra et al. Application of the Producer-Based Capacitance Resistance Model to Undersaturated Oil Reservoirs in Primary Recovery
CN107194067A (zh) 一种低渗透油藏数值模拟中启动压力梯度等效表征方法
CN116432366A (zh) 一种各向异性渗透率预测模型建立方法、装置及其应用
Zhao et al. Pore-scale hydraulic properties of virtual sandstone microstructures: spatial variations and voxel scale effects
CN109116440B (zh) 致密灰岩储层裂缝识别方法
CN113051518B (zh) 一种气油比定量计算方法、装置、计算机设备及存储介质
Schutjens et al. On the wellbore stress change caused by drawdown and depletion: an analytical model for a vertical well in a thin reservoir
CN113945992B (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