CN109459788A - 地层品质因子计算方法及*** - Google Patents

地层品质因子计算方法及*** Download PDF

Info

Publication number
CN109459788A
CN109459788A CN201710796553.6A CN201710796553A CN109459788A CN 109459788 A CN109459788 A CN 109459788A CN 201710796553 A CN201710796553 A CN 201710796553A CN 109459788 A CN109459788 A CN 109459788A
Authority
CN
China
Prior art keywords
wavelet
quality factors
power spectrum
decaying
interval quality
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.)
Granted
Application number
CN201710796553.6A
Other languages
English (en)
Other versions
CN109459788B (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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710796553.6A priority Critical patent/CN109459788B/zh
Publication of CN109459788A publication Critical patent/CN109459788A/zh
Application granted granted Critical
Publication of CN109459788B publication Critical patent/CN109459788B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/50—Corrections or adjustments related to wave propagation
    • G01V2210/58—Media-related
    • G01V2210/584—Attenuation
    • 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

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

地层品质因子计算方法及***
技术领域
本发明涉及地震勘探领域，更具体地，涉及一种地层品质因子计算方法及***。
背景技术
品质因子是描述地下介质衰减特征的重要参数，准确可靠的品质因子估计在地震处理解释中具有重要的意义。目前工业界常用的品质因子估计方法主要可分为两大类，即时间域法和频率域法，每一类方法又包括很多种不同的方法。其中频率域方法广泛应用于实际地震勘探中，该类方法以谱比法和质心频率法为代表。谱比法通过拟合衰减前后对数谱比值与频率的线性关系来求取品质因子，其在无噪声时的理论估计结果很精确，然而其对噪声很敏感，因此当地震资料存在噪声时其估计结果的稳定性和精度较差。质心频率法通过建立衰减前后地震波的质心频率变化量与品质因子Q值之间的关系式来反演地层的品质因子，由于利用的是振幅谱的统计特性，故该方法抗噪性能较好，然而其假设震源子波满足某种特定的形状(如满足高斯分布)，因此当实际子波与假设不一致时，估计的品质因子存在一定的误差。因此，有必要开发一种地层品质因子计算方法及***。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种地层品质因子计算方法及***，具有理论计算精度高、稳定性与抗噪性能好的特点，其计算结果可用于后续的地震数据处理及解释，如地震吸收衰减补偿、天然气检测等。
根据本发明的一方面，提出了一种地层品质因子计算方法。所述方法可以包括：根据衰减子波的振幅谱，计算所述衰减子波的功率谱；将所述衰减子波的功率谱转化为对数形式；对对数形式的衰减子波功率谱进行积分运算，获取地层品质因子。
优选地，所述衰减子波的振幅谱为：
其中，A2(f)为衰减子波的振幅谱，A1(f)为参考子波的振幅谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
优选地，所述衰减子波的功率谱为：
其中，P2(f)为衰减子波的功率谱，P1(f)为参考子波的功率谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
优选地，所述衰减子波功率谱的对数形式为：
其中，为衰减子波功率谱的对数形式，为参考子波的功率谱的对数形式。
优选地，所述地层品质因子为：
根据本发明的另一方面，提出了一种地层品质因子计算***，可以包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：根据衰减子波的振幅谱，计算所述衰减子波的功率谱；将所述衰减子波的功率谱转化为对数形式；对对数形式的衰减子波功率谱进行积分运算，获取地层品质因子。
优选地，所述衰减子波的振幅谱为：
其中，A2(f)为衰减子波的振幅谱，A1(f)为参考子波的振幅谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
优选地，所述衰减子波的功率谱为：
其中，P2(f)为衰减子波的功率谱，P1(f)为参考子波的功率谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
优选地，所述衰减子波功率谱的对数形式为：
其中，为衰减子波功率谱的对数形式，为参考子波的功率谱的对数形式。
优选地，所述地层品质因子为：
本发明的有益效果在于：(1)未假设子波频谱的形状，其对任何子波都是适用的，并且计算精度较高；(2)利用功率谱，相比振幅谱具有更好的抗噪性，积分项相当于求和，可以压制噪声，具有较高的稳定性和抗噪性；(3)本发明估计的地层品质因子可以作为地震反Q滤波的输入，以补偿地震波的吸收衰减，由于气体具有很强的吸收衰减作用，其通常对应小的地层品质因子，所以该方法的估计结果也可以用于天然气检测。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了根据本发明的地层品质因子计算方法的步骤的流程图。
图2示出了根据本发明的一个实施例的地层品质因子为20、40、80和120时的衰减地震记录的示意图。
图3示出了根据本发明、质心频率法和谱比法计算的地层品质因子的相对误差的示意图。
图4示出了根据本发明的一个实施例的含噪声地震记录的示意图。
图5a、图5b和图5c分别示出了根据本发明、质心频率法和谱比法计算的地层品质因子的概率分布的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的地层品质因子计算方法的步骤的流程图。
在该实施例中，根据本发明的地层品质因子计算方法可以包括：步骤101，根据衰减子波的振幅谱，计算衰减子波的功率谱；步骤102，将衰减子波的功率谱转化为对数形式；步骤103，对对数形式的衰减子波功率谱进行积分运算，获取地层品质因子。
在一个示例中，衰减子波的振幅谱为：
其中，A2(f)为衰减子波的振幅谱，A1(f)为参考子波的振幅谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
在一个示例中，衰减子波的功率谱为：
其中，P2(f)为衰减子波的功率谱，P1(f)为参考子波的功率谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
在一个示例中，衰减子波功率谱的对数形式为：
其中，为衰减子波功率谱的对数形式，为参考子波的功率谱的对数形式。
在一个示例中，地层品质因子为：
具体地，由于功率谱具有更高的抗噪性，因此用功率谱进行品质因子计算。根据衰减子波的振幅谱为公式(1)，计算衰减子波的功率谱为公式(2)；分别对公式(2)两边取对数，将公式(2)转化为对数形式为公式(3)；为了提高地层品质因子计算的稳定性和抗噪性，对公式(3)两边关于频率f求积分，获取地层品质因子为公式(4)。
本方法具有理论计算精度高、稳定性与抗噪性能好的特点，其计算结果可用于后续的地震数据处理及解释，如地震吸收衰减补偿、天然气检测等。
应用示例
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
图2示出了根据本发明的一个实施例的地层品质因子为20、40、80和120时的衰减地震记录的示意图，150ms处为初始地震子波，450ms处为衰减后的地震波，可以看到，地质品质因子越小，地震波衰减越大。
建立一个均匀介质模型，并合成对应的衰减地震记录，模型介质纵波速度为2000m/s，传播距离为600m，选取的雷克子波主频为40Hz。根据衰减子波的振幅谱为公式(1)，计算衰减子波的功率谱为公式(2)；分别对公式(2)两边取对数，将公式(2)转化为对数形式为公式(3)；为了提高地层品质因子计算的稳定性和抗噪性，对公式(3)两边关于频率f求积分，获取地层品质因子为公式(4)。
图3示出了根据本发明、质心频率法和谱比法计算的地层品质因子的相对误差的示意图，可以看出，本发明提出的方法计算的地质品质因子的精度最高。
图4示出了根据本发明的一个实施例的含噪声地震记录的示意图。
图5a、图5b和图5c分别示出了根据本发明、质心频率法和谱比法计算的地层品质因子的概率分布的示意图，分别是统计2000个品质因子估计结果的概率分布值，其真实品质因子为20，可以看出，本发明提出的方法估计结果最接近于真实值，误差小、稳定性好；而质心频率法和谱比法的估计结果分布范围较大，说明其估计结果的方差较大，意味着误差较大、稳定性较差。
综上所述，本发明具有理论计算精度高、稳定性与抗噪性能好的特点，其计算结果可用于后续的地震数据处理及解释，如地震吸收衰减补偿、天然气检测等。
根据本发明的实施例，提供了一种地层品质因子计算***，可以包括：存储器，存储有计算机可执行指令；处理器，处理器运行存储器中的计算机可执行指令，执行以下步骤：根据衰减子波的振幅谱，计算衰减子波的功率谱；将衰减子波的功率谱转化为对数形式；对对数形式的衰减子波功率谱进行积分运算，获取地层品质因子。
在一个示例中，衰减子波的振幅谱为：
其中，A2(f)为衰减子波的振幅谱，A1(f)为参考子波的振幅谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
在一个示例中，衰减子波的功率谱为：
其中，P2(f)为衰减子波的功率谱，P1(f)为参考子波的功率谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
在一个示例中，衰减子波功率谱的对数形式为：
其中，为衰减子波功率谱的对数形式，为参考子波的功率谱的对数形式。
在一个示例中，地层品质因子为：
本发明具有理论计算精度高、稳定性与抗噪性能好的特点，其计算结果可用于后续的地震数据处理及解释，如地震吸收衰减补偿、天然气检测等。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。

Claims (10)

1.一种地层品质因子计算方法，包括：
根据衰减子波的振幅谱，计算所述衰减子波的功率谱；
将所述衰减子波的功率谱转化为对数形式；
对对数形式的衰减子波功率谱进行积分运算，获取地层品质因子。
2.根据权利要求1所述的地层品质因子计算方法，其中，所述衰减子波的振幅谱为：
其中，A2(f)为衰减子波的振幅谱，A1(f)为参考子波的振幅谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
3.根据权利要求1所述的地层品质因子计算方法，其中，所述衰减子波的功率谱为：
其中，P2(f)为衰减子波的功率谱，P1(f)为参考子波的功率谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
4.根据权利要求3所述的地层品质因子计算方法，其中，所述衰减子波功率谱的对数形式为：
其中，为衰减子波功率谱的对数形式，为参考子波的功率谱的对数形式。
5.根据权利要求4所述的地层品质因子计算方法，其中，所述地层品质因子为：
6.一种地层品质因子计算***，其特征在于，该***包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
根据衰减子波的振幅谱，计算所述衰减子波的功率谱；
将所述衰减子波的功率谱转化为对数形式；
对对数形式的衰减子波功率谱进行积分运算，获取地层品质因子。
7.根据权利要求6所述的地层品质因子计算***，其中，所述衰减子波的振幅谱为：
其中，A2(f)为衰减子波的振幅谱，A1(f)为参考子波的振幅谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
8.根据权利要求6所述的地层品质因子计算***，其中，所述衰减子波的功率谱为：
其中，P2(f)为衰减子波的功率谱，P1(f)为参考子波的功率谱，f为频率，T为透射系数，t为两个子波间的旅行时间，Q为地层品质因子。
9.根据权利要求8所述的地层品质因子计算***，其中，所述衰减子波功率谱的对数形式为：
其中，为衰减子波功率谱的对数形式，为参考子波的功率谱的对数形式。
10.根据权利要求9所述的地层品质因子计算***，其中，所述地层品质因子为：
CN201710796553.6A 2017-09-06 2017-09-06 地层品质因子计算方法及*** Active CN109459788B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710796553.6A CN109459788B (zh) 2017-09-06 2017-09-06 地层品质因子计算方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710796553.6A CN109459788B (zh) 2017-09-06 2017-09-06 地层品质因子计算方法及***

Publications (2)

Publication Number Publication Date
CN109459788A true CN109459788A (zh) 2019-03-12
CN109459788B CN109459788B (zh) 2020-05-15

Family

ID=65606000

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710796553.6A Active CN109459788B (zh) 2017-09-06 2017-09-06 地层品质因子计算方法及***

Country Status (1)

Country Link
CN (1) CN109459788B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117724165A (zh) * 2024-02-07 2024-03-19 中国石油大学(华东) 一种基于时变子波的品质因子估计方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4884247A (en) * 1987-03-09 1989-11-28 Mobil Oil Company Method of processing geophysical data to compensate for earth filter attenuation
US5555218A (en) * 1995-02-27 1996-09-10 Western Atlas International, Inc. Computation of Q-derived static corrections from seismic data
US20020141287A1 (en) * 2001-02-13 2002-10-03 Lazaratos Spyridon K. Method for spectral balancing of near-and far-offset seismic data
US20050086005A1 (en) * 2003-10-16 2005-04-21 Taner M. T. Method for determining formation quality factor from seismic data
US20060265132A1 (en) * 2005-05-13 2006-11-23 Chevron U.S.A. Inc. Method for estimation of interval seismic quality factor
CN102023311A (zh) * 2010-08-10 2011-04-20 中国石油大学(华东) 地层的品质因子谱及其求取方法
CN103728662A (zh) * 2014-01-03 2014-04-16 中国海洋石油总公司 一种基于地震信号包络峰值的地层介质品质因子估计方法
CN104570108A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 估算等效品质因子方法及用其估算地层品质因子的方法
CN105388518A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种质心频率与频谱比联合的井中地震品质因子反演方法
CN105388523A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种高精度的品质因子提取方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4884247A (en) * 1987-03-09 1989-11-28 Mobil Oil Company Method of processing geophysical data to compensate for earth filter attenuation
US5555218A (en) * 1995-02-27 1996-09-10 Western Atlas International, Inc. Computation of Q-derived static corrections from seismic data
US20020141287A1 (en) * 2001-02-13 2002-10-03 Lazaratos Spyridon K. Method for spectral balancing of near-and far-offset seismic data
US20050086005A1 (en) * 2003-10-16 2005-04-21 Taner M. T. Method for determining formation quality factor from seismic data
US20060265132A1 (en) * 2005-05-13 2006-11-23 Chevron U.S.A. Inc. Method for estimation of interval seismic quality factor
CN102023311A (zh) * 2010-08-10 2011-04-20 中国石油大学(华东) 地层的品质因子谱及其求取方法
CN104570108A (zh) * 2013-10-29 2015-04-29 中国石油化工股份有限公司 估算等效品质因子方法及用其估算地层品质因子的方法
CN103728662A (zh) * 2014-01-03 2014-04-16 中国海洋石油总公司 一种基于地震信号包络峰值的地层介质品质因子估计方法
CN105388518A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种质心频率与频谱比联合的井中地震品质因子反演方法
CN105388523A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种高精度的品质因子提取方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHANGJUN ZHANG，ET AL: "Estimation of quality factors from CMP records", 《GEOPHYSICS》 *
LUKE J. W. PINSON，等: "Estimating quality factor and mean grain size of sediments from high-resolution marine seismic data", 《GEOPHYSICS》 *
余连勇，等: "时间域质心频移法估算品质因子Q", 《西南石油大学学报（自然科学版）》 *
崔琴，等: "基于加权指数型子波的质心频移法估算品质因子", 《物探与化探》 *
曹思远，等: "对数谱根式法Q值反演", 《石油地球物理勘探》 *
杨登锋，等: "能量谱质心频移法Ｑ值估计", 《石油地球物理勘探》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117724165A (zh) * 2024-02-07 2024-03-19 中国石油大学(华东) 一种基于时变子波的品质因子估计方法
CN117724165B (zh) * 2024-02-07 2024-05-03 中国石油大学(华东) 一种基于时变子波的品质因子估计方法

Also Published As

Publication number Publication date
CN109459788B (zh) 2020-05-15

Similar Documents

Publication Publication Date Title
Yang et al. Application of optimal transport and the quadratic Wasserstein metric to full-waveform inversion
EP3073296B1 (en) Full waveform inversion method for seismic data processing using preserved amplitude reverse time migration
US11333783B2 (en) Integrated method for estimation of seismic wavelets and synthesis of seismic records in depth domain
CN107843925B (zh) 一种基于修正相位的反射波波形反演方法
CN111723329B (zh) 一种基于全卷积神经网络的震相特征识别波形反演方法
CN107272062B (zh) 一种数据驱动的地下介质q场估计方法
US9612351B2 (en) System and method for estimating and attenuating noise in seismic data
CN110914718B (zh) 计算频域中地震速度反演的与振幅无关的梯度
AU2015383134B2 (en) Multistage full wavefield inversion process that generates a multiple free data set
CN109669212B (zh) 地震数据处理方法、地层品质因子估算方法与装置
EA032186B1 (ru) Сейсмическая адаптивная фокусировка
CN105353405B (zh) 一种全波形反演方法和***
CN109143331B (zh) 地震子波提取方法
Wang et al. Seismic data denoising for complex structure using BM3D and local similarity
AU2013201072B2 (en) Method and apparatus for automated noise removal from seismic data
CN104391324A (zh) 依赖频率的avo反演前的地震道集动校拉伸校正预处理技术
KR101693310B1 (ko) 평면파를 이용한 이방성 매질에서의 전파형 역산방법
Zhu et al. A multi-channel approach for automatic microseismic event association using RANSAC-based arrival time event clustering (RATEC)
Faucher et al. Eigenvector models for solving the seismic inverse problem for the Helmholtz equation
Shao et al. Simultaneous inversion of Q and reflectivity using dictionary learning
CN109459788A (zh) 地层品质因子计算方法及***
CN109581506B (zh) 一种数据驱动的无拉伸动校正方法
CN108680957B (zh) 基于加权的局部互相关时频域相位反演方法
CN102353991B (zh) 基于匹配地震子波的物理小波的地震瞬时频率分析方法
CN115170428A (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