0 引言

地震计传递函数是描述地震计输入、输出关系的表达式，是地震观测、地震学研究的第一手关键资料。早期机械地震计可以通过测试地震计自振周期和阻尼获得二阶传递函数。电子反馈型地震计的传递函数受电子线路影响，由于机械部分和电子器件性能离散性，通过正向推导获得准确的传递函数存在困难。地震计生产厂家一般采用实际测试方法确定地震计传递函数，主要采用振动台方法、电标定方法等。

振动台方法是对地震计整体施加振动，与地震计工作模式一致，是一种客观标定方式。目前该方法主要用于地震计幅频特性测试，由于缺少同步振动台记录波形，不能测定相频特性。

不仅如此，根据不同的媒介环境及时代背景因素，新闻叙事者的内涵也会随之发生变化。在传统新闻业中，记者、编辑及其所属机构占据于整个新闻叙事的主导地位，它们便是叙事者。而在大数据背景下，依托互联网及计算机技术产生的数据新闻中，叙事者则变得更为复杂。

电标定方法给地震计磁缸内的标定线圈输入一个电流（Peterson J et al，1980；王广福，1986； Rodgers P W et al，1995；李海亮，2000； Bormann P，2002），产生驱动力驱动地震计摆锤运动以等效成地面运动。目前电标定可以远程控制，标定时地震计放置在稳定地面，保持环境安静，干扰较小，可以更加准确地测定地震计响应。电标定方法根据信号不同主要可分为阶跃标定、伪随机码标定和正弦波标定。

河流的生态水利工程就是为了能使河流保持原有的自然属性，同时满足人类社会经济发展需要及河流生态系统本身健康发展的需求[8]。

阶跃标定信号自身频谱特性决定其低频信号强，高频成分随频率的平方衰减，有利于测定地震计低频特性。事实上，基于地震计阶跃标定响应特性推导的时域波形方程和响应频谱公式，准确给出了地震计自振周期和阻尼系数的关系。基于此，阶跃标定主要用于确定地震计低频特性。

（3）名义电流和实际电流误差对测试的影响。表8给出本实验中名义电流与实际电流的偏差，其平均相对偏差为-1.533%。在低频段，偏差与电流大小有关，主要是因为电流档位变化，改用不同电路，电路之间存在差异。在电流达到满刻度值后，电流明显有规律的变小，主要是标定信号发生器D/A转换采样率固定，当频率升高时，每个周期内样点减少所致。

在取得地震计的频率特性和相位特性后，采用拟合方法得到地震计传递函数，是一种通行方法（Bolduc P M et al，1972；王广福，1986； Kim W Y et al，1996）。正弦波标定使用广泛，不同频率可以设定相应激励电流，保证在测试点具有足够信噪比的信号用于确定频率特性。传统正弦标定方法限于测定地震计的幅频特性，而不测定相频特性。本文提出改进的技术措施，同时获取幅频频特性和相频特性，并利用阶跃标定特点，准确测定地震计自振周期和阻尼，确定地震计低频特性；在拟合地震计传递函数时，根据地震计频率特性的特点构造传递函数的形式，并约束低频传递函数参数，保证结果的合理性和可靠性。

1 测试原理与方法

1.1 反馈式地震计传递函数表述

任意一个系统的传递函数均可表示一个滤波器，其形式是2个多项式的商。本文定义反馈式地震计的传递函数为在频域经过归一化的地震计响应，由2个部分组成，其中低频部分取决于等效机械摆的自振角频率和阻尼系数，称为低频传递函数，记为H0(s)；另一部分主要影响地震计的高频特性，称为高频传递函数，记为H1(s)。地震计的总体频率特性、低频和高频特性可分别用式（1）、（2）、（3）来描述。

