非接触式城市洪涝积水监测异常数据修正方法技术

本发明专利技术的目的是提供一种非接触式城市洪涝积水监测异常数据修正方法，首先判断是积水监测设备数据传输故障造成的短时监测数据缺测，还是硬件故障导致无法监测；采用插补算法对短时积水监测数据缺测值进行数据插补；基于其他相近地区的积水监测设备值对因设备损坏导致长时间无法监测的积水点数据进行生成；然后对车辆经过或停车引起的异常波动数据进行修正；截取积水时段监测数据；最后对汽车行驶引起的振荡过程使用集合卡尔曼滤波算法进行修正。本发明专利技术解决了现有技术中存在的积水监测由于车辆或行人导致积水数据异常的问题。由于车辆或行人导致积水数据异常的问题。由于车辆或行人导致积水数据异常的问题。

【技术实现步骤摘要】
非接触式城市洪涝积水监测异常数据修正方法


[0001]本专利技术属于城市洪涝积水监测数据的质量
，具体涉及一种非接触式城市洪涝积水监测异常数据修正方法。

技术介绍

[0002]受全球气候变化导致的极端暴雨、城市“雨岛效应”形成的局部强降雨以及快速城市化引起的下垫面特性改变造成的影响，城市洪涝灾害频发，已成为制约城市发展与运行管理的一个非常突出的问题。目前新型城镇化快速推进，可以预见，城市洪涝灾害将在未来一段时期内长期存在。
[0003]积水监测城市洪涝防灾减灾体系最基本的组成部分之一，是洪涝预警预报、防洪预案编制及风险评估的重要基础。随着城市洪涝灾害形势越发严峻，许多城市已建立起了积水监测系统，实时监测区域内主要积水点的水深变化，并基于监测数据发布预警信息。随着监测数据的积累，逐步发现监测数据存在大量异常现象。对于非接触式积水监测设备，如超声波液位计，由于设备不与路面发生接触，安全性较高，但极易受到路面车辆行驶导致积水数据振荡，或车辆在液位计下停车和行人通过导致数据异常波动的问题，数据质量难以有效支撑城市洪涝灾害防治。

技术实现思路

