1.一种智能变电站保护系统的检修周期优化方法，其特征在于，包括：S1、根据各保护设备的尺度参数、形状参数计算得到满足威布尔分布的故障率；

S2、计算在不同检修类型下各保护设备的改善因子；

S3、根据所述各保护设备的满足威布尔分布的故障率和改善因子，计算得到计及计划检修的各设备的故障率；

S4、根据可靠性框图法，计算得到未考虑计划检修的各采样跳闸模式下单套配置或双重化配置的保护系统的故障率，得到计及计划检修的各采样跳闸模式下单套配置或双重化配置的保护系统的故障率；

S5、构建各保护系统的稳定运行期与损耗期的时间分界点的计算模型；

S6、计算智能变电站保护系统的分阶段变周期的故障率；

S7、构建智能变电站保护系统的可靠运行时间的计算模型；

S8、构建智能变电站保护系统的运行成本的计算模型；

S9、根据所述运行成本的计算模型，采用试值法，计算得到保护系统在可靠运行时间内年均运行成本最小情况下保护系统的稳定运行期的最优计划检修周期、损耗期的最优计划检修周期。

2.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S1中计算得到满足威布尔分布的故障率的方法为：根据各保护设备的尺度参数α、形状参数β，计算各保护设备满足威布尔分布的故障率λ(t)：其中，t为设备的实际年龄。

3.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S2中计算在不同检修类型下各保护设备的改善因子的方法为：基于不完全检修的原则，设备在检修之后并非修复如新，采用有效工龄、改善因子来表示不完全检修对设备运行状态的改善作用，设τ为设备运行的有效工龄，则设备在第k次维修后的有效工龄τk为：其中，T为计划检修周期，t为设备的实际年龄，qi为第i次维修时的改善因子；

计算第k次维修时的改善因子qk：

其中，li(i＝1,2,3,4)为与设备相关的常数，由专家根据设备的历史统计数据及其运行特性进行赋值，q1为第一次检修的改善因子，改善因子随着检修次数的增加逐渐减小；

为了保障智能变电站的安全稳定运行，对保护系统进行定期检修，根据检修程度分为A、B、C、D四类检修；A类检修为大修，对所有设备进行全面彻底的检查和修理；B类检修为故障检修，针对设备存在的问题，对其部分部件进行解体检查和修理；C类检修根据设备的老化规律进行的维护；D类检修是在设备总体运行良好的情况下进行的定期消缺检修；

由于大修时的检修成本投入较大，其改善因子比一般性维修的改善因子大Δq：Δq≤1-qi    (4)

根据式(3)-(4)得到A类大修时的改善因子qk：

B类故障检修只是将设备从故障状态变成正常运行状态，暂不考虑其对设备的改善作用，故设备的各类故障检修的改善因子为：其中，qk为第k次维修时各类故障检修的改善因子。

4.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S3中计算得到计及计划检修的设备的故障率的方法为：其中，λ(t)为计及计划检修的设备的故障率。

5.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S4的具体方法包括：利用可靠性框图法，得到直采直跳、直采网跳、网采网跳的采样跳闸模式下单套配置或双重化配置的各保护系统的可靠性框图，并根据可靠性框图中各设备在功能上的逻辑关系，得到未考虑计划检修的各采样跳闸模式下单套配置或双重化配置的保护系统故障率λ′s(t)：其中，m为保护设备的种类数；nj为第j类保护设备的等效台数；αj为第j类保护设备的尺度参数；βj为第j类保护设备的形状参数；

根据可靠性框图中各设备在功能上的逻辑关系，得到计及计划检修的各采样跳闸模式下单套配置或双重化配置的保护系统的故障率λs(t)：由于不同电压等级下各保护系统采用不同的配置方式，即保护单套配置或保护双重化配置；单套保护相关设备在功能逻辑上属于串联系统，因此单套保护的可靠性等于各类设备等效可靠运行概率P(t)的乘积；双重化保护的各类双重化配置设备故障率为单套保护对应部分设备故障率的一半，为此采用各类设备的等效台数nj来表示，第j类双重化配置设备的等效台数nj为该类设备实际台数的1/4。

6.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S5中构建各保护系统的稳定运行期与损耗期的时间分界点计算模型的方法为：设置故障率的变化量Δλ，考虑计划检修时各保护系统稳定运行期与损耗期的时间分界点的t1是由步骤S4中获得的两个相邻时刻的故障率满足式(10)的时间最小值：t1＝min(t|Δλ≤λs(t)-λs(t-Δt))    (10)其中，t1为考虑计划检修的情况下保护系统的稳定运行期与损耗期的分界时刻；Δλ为经验值，其取值与Δt的大小呈正相关；Δt为仿真时间步长。

