考虑供暖约束的时序蒙特卡洛模拟可靠性评估方法及系统的制作方法
【专利摘要】本公开涉及考虑供暖约束的时序蒙特卡洛模拟可靠性评估方法及系统，所述方法在传统时序蒙特卡洛模拟可靠性评估的框架下，对供暖约束进行了建模，在可靠性评估的框架中考虑了冬季供暖以及风电的交互作用对电力系统可靠性的影响，所述方法主要是通过考虑那些承担供暖约束的热电机组的出力约束条件和对优先接纳风电原则的合理建模来反映冬季供暖对电力系统可靠性的影响。所述系统基于所述方法实现。本公开能够更合理、更准确地对我国北方地区的可靠性水平进行评估。
【专利说明】
考虑供暖约束的时序蒙特卡洛模拟可靠性评估方法及系统
技术领域
[0001] 本公开设及电力系统规划评估领域，具体设及一种考虑了我国冬季供暖影响的使 用时序蒙特卡洛模拟进行可靠性评估方法及系统。
【背景技术】
[0002] 由于我国北方地区冬季需要供暖，而我国一般实行热定电"的方式，运样就要 求电力供应必须满足一定的要求，才能保证系统满足供暖的需求。目前热力负荷一般是由 火电机组承担，运样一般运些机组出力必须保证在一定范围才能保证供暖的影响。另一方 面，我国北方地区风电装机容量已初见规模，然而风电出力与热力需求有很强的时空相关 性，冬季常常是风电出力较大的时候，而供热的约束，导致系统调峰容量减小，运进一步导 致了冬季弃风事件的出现。传统的可靠性评估并未考虑运些影响，并未考虑供暖对系统可 靠性的影响。
[0003] 更加合理的可靠性评估方法，可W发现和识别与系统可靠性关系最为密切的环 节，更加准确的识别系统的薄弱环节，为保证系统的稳定运行和提高系统的可靠性提供合 理的建议。
[0004] 现有刻画系统可靠性最为准确的评估框架是时序蒙特卡洛可靠性评估方法，该框 架能够考刻画新能源出力的波动性，因此在含新能源区域电网的可靠性评估中得到了广泛 的应用。但是还未有可靠性评估框架考虑供暖约束的影响，W及供暖与风电之间的交互影 响。
【发明内容】
[0005] 针对上述问题，本公开提供了考虑供暖约束的时序蒙特卡洛模拟可靠性评估方法 及系统，所述方法及系统能够更合理、更准确地对我国北方地区的可靠性水平进行评估。
[0006] -种考虑供暖约束的时序蒙特卡洛模拟可靠性评估方法，所述方法包括下述步 骤：
[0007] S100、初始化元件状态:输入电力系统数据，进行初始潮流和初始连通性计算，设 定所有元件状态正常；
[000引S200、获取电力系统状态的持续时间Tk:对电力系统中每个元件的状态持续时间 进行抽样，得到电力系统状态的时间序列耗}/二耗，/2，...，^^+,,。}/，进而得到电 力系统状态的持续时间Tk,化二
[0009] S300、安排机组出力：利用S300中得到的Tk更新当前时间；如果当前属于供暖期， 先安排热电机组出力保证供热约束，然后再接纳风电机组和/或光伏机组的出力；否则按照 优先接纳风电机组和/或光伏机组的出力，在安排火电机组出力；
[0010] S400、进行状态评估:对化时间段内的电力系统状态进行可靠性评估。
[0011] 基于所述方法，实现了一种考虑供暖约束的时序蒙特卡洛模拟可靠性评估系统， 所述系统包括下述模块：
[0012] MlOO、元件状态初始化模块:输入电力系统数据，进行初始潮流和初始连通性计 算，设定所有元件状态正常；
[0013] M200、电力系统状态的持续时间化获取模块:对电力系统中每个元件的状态持续 时间进行抽样，得到电力系统状态的时间序列Ui}/二Ui，^2，一.，f,,t+?e}/，进而得到 电力系统状态的持续时间Tk，化=min {ti} k，输出给M300和M400;
[0014] M300、机组出力安排模块:利用M200中得到的化更新当前时间；如果当前属于供暖 期，先安排热电机组出力保证供热约束，然后再接纳风电机组和/或光伏机组的出力；否则 按照优先接纳风电机组和/或光伏机组的出力，在安排火电机组出力；
[0015] M400、状态评估模块:对化时间段内的电力系统状态进行可靠性评估。
【附图说明】
[0016] 图1是本公开一个实施例中的考虑供热约束的时序蒙特卡洛可靠性评估框架示意 图；
[0017] 图2是本公开一个实施例中的机组出力安排示意图；
[0018] 图3是本公开一个实施例中的时序蒙特卡洛模拟法采样示意图。
【具体实施方式】
[0019] 下面结合附图1-3对其进行说明，在一个基础实施例中，提供了一种考虑供暖约束 的时序蒙特卡洛模拟可靠性评估方法，所述方法包括下述步骤：
[0020] S100、初始化元件状态:输入电力系统数据，进行初始潮流和初始连通性计算，设 定所有元件状态正常；
[0021] S200、获取电力系统状态的持续时间化:对电力系统中每个元件的状态持续时间 进行抽样，得到电力系统状态的时间序列= Ul,h，进而得到电 力系统状态的持续时间化，化=min{ti}k;
[0022] S300、安排机组出力：利用S300中得到的化更新当前时间；如果当前属于供暖期， 先安排热电机组出力保证供热约束，然后再接纳风电机组和/或光伏机组的出力；否则按照 优先接纳风电机组和/或光伏机组的出力，在安排火电机组出力；
[0023] S400、进行状态评估:对化时间段内的电力系统状态进行可靠性评估。
[0024] 在运个实施例中，所述输入的电力系统数据包括：
[0025] (1)该电力系统各节点年度每小时的负荷水平，W及基于该数据建立的时序负荷 模型；
[0026] (2)该电力系统的线路容量、线路原始强迫停运率；
[0027] (3)该电力系统的各种发电机组的发电额定容量、发电机组最小出力等。
[0028] 在运个实施例中，S100中在进行初始潮流和初始连通性计算时，若计算结果为联 通，则执行S200，否则检查输入数据。
[0029] 所述S300中，所述风电机组和/或光伏机组的出力的数据来源于历史数据，即某一 年的风电和/或光伏机组的出力数据，那么在第k个电力系统状态下的风风电机组和光伏发 电机组的输出功率，即该年份中（Time，Time巧k)时段风电机组和光伏机组的出力，其中 Time为模拟的电力系统所处的一年中的时间。若电力系统此时工作的时间段(Τ ime，Time+ Tk)属于供暖期，如图2安排机组的出力，如果属于供暖期，先安排热电机组出力保证供热约 束，然后再保证优先接纳风电；否则按照优先接纳风电原则，安排机组出力。比如，一年有 8736小时，供热期的时段属于(Thi，化2)，如果时间属于运一区间即化i<Time<&2,则认为该时 段属于供暖期。
[0030]在运个实施例中，所述S400中对化时间段内的电力系统状态进行可靠性评估，所 述可靠性评估内容包括化时间段内总停电时间Tloss，Tk时间段内总损失电量Eloss，年度停电 时间TYloss，年度停电损失电量EYloss。若化时间段内总停电时间Tloss越短，化时间段内总损失 电量Eims越少，则电力系统在化时间段内可靠性越强;若年度停电时间TYiMs越短，年度停电 损失电量EYIdss越少，则电力系统可靠性越强，该方法可W通过比较可靠性指标的大小衡量 系统考虑供热约束后对于系统可靠性的影响。如果在T时刻，电力系统处于第k个状态，该状 态电力系统切负荷A Pk，那么电力系统从第1个状态到第k个状态总停电时间TiDss，总损失电 量EIdss，年度停电时间TYlDSS，年度停电损失电量EYIdss按照下面的式子计算：
[0035] 在一个实施例中，给出了在考虑供暖约束情况下的电力系统切负荷APk的求解方 法。如图1所示的状态评估部分，所述状态评估主要是对抽样的系统状态进行可靠性评估， 研究其是否会因为元件的停运而导致系统失负荷。具体流程包括：
[0036] 对在T时刻，处于第k状态的电力系统的网络拓扑分析，判断电力系统是否解列；
[0037] S100、若电力系统存在解列，则对电力系统各个解列块进行功率平衡分析，统计各 解列块内的失负荷量；
[0038] S200、若没有解列，则判断电力系统功率是否平衡；
[0039] S201、若电力系统功率平衡，则进入潮流计算;否则，执行步骤S203;
[0040] 所述潮流计算首先判断系统是否有孤立节点，有孤立节点去除孤立节点再形成导 纳矩阵，然后进行线路潮流计算；
[0041] S202、通过潮流计算结果判断电力系统是否存在线路过载;若无过负荷，则结束评 估;否则，执行步骤S203;
[0042] S203、进入基于直流潮流的线性规划模型进行最小切负荷Δ Pk计算：
[0043]
[0044] 式中；
[0045] Nd--负荷节点数目；
[0046] r--节点削减负荷量；
[0047] PgWin一一发电机最小有功出力；
[0048 ] Pg--发电机有功出力向量；
[0049] Pgmax一一发电机最大有功出力向量.
[0050] Pi--支路潮流向量；
[0化1 ] If*一一直流潮流极限有功功率向量；
[0化2] Pd--节点负荷向量；
[0化3] B一一节点导纳矩阵；
[0化4] Θ--节点相角向量。
[0055]其中，新能源机组约束如下：
[0化6]
[0057]式中；
[0化引 Pwind--风电机组此时的实际出力；
[0059] Psolar--光伏机组此时的实际出力；
[0060] 上面的约束表明新能源机组的出力在发电再调度和负荷调整模型中保持不变，不 发生调整，即新能源机组的出力不可调。
[0061] 相较于传统的可靠性评估，运里的可靠性评估过程更多地考虑了供热机组的约 束，其中最重要的是下式中不等式约束的下限，该约束考虑了若要保证当地居民的供暖需 求，必须要提供的最小火电出力。
[0062] 其中，供热机组约束如下：
[0063]
[0064] 式中；
[0065] 扣eat-一承担供热的机组此时的实际出力；
[0066] -一承担供热的机组出力的下限，其受到供热约束的影响，必须保证基本的 heat 供热需求得到满足；
[0067] --承担供热的机组出力的上限。
[0068] 优选的，所述S200中所述状态持续时间服从指数分布。
[0069] 优选的，按照时序蒙特卡洛进行抽样，其采样全过程如图3所示，因此所述S200中 所述电力系统状态的时间序列化}/二扣1，重2，.*.^,.+/,。}/满足：
[0070]
[0071] 式中；
[0072] 分别表示序列和Hi}冲的元素；
[007引山-康示第k-1时段的持续时间Tk-i。
[0074] 对于初始状态W及某个电力系统状态中需要抽样产生新的状态持续时间的情况， 在一个实施例中给出了通过蒙特卡洛生成电力系统的状态时间的方法，即按照下式生成：
[0075]
[0076] 式中；
[0077] MTTF，MTTR-一分别表示平均无故障时间，平均修复时间。
[0078] rand( ·) -一表示0-1之间的随机数。
[0079] 优选的，在S400之后，所述方法还包括下述步骤：
[0080] S500、重复S200~S400W通过对该电力系统状态的进行多年的可靠性评估来判断 该电力系统的可靠性。
[0081] 通过模拟N年的电力系统状态，获取电力系统状态总的持续时间T，总停电时间 TlDss，总损失电量Eidss，年度停电时间TYlDss，年度停电损失电量EYIdss，进而获得一个该电力 系统的年度可靠性指标，所述年度可靠性指标包括年电力不足概率LOLPn与年电力电量不 足期望值EENSn，进而可W获得该电力系统一个趋于常数的电力不足概率L0LP与电力电量 不足期望值EENS，运两个值可W很好的反映该电力系统在考虑供暖情况下的运行可靠性。 其中：
[0082]
[0083] EENSn=EYi〇ss
[0087] 进一步地，所述S500中所述电力不足概率的获得包括下述步骤：
[0088] S501、按照下式计算连续若干年的电力电量不足期望构成的数列的收敛方差系 数：
[0089]
[0090] 式中；
[0091] var化ENS(1:N))代表求第1年到第N年的电力电量不足期望值邸NS的方差，mean代 表求第1年到第N年的电力电量不足期望值EENS的均值；
[0092] S502、若化ENS<ft?AX，则判定已获取该电力系统的电力不足概率;其中，βΜΑΧ为设定 值。
[0093] 在运里，βΜΑΧ设置得越小，获取的电力不足概率越能反映于该电力系统运行时的状 态。
[0094] 更进一步地，若所模拟的电力系统在设定的收敛条件下不收敛，还可W设定最大 模拟的年数为Nmax。
[00M]基于所述方法，在一个实施例中给出了一种考虑供暖约束的时序蒙特卡洛模拟可 靠性评估系统，所述系统包括下述模块：
[0096] M100、元件状态初始化模块:输入电力系统数据，进行初始潮流和初始连通性计 算，设定所有元件状态正常；
[0097] M200、电力系统状态的持续时间化获取模块:对电力系统中每个元件的状态持续 时间进行抽样，得到电力系统状态的时间序列扭i }:/二托，，…，+,:c }/，进而得 到电力系统状态的持续时间Tk,化=min {ti} k，输出给M300和M400;
[0098] M300、机组出力安排模块:利用M200中得到的化更新当前时间；如果当前属于供暖 期，先安排热电机组出力保证供热约束，然后再接纳风电机组和/或光伏机组的出力；否则 按照优先接纳风电机组和/或光伏机组的出力，在安排火电机组出力；
[0099] M400、状态评估模块:对化时间段内的电力系统状态进行可靠性评估。
[0100] 在运个实施例中，所述输入的电力系统数据包括：
[0101] (1)该电力系统各节点年度每小时的负荷水平，W及基于该数据建立的时序负荷 模型；
[0102] (2)该电力系统的线路容量、线路原始强迫停运率；
[0103] (3)该电力系统的各种发电机组的发电额定容量、发电机组最小出力等。
[0104] 在运个实施例中，在进行初始潮流和初始连通性计算时，若计算结果为联通，则继 续计算;否则，检查输入数据。所述风电机组和/或光伏机组的出力的数据来源于历史数据， 即某一年的风电和/或光伏机组的出力数据，那么在第k个电力系统状态下的风风电机组和 光伏发电机组的输出功率，即该年份中（Time，Time巧k)时段风电机组和光伏机组的出力， 其中Time为模拟的电力系统所处的一年中的时间。若电力系统此时工作的时间段(Time, Time巧k)属于供暖期，如图2安排机组的出力，如果属于供暖期，先安排热电机组出力保证 供热约束，然后再保证优先接纳风电；否则按照优先接纳风电原则，安排机组出力。比如，一 年有8736小时，供热期的时段属于(Trl，Tr2)，如果时间属于运一区间即Trl<Time<Tr2,则认为 该时段属于供暖期。M400中所述对化时间段内的电力系统状态进行可靠性评估，所述可靠 性评估内容包括化时间段内总停电时间Tloss，Tk时间段内总损失电量Eloss，年度停电时间 了￥1。33，年度停电损失电量6￥1。33。若化时间段内总停电时间1'1。33越短，了1^寸间段内总损失电量 Eloss越少，则电力系统在化时间段内可靠性越强;若年度停电时间TYlMS越短，年度停电损失 电量EYIdss越少，则电力系统可靠性越强。如果在T时刻，系统处于第k个状态，该状态电力系 统切负荷A P，则所述总停电时间Tioss，总损失电量Eioss，年度停电时间TYioss，年度停电损失 电量EYioss按照下面的式子计算：
[0109] 优选的，所述M200中所述状态持续时间服从指数分布。
[0110] 优选的，12〇0中所述电力系统状态的时间序列捧}/二{《1，12，...，^,,,,+。<,}/ 满足：
[0111] 式中：
[0112]
[011引ti,k-i，ti,k分别表示序列{ti}k-i和Hi}冲的元素；
[0114] 山-康示第k-1时段的持续时间Tk-1。
[0115] 可选的，对于ti,k-1 = dk-i或者k=l的情形，在一个实施例中给出了通过蒙特卡洛生 成电力系统的状态时间的方法，即按照下式生成：
[0116]
[0117]式中；
[011引 MTTF，MTTR-一分别表示平均无故障时间，平均修复时间。
[0119] rand( ·)--表示0-1之间的随机数。
[0120] 优选的，所述系统还包括下述模块：
[0121] M500、电力不足概率获取模块:循环调用M200~M400 W通过对该电力系统状态的 进行多年的可靠性评估来判断该电力系统的可靠性。
[0122] 通过模拟N年的电力系统状态，获取电力系统状态总的持续时间T，总停电时间 TlDss，总损失电量EIdss，年度停电时间TYlDSS，年度停电损失电量EYIdss，进而获得一个该电力 系统的年度可靠性指标，所述年度可靠性指标包括年电力不足概率LOLPn与年电力电量不 足期望值EENSn，进而可W获得该电力系统一个趋于常数的电力不足概率L0LP与电力电量 不足期望值EENS，运两个值可W很好的反映该电力系统在考虑供暖情况下的运行可靠性。 其中：
[0128] 更进一步地，若所模拟的电力系统在设定的收敛条件下不收敛，还可W在可靠性 评估系统中设定最大模拟的年数为Nmax。
[0129] 进一步地，所述M500中所述电力不足概率通过下述单元获得：
[0130] U501、计算单元:按照下式计算连续若干年的电力电量不足期望构成的数列的收 敛方差系数，并输出给判断单元U502:
[0131]
[0132] 式中；
[0133] var化ENS(1:N))代表求第1年到第N年的电力电量不足期望值邸NS的方差，mean代 表求第1年到第N年的电力电量不足期望值EENS的均值；
[0134] 肌02、判断单元：若feENS<eMAX，则判定已获取该电力系统的电力不足概率；其中， fiwAX为设定值。
[0135] 综上，本公开的方法及系统能够更合理、更准确地对我国北方地区的可靠性水平 进行评估。
[0136] W上对本公开进行了详细介绍，本文中应用了具体个例对本公开的原理及实施方 式进行了阐述，W上实施例的说明只是用于帮助理解本公开的方法及其核屯、思想；同时，对 于本领域技术人员，依据本公开的思想，在【具体实施方式】及应用范围上均会有改变之处，综 上所述，本说明书内容不应理解为对本公开的限制。
【主权项】
1. 一种考虑供暖约束的时序蒙特卡洛模拟可靠性评估方法，其特征在于，所述方法包 括下述步骤： S100、初始化元件状态:输入电力系统数据，进行初始潮流和初始连通性计算，设定所 有元件状态正常； S200、获取电力系统状态的持续时间Tk:对电力系统中每个元件的状态持续时间进行抽 样，得到该电力系统状态的时间序列，进而得到该电力 系统状态的持续时间Tk，Tk=min{ti}k; S300、安排机组出力:利用S300中得到的Tk更新当前时间；如果当前属于供暖期，先安排 热电机组出力保证供热约束，然后再接纳风电机组和/或光伏机组的出力；否则按照优先接 纳风电机组和/或光伏机组的出力，在安排火电机组出力； S400、进行状态评估:对Tk时间段内的电力系统状态进行可靠性评估。2. 根据权利要求1所述的方法，其特征在于： 优选的，所述S200中所述状态持续时间服从指数分布。3. 根据权利要求1所述的方法，其特征在于，所述S200中所述电力系统状态的时间序列 鉍}/ = "...，人￡+%}/满足：式中： ti,k-1，ti,k分别表不序列{ti}k-1和{ti}k中的兀素； 表示第k-Ι时段的持续时间Τη。4. 根据权利要求1所述的方法，其特征在于，在S400之后，所述方法还包括下述步骤： 5500、 重复S200~S400以通过对该电力系统状态进行多年的可靠性评估来判断该电力 系统的可靠性。5. 根据权利要求4所述的方法，其特征在于，所述S500中所述电力不足概率的获得包括 下述步骤： 5501、 按照下式计算连续若干年的电力电量不足期望构成的数列的收敛方差系数：式中： var(EENS(l:N))代表求第1年到第Ν年的电力电量不足期望值EENS的方差，mean代表求 第1年到第N年的电力电量不足期望值EENS的均值； 5502、 若&ENS <βΜΑΧ，则判定已获取该电力系统的电力不足概率;其中，βΜΑΧ为设定值。6. -种考虑供暖约束的时序蒙特卡洛模拟可靠性评估系统，其特征在于，所述系统包 括下述模块： Μ100、元件状态初始化模块:输入电力系统数据，进行初始潮流和初始连通性计算，设 定所有元件状态正常； M200、电力系统状态的持续时间Tk获取模块:对电力系统中每个元件的状态持续时间进 行抽样，得到电力系统状态的时间序列Uih = 进而得到电力 系统状态的持续时间Tk，Tk=min {ti} k，输出给M300和M400; M300、机组出力安排模块:利用M200中得到的Tk更新当前时间；如果当前属于供暖期，先 安排热电机组出力保证供热约束，然后再接纳风电机组和/或光伏机组的出力；否则按照优 先接纳风电机组和/或光伏机组的出力，在安排火电机组出力； M400、状态评估模块:对Tk时间段内的电力系统状态进行可靠性评估。7. 根据权利要求6所述的系统，其特征在于： 所述M200中所述状态持续时间服从指数分布。8. 根据权利要求6所述的系统，其特征在于，M200中所述电力系统状态的时间序列 (广i } /. _ V1 .，# 2，·. · .，（/满足： 式中：ti,k-i，ti,k分别表不序列{ti}k-1和{ti}k中的兀素； 表示第k-1时段的持续时间IW。9. 根据权利要求6所述的系统，其特征在于，所述系统还包括下述模块： M500、电力不足概率获取模块:循环调用M200~M400以通过对该电力系统状态进行多 年的可靠性评估来判断该电力系统的可靠性。10. 根据权利要求9所述的系统，其特征在于，所述M500中所述电力不足概率通过下述 单兀获得：U501、计算单元:按照下式计算连续若干年的电力电量不足期望构成的数列的收敛方 差系数，并输出给判断单元 式中： var(EENS(l:N))代表求第1年到第N年的电力电量不足期望值EENS的方差，mean代表求 第1年到第N年的电力电量不足期望值EENS的均值； U502、判断单元:若&ENS<βΜΑΧ，则判定已获取该电力系统的电力不足概率;其中，ftlAX为 设定值。
【文档编号】H02J3/00GK105870913SQ201610170230
【公开日】2016年8月17日
【申请日】2016年3月23日
【发明人】刘慧超, 严超, 别朝红, 李更丰, 米渊, 杨扬
【申请人】国网山西省电力公司大同供电公司, 西安交通大学