共查询到19条相似文献,搜索用时 208 毫秒
1.
2.
《应用地球物理》2017,(2)
地震波场正演模拟是地震资料处理、解释中最为重要的技术之一。地震波场正演模拟在大时间步长、长时程的波场延拓中,存在计算不稳定的问题。本文基于声波方程的Hamilton表述,在波动方程求解中用辛差分格式进行时间网格离散,用傅里叶有限差分进行空间网格离散,提出一种新的保结构地震波场正演模拟方法一辛格式傅里叶有限差分法,在保证计算精度的同时提高计算的稳定性。利用声学近似处理空间-波数混合域的积分算子,将该方法推广至各向异性介质。给出各向同性和各向异性条件下的地震正演模拟的计算流程,并将本文方法用于BP盐丘、BP TTI等模型的波场正演模拟。数值算例表明本文开发的方法适用于速度变化剧烈的复杂介质地震波场正演模拟,计算精度高,数值频散小,在各向异性介质正演中能够有效避免qSV波残余,在大时间步长的迭代计算中稳定性好。本文为在辛算法的框架下实现高精度地震正演模拟提供了一种新的选择。 相似文献
3.
间断有限元(Discontinuous Galerkin:DG)方法具有低数值频散、网格剖分灵活、能模拟地震波在复杂介质中传播等优点.因此,本文将一种新的DG方法推广到双相和黏弹性等复杂介质的地震波场模拟,发展了求解Biot弹性波方程和D'Alembert介质波动方程的DG方法.首先通过引入辅助变量将Biot双相介质弹性波方程和D'Alembert介质波动方程转化为关于时间-空间的一阶偏微分方程组,然后对该方程组进行DG空间离散,得到半离散化的常微分方程组.最后,对此常微分方程组,应用加权的Runge-Kutta格式进行时间推进计算.数值结果表明,DG方法可以有效地求解Biot双相介质弹性波方程和D'Alembert介质波动方程,并能很好地压制因离散求解波动方程而产生的数值频散,获得清晰的各种地震波震相. 相似文献
4.
地震波传播过程本质上是能量在传播过程中逐步损耗直至殆尽的过程,而在实际应用中,常在无能量损耗假设下,用弹性波动方程或标量波动方程描述它.在哈密顿(Hamilton)体系表述下,地震波传播过程即为一个无限维的哈密顿系统随时间的演化过程.若不计能量损耗,波场演化过程实质上为一个单参数连续的辛变换,因而对应的数值算法应为辛几何算法.本文首先从地震波标量方程出发,给出哈密顿体系下地震波传播的表述,即任意两个时刻的波场是通过辛变换联系起来的.随后,把波场在时间和相空间离散化后,给出了用于波场计算的一些辛格式,如显式辛格式、隐式辛格式和蛙跳辛格式.并进一步讨论了有限差分格式和辛格式的异同.然后,应用显式辛格式和同阶的有限差分方法给出了同一理论速度模型下的波场和Marmousi速度模型下的单炮记录.数值结果表明,辛算法是一类可行的波场模拟的数值算法.在时间步长较小时,有限差分方法是辛算法的一个很好近似.文中的理论和方法,为地震波传播理论及实际应用研究提供了新的途径. 相似文献
5.
三角网格有限元法具有网格剖分的灵活性,能有效模拟地震波在复杂介质中的传播.但传统有限元法用于地震波场模拟时计算效率较低,消耗较大计算资源.本文采用改进的核矩阵存储(IKMS)策略以提高有限元法的计算效率,该方法不用组合总体刚度矩阵,且相比于常规有限元法节省成倍的内存.对于时间离散,将有限元离散后的地震波运动方程变换至Hamilton体系,在显式二阶辛Runge-Kutta-Nystr9m(RKN)格式的基础之上加入额外空间离散算子构造修正辛差分格式,通过Taylor展开式得到具有四阶时间精度时间格式,且辛系数全为正数.本文从理论上分析了时空改进方法相比传统辛-有限元方法在频散压制、稳定性提升等方面的优势.数值算例进一步证实本方法具有内存消耗少、稳定性强和数值频散弱等优点. 相似文献
6.
地震波场数值模拟是理解地震波在地下介质中的传播特点,帮助解释观测数据的有效手段,而提高计算精度和运算效率是所有波场数值模拟方法研究所追求的目标.有限差分技术是求解波动方程计算效率最高、应用最为广泛的方法之一.但传统的有限差分技术计算过程中的数值频散问题影响了该技术的计算精度与计算效率.本文通过交错网格高阶有限差分技术与通量校正传输方法(Flux|corrected transport method,FCT)相结合, 对横向各向同性介质(Transverse isotropic medium,TI)一阶速度|应力弹性波动方程组进行了数值求解研究.波场快照数值模拟结果表明,本文研究的数值模拟方法与波动方程二阶有限差分方法、交错网格四阶有限差分方法相比,在压制网格数值频散方面有明显的优势,计算精度提高,而且可以利用较大的空间步长,提高计算效率. 相似文献
7.
早期的褶积微分算子法都是基于正反傅立叶变换而实现的,其精度比四阶有限差分稍高。本文将计算数学中的Forsyte广义正交多项式微分算子与褶积算子相结合,构建了一个新的快速、高精度褶积微分算子,其计算结果非常接近实验函数微分的精确值,精度与l6阶有限差分相当。粘弹性波动方程更真实地描述了实际地下介质中弹性波的传播规律及其波场特征。本文以二维粘弹性波动方程为例,推导了粘弹性介质波动方程的离散格式,用迭积微分算子法实现了粘弹性介质的地震波场正演模拟,并对其波传播特征进行了分析。计算结果表明该算法能正确模拟粘弹性介质中的地震波,正确地反映粘弹性介质中波场的传播规律。 相似文献
8.
波动方程数值模拟方法是研究地震波场传播的一种重要手段,本文采用交错网格高阶有限差分方法分别对双程声波方程和双程弹性波方程进行了波场数值模拟,并且根据定位原理采用傅立叶有限差分算子进行了单程波方程数值模拟,在分析定位原理的基础上,对其计算过程稍作修改,将延拓到地面的波场直接由每个检波点接收,无需横向叠加过程,得到了单程声波方程共炮记录.基于不同波动方程的数值模拟结果表明,双程波方程结果包含直达波、多次波等干扰波,信噪比低;单程波数值模拟结果只包含了介质分界面的一次反射波,信噪比高,但对于大角度入射波误差较大,并且对于同一个地质模型而言,双程弹性波方程计算速度最慢,双程声波方程次之,单程声波方程计算速度最快.因此对于复杂地质模型,三种模拟方法可以取长补短,综合应用. 相似文献
9.
10.
波动方程数值模拟是深入研究地震波传播规律的有效方法.有限差分法因其方法简单、精度高而得到了广泛的应用.但其缺点是不能准确模拟具有复杂几何形态的物性界面.因而当遇到起伏地表或复杂构造时,求解精度低.为了准确模拟起伏地形、复杂构造和复杂介质条件下的地震波场,本文采用有限元法模拟二维声波方程.用三角形单元模拟地形和速度界面;把单元内的场和波速均看作单元上的线性函数,以适应复杂介质压制边角散射;采用吸收边界条件去除来自截断边界上的反射;采用集中质量矩阵和集中阻尼矩阵使得显式时间递推无需对矩阵求逆,提高了计算效率.对模型的计算表明该方法正确有效. 相似文献
11.
针对目前广泛用于解决工程地质问题的浅层地震勘探特点,通过对浅层小尺度介质进行地震波数值模拟,揭示地震波在不同浅层介质模型中的波场特征和传播规律.为浅层地震勘探野外数据采集前合理选择观测系统提供指导性的方案,并能对地震勘探成果解释的合理性给予客观评价.在分析各种数值模拟方法特点的基础上,选择具有精度高、易模拟复杂结构的有... 相似文献
12.
13.
基于精细积分法的三维弹性波数值模拟(英文) 总被引:1,自引:0,他引:1
波动方程有限差分法是地震数值模拟中的一种重要的方法,对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想,产生了不稳定性。精细积分法在有限差分法的基础上,在时间域采用解析解的表达形式,在空间域保留任意差分格式,发展成为半解析的数值方法。本文结合并发展了以往学者的成果,推导了任意精细积分法的三维弹性波正演模拟计算公式,并对其稳定性进行了数值分析。在计算实例中,实现了精细积分法二维和三维弹性波模型的地震正演模拟,对计算结果的分析表明,精细积分法反射信号走时准确,稳定性好,弹性波场相较于声波波场,弹性波波场成分更为丰富,包含了更多波型成分(PP-和PS-反射波、透射波和绕射波),这对实际地震资料的解释和储层分析有重要的意义。实践证明,该方法可直接应用到弹性波的地质模型的数值模拟中。 相似文献
14.
简单论述了横波的基本特点,指出在松散层勘探中,横波勘探的优势主要在于探测精度高且不受饱水度的影响.因此,可以选择地震横波勘探进行第四纪松散层层序划分及厚度探测.并以广州市城市地质调查项目为例,介绍了地震横波勘探在浅层地层勘测中的数据采集方法,分析了地质解释依据,并结合波速测井方法对地层层面进行了划分.分析结果表明了多数岩土层与反射界面有较好的对应关系,用地震横波勘探结合波速测井进行第四纪松散层层序划分及厚度探测是可行的. 相似文献
15.
本文用一个纵波信息,对一维声波方程的速度和源函数进行联合反演,并考虑到声波方程的反问题是一个不适应问题,对源函数和波速分别和正则化法分步迭代求解,减少反问题的计算工作量,改善该问题的计算稳定性,为计算实际工程和岩性学问题供了一种方法。文中给出只用一个反问题补充条件同时进行多参数反演的公式,并对相应的数值算例进行分析和比较。 相似文献
16.
将波动方程变换至Hamilton体系,构造了一种新的保结构算法,即最优化辛格式广义褶积微分算子(OSGCD). 在时间离散上,首先引入了Lie算子设计二级二阶辛格式,基于最小误差原理得到了优化的辛格式. 在空间离散上,引入广义离散奇异核褶积微分算子计算空间微分,提出了一种有效方法优化GCD并得到了稳定的算子系数. 针对本文发展的新方法,给出了OSGCD稳定性条件. 在数值实验中,将OSGCD与多种方法比较,从精度和计算效率两方面分析了OSGCD的计算优势,计算结果也表明OSGCD长时程以及非均匀介质中地震波模拟亦具有较强能力. 相似文献
17.
18.
地震波走时和射线的有限差分计算 总被引:5,自引:0,他引:5
以往都是采用射线追踪的方法计算地震波的走时和射线,但是当速度模型复杂时这种方法存在一些问题。本文提出另一种计算地震波走时和射线的方法。该方法从程函方程出发,利用互换原理和Fermat原理计算出各种波的到时和射线。解决了射线追踪方法存在的问题。为地震波走时和射线的计算以及地震波走时反演开辟了一条新途径。 相似文献
19.
初步探讨小波多分辨分析理论在地震波动方程正演模拟中的应用,给出一种基于地震波场局部变化性质而自适应调空间网格点的波动方程数值算法,目前求解波动方程所用的有限差分法、有限元法以及伪谱法都不能根据波场的局部变化性质而动态选择空间网格点的大小、小波基函数在空间域和频域中都具有局部性特征,它的性质优于有限差分法,有限元法和伪谱法中所用的基函数,通过阀值运算,地震波场失发辨表示变得非常稀疏,同时地下介质中的主要信息又不会受到损害,本文将声波场的多分辨表示变得非常稀疏,同时地下介质中的主要信息又不会受到损害,本文将声波方程的矩阵表示形式小波多分辨分析的框架下进行了展开,通过对算子长阵的地震波场矩阵进行多分辨分解和压缩,得到了小波域中地震波场正演模拟算法。 相似文献