CN111427086B - 利用纵横波振幅频率分量预测地层最大衰减参数的方法 - Google Patents

利用纵横波振幅频率分量预测地层最大衰减参数的方法 Download PDF

Info

Publication number
CN111427086B
CN111427086B CN202010255392.1A CN202010255392A CN111427086B CN 111427086 B CN111427086 B CN 111427086B CN 202010255392 A CN202010255392 A CN 202010255392A CN 111427086 B CN111427086 B CN 111427086B
Authority
CN
China
Prior art keywords
wave
inelastic
impedance
longitudinal
transverse
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
CN202010255392.1A
Other languages
English (en)
Other versions
CN111427086A (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.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN202010255392.1A priority Critical patent/CN111427086B/zh
Publication of CN111427086A publication Critical patent/CN111427086A/zh
Application granted granted Critical
Publication of CN111427086B publication Critical patent/CN111427086B/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
    • 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

本发明涉及一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，包括：S1、以地层的最大衰减参数表征出随频率变化的纵波和横波反射系数近似表达式，以纵波和横波反射系数近似表达式为基础，建立纵波非弹性阻抗和横波非弹性阻抗的表达式；S2、将提取的纵波和横波的反射振幅频率分量，分别代入纵波非弹性阻抗和横波非弹性阻抗的表达式，得到预测的纵波非弹性阻抗和横波非弹性阻抗；S3、以初始模型为约束，利用纵波非弹性阻抗和横波非弹性阻抗联合反演得到地层最大衰减参数。与现有技术相比，本发明不仅可用于计算地层的最大衰减参数，进行地震波传播过程中衰减补偿，同时也可以用于识别含油气储层。

Description

利用纵横波振幅频率分量预测地层最大衰减参数的方法
技术领域
本发明涉及地球物理勘探技术领域，尤其是涉及一种利用纵横波振幅频率分量预测地层最大衰减参数的方法。
背景技术
当地震波在地下介质中传播的时候，地下介质的粘滞性和非均匀性将会引起地震波的能量被介质吸收，同时不同频率的地震波以不同的速度传播。地层非完全弹性引起的衰减已成为研究的重点，而且这种衰减通常用地层逆品质因子(1/Q)值来描述。
研究表明，地下介质的品质因子Q值与其内部结构(孔隙度、渗透率等)及充填流体性质(饱和度、粘度等)等密切相关。同时，在反映介质的渗透率及流体粘度性质上，地震波的衰减特性要比地震波速度变化特性更明显。因此，研究地震波的吸收衰减进而预测地层的衰减参数，不仅有利于地震数据的吸收衰减补偿，同时能够通过地层Q值识别油气藏位置，进而预测油气藏充填流体特征。
现有技术当前大都利用信号分解方法，得到各频率成分的能量随时间的变化。由于高频信号比低频信号衰减快，因此可以利用高频与低频信号变换的差异提取品质因子。预测品质因子的方法包括：谱比法，质心频率移动法以及基于地震波频谱变化的方法。譬如，Li等(Li Y,Zheng X,Zhang Y.High-frequency anomalies in carbonate reservoircharacterization using spectral decomposition.Geophysics,2011,76(3):V47-V57.)通过地震信号频谱分析预测碳酸盐岩的高频衰减特征。以上方法都是针对地震波传播过程中的变化来求取Q值，主要利用地震信号自身的时间和频率变化特征，并未将油气储层的岩石物理性质和地震信号变化结合起来，而且忽略了叠前地震振幅所蕴藏的随角度和频率同时变化的有效信息。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种利用纵横波振幅频率分量预测地层最大衰减参数的方法。
本发明的目的可以通过以下技术方案来实现：
一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，包括：
S1、以地层的最大衰减参数表征出随频率变化的纵波和横波反射系数近似表达式，以纵波和横波反射系数近似表达式为基础，建立纵波非弹性阻抗和横波非弹性阻抗的表达式；
S2、将提取的纵波和横波的反射振幅频率分量，分别代入纵波非弹性阻抗和横波非弹性阻抗的表达式，得到预测的纵波非弹性阻抗和横波非弹性阻抗；
S3、以初始模型为约束，利用纵波非弹性阻抗和横波非弹性阻抗联合反演得到地层最大衰减参数。
优选的，所述S1中纵波和横波反射系数近似表达式由粘弹性介质中的地震纵波传播速度和横波传播速度得到。
优选的，所述纵波和横波反射系数近似表达式包括：
纵波-纵波反射系数近似式：
Figure BDA0002437106920000021
纵波-转换横波反射系数近似式：
Figure BDA0002437106920000022
其中，f为频率，M为纵波模量，μ为横波模量，ρ为密度，QPm为地层最大衰减参数，θP为纵波入射角，与θP有关的各参数具体为：
Figure BDA0002437106920000023
aμ(θP)＝-2γsin2θP
Figure BDA0002437106920000024
Figure BDA0002437106920000025
Figure BDA0002437106920000026
Figure BDA0002437106920000027
Figure BDA0002437106920000028
其中，γ为横纵波模量比值，
Figure BDA0002437106920000031
Figure BDA0002437106920000032
其中，fr为参考频率。
优选的，所述纵波非弹性阻抗的表达式为：
Figure BDA0002437106920000033
所述横波非弹性阻抗的表达式为：
Figure BDA0002437106920000034
其中，t为时间，M(t)、μ(t)、ρ(t)分别为纵波模量、横波模量及密度的时间采样点。
优选的，所述S2具体包括：采集不同角度、不同频率的纵波和横波反射振幅，输入纵波非弹性阻抗和横波非弹性阻抗的表达式中，利用最小二乘反演算法得到预测的纵波非弹性阻抗和横波非弹性阻抗。
优选的，所述纵波非弹性阻抗、横波非弹性阻抗和待反演的地层最大衰减参数之间的关系为：
d＝G(m)
其中，
Figure BDA0002437106920000035
m＝[M μ ρ 1/QPm]T
其中，G(m)为与入射角度和频率有关的正演算子，
Figure BDA0002437106920000036
Figure BDA0002437106920000037
分别为纵波非弹性阻抗、横波非弹性阻抗，M为纵波模量，μ为横波模量，ρ为密度，QPm为地层最大衰减参数。
优选的，所述S3中通过牛顿迭代法求解所述d＝G(m)的未知数变量m得到地层最大衰减参数，m的解为：
m＝m0+ζΔm，
其中，m0为利用测井数据和岩石物理分析得到的纵波模量、横波模量、密度和最大衰减参数的初始模型，ζ为迭代步长，Δm为模型更新值。
优选的，所述模型更新值为：
Δm＝-(YYT)-1Y
其中，
Figure BDA0002437106920000041
dobs为利用不同频率不同角度地震数据估测的非弹性阻抗数据体，dmod为已知初始模型值通过
Figure BDA0002437106920000042
和
Figure BDA0002437106920000043
的表达式计算得到的纵波和横波非弹性阻抗模型值。
与现有技术相比，本发明首先提出了线性化的包含地层最大衰减参数值的频变的纵波和横波反射系数近似公式，进而开展利用纵波和横波反射振幅频率分量联合反演预测地层最大衰减值的方法，不仅可以充分利用多波多分量地震资料包含的丰富的有效信息，而且对地层最大衰减参数值的预测精度会有较好地提高，具有较高的鲁棒性，对油气储层的预测以及孔隙流体类型的识别具有十分重要的意义。
附图说明
图1为本发明的方法流程图；
图2为本发明应用于含噪声合成数据得到的反演值与真实值之间的对比，其中(a)为信噪比为5的情况，(b)为信噪比为2的情况；
图3为本发明应用于实际地震数据联合反演得到的结果，其中(a)为纵波模量M的结果，(b)为横波模量μ的结果，(c)为密度ρ的结果，(d)为底层最大衰减参数值
Figure BDA0002437106920000044
的结果；
图4为抽取图3反演结果当中的过井道反演结果和测井数据对比。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
本申请提出一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，为地震波传播过程模拟的衰减补偿提供品质因子Q值，同时为地下油气储层的预测和流体识别提供衰减参数依据。该方法包括：
S1、以地层的最大衰减参数表征出随频率变化的纵波和横波反射系数近似表达式，以纵波和横波反射系数近似表达式为基础，建立纵波非弹性阻抗和横波非弹性阻抗的表达式；
S2、将提取的纵波和横波的反射振幅频率分量，分别代入纵波非弹性阻抗和横波非弹性阻抗的表达式，得到预测的纵波非弹性阻抗和横波非弹性阻抗；
S3、以初始模型为约束，利用纵波非弹性阻抗和横波非弹性阻抗联合反演得到地层最大衰减参数。
步骤S1中建立了地震纵波和横波速度与最大衰减参数的关系，用于直接模拟衰减参数变化对频变地震波速度的影响。纵波和横波反射系数近似表达式由粘弹性介质中的地震纵波传播速度和横波传播速度得到。
粘弹性介质中的地震纵波传播速度
Figure BDA0002437106920000051
表示为：
Figure BDA0002437106920000052
其中，α为参考频率下的纵波速度，QP(f为频变的纵波衰减参数，f为频率，fr为参考频率。
频变的纵波衰减参数QP(f可以表示为最大衰减值的函数：
Figure BDA0002437106920000053
其中，QPm为地层最大衰减参数。
将公式(2)代入公式(1)可以改写纵波速度(舍去虚数项，仅保留实部)：
Figure BDA0002437106920000054
其中
Figure BDA0002437106920000055
类比于纵波速度的表征，接下来给出横波速度的表达式为：
Figure BDA0002437106920000056
其中，β为参考频率下的横波速度，
Figure BDA0002437106920000057
σ是泊松比。可以看出，横波速度已经利用纵波最大衰减参数进行表征。
基于已经表征的地震纵波和横波速度表达式，通过求解粘弹性介质的扩展的Zoeppritz方程求解，进而获得纵波和横波反射系数近似表达式，为模拟衰减对地震纵波和横波反射振幅影响奠定了理论基础，具体如下：
纵波-纵波反射系数近似式：
Figure BDA0002437106920000058
纵波-转换横波反射系数近似式：
Figure BDA0002437106920000059
其中，M为纵波模量，μ为横波模量，ρ为密度，θP为纵波入射角，与θP有关的各参数具体为：
Figure BDA0002437106920000061
aμ(θP)＝-2γsin2θP
Figure BDA0002437106920000062
Figure BDA0002437106920000063
Figure BDA0002437106920000064
Figure BDA0002437106920000065
Figure BDA0002437106920000066
其中，γ为横纵波模量比值，
Figure BDA0002437106920000067
Figure BDA0002437106920000068
其中，fr为参考频率。
基于已推导的纵波-纵波反射系数和纵波-转换横波反射系数近似式，可以得到对应的纵波非弹性阻抗
Figure BDA0002437106920000069
和横波非弹性阻抗
Figure BDA00024371069200000610
的表达式。
纵波非弹性阻抗的表达式为：
Figure BDA00024371069200000611
横波非弹性阻抗的表达式为：
Figure BDA00024371069200000612
其中，t为时间，M(t)、μ(t)、ρ(t)分别为纵波模量、横波模量及密度的时间采样点。
步骤S2具体包括：采集不同角度、不同频率的纵波和横波反射振幅，输入纵波非弹性阻抗和横波非弹性阻抗的表达式中，利用最小二乘反演算法得到预测的纵波非弹性阻抗和横波非弹性阻抗。得到的纵波非弹性阻抗、横波非弹性阻抗和待反演的地层最大衰减参数之间的关系为：
d＝G(m) (10)
其中，
Figure BDA0002437106920000071
m＝[M μ ρ 1/QPm]T (12)
其中，
Figure BDA0002437106920000072
Figure BDA0002437106920000073
分别为纵波非弹性阻抗、横波非弹性阻抗，M为纵波模量，μ为横波模量，ρ为密度，QPm为地层最大衰减参数，G(m)为与入射角和频率有关的正演算子。
步骤S3中通过牛顿迭代法求解d＝G(m)得到地层最大衰减参数，未知数变量m的解为：
m＝m0+ζΔm (13)
其中，m0为利用测井数据和岩石物理分析得到的纵波模量、横波模量、密度和最大衰减参数的初始模型，ζ为迭代步长，Δm为模型更新值，其表示为：
Δm＝-(YYT)-1Y (14)
其中，
Figure BDA0002437106920000074
dobs为利用不同频率不同角度地震数据估测的非弹性阻抗数据体，dmod为已知初始模型值通过
Figure BDA0002437106920000075
Figure BDA0002437106920000076
Figure BDA0002437106920000077
的表达式计算得到的纵波和横波非弹性阻抗模型值。
实施例
为了验证本方法在预测地层最大衰减参数值的可靠性，将本方法分别应用于合成地震数据和实际采集地震数据处理分析。
图2中(a)和(b)分别为将本方法应用于含随机噪声的合成地震数据(信噪比为5和2)时，反演得到的模型真实值和预测值的对比。可以看出，即使在信噪比为2时，该方法得到的反演结果和真实值之间仍能较好吻合。
图3为将该方法应用实际地震数据得到的纵波模量、横波模量、密度和最大衰减值的反演结果。图4为抽取反演结果当中的过井道和真实测井数据的对比。可以看出，利用本方法纵横波联合反演得到的参数值与测井信息之间能够较好匹配。
图2～图4的数据处理效果验证了本方法中纵横波联合反演方法预测最大衰减值的可行性与有效性，具有较大的应用潜力。

Claims (5)

1.一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，其特征在于，包括：
S1、以地层的最大衰减参数表征出随频率变化的纵波和横波反射系数近似表达式，以纵波和横波反射系数近似表达式为基础，建立纵波非弹性阻抗和横波非弹性阻抗的表达式，
S2、将提取的纵波和横波的反射振幅频率分量，分别代入纵波非弹性阻抗和横波非弹性阻抗的表达式，得到预测的纵波非弹性阻抗和横波非弹性阻抗，
S3、以初始模型为约束，利用纵波非弹性阻抗和横波非弹性阻抗联合反演得到地层最大衰减参数；
所述S1中纵波和横波反射系数近似表达式由粘弹性介质中的地震纵波传播速度和横波传播速度得到；
所述纵波和横波反射系数近似表达式包括：
纵波-纵波反射系数近似式：
Figure FDA0002864253400000011
纵波-转换横波反射系数近似式：
Figure FDA0002864253400000012
其中，f为频率，M为纵波模量，μ为横波模量，ρ为密度，QPm为地层最大衰减参数，θP为纵波入射角，与θP有关的各参数具体为：
Figure FDA0002864253400000013
aμ(θP)＝-2γsin2θP
Figure FDA0002864253400000014
Figure FDA0002864253400000015
Figure FDA0002864253400000016
Figure FDA0002864253400000017
Figure FDA0002864253400000021
其中，γ为横纵波模量比值，
Figure FDA0002864253400000022
Figure FDA0002864253400000023
其中，fr为参考频率；
所述纵波非弹性阻抗的表达式为：
Figure FDA0002864253400000024
所述横波非弹性阻抗的表达式为：
Figure FDA0002864253400000025
其中，t为时间，M(t)、μ(t)、ρ(t)分别为纵波模量、横波模量及密度的时间采样点。
2.根据权利要求1所述的一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，其特征在于，所述S2具体包括：采集不同角度、不同频率的纵波和横波反射振幅，输入纵波非弹性阻抗和横波非弹性阻抗的表达式中，利用最小二乘反演算法得到预测的纵波非弹性阻抗和横波非弹性阻抗。
3.根据权利要求1所述的一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，其特征在于，所述纵波非弹性阻抗、横波非弹性阻抗和待反演的地层最大衰减参数之间的关系为：
d＝G(m)
其中，
Figure FDA0002864253400000026
m＝[M μ ρ 1/QPm]T
其中，G(m)为与入射角度和频率有关的正演算子，
Figure FDA0002864253400000027
分别为纵波非弹性阻抗、横波非弹性阻抗，M为纵波模量，μ为横波模量，ρ为密度，QPm为地层最大衰减参数。
4.根据权利要求3所述的一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，其特征在于，所述S3中通过牛顿迭代法求解所述d＝G(m)的未知数变量m得到地层最大衰减参数，m的解为：
m＝m0+ζΔm
其中，m0为利用测井数据和岩石物理分析得到的纵波模量、横波模量、密度和最大衰减参数的初始模型，ζ为迭代步长，Δm为模型更新值。
5.根据权利要求4所述的一种利用纵横波振幅频率分量预测地层最大衰减参数的方法，其特征在于，所述模型更新值为：
Δm＝-(YYT)-1Y
其中，
Figure FDA0002864253400000031
dobs为利用不同频率不同角度地震数据估测的非弹性阻抗数据体，dmod为已知初始模型值通过
Figure FDA0002864253400000032
和
Figure FDA0002864253400000033
的表达式计算得到的纵波和横波非弹性阻抗模型值。
CN202010255392.1A 2020-04-02 2020-04-02 利用纵横波振幅频率分量预测地层最大衰减参数的方法 Active CN111427086B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010255392.1A CN111427086B (zh) 2020-04-02 2020-04-02 利用纵横波振幅频率分量预测地层最大衰减参数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010255392.1A CN111427086B (zh) 2020-04-02 2020-04-02 利用纵横波振幅频率分量预测地层最大衰减参数的方法

Publications (2)

Publication Number Publication Date
CN111427086A CN111427086A (zh) 2020-07-17
CN111427086B true CN111427086B (zh) 2021-03-26

Family

ID=71557417

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010255392.1A Active CN111427086B (zh) 2020-04-02 2020-04-02 利用纵横波振幅频率分量预测地层最大衰减参数的方法

Country Status (1)

Country Link
CN (1) CN111427086B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113064204B (zh) * 2021-03-16 2022-07-22 中油奥博(成都)科技有限公司 一种井中地震吸收参数提取及地层岩性检测方法
CN113219536B (zh) * 2021-06-25 2022-03-01 成都理工大学 依赖频率的纵横波衰减参数的叠前地震反演方法
CN115267893A (zh) * 2022-07-26 2022-11-01 同济大学 泥质含量、孔隙度及流体模量的非线性反演方法和系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102109612A (zh) * 2009-12-23 2011-06-29 中国石油天然气集团公司 一种地震波吸收衰减补偿方法
CN102288994A (zh) * 2011-07-01 2011-12-21 中国海洋石油总公司 Radon谱约束下高维地震数据规则化方法
CN106249294A (zh) * 2015-06-12 2016-12-21 中国石油化工股份有限公司 一种储层烃类检测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102650701B (zh) * 2012-04-27 2014-02-05 中国石油天然气股份有限公司 一种基于偶极子波进行储层波阻抗预测的方法
US9897710B2 (en) * 2012-06-26 2018-02-20 Schlumberger Technology Corporation Borehole seismic acquisition tools, systems and methods

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102109612A (zh) * 2009-12-23 2011-06-29 中国石油天然气集团公司 一种地震波吸收衰减补偿方法
CN102288994A (zh) * 2011-07-01 2011-12-21 中国海洋石油总公司 Radon谱约束下高维地震数据规则化方法
CN106249294A (zh) * 2015-06-12 2016-12-21 中国石油化工股份有限公司 一种储层烃类检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冀中坳陷深潜山及潜山内幕地震采集关键技术及应用;陈敬国，等;《中国石油勘探》;20191130;第24卷(第6期);807-814 *

Also Published As

Publication number Publication date
CN111427086A (zh) 2020-07-17

Similar Documents

Publication Publication Date Title
CN111427086B (zh) 利用纵横波振幅频率分量预测地层最大衰减参数的方法
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
Park et al. Multichannel analysis of surface waves
CN102466816B (zh) 一种叠前地震数据地层弹性常数参数反演的方法
CN104570072B (zh) 一种粘弹性介质中的球面pp波反射系数建模方法
US20090132169A1 (en) Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements
WO2014191011A1 (en) High resolution estimation of attenuation from vertical seismic profiles
GB2382651A (en) Absolute preserved amplitude processing of seismic data
CN107065013B (zh) 一种地震尺度下的层速度确定方法及装置
CN103487831B (zh) Avo地震正演计算方法
Kamath et al. Elastic full-waveform inversion for VTI media: A synthetic parameterization study
Germán Rubino et al. Seismic characterization of thin beds containing patchy carbon dioxide-brine distributions: A study based on numerical simulations
Stovas et al. AVO attribute inversion for finely layered reservoirs
Takam Takougang et al. Imaging high-resolution velocity and attenuation structures from walkaway vertical seismic profile data in a carbonate reservoir using visco-acoustic waveform inversion
CN115061202B (zh) 一种页岩含气性地震储层直接检测方法
CN106249294A (zh) 一种储层烃类检测方法
CN116500679A (zh) 一种井震标定方法、装置和相关设备
CN115201915A (zh) 海洋宽方位数据深层甜点预测方法、存储介质和设备
CN109239777B (zh) 一种利用联合反演方法检测构造煤发育的方法
CN112558151A (zh) 一种衰减介质中ava波形模拟方法及系统
Schwenk Constrained parameterization of the multichannel analysis of surface waves approach with application at Yuma Proving Ground, Arizona
CN113806674A (zh) 古河道纵向尺度的量化方法、装置、电子设备及存储介质
Hu et al. Estimation of rock physics properties via full waveform inversion of vertical seismic profile data recorded by accelerometer and fibre-optic sensors
Li et al. Seismic attenuation and velocities of P-and S-waves in the German KTB area
Wang et al. Effects of the outermost boundary on acoustic waves in an artificial cased borehole

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