共查询到17条相似文献，搜索用时 156 毫秒
1.
根据瞬变电磁场所满足的扩散方程与波动场所满足的波动方程之间的数学关系,运用扫时波场变换技术,得到稳定的虚拟波场.在虚拟波场的基础上,借用地震勘探中的Born近似逆散射技术实现对电性界面进行成像的目的.在理论推导中,直接从时间域波动方程出发,运用格林定理、引入均匀半空间的Green函数建立散射场关于速度扰动量的表达式;在数值求解中,运用奇异值分解的方法实现对第一类Fredholm积分方程的求解.通过对理论模型的计算以及实例应用,均表明瞬变电磁虚拟波场Born近似算法对地下电性界面具有良好的分辨能力. 相似文献
2.
目前,瞬变电磁数据的主要处理结果是给出地下介质的(视)电阻率信息,较难给出清晰的地质界面信息.本文基于有限差分偏移成像原理,实现了电性源瞬变电磁虚拟波场的偏移成像技术.根据瞬变电磁扩散场与虚拟波场之间的数学关系,本文使用精细积分法对第一类Fredholm积分方程求解,然后消除直达波,得到只含有地层反射波的虚拟波场.在此基础上,本文从时间域波动方程出发,对上行波的一级近似方程反向外推,得到时间-空间域二维与三维的波场递推公式.通过求解地下空间的波场,并使用反褶积方法削弱虚拟波场的波形展宽效应,从而提高对地质界面的纵向分辨率,获得较为准确的地质界面信息.本文分别采用均匀半空间模型、D型模型、HK型模型以及两层含水采空区模型对本文方法进行验证,并对实测数据进行处理,偏移成像的结果与地质资料相吻合,说明瞬变电磁有限差分偏移成像方法在识别电性界面以及地层特征方面具有可行性. 相似文献
3.
在给出真实模型和相应光滑背景模型的情况下,如何计算扰动模型(散射体)产生的散射波场是一个有实际意义的正演问题.在Gabor变换域描述散射体,且入射波场为短时宽带信号时,散射波场可以在频率域用高斯束或时间域用高斯波包描述.相对于波动方程方法,高斯束和高斯波包的计算效率更高;背景模型光滑时,高斯束和高斯波包方法的精度也接近波动方程方法.文中导出了声波假设下应用高斯束和高斯波包计算散射波的方法.测试分析了高斯波包的计算精度.给出了一般散射体的散射波模拟策略.同时针对一个理论模型完成了本文方法计算散射波的实验,实验结果表明高斯波包散射波计算方法是有效可行的. 相似文献
4.
扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果. 相似文献
5.
本文提出了一种基于双平方根算子的共炮检距道集波动方程保幅叠前深度偏移方法,将振幅误差补偿作为偏移的一部分与“运动学偏移”一起在偏移过程中实现.其基本内容包括：(1)从保幅的单平方根算子方程出发,推导出由双平方根算子定义的保幅单程波方程;(2)根据地震波摄动理论把速度场分裂为层内常速背景和变速扰动,分别在频率－波数域和频率－空间域求得波场深度延拓的偏移时移量及振幅校正系数,从而得到最终的DSR保幅波场延拓算子;(3)在高频假设条件下,把DSR保幅波场延拓公式中的积分运算进行稳相近似,得到保幅波场延拓的相移公式.理论分析和模型数值试验表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析（如AVO/AVA）提供更真实的地震信息. 相似文献
6.
通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异. 相似文献
7.
传统的剩余校正（RMO）偏移速度分析方法基于走时原理,在陡倾角和欠照明地区,因为不能得到充分的角度域信息而失效.本文将展示一种基于波场延拓理论的偏移速度分析方法,即波动方程偏移速度分析（WEMVA）.这种方法先利用成像优化方法获得剩余成像,再利用剩余成像反演剩余速度.此类方法继承了波动方程偏移方法的优点和缺点.波动方程偏移速度分析是一种线性反演方法,它要求对Born近似的展开序列作一阶截断.高阶部分的丢失必然带来巨大的截断误差,因此剩余成像必须也进行线性化,以适应大速度扰动和大延拓步长.因此,在此类算法中,剩余成像的获取和线性化是偏移速度分析的关键.在叠前偏移算子中,因为双平方根算子的数学表达式更为简洁,所以本文基于对波动方程偏移速度分析初步讨论,并通过模型验证其原理. 相似文献
8.
在背景模型基础上,求解模型扰动后的地震波散射场,这是目前地震反演中的一个关键步骤.本文将计算数学中求解非线性积分方程的Adomian分解方法,应用到求解标量波散射场的Lippmann-Schwinger积分方程和Ricatti积分方程中,分别得到了散射场的Born序列解和Rytov序列解.通过一维和二维数值算例说明:在满足一定的条件下,散射场的这两种序列解稳定收敛,与传统的Born和Rytov近似解相比,引入散射序列中的高阶项可以更精确地描述地震波散射场. 相似文献
9.
10.
通过波场反变换从瞬变电磁信号中提取虚拟波场是定位电性界面和刻画地下导电目标形态的有效手段.波场反变换过程中,降低虚拟波场速度有利于改善虚拟波场对电性界面的分辨能力.然而,波场变换核函数随时间、虚拟时间以及虚拟波场速度变化的动态范围极大,且随虚拟波场降速而进一步扩大,造成波场反变换问题高度不适定,制约了虚拟波场提取的精度和效果.对此,在虚拟波场降速的条件下,本文以核函数在无穷区间上积分的解析解为准绳,通过最优化与数据拟合得到了经验公式,确定了波场反变换所需的最优积分区间;在适应采集时间与虚拟波场速度的同时,最大限度缩短了积分区间长度,降低了核函数变化的动态范围,压制了波场反变换方程的病态程度.此外,通过引入精细积分方法求解波场反变换方程,提高了波场反变换方程的求解精度.为验证方法的有效性与可靠性,对正演模拟的电性源瞬变电磁垂直磁场分量及其时间导数分量进行了波场反变换;结果回代后的拟合数据与正演数据的相对误差小于5%.在所提取的虚拟波场记录中观察到了合理的运动学特征,如实反映了低阻目标层深度、厚度和二次场观测时间范围的变化.在计算条件相同的前提下,通过降低虚拟波场速度,改善了波场反变换的精... 相似文献
11.
有限频率层析成像考虑了非均匀介质中波的散射、衍射、波前愈合等物理性质,使得其对速度异常体的分辨能力远大于射线层析成像.推导和计算有限频率敏感核是进行有限频率层析成像的关键,当前推导有限频率敏感核多借助一阶Born近似,但这只适用于弱散射介质的情况.本文基于二阶Born近似并利用傅里叶变换推导了三维均匀介质情况下有限频率敏感核的解析表达式,并将其推广到非均匀介质中得到了三维非均匀介质中有限频率敏感核.研究表明:当介质中速度扰动小于2%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核差别很小,可近似认为相同;当介质中速度扰动大于5%时,基于二阶Born近似的有限频率敏感核与基于一阶Born近似的有限频率敏感核有较大不同,表明此时已不能忽略二次散射. 相似文献
12.
最小二乘偏移是一种基于反射地震数据与地下反射率间线性关系而建立起来的地震数据线性反演方法,相比常规偏移成像具有更好的保幅性能.本文提出了一种基于照明补偿的单程波最小二乘偏移方法,首先利用单程波方程的稳定Born近似广义屏波场传播算子构建反射地震数据与地下反射率间的线性算子,然后再应用线性最优化方法求解最小二乘偏移所对应的线性反问题.在迭代求解最优化问题的过程中,以地震波场的地下照明强度作为迭代反演的预条件算子加快迭代的收敛速度.单程波传播过程中考虑了速度分界面产生的透射效应,并用单极震源代替常规偏移中的偶极震源.把本文提出的方法应用于层状理论模型和Marmosi模型地震数据的数值试验中均取得了理想的结果. 相似文献
13.
城市地下空间的高精度、多分辨探测是城市地下合理开发的前提.由于城市探测环境的特殊性,需要探测方法具备抗干扰能力强、分辨能力好、分异多尺度地下目标等特性,现有的物探方法难以兼顾城市地下空间的探测需要.本文利用三维矢量有限元方法为正演手段,在瞬变电磁高性能辐射源进行微分脉冲扫描的基础上,对探测数据进行多时窗的扫时波场变换,将微分脉冲扫描后的多分辨响应信息,进行多分辨信息提取.同时通过地震探测中的多次覆盖处理,提高探测场对地下目标的分辨能力,最后对多次叠加后的虚拟波场进行拟地震偏移成像,最终实现城市地下空间的高精度探测. 相似文献
14.
声波方程的逆散射反演乃是求解双曲型偏微分方程系数项反问题的一种解析方法,一般利用Born近似把这一非线性反问题线性化,并给出了恒参考波速介质中反问题解的解析表达式.由于Born近似假定波速扰动为一级无穷小,因此,在大多数情况下,恒参考波速介质模型的反问题的解无法得以应用.本文研究介质参考波速沿某个方向线性变化时的声散射理论,导出了声波方程逆散射问题解的解析表达式,从而既可使Born近似的假定在大多数情况下能得以满足,又可利用快速Fourier变换快速实现介质波速扰动的反演成象. 相似文献
15.
Seismic inversion with generalized Radon transform based on local second-order approximation of scattered field in acoustic media 
下载免费PDF全文
下载免费PDF全文 Sound velocity inversion problem based on scattering theory is formulated in terms of a nonlinear integral equation associated with scattered field. Because of its nonlinearity, in practice, linearization algorisms (Born/single scattering approximation) are widely used to obtain an approximate inversion solution. However, the linearized strategy is not congruent with seismic wave propagation mechanics in strong perturbation (heterogeneous) medium. In order to partially dispense with the weak perturbation assumption of the Born approximation, we present a new approach from the following two steps: firstly, to handle the forward scattering by taking into account the second-order Born approximation, which is related to generalized Radon transform (GRT) about quadratic scattering potential; then to derive a nonlinear quadratic inversion formula by resorting to inverse GRT. In our formulation, there is a significant quadratic term regarding scattering potential, and it can provide an amplitude correction for inversion results beyond standard linear inversion. The numerical experiments demonstrate that the linear single scattering inversion is only good in amplitude for relative velocity perturbation ( \( \delta_{c}/c_{0} \) ) of background media up to 10 %, and its inversion errors are unacceptable for the perturbation beyond 10 %. In contrast, the quadratic inversion can give more accurate amplitude-preserved recovery for the perturbation up to 40 %. Our inversion scheme is able to manage double scattering effects by estimating a transmission factor from an integral over a small area, and therefore, only a small portion of computational time is added to the original linear migration/inversion process. 相似文献
16.
在矿井瞬变电磁超前探测数据处理和解释中,为了更好地识别巷道前方含水体电性分界面,利用扩散场与波动场间的函数关系,将矿井瞬变电磁法超前探测中瞬变电磁波场转换为拟地震波场,实现波场转换.同时,为了解决波场转换后信号幅值偏弱,以增加信噪比利于解释,利用雷达的合成孔径方法对波场转换后的拟地震信号进行相关叠加处理.将矿井瞬变电磁法巷道超前探测的同一测点不同探测方向的数据作为孔径数据进行合成孔径成像,以突出弱异常而提高信噪比,结合应用实例验证了方法的效果.研究表明,矿井瞬变电磁法合成孔径成像方法能够提高电性界面信息的信噪比和界面几何分辨率,对矿井瞬变电磁法精细超前探测含水构造具有重要的理论意义和应用价值. 相似文献
17.
本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率. 相似文献