专利名称：一种桥梁结构安全的综合评价方法
技术领域：
本发明涉及一种桥梁结构安全预警技术，尤其涉及一种桥梁结构安全的综合评价方法。
背景技术：
由于桥梁结构本身和所处环境的复杂性以及结构安全评价的特殊性，现有的结构安全评 价方法存在较大的局限，在实际应用上未能取得令人满意的进展，主要体现在如下几点
(1) 现有的基于模型的桥梁结构安全评价理论均来源于力学学科，其评价过程严重依赖 于结构的精确理论模型和确定的系统激励；然而旧桥资料的严重缺失和桥梁结构参数在运行 过程中的不断变化使得结构精确模型的建立极为困难，加上桥梁结构健康监测中随机环境激 励的未知性，使得基于模型的结构安全评价遇到难以逾越的障碍。
(2) 基于统计分析的评价方法虽已得到逐歩的重视，但其研究和应用尚处在初期；桥梁尺 寸巨大，结构复杂，振动幅值很低，统计分析通常所采用的动力特性指标易被环境噪声淹没 或对局部损伤不敏感，而且没有将桥梁结构安全评价与环境特性(如温度)有效地联系起来， 其评价效果还无法满足工程要求。
(3) 神经网络虽然在处理非线性问题时具有优良的性能，但由于缺乏完备的训练样本，无 法在结构响应和安全状态之间建立有效的因果联系。
(4) 目前的安全评价大多采用单一的方法或单一的参数进行分析，环境影响的不确定性以 及测试数据的不完整使得监测变量未能全面反映结构的安全状况。
以上问题的存在，使得目前还没有一种令人满意的桥梁安全综合评价体系。
从以上分析可知，由于桥梁结构本身的构成极为复杂，其性能容易受到结构所处环境的 各种因素的影响，使得结构安全的评价工作也十分的复杂和困难。尽管目前针对桥梁结构安 全的评价方法有很多，但由于各种原因，这些评价方法都有其本身固有的缺陷，在实际应用 上尚未取得令人满意的结果。
D—S证据理论由Dempster提出，利用概率的上下限表示不确定性，后经Shafer进一歩 发展和完善，使之成为一种系统化的不确定性推理理论。常用于信息融合的决策层。
在实际应用中，不同的证据(本文中指评价结果)代表不同的信息来源，采集的手段不 同、时间不同，或采用的评估方法不同，都可能导致评价结果权重的不同分布。在进行融合 时，采用D-S融合规则，即基于直接和(正交和)的融合方法。这样，融合的结果可以大大 减小系统输出的不确定度，进一歩提高状态评估的准确性。
发明内容
本发明公开了一种桥梁结构安全的综合评价方法，它包括1)采用『i /指标和/ S/指标
对桥梁结构安全进行综合评价并形成评价结果；2)利用EWMA控制图的初级预警限和安全 告警限对桥梁结构安全进行综合评价并形成评价结果；3)采用Dempster融合规则对『及/指 标和及S/指标形成的评价结果和EWMA控制图的初级预警限和安全告警限形成的评价结果 进行融合并形成融合结果；根据融合结果来确定最终是否发出告警信息，或者，直接根据歩 骤l)或步骤2)中的评价结果来确定最终是否发出告警信息。
歩骤1)中，采用『及/指标和及S/指标对桥梁结构安全进行综合评价并形成评价结果， 包括
1) i S/指标是指实时监测信息的数据在某一周期内的单向波动量占总波动量的相对百分 数，具体计算方法为<formula>formula see original document page 8</formula>式中A为实时监测信息的数据在一个周期内正方向增加量的总和；
B为实时监测信息的数据在一个周期内负方向增加量绝对值的总和；
其中<formula>formula see original document page 8</formula>式中Z,为实时监测信息的数据在一个周期内的第z'个值；
Z^为实时监测信息的数据在一个周期呢的第/-1个值；
/为实时监测信息的数据在一个周期内的取值数目；
/为正整数。
2) 『及/指标是指实时监测信息的数据当前波动幅度占其前一周期最大波动幅度的百分 数，具体计算方法为
<formula>formula see original document page 8</formula>
式中Zmax为实时监测信息的数据在前一周期内的最大值；
为实时监测信息的数据在前一周期内的的最小值； Z为实时监测信息的数据的当前数值；
n为实际测量数据结果的周期长度内的测点数量；
其中，实时监测信息的数据在前一周期内的最大值Z^和最小值Z^分别为前一周期内 的测点数W减去距当前检测点W个测点后振幅的最大值和最小值，W取值可为3~6; 『 /指标和i^/指标形成的评价结果，包括
1) - 25 125时，不发出告警请求；反之则发出告警请求；或
2)同时满足-25<WRI/n)<125和45<及RSI/(n)<555时，不发出告警请求；任一条件不满
足时，则发出告警请求；
歩骤2)中，所述的利用ewma控制图的初级预警限和安全告警限对桥梁结构安全进行
综合评价并形成评价结果，包括ewma控制图的内控制限[zcz， t/o:]为初级预警限，外
控制限[VCL， HCL]为安全告警限，其评价结果为当桥梁结构的实时监测信息控制图统计量
的值超过内控制限[VCL, UCL或外控制限[VCL, HCL]时，发出告警请求，否则不发出告警 请求。
前文所述的融合，其步骤包括
1) 根据冲突因子《的值，确定是否存在评价结果冲突的情况如果没有冲突，则直接
用Dempster融合规则对评价结果进行融合；若有冲突，则进行下一步；
2) 采用专家评分法或相关系数法计算平均评价结果的基本概率赋值5iM;
3) 根据融合的具体要求及评价结果的属性，计算各评价结果的权值；
4) 用平均评价结果和没有冲突的评价结果建立新的评价结果源，并计算新评价结果的 se4，平均评价结果的权值由冲突评价结果的权值代替；
5) 用Dempster融合规则对新评价结果进行融合，得到融合结果。
融合歩骤2)中，采用专家评分法或相关系数法计算平均评价结果的基本概率赋值5iM 的步骤，包括
采用专家评分法或相关系数法对评价结果进行预处理的步骤，包括
(1)采用专家评分法构造融合评价变量的基本概率分布的方法，包括 根据下式确定专家评分法的基本概率赋值，即Bi^:
附'
式中
m (P表示由评价结果A确定的结构安全状态基本概率赋值^^4; m(/0表示由评价结果A不能确定的结构状态基本概率赋值5尸J; m为桥梁结构某一单元的安全状态的种数；
为根据某参数的评价结果A对结构安全状态进行评分专家人数；
C,G)为专家!'根据某参数的评价结果A判断结构安全状态为y的基本概率赋值5iM;
C,(//)为专家/根据某参数的评价结果A不能确定状态的基本概率赋值
其中，A为第/个监测变量的评价结果；
(2)采用相关系数法构造融合监测变量的基本概率分布，包括-
根据下式确定相关系数的基本概率分布
<formula>formula see original document page 10</formula>式中-
/、y均为自然数；
C,(;')为传感器/与目标状态_/的相关系数； m, (P为传感器!'对目标状态y的基本概率分布； 其中<formula>formula see original document page 10</formula>式中《,为系统总的不确定系数； M为传感器数量；
",为传感器/与目标状态y的最大相关系数; A为传感器/与各相关系数的分布系数； 及,为传感器/的可靠系数； l为环境性能系数，其值位于
内； 各参数表示如下
<formula>formula see original document page 10</formula>
传感器/对辨识框架的基本概率分布为
<formula>formula see original document page 10</formula>
式中m,(//)为传感器/对辨识框架a的基本概率分布。
融合步骤5)中，用Dempster融合规则对新评价结果进行融合，得到融合结果，包括 辨识框架// = {4，4}，
其中，A表示结构处于安全状态，A表示结构处于不安全状态，它们的基本概率分布分 别为和m(X2)， w(/i)表示不确定的状态， 基于基本概率赋值的决策
对任意<formula>formula see original document page 10</formula>,满足
<formula>formula see original document page 10</formula>
若有<formula>formula see original document page 11</formula> 则不发出安全告警信息；不满足上述条件时，发出安全告警信息；
e,、 ^为预先设定的阈值。
直接根据步骤l)或歩骤2)中的评价结果来确定最终是否发出告警信息，其具体内容如
下
直接根据『及/指标和i S/指标形成的评价结果来确定最终是否发出告警信息，包括1) <formula>formula see original document page 11</formula>时，不发出告警信息；反之则直接发出告警信息；或
2)同时满足-25 125和45S及S/(w)S55时，不发出告警信息；任一条件不满
足时，则直接发出告警信息；
直接根据EWMA控制图形成的评价结果来确定最终是否发出告警信息，包括当桥梁
结构的实时监测信息控制图统计量的值超过内控制限<formula>formula see original document page 11</formula>或外控制限<formula>formula see original document page 11</formula>时，
直接发出告警信息，否则不发出告警信息。
ewma控制图的内控制限[zo:， t/o:]的计算方法，包括
1) 将桥梁结构的任一实时监测参数的历史信息建立一个随机过程，表述为
式中/为正整数；
Z,为该监测信息的数据在周期内的第/个值； Z为该监测信息的数据的过程均值； 《为随机变化的参量；
2) 确定(T、 d、和^ "的值
对监测信息的数据进行分析，计算该数据的标准差tr以及最小偏移量"=//w (其中 hlSl);再根据桥梁结构的实际情况和评价的具体要求，选取数据变化的基本时间周期的 10倍~20倍，作为在某时间段内不发生误报警的平均运行时长Ji ^ (Average Run Length)。
3) 确定义和A:的值
通过査表法，从EWMA控制图最优义曲线图和EWMA控制图、(/1， 0组合曲线图确定/1 和A:的值；
4) 优化A和A:的值
根据桥梁结构监测的实际情况，建立一个迭代过程，通过d值的变化搜索最优的(A》)组 合，使被监测过程的误发警报的概率与预先确定的概率的误差最小；
5) 根据下式求得EWMA控制图的内控制限的上、下限
<formula>formula see original document page 12</formula>
<formula>formula see original document page 12</formula>
式中"为样本子集的大小，"=^及"/水，A为采样间隔；
ewma控制图的外控制限[rcz, //o:]计算方法，包括
1) 根据桥梁结构可靠度规范给出的可靠度指标，即可获得结构失效概率的允许值Pf;
2) 根据下式计算EWMA控制图的外控制限系数[-x,x]，(x>o);<formula>formula see original document page 12</formula>
3) 根据下式确定外控制限:<formula>formula see original document page 12</formula>
式中Z为监测信息的数据的过程均值； ^T为监测信息的数据的标准差。
本发明的有益技术效果大大减小系统输出的不确定度，提高状态评估的准确性，评价 结果可靠度高、适用面广。
附图1，桥梁结构状态示意附图2， Ji Z0=50， 100,250,370时，EWMA控制图最优X曲线；
附图3， ^ i:0=500， 1000， 1500,2000时，EWMA控制图最优X曲线；
附图4， ^ Z(t50， 100,250,370时，EWMA控制图(/1, it)组合曲线；
附图5， ^ "=500,1000,1500,2000时，EWMA控制图(;i， A:)组合曲线；
附图6，最优( ^)组合的搜索过程框附图7，挠度n5的活载效应EWMA控制图评价附图8，挠度s5瞬变的EWMA控制图评价附图9，挠度s5瞬变信息失控的模拟及评价附图IO，应变20瞬变信息的EWMA控制图评价附图ll，应变30瞬变信息的EWMA控制图评价附图12，应变30失控的模拟及评价附图13，挠度s5劣化效应信息的评价附图14，应变30劣化效应信息的评价附图中，图2至图5为(/l，A)组合曲线图，均为引用图，其来源为CrowderSV. Design of exponentially weighted moving average schemes. Journal of Quality Technology 1989,
1具体实施例方式
(一)采用融合规则进行评价 采用RSI指标和及S/指标对桥梁结构安全进行综合评价后形成证据mi、 m2 (即评价结 果)，利用EWMA控制图的初级预警限和安全告警限对桥梁结构安全进行综合评价并形成证据m3 (也即评价结果)，如果m,、 m2、 m3之间不存在冲突，则直接用Dempster融合规则对证据(评价结果)进行融合，得到融合结果；如果m卜m2、 rri3之间存在冲突，则按如下歩 骤进行处理
1 )采用专家评分法或相关系数法计算平均评价结果的基本概率赋值S户A
2) 根据融合的具体要求及评价结果的属性，计算各评价结果的权值；
3) 用平均评价结果和没有冲突的评价结果建立新的评价结果源，并计算新评价结果的SiM，平均评价结果的权值由冲突评价结果的权值代替；
4) 用Dempster融合规则对新评价结果进行融合，得到融合结果。
歩骤l)中，采用专家评分法或相关系数法计算平均评价结果的基本概率赋值5i^的步 骤，包括
采用专家评分法或相关系数法对评价结果进行预处理的歩骤，包括
(1) 采用专家评分法构造融合评价变量的基本概率分布的方法，包括 根据下式确定专家评分法的基本概率赋值，即BPA:
<formula>formula see original document page 13</formula>
式中
m (j)表示由评价结果D,确定的结构安全状态基本概率赋值BPA; m(U)表示由评价结果A不能确定的结构状态基本概率赋值BPA; m为桥梁结构某一单元的安全状态的种数；
n为根据某参数的评价结果A对结构安全状态进行评分专家人数；
C,C(j)为专家/根据某参数的评价结果Z),判断结构安全状态为j的基本概率赋值BPA; C,(U)为专家i根据某参数的评价结果Di不能确定状态的基本概率赋值BPA; 其中，Di为第i个监测变量的评价结果；
(2) 采用相关系数法构造融合监测变量的基本概率分布，包括 根据下式确定相关系数的基本概率分布<formula>formula see original document page 14</formula><formula>formula see original document page 14</formula><formula>formula see original document page 14</formula>式中
i, j均为自然数；
Ci(i)为传感器i与目标状态j的相关系数；
mi(j)为传感器i对目标状态j的基本概率分布；
其中<formula>formula see original document page 14</formula> 式中: Ki为系统总的不确定系数；
Ns为传感器数量；
alphai为传感器i与目标状态j的最大相关系数；
betai为传感器i与各相关系数的分布系数；
i ,为传感器/的可靠系数；
Yi为环境性能系数，其值位于
内； 各参数表示如下
<formula>formula see original document page 14</formula>
传感器i对辨识框架u的基本概率分布为<formula>formula see original document page 14</formula>式中: mi(u)为传感器i对辨识框架u的基本概率分布。
步骤4)中，用Dempster融合规则对新评价结果进行融合，得到融合结果，包括 辨识框架<formula>formula see original document page 14</formula>
其中，A1表示结构处于安全状态，A2表示结构处于不安全状态，它们的基本概率分布分
别为和m(A1)和m(A2)表示不确定的状态， 基于基本概率赋值的决策
对任意<formula>formula see original document page 14</formula> 满足
<formula>formula see original document page 14</formula>若有<formula>formula see original document page 15</formula>
则不发出安全告警信息；不满足上述条件时，发出安全告警信息；
<formula>formula see original document page 15</formula>为预先设定的阈值。
(二 )采用『i /指标和指标进行评价
鉴于活载效应和劣化效应信息的变化特点，可通过监测一个时间段内监测变量的单向变 化程度来确定过程是否失控，从而判断结构的安全情况。显然，监测变量在某一时间段内的 单向变化程度需要一个定量的表示，这其实就是一种安全评价指标。
桥梁结构安全评价的统计指标『及/指标由证券市场技术分析指标中的威廉指标演化得 到。 .
威廉指标反映了计算日收盘价在过去一段时间里的全部价格范围内所处的相对位置。其 计算公式如下
<formula>formula see original document page 15</formula>
式中Zmax为评价变量前一周期内的最大值；
Zmm为评价变量前一周期内的的最小值；
Z为评价变量的当前数值；
W为测点数，即周期长度；
其中，评价变量前一周期内的最大值Zmax和最小值Zmin分别为前一周期内的测点数w减 去距当前检测点m个测点后振幅的最大值和最小值，w取值大于等于144， m取值为3 6。
综上所述，用WRI指标进行桥梁结构安全评价的判断方法如下
① 如果-25<WRI(n)125，则结构处于安全状态。
② 如果WRI(n)》125或WRI(n)《-25，结构可能出现问题，系统给出预警信息。 桥梁结构安全评价的统计指标及S/指标反映监测变量在某一周期内的单向波动量占总波
动量的相对百分数。以结构活载效应或劣化效应信息为监测变量时，其计算公式如下:
<formula>formula see original document page 15</formula>
其中A为该效应在一个周期内正方向增加量的总和，S为一个周期内负方向增加量的总
和(绝对值)，n为周期长度。即
<formula>formula see original document page 15</formula>
式中A为该效应信息在周期内的第i个值(J l， 2，…，77)。
S/指标其实就是反映了监测变量在某一个周期内正方向增量所产生的波动占总的波动
的百分数，其值越大，表示其向上运动的趋势越强；其值越小，则表示向下运动的趋势越明 显。显然，在周期长度为无限长的情况下，如果结构处于安全状态，则i S/指标的值必然为 50，而在选定足够长度的周期内，则其值在50附近上下小幅波动，上、下幅度不会超过10%。 因此，以此指标对桥梁结构进行安全评价时， 一般而言，如果7 S/指标的值大于55，则表明 在此周期内，监测变量(挠度)向正方向波动的趋势很强烈，此时需要发出预警信号；反之， 如果i^/指标的值小于45，则表明监测变量朝负方向波动的趋势很强，也需要发出预警信号； 如果i S/指标的值位于45至55之间，可以认为监测变量的波动趋势不明显，桥梁结构处于 正常范围内。
(三)采用EWMA控制图进行评价 桥梁结构的活载效应和劣化效应的历史信息分别构成一个随机过程^k。， &，
不妨表示为
式中，z为过程均值，《为随机变化的参量。
参见附图1，当结构处于正常状态时，z不变，《呈随机变化，整个过程没有明显的趋
势。当结构出现损伤或有安全问题时，z,的变化幅度持续增大，过程均值出现持续的不可恢
复的单向变化趋势。因此，通过对这个过程的统计分析，监测这个过程均值是否出现不可恢 复的单向变化趋势，即可评价桥梁结构的安全情况。
EWMA控制图的具体设计步骤如下
1) EWMA控制图的内控制限[Zr丄，(/C丄]的计算方法，包括
(1) 将桥梁结构的任一实时监测参数的历史信息建立一个随机过程，表述为
式中/为正整数；
Z为该监测信息的数据在周期内的第/个值； Z为该监测信息的数据的过程均值； 《为随机变化的参量；
(2) 确定(T、 d、和Ji "的值
对监测信息的数据进行分析，计算该数据的标准差<7以及最小偏移量d = ;/w (其中 I il 1^1);再根据桥梁结构的实际情况和评价的具体要求，选取数据变化的基本时间周期的
10倍~20倍，作为在某时间段内不发生误报警的平均运行时长ARL0(Average Run Length)。
ARL0的确定是根据桥梁结构的实际情况和评价的具体要求进行的。对桥梁这种大型结构 而言，由于其监测系统误发警报造成的成本相对较高，因此在结构处于正常状态时，希望在 较长时间内不发生误发警报的情况。以重庆马桑溪长江大桥的监测为例，结构测量信息是每 十分钟采集一次，对结构的某一测量点而言， 一天之内便有144数据点。为此, ARL0可选择 为2000点，即2周内不发生误发警报。
(3) 确定/l和k的值
通过査表法，从EWMA控制图最优/l曲线图和EWMA控制图、(义，K)组合曲线图确定;i 和A的值；
ARL0确定后，为了选择;i，需要先确定过程所能监测的最小偏移量J的大小。"的确定 需要综合考虑系统安全要求所允许的最大偏移量，最后将d的大小表示为正常观测样本标准 差(7的倍数。 一旦d的大小确定，即可根据它以及ARL0^的值，利用图2至图5插值确定X和 A:的值。
(4) 优化;i和k的值
根据桥梁结构监测的实际情况，建立一个迭代过程，通过"值的变化搜索最优的(/l力组 合，使被监测过程的误发警报的概率与预先确定的概率的误差最小；
因为^的值位于0-4之间，故最优( a)组合可根据桥梁结构处于正常状态的实际监测过 程的分析来确定。具体过程为根据桥梁结构响应监测过程的实际情况，建立一个迭代过程， 通过d值的变化搜索最优的( gt)组合，使被监测过程的误发警报的概率与预先确定的的概率 的误差最小。其搜索过程如图6，搜索所得的控制限记为[ZTZ， t/CZ]。
(5) 根据下式求得EWMA控制图的内控制限的上、下限
式中M为样本子集的大小，W=」i iVz^， A为采样间隔。 2) EWMA控制图的外控制限[rC乙//C丄]计算方法，包括
(1) 根据桥梁结构可靠度规范给出的可靠度指标，即可获得结构失效概率的允许值尸/;
(2) 根据下式计算EWMA控制图的外控制限系数[-x,x],(x>0);<formula>formula see original document page 17</formula>
(3)根据下式确定外控制限:<formula>formula see original document page 17</formula>
式中Z为监测信息的数据的过程均值； O"为监测信息的数据的标准差。
3)外控制限[FCX， //CZ]与内控制限[iTZ， t/CX]相结合形成"双控制限"，以此为判据 对桥梁结构进行评估并报警
内、外控制限都确定后，控制图即设计完毕，可用来监测桥梁各点结构响应的活载效应
和劣化效应信息的变化情况。 一旦控制图统计量的值超出内控制限[丄cz， c/cz]的范围，即表
示过程有失控的征兆，可进行初级预警；而一旦统计量的值超越内控制限[ZCX， f/CX]后，又 快速超越外控制限[KCi:， //C￡]，则反映了桥梁结构的安全状态打破了平衡，这时需要给出安 全告警信息。高度关注桥梁的运行情况，并采取进一步的应急措施。
上述方法在工程实际中进行了验证，评价结果如下
①瞬变信息控制图的评价结果
以下是对某座大桥部分测点挠度和应变瞬变信息用EWMA控制图进行评价的结果。对 跨中烧度测点n5和s5的瞬变信息建立控制图时，通过对其历史信息构成的过程的分析和迭 代计算，可得到最优的最小偏移量为^1.5cr，从而根据图2至图5确定A的值为0.21， A的 值为3.4，而根据桥梁结构可靠度规范计算所得到的[//<^:,^0:]为[-4.32， 4.32]。至此，用于 桥梁结构安全评价的控制图设计完毕，可将之用于实际应用。
为了模拟桥梁结构出现安全问题的情况，将挠度s5在7月1日后面的的信息加上一个持 续增大的随机变量，再用EWMA控制图进行评价。对于其余各点的控制限，通过计算也可 得到类似跨中测点的结果，只是由于正常状态的标准差o"不一样，而具体的控制限数值不一 样而已。
参见附图7至9，图中分别为跨中测点n5、 s5的挠度瞬变信息的控制图评价情况。从评 价结果可以看到，当结构处于正常状态时，挠度的瞬变信息控制图统计量的值均落在控制限 之间。图9为模拟测点s5在7月1日后挠度瞬变信息的波动幅度出现持续增大的情况评价， 可见当结构出现安全问题时，其控制图的统计量的波动幅度也出现持续的增大，并落在两级 控制限之外。
应变瞬变信息的控制图的设计方法和步骤与挠度信息的相似，只是由于原始信息不一样， 其波动幅度、标准差、最大允许偏移等参数也不同，从而使控制图的相应的控制限的宽度不 一样而已。
由于应变的影响因素繁多，而且影响方式极为复杂，而应变又是局部参量，在动荷载和 其它因素的的作用下，各个截面测点应变的变化情况不一样，但相对而言，在同一个截面内 的同一层面(上缘或下缘)应变的变化情况较为一致，故该大桥结构应变的控制图设计应根 据各个截面的具体情况进行，分别考虑各截面上下缘测点的布置情况，并按控制图的设计步 骤完成设计。
以下是对上述大桥主梁截面IV'的下缘测点20和O截面上缘测点30在6月25日开始到 7月2日中午的瞬变信息设计了相应的控制图，并进行安全评价的情况,参见附图10至12。 各控制图的X).18, )H3.4。
从图中看出，在同一截面的同一层面，两个测点的控制图的控制量变化相似，而截面不 同时，控制量的波动情况明显不一样。在结构处于正常状态时，控制量的值均落在控制限之 间。
同样，为了模拟结构出现安全问题的情况，在测点30在7月1闩下午20时30分丌始加 上一个持续增大的随机变量，其评价结果参见附图12。从图中可以看到，随着外加随机变量 的加入，控制量的值波动幅度不断增大，并迅速超出控制限。此时离出现设定的最小偏移量 的间隔点数只有6点。可见，控制图具有较好的性能。
②劣化效应控制图的评价结果
劣化效应信息的控制图设计与前述方法相似。以下是对跨中挠度测点s5和IV'截面应变 测点20的劣化效应信息的评价结果，参见附图13、 14。经搜索可得到控制图的内控制限的A: 值为[-3.1,3.1]，而X的值为0.25，外控制限的A:值仍然为[-4.32,4.32]。其中，横坐标是按效应 差中基准效应(减数)的采样时刻。A"m和AS^分别表示12小时内的劣化挠度差和应变差。 从图中可以看出，在结构处于正常状态时，统计量的值均落在内控制限[iXX,C/Ci:]内。
为进一歩提高本发明桥梁结构安全预警的ewma控制图方法安全预警的准确性，如结
构实时监测信息控制图统计量的值达到内控制限[zcz， c/o:]或外控制限[r:z， //cz]的80%
以上时，通过实时监测信息的变动趋势特性『及/指标和i s/指标进行二次或三次评估，并根 据评估结果判断是否发出初级预警或安全告警(所述变动趋势特性『及/指标和i s/指标实质 上是一个时间段内监测变量的单向变化程度，通过定量分析同样可以对桥梁结构的安全进行 评价)。
具体操作如下
如结构实时监测信息控制图统计量的值直接达到或超越内控制限[zcx， c/cz]或外控制限 [rcz， //c丄]时，直接发出发出初级预警或安全告警。
如结构实时监测信息控制图统计量的值达到内控制限[iXZ， C/CZ]或外控制限ITCi:， //GL] 的80%以上时，通过实时监测信息的变动趋势特性『i /指标和i S/指标进行二次评估，如果-25<WRI(n)<125，则不发出初级预警或安全告警，反之，则发出初级预警或安全告警。
如结构实时监测信息控制图统计量的值达到内控制限[zo:， t/cz]的80%以上时，通过实
时监测信息的变动趋势特性『R/指标和i s/指标进行三次评估，当评估结果同时满足条件
-25<WRI(n)<125和455及5/( )￡55时，则不发出初级预警；如果其中任一条件不满足，则 发出初级预警。
综上所述，在实际工程中，(一)、(二)和(三)中所述方法均可单独用于桥梁安全告警
19
系统，但如果采用(一)中所述方法，可以大大提高安全告警的准确性。实际工程中会出现 这样的情况根据(二)方法和(三)方法进行判断，两种方法都发出告警请求、两种方法 都不发出告警请求、其中一个发出告警请求另一个不发出；针对上述四种情况，根据(一) 方法仅可能作出两种判断，发出告警或不发出告警，这样就大大减小系统输出的不确定度， 提高了状态评估的准确性。
权利要求
1、一种桥梁结构安全的综合评价方法，其特征在于1)采用WRI指标和RSI指标对桥梁结构安全进行综合评价并形成评价结果；2)利用EWMA控制图的初级预警限和安全告警限对桥梁结构安全进行综合评价并形成评价结果；3)采用Dempster融合规则对WRI指标和RSI指标形成的评价结果和EWMA控制图的初级预警限和安全告警限形成的评价结果进行融合并形成融合结果；根据融合结果来确定最终是否发出告警信息，或者，直接根据步骤1)或步骤2)中的评价结果来确定最终是否发出告警信息。
2、 根据权利要求1所述的桥梁结构安全的综合评价方法，其特征在于融合的步骤包括1) 根据冲突因子《的值，确定是否存在评价结果冲突的情况如果没有冲突，则直接 用Dempster融合规则对评价结果进行融合；若有冲突，则进行下一步；2) 采用专家评分法或相关系数法计算平均评价结果的基本概率赋值5PJ;3) 根据融合的具体要求及评价结果的属性，计算各评价结果的权值；4) 用平均评价结果和没有冲突的评价结果建立新的评价结果源，并计算新评价结果的BiM，平均评价结果的权值由冲突评价结果的权值代替；5) 用Dempster融合规则对新评价结果进行融合，得到融合结果。
3、 根据权利要求2所述的桥梁结构安全的综合评价方法，其特征在于步骤2)中，采 用专家评分法或相关系数法计算平均评价结果的基本概率赋值万iM的步骤，包括采用专家评分法或相关系数法对评价结果进行预处理的步骤，包括(1)采用专家评分法构造融合评价变量的基本概率分布的方法，包括-根据下式确定专家评分法的基本概率赋值，即SiM:<formula>formula see original document page 2</formula>式中m (P表示由评价结果A确定的结构安全状态基本概率赋值^/M; m(//)表示由评价结果A不能确定的结构状态基本概率赋值5尸A w为桥梁结构某一单元的安全状态的种数；w为根据某参数的评价结果A对结构安全状态进行评分专家人数； C,G)为专家z'根据某参数的评价结果A判断结构安全状态为J'的基本概率赋值BiM; C,(//)为专家/根据某参数的评价结果A不能确定状态的基本概率赋值WP」； 其中，A为第i个监测变量的评价结果； (2)采用相关系数法构造融合监测变量的基本概率分布，包括: 根据下式确定相关系数的基本概率分布式中/、 /均为自然数；C,(/)为传感器/与目标状态7'的相关系数； w, (p为传感器/对目标状态/的基本概率分布； 其中尺，=W (1 — i , )(l —"力,) 式中《为系统总的不确定系数； M为传感器数量；",为传感器/与目标状态/的最大相关系数；A为传感器/与各相关系数的分布系数； 及,为传感器/的可靠系数； Yz为环境性能系数，其值位于
内； 各参数表示如下<formula>formula see original document page 3</formula>传感器/对辨识框架//的基本概率分布为<formula>formula see original document page 3</formula>式中附,(//)为传感器/对辨识框架//的基本概率分布。
4、根据权利要求2所述的桥梁结构安全的综合评价方法，其特征在于:步骤5)用Dempster 融合规则对新评价结果进行融合，得到融合结果，包括 辨识框架// = ^41,」2}，其中，^表示结构处于安全状态，A表示结构处于不安全状态，它们的基本概率分布分 别为和m(J2)，w(/z)表示不确定的状态， 基于基本概率赋值的决策对任意4，4G/Z,满足<formula>formula see original document page 4</formula>若有<formula>formula see original document page 4</formula>则不发出安全告警信息；不满足上述条件时，发出安全告警信息； S,、 ^为预先设定的阈值。
5、 根据权利要求1所述的桥梁结构安全的综合评价方法，其特征在于步骤2)中，所 述的利用EWMA控制图的初级预警限和安全告警限对桥梁结构安全进行综合评价并形成评价结果，包括EWMA控制图的内控制限[丄CZ， f/C丄]为初级预警限，外控制限[KC丄，州]为安全告警限，其评价结果为当桥梁结构的实时监测信息控制图统计量的值超过内控制限[ixx， c/o:]或外控制限[Fa:， //cz]时，发出告警请求，否则不发出告警请求。
6、 根据权利要求1所述的桥梁结构安全的综合评价方法，其特征在于步骤l)中，采 用『及/指标和i^/指标对桥梁结构安全进行综合评价并形成评价结果，包括1) i s/指标是指实时监测信息的数据在某一周期内的单向波动量占总波动量的相对百分 数，具体计算方法为-爿+ B式中A为实时监测信息的数据在一个周期内正方向增加量的总和；B为实时监测信息的数据在一个周期内负方向增加量绝对值的总和；其中<formula>formula see original document page 4</formula>式中Z,为实时监测信息的数据在一个周期内的第f个值；Z,爿为实时监测信息的数据在一个周期呢的第个值； /为实时监测信息的数据在一个周期内的取值数目； /为正整数。2) ^T /指标是指实时监测信息的数据当前波动幅度占其前一周期最大波动幅度的百分 数，具体计算方法为 證<formula>formula see original document page 5</formula>式中ZmM为实时监测信息的数据在前一周期内的最大值；Zmin为实时监测信息的数据在前一周期内的的最小值；Z为实时监测信息的数据的当前数值；n为实际测量数据结果的周期长度内的测点数量； 其中，实时监测信息的数据在前一周期内的最大值Z^和最小值Z皿分别为前一周期内 的测点数W减去距当前检测点W个测点后振幅的最大值和最小值，m取值可为3~6; 『i /指标和i S/指标形成的评价结果，包括1) - 25<『^ /(")<125时，不发出告警请求；反之则发出告警请求；或2) 同时满足-25<『7 /(")<125和455i S/(")555时，不发出告警请求；任一条件不满 足时，则发出告警请求；
7、 根据权利要求1或5或6所述的桥梁结构安全的综合评价方法，其特征在于 直接根据『及/指标和i S/指标形成的评价结果来确定最终是否发出告警信息，包括1)-25<『六/(")<125时，不发出告警信息；反之则直接发出告警信息；或2)同时满足-25<^7 /(")<125和45Si 57(")^55时，不发出告警信息；任一条件不满足时，则直接发出告警信息；直接根据EWMA控制图形成的评价结果来确定最终是否发出告警信息，包括当桥梁结构的实时监测信息控制图统计量的值超过内控制限[丄o:, t/o:]或外控制限[ra:， //ci]时， 直接发出告警信息，否则不发出告警信息。
8、 根据权利要求4所述的桥梁结构安全的综合评价方法，其特征在于EWMA控制图的内控制限[丄o:， [/ci]的计算方法，包括1) 将桥梁结构的任一实时监测参数的历史信息建立一个随机过程，表述为<formula>formula see original document page 5</formula>式中/为正整数；Z,为该监测信息的数据在周期内的第/个值； 'Z为该监测信息的数据的过程均值；《为随机变化的参量；2) 确定ff、 d、和^ ^的值对监测信息的数据进行分析，计算该数据的标准差ff以及最小偏移量d = (其中 hl^l);再根据桥梁结构的实际情况和评价的具体要求，选取数据变化的基本时间周期的10倍~20倍，作为在某时间段内不发生误报警的平均运行时长(Average Run Length)。3) 确定A和A:的值通过查表法，从EWMA控制图最优/l曲线图和EWMA控制图、(义，A)组合曲线图确定;i 和A的值；4) 优化义和A的值根据桥梁结构监测的实际情况，建立一个迭代过程，通过d值的变化搜索最优的(/i，;t)组合，使被监测过程的误发警报的概率与预先确定的概率的误差最小；5) 根据下式求得EWMA控制图的内控制限的上、下限<formula>formula see original document page 6</formula>式中W为样本子集的大小，"=^ """ A为采样间隔；ewma控制图的外控制限[ro:， i/C丄]计算方法，包括1) 根据桥梁结构可靠度规范给出的可靠度指标，即可获得结构失效概率的允许值尸/;2) 根据下式计算EWMA控制图的外控制限系数<formula>formula see original document page 6</formula>;3) 根据下式确定外控制限<formula>formula see original document page 6</formula>式中Z为监测信息的数据的过程均值； C7为监测信息的数据的标准差。
全文摘要
本发明公开了一种桥梁结构安全的综合评价方法，包括1)采用WRI指标和RSI指标形成评价结果；2)利用EWMA控制图形成评价结果；3)采用Dempster融合规则对WRI指标和RSI指标形成的评价结果和EWMA控制图的初级预警限和安全告警限形成的评价结果进行融合并形成融合结果；根据融合结果来确定最终是否发出告警信息，或者，直接根据步骤1)或步骤2)中的评价结果来确定最终是否发出告警信息；本发明的有益技术效果是大大减小系统输出的不确定度，提高状态评估的准确性，评价结果可靠度高、适用面广。
文档编号G01M99/00GK101382474SQ20081023282
公开日2009年3月11日 申请日期2008年10月8日 优先权日2008年10月8日
发明者梁宗保, 鹏 章, 胡顺仁, 伟 郑, 陈伟民, 黄晓微 申请人:重庆大学