P波入射反射、透射系数推导

P波入射Zoeppritz方程的推导根据弹性力学的假设，介质是均匀各向同性的无限大介质，平面波是一种最简单的波动形式，其以波面为平面的形式在介质中传播，即平面波在垂直于波传播的任一平面上，各点的振动是同相的，实际上并不存在激发平面波的震源，所以它是一个数学抽象了的波动过程。点震源激发的球面波向四面八方传播，当其距震源足够远时，在这个地方研究一个局部的等相位面，可以将其看成一个平面波。在理论上，任何类型的波都可以用平面波的合成形式来表示，所以平面波是波动现象中最基本的形式，也是理论研究和实际应用的基础。在地震勘探中，讨论在两种不同的介质分界面上的波的传播现象是十分重要的。一般分为两种情况进行讨论，第一种，我们所研究的地球介质按其物性变化是分层的，具有层装结构。因此，讨论两种弹性性质不同的介质分界面上波的传播情况。第二种，地球表面是一个特殊的分界面，它将无限介质划分为两个半空间。地面以上的空气介质，其密度与地面以下的岩石或海平面以下的海水层及岩石层的密度相比可以忽略。因此，地球表面可以看成是一个弹性半空间表面，称为自由面，其上的应力作用为零。根据本文所讨论的地质模型所涉及到的地质灾害，我们只讨论波在第一种介质分界面情况下波的传播，即平面波在弹性分界面上的反射与透射。1.1波函数设有一平面谐纵波入射到两种半无限弹性介质的分界面上。在这种情况下，波不仅会折回到入射介质中传播，而且会透射到另一种介质中传播；即同时存在反射波和透射射波。反射波和透射波中都包含纵波和横波两种成份。P波在介质分界面上的反射和透射情况如图所示：关于位函数我们首先看：沿任意方向传播的平面波。设N是一个任意取定的单位方向矢量。Nlimjnk（1）下面来看沿N方向的平面波，或称三维平面波的波函数形式。三维平面波的波函数f满足三维波动方程，即：2222222221ffffxyzVt（2）这里我们通过和一维平面波函数类比，可以得出三维平面波函数的形式。我们知道，在一维平面波的情况下，空间任意一点,,xyz上的波函数值只取决于x。于是沿x正方向传播的平面波的波函数为1(,)()fxtfxVt。其中的x实际上是从原点至,,xyz点所在波面的垂直距离，即00dxyz（一维平面波的传播方向的单位矢量为Ni。在三维平面波情况下，这一距离应为dlxmynz。因此，将一维平面波函数中的x以lxmynz代替应该可以得到三维平面波的波函数）即：1(,,,)()fxyztflxmynzVt（3）同一维平面波一样，式中的t为波沿N方向的传播时间。1()flxmynzVt代表一个沿N的正方向传播的平面波。同理，1(,,,)()fxyztflxmynzVt代表一个沿N的负方向传播的平面波，在一般情况下，沿任意方向N传播的平面波的波函数可写成：11(,,,)()()fxyztflxmynzVtflxmynzVt（4）1.2平面简谐波：平面简谐波是是波函数为简谐形式的平面波，也是数学上最容易处理的一种波。因此，在研究波的传播问题时经常使用简谐波假定。沿x正方向传播的平面简谐波的波函数可写成：0(,)cos()fxtfkxVt（5）或0(,)sin()fxtfkxVt（6）上面两式分别代表的是余弦形式和正弦形式的平面简谐波。我们最常使用的是指数形式的平面简谐波()0(,)jkxVtfxtfe（7）通过取上式的实部或虚部即可得到余弦形式或正弦形式的平面简谐波的波函数。上面各波函数中的0f称为波的振幅，因为波函数值总是在0f和0f之间变化。下面讨论波函数中其他各量的意义及它们之间的关系。为此，首先“固定”时间变量t以考查波剖面的情况。不难验证，2(,)(,)fxtfxtk（8）这表明，波剖面的值每隔2k距离重复一次。因此我们将这个量称为波长，记为，2k同时，把2k称为波数。可见波数就是2距离内所含的波长个数。再“固定”空间变量x以考查振动图的情况。容易看出，2(,)(,)fxtfxtkV（9）这说明，振动图的值每隔2kV时间重复一次。因此将这个量称作周期，记为T，2TkVV由此可见，周期即为波传播一个波长距离所用的时间。另外，22kVTVV其中1T和2分别为频率和圆频率。利用上面得到的各量之间的关系，可将平面简谐波的波函数写成如下等价形式：2()()00(,)jxVtjkxVtfxtfefe()()(2)000xjtjkxtjkxtVfefefe2()0xjtfe（10）沿任意方向Nlimjnk传播的平面简谐波的波函数可写为()0(,,,)jklxmynzVtfxyztfe()0xyzjkxkykzVtfe（11）因此二维平面波的波函数可以写成：(,,,)fxyzt=()xyjkxkytAe（12）我们可以写出入射P波、反射波P波、反射SV波、透射P波和透射SV波的位函数：ZXP(A1)SV(A3)P(A2)SV(A5)P(A4)ββαβ'α’介质1（λ1μ1ρ1Pv1Ps1）介质2（λ2μ2ρ2Pv2Ps2）α(1)(1)()(1)1xzjkxkzwtAe（13）(2)(2)()(2)2xzjkxkzwtAe（14）(3)(3)()(3)3xzjkxkzwtAe(15)(4)(4)()(4)4xzjkxkzwtAe(16)(5)(5)()(5)5xzjkxkzwtAe(17)上式中1(1)(2)/sinxxpkkwv，1(3)/sinxskwv，2(4)'/sinxpkxwv，2(5)'/sinxskwv（18）且有(1)(2)(3)(4)(5)xxxxxkkkkk（19）由此可得反射和透射定律（斯奈尔定律）如下：1122''/sin/sin/sin/sinpspsvvvv（20）另外，由图可见：1(1)/coszpkwv，1(2)/coszpkwv，1(3)/coszskwv，2(4)'/coszpkwv，2(5)'/coszskwv在介质I中，总的位函数为(1)(1)(2)(2)()()(1)(2)112xzxzjkxkzwtjkxkzwtAeAe（21）(3)(3)()(3)13xzjkxkzwtAe（22）在介质中，总的位函数为(4)(4)()(4)24xzjkxkzwtAe（23）(5)(5)()(5)25xzjkxkzwtAe（24）1.3边界条件我们知道，介质分界面处的边界条件为位移连续和应力连续。因此，可写出本问题的边界条件如下：在Z=0处12121212()()()()zzzzzxzxuuww（25）（1）位移连续：地震波在传播过程中质点振动的位移u可以分解为其标量位的梯度与与其矢量位的旋度之和的形式，有：ugradrot（26）同时uuivjk（27）设xyzijk（28）将式（26）按梯度和旋度公式展开，得到u的3个分量为：yzxzyxuxyzvyzxwzxy（29）研究空间传播的平面波时，一般情况下选择直角坐标系，可使得波前面与一个坐标轴（如y轴）平行，此时方向余弦cos0。这样，波前面在y轴方向上无限延伸，波函数与坐标y无关，于是有0y此时，式（29）中对y的导数项变为0，则式（29）变为：yxzyuxzvzxwzx（30）这说明位移分量可以分为两部分其中一部分时位于xz平面内的位移分量u和w，它们只与和y有关，含有P波和SV波成份；另一部分是垂直于xz平面的位移分量v，它只与x和z有关。且只含有SH波成份。这一结果表明，可将P波和SV波作为一组与SH波分开来处理。我们在讨论P波和SV波时使用位函数和z然后由（30）式过渡到位移。为简单起见，记y。uxzwzx（31）和满足下面的波动方程：2222222211psVtVt（32）（2）应力连续首先由虎克定律有：2zzzze（33）zxzxe（34）虎克定律阐述了应力和应变的关系。再看应变的定义式：zzwez（35）zxxzwueexz（36）应变的定义式阐述了应变和位移的关系。再由位移和位移位的关系式：xzuxzvzxwzx（37）体应变的关系式：xxyyzzuvweeexyz（38）222psvv（39）2sv（40）由以上各式可得到：222222zzzzpssuvwwevvvxyzz将（37）式代入上式得到：22222222222xzzzpsvvxxzyzxyzxz22222svzzx式中：223()xyzyzxxyz，223()zxyxyzxyz故222222222222zzpssvvvxzzzx2222222xyz而0y故22222xz所以222222222zzpssvvvzzx而22221pVt（波函数满足波动方程）故22222222222pszzspvvvvtzzx（41）()zxzxwuueuxz2222222()svxzxxzz=222222(2)svxzxz（42）1.4反射系数和透射系数以下的工作是使波函数满足上面的边界条件，为此将（21）~（24）式代入（25）式，并整理。首先代（25）式的第一式有：(1)(1)(2)(2)(3)(3)()()()(1)(2)(3)12x3xzxzxzjkxkztjkxkztjkxkztxzAekAekAek(4)(4)(5)(5)()()(4)(5)45xzxzjkxkztjkxkztxzAekAek由于(1)(2)(3)(4)(5)xxxxxkkkkk故上式变为：(1)(2)(3)(4)(5)(1)(2)(3)(4)(5)12x345zzzzzjkzjkzjkzjkzjkzxzxzAekAekAekAekAek将1(1)(2)/sinxxpkkv，2(4)'/sinxpkv，1(3)/coszskv，2(5)'/coszskv，并且z=0代入上式：111111cos.cos.cos.123sinsincospps''2222cos.cos.''45sincosspjzjzvvpsAeAevv1122''12345sin111()cossincospspsAAAAAvvvv（43）代入（25）式的二式有：(1)(1)(2)(2)(3)(3)()()()(1)(2)(3)123xzxzxzjkxkztjkxkztjkxkztzzxAekA