共查询到13条相似文献,搜索用时 46 毫秒
1.
任意广角波动方程(Arbitrarily Wide-angle Wave Equation,AWWE)是一种高精度的时空域单程波动方程,与其他时空域单程波方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文改进了现有的AWWE有限差分计算方案,避免了原计算方案中对矩阵反复求逆的过程,同时减少了矩阵相乘的次数,显著地提高了AWWE数值计算的效率.数值实验证实,采用改进之后的计算方案,其计算结果与原方案完全一致.此外,本文采用了一种方法,在有倏逝波干扰的情况下,能在频率-波数(f-k)域区分出倏逝波与非倏逝波,为设计f-k滤波器提供合理的视速度门槛信息.偏移实验证实,在深度延拓时仅仅进行数次f-k滤波就能很好地抑制倏逝波干扰,提高偏移成像质量. 相似文献
2.
传统的有限差分法偏移是建立在单程波波动方程的不同阶数的近似式基础上的。因此,各阶近似式的偏移方法都存在一个倾角极限问题。克服倾角极限和提高极限倾角是八十年代以来有限差分偏移技术的研究目标。在这方面取得了显著的进展。本文从标量波动方程出发,通过函数替代,将它化为一个等价的方程组,用有限差分法解此方程组进行波场外推,实现全倾角偏移成象的目的。与国外已有的使用波动方程全式进行偏移的方法相比较,除方法原理本身不同外,该方法具有稳定性好,计算工作量较少的优点。 相似文献
3.
4.
波动方程Maslov射线解的有限差分计算 总被引:2,自引:0,他引:2
借助于守恒型程函方程的有限差分解法,构造出一种地震波场Maslov渐近解的数值计算方法,并针对不同坐标系进行了讨论。文中通过不同模型分别计算了透射波和反射波的射线路径、走时及相应的合成地震记录,结果表明本算法具有快速、准确的特点。本文算法可用于叠前Kirchhoff深度偏移及非均匀介质逆散射问题中Green函数的计算,并就算法本身的局限性作了讨论。 相似文献
5.
6.
7.
8.
强横向变速和陡倾角介质的成像问题是地震偏移成像的难点。本文在频率-空间域利用有限差分高阶分裂法求解任意广角波动方程实现了地震偏移成像。试验表明,参考速度的选取是影响成像精度的关键因素。脉冲响应测试表明通过优化参数任意广角波动方程能够以较低阶的方程获得较大的偏移角度。Marmousi模型数据叠前深度偏移试验表明,该方法能够适应强横向变速和陡倾角地质结构的成像,参数优化后能够取得更好的成像效果。 相似文献
9.
基于任意广角波动方程的频率-空间域深度偏移方法研究 总被引:1,自引:1,他引:0
孙歧峰 《CT理论与应用研究》2010,19(2):17-24
传统的单程波动方程偏移算法对大倾角成像困难,本文基于时空域任意广角单程声波方程,通过对时间变量进行傅立叶变换得到其频率-空间域形式,利用有限差分进行离散化,设计并实现了频率-空间域有限差分叠后偏移成像算法。模型试算表明,该方法能够通过参数优化使得低阶数算子适应较大倾角地层的偏移成像。 相似文献
10.
基于精细积分法的三维弹性波数值模拟(英文) 总被引:1,自引:0,他引:1
波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。 相似文献
11.
Jia Xiaofeng 《中国地震研究》2004,18(1):36-41
Wave equation migration is often applied to solve seismic imaging problems. Usually, the finite difference method is used to obtain the numerical solution of the wave equation. In this paper, the arbitrary difference precise integration (ADPI) method is discussed and applied in seismic migration. The ADPI method has its own distinctive idea. When dispersing coordinates in the space domain, it employs a relatively unrestrained form instead of the one used by the conventional finite difference method. Moreover, in the time domain it adopts the sub-domain precise integration method. As a result, it not only takes the merits of high precision and narrow bandwidth, but also can process various boundary conditions and describe the feature of an inhomogeneous medium better. Numerical results show the benefit of the presented algorithm using the ADPI method. 相似文献
12.