基于树莓派的VLF数据采集系统

基于树莓派的甚低频数据采集系统开发

摘要

空间研究已成为新时代最具信息量和最受欢迎的领域之一。对空间的监测及其现象的研究早已开始。为了使研究成本更低且更具便携性，斯坦福太阳中心首次提出了名为突然电离层扰动监测器（SID）的空间天气监测项目。该监测系统用于监测影响甚低频（VLF）无线电传播的地球电离层变化。这种低成本但高灵敏度的探测系统已通过各种改进在全球范围内广泛应用。本文提出了一种基于树莓派的监测系统设计与实现方案，具有体积小、功耗低、灵敏度高的特点。树莓派作为基础平台，连接环形天线和前置放大器。数据采集系统采用 Python平台编写，并在Linux系统上运行。最终成功开发出一种高效且更便携的基于树莓派的甚低频数据采集系统，用于检测突然电离层扰动。该系统能够接收并检测来自澳大利亚 NWC（19.8 kHz）和印度VTX3（18.2 kHz）发射的监测频率。

关键词 —甚低频，电离层，Python，突发电离层监测仪，树莓派

一、引言

各种太阳活动，即日冕物质抛射、太阳耀斑、太阳风以及空间中的其他有害事件，引起了研究人员的关注，亟需解决这些可能对宇航员及地表人类造成威胁的现象。当电离层中的电子密度发生变化时，就会发生突然电离层扰动（SID）。电离层中电子密度的变化主要由太阳耀斑期间产生的X射线辐射引起。这一现象会影响甚低频（VLF）波的传播[1]。

斯坦福太阳中心曾提出创建并分发一种低成本的突然电离层扰动（SID）监测系统，用于探测太阳耀斑及地球电离层的其他扰动，此类尝试已取得一定进展。目前，全球大多数人正在使用由斯坦福太阳中心开发的低成本但高灵敏度的空间天气监测设备“超级SID”（Super SID）来进行SID探测[1]。在现代技术和先进设备发展的今天，以往研究中的局限性可以进一步改进，以提高有效性与便携性。十多年来，通过甚低频（VLF）进行空间天气研究已引起研究人员的广泛关注。研究人员相信，空间天气正变得越来越重要。

由于最近对太阳对地球环境及其他通信系统影响的研究，分析所需的区域变得至关重要[2, 3]。

一些典型因素威胁着便携式甚低频数据采集的理念。例如，甚低频无线电波传播通常被视为一个波导问题，其中电离层和地球表面作为边界表面[4, 5]。其他问题，如天空低层结构中的昼夜、季节性和瞬态变化，会在特定位置的射频和电磁场中产生巨大的幅度变化，最终影响通信系统的性能[4]。大量研究表明，射频传播参数的分布在实时、空间和频率的多个区域具有显著差异。因此，由于这些因素动态变化，实现完美的检测并不总是可行的。

利用甚低频获取的数据可以增强对空间及其他大气问题的研究能力。地面科学家可以利用这种增强能力，通过监测远处甚低频发射器的信号强度来探测太阳耀斑。当发生太阳耀斑时，会记录并显示太阳耀斑引起的突然电离层扰动[6‐8]。空间天气是一个广泛的研究领域，尤其是发生甚低频波反射的电离层区域。

使用低成本设备测量甚低频无线电传播是一项巨大的挑战。多年来，研究人员一直使用体积庞大且敏感的空间监测工具，其便携性和有效性较差。该研究领域的许多方面在精度、处理时间、便携性和成本效益方面仍有待提高。因此，本文旨在利用配备环形天线、高质量模数转换器（ADC）和液晶显示器输出的树莓派，构建一种更加便携且高效的甚低频检测系统。首先，设计了用于捕获甚低频信号的环形天线及其前置放大器。接着，使用高质量USB声卡作为模数转换器连接到树莓派。采集的数据在树莓派中进行处理并存储，检测结果通过液晶显示器显示，或通过计算机进一步分析。获取并利用频谱来识别来自附近发射器的甚低频信号。

II. 数据采集系统开发

传统的甚低频数据采集流程如图1所示。对于该频率处理系统，甚低频频率提取标志着甚低频数据采集的开始。该过程包括天线选择、突然电离层扰动监测以及突然电离层扰动分析（即超级SID）的开发。数据记录以及数据的处理与分析均在单板计算机树莓派上进行。

A. 甚低频提取

