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

1、这里我们通过和一维平面波函数类比,可以得出三维平面波函数的形式。我们知道，在一维P波入射Zoeppritz方程的推导根据弹性力学的假设，介质是均匀各向同性的无限大介质，平面波是一种最简单的波动形式，其以波面为平面的形式在介质中传播，即平面波在垂直于波传播的任一平面上，各点的振动是同相的，实际上并不存在激发平面波的震源，所以它是一个数学抽象了的波动过程。点震源激发的球面波向四面八方传播，当其距震源足够远时，在这个地方研究一个局部的等相位面，可以将其看成一个平面波。在理论上，任何类型的波都可以用平面波的合成形式来 表示，所以平面波是波动现象中最基本的形式，也是理论研究和实际应用的基础。在地震勘探中

2、，讨论在两种不同的介质分界面上的波的传播现象是十分重要的。一般分为两种情况进行讨论， 第一种，我们所研究的地球介质按其物性变化是分层的，具有层装结构。因此，讨论两种弹性性质不同的介质分界面上波的传播情况。第二种，地球表面是一个特殊的分界面，它将无限介质划分为两个半空间。地面以上的空气介质， 其密度与地面以下的岩石或海平面以下的海水层及岩石层的密度相比可以忽略。因此，地球表面可以看成是一个弹性半空间表面，称为自由面，其上的应力作用为零。 根据本文所讨论的地质模型所涉及 到的地质灾害，我们只讨论波在第一种介质分界面情况下波的传播，即平面波在弹性分界面上的反射与透射。i.i波函数设有一平面谐纵波入射

3、到两种半无限弹性介质的分界面上。在这种情况下，波不仅会折回到入射介质中传播， 而且会透射到另一种介质中传播；即同时存在反射波和透射射波。反射波和透射波中都包含纵波和横波两种成份。p波在介质分界面上的反射和透射情况如图所示：4关于位函数我们首先看：沿任意方向传播的平面波。设N是一个任意取定的单位方向矢量。N = li mj nk（ 1）下面来看沿N方向的平面波，或称三维平面波的波函数形式。三维平面波的波函数f满足三维波动方程，即：(2)平面波的情况下，空间任意一点x,y,z上的波函数值只取决于 x。于是沿x正方向传播的平面波的波函数为f (x,t)二fi(x -Vt)。其中的x实际上是从原点至x

