本发明属于配电网可靠性评估方法，尤其涉及一种基于马尔科夫链的配电网序贯可靠性评估方法。
背景技术
0、
背景技术：
1、据电力公司统计，大部分停电事故是由配电网故障引起的。配电网的分析与评价具有重要意义。然而，配电网的结构往往复杂多样，难以进行此类分析和评估。例如，现有的规划中经常采用环网或多分段多联络结构，利用多个联络开关连接其他部分。另外，配电网中元件种类繁多，因此配电网可靠性评估是一个非常复杂的问题。
2、现有的可靠性评估方法基本上考虑了n-1原则，没有考虑部件的概率因素。如果n元件配电系统的任何一元件(馈线、变压器等)因故障而被切断，则不应造成潮流破坏或负荷削减。n-1准则忽略了系统失效状态的可能性，因此仅考虑n-1准则的可靠性评估是不准确的。同时n-1原理忽略了多重故障事件发生的可能性，从而导致对配电系统可靠性的高估。
3、配电网常用的考虑系统状态可能性的可靠性评估方法可分为两类：解析法和模拟法。解析法包含网络法和系统状态枚举法。网络法是基于配电网的拓扑结构对配电进行分析，具体包括故障模式影响分析法、网络等效法、最短路径法、最小割集法和故障扩散法等方法。状态枚举方法直接枚举所有的系统状态，同时忽略状态之间的转换关系。所有的解析法都以元件可靠性模型为基础，通过列举系统故障进行分析。但系统的失效状态数直接随系统构件数和需要考虑的故障阶数的增加而增加，当系统较大且需要考虑的故障阶数较高时，系统会出现大量的失效状态，分析方法的计算量也会很大。此外，解析法是一种非序贯方法，频率指标仅通过元件故障率进行计算，在考虑了高阶故障的情况下，其频率相关可靠性指标会与真实值出现较大误差。
4、典型的模拟方法是蒙特卡洛模拟法，同样可分为两类：非序贯蒙特卡洛法(nmcs)和序贯蒙特卡洛法(smcs)。其中nmcs只对系统状态进行采样，因此可以提供系统概率和期望指数。smcs不仅对系统状态进行采样，还对其持续时间进行采样。因此，可以考虑时间因素，并获得综合指数，包括指数的概率分布。显然，smc比nmc需要更多的运行时间。mcs的一个重要优点是它对系统维度与结果精度之间的关联较小。因此，相较于解析法，蒙特卡洛法风适合应用于大型系统。为了提高其收敛速度，提出了各种改进方法，如控制变量、简化状态空间等。然而，对于可靠性较高的系统，异常状态的采样依旧十分困难。
5、电力系统中可修复部件的运行通常被视为连续时间马尔可夫链(ctmc)。ctmc模型是一种通过建立状态空间图求解频率，持续时间和稳态概率等可靠性指标。该方法已用于分析单个部件和简单电力系统，因此有可能将该方法与其他解析法相结合发展为序贯分析方法，提高频率相关指标的计算精度。但ctmc模型状态空间的维数是一个指数函数，同样受系统维度较大影响，无法直接在大规模系统中使用。虽然已有一些基于状态融合的方法对状态空间图进行了简化，但是在大规模系统中的应用依然不甚理想，仍需进一步改进。
6、参考文献
7、[1]r.billinton and p.wang.“reliability-network-equivalent approach todistribution-system-reliability evaluation,”.ieee proceedings-generation,transmission and distribution,vol 145,no.2,pp.149-153,mar.1998.
8、[2]g.t.heydtandt.j.graf.“distribution system reliability evaluationusing enhanced samples in a monte carlo approach”.ieee transactions on powersystems,vol 25,no.4,pp.2006-2008,nov.2010.
9、[3]r.billinton andw.li.reliability assessment of electric powersystems using monte carlo methods.new york:springer us,1994
10、[4]y.tao,t wei,w.yue,et al.“reliability assessment of a multi-statedistribution system with microgrids based on an accelerated monte-carlomethod,”.iet generation,transmission&distribution,vol.12,no.13,pp.3221-3229,jul.2018.
11、[5]g.li,z.bie,h.xie,et al.“customer satisfaction based reliabilityevaluation of active distribution networks,”applied energy,vol.162,no.1,pp.1571-1578,jan.2015.
12、[6]k.hou,z.he,y.wang,et al.“reliability modeling for integratedcommunity energy system considering dynamic process of thermal loads”.ietenergy systems integration,vol.1,no.3,apr.2019.
13、[7]k.hou,h.jia,x.xu,et al.“a continuous time markov chain basedsequential analytical approach for composite power system reliabilityassessment,”in power&energy society general meeting.boston,usa,jul.2016,pp.1-5.
14、[8]j.huang,z.li,and h.c.liu.“new approach for failure mode and effectanalysis using linguistic distribution assessments and todim method,”reliability engineering&systemsafety,vol.167,no.1,pp.302-309,nov.2017.
15、[9]c.juanwei,y.tao,x.yue,et al.“fast analytical method forreliability evaluation of electricity-gas integrated energy systemconsidering dispatch strategies,”applied energy,vol.242,no.15,pp.260-272.mar.2019.
16、[10]y.liu and c.singh.“reliability evaluation of composite powersystems using markov cut-set method,”ieee transactions on power systems,vol.25,no.2,pp.777-785,may.2010.
技术实现思路
1、针对现有技术存在技术问题，本发明提出了一种基于马尔科夫链的配电网序贯可靠性评估方法。在引入马尔科夫链模型后，由于需要对节点进行分析，序贯可靠性分析的计算量急剧增加，因此需要新的简化计算方法。对于负荷节点，由于配电网的结构较大，存在较多分枝，因此会存在一些故障与否都不影响负荷节点状态的元件，也会存在一些元件故障将直接导致负荷削减的元件。元件分组就是将涉及这些元件的状态挑选出来并进行相应的简化，减少状态空间图中的状态数量，减少计算量。在完成简化后，使用新的频率计算公式计算相应可靠性指标。
2、该发明核心内容可归结如下：
3、一种基于马尔科夫链的配电网序贯可靠性评估方法，包括如下步骤：
4、步骤1：获取配电网系统中节点结构、故障率和故障阶数阀值；
5、步骤2：使用影响增量法对配电网系统的所有节点状态的负荷削减情况进行计算，并将负荷节点在配电网系统状态下的负荷削减情况和在配电网系统状态下各元件之间相关性按照矩阵方式储存；
6、步骤3：对第一个负荷节点遍历在所有配电网系统状态下负荷削减情况划分；
7、步骤4：将步骤3得到的节点状态按照不相关元件、单一相关元件的状态进行分组获得第一负荷状态削减配电网系统；
8、步骤5：将步骤3中未分组相关元件的正常状态进行融合获得第二负荷状态削减配电网系统；
9、步骤6：根据第二负荷状态削减配电网系统构建融合状态空间图和配电网系统转移率矩阵；
10、步骤7：按照如下公式对所述第三负荷削减配电网系统计算节点的年平均停电次数；
11、
12、其中：i为负荷节点编号，m为分组编号，mmax为分组上限；
13、步骤8：重复步骤3至7，直到将配电网系统中所涉及到的所有负荷节点遍历完成；
14、步骤9：按照如下公式获得配电网系统的可靠性指标saif；
15、
16、其中：fi是负荷节点i的停电频率，ni是负荷节点i的用户数，imax为系统中的负荷节点数。
17、进一步地，所述步骤2使用影响增量法对配电网系统的所有节点状态的负荷削减情况进行计算过程，包括：
18、令某一元件s1为配电网系统负荷节点lpm的不相关元件，其影响为0，即：is1，lpm＝0；
19、当s1与配电网系统内另一元件同时故障时的系统削减情况为：
20、
21、其中，plpm为负荷节点lpm的负荷数值，根据影响增量法高阶故障独立性可知其影响增量为
22、
23、其中，δi{s1.sn},lpm为负荷节点lpm在元件s1和sn同时故障时的影响增量；
24、当元件sn满足以下判据时为lpm的不相关元件：
25、
26、其中，sr为除去系统中除去sn的其余元件，该判据的含义为如果元件sn发生故障时不会导致该负荷节点lpm出现负荷削减。
27、进一步地，所述步骤4按照不相关元件、单一相关元件的状态进行分组获得第一负荷状态削减配电网系统过程，包括：
28、建立不相关元件分组系统状态的数学模型
29、πu＝1
30、qu＝0
31、其中，u是所有负荷节点y包含元件x故障的状态集合，u的定义如下：
32、u(x,y)＝{s|x∈s}
33、建立独立元件分组系统状态的数学模型
34、πu＝ax
35、qu＝λu
36、其中，x为系统元件，s为由故障元件组成的系统状态；
37、根据独立元件分组系统状态的数学模型构建负荷节点的单一相关元件的判定条件为：
38、
39、式中，u是所有包含元件x故障的状态集合，该集合内的所有状态均会导致节点y产生负荷削减，u的定义如下
40、
41、通过如下公式将配电网系统中已经分组的不相关元件、单一相关元件的概率除去获得第一负荷状态削减配电网系统:
42、
43、其中，px为分组元件x的概率。
44、进一步地，所述步骤5按照未分组相关元件的正常状态进行融合获得第二负荷状态削减配电网系统过程，包括：
45、建立未分组相关元件的数学模型：
46、
47、按照如下公式进行融合获得第二负荷状态削减配电网系统：
48、
49、
50、
51、
52、
53、其中：n为被融合的状态总数。
54、有益效果
55、本方法在改进后的rbtsbus6算例中进行了测试，其中包括30个馈线段、26个节点、23个配电变压器、23个负载点、断路器和断开装置。如图6所示，节点40通过联络开关连接到另一个馈线，其外部特性可以用6kw容量的电源来模拟。线路和变压器的平均故障率设为0.014occ/年，线路mttr为8h，变压器平均修复时间(mttr)为15h。以序贯蒙特卡洛法5000年的结果为基准，其中eens为11.76mwh/年，saifi为0.1947次/年。
56、从表1中可以看出，随着枚举阶数的增加，传统枚举法和马尔科夫链都能有效提升所得结果的精度，但是传统枚举法的提升效果不如马尔科夫链法的效果明显。在仅考虑n-1情况时，传统枚举法的saifi精度要略高于ctmc，这是因为根据计算公式，在n-1情况下，ctmc的每一项为传统解析法的每一项乘以当前非故障状态的概率，因此n-1情况下的ctmc的saifi误差大于传统枚举法。但是在n-2情况下，ctmc对于saifi的计算精度要高于传统解析法，且在影响增量法能够有效减少高阶故障计算时间这一优势的帮助下，ctmc的运行时间也明显小于传统解析法。在n-2情况下，ctmc和1000年smcs相比，ctmc仅有0.30％的误差，小于smcs0.42％的误差，此外，在运行时间上，ctmc也仅为smcs的2/3左右。
57、表1低故障率时系统可靠性指标对比
58、
59、
60、图7是smc的运行时间和结果误差曲线，图上的点是考虑n-2故障情况下的ctmc-iise和传统解析法的运行时间及误差。从图中能够看出ctmc-iise在计算效率上的优势，ctmc-iise在100s以下就能达到小于1％的误差，而smc需要接近1000s才能稳定在同样的水准。
1.一种基于马尔科夫链的配电网序贯可靠性评估方法，其特征在于：
2.根据权利要求1所述的一种基于马尔科夫链的配电网序贯可靠性评估方法，其特征在于：所述步骤2使用影响增量法对配电网系统的所有节点状态的负荷削减情况进行计算过程，包括：
3.根据权利要求1所述的一种基于马尔科夫链的配电网序贯可靠性评估方法，其特征在于：所述步骤4按照不相关元件、单一相关元件的状态进行分组获得第一负荷状态削减配电网系统过程，包括：
4.根据权利要求1所述的一种基于马尔科夫链的配电网序贯可靠性评估方法，其特征在于：所述步骤5按照未分组相关元件的正常状态进行融合获得第二负荷状态削减配电网系统过程，包括：
