地球物理勘探 1-5波动方程的解及地震波的特点

V P2 2 0
显 然 ，r
x2 y2 z2
轴之间的夹角
为矢量r和z轴之间的夹角，
为矢量r在 xoy平面上的投影与x
14
3.2 无限大、均匀各向同性介质中的球面波
各种算子在球坐标系中的表达式为：
gradu u r er 1 u r e 1 u r sin e u u 0
0 ( t ) ( t ) 0
t0 0 t t t t
12
（初始条件）
3.2 无限大、均匀各向同性介质中的球面波

当t<0时，整个空间位函数有：
x , y , z ,t
0, ( x , y , z , t ) t 0

在t=0时，点震源开始作用，作用时间为Δ t；
2
1 2u ) 2 2 r sin
2u r 2
)
2u r 2

2 u r r
15
3.2 无限大、均匀各向同性介质中的球面波
将各种算子带入纵波的波动传播方程，得到著名的弦方程：
2 1 t
2
V
2 P
2 1 r
2
0
1 r
其解为
r 0 W r 0 S
将左端一、二项带回可得 ： lim grad ndS
r 0 S
1 V
2 r 0 p
lim ( t )dW
W
1 V
2 p
1( t ) r
因为是纵波，只存在r方向的分量，即grad lim
r 0 S
, 则上式变为：


t>△t时，点震源作用完毕。
边界条件：因已假设弹性介质的空间是无限的，其内不 存在任何弹性分界面，故无边界条件。

球面纵波：在在均匀各向同性介质中激发点源，所产生 的胀缩力作用面具有球对称性，所产生的波前面是一个 球面，当研究任意一球半径r方向上的纵波的传播特点， 就可以代表其他方向的传播特点，称此为球面纵波。
A2 V 2 A2 0 A3 V 2 A3 0
10
无限大、均匀各向同性介质中的波动方程的解有两组。 第一组解：当 V V p ( 2 ) / 时，
U A1 exp( vw0 i V ( x V p t ))
沿x方向的位移分量不为零，其他方向的位移为零，即波的传播方向与位 移方向一致，所以称为平面纵波，也称为胀缩波，通常简称为P波。 第二组解：当 V V s / 时，
3 波动方程的解及地震波的特点
本章包括：

无限大、均匀各向同性介质中的平面波 无限大、均匀各向同性介质中的球面波 地震波的动力学特点 地震波的运动学特点
1
3 波动方程的解及地震波的特点
P波—波动方程 S波—波动方程 SV波 SH波
无限大、均匀各向同 性介质中的平面波
2
3 波动方程的解及地震波的特点
对于球面纵波，位移u 只存在 r方向上，即u 只是 ( r , t )的函数，则 u r er u r r r
拉普 拉斯 算子 ： 1 u 1 1 u 2u 2 ( r 2 ) (sin r r r r sin r u u 0 1 r2 ( 2r u r r
u0 v A2 exp[ i V i V ( x V s t )] ( x V s t )]
w A3 exp[
其位移方向与波的传播方向垂直，所以称为平面横波，也称为剪 切波，通常简称为S波。S波有两个质点振动方向：沿Z轴振动的S波分 量为垂直偏振剪切波，称为SV波，沿Y轴振动的S波为水平偏振剪切波， 11 称为SH波。
位移方程(P74)
胀缩点震源——球面纵波 震 胀缩力 物理含义
无限大、各向同性 介质中的球面波
源 性 旋转力
质
旋转点震源——球面横波
位移方程
物理含义
3
3 波动方程的解及地震波的特点
球面纵波的传播特点 视波长λ 波剖面 视波数k 振动图（实际记录） 视周期T 视频率f 地震波的动力学特点 能量密度 能量和球面扩散 能流密度 球面扩散 地震波的谱分析（傅立叶变换 ） 应用 识别不同的地震波 识别岩性 关 系
8
3.1
无限大、均匀各向同性介质中的平面波
若 k1 x k 2 y k 3 z Vt c 为常数，t固定，该方程代表一个 以 k k1 , k 2 , k 3 为法向量的平面，波在每个这样的平面上 必然有相同的相位，即平面波是垂直于 k1 x k 2 y k 3 z 平面传播的。不同的t，有不同的波前面。平面波的波 前面是平行的。 k1 , k 2 , k 3 是平面的法线方向数。有 k12 k 22 k 32 1 取负号时，表示随时间t的增加，波沿k方向前进，即延 迟一个时间。 取正号时，表示随时间t的增加，波沿-k方向前进，即 提前一个时间 当K是任意矢量的，则平面法向量为任意方向的。即表 示沿任意方向传播的平面简谐波。
r c ( t
1 1
r VP
) c2 ( t
r VP
)
如果使 t
r VP
为常数，则t 随 r 增大而增大，( t c
1
r VP
)代表了沿
r向“外”传播的波，称为发散波。t c (
2
r VP
)
代表了沿-r向“内”传
播的波，即向震源方向传播的波，称为聚会波。聚会波只存在于t为
3.2 无限大、均匀各向同性介质中的球面波
下面进一步讨论在地震勘探的初始和边界条件下，胀缩力divF和旋 转力rotF的作用下，求解波函数，并分析其性质（以纵波为主）
2.2.1

胀缩点震源条件下的球面纵波
1、初始和边界条件
初始条件：在均匀各向同性介质中，炸药爆炸后产生一个 均匀的力垂直作用在半径为a的球形空腔壁上。当 a 0 或 相对无限大空间而言，这个震源可以看成是点震源，其力 位函数或者震源函数可以表示为
) dS S
r ' lim 4 ( C 1 C 1 r 0 V 4C 1 带回原式，则得： C1( t r Vp ) 1 4V
2 p
)
1( t
r Vp 1
)
再将 C 1 带到特解式中，则得
（ r , t )
负值的情况，这与实际不合，则该波是不存在的。
16
因此，上式又可写为：


