Rdseed学习笔记

Posted on | In

Rdseed软件是IRIS提供，用于读取SEED格式文件、并转化为SAC、ASCII、SEGY等
常用格式的转化软件。

由于Rdseed安装较为繁琐，在本电脑上未能成功实现安装。若安装成功，将在未来几个月内进行更新。

参考

[1] https://github.com/iris-edu-legacy/rdseed

[2] https://blog.seisman.info/rdseed-notes/

说明：
()内文字为[2]中的中文内容，[]中文字为[1]中README文件的中文直译。语法说明中其他中文补充内容为[1]中README文件中文直译。

语法说明

使用rdseed -h或者rdseed -u可查看命令选项和语法说明。

The following options may be selected:

  • -a retrieve the abbreviation dictionaries (提取缩略词词典)

  • -b n set the seismogram input buffer size (bytes) (设置输出波形数据的最大数据点数
    ，默认为2*E7)

  • -c retrieve volume table of contents (提取文件内容的目录信息)[检索卷目录]

  • -C arg retrieve the comments where ‘arg’ is either STN or CHN
    (提取台站或分量的注释信息) [检索 “arg”为STN或CHN的注释]

  • -d output data records (follow with -o to select a format or it defaults to SAC)
    (从seed数据中提取波形数据)[输出数据记录(后跟-o以选择格式，否则默认为SAC)]

  • -e extract event/station data
    (提取事件或台站数据到文件)[提取事件/站数据]

  • -E output filename will include endtime
    (输出的波形数据的文件名中包含结束时间)[输出文件名将包括结束时间]

  • -f file input file name
    (输入文件名)[输入文件名]

  • -g file specify alternate SEED volume for metadata (see ALT_RESPONSE_FILE)

    (为 SEED 或 miniSEED 数据单独指定响应文件。响应文件可以是 SEED 格式也可以是 dateless SEED 格式，也可以通过设置环境变量 ALT_RESPONSE_FILE 指定响应文件，这样 做的好处在于可以多个 SEED 文件共用一套响应文件。)

    [为元数据指定备用SEED卷(参见ALT_RESPONSE_FILE)]

  • -h help (this list)
    (显示语法)[显示这个用法列表]

  • -i ignore location codes
    (忽略位置码)[忽略位置代码]

  • -k strip records from output which have a zero sample count (miniSEED and SEED output only)
    (跳过数据点为0的记录)[采样计数为零的输出条带记录(仅限miniSEED和SEED输出)]

  • -l list contents of each record in the volume
    (列出每个 block 的内容)[列出卷中每条记录的内容]

  • -o n specify output format (numeric value, see examples below). Use -d with this option.
    (输出波形数据的格式。默认为SAC格式。n取值范围为1-9。其中 SAC（1）、 AH（2）、CSS（3）、miniSEED（4）、SEED（5）、SAC ASCII（6）、SEGY（7）、 Simple ASCII(SLIST)（8）、Simple ASCII(TSPAIR)（9）)

    [指定输出格式(数值，参见下面的示例)。在此选项中使用-d。]

  • -p output poles and zeros
    (输出 SAC PZ 格式的仪器响应文件)[输出极点和零点]

  • -q dir specify the output directory
    (指定输出目录，该目录必须已经存在，默认输出到当前目录)[指定输出目录为dir]

  • -Q qual specify data quality [E,D,M,Q,R] (see Input Options below for details)
    (选择波形数据质量)[qual指定数据质量E,D,M,Q,R]

    • E： 全部质量的波形数据
    • D： 不确定状态的数据
    • M： 已合并的数据
    • R： 原始波形数据
    • Q： 经过质量控制的数据

  • -R print response data
    (输出RESP格式的仪器响应文件)[打印响应数据]

  • -s retrieve all station header information
    (输出全部台站的RESP格式一起响应文件到终端)[检索所有站标头信息]

  • -S retrieve station summary information
    (提取台站汇总信息到文件rdseed.stations，内容包括台站名、台网名、经纬度、海拔、分量、台站开始时间和结束时间)

    [检索站摘要信息]

  • -t show start/stop times, start blocks of events
    (输出波形相关信息到终端，包括：台站名、分量名、台网名、位置码、质量控制符、波形开始时间、结束时间、采样率、数据点数。)

    [显示开始/停止时间，开始事件块]

  • -u show usage (this list)
    (显示语法)[显示使用方式]

  • -v n select volume number
    (选择卷号，默认值为1.对于seed文件n只能取1)[选择卷号]

    • 可以使用(-v选项)选择文件或磁带上所需的顺序卷。
      第一卷是“1”，第二卷是“2”，以此类推。备注:该选项仅适用于包含多个逻辑卷的物理卷。

  • -x use summary file
    (使用JWEED生成的summary文件，根据summary文件提取指定台站、分量和时间窗内的波形数据)[使用summary文件]

  • -z n check for reversals and apply sign change correction based on the number n=[0-3]
    (see Check Reversal under Options section below)
    (检查并校正数据极性)

    [检查反转情况，并根据数字n=[0-3]应用符号更改纠正

(请参阅下面选项部分的逆转检查)]

极性及其校正

参考于seisman博客

地震仪的传感器记录三分量值，若地面运动与仪器的敏感轴运动方向一致时，输出为正值。

SAC头段中cmpaz和cmpinc是用于描述正方向的方法。

cmpaz：敏感轴的方位角，从正北方向N开始顺时针旋转。

cmpinc：敏感轴的入射角，垂直向上方向开始旋转的角度。

下列方向对应的方位角和入射角参数：

方向 cmpaz（方位角） cmpinc（入射角） 方位码 极性
垂直向上 >0 >0 >Z >正
垂直向下 >0 >180 >Z >负
正北 >0 >90 >N >正
正南 >180 >90 >N >负
正东 >90 >90 >E >正
正西 >270 >90 >E >负

-z n选项可以用于检测和校正负极性。

  • n=0 不检测极性
  • n=1 只检查两个参数； 如果是负极性，就反转所有数据点的正负号并修改这两个参数的值。
  • n=2 只检查总增益的正负值；如果是负值表示负极性，反转所有数据的正负号，不修改这两个参数的值。
  • n=3 同时检查两个参数和总增益，只有其中之一符合负极性要求时做校正。

注意：两个参数指cmpaz和cmpnic

由于正负极性概念仅适用于6个标准分量方向，对于垂向分量而言，通常需要校正极性，负责查看Z分量波形数据时，可能出现某个台站波形极性错误。对于水平分量不做校正也不会有问题。
seisman建议使用n=1做极性校正。