CN111596349B - 基岩含气储层的预测方法及系统 - Google Patents

基岩含气储层的预测方法及系统 Download PDF

Info

Publication number
CN111596349B
CN111596349B CN201910129545.5A CN201910129545A CN111596349B CN 111596349 B CN111596349 B CN 111596349B CN 201910129545 A CN201910129545 A CN 201910129545A CN 111596349 B CN111596349 B CN 111596349B
Authority
CN
China
Prior art keywords
frequency
well
spectrum
gas
seismic data
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
Application number
CN201910129545.5A
Other languages
English (en)
Other versions
CN111596349A (zh
Inventor
于雪娇
张军勇
陈国文
周川江
王光华
郗小萍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201910129545.5A priority Critical patent/CN111596349B/zh
Publication of CN111596349A publication Critical patent/CN111596349A/zh
Application granted granted Critical
Publication of CN111596349B publication Critical patent/CN111596349B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • 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
    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (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

基岩含气储层的预测方法及系统
技术领域
本发明涉及地球物理勘探领域，具体地，涉及一种基岩含气储层的预测方法及系统。
背景技术
随着基岩勘探开发的不断深入，对基岩储层预测的需求也越来越高。目前，基岩储层预测的方法很多，但基岩含气性预测还处于探索阶段。基岩油气富集程度只依赖于裂缝和溶孔的发育程度，而基岩由于地层埋藏较深，又为塑性岩体，裂缝和溶孔并不是很发育，依据目前钻至基岩的测井解释成果，基岩孔隙度基本上都在10％以下，再加上并不是所有裂缝和溶孔都具有储集油气的能力，因此，基岩含气性预测存在一定难度。
发明内容
本发明实施例的主要目的在于提供一种基岩含气储层的预测方法及系统，以有效预测含气储层的平面分布范围。
为了实现上述目的，本发明实施例提供一种基岩含气储层的预测方法，包括：
采集宽频地震数据体；
分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
分析钻井的测井响应特征，得到总含气段位置；
在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据频谱能量体差值预测含气储层的平面分布范围。
本发明实施例还提供一种基岩含气储层的预测系统，包括：
采集单元，用于采集宽频地震数据体；
有效频带宽度单元，用于分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
总含气段位置单元，用于分析钻井的测井响应特征，得到总含气段位置；
频谱分析单元，用于在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
频率单元，用于确定每个钻井的频谱分析曲线的波峰对应的频率；
选取单元，用于从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
频谱分解单元，用于利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
频谱能量体差值单元，用于将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
预测单元，用于根据频谱能量体差值预测含气储层的平面分布范围。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：
采集宽频地震数据体；
分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
分析钻井的测井响应特征，得到总含气段位置；
在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据频谱能量体差值预测含气储层的平面分布范围。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
采集宽频地震数据体；
分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
分析钻井的测井响应特征，得到总含气段位置；
在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据频谱能量体差值预测含气储层的平面分布范围。
本发明实施例的基岩含气储层的预测方法及系统，先采集宽频地震数据体，分析宽频地震数据体得到有效频带宽度，再分析钻井的测井响应特征得到总含气段位置，并在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线，以确定每个钻井的频谱分析曲线的波峰对应的频率；然后从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；接着利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；最后将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值，根据频谱能量体差值预测含气储层的平面分布范围，可以有效预测含气储层的平面分布范围。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例中基岩含气储层的预测方法的流程图；
图2是本发明实施例中宽频地震数据体的频带宽度示意图；
图3是本发明实施例中井的测井响应特征示意图；
图4是本发明实施例中总含气段位置在地震剖面上的标定示意图；
图5是本发明实施例中井的频谱分析曲线示意图；
图6是本发明实施例中最小频率对应的单频频谱能量体的剖面示意图；
图7是本发明实施例中最大频率对应的单频频谱能量体的剖面示意图；
图8是本发明实施例中频谱能量体差值的剖面示意图；
图9是本发明实施例中频谱能量体差值的平面示意图；
图10是本发明实施例中基岩含气储层的预测系统的结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本领域技术技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。
鉴于现有技术无法预测基岩含气性，本发明实施例提供一种基岩含气储层的预测方法，以有效预测含气储层的平面分布范围。以下结合附图对本发明进行详细说明。
图1是本发明实施例中基岩含气储层的预测方法的流程图。如图1所示，基岩含气储层的预测方法包括：
S101：采集宽频地震数据体。其中，宽频地震数据体的格式为SEG-Y格式。数据采集时的最低激发频率可以达到1.5Hz。
S102：分析宽频地震数据体，得到宽频地震数据体的有效频带宽度。
图2是本发明实施例中宽频地震数据体的频带宽度示意图。图2的横轴为距离，单位为km；纵轴为频率，单位为赫兹(Hz)；JT2为JT2井，JT3为JT3井。如图2所示，宽频地震数据体的有效频带宽度为3Hz至60Hz。
S103：分析钻井的测井响应特征，得到总含气段位置。
S104：在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线。
S105：确定每个钻井的频谱分析曲线的波峰对应的频率。
S106：从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率。
S107：利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体。
S108：将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值。
S109：根据频谱能量体差值预测含气储层的平面分布范围。
图1所示的基岩含气储层的预测方法的执行主体可以为计算机。由图1所示的流程可知，本发明实施例的基岩含气储层的预测方法，先采集宽频地震数据体，分析宽频地震数据体得到有效频带宽度，再分析钻井的测井响应特征得到总含气段位置，并在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线，以确定每个钻井的频谱分析曲线的波峰对应的频率；然后从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；接着利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；最后将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值，根据频谱能量体差值预测含气储层的平面分布范围，可以有效预测含气储层的平面分布范围。
一实施例中，S103包括：分析每个钻井的测井响应特征，得到每个钻井的含气段位置；根据每个钻井的含气段位置求取平均值，将平均值作为总含气段位置。
一实施例中，测井响应特征包括：声波时差、补偿密度、纵波阻抗、自然伽马、孔隙度和含气饱和度。
图3是本发明实施例中井的测井响应特征示意图。如图3所示，分析JT1井的声波时差、补偿密度、纵波阻抗、自然伽马、孔隙度以及含气饱和度等测井响应特征，可以得出井的含气段位置为基岩顶面至基岩顶面下方的某一位置。本领域技术人员还可以分析JT2井、JT3井的测井响应特征，分别得到JT2井、JT3井的含气段位置。求取上述三个井的含气段位置的平均值，可以得到总含气段位置，为基岩顶面至基岩顶面下方的100米。
执行S103之后，可以利用浅层的标志层完成钻井的总含气段位置的标定。图4是本发明实施例中总含气段位置在地震剖面上的标定示意图。如图4所示，图4的横轴为地震道号；纵轴为深度，单位为米(m)；JT2为JT2井，JT3为JT3井。图4中的虚线为基岩顶面，实线为断层。其中，JT2井为气井，平均日产气5.1万立方米，JT3井为水井，累计产水976立方米。
图5是本发明实施例中井的频谱分析曲线示意图。如图5所示，图5的横轴为频率，单位为赫兹(Hz)；纵轴为振幅。JB101井为含气水井，该井在总含气段位置的频谱特征为双峰，且第一峰的低频振幅能量强于第二峰的高频振幅能量；JT2井为工业气流井，其频谱特征为单峰，且低频振幅能量较强；JT3井为水井，频谱特征为双峰，第一峰的振幅能量较低，第二峰的振幅能量较强。从图5中可以得到5个波峰对应的频率，而上述频率中的最小频率为10Hz，最大频率为22Hz，因此本领域技术人员可以利用宽频地震数据体进行频谱分解，分别得到10Hz对应的单频频谱能量体和22Hz对应的单频频谱能量体。
其中，可以应用小波变换算法对宽频地震数据体进行频谱分解，分别得到10Hz对应的单频频谱能量体和22Hz对应的单频频谱能量体。
本发明的其中一个具体实施例如下：
1、采集宽频地震数据体。其中，宽频地震数据体的格式为SEG-Y格式。
2、分析宽频地震数据体，得到宽频地震数据体的有效频带宽度。如图2所示，宽频地震数据体的有效频带宽度为3Hz至60Hz。
3、分析每个钻井的声波时差、补偿密度、纵波阻抗、自然伽马、孔隙度和含气饱和度，得到每个钻井的含气段位置；根据每个钻井的含气段位置求取平均值，将平均值作为总含气段位置。
4、在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线。
5、确定每个钻井的频谱分析曲线的波峰对应的频率，从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率。如图5所示，5个波峰对应的频率中的最小频率为10Hz，最大频率为22Hz。
6、应用小波变换算法对宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体。
7、将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值。
8、根据频谱能量体差值预测含气储层的平面分布范围。图6是本发明实施例中最小频率对应的单频频谱能量体的剖面示意图。图7是本发明实施例中最大频率对应的单频频谱能量体的剖面示意图。图8是本发明实施例中频谱能量体差值的剖面示意图。图9是本发明实施例中频谱能量体差值的平面示意图。如图6至图9所示，图6至图8的横轴为地震道号；纵轴为深度，单位为米。图9的横轴为东坐标，单位为米；纵轴为北坐标，单位为米。图6至图9中的颜色代表频谱能量体差值，颜色越浅表示频谱能量体差值越大。图9中的JB1-4井、JT1井均为基岩气井，其平面分布规律与含气预测结果相吻合，运算所得结果与钻井吻合率达90％以上。
综上，本发明实施例的基岩含气储层的预测方法先采集宽频地震数据体，分析宽频地震数据体得到有效频带宽度，再分析钻井的测井响应特征得到总含气段位置，并在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线，以确定每个钻井的频谱分析曲线的波峰对应的频率；然后从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；接着利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；最后将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值，根据频谱能量体差值预测含气储层的平面分布范围，可以有效预测含气储层的平面分布范围。
基于同一发明构思，本发明实施例还提供了一种基岩含气储层的预测系统，由于该系统解决问题的原理与基岩含气储层的预测方法相似，因此该系统的实施可以参见方法的实施，重复之处不再赘述。
图10是本发明实施例中基岩含气储层的预测系统的结构框图。如图10所示，基岩含气储层的预测系统包括：
采集单元，用于采集宽频地震数据体；
有效频带宽度单元，用于分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
总含气段位置单元，用于分析钻井的测井响应特征，得到总含气段位置；
频谱分析单元，用于在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
频率单元，用于确定每个钻井的频谱分析曲线的波峰对应的频率；
选取单元，用于从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
频谱分解单元，用于利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
频谱能量体差值单元，用于将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
预测单元，用于根据频谱能量体差值预测含气储层的平面分布范围。
在其中一种实施例中，总含气段位置单元具体用于：
分析每个钻井的测井响应特征，得到每个钻井的含气段位置；
根据每个钻井的含气段位置求取平均值，将平均值作为总含气段位置。
在其中一种实施例中，测井响应特征包括：声波时差、补偿密度、纵波阻抗、自然伽马、孔隙度和含气饱和度。
在其中一种实施例中，频谱分解单元具体用于：
应用小波变换算法对宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体。
综上，本发明实施例的基岩含气储层的预测系统先采集宽频地震数据体，分析宽频地震数据体得到有效频带宽度，再分析钻井的测井响应特征得到总含气段位置，并在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线，以确定每个钻井的频谱分析曲线的波峰对应的频率；然后从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；接着利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；最后将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值，根据频谱能量体差值预测含气储层的平面分布范围，可以有效预测含气储层的平面分布范围。
本发明实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现以下步骤：
采集宽频地震数据体；
分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
分析钻井的测井响应特征，得到总含气段位置；
在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据频谱能量体差值预测含气储层的平面分布范围。
综上，本发明实施例的计算机设备先采集宽频地震数据体，分析宽频地震数据体得到有效频带宽度，再分析钻井的测井响应特征得到总含气段位置，并在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线，以确定每个钻井的频谱分析曲线的波峰对应的频率；然后从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；接着利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；最后将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值，根据频谱能量体差值预测含气储层的平面分布范围，可以有效预测含气储层的平面分布范围。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
采集宽频地震数据体；
分析宽频地震数据体，得到宽频地震数据体的有效频带宽度；
分析钻井的测井响应特征，得到总含气段位置；
在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；
将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据频谱能量体差值预测含气储层的平面分布范围。
综上，本发明实施例的计算机可读存储介质先采集宽频地震数据体，分析宽频地震数据体得到有效频带宽度，再分析钻井的测井响应特征得到总含气段位置，并在有效频带宽度内对每个钻井在总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线，以确定每个钻井的频谱分析曲线的波峰对应的频率；然后从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；接着利用宽频地震数据体进行频谱分解，分别得到最小频率对应的单频频谱能量体和最大频率对应的单频频谱能量体；最后将最小频率对应的单频频谱能量体与最大频率对应的单频频谱能量体相减，得到频谱能量体差值，根据频谱能量体差值预测含气储层的平面分布范围，可以有效预测含气储层的平面分布范围。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (8)

1.一种基岩含气储层的预测方法，其特征在于，包括：
采集宽频地震数据体；
分析所述宽频地震数据体，得到所述宽频地震数据体的有效频带宽度；
分析每个钻井的测井响应特征，得到每个钻井的含气段位置；
根据每个钻井的含气段位置求取平均值，将所述平均值作为总含气段位置；
在所述有效频带宽度内对每个钻井在所述总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体；
将所述最小频率对应的单频频谱能量体与所述最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据所述频谱能量体差值预测含气储层的平面分布范围。
2.根据权利要求1所述的基岩含气储层的预测方法，其特征在于，所述测井响应特征包括：声波时差、补偿密度、纵波阻抗、自然伽马、孔隙度和含气饱和度。
3.根据权利要求1所述的基岩含气储层的预测方法，其特征在于，利用所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体，包括：
应用小波变换算法对所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体。
4.一种基岩含气储层的预测系统，其特征在于，包括：
采集单元，用于采集宽频地震数据体；
有效频带宽度单元，用于分析所述宽频地震数据体，得到所述宽频地震数据体的有效频带宽度；
总含气段位置单元，用于分析每个钻井的测井响应特征，得到每个钻井的含气段位置；根据每个钻井的含气段位置求取平均值，将所述平均值作为总含气段位置；
频谱分析单元，用于在所述有效频带宽度内对每个钻井在所述总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
频率单元，用于确定每个钻井的频谱分析曲线的波峰对应的频率；
选取单元，用于从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
频谱分解单元，用于利用所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体；
频谱能量体差值单元，用于将所述最小频率对应的单频频谱能量体与所述最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
预测单元，用于根据所述频谱能量体差值预测含气储层的平面分布范围。
5.根据权利要求4所述的基岩含气储层的预测系统，其特征在于，所述测井响应特征包括：声波时差、补偿密度、纵波阻抗、自然伽马、孔隙度和含气饱和度。
6.根据权利要求4所述的基岩含气储层的预测系统，其特征在于，所述频谱分解单元具体用于：
应用小波变换算法对所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体。
7.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现以下步骤：
采集宽频地震数据体；
分析所述宽频地震数据体，得到所述宽频地震数据体的有效频带宽度；
分析每个钻井的测井响应特征，得到每个钻井的含气段位置；
根据每个钻井的含气段位置求取平均值，将所述平均值作为总含气段位置；
在所述有效频带宽度内对每个钻井在所述总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体；
将所述最小频率对应的单频频谱能量体与所述最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据所述频谱能量体差值预测含气储层的平面分布范围。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现以下步骤：
采集宽频地震数据体；
分析所述宽频地震数据体，得到所述宽频地震数据体的有效频带宽度；
分析每个钻井的测井响应特征，得到每个钻井的含气段位置；
根据每个钻井的含气段位置求取平均值，将所述平均值作为总含气段位置；
在所述有效频带宽度内对每个钻井在所述总含气段位置的井旁道进行频谱分析，得到每个钻井的频谱分析曲线；
确定每个钻井的频谱分析曲线的波峰对应的频率；
从每个钻井的频谱分析曲线的波峰对应的频率中选取最小频率和最大频率；
利用所述宽频地震数据体进行频谱分解，分别得到所述最小频率对应的单频频谱能量体和所述最大频率对应的单频频谱能量体；
将所述最小频率对应的单频频谱能量体与所述最大频率对应的单频频谱能量体相减，得到频谱能量体差值；
根据所述频谱能量体差值预测含气储层的平面分布范围。
CN201910129545.5A 2019-02-21 2019-02-21 基岩含气储层的预测方法及系统 Active CN111596349B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910129545.5A CN111596349B (zh) 2019-02-21 2019-02-21 基岩含气储层的预测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910129545.5A CN111596349B (zh) 2019-02-21 2019-02-21 基岩含气储层的预测方法及系统

Publications (2)

Publication Number Publication Date
CN111596349A CN111596349A (zh) 2020-08-28
CN111596349B true CN111596349B (zh) 2023-05-26

Family

ID=72184799

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910129545.5A Active CN111596349B (zh) 2019-02-21 2019-02-21 基岩含气储层的预测方法及系统

Country Status (1)

Country Link
CN (1) CN111596349B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107272064A (zh) * 2017-07-18 2017-10-20 中国石油化工股份有限公司 碳酸盐岩缝洞体内部结构的刻画方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7729862B2 (en) * 2006-06-09 2010-06-01 Spectraseis Ag Frequency autonormalisation reservoir mapping
US8467266B2 (en) * 2006-06-13 2013-06-18 Seispec, L.L.C. Exploring a subsurface region that contains a target sector of interest
CA2818641C (en) * 2010-11-23 2017-07-04 Apex Spectral Technology, Inc. Methods of hydrocarbon detection using spectra dominant frequency and measures of energy decay on the low side and high side of spectra dominant frequency
CN102692647B (zh) * 2011-03-23 2014-12-03 中国石油天然气集团公司 一种高时间分辨率的地层含油气性预测方法
CN104516019B (zh) * 2013-09-27 2017-05-31 中国石油天然气集团公司 一种基于频谱形状的油气预测方法
CN104318106B (zh) * 2014-10-24 2020-07-21 西北大学 一种利用分频能量的衰减预测含气性和产能的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107272064A (zh) * 2017-07-18 2017-10-20 中国石油化工股份有限公司 碳酸盐岩缝洞体内部结构的刻画方法

Also Published As

Publication number Publication date
CN111596349A (zh) 2020-08-28

Similar Documents

Publication Publication Date Title
CN101575970B (zh) 一种随钻岩性和油藏特性识别方法
CN110847901B (zh) 一种变矿化度地层水下的致密砂岩储层流体识别方法
CN103698811B (zh) 一种碳酸盐岩岩石结构组分测井定量识别方法及其用途
CN113685172B (zh) 一种随钻声波固井质量评价方法及处理装置
CN108694264A (zh) 一种确定页岩气储层渗透率的方法
CN105064986A (zh) 利用常规测录井资料构建储层四性关系谱的方法
CN113887078B (zh) 深层-超深层海相碳酸盐岩油藏赋存深度下限预测方法
CN116122801A (zh) 一种页岩油水平井体积压裂可压性综合评价方法
CN111381292B (zh) 一种预测砂岩含烃储层的测井解释方法与装置
CN112363242A (zh) 基于测录井融合的储层流体识别方法与装置
CN106501146B (zh) 一种致密油藏物性上限确定方法及装置
CN111596349B (zh) 基岩含气储层的预测方法及系统
CN106934725B (zh) 岩石储层中值半径预测模型建立方法、装置及系统
CN110322363A (zh) 页岩气储层改造体积计算方法及系统
CN106568918B (zh) 页岩有机碳含量toc预测方法
US9989667B2 (en) Pore size classification in subterranean formations based on nuclear magnetic resonance (NMR) relaxation distributions
CN113530533B (zh) 基于三品质参数的甜点区预测方法及系统
CN113327070B (zh) 智能勘查煤系气的方法、装置及电子设备
CN108756867B (zh) 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法
CN111596351A (zh) 碳酸盐岩输导体系定量评价方法、系统、装置及存储介质
RU2259575C1 (ru) Способ определения продуктивности нефтяного пласта в трехмерном межскважинном пространстве
CN112415588A (zh) 基于三维地震网格的储层参数油气储量计算方法及系统
CN104699975A (zh) 从声电效应井下探测器测量数据中提取参数的方法
CN113900143B (zh) 天然气水合物的属性特征的确定方法、装置和服务器
CN113376692B (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