Page 85 - 中国仿真学会通讯
P. 85
量的百分比［１５］ ， 可表示为
ｍ
ζＰＶ％ ＝ ｎ ＰＮ， ＰＶ， ｉ （６）
ｉ＝１
ｎ
ＰＮ， ＰＶ， ｉ ＋ ＰＮ， Ｇ， ｊ
ｉ＝１ ｊ＝１
式中： ｍ、 ｎ 分别为光伏发电机组和传统发电机组的个数； ＰＮ，ＰＶ，ｉ为第 ｉ 台光伏发电机组的额定容
量； ＰＮ，Ｇ，ｊ为第 ｊ 台传统发电机组的额定容量。
１ ４ 基于非序贯蒙特卡罗模拟法的电力系统可靠性评估过程
（１） 建立多级水平负荷模型， 对每一级负荷水平进行系统状态抽样。
（２） 利用蒙特卡洛模拟选择系统状态。
１） 利用式（７） 所示的正态分布随机变量以计及母线负荷的不确定性。
Ｍσｉｊ ＝ （ ｚｉｊ ＋ １） Ｍｉｊ （７）
式中： Ｍσｉｊ为遵循曲线 ｊ 的母线负荷在第 ｉ 级负荷水平下的抽样值； Ｍｉｊ 为聚类法建立的多级负荷
模型中曲线 ｊ 的第 ｉ 级负荷水平， Ｍｉｊ和 Ｍσｉｊ都要用峰荷为基准的标幺值表示； ｚｉｊ为对应与曲线 ｊ 的
母线负荷在第 ｉ 级负荷水平的标准正态分布随机数。
遵循负荷曲线 ｊ 的母线负荷抽样值由 Ｍσｉｊ乘以母线峰荷值而得。 如果要考虑母线负荷之间的
相关性， 则应使用相关抽样法进行抽样［１８］ 。
２） 元件状态（运行、 停运或降额） 利用均匀分布随机变量来模拟。
ｓｋ ＝ ïïíì１０，， 如果 Ｒｋ ＞ ｐＰｋ ＋ ｐＦｋ（ 运行） （８）
ïïî２， 如果 ｐＰｋ ≤ Ｒｋ ≤ ｐＦｋ（ 停运）
如果 ０ ≤ Ｒｋ ≤ ｐＰｋ（降额）
式中： ｓｋ 为元件 ｋ 的抽样状态； Ｒｋ 为第 ｋ 个元件在［０， １］ 区间的均匀分布随机数； ｐＦｋ、 ｐＰｋ 分别
为第 ｋ 个元件处于停运和降额状态的概率［１８］ 。
（３） 进行预想故障分析， 如果必要的话， 则还需要进行最优潮流分析， 估计所选择的状态
需要削减的负荷量。 如果负荷削减不为 ０， 则这个被选择了的状态是一个失效状态， 记录系统失
效状态的负荷削减量。
（４） 计算风险指标。
１） 负荷削减频率 σＰＬＣ 。
σＰＬＣ ＝ ＮＬ êêéë ｓ∈Ｆｉ ｎ（ ｓ） úúûù Ｔｉ （９）
ｉ＝１ Ｎｉ Ｔ
式中： ＮＬ 为负荷水平分级数； Ｆｉ 为多级负荷模型中第 ｉ 级负荷水平下系统失效状态的集合； ｎ
（ｓ） 为抽样中 ｓ 状态的发生数； Ｎｉ 为抽样总数； Ｔｉ 为第 ｉ 级负荷水平的时间长度； Ｔ 为以 ３６５ 天
计算的 １ 年的时间长度。
２） 期望缺供电量 ＥＥＮＳ。
ＥＥＮＳ ＝ ＮＬ êëéê ｓ∈Ｆｉ ｎ（ ｓ） Ｃ（ ｓ） úùúû Ｔｉ （ １０）
ｉ＝１ Ｎｉ
式中 Ｃ（ｓ）为状态 ｓ 的负荷削减量。
为了更好地研究 ＣＨＰ 渗透率对 ＩＣＥＳ 可靠性的影响， 本文定义， 在研究 ＩＣＥＳ 中 ＣＨＰ 渗透率
发生变化的时候， 利用 ＥＥＮＳ变化率 ΔＥＥＮＳ来描述系统可靠性指标的变化情况。
ΔＥＥＮＳ ＝ ＥＥＮＳ × １００％ （ １１）
Ｅ ＥＥＮＳ０
８１