[0004]本专利技术的目的是提供一种非接触式城市洪涝积水监测异常数据修正方法，解决了现有技术中存在的积水监测由于车辆或行人导致积水数据异常的问题。
[0005]本专利技术所采用的技术方案是，非接触式城市洪涝积水监测异常数据修正方法，具体按照以下步骤实施：
[0006]步骤1、判断是积水监测设备数据传输故障造成的短时监测数据缺测，还是硬件故障导致无法监测；
[0007]步骤2、采用插补算法对短时积水监测数据缺测值进行数据插补；
[0008]步骤3、基于其他相近地区的积水监测设备值对因设备损坏导致长时间无法监测的积水点数据进行生成；
[0009]步骤4、对车辆经过或停车引起的异常波动数据进行修正；
[0010]步骤5、截取积水时段监测数据；
[0011]步骤6、对汽车行驶引起的振荡过程使用集合卡尔曼滤波算法进行修正。
[0012]本专利技术的特点还在于，
[0013]步骤1具体按照以下步骤实施：
[0014]计算积水深度缺测数据的时间长短，计算公式为：
[0015]t
q
＝t
j
‑
t
i
[0016]其中，t
q
为缺测数据的时间间隔，t
j
为最新的缺测数据时间，t
i
为第一个缺测数据时间的前一时刻，当t
q
＜α时，α为阈值，则是传输故障造成的短时监测数据缺测；当t
q
≥α时，
则是硬件故障导致设备无法监测；i为第一个缺测数据的前一时刻数据序列编号；j为最新的缺测数据序列编号；q表示是缺测数据的时间间隔。
[0017]步骤2具体如下：
[0018][0019]式中，为缺测数据；i为缺测数据的时间顺序，i＝1,2,3
……
；d
n
为缺测前最后一位实测数据，d
m
为缺测后第一位实测数据，Δt为监测数据的时间间隔，Δt
m
‑
n
为监测数据d
m
和d
n
的时间间隔；n为缺测前最后一位实测数据的数据序列编号；m为缺测后第一位实测数据的数据序列编号；q表示当前计算的水位数据为缺测数据。
[0020]步骤3具体按照以下步骤实施：
[0021]以故障设备为圆心，计算各邻近设备与故障设备的距离，进一步获得最短距离，确定距离故障设备最近的1个正常设备；
[0022]计算各正常运行设备与已故障设备的距离l
b
：
[0023][0024]l
min
＝min(l1,l2,l3,...,l
b
)
[0025]式中，l
b
为正常设备与故障设备的距离；b为设备编号，b＝1、2、3
…
；x
f
、y
f
分别为故障设备在x和y方向相对坐标；x
b
、y
b
分别为各正常设备在x和y方向相对坐标；l
min
为最小距离；f表示坐标为故障设备的坐标值。
[0026]计算最近一次积水事件中两个设备的监测数据倍比关系：
[0027][0028][0029]式中，m
z
为故障设备数据与正常设备数据的比值，z为积水数据的时间顺序，z＝1,2,3
……
；m
average
为所有积水数据比值的平均比值；为正常设备的监测值；为已故障设备的监测值；
[0030]根据比值与正常设备的监测数据生成故障设备的监测数据：
[0031][0032]式中，为缺测数据；i为缺测数据的时间顺序，i＝1,2,3
……
；为正常设备的监测数据。
[0033]步骤4具体按照以下步骤实施：
[0034]首先对车辆经过引起的异常波动进行识别与处理：
[0035]对于i时刻水深d
i
≠0，且i时刻前一时刻的水深d
i
‑1，i时刻后一时刻的水深d
i+1
均等于0的情况，令d
i
＝0；
[0036]其次对车辆停车引起的异常波动数据进行识别与处理：
[0037]对于i时刻水深d
i
≥γ，γ为事先设定的水深阈值，且i时刻前一时刻的水深d
i
‑1＝0，i时刻后一时刻的水深d
i+1
＝d
i
，且令d
i
＝d
i+1
＝0；若d
c
，c＝i+2,i+3,
…
都等于未修正前的
d
i
，则令d
c
＝0，直到监测值为0时，结束该项循环修正。
[0038]步骤5具体按照以下步骤实施：
[0039]确定积水时段的开始时刻和结束时刻，非零值前五分钟无积水，下一时刻有积水，则非零值前一时刻视为这一积水的开始时刻t
star
；最后一个非零值后连续五分钟无积水，非零值后一时刻视为这一积水时段的结束时刻t
end
，
[0040]确定了时段t
star
～t
end
后，从数据存储数据库筛选出该时段内的积水监测数据值d
i
，i＝star,stari+1,
…
,end。
[0041]步骤6具体按照以下步骤实施：
[0042]使用集合卡尔曼滤波算法EnKF修正因汽车行驶引起的异常震荡数据；
[0043]滤波算法具体流程如下：
[0044]步骤6.1、初始化背景场，给出N个符合高斯分布的随机变本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.非接触式城市洪涝积水监测异常数据修正方法，其特征在于，具体按照以下步骤实施：步骤1、判断是积水监测设备数据传输故障造成的短时监测数据缺测，还是硬件故障导致无法监测；步骤2、采用插补算法对短时积水监测数据缺测值进行数据插补；步骤3、基于其他相近地区的积水监测设备值对因设备损坏导致长时间无法监测的积水点数据进行生成；步骤4、对车辆经过或停车引起的异常波动数据进行修正；步骤5、截取积水时段监测数据；步骤6、对汽车行驶引起的振荡过程使用集合卡尔曼滤波算法进行修正。2.根据权利要求1所述的非接触式城市洪涝积水监测异常数据修正方法，其特征在于，所述步骤1具体按照以下步骤实施：计算积水深度缺测数据的时间长短，计算公式为：t
q
＝t
j
‑
t
i
其中，t
q
为缺测数据的时间间隔，t
j
为最新的缺测数据时间，t
i
为第一个缺测数据时间的前一时刻，当t
q
＜α时，α为阈值，则是传输故障造成的短时监测数据缺测；当t
q
≥α时，则是硬件故障导致设备无法监测；i为第一个缺测数据的前一时刻数据序列编号；j为最新的缺测数据序列编号；q表示是缺测数据的时间间隔。3.根据权利要求2所述的非接触式城市洪涝积水监测异常数据修正方法，其特征在于，所述步骤2具体如下：式中，为缺测数据；i为缺测数据的时间顺序，i＝1,2,3
……
；d
n
为缺测前最后一位实测数据，d
m
为缺测后第一位实测数据，Δt为监测数据的时间间隔，Δt
m
‑
n
为监测数据d
m
和d
n
的时间间隔；n为缺测前最后一位实测数据的数据序列编号；m为缺测后第一位实测数据的数据序列编号；q表示当前计算的水位数据为缺测数据。4.根据权利要求3所述的非接触式城市洪涝积水监测异常数据修正方法，其特征在于，所述步骤3具体按照以下步骤实施：以故障设备为圆心，计算各邻近设备与故障设备的距离，进一步获得最短距离，确定距离故障设备最近的1个正常设备；计算各正常运行设备与已故障设备的距离l
b
：l
min
＝min(l1,l2,l3,...,l
b
)式中，l
b
为正常设备与故障设备的距离；b为设备编号，b＝1、2、3
…
；x
f
、y
f
分别为故障设备在x和y方向相对坐标；x
b
、y
b
分别为各正常设备在x和y方向相对坐标；l
min
为最小距离；f表示坐标为故障设备的坐标值；计算最近一次积水事件中两个设备的监测数据倍比关系：
式中，m
z
为故障设备数据与正常设备数据的比值，z为积水数据的时间顺序，z＝1,2,3
……
；m
average
为所有积水数据比值的平均比值；为正常设备的监测值；为已故障设备的监测值；根据比值与正常设备的监测数据生成故障设备的监测数据：式中，为缺测数据；i为缺测数据的时间顺序，i＝1,2,3
……
；为正常设备的监测数据。5.根据权利要求4所述的非接触式城市洪涝积水监测异常数据修正方法，其特征在于，所述步骤4具体按照以下步骤实施：首先对车辆经过引起的异常波动进行识别与处理：对于i时刻水深d
i
≠0，且i时刻前一时刻的水深d
i
‑1，i时刻后一时刻的水深d
i+1

【专利技术属性】
技术研发人员：陈光照，侯精明，潘鑫鑫，吕佳豪，景静，李轩，梁鑫，王兴桦，范臣臣，高徐军，刘园，薛树红，吴豪，毛叶丰，马红丽，
申请(专利权)人：中国电建集团西北勘测设计研究院有限公司中智水科宁波科技有限公司，
类型：发明
国别省市：