地震预警数据流监控与分析

王红蕾; 朱宏; 毛国良; 李冬圣; 李振波; 蒋宏毅

doi:10.11899/zzfy20230220

地震预警数据流监控与分析

doi: 10.11899/zzfy20230220

王红蕾^{1, 2,},
朱宏³,
毛国良^{1, 2},
李冬圣^{1, 2},
李振波^{1, 2},
蒋宏毅^{1, 2}

1.
河北省地震局, 石家庄 050021
2.
河北红山巨厚沉积与地震灾害国家野外科学观测研究站, 河北邢台 055350
3.
天津市地震局, 天津 300201

基金项目: 中国地震局地震科技星火计划专项（XH23005A）；河北省地震科技星火计划攻关项目（DZ20200326044）；国家地震烈度速报与预警工程（河北子项目）（152004000000163041）

详细信息

作者简介:
王红蕾，女，生于1981年。高级工程师。主要从事地震预警及地震信息网络相关研究。E-mail：wanghlcea@126.com

计量
- 文章访问数: 256
- HTML全文浏览量: 24
- PDF下载量: 31
- 被引次数: 2
出版历程
- 收稿日期: 2022-06-24
- 刊出日期: 2023-06-30

Data Flow Monitoring and Analysis of Earthquake Early Warning System

Wang Honglei^{1, 2
,},
Zhu Hong³,
Mao Guoliang^{1, 2},
Li Dongsheng^{1, 2},
Li Zhenbo^{1, 2},
Jiang Hongyi^{1, 2}

1.
Hebei Earthquake Agency, Shijiazhuang 050021, China
2.
Hebei Hongshan National Observatory on Thick Sediments and Seismic Hazards, Xingtai 055350, Hebei, China
3.
Tianjin Earthquake Agency, Tianjin 300201, China

摘要

摘要: 为实现地震预警网络环境实时监控与分析，明确各业务系统之间的数据流向和范围，根据地震预警系统业务分区梳理关键业务、网段和终端，并通过数据流分析实现地震预警系统网络流量监控。通过数据包采集与解码技术，实现对主要业务系统关键指标的监控，从而对地震预警系统的网络访问性能、系统服务性能、应用响应性能等关键指标进行综合监控与分析，实现以业务为核心的网络环境监控与支撑。
- 地震预警 /
- 流量分析 /
- 数据包检测 /
- 运维监控
Abstract: This paper studies the data flow monitoring of the Earthquake Early Warning system based on the connection of the core components. We sort out the key servers, network segments and terminals of the EEW system, and then analysis the data flow, extract key indicators by data packet acquisition and decoding technology, so as to comprehensively monitor and analyze the key performance indicators such as network access performance, system service performance and application response performance of the full life cycle of EEW system.
- Earthquake early warning /
- Network traffic analysis /
- Packet analysis /
- Operation and maintenance

HTML全文

引言

地震动高频衰减在工程地震领域有着重要的应用，常被用于震源、地震动模拟和衰减关系研究等工作中。kappa值是描述地震动高频衰减参数之一，多位学者对其进行了研究，如Sun等（2013）选用汶川地震52个主震加速度数据，得到kappa值变化范围，认为kappa值与震源和传播路径有关。Anderson等（1984，1986，1991）发现在美国加州地区强震数据中，kappa值与震中距表现出正相关关系。李文倩（2014）基于华北地区近2000条数据，利用Anderson等（1984）的方法得到kappa值，认为kappa值与震级和震源深度无关，但离散性较大；朱百慧（2016）采用Anderson等（1984）经典计算方法、固定应力降计算方法及宽带反演计算方法，以日本东北地区为研究区域，利用K-NET和KiK-net强震记录，对kappa值与传播途径、震级及局部场地特征的关系进行分析，发现kappa值与震中距具有一定线性关系，随着震中距的增大kappa值呈增长趋势，kappa值与震级的线性关系不明显，随着V_S30的增大kappa值呈对数下降趋势。

鄂尔多斯块体是古克拉通块体，其内部结构稳定，但其外部地震频发，尤其是鄂尔多斯块体西缘地区，该区域属于亚洲东部“大三角”地震构造区东边缘（高祥林等，2010），“大三角”地震构造沿华南、鄂尔多斯西边缘向北延伸至贝加尔湖附近，曾被称为中蒙大陆中轴构造带，研究该区域具有深远意义。鄂尔多斯西缘地区由内蒙古强震动台网、宁夏强震动台网属地监控，该区域共有14个强震动台站，台站积累了真实可靠的地震数据，为研究提供了数据支撑。

