CN111814110B - 一种桥梁健康监测数据控制图分析方法 - Google Patents

一种桥梁健康监测数据控制图分析方法 Download PDF

Info

Publication number
CN111814110B
CN111814110B CN202010440795.3A CN202010440795A CN111814110B CN 111814110 B CN111814110 B CN 111814110B CN 202010440795 A CN202010440795 A CN 202010440795A CN 111814110 B CN111814110 B CN 111814110B
Authority
CN
China
Prior art keywords
data
bridge
probability
value
interval
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.)
Active
Application number
CN202010440795.3A
Other languages
English (en)
Other versions
CN111814110A (zh
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.)
Guangdong Jianke Innovation Technology Research Institute Co ltd
Original Assignee
Guangdong Jianke Innovation Technology Research Institute Co ltd
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 Guangdong Jianke Innovation Technology Research Institute Co ltd filed Critical Guangdong Jianke Innovation Technology Research Institute Co ltd
Priority to CN202010440795.3A priority Critical patent/CN111814110B/zh
Publication of CN111814110A publication Critical patent/CN111814110A/zh
Application granted granted Critical
Publication of CN111814110B publication Critical patent/CN111814110B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/206Drawing of charts or graphs

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种桥梁健康监测数据控制图分析方法,其包括如下步骤:步骤1)将数据按照时间周期分成不同的数据块区间;步骤2)基于该数据块区间绘制控制图,通过保存形成该控制图的相关数据来保存桥梁健康监测历史数据;所述控制图包括如下信息:测点、周期时间、最大值,最小值,基于设定置信区间的频率直方图,最大概率区间中程数。本发明采用控制图方法对桥梁的这些重要数据进行提炼、保存和展示,提炼的数据已经可以用于分析桥梁结构长期运营的过程,所以,不仅保存了计算周期中的关键数据,有效地避免了海量数据的存储问题,而且通过直观的数学分析模型更好地表达了桥梁的周期变化。

Description