其中：H0(s)为等效固有周期机械摆的二阶传递函数，即低频传递函数，s为Laplace算子，D0为等效机械摆的阻尼系数，ω0为等效机械摆自振周期，(an，an-1，...，a1，a0)、(bm，bm-1，...，b1，b0)分别为其分子分母多项式的系数。式（2）中的H0(s)由阶跃标定获得的自振周期和阻尼系数确定。式（3）中的H1(s)描述地震计传递函数的高频段特性，取决于地震计的电路特性。(an，an-1，...，a1，a0)、(bm，bm-1，...，b1，b0)系数用最优化拟合确定，是本文研究重点。

研究发现石渠县南部强度负向人为干扰最为集中；中度负向人为干扰区域集中在甘孜县、色达县以及壤塘县；轻度负向人为干扰区域在整个研究区分布较为均匀；甘孜县与色达县无明显人为干扰区域占比较大；轻度正向人为干扰区域则主要分布在石渠县、若尔盖县以及红原县；中度正向人为干扰区域主要集中在石渠县，其余各县均有分布；强度正向人为干扰区域，仅在甘孜最南端有分布。

再者，微课的实行满足了教育公平。由于高职院校所在地不同，导致一些学校在师资力量、教学设备方面都存在差异性，学生受教程度各有不同。微课的出现最开始主要是因为人们对一线教师授课方式共享的渴望，主要作用是可以实现优质的教学资源和学习资源共享。随着时代不断发展，微课逐渐被引入教学，对高职院校的学生来说，大大减轻了一些中低收入家庭的学生对教育费用的负担。同时，高职学生还可以利用微课提供的学习情境，将学习的知识学以致用。

1.2 地震计自振周期与阻尼系数确定

阶跃标定是确定地震计自振周期和阻尼的有效可靠的方式（李海亮，2000）。地震计阶跃响应的记录波形在时间域和频率域内可分别由式（4）、（5）描述。

阶跃响应由K、ω0、D0参数唯一描述，其中的K为阶跃等效加速度a与地震计速度灵敏度S的乘积，ω0为固有角频率，D0为阻尼系数，从式（4）可见，当D0≤1时，r(t)的定义在实数域内将失去意义，式（4）只在D0 ＜1时成立。

使用实际记录的地震计阶跃响应波形计算其响应频谱，通过式（2），采用最小二乘法在频域内拟合，可以方便、可靠地求得K、ω0、D0，其中ω0=2π/T0。本文关注重点是T0、D0，它们决定了地震计的低频特性。

1.3 地震计响应特性精确测定及H1(s)求取

地震计传递函数包括幅频信息和相频信息。在普通的正弦标定中，无参考信号记录，只能测定幅频特性。利用本文方法进行正弦标定时，记录地震信号输出时，同步记录标定信号在测流电阻上产生的电压，处理正弦波标定信号可以同时获得地震计的振幅信息和相移信息。

为验证所提出方法的可行性，首先测试标准地震计的传递函数。选定北京港震仪器设备有限公司生产的BBVS-120甚宽带地震计为标准地震计，其中-3 dB频带为120 s—50 Hz，标称灵敏度为2 000 V·s/m。该地震计灵敏度在振动台经过精确测试，并校正其标定灵敏度，数值分别为 10.986 70、10.785 00、10.509 10。

在使用式（6）进行正弦波标定时，需要设定正弦波的频率和电流参数。电流参数可能受到恒流源精度和分流电路的影响，但频率只与采集系统的时钟频率有关，是高度精确的，因此可以把频率当成已知量来计算正弦波的振幅和周期。本文采用最小二乘法拟合求得精确振幅和相位值。在拟合中设定振幅精度目标为10-12，保证所求振幅和相位高度精确，满足正弦波标定数据处理的需要。地震计传递函数测定数据处理步骤如下。

向地震计发出如表2所示标定参数的正弦波标定信号，信号参数在通带内尽量按5 mm/s的速度设计，在通带外适当增加，但在高频部分受采集系统标定信号幅度上限20 000 µA的限制，速度小于5 mm/s。在标定过程中，数据采集器自动波形使用记录波形，计算得到地震计3个分向输出和测流电压信号在各同频率点的幅度和相位，使用地震计输出电压和测流电压，并结合标定灵敏度，计算地震计的测试幅频响应和相频响应（图3），实际测得的幅频响应值和相频响应值分别见表3和表4。

权筝把手镯又塞给何东：“又不是你媳妇我留着干嘛？”她其实是在试探何东，想听何东说，你先留着，我这劲儿一过，咱们过几天就登记去。这也是她今天来爷爷家的目的，不把事儿弄僵，一切就还有挽回的余地。

（2）计算地震计输出各频点的3个方向的相位响应。电流流过标定线圈产生驱动力推定摆锤运动，力与加速度成正比，地震计输出与速度量成比例的电压。因此，地震计响应在指定频率点相位响应为

其中：φr为相位响应，φs为地震计输出信号相位，φI为标定电流相位。

（3）计算实际标定电流。用获得的测流电压除以测流电阻值即得到实际标定电流。

（4）计算地震计幅频特性。使用已有地震计输出3个分向各频率的振幅、标定灵敏度Sc和实测标定电流值，按式（8）计算地震计幅频特性。

（5）计算地震计响应函数 R(f)。将式（8）计算的地震计幅频特性带入式（9），可得

（6）计算所得地震计传递函数在各频点的实测值Hm(f )。

（7）拟合确定H1(s)。采用最优化方法拟合求H1(s)，其目标函数为

其中，W(i)为对应频点的权重。在实际计算中，可对频带内数据优先保证拟合水平指定较高的权重，对频带外数据指定较小的权重；如无特定要求，可指定相同的非零权重。D0、ω0阶跃标定已确定，拟合的目标是确定 A∈(an，an-1，...，a1，a0)、B∈(bm，bm-1，...，b1，b0)。

图1 精确测定地震计传递函数电路示意 Fig.1 The schematic diagram of measuring seismometer transfer function

1.4 测试电路设计

与普通地震计电标定不同，本文提出的测试方法引入测流电阻R，用于同步测试标定回路中的电流。实验中采用的电路见图1，由标准EDAS-24GN地震数据采集器输出的电流标定信号依次流过测流电阻、地震计UD、EW、NS向的标定线圈后返回。地震计3个通道的输出信号输出到标准地震计数据采集器的A口3个主通道，测流电阻两端电压输出到地震数据采集器B口的某一通道（实测时选用第2通道）。当正弦标定电流流过标定线圈时，产生一个力作用在地震计摆锤上，地震计的响应波形由地震数据采集器记录下来，并同步记录测流电阻两端的电压信号。测流电阻两端的电压信号有2个作用，即精确计算实际流过标定回路的电流及作为测试地震计相位特性的参考信号。

2 测试实验及数据处理

任意一段正弦波形均可表示为

2.1 实验仪器

测试使用的标准数据采集器为一台6通道EDAS-24GN24高分地震数据采集器，其采集灵敏度经过标准信源校正，误差小于1/10 000。标准地震数据采集器共有6个信号输入通道，分别位于A、B端口。地震计及测试信号连接如图1所示，2个端口的输入范围均设置为±20 V，采样率均设置为500 sps，采集器采用28位字长编码，电压分辨率在该设置下为149 nV。来自地震计和测流电阻端的电压信号经相同处理并记录下来，地震计输出信号与测流信号的差别与采集记录下来的信号差别相同。EDAS-24GN的6个通道同步采集数据，进行同步数据分析即可取得准确的仪器振幅和相位响应。

如果生产工艺中对电机降速时间没有要求，可以手动更改电机减速时间参数，在保证不引起超电压的前提条件下，避免出现负载减速太快情况。如果负载减速时间收到生产工艺流程的制约，为预防变频器在减速时间内出现过电压跳停，可降低变频器的频率值，通过降低频率控制降低设备运行速度。

2.2 阶跃标定确定周期和阻尼

向地震计发出阶跃响应，通过阶跃波形拟合获得地震计的周期和阻尼。阶跃波形上升侧频域、时间域的拟合见图2，下降侧拟合与上升侧相似。阶跃拟合获得的地震计等校机械摆参数见表1，由拟合结果可见，上升侧和下降侧拟合获得的参数高度一致。