本文以鄂尔多斯西缘地区（34°N—42°N，103°E—109°E）14个固定强震台站（包括8个内蒙古强震台站和6个宁夏强震台站）2007—2020年116条强震记录（1≤M_S≤8）为研究对象，利用Anderson等（1984）提出的经典算法，通过MATLAB软件计算kappa0值，讨论kappa值与震中距、kappa0值与V_S20、V_S30的关系。

1. 研究区域台站分布与地质构造特征

鄂尔多斯西缘地区地质条件复杂，北部临河断陷盆地西南方向与潮水—吉兰泰断陷盆地相连。南部有银川地堑、干盐池盆地，东部鄂尔多斯盆地是古老的克拉通地块，地块相对稳定。西北部地区为阿拉善块体，近年中小震频发。区域内有河套平原和宁夏平原，均为第四纪泥沙沉积覆盖的冲击平原。

临河断陷盆地自早第三纪晚期开始发育，新生界沉积厚度最大达10000 m以上，第四纪以来断陷活动依然保持了第三纪晚期的活动特征，周边山地继续抬升，断陷继承性沉降，第四纪厚度2400 m。吉兰泰盆地最初发育于早白垩纪，但晚白垩纪-早第三纪古新世时期抬生为陆地，始新世又重新开始发育，直到现在。该断陷盆地内新生界厚度约2 000 m，其中第四纪最大厚度约400 m，沉降中心靠近巴彦乌拉山东麓断裂。银川地堑是由喜马拉雅造山运动时期构造活跃的贺兰山褶皱带与鄂尔多斯盆地相对上升形成的，其为本区域最大的拉张型活动盆地，北起石嘴山，南至青铜峡，长160 km，最宽55 km，地堑形成于渐新世，其后持续下陷，沉积了巨厚的新生代地层，总厚度达7000 m，其中古近纪最大厚度3600 m，新近系最大厚度2500 m，第四纪最大厚度>1609 m。干盐池盆地为青藏高原东北边界海原断裂带内最大的拉分盆地，沉积物粒度呈向上逐渐变细的趋势，湖泊沉积为主，干盐池盆地主要受南—西华山与唐家坡—干盐池断层的控制，最深处可能>550 m。第四纪以来，构造运动主要以大面积的周期节奏性震荡上升为主，在西吉县形成了剥蚀构造山地、构造剥蚀红岩丘陵，剥蚀堆积黄土丘陵，侵蚀堆积河谷平原。西吉县大面积被第四纪覆盖。河套平原是冲击平原，河套地区北部、南部断层证明该区域存在断层下陷，黄河流经低洼地区带来泥沙淤积，形成河套平原。宁夏平原也是冲击平原，由于贺兰山与鄂尔多斯高原的抬升，宁夏平原断裂下陷形成地堑带，后由于黄河携带大量泥沙堆积形成宁夏平原，可知宁夏平原是由地壳运动、断裂下陷、河流流水搬运堆积而成。

2. 计算分析

2.1 数据选取

研究数据来自于内蒙古自治区地震局和宁夏自治区地震局提供的地震记录数据、钻孔资料、剪切波速测试资料、地质条件资料等，研究区域地震台站信息如表1所示，台站分布及地质构造如图1所示。

表 1 研究区域地震台站信息

Table 1. The research area seismic station information

序号	台站名称	经度	纬度	记录器型号	场地类型
1	巴彦浩特	105.7°E	38.8°N	ETNA	土层
2	巴彦木仁	106.7°E	39.9°N	ETNA	土层
3	宝丰	106.3°E	39.0°N	MR-2002	土层
4	磴口	106.9°E	40.3°N	ETNA	土层
5	干盐池	105.3°E	36.6°N	MR-2002	土层
6	公地	106.8°E	40.6°N	ETNA	土层
7	吉兰泰	105.7°E	39.7°N	ETNA	土层
8	临河	107.6°E	40.8°N	ETNA	土层
9	灵武	106.3°E	38.1°N	MR-2002	土层
10	前进农场	106.4°E	38.8°N	MR-2002	土层
11	青铜峡	106.0°E	38.0°N	MR-2002	基岩
12	沙海	106.9°E	40.9°N	ETNA	土层
13	乌海	106.8°E	39.4°N	ETNA	土层
14	西吉	105.4°E	35.5°N	MR-2002	土层

