CN110175412B - 油藏数据拟合方法、油藏数据分析方法及装置 - Google Patents

Abstract

本申请提供了一种油藏数据拟合方法、油藏数据分析方法及装置，其中，该油藏数据拟合方法包括：根据多组实测油藏数据获取压力数据及压力相关数据；根据压力相关数据的变化趋势确定出压力相关数据中的关键数据节点；根据关键数据节点将压力数据形成的第一变化趋势信息及压力相关数据形成的第二变化趋势信息分成多段数据段；从多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段；将每个目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出目标油藏模型。根据上述方法可以提高目标油藏模型拟合的收敛速度。

Description

油藏数据拟合方法、油藏数据分析方法及装置

技术领域

本申请涉及数据处理领域，具体而言，涉及一种油藏数据拟合方法、油藏数据分析方法及装置。

背景技术

随着计算机技术的进步，自动拟合技术被引入到试井分析过程，其主要思想是将试井分析问题转化为数学意义下的非线性最小二乘问题。但是对上述的试井分析过程中的非线性最小二乘问题的解决的手段存在收敛速度较慢，算法需要耗费大量时间等不足。

发明内容

有鉴于此，本申请实施例的目的在于提供一种油藏数据拟合方法、油藏数据分析方法及装置。

第一方面，本申请实施例提供了一种油藏数据拟合方法，包括：

根据多组实测油藏数据获取压力数据及压力相关数据；

根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点；

根据所述关键数据节点将所述压力数据形成的第一变化趋势信息及所述压力相关数据形成的第二变化趋势信息分成多段数据段；

从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段；

将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型。

本申请实施例提供的油藏数据拟合方法，采用两种变化趋势信息中的多段数据段为不同的目标参数进行校准，可以有针对性的拟合出目标参数，从而可以提高拟合的收敛速度。由于使用不同的数据段对不同的目标参数进行拟合校准，可以降低算法对参数的初值的敏感度，从而可以提高油藏模型的拟合速度。

结合第一方面，本申请实施例提供了第一方面的第一种可能的实施方式，其中：所述目标参数包括：渗透率表皮系数及井储参数，所述目标数据段包括用于验证所述渗透率的第一目标数据段、用于验证所述表皮系数的第二目标数据段及用于验证所述井储参数的第三目标数据段；所述将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型的步骤，包括：使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到误差值小于第一设定值；使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到误差值小于第二设定值；使用第三设定拟合方式对所述井储参数进行拟合，并使用所述第三目标数据段对拟合所述井储参数后对应的所述目标油藏模型进行验证，直到误差值小于第三设定值；重复上述的步骤，直到所述目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差小于目标设定值，确定出目标参数，以得到所述目标油藏模型。

进一步，本申请实施例提供的油藏数据拟合方法，还可以使用三段数据段分别对试井分析中的参数：井储参数、表皮系数及渗透率进行拟合校准，从而可以有效拟合出有用于对气田进行分析的目标油藏模型。

结合第一方面第一种可能的实施方式，本申请实施例提供了第一方面的第二种可能的实施方式，其中：所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述第一目标数据段属于所述第二变化趋势信息中的一数据段；所述使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到误差值小于第一设定值的步骤，包括：使用第一设定拟合方式确定出所述渗透率的当前值，形成当前的所述导数函数；使用所述第一目标数据段对所述导数函数进行验证；判断验证得到的所述导数函数与所述第一目标数据段的误差值是否小于第一设定值；若所述误差值不小于所述第一设定值，则重复上述的步骤，直到所述误差值小于所述第一设定值。

进一步，本申请实施例提供的油藏数据拟合方法，上述的目标油藏模型包括压力函数及所述压力函数对应的导数函数，通过对其中的导数函数中的渗透率进行拟合，使最后拟合目标油藏模型时的速度更快，从而使目标油藏模型与实测数据更匹配。

结合第一方面第一种可能的实施方式，本申请实施例提供了第一方面的第三种可能的实施方式，其中：所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述第二目标数据段属于所述第一变化趋势信息中的一数据段，所述使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到误差值小于第二设定值的步骤，包括：使用第二设定拟合方式确定出所述表皮系数的当前值，形成当前的所述压力函数；使用所述第二目标数据段对所述压力函数进行验证；判断验证得到所述压力函数与所述第二目标数据段的误差值是否小于第二设定值；若所述误差值不小于所述第二设定值，则重复上述的步骤，直到所述误差值小于所述第二设定值。

进一步，本申请实施例提供的油藏数据拟合方法，上述的目标油藏模型包括压力函数及所述压力函数对应的导数函数，通过对其中的压力函数中的表皮系数进行拟合，使最后拟合目标油藏模型时的速度更快，从而使目标油藏模型与实测数据更接近。

结合第一方面第一种可能的实施方式，本申请实施例提供了第一方面的第四种可能的实施方式，其中：所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述第三目标数据段属于所述第二变化趋势信息中的一数据段，所述使用第三设定拟合方式对所述井储参数进行拟合，并使用所述第三目标数据段对拟合所述井储参数后对应的所述目标油藏模型进行验证，直到误差值小于第三设定值的步骤，包括：使用第三设定拟合方式确定出所述井储参数的当前值，形成当前的所述导数函数；使用所述第三目标数据段对所述导数函数进行验证；判断验证得到所述导数函数与所述第三目标数据段的误差值是否小于第三设定值；若所述误差值不小于所述第三设定值，则重复上述的步骤，直到所述误差值小于所述第三设定值。

