版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、第九章第九章 平面波的反射与透射平面波的反射与透射 在无限弹性介质中有无旋波(纵波)和等容波(横波)这两种弹性波的传播。 但是实际介质内部,一般都会存在分界面。如地表表面、岩层的分界面。 通常把实际地层称为层状介质。本章研究弹性波在传播的过程中遇到分界面的情况 纵波和横波 地震波到达之处,介质就产生形变。由力学定律知道,任何小的形变都可以分解为两部分:一部分表示胀缩,即变体积而不变形状;另一部分表示畸变,即变形状而不变体积。 形变传播时,两部分的传播速度不同。在震源附近,两部分还未分开,所以波经过处的形变是复杂的。在较远的地方,波阵面就分成两个。胀缩波传播较快,波阵面上的质点位移和传播方向一致
2、,所以叫做纵波,一般用字母P表示。较慢的叫畸变波,质点位移和传播方向垂直,所以叫做横波,一般用字母S表示。在均匀的介质中,波阵面在震源附近是曲面,但在相当距离后就趋近于平面。为简单起见,现只讨论平面波。 1;.2地震地震P波波(纵波纵波)和和S波波(横波横波)运行时弹性岩石运动的形态运行时弹性岩石运动的形态 34纵波:间隔形成压缩带(密集带)和膨胀带(稀疏带),传播方向与振动方向一致,波速 Vp横波:传播方向与振动方向相垂直,波速Vs水平面内分量:称SH波 垂直面内分量:称SV波 59-1平面波在自由界面上的反射平面波在自由界面上的反射一、平面纵波在自由界面上的反射一、平面纵波在自由界面上的反
3、射 设半无限弹性介质的自由表面为yoz面,z轴与图面垂直。假定yoz面的左边为真空,由于没有传播振动的介质,故不会产生透射问题,全部入射波都在界面上被反射。 因为地表面可以近似的看作为自由界面,任意平面波又可以用频谱分析的方法分解为许多简谐平面波的线性叠加,于是我们现在讨论简谐平面波在自由界面上的反射就可以了 。P1入射纵波P1S1反射横波P11反射纵波6设一平面纵波与x轴夹角为1的方向入射到自由界面;设入射纵波中质点的位移函数(即波动方程的解)为:22222pUUVtr1212()()ppUUUf rV tfrV t1111cossinsin()pxyUAtV1111sin()UAtf xg
4、 y改写为1111cossin;ppfgVV7 设入射波为拉伸波,即质点的运动方向与波的前进方向相反,于是可以知道在入射波中质点的位移分量为:111111cos,sinuUvU 假设与自由界面作用后,只有纵波反射,且反射纵波的位移函数为:22221sin()UAtf xg y2222cossin;ppfgVVF2前面的负号原因:反射波相对于x轴而言,是正向传播。1为常数,表示由反射所引起的相位改变。仍假设反射波为拉伸波,在此波动中质点的位移分量为:222222cos,sinuUvU 8问题:确定反射波的各参数问题:确定反射波的各参数方法:由边界条件可以确定方法:由边界条件可以确定边界条件:在自
5、由界面上位移不受任何限制,即应力等于零;即在x=0的平面上,对于y与t的任意值有0 xxyxz由广义胡克定理和几何方程可以得到:0000002000 xxtxxyxxxyxxuxuvyxuwzx自然满足在此波动中,质点的位移函数与z无关,且位移在z方向的分量w=09 因为在自由界面上既有入射波带来的位移,又有反射波带来的位移。所以在假定只有反射纵波的时候,有1212;uuu vvv1111sin()UAtf xg y111111cos,sinuUvU22221sin()UAtf xg y222222cos,sinuUvU 代入边界条件可以得到:2111222211112221(2 cos)co
6、s()(2 cos)cos()0sin2cos()sin2cos()0pAtg yAtg yAtg yAtg yV10对于y的任意值,为了使第一式得到满足,必须有1211112112(1)();(2)0,ggAAAA 即或上面两式是等价的,说明经过自由界面的反射,位移的位相改变了,但是将其代入第二式却不能使第二式得到满足。于是我们可以知道仅有一反射的纵波,是不能同时满足边界上无正应力和无剪应力的边界条件的。现在我们证明还有一反射的横波,入射纵波、反射纵波、反射横波代入可以同时满足边界上的边界条件。11 设反射横波的传播方向与x轴的夹角为2,相应于此波质点的位移函数为:33332sin()UAt
7、f xg y2233cossin;ssfgVV 考虑到横波中质点的位移与传播方向垂直,且假定z方向无运动,因此振动必然发生在oxy平面内,这样,在次波动中质点的位移分量为:332332sin,cosuUvU12 入射纵波、反射纵波、反射横波同时存在,在自由界面上质点的位移分量为:123123;uuuu vvvv1111sin()UAtf xg y111111cos,sinuUvU22221sin()UAtf xg y222222cos,sinuUvU 33332sin()UAtf xg y332332sin,cosuUvU代入边界条件可以得到:0000002000 xxtxxyxxxyxxux
8、uvyxuwzx自然满足13整理可以得到:211122221323212112213232(2 cos)cos()(2 cos)cos()( sin2)cos()0sin2cos()sin2cos()cos2cos()0ppsppsAtg yVAtg yVAtg yVAAtg ytg yVVAtg yV 14由此可以得到:1231221122,sinsinsinsinsinppspsgggVVVVV即由此可以得到:,一、二、111221322121232和必须等于零或者等于如果 ,带入边界条件可以得到(A +A )cos2sin-A sinsin2=02(A -A )cossin-A cos2
9、=02221122222sinsin2,2sinsinpsVV即入射波、反射波振幅间的关系15结论:平面纵波在自由界面上反射的特点(1)、纵波入射到自由界面上时,将产生两种反射波,即反射纵波和反射横波。(2)、入射波与反射波的传播方向存在关系。当入射波为纵波的时候,反射纵波的射线与自由界面法线的夹角等于入射角,反射横波的射线与自由界面法线的夹角与入射角的夹角,类似与光的折射的关系。(3)如果取反射波的位相改变量为零,入射波和反射波的振幅也存在相应的关系。当入射波的振幅、入射角以及材料的弹性常数已知时,可以求出反射纵波、反射横波的反射角及振幅。16反射时波的能量分配关系波的能量与振幅的平方成正比
10、,引入振幅比为纵波的反射系数,为横波的反射系数振幅比反映了反射波相对于入射波振幅的相对强度。21222112212212212221122122212cos2sin-sinsin2=cos2sincossincos2=cossin解得:cos2sin-sinsin2cossincos2cos2sin-sinsin2cossincos2cos2sincos2sin1121221122122cossincossincos2sin-sinsin2cossincos21722222(1)(1 2 )psVKV12212212122122展开整理可以得到:sinsin2 cos 2sinsin2 cos
11、 2sincos2sinsin2 cos 2sinsin18研究入射纵波的应力与反射波的应力之间的关系:入射纵波在传播方向引起的应力为:12(2 )(2 )rtrrUr 111111111cossinsin()sin()cossinsin()pppxyrUAtAtVVrxyrUAtV1(2 )cos()rppArtVV19反射纵波的位移方程为:22sin()prUAtV ,取 (2 )(2 )cos()rppUArtrVV反射横波的位移方程为:sin()rUAtV 2cos()rssUArtrVV3320结论:入射纵波引起的应力和两个反射波引起的应力都是时间和位置的函数。在自由界面上存在:22
12、2222(1)(1 2 )rrsrrppsAAAVAVVKV 12211223121122122sinsin2 cos 2-sinsin2 cos 2sincos2sinsin2 cos 2sinsin入射纵波的应力与反射波应力的关系。21二、平面横波在自由界面上的反射二、平面横波在自由界面上的反射横波:质点的运动方向与波的传播方向垂直,但是不限定于那个面内。故考虑任一横波所引起的位移可以是振动方向互相垂直的二波叠加而成。无论质点在什么方向振动,它与波的传播方向垂直,可以分解为一个平行于z轴和垂直于z轴(即在xoy面)的运动。质点平行于z轴振动的横波,我们称为SH波;质点垂直于z轴振动的横波,
13、我们称为SV波。22在研究中,通常把横波看作是由两个方向的振动所组成,一个是质点振动在水平平面内的横波分量,称为SH波,一个是质点振动在垂直平面内的横波分量,称为SV波,SH 波SV波23(1)质点平行于z轴振动的横波(SH波) ,到达自由界面时应该满足的边界条件。边界条件为:在x=0的平面上0 xxyxz由于振动平行于z轴,故在此波动中质点没有x和y方向的位移,即u=v=0,由于w只是x、y、t的函数,与z没有关系。00 xxyxz自然满足,现在考虑的条件就可以了。由这个唯一的边界条件可以分析得到以下结论: 质点平行于z轴振动的横波(SH波),入射到自由界面时,只产生一个质点是平行于z轴振动
14、的反射横波(SH波);反射角等于入射角;反射波的振幅与入射波的振幅数值相等,但是位相改变了 ;没有纵波反射。24(2)对于质点垂直于z轴(即在xoy面)振动的横波(sv波),质点在z方向上没有运动。分析问题的方式和入射纵波在自由界面上反射的问题相似边界条件:在自由界面上位移不受任何限制,即应力等于零;即在x=0的平面上,对于y与t的任意值有0 xxyxz0 xz因为质点没有z方向的运动,恒成立由边界条件得到结论: 质点垂直于z轴(即在xoy面)振动的横波(sv1波)入射到自由界面上时,将产生两种反射波,即反射横波sv11和反射纵波sv1 P1。25 入射波与反射波的传播方向存在关系。当入射波为
15、横波的时候,反射横波的射线与自由界面法线的夹角等于入射角,反射纵波的射线与自由界面法线的夹角与入射角有如下的关系。1222121sinsinsinsinsinssppsVVVVV由此可以得到:, 如果取反映反射时位相改变的量等于零,则反射横波的振幅和反射纵波的振幅与入射横波的振幅之间的关系为111211321121312和必须等于零或者等于如果 ,带入边界条件可以得到(B +B )sin2 sin-B sincos2=0(B -B )cos2-2B sincos=0269-2平面波在两种介质分界面上的反射和透射平面波在两种介质分界面上的反射和透射 实际地球内部是成层结构的,我们称为层状介质。地
16、震波在层与层之间的分界面上将会产生反射和透射。 可以证明,当任何一种弹性波到达紧密接触的两种介质的分界面时,一般来讲,都会产生四种波,其中两种波透射到第二种介质中去;另外两种波反射回原来的介质。 紧密接触的两种介质:对于地球内部的情况来说,压缩很大,对于相当微弱的弹性振动,可以认为分界面两边的介质质点是紧密接触的,界面上的介质是连续的。 在两种介质的分界面上的位移和应力都是连续的。 于是可以得到分界面上的边界条件。27边界条件:在分界面上有力的边界条件:分界面两边的应力相等;在分界面上有位移的边界条件:分界面两边的位移相等。即:下述四个量应该相等1、正应力2、剪应力3、质点的法向位移4、质点的
17、切向位移上述每个量都由五个位移分量组合而成。其中三个分量来自原来介质的入射波和两个反射波,另外两个来自第二种介质的两个透射波。一、入射波为纵波。yoz面为两种介质a介质和b介质的分界面,x轴为界面的法线,一平面纵波平行于xoy面与x轴夹角为1的方向入射到自由界面。28P1入射纵波P1S1反射横波P12透射纵波P1S2透射横波设入射纵波中质点的位移函数为:1111sin()UAtf xg y1111cossin;papafgVV相应的位移分量为:111111cos,sinuUvU29设反射纵波中质点的位移函数为:2222sin()UAtf xg y2222cossin;papafgVV相应的位移
18、分量为:222222cos,sinuUvU 3333sin()UAtf xg y2233cossin;sasafgVV相应的位移分量为:332332sin,cosuUvU设反射横波中质点的位移函数为:30设透射纵波中质点的位移函数为:4444sin()UAtf xg y3344cossin;pbpbfgVV相应的位移分量为:443443cos,sinuUvU5555sin()UAtf xg y3355cossin;sbsbfgVV相应的位移分量为:553553sin,cosuUvU设透射横波中质点的位移函数为:31在a介质中质点的总位移分量为:123123;aauuuu vvvv在b介质中质点
19、的总位移分量为: 0000abxxabxxuuvv 设入射纵波的各个参数为已知,于是可以由边界条件确定反射波和投射波的各参数。1、在分界面上位移连续,有4545;bbuuu vvv32代入可得:111222323434535111222323434535cossin()cossin()sinsin()cossin()sinsin()sinsin()sinsin()cossin()sinsin()cossin()Atg yAtg yAtg yAtg yAtg yAtg yAtg yAtg yAtg yAtg y1234533122sinsinsinsinsinpapasapbsbgggggVVV
20、VV121324353121324353cossincossin0sincossincos0AAAAAAAAAA330000 xxxxabxyxyxxab2、在分界面上应力连续,有122324353()cos2sin2cos2sin20pasabbpbsbaaAA VAVAVAV21213224353()sin2cos2sin2cos20paasasapapabsbpasbVVAAAVVVVAAVV34结论:1、入射纵波到达两种介质的分界面上时,反射两种波,即反射纵波和反射横波;透射两种波,即透射纵波和透射横波。2、入射波、反射波及透射波的传播方向之间存在关系:33122sinsinsinsi
21、nsinpapasapbsbVVVVV3、入射波、反射波及透射波的振幅之间存在关系:35121211212121221212121212211121121121sincossincossincossincossincossin2cos2sin2cos2sin2cos2sin2pppspppspppspppspsppspppspppssspspspppspppAABBAABBVV VV VAABBVV VVVVAABV 2221211111cos2sin2cos2spsppVBVV 121212122122121212122211121122122121211111sincossincoscos
22、sincossinsin2cos2sin2cos2cos2sin2cos2sin2pppsppspssspspppspsspppAVV VV VAVV VVBBVVVVVV1111sincossin2cos2写成矩阵形式:写成矩阵形式:这就是著名的佐普瑞兹这就是著名的佐普瑞兹 Zoepritz方程。方程。ppApsAppBpsB36111222pspsVVVV、ppApsA ppBpsB 0pspsBA371121ppVV21ppVV 12ppVV12ppVV1pspppsppBBAA根据能量守恒原则有:根据能量守恒原则有:3839221 1221 11 1221 12RRITTIAVVKAV
23、VAVKAVV垂直入射(或法向入射)时波的反射和透射垂直入射(或法向入射)时波的反射和透射假设:A.各向同性弹性介质; B.波通过分界面时满足两个连续性条件: 1. 弹性位移的法向分量和切向分量连续; 2. 应力的法向分量与切向分量连续。依据:A.斯奈尔定律; B.地震波可分解为无限个简谐振动由此可推导出用于揭示入射、反射、透射波能量分配关系的佐普瑞兹由此可推导出用于揭示入射、反射、透射波能量分配关系的佐普瑞兹(Zoepritz)(Zoepritz)方方程。在该方程中,令入射角程。在该方程中,令入射角0求解得:求解得: 401IRTTRAAAKKTKRK41TK42要把质点的振动分解为平行于z
24、轴和在oxy平面内的两个分量。1、对于质点平行于z轴的振动;即SH波它没有垂直于分界面的运动,因此不产生反射和透射的纵波,即存在如下关系:312112511253sin;sin0sin2()sin20sbsaabVVBBBBBB二、入射波为横波二、入射波为横波432、对于质点垂直于z轴(即在oxy平面内)的振动;即SV波SV1:入射横波SV11:反射横波SV1P1:反射纵波SV1P2:透射纵波SV12:透射横波44入射波、反射波即透射波相应的各参数关系如下:33122sinsinsinsinsinsapasapbsbVVVVV121324353121324353sincoscossin0cos
25、sinsincos0BBBBBBBBBB122314353()sin2cos2cos2sin20sapabbpbsbaaBB VBVB VBV121324353()cos2sin2sin2cos20saasapasbbsbpbVVBBBVVVBBV45全反射问题: 当反射波和透射波的速度大于入射波的速度时,就存在入射的临界角,在这种角度之下,反射角和透射角成为90。如果入射角比它还大,关系式就不对。从而出现了类似在光学中所谓的全反射的情况。全反射全反射:;1212 vv46当当 到一定程度,但还未到到一定程度,但还未到 时,时, 已增大到已增大到 ,这时透射波在第,这时透射波在第二种介质中沿界面二种介质中沿界面“滑行滑行”,出现,出现“全反射全反射”现象。现象。1。902。9012sincvv开始出现开始出现“全反射全反射”时的入射角叫临界角时的入射角叫临界角 ,c47对于入射角大于临界角的问题,必须借助复变函数来解决。 可以证明。代替一部分反射或透射的平面波,将产生一种波动,它是随着离开分界面距离的增加而按指数规律减弱的波,这个波并不从边界带走能量,入射波的能量是分给了剩下来的反射波和透射波,但是,这种衰减波的存在,使所产生的其他波的位相有所改变。横波SV1入射于自由界面上时,临界角的求法:
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- WPS 办公应用-教学大纲、授课计划
- 2024年汽车热交换器项目资金申请报告代可行性研究报告
- 有关追梦演讲稿(17篇)
- 文明礼仪伴我行演讲稿400(34篇)
- 学校表彰大会校长致辞
- 河西走廊观后感600字范文(6篇)
- 珍惜粮食学生个人倡议书
- 理疗师劳务合同范本
- 疫情期间幼儿工作总结5篇
- 新教材高考地理二轮专题复习单元综合提升练3地球上的水含答案
- 安全总监安全职责
- 云南白族课件
- 消防应急预案组织结构图
- 油站使用说明书
- 小学班主任工作经验交流ppt
- 如何识别真假币(共34张PPT)
- 2023届高考数学复习微难点7 三角函数中ω的范围问题(共11张PPT)
- 计算机科学与技术本科专业自评报告(共64页)
- 工程建设情况汇报PPT课件
- GB∕T 39116-2020 智能制造能力成熟度模型
- 数学-九宫数独100题(附答案)
评论
0/150
提交评论