图2 地震计阶跃标定波形拟合结果 Fig.2 The waveform fitting result of step calibration of the seismometer

表1 阶跃标定拟合得到等效机械摆参数 Table 1 The equivalent mechanical pendulum parameters from step calibration

参数 上升侧阶跃 下降侧阶跃UD EW NS UD EW NS灵敏度/V·s·m-1 2 016.8 2 027.4 1 930.9 2 017.3 2 027.2 1 930.5自振周期/s 120.5 120.4 120.1 120.5 120.4 120.1阻尼系数 0.703 34 0.703 59 0.701 22 0.703 14 0.703 40 0.700 80

2.3 正弦波标定获得地震计幅频特性和相频特性

（1）对每一个频率、每一个通道和测量电压的波形取出一段同步数据，使用正弦波拟合方法，得到同一段地震计3个分向和测流电压的振幅和相位。

（2）拟合结果的多解性及最终结果选取。根据选择的参数不同，同一频率特性可以有多个传递函数表达式，这些表达式在形式上是有差别的，是实际测量数据拟合得到的结果，只要与实际数据吻合，均可作为地震计传递函数。在实际工作中应当选用与实际数据拟合残差小于一定水平范围内最短数据长度的拟合结果。

表2 正弦波标定参数设计 Table 2 Sine wave calibration parameters

注：T/f表示频率或周期，以负数表示周期，绝对值为周期值，单位：s；正数表示频率，单位：Hz。 下文相同

重复周期数 T/f/ (s/Hz)名义电流/µA 设计激励速度/mm ·s-1重复周期数 T/f/ (s/Hz)名义电流/µA 设计激励速度/mm ·s-1 2-1 000 5 7.957 7 60 2 6 283 4.999 9 2-720 6 6.875 5 150 5 15 707 4.999 9 2-360 13 7.448 5 300 10 19 999 3.183 0 2-240 20 7.639 4 600 20 19 999 1.591 5 2-120 41 7.830 4 900 30 19 999 1.061 0 3-60 83 7.925 9 1 200 40 19 999 0.795 7 4-30 167 7.973 7 1 500 50 19 999 0.636 6 5-20 251 7.989 6 1 800 60 19 999 0.530 5 6-10 314 4.997 5 2 100 70 19 999 0.454 7 6-5 628 4.997 5 2 400 80 19 999 0.397 9 15 -2 1 570 4.997 5 2 700 90 19 999 0.353 7 30 1 3 141 4.999 1 3 000 100 19 999 0.318 3

表3 实际测得的地震计幅频响应值 Table 3 The actual measured amplitude frequency response of the seismometer