下载: 导出CSV

| 显示表格

图 1 台站分布及地质构造

Figure 1. The distribution of seismograph station and geologic structure

下载: 全尺寸图片幻灯片

2.2 计算方法

加速度反应谱法由Anderson等（1984）提出，是目前工程地震学领域中使用最广泛的方法，本文采用该方法计算kappa值，开展区域高频衰减特性研究。频率高于特定频率f_E的单对数坐标系加速度傅立叶谱幅值衰减规律可表示为：

$$ {{A}}\left( {{f}} \right)={{{A}}_0}\cdot{\rm{exp}}\left( { - {{{\text{π}} kf}}} \right),\;f > {f_{\rm{E}}} $$

(1)

式中，A₀为依赖于震源特征和传播距离的因子，k为谱衰减参数，f为地震动频率，f_E为谱幅值对数值近似直线下降的起始点。

一般来说，当频率达f_E后，谱幅值随着频率的增加呈指数下降。在线性-对数坐标系（即单对数坐标系）下，高频衰减参数k可用最简单的线性关系式近似得到：

$$ {{k}}= - {\rm{\lambda }}/{\rm{{\text{π}} }} $$

(2)

$$ {\rm{\lambda }}=\frac{{\ln \left( {{A_2}\left( f \right)} \right) - \ln \left( {{A_1}\left( f \right)} \right)}}{{{f_2} - {f_1}}}=\frac{{\ln \left( {\varDelta a} \right)}}{{\varDelta f}} $$

(3)

kappa值计算截止点为傅立叶谱与噪声谱重合的点或幅值平行起始点，利用最小二乘法进行线性拟合，得到的斜率即为kappa值（图2）。

图 2 Anderson等（1984）经典计算方法计算的kappa值

Figure 2. Anderson&Hough's the classical method calculates k

下载: 全尺寸图片幻灯片

为计算和确定kappa值与震中距R的关系，利用最小二乘法得到线性回归方程：

$$ {{K=}}{{{K}}_{{0}}}{+}{{{K}}_{\rm{R}}}\cdot{{R}}\left( {{s}} \right) $$

(4)

$$ {{\bar{{K}}}_{\rm{0}}}{\rm{=[kappa0}}\left( {{\rm{N - S}}} \right){\rm{+kappa0}}\left( {{\rm{E - W}}} \right){\rm{]/2}} $$

(5)

式中，R(s)为震中距，K₀和K_R为常数量，K₀为由台站附近近地表地下地质构造引起的高频衰减参数值，K_R为与区域相关的衰减参数值（Anderson等，1984）kappa（N-S）和kappa（E-W）分别为kappa南北向和东西向分量，$ {{\bar{\rm{K}}}_{\rm{0}}}$表示强震台南北向、东西向kappa0均值。

2.3 数据处理

利用MATLAB语言编辑计算处理全部程序，包括kappa值、kappa值与震中距关系、kappa0值、kappa0值与V_S20、V_S30关系计算程序。

（1）通过编程实现对强震观测记录数据的处理，包括基线调整、采用0.1～20 Hz的带通滤波（白立新等，2019）、数据初筛及人工复筛等工作（崔建文等，2007；万永革，2012）。

（2）分析单个台站多次地震事件，利用计算机程序获得每次地震每个分量的kappa值，进而利用每个分量的kappa值结合震中距，采用最小二乘法进行拟合，得到kappa值与震中距的线性回归方程（陶正如等，2015）。

（3）分别求出116条强震记录kappa0（N-S）与kappa0（E-W）值与震中距的回归方程，并计算拟合优度R²。计算kappa0（N-S）、kappa0（E-W）、$ {{\bar{{K}}}_{\rm{0}}}$值，对$ {{\bar{{K}}}_{\rm{0}}}$值与V_S20、V_S30进行拟合，得到拟合关系曲线。

2.3.1 kappa值与震中距的关系

kappa值与震中距的关系计算结果如表2、图3所示，由表2、图3可知，kappa值与震中距存在线性回归关系；kappa值在鄂尔多斯西缘大部分地区拟合优度值相对较高，为0.6～0.9，而西吉台站地区南北向、东西向拟合优度值分别为0.47、0.49，可能是由于台站记录地震数据较少、仪器噪声过高等原因导致kappa值相对离散，拟合优度较低。

