CN113866826A - 一种混合域地震偏移海森矩阵估计方法 - Google Patents
一种混合域地震偏移海森矩阵估计方法 Download PDFInfo
- Publication number
- CN113866826A CN113866826A CN202110970271.XA CN202110970271A CN113866826A CN 113866826 A CN113866826 A CN 113866826A CN 202110970271 A CN202110970271 A CN 202110970271A CN 113866826 A CN113866826 A CN 113866826A
- Authority
- CN
- China
- Prior art keywords
- migration
- domain
- mixed
- seismic
- imaging result
- 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
- 230000005012 migration Effects 0.000 title claims abstract description 54
- 238000013508 migration Methods 0.000 title claims abstract description 54
- 239000011159 matrix material Substances 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000003384 imaging method Methods 0.000 claims abstract description 69
- 238000002310 reflectometry Methods 0.000 claims abstract description 11
- 230000005526 G1 to G0 transition Effects 0.000 claims abstract description 5
- 238000001514 detection method Methods 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 3
- 230000003094 perturbing effect Effects 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 10
- 150000003839 salts Chemical class 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005286 illumination Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000059 patterning Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction 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. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- 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
本发明公开了一种混合域地震偏移海森矩阵估计方法,涉及勘探地球物理技术领域,该方法首先利用常规偏移方法获得初始成像结果,然后将该偏移结果作为反射率模型,进行线性正演和偏移获得二次成像结果,使用这两个偏移结果在空间‑波数域估计非稳相滤波器,用于近似汉森矩阵,并利用混合域反褶积计算汉森矩阵的逆,再使用估计的汉森矩阵的逆在最小二乘偏移中作为泛函的梯度,以提高反演的收敛速率,使得最小二乘反演在4‑5次迭代内收敛至准确解。本发明的有益效果是,该混合域地震偏移汉森矩阵估计方法既具有较高的计算效率,又可加快最小二乘偏移收敛速度,降低最小二乘反射率反演的计算成本,具有较好的应用前景。
Description
技术领域
本发明涉及勘探地球物理技术领域,尤其涉及一种混合域地震偏移海森矩阵估计方法。
背景技术
目前,高精度地震偏移成像是油气勘探中重要的关键技术之一,常规地震偏移方法,包括:Kirchhoff偏移、高斯束偏移、单程波偏移以及逆时偏移等,在数学意义上均可看成线性正演模拟的伴随算子,而非真正意义的逆算子。因此,当地震采集不规则、震源及检波器频带较窄以及地下照明不均匀时,基于伴随算子的常规偏移方法很难获得高质量成像结果。基于反演框架的最小二乘偏移通过数据拟合可以逐渐逼近正演算子的逆,进而大幅提高偏移结果的振幅均衡性和空间分辨率,是目前地震成像领域研究的热点。
最小二乘偏移成像的核心在于如何高效计算汉森矩阵并求其逆算子。常规数据域最小二乘偏移首先构建观测数据和预测数据的误差泛函,然后使用线性正演计算预测数据,使用伴随偏移计算泛函梯度,并通过迭代的方式逐步逼近汉森矩阵的逆。该方法通常需要数十次至上百次正演模拟和伴随偏移,计算效率较低,往往很难应用于大规模实际生产处理中。通过设计合理的预条件算子,可以加快最小二乘偏移的收敛速率,但目前常用的震源照明补偿只能改善振幅的相对强弱关系,不能够提高空间分辨率。
为使最小二乘偏移走向生产化应用,急需研发精确高效的预条件算子,使其使用较少次数的迭代就可以收敛至真实解附近。
发明内容
为解决上述最小二乘偏移成像中的收敛速度较慢、计算效率较低的技术问题,本发明公开了一种混合域地震偏移海森矩阵估计方法。
为实现上述目的,本发明采用下述技术方案:
一种混合域地震偏移海森矩阵估计方法,具体包括以下步骤:
步骤一,获取输入数据;
步骤二,通过求解波动方程计算震源波场,计算检波点波场,再采用零延迟互相关成像条件计算初步成像结果;
步骤三,将初步成像结果作为反射率模型,使用线性正演模拟计算合成地震数据;
步骤四,使用步骤三合成的地震数据和偏移速度模型,采用步骤二的成像处理过程进行再次偏移,获得二次成像结果;
步骤五,将初步成像结果和二次成像结果应用空间-波数域变换,获得混合域成像结果;
步骤六,计算近似汉森矩阵的非稳相混合域滤波器;
步骤七,使用混合域反褶积计算汉森矩阵的逆;
步骤八,将估计的汉森矩阵的逆作为预条件算子应用于最小二乘偏移的泛函梯度中,以提高反演的收敛速率;
步骤九,在每次最小二乘迭代中都应用步骤八中的预条件算子,使得最小二乘偏移快速收敛于准确解。
进一步地,步骤一中,输入数据包括:偏移速度模型v(x)、震源函数子波f(t)、野外观测数据dobs(xS,xR,t),混合域滤波器窗口长度l,xR为检波点位置。
进一步地,步骤二包括:
(2.1)使用震源函数子波f(t)和偏移速度模型v(x),通过求解波动方程计算震源波场US(x,t),波动方程如下:
其中,x为地下成像空间位置,xS代表了震源位置,t代表时间变量,δ(x-xs)为狄拉克delta函数。
(2.2)使用野外观测数据dobs(xS,xr,t)和偏移速度模型v(x),计算检波点波场UR(x,t),具体表达式如下:
(2.3)采用零延迟互相关成像条件计算初步成像结果m0(x),且互相关成像条件为:
其中,T为地震记录的长度,US(x,t)为震源波场,为UR(x,t)检波点波场。
进一步地,步骤三中,将初步成像结果m0(x)作为反射率模型,计算合成地震数据dsyn(xS,xr,t)所采用的线性正演模拟表达式如下:
其中,U0(x,t)为背景波场,U1(x,t)为扰动波场。
进一步地,步骤五中,对初次成像结果m0(x)和二次成像结果m1(x)应用空间-波数域变换,该空间-波数域变换表达式为:
其中,x=[x,y,z]T为空间矢量,k=[kx,ky,kz]T为波数矢量,x'=[x',y',z']T为空间-波数域变换后的空间矢量,M0(x',k)、M1(x',k)分别为m0(x)、m1(x)的空间-波数域变换结果,σx、σy和σz分别为在x、y和z方向的窗口长度,设置σx、σy和σz的值为输入混合域滤波器窗口长度l。
进一步地,步骤六中,使用M0(x',k)和M1(x',k),计算近似汉森矩阵的非稳相混合域滤波器F(x',k),其表达式如下:
其中,ε为避免除以零的小值,取ε=0.001*max|M1(x',k)|。
进一步地,步骤七中,汉森矩阵的逆Hinv(x',k)为:
进一步地,步骤八中,预条件算子如下:
其中,g(x',k)为泛函梯度的空间-波数域谱,gpre(x)为应用预条件后的梯度。
本发明的有益效果是,首先计算了均匀介质中解析的汉森矩阵,并通过分析发现当把该矩阵的列向量投影至成像空间时所形成的点拓展函数具有局部化特征,基于该分析,本发明专利提出了一种混合域地震偏移汉森矩阵估计方法,即根据地震偏移的结果构建一系列空间-波数域的滤波器,用于近似局部化的点拓展函数,进而可以通过混合域反褶积计算汉森矩阵的逆。再将该汉森矩阵的逆作为最小二乘偏移的预条件算子,可大幅提高反演的收敛速率,使得最小二乘反演在4-5次迭代内收敛至准确解。该混合域地震偏移汉森矩阵估计方法既具有较高的计算效率,又可加快最小二乘偏移收敛速度,降低最小二乘反射率反演的计算成本,具有较好的应用前景。
附图说明
图1为部分Sigsbee速度模型图;
图2为由传统偏移方法获得的成像结果m0(x);
图3为由m0(x)作为反射率模型,进行线性正演和再次偏移所获得的成像结果m1(x);
图4为对m0(x)应用汉森矩阵的逆获得的校正后的偏移结果;
图5为常规最小二乘偏移迭代更新5次后,获得的反演成像结果
图6为使用本发明汉森矩阵作为预条件算子,优化的最小二乘偏移迭代更新5次后,获得的反演成像结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
基于线性反演框架,最小二乘偏移可以消除由于采集不规则造成的偏移假象,并可以改善深部振幅均衡性和整体的空间分辨率。,但常规数据域最小二乘偏移通过迭代数据拟合逐步逼近汉森矩阵的逆,在实际应用中,往往需要数十次正演模拟和伴随偏移,计算量较大,在当前计算条件下仍较难走向实用化生产应用。
针对最小二乘偏移成像中的收敛速度较慢、计算效率较低的问题,本发明公开了一种混合域地震偏移海森矩阵估计方法,具体包括以下步骤:
(1)获取输入数据。
输入数据包括:偏移速度模型v(x)、震源函数子波f(t)、野外观测数据dobs(xS,xR,t),混合域滤波器窗口长度l,xR为检波点位置。
(2)通过求解波动方程计算震源波场,计算检波点波场,再采用零延迟互相关成像条件计算初步成像结果。
具体地:
(2.1)使用震源函数子波f(t)和偏移速度模型v(x),通过求解波动方程计算震源波场US(x,t),波动方程如下:
其中,x为地下成像空间位置,xS代表了震源位置,t代表时间变量,δ(x-xs)为狄拉克delta函数。
(2.2)使用野外观测数据dobs(xS,xr,t)和偏移速度模型v(x),计算检波点波场UR(x,t),具体表达式如下:
(2.3)采用零延迟互相关成像条件计算初步成像结果m0(x),且互相关成像条件为:
其中,T为地震记录的长度,US(x,t)为震源波场,为UR(x,t)检波点波场。
(3)将初步成像结果作为反射率模型,使用线性正演模拟计算合成地震数据。
将初步成像结果m0(x)作为反射率模型,计算合成地震数据dsyn(xS,xr,t)所采用的线性正演模拟表达式如下:
其中,U0(x,t)为背景波场,U1(x,t)为扰动波场。
(4)使用步骤(3)合成的地震数据dsyn(xS,xr,t)和偏移速度模型,采用步骤(2)的成像处理过程进行再次偏移,获得二次成像结果m1(x);
(5)将初步成像结果和二次成像结果应用空间-波数域变换,获得混合域成像结果M0(x',k)和M1(x',k)。
对初次成像结果m0(x)和二次成像结果m1(x)应用空间-波数域变换,该空间-波数域变换表达式为:
其中,x=[x,y,z]T为空间矢量,k=[kx,ky,kz]T为波数矢量,x'=[x',y',z']T为空间-波数域变换后的空间矢量,M0(x',k)、M1(x',k)分别为m0(x)、m1(x)的空间-波数域变换结果,σx、σy和σz分别为在x、y和z方向的窗口长度,本方法中设置σx、σy和σz的值为输入混合域滤波器窗口长度l。
(6)计算近似汉森矩阵的非稳相混合域滤波器。
使用M0(x',k)和M1(x',k),计算近似汉森矩阵的非稳相混合域滤波器F(x',k),其表达式如下:
其中,ε为避免除以零的小值,取ε=0.001*max|M1(x',k)|。
(7)使用混合域反褶积计算汉森矩阵的逆。
汉森矩阵的逆Hinv(x',k)为:
(8)将估计的汉森矩阵的逆作为预条件算子应用于最小二乘偏移的泛函梯度中,以提高反演的收敛速率。
预条件算子如下:
其中,g(x',k)为泛函梯度的空间-波数域谱,gpre(x)为应用预条件后的梯度。
(9)在每次最小二乘迭代中都应用步骤八中的预条件算子,使得最小二乘偏移快速收敛于准确解。
将本方法应用于国际标准Sigsbee模型地震偏移成像中,取得了较好的成像结果。
图1为部分Sigsbee速度模型图,观测数据的震源函数为15Hz的雷克子波,观测数据集为50个炮集,均匀分布在地表,每个炮集包含500个检波点,间隔为10m。图2为使用观测数据获得的偏移成像结果m0(x),可以看出,常规偏移对于浅部反射层刻画较好,但对于盐下深部沉积层,受盐丘能量屏蔽作用影响,很难获得高质量成像结果。图3为使用图2所示的成像结果m0(x)作为反射率模型,进行线性正演和偏移获得的二次偏移结果m1(x),受汉森矩阵影响,二次偏移结果相对初始偏移结果深部振幅更弱、分辨率更差。利用图2和图3所示的偏移结果在空间-波数域估计汉森矩阵,并应用于初始结果上,获得了图4的成像结果,可以看出深部能量得到了显著补偿;此外,整体的空间分辨率也得到了很大提高。图5和图6对比了常规最小二乘和使用本发明专利提出的汉森矩阵作为预条件的最小二乘偏移在5次迭代更新后获得的成像结果。可以看出,常规最小二乘由于收敛较慢,在经过5次迭代后成像结果相对图2中的成像结果改善不大。但是,在使用本发明所述的汉森矩阵作为预条件算子后,经过5次迭代,盐丘和盐下沉积层的成像质量均得到了大幅度改善。
本发明专利提出了一种混合域地震偏移汉森矩阵估计方法,主要解决常规最小二乘偏移中收敛速率较慢、计算效率较低问题,通过将估计的汉森矩阵作为最小二乘偏移中泛函梯度的预条件算子,可使得最小二乘偏移在4-5次迭代后收敛至准确解。
本发明首先利用常规偏移方法获得初始成像结果,然后将该偏移结果作为反射率模型,进行线性正演和偏移获得二次成像结果,使用这两个偏移结果在空间-波数域估计非稳相滤波器,用于近似汉森矩阵,并利用混合域反褶积计算汉森矩阵的逆,再使用估计的汉森矩阵的逆在最小二乘偏移中作为泛函的梯度,可大题提高收敛速率,使其在4-5次迭代内收敛至准确解,最终推进最小二乘偏移的实用化进程。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。
Claims (8)
1.一种混合域地震偏移海森矩阵估计方法,其特征在于,具体包括以下步骤:
步骤一,获取输入数据;
步骤二,通过求解波动方程计算震源波场,计算检波点波场,再采用零延迟互相关成像条件计算初步成像结果;
步骤三,将初步成像结果作为反射率模型,使用线性正演模拟计算合成地震数据;
步骤四,使用步骤三合成的地震数据和偏移速度模型,采用步骤二的成像处理过程进行再次偏移,获得二次成像结果;
步骤五,将初步成像结果和二次成像结果应用空间-波数域变换,获得混合域成像结果;
步骤六,计算近似汉森矩阵的非稳相混合域滤波器;
步骤七,使用混合域反褶积计算汉森矩阵的逆;
步骤八,将估计的汉森矩阵的逆作为预条件算子,应用于最小二乘偏移的泛函梯度中,以提高反演的收敛速率;
步骤九,在每次最小二乘迭代中都应用步骤八中的预条件算子,使得最小二乘偏移快速收敛于准确解。
2.如权利要求1所述的一种混合域地震偏移海森矩阵估计方法,其特征在于,步骤一中,输入数据包括:偏移速度模型v(x)、震源函数子波f(t)、野外观测数据dobs(xS,xR,t),混合域滤波器窗口长度l,xR为检波点位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110970271.XA CN113866826B (zh) | 2021-08-23 | 2021-08-23 | 一种混合域地震偏移海森矩阵估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110970271.XA CN113866826B (zh) | 2021-08-23 | 2021-08-23 | 一种混合域地震偏移海森矩阵估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113866826A true CN113866826A (zh) | 2021-12-31 |
CN113866826B CN113866826B (zh) | 2024-03-15 |
Family
ID=78987982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110970271.XA Active CN113866826B (zh) | 2021-08-23 | 2021-08-23 | 一种混合域地震偏移海森矩阵估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113866826B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166827A (zh) * | 2022-07-15 | 2022-10-11 | 中山大学 | 基于反褶积成像条件的最小二乘偏移成像方法、设备及存储介质 |
CN115951401A (zh) * | 2022-07-19 | 2023-04-11 | 中山大学 | 成像条件驱动的最小二乘逆时偏移成像方法、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120075954A1 (en) * | 2010-09-24 | 2012-03-29 | CGGVeritas Services (U.S.) Inc. | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration |
CN108369289A (zh) * | 2015-12-18 | 2018-08-03 | 埃克森美孚上游研究公司 | 使用全波场反演点扩展函数分析设计地球物理勘测的方法 |
-
2021
- 2021-08-23 CN CN202110970271.XA patent/CN113866826B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120075954A1 (en) * | 2010-09-24 | 2012-03-29 | CGGVeritas Services (U.S.) Inc. | Device and Method for Calculating 3D Angle Gathers from Reverse Time Migration |
CN108369289A (zh) * | 2015-12-18 | 2018-08-03 | 埃克森美孚上游研究公司 | 使用全波场反演点扩展函数分析设计地球物理勘测的方法 |
Non-Patent Citations (2)
Title |
---|
JIDONG YANG 等: "Approximating the Gauss–Newton Hessian using a space-wavenumber filter and its applications in least-squares seismic imaging", TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, vol. 60, pages 1 - 13, XP011892348, DOI: 10.1109/TGRS.2021.3064852 * |
李闯 等: "预条件最小二乘逆时偏移方法", 石油地球物理勘探, vol. 51, no. 03, pages 513 - 520 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166827A (zh) * | 2022-07-15 | 2022-10-11 | 中山大学 | 基于反褶积成像条件的最小二乘偏移成像方法、设备及存储介质 |
CN115166827B (zh) * | 2022-07-15 | 2023-04-28 | 中山大学 | 基于反褶积成像条件的最小二乘偏移成像方法、设备及存储介质 |
CN115951401A (zh) * | 2022-07-19 | 2023-04-11 | 中山大学 | 成像条件驱动的最小二乘逆时偏移成像方法、设备及存储介质 |
CN115951401B (zh) * | 2022-07-19 | 2023-09-15 | 中山大学 | 成像条件驱动的最小二乘逆时偏移成像方法、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113866826B (zh) | 2024-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Aghamiry et al. | Compound regularization of full-waveform inversion for imaging piecewise media | |
CN105467444B (zh) | 一种弹性波全波形反演方法及装置 | |
CN113866826A (zh) | 一种混合域地震偏移海森矩阵估计方法 | |
CN107894612B (zh) | 一种q吸收衰减补偿的声波阻抗反演方法及系统 | |
CN106483559B (zh) | 一种地下速度模型的构建方法 | |
Maharramov* et al. | Joint full-waveform inversion of time-lapse seismic data sets | |
CN109765616B (zh) | 一种保幅波场延拓校正方法及系统 | |
CN105353405B (zh) | 一种全波形反演方法和系统 | |
CN110488354B (zh) | 一种q补偿的起伏地表棱柱波与一次波联合最小二乘逆时偏移成像方法 | |
CN110187382B (zh) | 一种回折波和反射波波动方程旅行时反演方法 | |
CN111505718A (zh) | 一种高分辨率地下结构保幅成像方法 | |
CN110857998B (zh) | 一种基于lowrank有限差分的弹性逆时偏移方法及系统 | |
CN108983300B (zh) | 一种隧道掘进机施工条件下的瞬变电磁隧道超前预报方法 | |
Zand et al. | Consensus optimization of total variation–based reverse time migration | |
Maharramov et al. | Robust joint full-waveform inversion of time-lapse seismic data sets with total-variation regularization | |
Fang et al. | Source-independent elastic least-squares reverse time migration | |
CN104459770A (zh) | 一种高维地震数据规则化方法 | |
Wason et al. | Source separation via SVD-free rank minimization in the hierarchical semi-separable representation | |
CN112965102A (zh) | 一种基于阻抗敏感核函数的最小二乘逆时偏移成像方法 | |
CN116088048A (zh) | 包含不确定性分析的各向异性介质多参数全波形反演方法 | |
CN116626765A (zh) | 一种基于二维k-svd和卷积稀疏编码的多道地震反褶积方法 | |
Zhang et al. | Multitrace impedance inversion based on structure-oriented regularization | |
Wang et al. | Robust nonstationary local slope estimation | |
Guan et al. | Love wave full-waveform inversion for archaeogeophysics: From synthesis tests to a field case | |
Tu et al. | Fast least-squares migration with multiples and source estimation |
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 |