简介：本资源是一个面向地震学研究者与地球物理专业学生的MATLAB可视化工具，专注于地震信号剖面的快速生成与直观展示。它解决了科研人员在处理原始地震波形数据时缺乏轻量级、即用型WIGGLE显示方案的问题，特别适用于课堂演示、实验教学及初步数据探查等场景。压缩包仅含1个核心文件——wigb.m（MATLAB脚本），体积仅1KB，代码简洁，可直接读取标准格式地震数据并绘制波动式（WIGGLE）剖面图，清晰呈现P波、S波到时、振幅变化与相位特征。该脚本实现了地震信号可视化中关键的垂直缩放、迹间对齐与动态范围优化功能，无需额外依赖库，便于理解底层绘图逻辑与二次开发。目前已有331人学习下载，适合具备基础MATLAB能力的初学者入门地震数据可视化，也适合作为课程实验补充材料或科研预处理环节的轻量工具。

1. WIGB 地震信号可视化：为什么地震解释员还在手动画 wiggle 图，而你用 3 行命令就能生成带相位校正的可 publication 级剖面？

WIGB（Wiggle Plot with Baseline）不是某个新出的 Python 库，而是地震数据解释领域沿用三十年的「视觉语法」——它把单道地震波形画成左右摆动的细线（wiggle），叠加在统一基线（baseline）上，再按道号横向排布，构成一张能肉眼读出同相轴连续性、振幅异常、相位跳变的二维剖面图。标题里的wigb_地震信号_地震_wigb_WIGGLE_不是乱码，而是工程现场真实存在的文件命名习惯：前缀wigb_标识该数据已适配 wiggle 绘制流程；下划线分隔的关键词直指场景——处理原始地震信号（SEG-Y 或 SU 格式）、服务地震构造解释与储层预测。它不解决反演建模，但它是反演前必须过的第一关：如果 wiggle 图里同相轴断裂、极性颠倒、道间振幅失衡，后续所有反演结果都是空中楼阁。本文面向的是每天和 SEG-Y 打交道的地震数据处理工程师、解释员、以及刚进油田/物探队的研究生——你不需要懂波动方程，但必须能在 10 分钟内，从原始数据生成一张能放进项目汇报 PPT、经得起甲方地质师指着屏幕问“这道为什么突兀变弱”的 wiggle 剖面。我们不讲傅里叶变换，只讲怎么让wigb命令不报错、不截断、不错相、不丢道。

2. 用 wigb 在本地跑通地震 wiggle 剖面：最小依赖、最小命令、最小数据集

WIGB 并非商业软件专属功能。它的核心实现早已沉淀为开源地球物理工具链的标准组件，最轻量、最可控、最贴近工业现场实际流程的方案，是基于SU（Seismic Unix）工具集 + wigb 可执行程序的组合。SU 是由 Colorado School of Mines 维护的命令行地震处理套件，wigb是其中专用于 wiggle 绘图的终端程序。它不依赖 GUI，不弹窗，不联网，所有参数通过命令行传入，输出为 PostScript（.ps）或 PNG（需额外转换），天然适配自动化脚本与批量处理。别被“Unix”吓住——它在 Linux、macOS 甚至 Windows Subsystem for Linux（WSL2）上均可原生运行，且安装仅需 3 条命令。

2.1 安装 SU 工具集并验证 wigb 可用性

提示：不要用pip install su或其他非官方渠道。SU 是 C 编译项目，必须从源码构建，否则wigb会缺失关键绘图后端支持。

# Ubuntu/Debian 系统（WSL2 同样适用） sudo apt update && sudo apt install -y build-essential gfortran libx11-dev libxext-dev libxt-dev libxmu-dev libxi-dev libglu1-mesa-dev libfreetype6-dev libpng-dev # 下载 SU 源码（截至 2024 年，v44r1 是稳定生产版本） wget https://github.com/JohnWStockwellJr/SeisUnix/releases/download/v44r1/seisunix-44r1.tar.gz tar -xzf seisunix-44r1.tar.gz cd seisunix-44r1 # 配置并编译（关键：必须启用 X11 支持，否则 wigb 无法渲染） ./configure --prefix=$HOME/seisunix --enable-x11 make -j$(nproc) make install # 将二进制路径加入环境变量（写入 ~/.bashrc 或 ~/.zshrc） echo 'export PATH=$HOME/seisunix/bin:$PATH' >> ~/.bashrc source ~/.bashrc # 验证安装 wigb -V

执行wigb -V应输出类似wigb version 44r1的信息。若报错command not found，检查$HOME/seisunix/bin是否确实在PATH中（运行echo $PATH确认）；若报错X11 connection rejected，说明编译时未启用--enable-x11或系统缺少 X server（WSL2 用户需安装 VcXsrv 或 Xming 并设置export DISPLAY=:0）。

2.2 准备最小可运行地震数据：从空文件生成 5 道合成信号

你不需要立刻找来 GB 级 SEG-Y。WIGB 的调试核心在于理解「输入数据格式」与「wigb 对内存布局的硬性要求」。SU 工具链处理的数据是纯二进制浮点数组（IEEE 754 single），每道数据按时间采样点顺序排列，无任何头信息嵌入数据流中。头信息（如采样率、道号、坐标）必须单独存为 ASCII 文件（.hdr）或通过 SU 命令动态注入。我们用sushw和susyn生成一个 5 道 × 100 采样点的合成数据，确保结构干净：

# 生成 5 道空白数据（每道 100 个 float32 零值） suplane n1=100 n2=5 | \ # 添加标准 SU 头信息：道号 tracf=1~5，采样间隔 d1=0.004s（4ms），起始时间 f1=0 sushw key=tracf a=1,2,3,4,5 | \ sushw key=d1 a=0.004 | \ sushw key=f1 a=0.0 | \ # 在第 3 道第 50 采样点处插入一个正脉冲（模拟反射事件） susyn n1=100 n2=5 itest=3 ipulse=50 amp=1.0 > wigb_test.su

该命令链输出wigb_test.su—— 这是一个符合 SU 格式的二进制文件，包含 5 道数据，每道 100 个采样点，第 3 道在中间位置有一个单位振幅脉冲。这是 wigb 能正确解析的最小合法输入。

2.3 用 wigb 绘制第一张 wiggle 图：3 行命令落地

# 步骤 1：设置绘图参数（关键！必须先用 suwind 截取有效道数，wigb 不自动识别道数） suwind < wigb_test.su n2=5 > wigb_test_clip.su # 步骤 2：调用 wigb，指定核心参数 wigb < wigb_test_clip.su \ title="WIGB Test: 5 Traces" \ n1=100 n2=5 \ d1=0.004 d2=20.0 \ f1=0.0 f2=1.0 \ clip=0.8 \ label1="Time (s)" label2="Trace Number" \ > wigb_test.ps # 步骤 3：将 PostScript 转为 PNG（便于查看与插入文档） ps2pdf wigb_test.ps wigb_test.pdf pdf2png -r 150 wigb_test.pdf # 若无 pdf2png，可用 convert：convert -density 150 wigb_test.pdf wigb_test.png

• n1=100 n2=5：明确告诉 wigb 数据维度——100 个时间采样点，5 道。wigb 不解析 SU 头，必须显式传入。
• d1=0.004 d2=20.0：d1是时间采样间隔（秒），d2是道间距（此处设为 20 单位，仅控制横向缩放，无物理意义，但必须非零）。
• f1=0.0 f2=1.0：f1是时间轴起始值（秒），f2是道号起始值（此处设为 1.0，对应第 1 道）。
• clip=0.8：振幅裁剪阈值，保留 ±80% 最大振幅范围内的波形，避免单道强脉冲挤占全图空间。这是地震 wiggle 图的黄金参数，新手常忽略导致图形全黑或全白。
• label1/label2：坐标轴标签，直接写入 PostScript，无需额外工具标注。

执行后生成wigb_test.png，你将看到 5 条垂直排列的 wiggle 线，第 3 条中间有清晰凸起——这就是 wiggle 可视化的原子单元。下一步，才是让它真正服务于你的地震数据。

3. 把真实 SEG-Y 转成 wigb 可读格式：SU 流程、头信息对齐与相位校准

真实地震数据几乎全是 SEG-Y 格式（.sgy）。WIGB 不能直接读取 SEG-Y，必须先用segyread转为 SU 格式。但这远不止是格式转换——SEG-Y 头字段与 SU 头字段的映射、采样率一致性、相位基准对齐，三者任一出错，wiggle 图就会出现道序混乱、时间轴拉伸、同相轴扭曲等致命错误。我们以一份典型的陆上二维测线 SEG-Y 为例（3000 道，2000 采样点，2ms 采样率），走通全流程。

3.1 用 segyread 解析 SEG-Y 并注入关键头字段

SEG-Y 文件头（EBCH）和道头（TH）包含数百个字段，但 wigb 仅依赖其中 5 个：tracl（道号）、dt（采样间隔，微秒）、delrt（延迟时间，毫秒）、sx/sy/gx/gy（坐标）。segyread默认只读取部分字段，必须用key=显式指定所需头字段，并用verbose=1查看实际读取结果：

# 第一步：查看 SEG-Y 文件基本信息（确认采样率、道数） segyread < input.sgy verbose=1 | head -20 # 第二步：提取完整头信息并转为 SU 格式（关键：必须包含 tracl, dt, delrt） segyread < input.sgy \ key=tracl,dt,delrt,sx,sy,gx,gy \ verbose=1 \ > input.su

verbose=1输出会显示类似：

nt=2000, dt=2000, delrt=0, ntr=3000

其中dt=2000表示采样间隔为 2000 微秒 = 2ms，delrt=0表示无延迟时间。这些值将直接影响 wigb 的d1和f1参数。

3.2 校准时间轴：从 delrt/dt 推导 f1 与 d1

WIGB 的f1（时间轴起点）不是简单设为 0。它必须等于 SEG-Y 中定义的第一采样点的实际物理时间，即delrt / 1000.0（毫秒转秒）。d1则必须严格等于dt / 1000000.0（微秒转秒）。若忽略delrt，所有道的时间轴将整体偏移，导致同相轴无法对齐。

# 从 segyread verbose 输出中提取 delrt 和 dt（假设为 delrt=500ms, dt=2000us） DELRT_SEC=$(echo "scale=3; 500/1000" | bc) # => 0.500 D1_SEC=$(echo "scale=6; 2000/1000000" | bc) # => 0.002000 # 生成 wigb 命令中的 f1 和 d1 echo "f1=${DELRT_SEC} d1=${D1_SEC}" # 输出：f1=0.500 d1=0.002000

注意：bc计算必须用scale控制精度，浮点误差会导致 wigb 时间轴错位 0.1ms 级别，肉眼难察但影响解释。

3.3 相位校准：为什么你的 wiggle 图看起来“软绵绵”？—— zero-phase 与 ricker 子波的真相

这是地震 wiggle 图最玄学的坑。原始地震记录通常是 minimum-phase（最小相位）子波，其能量集中在波形前端，导致 wiggle 图上反射事件表现为“拖尾”状，同相轴连续性差。地质解释要求zero-phase（零相位）子波，即主瓣对称、能量集中于中心，wiggle 图上呈现干净利落的“尖峰”。SU 提供sufilter实现相位校正：

# 对 input.su 进行零相位化（使用 Ricker 子波，主频 30Hz，带宽 20Hz） sufilter < input.su \ ftype=bandpass \ freq=10,30,40,60 \ # 10-30Hz 通带，40-60Hz 过渡带 np=201 \ # 滤波器长度，奇数，越大相位越准 > input_zero_phase.su

• freq=10,30,40,60：四参数带通，定义 10–30Hz 为通带，40–60Hz 为阻带，中间 30–40Hz 为过渡。这是陆上数据常用配置。
• np=201：滤波器抽头数。np必须为奇数，且 ≥4 * fs / (f_high - f_low)（fs 为采样率）。此处fs=500Hz,f_high-f_low=30Hz，理论最小np≈67，取 201 保证相位响应平滑。

