CN111505718B - 一种高分辨率地下结构保幅成像方法 - Google Patents
一种高分辨率地下结构保幅成像方法 Download PDFInfo
- Publication number
- CN111505718B CN111505718B CN202010352142.XA CN202010352142A CN111505718B CN 111505718 B CN111505718 B CN 111505718B CN 202010352142 A CN202010352142 A CN 202010352142A CN 111505718 B CN111505718 B CN 111505718B
- Authority
- CN
- China
- Prior art keywords
- imaging
- imaging result
- result
- curvelet
- migration
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 98
- 238000013508 migration Methods 0.000 claims abstract description 31
- 230000005012 migration Effects 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 28
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 238000004088 simulation Methods 0.000 claims description 8
- 230000000644 propagated effect Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000005286 illumination Methods 0.000 description 2
- 238000002310 reflectometry Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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. for interpretation or for event detection
- G01V1/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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
- G01V1/364—Seismic filtering
-
- 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
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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
-
- 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/21—Frequency-domain filtering, e.g. band pass
-
- 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/51—Migration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/70—Other details related to processing
- G01V2210/74—Visualisation of seismic data
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
本发明公开了一种高分辨率地下结构保幅成像方法,首先由逆时偏移得到初始成像结果,在初始成像结果的基础上进行Born正演得到仿真地震数据,将仿真数据逆时偏移得到再次偏移成像的结果;然后对两次成像结果进行曲波变换,在曲波域逐点估计,将匹配两组曲波系数的维纳解作为匹配滤波器的解;最后将估计得到的匹配滤波器作用于初始成像结果得到高分辨率的保幅成像结果。该方法针对由于逆时偏移成像算子实际为正演算子的伴随算子而导致的成像结果模糊、振幅不均衡的问题,在初始偏移成像的结果上进行再次偏移成像,利用曲波变换在变换域寻找两次偏移成像结果的匹配滤波器作为逆Hessian算子的近似,将匹配滤波器作用于初始成像结果以达到改善成像质量的效果。
Description
技术领域
本发明属于地震勘探技术领域,涉及一种地下结构成像方法,具体涉及一种高分辨率地下结构保幅成像方法。
背景技术
在地震勘探技术领域,逆时偏移(RTM)是当前最先进的地下结构成像技术。逆时偏移采用双程波波动方程描述地震波在地下的传播规律,能够利用各种类型的波进行成像,如棱柱波、回转波等,可以对任意结构成像,无倾角限制,比如倒挂体的成像。然而,逆时偏移仍然是正演算子的伴随算子,其只能生成一个地下反射率模型的近似。结果是,逆时偏移生成的地震成像剖面是模糊的,并且具有不平衡的振幅。为了得到真振幅的成像结果,可以使用最小二乘方法得到正演算子的逆算子的近似。
最小二乘逆时偏移(LSRTM)的目标是找到最匹配观测数据的地下反射率模型。也就是说,它通过最小化由一个Born正演算子得到的仿真数据和真实观测数据之间的差异以得到一个最优的成像剖面。LSRTM通常以一种迭代的方式执行,计算量庞大,对于规模较大的数据集来说计算成本十分昂贵。为了减少计算成本,可以通过近似最小二乘泛函的逆Hessian算子的方法改善成像质量。在地震波的传播过程中,地震子波是带限和时变的,地下照明不均衡。最小二乘泛函的Hessian算子可用于描述这些波传播的效应。将Hessian算子的逆算子作用于原始成像结果即可校正这些波传播的效应,达到改善成像质量的目的。
然而,Hessian算子的规模庞大,难以在实际应用中得到,更不要提求得Hessian算子的逆算子。如何估计得到Hessian算子的逆算子是一个热点关注且未得到妥善解决的问题。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种高分辨率地下结构保幅成像方法,能够提高地震成像结果的分辨率,并且使成像结果拥有更加均衡的幅度,且具有较高的计算效率,从而能够针对性地消除地震成像剖面中模糊和振幅不均衡的效应。
为了达到上述目的,本发明采用以下技术方案予以实现:
本发明公开了一种高分辨率地下结构保幅成像方法,包括:
首先,由逆时偏移得到初始成像结果,在初始成像结果的基础上进行Born正演得到仿真地震数据,将仿真数据逆时偏移得到再次偏移成像的结果;
然后,对两次成像结果进行曲波变换,在曲波域逐点估计,将匹配两组曲波系数的维纳解作为匹配滤波器的解;
最后,将估计得到的匹配滤波器作用于初始成像结果得到高分辨率的保幅成像结果。
优选地,上述的高分辨率地下结构保幅成像方法,包括以下步骤:
1)将观测数据进行逆时偏移,得到初始成像结果m0(x);
2)将初始成像结果m0(x)带入一个Born正演算子得到仿真地震数据;
3)将仿真数据进行逆时偏移,得到再次成像的结果m'(x);
优选地,步骤1)中,将观测数据进行逆时偏移采用有限差分方法求解波动方程,其包含三个步骤,即震源波场的正向传播、检波器波场的反向传播和应用一个互相关成像条件。
进一步优选地,震源波场沿时间正向传播,检波器波场从最大记录时刻开始沿时间反向传播;观测数据是以真实速度模型作为输入,通过波动方程的有限差分正演模拟得到的数据。
优选地,步骤2)中,Born正演是通过已经得到的地震成像剖面得到仿真地震数据的过程,具体包括:
2.1)震源波场的正向传播得到背景波场ps(x),其中,背景波场是以偏移速度模型作为输入通过波动方程的有限差分正演模拟得到;
2.3)在检波点位置记录不同时刻的散射波场得到仿真地震数据。
与现有技术相比,本发明具有以下有益效果:
(1)本发明能够改善成像结果中由于震源子波频带有限引入的模糊效应,达到提高分辨率的效果。
(2)本发明能校正成像结果中反射体的振幅,改善波在地下传播时不同位置处照明不均衡导致的振幅不均衡的问题。
(3)与基于迭代算法的LSRTM相比,本发明的计算效率显著提高,本发明生成的成像结果与传统LSRTM迭代十次后的成像结果质量相当。然而,本发明只需要执行两次偏移成像和一次Born正演,而传统LSRTM迭代十次需要10次RTM偏移成像和9次Born正演,曲波变换的计算量与一次波传播过程的计算量相比微不足道。
附图说明
图1是本发明用于测试的真实速度模型;
图2是本发明用于测试的偏移速度模型;
图3是观测地震数据;
图4是通过Born正演得到的仿真地震数据;
图5是由观测地震数据反演得到的RTM初始成像结果;
图6是由仿真地震数据反演得到的RTM再次成像结果;
图7是由本发明方法得到的改善后的成像结果;
图8是由传统LSRTM迭代10此后得到的成像结果;
图9是真实反射系数、RTM结果和本发明结果的单道对比;
图10是RTM结果、LSRTM结果、本发明结果的波数谱比较。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
下面结合附图对本发明做进一步详细描述:
步骤1、将观测数据进行逆时偏移得到初始成像结果m0(x)
使用有限差分数值方法求解波动方程实现波动方程的模拟。RTM的成像过程可以描述为震源波场的正向传播,检波器波场的反向传播,和应用一个互相关成像条件。对于一个观测记录D(xr;t;xs),其震源位于xs=(xs,ys,zs=0),检波器位于xr=(xr,yr,zr=0)。首先,震源波场的正向传播可以描述为:
其中v(x)表示速度模型,x表示地下位置坐标,ps是震源波场,f(t)是震源子波。本发明用于测试的真实速度模型和偏移速度模型分别如图1和图2所示。
其次,检波器波场的反向传播可以描述为:
其中,pr表示检波器波场。
最后应用互相关成像条件得到地下结构的像:
其中,m(x)表示叠加成像剖面,并且Tmax是地震记录的总接收时间。将图3所示的观测数据用于RTM偏移成像得到的成像结果如图5所示。
步骤2、通过Born正演得到仿真地震数据。首先,通过震源波场的正向传播得到背景波场ps,描述为:
最后,记录检波点位置的散射波场即得到仿真地震数据:
d(xr;t;xs)=pr(xr;t;xs) (6)
结果如图4所示。
步骤3、将仿真数据再次逆时偏移成像得到再次偏移成像结果m'(x)
参见图6,具体实施方法与1中相同,只需要将观测数据D(xr;t;xs)替换为仿真数据d(xr;t;xs)即可。
和
其中,*表示复共轭。
步骤5、求解曲波域的匹配滤波器作为逆Hessian算子的近似。
其中F表示匹配滤波器。这个均方误差最小化问题的维纳解可以通过下面的逐点估计得到:
其中,k表示空间频率域的波数向量,a表示尺度,θ表示角度,ε是一个正则化参数,以保证除法的稳定性。
其中Ψ-1表示逆曲波变换,F表示在步骤5中估计得到的的匹配滤波器,表示初始成像结果m0(x)的曲波变换系数,表示本方法最终得到的高分辨率成像结果。最终得到的成像结果如图7所示,与图8所示的传统LSRTM迭代10次后的成像结果质量相当,分辨率得以提高,并且幅度更加均衡。
本发明的方法只需两次RTM偏移成像和一次Born正演,而传统LSRTM迭代十次需要10次RTM偏移成像和9次Born正演,所以本发明的方法大大提高了计算效率。其中,曲波变换及其逆变换的计算量与一次波传播过程的计算量相比可以忽略不计。图9中的单道比较说明了本发明方法得到的成像结果与RTM得到的结果相比,其与真实反射系数更加接近,保持了反射系数相对幅度关系。图10中的波数谱比较说明了本发明方法与LSRTM都具有通过恢复高波数分量以提高地震图像分辨率的的效果。
综上所述,本发明方法针对由于逆时偏移成像算子实际为正演算子的伴随算子而导致的成像结果模糊、振幅不均衡的问题,在初始偏移成像的结果上进行再次偏移成像,利用曲波变换在变换域寻找两次偏移成像结果的匹配滤波器作为逆Hessian算子的近似,将匹配滤波器作用于初始成像结果以达到改善成像质量的效果。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。
Claims (1)
1.一种高分辨率地下结构保幅成像方法,其特征在于,包括:
首先,由逆时偏移得到初始成像结果,在初始成像结果的基础上进行Born正演得到仿真地震数据,将仿真数据逆时偏移得到再次偏移成像的结果;
然后,对两次成像结果进行曲波变换,在曲波域逐点估计,将匹配两组曲波系数的维纳解作为匹配滤波器的解;
最后,将估计得到的匹配滤波器作用于初始成像结果得到高分辨率的保幅成像结果;
包括以下步骤:
1)将观测数据进行逆时偏移,得到初始成像结果m0(x);
2)将初始成像结果m0(x)带入一个Born正演算子得到仿真地震数据;
3)将仿真数据进行逆时偏移,得到再次成像的结果m'(x);
步骤1)中,将观测数据进行逆时偏移采用有限差分方法求解波动方程,其包含三个步骤,即震源波场的正向传播、检波器波场的反向传播和应用一个互相关成像条件;
震源波场沿时间正向传播,检波器波场从最大记录时刻开始沿时间反向传播;观测数据是以真实速度模型作为输入,通过波动方程的有限差分正演模拟得到的数据;
步骤2)中,Born正演是通过已经得到的地震成像剖面得到仿真地震数据的过程,具体包括:
2.1)震源波场的正向传播得到背景波场ps(x),其中,背景波场是以偏移速度模型作为输入通过波动方程的有限差分正演模拟得到;
2.3)在检波点位置记录不同时刻的散射波场得到仿真地震数据;
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010352142.XA CN111505718B (zh) | 2020-04-28 | 2020-04-28 | 一种高分辨率地下结构保幅成像方法 |
US17/242,373 US11294086B2 (en) | 2020-04-28 | 2021-04-28 | Method of high-resolution amplitude-preserving seismic imaging for subsurface reflectivity model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010352142.XA CN111505718B (zh) | 2020-04-28 | 2020-04-28 | 一种高分辨率地下结构保幅成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111505718A CN111505718A (zh) | 2020-08-07 |
CN111505718B true CN111505718B (zh) | 2021-11-19 |
Family
ID=71864399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010352142.XA Active CN111505718B (zh) | 2020-04-28 | 2020-04-28 | 一种高分辨率地下结构保幅成像方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11294086B2 (zh) |
CN (1) | CN111505718B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111505718B (zh) * | 2020-04-28 | 2021-11-19 | 西安交通大学 | 一种高分辨率地下结构保幅成像方法 |
CN115166827B (zh) * | 2022-07-15 | 2023-04-28 | 中山大学 | 基于反褶积成像条件的最小二乘偏移成像方法、设备及存储介质 |
CN115951401B (zh) * | 2022-07-19 | 2023-09-15 | 中山大学 | 成像条件驱动的最小二乘逆时偏移成像方法、设备及存储介质 |
CN115993650B (zh) * | 2023-03-22 | 2023-06-06 | 中国石油大学(华东) | 一种基于棱柱波的地震干涉成像方法 |
CN117647838B (zh) * | 2024-01-29 | 2024-04-12 | 东北石油大学三亚海洋油气研究院 | 一种角度约束的逆时偏移成像方法 |
CN118131333B (zh) * | 2024-05-06 | 2024-08-27 | 南方科技大学 | 一种最小二乘逆时偏移成像方法、系统、终端及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016001750A3 (en) * | 2014-06-30 | 2016-03-31 | Cgg Services Sa | Seismic data processing using matching filter based cost function optimization |
CN107102355A (zh) * | 2017-04-27 | 2017-08-29 | 吉林大学 | 低频重构并行Marchenko成像方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3163328B1 (en) * | 2015-11-02 | 2021-09-29 | CGG Services SAS | Seismic data least-square migration method and device |
US10191168B2 (en) * | 2016-01-12 | 2019-01-29 | Cgg Services Sas | AVA compliant pre-stack frequency spectrum enhancement of seismic data |
WO2019182564A1 (en) * | 2018-03-20 | 2019-09-26 | Equinor Us Operations Llc. | Sub-surface structure determination using seismic migration |
US20210215824A1 (en) * | 2020-01-14 | 2021-07-15 | Emerson Paradigm Holding Llc | Seismic migration techniques for improved image accuracy |
CN111505718B (zh) * | 2020-04-28 | 2021-11-19 | 西安交通大学 | 一种高分辨率地下结构保幅成像方法 |
-
2020
- 2020-04-28 CN CN202010352142.XA patent/CN111505718B/zh active Active
-
2021
- 2021-04-28 US US17/242,373 patent/US11294086B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016001750A3 (en) * | 2014-06-30 | 2016-03-31 | Cgg Services Sa | Seismic data processing using matching filter based cost function optimization |
CN107102355A (zh) * | 2017-04-27 | 2017-08-29 | 吉林大学 | 低频重构并行Marchenko成像方法 |
Non-Patent Citations (5)
Title |
---|
A new correlation based least-squares reverse time migration;Feipeng Li et al.;《2018 SEG International Exposition and 88th Annual Meeting》;20181231;第4237页 * |
Improved iterative least-squares migration using curvelet-domain Hessian filters;Ming Wang et al.;《2017 SEG International Exposition and 87th Annual Meeting》;20171231;第4555-4556页 * |
Least-squares RTM: Reality and possibilities for subsalt imaging;Ping Wang et al.;《2016 SEG International Exposition and 86th Annual Meeting》;20161231;第4204-4206页 * |
Ming Wang et al..Improved iterative least-squares migration using curvelet-domain Hessian filters.《2017 SEG International Exposition and 87th Annual Meeting》.2017, * |
Ping Wang et al..Least-squares RTM: Reality and possibilities for subsalt imaging.《2016 SEG International Exposition and 86th Annual Meeting》.2016, * |
Also Published As
Publication number | Publication date |
---|---|
CN111505718A (zh) | 2020-08-07 |
US20210341635A1 (en) | 2021-11-04 |
US11294086B2 (en) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111505718B (zh) | 一种高分辨率地下结构保幅成像方法 | |
US10605941B2 (en) | Methods for simultaneous source separation | |
CN110687600B (zh) | 一种基于声-弹耦合方程的弹性波最小二乘逆时偏移方法 | |
CN108983287B (zh) | 一种基于凸集投影算法的曲波变换反假频地震数据重建方法 | |
CN112882115B (zh) | 基于gwo优化小波阈值的大地电磁信号去噪方法及系统 | |
CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
CN111077567B (zh) | 一种基于矩阵乘法的双程波叠前深度偏移的方法 | |
CN112147236A (zh) | 一种基于稀疏盲解卷积的超声信号分辨率提升方法 | |
Métivier et al. | A review of the use of optimal transport distances for high resolution seismic imaging based on the full waveform | |
CN110749925A (zh) | 一种宽频逆时偏移成像处理方法 | |
CN114325821A (zh) | 基于3d-snacnn网络的叠前地震资料中强散射噪声压制方法及系统 | |
CN117724154A (zh) | 一种考虑强各向异性特征的vti介质参数预测方法 | |
CN113866826A (zh) | 一种混合域地震偏移海森矩阵估计方法 | |
CN110967734A (zh) | 基于快速傅立叶变换的虚源重构方法及系统 | |
CN115598704A (zh) | 一种基于最小二乘逆时偏移生成保幅角道集的方法、设备及可读存储介质 | |
CN110749923A (zh) | 一种基于范数方程提高分辨率的反褶积方法 | |
CN114721044A (zh) | 一种多频率接收函数和振幅比联合反演地壳结构的方法及系统 | |
CN111077566B (zh) | 一种基于矩阵分解的双程波叠前深度偏移的方法 | |
CN110596756B (zh) | 基于自适应混合复扩散模型的沙漠地震勘探噪声压制方法 | |
CN113866821A (zh) | 一种基于照明方向约束的被动源干涉偏移成像方法和系统 | |
CN107589456B (zh) | 一种获取地震数据的方法和装置及计算机可读存储介质 | |
CN114779335B (zh) | 基于各向异性全变分约束的弹性波直接包络反演方法 | |
CN113064205B (zh) | 菲涅尔带约束的浅水多次波衰减方法 | |
Henley | More coherent noise attenuation in the radial trace domain | |
Wu et al. | An inverse Q filtering method with efficient Q estimation using a shallow neural network |
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 |