表 2 kappa值与震中距拟合结果

Table 2. The relationship between kappa and epicentral distance

台站名称	方向	回归方程	K₀/s	拟合优度值/R²	$ {{\bar{{K}}}_{\rm{0}}}$/s	V_S20/m·s⁻¹	V_S30/m·s⁻¹
巴彦浩特	N-S	kappa=0.000 3R+0.0264	0.0264	0.9803	0.02530	343.05	370.2
巴彦浩特	E-W	kappa=0.000 3R+0.0242	0.0242	0.9854	0.02530	343.05	370.2
巴彦木仁	N-S	kappa=0.000 2R+0.0363	0.0363	0.8691	0.03485	353.4	389.4
巴彦木仁	E-W	kappa=0.000 2R+0.0334	0.0334	0.9280	0.03485	353.4	389.4
宝丰	N-S	kappa=0.000 1R+0.0687	0.0687	0.6055	0.06560	180.4	222.1
宝丰	E-W	kappa=0.000 1R+0.0625	0.0625	0.7876	0.06560	180.4	222.1
磴口	N-S	kappa=0.000 1R+0.0605	0.0605	0.8971	0.05845	231.7	261.7
磴口	E-W	kappa=0.000 1R+0.0564	0.0564	0.9425	0.05845	231.7	261.7
干盐池	N-S	kappa=0.000 1R+0.0531	0.0531	0.7542	0.05245	317.4	357.3
干盐池	E-W	kappa=0.000 1R+0.0518	0.0518	0.7194	0.05245	317.4	357.3
公地	N-S	kappa=0.000 2R+0.0541	0.0541	0.9790	0.05450	222.7	251.2
公地	E-W	kappa=0.000 2R+0.0549	0.0549	0.9730	0.05450	222.7	251.2
吉兰泰	N-S	kappa=0.000 2R+0.0526	0.0526	0.8812	0.05525	327.7	335.8
吉兰泰	E-W	kappa=0.000 2R+0.0579	0.0579	0.7826	0.05525	327.7	335.8
临河	N-S	kappa=0.000 1R+0.0530	0.0530	0.9977	0.05175	232.3	249.7
临河	E-W	kappa=0.000 1R+0.0505	0.0505	0.9960	0.05175	232.3	249.7
灵武	N-S	kappa=0.000 3R+0.0449	0.0449	0.9996	0.04415	292.6	325.1
灵武	E-W	kappa=0.000 3R+0.0434	0.0434	0.9957	0.04415	292.6	325.1
前进农场	N-S	kappa=0.000 3R+0.0637	0.0637	0.7536	0.06555	209.5	248.9
前进农场	E-W	kappa=0.000 3R+0.0674	0.0674	0.6815	0.06555	209.5	248.9
青铜峡	N-S	kappa=0.000 1R+0.0428	0.0428	0.4987	0.04315	219.0	260.0
青铜峡	E-W	kappa=0.000 1R+0.0435	0.0435	0.6287	0.04315	219.0	260.0
沙海	N-S	kappa=0.000 2R+0.0489	0.0489	0.9743	0.04870	188.3	216.5
沙海	E-W	kappa=0.000 2R+0.0485	0.0485	0.9710	0.04870	188.3	216.5
乌海	N-S	kappa=0.000 5R+0.0155	0.0155	0.7320	0.01545	352.7	390.5
乌海	E-W	kappa=0.000 5R+0.0154	0.0154	0.6702	0.01545	352.7	390.5
西吉	N-S	kappa=0.000 1R+0.0429	0.0429	0.4713	0.04240	195.1	228.7
西吉	E-W	kappa=0.000 1R+0.0419	0.0419	0.4909	0.04240	195.1	228.7

下载: 导出CSV

| 显示表格

下载: 全尺寸图片幻灯片

图 3 鄂尔多斯西缘地区14个台站kappa值与震中距拟合关系

Figure 3. The fitting relationship between kappa and epicentral distance of 14 stations in the western margin of Erdos region

下载: 全尺寸图片幻灯片

2.3.2 计算结果不确定性分析

