Characteristics and Background of Shandong Rushan Earthquake Sequences
-
摘要: 2013年10月1日在山东省威海市乳山市发生M3.2级地震，之后发生了一系列震群活动。截至2016年5月，山东台网已经记录到了1万多次余震，其中3级以上地震9次，4级以上地震3次。频繁的地震构造活动引起了乳山市及周边地区强烈震感。为研究乳山震群的发震机理，本文利用山东台网数字化地震波资料和新建的乳山台阵资料，通过双差精定位方法重新确定了震中位置。研究结果表明：余震序列呈现出NW向的条带分布；采用CAP方法（Cut and Paste）反演震群中9次3级以上地震的震源机制解，结果显示几次较大地震的震源深度平均约为5km，与台网编目定位的结果基本相同。从得到的精定位结果并结合震源机制解的结果来看，震群的走向是NW向，倾角是NE向，与最近的乳山断裂有一定距离。由此推断该区域可能是乳山断裂的分支，或者有一条或多条隐伏断裂。Abstract: An earthquake with MS3.2 occured in Weihai city of Shandong province on October 1st, 2013, followed by a series of earthquake sequences. Until now (May, 2016), the Shandong Seismic Network has recorded more than ten thousand aftershocks, among which there are nine events with magnitude over MS3.0 and three over MS4.0. Since then frequent seismotectonic activities have been perceived in Rushan and the surrounding areas. To understand the mechanism of Rushan earthquake swarms better, we collected and analyzed the digital seismic wave data from Shandong Seismic Network and the seismic array data from the newly-built Rushan station, and relocated the epicenter location and the sequence distribution with the Double-Difference method, and got a N-W stripe distribution from the aftershock sequence. The CAP (Cut and Paste) method was applied to inverse the focal mechanism solution three times of the earthquakes over MS4.0 and 6.0 times of the earthquakes over MS3.0-4.0. The results showed that the average focal depth of several great earthquakes was 5km, about the same with that of the seismic network. The MS3.2 earthquake in Rushan happened on January 7, 2014, was analyzed with the detailed result of nodal planeⅠ(strike angle 299.5°, dip angle 64ånd rake angle 16.7°) and nodal planⅡ(strike angle 202°, dip angle 75° and rake angle 162°). The strike distribution of the earthquake swarms is NW while the dip distribution NE. Because the certain distance between the distribution and the Rushan fault, it can be deduced that this area may belong to the part of Rushan fault, or there may be some undiscovered faults.
-
Key words:
- Rushan /
- Accurate relocation /
- Focal mechanism
-
引言
随着国家地震烈度速报与预警项目工程的实施，将在全国建设大量烈度计观测台站，弥补现有地震台网能力的不足，有效提高减灾和社会服务能力。其中，天津地区已在京津冀简易烈度计预警示范工程中先行建设了80个简易烈度计台站，与天津行政区内具备实时传输能力的测震台站和强震动台站共同组建天津地震预警观测系统（许可等，2019）。现有观测系统中缺少台站设备状态监控，台站各种设备基本处于未知状态，台网中心设备维护人员对台站设备状态的判断仅限于烈度计信号通断与否，台站出现故障后（如市电故障、电压不稳、网络故障、设备死机等），运维人员无法远程判断故障原因并进行有效处理，须到现场进行排查及维修，增加了运维成本，且效率较低。针对上述问题，设计烈度计台站远程监控系统，使台网和台站值班人员可在线实时查看仪器运行状态，及时发现各类设备故障，并对故障进行远程处理，提高台网管理与维护能力。
1. 监控系统的构成
烈度计台站远程监控系统物理架构如图 1所示，分为硬件设备和软件平台。硬件设备主要包括信息采集设备、服务器、PC和手机，其中，信息采集设备部署在烈度计台站，是整个监控系统的数据支撑，也是逻辑架构中的数据采集层。软件平台部署在监控中心，涉及业务层的使能平台和大数据页面、数据仓库层中的数据库、服务器中的MQTT中间件、手机APP等。
信息采集设备采集台站供电、网络、仪器状态等信息，将采集到的原始数据上传至使能平台进行存储，将原始数据传至解析器，解析器对电源数据、专业仪器数据、网络数据等进行类型分析，同时判断数据是否正常，并进行数据分类存储与统计分析，大数据界面通过GIS地图显示台站状态，对异常数据台站进行报警。运维人员查看报警台站实时数据信息，通过使能平台或手机APP向信息采集设备发送控制指令，对观测仪器等设备进行远程维护。
2. 信息采集设备
2.1 设备功能
烈度计台站信息采集设备通过对台站设备运行状态、供电状态、网络状态等进行全方位监测，及时发现并处理存在的问题，最大程度地减小设备离线率和故障率，确保台站设备安全可靠地运行。根据实际需求，本系统具有以下功能：
（1）基于NB-iot通信具有功耗低和费用低的特点，信息采集设备上行数据接口采用单独的NB-iot通信链路，信息采集设备采集台站监控信息，并通过NB-iot网络回传至台网监控中心的监控平台。
（2）信息采集设备可通过以太网口对台站现有连接网络设备（路由器和烈度计）进行网络通讯，发现问题及时报警。设备可通过以太网读取烈度计状态信息及实时数据信息，进而对专业设备状态进行监控（王建国等，2010）。
（3）信息采集设备支持接入、接出市电，并对市电电压、电流、功率、功率因数、频率、用电量等参数进行监测。同时信息采集设备不影响其他设备供电，即使信息采集设备出现故障，也不会影响其他设备的正常供电。
（4）设备具有多路I/O输入及1路RS485输入接口，支持其他设备接口接入，对设备运行状态进行监测，从而判断设备是否正常运行。
（5）信息采集设备设有1路继电器输出，将需控制设备的供电线路接入继电器常闭端，即可通过终端发送特定指令完成相应设备的关闭、打开和重启（陈吉锋等，2012）。
（6）信息采集设备内置电容，能保证台站在断电情况下短时间工作，将故障信息回传。
2.2 功能实现
设备核心功能模块主要包括核心嵌入式控制器（MCU）模块、NB通信模组模块、能耗测量模块、网络通讯模块。信息采集设备功能模块逻辑图如图 2所示，其中，主控模块是设备的核心，通过C语言和汇编语言实现设备核心算法与逻辑控制，完成对外围硬件的控制功能；电源模块为控制电路提供电源；数据采集模块实现不同通讯协议设备的接入，包括电能计量模块、I/O接口输入、RS485通讯接口等；以太网模块实现与本地路由器之间的通讯，通过TCP传输协议查询台站设备间的网络连通情况，包括烈度计、路由器等网络联通信息（宁晓青等，2019）；控制模块通过接收主控模块的指令实现设备断电、重启等操作；NB模块为通信模块，可使用3家运营商的NB服务，实现监控信息及控制指令的传输。
2.3 设备配置
上位机设备配置软件通过RS485接口与设备进行通信，实现NB-iot网络、以太网、测量量、烈度计等参数的配置，上位机配置软件界面如图 3所示。NB-iot网络参数配置实现设备与监控中心接收数据服务器的通信，需配置的主要参数包括MQTT服务器地址和端口、登陆ID及订阅主题；以太网参数配置实现设备与烈度计和路由器间的通信，需配置的主要参数包括设备网卡IP地址、网关地址和掩码地址；测量量参数配置实现台站网络与供电的监控测量，需配置的主要参数包括台站设备IP地址、环境参数（如电压、功率、温度等）；地震烈度计监测参数配置实现烈度计状态信息的监控测量，需配置的参数主要包括烈度计类型和IP地址。
3. 监控平台
监控平台是基于物联网侧设备接入使能的云化平台系统，能有效监控烈度计台站设备数据的采集、存储、分析、数据展现及发布、智能管控等，平台具备丰富的对外数据接口，可简单灵活地通过插件编程实现不停机对接收数据保存、解析、加解密、格式转换等。
3.1 逻辑架构
监控平台逻辑架构从功能层上分为数据仓库层、功能层和业务层（图 4）。
数据仓库层实现对数据存储表的管理，包括信息采集设备上传的原始数据表、解析数据表、用户管理表、业务报表、系统监控表、参数设备表。功能层包括监控平台实现的功能，如身份认证、数据解析、数据查询和报表分析等。业务层包括使能平台和大数据页面，使能平台主要完成数据查询统计和下行，大数据页面完成数据统计分析及展示。
3.2 功能模块
监控平台从使用上分为用户模块、设备模块、应用模块和解析器模块，其中，用户模块为基础，设备模块为纽带，应用模块为主干，解析器模块为重点，各模块功能见表 1。
表 1 监控平台各模块功能Table 1. Functions of monitoring platform modules模块 功能 用户模块 存储用户的基本信息，按登录用户信息显示不同的平台信息 设备模块 对上行数据进行分类，将数据与设备对应，对设备基本信息进行查询与维护 应用模块 作为设备与解析器的依托，控制设备离线监测、数据解析、设备报警等功能是否开启 解析器模块 将所有上行数据解析后展现在监控平台的大数据界面上 4. 系统应用
烈度计台站信息采集设备已在80个简易烈度计台站安装部署，监控平台在台网中心服务器进行部署，监控平台基于B/S架构进行设计，方便用户操作。值班人员通过监控平台大数据界面（图 5）查看台站运行状态，发现报警及时远程维护。监控系统在实际运行过程中多次监测到台站供电中断、网络故障和烈度计数据异常。当监测到台站供电中断时，运维人员第一时间给烈度计台站看护人员打电话确认供电故障的具体原因，确保供电故障及时修复；当监测到台站网络中断时，通过远程控制路由器重启解决由于路由器死机导致的网络故障；当监测到烈度计数据异常时，通过远程控制烈度计重启解决由于烈度计死机导致的数据异常问题。
监控平台还具有大数据统计分析功能，如可对台站报警类型及报警次数日排名、台站通信流量排名、台站电压日统计报表、报警次数月统计报表、台站报警类型占比、台站报警状态占比进行展示。运维人员可根据相关统计报表有针对性地对台站各类设备进行定向优化升级，保证观测数据的稳定可靠。
5. 结语
烈度计台站远程监控系统可实时监控烈度计台站各设备运行状态，设备如果出现故障可判断具体故障原因，并及时进行远程维护，减少运维成本，有效提高运维人员工作效率和监测数据的连续率，具有应用与推广价值。
-
图 1 乳山序列附近断层、发生位置和乳山台阵台站分布
F1乳山断裂；F2海阳断裂；F3朱吴断裂；F4米山断裂
Figure 1. Location of nearby faults, Rushan earthquake sequence and seismic stations
图 4 乳山震群地震序列台网原始震中分布图
Figure 4. Epicenter distribution of Rushan catalog earthquake sequences
图 5 乳山震群地震序列重新精定位后震中分布图
Figure 5. Epicenter distribution after relocation of Rushan earthquake sequences
图 6 乳山震群精定位深度剖面分布图
Figure 6. The depth profiles of Rushan earthquake sequences after relocation
图 7 乳山震群较大地震的震源机制解
蓝色沙滩球为3次4级以上震源机制解，红色填充的沙滩球为6次3—4级地震的震源机制解，白色区域的蓝点代表压缩P轴，红色和蓝色区域的黄点代表拉张T轴
Figure 7. Focal mechanisms of earthquakes with relatively large magnitude of Rushan earthquake sequences
图 8 CAP方法的最佳拟合深度和波形拟合
Figure 8. Tendency chart of focal depth variation with different depth errors and focal mechanism
图 9 乳山4.2级地震的理论地震图和观测地震图的对比
Figure 9. The comparison of theoretical seismogram and observed seismogram of Rushan MS4.2 earthquake
表 1 地壳速度模型
Table 1. Crustal velocity model
深度/km Vp/km•s-1 Vp/VS 0 4.00 1.73 5 6.10 1.73 10 6.20 1.73 20 6.40 1.73 30 6.78 1.73 40 8.20 1.73 
下载: 导出CSV
表 2 乳山震群9次3级以上地震震源机制解
Table 2. Focal mechanisms of nine events over MS3.0 of Rushan earthquake sequence
发震时刻 经度/°E 纬度/°N 震级/M 节面Ⅰ/° 节面Ⅱ/° 最佳拟合震源深度/km 走向 倾角 滑动角 走向 倾角 滑动角 2013-10-01 12：07 121.70 36.83 3.2 297 68 －13 32 78 －158 5.5 2014-01-07 22：24 121.67 36.82 4.2 299.5 64 16.7 202 75 153 8.6 2014-04-04 00：12 121.70 36.83 4.0 289 90 27 199 63 180 8.3 2014-07-16 00：40 121.69 36.84 3.0 260.1 55.1 －41.6 17 57 －137 4.9 2014-09-16 14：42 121.65 36.83 3.5 317.6 61 －61.2 89 40 －131 3.8 2014-09-16 14：43 121.67 36.83 3.3 270 69 －52 24.6 42.6 －148.1 4.1 2015-05-22 00：05 121.69 36.83 4.6 295 85 －19 27 71 －175 10.2 2015-06-09 22：29 121.69 36.84 3.1 290 90 22 200 68 －180 7.6 2015-07-01 19：31 121.67 36.80 3.0 281 71 36 178 56 157 6.3
下载: 导出CSV
-
蔡一川, 戴仕贵, 魏娅玲, 2015.四川木里小震群活动特征研究.震灾防御技术, 10(4):933-940. doi: 10.11899/zzfy20150411 胡幸平, 崔效锋, 2013.华北地区中部地震精定位与构造应力场研究.震灾防御技术, 8(4):351-360. doi: 10.11899/zzfy20130402 黄媛, 2008.结合波形互相关技术的双差算法在地震定位中的应用探讨.国际地震动态, (4):29-34. http://www.cnki.com.cn/Article/CJFDTOTAL-GJZT200804005.htm 李霞, 刘希强, 李亚军等, 2012.山东及邻区地震的重新定位及其与活动构造的关系.中国地震, 28(4):381-392. http://www.cnki.com.cn/Article/CJFDTOTAL-ZGZD201204005.htm 刘建达, 李清河, 张元生等, 2009.江苏中、南部及邻区地震重新定位及在断层研究中的应用.地球物理学进展, 24(1):58-66. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ200901005.htm 曲均浩, 刘瑞峰, 李金等, 2014.CAP方法反演2014年山东乳山M4.2、M4.0地震震源机制解.地震工程学报, 36(4):1076-1080. http://www.cnki.com.cn/Article/CJFDTOTAL-ZBDZ201404051.htm 杨歧焱, 孟立朋, 彭远黔等, 2015.河北省张家口市涿鹿MS4.3级地震震害调查及地震特征分析.震灾防御技术, 10(2):353-360. doi: 10.11899/zzfy20150215 张岭, 刘劲松, 郑天珧, 2005.利用Ground-Truth资料和Messy GA方法反演渤海湾盆地及其邻区地壳P波速度模型.地球物理学进展, 20(1):186-197. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ20050100W.htm 郑建常, 林眉, 王鹏等, 2015.CAP方法反演震源机制的误差分析:以胶东半岛两次显著中等地震为例.地球物理学报, 58(2):453-462. doi: 10.6038/cjg20150209 Helmberger D. V., Engen G. R., 1980. Modeling the long-period body waves from shallow earthquakes at regional ranges. Bulletin of the Seismological Society of America, 70(5):1699-1714. http://www.bssaonline.org/content/70/5/1699.abstract Waldhauser F., Ellsworth W. L., 2000. A double-difference earthquake location algorithm:method and application to the northern Hayward Fault, California. Bulletin of the Seismological Society of America, 90(6):1353-1368. doi: 10.1785/0120000006 Zhang H. J., Thurber C. H., 2003. Double-difference tomography:the method and its application to the Hayward Fault, California. Bulletin of the Seismological Society of America, 93(5):1875-1889. doi: 10.1785/0120020190 Zhu L. P., Helmberger D. V., 1996. Advancement in source estimation techniques using broadband regional seismograms. Bulletin of the Seismological Society of America, 86(5):1634-1641. http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.482.6036 -
下载:
点击查看大图
计量
- 文章访问数: 162
- HTML全文浏览量: 60
- PDF下载量: 11
- 被引次数: 0
