CN109459788A - 地层品质因子计算方法及系统 - Google Patents
地层品质因子计算方法及系统 Download PDFInfo
- 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
- power spectrum
- quality factor
- formation quality
- attenuation
- 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
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 32
- 238000001228 spectrum Methods 0.000 claims abstract description 94
- 230000015572 biosynthetic process Effects 0.000 claims description 53
- 230000002238 attenuated effect Effects 0.000 claims description 15
- 230000005540 biological transmission Effects 0.000 claims description 12
- 238000000034 method Methods 0.000 abstract description 35
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 12
- 238000010521 absorption reaction Methods 0.000 abstract description 7
- 238000001514 detection method Methods 0.000 abstract description 6
- 239000003345 natural gas Substances 0.000 abstract description 6
- 230000003595 spectral effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 5
- 230000036039 immunity Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 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/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所述的地层品质因子计算系统,其中,所述地层品质因子为:
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117724165A (zh) * | 2024-02-07 | 2024-03-19 | 中国石油大学(华东) | 一种基于时变子波的品质因子估计方法 |
Citations (10)
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 | 中国石油化工股份有限公司 | 一种高精度的品质因子提取方法 |
-
2017
- 2017-09-06 CN CN201710796553.6A patent/CN109459788B/zh active Active
Patent Citations (10)
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)
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值反演", 《石油地球物理勘探》 * |
杨登锋,等: "能量谱质心频移法Q值估计", 《石油地球物理勘探》 * |
Cited By (2)
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 |
---|---|---|
US6556922B2 (en) | Non-stationary filter design and application for seismic signals | |
Sun et al. | Cross-correlation analysis and time delay estimation of a homologous micro-seismic signal based on the Hilbert–Huang transform | |
CN108535775B (zh) | 非平稳地震资料声波阻抗反演方法及装置 | |
CN107132579B (zh) | 一种保地层结构的地震波衰减补偿方法 | |
CN112835103B (zh) | 自适应去鬼波与宽频准零相位反褶积联合处理方法及系统 | |
WO2013090713A2 (en) | Iterative dip-steering median filter for seismic data processing | |
CN110646841B (zh) | 时变稀疏反褶积方法及系统 | |
CN110658558A (zh) | 吸收衰减介质叠前深度逆时偏移成像方法及系统 | |
CN109459788B (zh) | 地层品质因子计算方法及系统 | |
CN110487916A (zh) | 一种适用于换热污垢超声检测回波信号的降噪方法 | |
CN117574062A (zh) | 基于vmd-dnn模型的小回线瞬变电磁信号去噪方法 | |
Oboué et al. | An advanced median filter for improving the signal-to-noise ratio of seismological datasets | |
Atefi et al. | Rapid estimation of earthquake magnitude by a new wavelet‐based proxy | |
CN110837119B (zh) | 一种增强反q值补偿稳定性的地震资料处理方法及系统 | |
CN105866838B (zh) | 一种地震数据低频信息补偿方法及装置 | |
Zhang et al. | Q estimation via the discriminant method based on error modeling | |
CN113341463A (zh) | 一种叠前地震资料非平稳盲反褶积方法及相关组件 | |
KR101355107B1 (ko) | 잡음 제거를 통한 지하 매질구조 추정방법 및 그 장치 | |
CN109212590B (zh) | 一种保幅的逆时偏移低频噪音压制方法及系统 | |
CN118091753B (zh) | 基于稀疏反褶积与时频分析的强干扰环境震源初至拾取方法、装置、设备及介质 | |
Angelopoulos et al. | Nonparametric spectral estimation-an overview | |
Wang et al. | A common-reflection-point gather random noise attenuation method based on the synchrosqueezing wavelet transform | |
CN117890979B (zh) | 地震数据自适应弱反射信号补偿方法、系统、设备及介质 | |
CN110794455B (zh) | 一种地震波传播能量衰减补偿方法 | |
CN113917526B (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 |