基于阶梯波参考信号的改进型正交锁定放大
周逢道1,2, 丁凯来1, 曾新森1, 薛开昶1, 孙彩堂1,2
1.吉林大学 仪器科学与电气工程学院,长春 130061
2.吉林大学 国家地球物理探测仪器工程技术研究中心,长春 130061
通讯作者:孙彩堂(1975-),男,副教授,博士.研究方向:数据处理与可视化.E-mail:sunct@jlu.edu.cn

作者简介:周逢道(1970-),男,教授,博士.研究方向:电磁探测技术及仪器.E-mail:zhoufd@jlu.edu.cn

摘要

对以阶梯波作为参考信号的正交锁定放大的探测精度进行了建模分析,基于硬件简单原则给出了以阶梯波作为参考信号的正交锁定放大信号检测方法的实现策略。综合考虑参考波形的失真度与采样率之间的关系,对比分析了基于面积等效法和特定谐波消去法的阶梯波合成方法。根据阶梯波参考信号对采集电路相位精度的要求及待检测信号的周期性,提出了匹配采样技术。实测结果表明:采用现场可编程门阵列(Field programmable gate array, FPGA)实现所提出的正交锁定放大,方法简单、硬件要求低、可用于识别不同异常体。改进后方法的误差为0.7%~7%,与以方波为参考信号的正交锁定放大相比,提高了探测精度。

关键词: 信息处理技术; 近地表; 正交锁定放大; 阶梯波
中图分类号:TH763.1 文献标志码:A 文章编号:1671-5497(2016)03-0996-08
Improved orthogonal lock-in amplification based on step wave reference signal
ZHOU Feng-dao1,2, DING Kai-lai1, ZENG Xin-sen1, XUE Kai-chang1, SUN Cai-tang1,2
1.College of Instrumentation and Electrical Engineering, Jilin University, Changchun 130061, China
2.National Geophysical Exploration Equipment Engineering Research Center, Jilin University,Changchun 130061, China
Abstract

The detection accuracy of the step wave reference orthogonal lock-in amplification is modeled and analyzed, and the implementation strategy based the principle of using simple hardware is presented. Considering the relationship between the frequency waveform distortion and the sampling rate, the synthesis strategies of step wave based on equivalent area method and selected harmonic elimination method are comparatively analyzed. According to the phase accuracy requirement of step wave reference signal acquisition and the periodicity of the signal to be tested, a matching sampling technique is proposed. Experiments results show that the proposed orthogonal lock-in amplification can achieved by Field Programmable Gate Array (FPGA), and the method is simple and has low hardware requirement, can be used to identify different anomaly bodies. The error of the improved method is 0.7% to 7.0%. Compared with the square wave reference orthogonal lock-in amplification, the detection accuracy is improved by the new method.

Keyword: information processing technology; near-surface; orthogonal lock-in amplification; step wave
0 引 言

近地表电磁探测是指在电磁波一次场的激励下, 通过采集由地下异常所产生的二次场, 获取地下介质电磁特性的方法, 被广泛应用于地下管线检测、工程质量监测和考古勘测等领域[1]。近地表电磁探测中, 低频部分数据反映深层地下结构信息, 而高频部分数据反映浅层地下结构信息, 一般频率范围为1~40 kHz, 频率近似按照对数等间隔分布, 每10倍频程中分布着3~10个频率点。同时, 为了提高野外工作效率, 仪器应尽可能轻便、低功耗并具有一定的实时处理能力[2]。因此, 对能满足上述要求的数据检测方法进行研究具有较高的实际应用价值。

在近地表电磁探测中, 方波激励源与正弦波激励源相比具有如下优点[3, 4, 5, 6, 7]:①正弦波发射一般采用正弦脉冲宽度调制(Sinusoidal pulse width modulation, SPWM)技术, 调制过程引入交调失真, 导致发射波形随机度大、探测精度低, 而方波发射时不需要调制; ②产生方波时发射电路只需要工作在发射频率, 而产生正弦波时一般需要工作在发射频率的数倍以上, 发射方波具有低的功率管损耗; ③在相同的母线电压等级下, 方波的基波幅度是合成正弦波的4/π 倍, 方波能量更高。

