提高任意广角波动方程计算效率及压制倏逝波干扰方法

任意广角波动方程（Arbitrarily Wide-angle Wave Equation,AWWE）是一种高精度的时空域单程波动方程,与其他时空域单程波方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文改进了现有的AWWE有限差分计算方案,避免了原计算方案中对矩阵反复求逆的过程,同时减少了矩阵相乘的次数,显著地提高了AWWE数值计算的效率.数值实验证实,采用改进之后的计算方案,其计算结果与原方案完全一致.此外,本文采用了一种方法,在有倏逝波干扰的情况下,能在频率-波数（f-k）域区分出倏逝波与非倏逝波,为设计f-k滤波器提供合理的视速度门槛信息.偏移实验证实,在深度延拓时仅仅进行数次f-k滤波就能很好地抑制倏逝波干扰,提高偏移成像质量.     

标量波动方程全倾角有限差分偏移

传统的有限差分法偏移是建立在单程波波动方程的不同阶数的近似式基础上的。因此,各阶近似式的偏移方法都存在一个倾角极限问题。克服倾角极限和提高极限倾角是八十年代以来有限差分偏移技术的研究目标。在这方面取得了显著的进展。本文从标量波动方程出发,通过函数替代,将它化为一个等价的方程组,用有限差分法解此方程组进行波场外推,实现全倾角偏移成象的目的。与国外已有的使用波动方程全式进行偏移的方法相比较,除方法原理本身不同外,该方法具有稳定性好,计算工作量较少的优点。     

标量波动方程全倾角有限差分偏移

传统的有限差分法偏移是建立在单程波波动方程的不同阶数的近似式基础上的。因此,各阶近似式的偏移方法都存在一个倾角极限问题。克服倾角极限和提高极限倾角是八十年代以来有限差分偏移技术的研究目标。在这方面取得了显著的进展。本文从标量波动方程出发,通过函数替代,将它化为一个等价的方程组,用有限差分法解此方程组进行波场外推,实现全倾角偏移成象的目的。与国外已有的使用波动方程全式进行偏移的方法相比较,除方法原理本身不同外,该方法具有稳定性好,计算工作量较少的优点。     

波动方程Maslov射线解的有限差分计算

借助于守恒型程函方程的有限差分解法，构造出一种地震波场Ｍａｓｌｏｖ渐近解的数值计算方法，并针对不同坐标系进行了讨论。文中通过不同模型分别计算了透射波和反射波的射线路径、走时及相应的合成地震记录，结果表明本算法具有快速、准确的特点。本文算法可用于叠前Ｋｉｒｃｈｈｏｆｆ深度偏移及非均匀介质逆散射问题中Ｇｒｅｅｎ函数的计算，并就算法本身的局限性作了讨论。     

高精度有限差分地震波正演方法

三维地震模拟不仅可以更精确地研究地震波的传播规律,而且是三维地震资料处理和解释的工具。本文采用精细积分法用于求解波动方程.对波动方程在空间采用差分格式,时间域采用积分法求解析的方法。文中详细论述了精细积分法的数值方法,并给出了计算公式。理论分析和数值算例的结果表明了用这种混合方法得到的解与精确解十分吻合,比有限差分法具有更高的精度。文中给出的地震波正演模拟算例说明了该方法适用于复杂地表和复杂构造地质体。     

三维波动方程有限差分正演方法

三维地震资料的处理和解释都需要有效的三维正演模型予以验证．本文提出一种在ｘ－ｔ域实现快速、高精度的有限差分正演方法,采用了独特的“平行四边形网格”,并用Ｐ－Ｒ交替差分格式使三维波动方程可以局部地分裂成二维求解方程,从而有效地减少运算量．     

三维波动方程有限差分正演方法

三维地震资料的处理和解释都需要有效的三维正演模型予以验证．本文提出一种在ｘ－ｔ域实现快速、高精度的有限差分正演方法，采用了独特的“平行四边形网格”，并用Ｐ－Ｒ交替差分格式使三维波动方程可以局部地分裂成二维求解方程，从而有效地减少运算量．     

任意广角波动方程深度偏移参数优化分析

强横向变速和陡倾角介质的成像问题是地震偏移成像的难点。本文在频率-空间域利用有限差分高阶分裂法求解任意广角波动方程实现了地震偏移成像。试验表明,参考速度的选取是影响成像精度的关键因素。脉冲响应测试表明通过优化参数任意广角波动方程能够以较低阶的方程获得较大的偏移角度。Marmousi模型数据叠前深度偏移试验表明,该方法能够适应强横向变速和陡倾角地质结构的成像,参数优化后能够取得更好的成像效果。     

基于任意广角波动方程的频率-空间域深度偏移方法研究

传统的单程波动方程偏移算法对大倾角成像困难,本文基于时空域任意广角单程声波方程,通过对时间变量进行傅立叶变换得到其频率-空间域形式,利用有限差分进行离散化,设计并实现了频率-空间域有限差分叠后偏移成像算法。模型试算表明,该方法能够通过参数优化使得低阶数算子适应较大倾角地层的偏移成像。     

基于精细积分法的三维弹性波数值模拟(英文)

波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。     

Arbitrary Difference Precise Integration Method for Solving the Seismic Wave Equation

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.     

基于双差定位法的鄂尔多斯块体西北缘地震精定位研究

针对鄂尔多斯块体西北缘地震活动的复杂性及目前台网定位方法存在偏差的问题,采用双差定位法对2009—2019年发生在鄂尔多斯块体西北缘的地震事件进行重新定位.经研究表明,双差定位之后地震分布更集中,鄂尔多斯块体西缘的地震沿银川吉兰泰断陷带分布,北缘的地震沿河套断陷带分布,地震定位精度明显提高,这与块体周缘复杂的地质构造背...     

频率波数域内层状地基动力柔度的精细求解

精细积分法既可得到在计算机精度意义下的精确解,又能够保持哈密顿体系的辛结构。其是求解一阶线性常微分方程组的精确数值方法,既可以用于时间域的初值问题,又可以应用于空间域的两点边值问题。运用精细积分法求解微层区段矩阵,并对微层区段矩阵合并得到整个层状地基的区段矩阵,最终得到层状地基的动力柔度值。运用数值算例验证了本文方法的计算精度。