平面波反射与透射ppt课件.pptx

《平面波反射与透射ppt课件.pptx》由会员分享，可在线阅读，更多相关《平面波反射与透射ppt课件.pptx（69页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第七章第七章 平面波的反射与透射平面波的反射与透射 在无限弹性介质中有无旋波（纵波）和等容波（横波）这两种弹性波的传播。 但是实际介质内部，一般都会存在分界面。如地表表面、岩层的分界面。 通常把实际地层称为层状介质。本章研究弹性波在传播的过程中遇到分界面的情况。 纵波和横波 地震波到达之处,介质就产生形变。由力学定律知道，任何小的形变都可以分解为两部分：一部分表示胀缩，即变体积而不变形状；另一部分表示畸变，即变形状而不变体积。 形变传播时，两部分的传播速度不同。在震源附近，两部分还未分开，所以波经过处的形变是复杂的。在较远的地方，波阵面就分成两个。胀缩波传播较快，波阵面上的质点位移和传播方向一

2、致，所以叫做纵波，一般用字母P(Primary)表示。较慢的叫畸变波,质点位移和传播方向垂直，所以叫做横波,一般用字母S(shear)表示。在均匀的介质中，波阵面在震源附近是曲面，但在相当距离后就趋近于平面。为简单起见，现只讨论平面波。 地震地震P波波(纵波纵波)和和S波波(横波横波)运行时弹性岩石运动的形态运行时弹性岩石运动的形态 纵波：间隔形成压缩带（密集带）和膨胀带（稀疏带），传播方向与振动方向一致，波速 Vp横波：传播方向与振动方向相垂直，波速Vs水平面内分量：称SV波 垂直面内分量：称SH波 7-1平面波在自由界面上的反射平面波在自由界面上的反射自由界面：指地表应力为零的界面，半无限

3、弹性体的界面就是自由界面。由于地球表面大气压相对于地球内部压力来说是十分小的，在讨论中可把大气压忽略不计，于是地球表面可以看作自由界面。三点说明 因为地球介质在短暂力（如爆炸）的作用下，在离开震源稍远的大部分地区 可看成弹性体，并且地球半径比地震波波长大得多，所以可将地球看作半无限大弹性介质，同时，可近似将地震波视为平面波。 任何复杂的波都可看成一系列不同振幅，不同频率及波长谐波的叠加，因此仅讨论一个平面谐波入射到自由界面的情形即可。 如果平面波的传播方向与z轴垂直（即在xoy面或平行于xoy的平面内），弹性动力学问题中的场变量都依赖于x和y两维，此时弹性动力学问题为二维问题，这 时讨论分层介

4、质波传播问题就是讨论这样二维弹性动力学问题。7-1平面波在自由界面上的反射平面波在自由界面上的反射一、平面纵波在自由界面上的反射一、平面纵波在自由界面上的反射 设半无限弹性介质的自由表面为yoz面，z轴与图面垂直。假定yoz面的左边为真空，由于没有传播振动的介质，故不会产生透射问题，全部入射波都在界面上被反射。P1入射纵波P1S1反射横波P11反射纵波设一平面纵波与x轴夹角为1的方向入射到自由界面；设入射纵波中质点的位移函数（即波动方程的解）为：22222pUUVtr1111cossinsin()pxyUAtV1111sin()UAtf xg y改写为1111cossin;ppfgVV 设入射

5、波为拉伸波，即质点的运动方向与波的前进方向相反，于是可以知道在入射波中质点的位移分量为：111111cos,sinuUvU 假设与自由界面作用后，只有纵波反射，且反射纵波的位移函数为：22221sin()UAtf xg y2222cossin;ppfgVVf2前面的负号原因：反射波相对于x轴而言，是正向传播。1为常数，表示由反射所引起的相位改变。仍假设反射波为拉伸波，在此波动中质点的位移分量为：222222cos,sinuUvU 问题：确定反射波的各参数问题：确定反射波的各参数方法：由边界条件可以确定方法：由边界条件可以确定边界条件：在自由界面上位移不受任何限制，即应力等于零；即在x=0的平面

