The Status of Ion, Gas and Other Items of Groundwater Observation in Seismological Geochemistry
-
摘要: 地震地球化学离子、气体成分等其它观测项目是地震地下流体观测网的重要组成部分,在地震监测预报工作中发挥了重要作用。本文对观测项目的规模、观测数据质量和数据预报应用3方面进行分析,梳理了地球化学离子、气体成分等观测项目的背景变化规律;指出近年来观测规模急剧缩小、观测数据评价体系不完善、台网管理缺失严重等现状;并提出需清理观测项目、完善仪器更新机制、加强台网管理和建立人员培训机制等建议。Abstract: Ion, gas and other items of groundwater observation are important parts of seismic underground fluid observation network, which play the important roles in earthquake monitoring and prediction. In this article we review the three aspects of the scale, the quality of observation data and the application of earthquake monitoring and prediction in China. The background dynamic changes of ion, gas and other items is also analyzed. Our results show some problems in groundwater observation in recent years such as insufficient projects, lack of data quality evaluation system and lack of effective management of network. In the future, we need to reorganize the observation items better, to improve the instrument renewal mechanism, to strengthen the management of the networks, and to set up personnel training mechanism.
-
Key words:
- Ion /
- Gas /
- Seismological geochemistry
-
引言
力平衡加速度传感器和强震动数据采集器是进行强震动观测的2种重要仪器。通常2种仪器各自独立,但是部分厂商已经开始将其合成为1种仪器,即将强震动数据采集器嵌入到力平衡加速度传感器内部,使得2种独立仪器合成为强震动测量、数据转换及记录一体机。这种强震动测量一体机具有体积小、便于携带和联网等特点,尤其适用于流动观测。同时,为了保证获得高质量的强震动观测数据,要求强震动数据采集器的信噪比和动态范围足够大、采样率高且具有足够的存储空间;随着高速、大容量存储介质的出现,高频率采样成为强震仪器的发展趋势。本文基于传统三分量力平衡加速度传感器技术,设计了三分量一体化加速度传感器机械底座,并以低噪音宽动态范围的信号调制解调电路完成加速度传感器模拟电压信号输出;在传感器内部嵌入以24位Δ-Σ型模数转换器ADS1294为基础的高精度模数转换电路,嵌入扩展多种接口的1.4GHz×4微控制器主板,形成1款具有高精度模数转换、具备串口、USB口、网口等多种接口的低功耗、24位数字输出的力平衡加速度传感器。
该数字化力平衡加速度传感器共包含4部分,分别为机械振动传感器底座、模拟信号调制解调电路、三通道同步高精度模数转换电路和高速ARM微控制器主板及数据处理电路,其基本结构见图 1。其中前2部分是传统模拟式力平衡加速度传感器,而后2部分组成了强震动数据采集器或记录器。其中,机械振动传感器底座负责对震动信号进行检测,模拟信号调制解调电路部分负责将传感器机械底座检测到的电荷变化转变为便于测量的电压信号,三通道同步高精度模数转换电路部分完成三通道信号的同步模数转换、数据处理,并提供精准时钟及时钟校正功能;高速ARM微控制器主板完成模数转换电路的控制与数据通讯、数据存储、RS232串口、USB口、网络、电源管理模块等多接口操作。
1. 传感器原理及设计
传感器机械振动底座负责对振动信号进行检测。一般单分量的力平衡加速度传感器由磁缸、磁铁、反馈线圈和中间极板、上极板、下极板等组成部件(李彩华等,2014)。反馈线圈与中间极板共同构成传感器的可动部件,其余部分均为不可动部件。3个正交的单分量力平衡加速度传感器组成三分向力平衡加速度传感器底座。本文中的传感器机械底座具有一体化加工的三分量正交磁缸,采用先进的数字加工技术,在1块纯铁上同时加工出3个垂直正交的机械磁缸,减小了正交误差,解决了传统分体组装传感器横向灵敏度较大的问题,同时缩小传感器的体积。并且通过增强传感器内部磁场强度、调整传感器质量块重量方式将传感器测量频带范围扩展到0—200Hz。
传感器中调制解调电路部分负责将传感器机械底座的电荷变化转变为便于测量的电压信号。此部分电路主要包含桥式振荡电路、解调电路、放大输出电路和反馈电路。桥式振荡电路自激产生一定频率的正弦波,分别被处理成幅值相同、极性相反的正弦波信号后被分别施加到机械底座的上、下极板,因此中间极板输出信号由低频信号调整成了高频信号。该高频信号经由场效应管组成的解调电路后,代表地面振动的电压信号被提取出来,最后由运算放大电路调理成模拟差动电压信号进行输出。为了调整传感器频带特性,输出电压信号经过反馈电路将重新输入到传感器线圈中。
2. 高精度数据采集器原理及设计
高精度数据采集器电路负责将传感器输出的3个模拟电压信号转换为数字信号,完成计算处理操作后,由本地存储介质存储,经串口、USB口、网口等通讯端口发送出去(吴海超,2013;Hu等,2014)。本文的采集器电路从功能上主要分为高精度数据采集电路和高速ARM核心处理电路。
2.1 高精度数据采集电路设计
高精度数据采集电路板主要完成三通道数据的同步模数转换、时钟读取,其基本设计原理见图 2。主要分为前置调理电路、模数转换电路、低速ARM处理器电路和时钟管理电路4部分。
2.1.1 前置调理电路
强震信号尤其远场强震信号是极其微弱的低频信号,其有用频率在0.1—20Hz频带内,从抑制基线漂移和带外噪声以及保护后续电路等方面考虑,在采集器前置处理电路中需要对输入信号进行滤波、调理及限幅。本文设计的前置调理电路为3路完全相同的差动信号调理电路,将加速度传感器的输出电压信号调理成与模数转换器ADS1294输入范围一致的电压信号。图 3为单通道前置调理电路,其中电阻对R1与R3、R2与R4分别完成正负向电压信号的比例分压,阻容电路R5与C2、R6与C3完成正负向电压信号的低通滤波,并经由双二极管BAT54S并将AINP、AINN端输出电压信号钳制在-2.5—2.5V范围内,将后续模数转换电路输入电压钳制在合理输入范围内,可以避免后续电路遭高电压破坏。
2.1.2 模数转换电路
电路选用四通道的同步24位Δ-Σ模数转换器ADS1294。该元件内置多通道可变增益差分输入放大器,还具有电源参考源和晶振。该芯片每个通道的最高采样率为32ksps,是1款基于SPI接口的可编程芯片,通过内部寄存器可以选择内部电路开关,设置增益、采样率等,其24位A/D转换精度可以达到μV级,完全满足强震动观测仪器数据采集的需要。高分辨率模式下功耗为10.1mW,低功耗模式下为8.3mW。在采样率不超过8ksps时仍能具有24位的转换精度。内部集成了4个程控差分输入放大器、Wilson电阻网络等,结合高精度的模数转换,将基线漂移、工频干扰等噪声处理设置在模数转换后的数字处理部分,可以大大简化数据采集前端的设计,实现强震信号的采集。
2.1.3 低速ARM控制器及时钟管理电路
为了完成500Hz、1000Hz等高频率采样时大量数据的处理工作,本研究在该模数转换电路板上内置了低速ARM芯片STM32F407,由该芯片及其相关电路与GPS模块、高精度晶振电路、高速ARM核心微控制器主板共同配合完成模数转换电路的各种时序控制。其中GPS模块产生3个控制信号,即GPS_TX、GPS_RX和TIMEPUSLE,分别实现GPS信息发送、接收和GPS秒脉冲传送,并经通用数据接口PB5口输入到STM32F407内部参与逻辑控制和时序控制。高精度晶振电路负责生成高精度时钟信号,由STM32F407根据GPS秒脉冲数据对高精度晶振进行时钟校准。具体实现方法为STM32F407同时读入GPS秒脉冲信号TIMEPUSLE、高精度时钟晶振输出信号OSC_IN,计算两者差值,并根据差值大小由STM32F407通用数据输出口输出电压调节信号DAC到运放OPA333的正向输入端,该运放的输出端信号再连接到高精度晶振的电压控制端进行调压,实现调节高精度晶振的输出信号OSC_IN的频率偏差,从而利用GPS秒脉冲信号完成了校正实时时钟功能。
2.2 高速ARM核心电路
高速ARM核心电路板需要完成与高精度数据采集电路板的通讯,数据处理与SD存储卡的数据存储,串口模块、USB模块、有线网络模块和无线Wi-Fi模块4种通讯模块管理与数据通讯,以及各种模块的电源管理控制。以上功能对ARM核心电路的CPU要求较高,不仅要求芯片速度快、可管理多种接口通讯,而且必须具备低功耗特性。因此,本设计中选用三星公司相对成熟的型号为S5P4418的ARM芯片做为核心电路主CPU,且为本仪器定制1款专用ARM核心电路板。
2.2.1 S5P4418型ARM芯片简介
ARM芯片S5P4418采用4核处理器Cortex-A9,主频为1.4GHz×4,片内配置32KB×4 I/D一级缓存,1MB二级缓存,单通道32位数据总线,内存DDR3 800MHz存储带宽,板载4个USB Host 2.0主设备接口和1个USB Otg 2.0从设备接口,板上2路SDIO 2.0,引出1个SD/MMC卡槽,支持4个UART串行口,波特率115200bps,用于GPS通信、普通串口、调试信息输出等。采用此ARM芯片的电路板可直接用作高性能强震记录器的核心CPU板,从而节约强震仪器的开发时间。
2.2.2 高速ARM核心电路板设计
高速ARM核心电路板有数据存储、数据与命令显示、数据传输与命令通讯和电源管理4大功能模块,模块与ARM芯片S5P4418的接口设计见图 5。对数据存储模块设计,首先S5P4418通过GPIO端口从高精度数据采集板卡读出模数转换数据和发送命令信息,同时也通过SPI端口读取GPS信息及秒脉冲数据,然后S5P4418把读出的数据信息和GPS时间信息通过SDIO端口存储到SD卡中,完成永久固态存储强震记录;在数据与命令显示模块的设计中,S5P4418通过1个I2C总线和LCD-TFT控制器管理控制触摸屏显示,实现发送显示数据和命令数据读取的交互操作;本设计中数据传输与命令通讯模块共支持4种通讯端口,分别为串口、USB口、有线网络端口和无线Wi-Fi口。其中串口、USB口、有线网络由S5P4418芯片内置相应功能端口直接管理,而无线Wi-Fi口则由S5P4418芯片的SPI和GPIO口通过无线Wi-Fi模块GS2011扩展实现,由于GS2011内部具有CPU模块和存储器,当需要处理数据时,GS2011与S5P4418通过中断进行程序执行,实现高速数据通讯;在电源管理模块设计中,S5P4418通过线性电源模块输出高精度电源,并且由ARM芯片判断数据存储、显示、通讯等各个模块是否需要工作,来开关各个模块的线性电源供电操作,从而实现整个传感器电路的低功耗管理。
3. 传感器测试
数字化力平衡加速度传感器整机达到了如下技术指标:
(1)加速度测量频率范围为0—200Hz;
(2)加速度测量范围±4g;
(3)输出灵敏度为2.5V/g;
(4)动态范围≥130dB;
(5)线性度优于1%;
(6)2个水平向、1个垂直向,三分向一体;
(7)具有防雷保护电路;
(8)采集器采样率支持50Hz、100Hz、200Hz、400Hz、800Hz、1000Hz;
(9)采集器量程支持±1.25V、±2.5V、±5V、±10V;
(10)支持USB口、网络端口通讯;
(11)单+12V DC供电,整机功耗1.1W;
为了验证该数字化加速度传感器各项技术指标及整机的可靠性,对传感器样机分别进行了传感器的标准振动台检定试验、传感器运行观测实验以及结构模型的振动台观测实验,以上测试实验中采集器采样率均为100Hz,5V量程。
首先,在中航工业北京长城测试计量研究院的标准振动台对传感器的幅频特性、灵敏度等进行测试,表 1给出了该传感器的幅频特性检定结果,从中可以看出传感器频带可以达到200Hz。
表 1 加速度传感器幅频特性标定数据Table 1. Calibration data for amplitude and frequency characteristics of accelerometer频率/Hz 输出灵敏度/V·g-1 1 2.549 5 2.548 10 2.541 20 2.541 40 2.432 60 2.348 80 2.369 100 2.348 110 2.329 120 2.366 140 2.382 160 2.330 180 2.124 200 1.748 220 1.337 其次,在北京白家疃地震台、安徽金寨县地震台进行了多天传感器运行观测实验,取得了大量的强震动噪声观测数据,并且在安徽省金寨县地震台获取了2017年8月8日四川省九寨沟强震记录。随机选取了2017年2月8日3时的北京白家疃地震台噪声观测数据、2017年8月3日3时的安徽省金寨县地震台噪声观测数据进行分析。图 6、7分别为北京白家疃地震台、安徽金寨县地震台的1000点噪声观测数据曲线,图 8、9分别为以上2组数据剔除基线值的噪声波形曲线。从噪声曲线波形来看,该传感器噪声峰值小于5μg,对上述2时段的噪声数据计算得到噪声均方根值和动态范围值(表 2),从计算结果可以看出传感器噪声均方根值最低为0.7235μg,动态范围最大可达134.85dB。为了分析该传感器的频域特性,对2个台站噪声观测数据进行自功率谱分析,图 10(a)、(b)、(c)分别为白家疃台站传感器东西向、南北向和垂直向的自功率谱曲线,图 10(d)、(e)、(f)分别为金寨县台站传感器东西向、南北向和垂直向自功率谱曲线。从图中可以看出,该传感器0.1Hz处功率谱约为-120dB,1Hz处功率谱约为-125dB。同时,该传感器在安徽省金寨县地震台进行观测实验中,获得了2018年8月8日四川省九寨沟强震动记录,其强震动波形见图 11。
表 2 传感器噪声及动态范围Table 2. Noise result and dynamic range of the force balanced accelerometer计算项 北京市白家疃地震台 安徽省金寨县地震台 东西向 南北向 垂直向 东西向 南北向 垂直向 噪声值RMS/μg 1.2051 1.1260 0.7235 1.3720 1.5805 1.2861 动态范围/dB 130.42 131.01 134.85 129.29 128.01 129.86 最后,为了进一步测试该传感器性能,将传感器与美国凯尼公司的TSA-SMA型强震动记录器同时安装在曲面桥结构模型实验的振动台台面上,测试实验在中国地震局工程力学研究所的燕郊园区振动台进行,共安排了50组不同参数的振动波形实验,2套设备都完整记录了50组震动波形,此处仅展示其中1组振动波形曲线,该振动波形发生于2017年12月14日16时29分,图 12为本文研制设备所记录的地震波曲线,图 13为凯尼公司设备所记录的地震波曲线,从上至下依次为东西向、南北向和垂直向震动波形,其中振动台东西向输出波形加速度峰值为0.7g,从图中可以看出2种设备记录的强震动峰值基本相同。
通过以上3种不同类型的测试实验证明,该数字化加速度传感器技术指标达到预期,且质量稳定、性能可靠,该类型仪器完全可以应用于各种低频振动的测量。
4. 结论
本文选用数字加工工厂技术设计了传感器三分量一体化机械底座,并通过改进机械参数、调制电路参数,实现扩展传感器频带范围的目标。在此技术基础上,通过内嵌高精度三通道数据采集器,完成3路同步24位的模数转换,使该传感器成为三分量加速度数字化测量一体机。测试实验数据分析表明,该传感器不仅达到设计的技术指标,并且能够稳定可靠地测量微弱远场强震动,可以普遍应用于强震动观测和各种结构振动监测。
-
表 1 地球化学离子、气体成分等其它观测项目一览表(截至2017年底)
Table 1. List of ion, gas and other items of groundwater observation in seismological geochemistry
区域 观测数量 台站 测项 离子 气体 电导率 其它(PH、硬度) 安徽 1 7 6 1 北京 4 14 4 10 地质所 2 2 2 福建 3 4 1 3 甘肃 3 16 9 7 广东 2 12 4 6 1 1(pH) 河北 15 18 18 河南 4 18 10 6 1 1(硬度) 吉林 1 4 4 辽宁 1 9 4 5 宁夏 1 7 3 3 1 山东 5 22 8 11 2 1(硬度) 山西 1 15 5 10 陕西 2 4 4 四川 8 53 33 9 7 4(pH) 天津 9 15 8 7 新疆 16 36 9 21 4 云南 4 8 1 4 1 2(pH) 合计 82 262 103 129 21 9 表 2 地震地球化学离子、气体成分等其它测项观测仪器一览表
Table 2. Instrument of ion, gas and other items of groundwater observation in seismological geochemistry observation
观测内容 仪器名称 型号 采样率 观测方式 启用时间 离子、硬度 水质观测仪 PXJ-1C 日 人工 1980 电位滴定仪 ZD-2 日 人工 1995 离子 离子色谱仪 CIC-200 日 人工 2014 滴定管 日 人工 1980 电导率 电导率仪 CM-30/DDS-307/DDS-11A 日 人工 2008/2009/2010 pH值 酸度计 PHS-2/PHS-3C 日 人工 1990/2002 气体量、氮气浓度、甲烷浓度、氩气浓度、氧气浓度、二氧化碳浓度、氢气浓度、氦气浓度 气相色谱仪 SP-3400/TP-2060T/SP-2304A 日 人工 1994/1994/1981 二氧化碳 二氧化碳测量管 EY-2 日 人工 1990 氢气 痕量氢在线自动分析仪 ATG-6118H 小时 数字 2014 氦气 测氦仪 WGK-1 小时 数字 1994 表 3 地球化学离子、气体成分等其它测项异常统计表(1966—2012)
Table 3. Statistical anomaly for ion, gas and other items of groundwater observation in seismological geochemistry
测项目 异常次数 气体观测 气体总量 硫化氢 二氧化碳 氢气 氮气 氩气 氧气 甲烷 3 4 15 4 2 3 2 2 离子观测 氟离子 钙离子 碳酸氢根离子 氯离子 镁离子 硫酸根离子 碳酸根离子 钾离子 37 15 16 12 11 4 2 1 其它观测 电导率 pH值 5 7 -
车时, 蒋海昆, 付虹等, 2014.中国震例(2003-2006).北京:地震出版社. 车用太, 王吉易, 李一兵等, 2004.首都圈地下流体监测与地震预测.北京:气象出版社. 杜方, 蒋海昆, 杨马陵等, 2018.中国震例(2008年5月12日四川汶川8.0级地震).北京:地震出版社. 高文学, 蒋凤亮, 高庆华等, 2000.地球化学异常(地震预测整体观的探索).北京:地震出版社. 高小其, 迪里夏提, 许秋龙等, 2002.乌鲁木齐9、10号泉F-映震特征的研究.地震学刊, 22(1):5-11. doi: 10.3969/j.issn.1672-2132.2002.01.002 蒋海昆, 杨马陵, 付虹, 2018a.中国震例(2011-2012).北京:地震出版社. 蒋海昆, 付虹, 杨马陵, 2018b.中国震例(2007-2010).北京:地震出版社. 林国元, 朱继承, 黄其昌等, 2015. 2013年仙游MS 4.8地震福建水化学前兆异常.地震地磁观测与研究, 36(2):64-70. doi: 10.3969/j.issn.1003-3246.2015.02.012 吕超甲, 刘雷, 周晓成等, 2017.甘肃东南地区温泉流体地球化学特征.地震, 37(1):52-60. http://d.old.wanfangdata.com.cn/Periodical/diz201701006 宋臣田, 宋彦云, 唐久安等, 2008.地震监测仪器大全.北京:地震出版社. 张炜, 王吉易, 鄂秀满等, 1988.水文地球化学预报地震的原理与方法.北京:教育科学出版社. 张肇诚, 罗兰格, 李海华等, 1988a.中国震例(1976-1980).北京:地震出版社. 张肇诚, 罗兰格, 李海华等, 1988b.中国震例(1966-1975).北京:地震出版社. 张肇诚, 罗兰格, 李海华等, 1990.中国震例(1981-1985).北京:地震出版社. 中国地震局, 2011. DB/T 3-2011地震测项分类与代码.北京:中国地震局. 中国地震局监测预报司, 2007.地震地下流体理论基础与观测技术.北京:地震出版社. 周晓成, 孙凤霞, 陈志等, 2017a.汶川8.0地震破裂带CO2、CH4、Rn和Hg脱气强度.岩石学报, 33(1):291-303. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ysxb98201701023 周晓成, 石宏宇, 陈超等, 2017b.汶川MS 8.0地震破裂带土壤气中H2浓度时空变化.地球科学进展, 32(8):818—827. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqkxjz201708004 -