CN106526674A - 一种三维全波形反演能量加权梯度预处理方法 - Google Patents
一种三维全波形反演能量加权梯度预处理方法 Download PDFInfo
- Publication number
- CN106526674A CN106526674A CN201610999220.9A CN201610999220A CN106526674A CN 106526674 A CN106526674 A CN 106526674A CN 201610999220 A CN201610999220 A CN 201610999220A CN 106526674 A CN106526674 A CN 106526674A
- Authority
- CN
- China
- Prior art keywords
- wave
- wave field
- gradient
- energy
- field
- 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
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000007781 pre-processing Methods 0.000 title abstract 4
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- FNMKZDDKPDBYJM-UHFFFAOYSA-N 3-(1,3-benzodioxol-5-yl)-7-(3-methylbut-2-enoxy)chromen-4-one Chemical compound C1=C2OCOC2=CC(C2=COC=3C(C2=O)=CC=C(C=3)OCC=C(C)C)=C1 FNMKZDDKPDBYJM-UHFFFAOYSA-N 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims 1
- 238000003384 imaging method Methods 0.000 abstract description 5
- 238000009792 diffusion process Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 7
- 230000007423 decrease Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 3
- 230000006854 communication Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011084 recovery Methods 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/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data 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/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
本方法涉及一种三维全波形反演能量加权梯度预处理方法,其主要步骤为:(1)时间二阶空间有限差分方法三维正演模拟,存储边界波场同时计算并存储到达每个速度网格点的振幅最大值;(2)同时进行源波场重建和残差波场逆时传播,对应时刻互相关求取梯度;(3)读取步骤(1)中存储的信息求取能量加权算子对梯度进行预处理;(4)求取合适步长作用在预处理后的梯度算子求取速度更新量迭代更新速度。本方法在不增加计算量的前提下解决了由于球面波传播几何扩散等原因造成的梯度算子深浅层能量不均衡的问题,有效的提高了全波形反演对深层速度成像的精度。
Description
技术领域
本发明涉及地球物理勘查技术领域,主要用于提高三维全波形反演纵波速度场的精度。
背景技术
速度是描述地下介质情况的重要参数,地球物理勘探的重点在于如何恢复地下介质所有尺度的速度信息。速度分析的方法有很多,但是始终达不到理想的效果。射线层析、波动方程层析等层析成像的方法可恢复低波数的速度信息,偏移的方法可提供高波数的反射率信息,然而这些方法都不能同时恢复所有波数的速度参数。全波形反演利用叠前地震数据的全部信息,可提供地下介质的高分辨率成像。常规全波形反演是非线性梯度类最优化方法,以观测数据与模拟数据的最小二乘目标函数值最小为标准更新地下速度模型。三维全波形反演主要存在以下难点:(1)巨大的计算量问题:三维全波形反演基于三维波动方程正演模拟的迭代反演方法,在一次迭代过程中,需要计算多炮的波动方程正演模拟,误差波场反向传播,这两部分每个时间步都需要进行差分运算,伴随着巨大的计算量。(2)海量存储的问题:基于伴随状态梯度算子的全波形反演方法要求存储整个正演波场,三维情况下的存储量是不可容忍的,可以说大部分的时间消耗在数据的I/O读写上。计算量大可以通过多种并行算法进行加速,而I/O读写速度只能靠提高硬件性能来改善。(3)梯度算子的求取效率与精度:常规全波形反演梯度算子是由震源正传波场的时间二阶偏导数与残差波场反传波场互相关所得,震源正传波场包含着球面波传播过程中几何扩散的能量损失,震源波场与残差反传波场互相关就导致了梯度深、浅层能量更加不均衡,会引起反演速度深层精度不足。(4)海森矩阵的求取:对梯度预处理同样是全波形反演中的一个重要环节,可有效的提高反演的精度。然而海森矩阵的求取伴随着巨大的存储量和计算量,如何避开海森矩阵的求取或是对海森矩阵进一步的近似都是需要考虑的问题。如何在提高效率、降低耗费且保证精度是三维全波形反演的首要问题。因此需要发展一种三维全波形反演能量加权梯度预处理方法。
发明内容
本发明的目的正是为了解决现有技术存在的问题,提供一种高效高精度的三维全波形反演能量加权梯度预处理方法。
本发明的技术方案为:
一种三维全波形反演能量加权梯度预处理方法,具体步骤如下:
(1)三维高精度有限差分正演模拟
第一,根据初始速度的道头信息确定三维正演模拟的观测系统,主要确定炮检波点的位置;第二,根据初始速度的最大最小值求取满足有限差分数值模拟差分稳定性和频散关系的三维正演模拟参数;第三,引入完全匹配层边界条件,确定其需要的参数,用于消除正演模拟的边界反射干扰;第四,用时间二阶空间十阶有限差分方法进行三维正演模拟,存储边界波场用于波场重建,并计算存储到达每个速度网格点振幅最大值。
(2)源波场重建及误差波场逆时传播求取梯度
正演模拟得到的模拟炮集与实际观测炮集对应做差求取残差波场。读取步骤(1)中存储的边界波场作为边界条件用时间二阶空间十阶有限差分方法进行源波场重建。源波场重建的同时进行残差波场逆时传播。根据式(1)的伴随状态法进行梯度计算,将源波场重建的波场与残差波场逆时传播得到的波场对应时刻进行零延迟互相关得到梯度算子。
其中E为目标函数,表示梯度,m为模型参数,v表示各网格点的速度值,xs表示震源点的位置,x表示各网格点的位置,t表示每个时间步,T表示最大时间步,为源波场关于时间的二阶偏导数,q为以残差波场为震源的逆时反传波场。
(3)能量加权梯度预处理
取存储的波场到达每个速度网格点能量的最大值,即进而求取每点能量的最大值为波场到达每个速度网格点能量的最大值即初至波的能量值,表征的是波传播球面波几何扩散的过程,用此对梯度进行预处理,得到能量加权梯度算子。
其中,表示能量加权梯度算子。
(4)求取合适步长迭代更新速度
首先给一个试探步长,再用Armijo条件的一维线搜索方法求取合适的步长作用与能量加权梯度算子对速度进行迭代更新。
本发明的技术效果体现在:
常规全波形反演梯度算子是由震源正传波场的时间二阶偏导数与残差波场反传波场互相关所得,震源正传波场包含着球面波传播过程中几何扩散的能量损失,震源波场与残差反传波场互相关就导致了梯度深、浅层能量更加不均衡,会引起反演速度深层精度不足。利用海森矩阵对梯度进行预处理可有效的均衡梯度的能量,但是三维情况下的海森矩阵构建困难,计算量与存储量巨大,因此我们只提取海森矩阵中球面波传播几何扩散的信息对梯度进行预处理达到均衡梯度能量的效果,提高速度深层的反演精度。
该方法基于GPU集群搭建多级异构并行算法解决三维全波形反演的运算效率问题;利用源波场重建策略解决三维情况下的海量存储问题;在保证运算效率的情况下提高梯度算子构建精度。
附图说明
图1为截取的部分二维Marmousi模型并重采样作为真实速度
图2为对图1的真实速度平滑作为初始速度
图3为不同全波形反演算法迭代60次速度的对比,其中图a为常规算法反演速度;图b为本文算法反演速度
图4为常规全波形算法迭代60次与本文算法迭代60次误差下降曲线的对比
图5为不同全波形反演算法迭代60次单道速度的对比,其中图a为常规算法反演速度;图b为本文算法反演速度
图6为全波形反演能量加权梯度预处理方法流程图
图7为三维SEG/EAGE推覆体速度模型,其中图a为真实速度;图b为初始速度
图8为三维SEG/EAGE推覆体模型不同迭代次数全波形反演速度,其中图a为第1次迭代的三维速度体与三个不同方向的剖面;图b为第10次迭代的三维速度体与三个不同方向的剖面;图c为第40次迭代的三维速度体与三个不同方向的剖面
图9为三维SEG/EAGE推覆体模型不同迭代次数的误差下降曲线
图10为三维SEG/EAGE推覆体模型最终反演速度与真实速度以及初始速度的对比,其中图a为真实速度;图b为初始速度;图c为反演速度
具体实施方式
下面结合附图说明本发明的具体实施方式:
通过模型测试说明具体的技术方案:
第一步:三维高精度有限差分正演模拟
第一,根据初始速度的道头信息确定三维正演模拟的观测系统,主要确定炮检波点的位置;第二,根据初始速度的最大最小值求取满足有限差分数值模拟差分稳定性和频散关系的三维正演模拟参数;第三,引入完全匹配层边界条件,确定边界厚度以及衰减吸收系数,波场能量到达边界处逐渐衰减为零,从而消除三维正演模拟的边界反射干扰;第四,用时间二阶空间十阶有限差分方法对速度模型进行三维正演模拟,存储边界波场用于波场重建,并计算波场到达每个速度网格点振幅最大值。
第二步:源波场重建及误差波场逆时传播求取梯度
第一步中三维正演模拟得到的模拟炮集与实际观测炮集对应做差求取残差波场。读取第一步中存储的边界波场作为边界条件,用时间二阶空间十阶有限差分方法进行源波场重建。源波场重建的同时进行残差波场逆时传播。根据下式(1)的伴随状态法进行梯度计算,将源波场重建的波场与残差波场逆时传播得到的波场对应时刻进行零延迟互相关得到梯度算子。
式中E为目标函数,表示梯度,m为模型参数,v表示各网格点的速度值,xs表示震源点的位置,x表示各网格点的位置,t表示每个时间步,T表示最大时间步,为源波场关于时间的二阶偏导数,q为以残差波场为震源的逆时反传波场。
第三步:能量加权梯度预处理
读取存储的源波场到达每个速度网格点振幅的最大值,即进而求取每点能量的最大值为(前一个公式是最大振幅值,后一个公式是最大能量值)。源波场到达每个速度网格点能量的最大值即初至波的能量值,表征的是波传播球面波几何扩散的过程,用此对梯度进行预处理,得到能量加权梯度算子,
其中,表示能量加权梯度算子。
第四步:求取合适步长迭代更新速度
每次迭代,首先给一个试探步长,再用一维线搜索方法求取新的步长,此步长满足Armijo条件使误差下降,然后将求得的步长作用于能量加权梯度算子对速度进行迭代更新。图3到图5展示的是二维Marmousi模型反演测试的效果图。从图3可以看出图b中的速度成像浅层与深层的能量较图a更加均衡。而从图4误差曲线下降图的也可以看出文中基于能量加权梯度算子的全波形反演算法误差收敛更快,且收敛的误差值更小。从图5为单道速度反演结果的对比可以看出基于能量加权梯度算子的全波形反演算法深部速度更新更加准确。
图8到图10展示的是按照图6流程进行三维SEG/EAGE推覆体模型反演测试的效果图。图8展示的是不同迭代次数三维SEG/EAGE推覆体模型全波形反演的速度结果。其中图a为第1次迭代的三维速度体与三个不同方向的剖面;图b为第10次迭代的三维速度体与三个不同方向的剖面;图c为第40次迭代的三维速度体与三个不同方向的剖面。三个不同方向分别为横纵测线方向8000m处的剖面与深度3000m处的切片。从图7中可以看出随着迭代次数的增加,构造细节成像越来越清晰,特别对于3000m处目的层切片的河道处的细节恢复明显。图9是迭代过程中的误差下降曲线,误差做过归一化。从图中可以看出,17次之前误差下降迅速,17次之后逐渐趋向收敛。图10是三维SEG/EAGE推覆体模型全波形反演迭代第40次的最终反演速度。图a是初始速度,图b是真实速度,图c是全波形反演速度。从图中可以看出速度深浅层成像效果均衡,并可以看出全波形反演在构造细节上优势明显,可以对类似于河道等的构造细节精细刻画。
Claims (3)
1.一种三维全波形反演能量加权梯度预处理方法,其特征在于包括:
(1)三维高精度有限差分正演模拟
第一,根据初始速度的道头信息确定三维正演模拟的观测系统,主要确定炮检波点的位置;第二,根据初始速度的最大最小值求取满足有限差分数值模拟差分稳定性和频散关系的三维正演模拟参数;第三,引入完全匹配层边界条件,确定其需要的参数,用于消除正演模拟的边界反射干扰;第四,用时间二阶空间十阶有限差分方法进行三维正演模拟,存储边界波场用于波场重建,并计算存储到达每个速度网格点振幅最大值;
(2)源波场重建及误差波场逆时传播求取梯度
正演模拟得到的模拟炮集与实际观测炮集对应做差求取残差波场,读取步骤(1)中存储的边界波场作为边界条件用时间二阶空间十阶有限差分方法进行源波场重建,源波场重建的同时进行残差波场逆时传播;根据式(1)的伴随状态法进行梯度计算,将源波场重建的波场与残差波场逆时传播得到的波场对应时刻进行零延迟互相关得到梯度算子;
其中E为目标函数,表示梯度,m为模型参数,v表示各网格点的速度值,xs表示震源点的位置,x表示各网格点的位置,t表示每个时间步,T表示最大时间步,为源波场关于时间的二阶偏导数,q为以残差波场为震源的逆时反传波场;
(3)能量加权梯度预处理
读取存储的波场到达每个速度网格点能量的最大值,进而求取每点能量的最大值,波场到达每个速度网格点能量的最大值即初至波的能量值,表征的是波传播球面波几何扩散的过程,用此对梯度进行预处理,得到能量加权梯度算子;
(4)求取合适步长迭代更新速度
首先给一个试探步长,再用Armijo条件的一维线搜索方法求取合适的步长作用于能量加权梯度算子对速度进行迭代更新。
2.如权利要求1所述的一种三维全波形反演能量加权梯度预处理方法,其特征在于:
步骤(2)中源波场重建的同时进行残差波场逆时传播,根据式(1)的伴随状态法进行梯度计算,将源波场重建的波场与残差波场逆时传播得到的波场对应时刻进行零延迟互相关得到梯度算子;
其中E为目标函数,表示梯度,m为模型参数,v表示各网格点的速度值,xs表示震源点的位置,x表示各网格点的位置,t表示每个时间步,T表示最大时间步,为源波场关于时间的二阶偏导数,q为以残差波场为震源的逆时反传波场。
3.如权利要求2所述的一种三维全波形反演能量加权梯度预处理方法,其特征在于:
步骤(3)能量加权梯度预处理
读取存储的波场到达每个速度网格点能量的最大值,即进而求取每点能量的最大值为得到能量加权梯度算子:
其中,表示能量加权梯度算子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610999220.9A CN106526674B (zh) | 2016-11-14 | 2016-11-14 | 一种三维全波形反演能量加权梯度预处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610999220.9A CN106526674B (zh) | 2016-11-14 | 2016-11-14 | 一种三维全波形反演能量加权梯度预处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106526674A true CN106526674A (zh) | 2017-03-22 |
CN106526674B CN106526674B (zh) | 2020-03-31 |
Family
ID=58351909
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610999220.9A Active CN106526674B (zh) | 2016-11-14 | 2016-11-14 | 一种三维全波形反演能量加权梯度预处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106526674B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950596A (zh) * | 2017-04-11 | 2017-07-14 | 中国石油大学(华东) | 一种基于子波迭代估计的有限差分对比源全波形反演方法 |
CN107505654A (zh) * | 2017-06-23 | 2017-12-22 | 中国海洋大学 | 基于地震记录积分的全波形反演方法 |
CN107894618A (zh) * | 2017-11-10 | 2018-04-10 | 中国海洋大学 | 一种基于模型平滑算法的全波形反演梯度预处理方法 |
CN108181650A (zh) * | 2017-12-19 | 2018-06-19 | 中国石油天然气集团公司 | 一种数据处理方法和装置 |
CN108693560A (zh) * | 2017-04-12 | 2018-10-23 | 中国石油化工股份有限公司 | 一种基于互相关道的散射波成像方法及系统 |
CN108802818A (zh) * | 2018-06-11 | 2018-11-13 | 中国石油大学(北京) | 一种全波形反演中原始梯度的层析分量提取方法 |
CN110058302A (zh) * | 2019-05-05 | 2019-07-26 | 四川省地质工程勘察院 | 一种基于预条件共轭梯度加速算法的全波形反演方法 |
CN110531410A (zh) * | 2019-09-23 | 2019-12-03 | 中国海洋大学 | 一种基于直达波场的最小二乘逆时偏移梯度预条件方法 |
CN111190224A (zh) * | 2020-01-09 | 2020-05-22 | 中国石油大学(华东) | 一种基于三维地震波反向照明的动态采样全波形反演系统及方法 |
CN112578431A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 一种限存状态全波形反演波场最优化存储方法及系统 |
CN113296146A (zh) * | 2021-05-19 | 2021-08-24 | 中国海洋大学 | 一种基于梯度道集相关加权的全波形反演梯度预处理方法 |
CN113341456A (zh) * | 2021-07-08 | 2021-09-03 | 清华大学 | 一种偏移成像方法、装置及电子设备 |
CN113376695A (zh) * | 2021-06-11 | 2021-09-10 | 中国矿业大学 | 一种适用于煤层底板复杂陷落柱的全波形反演方法 |
CN114624766A (zh) * | 2022-05-16 | 2022-06-14 | 中国海洋大学 | 基于行波分离的弹性波最小二乘逆时偏移梯度求取方法 |
CN117849879A (zh) * | 2023-12-20 | 2024-04-09 | 京全品质能源技术(北京)有限公司 | 波动方程反演方法及装置、电子设备以及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100322032A1 (en) * | 2009-06-17 | 2010-12-23 | Shin Chang-Soo | Apparatus and method for imaging subsurface structure of target area by using waveform inversion |
CN103270430A (zh) * | 2011-06-08 | 2013-08-28 | 雪佛龙美国公司 | 用于地震数据反演的系统和方法 |
-
2016
- 2016-11-14 CN CN201610999220.9A patent/CN106526674B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100322032A1 (en) * | 2009-06-17 | 2010-12-23 | Shin Chang-Soo | Apparatus and method for imaging subsurface structure of target area by using waveform inversion |
CN103270430A (zh) * | 2011-06-08 | 2013-08-28 | 雪佛龙美国公司 | 用于地震数据反演的系统和方法 |
Non-Patent Citations (3)
Title |
---|
ZHIGANG ZHANG,等: "DOUBLE-DIFFERENCE ELASTIC-WAVEFORM INVERSION WITH WEIGHTED GRADIENTS FOR MONITORING EGS RESERVOIRS", 《PROCEEDINGS, THIRTY-SEVENTH WORKSHOP ON GEOTHERMAL RESERVOIR ENGINEERING 》 * |
张凯,等: "一种新的用于全波形反演的能量加权梯度方法", 《石油物探》 * |
曲英铭,等: "基于多尺度双变网格的时间域全波形反演", 《石油物探》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106950596B (zh) * | 2017-04-11 | 2019-02-26 | 中国石油大学(华东) | 一种基于子波迭代估计的有限差分对比源全波形反演方法 |
CN106950596A (zh) * | 2017-04-11 | 2017-07-14 | 中国石油大学(华东) | 一种基于子波迭代估计的有限差分对比源全波形反演方法 |
CN108693560B (zh) * | 2017-04-12 | 2020-04-07 | 中国石油化工股份有限公司 | 一种基于互相关道的散射波成像方法及系统 |
CN108693560A (zh) * | 2017-04-12 | 2018-10-23 | 中国石油化工股份有限公司 | 一种基于互相关道的散射波成像方法及系统 |
CN107505654A (zh) * | 2017-06-23 | 2017-12-22 | 中国海洋大学 | 基于地震记录积分的全波形反演方法 |
CN107505654B (zh) * | 2017-06-23 | 2019-01-29 | 中国海洋大学 | 基于地震记录积分的全波形反演方法 |
CN107894618A (zh) * | 2017-11-10 | 2018-04-10 | 中国海洋大学 | 一种基于模型平滑算法的全波形反演梯度预处理方法 |
CN107894618B (zh) * | 2017-11-10 | 2018-08-21 | 中国海洋大学 | 一种基于模型平滑算法的全波形反演梯度预处理方法 |
CN108181650A (zh) * | 2017-12-19 | 2018-06-19 | 中国石油天然气集团公司 | 一种数据处理方法和装置 |
CN108802818A (zh) * | 2018-06-11 | 2018-11-13 | 中国石油大学(北京) | 一种全波形反演中原始梯度的层析分量提取方法 |
CN110058302A (zh) * | 2019-05-05 | 2019-07-26 | 四川省地质工程勘察院 | 一种基于预条件共轭梯度加速算法的全波形反演方法 |
CN110531410A (zh) * | 2019-09-23 | 2019-12-03 | 中国海洋大学 | 一种基于直达波场的最小二乘逆时偏移梯度预条件方法 |
CN112578431A (zh) * | 2019-09-27 | 2021-03-30 | 中国石油化工股份有限公司 | 一种限存状态全波形反演波场最优化存储方法及系统 |
CN112578431B (zh) * | 2019-09-27 | 2024-04-09 | 中国石油化工股份有限公司 | 一种限存状态全波形反演波场最优化存储方法及系统 |
CN111190224A (zh) * | 2020-01-09 | 2020-05-22 | 中国石油大学(华东) | 一种基于三维地震波反向照明的动态采样全波形反演系统及方法 |
CN113296146A (zh) * | 2021-05-19 | 2021-08-24 | 中国海洋大学 | 一种基于梯度道集相关加权的全波形反演梯度预处理方法 |
CN113376695A (zh) * | 2021-06-11 | 2021-09-10 | 中国矿业大学 | 一种适用于煤层底板复杂陷落柱的全波形反演方法 |
CN113376695B (zh) * | 2021-06-11 | 2022-07-05 | 中国矿业大学 | 一种适用于煤层底板复杂陷落柱的全波形反演方法 |
CN113341456A (zh) * | 2021-07-08 | 2021-09-03 | 清华大学 | 一种偏移成像方法、装置及电子设备 |
CN113341456B (zh) * | 2021-07-08 | 2024-02-02 | 清华大学 | 一种基于非规则网格的地震偏移方法、装置及电子设备 |
CN114624766A (zh) * | 2022-05-16 | 2022-06-14 | 中国海洋大学 | 基于行波分离的弹性波最小二乘逆时偏移梯度求取方法 |
CN117849879A (zh) * | 2023-12-20 | 2024-04-09 | 京全品质能源技术(北京)有限公司 | 波动方程反演方法及装置、电子设备以及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106526674B (zh) | 2020-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106526674A (zh) | 一种三维全波形反演能量加权梯度预处理方法 | |
CN111239802B (zh) | 基于地震反射波形和速度谱的深度学习速度建模方法 | |
CN108345031B (zh) | 弹性介质主动源和被动源混采地震数据全波形反演方法 | |
CN103713315B (zh) | 一种地震各向异性参数全波形反演方法及装置 | |
CN110058302A (zh) | 一种基于预条件共轭梯度加速算法的全波形反演方法 | |
CN101334483B (zh) | 一种在地震数据处理中衰减瑞雷波散射噪声的方法 | |
Pan et al. | Love-wave waveform inversion in time domain for shallow shear-wave velocity | |
CN107894618B (zh) | 一种基于模型平滑算法的全波形反演梯度预处理方法 | |
CN107765308B (zh) | 基于褶积思想与精确震源的重构低频数据频域全波形反演方法 | |
CN104570082B (zh) | 一种基于格林函数表征的全波形反演梯度算子的提取方法 | |
CN105093278B (zh) | 基于激发主能量优化算法的全波形反演梯度算子提取方法 | |
CN106842295A (zh) | 测井信息约束的波形反演方法 | |
CN104237937B (zh) | 叠前地震反演方法及其系统 | |
CN103499835A (zh) | 初至波波形反演近地表速度模型方法 | |
CN107505654A (zh) | 基于地震记录积分的全波形反演方法 | |
CN113740901B (zh) | 基于复杂起伏地表的陆上地震数据全波形反演方法及装置 | |
CN113552625B (zh) | 一种用于常规陆域地震数据的多尺度全波形反演方法 | |
CN105319581A (zh) | 一种高效的时间域全波形反演方法 | |
CN104570101A (zh) | 一种基于粒子群算法的avo三参数反演方法 | |
CN106569262B (zh) | 低频地震数据缺失下的背景速度模型重构方法 | |
CN109459787B (zh) | 基于地震槽波全波形反演的煤矿井下构造成像方法及系统 | |
CN104965222B (zh) | 三维纵波阻抗全波形反演方法及装置 | |
BR112012022909B1 (pt) | Métodos de exploração de hidrocarbonetos dentro de volume predeterminado de terra | |
CN107462924A (zh) | 一种不依赖于测井资料的绝对波阻抗反演方法 | |
CN103777236A (zh) | 一种可控震源同步扫描波场正演数值模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |