基金项目:国家重点研发计划(2019YFC1509502).
第一作者简介:陈卓琳(1998-),硕士研究生在读,主要从事地球探测与信息技术研究.E-mail:chenzl1692948775@163.com.
通信作者简介:胡星星(1978-),研究员,主要从事地磁观测技术、高精度大动态范围地震数据采集技术、MEMS地震传感技术、结构健康监测等方面研究.E-mail:huxx@cea-igp.ac.cn.
(1.中国地震局地球物理研究所,北京 100081; 2.中国地震台网中心,北京 100045; 3.苏州工业职业技术学院,江苏 苏州 215000)
(1.Institute of Geophysics,China Earthquake Administration,Beijing 100081,China;2.China Earthquake Networks Center,Beijing 100045,China;3.Suzhou Institute of Industrial Technology,Suzhou 21500,Jiangsu,China)
DOI: 10.20015/j.cnki.ISSN1000-0666.2025.0016
研究地表处的磁场有助于地震学家了解孕震区地下介质电导率在地震发生前的变化。地震预报专家基于震磁异常效应的预测分析方法,在以往许多中强震前给出了有效的预测意见(倪晓寅,洪旭瑜,2022; 翟世龙等,2024)。因此,开展高分辨率的地磁场观测在地震科学研究和地震监测预报等方面具有重要应用价值。无人机磁测是近年来迅速发展起来的一种地磁观测方式(崔志强等,2016)。通过无人机搭载磁测系统,可以探测地磁监测空白区的磁场变化情况,能够在很大程度上弥补地面地磁台站观测的不足(王猛等,2022; 张永旺等,2022)。但由于无人机机体铁磁性材料及电子设备的干扰,无人机矢量磁测数据经姿态校正后仍无法与地面台站观测数据质量相比(郭华等,2022; 杜俊岐,2021),因此,对机体干扰磁场补偿是提高航磁测量数据质量的关键(高全明,2020)。现有的采用多传感器的差分补偿方法,即在机头、两侧机翼分别放置一个传感器,能够克服机体自身剩磁干扰(朱学军,2014; 赵琳等,2018),提高磁测数据质量。传感器一般采用低功耗、小体积的磁通门矢量传感器,而数据采集系统需要根据实际应用需求进行定制化研制。为了采集低噪声、高精度的无人机磁测数据,本文针对3个磁通门传感器补偿系统,研制了能够接入3个三分量传感器的多通道数据同步采集系统。该系统由模拟采集和控制两部分功能模块组成,分开设计成相应的两块体积较小的电路板,以便适应空间有限的安装应用环境。其中,模拟采集部分采用高性能低功耗的32位微处理器STM32F103C8T6,完成9通道数据的同步采集、温度数据和倾角数据的获取,可满足无人机磁测系统的重量要求和姿态校正需要(Funaki et al,2014); 控制部分基于ARM+LINUX嵌入式系统为程序运行平台,实现GNSS精确授时和本地大容量数据存储及网络传输接口,方便存储数据文件的下载和后接无线模块进行实时传输。该采集系统体积小、功耗低,能够满足中小型无人机对载荷重量和搭载空间有限的需求(周普志等,2022)。
本系统包括模拟采集部分和控制部分。系统框图如图1所示。模拟采集部分需要实现的功能有:将3个三分量磁通门传感器的输入模拟信号连续转换为24位高分辨率的数字信号,通过SPI串口同步发送给微处理器STM32,微处理器STM32将9通道AD数据、倾角数据和温度数据打包并按照规定的通信协议发送给控制部分。控制部分需要实现的功能有:①与模拟采集部分通信,对原始采集数据进行解析、格式转换、整理等得到最终的数据结果,并进行本地存储管理; ②实现远程FTP存储数据下载和传输等; ③为系统的实时时钟提供时钟驱动,并完成GNSS的同步校时和授时功能; ④本系统提供稳定、低噪声的供电电压。
模拟采集部分的电压由控制部分提供,输入电压为12 V左右,首先经过稳压芯片WRB1206S-3WR2降为6 V,其体积小、效率高,拥有可持续短路保护功能,电源模块电路图如图2所示。由图2可见,输入电压又可分为两路,一路经功耗较低的稳压芯片K7803MT-1000R4转为5 V供给模拟采集电路板中的数字部分,另一路经低压差线性稳压器LDO转为5 V后又经高精度低温漂低噪声电压源ADR4525CRZ降低至2.5 V,供给模拟采集电路板中的模拟部分。为避免模拟芯片受到高频噪声干扰,对模拟部分和数字部分进行了隔离供电,并采用一点接地的方式。
多传感器差分补偿的实现要求传感器的数据采集同步。为此,9通道AD芯片采用同一高精度外部时钟驱动,能够使得各通道的采集具有严格同步的时钟周期,模块设计如图3所示。
为满足传感器输出信号的大动态范围,本系统采用24位高精度模数转换芯片ADS1255。这是一款内带四阶ΣΔ调制器的高精度24位AD,具有可达23位无噪声分辨率、最大非线性±0.001%的优良特性,且外围电路非常简洁,有利于提高系统的可靠性和信噪比。采集部分采用低功耗单片机处理器STM32F103C8T6,其与9个模数转换芯片ADS1255之间采用SPI高速串行通讯接口,各通道AD转换完成后输出的/DRDY低电平中断信号经或门电路合并为一路,作为采集处理器读取AD转换数据的总中断信号,因而能够实现9通道AD转换数据的同步读取。
采集数据的绝对时间戳的准确性由控制部分的上位机采用GNSS校时实现。上位机接收到采集部分的实时数据后,通过读取系统RTC获得绝对时间。RTC时间由GNSS接收模块进行校时,并为保证校时精确性,系统采用GNSS模块的秒脉冲信号PPS作为校时的中断信号。
采集程序首先对串口等外围接口及外设进行初始化,配置中断优先级、设置采集参数等。采集部分控制器收到采集中断信号后,进入中断服务函数,读取转换结果。
模拟采集和控制部分之间采用串口通信,为保证传输数据的完整性,设计通信协议为4 Byte帧头+48 Byte有效数据+4 Byte帧尾,数据的接收和发送均按照通信协议进行,如图4所示。
为满足控制部分对系统的软、硬件要求,采用低功耗高性能的32位嵌入式ARM处理器S3C2416作为主控制器,以裁剪后的嵌入式LINUX操作系统作为程序运行平台。采用RS232串口与模拟采集部分进行通信,外接32 GB存储卡作为采集数据存储器; 采用10/100 M自适应网卡作为对外传输通信,并用以查询和设置采集器参数,如图5所示。
数据存储模块采用容量为32 GB的TF卡,其体积小、容量大、存取速度快。在9通道50 Hz采样时,可存储大约5个月的连续采集数据。控制器采用SDIO(安全数字输入、输出)接口连接TF卡。
GNSS模块采用灵敏度高、体积小、功耗低的NEO-6M芯片,其内置NMEA-0183标准协议。NMEA-0183协议是为了在不同的GNSS导航设备中建立统一的RTCM标准。该协议定义了多种命令,最常用的6种命令分别为$GPGGA、$GPGLL、$GPGSA、$GPGSV、$GPRMC和$GPVTG,其中$GPRMC包含了最小定位信息,如时间、经纬度等,满足本采集系统校时和定位的设计需求,因此本模块仅解析该命令。
控制部分程序基于LINUX系统设计,分为应用层和驱动层。驱动层通过响应硬件中断,实现实时处理。应用层主要负责数据的接收和处理,并响应各类控制命令,实现人机接口。
驱动程序设计如图6所示。设备驱动程序响应硬件请求,是内核与硬件之间的接口。为使采集数据有精确的采样时刻,采集系统在驱动程序中通过中断方式响应来自GNSS模块的秒脉冲(Pulse Per Second,PPS)信号,并通过SIGNAL信号机制通知应用程序。
应用程序完成数据采集、数据校时、数据存储、数据处理和网络传输通讯等任务,主要有实时数据接收模块、GNSS高精度时间服务模块、本地数据存储模块、远程数据服务模块,如图7所示。为了实现多任务的“并行”处理,应用程序采用多线程设计。在采集器应用程序设计中,以全局数组变量作为公共的数据缓冲池,并创建4个处理线程:线程1是程序主线程,循环监控网络连接请求。当有远程用户试图连接时,主线程创建一个新的线程去处理,然后继续进行循坏监控; 线程2用于与下位机通信获取采集数据,同时处理GNSS秒响应信号,为采集数据标上精确时间戳; 线程3用于数据的本地存储管理,根据采集器的设置要求按每小时的时间长度保存为连续数据文件,并监控存储空间容量变化情况,根据设置是否要删除过期数据文件等处理; 线程4对采集的原始数据进行解析处理,如数据格式转换等。
设置该数据采集系统采样率为50 SPS,使用高精度直流电源输出直流电压信号,测试信号幅度分别选择为-5 V、-4 V、-3 V、-2 V、-1 V、0 V、1 V、2 V、3 V、4 V、5 V,每个测试点记录30 s长的数据,线性测试数据记录见表1。由最小二乘法(地震观测仪器进网技术要求 地震烈度仪,DB/T 59—2015)计算得到线性度误差为5.56×10-6,线性测试结果如图8所示。
采用信号发生器输出正弦信号作为采集系统的测试信号。测试时,使输入正弦测试信号幅值保持2 V(峰值)不变,输入频率为0.01~25 Hz,采集系统采样率为50 SPS,记录采集系统的输出并通过快速傅里叶变换计算响应信号幅值。以0.01 Hz时的输出信号幅度作为参考,输出信号幅值降至0.707倍即-3 dB时的输入信号频率即为该采集系统的响应截止频率。幅频特性测试记录见表2,表中分贝数计算如下:
dB=20×lg(Ni/N0) (1)
式中:Ni为测量值,i= 0、1、2……24。
幅频特性测试结果如图9所示,结合图中数据可以看出,频带范围为DC~22 Hz。
在北京国家地球观象台的地磁观测室内进行数采连接磁通门磁力仪的观测实验。在数采的输
入端接入3台磁通门磁力仪,等待磁通门磁力仪运行稳定后,选取一天的测试数据(时间为国际标准时)与台站仪器记录的数据进行对比分析,数采所记录的3台磁通门磁力仪测试数据与台站磁力仪观测数据经去平均值后其日变曲线如图 10所示。从图中可以看出,3台测试仪器与台站仪器记录到的磁日变曲线大致重合。引入皮尔逊互相关系数进一步分析测试数据的一致性,皮尔逊法则是一种经典的相关系数计算方法,主要用于表征线性相关性能,皮尔逊相关系数为-1~1,绝对值越接近于1,则表明两个变量的相关程度越高(刘若男等,2022)。因此,该方法常用于比较地磁矢量观测数据的一致性(刘成等,2019; 周思远等,2023)。3台磁力仪测试数据与台站磁力仪观测数据的皮尔逊互相关系数分析见表3,从表中可以看出,3台测试仪器之间以及3台测试仪器与台站仪器之间的数据互相关系数均为0.9以上,表明其一致性较高。
图 10 3台磁通门磁力仪(50 Hz)与台站磁力仪(1 Hz)日变曲线记录
Fig.10 Diurnal curves of three fluxgate magnetometers(50 Hz)and the station magnetometer(1 Hz)
控制部分单独工作时,功耗为1.3 W。连接上模拟采集部分,等待整个采集系统稳定之后测得的功耗为4.5 W。
基于嵌入式技术设计的多通道数据采集系统,采用模拟采集和控制部分分体式的设计,模拟采集电路板大小为18 cm×11 cm,控制部分电路板大小为11 cm×8 cm,较小的尺寸使得其在空间有限的无人机上的安装和固定更具灵活性,并可以使易受干扰的传感器精密模拟信号处理系统远离噪声较大的数字系统。整个系统稳定工作后的功耗仅为4.5 W,较低的功耗可以提高采集效率,延长续航时间。对于地面地磁观测,1 SPS的采样率已能满足观测准静态地磁场的要求。对于航磁,为了获得较高的空间分辨率则需要提高采样率。在不高于100 km/h的航速时,采用50 SPS的采样率能实现每米距离得到1个采样点的精度。幅频特性测试结果表明频带满足采样率的要求。较高采样率对于通道间的同步提出了更高的要求。本文对模拟采集部分的9个AD芯片使用同一外部时钟驱动,其采集信号经或门电路合为一路作为中断信号,使处理器可以同步采集9个AD芯片的数据,保证了通道间数据具有高精度的相对同步。控制部分的GNSS模块为整个系统提供时间信息和PPS秒信号,从而保证本地时钟具有高精度绝对准确性。从静态测试结果来看,该数据采集系统能够满足多通道航磁测量的需求。但是在飞行情况下,该系统的采样精度、稳定性等是否会受影响,无人机自身以及环境的磁干扰是否会降低采集质量等问题都需要开展无人机动态测试验证,但因北京地区飞控管制以及无人机设备借用不便等问题,机载动态测试目前尚未进行。此外,动态测试需要选择机身材料磁屏蔽能力较高、续航时间长、可控性较强的无人机,在磁干扰较小的环境内以重复路线多次飞行的方式对系统的噪声、稳定性等进行测试。
本文针对多传感器磁测系统,研制了一套可适用于无人机的多通道航磁同步数据采集系统,该系统基于嵌入式LINUX系统,提供具有高精度的本地时钟,并通过同一高精度外部时钟驱动9通道AD同步采集模块,设置系统采样率为50 Hz以提高采集精度。对该系统进行测试,主要得到以下结论:
(1)该数据采集系统体积小、质量轻、功耗低,适用于无人机。
(2)该系统在50 SPS的采样率下,线性度误差为5.56×10-6,-3 dB频带范围为DC~22 Hz,满足航磁的高分辨率采集需求。
(3)3个磁力仪之间以及3个磁力仪与台站比对仪器之间的采集数据一致性良好,表明该数据采集系统各通道之间的同步性较好、数据质量较高。