CN113050163B - 一种振幅、相位信息可调节的全波形反演方法及装置 - Google Patents
一种振幅、相位信息可调节的全波形反演方法及装置 Download PDFInfo
- Publication number
- CN113050163B CN113050163B CN202110489838.1A CN202110489838A CN113050163B CN 113050163 B CN113050163 B CN 113050163B CN 202110489838 A CN202110489838 A CN 202110489838A CN 113050163 B CN113050163 B CN 113050163B
- Authority
- CN
- China
- Prior art keywords
- inversion
- seismic source
- wave field
- gradient
- representing
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 208000037516 chromosome inversion disease Diseases 0.000 claims abstract description 114
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 238000004088 simulation Methods 0.000 claims abstract description 31
- 230000014509 gene expression Effects 0.000 claims abstract description 28
- 238000011896 sensitive detection Methods 0.000 claims abstract description 8
- 230000005540 biological transmission Effects 0.000 claims abstract description 6
- 238000001514 detection method Methods 0.000 claims description 10
- 238000002939 conjugate gradient method Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims description 3
- 230000010354 integration Effects 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 7
- 238000007796 conventional method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000011541 reaction mixture Substances 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
- 238000012795 verification 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/282—Application of seismic models, synthetic seismograms
-
- 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
- G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
- G01V2200/10—Miscellaneous details
- G01V2200/14—Quality control
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
本发明公开了一种振幅、相位信息可调节的全波形反演方法及装置,通过多个单频编码信号同时激发作为震源;采用一阶速度—应力声波动方程,在一次正演模拟中完成多个震源子波的波场模拟,通过相敏检测方法分别对正传和反传波场实施波场解耦计算,求取无串扰的梯度结果,在目标函数关于速度场模型参数的梯度计算表达式计算中,采用在不同的反演阶段使用不同的计算表达式,即在反演初期主要以相位信息反演为主,在反演的中后阶段,采用全波形反演,可以减弱反演对初始模型的依赖性;通过使用相位反演和全波形反演相结合的策略,可以降低目标函数的复杂度,从而减少反演中的多解性,最终实现高精度的反演。
Description
技术领域
本发明涉及一种振幅、相位信息可调节的全波形反演方法及装置,属于地震勘探速度建模技术领域。
背景技术
随着我国石油勘探向深层及超深层转移,深部目标层高陡复杂构造的高精度成像问题已成为当前制约我国石油勘探和开发迫在眉睫的瓶颈。速度建模是实现复杂构造成像的灵魂,而诸如速度层析反演等传统速度建模技术虽在过去的数十年中取得了巨大的成功,但已无法满足高分辨率地震成像的需求。因此,亟需发展新一代速度建模技术,以满足复杂构造成像的需求。
目前,全波形反演(Full Waveform Inversion,FWI)是公认的最具潜力的新一代速度建模技术,在准确的初始模型下能够恢复出高分辨率的速度模型,可用于提高深层目标层偏移成像的质量,具有广阔的应用前景和潜在的经济效益。全波形反演以层析反演或者速度分析结果为初始模型,可实现高分辨率的速度建模,结合保幅的偏移成像方法,获得高质量的成像剖面,为地质解释提供支撑和依据,因此全波形反演成为近些年来的研究热点之一。然而,计算效率低下是制约全波形反演进行大规模计算以及实际资料处理的主要因素,现有的一种解决方法是使用震源编码策略,其思想是通过某种震源组合方式将很多单炮形成一个超级震源,从而减少波场模拟的次数。其中被广泛应用的一种为基于频率选择的同时源全波形反演方法,它可以获得无串扰的反演结果,但是这种方法如果在初始模型不是很理想时,由于波场振幅信息对初始模型比较敏感,使用全波场信息匹配很容易使得全波形反演陷入局部极小值,所以开展数据子集或者波场属性反演就显得十分需要。由于振幅信息对初始模型的敏感度比相位信息高,所以当初始模型不准确时,直接使用全波形反演,可能导致全波形反演低波数更新不准确,从而导致高波数更新的界面位置不对,降低反演的精度。
因此,如何提供一种新的同时源反演方法,能够克服现有技术存在的反演精度低的问题是本领域技术人员亟需解决的技术难题。
发明内容
针对上述现有技术存在的问题,本发明提供一种振幅、相位信息可调节的全波形反演方法及装置,通过一系列单频子波同时激发形成混叠波场;并采用相敏检测法从混叠波场中解耦出每一个编码信号对应的波场变量,实现无串扰的梯度计算,从而有效提高反演精度。
为了实现上述目的,本发明采用的技术方案是:一种振幅、相位信息可调节的全波形反演方法,具体步骤为:
A、通过多个位于不同位置的单频子波同时激发作为震源,采用多个检波点组成观测系统,通过观测系统获得观测数据;建立速度场模型(该速度场模型从现有模型中选择一个),并根据震源情况计算得出模拟数据;
B、基于观测数据和模拟数据,构建表征所述观测数据与所述模拟数据之间数据残差的目标函数;
C、根据震源与检波点之间的位置,得出虚拟震源的信息,然后根据观测数据获取震源的正传波场;并根据目标函数获得虚拟震源的反传波场;
D、确定目标函数关于速度场模型参数的梯度计算表达式,在不同的反演阶段具有两种计算表达式,其中在反演初期即低频反演阶段,采用具有相位信息反演的梯度表达式,在反演的中后阶段,采用兼顾波场相位和振幅信息的全波形反演的梯度表达式;先对震源的正传波场和虚拟震源的反传波场分别进行无串扰解耦,获得解耦的波场变量;然后根据得到的波场变量和梯度计算表达式,计算得到速度场模型参数的梯度;反演初期和反演中后期具体划分为:多尺度全波形反演一般反演3至5个频带,通常前一到两个频带被认为是低波数恢复,也就是低频反演阶段,中间频带的反演恢复中波数信息,后续频带为高频反演。
E、采用共轭梯度方法求取更新方向,根据得到的梯度对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。
进一步,所述步骤B中的目标函数采用最小化位于检波点Xr的观测数据dobs(Xr,Ω,t)和模拟数据u(Xr,Ω,t;m)的残差,具体为:
其中,x表示计算区域,m(x)表示计算区域x处的速度场模型参数;Xr为检波点的位置信息;Ω表示角频率,t表示时间;m表示速度场模型参数,u(Xr,Ω,t;m)表示在当前速度场模型参数m下获得的模拟数据。
进一步,所述步骤A中的震源为同时源模拟波场:
给定一系列位于Xs的单频随机信号s(Xs,Ω),则震源能表达为:
进一步,所述步骤C中震源的正传波场由求解模拟方程得到,具体模拟方程为:
震源对应的虚拟震源为:
根据确定的虚拟震源,获取伴随方程,所述伴随方程为:
求解伴随方程,获得虚拟震源的反传波场。
进一步,所述步骤D中无串扰解耦具体为:
采用相敏检测法先获取第一参考信号和第二参考信号,所述第一参考信号的振幅和频率均与所述第二参考信号的振幅和频率相同,所述第一参考信号的相位与所述第二参考信号的相位相差90°;然后根据第一参考信号和第二参考信号,采用三角函数正交性对震源的正传波场和虚拟震源的反传波场分别进行无串扰解耦。
一种振幅、相位信息可调节的全波形反演装置,包括:
目标函数构建模块,用于基于观测数据和模拟数据,构建表征所述观测数据与所述模拟数据之间数据残差的目标函数;
正传波场获取模块,用于获取震源的正传波场;
反传波场获取模块,用于获取虚拟震源的反传波场;
梯度计算模块,用于根据正传波场、反传波场和目标函数关于速度场模型参数的相位反演、全波场反演梯度计算表达式,计算速度场模型参数的梯度;
迭代更新模块,用于根据计算获得的梯度对当前速度场模型参数进行迭代更新,直至所述数据残差达到设定阈值或者迭代更新次数达到预设值,确定速度场模型参数。
与现有技术相比,本发明采用一系列单频子波同时激发形成混叠波场;并采用了相敏检测法从混叠波场中解耦出每一个编码信号对应的波场变量,实现了无串扰的梯度计算,从而有利于获得准确的速度场模型参数的更新方向。另外在目标函数关于速度场模型参数的梯度计算表达式计算中,采用在不同的反演阶段使用不同的计算表达式,即在反演初期(即低频反演阶段)主要以相位信息反演为主,在反演的中后阶段,采用全波形反演,可以减弱反演对初始模型的依赖性。
综上,本发明通过多个单频编码信号同时激发实现同时源模拟;采用一阶速度—应力声波动方程,在一次正演模拟中完成多个震源子波的波场模拟,通过相敏检测方法分别对正传和反传波场实施波场解耦计算,求取无串扰的梯度结果,通过使用相位反演和全波形反演相结合的策略,可以降低目标函数的复杂度,从而减少反演中的多解性,实现高精度的反演。
附图说明
图1是本发明实施例的步骤示意图;
图2是Overthrust模型下真实速度场模型参数和初始速度场模型参数的示意图;
其中,图2a为真实的速度模型,图2b为反演所用的初始速度模型;
图3是本发明中采用相敏检测方法对反传源波场进行解耦得到的其对应的解耦波场的振幅和相位信息;
其中,图3a表示解耦波场的振幅信息,图3b表示解耦波场的相位信息;
图4是本发明实施例运用到Overthrust模型得到的反演结果与常规方法得到的反演结果的对比示意图;
其中,图4a是仅使用波场振幅信息的同时源反演结果;图4b是仅使用波场相位信息的同时源反演结果;图4c是同时使用波场振幅信息的同时源反演结果,即全波形反演方法;图4d是采用本发明实施例得到的结果;
图5是本发明实施例运用到Overthrust模型中得到的反演的目标函数下降曲线对比示意图。
具体实施方式
下面将结合本文实施例中的附图,对本文实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本文一部分实施例,而不是全部的实施例。基于本文中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本文保护的范围。
如图1所示,本发明实施例的反演方法,包括:
S1:通过多个位于不同位置的单频子波同时激发作为震源,采用多个检波点组成观测系统,通过观测系统获得观测数据;建立已知的速度场模型(采用Overthrust模型),并根据震源情况计算得出模拟数据;
S2:建立表征观测数据与模拟数据之间数据残差的目标函数;
所述目标函数采用最小化位于检波点Xr的观测数据dobs(Xr,Ω,t)和模拟数据u(Xr,Ω,t;m)的残差,具体为:
其中,x表示计算区域,m(x)表示计算区域x处的速度场模型参数;Xr为检波点的位置信息;Ω表示角频率,t表示时间;m表示速度场模型参数,u(Xr,Ω,t;m)表示在当前速度场模型参数m下获得的模拟数据。
S3:根据震源与检波点之间的位置,得出虚拟震源的信息,然后根据观测数据获取震源的正传波场;并根据目标函数获得虚拟震源的反传波场;
给定一系列位于Xs的单频随机信号s(Xs,Ω),则所述震源可表达为:
相应地,伴随方程为:
通常情况下,由于观测数据和模拟数据的观测系统不一致,所以需要分别从观测和模拟数据对应的单个信号数据中提取相同频率的信号,在观测数据的观测系统约束下求取数据的残差,以此解决观测系统不一致的问题。
S4:使用正传波场和反传波场进行梯度计算,迭代计算目标函数数据残差以更新速度场模型参数:
具体地,包括如下步骤:
S41:获得目标函数对速度场模型参数的梯度计算表达式;
S42:根据同时源模拟进行正演和反传计算,获得解耦的波场变量;
本实施例中,梯度计算表达式为正传波场与反传波场的特定形式的互相关,可简写为:
梯度简写形式G(x,m)表达式中,ψs(x,Ω)和ψr(x,Ω)分别表示正传波场与反传波场。从而,由(6)式的第二个等式可以看到,等式右端的第一项表示每一个编码信号对应波场数据的自相关运算;第二项表示当前编码信号与其它编码信号的波场互相关运算,即串扰。
由(6)式可知,只有从混叠的波场中解耦出所有编码信号的波场,然后只计算当前编码信号对应的正反传波场的自相关,才能从根本上消除串扰噪声。
即,有效地解耦混叠波场是解决串扰噪声的关键,本实施例引入相敏检测(PSD)法来解决这个问题。PSD的主要思想是:利用两个与未知信号具有相同频率但彼此相位相差90°的参考信号进行时间积分来提取振幅为sin(ωit)、相位为cos(ωit)的未知信号(即目标信号)。
本实施例中,取Erefsin(ωit+θref)为第一参考信号,取Erefcos(ωit+θref)为第二参考信号;其中,Eref为第一参考信号和第二参考信号的振幅;θref为所述第一参考信号和第二参考信号的相位。由此可以看出,第一参考信号的振幅和频率均与第二参考信号的振幅和频率相同,第一参考信号的相位与第二参考信号的相位相差90°。
由于空间某一点的混叠波场,均可以表示为:
使|ωi-ωj|≥pΔω(i,j∈ns,p∈Z),且Δω=2π/T。这里的T是地震记录的长度,p是最小的整数,pΔω表示编码信号之间的最小频率间隔。当定义最小的积分区间为Tp=T/p,对于任何q(p∈[1,p])倍的Tp区间qTp都作为合理的积分区间。通常,需要一定的时间等待波场达到稳态状况,所以p比q要小。
在积分区间qTp上,将所述第一参考信号与混叠波场的信号相乘并积分,得到:
Tp为最小积分区间;
以及,将所述第二参考信号与所述混叠波场的信号相乘并积分,得到:
从而,与第一参考信号和第二参考信号频率相同的编码信号的对应波场便能解耦出来,其计算公式为:
上述积分过程主要利用了三角函数正交性,只有与第一参考信号和第二参考信号频率相同的信号其得到积分结果不为零,其他与第一参考信号和第二参考信号频率不同的信号其得到的积分结果为零,从而实现对空间中任意一混叠波场中各频率编码信号其对应波场的解耦。
S43:利用目标函数对速度场模型参数的梯度计算表达式和解耦的波场变量,获取模型参数梯度;
其中速度模型的梯度计算式为:
其中,Re表示取向量的实数部分;u(x,ωi,t;m)为正传波场中,频率为ωi的编码信号解耦得到的其对应的波场变量,u′(x,ωi,t;m)为反传波场中,频率为ωi的编码信号经解耦得到其对应的波场变量。式(11)即是对频率均为ωi的两个波场变量进行正相关运算,从而实现无串扰的梯度求解。
本实施例采用五个频带的多尺度反演,在前两个频带反演(即反演的初始阶段),采用具有相位信息反演的梯度表达式,在后面三个频带反演(即反演的中后阶段),采用兼顾波场相位和振幅信息的全波形反演的梯度表达式。
基于上述梯度表达式,拓展出具有相位信息反演的梯度表达式为:
与之对应的兼顾波场相位和振幅信息的全波形反演的梯度表达式为
S5:采用共轭梯度方法求取更新方向,对模型参数梯度进行迭代更新,直至残差达到设定值或者迭代更新次数达到预设值,停止迭代更新输出速度场模型参数,此时采用该参数的速度场模型进行全波形反演。
采用共轭梯度方法求取更新方向Δmk+1,
式(14)中的下标k表示迭代次数,sk是用于计算共轭方向的一个中间变量,上标T表示矩阵转置,更新后的速度场模型为:
mk+1=mk+tkΔmk+1 (15)
即将式(14)中计算得到的Δmk+1带入到式(15)中,式(15)中的tk表示选取的更新步长,tkΔmk+1即为此次迭代过程中速度场模型参数的更新量。
上述反演方法所需的反演装置,包括:
目标函数构建模块,用于基于观测数据和模拟数据,构建表征所述观测数据与所述模拟数据之间数据残差的目标函数;
正传波场获取模块,用于获取震源的正传波场;
反传波场获取模块,用于获取虚拟震源的反传波场;
梯度计算模块,用于根据正传波场、反传波场和目标函数关于速度场模型参数的相位反演、全波场反演梯度计算表达式,计算速度场模型参数的梯度;
迭代更新模块,用于根据计算获得的梯度对当前速度场模型参数进行迭代更新,直至所述数据残差达到设定阈值或者迭代更新次数达到预设值,确定速度场模型参数。
效果验证:
将本发明实施例运用到Overthrust模型得到的反演结果与常规方法得到的反演结果进行对比;由图4b与图4a对比,能看出相位反演相对于振幅反演的精度更高,说明相位反演对初始模型的依赖性较小,可以恢复出有效的低波数信息;图4c与图4d相比,可以看出,采用本发明实施例全波形反演方法其结果比常规反演在界面的清晰度,高波数的连续性等方面更好。
如图5所示,采用本发明实施例的全波形反演方法与常规方法对应的目标函数下降曲线对比,由图可知,相比于其他几种反演方法,振幅反演的目标函数在低频下降缓慢,这说明振幅反演过早陷入局部极值;虽然同时使用相位和振幅的全波形反演在低频有着最大的下降率,但是随着迭代次数的增加,其目标函数的下降曲线逐渐与相位反演相一致;对于本发明实施例的方案,在开始阶段与相位反演吻合度很高,随着迭代次数的增加,其目标函数值逐渐低于相位反演,直至最后的频带、组合反演的目标函数值是四种反演算法中最低的,这说明本发明实施例的反演方法可以降低目标函数的复杂度,使得目标函数下降的更多。
Claims (5)
1.一种振幅、相位信息可调节的全波形反演方法,其特征在于,该方法的具体步骤为:
A、通过多个位于不同位置的单频子波同时激发作为震源,采用多个检波点组成观测系统,通过观测系统获得观测数据;建立速度场模型,并根据震源情况计算得出模拟数据;
B、基于观测数据和模拟数据,构建表征所述观测数据与所述模拟数据之间数据残差的目标函数;
C、根据震源与检波点之间的位置,得出虚拟震源的信息,然后根据观测数据获取震源的正传波场;并根据目标函数获得虚拟震源的反传波场;
D、确定目标函数关于速度场模型参数的梯度计算表达式,在不同的反演阶段具有两种计算表达式,其中在反演初期阶段,采用具有相位信息反演的梯度表达式,在反演的中后阶段,采用兼顾波场相位和振幅信息的全波形反演的梯度表达式;先对震源的正传波场和虚拟震源的反传波场分别进行无串扰解耦,获得解耦的波场变量;然后根据得到的波场变量和梯度计算表达式,计算得到速度场模型参数的梯度;其中无串扰解耦具体为:
采用相敏检测法先获取第一参考信号和第二参考信号,所述第一参考信号的振幅和频率均与所述第二参考信号的振幅和频率相同,所述第一参考信号的相位与所述第二参考信号的相位相差90°;然后根据第一参考信号和第二参考信号,采用三角函数正交性对震源的正传波场和虚拟震源的反传波场分别进行无串扰解耦;
基于解耦的波场,构建具有相位反演的梯度表达式
兼顾波场相位和振幅信息的梯度表达式为
E、采用共轭梯度方法求取更新方向,根据得到的梯度对速度场模型参数进行迭代更新,直至数据残差达到设定阈值或迭代更新次数达到预设值,停止迭代更新确定速度场模型参数,此时采用该参数的速度场模型进行全波形反演。
5.一种实现权利要求1至4任一项所述振幅、相位信息可调节的全波形反演方法的反演装置,其特征在于,包括:
目标函数构建模块,用于基于观测数据和模拟数据,构建表征所述观测数据与所述模拟数据之间数据残差的目标函数;
正传波场获取模块,用于获取震源的正传波场;
反传波场获取模块,用于获取虚拟震源的反传波场;
梯度计算模块,用于根据正传波场、反传波场和目标函数关于速度场模型参数的相位反演、全波场反演梯度计算表达式,计算速度场模型参数的梯度;
迭代更新模块,用于根据计算获得的梯度对当前速度场模型参数进行迭代更新,直至所述数据残差达到设定阈值或者迭代更新次数达到预设值,确定速度场模型参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110489838.1A CN113050163B (zh) | 2021-05-06 | 2021-05-06 | 一种振幅、相位信息可调节的全波形反演方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110489838.1A CN113050163B (zh) | 2021-05-06 | 2021-05-06 | 一种振幅、相位信息可调节的全波形反演方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113050163A CN113050163A (zh) | 2021-06-29 |
CN113050163B true CN113050163B (zh) | 2022-05-03 |
Family
ID=76518425
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110489838.1A Active CN113050163B (zh) | 2021-05-06 | 2021-05-06 | 一种振幅、相位信息可调节的全波形反演方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113050163B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116819602B (zh) * | 2023-07-12 | 2024-02-09 | 中国矿业大学 | 一种深度学习优化的变密度声波方程全波形反演方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353409A (zh) * | 2014-08-21 | 2016-02-24 | 中国石油化工股份有限公司 | 一种用于抑制全波形反演震源编码串扰噪音的方法和系统 |
CN105549079A (zh) * | 2016-01-12 | 2016-05-04 | 中国矿业大学(北京) | 一种地球物理参数的全波形反演模型的建立方法和装置 |
CN107843925A (zh) * | 2017-09-29 | 2018-03-27 | 中国石油化工股份有限公司 | 一种基于修正相位的反射波波形反演方法 |
CN110441816A (zh) * | 2019-09-20 | 2019-11-12 | 中国科学院测量与地球物理研究所 | 不依赖子波的无串扰多震源全波形反演方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9910174B2 (en) * | 2014-07-25 | 2018-03-06 | Seoul National University R&Db Foundation | Seismic imaging apparatus and method for performing iterative application of direct waveform inversion |
-
2021
- 2021-05-06 CN CN202110489838.1A patent/CN113050163B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105353409A (zh) * | 2014-08-21 | 2016-02-24 | 中国石油化工股份有限公司 | 一种用于抑制全波形反演震源编码串扰噪音的方法和系统 |
CN105549079A (zh) * | 2016-01-12 | 2016-05-04 | 中国矿业大学(北京) | 一种地球物理参数的全波形反演模型的建立方法和装置 |
CN107843925A (zh) * | 2017-09-29 | 2018-03-27 | 中国石油化工股份有限公司 | 一种基于修正相位的反射波波形反演方法 |
CN110441816A (zh) * | 2019-09-20 | 2019-11-12 | 中国科学院测量与地球物理研究所 | 不依赖子波的无串扰多震源全波形反演方法及装置 |
Non-Patent Citations (2)
Title |
---|
基于分频编码的弹性波全波形反演;劭祥奇等;《石油地球物理勘探》;20201231;第55卷(第6期);第1321-1329页 * |
瞬时相位信息用于时域弹性波全波形反演初始模型建立;罗静蕊等;《地球物理学进展》;20181231;第33卷(第6期);第2435-2440页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113050163A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2690373C (en) | Method for velocity analysis using waveform inversion in laplace domain for geophysical imaging | |
US11828894B2 (en) | Multi-scale unsupervised seismic velocity inversion method based on autoencoder for observation data | |
CN113156493B (zh) | 一种使用归一化震源的时频域全波形反演方法及装置 | |
CN106324662B (zh) | 一种针对目标层的全波形反演方法及系统 | |
CN103424777A (zh) | 一种提高地震成像分辨率的方法 | |
CN109541681A (zh) | 一种拖缆地震数据和少量obs数据联合的波形反演方法 | |
CN105242328B (zh) | 古热岩石圈厚度的确定方法及装置 | |
CN105093278A (zh) | 基于激发主能量优化算法的全波形反演梯度算子的提取方法 | |
CN113050163B (zh) | 一种振幅、相位信息可调节的全波形反演方法及装置 | |
CN112462427B (zh) | 多分量地震资料保幅角度域共成像点道集提取方法及系统 | |
CN113050160B (zh) | 一种数据阻尼全波形反演方法、装置和计算机设备 | |
GB2503640A (en) | Quality Assurance in a Full Waveform Inversion Process | |
WO2022153984A1 (ja) | 学習データ生成方法、モデル生成方法および学習データ生成装置 | |
CN114966831A (zh) | 一种基于速度-衰减解耦的粘声全波形反演方法 | |
CN111208568B (zh) | 一种时间域多尺度全波形反演方法及系统 | |
CN115373022A (zh) | 一种基于振幅相位校正的弹性波场Helmholtz分解方法 | |
CN111123373A (zh) | 一种基于波场扩展重构的探地雷达全波形反演方法 | |
CN110764146A (zh) | 一种基于声波算子的空间互相关弹性波反射波形反演方法 | |
EP3830608A1 (en) | Fast plane-wave reverse time migration | |
CN114488286B (zh) | 基于振幅加权的拖缆与海底地震资料联合波形反演方法 | |
CN118244355B (zh) | 一种基于重建观测地震数据的反射波走时反演方法 | |
CN115480306B (zh) | 一种基于深度学习的隧道全波形反演优化方法及系统 | |
CN111722287B (zh) | 一种基于渐进数据同化方法的震相特征识别波形反演方法 | |
CN118226510A (zh) | 一种高分辨率的多尺度波动方程反演方法及系统 | |
CN118033743A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |