简介本资源为地震学分析软件SACSeismic Analysis Code的Linux x86_64版本安装包主要面向从事地震波形数据处理与格式转换的科研人员及学生。包内包含SAC 101.5c及101.6a等版本可满足SAC、MiniSEED与SEED三种地震数据格式之间的相互转换需求便于将数据接入不同分析平台或数据库。压缩包共550个文件以m、z、e、n、c、f等源码与头文件为主辅以csh、sh脚本、readme说明、tbl表格及py辅助脚本整体约3.59MB目录结构清晰便于按模块查阅与编译安装。目前已有689人学习下载。借助该工具读者可完成地震数据的读取、处理、可视化与格式互转解决跨软件、跨平台的数据兼容问题适合需要灵活处理地震波形数据的中高级用户参考使用。1. 拿到 sac-101.5c-linux_x86_64.tar.gz 之后它是什么能解决什么如果你在 Linux 上做地震信号处理、台站波形分析或者 SAC 格式数据的批处理大概率绕不开 SACSeismic Analysis Code这套老牌工具。sac-101.5c-linux_x86_64.tar.gz就是 SAC 在 Linux x86_64 平台上的一个预编译发行包解压后直接得到可执行文件不需要从源码一步步编译。同系列的sac-101.6a-linux属于更新的版本分支而seed_sac通常指配套的 SEED 格式与 SAC 格式互转工具链。这套资源解决的核心问题是让你在一台干净的 Linux 机器上用最短路径跑起 SAC 命令行环境完成波形读取、滤波、去均值、写 SAC 文件、批量脚本处理等操作。适合地震学方向的研究生、台网运维人员以及需要把 SAC 集成进自动化流程的工程师。它不依赖图形界面纯命令行适合服务器和嵌入式 Linux 场景。2. 解压、环境变量与首次运行把 tar.gz 变成能用的命令2.1 解压前先确认三件事拿到sac-101.5c-linux_x86_64.tar.gz之后别急着tar -zxvf。先确认三件事第一你的系统是 x86_64 架构uname -m输出x86_64才匹配第二确认动态链接库依赖SAC 预编译包通常依赖libc、libm、libX11即使命令行模式也可能链接到 X11 库用ldd检查解压后的二进制第三确认磁盘空间SAC 本体加示例数据通常几百 MB但如果你要挂载 NAS 存储跑大批量数据预留空间要按波形数据量算。常见做法是先在/opt或用户目录下建一个独立目录避免和系统包管理器冲突。我一般会这样操作# 创建安装目录避免污染系统路径 mkdir -p /opt/sac-101.5c # 解压到指定目录注意 tar.gz 内部通常有一层顶层目录 tar -zxvf sac-101.5c-linux_x86_64.tar.gz -C /opt/sac-101.5c --strip-components1 # 查看解压后的目录结构 ls -la /opt/sac-101.5c--strip-components1的作用是去掉压缩包内可能存在的顶层目录直接把内容铺到目标目录。如果解压后发现多了一层sac-101.5c/说明压缩包本身没有顶层目录去掉这个参数重新解压即可。解压后你应该能看到bin/、lib/、aux/这类目录bin/里就是sac可执行文件。2.2 环境变量配置与验证SAC 运行依赖两个关键环境变量SACAUX指向辅助数据目录包含仪器响应、台站信息等PATH里加入bin目录。如果你用的是 bash在~/.bashrc或/etc/profile.d/sac.sh里写# SAC 安装根目录 export SAC_ROOT/opt/sac-101.5c # 辅助数据目录SAC 启动时会读取这里的台站和响应文件 export SACAUX$SAC_ROOT/aux # 把 sac 可执行文件加入 PATH export PATH$SAC_ROOT/bin:$PATH # 可选设置默认的 SAC 宏文件搜索路径 export SACMACRO$HOME/.sacmacros配置完成后执行source ~/.bashrc然后运行sac进入交互式命令行。如果看到SAC提示符说明安装成功。输入quit退出。如果报错error while loading shared libraries: libX11.so.6说明缺少 X11 运行库在 CentOS 7 上可以yum install libX11在 Debian/Ubuntu 上apt install libx11-6。这是预编译包最常见的翻车点因为很多服务器最小化安装不带 X11 库。2.3 用 seed_sac 做格式转换的典型流程seed_sac相关工具通常用于把 SEED 格式的连续波形切成 SAC 文件。假设你有一个 SEED 文件data.seed和对应的台站信息stations.txt常见做法是先用rdseed或配套工具解出 SAC 文件再用 SAC 做后续处理。如果你拿到的资源里包含seed_sac脚本或二进制典型调用方式如下# 从 SEED 文件提取 SAC 波形-f 指定输出格式为 SAC seed_sac -f SAC -o ./sac_out data.seed # 检查输出目录确认每个台站通道生成了独立的 SAC 文件 ls ./sac_out/*.sac | head -20参数-f SAC指定输出格式-o指定输出目录。不同版本的seed_sac参数名可能略有差异用seed_sac -h查看实际支持的选项。转换完成后用 SAC 读取一个文件验证# 进入 SAC 交互环境 sac # 读取一个 SAC 文件 SAC r ./sac_out/BJN.BHZ.sac # 查看头段信息确认台站、通道、采样率正确 SAC lh # 画波形图需要图形终端或 X11 转发 SAC plot # 退出 SAC quitlh列出头段变量重点看kstnm台站名、kcmpnm通道名、delta采样间隔、npts采样点数。如果delta为 0 或npts异常说明转换过程中头段信息丢失需要检查 SEED 文件的完整性。3. SAC 命令行核心操作从读波形到批量处理脚本3.1 读、写、滤波、去均值的基本命令链SAC 的命令行交互模式适合快速验证但真正干活要靠脚本。先掌握基本命令链读文件用r写文件用w滤波用filter去均值用rmean去趋势用rtrend。一个典型的预处理流程如下# 启动 SAC sac # 读取三个分量的波形 SAC r ./data/*.BHZ ./data/*.BHN ./data/*.BHE # 去均值 SAC rmean # 去线性趋势 SAC rtrend # 带通滤波频率范围 0.5-10 Hz SAC filter bandpass 0.5 10 # 写回 SAC 文件文件名加 _proc 后缀 SAC w over ./proc/*.sac SAC quitrmean去掉直流偏移rtrend去掉线性趋势这两个操作在波形预处理里几乎是标配。filter bandpass 0.5 10的四个参数含义是滤波类型bandpass低频截止 0.5 Hz高频截止 10 Hz。注意 SAC 的滤波默认是 Butterworth 滤波器阶数可以通过filter bandpass 0.5 10 npasses 2调整npasses表示滤波次数次数越多边缘效应越明显一般 2 次足够。3.2 用宏文件把重复操作固化下来每次手动敲命令不现实SAC 支持宏文件macro把命令序列写进.mac文件用m命令调用。比如建一个preprocess.mac# preprocess.mac - SAC 预处理宏 # 用法SAC m preprocess.mac r %1 rmean rtrend filter bandpass 0.5 10 w over %2调用时传入输入和输出文件模式SAC m preprocess.mac ./data/*.sac ./proc/*.sac%1和%2是宏参数占位符。这种写法适合批量处理但要注意 SAC 的宏参数替换是文本级的路径里有空格会出问题。我一般会在脚本里用for循环逐个文件处理而不是依赖通配符一次性读入因为一次性读入多个文件时w over的命名规则容易乱。3.3 用 shell 脚本驱动 SAC 批处理更可控的方式是用 shell 脚本生成 SAC 命令流然后管道传给sac。比如批量处理一个目录下所有 SAC 文件#!/bin/bash # batch_sac.sh - 批量预处理 SAC 文件 INPUT_DIR./data OUTPUT_DIR./proc mkdir -p $OUTPUT_DIR for sacfile in $INPUT_DIR/*.sac; do basename$(basename $sacfile) # 生成 SAC 命令流通过管道传给 sac cat EOF | sac r $sacfile rmean rtrend filter bandpass 0.5 10 w $OUTPUT_DIR/$basename quit EOF done这种写法的好处是每个文件独立处理不会因为某个文件头段异常导致整批失败。cat EOF | sac是 here-document 语法把多行命令喂给 SAC 的标准输入。注意w命令后面直接跟输出路径不要加over因为这里是逐个文件写不需要覆盖模式。3.4 头段变量修改与时间窗截取实际工作中经常需要修改头段变量比如统一台站名、调整参考时间、截取特定时间窗。SAC 提供chchange header命令SAC r ./data/BJN.BHZ.sac # 修改台站名为 BJN通道为 BHZ SAC ch kstnm BJN kcmpnm BHZ # 截取从参考时间后 10 秒到 60 秒的数据 SAC cut 10 60 # 查看截取后的采样点数 SAC lh npts SAC w ./proc/BJN.BHZ.cut.sac SAC quitcut 10 60的参数是相对于参考时间o标记的秒数。如果参考时间不是 0需要先用ch o 0把参考时间归零否则截取窗口会偏移。这是新手最容易踩的坑之一读进来的 SAC 文件参考时间可能是发震时刻直接cut会得到意料之外的结果。4. 避坑与排查SAC 在 Linux 上最容易翻车的五个场景4.1 现象运行 sac 报 libX11.so.6 找不到原因预编译包在链接时依赖 X11 库即使你只用命令行模式动态链接器仍然会加载它。最小化安装的 CentOS 7 或 Debian 默认不带 X11 运行库。解决安装 X11 运行库即可不需要装完整桌面环境。CentOS 7 执行yum install libX11 libXextDebian/Ubuntu 执行apt install libx11-6 libxext6。如果服务器完全不允许装图形库可以考虑从源码编译 SAC 并禁用 X11 支持但预编译包没有这个选项。4.2 现象SAC 启动后 SACAUX 路径为空读仪器响应文件失败原因SACAUX环境变量没有设置或者指向的目录不存在。SAC 启动时不会自动推断辅助数据目录必须显式指定。解决确认echo $SACAUX输出非空且目录存在。如果目录存在但文件缺失检查解压时是否完整解出了aux/目录。有些压缩包为了减小体积会省略aux/下的部分文件需要单独下载补充。4.3 现象批量处理时部分文件写出的 SAC 头段 delta 为 0原因输入文件的采样间隔头段本身有问题或者r命令读取时因为文件格式不标准导致头段解析失败。SEED 转 SAC 过程中如果台站信息不完整delta可能丢失。解决在预处理脚本里加一步检查用lh delta判断如果为 0 则跳过或报错。更根本的办法是回到seed_sac转换环节确认台站响应文件完整。我一般会在批处理脚本开头加一个sac命令流专门检查头段# 检查 delta 是否为 0 cat EOF | sac r $sacfile lh delta quit EOF4.4 现象filter bandpass 后波形两端出现明显畸变原因Butterworth 滤波器在数据边缘会产生瞬态响应滤波次数越多、频带越窄边缘畸变越严重。这是数字滤波的固有特性不是 bug。解决在滤波前先截取比目标时间窗更长的数据滤波后再截取目标窗。或者使用 SAC 的taper命令在数据两端加窗减少边缘突变。常见做法是taper后再filter代价是两端数据会被压低需要根据实际需求权衡。4.5 现象在 NAS 挂载目录上跑 SAC 批处理速度极慢甚至卡死原因SAC 的r和w命令对文件系统的随机读写性能敏感NAS 的网络延迟和锁机制会导致大量小文件读写时性能急剧下降。解决先把数据从 NAS 拷贝到本地磁盘再处理处理完再拷回。如果数据量太大无法全部拷贝至少把当前批次的文件预取到本地/tmp或 SSD 目录。我一般会在脚本里加一个rsync预取步骤处理完再rsync回 NAS。5. 进阶技巧用 SAC 做波形质量检查与自动化验证5.1 用 SAC 命令流做批量质量检查批量处理最怕的是“跑完了但不知道结果对不对”。一个实用技巧是写一个质量检查宏对每个处理后的文件输出关键指标最大振幅、均方根、过零率。这些指标可以帮你快速定位异常波形。# qc.mac - 波形质量检查宏 # 输出文件名 最大振幅 RMS 过零率 r %1 # 计算最大振幅 mt max # 计算均方根 mt rms # 计算过零率 mt zcross quitmt是 SAC 的测量命令max返回最大值rms返回均方根zcross返回过零率。这些值会打印到标准输出重定向到日志文件后可以用awk做统计分析。比如找出 RMS 超过阈值 3 倍标准差的文件# 批量运行质量检查输出到日志 for f in ./proc/*.sac; do echo $f qc.log cat EOF | sac qc.log 21 r $f mt max mt rms mt zcross quit EOF done # 用 awk 提取 RMS 值并标记异常 grep -A1 rms qc.log | grep -v ^-- | awk {print $2} | sort -n | tail -55.2 用 SAC 的 cut 和 interpolate 做时间对齐多台站数据对比时时间对齐是关键。SAC 提供cut和interpolate命令可以把不同采样率的波形统一到相同时间基准。常见做法是先用cut截取公共时间窗再用interpolate重采样到统一采样率SAC r ./data/*.sac # 截取所有波形的公共时间窗 SAC cut 0 120 # 统一重采样到 100 Hz SAC interpolate delta 0.01 # 写出一致采样率的文件 SAC w over ./aligned/*.sac SAC quitinterpolate delta 0.01表示重采样到采样间隔 0.01 秒即 100 Hz。注意interpolate默认使用线性插值如果对精度要求高可以指定interpolate delta 0.01 akima使用 Akima 插值。重采样后头段delta会自动更新但npts会变化写文件前最好用lh确认一下。5.3 一个我常用的验证习惯每次跑完批处理我不会直接看波形图而是先跑一遍质量检查脚本把 RMS 和过零率排序挑出偏离中位数最远的几个文件单独看。这个习惯帮我省了很多“后悔药”——有次一批数据里混进了几个仪器响应文件没去掉的波形振幅比正常值大两个数量级质量检查脚本直接把它们排到了最前面。从那以后我每次批量处理 SAC 数据都强制走一遍mt max和mt rms检查确认没有异常值再进入下一步分析。希望帮到你。本文还有配套的精品资源点击获取
