CN113887110A - 油井产能预测方法、装置、电子设备及计算机存储介质 - Google Patents

油井产能预测方法、装置、电子设备及计算机存储介质 Download PDF

Info

Publication number
CN113887110A
CN113887110A CN202111219918.1A CN202111219918A CN113887110A CN 113887110 A CN113887110 A CN 113887110A CN 202111219918 A CN202111219918 A CN 202111219918A CN 113887110 A CN113887110 A CN 113887110A
Authority
CN
China
Prior art keywords
reservoir
finite element
determining
element integral
position point
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.)
Pending
Application number
CN202111219918.1A
Other languages
English (en)
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 Oilfield Services Ltd
Original Assignee
China Oilfield Services 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 China Oilfield Services Ltd filed Critical China Oilfield Services Ltd
Priority to CN202111219918.1A priority Critical patent/CN113887110A/zh
Publication of CN113887110A publication Critical patent/CN113887110A/zh
Pending legal-status Critical Current

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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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"

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Development Economics (AREA)
  • Geometry (AREA)
  • Game Theory and Decision Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种油井产能预测方法、装置、电子设备及计算机存储介质。其中，所述方法包括：根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型；对渗流数学模型进行无因次化分析，确定污染区和非污染区每个位置点的无因次压力；根据每个位置点的无因次压力对渗流数学模型进行推导，确定有限元积分方程组；根据有限元积分方程组预测油田开发过程中储层受到不同污染程度时的弹性产量。本方案通过对污染区和非污染区的渗流数学模型进行无因次化分析，确定有限元积分方程组，能够准确计算在储层污染区边界不规则情况下的油井产能，更贴近地层的实际情况，客观反映污染区和未污染区之间的关联。

Description