T/f/(s/Hz)地震计输出电压/V 实测标定 电流/μA幅频响应值/V ·s·m-1 UD EW NS UD EW NS-1000 0.123 33 0.122 94 0.116 89 4.913 60 28.71 29.15 28.44 -720 0.203 97 0.202 36 0.194 29 5.886 98 55.04 55.63 54.81 -360 0.874 68 0.868 33 0.835 80 12.743 5 218.07 220.54 217.85 -240 2.006 44 1.992 66 1.916 28 19.745 8 484.27 489.93 483.52 -120 5.790 55 5.748 78 5.551 88 40.416 5 1 365.60 1 381.10 1 368.81 -60 8.208 09 8.152 41 7.878 22 81.891 4 1 910.71 1 933.24 1 917.27 -30 8.502 21 8.445 14 8.153 38 164.732 1 967.77 1 991.12 1 972.80 -20 8.504 79 8.447 84 8.154 00 247.276 1 966.96 1 990.33 1 971.53 -10 5.322 65 5.286 91 5.102 28 309.361 1 967.90 1 991.25 1 972.16 -5 5.319 67 5.283 81 5.098 72 619.055 1 965.75 1 989.01 1 969.73 -2 5.336 78 5.298 90 5.111 59 1 548.39 1 971.12 1 993.72 1 973.74 1 5.369 38 5.325 51 5.132 20 3 091.02 1 986.85 2 007.47 1 985.39 2 5.463 37 5.406 64 5.196 52 6 185.91 2 020.37 2 036.78 2 009.02 5 5.451 90 5.409 06 5.180 21 15 471.4 2 015.26 2 036.82 2 001.85 10 3.390 72 3.392 88 3.259 30 19 703.0 1 968.35 2 006.44 1 978.05 20 1.767 06 1.729 56 1.701 06 19 711.6 2 050.70 2 044.71 2 063.82 30 1.159 59 1.112 03 1.139 61 19 685.8 2 021.22 1 974.58 2 076.67 40 0.747 18 0.721 66 0.761 57 19 649.5 1 739.71 1 711.73 1 853.79 50 0.473 50 0.470 32 0.495 22 19 614.9 1 380.52 1 396.89 1 509.46 60 0.326 30 0.322 53 0.353 75 19 590.4 1 143.07 1 150.98 1 295.54 70 0.206 56 0.217 63 0.238 60 19 570.5 845.07 907.00 1 020.47 80 0.137 16 0.148 68 0.160 08 19 520.8 642.93 709.96 784.45 90 0.086 09 0.147 17 0.150 58 19 448.7 455.66 793.54 833.23 100 0.054 17 0.069 91 0.128 75 19 393.6 319.49 420.02 793.83

表4 实际测得的地震计相频响应值 Table 4 The actual measured phase frequency response of the seismometer

T/f/(s/Hz)实测地震计相位响应/rad 实测标定 电流/μA实测地震计相位响应/rad UD EW NS UD EW NS-1 000 2.970 356 2.968 107 2.9715 25 2 -0.063 51 -0.047 69 -0.045 79-720 2.903 391 2.904 945 2.903 267 5 -0.220 74 -0.173 85 -0.159 63-360 2.653 914 2.653 044 2.657 528 10 -0.417 86 -0.336 64 -0.301 20-240 2.386 593 2.385 801 2.390 558 20 -0.840 56 -0.687 79 -0.615 24-120 1.549 783 1.549 258 1.552 565 30 -1.378 40 -1.114 03 -1.025 66-60 0.740 88 0.741 099 0.741 412 40 -1.918 64 -1.531 17 -1.456 51-30 0.358 466 0.358 487 0.358 309 50 -2.362 35 -1.871 90 -1.791 56-20 0.236 734 0.236 643 0.236 446 60 -2.814 13 -2.170 19 -2.106 65-10 0.116 068 0.115 745 0.115 597 70 -3.227 06 -2.431 84 -2.386 39-5 0.056 371 0.055 687 0.055 529 80 -3.591 67 -2.503 47 -2.561 57-2 0.017 362 0.015 685 0.015 453 90 -3.926 57 -2.902 90 -2.814 15 1 -0.004 26 -0.007 1 -0.007 16 100 -4.132 96 -3.382 88 -3.121 67

图3 实测标准地震计频率特性 (a)地震计幅频特性；(b)地震计相频特性 Fig.3 The actual measured frequency characteristics of standard seismometer

2.4 拟合求地震计传递函数

坡度空间数据由DEM数据在ArcGIS 9.3中，通过Slope功能生成。曼宁系数空间数据利用ArcGIS 9.3将查阅文献获得的曼宁系数属性数据（表1）与土地覆盖类型空间数据相关联生成。土壤饱和导水率与土壤储水能力，通过结合土壤类型组成及其土壤剖面等属性数据，借助于土壤水分运动参数模型RETC推导获得，空间数据在ArcGIS 9.3下通过建立土壤水分运动参数与土壤类型空间数据之间的关联生成。

表5 标准地震计UD向传递函数拟合结果 Table 5 The fitted result of the UD component of the seismometer transfer function