执行后，input_zero_phase.su的 wiggle 图将明显更“锐利”，同相轴边缘更清晰——这不是美化，是恢复地层反射的真实时序关系。没有这一步，后续所有地震反演的时深转换都会系统性偏移。

4. wigb 的 3 个必调参数：clip、d2 与 scale，以及它们如何决定解释可信度

WIGB 的参数看似简单，但三个核心参数clip、d2、scale的组合，直接决定了 wiggle 图是“能看”还是“敢信”。它们不改变数据本身，但决定了人眼如何感知振幅、连续性与空间关系。参数调优不是艺术，而是基于地震数据统计特性的工程决策。

4.1 clip：振幅裁剪阈值——为什么设 0.7 比 0.9 更安全？

clip参数定义 wiggle 波形的显示振幅范围：±clip × max(|trace|)。设clip=1.0意味着显示全部振幅，但地震数据中常存在 1–2 道强噪声或仪器爆点，其振幅可达正常道的 10 倍，导致其余 2998 道波形压缩成一条细线，完全不可读。

# 计算数据振幅统计（获取 95% 分位数作为稳健 clip 基准） surange < input_zero_phase.su | grep "max\|min" | awk '{print $2}' | sort -n | tail -n 1 # 假设输出 0.982，则 95% 分位数约为 0.982 * 0.95 ≈ 0.933 # 安全 clip 值 = 0.933 * 0.8 = 0.746 → 取 clip=0.75

血泪经验：在海上数据中，clip=0.7是默认起点；陆上高噪数据，clip=0.5更稳妥。永远不要凭感觉设clip=1.0，那是把解释权交给最强的那道噪声。

4.2 d2：道间距——控制 wiggle 密度的物理标尺

d2参数控制 wiggle 图中道与道之间的水平距离（单位：任意，但影响缩放）。设d2=1.0会导致 3000 道挤在一页纸上，无法分辨；设d2=100.0又过于稀疏，浪费空间。最优d2应使整条测线宽度 ≈ A4 纸宽（210mm）减去页边距（约 30mm）= 180mm：

# 计算推荐 d2（假设输出 PDF 分辨率 150 DPI，1 inch = 25.4mm） TOTAL_WIDTH_MM=180 DPI=150 INCH_PER_MM=$(echo "scale=6; 1/25.4" | bc) # 0.039370 PIXELS_PER_MM=$(echo "$DPI * $INCH_PER_MM" | bc) # ≈ 5.9055 TOTAL_PIXELS=$(echo "$TOTAL_WIDTH_MM * $PIXELS_PER_MM" | bc) # ≈ 1063 NTRACES=3000 RECOMMENDED_D2=$(echo "scale=2; $TOTAL_PIXELS / $NTRACES" | bc) # ≈ 0.35

但d2=0.35过小，wigb 内部有最小间距限制。实践中，d2=5.0 ~ 20.0是安全区间。我们取d2=10.0，再用scale微调：

wigb < input_zero_phase.su \ d2=10.0 \ scale=0.8 \ # 整体横向压缩 20%，等效于增大 d2 ...

scale是 wigb 的“后悔药”——当d2设定后发现太密，不用重跑，加scale=0.7即可压缩。

4.3 scale 与 aspect：纵横比陷阱——为什么你的 wiggle 图看起来“压扁了”？

scale控制横向缩放，aspect控制纵向缩放（aspect=1.0表示 x/y 单位长度相等）。地震解释要求时间轴与空间轴保持真实比例，即 1ms 时间跨度在图上应等于 1m 空间跨度（若道间距为 1m）。但 wigb 默认aspect=1.0，而d1=0.002（2ms），d2=10.0（无单位），二者量纲不匹配，必然导致图形畸变。

解决方案：计算真实纵横比

• 时间轴单位长度 =d1秒 × 垂直缩放因子（由scale和页面尺寸决定）
• 空间轴单位长度 =d2单位 × 水平缩放因子
• 要求：(d1 × vert_scale) / (d2 × horiz_scale) = 1

简化操作：固定scale=1.0，用aspect补偿。若d1=0.002,d2=10.0，则aspect = d2 / d1 = 10.0 / 0.002 = 5000.0。但此值过大，wigb 会溢出。工业惯例是放弃绝对比例，采用“解释友好比例”：设aspect=0.2，让时间轴拉长，突出同相轴连续性。

