CN113534255A - 一种任意形态不连续面自适应表达的方法 - Google Patents
一种任意形态不连续面自适应表达的方法 Download PDFInfo
- Publication number
- CN113534255A CN113534255A CN202110768611.0A CN202110768611A CN113534255A CN 113534255 A CN113534255 A CN 113534255A CN 202110768611 A CN202110768611 A CN 202110768611A CN 113534255 A CN113534255 A CN 113534255A
- Authority
- CN
- China
- Prior art keywords
- free surface
- model
- constitutive relation
- parameters
- free
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 239000006185 dispersion Substances 0.000 claims abstract description 5
- 239000007787 solid Substances 0.000 claims description 27
- 239000007788 liquid Substances 0.000 claims description 19
- 238000012986 modification Methods 0.000 claims description 10
- 230000004048 modification Effects 0.000 claims description 10
- 238000012937 correction Methods 0.000 claims description 8
- 230000003044 adaptive effect Effects 0.000 claims description 5
- 230000001788 irregular Effects 0.000 claims description 4
- 238000004088 simulation Methods 0.000 abstract description 7
- 238000006073 displacement reaction Methods 0.000 abstract description 3
- 230000002349 favourable effect Effects 0.000 abstract 1
- 238000004364 calculation method Methods 0.000 description 12
- 239000011159 matrix material Substances 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 241001344923 Aulorhynchidae Species 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000011438 discrete method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
-
- 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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G01V1/01—
-
- 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/626—Physical property of subsurface with anisotropy
-
- 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/64—Geostructures, e.g. in 3D data cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
- G01V2210/641—Continuity of geobodies
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Abstract
本发明公开了一种任意形态不连续面自适应表达的方法,包括如下步骤:导入正演的模型,包括纵横波速度以及密度,对于各向异性介质,导入各向异性参数,并根据模型参数设置空间步长和时间步长;对自由表面进行阶梯化离散;对自由表面上不同点进行本构关系修正;当模型地表起伏时,利用修正后的本构关系更改模型的第一层参数,这样在不对原始的有限差分法代码进行修改的情况下就可以引入自由表面的影响;将修正后的本构关系带入到位移应力方程中经过一系列的运算可在各异性介质的情况下引入自由表面的影响。本发明所述方法可以准确地对不连续面附近的波场进行数值模拟,有利于对此类地震数据的信息进行提取和分析。
Description
技术领域
本发明属于地震勘探技术领域,尤其涉及一种任意形态不连续面自适应表达的方法。
背景技术
在地震学中,通常把地震传播过程中突然发生变化的地方成为不连续面。根据不连续面两侧介质物理的物理性质(固体,液体和气体),可以将不连续面划分为以下四类:1)固体/固体不连续面;2)固体/液体不连续面;3)固体自由表面(陆地表面),也称固体/真空不连续面;4)液体自由表面(海平面),也称液体/真空不连续面。这类不连续面附近均会发生复杂的波转换现象和强能量的界面波生成。例如,固体/液体不连续面往往用于描述海底界面或者地球内外核边界,在其附近会产生一种沿该界面传播的Scholte wave,以及发生P-S波转换,会产生强振幅面波(Rayleigh wave或love wave),其频散特征在近地表面地震勘探中有重要的意义和用途;液体自由表面用于描述海洋表面,是海洋勘探中多次波产生的原因,其可以被视为噪音进行除去(多次波去噪技术),也可以用于照明角度和区域补偿(多次波成像技术)。
在复杂的介质地震波数值模拟中,往往会遇到以上不连续面的组合,对于这些不连续面的正确数值建模表达会直接影响数值模拟的准确性和精度,特别是对于地表观测地震数据和海底地震仪观测地震数据,其检波器和地震仪直接位于不连续面上或附附近,因此受不连续面的影响更大。能否正确的对不连续面附近的波长进行数值模拟会直接对这类地震数据中提取和分析信息产生影响。
此外,这类不连续面还会破坏一些数值离散方法的完整性,例如有限差分法,伪谱法等基于偏微分方程强解形式的算法。目前工业界和科研界最常用的地震数值模拟技术是有限差分法,其适用范围为参数均匀或光滑变化介质,直接应用于含强阻抗对比的不连续面时会导致计算不稳定,数值假象产生和计算精度降低,进而影响最终模拟结果的准确性。
目前,一般采用直接离散化边界条件,以差分形式来近似边界条件中的控件导数,以固体自由表面的实现为例,有中心差分显示近似格式(Alterman和Karal,1968)、单边显示近似格式(Vidale和Clayton,1986)以及复合近似等格式(Ilan等,1975;Ilan和Loewenthal,1976;Lan和Zhang,2011);另外,在一阶速度-应力方程表达形式下对应力张量和速度矢量分别操作,以固体自由表面的实现为例,有真空格式(Zahradnik等,1993;Bohlen和Saenger,2006)、应力镜像法(Levander,1988;Graves,1996;Kristek等,2002)和一些介质平均方法(Mittet,2002;Xu等,2007);更有,基于Mimetic算子的离散方法,例如dela Puente等(2014),Shragge和Tapley(2017),Konuk和Shragge(2020),Shthi等(2021);除以上三种方法,还有有限元类方法,包含谱yuan法(Komatitsch和Tromp,1999;Komatitsch等,1999;Tromp等,2008;Peter等,2011),这类基于偏微分弱解形式的数值模拟方法,可参考现有的普元发开源程序包:SPECFEM2D和SPECFEM3D。
但是,现有技术存在以下问题:1)应用于高泊松比介质的正演模拟会产生计算不稳定的现象,并且最高只能到达二阶的计算精度。在截面起伏情况下其计算稳定性还会受到截面最大曲率的影响,曲率越大越容易不稳定;2)计算精度低,且受泊松比变化影响;在截面起伏情况下,截面两侧物理场所需的关于界面的对称操作变得及其复杂;3)算法实现复杂,需额外开辟内存控件来存储集合拓扑变量,且计算稳定性还会收到截面最大曲率的影响,曲率越大越容易不稳定;4)相比较于基于有限差分法的数值模拟技术,算法理论和实现复杂,计算量大,内存消耗费多。
发明内容
针对上述背景技术中提到的技术缺陷,为了解决上述题述问题,本发明提供了一种任意形态不连续面自适应表达的方法,其特征在于,包括如下步骤:
S1:导入初始正演模型,包括纵横波的速度和密度,设置震源子波的类型和主频的参数,并根据各模型的参数设置空间步长和时间步长;
S2:根据模型参数选取时间步长和空间步长,设置差分的阶数和震源子波;
S3:当模型表面起伏时,对模型的自由表面进行阶梯化离散,将非规则的起伏地表变成规则的阶梯状网格;
S4:运用介质参数修正法对自由表面上非水平的点进行本构关系和密度的修正;
S5:利用修正后的本构关系来替换自由表面处原来的本构关系,在不对原始的有限差分法代码进行修改的情况下,即可引入自由表面的影响;
S6:对于各向异性介质,将修正后的本构关系带入到有限差分方法的代码中,即可在各向异性介质的情况下引入自由表面的影响。
进一步的,所述自由表面包括液体自由表面和固体自由表面。
进一步的,所述不连续面包括固体/液体不连续面和液体/真空不连续面。
进一步的,对于各异性介质,需要导入各向异性参数。
进一步的,所述异性介质适用于固体自由表面部分。
本发明所述通过介质参数修正实现任意形态不连续面自适应表达的方法,可直接继承与应用于现有的工业界和学术界主流的速度-应力交错网格有限分数值程序,精确贴合不连续面的空间位置,没有传统方法中的版网格偏移误差;通过简单的修正不连续面附近网格点的介质参数来实现,且该修正过程只需一次,计算高效、操作简单;可应用高阶有限差分算子,实现与基于偏微分方程弱解类精确实现方法相当的计算精度;适用于各向异性介质附近的不连续面,且具有泊松比自适应的特点。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面对实施例描述中所需要使用的附图作一个简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例的介质参数修正实现任意形态不连续面自适应表达的方法的流程图;
图2为本发明实施例的自由表面离散化示意图;
图3为本发明实施例的固体与空气的位置关系示意图,图a)指空气在自由表面左侧,图b)指空气在自由表面前面,图c)空气在自由表面右侧,图d)指空气在自由表面后面,图e)指空气在自由表面上面。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提供了一种任意形态不连续面自适应表达的方法,其特征在于,包括如下步骤:
S1:导入初始正演模型,包括纵横波的速度和密度,设置震源子波的类型和主频的参数,并根据各模型的参数设置空间步长和时间步长;
S2:根据模型参数选取时间步长和空间步长,设置差分的阶数和震源子波;
S3:当模型表面起伏时,对模型的自由表面进行阶梯化离散,将非规则的起伏地表变成规则的阶梯状网格;
S4:运用介质参数修正法对自由表面上非水平的点进行本构关系和密度的修正;
S5:利用修正后的本构关系来替换自由表面处原来的本构关系,在不对原始的有限差分法代码进行修改的情况下,即可引入自由表面的影响,所述影响包括但不限于自由表面引起的地震波的反射;
S6:对于各向异性介质,将修正后的本构关系带入到有限差分方法的代码中,即可在各向异性介质的情况下引入自由表面的影响。
进一步的,所述自由表面包括液体自由表面和固体自由表面;所述不连续面包括固体/液体不连续面和液体/真空不连续面;所述异性介质适用于固体自由表面部分。
具体地,对固体/液体不连续面假设固体介质1对应的本构关系为固体介质2对应的本构关系为根据平均介质思想,可将两个连接的介质等效成一个等效平均介质,其本构关系为同理,对于固体/液体不连续面,假设液体介质1对应的本构关系为液体介质2对应的本构关系为根据平均介质思想,可将两个连接的介质等效成一个等效平均介质,其本构关系为
上述式中:τ指应力张量矩阵;M1指弹性介质1,M2指弹性介质2,A代表平均介质;E指等效的弹性系数矩阵;λ1、λ2和μ1、μ2是拉梅常数;ε为应变张量矩阵。
上述公式中c11,c13,c33,c44为各向异性介质的弹性参数。
上述式中:τ指应力张量矩阵;C代表弹性系数;M1指弹性介质1,M2指弹性介质2,A代表平均介质;E指等效的弹性系数矩阵;λ1、λ2和μ1、μ2是拉梅常数;ε为应变张量矩阵。
具体的,在步骤S1中,对于各向异性介质,需要导入各向异性参数c11,c13,c33,c44。同时需要确定差分的阶数、震源类型和震源主频等参数,根据模型参数设置空间步长和时间步长。
进一步的,当模型的地表面起伏时,需要对模型的自由面进行阶梯化离散,离散如图1所述,将非规则的欺负地面变成规则的阶梯状网格,如图2,图2中0L、0R指网格点的外角点,IL、IR代表网格点的内角点,图2显示了起伏地表的实现方式;对自由表面网格点进行划分,并对不同网格点根据表1进行本构关系修正。
在步骤S4中,运用介质参数法对自由表面上非水平的点进行本构关系修正,修正关系如下表所示,其中VL代表自由表面左边为空气的情况,VR代表自由表面右边为空气的情况,VB代表自由表面后边为空气的情况。H代表自由表面上方为空气的情况,具体可参考图3,图3中A指空气,S指地下介质。
表1:起伏地表自由表面处本构关系
上述表中:τij指应力张量的第ij分量;εij指应变张量的第ij分量;λ和μ指拉梅常数;ρ指介质的密度;α=(2μ(λ+μ)/(λ+2μ)),β=(μλ)/(λ+2μ)。
进一步的,利用修正后的本构关系,只需更改模型的第一层参数,即自由表面处的参数,这样在不对原始的有限差分法代码进行修改的情况下就可以引入自由表面的影响。在自由表面处的本构关系如下所示:
τzz|z=0=0,
τxy|z=0=τyx|z=0=μεyx=μεxy,
τyz|z=0=τzy|z=0=0,
τxz|z=0=τzx|z=0=0,
ui是位移的第i个分量,τij指应力张量的第ij分量;εij指应变张量的第ij分量;λ和μ指拉梅常数;ρ指介质的密度。
进一步的,对于各向异性介质,以VTI介质为例,VTI为一种各向异性介质,其修正后的本构关系如下式所示:
修正后的本构关系带入到位移应力方程中可得:
由此即可在各向异性介质的情况下引入自由表面的影响。
本发明所述通过介质参数修正实现任意形态不连续面自适应表达的方法,可直接继承与应用于现有的工业界和学术界主流的速度-应力交错网格有限分数值程序,精确贴合不连续面的空间位置,没有传统方法中的版网格偏移误差;通过简单的修正不连续面附近网格点的介质参数来实现,且该修正过程只需一次,计算高效、操作简单;可应用高阶有限差分算子,实现与基于偏微分方程弱解类精确实现方法相当的计算精度;适用于各向异性介质附近的不连续面,且具有泊松比自适应的特点。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例和落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (5)
1.一种任意形态不连续面自适应表达的方法,其特征在于,包括如下步骤:
S1:导入初始正演模型,包括纵横波的速度和密度,设置震源子波的类型和主频的参数,并根据各模型的参数设置空间步长和时间步长;
S2:根据模型参数选取时间步长和空间步长,设置差分的阶数和震源子波;
S3:当模型表面起伏时,对模型的自由表面进行阶梯化离散,将非规则的起伏地表变成规则的阶梯状网格;
S4:运用介质参数修正法对自由表面上非水平的点进行本构关系和密度的修正;
S5:利用修正后的本构关系来替换自由表面处原来的本构关系,在不对原始的有限差分法代码进行修改的情况下,即可引入自由表面的影响;
S6:对于各向异性介质,将修正后的本构关系带入到有限差分方法的代码中,即可在各向异性介质的情况下引入自由表面的影响。
2.根据权利要求1所述通过介质参数修正实现任意形态不连续面自适应表达的方法,其特征在于:所述自由表面包括液体自由表面和固体自由表面。
3.根据权利要求1所述通过介质参数修正实现任意形态不连续面自适应表达的方法,其特征在于:所述不连续面包括固体/液体不连续面和液体/真空不连续面。
4.根据权利要求1所述通过介质参数修正实现任意形态不连续面自适应表达的方法,其特征在于:对于各异性介质,需要导入各向异性参数。
5.根据权利要求1所述通过介质参数修正实现任意形态不连续面自适应表达的方法,其特征在于:所述异性介质适用于固体自由表面部分。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110768611.0A CN113534255A (zh) | 2021-07-07 | 2021-07-07 | 一种任意形态不连续面自适应表达的方法 |
US17/393,317 US20230020158A1 (en) | 2021-07-07 | 2021-08-03 | Media parameter-modified method for realizing an adaptive expression of an arbitrary discontinuous surface |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110768611.0A CN113534255A (zh) | 2021-07-07 | 2021-07-07 | 一种任意形态不连续面自适应表达的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113534255A true CN113534255A (zh) | 2021-10-22 |
Family
ID=78097967
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110768611.0A Pending CN113534255A (zh) | 2021-07-07 | 2021-07-07 | 一种任意形态不连续面自适应表达的方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20230020158A1 (zh) |
CN (1) | CN113534255A (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749643A (zh) * | 2011-04-22 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种波动方程正演的瑞利面波频散响应计算方法及其装置 |
CN103562749A (zh) * | 2011-05-23 | 2014-02-05 | 帝国创新有限公司 | 可用于金属矿物和/或金刚石矿床的勘探、矿井设计、评估和/或开采的方法 |
-
2021
- 2021-07-07 CN CN202110768611.0A patent/CN113534255A/zh active Pending
- 2021-08-03 US US17/393,317 patent/US20230020158A1/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749643A (zh) * | 2011-04-22 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种波动方程正演的瑞利面波频散响应计算方法及其装置 |
CN103562749A (zh) * | 2011-05-23 | 2014-02-05 | 帝国创新有限公司 | 可用于金属矿物和/或金刚石矿床的勘探、矿井设计、评估和/或开采的方法 |
Non-Patent Citations (1)
Title |
---|
曹健: "弹性波有限差分模拟中自由表面的自适应表达", 《石油物探》, 25 July 2018 (2018-07-25), pages 1 - 3 * |
Also Published As
Publication number | Publication date |
---|---|
US20230020158A1 (en) | 2023-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Downton et al. | Azimuthal simultaneous elastic inversion for fracture detection | |
US10459117B2 (en) | Extended subspace method for cross-talk mitigation in multi-parameter inversion | |
CA2920008C (en) | Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media | |
US7675818B2 (en) | Method for tomographic inversion by matrix transformation | |
US20190018155A1 (en) | Visco-Pseudo-Elastic TTI FWI/RTM Formulation and Implementation | |
CN109143351B (zh) | 叠前各向异性特征参数反演方法及计算机可读存储介质 | |
Hu | An improved immersed boundary finite-difference method for seismic wave propagation modeling with arbitrary surface topography | |
Guitton et al. | A preconditioning scheme for full waveform inversion | |
Fu et al. | Multiscale phase inversion for 3D ocean‐bottom cable data | |
CA2615125A1 (en) | Method for predicting the best and worst in a set of non-unique solutions | |
NO20190489A1 (en) | Seismic modeling | |
Zhang et al. | Attenuation compensation for wavefield‐separation‐based least‐squares reverse time migration in viscoelastic media | |
US8437219B2 (en) | Correcting an acoustic simulation for elastic effects | |
CN113534255A (zh) | 一种任意形态不连续面自适应表达的方法 | |
Jiang et al. | Full waveform inversion based on inversion network reparameterized velocity | |
CN114895351A (zh) | 模拟地震波在任意不连续界面处传播的介质模型化的方法及装置 | |
AlSalem et al. | Embedded boundary methods for modeling 3D finite-difference Laplace-Fourier domain acoustic-wave equation with free-surface topography | |
WO2022069441A1 (en) | Method and system for processing seismic images to obtain seismic horizon surfaces for a geological formation | |
Kim et al. | Least-squares reverse time migration using analytic-signal-based wavefield decomposition | |
Chen* | Deblending by iterative orthogonalization and seislet thresholding | |
CN112099079B (zh) | 自适应分频串联反射率反演方法及系统 | |
Santos et al. | Evaluation of L-Curve and Theta-Curve approaches for the selection of regularization parameters in anisotropic traveltime tomography | |
Nikitin et al. | Simulation of seismic responses from the 3D non-linear model of the Bazhenov formation | |
Li et al. | 3D anisotropic full-waveform inversion for complex salt provinces | |
CN112684505B (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 |