拟合系数零极点数Nz=3，Np =6 Nz=4，Np =7 Nz=5，Np=8 A 1 1 1 1 587.857 448 668 78 1 588.702 902 536 54 1 647.690 311 319 30 542 438.583 793 854 542 707.687 302 222 574 489.622 336 054 108 675 514.140 093 108 706 728.312 693 119 200 934.216 165 413 920 092.683 199 407 533 343.868 138 1 794 053 227.686 39-18 428 207.021 333 3 -39 407 907 373.359 8-186 708 332 709.790 B 107 080 966.507 916 107 170 107.886 847 109 245 360.053 649 407 918 946.532 043 402 305 252.224 082 2 015 721 360.077 90-17 789 371.902 507 2 -39 236 882 637.169 3-184 069 429 964.200残差 0.013 855 0.013 518 0.013 189

表6 标准地震计EW向传递函数拟合结果 Table 6 The fitted result of the EW component of the seismometer transfer function

拟合系数零极点数Nz =3，Np =6 Nz=4，Np =7 Nz =5，Np =8 A 1 1 1 1 498.689 105 893 48 1 232.246 518 234 47 1 499.763 192 190 55 557 855.282 037 954 932 333.152 148 272 559 768.018 218 627 107 097 923.765 168 264 074 576.270 363 107 721 394.256 087 659 740 219.695 601 39 402 688 659.045 2 779 108 504.792 030 112 581 804 701.395 286 491 664.840 142-14 799 496.641 555 8 B 105 463 984.554 566 63764062.0550742 105 515 797.520 385 651 887 013.348 565 39 007 480 880.486 2 772 019 916.409 103 110 986 887 496.447 283 938 545.640 724-14 278 864.590 389 9残差 0.016 648 0.016 491 0.016 497

表7 标准地震计NS向传递函数拟合结果 Table 7 The fitted result of the NS component of the seismometer transfer function

拟合系数零极点数Nz=3，Np =6 Nz =4，Np=7 Nz=5，Np=8 A 1 1 1 1 642.974 777 532 25 1 655.440 785 993 71 1 671.457 496 423 85 606 675.498 387 427 610 955.025 177 009 605 549.630 542 702 129 492 940.176 938 130 525 305.852 803 128 794 261.650 657 772 209 380.070 837 722 541 333.878 945 61 556 478.901 580 6-49 914 645.046 524 6 -12 415 574.962 803 8-248 016 599.777 630 B 127 719 713.159 972 128 997 760.622 664 129 775 980.995 086 764 319 935.286 854 716 756 127.299 418 60 055 004.663 173 8-48 612 286.841 825 7 -8 970 368.983 448 19-245 388 893.059 357残差 0.019 351 0.020 430 0.020 043

图4 幅频特性和相频特性传递函数拟合结果 Fig.4 The fitted frequency characteristics of the seismometer transfer function

3 结果分析

（1）地震计传递函数低频、高频部分的关系及各自作用。图5清楚显示地震计传递函数可以分为低频段和高频段2个部分，其中低频段由其等效机械摆的周期和阻尼决定，高频段通过拟合求得。低频段和高频段综合在一起形成地震计的整体传递函数。

按前文所述，标准传递函数表达式见式（1）、式（2）、式（3）所示，拟合传递函数的幅频特性和相频特性见图4，标准地震计UD向、EW向和NS向传递函数拟合系数及残差分别见表5—表7。

陈主任比吕凌子先到一步。两人找到经销老板，老板姓莫，是位四十岁左右的女姓，身材有点富态，但穿着时髦，脸上白得很不自然。

观察组:腹腔镜手术。患者气管插管全身麻醉,于脐轮上缘或下缘做横切口约1.5cm,气腹针穿刺成功,建立气腹并维持12-14mm Hg气腹压,常规三点作穿刺点、置镜点及探查点。腹腔镜下探查肿瘤位置、体积和性质,观察卵巢周围与盆腔关系,充分暴露肿瘤表面,在卵巢最薄弱处,电凝切开卵巢皮质,卵巢皮质及卵巢壁间间钝性或锐性分离,完整剥除肿瘤,切除附件。2-0可吸收缝线缝合剩余正常的卵巢及包膜组织。待腹腔无出血后,关闭切口。

