PZ文件的几点说明
- 完整的
PZ文件包含注释和零点、极点、放大系数两部分
* **********************************
* NETWORK (KNETWK): IU
* STATION (KSTNM): COLA
* LOCATION (KHOLE): 00
* CHANNEL (KCMPNM): BHZ
* CREATED : 2013-06-22T14:12:09
* START : 2012-09-14T04:00:00
* END : 2599-12-31T23:59:59
* DESCRIPTION : College Outpost, Alaska, USA
* LATITUDE : 64.873599
* LONGITUDE : -147.861600
* ELEVATION : 84.0
* DEPTH : 116.0
* DIP : 0.0
* AZIMUTH : 0.0
* SAMPLE RATE : 20.0
* INPUT UNIT : M
* OUTPUT UNIT : COUNTS
* INSTTYPE : Geotech KS-54000 Borehole Seismometer
* INSTGAIN : 2.013040e+03 (M/S)
* COMMENT : N/A
* SENSITIVITY : 3.377320e+09 (M/S)
* A0 : 8.627050e+04
* **********************************
ZEROS 3
+0.000000e+00 +0.000000e+00
+0.000000e+00 +0.000000e+00
+0.000000e+00 +0.000000e+00
POLES 5
-5.943130e+01 +0.000000e+00
-2.271210e+01 +2.710650e+01
-2.271210e+01 -2.710650e+01
-4.800400e-03 +0.000000e+00
-7.384400e-02 +0.000000e+00
CONSTANT +2.913631e+14
-
*开头的行为注释行，并不要求这些注释行全部存在，即可以删除部分注释行。虽然我们称之为注释行，PZ文件中一旦包含这些注释行，在去仪器响应过程中，sac会根据这些注释行的内容检查数据的头段变量，如果不匹配，sac会报错。 - 一个PZ文件中可包含多个台站，多个分量的仪器响应信息，但必须包含注释行（同样，可以省略一部分注释行，但至少要包含STATION和CHANNEL 信息，以便区分不同台站，不同分量的仪器响应信息）
如：
* **********************************
* STATION (KSTNM): COLA
* CHANNEL (KCMPNM): BHZ
* **********************************
ZEROS 3
....
POLES 4
....
* **********************************
* STATION (KSTNM): COAA
* CHANNEL (KCMPNM): BHZ
* **********************************
ZEROS 3
....
POLES 4
sac中使用PZ文件去除仪器响应文件包含两种方式：
指定PZ文件
trans FROM polezero subtype mypz.pz TO none
使用整合的PZ文件（多台站、多分量，多时段的[如中途更换过仪器]PZ信息）
trans FROM polezero s all.pz TO none
使用该方式，可以一次性读入较多的数据，sac自动去匹配每个文件对应的PZ信息。
- 严格来说，保留PZ文件中的全部注释信息是最好的，但为防止SAC数据道头中信息缺失而导致匹配PZ信息错误，我们可以删除PZ文件中的全部注释信息（此时，只能使用指定PZ文件的方式
subtype），或者仅保留部分注释信息
* **********************************
* NETWORK (KNETWK): IU
* STATION (KSTNM): COLA
* LOCATION (KHOLE): 00
* CHANNEL (KCMPNM): BHZ
* CREATED : 2013-06-22T14:12:09
* START : 2012-09-14T04:00:00
* END : 2599-12-31T23:59:59
* DESCRIPTION : College Outpost, Alaska, USA
* INSTTYPE : Geotech KS-54000 Borehole Seismometer
* **********************************
为便于检查，建议一个台站制作一个PZ文件，同时使用部分注释信息。这样就可应对三个分量（当然，现在大部分地震计的三分量仪器响应是一致的），中途更换仪器这些复杂情况。DESCRIPTION与INSTTYPE是很有用的两行注释，可以帮助我们检查制作PZ文件过程是否搞混仪器型号。
参考：
http://seisman.info
SAC 参考手册 3.4
网友评论