进一步，本申请实施例提供的油藏数据拟合方法，上述的目标油藏模型包括压力函数及所述压力函数对应的导数函数，通过对其中的导数函数中的井储参数进行拟合，使最后拟合目标油藏模型时的速度更快，从而使目标油藏模型与实测数据更接近。

结合第一方面第一种可能的实施方式，本申请实施例提供了第一方面的第五种可能的实施方式，其中：所述压力相关数据包括所述压力数据的导数数据，所述根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点的步骤，包括：从所述导数数据中获取导数最大值点、平稳起点、平稳终点作为所述关键数据点；

所述从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段的步骤，包括：将所述第二变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第一目标数据段；将所述第一变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第二目标数据段；将所述第二变化趋势信息中对应所述导数最大值点至所述平稳起点之间的数据段作为所述第三目标数据段。

进一步，本申请实施例提供的油藏数据拟合方法，还可以通过将能够代表导数数据的特性的导数最大值点、平稳起点、平稳终点作为关键数据节点，由于井储参数的变化对第二变化趋势信息最大值点至所述平稳起点之间的数据段影响较明显，表皮系数的变化对第一变化趋势信息中对应所述平稳起点至所述平稳终点之间的影响较明显，渗透率对第二变化趋势信息中对应所述平稳起点至所述平稳终点之间影响明显，基于此确定出用作校准数据的目标数据段，使最后拟合目标油藏模型时的速度更快，从而使目标油藏模型与实测数据更接近。

结合第一方面第五种可能的实施方式，本申请实施例提供了第一方面的第六种可能的实施方式，其中：在所述将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型的步骤之前，所述方法还包括：任意设定所述目标参数的值，以得到当前的所述目标油藏模型；通过所述第二变化趋势信息中所述导数最大值点对应的数据点，对所述井储参数进行调整，以使所述目标油藏模型中的导数函数的最大值点与所述导数最大值点之间的误差小于第四设定值；将当前的所述目标参数的值作为拟合每个待确定的所述目标参数的初始值。

进一步，本申请实施例提供的油藏数据拟合方法，还可以先初步对井储参数进行调整，可以使具有确定的初始值的目标油藏模型中的导数函数的最大值点与所述导数最大值点距离接近，从而使后续步骤中再对渗透率进行拟合时选定的目标数据段偏差相对较小，从而可以更快地拟合出渗透率参数。

结合第一方面第一种至第六种可能的实施方式，本申请实施例提供了第一方面的第七种可能的实施方式，其中：所述根据多组实测油藏数据获取压力数据及压力相关数据的步骤，包括：从所述多组实测油藏数据确定出压力数据；使用移动窗口平滑法对所述压力数据进行计算，得到压力导数数据作为压力相关数据。

进一步，本申请实施例提供的油藏数据拟合方法，还可以使用拉氏变换反演算法计算出实测数据对应的压力数据及导数数据，可以使实测的数据与压力数据及导数数据对应，以用于后续拟合模型时候的校准数据。

第二方面，本申请实施例还提供一种油藏数据分析方法，包括：

使用上述的油藏数据拟合方法拟合得到的目标油藏模型确定目标气田的油藏动态数据。

本申请实施例提供的油藏数据分析方法，采用拟合出的目标油藏模型对目标气田进行数据分析，从而了解目标气田的状态。

第三方面，本申请实施例还提供一种油藏数据拟合装置，包括：

获取模块，用于根据多组实测油藏数据获取压力数据及压力相关数据；

确定模块，用于根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点；

分段模块，用于根据所述关键数据节点将所述压力数据形成的第一变化趋势信息及所述压力相关数据形成的第二变化趋势信息分成多段数据段；

匹配模块，用于从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段；

拟合模块，用于将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型。

第四方面，本申请实施例还提供一种油藏数据分析装置，包括：

确定模块，用于使用上述的油藏数据拟合方法拟合得到的目标油藏模型确定目标气田的油藏动态数据。

第五方面，本申请实施例还提供一种电子设备，包括：处理器、存储器和总线，所述存储器存储有所述处理器可执行的机器可读指令，当电子设备运行时，所述处理器与所述存储器之间通过总线通信，所述机器可读指令被所述处理器执行时执行上述第一方面或第二方面的方法的步骤。

第六方面，本申请实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述第一方面或第二方面中任一种可能的实施方式中的方法的步骤。

为使本申请的上述目的、特征和优点能更明显易懂，下文特举实施例，并配合所附附图，作详细说明如下。

附图说明

为了更清楚地说明本申请实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

图1为本申请实施例提供的电子设备的方框示意图。

图2为本申请实施例提供的油藏数据拟合方法的流程图。

图3为本申请实施例提供的压力数据和压力相关数据拟合成曲线的示意图。

图4为本申请实施例提供的曲线拟合过程的曲线示意图。

图5为本申请实施例提供的另一曲线拟合过程的曲线示意图。

图6为本申请实施例提供的再一曲线拟合过程的曲线示意图。

图7为本申请实施例提供的油藏数据拟合方法的步骤S205的详细流程图。

图8为本申请实施例提供的一参数的变化对目标油藏模型影响示意图。

图9为本申请实施例提供的另一参数的变化对目标油藏模型影响示意图。

图10为本申请实施例提供的再一参数的变化对目标油藏模型影响示意图。