伪随机码标定方法（Berg E et al，1976；Berger J et al，1979；林湛等，2015）是向地震计发送精心设计的伪随机二进制码信号，通过处理地震计响应来获得地震计幅频特性。本质上二进制码信号是一系列不同宽度的阶跃信号组合，其能量在频率分布上具有不均匀性，且高频信号强度不足，通过谱除法获得频率特性有较大波动，信噪比变化较大，经过拟合也可求得传递函数，置信程度有一定不足。

图5 拟合结果显示的低频段、高频段和综合频率特性 (a)低频段传递函数的幅频特性和相频特性；(b)传递函数高频段的幅频特性和相频特性；(c)综合幅频特性和相频特性 Fig.5 The fitted result of the low frequency band,high frequency band and full frequency band

表8 名义标定电流和实际标定电流偏差 Table 8 The deviation of nominal calibration current and actual calibration current

T/f/ (s/Hz)名义电流/μA 实际电流/μA 相对误差（%）T/f/ (s/Hz)名义电流/μA 实际电流/μA 相对误差(%)-1000 5 4.913 60 -1.728 01 2 6 283 6 185.91 -1.545 35-720 6 5.886 98 -1.883 6 5 15 707 15 471.4 -1.499 93-360 13 12.743 5 -1.972 76 10 19 999 19 703.0 -1.480 22-240 20 19.745 8 -1.271 19 20 19 999 19 711.6 -1.436 87-120 41 40.416 5 -1.423 15 30 19 999 19 685.8 -1.566 18-60 83 81.891 4 -1.335 63 40 19 999 19 649.5 -1.747 81-30 167 164.732 -1.357 9 50 19 999 19 614.9 -1.920 53-20 251 247.276 -1.483 84 60 19 999 19 590.4 -2.043 05-10 314 309.361 -1.477 23 70 19 999 19 570.5 -2.142 48-5 628 619.055 -1.424 32 80 19 999 19 520.8 -2.391 03-2 1 570 1 548.39 -1.376 54 90 19 999 19 448.7 -2.751 56 1 3 141 3 091.02 -1.591 19 100 19 999 19 393.6 -3.026 94

表9 垂直向4零点、7极点模式下的零极点 Table 9 The zero-pole parameters of the UD component

?

（4）频率特性的几种形式及转换。地震计传递函数有3种形式：①幅频、相频特性数据列表；②s算子的多项式商；③零极点形式。本文已明确证明可以通过拟合把幅频、相频特性列表数据转换为s算子的多项式商，进一步假定分母多项式和分子多项式为0，求出各自的解，得到传递函数的零极点。表9为利用该方法求得的1组垂直向零极点。

若必要，也可按式（13）把零极点转换成传递函数，进一步代入式（14）计算各频点的复数响应，计算幅频响应和相频响应。

世袭制作为“其父善游”的思想变种，其实并不如看上去那般美好。比如南唐后主李煜，明明可以成为诗词大家，偏偏做了皇帝，结果沦为阶下囚，被赐毒酒而死；再比如赵佶，本来可以发扬光大花鸟画和瘦金体，做个举世无双的大艺术家，偏偏做了宋徽宗，落得个“诸事皆能，独不能为君”的盖棺论定。此二人便是被“其父善游”所害。他们就像那即将被抛下江的小儿，只能被历史的洪流所裹挟，绝无翻身逃脱的命运。只要那个扔孩子的坚持奉“其父善游”为圭臬，世上就永远不缺冤死鬼。

（5）其他类型地震计测试结果。为验证方法的普适性，选取1台北京港震仪器设备有限公司开发的GL-CS60宽频带地震计和1台FSS-3M短周期地震计进行传递函数对比测试。选择3零点、6极点模式，测试结果见表10，可见拟合结果较好。GL-CS60是斜轴悬挂三分向地震计，通过变换后得到UD、EW、NS分向的输出，标定时以斜轴方式输出，得到斜轴地震计传递函数。

