CN108664679A - 一种油气井生产数据分析方法 - Google Patents

Abstract

一种油气井生产数据分析方法，其包括：根据原始地层压力、井底流压和日产量构建待分析井的压力规整化产量函数；根据压力规整化产量函数，基于累产量确定双对数图版绘制函数，根据双对数图版绘制函数绘制压力规整化产量双对数图版；根据获取到的压力规整化产量双对数图版中的实测数据对压力规整化产量双对数图版进行拟合，并根据双对数图版拟合结果进行产量与压力历史拟合，根据产量与压力历史拟合结果确定出待分析井的相关参数。本方法能够有效避开物质物质平衡时间随真实时间不单调变化的问题，消除了因生产数据间断、大的波动给生产数据解释带来的负面影响，同时能实现双对数图版数据和生产历史数据的同步拟合。

Description

一种油气井生产数据分析方法

技术领域

本发明涉及油气井勘探开发技术领域，具体地说，涉及一种油气井生产数据分析方法。

背景技术

试井是以渗流力学为理论基础，以测试仪表为手段，通过测量油气、水、井的压力以及产量等生产动态数据，研究和确定测试井和油、气、水层的生产能力、物性参数以及辨别井间或层间连通关系的技术，其对油气田新区勘探开发、老区动态调整都具有极为重要的作用。

为了解某特定井及对应储层的参数，根据所利用的资料来源不同，主要有压力恢复试井(简称压恢试井)、生产数据试井分析两种常用方法。然而，现有的生产数据分析方法对部分油气井并不适用，这使得生产数据分析的解释结果具有很强的不确定性和多解性。

发明内容

为解决上述问题，本发明提供了一种油气井生产数据分析方法，所述方法包括：

步骤一、根据原始地层压力、井底流压和日产量构建待分析井的压力规整化产量函数或拟压力规整化产量函数；

步骤二、根据所述压力规整化产量函数或拟压力规整化产量函数，基于累产量确定双对数图版绘制函数，根据所述双对数图版绘制函数绘制压力规整化产量双对数图版或拟压力规整化产量双对数图版；

步骤三、根据获取到的压力规整化产量双对数图版或拟压力规整化产量双对数图版中的实测数据对所述压力规整化产量双对数图版或拟压力规整化产量双对数图版进行拟合，并根据双对数图版拟合结果进行产量与压力历史拟合，根据产量与压力历史拟合结果确定出待分析井的相关参数。

根据本发明的一个实施例，在所述步骤一中，

如果所述待分析井为油井，则根据如下表达式构建所述压力规整化产量函数：

如果所述待分析井为气井，则根据如下表达式计算所述拟压力规整化产量函数：

其中，PNR表示压力规整化产量或拟压力规整化产量，p_i表示原始地层压力，p_wf表示井底压力，t表示时间，q表示产量，ψ_i表示气井原始地层拟压力，ψ_wf表示气井井底拟压力。

根据本发明的一个实施例，在所述步骤二中，

根据所述压力规整化产量函数确定压力规整化产量积分函数，并根据所述压力规整化产量积分函数确定压力规整化产量积分导数函数，根据所述压力规整化产量函数、压力规整化产量积分函数和压力规整化产量积分导数函数构建所述压力规整化产量双对数图版。

根据本发明的一个实施例，根据如下表达式确定所述压力规整化产量积分函数：

其中，PNR_i表示压力规整化产量积分，V_e表示在物质平衡时间的累产量，PNR表示压力规整化产量，V表示累产量。

根据本发明的一个实施例，根据如下表达式确定所述压力规整化产量积分导数函数：

其中，PNR_i表示压力规整化产量积分，PNR_id表示压力规整化产量积分导数，V表示累产量。

根据本发明的一个实施例，在所述步骤三中，通过调整累产拟合值和压力规整化产量拟合值来调整压力规整化产量双对数图版的拟合效果。

