Page 85 - 中国仿真学会通讯

量的百分比［１５］ ， 可表示为

ｍ

ζＰＶ％ ＝ ｎ ＰＮ， ＰＶ， ｉ （６）

ｉ＝１
ｎ

ＰＮ， ＰＶ， ｉ ＋ ＰＮ， Ｇ， ｊ
ｉ＝１ ｊ＝１

式中： ｍ、 ｎ 分别为光伏发电机组和传统发电机组的个数； ＰＮ，ＰＶ，ｉ为第 ｉ 台光伏发电机组的额定容

量； ＰＮ，Ｇ，ｊ为第 ｊ 台传统发电机组的额定容量。

１ ４ 基于非序贯蒙特卡罗模拟法的电力系统可靠性评估过程

（１） 建立多级水平负荷模型， 对每一级负荷水平进行系统状态抽样。

（２） 利用蒙特卡洛模拟选择系统状态。

１） 利用式（７） 所示的正态分布随机变量以计及母线负荷的不确定性。

Ｍσｉｊ ＝ （ ｚｉｊ ＋ １） Ｍｉｊ （７）

式中： Ｍσｉｊ为遵循曲线 ｊ 的母线负荷在第 ｉ 级负荷水平下的抽样值； Ｍｉｊ 为聚类法建立的多级负荷

模型中曲线 ｊ 的第 ｉ 级负荷水平， Ｍｉｊ和 Ｍσｉｊ都要用峰荷为基准的标幺值表示； ｚｉｊ为对应与曲线 ｊ 的

母线负荷在第 ｉ 级负荷水平的标准正态分布随机数。

遵循负荷曲线 ｊ 的母线负荷抽样值由 Ｍσｉｊ乘以母线峰荷值而得。 如果要考虑母线负荷之间的
相关性， 则应使用相关抽样法进行抽样［１８］ 。

２） 元件状态（运行、 停运或降额） 利用均匀分布随机变量来模拟。

ｓｋ ＝ ïïíì１０，， 如果 Ｒｋ ＞ ｐＰｋ ＋ ｐＦｋ（ 运行） （８）
ïïî２， 如果 ｐＰｋ ≤ Ｒｋ ≤ ｐＦｋ（ 停运）
如果 ０ ≤ Ｒｋ ≤ ｐＰｋ（降额）

式中： ｓｋ 为元件 ｋ 的抽样状态； Ｒｋ 为第 ｋ 个元件在［０， １］ 区间的均匀分布随机数； ｐＦｋ、 ｐＰｋ 分别

为第 ｋ 个元件处于停运和降额状态的概率［１８］ 。

（３） 进行预想故障分析， 如果必要的话， 则还需要进行最优潮流分析， 估计所选择的状态

需要削减的负荷量。 如果负荷削减不为 ０， 则这个被选择了的状态是一个失效状态， 记录系统失

效状态的负荷削减量。

（４） 计算风险指标。

１） 负荷削减频率 σＰＬＣ 。

σＰＬＣ ＝ ＮＬ êêéë ｓ∈Ｆｉ ｎ（ ｓ） úúûù Ｔｉ （９）
ｉ＝１ Ｎｉ Ｔ

式中： ＮＬ 为负荷水平分级数； Ｆｉ 为多级负荷模型中第 ｉ 级负荷水平下系统失效状态的集合； ｎ

（ｓ） 为抽样中 ｓ 状态的发生数； Ｎｉ 为抽样总数； Ｔｉ 为第 ｉ 级负荷水平的时间长度； Ｔ 为以 ３６５ 天

计算的 １ 年的时间长度。

２） 期望缺供电量 ＥＥＮＳ。

ＥＥＮＳ ＝ ＮＬ êëéê ｓ∈Ｆｉ ｎ（ ｓ） Ｃ（ ｓ） úùúû Ｔｉ （ １０）
ｉ＝１ Ｎｉ

式中 Ｃ（ｓ）为状态 ｓ 的负荷削减量。

为了更好地研究 ＣＨＰ 渗透率对 ＩＣＥＳ 可靠性的影响， 本文定义， 在研究 ＩＣＥＳ 中 ＣＨＰ 渗透率

发生变化的时候， 利用 ＥＥＮＳ变化率 ΔＥＥＮＳ来描述系统可靠性指标的变化情况。

ΔＥＥＮＳ ＝ ＥＥＮＳ × １００％ （ １１）
Ｅ ＥＥＮＳ０

８１