4、, y, z点所在波面的垂直距离，即d = x 0y 0z (一维平面波的传播方向的单位矢量为N =?。在三维平面波情况下，这一距离应为d = lx my nz。因此，将一维平面波函数中的 x以lx my nz 代替应该可以得到三维平面波的波函数)即：f(x,y,z,t) = f|(lx my nz-Vt)(3)同一维平面波一样，式中的t为波沿N方向的传播时间。f1 (lx my nz -Vt)代表一个沿N的正方向传播的平面波。同理，f ( x, y, z, 4)1 f ( l x m代表一个沿N的负方向传播的平面波，在一般情况下，沿任意方向N传播的平面波的波函数可写成：f (x,y, z,t

5、)二 (lx my nz -Vt) f1(lx my nz Vt)(4)1.2平面简谐波：平面简谐波是是波函数为简谐形式的平面波，也是数学上最容易处理的一种波。因此，在研究波的传播问题时经常使用简谐波假定。沿x正方向传播的平面简谐波的波函数可写成：f ( x, t> of cosk -(x Vt)( 5)或f(x,t)二 f°sink(x-Vt)( 6)上面两式分别代表的是余弦形式和正弦形式的平面简谐波。我们最常使用的是指数形式的平面简谐波f(X,t)二f0ejk(xt)( 7)通过取上式的实部或虚部即可得到余弦形式或正弦形式的平面简谐波的波函数。上面各波函数中的f°

6、称为波的振幅，因为波函数值总是在-f。和 f。之间变化。下面讨论波函数中其他各量的意义及它们之间的关系。为此，首先“固定”时间变量t以考查波剖面的情况。不难验证，2兀f (x ,t)二 f (x,t)( 8)k这表明，波剖面的值每隔距离重复一次。因此我们将这个量称为波长，记为k同时，把称为波数。可见波数就是2二距离内所含的波长个数。再“固定”空间变量x以考查振动图的情况。容易看出,f(x,烧f( x, t)(9)这说明，振动图的值每隔2 :时间重复一次。因此将这个量称作周期，记为kV2 二 T 二kV V由此可见，周期即为波传播一个波长距离所用的时间。另外,kZVT V V1其中和小=2二分别

7、为频率和圆频率。T利用上面得到的各量之间的关系，可将平面简谐波的波函数写成如下等价形式:j 丝(x_Vt)f(x,t) = f0ejWf0j()(10)沿任意方向N=link传播的平面简谐波的波函数可写为jk (lx my nz _Vt)f(x, y,z,t)二 f°ej (kxx kyy z-Vt)(11)因此二维平面波的波函数可以写成:jg 初一 t)(12)f (x,y, z,t)= Ae我们可以写出入射 P波、反射波P波、反射SV波、透射P波和透射SV波的位函数:：(1) _Aej(kxxkjzvt)(13): _ Aej(kx2) x 好旳(14),；(3) _ Aej(k

8、x3)x kz3)z_wt)(15):了._ Aej(kx4) x k!4) z _wt)(16). -I f (5) _ Aej(k'5) x kz5) z _wt)(17)上式中k：kx(W psi n，k：3= w/vs sin ：，kx4) x = w / Vp2 sina，kx5)=w/Vq sin 0(18)且有x(19)由此可得反射和透射定律(斯奈尔定律)如下：v/sin ：一 Vsi/Sin 1 二 v/sin :丄沧/sin 1'(20)另外，由图可见：kZ1)=w/vpcos，k；2)=w/vpcos二,k：3)=w/vs cos ：,1 11kZ4) =

9、w/vp cos ， kZ5) =w/vscos：(21)在介质I中，总的位函数为 =(1)fe曲"1亠旳 A2ej(kx2)x kz(2)wt)=A3ej(kx3)x：；kZ3)z_wt)(22)在介质丨丨中，总的位函数为-2 =Aej(k)(4)x kr)zjwt)(23)1.3边界条件_ Aej(kx5)x kZ5) z_wt)(24)我们知道，介质分界面处的边界条件为位移连续和应力连续。因此，可写出本问题的边界条U1 = U2件如下：在Z=0处W1 -W2f zz)1 =(Jz)2(zx)1 = ( zx)2(1)位移连续：地震波在传播过程中质点振动的位移度之和的形式，有：u

10、可以分解为其标量位的梯度与与其矢量位的旋(25)H = grad 仁亠 rot'，(26)4 呻 耳 耳同时 u = ui vj k设T'zk将式（26）按梯度和旋度公式展开，得到U的3个分量为:丹評z 評yu-cXdy沖、汀X Lz：z(27)(28)(29)wy xGZCX研究空间传播的平面波时，一般情况下选择直角坐标系，可使得波前面与一个坐标轴（如y轴）平行，此时方向余弦 cos ： = 0。这样，波前面在 y轴方向上无限延伸，波函数与坐标y无关，于是有-=0 y此时，式（29）中对y的导数项变为0,则式（29）变为:刖評yU 二.X；z(30)严XZ.Z .X汐丄訓yw

11、:Z:X这说明位移分量可以分为两部分其中一部分时位于x-z平面内的位移分量U和w，它们只与和y有关，含有P波和SV波成份；另一部分是垂直于 x - z平面的位移分量V，它只与-：x和有关。且只含有SH波成份。这一结果表明，可将P波和SV波作为一组与SH波分开来处理。我们在讨论P波和SV波时使用位函数;：和 z然后由（30）式过渡到位移。为简单起见，记'。(31)和t满足下面的波动方程:V；1Vs2?2 :(2)应力连续首先由虎克定律有:ZX = 'Qx虎克定律阐述了应力和应变的关系。再看应变的定义式:wdx cz应变的定义式阐述了应变和位移的关系。再由位移和位移位的关系式:X.