根据本发明的一个实施例，根据如下表达式确定所述累产拟合值：

根据如下表达式确定所述压力规整化产量拟合值：

其中，V_M和PNR_M分别表示累产拟合值和压力拟合值，V_TB和V_SJ分别表示双对数曲线中拟合点的图版累产量数据和实测累产量数据，PNR_TB和PNR_SJ分别表示双对数曲线中拟合点的图版压力规整化产量数据和实测压力规整化产量数据。

根据本发明的一个实施例，在所述步骤三中，还基于双对数图版拟合结果进行产量和压力数据的历史拟合，通过历史拟合来对所述双对数图版拟合结果进行优化，实现图版数据和生产历史数据的同步拟合，从而得到优化后的待分析井的相关参数。

根据本发明的一个实施例，在所述步骤三中，

步骤a、基于本次双对数图版拟合结果，根据各个时刻的实测累产量以及累产量拟合值计算双对数图版上对应的累产量图版值；

步骤b、根据所述累产量图版值在所述双对数图版上查找该累产量图版值所对应的压力规整化产量图版值，根据所述压力规整化产量图版值和压力规整化产量拟合值计算对应的压力规整化产量实测值；

步骤c、根据所述压力规整化产量实测值计算井底压力拟合数据；

步骤d、根据实测井底压力数据和步骤b中得到的压力规整化产量实测值计算日产量拟合数据，根据所述日产量拟合数据计算各个时刻的累产量拟合数据；

步骤e、根据各个时刻的井底压力拟合数据、日产量拟合数据以及累产量拟合数据与实测压力井底压力数据、实测日产量数据以及累产量实测数据进行拟合，并根据拟合结果对所述双对数图版拟合结果进行优化。

本发明所提供的油气井生产数据分析方法基于(拟)压力规整化产量与累产量的双对数图版，利用油气井生产数据通过图版拟合、产量和压力数据的历史拟合进行储层参数(如渗透率、表皮系数)和井筒及改造参数(如压裂裂缝半长、裂缝导流能力、SRV体积)解释，这样能够有效避开物质物质平衡时间随真实时间不单调变化的问题，消除了因生产数据间断、大的波动给生产数据解释带来的负面影响，从而提高了生产数据(尤其是日产量)间断、波动较大的油气井的生产数据解释精度，提高了产量和压力历史拟合效果，大大提高了解释精度，降低了多解性。

本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要的附图做简单的介绍：

图1是根据本发明一个实施例的油气井生产数据分析方法的实现流程示意图；

图2是根据本发明一个实施例的双对数图版示意图；

图3是根据本发明一个实施例的生产数据双对数图版拟合示意图；

图4是现有方法的生产数据双对数图版拟合示意图；

图5是根据本发明一个实施例的历史拟合示意图；

图6是现有方法的历史拟合示意图。

具体实施方式

以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例以及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。

同时，在以下说明中，出于解释的目的而阐述了许多具体细节，以提供对本发明实施例的彻底理解。然而，对本领域的技术人员来说显而易见的是，本发明可以不用这里的具体细节或者所描述的特定方式来实施。

另外，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。

为了解某特定井及对应储层的参数，根据所利用的资料来源不同，主要有压力恢复试井(简称压恢试井)、生产数据试井分析两种常用方法。

压恢试井是一种传统试井方法，其应用时在井底或井口放置压力井，关井后记录压力变化。该方法的优点是压力计采点密精度高，解释结果能反映油气井开井生产早期信息。但由于页岩气藏储层渗透率极低，压恢试井只能记录短期的部分数据，无法反映页岩气井诸如裂缝干扰、裂缝边界等远井信息。因此压恢试井在页岩气井中能提供的信息非常有限。

生产数据试井分析是利用开井生产的动态数据进行试井解释，主要方法有PNR方法(压力归一化流量方法，也称为Blasingame方法)、Agarwal-Gardner方法、RNP方法(流量归一化压力方法)。生产数据试井方法优点是不需要关井，不影响正常生产，通过记录的长期生产数据来解释井和储层特征参数，能反映油气井远井信息，因此这种方法非常适用于低渗透储层的油气井。

