地震波场模拟的伪谱和高阶有限差分混合方法

来源 :北京大学 | 被引量 : 0次 | 上传用户:scg5252
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文首先系统回顾了波动方程地震数值模拟方法的发展过程和研究现状,主要介绍了地震波场数值模拟的两种常用方法:交错网格的伪谱法和有限差分法,并详细比较了二者的优势和劣势。伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值方法,但是由于它的微分算子的全局性,使得该方法不适用于分散内存的并行计算。 本文将有限差分算子的局部性和伪谱法算子的高效、高精度相结合,发展基于两种方法的PSM/FDM混合方法,该方法在一个空间坐标方向上利用交错网格高阶有限差分算子,在另外的空间坐标方向上利用交错网格伪谱法算子。我们给出了混合方法的数值实现过程,并对均匀弹性半空间进行计算,与伪谱法的计算结果做比较,分析了计算精度和效率,显示了混合方法的优越性。最后,应用交错网格混合方法模拟了二维椭圆形低速盆地模型的P-SV波传播,计算结果充分显示了低速盆地的场地效应。 本文的研究结果表明,基于伪谱法和有限差分方法的数值模拟混合方法,结合了有限差分法和伪谱法的优点,弥补了二者的不足,既能有效的较好地处理介质不连续面的计算,又能在保证和伪谱法相当精度的前提下,提高计算效率,节省计算时间。由于在垂直方向上采用局部有限差分算子,混合方法很适宜于在分布内存的PC集群上实现并行计算,在大规模三维非均匀介质模型小波长地震波场模拟方面有应用价值。
其他文献
磁层亚暴活动反映了磁层和太阳风相互耦合作用的过程,对亚暴活动的研究是日地物理研究中的重要课题。磁层亚暴的预报是人类空间环境预报的一个重要组成部分。通常用地磁活动
哈密位于新疆东部,是新疆通向祖国内地的要道,有“西域襟喉,中华拱卫”和“新疆门户”之称。地区矿产资源种类多、品位高、储量大,目前已探明各类矿种76种,占新疆已发现矿产
本文对井中激电地—井方式正反演进行了研究。文章以井中激电地一井方式为例,讨论了在直流点电流源的作用下,地下半空间存在各种体极化异常体时,不同电极方位测量井中二次场的异
本文的核心就是对分钟值的潮汐分析方法的研究。在此背景之下,从以下几个方面进行了探讨: PET的观测输出数据的是秒采样的数据,因此首先建立了以切比雪夫多项式和最小二乘原
随着找矿程度的不断深入,浅部易于发现和利用的矿产日趋减少,对寻找深部矿藏的快速定位预测方法的研究已构成当今整个地质界的重要研究课题之一。重磁勘探法是矿产资源勘查及
在领导实践中,几乎每一个新任领导者都会遇到这样的问题:如何对待前任领导干部及他们的政绩,怎样在前任的基础上,创造性地开展工作。我认为,这既是思想方法问题,也是从政道德
各位领导,各位代表,同志们,中国计算机自动测量与控制技术协会第五次代表大会,回顾总结了协会五年来的工作,特别是总结了第四届理事会的工作,确定了协会今后的工作目标和任务
本文利用国际地震中心(ISC)提供的高精度地震资料,相应的哈佛大学CMT机制解及其聚类分析结果,用三维空间图象和二维剖面相结合的方法,仔细研究帕米尔一兴都库什地区地震的空间分
随着近年来空间探测技术的快速发展,借助科学卫星、地面大功率雷达等手段,通过国际间的广泛合作,人们可以获得更为全面、准确的空间科学数据,在高纬电离层形态学,磁层粒子沉
有关青藏高原的隆升和变形机制仍存在诸多争议,作为青藏高原生长变形前缘带,青藏高原东南缘是研究高原隆升与变形的理想地区。为此,我们利用ChinArray台阵项目在青藏高原东南缘3