基于支持向量机的高速铁路地震预警震级连续预测方法
地震是1 种对铁路行车安全危害极大的突发性自然灾害，即使是较小震级的地震，对路基、轨道和桥梁等造成的冲击都可能导致危害旅客生命安全的重大事故［1-4］。高速铁路地震预警系统是减轻地震带来重大损失的有效手段之一。如果能够在地震的破坏性震动到来前，对高速铁路提早哪怕是短短的几十秒甚至几秒实施报警并开展紧急处置，都能够大大降低旅客生命、财产损失的发生概率［5］。随着我国高速铁路的迅猛发展，建设高速铁路地震预警系统并提升地震预警能力迫在眉睫，同时，相关研究也将为穿越多个地震断裂带的川藏铁路设计建设提供支持参考。
综上所述，GDM双胎妊娠属于高危妊娠。加强GDM双胎妊娠的管理，控制血糖，选择合适的分娩方式，做到早发现、早应对，进而改善围产结局。
震级预测是高速铁路地震预警的重要环节，地震预警信息发布、地震影响范围判断、紧急处置控车范围确定等都依赖震级预测的结果。震级预测通常利用P 波初期特征参数与震级的线性比例关系，并依据这个比例关系建立预测震级的模型。预测模型中，采用的特征参数主要分为2 大类，即周期类参数［6-10］与幅值类参数［11-13］，基于这2 类参数的常用震级预测模型分别是τc方法［6］与Pd方法［11］。
在基于网状密集布设地震监测台网的地震预警系统中，震级估计可以采用加权平均或多台信息协同处理等方式。与地震行业不同的是，因铁路具有线性特征，高速铁路地震监测台站只能在铁路沿线呈稀疏线性布设，受台站数量、分布等方面的限制，高铁地震预警只能采用单台震级估计的模式，因此需要建立准确性更高的震级预测模型。
既有高速铁路地震预警通常利用地震P波到达后3.0 s 的固定时间窗进行震级预测，而常用的震级预测模型τc方法与Pd方法，都只用到地震波初期的单一参数特征，导致预测模型的泛化能力低，预测结果的离散性大，小震高估与大震低估现象明显，且建立特征参数与震级的线性比例关系时需要筛选震中距与信噪比。随着大数据分析技术的快速发展，机器学习方法在地震预警领域的应用逐步扩大，一些学者、专家将机器学习方法应用于地震预警震级预测的研究中，例如，马强［7］基于人工神经网络建立了多个参数与震级的关系，Reddy等［14］使用支持向量机的方法建立了小波系数与震级的关系，这些工作为利用多参数输入与人工智能方法进行震级预测提供了可行性参考。
“天命之谓性，率性之谓道，修道之谓教。道也者，不可须臾离也……”在全校师生齐诵《中庸》中纪念孔子诞辰暨新生开笔礼活动圆满落幕。
为尽早发出准确的地震预警警报，满足高速铁路地震预警警报随时间变化而不断更新的要求，实现“在3.0 s 以内提高震级预测的准确性、3.0 s 以后提高震级预测的连续性”的目标，本文选取P波到达后的0.5～10.0 s 范围内，以0.5 s 为间隔时间建立预测时间窗，利用日本K-net 强震动数据，基于人工智能机器学习领域经典的支持向量机（Sup⁃port Vector Machine，SVM）方法，以幅值参数、周期参数、能量参数、衍生参数这4 大类共计12 个P 波特征参数作为SVM 的输入参数，构建支持向量机的高速铁路地震预警震级预测模型（SVM-based High-speed Railway Magnitude Pre⁃diction Model，SVM-HRM），将3.0 s 时间窗下的震级预测结果与τc方法与Pd方法分别对比，并依据现行的《高速铁路地震预警监测系统试验方法》的相关条文［15］，将统计得到的震级预测实现率与实现率标准进行对比，旨在为中国高速铁路地震预警系统的建设与完善提供借鉴与参考。
语文学习是一个知识积累的过程，只有通过知识的长期积累，学生才能有语感，才能做到厚积薄发。通过课内外的阅读，学生才能积累到大量的语文知识，久而久之，语文素养也能得到有效提升。所谓“熟读唐诗三百首，不会作诗也会吟”，就是知识积累的写照。
1 数据及处理
机器学习算法需建立在大数据统计的基础上。中国现有的高铁地震监测台站记录到的数据较为有限，而日本的K-net 强震观测台网采用了与前者类型相同的强震观测台站，且记录有海量的高质量强震数据，因此本文的研究主要选用日本K-net 强震观测台网［16］数据，并基于以下原则选取强震动数据［17］。
（1） 发震时刻：2007年10月1日至2017年9月1日。
（2） 发震区域：日本岛内及周边海域。
（3） 震源深度：10 km范围以内。
（4） 地震震级：范围在3～8级之间。
（5） 震中距：未做筛选。
对选取的数据进行如下处理。
第四，重视史料教学,增强学科教学能力,应立足于历史学科的主要知识,充分挖掘史料实物、史象、新文本等新的史料信息库。这些新的历史资料用于历史教学实践,突出了历史变化的本质,提高了学生获得、整合和解读历史事迹的能力、解读历史现象、概括历史问题的能力,同时还可以提升高考命题考试的准确目的,结合高考问题中进一步把握新课程的方向和脉搏,从而优化我们的历史教学设置,创新一种面临高考历史自主教学的方法模式。
（1）采取马强等［18］与王子珺等［19］提出的方法，对加速度数据进行P波到时自动捡拾，并做人工核对。
（2）对加速度数据进行一次积分得到速度记录，对速度记录进行一次积分得到位移记录，对积分后的记录进行0.075 Hz 巴特沃斯高通滤波，消除积分带来的低频漂移影响。
（3）考虑到特征参数值巨大变化产生的数值失真，以及模型的训练效率提升等问题，取P波到达后0.5～10.0 s 的时间范围内，以0.5 s 为间隔计算特征参数，对数据依次进行归一化处理，归一化方法可表示为
其中，
式中：bg为各个特征参数归一化的结果；x 为各个特征参数对应的数据；xmax和xmin分别为相应P波特征参数的最大值和最小值。
经过处理，共筛选出地震1 837 次，强震动数据19 263组、57 789条，记录到强震动数据的台站922 个，筛选得到的地震震中与记录台站分布如图1所示。图中，红色空心圆表示地震震中位置，圆的直径与震级大小成正比；绿色三角形表示记录到数据的台站。
图1 选取的日本地震震中及K-net台站分布
将筛选得到的强震动数据随机划分为互不重复的2 组，即训练数据集和测试数据集，其中训练数据集占全部数据的80%，共有数据15 410 组、46 230 条；测试数据集占全部数据的20%，共有数据3 853 组、11 559 条。筛选数据的震级、震中距和记录数量关系如图2所示，图中绿色圆点表示用于建立SVM-HRM 预测模型的训练数据集，红色圆点表示用于测试SVM-HRM 预测模型的测试数据集。
图2 震级与震中距、记录数量的关系
2 模型特征参数
高铁地震预警只能采用单台震级估计的形式，对准确度的要求更高，因此在分析震级估计结果时需要融合多种参数。本文以获得最优预测结果为目标进行参数设置，选取幅值参数、周期参数、能量参数、衍生参数这4 大类共12 个特征参数作为支持向量机的高速铁路地震预警震级预测模型SVMHRM 的输入参数，其中，幅值参数、能量参数、衍生参数均进行震源距修正，统一校正到参考震源距10 km［12，20-21］。各特征参数的定义分别如下。
1）幅值参数
包括峰值位移Pd［11］、峰值速度Pv和峰值加速度Pa，计算式分别为
式中：0 为P 波的到达时刻；τ0为P 波到达后的时间窗长度；s（t）为竖向位移时程；v（t）为竖向速度时程；a（t）为竖向加速度时程。
2）周期参数
包括特征周期τc［5］、峰值比Tva［6］和构造参数PP［8］，计算式分别为
其中，
式中：r为依据帕塞瓦尔定理取的中间变量。
3）能量参数
6个月至2岁，如果最近没有给宝宝打过疫苗，孩子只是发烧，精神还好，没有其他症状，并且是第一次发烧，首先要怀疑幼儿急疹，到医院化验血常规和C反应蛋白。确认是病毒性的，这个月龄的孩子基本就是幼儿急疹了。
包括速度平方积分
累积绝对速度
和累计能量变化率
计算式分别为
式中：a3（t）为三分向合成的加速度。
4）衍生参数
包括竖向累积绝对位移cad、竖向累积绝对速度cav和竖向累积绝对加速度caa，计算式分别为
3 预测模型算法
支持向量机是人工智能领域中1 种基于统计学习的机器学习方法，可利用多个参数来进行模式分类以及非线性回归分析［21］。本文利用前述筛选得到的训练数据集，在获得其拟合函数的基础上，建立支持向量机模型算法，算法的关键是确立高斯核参数与训练参数。
3.1 支持向量机算法
定义参数：f（X）为预测震级；W 为各特征参数的权重向量；X 为特征参数组成的向量；b 为截距；n 为训练集数据数量，i 为其中第i 条数据；yi为与特征参数对应的震级；E 为容忍误差。利用线性回归函数f（X）=WT ∙X+b，对基于前述训练数据集计算得到的4大类共12个参数进行拟合。
假设经过拟合，所有的样本数据都可以在［-E，E］范围内通过线性函数f（X）表示［22］，其数学表达式为
记者看养老服务中心条件这么好，收费一定不低吧，张振美表示，企事业单位退休下来的老人完全能够负担得起，能自理老人收费标准在2000元左右，半自理老人收费在3000元左右，完全不能自理老人收费在4000元以上。
在此基础上，通过如下步骤建立支持向量机模型。
修正评价系统包括水资源、生态条件与环境质量两项要素，具体包含8项指标。区别于基础评价，修正评价各指标属于较宏观指标，一般不能利用小单元的土地利用图斑数据统计得到，而应该从区县等大范围数据汇总实现。该文根据《山东省土地资源承载状况评估研究》报告，获取荣成市各项修正评价指标。由于各项指标值之间单位不统一，不便于比较和分析，因此采用归一化法将各项修正评价指标进行归一化处理。归一化处理后，所有指标的取值范围为0～1，修正评价指标计算结果及归一化结果如表3所示。
1）问题转化
优化数据点与线性回归函数之间的距离，并将计算距离问题转化为极值问题，则有
式中：‖W‖为权重向量的模。
2）引入惩罚参数
不同于最小二乘拟合法，支持向量机算法允许一定范围内的拟合误差，所以可引入惩罚参数C，表示超出允许误差时对样本的惩罚程度，并将式（16）的目标函数（极值计算式）转化为
其中，
式中：
为松弛变量，也记为损失函数，即回归函数上侧和下侧的训练误差。
支持向量机的线性回归函数与损失函数之间的关系如图3所示。图中的实线表示线性回归函数
，虚线范围［-E，E］称为容忍误差范围，表示支持向量机的超平面。如果落在虚线范围内的样本点（空心圆）的误差可以忽略不计，那么落在这2 条虚线上的样本点（红色实心圆）则记为支持向量；落在虚线范围外的样本点（空心圆），即表示超过［-E，E］范围，就可记为损失函数S。
图3 支持向量机的线性回归函数与损失函数
3）目标函数转化
为了解决第2 步中的极值问题，将上述的目标函数转化为拉格朗日函数，从而解决有约束下的极值问题，由此，式（17）可转化为
（5）撰写技术指导原则有适用于各种抗菌药的下列描述，“在选择或调整抗菌药物治疗方案时，应进行细菌培养和药敏试验以分离并鉴定感染病原菌，确定其对该抗菌药物的敏感性。如果没有这些试验的数据做参考，则应根据当地细菌耐药性和抗菌药物敏感性等流行病学情况进行经验性治疗。在获得以上药敏结果之前可以先使用该抗菌药物进行治疗，获得药敏结果后再选择进行针对性的病原治疗。”中国说明书没有完整的这种描述，而英文说明书则有（见适应证和应用项目的应用部分）。
式中：α，β，γ为拉格朗日因子。
4）求拉格朗日函数
对式（18）求导，分别得到
急性化脓性中耳炎是耳鼻喉科的一种常见疾病，其发生频率较高，多是因细菌感染所导致的急性中耳粘膜化脓性病变，鼓室及其它部位均可受到累及，临床表现有耳痛、耳中流脓、鼓膜充血、耳道穿孔、听力下降等[1] 。急性化脓性中耳炎多继发于急性上呼吸道感染，各个年龄段均可发病，但以儿童其青年发病率为高[2] ，发病后需要立即给予有效的治疗。近年来笔者采用盐酸洛美沙星滴耳液治疗急性化脓性中耳炎患者，效果理想，现报告如下。
式中：αi，βi不同时为0，且αi＞0；K（Xi，X）为支持向量机的核函数。
通过以上4 步即可以得到支持向量机的线性回归函数为
3.2 支持向量机的高斯核参数
支持向量机与线性回归最重要的区别在于，前者通过核函数运算，将样本数据映射到高维空间，即：基于非线性函数的线性组合，其网络结构如图4所示。图中的样本点X 基于i 个支持向量进行非线性变换K（Xi，X），再基于支持向量机的线性回归函数即式（21），便可求得预测震级。
图4 支持向量机网络结构
核函数有多种形式，本文中核函数采用的是高斯核函数，该函数应用范围广、差值能力强，能够更好地提取样本数据的局部特征［27］，计算方法为
式中：λ为高斯核参数。
3.3 支持向量机的训练参数
合理的支持向量机模型需要选择合适的训练参数，以使预测结果与真实值的误差尽可能小。在支持向量机中与误差相关的参数有惩罚参数C、容忍误差E、高斯核参数λ，本文计算依据的是Cher⁃kassky和Ma［28-29］给出的经验计算，分别如下
式中：μ 为训练集数据输出结果的平均值；γ 为训练集数据输出结果的标准差；n 为训练集中的数据数量；η 为训练得到的预测值与真实值误差的标准差；q 为输入特征参数值的范围；m 为输入特征参数的数量。
根据式（23）—式（25）给出的训练参数计算方法，经过6 次交叉验证，确定不同时间窗下的惩罚参数C、容忍误差E 和高斯核参数λ，进而建立不同预测时间窗下基于支持向量机的高速铁路地震预警震级预测模型SVM-HRM。
4 SVM-HRM 预测模型的震级预测结果对比
基于测试数据集，对SVM-HRM 预测模型的震级预测结果进行如下对比。
（1）在3.0 s的时间窗下，将SVM-HRM 预测模型的震级预测结果与传统的震级预测模型τc方法与Pd方法进行对比。
（2）依据现行《高速铁路地震预警监测系统试验方法》相关条文，计算SVM-HRM 预测模型震级预测的单台实现率，并与《高速铁路地震预警监测系统试验方法》要求的震级预测实现率标准进行对比。
4.1 与传统震级预测模型的对比
定义SVM-HRM 预测模型的预测震级为预测值，定义地震事件的编目震级为真实值，如果取预测值与真实值之差为误差ωi，则ωi的标准差σ为
式中：N 为测试集中的数据数量，k 为其中第k 条数据。
在P 波到达后3.0 s 的时间窗下，将SVMHRM 预测模型与传统震级预测模型τc方法与Pd方法的震级预测结果进行对比，并按式（26）计算3种模型的标准差σ 值，结果如图5所示。图中的黑色实线表示预测值与真实值的1∶1 线性比例关系，红色虚线表示预测值与真实值误差的1 倍标准差±σ。
由计算并结合图5 可知：τc 方法、Pd 方法和SVM-HRM 预测模型的震级预测误差标准差分别是1.64，0.43 和0.30 个震级单位，即SVM-HRM预测模型远小于τc方法，也小于Pd方法。
《中国造纸学报》为大16开本。国内定价：纸质版30元/期，电子版30元/期，纸质版+电子版50元/期；国外及港澳台地区定价：纸质版30美元/期，电子版30美元/期，纸质版+电子版50美元/期。
图5 τc方法、Pd方法和SVM-HRM 的震级预测结果对比
需要注意的是，τc方法对5 级以下地震的震级预测结果存在明显高估，即存在“小震高估”现象，这是因为τc方法需要对地震数据进行震中距与信噪比筛选才能构建合适的预测模型。本文构建的SVM-HRM 预测模型没有这一步骤，大幅度提高了模型的普适性。同时，对比Pd方法的震级预测结果能清晰发现，SVM-HRM 预测模型的“小震高估”现象也得到了改善。
显然，对比τc方法与Pd方法，SVM-HRM 预测模型得到的震级预测准确性明显提升。
4.2 与高速铁路震级预测实现率标准的对比
SVM-HRM 预测模型基于训练数据集的震级预测误差标准差，以及基于测试数据集的震级预测误差标准差随预测时间窗的变化如图6所示。图中，横坐标为时间窗，0 s为P波的到达时刻。
由图6 可知：2 条曲线趋于重合、且在相同时间窗所对应误差标准差的最大差值不超过0.02，这表明在P 波到达后，所有预测时间窗所构建的SVM-HRM 预测模型均有极强的泛化能力，即SVM-HRM 预测模型在新鲜数据样本下的适应能力强，对具有同一规律的训练数据集以外的数据，经过训练的网络也能给出合适的输出；随着时间窗的增加，震级预测的误差标准差显著减小，这表明SVM-HRM 预测模型具备震级预测的连续性，且随着P波到达后时间窗的逐步增加，预测震级准确性显著增强。
图6 SVM-HRM 震级预测误差标准差随时间窗的变化
《高速铁路地震预警监测系统试验方法》对地震预警首报（第一次发布警报）预测震级实现率做出了规定与要求［15］，并以实现率作为衡量高速铁路地震预警系统震级预测准确性的指标参数。实现率根据预测震级误差绝对值小于等于1 的数量所占百分比来定义，计算式为
式中：r 为实现率；h 为测试数据中预测震级误差ωi≤1的数量；H为测试数据的总数量。
《高速铁路地震预警监测系统试验方法》中所规定的震级预测实现率标准见表1。本文在此研究实现率随预测时间窗的变化规律，旨在分析SVMHRM 预测模型在P 波到达后多久可以给出准确可靠的首报震级。
结合京港澳高速公路涿州(京冀界)—石家庄段改扩建工程，分析了地基处理前后，新老路堤表面工后沉降和横坡度的变化规律，研究了工后沉降作用下的路面结构层附加应力响应，计算了路面结构层的强度发挥率，并确定了应力控制层，制定了合理的差异沉降控制标准和经济有效的地基处理措施。主要结论如下：
近年来有学者提出，肠道菌群的改变可以影响NAFLD的疾病进展[60]。肠道微生物群破坏(生态失调)诱导的炎症可能有助于改变肠道的屏障功能和细菌移位。最近的证据表明，炎症小体是这一过程的核心，通过调节肠道微生物群从而调节由微生物群驱动的NAFLD表型。炎症小体是外源性PAMPs和DAMPs的传感器，其调节caspase-1介导对炎性细胞因子pro-IL-1β和pro-IL-18的前体的切割。NLRP3和NLRP6炎症小体敲除后，小鼠体内普雷沃菌科(Prevotellaceae)和细菌科TM7等主要成分显著增加，进而促进NAFLD疾病进展。
表1 现行震级预测实现率标准
震级预测模式单台多台震级范围/级3～8 3～5 5～7 7～8首报预测震级与实际震级偏差≤1的实现率不小于50%不小于30%不小于90%不小于60%
需要指出的是，在高速铁路地震预警中，多台的震级预测结果需要依据单台震级预测结果加权平均后得到，因此在多台实现率的对比分析中，只要单台分震级范围实现率统计结果优于多台实现率标准，就表明多台实现率结果优于标准，故只需要统计不同震级范围下的单台实现率与《高速铁路地震预警监测系统试验方法》中相应的震级范围实现率标准进行对比即可。
不同时间窗下，SVM-HRM 预测模型开展震级预测的单台实现率计算结果以及与《高速铁路地震预警监测系统试验方法》规定的实现率标准的对比如图7所示。图中的虚线分别为现行《高速铁路地震预警监测系统试验方法》规定的单台（红色虚线）与多台（蓝色虚线）实现率标准［15］。
图7 不同时间窗下SVM-HRM 预测模型的单台实现率与现行实现率标准对比
由图7（a）可知：对于震级范围在3～8 级之间的地震事件，SVM-HRM 预测模型得到单台预测震级实现率在P 波到达后的0.5 s 达到95%，即优于《高速铁路地震预警监测系统试验方法》规定的单台实现率标准；随着预测时间窗的增加，实现率逐渐增大，表明震级预测的准确性持续增加；当时间窗达到2.0 s 时，SVM-HRM 预测模型的预测震级实现率接近100%。
由图7（b）可知：对于震级范围在3～5 级之间的地震事件，SVM-HRM 预测模型在P 波到达后的0.5 s 实现率达到98%，即优于《高速铁路地震预警监测系统试验方法》规定的多台实现率标准；随着预测时间窗的增加，实现率逐渐增大，表明震级预测的准确性持续增加；当时间窗达到1.0 s时，SVM-HRM 预测模型的震级预测实现率接近100%。
对于进入军品市场的民营企业，定期对其法人资格、专业队伍、质量管理、财务资金状况、保密资格等要素实施连续性的承制资格例行监督。实行准入准出并行制度，对监督中发现的问题，应加大整改督促力度，充分利用警告、暂停资格等手段，促进民营企业整改。对于整体承制能力、质量管理水平下滑或出现重大问题的民营企业，依法取缔军品资格。
由图7（c）可知：对于震级范围在5～7 级之间的地震事件，SVM-HRM 预测模型在P 波到达后的1.5 s 实现率达到92%，即优于《高速铁路地震预警监测系统试验方法》规定的多台实现率标准；随着预测时间窗的增加，实现率逐渐增大，表明震级预测的准确性持续增加；当时间窗达到4.5 s时，SVM-HRM 预测模型的震级预测实现率接近100%。
由图7（d）可知：对于震级范围在7～8 级之间的地震事件，SVM-HRM 预测模型在P 波到达后的0.5 s 实现率为67%，即优于《高速铁路地震预警监测系统试验方法》规定的多台实现率标准；随着预测时间窗的增加，实现率逐渐增大，表明震级预测的准确性继续增加；当时间窗达到2.5 s时，SVM-HRM 预测模型的震级预测实现率达到95%并进入平台阶段；当时间窗达到6.0 s 时，实现率又继续增加，且在时间窗达到7.0 s 时，震级预测实现率达到100%。
通过上述分析可以发现，本文建立的支持向量机震级预测模型SVM-HRM 的震级预测准确性和连续性都得到极大提升，满足《高速铁路地震预警监测系统试验方法》中对震级预测的相关条文要求，可用于高速铁路地震预警系统震级预测。
5 结论
（1）在P 波到达后3.0 s 时间窗下，与传统的地震预警震级预测τc方法与Pd方法相比，SVMHRM 预测模型的震级预测误差明显减小，且小震高估的现象得到明显改善，准确性得到极大提升。
（2）SVM-HRM 预测模型的训练数据集可不做震中距与信噪比筛选，表明该模型具备极强的普适性；模型的训练数据集与测试数据集震级预测结果误差标准差趋于一致，且随着时间窗的增长而显著减小，表明该模型具备极强的泛化性能，震级预测连续性得到极大提升。
（3）将SVM-HRM 预测模型的震级预测结果对比《高速铁路地震预警监测系统试验方法》要求的首报震级预测实现率标准，可发现对于震级范围在3～8 级的地震事件，SVM-HRM 预测模型的单台震级预测实现率在P波到达后的0.5 s达到95%，且优于《高速铁路地震预警监测系统试验方法》要求的首报震级预测实现率标准；对于3～5，5～7，7～8这3 个震级范围下的地震事件，震级预测单台实现率分别在P 波到达后的0.5，1.5，0.5 s 优于多台实现率标准。这表明利用本文构建的SVMHRM 预测模型，最迟在P 波到达后的1.5 s可以发出满足《高速铁路地震预警监测系统试验方法》震级预测实现率标准的地震预警首报，准确性得到极大提升。
参考文献
［1］孙汉武，王澜，戴贤春，等.高速铁路地震紧急自动处置系统的研究［J］.中国铁道科学，2007，28（5）：121-127.
（SUN Hanwu，WANG Lan，DAI Xianchun，et al.Study on the Earthquake Urgent Automatic Alarm System of High Speed Railway［J］.China Railway Science，2007，28（5）：121-127.in Chinese）
［2］郭奇园.高速铁路地震监测预警信号接口单元的研发［J］.中国铁道科学，2016，37（3）：138-144.
（GUO Qiyuan.Research and Development of Signal Interface Unit in Earthquake Early Warning System for High Speed Railway［J］.China Railway Science，2016，37（3）：138-144.in Chinese）
［3］王澜，戴贤春，习年生，等.基于概率模型的高速铁路地震预警实时性与准确性分析［J］.中国铁道科学，2018，39（1）：131-137.
（WANG Lan，DAI Xianchun，XI Niansheng，et al.Real Time and Accuracy Analysis of Earthquake Early Warning for High Speed Railway Based on Probability Model［J］.China Railway Science，2018，39（1）：131-137.in Chinese）
［4］杨林.高速铁路地震预警系统紧急处置信息快速生成算法［J］.中国铁道科学，2018，39（3）：125-130.
（YANG Lin.Fast Generation Algorithm for Emergency Disposal Information in Earthquake Early Warning System for High Speed Railway［J］.China Railway Science，2018，39（3）：125-130.in Chinese）
［5］宋晋东.高速铁路运行控制用地震动参数及单台地震预警技术研究［D］.哈尔滨：中国地震局工程力学研究所，2013.
（SONG Jindong.Research on Seismic Ground Motion Indices for Operation Control and Single Station Earthquake Early Warning Applied for High-Speed Railway［D］.Harbin：Institute of Engineering Mechanics，China Earthquake Administration，2013.in Chinese）
［6］KANAMORI H.Real-Time Seismology and Earthquake Damage Mitigation［J］.Annual Review of Earth and Planetary Sciences，2005，33（1）：195-214.
［7］马强.地震预警技术研究及应用［D］.哈尔滨：中国地震局工程力学研究所，2008.
（MA Qiang.Study and Application on Earthquake Early Warning Technology［D］.Harbin：Institute of Engineering Mechanics，China Earthquake Administration，2008.in Chinese）
［8］HUANG P L，LIN T L，WU Y M.Application of τc*Pd in Earthquake Early Warning［J］.Geophysical Research Letters，2015，42（5）：1403-1410.
［9］ZIV A.New Frequency-Based Real-Time Magnitude Proxy for Earthquake Early Warning［J］.Geophysical Research Letters，2014，41（20）：7035-7040.
［10］李鸿杰，张建经.基于汶川主震及余震的一种新型地震预警震级估计参数τlog研究［J］.地球物理学进展，2018，33（1）：45-50.
（LI Hongjie，ZHANG Jianjing.Research on a New Magnitude Estimation Parameter τlog in Earthquake Early Warning System Based on the Wenchuan Sequences［J］.Progress in Geophysics，2018，33（1）：45-50.in Chinese）
［11］WU Y M，ZHAO L.Magnitude Estimation Using the First Three Seconds P-Wave Amplitude in Earthquake Early Warning［J］.Geophysical Research Letters，2006，33（16）：L16312.
［12］ZOLLO A，LANCIERI M，NIELSEN S.Earthquake Magnitude Estimation from Peak Amplitudes of Very Early Seismic Signals on Strong Motion Records［J］.Geophysical Research Letters，2006，33（23）：L23312.
［13］王子珺，赵伯明.高速铁路地震预警震级快速准确预测方法［J］.中国铁道科学，2017，38（2）：127-132.
（WANG Zijun，ZHAO Boming.Rapid and Accurate Prediction Method for Magnitude in Earthquake Early Warning for High Speed Railway［J］.China Railway Science，2017，38（2）：127-132.in Chinese）
［14］REDDY R，NAIR R R.The Efficacy of Support Vector Machines（SVM）in Robust Determination of Earthquake Early Warning Magnitudes in Central Japan［J］.Journal of Earth System Science，2013，122（5）：1423-1434.
［15］中国铁路总公司.Q/CR 634—2018高速铁路地震预警监测系统试验方法［S］.北京：中国铁道出版社，2018.
（China Railway Corporation.Q/CR 634—2018 Test Method for High-Speed Railway Earthquake Early-Warning and Monitoring System［S］.Beijing：China Railway Publishing House，2018.in Chinese）
［16］National Research Institute for Earth Science and Disaster Resilience（NIED）.Strong-Motion Seismograph Networks［DB/OL］.［2018-12-31］.http：//www.kyoshin.bosai.go.jp/.
［17］宋晋东，教聪聪，李山有，等.基于地震P 波双参数阈值的高速铁路Ⅰ级地震警报预测方法［J］.中国铁道科学，2018，39（1）：138-144.
（SONG Jindong，JIAO Congcong，LI Shanyou，et al.Prediction Method of First-Level Earthquake Warning for High Speed Railway Based on Two-Parameter Threshold of Seismic P-Wave［J］.China Railway Science，2018，39（1）：138-144.in Chinese）
［18］马强，金星，李山有，等.用于地震预警的P波震相到时自动拾取［J］.地球物理学报，2013，56（7）：2313-2321.
（MA Qiang，JIN Xing，LI Shanyou，et al.Automatic P-Arrival Detection for Earthquake Early Warning［J］.Chinese Journal of Geophysics，2013，56（7）：2313-2321.in Chinese）
［19］王子珺，赵伯明.高速铁路地震预警P 波与S 波复合自动快速识别的理论方法与应用［J］.中国铁道科学，2016，37（4）：121-127.
（WANG Zijun，ZHAO Boming.Theatrical Method and Application of Auto-Rapid P and S Waves Recognition in Earthquake Early Warning for High Speed Railway［J］.China Railway Science，2016，37（4）：121-127.in Chinese）
［20］彭朝勇，杨建思，薛兵，等.基于汶川主震及余震的预警参数与震级相关性研究［J］.地球物理学报，2013，56（10）：3404-3415.
（PENG Chaoyong，YANG Jiansi，XUE Bing，et al.Research on Correlation between Early-Warning Parameters and Magnitude for the Wenchuan Earthquake and Its Aftershocks［J］.Chinese Journal of Geophysics，2013，56（10）：3404-3415.in Chinese）
［21］YAMADA M，HEATON T，BECK J.Real-Time Estimation of Fault Rupture Extent Using Near-Source Versus Far-Source Classification［J］.Bulletin of the Seismological Society of America，2007，97（6）：1890-1910.
［22］FESTA G，ZOLLO A，LANCIERI M.Earthquake Magnitude Estimation from Early Radiated Energy［J］.Geophysical Research Letters，2008，35（22）：L22307.
［23］REED J W，KASSAWARA R P.A Criterion for Determining Exceedance of the Operating Basis Earthquake［J］.Nuclear Engineering and Design，1990，123（2/3）：387-396.
［24］NAKAMURA Y.A New Concept for the Earthquake Vulnerability Estimation and Its Application to the Early Warning System［M］.Berlin Heidelberg：Springer Berlin Heidelberg，2003：693-699.
［25］SAUNDERS C，STITSON M O，WESTON J，et al.Support Vector Machine［J］.Computer Science，2002，1（4）：1-28.
［26］李娇.支持向量机参数优化研究［D］.武汉：华中师范大学，2011.
（LI Jiao.Research of Parameter Optimization for Support Vector Machine［D］.Wuhan：Central China Normal University，2011.in Chinese）
［27］STEINWART I，HUSH D，SCOVEL C.Learning from Dependent Observations［J］.Journal of Multivariate Analysis，2009，100（1）：175-194.
［28］CHERKASSKY V，MA Y Q.Selection of Meta-Parameters for Support Vector Regression［C］// Artificial Neural Networks—ICANN，2002，2002（10）：687-693.
［29］CHERKASSKY V，MA Y Q.Practical Selection of SVM Parameters and Noise Estimation for SVM Regression［J］.Neural Networks，2004，17（1）：113-126.