但上述现有的生产数据分析方法仍存在两个问题。其一，现有的生产数据分析方法主要通过(拟)压力和(拟)压力导数与物质平衡时间双对数图版的拟合来进行试井分析，是在分析(拟)压力和(拟)压力导数与物质平衡时间的关系。当物质平衡时间随真实生产时间单调递增时，双对数图版拟合较好。但由于气井在生产过程中经常伴随频繁的生产制度改变(如开关井、排液、更换油嘴尺寸、井间干扰)或人为原因导致生产状态突然变化，因此生产数据质量通常较差，压力和流量数据经常存在间断、波动甚至缺失，导致物质平衡时间随真实生产时间波动幅度很大以至于无法单调递增。在这种情况下压力和压力导数与物质平衡时间双对数图版的拟合效果则无法保证，则解释结果具有很强的不确定性和多解性。

其二，现有的生产数据分析方法主要通过(拟)压力和(拟)压力导数与物质平衡时间双对数图版的拟合来进行试井分析，双对数图版的拟合与产量与压力数据的历史拟合是独立进行的，无法保证同步保证好的图版拟合和历史拟合效果。

上述两个问题使得现有的生产数据分析方法对部分油气井并不适用，双对数图版拟合和产量与压力历史拟合无法保证同步拟合，解释结果具有很强的不确定性和多解性。

针对现有技术中所存在的上述问题，本发明提供了一种新的油气井生产数据分析方法，该方法采用(拟)压力规整化产量以及(拟)压力规整化产量积分、(拟)压力规整化产量积分导数与累产量Blasingame双对数关系图版进行分析，保证了双对数关系图版中的横坐标(累产量)随真实时间的单调递增，从而解决了因频繁的生产制度改变(如开关井、排液、更换油嘴尺寸、井间干扰)或人为原因导致的生产数据间断、大的波动而现有生产数据试井解释方法难以准确解释的难题，并保证了双对数图版和产量与压力历史的同步拟合效果。

图1示出了本实施例中油气井生产数据分析方法的实现流程示意图。

本实施例所提供的油气井生产数据分析方法在实现过程中需要利用到井底压力与生产时间内的累产量的关系模型，因此首先需要构建得到井底压力与生产时间内的累产量的关系模型。

对于随时间变化的变流量和变压力数据，可通过如下表达式可以计算生产到t时刻的累产量：

其中，V表示累产量，t表示生产时间，q表示日产量。

以无限大均质地层中一口直井油井为例，变产量生产时地层压力分布可表示为：

其中，p表示压力，p_i表示原始地层压力，r表示地层中的点与井筒的距离，B表示体积系数，h表示地层厚度，k表示渗透率，μ表示气体粘度。

χ表示导压系数，其可以采用如下表达式计算得到：

其中，φ表示孔隙度，C_t表示压缩系数。

当时间较大时即存在：

此时表达式(2)可以表示成：

其中，t₁表示一趋近于0的时间点。

将表达式(1)代入表达式(5)中，井底压力可以表示成：

其中，p_wf表示井底压力，表示平均产量，a、b和c均表示方程系数。

从表达式(6)可以看出，对于无限大均质地层中一口直井油井，当流动达到径向流时，井底压力与产量的对数之间为线性关系。

以无限大均质地层中一口垂直裂缝井为例，变产量生产时地层压力分布可表示为：

其中，x_f表示裂缝半长，x和y表示地层中任意一点与井的距离。erf(x)表示误差函数，其可以采用如下表达式表示：

其井底压力可以表示成：

根据误差函数的性质，当时间较小，即时，此时存在：

此时，表达式(9)可以近似为：

将表达式(1)代入表达式(11)，其井底压力可以表示成：

