人工合成地震记录程序设计
（一）、人工合成地震记录原理：
地震记录上看到的反射波波形是地震子波在地下各反射界面上发生反射时形成的。反射波的振幅有大有小（决定于界面反射系数的绝对值）、极性有正有负（取决于反射系数的正负）、到达时间有先有后（取决于反射界面的深度）的地震反射子波叠加的结果。
如果地震子波的波形用S(t)表示，地震剖面的反射系数为双程垂直反射时间t的函数，用R(t)表示，那么反射波地震记录形成的物理过程在数学上就可以用S(t)的R(t)的褶积表示，即某一时刻的反射波地震记录f(t)是:
f(t)?S(t)?R(t)
其离散形式为：
Mf(n?t)??S(m?t)?R((n?m)?t)
m?1如果大地为多层介质，在地面记录长度内可接收的反射波地震记录为：
f(n?t)??n?1N?S(m?t)?R((n?m)?t)
m?1M式中，n为合成地震记录的采样序号，n=1,2,3...N；N为合成一道地震记录的采样点数；m=1,2,3...M，为离散子波的采样点数；△t为采样间隔。
这种褶积模型将地震波的实际传播过程进行了简化：
1、在合成地震记录的过程中没有考虑大地的吸收作用，所有薄层的反射波都与地震子波的形式相同，只是振幅和符号不同。
2、假设地震波垂直入射到界面上，并原路径返回。
3、假设地层横向是均匀的，在深度（纵向）方向上假设密度为常数，只是速度发生变化。 4、不考虑地震波在传播过程中的透射损失。 （二）、人工合成地震记录的方法
1、 反射系数序列
在有速度测井资料的情况下，可以用速度曲线代替波阻抗曲线，计算反射系数序列。在没有速度资料的情况下，可根据干扰波调查剖面分析的结果设计地质模型。
如设计的地质模型如图a所示，图中H为层厚度，V为层速度，根据下式计算反射系数：
R(HN)?VN?VN?1
VN?VN?1式中H为反射界面的深度，N为反射层序号，随深度变化的反射系数序列如图b所示。但褶积计算中需要与时间有关的反射系数，深度与时间的转换可用下列公式计算：
R(H1)?R(n1),n1?2H1 V1?t2H2?n1 V2?tR(H2)?R(n2),n2?R(HN)?R(nN),nN?2HN?n1?n2???nN?1 VN?t式中n1、n2…为离散形式的时间序号，因序号需要用整数表示，计算出n1、n2…后需要取整处理。随时间变化的反射系数序列如图c所示。
V1=300m/s,H1=10m
V2=1000m/s,H210m V3=500m/s,H3=10m V4=3000m/s
H t 图a 地质模型 反射系数（H） 反射系数（t）t=n△t 2.地震子波
地震子波的选择是制作合成记录中的重要部分，子波选择是否合适，直接影响合成记录的质量。一般可用以下二种方法：
（1）根据已总结出的地震子波的特点，地震子波可以用一些函数表示，如最小相位子波：
S(m?t)?Asin(2?fm?t)?e?b?t
或简单的正弦波：
S(m?t)?sin(?2?fm?t)
（2）在地震记录上识别波形没有被其它波干涉的、形态完整的反射波。并在地震记录上读出离散数据，作为地震子波。
3.编程参数及要求
（1）程序设计的采样间隔△t=0.2ms，记录长度N为1024。
（2）选用正弦子波时，子波采样点数M=60个；选用最小相位子波时，子波采样点数为120，频率f=100，衰减系数b=120，振幅极值A自定。
（3）先用地质模型上给出的层参数编程，并绘图输出。