Signal detection method for AC geo-resistivity observation
-
摘要: 目前我国地电阻率观测台站受地铁、 轻轨等干扰的情况越来越多， 严重影响地电阻率观测数据质量. 地电阻率交流观测方法被认为是一种可以解决此类干扰问题的有效方法. 本文首先介绍该观测方法的理论基础； 然后采用数值模拟方法对实际观测过程进行模拟， 给出在不同参数条件下， 分别采用数字滤波器法、 频谱分析法和相关检测法进行信号检测所得到的相对误差和相对均方差结果； 最后， 通过对3种信号检测方法进行比较和分析， 给出频谱分析法在地电阻率交流观测中可以作为最佳信号检测方法的建议， 为地电阻率交流观测系统的研制提供理论参考.
-
关键词:
- 地电阻率观测 /
- 地电阻率交流观测方法 /
- 选频接收 /
- 数值模拟
Abstract: In the past ten or twenty years, more and more geo-resistivity observation stations are being faced with increasingly electromagnetic interference, such as subway and rail traffic system in urban area, which seriously affects the quality of observation data. The AC geo-resistivity observation method has been found to be effective to solve this problem and can be used for earthquake monitoring and prediction. This paper introduced firstly the theory of the AC observation method, and then three different signal detection methods, i.e., the digital filter method, the spectrum analysis method and the correlation detection method, are used in the simulation. The analysis results are obtained by simulating an actual work process of AC geo-resistivity observation method. Comparison of the relative error and relative standard deviation on the condition of different signal-noise ratio, sampling rate or sampling time shows that frequency spectra method can be taken as the best method for the signal observation, providing theoretical basis for developing the AC geo-resistivity observation system. -
引言
目前，我国地电前兆台网广泛使用的地电阻率观测方法是直流观测方法，其观测装置采用对称四电极装置(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 methodA and B are current-emitting electrodes, and M and N are measuring electrodes地电阻率交流测量方法采用低频电流作为激励源，通过供电电极A和B向大地供入低频电流IAB，同时通过测量电极M和N测量MN之间的同频率响应信号UMN，通过计算得出地电阻率ρ的具体表达式为(张国民等，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)与参考信号r1(t)和r2(t)相乘，分别为y1(t)和y2(t)，即
式中，x(t)为被检测信号，n(t)为噪声信号，r1(t)和r2(t)为参考信号，A为被测信号振幅，ω0为被测信号频率，φ为相对相位. 上式中由于信号与噪声是不相关的，因此sin(ω0t+φ)n(t)和cos(ω0t+φ)n(t)两项的结果为零，令
这样，Y1(t)与Y2(t)经过低通滤波器后，频率为2ω0的高频项被滤波处理，可以得到 z1=AB/2cosφ和z2=AB/2sinφ两个结果项，被检测信号的振幅值则为
![]()
图 2 相关检测法的测量原理x(t)为含有噪声的信号，r1(t)和r2(t)为与x(t) 频率相同而相位不同的参考信号Figure 2. Measuring principle of correlation detection methodx(t) is the signals including the noise, r1(t) and r2(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πf0t+φ0)，式中f0为信号频率，φ0为初始相位，A为被检测信号振幅，正弦信号平均功率值Ps=A2/2．根据信噪比计算公式SNR=10lg(Ps/Pn)保持噪声信号不变，由不同信噪比得到不同正弦信号振幅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 输入信号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:
error on line 1 at column 1: Start tag expected, '<' not foundBelow 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 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% 表 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% 表 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% ![]()
图 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之间波动．假设输入信号频率为f0，振幅为A0，在频点f0处噪声的影响幅度为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% 表 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% 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% 表 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% 从上述仿真结果中可以得出以下结论：
1)相同条件下，信噪比越高，信号检测的相对误差和相对均方差越小，检测效果越好．信噪比为-10—50 dB时，频率分析法和相关检测法的信号检测结果要优于数字滤波器法；当信噪比高于20 dB时，3种信号检测方法的检测效果差别不大，且随着信噪比的提高，相对误差和相对均方差的变化趋于平稳．
2)相同条件下，采样率为10—30 Hz时，频谱分析法的相对误差和相对均方差最小，信号检测效果优于相关检测法和数字滤波器法．
3)相同条件下，随着采样点数的增加，信号检测的均方差逐渐减小，相关检测法和频谱分析法的差值变化不大．从相对误差结果上看，频谱分析法和相关检测法的检测效果要优于数字滤波器法．
4)信号频率一定时，对于同一噪声序列，相关检测法和频谱分析法的检测结果仅与噪声序列的频谱相关，与正弦信号的振幅无关．
4. 讨论与结论
本文利用数值方法模拟频谱分析法、 相关检测法和数字滤波器法等3种信号检测方法，得出在不同参数条件下对单一频率输入信号的检测精度. 通过对比分析检测结果，并结合3种方法在实际系统中实现的难易程度，给出适用于地电阻率交流观测系统的最优信号检测方法.
从仿真结果上看，频谱分析法和相关检测法的相对误差和相对均方差结果均优于数字滤波器法．其主要原因在于，数字滤波器法中所使用的带通滤波器有一定的通带范围，所受到的噪声影响较大；另外，受滤波器阶数和仿真参数的限制，带通滤波器的纹波因数不能达到理想状态，对被测信号会产生一定程度的抑制. 两种影响的同时作用使得数字滤波器法的检测效果不够理想．
相关检测法的检测效果与频谱分析法差别不大．但是，相关检测法在使用时需要满足两个条件：① 该检测方法中的参考信号需要与被测信号的频率完全一致，这就要求在实际应用时低频电源发射的供电信号需要非常稳定；② 该检测方法中需要进行低通滤波，低通滤波器的通带范围应该在直流附近，通带越窄越好． 若采用有限冲激响应滤波器设计低通滤波器的话，会需要很高阶数的滤波器，从而导致相关检测法需要占用更多的系统资源．并且由于滤波器性能的限制，相关检测法的检测效果也会受到影响．因此，在实际应用中，该检测方法的实现难度较大，占用的系统资源也较多．
与相关检测法相比，频谱分析法的处理方法更简单，易于实现．但是，该检测方法的前提是噪声频率与被提取信号的频率不在同一个频段．因此，在实际观测时，需要提前对场地进行测试，找到干扰比较小的频段，选择合适的信号频率．
目前直流地电阻率观测的相对误差要求不大于0.3%(中国地震局，2009)，因此采用地电阻率交流测量方法时，在采样率不变的前提下，通过提高采样时间长度，以增加频谱分析时的数据长度，提高频率分辨力；或者采用滑动多次平均的方法来提高频谱分析法对数据处理的准确度，以提高检测精度．
综上所述，频谱分析法是适用于地电阻率交流观测的一种较好的信号检测技术，希望能够为地电阻率交流观测系统的研制提供借鉴．
-
![]()
图 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
![]()
图 2 相关检测法的测量原理
x(t)为含有噪声的信号，r1(t)和r2(t)为与x(t) 频率相同而相位不同的参考信号
Figure 2. Measuring principle of correlation detection method
x(t) is the signals including the noise, r1(t) and r2(t) are the reference signals whose frequency are same as that of x(t) and whose phase are different from that of x(t)
![]()
图 3 10万点高斯白噪声序列(均值为0.0071，方差为1.0051)
Figure 3. 100000-point Gaussian white noise with mean 0.0071 and variance 1.0051
![]()
图 4 信噪比为0 dB、信号频率为1 Hz的输入信号x(t)的波形
Figure 4. Signal with Gaussian white noise(SNR is 0 dB, signal frequency is 1 Hz)
![]()
图 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)
表 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
表 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 采样率为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
表 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
表 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
-
门爱东, 苏菲, 王雷, 王海婴, 李江军. 2009. 数字信号处理[M]. 第二版. 北京: 科学出版社: 75-136. Men A D, Su F, Wang L, Wang H Y, Li J J. 2009. Digital Signal Processing[M]. 2nd ed. Beijing: Science Press: 75-136 (in Chinese).
万永革. 2012. 数字信号处理的MATLAB实现[M]. 第二版. 北京: 科学出版社: 86-132. Wan Y G. 2012. Implementation of Digital Signal Processing Based on MATLAB[M]. 2nd ed. Beijing: Science Press: 86-132 (in Chinese).
张国民, 傅征祥, 桂燮泰. 2001. 地震预报引论[M]. 北京: 科学出版社: 241-246. Zhang G M, Fu Z X, Gui X T. 2001. Introduction to Earthquake Prediction[M]. Beijing: Science Press: 241-246 (in Chinese).
中国地震局. 2009. 地震标准汇编[M]. 北京: 地震出版社: 1214-1220. China Earthquake Administration. 2009. Seismic Standards Compilation[M]. Beijing: Seismological Press: 1214-1220 (in Chinese).
计量
- 文章访问数: 547
- HTML全文浏览量: 264
- PDF下载量: 16