通过大量计算得到研究区域116条地震记录地震动峰值加速度PGA及震中距，震级与震中距关系分布如图4所示，震级与峰值加速度PGA关系分布如图5所示。由图4、5可知，震级为M_s2.0～M_s8.0，震中距为0～1200 km，峰值加速度PGA为0～180 cm/s²，多数集中在0～60 cm/s²。

图 4 震级与震中距关系分布图

Figure 4. Distribution of earthquake magnitude and epicenter distance

下载: 全尺寸图片幻灯片

图 5 震级与峰值加速度PGA关系分布图

Figure 5. Distribution of earthquake magnitude and PGA(peak acceleration value）

下载: 全尺寸图片幻灯片

由于研究区域地震较少，覆盖层较厚，台站有效记录地震数据不足，且kappa值与震中距拟合关系受远震（汶川地震等）影响，故拟合结果存在一定不确定性。针对此问题，剔除了震中距500 km以上的地震。部分台站拟合结果偏差较大，以巴彦木仁、磴口台站为例进行分析，结果如表3、图6所示。由表3、图6可知，剔除汶川地震等远震后，拟合优度值明显下降，巴彦木仁台南北向拟合优度值由0.86911降至0.6186，东西向拟合优度值由0.9280降至0.7081；磴口台南北向拟合优度值由0.8971降至0.2447，东西向拟合优度值由0.9425降至0.6064。综上所述，kappa值与震中距线性拟合受远震的影响较大，因此计算kappa值时应尽量考虑震中距为500 km以内的地震数据。

表 3 巴彦木仁、磴口台站kappa值与震中距拟合结果

Table 3. Fitting results of kappa value and epicenter distance at Bayanmuren and Dengkou station

台站名称	方向	回归方程	K₀/s	拟合优度值/R²
巴彦木仁	N-S	kappa=0.000 3R+0.0275	0.0275	0.6186
巴彦木仁	E-W	kappa=0.000 3R+0.0266	0.0266	0.7081
磴口	NS	kappa=0.000 1R+0.0650	0.0650	0.2447
磴口	EW	kappa=0.000 1R+0.0562	0.0562	0.6064

下载: 导出CSV

| 显示表格

图 6 巴彦木仁、磴口台站kappa值与震中距拟合关系

Figure 6. Fitting relationship between kappa and epicenter distance at Banyanmuren and Dengkou station

下载: 全尺寸图片幻灯片

2.3.3 kappa0值与剪切波速关系

研究区域有14个强震动台站，内蒙古强震动台网8个台站分别为：巴彦木仁台、巴彦浩特台、磴口台、公地台、吉兰泰台、临河台、沙海台、乌海台，宁夏强震动台网6个台站分别为：宝丰台、干盐池台、灵武台、前进农场台、青铜峡台、西吉台。计算得到研究区域内14个台站的kappa0值，结合合站提供的钻孔资料，得到剪切波速V_S20、V_S30值，进而得出台站kappa0值与剪切波速V_S20拟合关系式为：

$$ {\rm{kappa}}0 = - 0.035\;89 \times {\rm{ln}}\;{V_{{\rm{S}}20}} + 0.245\;8,\;{{{R}}^2} = 0.392\;9 $$

(6)

kappa0值与剪切波速V_S30拟合关系式为：

$$ {\rm{kappa}}0 = - 0.042\;38 \times {\rm{ln}}\;{V_{{\rm{S}}30}} + 0.286\;8,\;{{{R}}^2} = 0.414\;4 $$

(7)

kappa0值与剪切波速V_S20、V_S30拟合关系曲线如图7所示。由式（6）、式（7）、图7可知，kappa0值与剪切波速V_S20、V_S30存在对数关系，拟合优度值约为0.4，拟合效果不理想，结合朱百慧（2016)的研究成果，笔者认为kappa0值随着V_S20、V_S30的增大逐渐减至0（Ktenidou等，2013）。

图 7 kappa0值与剪切波速V_S20、V_S30拟合关系曲线

Figure 7. The fitting relationship between kappa0 and shear-wave（V_S20、V_S30）velocity

下载: 全尺寸图片幻灯片

3. 结论

本文利用Anderson等（1984）提出的经典方法，通过计算得出以下结论：

（1）鄂尔多斯块体西缘地区强震台站kappa值与震中距呈线性相关，且拟合优度值较高；

（2）kappa0值为0.01545～0.06560 s；

（3）kappa0值与V_S20、V_S30呈指数相关关系，离散性较大，拟合优度值约为0.4；