一种桥梁健康监测数据控制图分析方法
技术领域
本发明属于数据统计分析技术领域,尤其涉及一种桥梁监测数据的统计分析方法。
背景技术
桥梁在运营过程中,由于受到各种荷载、材料本身的退化以及自然条件因素的影响,桥梁的使用寿命以及安全受到了严峻的考验。传统的养护管理方法已经不能满足及时了解桥梁安全状况的要求,因此,桥梁结构健康监测系统应运而生。该系统运用现代传感技术与通信技术,可实时获取桥梁结构状态和环境信息的各种数据,桥梁工作者可利用数据分析评估桥梁结构的安全状态。这些监测数据是桥梁结构健康监测系统的核心和关键。
目前,桥梁结构健康监测系统中桥梁测点多,监测频次高,每天产生的数据呈爆炸性趋势增长。海量数据的累加,在数据储存、管理、访问和分析方面都面临着巨大的挑战。由于数据信息体量庞大,几乎很少会有人回看桥梁的历史数据,造成大量有效资源的浪费。往往都是新数据取代旧数据,旧数据丧失了参考价值,从而影响到桥梁的安全评估。
如何高效长期地保存桥梁运营过程中的重要数据,解决数据泛滥导致历史数据应用率低的问题,成为桥梁结构健康监测系统中面临最紧迫的问题之一。数据规模越庞大,越需要有效统一的分析平台,针对桥梁健康监测数据,该平台需要满足既能够看到桥梁目前的状况,又可以了解桥梁的发展变化,让数据充分有效地利用,展现其生命力与价值。
发明内容
前面已经提到,桥梁健康监测历史数据的泛滥导致数据利用率低,而且带来存储、管理、访问等各方面的问题,针对这些问题,本发明提供一种适于长期保存桥梁健康监测数据的方法。
本发明技术案如下:一种桥梁健康监测数据控制图分析方法,包括如下步骤:
步骤1)将数据按照时间周期(如日、周、月、年)分成不同的数据块区间;
步骤2)基于该数据块区间绘制控制图,通过保存形成该控制图的相关数据来保存桥梁健康监测历史数据;
所述控制图包括如下信息:测点、周期时间、最大值,最小值,基于设定置信区间的频率直方图,最大概率区间中程数。
最大值,最小值,以及最大概率区间中程数为桥梁健康监测数据中最关键的三个数据。本发明采用控制图方法对桥梁的这些重要数据进行提炼、保存和展示,提炼的数据已经可以用于分析桥梁结构长期运营的过程,所以,不仅保存了计算周期中的关键数据,有效地避免了海量数据的存储问题,而且通过直观的数学分析模型更好地表达了桥梁的周期变化。本发明方法可高效地长期保留桥梁运营中的重要数据,了解桥梁的变化趋势,方便数据管理以及查看。
保存的数据还包括在所述设定置信区间内的数据的均值和标准差,所述控制图展示的信息还包括根据所述均值与标准差推断的桥梁荷载超限概率以及破坏概率。
按照统计分析方法推断桥梁荷载超限概率以及破坏概率,从而使上述控制图方法不仅可以用于分析桥梁结构长期运营的过程,还可为桥梁管养人员了解桥梁的变化趋势提供相应信息。桥梁监测系统中,实时数据库中记录的信息通常包括:测点、时间、实测值。采用本发明控制图方法后,控制图数据库中记录的信息包括:测点、周期时间、最大值、最小值、最大概率区间中程数、均值、标准差、直方图频数(以划分10个小区间为例)。按一天时间来计算,一般桥梁健康监测测点的采样频率大约为30s/次,即一天时间内一个测点实时数据库要记录2880个实测值,而控制图数据库只需要记录15个值,控制图通过这15个值即可绘制完成。可见,本发明方法能极大缩小数据量,有利于数据的长期保存。
所述最大值,最小值,基于设定置信区间的频率直方图,最大概率区间中程数,在同一坐标系中标识。
所述控制图包括日控制图,所述日控制图排列在一个坐标系中,其中,横坐标为日期,纵坐标为监测值的大小。
所述桥梁荷载超限概率以及破坏概率计算方式如下:
正态概率密度函数的数学表达为:
式中,
x为随机变量;
P(x)为特定值的概率密度;
为均值,/>
σ为标准差,
根据每个计算周期的数据,通过上式得到相应的正态概率密度函数,针对荷载超限概率以及破坏概率通过下式积分进行计算:
当t=设计值时,P所求得的值为荷载超限概率;
当t=极限值时,P所求得的值为破坏概率;
设计值和极限值通过设计规范以及有限元仿真计算得到。
相比于现有技术,本发明具有如下技术效果:
本发明方法可有效提取桥梁健康监测历史数据中的重要数据,极大的缩小了数据量,有利于桥梁监测数据的长期保留,而且,本发明方法通过直观的数学分析模型表达桥梁的周期变化,方便快速了解桥梁结构参数的变化情况,分析桥梁结构长期运营的过程和预测发展趋势,让数据可以得以充分有效地利用,从而更好地为桥梁管理者提供决策服务。
附图说明
图1为形成的日控制图实例;
图2为图1中某一日的日控制图;
图3为绘制的频率直方图;
图4为绘制的概率密度函数曲线。
具体实施方式
下面通过具体实施例对发明桥梁健康监测数据控制图分析方法进行详细介绍。
本发明控制图分析方法主要包括两个步骤,步骤1)将数据按照时间周期分成不同的数据块区间;
步骤2)基于该数据块区间绘制控制图,通过保存形成该控制图的相关数据来保存桥梁健康监测历史数据;这里的控制图包括如下信息:测点、周期时间、最大值,最小值,基于设定置信区间的频率直方图,最大概率区间中程数。具体介绍如下:
控制图的计算周期根据时间可分为日控制图、周控制图、月控制图和年控制图等。控制图的形态取决于每个计算周期中的基本数据,即:最大值、最小值、直方图区间、最大概率区间中程数、置信区间以及概率密度分布函数。从控制图的形态中可以判断出桥梁随时间(日、周、月、年)周期变化的情况。
以日控制图为例,直方图是将每日的数据,以采样频率30s/次为例,2880个数据筛选至95%置信区间后,按照数据所处的区间进行数量的统计然后转化的频率直方图。
图1为日控图,日控制图排列在一个坐标系中,其中,横坐标为日期,纵坐标为监测值的大小。本发明日控制图可以起到一个统计作用,方便查看测点每日桥梁所处范围,是否在正常区间,状态是否平稳等。图2为图1中某一日的日控制图。如图2所示,各个位置的名称以及代表的含义如下:
1——当天数据最大值dmaxWeight;
2——当天数据最小值dminWeight;
3——95%置信区间上限fmaxWeight;
4——95%置信区间下限fminWeight;
5——最大概率区间上限sectionMax;
6——最大概率区间下限sectionMin;
7——最大概率区间中程数(sectionMax+sectionMin)/2。
它们都标识在同一坐标系中。
同样以日控图为例,当晚上24点某类传感器监测完毕后,对数据库中所存储的当日数据x1,x2……xn进行统计分析,记当天数据最大值dmaxWeight,当天数据最小值dminWeight;
置信区间通俗而言即所要求达到的可信度所跨度的范围,样本的置信区间可作为总体均值的一个区间估计。置信区间的计算方法:
第一步:求得数据的平均值及标准误差/>
第二步:确定置信区间,以置信度等于95%为例,当置信度等于95%时,Z=1.96;
确定置信区间
根据配置的等份Y,将上述置信区间分成Y组,记组距s=(fmaxWeight-fminWeight)/Y。根据采样频率可对数据量作出合理的配置等份。监测数据越多,就可以不断减小组距。直方图的组距越小,组数越多,数据分布的效果越明显。直方图的高度表示频数,即该变量值在该区间内出现的次数。
第三步:按照组距统计落入每个小区间的数据的数量,根据该数据量除以当天总数据量,计算落入该小区间的概率,进行排序后取出概率最大值,作为sectionProportion,该区间的上下限分别作为sectionMax、sectionMin,(sectionMax+sectionMin)/2即为最大概率区间中程数。所以,形成上述控制图,还要保留直方图的频数数据。
本发明控制图不仅保留了桥梁监测数据中最关键的三个数据:最大值,最小值,以及最大概率区间中程数,而且,还保留了形成基于置信区间的频率直方图数据。根据频率直方图可直观判断数据的分布状态,数据是用于对桥梁进行状态评估的,基于置信区间是为了提高数据的可信度和精度。上述数据是桥梁健康监测数据中真正重要的数据,如最大概率区间中程数可定性判断桥梁状态的发展趋势。本发明提炼出的这些数据,占有储存空间小,方便长期保留,形成控制图,直观方便查看,让这些历史数据能充分展现其生命力与价值。
此外本发明对该监测参数的超限概率以及破坏概率通过概率统计方法进行分析并列入上述控制图中,供桥梁管养人员参考。
趋势预测(超限概率,破坏概率)
根据科学试验和工程试验的结果,数据几乎都遵循正态分布。假设每日传感器的数据符合正态分布,取置信区间为95%的区间数据,采用正态概率密度函数对荷载超限概率以及破坏概率进行分析。正态概率密度函数的数学表达为:
式中,
x——随机变量
P(x)——特定值的概率密度
——均值,/>
σ——标准差,
根据每个计算周期的数据,通过上式得到相应的正态概率密度函数,针对荷载超限概率以及破坏概率通过下式积分进行计算:
当t=设计值时,P所求得的值为荷载超限概率。
当t=极限值时,P所求得的值为破坏概率。
设计值和极限值通过设计规范以及有限元仿真计算得到。
利用该指标模型可以对桥梁测点位置状态的发展趋势进行评价,可适用于多种不同类型的监测参数。
数据实例
给出一组测试数据,共120个。如表1:
表1
step1:
当天数据最大值dmaxWeight=1.2286
当天数据最小值dminWeight=-0.4060
平均值
标准误差SE=0.2922。
step2:
95%置信区间上限fmaxWeight=1.1337
95%置信区间下限fminWeight=-0.0117
从120个数据中筛选出置信区间的数据,共116个数据符合要求。如表2:
表2
Step3:
按照上面第三步绘制上述频率直方图,如图3所示(本数据为测试数据,工程中的数据有正有负,如正代表受拉,负代表受压)。其中,最大概率区间上限sectionMax=0.6;最大概率区间下限sectionMax=0.4;最大概率区间中程数为0.5。
趋势预测数据处理:
取置信区间为95%的区间数据,该区间内的
平均值
标准差σ=0.2583。
根据公式(2),可绘制相应的概率密度函数曲线,如图4所示。
设t1=1.3设计值时,P所求得的值为荷载超限概率。结合(3)可得
荷载超限概率P=0.25%。
设t2=1.5极限值时,P所求得的值为破坏概率。结合(3)可得
破坏概率为P=0.0174%。
本发明控制图数据库只需要记录如下信息:测点、周期时间、最大值、最小值、最大概率区间中程数、均值、标准差、直方图频数。基于原始数据形成相应的周控制图、月控制图以及一年周期控制图后(形成过程同日控制图),如果桥梁运营状况良好,可以适时摒弃数据库中的原始数据。所以,本发明能极大的缩小数据量,解决历史数据泛滥引起的各种问题。本发明控制图还呈现了荷载超限概率和破坏概率。本发明在高效长期保留桥梁运营中的重要数据的同时,方便了数据的管理和查看,而且可给桥梁管养人员了解桥梁的变化趋势提供相应信息。
桥梁的应力、应变、挠度、风力风速、位移、倾角等结构参数,均可以用本发明控制图分析方法来进行描述,本发明方法具有优良的包容性。

