CN110988986B - 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法 - Google Patents

改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法 Download PDF

Info

Publication number
CN110988986B
CN110988986B CN201911353724.3A CN201911353724A CN110988986B CN 110988986 B CN110988986 B CN 110988986B CN 201911353724 A CN201911353724 A CN 201911353724A CN 110988986 B CN110988986 B CN 110988986B
Authority
CN
China
Prior art keywords
seismic
low
seismic data
frequency
frequency information
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
CN201911353724.3A
Other languages
English (en)
Other versions
CN110988986A (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.)
Chengdu Univeristy of Technology
Original Assignee
Chengdu Univeristy of Technology
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 Chengdu Univeristy of Technology filed Critical Chengdu Univeristy of Technology
Priority to CN201911353724.3A priority Critical patent/CN110988986B/zh
Publication of CN110988986A publication Critical patent/CN110988986A/zh
Priority to US17/133,649 priority patent/US11262470B2/en
Application granted granted Critical
Publication of CN110988986B publication Critical patent/CN110988986B/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
    • 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • 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/282—Application of seismic models, synthetic seismograms
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
    • G01V2210/23—Wavelet filtering

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

改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法
技术领域
本发明属于石油地震勘探领域，涉及一种改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法，对于深层碳酸盐岩缝洞型储层和礁滩储层的地震反射资料，通过本专利的方法处理后，可以显著改善深层碳酸盐储层中各种有效地质特征的成像精度。
背景技术
目前，采用常规震源和检波器采集的地震数据，由于震源和接收器的硬件功能限制，导致采集到的地震数据往往会缺失低频成分。尤其是在我国的很多地区，由于碳酸盐岩储层埋藏很深、其地质构造和空间分布十分复杂，且内部结构非均质性和各向异性强，使地震勘探信号的吸收衰减问题十分严重，深层碳酸盐储层的地震反射信号的振幅弱、主频低、频带窄、信噪比低、从而造成深层碳酸盐岩储层的反射结构不清，地层横向展布和可追踪性差，断裂系统的刻画和储层裂缝发育区的预测存在很大的多解性和不确定性。
充分可靠的地震低频信息对于深层油气储层勘探具有重要意义，对于深层碳酸盐岩储层描述，当地震资料中的低频分量丰富、低频端的能量较强时，可显著改善断层和裂缝系统等地质异常信息的成像质量。地震波在黏弹性介质中传播时，由于地层吸收、能量衰减因素导致地震信号的高频成分更易被吸收，而相对于高频分量来说，地震信号的低频成分衰减程度相对较小、传播更深，因此，充分利用地震资料中低频端的有效信息，有利于改善深部地层地震反射波的分辨率和深层地质体的成像质量。另外，地震信号低频端的成分越完整、能量越强，则地震反射波中的子波旁瓣越少，且旁瓣相对于主瓣的能量越弱，也有利于提高地震资料的分辨率。目前，一些常规的增强低频信息的方法通常是基于已有的地震数据构建补偿算子(亦称拓频算子)，这其中涉及的控制参数较多，使得这类方法在应用于实际资料处理时，调整参数的过程非常繁琐，难以获得最佳的处理效果。
发明内容
本发明的目的在于提供一种有效地改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法，所述方法包括以下主要步骤：
⑴输入地震数据集，首先利用地震数据的统计信息反演出地震子波，再用最优化稀疏正则化方法反演出反射系数；
⑵按如下步骤获得低频信息丰富的地震子波wb：
①基于步骤⑴中反演出的地震子波，计算输入地震数据集的平均地震子波wa；
②对平均地震子波wa做N点傅里叶变换，并求得其归一化的振幅谱Sa，N为地震道包含的采样点数量；
③确定参考位置n＝N/2，给定控制参数Pa和Pb；
④按照下式计算新的振幅谱Sb：
Figure BDA0002335339560000021
式中，j为采样点号，h是归一化振幅谱Sa中最大值对应的采样点号；
⑤按照下式计算转换系数C：
Figure BDA0002335339560000022
其中，
Figure BDA0002335339560000023
表示基本积运算；
⑥按照下式计算中间地震子波w：
Figure BDA0002335339560000024
式中，ft(·)表示傅里叶变换，ift(·)表示反傅里叶变换，real(·)表示取实部运算；
⑦按照下式计算比例变换系数λ：
Figure BDA0002335339560000025
式中，sum(·)表示求和运算；
⑧按照下式计算得到低频信息丰富的地震子波wb：
Figure BDA0002335339560000026
⑶将低频信息丰富的地震子波wb与步骤⑴中反演得到的反射系数进行褶积运算，最终得到低频信息丰富和低频能量增强的地震数据集。
还需要说明的是，步骤⑵中控制参数Pa的取值范围为
Figure BDA0002335339560000027
值越小，增强的低频成分越多，值越大，增强的低频成分越少。步骤⑵中控制参数Pb的取值范围为[1，3]，值越小，所增强低频成分的能量越小，值越大，所增强低频成分的能量越大。
附图说明
图1是本发明实施例的某工区地震数据的平均地震子波wa和对应的低频信息丰富地震子波wb。其中，图1(a)是输入地震数据集的平均地震子波wa，纵坐标为振幅，横坐标为时间，单位是毫秒(ms)。图1(b)是平均地震子波wa的振幅谱，纵坐标为振幅，横坐标为频率，单位是赫兹(Hz)。图1(c)是低频信息丰富的地震子波wb，纵坐标为振幅，横坐标为时间，单位是毫秒(ms)。图1(d)低频信息丰富的地震子波wb的振幅谱，纵坐标为振幅，横坐标为频率，单位是赫兹(Hz)。
图2是本发明实施例对某工区地震数据集用图1(c)中的低频信息丰富的地震子波wb低频增强处理前后的地震剖面对比。其中，图2(a)是低频信息增强处理前的原始地震剖面，横坐标为道号，纵坐标为时间，单位是秒(s)。图2(b)是利用本发明的方法在低频信息增强处理后的地震剖面，横坐标为道号，纵坐标为时间，单位是秒(s)。
图3是本发明实施例对某工区低频信息增强处理前后的地震数据集，分别计算方差属性后的时间切片对比，用于说明本发明的方法在处理前后，地震数据中深层地质异常信息的成像质量的差异。其中，图3(a)是低频信息增强处理前，原始地震数据的方差属性在4452毫秒处的时间切片，横坐标为InLine线号，纵坐标为XLine线号。图3(b)是利用本发明的方法进行低频信息增强处理后，地震数据的方差属性在4452毫秒处的时间切片，横坐标为InLine线号，纵坐标为XLine线号。
具体实施方式
⑴输入地震数据集，本实施例中首先用地震数据的统计信息反演得到地震子波，再用最优化稀疏正则化方法反演出反射系数；
⑵按以下步骤构建低频信息丰富的地震子波wb：
①基于步骤⑴中反演出的地震子波，计算输入地震数据集的平均地震子波wa；
②对平均地震子波wa做N点傅里叶变换，并求得其归一化的振幅谱Sa，N为地震道包含的采样点数量，本实施例中N＝1750；
③确定参考位置n＝N/2，给定控制参数Pa＝0.07和Pb＝1.5；
④按照下式计算新的振幅谱Sb：
Figure BDA0002335339560000041
式中，j为采样点号，h是归一化振幅谱Sa中最大值对应的采样点号；
⑤按照下式计算转换系数C：
Figure BDA0002335339560000042
式中，
Figure BDA0002335339560000043
表示基本积运算；
⑥按照下式计算中间地震子波w：
Figure BDA0002335339560000044
式中，ft(·)表示傅里叶变换，ift(·)表示反傅里叶变换，real(·)表示取实部运算；
⑦按照下式计算比例变换系数λ：
Figure BDA0002335339560000045
式中，sum(·)表示求和运算；
⑧按照下式计算得到低频信息丰富的地震子波wb：
Figure BDA0002335339560000046
⑶将低频信息丰富的地震子波wb与步骤⑴中反演得到的反射系数进行褶积运算，最终得到低频信息丰富和低频能量增强的地震数据集。
图1(c)和1(d)为本发明实施例基于某工区地震数据集的平均地震子波wa，控制参数给为Pa＝0.07、Pb＝1.5时对应的低频信息丰富地震子波wb及其归一化振幅谱。
图2为本发明实施例对某工区地震数据集用图1(c)中的低频信息丰富地震子波wb，执行低频补偿处理前后地震剖面的对比。从图中可以看出，经过低频信息增强处理后的图2(b)对于深层碳酸盐岩断裂系统的成像质量优于图2(a)的原始地震剖面(如图2中白色虚线矩形框所示区域)。
图3为本发明实施例对某工区原始地震数据集，用本发明的方法进行低频信息增强处理后得到低频信息增强的地震数据集，分别对处理前的原始地震数据集和处理后的地震数据集计算方差属性，然后从它们的方差属性数据集中分别抽取4452毫秒处的时间切片作对比。对比图3(b)和图3(a)可见，本发明的方法进行低频信息增强处理后，图3(b)的方差属性时间切片中更清晰地显示了更多、更丰富的地质异常信息。
本发明的优点在于：(1)利用输入的整个地震数据集的平均地震子波来构建低频信息丰富的地震子波，考虑了工区内地震数据在横向和纵向上的变化；(2)仅包含2个控制参数，通过参数Pa可以控制增强多少低频成分，通过参数Pb可以控制低频成分能量增强程度的大小，在实际地震资料处理的应用中，可以快速根据工区资料的实际情况确定处理参数，获得最佳的处理效果。
上述各实施例仅用于说明本发明，其中方法的各实施步骤等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (1)