（4）随着V_S20、V_S30的增大，kappa0值逐渐减至0。

（5）本研究采用的地震数据较少，分析得到的kappa值与震中距线性拟合受远震的影响较大，因此计算kappa值时应尽量考虑震中距500 km以内的地震数据。

致谢本文在撰写过程中得到中国地震局工程力学研究所研究员陶正如、新疆地震局工程师李文倩、广东地震局高级工程师吴华灯、云南地震局研究员崔建文、防灾学院教授万永革、宁夏地震局高级工程师吕俊强及工程师李鸿廷的悉心帮助，在此向他们表示衷心感谢，同时感谢评审专家提出的宝贵意见。

图 1 省中心地震预警网拓扑

Figure 1. Topology of earthquake early warning network for provincial center

下载: 全尺寸图片幻灯片

图 2 地震预警业务流程

Figure 2. Diagram of earthquake early warning data flow

下载: 全尺寸图片幻灯片

图 3 整体业务流程

Figure 3. EEW system flow chart

下载: 全尺寸图片幻灯片

图 4 决策平台数据流交互

Figure 4. Decision modular data flow diagram

下载: 全尺寸图片幻灯片

图 5 发布平台数据流交互

Figure 5. Message sending modular data flow diagram

下载: 全尺寸图片幻灯片

图 6 数据流分发告警规则

Figure 6. Alarm rules diagram of data flow distribution

下载: 全尺寸图片幻灯片

表 1 核心应用

Table 1. List of core applications

业务分区	分系统部署单元	部署方式	IP地址
数据接入区	波形流服务（分流服务）	物理机	10.×.×.5（国家项目台站流）
			10.×.×.6（省项目台站流1）
			10.×.×.7（省项目台站流2）
			10.×.×.8（唐山台站流）
			10.×.×.9（张家口台站流）
			10.×.×.10（背景场台站流）
			10.×.×.11（一般站）
汇集交换区	波形流服务（综合）	物理机	10.×.×.12
	波形流服务（触发流）	物理机	10.×.×.13
	波形流服务（台网中心）	未知	10.×.×.14
	波形流服务（备份中心）	未知	10.×.×.15
	波形流服务（北京局）	未知	10.×.×.16
	波形流服务（天津局）	未知	10.×.×.17
业务处理区	预警处理软件（JEEW）	物理机	10.×.×.18
	预警处理软件（EEW）	物理机	10.×.×.19
	烈度速报	物理机	10.×.×.20
	决策平台	物理机	10.×.×.21
	决策平台（台网中心）	未知	10.×.×.22
数据存储区	数据存储AWS	物理机	10.×.×.23
信息发布区	信息发布子系统	集群	10.×.×.24；10.×.×.25
	消息中间件	主备	10.×.×.26
	消息通道：预警终端政府	集群	10.×.×.27
	消息通道：预警终端学校	集群	10.×.×.28
	消息通道：PC端、通信运营商、新闻媒体通道、特殊行业通道	集群	10.×.×.29
	消息通道：市级转发通道	集群	10.×.×.30；10.×.×.31
	Nginx（负载均衡）	主备	10.×.×.32
	MYSQL（紧急地震信息服务）	主备	10.×.×.12
	预警转发平台（IPTV、广电等预警信息发布）	集群	10.×.×.33
运维监控区	站点监控	集群（负载均衡）	10.×.×.34
	业务系统监控/预警终端	集群（负载均衡）	10.×.×.35
	观测数据监控	主备	10.×.×.36
	NTP	单机	10.×.×.37

下载: 导出CSV

表 2 预警业务网段划分

Table 2. Network segments of EEW system

地址段划分	业务分区	VLAN	网关
10.×.0.0/24	专线预警发布区（接入政府专线）	2100	10.×.0.254
10.×.1.0/24	互联网对外发布区	2101	10.×.1.254
10.×.2.0/24	DMZ区	2102	10.×.2.254
10.×.7.0/24	预警数据交换区	2107	10.×.7.254
10.×.8.0/24	预警数据处理区	2108	10.×.8.254
10.×.9.0/24	预警数据存储与备份区	2109	10.×.9.254
10.×.10.0/24	预警运行监控与展示区	2110	10.×.10.254
10.×.11.0/24	预警数据信息服务区	2111	10.×.11.254
10.×.13.0/24	预警系统安全管理区	2113	10.×.13.254
10.×.63.0/24	测试区	2163	10.×.63.254

下载: 导出CSV

表 3 预警台站网段划分

Table 3. Network segments of EEW stations

地市	台站数	网段	总带宽/Mbps
廊坊	9	15.×.4.0/22	18
保定	17	15.×.8.0/22	34
邢台	16	15.×.12.0/22	32
唐山	13	15.×.16.0/22	26
沧州	13	15.×.20.0/22	26
邯郸	15	15.×.24.0/22	30
衡水	11	15.×.28.0/22	22
承德	5	15.×.32.0/22	10
石家	4	15.×.36.0/22	8
张家口	34	15.×.40.0/22	68

下载: 导出CSV

表 4 监控指标

Table 4. Network monitoring metrics

监控类别	监控指标名称
网络性能	数据流量（上行和下行）	3次握手平均RTT （客户端、服务端）	ACK时延（客户端、服务端）	丢包率	重传率	丢包重传比
应用性能	响应时间	超时比例	交易次数	交易响应率
主机性能	会话数量	连接状态	服务器窗口尺寸	客户端窗口尺寸

下载: 导出CSV

表 5 主要应用监控指标

Table 5. Monitoring metrics of core applications

类别	监控指标	波形流服务（综合）	台网中心流交换	广东流服务（备份）	北京流交换	波形流服务（分流）	流服务（一般站）	决策系统
连接类型		长连接	长连接	长连接	长连接	长连接	短连接	短连接
网络性能	比特率/Mbps	142.84	21.85	22.49	19.30	6.20	9.54	0.064
	上行比特率/Mbps	1.79	20.07	20.45	16.28	2.14	8.37	0.032
	下行比特率/Mbps	141.05	1.79	2.05	3.02	4.06	1.18	0.032
	客户端平均ACK时延/ms	2.85	1 693.50	0.00	1.86	1.27	282.38	29.08
	服务器平均ACK时延/ms	1.18	15.95	39.41	0.00	26.04	34.30	33.10
	TCP重传率/%	0.09	0	0.03	0	0	0.39	0
应用性能	平均响应时间/ms	0.39	0	0	0	0	1.48	0.31
应用性能	服务器响应平均传输时间/ms	0.05	0	0	0	0	137.73	1.13
主机性能	活动会话数	10	1	1	1	134	552	10
	客户端平均窗口大小/KB	662.78	1.00	1.00	164.92	9.38	18.82	382.75
	服务器平均窗口大小/KB	37.33	5.45	25.48	64.19	8.32	62.77	56.84
注：监控时间段为2022年3月15日08:05—2022年3月15日08:10。

下载: 导出CSV

表 6 台站链路监控指标

Table 6. Monitoring metrics of station links

地市	台站数	数据流量/Kbps	上行流量/Kbps	上行丢包数/个	TCP重传率/%	客户端3次握手平均RTT/ms	服务器3次握手平均RTT/ms	平均ACK 时延/ms
廊坊	9	161.78	141.77	1 386	0	13.33	0.01	20.73
张家口	34	492.07	453.99	44	0	16.70	75.95	38.07
保定	17	324.12	288.55	2 278	0.01	71.19	1.29	19.47
邢台	16	285.91	250.55	2 856	0	8.11	0.17	20.54
唐山	13	251.64	231.60	1 862	0	165.34	1.29	25.08
沧州	13	229.89	203.44	9 675	0.01	45.67	0.59	20.24
邯郸	15	218.83	201.43	0	0	82.60	0.22	39.12
衡水	11	207.44	184.05	1 539	0	4.24	0.43	19.72
承德	5	89.99	83.03	217	0.01	11.66	4.62	27.13
石家庄	4	81.67	73.84	549	0.02	1.27	0.28	18.11

下载: 导出CSV

参考文献(16)

陈军华, 王忠民, 2006. BGP/MPLS VPN实现原理. 计算机工程, 32（23）: 124—126

Chen J. H. , Wang Z. M. , 2006. Theory and implementation of BGP/MPLS VPN. Computer Engineering, 32(23): 124—126. (in Chinese)

陈琳, 南洋, 2017. 一种面向大规模系统域网络性能管理系统. 计算机工程与科学, 39（9）: 1588—1593

Chen L. , Nan Y. , 2017. A performance management system for large-scale system area networks. Computer Engineering & Science, 39(9): 1588—1593. (in Chinese)

凯文 R. 福尔, W. 理查德·史蒂文斯, 2016. TCP/IP详解-卷1: 协议. 吴英, 张玉, 许昱伟, 译. 北京: 机械工业出版社, 423—427

Fall K. R., Stevens W. R., 2016. TCP/IP illustrated-volume 1: the protocols. Wu Y., Zhang Y., Xu Y. H., trans. Beijing: China Machine Press, 423—427. (in Chinese)

刘恩海, 刘微, 樊世燕等, 2015. 增强的NetFlow数据采集系统设计与实现. 计算机工程与设计, 36（5）: 1150—1155

Liu E. H. , Liu W. , Fan S. Y. , et al. , 2015. Design and implementation of enhanced Netflow data collection system. Computer Engineering and Design, 36(5): 1150—1155. (in Chinese)

孙路强, 张春莉, 徐小远等, 2020. 地震综合信息多渠道发布系统设计与实现. 震灾防御技术, 15（3）: 563—570

Sun L. Q. , Zhang C. L. , Xu X. Y. , et al. , 2020. Design and implementation of multi-channel release system of seismic comprehensive information. Technology for Earthquake Disaster Prevention, 15(3): 563—570. (in Chinese).

王红蕾, Walter M. , 高景春等, 2019. 美国地震预警系统ShakeAlert测试与评估. 国际地震动态, （5）: 17—25

Wang H. L. , Walter M. , Gao J. C. , et al. , 2019. Instruction of earthquake early warning ShakeAlert testing platform. Recent Developments in World Seismology, (5): 17—25. (in Chinese)

王俊, 刘红桂, 周昱辰, 2021. 地震预警技术的应用与展望. 防灾减灾工程学报, 41（4）: 874—882

Wang J. , Liu H. G. , Zhou Y. C. , 2021. Application and prospect of earthquake early warning technology. Journal of Disaster Prevention and Mitigation Engineering, 41(4): 874—882. (in Chinese)

王颀, 刘丰, 张彬等, 2021. SNMP与NetStream流量监控方法的对比分析. 海洋信息, 36（3）: 18—25

Wang Q. , Liu F. , Zhang B. , et al. , 2021. Comparative analysis of SNMP and netstream traffic monitoring methods. Marine Information, 36(3): 18—25. (in Chinese)

许可, 郭巍, 高也等, 2019. 天津简易烈度计地震预警试验区建设. 震灾防御技术, 14（2）: 456—463

Xu K. , Guo W. , Gao Y. , et al. , 2019. Development of simple intensity meter in Tianjin early warn experimental site. Technology for Earthquake Disaster Prevention, 14(2): 456—463. (in Chinese)

张红才, 金星, 李军等, 2013. 地震预警系统研究及应用进展. 地球物理学进展, 28（2）: 706—719

Zhang H. C. , Jin X. , Li J. , et al. , 2013. Progress of research and application of earthquake early warning system (EEWs). Progress in Geophysics, 28(2): 706—719. (in Chinese)

Allen R. M. , Melgar D. , 2019. Earthquake early warning: advances, scientific challenges, and societal needs. Annual Review of Earth and Planetary Sciences, 47: 361—388. doi: 10.1146/annurev-earth-053018-060457

Ganguli S., Bhalla V., Prasad P., 2018. Gartner magic quadrant for network performance monitoring and diagnostics. (2018-02-21)[2022-04-21].https://www.gartner.com/en/documents/3859467.

Given D. D., Allen R. M., Baltay A. S., et al., 2018. Revised technical implementation plan for the ShakeAlert system—An earthquake early warning system for the West Coast of the United States. Reston: USGS, 15—17.

Peng C. Y. , Jiang P. , Ma Q. , et al. , 2022. Chinese nationwide earthquake early warning system and its performance in the 2022 Lushan M6.1 earthquake. Remote Sensing, 14(17): 4269. doi: 10.3390/rs14174269

Saito H. , Chusho T. , 2005. Development of performance evaluation system for network systems by client observation. Electronics and Communications in Japan (Part I: Communications), 88(2): 71—80. doi: 10.1002/ecja.20083

Siswanto A., Syukur A., Kadir E. A., et al., 2019. Network traffic monitoring and analysis using packet sniffer. In: 2019 International Conference on Advanced Communication Technologies and Networking (CommNet). Rabat: IEEE, 1—4.

施引文献