CN111400911A - 一种基于ewma控制图的gnss变形信息识别与预警方法 - Google Patents

一种基于ewma控制图的gnss变形信息识别与预警方法 Download PDF

Info

Publication number
CN111400911A
CN111400911A CN202010183171.8A CN202010183171A CN111400911A CN 111400911 A CN111400911 A CN 111400911A CN 202010183171 A CN202010183171 A CN 202010183171A CN 111400911 A CN111400911 A CN 111400911A
Authority
CN
China
Prior art keywords
ewma
control chart
gnss
early warning
data
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010183171.8A
Other languages
English (en)
Inventor
代阳
吴昊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Anhui University of Science and Technology
Original Assignee
Anhui University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Anhui University of Science and Technology filed Critical Anhui University of Science and Technology
Priority to CN202010183171.8A priority Critical patent/CN111400911A/zh
Publication of CN111400911A publication Critical patent/CN111400911A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
    • G01B7/16—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring the deformation in a solid, e.g. by resistance strain gauge
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02—Alarms for ensuring the safety of persons
    • G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明公开了一种基于EWMA控制图的GNSS变形信息识别与预警，包括以下步骤：(1)获取监测体的GNSS变形监测数据，并计算其计算均值μ和标准差σ；(2)根据监测数据的特点，设置EWMA控制图的权重参数λ和控制限限宽L；(3)将含有变形数据的监测数据，构建EWMA控制图的偏移统计量Z；(4)确定异常值识别模型的中心线和上下参考限；(5)对比判断EWMA偏移统计量是否超过控制限，对于超过控制限的数据定义为异常数据并进行预警。本发明首次将EWMA控制图运用于GNSS变形信息的识别与预警，能够准确的分析监测数据的变化趋势和变形信息，具有较好的变形信息检验能力和误报率较低的优点。

Description

