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

本科生实验报告实验课程 油气勘探新方法 学院名称 地球物理学院 专业名称 勘查技术与工程（石油物探） 学生姓名 学生学号 指导教师 熊高君 实验地点 5417 实验成绩 2015年12月 成都理工大学《油气勘探新方法》实验报告实验时间2015年 12月开课单位地球物理学院指导教师熊高君实验题目：叠前地震记录的相移波动方程正演模拟实验姓名学号班级专业勘测技术与工程（石油物探）院（系）地球物理学院 地球探测与信息技术系单项成绩内容理解写作结构程序设计模型设计计算结果结果分析总成绩实验报告实验题目叠前地震记录的相移波动方程正演模拟实验二、实验目的掌握各向同性介质任意构造、水平层状速度结构地质模型的叠前地震记录相移波动方程正演模拟基本理论、实现方法与程序编制，由正演记录初步分析地震信号的分辨率。三、原理公式1、地震波传播的波动方程设（x,z）为空间坐标，t为时间，地震波传播速度为v(x,z)，则二维介质中任意位置、任意时刻的地震波场为p(z,x,t)：压缩波——纵波。则二维各向同性均匀介质中地震波传播遵循的声波方程：2、傅里叶变换的微分性质p(t)与其傅里叶变换的P()的关系：则有时间微分性质：ω为频率，，T为周期。同理有空间微分性质：k为波数，，为波长。3、地震波传播的相移外推公式令速度v不随x变化，只随z变化，则利用傅里叶变换微分性质（3）和（4）式，把波动方程(1)式变换到频率‐波数域，得：或：令：则(6)式的解为： 包括上行波和下行波两项，正演模拟取上行波： 若和间隔为，速度v(z)为在此间隔内不随Z变的常数，(8)式实现波场从到的延拓，即： 在深度开始向上延拓到，若延拓深度为零，即：，则 对于任意深度到的延拓，可得正演模拟中地震波的传播方程（延拓公式） 4、初始条件和边界条件按照爆炸界面理论，反射界面震源在 t=0 时刻同时起爆，此时刻的波场就是震源。根据不同情况，可直接使用反射系数脉冲或子波作震源。如果直接使用反射系数作震源脉冲，则初始条件可表示为： 对时间t和空间x做二维傅立叶变换，则得频率-波数域的初始波场。边界条件： 其他参数都是在范围内定义的。5、边界处理（1）边界反射问题把实际无穷空间区域中求解波场的问题化为有穷区域求解时，左右两边使用零边界条件。物理上假设探区距与两个端点很远，在两个端点上收到的反射波很弱。但是，上述条件在实际中不能成立，造成零边界条件反而成为绝对阻止波通过的强反射面。在正演模拟的剖面上出现了边界假反射干涉正常界面的反射。（2）边界强反射的处理镶边法、削波法、吸收边界都能有效消除边界强反射。削波法就是在波场延拓过程中，每延拓一次，在其两侧均匀衰减到零，从而消除边界强反射的影响。假设横向总长度为NX，以两边Lx道吸波为例，有以下吸波公式：6、定位原理某一震源在某一固定检波点产生的记录， 就等于该固定检波点主动到所有可能反射点处接收的由这一震源激发的反射波的总和。形象地称为检波点下延记录原理。7、反射记录的形成 1）待激发的二次源的描述：介质的离散化把整个介质空间离散化，每个离散点都可以看成一个绕射点，某深度的反射界面与这一深度线的交点为这一反射界面的分布点。 每个离散点就是待激发的二次源 2）反射界面的分布某一深度层的反射界面分布在这一深度的离散点上，即反射界面与这一深度线的交点为这一反射界面的分布点。 3）反射界面二次源的形成 ：绕射波的产生震源从炸点开始向下传播，遍历介质的每一个离散点绕射点。每遇到一个绕射点，都会激发该点的绕射波。反射界面与某深度离散点的 交叉点的反射系数不为零， 产生的绕射波也不为零， 而其他没有反射界面的地方，反 反射系数为零，是均匀介质， 产生的绕射波为零。 4）反射记录的形成把检波器的检波机制看成一个数学算子，或检波因子，可用脉冲或雷克子波。在延拓震源的同时，把检波器即检波因子也做完全相同的向下延拓，同样遍历每个离散点。震源和检波因子每延拓一个深度步长，震源、检波因子及反射系数三者做相关运算，则检波点就记录到了该深度每个绕射点的绕射波。所有二次源地震波都传播到测线某检波点，叠加而成该点的记录——炮集中的某一记录道。测线排列多少道检波器，则都可记录道所在位置的地震模拟记录。8、数字化根据数字信号处理的采样定理，把连续的信号变为计算机能处理的数字信号，使相移法正演模拟得以实现。频域抽样定理：一个频谱受限信号，如果时间只占据 的范围，若在频域以不大于频率间隔对信号的频谱采样，则抽样到的离散信号可以唯一表示原信号。时域抽样定理：一个时间受限信号，如果频谱只占据的范围，则信号可以用等间隔的抽样值唯一表示出来，而时间抽样间隔必须不大于，，。四、实验内容削波法相位移正演模拟（1）点绕射构造和水平层状速度模型（参数如图1所示）的正演数值模拟；