1、第 卷 第 期 年 月人 民 长 江 ，收稿日期：基金项目：国家工信部智能制造综合标准化与新模式应用项目（）；河南省水利部水利科技攻关项目（）作者简介：杨 涛，男，硕士研究生，主要从事水工金属结构研究。：通信作者：赵华东，男，教授，博士，主要从事水工金属结构、智能制造研究。：文章编号：（）引用本文：杨涛，张钰奇，李凯旋，等 基于二元维纳过程的水工闸门剩余寿命预测 人民长江，（）：基于二元维纳过程的水工闸门剩余寿命预测杨涛，张 钰 奇，李 凯 旋，韩 明 海，赵 华 东（郑州大学 机械与动力工程学院，河南 郑州；河南省陆浑水库管理局，河南 洛阳）摘要：针对水工闸门运行过程中因性能退化而引起的剩余
2、寿命预测问题，提出了一种基于二元维纳过程的水工闸门剩余寿命预测方法。利用二元维纳过程建立闸门因腐蚀而引起的变形、应力退化过程模型，并基于 准则选择出合适的 函数分析其相关特性；通过分步极大似然估计法估计出模型中的参数，得到闸门剩余寿命概率密度函数，同时基于剩余寿命预测结果确定闸门的维修时间。最后，以河南省陆浑水库的弧形钢闸门为例进行数值仿真，验证了结果的合理性。结果表明：采用多指标方法更能全面、准确地反映出闸门的健康状态，基于二元维纳过程预测的闸门剩余寿命其 为 ，小于一元维纳过程的 ，具有较高的预测精度。研究成果可为合理评估闸门的寿命状态，确定合适的维修时间以延长闸门寿命提供依据。关 键 词
3、：闸门腐蚀；剩余寿命预测；二元维纳过程；函数中图法分类号：文献标志码：引 言水工闸门随着运行年限的增加，会出现腐蚀、磨损等情况，从而严重影响闸门的运行安全。一旦闸门出现破坏，会带来巨大的经济损失，甚至会造成人员伤亡。目前，随着数字化、信息化技术的发展，对闸门进行健康状态评估越来越重要。剩余寿命预测作为健康监测的核心部分，对于指导闸门进行预测性维护以保证闸门安全运行具有重要意义。近年来，闸门寿命预测方法越来越受到学者们的重视。赵林章等利用分形维数法对闸门进行了寿命预测；等计算了闸门侧轮和轨道的疲劳寿命；徐卫敏等利用 随机过程对钢板锈蚀量进行随机模拟，并采用 法对闸门上某构件进行了时变可靠度计算；
4、李伟康等考虑到闸门锈蚀对材料性能会产生影响，基于 曲线推导出了闸门构件寿命预测的公式。但以上研究还存在着一些不足之处，主要体现在寿命预测的精度不高，用某个构件代替整个闸门，缺乏合理性。基于随机过程的寿命预测方法可以很好地描述装备的不确定性，现已逐渐被应用于长寿命、高可靠度的装备寿命预测。等基于逆高斯过程对刀具寿命进行预测，达到了较高的精度；等基于 过程对锂电池剩余寿命进行了预测；等利用维纳过程建立了状态空间模型，并对轴承进行了寿命预测；等基于维纳过程对复杂机电系统进行了寿命预测。但是可以发现，以往研究未充分考虑闸门在运行过程中的复杂工况，且单个指标难以准确反映出其健康状态。函数可以合理刻画出不
5、同性能指标间的相关特性。等运用 函数表征了腐蚀管道不同缺陷间的相关关系；等利用 函数分析了轴承多种失效模式的相关性，并进行了寿命预测。因此本文采用多个指标进行寿命预测。第 期 杨 涛，等：基于二元维纳过程的水工闸门剩余寿命预测基于上述分析，本文选取闸门在腐蚀情况下变形和应力两个性能指标来表征其健康状态。首先，基于维纳过程对闸门变形和应力进行退化性建模；其次利用 函数分析其相关特性，得到闸门剩余寿命概率密度函数；通过分步极大似然估计法对维纳过程和 函数中的参数进行估计；最后对比分析了基于二元维纳过程和只考虑应力影响下的一元维纳过程进行寿命预测的精度，同时利用剩余寿命预测结果指导闸门确定了合理的维
6、修时间。基于二元维纳过程的退化建模 维纳过程当采用维纳过程进行闸门退化建模时，第 个性能指标的退化过程可表示为（）（）（）（）（）（）（）（）式中：（）（）为 时刻第 个指标退化值；（）（）为退化初始值；（）为漂移系数；（）为扩散系数；（）为标准维纳过程。基于 函数的相关性分析由于闸门的变形和应力两个指标之间不是相互独立，也不是简单的线性关系，所以需要选择出合理的方法刻画不同指标间的关系。函数能有效地将联合分布函数与边缘分布函数联合起来，常用于分析不同指标之间的相关特性。因此，本文利用 函数来分析闸门上多元性能指标之间的相关性，如式（）所示：（，）（），（），（）（，；）（）式中：（，）为各个
7、变量联合分布函数；（），（），（）为各个变量的边缘分布函数；为 函数；为 中的参数。函数常见的 种类型如表 所列，其中，表示边缘分布函数。表 常用的 函数 函数（，；）（）（）（）（）（）()（）（），由 定理得，闸门变形和应力的相关特性可表示为（，）（），（）；（）由于不同的 函数对于描述变形和应力两个退化指标的相关特性具有较大影响，所以应结合闸门实际退化情况选择合适的 函数。信息准则鼓励数据拟合的优良性，被广泛运用于评价模型拟合数据效果的优劣。值越小，则说明拟合的效果越好。本文采用 准则选择 函数：（）式中：为模型对应的似然函数；为模型中参数的个数。剩余寿命预测模型 剩余寿命分布当闸门两个
8、性能指标均服从维纳过程时，其剩余寿命服从逆高斯分布，剩余寿命的概率密度函数为（，）（）（）（）（）式中：为第 个性能指标的失效阈值；为第 个性能指标在 时刻的性能退化值。当获得闸门上变形和应力各自的剩余寿命概率密度函数后，由 函数可以得到应力和变形剩余寿命分布函数的联合概率密度函数：（）（），（）；）（，）（）式中：（）、（）为变形和应力剩余寿命边缘分布函数；（），（）；）为（），（）；）的概率密度函数。定义变形和应力中任一性能指标超过失效阈值，便认为装备失效。则闸门的剩余寿命 （，）。因此，闸门剩余寿命的分布函数为（，）（）（）（），（）（）闸门剩余寿命的概率密度函数为（）（，）（，）（），
9、（）；）（，）（）参数估计由于要估计的参数较多，分步极大似然估计法具有计算简单、流程清晰等优点，所以采用分步极大似然估计法更新模型中的参数。（）由维纳过程的性质可知，增量服从正态分布：人 民 长 江 年（）（）（，（）（）式中：为第 个性能指标退化量增量。所以参数（，）的极大似然函数为（，）（）（）（）（）对，求偏导，可得，的极大似然估计值：（）（）（）（）将变形和应力的剩余寿命边缘分布函数（，）代入 函数中，利用极大似然估计法求出 函数中参数 的估计值：（，），（，）；）（）工程实例分析陆浑水库溢洪道露顶式弧形钢闸门构件采用 钢，闸门面板中心弧面半径为 ，面板、纵梁厚度为 ，支臂厚度为 ，主
10、横梁厚度为，作用水头为 。闸门变形和应力作为闸门上重要的性能指标，在运行过程中会发生退化，可反映出闸门的健康状态状况。变形会影响闸门安全运行，本文定义变形阈值为 。根据 水利水电工程钢闸门设计规范，容许应力与闸门运行情况和钢材的厚度有关。该弧形闸门最大应力处构件厚度均不大于 ，属于钢材尺寸分组中的第一组，同时取修正系数 ，所以定义闸门上最大应力阈值为 。本文采用文献中所论述的 锈蚀深度非线性变化公式，如式（）所示，并根据文献中的统计数据确定参数，计算时段取 ，同时利用有限元仿真得到闸门不同年份上的变形和应力变化。（），（），（）式中：为腐蚀量的长期厚度，全面腐蚀取 ，局部腐蚀取 ；（）为 时刻
11、的腐蚀量厚度；为涂层寿命，取 ；为过渡时间，取 。由于 中的 模块有着简单的格式，可以方便对有限元模型进行计算，所以利用 对弧形钢闸门进行参数化建模和求解。由于闸门属于空间薄壁结构，因此采用 单元进行建模。闸门上锈蚀主要分为全面锈蚀和局部锈蚀：对于全面锈蚀，采用平均蚀余厚度法进行模拟；对于局部锈蚀，通过改变锈蚀坑所在位置节点的厚度来实现对锈蚀坑的模拟。根据闸门实际运行情况，局部锈蚀一般情况下出现在闸门的面板右下角漏水的部位、右下角边梁处以及下主梁经常性积水区域，所以通过改变上述区域的节点厚度实现对局部锈蚀坑的模拟。闸门的局部锈蚀区域如图 所示。图 闸门局部锈蚀区域 利用 对闸门进行有限元分析，
12、可以得到闸门的变形和等效应力。闸门初始状态的变形和等效应力云图如图 所示。图 闸门变形云图 图 闸门应力云图 根据腐蚀非线性规律改变闸门对应部分的厚度，同时利用 仿真可以得出闸门在退化过程中不同运行年限的变形和等效应力，将不同运行年限的闸门 第 期 杨 涛，等：基于二元维纳过程的水工闸门剩余寿命预测变形和应力值连接成线，可以得到闸门上变形和应力退化曲线，如图 所示。图 闸门性能指标退化曲线 为了验证使用维纳过程模型对闸门因腐蚀而产生的退化过程建模是否合理，需要对变形和应力退化数据进行检验。由于 检验法能够利用样本数据推断样本的总体是否服从某一理论分布，因此利用 检验法对闸门退化数据进行检验，检
13、验结果如表 所列。其中：表示在置信水平为 时接受退化量增量为正态分布的原假设，时表示拒绝原假设；当 时，表示不拒绝退化增量为正态分布的原假设。通过分析可以看出闸门上的变形和应力退化量均满足正态分布假设，故可以采用维纳过程模型来描述闸门退化过程。表 检验结果 退化性能指标变形 应力 通过式（）得到变形和应力的剩余寿命边缘概率密度函数后，利用 准则选择合适的 函数，不同 函数得到的 值如表 所列。由表 可以看出：函数的 值最小，所以选择 函数来分析变形和应力的相关特性。由 函数得到变形和应力的剩余寿命联合概率密度函数图如图 所示。由图 可以看出 函数具有不对称的尾部结构，上尾高，下尾低，这说明在分
14、布的上尾部变形和应力之间具有较强的相关性，在分布的下尾部两者之间则为渐进独立的。表 种 函数的 值 函数 值图 二元 概率密度函数 在选择好合适的 函数之后，由式（）（）可得到不同运行年限的参数值，如表 所列。之后结合式（）可以得到不同时刻的闸门剩余寿命概率密度函数，同时计算出基于应力退化的一元维纳过程预测的剩余寿命概率密度函数，将闸门概率密度函数最大值对应的时刻作为剩余寿命的预测值，如图 所示。由图 可以看出，随着闸门运行年限的增加，数据的增多，剩余寿命预测结果的分布范围越来越小，说明其预测结果的精确度越来越高，闸门发生失效的概率越来越高，闸门寿命减少，符合工程实际。表 不同运行年限的参数估
15、计值 运行年限（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）（，）为了能够清晰地看出基于一元维纳过程和二元维纳过程预测闸门剩余寿命的变化趋势，将闸门实际寿命取 水利水电金属结构报废标准中所规定的大型闸门折旧年限 ，将两种方法预测的闸门剩余寿命与真实剩余寿命相比较，如图 所示。从图 可以看出随着闸门运行年限的增加，基于二元 人 民 长 江 年维纳过程的闸门剩余寿命预测与真实寿命差距较小。图 闸门剩余寿命概率密度 图 不同方法预测寿命变化趋势 为了客观比较两种方法预测闸门剩余寿命的精度，引入均方根误差（）和相对误差指标进行评价。定义相对误差指标：（）式中
16、：为闸门剩余寿命预测值；为闸门实际剩余寿命。一元维纳过程和二元维纳过程预测结果的 分别为 和 ，这说明相较于单一指标，两个指标更能全面准确反映出闸门的寿命状态。不同预测方法的相对误差 如表 所列，由表 可以看出二元维纳过程的 占比达到 ，而一元维纳过程 仅占比 ，同时，基于二元维纳过程的，占比为，这说明基于二元维纳过程进行闸门寿命预测的方法可信度较高。表 不同预测方法的相对误差指标 预测方法 一元维纳过程二元维纳过程 通过对闸门剩余寿命预测可以指导闸门确定合理的维修时间，既避免了经常性维修造成的资金浪费，又避免了闸门因过久不维修而发生破坏，造成巨大的安全事故。基于随机过程预测闸门维修时间，可以根据闸门的重要性程度和维修的费用，确定闸门上的维修阈值。本文假设闸门上的变形维修阈值为 ，应力维修阈值为 ，当闸门到达维修阈值后，便进行一次全面大修，以延长闸门的剩余使用寿命。由于数据越多，预测越精确，因此根据 的数据，将维修阈值代入式（）（）可以得到闸门的维修时间，如图 所示。由图 可知，维修阈值下剩余寿命还剩 ，说明至少应在闸门运行 的时候进行一次大修，以确保闸门的运行安全。图 闸门维修时间预