地震波的吸收.docxVIP

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