表10 GL-CS60宽频带、FSS-3M短周期地震计实测传递函数 Table 10 The actual measured transfer function of GL-CS60 and FSS-3M

参数 GL-CS60 FSS-3M 参数 GL-CS60 FSS-3M A 11 B 583 563 751.932 584 -91 760 221.279 184 4 2 783.947 788 532 72 1 085.774 591 695 27 19 029 361 549.143 4 65 420 175 090.852 4 2 620 150.896 956 19 859 621.969 217 559 自振周期/s 60.430 91 2.069 0 642 178 684.054 033 333 371 381.594 190 阻尼系数 0.707 351 0.694 1 19 139 236 259.644 3 64 415 601 944.050 0 残差 0.004 915 0.040 189

4 结论

综上所述，可知本文提出的使用阶跃标定和正弦标定测定地震计传递函数的方法，具有以下特点。

（1）本文提出联合使用阶跃标定和正弦波标定测定地震计传递函数的方法，只需要1台6通道地震计采集器和1条经过修改的摆线，采用数据采集器固有功能即可完成地震计传递函数测试，具有操作简单、容易实现、结果置信度、精度高的特点。

（2）本文提出的地震计传递函数测量方法与伪随机码标定方法相比，具有信号频点信噪比高、结果可靠的特点，拟合结果准确，直接反映地震计的频率特性高。

（3）本文提出的方法采用的目标函数和传递函数形式考虑到地震计传递函数的特点，拟合结果物理意义明确，测试结果说明低频段受等效机械摆自振周期和阻尼控制，高频段取决于地震计滤波器特性。

林业资源管理中机构配置不相符主要体现在基层林业站中。众所周知，基层林业站主要负责的工作除管理林木、爱护林木和加大绿化之外，还要宣传和推广相关的先进技术和林业政策等。由于我国基层林业站的人员有限，相关单位的工作人员承受着很大的工作压力，要高效完成林业资源管理工作，难度较大。除此之外，基层工作站工作人员短缺，导致一人多职的现象较为严重。身兼多职在一定程度上降低了管理人员的工作效率，容易在工作中出现纰漏。

（4）本文提出的方法具有普适性，可以用于甚宽带地震计、宽带地震计、短周期地震计的传递函数测试。

参考文献

李海亮.地震计阻尼和自振频率的频域测定[J].地震地磁观测与研究，2000，21（2）：7-12.

林湛，薛兵，朱小毅，等.地震计的改进型伪随机二进制标定方法[J].地球物理学进展，2015，30（3）：1 151-1 158.

王广福.地震观测系统的标定与检查[M].北京：地震出版社，1986：165-170.

Berg E and Chesley D M.Automated High-Procision Amplitude and Phase Calibration of Seismic System[J].BSSA,1976,66（4）:1 413-1 424.

Berger J,Agnew D C,Parker R L,Farrell W E.Seismic system calibration: 2.cross-spectral calibration using random binary signal [J].BSSA,1979,69（1）: 271-288.

Bolduc P M,Ellis R M and Russell R D.Determination of the seismograph phase response from the amplitude response[J].BSSA,1972,62（6）:1 665-1 672.

Bormann P.New Manual of Seismological Observatory Practice（NMSOP-2）[M].Potsdam: GeoForschungs Zentrum,2002: 34-40.

Kim W Y and Ekström G.Instrument Responses of Digital Seismographs at Borovoye,Kazakhstan,by inversion of transient calibration pulses [J].BSSA,1996,86（1A）: 191-203.

Peterson J,Hutt C R and Holcomb L G.Test and Calibration of the Seismic Research Observatory[M].USGS Open-File Report,1980: 80-187.

Rodgers P W,Martin A,Haver M R,et al.Signal-coil Calibration of Electromagnetic Seismometers[J].BSSA,1995,85（3）: 845-850.