引言

目前，我国地电前兆台网广泛使用的地电阻率观测方法是直流观测方法，其观测装置采用对称四电极装置(Zhao et al，2011)．随着我国国民经济和社会的发展，地电阻率观测受到越来越多的电磁干扰，特别是城市轨道交通(地铁、 轻轨)的干扰，其干扰频段主要集中在直流附近，且影响范围可达数十千米，严重影响了城市附近地电阻率观测台站的观测(张宇等，2014)．张世中等(2013)在北京西集、 天津青光两个台站进行场地干扰测试，发现受地铁干扰的影响，两台站每天有2/3以上的小时观测值的相对方差大于0.3%，有时甚至超过3%．

采用交流供电的地电阻率方法是消除此类干扰的一个有效方法．该方法采用低频电流作为激励源，接收低频电流经大地传输后的响应信号，可以采用频率域测量技术进行选频接收(罗维斌等，2004)，以获得较高的信噪比．

20世纪80年代，中国地震局曾引进美国研制的MARK型低频交流地电仪，在河北开滦马家沟地震台开展地电阻率交流法的监测试验，其结果证明了地电阻率交流观测方法的可行性(桂燮泰等，1988；马希融，1989)．

本文将首先阐述地电阻率交流观测方法的原理；然后基于地电阻率的实际测量过程，采用数值模拟方法分析不同参数条件下，数字滤波器法、 频谱分析法和相关检测法等3种信号检测方法的检测效果差异；最后根据仿真结果给出适用于地电阻率交流观测方法的信号检测方法，以期为地电阻率交流观测系统的研制提供理论参考．

1.   理论基础

1.1   地电阻率交流观测原理

在0.01—1000 Hz频段，交流电场分布可近似地遵循欧姆定律，电阻率交流测量方法的原理与直流测量方法相同(桂燮泰等，1988；马希融，1989；张国民等，2001)．低频交流电法测量示意图如图 1所示．

图  1  地电阻率交流观测方法工作原理示意图

A和B为供电电极，M和N为测量电极

Figure  1.  Schematic diagram of principle for AC geo-resistivity method

A and B are current-emitting electrodes, and M and N are measuring electrodes

下载: 全尺寸图片 幻灯片

地电阻率交流测量方法采用低频电流作为激励源，通过供电电极A和B向大地供入低频电流I_AB，同时通过测量电极M和N测量MN之间的同频率响应信号U_MN，通过计算得出地电阻率ρ的具体表达式为(张国民等，2001；刘义国等，2010)

式中K为装置系数，其计算公式为

1.2   单频信号检测法

地电阻率交流观测方法采用低频电流源发射单一频率的正弦信号，利用高采样率数据采集系统同步测量人工供电电流和人工供电电位差. 要得到高精度的地电阻率观测数据，就意味着从含有噪声的测量数据中，准确检测出已知频率信号的振幅．

目前单一频率信号振幅检测方法很多，通常采用的有数字滤波器法、 频谱分析法和相关检测法等3种方法．

1.2.1   数字滤波器法

数字滤波是信号处理的一种基本而重要的技术，利用数字滤波器的频率选择特性，可把滤波器的通带设置在能够覆盖有用信号的频谱范围处，使得通带内信号振幅不变，而通带外的噪声受到大幅度衰减，从而提高信噪比，得出已知频率信号的振幅值(高晋占，2007)．

数字滤波器按照设计方法分为有限冲激响应滤波器和无限冲激响应滤波器，两种滤波器理论上都能够实现高精度的滤波效果．与无限冲激响应滤波器相比，有限冲激响应滤波器有突出的优点: 第一，系统总是稳定的，第二，线性相位特性；但其也有明显的缺点，即滤波器的阶数较高，在资源有限的情况下，不易满足实时需求．

常用的有限冲激响应滤波器设计方法有窗函数法、 频率抽取法和最佳一致逼近法，其中窗函数设计法最为常用. 对于各种窗函数来说，在滤波器阶数相同的情况下，凯塞(Kaiser)窗函数具有更好的性能(杜勇等，2012)．因此，本文选择凯塞窗来设计带通数字滤波器进行仿真．

1.2.2   频谱分析法

