CN109025961B - 页岩储层含气饱和度的计算方法、装置和电子设备 - Google Patents
页岩储层含气饱和度的计算方法、装置和电子设备 Download PDFInfo
- Publication number
- CN109025961B CN109025961B CN201810841567.XA CN201810841567A CN109025961B CN 109025961 B CN109025961 B CN 109025961B CN 201810841567 A CN201810841567 A CN 201810841567A CN 109025961 B CN109025961 B CN 109025961B
- Authority
- CN
- China
- Prior art keywords
- neutron
- shale reservoir
- density
- curve
- logging
- 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.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 68
- 230000004044 response Effects 0.000 claims abstract description 61
- 238000000034 method Methods 0.000 claims abstract description 51
- 239000002734 clay mineral Substances 0.000 claims abstract description 46
- 229910052770 Uranium Inorganic materials 0.000 claims abstract description 15
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 claims abstract description 15
- 238000010219 correlation analysis Methods 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 15
- 239000004927 clay Substances 0.000 claims description 11
- 239000011148 porous material Substances 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 9
- 239000012530 fluid Substances 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 238000012216 screening Methods 0.000 claims description 6
- 239000011435 rock Substances 0.000 description 13
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000002474 experimental method Methods 0.000 description 5
- 229930195733 hydrocarbon Natural products 0.000 description 5
- 150000002430 hydrocarbons Chemical class 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000004215 Carbon black (E152) Substances 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 239000005416 organic matter Substances 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 229910052683 pyrite Inorganic materials 0.000 description 2
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 description 2
- 239000011028 pyrite Substances 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000000802 evaporation-induced self-assembly Methods 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000011425 standardization method Methods 0.000 description 1
- 238000006467 substitution reaction 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
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
- E21B49/087—Well testing, e.g. testing for reservoir productivity or formation parameters
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种页岩储层含气饱和度的计算方法、装置和电子设备，涉及页岩气勘查技术领域，该方法包括选取参与目标页岩储层含气饱和度计算的的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；标准化处理原始测井曲线得到标准化测井曲线；根据该标准化测井曲线计算粘土矿物含量，以及获取区域致密页岩的密度测井响应值和中子测井响应值，并计算目标页岩储层的密度视孔隙度和中子视孔隙度；根据计算的粘土矿物含量、密度视孔隙度和中子视孔隙度计算目标页岩储层的含气饱和度。本发明实施例提供的一种页岩储层含气饱和度的计算方法、装置和电子设备，可以减小页岩储层含气饱和度的计算误差，增加计算结果的可靠性。
Description
技术领域
本发明涉及页岩气勘查技术领域，尤其是涉及一种页岩储层含气饱和度的计算方法、装置和电子设备。
背景技术
传统的含气饱和度计算模型主要以Archie公式为主。但该方法尚存在以下几点不足：首先，Archie公式适用于孔隙结构简单的中高孔隙度、高渗透率的纯砂岩地层，而页岩储层岩性较复杂，大多为低孔隙度、特低渗透率地层；其次，Archie公式中的重要参数包括a(与岩性有关的系数)、b(与岩性有关的系数)、m(胶结指数)、n(饱和度指数)，这四个参数通常根据岩电实验确定，但页岩储层样品很难满足岩电实验要求，故而影响四个参数的求取，对饱和度计算的准确性产生直接影响；另外，Archie公式中采用了地层电阻率，前提是认为地层中仅有地层水导电，但页岩储层层理发育、且普遍发育黄铁矿，层间缝和黄铁矿使电阻率测井响应值降低，影响含气饱和度解释结果的可靠性。
目前，国内外大多采用Total-shale公式计算页岩储层含气饱和度。Total-shale公式实质是在Archie公式基础上考虑了泥质的影响，但含气饱和度的计算结果仍然与岩电实验参数、地层电阻率数值密切相关，仍不能避开Archie公式在页岩储层含气饱和度计算中的缺陷。
发明内容
有鉴于此，本发明的目的在于提供一种页岩储层含气饱和度的计算方法、装置和电子设备，可以减小页岩储层含气饱和度的计算误差，增加计算结果的可靠性。
第一方面，本发明实施例提供了一种页岩储层含气饱和度的计算方法，包括：选取参与目标页岩储层含气饱和度计算的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；标准化处理该原始测井曲线得到标准化测井曲线；根据该标准化测井曲线计算目标页岩储层的粘土矿物含量；获取区域致密页岩的密度测井响应值和中子测井响应值，并根据该密度测井响应值和密度曲线、该中子测井响应值和中子曲线计算目标页岩储层的密度视孔隙度和中子视孔隙度；根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算目标页岩储层的含气饱和度。
结合第一方面，本发明实施例提供了第一方面的第一种可能的实施方式，其中，在上述标准化处理该原始测井曲线得到标准化测井曲线的步骤之前，还包括：获取区域页岩储层岩心的含水饱和度数据；对该含水饱和度数据和原始测井曲线进行相关性分析得到相关性分析结果；根据该相关性分析结果从该原始测井曲线中筛选参与计算含气饱和度的测井曲线。
结合第一方面的第一种可能的实施方式，本发明实施例提供了第一方面的第二种可能的实施方式，其中，上述根据该相关性分析结果从该原始测井曲线中筛选参与计算含气饱和度的测井曲线的步骤，包括：根据该相关性分析结果，从该原始测井曲线中选取相关性好的测井曲线参与计算含气饱和度。
结合第一方面，本发明实施例提供了第一方面的第三种可能的实施方式，其中，上述根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算目标页岩储层的含气饱和度的步骤，包括：根据该密度视孔隙度和该中子视孔隙度计算密度-中子视孔隙度差；根据该密度-中子视孔隙度差和该粘土矿物含量计算目标页岩储层的含气饱和度。
结合第一方面的第三种可能的实施方式，本发明实施例提供了第一方面的第四种可能的实施方式，其中，上述根据该密度-中子视孔隙度差和该粘土矿物含量计算目标页岩储层的含气饱和度的公式为：Sg＝1-(a×Vclay+b×ΦD-N+c)，ΦD-N＝ΦD-ΦN，式中，Sg为页岩储层含气饱和度，Vclay为页岩储层的粘土矿物含量，ΦD为密度视孔隙度，ΦN中子视孔隙度，ΦD-N为密度-中子视孔隙度差，a、b、c为区域经验回归系数。
结合第一方面，本发明实施例提供了第一方面的第五种可能的实施方式，其中，上述根据该密度测井响应值和该中子测井响应值计算目标页岩储层的密度视孔隙度和中子视孔隙度的计算公式为：
式中，ΦD为密度视孔隙度，ΦN中子视孔隙度，ρb为页岩储层段密度测井响应值，ρma为区域致密页岩的密度测井响应值，ρf为地层孔隙流体密度值，CNL为页岩储层段中子测井响应值，CNLma为区域致密页岩的中子测井响应值，CNLf为地层孔隙流体中子值。
结合第一方面，本发明实施例提供了第一方面的第六种可能的实施方式，其中，上述根据该标准化测井曲线计算该目标页岩储层的粘土矿物含量的步骤，包括：获取区域页岩储层岩心的粘土矿物含量测试数据；将该标准化测井曲线与该粘土矿物含量测试数据进行拟合，得到拟合公式；根据该拟合公式计算该目标页岩储层的粘土矿物含量。
第二方面，本发明实施例还提供了一种页岩储层含气饱和度的计算装置，该装置包括：原始测井曲线选取模块，用于选取参与目标页岩储层含气饱和度计算的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；标准化处理模块，用于标准化处理该原始测井曲线得到标准化测井曲线；粘土矿物含量计算模块，用于根据该标准化测井曲线计算该目标页岩储层的粘土矿物含量；视孔隙度计算模块，用于获取区域致密页岩的密度测井响应值和中子测井响应值，并根据该密度测井响应值和密度曲线、该中子测井响应值和中子曲线计算该目标页岩储层的密度视孔隙度和中子视孔隙度；含气饱和度计算模块，用于根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算该目标页岩储层的含气饱和度。
第三方面，本发明实施例还提供了一种页岩储层含气饱和度的计算装置，该装置包括处理器，存储器，总线和通信接口，该处理器、通信接口和存储器通过该总线连接；该存储器用于存储程序；该处理器，用于通过该总线调用存储在该存储器中的程序，执行上述第一方面及其可能的实施方式之一提供的页岩储层含气饱和度的计算方法。
第四方面，本发明实施例还提供了一种电子设备，包括存储器、处理器，该存储器中存储有可在该处理器上运行的计算机程序，该处理器执行该计算机程序时实现上述第一方面及其可能的实施方式之一提供的页岩储层含气饱和度的计算方法的步骤。
本发明实施例带来了以下有益效果：
本发明实施例提供的一种页岩储层含气饱和度的计算方法、装置和电子设备，该方法包括选取参与目标页岩储层含气饱和度计算的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；标准化处理该原始测井曲线得到标准化测井曲线；根据该标准化测井曲线计算目标页岩储层的粘土矿物含量；获取区域致密页岩的密度测井响应值和中子测井响应值，并根据该密度测井响应值和密度曲线、该中子测井响应值和中子曲线计算目标页岩储层的密度视孔隙度和中子视孔隙度；根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算目标页岩储层的含气饱和度；可以减小页岩储层含气饱和度的计算误差，增加计算结果的可靠性。
本公开的其他特征和优点将在随后的说明书中阐述，或者，部分特征和优点可以从说明书推知或毫无疑义地确定，或者通过实施本公开的上述技术即可得知。
为使本公开的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种页岩储层含气饱和度的计算方法的流程图；
图2为本发明实施例提供的另一种页岩储层含气饱和度的计算方法的流程图；
图3为本发明实施例提供的一种页岩储层含气饱和度计算的多种方法处理结果对比图；
图4为本发明实施例提供的一种页岩储层含气饱和度的计算装置的结构示意图；
图5为本发明实施例提供的另一种页岩储层含气饱和度的计算装置的结构示意图。
图标：
41-原始测井曲线选取模块；42-标准化处理模块；43-粘土矿物含量计算模块；44-视孔隙度计算模块；45-含气饱和度计算模块；50-处理器；51-存储器；52-总线；53-通信接口。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
目前，计算页岩气饱和度的方法包括使用Archie公式和使用Total-shale公式，但是，这两种方法都存在因计算过程受岩电实验参数、地层电阻率数值影响而计算结果可靠性差的问题。基于此，本发明实施例提供的一种页岩储层含气饱和度的计算方法、装置和电子设备，可以减小页岩储层含气饱和度的计算误差，增加计算结果的可靠性。
为便于对本实施例进行理解，首先对本发明实施例所公开的一种页岩储层含气饱和度的计算方法进行详细介绍。
实施例一
本发明实施例提供了一种页岩储层含气饱和度的计算方法，如图1所示，为该方法的流程示意图，由图1可见，该方法包括以下步骤：
步骤S102：选取参与目标页岩储层含气饱和度计算的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线。
原始测井曲线是实际测井之后得到的未经处理的测井曲线，其中，测井的方法有很多，电、声、放射性是三种基本测井方法，特殊测井方法有电缆地层测试、地层倾角测井、成像测井、核磁共振测井等，其他测井方式还有随钻测井等。各种测井方法基本上是间接地、有条件地反映岩层地质特性的某一侧面，要全面认识地下地质面貌，发现和评价油气层，需要综合使用多种测井方法。
目标页岩储层是经分析后认为有利含气的页岩储层。在确定页岩储层是否为有利含气储层时，首先，需要分析测井曲线的响应特征，结合现场录井岩性描述、气测全烃和甲烷显示资料，确定含有机质的页岩层段；其次，在含有机质的页岩层段中，根据自然伽马能谱、声波、电阻率、密度、中子等资料，划分出高自然伽马、高铀、高声波时差、中等电阻率、低密度、相对低中子的有利含气页岩储层。
在本实施例中，获取的原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线。
步骤S104：标准化处理该原始测井曲线得到标准化测井曲线。
这里，根据区域标志层的测井曲线响应值的统计分析结果，对参与含气饱和度计算的原始测井曲线进行标准化处理。测井曲线标准化是为测井解释储集层参数，进行油藏描述的前期准备工作。
步骤S106：根据该标准化测井曲线计算目标页岩储层的粘土矿物含量。
其中，计算目标页岩储层的粘土矿物含量的步骤包括：
(a1)获取区域页岩气田岩心的粘土矿物含量测试数据；
(a2)将该标准化测井曲线与该粘土矿物含量测试数据进行拟合，得到拟合公式；
(a3)根据该拟合公式计算该目标页岩储层的粘土矿物含量。
步骤S108：获取区域致密页岩的密度测井响应值和中子测井响应值，并根据该密度测井响应值和密度曲线、该中子测井响应值和中子曲线计算目标页岩储层的密度视孔隙度和中子视孔隙度。
根据区域致密页岩的测井曲线响应值的统计分析结果，获取区域致密页岩的密度测井响应值和中子测井响应值。
根据标准化密度曲线和中子曲线，结合得到的区域致密页岩的密度测井响应值和中子测井响应值，计算目标页岩储层的密度视孔隙度和中子视孔隙度，其计算公式如下：
式中，ΦD为密度视孔隙度，ΦN中子视孔隙度，ρb为页岩储层段密度测井响应值，ρma为区域致密页岩的密度测井响应值，ρf为地层孔隙流体密度值，CNL为页岩储层段中子测井响应值，CNLma为区域致密页岩的中子测井响应值，CNLf为地层孔隙流体中子值。
步骤S110：根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算目标页岩储层的含气饱和度。
首先，根据该密度视孔隙度和该中子视孔隙度计算密度-中子视孔隙度差；然后，根据该密度-中子视孔隙度差和该粘土矿物含量计算目标页岩储层的含气饱和度。其中，目标页岩储层的含气饱和度的计算公式为：
Sg＝1-(a×Vclay+b×ΦD-N+c)
ΦD-N＝ΦD-ΦN
式中，Sg为页岩储层含气饱和度，Vclay为页岩储层的粘土矿物含量，ΦD为密度视孔隙度，ΦN中子视孔隙度，ΦD-N为密度-中子视孔隙度差，a、b、c为区域经验回归系数。
这样，本发明实施例，通过计算页岩储层中的粘土矿物含量以及密度视孔隙度和中子视孔隙度，从而求得页岩储层的含气饱和度，避免了传统的页岩储层的含气饱和度计算方法受岩电实验参数、地层电阻率数值测量不准确的影响，该计算方法的整个过程计算误差更小，计算结果更加可靠。
本发明实施例提供的一种页岩储层含气饱和度的计算方法，该方法包括选取参与目标页岩储层含气饱和度计算的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；标准化处理该原始测井曲线得到标准化测井曲线；根据该标准化测井曲线计算目标页岩储层的粘土矿物含量；获取区域致密页岩的密度测井响应值和中子测井响应值，并根据该密度测井响应值和密度曲线、该中子测井响应值和中子曲线计算目标页岩储层的密度视孔隙度和中子视孔隙度；根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算目标页岩储层的含气饱和度；可以减小页岩储层含气饱和度的计算误差，增加计算结果的可靠性。
实施例二
在图1的基础上，参见图2，为本发明实施例提供的另一种页岩储层含气饱和度的计算方法的流程图，图2示出的实施方式中，该方法在上述标准化处理原始测井曲线得到标准化测井曲线的步骤之前，还包括以下步骤：
步骤S202：获取区域页岩储层岩心的含水饱和度数据。
含水饱和度是指在油层中水所占的孔隙的体积与岩石孔隙体积之比。流体饱和度是最重要的储层参数之一。
步骤S204：对该含水饱和度数据和原始测井曲线进行相关性分析得到相关性分析结果。
这里，相关性分析是指对两个或多个具备相关性的变量元素进行分析，从而衡量两个变量因素的相关密切程度。
步骤S206：根据该相关性分析结果从该原始测井曲线中筛选参与计算含气饱和度的测井曲线。
根据得到的相关性分析结果，从上述原始测井曲线中选取相关性好的测井曲线参与计算含气饱和度。
这样，本实施例在上述实施例一的基础上，在对获得的原始测井曲线进行标准化处理之前，还进行了筛选的工作，通过含水饱和度数据和原始测井曲线的相关性分析，筛选出相关性好的测井曲线进行实际的含气饱和度计算，对于最终计算得到的含气饱和度数据，提高了其准确度和可信度。
实施例三
本发明实施例介绍了一个计算页岩储层含气饱和度的应用实例，在该实例中，分别利用本发明的页岩储层含气饱和度的计算方法、传统的Archie公式以及Total-shale公式对同一实验场地的测井资料进行了处理，其中，具体介绍了利用本发明的页岩储层含气饱和度的计算方法的计算过程，内容如下。
首先，分析待计算井的各条测井曲线响应特征，寻找现场录井岩性描述为页岩的层段，分析现场气测资料的全烃和甲烷异常显示情况，综合确定含有机质的页岩层段。在本实施例中，3559.6m-3649.2m为含有机质的页岩层段。
其次，在已划分为含有机质的页岩层段中，根据高自然伽马、相对低去铀伽马、中等电阻率(本实施例中为10～90Ωm)、高声波时差、低密度、相对低中子的特征，划分有利含气页岩储层。在本实施例中，3559.6m-3649.2m含有机质的页岩层段同时也为有利含气页岩储层段，并综合各类资料自上而下细分为Ⅲ类页岩气层、Ⅱ类页岩气层、Ⅰ类页岩气层，参见图3所示。
然后，选取参与含气饱和度计算的原始测井曲线，本实施例中，根据涪陵页岩气田页岩储层岩心实测含水饱和度和测井曲线的相关性分析，选取去铀伽马(KTH)、密度(DEN)、中子(CNL)为计算含气饱和度的原始测井曲线，采用逐点采样方式(一般为8点/1m)。
接着，根据区域标志层和标准化方法，对KTH、DEN、CNL测井值进行标准化处理，并输出标准化后的KTH、DEN、CNL曲线。
并且，计算页岩气储层段的粘土矿物含量Vclay，本实施例采用的是经涪陵页岩气田岩心粘土矿物含量测试数据与测井曲线进行拟合的公式：
Vclay＝0.157×KTH+34.51×DEN+1.924×CNL-95.773
上式中，Vclay的单位为百分数。
选取区域致密页岩的密度测井响应值ρma、中子测井响应值CNLma，其中，ρma＝2.7g/cm3、CNLma＝16％。
然后，计算页岩气储层段的密度、中子视孔隙度，采用经验公式：
密度视孔隙度：
中子视孔隙度：
式中，ρb为页岩储层段标准化后的密度值，ρf为地层孔隙流体密度值、取值1.0g/cm3；CNL为页岩储层段标准化后的中子值，CNLf为地层孔隙流体中子值，取值100％。
计算密度—中子视孔隙度差ΦD-N＝ΦD-ΦN，在本实施例中，ΦD、ΦN都转化为百分数形式，ΦD-N也为百分数。
计算页岩储层含气饱和度，采用公式：
Sg＝1-(a×Vclay+b×ΦD-N+c)
在本实施例中，
Sg＝100-(0.236×Vclay-0.938×ΦD-N+38.766)
上式中，Sg、Vclay、ΦD-N都采用百分数形式。
最后，输出含气饱和度计算结果。
如图3所示，为本发明实施例提供的一种页岩储层含气饱和度计算的多种方法处理结果对比图，在图3中，
第一道为自然伽马、去铀伽马曲线，自然伽马为实线、去铀伽马为虚线；
第二道为深度道；
第三道为录井岩性剖面；
第四道为气测全烃、甲烷显示值，全烃为实线、甲烷为虚线；
第五道为标准化后的电阻率曲线；
第六道为标准化后的声波时差曲线；
第七道为标准化后的密度、中子曲线，密度为实线、中子为虚线；
第八道为页岩储层划分道；
第九道为岩心实测含气饱和度、Archie公式计算含气饱和度的对比显示，岩心实测含气饱和度为散点，Archie公式计算的含气饱和度为实线；
第十道为岩心实测含气饱和度、Total-shale公式计算含气饱和度的对比显示，岩心实测含气饱和度为散点，Total-shale公式计算的含气饱和度为实线；
第十一道为岩心实测含气饱和度、本发明计算含气饱和度的对比显示，岩心实测含气饱和度为散点，本发明计算的含气饱和度为实线。
通过本发明提供的页岩储层含气饱和度计算方法，得到本实例中3559.6-3610.6m井段含气饱和度平均为54.6％，解释为Ⅲ类页岩气层；3610.6-3630.8m井段含气饱和度平均为60.2％，解释为Ⅱ类页岩气层；3630.8-3649.2m井段含气饱和度平均为66.1％，解释为Ⅰ类页岩气层。
由图3可见，与传统含气饱和度算法相比，第九道中Archie公式计算含气饱和度与岩心实测含气饱和度吻合差、误差大；第十道中Total-shale公式计算含气饱和度与岩心实测含气饱和度吻合性比Archie算法略好，但误差仍较大；第十一道中本发明计算含气饱和度与岩心实测含气饱和度吻合好、误差较小，说明本发明计算的含气饱和度更符合页岩储层实际情况。
实施例四
如图4所示，为本发明实施例提供的一种页岩储层含气饱和度的计算装置的结构示意图，由图4可见，该装置包括原始测井曲线选取模块41、标准化处理模块42、粘土矿物含量计算模块43、视孔隙度计算模块44和含气饱和度计算模块45，其中，各个模块的功能如下：
原始测井曲线选取模块41，用于选取参与目标页岩储层含气饱和度计算的原始测井曲线；该原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；
标准化处理模块42，用于标准化处理该原始测井曲线得到标准化测井曲线；
粘土矿物含量计算模块43，用于根据该标准化测井曲线计算该目标页岩储层的粘土矿物含量；
视孔隙度计算模块44，用于获取区域致密页岩的密度测井响应值和中子测井响应值，并根据该密度测井响应值和密度曲线、该中子测井响应值和中子曲线计算该目标页岩储层的密度视孔隙度和中子视孔隙度；
含气饱和度计算模块45，用于根据该粘土矿物含量、密度视孔隙度和中子视孔隙度计算该目标页岩储层的含气饱和度。
本发明实施例所提供的页岩储层含气饱和度的计算装置，其实现原理及产生的技术效果和前述页岩储层含气饱和度的计算方法实施例相同，为简要描述，装置实施例部分未提及之处，可参考前述方法实施例中相应内容。
实施例五
参见图5，本发明实施例还提供一种页岩储层含气饱和度的计算装置，包括：处理器50，存储器51，总线52和通信接口53，处理器50、通信接口53和存储器51通过总线52连接；处理器50用于执行存储器51中存储的可执行模块，例如计算机程序。
其中，存储器51可能包含高速随机存取存储器(RAM，Random Access Memory)，也可能还包括非易失性存储器(non-volatile memory)，例如至少一个磁盘存储器。通过至少一个通信接口53(可以是有线或者无线)实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。
总线52可以是ISA总线、PCI总线或EISA总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，图5中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
其中，存储器51用于存储程序，处理器50在接收到执行指令后，执行程序，前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器50中，或者由处理器50实现。
处理器50可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器50中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器50可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(Digital SignalProcessing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器51，处理器50读取存储器51中的信息，结合其硬件完成上述方法的步骤。
本发明实施例提供的页岩储层含气饱和度的计算装置，与上述实施例一、实施例二及实施例三提供的页岩储层含气饱和度的计算方法具有相同的技术特征，所以也能解决相同的技术问题，达到相同的技术效果。
实施例六
本发明实施例还提供了一种电子设备，包括存储器、处理器，该存储器中存储有可在该处理器上运行的计算机程序，该处理器执行该计算机程序时实现上述实施例一、实施例二、实施例三及其可能的实施方式之一提供的页岩储层含气饱和度的计算方法的步骤。
除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本发明的范围。
在这里示出和描述的所有示例中，任何具体值应被解释为仅仅是示例性的，而不是作为限制，因此，示例性实施例的其他示例可以具有不同的值。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种页岩储层含气饱和度的计算方法，其特征在于，包括：
选取参与目标页岩储层含气饱和度计算的原始测井曲线；所述原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；
标准化处理所述原始测井曲线得到标准化测井曲线；
根据所述标准化测井曲线计算所述目标页岩储层的粘土矿物含量；
获取区域致密页岩的密度测井响应值和中子测井响应值，并根据所述密度测井响应值和所述密度曲线、所述中子测井响应值和所述中子曲线计算所述目标页岩储层的密度视孔隙度和中子视孔隙度；
根据所述密度视孔隙度和所述中子视孔隙度计算密度-中子视孔隙度差；
根据所述密度-中子视孔隙度差和所述粘土矿物含量计算所述目标页岩储层的含气饱和度；
所述根据所述密度-中子视孔隙度差和所述粘土矿物含量计算所述目标页岩储层的含气饱和度的公式为：
Sg＝1-(a×Vclay+b×ΦD-N+c)
ΦD-N＝ΦD-ΦN
式中，Sg为页岩储层含气饱和度，Vclay为页岩储层的粘土矿物含量，ΦD为密度视孔隙度，ΦN中子视孔隙度，ΦD-N为密度-中子视孔隙度差，a、b、c为区域经验回归系数。
2.根据权利要求1所述的页岩储层含气饱和度的计算方法，其特征在于，在所述标准化处理所述原始测井曲线得到标准化测井曲线的步骤之前，还包括：
获取区域页岩储层岩心的含水饱和度数据；
对所述含水饱和度数据和所有原始测井曲线进行相关性分析得到相关性分析结果；
根据所述相关性分析结果从所述原始测井曲线中筛选参与计算含气饱和度的测井曲线。
3.根据权利要求2所述的页岩储层含气饱和度的计算方法，其特征在于，所述根据所述相关性分析结果从所述原始测井曲线中筛选参与计算含气饱和度的测井曲线的步骤，包括：
根据所述相关性分析结果，从所述原始测井曲线中选取相关性好的测井曲线参与计算含气饱和度。
4.根据权利要求1所述的页岩储层含气饱和度的计算方法，其特征在于，所述根据根据所述密度测井响应值和所述密度曲线、所述中子测井响应值和所述中子曲线计算所述目标页岩储层的密度视孔隙度和中子视孔隙度的计算公式为：
式中，ΦD为密度视孔隙度，ΦN中子视孔隙度，ρb为页岩储层段密度测井响应值，ρma为区域致密页岩的密度测井响应值，ρf为地层孔隙流体密度值，CNL为页岩储层段中子测井响应值，CNLma为区域致密页岩的中子测井响应值，CNLf为地层孔隙流体中子值。
5.根据权利要求1所述的页岩储层含气饱和度的计算方法，其特征在于，所述根据所述标准化测井曲线计算所述目标页岩储层的粘土矿物含量的步骤，包括：
获取区域页岩储层岩心的粘土矿物含量测试数据；
将所述标准化测井曲线与所述粘土矿物含量测试数据进行拟合，得到拟合公式；
根据所述拟合公式计算所述目标页岩储层的粘土矿物含量。
6.一种页岩储层含气饱和度的计算装置，其特征在于，包括：
原始测井曲线选取模块，用于选取参与目标页岩储层含气饱和度计算的原始测井曲线；所述原始测井曲线包括去铀伽马曲线、密度曲线和中子曲线；
标准化处理模块，用于标准化处理所述原始测井曲线得到标准化测井曲线；
粘土矿物含量计算模块，用于根据所述标准化测井曲线计算所述目标页岩储层的粘土矿物含量；
视孔隙度计算模块，用于获取区域致密页岩的密度测井响应值和中子测井响应值，并根据所述密度测井响应值和所述密度曲线、所述中子测井响应值和所述中子曲线计算所述目标页岩储层的密度视孔隙度和中子视孔隙度；
含气饱和度计算模块，用于根据所述粘土矿物含量、所述密度视孔隙度和所述中子视孔隙度计算所述目标页岩储层的含气饱和度。
7.一种页岩储层含气饱和度的计算装置，其特征在于，所述装置包括处理器，存储器，总线和通信接口，所述处理器、通信接口和存储器通过所述总线连接；
所述存储器用于存储程序；
所述处理器，用于通过所述总线调用存储在所述存储器中的程序，执行所述权利要求1-5任一所述方法。
8.一种电子设备，包括存储器、处理器，所述存储器中存储有可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至5任一项所述的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810841567.XA CN109025961B (zh) | 2018-07-26 | 2018-07-26 | 页岩储层含气饱和度的计算方法、装置和电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810841567.XA CN109025961B (zh) | 2018-07-26 | 2018-07-26 | 页岩储层含气饱和度的计算方法、装置和电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109025961A CN109025961A (zh) | 2018-12-18 |
| CN109025961B true CN109025961B (zh) | 2020-03-06 |
Family
ID=64647030
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810841567.XA Active CN109025961B (zh) | 2018-07-26 | 2018-07-26 | 页岩储层含气饱和度的计算方法、装置和电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109025961B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112922587B (zh) * | 2019-11-20 | 2023-10-27 | 中国石油天然气股份有限公司 | 盐霜位置的识别方法及装置 |
| CN112987121A (zh) * | 2019-12-13 | 2021-06-18 | 北京国双科技有限公司 | 一种孔隙度预测方法、装置、电子设备及存储介质 |
| CN113139877B (zh) * | 2020-01-17 | 2024-03-26 | 中国石油天然气股份有限公司 | 页岩气井的筛选方法及装置 |
| CN111562207B (zh) * | 2020-05-13 | 2022-06-14 | 南京宏创地质勘查技术服务有限公司 | 一种页岩含油饱和度的计算方法 |
| CN114135269B (zh) * | 2020-08-12 | 2024-04-12 | 中国石油化工股份有限公司 | 一种致密砂岩油层识别方法及装置 |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015021088A1 (en) * | 2013-08-06 | 2015-02-12 | Schlumberger Canada Limited | Methods for determining a saturation-height function in oil and gas reservoirs |
| CN103912268B (zh) * | 2014-03-28 | 2016-09-28 | 中石化江汉石油工程有限公司测录井公司 | 基于toc的页岩储层含气饱和度确定方法 |
| CN104500049A (zh) * | 2014-10-20 | 2015-04-08 | 成都创源油气技术开发有限公司 | 页岩气地球物理快速评价方法 |
| CN104573344B (zh) * | 2014-12-25 | 2017-10-13 | 中国海洋石油总公司 | 一种通过测井数据获取页岩储层含气量的方法 |
| CN104712329B (zh) * | 2015-01-29 | 2018-06-05 | 中国石油大学(华东) | 一种泥页岩油气饱和度的计算模型 |
| US10422917B2 (en) * | 2016-04-19 | 2019-09-24 | Halliburton Energy Services, Inc. | Sand-resistivity/conductivity and saturation evaluation in laminated formations with biaxial anisotropy |
| CN207261006U (zh) * | 2017-08-23 | 2018-04-20 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 致密油气藏油水饱和度测试仪 |
| CN107808032B (zh) * | 2017-09-28 | 2021-07-16 | 中国石油化工股份有限公司 | 页岩气储层含水饱和度计算方法及系统 |
| CN107808055B (zh) * | 2017-10-30 | 2021-05-14 | 中石化石油工程技术服务有限公司 | 一种基于双重影响的页岩气饱和度测井计算方法 |
-
2018
- 2018-07-26 CN CN201810841567.XA patent/CN109025961B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109025961A (zh) | 2018-12-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109025961B (zh) | 页岩储层含气饱和度的计算方法、装置和电子设备 | |
| US10324229B2 (en) | System and method of pore type classification for petrophysical rock typing | |
| US9372162B2 (en) | Characterization of subterranean formation properties derived from quantitative X-Ray CT scans of drill cuttings | |
| CN111077568B (zh) | 致密油气储层流体因子检测油气储层的方法及设备 | |
| CN107402176B (zh) | 裂缝孔隙度的确定方法和装置 | |
| US20150234069A1 (en) | System and Method for Quantifying Vug Porosity | |
| US10422917B2 (en) | Sand-resistivity/conductivity and saturation evaluation in laminated formations with biaxial anisotropy | |
| CN107450108B (zh) | 甜点区的确定方法和装置 | |
| CN106355571B (zh) | 一种白云岩储层质量的确定方法及装置 | |
| CN107829731B (zh) | 一种黏土蚀变的火山岩孔隙度校正方法 | |
| CN112363242B (zh) | 基于测录井融合的储层流体识别方法与装置 | |
| CN107688037A (zh) | 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法 | |
| CN105064986A (zh) | 利用常规测录井资料构建储层四性关系谱的方法 | |
| CN111625750B (zh) | 一种针对致密砂岩储层气层与气水同层的识别方法 | |
| CN110344822B (zh) | 含钙质致密储层含水饱和度的确定方法 | |
| CN114427457A (zh) | 潮坪相碳酸盐岩储层测井五性关系确定、及测井评价方法 | |
| CN106934725B (zh) | 岩石储层中值半径预测模型建立方法、装置及系统 | |
| CN114183121A (zh) | 裂缝有效性定量评价方法、装置、电子设备及存储介质 | |
| AU2017279838B1 (en) | Method for classifying deep rock geofacies based on data mining | |
| CN116821786A (zh) | 一种低渗储层分类方法、装置及介质 | |
| CN112130227B (zh) | 一种地表水侵入型储集层中油水层识别方法 | |
| US20160146970A1 (en) | Automated method for selecting positions within formations from which to extract samples thereof | |
| CN114428365A (zh) | 一种基于含气特征的致密砂岩气层测井识别方法及装置 | |
| CN117328852B (zh) | 海陆过渡相页岩气有利区确定方法、终端及介质 | |
| CN113960666B (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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200414 Address after: 100027 Beijing, Chaoyangmen, North Street, No. 22, No. Co-patentee after: EXPLORATION AND DEVELOPMENT RESEARCH INSTITUTE OF SINOPEC JIANGHAN OILFIELD Co. Patentee after: China Petroleum & Chemical Corp. Address before: 430000 18 University Garden Road, East Lake hi tech Development Zone, Wuhan, Hubei (China Science and Technology Park) Patentee before: EXPLORATION AND DEVELOPMENT RESEARCH INSTITUTE OF SINOPEC JIANGHAN OILFIELD Co. |
|
| TR01 | Transfer of patent right |