6、上，对于y与t的任意值有0 xxyxz由广义胡克定理和几何方程可以得到：0000002000 xxtxxyxxxzxxuxuvyxuwzx自然满足在此波动中，质点的位移函数与z无关，且位移在z方向的分量w=0 因为在自由界面上既有入射波带来的位移，又有反射波带来的位移。所以在假定只有反射纵波的时候，有1212;uuu vvv1111sin()UAtf xg y111111cos,sinuUvU22221sin()UAtf xg y222222cos,sinuUvU 代入边界条件可以得到：2211122211112221(2 cos)cos()(2 cos)cos()0sin2cos()sin2

7、cos()0pAtg yAtg yAtg yAtg yV对于y的任意值，为了使第一式得到满足，必须有1211112112(1)();(2)0,ggAAAA 即或上面两式是等价的，说明经过自由界面的反射，位移的位相改变了，但是将其代入第二式却不能使第二式得到满足。于是我们可以知道仅有一反射的纵波，是不能同时满足边界上无正应力和无剪应力的边界条件的。现在我们证明还有一反射的横波，入射纵波、反射纵波、反射横波代入可以同时满足边界上的边界条件。 设反射横波的传播方向与x轴的夹角为2，相应于此波质点的位移函数为：33332sin()UAtf xg y2233cossin;ssfgVV 考虑到横波中质点的

8、位移与传播方向垂直，且假定z方向无运动，因此振动必然发生在oxy平面内，这样，在次波动中质点的位移分量为：332332sin,cosuUvU 入射纵波、反射纵波、反射横波同时存在，在自由界面上质点的位移分量为：123123;uuuu vvvv1111sin()UAtf xg y111111cos,sinuUvU22221sin()UAtf xg y222222cos,sinuUvU 33332sin()UAtf xg y332332sin,cosuUvU代入边界条件可以得到：0000002000 xxtxxyxxxzxxuxuvyxuwzx自然满足整理可以得到：2111222213232121

9、12213232(2 cos)cos()(2 cos)cos()( sin2)cos()0sin2cos()sin2cos()cos2cos()0ppsppsAtg yVAtg yVAtg yVAAtg ytg yVVAtg yV 由此可以得到:1231221122,sinsinsinsinsinppspsgggVVVVV即由此可以得到：，一、二、111221322121232和必须等于零或者等于如果 ，带入边界条件可以得到(A +A )cos2sin-A sinsin2=02(A -A )cossin-A cos2=02221122222sinsin2,2sinsinpsVV即入射波、反射波

10、振幅间的关系结论：平面纵波在自由界面上反射的特点（1）、纵波入射到自由界面上时，将产生两种反射波，即反射纵波和反射横波。（2）、入射波与反射波的传播方向存在关系。当入射波为纵波的时候，反射纵波的射线与自由界面法线的夹角等于入射角，反射横波的射线与自由界面法线的夹角与入射角的夹角，类似与光的折射的关系。（3）如果取反射波的位相改变量为零，入射波和反射波的振幅也存在相应的关系。当入射波的振幅、入射角以及材料的弹性常数已知时，可以求出反射纵波、反射横波的反射角及振幅。反射时波的能量分配关系波的能量与振幅的平方成正比，引入振幅比为纵波的反射系数，为横波的反射系数振幅比反映了反射波相对于入射波振幅的相对

11、强度。212221122122122122211221222121cos2sin-sinsin2=cos2sincossincos2=cossin解得：cos2sin-sinsin2cossincos2cos2sin-sinsin2cossincos2cos2sincos2sin121221122122cossincossincos2sin-sinsin2cossincos2应用线性代数的克莱姆法则22222(1)(1 2 )psVKV12212212122122展开整理可以得到：sinsin2cos 2sinsin2cos 2sincos2sinsin2cos 2sinsin不同材料反射系数