由于采用方波作为一次场激励源, 给常规正交锁定放大技术带来了不便[8, 9, 10, 11]。常规正交锁定放大中的参考信号一般采用正弦波或方波。若选用正弦波作为参考信号, 不可避免地要用乘法运算实现正交锁定放大中的相敏检测(Phase sensitive detection, PSD), 通常需要基于电脑进行运算, 具有高的硬件功耗及成本; 若选用方波作为参考信号, PSD可以通过电子开关实现[12], 故功耗低、成本低。但在方波激励源的作用下, 因接收信号含有与参考信号频率一致的谐波成份, 故具有较大的检测误差, 虽然可以采用选频带通滤波器抑制接收信号中的谐波成份, 但带通滤波器的中心频率要时刻匹配检测频率, 难以实现。

本文提出了一种采用阶梯波逼近正弦参考信号的方法实现了正交锁定放大。因阶梯参考信号与方波参考信号相比, 具有低的谐波成份, 故精度高于方波参考的正交锁定放大。可以采用在现场可编程门阵列(Field programmable gate array, FPGA)内部对数据进行移位的方式实现阶梯波参考PSD运算, 故不需要做乘法运算, 具有较低的功耗及硬件成本。

1 正交锁定放大模型
1.1 正交锁定放大原理

图1为近地表电磁探测的正交锁定放大原理。发射电路向发射线圈发送方波激励信号vT:

vT=k=1, 3, 5, 4Asin(kωt)(1)

式中:A表示双极性方波的幅值。

发射线圈产生一次激励场H1并作用于地下异常体, 感应产生二次场H2, 因此, H2包含了地下信息。由发射线圈直接耦合至接收线圈的一次场H1和感应二次场H2在接收线圈上产生感应电压vR, 因为发射线圈、地下介质和接收线圈可以认为是线性时不变系统, 故vR可以用式(2)表示:

vR=k=1, 3, 5, aRksinkωt+θRk(2)

式中:aRk为第k次频率成份对应的幅度; θ Rk为相对于发射波形, 第k次频率成份的相移。

vR经正交锁定放大得到同相分量I和正交分量Q。正交参考信号vIvQ通过同步控制与vT保持同步, 对于周期性的参考信号, vIvQ分别如式(3)(4)所示:

vI=i=1, 3, 5, aIQisiniωt(3)

vQ=i=1, 3, 5, aIQisiniωt-π2i(4)

式中:aIQi表示vIvQ中第i次频率成份对应的幅度。

vIvT同相, vQvI滞后π /2。正交锁定放大中PSD等价于乘法运算; 低通滤波器(Low pass filter, LPF), 用于获取直流分量。IQ分别为vRvIvQ乘积的直流分量。IQ分别如式(5)和(6)所示:

I=aIQ1aR12sinθR1+i=k=1, 3, 5, aIQiaRk2sinθRk(5)

Q=aIQ1aR12cosθR1+i=k=1, 3, 5, aIQiaRk2sinθRk+π2i(6)

当按式(7)(8)取基波幅度aR1和相移θ R1时, 式(5)(6)中等号右侧第1项即为精确成份, 第2项为误差成份。

aR1=2aIQ1·I2+Q2(7)

θR1=arctanIQ(8)

因此, 减小式(3)(4)中i≠ 1的频率成分能减小式(5)(6)中的误差项, 进而提高正交锁定放大的精度。减小式(3)(4)中i≠ 1的频率成分可以转换为降低参考信号的失真度。

图1 电磁探测的正交锁定放大原理Fig.1 Orthogonal lock-in amplification principle of electromagnetic prospecting

1.2 阶梯波参考的PSD

为了降低参考信号vIvQ的失真度, 可采用图2所示的阶梯波x(ω t)逼近正弦波。

图2 阶梯波的正弦逼近原理Fig.2 Sinusoidal approximation principle of step wave

图3 阶梯波的PSD原理Fig.3 PSD principle of step wave

