第２２卷第４期 ２０１０年ｌ２月 高原地震 ＰＬＡＴＥＡｕ ＥＡＲＴＨＱＵＡＫＥ ＲＥＳＥＡＲＣＨ Ｖ０１．２２ Ｎｏ．４ Ｄｅｃ．２０１０ 地电阻率资料日常处理软件研制 狄国荣 ，汤 吉。，高曙德 ，苏永刚 ，叶媛媛 ，陈彦平 （１．甘肃省地震局陇南中心地震台，甘肃陇南７４６０００； ２．中国地震局地质研究所，北京１０００２９） 摘要：地电阻率观测在地震前兆观测中具有重要的地位，大量的观测资料表明，多数中强地震前地 电阻率出现了异常变化，而快速地发现、甄别地电阻率的变化是地震前兆信息还是干扰，高效率地做好 地电阻率资料日常处理工作是前提。由此，依据地电阻率理论和有关规定，基于ＶｉｓｕａｌＣ＋４－６．０平台， 开发了一款关于ＺＤＳＢ或ＺＤ８ＢＩ地电仪产出的地电阻率资料日常处理软件，旨在介绍本软件算法的理 论原理、设计思路及实现方法，并结合陇南中心地震台地电阻率观测资料给出了相应的应用实例。应用 表明，研发的该软件可以将地电阻率资料日常处理工作系统化，将日常处Ａ．Ｚ－作中大量繁杂的手工操作 自动化，能及时地对资料的异常变化做出初步判定，快速地为地震预测部门提供依据。 关键词：地电阻率；前兆异常；地电仪；软件 中图分类号：Ｐ３１９．３ 文献标识码：Ａ ．文章编号：１００５—５８６Ｘ（２０１０）Ｏ４一ｏｏ５２—１Ｏ ０引 言 地电阻率观测手段在我国已经进行了４０多 年的应用，积累了相当丰富的第一手观测资料，大 量的观测资料表明，多数中强地震前地电阻率发 生了异常变化。例如：在１９７６年８月１６日四川 日常处理方法的系统化、自动化程度不高，手工计 算量大，过程繁琐。为了高效率地做好原始资料 的日常处理，快速对地电阻率资料异常变化作出 初步判定，并及时为进一步深入研究提供日常处 理的规范化基础数据，作者基于ＶｉｓｕａｌＣ＋＋６．０ 平台，依据中国地震局颁发的《数字地震前兆观 松潘—平武 ｓ７．２（Ａ：１０５ｋｍ）地震前，甘肃武都 地震台地电阻率北７３。西测道的观测资料出现明 ｉ贝０技术规范》（电磁分册） 和地电台管组的相关 规定开发了本软件。 显的阶段性负异常，异常幅度达一１４％，该异常特 征被作为讨论地电阻率预报方法的典型震例 Ｊ。 唐山地震前，在震中周围１００～２００ｋｉｎ范围的多 个地电阻率观测台站，也出现了明显的异常特 征 。还有很多地电阻率观测资料都在一些地 震前有明显的异常反应 －４］。２００８年５月１２ Ｅｔ １ 地电阻率观澳４的基本原理与观测 数据 图１中，Ａ、Ｂ为２个供电电极，Ｍ、Ｎ是２个 测量电极。用测量得到的电流，和电位差△ 计 算大地电阻率： ＾Ｔ， 汶川８．０级地震前，距震中３２０ｋｍ的甘肃陇南武 都汉王台观测到了２—３年尺度的长趋势地电阻 率上升和下降异常 Ｊ。 随着地震预测研究的不断发展，对地电阻率 Ｐ ：Ｋ 』 （１） 其中，Ｐ 为视电阻率，只在地下为一种 均匀介质时等于其真电阻率。通常地下为不均 匀介质，由（１）式计算得到的大地电阻率一般不 等于某种介质的真电阻率，而是地下各种介质导 资料分析和研究也逐渐深入，对原始资料的日常 处理工作量急剧增加，但是各台站地电阻率资料 收稿日期：２０１０—０８—１２ 基金项目：科技支撑计划项目（２００８ＢＡＣ３５Ｂ０１一１．２００８ＢＡＣ３５Ｂ０１—８）。 作者简介：狄国荣（１９８１～），男，甘肃白银人，工学学士，助理工程师，主要从事地震监测与资料分析处理工作。Ｅ—ｍａｉｌ：ｘｕ￣－ ｈｕ９１０７＠ｓｉｎａ．ｃｏｒｎ。 第４期 狄国荣，等：地电阻率资料日常处理软件研制 ５３ 电性的综合反映，即视电阻率，在地震地电学中通 ｒ 是图１中Ａ、Ｍ间的距离，其余类推。目前， 常视电阻率即指地电阻率。Ｋ是一个常数，称为装 我国开展的地震前兆监测的地电阻率测量，绝大 置系数，其大小取决于４个电极的排列方式，由 多数采用的是四极对称装置，即４个电极排列在１ （２）式确定： 条直线上，且相对中心点是对称排列的，观测仪器 （２） 大多采用国产的ＺＤ８Ｂ或ＺＤ８ＢＩ地电仪，该类仪 器采样率１次／小时，２４次／天。 ＦＡＭ ｒ＾Ｎ ＴＢＭ ＴＢＮ 图１ 视电阻率测量示意图 ＺＤ８Ｂ或ＺＤ８ＢＩ地电仪每天将采集的数据按 ２．２制作月报表 照《数字地震前兆观测技术规范》（电磁分册）规 每月观测结束时，将每天 日均值，Ｐ 日均 定处理后，用十六进制补码的形式，以原始数据格 值及其均方根误差，小时值均方根误差Ｏ－ 一。及日 式“ ．ｏｒｇ”保存，其中包括自然电位差 的小时 均值提取出来，再分别计算出这４项的５日均值、 值、日均值；视电阻率Ｐ 的小时值及其均方根误 月均值和Ｐ 月均值的均方根误差，导入地电台管 差 一，和相对均方根误差 ；ｐ 的初始日均值 组设计的Ｅｘｃｅｌ表中上报台管组。 Ｐ 初及其均方根误差（ ）初和相对均方根误差 。、Ｐ 日均值、５日均值和月均值及小时值均 （ ）初；将小时值中Ｉ Ｐ —Ｐ 初ｆ＞２（ ）初的测 方根误差的日均值计算用（３）式： 值剔除后，重新计算所得的Ｐ 平均值ｐ 及其均方 根误差（ 一 ）日和相对均方根误差（ ）目；．Ｄ 小 一 ∑置 Ｘ＝生 （３） 时值均方根误差 一 的Ｆｔ均值（ 川）日；被剔除 ｎ 的测值个数即坏点数。其中， 。精度保留～位小 其中， 表示被计算项，ｎ表示计算项的数目。 数，视电阻率各种均值及其均方根误差和相对均 日均值用小时值数据计算，５日均值和月均 方根误差精度都保留２位小数。 值分别用日均值数据计算。 ２ 地电阻率资料日常处理 Ｐ 日均值、月均值的均方根误差（ｏｒ 一 ） 、 （ 一。）月计算用（４）式： ． 依据《数字地震前兆观测技术规范》和地电 台管组规定，对“ ．ｏｒｇ”数据文件资料现阶段所 ｎ一１ 作的日常处理包括以下几项。 （４） ２．１ 数据查看和表示 其中，ｎ为观测次数 Ｐ为第ｉ次观测值，Ｐ 为 将每天十六进制补码形式的数据资料转换成 ，ｚ个Ｐ 算术平均值。ｏｒ 币 最低有效位一致。 十进制形式查看并保存。每月观测结束时，将本月 ２．３ 外线路漏电影响系数８ 、￡ 和外负载电阻 各测道分量每天的小时值以图形的方式表示出来 计算 并保存、上报地电台管组。每年观测结束时，将本 每月定期巡检外线路，并计算外线路漏电影 年度 日均值、５日均值和月均值以图形方式 响系数 、 ：和外负载电阻 ，对于外线路漏电影 表示出来并上报台管组。 响系数和外负载电阻的计算分别为： ５４ 高 原 １＝ ×１００％ （５） ｚ ×１００％ （６） Ｒ＝ （７） 』ＡＢ 按照科学测量数值计算的习惯方法：对于ｓ． 来说，当 ≥０．０１ｍＡ时， 用（５）式计算，取位与 的有效位数保持一致；当ＩＬ＝０．００ｍＡ时，根据 误差理论中关于有效数字的规定， 的最大误差 ＩＬ＝ ＜０．０１ｍＡ，因此 ｓｌ＜（半）×１００％ （８） 最后的结果，不同０的表示形式代表了不 同的意义。如果 １＝０％，则表示ｓＩ＜１％；ｓ１＝ ０．０％，则表示ｓ１＜０．１％， １＝０．００％，则表示ｓｌ ＜０．０１％，依此类推。 计算利用（６）式，取位方 法与 相同。 ２．４ 月数据连续率 月、完整率 月 月数据连续率和完整率的计算分别由下式计 算： Ｌ月： ＩＪ月一 ｇｔ ×１００％＾ （９）＼７， ／Ｊ应测个数 Ｗｎ： 二鲁 ／］应测个数 ×１００％ （１Ｏ） 其中，“月应测个数”＝２４×ｒ／，，ｎ为当月天数， “月缺测个数”为当月无测量和粗差总数，“月舍弃个数” 是当月在计算该测道Ｐ 日均值时，将Ｉ Ｐ Ｐ 初Ｉ＞ ２（ｏ－．一 ）初的测点值舍去的个数。 ２．５ 年度观测资料精度系数厂值 每年观测结束后，计算出本年度观测资料精 度系数 值，该值是现阶段衡量地电祖率资料内 在质量的一个重要指标，定义如下： ，－１＿（ １ ｋｏｑ＋ ）（１１） 其中： ． ２ １２ ． 后 。＝了１∑∑［（ｏ－ ）月 ／（ ）月 ］ （１２） ：１ ｉ＝１ 。 １ ２ ｘ＋４ ＝了１∑∑［（ ）月 ／（ｐ ）月 ］ （１３） 一 １ ｌ （ ）月 为第 道第ｉ月的ｏｒ 的月均值， 地 震 第２２卷 （Ｏ－ 一 ）月 为第．『道第 月由（ ）Ｈ算出的月均 值；（Ｐ ）月 为第 道第ｉ月以Ｐ 日均值计算的月均 值， 为起始月份。计算精度要求达到０．０００１，并 限定：当（ ｌ／１２）＜０．００２和（｜ｉ｝以／５）＜０．００１时， 均按０对待。 ２．６ 观测资料异常的初步判定 视电阻率正常动态指标值是判定是否存在异 常的指标，它定义为： 均值被动幅度一 。目 ＝ ＿｜ｎ — ×１００％ （１４）＼ ， ∑Ｐ 』三 —一 ／／， 其中，／７，为天数，Ｐ 为第ｉ天的Ｅｔ均值，Ｐ… 为 这些天Ｐ 中的最大者，Ｐ 为这些天Ｐ，中的最小 者，如果波动幅度Ａ 均值波动幅度＜１．５％，并且在这 一时间段内无地震发生，那么可以将这一电阻率 的平均值视为正常动态指标值；如果当电阻率均 值偏离正常动态值的相对变化幅度大于１．５％， 且超过正常动态值时段标准偏差３倍时，即判定 为异常变化。异常的起始时间应向前追溯到异常 开始偏正常值的时间算起。在异常没有恢复到正 常值的情况下，从异常起始时间到所分析的数据 截止之日的时段为异常持续时间‘；当异常已经恢 复到正常值的情况时，从异常起始时间到恢复正 常动态值的时段为异常持续时问；在异常持续时 间内，地电阻率偏离正常值的最大相对变化值为 异常幅度。 ３ 软件设计与实现 实用软件的设计，应考虑如下３个因素：①软 件要便于使用，要有方便灵活的数据接口；②有 坚实的理论作为基础，并结合实际观测数据处理 的需求，结果输出要便利化、可视化；③软件要有 不断发展的生命力。据此，本软件设计了２大模 块。 （１）地电阻率资料日常处理模块。主要是将 日常处理工作避开现阶段的繁琐处理过程，将其 系统化、自动化。其基本功能包括数据查看、制作 月报表、漏电检查计算、月数据连续率和完整率计 算、年度观测资料精度系数．厂值计算、每月各测道 小时值绘图、全年观测值绘图（包括 。和Ｐ 的年 度日均值图、５日均值图、月均值图， 年度日均 第４期 狄国荣，等：地电阻率资料日常处理软件研制 ５５ 值图）等。 ３．１ 软件的基本结构 （２）功能模块。主要是增强软件的生命力，可 根据地电阻率的观测规范和行业技术要求， 以根据用户的需求添加功能。现阶段其主要功能 考虑到台站实际应用和日常数据管理以及数据上 包括资料提取与编辑，正常动态指标计算等。 报的需要，软件总体结构设计如图２所示。 图２软件总体结构图 ３．２ 软件数据接口 ＣｈｅｃｋＢｏｘ控件。 为了提高软件数据输入的准确性与效率，软 对资料处理后产生的数据及图像，软件提供 件的各个功能直接使用“ ．ｏｒｇ”原始数据文件， 了即时存储功能。处理数据以文本形式保存，便于 避免了手工录入误差。在设计各个功能模块与用 对数据结果的查看与后续处理。图像以“ ．ｂｍｐ” 户交互的数据接口时利用了ＭＦＣ提供的文件管 格式存储。 理对话框ＣＦｉｌｅＤｉａｌｏｇ类，很直观、方便地让用户 ３．３各功能实现 直接选择所要处理的数据文件，极大地简化了操 ３．３．１数据查看功能 作。 首先，利用ＣＦｉｌｅ类的Ｒｅａｄ方法将“￥．ｏｒｇ” 由于中国地震局地电前兆监测网的“九五” 数据文件十六进制补码形式的数据读出，然后，利 和“十五”观测数据的文件命名方式不一样，软件 用码制转换规则将其转换成十进制形式并利用 在处理连续一段时间内的资料时，具有自动搜索 ＥｄｉｔＢｏｘ控件将其可视化。 下一日文件并继续处理的功能，所以，为方便软件 ３．３．２制作月报表功能 按照相对应的规则来自动搜索处理文件，作者在 利用ＣＦｉｌｅ类的Ｒｅａｄ方法将每月每天的 设计各个功能时给出了选择文件类型的 “ ．ｏｒｇ”数据文件里的 日均值，Ｐ 日及其均方 ５６ 高 原 地 震 第２２卷 根误差，还有小时值均方根误差的日均值 （ ） 提取出来并转换成～卜进制形式，再依据 文件中满足要求的数据提取出来并转换成十进制 形式在ＥｄｉｔＢｏｘ控件中显示，利用ＣＦｉｌｅ类的 Ｗｒｉｔｅ方法实现保存功能。 （３）式分别计算出各自的５ Ｅｌ均值和月均值，依 据（４）式计算出Ｐ 月均值的均方根误差，最后利 用Ｏｆｉｆｃｅ产品的ＯＬＥ Ａｕｔｏｍａｔｉｏｎ自动化程序接口 将提取出的数据和计算后的数据导人到台管组设 计的Ｅｘｃｅｌ表格中生成当月的报表。 ３．３．３ 月数据连续率、完整率检测 首先，计算出当月的天数ｎ，即当月“ ．ｏｒｇ” 数据文件个数，依据ｎ计算出各测道“月直测个数” ＝２４ｎ（采样率１ Ｒ／／｜，Ｕ￣，２４次／天）。其次，利用 ＣＦｉｌｅ类的Ｒｅａｄ方法从当月“ ．ｏｒｇ”数据文件中 分别提取出每天各测道舍点数和缺测数，把各测 道每天的舍点数分别求和，得当月各测道的 “月舍弃个数”；把各测道每天缺测数分别求和，得当 月各测道缺测数总和。如果当月有粗差，将各测 道粗差总数利用ＥｄｉｔＢｏｘ控件由用户输入，然后 由缺测数总和以及粗差数分别计算出当月各测道 “月缺测个数”。最后，利用（９）、（１０）式分别计算出 当月各测道数据连续率、完整率，并用ＥｄｉｔＢｏｘ控 件将这些数据可视化。 ３．３．４观测资料精度系数．厂值计算 用１，２，３分别代表第１测道、第Ⅱ测道、第Ⅲ 测道，１…ｌ２分别表示第一月…第十二月，由Ｅｄｉｔ— Ｂｏｘ控件读入用户的要求，利用ＣＦｉｌｅ类的Ｒｅａｄ 方法从本年度的“ ．ｏｒｇ”数据文件中提取出满足 条件的数据并转换成十进制形式，而后将这些数 据代入（１１）、（１２）、（１３）式计算出观测资料本年 度的精度系数，并用ＥｄｉｔＢｏｘ控件将其可视化。 ３．３．５图形绘制 利用ＣＦｉｌｅ类的Ｒｅａｄ方法将所需的数据从 “ ．ｏｒｇ”数据文件中提取出来并转换成十进制形 式，然后利用ＣＣｌｉｅｎｔＤＣ类的ＭｏｖｅＴｏ和ＬｉｎｅＴｏ方 法，在对话框窗口中绘制图形并显示。为了让图 形具有保持功能，在对话框窗口的ＰｒｅＴｒａｎｓ— ｌａｔｅＭｅｓｓａｇｅ方法中调用ＣＣｌｉｅｎｔＤＣ类的ＭｏｖｅＴｏ 和ＬｉｎｅＴｏ方法ｏ ３．３．６资料提取与编辑 资料提取是依据“ ．ｏｒｇ”数据文件的规定以 及用户的要求，包括时间段，需要的数据类型（整 点值数据，日均值数据）及数据项（ ，Ｐ ， ， ），利用ＣＦｉｌｅ类的Ｒｅａｄ方法将“ ．ｏｒｇ”数据 资料编辑是利用ＣＦｉｌｅ类的Ｒｅａｄ方法将用 户要求的资料读人到ＥｄｉｔＢｏｘ控件中，然后利用 ＣＥｄｉｔＢｏｘ类的字符处理方法编辑这些数据，编辑 完毕后保存。 ３．３．７动态指标值的计算 首先，利用ＣＦｉｌｅ类的Ｒｅａｄ方法将一段时间 “ ．ｏｒｇ”数据文件里的Ｐ 日均值提取出来，根据 冒泡排序法求出这个时间段内ｐ，的最大值、最小 值，然后用（３）和（４）式计算出这段时间内Ｐ 的 均值和该均值的均方根误差，再利用（１４）式计算 出这个时段Ｐ 的波动幅度 均值渡动幅度均在Ｅｄｉｔ— Ｂｏｘ控件中显示以及根据用户需要利用ＣＦｉｌｅ类 的Ｗｒｉｔｅ方法保存。用户可以依据这些数据和当 地在这段时间内的地震活动情况判断Ｐ 这个均 值是否可以充当正常动态指标值。 ４软件应用及对比处理实例 采用陇南中心地震台２００７年地电阻率观测 资料，将传统的资料１３常处理方法和本软件对比 来做如下几项日常处理工作。 ４．１数据查看检验 查看２００７—０１—０１～２００７—０１—３１这一时 间段的数据。点击“资料日常处理”菜单下的“查 看”功能，然后在弹出的页面下选定“开始日期” 和“结束日期”，接着点击“查看”按钮后立即可得 到如图３所示数据（由于ＥｄｉｔＢｏｘ使用滚屏显示， 所以本图只能显示２００７年１月１日００—２３点整 点值数据）。 由此，可以看出： （１）本软件可以快速、简洁、直观的查看任意 时间段的相关数据； （２）可以很方便的通过查看仪器每天的零 输入值和坏点数来简单判断仪器的工作状况以及 资料受影响的程度，为及时发现资料异常提供方 便。 ４．２制作月报表 传统制作方法是将这一个月的“ ．ｏｒｇ”数据 文件转换后，计算出报表所需的每天的各个测向 的日均值观测数据，然后以人工录人的方式制作。 第４期 狄国荣，等：地电阻率资料日常处理软件研制 ５７ 图３ 由本软件转换所得陇南台２００７年１月地电阻率数据 而用本软件制作，只需点击“资料日常处理” 动在后台将报表所需的这个月的所有数据计算 出，直接导人到台管组设计好的Ｅｘｃｅｌ模板中，即 可得所需报表，如图４所示。 菜单下的“制作月报表”功能后，在该功能页面下 选定２００７—０１—０１的“ ．ｏｒｇ”数据文件，软件自 图４陇南台２００７年１月地电阻率月报表 通过对比制作可以看出，将繁杂的数据检索、 和２００７年度厂值。 计算工作转为软件处理，在很大程度上减轻了工作 人员的计算量，避免了手工录入数据误差的出现。 ４．３月数据连续率、完整率及观测资料精度系数 ．传统处理方法计算月数据连续率、完整率，首 先将这个月所有的十六进制“ ．ｏｒｇ”数据文件转 换后，手工检索出公式（１０）、（１１）所需数据，而后 厂值计算 人工计算。计算．厂值，首先从２００７年全年１～１２ 月的报表中手工检索出公式（１３）、（１４）所需的数 例如：计算２００７年２月数据连续率、完整率 ５８ 高 原 地 震 第２２卷 据，然后利用ｆｃａｌｃｕｌａｔｅ软件计算。 下的“观测精度（ 值）”功能后，在弹出的页面下， 用本软件计算连续率、完整率，在点击“资料 按照提示输入相应的数据代号，点击“计算”即可 日常处理”菜单下的“连续率Ⅱ，完整率”功能后， 得如图６的结果。 在弹出的页面下，选定这个月的“开始日期”和 通过对比计算可以看出，本软件在计算月数 “结束日期”后再点击“计算”按钮即可得结果，如 据连续率、完整率以及观测资料精度系数 值时 图５所示。计算厂值，点击“资料日常处理”菜单 较传统计算方法简捷、方便。 图５陇南台２００７年２月数据连续率、完整率 图６陇南台２００７年度／值 ４．４图形输出 “ＭａｐＳＩＳ软件”绘制。对于月均值图，按照《数字 例如：绘制２００７年度视电阻率日均值和月均 地震前兆观测技术规范》（电磁分册）要求，要绘 值图。 制出每月Ｐ 均值的误差限线段（以Ｐ 月均值为中 传统绘制方法：首先，利用软件将全年的资料 心，沿纵坐标方向画一对称的细线段，线段的总尺 转换，如：“ＥＰＣＣ一通讯控制软件”等，而后利用 度为该月Ｐ 月均方根误差值的２倍），而利用 第４期 狄国荣，等：地电阻率资料日常处理软件研制 ５９ ＭａｐＳＩＳ并不能绘出误差限线段，该线段要在绘出 的月均值图上手工测量换算绘制。 图７陇南台２００７年度Ｐ 日均值曲线 图８陇南台２００７年度Ｐ 月均值曲线 用本软件绘制，点击“资料日常处理”菜单下 择日期页面，在该页面下选定２００７年１月１ １３的 的“年度图一ｐ １３均值年度图”功能会出现选择 “ ．ｏｒｇ”数据文件后，弹出“是否要绘制误差限 ３１期页面，在该页面下选定“２００７年１月１日”的 线”页面，如果要绘制误差限线，选择“是”按钮即 “ ．ｏｒｇ”数据文件后即可得到如图７所示的２００７ 可得如图８所示的２００７年度带误差限线的月均 年度１３均值图形。点击“资料１３常处理”菜单下 值图形，图形绘制利用ＣＦｉｌｅ类的Ｒｅａｄ方法将所 的“年度图一ｐ 月均值年度图”功能后会山现选 需的数据从“ ．ｏｒｇ”数据文件中提取出来并转换 高 原 地 震 第２２卷 成十进制形式，然后利用ＣＣｌｉｅｎｔＤＣ类的ＭｏｖｅＴｏ 员，减小了误差。 和ＬｉｎｅＴｏ方法，在对话框窗口中绘制图形并显 ４．５查看日均值波动幅度 示，为了让图形具有保持功能，在对话框窗口的 例如查看２００７年度的日均值波动幅度。 ＰｒｅＴｒａｎｓｌａｔｅＭｅｓｓａｇｅ方法中调用ＣＣｈｅｎｔＤＣ类的 传统处理方法是通过转换２００７年原始 ＭｏｖｅＴｏ和ＬｉｖｅＴｏ方法。 “ ．ｏｒｇ”数据文件后，手工检索出公式（１５）所需 通过对比绘制可以看出，本软件较传统绘制 数据并且计算。 方法简洁、快速、方便。在绘制日均值图时直接利 用本软件查看，只需点击“工具”菜单栏下的 用“ ．ｏｒｇ”原始数据文件，运用自身的绘图功能 “正常动态指标计算”，然后在出现页面下的“开 绘制，不用像传统处理方法利用多个软件协同绘 始日期”和“结束日期”栏分别输入２００７年１月１ 制；在绘制带有误差限线段的月均值图时，直接将 日和２００７年１２月３１日的“％．ｏｒｇ”数据文件，接 手工换算绘制自动化，很大程度地方便了工作人 着点击“计算”按钮，即可得结果，如图９所示。 功青色区 ～ 一 一 一 ’ ～～ 九五ｏｒｇ文件 季孽 Ｊ芋试 盒 ｇ、晕００７ｏｒｇ、碡 ８建ｌ０１．ｏｒｇ穗 雾。 》 塑里苎 ｊ Ｉ 程 骚 ＼ －－－￣－ ＿ …结柬日期 ｌ… 一 Ｊ匡塞 团 ．堡 ～堡Ｊ—堡 ｊ 显示区一 一 一 一 …－ ・ …一 ……一…～…一 ． ，＾ ～—． ………－ 一 …… ～ 图９陇南台２００７年正常动态相关数据 通过对比可以看出，本软件可以很方便地利 地方便了工作人员，也可以为更进一步深入研究 用原始数据文件检索出一段时问内的视电阻率最 提供基础数据。 大值、最小值并给出这段时间视电阻率的均值以 （３）为基层台站提供了一个初步的地电阻率 及均值的均方根误差和视电阻率波动幅度。这对 资料质量自评价和快速判断资料变化的工具，为 基层台站初步评定一段时间的资料，及时发现资 地震预测部门可及时提供基础数据。 料异常提供了很大的便利。 本软件由于使用时间较短、范围有限，在全国 ５结论与建议 地电台站推广使用时出现的不足恳请批评指正， 有好的建议欢迎指导。 通过软件的试用和与传统处理方法的对比， 感谢甘肃省地震局嘉峪关中心地震台在软件 可以看出本软件具有如下特点： 的测试过程中提供的无私帮助！ （１）将地电阻率资料日常处理工作系统化， 将传统处理方法中大量繁琐的手工计算自动化， 参考文献： 在很大程度上提高了地电阻率资料日常处理工作 ［１］ 赵玉林，钱复业，薛顺章，等．地电阻率前兆场与应 的时效性和准确性。 变场关系及其预报地震方法的研究［Ａ］．石特临． （２）数据输入与结果输出灵恬方便，可以直 地震地电学发展与展望［Ｍ］．兰州：兰州大学出版 接产出地电台管组所需日常处理数据结果，极大 社．１９９８：４１—５０． １Ｊ］Ｊ １Ｊ １ｊ］Ｊ １Ｊ 第４期 狄国荣，等：地电阻率资料日常处理软件研制 ［８］ 钱能．ｃ＋＋程序设计教程（第二版）［Ｍ］．北京： 钱家栋．与大震孕育过程有关的地电阻率变化研究 ［Ｊ］．中国地震，１９９３，９（４）：５５—６４． 清华大学出版社，２００６：１０—２０． 钱复业，赵玉林．地震前地电阻率变化十例［Ｊ］． 地震学报，１９８０，２（２）：１８６—１９７． ［９］ 白中英．计算机组成原理［Ｍ］．北京：科学出版社， ２００３：１—１０． 陈彦平，狄国荣，梅东林，等．武都地电阻率资料在 文县一临江Ｍｓ５．０地震的异常反应［Ｊ］．地震地磁 观测与研究，２００９，３０（增刊）：１６一ｌ９． ［１０］ Ｂｊａｒｎｅ Ｓｔｒｏｕｓｔｒｕｐ．Ｔｈｅ Ｃ＋＋Ｐｒｏｇｒａｍｍｉｎｇ Ｌａｎｇｕａｇｅ （Ｓｐｅｃｉａｌ Ｅｄｉｔｉｏｎ）［Ｍ］．北京：Ａｄｄｉｓｏｎ Ｗｅｓｌｅｙ／Ｐｅａｒ－ ＳＯＩｌ，２００１：１—８． 高曙德，汤 吉，杜学彬，等．汶川８．０级地震前后 电磁场的变化特征［Ｊ］．地球物理学报，２０１０，５３ （３）：５１２—５２５． Ｒｏｂｅｒｔ Ｌ．Ｋｒｕｓｅ．Ａｌｅｘａｎｄｅｒ Ｊ．Ｒｙｂａ．Ｄａｔａ Ｓｔｒｕｃｔｕｒｅｓ Ａｎｄ Ｐｒｏｇｒａｍｄｅｓｉｇｎ Ｉｎ Ｃ＋＋［Ｍ］．北京：高等教育 出版社，Ｐｅａｍｏｎ Ｅｄｕｃａｔｉｏｎ出版集团，２００１：１５— ３０． 中国地震局．地震及前兆数字观测技术规范［Ｍ］． 北京：地震出版社，２００ｌ：ｌ０—３Ｏ． ［１２］ 侯俊杰．深入浅出ＭＦＣ（第二版）［Ｍ］．武汉：华 中科技大学出版社，２００１：１—２０． 谭浩强．ｃ程序设计（第二版）［Ｍ］．北京：清华大 学出版社，２ｏ００：１—６． ＴＨＥ ＳｏＦＴＷＡＲＥ ＤＥＶＥＬｏＰＭＥＮＴ ＦｏＲ ＤＡＩＬＹ ＰＲＯＣＥＳＳＩＮＧ ｏＦ ＴＨＥ ＤＣ ＲＥＳＩＳＴＩＶＩＴＹ ＭＯＮＩＴｏＩＵＮＧ ＤＡＴＡ ＤＩ Ｇｕｏｒｏｎｇ ，ＴＡＮＧ Ｊｉ ，ＧＡＯ Ｓｈｕｄｅ ，ＳＵ Ｙｏｎｇｇａｎｇ ，ＹＥ Ｙｕａｎｙｕａｎ ，ＣＨＥＮ Ｙａｎｐｉｎｇ （１．Ｌｏｎｇｎａｎ Ｃｅｎｔｅｒ Ｓｅｉｓｍｉｃ Ｓｔａｔｉｏｎ，Ｇａｎｓｕ Ｐｒｏｖｉｎｃｅ，Ｌｏｎｇｎａｎ ７４６０００，Ｃｈｉｎａ； ２．Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｇｅｏｌｏｇｙ，Ｃｈｉｎａ Ｅａｒｔｈｑｕａｋｅ Ａｄｍｉｎｉｓｔｒａｔｉｏｎ，Ｂｅｉｊｉｎｇ １ ０００２９，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：ＤＣ ｒｅｓｉｓｔｉｖｉｔｙ ｍｅａｓｕｒｅｍｅｎｔｓ ｐｌａｙ ａｎ ｉｍｐｏｒｔａｎｔ ｒｏｌｅ ｉｎ ｔｈｅ ｍｏｎｉｔｏｒｉｎｇ ｏｆ ｅａｒｔｈｑｕａｋｅ ｐｒｅｃｕｒｓｏｒｓ．Ｍａｎｙ ｏｂｓｅｒｖａｔｉｏｎｓ ｏｆ ＤＣ ｒｅｓｉｓｔｉｖｉｔｙ ｓｈｏｗ ｔｈａｔ ｔｈｅ ｒｅｓｉｓｔｉｖｉｔｙ ｉｓ ｃｈａｎｇｅｄ ｂｅｆｏｒｅ ｅａｒｔｈｑｕａｋｅ．Ａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ｔｈｅｏｒｉｅｓ ａｎｄ ｒｕｌｅｓ ｏｆ ｈｅ ｅｌｔｅｃｔｒｏｍａｇｎｅｔｉｃ ｏｂｓｅｒｖａｔｉｏｎ．ａ ｓｏｆｔｗａｒｅ ｕｓｉｎｇ Ｖｉｓｕａｌ Ｃ＋＋６．０ ｗａｓ ｄｅｖｅｌｏｐｅｄ ｆｏｒ ｔｈｅ ｄａｉｌｙ ｐｒｏｃｅｓｓｉｎｇ ｏｆ ＤＣ ｒｅｓｉｓｔｉｖｉｔｙ ｄａｔａ ｒｅｃｏｒｄｅｄ ｂｙ ＺＤ８ Ｂ ｏｒ ＺＤ８ ＢＩ ｉｎｓｔｒｕｍｅｎｔｓ．Ｔｈｅ ｂａｃｋｇｒｏｕｎｄ ｏｆ ｔｈｅ ｓｏｆｔｗａｒｅ， ａｎｄ ｉｔｓ ｄｅｓｉｇｎ ａｎｄ ｒｅａｌｉｚａｔｉｏｎ ｍｅｔｈｏｄｓ ａｒｅ ｉｎｔｒｏｄｕｃｅｄ ｉｎ ｔｈｅ ｐａｐｅｒ．Ｔｈｅ ａｕｔｈｏｒｓ ａｌｓｏ ｇｉｖｅ ｓｅｖｅｒａｌ ｅｘａｍｐｌｅｓ ｏｆ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎ ｉｎ ｔｈｅ ｐｒｏｃｅｓｓｉｎｇ ｏｆ ｔｈｅ ｄａｔａ ｍｅａｓｕｒｅｄ ａｔ Ｌｏｎｇｎａｎ ｃｅｎｔｅｒ ｓｅｉｓｍｉｃ ｓｔａｔｉｏｎ，Ｇａｎｓｕ ｐｒｏｖｉｎｃｅ．Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗ ｔｈａｔ ｉｔ ｉｓ ｖｅｒｙ ａｖａｉｌａｂｌｅ ａｎｄ ｃｏｎｖｅｎｉｅｎｔ ｆｏｒ ｔｈｅ ｄａｉｌｙ ｐｒｏｃｅｓｓｉｎｇ ｗｏｒｋ，ｉｔ ｃｏｕｌｄ ａｓｓｅｓｓ ｔｈｅ ｐｒｅｃｕｒｓｏｒ ａｎｏｍａｌｉｅｓ ｓｉｍｐｌｙ ａｎｄ ｐｒｏｖｉｄｅ ｔｈｅ ｄａｔａ ｔｏ ｔｈｅ ｓｅｃｔｉｏｎｓ ｏｆ ｅａｒｔｈｑｕａｋｅ ｐｒｅｄｉｃｔｉｏｎ ｑｕｉｃｋｌｙ． Ｋｅｙ ｗｏｒｄｓ：Ｅａｒｔｈ ｒｅｓｉｓｔｉｖｉｔｙ；Ｐｒｅｃｕｒｓｏｒｙ ａｎｏｍａｌｙ；Ｇｅｏ—ｅｌｅｃｔｒｉｃ ｉｎｓｔｒｕｍｅｎｔ；Ｓｏｆｔｗａｒｅ