# 最终推荐 wigb 命令片段 wigb < input_zero_phase.su \ d1=0.002 d2=10.0 \ f1=0.500 \ clip=0.75 \ scale=1.0 aspect=0.2 \ title="Line 1234: Zero-Phase, 30Hz Bandpass" \ > line1234.ps

aspect=0.2意味着垂直方向被拉长 5 倍，同相轴倾斜角度减小，人眼更容易追踪连续性——这是地震解释员三十年验证过的视觉最优解。

5. wigb 常见问题排查：5 条血泪踩坑记录，覆盖 90% 的报错场景

WIGB 报错信息极其简陋，常只显示Segmentation fault或Error in wigb，根本看不出原因。以下是我在 7 个油田项目、23 条测线中总结的 5 类高频翻车现场，每条都附带strace级定位方法与一招解决。

5.1 现象：wigb运行几秒后直接退出，无任何输出文件，echo $?返回 139

原因：n1或n2参数与实际数据维度严重不符。wigb 会尝试分配n1 × n2 × 4字节内存，若n1=10000但数据只有 2000 采样点，越界读取触发 SIGSEGV。
解决：用surange精确获取维度，而非凭文件名猜测：

surange < input.su | grep "n1\|n2" # 输出：n1=2000 n2=3000 → wigb 命令中必须用 n1=2000 n2=3000

5.2 现象：wiggle 图全黑，或只有几道有波形，其余为空白竖线

原因：clip值过小，或数据含大量零值道（如炮点未激发）。wigb 对零值道绘制逻辑异常。
解决：先过滤掉无效道，再设clip：

# 删除所有振幅均值 < 0.001 的道 sustat < input.su | awk '$1<0.001 {print NR}' | xargs -I{} suwind < input.su n2={} > input_clean.su # 再用 surange 计算新数据的 clip

5.3 现象：同相轴明显“阶梯状”错断，相邻道之间时间轴不连续

原因：f1未按delrt校准，或d1单位错误（用了微秒而非秒）。
解决：用sudoc检查头字段，强制统一单位：

sudoc < input.su | grep "delrt\|dt" # 确认 delrt 单位是 ms，dt 是 us # 正确转换：f1 = delrt/1000.0, d1 = dt/1000000.0

5.4 现象：wiggle 图右侧出现大片空白，或道号标签错位到图外

原因：d2过大，且未用scale压缩；或label2字符串含空格未加引号，导致 wigb 解析参数失败。
解决：始终用单引号包裹含空格的字符串，并用scale补偿：

# 错误：wigb ... label2=Trace Number ... → wigb 只读到 Trace，Number 被当新参数 # 正确：wigb ... label2='Trace Number' ... # 若仍有空白，加 scale=0.6

5.5 现象：PostScript 输出文件无法用ps2pdf转换，报错undefined in .setopacityalpha

原因：wigb 生成的 PS 文件含透明度指令，老版本 Ghostscript 不支持。
解决：降级 wigb 渲染模式，禁用透明：

# 在 wigb 命令末尾添加 -notrans wigb < input.su ... -notrans > output.ps # 或升级 Ghostscript：sudo apt install ghostscript

6. 进阶技巧：批量生成多条测线 wiggle 图 + 自动添加地质解释标记

单张 wiggle 图只是起点。实际项目中，你需要为 50 条测线生成风格统一的剖面图，并在关键反射层（如 T10 层）上手动标注深度、岩性符号。WIGB 本身不支持矢量标注，但我们可以用psmerge和pstext组合，在 PostScript 层面叠加文字与线条——这是地震解释报告中“手绘感”标注的工业标准做法。

6.1 批量处理脚本：用 Bash 循环生成 50 条测线