采用阶梯波参考作PSD的原理图如图3所示。图3中vR(m)为经模数转换器(Analog-to-digital converter, ADC)采样后的离散化接收信号, 数据以二进制的补码形式表示; x(ω t)为阶梯波参考信号, 设阶梯级数为n, 阶梯幅值步进为1。当一个接收信号的采样点vR(m)到来时, 首先判断其对应阶梯波的哪一段, 并取出相应段的幅值绝对值|x(ω t)|, 以二进制数Dh, …, D1D0表示。当Dh, …, D1D0中的第l(l=0, 1, …, n)位为1时, 将vR(m)左移l位实现vR(m)与2l的乘积, 并将移位结果保存。按照相同的方法处理所有Dh, …, D1D0中为1的位, 将移位后保存的结果全部累加, 当x(ω t)为负时, 可通过将移位后的结果逐位取反, 然后加1的方式实现取相反数的运算; 当x(ω t)为0时, 将vR(m)置零。通过这种方法即可完成一个采样点的阶梯波PSD计算。

由于阶梯波参考信号通过拟合正弦波的方式实现, 故失真度低于方波, 精度高于方波。由于阶梯波参考信号与接收信号进行PSD运算时, 可以通过移位累加的方式实现, 故硬件要求及功耗低于正弦波参考信号。

2 阶梯波参考信号的合成

阶梯波参考信号的失真度与电平切换的相角有直接关系, 通常采用特定谐波消去法[13, 14]和等效面积法[15, 16]确定各相角。由于这两种方法主要应用于电力电子领域, 与本文的应用背景及关注点有较大的不同, 故这一节将针对正交锁定放大这一特殊目的进行分析。

2.1 阶梯波失真度分析模型

图4为阶梯波的合成原理, 设所需合成的阶梯波x(ω t)幅值为n, 共有n级阶梯, 可由图4(b)(c)(d)所示的幅值为1的矩形波叠加得到。为了分析各级阶梯波x1(ω t), x2(ω t), …, xn(ω t)的频率成分, 将它们进行傅里叶级数展开, 展开式如(9)所示:

xjωt=i=1, 3, 5, aijsin(iωt)aij=4cosiθj, j=1, 2, , n(9)

根据叠加定理, 所合成的阶梯波x(ω t)的频率成分可由各级阶梯波的展开式相对应的频率成分叠加得到, 如式(10)所示:

xωt=j=1nxjωt=i=1, 3, 5, Aisin(iωt)Ai=4j=1ncosiθj(10)

图4 阶梯波的合成原理Fig.4 Synthesis principle of step wave

由式(10)可知各次频率成分的幅度Ai。结合式(10)及帕斯瓦尔定理, 总谐波失真度(Total harmonic distortion, THD)可用式(11)表示:

THD=i=3, 5, 7, Ai/22A1/2=RMS2-A1/22A1/2(11)

式中:RMS为x(ω t)的有效值:

RMS=12π-ππx2(ωt)dωt

由式(9)~(11)可知, x1(ω t), x2(ω t), …, xn(ω t)的相位角θ 1, θ 2, …, θ n直接影响THD。

2.2 特定谐波消去法合成的阶梯波

由式(10)可知, 当i=3, 5, 7, …时, 各次谐波成分的幅度Aii成反比。因此, 对于具有n个阶梯的阶梯波, 可以通过控制θ 1, θ 2, …, θ n使如式(10)所示的谐波分量A3, A5, …, A2n-1为零。根据这一原理, 可以采用式(12)所示方程求取θ 1, θ 2, …, θ n:

A1=4πj=1ncosθj=nA3=43πj=1ncos3θj=0A5=45πj=1ncos5θj=0    A2n-1=42n-1πj=1ncos2n-1θj=00< θ1< θ2< < θn< π/2(12)

由式(12)可知, 虽然方程中等式个数与未知数的个数相等, 但θ 1, θ 2, …, θ n是与三角函数相关的隐函数, 难以精确求解, 故采用计算机进行迭代计算。表1表2分别为n=1~4时的相位角、频率成份值及THD。表2中, f0表示基波的频率, 幅值做了归一化处理, 即除以阶数n表4~表7中幅度值也做了归一化。

表1 特定谐波消去法的相位角(° ) Table 1 Phase angles of selected harmonic elimination(° )
表2 特定谐波消去法的频率成分及THD Table 2 Frequency components and THD of selected harmonic elimination