1.一种改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法，其包括以下主要步骤：
⑴输入地震数据集，首先利用地震数据的统计信息反演出地震子波，再用最优化稀疏正则化方法反演出地震反射系数；
⑵按以下步骤构建低频信息丰富的地震子波wb：
①基于步骤⑴中反演出的地震子波，计算输入地震数据集的平均地震子波wa；
②对平均地震子波wa做N点傅里叶变换，并求得其归一化的振幅谱Sa，N为地震道包含的采样点数量；
③确定参考位置n＝N/2，给定控制参数Pa和Pb；
④按照下式计算新的振幅谱Sb：
Figure FDA0002734894110000011
式中，j为采样点号，h是归一化振幅谱Sa中最大值对应的采样点号；控制参数Pa的取值范围为
Figure FDA0002734894110000012
Pa的值越小，增强的低频成分越多，Pa的值越大，增强的低频成分越少；控制参数Pb的取值范围为[1,3]，Pb的值越小，所增强低频成分的能量越小，Pb的值越大，所增强低频成分的能量越大；
⑤按照下式计算转换系数C：
Figure FDA0002734894110000013
其中，
Figure FDA0002734894110000015
表示基本积运算；
⑥按照下式计算中间地震子波w：
Figure FDA0002734894110000016
式中，ft(·)表示傅里叶变换，ift(·)表示反傅里叶变换，real(·)表示取实部运算；
⑦按照下式计算比例变换系数λ：
Figure FDA0002734894110000014
式中，sum(·)表示求和运算；
⑧按照下式计算得到低频信息丰富的地震子波wb：
Figure FDA0002734894110000021
⑶将低频信息丰富的地震子波wb与步骤⑴中反演得到的地震反射系数进行褶积运算，得到低频信息丰富和低频能量增强的地震数据集。
CN201911353724.3A 2019-12-25 2019-12-25 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法 Active CN110988986B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201911353724.3A CN110988986B (zh) 2019-12-25 2019-12-25 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法
US17/133,649 US11262470B2 (en) 2019-12-25 2020-12-24 Method of low-frequency seismic data enhancement for improving characterization precision of deep carbonate reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911353724.3A CN110988986B (zh) 2019-12-25 2019-12-25 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法

Publications (2)

Publication Number Publication Date
CN110988986A CN110988986A (zh) 2020-04-10
CN110988986B true CN110988986B (zh) 2021-01-01

Family

ID=70076414

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911353724.3A Active CN110988986B (zh) 2019-12-25 2019-12-25 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法

Country Status (2)

Country Link
US (1) US11262470B2 (zh)
CN (1) CN110988986B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113960671B (zh) * 2020-07-20 2024-04-16 中国石油化工股份有限公司 频变子波压缩处理方法、装置、计算机设备和存储介质
CN112433248B (zh) * 2020-11-10 2022-09-30 成都理工大学 一种碳酸盐岩沉积环境下隐蔽储层的探测方法
CN113391352B (zh) * 2021-06-11 2022-03-25 成都理工大学 一种突出含气储层低频地震异常的地震信号时频分析方法
CN115680593A (zh) * 2021-07-30 2023-02-03 中国石油天然气股份有限公司 火山油藏裂缝集中发育区的预测方法
CN114280671A (zh) * 2021-09-17 2022-04-05 成都理工大学 一种广义w变换地震含气储层表征方法
CN114002738B (zh) * 2021-11-02 2023-08-15 中海石油(中国)有限公司 高角度断裂约束的叠前裂缝预测方法、装置及存储介质
CN114663453B (zh) * 2022-03-11 2023-04-18 成都理工大学 储层多尺度裂缝网络定量描述方法及装置
CN115327624B (zh) * 2022-08-02 2023-07-14 西安交通大学 一种地震子波和反射系数的反演方法及反演系统
CN115639603B (zh) * 2022-10-25 2024-02-23 电子科技大学 基于采样点叠加去干涉的地震拓频方法、设备及存储介质
CN116840916B (zh) * 2023-07-04 2024-03-26 成都理工大学 一种地震速度信号和加速度信号联合子波提取方法
CN117148432B (zh) * 2023-10-27 2024-03-19 胜利信科(山东)勘察测绘有限公司 基于多尺度分量提取的浅剖数据空间插值方法
CN117890979B (zh) * 2024-03-14 2024-05-24 山东科技大学 地震数据自适应弱反射信号补偿方法、系统、设备及介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012134621A1 (en) * 2011-03-30 2012-10-04 Exxonmobil Upstream Research Company Convergence rate of full wavefield inversion using spectral shaping
CN102928879A (zh) * 2012-10-31 2013-02-13 中国石油天然气股份有限公司 一种弱地震反射特征的碳酸盐岩勘探方法及装置
CN104280765A (zh) * 2013-07-11 2015-01-14 中国石油化工股份有限公司 基于变子波反射系数反演的地震高分辨处理方法
CN106873038A (zh) * 2017-03-15 2017-06-20 成都理工大学 一种从深度域地震数据中提取深度域地震子波的方法
CN108020863A (zh) * 2017-11-28 2018-05-11 河海大学 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法
CN109633744A (zh) * 2019-01-18 2019-04-16 中海石油(中国)有限公司上海分公司 地震子波的提取方法、装置、设备及存储介质
CN109738951A (zh) * 2019-01-03 2019-05-10 国家深海基地管理中心 一种基于地震同相轴子波谱的时变反褶积方法
CN110146923A (zh) * 2019-07-03 2019-08-20 成都理工大学 一种高效的高精度深度域地震子波提取方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6131071A (en) * 1996-12-06 2000-10-10 Bp Amoco Corporation Spectral decomposition for seismic interpretation
US6745129B1 (en) * 2002-10-29 2004-06-01 The University Of Tulsa Wavelet-based analysis of singularities in seismic data
US7525873B1 (en) * 2004-07-23 2009-04-28 Bf Corporation North America Inc. Seismic inversion of conditioned amplitude spectra
WO2009088602A1 (en) * 2008-01-08 2009-07-16 Exxonmobil Upstream Research Company Spectral shaping inversion and migration of seismic data
CN103376464B (zh) * 2012-04-13 2016-04-06 中国石油天然气集团公司 一种地层品质因子反演方法
US10353099B2 (en) * 2014-10-24 2019-07-16 Ion Geophysical Corporation Methods and systems for seismic inversion and related seismic data processing
US10191168B2 (en) * 2016-01-12 2019-01-29 Cgg Services Sas AVA compliant pre-stack frequency spectrum enhancement of seismic data
US20180017692A1 (en) * 2016-07-18 2018-01-18 Cgg Services Sas Device and method for estimating pre-stack wavelet model from seismic gathers

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012134621A1 (en) * 2011-03-30 2012-10-04 Exxonmobil Upstream Research Company Convergence rate of full wavefield inversion using spectral shaping
CN102928879A (zh) * 2012-10-31 2013-02-13 中国石油天然气股份有限公司 一种弱地震反射特征的碳酸盐岩勘探方法及装置
CN104280765A (zh) * 2013-07-11 2015-01-14 中国石油化工股份有限公司 基于变子波反射系数反演的地震高分辨处理方法
CN106873038A (zh) * 2017-03-15 2017-06-20 成都理工大学 一种从深度域地震数据中提取深度域地震子波的方法
CN108020863A (zh) * 2017-11-28 2018-05-11 河海大学 一种基于地震奇偶函数的碳酸盐岩薄储层孔隙度预测方法
CN109738951A (zh) * 2019-01-03 2019-05-10 国家深海基地管理中心 一种基于地震同相轴子波谱的时变反褶积方法
CN109633744A (zh) * 2019-01-18 2019-04-16 中海石油(中国)有限公司上海分公司 地震子波的提取方法、装置、设备及存储介质
CN110146923A (zh) * 2019-07-03 2019-08-20 成都理工大学 一种高效的高精度深度域地震子波提取方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A procedure to reduce side lobes of reflectionwavelets: A contribution to low frequency information;Hakan Karslı，等;《Journal of Applied Geophysics》;20131231;第107-118页 *
Seismic low-frequency-based calculation of reservoir fluid mobility and its applications;chen xuehua，等;《APPLIED GEOPHYSICS》;20120930;第326-332页 *
基于时频谱模拟的时变混合相位子波提取;戴永寿，等;《石油地球物理勘探》;20151031;第830-838页 *