来自地球各处的无线电信号到达电离层，电离层会反射甚低频无线电信号。反射的无线电信号会随着电离层电离强度的变化而变化。电离层根据高度层次被划分为多个区域，即D层、E层和F层。电子密度或X射线辐射也随高度变化而变化。电离会引起所有传播参数的改变。因此，传播信号的波长变化为3‐30 kHz频率范围及大约70至90 km高度层次的相关研究提供了途径。因此，为了定位空间天气中的这种突然变化，引入了一种名为突然电离层扰动监测的新监测系统。

为了捕获发射器产生的频率，天线是一个关键组件。接收端应能够完美地捕获所需频率而没有任何缺陷。环形天线是一种射频天线，由环路、线圈或其他电气元件组成，通常由平衡源供电或向平衡负载馈电，以在捕获目标频率时获得完全精度。环形天线具有双瓣辐射模式，对相反方向上的两个宽瓣内的无线电波敏感，二者相距180°。因此，它们可以轻松区分方向图以定位发射器的位置。该准则将有助于定位太阳扰动识别。

B. 频率分类器

在检测到突然电离层扰动（SID）和电离现象后，相关信息会被传递给一组分类器。为此目的存在多种分类器模型，例如有限差分频率（FDFD）模型、时域有限差分（FDTD）模型以及最小频移键控（MSK）模型。分类过程更侧重于找出所捕获的甚低频信号的相位与电离层状态之间的关系。FDFD模型描述的是频域有限差分，但通常用于解决散射问题。FDTD模型基于时域，该模型能够通过自然处理非线性材料特性，在一次仿真中覆盖较宽的频率范围。MSK模型通常是连续相位频移键控过程。

C. 数据记录与处理

数据记录机制主要强调对前置放大器捕获的数据进行处理。然后利用先进的软件工具和处理方法，获得实际结果，从而捕捉太阳耀斑事件。通常有一些方法用于此目的，其中超级SID监测软件的效果最为显著。数据处理操作包括数据转换、模拟到数字转换、采样、本地时间转换、峰值检测过程以及存储。在此过程中，使用声卡和树莓派作为数据记录器。该项目的采样率设置为48千赫，因此可获取的最大频率为24千赫。

III. 设计与实现

该项目分为硬件和软件开发两个主要部分。在硬件部分，开发并使用了环形天线、前置放大器、USB声卡和液晶显示器。在软件方面，已在树莓派上安装用Python编写的SuperSID软件[9]。此外，还针对所使用的条件和特定硬件对其进行了适配和修改。

图2显示了所提出系统的框图。首先，环形天线在接收端接收甚低频信号。接着，前置放大器将对接收到的信号进行放大。此外，模数转换器将把模拟信号转换为数字信号，经过处理后保存到SD卡中。最后，处理后的数据将在液晶显示器屏幕上显示。

A. 环形天线设计

该项目的主要标准之一是使用环形天线捕获来自多个站点的甚低频信号。为了获得有效结果并实现成本效益高，天线设计的主体结构采用PVC管。框架使用22或更细的AWG #26导线缠绕，尺寸为半米到一米方形空间。此外，使用RG58同轴电缆通过接线端子将前置放大器与天线连接。完整的环形天线设计如图3所示。

B. 前置放大器设计

在本研究中，前置放大器将引导环形天线捕获的频率信号，并将其直接存储到与树莓派相连的SD卡中。该前置放大器的设计灵感来源于斯坦福太阳中心在此领域获取数据的首次尝试。该设计采用一块单面条状电路板（10厘米×12.5厘米），由于电路板尺寸较小，选用了表面贴装器件（SMD）进行印刷电路板的焊接。使用函数发生器提供2‐200千赫兹的频率信号，并保持20毫伏的幅度。通过3.5毫米TRS麦克风接口将前置放大器连接至与树莓派相连的声卡。整个电路由一个9伏适配器供电。电路原理图及所搭建的电路如图4所示。

C. 树莓派

作为本研究基础的单板计算机树莓派负责运行所有的软件操作。Raspbian Wheezy操作系统将运行所有分析软件以检测频率峰值。使用了一张Class 10 SD卡，容量为32GB，碎片化顺序写入速度为10MB/s[10, 11]。

项目选用了一块高质量的声卡，采样率为44100千赫兹。通过HDMI连接，使用液晶显示器在屏幕上显示必要信息。树莓派连接了无线网络，并通过5伏适配器为单板计算机供电。完整的实验装置如图5所示。

D. 软件设计

超级SID监测软件是分析太阳耀斑事件的主要软件。该开源软件可从GitHub网站[9]下载。由于该软件基于Python语言编写，因此必须在树莓派上安装一些模块才能正常运行SID软件[11, 12]。这些模块包括AlsaAudio、wxPython（wxgtk2.8）、NumPy和Matplotlib。

Matplotlib模块用于绘制SID数据，而AlsaAudio模块使用高级Linux声音架构来录制和回放保存的音频信号。另一方面，NumPy模块是一个用于控制大规模多维数组集的通用数组处理包。对于图形用户界面，我们使用了wxPython模块工具包为本项目创建了一个功能强大的图形界面。

在执行SuperSID软件程序之前，需要正确配置配置文件（supersid.cfg），以显示来自目标站点的甚低频指示器。通过写入命令“viewer=wx”来启用图形用户界面（GUI）模式。与通过命令“viewer = textual content”启用的文本模式相比，GUI模式会从树莓派调用更多资源。在（sampler.py）和（supersid.cfg）文件中的命令行“audio_sampling_rate”应根据项目中所用声卡支持的采样频率更改为48000赫兹。采集到的记录可保存在预配置的“data_path”目录中。USB声卡已在配置文件中使用“Card = Pro”输入进行添加。最后，选择了两个发射器，包括NWC（澳大利亚）和VTX3（印度），作为主要发射器位置。随后，接收到的信号被转换为本地时间，并通过软件分析给出绝对峰值频率，以检测大气扰动。

IV. 实验结果与分析

根据软件参数，VTX3（印度）在18200赫兹、NWC（澳大利亚）在19800赫兹的发射器提供了垂直红线以指示截止频率水平。具体数据采集第一天从2:00至5:00持续运行6小时，第二天再运行2小时。

图6显示了其y轴上的信号强度，其中展现了突然的峰值。以同样的方式，进行了每小时模拟，然后在编译框架中与其他文件合并，以了解来自发射站的信号强度，因为图表中的峰值表明了改进的结果。

这六小时的数据编译后自动生成了“.csv”文件。现在通过使用plot.py软件，数据采集结果如图7所示。编译图中的24小时时长明确表明，除了特定的几个小时外，并未捕获到其他时间的数据，因此数值保持为0。另一方面，两个不同站点根据其信号强度产生了不同的峰值。然而，图8显示了从原型数据采集中获得的合并的编译文件具有更好的性能。

在仿真过程的初期，从原型设备获取的数据显示出信号强度的最小范围处于负值区域。然而，随着流程改进的推进，结果得到了进一步改善和优化。图8清晰地显示了系统在早期数据采集阶段的极低信号强度。该图还表明，在8:00至10:00之间，VTX3站传输的信号强度最高达到了40分贝。尽管如此，在仿真过程中，低强度信号水平的峰值更为明显，因为数据采集环境发射出低频噪声，这些噪声反映在图中的其他部分。图9展示了实时采集的一个示例。

图10中绘制的低频部分表示低频峰值。监测系统在5秒时间帧内捕捉到了一次突然出现的电离层信号的信号足迹，在接近18200赫兹处产生了一个突然的峰值。此外，监测系统同时也在接收19800赫兹的发射器信号，声卡在树莓派系统中捕获了发射的数据。图11显示了接近发射频率处高频强度的突然峰值。

收集的原始数据来自空间研究站。本次操作选定的频率为18200赫兹和19800赫兹。另一方面，软件和硬件均进行了增强，以检测所需数据，通过在白天两小时的时间边界内出现明显尖峰来确认太阳耀斑事件。然而，由于整个过程是在较短的时间安排和封闭环境中完成的，因此未能生成精确的结果。其主要原因是国际伊斯兰大学马来西亚分校的仪器与测量实验室中，所有其他电气设备同时运行，导致监测系统受到噪声和干扰。

此外，如果环形天线能够架设在高于地面的开阔区域，将能够接收到更精确的信号。另外，需要更长时间的数据采集过程以提高数据的精度和有效性。在我们的仿真试验中，大多数日子都是雨天和阴天天气。因此，将会更加如果能够在晴朗的白天运行模拟，从而无缝检测到突然的电离，则效率更高。

V. 结论与未来工作

本文介绍了使用树莓派进行甚低频数据采集的设计与实现。首先，设计了环形天线以有效捕获发射器频率。然后，设计了前置放大器电路，并将其连接到天线和高质量的USB声卡。数字化信号随后通过配备相关软件库的树莓派进行处理、记录和显示。实验结果表明，所提出的系统能够检测大气扰动。然而，还需在不同天气条件下进行更多实验以及更长时间的实验，以捕捉各种突发性电离层扰动。