一种基于EWMA控制图的GNSS变形信息识别与预警方法
技术领域
本发明涉及测绘变形监测的预警领域，具体是一种基于EWMA控制图的GNSS变形信息识别与预警方法。
背景技术
随着全球经济的快速发展，人类的活动范围不断扩张，高层建筑、大型桥梁和高速公路不断地建造，由建筑物变形所带来的灾害事件层出不穷。为了保障人民的生命财产安全，有必要对建筑物进行有效地变形监测，而GNSS定位技术具有实时动态高精度的优势，被广泛地应用于变形监测领域。将GNSS监测数据，采用一定的数学方法对所获取的变形信息进行数据处理与分析，及时进行预警，以期降低灾害发生的概率并减少其影响范围。
目前在对GNSS监测数据进行检验和预警常用的方法为CUSUM控制图方法，但其存在两个问题：一是对连续大偏移变形信息的识别结果较差，二是随着变形量的不断增加其误报率也随之增大。鉴于此，本文提出采用EWMA控制图进行GNSS灾害信息的识别与预警，并结合GNSS坐标序列的相关特征，提出一种新的GNSS变形信息识别与预警方法。该方法对于变形信息的识别精度高于经典的CUSUM控制图方法，且误警率得到了较大改善。
发明内容
针对现有技术存在的问题，本发明所要解决的是提供一种基于EWMA控制图的GNSS变形信息识别与预警方法。本发明首次将EWMA控制图用于GNSS变形信息的识别与预警，能够准确的分析监测数据的变化趋势和变形信息，具有较好的变形信息检验能力和误报率较低的优点。
本发明实现发明目的采用如下技术方案：
一种基于EWMA控制图的GNSS变形信息识别与预警方法，包括以下步骤：
步骤1，获取监测体的GNSS变形监测数据，并计算其计算均值μ和标准差σ；
步骤2，设置EWMA控制图的权重参数λ和控制限限宽L；
步骤3，将获取到的监测数据，构建EWMA控制图的偏移统计量Z；
步骤4，确定异常值识别模型的中心线和上下参考限；
步骤5，对比判断EWMA偏移统计量是否超过控制限，对于超过控制限的数据定义为异常数据并进行预警。
作为优选，根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤1中，采集对监测体(构筑物或建筑物)进行形变监测的 GNSS时间序列，对时间序列进行初步分析，计算均值μ和标准差σ。
作为优选，根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤2中，EMWA控制图的权重参数λ一般取值范围为[0.05, 0.25]，为突出控制图检验性能，λ取值为0.1，而控制限限宽取值为3。
作为优选，根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤3中，构建EWMA控制图的偏移统计量Z。
对于长时间观测的GNSS监测数据X(i),i＝1,2,…,n，n为样本大小。EWMA控制图的偏移统计量Zi定义为：
Zi＝λXi+(1-λ)Zi-1 (1)
式中λ是EWMA控制图的权值参数，是属于[0,1]中的常数，λ取值为0.1；初始值Z0(i＝1处)是过程目标值，用观测值X(i)的均值取代，Z0＝μ0。
因为Zi是EWMA所有先前样本均值的加权平均值，通过替换方程右侧的Zi-1来获得：
Zi＝λXi+(1-λ)[λXi-1+(1-λ)Zi-2]＝λXi+λ(1-λ)Xi-1+(1-λ)2Zi-2 (2)
继续以递归代替Zi-j，j＝2,3,…,t，我们得到
Figure DEST_PATH_1
作为优选，根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤4中，确定异常值识别模型的中心线CL和上参考限UCL 与下参考线LCL。
GNSS观测数据X(t)是具有方差σ2的独立随机变量，则构建的EWMA统计量Zi的方差如下所示：
Figure BDA0002413266670000022
因此，通过将GNSS监测数据构建为EWMA控制图的统计量Zi，与时间(历元) 共同组成EWMA控制图，其变形信息的检验准则如下：
Figure BDA0002413266670000031
CL＝μ0 (6)
Figure BDA0002413266670000032
式中L为构建检验阈值时控制限的宽度，一般取值为3；异常值识别模型的上参考限UCL、中心线CL和与下参考线LCL计算公式如式(5)，式(6)与式(7)所示；当构建的偏移统计量超过UCL或-LCL时，视为变形信息，并发出预警。
作为优选，根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤5中，对比判断由步骤3得到的EWMA偏移统计量Z是否超过上参考限UCL与下参考线LCL构成的控制限范围，对于超过控制限的数据定义为异常数据并进行预警
与现有技术相比，本发明的有益效果是：本发明首次将EWMA控制图用于GNSS 变形信息的识别与预警，能够准确的分析监测数据的变化趋势和变形信息，具有较好的变形信息检验能力和误报率较低的优点。
附图说明
图1是本发明一种基于EWMA控制图的GNSS变形信息识别与预警方法的较佳实施例的流程图；
图2是本发明案例中的实验采用的无变形信息的原始GNSS坐标序列；
图3是本发明案例中的分别加入1-4倍标准差变形信息的数据；
图4是本发明案例中的EWMA控制图对不同标准差变形数据的识别结果；
图5是本发明案例中的CUSUM控制图对不同标准差变形数据的识别结果；
具体实施方式
下面结合附图对本发明的较佳实施例进行详细阐述，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。
请参阅图1，本发明实施例包括：
一种基于EWMA控制图的GNSS变形信息识别与预警方法，包括以下步骤：
步骤1，获取监测体的GNSS变形监测数据，并计算其计算均值μ和标准差σ；
步骤2，设置EWMA控制图的权重参数λ和控制限限宽L；
步骤3，将获取到的监测数据，构建EWMA控制图的偏移统计量Z；
步骤4，确定异常值识别模型的中心线和上下参考限；
步骤5，对比判断EWMA偏移统计量是否超过控制限，对于超过控制限的数据定义为异常数据并进行预警。
以某段建筑物的GNSS监测时间序列为例，具体实施步骤如下：
第一步，采用一组含有模拟变形信息的GNSS监测数据进行试验。
原始的GNSS数据来源于安徽理工大学的某楼顶的GNSS连续跟踪站，站点相对稳固，并未有变形迹象。截取站点中一段时间的X方向坐标序列(采样时间是1500s，采样间隔是1s，采样频率1Hz)，并减去其真实值(由长时间观测平均得到)的数据作为原始GNSS监测数据，如图2所示。
在原始GNSS坐标序列的701-900历元处分别加入1-4倍的标准差的上偏移，形成含有四组突变型变形的GNSS监测数据，如图3所示。并对含有变形信息的GNSS监测序列进行分析，计算其计算均值μ和标准差σ。
第二步，设置EWMA控制图的权重参数λ和控制限限宽L。EMWA控制图的权重参数λ一般取值范围为[0.05,0.25]，为突出控制图检验性能，λ取值为0.1，而控制限限宽取值为3。
第三步，根据含有变形信息的GNSS监测数据，构建EWMA控制图的偏移统计量 Z，公式如下：
对于长时间观测的GNSS监测数据X(i),i＝1,2,…,n，n为样本大小。EWMA控制图的偏移统计量Zi定义为：
Zi＝λXi+(1-λ)Zi-1 (1)
式中λ是EWMA控制图的权值参数，是属于[0,1]中的常数，λ取值为0.1；初始值Z0(i＝1处)是过程目标值，用观测值X(i)的均值取代，Z0＝μ0。
因为Zi是EWMA所有先前样本均值的加权平均值，通过替换方程右侧的Zi-1来获得：
Zi＝λXi+(1-λ)[λXi-1+(1-λ)Zi-2]＝λXi+λ(1-λ)Xi-1+(1-λ)2Zi-2 (2)
继续以递归代替Zi-j，j＝2,3,…,t，我们得到
Figure 456021DEST_PATH_1
第四步，确定基于EWMA控制图的异常值识别模型的中心线CL和上参考限UCL 与下参考线LCL。
GNSS观测数据X(t)是具有方差σ2的独立随机变量，则构建的EWMA统计量Zi的方差如下所示：
Figure BDA0002413266670000051
因此，通过将GNSS监测数据构建为EWMA控制图的统计量Zi，与时间(历元) 共同组成EWMA控制图，其变形信息的检验准则如下：
Figure BDA0002413266670000052
CL＝μ0 (6)
Figure BDA0002413266670000053
式中L为构建检验阈值时控制限的宽度，一般取值为3；异常值识别模型的上参考限UCL中心线CL和与下参考线LCL计算公式如式(5)，式(6)与式(7)所示；当构建的偏移统计量超过UCL或-LCL时，视为变形信息，并发出预警。
第五步，基于上述步骤，采用EWMA控制图对上述加入1-4倍标准差的四组含有变形信息的GNSS监测数据，进行异常值检验分析，如图4所示。采用通常方法CUSUM 控制图对上述加入1-4倍标准差的四组含有变形信息的GNSS监测数据，进行异常值检验分析，如图5所示。将两种方法得到的结果进行分析，如表1与表2所示。
表1不同标准差变形数据的CUSUM控制图识别结果
Figure BDA0002413266670000061
表2不同标准差变形数据的EWMA控制图识别结果
Figure BDA0002413266670000062
Figure BDA0002413266670000071
由图4、图5与表1、表2所知：对于在701-900历元加入1-4倍标准差上偏移的变形信息中，CUSUM控制图的上偏移的预警数量分别为327、471、492和518，而下偏移的预警数量为324、408、427和438；EWMA控制图的上偏移的预警数量分别为 156、213、215和267，下偏移的预警数量分别为117、124、123和367。总的来说，无论是加入到上偏移变形信息，还是实测数据中的下偏移变形数据，EWMA控制图检验与预警的准确度均高于CUSUM控制图，其结果具有一定的可靠度。在对2倍以下的连续小偏移变形信息进行检验与预警时，CUSUM控制图检验与预警的结果一般，并且在2倍时上偏移变形信息的误警数量已经超过了实际所加入变形信息的真实范围，因此其检验与预警结果的准确度较低；与CUSUM控制图相比，EWMA控制图所检验出的变形信息多围绕在所加入的真实变形信息范围内，其结果优于CUSUM控制图。此外，本研究在测试数据中并没有加入任何下偏移的变形信息，但两种控制图算法均检验到了一定的变形数据，这可能由实测数据中的粗差所引起。其中，随着标准差变形信息的不断增大，EWMA控制图与CUSUM控制图的误警率均不断增加，从表1-2和图4- 5中可以看出，EWMA控制图在误警率方面得到了较大改善，其误警率下降了50％以上。同时，针对下偏移变形信息的检验与预警，与CUSUM控制图相比，EWMA控制图对粗差有一定的抑制作用，使其检验预警的结果具有一定的可靠性和稳定性。
通过对结果分析，可得到以下结论：
(1)针对不同标准差变形信息的检验与预警，EWMA控制图与经典的CUSUM控制图均能有效地检验出变形信息。
(2)在对2倍标准差以上的连续大偏移变形信息的检验与预警时，EWMA控制图的检测预警能力明显优于CUSUM控制图，同时误警率也得到了较大改善
(3)虽然在测试数据中没有加入任何下偏移的变形信息，但两种控制图算法均检验到了一定的变形数据，这可能由实测数据中的粗差所引起。但随着标准差的不断增大，检验出的变形信息范围也在不断增加，与CUSUM算法相比，EWMA控制图对粗差有一定的抑制效果，其检验与预警的结果具有一定的可靠性和稳定性。
由本发明所要解决的是提供一种基于EWMA控制图的GNSS变形信息识别与预警方法。本发明首次将EWMA控制图运用于GNSS变形信息的识别与预警，能够准确的分析监测数据的变化趋势和变形信息，具有较好的变形信息检验能力和误报率较低的优点。
以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围，因此，对于本技术领域的普通技术人员来说，凡在本发明的精神和原则之内所作的任何修改、等同替换、润饰、改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于，包括以下步骤：
步骤1，获取监测体的GNSS变形监测数据，并计算其计算均值μ和标准差σ；
步骤2，设置EWMA控制图的权重参数λ和控制限限宽L；
步骤3，将获取到的监测数据，构建EWMA控制图的偏移统计量Z；
步骤4，确定异常值识别模型的中心线和上下参考限；
步骤5，对比判断EWMA偏移统计量是否超过控制限，对于超过控制限的数据定义为异常数据并进行预警。
2.根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤1中，采集对监测体(构筑物或建筑物)进行形变监测的GNSS时间序列，对时间序列进行初步分析，计算均值μ和标准差σ。
3.根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤2中，EMWA控制图的权重参数λ一般取值范围为[0.05,0.25]，为突出控制图检验性能，λ取值为0.1，而控制限限宽取值为3。
4.根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤3中，构建EWMA控制图的偏移统计量Z。
对于长时间观测的GNSS监测数据X(i),i＝1,2,…,n，n为样本大小。EWMA控制图的偏移统计量Zi定义为：
Zi＝λXi+(1-λ)Zi-1 (1)
式中λ是EWMA控制图的权值参数，是属于[0,1]中的常数，λ取值为0.1；初始值Z0(i＝1处)是过程目标值，用观测值X(i)的均值取代，Z0＝μ0。
因为Zi是EWMA所有先前样本均值的加权平均值，通过替换方程右侧的Zi-1来获得：
Zi＝λXi+(1-λ)[λXi-1+(1-λ)Zi-2]＝λXi+λ(1-λ)Xi-1+(1-λ)2Zi-2 (2)
继续以递归代替Zi-j，j＝2,3,…,t，我们得到
Figure 1
。
5.根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤4中，确定异常值识别模型的中心线CL和上参考限UCL与下参考线LCL。
GNSS观测数据X(t)是具有方差σ2的独立随机变量，则构建的EWMA统计量Zi的方差如下所示：
Figure FDA0002413266660000021
因此，通过将GNSS监测数据构建为EWMA控制图的统计量Zi，与时间(历元)共同组成EWMA控制图，其变形信息的检验准则如下：
Figure FDA0002413266660000022
CL＝μ0 (6)
Figure FDA0002413266660000023
式中L为构建检验阈值时控制限的宽度，一般取值为3；异常值识别模型的上参考限UCL、中心线CL和与下参考线LCL计算公式如式(5)，式(6)与式(7)所示；当构建的偏移统计量超过UCL或-LCL时，视为变形信息，并发出预警。
6.根据权利要求1所述的一种基于EWMA控制图的GNSS变形信息识别与预警方法，其特征在于步骤5中，对比判断由步骤3得到的EWMA偏移统计量Z是否超过上参考限UCL与下参考线LCL构成的控制限范围，对于超过控制限的数据定义为异常数据并进行预警。
CN202010183171.8A 2020-03-16 2020-03-16 一种基于ewma控制图的gnss变形信息识别与预警方法 Pending CN111400911A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010183171.8A CN111400911A (zh) 2020-03-16 2020-03-16 一种基于ewma控制图的gnss变形信息识别与预警方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010183171.8A CN111400911A (zh) 2020-03-16 2020-03-16 一种基于ewma控制图的gnss变形信息识别与预警方法

