下载文档原格式
地震波数值模拟方法研究综述在地学领域，对于许多地球物理问题，人们已经得到了它应遵循的基本方程(常微分方程或偏微分方程)和相应的定解条件，但能用解析方法求得精确解的只是少数方程性质比较简单，且几何形状相当规则的问题。
对于大多数问题，由于方程的非线性性质，或由于求解区域的几何形状比较复杂，则不能得到解析解。
这类问题的解决通常有两种途径。
一是引入简化假设，将方程和几何边界简化为能够处理的情况，从而得到问题在简化状态下的解答。
但这种方法只是在有限的情况下是可行的，过多的简化可能导致很大的误差甚至错误的解答。
因此人们多年来寻找和发展了另一种求解方法——数值模拟方法。
地震数值模拟(SeismicNumericalModeling)是地震勘探和地震学的基础，同时也是地震反演的基础。
所谓地震数值模拟，就是在假定地下介质结构模型和相应的物理参数已知的情况下，模拟研究地震波在地下各种介质中的传播规律，并计算在地面或地下各观测点所观测到的数值地震记录的一种地震模拟方法。
地震波场数值模拟是研究复杂地区地震资料采集、处理和解释的有效辅助手段，这种地震数值模拟方法已经在地震勘探和天然地震领域中得到广泛应用。
地震数值模拟的发展非常迅速，现在已经有各种各样的地震数值模拟方法在地震勘探和地震学中得到广泛而有效的应用。
这些地震波场数值模拟方法可以归纳为三大类，即几何射线法、积分方程法和波动方程法。
波动方程数值模拟方法实质上是求解地震波动方程，因此模拟的地震波场包含了地震波传播的所有信息，但其计算速度相对于几何射线法要慢。
几何射线法也就是射线追踪法，属于几何地震学方法，由于它将地震波波动理论简化为射线理论，主要考虑的是地震波传播的运动学特征，缺少地震波的动力学信息，因此该方法计算速度快。
因为波动方程模拟包含了丰富的波动信息，为研究地震波的传播机理和复杂地层的解释提供了更多的佐证，所以波动方程数值模拟方法一直在地震模拟中占有重要地位。
数值地球物理学中的地震波传播模拟地震波传播模拟是数值地球物理学中的重要研究领域，它对于理解地震的机制、预测地震的危险性以及构建地震防灾体系具有重要意义。
本文将从数值地球物理学的角度，探讨地震波传播模拟的原理、方法和应用。
地震波传播模拟是通过计算机模拟地震波在地下介质中的传播过程，以获取地震波的传播路径、传播速度和传播强度等信息。
在地震波传播模拟中，地下介质被划分为离散的网格单元，每个网格单元的物理性质（如密度、速度等）被赋予数值，通过求解弹性波方程，可以模拟地震波在地下介质中的传播。
地震波传播模拟的方法主要有有限差分法、有限元法和谱元法等。
其中，有限差分法是最常用的方法之一。
它将地震波方程离散化为差分方程，通过迭代计算，逐步推进地震波的传播。
有限差分法具有计算效率高、适用范围广等优点，因此在地震波传播模拟中得到广泛应用。
地震波传播模拟在地震学研究中有着广泛的应用。
首先，地震波传播模拟可以帮助我们理解地震的机制。
通过模拟地震波在地下介质中的传播过程，可以揭示地震波与地下介质的相互作用，进而深入研究地震的发生机制和演化规律。
其次，地震波传播模拟对于预测地震的危险性具有重要意义。
通过模拟地震波在地下介质中的传播，可以预测地震波在不同地点的传播路径和传播速度，进而评估地震对人类和建筑物的影响。
这对于制定地震防灾措施、规划城市发展具有重要的指导意义。
此外，地震波传播模拟还可以用于地震勘探和地震监测。
在地震勘探中，通过模拟地震波在地下介质中的传播，可以预测地下构造和地质条件，为油气勘探和矿产资源勘探提供重要信息。
在地震监测中，通过模拟地震波在地下介质中的传播，可以解释地震观测数据，帮助我们更好地理解地震活动的特征和规律。
然而，地震波传播模拟也面临一些挑战和困难。
首先，地震波传播模拟需要大量的计算资源和时间。
由于地下介质的复杂性和地震波的多尺度特性，地震波传播模拟需要高分辨率的模型和大规模的计算，这对计算机的性能和存储能力提出了很高的要求。
地震波波动方程数值模拟方法地震波波动方程数值模拟方法主要包括克希霍夫积分法、傅里叶变换法、有限元法和有限差分法等。
克希霍夫积分法引入射线追踪过程，本质上是波动方程积分解的一个数值计算，在某种程度上相当于绕射叠加。
该方法计算速度较快，但由于射线追踪中存在着诸如焦散、多重路径等问题，故其一般只能适合于较简单的模型，难以模拟复杂地层的波场信息。
傅里叶变换法是利用空间的全部信息对波场函数进行三角函数插值，能更加精确地模拟地震波的传播规律，同时，利用快速傅里叶变换(FFT)进行计算，还可以提高运算效率，其主要优点是精度高，占用内存小，但缺点是计算速度较慢，对模型的适用性差，尤其是不适应于速度横向变化剧烈的模型.波动方程有限元法的做法是：将变分法用于单元分析，得到单元矩阵，然后将单元矩阵总体求和得到总体矩阵，最后求解总体矩阵得到波动方程的数值解；其主要优点是理论上可适宜于任意地质体形态的模型，保证复杂地层形态模拟的逼真性，达到很高的计算精度，但有限元法的主要问题是占用内存和运算量均较大，不适用于大规模模拟，因此该方法在地震波勘探中尚未得到广泛地应用。
相对于上述几种方法，有限差分法是一种更为快速有效的方法。
虽然其精度比不上有限元法，但因其具有计算速度快，占用内存较小的优点，在地震学界受到广泛的重视与应用。
声波方程的有限差分法数值模拟对于二维速度-深度模型，地下介质中地震波的传播规律可以近似地用声波方程描述：)()(2222222t S zu x u v t u +∂∂+∂∂=∂∂ （4-1） (,)v x z 是介质在点（x , z ）处的纵波速度，u 为描述速度位或者压力的波场，)(t s 为震源函数。
为求式（4-1）的数值解，必须将此式离散化，即用有限差分来逼近导数，用差商代替微商。
为此，先把空间模型网格化（如图4-1所示）。
设x 、z 方向的网格间隔长度为h ∆，t ∆为时间采样步长，则有：z∆,i j1,i j +2,i j+1,i j-h i x ∆= （i 为正整数）h j z ∆= (j 为正整数)t n t =∆ (n 为正整数)k j i u , 表示在(i,j)点，k 时刻的波场值。
地震波传播模拟中的数值方法一、引言对地球上发生的自然灾害进行研究和预测一直是人类所探究的课题之一。
其中，地震是一种造成极大灾害的自然现象，它的预测和探测对减轻地震对社会影响，提高人类对灾害的应对能力，具有重要意义。
地震波传播模拟是地震研究领域的重要课题，为了更好地预测地震和应对地震灾害，需要对地震波传播的数值模拟方法进行深入研究。
二、地震波传播数值模拟的方法1. 有限差分法（FDTD）有限差分法，英文全称为Finite Difference Time Domain，是一种常用的求解电磁场和声场传播问题的数值方法。
FDTD方法利用有限差分逼近微分算符，将偏微分方程离散化，然后通过差分方程组求解离散化问题。
FDTD方法的优点是较为简便和直观，对于一些基础场问题可以精确求解，但是FDTD方法在离散化问题域时会导致误差，对于具有复杂形状、边界不规则和含有多个介质的问题，其求解需要繁琐的预处理工作和较为复杂的网格划分，求解过程也较为复杂。
2. 有限元法（FEM）有限元法，英文全称为Finite Element Method，是一种广泛应用于工程和科学计算领域的数值方法。
它是通过将一个复杂的问题域分解成多个小问题域，用简单的数学公式在每个小问题域内求解，通过对这些小问题域的求解累加得到整个问题域的解。
FEM方法的特点是能够对不规则的计算域进行处理，求解过程较为直观和简单，对于多介质、弹性、非线性等问题也有很好的处理能力。
但FEM方法对于较为复杂的问题各向异性和自由面的处理比较困难。
3. 间接边界积分法（BEM）边界积分法，英文全称为Boundary Element Method，是近年来发展起来的一种求解偏微分方程的数值方法。
BEM方法将待求解的域分为界面和域外两部分，通过界面上的边界积分求解内部问题。
BEM方法对于不规则和异形问题的边界条件求解有很好的处理能力，并且具有较高的精度和较低的计算量。
但是对于非线性问题处理不够准确，对纯内部问题的求解效果不如其他方法。
汶川地震的发生对周围断层影响的数值模拟研究李玉江，陈连旺，陆远忠，詹自敏中国地震局地壳应力研究所，100085[ 摘要 ] 基于川西-藏东地区三维粘弹性有限元模型，以GPS和震源机制解作为约束条件，利用汶川地震同震静态位移量结果，在考虑地形差异和粘弹性松弛等因素的影响下，研究汶川地震的发生造成周围主要断层库仑破裂应力变化。
研究结果表明：1）汶川地震的发生除造成震源区同震库仑破裂应力减小外，还造成断裂带两端不同程度的应力增加，尤其是断裂带北东段，这与震后较强的余震分布基本吻合。
断层面上库仑破裂应力的变化与静态位移量的分布存在一定的关系，静态位移量的高值区对应库仑破裂应力减小的区域。
2）汶川地震的发生显著增加了鲜水河断裂北西段、东昆仑断裂、龙日坝断裂、岷江断裂以及虎牙断裂库仑应力水平，将提升这些断层面上发震概率。
3）鲜水河断裂带作为菱形块体北部的一条强震活动频发区域，受汶川地震影响较为显著，地震的发生造成断裂带强震复发间隔缩短约52~104a，是未来值得重点关注的强震危险区。
[ 关键词]川西藏东地区，龙门山断裂带，汶川地震，库仑应力变化，数值模拟Numerical Simulation Research on effect of the surrounding faults induced by the Ms8.0 Wenchuan earthquake，Sichuan，ChinaLI Yu-Jiang，CHEN Lian-Wang，Lu Yuan-Zhong，ZHAN Zi-MinInstitute of Crustal Dynamics，CEA，Beijing 100085，China[ Abstract ] In this paper，with GPS and focal mechanism constraints，utilizing the findings of co-seismic static slip distribution，and taking the terrain variance and post-seismicrelaxation effect into consideration，we employed the three dimensional viscoelasticmodel of western Sichuan and eastern Tibet regions to calculate the Coulomb failure stressaccumulation rate on the major active faults and the dynamic stress change caused by theWenchuan earthquake. Our results indicated that：1) The Wenchuan earthquake not onlycaused the Coulomb stress decrease in the source region，but also the stress increase inthe two-terminal，easpecially the northeastern segment，and it’s comparatively consistentwith the aftershock distrbution.Meanwhile，the high concentration area of the static slipdistrbution was corresponding to the Coulomb stress reduction.2)The Coulomb stresschange caused by Wenchuan earthquake showed significant increase on five major faultsegments，which were northwestern segment of Xianshui fault、Eastern Kunlun fault、Longriba fault、Minjiang fault and Huya fault respectively，and these faults need to bepreferentially concentrated in the future.3) As the northern boundary of theSichuan-Yunnan rhombic block，there were frequently strong earthquakes along theXianshuihe fault，also the effect derived from Wenchuan earthquake was distinct，andthe characteristic earthquake recurrence interval of this fault was delayed 52~104a，so itwas ranked the high seismic hazard region in the future.[ Keyword ] Western Sichuan and Eastern Tibet regions，Longmenshan fault zone，Wenchuan earthquake，Coulomb stress change，Numerical simulation1前言2008年5月12日汶川Ms8.0级地震发生在南北地震带中段、青藏高原东缘的NE向龙门山断裂带上。
地震波数值模拟与分析地震波是地震活动中最重要的研究对象之一。
而地震波数值模拟和分析则是地震学领域中的重要研究方向之一。
在地震波数值模拟和分析的过程中，人们可以通过计算机模拟地震波的传播过程，并从中获取有关地震特征及其引起的地表破坏和建筑物结构变形等各种信息。
这对于地震灾害的预防、预测和减轻有着重要的意义。
地震波的数值模拟方法主要有有限差分法、有限元法、边界元法和谱元法等。
其中，有限差分法是目前地震波数值模拟中应用最为广泛的一种方法。
有限差分法在解决非线性、多维度和非静态问题方面表现尤为出色。
其基本思想是将地震波场离散成网格，并利用二阶精度差分公式计算各个时刻在网格点处的地震波场值。
有限差分法的优点在于精度高、计算速度快，同时可以对复杂地质构造及其他复杂条件进行模拟分析。
地震波的数值分析方法主要有PTA和TFI等。
其中，PTA是计算地震波传播中频谱组成的一种方法。
PTA方法基于傅里叶变换，将地震波在频域中进行分析，主要考虑波振幅和频率之间的关系。
通过对地震波的频谱进行分析，可以得出波传播路径、应变速率及层间的速度等信息。
而TFI则是通过时间域内的雷克子波分析地震波的能量分布，从而得出地表加速度和地震破坏信息。
当我们研究地震波数值模拟的同时，还要重视地震波分析的意义。
地震波的分析能够帮助我们对地震发生的原因、机制及它们对地表的影响进行研究。
同时，地震波分析也可以帮助我们评估地震对建筑物和基础设施的破坏。
这项工作通常涉及结构动力学模拟、震害评估、震害预测等研究领域。
此外，通过地震波分析，我们也可以了解地震所带来的生态影响和异常现象（如水波、地陷等）。
在地震波数值模拟和分析过程中，实际数据采集十分必要。
地震数据采集主要分为地震观测和近场强动观测两种方法。
地震观测是通过装置地震仪器等方法获得的数据。
而近场强动观测则是通过现场安装观测设备，获取地震波传播的信息。
同时，人工模拟地震波也是一种可行的方法，但其对于地震波的形态和波速等方面需进行较为精确的估计。