12、随入射角的变化研究入射纵波的应力与反射波的应力之间的关系：入射纵波在传播方向引起的应力为：12(2 )(2 )rtrrUr 111111111cossinsin()sin()cossinsin()pppxyrUAtAtVVrxyrUAtV1(2 )cos()rppArtVV反射纵波的位移方程为：22sin()prUAtV ，取 (2 )(2 )cos()rppUArtrVV反射横波的位移方程为：sin()rUAtV 2cos()rssUArtrVV33结论：入射纵波引起的应力和两个反射波引起的应力都是时间和位置的函数。在自由界面上存在：222222(1)(1 2 )rrsrrppsAAAVAV

13、VKV 12211223121122122sinsin2cos 2-sinsin2cos 2sincos2sinsin2cos 2sinsin入射纵波的应力与反射波应力的关系。二、平面横波在自由界面上的反射二、平面横波在自由界面上的反射横波：质点的运动方向与波的传播方向垂直，但是不限定于那个面内。故考虑任一横波所引起的位移可以是振动方向互相垂直的二波叠加而成。无论质点在什么方向振动，它与波的传播方向垂直，可以分解为一个平行于z轴和垂直于z轴（即在xoy面）的运动。质点平行于z轴振动的横波，我们称为SH波；质点垂直于z轴振动的横波，我们称为SV波。在研究中，通常把横波看作是由两个方向的振动所组成

14、，一个是质点振动在水平平面内的横波分量，称为SV波，一个是质点振动在垂直平面内的横波分量，称为SH波，SV 波SH波（1）质点平行于z轴振动的横波（SH波） ，到达自由界面时应该满足的边界条件。边界条件为：在x=0的平面上0 xxyxz由于振动平行于z轴，故在此波动中质点没有x和y方向的位移，即u=v=0，由于w只是x、y、t的函数，与z没有关系。00 xxyxz自然满足，现在考虑的条件就可以了。由这个唯一的边界条件可以分析得到以下结论： 质点平行于z轴振动的横波（SH波），入射到自由界面时，只产生一个质点是平行于z轴振动的反射横波（SH波）；反射角等于入射角；反射波的振幅与入射波的振幅数值相

15、等，但是位相改变了 ；没有纵波反射。SH 波入射到自由界面SH波11112222cossin;cossin;ssssfgVVfgVV1111sin()wBtf xg y22221sin()wBtf xg y12www1 112221sin()sin()0wB ftg yB ftg yx1212BB 由这个唯一的边界条件可以分析得到以下结论： 质点平行于z轴振动的横波（SH波），入射到自由界面时，只产生一个质点是平行于z轴振动的反射横波（SH波）；反射角等于入射角；反射波的振幅与入射波的振幅数值相等，但是位相改变了 ；没有纵波反射。SH波（2）对于质点垂直于z轴（即在xoy面）振动的横波（sv波

16、），质点在z方向上没有运动。分析问题的方式和入射纵波在自由界面上反射的问题相似边界条件：在自由界面上位移不受任何限制，即应力等于零；即在x=0的平面上，对于y与t的任意值有0 xxyxz0 xz因为质点没有z方向的运动，恒成立由边界条件得到结论： 质点垂直于z轴（即在xoy面）振动的横波（sv1波）入射到自由界面上时，将产生两种反射波，即反射横波sv11和反射纵波sv1 P1。SV波入射到自由界面 入射波与反射波的传播方向存在关系。当入射波为横波的时候，反射横波的射线与自由界面法线的夹角等于入射角，反射纵波的射线与自由界面法线的夹角与入射角有如下的关系。1222121sinsinsinsins

17、inssppsVVVVV由此可以得到：， 如果取反映反射时位相改变的量等于零，则反射横波的振幅和反射纵波的振幅与入射横波的振幅之间的关系为111211321121312和必须等于零或者等于如果 ，带入边界条件可以得到(B +B )sin2 sin-B sincos2=0(B -B )cos2-2B sincos=0用位移场来讨论 取自由界面为xoy平面，y为波前与自由界面的交线，z轴指向介质内部。 波前面垂直于xoz面， 因此 和 都与y无关。(,)xyz u vU rr0rP波入射SV波入射直角坐标上进行投影( )( )yxzuxzawzxvbzx令yyzxxzzyyxuxyzxzvyzxz