Publications (1)

Publication Number Publication Date
CN111400911A true CN111400911A (zh) 2020-07-10

Family

ID=71431042

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010183171.8A Pending CN111400911A (zh) 2020-03-16 2020-03-16 一种基于ewma控制图的gnss变形信息识别与预警方法

Country Status (1)

Country Link
CN (1) CN111400911A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112487364A (zh) * 2020-12-10 2021-03-12 北部湾大学 一种基于指数加权移动平均值算法的小泄漏检测方法
CN112581727A (zh) * 2020-11-16 2021-03-30 西人马联合测控（泉州）科技有限公司 桥梁的位移漂移预警方法、装置、设备及存储介质
CN112734977A (zh) * 2020-12-25 2021-04-30 安徽省安泰科技股份有限公司 一种基于物联网的设备风险预警系统及算法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112581727A (zh) * 2020-11-16 2021-03-30 西人马联合测控（泉州）科技有限公司 桥梁的位移漂移预警方法、装置、设备及存储介质
CN112581727B (zh) * 2020-11-16 2022-08-19 西人马联合测控（泉州）科技有限公司 桥梁的位移漂移预警方法、装置、设备及存储介质
CN112487364A (zh) * 2020-12-10 2021-03-12 北部湾大学 一种基于指数加权移动平均值算法的小泄漏检测方法
CN112487364B (zh) * 2020-12-10 2022-03-11 北部湾大学 一种基于指数加权移动平均值算法的小泄漏检测方法
CN112734977A (zh) * 2020-12-25 2021-04-30 安徽省安泰科技股份有限公司 一种基于物联网的设备风险预警系统及算法
CN112734977B (zh) * 2020-12-25 2022-07-05 安徽省安泰科技股份有限公司 一种基于物联网的设备风险预警系统及算法

