共查询到20条相似文献，搜索用时 168 毫秒
1.
本文详细地给出了基于非均匀介质的体散射广义散射层析成像反演的基本理论.广义散射层析成像反演可以描述为波场的反向传播和对成像场进行局部波数域滤波的过程.在数值算例中,利用背景速度沿深度方向均匀变化的v(z)介质中的简单的方块作为速度异常体的模型,通过对该模型产生的低频的Born数据和声波的正演数据的测试,在对采集系统进行有限频率带宽和空间孔径的校正来进行局部成像矩阵谱的恢复中,可以看出模型中各点的谱在恢复后的质量无论从覆盖的面积范围还是幅值的均一性上都有着明显的提高;在对速度模型的重建中,广义散射层析成像反演能够很好地恢复速度模型的低频分量,即便是方块速度异常体相对于背景速度的平均速度扰动是23%也能很好地重建模型中的速度,且对于不同的背景速度模型基本上都能很好地恢复Marmousi速度模型的低频分量.所以该方法将基于Born模型的层析成像反演适应范围进行了一定程度的扩展. 相似文献
2.
《地球物理学进展》2017,(2)
传统的地震逆散射广义Radon变换(GRT)保幅反演方法是建立在散射场一阶Born近似(单散射)的基础上,仅仅适用于弱扰动介质模型.本文从散射场积分方程出发,通过研究二次散射的特征,讨论和验证了基于局部二阶Born近似的GRT非线性保幅反演方法,将传统GRT线性保幅反演算子的适用范围扩展至非均匀强扰动介质.数值测试结果表明:在散射场近似模拟方面,二阶Born近似比一阶Born近似更为准确,二次散射效应主要集中在主散射点周围的局部区域内,超过这一范围,二次散射强度趋于稳定;在保幅反演方面,本文基于局部二阶Born近似的GRT非线性反演算法,明显优于传统的GRT线性反演算法,可以准确重构强扰动介质模型,而计算效率与线性反演方法相当. 相似文献
3.
有限频率层析成像考虑了非均匀介质中波的散射、衍射、波前愈合等物理性质,使得其对速度异常体的分辨能力远大于射线层析成像.推导和计算有限频率敏感核是进行有限频率层析成像的关键,当前推导有限频率敏感核多借助一阶Born近似,但这只适用于弱散射介质的情况.本文基于二阶Born近似并利用傅里叶变换推导了三维均匀介质情况下有限频率敏感核的解析表达式,并将其推广到非均匀介质中得到了三维非均匀介质中有限频率敏感核.研究表明:当介质中速度扰动小于2%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核差别很小,可近似认为相同;当介质中速度扰动大于5%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核有较大不同,表明此时已不能忽略二次散射. 相似文献
4.
5.
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异. 相似文献
6.
《应用地球物理》2015,(3)
多分量地震数据中低频缺失是弹性波全波形反演中的一大难题,低频的缺失导致全波形反演无法有效恢复介质的长波长成分进而使反演陷入局部极值。为此,本文提出了反演介质纵横波速度长波长分量的弹性波包络反演方法。该方法利用包络算子具有的解调多分量数据中隐含的低频信息的能力,构造多分量地震数据的包络目标函数进行反演,用以恢复地下介质纵横波速度的长波长成分。一系列数值试验表明,即使在多分量地震数据中缺失低频信息、并且初始模型缺少先验信息的情况下,这种弹性波包络反演方法能够有效降低波形反演的非线性,可以为后续的常规弹性波全波形反演或者深度偏移提供足够精确的初始模型,且该方法对横波速度长波长分量的重建尤为有效。Mamousi-2模型的高精度纵横波速度的反演结果表明,利用该方法反演的纵横波速度作为常规弹性波全波形反演的初始模型,可以显著提高反演结果的精度。此外,本文对弹性波包络反演方法的适用性也进行了初步的研究与讨论。 相似文献
7.
地震波在非弹性介质中的衰减效应常用品质因子Q度量.相对准确的Q模型对提高强衰减介质中地震波成像的质量至关重要.本文提出了黏声介质反射波形反演(QRWI)方法来重建地下宏观Q模型.在缺乏大偏移距和低频地震数据时,该方法以黏声波方程为波场传播引擎,利用反射波核函数对模型中深部的敏感性去提取背景Q值.当速度高、低波数成分均已知时,基于波形拟合的QRWI可以获得较高分辨率的反演结果.由于地下介质速度的高波数扰动很难准确估计,本文通过引入峰值频移目标函数,极大地降低了QRWI对速度高波数成分的依赖.理论合成数据实验结果表明,本文方法反演得到的宏观Q模型可以满足衰减补偿逆时偏移成像的要求. 相似文献
8.
根据地震波传播的有限频理论,对于某个特定震相的观测信息,不仅射线路径上的点对该信息具有影响,射线领域上的其他点对接收信息也具有影响,这种影响可以用核函数来表达.本文基于波动方程的Born近似与Rytov近似,给出了非均匀介质情况下初至地震波振幅与走时菲涅尔体层析成像单频、带限层析核函数的计算方法.通过对均匀介质情况下初至地震波菲涅尔体层析成像核函数解析表达式的理论模型实验与分析,给出了不同维度振幅、走时单频菲涅尔体的空间分布范围,进而给出了带限菲涅尔体边界的确定方法.将本文的走时菲涅尔体层析成像理论应用于表层速度结构反演中,理论模型试验与实际资料处理结果表明,初至波菲涅尔体地震层析成像方法比传统的初至波射线层析成像理论具有更高的反演精度. 相似文献
9.
高频假设下的地震射线理论以及相应的地震成像理论表明,在射线稀疏条件下,不可能得到较高分辨率的构造成像;而有限频射线理论更符合实际地震的传播规律,即地震波的走时不仅与中心射线(传统的几何射线)上的速度分布有关,而且与中心射线附近一定范围(称其为第一菲涅耳体)内的速度异常分布有关.鉴于此,本文提出了计算多震相地震波菲涅耳体有限频射线的方法,并定义了走时敏感核函数,同时给出了利用多震相菲涅耳体有限频射线进行速度模型和反射界面同时反演成像的公式.利用多震相走时资料,使用传统射线层析成像方法与有限频射线层析成像方法进行了速度和界面的同时反演成像.结果表明,当射线密度较小时,无论是对速度模型的重建还是对反射界面几何形状的更新,有限频射线层析成像方法均优于传统射线层析成像方法, 而变频有限频射线层析成像则是实际地震层析成像的首选反演算法. 相似文献
10.
《中国科学:地球科学》2017,(10)
得益于结合最优化算法和完整的波动方程,全波形反演已成为地震勘探最前沿的研究方向,正逐渐发展成为获取地球内部信息的重要工具.但是,在初始速度模型不够精准的条件下,利用纯反射波进行波形反演来更新光滑的背景速度模型是非常困难的.本文展示了一种利用纯反射波通过交替进行最小二乘逆时偏移和波形反演的方法来成功地同时反演反射系数和背景速度的方法.其波形反演中用于更新背景速度的层析梯度是通过基于Born近似的正演模拟来分离反射波波场得到的,这样可以很好地更新光滑的背景速度模型.之后,常规全波形反演把上一步得到的光滑背景速度模型作为初始模型,利用纯反射波获得最终的高保真、高分辨率的速度模型.进一步分析揭示了高分辨率的速度模型是由于在全波形反演中应用了反褶积成像条件.相对于其他类型的地震波,反射波带有地下更深部的信息,因此,反射波波形反演能够进一步提高我们对地球内部的成像精度. 相似文献
11.
Migration velocity analysis is a method devoted to the evaluation of both reflectivity and background velocity models, associated with the high and low wavenumber components of the model, respectively. Inversion velocity analysis is one of its improved versions, leading to more stable background velocity updates. Still, the impact of the user parameters should be understood for an optimal update of the background velocity. We show that a sign reversal of the background velocity gradient could occur when the selected surface offset range or the space lag range is too small. We derive the theoretical limits and check their consistency through simulations in a simple model with a single interface. These guidelines determine the necessary ranges of surface offsets and space lags for a proper update of the background velocity model. We discuss their applicability on the Marmousi model. Artefacts in the retrieved background velocity model are observed when the guidelines are not satisfied. 相似文献
12.
A local plane-wave approach of generalized diffraction tomography in heterogeneous backgrounds, equivalent to Kirchhoff summation techniques when applied in seismic reflection, is re-programmed to act as repeated synthetic aperture radar (SAR) imaging for seismic prestack depth migration. Spotlight-mode SAR imaging quickly provides good images of the electromagnetic reflectivity of the ground via fast Fourier transform (FFT)-based signal processing. By calculating only the Green's functions connecting the aircraft to the centre of the illuminated patch, scattering structures around that centre are also recovered. SAR technology requires us to examine seismic imaging from the local point of view, where the quantity and quality of the available information at each image point are what are important, regardless of the survey geometry. When adapted to seismics, a local image of arbitrary size and sampling is obtained by FFT of seismic energy maps in the scattering wavenumber domain around each node of a pre-calculated grid of Green's functions. These local images can be used to generate a classic prestack depth-migrated section by collecting only their centres. However, the local images also provide valuable information around the centre, as in SAR. They can therefore help to pre-analyse prestack depth migration efficiently, and to perform velocity analysis at a very low cost. The FFT-based signal-processing approach allows local, efficient and automatic control of anti-aliasing, noise and resolution, including optimized Jacobian weights. Repeated local imaging could also be used to speed up migration, with interpolation between local images associated with a coarse grid of Green's functions, as an alternative to interpolation of Green's functions. The local images may, however, show distortions due to the local plane-wave approximation, and the velocity variations across their frame. Such effects, which are not necessarily a problem in SAR, should be controlled and corrected to further enhance seismic imaging. Applications to realistic models and to real data show that, despite the distortion effects, the local images can yield similar information to prestack depth migration, including common-image-point gathers for velocity analyses and AVO/AVA effects, at a much lower cost when a small target is considered. 相似文献
13.
Wave equation least square imaging using the local angular Hessian for amplitude correction 总被引：4，自引：0，他引：4
Local angular Hessian can be used to improve wave equation least square migration images. By decomposing the original Hessian operator into the local wavenumber domain or the local angle domain, the least square migration image is obtained as the solution of a linearized least‐squares inversion in the frequency and local angle domains. The local angular Hessian contains information about the acquisition geometry and the propagation effects based on the given velocity model. The inversion scheme based on the local angular Hessian avoids huge computation on the exact inverse Hessian matrix. To reduce the instability in the inversion, damping factors are introduced into the deconvolution filter in the local wavenumber domain and the local angle domain. The algorithms are tested using the SEG/EAGE salt2D model and the Sigsbee2A model. Results show improved image quality and amplitudes. 相似文献
14.
全波形反演方法是一种数据域高精度反演方法,该方法通过匹配观测数据与模拟数据的地震波形,利用梯度法准确反演地下介质参数的分布情况.由于观测数据普遍缺少低频信息,该方法易受周期跳跃现象影响.特别是当地下存在大尺度强反射界面的构造时,地下介质的反演转化为强非线性问题求解.该情形下,即使观测数据包含充足的低频信息,全波形反演也难以给出准确的反演结果.一般可以通过减弱反演对初始模型参数的依赖性来克服上述问题,具体表现为使用新变量(例如瞬时相位、包络等)代替目标函数中的采样后波场,以增强新目标函数的凸性.但是,对该新目标函数进行反演时,伴随状态方程中存在关于新变量和波场的一个链式微分项,该项保留了反演问题的非线性,导致新的反演方法难以处理包含大尺度构造的强非线性反演问题.此外,基于新变量的反演问题依然在波场空间中计算模型梯度,难以充分利用新变量与模型参数之间的弱非线性关系.因此,本文提出用频率域波动方程的相位形式代替传统的波动方程来消除伴随状态方程中的链式微分项,用解缠绕的相位代替目标函数中采样前波场并在相位空间进行反演.该方法可以最大程度地利用地下介质参数和解缠绕相位之间的弱非线性关系,从而削弱反演的非线性性.由于基于频率域波场计算得到相位有严重的缠绕问题,本文采用基于振幅排序的多聚类算法来对相位进行解缠绕.虽然将介质参数到波场的映射替换为介质参数与解缠绕相位的映射,会导致反演结果的分辨率有所下降,但该方法可以在相位空间恢复介质参数的大尺度低波数分量.Marmousi模型测试证明了该方法的有效性和准确性,针对部分BP模型的测试也证明了该方法处理强非线性问题的能力. 相似文献
15.
Liangguo Dong Zhongyi Fan Hongzhi Wang Benxin Chi Yuzhu Liu 《Geophysical Prospecting》2018,66(8):1503-1520
Reflection full waveform inversion can update subsurface velocity structure of the deeper part, but tends to get stuck in the local minima associated with the waveform misfit function. These local minima cause cycle skipping if the initial background velocity model is far from the true model. Since conventional reflection full waveform inversion using two‐way wave equation in time domain is computationally expensive and consumes a large amount of memory, we implement a correlation‐based reflection waveform inversion using one‐way wave equations to retrieve the background velocity. In this method, one‐way wave equations are used for the seismic wave forward modelling, migration/de‐migration and the gradient computation of objective function in frequency domain. Compared with the method using two‐way wave equation, the proposed method benefits from the lower computational cost of one‐way wave equations without significant accuracy reduction in the cases without steep dips. It also largely reduces the memory requirement by an order of magnitude than implementation using two‐way wave equation both for two‐ and three‐dimensional situations. Through numerical analysis, we also find that one‐way wave equations can better construct the low wavenumber reflection wavepath without producing high‐amplitude short‐wavelength components near the image points in the reflection full waveform inversion gradient. Synthetic test and real data application show that the proposed method efficiently updates the background velocity model. 相似文献
16.
将局部余弦基小波束波场分解、传播与观测系统沉降法叠前深度偏移相结合,推导了源-检波器观测系统沉降法传播算子.本算法中,先对频率域的共点源和共点检波器道集做局部余弦小波束分解,然后分别沿共小波束源和共小波束检波器在深度方向延拓得到下一层波场.每个深度的波场,都等效于把源和检波器放在该层后所能接收到的地震记录,每点的像值由炮点和检波点重合时的零时刻波场值给出.通过二维SEG/EAGE盐丘模型和Marmousi模型的偏移成像结果验证该方法理论推导的正确性.另外,结果显示该方法继承了小波束域波场延拓在速度扰动较大情况下波传播及成像精度高的优点. 相似文献
17.
Diffractions play a vital role in seismic processing as they can be utilized for high‐resolution imaging applications and analysis of subsurface medium properties like velocity. They are particularly valuable for anisotropic media as they inherently possess a wide range of dips necessary to resolve the angular dependence of velocity. However, until recently, the focus of diffraction imaging or inversion algorithms have been only on the isotropic approximation of the subsurface. Using diffracted waves, we develop a framework to invert for the effective η model. This effective model is obtained through scanning over possible effective η values and selecting the one that best fits the observed moveout curve for each diffractor location. The obtained effective η model is then converted to an interval η model using a Dix‐type inversion formula. The inversion methodology holds the potential to reconstruct the true η model with sufficiently high accuracy and resolution properties. However, it relies on an accurate estimation of diffractor locations, which in turn requires good knowledge of the background velocity model. We test the effectiveness and applicability of our method on the vertical transverse isotropic Marmousi model. The inversion results yield a reasonable match even for the complex Marmousi model. 相似文献
18.
Zhang Zhifu Meng Xiaohong Liu Hong Chen Jingbo Li YouMing 《应用地球物理》2006,3(4):218-225
An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure. 相似文献
19.
20.
层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行走时层析速度反演是当前比较成熟有效且广泛应用的技术.本文从高斯束偏移成像条件出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域走时扰动与速度扰动的线性关系,建立了成像域走时层析方程及其显式表达的层析核函数.该核函数的本质是有限频层析核函数,利用该核函数替换常规射线层析核函数可以明显提高层析反演精度.该核函数的计算关键是背景波场格林函数的计算,本文利用高斯束传播算子计算格林函数进而得到走时层析核函数,实现方式灵活高效且计算精度较高.基于高斯束传播算子的偏移成像与层析成像相结合进行深度域建模迭代,体现了速度建模与偏移成像一体化的思想.数值计算及实际数据应用证明了基于高斯束传播算子的成像域走时层析方法的有效性. 相似文献