mcluster-微地震波形相关及震群自动识别软件
地震资料处理集群中CMU管理软件的应用
47SYS MANAGEMENT 系统管理摘要：在分析集群管理软件CMU功能及在地震资料处理集群中管理集群方法的基础上，通过实例介绍CMU在石油地震资料处理集群部署、监控、管理取得的效果，指出CMU在实际应用中的优点和存在的不足。
关键词：CMU；集群；监控一、CMU 简介CMU 是集群控制单元（Cluster Management Unit）的简称，与仅在集群基础层监控基础架构的工具不同，它可以对集群的工作负载和资源进行全面监控，包括节点硬件信息、运行状态、报警信息等，对集群系统各个节点状态一目了然。
单一的界面中能够同时动态监控1024个节点的状态，同时利用克隆技术对集群节点操作系统进行批量分发（见图1）。
图1 CMU 监控主界面CMU 的监控功能是通过收集每个集群节点上CMU 守护进程所采集的信息实现的。
将集群节点物理位置信息和节点的ip 地址等信息输入到CMU 系统的控制台中，能够自动扫描收集节点的硬件信息（CPU、内存、硬盘等）把数据集合和显示图形分成不同对象后重新处理、重新组织，以流式网络图形格式PNG（Portable Network Graphic Format）来显示各种不同的状态图，即可实现对集群的动态实时监控[1]。
主要组成部分如下所述：（一）守护进程CMU 守护进程部署在集群计算节点上的信息采集守护程序，通过SmallMonitoringDaemon 进程实现。
采集的信息包括计算节点的物理信息（如CPU，内存使用量）等监控信息，守护进程查看方法如下：[root@hc18n01a ~]# ps -ef |grep cmu/o p t /c m u /b i n /S m a l l M o n i t o r i n g D a e m o n -h 192.168.14.10 -o 48560 -O 49074 -i 48557 -a /opt/cmu/etc/ActionAndAlertsFile.txt -t 5000000 -M 192.168.14.254 -f 1 -s 1 -L 0 -r 0（二）CMU GUI 功能CMU GUI 功能是集群管理监控系统的核心，通过管理服务器中的JAVA 实现控制界面的呈现，系统管理员可以在控制界面中操作所有功能，如在node management 选项中添加需要监控的集群的名称、节点名称、ip 地址等，CMU 将实时监控的结果整理后通过图表方式分类显示给用户[1]。
Omega地震数据处理系统介绍
• •
同相轴拾取控制参数窗口
数据显示参数窗口
辅助实用程序
JOB_VU:Omega作业监测器---监视提交作业的状态，运行的环境，可对正在运行的作业进
行删除和查看宽行
TraceDump:Omega数据体道头信息卸载---卸载某一数据体的道头信息和查询数据体
处理历史
Xfedit:Omega文本编辑器---提供方便的文本编辑工具 Xshow:Omega的文件管理器---可以对Omega工区内的各种文件进行拷贝、
是承接野外与提供可应用的地震成果资料的 重要中间环节。
1、对野外采集的地震数据进行检查、分析 2、为下一步的地震资料解释提供可靠的基础资料
地震 资料 解释
进行构造解释和岩性等综合解释， 为地质综合研究 提供各种地震成果资料。
2、主要处理软件介绍：
Omega Grisys Focus Promax CGG 西方地球物理公司地震处理软件 中国东方地球物理公司处理软件 帕拉代姆地球物理公司处理软件 LANDMARK 处理软件 法国CGG公司地震处理软件
工区管理结构示意
IVP（交互速度分析）— 交互速度分析的图形界面（2）
速度工区列表
速 度 拾 取 、 编 辑 点 显 示
组名
交互速度分析结构示意
线号
主要功能：
Cmp点号
1、交互速度拾取、编辑 2、多种辅助数据显示，包括速度谱、CMP道集 、交互叠加剖面，MVF叠加剖面、层位时间、 邻近速度曲线、速度组、层速度显示等。 3、层位拾取、编辑功能 4、多种速度、层位QC手段
切除前后道集对比
DDI（数据驱动式交互处理）— 交互处理的图形界面
时窗定义
显示信息定义
起始时间定义
IMS微震系统介绍
详细产品概览IMS提供了数字化，智能化，高分辨率的地震监测和控制系统，具有在线地震信息处理，分析和可视化功能。
该系统易于使用，可在Microsoft Windows或Linux操作系统下运行。
除地震方面，许多非地震岩土工程传感器也可以用于监测。
当信号或某些参数超过阈值时，具有报警、控制和（或）停机功能。
该系统基于模块化设计，易于扩展，可从自记式监测单元扩展成连接数个台站的复杂网络。
并提供全天候24小时技术支持。
硬件概览IMS微震系统的硬件主要分为三个部分，即传感器，数据采集器和数据通信部分。
∙传感器将地面运动（地面速度或加速度）转换成一个可衡量的电子信号。
非地震传感器也可以用于IMS地震网络。
∙数据采集器负责将来自传感器的模拟信号转换成数字格式。
数据可以被连续记录采集，或采用触发模式，通过特殊算法来确定是否记录微震事件发生的数据。
∙地震数据同时被传输到一个中央计算机或本地磁盘以待储存或处理。
系统可以采用多种数据通讯手段，以适应不同的系统环境需要。
微震传感器微震传感器通过将地面运动（地面速度或加速度）转换成一个可衡量的电子信号来衡量微震活动。
由于信号在本质上是模拟信号，传感器必须被连接到一个数据采集装置，将其转换成数字格式以待被计算机读取输出。
所有IMS传感器都包含智能电子元件以提供传感器类型，序列号和方向标识。
此外，智能传感器能够产生内部的振动，以验证传感器的操作和检测安装后电缆布线是否正确。
传感器类型微震传感器的分类取决于所要监测的地面运动类型，即地面速度（检波器）或地面加速度（加速度计和FBA）；传感器的传感轴数量和传感器是否部署在岩石钻孔里或岩石表面。
每个传感器类型在幅度范围，频率范围，可靠性和成本方面等有不同的优势。
一个IMS微震监测系统可基于检波器，加速计和力平衡加速计的任意组合，并同时搭配单分量和三分量传感器。
三分量传感器能够提供最精确的信息数据采集单元NETADCIMS的netADC是24位，4 或8通道，低噪声的模拟-数字转换器（ADC），以太网接口。
MSC系列软件介绍
基于有限元技术的 CAE 软件，无论在数量、规模上，还是应用范围上都处于主要 地位，作为 20 世纪中期兴起的技术手段，有限元技术随着计算机技术的迅猛发展，得 到了飞速发展和广泛的应用。
目前，CAE 软件主要应用于汽车、航空、航天、电子、土木工程、通用机械、兵器、 核能、医疗器械、铁道、石油和化工等行业。在国外，汽车行业是有限元软件的主要应 用行业，其所涉及的专业领域相当广泛，并且应用历史长、应用成熟度高。
2.1 MSC 公司介绍
MSC.Software Corporation 创建于 19863 年，总部设在美国洛杉矶，是享誉全球最 大的工程校验、有限元分析和计算机仿真预测应用软件(CAE，computer aided computer) 供 应 商 ， 其 最 早 最 著 名 的 产 品 是 大 型 通 用 结 构 有 限 元 分 析 软 件 MSC.Nastran 。 MSC.Nastran 软件始终作为美国联邦航空管理局(FAA)飞行器适航证领取的唯一验证软 件。
-3-
所有这些都是可以实现的。 VPD 建立在企业现有的投资基础上，包括设计和分析/仿真软件、计算机硬件、样
微地震数据处理方法
微地震数据处理方法
微地震数据处理是指对微地震（M≤4.0）的地震数据进行分析
和处理的过程。
以下是常见的微地震数据处理方法：
1. 数据预处理：主要包括数据采集、去噪和标定。
数据采集是通过地震仪器记录地震波形数据，去噪可以去除测量中的噪声干扰，标定是对地震仪器进行校准，保证数据的准确性。
2. 地震事件检测：通过对地震数据进行分析，检测出地震事件发生的时间和位置。
常用的方法包括振幅和能量的阈值检测、频率域分析等。
3. 地震定位：将地震事件的震源位置精确确定。
常见的方法包括到时差定位和波形反演等。
到时差定位是根据地震波在不同测点的到时差来确定震源位置；波形反演则是通过对地震波形数据进行逆过程反演，得到地震事件的震源位置。
4. 地震分析：分析地震事件的震级、震源深度、震中距等参数。
常见的方法包括振幅衰减分析、震级估计、地震波形拟合等。
5. 数据解释：对处理后的地震数据进行解释，分析地震活动的规律和构造特征。
常见的方法包括地震活动性分析、应力场分析等。
以上是常见的微地震数据处理方法，根据具体情况和需求，可以选择适合的方法来处理微地震数据。
多道法横波分裂各向异性参数反演软件V10使用说明
多道法横波分裂各向异性参数反演软件V1.0使用说明软件著作权登记号：2013SR057960刘劲松中国科学院地质与地球物理研究所1.功能简介本软件采用多道法(multi-channel method)对横波分裂各向异性参数进行反演。
软件包含5个模块，lev模块根据给定的参数对地震事件进行筛选，可按照震中经纬度范围、震级大小、震中距范围、震源深度范围、射线参数范围等条件筛选数据。
mafsac模块根据震源和台站位置信息计算震相的走时，并标注到sac 数据头中；tramp模块根据指定的时窗，计算径向横向分量数据的振幅比（即分裂强度si，splitting intensity）; azwin模块按方位窗对数据进行平均处理； fitsin 模块根据方位角/分裂强度数据反演各向异性参数。
该软件可用于SKS/SKKS等震相的横波分裂分析，研究地幔上地幔各向异性。
采用多道分析法的处理流程如图1所示。
输入数据为sac格式的三分量波形数据。
首先用mafsac软件标注sac数据的震相到时，再用可视化软件筛选并调整所需数据的起止时窗(sac数据头中的A、F字段)，然后用tramp模块计算分裂强度si，其输出结果用azwin按时窗进行平均后，用fitsin模块反演得到各向异性参数。
tramp模块也可直接对计算得到的方位角、分裂强度数据进行拟合得到各向异性参数。
图1.反演处理流程示意图。
2．模块用法详解2.1 lev模块lev模块根据给定的参数对地震事件目录进行筛选，可按照震中经纬度范围、震级大小、震中距范围、震源深度范围、射线参数范围等条件对地震事件进行筛选。
模块通过命令行变量输入参数，无变量执行lev会显示程序帮助信息，内容如下：******************************************************************** List Event catalogue according to given parameters* Version 2009.5.21********************************************************************usage: lev elist [ll=] [m=] [ymd=] [p=] [d=] [dp=] [baz=] [rp=]* [ck=] [-v] [-oa|-ob|-on|-ot] [-rp] [rq= net=] [-eN]* [ts=] [off=]* elist event catalogue file.* -eN input catalogue format.* N=0 anss format(default)* N=1 hypodd phase head fmt.* N=2 zxf catalogue fmt.* N=3 NEIC catalogue fmt.* N=4 BKC catalogue fmt.* N=5 ChaiXi fmt.* ll=lat1,lat2,lon1,lon2 output events in the region.* m=m1,m2 output events that magnitude=[m1,m2].* d=d1,d2 output events that delta=[d1,d2].* baz=bz1,bz2 output events that baz=[bz1,bz2].* dp=dp1,dp2 output events that depth=[dp1,dp2].(km)* rp=rp1,rp2 output events that ray prameter=[rp1,rp2](s/deg)* ymd=yb,mb,db,ye,me,de output events within the time window.* yjd=yb,jdb,ye,jde output events within the time window.* p=lat,lon compute delta between each epicenter and* the specified point* ck=pfx check if the file pfx-froot.BHZ exist,output* a X infront of each line if not exist* -v verbose mode output additional info of epicenter.* -oa output in original catalogue format.* -ob output in BKC catalogue format.* -on output in navigation format.* -ot output in reftek time range format.* -od output begining time of each day,in NEIC format.* for use with ymd & yjd options only.* off=before,after specify time range in seconds(for -ot opt).* -rp compute ray pra.* rq=stn output in IRIS request format,this gives station* net=net this gives networkname* ts=0 time shift in hour.****************************************************************以下详细解释每个命令行参数的意义：elist 地震事件目录文件。
微震监测技术发展状况及各提供商技术特点简介
微震监测技术发展状况及各提供商技术特点简介周乐1. 题引：微震技术简介 (3)2. 国际微震监测技术主要提供商 (4)3. 各大公司软硬技术特点及其主要业务领域 (5)3.1. 矿震研究院INSTITUTE OF MINE SEISMOLOGY (IMS) (5)3.2. 加拿大ESG Solutions 公司(ESG) (6)3.3. 美国MicroSeismic公司 (7)3.4. 英国Semore Seismic公司 (8)3.5. 加拿大Microseismic Industry Consortium (mu-SIC) (9)3.6. 澳大利亚联邦科学与工业研究组织CSIRO (10)3.7. Schlumberger 公司 (11)3.8. 英国Applied Seismology Consultants公司 (12)4. 各大提供商硬件主要技术指标对比 (13)5. 各大主要厂商用户群区域及市场占有率分析 (14)6. 主要信息来源 (15)1.题引：微震技术简介微震监测技术用于监测岩体在变形和破坏过程中，裂纹产生、扩展、摩擦时内部积聚的能量以应力波的形式释放，产生微震事件。
1990年代以来，由于引入了现代计算机技术、现代通讯技术、GPS 授时定位技术、地震学相关知识，该技术取得了突破性进展。
借助专业化的数据处理软件，能够实现在三维空间中实时准确地确定岩体中微震事件发生的位置、量级，从而对岩体的变形破坏的活动范围、稳定性及其发展趋势做出定性、定量评价。
微震监测技术在工程上的应用领域主要是矿山微震监测，监测并定位不稳定岩体，做出预警，还可以预测岩体对工程开展的反应，进行灾害等级评估等。
在石油开采领域，通过监测水制压裂过程中储油或储气层产生的微震信号，可以帮助操作人员得到压裂过程信息，以指导并优化工程参数设置。
在土木工程领域，利用微震监测技术，还可以进行隧道稳定性监测、边坡稳定性和大坝稳定性监测等。
Mshale基本参数设置
Meyer软件Mshale基本参数设置模块使用方法目录一、Data Options 数据选择 (2)1.1 General 概况 (2)1.2 Fracture 裂缝 (3)1.3 Proppant 支撑剂 (3)二、Data Description 数据描述 (4)三、Wellbore Hydraulics 水力井筒数据 (4)3.1 General 概况 (4)3.2 Deviation 井斜 (5)3.3 Casing 套管 (5)3.4 Tubing 油管 (5)3.5 Restrictions 油管变径 (6)3.6 BHTP References 深度参考 (6)3.7 Profile 井筒结构示意图 (6)四、Zones 射孔段 (7)4.1 Perforations 射孔参数 (7)4.2 Pay Zone 射孔层 (7)4.3 Fracture Network Options 裂缝网络模型 (8)4.4 Interaction 裂缝相互影响 (11)4.5 Proppant Distribution 支撑剂分布模式 (12)4.6 Near Wellbore 近井筒效应 (12)五、Treatment Schedule 泵注程序 (12)5.1 Auto Design自动设计泵注程序 (12)5.2 Input人工输入泵注程序 (13)六、Rock Properties 岩石特性 (14)6.1 Insert from Database 从数据库选择 (14)6.2 Import Log 测井数据导入/计算 (15)七、Fluid Loss Data 液体漏失数据 (18)八、Proppant Criteria 支撑剂标准 (22)九、Heat Transfer 热传导效应 (22)打开Mshale模块，菜单栏选择“Data”进入数据编辑模块，该模块共分：Data Options（基础数据）、Data Description（概况描述）、Wellbore Hydraulics （水力井筒）、Zunes（目的层-射孔）、Treatment Schedule（泵注程序）、Rock Properties（岩石特性-产层）、Fluid Loss（液体漏失）、Proppant Criteria （支撑剂性能）、Heat Transfer（热传递效应）九个模块。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
微地震波形相关及震群自动识别软件V1.0使用说明软件著作权登记号：2012SR079382刘劲松中国科学院地质与地球物理研究所1.功能简介微地震波形相关及震群识别软件，可对微地震事件的数字波形计算互相关函数，根据波形的相似性对微地震事件进行分类，找出属于同一震群的地震事件，软件可同时计算两个地震事件在同一个台站上的走时差用于双差法精确定位。
软件由三个模块组成：mecor模块计算微地震波形的互相关函数；corcluster模块根据相关函数的最大值自动识别震群；dtcorcc 模块根据双差定位的事件对数据和mecor的输出结果计算地震事件的走时差。
软件采用快速付氏变换法计算互相关函数，用抛物线3点插值法精确确定相关最大值点的位置，时差计算精度可达到一个采样间隔以内。
2．用法详解2.1 mecor模块该程序利用快速付氏变换计算两个地震波形的互相关函数，输出相关函数的最大值以及最大值的时移位置，最大值位置采用抛物线3点插值法精确求解，精度可达到一个采样间隔以内。
程序的输入数据是sac格式的地震波形数据。
程序会根据给定的地震波形，自动完成波形的配对并计算互相关函数，配对是指各个事件所有可能的组合。
例如，如果给定了同一个台站上10个事件的波形记录，程序将自动对这10个事件进行组合配对，最终计算得到的互相关记录为10*9/2=45个。
输入数据的sac头中的O字段(发震时刻)必须有赋值，否则计算出的相关时移值没有意义。
程序通过命令行变量输入参数，无变量执行mecor会显示程序帮助信息，内容如下：**** Help information ****************************************************** mecor [-bo] [-oc] [-cs] [-poc] [tl=] [-pz|-sz] sfn1 sfn2 ....sfnN* -bo start from the origin time of sac header,default is 1st pt.* -oc output correlation function.* -cs only compute cross-station record.* -poc phase only correlation.* tl=0. time window length in seconds,0 means the whole trace.* -pz set P wave area to zero(before s arrival)* -sz set S wave area to zero(after s arrival)* sfnx event SAC trace filenames.***************************************************************************以下详细解释每个命令行变量的意义：-bo 波形起始位置从sac头的o字段所指的时间开始(即发震时刻)，缺省从第1个采样点开始。
-oc 输出相关函数。
输出格式为sac格式，文件名为cor.sfnx.sfny。
缺省不输出相关函数。
-cs 只计算不同台站间波形记录的互相关函数。
缺省只计算相同台站波形记录的相关函数。
-poc 采用相位相关算法。
tl=0. 波形时窗长度，单位秒，0表示从起始位置到波形结束。
-pz 将S波到时前的数据清零，要求sac数据头中t2标记为横波到时。
-sz 将S波到时后的数据清零。
sfnX sac数据文件名。
2.2 corcluster模块该模块根据mecor计算的地震波形之间的互相关系数最大值和给定的阀值，。
程序通过命令行变量输入参数，无变量执行corcluster会显示程序帮助信息，内容如下：**** Help information ****************************************************** corcluster [cth=] [efn=] cfn1 cfn2 ....* cth=0.6 correlation coefficient threshold.* efn event_name epi info file.* cfnx correlation data from mecor.***************************************************************************以下详细解释每个命令行变量的意义：cth=0.6 相关系数阀值，缺省为0.6。
efn= 地震事件名-事件id映射表文件。
cfn 包含事件相关系数信息的mecor输出文件，可以有多个。
Corcluster模块的输出直接写到标准输出。
2.3 dtcorcc模块该模块根据mecor计算得到的地震事件时差，替换双差定位软件生成的双差数据中的时差，用于精确定位。
程序通过命令行变量输入参数，无变量执行dtcorcc会显示程序帮助信息，内容如下：**** Help information ****************************************************** dtcorcc dt.ct corfn emfn [-p|-s]* dt.ct dt.ct file of hypodd program.* corfn correlation file(output of mecor program).* emfn event-name epi info mapping file.* -p assign P phases(default)* -s assign S phases***************************************************************************下面详细解释每个命令行变量的意义：dt.ct 双差定位软件生成的双差数据文件。
corfn mecor模块生成事件相关系数数据。
emfn 地震事件名-地震事件id映射文件。
-p 只赋值P波震相双差。
-s 只赋值S波震相双差。
3．用法举例及数据格式3.1 mecor用法举例下面的命令对当前目录下所有垂直分量波形记录计算相关函数，波形起始点从发震时刻开始，波形长度10秒，计算结果输出到文件cor.dat中。
mecor tl=10 -bo IG*.BHZ >cor.datcor.dat的部分内容如下：0.149 -1.135 5.649 90 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_298_18_10_26-m-00.BHZ0.268 -0.265 0.516 3 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_300_04_19_30-m-00.BHZ0.160 0.090 4.542 70 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_301_12_39_58-m-01.BHZ0.157 -0.650 4.580 73 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_301_16_24_34-m002.BHZ0.633 -0.005 0.282 5 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_305_22_12_54-m-01.BHZ0.218 0.055 5.143 95 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_306_05_50_16-m003.BHZ0.126 -0.155 4.714 75 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_308_22_20_30-m-02.BHZ0.142 0.860 4.938 74 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_309_22_29_22-m002.BHZ0.122 -1.580 4.790 74 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_310_22_29_09-m000.BHZ0.133 -1.615 4.726 74 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_310_22_30_09-m010.BHZ0.137 0.590 4.799 75 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_311_03_22_20-m001.BHZ0.747 -0.090 0.238 3 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_312_16_07_59-m001.BHZ0.203 -0.330 4.180 84 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_315_00_37_11-m005.BHZ0.129 1.440 4.175 65 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_315_19_47_35-m-03.BHZ0.165 -0.140 4.532 70 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_316_08_05_44-m-01.BHZ0.268 -0.090 1.108 7 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_316_22_14_00-m001.BHZ0.187 0.825 4.411 70 IG01-2004_297_11_07_50-m-02.BHZ IG01-2004_317_05_14_07-m-00.BHZ数据的第1列为互相关系数最大值；第2列为相关系数最大值的时移位置，单位为秒；第3列为两地震事件震中的距离；第4列为两事件的方位角之差。
第5-6列分别为两个地震波形记录文件名。
在用于双差定位时，第2列的时差为第2个事件的到时减第1个事件的到时。
3.2 corcluster用法举例下面的命令根据mecor计算得到互相关记录数据cor.dat，对地震波形的关联性进行分析，自动识别出属于同一震群的地震。
互相关系数阀值取0.9，输出结果定向到文件clust.09中。
corcluster efn=eventnm_newepi.dat cth=0.9 cor.dat >clust.09输出的数据按事件数的多少从小到大排序，下面是部分内容，包含了3个震群：cluster 72 6 events maxrd: 0.926 0.6702006_170_07_44_59-m003 cid:72 38.4445 91.0028 4.5572006_170_14_32_16-m003 cid:72 38.4455 91.0032 4.6162006_163_17_45_56-m-01 cid:72 38.4406 91.0050 5.2272006_170_07_53_34-m-03 cid:72 38.4375 91.0066 5.0762006_172_18_59_04-m001 cid:72 38.4406 91.0047 4.9702006_172_18_41_20-m-02 cid:72 38.4403 91.0046 4.839cluster 73 7 events maxrd: 0.629 0.6932006_252_06_23_51-m004 cid:73 38.3582 90.8311 5.2472006_252_06_40_37-m003 cid:73 38.3603 90.8376 5.4842006_252_06_43_24-m003 cid:73 38.3573 90.8336 5.5082006_252_07_08_52-m003 cid:73 38.3572 90.8382 4.9492006_252_07_11_15-m004 cid:73 38.3577 90.8371 5.5152006_252_06_57_37-m005 cid:73 38.3574 90.8380 5.6422006_252_06_32_39-m003 cid:73 38.3591 90.8374 5.488cluster 74 12 events maxrd: 0.854 0.5202006_199_01_45_17-m002 cid:74 38.3527 90.8589 7.5352006_199_06_43_56-m005 cid:74 38.3585 90.8548 7.3052006_199_06_47_50-m007 cid:74 38.3578 90.8554 7.3342006_199_13_09_14-m003 cid:74 38.3585 90.8530 7.5422006_199_22_01_50-m009 cid:74 38.3577 90.8552 7.3652006_199_09_10_02-m003 cid:74 38.3554 90.8563 7.6432006_199_09_11_22-m001 cid:74 38.3535 90.8573 7.6852006_199_10_14_04-m004 cid:74 38.3576 90.8550 7.5002006_199_12_50_09-m002 cid:74 38.3552 90.8570 7.5942006_199_13_55_29-m004 cid:74 38.3582 90.8565 7.2892006_199_22_26_11-m000 cid:74 38.3577 90.8571 7.8092006_199_10_28_12-m-01 cid:74 38.3535 90.8604 7.784输出结果中以cluster打头的行为震群头，后面的第1个数字为震群序号，第2个数字为该震群内包含的地震事件数目，maxrd后面的2个数值分别为震群内地震震源位置间的最大水平间距和最大深度间距。