1、姓名： 學(xué)號(hào)： 專(zhuān)業(yè)：地球物理勘察技術(shù) 2012級(jí)一、實(shí)驗(yàn)?zāi)康? .利用密度、上下界面的縱橫波速度通過(guò)斯奈爾定律求出該界面的00 90°入射角的反射角度和透射角度。sin sin 1 sin 1 sin 2 sin 2 pVp1Vp1 Vs1 Vp2 Vs22 .利用Zoeppritz方程繪制反射系數(shù)和透射系數(shù)曲線。二、實(shí)驗(yàn)步驟1、模型：3Vp1=3300 m/s V s1=1585 m/s p 1 = g/cm3Vp1=3100 m/s V s1=1989 m/s p 1 = g/cm計(jì)算從090度的反射透射系數(shù)曲線 如圖1-1-4三、實(shí)驗(yàn)結(jié)果P波入射的反射透射角度：圖1P波入射的

2、反射透射系數(shù)曲線：圖2SV波入射的反射透射角度正弦值：圖3SV波入射的反射透射系數(shù)曲線：圖4SH波入射的反射透射角度正弦值：圖5SH波入射的反射透射系數(shù)曲線：四、實(shí)驗(yàn)分析因?yàn)楦鞑ǚ瓷浜屯干洳ㄕ穹禂?shù)與其能量成正比，由此可以看出其能量的變化以及入射波能量的分配。當(dāng)上層介質(zhì)為密介質(zhì)，下層介質(zhì)為疏介質(zhì)時(shí)：由圖2知，P波從上層密介質(zhì)入射到界面時(shí)：隨著入射角增大，P波反射系數(shù)、P波透射系數(shù)在減小，即隨著入射角的增大，反射 P波、透射P波的能量 在減小，而S波反射系數(shù)、S波透射系數(shù)在增大，即隨著入射角的增大，反射 S 波、透射S 波的能量在增大。說(shuō)明當(dāng)入射角發(fā)生變化時(shí)，入射波的能量分配在改變。由圖4可知，

3、SV波從上層密介質(zhì)入射到界面時(shí)：入射角小于29度時(shí)，隨著入射角增大，SV波反射系數(shù)在減小，P波反射系 數(shù)、SV波透射系數(shù)、P波透射系數(shù)在增大，SV波透射系數(shù)最大，即反射 SV波 能量在減小，反射P波、透射SV波、透射P波能量在增大，透射SV波的能量 最大；入射角大于等于29度且小于31度時(shí)，SV波反射系數(shù)、SV波透射系數(shù)減 小，P波反射系數(shù)、P波透射系數(shù)在增大，SV波透射系數(shù)最大，即隨著入射角 的增大，反射SV波、透射SV波能量在減小，反射P波、透射P波能量在增 大，透射SV波能量最大；入射角大于等于31度且小于53度時(shí)，SV波透射系數(shù)增大，SV波反射系 數(shù)、P波反射系數(shù)、P波透射系數(shù)先減小后

4、增大，SV波透射系數(shù)仍最大，即隨 著入射角的增大，透射SV波能量增大，反射SV波、反射P波、透射P波能量 先減小后增大，透射SV波能量最大；入射角大于53度時(shí)，SV波反射系數(shù)為1, SV波透射系數(shù)、P波反射系 數(shù)、P波透射系數(shù)減小，即隨著入射角的增大，反射 SV波能量不變，透射SV 波、反射P 波、透射P 波能量減小。由圖6可知，SH波從上層密介質(zhì)入射到界面時(shí)：不產(chǎn)生轉(zhuǎn)換波，入射角小 于臨界角時(shí)，SH波的反射系數(shù)、透射系數(shù)均隨著入射角的增大而增大，即反 射、透射SH波能量增大；入射角大于等于臨界角時(shí)，隨著入射角的增大，SH波的反射系數(shù)幾乎不變，透射系數(shù)減小，即反射 SH波的能量減小，透射SH波

5、 的能量幾乎不變。五、附：源程序代碼P 波入射時(shí)，程序：#include<> #include<>#include"" #define PIvoid main()FILE *fp1,*fp2;int i,n91;doubleipp,x1,x2,y1,y2,pr,sr,pt,st,a195,a295,b195,b295,vp1=3300, vp2=3100,vs1=1585,vs2=1989,den1=,den2=,k=den2/den1;static double a44=,b4=;fp1=fopen("","w&quo

6、t;);fp2=fopen("P 波入射反射透射系數(shù).csv","w");for(i=0;i<=90;i+)ni=i;ipp=i*PI/180;x1=ipp;a1i=x1*180/PI;x2=sin(ipp)*vs1/vp1;a2i=asin(x2)*180/PI;y1=sin(ipp)*vp2/vp1;b1i=asin(y1)*180/PI;y2=sin(ipp)*vs2/vp1;b2i=asin(y2)*180/PI;fprintf(fp1,"%d,%f,%f,%f,%fn",ni,a1i,a2i,b1i,b2i);sv&q

7、uot;,"w");ipp1=asin(vs1/vp1);ipp2=asin(vs1/vs2);ipp3=asin(vs1/vp2);sv","w");ipp=asin(vs1/vs2);/* 臨界角for(i=0;i<=90;i+)ni=i;in=i*PI/180;x=sin(in);y=sin(in)*vs2/vs1;/*snell 定律ri=x;ti=y;fprintf(fp1,"%d,%f,%fn",ni,ri,ti);/*輸出SH波入射反射透射角度正弦值if(in<ipp)st=asin(y);p=vs1*cos(in);q=vs2*cos(st);rsh=(p-k*q)/(p+q);tsh=(2*p+(1-k)*q)/(p+q);fprintf(fp2,"%d,%f,%fn",ni,rsh,tsh);/*輸出 SH波入射反射透射系數(shù)elsep=vs1*cos(in);q=vs2*sqrt(y*y-1);rsr=(p*p-k*q*q)/(p*p+q*q);rsi=-(1+k)*p*q/(p*p+q*q);tsr=(2*p*p+(1-k)*q*q)/(p*p+q*q);tsi=-(1+k)*p*q/(p*p+q*q);rsh=sqrt(rsr*