从表达式(12)可以看出，对于无限大均质地层中垂直裂缝井，当流动达到线性流时，井底压力与累产量的平方根之间为线性关系。

对于气井，利用表达式(13)的拟压力代替表达式(6)和(12)中的压力，也可得出如表达式(6)和表达式(12)的方程，以表示拟压力与累产量之间的关系。

其中，ψ表示气井拟压力，z表示气体压缩因子，p₀表示参考压力(优选地为大气压)。

如图1所示，本实施例所提供的油气井生成数据分析方法首先在步骤S101中根据原始地层压力、井底流压和日产量构建待分析井的压力规整化产量函数或拟压力规整化产量函数。

本实施例中，对于日产量和压力随时间变化的油气井，本方法根据获得的原始地层压力、井底流压、日产量来计算压力规整化产量或拟压力规整化产量。若未测量井底流压数据，该方法则可以通过压力折算将井口测量的套压或油压数据折算至井底流压。

具体地，如果待分析井为油井，那么本方法优选地根据如下表达式来构建压力规整化产量函数：

如果待分析井为气井，类似地，本方法则优选地根据如下表达式来构建拟压力规整化产量函数：

其中，p_i表示原始地层压力，p_wf表示井底压力，t表示时间，q表示产量，ψ_i表示气井原始地层拟压力，ψ_wf表示气井井底拟压力。对于油井而言，PNR表示压力规整化产量；对于气井而言，PNR则表示拟压力规整化产量。

由于对于油井和气井而言，二者的生产数据分析过程类似，故本实施例中以油井为例来进行后续的说明。

如图1所示，本实施例中，在得到待分析井的压力规整化产量函数后，该方法会在步骤S102中根据上述压力规整化产量函数确定双对数图版绘制函数，并在步骤S103中根据上述双对数图版绘制函数构建压力规整化产量双对数图版。具体地，本实施例中，该方法所确定出的双对数图版绘制函数优选地包括：压力规整化产量积分函数、压力规整化产量积分导数函数以及压力规整化产量函数本身。

其中，本实施例中，压力规整化产量积分函数优选地可以采用如下表达式进行表示：

压力规整化产量积分导数函数优选地可以采用如下表达式进行表示：

其中，PNR_i表示压力规整化产量积分，PNR_id表示压力规整化产量积分导数，V表示累产量，V_e表示在物质平衡时间的累产量。

本实施例中，对于油井来说，物质平衡时间t_e可以采用如下表达式计算得到：

而对于气井来说，物质平衡时间t_e则可以采用如下表达式计算得到：

在物质平衡时间的累产量V_e则可以通过将上述物质平衡时间t_e代入表达式(1)中计算得到。

对于给定的油气井类型、油气藏类型、内边界类型、外边界类型、流体类型，通过上述部分中所建立的压力与累产量的关系模型以及得到的双对数图版绘制函数，以累产量为横坐标，在双对数坐标上可以绘制3条曲线，即：压力规整化产量曲线、压力规整化产量积分曲线以及压力规整化产量积分导数曲线，从而得到如图2所示的以累产量为横坐标并以压力规整化产量、压力规整化产量积分以及压力规整化产量积分导数为坐标的Blasingame双对数图版示意图。

利用生产数据进行油气井试井分析时，需要进行图版拟合以及产量、压力的历史拟合。如图1所示，本实施例所提供的方法在步骤S104中根据获取到的压力规整化产量双对数图版的实测数据对产量规整化双对数图版进行拟合。

具体地，对于油气井的压力规整化产量双对数图版，拖动双对数图版上的理论模型曲线，使之与实测数据的双对数曲线拟合。在图版拟合过程中，将图版拖动的水平距离称为累产量拟合值，垂直距离称为压力规整化产量拟合值。

其中，本实施例中，累产拟合值定义如下：

压力规整化产量拟合值定义如下：