图11为本申请实施例提供的其中一阶段的拟合参数过程的示意图。

图12为本申请实施例提供的油藏数据拟合装置的功能模块示意图。

具体实施方式

下面将结合本申请实施例中附图，对本申请实施例中的技术方案进行描述。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本申请的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

目前，试井分析利用瞬变压力分析储层特征，是油藏工程的重要分支。其主要技术要求是能够快速准确的估计出渗透率、表皮系数、井储系数等相关参数，为管理和预测油藏动态提供依据。早期试井分析一般采用手动图版拟合方法，近年来，随着计算机技术的进步，自动拟合技术被引入到试井分析过程，其主要思想是将问题视为数学意义下的非线性最小二乘问题。其求解方法主要分为两类：1.以牛顿法为基础的相关算法，如LM(Levenberg-Marquard，中文称：列文伯格-马夸尔特)算法、信赖域算法等，这类方法属于局部搜索算法，对初值一般比较敏感，需要对初值有比较好的估算。2.基于智能算法的全局优化算法，如遗传算法、差异进化算法、粒子群优化算法等。这类算法可能存在收敛速度较慢、算法需要耗费大量时间等不足。

具体地，试井参数自动拟合算法的实质是寻找最优参数下的理论曲线和实测曲线的最佳拟合，使其偏差达到最小。因此，可以建立如下数学模型：

其中，E(α)表示目标函数，该目标函数是关于α的函数；α表示所求油藏参数组成的向量；y_i表示实测压力数据；x_i表示采样节点；n表示实测数据的数量；F(α,x_i)表示试井分析模型的响应，可以是α的非线性函数。其中，F(α,x_i)可以用于计算拟合的气田的压力数据。

试井参数自动拟合方法可以基于上述模型的，其主要问题是把试井参数的拟合视为一个纯数学问题进行求解。因此，上述的方法可能会存在以下缺点：1.计算复杂度大，算法需要大量的运算时间。2.采用牛顿类方法时，对初值要求高，具有局部收敛缺点，算法对初值极为敏感，经常无法获得合理的参数估计。3.采用智能算法时，收敛速度慢，因为需要计算大量F(α,x_i)值，使得算法的计算量很大，而且所获得的结果可能距离最优值较远。

基于上述研究，本申请的发明人进行进一步地研究，考虑试井曲线本身的特点考虑到拟合过程中。

基于此，本申请提供的一种油藏数据拟合方法、油藏数据分析方法及装置，可以通过对不同的参数拟合过程中使用不同的数据段进行校准，可以提高模型拟合的速度。

为便于对本实施例进行理解，首先对执行本申请实施例所公开的油藏数据拟合方法或油藏数据分析方法的电子设备进行详细介绍。

实施例一

如图1所示，是电子设备的方框示意图。本实施例中的电子设备包括：存储器111、存储控制器112、处理器113、外设接口114、输入输出单元115、显示单元116。本领域普通技术人员可以理解，图1所示的结构仅为示意，其并不对电子设备100的结构造成限定。例如，电子设备100还可包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。

上述的存储器111、存储控制器112、处理器113、外设接口114、输入输出单元115及显示单元116各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。上述的处理器113用于执行存储器中存储的可执行模块。

其中，存储器111用于存储计算机程序，所述处理器113在接收到执行指令后，执行所述程序，本申请实施例任一实施例揭示的过程定义的电子设备100所执行的方法可以应用于处理器113中，或者由处理器113实现。

处理器113可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

下面通过方法实施例对上述的电子设备100可以执行的方法步骤进行详细描述。

实施例二

请参阅图2，是本申请实施例提供的油藏数据拟合方法的流程图。下面将对图2所示的具体流程进行详细阐述。

步骤S201，根据多组实测油藏数据获取压力数据及压力相关数据。

上述的实测油藏数据可以包括：压力数据、流量数据、粘度、渗透率、厚度、井半径、井储系数、孔隙度等数据。

在一个实例中，上述的实测油藏数据中的压力数据可以由油气田井的井底的压力计采集得到。一般油气井打开后，可以放一个永久压力计到井底，井底的压力数据会以固定间隔传输给一电子设备。可选地，可以将井底的压力数据描绘成压力曲线。由于会有关井和开井操作，开井压力下降，关井压力上升，因此，井底的压力数据对应的压力曲线呈分段状态。其中，拟合用的压力数据可以是一段关井状态下的压力数据；或，拟合用的压力数据可以是一段开井状态下的压力数据。

在一个可选的实施方式中：步骤S201可以被实现为：对压力数据进行处理，得到压力相关数据。

下面就得到压力相关数据的相关方法进行介绍。

在实际的试井分析中，为了使用典型曲线进行拟合，需要对实测压力差ΔP求导，导数公式可以表示为：

若/>

存在误差，则/>

的误差将增大到Δt倍。因此，当Δt较大时，可以选用一种较好的求导方法来求实测压力数据的压力导数，从而减少由此造成的误差。

求压力导数的常用方法有：三点差分加权法，移动窗口平滑法[25,26]。

三点差分加权法的公式为：

其中，1和2分别是第i点的前后相邻点。X表示时间的函数，对压降数据，X＝ln(Δt)。对压恢数据，还需在求导后做多流量修正。

移动窗口平滑法是三点法的推广，通过增加点1，2与当前点i的距离L，可减少压力的相对变差，达到平滑的目的。窗口长度2L与数据全长的比例称为“光滑系数”，记为w，表示相对于时间对数ln(Δt)的比例。w的取值范围可以为0～0.5，在一个实例中，可以令w＝0.2，当w＝0时，即为三点加权法。试井分析软件可以定义“光滑系数w”为左右半窗在对数坐标下的长度，其中w＝0.1。

对左右两端的抽样点，若左右半窗长度可能小于给定的光滑系数，则可以缩减半窗长度。例如，当左端的点左边半窗长度不足时，缩减左半窗口长度，右半窗取同样的比例。当右端的点右边半窗长度不足时，缩减右半窗口长度，左半窗取同样的比例。而对首末两个端点，左半窗或右半窗长度为0，其压力导数可以使用三点拉格朗日插值求导，若取值小于0，则使用右半窗或左半窗的导数代替。

步骤S202，根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点。

在一个可选的实施方式中，从所述导数数据中获取导数最大值点、平稳起点、平稳终点作为所述关键数据点。

示例性的，可以将上述的压力相关数据可以是压力数据对应的导数数据。可选地，可以将导数数据拟合成曲线。如图3所示，示出了将压力数据和压力相关数据拟合成曲线的示意图。其中，图3所示的两条曲线中上边的曲线表示压力数据拟合的压力曲线，下边的曲线表示压力相关数据拟合的导数曲线。导数曲线上的三个圈出的点分别表示：导数最大值点、平稳起点、平稳终点。

步骤S203，根据所述关键数据节点将所述压力数据形成的第一变化趋势信息及所述压力相关数据形成的第二变化趋势信息分成多段数据段。

在一实施方式中，上述的第一变化趋势信息可以是压力曲线，上述的第二变化趋势信息可以表示导数曲线。示例性的，上述的第一变化趋势信息及第二变化趋势信息可以表示其它能够表示压力数据及压力相关数据态势的信息，例如，离散点图、三维实景拟合图等。

可再次参考图3所示的示意图，通过上述的导数最大值点、平稳起点、平稳终点三个点可以将导数曲线划分成三个数据段，分别是：导数曲线上升至导数最大值点导数曲线段对应的数据段；导数最大值点下降至平稳起点导数曲线段对应的数据段；平稳起点至平稳终点导数曲线段对应的数据段。通过上述的导数最大值点、平稳起点、平稳终点三个点可以将压力曲线划分成三个数据段，分别是：压力曲线上升至导数最大值点的横坐标位置的压力曲线段对应的数据段；导数最大值点的横坐标处下降至平稳起点的横坐标位置处的压力曲线段对应的数据段；平稳起点的横坐标处至平稳终点的横坐标处的压力曲线段对应的数据段。

步骤S204，从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段。

可选地，可以根据拉氏变换反演算法得到待拟合的目标油藏模型。

本申请实施例可以使用如下模型对应的油井的参数拟合问题，井筒储集和表皮效应影响的均质无限大油藏无因次数学模型表示为：

在一可选的实例中，上述的参数：p_D、t_D、r_D、C_D可以表示为：

其中，p_D表示无因次压力；t_D表示无因次时间；r_D表示无因次井半径；C_D表示无因次井储系数；p_wD表示井底压力；S表示表皮系数；q表示流量；p表示压力；μ表示粘度；K表示渗透率；h表示厚度；r井半径；t表示时间；C表示井储系数；φ表示孔隙度；r_w表示井底半径。

拉氏空间上的井底压力

的表示式为：

其中，K₀、K₁分别表示零阶虚宗量贝塞尔(Bessel)函数、一阶虚宗量Bessel函数；S表示表皮系数；μ表示粘度；

表示拉氏空间上的井底压力。

采用Gaver–Stehfest(中文称：盖弗-施泰福斯特，英文简称：G-S)中的拉氏变换反演算法获得实空间解，也就是拟合过程中的响应函数F(α,x_i)，也就是上述的待拟合的目标油藏模型。

首先，通过一些实例调整各个目标参数分参数对曲线影响的分析。

在一个实例中，均质无限大模型可以拟合渗透率(K)、表皮系数(S)、井储参数(C)三个参数，下面图4至图6分别示出了三个参数变化对压力曲线和相应导数曲线的影响。其中，示意图中的上面一组曲线表示随参数变化对应的压力曲线，示意图中下面一组曲线表示随参数变化对应的导数曲线。

如图4所示，参数C对压力曲线的上升段，及压力曲线的上升段对应在导数曲线的区域影响显著。随着C的变化图像是左右平移效果，图像随着参数C的增大向右平移。如图5所示，参数S对压力曲线的平稳段影响显著，随着参数S的变化压力曲线是上下平移效果，随着参数S的增大，压力曲线平稳段曲线上移。如图6所示，参数K对压力曲线及导数曲线的平稳段影响显著，随着参数K变化压力曲线及导数曲线是上下平移效果，随着参数K的增大压力曲线及导数曲线下移。

基于上述的各个目标参数对压力曲线和导数曲线的影响研究，上述的步骤S204可以被实施为：将第二变化趋势信息中对应平稳起点至平稳终点之间的数据段作为所述第一目标数据段；将所述第一变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第二目标数据段；将所述第二变化趋势信息中对应所述导数最大值点至所述平稳起点之间的数据段作为所述第三目标数据段。

步骤S205，将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型。

其中，上述的目标参数可以包括：井储参数、表皮系数及渗透率。与目标参数对应，上述的目标数据段也可以包括：用于验证渗透率的第一目标数据段、用于验证表皮系数的第二目标数据段及用于验证井储参数的第三目标数据段。

如图7所示，上述的步骤S205可以包括以下步骤。

步骤S2051，使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到误差值小于第一设定值。

其中，目标油藏模型包括压力函数及所述压力函数对应的导数函数。上述的第一目标数据段属于所述第二变化趋势信息中的一数据段。

可选地，步骤S2051可以包括：使用第一设定拟合方式确定出所述渗透率的当前值，形成当前的所述导数函数；使用所述第一目标数据段对所述导数函数进行验证；判断验证得到的所述导数函数与所述第一目标数据段的误差值是否小于第一设定值；若所述误差值不小于所述第一设定值，则重复上述的步骤，直到所述误差值小于所述第一设定值。

当导数函数与第二目标数据段的误差值小于第一设定值时，对应的渗透率的值，可以作为一次拟合过程中渗透率的拟合结果。

上述的导数函数与所述第二目标数据段的误差值的计算方式可以是：计算导数函数与第二目标数据段中每个数据的差值，然后计算所有数据差值的平均值、方差等能表示误差的数值作为上述的误差值。

可选地，上述的第一设定值可以是一较小的数值，具体可以按照精度要求选择合适的数值，精度要求高时，第一设定值可以选择较小数值，精度要求低时，第一设定值可以选择较大数值。例如，上述的第一设定值可以是0.1、0.01、0.03、0.05等数值。

上述的第一设定拟合方式可以是二分拟合算法等。

在一个可选的实施方式中，可以将第二变化趋势信息中对应平稳起点至平稳终点之间的数据段作为所述第一目标数据段，以作为参数K的拟合校准数据，采用二分方式对参数K进行拟合。如图8所示，图示中多条上下叠置的线显示了随着参数K的变化平稳段导数拟合曲线的变化过程。可以选择拟合曲线与真解导数曲线基本重叠时，对应的参数K值作为参数K的当前拟合结果。

可选地，当上述的第一设定拟合方式是二分拟合算法时，可以计算最后两次的导数函数之间的误差，作为是否退出计算的标准，当最后两次的导数函数之间的误差小于第一设定值时，则可以退出计算，得到渗透率的当前值。其中，最后两次计算得到的渗透率的当前值，形成最后两次的导数函数。

步骤S2052，使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到误差值小于第二设定值。

上述的第二目标数据段属于所述第一变化趋势信息中的一数据段。

可选地，步骤S2052可以包括：使用第二设定拟合方式确定出所述表皮系数的当前值，形成当前的所述压力函数；使用所述第二目标数据段对所述压力函数进行验证；判断验证得到所述压力函数与所述第二目标数据段的误差值是否小于第二设定值；若所述误差值不小于所述第二设定值，则重复上述的步骤，直到所述误差值小于所述第二设定值。

当上述的压力函数与第二目标数据段的误差值小于上述的第二设定值时，对应的表皮系数的值，可以作为一次拟合过程中表皮系数的拟合结果。

可选地，上述的第二设定值可以是一较小的数值，例如，上述的第二设定值可以是0.1、0.01、0.03、0.05、0.15等数值。

上述的第二设定拟合方式可以是与上述的第一设定拟合方式相同的算法，也可以是不同的算法。例如，上述的第二设定拟合方式也可以是二分拟合算法。

在一个可选的实施方式中，将所述第一变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第二目标数据段，以作为参数S的校准数据，以参数S进行校准。如图9所示，上述的压力曲线对应的水平段处上下叠置的多条线表示随着参数S的变化平稳段拟合曲线的变化过程。可以选择拟合曲线与真解曲线基本重叠时，对应的参数S值作为参数S的当前拟合结果。

可选地，当上述的第二设定拟合方式是二分拟合算法时，可以计算最后两次的导数函数之间的误差，作为是否退出计算的标准，当最后两次的导数函数之间的误差小于第二设定值时，则可以退出计算，得到表皮系数的当前值。其中，最后两次计算得到的表皮系数的当前值，形成最后两次的导数函数。

步骤S2053，使用第三设定拟合方式对所述井储参数进行拟合，并使用所述第三目标数据段对拟合所述井储参数后对应的所述目标油藏模型进行验证，直到误差值小于第三设定值。

可选地，上述的第三设定值可以是一较小的数值，例如，上述的第三设定值可以是0.1、0.01、0.03、0.05、0.2等数值。

上述的第三设定拟合方式可以是与上述的第一设定拟合方式和第二设定拟合方式相同的算法，也可以是不同的算法。例如，上述的第三设定拟合方式也可以是二分拟合算法。

上述的第三目标数据段属于所述第二变化趋势信息中的一数据段。

可选地，上述的步骤S2053可以包括：使用第三设定拟合方式确定出所述井储参数的当前值，形成当前的所述导数函数；使用所述第三目标数据段对所述导数函数进行验证；判断验证得到所述导数函数与所述第三目标数据段的误差值是否小于第三设定值；若所述误差值不小于所述第三设定值，则重复上述的步骤，直到所述误差值小于所述第三设定值。

当上述的导数函数与第三目标数据段的误差值小于上述的第三设定值时，对应的井储参数的值，可以作为一次拟合过程中井储参数的拟合结果。

