汶川大地震的动力学成因
Hu, C., Y. Cai, and Z. Wang (2012), Effects of large historical earthquakes, viscous relaxation, and tectonic loading on the 2008 Wenchuan earthquake, Journal of Geophysical Research, 117, B06410, doi:10.1029/2011JB009046. (SCI, IF: 3.303)
d3 dx3
f (x) (x)4 4!
d4 dx4
f (x)
f (x x)
f (x) x d dx
f (x) (x)2 2!
d2 dx2
f (x) (x)3 3!
d3 dx3
f (x) (x)4 4!
d4 dx4
f (x)
（1） （2）
对于函数f(x)，通常意义下的导数（微商）定义为：
x
lim
f
(x dx) dx
f (x)
x
lim
f
(x) f (x dx) dx
x
lim
f
(x dx) f (x dx) 2dx
当dx→0时，以上三种形式都是微商的正确定义。 如果dx是有限的，如何给出微商的近似定义？
如果Δx很小，f(x)可微，则以上级数收敛。 次数越高，收敛级数的项的绝对值越小。
由（1）得到，
f (x x) f (x) x d f (x) (x)2 d 2 f (x) (x)3 d 3 f (x) (x)4 d 4 f (x)
dx
2! dx2
3! dx3
• 有限差分法以变量离散取值后对应的函数值来近似微分方程中 独立变量的连续取值。
• 我们放弃了微分方程中独立变量可以取连续值的特征，而关注 独立变量离散取值后对应的函数值。
• 有限差分法的具体操作分为两个部分： • (1)用差分代替微分方程中的微分，将连续变化的变量离散化，
从而得到差分方程组的数学形式； • (2)求解差分方程组。
由（4）给出导数的一阶精度（first order accurate）近似为：
d f (x) f (x x) f (x)
dx
x
（5）
（5）式称为向前差分格式（forward-difference formula）
由（2）式得到
f (x)
f (x-x) x d dx
f (x)- (x)2 2!
3. 如何数值求解差分方程组
2.2 网格剖分
• 网格剖分就是研究区域和边界的离散化 • 1.矩形分割 • 2.三角形分割 • 3.极网格分割
对地球物理问题的连续求解区域通过网格划分离散为空间上得一系 列网格点，接下来需要利用一定的差分格式对偏微分方程组中的导 数用差商进行近似，从而将偏微分方程组离散化为差分方程组。
d2 dx2
f (x) (x)3 3!
d3 dx3
f (x)- (x)4 4!
d4 dx4
f (x)
d f (x) f (x) f (x x) +O(x)
dx
x
由（7）式得到导数的另一个一阶精度近似：
d f (x) f (x) f (x-x)
dx
x
（8）
有限差分法的主要内容
1. 建立地球物理问题的离散有限差分模型 （1）如何根据问题的特点将定解区域做网格划分； （2）如何在所有网格节点上用有限差分格式对导数求近似， 对函数、初始条件和边界条件求近似； （3）如何把原方程离散化为代数方程组，即有限差分方程组。
2.从理论上研究有限差分模型的形态，以保证计算过程的可行性和计算结果的正确性 （1）解的相容性； （2）解的稳定性； （3）解的收敛性。
有限差分方法的基本特点 该方法是一种直接将微分问题变为代数问题的近似数值解法，数学 概念直观，表达简单，是发展较早且比较成熟的数值方法。
2.1 有限差分法基础
有限差分方法的基本原理
该方法将求解域划分为差分网格，用有限个网格节点代替连续的求 解域。有限差分方法以Taylor级数展开等方法，把控制方程中的导数 用网格节点上的函数值的差商代替进行离散，从而建立以网格节点 上的函数值为未知数的代数方程组。
4! dx4
d f (x) f (x x) f (x) O(x)
dx
x
（3） （4）
d f (x) f (x x) f (x) O(x)
dx
x
（4）
式中的O(x)项表示忽略掉的所有项中的最大项的量级 是Δx，也就是说，忽略掉这些项带来的误差中的最大 项和Δx成正比。
2.3 差分格式
用Taylor级数展开可以给出微商的近似形式。
ห้องสมุดไป่ตู้
dx变为Δx
对于连续函数f(x)，它在相邻点上的值f(x+Δx)和f(x- Δx)可以用Taylor 级数展开为
f (x x)
f (x) x d dx
f (x) (x)2 2!
d2 dx2
f (x) (x)3 3!
第二章 有限差分法
主讲人：胡才博 中国科学院大学地球科学学院 中国科学院计算地球动力学重点实验室
第二章 有限差分法
• 2.1 有限差分法基础 • 2.2 网格剖分 • 2.3 差分格式 • 2.4 差分方程 • 2.5 应用实例
解析方法的局限性 1. 地球内部介质，不仅存在纵向非均匀结构（一维地球模型）， 也存在横向非均匀结构（不同块体、断层系统）； 2. 几何模型也呈现出相当的复杂性； 3. 另外，边界条件和初始条件对于不同问题具有特殊性。
（8）式称为向后差分形式（backward-difference formula）。
（6） （7）
f (x x)
f (x) x d dx
2.1 有限差分法基础
对于存在复杂介质和几何、特殊边界条件和初始条件的实际地质问题， 一般不存在解析解，需要近似的数值求解方法。 有限差分方法是地球物理方法中最常见的一种。
有限差分方法（Finite Difference Method, FDM）是计算机数值模拟最 早采用的方法，至今仍被广泛使用。