r
1

1 r
c1 ( t
r VP
)
该式是齐次方程的解，只反映了波的传播特点。当力位 函数不为零时，需求非齐次方程的解，即达朗贝尔解。
2 t 2 V p2 2 2 t 2 V p2 divgrad ( t )
Vt d Vt Vt Vt

x
将上式代入波的Navier方程
U ( ) grad F
2
2U t 2
整理简化，并令体力F=0，可得
2 A1 V 2 A1 0
r
dS
1 V
2 p
1( t )
再将特解带入左端项，则 有 ：
19
lim
r 0 S

' C1 C1 1 r ' dS lim 2 C1 dS lim 2 ( C 1 r0 r0 r rV p Vp r S r
7
3ຫໍສະໝຸດ Baidu1
无限大、均匀各向同性介质中的平面波
一、沿任意方向传播的平面波
直接用位移向量所表示的波动方程式求解
2i k1 x k 2 y k 3 z Vt d U A exp Ae i Acos i sin
式中：A为振幅，决定位移的大小，ψ 为波的相位. 2π f/V = w/V为简谐波参数，f频率，w圆频率，V波速。 i为虚数符号，仅考虑实数时为简谐波。 k1 x k 2 y k 3 z Vt c 为传播项。 此式表达的波函数为沿k方向传播的平面简谐波。
21
3.2 无限大、均匀各向同性介质中的球面波
3、球面纵波的位移解
球面纵波的位移为：
U p grad 1 4V
2 p
[
1 r
2
1( t
r Vp
)
1 rV p
(t
' 1
r Vp
)]
r r
)



其物理含义为： 球面纵波以 V P 的速度沿r方向向外传播 纵波质点的位移与震源强度及其一阶导数成正比 质点位移与波的传播方向一致 纵波质点位移与r及其平方成反比 t r / V p 表示延迟位 质点位移在一维空间振动，为线形极化波
2U t 2
两边取分别散度和旋度，并且令
VP2 ( 2 ) / VS2 /
则可得纵波方程和横波方程
2 t
2
VP2 2
2 t 2
V 2
6
3 波动方程的解及地震波的特点
波动方程反映了物体波动过程的普遍 规律。波动方程的求解通常是和定解问题 联系起来考虑。波动方程的解就是波函数。 在不同的情况下可以得到不同的解， 即波函数有不同的形式。
9
二、沿X轴方向传播的平面波（即
k1 1, k 2 0 , k3 0 2 i U A exp 2 i u A1 exp 2 i v A2 exp 2 i w A3 exp
k x
）
x x x
左端的第二项，按奥斯特洛斯公式：任意一个矢量场，若 在空间域W中该场的一阶导数是存在，则该场边界S上的 通量等于在W域中散度的体积分。即：
18
divgrad dW grad ndS
W S
由此可得出： lim divgrad dW lim grad ndS
W
17
左端的第一项，其特解带入则得，
lim
r0 W
2 t
2
dW lim
r0 W
1 2C 1 r t
2
dW
lim
r0
1 2C 1 r t
2
dW
W
dW
W
4 3
r 3
2C 1 4 3 lim ( r ) 2 r 0 t 3 0
将点震源用半径r=a的小球代替，小球体积为W。对上式 求体积分，并令r->0,其极限情况就是点震源的达朗贝 尔解。
lim
r 0 W
2 t 2
dW V p2 lim divgrad dW lim ( t )dW
r0 W r0 W r0
若令： 1 ( t ) lim ( t )dW
13
3.2 无限大、均匀各向同性介质中的球面波
2、坐标变换和球坐标下球面纵波的传播方程解
已知球面纵波传播波动方程如下：
2 t 2 此式是直角坐标系中的波动方程，需转换到球坐标系中，即
x r sin cos y r sin sin z r cos ( 0 r ,0 ,0 2 )
4
3 波动方程的解及地震波的特点
地震波的运动特点 惠更斯-夫列涅尔原理 射线积分理论-克希霍夫积分

费马原理和波的射线 时间场和视速度定理
5
3 波动方程的解及地震波的特点
在均匀、各向同性、理想弹性介质中的弹性波方程为：
U ( ) grad F
2
1 r
C1( t
r Vp
)
4 rV
2 p
1( t
r Vp
)
20
3.2 无限大、均匀各向同性介质中的球面波
力位函数不为零的波动方程的达郎贝尔解为： 1 ( r ,t ) 1( t r / Vp ) 2 4 rV p 该式为用震源函数表示的波动方程的位移位解。 在实际工作中，人们不可能接收到质点的位移 位，而只能接收到质点的位移。 地震记录上地震波的振幅A值就是反映质点的位 移。所以必须把位移位转换成位移。