其中，V_M和PNR_M分别表示累产量拟合值和压力规整化产量拟合值，V_TB和V_SJ分别表示双对数曲线中拟合点的图版累产量数据和实测累产量数据，PNR_TB和PNR_SJ分别表示双对数曲线中拟合点的图版压力规整化产量数据(即压力规整化产量图版值)和实测压力规整化产量数据(即压力规整化产量实测值)。

本实施例中，压力规整化产量双对数图版的拟合过程是通过调整累产量拟合值和压力规整化产量拟合值来调整拟合效果来使得拟合效果达到较好效果的过程。当拟合效果达到预设要求时，双对数图版的拟合过程也即停止。

为了使得最终得到的待分析井的相关参数(例如储层、井筒和/或改造参数)的结果更加准确可靠，如图1所示，本实施例所提供的方法还会在步骤S105中根据双对数图版拟合结果进行历史数据拟合，以实现对双对数图版和历史数据的同步拟合，从而得到优化后的待分析井的相关参数。

具体地，本实施例中，该方法通过本次图版拟合过程，可以得到包括一组累产量拟合值V_M、压力规整化产量拟合值PNR_M以及待分析井的储层、井筒、改造参数的图版拟合结果。对于任一时刻的生产数据，利用表达式(1)可以计算得到对应该时刻的实测累产量V_SJ，利用表达式(20)可以根据上述实测累产量V_SJ以及累产量拟合值V_M计算得到双对数图版上对应的累产量图版值V_TB。

在得到上述累产量图版值V_TB后，可以在双对数图版上查找到与该累产量图版值V_TB相对应的压力规整化产量的图版值PNR_TB(即图版压力规整化产量数据)。利用表达式(21)，根据上述压力规整化产量的图版值PNR_TB和压力规整化产量拟合值PNR_M可以计算得到对应的压力规整化产量的实测值PNR_SJ(即实测常量规整话压力数据)。

由于对应时刻的日产量已知，利用表达式(14)便可以根据压力规整化产量的实测值PNR_SJ计算得到井底压力的拟合数据。同时，由于对应时刻的实测井底压力已知，利用表达式(14)可以根据实测得到的井底生产压差(即原始压力减去井底压力)和压力规整化产量的实测值PNR_SJ计算得到日产量拟合数据。根据该时刻的日产量拟合数据以及上一时刻的累产量，也可以计算得到该时刻的累产量拟合数据。

基于相同原理，该方法可以得到从时间0到生产时间最末点中各个时刻的压力拟合数据、日产量拟合数据以及累产量拟合数据。在得到上述拟合数据后，该方法会将上述压力拟合数据、日产量拟合数据以及累产量拟合数据与实测压力数据、实测日产量数据以及实测累产量数据分别进行拟合。其中，如果能够满足预设拟合精度要求(例如拟合数据与实测数据之差的绝对值小于或等于预设差值阈值)，那么则停止拟合，并将本次双对数图版拟合所得到的相关参数(例如储层参数、井筒参数和/或改造参数)作为最终的解释结果。而如果无法满足预设拟合精度要求，那么则调整双对数图版拟合参数并重新进行图版拟合以及历史拟合过程，直至满足预设拟合精度要求。

需要指出的是，本实施例中，根据实际需要，本方法所采用的双对数图版拟合过程以及历史数据拟合过程也可以与现有技术相同，故在此不再对步骤S104以及步骤S105的相关内容进行赘述。

以某油田一口气井为例说明本发明的应用。该井是水平井，完钻井深4168.0m，水平段长1533m，压裂22段，储层原始压力为36.7MPa，地层温度81℃，孔隙度4％。该井生产了17184小时，日产量和压力数据为每24小时1个点，共711组数据。

应用本发明的方法，对该井的生产数据进行试井分析，如图3所示将生产数据进行处理并与理论图版进行拟合，拟合结果为裂缝平均半长为77m，地层渗透率为0.024md，井筒表皮因子0.1，地层边界为1184m×307m，累产拟合值V_M为9.46，压力拟合值P_M为0.054MPa^-1。