Similar Documents

Publication Publication Date Title
CN111651941B (zh) 一种全球电离层电子总含量预测的算法
CN111400911A (zh) 一种基于ewma控制图的gnss变形信息识别与预警方法
Rhoades et al. Adjustment of temperature and rainfall records for site changes
CN102944583B (zh) 基于漂移补偿的金属氧化物气体传感器阵列浓度检测方法
CN110260774B (zh) 一种基于Pettitt算法的GNSS变形信息检验与预警方法
CN110441305B (zh) 一种二次喷丸的表面覆盖率检验方法
CN106021710A (zh) 基于大气电离层参数的震前卫星轨道异常识别方法
CN112733419B (zh) 一种突发大气污染事故源反演准确性提升方法
CN116910041B (zh) 一种基于尺度分析的遥感降水产品的逐日订正方法
CN106055885A (zh) 基于过采样投影近似基追踪无人机飞行数据异常检测方法
CN111367959A (zh) 一种零时滞非线性扩展Granger因果分析方法
CN107273659B (zh) 一种基于ransac算法改进的用于空间碎片光电跟踪的轨迹预测方法
Liu et al. Research on data correction method of micro air quality detector based on combination of partial least squares and random forest regression
CN116985183A (zh) 一种近红外光谱分析仪质量监测管理方法及系统
CN112949735A (zh) 一种基于离群数据挖掘的液态危化品挥发浓度异常发现方法
CN113592308A (zh) 一种基于常态模型的监测数据告警阈值提取方法
CN109614758A (zh) 具有空间相关性的圆形形状误差的监控方法
CN115600067A (zh) 一种基于物联网的建筑结构形变监测系统及方法
CN114001887B (zh) 一种基于挠度监测的桥梁损伤评定方法
KR101194677B1 (ko) 도시 지역 고농도 미세먼지의 발원지 추정 방법
CN114943493A (zh) 一种耕地质量等别监测评价系统及方法
CN114563771A (zh) 基于聚类分析的双阈值激光雷达云层检测算法
CN112116014A (zh) 一种配电自动化设备测试数据离群值检测方法
CN114019139A (zh) 一种农用地土壤重金属异常数据的探测方法
CN114299693A (zh) 一种基于gnss的边坡监测预警方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200710

WD01 Invention patent application deemed withdrawn after publication