c语言地震子波6,地震子波显示及合成地震记录

一、实验目的

通过编制和运行相关C或matlab程序，理解地震子波的概念和特点，掌握地震合成记录的制作流程。认识子波，对子波的波形有直观的认识。(零相位子波，混合相位子波，最小相位子波；了解子波的分辨率与频宽的关系；)

利用褶积公式合成一维地震记录

二、主要内容和原理

1、理论子波的计算及显示

(1)零相位子波、最小相位子波、钟型子波的计算

零相位子波的表达式

最小相位子波的表达式为

；

n=m1/m2，为最大波峰m1和最大波谷m2之比，取2.0或1.5，本次试验采用2.0.；

代表子波的中心频率，t=i*dt，dt为时间采样间隔，i为时间离散点序号;

本次试验通过零相位子波和最小相位子波的公式，利用0.002s的采样间隔，对零相位子波和最小相位子波进行离散采样，得到一组数据。利用excel绘制图像，得到中心频率分别为10,25,40,100hz的8张图像。

2、利用声波测井资料计算层速度，密度和反射系数序列

本实验通过C语言编程序处理获得的声波测井资料，把资料中的数据用于层速的计算，公式为

和密度计算，当速度与密度呈线性关系时，计算公式为：

本次试验假设速度和密度为线性关系，采用这个公式参与计算。

3、 制作合成地震记录

地震合成记录是在一定的假设条件下进行的，如果地下实际介质有n层面都是反射界面的话，地面可以接收到每一层界面的反射波，于是一个实际地震记录道上会记录n个反射波。忽略介质吸收、透射损失影响，使用地震子波和反射系数离散卷积的方法得到。因为离散波形卷积需要采样间隔相同，所以使用0.002s的间隔对第i个反射面到地表的垂直往返时间进行离散采样.计算公式为：