7.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S6中计算智能变电站保护系统的分阶段变周期的故障率的方法为：当考虑计划检修情况下的运行时间等于t1时，稳定运行期已按照检修周期T1做了k1次检修，则稳定运行期的计划检修累积效果Q1为：其中，floor[]为向下取整函数；

当运行时间为t1时保护系统进入损耗期，计划检修周期变为T2，此时保护系统的故障率为λs(t1)；将损耗期的初始时刻t1看作一个新阶段的开始，则与t1时刻的故障率λs(t1)具有相同的故障率、且未考虑计划检修时系统故障率曲线上的等效时刻为t1′，根据步骤S4和式(12)得到t1'：λ′s(t1')＝λs(t1)    (12)

保护系统在损耗期故障率的计算方式与其在稳定运行期的类似，则保护系统在损耗期第一个检修周期内的故障率为λs(t)：损耗期检修作用的改善因子应从第k1+1次开始，则损耗期的第k次计划检修后的累积效果为Q2：考虑计划检修的损耗期的故障率λs(t)：

将式(11)和(14)代入(15)得到λs(t)：

智能变电站保护系统的分阶段变周期的故障率λs(t)：

由于检修周期的变化及计划检修和设备老化的累积效果，损耗期的故障率计算方式与稳定运行期有所不同，设T1为保护系统的稳定运行期的计划检修周期，T2为保护系统的损耗期的计划检修周期。

8.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S7中构建智能变电站保护系统的可靠运行时间计算模型的方法为：根据可靠运行概率与故障率的关系式，得到分段检修的保护系统可靠运行概率Ps(t)：假设R为保护系统的最低可靠运行概率，计算得到线路保护系统的可靠运行时间ta：ta＝max(t|Ps(t)≥R)    (19)

其中，在工程实际中，当系统的可靠性低于一定数值时，则认为该系统失效，此时对应的最大时刻为系统的可靠运行时间。

9.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S8中构建智能变电站保护系统的运行成本的计算模型的方法为：收集来自设备厂家的保护设备的故障维修成本、设备维护成本及设备购置成本；假设各设备在稳定运行期的维护成本为CD-pm，损耗期的检修成本为CC-pm，大修的检修成本为CA-pm；保护系统在可靠运行时间ta内的总成本CS-total由计划检修的维护总成本Cpm-total、故障检修总成本Cmr-total、设备购置总成本Cpr-total三部分组成；

各保护系统在可靠运行时间ta内维护总成本Cpm-total等于该保护系统中每类设备在不同阶段不同类型的维护总成本Cj-pm之和；

其中，第j类设备的维护总成本Cj-pm等于第j类设备的实际台数hj、检修次数N、第j类设备的单台维护成本 的三者乘积；

各保护系统在可靠运行时间内故障维修总成本Cmr-total等于该保护系统中每类设备的故障维修总成本Cj-mr之和；

其中，第j类设备的故障维修总成本Cj-mr等于第j类设备的实际台数hj、第j类设备的单台故障维修成本 第j类设备的单台故障概率 的三者乘积；

各保护系统在可靠运行时间内设备购置总成本Cpr-total等于该保护系统中每类设备的设备购置总成本Cj-pr之和，第j类设备的购置总成本Cj-pr等于第j类设备的实际台数hj、第j类设备的单台购置成本 的两者乘积：其中，m为系统内设备的种类数；λj(t)为第j类设备的等效故障率函数；hj为第j类设备的实际台数； 为第j类设备的故障维修成本； 为第j类设备的购置成本， 为稳定运行期的维护成本， 为第j类设备在损耗期的检修成本， 为第j类设备在修的检修成本；TM为大修周期；

可得到总成本CS-total：

CS-total＝Cpm-total+Cmr-total+Cpr-total    (23)保护系统在可靠运行时间内的年均运行成本为：

其中，c(T1,T2)为保护系统在可靠运行时间内的年均运行成本。

10.根据权利要求1所述的智能变电站保护系统的检修周期优化方法，其特征在于，所述步骤S9的具体方法包括：采用试值法，分别针对单套和双套的线路保护系统、母线保护系统、主变保护系统，先列出各保护系统在直采直跳、直采网跳、网采网跳的采样跳闸模式下可能存在的T1、T2组合，再将T1、T2的各组合分别代入运行成本的计算模型，计算得到各组相应的c(T1,T2)；最后将其中使c(T1,T2)最小的T1、T2组合作为计划检修周期的最优解，即求解得到c(T1,T2)最小情况下保护系统的稳定运行期的最优计划检修周期T1、损耗期的最优计划检修周期T2。