Claims (3)

1.一种桥梁健康监测数据控制图分析方法,其特征在于,其包括如下步骤:
步骤1)将数据按照时间周期分成不同的数据块区间;
步骤2)基于该数据块区间绘制控制图,通过保存形成该控制图的相关数据来保存桥梁健康监测历史数据;
所述控制图包括如下信息:测点、周期时间、最大值,最小值,基于设定置信区间的频率直方图,最大概率区间中程数;
保存的数据还包括在所述设定置信区间内的数据的均值和标准差,所述控制图展示的信息还包括根据所述均值与标准差推断的桥梁荷载超限概率以及破坏概率;
所述桥梁荷载超限概率以及破坏概率计算方式如下:
正态概率密度函数的数学表达为:
式中,
x为随机变量;
P(x)为特定值的概率密度;
为均值,/>
σ为标准差,
根据每个计算周期的数据,通过上式得到相应的正态概率密度函数,荷载超限概率以及破坏概率具体通过下式积分进行计算:
当t=设计值时,P所求得的值为荷载超限概率;
当t=极限值时,P所求得的值为破坏概率;
设计值和极限值通过设计规范以及有限元仿真计算得到。
2.根据权利要求1所述的桥梁健康监测数据控制图分析方法,其特征在于,所述最大值,最小值,基于设定置信区间的频率直方图,最大概率区间中程数,在同一坐标系中标识。
3.根据权利要求2所述的桥梁健康监测数据控制图分析方法,其特征在于,所述控制图包括日控制图,所述日控制图排列在一个坐标系中,其中,横坐标为日期,纵坐标为监测值的大小。
CN202010440795.3A 2020-05-22 2020-05-22 一种桥梁健康监测数据控制图分析方法 Active CN111814110B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010440795.3A CN111814110B (zh) 2020-05-22 2020-05-22 一种桥梁健康监测数据控制图分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010440795.3A CN111814110B (zh) 2020-05-22 2020-05-22 一种桥梁健康监测数据控制图分析方法

