高渗透率下分布式决策的负荷预测模型
2020-05-02
来源:榕意旅游网
第38卷第5期 2016年1O月 黑龙江电力 Vo1.38 No.5 Oct.2016 Heilongjiang Electirc Power 高渗透率下分布式决策的负荷预测模型 郭国栋 (天津科技大学,天津300222) 摘要:为提高分布式新能源在电网中的渗透率和寻找有效的负荷预测方法,提出一种基于动态贝叶斯网络的分布式决策模 型,对电网负荷进行精确预测,即先用分类器按照负荷性质将大电网进行聚类划分,通过考虑影响负荷的各类因素以及这些 因素之间的关系,对各个分布式网络进行动态贝叶斯网络建模,然后用贝叶斯优化的方法以及转移网络对模型进行优化,再 用决策网络对其进行总体预测。最后,通过实际算例计算、分析结果表明,这种方法在渗透率小的一般预测时,具有较高的精 度,且学习时间短,适合在线预测。 关键词:动态贝叶斯;分布式;负荷预测;决策网络;渗透率 中图分类号:TP393 文献标志码:A 文章编号:2095—6843(2016)05—0432—07 Load forecasting model based on distributed decision—-making with high penetration GU0 Guodong (Tianjin University of Science and Technology,Tianjin 300222,China) Abstract:In order to enhance the penetration rate of distirbuted new energy in the grid and to find the effective way to forecast load,this paper proposed the distirbuted decision—making model based on dynamic Bayesian network to accurately predict the load on the grid.It is a method that ifrstly uses the classifier to divide and cluster the large grid according to the nature of the load.Then by considering all kinds of factors that influence the load and the the relationship between these factors,it establishes the model of distributed networks by dynamic Bayesian network, optimizes the model by the optimization of Bayesian method and the transfer network,and forecasts the overall load by decision—making network.Finally,the calculation and the analysis of the example shows that this method,of which the learning time is short,is effective and accurate in the system with small penetration and,which is suit— able for on—line load forecasting. Key words:dynamic Bayesian;distirbuted;load forecasting;decision—making network;penetration 智能电网是新型现代化电网,可以最大限度接 纳清洁能源,实现低碳绿色生活,提高电网安全稳 定性…。在这个框架下,随着大规模可再生能源 对大规模训练样本难以实施且存在分类问题,神经 网络算法收敛速度慢且容易失败,灰色模型误差偏 大。然而贝叶斯网络(Bayesian network,BN)为分 布式电网负荷预测提供了一种可行的方法。BN理 的并网,其对电网的渗透率不断提高,分布式电源 输出功率的波动也会在电网系统的安全性与经济 性上产生重大影响_2 J。精准的电力负荷预测可以 有效保障电力系统的安全,并且经济合理地安排 电力发电系统的启停 J,最大程度地利用可再生 能源 ]。 论是采用概率理论在网络节点上进行概率推理,可 以由一些已知节点的概率推算出另外一些节点的 概率。动态贝叶斯网络(dynamic Bayesian network, 传统的负荷预测算法包括SVM预测算法 6 J、 DBN)在普通的贝叶斯网络基础上加入了时间信 息,把原来的静态贝叶斯网络与时间信息结合,形 成具有处理时序数据的随机模型。在分布式新能 源发电负荷并网的电力系统中,风能发电的功率与 神经网络算法p J、灰色模型预测方法 J。SVM算法 收稿日期:2016—04—12。 某一时刻的风速、风向有关,光伏发电的输出功率 与太阳辐射度、环境因素有关,而一般情况下的负 荷预测与天气因素(最高温度、最低温度、降雨量 作者简介:郭国栋(1995一),男,研究方向为电力系统及自动化。 第5期 郭国栋:高渗透率下分布式决策的负荷预测模型 等)和可中断负荷因素有关。这些因素都是随机 的,且同一时刻不同因素之间与不同时刻同一因素 之间是相互影响的。基于此,本文根据电网中电源 性质的相近性利用构造贝叶斯网络分类器的方法 对整个电力系统网络进行划分,分析影响每个小网 络负荷因素之间的关系,构造动态贝叶斯预测模 型,然后探讨电网总负荷与各个子电网负荷之间的 关系,从而对整个电力网络进行精准预测,最后通 过算例计算、分析结果验证其正确性。 1电力网络的划分 1.1根据负荷性质划分 因为分布式电源的渗透率不断提高 ,对电 网来说是一个不可控源,所以在实际应用中将分布 式电源结合在微网中去。微网一般通过单点接人 电网,在电网端被视作可控发电源或负荷,微网可 以并人电网或在孤岛模式下运行,同一微网地理位 置相近,条件也比较相似,往往可以根据性质 划 分为一类电网。在电网中,商业用电、居民用电、工 厂用电是不同性质的用电,体现为不同的负荷特 性,因此在电网划分中也可以划分为几块电网。 1.2 Multi Agent技术的借鉴 划分好的电网与上级电网之间是一个分散式 的控制(如微网与上级电网之间),Multi Agent技术 是将一个大的系统划分为几个小的系统,整个系统 的预测任务可以由几个小的系统相互协作完成,并 具有以下几个方面的优势: 1)多代理技术符合电力系统分布化的趋势,能 够达到分布电源精确预测的目的。 2)除了预测之外,这种划分方式的应用有利于 上级电网与下级电网以及下级电网之间的协调 控制。 3)通过这种划分可以找到影响负荷的主要因 素,有利于制定滚动计划或调峰计划6 4)多代理系统可以加强下级电网对上级电网 的支撑,改善预测结构。借鉴多代理系统的结构, 可以根据电网负荷性质的相似性构造贝叶斯网络 分类器对整个电网进行划分。 1.3贝叶斯分类器的构造 贝叶斯分类实现的核心方法是构造贝叶斯网 络。其主要思想是:在给定的待分样本的条件下计 算类别的后验概率,确定每个小部分的后验概率。 ]/x r : 式中: [t]={ 。[t],X2[t],…, [ ]}为在t时刻属 性变量的集合,0≤ ≤T,而 t]={ 。[t], :[ ],…, [ ]}是来描述n个属性节点 。[t], [£],…, [t]在t时刻的rt个度量, [t](i=1,2,…,n)表 示属性置可能的取值;C[t]表示t时刻片段的类节 点,c [t]则表示相对应的类节点c的取值。对于类 来说,P( [t])为常数,P(置[t])表示类的先验 概率。 ’考虑到时间因素,动态贝叶斯分类器是由许多时 间片段构成的。在时间序列下,C;{c[o],C[1],…, c[ ]}为类变量序列,X={X[0],X[1],…, 71]} 则是属性变量序列,其中X[t]={ [t], [t],…, [t]},此时0≤t≤T。动态贝叶斯分类器的结构 如图1所示。 图1动态贝叶斯分类器结构 Fig.1 Structure of dynamic Bayesian classifier 通过网络分类,将各类信息归类到各个数据 库,对结构进行分析,调整划分方向。 2动态贝叶斯网络的分布式模型 2.1动态贝叶斯预测网络的构造 通过将较大的电网进行划分,可以得到多个较 小的网络。因为小网络性质相似,所以对负荷预测 更为方便,结果也更为精确。各个分布式网络的预 测遵循相似的原则。 首先,为了能对负荷进行准确预测,这个网络 应包含对分布式负荷起重要作用的关键因素 。 F( )=1一e一(∥ ) P b(0)=F( )+[1一F( 3)]= 1一e一( l ) +e一( 3n) P b(P )=F( )一F( 2)= e一(v2/A) e一(v3/A) 一P b(P +P…)/2)=F( m+1)+F( )= e一( e一(口 1/A) —式中:t)为风速,A为尺度参数,k为形状参数, 为 区间[ , ]经过离散化之后的第 段风速,F( ) 为风速的累积分布函数,P (0)为当风机的有功功 率为零时的概率,P (P )为风机的有功功率为P 时的概率,P (P +P…)/2为风机在第i段内为 黑龙江电力 第38卷 (P +P…)/2时的概率。 对于光伏来说,其负荷与光照强度G和组件温 度 有关,而组件温度又与光照强度和环境温度有 关,有如下关系: T:ro+后G 式中: 为环境温度,k为相关系数。 上述各种因素之间在同一时刻是相互影响的, 相同的因素在相邻时刻也是相互影响的,同一因素 在时间序列下具有一定的自相关性。 这里先选定风速 、光照强度G、组件温度 、环 境温度以及负荷功率P,其他因素可以用同样的方 式加人到网络中。 因为某一时刻负荷功率受到这一时刻风速、光 照强度、组件温度和环境温度的影响,同时在时序 上也与上一时刻的负荷功率有着紧密联系,因此有 如下因果关系: P(P…lP , + ,G…, + ) 式中:P为负荷功率,下标t和t+1分别表示当前时 刻和预测时刻。 组件温度、光照强度和环境温度有如下因果 关系: P( +1 l G…,ro ) 再考虑时序上的关联性,风速、组件温度、环境 温度和光照强度还有如下因果关系: P( + l ),P( + l ), P(ro +1 l ),P(G…IG ) 上述因果关系可以用局部贝叶斯网络表示,如 图2所示。 \—/ (a)P( +-1只,1,件-,G“t, “) (b)P(B+IIG t,7 。) (c)P( IV ) (d)P( 。I ) (e)P(z -l )(f)PiG ・IG ) 图2局部贝叶斯网络图 Fig.2 Diagram of local Bayesian network 将以上网络进行组合,得到静态贝叶斯网络, 如图3所示。 图3静态贝叶斯网络图 Fig.3 Diagram of static Bayesian network 再结合时间序列信息,可以得到动态贝叶斯网 络,如图4所示。 图4动态贝叶斯网络图 Fig.4 Diagram of dynamic Bayesian network 这样,动态贝叶斯网络就形成了。 2.2动态贝叶斯网络参数节点的计算 当数据比较完备时,根据EMⅢ1 算法,可以算 出参数的极大似然估计,完成参数一个周期的更 新。似然函数为 ( :D),X={ , ,…, }为动 态贝叶斯网络的变量集,置[m]是节点在第m个样 本中的数据,P [m]是父节点集的配置, 是待估参 数,D是样本集。则公式为 LiO:D)=n Pix [m ・,X m]: )= n Ⅱ e(x m]I P。 [ml:Oi)= n n pix [m]I P [m 3: )= n Li( :D) 似然函数就被分解成为17,个独立的参数极大似 然估计问题。令NiX ,P )为变量的充分统计因子, 可以得到 ( :D)=FI P( [m 3 I P [m]: = n n P(五[m] [m]:oi= 兀 n xiP(X I P Oi)州鼍 ). 第5期 郭国栋:高渗透率下分布式决策的负荷预测模型 兀 n 0剐 。虑有向边出现的次数。 的最大似然估计为 Ⅳ(X ,P。) .3决策网络的总体预测 (1) 3.1 BN决策网络的构建 .—N(P ) 2.3动态贝叶斯网络的优化 前面已经将一个大电网进行划分,又对各个划 为了能使网络在负荷预测中实时地解决问题, 实现在线预测,本文提出了贝叶斯优化(Bayesian Optimization Algorithm,BOA)与动态贝叶斯(DBN) 分的电网进行了分布式负荷预测。除了考虑各个 分布式电网外,还得结合上级电网和用户测的数 据,从而对整个电网进行一个预测,决策网络的工 网络结合的算法。 在网络图形模型智能优化的基础上,将DBN网 络作为转移网络,在整个动态优化的过程中,转移 网络的作用非常重要。当样本信息在时序轴上的t 到t+1轮优化条件变化时,转移网络会感知到所处 环境的变化,基于对当前状态的理解从而进行实时 负荷预测。其过程如图5所示。 图5 BOA与DBN的混合优化机制 Fig.5 Hyb ̄d optimization mechanism of DBN and BOA DBN网络感知到网络节点集合发生变化后,适 时改变基本条件,重新确立了新的种群以及优化方 向。大大提高了负荷智能预测的自主能力。 为了简化计算,一般往往可以去掉网络一些不 重要的节点,继承上一代转移网络的优点,重新计 量一些边。这样就需要进行网络搜索,而一般的搜 索算法往往要搜索大量的节点空间,为了提高运算 速度,采用一种启发式算法。 首先将时间 分成n个时间段,每个时间段长 度为At=T/n,将每段进行计算再综合。设定评分 函数为 SCOREsID=∑ n H P(Xf,YI Dj, )(2) 从0时刻起的△ 时间内,测得一批动态数据 D ,基于前面的BOA方法,从所有的DBN结构中搜 索到的第一个最高评分SCORE 对应网络DBN ,再 经过一个 ,得到SCORE2和DBN2,最后在nat即 时间内得到SCORE 对应的DBN 。这里所进 行的综合指的是以网络结构出现的概率为基础,考 作机制如图6所示。 图6决策网络的工作机制 Fig.6 Mechanism of decision—making network 为了精确、简洁、快速地对大电网进行总体的 预测,这里构造一个决策网络。 由于分布式网络通过感知环境后所形成的网 络节点集合的空间是很庞大的,实际应用中需要在 较短的时间内找到决策网络的结构,而以往的遗传算 法和进化算法并不具有处理这种结构的效率性,因此 可以采用改进进化算法的快速BN构图法进行寻找。 设所有变量的集合,={ , ,…, },其基本 流程如下: 1)选择初始节点 其原则是使变量取值分 别为0和1时所对应的染色体个数相差的绝对值最大。 2),in表示之前已经生成的部分树形模型中的 变量集,Io =,一,in。对于集合中的X ∈In和墨∈ ,mn,搜索到最大的BD(Xi 一 — 。 一 ),也就是寻找: BD( i 一 。 。一 )=max{BD( )I Xi∈ ,in, ∈, }。 3)把节点Xout-j加入Ii 中去,使之成为 i 一 的 子节点,这样局部就增长了一个树叶,于是 ,i ・一,i u{X。 }, 。 +一,0 u{ 。 一 } ・436・ 黑龙江电力 第38卷 4)返回2),当集合,与 相等时结束。 用上述过程构造决策网络时,为了进一步简化 运算,应按照以下原则进行:节点的独立性越弱越 先添加一个边指向它,一对节点相似性越强越先连 接它们的边,不同节点连接边的添加具有可加性, 而同一节点连接边的添加并不具有可加性。 在整个动态时序过程中,某一时刻的预测方向 荷的时间性质和空间性质,在matlab环境下做出相 对负荷在时间维度与空间维度的三维仿真图,这里 先将样本进行归一化处理。将负荷处理到[0,1]区 间内,时间维度本身就已限制在[0,24]区间内,无 需再进行处理。在空间位置上选取几何中心点,由 远及近将空间维度处理在[0,10]区间内。归一化 公式为 一 和优化点的集合都可以遗传给下一代的BN决策网 络做负荷预测使用。 3.2 BN决策网络预测的优点 1, =———一 nlin _I m“tnin 将所有数据进行处理,得到负荷的三维特性如 图7所示。 BN决策网络的构造及前文转移网络共同组成 了整个模型的双核心结构。 将大电网进行划分,再用动态贝叶斯网络进行 学习和预测,大大提高了预测的鲁棒性和精确性。 BN决策网络的使用使算法有了较快的收敛速度, 能够实现实时自主控制,即高度智能化。 BN决策网络不仅考虑了每个分布式电网,还 萎 量一。 将客户端与上级电网考虑了进来。此外,它还具有 良好的扩展性,体现在以下几个方面: 1)当决策网络与专家系统结合时,受益于决策 网络的高度智能化,将会大幅度提高电网的运行 效率。 图7负荷的三维特性 Fig.7 Three—dimensional characteristics of the load 从图7可以看出,颜色的深浅代表相对负荷的 大小,在时间维度l0到ll,空间维度4和7的位置 有2个高峰,而在时间维度20到22时,在整个空间 2)将决策网络的优化方向进行改变,则可以将 其功能扩展到负荷预测之外的广阔领域。如果将 上延绵一片小高峰,在时间维度0到4时又出现成 片蓝色低谷,在空间上也有较好的集中性。为了给 其优化方向改为改善不间断供电,那么结合新能源 发电技术和新型储能技术,可以使功率输出平滑, 提高其稳定性。此外,在故障检测和用户交互方面 也有天然的优势。 贝叶斯处理提供条件,这里采用非负矩阵分解法将 其进行降维处理,如图8所示。 4算例分析 冒 为验证本文提出的负荷预测模型,使用Matlab 捌 中的FULLBNT工具箱进行网络构建。 以某地区夏季的日负荷数据为研究对象,数据 采集时间间隔为60 min。选取天气因素相似的 2013年5月31日-20l3年5月20日的电力数据 和气候数据进行计算,来预测5月21日0—24时的 负荷,其中3个捕捉到的时间点的数据如表1所示。 表1 三个时间点的历史数据 Tablel Historical data of three time points 图8负荷特性的降维分类 Fig.8 Classiicatifon of Dimension Reduction of Load characteristis c可以利用贝叶斯分类器来处理这个二维图像 减少计算量。再利用式(2)计算各分类网络在各个 时间段内的评分。 SCORE . =∑ n兀 P(X ,Yl , ) 分析各个分布式网络与风速、温度、湿度关系 的强弱,找到强属性集G,和弱属性集G 。强属性 按照负荷性质对样本进行划分。为了了解负 集G。={A。, },弱属性集G ={A。},在要求计算 速度的前提下可以忽略A ,如表2所示。 第5期 郭国栋:高渗透率下分布式决策的负荷预测模型 表2影响因素的关联性结果 Table 2 Related results of influencing factors 在对分类后的每一部分电网中,利用本文所提 出转移矩阵的方法,以式(1)作为目标函数,利用分 布式协同进化算法将其在动态贝叶斯决策网络中 进行迭代。 决策网络中,为寻找BD(置 一 一 )= max{BD(X 置)IX ∈,in,Xi∈Io },将风速、湿度、 温度的值处理为网络点之间的概率值,经过多次学 习得到最优参数,如图9所示。 蠹 -2.0 辎一2.5 图9决策网络参数学习曲线 Fig.9 Decision network parameter learning curve 在图9中,实际运行65代就已能得到最优解和 最优参数,当前后代数参数学习极大似然值的差小 于0.001时就终止迭代了,但因为在线运行时继续 调整参数,故65代后还在继续。 通过动态贝叶斯方法可以得到各个时刻的状 态区间,即各负荷出现在某一个值上的概率。本文 采取某一时刻概率最大值所对应负荷区间的中位 数作为这一时刻的预测值。如5月21日8时的状 态区间,如图10所示。这一刻的预测值就为 270 MW 0.25 0.2O O.15 0.1O O.O5 0 O 1oo 300 40o 图1O负荷预测的概率分布 Fig.10 Probability distributions of load forecasting 将本文的算法与一般的灰色神经算法的预测 结果进行对比,如图11所示。 450 400 350 3o0 善 250 挺 2o0 一实际负荷 l5O 蕾DBN预测 l0o 咿灰色神经预测 5O O O 5 10 15 20 25 t|h 图l1负荷预测对比图 iFg.11 Load forecasting contrast curve ’ 从图11可以看到,本文所提出的模型可以有 效地跟上实际负荷的变化,而灰色神经预测模型在 中间段有明显的滞后性。 采取均方根误差作为判断负荷预测结果准确 椭酱} 嚣 性的指标,负荷预测的均方根误差曲线如图12 所示。 灰色神经预测 ・一DBN预测 图12负荷预测的均方根误差曲线 iFg.12 Mean square error curve of load forecasting 因为实际选址条件限制,本例中分布式能源的 渗透率并不高,但通过图9可以看出,DBN预测的 结果还是较为精确的,误差可控制在5%以内,符合 负荷预测的误差标准,证明基于DBN的分布式决策 网络预测方法是可行的。而灰色预测方法较为粗 糙,误差达到了10%。由此可见,在一般负荷预测 中,本文所述方法也是适用的。当新能源的渗透率 更高时,本文方法的优越性将会更加明显。为了分 析本文模型是否适用于在线预测,本文采用整个学 习过程的收敛时间作为判断指标。普通动态贝叶 斯网络与分布式决策贝叶斯网络计算性能的对比 如图13所示。 从图13可以看到,当样本数据较小时,两者的 差别并不大,但随着数据量的增大,预测所需的学 习时间有了明显的不同,这说明分布式决策比贝叶 斯更适合在线预测。 ・438・ 黑龙江电力 35(8):1829—1836. 第38卷 [5]张东霞,苗新,刘丽平,等.智能电网大数据技术发展研究[J]. 中国电机工程学报,2015,35(1):2一l2. ∞∞舳加∞∞∞∞加m O 咿普通动态贝叶斯 喈分布式决策贝叶斯 ZHANG Dongxia,MIAO Xin,LIU Liping,et a1.Research on de- velopment strategy or fsmart d big data[J].Proceedings of the CSEE,2015,35(1):2—12. [6]胡时雨,罗滇生,阳霜,等.基于多变量LS—SVM和模糊循环 推理系统的负荷预测[J].计算机应用,2015,35(2):595— 6oo. 学习数据集大d ̄/GB HU Shiyu,LUO Diansheng,YANG Shuang.et a1.L0ad forecas- ring based on multi—-variable IS——SVM and fuzzy reeursive infer- 图13算法所需时间对比 Fig.13 Comparison of computational performance of general dynamic Bayesian networks and distributed decision making Bayesina networks 5 结 论 1)利用动态贝叶斯网络方法,考虑时序因素, 充分利用历史大数据进行学习,能提高预测精度。 2)按照负荷性质,通过分类器对大电网进行聚 类划分,并进行分布式计算,能有效降低了网络的 复杂度。 3)提出了决策网络方法,结合分布式网络,可 大大降低计算复杂度。而且决策网络扩展性好,可 以用于电网优化、故障检测等方面。 4)这种方法在渗透率较小的一般预测时,具有 较高的精度,且学习时间短,适合在线预测。 参考文献 [1]张东霞,姚良忠,马文媛.中外智能电网发展战略[J].中国 电机工程学报,2013,33(31):1—14. ZHANG Dongxia,YAO Liangzhong,MA Wenyuan.Development strategies of smart d in China and abroad[J].Proceedings of the CSEE,2013,33(31):1—14. [2]鲍薇,胡学浩,何国庆,等.分布式电源并网标准研究[J]. 电网技术,2012,36(11):46—52. BAO Wei,HU Xuehao,HE Guoqing,et a1.Study on standard f0r grid—integration of distirbuted resources[J].Power System Technology,2012,36(11):46—52. [3]黄海涛,胡学英,李翔,等.实用化的激励性可中断负荷最优 补偿定价模型[J].电网技术,2014,38(8):2149—2154. HUANG Haitao,HU Xueying,LI Xiang,et a1.A practicable。p— timal compensation pricing model for incentive interruptible load contract[J].Power System Technology,2014,38(8):2149— 2154. [4]王继业,季知祥,史梦洁,等.智能配用电大数据需求分析与 应用研究[J].中国电机工程学报,2015,35(8):1829— 1836. WANG Jiye,J1 Zhixiang,SHI Mengjie,et a1.Scenario analysis and application research on big data in smart power distlibution and consumption systems[J].Proceedings of the CSEE,2015. ence system[J].Journal of Computer Applications,2015,35(2): 595—600. [7]庞清乐.基于粗糙集理论的神经网络预测算法及其在短期负 荷预测中的应用[J].电网技术,2010,34(12):168—173. PANG Qingle.A rough set—based neurla network load forecasting algoirthm and its application in sh0n—term load forecasting[J]. Power System Technology,2010,34(12):168—173. [8]葛少云,贾鸥莎,刘洪.基于遗传灰色神经网络模型的实时 电价条件下短期电力负荷预测[J].电网技术,2012,36(1): 224—229. GE Shaoyun,JIA Ousha,LIU Hong.A rgay neural network model improved by genetic algorithm for short—Term load forecasting in price—sensitive environment[J].Power System Technology, 2012,36(1):224—229. [9]王深哲,高山,李海峰,等.含高渗透率间歇性电源的电网规 划方案优选方法[J].电网技术,2013,37(8):2129—2135. WANG Shenzhe,GAO Shan,LI Haifeng,et a1.Optimal selection of planning schemes for power gdd with high penetration of inter- mittent generation[J],Power System Technology,2013,37(8): 2129—2135. [10]KERIN U.Dynamic security indication in power systems with large amount of renewanles[C]//Power and Energy Society General Meeting.San Diego:IEEE,2012:1—6. [11]施琳,罗毅,施念,等.高渗透率风电一储能孤立电网控制策 略[J].中国电机工程学报,2013,33(16):78—85. SHI Lin,LUO Yi,SHI Nian,et a1.A contorl strategy of isolated grid with high penetration of wind nad energy storage systems[J]. Prcoeedings ofthe CSEE,2013,33(16):78—85. [12]田书,王丽.考虑负荷类型的含分布式电源系统对配电网电 压的影响[J].电网与清洁能源,2014,30(2):43—47. TIAN Shu,WANG Li.Influence of DG on the distribution net. work voltgae considering the load type[J].Power System and Clean Eneryg,2014,30(2):43—47. [13]高元海,王淳.基于全概率公式的含风电配电系统概率潮流 计算[J].中国电机工程学报,2015,35(2):327—334. GAO Yuanhai,WANG Chun.Probabilistic load flow calculation of distribution system including wind farms based on total probabil— ity formula[J].Proceedings ofthe CSEE,2015,35(2):327—334. [14]俞奎,王浩,吴信东,等.贝叶斯网络的并行EM学习算 法[J].模式识别与人工智能,2008,21(5):670—675. YU Kui,WANG Hao,WU Xindong,et a1.Leaning Bayesina net— works using a parlalel EM apporach[J].PR&A1,2008,21(5): 670—675. (责任编辑张兴业)