12、Z.:x(32)(33)(34)(35)(36)(37)体应变的关系式:(38)vp -2v；(39)(40)由以上各式可得到:ezZ将(37)式代入上式得到:(丹% = P(v：-2v；)rCX式中:(兰)y. z； y z :x:x：y：z/亠 2匕2+厂+ 2叫故二 ZZ 二Vp -2VS2.z2宀丄二二2 2 2.X2鋼-ZT矿0所以zz = P(vp -2v；严那 +2PVs2I2敖' 2 + _ _ l czczcx 丿1 戸2（p而12 =,（波函数满足波动方程）故二ZZ2Vp.:t2.Z2(41)zx uezx,:w=u(一:X.x :z2 2:X:X：z :Z(42)

13、=讼"GXCz CX1.4反射系数和透射系数以下的工作是使波函数满足上面的边界条件,为此将（21）（ 24）式代入（25）式,并整理。首先代（25）式的第一式有:Aej(kX°x 必亠.乍(1) . Aej(kX2)x kzz-t.*(2) _ Aej(kX3)X 1乍(3)_ Aej(kX4) X kz4) z ；t)k _ Aej(kX5)x kz5) z ；：；t)k(5)由于 kXkX2) *x3)二 kX4)二 kX5)故上式变为：从灣条 Aejkz2)zk 一 Aejkz3)zk=Aejk儿k- Aejkz5) zk(5)ix2x3z4x5z将kJ二 k：2)=

14、 /vp sin ：，k：4)=/vpsin ： '，k；3)= -/vscos k；5） /vs cos :'，并且z=0代入上式:.w. wjCOSQZj COSQZVpiVp.=A.e P1sin 二亠 A>eP1sin 二亠 "e上 cos|.”zVsi cos -VjO二 Aesin 二VP2j cos 典Vs2f A5e.z一COS：Vs?=_(A A2)丄cos ： A =VPiVsi1 ' 1sin： A -VP2V(43)VA代入(25)式的二式有：Ae j(kX°x kZ1* - t)k(1) . Aej(kX2)x kZ2

15、)z_ t)k(2) _ Aej(kX3)x kZ3)，t)k(3) 二 Aej(kx4)x kz4)z_. °k一 Aej(kx5)x kz5)- t)k(5)由于 kx1 kx2）二 kx3）二 kx4）二 kx5）故上式变为:Ae好k德尹诗-&宀kJk：3) = /vs sin ： , k：5) =/vs sin :,kJ(2)=/vpcos： , kz/vpcos：,kz3)-/vs1cos： , k；4) = /vP2cos：,kz5)F：/Vs2COS：，且 z=0 故:©©«fycos-： A?(cos： )-A sin-Vp1Vp

