叠前地震记录的相移波动方程正演模拟实验

本科生实验报告实验课程油气勘探新方法学院名称地球物理学院专业名称勘查技术与工程（石油物探）学生姓名学生学号指导教师熊高君实验地点5417实验成绩2015年12月1成都理工大学《油气勘探新方法》实验报告实验时间2015年12月开课单位地球物理学院指导教师熊高君实验题目：叠前地震记录的相移波动方程正演模拟实验姓名学号班级专业勘测技术与工程（石油物探）院（系）地球物理学院地球探测与信息技术系单项成绩内容理解写作结构程序设计模型设计计算结果结果分析总成绩2实验报告一、实验题目叠前地震记录的相移波动方程正演模拟实验二、实验目的掌握各向同性介质任意构造、水平层状速度结构地质模型的叠前地震记录相移波动方程正演模拟基本理论、实现方法与程序编制，由正演记录初步分析地震信号的分辨率。三、原理公式1、地震波传播的波动方程设（x,z）为空间坐标，t为时间，地震波传播速度为v(x,z)，则二维介质中任意位置、任意时刻的地震波场为p(z,x,t)：压缩波——纵波。则二维各向同性均匀介质中地震波传播遵循的声波方程：2、傅里叶变换的微分性质p(t)与其傅里叶变换的P()的关系：{∫正傅里叶变换∫逆傅里叶变换则有时间微分性质：3{∫一阶微分∫二阶微分ω为频率，⁄，T为周期。同理有空间微分性质：{∫一阶微分∫二阶微分k为波数，⁄，为波长。3、地震波传播的相移外推公式令速度v不随x变化，只随z变化，则利用傅里叶变换微分性质（3）和（4）式，把波动方程(1)式变换到频率‐波数域，得：（）或：()令：则(6)式的解为：包括上行波和下行波两项，正演模拟取上行波：4若和间隔为，速度v(z)为在此间隔内不随Z变的常数，(8)式实现波场从到的延拓，即：()（）在深度开始向上延拓到，若延拓深度为零，即：，则()对于任意深度到的延拓，可得正演模拟中地震波的传播方程（延拓公式）()()4、初始条件和边界条件按照爆炸界面理论，反射界面震源在t=0时刻同时起爆，此时刻的波场就是震源。根据不同情况，可直接使用反射系数脉冲或子波作震源。如果直接使用反射系数作震源脉冲，则初始条件可表示为：{其他对时间t和空间x做二维傅立叶变换，则得频率-波数域的初始波场。边界条件：{其他其他其他其他参数都是在范围内定义的。5、边界处理（1）边界反射问题把实际无穷空间区域中求解波场的问题化为有穷区域求解时，左右两边使用零边界条件。物理上假设探区距与两个端点很远，在两个端点上收到5的反射波很弱。但是，上述条件在实际中不能成立，造成零边界条件反而成为绝对阻止波通过的强反射面。在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。（2）边界强反射的处理镶边法、削波法、吸收边界都能有效消除边界强反射。削波法就是在波场延拓过程中，每延拓一次，在其两侧均匀衰减到零，从而消除边界强反射的影响。假设横向总长度为NX，以两边Lx道吸波为例，有以下吸波公式：{√其它6、定位原理某一震源在某一固定检波点产生的记录，就等于该固定检波点主动到所有可能反射点处接收的由这一震源激发的反射波的总和。形象地称为检波点下延记录原理。7、反射记录的形成1）待激发的二次源的描述：介质的离散化把整个介质空间离散化，每个离散点都可以看成一个绕射点，某深度的反射界面与这一深度线的交点为这一反射界面的分布点。每个离散点就是待激发的二次源2）反射界面的分布某一深度层的反射界面分布在这一深度的离散点上，即反射界面与这一深度线的交点为这一反射界面的分布点。3）反射界面二次源的形成：绕射波的产生震源从炸点开始向下传播，遍历介质的每一个离散点绕射点。每遇到一个绕射点，都会激发该点的绕射波。反射界面与某深度离散点的交叉点的反射系数6不为零，产生的绕射波也不为零，而其他没有反射界面的地方，反反射系数为零，是均匀介质，产生的绕射波为零。4）反射记录的形成把检波器的检波机制看成一个数学算子，或检波因子，可用脉冲或雷克子波。在延拓震源的同时，把检波器即检波因子也做完全相同的向下延拓，同样遍历每个离散点。震源和检波因子每延拓一个深度步长，震源、检波因子及反射系数三者做相关运算，则检波点就记录到了该深度每个绕射点的绕射波。所有二次源地震波都传播到测线某检波点，叠加而成该点的记录——炮集中的某一记录道。测线排列多少道检波器，则都可记录道所在位置的地震模拟记录。8、数字化根据数字信号处理的采样定理，把连续的信号变为计算机能处理的数字信号，使相移法正演模拟得以实现。频域抽样定理：一个频谱受限信号，如果时间只占据的范围，若在频域以不大于频率间隔对信号的频谱采样，则抽样到的离散信号可以唯一表示原信号。时域抽样定理：一个时间受限信号，如果频谱只占据的范围，则信号可以用等间隔的抽样值唯一表示出来，而时间抽样间隔必须不大于，，。四、实验内容削波法相位移正演模拟（1）点绕射构造和水平层状速度模型（参数如图1所示）的正演数值模拟；1）削波的正演；2）无削波的正演；（2）计算中点和两个边界的信号位置，分析实验结果的正确性；（3）做同样模型的褶积模型数值模拟，对比分析两者的异同；7（4）改变绕射点位置、速度，再做正演模拟；图1点绕射构造与水平速度模型五、方法路线1、参数初始化；2、二次源的形成：1）震源波场定义：Scr(x,z,t){其它（）2）界面待激发爆炸源定义：Rsht(x,z,t){其它（）3）形成削波数据：Absb(x)在所有深度都用同样的削波数据：{√其它（）84）从测线深度开始向下做波场延拓在所有深度都用同样的削波数据①震源波场做削波处理()()（）②震源做x和t的傅里叶变换到频率-波数域：()→()（）③计算相移延拓算子：()(√())（）④震源波场延拓计算()()()（）⑤延拓后的震源波场变换到时间-空间域：()→()（）⑥激发当前深度的二次源：震源波场与待激发波场相关()∫()()（）⑦将二次源变换到频率域存储：()→()⑧从目前深度开始，重①到⑦步，震源波场继续向下延拓,直至介质最深处，生成介质所有点的绕射波--整个介质的二次源形成。3、正演记录产生的数学过程1）定义频率域延拓波场:；2）调入削波数据：；3）在频率域延拓波场：9①每给一个频率，延拓波场初始化：()；②从最深处开始，波场向上延拓：a、取出当前频率和深度的二次源：；b、形成新延拓波场：()()()（）c、新延拓波场做削波处理()()（）d、新延拓波场变换到波数域：()→()（）e、计算相移延拓算子：()(√())（）f、波场延拓计算()()()（）g、波场变换到空间域：()→()（）h、重复a到g的计算，直至波场延拓到测线深度，ZJ=Z1，得到测线上的频率域波场：()，③存储测线上的波场:();④重复①到③的计算，直到完成所有频率的循环；4)把波场变换到时间域，得一炮的叠前正演模拟结果：10→（）4、频率-空间域波场对频率做反傅里叶变换，得时间-空间波场；5、使用Fimage软件显示所得结果。六、实验结果1、结果显示及分析（点绕射结构）图2削波（左）正演模拟结果与未削波正演模拟结果（右）由图2分析可知：削波前，由于在两个端点上收到的反射波很弱，造成零边界条件反而成为绝对阻止波通过的强反射面，在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。削波后，边界假反射的影响消除了，且曲线变得更光滑了。a、速度V1=5000、V2=5500，深度h=2000时，改变绕射点x的位置；11图3x=Nx/4-1时，反射系数（左）与削波正演模拟结果（右）图4x=Nx/2-1时，反射系数（左）与削波正演模拟结果（右）由图3与图4对比分析可知：当速度，深度一定时，改变绕射点的位置，曲线随绕射点位置的改变而左右移动，且移动趋势与绕射点位置的移动一致。b、绕射点x=Nx/2-1，深度h=2000，速度V2=5500时，改变速度V1；12图5V1=4000时，速度模型（左）与削波正演模拟结果（右）图6V1=6000时，速度模型（左）与削波正演模拟结果（右）由图5与图6对比分析可知：当绕射点位置固定，深度一定时，随着介质速度增大，曲线变得越平缓。这是由于速度增大时，波的传播时间减小所致。c、速度V1=5000、V2=5500，绕射点x=Nx/2-1，改变深度h；V1=4000V2=5500V2=5500V1=600013图7h=1000时，反射系数（左）与削波正演模拟结果（右）图8h=3000时，反射系数（左）与削波正演模拟结果（右）由图7与图8对比分析可知：当介质速度一定，绕射点位置不变时，随着绕射点深度的增加，曲线变得越来越平滑。显然是由于深度变大，波的传播时间变小所致。2、结果分析（水平层速度模型）14图9未削波（左）正演模拟结果与削波正演模拟结果（右）由图9分析可知：削波前，由于在两个端点上收到的反射波很弱，造成零边界条件反而成为绝对阻止波通过的强反射面，在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。削波后，边界假反射的影响消除了，且曲线变得更光滑了。a、深度h=2000，速度V2=5500时，改变速度V1；图10V1=4000时，速度模型（左）与削波正演模拟结果（右）V1=4000V2=550015图11V1=6000时，速度模型（左）与削波正演模拟结果（右）由图10与图11对比分析可知：水平层深度不变时，随着速度的增加，结果显示的水平层向上移动。这是由于速度增大时，波的传播时间减小的。b、速度V1=5000、V2=5500，改变深度h；图12h=1000时，反射系数（左）与削波正演模拟结果（右）V2=5500V1=600016图13h=3000时，反射系数（左）与削波正演模拟结果（右）由图12与图13对比分析可知：当波的传播速度不改变时，随着水平层深度的增加，结果显示的水平层向下移动。这是由于深度增加，波的传播时间增大。七、讨论建议1、实验收获通过本次实验，对叠后相移波动方程正演模拟的过程有了更深入的理解，以及对其公式和意义都有了进一步的认识。编程能力又有了一定提高。2、存在问题图14x=3*Nx/4-1时，反射系数（左）与削波正演模拟结果（右）17如图14，当速度，深度一定时，绕射点位置x=3*Nx/4-1，其反射系数不能用Fimage正确显示，由显示的削波正演模拟结果可知，使用削波函数并未起到削波的作用。3、心得体会本次实验，并不容易。在实验中，必须要理解理论知识，才能够正确的编写程序。在实验过程中，遇到许多问题，通过和同学讨论，向老师求解，最终才得出实验结果。附程序代码：//Include#includestdlib.h#includeconio.h#includemath.h#includestdio.h#includestring.h//Parametersset#defineNx128//TraceNumber#defineNt256//RecordNumber#defineNz100//DepthNumber#defineLabs15//LengthOfBoundaryAbsorbing#defineDx20//TraceInterval#defineDt0.004//RecordInterval#defineDz20//DepthInterval#definePai3.1415926//function:Judgethepowerof2intPo2Judge(intN){intk=0;longLn=0;for(k=0;N-Ln0;k++){Ln=(long)pow(2,k);}Ln=(long)pow(2,k-1);if(fabs(Ln-N)=1)return(0);return(1);}18//function:FormAbsorbBounderyintAbsorb(){FILE*fp_Abs;intIx;doubleAbs[Nx];if((fp_Abs=fopen(Absb.dat,wb))==NULL)printf(ConnotopenfileAbsb);for(Ix=0;IxNx;Ix++)Abs[Ix]=1.;for(Ix=0;Ix=Labs-1;Ix++