频谱分析法实现某一特定频率信号振幅的检测是基于任意一个函数x(t)都可以分解为无穷多个不同频率的正弦信号之和的理论，而傅立叶变换就是这样一个工具，它能通过傅里叶分析得到信号的频谱．数字信号经过傅立叶变换后可成为一个N点的复数(N为采样点数)，每一个点对应着一个频率点，这个点的模值就是该频率下振幅谱的相对值．真正的振幅值为原信号相应频率下振幅的N/2倍．因此，通过傅立叶变换，可以很容易地找到指定频率信号的振幅(万永革，2012)．

离散傅立叶变换定义如下：

式中，x(n)为一有限长序列，其长度为N(门爱东等，2009).

1.2.3   相关检测法

相关检测法是基于信号和噪声的统计特性进行信号检测，由于确定性信号不同时刻的取值之间一般都具有较强的相关性，而干扰噪声由于随机性较强，其不同时刻的取值之间相关性一般较差，利用相关性差别可以把确定性信号与干扰信号区别开来(王兰炜等，2004； 张宇等，2014)．因此，从本质上说，相关函数是两个时域信号(或空间域信号)相似性的一种度量(高晋占，2007)．

对于一个已知频率的信号，可在其接收端给出与被检测信号频率相同而相位不同的正弦信号作为参考信号，将参考信号与混有噪声的输入信号进行互相关计算，从而检测出被检测信号的振幅和相位信息． 相关检测法的测量原理如图 2所示． 含有噪声的信号x(t)与参考信号r₁(t)和r₂(t)相乘，分别为y₁(t)和y₂(t)，即

式中，x(t)为被检测信号，n(t)为噪声信号，r₁(t)和r₂(t)为参考信号，A为被测信号振幅，ω₀为被测信号频率，φ为相对相位. 上式中由于信号与噪声是不相关的，因此sin(ω₀t+φ)n(t)和cos(ω₀t+φ)n(t)两项的结果为零，令

这样，Y₁(t)与Y₂(t)经过低通滤波器后，频率为2ω₀的高频项被滤波处理，可以得到 z₁=AB/2cosφ和z₂=AB/2sinφ两个结果项，被检测信号的振幅值则为

图  2  相关检测法的测量原理

x(t)为含有噪声的信号，r₁(t)和r₂(t)为与x(t) 频率相同而相位不同的参考信号

Figure  2.  Measuring principle of correlation detection method

x(t) is the signals including the noise, r₁(t) and r₂(t) are the reference signals whose frequency are same as that of x(t) and whose phase are different from that of x(t)

下载: 全尺寸图片 幻灯片

2.   数值模拟

2.1   输入信号的生成

采用高斯白噪声信号作为数值模拟的噪声信号，为保证模拟仿真中生成的高斯白噪声符合理想条件(均值为0，方差为1，否则会因噪声序列的不均匀而引入误差)，在Matlab软件下生成长度为10万点的固定序列的高斯白噪声n(t)，作为仿真的噪声序列．图 3所示为10万点高斯白噪声序列，其均值为0.0071，方差为1.0051．

图  3  10万点高斯白噪声序列(均值为0.0071，方差为1.0051)

Figure  3.  100000-point Gaussian white noise with mean 0.0071 and variance 1.0051

下载: 全尺寸图片 幻灯片

采用标准正弦信号s(t)作为被检测信号，s(t)=Asin(2πf₀t+φ₀)，式中f₀为信号频率，φ₀为初始相位，A为被检测信号振幅，正弦信号平均功率值P_s=A²/2．根据信噪比计算公式SNR=10lg(P_s/P_n)保持噪声信号不变，由不同信噪比得到不同正弦信号振幅A(表 1)．

表  1  信噪比与正弦信号振幅对应表

Table  1.  SNRs corresponding to sine signal amplitudes

信噪比/dB	正弦信号振幅/mV
-10	0.5
0	1.5
10	4.5
20	15
30	45
40	141.5
50	447

下载: 导出CSV 

| 显示表格

输入信号x(t)由标准正弦信号和噪声两部分组成，即x(t)=s(t)+n(t)． 图 4给出了信噪比为0 dB、 正弦信号频率为1 Hz的输入信号x(t)的波形．