表2可知, 特定谐波消去法能将3, 5, …, 2n-1次谐波精确消去。

2.3 等效面积法合成的阶梯波

在电力电子领域的SPWM中, 通过使矩形单元的面积与正弦波单元的面积相等来使矩形阶梯波的低频成份与正弦波接近。借鉴这一理论, 可以通过图5所示方法来合成阶梯波。图5(a)为正弦波单元S1, S2, …, Sn的划分, S1S2, …, Snθ S1, θ S2, …, θ Sn作为各正弦波单元的划分, θ Sj表示nsinθ Sj=j(j=1, 2, …, n)时对应的相角; 图5(b)为矩形单元Sq1Sq2, …, Sqn的划分。令Sj=Sqj相当于:

0θ1nsinωtdωt=θ1θS11-nsinωtdωt(13a)

θS(j-1)θjnsinωt-(j-1)dωt=θjθSjj-nsinωtdωt(13b)

θS(n-1)θnnsinωt-(n-1)dωt=θnπ2n-nsinωtdωt(13c)

式中:j=2, 3, …, n-1; θ S1, θ S2, …, θ S(n-1)=sin-1(1/n), sin-1(2/n), …, sin-1[(n-1)/n]。

图5 等效面积法原理Fig.5 Principle of equivalent area

根据式(13a)(13b)(13c)获取阶梯波的相位角θ 1, θ 2, …, θ n

表3表4分别为n=1~5时的相位角、频率成份值及THD。

表3 等效面积法的相位角 Table 3 Phase angles of equivalent area
表4 等效面积法的频率成分及THD Table 4 Frequency components and THD of equivalent area
2.4 阶梯波合成方法的比较

表5为单位幅值方波的主要频率成分及THD, 将表5表2表4对比可知, 阶梯波的谐波及THD显著低于方波。

表5 方波的频率成分及THD Table 5 Frequency components and THD of square wave

对比表2表4可知, 在相同阶梯数n下, 虽然特定谐波消去能精确消去奇次谐波, 但其THD仍低于等效面积法。同时, 特定谐波消去的相角有分布极近的情况, 如n=3时θ 1=19.86° θ 2=20.24° , 其对采样频率有过高的要求。

在实际应用中, 受采样频率的限制, 实际可实现的相位精度一般比表1表3中数据低得多。因此, 表6给出了n=4的阶梯波, 单个相角依次偏差+3° 时两种方法的THD。

表6 相位偏差对特定谐波消去法和等效面积法THD的影响 Table 6 Effect of phase deviation in selected harmonic elimination and equivalent area

表6中, THD1、THD2分别为特定谐波消去法和等效面积法在相位偏差时的THD。由表6表2表4对比可知:n=4, 某相位角偏差为+3° 时, THD1误差为-0.0091~+0.0110, THD2误差为+0.0048~+0.0053, 等效面积法受相位精度的影响较低。

综合表2表4表6的比较结果, 等效面积法合成的阶梯波更适合所提出的正交锁定放大。

综合计算结果及工程应用的可行性, 当n=3~5时, 阶梯波的相位角可近似取表7所示的值, 其THD与表4相比差别不大。需特别注意:获取IQ应除以阶梯波的基波幅值, 阶梯波的基波幅值(即表7中1f0)不是1。

表7 工程近似阶梯波参数 Table 7 Parameters of engineering approximation step wave
3 匹配采样技术

在数据采集过程中, 采样频率越高, 数据对二次场信号描述越细致, 与阶梯波相角匹配得越准确, 但是采用过高的采样频率会对硬件电路提出更高要求。当采用图6所示的降采样频率原理后, 可以用相对低的采样频率实现较高采样频率的效果。对于图6(a)所示的时域信号, 它的周期为T, 当要保证一个周期内采集到8个点时, 对于常规技术, 采样间隔Ts=(1/8)T, 即采样频率为信号频率的8倍。然而, 当取Ts=(3/8)T时, 按图6(b)所示方法, 则可以在3个信号周期T内重构一个周期的8个点。显然, 图6(b)用常规技术1/3倍的采样频率重构了相同相位精度的周期信号。

该原理可以定义为:设希望将一个周期为T的信号进行N点采样, 当取K为一个不能被N整除的质数, 采样周期Ts=(K/N)T时, 在K个信号周期T内可重构N个相位不同的采样点, 采样频率可以比常规技术低K倍。

图6 降采样频率原理Fig.6 Principle of sampling rate reduction

4 试验结果

在图7所示的硬件平台上进行试验。采用ALTERA公司的EP3C40Q240C8系列FPGA, 用于实现对发射电路和ADC采集的同步控制以及n=4的阶梯波正交锁定放大, 在采样间隔即可完成运算, 节约了探测时间。上位机用于进行传统的正弦波和方波参考的正交锁定放大。USB模块CY7C68013用于将FPGA运算结果传输到上位机, 便于在上位机中进行本文方法与传统方法的结果对比。

图7 硬件平台的结构框图Fig.7 Structure diagram of hardware platform

4.1 采样结果

将接收线圈上的感应信号vR适当放大后, 采用泰克TDS1012C示波器获取的波形图如图8所示。在相同采样率下, 采用本文所述采样技术和常规采样技术对vR进行采样, 获得的波形分别如图9中实线和虚线所示。将图9中波形与图8中示波器波形比较可知, 本文所提出的降采样频率技术能更好地体现波形细节, 具有更高的相位精度。

图8 示波器波形Fig.8 Waveform in oscilloscope

图9 本文方法和常规方法的采样波形Fig.9 Sampling waveforms using the proposed method and conventional method

4.2 金属异常探测结果

对两个大小相同的异常体进行探测实验, 一个为圆形铜片, 属于非铁磁性材料, 另一个为圆形铁片, 属于磁性材料。采用正弦波、方波及本文所述的阶梯波作为参考信号进行正交锁定放大, 对于铜片和铁片, 获取IQ曲线分别如图10(a)(b)所示。

图10 铜片和铁片的正交锁定放大结果Fig.10 Orthogonal lock-in amplification results of copper and iron plate

将方波、阶梯波、正弦波的测试结果作对比可知, 采用方波参考的探测结果在正弦波参考的探测结果两侧大幅跳动, 而以阶梯波参考的结果较为接近于以正弦波参考的探测结果。对比图10(a)和图10(b)可知, 铜片和铁片具有不同的IQ分量特性, 阶梯波的IQ与正弦波的IQ走势一致, 故本文所提出的正交锁定放大技术可用于近地表电磁探测。

表8给出了图10(a)(b)中阶梯波探测结果误差最大的一条曲线中的数值, 即图10(a)中铜片的I分量值, 用于定量分析方波和阶梯波参考下正交锁定放大的误差。表中的方波误差和阶梯波误差均以正弦波参考下的数值为标准求取。

表8 铜片的I分量 Table 8 I of copper plate

表8可知:在本文所述阶梯波参考下, I分量的误差为0.7%~7.0%, 远低于方波参考下的误差-14.3%~+18.2%。

5 结束语

提出了基于阶梯波参考的正交检测方法。论述了阶梯波合成原理, 推导了阶梯波谐波表达式, 用特定谐波消去法和等面积法分别确定了阶梯波各级阶梯的相角, 并对比了两种方法所合成的阶梯波的频率成分及THD。采用匹配采样方法更细致地还原了信号波形, 以便与阶梯波进行更准确的运算。试验结果表明:所提出的正交锁定放大方法中的阶梯波相角更适合用等面积法求出, 阶梯波与接收信号的PSD运算可通过在FPGA内部用移位累加来实现, 无需通过乘法器进行运算, 实现方法简单、具有低的硬件要求、精度高于方波正交锁定放大, 可用于识别不同异常体。

The authors have declared that no competing interests exist.

参考文献
[1] 周逢道, 王金玉, 唐红忠, . 近地表电磁探测多频数字驱动信号产生技术[J]. 吉林大学学报: 工学版, 2013, 43(3): 682-687.
Zhou Feng-dao, Wang Jin-yu, Tang Hong-zhong, et al. Multi-frequency digital drive signal generation technology in near surface electromagnetic detection domain[J]. Journal of Jilin University (Engineering and Technology Edition), 2013, 43(3): 682-687. [本文引用:1]
[2] 张文秀, 林君, 刘立超, . 分布式电磁探测宽频数据采集系统设计与实现[J]. 吉林大学学报: 工学版, 2012, 42(6): 1426-1431.
Zhang Wen-xiu, Lin Jun, Liu Li-chao, et al. Design and implementation of broadband data acquisition system for distributed electromagnetic exploration[J]. Journal of Jilin University(Engineering and Technology Edition), 2012, 42(6): 1426-1431. [本文引用:1]
[3] 薛开昶, 董威, 林君, . 逆m伪随机相位PWM逆变器的EMI峰值抑制方法[J]. 吉林大学学报: 工学版, 2013, 43(6): 1589-1595.
Xue Kai-chang, Dong Wei, Lin Jun, et al. Suppression method of inverse m sequence rand om phase PWM inverter’s EMI peak[J]. Journal of Jilin University (Engineering and Technology Edition), 2013, 43(6): 1589-1595. [本文引用:1]
[4] Zitoun A, Beheim L, Huez R, et al. Smart electrom-agnetic sensor for buried conductive targets identi-cation[J]. IEEE Sensors Journal, 2006, 6(6): 1580-1591. [本文引用:1]
[5] Svatos J, Vedral J, Fexa P. Metal detector excited by frequency-swept signal[J]. Metrology and Measurement Systems, 2011, 18(1): 57-67. [本文引用:1]
[6] Ramachand ran G, Gader P D, Wilson J N. GRANMA: Gradient angle model algorithm on wideband EMI data for land -mine detection[J]. Geoscience and Remote Sensing Letters, 2010, 7(3): 535-539. [本文引用:1]
[7] Shubitidze F, Fernand ez J P, Barrowes B E, et al. The orthonormalized volume magnetic source model for discrimination of unexploded ordnance[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(8): 4658-4670. [本文引用:1]
[8] 罗海军, 何为, 徐征, . 基于同步检波的单通道磁感应成像技术研究[J]. 仪器仪表学报, 2012, 44(4): 899-904.
Luo Hai-jun, He Wei, Xu Zheng, et al. Study on single-channel magnetic induction tomography technique based on synchronous detection[J]. Chinese Journal of Scientific Instrument, 2012, 44(4): 899-904. [本文引用:1]
[9] Aguirre J, Garcia R D, Medrano N, et al. Square-signal-based algorithm for analog lock-in amplifiers[J]. IEEE Transactions on Industrial Electronics, 2014, 61(10): 5590-5598. [本文引用:1]
[10] Qin J H, Huang Z M, Ge Y J, et al. Tand em demodulation lock-in amplifier based on digital signal processor for dualmodulated spectroscopy[J]. Review of Scientific Instruments, 2009, 80(3): 033112. [本文引用:1]
[11] Yang Y, Huang L, Lu Z L, et al. AC voltage reference using the fundamental of DAC stepwise-approximated sine waves[J]. IEEE Transactions on Instrumentation and Measurement, 2014, 63(2): 473-478. [本文引用:1]
[12] 高晋占. 微弱信号检测[M]. 2版, 北京: 清华大学出版社, 2011: 171-214. [本文引用:1]
[13] Suresh Y, Pand a A K. Research on a cascaded multi-level inverter by employing three-phase transformers[J]. IET Power Electronics, 2012, 5(5): 561-570. [本文引用:1]
[14] Young C M, Wu S F. Selective harmonic elimination in multi-level inverter with zig-zag connection transformers[J]. IET Power Electronics, 2014, 7(4): 876-885. [本文引用:1]
[15] Zhang Y, Zou Y P, Wang C Z, et al. A novel modulation technology for multilevel inverter based on equivalent area[C]//The 33rd Annual Conference of the IEEE Industrial Electronics Society, Taipei, Taiwan, 2007: 2115-2118. [本文引用:1]
[16] Hua C C, Wu C W, Chuang C W. Control of low-distortion 27-level cascade inverter with three H-bridge inverter modules[C]//IEEE International Conference on Industrial Technology, Bombay, India, 2006: 1844-1849. [本文引用:1]