地震波的吸收.docVIP

中国地质大学（北京） 应用地震学 — §1.3(C) 备课教师：段云卿 五、 地震波的吸收 对吸收的概括说明 弹性介质是实际介质的近似，实际介质对地震波有吸收。使波传播时能量 很快衰减。地震波的能量不可逆地转化成了热能散发掉叫做吸收。 目前公认岩石颗粒之间的内摩擦力是吸收的主要原因，这种内摩擦力也叫粘滞力。 不同的介质中，内摩擦力所遵守的规律不同，其中沃尹特（ Voigt ）的假设比较有代表性。 2．沃尹特（ Voigt ）假设下的波动方程 （ 1）沃尹特假设 前面讨论的完全弹性介质中的波动方程是在应力与应变成正比 （虎克定律） 的条件下得到的。 沃尹特假设：应力与应变的关系包括两部分，一部分是应力与应变成正比的满足虎克定律的应变，另一部分是应力与应变的时间变化率成比例的粘滞效应。 （ 2）沃尹特假设的粘弹介质中的波动方程 2 u ( ) 2 u 1 t 2 u (6.1-45) t 2 3 t grad x i j k y z 是粘滞系数。 对（ 6.1-45 ）求散度，可得纵波的波动方程为： 2 4 ] 2 t 2 [( 2 ) (6.1-46) 3 t 对（ 6.1-45 ）求旋度，可得横波的波动方程为： 2 2 rot u t 2 rot u ( ) (6.1-47) t 沃尹特粘弹介质中一维谐波方程的求解 1 中国地质大学（北京） 应用地震学 — §1.3(C) 备课教师：段云卿 （ 1）求解 分析平面谐波沿 X 方向传播，以纵波为例。 假设纵波的位移位是： j x ) (t j ( t Kx) (x, t) 0e v 0e (6.1-48) u v w u （注：波沿 x 方向传播， v =0 ， w =0） x y z x y z 又 u ( 注： u grad i 0 j 0k ui v j wk ) x x u 2 K 2 x2 (6.1-49) x 将(6.1-49) 代入 (6.1-46) ，得 2 ( 2 )K K 2 ( 2 )  2 2  j 4 K 2 3 2 4 4 （ ） (6.1-50) (2 ) j 3 j 3 将上式有理化，得 K 2 2 2 2 )2 1 { 1 j 2 )2 1 } [( 2 2 ] 2 [( 2 )2 2 2 ] 2 [( 2 2 ] 2 (6.1-51) ( 2 ) 2 2 2 2 [( 2 ) 2 22 1 ] 2  β 2 cos [( 2 ) 2 2 2 1 sin ] 2 tan 2 2 中国地质大学（北京） 应用地震学 — §1.3(C) 备课教师：段云卿 1 tan ( ) 这样， (6.1-51) 可写成 K 2 2 j tan 1 ( [( 2 ) 2 2 2 ] 1 2 exp[ 2 )] （ e j cos j sin ej cos j sin ） 2 4 1 4 j 1 tan 1( K [ ] exp[ )] 2 ) 2 2 2 ( 2 2 或 K [ 2 4 2 ] {cos[ 1 tan 1( )] j sin[ 1 tan 1( )]} (6.1-52) 1 4 ( 2 ) 2 2 2 2 令 R 表示上式中的实数部分， α表示虚数部分， 2 4 1 4 1 1 则 R [ ( )] (6.1-53) 2 ) 2 2 2 ] cos[ tan ( 2 2 2 4 1 4 1 tan 1 ( [ 2 ) 2 2 2 ] sin[ )] (6.1-54) ( 2 2 于是 K R j (6.1-55) 将(6.1-55) 代入 (6.1-48) ，得 (x, t) 0 e xe j ( t Rx) (6.1-56) 对上述求解过程的大致说明 假设的解为 (x, t) 0 e j ( t Kx ) （6.1-48 ） 真正的解为 ( x, t) 0 e xe j ( t Rx) （ 6.1-56 ） 二者的差别在于衰减项 e x ，这正是吸收项。 吸收作用的实测结果是随传播距离的增加振幅呈指数衰减，我们正希望最 终的解当中出现 e x 。 分析（ 6.1-48 ）只有 K 变成复数才能出现 e x 项，可以认为沃尹特已经按 这套思路去处理了。 另一点要说明的就是：沃尹特假设下的波动方程（ 6.1-46 ）无法求解。所 以沃尹特先假设解为（ 6.1-48 ）的形式，然后代入（ 6.1-46 ）来确定 K 值。最 3 中国地质大学（北京） 应用地震学 — §1.3(C) 备课教师：段云卿 终确定的 k R j 达到了开始的设想。 3） 沃尹特介质中的波动方程解（ 6.1-56 ）的物理意义①有频散 波在沃尹特介质中传