一种基于增量判别的水电工程安全监测数据粗差处理方法
【专利摘要】本发明涉及一种基于增量判别的水电工程安全监测数据粗差处理方法，属于水电工程安全监测技术领域。该方法对监测数据求取两相邻测值的增量形成新的增量数据序列，在对该增量序列依次进行莱茵达准则判别，若满足增量序列测值和下一增量序列测值符合莱茵达准则、下下一增量序列测值不符合莱茵达准则的条件，则判断该增量序列测值所对应的原始测值为尖点异常值，应予以剔除并重新插值补充。本方法方法可有效识别并剔除尖点异常值，解决传统人工识别效率低下易错的问题，并提高了生产效率，节约了生产成本，为监测报表等后续工作提供了坚实基础。
【专利说明】
-种基于増量判别的水电工程安全监测数据粗差处理方法
技术领域
[0001] 本发明属于水电工程安全监测技术领域，具体设及一种基于增量判别的水电工程 安全监测数据粗差处理方法。
【背景技术】
[0002] 近年来我国水电发展尤为迅速，一大批水电枢纽已完建开始投入运行，由于水电 工程规模庞大，关系到库区和下游人民生命财产安全，因而无论是在建设期还是运行期，都 需定期评估大巧的运行性态。生产实际中，由于环境扰动、信号干扰、测量误差等各种原因 导致大巧安全监测数据出现异常，为了有效识别并剔除异常值W便为下一步的建筑物预警 及报表生成提供依据，因而有必要研究一种可W针对水电工程安全监测数据粗差处理的方 法。
【发明内容】
[0003] 本发明的目的是为了解决现有技术的不足，提供一种基于增量判别的水电工程安 全监测数据粗差处理方法，W有效识别并剔除异常值，为后续工作提供扎实基础。
[0004] 为实现上述目的，本发明采用的技术方案如下：
[0005] -种基于增量判别的水电工程安全监测数据粗差处理方法，包括如下步骤：
[0006] 步骤（1)，对于工程安全监测中的某一测点的监测数据序列xi，X2,…，xn(n>10)， 依次求取两相邻测值的增量A Xi = XW-XiU = 1，2,…，n-1)，得到增量数据序列A XI，A X2 , ... , A Xn-I;
[0007] 步骤（2 )，求取增量数据序列的平均值
和标准差
，并对增量数据A Xi依次进行莱茵达准则判别（亦称30准则），具体判
别万法如h :
[000引a.若增量数据序列值A Xi不满 ，则进行下一增 量数据序列值的检查；
[0009] b.若增量数据序列值A Xi满足 i下一增量数据 序列值A Xi+地满足该不等式，且下下一测值A XW不满足该不等式，则增量序列测值A Xi 所对应的原测值Xi + l为尖点异常测值，予W剔除，并将在原测值位置替换成新测值
[0010] C .若增量数据序列值A Xi满足
但下一测值A XW不满足该不等式，则进行下一增量数据序列值的检查；
[0011] 按照W上判别方法依次对所有增量数据序列值进行判别，得到新的监测数据序 列。
[0012] 通过对测值两相邻增量形成的新序列进行莱茵达法则判别，根据判别结果确定某 测值是否为异常值，从而达到对监测数据粗差处理的目的。对测值而言，所述的测点可W是 测缝计、应变计、多点位移计、测斜仪、压应力计等的监测数据，但不仅限于上述监测仪器。 [OOU]监测数据为测缝计、应变计、多点位移计、测斜仪或压应力计监测得到的数据。
[0014] 常规30准则为针对若干次独立等精度测量结果X1，X2,…，Xn计算其算术平均值
，若测量结果剩余误差|x, - ^大于3倍标准误差，良
则表示该测量结果为异常数据，由于常规30准则基于测值属于正态分布类型，而从实际情 况来看，常规30准则判断相关联性差，容易导致误判，误判率较高。
[0015] 本发明与现有技术相比，其有益效果为：
[0016] 本发明通过建立一种基于增量判别的水电工程安全监测数据粗差处理方法，有效 识别并剔除了水电工程安全监测数据中的尖点异常值，解决了传统人工识别效率低下易错 的问题，为监测数据的后续处理提供了坚实基础。
[0017] 本发明针对海量的监测数据，只需输入时间即可自适应捜索最匹配的原始记录， 并对原始记录进行粗差处理，识别并剔除尖点异常值，本发明方法将大大提高水电工程在 施工期、运行期内的安全监测数据粗差的识别效率，并极大提高了生产效率，节约了生产成 本，为监测报表等后续工作提供了坚实基础。同时，对于中国西部完建、在建、拟建万余座水 库大巧，本发明的市场需求和经济效益均较为显著。
[0018] 使用传统人工方法对100万条测缝计开合度数据导入到Excel进行粗差处理大致 需要1人工日，而使用前述方法编制计算机程序后，处理相同数据仅需要93秒，加上后续人 工检验也只需0.2人工日，提高生产效率80% W上，且处理准确性也大大提高，降低误判率 80% W上。由于大巧施工期及运行期每天都会有大量的监测数据产生，每周都需要对数据 进行粗差处理，按每次减少工作量80%估计，按每人工日400元计，对一座水电站而言，每月 节省生产成本1280元，年节约成本15360元。
【附图说明】
[0019] 图1是本发明实施例中测缝计开合度测值过程线。
[0020] 图2是本发明实施例中经过粗差处理后的测缝计测值过程线。
【具体实施方式】
[0021 ]下面结合实施例对本发明作进一步的详细描述。
[0022] 本领域技术人员将会理解，下列实施例仅用于说明本发明，而不应视为限定本发 明的范围。实施例中未注明具体技术或条件者，按照本领域内的文献所描述的技术或条件 或者按照产品说明书进行。所用试剂或仪器未注明生产厂商者，均为可W通过购买获得的 常规产品。
[0023] 表1是本发明实施例中某水电站22#巧段某测点（C4-A22-J-17)的部分监测数据 表。
[0024] 表2是本发明实施例中增量序列值及其偏差值表。
[0025] 表3是本发明实施例中经过粗差处理后的测值表。
[0026] 本发明的具体实例如下，
[0027] (A)针对给定的时间范围，自适应捜索并导出数据库中最接近给定时间范围的原 始记录；
[0028] 某水电站2?!巧段某测点（C4-A22-J-17)的部分监测数据如表1所示。W下需要导 出2013年10月5日7:00~2013年10月26日9:00的数据。
[0029] 表1某水电站2?巧段某测点编号为C4-A22-J-17的部分监测数据
[0030]
[003。 捜索表1的原始数据，与给定时间段2013年10月5日7:00~2013年10月26日9:00最 匹配的数据为序号2~25的数据，将其导出，绘制如图1所示的测值过程线。
[0032] (B)针对给出的数据序列xi，X2,…，xn(n>10)依次进行判断，分别计算相邻两测值 的增量 A Xi = Xi+1-Xi(i = 1，2，…，n-1)，得到增量数据序列 A Xi，A X2，…，A Xn-I，即（X2-X1)， (X3-X2)，…，（Xi-Xi-I)，…，（Xn-Xn-I)，求取该增量序列的平均{1
和标准偏差
[0033] 对表1中的序号2~25的数据形成新测值序列xi，X2,…，xn(n = 24)依次求取相邻两 测值的增量A xi = Xi+：L-Xi(i = l，2,…，n-l)，形成新的增量数据序列A XI，A X2，。'，A xn-i(n = 24)如下：
[0034] -0.01，0.01，-0.01，-0.01，0.02,0,0,0,0.03，-0.04,0,41.57，-41.55，-0.02， 0?03，0，-〇?Ol，0，0?03，-〇?03，0，-〇?02，0?02〇
[0035] 上述增量数据的平均值巧
[0036] H术说吾狱报席别的3倍掠准偏差责
[0037] (C)对该增量序列的每个测值依次进行莱茵达准则判别（亦称30准则），即
，式中A Xi为当前测值；
[0038] 具体判别方法如下：
[0039] a.若增量数据序列值A Xi不满 ，则进行下一增 量数据序列值的检查；
[0040] b.若增量数据序列值A Xi满盾 i下一增量数据 序列值A Xi+地满足该不等式，且下下一测值A XW不满足该不等式，则增量序列测值A Xi 所对应的原测值Xi + l为尖点异常测值，予W剔除，并将在原测值位置替换成新测值
[0041 ] C.若增量数据序列值A Xi满^
,但下一测值A XW不满足该不等式，则进行下一增量数据序列值的检查；
[0042] 按照W上判别方法依次对所有增量数据序列值进行判别，得到新的监测数据序 列。
[0043] 依次计算上述增量序列中每个测值及其偏差值- A￥|，如表2所示。
[0044] 表2增量数据序列值及相应的偏差值表
[0045]
[0046] (D)通过比较每行
k小可知，连续两序列值41.57 和-41.55偏差值超过了 3倍标准差，紧后序列值-0.02未超过上述标准，属于前述判别方 法中的b情况，属于尖点异常数据，应剔除原始测值序列Xi中的第14行数据，即表1中 2 0 1 3 / 1 0/ 1 6 8 : 00所对应的测值42mm。并将在原测值位置替换成新测值
[0047] 测值-41.55虽满足偏差超过3倍标准差，但其紧后测值-0.02不满足，则不属于前 述b情况，属于前述C情况，不予剔除。
[0048] 化)经过粗差处理后的测缝计测值数据如表3所示，其开合度过程线如图2所示。
[0049] 表3经过粗差处理后的测值数据
[(K)加 ]
[0051]对大巧安全监测而言，每种仪器的监测数据（测缝计、应变计、多点位移计、测斜 仪、压应力计等等)在生成报表之前都需要对其进行粗差处理，筛选出明显的异常值方能进 行下一步工作。而大巧安全监测规范要求大巧安全监测周报、月报W及年报均需要对该时 间段内的监测数据进行粗差处理、整编后方可用于生成报表。前述针对测缝计开合度的监 测数据进行粗差处理仅用于说明该方法的步骤，并不局限于测缝计。因而本发明的处理方 法可对大巧施工期和运行期的大巧安全监测发挥重要作用。
[0052] 使用传统人工方法对100万条测缝计开合度数据导入到Excel进行粗差处理大致 需要1人工日，而使用前述方法编制计算机程序后，处理相同数据仅需要93秒，加上后续人 工检验也只需0.2人工日，提高生产效率80% W上，且处理准确性也大大提高，降低误判率 80% W上。由于大巧施工期及运行期每天都会有大量的监测数据产生，每周都需要对数据 进行粗差处理，按每次减少工作量80%估计，按每人工日400元计，对一座水电站而言，每月 节省生产成本1280元，年节约成本15360元。
[0053] W上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术 人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本 发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，运些变 化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其 等效物界定。
【主权项】
1. 一种基于增量判别的水电工程安全监测数据粗差处理方法，其特征在于，包括如下 步骤：步骤（1)，对于工程安全监测中的某一测点的监测数据序列X1，X2,-SXn(H)IO),依次 求取两相邻测值的增量Δ Xi = xi+1-Xi(i = l，2,…，n-l)，得到增量数据序列Δ X1, Δ X2，···， A Xn-I ； 步骤(2)，求取增量数据序列的平均_然 后对增量数据A X1依次进行莱茵达准则判别;具体判别方法如下： a. 若增量数据序列值△ X1不满>，则进行下一增量数 据序列值的检查；值Δ Xl+1也满足该不等式，且下下一测值△ Xl+2不满足该不等式，则增量序列测值△ Xl所对 应的原测值Xl + 1为尖点异常测值，予以剔除，并将在原测值位置替换成新测值 b. 若增量数据序列值△ Xi满， ，且下一增量数据序列c .若增量数据序列值Δ Xi满足 ，但下一测值Δ xi+1f 满足该不等式，则进行下一增量数fe/予列俚的恨宜； 按照以上判别方法依次对所有增量数据序列值进行判别，得到新的监测数据序列。2. 根据权利要求1基于增量判别的水电工程安全监测数据粗差处理方法，其特征在于， 监测数据为测缝计、应变计、多点位移计、测斜仪或压应力计监测得到的数据。
【文档编号】G06Q50/06GK106022966SQ201610372702
【公开日】2016年10月12日
【申请日】2016年5月31日
【发明人】张宗亮, 艾永平, 易魁, 陈亚军, 赵志勇, 王永祥, 张文学, 鲁俊兵, 邱小弟, 张礼兵, 陈豪, 李然, 张鹏, 胡灵芝, 张帅, 廖贵能, 彭欣欣, 庞博慧, 陈鸿杰
【申请人】中国电建集团昆明勘测设计研究院有限公司, 华能澜沧江水电股份有限公司