可选地，当上述的第三设定拟合方式是二分拟合算法时，可以计算最后两次的导数函数之间的误差，作为是否退出计算的标准，当最后两次的导数函数之间的误差小于第二设定值时，则可以退出计算，得到井储参数的当前值。其中，最后两次计算得到的井储参数的当前值，形成最后两次的导数函数。

在一实施方式中，将第二变化趋势信息中对应导数最大值点至平稳起点之间的数据段作为所述第三目标数据段，以用于对参数C进行拟合。如图10所示，上述的导数曲线对应的下降段处上下叠置的多条线表示随着参数C的变化导数曲线对应的下降段拟合曲线的变化过程。可以选择拟合曲线与真解曲线基本重叠时，对应的参数C值作为参数C的当前拟合结果。

可选地，也可以使用第二变化趋势信息中的起始数据点至导数最大值点之间的数据段作为参数C的校准数据。

步骤S2054，判断目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差是否小于目标设定值。

上述的目标设定值也可以是一较小的数值，例如，0.1、0.01、0.03、0.001、0.003等数值。

可选地，拟合出参数K、参数S及参数C后，可确定出目标油藏模型。可随机选取目标油藏模型中的数据点与压力数据及压力相关数据计算误差平均值或误差方差作为目标油藏模型与压力数据及压力相关数据的误差。

若目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差大于目标设定值，则返回步骤S2051；若目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差小于目标设定值，则流程结束，输出目标油藏模型。

在一个实例中，图8、9、10所示的拟合过程的示意图，采用的真解为C＝5，S＝2，K＝300。随机选择拟合的初值：C＝20，S＝50，K＝800，二分拟合过程参数的控制误差为1e^-4，迭代相对误差控制阈值为1％。拟合的最终参数值为C＝4.9396，S＝1.9935，K＝299.8439，与真解C＝5，S＝2，K＝300相比较数值误差较小。

在其它可选的实施方式中，在对参数K进行拟合之前，可以先进行参数的初步调整，可以使参数K的拟合过程中，选定数据段的位置不会有太大偏差。在步骤S205之前，油藏数据拟合方法方法还可以包括：任意设定所述目标参数的值，以得到当前的所述目标油藏模型；通过所述第二变化趋势信息中所述导数最大值点对应的数据点，对所述井储参数进行调整，以使所述目标油藏模型中的导数函数的最大值点与所述导数最大值点之间的误差小于第四设定值；将当前的所述目标参数的值作为拟合每个待确定的所述目标参数的初始值。

在一个实例中，任意设定的目标参数的值可以是：C＝20，S＝50，K＝800。

上述的第四设定值可以是一较小的值，例如，0.1、0.01、0.02、0.001、0.0001等数值。

给定拟合参数的任意初值之后，首先通过导数最大值点的位置坐标，对参数C进行二分调整，使用导数最大值点与目标油藏模型中的导数函数的最大值点距离差小于第四设定值。在一个实例中，可以使导数最大值点与目标油藏模型中的导数函数的最大值点重合。如图11所示，图示中的最上边的一条曲线表示实测压力曲线，下面的多条曲线分别表示拟合曲线及实测导数曲线。随着参数C的变化导数曲线的随着变化，最终拟合的导数曲线与实测的导数曲线的导数最大值点重叠。

本申请实施例提供的油藏数据拟合方法，还可以先初步对井储参数进行调整，可以使具有确定的初始值的目标油藏模型中的导数函数的最大值点与所述导数最大值点距离接近，从而使后续步骤中再对渗透率进行拟合时选定的目标数据段偏差相对较小，从而可以更快地拟合出渗透率参数。

本本申请实施例提供的油藏数据拟合方法，可以至少包括以下有益效果：1.由于本申请是从试井曲线的特征出发，深入分析了各参数对曲线形态的影响，指出了各参数的相对独立性，通过曲线的分段，分部对参数进行拟合。相对于现有技术中将参数拟合直接转化成纯数学问题，从非线性最小二乘问题求解方法入手来实现参数拟合而言，算法具有更好的强健性，对初值不敏感。2.各目标参数的拟合可选地采用二分法进行，二分法的执行效率高，而且对曲线进行分段后，我们不必把分段内所有的节点列入计算，这可以大幅度提高算法的执行效率。3.本申请中的算法的停止准则直接是针对数据误差的，这避免了智能算法会出现距离真解有较大偏差问题。

实施例三

本申请实施例提供一种油藏数据分析方法，本实施例中的方法包括：使用实施例二中的方法拟合得到的目标油藏模型确定目标气田的油藏动态数据。可选地，可以将采样节点数据输入拟合得到的目标油藏模型进行计算可以得到油藏动态数据。上述的油藏动态数据可以表示目标气田的压力。

实施例四

基于同一申请构思，本申请实施例中还提供了与油藏数据拟合方法对应的油藏数据拟合装置，由于本申请实施例中的装置解决问题的原理与本申请实施例上述油藏数据拟合方法相似，因此装置的实施可以参见方法的实施，重复之处不再赘述。

请参阅图12，是本申请实施例提供的油藏数据拟合装置的功能模块示意图。油藏数据拟合装置包括：获取模块301、确定模块302、分段模块303、匹配模块304，以及拟合模块305；其中：

获取模块301，用于根据多组实测油藏数据获取压力数据及压力相关数据；

确定模块302，用于根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点；

分段模块303，用于根据所述关键数据节点将所述压力数据形成的第一变化趋势信息及所述压力相关数据形成的第二变化趋势信息分成多段数据段；

匹配模块304，用于从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段；

拟合模块305，用于将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型。

一种可能的实施方式中，所述目标参数包括：渗透率表皮系数及井储参数，所述目标数据段包括用于验证所述渗透率的第一目标数据段、用于验证所述表皮系数的第二目标数据段及用于验证所述井储参数的第三目标数据段；拟合模块305，还用于：使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到误差值小于第一设定值；使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到误差值小于第二设定值；使用第三设定拟合方式对井储参数进行拟合，并使用所述第三目标数据段对拟合井储参数后对应的所述目标油藏模型进行验证，直到误差值小于第三设定值；重复上述的步骤，直到所述目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差小于目标设定值，确定出目标参数，以得到所述目标油藏模型。

一种可能的实施方式中，所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述第一目标数据段属于所述第二变化趋势信息中的一数据段；所述拟合模块305，还用于：使用第一设定拟合方式确定出所述渗透率的当前值，形成当前的所述导数函数；使用所述第一目标数据段对所述导数函数进行验证；判断验证得到的所述导数函数与所述第一目标数据段的误差值是否小于第一设定值；若所述误差值不小于所述第一设定值，则重复上述的步骤，直到所述误差值小于所述第一设定值。

一种可能的实施方式中，所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述第二目标数据段属于所述第一变化趋势信息中的一数据段，上述的拟合模块305，还可以用于：使用第二设定拟合方式确定出所述表皮系数的当前值，形成当前的所述压力函数；使用所述第二目标数据段对所述压力函数进行验证；判断验证得到所述压力函数与所述第二目标数据段的误差值是否小于第二设定值；若所述误差值不小于所述第二设定值，则重复上述的步骤，直到所述误差值小于所述第二设定值。

一种可能的实施方式中，所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述第三目标数据段属于所述第二变化趋势信息中的一数据段，上述的拟合模块305，还可以用于：使用第三设定拟合方式确定出所述井储参数的当前值，形成当前的所述导数函数；使用所述第三目标数据段对所述导数函数进行验证；判断验证得到所述导数函数与所述第三目标数据段的误差值是否小于第三设定值；若所述误差值不小于所述第三设定值，则重复上述的步骤，直到所述误差值小于所述第三设定值。

一种可能的实施方式中，所述压力相关数据包括所述压力数据的导数数据，上述的确定模块302，还可以用于：从所述导数数据中获取导数最大值点、平稳起点、平稳终点作为所述关键数据点；

上述的匹配模块304，还可以用于：将所述第二变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第一目标数据段；将所述第一变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第二目标数据段；将所述第二变化趋势信息中对应所述导数最大值点至所述平稳起点之间的数据段作为所述第三目标数据段。

一种可能的实施方式中，上述的油藏数据拟合装置还可以包括：校准模块306，用于：任意设定所述目标参数的值，以得到当前的所述目标油藏模型；通过所述第二变化趋势信息中所述导数最大值点对应的数据点，对所述井储参数进行调整，以使所述目标油藏模型中的导数函数的最大值点与所述导数最大值点之间的误差小于第四设定值；将当前的所述目标参数的值作为拟合每个待确定的所述目标参数的初始值。

一种可能的实施方式中，上述的获取模块301，还可以用于：从所述多组实测油藏数据确定出压力数据；

使用移动窗口平滑法对所述压力数据进行计算，得到压力导数数据作为压力相关数据。

此外，本申请实施例还提供一种计算机可读存储介质，该计算机可读存储介质上存储有计算机程序，该计算机程序被处理器运行时执行上述方法实施例中所述的油藏数据拟合方法方法或油藏数据分析方法的步骤。

本申请实施例所提供的油藏数据拟合方法的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行上述方法实施例中所述的油藏数据拟合方法的步骤，具体可参见上述方法实施例，在此不再赘述。

本申请实施例所提供的油藏数据分析方法的计算机程序产品，包括存储了程序代码的计算机可读存储介质，所述程序代码包括的指令可用于执行上述方法实施例中所述的油藏数据分析方法的步骤，具体可参见上述方法实施例，在此不再赘述。

在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

另外，在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。

所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。

以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以权利要求的保护范围为准。

Claims (9)

1.一种油藏数据拟合方法，其特征在于，包括：

根据多组实测油藏数据获取压力数据及压力相关数据；

根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点；

根据所述关键数据节点将所述压力数据形成的第一变化趋势信息及所述压力相关数据形成的第二变化趋势信息分成多段数据段；

从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段；

将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型；

所述目标参数包括：渗透率、表皮系数及井储参数，所述目标数据段包括用于验证所述渗透率的第一目标数据段、用于验证所述表皮系数的第二目标数据段及用于验证所述井储参数的第三目标数据段；所述将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型的步骤，包括：

使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到第一误差值小于第一设定值；

使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到第二误差值小于第二设定值；

使用第三设定拟合方式对所述井储参数进行拟合，并使用所述第三目标数据段对拟合所述井储参数后对应的所述目标油藏模型进行验证，直到第三误差值小于第三设定值；

重复上述对所述目标参数进行拟合，并使用对应的所述目标数据段对所述目标油藏模型进行验证，直到相应的误差值小于相应的设定值的步骤，直到所述目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差小于目标设定值，确定出目标参数，以得到所述目标油藏模型；

其中，所述第一目标数据段属于所述第二变化趋势信息中的一数据段；所述第二目标数据段属于所述第一变化趋势信息中的一数据；所述第三目标数据段属于所述第二变化趋势信息中的一数据。

2.根据权利要求1所述的方法，其特征在于，所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到第一误差值小于第一设定值的步骤，包括：

使用第一设定拟合方式确定出所述渗透率的当前值，形成当前的所述导数函数；

使用所述第一目标数据段对所述导数函数进行验证；

判断验证得到的所述导数函数与所述第一目标数据段的第一误差值是否小于第一设定值；

若所述第一误差值不小于所述第一设定值，则重复上述的步骤，直到所述第一误差值小于所述第一设定值。

3.根据权利要求1所述的方法，其特征在于，所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到第二误差值小于第二设定值的步骤，包括：

使用第二设定拟合方式确定出所述表皮系数的当前值，形成当前的所述压力函数；

使用所述第二目标数据段对所述压力函数进行验证；

判断验证得到所述压力函数与所述第二目标数据段的第二误差值是否小于第二设定值；

若所述第二误差值不小于所述第二设定值，则重复上述的步骤，直到所述第二误差值小于所述第二设定值。

4.根据权利要求1所述的方法，其特征在于，所述目标油藏模型包括压力函数及所述压力函数对应的导数函数；所述使用第三设定拟合方式对所述井储参数进行拟合，并使用所述第三目标数据段对拟合所述井储参数后对应的所述目标油藏模型进行验证，直到第三误差值小于第三设定值的步骤，包括：

使用第三设定拟合方式确定出所述井储参数的当前值，形成当前的所述导数函数；

使用所述第三目标数据段对所述导数函数进行验证；

判断验证得到所述导数函数与所述第三目标数据段的第三误差值是否小于第三设定值；

若所述第三误差值不小于所述第三设定值，则重复上述的步骤，直到所述第三误差值小于所述第三设定值。

5.根据权利要求1所述的方法，其特征在于，所述压力相关数据包括所述压力数据的导数数据，所述根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点的步骤，包括：

从所述导数数据中获取导数最大值点、平稳起点、平稳终点作为关键数据点；

所述从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段的步骤，包括：

将所述第二变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第一目标数据段；

将所述第一变化趋势信息中对应所述平稳起点至所述平稳终点之间的数据段作为所述第二目标数据段；

将所述第二变化趋势信息中对应所述导数最大值点至所述平稳起点之间的数据段作为所述第三目标数据段。

6.根据权利要求5所述的方法，其特征在于，在所述将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型的步骤之前，所述方法还包括：

任意设定所述目标参数的值，以得到当前的所述目标油藏模型；

通过所述第二变化趋势信息中所述导数最大值点对应的数据点，对所述井储参数进行调整，以使所述目标油藏模型中的导数函数的最大值点与所述导数最大值点之间的误差小于第四设定值；

将当前的所述目标参数的值作为拟合每个待确定的所述目标参数的初始值。

7.根据权利要求1-6任意一项所述的方法，其特征在于，所述根据多组实测油藏数据获取压力数据及压力相关数据的步骤，包括：

从所述多组实测油藏数据确定出压力数据；

使用移动窗口平滑法对所述压力数据进行计算，得到压力导数数据作为压力相关数据。

8.一种油藏数据分析方法，其特征在于，包括：

使用权利要求1-7任意一项所述的方法拟合得到的目标油藏模型确定目标气田的油藏动态数据。

9.一种油藏数据拟合装置，其特征在于，包括：

获取模块，用于根据多组实测油藏数据获取压力数据及压力相关数据；

确定模块，用于根据所述压力相关数据的变化趋势确定出所述压力相关数据中的关键数据节点；

分段模块，用于根据所述关键数据节点将所述压力数据形成的第一变化趋势信息及所述压力相关数据形成的第二变化趋势信息分成多段数据段；

匹配模块，用于从所述多段数据段中为待拟合的目标油藏模型中的每个待确定的目标参数匹配对应的目标数据段；

拟合模块，用于将每个所述目标参数对应的目标数据段作为校准数据，拟合出每个待确定的目标参数，以确定出所述目标油藏模型；

所述目标参数包括：渗透率表皮系数及井储参数，所述目标数据段包括用于验证所述渗透率的第一目标数据段、用于验证所述表皮系数的第二目标数据段及用于验证所述井储参数的第三目标数据段；所述拟合模块，还用于：使用第一设定拟合方式对所述渗透率进行拟合，并使用所述第一目标数据段对拟合所述渗透率后对应的所述目标油藏模型进行验证，直到第一误差值小于第一设定值；使用第二设定拟合方式对所述表皮系数进行拟合，并使用所述第二目标数据段对拟合所述表皮系数后对应的所述目标油藏模型进行验证，直到第二误差值小于第二设定值；使用第三设定拟合方式对所述井储参数进行拟合，并使用所述第三目标数据段对拟合所述井储参数后对应的所述目标油藏模型进行验证，直到第三误差值小于第三设定值；重复上述对所述目标参数进行拟合，并使用对应的所述目标数据段对所述目标油藏模型进行验证，直到相应的误差值小于相应的设定值的步骤，直到所述目标油藏模型的计算数值与所述压力数据及所述压力相关数据的误差小于目标设定值，确定出目标参数，以得到所述目标油藏模型；

其中，所述第一目标数据段属于所述第二变化趋势信息中的一数据段；所述第二目标数据段属于所述第一变化趋势信息中的一数据；所述第三目标数据段属于所述第二变化趋势信息中的一数据。

Images