地电观测干扰数据检查仪研制.pdf

郅红魁，叶青，孙召华，等.地电观测干扰数据检查仪研制 J.华北地震科学，2023,41（3）：52-59.doi:10.3969/j.issn.10031375.2023.03.008.ZHI Hongkui，YE Qing，SUN Shaohua，et al.Development and Application of Interference Data Anomaly Instrument for GeoelectricObservationJ.North China Earthquake Sciences，2023,41（3）：52-59.doi:10.3969/j.issn.10031375.2023.03.008.地电观测干扰数据检查仪研制郅红魁1，叶 青2*，孙召华1，赵彦旭1，谢佳兴1，郜冰莹1，郭少峰1（1.河南省地震局,郑州 450016；2.中国地震台网中心，北京 100045）摘要：地电观测过程中，由于受漏电和观测系统故障会引起干扰数据，如突跳、台阶、超差等，这些干扰数据不易被快速判断，需要做大量检查工作才能完成。为解决此问题，研制了地电观测干扰数据检查仪，能实时记录仪器工作状态和地电场变化，当出现干扰时调取相应干扰时段的秒数据或分数据，通过对干扰数值波形特征分析，进而快速判别并排除干扰。实践应用表明，该仪器对干扰异常能够较好识别，提高了排查效率。关键词：地电；故障干扰；供电波形；数据分析；干扰源排查中图分类号：P315.72 文献标志码：A 文章编号：10031375（2023）03005208doi:10.3969/j.issn.10031375.2023.03.008 0 引言地震地电阻率观测在国内地震监测、预报中发挥着重要的作用，由于城市化进程的加快，地表电阻率台站观测环境干扰严重1。各种漏电干扰不断增加，同时很多台站使用的“十五”“背景场”项目的仪器日益老化，容易引起干扰数据。造成原因主要有 2 种：一是观测系统故障，主要由外线路老化、避雷装置故障、接头接触不良等造成；二是场地环境变化，主要由观测场地附近新增或拆除供电设备、金属管道和交直流漏电造成数据异常等。对于上述干扰的检查方法，主要是通过仪器标定、替换、外环境巡查、线路绝缘等。但上述检查方法对测区内漏电无法准确判定，存在排查难度大、耗时长，导致无法快速、准确地判断干扰原因。目前，台站没有地电专用仪器故障、漏电检查仪，不定时出现的异常判别对工作人员地电观测技术要求较高。本文针对地电仪器的测量原理和测区漏电特性，研发了地电观测干扰数据检查仪，为地电台站提供一个新的排查依据，从而提高工作效率。1 仪器设计思路、工作原理1.1设计思路地电观测仪器观测时间间隔较长，中间无法记录到地电场变化和仪器工作状态，出现数据异常时不能准确判断干扰原因，因此需要研制一台高精度、高分辨率、较高采样率的的仪器，接入地电测量极获得实时数据，再对这些干扰数据变化特征进行分析，确定干扰原因。1.2结构、工作原理该仪器由 ARM 工控板、PIC18F4520 8 位单片机、U 北斗/GPS 授时模块+天线、ADS12568 通道模数转换模块、锂电池/12 V 线性电源、数据线、测量线、4 根铜针（野外漏电检查）等组成。仪器采用32 位的 ARM9 工控板，内含 WINCE6.0 操作系统，系统内有 c#开发的地电观测环境干扰数据采集波形显示软件。该软件通过 RS232 串口向 PIC 单片机发送读取数据指令，单片机处理通过 SPI 接口与24 位高精度 ADS156 模数采集模块通信，发送读取测量电压数据命令；采集模块对地电测量极 M、N 电极或观测场地电压测量1，将测量数据通过SPI 接口传给单片机；单片机将数据进行计算处理、收稿日期：2022-12-21基金项目：北京市自然科学基金（8212045）第一作者简介：郅红魁（1974），男，河南洛阳人，高级工程师，主要从事地震监测预报、仪器运维等研究.Email：*通讯作者简介：叶青（1977），女，河南南阳人，高级工程师，主要从事地震电磁学方面研究.E-mail：qing_第 41 卷第 3 期华 北 地 震 科 学Vol.41，No.32023 年 7 月NORTH CHINA EARTHQUAKE SCIENCESJul.，2023格式转换后，发送给工控板上的采集软件；采集软件再按指定格式生成秒数据和分钟值文件。当地电数据出现异常时，用 FTP 软件通过本地或远程下载数据进行分析，识别出观测干扰原因。2 仪器的研制2.1硬件的实现仪器由工控板、电源、GPS 授时、数据处理芯片、AD 模数转换模块、数据线及 4 根铜针组成。工控板采用 LJD-eWinV5-LT3 开发板，V5-3.5 寸触控屏 wince6.0 中文系统、3 路 RS232/RS485 通信串口、TF 卡最大支持 64 G 存储卡、百兆网口，主要功能是完成数据存储、实时波形显示、网络通信、数据传输等。电源部分由 12 V、4 A 的 DC 外接电源、内部 LT3045 单电源低噪声线性模块和 12 V 锂电池组成，主要作用是给仪器供电，同时给锂电池充电，保证在野外查漏电干扰时锂电池可以为仪器工作8 小时以上供电。单电源线性模块是给数采 AD 模数转换模块和授时模块提供电压。GPS 授时由高精度 GPS/北斗一体化模块 NEO-M8M 陶瓷天线组成，主要功能是给仪器授时，保证仪器时间与工作仪器时间同步。数据通信处理芯片是由8 位PCI18F4520单片机组成，主要功能是通过 SPI 接口读取 AD 模块数据进行计算处理和格式转换成字符串形式，通过单片机串口与工控板 RS232 串口 2 通信，将数据发送到工控板的数据采集软件中。AD 模数转换是由 24 位高精度的 ADS1256 模块组成，内含有 8 路输入通道，采用差分输入，分辨率达 1 V，低噪声低温漂 Vishay 晶圆金属膜电阻，测量电压时确保不受温度影响，数字信号高频衰减电阻，隔离外部高噪 数据下载分析观测系统是否正常故障位置漏电位置锂电池电源滤波液晶显示屏ARM 工控板北斗/GPS授时模块中央处理器RS232 端口 1RS232 端口 2格式转换文件存储8 位单片机数据通信处理、SPI 接口ADS1256 电压采集模块地点观测系统各测量极信号检测地电场观测系统各测量极信号检测地电和地电场观测场地漏电检查网口FTP协议观测场地是否漏电通道 1、2 图 1 仪器设计结构图Fig.1 Instrument design structure diagram 图 2 工控板Fig.2 Industrial control boards 图 3 授时模块Fig.3 Timing modules3 期郅红魁，等：地电观测干扰数据检查仪研制53 声。数据线采用的是 4 芯信号屏蔽线，一端接入仪器，另一端为 4 个鳄鱼夹子，分别代表示 A1（北+）、B1（南-）、A2（东+）、B2（西-）；当地电数据出现干扰时，将 4 个夹子，A1 接北南向测量极 M1、A2 接测量 极 N1，A2 接 东 西 向 测 量 极 M2、B2 接 测 量 极N2。如果是室外查漏电干扰源时，将 4 根铜针正十字交叉间距 5 m 插在测区干扰附近，记录漏电干扰源数据，根据干扰方向和幅度确定干扰源位置。2.2软件系统的实现软件系统由数据采集波形实时显示软件、授时软件、PIC 单片机数据处理等组成。波形实时显示软件是基于 Visual Studio 2008 运行环境来开发的，使用 c#语言进行编程，可根据自己的需求进行编辑，移植编译运行在基于 wince 系统的 arm 体系结构操作平台上。软件系统的主要功能是显示记录波形和数据备份。具体实现方法：首先对两个串口初 始 化，通 过 串 口 控 件 comport 中 BaudRate 和PortName 属性设置波特率和 COM 端口号，用文件读取函数 StreamReader 获得所有参数字符串，用Split()分割字符串获得台站每个参数，File.Exists查看当天有无数据文件，无则通过 File.AppendText函数新建空文件，将所有缺测数据置为 NULL。利用 timer1 定时器每秒钟触发一次，利用 comportbf.Write(zhkwhx)函数向单片机发送电压采集指令，并通过 comportbf.ReadExisting()读取单片机发送的数据，用 ByteArrayToHexString(byte data)函数将十六进制数据转换成字符串，然后每分钟对秒数据进行平均计算得到分钟值数据，同时按照 TXT 文本格式“年月日时分（秒）+制表符+数据 1+制表符+数据 2+换行符”，方便直接复制粘贴到 EXECL 表格绘图形，通过 File.AppendText 函数实时保存，并按时间序列存入秒、分钟数据日文件，文件名分别为YYYYMMDD.txt 分和 YYYYMMDD 分钟.txt。仪器内设有两个选项，“地电自然电位差”和“地电场”，根据干扰异常数据选择相应测量方式，当选择“地电场”，软件通过地电场公式换算自动计算出电场值，按上述方法存入 DC 文件夹中。软件实时绘图功能，利用 Graphics.FromImage(bMap)函数创建一个画板，在画板上用 gph.DrawLine 函数用循环语句绘制 2 分钟两个测向的秒数据曲线。3 数据波形分析方法s连接该检查仪通过 LeapFTP 软件登陆下载文件，将数据复制粘贴到 EXECL 表中绘制曲线图形，通过地电仪器在工作前、中、后的数据变化形态特征，分析出地电干扰数据异常原因和干扰位置。首先，需要了解地电阻率供电原理、测量时序，有助于快速数据异常分析。地电阻率采用四极对称装置进行测量，工作过程中采用 A、B 电极正反向供电，利用测量极在供电期间测量得人工供电电位差V，通过视电阻率计算公式=K*V/I 得到地电阻率（图 6）。从式中可以看出影响地电阻率观测值主要是V 和 I，而V 则是由 N 次正反向供电中测得，如果 B、C 处任何一处出现方波畸变，都会影响地电阻率观测值。AVAMNB 图 6 地电阻率四极对称装置图Fig.6 Four-pole symmetrical device diagram of geoelectricalresistivity 地电阻率仪器观测的供电方法和测量时序理论上为 M、N 之间人工电位差测值变化时序2，通过了解观测供电波形特征，有助于分析仪器工作状态 图 4 AD 转换模块Fig.4 AD conversion module 图 5 单片机处理器Fig.5 Microcontroller processor54华 北 地 震 科 学41 卷 和观测环境变化情况。仪器测量时序如图 7，5 分时开始先测 M、N 之间自然电位差 A，之后进行 2 次正反向供电测量供电电流 B，接着进行 N 次正反向供电测得人工供电电位差 C，最后测量一次正反向供电测量供电电流 B；然后开始计算 N 次的地电阻率值，每小时地电阻率观测值可视为由 N 次观测的地电阻率的平均值。AI/At/sTB BBC C C C CCB BBC C C C CCN+1 图 7 地电阻率仪器供电测量时序图Fig.7 Timing diagram of geoelectrical resistivity for instrumentmeasurement 以洛阳台地电阻率供电波形为例（图 8），第一条曲线 A1 为观测前 NS 向自然电位差变化形态，B1 段为两次正反供电，用于测量电流大小，C1 段为4 次正反向供电测量人工供电电位差，为了确保电流稳定可靠 C1 段后 B1 段再进行一次正反向电流测量；然后仪器自动计算自然电位差、地电阻率、均方差，完成一次 NS 向地电阻率观测，D1 段为EW 向在观测过程中 NS 向记录到的供电波形。第二条曲线 EW 向和第一条 SN 向供电波形是一样的，但位置左右相反。从图