CN106849176B - 一种计及风电不确定性的随机序列优化模型的建模方法 - Google Patents
一种计及风电不确定性的随机序列优化模型的建模方法 Download PDFInfo
- Publication number
- CN106849176B CN106849176B CN201710229367.4A CN201710229367A CN106849176B CN 106849176 B CN106849176 B CN 106849176B CN 201710229367 A CN201710229367 A CN 201710229367A CN 106849176 B CN106849176 B CN 106849176B
- Authority
- CN
- China
- Prior art keywords
- wind power
- output
- line
- representing
- probability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 238000005457 optimization Methods 0.000 title claims abstract description 13
- 238000009826 distribution Methods 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- UNIDAFCQFPGYJJ-UHFFFAOYSA-N 2-amino-2-(2-chloro-5-hydroxyphenyl)acetic acid Chemical compound OC(=O)C(N)C1=CC(O)=CC=C1Cl UNIDAFCQFPGYJJ-UHFFFAOYSA-N 0.000 claims description 2
- 230000010354 integration Effects 0.000 abstract description 3
- 238000010438 heat treatment Methods 0.000 description 5
- 230000033228 biological regulation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000010248 power generation Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 2
- 239000003245 coal Substances 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
Images
Classifications
-
- H02J3/386—
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种计及风电不确定性的随机序列优化模型的建模方法,该建模方法包括以下步骤:利用概率性序列描述风电出力的不确定性;建立基于概率性序列的线路潮流的约束模型;建立基于概率性序列的系统旋转储备约束模型。该模型可以很好地处理风电的不确定性,并可以方便的分析风电并网给电力系统安全带来的影响。同时,该模型可以应用在日前电力系统的优化调度中,能有效的提升电力系统运行经济性和安全性。
Description
技术领域
本发明涉及一种计及风电不确定性的建模方法,尤其是一种计及风电不确定性的随机序列优化模型的建模方法。
背景技术
风能是当前世界上最具大规模商业化开发潜力的可再生能源。大规模开发利用风能发电,已成为世界各国解决能源问题和环境问题,改善能源结构,保证国民经济可持续发展的有效措施。吉林省可再生能源储量丰富,具备建设国家级清洁能源基地的条件。其中风电可装机容量5400万千瓦,是国家确定的9个千万千瓦风电风电基地之一。
我国“三北”地区风能资源丰富,但能源结构不合理,例如吉林省的热电矛盾异常突出。全省热电联产机组装机容量达到1314.84万千瓦,占燃煤火电装机的74%,比全省统调最小负荷多800万千瓦,比统调最大负荷多400万千瓦。进入供暖期后,电网调峰难度极大,即使热电联产机组全部按照最小方式运行,在夜晚低负荷期,风电仍被迫大量参与调峰。在春节等极端低负荷期,风电机组和纯凝机组必须全停,仍有300多万千瓦电力无法消纳,必须采取紧急措施并得到东北电网联络线支持才能保证电网安全。受此影响,吉林省弃风问题非常严重,近几年弃风率连续位列全国前列,供热中期风电大面积弃风,夜间风电全停成为常态,累计弃风率达到31.1%。
发明内容
本发明要解决的技术问题是由于风电的不确定性导致风电难以并网,造成大面积弃风。
为了解决上述技术问题,本发明提供了一种计及风电不确定性的随机序列优化模型的建模方法,其特征在于,包括如下步骤:
步骤1,利用概率性序列理论描述风电出力情况;
步骤2,建立基于概率性序列的线路潮流的约束模型;
步骤3,建立基于概率性序列的系统旋转储备约束模型。
通过建立线路潮流的约束模型以及系统旋转储备约束模型描述风电的不确定性以及风电并网对电力系统的影响,根据风电预测信息,实现日前热电联合经济调度,促进风电的消纳。
作为本发明的进一步限定方案,步骤1中利用概率性序列理论描述风电出力情况具体是:将风电出力离散化并赋予与之对应的概率,从而描述整个风电场的出力状况,同时定义并求解出风电场出力的期望值。
作为本发明的进一步限定方案,步骤2中建立基于概率性序列的线路潮流的约束模型具体是:计算出不同风电出力情况下的潮流分布,并赋予出现不同潮流分布时的对应概率,将当前运行方式下满足潮流约束的解所对应的概率相加,如果累加值大于设定的置信水平,则当前运行方式满足电力系统的安全稳定运行要求。
作为本发明的进一步限定方案,步骤3中建立基于概率性序列的系统旋转储备约束模型具体是:将当前运行方式下满足系统旋转储备要求的解所对应的概率相加,如果累加值大于设定的置信水平,则当前运行方式满足电力系统的安全稳定运行。
本发明的有益效果在于:(1)采用概率性序列理论,对风电的不确定性加以描述,方便计算机加以分析计算,与传统的蒙特卡洛法相比,在保证相同精度的前提下,大大节省了计算时间,不仅可以应用在日前优化调度中,也可以应用在调度时间间隔大于等于15分钟的短期经济调度中;(2)可以根据电网实际运行要求,对不同的线路设定不同的置信水平,不再统一的对所有线路设定相同的要求,可以单独分析风电并网对某一条线路的影响,从而帮助系统规划人员设计合理的线路容量;(3)可以直观的分析旋转储备容量对风电消纳的影响,为系统的安全稳定运行提供保障,减少基础建设成本,提高系统对风电的消纳能力,具有一定的经济和社会效益。
附图说明
图1为本发明的方法流程图;
图2为本发明的风电的离散化概率分布图;
图3为本发明的IEEE-30节点拓扑结构图;
图4为本发明的旋转储备约束与热负荷的关系柱状图。
具体实施方式
如图1所示,本发明公开的计及风电不确定性的随机序列优化模型的建模方法,包括如下步骤:
步骤1,利用概率性序列理论描述风电出力情况;
步骤2,建立基于概率性序列的线路潮流的约束模型;
步骤3,建立基于概率性序列的系统旋转储备约束模型。
其中,步骤1中利用概率性序列理论描述风电出力情况具体是:将风电出力离散化并赋予与之对应的概率,从而描述整个风电场的出力状况,同时定义并求解出风电场出力的期望值;具体步骤为:
利用风电的历史数据计算得到一段时间内的风电平均相对误差ε为:
由于风电预测误差服从正态分布,于是推导出风电平均相对误差ε与标准差σ的关系为:
式中,σ为风电预测的标准差,再将风电预测误差作为期望值,得到风电出力的概率密度函数为:
由于风电预测误差服从正态分布,所以风电出力在负半轴有一定的分布,然而负值是没有意义的,所以近似认为风电出力为零的概率是正态分布中小于等于0的部分对应的概率之和,将分布在负半轴上的概率归并至原点上,风电的离散化概率分布为:
风电的离散化概率分布图如附图2所示,设风电功率的离散化步长为ΔP,Nwp(t)是t时刻风电概率性序列的长度,由于在正态分布中,其变量离中心位置的距离超过三倍标准差的概率不到千分之三,所以Nwp(t)的取值可以选择覆盖99.7%的出力情况。
步骤2中建立基于概率性序列的线路潮流的约束模型具体是:计算出不同风电出力情况下的潮流分布,并赋予出现不同潮流分布时的对应概率,将当前运行方式下满足潮流约束的解所对应的概率相加,如果累加值大于设定的置信水平,则当前运行方式满足电力系统的安全稳定运行要求;具体步骤为:
首先建立供热平衡及有功平衡约束:
式中,COP表示电热泵的性能系数,PLoadt和HLoadt分别代表电负荷和热负荷,CG代表常规机组的集合,CHPG代表热电联产机组的集合,pi,t代表常规机组的电出力,代表热电联产机组的电出力,hj,t代表热电联产机组的热出力,代表热泵的热出力,代表热泵消耗的电功率,iwp代表风电出力,WPt(iwp)代表当风电出力为iwp时对应的概率,Ewp(t)代表风电出力的期望值,线路潮流约束如下:
Pline(t)≥βline (9)
式中,b1,b2,…,bN∈Bus,Bus是电网所有母线的集合,Pline(t)是用来计算满足潮流约束的概率之和的矩阵,是各线路所能承受的最大有功矩阵,βline是保障各线路运行安全的置信水平构成的矩阵,当某一条线路的有功i小于时,把他们对应的概率累加,得到该线路在这种运行方式下的安全概率,如果此概率大于事先设定好的置信水平βline,则该运行方式满足系统要求,当两个条母线不相连时,则对应元素的赋值如下:
这样的幅值可以保证不相连的母线满足下式:
由于式(10)-(12)三个矩阵都是稀疏矩阵,所以可以采用压缩行储存或者是压缩列储存,以节省存储空间,加快计算速度。
步骤3中建立基于概率性序列的系统旋转储备约束模型具体是:将当前运行方式下满足系统旋转储备要求的解所对应的概率相加,如果累加值大于设定的置信水平,则当前运行方式满足电力系统的安全稳定运行;具体步骤为:
式(15)-(20)中,pi,t代表常规机组i在t时刻的出力,和表示常规机组i的最大和最小出力,ΔT是调度时间间隔,rupi是机组i的向上爬坡速率,rdowni是机组i的向下爬坡速率,iwp代表风电场的实际出力,Ewp(t)代表风电场在t时刻的实际出力,式(17)和(18)是用来计算每台机组可提供的向上旋转储备容量和向下旋转储备容量,式(19)和(20)是用来挑选出满足系统旋转储备要求的储备容量,把它们的值记录在USR和DSR中,Rut(ru)和Rdt(rd)是t时刻的向上和向下旋转储备约束的概率性序列,把他们中满足系统要求的值所对应的概率累加,如果大于置信水平βRU和βRD,则此种运行方式满足系统的安全稳定要求。
附图3为IEEE-30节点的系统拓扑结构图,假设装机容量为40MW和90MW的风电场分别接在22号和27号母线上。热电联产机组接在1号和2号母线上。假设εwp等于该风电场10%的装机容量,置信水平βRU=βRD=95%,minβline≥85%。设置调度时间间隔为15分钟。
表1展示了随着系统增加对风电的消纳,风电对系统稳定性的影响。EWP表示该时刻风电的期望值,PUSR表示满足向上旋转储备约束的概率,PDSR表示满足向下旋转储备约束的概率。经潮流计算可知,线路L25-27承担着风电功率外送的任务,所以主要检查该线路上是否会出现过载的情况。表中L25-27表示该线路不过载的概率。随着风电期望值的增加,系统的向下旋转储备容量愈加不足,线路L25-27安全性降低,这也符合我们对风电消纳通常的认知。
表1实际风电功率对系统稳定性的影响
表2具体优化调度结果
表3与蒙特卡洛算法的比较
表2展示了具体的各个时刻发电机的出力和总燃煤成本。表3将本研究提出的算法与蒙特卡洛算法(MC)进行对比,展示了该算法和模型的有效性。在相同的计算精度条件下,本文提出的方法可以大幅减小计算时间。如表1示,线路潮流约束逐渐成为限制风电消纳的主要原因。所以,假设接下来给线路L25-27扩容,扩容后的优化结果如表4所示。CWP表示实际消纳风电的功率。结果显示,扩容后该线路不再是阻碍风电消纳的障碍,但是系统的旋转储备容量逐渐降低,系统的可靠性也随之降低。表5给出了旋转储备约束无法被满足的原因。
表4线路扩容后的优化结果
表5供热机组的优化结果
由表5可知,当消纳13MW风电时,供热1号机组电出力155.1MW,其电出力下限是105.7MW。然而,当消纳65MW风电时,供热1号机组电出力110MW,其电出力下限是109MW,已经没有办法进一步下调出力,所以导致了旋转储备容量的不足。可以看出,此时热电耦合约束限制了风电的消纳。以上展示的优化结果都是假设热负荷不变的情况下,随着实际风电功率的变化,各个参数的变化情况。附图4展示热负荷增加对风电消纳的影响。
当风电消纳量不变时,随着热负荷的增加,旋转储备约束越来越难以得到满足。系统消纳风电越多,满足向下旋转储备容量约束的概率下降得越快。例如当消纳52MW的风电时,满旋转储备约束的概率只有46%;但当风电消纳量较小时,例如仅消纳26MW的风电时,即使热负荷达到240MW也可以有较高的旋转储备容量。这说明风电的消纳给系统带来了更多的不确定性因素。表6展示了各个供热机组的热负荷的分配情况。H1是1号供热机组的热出力,H2是2号供热机组的热出力,由于1号机组的容量较大,且具有较好的经济性,所以通常1号机组会承担更多的热负荷。在风电可以被完全消纳的情况下,让2号机组承担更多的热负荷可以提高系统的向下旋转储备容量,但是会增加系统的煤耗成本。
表6供热机组的热出力情况(风电完全消纳)
表7供热机组的热出力情况(风电不完全消纳)
以上分析结果都是在假设风电可以被全部消纳的前提下得出的。但是,当风电消纳不完全时,可以适当的多分配热负荷给2号机组,提高系统的旋转储备容量,以促进风电的进一步消纳。利用多消纳的风电解决2号供热机组经济性较差的问题。表7展示了具体的优化结果,当风能资源较为充足的情况下,适当的分配一定的热负荷给经济性较差的机组,不仅可以提高系统的旋转储备容量,还可以促进风电的消纳,减少化石能源的燃烧。
Claims (1)
1.一种计及风电不确定性的随机序列优化模型的建模方法,其特征在于,包括如下步骤:
步骤1,利用概率性序列理论描述风电出力情况;
步骤2,建立基于概率性序列的线路潮流的约束模型;
步骤3,建立基于概率性序列的系统旋转储备约束模型;
步骤1中利用概率性序列理论描述风电出力情况具体是:将风电出力离散化并赋予与之对应的概率,从而描述整个风电场的出力状况,同时定义并求解出风电场出力的期望值;
将风电出力离散化并赋予与之对应的概率的具体步骤为:
利用风电的历史数据计算得到一段时间内的风电平均相对误差ε为:
由于风电预测误差服从正态分布,于是推导出风电平均相对误差ε与标准差σ的关系为:
式中,σ为风电预测的标准差,再将风电预测误差作为期望值,得到风电出力的概率密度函数为:
由于风电预测误差服从正态分布,所以风电出力在负半轴有一定的分布,风电出力为零的概率是正态分布中小于等于0的部分对应的概率之和,将分布在负半轴上的概率归并至原点上,则风电的离散化概率分布为:
式中,ΔP为风电功率的离散化步长,Nwp(t)是t时刻风电概率性序列的长度;
步骤2中建立基于概率性序列的线路潮流的约束模型具体是:计算出不同风电出力情况下的潮流分布,并赋予出现不同潮流分布时的对应概率,将当前运行方式下满足潮流约束的解所对应的概率相加,如果累加值大于设定的置信水平,则当前运行方式满足电力系统的安全稳定运行要求;
建立线路潮流的约束模型的具体步骤为:
建立供热平衡及有功平衡约束为:
式中,COP表示电热泵的性能系数,PLoadt和HLoadt分别代表电负荷和热负荷,CG代表常规机组的集合,CHPG代表热电联产机组的集合,pi,t代表常规机组的电出力,代表热电联产机组的电出力,hj,t代表热电联产机组的热出力,代表热泵的热出力,代表热泵消耗的电功率,iwp代表风电出力,WPt(iwp)代表当风电出力为iwp时对应的概率,Ewp(t)代表风电出力的期望值,则线路潮流约束如下:
Pline(t)≥βline (9)
式中,b1,b2,L,bN∈Bus,Bus是电网所有母线的集合,Pline(t)是用来计算满足潮流约束的概率之和的矩阵,是各线路所能承受的最大有功矩阵,βline是保障各线路运行安全的置信水平构成的矩阵,当某一条线路的有功i小于时,把他们对应的概率累加,得到该线路在这种运行方式下的安全概率,如果此概率大于事先设定好的置信水平βline,则该运行方式满足系统要求;
步骤3中建立基于概率性序列的系统旋转储备约束模型具体是:将当前运行方式下满足系统旋转储备要求的解所对应的概率相加,如果累加值大于设定的置信水平,则当前运行方式满足电力系统的安全稳定运行;
建立系统旋转储备约束模型的具体计算模型为:
式(15)-(20)中,pi,t代表常规机组i在t时刻的出力,和表示常规机组i的最大和最小出力,ΔT是调度时间间隔,rupi是机组i的向上爬坡速率,rdowni是机组i的向下爬坡速率,iwp代表风电场的实际出力,Ewp(t)代表风电场在t时刻的实际出力,式(17)和(18)是用来计算每台机组可提供的向上旋转储备容量和向下旋转储备容量,式(19)和(20)是用来挑选出满足系统旋转储备要求的储备容量,把他们的值记录在USR和DSR中,Rut(ru)和Rdt(rd)是t时刻的向上和向下旋转储备约束的概率性序列,把他们中满足系统要求的值所对应的概率累加,如果大于置信水平βRU和βRD,则此种运行方式满足系统的安全稳定要求。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710229367.4A CN106849176B (zh) | 2017-04-10 | 2017-04-10 | 一种计及风电不确定性的随机序列优化模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710229367.4A CN106849176B (zh) | 2017-04-10 | 2017-04-10 | 一种计及风电不确定性的随机序列优化模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106849176A CN106849176A (zh) | 2017-06-13 |
CN106849176B true CN106849176B (zh) | 2023-02-17 |
Family
ID=59148052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710229367.4A Active CN106849176B (zh) | 2017-04-10 | 2017-04-10 | 一种计及风电不确定性的随机序列优化模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106849176B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107508327B (zh) * | 2017-10-20 | 2020-07-24 | 云南电网有限责任公司 | 基于模型预测控制的配电网自治区域并网功率控制方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102880989B (zh) * | 2012-09-05 | 2016-04-20 | 中国电力科学研究院 | 一种风电出力时间序列建模方法 |
CN103997039B (zh) * | 2014-05-30 | 2017-01-18 | 国网甘肃省电力公司 | 基于概率区间预测的计及风电接纳的旋转备用区间预测方法 |
CN106410859B (zh) * | 2015-07-28 | 2019-04-05 | 中国电力科学研究院 | 一种基于随机生产模拟的新能源并网接纳能力计算方法 |
-
2017
- 2017-04-10 CN CN201710229367.4A patent/CN106849176B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106849176A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109449973B (zh) | 一种含光伏和储能的轨道交通供电系统能量优化方法 | |
CN104361416B (zh) | 一种考虑大规模电动汽车接入的电网双层优化调度方法 | |
CN103248048B (zh) | 一种弃风电量预估方法 | |
CN106712033A (zh) | 一种热电厂消纳弃风的方法 | |
CN114069688A (zh) | 一种基于时序生产模拟的多电源容量布局规划方法 | |
Chen et al. | Scheduling strategy of hybrid wind-photovoltaic-hydro power generation system | |
CN115081700A (zh) | 基于综合储能技术的数据中心多能协同优化方法及系统 | |
CN103441494B (zh) | 基于非并网多能源协同供电的智能电网系统 | |
CN113078684B (zh) | 一种基于双层优化的区域能源社区规划方法 | |
CN106849176B (zh) | 一种计及风电不确定性的随机序列优化模型的建模方法 | |
CN112653137A (zh) | 计及碳交易的光热电站与风电系统、低碳调度方法及系统 | |
CN112736969A (zh) | 基于新能源经济调度的分布式光伏数据处理方法及系统 | |
CN105514985B (zh) | 一种电网聚合模型的构建方法 | |
CN112767189B (zh) | 一种计及风电并网下电网灵活性评估方法 | |
Julianto et al. | Confronting the Duck Curve Problem Using Dynamic Economic Emission Dispatch with CAES. | |
Hao et al. | PSO applied to optimal operation of a micro-grid with wind power | |
Zhou et al. | Research on Pumped Storage Capacity Allocation of Cascade Hydro-Wind-Solar-Pumped Storage Hybrid System Considering Economy and Operational Stability | |
Chen et al. | Study on Wind Power Accommodation Method Based on Controllable Electric Heating | |
Han et al. | Research on the Optimal Scheduling Strategy of Energy Storage Plants for Peak-shaving and Valley-filling | |
Deng et al. | Analysis of Renewable Energy Accommodation Capability of Shanxi Power Grid Based on Operation Simulation Method | |
Du et al. | Configuration of thermal storage tank of microgrid clusters considering thermal interaction | |
Pan et al. | A multi-source coordinated spinning reserve capacity optimization considering wind and photovoltaic power uncertainty | |
Xu et al. | Optimal Scheduling Method for Cogeneration System Based on Heat Storage Device | |
Zhang | Study on the Effects of Different Measures in Promoting Renewable Energy Consumption | |
Xiao et al. | Complementary coordinated dispatch of multi-energy systems considering power flexibility margin |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |