高速公路雾灾风险等级预报技术实现方法
【专利摘要】本发明涉及一种高速公路雾灾风险等级预报技术实现方法。本发明是按照气象灾害理论的概念，从致灾因子危险性、承灾体暴露性等方面入手，将高速公路雾灾风险指标量化、按照等级分别赋值，易计算、易理解，达到可预报性，实现了逐日、甚至是更短时效内计算和发布高速公路沿线相邻两个收费站未来一定时段内雾灾风险等级预报，改变了以往气象预报预警服务产品仅仅基于气象条件并且无法细分到局部路段的局面。本发明不仅基于气象因素，还综合考虑到了公路路网密度、车流量等因素，所建立的高速公路雾灾风险计算模型，可计算出高速公路雾灾风险等级预报结果，实现对高速公路雾灾综合风险的预报和预警。
【专利说明】
高速公路雾灾风险等级预报技术实现方法
技术领域
[0001] 本发明设及一种应用气象技术，具体地说是一种高速公路雾灾风险等级预报技术 实现方法。
【背景技术】
[0002] 目前，我国高速公路的通车里程为世界之首。统计发现，高速公路通行受阻（即采 取封闭管制措施)事件有38.3%是因气象因素(包括雾、路面结冰、降雨、大风等)造成的。其 中W雾造成的封闭事件局达80.2%，同时雾造成的通行雙阻时间也最长。可见，雾是对局速 公路通行影响最为严重的灾害性天气。
[0003] 然而，W往的灾害性天气预报预警几乎完全取决于气象条件，依照气象行业的标 准和指标，而并没有考虑到公路交通方面的综合因素，也没有实现基于影响和风险的预报 预警服务。
【发明内容】
[0004] 本发明的目的就是提供一种高速公路雾灾风险等级预报技术实现方法，W解决高 速公路交通服务中缺少基于影响和风险的预报预警服务的问题。
[0005] 本发明是运样实现的：一种高速公路雾灾风险等级预报技术实现方法，包括W下 步骤：
[0006] ①、巧吿雾巧风盼预根横巧为：
[0007]
[000引其中，R为雾灾风险等级指数，P代表雾的强度分指标(运里用能见度最小值Vmin表 示），T代表雾的持续时间分指标，D代表路网密度分指标，N代表单项车流量分指标;P、t、d、n 分别为各分指标对应的权重系数；
[0009] ②、按下式确定雾的强度分指标P:
[0010]
[00川其中，Vmin为能见度最小值；
[0012] ③、确定雾的持续时间分指标Tj:为高速公路每个收费站匹配雾的持续时间预报， W小时为单位，能见度《50m的持续时间记为Tl, 50m<能见度《200m的持续时间记为T2,200m <能见度《500m的持续时间记为T3;tj分别表示它们的权重系数；
[0013] ④、按下式确定路网密度分指标D:
[0014]
[001引其中，P为路网密度值，单位为km/km2;
[0016]⑤、按下式确定单项车流量分指标N:
[0017]
[001引其中，1为高速公路某收费站预报时段内单项车流量(含上、下站口和过路车辆）， 单位为万辆；
[0019] ⑥、各分指标按下式进行无量纲化处理后参与计算：
[0020]
[0021] ⑦、雾灾风险预报等级指数计算模型的权重系数的确定：
[0022] 对因子进行两两比较建立成对比较矩阵，即每次取两个因子Xi和Xj ,Wau表示Xi和 Xj对某个指标Z的影响大小之比，引用数字1-9及其倒数作为影响程度的标度：
[0023]
[0024] 建立如下矩阵：
[0025：
[0026] 计算矩阵A的最大特征值和特征向量，通过一致性检验，作为最终的权重系数；
[0027] ⑧、雾灾风险预报等级划分：
[0028] 通过雾灾风险预报模型计算出R的数值序列，采用标准差分类法，W雾灾风险预报 指数的平均值n和标准差m为依据，将雾灾风险预报指数分为W下四个等级：
[0029] 当3<]1-1]1时，高速公路雾灾风险预报等级为1级；
[0030] 当n-m《R<^4，高速公路雾灾风险预报等级为2级；
[0031] 当n《R<n+m时，高速公路雾灾风险预报等级为3级；
[0032] 当3>11+111时，高速公路雾灾风险预报等级为4级。
[0033] 本发明还包括W下步骤：
[0034] ⑨、根据每条高速公路的各个收费站计算得出的雾灾风险预报指数等级，绘制实 时的雾灾风险等级预报图。
[0035] 本发明中的雾灾风险预报模型是综合气象因素、公路交通因素，从致灾因子危险 性和承灾体的空间脆弱性角度考虑，W雾的强度、持续时间（由于雾灾风险等级预报是针对 高速公路沿线某个收费站或气象站具体点的，因此不需要将雾的空间影响范围考虑到致灾 因子危险性指标)作为致灾因子危险性指标，W公路路网密度、逐日车流量作为承载体空间 脆弱性指标，因子建立高速公路雾灾风险预报模型，计算出雾灾风险预报指数序列后，再进 行雾灾风险预报等级划分。
[0036] 本发明是按照气象灾害理论的概念，从致灾因子危险性、承灾体暴露性等方面入 手，将高速公路雾灾风险指标量化、按照等级分别赋值，易计算、易理解，达到可预报性，实 现了逐日、甚至是更短时效内计算和发布高速公路沿线相邻两个收费站未来一定时段内雾 灾风险等级预报，改变了 W往气象预报预警服务产品仅仅基于气象条件并且无法细分到局 部路段的局面，解决了高速公路交通服务中缺少基于影响和风险的预报预警服务的问题。
[0037] 本发明高速公路雾灾风险等级预报技术实现方法不仅基于气象因素，还综合考虑 到了公路路网密度、车流量等因素，所建立的高速公路雾灾风险计算模型，根据模型计算出 高速公路雾灾风险等级预报结果，实现了对高速公路雾灾综合风险的预报和预警。
【附图说明】
[0038] 图1是河北省各高速公路收费站的某一时刻雾灾风险等级图表。
[0039] 图2是河北省各高速公路某一时刻雾灾风险等级预报示意图。
【具体实施方式】
[0040] 本发明W河北省2013年12月8日的大范围降雾天为例，进行雾灾风险等级预报，步 骤如下：
[0041 ]第一步，确定雾灾风险预报模型：
[0042]
[0043] 其中，R为雾灾风险等级指数，P代表雾的强度分指标(运里用能见度最小值Vmin表 示），T代表雾的持续时间分指标，D代表路网密度分指标，N代表单项车流量分指标;P、t、d、n 分别为各分指标对应的权重系数。
[0044] 第二步，确定雾的强度分指标P:通过前一日未来24小时内能见度预报，为高速公 路每个收费站匹配能见度最小值预报，能见度最小值Vmin巧Om的赋值为5，50m《 Vmin<200m的 赋值为4，200m《 vmin<500m的赋值为3，500m《 VminCOOOm的赋值为2，Vmin^ 1000m的赋值为1。
[0045] 第=步，确定雾的持续时间分指标L :为高速公路每个收费站匹配雾的持续时间 预报，W小时为单位;能见度《50m的持续时间记为Tl,50m<能见度《200m的持续时间记为 T2,200m<能见度《500m的持续时间记为T3。
[0046] 第四步，确定路网密度分指标D:利用Arc gis处理好高速公路各收费站处的路网 密度数值(一定时段内是固定值，单位为km/km2)，对各收费站的路网密度分指标进行赋值， 路网密度<4.33的赋值为3,4.33《路网密度<11.34的赋值为2,路网密度>11.34的赋值为 Io
[0047] 第五步，确定单项车流量分指标N:对高速公路各收费站一天内单项车流量(含上、 下站口和过路车辆)进行统计分析，得出规律，便于做车流量预报;利用各收费站八日单向 车流量(含上、下站口和过路车辆，取单向最大值)对车流量分指标进行赋值，车流量在2万 辆W上的赋值为3，车流量在0.5-2万辆之间的赋值为2，车流量小于0.5万辆的赋值为1。
[004引 笛六击.刑田W下公立；
[0049]
[0050] 对上述四个分指标(P、Tj、D、N)进行无量纲化处理，maxi《i《n{xi}对于各分指标取最 大等级对应的数值或取最大值，即雾的强度分指标P取5,雾的持续时间分指标L取24h，路 网密度分指标D取3，单向车流量分指标N取3。
[0051 ]第屯步，确定权重系数：
[0052] 1、确定雾的持续时间分指标L中的权重系数ti，根据ti对该分指标的贡献，建立如
下矩r
[0053:
[0化4] 计算得出该矩阵的最大特征值为3.0649;特征向量为0.6491，0.2790，0.0719; - 致性指标为0.0324，平均随机一致性指标取0.58，计算得出一致性比例为0.0559<0.1，因此 =个特征向量可W作为ti的权重系数。
[0055] 2、雾灾风险预报模型的权重系数p、t、d、n，针对四个分指标对雾灾风险等级的贡 献，建立如下矩阵：
[0化6]
[0057]计算得出该矩阵的最大特征值为4.2272;特征向量为0.5255,0.3128,0.1000， 0.0617; -致性指标为0.0757，平均随机一致性指标取0.9，计算得出一致性比例为0.0848< 0.1，因此四个特征向量可W作为模型四个分指标的权重系数。
[005引第八步，根据第一步确定的雾灾风险预报模型，计算出河北省各高速公路收费站 的雾灾风险等级，并绘制成图表(见图1)，河北省全省37条高速公路沿线345个收费站的雾 灾风险等级指数在0-0.7之间。
[0059]第九步，划分雾灾风险:利用标准差分类法，W雾灾风险预报指数的平均值n和标 准差111为依据，当1?<]1-1]1，]1-1]1《1?<]1，]1《1?<]1+1]1，1?>]1+1]1时，高速公路雾灾风险预报等级分 别对应I级、2级、3级、4级；
[0060]
[DD61」 呆下巧，化化母巧同化公蹄合/ r ^1义巧化TT异巧m的务火W陶TJM化巧双寺级，5会制 雾灾风险等级预报图（见图2)。
【主权项】
1. 一种高速公路雾灾风险等级预报技术实现方法，其特征是，包括以下步骤： ① 、设定雾灾风险预报模型为：其中，R为雾灾风险等级指数，P代表雾的强度分指标(这里用能见度最小值Vmin表示），τ 代表雾的持续时间分指标，D代表路网密度分指标，Ν代表单项车流量分指标;p、t、d、n分别 为各分指标对应的权重系数； ② 、按下式确定雾的强度分指标P:其中，Vmin为能见度最小值； ③ 、确定雾的持续时间分指标为高速公路每个收费站匹配雾的持续时间预报，以小 时为单位，能见度彡50m的持续时间记为T!，50m〈能见度彡200m的持续时间记为T 2，200m〈能 见度<500m的持续时间记为T3;tj分别表示它们的权重系数； ④ 、按下式确定路网密度分指标D:其中，P为路网密度值，单位为km/km2; ⑤ 、按下式确定单项车流量分指标N:其中，1为高速公路某收费站预报时段内单项车流量(含上、下站口和过路车辆），单位 为万辆； ⑥ 、各分指标按下式进行无量纲化处理后参与计算： ⑦ 、雾灾风险预报等级指数计IfT关乂里 亦:mmjfl用疋： 对因子进行两两比较建立成对比较矩阵，即每次取两个因子，以表示xdPxWi 某个指标Z的影响大小之比，引用数字1-9及其倒数作为影响程度的标度：建立如下矩阵：计算矩阵A的最大特征值和特征向量，通过一致性检验，作为最终的权重系数； ⑧ 、雾灾风险预报等级划分： 通过雾灾风险预报模型计算出R的数值序列，采用标准差分类法，以雾灾风险预报指数 的平均值η和标准差m为依据，将雾灾风险预报指数分为以下四个等级： 当R<n-m时，高速公路雾灾风险预报等级为1级； 当n-m彡R<n时，高速公路雾灾风险预报等级为2级； 当η彡R<n+m时，高速公路雾灾风险预报等级为3级； 当R>n+m时，高速公路雾灾风险预报等级为4级。2.根据权利要求1所述的高速公路雾灾风险等级预报技术实现方法，其特征是，还包括 以下步骤： ⑨ 、根据每条高速公路计算得出的雾灾风险预报指数等级，绘制实时的雾灾风险等级 预报图。
【文档编号】G06Q50/26GK105956986SQ201610541921
【公开日】2016年9月21日
【申请日】2016年7月11日
【发明人】曲晓黎, 张娣, 郭蕊, 赵增保, 武辉芹, 范俊红, 张成伟
【申请人】河北省气象服务中心