2. 地震子波的基本概念与数学描述

在地震勘探中，震源激发后产生的并不是一个理想的尖脉冲，而是一段具有有限持续时间、特定频率成分和相位特征的振荡波形。这个从震源出发、经地层传播后到达检波器的基本波形单元，就是地震子波。理解子波的定义、分类和数学描述，是后续进行子波提取和合成地震记录制作的理论基础。

子波的定义与物理意义

地震子波可以定义为：震源激发产生、经大地滤波作用后，具有确定起始时间、有限能量和特定形态的地震波形。它是构成地震道的基本单元，地震道可视为地下各反射界面反射系数序列与子波的褶积结果。

从物理意义上看，子波包含了三个关键信息：振幅反映反射界面的波阻抗差异大小；极性反映波阻抗变化的方向（增大或减小）；时间位置反映反射界面的深度（双程旅行时）。子波的形态则受震源类型、地层吸收衰减、检波器响应等多种因素共同影响。理想情况下，我们希望子波能量集中、旁瓣小、分辨率高，但实际子波往往存在一定的延续长度和振荡拖尾。

子波的分类

按相位特征，地震子波主要分为三类：

• 零相位子波：能量集中在时间原点，波形关于中心对称，旁瓣最小，分辨率最高。其频谱为实数，相位谱恒为零。零相位子波在地震资料处理中最为理想，因为它的峰值直接对应反射界面位置，便于解释。
• 最小相位子波：能量集中在子波前部，相位谱为最小相位，即其z变换的零极点全部位于单位圆内。最小相位子波具有因果性，物理上最易实现，实际震源（如炸药震源）产生的子波往往接近最小相位。
• 混合相位子波：介于零相位和最小相位之间，能量分布和相位特征兼具两者特点。实际地震资料中的子波多为混合相位，这也是子波提取中需要估计相位谱的原因。

此外，按时间域形态还可分为雷克子波（Ricker wavelet）、带通子波等；按频率成分可分为宽带子波和窄带子波。雷克子波是最常用的理论子波模型，其表达式为：

w(t) = (1 - 2π²f₀²t²) · exp(-π²f₀²t²)

其中 f₀ 为子波主频，t 为时间。该子波为零相位，峰值在 t=0 处，主瓣宽度与主频成反比。

子波的数学描述

子波的数学描述主要从时间域和频率域两个角度展开。时间域上，子波 w(t) 是一个有限能量的振荡波形；频率域上，其频谱 W(f) 可表示为振幅谱 A(f) 和相位谱 φ(f) 的组合：

W(f) = A(f) · exp[jφ(f)]

振幅谱决定子波的频率成分和频带宽度，相位谱决定子波的形态和能量分布。对于零相位子波，φ(f)=0；对于最小相位子波，φ(f) 满足最小相位条件；混合相位子波则具有更复杂的相位谱。

地震道 s(t) 与反射系数 r(t)、子波 w(t) 的关系为褶积模型：

s(t) = r(t) * w(t) + n(t)

其中 n(t) 为噪声。在频率域中，该关系变为乘积：S(f) = R(f)·W(f) + N(f)。这一模型是子波提取和合成地震记录制作的核心依据——已知子波和反射系数可合成地震道，反之从地震道中估计子波即为子波提取。

子波参数与分辨率的关系

子波的主频和频带宽度直接决定地震分辨率。以雷克子波为例，主频 f₀ 越高，子波主瓣越窄，纵向分辨率越高；但高频成分在地层中衰减更快，深层资料往往主频偏低。频带宽度则影响子波的旁瓣水平，频带越宽，旁瓣越小，波形越接近理想脉冲。

在实际工作中，常通过以下参数描述子波特征：主频（峰值频率）、频带范围（如 10–60 Hz）、视周期、延续时间、旁瓣比等。例如，一个主频 30 Hz 的雷克子波，其主瓣宽度约为 1/(2f₀) ≈ 16.7 ms，这意味着在速度 3000 m/s 的地层中，理论纵向分辨率约为四分之一波长，即约 25 m。这些估算值属于教学示例，实际应用中需结合具体资料标定。

从概念到实践：子波描述在合成记录中的作用

在制作合成地震记录时，通常先由声波和密度测井曲线计算波阻抗，再求取反射系数序列，然后与选定的子波进行褶积。子波的选择直接影响合成记录的波形形态和与实测地震道的匹配效果。若子波主频偏高，合成记录会出现过多高频细节，与实际资料频带不符；若相位选择错误，则合成记录与地震道的同相轴会出现错位。