本发明专利技术公开了一种广域电磁法工频滤波与实时计算技术，包括：基于现场可编程门阵列FPGA芯片实现，由FIR滤波器、巴特沃斯带阻滤波器、插值抽取模块、FFT计算模块组成，并基于凯瑟窗函数设计，采用凯塞窗函数设计FIR滤波器用于滤除高频干扰信号；采用巴特沃斯带阻滤波器用于滤除50Hz工频干扰噪声；采用插值与抽取方法将广域电磁法电场信号同一变换为2048点后，传入FFT IP核进行FFT计算。其特征在于，通过分析广域电磁法电场信号频率特征，设置滤波器参数和插值与抽取因子，同时对滤波器参数进行量化，利用FPGA可编程与并行处理的优势，对接收到的广域电磁法电场信号进行滤波与实时计算。解决了广域电磁法电场信号采集数据中工频噪声影响大和不能实时FFT计算的问题。频噪声影响大和不能实时FFT计算的问题。频噪声影响大和不能实时FFT计算的问题。
【技术实现步骤摘要】
一种广域电磁法工频滤波与实时计算技术
[0001]本专利技术涉及地球物理数据处理领域，更具体地涉及广域电磁法接收机采集数据的工频滤波和实时FFT计算。
技术介绍
[0002]传统电磁法勘探主要缺陷有：理论方法上采用平面波近似公式，探测深部信息不准确；探测技术上单次探测数据量少，大量谐波数据被丢弃，勘测效率不高、分辨率较差等；仪器装备上采集信号信噪比不高，转置笨重成本高等。而广域电磁法建立了全新的曲面波电磁勘探理论，以严格的电磁波方程为基础，不做简化，将电流
‑
频率
‑
位置
‑
大地电阻率
‑
磁导率等全部信息考虑在内，只用测量电分量或磁分量中的一个，采用迭代方式计算出广域视电阻率，实现全息电磁法勘探；探测技术上广域电磁法与伪随机电法相结合，实现了一次可发送和接收多个主频信号，能够反映不同深度的地电信息，提高了勘测效率；目前何继善院士团队基于广域电磁法理论，已经设计研发出高精度广域电磁法接收机，测量精度在1％以内，模数转换分辨率达到了24bit。
[0003]目前广域电磁法接收机主控芯片采用的是XILINX的ZYNQ系列芯片，主要包含ARM和FPGA两部分，资源丰富，功能强大。广域电磁法接收基于ZYNQ芯片设计，实现了对广域电磁法电场信号的高精度采集，其主要思路为：传感器输入模拟信号后经过前置放大滤波和噪声压制后，经过A/D高精度模数转换输入到ZYNQ芯片的FPGA端，在FPGA内实现对信号采集的控制以及输入信号的串并转换和缓冲，通过AXI总线传输至ARM的DDR存储器中储存。数据的处理一般将数据读入到PC电脑中进行后续处理。
[0004]FPGA仅仅完成了广域电磁法电场信号的接收控制与传输，不能完全发挥FPGA内部更多逻辑资源的作用。FPGA作为可编程逻辑器件，基于RAM实现，非常适合并行计算。FPGA在数字信号处理领域也有着广泛的应用。相比于电脑上的软件处理，基于FPGA设计广域电磁法电场信号的工频滤波和FFT计算，延迟能够控制在us级，达到实时计算的效果。
技术实现思路
[0005]本专利技术提出一种广域电磁法工频滤波与实时计算技术，解决了广域电磁法电场信号采集数据中工频噪声影响大和不能实时FFT计算的问题。
[0006]为实现上述目的，本专利技术通过以下技术方案予以实现：本专利技术提出一种广域电磁法工频滤波与实时计算技术，包括：基于现场可编程门阵列FPGA芯片实现，由FIR滤波器、巴特沃斯带阻滤波器、插值抽取模块、FFT计算模块组成，包括：
[0007]一套FIR高通滤波器参数和多套FIR低通抗混叠滤波器参数，根据带计算广域电磁法电场信号频率特征，并基于凯瑟窗函数设计，由MATLAB的filter designer工具计算得出滤波器系数，计算出的FIR滤波器系数数符合对称性，减小设计复杂度。
[0008]一套巴特沃斯50Hz工频带阻滤波器参数，在50Hz附近设置滤波器的通带低频截止频率、通带高频截止频率、阻带低频截止频率、阻带高频截止频率、通带最大衰减、阻带最小
衰减等参数，在MATLAB上由buttord算法计算出高阶巴特沃斯带阻滤波器的阶数和截止频率，接着由butter函数求出巴特沃斯带阻滤波器系数，并转换为二阶级联形式，更容易在FPGA上实现。
[0009]插值和抽取模块根据广域电磁法电场信号频率特征和A/D采样芯片采样率特征，设计插值与抽取倍数。为减少信号损失，先进行插值，FIR低通滤波除去镜像后，再进行抽取，再通过FIR低通滤波去除高频噪声混叠。经过插值抽取后的信号变换采样率同样满足奈奎斯特采样定理。
[0010]经过插值和抽取后的广域电磁法电场信号每次计算的FFT点数统一变换为2048点，通过FPGA上的FFT IP核实时计算。
[0011]凯塞窗函数实现FIR滤波器的方法如下：
[0012]理想低通滤波器可表示为：
[0013][0014]其中w
c
为滤波器低通截止频率，α为群延时。(1)式所对应的单位冲激响应为：
[0015][0016]由于上式冲激响应函数无限长，且非因果，需要使用窗函数进行截断。凯塞窗函数具有阻带衰减快，通带起伏小的优势，所以这里采用凯瑟窗函数设计FIR滤波器，凯塞窗函数为：
[0017][0018]其中
[0019][0020]其中I0(x)表示零阶第一类贝塞尔函数。同样在MATLAB上的Filter Desiner工具上设置好滤波器参数后，便能直接求出基于凯塞窗的FIR滤波器系数。
[0021]FIR滤波器在FPGA上实现，需要先求解出其滤波器系数，而MATLAB的Filter Designer工具上，选择基于凯塞窗的FIR低通滤波器，设置好滤波器采样频率38400Hz，可以根据截止频率8500Hz和9500Hz，通带最大衰减3dB，阻带最小衰减50dB等参数设计FIR滤波器并自动计算出滤波器系数，对滤波器系数进行指定16位宽的量化，最后生成coe文件导入到FIR IP核中，在FPGA上实现FIR滤波。
[0022]设计的工频滤波器滤除的频率在50Hz左右，现将受工频影响较大的第2、3频组采样率降低为400Hz，再设计7个参数：滤波器采样频率fs为400Hz,通带低频截止信号频率f
pl
为48Hz、通带高频截止信号频率f
ph
为52Hz、阻带低频截止信号频率f
sl
为49Hz、阻带高频截止信号频率f
sh
为51Hz、通带最大衰减A
p
为2dB、阻带最小衰减A
s
为16dB。将截止信号频率转换为归一化角频率：
[0023][0024]将以上参数代入MATLAB中的buttord函数求出巴特沃斯带阻滤波器的阶数N为8和截止频率w
n
为[0.2421,0.2580]。利用N和w
n
求出巴特沃率零点z、极点p、增益g模式的传递函数H(z)：
[0025][0026]将上式分解为二阶级联形式：
[0027][0028]求出的系数a和b组成一个4行6列矩阵，增益g为0.9370。8阶巴特沃斯带阻滤波器可以拆分为4个二阶巴特沃斯带阻滤波器级联形式，其中每个二阶传递函数可看成是输入信号移位相乘的结果除去递归输出移位相乘的结果，存储系数g，a，b，在FPGA实现二阶巴特沃斯带阻滤波器，通过级联多个二阶巴特沃斯带阻滤波器，实现高阶滤波器工频滤波。
[0029]插值与抽取模块，广域电磁法电场信号频点范围在0.01Hz到8192Hz频带上分布，共分为0频组、1频组、2频组、3频组、4频组、5频组共六个频组，每个频组包含7个频点，其中0频组、2频组、4频组主频成分按2
n
Hz(
‑
6≤n≤13)分布，1频组、3频组、5频组主频成分按0.75*2
n
Hz(
‑
6≤n≤13)分布本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种广域电磁法工频滤波与实时计算技术，其特征在于，包括：FIR高通滤波器参数和FIR低通抗混叠滤波器参数，根据带计算广域电磁法电场信号频率特征，并基于凯瑟窗函数设计，由MATLAB的filter designer工具计算得出滤波器系数，计算出的FIR滤波器系数数符合对称性，减小设计复杂度；基于现场可编程门阵列FPGA芯片实现，由FIR滤波器、巴特沃斯带阻滤波器、插值抽取模块、FFT计算模块组成，巴特沃斯50Hz工频带阻滤波器参数，在50Hz附近设置滤波器的通带低频截止频率、通带高频截止频率、阻带低频截止频率、阻带高频截止频率、通带最大衰减、阻带最小衰减等参数，在MATLAB上由buttord算法计算出高阶巴特沃斯带阻滤波器的阶数和截止频率，由butter函数求出巴特沃斯带阻滤波器系数，并转换为二阶级联形式，更容易在FPGA上实现；插值和抽取模块根据广域电磁法电场信号频率特征和A/D采样芯片采样率特征，设计插值与抽取倍数；为减少信号损失，先进行插值，FIR低通滤波除去镜像后，再进行抽取，再通过FIR低通滤波去除高频噪声混叠；经过插值抽取后的信号变换采样率同样满足奈奎斯特采样定理；经过插值和抽取后的广域电磁法电场信号每次计算的FFT点数统一变换为2048点，通过FPGA上的FFT IP核实时计算。2.根据权利要求1所述的一种广域电磁法工频滤波与实时计算技术，其特征在于，凯塞窗函数实现FIR滤波器的方法如下：理想低通滤波器可表示为：其中w
c
为滤波器低通截止频率，α为群延时；(1)式所对应的单位冲激响应为：由于上式冲激响应函数无限长，且非因果，需要使用窗函数进行截断；凯塞窗函数具有阻带衰减快，通带起伏小的优势，所以这里采用凯瑟窗函数设计FIR滤波器，凯塞窗函数为：其中其中其中I0(x)表示零阶第一类贝塞尔函数；同样在MATLAB上的Filter Desiner工具上设置好滤波器参数后，便能直接求出基于凯塞窗的FIR滤波器系数。3.根据权利要求2所述的一种广域电磁法工频滤波与实时计算技术，其特征在于，FIR滤波器在FPGA上实现，需要先求解出其滤波器系数，而MATLAB的FilterDesigner工具上，选择基于凯塞窗的FIR低通滤波器，设置好滤波器采样频率38400Hz，可以根据截止频率
8500Hz和...
【专利技术属性】
技术研发人员：唐红非，蒋奇云，裴婧，
申请(专利权)人：中南大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。