CN112832737A - 页岩气井eur的确定方法、装置、设备及存储介质 - Google Patents
页岩气井eur的确定方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN112832737A CN112832737A CN202110002077.2A CN202110002077A CN112832737A CN 112832737 A CN112832737 A CN 112832737A CN 202110002077 A CN202110002077 A CN 202110002077A CN 112832737 A CN112832737 A CN 112832737A
- Authority
- CN
- China
- Prior art keywords
- gas well
- production
- yield
- sample gas
- eur
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000003860 storage Methods 0.000 title claims abstract description 17
- 238000013507 mapping Methods 0.000 claims abstract description 89
- 238000004519 manufacturing process Methods 0.000 claims description 292
- 230000001186 cumulative effect Effects 0.000 claims description 82
- 230000007423 decrease Effects 0.000 claims description 63
- 230000015654 memory Effects 0.000 claims description 18
- 238000011156 evaluation Methods 0.000 abstract description 2
- 206010017076 Fracture Diseases 0.000 description 17
- 208000010392 Bone Fractures Diseases 0.000 description 13
- 230000006870 function Effects 0.000 description 13
- 230000009467 reduction Effects 0.000 description 12
- 230000008859 change Effects 0.000 description 11
- 230000002093 peripheral effect Effects 0.000 description 10
- 230000001133 acceleration Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 9
- 239000011148 porous material Substances 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 5
- 238000012795 verification Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000000919 ceramic Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 208000013201 Stress fracture Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- 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
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Sampling And Sample Adjustment (AREA)
Abstract
本申请提供了一种页岩气井EUR的确定方法、装置、设备及存储介质，属于油气储量评估领域。所述方法包括：通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型；通过每个样本气井的产量预测模型，预测每个样本气井在之后的目标时间段内的页岩气产量；将每个样本气井的总历史产量与目标时间段内的页岩气产量之和，确定为样本气井的EUR；基于多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系；对于任一未达到边界流形态的目标气井，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR。上述方法能够保证确定的页岩气井EUR的准确度。
Description
技术领域
本申请涉及油气储量评估领域，特别涉及一种页岩气井EUR的确定方法、装置、设备及存储介质。
背景技术
EUR(Estimated ultimate recovery，最终可采储量)是指从页岩气井投产开始，以某一指定日期、某一废弃压力界限或某一最低极限气量为截止条件所确定的可采出气量，其直接反映油气勘探项目的技术风险和经济价值，对油气储量的勘探和开发有决定性作用，因此，如何确定页岩气的EUR是目前研究的热点问题。
相关技术中，一般根据页岩气井的历史产量来预测之后的页岩气产量，进而得到该页岩气井的EUR，然而对于处于投产初期的页岩气井来讲，其产量变化趋势还不稳定，因此，获取的EUR不准确。
发明内容
本申请实施例提供了一种页岩气井EUR的确定方法、装置、设备及存储介质，能够保证确定的页岩气井EUR的准确度。所述技术方案如下：
一方面，提供了一种页岩气井EUR的确定方法，所述方法包括：
通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，所述第一历史产量数据表示所述样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量；
通过所述每个样本气井的产量预测模型，预测所述每个样本气井在之后的目标时间段内的页岩气产量，所述目标时间段包括多个单位时间段，所述样本气井在所述目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量；
将所述每个样本气井的总历史产量与所述目标时间段内的页岩气产量之和，确定为所述样本气井的EUR；
基于所述多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系，所述推定累计产量是根据所述样本气井在达到边界流形态之前的产量下降速率确定的；
对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR。
在一种可能的实现方式中，所述基于所述多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系之前，所述方法还包括：
对于所述每个样本气井，基于所述样本气井在达到边界流形态之前的产量下降速率，推定所述样本气井的产量下降为零的截止时间段；
基于所述产量下降速率，确定所述样本气井在所述截止时间段之前的推定累计产量。
在另一种可能的实现方式中，所述基于所述样本气井在达到边界流形态之前的产量下降速率，推定所述样本气井的产量下降为零的截止时间段，包括：
在坐标系中标注所述样本气井的第二历史产量数据，所述第二历史产量数据表示所述样本气井在达到边界流形态之前的多个第二历史时间段内的页岩气产量，所述坐标系的横坐标表示时间段，所述坐标系的纵坐标表示页岩气产量；
基于所述样本气井的产量达到峰值之后且所述样本气井达到边界流形态之前的多个第二历史时间段内的页岩气产量拟合出至少一条线段，延长所述至少一条线段中斜率最小的线段至所述坐标系的横轴；
将所述斜率最小的线段与所述横轴的交点的横坐标确定为所述截止时间段。
在另一种可能的实现方式中，所述通过所述每个样本气井的产量预测模型，预测所述每个样本气井在之后的目标时间段内的页岩气产量，包括：
对于所述每个样本气井，通过所述样本气井对应的产量预测模型预测所述样本气井在之后的多个单位时间段内的页岩气产量；
确定所述最低极限产量对应的第一单位时间段，将所述第一单位时间段以及所述第一单位时间段之前的多个单位时间段确定为所述目标时间段；
统计所述样本气井在所述目标时间段内的页岩气产量。
在另一种可能的实现方式中，所述基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR之前，所述方法还包括：
将验证气井的总历史产量与所述验证气井在所述目标时间段内的页岩气产量之和，确定为所述验证气井的第一EUR，所述验证气井的页岩气的流动形态达到边界流形态；
基于所述验证气井的推定累计产量和所述映射关系，确定所述验证气井的第二EUR；
将所述第一EUR和所述第二EUR的绝对差值与所述第一EUR的商，确定为所述映射关系的误差；
在所述误差小于参考误差的前提下，执行对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR的步骤。
在另一种可能的实现方式中，所述通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型之前，所述方法还包括：
从所述多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井。
在另一种可能的实现方式中，所述从所述多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井，包括：
在坐标系中标注生产时长不小于所述第一参考时长的每个页岩气井的历史产量数据，所述历史产量数据表示所述页岩气井在多个历史时间段内的页岩气产量，所述坐标系的横坐标表示时间段，所述坐标系的纵坐标表示页岩气产量；
确定所述页岩气井的产量趋势，若所述页岩气井的产量趋势包括至少两个产量下降速率不同的产量下降阶段，则确定产量下降速率最小的产量下降阶段的持续时长；
若所述持续时长不小于所述第二参考时长，则将所述页岩气井确定为所述样本气井。
另一方面，提供了一种页岩气井EUR的确定装置，所述装置包括：
模型建立模块，被配置为通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，所述第一历史产量数据表示所述样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量；
产量预测模块，被配置为通过所述每个样本气井的产量预测模型，预测所述每个样本气井在之后的目标时间段内的页岩气产量，所述目标时间段包括多个单位时间段，所述样本气井在所述目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量；
第一EUR确定模块，被配置为将所述每个样本气井的总历史产量与所述目标时间段内的页岩气产量之和，确定为所述样本气井的EUR；
关系建立模块，被配置为基于所述多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系，所述推定累计产量是根据所述样本气井在达到边界流形态之前的产量下降速率确定的；
第二EUR确定模块，被配置为对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR。
在一种可能的实现方式中，所述装置还包括：
时间推定模块，被配置为对于所述每个样本气井，基于所述样本气井在达到边界流形态之前的产量下降速率，推定所述样本气井的产量下降为零的截止时间段；
产量确定模块，被配置为基于所述产量下降速率，确定所述样本气井在所述截止时间段之前的推定累计产量。
在另一种可能的实现方式中，所述时间推定模块，被配置为在坐标系中标注所述样本气井的第二历史产量数据，所述第二历史产量数据表示所述样本气井在达到边界流形态之前的多个第二历史时间段内的页岩气产量，所述坐标系的横坐标表示时间段，所述坐标系的纵坐标表示页岩气产量；基于所述样本气井的产量达到峰值之后且所述样本气井达到边界流形态之前的多个第二历史时间段内的页岩气产量拟合出至少一条线段，延长所述至少一条线段中斜率最小的线段至所述坐标系的横轴；将所述斜率最小的线段与所述横轴的交点的横坐标确定为所述截止时间段。
在另一种可能的实现方式中，所述产量预测模块，被配置为对于所述每个样本气井，通过所述样本气井对应的产量预测模型预测所述样本气井在之后的多个单位时间段内的页岩气产量；确定所述最低极限产量对应的第一单位时间段，将所述第一单位时间段以及所述第一单位时间段之前的多个单位时间段确定为所述目标时间段；统计所述样本气井在所述目标时间段内的页岩气产量。
在另一种可能的实现方式中，所述装置还包括：
误差确定模块，被配置为将验证气井的总历史产量与所述验证气井在所述目标时间段内的页岩气产量之和，确定为所述验证气井的第一EUR，所述验证气井的页岩气的流动形态达到边界流形态；基于所述验证气井的推定累计产量和所述映射关系，确定所述验证气井的第二EUR；将所述第一EUR和所述第二EUR的绝对差值与所述第一EUR的商，确定为所述映射关系的误差；
所述第二EUR确定模块，被配置为在所述误差小于参考误差的前提下，执行对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR的步骤。
在另一种可能的实现方式中，所述装置还包括：
样本气井确定模块，被配置为从所述多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井。
在另一种可能的实现方式中，所述样本气井确定模块，被配置为在坐标系中标注生产时长不小于所述第一参考时长的每个页岩气井的历史产量数据，所述历史产量数据表示所述页岩气井在多个历史时间段内的页岩气产量，
所述坐标系的横坐标表示时间段，所述坐标系的纵坐标表示页岩气产量；确定所述页岩气井的产量趋势，若所述页岩气井的产量趋势包括至少两个产量下降速率不同的产量下降阶段，则确定产量下降速率最小的产量下降阶段的持续时长；若所述持续时长不小于所述第二参考时长，则将所述页岩气井确定为所述样本气井。
另一方面，提供了一种电子设备，所述电子设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述程序代码由所述处理器加载并执行以实现上述任一种可能实现方式中的页岩气井EUR的确定方法中执行的操作。
另一方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有至少一条程序代码，所述程序代码由处理器加载并执行以实现上述任一种可能实现方式中的页岩气井EUR的确定方法中执行的操作。
另一方面，提供了一种计算机程序产品，所述计算机程序产品中包括至少一条程序代码，所述程序代码由处理器加载并执行以实现上述任一种可能实现方式中的页岩气井EUR的确定方法中执行的操作。
本申请实施例提供的技术方案带来的有益效果至少包括：
本申请实施例提供的技术方案，考虑到页岩气流动形态达到边界流形态的页岩气井的产量变化趋势已经趋于稳定，且达到边界流形态之后的产量变化趋势更能体现之后的页岩气产量，因此，基于样本气井在达到边界流形态之后的历史时间段内的页岩气产量来预测之后的页岩气产量，保证了确定的样本气井EUR的准确度。在此基础上，考虑到处于投产初期的目标气井的推定累计产量能够准确确定出来，因此基于样本气井的推定累计产量和该样本气井的准确的EUR，建立推定累计产量与EUR之间的映射关系，使得能够根据目标气井的推定累计产量和该映射关系来确定目标气井的EUR，保证了确定的目标气井EUR的准确度。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种实施环境的示意图；
图2是本申请实施例提供的一种页岩气井EUR的确定方法的流程图；
图3是本申请实施例提供的一种页岩气井EUR的确定方法的流程图；
图4是本申请实施例提供的一种页岩气井的历史产量数据的示意图；
图5是本申请实施例提供的一种页岩气井的历史产量数据的示意图；
图6是本申请实施例提供的一种关系数据的示意图；
图7是本申请实施例提供的一种页岩气井EUR的确定方法的流程图；
图8是本申请实施例提供的一种页岩气井EUR的确定装置的框图；
图9是本申请实施例提供的一种终端的结构示意图；
图10是本申请实施例提供的一种服务器的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施方式作进一步地详细描述。
本申请所使用的术语“第一”、“第二”、“第三”、“第四”等可在本文中用于描述各种概念，但除非特别说明，这些概念不受这些术语限制。这些术语仅用于将一个概念与另一个概念区分。举例来说，在不脱离本申请的范围的情况下，可以将第一历史产量数据称为历史产量数据，且类似地，可将第二历史产量数据称为第一历史产量数据。
本申请所使用的术语“至少一个”、“多个”、“每个”、“任一”，至少一个包括一个、两个或两个以上，多个包括两个或两个以上，而每个是指对应的多个中的每一个，任一是指多个中的任意一个。举例来说，多个样本气井包括3个样本气井，而每个是指这3个样本气井中的每一个样本气井，任一是指这3个样本气井中的任意一个，可以是第一个，可以是第二个、也可以是第三个。
图1是本公开实施例提供的一种实施环境的示意图。参见图1，该实施环境包括终端101和服务器102。终端101和服务器102之间通过无线或者有线网络连接。可选地，终端101为电脑、手机、平板电脑或者其他终端。可选地，服务器102为该目标应用的后台服务器或者为提供云计算以及云存储等服务的云服务器。
可选地，终端101上安装由服务器102提供服务的目标应用，终端101能够通过该目标应用实现例如数据传输、消息交互等功能。可选地，目标应用为终端101操作系统中的目标应用，或者为第三方提供的目标应用。该目标应用具有确定页岩气井EUR的功能，当然，该目标应用还能够具有其他功能，本申请对此不做限制。
本申请实施例中，服务器102用于确定页岩气井的EUR，将EUR发送给终端101，而终端101用于展示该页岩气井的EUR。
本申请提供的页岩气井EUR的确定方法能够应用于确定任何页岩气井的EUR，例如，能够应用于确定处于投产初期的页岩气井的EUR，也能够应用于确定投产时间较长的页岩气井的EUR。
图2是本申请实施例提供的一种页岩气井EUR的确定方法的流程图。该方法实施例的执行主体为电子设备。参见图2，该方法包括：
201：通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，第一历史产量数据表示样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量。
202：通过每个样本气井的产量预测模型，预测每个样本气井在之后的目标时间段内的页岩气产量，目标时间段包括多个单位时间段，样本气井在目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量。
203：将每个样本气井的总历史产量与目标时间段内的页岩气产量之和，确定为样本气井的EUR。
204：基于多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系，推定累计产量是根据样本气井在达到边界流形态之前的产量下降速率确定的。
205：对于任一未达到边界流形态的目标气井，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR。
本申请实施例提供的技术方案，考虑到页岩气流动形态达到边界流形态的页岩气井的产量变化趋势已经趋于稳定，且达到边界流形态之后的产量变化趋势更能体现之后的页岩气产量，因此，基于样本气井在达到边界流形态之后的历史时间段内的页岩气产量来预测之后的页岩气产量，保证了确定的样本气井EUR的准确度。在此基础上，考虑到处于投产初期的目标气井的推定累计产量能够准确确定出来，因此基于样本气井的推定累计产量和该样本气井的准确的EUR，建立推定累计产量与EUR之间的映射关系，使得能够根据目标气井的推定累计产量和该映射关系来确定目标气井的EUR，保证了确定的目标气井EUR的准确度。
在一种可能的实现方式中，基于多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系之前，方法还包括：
对于每个样本气井，基于样本气井在达到边界流形态之前的产量下降速率，推定样本气井的产量下降为零的截止时间段；
基于产量下降速率，确定样本气井在截止时间段之前的推定累计产量。
在另一种可能的实现方式中，基于样本气井在达到边界流形态之前的产量下降速率，推定样本气井的产量下降为零的截止时间段，包括：
在坐标系中标注样本气井的第二历史产量数据，第二历史产量数据表示样本气井在达到边界流形态之前的多个第二历史时间段内的页岩气产量，坐标系的横坐标表示时间段，坐标系的纵坐标表示页岩气产量；
基于样本气井的产量达到峰值之后且样本气井达到边界流形态之前的多个第二历史时间段内的页岩气产量拟合出至少一条线段，延长至少一条线段中斜率最小的线段至坐标系的横轴；
将斜率最小的线段与横轴的交点的横坐标确定为截止时间段。
在另一种可能的实现方式中，通过每个样本气井的产量预测模型，预测每个样本气井在之后的目标时间段内的页岩气产量，包括：
对于每个样本气井，通过样本气井对应的产量预测模型预测样本气井在之后的多个单位时间段内的页岩气产量；
确定最低极限产量对应的第一单位时间段，将第一单位时间段以及第一单位时间段之前的多个单位时间段确定为目标时间段；
统计样本气井在目标时间段内的页岩气产量。
在另一种可能的实现方式中，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR之前，方法还包括：
将验证气井的总历史产量与验证气井在目标时间段内的页岩气产量之和，确定为验证气井的第一EUR，所述验证气井的页岩气的流动形态达到边界流形态；
基于验证气井的推定累计产量和映射关系，确定验证气井的第二EUR；
将第一EUR和第二EUR的绝对差值与第一EUR的商，确定为映射关系的误差；
在误差小于参考误差的前提下，执行对于任一未达到边界流形态的目标气井，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR的步骤。
在另一种可能的实现方式中，通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型之前，方法还包括：
从多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井。
在另一种可能的实现方式中，从多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井，包括：
在坐标系中标注生产时长不小于第一参考时长的每个页岩气井的历史产量数据，历史产量数据表示页岩气井在多个历史时间段内的页岩气产量，坐标系的横坐标表示时间段，坐标系的纵坐标表示页岩气产量；
确定页岩气井的产量趋势，若页岩气井的产量趋势包括至少两个产量下降速率不同的产量下降阶段，则确定产量下降速率最小的产量下降阶段的持续时长；
若持续时长不小于第二参考时长，则将页岩气井确定为样本气井。
上述所有可选技术方案，可以采用任意结合形成本申请的可选实施例，在此不再一一赘述。
图3是本申请实施例提供的一种页岩气井EUR的确定方法的流程图。参见图3，该方法包括：
301：服务器从多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井。
页岩气井在投产后，随着页岩气储层中的页岩气不断产出，页岩气井井底的页岩气流动形态会产生阶段性的变化，由一开始的线性流形态变为边界流形态。在页岩气井处于线性流形态的情况下，页岩气井所产出的页岩气主要来源于水力压裂缝，此时，页岩气井的产量下降速率较快。而在页岩气井处于边界流形态的情况下，页岩气井所产出的页岩气则主要来源于水力压裂缝周围的页岩基质孔隙和天然微裂缝孔隙，此时，页岩气井的产量下降速率较慢。
生产时长是指页岩气井从投产时刻开始到当前时刻所持续的总时长。处于边界流形态的持续时长是指：从页岩气井进入到边界流形态的开始时刻到当前时刻的持续时长。
在一种可能的实现方式中，服务器从多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井，包括下述步骤(1)-(3)：
(1)服务器在坐标系中标注生产时长不小于第一参考时长的每个页岩气井的历史产量数据。
其中，历史产量数据表示页岩气井在多个历史时间段内的页岩气产量。该多个历史时间段的时长相同。可选地，该历史时间段的时长为一天，当然，该历史时间段的时长为还能够为其他时长，例如一个月，本申请实施例对此不做限制。其中，坐标系的横坐标表示时间段，坐标系的纵坐标表示页岩气产量。图4为在坐标系中标注出一个页岩气井的历史产量数据所得到的示意图。参考图4，每个点的横坐标表示时间段，该时间段的时长为一天。每个点的纵坐标表示日产气量，即在一天内的页岩气产量。
其中，第一参考时长为页岩气井是否达到边界流形态的生产时长分界线，若页岩气井的生产时长不小于第一参考时长，则认为该页岩气井已经处于边界流形态。若页岩气井的生产时长小于第一参考时长，则认为该页岩气井还未达到边界流形态。
第一参考时长能够根据页岩气井所在的页岩气田的实际情况设置，例如，将第一参考时长设置为1年，当然，第一参考时长还能够设置为其他数值，本申请实施例对此不做限制。
(2)服务器确定页岩气井的产量趋势，若页岩气井的产量趋势包括至少两个产量下降速率不同的产量下降阶段，则确定产量下降速率最小的产量下降阶段的持续时长。
经过步骤(1)，服务器在坐标系中标注出了已经处于边界流形态的页岩气井的历史产量数据。由于页岩气井在达到边界流形态之前的产量下降速率大于处于边界流形态的产量下降速率，且页岩气储层中水力压裂缝的变化可能导致页岩气井的产量下降速率发生变化，因此，页岩气井的产量趋势会呈现至少两个产量下降速率不同的产量下降阶段。并且，由于页岩气储层中水力压裂缝的变化往往在页岩气井达到边界流形态之前，因此，在至少两个产量下降阶段中，产量下降速率最小的产量下降阶段对应页岩气井的边界流形态，而其他产量下降阶段对应页岩气井的线性流形态，即未达到边界流形态的页岩气流动形态。因此，产量下降速率最小的产量下降阶段的持续时长即表示页岩气井处于边界流形态的持续时长。
(3)若该持续时长不小于第二参考时长，服务器则将页岩气井确定为样本气井。
其中，第二参考时长能够根据实际情况设置，例如，第二参考时长能够设置为1个月，本申请实施例对此不做限制。
在本申请实施例中，通过在坐标系中标注历史产量数据并确定产量趋势的方法，使得能够根据页岩气井的产量下降规律直观地确定出页岩气井处于边界流形态的持续时长，进而快速地确定出样本气井。
可选地，服务器根据流动形态诊断原理对多个生产时长不小于第一参考时长的页岩气井进行流动形态诊断，确定每个页岩气井是否达到边界流形态，以及处于边界流形态的持续时长，然后选择出处于边界流形态的持续时长不小于第二参考时长的多个样本气井。或者服务器根据试井解释成果资料确定每个页岩气井是否达到边界流形态，以及处于边界流形态的持续时长，然后选择出处于边界流形态的持续时长不小于第二参考时长的多个样本气井。其中，试井解释成果资料中包括页岩气井的产量、压力、储层结构、储层性质等多项数据，服务器能够根据试井解释成果资料确定每个页岩气井是否达到边界流形态，以及处于边界流形态的持续时长。由于页岩气井的历史产量数据可能发生波动，通过产量趋势可能无法准确判断页岩气井达到边界流形态的时间，而通过流动形态诊断和试井解释成果资料则能够准确确定页岩气井达到边界流形态的时间，进而能够准确确定出样本气井。
需要说明的一点是，上述多种确定样本气井的方式能够以任意方式结合，本申请实施例对此不做限制。
302：服务器通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，第一历史产量数据表示样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量。
其中，多个第一历史时间段的时长相同。可选地，第一历史时间段的时长为一天，当然，第一历史时间段的时长为还能够为其他时长，例如一个月，本申请实施例对此不做限制。
可选地，服务器通过双曲递减模型来拟合多个样本气井的第一历史产量数据，得到产量预测模型，当然，服务器还能够通过其他模型来拟合第一历史产量数据，例如扩展指数模型等，本申请实施例对此不做限制。
303：服务器通过每个样本气井的产量预测模型，预测每个样本气井在之后的目标时间段内的页岩气产量。
其中，目标时间段包括多个单位时间段，并且样本气井在目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量。
其中，多个单位时间段的时长相同。可选地，单位时间段的时长为一天，当然，单位时间段的时长为还能够为其他时长，例如一个月，本申请实施例对此不做限制。样本气井在多个单位时间段内的页岩气产量是逐渐下降的，并且，在最后一个单位时间段内的页岩气产量达到最低极限产量。
最低极限产量为技术极限产量或经济极限产量。其中，技术极限产量是指开采技术上能够支持的最低产量，页岩气的产量如果低于技术极限产量，在当前的开采技术下，剩余的页岩气无法开采出来。可选地，技术极限产量能够根据废弃的页岩气井的最低产量确定。经济极限产量是指能够支持生产成本的最低产量。经济极限产量能够利用油气价格，操作成本等经济参数按照行业标准确定。
在一种可能的实现方式中，服务器通过每个样本气井的产量预测模型，预测每个样本气井在之后的目标时间段内的页岩气产量，包括：对于每个样本气井，服务器通过样本气井对应的产量预测模型预测样本气井在之后的多个单位时间段内的页岩气产量；确定最低极限产量对应的第一单位时间段，将第一单位时间段以及第一单位时间段之前的多个单位时间段确定为目标时间段；统计样本气井在目标时间段内的页岩气产量。
例如，服务器通过样本气井对应的产量预测模型预测出样本气井在之后的100天内每天的页岩气产量，其中，第80天内的页岩气产量为最低极限产量，则第80天即为第一单位时间段，目标时间段即为这100天内的前80天。然后，将这80天的页岩气产量相加，则统计出样本气井在目标时间段内的页岩气产量。
在本申请实施例中，考虑到样本气井在达到边界流形态之后，其产量变化趋势已经趋于稳定，因此，利用样本气井在达到边界流形态之后的历史产量数据来建立产量预测模型，则该产量预测模型的准确度高，预测的之后的页岩气产量更加准确，从而能够保证确定的样本气井的EUR的准确度。
304：服务器将每个样本气井的总历史产量与目标时间段内的页岩气产量之和，确定为该样本气井的EUR。
其中，样本气井的总历史产量包括样本气井在达到边界流形态之前的历史产量和达到边界流形态之后的历史产量。由于样本气井在之前的多个历史时间段内的页岩气产量能够从页岩气井的生产数据中直接确定，因此，只需要将该多个历史时间段内的页岩气产量相加，则能够得到样本气井的总历史产量。然后，再将该总历史产量与预测的样本气井在之后的目标时间段内的页岩气产量相加，则得到该样本气井的EUR。
305：服务器根据样本气井在达到边界流形态之前的产量下降速率确定推定累计产量。
页岩气井在投产后，在未达到边界流形态之前，页岩气井所产出的页岩气主要来源于水力压裂缝，此时，页岩气井的产量下降速率较快。而在页岩气井达到边界流形态之后，由于水力压裂缝周围的页岩基质孔隙和天然裂缝孔隙气体的补给，使得页岩气的产量下降速率变慢。而服务器根据样本气井在达到边界流形态之前的产量下降速率确定的推定累计产量则表示：样本气井的水力压裂缝中总共的页岩气产量，该页岩气产量不包括水力压裂缝周围的页岩基质孔隙和天然裂缝孔隙中的页岩气产量。
在一种可能的实现方式中，服务器根据样本气井在达到边界流形态之前的产量下降速率确定推定累计产量，包括：对于每个样本气井，服务器基于样本气井在达到边界流形态之前的产量下降速率，推定样本气井的产量下降为零的截止时间段；基于产量下降速率，确定样本气井在截止时间段之前的推定累计产量。
在本申请实施例中，考虑到不存在周围补给区域页岩基质孔隙和天然裂缝孔隙的气体补给且开采技术允许的情况下，样本气井的页岩气产量会按照该产量下降速率持续下降，直至产量为零。因此，基于该产量下降速率推定样本气井的产量下降为零的截止时间段，并且基于该产量下降速率，确定样本气井在截止时间段之前的推定累计产量，则该推定累计产量能够表示样本气井的水力压裂缝中总共的页岩气产量。另外，由于样本气井的水力压裂缝的体积大小能够决定周围补给区域的体积大小，因此，样本气井的水力压裂缝中的页岩气产量决定了补给区域中的页岩气产量，而样本气井的EUR又由样本气井的水力压裂缝和补给区域中的页岩气产量决定，因此实际上，样本气井的水力压裂缝中的页岩气产量决定了样本气井的EUR，也即是，样本气井的水力压裂缝中的页岩气产量与EUR之间存在确定的关系。因此，在本申请实施例中，确定出能够表示水力压裂缝中的页岩气产量的推定累计产量后，就能够确定出该推定累计产量与EUR之间的映射关系。
在一种可能的实现方式中，服务器基于样本气井在达到边界流形态之前的产量下降速率，推定样本气井的产量下降为零的截止时间段，包括：服务器在坐标系中标注样本气井的第二历史产量数据，第二历史产量数据表示样本气井在达到边界流形态之前的多个第二历史时间段内的页岩气产量；服务器基于样本气井的产量达到峰值之后且样本气井达到边界流形态之前的多个第二历史时间段内的页岩气产量拟合出至少一条线段，延长至少一条线段中斜率最小的线段至坐标系的横轴；服务器将斜率最小的线段与横轴的交点的横坐标确定为截止时间段。
其中，多个第二历史时间段的时长相同。可选地，第二历史时间段的时长为一天，当然，第二历史时间段的时长为还能够为其他时长，例如一个月，本申请实施例对此不做限制。其中，坐标系的横坐标表示时间段，坐标系的纵坐标表示页岩气产量。
继续参考图4，其中分隔线左侧的多个点表示在达到边界流形态之前的多个第二历史时间段内的页岩气产量，q0表示样本气井的产量峰值，q0至t之间的线段即为斜率最小的线段，t表示截止时间段。在此情况下，样本气井在截止时间段之前的推定累计产量能够通过下述公式(1)和(2)确定。
S＝S0+S1 (1)
其中，S表示推定累计产量，S0表示样本气井达到峰值时间段t0前的累计产量，S1表示从峰值时间段t0到截止时间段t的累计产量，q0表示样本气井的产量峰值。
图4对应样本气井在达到边界流形态之前只有1个产量下降阶段的情况，图5对应样本气井在达到边界流形态之前有2个产量下降阶段的情况。参考图5，q0至q1之间的线段以及q1至t2之间的线段均为拟合出的线段，且q1至t2之间的线段为斜率最小的线段，t2表示截止时间段，在此情况下，样本气井在截止时间段之前的推定累计产量能够通过下述公式(3)、(4)、(5)来确定。
S＝S0+S11+S12 (3)
其中，S表示推定累计产量，S0表示样本气井达到峰值时间段t0前的累计产量，S11表示峰值时间段t0到时间段t1的累计产量，S12表示时间段t1到截止时间段t2的累计产量，q0表示样本气井的产量峰值，q1表示产量下降速率发生变化的时间段t1对应的页岩气产量。
需要说明的一点是，对于样本气井在达到边界流形态之前包括2个以上产量下降阶段的情况，确定推定累计产量的方式与确定包括2个产量下降阶段的样本气井的推定累计产量的方式同理，此处不再赘述。当然，还能够通过其他方式确定在达到边界流形态之前，包括2个以上产量下降阶段的样本气井的推定累计产量，例如，按照产量趋势将2个以上产量下降阶段等效为2个产量下降阶段来处理。
306：服务器基于多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系。
可选地，服务器将一个样本气井的推定累计产量和EUR作为一个关系数据，通过参考映射关系拟合多个关系数据，得到推定累计产量与EUR之间的映射关系。
可选地，参考映射关系为线性关系。参考映射关系能够根据实际情况设置，例如，参考映射关系设置为下述公式(6)。
EUR＝AS+B (6)
其中，A和B表示关系参数，S表示推定累计产量。
可选地，服务器通过参考映射关系拟合多个关系数据的实现方式为：服务器通过参考映射关系和最小二乘算法来拟合多个关系数据。当然，服务器也能够通过最小二乘算法以外的其他数据拟合算法来拟合多个关系数据，本申请实施例对此不做限制。图6为在坐标系中标注出多个关系数据的示意图。参考图6，其中，每个点表示一个关系数据，该点的横坐标表示推定累计产量，纵坐标表示EUR。示例性的，该EUR为TEUR。图6中的线段为拟合该多个关系数据所得到的，该线段对应的映射关系为y＝0.0002x+0.4727，其中，x表示推定累计产量，y表示EUR。
在本申请实施例中，不是根据一个样本气井的关系数据来确定推定累计产量与EUR之间的映射关系，而是根据多个样本气井的关系数据来确定推定累计产量与EUR之间的映射关系，如此能够保证确定的映射关系客观准确。
可选地，在确定出推定累计产量与EUR之间的映射关系前，服务器先根据多个关系数据，确定推定累计产量与EUR之间的相关度，在相关度不小于参考阈值的情况下，基于该多个关系数据确定推定累计产量与EUR之间的相关度。在确定的相关度小于参考阈值的情况下，从多个关系数据中筛选多个关系数据，根据筛选出的多个关系数据确定推定累计产量与EUR之间的相关度，直至确定的相关度不小于参考阈值，然后服务器根据筛选出的多个关系数据确定推定累计产量与EUR之间的映射关系。如此，能够排除部分异常的关系数据对映射关系的影响，使得确定的映射关系更加准确。
可选地，该相关度为皮尔逊相关系数的平方，当然，该相关度也能够为其他相关度，本申请实施例对此不做限制。可选地，相关度的参考阈值能够根据实际情况设置，例如，将参考阈值设置为0.7，本申请实施例对此不做限制。继续参考图6，其中，R2表示根据其中的多个关系数据确定的推定累计产量与EUR之间的相关度。
307：对于任一未达到边界流形态的目标气井，服务器基于目标气井的推定累计产量和映射关系，确定目标气井的EUR。
服务器在基于目标气井的推定累计产量和映射关系，确定目标气井的EUR之前，先要确定目标气井的推定累计产量。服务器确定目标气井的推定累计产量的实现方式与确定样本气井的推定累计产量的实现方式同理，此处不再赘述。
确定出目标气井的推定累计产量后，将该推定累计产量代入该映射关系，则能够得到目标气井的EUR。
需要说明的一点是，在投产初期，页岩气井的页岩气流动形态还未达到边界流形态，但页岩气井的产量下降速率能够准确确定，因此，能够根据页岩气井的产量下降速率准确确定出页岩气井的推定累计产量，进而能够根据该推定累计产量和映射关系准确确定该页岩气井的EUR。
需要说明的另一点是，在本申请实施例中，EUR包括TEUR(Technique estimatedultimate recovery，技术能够实现的最终可采储量)和EEUR(Economic estimatedultimate recovery，具有经济效益的最终可采储量)。上述在确定样本气井在目标时间段内的页岩气产量时，若样本气井在目标时间段内的最后一个单位时间段内的页岩气产量达到技术极限产量，则映射关系表示推定累计产量与TEUR之间的映射关系，若样本气井在目标时间段内的最后一个单位时间段内的页岩气产量达到经济极限产量，则映射关系表示推定累计产量与EEUR之间的映射关系。相应的，在确定目标气井的EUR时，若采用的映射关系为推定累计产量与TEUR之间的映射关系，则确定的EUR为TEUR，若采用的映射关系为推定累计产量与EEUR之间的映射关系，则确定的EUR为EEUR。
在一种可能的实现方式中，服务器基于目标气井的推定累计产量和映射关系，确定目标气井的EUR之前，还要确定映射关系的误差，在误差小于参考误差的前提下，才执行对于任一未达到边界流形态的目标气井，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR的步骤。如此，能够进一步保证确定的目标气井的EUR的准确度。其中，参考误差能够根据实际情况设置，本申请实施例对此不做限制。
在一种可能的实现方式中，服务器确定映射关系的误差，包括：服务器将验证气井的总历史产量与验证气井在目标时间段内的页岩气产量之和，确定为验证气井的第一EUR；基于验证气井的推定累计产量和映射关系，确定验证气井的第二EUR；将第一EUR和第二EUR的绝对差值与第一EUR的商，确定为映射关系的误差。其中，验证气井的页岩气的流动形态达到边界流形态。
可选地，验证气井的数量为多个，相应的，对于每个验证气井，确定映射关系的误差，将多个映射关系的误差的平均值作为最终的误差。如此，极大地提高了确定的误差的准确度。
其中，验证气井的确定方法与样本气井的确定方法同理。验证气井的总历史产量与样本气井的总历史产量的确定方法同理，验证气井在目标时间段内的页岩气产量与样本气井在目标时间段内的页岩气产量的确定方法同理，此处不再赘述。
由于验证气井的页岩气流动形态已经达到边界流形态，其页岩气产量趋于稳定，因此，服务器所预测的验证气井在目标时间段内的页岩气产量准确度高，从而获取的第一EUR的准确度高，能够表示该验证气井的真实的EUR。因此，第一EUR和第二EUR之间的绝对差值与第一EUR的商能够客观地表示映射关系的误差。
需要说明的再一点是，上述方法仅以服务器为执行主体为例进行说明，实际上，该方法的执行主体能够为除了服务器之外的其他电子设备，本申请对此不做限制。
以下为该页岩气井EUR的确定方法的具体应用：
某页岩气田共有99个页岩气井，通过上述方法确定其中70个页岩气井达到了边界流形态，其中55个页岩气井的产量趋势与图4类似，另外15个页岩气井的产量趋势与图5类似。通过双曲递减模型拟合这70个页岩气井在达到边界流形态之后的历史产量数据，得到每个页岩气井的产量预测模型后，按照最低极限产量为1000m3/day预测了每个页岩气井在之后的目标时间段内的页岩气产量，并确定出了这70个页岩气井的TEUR。选择其中前60个页岩气井作为样本气井，建立了推定累计产量与EUR之间的映射关系，剩余10个页岩气井作为验证气井，来验证映射关系的准确度。其中，建立的映射关系为TEUR＝0.0002S+0.4727。
下述表1为这70个页岩气井的相关数据。其中，S0表示样本气井达到峰值时间段t0前的累计产量，S1表示从峰值时间段t0到截止时间段t的累计产量，S11表示峰值时间段t0到时间段t1的累计产量，S12表示时间段t1到截止时间段t2的累计产量，S表示推定累计产量。
表1
在通过10个验证气井验证映射关系的误差后，确定出每个验证气井对应的误差都在15％以内，且平均误差为9.764％。下述表2为验证气井的误差数据。
表2
图7为确定EUR的流程图。参考图7，首先获取目标页岩气田中的多个已投产的页岩气井的历史产量数据。然后对多个历史产量数据进行预处理，以去除历史产量数据中的异常数据，其中，异常数据是指错误的产量数据。根据多个历史产量数据对应的产量-时间曲线筛选多个样本气井。接着确定每个样本气井的推定累计产量和EUR。然后通过线性拟合的方式，对多个样本气井的推定累计产量和EUR进行拟合，得到推定累计产量与EUR之间的映射关系。然后通过验证气井确定映射关系的误差，在误差小于参考误差的前提下，根据该映射关系，确定目标气井的EUR。
需要说明的一点是，本申请实际上是将生产时间较长且EUR计算结果可靠的页岩气井作为样本气井，建立样本气井的EUR与投产初期的推定累计产量的映射关系，以此映射关系来计算生产时间短，生产规律尚不明确的页岩气井的EUR。
本申请实施例提供的技术方案，考虑到页岩气流动形态达到边界流形态的页岩气井的产量变化趋势已经趋于稳定，且达到边界流形态之后的产量变化趋势更能体现之后的页岩气产量，因此，基于样本气井在达到边界流形态之后的历史时间段内的页岩气产量来预测之后的页岩气产量，保证了确定的样本气井EUR的准确度。在此基础上，考虑到处于投产初期的目标气井的推定累计产量能够准确确定出来，因此基于样本气井的推定累计产量和该样本气井的准确的EUR，建立推定累计产量与EUR之间的映射关系，使得能够根据目标气井的推定累计产量和该映射关系来确定目标气井的EUR，保证了确定的目标气井EUR的准确度。
图8是本申请实施例提供的一种页岩气井EUR的确定装置的框图。参见图8，该装置包括：
模型建立模块801，被配置为通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，第一历史产量数据表示样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量；
产量预测模块802，被配置为通过每个样本气井的产量预测模型，预测每个样本气井在之后的目标时间段内的页岩气产量，目标时间段包括多个单位时间段，样本气井在目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量；
第一EUR确定模块803，被配置为将每个样本气井的总历史产量与目标时间段内的页岩气产量之和，确定为样本气井的EUR；
关系建立模块804，被配置为基于多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系，推定累计产量是根据样本气井在达到边界流形态之前的产量下降速率确定的；
第二EUR确定模块805，被配置为对于任一未达到边界流形态的目标气井，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR。
在一种可能的实现方式中，装置还包括：
时间推定模块，被配置为对于每个样本气井，基于样本气井在达到边界流形态之前的产量下降速率，推定样本气井的产量下降为零的截止时间段；
产量确定模块，被配置为基于产量下降速率，确定样本气井在截止时间段之前的推定累计产量。
在另一种可能的实现方式中，时间推定模块，被配置为在坐标系中标注样本气井的第二历史产量数据，第二历史产量数据表示样本气井在达到边界流形态之前的多个第二历史时间段内的页岩气产量，坐标系的横坐标表示时间段，坐标系的纵坐标表示页岩气产量；基于样本气井的产量达到峰值之后且样本气井达到边界流形态之前的多个第二历史时间段内的页岩气产量拟合出至少一条线段，延长至少一条线段中斜率最小的线段至坐标系的横轴；将斜率最小的线段与横轴的交点的横坐标确定为截止时间段。
在另一种可能的实现方式中，产量预测模块802，被配置为对于每个样本气井，通过样本气井对应的产量预测模型预测样本气井在之后的多个单位时间段内的页岩气产量；确定最低极限产量对应的第一单位时间段，将第一单位时间段以及第一单位时间段之前的多个单位时间段确定为目标时间段；统计样本气井在目标时间段内的页岩气产量。
在另一种可能的实现方式中，装置还包括：
误差确定模块，被配置为将验证气井的总历史产量与验证气井在目标时间段内的页岩气产量之和，确定为验证气井的第一EUR，验证气井的页岩气的流动形态达到边界流形态；基于验证气井的推定累计产量和映射关系，确定验证气井的第二EUR；将第一EUR和第二EUR的绝对差值与第一EUR的商，确定为映射关系的误差；
第二EUR确定模块805，被配置为在误差小于参考误差的前提下，执行对于任一未达到边界流形态的目标气井，基于目标气井的推定累计产量和映射关系，确定目标气井的EUR的步骤。
在另一种可能的实现方式中，装置还包括：
样本气井确定模块，被配置为从多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井。
在另一种可能的实现方式中，样本气井确定模块，被配置为在坐标系中标注生产时长不小于第一参考时长的每个页岩气井的历史产量数据，历史产量数据表示页岩气井在多个历史时间段内的页岩气产量，
坐标系的横坐标表示时间段，坐标系的纵坐标表示页岩气产量；确定页岩气井的产量趋势，若页岩气井的产量趋势包括至少两个产量下降速率不同的产量下降阶段，则确定产量下降速率最小的产量下降阶段的持续时长；若持续时长不小于第二参考时长，则将页岩气井确定为样本气井。
本申请实施例提供的技术方案，考虑到页岩气流动形态达到边界流形态的页岩气井的产量变化趋势已经趋于稳定，且达到边界流形态之后的产量变化趋势更能体现之后的页岩气产量，因此，基于样本气井在达到边界流形态之后的历史时间段内的页岩气产量来预测之后的页岩气产量，保证了确定的样本气井EUR的准确度。在此基础上，考虑到处于投产初期的目标气井的推定累计产量能够准确确定出来，因此基于样本气井的推定累计产量和该样本气井的准确的EUR，建立推定累计产量与EUR之间的映射关系，使得能够根据目标气井的推定累计产量和该映射关系来确定目标气井的EUR，保证了确定的目标气井EUR的准确度。
需要说明的是：上述实施例提供的页岩气井EUR的确定装置在确定页岩气井的EUR时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将电子设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的页岩气井EUR的确定装置与页岩气井EUR的确定方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本申请实施例还提供了一种电子设备，该电子设备包括处理器和存储器，存储器中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以实现上述实施例的页岩气井EUR的确定方法中执行的操作。
可选地，该电子设备提供为终端。图9示出了本申请一个示例性实施例提供的终端900的结构框图。该终端900可以是：智能手机、平板电脑、MP3播放器(Moving PictureExperts Group Audio Layer III，动态影像专家压缩标准音频层面3)、MP4(MovingPicture Experts Group Audio Layer IV，动态影像专家压缩标准音频层面4)播放器、笔记本电脑或台式电脑。终端900还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。
终端900包括有：处理器901和存储器902。
处理器901可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器901可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器901也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器901可以集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器901还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。
存储器902可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器902还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器902中的非暂态的计算机可读存储介质用于存储至少一个程序代码，该至少一个程序代码用于被处理器901所执行以实现本申请中方法实施例提供的页岩气井EUR的确定方法。
在一些实施例中，终端900还可选包括有：外围设备接口903和至少一个外围设备。处理器901、存储器902和外围设备接口903之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口903相连。具体地，外围设备包括：射频电路904、显示屏905、摄像头组件906、音频电路907、定位组件908和电源90中的至少一种。
外围设备接口903可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器901和存储器902。在一些实施例中，处理器901、存储器902和外围设备接口903被集成在同一芯片或电路板上；在一些其他实施例中，处理器901、存储器902和外围设备接口903中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
射频电路904用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路904通过电磁信号与通信网络以及其他通信设备进行通信。射频电路904将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路904包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路904可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：城域网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路904还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。
显示屏905用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏905是触摸显示屏时，显示屏905还具有采集在显示屏905的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器901进行处理。此时，显示屏905还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏905可以为一个，设置终端900的前面板；在另一些实施例中，显示屏905可以为至少两个，分别设置在终端900的不同表面或呈折叠设计；在另一些实施例中，显示屏905可以是柔性显示屏，设置在终端900的弯曲表面上或折叠面上。甚至，显示屏905还可以设置成非矩形的不规则图形，也即异形屏。显示屏905可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode，有机发光二极管)等材质制备。
摄像头组件906用于采集图像或视频。可选地，摄像头组件906包括前置摄像头和后置摄像头。前置摄像头设置在终端的前面板，后置摄像头设置在终端的背面。在一些实施例中，后置摄像头为至少两个，分别为主摄像头、景深摄像头、广角摄像头、长焦摄像头中的任意一种，以实现主摄像头和景深摄像头融合实现背景虚化功能、主摄像头和广角摄像头融合实现全景拍摄以及VR(Virtual Reality，虚拟现实)拍摄功能或者其它融合拍摄功能。在一些实施例中，摄像头组件906还可以包括闪光灯。闪光灯可以是单色温闪光灯，也可以是双色温闪光灯。双色温闪光灯是指暖光闪光灯和冷光闪光灯的组合，可以用于不同色温下的光线补偿。
音频电路907可以包括麦克风和扬声器。麦克风用于采集用户及环境的声波，并将声波转换为电信号输入至处理器901进行处理，或者输入至射频电路904以实现语音通信。出于立体声采集或降噪的目的，麦克风可以为多个，分别设置在终端900的不同部位。麦克风还可以是阵列麦克风或全向采集型麦克风。扬声器则用于将来自处理器901或射频电路904的电信号转换为声波。扬声器可以是传统的薄膜扬声器，也可以是压电陶瓷扬声器。当扬声器是压电陶瓷扬声器时，不仅可以将电信号转换为人类可听见的声波，也可以将电信号转换为人类听不见的声波以进行测距等用途。在一些实施例中，音频电路907还可以包括耳机插孔。
定位组件908用于定位终端900的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件908可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统、俄罗斯的格雷纳斯系统或欧盟的伽利略系统的定位组件。
电源90用于为终端900中的各个组件进行供电。电源90可以是交流电、直流电、一次性电池或可充电电池。当电源90包括可充电电池时，该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。
在一些实施例中，终端900还包括有一个或多个传感器910。该一个或多个传感器910包括但不限于：加速度传感器911、陀螺仪传感器912、压力传感器913、指纹传感器914、光学传感器915以及接近传感器916。
加速度传感器911可以检测以终端900建立的坐标系的三个坐标轴上的加速度大小。比如，加速度传感器911可以用于检测重力加速度在三个坐标轴上的分量。处理器901可以根据加速度传感器911采集的重力加速度信号，控制显示屏905以横向视图或纵向视图进行用户界面的显示。加速度传感器911还可以用于游戏或者用户的运动数据的采集。
陀螺仪传感器912可以检测终端900的机体方向及转动角度，陀螺仪传感器912可以与加速度传感器911协同采集用户对终端900的3D动作。处理器901根据陀螺仪传感器912采集的数据，可以实现如下功能：动作感应(比如根据用户的倾斜操作来改变UI)、拍摄时的图像稳定、游戏控制以及惯性导航。
压力传感器913可以设置在终端900的侧边框和/或显示屏905的下层。当压力传感器913设置在终端900的侧边框时，可以检测用户对终端900的握持信号，由处理器901根据压力传感器913采集的握持信号进行左右手识别或快捷操作。当压力传感器913设置在显示屏905的下层时，由处理器901根据用户对显示屏905的压力操作，实现对UI界面上的可操作性控件进行控制。可操作性控件包括按钮控件、滚动条控件、图标控件、菜单控件中的至少一种。
指纹传感器914用于采集用户的指纹，由处理器901根据指纹传感器914采集到的指纹识别用户的身份，或者，由指纹传感器914根据采集到的指纹识别用户的身份。在识别出用户的身份为可信身份时，由处理器901授权该用户执行相关的敏感操作，该敏感操作包括解锁屏幕、查看加密信息、下载软件、支付及更改设置等。指纹传感器914可以被设置终端900的正面、背面或侧面。当终端900上设置有物理按键或厂商Logo时，指纹传感器914可以与物理按键或厂商Logo集成在一起。
光学传感器915用于采集环境光强度。在一个实施例中，处理器901可以根据光学传感器915采集的环境光强度，控制显示屏905的显示亮度。具体地，当环境光强度较高时，调高显示屏905的显示亮度；当环境光强度较低时，调低显示屏905的显示亮度。在另一个实施例中，处理器901还可以根据光学传感器915采集的环境光强度，动态调整摄像头组件906的拍摄参数。
接近传感器916，也称距离传感器，设置在终端900的前面板。接近传感器916用于采集用户与终端900的正面之间的距离。在一个实施例中，当接近传感器916检测到用户与终端900的正面之间的距离逐渐变小时，由处理器901控制显示屏905从亮屏状态切换为息屏状态；当接近传感器916检测到用户与终端900的正面之间的距离逐渐变大时，由处理器901控制显示屏905从息屏状态切换为亮屏状态。
本领域技术人员可以理解，图9中示出的结构并不构成对终端900的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
可选地，该电子设备提供为服务器。图10是本申请实施例提供的一种服务器的结构示意图，该服务器1000可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(central processing units，CPU)1001和一个或一个以上的存储器1002，其中，所述存储器1002中存储有至少一条程序代码，所述至少一条程序代码由所述处理器1001加载并执行以实现上述各个方法实施例提供的页岩气井EUR的确定方法。当然，该服务器还可以具有有线或无线网络接口、键盘以及输入输出接口等部件，以便进行输入输出，该服务器还可以包括其他用于实现设备功能的部件，在此不做赘述。
本申请实施例还提供了一种计算机可读存储介质，该计算机可读存储介质中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以实现上述实施例的页岩气井EUR的确定方法中执行的操作。
本申请实施例还提供了一种计算机程序，该计算机程序中存储有至少一条程序代码，该至少一条程序代码由处理器加载并执行，以实现上述实施例的页岩气井EUR的确定方法中执行的操作。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本申请的可选实施例，并不用以限制本申请，凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。
Claims (10)
1.一种页岩气井EUR的确定方法，其特征在于，所述方法包括：
通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，所述第一历史产量数据表示所述样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量；
通过所述每个样本气井的产量预测模型，预测所述每个样本气井在之后的目标时间段内的页岩气产量，所述目标时间段包括多个单位时间段，所述样本气井在所述目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量；
将所述每个样本气井的总历史产量与所述目标时间段内的页岩气产量之和，确定为所述样本气井的EUR；
基于所述多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系，所述推定累计产量是根据所述样本气井在达到边界流形态之前的产量下降速率确定的；
对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR。
2.根据权利要求1所述的方法，其特征在于，所述基于所述多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系之前，所述方法还包括：
对于所述每个样本气井，基于所述样本气井在达到边界流形态之前的产量下降速率，推定所述样本气井的产量下降为零的截止时间段；
基于所述产量下降速率，确定所述样本气井在所述截止时间段之前的推定累计产量。
3.根据权利要求2所述的方法，其特征在于，所述基于所述样本气井在达到边界流形态之前的产量下降速率，推定所述样本气井的产量下降为零的截止时间段，包括：
在坐标系中标注所述样本气井的第二历史产量数据，所述第二历史产量数据表示所述样本气井在达到边界流形态之前的多个第二历史时间段内的页岩气产量，所述坐标系的横坐标表示时间段，所述坐标系的纵坐标表示页岩气产量；
基于所述样本气井的产量达到峰值之后且所述样本气井达到边界流形态之前的多个第二历史时间段内的页岩气产量拟合出至少一条线段，延长所述至少一条线段中斜率最小的线段至所述坐标系的横轴；
将所述斜率最小的线段与所述横轴的交点的横坐标确定为所述截止时间段。
4.根据权利要求1所述的方法，其特征在于，所述通过所述每个样本气井的产量预测模型，预测所述每个样本气井在之后的目标时间段内的页岩气产量，包括：
对于所述每个样本气井，通过所述样本气井对应的产量预测模型预测所述样本气井在之后的多个单位时间段内的页岩气产量；
确定所述最低极限产量对应的第一单位时间段，将所述第一单位时间段以及所述第一单位时间段之前的多个单位时间段确定为所述目标时间段；
统计所述样本气井在所述目标时间段内的页岩气产量。
5.根据权利要求1所述的方法，其特征在于，所述基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR之前，所述方法还包括：
将验证气井的总历史产量与所述验证气井在所述目标时间段内的页岩气产量之和，确定为所述验证气井的第一EUR，所述验证气井的页岩气的流动形态达到边界流形态；
基于所述验证气井的推定累计产量和所述映射关系，确定所述验证气井的第二EUR；
将所述第一EUR和所述第二EUR的绝对差值与所述第一EUR的商，确定为所述映射关系的误差；
在所述误差小于参考误差的前提下，执行对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR的步骤。
6.根据权利要求1所述的方法，其特征在于，所述通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型之前，所述方法还包括：
从所述多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井。
7.根据权利要求6所述的方法，其特征在于，所述从所述多个页岩气井中确定生产时长不小于第一参考时长，且处于边界流形态的持续时长不小于第二参考时长的多个样本气井，包括：
在坐标系中标注生产时长不小于所述第一参考时长的每个页岩气井的历史产量数据，所述历史产量数据表示所述页岩气井在多个历史时间段内的页岩气产量，所述坐标系的横坐标表示时间段，所述坐标系的纵坐标表示页岩气产量；
确定所述页岩气井的产量趋势，若所述页岩气井的产量趋势包括至少两个产量下降速率不同的产量下降阶段，则确定产量下降速率最小的产量下降阶段的持续时长；
若所述持续时长不小于所述第二参考时长，则将所述页岩气井确定为所述样本气井。
8.一种页岩气井EUR的确定装置，其特征在于，所述装置包括：
模型建立模块，被配置为通过拟合多个样本气井的第一历史产量数据，建立每个样本气井的产量预测模型，所述第一历史产量数据表示所述样本气井在达到边界流形态之后的多个第一历史时间段内的页岩气产量；
产量预测模块，被配置为通过所述每个样本气井的产量预测模型，预测所述每个样本气井在之后的目标时间段内的页岩气产量，所述目标时间段包括多个单位时间段，所述样本气井在所述目标时间段内的最后一个单位时间段内的页岩气产量达到最低极限产量；
第一EUR确定模块，被配置为将所述每个样本气井的总历史产量与所述目标时间段内的页岩气产量之和，确定为所述样本气井的EUR；
关系建立模块，被配置为基于所述多个样本气井的推定累计产量和EUR，建立推定累计产量与EUR之间的映射关系，所述推定累计产量是根据所述样本气井在达到边界流形态之前的产量下降速率确定的；
第二EUR确定模块，被配置为对于任一未达到边界流形态的目标气井，基于所述目标气井的推定累计产量和所述映射关系，确定所述目标气井的EUR。
9.一种电子设备，其特征在于，所述电子设备包括处理器和存储器，所述存储器中存储有至少一条程序代码，所述程序代码由所述处理器加载并执行以实现如权利要求1至权利要求7任一项所述的页岩气井EUR的确定方法所执行的操作。
10.一种计算机可读存储介质，其特征在于，所述存储介质中存储有至少一条程序代码，所述程序代码由处理器加载并执行以实现如权利要求1至权利要求7任一项所述的页岩气井EUR的确定方法所执行的操作。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110002077.2A CN112832737B (zh) | 2021-01-04 | 2021-01-04 | 页岩气井eur的确定方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110002077.2A CN112832737B (zh) | 2021-01-04 | 2021-01-04 | 页岩气井eur的确定方法、装置、设备及存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112832737A true CN112832737A (zh) | 2021-05-25 |
| CN112832737B CN112832737B (zh) | 2023-10-31 |
Family
ID=75927333
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110002077.2A Active CN112832737B (zh) | 2021-01-04 | 2021-01-04 | 页岩气井eur的确定方法、装置、设备及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112832737B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114358427A (zh) * | 2022-01-07 | 2022-04-15 | 西南石油大学 | 一种预测页岩气井最终可采储量的方法 |
| CN115199240A (zh) * | 2022-08-25 | 2022-10-18 | 西南石油大学 | 一种页岩气井产量预测方法、装置及存储介质 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130346040A1 (en) * | 2012-06-20 | 2013-12-26 | Schlumberger Technology Corporation | Shale gas production forecasting |
| CN106351651A (zh) * | 2016-08-26 | 2017-01-25 | 中国石油天然气股份有限公司 | 页岩气井产能的预测方法及装置 |
| US20180202264A1 (en) * | 2016-02-08 | 2018-07-19 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
| US20180293336A1 (en) * | 2017-04-06 | 2018-10-11 | Qingfeng TAO | Forecasting ultimate recovery of oil and oil production for a multiply-fractured horizontal well |
| CN108952676A (zh) * | 2018-07-10 | 2018-12-07 | 中国石油天然气股份有限公司 | 一种页岩气藏非均质性评价方法及其装置 |
| CN109057786A (zh) * | 2018-07-27 | 2018-12-21 | 中国石油大学(华东) | 一种用于非常规油气藏可采储量估算的自适应模型 |
| CN111878072A (zh) * | 2020-07-31 | 2020-11-03 | 中国石油天然气集团有限公司 | 一种基于流态识别外推法的页岩气井eur计算方法 |
-
2021
- 2021-01-04 CN CN202110002077.2A patent/CN112832737B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130346040A1 (en) * | 2012-06-20 | 2013-12-26 | Schlumberger Technology Corporation | Shale gas production forecasting |
| US20180202264A1 (en) * | 2016-02-08 | 2018-07-19 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
| CN106351651A (zh) * | 2016-08-26 | 2017-01-25 | 中国石油天然气股份有限公司 | 页岩气井产能的预测方法及装置 |
| US20180293336A1 (en) * | 2017-04-06 | 2018-10-11 | Qingfeng TAO | Forecasting ultimate recovery of oil and oil production for a multiply-fractured horizontal well |
| CN108952676A (zh) * | 2018-07-10 | 2018-12-07 | 中国石油天然气股份有限公司 | 一种页岩气藏非均质性评价方法及其装置 |
| CN109057786A (zh) * | 2018-07-27 | 2018-12-21 | 中国石油大学(华东) | 一种用于非常规油气藏可采储量估算的自适应模型 |
| CN111878072A (zh) * | 2020-07-31 | 2020-11-03 | 中国石油天然气集团有限公司 | 一种基于流态识别外推法的页岩气井eur计算方法 |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114358427A (zh) * | 2022-01-07 | 2022-04-15 | 西南石油大学 | 一种预测页岩气井最终可采储量的方法 |
| CN114358427B (zh) * | 2022-01-07 | 2024-05-07 | 西南石油大学 | 一种预测页岩气井最终可采储量的方法 |
| CN115199240A (zh) * | 2022-08-25 | 2022-10-18 | 西南石油大学 | 一种页岩气井产量预测方法、装置及存储介质 |
| CN115199240B (zh) * | 2022-08-25 | 2023-05-19 | 西南石油大学 | 一种页岩气井产量预测方法、装置及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112832737B (zh) | 2023-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108304265B (zh) | 内存管理方法、装置及存储介质 | |
| CN113204298B (zh) | 显示发布进度的方法、装置、电子设备及存储介质 | |
| CN111127509B (zh) | 目标跟踪方法、装置和计算机可读存储介质 | |
| CN111005715B (zh) | 确定气井产量的方法、装置及存储介质 | |
| CN108764530B (zh) | 油井抽油机的工作参数配置方法及装置 | |
| CN112832737B (zh) | 页岩气井eur的确定方法、装置、设备及存储介质 | |
| CN111255434B (zh) | 气井的试井方法、装置和计算机存储介质 | |
| CN110909264A (zh) | 信息处理方法、装置、设备及存储介质 | |
| CN110166275B (zh) | 信息处理方法、装置及存储介质 | |
| CN110990728B (zh) | 兴趣点信息的管理方法、装置、设备及存储介质 | |
| CN109281648B (zh) | 确定油藏的合理井网密度的方法和装置 | |
| CN109036463B (zh) | 获取歌曲的难度信息的方法、装置及存储介质 | |
| CN109899051B (zh) | 油井设备的评价标准确定方法、装置及存储介质 | |
| CN111063372B (zh) | 确定音高特征的方法、装置、设备及存储介质 | |
| CN109343113B (zh) | 一种预测油气藏位置的方法和装置 | |
| CN111080033B (zh) | 预测服务容量的方法、装置、计算机设备及存储介质 | |
| CN109902844B (zh) | 注水系统的优化信息确定方法、装置及存储介质 | |
| CN112132472A (zh) | 资源管理方法、装置、电子设备及计算机可读存储介质 | |
| CN111274662B (zh) | 对流连通老腔的筛选方法及装置 | |
| CN111219206B (zh) | 盐穴储气库的造腔方法及装置 | |
| CN111736215A (zh) | 断层断距确定方法及装置 | |
| CN111861780A (zh) | 生产井中生产层的产量确定方法、装置及存储介质 | |
| CN115263244B (zh) | 气水层的排水控制方法、装置及计算机存储介质 | |
| CN111798358B (zh) | 确定算路时间的方法、装置、电子设备及可读存储介质 | |
| CN112308587B (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 |