CN108594305A - 页岩储层含气性分析方法及装置 - Google Patents
页岩储层含气性分析方法及装置 Download PDFInfo
- Publication number
- CN108594305A CN108594305A CN201810917461.3A CN201810917461A CN108594305A CN 108594305 A CN108594305 A CN 108594305A CN 201810917461 A CN201810917461 A CN 201810917461A CN 108594305 A CN108594305 A CN 108594305A
- Authority
- CN
- China
- Prior art keywords
- data
- shale
- parameter
- curve
- shale reservoir
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 25
- 239000011435 rock Substances 0.000 claims abstract description 48
- 239000011148 porous material Substances 0.000 claims abstract description 34
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 28
- 238000011156 evaluation Methods 0.000 claims abstract description 21
- 238000005056 compaction Methods 0.000 claims abstract description 11
- 238000005516 engineering process Methods 0.000 claims abstract description 10
- 230000008859 change Effects 0.000 claims description 18
- 230000000153 supplemental effect Effects 0.000 claims description 15
- 238000001615 p wave Methods 0.000 claims description 12
- 238000010606 normalization Methods 0.000 claims description 11
- 230000002706 hydrostatic effect Effects 0.000 claims description 10
- 230000001419 dependent effect Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 241000208340 Araliaceae Species 0.000 claims 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 2
- 235000008434 ginseng Nutrition 0.000 claims 2
- 238000012098 association analyses Methods 0.000 claims 1
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 abstract 2
- 229910052799 carbon Inorganic materials 0.000 abstract 2
- 238000011158 quantitative evaluation Methods 0.000 abstract 1
- 230000015654 memory Effects 0.000 description 15
- 238000012545 processing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000002093 peripheral effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 1
- 235000011613 Pinus brutia Nutrition 0.000 description 1
- 241000018646 Pinus brutia Species 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 125000001183 hydrocarbyl group Chemical group 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229910052711 selenium Inorganic materials 0.000 description 1
- 239000011669 selenium Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种页岩储层含气性分析方法及装置，涉及页岩气勘探开发技术领域。该方法在明确页岩含气性定量评价参数的基础上，求取岩心实测散点数据与弹性参数之间的关联关系，优选出页岩的敏感弹性参数；利用敏感弹性参数，建立页岩有机碳含量、孔隙度、含气量的多元回归方程；通过地震叠前反演技术，依据多元回归方程得到页岩层段的有机碳含量、孔隙度和含气量的三维空间分布；依据单井正常压实趋势线和区域伊顿指数，利用地震层速度体，采用伊顿法得到地层孔隙压力数据体；综合页岩含气性定量评价参数，求得页岩含气性综合评价指数。本发明提供的方法及装置可准确预测出页岩含气性综合评价指数，实现对页岩储层的含气性定量化评价。
Description
技术领域
本发明涉及页岩气勘探开发技术领域，具体而言，涉及一种页岩储层含气性分析方法及装置。
背景技术
现有常规储层含气性评价常基于含油气地层对地震波高频能量吸收强的特点，利用叠后地震资料的频谱衰减技术进行含油气平面预测。
然而，采用频谱衰减技术只能定性识别出含油气的大体区域，无法量化评价含气性的平面差异。同时，造成地震波能量和频率衰减的因素是多方面的，除了页岩储层含油气的影响外，还可能是由于页岩储层岩性、岩相变化以及薄层调谐等引起的，而泥页岩本身属于高吸收岩性页岩储层，且纵向非均质性强，频谱衰减属性的预测结果包含多种因素的综合反映，不能完全体现出页岩含气性的差异。另外，频率衰减属性的预测精度受地震资料品质、有效带宽以及主频等因素的制约，预测精度有限。
发明内容
本发明的目的在于提供一种页岩储层含气性分析方法及装置，以改善上述的问题。
本发明是这样实现的：
第一方面，本发明实施例提供了一种页岩储层含气性分析方法，包括：
依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线，所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨式模量曲线、泊松比曲线和纵横波速度比曲线；
依据积差方法求取通过岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数，所述岩心实测散点数据包括有机含碳量、孔隙度和含气量；
选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数；
以敏感弹性参数为自变量、与所述敏感弹性参数对应的岩心实测散点数据为因变量，通过最小二乘法确定与所述敏感弹性参数对应的岩心实测散点数据最佳匹配的多元回归方程；
基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到弹性参数体数据，所述弹性参数数据体包括纵波阻抗三维数据体、横波阻抗三维数据体和密度三维数据体；
依据所述弹性参数数据体得到对应的敏感弹性参数数据体；
依据所述敏感弹性参数数据体和所述多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体；
利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线；
依据所述声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数；
依据所述声波速度趋势线和所述区域伊顿指数，得到地层孔隙压力数据体；
依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体；
采用极值法对对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理；
对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。
可选地，所述岩心实测散点数据与弹性参数曲线之间的关联系数为其中，R为关联系数，n为进行关联分析的样本个数，xi为有机含碳量、孔隙度和含气量中的其中一个，yi为弹性参数曲线中的对应样本值。
可选地，所述多元回归方程为Y＝β0+β1χ1+β2χ2+...+βpχp+ε，其中，Y为有机含碳量、孔隙度和含气量中的其中一个，χ1、χ2、χp为敏感弹性参数，β0、β1、β2、βp为建立回归方程的回归系数，ε为随机误差项，当ε为0时得到最佳匹配的多元回归方程。
可选地，所述声波速度趋势线为其中，VP为泥岩地层的声波速度，VPMatrix为页岩储层骨架的声波速度，VPTop为深度为0时声波的初始速度，b为页岩储层实压系数，Zml为垂直深度。
可选地，所述单井地层孔隙压力为Pp＝Po-(Po-Pw)×(νinst//νn)N，其中，Pp为预测的页岩储层孔隙压力，Po为静岩压力，Pw为静水压力，νn为正常压实页岩储层的声波速度，νinst为页岩储层实测的声波速度，N为幂指数。
可选地，所述地层压力系数体为
可选地，归一化处理的公式为其中，XN为归一化处理后的参数，X为需要进行归一化处理的参数值；Xmin为归一化处理参数的最小值，Xmax为归一化处理参数的最大值。
可选地，所述页岩含气性综合评价指数为C＝(X1，X2，X3，X4)(A1，A2，A3，A4)T，其中，X1为归一化处理后有机含碳量数据体所对应的归一化值，X2为归一化处理后孔隙度数据体所对应的归一化值，X3为归一化处理后含气量数据体所对应的归一化值，X4为归一化处理后地层压力系数体所对应的归一化值，A1为有机含碳量所对应的权重，A2为孔隙度所对应的权重，A3为含气量所对应的权重，A4为地层压力系数体所对应的权重。
可选地，所述选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数，包括：
选取与对应的岩心实测散点数据之间的关联系数大于0.5的弹性参数曲线作为所述敏感弹性参数曲线。
第二方面，本发明实施例提供了一种页岩储层含气性分析装置，应用于用户终端，页岩储层含气性分析装置包括：运算模块、选取模块和归一化模块；
所述运算模块用于依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线，所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨式模量曲线、泊松比曲线和纵横波速度比曲线；
所述运算模块还用于依据积差方法求取通过岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数，所述岩心实测散点数据包括有机含碳量、孔隙度和含气量；
所述选取模块用于选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数；
所述运算模块还用于以敏感弹性参数为自变量、与所述敏感弹性参数对应的岩心实测散点数据为因变量，通过最小二乘法确定与所述敏感弹性参数对应的岩心实测散点数据最佳匹配的多元回归方程；
所述运算模块还用于基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到弹性参数体数据，所述弹性参数数据体包括纵波阻抗三维数据体、横波阻抗三维数据体和密度三维数据体；
所述运算模块还用于依据所述弹性参数数据体得到对应的敏感弹性参数数据体；
所述运算模块还用于依据所述敏感弹性参数数据体和所述多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体；
所述运算模块还用于利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线；
所述运算模块还用于依据所述声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数；
所述运算模块还用于依据所述声波速度趋势线和所述区域伊顿指数，得到页岩储层压力数据体；
所述运算模块还用于依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体；
所述归一化模块用于采用极值法对对对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理，得到含气性评价参数；
所述运算模块还用于对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。
对于现有技术，本发明提供的页岩储层含气性分析方法及装置具有如下的有益效果：
本发明提供的页岩储层含气性分析方法及装置能够根据地震资料并综合多种评价参数对页岩储层含气性进行综合评价分析，能够准确预测出页岩储层含气性的综合评价指数，提升预测的准确性。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本发明较佳实施例提供的用户终端的方框示意图。
图2为本发明较佳实施例提供的岩层页岩储层含气性分析方法的流程图。
图3为本发明较佳实施例提供的页岩储层含气性分析装置的功能模块图。
图标：100-用户终端；110-页岩储层含气性分析装置；120-存储器；130-存储控制器；140-处理器；150-外设接口；160-输入输出单元；170-显示单元。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，或者是该发明产品使用时惯常摆放的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
此外，术语“水平”、“竖直”、“悬垂”等术语并不表示要求部件绝对水平或悬垂，而是可以稍微倾斜。如“水平”仅仅是指其方向相对“竖直”而言更加水平，并不是表示该结构一定要完全水平，而是可以稍微倾斜。
在本发明的描述中，还需要说明的是，除非另有明确的规定和限定，术语“设置”、“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
请参阅图1，是本发明较佳实施例提供的用户终端100的方框示意图，所述用户终端100包括有包括页岩储层含气性分析装置110、存储器120、存储控制器130、处理器140、外设接口150、输入输出单元160、显示单元170。
所述存储器120、存储控制器130、处理器140、外设接口150、输入输出单元160、显示单元170各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述页岩储层含气性分析装置110包括至少一个可以软件或固件(firmware)的形式存储于所述存储器120中或固化在所述用户终端100的操作系统(operating system，OS)中的软件功能模块。所述处理器140用于执行存储器120中存储的可执行模块，例如所述页岩储层含气性分析装置110包括的软件功能模块或计算机程序。
其中，存储器120可以是，但不限于，随机存取存储器(Random Access Memory，RAM)，只读存储器(Read Only Memory，ROM)，可编程只读存储器(Programmable Read-OnlyMemory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。其中，存储器120用于存储程序，所述处理器140在接收到执行指令后，执行所述程序，本发明实施例任一实施例揭示的流过程定义的用户终端100所执行的方法可以应用于处理器140中，或者由处理器140实现。
处理器140可能是一种集成电路芯片，具有信号的处理能力。上述的处理器140可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述外设接口150将各种输入/输入装置耦合至处理器140以及存储器120。在一些实施例中，外设接口150，处理器140以及存储控制器130可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。
输入输出单元160用于提供给用户输入数据实现用户与所述用户终端100的交互。所述输入输出单元160可以是，但不限于，鼠标和键盘等。
显示单元170在所述用户终端100与用户之间提供一个交互界面(例如用户操作界面)或用于显示图像数据给用户参考。在本实施例中，所述显示单元170可以是液晶显示器或触控显示器。若为触控显示器，其可为支持单点和多点触控操作的电容式触控屏或电阻式触控屏等。支持单点和多点触控操作是指触控显示器能感应到来自该触控显示器上一个或多个位置处同时产生的触控操作，并将该感应到的触控操作交由处理器140进行计算和处理。
请参阅图2，是本发明较佳实施例提供的页岩储层含气性分析方法的流程图，下面将对图2所示的流程进行阐述。
步骤S101，依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线。
本发明实施例提供的方法可应用于上述所述用户终端100。在对页岩含气性进行分析时，用户终端100首先根据预先测得的页岩气田的横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线。其中，横波测井曲线包括有声波在页岩储层中的横波速度、声波在页岩储层中的纵波速度以及页岩储层的密度。所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨氏模量曲线、泊松比曲线和纵横波速度比曲线等一系列弹性参数曲线。
步骤S102，求取岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数。
在对页岩含气性进行分析之前，可通过岩心测试得到页岩气田页岩储层的岩心实测散点数据，所述岩心实测散点数据包括页岩储层的有机含碳量、孔隙度和含气量。
在得到页岩储层随深度变化的多个弹性参数曲线后，用户终端100依据积差方法求取通过岩心测试得到的每个岩心实测散点数据与各个弹性参数曲线上对应弹性参数之间的关联系数。所述岩心实测散点数据与弹性参数曲线之间的关联系数为其中，R为关联系数，n为进行关联分析的样本个数，xi为有机含碳量、孔隙度和含气量中的其中一个，yi为弹性参数曲线中的对应样本值(即弹性参数曲线中对应弹性参数)。弹性参数曲线中的对应样本值是指，弹性参数曲线中的样子值所对应的地底深度与对应岩心实测散点数据所对应的地底深度相同。例如，xi为有机含碳量，yi为纵波阻抗曲线，在求取有机含碳量岩心实测散点数据xi与纵波阻抗曲线yi上的样本值之间的关联系数时，若该有机含碳量岩心实测散点数据xi为地底深度100米处的值，则纵波阻抗曲线yi的样本值也为地底深度100米处的取样值。
步骤S103，选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数。
求取岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数后，用户终端100选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数。本发明实施例中所述预设阈值为0.5，可以理解的，在其他的一些实施例中所述预设阈值也可以根据实际情形进行调整。
步骤S104，基于敏感弹性参数，确定岩心实测散点数据最佳匹配的多元回归方程。
得到敏感弹性参数后，用户终端100以每个敏感弹性参数为自变量、与敏感弹性参数对应的岩心实测散点数据为因变量，通过最小二乘法确定与每个敏感弹性参数对应的岩心实测散点数据(页岩的有机含碳量、孔隙度、含气量)最佳匹配的多元回归方程。所述多元回归方程为Y＝β0+β1χ1+β2χ2+...+βpχp+ε，其中，Y为有机含碳量、孔隙度和含气量中的其中一个，χ1、χ2、χp为敏感弹性参数，β0、β1、β2、βp为建立回归方程的回归系数，ε为随机误差项，当ε为0时得到最佳匹配的多元回归方程。
步骤S105，基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到对应的敏感弹性参数数据体。
用户终端100可预先记录有页岩气田三维叠前道集数据，所述三维叠前道集数据是根据预先测得的地震数据运算得到的。与此同时，用户终端100基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到弹性参数体数据，所述弹性参数数据体包括纵波阻抗三维数据体、横波阻抗三维数据体和密度三维数据体。
基于与步骤S102以及步骤S103相同的原理，用户终端100依据所述弹性参数数据体得到对应的敏感弹性参数数据体。
步骤S106，依据敏感弹性参数数据体和多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体。
得到敏感弹性参数数据体后，用户终端100依据得到的敏感弹性参数数据体和步骤S104中的多元回归方程，进行运算，得的页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体。
步骤S107，利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线。
与此同时，用户终端100利用预先测得的泥质含量对测得的页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线。泥质是指颗粒直径小于0.01mm的碎屑物质，泥质含量是指泥质的体积占岩石总体积的比。所述声波速度趋势线为其中，VP为泥岩地层的声波速度，VPMatrix为页岩储层骨架的声波速度，VPTop为深度为0时声波的初始速度，b为页岩储层实压系数，Zml为垂直深度。
步骤S108，依据声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数。
得到声波速度趋势线后，用户终端100依据得到的声波速度趋势线进行运算得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据计算得到区域伊顿指数。所述单井地层孔隙压力为Pp＝Po-(Po-Pw)×(νinst//νn)N，其中，Pp为预测的页岩储层孔隙压力，Po为静岩压力，Pw为静水压力，νn为正常压实页岩储层的声波速度，νinst为页岩储层实测的声波速度，N为幂指数。
步骤S109，依据声波速度趋势线和区域伊顿指数，得到地层孔隙压力数据体。
得到声波速度趋势线和区域伊顿指数后，用户终端100依据声波速度趋势线和区域伊顿指数，利用地震层速度提，通过Eaton算法即可得到地层孔隙压力数据体。
步骤S110，依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体。
然后，用户终端100依据地层孔隙压力数据体和测得的静水压力数据体进行运算，得到地层压力系数体，所述地层压力系数体为
步骤S111，采用极值法对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理。
有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体等参数是影响页岩储层的含气量的主控因素，基于此用户终端100预先定义有归一化处理公式，在得到有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体后，用户终端100采用极值法对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理。归一化处理的公式为其中，XN为归一化处理后的参数，X为需要进行归一化处理的参数值；Xmin为归一化处理参数的最小值，Xmax为归一化处理参数的最大值。
例如，对于有机含碳量数据体VTOC，设定当VTOC<0.5时，X1＝0，当VTOC>6时，X1＝1；则当0.5<VTOC<6时，对于孔隙度数据体φ，设定当φ<1时，X2＝0，当φ>7时，X2＝1；则当1<φ<7时，对于含气量VGAS，设定当VGAS<1时，X3＝0；当VGAS>10时，X3＝1；则当1<VGAS<10时，对于地层压力系数体αp，设定当αp<1时，X4＝0；当αp>10时，X4＝1；则当1<αp<1.6时，其中，X1为归一化处理后有机含碳量数据体所对应的归一化值，X2为归一化处理后孔隙度数据体所对应的归一化值，X3为归一化处理后含气量数据体所对应的归一化值，X4为归一化处理后地层压力系数体所对应的归一化值，X1、X2、X3和X4即为含气性评价参数。
步骤S112，对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。
归一化处理后，用户终端对不同的含气性评价参数赋予不同的权重，采用权重分配法，计算页岩含气性综合评价指数。所述页岩含气性综合评价指数为C＝(X1，X2，X3，X4)(A1，A2，A3，A4)T，其中，A1为有机含碳量所对应的权重，A2为孔隙度所对应的权重，A3为含气量所对应的权重，A4为地层压力系数体所对应的权重。
含气性评价参数权重系数的分配主要考虑两点：一是对页岩含气性变化反映的敏感程度，二是考虑参数获取的可靠程度。本发明实施例中，地层压力系数体作为评价页岩储层保存条件的重要参数，孔隙度数据体作为页岩储集性能的重要评价参数，是反映页岩含气性差异的两个重要因素，权重赋予0.3，含气量数据体本身虽能较好反映页岩含气性变化，但现有页岩储层含气量测量方法存在一定局限性，导致含气量预测可靠性略有降低，权重赋予0.25，页岩的有机含碳量数据体作为页岩品质评价的关键参数，反映沉积时期的气源物质基础，在页岩气田整体分布相对稳定，权重给予0.15。
请参阅图3，是本发明较佳实施例提供的页岩储层含气性分析装置110的功能模块示意图，所述页岩储层含气性分析装置110包括有运算模块111、选取模块112和归一化模块113。
所述运算模块111用于依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线，所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨式模量曲线、泊松比曲线和纵横波速度比曲线。
可以理解的，所述运算模块111可以用于执行上述的步骤S101。
所述运算模块111还用于依据积差方法求取通过岩心测试岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数，所述岩心实测散点数据包括有机含碳量、孔隙度和含气量。
可以理解的，所述运算模块111还可以用于执行上述的步骤S102。
所述选取模块112用于选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数。
可以理解的，所述选取模块112可以用于执行上述的步骤S103。
所述运算模块111还用于基于敏感弹性参数，确定岩心实测散点数据最佳匹配的多元回归方程。
可以理解的，所述运算模块111还可以用于执行上述的步骤S104。
所述运算模块111还用于基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到对应的敏感弹性参数数据体。
可以理解的，所述运算模块111还可以用于执行上述的步骤S105。
所述运算模块111还用于依据所述敏感弹性参数数据体和所述多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体。
可以理解的，所述运算模块111还可以用于执行上述的步骤S106。
所述运算模块111还用于利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线。
可以理解的，所述运算模块111还可以用于执行上述的步骤S107。
所述运算模块111还用于依据所述声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数。
可以理解的，所述运算模块111还可以用于执行上述的步骤S108。
所述运算模块111还用于依据所述声波速度趋势线和所述区域伊顿指数，得到地层孔隙压力数据体。
可以理解的，所述运算模块111还可以用于执行上述的步骤S109。
所述运算模块111还用于依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体。
可以理解的，所述运算模块111还可以用于执行上述的步骤S110。
所述归一化模块113用于采用极值法对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理。
可以理解的，所述归一化模块113可以用于执行上述的步骤S111。
所述运算模块111还用于对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。
可以理解的，所述运算模块111还可以用于执行上述的步骤S112。
本发明提供的页岩储层含气性分析方法及装置能够根据地震资料并综合多种评价参数对页岩储层含气性进行综合评价分析，能够准确预测出页岩储层含气性的综合评价指数，提升预测的准确性。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种页岩储层含气性分析方法，其特征在于，包括：
依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线，所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨式模量曲线、泊松比曲线和纵横波速度比曲线；
依据积差方法求取通过岩心测试岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数，所述岩心实测散点数据包括有机含碳量、孔隙度和含气量；
选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数；
以敏感弹性参数为自变量、与所述敏感弹性参数对应的岩心实测散点数据为因变量，通过最小二乘法确定与所述敏感弹性参数对应的岩心实测散点数据最佳匹配的多元回归方程；
基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到弹性参数体数据，所述弹性参数数据体包括纵波阻抗三维数据体、横波阻抗三维数据体和密度三维数据体；
依据所述弹性参数数据体得到对应的敏感弹性参数数据体；
依据所述敏感弹性参数数据体和所述多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体；
利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线；
依据所述声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数；
依据所述声波速度趋势线和所述区域伊顿指数，得到地层孔隙压力数据体；
依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体；
采用极值法对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理；
对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。
2.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述岩心实测散点数据与弹性参数曲线之间的关联系数为其中，R为关联系数，n为进行关联分析的样本个数，xi为有机含碳量、孔隙度和含气量中的其中一个，yi为弹性参数曲线中的对应样本值。
3.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述多元回归方程为Y＝β0+β1χ1+β2χ2+...+βpχp+ε，其中，Y为有机含碳量、孔隙度和含气量中的其中一个，χ1、χ2、χp为敏感弹性参数，β0、β1、β2、βp为建立回归方程的回归系数，ε为随机误差项，当ε为0时得到最佳匹配的多元回归方程。
4.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述声波速度趋势线为其中，VP为泥岩地层的声波速度，VPMatrix为页岩储层骨架的声波速度，VPTop为深度为0时声波的初始速度，b为页岩储层实压系数，Zml为垂直深度。
5.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述单井地层孔隙压力为Pp＝Po-(Po-Pw)×(νinst/νn)N，其中，Pp为预测的页岩储层孔隙压力，Po为静岩压力，Pw为静水压力，νn为正常压实页岩储层的声波速度，νinst为页岩储层实测的声波速度，N为幂指数。
6.根据权利要求5所述的页岩储层含气性分析方法，其特征在于，所述地层压力系数体为
7.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，归一化处理的公式为其中，XN为归一化处理后的参数，X为需要进行归一化处理的参数值；Xmin为归一化处理参数的最小值，Xmax为归一化处理参数的最大值。
8.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述页岩含气性综合评价指数为C＝(X1，X2，X3，X4)(A1，A2，A3，A4)T，其中，X1为归一化处理后有机含碳量数据体所对应的归一化值，X2为归一化处理后孔隙度数据体所对应的归一化值，X3为归一化处理后含气量数据体所对应的归一化值，X4为归一化处理后地层压力系数体所对应的归一化值，A1为有机含碳量所对应的权重，A2为孔隙度所对应的权重，A3为含气量所对应的权重，A4为地层压力系数体所对应的权重。
9.根据权利要求1所述的页岩储层含气性分析方法，其特征在于，所述选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数，包括：
选取与对应的岩心实测散点数据之间的关联系数大于0.5的弹性参数曲线作为所述敏感弹性参数曲线。
10.一种页岩储层含气性分析装置，应用于用户终端，其特征在于，包括：运算模块、选取模块和归一化模块；
所述运算模块用于依据横波测井曲线计算得到页岩储层随深度变化的多个弹性参数曲线，所述多个弹性参数曲线包括纵波阻抗曲线、横波阻抗曲线、杨式模量曲线、泊松比曲线和纵横波速度比曲线；
所述运算模块还用于依据积差方法求取通过岩心实测散点数据与弹性参数曲线上对应弹性参数之间的关联系数，所述岩心实测散点数据包括有机含碳量、孔隙度和含气量；
所述选取模块用于选取与岩心实测散点数据之间的关联系数大于预设阈值的弹性参数作为敏感弹性参数；
所述运算模块还用于以敏感弹性参数为自变量、与所述敏感弹性参数对应的岩心实测散点数据为因变量，通过最小二乘法确定与所述敏感弹性参数对应的岩心实测散点数据最佳匹配的多元回归方程；
所述运算模块还用于基于横波测井曲线和三维叠前道集数据通过叠前同时反演技术得到弹性参数体数据，所述弹性参数数据体包括纵波阻抗三维数据体、横波阻抗三维数据体和密度三维数据体；
所述运算模块还用于依据所述弹性参数数据体得到对应的敏感弹性参数数据体；
所述运算模块还用于依据所述敏感弹性参数数据体和所述多元回归方程，得到页岩储层的有机含碳量数据体、孔隙度数据体和含气量数据体；
所述运算模块还用于利用测得的泥质含量对页岩储层的声波速度曲线进行过滤，并依据页岩储层随深度变化的声波速度趋势，建立单井正常压实地层的声波速度趋势线；
所述运算模块还用于依据所述声波速度趋势线得到单井地层孔隙压力，并依据实际测得的地层孔隙压力数据得到区域伊顿指数；
所述运算模块还用于依据所述声波速度趋势线和所述区域伊顿指数，得到地层孔隙压力数据体；
所述运算模块还用于依据地层孔隙压力数据体和静水压力数据体得到地层压力系数体；
所述归一化模块用于采用极值法对对有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行归一化处理；
所述运算模块还用于对归一化处理后的有机含碳量数据体、孔隙度数据体、含气量数据体和地层压力系数体进行加权运算，得到页岩储层含气性综合评价指数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810917461.3A CN108594305A (zh) | 2018-08-13 | 2018-08-13 | 页岩储层含气性分析方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810917461.3A CN108594305A (zh) | 2018-08-13 | 2018-08-13 | 页岩储层含气性分析方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108594305A true CN108594305A (zh) | 2018-09-28 |
Family
ID=63618968
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810917461.3A Pending CN108594305A (zh) | 2018-08-13 | 2018-08-13 | 页岩储层含气性分析方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108594305A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114462237A (zh) * | 2022-01-27 | 2022-05-10 | 中国石油化工股份有限公司 | 一种泥页岩油气甜点预测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105626056A (zh) * | 2015-12-22 | 2016-06-01 | 中国石油天然气集团公司 | 一种确定目标区三维地层孔隙压力的方法和装置 |
| CN107605469A (zh) * | 2016-07-12 | 2018-01-19 | 中国石油化工股份有限公司 | 预测地层孔隙压力的方法 |
| CN108303510A (zh) * | 2017-12-25 | 2018-07-20 | 中国石油天然气股份有限公司 | 页岩气储层性能的评价方法、装置及计算机存储介质 |
-
2018
- 2018-08-13 CN CN201810917461.3A patent/CN108594305A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| CN105626056A (zh) * | 2015-12-22 | 2016-06-01 | 中国石油天然气集团公司 | 一种确定目标区三维地层孔隙压力的方法和装置 |
| CN107605469A (zh) * | 2016-07-12 | 2018-01-19 | 中国石油化工股份有限公司 | 预测地层孔隙压力的方法 |
| CN108303510A (zh) * | 2017-12-25 | 2018-07-20 | 中国石油天然气股份有限公司 | 页岩气储层性能的评价方法、装置及计算机存储介质 |
Non-Patent Citations (7)
| Title |
|---|
| 《地质监督与录井手册》编辑委员会: "《地质监督与录井手册》", 31 August 2001, 石油工业出版社 * |
| 《首届中国石油工业录井技术交流会论文集》编委会: "《首届中国石油工业录井技术交流会论文集2011》", 30 September 2011, 中国石油大学出版社 * |
| 罗兵等: "四川盆地涪陵地区页岩气层构造特征与保存评价", 《石油实验地质》 * |
| 胜亚楠等: "基于钻前风险预测的井身结构优化方法", 《石油钻采工艺》 * |
| 贺振华等: "《复杂油气藏地震波场特征方法理论及应用》", 31 March 1999, 四川科学技术出版社 * |
| 郝瑞: "《钻井工程》", 31 December 1989, 石油工业出版社 * |
| 陈亚琳: "涪陵页岩气田焦石坝地区含气量定量预测技术", 《海洋地质前沿》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114462237A (zh) * | 2022-01-27 | 2022-05-10 | 中国石油化工股份有限公司 | 一种泥页岩油气甜点预测方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103884633B (zh) | 一种确定岩石渗透率的方法及装置 | |
| CN103993871B (zh) | 针对薄互层地层的测井资料标准化处理方法及装置 | |
| US20140149041A1 (en) | Rock facies prediction in non-cored wells from cored wells | |
| CN105863628B (zh) | 一种油田开发阶段储层精细预测方法 | |
| CN105221133A (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
| CN108897066A (zh) | 碳酸盐岩裂缝密度定量预测方法及装置 | |
| CN102232198A (zh) | 处理井眼核磁共振测井记录以增强t2分布的连续性的方法 | |
| CN108680955B (zh) | 一种岩性识别方法及识别系统 | |
| CN108181649B (zh) | 一种介质分界面振幅补偿方法及装置 | |
| CN107450108A (zh) | 甜点区的确定方法和装置 | |
| CN114254505A (zh) | 渗透率模型的构建方法、装置、电子设备及存储介质 | |
| CN108594305A (zh) | 页岩储层含气性分析方法及装置 | |
| US8244472B2 (en) | Enabling well-logging instrument development on a common, realistic earth model | |
| US11346968B2 (en) | Estimation of reservoir flow properties from seismic data | |
| CN108693557A (zh) | 一种地震数据处理方法及装置 | |
| Christensen et al. | Seismically driven reservoir characterization using an innovative integrated approach: Syd Arne Field | |
| Xing et al. | Development and validation of an acoustic-electrical joint testing system for hydrate-bearing porous media | |
| CN110441822A (zh) | 一种地震数据去噪方法、装置、电子设备及存储介质 | |
| CN108490490A (zh) | 裂缝带表征方法、装置、电子设备及计算机存储介质 | |
| CN107991705A (zh) | 基于二维统计特征的测井曲线校正方法和装置 | |
| CN110119587B (zh) | 一种地层分层方法及装置 | |
| CN109100802B (zh) | 储层确定方法和装置 | |
| CN108761533B (zh) | 一种确定纵横波速度比的方法、装置及系统 | |
| CN107976711A (zh) | 一种地震信息定量预测烃源岩有机碳质量分数方法及装置 | |
| CN112666610A (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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20200416 Address after: 100027 Beijing, Chaoyangmen, North Street, No. 22, No. Applicant after: China Petroleum & Chemical Corp. Applicant after: EXPLORATION AND DEVELOPMENT RESEARCH INSTITUTE OF SINOPEC JIANGHAN OILFIELD Co. Address before: 430000 18 University Garden Road, East Lake hi tech Development Zone, Wuhan, Hubei (China Science and Technology Park) Applicant before: EXPLORATION AND DEVELOPMENT RESEARCH INSTITUTE OF SINOPEC JIANGHAN OILFIELD Co. |
|
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180928 |