基于历史故障信息和运行状态的电力电缆运行维护方法
基于历史故障信息和运行状态的电力电缆运行维护方法
【技术领域】
[0001] 本发明涉及一种基于历史故障信息和运行状态的电力电缆运行维护方法，其属于 电力设备资产管理领域，适用于预测电力电缆的故障率，并通过提前制定更换策略来更换 一些运行年份较长的电缆，减缓故障率的增长趋势，提高电缆线路的可靠性，从而提高电网 运行的稳定性。
【背景技术】
[0002] 国内外主要通过建立统计模型分析电缆的故障模式，得出电缆整体的故障概率密 度分布，预测未来的电缆故障数；同时结合电缆的状态以及电缆在系统中的地位，制定电缆 的维修更换计划。通常，电缆的故障分为早期、中期和晚期三个阶段，在进行故障分布的拟 合时，Weibull分布和Crow-AMSAA模型（可靠性增长模型）能够很好的反映这三个阶段，因 此在可靠性领域使用较多，但Weibul1分布和Crow-AMSAA模型存在其局限性。使用Weibul1 分布的不足：
[0003] l)Weibull分布分析的是故障前运行时间，对故障数据的质量要求相对较高，需要 知道设备详细的投运日期和故障日期，当故障前运行时间不满足Weibull分布时，无法使 用Weibul1分布进行分析；
[0004] 2)Weibull分布通常假设电缆为一种不可修复的元件，这显然与电力电缆的实际 运行情况不相符；
[0005] 3)Weibull分布在国外的研宄工作通常都基于大量的老化故障数据，但国内的 电缆故障主要为早期故障，由于早期故障数据的特点与老化故障数据不同，不能简单套用 Weibull分布模型。
[0006] 使用Crow-AMSAA模型的不足：
[0007] 1)Crow-AMSAA模型适用于分析多种混合故障模式的故障数据，对于分析单一故障 模式的故障数据适用性差；
[0008] 2)当分析早期电缆故障数据特点时，由于故障数据样本少，Crow-AMSAA模型抗扰 动性差，分析结果可靠性较差；
[0009] 此外，Weibull分布和Crow-AMSAA模型仅基于电缆历史故障信息来预测电力电缆 故障率也存在明显不足。实际上，预测电缆故障率不仅与电缆历史故障信息有关，还与电缆 实际运行状态有关。电缆运行状态随负荷水平、安装方法以及安装位置的不同而异，从而导 致电缆故障率随运行状态而变化。且进行电缆故障预测的目的就是了解电力电缆运行情 况，并通过提前制定维修更换电缆策略，减缓故障率的增长趋势，提高电缆线路的可靠性。 因此，提出一种能够同时考虑电缆历史故障信息和运行状态的电力电缆运维策略，具有十 分重要的现实意义。
【发明内容】
[0010] 本发明的目的是提供一种基于历史故障信息和运行状态来的电力电缆运维方法， 用来预测电缆剩余寿命，同时提出电缆更换策略，提高整体电缆的可靠性。
[0011] 为达到上述目的，本发明采用的技术方案是：
[0012] 一种基于历史故障信息和运行状态的电力电缆运行维护方法，用于预测电力电缆 的剩余寿命，以及制定电力电缆的更换策略，该方法包括以下步骤：
[0013] (l)weibull比例风险模型的建立：
[0014] 所述的weibull比例风险模型基于Cox比例风险模型和weibull分布而建立，当 所述的Cox比例风险模型的基础风险函数服从所述的weibull分布时，所述的Cox比例风 险模型的函数表达式用作分析所述的电力电缆的故障数据和状态量的weibull比例风险 丰旲型；
[0015] (2)weibull比例风险模型中参数的求解：
[0016] 选取所述的weibull比例风险模型中各个参数的初值，从而求解所述的weibull 比例风险模型中的各个参数；
[0017] (3)电力电缆状态决策曲线的求取：
[0018] 将求解得的所述的weibull比例风险模型中的各个参数代入所述的weibull比例 风险模型，得出所述的weibull比例风险模型的表达式，再根据所述的weibull比例风险模 型的表达式绘制所述的电力电缆的状态决策曲线；
[0019] (4)电力电缆剩余寿命的预测：
[0020] 获取所述的电力电缆的历史状态数据并拟合得到所述的电力电缆的状态信息随 时间变化的关系曲线，获取所述的关系曲线与所述的状态决策曲线的相交点，则所述的电 力电缆当前运行点与所述的相交点之间的时间差即为所述的电力电缆的剩余寿命；
[0021] (5)电力电缆维修更换策略的制定：
[0022] 所述的电力电缆维修更换策略包括最小成本法和最大可用度法：
[0023] ①最小成本法：建立所述的weibull比例风险模型下所述的电力电缆发生故障而 进行更换时的成本率模型，并求解所述的成本率最小时的时间，从而得到所述的电力电缆 更换成本率最小的更换周期；
[0024] ②最大可用度法：建立所述的weibull比例风险模型下所述的电力电缆的可用度 模型，并求解所述的可用度最大时的时间，从而得到所述的电力电缆可用度最大的预防性 维护周期。
[0025] 所述的步骤（1)中，对不同类型的所述的状态量分别进行分析。
[0026] 所述的步骤（2)中，利用迭代法求解所述的weibull比例风险模型。
[0027] 所述的步骤（3)中，根据所述的weibull比例风险模型的表达式获得所述的电力 电缆的可靠度，使所述的电力电缆的可靠度大于满足要求的最小可靠度并分别取自然对 数，以时间为横坐标、以取自然对数后的电力电缆的可靠度为纵坐标，从而绘制所述的电力 电缆的状态决策曲线。
[0028] 所述的步骤（3)中，所述的电力电缆的可靠度取两个不同值分别绘制所述的电力 电缆的状态决策曲线。
[0029] 所述的电力电缆运行维护方法通过电力电缆状态监测和智能评估系统自动实现。
[0030] 由于上述技术方案运用，本发明与现有技术相比具有下列优点：
[0031] l)Weibull比例风险模型同时分析电缆历史故障数据和电缆运行状态（如护层环 流、最大局部放电量变化率与介损变化率），进行电缆故障预测，更加接近电缆的实际情况， 提高了预测的合理性、可靠性；
[0032] 2)根据Weibull比例风险模型可绘制状态决策曲线，直观明了的显示出电缆的运 行情况，可为电缆是否需要维修更换提供参考；
[0033] 3)如果已知某电缆的历史状态信息，Weibull比例风险模型通过拟合可得到电缆 的状态信息随时间变化的关系曲线，预测电缆的剩余寿命；
[0034] 4)Weibull比例风险模型的基础上可提出最小成本法和最大可用度法，结合考虑 经济因素，制定电缆维修更换策略。
【附图说明】
[0035] 附图1为本发明的电力电缆运行维护方法的流程示意图。
[0036] 附图2为电力电缆的Weibull比例风险模型h(t，X)示意图。
[0037] 附图3为电力电缆的状态决策曲线示意图。
[0038] 附图4为电力电缆的状态信息随时间变化的关系曲线与状态决策曲线的相交点 示意图。
[0039] 附图5为电力电缆更换成本率最小决策图。
[0040] 附图6为电力电缆预防性维修最佳时间决策图。
【具体实施方式】
[0041] 下面结合实施例对本发明作进一步描述。
[0042] 实施例一：一种用于预测电力电缆的剩余寿命，以及制定电力电缆的更换策略的 基于历史故障信息和运行状态的电力电缆运行维护方法，如附图1所示，其在电力电缆状 态监测和智能评估系统中通过weibull比例风险模型的建立、weibull比例风险模型中参 数的求解、电力电缆状态决策曲线的求取、电力电缆剩余寿命的预测、电力电缆维修更换策 略的制定这几个步骤实现。
[0043] (1)weibull比例风险模型的建立
[0044]weibull比例风险模型基于Cox比例风险模型和weibull分布而建立，当Cox比例 风险模型的基础风险函数服从weibull分布时，Cox比例风险模型的函数表达式用作分析 电力电缆的故障数据和状态量的 weibull比例风险模型。
[0045]Cox比例风险模型的函数表达式如式（1)
[0047] 其中，X#依时协变量，对应为电力电缆的运行状态信息，Y对应的回归系 数，n为依时协变量的个数，0为Weibull分布形状参数，t为时间，n为Weibull尺度参 数。此时模型称为Weibull比例风险模型，其能够同时分析故障数据和在线监测数据。
[0048] 采用上述Weibull比例风险模型分析电力电缆时，对不同类型的状态量分别进行 分析。在应用Weibull比例风险模型前，首先对纳入模型的状态量进行选择。比如，局部放 电和介损反映的是电缆绝缘内部缺陷或老化的状态，而护层环流反映的是绝缘外部的缺陷 或故障，因此在对电缆的状态进行评估时，将局部放电和介损归为一类，将护层环流归为另 一类，分别进行分析。由于局放特征量与电缆绝缘材料、绝缘缺陷形式、大小、电缆所受应力 等均有关系，局放特征量与绝缘老化程度的关联度尚不清晰，同时介损的测量也受到环境 因素以及测量误差的影响，因此将这两个特征量纳入模型进行分析时，以单位时间内的介 损的变化率和单位时间内的局部放电能量作为变量。
[0049] (2)weibull比例风险模型中参数的求解
[0050] 求解方法为：选取weibull比例风险模型中各个参数的初值，从而利用迭代法求 解weibull比例风险模型中的各个参数。weibull比例风险模型中参数的求解包括求解含 一个协变量的weibull比例风险模型和求解两个协变量的weibull比例风险模型。
[0051] ①求解含一个协变量的weibull比例风险模型
[0052] 含有一个协变量的weibull比例风险模型的函数表达式如式（2)
[0054] 其中，X为依时协变量，a为X对应的回归系数。
[0055] 可靠度表示电力电缆在时间t内正常工作的概率，Weibull比例风险模型对应的 可靠度函数为：
[0057] 故障概率密度函数为
[0058]f(t,X) =h(t,X) ?R(t,X) (4)
[0059] 构造故障概率密度函数的似然函数：
[0061] 其中n为数据的总数，m为故障数据个数，n-m表示截尾数据个数。
[0062] 两边取自然对数得
[0064] 为了求解待求参数为0，n和a，对上式分别求偏导数和二阶偏导数
[0066]
[0076]
[0077] 则对数似然函数的二阶求导矩阵
[0079]然后通过Newton-Raphson迭代，直至收敛，贝丨J :
[0081] 当k = 0时，代入0，n和a的初值0。，qQ和a。进行求解。
[0082] ②求解含两个协变量的weibull比例风险模型
[0083] 当weibull比例风险模型中有两个协变量时，求解模型参数时，需对上述求解方 法重新构造迭代式子。其中两个协变量时，Weibull比例风险模型函数为：
[0085] 故障概率密度函数为
[0086] f (t, X) = h(t, X) ? R(t, X) (22)
[0087] 构造故障概率密度函数的似然函数
[0089] 其中n为数据的总数，m为故障数据个数，n-m表示截尾数据个数。
[0090] 两边取自然对数，得：
[0092] 为了求解待求参数为0，n，a种a2，对上式分别求0，q，a种a2的偏导数 和二阶偏导数。
[0093]
[0103]
[0113] 则对数似然函数的二阶求导矩阵：
[0114]
[0115] 然后通过Newton-Raphson迭代，贝丨J:
[0117] 当k= 0时，代入0，q，a兩a2的初值0 〇，q〇，a1〇和a2〇进行求解。
[0118] 现阶段由于受条件限制，目前很难获得电缆的全部状态数据，因此对状态数据进 行了一些假设，如表1所示，其中电缆状态为1表示发生过故障，电缆状态为0表示未发生 故障。对于故障电缆，假设一些整体处于晚期故障模式的数据。
[0119] 表1电缆状态数据
[0121]
[0122] (3)电力电缆状态决策曲线的求取
[0123] 将迭代求解得的weibull比例风险模型中的各个参数代入weibull比例风险模 型，得出weibull比例风险模型的表达式；再根据weibull比例风险模型的表达式绘制电力 电缆的状态决策曲线。绘制电力电缆的状态决策曲线时，根据weibull比例风险模型的表 达式获得电力电缆的可靠度，使电力电缆的可靠度大于满足要求的最小可靠度并分别取自 然对数，以时间为横坐标、以取自然对数后的电力电缆的可靠度为纵坐标，从而绘制电力电 缆的状态决策曲线。
[0124] 以上述含两个协变量的weibull比例风险模型为例，选取初值|3 Q= 1. 8，n。= 10000,a1(|= 〇.〇〇1，a2。= 〇.〇9,代入程序进行迭代，得到 0 = 3.246,n= 114252,ai =0. 078,a2= 〇. 022,Weibull比例风险模型的函数形式如下：
[0126] 则电力电缆的Weibull比例风险模型h(t，X)如附图2所示。
[0127] 假如电缆系统所要求的最小可靠度为心，为了保证电缆的安全运行，电缆的可靠 度必须大于&，则有：
[0129] 两边分别取自然对数有
[0130] aj ?Xi<ln(-lnR〇) - 0 (lnt-lnn) (49)
[0131] 以t为横坐标，ln(-lnRQ)-f3 (lnt-lnq)为纵坐标绘制曲线即为决策曲线。
[0132] 可以对电力电缆的可靠度取两个不同值分别绘制电力电缆的状态决策曲线，如分 别绘制系统最小可靠度为0. 95和0. 99时的决策曲线，如附图3所示，代入某时刻t电力电 缆的状态量X，计算a*X，如果点（t，a 落在虚线上方，表明此时需要对电力电缆进行 维修更换；如果点（t，a 落在实线下方，表明电力电缆仍可继续工作；如果点（t，a*X) 落在实线与虚线之间，这部分区域为临界区，表明需要对电力电缆状态进行密切关注，重点 监测。
[0133] (4)电力电缆剩余寿命的预测
[0134] 获取电力电缆的历史状态数据并通过数学拟合得到电力电缆的状态信息随时间 变化的关系曲线，获取关系曲线与状态决策曲线的相交点，则电力电缆当前运行点与相交 点之间的时间差即为电力电缆的剩余寿命，如附图4所示。假设电力电缆的状态信息随时 间变化的关系曲线（图中圆圈表示）与状态决策曲线相交于点B，计算点A(当前运行点） 和B的时间差h即为电缆的剩余寿命。如果电缆的状态量从点A开始变化缓慢，则该关系 曲线（图中三角表示）与决策曲线相交于点C，此时的剩余寿命&大于h。为了提高电缆 寿命预测的精度，需要对电缆的状态量进行连续监测，收集足够多的状态信息，从而实现对 电缆寿命的预测。
[0135] (5)电力电缆维修更换策略的制定：
[0136] 电力电缆维修更换策略包括最小成本法和最大可用度法。
[0137] ①最小成本法：建立weibull比例风险模型下电力电缆发生故障而进行更换时的 成本率模型，并求解成本率最小时的时间，从而得到电力电缆更换成本率最小的更换周期。
[0138] 电缆发生故障后，目前一般采取更换措施。假设购买电缆的费用为，电缆故障后 的损失为C2，电缆的更换周期为t，即电缆运行时间达到t时，即使电缆没有发生故障也对 电缆进行更换。则电缆更换的成本率为：
[0140] 对成本率求导可得：
[0142] 其中
h(t)为故障率函数，f(t)为故障概率密 度函数，F(t)为累积故障概率密度函数，R(u)为残存函数，X表示状态量，则当
[0144] 时成本率最小，求出此时的t即为电缆更换成本率最小的更换时间。
[0145] 对上式进一步变换可得
[0147]根据已求出的h(t，X)的函数表达式，可以求出F(t，X)和R(t，X) 。由于
积分困难，首先在x轴（0~5000)范围内绘制R(u，X)的图形，然后通过Origin 软件对该图形进行数值积分。在进行数值积分时，Origin对x轴（0~5000)范围内均勾 划分为99个区间Xl= [50. 5,101，…，5000]，求得数值积分的值yi= [50. 5,100. 9，…， 2333. 2]。然后以Xl为横坐标，y 纵坐标，绘制并拟合曲线，该拟合的曲线方程即为
的积分表达式，以上式中的参数为例时
[0149] 假设(VC2= 1，当X:= 150,X:= 30时，分别绘制h(t,X)和
这两条 曲线，如图5所示，两者相交于点A，点A对应的时间为2018，表明当电缆的状态量分别为Xi=150,X2= 30时，最佳的更换时间为电缆运行的第2018天。
[0150] 如表2和3所示，讨论了当协变量值不变，（VC2比值发生变化时对最佳更换时间 的影响，发现随着(Vc2比值的增大，最佳更换时间也在增大，即如果故障损失与购买电缆 成本相比很小，那么电缆最佳更换时间会增大，电缆可用时间增大；当Q/^比值不变，发现 随着协变量&值的增大，最佳更换时间在减小，这表明当表征电缆的状态变化较快时，应尽 快更换电缆。
[0151] 表2 (VC2比值变化时对最佳更换时间的影响（X150,X2= 30)
[0153] 表3 值变化时对最佳更换时间的影响（X2= 30,C^4= 1)
[0154]
[0155] ②最大可用度法：建立weibull比例风险模型下电力电缆的可用度模型，并求解 可用度最大时的时间，从而得到电力电缆可用度最大的预防性维护周期。
[0156] 考察可用度时，以一个预防性维修周期为例进行说明。假设在每一个预防性维修 周期内，电缆的可用时间为Ta，电缆的不可用时间为Tb，电缆的预防性维修的时间为，故 障后维修时间为T2。因此在每一个预防性维修周期内电缆不可用时间为
[0157]Tb=R(t,X) ?Ti+d-Ra.X)) ?T2 (55)
[0158] 根据可用度的定义
[0160]则Weibull比例风险模型下的可用度为
[0162] 令
为使可用度最大，贝U只需B最小，由于积分
不易求得，另1(1天），T2= 2，X1= 150，X2= 30,绘制曲线B，如附图6所 示，得到A点为曲线的最低点，进而得出A点对应的x轴值为1832,表明使可用度最大化的 预防性维修的时间间隔为1832天。
[0163] 如表4和5所示，讨论了当协变量值和T2值不变，Ti值发生变化时对最佳预防性 维修时间的影响，发现随着1\值的增大，最佳预防性维修时间也在增大；当T:和T2值以及 &值不变时，发现随着协变量值的增大，最佳预防性维修时间在减小，这表明当表征电缆的 状态变化较快时，应尽快对电缆进行维修。
[0164] 表4 值变化时对最佳预防性维修时间的影响（X150,X2= 30,T2= 2)
[0165]
[0166]表5Xi值变化时对最佳预防性维修时间的影响（X2= 30,Tf1，T2= 2)
[0168] 本发明的基于历史故障信息和运行状态的电力电缆运维策略，基于Cox比例风险 模型和Weibull分布，建立Weibull比例风险模型。Weibull比例风险模型能够同时分析故 障数据并在线监测状态数据，从而能够根据Weibull比例风险模型对电缆历史故障信息和 运行状态的分析，预测电缆剩余寿命，提出最小成本法和最大可用度法，用来指导制定电力 电缆的更换策略。
[0169] 上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人 士能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明 精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种基于历史故障信息和运行状态的电力电缆运行维护方法，用于预测电力电缆的 剩余寿命，以及制定电力电缆的更换策略，其特征在于：该方法包括以下步骤： (Dweibull比例风险模型的建立： 所述的weibull比例风险模型基于Cox比例风险模型和weibull分布而建立，当所述 的Cox比例风险模型的基础风险函数服从所述的weibull分布时，所述的Cox比例风险模 型的函数表达式用作分析所述的电力电缆的故障数据和状态量的weibull比例风险模型； (2) weibull比例风险模型中参数的求解： 选取所述的weibull比例风险模型中各个参数的初值，从而求解所述的weibull比例 风险模型中的各个参数； (3) 电力电缆状态决策曲线的求取： 将求解得的所述的weibull比例风险模型中的各个参数代入所述的weibull比例风险 模型，得出所述的weibull比例风险模型的表达式，再根据所述的weibull比例风险模型的 表达式绘制所述的电力电缆的状态决策曲线； (4) 电力电缆剩余寿命的预测： 获取所述的电力电缆的历史状态数据并拟合得到所述的电力电缆的状态信息随时间 变化的关系曲线，获取所述的关系曲线与所述的状态决策曲线的相交点，则所述的电力电 缆当前运行点与所述的相交点之间的时间差即为所述的电力电缆的剩余寿命； (5) 电力电缆维修更换策略的制定： 所述的电力电缆维修更换策略包括最小成本法和最大可用度法： ① 最小成本法：建立所述的weibull比例风险模型下所述的电力电缆发生故障而进行 更换时的成本率模型，并求解所述的成本率最小时的时间，从而得到所述的电力电缆更换 成本率最小的更换周期； ② 最大可用度法：建立所述的weibull比例风险模型下所述的电力电缆的可用度模 型，并求解所述的可用度最大时的时间，从而得到所述的电力电缆可用度最大的预防性维 护周期。2. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（1)中，对不同类型的所述的状态量分别进行分析。3. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（2)中，利用迭代法求解所述的weibull比例风险模型。4. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（3)中，根据所述的weibull比例风险模型的表达式获得所述的电 力电缆的可靠度，使所述的电力电缆的可靠度大于满足要求的最小可靠度并分别取自然对 数，以时间为横坐标、以取自然对数后的电力电缆的可靠度为纵坐标，从而绘制所述的电力 电缆的状态决策曲线。5. 根据权利要求4所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（3)中，所述的电力电缆的可靠度取两个不同值分别绘制所述的电 力电缆的状态决策曲线。6. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：其通过电力电缆状态监测和智能评估系统自动实现。
【专利摘要】本发明涉及一种用于预测电力电缆的剩余寿命，以及制定电力电缆的更换策略的基于历史故障信息和运行状态的电力电缆运行维护方法，其通过weibull比例风险模型的建立、weibull比例风险模型中参数的求解、电力电缆状态决策曲线的求取、电力电缆剩余寿命的预测、电力电缆维修更换策略的制定这几个步骤实现。该方法提高了电缆故障预测的合理性、可靠性，直观明了的显示出电缆的运行情况，可为电缆是否需要维修更换提供参考，在Weibull比例风险模型的基础上可提出最小成本法和最大可用度法，结合考虑了经济因素，能够制定电缆维修更换的合适策略。
【IPC分类】G06Q50/06, G06Q10/06
【公开号】CN104899690
【申请号】CN201510308384
【发明人】张梁, 张俊, 郑勇, 黄晓军, 周承科, 周文俊, 王航, 付光攀, 周灏, 史筱川, 宋晨杰, 陆文斌, 诸葛烽, 黄兴, 吴菲菲
【申请人】江苏省电力公司苏州供电公司, 国家电网公司, 武汉大学, 北京伏佳安达电气技术有限公司
【公开日】2015年9月9日
【申请日】2015年6月8日
【技术领域】
[0001] 本发明涉及一种基于历史故障信息和运行状态的电力电缆运行维护方法，其属于 电力设备资产管理领域，适用于预测电力电缆的故障率，并通过提前制定更换策略来更换 一些运行年份较长的电缆，减缓故障率的增长趋势，提高电缆线路的可靠性，从而提高电网 运行的稳定性。
【背景技术】
[0002] 国内外主要通过建立统计模型分析电缆的故障模式，得出电缆整体的故障概率密 度分布，预测未来的电缆故障数；同时结合电缆的状态以及电缆在系统中的地位，制定电缆 的维修更换计划。通常，电缆的故障分为早期、中期和晚期三个阶段，在进行故障分布的拟 合时，Weibull分布和Crow-AMSAA模型（可靠性增长模型）能够很好的反映这三个阶段，因 此在可靠性领域使用较多，但Weibul1分布和Crow-AMSAA模型存在其局限性。使用Weibul1 分布的不足：
[0003] l)Weibull分布分析的是故障前运行时间，对故障数据的质量要求相对较高，需要 知道设备详细的投运日期和故障日期，当故障前运行时间不满足Weibull分布时，无法使 用Weibul1分布进行分析；
[0004] 2)Weibull分布通常假设电缆为一种不可修复的元件，这显然与电力电缆的实际 运行情况不相符；
[0005] 3)Weibull分布在国外的研宄工作通常都基于大量的老化故障数据，但国内的 电缆故障主要为早期故障，由于早期故障数据的特点与老化故障数据不同，不能简单套用 Weibull分布模型。
[0006] 使用Crow-AMSAA模型的不足：
[0007] 1)Crow-AMSAA模型适用于分析多种混合故障模式的故障数据，对于分析单一故障 模式的故障数据适用性差；
[0008] 2)当分析早期电缆故障数据特点时，由于故障数据样本少，Crow-AMSAA模型抗扰 动性差，分析结果可靠性较差；
[0009] 此外，Weibull分布和Crow-AMSAA模型仅基于电缆历史故障信息来预测电力电缆 故障率也存在明显不足。实际上，预测电缆故障率不仅与电缆历史故障信息有关，还与电缆 实际运行状态有关。电缆运行状态随负荷水平、安装方法以及安装位置的不同而异，从而导 致电缆故障率随运行状态而变化。且进行电缆故障预测的目的就是了解电力电缆运行情 况，并通过提前制定维修更换电缆策略，减缓故障率的增长趋势，提高电缆线路的可靠性。 因此，提出一种能够同时考虑电缆历史故障信息和运行状态的电力电缆运维策略，具有十 分重要的现实意义。
【发明内容】
[0010] 本发明的目的是提供一种基于历史故障信息和运行状态来的电力电缆运维方法， 用来预测电缆剩余寿命，同时提出电缆更换策略，提高整体电缆的可靠性。
[0011] 为达到上述目的，本发明采用的技术方案是：
[0012] 一种基于历史故障信息和运行状态的电力电缆运行维护方法，用于预测电力电缆 的剩余寿命，以及制定电力电缆的更换策略，该方法包括以下步骤：
[0013] (l)weibull比例风险模型的建立：
[0014] 所述的weibull比例风险模型基于Cox比例风险模型和weibull分布而建立，当 所述的Cox比例风险模型的基础风险函数服从所述的weibull分布时，所述的Cox比例风 险模型的函数表达式用作分析所述的电力电缆的故障数据和状态量的weibull比例风险 丰旲型；
[0015] (2)weibull比例风险模型中参数的求解：
[0016] 选取所述的weibull比例风险模型中各个参数的初值，从而求解所述的weibull 比例风险模型中的各个参数；
[0017] (3)电力电缆状态决策曲线的求取：
[0018] 将求解得的所述的weibull比例风险模型中的各个参数代入所述的weibull比例 风险模型，得出所述的weibull比例风险模型的表达式，再根据所述的weibull比例风险模 型的表达式绘制所述的电力电缆的状态决策曲线；
[0019] (4)电力电缆剩余寿命的预测：
[0020] 获取所述的电力电缆的历史状态数据并拟合得到所述的电力电缆的状态信息随 时间变化的关系曲线，获取所述的关系曲线与所述的状态决策曲线的相交点，则所述的电 力电缆当前运行点与所述的相交点之间的时间差即为所述的电力电缆的剩余寿命；
[0021] (5)电力电缆维修更换策略的制定：
[0022] 所述的电力电缆维修更换策略包括最小成本法和最大可用度法：
[0023] ①最小成本法：建立所述的weibull比例风险模型下所述的电力电缆发生故障而 进行更换时的成本率模型，并求解所述的成本率最小时的时间，从而得到所述的电力电缆 更换成本率最小的更换周期；
[0024] ②最大可用度法：建立所述的weibull比例风险模型下所述的电力电缆的可用度 模型，并求解所述的可用度最大时的时间，从而得到所述的电力电缆可用度最大的预防性 维护周期。
[0025] 所述的步骤（1)中，对不同类型的所述的状态量分别进行分析。
[0026] 所述的步骤（2)中，利用迭代法求解所述的weibull比例风险模型。
[0027] 所述的步骤（3)中，根据所述的weibull比例风险模型的表达式获得所述的电力 电缆的可靠度，使所述的电力电缆的可靠度大于满足要求的最小可靠度并分别取自然对 数，以时间为横坐标、以取自然对数后的电力电缆的可靠度为纵坐标，从而绘制所述的电力 电缆的状态决策曲线。
[0028] 所述的步骤（3)中，所述的电力电缆的可靠度取两个不同值分别绘制所述的电力 电缆的状态决策曲线。
[0029] 所述的电力电缆运行维护方法通过电力电缆状态监测和智能评估系统自动实现。
[0030] 由于上述技术方案运用，本发明与现有技术相比具有下列优点：
[0031] l)Weibull比例风险模型同时分析电缆历史故障数据和电缆运行状态（如护层环 流、最大局部放电量变化率与介损变化率），进行电缆故障预测，更加接近电缆的实际情况， 提高了预测的合理性、可靠性；
[0032] 2)根据Weibull比例风险模型可绘制状态决策曲线，直观明了的显示出电缆的运 行情况，可为电缆是否需要维修更换提供参考；
[0033] 3)如果已知某电缆的历史状态信息，Weibull比例风险模型通过拟合可得到电缆 的状态信息随时间变化的关系曲线，预测电缆的剩余寿命；
[0034] 4)Weibull比例风险模型的基础上可提出最小成本法和最大可用度法，结合考虑 经济因素，制定电缆维修更换策略。
【附图说明】
[0035] 附图1为本发明的电力电缆运行维护方法的流程示意图。
[0036] 附图2为电力电缆的Weibull比例风险模型h(t，X)示意图。
[0037] 附图3为电力电缆的状态决策曲线示意图。
[0038] 附图4为电力电缆的状态信息随时间变化的关系曲线与状态决策曲线的相交点 示意图。
[0039] 附图5为电力电缆更换成本率最小决策图。
[0040] 附图6为电力电缆预防性维修最佳时间决策图。
【具体实施方式】
[0041] 下面结合实施例对本发明作进一步描述。
[0042] 实施例一：一种用于预测电力电缆的剩余寿命，以及制定电力电缆的更换策略的 基于历史故障信息和运行状态的电力电缆运行维护方法，如附图1所示，其在电力电缆状 态监测和智能评估系统中通过weibull比例风险模型的建立、weibull比例风险模型中参 数的求解、电力电缆状态决策曲线的求取、电力电缆剩余寿命的预测、电力电缆维修更换策 略的制定这几个步骤实现。
[0043] (1)weibull比例风险模型的建立
[0044]weibull比例风险模型基于Cox比例风险模型和weibull分布而建立，当Cox比例 风险模型的基础风险函数服从weibull分布时，Cox比例风险模型的函数表达式用作分析 电力电缆的故障数据和状态量的 weibull比例风险模型。
[0045]Cox比例风险模型的函数表达式如式（1)
[0047] 其中，X#依时协变量，对应为电力电缆的运行状态信息，Y对应的回归系 数，n为依时协变量的个数，0为Weibull分布形状参数，t为时间，n为Weibull尺度参 数。此时模型称为Weibull比例风险模型，其能够同时分析故障数据和在线监测数据。
[0048] 采用上述Weibull比例风险模型分析电力电缆时，对不同类型的状态量分别进行 分析。在应用Weibull比例风险模型前，首先对纳入模型的状态量进行选择。比如，局部放 电和介损反映的是电缆绝缘内部缺陷或老化的状态，而护层环流反映的是绝缘外部的缺陷 或故障，因此在对电缆的状态进行评估时，将局部放电和介损归为一类，将护层环流归为另 一类，分别进行分析。由于局放特征量与电缆绝缘材料、绝缘缺陷形式、大小、电缆所受应力 等均有关系，局放特征量与绝缘老化程度的关联度尚不清晰，同时介损的测量也受到环境 因素以及测量误差的影响，因此将这两个特征量纳入模型进行分析时，以单位时间内的介 损的变化率和单位时间内的局部放电能量作为变量。
[0049] (2)weibull比例风险模型中参数的求解
[0050] 求解方法为：选取weibull比例风险模型中各个参数的初值，从而利用迭代法求 解weibull比例风险模型中的各个参数。weibull比例风险模型中参数的求解包括求解含 一个协变量的weibull比例风险模型和求解两个协变量的weibull比例风险模型。
[0051] ①求解含一个协变量的weibull比例风险模型
[0052] 含有一个协变量的weibull比例风险模型的函数表达式如式（2)
[0054] 其中，X为依时协变量，a为X对应的回归系数。
[0055] 可靠度表示电力电缆在时间t内正常工作的概率，Weibull比例风险模型对应的 可靠度函数为：
[0057] 故障概率密度函数为
[0058]f(t,X) =h(t,X) ?R(t,X) (4)
[0059] 构造故障概率密度函数的似然函数：
[0061] 其中n为数据的总数，m为故障数据个数，n-m表示截尾数据个数。
[0062] 两边取自然对数得
[0064] 为了求解待求参数为0，n和a，对上式分别求偏导数和二阶偏导数
[0066]
[0076]
[0077] 则对数似然函数的二阶求导矩阵
[0079]然后通过Newton-Raphson迭代，直至收敛，贝丨J :
[0081] 当k = 0时，代入0，n和a的初值0。，qQ和a。进行求解。
[0082] ②求解含两个协变量的weibull比例风险模型
[0083] 当weibull比例风险模型中有两个协变量时，求解模型参数时，需对上述求解方 法重新构造迭代式子。其中两个协变量时，Weibull比例风险模型函数为：
[0085] 故障概率密度函数为
[0086] f (t, X) = h(t, X) ? R(t, X) (22)
[0087] 构造故障概率密度函数的似然函数
[0089] 其中n为数据的总数，m为故障数据个数，n-m表示截尾数据个数。
[0090] 两边取自然对数，得：
[0092] 为了求解待求参数为0，n，a种a2，对上式分别求0，q，a种a2的偏导数 和二阶偏导数。
[0093]
[0103]
[0113] 则对数似然函数的二阶求导矩阵：
[0114]
[0115] 然后通过Newton-Raphson迭代，贝丨J:
[0117] 当k= 0时，代入0，q，a兩a2的初值0 〇，q〇，a1〇和a2〇进行求解。
[0118] 现阶段由于受条件限制，目前很难获得电缆的全部状态数据，因此对状态数据进 行了一些假设，如表1所示，其中电缆状态为1表示发生过故障，电缆状态为0表示未发生 故障。对于故障电缆，假设一些整体处于晚期故障模式的数据。
[0119] 表1电缆状态数据
[0121]
[0122] (3)电力电缆状态决策曲线的求取
[0123] 将迭代求解得的weibull比例风险模型中的各个参数代入weibull比例风险模 型，得出weibull比例风险模型的表达式；再根据weibull比例风险模型的表达式绘制电力 电缆的状态决策曲线。绘制电力电缆的状态决策曲线时，根据weibull比例风险模型的表 达式获得电力电缆的可靠度，使电力电缆的可靠度大于满足要求的最小可靠度并分别取自 然对数，以时间为横坐标、以取自然对数后的电力电缆的可靠度为纵坐标，从而绘制电力电 缆的状态决策曲线。
[0124] 以上述含两个协变量的weibull比例风险模型为例，选取初值|3 Q= 1. 8，n。= 10000,a1(|= 〇.〇〇1，a2。= 〇.〇9,代入程序进行迭代，得到 0 = 3.246,n= 114252,ai =0. 078,a2= 〇. 022,Weibull比例风险模型的函数形式如下：
[0126] 则电力电缆的Weibull比例风险模型h(t，X)如附图2所示。
[0127] 假如电缆系统所要求的最小可靠度为心，为了保证电缆的安全运行，电缆的可靠 度必须大于&，则有：
[0129] 两边分别取自然对数有
[0130] aj ?Xi<ln(-lnR〇) - 0 (lnt-lnn) (49)
[0131] 以t为横坐标，ln(-lnRQ)-f3 (lnt-lnq)为纵坐标绘制曲线即为决策曲线。
[0132] 可以对电力电缆的可靠度取两个不同值分别绘制电力电缆的状态决策曲线，如分 别绘制系统最小可靠度为0. 95和0. 99时的决策曲线，如附图3所示，代入某时刻t电力电 缆的状态量X，计算a*X，如果点（t，a 落在虚线上方，表明此时需要对电力电缆进行 维修更换；如果点（t，a 落在实线下方，表明电力电缆仍可继续工作；如果点（t，a*X) 落在实线与虚线之间，这部分区域为临界区，表明需要对电力电缆状态进行密切关注，重点 监测。
[0133] (4)电力电缆剩余寿命的预测
[0134] 获取电力电缆的历史状态数据并通过数学拟合得到电力电缆的状态信息随时间 变化的关系曲线，获取关系曲线与状态决策曲线的相交点，则电力电缆当前运行点与相交 点之间的时间差即为电力电缆的剩余寿命，如附图4所示。假设电力电缆的状态信息随时 间变化的关系曲线（图中圆圈表示）与状态决策曲线相交于点B，计算点A(当前运行点） 和B的时间差h即为电缆的剩余寿命。如果电缆的状态量从点A开始变化缓慢，则该关系 曲线（图中三角表示）与决策曲线相交于点C，此时的剩余寿命&大于h。为了提高电缆 寿命预测的精度，需要对电缆的状态量进行连续监测，收集足够多的状态信息，从而实现对 电缆寿命的预测。
[0135] (5)电力电缆维修更换策略的制定：
[0136] 电力电缆维修更换策略包括最小成本法和最大可用度法。
[0137] ①最小成本法：建立weibull比例风险模型下电力电缆发生故障而进行更换时的 成本率模型，并求解成本率最小时的时间，从而得到电力电缆更换成本率最小的更换周期。
[0138] 电缆发生故障后，目前一般采取更换措施。假设购买电缆的费用为，电缆故障后 的损失为C2，电缆的更换周期为t，即电缆运行时间达到t时，即使电缆没有发生故障也对 电缆进行更换。则电缆更换的成本率为：
[0140] 对成本率求导可得：
[0142] 其中
h(t)为故障率函数，f(t)为故障概率密 度函数，F(t)为累积故障概率密度函数，R(u)为残存函数，X表示状态量，则当
[0144] 时成本率最小，求出此时的t即为电缆更换成本率最小的更换时间。
[0145] 对上式进一步变换可得
[0147]根据已求出的h(t，X)的函数表达式，可以求出F(t，X)和R(t，X) 。由于
积分困难，首先在x轴（0~5000)范围内绘制R(u，X)的图形，然后通过Origin 软件对该图形进行数值积分。在进行数值积分时，Origin对x轴（0~5000)范围内均勾 划分为99个区间Xl= [50. 5,101，…，5000]，求得数值积分的值yi= [50. 5,100. 9，…， 2333. 2]。然后以Xl为横坐标，y 纵坐标，绘制并拟合曲线，该拟合的曲线方程即为
的积分表达式，以上式中的参数为例时
[0149] 假设(VC2= 1，当X:= 150,X:= 30时，分别绘制h(t,X)和
这两条 曲线，如图5所示，两者相交于点A，点A对应的时间为2018，表明当电缆的状态量分别为Xi=150,X2= 30时，最佳的更换时间为电缆运行的第2018天。
[0150] 如表2和3所示，讨论了当协变量值不变，（VC2比值发生变化时对最佳更换时间 的影响，发现随着(Vc2比值的增大，最佳更换时间也在增大，即如果故障损失与购买电缆 成本相比很小，那么电缆最佳更换时间会增大，电缆可用时间增大；当Q/^比值不变，发现 随着协变量&值的增大，最佳更换时间在减小，这表明当表征电缆的状态变化较快时，应尽 快更换电缆。
[0151] 表2 (VC2比值变化时对最佳更换时间的影响（X150,X2= 30)
[0153] 表3 值变化时对最佳更换时间的影响（X2= 30,C^4= 1)
[0154]
[0155] ②最大可用度法：建立weibull比例风险模型下电力电缆的可用度模型，并求解 可用度最大时的时间，从而得到电力电缆可用度最大的预防性维护周期。
[0156] 考察可用度时，以一个预防性维修周期为例进行说明。假设在每一个预防性维修 周期内，电缆的可用时间为Ta，电缆的不可用时间为Tb，电缆的预防性维修的时间为，故 障后维修时间为T2。因此在每一个预防性维修周期内电缆不可用时间为
[0157]Tb=R(t,X) ?Ti+d-Ra.X)) ?T2 (55)
[0158] 根据可用度的定义
[0160]则Weibull比例风险模型下的可用度为
[0162] 令
为使可用度最大，贝U只需B最小，由于积分
不易求得，另1(1天），T2= 2，X1= 150，X2= 30,绘制曲线B，如附图6所 示，得到A点为曲线的最低点，进而得出A点对应的x轴值为1832,表明使可用度最大化的 预防性维修的时间间隔为1832天。
[0163] 如表4和5所示，讨论了当协变量值和T2值不变，Ti值发生变化时对最佳预防性 维修时间的影响，发现随着1\值的增大，最佳预防性维修时间也在增大；当T:和T2值以及 &值不变时，发现随着协变量值的增大，最佳预防性维修时间在减小，这表明当表征电缆的 状态变化较快时，应尽快对电缆进行维修。
[0164] 表4 值变化时对最佳预防性维修时间的影响（X150,X2= 30,T2= 2)
[0165]
[0166]表5Xi值变化时对最佳预防性维修时间的影响（X2= 30,Tf1，T2= 2)
[0168] 本发明的基于历史故障信息和运行状态的电力电缆运维策略，基于Cox比例风险 模型和Weibull分布，建立Weibull比例风险模型。Weibull比例风险模型能够同时分析故 障数据并在线监测状态数据，从而能够根据Weibull比例风险模型对电缆历史故障信息和 运行状态的分析，预测电缆剩余寿命，提出最小成本法和最大可用度法，用来指导制定电力 电缆的更换策略。
[0169] 上述实施例只为说明本发明的技术构思及特点，其目的在于让熟悉此项技术的人 士能够了解本发明的内容并据以实施，并不能以此限制本发明的保护范围。凡根据本发明 精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种基于历史故障信息和运行状态的电力电缆运行维护方法，用于预测电力电缆的 剩余寿命，以及制定电力电缆的更换策略，其特征在于：该方法包括以下步骤： (Dweibull比例风险模型的建立： 所述的weibull比例风险模型基于Cox比例风险模型和weibull分布而建立，当所述 的Cox比例风险模型的基础风险函数服从所述的weibull分布时，所述的Cox比例风险模 型的函数表达式用作分析所述的电力电缆的故障数据和状态量的weibull比例风险模型； (2) weibull比例风险模型中参数的求解： 选取所述的weibull比例风险模型中各个参数的初值，从而求解所述的weibull比例 风险模型中的各个参数； (3) 电力电缆状态决策曲线的求取： 将求解得的所述的weibull比例风险模型中的各个参数代入所述的weibull比例风险 模型，得出所述的weibull比例风险模型的表达式，再根据所述的weibull比例风险模型的 表达式绘制所述的电力电缆的状态决策曲线； (4) 电力电缆剩余寿命的预测： 获取所述的电力电缆的历史状态数据并拟合得到所述的电力电缆的状态信息随时间 变化的关系曲线，获取所述的关系曲线与所述的状态决策曲线的相交点，则所述的电力电 缆当前运行点与所述的相交点之间的时间差即为所述的电力电缆的剩余寿命； (5) 电力电缆维修更换策略的制定： 所述的电力电缆维修更换策略包括最小成本法和最大可用度法： ① 最小成本法：建立所述的weibull比例风险模型下所述的电力电缆发生故障而进行 更换时的成本率模型，并求解所述的成本率最小时的时间，从而得到所述的电力电缆更换 成本率最小的更换周期； ② 最大可用度法：建立所述的weibull比例风险模型下所述的电力电缆的可用度模 型，并求解所述的可用度最大时的时间，从而得到所述的电力电缆可用度最大的预防性维 护周期。2. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（1)中，对不同类型的所述的状态量分别进行分析。3. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（2)中，利用迭代法求解所述的weibull比例风险模型。4. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（3)中，根据所述的weibull比例风险模型的表达式获得所述的电 力电缆的可靠度，使所述的电力电缆的可靠度大于满足要求的最小可靠度并分别取自然对 数，以时间为横坐标、以取自然对数后的电力电缆的可靠度为纵坐标，从而绘制所述的电力 电缆的状态决策曲线。5. 根据权利要求4所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：所述的步骤（3)中，所述的电力电缆的可靠度取两个不同值分别绘制所述的电 力电缆的状态决策曲线。6. 根据权利要求1所述的基于历史故障信息和运行状态的电力电缆运行维护方法，其 特征在于：其通过电力电缆状态监测和智能评估系统自动实现。
【专利摘要】本发明涉及一种用于预测电力电缆的剩余寿命，以及制定电力电缆的更换策略的基于历史故障信息和运行状态的电力电缆运行维护方法，其通过weibull比例风险模型的建立、weibull比例风险模型中参数的求解、电力电缆状态决策曲线的求取、电力电缆剩余寿命的预测、电力电缆维修更换策略的制定这几个步骤实现。该方法提高了电缆故障预测的合理性、可靠性，直观明了的显示出电缆的运行情况，可为电缆是否需要维修更换提供参考，在Weibull比例风险模型的基础上可提出最小成本法和最大可用度法，结合考虑了经济因素，能够制定电缆维修更换的合适策略。
【IPC分类】G06Q50/06, G06Q10/06
【公开号】CN104899690
【申请号】CN201510308384
【发明人】张梁, 张俊, 郑勇, 黄晓军, 周承科, 周文俊, 王航, 付光攀, 周灏, 史筱川, 宋晨杰, 陆文斌, 诸葛烽, 黄兴, 吴菲菲
【申请人】江苏省电力公司苏州供电公司, 国家电网公司, 武汉大学, 北京伏佳安达电气技术有限公司
【公开日】2015年9月9日
【申请日】2015年6月8日
最新回复(0)