油井产能预测方法、装置、电子设备及计算机存储介质
技术领域
本发明涉及油井产能预测技术领域，具体涉及一种油井产能预测方法、装置、电子设备及计算机存储介质。
背景技术
油井动态分析主要包括不稳定压力响应(TPR)和产量递减分析(RDA)，已有30多年的研究历史。学者提出了多种水平井模型来绘制典型曲线。同时，还推导出了大量可用于现场应用的简化公式。1973年，Gringarten和Ramey将源函数和格林函数引入了石油文献，该油井产能预测方法在解决水平井等非稳定流动问题方面具有很强的优越性。
之后，为了分析水平井的性动态分析，开展了多项研究工作。在传统方法中，水平井通常被简化为线源；然后，Ezulike和Igbokoyi研究了有限半径作为杆源的方法，使其更适合于水平井的动态分析。Ozkan and Raghavan提出了一种新的源函数方法来研究双区复合油藏水平井的不稳定压力特性，该方法一直是动态分析领域的主要求解方法之一。
国外方面，Chen等给出了复合油藏污染区、非污染区的点源解和线源解，在此基础上给出了复合油藏水平井压力动态。Ambastha等提出了利用矩形复合油藏模拟稠油注蒸汽热采的模型，采用数值方法分析了水平井完全在矩形复合油藏污染区时的压力动态，论证了当内非污染区存储比和流度比很大时会出现反映污染区边界的拟稳态阶段。
国内方面，王晓冬等运用汇源叠加的方法给出了不同边界条件下水平井完全在复合油藏污染区时的压力动态特征并绘制了典型曲线。石国新等运用拉普拉斯变换、分离变量法及特征函数和特征值法对双区复合油藏不稳定试井问题进行了研究，给出了不同边界条件下的井底压力解，并对实例进行了分析。
但是，上述方法均不能很好的反映油田开发的实际产量。
发明内容
鉴于上述问题，提出了本发明以便提供一种克服上述问题或者至少部分地解决上述问题的油井产能预测方法、装置、电子设备及计算机存储介质。
根据本发明的一个方面，提供了一种油井产能预测方法，包括：
根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型；
对所述渗流数学模型进行无因次化分析，确定所述污染区和非污染区每个位置点的无因次压力；
根据所述每个位置点的无因次压力对所述渗流数学模型进行推导，确定有限元积分方程组；
根据所述有限元积分方程组预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
根据本发明的另一个方面，提供了一种油井产能预测装置，包括：
渗流数学模型确定模块，用于根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型；
无因次化分析模块，用于对所述渗流数学模型进行无因次化分析，确定所述污染区和非污染区每个位置点的无因次压力；
方程推导模块，用于根据所述每个位置点的无因次压力对所述渗流数学模型进行推导，确定有限元积分方程组；
产量预测模块，用于根据所述有限元积分方程组预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
根据本发明的另一方面，提供了一种电子设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行本发明所述的油井产能预测方法对应的操作。
根据本发明的另一方面，提供了一种计算机存储介质，所述存储介质中存储有至少一可执行指令，所述可执行指令使处理器执行本发明所述的油井产能预测方法对应的操作。
根据本发明公开的油井产能预测方法、电子设备及计算机存储介质，通过对污染区和非污染区的渗流数学模型进行无因次化分析，确定有限元积分方程组，根据有限元积分方程组预测油田开发过程中储层受到不同污染程度时的弹性产量，能够准确计算在储层污染区边界不规则情况下的油井产能，更贴近地层的实际情况，客观反映污染区和未污染区之间的关联。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1示出了根据本发明实施例一提供的一种油井产能预测方法的流程示意图；
图2示出了根据本发明实施例一提供的一种油井产能预测方法中直井污染区域示意图；
图3示出了根据本发明实施例一提供的一种油井产能预测方法中水平井完全在污染区示意图；
图4示出了根据本发明实施例一提供的一种油井产能预测方法中水平井穿过污染区示意图；
图5示出了根据本发明实施例二提供的一种油井产能预测方法的流程示意图；
图6示出了根据本发明实施例二提供的一种油井产能预测方法中直井在污染区情况下的径向流流动状态的主视图；
图7示出了根据本发明实施例二提供的一种油井产能预测方法中直井在污染区情况下的径向流流动状态的俯视图；
图8示出了根据本发明实施例二提供的一种油井产能预测方法中直井在非污染区情况下的径向流流动状态的俯视图；
图9示出了根据本发明实施例二提供的一种油井产能预测方法中直井污染区半径的影响示意图；
图10示出了根据本发明实施例二提供的一种油井产能预测方法中直井污染区渗透率的影响示意图；
图11示出了根据本发明实施三提供的一种油井产能预测方法的流程示意图；
图12示出了根据本发明实施例三提供的一种油井产能预测方法中水平井完全在污染区情况下的径向流流动状态的主视图；
图13示出了根据本发明实施例三提供的一种油井产能预测方法中水平井完全在污染区情况下的径向流流动状态的俯视图；
图14示出了根据本发明实施例三提供的一种油井产能预测方法中水平井完全在污染区情况下的径向流流动状态的俯视图；
图15示出了根据本发明实施例三提供的一种油井产能预测方法中水平井完全在污染区情况下污染区半径的影响示意图；
图16示出了根据本发明实施例三提供的一种油井产能预测方法中水平井完全在污染区情况下污染区渗透率的影响示意图；
图17示出了根据本发明实施四提供的一种油井产能预测方法的流程示意图。
图18示出了根据本发明实施例四提供的一种油井产能预测方法中水平井穿过污染区情况下的径向流流动状态的主视图；
图19示出了根据本发明实施例四提供的一种油井产能预测方法中水平井穿过污染区情况下的径向流流动状态的俯视图；
图20示出了根据本发明实施例四提供的一种油井产能预测方法中水平井穿过污染区情况下污染区半径的影响示意图；
图21示出了根据本发明实施例四提供的一种油井产能预测方法中水平井穿过污染区情况下污染区渗透率的影响示意图；
图22示出了根据本发明实施例五提供的一种油井产能预测装置的结构示意图；
图23示出了根据本发明实施例七提供的一种电子设备的结构示意图。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
实施例一
图1示出了根据本发明实施例一提供的一种油井产能预测方法的流程示意图。如图1所示，该方法包括：
步骤S11，根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型。
其中，污染程度包括污染区的渗透率、污染区的伤害范围等。具体的，通过分析油田开发储层的污染程度可以确定储层污染区的伤害范围即半径，并对污染区和非污染区进行刻画，建立考虑储层污染区和非污染区的双区复合油藏直井和水平井的渗流数学模型。
其中，污染区和非污染区具有不同的储层性质。建立渗流数学模型之前有以下假设：(1)一个圆形的非污染区储层是无限大的，污染区半径为r1；(2)储层水平分布，厚度h均匀，原始压力Pi均匀；(3)对于污染区，水平渗透率为kh1、垂直渗透率为kv1、压力系数为ct1和孔隙度为φ1；对于非污染区，它们分别为kh2、kv2、ct2、φ2；(4)忽略重力和毛管力的影响，考虑井筒储存效应和地层伤害。基于以上假设分别建立污染区和非污染区的流动方程。
例如，对于不同的污染区，考虑如图2、图3和图4所示的三种情况。其中，图2为直井示意图，图3和图4为水平井示意图，对于图3，污染区较大，水平井完全在污染区内，对于图4，污染区较小或者水平井较长，水平井穿过污染区。
以水平井为例，分别考虑两种情况：(1)水平井完全在污染区内；(2)水平井穿过污染区。
在笛卡尔坐标系下，分别建立污染区和非污染区的流动方程，其中，污染区的流动方程可以表示为：
Figure BDA0003312189620000061
非非污染区的流动方程为：
Figure BDA0003312189620000062
初始条件为：
p(x,y,z,t)|t＝0＝pi (4)
外边界条件为：
Figure BDA0003312189620000063
Figure BDA0003312189620000064
交界条件为：
p1＝p2,x＝y＝r1 (7)
Figure BDA0003312189620000071
Figure BDA0003312189620000072
其中，x、y、z分别为笛卡尔坐标系的三个坐标方向；t为时间，单位s；kh1、kv1分别为污染区的水平渗透率和垂直渗透率，单位μm2；kh2、kv2分别为污染区的水平渗透率和垂直渗透率，单位μm2；μ为地层流体粘度，单位mPa·s；p1、p2分别为污染区的压力和非污染区的压力，单位MPa；φ1、φ2分别为污染区和非污染区的孔隙度，为小数；ct1、ct2分别为污染区和非污染区的地层流体综合压缩系数，单位为1/MPa；Pi为储层原始压力，单位为MPa；h为储层厚度，单位m；r1为污染区半径，单位m。
步骤S12，对渗流数学模型进行无因次化分析，确定污染区和非污染区每个位置点的无因次压力。
具体的，对污染区和非污染区进行网格化处理，针对每个网格即位置点对渗流数学模型进行无因次化分析，得到每个位置点的无因次压力。
步骤S13，根据每个位置点的无因次压力对渗流数学模型进行推导，确定有限元积分方程组。
其中，针对每个位置点的无因次压力对渗流数学模型进行推导，得到每个位置点的有限元积分方程，所有位置点的有限元积分方程组成有限元积分方程组。
步骤S14，根据有限元积分方程组预测油田开发过程中储层受到不同污染程度时的弹性产量。
具体的，可以通过有限元数值求解该有限元积分方程组，进而根据求解结果预测油田开发过程中储层受到不同污染程度时的弹性产量。
由此可见，本实施例通过对污染区和非污染区的渗流数学模型进行无因次化分析，确定有限元积分方程组，根据有限元积分方程组预测油田开发过程中储层受到不同污染程度时的弹性产量，能够准确计算在储层污染区边界不规则情况下的油井产能，更贴近地层的实际情况，客观反映污染区和未污染区之间的关联。
在一个可选实施例中，步骤S13具体包括：
步骤S131，确定每个位置点的插值函数。
其中，插值函数可以为假设的形函数或基本函数，如下所示：
N＝(N1,N2,...,Nn) (10)
其中，Nn为第n个位置点对应的差值函数，N为n个位置点的差值函数集合。
其中，为了提高插值函数的计算精度，插值多项式可以为完全多项式，但是完整多项式的项数和单元(计算区域)的自由度数有时并不相同，可以采用减少多项式项数的方法，同时保持多项式的对称性。
步骤S132，根据每个位置点的无因次压力和插值函数确定每个位置点的位移函数。
其中，位移函数为下式：
Figure BDA0003312189620000081
其中，n为位置点数，pDj为位置点j处的无因次压力，Nj为位置点j处的差值函数。
步骤S133，根据每个位置点的位移函数和渗流数学模型得到每个位置点的体积积分。
例如，根据每个位置点的位移函数和渗流数学模型对单元Ωe的体积积分，得到污染区的体积积分：：
Figure BDA0003312189620000082
非污染区的体积积分：
Figure BDA0003312189620000091
其中，piD为污染区的位移函数，LiD为污染区的无因次水平井半长，hiD为污染区的无因次储层厚度，poD为污染区的位移函数，LoD为非污染区的无因次水平井半长，hoD为非污染区的无因次储层厚度，t为时间，单元Ωe为体积积分区域。
步骤S134，根据每个位置点的体积积分确定有限元积分方程组。
具体的，通过对每个位置点的体积积分进行进一步处理，可以得到有限元积分方程组。
在一个可选实施例中，步骤S134具体包括：
步骤S1341，采用格林函数对每个位置点的体积积分进行处理，并去除污染区的零界点和边界点，得到每个位置点的有限元积分方程。
具体的，以污染区为例，采用格林函数对每个位置点的体积积分进行处理，可以得到：
Figure BDA0003312189620000092
其中，
Figure BDA0003312189620000093
其中，Γe为面积积分区域，A为面积。
然后，去除污染区半径为0和r1的区域，得到污染区内部区域的有限元积分方程为：
Figure BDA0003312189620000101
步骤S1342，根据每个位置点的有限元积分方程确定有限元积分方程组。
具体的，由每个位置点的有限元积分方程组成有限元积分方程组。
在一个可选实施例中，步骤S1342具体包括：对每个位置点的有限元积分方程进行形式变换，得到简化的有限元积分方程组。
具体的，根据每个位置点的有限元积分方程，确定如下矩阵：
Figure BDA0003312189620000102
并定义以下参数：
Figure BDA0003312189620000103
则有限元方程变为：
Figure BDA0003312189620000104
为了方便产量计算软件计算，则进一步简化为：
KePm＝Fe (19)
其中，
Figure BDA0003312189620000111
Figure BDA0003312189620000112
在一个可选实施例中，步骤S15具体包括：
步骤S151，对有限元积分方程组进行求解，得到各个时刻储层内各个位置的压力值。
具体的，通过组装整体矩阵和列阵，建立相应的整体有限元方程组，选择二维三角形网格进行计算域的离散，通过有限元法数值求解方程组，选取任意网格单元进行单元特性分析，便能得到各个时刻复合油藏内各个位置的压力值，即p1(x，y，z，t)和p2(x，y，z，t)。
步骤S152，根据各个时刻各个位置的压力值预测油田开发过程中储层受到不同污染程度时的弹性产量。
具体的，可结合弹性产量计算公式，利用有限元数值方法对油田开发过程中储层受到不同污染程度时的弹性产量进行数值求解和分析。
在一个可选实施例中，步骤S152包括：
步骤S1521，根据各个时刻的综合压缩系数、复合油藏体积、油藏原始地层压力和油藏饱和压力确定弹性储量的计算方式。
具体的，根据各个时刻的综合压缩系数、复合油藏体积、油藏原始地层压力和油藏饱和压力可得出靠弹性能驱动的弹性储量的计算方式：
ΔV＝CtVf(Pi-Pb) (22)
其中，ΔV为弹性储量，Ct为t时刻的综合压缩系数，Vf为复合油藏体积，Pi为油藏原始地层压力，Pb为油藏饱和压力。
步骤S1522，根据弹性储量的计算方式和各个时刻各个位置的压力值确定储层受到不同污染程度时的弹性产量的计算方式。
具体的，可根据弹性储量的计算方式确定储层受到不同污染程度时的弹性产量的计算方式。例如，根据各个时刻的综合压缩系数、复合油藏体积、各个时刻各个位置的压力和井底压力确定储层受到不同污染程度时的弹性产量的计算方式。
具体的，通过对整个油藏进行体积积分的方法，求得考虑储层受到不同污染程度时的弹性产量：
Figure BDA0003312189620000121
其中，Np为弹性产量，Pt为t时刻某位置的压力，Pw为井底压力。
实施例二
图5示出了根据本发明实施例二提供的一种油井产能预测方法的流程示意图。其中，本实施例为考虑储层污染区的双区复合油藏直井对应的应用场景实施例。如图5所示，该方法包括：
步骤S21，确定油田开发储层的范围和相关性质，并对直井污染区进行刻画，定义污染区半径、污染区渗透率等参数。
如图2为本实施例中考虑储层污染区的双区复合油藏直井模型的二维示意图，其中直井贯穿整个油层，基本参数如表一所示：
表一
Figure BDA0003312189620000122
Figure BDA0003312189620000131
步骤S22，建立考虑储层污染区和非污染区的双区复合油藏直井渗流数学模型。
其中，直井渗流数学模型参考公式(1)和公式(2)。图6为直井在污染区情况下的径向流流动状态的主视图，图7为直井在污染区情况下的径向流流动状态的俯视图，图8为直井在非污染区情况下的径向流流动状态的俯视图。
步骤S23，对直井渗流数学模型进行无因次化分析，确定污染区和非污染区每个位置点的无因次压力。
步骤S24，根据每个位置点的无因次压力对渗流数学模型进行推导，建立考虑储层污染区和非污染区的双区复合油藏直井的有限元积分方程组。
步骤S25，对有限元积分方程组进行数值求解。
其中，数值求解方法参考公式(10)-(21)。
步骤S26，结合弹性产量计算公式，利用有限元数值方法对油田开发过程中储层受到不同污染程度时的弹性产量进行数值求解和分析。
其中，求解方法参考公式(22)-(23)。
根据上述弹性产量计算结果，分析储层污染区的双区复合油藏中直井污染区大小和渗透率对油井产量的影响。其中，图9为污染区半径的影响示意图，如图9所示，可以看出产量递减动态随污染区大小变化的差异，其中，污染区渗透率定为0.002μm^2。从产量递减曲线可以看出，污染区半径越小，及污染区越小，油藏依靠弹性储量开采的初开井日产油量越高，但是其产量递减的最快。弹性产量处于缓慢递减阶段时，污染区半径越大，产量越高。后期污染区大小对产量几乎没有影响。
图10为污染区渗透率的影响示意图，如图10所示，可以看出产量递减动态随污染区渗透率变化的差异，其中，污染区半径定为300m。从产量递减曲线可以看出，污染区渗透率越小，及污染越严重，油藏依靠弹性能量开采的初开井日产油量越高，但是其产量递减幅度最小。后期产量趋于不变时，污染区渗透率越大，产量越高。
实施例三
图11示出了根据本发明实施三提供的一种油井产能预测方法的流程示意图。本发明实施例考虑到水平井完全在污染区的情景，如图11所示，该方法包括：
步骤S31，确定油田开发储层的范围和相关性质，并对完全在污染区的水平井的污染区进行刻画，定义污染区半径、污染区渗透率等参数。
如图3为本实施例中考虑储层污染区的双区复合油藏直井模型的二维示意图，其中，水平井完全在污染区中，基本参数如表二所示：
表二
Figure BDA0003312189620000141
Figure BDA0003312189620000151
步骤S32，建立考虑储层污染区和非污染区的双区复合油藏完全在污染区的水平井渗流数学模型。
其中，完全在污染区的水平井渗流数学模型参考公式(1)和公式(2)。图12为水平井直井完全在污染区情况下的径向流流动状态的主视图，图13为水平井直井完全在污染区情况下的径向流流动状态的俯视图，图14为水平井直井完全在污染区情况下的径向流流动状态的俯视图。
步骤S33，对渗流数学模型进行无因次化分析，确定污染区和非污染区每个位置点的无因次压力。
步骤S34，根据每个位置点的无因次压力对渗流数学模型进行推导，建立考虑储层污染区和非污染区的双区复合油藏完全在污染区水平井的有限元积分方程组。
步骤S35，对有限元积分方程组进行数值求解。
其中，数值求解方法参考公式(10)-(21)。
步骤S36，结合弹性产量计算公式，利用有限元数值方法对油田开发过程中储层受到不同污染程度时的弹性产量进行数值求解和分析。
其中，求解方法参考公式(22)-(23)。
根据上述弹性产量计算结果，分析储层污染区的双区复合油藏中完全在污染区的水平井的污染区大小和渗透率对油井产量的影响。其中，图15为污染区半径的影响示意图，如图15所示，可以看出产量递减动态随污染区大小变化的差异，其中，污染区渗透率定为0.002μm2。从产量递减曲线可以看出，产量处于递减阶段时，污染区半径越小，产量越高。污染区半径达到一定大小对产量的影响不明显。后期污染区大小对产量几乎没有影响。
图16为污染区渗透率的影响示意图，如图16所示，可以看出产量递减动态随污染区渗透率变化的差异，其中，污染区半径定为300m。从产量递减曲线可以看出，产量处于递减阶段时，污染区渗透率越大，产量越高。后期产量趋于不变时，污染区渗透率对产量几乎没有影响。
实施例四
图17示出了根据本发明实施四提供的一种油井产能预测方法的流程示意图。本发明实施例考虑到水平井穿过污染区的情景，如图17所示，该方法包括：
步骤S41，确定油田开发储层的范围和相关性质，并对穿过污染区的水平井的污染区进行刻画，定义污染区半径、污染区渗透率等参数。
如图4为本实施例中考虑储层污染区的双区复合油藏直井模型的二维示意图，其中，水平井穿过污染区，基本参数如表三所示：
表三
Figure BDA0003312189620000161
步骤S42，建立考虑储层污染区和非污染区的双区复合油藏穿过污染区的水平井渗流数学模型。
其中，穿过污染区的水平井渗流数学模型参考公式(1)和公式(2)。图18为水平井直穿过污染区情况下的径向流流动状态的主视图，图19为水平井直井穿过污染区情况下的径向流流动状态的俯视图。
步骤S43，对渗流数学模型进行无因次化分析，确定污染区和非污染区每个位置点的无因次压力。
步骤S44，根据每个位置点的无因次压力对渗流数学模型进行推导，建立考虑储层污染区和非污染区的双区复合油藏穿过污染区水平井的有限元积分方程组。
步骤S45，对有限元积分方程组进行数值求解。
其中，数值求解方法参考公式(10)-(21)。
步骤S46，结合弹性产量计算公式，利用有限元数值方法对油田开发过程中储层受到不同污染程度时的弹性产量进行数值求解和分析。
其中，求解方法参考公式(22)-(23)。
根据上述弹性产量计算结果，分析储层污染区的双区复合油藏中穿过污染区的水平井的污染区大小和渗透率对油井产量的影响。其中，图20为污染区半径的影响示意图，如图20所示，可以看出产量递减动态随污染区大小变化的差异，其中，污染区渗透率定为0.002μm^2。从产量递减曲线可以看出，污染区半径越大，及污染区越大，产量越低。后期污染区大小对产量几乎没有影响。
图21为污染区渗透率的影响示意图，如图21所示，可以看出产量递减动态随污染区渗透率变化的差异，污染区半径定为150m。从产量递减曲线可以看出，污染区渗透率越小，及污染越严重，产量越低。后期污染区渗透率对产量几乎没有影响。
实施例五
图22示出了根据本发明实施五提供的一种油井产能预测装置的结构示意图。如图22所示，该装置包括：渗流数学模型确定模块51、无因次化分析模块52、方程推导模块53和产量预测模块54；其中，
渗流数学模型确定模块51用于根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型；
无因次化分析模块52用于对所述渗流数学模型进行无因次化分析，确定所述污染区和非污染区每个位置点的无因次压力；
方程推导模块53用于根据所述每个位置点的无因次压力对所述渗流数学模型进行推导，确定有限元积分方程组；
产量预测模块54用于根据所述有限元积分方程组预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
进一步的，所述方程推导模块53包括：
插值函数确定单元531用于确定所述每个位置点的插值函数；
位移函数确定单元532用于根据所述每个位置点的无因次压力和插值函数确定所述每个位置点的位移函数；
积分单元533用于根据所述每个位置点的位移函数和所述渗流数学模型得到所述每个位置点的体积积分；
方程确定单元534用于根据所述每个位置点的体积积分确定有限元积分方程组。
进一步的，所述方程确定单元534具体用于采用格林函数对所述每个位置点的体积积分进行处理，并去除所述污染区的零界点和边界点，得到所述每个位置点的有限元积分方程；根据所述每个位置点的有限元积分方程确定有限元积分方程组。
进一步的，所述方程确定单元534具体用于对所述每个位置点的有限元积分方程进行形式变换，得到简化的有限元积分方程组。
所述产量预测模块54包括：
求解单元541用于对所述有限元积分方程组进行求解，得到各个时刻所述储层内各个位置的压力值；
产量预测单元542用于根据所述各个时刻各个位置的压力值预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
进一步的，所述产量预测单元542具体用于根据各个时刻的综合压缩系数、复合油藏体积、油藏原始地层压力和油藏饱和压力确定弹性储量的计算方式；根据所述弹性储量的计算方式和所述各个时刻各个位置的压力值确定所述储层受到不同污染程度时的弹性产量的计算方式。
进一步的，所述产量预测单元542具体用于根据各个时刻的综合压缩系数、复合油藏体积、各个时刻各个位置的压力和井底压力确定所述储层受到不同污染程度时的弹性产量的计算方式。
本实施例所述的油井产能预测装置用于执行上述实施例一所述的油井产能预测方法，其工作原理与技术效果类似，这里不再赘述。
实施例六
本发明实施例六提供了一种非易失性计算机存储介质，所述计算机存储介质存储有至少一可执行指令，该计算机可执行指令可执行上述任意方法实施例中的油井产能预测方法。
实施例七
图23示出了根据本发明实施例八提供的一种电子设备的结构示意图。本发明具体实施例并不对电子设备的具体实现做限定。
如图23所示，该电子设备可以包括：处理器(processor)602、通信接口(Communications Interface)604、存储器(memory)606、以及通信总线608。
其中：处理器602、通信接口604、以及存储器606通过通信总线608完成相互间的通信。通信接口604，用于与其它设备比如客户端或其它服务器等的网元通信。处理器602，用于执行程序610，具体可以执行上述方法实施例中的相关步骤。
具体地，程序610可以包括程序代码，该程序代码包括计算机操作指令。
处理器602可能是中央处理器CPU，或者是特定集成电路ASIC(ApplicationSpecific Integrated Circuit)，或者是被配置成实施本发明实施例的一个或多个集成电路。电子设备包括的一个或多个处理器，可以是同一类型的处理器，如一个或多个CPU；也可以是不同类型的处理器，如一个或多个CPU以及一个或多个ASIC。
存储器606，用于存放程序610。存储器606可能包含高速RAM存储器，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。
程序610具体可以用于使得处理器602执行上述任意方法实施例中的油井产能预测方法。
在此提供的算法或显示不与任何特定计算机、虚拟系统或者其它设备固有相关。各种通用系统也可以与基于在此的示教一起使用。根据上面的描述，构造这类系统所要求的结构是显而易见的。此外，本发明实施例也不针对任何特定编程语言。应当明白，可以利用各种编程语言实现在此描述的本发明的内容，并且上面对特定语言所做的描述是为了披露本发明的最佳实施方式。
在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
类似地，应当理解，为了精简本发明并帮助理解各个发明方面中的一个或多个，在上面对本发明的示例性实施例的描述中，本发明实施例的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该公开的方法解释成反映如下意图：即所要求保护的本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如下面的权利要求书所反映的那样，发明方面在于少于前面公开的单个实施例的所有特征。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。
本领域那些技术人员可以理解，可以对实施例中的设备中的模块进行自适应性地改变并且把它们设置在与该实施例不同的一个或多个设备中。可以把实施例中的模块或单元或组件组合成一个模块或单元或组件，以及此外可以把它们分成多个子模块或子单元或子组件。除了这样的特征和/或过程或者单元中的至少一些是相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的的替代特征来代替。
此外，本领域的技术人员能够理解，尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在下面的权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
本发明的各个部件实施例可以以硬件实现，或者以在一个或者多个处理器上运行的软件模块实现，或者以它们的组合实现。本领域的技术人员应当理解，可以在实践中使用微处理器或者数字信号处理器(DSP)来实现根据本发明实施例的一些或者全部部件的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的设备或者装置程序(例如，计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上，或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到，或者在载体信号上提供，或者以任何其他形式提供。
应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。单词“包含”不排除存在未列在权利要求中的元件或步骤。位于元件之前的单词“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。上述实施例中的步骤，除有特殊说明外，不应理解为对执行顺序的限定。

Claims (10)

1.一种油井产能预测方法，其特征在于，包括：
根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型；
对所述渗流数学模型进行无因次化分析，确定所述污染区和非污染区每个位置点的无因次压力；
根据所述每个位置点的无因次压力对所述渗流数学模型进行推导，确定有限元积分方程组；
根据所述有限元积分方程组预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
2.根据权利要求1所述的方法，其特征在于，所述根据所述每个位置点的无因次压力对所述渗流数学模型进行推导，确定有限元积分方程组，包括：
确定所述每个位置点的插值函数；
根据所述每个位置点的无因次压力和插值函数确定所述每个位置点的位移函数；
根据所述每个位置点的位移函数和所述渗流数学模型得到所述每个位置点的体积积分；
根据所述每个位置点的体积积分确定有限元积分方程组。
3.根据权利要求2所述的方法，其特征在于，所述根据所述每个位置点的体积积分确定有限元积分方程组，包括：
采用格林函数对所述每个位置点的体积积分进行处理，并去除所述污染区的零界点和边界点，得到所述每个位置点的有限元积分方程；
根据所述每个位置点的有限元积分方程确定有限元积分方程组。
4.根据权利要求3所述的方法，其特征在于，所述根据所述每个位置点的有限元积分方程确定有限元积分方程组，包括：
对所述每个位置点的有限元积分方程进行形式变换，得到简化的有限元积分方程组。
5.根据权利要求1所述的方法，其特征在于，所述根据所述有限元积分方程组预测油田开发过程中所述储层受到不同污染程度时的弹性产量，包括：
对所述有限元积分方程组进行求解，得到各个时刻所述储层内各个位置的压力值；
根据所述各个时刻各个位置的压力值预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
6.根据权利要求5所述的方法，其特征在于，所述根据所述各个时刻各个位置的压力值预测油田开发过程中所述储层受到不同污染程度时的弹性产量，包括：
根据各个时刻的综合压缩系数、复合油藏体积、油藏原始地层压力和油藏饱和压力确定弹性储量的计算方式；
根据所述弹性储量的计算方式和所述各个时刻各个位置的压力值确定所述储层受到不同污染程度时的弹性产量的计算方式。
7.根据权利要求6所述的方法，其特征在于，所述根据所述弹性储量的计算方式和所述各个时刻各个位置的压力值确定所述储层受到不同污染程度时的弹性产量的计算方式，包括：
根据各个时刻的综合压缩系数、复合油藏体积、各个时刻各个位置的压力和井底压力确定所述储层受到不同污染程度时的弹性产量的计算方式。
8.一种油井产能预测装置，其特征在于，包括：
渗流数学模型确定模块，用于根据油田开发储层的污染程度确定包含污染区和非污染区的渗流数学模型；
无因次化分析模块，用于对所述渗流数学模型进行无因次化分析，确定所述污染区和非污染区每个位置点的无因次压力；
方程推导模块，用于根据所述每个位置点的无因次压力对所述渗流数学模型进行推导，确定有限元积分方程组；
产量预测模块，用于根据所述有限元积分方程组预测油田开发过程中所述储层受到不同污染程度时的弹性产量。
9.一种电子设备，包括：处理器、存储器、通信接口和通信总线，所述处理器、所述存储器和所述通信接口通过所述通信总线完成相互间的通信；
所述存储器用于存放至少一可执行指令，所述可执行指令使所述处理器执行如权利要求1-7中任一项所述的油井产能预测方法对应的操作。
10.一种计算机存储介质，所述存储介质中存储有至少一可执行指令，所述可执行指令使处理器执行如权利要求1-7中任一项所述的油井产能预测方法对应的操作。
CN202111219918.1A 2021-10-20 2021-10-20 油井产能预测方法、装置、电子设备及计算机存储介质 Pending CN113887110A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111219918.1A CN113887110A (zh) 2021-10-20 2021-10-20 油井产能预测方法、装置、电子设备及计算机存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111219918.1A CN113887110A (zh) 2021-10-20 2021-10-20 油井产能预测方法、装置、电子设备及计算机存储介质

Publications (1)

Publication Number Publication Date
CN113887110A true CN113887110A (zh) 2022-01-04

Family

ID=79003810

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111219918.1A Pending CN113887110A (zh) 2021-10-20 2021-10-20 油井产能预测方法、装置、电子设备及计算机存储介质

Country Status (1)

Country Link
CN (1) CN113887110A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115492561A (zh) * 2022-09-16 2022-12-20 重庆地质矿产研究院 一种页岩气井注入二氧化碳增产增效与地质封存一体化方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110656923A (zh) * 2018-06-28 2020-01-07 中国石油化工股份有限公司 一种碳酸盐岩缝洞型油藏生产井产量预测方法
CN111980654A (zh) * 2020-10-12 2020-11-24 西南石油大学 一种非均匀页岩油藏分段压裂水平井产能计算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110656923A (zh) * 2018-06-28 2020-01-07 中国石油化工股份有限公司 一种碳酸盐岩缝洞型油藏生产井产量预测方法
CN111980654A (zh) * 2020-10-12 2020-11-24 西南石油大学 一种非均匀页岩油藏分段压裂水平井产能计算方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
向开理, 郭建春, 李允: "侧钻井压力分析的边界元法及储层伤害的评价", 计算物理, no. 05, 20 September 2000 (2000-09-20), pages 579 - 587 *
孙建孟;运华云;冯春珍;: "测井产能预测方法与实例", 测井技术, no. 06, 20 December 2012 (2012-12-20), pages 628 - 634 *
岳建伟: "低渗透气藏压裂井产能研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》, no. 2006, 15 April 2006 (2006-04-15), pages 019 - 45 *
李成勇: "格林函数与源函数方法在试井分析中的应用研究", 《中国优秀硕士学位论文全文数据库 工程科技I辑》, no. 2006, 15 April 2006 (2006-04-15), pages 019 - 32 *
杨悦;李相方;卢巍;李异才;吴克柳;羊新州;龚崛;: "浅层低渗油藏压裂水平裂缝直井产能方程推导及应用", 科学技术与工程, vol. 13, no. 24, 28 August 2013 (2013-08-28), pages 7015 - 7020 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115492561A (zh) * 2022-09-16 2022-12-20 重庆地质矿产研究院 一种页岩气井注入二氧化碳增产增效与地质封存一体化方法
CN115492561B (zh) * 2022-09-16 2023-03-21 重庆地质矿产研究院 一种页岩气井注入二氧化碳增产增效与地质封存一体化方法

Similar Documents

Publication Publication Date Title
EP2253797B1 (fr) Méthode d'exploitation de milieu poreux au moyen d'une modélisation d'écoulements de fluide
CN108798634A (zh) 一种缝洞型储层单洞-单层均值试井解释方法
CN101446196A (zh) 三重介质油藏分支水平井的试井分析方法及装置
CN114169263B (zh) 基于压裂液返排数据的页岩油气藏裂缝参数反演方法
FR2961844A1 (fr) Ecoulement multiphase dans un puits de forage et fracture hydraulique connectee
CN109684685B (zh) 一种多分支井水合物降压开采条件下的产能及储层稳定性分析方法
CN114818533B (zh) 基于排采数据的页岩油气藏裂缝参数确定方法及装置
CN116306385B (zh) 一种油藏压裂渗吸增能数值模拟方法、系统、设备及介质
CN113887110A (zh) 油井产能预测方法、装置、电子设备及计算机存储介质
CN113255092A (zh) 体积压裂水平井特征点拟合反演方法、电子设备及介质
Nakashima et al. Near-well upscaling for three-phase flows
CN111577264A (zh) 裂缝孔隙型油藏水平井产能预测方法及装置
CN112377182A (zh) 大型孔洞裂缝型碳酸盐岩储集体参数确定方法及装置
Ding Modeling formation damage for flow simulations at reservoir scale
Milliotte et al. Well-data-based discrete fracture and matrix modelling and flow-based upscaling of multilayer carbonate reservoir horizons
WO2017039622A1 (en) Integrated workflow for feasibility study of cuttings reinjection based on 3-d geomechanics analysis
CN112377184A (zh) 串珠状缝洞型碳酸盐岩储层物性参数分析方法及装置
FR3046810A1 (fr) Procede de production d'hydrocarbures comportant un index de productivite des puits sous effet thermique
CN105649610B (zh) 一种获取油藏压力的方法及装置
CN112419493A (zh) 页岩储层三维属性模型建立方法及装置
CN114086933B (zh) 一种试井解释方法、装置、存储介质以及计算机设备
Khramchenkov et al. Mathematical modeling and experimental study of erosion-deposition process in deformable porous media
CN112682033B (zh) 孔洞裂缝型碳酸盐岩储集体参数确定方法及装置
CN112377185B (zh) 单洞与多重介质复合型储层参数敏感性分析方法及装置
CN113919242B (zh) 微粒运移损害油气层建模方法、损害程度时空演化4d定量与智能诊断方法及其系统

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