图  4  信噪比为0 dB、信号频率为1 Hz的输入信号x(t)的波形

Figure  4.  Signal with Gaussian white noise(SNR is 0 dB, signal frequency is 1 Hz)

下载: 全尺寸图片 幻灯片

2.2   数值模拟方法说明

数值模拟时，使用系统相对误差ε和相对均方差δ作为评判标准．理论上，相对误差ε越小，说明检测结果与真值越接近，准确度越高；相对均方差δ越小，说明检测结果精密度越好．因此，应该综合考虑相对误差和相对均方差结果，从数字滤波器法、 相关检测法和频谱分析法中选取相对误差和相对均方差结果均较小的方法作为最优方法．

This page contains the following errors:

Below is a rendering of the page up to the first error.

相对误差ε计算公式为

均方差σ的计算公式为

相对均方差结果为均方差与检测均值的比值，计算公式为

数值模拟时，数字滤波器法仿真中设计的带通滤波器性能如表 2所示．

表  2  本文仿真中使用的带通滤波器的性能参数

Table  2.  Bandpass filter performance parameters used in simulation

滤波器序号	中心频率/Hz	采样率/Hz	带宽/Hz	通带波动/dB	阻带衰减/dB	阶数
1	0.1	10	0.02	0.01	80	837
2	0.5	10	0.02	0.01	80	837
3	1.0	10	0.002	0.01	80	541
4	0.1	20	0.002	0.01	80	1068
5	0.5	20	0.002	0.01	80	1068
6	1.0	20	0.002	0.01	80	1014
7	0.1	30	0.002	0.01	80	1602
8	0.5	30	0.002	0.01	80	1602
9	1.0	30	0.002	0.01	80	1521

下载: 导出CSV 

| 显示表格

3.   仿真结果

3.1   不同信噪比的仿真结果比较

通过模拟仿真给出一定采样率和采样点数(采样时间)的条件下，不同信噪比时数字滤波器法、 相关检测法和频谱分析法的检测效果．

假定采样率为10 Hz，采样点数为2000点，在信噪比为-10—50 dB的情况下，表 3—5分别给出了信号频率为0.1，0.5和1 Hz时的仿真结果．图 5给出了不同频率时的相对误差和相对均方差随信噪比变化的结果．可以看出，随着信噪比的增大，3种信号检测方法得到的相对误差和相对均方差均呈下降趋势，表明信噪比越大，信号检测效果越好．当信噪比大于20 dB时，相对误差和相对均方差均趋于平稳，且不同方法的检测结果相近；当信噪比小于20 dB时，频谱分析法的相对误差和相对均方差结果最小，数字滤波器法的结果最大．

表  3  采样率为10 Hz，采样点数为2000，信号频率为0.1 Hz时的仿真结果

Table  3.  Simulation results of three signal detection methods on the condition of different SNR (sampling rate is 10 Hz, number of samples is 2000, signal frequency is 0.1 Hz)

信噪比/dB	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.537	0.037	7.342%	0.044	8.138%
-10	相关检测法	0.5	0.507	0.007	1.445%	0.045	8.837%
	频谱分析法	0.5	0.503	0.003	0.643%	0.034	6.789%
	数字滤波器法	1.5	1.518	0.018	1.179%	0.044	2.930%
0	相关检测法	1.5	1.506	0.006	0.414%	0.045	2.988%
	频谱分析法	1.5	1.503	0.003	0.175%	0.034	2.269%
	数字滤波器法	4.5	4.507	0.007	0.156%	0.046	1.012%
10	相关检测法	4.5	4.506	0.006	0.131%	0.045	1.000%
	频谱分析法	4.5	4.502	0.002	0.054%	0.034	0.757%
	数字滤波器法	15	15.001	0.001	0.009%	0.046	0.305%
20	相关检测法	15	15.006	0.006	0.039%	0.045	0.300%
	频谱分析法	15	15.002	0.002	0.016%	0.034	0.227%
	数字滤波器法	45	44.986	-0.014	0.032%	0.046	0.102%
30	相关检测法	45	45.006	0.006	0.014%	0.045	0.100%
	频谱分析法	45	45.002	0.002	0.005%	0.034	0.076%
	数字滤波器法	141.5	141.435	-0.065	0.046%	0.046	0.032%
40	相关检测法	141.5	141.507	0.007	0.005%	0.045	0.032%
	频谱分析法	141.5	141.502	0.002	0.002%	0.034	0.024%
	数字滤波器法	447	446.776	-0.223	0.050%	0.046	0.010%
50	相关检测法	447	447.010	0.010	0.002%	0.045	0.010%
	频谱分析法	447	447.002	0.002	0.001%	0.034	0.008%

下载: 导出CSV 

| 显示表格

表  4  采样率为10 Hz, 采样点数为2000, 信号频率为0.5 Hz时的仿真结果

Table  4.  Simulation results of three signal detection methods on the condition of different SNR(sampling rate is 10 Hz, number of samples is 2000, signal frequency is 0.5 Hz)

信噪比/dB	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.507	0.007	1.371%	0.035	7.000%
-10	相关检测法	0.5	0.497	-0.003	0.502%	0.036	7.202%
	频谱分析法	0.5	0.498	-0.002	0.325%	0.028	5.629%
	数字滤波器法	1.5	1.496	-0.004	0.236%	0.036	2.434%
0	相关检测法	1.5	1.496	-0.004	0.240%	0.036	2.409%
	频谱分析法	1.5	1.498	-0.002	0.152%	0.028	1.864%
	数字滤波器法	4.5	4.496	-0.004	0.094%	0.036	0.809%
10	相关检测法	4.5	4.496	-0.004	0.087%	0.036	0.803%
	频谱分析法	4.5	4.497	-0.003	0.056%	0.028	0.619%
	数字滤波器法	15	14.995	-0.005	0.034%	0.036	0.243%
20	相关检测法	15	14.996	-0.004	0.026%	0.036	0.241%
	频谱分析法	15	14.997	-0.003	0.017%	0.028	0.186%
	数字滤波器法	45	44.993	-0.007	0.016%	0.036	0.081%
30	相关检测法	45	44.996	-0.004	0.008%	0.036	0.080%
	频谱分析法	45	44.997	-0.003	0.006%	0.028	0.062%
	数字滤波器法	141.5	141.486	-0.014	0.010%	0.036	0.026%
40	相关检测法	141.5	141.497	-0.003	0.002%	0.036	0.026%
	频谱分析法	141.5	141.497	-0.003	0.002%	0.028	0.020%
	数字滤波器法	447	446.963	-0.037	0.008%	0.036	0.008%
50	相关检测法	447	447.000	-0.000	0.000%	0.036	0.008%
	频谱分析法	447	446.997	-0.003	0.001%	0.028	0.006%

下载: 导出CSV 

| 显示表格

表  5  采样率为10 Hz, 采样点数为2000, 信号频率为1 Hz时的仿真结果

Table  5.  Simulation results of three signal detection methods on the condition of different SNR(sampling rate is 10 Hz, number of samples is 2000, signal frequency is 1 Hz)

信噪比/dB	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.514	0.014	2.743%	0.048	9.317%
-10	相关检测法	0.5	0.502	0.002	0.489%	0.053	10.529%
	频谱分析法	0.5	0.500	0.000	0.037%	0.033	6.576%
	数字滤波器法	1.5	1.505	0.005	0.334%	0.049	3.278%
0	相关检测法	1.5	1.501	0.001	0.059%	0.053	3.527%
	频谱分析法	1.5	1.500	0.000	0.031%	0.033	2.189%
	数字滤波器法	4.5	4.505	0.005	0.111%	0.049	1.095%
10	相关检测法	4.5	4.500	0.000	0.009%	0.053	1.176%
	频谱分析法	4.5	4.499	-0.001	0.015%	0.033	0.729%
	数字滤波器法	15	15.005	0.005	0.033%	0.049	0.329%
20	相关检测法	15	15.000	0.000	0.002%	0.053	0.353%
	频谱分析法	15	14.999	-0.001	0.005%	0.033	0.219%
	数字滤波器法	45	45.005	0.005	0.011%	0.049	0.110%
30	相关检测法	45	45.001	0.001	0.001%	0.053	0.118%
	频谱分析法	45	44.999	-0.001	0.002%	0.033	0.073%
	数字滤波器法	141.5	141.505	0.005	0.004%	0.049	0.035%
40	相关检测法	141.5	141.501	0.001	0.001%	0.053	0.037%
	频谱分析法	141.5	141.499	-0.001	0.001%	0.033	0.023%
	数字滤波器法	447	447.005	0.005	0.001%	0.049	0.011%
50	相关检测法	447	447.004	0.004	0.001%	0.053	0.012%
	频谱分析法	447	446.999	-0.001	0.000%	0.033	0.007%

下载: 导出CSV 

| 显示表格

图  5  3种方法在不同频率条件下的相对误差(左)和相对均方差(右)随 信噪比的变化结果比较(采样率为10 Hz, 采样点数为2000)

Figure  5.  Comparison of relative error (left) and relative standard deviation (right) with the variation of SNR on the condition of different frequency for three signal detection methods (sampling rate is 10 Hz, number of samples is 2000)

下载: 全尺寸图片 幻灯片

从表 3—5及图 5中可以看出，当信号频率一定时，相关检测法与频谱分析法的差值和均方差均不随信噪比的变化而变化，说明当噪声序列固定时，相关检测法和频谱分析法的检测结果仅受到噪声序列的影响，与正弦信号的振幅大小无关；当信号频率一定时，数字滤波器法的检测差值随信噪比的增大逐渐变小，这是因为数字滤波器法采用带通滤波器对输入信号进行滤波. 由于数字滤波器阶数的限制，数值模拟中滤波器的通带波动不能够达到理想状态，因此在不同频点处数字滤波器对输入信号会存在不同程度的抑制作用．根据本文仿真中选用的滤波器，其通带波动为0.01 dB，说明滤波器对不同频点的抑制在0—0.01 dB之间波动．假设输入信号频率为f₀，振幅为A₀，在频点f₀处噪声的影响幅度为B，滤波器的抑制为δ，那么最终数字滤波器法检测出的信号振幅为

那么，检测得到的误差结果为

因此，当信噪比变化，即被检测信号振幅变化时，由于所使用的带通滤波器没有变，通带范围内噪声的影响幅度不变，滤波器对已知频点的抑制作用不变，由式(9)可知，数字滤波器法的检测误差将随信噪比的增大而减小．

3.2   不同采样率的仿真结果比较

通过模拟仿真，给出了在一定采样时间和信噪比的条件下，当采样率不同时，数字滤波器法、 相关检测法和频谱分析法的检测效果．

假定采样时间为200 s，信噪比为-10 dB，采样率为10—30 Hz，表 6和表 7分别给出了信号频率为0.1 Hz和1 Hz时的相对误差和相对均方差的仿真结果． 可以看出，在采样率为10—30 Hz、 信号频率分别为0.1 Hz和1 Hz时，频谱分析法的相对误差结果和相对均方差结果均最小，信号检测效果要优于相关检测法和数字滤波器法．

表  6  采样时间为200 s，信噪比为-10 dB，信号频率为0.1 Hz时的仿真结果

Table  6.  Simulation results of three signal detection methods on the condition of different sampling rate (SNR is -10 dB, sampling time is 200 s, signal frequency is 0.1 Hz)

采样率/Hz	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.516	0.016	3.269%	0.040	7.702%
10	相关检测法	0.5	0.509	0.009	1.832%	0.041	8.052%
	频谱分析法	0.5	0.501	0.001	0.170%	0.025	5.004%
	数字滤波器法	0.5	0.509	0.009	1.716%	0.025	4.849%
20	相关检测法	0.5	0.503	0.003	0.544%	0.026	5.150%
	频谱分析法	0.5	0.499	-0.001	0.153%	0.025	4.994%
	数字滤波器法	0.5	0.507	0.007	1.374%	0.020	4.010%
30	相关检测法	0.5	0.503	0.003	0.696%	0.020	4.057%
	频谱分析法	0.5	0.502	0.002	0.392%	0.020	3.982%

下载: 导出CSV 

| 显示表格

表  7  采样时间为200 s，信噪比为-10 dB，信号频率为1 Hz时的仿真结果

Table  7.  Simulation results of three signal detection methods on the condition of different sampling rate (SNR is -10 dB, sampling time is 200 s, signal frequency is 1 Hz)

采样率/Hz	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.537	0.037	7.342%	0.044	8.138%
10	相关检测法	0.5	0.507	0.007	1.445%	0.045	8.837%
	频谱分析法	0.5	0.503	0.003	0.643%	0.034	6.789%
	数字滤波器法	0.5	0.517	0.017	3.440%	0.028	5.491%
20	相关检测法	0.5	0.500	0.000	0.040%	0.027	5.400%
	频谱分析法	0.5	0.500	0.000	0.029%	0.025	5.011%
	数字滤波器法	0.5	0.516	0.016	3.173%	0.022	4.329%
30	相关检测法	0.5	0.506	0.006	1.276%	0.025	4.891%
	频谱分析法	0.5	0.504	0.004	0.854%	0.025	5.016%

下载: 导出CSV 

| 显示表格

3.3   不同采样点数的仿真结果比较

通过模拟仿真，给出了在一定的采样率和信噪比条件下，当采样点数(采样时间)和信号频率不同时，数字滤波器法、 相关检测法和频谱分析法的检测效果．

假定采样率为10 Hz，信噪比为-10 dB，采样点数为2000—10000点，表 8和表 9分别给出了信号频率为0.1 Hz和1 Hz时的相对误差和相对均方差的仿真结果． 可以看出： 随着采样点数的增加，相关检测法和频谱分析法的差值变化不大，频谱分析法的相对误差结果要优于其它两种方法；随着采样点数的增加，信号检测的均方差结果逐渐减小，其原因是采样点数的增加，能够提高频谱的分辨能力．

表  8  采样率为10 Hz，信噪比为-10 dB, 信号频率为0.1 Hz时的仿真结果

Table  8.  Simulation results of three signal detection methods on the condition of different sampling number (SNR is -10 dB, sampling rate is 10 Hz, signal frequency is 0.1 Hz)

采样点数	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.537	0.037	7.342%	0.044	8.138%
2000	相关检测法	0.5	0.507	0.007	1.445%	0.045	8.837%
	频谱分析法	0.5	0.503	0.003	0.643%	0.034	6.789%
	数字滤波器法	0.5	0.532	0.032	6.371%	0.023	4.345%
4000	相关检测法	0.5	0.501	0.001	0.211%	0.023	4.495%
	频谱分析法	0.5	0.503	0.003	0.539%	0.021	4.128%
	数字滤波器法	0.5	0.534	0.034	6.857%	0.022	4.114%
5000	相关检测法	0.5	0.504	0.004	0.739%	0.021	4.202%
	频谱分析法	0.5	0.503	0.003	0.516%	0.019	3.688%
	数字滤波器法	0.5	0.531	0.031	6.227%	0.018	3.384%
8000	相关检测法	0.5	0.501	0.001	0.103%	0.017	3.383%
	频谱分析法	0.5	0.502	0.002	0.429%	0.017	3.332%
	数字滤波器法	0.5	0.534	0.034	6.746%	0.016	2.996%
10000	相关检测法	0.5	0.503	0.003	0.560%	0.015	2.940%
	频谱分析法	0.5	0.502	0.002	0.490%	0.014	2.839%

下载: 导出CSV 

| 显示表格

表  9  样率为10 Hz，信噪比为-10 dB、 信号频率为1 Hz时的仿真结果

Table  9.  Simulation results of three signal detection methods on the condition of different sampling number (SNR is -10 dB, sampling rate is 10 Hz, signal frequency is 1 Hz)

采样点数	方法	信号振幅/mV	检测均值/mV	差值/mV	相对误差	均方差	相对均方差
	数字滤波器法	0.5	0.514	0.014	2.743%	0.048	9.317%
2000	相关检测法	0.5	0.502	0.002	0.489%	0.053	10.529%
	频谱分析法	0.5	0.5	0.000	0.037%	0.033	6.576%
	数字滤波器法	0.5	0.508	0.008	1.695%	0.027	5.363%
4000	相关检测法	0.5	0.499	-0.001	0.269%	0.029	5.789%
	频谱分析法	0.5	0.5	0.000	0.060%	0.025	4.965%
	数字滤波器法	0.5	0.507	0.007	1.323%	0.021	4.166%
5000	相关检测法	0.5	0.499	-0.001	0.263%	0.020	4.097%
	频谱分析法	0.5	0.5	0.000	0.072%	0.020	3.903%
	数字滤波器法	0.5	0.507	0.007	1.309%	0.024	4.832%
8000	相关检测法	0.5	0.497	-0.003	0.541%	0.025	5.024%
	频谱分析法	0.5	0.498	-0.002	0.319%	0.020	4.104%
	数字滤波器法	0.5	0.509	0.009	1.750%	0.014	2.748%
10000	相关检测法	0.5	0.5	0.000	0.058%	0.014	2.791%
	频谱分析法	0.5	0.499	-0.001	0.123%	0.014	2.713%

下载: 导出CSV 

| 显示表格

从上述仿真结果中可以得出以下结论：

1)相同条件下，信噪比越高，信号检测的相对误差和相对均方差越小，检测效果越好．信噪比为-10—50 dB时，频率分析法和相关检测法的信号检测结果要优于数字滤波器法；当信噪比高于20 dB时，3种信号检测方法的检测效果差别不大，且随着信噪比的提高，相对误差和相对均方差的变化趋于平稳．

2)相同条件下，采样率为10—30 Hz时，频谱分析法的相对误差和相对均方差最小，信号检测效果优于相关检测法和数字滤波器法．

3)相同条件下，随着采样点数的增加，信号检测的均方差逐渐减小，相关检测法和频谱分析法的差值变化不大．从相对误差结果上看，频谱分析法和相关检测法的检测效果要优于数字滤波器法．

4)信号频率一定时，对于同一噪声序列，相关检测法和频谱分析法的检测结果仅与噪声序列的频谱相关，与正弦信号的振幅无关．

4.   讨论与结论

本文利用数值方法模拟频谱分析法、 相关检测法和数字滤波器法等3种信号检测方法，得出在不同参数条件下对单一频率输入信号的检测精度. 通过对比分析检测结果，并结合3种方法在实际系统中实现的难易程度，给出适用于地电阻率交流观测系统的最优信号检测方法.

从仿真结果上看，频谱分析法和相关检测法的相对误差和相对均方差结果均优于数字滤波器法．其主要原因在于，数字滤波器法中所使用的带通滤波器有一定的通带范围，所受到的噪声影响较大；另外，受滤波器阶数和仿真参数的限制，带通滤波器的纹波因数不能达到理想状态，对被测信号会产生一定程度的抑制. 两种影响的同时作用使得数字滤波器法的检测效果不够理想．

相关检测法的检测效果与频谱分析法差别不大．但是，相关检测法在使用时需要满足两个条件：① 该检测方法中的参考信号需要与被测信号的频率完全一致，这就要求在实际应用时低频电源发射的供电信号需要非常稳定；② 该检测方法中需要进行低通滤波，低通滤波器的通带范围应该在直流附近，通带越窄越好． 若采用有限冲激响应滤波器设计低通滤波器的话，会需要很高阶数的滤波器，从而导致相关检测法需要占用更多的系统资源．并且由于滤波器性能的限制，相关检测法的检测效果也会受到影响．因此，在实际应用中，该检测方法的实现难度较大，占用的系统资源也较多．

与相关检测法相比，频谱分析法的处理方法更简单，易于实现．但是，该检测方法的前提是噪声频率与被提取信号的频率不在同一个频段．因此，在实际观测时，需要提前对场地进行测试，找到干扰比较小的频段，选择合适的信号频率．

目前直流地电阻率观测的相对误差要求不大于0.3%(中国地震局，2009)，因此采用地电阻率交流测量方法时，在采样率不变的前提下，通过提高采样时间长度，以增加频谱分析时的数据长度，提高频率分辨力；或者采用滑动多次平均的方法来提高频谱分析法对数据处理的准确度，以提高检测精度．

综上所述，频谱分析法是适用于地电阻率交流观测的一种较好的信号检测技术，希望能够为地电阻率交流观测系统的研制提供借鉴．