共查询到19条相似文献，搜索用时 156 毫秒
1.
为更好地研究复杂构造和速度分布条件下地震资料的精确偏移.在波场模拟方面,本文采用了高阶交错网格有限差分[1]推导了二维各向同性介质声波方程的数值模拟公式.并给出了声波方程完全匹配层吸收边界[2-5](PML).在偏移方面,本文介绍了一种能适应介质速度横向变化的非稳态相移算子(NSPS)[6]及其叠后深度偏移方法.并对简单的四层模型和Marmousi模型进行试算.偏移结果表明,该方法可以有效地处理复杂条件下的各向同性介质声波偏移问题. 相似文献
2.
从含Thomsen各向异性参数的qP波相速度表示式出发,建立并求解三维VTI介质中的频散方程,得到三维VTI介质中的相移算子,进而将以相移算子为基础的最优分裂Fourier方法推广到三维VTI介质,发展了一个三维VTI介质的深度偏移方法.文中使用的各向异性介质的速度模型与现行的各向异性构造的速度估计方法一致,将各向同性、弱各向异性及强各向异性统一在一个模型中.文中提出的偏移算法对相移法引入了高阶校正项来补偿介质横向变化的影响,使该方法可应用于横向非均匀VTI介质的陡角度成像,文中给出的偏移脉冲响应很好地证明了这一点. 相似文献
3.
横向各向同性介质是地球内部广泛分布的一种各向异性介质.针对这种介质，我们对各向同性介质的最小走时树走时模拟方法进行了推广，推广后的方法可适用于非均匀、对称轴任意倾斜的横向各向同性介质模型.为保证计算效率，最小走时树的构建采用了一种子波传播区域随地震波传播动态变化的改进算法.对于弱各向异性介质，我们使用了一种新的地震波群速度近似表示方法，该方法基于用射线角近似表示相角的思想，对3种地震波（qP, qSV和qSH）均有较好的精度.应用本文地震波走时模拟方法对均匀介质、横向非均匀介质模型进行了计算，并将后者结果与弹性波方程有限元方法的模拟结果进行了对比，结果表明两者符合得很好.本文方法可用于横向各向同性介质的深度偏移及地震层析成像的深入研究. 相似文献
4.
根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果． 相似文献
5.
为了处理横向强变速介质中的深度成像问题,本文提出一种基于共炮道集的优化系数的傍轴近似方程叠前深度偏移算子,并在基于反射系数估算的成像条件下,可实现叠前深度偏移成像.该算子具有方程阶数低且能对陡倾角成像的特征,并采用有限差分法波场延拓,能适应速度场的任意变化.当在频率-空间域进行计算时,相对于纯粹的时间-空间域有限差分算法有计算效率高、成像方便的优点.脉冲响应测试和对Marmousi模型进行的叠前深度偏移结果表明,该偏移方法在强横向变速情况下具有非常好的成像效果. 相似文献
6.
非稳态相移法叠前深度偏移 总被引：5，自引：3，他引：2
介绍一种能够适应介质速度横向变化的非稳态相移算子及其叠前深度偏移方法．为了克服常规相移偏移算法中要求速度横向不变的缺点,出现了基于非稳态滤波器理论的非稳态相移算子,即PSPI算子、NSPS算子和SNPS算子,其中SNPS算子是将前二者结合起来的一种对称的非稳态相移算子,它比前二者具有更高的精度和稳定性．为了提高运算速度,基于非稳态相移算子的叠前深度偏移算法采取了分片均匀近似的策略,Marmousi模型的叠前深度偏移结果证明了该算法的可行性和有效性。 相似文献
7.
利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明:基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性. 相似文献
8.
波动方程Maslov射线解的有限差分计算 总被引：2，自引：0，他引：2
借助于守恒型程函方程的有限差分解法，构造出一种地震波场Ｍａｓｌｏｖ渐近解的数值计算方法，并针对不同坐标系进行了讨论。文中通过不同模型分别计算了透射波和反射波的射线路径、走时及相应的合成地震记录，结果表明本算法具有快速、准确的特点。本文算法可用于叠前Ｋｉｒｃｈｈｏｆｆ深度偏移及非均匀介质逆散射问题中Ｇｒｅｅｎ函数的计算，并就算法本身的局限性作了讨论。 相似文献
9.
合成三维横向非均匀介质远震体波接收函数的Maslov方法 总被引：5，自引：3，他引：2
三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要．本文发展了基于Ｍａｓｌｏｖ理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪,并采用递归算法组织射线的追踪过程,从而在保证追踪精度的同时提高了射线追踪的效率,并可完全避免传统射线追踪过程中的射线编码问题．正演计算结果表明,此法是一种高效、实用的三维横向非均匀介质接收函数计算方法． 相似文献
10.
11.
阐述了矩阵多分辨分析理论中的标准形式与非标准形式,并以Hilbert算子为例,说明了算子多分辨表示的压缩效果,为小波域偏移算法奠定了理论基础.从三维雷达波动方程出发,利用爆炸反射原理和浮动坐标变换,推导出三维探地雷达波动方程差分格式,并通过方程分裂算法及小波多分辨算法,在小波域求解波场外推矩阵,进而得到探地雷达小波域三维波动方程偏移算法,在此基础上,开发了探地雷达小波域偏移处理程序,并把该程序应用于三个球体空洞的3-D正演结果及实际的雷达数据中,通过对比偏移处理前后的雷达资料,得知该三维偏移算法能使3-D正演剖面中的反射波归位、绕射波收敛,极大地提高了雷达剖面的分辨率,有利于探地雷达资料的地质解释. 相似文献
12.
探讨了井巷隧道工程反射波法超前地质预报成像问题.从一阶速度－应力弹性波方程出发,推导了二维各向同性介质情况下弹性波逆时传播的高阶差分格式,实现了弹性波在数值空间中的逆时延拓.构建反射波法隧道超前探测中的断层、软弱夹层等介质模型,以反射波探测的正演记录作为初始条件,并从程函方程出发,采用逆时差分格式求取介质模型网格空间中各点的直达波旅行时作为弹性波逆时偏移的成像条件,实现多波多分量资料的逆时偏移.偏移结果表明,逆时偏移能够使隧道壁接收到的波场准确归位,提高隧道反射波超前探测的资料处理的精度. 相似文献
13.
利用地震波正向传播方程对属于波形线性反演问题近似求解方法的地震数据偏移成像进行重新推导,得到了适合散射地震数据的散射偏移成像方法和适合反射地震数据的反射偏移成像方法.以地震波传播的散射理论为出发点,首先根据描述一次散射波正向传播的线性方程研究建立散射地震数据的偏移成像方法理论;利用高频近似对产生散射波场的地下速度扰动函数的空间变化进行近似,推导出地下反射率函数,再由散射波传播方程推导出基于反射率函数的反射波传播方程,然后根据描述一次反射波正向传播的线性方程研究建立反射地震数据的偏移成像方法理论.本文指出和修正了Claerbout偏移成像方法中的不足,提出的地震数据偏移成像方法是对当前偏移成像方法理论的完善,使反射地震数据偏移成像具有了更坚实的数学物理理论基础,得到的偏移成像结果相位正确、位置准确、分辨率提高. 相似文献
14.
最小二乘偏移是一种基于反射地震数据与地下反射率间线性关系而建立起来的地震数据线性反演方法,相比常规偏移成像具有更好的保幅性能.本文提出了一种基于照明补偿的单程波最小二乘偏移方法,首先利用单程波方程的稳定Born近似广义屏波场传播算子构建反射地震数据与地下反射率间的线性算子,然后再应用线性最优化方法求解最小二乘偏移所对应的线性反问题.在迭代求解最优化问题的过程中,以地震波场的地下照明强度作为迭代反演的预条件算子加快迭代的收敛速度.单程波传播过程中考虑了速度分界面产生的透射效应,并用单极震源代替常规偏移中的偶极震源.把本文提出的方法应用于层状理论模型和Marmosi模型地震数据的数值试验中均取得了理想的结果. 相似文献
15.
本文基于拟微分算子理论和李代数积分法,根据程函方程和波场坐标变换,提出一种新的适于横向变速介质Kirchhoff叠前深度偏移的地震波走时算法.该算法与Kirchhoff叠前时间偏移所用李代数时间积分表达相比,差异在于增加了波数一次项,且二次项的系数在求积时亦需进行修正.针对单平方根算子象征、李代数积分、指数映射和走时多项式的求解而言,皆需对以往Kirchhoff叠前时间偏移中所用算法进行深化调整.文中数值算例对比了本文李代数积分表达与时间积分的区别,本算法计算结果与线性横向变速介质中的理论值相当吻合.通过走时多项式中各项对结果的影响分析,可知非对称项使计算精度得到了进一步提高.数值试验表明,本算法对横向变速介质中走时求取是可行的,且不需要存储海量走时表,有利于提高Kirchhof叠前深度偏移的精度和效率. 相似文献
16.
基于任意广角波动方程的频率-空间域深度偏移方法研究 总被引：1，自引：1，他引：0
孙歧峰 《CT理论与应用研究》2010,19(2):17-24
传统的单程波动方程偏移算法对大倾角成像困难,本文基于时空域任意广角单程声波方程,通过对时间变量进行傅立叶变换得到其频率-空间域形式,利用有限差分进行离散化,设计并实现了频率-空间域有限差分叠后偏移成像算法。模型试算表明,该方法能够通过参数优化使得低阶数算子适应较大倾角地层的偏移成像。 相似文献
17.
有限差分法广泛应用于地震波场的数值延拓,确定合适的有限差分算子以减小数值频散是有限差分法的一个重要研究内容。近年来为了进一步抑制数值频散和增加时间步长,新的有限差分模板得到了应用,对于此,前人使用泰勒展开方法和最小二乘方法确定有限差分算子系数。本文在以前工作的基础上,使用改进的线性方法确定新模板的有限差分系数,并与传统模板线性方法进行对比;通过频散分析和正演模拟验证出新模板线性方法能够更好地保持频散关系,在相同的精度下效率提高了一倍,从而说明了改进的线性方法的有效性。 相似文献
18.
19.
叠前深度偏移理论及方法一直是地震数据成像中研究的热点问题.业界对单程波叠前深度偏移方法和逆时深度偏移开展了深入的研究,但对双程波方程波场深度延拓理论及成像方法的研究还鲜有报道.本文以地表记录的波场值为基础,利用单程波传播算子估计波场对深度的偏导数,为在深度域求解双程波方程提供充分的边界条件,并提出利用矩阵分解理论实现双程波方程的波场深度外推.通过对强速度变化介质中传播波场的计算,与传统的单程波偏移方法相比,本文提出的偏移方法计算的波场与常规有限差分技术计算的波场相一致,证明了本方法计算的准确性.通过对SEAM模型的成像,在相同的成像参数下,与传统的单程波偏移算法和逆时深度偏移算法方法相比,本文提出的偏移方法能够提供更少的虚假成像和更清晰的成像结果.本文所提偏移算法具有深度偏移和双程波偏移的双重特色,推动和发展了双程波叠前深度偏移的理论和实践. 相似文献