为了对比，应用现有的方法，对该井的生产数据进行试井分析，如图5所示将生产数据进行处理并与理论图版进行拟合，拟合结果为裂缝平均半长为67.6m，地层渗透率为0.028md，井筒表皮因子0.1，地层边界为1172m×275m，时间拟合值T_M为14.19(1/hr)，压力拟合值P_M为0.0625MPa^-1。

根据图版拟合所得到的结果，分别应用本发明的方法和现有方法，进行日产量和压力数据的历史拟合，结果分别如图4和图6所示。其中，由本发明方法得到的历史拟合如果如图4所示，应用现有方法的历史拟合结果如图6所示。

通过图4至图6可以看出，由于该井在4000和8000小时时日产量发生突变，导致其物质平衡时间随真实时间产生大的波动，破坏了物质平衡时间与真实时间之间的单调性，因此用现有方法进行解释时，在拟压力规整化产量与物质平衡时间的Blasingame双对数理论曲线上出现了明显的拐点(对应约4000小时物质平衡时间处)。由于该拐点前后是不同的趋势，所以在图5和图6上都只能在拐点之前或拐点之后得到较好的拟合效果，而无法得到整条曲线的好的拟合效果。因此通过现有方法得到的解释结果可信度低，且结果的多解性较强。

而本发明的方法避开了物质平衡时间随真实时间不单调变化的问题，因此在图3的拟压力规整化产量与累产量的双对数图版上拟合效果非常好，也改进了图4所示的日产量、累产量和压力历史拟合效果。因此通过本发明的方法得到的结果可信度较高，降低了结果的多解性。

通过对比可以看出，本发明提出的基于(拟)压力规整化产量与累产量的双对数图版的生产数据试井分析方法，避开了物质平衡时间随真实时间不单调变化的问题，提高了油气井产量和压力历史拟合效果，大大提高了解释精度，降低了多解性，从而验证了本发明的正确性和实用性。本发明尤其提高了生产数据(尤其是日产量)波动较大的油气井的生产数据解释精度，当生产数据(尤其是日产量)变化较小时，该方法同样适用，因此本发明具有广泛的实用性。

从上述描述中可以看出，本实施例所提供的油气井生产数据分析方法基于(拟)压力规整化产量与累产量的双对数图版，利用油气井生产数据通过图版拟合、产量和压力数据的历史拟合进行储层参数(如渗透率、表皮系数)和井及改造参数(如压裂裂缝半长、裂缝导流能力、SRV体积)解释，这样能够有效避开物质物质平衡时间随真实时间不单调变化的问题，消除了因生产数据间断、大的波动给生产数据解释带来的负面影响，从而提高了生产数据(尤其是日产量)间断、波动较大的油气井的生产数据解释精度，提高了产量和压力历史拟合效果，大大提高了解释精度，降低了多解性。

应该理解的是，本发明所公开的实施例不限于这里所公开的特定结构或处理步骤，而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是，在此使用的术语仅用于描述特定实施例的目的，而并不意味着限制。

说明书中提到的“一个实施例”或“实施例”意指结合实施例描述的特定特征、结构或特性包括在本发明的至少一个实施例中。因此，说明书通篇各个地方出现的短语“一个实施例”或“实施例”并不一定均指同一个实施例。

虽然上述示例用于说明本发明在一个或多个应用中的原理，但对于本领域的技术人员来说，在不背离本发明的原理和思想的情况下，明显可以在形式上、用法及实施的细节上作各种修改而不用付出创造性劳动。因此，本发明由所附的权利要求书来限定。

Claims (9)

1.一种油气井生产数据分析方法，其特征在于，所述方法包括：

步骤一、根据原始地层压力、井底流压和日产量构建待分析井的压力规整化产量函数或拟压力规整化产量函数；

