CN115263284A - 页岩气层的筛选方法、装置及存储介质 - Google Patents
页岩气层的筛选方法、装置及存储介质 Download PDFInfo
- Publication number
- CN115263284A CN115263284A CN202110484937.0A CN202110484937A CN115263284A CN 115263284 A CN115263284 A CN 115263284A CN 202110484937 A CN202110484937 A CN 202110484937A CN 115263284 A CN115263284 A CN 115263284A
- Authority
- CN
- China
- Prior art keywords
- shale gas
- gas layer
- porosity
- target shale
- formation
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000012216 screening Methods 0.000 title claims abstract description 19
- 238000003860 storage Methods 0.000 title claims abstract description 15
- 239000011148 porous material Substances 0.000 claims abstract description 128
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 125
- 238000006243 chemical reaction Methods 0.000 claims description 67
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 64
- 229910052799 carbon Inorganic materials 0.000 claims description 64
- 238000005070 sampling Methods 0.000 claims description 33
- 238000012545 processing Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 2
- 239000004215 Carbon black (E152) Substances 0.000 abstract description 9
- 229930195733 hydrocarbon Natural products 0.000 abstract description 9
- 150000002430 hydrocarbons Chemical class 0.000 abstract description 9
- 230000002159 abnormal effect Effects 0.000 abstract description 5
- 238000005755 formation reaction Methods 0.000 description 99
- 230000002093 peripheral effect Effects 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 238000005056 compaction Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 6
- 229910052770 Uranium Inorganic materials 0.000 description 5
- 238000004891 communication Methods 0.000 description 5
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000013473 artificial intelligence Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000002706 hydrostatic effect 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
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 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
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000011160 research Methods 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- 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
-
- 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
- E21B47/06—Measuring temperature or pressure
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请实施例公开了一种页岩气层的筛选方法、装置及存储介质，属于勘探技术领域。该方法包括：根据目标页岩气层对应的测井资料，确定目标页岩气层的修正孔隙度；基于目标页岩气层的修正孔隙度，确定目标页岩气层的地层孔隙压力系数；当地层孔隙压力系数大于参考数值时，将目标页岩气层确定为待开发的页岩气层。本申请实施例中，由于确定的目标页岩气层的修正孔隙度是指除了有机孔之外的其他孔孔隙度，从而能够消除目标页岩气层中因有机质生烃排烃而产生的异常高压的影响。这样通过修正孔隙度确定的目标页岩气层的地层孔隙压力系数会更符合实际情况，从而根据地层孔隙压力系数才能更准确的确定的目标页岩气层是否为待开发层。
Description
技术领域
本申请实施例涉及勘探技术领域，特别涉及一种页岩气层的筛选方法、装置及存储介质。
背景技术
地层孔隙压力系数是反映页岩气藏保存条件好坏的重要参数，也是评价一个地区页岩气资源前景的重要依据，因此，研究页岩气层的地层孔隙压力系数对页岩气勘探开发具有十分重要的意义。示例地，可根据页岩气层的地层孔隙压力系数选择待开发的页岩气层。其中，地层孔隙压力系数是指地层孔隙中流体所承受的应力与静水压力的比值。
相关技术中，通过测井曲线得到正常压实趋势线，进而采用正常压实趋势线法确定页岩气层的地层孔隙压力系数。在正常压实情况下，泥页岩的垂直有效应力随着埋深的增加而逐渐增大，孔隙度减小。而当泥页岩地层发生不平衡压实时，孔隙压力增大，同时有效应力降低，孔隙度增加。这样，可根据泥页岩地层孔隙度的异常变化预测地层孔隙压力系数，进而根据地层孔隙压力系数指导页岩气层的开发。
然而，一方面，页岩气层中矿物组分的差异很容易造成测井曲线出现偏差，进而造成正产压实趋势线出现偏差；另一方面，部分泥页岩的异常高压不再以不平衡压实为主，还可能因泥页岩中有机质生烃排烃形成的异常高而产生影响，因此通过上述方法确定的地层孔隙压力系数存在偏差。
发明内容
本申请实施例提供了一种页岩气层的筛选方法、装置及存储介质，能够有效消除有机质生烃排烃产生的影响，提高待开发页岩气层确定的准确性。所述技术方案如下：
第一方面，提供了一种页岩气层的筛选方法，所述方法包括：
根据目标页岩气层对应的测井资料，确定所述目标页岩气层的修正孔隙度，所述修正孔隙度是指所述目标页岩气层中除有机孔之外的其他孔的孔隙度；
基于所述目标页岩气层的修正孔隙度，确定所述目标页岩气层的地层孔隙压力系数；
当所述地层孔隙压力系数大于参考数值时，将所述目标页岩气层确定为待开发的页岩气层。
可选地，所述根据目标页岩气层对应的测井资料，确定所述目标页岩气层的修正孔隙度，包括：
根据所述测井资料确定所述目标页岩气层的总有机碳含量和总孔隙度；
根据所述总有机碳含量对所述总孔隙度进行修正，得到所述目标页岩气层的修正孔隙度。
可选地，所述根据所述总有机碳含量对所述总孔隙度进行修正，得到所述目标页岩气层的修正孔隙度，包括：
根据所述总有机碳含量，确定所述目标页岩气层中有机孔的总孔隙度；
将所述总孔隙度和所述有机孔的总孔隙度的差值确定为所述目标页岩气层的修正孔隙度。
可选地，所述根据所述总有机碳含量，确定所述目标页岩气层中有机孔的总孔隙度，包括：
基于所述总有机碳含量，根据第一转换公式确定所述目标页岩气层中有机孔的总孔隙度；
所述第一转换公式：PORO＝a1×TOC
其中，在所述第一转换公式中，PORO是指所述目标页岩气层中有机孔的总孔隙度；TOC是指所述总有机碳含量；a1是指第一参数。
可选地，所述基于所述总有机碳含量，根据第一转换公式确定所述目标页岩气层中有机孔的总孔隙度之前，还包括：
获取多个岩心样品中每个岩心样品的有机碳含量和有机孔孔隙度，所述多个岩心样品是在不同深度的页岩气层进行取样得到的，且所述多个岩心样品的取样深度均不同；
对所述多个岩心样品的有机碳含量和有机孔孔隙度进行拟合，得到所述第一转换公式。
可选地，所述基于所述目标页岩气层的修正孔隙度，确定所述目标页岩气层的地层孔隙压力系数，包括：
基于所述目标页岩气层修正孔隙度确定所述目标页岩气层的地层有效应力；
对所述目标页岩气层的地层有效应力进行转换，得到所述目标页岩气层的地层孔隙压力系数。
可选地，所述基于所述目标页岩气层修正孔隙度确定所述目标页岩气层的地层有效应力，包括：
基于所述目标页岩气层修正孔隙度，根据第二转换公式确定所述目标页岩气层的地层有效应力；
所述第二转换公式：
其中，在所述第二转换公式中，σ是指所述目标页岩气层的地层有效应力；POR是指所述目标页岩气层的修正孔隙度；a2是指第二参数；a3是指第三参数。
可选地，所述基于所述目标页岩气层修正孔隙度，根据第二转换公式确定所述目标页岩气层的地层有效应力之前，还包括：
获取多个岩心样品的取样深度，所述多个岩心样品是在不同深度的页岩气层进行取样得到的，且所述多个岩心样品的取样深度均不同；
基于所述多个岩心样品的取样深度，确定每个岩心样品的上覆岩层压力，得到多个上覆岩层压力；
从零压力到最大上覆岩层压力的一半之间的范围内选择多个有效应力；
基于所述多个有效应力获取每个岩心样品对应的多个孔隙度；
对所述多个岩心样品的孔隙度进行归一化处理，以及对每个岩心样品的多个归一化孔隙度进行均值处理，得到与所述多个有效应力一一对应的多个平均孔隙度；
对所述多个有效应力和所述多个平均孔隙度进行指数曲线拟合，得到所述第二转换公式。
第二方面，提供了一种页岩气层的筛选装置，所述装置包括：
第一确定模块，用于根据目标页岩气层对应的测井资料，确定所述目标页岩气层的修正孔隙度，所述修正孔隙度是指所述目标页岩气层中除有机孔之外的其他孔的孔隙度；
第二确定模块，用于基于所述目标页岩气层的修正孔隙度，确定所述目标页岩气层的地层孔隙压力系数；
第三确定模块，用于当所述地层孔隙压力系数大于参考数值时，将所述目标页岩气层确定为待开发的页岩气层。
可选地，所述第一确定模块包括：
第一确定单元，用于根据所述测井资料确定所述目标页岩气层的总有机碳含量和总孔隙度；
修正单元，用于根据所述总有机碳含量对所述总孔隙度进行修正，得到所述目标页岩气层的修正孔隙度。
可选地，所述修正单元主要用于：
根据所述总有机碳含量，确定所述目标页岩气层中有机孔的总孔隙度；
将所述总孔隙度和所述有机孔的总孔隙度的差值确定为所述目标页岩气层的修正孔隙度。
可选地，所述修正单元还用于：
基于所述总有机碳含量，根据第一转换公式确定所述目标页岩气层中有机孔的总孔隙度；
所述第一转换公式：PORO＝a1×TOC
其中，在所述第一转换公式中，PORO是指所述目标页岩气层中有机孔的总孔隙度；TOC是指所述总有机碳含量；a1是指第一参数。
可选地，所述修正单元还用于：
获取多个岩心样品中每个岩心样品的有机碳含量和有机孔孔隙度，所述多个岩心样品是在不同深度的页岩气层进行取样得到的，且所述多个岩心样品的取样深度均不同；
对所述多个岩心样品的有机碳含量和有机孔孔隙度进行拟合，得到所述第一转换公式。
可选地，所述第二确定模块包括：
第二确定单元，用于基于所述目标页岩气层修正孔隙度确定所述目标页岩气层的地层有效应力；
转换单元，用于对所述目标页岩气层的地层有效应力进行转换，得到所述目标页岩气层的地层孔隙压力系数。
可选地，所述第二确定单元主要用于：
基于所述目标页岩气层修正孔隙度，根据第二转换公式确定所述目标页岩气层的地层有效应力；
所述第二转换公式：
其中，在所述第二转换公式中，σ是指所述目标页岩气层的地层有效应力；POR是指所述目标页岩气层的修正孔隙度；a2是指第二参数；a3是指第三参数。
可选地，所述第二确定单元主要用于：
获取多个岩心样品的取样深度，所述多个岩心样品是在不同深度的页岩气层进行取样得到的，且所述多个岩心样品的取样深度均不同；
基于所述多个岩心样品的取样深度，确定每个岩心样品的上覆岩层压力，得到多个上覆岩层压力；
从零压力到最大上覆岩层压力的一半之间的范围内选择多个有效应力；
基于所述多个有效应力获取每个岩心样品对应的多个孔隙度；
对所述多个岩心样品的孔隙度进行归一化处理，以及对每个岩心样品的多个归一化孔隙度进行均值处理，得到与所述多个有效应力一一对应的多个平均孔隙度；
对所述多个有效应力和所述多个平均孔隙度进行指数曲线拟合，得到所述第二转换公式。
第三方面，提供了一种计算机可读存储介质，所述存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现上述第一方面提供的任一所述的方法。
本申请实施例提供的技术方案带来的有益效果至少包括：
本申请实施例中，由于确定的目标页岩气层的修正孔隙度是指除了有机孔之外的其他孔孔隙度，从而能够消除目标页岩气层中有机质生烃排烃而产生的异常高压的影响。这样通过修正孔隙度确定的目标页岩气层的地层孔隙压力系数会更符合实际情况，避免了采用正常压实趋势线法时页岩气层中矿物组分差异引起的测井曲线不能准确评价页岩气储层的地层孔隙压力，进而根据地层孔隙压力系数才能更准确的确定的目标页岩气层是否为待开发层。
附图说明
为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请实施例的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的一种页岩气层的筛选方法的流程示意图；
图2是本申请实施例提供的一种总有机碳含量和有机孔孔隙度的关系示意图；
图3是本申请实施例提供的一种有效应力和平均孔隙度的关系示意图；
图4是本申请实施例提供的一种两种方式确定的地层孔隙压力系数的关系示意图；
图5是本申请实施例提供的一种页岩气层的筛选装置的结构示意图；
图6是本申请实施例提供的另一种页岩气层的筛选装置的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合附图对本申请实施例的实施方式作进一步地详细描述。
图1是本申请实施例提供的一种页岩气层的筛选方法的流程示意图，该方法应用于终端，该终端为智能手机、笔记本电脑、台式电脑等。参见图1，该方法包括如下步骤。
本申请实施例中，由于地层孔隙压力系数是反映页岩气藏保存条件好坏的重要参数，可根据目标页岩气层的地层孔隙压力系数确定目标页岩气层是否适合开发，也即是确定目标页岩气层是否具有开发价值。因此，在确定目标页岩气层是否适合开发之前需要确定目标页岩气层的地层孔隙压力系数。由于目标页岩气层的地层孔隙压力系数取决于目标页岩气层的属性信息和组分信息，而目标页岩气层的属性信息和组分信息可根据目标页岩气层对应的测井资料获取。由此，可以先获取目标页岩气层对应的测井资料，之后按照如下步骤101-步骤102确定目标页岩气层的地层孔隙压力系数。
其中，目标页岩气层的属性信息包括目标页岩气层的密度值和声波时差，目标页岩气层的组分信息包括铀含量值。
步骤101：根据目标页岩气层对应的测井资料，确定目标页岩气层的修正孔隙度。
其中，修正孔隙度是指目标页岩气层中除有机孔之外的其他孔的孔隙度。
在一些实施例中，根据测井资料确定目标页岩气层的总有机碳含量和总孔隙度；根据总有机碳含量对总孔隙度进行修正，得到目标页岩气层的修正孔隙度。
对于目标页岩气层的总有机碳含量，先获取目标页岩气层所在气藏区块的总有机碳含量和测井曲线拟合建立的第一经验公式，之后根据测井资料获取目标页岩气层的密度值和铀含量值，再将目标页岩气层的密度值和铀含量值带入第一经验公式即可得到目标页岩气层的总有机碳含量。
其中，第一经验公式如下：TOC＝a×DEN+b×log(U)+c
其中，在第一经验公式中，TOC是指目标页岩气层的总有机碳含量；DEN是指目标页岩气层的密度值；U是指目标页岩气层中油含量值；a是指第一计算参数；b是指第二计算参数；c是指第三计算参数。不同气藏区块a、b、c的取值不同，示例地，目标页岩气层所在气藏区块的a的取值为-16.89，b的取值为0.294，c的取值为45.367，确定的总有机碳含量为3.5。
需要说明的是，除了通过上述方式确定目标页岩气层的总有机碳含量之外，还能够通过其他方式确定目标页岩气层的总有机碳含量。示例地，目标页岩气层的测井资料还包括通过元素俘获测井的方法得到的非弹性伽马谱，此时对非弹性伽马谱进行解谱，即可得到目标页岩气层的总有机碳含量。其中，通过元素俘获测井资料得到目标页岩气层中的总有机碳含量的具体方式可参考相关技术，本申请实施例对此不做限定。
对于目标页岩气层的总孔隙度，先获取目标页岩气层所在气藏区块的总孔隙度和测井曲线拟合建立的第二经验公式，之后根据测井资料获取目标页岩气层的密度值、声波时差和铀含量值，再将目标页岩气层的密度值、声波时差和铀含量值带入第二经验公式即可得到目标页岩气层的总孔隙度。
其中，第二经验公式如下：POR＝d×DEN+f×AC+g×log(U)+h
其中，在第二经验公式中，POR是指目标页岩气层的总孔隙度；DEN是指目标页岩气层的密度值；AC是指目标页岩气层的声波时差；U是指目标页岩气层中油含量值；d是指第四计算参数；f是指第五计算参数；g是指第六计算参数；h是指第七计算参数。不同气藏区块d、f、g、h的取值不同，示例地，目标页岩气层所在气藏区块的d的取值为-6.51，f的取值为0.124，g的取值为0.093，h的取值为12.98，确定的总孔孔隙度为5.5。
需要说明的是，除了通过上述方式确定目标页岩气层的总孔隙度之外，还能够通过其他方式确定目标页岩气层的总孔隙度。示例地，根据岩石体积模型采用最优化求解目标页岩气层的总孔隙度。其中，采用最优化求解目标页岩气层的总孔隙度的方式可参考相关技术的描述，本申请实施例对此不在赘述。
对于目标页岩气层的修正孔隙度，则可先根据总有机碳含量，确定目标页岩气层中有机孔的总孔隙度；再将总孔隙度和有机孔的总孔隙度的差值确定为目标页岩气层的修正孔隙度。
本申请实施例中，有机质生烃排烃后气体会聚集在有机孔内，而有机质的含量可通过有机碳的含量表示，因此，有机碳的含量和有机孔的总孔隙度之间存在相关性，也即是由于有机孔形成的多少取决于有机碳含量的多少。
因此在一些实施例中，可以先确定总有机碳含量和有机孔的总孔隙度之间的第一转换公式，之后基于总有机碳含量，根据第一转换公式确定目标页岩气层中有机孔的总孔隙度。
在确定了目标页岩气层中有机孔的总孔隙度之后，则可将总孔隙度和有机孔的总孔隙度的差值确定为目标页岩气层的修正孔隙度，以消除目标页岩气层中因有机质生烃排烃而产生的异常高压的影响。
需要说明的是，可通过如下步骤(1)-(2)确定总有机碳含量和有机孔的总孔隙度之间的第一转换公式。
(1)、获取多个岩心样品中每个岩心样品的有机碳含量和有机孔孔隙度。
其中，多个岩心样品是在不同深度的页岩气层进行取样得到的，且多个岩心样品的取样深度均不同。另外，为了保证确定的第一转换公式的准确性，获取的多个岩心样品的总有机碳含量均不同，且多个岩心样品的有机孔孔隙度也均不同。
可选地，对于不同深度的多个页岩气层，可在每个页岩气层获取至少一个岩心样品，之后获取每个岩心样品的有机碳含量和有机孔孔隙度，之后从在多个页岩气层获取的岩心样品中选择有机碳含量不同，且有机孔孔隙度不同的多个岩心样品。
其中，每个岩心样品的有机碳含量和有机孔孔隙度可在实验条件下获取，示例地，可在场发射扫描电子显微镜下进行切片观察获得，具体获得的方式可参考相关技术，本申请实施例对此不做限定。
(2)、对多个岩心样品的有机碳含量和有机孔孔隙度进行拟合，得到第一转换公式。
其中，第一转换公式：PORO＝a1×TOC
在第一转换公式中，PORO是指目标页岩气层中有机孔的总孔隙度；TOC是指目标页岩气层的总有机碳含量；a1是指第一参数。不同气藏区块a1的取值不同，示例地，取样得到24个岩心样品，建立直角坐标系，并以有机碳含量作为横轴，以有机孔孔隙度作为纵轴，如图2所示，将每个岩心样品的有机碳含量和有机孔孔隙度标注在直角坐标系中，得到24个坐标点，之后以最小二乘法对这24个点进行拟合，得到a1的取值为0.2068的第一转换公式，根据第一转换公式确定的目标页岩气层的有机孔的总孔隙度为0.724，此时对目标页岩气层的总孔隙度进行修正后得到的修正孔隙度为4.776。
步骤102：基于目标页岩气层的修正孔隙度，确定目标页岩气层的地层孔隙压力系数。
可选地，基于目标页岩气层修正孔隙度确定目标页岩气层的地层有效应力；对目标页岩气层的地层有效应力进行转换，得到目标页岩气层的地层孔隙压力系数。
由于孔隙度与地层有效应力之间存在相关性，因此可以先确定孔隙度与地层有效应力之间的第二转换公式，之后基于目标页岩气层修正孔隙度，根据第二转换公式确定目标页岩气层的地层有效应力。
在得到目标页岩气层的地层有效应力后，可通过第三转换公式将目标页岩气层的地层有效应力转换为目标地层的孔隙压力系数。
其中，第三转换公式：
在第三转换公式中，PP是指目标页岩气层的地层孔隙压力系数；σ是指目标页岩气层的地层有效应力；H是指目标页岩气层的深度；ρ是指不同深度地层的密度值；g是指重力单位，取常数9.8牛顿/千克。
需要说明的是，可通过如下方式确定总有机碳含量和有机孔的总孔隙度之间的第二转换公式。
获取多个岩心样品的取样深度，多个岩心样品是在不同深度的页岩气层进行取样得到的，且多个岩心样品的取样深度均不同；基于多个岩心样品的取样深度，确定每个岩心样品的上覆岩层压力，得到多个上覆岩层压力；从零压力到最大上覆岩层压力的一半之间的范围内选择多个有效应力；基于多个有效应力获取每个岩心样品对应的多个孔隙度；对多个岩心样品的孔隙度进行归一化处理，以及对每个岩心样品的多个归一化孔隙度进行均值处理，得到与多个有效应力一一对应的多个平均孔隙度；对多个有效应力和多个平均孔隙度进行指数曲线拟合，得到第二转换公式。
其中，对于每个岩心样品的上覆岩层压力，可根据该岩心样品的取样深度在测井资料中获取对应深度的测井数据，进而基于测井数据确定该岩心样品的上覆岩层压力。基于测井数据获取上覆岩层压力的具体方式可参考相关技术，本申请实施例对此不做限定。
其中，从零压力到最大上覆岩层压力的一半之间的范围内可均匀的选择4～8个压力值作为选取的有效应力。之后对于每个岩心样品在每个有效应力下的孔隙度的获取方法，可参考《SY/T6385-2016覆压下岩石孔隙度和渗透率的测定方法》，本申请实施例对此不做限定。
其中，第二转换公式：
在第二转换公式中，σ是指目标页岩气层的地层有效应力；POR是指目标页岩气层的修正孔隙度；a2是指第二参数；a3是指第三参数。不同气藏区块a2和a3的取值均不同，示例地，取样得到24个岩心样品，且获取每个岩心样品的8个有效应力，以及每个有效应力对应的平均孔隙度，之后建立直角坐标系，并以平均孔隙度作为横轴，以有效应力作为纵轴，将8个有效应力和分别对应的平均孔隙度标注在直角坐标系中，如图3所示，得到8个坐标点，之后进行拟合，得到a2的取值为0.214.61、a3的取值为-0.482的第二转换公式；目标页岩气层的深度为2505米，通过第三转换公式对目标页岩气层的地层有效应力进行转换后，得到的目标页岩气层的地层孔隙压力系数为1.685。
上述第一转换公式和第二转换公式可同时获取，也可先获取第一转换公式再获取第二转换公式；或者先获取第二转换公式再获取第一转换公式。由于在获取第一转换公式和第二转换公式时，均需要获取多个岩心样品，由此为了避免岩心样品取样的重复性，可一次性获取多个岩心样品，之后根据获取的多个岩心样品同时获取第一转换公式和第二转换公式。
实际使用过程中，对于24个不同的页岩气层，通过本申请实施例提供的方式确定得到24个地层孔隙压力系数；而通过电缆地层测试技术(XPT)获取到24个地层孔隙压力系数，以得到24个页岩气层中每个页岩气层的两个地层孔隙压力系数。建立坐标系，并以通过电缆地层测试技术获取的地层孔隙压力系数为横坐标，以本申请实施例提供的方式确定的地层孔隙压力系数为纵坐标，并将每个页岩气层的地层孔隙压力系数按照获取方式标注在直角坐标系中，如图4所示。此时通过本申请实施例提供的方式确定得到24个地层孔隙压力系数与通过电缆地层测试技术(XPT)获取到24个地层孔隙压力系数之间的相关系数近似等于1，从而表明通过本申请实施例确定的地层孔隙压力系数的准确性，且相对于通过电缆地层测试技术，本申请实施例提供的方式的方便性和快捷性。
步骤103：当地层孔隙压力系数大于参考数值时，将目标页岩气层确定为待开发的页岩气层。
结合上述步骤102确定了目标页岩气层的地层孔隙压力系数后，由于地层孔隙压力系数越大，则表明地层孔隙中所含有的气体越多，由此可根据地层孔隙压力系数的大小确定目标页岩气层是否为待开发的页岩气层。
可选地，将目标页岩气层的地层孔隙压力系数与参考数值进行比较，若地层孔隙压力系数大于参考数值，则将目标页岩气层确定为待开发的页岩气层；相反，若地层孔隙压力系数小于或等于参考数值，则确定目标页岩气层不适合开发。
示例地，参考数值为1.2或1.5。1.2或1.5只是参考数值的一种实例，参考数值也可以为其他数值，本申请实施例对此不做限定。
本申请实施例中，由于确定的目标页岩气层的修正孔隙度是指除了有机孔之外的其他孔孔隙度，从而能够消除目标页岩气层中因有机质生烃排烃而产生的异常高压的影响。这样通过修正孔隙度确定的目标页岩气层的地层孔隙压力系数会更符合实际情况，从而根据地层孔隙压力系数才能更准确的确定的目标页岩气层是否为待开发层。另外，本申请实施例中不需要特殊仪器测量或者开展小型压裂测试，使得目标页岩气层的地层孔隙压力系数的确定更为简化且高效。
图5是本申请实施例提供的一种页岩气层的筛选装置的结构示意图，该装置集成在智能手机、笔记本电脑、台式电脑等中。参见图5，该装置包括：
第一确定模块501，用于根据目标页岩气层对应的测井资料，确定目标页岩气层的修正孔隙度，修正孔隙度是指目标页岩气层中除有机孔之外的其他孔的孔隙度；
第二确定模块502，用于基于目标页岩气层的修正孔隙度，确定目标页岩气层的地层孔隙压力系数；
第三确定模块503，用于当地层孔隙压力系数大于参考数值时，将目标页岩气层确定为待开发的页岩气层。
可选地，第一确定模块501包括：
第一确定单元，用于根据测井资料确定目标页岩气层的总有机碳含量和总孔隙度；
修正单元，用于根据总有机碳含量对总孔隙度进行修正，得到目标页岩气层的修正孔隙度。
可选地，修正单元主要用于：
根据总有机碳含量，确定目标页岩气层中有机孔的总孔隙度；
将总孔隙度和有机孔的总孔隙度的差值确定为目标页岩气层的修正孔隙度。
可选地，修正单元还用于：
基于总有机碳含量，根据第一转换公式确定目标页岩气层中有机孔的总孔隙度；
第一转换公式：PORO＝a1×TOC
其中，在第一转换公式中，PORO是指目标页岩气层中有机孔的总孔隙度；TOC是指目标页岩气层的总有机碳含量；a1是指第一参数。
可选地，修正单元还用于：
获取多个岩心样品中每个岩心样品的有机碳含量和有机孔孔隙度，多个岩心样品是在不同深度的页岩气层进行取样得到的，且多个岩心样品的取样深度均不同；
对多个岩心样品的有机碳含量和有机孔孔隙度进行拟合，得到第一转换公式。
可选地，第二确定模块502包括：
第二确定单元，用于基于目标页岩气层修正孔隙度确定目标页岩气层的地层有效应力；
转换单元，用于对目标页岩气层的地层有效应力进行转换，得到目标页岩气层的地层孔隙压力系数。
可选地，第二确定单元主要用于：
基于目标页岩气层修正孔隙度，根据第二转换公式确定目标页岩气层的地层有效应力；
第二转换公式：
其中，在第二转换公式中，σ是指目标页岩气层的地层有效应力；POR是指目标页岩气层的修正孔隙度；a2是指第二参数；a3是指第三参数。
可选地，第二确定单元主要用于：
获取多个岩心样品的取样深度，多个岩心样品是在不同深度的页岩气层进行取样得到的，且多个岩心样品的取样深度均不同；
基于多个岩心样品的取样深度，确定每个岩心样品的上覆岩层压力，得到多个上覆岩层压力；
从零压力到最大上覆岩层压力的一半之间的范围内选择多个有效应力；
基于多个有效应力获取每个岩心样品对应的多个孔隙度；
对多个岩心样品的孔隙度进行归一化处理，以及对每个岩心样品的多个归一化孔隙度进行均值处理，得到与多个有效应力一一对应的多个平均孔隙度；
对多个有效应力和多个平均孔隙度进行指数曲线拟合，得到第二转换公式。
本申请实施例中，由于确定的目标页岩气层的修正孔隙度是指除了有机孔之外的其他孔孔隙度，从而能够消除目标页岩气层中因有机质生烃排烃而产生的异常高压的影响。这样通过修正孔隙度确定的目标页岩气层的地层孔隙压力系数会更符合实际情况，从而根据地层孔隙压力系数才能更准确的确定的目标页岩气层是否为待开发层。
需要说明的是：上述实施例提供的页岩气层的筛选装置在筛选待开发页岩气层时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将设备的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的页岩气层的筛选装置与页岩气层的筛选方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
图6示出了本申请实施例一个示例性实施例提供的终端600的结构框图。参见图6，该终端600可以是：智能手机、平板电脑、笔记本电脑或台式电脑。终端600还可能被称为用户设备、便携式终端、膝上型终端、台式终端等其他名称。参见图6，终端600可以包括处理器601和存储器602。
处理器601可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器601可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field－Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器601也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中，处理器601可以在集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器601还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。
存储器602可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器602还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器602中的非暂态的计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器601所执行以实现本申请中方法实施例提供的一种页岩气层的筛选方法。
在一些实施例中，终端600还可选包括有：外围设备接口603和至少一个外围设备。处理器601、存储器602和外围设备接口603之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口603相连。具体地，外围设备包括：射频电路604、显示屏605、定位组件606和电源607中的至少一种。
外围设备接口603可被用于将I/O(Input/Output，输入/输出)相关的至少一个外围设备连接到处理器601和存储器602。在一些实施例中，处理器601、存储器602和外围设备接口603被集成在同一芯片或电路板上；在一些其他实施例中，处理器601、存储器602和外围设备接口603中的任意一个或两个可以在单独的芯片或电路板上实现，本实施例对此不加以限定。
射频电路604用于接收和发射RF(Radio Frequency，射频)信号，也称电磁信号。射频电路604通过电磁信号与通信网络以及其他通信设备进行通信。射频电路604将电信号转换为电磁信号进行发送，或者，将接收到的电磁信号转换为电信号。可选地，射频电路604包括：天线系统、RF收发器、一个或多个放大器、调谐器、振荡器、数字信号处理器、编解码芯片组、用户身份模块卡等等。射频电路604可以通过至少一种无线通信协议来与其它终端进行通信。该无线通信协议包括但不限于：万维网、城域网、内联网、各代移动通信网络(2G、3G、4G及5G)、无线局域网和/或WiFi(Wireless Fidelity，无线保真)网络。在一些实施例中，射频电路604还可以包括NFC(Near Field Communication，近距离无线通信)有关的电路，本申请对此不加以限定。
显示屏605用于显示UI(User Interface，用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏605是显示屏时，显示屏605还具有采集在显示屏605的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器601进行处理。此时，显示屏605还可以用于提供虚拟按钮和/或虚拟键盘，也称软按钮和/或软键盘。在一些实施例中，显示屏605可以为一个，设置终端600的前面板；在另一些实施例中，显示屏605可以为至少两个，分别设置在终端600的不同表面或呈折叠设计；在再一些实施例中，显示屏605可以是柔性显示屏，设置在终端600的弯曲表面上或折叠面上。甚至，显示屏605还可以设置成非矩形的不规则图形，也即异形屏。显示屏605可以采用LCD(LiquidCrystal Display，液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
定位组件606用于定位终端600的当前地理位置，以实现导航或LBS(LocationBased Service，基于位置的服务)。定位组件606可以是基于美国的GPS(GlobalPositioning System，全球定位系统)、中国的北斗系统或俄罗斯的伽利略系统的定位组件。
电源607用于为终端600中的各个组件进行供电。电源607可以是交流电、直流电、一次性电池或可充电电池。当电源607包括可充电电池时，该可充电电池可以是有线充电电池或无线充电电池。有线充电电池是通过有线线路充电的电池，无线充电电池是通过无线线圈充电的电池。该可充电电池还可以用于支持快充技术。
本领域技术人员可以理解，图6中示出的结构并不构成对终端600的限定，可以包括比图示更多或更少的组件，或者组合某些组件，或者采用不同的组件布置。
在上述实施例中，还提供了一种包括指令的非暂态的计算机可读存储介质，该计算机可读存储介质用于存储至少一个指令，该至少一个指令用于被处理器所执行以实现上述图1所示实施例提供的方法。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本申请实施例的可选实施例，并不用以限制本申请实施例，凡在本申请实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请实施例的保护范围之内。
Claims (10)
1.一种页岩气层的筛选方法，其特征在于，所述方法包括：
根据目标页岩气层对应的测井资料，确定所述目标页岩气层的修正孔隙度，所述修正孔隙度是指所述目标页岩气层中除有机孔之外的其他孔的孔隙度；
基于所述目标页岩气层的修正孔隙度，确定所述目标页岩气层的地层孔隙压力系数；
当所述地层孔隙压力系数大于参考数值时，将所述目标页岩气层确定为待开发的页岩气层。
2.如权利要求1所述的方法，其特征在于，所述根据目标页岩气层对应的测井资料，确定所述目标页岩气层的修正孔隙度，包括：
根据所述测井资料确定所述目标页岩气层的总有机碳含量和总孔隙度；
根据所述总有机碳含量对所述总孔隙度进行修正，得到所述目标页岩气层的修正孔隙度。
3.如权利要求2所述的方法，其特征在于，所述根据所述总有机碳含量对所述总孔隙度进行修正，得到所述目标页岩气层的修正孔隙度，包括：
根据所述总有机碳含量，确定所述目标页岩气层中有机孔的总孔隙度；
将所述总孔隙度和所述有机孔的总孔隙度的差值确定为所述目标页岩气层的修正孔隙度。
4.如权利要求3所述的方法，其特征在于，所述根据所述总有机碳含量，确定所述目标页岩气层中有机孔的总孔隙度，包括：
基于所述总有机碳含量，根据第一转换公式确定所述目标页岩气层中有机孔的总孔隙度；
所述第一转换公式：PORO＝a1×TOC
其中，在所述第一转换公式中，PORO是指所述目标页岩气层中有机孔的总孔隙度；TOC是指所述总有机碳含量；a1是指第一参数。
5.如权利要求4所述的方法，其特征在于，所述基于所述总有机碳含量，根据第一转换公式确定所述目标页岩气层中有机孔的总孔隙度之前，还包括：
获取多个岩心样品中每个岩心样品的有机碳含量和有机孔孔隙度，所述多个岩心样品是在不同深度的页岩气层进行取样得到的，且所述多个岩心样品的取样深度均不同；
对所述多个岩心样品的有机碳含量和有机孔孔隙度进行拟合，得到所述第一转换公式。
6.如权利要求1所述的方法，其特征在于，所述基于所述目标页岩气层的修正孔隙度，确定所述目标页岩气层的地层孔隙压力系数，包括：
基于所述目标页岩气层修正孔隙度确定所述目标页岩气层的地层有效应力；
对所述目标页岩气层的地层有效应力进行转换，得到所述目标页岩气层的地层孔隙压力系数。
7.如权利要求6所述的方法，其特征在于，所述基于所述目标页岩气层修正孔隙度确定所述目标页岩气层的地层有效应力，包括：
基于所述目标页岩气层修正孔隙度，根据第二转换公式确定所述目标页岩气层的地层有效应力；
所述第二转换公式：
其中，在所述第二转换公式中，σ是指所述目标页岩气层的地层有效应力；POR是指所述目标页岩气层的修正孔隙度；a2是指第二参数；a3是指第三参数。
8.如权利要求7所述的方法，其特征在于，所述基于所述目标页岩气层修正孔隙度，根据第二转换公式确定所述目标页岩气层的地层有效应力之前，还包括：
获取多个岩心样品的取样深度，所述多个岩心样品是在不同深度的页岩气层进行取样得到的，且所述多个岩心样品的取样深度均不同；
基于所述多个岩心样品的取样深度，确定每个岩心样品的上覆岩层压力，得到多个上覆岩层压力；
从零压力到最大上覆岩层压力的一半之间的范围内选择多个有效应力；
基于所述多个有效应力获取每个岩心样品对应的多个孔隙度；
对所述多个岩心样品的孔隙度进行归一化处理，以及对每个岩心样品的多个归一化孔隙度进行均值处理，得到与所述多个有效应力一一对应的多个平均孔隙度；
对所述多个有效应力和所述多个平均孔隙度进行指数曲线拟合，得到所述第二转换公式。
9.一种页岩气层的筛选装置，其特征在于，所述装置包括：
第一确定模块，用于根据目标页岩气层对应的测井资料，确定所述目标页岩气层的修正孔隙度，所述修正孔隙度是指所述目标页岩气层中除有机孔之外的其他孔的孔隙度；
第二确定模块，用于基于所述目标页岩气层的修正孔隙度，确定所述目标页岩气层的地层孔隙压力系数；
第三确定模块，用于当所述地层孔隙压力系数大于参考数值时，将所述目标页岩气层确定为待开发的页岩气层。
10.一种计算机可读存储介质，其特征在于，所述存储介质内存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1-8任一所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110484937.0A CN115263284A (zh) | 2021-04-30 | 2021-04-30 | 页岩气层的筛选方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110484937.0A CN115263284A (zh) | 2021-04-30 | 2021-04-30 | 页岩气层的筛选方法、装置及存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115263284A true CN115263284A (zh) | 2022-11-01 |
Family
ID=83745462
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110484937.0A Pending CN115263284A (zh) | 2021-04-30 | 2021-04-30 | 页岩气层的筛选方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115263284A (zh) |
-
2021
- 2021-04-30 CN CN202110484937.0A patent/CN115263284A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111396013B (zh) | 页岩气井压裂改造方案的确定方法、装置及存储介质 | |
| CN109191373A (zh) | 导航坐标转换高斯投影坐标的方法、系统、设备及介质 | |
| CN111485850B (zh) | 基于大孔道数据的油井堵水方法及装置 | |
| CN111255434B (zh) | 气井的试井方法、装置和计算机存储介质 | |
| CN109002420A (zh) | 影响程度的确定方法、装置及存储介质 | |
| CN110109178B (zh) | 优势砂体分布区域的确定方法、装置及存储介质 | |
| CN110083851B (zh) | 气井井底压力的确定方法、装置及存储介质 | |
| CN115263284A (zh) | 页岩气层的筛选方法、装置及存储介质 | |
| CN113216926A (zh) | 页岩气井的筛选方法、装置、终端及存储介质 | |
| CN111999140B (zh) | 制备声发射样品、确定压裂泵的方法、装置及存储介质 | |
| CN113495046A (zh) | 储层类型的确定方法、装置及存储介质 | |
| CN111188608B (zh) | 页岩气井压裂改造方案的确定方法、装置及存储介质 | |
| CN111852453B (zh) | 渗流通道的预警方法、装置及存储介质 | |
| CN110824566B (zh) | 地震属性融合方法、装置及存储介质 | |
| CN115217466A (zh) | 驱油体系的乳化参数确定方法、装置、设备及存储介质 | |
| CN110894782B (zh) | 确定油藏存气量的方法、装置和存储介质 | |
| CN108487902B (zh) | 井底流压的确定方法、装置及存储介质 | |
| CN111426813B (zh) | 砂体连通质量的确定方法、装置及存储介质 | |
| CN114864008A (zh) | 确定储层中tsr的反应强度的方法和装置 | |
| CN113027426B (zh) | 漏失压力确定方法、装置及存储介质 | |
| CN114458306A (zh) | 基于噪声测井的流体流量的确定方法、装置、设备及介质 | |
| CN111425189B (zh) | 超深层缝洞碳酸盐岩储层品质评价方法、装置及存储介质 | |
| CN111812744B (zh) | 地层可压性确定方法、装置及计算机存储介质 | |
| CN110320561B (zh) | 火成岩储层边界的刻画方法、装置及存储介质 | |
| CN114165219B (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 |