射线追 踪及动校正编程报告_中国石油大学.docVIP

射线追踪及动校正编程报告 原理 报告中使用的地震子波为瑞克子波。射线追踪使用的二分法的方式，确定了极大值和极小值后，利用中间值进行试射逼近，然后以中间值为一个极值再次进行取中间值进行试射，循环多次，当射线得到的点与我们理想点之间的差距在给定的误差范围内时，记录此时射线经过各层的坐标和角度，对记录的坐标进行连接，得到射线的追踪路径。射线在各层之间的透射满足斯奈尔定理（公式1）。 Sin(a)/v1=sin(b)/v2 （1） 其中，a是入射角，b是透射角，v1为入射层的层速度，v2为透射层的层速度。 在反射系数的计算中没有考虑每层的透射损失只考虑了最终层位的反射，反射系数使用的计算方法是做不离子方程的简化形式Bortfeld近似式（公式2）。 （2） 其中Rpp表示纵波反射系数。vp1为入射层的纵波速度，vp2为透射层的纵波速度。θ1为入射角，θ2为透射角。ρ1为入射层密度，ρ2为折射层密度。Vs1为入射层横波速度，Vs2为折射层的横波速度。在本文中由于在临界角以外只有入射角，此时的反射系数可能会存在误差。在第六层与第七层的反射系数过小造成在最后的合成单炮记录上几乎肉眼不能看到。 动校正的方法主要使用了常规的双曲动校正和无拉伸动校正，常规双曲动校正的原理是公式（3），在炮检距较大时动校正会有明显的动校拉伸，它的拉伸系数为公式（4）。 （3） 是自激自收的时间，x为炮检距，v为动校正的速度。在实际中动校正的速度应该是变化的，但是为了计算的方便一般是设定成订制，这就导致了动校正在大偏移距时发生动校误差。 （4） △t为反射点的时间与校正到点的时间的差，t为校正的点的时间。 而无拉伸动校正的原理是在原始地震剖面上实现反射层位的部分实现整体搬离。本文中由于动校速度为定值，在大偏移距时存在比较明显的误差。在反射见面附近的动校速度为给定的动校速度，在该区域外运用的速度为插值得到的动校速度进行常规动校，该动校的好坏有直接关系的是双曲线反射界面的追踪的准确程度。 流程图 图1合成地震记录的流程图 图2 常规动校正的流程图  图3 无拉伸动校的流程图 该动校的好坏与中心线的确定有直接关系。 图4 追踪的射线路径  图5时间与炮检距的变化 图4中在起始位置的由大到小为层位由小到大。 图6 反射系数 由于在临界角以后求得的反射系数为复数，在图中取的模长，反射系数第一层界面为深蓝线，第二层反射界面为红线，第三层反射系数为黄线，第四层的反射系数为紫线，第五层的反射系数为绿线，第六层的反射系数为浅蓝线。 图7反射系数在时间域 图8 单炮记录 图9 均方根速度的动校正 图10均方根速度存在校正过量调整 图11 切除拉伸系数30% 图12基于图9的无拉伸动校正 程序 主程序部分 clear all;%清除内存 close all;clc;%，关闭所有窗口 v=[2000;2500;2900;4000;3600;4000;4200];%纵波速度 vs=[700;1050;1350;2000;1900;2200;2400];%横波速度 mi=[2.0;2.1;2.2;2.4;2.3;2.5;2.6];%密度 h=[500;400;200;400;300;200];%各层厚度 nc=6; %层数目 nd=101; %接收道数目 os=0; %最小炮检距； dj=40; %道间距 offset=zeros(1,nd);%偏移距的计算 for i=1:nd; offset(i)=i*dj; end tnmo=zeros(1,nc);%各层垂直入射时间 tnmo(1)=0.5; for i=2:nc; tnmo(i)=tnmo(i-1)+2*h(i)/v(i); end % vnm=zeros(1,nc);%均方根速度 vnmo=zeros(1,nc);%动校正速度 mmm=0; nnn=0; for i=1:nc; for j=1:i; mmm=mmm+tnmo(j)*(v(j)^2); nnn=nnn+tnmo(j); end vnmo(i)=sqrt(mmm/nnn); end vnmo(1)=2040;%由于均方根速度把存在明显地校正过量对速度调整 vnmo(2)=2265; vnmo(3)=2400; vnmo(4)=2760; vnmo(5)=2976; vnmo(6)=3143