步骤二、根据所述压力规整化产量函数或拟压力规整化产量函数，基于累产量确定双对数图版绘制函数，根据所述双对数图版绘制函数绘制压力规整化产量双对数图版或拟压力规整化产量双对数图版；

步骤三、根据获取到的压力规整化产量双对数图版或拟压力规整化产量双对数图版中的实测数据对所述压力规整化产量双对数图版或拟压力规整化产量双对数图版进行拟合，并根据双对数图版拟合结果进行产量与压力历史拟合，根据产量与压力历史拟合结果确定出待分析井的相关参数。

2.如权利要求1所述的方法，其特征在于，在所述步骤一中，

如果所述待分析井为油井，则根据如下表达式构建所述压力规整化产量函数：

如果所述待分析井为气井，则根据如下表达式计算所述拟压力规整化产量函数：

其中，PNR表示压力规整化产量或拟压力规整化产量，p_i表示原始地层压力，p_wf表示井底压力，t表示时间，q表示产量，ψ_i表示气井原始地层拟压力，ψ_wf表示气井井底拟压力。

3.如权利要求1或2所述的方法，其特征在于，在所述步骤二中，

根据所述压力规整化产量函数确定压力规整化产量积分函数，并根据所述压力规整化产量积分函数确定压力规整化产量积分导数函数，根据所述压力规整化产量函数、压力规整化产量积分函数和压力规整化产量积分导数函数构建所述压力规整化产量双对数图版。

4.如权利要求3所述的方法，其特征在于，根据如下表达式确定所述压力规整化产量积分函数：

其中，PNR_i表示压力规整化产量积分，V_e表示在物质平衡时间的累产量，PNR表示压力规整化产量，V表示累产量。

5.如权利要求3或4所述的方法，其特征在于，根据如下表达式确定所述压力规整化产量积分导数函数：

其中，PNR_i表示压力规整化产量积分，PNR_id表示压力规整化产量积分导数，V表示累产量。

6.如权利要求1～5中任一项所述的方法，其特征在于，在所述步骤三中，通过调整累产拟合值和压力规整化产量拟合值来调整压力规整化产量双对数图版的拟合效果。

7.如权利要求6所述的方法，其特征在于，根据如下表达式确定所述累产拟合值：

根据如下表达式确定所述压力规整化产量拟合值：

其中，V_M和PNR_M分别表示累产拟合值和压力规整化产量拟合值，V_TB和V_SJ分别表示双对数曲线中的图版累产量数据和实测累产量数据，PNR_TB和PNR_SJ分别表示双对数曲线中的图版压力规整化产量数据和实测压力规整化产量数据。

8.如权利要求1～7中任一项所述的方法，其特征在于，在所述步骤三中，还基于双对数图版拟合结果进行产量和压力数据的历史拟合，通过历史拟合来对所述双对数图版拟合结果进行优化，实现图版数据和生产历史数据的同步拟合，从而得到优化后的待分析井的相关参数。

9.如权利要求8所述的方法，其特征在于，在所述步骤三中，

步骤a、基于本次双对数图版拟合结果，根据各个时刻的实测累产量以及累产量拟合值计算双对数图版上对应的累产量图版值；

步骤b、根据所述累产量图版值在所述双对数图版上查找该累产量图版值所对应的压力规整化产量图版值，根据所述压力规整化产量图版值和压力规整化产量拟合值计算对应的压力规整化产量实测值；

步骤c、根据所述压力规整化产量实测值计算井底压力拟合数据；

步骤d、根据实测井底压力数据和步骤b中得到的压力规整化产量实测值计算日产量拟合数据，根据所述日产量拟合数据计算各个时刻的累产量拟合数据；

步骤e、根据各个时刻的井底压力拟合数据、日产量拟合数据以及累产量拟合数据与实测压力井底压力数据、实测日产量数据以及累产量实测数据进行拟合，并根据拟合结果对所述双对数图版拟合结果进行优化。