16、1Vs1二 A cos：Vp2*v;sin1 11Acosj A2cos：A sinVP1VP2乂 1人丄cos：VP2人丄si nVs?叱(A_A)+泄 A3 =丄 cos。VP2VP1(44)应力连续故代入（25）式第三式有:p1 严牙 2vs12e-21Vp1t: z2p22Vsvf空2(空2-2 22Vs2 ( - 2p2 : t: z-2 2丄）；XZ z£2 2 p/p1 _2Vs1 (Aej(kxF1)z%2 + 民丿炉川日材 2) +Vp12V2(AeJ(ki1)xkZn t)k2(1). Aej(kJ2)x-kZ2)t)k2(2).4泌3.必3*，t)k(3)r(3

17、)2 2 =订严 半(代ej(k'4) V2)2vs22(A4ej(k(4)x'kZ4)"t)k2Z4) +Vp2Aej（kX5）xkZ5）z-t）k（5）k（5）因为kX1 kX2)二 kX3)kX3)F：/VqSin ：,kX5)IZd =()/vs2 sin P , kz =co/vPlcosa ,kZ2)/VpCOS：,kZ3)-/vs cos 1S1kZ4) = /Vp2cos：',k；5) = - /vs cos J ,z=0 ,故上式变为:2 2Vp1 - 2vs122222(A '宀 ) 2vs1 (A cos :-Vp12 2pl -

18、 2Vs12VP12A $ 2A2 cos :-vpi2Vs2sin : cos :)2 2p2 - 2Vs22Vp22(A 2) 2vs22(A4cos2 ：Vp22A sin :' cos :')Vs222v 2vpis2Vp1(A A)2p p (A AOcos2 : -sin2 I'A?=-2 2 2 、 2V2Vs1sin ：(A A2)s2I± 卩12Vp1-2v： sin2 av¥ A cos.篇 sin 2IAsVp2A psin 2 ： A (45)代入（25）式的四式:汕2二厶L、 J '、L、JZ:xz Xvp1(2Ae

19、j(kX1)x'kZ1)t)kX1)kZ1) PA2ej(kX2)x kZ2)t)kX2)kZ2)-AekXx kz- t)k2(3) _ AeHkX3' kZ3)Z-，t)k2(3) )_r v? (2Aej(kX4)x kZ4)zt)kk Aej(kX5)x kZ5)z-丄乍?® _ Aej(kX5)x席*-丄乍2(5) 由于 kX1 - kx ( 2 kx- 3k)x= k： kJ= kf=/vp sin ：, k：3)=/vs sin ：,kX4) x -/ vp sin ： ,-/ vs sin : , kJ-/ vp cos：, k；2 -/ vp cos

20、：,k；3- - /Vscos： , k；4) = /vppcos , k；5) - - /vsp cos :,且 z=0，故上式变为：、2 - -2-'2* '2 2 - ''2 2 v (2A|p sin ： cos：-2APpsin: cos：A psin - A cos -)Vp1Vp1Vs1Vs1=*：(2人2sin ： cos：VppAr s in2Vsp:- A p cos :')VspHkA - A)丄Sin21 A2Vpi11sin2 ： - A 2 cos -)VslVs!2vS1(2A-2si n2：2Vp21 2 一A 2 sin

21、AsVs2-yCOS2 :)Vs2v：(AVP1sin2: -A2 2psin2工"Afe-ysir2ViV2 2 2 2P -人卑 c°s B)=v：(Asin2o( +A 罕siB -A第 c°sP )ViRVP2VVs222VPisin2：(几 - A2) cos2 : Asin 2 二 A12；1cos 2 : A5(46)联立(43) ( 44) ( 45) (46)有1 1 . sin ： (A A2)cos : A3VnV11.cos： (A -ADsin : APlVS12 2 . 2Vp - 2vs sin :-P1|(A A2) sin2 ：人

22、二VS11 cos : AVP2V.1 ' 1 -'cos： A sin AVP2V2P11 ' sin ： A42p2-Vsin2o(A| -A2) +cos2P 人=2 %VP1S22c 22丄 Vp2-2VqSin :2'鳥.'sin 2一：仏cos2： AVP2：?1(47)2VP2由斯奈尔定律可得：v： -2v： sin2 ：二 v： -2v； sin2 ： = v； cos2 :vp? -2v： sin2 ： ' = vp? -2vp2 sin2 ' = v： cos2 -代入（47）式中的第三式，并将其方程组的各项同除A，得sin：企 亞 cos-s in :匕込丁公二飞 in：A VqAVp2A1Vs2AA2VP1. r A +VP1丹代Vs甘A5cossincossincos -A1Vs,A Vp2A1Vs2A1cos2 ：仏 -sin 2 虫-Wcos2 ： “ -Wsin 2 ：-cos2 -AA 叱A 叱 A2V； sin2G A2 +cos2 0 A> +VP1AA22VP2sin 2A121VP1sin 2此方程组