18、xwzxyzx质点平行于xoz面的位移由平面纵波P和平面SV波引起；质点垂直于xoz面的位移由平面SH、波引起。P波和平面SV波0uxzwzxv2222222200psvtvt()( )()( )iK ctcf z eiK ctxg z e22222222221010psd fck fdzvd gck gdzv上面公式的通解形式为：22222222(1)(1)12(1)(1)12( )( )ppsscciKziKzvvcciKziKzvvf zAeA eg zBeB eP波和平面SV波22222222(1)(1)12(1)(1)12ppsscciK ct zxiK ct zxvvcciK ct

19、 zxiK ct zxvvAeA eBeB e 222221( )pppcvcvvAKtgcAOpcossin ( )xxxxzd 是入射线和 的夹角，称为入射角，取 沿入射的 波方向P波和平面SV波sincos()()cos11()cos1()coscos222 coscosOBpzxiK ctiK ct ztgxKi Kctxi wt kxpAeAeAeAexwKccvKkKkK 与简谐波比较，它代表沿 方向传播的平面纵波。传播速度相应的波长同理第二项代表沿方向的 波222221ssscvcvvAKtgfAO()()12iK ct ztgfxiK ct ztgfxBeB e边界条件0000

20、002000zztxzyxxxzzzwzwvyzuwzx自然满足( )uxzawzx将代入上式得：22220222220(2 ()0(2)0zzzx zx zxz 求解2222222222112222coscospppssscvcvvcvcvvpsAKtgAOAKtgfAOvfv2222()()12()12122212222222121P,3,cos3cos42(-)(1)0()(1 3)2041 3=41 3psiK ct ztgxiK ct ztgxiK ct ztgfxvvfAeA eAeA A tgBtg fAAtgB tgftg tgftgAAtg tgftgBA当入射只有 波时：取

21、代入边界条件 得：求解：222141 3=41 3tgtgtg tgftg当只有SV波入射时2222222222112222coscospppssscvcvvcvcvvpsAKtgAOAKtgfAOvfv2222()2()()122122212222221211SV,3 ,cos3cos4()(1) 202(-)(1 3)041 3=41 34=psiK ct ztgxiK ct ztgfxiK ct ztgfxvvfAeBeB eBBtg fAtgB B tgfAtgtgftgABtg tgftgtBB当入射只有波时：取代入边界条件 得：求解：22221 341 3g tgftgtg tgf

22、tg当SH波入射时0000=0=-=,=-xzyxzyxzxzuwvzxxyzyxzGGzxr因为令222222222222121122,00 xyzsxyzssik ct c ztgfik ct c ztgfvtGvGztGvGxtGGGGBeGB e 7-2平面波在两种介质分界面上的平面波在两种介质分界面上的反射和透射反射和透射 实际地球内部是成层结构的，我们称为层状介质。地震波在层与层之间的分界面上将会产生反射和透射。 可以证明，当任何一种弹性波到达紧密接触的两种介质的分界面时，一般来讲，都会产生四种波，其中两种波透射到第二种介质中去；另外两种波反射回原来的介质。 紧密接触的两种介质：对

23、于地球内部的情况来说，压缩很大，对于相当微弱的弹性振动，可以认为分界面两边的介质质点是紧密接触的，界面上的介质是连续的。 在两种介质的分界面上的位移和应力都是连续的。 于是可以得到分界面上的边界条件。边界条件：在分界面上有力的边界条件：分界面两边的应力相等；在分界面上有位移的边界条件：分界面两边的位移相等。即：下述四个量应该相等1、正应力2、剪应力3、质点的法向位移4、质点的切向位移上述每个量都由五个位移分量组合而成。其中三个分量来自原来介质的入射波和两个反射波，另外两个来自第二种介质的两个透射波。一、入射波为纵波。yoz面为两种介质a介质和b介质的分界面，x轴为界面的法线，一平面纵波平行于x

24、oy面与x轴夹角为1的方向入射到自由界面。层状介质P1入射纵波P1S1反射横波P12透射纵波P1S2透射横波设入射纵波中质点的位移函数为：1111sin()UAtf xg y1111cossin;papafgVV相应的位移分量为：111111cos,sinuUvU设反射纵波中质点的位移函数为：2222sin()UAtf xg y2222cossin;papafgVV相应的位移分量为：222222cos,sinuUvU 3333sin()UAtf xg y2233cossin;sasafgVV相应的位移分量为：332332sin,cosuUvU设反射横波中质点的位移函数为：设透射纵波中质点的位移

25、函数为：4444sin()UAtf xg y3344cossin;pbpbfgVV相应的位移分量为：443443cos,sinuUvU5555sin()UAtf xg y3355cossin;sbsbfgVV相应的位移分量为：553553sin,cosuUvU设透射横波中质点的位移函数为：在a介质中质点的总位移分量为：123123;aauuuu vvvv在b介质中质点的总位移分量为： 0000abxxabxxuuvv 设入射纵波的各个参数为已知，于是可以由边界条件确定反射波和投射波的各参数。1、在分界面上位移连续，有4545;bbuuu vvv代入可得：11122232343453511122

26、2323434535cossin()cossin()sinsin()cossin()sinsin()sinsin()sinsin()cossin()sinsin()cossin()Atg yAtg yAtg yAtg yAtg yAtg yAtg yAtg yAtg yAtg y1234533122sinsinsinsinsinpapasapbsbgggggVVVVV121324353121324353cossincossin0sincossincos0AAAAAAAAAA0000 xxxxabxyxyxxab2、在分界面上应力连续，有122324353()cos2sin2cos2sin20p

27、asabbpbsbaaAA VAVAVAV21213224353()sin2cos2sin2cos20paasasapapabsbpasbVVAAAVVVVAAVV结论：1、入射纵波到达两种介质的分界面上时，反射两种波，即反射纵波和反射横波；透射两种波，即透射纵波和透射横波。2、入射波、反射波及透射波的传播方向之间存在关系：33122sinsinsinsinsinpapasapbsbVVVVV3、入射波、反射波及透射波的振幅之间存在关系：121211212121221212121212211121121121sincossincossincossincossincossin2cos2sin2c

28、os2sin2cos2sin2pppspppspppspppspsppspppspppssspspspppspppAABBAABBVV VV VAABBVV VVVVAABV 2221211111cos2sin2cos2spsppVBVV 121212122122121212122211121122122121211111sincossincoscossincossinsin2cos2sin2cos2cos2sin2cos2sin2pppsppspssspspppspsspppAVV VV VAVV VVBBVVVVVV1111sincossin2cos2写成矩阵形式：写成矩阵形式：这就是著名

29、的佐普瑞兹这就是著名的佐普瑞兹 Zoepritz方程。方程。11111212pppspppsAPAPSBPBPS设为的反射系数为的反射系数为的反射系数为的反射系数1121ppVV21ppVV 12ppVV12ppVV1pspppsppBBAA根据能量守恒原则有：根据能量守恒原则有：1352141=002bpbapabpbapaapabpbapaAAVVAAVVVAAVV当时，其余四个角都等于零，垂直入射（或法向入射）时波的反射和透射垂直入射（或法向入射）时波的反射和透射由此可推导出用于揭示入射、反射、透射波能量分配关系的佐由此可推导出用于揭示入射、反射、透射波能量分配关系的佐普瑞兹普瑞兹(Zo

30、epritz)(Zoepritz)方程。在该方程中，令入射角方程。在该方程中，令入射角0求解得：求解得： 要把质点的振动分解为平行于z轴和在oxy平面内的两个分量。1、对于质点平行于z轴的振动；即SH波它没有垂直于分界面的运动，因此不产生反射和透射的纵波，即存在如下关系：312112511253sin;sin0sin2()sin20sbsaabVVBBBBBB二、入射波为横波二、入射波为横波2、对于质点垂直于z轴（即在oxy平面内）的振动；即SV波SV1：入射横波SV11：反射横波SV1P1：反射纵波SV1P2：透射纵波SV12：透射横波入射波、反射波即透射波相应的各参数关系如下：33122s

31、insinsinsinsinsapasapbsbVVVVV121324353121324353sincoscossin0cossinsincos0BBBBBBBBBB122314353()sin2cos2cos2sin20sapabbpbsbaaBB VBVB VBV121324353()cos2sin2sin2cos20saasapasbbsbpbVVBBBVVVBBV1341251252=00,0()0=0asabsbasabsbBBBBBVBBV BVVB当时，其余四个角都等于零，无纵波当时无任何的反射波垂直入射（或法向入射）时波的反射和透射垂直入射（或法向入射）时波的反射和透射求解得：

32、求解得： 全反射问题： 当反射波和透射波的速度大于入射波的速度时，就存在入射的临界角，在这种角度之下，反射角和透射角成为90。如果入射角比它还大，关系式就不对。从而出现了类似在光学中所谓的全反射的情况。全反射全反射:;1212 vv当当 到一定程度，但还未到到一定程度，但还未到 时，时， 已增已增大到大到 ，这时透射波在第二种介质中沿界，这时透射波在第二种介质中沿界面面“滑行滑行”，出现，出现“全反射全反射”现象。现象。1。902。901122sin,arcsinccvvvv开始出现开始出现“全反射全反射”时的入射角叫临界角时的入射角叫临界角 ,c对于入射角大于临界角的问题，必须借助复变函数来

33、解决。 可以证明。代替一部分反射或透射的平面波，将产生一种波动，它是随着离开分界面距离的增加而按指数规律减弱的波，这个波并不从边界带走能量，入射波的能量是分给了剩下来的反射波和透射波，但是，这种衰减波的存在，使所产生的其他波的位相有所改变。纵波入射到两种介质的分界面上的透射纵波；SV波入射于两种介质的分界面上时的透射纵波和透射SV波；或者SH波入射于两种介质的分界面上时的透射SH波等一般也会出现全反射的情况。折射波的形成与传播折射波的形成与传播 折射波的形成条件：折射波的形成条件：V2V1 1221sinsin,9090,90SnellVVVV根据定律：若：则一定有个度时，1122sinsin

34、VVVV112sin ()VV波动传播原理波动传播原理 一一、惠更斯原理、惠更斯原理tttttrV t 惠更斯(Huygens)最早（1690 年）提出了由某时刻波前求另一时刻波前的原理，其基本点是：波前上所有点都可以看作广义绕射源，此源发出二次子波，下一时刻的波前即为这些二次子波波前的包络。二、费马原理二、费马原理(Fermat)（最小时间原理）（最小时间原理） 波动总是沿传播时间最小的路径传播，这些路径就是波动总是沿传播时间最小的路径传播，这些路径就是射线射线。 在均匀介质中，射线为直线；在均匀层状介质中，射线满在均匀介质中，射线为直线；在均匀层状介质中，射线满足足斯奈尔定律斯奈尔定律，即

35、波动沿满足斯奈尔定律的路径传播时所用时，即波动沿满足斯奈尔定律的路径传播时所用时间最小。同时，折射波的存在也有力地说明了费马原理。间最小。同时，折射波的存在也有力地说明了费马原理。 三、视速度定理三、视速度定理*VV真速度真速度V V 沿射线方向估计出的波的传播速度。沿射线方向估计出的波的传播速度。出射角出射角 射线与地面法线之间的夹角。射线与地面法线之间的夹角。*V视速度视速度 由于观测方向偏离射线方向，此时由观测数据由于观测方向偏离射线方向，此时由观测数据估计出的波的传播速度将不等于真速度，而称为视速度。估计出的波的传播速度将不等于真速度，而称为视速度。视速度定理：视速度定理： *sin( )cos( )VVVe 按波在传播过程中的传按波在传播过程中的传播路径：直达波，反射波，播路径：直达波，反射波，折射波，透射波。折射波，透射波。 由震源出发向外传播，没有遇到分界面直接到达接收点的波叫直达波。一个纵波入射到反射面时 ，即产生反射纵波和反射横波，也产生透射纵波和透射横波。与入射波类型相同的反射波或透射波称为同类波。改变了类型的反射波或透射波称为转换波。入射角不大，转换波很小，垂直入射不产生转换波。