技术领域
[0001] 本
发明涉及
天然气开发技术领域，特别涉及一种确定气井最终累积产量的方法及装置。
背景技术
[0002] 气井最终累积产量(EUR)即在废弃条件约束下气井由测试到废弃的过程中最终可采出的天然气总量，它是表征气井生产动态特征、评价气井开发效果的关键技术指标之一，准确、快速地预测该项指标对指导气井/气田的经济效益评价以及气田开发方案编制/调整具有重要意义。
[0003] 目前，气井EUR的评价方法众多，主要包括经验公式法、理论
解析法、数值模拟法等。经验公式法以气井产量数据为分析对象，通过数据最优拟合来建立气井产量表达式，进而预测气井EUR，该方法以Arps产量递减模型为典型代表，该经典模型在常规油气井EUR预测中的应用非常成功，而后，随着致密气、
页岩气等非常规资源的兴起，经验公式法得到了较大的发展，先后发展出了Valkó扩展指数模型、Ilk幂律指数模型、Duong模型等。理论解析法根据分析原理的不同，又可分为流动物质平衡法、典型曲线拟合法、线性流分析法、解析模型法等；流动物质平衡法是利用气井流动达到拟稳态时储
层压力匀速降落的特点，建立气井累积产量(Gp)与压力参数(Pwf/Z)的线性关系来预测气井最终累积产量；典型曲线拟合法是借鉴试井解释的思想，将生产数据(压力、产量)导入模型，并与理论曲线拟合，获得相应的动态解释参数，进而预测气井未来产量变化；线性流分析法是针对特定流动阶段内的动态数据，开展双线性流、线性流等分析，以获得裂缝、储层特性参数，作为未来产量预测的理论
基础；解析模型法是建立以分段压裂
水平井为对象的概念模型，通过全生命周期的生产数据获得确定性的动态参数，进而预测气井未来产量及EUR。数值模拟法是基于流动过程的生产动态模拟，通过物理-数学模型与实际数据的历史拟合，厘清影响生产的主控因素，掌握气井的生产动态规律，从而获得可靠的EUR，该方法计算量大，应用相对复杂。
[0004] 上述经验公式法和理论解析法的应用都是针对单井开展的，对于一个区
块或气田而言，随着开发
进程的推进，气井数量会快速增加，针对数百口甚至数千口的气井动态资料，如果逐井开展EUR评价，将会带来巨大的工作量和实施难度；数值模拟法虽可批量预测区块气井EUR，但由于模型建立、历史拟合等环节的工作量相对较大，且耗时较长，使得该方法应用同样存在不便。因而开展EUR批量快速评价方法研究具有重要的现实意义。
发明内容
[0005] 本发明
实施例提供了一种确定气井最终累积产量的方法，以解决
现有技术中在批量预测区块气井EUR时存在工作量大、耗时长的技术问题。该方法包括：
[0006] 对目标区域内的气井进行分类，不同类的气井满足不同的开采效果标准；
[0007] 确定所述目标区域内每口气井的产量递减的初始时间点；
[0008] 确定所述目标区域内每口气井在所述产量递减的初始时间点之前的累计产量，该累计产量为第一累计产量；
[0009] 确定所述目标区域内每口气井由所述产量递减的初始时间点至产量递减期内预设时间点的累计产量，该累计产量为第二累计产量；
[0010] 计算每类气井在产量递减期内的产量系数；
[0011] 针对每类气井，根据该类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，计算该类气井中每口气井的最终累积产量。
[0012] 本发明实施例还提供了一种计算机设备，包括
存储器、处理器及存储在存储器上并可在处理器上运行的
计算机程序，所述处理器执行所述计算机程序时实现上述任意的确定气井最终累积产量的方法。
[0013] 本发明实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述任意的确定气井最终累积产量的方法的计算机程序。
[0014] 本发明实施例还提供了一种确定气井最终累积产量的装置，以解决现有技术中在批量预测区块气井EUR时存在工作量大、耗时长的技术问题。该装置包括：
[0015] 气井分类模块，用于对目标区域内的气井进行分类，不同类的气井满足不同的开采效果标准；
[0016] 时间点确定模块，用于确定所述目标区域内每口气井的产量递减的初始时间点；
[0017] 第一产量确定模块，用于确定所述目标区域内每口气井在所述产量递减的初始时间点之前的累计产量，该累计产量为第一累计产量；
[0018] 第二产量确定模块，用于确定所述目标区域内每口气井由所述产量递减的初始时间点至产量递减期内预设时间点的累计产量，该累计产量为第二累计产量；
[0019] 系数计算模块，用于计算每类气井在产量递减期内的产量系数；
[0020] 最终累积产量确定模块，用于针对每类气井，根据该类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，计算该类气井中每口气井的最终累积产量。
[0021] 在本发明实施例中，通过对目标区域内的气井进行分类，进而得到目标区域内每口气井的第一累计产量和第二累计产量，再确定出每类气井在产量递减期内的产量系数，最后，就可以根据每类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，批量计算出每类气井中每口气井的最终累积产量。该方案主要是确定出了每类气井在产量递减期内的产量系数，使得可以实现根据每类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，批量计算出每类气井中每口气井的最终累积产量，相对现有技术而言避免了模型建立、历史拟合等环节，有利于在批量预测区块气井EUR的同时，减少工作量、减少耗时。
附图说明
[0022] 此处所说明的附图用来提供对本发明的进一步理解，构成本
申请的一部分，并不构成对本发明的限定。在附图中：
[0023] 图1是本发明实施例提供的一种确定气井最终累积产量的方法的
流程图；
[0024] 图2是本发明实施例提供的一种确定气井最终累积产量的示例流程图；
[0025] 图3是本发明实施例提供的一种气井产量递减起始点识别及递减期内阶段累积产量统计示意图；
[0026] 图4是本发明实施例提供的一种区块A的I类气井“递减240天的累积产量”与“递减180天内的累积产量”之间的线性关系图；
[0027] 图5是本发明实施例提供的一种区块A的I类气井“递减330天的累积产量”与“递减180天内的累积产量”之间的线性关系图；
[0028] 图6是本发明实施例提供的一种区块A的I类气井“递减期内的阶段累积产量”与“递减180天内的累积产量”之间的线性关系系列，线性斜率随递减时间的延长而增大的示意图；
[0029] 图7是本发明实施例提供的一种区块A的I类气井“递减期内的阶段累积产量”与“递减180天内的累积产量”之间的线性关系斜率与递减时间之间的函数关系拟合示意图；
[0030] 图8是本发明实施例提供的一种确定气井最终累积产量的装置的结构
框图。
具体实施方式
[0031] 为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施方式和附图，对本发明做进一步详细说明。在此，本发明的示意性实施方式及其说明用于解释本发明，但并不作为对本发明的限定。
[0032] 在本发明实施例中，提供了一种确定气井最终累积产量的方法，如图1所示，该方法包括：
[0033] 步骤101：对目标区域内的气井进行分类，不同类的气井满足不同的开采效果标准；
[0034] 步骤102：确定所述目标区域内每口气井的产量递减的初始时间点；
[0035] 步骤103：确定所述目标区域内每口气井在所述产量递减的初始时间点之前的累计产量，该累计产量为第一累计产量；
[0036] 步骤104：确定所述目标区域内每口气井由所述产量递减的初始时间点至产量递减期内预设时间点的累计产量，该累计产量为第二累计产量；
[0037] 步骤105：计算每类气井在产量递减期内的产量系数；
[0038] 步骤106：针对每类气井，根据该类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，计算该类气井中每口气井的最终累积产量。
[0039] 由图1所示的流程可知，在本发明实施例中，通过对目标区域内的气井进行分类，进而得到目标区域内每口气井的第一累计产量和第二累计产量，再确定出每类气井在产量递减期内的产量系数，最后，就可以根据每类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，批量计算出每类气井中每口气井的最终累积产量。该方案主要是确定出了每类气井在产量递减期内的产量系数，使得可以实现根据每类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，批量计算出每类气井中每口气井的最终累积产量，相对现有技术而言避免了模型建立、历史拟合等环节，有利于在批量预测区块气井EUR的同时，减少工作量、减少耗时。
[0040] 具体实施时，上述目标区域可以是任意待预测气井最终累积产量的气井所在的区域，例如，可以是一个气井区块或气田等。
[0041] 具体实施时，在本实施例中，提出根据综合分类系数、生产测试时的测试产量、投产首年的平均日产量以及气井内部收益率，对目标区域内的气井进行分类。
[0042] 具体实施时，在对目标区域内的气井进行分类之前，进行基础数据准备工作，采集目标区域内所有气井的动、静态资料，包括地质特征参数、工程施工资料以及生产动态数据；检查、核实气井的动、静态数据的合理性，尤其对于气井生产动态资料，删除或校正异常的产量和压力数据。
[0043] 具体实施时，依据每口气井的特征参数、工程施工资料等关键参数构建综合分类系数，具体地，从影响气井产量的诸多因素中筛选优质储层钻遇比例、压裂水平段长度、压裂改造段数以及压裂加砂量，构建气井综合分类系数的公式如下：
[0044]
[0045] 其中，Ψ为气井的综合分类系数，无量纲；L为气井压裂水平段长度，单位为：m；N为改造段数；S为压裂加砂量，单位为：t；β为优质储层钻遇比例，单位为：％；脚标avg代表相应的样本平均值。
[0046] 具体实施时，以气井早期的生产动态资料为基础，优选关键评价指标进行气井分类。具体地，对于钻完井后、正式投产前，气井普遍进行生产测试，测试产量是可信度较强的第一手动态资料，初步较客观地反映了气井的开发效果；对于正式投产早期阶段，首年平均日产量可有效反映气井产能。因而，优选生产测试时的测试产量qtest和投产首年的平均日产量qavg,1st为分类指标。
[0047] 进一步，以综合分类系数、生产测试时的测试产量、投产首年的平均日产量以及气井内部收益率为分类指标，综合分类系数、生产测试时的测试产量、投产首年的平均日产量以及气井内部收益率的具体分类之间的界限值可以视情况而定，本申请不做具体限定。例如，综合分类系数可以分别以1.06和1.28为分类参考界限，生产测试时的测试产量可以分别以13.9 104m3/d和15.3 104m3/d为分类参考界限，投产首年的平均日产量可以分别以7.2 4 3 4 3
10m/d和8.0 10m/d为分类参考界限，气井内部收益率可以分别以8％和12％为分类参考界限，来确定气井分类标准。
[0048] 具体实施时，对气井进行分类后，以气井分类标准为指导，获得每类气井的样本集合，基于每类气井的样本集合计算每类气井在产量递减期内的产量系数。
[0049] 具体的，通过以下步骤计算每类气井在产量递减期内的产量系数：
[0050] 确定该类气井中每口气井由所述产量递减的初始时间点至产量递减期内任意时间点的累计产量，该累计产量为递减期内的阶段累积量，其中，相邻时间点之间相差预设时间步长；例如，时间步长可以为“月”、“日”等。
[0051] 针对该类气井中所有气井，以第二累计产量为横坐标，以不同时间点对应的递减期内的阶段累积量为纵坐标，绘制散点图，得到一系列的线性关系，对各散点进行线性拟合，得到与时间点对应的线性斜率；
[0052] 以每个时间点对应的线性斜率为样本数据，绘制线性斜率与时间点的关系图，根据数据点分布趋势优选函数进行曲线拟合，得到线性斜率与时间点的函数表达式，该线性斜率为该类气井在产量递减期内的产量系数，依据此函数表达式，可计算任意产量递减期内时间点对应的气井总累积产量。
[0053] 具体的，线性斜率与时间点的函数表达式如下：
[0054] Gp(t)＝Gp0+k(t)·Gptn
[0055] 其中，Gp(t)为产量递减期内任意时间点的气井总累积产量；Gp0为第一累计产量；Gptn为第二累计产量；k(t)为线性斜率；t为产量递减期内的时间点。
[0056] 具体实施时，第二累计产量Gptn对应的产量递减期内预设时间点可以就具体情况而定，本申请不做限定，例如，该产量递减期内预设时间点可以是产量递减期内第180日，产量递减期内预设时间点还可以是产量递减期内第5个月等。
[0057] 具体实施时，得到每类气井的产量系数后，可以通过以下公式实现根据该类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，计算该类气井中每口气井的最终累积产量：
[0058] EUR＝Gp0+k(tabandon)·Gptn
[0059] 其中，EUR为每口气井的最终累积产量；Gp0为第一累计产量；Gptn为第二累计产量；k(tabandon)为该类气井在产量递减期内的产量系数；tabandon为气井最终寿命对应的时间点。
[0060] 具体的，通过以下公式确定气井最终寿命对应的时间点：
[0061] 或tabandon≥tlimit
[0062] 其中，qabandon为气井经济极限产量；Gp(tabandon)为气井最终寿命对应的时间点的气井总累积产量；Gp(tabandon-1)为气井最终寿命对应时间点的前一时间点的气井总累积产量，相邻时间点之间的时间步长为一个月；tlimit为气井最终寿命。
[0063] 以下基于具体示例来描述上述确定气井最终累积产量的方法，如图2所示，包括以下步骤：
[0064] 步骤1，采集区块A的页岩气井动、静态资料：气井动态数据主要包括单井日产气量qg、套压Pc、测试产量qtest、测试压力Ptest等，如图3所示，带平滑线的实心圆散点即单井日产气量；静态资料主要包括气井的优质储层钻遇比例β、压裂水平段长度L、压裂改造段数N、压裂加砂量S等，如表1所示。检查数据资料的合理性，尤其对于气井生产动态资料，删除或校正异常的产量和压力数据。
[0065] 表1
[0066]
[0067]
[0068] 步骤2，以开发效果为评价目标，对区块A的页岩气井进行分类。首先，以地质、工程资料为基础求取气井分类评价指标：基于每口气井的压裂水平段长度L、压裂改造段数N、压裂加砂量S和优质储层钻遇比例β，求取气井综合分类系数Ψ；其次，以气井早期的生产动态资料为基础求取气井分类评价指标：区块A气井测试产量qtest数据较少，因而只选用首年平均日产qavg,1st作为气井分类评价指标；具体数据如表1所示。
[0069] 依据区块A已经建立的分类评价标准(表2)对气井进行分类，获得每类气井的样本集，分类结果如表1所示。
[0070] 表2
[0071]
[0072] 步骤3，针对每类气井，逐井开展气井产量递减识别，图3中的实心三
角形即标记出了气井产量递减的起始点。
[0073] 步骤4，针对每类气井，计算每口气井在产量递减现象出现前的累积产量Gp0(即上述第一累计产量)，如图3中的实心正方形所示。
[0074] 步骤5，针对每类气井的产量递减期，自产量递减始点(t＝0)起，以“月”为时间步长，逐井统计不同递减时间点t对应的“递减期内的阶段累积产量Gp(t)”，如图3中的实心星形所示。
[0075] 步骤6，针对每类气井，将步骤5中统计数据进行汇总，如表3所示，基于这些数据，以“递减180天内的累积产量Gp180(即上述第二累计产量)”为横坐标，以不同递减时间点t对应的“递减期内的阶段累积产量Gp(t)”为纵坐标，绘制散点图，得到一系列的线性关系，对散点进行线性拟合并求取线性斜率K(t)，本实施例的图4～图6举例展示了I类气井的线性关系特征，后文的论述均只针对I类气井展开。
[0076] 表3
[0077]
[0078]
[0079] 步骤7：以每个时间点t对应的线性斜率K(t)为样本数据，如表4所示，绘制线性斜率K(t)与递减时间点t的关系图，如图7所示，根据数据点分布趋势优选函数进行曲线拟合，明确线性斜率K(t)与递减时间t之间的函数表达式：
[0080] K(t)＝1.292×(6.602t+23.756)0.2852-3.215
[0081] 依据此函数表达式，可计算任意递减时间点t对应的气井总累积产量Gp(t)。
[0082] 表4
[0083]
[0084]
[0085] 步骤8：以气井的“经济极限产量qabandon”取0.1万方/天和“根据油气田实际需求所设定的气井极限生产时间tlimit”取20年为例求取气井最终累积产量(EUR)，则约束条件如下：
[0086] 或t≥20
[0087] 上述约束条件满足任意其一，气井即废弃。
[0088] 计算结果表明，I类气井均因达到极限生产时间20年而废弃；而20年对应的线性斜率计算结果为7.3943，据此，便可批量快速计算I类气井的最终累积产量(EUR)，计算结果如表5所示。
[0089] 表5
[0090]
[0091]
[0092] 基于同一发明构思，本发明实施例中还提供了一种确定气井最终累积产量的装置，如下面的实施例所述。由于确定气井最终累积产量的装置解决问题的原理与确定气井最终累积产量的方法相似，因此确定气井最终累积产量的装置的实施可以参见确定气井最终累积产量的方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的
软件和/或
硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。
[0093] 图8是本发明实施例的确定气井最终累积产量的装置的一种结构框图，如图8所示，包括：
[0094] 气井分类模块801，用于对目标区域内的气井进行分类，不同类的气井满足不同的开采效果标准；
[0095] 时间点确定模块802，用于确定所述目标区域内每口气井的产量递减的初始时间点；
[0096] 第一产量确定模块803，用于确定所述目标区域内每口气井在所述产量递减的初始时间点之前的累计产量，该累计产量为第一累计产量；
[0097] 第二产量确定模块804，用于确定所述目标区域内每口气井由所述产量递减的初始时间点至产量递减期内预设时间点的累计产量，该累计产量为第二累计产量；
[0098] 系数计算模块805，用于计算每类气井在产量递减期内的产量系数；
[0099] 最终累积产量确定模块806，用于针对每类气井，根据该类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，计算该类气井中每口气井的最终累积产量。
[0100] 在一个实施例中，所述气井分类模块，具体用于根据综合分类系数、生产测试时的测试产量、投产首年的平均日产量以及气井内部收益率，对目标区域内的气井进行分类。
[0101] 在一个实施例中，所述气井分类模块还用于通过以下公式计算所述综合分类系数：
[0102]
[0103] 其中，Ψ为气井的综合分类系数，无量纲；L为气井压裂水平段长度；N为改造段数；S为压裂加砂量；β为优质储层钻遇比例；脚标avg代表相应的样本平均值。
[0104] 在一个实施例中，所述系数计算模块，包括：
[0105] 产量获取单元，用于确定该类气井中每口气井由所述产量递减的初始时间点至产量递减期内任意时间点的累计产量，该累计产量为递减期内的阶段累积量，其中，相邻时间点之间相差预设时间步长；
[0106] 斜率获取单元，用于针对该类气井中所有气井，以第二累计产量为横坐标，以递减期内的阶段累积量为纵坐标，绘制散点图，对各散点进行线性拟合，得到与时间点对应的线性斜率；
[0107] 产量系数确定单元，用于绘制线性斜率与时间点的关系图，根据数据点分布趋势优选函数进行曲线拟合，得到线性斜率与时间点的函数表达式，该线性斜率为该类气井在产量递减期内的产量系数。
[0108] 在一个实施例中，线性斜率与时间点的函数表达式如下：
[0109] Gp(t)＝Gp0+k(t)·Gptn
[0110] 其中，Gp(t)为产量递减期内任意时间点的气井总累积产量；Gp0为第一累计产量；Gptn为第二累计产量；k(t)为线性斜率；t为产量递减期内的时间点。
[0111] 在一个实施例中，所述最终累积产量确定模块通过以下公式计算该类气井中每口气井的最终累积产量：
[0112] EUR＝Gp0+k(tabandon)·Gptn
[0113] 其中，EUR为每口气井的最终累积产量；Gp0为第一累计产量；Gptn为第二累计产量；k(tabandon)为该类气井在产量递减期内的产量系数；tabandon为气井最终寿命对应的时间点。
[0114] 在一个实施例中，所述最终累积产量确定模块还用于通过以下公式确定气井最终寿命对应的时间点：
[0115] 或tabandon≥tlimit
[0116] 其中，qabandon为气井经济极限产量；Gp(tabandon)为气井最终寿命对应的时间点的气井总累积产量；Gp(tabandon-1)为气井最终寿命对应时间点的前一时间点的气井总累积产量，相邻时间点之间的时间步长为一个月；tlimit为气井最终寿命。
[0117] 在另外一个实施例中，还提供了一种软件，该软件用于执行上述实施例及优选实施方式中描述的技术方案。
[0118] 在另外一个实施例中，还提供了一种存储介质，该存储介质中存储有上述软件，该存储介质包括但不限于：光盘、
软盘、
硬盘、可擦写存储器等。
[0119] 本发明实施例实现了如下技术效果：通过对目标区域内的气井进行分类，进而得到目标区域内每口气井的第一累计产量和第二累计产量，再确定出每类气井在产量递减期内的产量系数，最后，就可以根据每类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，批量计算出每类气井中每口气井的最终累积产量。该方案主要是确定出了每类气井在产量递减期内的产量系数，使得可以实现根据每类气井的产量系数、每口气井的第一累计产量以及每口气井的第二累计产量，批量计算出每类气井中每口气井的最终累积产量，相对现有技术而言避免了模型建立、历史拟合等环节，有利于在批量预测区块气井EUR的同时，减少工作量、减少耗时。
[0120] 显然，本领域的技术人员应该明白，上述的本发明实施例的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成
电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明实施例不限制于任何特定的硬件和软件结合。
[0121] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何
修改、等同替换、改进等，均应包含在本发明的保护范围之内。