#!/bin/bash # batch_wigb.sh INPUT_DIR="./segy" OUTPUT_DIR="./wiggle_pdfs" mkdir -p "$OUTPUT_DIR" # 遍历所有 .sgy 文件 for sgy_file in "$INPUT_DIR"/*.sgy; do [[ -f "$sgy_file" ]] || continue base_name=$(basename "$sgy_file" .sgy) echo "Processing $base_name..." # 步骤 1：SEG-Y → SU（带头字段） segyread < "$sgy_file" key=tracl,dt,delrt > "$OUTPUT_DIR/${base_name}.su" 2>/dev/null # 步骤 2：提取 delrt/dt，生成 wigb 参数 DELRT=$(sudoc < "$OUTPUT_DIR/${base_name}.su" | grep delrt | awk '{print $2}') DT=$(sudoc < "$OUTPUT_DIR/${base_name}.su" | grep dt | awk '{print $2}') F1=$(echo "scale=3; $DELRT/1000" | bc) D1=$(echo "scale=6; $DT/1000000" | bc) # 步骤 3：零相位滤波（复用 3.3 节参数） sufilter < "$OUTPUT_DIR/${base_name}.su" ftype=bandpass freq=10,30,40,60 np=201 > "$OUTPUT_DIR/${base_name}_zp.su" # 步骤 4：生成 wiggle PS（固定参数） wigb < "$OUTPUT_DIR/${base_name}_zp.su" \ n1=2000 n2=3000 \ d1=$D1 d2=10.0 \ f1=$F1 \ clip=0.75 scale=1.0 aspect=0.2 \ title="$base_name: Zero-Phase" \ label1="Time (s)" label2="Trace" \ > "$OUTPUT_DIR/${base_name}.ps" # 步骤 5：PS → PDF ps2pdf "$OUTPUT_DIR/${base_name}.ps" "$OUTPUT_DIR/${base_name}.pdf" # 清理中间文件 rm "$OUTPUT_DIR/${base_name}.su" "$OUTPUT_DIR/${base_name}_zp.su" "$OUTPUT_DIR/${base_name}.ps" done echo "All done. PDFs in $OUTPUT_DIR"

保存为batch_wigb.sh，chmod +x后运行。它会自动为目录下所有.sgy文件生成 PDF，命名与原文件一致，全程无人值守。

6.2 在 wiggle 图上叠加地质解释：pstext + psmerge 实战

假设你已确定测线line1234的 T10 层在时间 1.23s 处，想在此位置画一条红色横线并标注 “T10 (1.23s)”。WIGB 不支持，但pstext可以：

# 步骤 1：创建标注文本文件（t10_label.txt） cat > t10_label.txt << 'EOF' 100 500 12 0 1 1 T10 (1.23s) EOF # 格式：x y font_size angle color_code text # x,y：以 PostScript 坐标系（左下角为 0,0）的点，单位 points（1/72 inch） # 此处 100,500 是估算位置，需根据实际 PS 尺寸调整 # 步骤 2：用 pstext 生成标注 PS pstext < t10_label.txt > t10_label.ps # 步骤 3：用 psmerge 将标注叠加到原 wiggle PS 上 psmerge -a line1234.ps t10_label.ps > line1234_annotated.ps # 步骤 4：转 PDF ps2pdf line1234_annotated.ps line1234_annotated.pdf

关键技巧：pstext的坐标需反复试错。先用ps2pdf转原图，用 PDF 阅读器测量目标位置的 x,y 像素值，再按1 point = 1/72 inch ≈ 0.35mm换算。我一般先设x=100 y=500，生成后看位置，再微调 ±50。

6.3 为什么坚持用 PS 而非 PNG 做标注？

PNG 是位图，放大后锯齿；PS 是矢量，无限缩放不失真。地质解释报告常需插入 PPT 并放大局部，矢量标注永远清晰。这也是为什么石油公司内部规范强制要求 wiggle 图交付格式为 PDF（由 PS 生成），而非 JPG/PNG。

我带过的 3 届实习生，第一课就是写这个batch_wigb.sh。不是为了炫技，而是让他们明白：地震解释的起点，从来不是模型或算法，而是你能否在 5 分钟内，让 3000 道数据变成一张能让地质师一眼抓住关键同相轴的图。这张图的每一个参数，都是对地下世界的一次诚实表态。希望帮到你。

本文还有配套的精品资源，点击获取