Also Published As

Publication number Publication date
US11262470B2 (en) 2022-03-01
US20210199830A1 (en) 2021-07-01
CN110988986A (zh) 2020-04-10

Similar Documents

Publication Publication Date Title
CN110988986B (zh) 改善深层碳酸盐岩储层刻画精度的地震资料低频增强方法
CN106597532B (zh) 一种结合井资料与层位资料的叠前地震数据频带拓展方法
Wu et al. Noise attenuation for 2-D seismic data by radial-trace time-frequency peak filtering
CN107422381B (zh) 一种基于eemd-ica的地震低频信息流体预测方法
CN101923176B (zh) 一种利用地震数据瞬时频率属性进行油气检测的方法
CN110554428A (zh) 一种基于变分模态分解的地震波低频能量变化率提取方法
CN107703546A (zh) 一种基于小波变换的新阈值函数地震资料去噪方法
CN111427088A (zh) 用于识别薄互储层的地震数据低频补偿方法
Zhou et al. Amplitude spectrum compensation and phase spectrum correction of seismic data based on the generalized S transform
CN107703547A (zh) 一种自适应小波阈值函数地震资料去噪方法
CN107728213A (zh) 一种小波新阈值函数地震资料去噪方法
Wang et al. Application of a new wavelet threshold method in unconventional oil and gas reservoir seismic data denoising
CN115390133A (zh) 一种基于压缩感知和统计学习的地震弱信号增强方法
CN109884705B (zh) 双重约束时频域子波谱提高地震分辨率处理方法
CN112684501B (zh) 一种基于谱比面积的q值估计方法与应用
CN110749923A (zh) 一种基于范数方程提高分辨率的反褶积方法
CN112180464B (zh) 一种储层物性的识别方法
CN114428280A (zh) 一种地震数据低频信息补偿方法及其应用
CN112147687B (zh) 一种储层含气性预测方法及预测系统
CN115639603B (zh) 基于采样点叠加去干涉的地震拓频方法、设备及存储介质
Erhan et al. Application of non-stationary iterative time-domain deconvolution
Kong et al. An amplitude suppression method based on the decibel criterion
CN110673211A (zh) 一种基于测井与地震数据的品质因子建模方法
CN114355442B (zh) 一种三参数w变换的地震储层识别时频分析方法
CN104237941A (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