Publications (2)

Publication Number Publication Date
CN111814110A CN111814110A (zh) 2020-10-23
CN111814110B true CN111814110B (zh) 2024-05-10

Family

ID=72848338

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010440795.3A Active CN111814110B (zh) 2020-05-22 2020-05-22 一种桥梁健康监测数据控制图分析方法

Country Status (1)

Country Link
CN (1) CN111814110B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112800602B (zh) * 2021-01-25 2023-05-23 国家能源集团新疆吉林台水电开发有限公司 一种安全监测数据的整体可视化分析方法
CN113626906B (zh) * 2021-07-07 2024-03-22 武汉轻工大学 基于数据挖掘与控制图的地下支护结构变形监控方法
CN114004467A (zh) * 2021-10-15 2022-02-01 河南工业大学 一种基于监测数据的预制装配式桥梁结构性能分析方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382473A (zh) * 2008-10-08 2009-03-11 重庆大学 桥梁结构安全预警的ewma控制图方法
CN101382474A (zh) * 2008-10-08 2009-03-11 重庆大学 一种桥梁结构安全的综合评价方法
CN101696912A (zh) * 2009-11-02 2010-04-21 重庆交通大学 基于统计指标的混凝土桥梁结构健康诊断的方法
CN103985415A (zh) * 2013-02-10 2014-08-13 Lsi公司 基于保留漂移历史的非易失性存储器读取阈值最优化
CN105184065A (zh) * 2015-08-27 2015-12-23 北京特希达交通设施顾问有限公司 基于常态均值的桥梁损伤识别方法
CN105488352A (zh) * 2015-12-11 2016-04-13 河南省交通科学技术研究院有限公司 基于长期挠度监测数据的混凝土桥梁刚度可靠度评估方法
CN110017929A (zh) * 2019-03-12 2019-07-16 华中科技大学 基于子结构灵敏度分析的船撞桥荷载与损伤同步识别方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI669617B (zh) * 2018-10-12 2019-08-21 財團法人工業技術研究院 設備健康狀態監控方法及其系統

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382473A (zh) * 2008-10-08 2009-03-11 重庆大学 桥梁结构安全预警的ewma控制图方法
CN101382474A (zh) * 2008-10-08 2009-03-11 重庆大学 一种桥梁结构安全的综合评价方法
CN101696912A (zh) * 2009-11-02 2010-04-21 重庆交通大学 基于统计指标的混凝土桥梁结构健康诊断的方法
CN103985415A (zh) * 2013-02-10 2014-08-13 Lsi公司 基于保留漂移历史的非易失性存储器读取阈值最优化
CN105184065A (zh) * 2015-08-27 2015-12-23 北京特希达交通设施顾问有限公司 基于常态均值的桥梁损伤识别方法
CN105488352A (zh) * 2015-12-11 2016-04-13 河南省交通科学技术研究院有限公司 基于长期挠度监测数据的混凝土桥梁刚度可靠度评估方法
CN110017929A (zh) * 2019-03-12 2019-07-16 华中科技大学 基于子结构灵敏度分析的船撞桥荷载与损伤同步识别方法

