郑淑梅

1 中国地震局第一监测中心，天津市耐火路7号，300180

线性度是地震数据采集器的关键性能指标之一，可反映地震数据采集器输出信号是否失真。评价线性度的常用方法有最小二乘法[1-3]和端基直线法[4]，地震数据采集器线性度计算多采用最小二乘法。

为评估最小二乘法和端基直线法对地震数据采集器线性度测量不确定度的影响，本文从理论上给出2种方法评价线性度的区别，通过实例展示不确定度的计算过程。最后通过对比多种型号的地震数采测量结果，分析2种方法的差别及不同输入点个数和位置对线性度不确定度的影响。

1 最小二乘法和端基直线法的线性度理论

对于线性系统，输入量x和输出量y之间的关系为：

f(x)=kx+b

(1)

式中，k为直线斜率，b为直线截距。

用xi、yi(i=1,2,…,n)分别表示输入值和测量值，则在信号点xi上测量值yi与线性回归值f(xi)的残差可表示为：

Δyi=yi-(kxi+b)

(2)

根据线性度定义，最大残差可反映线性度大小。线性度公式为：

(3)

图1为最小二乘法和端基直线法拟合示意图。可以看出，二者计算线性度的区别在于k和b的估计值不同导致残差Δy不同，残差的正负和大小可表征测量值在拟合直线两侧的分散程度。按照最小二乘法原理，当残差平方和最小时，k和b估计值最佳；端基直线法中，k和b只与首尾点有关。

图1 线性拟合示意图Fig.1 Schematic diagram of linear fitting

(4)

最小二乘法

(5)

(6)

2 测量不确定度模型

考虑到被测数据采集器的重复性和分辨力存在重复，合成不确定度时将二者中的最小值舍去。

根据不确定度理论公式[5]与文献[1,4]可得，线性度L的合成标准不确定度为：

(7)

其有效自由度veff=m-3。

3 测量数据处理及分析

3.1 测量数据获取及处理

采用标准源法，利用5730A输出标准直流电压信号作为地震数据采集器的输入。选用EDAS-27HR型号地震数据采集器，测量范围设为(-10 V, 10 V)，采样率设为100 Hz，其A/D位数为28 bits，分辨率为0.149 μV/count。根据地震监测专业设备入网测试规程[6]选择测试点，从-10～10 V每隔1 V选取1个点，首尾选取-9.5 V和9.5 V，共21个测试点。每个测试点测量2 min，截取20 s的数据计算其平均值，并将平均值折合到对应的输入端作为最终的测量结果，由5730A标准仪器引入的不确定度分量u(Ei)采用B类评定，Ei的允许误差极限由仪器说明书获得。如参照5730A技术指标，测量2 V直流电压时的扩展不确定度为5 ppm×Vout+0.7 μV，k=2，则其区间半宽为：

5.35×10-6V

在日常测量中，一般取2 000(20 s)个点的平均值作为测量结果，测量结果的实验方差由式(4)获得。最终的不确定度分量以及灵敏系数见表1。

对于最小二乘法，由式(2)和式(3)可得线性度估计值L=4.74×10-6，j=1。将表1中灵敏系数代入式(7)，计算得到最小二乘法线性度的合成标准不确定度为7.8×10-7V。

对于端基直线法，由式(2)和式(3)可得线性度估计值L=9.1×10-7，j=12。将表1中灵敏系数代入式(7)，计算得到端基直线法线性度的合成标准不确定度为6.5×10-7V，有效自由度均为veff=21-3=18。

表1 不确定度分量及灵敏系数Tab.1 Uncertainty components and sensitivity coefficients

3.2 测量数据分析

3.2.1 不同方法对计算结果的影响

分别采用最小二乘法和端基直线法对不同分辨力(数据位数)、不同型号的地震数据采集器的线性度结果进行不确定度评定，结果见表2。

从表2可以看出，最小二乘法与端基直线法线性度不确定度评定结果在同一数量级，与分辨力无关。因此，对于目前常用的地震数据采集器而言，无论采用最小二乘法还是端基直线法计算的线性度，其结果的离散度均很小。端基直线法线性度不确定度的评定结果较最小二乘法稳定，这主要与不确定度计算公式有关。

表2 最小二乘法与端基直线法线性度的标准不确定度Tab.2 Standard uncertainty of linearity of the least squares method and the terminal-based straight line method

3.2.2 数据点数和位置对计算结果的影响

为研究数据点数和位置对不确定度评定结果的影响，分别对19点、17点、15点、13点、11点、9点和7点数据在不同位置组合下线性度的不确定度进行评定。为不失一般性，仍然采用对称的选点原则。如图2所示，以19点数据为例，即在原21点数据的基础上从大到小依次去掉对称的一组数据，因此有9个位置组合。同理，17点、15点、13点、11点、9点和7点分别有8个、7个、6个、5个、4个和3个位置组合。2种方法在不同数据点及不同位置组合下的不确定度评定结果见图3。

图2 不同数据点数的位置组合示意图Fig.2 Schematic diagram of the location combination of different data points

图3 不同数据点数和位置对不确定度的影响Fig.3 The influence of different data points and locations on uncertainty

从图3可以看出，最小二乘法线性度不确定度随数据点数的增加而增加，其主要原因是不确定度计算公式为求和过程，数据点数越多其不确定度值越大。相对于最小二乘法，端基直线法线性度不确定度评定结果受数据点数影响较小，虽然不确定度的理论公式相同，但从灵敏系数可以看出，实际有效点数为3个值累加，因此其受数据点数的影响较小。

对于最小二乘法线性度不确定度而言，数据点数越多，对数据位置越不敏感；选取的数据越靠两侧，不确定度越小，说明靠边的数据对不确定度的影响较大。端基直线法线性度不确定度均在末位位置组合处发生突变，这是因为除末位位置组合外，其他位置组合的最大残差均在1 V点处出现，因此不确定度未发生任何变化，而末位位置组合已剔除1 V点，最大残差点发生变化，因此出现不确定度不同程度的突变。

4 结 语

本文对地震数据采集器线性度不确定度进行评定，分析目前常用的最小二乘法和端基直线法线性度不确定度的差别。多种型号的地震数据采集器线性度不确定度评定结果表明，2种方法线性度不确定度评定结果在同一数量级，可以满足地震数据采集器线性度不确定度评定的要求，但最小二乘法线性度不确定度受测试点数的影响较大，端基直线法线性度受测试点数的影响较小。