1. GReX仪器设计概述
GReX（Galactic Radio Explorer）是一套专为探测银河系内快速射电暴（FRBs）设计的全天空监测系统。作为STARE2项目的升级版，GReX通过全球部署的低成本终端网络实现了全天候、高灵敏度的射电监测能力。整套系统的核心创新在于其模块化设计——每个终端都集成在单个防水铝制箱体内，包含完整的信号接收和处理链。
关键设计理念：牺牲角分辨率换取视场覆盖，通过分布式部署实现全天区监测。这种"广撒网"策略特别适合捕捉随机出现的超亮瞬变事件。
仪器工作频段为1280-1530MHz（L波段），这个频段选择基于两点考量：一是避开地面通信的主要干扰频段；二是包含中性氢线（1420MHz）这一重要的校准源。系统时间分辨率达到32.768μs，并能存储原始8.192μs采样数据，为后续分析保留了充分的信息量。
2. 硬件架构解析
2.1 射频前端设计
射频信号链的起点是改进型"蛋糕盘"波导喇叭天线。这种轴向波纹喇叭天线具有70°±5°的波束宽度（1420MHz时），其宽波束特性虽然会降低点源灵敏度，但大幅提高了瞬时视场覆盖率。实测表明，该天线在1280-1530MHz范围内波束图案变化小于5%，确保了宽频带性能的一致性。
低噪声放大器（LNA）采用Caltech为DSA-110项目开发的新型设计：
- 噪声温度：7K（旧版30K）
- 采用Diramics超低噪晶体管
- 悬浮带状线输入匹配网络
- 工作带宽：1280-1530MHz
- 增益：40dB
这种LNA的创新之处在于在常温下实现了通常需要低温冷却才能达到的噪声性能。我们通过Y因子测量法验证，实际系统噪声温度控制在50K以下，比STARE2提高了约4倍灵敏度。
2.2 前端模块(FEM)
FEM是模拟信号处理的核心，集成在单个屏蔽盒内，包含：
- 下变频器：将RF信号转换到307.2MHz中频
- 带通滤波器：1280-1530MHz→250MHz带宽
- 可变衰减器：0-31.5dB，1dB步进
- 监控系统：基于RP2040微控制器
- 实时监测LNA工作状态
- 温度传感器（±0.5℃精度）
- 平均功率检测
- 串口控制接口
整个FEM采用Rust语言编写固件，使用RTIC框架实现实时任务调度。模块化设计使得单位成本控制在200美元以内，而性能指标媲美专业射电天文设备。
2.3 数字化处理单元
信号数字化由SNAP（Smart Network ADC Processor）板卡完成：
- ADC：HMCAD1511（3片）
- 采样率：500MSPS
- 分辨率：8bit
- FPGA：Xilinx Kintex-7 160T
- 数据传输：10G以太网（SFP+）
FPGA固件采用CASPER工具流开发，核心处理流程包括：
- 多相滤波器组(PFB)
- 8抽头FIR滤波器（Hamming窗）
- 2048点FFT
- 输出18+18bit复数
- 重量化：18bit→8bit
- 数据打包：8200字节/包（含64bit时间戳）
这种设计在保留信号特征的同时，将数据率从原始的8Gb/s（18bit×2048通道×500MSPS）压缩到1GB/s，适合通过10G网络传输。
3. 软件系统架构
3.1 实时处理流水线
T0程序构成数据处理的第一级：
- 线程1：网络包捕获
- 检查包序号连续性
- 统计丢包率（典型值10^-6%）
- 线程2：测试脉冲注入
- 插入模拟FRB信号
- SQLite记录注入参数
- 线程3：时间降采样
- 32μs→1ms
- 计算Stokes I参数
- 线程4：数据输出
- SIGPROC格式滤波器组文件
- PSRDADA环形缓冲区
程序采用Rust编写，利用线程间通道实现无锁通信。监控子系统通过Prometheus收集性能指标，OpenTelemetry处理日志，形成完整的可观测性体系。
3.2 射频干扰抑制
clean_rfi程序实现三级RFI消除：
- 静态掩码：标记已知干扰频点
- 零值过滤：丢弃异常采样点
- 带通校正：除以时间平均谱
- 方差截断：
- 第一轮：3σ剔除
- 第二轮：5σ剔除
实测表明，这套方法能有效抑制手机基站、雷达等常见干扰源，同时保留真实的射电瞬变信号特征。
3.3 实时检测算法
改进版HEIMDALL实现暴力去色散搜索：
- 去色散范围：0-5000 pc/cm³
- 盒car滤波宽度：1-256ms
- 输出候选体通过Socket传输
- 触发原始电压数据存盘（NetCDF格式）
检测流程后接T2/T3程序进行候选体聚类和可视化，最终结果通过Slack实时通知观测人员。整套系统采用bash脚本编排，支持多种运行模式切换。
4. 系统部署与性能
4.1 全球站点布局
目前已部署5个终端：
- 欧文斯谷（美国）：替换原STARE2系统
- 康奈尔大学（美国）：屋顶安装
- 哈佛大学（美国）：天文台楼顶
- Hat Creek（美国）：与Allen Telescope Array共址
- Birr（爱尔兰）：历史天文台遗址
每个站点安装流程标准化：
- 天线组装（30分钟）
- 光纤熔接（需专业设备）
- 系统加电测试（2小时）
- Y因子测量（1小时）
站点选择考虑四要素：
- 电力/网络接入
- 物理安全性
- 天空可视范围
- 射频环境质量
4.2 实测性能指标
通过氢线测量验证系统灵敏度：
- 系统等效通量密度：3.8kJy
- 最小可检测流量密度（5σ）：
- 32μs时标：580Jy
- 1ms时标：104Jy
- 动态范围：>60dB
这些指标意味着GReX能探测到类似FRB 200428（峰值约1.5MJy）的银河系内爆发，但对更弱的extragalactic FRBs灵敏度不足——这正是设计时的权衡选择。
5. 科学成果与银河FRB上限
基于2020年4月（FRB 200428探测）至2025年5月的观测数据，我们更新了银河系FRB发生率的上限：
| 时间窗口 | 曝光量(终端·天) | 上限(90%置信度) |
|---|---|---|
| 1天 | 892 | <3.4×10⁻³/天 |
| 1年 | 325,580 | <1.1/年 |
计算采用泊松统计：
上限 = -ln(1-CL)/(ε×T)其中CL=0.9（置信度），ε=0.7（检测效率），T为有效曝光时间。
这一结果暗示：
- 银河系磁星产生FRB 200428级别爆发的概率极低
- 需要更大规模的监测网络或更长时基观测
- 可能存在爆发活动的周期性或集群性
6. 技术挑战与解决方案
6.1 时间同步问题
初期遇到GPS 1PPS信号抖动（~100ns）导致时间戳错误。解决方案：
- 改用温补晶振(TCXO)作为基准
- FPGA固件增加数字锁相环(DPLL)
- 软件端实现NTP微调（<1ms精度）
6.2 数据丢包优化
10G网络持续传输时出现随机丢包。通过以下措施将丢包率降至<10⁻⁶：
- 调整网卡DMA缓冲区（增至128MB）
- 禁用CPU节能模式
- 设置实时进程优先级
- 采用DPDK用户态驱动（备用方案）
6.3 射频干扰抑制
城市站点受WiFi/5G干扰严重。除常规RFI消除外，还采取：
- 物理屏蔽：导电橡胶垫圈
- 空间滤波：利用多站点符合检测
- 频域标记：建立干扰数据库
7. 扩展应用与未来计划
7.1 技术签名搜索
GReX可用于搜寻地外文明信号：
- 窄带搜索（δν=1Hz）
- 最小可检测EIRP：1.2×10¹² W/Hz（比邻星）
- 人工色散信号（τ=1ms）
- 最小可检测EIRP：4.9×10¹² W
通过多站符合检测可有效排除地面干扰，相关算法已集成到SPANDAK处理管线中。
7.2 系统升级路线
下一步改进方向：
- 增加部署密度（目标20个站点）
- 开发第二代前端：
- 双极化接收
- 数字波束形成
- 光学频率梳校准
- 机器学习实时分类：
- 基于ResNet-50的候选体筛选
- 异常检测算法找新现象
当前系统已在GitHub开源（GReX-Telescope组织），包含完整的硬件设计图、FPGA固件和数据处理软件，方便研究机构复现。这种开放科学模式显著降低了射电天文的研究门槛，已有三所大学基于我们的设计建立了教学用射电望远镜。