Also Published As

Publication number Publication date
CN111814110A (zh) 2020-10-23

Similar Documents

Publication Publication Date Title
CN111814110B (zh) 一种桥梁健康监测数据控制图分析方法
CN114168906B (zh) 一种基于云计算的测绘地理信息数据采集系统
CN106529704A (zh) 月最大电力负荷预测方法及装置
CN110782153A (zh) 一种企业园区综合能效评估体系建模方法及系统
CN112084684B (zh) 一种基于物联网的桥梁健康可视化监测系统
CN109308589B (zh) 电网自动化数据质量监测方法、存储介质、终端设备和系统
CN117932501B (zh) 一种电能表运行状态管理方法和系统
CN113887898A (zh) 一种基于电力数据指标的区域产业发展雷达图评价方法
CN114254806A (zh) 配电网重过载预警方法、装置、计算机设备和存储介质
CN113657610A (zh) 一种基于随机森林的冰雹气候特征预测方法
CN116894687A (zh) 一种基于机器学习的用电量分析方法及系统、电子设备
CN115115470A (zh) 一种基于排放因子法的绿色数据中心碳排放管理方法
CN107045548B (zh) 一种计算风电能量利用率的系统及方法
CN115564093A (zh) 一种面向化工区工业的用水预测、预警方法
CN112434750B (zh) 基于卷积神经网络的大坝监测数据发展模式识别方法
CN112257233B (zh) 弹性电网恢复力评估方法、装置、计算机设备和介质
CN112001551B (zh) 一种基于大用户电量信息的地市电网售电量预测方法
CN111368257B (zh) 煤改电负荷特性的分析预测方法及装置
CN114298131A (zh) 一种基于水电设备状态指标特征值进行状态监测的方法
CN114169802A (zh) 电网用户需求响应潜力分析方法、系统及存储介质
CN111079069A (zh) 一种基于误差分布的预报难度计算方法及系统
CN113077105B (zh) 一种长假日负荷预测方法及装置
CN117955245B (zh) 电网的运行状态的确定方法、装置、存储介质和电子设备
CN113656238B (zh) 一种智能终端的ai工业应用能力测试方法和系统
CN118228510A (zh) 基于实时数据采集的材料腐蚀仿真建模方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20210318

Address after: 528437 room 409-17, 4 / F, building C, innovation center building, 34 Xiangshan Avenue, Cuiheng New District, Zhongshan City, Guangdong Province

Applicant after: Guangdong Jianke Innovation Technology Research Institute Co.,Ltd.

Address before: 8 / F, new office building, Jianke academy, 121 Xianlie East Road, Tianhe District, Guangzhou City, Guangdong Province, 510500

Applicant before: GUANGDONG PROVINCIAL ACADEMY OF BUILDING RESEARCH GROUP Co.,Ltd.

SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant