CN110764146B - 一种基于声波算子的空间互相关弹性波反射波形反演方法 - Google Patents
一种基于声波算子的空间互相关弹性波反射波形反演方法 Download PDFInfo
- Publication number
- CN110764146B CN110764146B CN201911016597.8A CN201911016597A CN110764146B CN 110764146 B CN110764146 B CN 110764146B CN 201911016597 A CN201911016597 A CN 201911016597A CN 110764146 B CN110764146 B CN 110764146B
- Authority
- CN
- China
- Prior art keywords
- wave
- calculating
- formula
- longitudinal
- wave 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000003384 imaging method Methods 0.000 claims abstract description 23
- 238000005314 correlation function Methods 0.000 claims abstract description 9
- 238000013508 migration Methods 0.000 claims description 26
- 230000005540 biological transmission Effects 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 12
- 230000005012 migration Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 150000001875 compounds Chemical class 0.000 claims description 4
- 238000001615 p wave Methods 0.000 claims description 3
- 230000001902 propagating effect Effects 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 230000001808 coupling effect Effects 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 14
- 230000002596 correlated effect Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000010354 integration Effects 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
- 230000002194 synthesizing effect 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/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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
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)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提出了一种基于声波算子的空间互相关弹性波反射波形反演方法,采用空间相关函数实现对多波多分量反射地震数据运动学信息的有效利用,降低在数据匹配过程中对完整地震波形的高度依赖,缓解反演的非线性问题,从而准确、高效地反演纵、横波速度的背景模型。在实现过程中,本发明利用声波算子分别构建反射PP和PS的波路径,进行纵、横波速度的分步反演,这种分布反演策略可以有效压制纵、横波速度之间的耦合效应,提高两种参数的反演精度,同时显著降低算法对计算量和内存的需求。本发明得到的纵、横波速度反演结果可以为高分辨率、高精度的地震成像提供了准确的背景模型,因此对解决实际多分量地震数据的处理需求具有重要意义。
Description
技术领域
本发明涉及一种弹性波多参数背景建模方法,属于地震波成像技术领域。
背景技术
弹性波全波形反演方法(EFWI)可以利用观测的多分量地震数据获取地下弹性参数模型,如纵、横波速度、密度、弹性阻抗、各向异性参数等,是探索地球内部构造、进行油气资源勘探开发等重要手段。然而,受实际观测孔径有限的影响,EFWI通常只能从观测的初至波中获取浅层速度模型的有效更新,而有关中深层模型信息则体现在反射波中,却很难利用。为此,有学者提出了弹性波反射波形反演方法(ERWI)。该方法通过提取较低波数的反射波路径,可以实现多种参数模型的中深层更新,是目前研究的热点。
然而,基于波形残差目标函数的ERWI仍然具有高度的反演非线性,导致在匹配反射波波形时出现周波跳跃现象,随着偏移距的增大,这种现象更为明显,从而产生错误的反演结果。另外,反射波形的精确模拟依赖于真振幅成像技术。事实上,这种保真成像本身就是目前成像问题中的一个巨大挑战。相较于利用反射数据的波形信息,反射数据的运动学信息更为稳健。其中,空间相关函数可以有效地衡量反射波的空间聚焦程度,可以很好地描述反射波的运动学特征,因此可以有效降低反演的非线性,同时规避对真振幅偏移的依赖。另外,弹性波场通常在时间序列上具有较长的延续,采用空间相关函数可以提高运动学信息提取的准确性,因而在复杂介质参数反演中具有更好的适应性。
这种空间相关目标函数目前仅在单分量声波RWI问题中取得了较好的实际应用案例,尚未应用于多分量多参数的弹性波参数反演问题中。值得注意的是,相对于单参数声波RWI反演问题,ERWI还需要考虑更为严重的参数耦合效应,同时负担指数增长的内存需求、计算代价,因此EFWI中常用的Hessian预条件、波场分解等技术并不适用于ERWI问题。可见,发展基于反射波运动学信息的弹性波反射波形反演方法,同时设计更为合理的多参数反演策略,对实现弹性参数背景建模至关重要。
发明内容
针对现有ERWI方法存在的周波跳跃问题和对真振幅偏移技术的依赖等问题,本发明引入空间相关函数,以充分利用地震反射数据的运动学信息,缓解周波跳跃问题和对真振幅偏移的过度依赖,获取纵横波速度的背景模型。同时,利用采用分步反演策略,基于声波算子分别构建反射PP和PS波路径,实现纵横波速度的有效解耦,提高参数反演精度,降低算法对计算量和内存的需求。
本发明的具体技术方案如下:
一种基于声波算子的空间互相关弹性波反射波形反演方法,该方法具体包括以下步骤:
(1)读取的纵、横波速度初始模型α0和β0以及地震观测的PP波和PS波数据dp,obs,ds,obs;
(2)利用声波算子基于在纵、横波速度的当前模型αk和βk下计算炮端正传P波场up,b;同时,反传多分量地震观测数据dp,obs,ds,obs,计算检端反传P波场wp,b和S波场ws,b;其中本申请涉及到的所有正、反传波场均是时间t的函数;
(3)利用步骤2中得到的炮、检端波场,获取PP波的成像结果IpiIpp=∫up,b·ψp,bdt,此处up,b和ψi,b均是时间t的函数;
同时,基于炮端正传P波场up,b和PP波的成像结果Ipp,分别获取炮端反偏移P波场up,r和反偏移S波场us,r,并在检波点处抽取合成PP波和PS波数据dp,cal,ds,cal;
(4)基于地震观测数据和合成数据以及空间相关函数,求取空间相关道集,其计算公式为Cp,h(t,λ)=∫dp,obs(t,h+λ)·dp,cal(t,h)dh(此处均以PP波为例,PS波也同样计算),其中,变量t代表时间序列,变量h代表观测道号,变量λ代表观测道的偏移量,h取值范围一般设置为[-50,+50];
根据上述公式得到PP波的二维空间相关剖面Cp,h(t,λ);对应t中每一个时刻,记录Cp,h(t,λ)取最大值时的空间偏移量λ,生成所有时刻的Δp(t),并求取t中每一时刻上,检端反传的震源函数,
同理,计算PS波的相关道集Cs,h(t,λ),最大相关性对应的偏移量Δs(t),以及PS波的伴随震源gs;
然后,根据空间相关目标泛函计算目标函数值,其中,s代表观测数据的炮信息;根据目标函数值,判断当前模型下反射波的空间聚焦性是否满足收敛准则,若不满足,进入步骤(5),迭代更新纵波速度和横波速度的当前模型,直至目标函数满足收敛准则,进入步骤(9);
(5)计算地下炮端正传波场up,b及其反偏移波场up,r,us,r,并根据步骤(4)中获得的PP波和PS波的伴随震源gp和gs分别计算检端反传波场ψp,r,ψs,b及其反偏移波场ψp,r;
(6)求取炮、检端波场的零延迟互相关,并沿时间积分,即构建单炮情况下反射PP波和PS波的波路径,
Pathpp=∫(up,b·ψp,r+up,r·ψp,b)dt
Pathps=∫(us,r·ψs,b)dt
(7)将多炮反射PP、PS波的波路径进行累加,计算纵波速度和横波速度的梯度方向gα和gβ:
(8)根据纵波和横波速度的梯度利用给定的步长,对纵波和横波速度的当前模型αk和βk进行更新,
αk+1=αk+γαgα
βk+1=βk+γβgβ
最后,返回步骤(2)进行下一次迭代,k=k+1;
(9)输出纵、横波速度的当前模型αk和βk,作为最终的反演结果。
优选地,步骤(2)中,计算炮端正传P波场up,b公式如下:
计算检端反传P波场wp,b和S波场ws,b公式如下:
其中,dp代表P波的数据残差,ds代表S波的数据残差,符号*代表共轭算子。由于本申请采用空间相关值作为衡量标准,因而对地震数据的振幅和相位差具有较高的容忍度,因此使用声波算子实现P波和S波的快速计算,节省反射波形反演中对内存和计算量的过度需求。
优选地,步骤(3)中,获取炮端反偏移P波场up,r和反偏移S波场us,r的公式如下:
优选地,步骤(5)中,
a)对于纵波速度的当前模型更新,计算炮端正传P波场up,b的公式为:
计算炮端反偏移波P波场up,r的公式为:
计算检端反传P波场ψp,b的公式为:
检端反偏移P波场ψp,r的公式为:
b)对于横波速度的当前模型更新,需要计算炮端反偏移S波场us,r的公式为:
和计算检端反传S波场ψs,b的公式为:
由于PS像通常具有更多的偏移假象(耦合脚印和非物理波现象),同时需要解决极性反转问题,因此,在计算反偏移的S波场时,本申请仍然采用PP波的像进行替代。
本发明首先利用观测的多分量地震数据和纵横波速度的初始模型获取地下界面信息;然后,基于成像结果利用反偏移技术预测反射数据,并根据空间相关函数估计运动学差异,求取反射PP波和PS波的伴随震源;最后,利用构建的反射PP和PS波路径合成梯度,更新参数模型。为了尽可能地压制纵横波速度之间的耦合效应,本发明需要优先反演纵波速度模型,当纵波速度较为准确后,再进行横波速度背景反演。另外,由于空间相关目标函数对地震数据的振幅、相位精度具有较高的容忍度,因此,本发明在计算弹性波场时均采用声波算子,已实现多分量地震数据的快速计算,以显著降低内存和计算量,同时便于分步反演策略的有序开展。
本发明相比现有技术具有如下有益效果:
本发明提的基于声波算子的空间互相关弹性波反射波形反演方法,采用空间相关函数突出对多分量反射数据运动学信息的有效利用,有效降低反演的非线性问题,从而准确、稳健地反演纵横波速度的背景模型。
在实现过程中,本发明采用声波算子分别构建反射PP和PS的波路径,实现纵、横波速度的分步反演,这种分布反演策略可以有效压制纵横波速度之间的耦合效应,提高参数的反演精度,同时显著降低算法对计算量和内存的需求。
本发明得到的反演结果可以为高分辨率地震成像和高精度的全波形反演提供了必要的参数建模,因此对实现参数高精度建模的实际应用具有重要意义。
附图说明
图1是本发明流程图(也作摘要附图);
图2是Sigbee2A实验中真实纵(a)、横波(b)速度模型;
图3是Sigbee2A实验中初始纵(a)、横波(b)速度模型;
图4是地震观测反射数据水平(a)和垂直(b)分量;
图5是初始模型下得到的PP波(a)和PS波(b)成像结果;
图6是真实模型下得到的PP波(a)和PS波(b)成像结果;
图7是初始模型得到的像生成的反射PP波(a)和PS波(b)数据;
图8是初始模型下得到的反射PP波(a)和PS波(b)空间相关道集;
图9是单炮情况下的反射PP波炮(a)、检(b)端反射波路径;
图10是反演得到的纵(a)、横波(b)速度模型;
图11是反演得到的PP波(a)和PS波(b)空间相关道集;
图12是反演结果作为偏移速度得到的PP波(a)和PS波(b)成像结果;
图13是反演结果作为偏移速度得到的PP波(a)和PS波(b)空间相关道集;
图14是反演结果作为EFWI初始模型得到的纵(a)、横波(b)速度反演结果;
图15是水平位置2km处纵(a)、横波(b)速度模型的垂向抽线结果;
图16是初始模型作为EFWI初始模型得到的纵(a)、横波(b)速度反演结果。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步说明,
实施例:
本发明基于声波算子的空间互相关弹性波反射波形反演方法,下面以Sigbee2A模型实验为例,其真实模型如图2所示,该方法具体包括以下步骤:
(1)读取的纵、横波速度初始模型α0和β0,如图3所示,以及多分量地震观测,包括反射PP波和PS波dp,obs,ds,obs,如图4所示;
(2)根据声波算子以及纵、横波速度的当前模型αk和βk,计算炮端正传P波场up,b;根据观测的反射PP波数据dp,obs,计算检端反传P波场wp,b。本申请涉及到的所有正、反传波场均是时间t的函数;
(3)根据炮端正传P波场up,b和检端反传P波场wp,b,以及成像公式Ipp=∫up,b·wp, bdt,获取PP波成像结果,如图5a所示,为了同时反应纵横波速度模型的准确性,图5也展示了反射PS波的成像结果,但在算法中仍然采用PP波成像结果来合成反射PS波;
根据真实模型得到的PP波成像结果,如图6所示。对比可以看出,在PP波的成像结果中,由于初始模型距离真实模型较远,主要层位定位不准确,同时底部散射点无法收敛。
根据炮端正传P波场up,b和PP波的成像结果Ipp,计算炮端反偏移P波场up,r和反偏移S波场us,r,并在检波点处抽取反偏移得到的合成反射PP波和PS波数据dp,cal,ds,cal,如图7所示;由于反射PS波数据在零偏位置处存在极性转换,这是声波算子无法正确模拟的,因此在实际匹配之前需要对单炮PS数据在近偏移距进行置零处理);
(4)根据地震观测P波数据和合成P波数据以及空间相关函数,其计算公式为Cp,h(t,λ)=∫dp,obs(t,h+λ)·dp,cal(t,h)dh,计算反射PP波空间相关道集,如图8所示,其中,变量t代表时间序列,变量h代表观测道号,变量λ代表空间偏移量,λ取值一般为[-50,+50]内的所有整数;
根据上述公式得到PP波的二维空间相关剖面Cp,h(t,λ);对应t中每一个时刻,记录Cp,h(t,λ)取最大值时的空间偏移量λ,生成所有时刻的Δp(t),并求取t中每一时刻上,检端反传的震源函数,
根据地震观测S波数据和合成S波数据,计算反射PS波的相关道集Cs,h(t,λ),最大相关性对应的偏移量Δs(t),以及PS波的伴随震源gs;
根据目标函数值,判断当前模型下反射波的空间聚焦性是否满足收敛准则,如:是否小于初始目标函数的1%,若不满足,进入步骤(5);若满足,进入步骤(9);
(5)根据纵横波速度当前模型,计算炮端正传波场up,b及其反偏移波场up,r,us,r,并根据反射PP波和PS波的伴随震源gp和gs计算检端反传波场ψp,r,ψs,b及其反偏移波场ψp,r;
(6)根据炮端正传波场up,b及其反偏移波场up,r,检端反传波场ψp,r,ψs,b及其反偏移波场ψp,r,计算单炮情况下反射PP波和PS波的波路径,如图9所示,其公式如下
Pathpp=∫(up,b·ψp,r+up,r·ψp,b)dt
Pathps=∫(us,r·ψs,b)dt
(7)根据单炮情况下反射PP波和PS波的波路径,计算计算纵波速度和横波速度的梯度gα和gβ
(8)根据纵波和横波速度的梯度利用给定的步长,对纵波和横波速度的当前模型αk和βk进行更新,
αk+1=αk+γαgα
βk+1=βk+γβgβ
最后,返回步骤(2)进行下一次迭代,k=k+1;
(9)输出纵、横波速度的当前模型αk和βk,作为最终的反演结果,如图10所示。
步骤(2)中,计算炮端正传P波场up,b公式如下:
计算检端反传P波场wp,b和S波场ws,b公式如下:
其中,dp代表P波的数据残差,ds代表S波的数据残差,符号*代表共轭算子。需要注意的是,由于本申请采用空间相关值作为衡量标准,因而对地震数据的振幅和相位差具有较高的容忍度,因此使用声波算子实现P波和S波的快速计算,节省反射波形反演中对内存和计算量的过度需求。
步骤(3)中,计算炮端反偏移P波场up,r和反偏移S波场us,r的公式如下:
步骤(5)中,
a)计算炮端正传P波场up,b的公式为:
计算炮端反偏移波P波场up,r的公式为:
计算检端反传P波场ψp,b的公式为:
检端反偏移P波场ψp,r的公式为:
b)计算炮端反偏移S波场us,r的公式为:
和计算检端反传S波场ψs,b的公式为:
从空间相关道集(如图11所示)可以看出,经过SC-ERWI方法反演后得到的空间聚焦性得到改善,偏离中心位置的程度明显降低。同时,基于反演得到的纵横波速度模型,分别求取PP波和PS波的像(如图12所示),并提取共成像点道集(如图13所示)。可见,经过多次迭代后,界面位置更加准确,底部散射点收敛性更好,且与真实模型得到的像非常接近,这说明了利用SC-ERWI方法反演得到的纵横波速度模型已经可以较好地描述参数的背景,使得预测的反射波可以很好的匹配观测数据的运动学特征。
同时,可以将本发明反演得到的纵横波速度模型作为高精度EFWI的初始模型,得到的反演结果如图14所示,其垂直抽线如图15所示。从图中可以看出,本发明所研究的方法可以准确地恢复纵横波速度的宏观背景,再经过EFWI反演得到的纵横波速度模型则非常准确,可以很好地拟合真实模型,而如果以图2模型作为EFWI的初始模型,EFWI则会陷入局部极值,得到错误的反演结果(如图16所示)。可见,本发明方法,可以正确利用反射数据的运动学信息,恢复纵横波速度的背景模型,可以满足EFWI对初始模型精度的需求。
最后,需要指出的是,由于利用声波算子构建反射波路径来替代常规ERWI中使用的标准弹性波算子,因此,在内存上,本发明研究方法只需要存储压力分量,而ERWI至少需要存储水平和位移压力分量,因此在二维情况下内存上可以节省50%(三维情况下为66.7%)。而在计算量方面,求解两次声波方程和弹性波方程并实现波场分解所需的计算量比约为4:12,因此SC-ERWI方法在计算量上可以节省67%(如表1所示)。
表1 SC-ERWI和ERWI方法资源消耗对比
Claims (4)
1.一种基于声波算子的空间互相关弹性波反射波形反演方法,该方法具体包括以下步骤:
(1)读取的纵、横波速度初始模型α0和β0以及多分量地震观测,包括反射PP波数据和PS波数据dp,obs,ds,obs;
(2)根据声波算子以及纵、横波速度的当前模型αk和βk,计算炮端正传P波场up,b;根据观测的反射PP波数据dp,obs,计算检端反传P波场wp,b,其中正、反传波场均是时间t的函数;
(3)根据炮端正传P波场up,b和检端反传P波场wp,b,以及成像公式Ipp=∫up,b·wp,bdt,获取PP波的成像结果;
根据炮端正传P波场up,b和PP波的成像结果Ipp,计算炮端反偏移P波场up,r和反偏移S波场us,r,并在检波点处抽取反偏移得到的合成反射PP波和PS波数据dp,cal,ds,cal;
(4)根据地震观测PP波数据和合成PP波数据以及空间相关函数,其计算公式为Cp,h(t,λ)=∫dp,obs(t,h+λ)·dp,cal(t,h)dh,计算反射PP波的空间相关道集,其中,变量t代表时间序列,变量h代表观测道号,变量λ代表空间偏移量,λ取值一般为[-50,+50]内的所有整数;
根据上述公式得到PP波的二维空间相关剖面Cp,h(t,λ);对应t中每一个时刻,记录Cp,h(t,λ)取最大值时的空间偏移量λ,生成所有时刻的反射PP波空间偏移量时间序列Δp(t),并求取t中每一时刻上,检端反传的震源函数,获得单炮情况下PP波的伴随震源gp;
根据地震观测PS波数据和合成PS波数据,计算反射PS波的相关道集Cs,h(t,λ),最大相关性对应的偏移量Δs(t),以及PS波的伴随震源gs;
根据目标函数值,判断当前模型下反射波的空间聚焦性是否满足收敛准则,即:是否小于初始目标函数的1%,不满足,则进入步骤(5);若满足条件,进入步骤(9);
(5)根据纵横波速度当前模型,计算炮端正传P波场up,b、反偏移P波场up,r和S波场us,r,并根据反射PP波和PS波的伴随震源gp和gs计算检端反传P波场ψp,b、S波场ψs,b和反偏移P波场ψp,r;
(6)根据炮端正传P波场up,b,炮端反偏移P波场up,r,检端反传P波场ψp,b,检端反传S波场ψs,b和检端反偏移P波场ψp,r,计算单炮情况下反射PP波和PS波的波路径Pathpp和Pathps,其公式如下
Pathpp=∫(up,b·ψp,r+up,r·ψp,b)dt
Pathps=∫(us,r·ψs,b)dt
(7)根据单炮情况下反射PP波和PS波的波路径,计算纵波速度和横波速度的梯度gα和gβ:
(8)根据纵波和横波速度的梯度利用给定的步长,对纵波和横波速度的当前模型αk和βk进行更新,其中γα和γβ分别代表纵波和横波速度的更新步长;
αk+1=αk+γαgα
βk+1=βk+γβgβ
最后,返回步骤(2)进行下一次迭代,k=k+1;
(9)输出纵、横波速度的当前模型αk和βk,作为最终的反演结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911016597.8A CN110764146B (zh) | 2019-10-24 | 2019-10-24 | 一种基于声波算子的空间互相关弹性波反射波形反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911016597.8A CN110764146B (zh) | 2019-10-24 | 2019-10-24 | 一种基于声波算子的空间互相关弹性波反射波形反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110764146A CN110764146A (zh) | 2020-02-07 |
CN110764146B true CN110764146B (zh) | 2020-05-19 |
Family
ID=69333523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911016597.8A Active CN110764146B (zh) | 2019-10-24 | 2019-10-24 | 一种基于声波算子的空间互相关弹性波反射波形反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110764146B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111781639B (zh) * | 2020-06-04 | 2021-06-04 | 同济大学 | 针对obs多分量数据的炮检互易弹性波全波形反演方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104391323B (zh) * | 2014-11-21 | 2015-11-18 | 中国石油大学(华东) | 一种利用反射波信息反演速度场中低波数成分的方法 |
US10353092B2 (en) * | 2015-12-10 | 2019-07-16 | Pgs Geophysical As | Velocity model update with an inversion gradient |
US11487036B2 (en) * | 2017-01-12 | 2022-11-01 | Cgg Services Sas | Reflection full waveform inversion methods with density and velocity models updated separately |
CN107843925B (zh) * | 2017-09-29 | 2019-03-08 | 中国石油化工股份有限公司 | 一种基于修正相位的反射波波形反演方法 |
CN108873066B (zh) * | 2018-06-26 | 2020-03-10 | 中国石油大学(华东) | 弹性介质波动方程反射波旅行时反演方法 |
CN110187382B (zh) * | 2019-03-05 | 2020-10-13 | 中国石油大学(华东) | 一种回折波和反射波波动方程旅行时反演方法 |
-
2019
- 2019-10-24 CN CN201911016597.8A patent/CN110764146B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110764146A (zh) | 2020-02-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108139499B (zh) | Q-补偿的全波场反演 | |
CN108873066B (zh) | 弹性介质波动方程反射波旅行时反演方法 | |
US11041971B2 (en) | Full wavefield inversion with an image-gather-flatness constraint | |
CN111025387B (zh) | 一种页岩储层的叠前地震多参数反演方法 | |
CN109188519B (zh) | 一种极坐标下的弹性波纵横波速度反演系统及方法 | |
EA032186B1 (ru) | Сейсмическая адаптивная фокусировка | |
CN111025386B (zh) | 一种无分离假象的纵横波分离方法 | |
RU2570827C2 (ru) | Гибридный способ для полноволновой инверсии с использованием способа одновременных и последовательных источников | |
CN111580163B (zh) | 一种基于非单调搜索技术的全波形反演方法及系统 | |
US20150362608A1 (en) | Combined interpolation and primary estimation | |
WO2021116800A1 (en) | System and method for using a neural network to formulate an optimization problem | |
CN110764146B (zh) | 一种基于声波算子的空间互相关弹性波反射波形反演方法 | |
CN111257930B (zh) | 一种黏弹各向异性双相介质区域变网格求解算子 | |
CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
CN116520418A (zh) | 一种弹性波角度域共成像点道集高效提取方法 | |
CN117388944A (zh) | 地质模型约束的多物性参数反演方法 | |
CN111308550B (zh) | 一种页岩vti储层的各向异性参数多波联合反演方法 | |
CN115373022A (zh) | 一种基于振幅相位校正的弹性波场Helmholtz分解方法 | |
CN110857999B (zh) | 一种基于全波形反演的高精度波阻抗反演方法及系统 | |
CN113156498B (zh) | 一种基于同伦延拓的叠前avo三参数反演方法和系统 | |
CN113791447B (zh) | 一种反射结构导引的反射波层析反演方法 | |
CN112649848A (zh) | 利用波动方程求解地震波阻抗的方法和装置 | |
Xue et al. | Joint inversion of location, excitation time, and amplitude of microseismic sources | |
CN113589362B (zh) | 三维陆上耦合波正演模拟方法 | |
Fang et al. | Deep learning optimization using the gradient of a custom objective function: An FWI example study on the convolutional objective function |
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 |