CN109347152A - 考虑多类型电源参与调峰的随机生产模拟方法及应用 - Google Patents
考虑多类型电源参与调峰的随机生产模拟方法及应用 Download PDFInfo
- Publication number
- CN109347152A CN109347152A CN201811457102.0A CN201811457102A CN109347152A CN 109347152 A CN109347152 A CN 109347152A CN 201811457102 A CN201811457102 A CN 201811457102A CN 109347152 A CN109347152 A CN 109347152A
- Authority
- CN
- China
- Prior art keywords
- load
- power
- unit
- peak
- load curve
- 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.)
- Granted
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 49
- 240000002853 Nelumbo nucifera Species 0.000 claims abstract description 58
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims abstract description 58
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims abstract description 58
- 238000004088 simulation Methods 0.000 claims abstract description 42
- 230000005611 electricity Effects 0.000 claims abstract description 40
- 230000002045 lasting effect Effects 0.000 claims abstract description 29
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 23
- 239000000446 fuel Substances 0.000 claims description 10
- 238000010248 power generation Methods 0.000 claims description 8
- 230000001105 regulatory effect Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000013461 design Methods 0.000 claims description 4
- 238000013439 planning Methods 0.000 claims description 4
- 230000007812 deficiency Effects 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 7
- 238000010586 diagram Methods 0.000 description 10
- 239000003245 coal Substances 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000005619 thermoelectricity Effects 0.000 description 5
- 238000003860 storage Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 description 2
- 229910052796 boron Inorganic materials 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000033001 locomotion Effects 0.000 description 2
- 230000005622 photoelectricity Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009257 reactivity Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000029087 digestion Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003758 nuclear fuel Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 230000007096 poisonous effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- 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
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Power Engineering (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种考虑多类型电源参与调峰的随机生产模拟方法及应用，所述随机生产模拟方法包括以下步骤：1)获取全年负荷曲线；2)从所述全年负荷曲线中剔除风电、光伏和核电全年出力，获得持续负荷曲线，并根据全年负荷水平将所述持续负荷曲线划分为基荷、腰荷、峰荷三个阶段；3)根据所述持续负荷曲线建立等效电量函数；4)根据基荷、腰荷和峰荷的等效电量，实现火电机组和水电机组在不同阶段的出力，获得电力系统随机生产模拟数据。与现有技术相比，本发明具有提高区域电网的调峰能力等优点。
Description
技术领域
本发明涉及电力系统规划与运行控制技术领域，尤其是涉及一种考虑多类型电源参与调峰的随机生产模拟方法及应用。
背景技术
电力系统随机生产模拟是电力系统规划与运行的基本分析工具之一，它通过优化发电机组的生产运行，考虑机组故障和负荷的不确定性，计算各电厂的发电量、系统生产成本及系统的可靠性指标。传统地，核电机组不参加电网调峰，这并不是由于机组设计方面的原因，而是因为核燃料成本较高，核电所占比重较小。随着核电、风电和光伏发电的发展建设，仅靠火电和抽水蓄能机组进行调峰将难以满足电网负荷变化的需要，核电机组有必要参与电网调峰。另一方面，在水电资源丰富的地区，水电调峰能力由于外部天气条件和水电站类型的不同而差别较大，水电机组可以根据调节能力由弱到强，依次细分为年调型水电机组、季调型水电机组和日调型水电机组，有必要研究并充分发挥不同类型机组在不同天气条件下的调峰能力。因此，研究包括核电、不同类型水电机组等多类型电源参与调峰的电力系统随机生产模拟方法对合理安排发电计划、降低成产成本、提高可靠性具有重要意义。
国内对随机生产模拟问题进行了广泛研究。陈亚博、盛戈皞、黎建等在《电力系统及其自动化学报》(2015,27(05):1-6)上发表的《含光伏和风电的电力系统随机生产模拟》基于随机生产模拟方法以及其中多状态机组的处理方法,建立了光伏的多状态计算模型。同时，结合该模型和风电等效多状态模型分析处理了光伏和风电两种间歇式能源同时并入电力系统后的随机生产模拟情况。张宏宇、方鑫、李碧辉等在《中国电力》(2012,45(06):73-767)上发表的《含大规模风光电源电力系统随机生产模拟》根据风电与光伏电源的自身电源特性及其与电网的相互关系，提出风电与光伏电源参与随机生产模拟的原则。根据所提原则对应提出风功率损失期望及弃风比率2个评价指标。冯长有、梁志峰在《电网技术》(2013,37(02):493-499)上发表的《考虑潮流断面约束的电力系统随机生产模拟》结合电网实际运行情况,根据有无电量交易和断面构成情况，划分了潮流断面类型并构建了相应数学描述模型。基于此，给出了等效电量函数法中不同类型断面的处理方法，设计了相应的随机生产模拟计算流程。姚力、黄镔、王秀丽等在《电网技术》(2015,39(05):1219-1225)上发表的《考虑风火联合外送的互联系统随机生产模拟》文章采用风电多状态模型，引入机组剩余容量的概念，提出了联合外送概率模型以模拟送端电网对受端电网的电力支援过程。结合各类型电源特点，设计了考虑风火联合外送的互联系统随机生产模拟算法框架。以上文献中，都着重考察风、光电波动性、间歇性、随机性对电力系统生产模拟的影响，而没有考虑到核电机组、多类型水电机组参与调峰的因素。专利中，范宏、左路浩、马莲等发明人申请的发明专利《能效电厂的随机生产模拟计算方法》针对能效电厂需求侧管理方面的创新与对风电光电机组出力的精确模拟，探究能效电厂与大型光伏风电站出力的对电力系统运行的影响。施涛、韩华玲、朱凌志等发明人申请的发明专利《一种含大规模光伏电站的电力系统随机生产模拟方法》确定不含光伏电站的系统机组等值有效容量各阶累积量；依次获得对应时段的光伏发电的时序有效容量概率分布Beta模型及参数；依次获得对应时段的含光伏电站的系统机组等值有效容量分布概率模型及参数及对应负荷水平下的可靠性指标。吴红斌、白雪、丁明等发明人申请的发明专利《一种用于风光储联合发电系统的随机生产模拟方法》采用K-Means算法进行气象模式划分将风力发电模块、光伏发电模块、储能装置及负荷的序列化，基于序列运算理论进行随机生产模拟，更好地描述了风光互补特性，并进行风光储联合发电系统的随机生产模拟，从而能为风光储联合发电系统的运行与调度提供参考。以上专利的随机生产模拟进一步考虑了抽水蓄能电站和负荷变化的影响，但是仍然没有考虑核电、多类型水电机组参与调峰的影响。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种考虑多类型电源参与调峰的随机生产模拟方法及应用。
本发明的目的可以通过以下技术方案来实现：
一种考虑多类型电源参与调峰的随机生产模拟方法，包括以下步骤：
1)获取全年负荷曲线；
2)从所述全年负荷曲线中剔除风电、光伏和核电全年出力，获得持续负荷曲线，并根据全年负荷水平将所述持续负荷曲线划分为基荷、腰荷、峰荷三个阶段；
3)根据所述持续负荷曲线建立等效电量函数；
4)根据基荷、腰荷和峰荷的等效电量，实现火电机组和水电机组在不同阶段的出力，获得电力系统随机生产模拟数据。
进一步地，所述持续负荷曲线的获取步骤包括：
201)获取全年风电和光伏出力预测曲线，将风电与光伏出力从所述全年负荷曲线中剔除；
202)考虑核电参与调峰，设计核电在峰谷差大于设定值的典型日以“12-3-6-3”调峰模式出力，其余时段以额定功率带基荷，获得核电全年出力时序曲线；
203)基于步骤201)，再将所述核电全年出力时序曲线从全年负荷曲线中剔除，获得持续负荷曲线。
进一步地，所述等效电量函数的建立过程包括：
301)在持续负荷曲线上，将横轴x按Δx分段，将所述持续负荷曲线转换为离散的初始电量函数：
式中，k＝＜x/Δx＞+1，＜x/Δx＞表示不大于x/Δx的整数，E(k)等于从x到x+Δx这一段负荷曲线对应的负荷电量，F(x)表示持续负荷曲线；
302)考虑发电机组随机停运的影响，以初始电量函数为基础进行修正，获得最终的等效电量函数：
E(i)(k)＝piE(i-1)(k)+qiE(i-1)(k-mi)
式中，E(i)(k)为运行第i台发电机组运行后的等效电量函数，qi为强迫停运率，pi＝1-qi，mi＝Ci/Δx，Ci为第i台发电机组的装机容量。
进一步地，所述实现火电机组和水电机组在不同阶段的出力具体包括以下步骤：
401)根据经济性对火电机组排序，优先运行所有火电机组按最小出力段带基荷；
402)根据调节能力对水电机组排序，优先运行所有日调型机组带基荷；
403)判断基荷是否分配完成，若是，则执行步骤404)，若否，则依次运行季调型机组、年调型机组带基荷；
404)优先运行所有火电机组带腰荷，然后依次运行日调型、季调型机组带腰荷；
405)判断腰荷是否分配完成，若是，则执行步骤406)，若否，则运行年调型机组带腰荷；
406)依次运行年调型、季调型、日调型机组带峰荷；
407)判断峰荷是否分配完成，若是，则结束，若否，则运行火电机组带峰荷。
进一步地，所述水电机组在生产模拟中满足以下条件：
式中，PHL为水电机组担任的最大负荷功率，CH为水电机组容量，EL为负荷电量，EA为水电机组给定电量。
进一步地，单台所述水电机组带负荷的位置通过以下方式确定：
在所述持续负荷曲线下构建一台水电机组的特征矩形，该特征矩形的宽度为水电机组容量，高度为水电机组在模拟周期内的利用小时数，特征矩形底部与特征矩形的横轴重合，使该特征矩形在持续负荷曲线下移动，当该特征矩形的面积等于持续负荷曲线在相应区间内的面积时，将所述相应区间作为该台水电机组带负荷的位置。
进一步地，该方法还包括步骤：
5)根据电力系统随机生产模拟数据计算发电量、生产成本和可靠性指标。
进一步地，所述可靠性指标包括系统电量不足期望值、电力不足概率和电力不足期望值。
进一步地，所述生产成本包括燃料费用、火电机组启动费用和停机费用。
本发明还提供一种电力系统规划与运行方法，该方法包括：
采用如所述的考虑多类型电源参与调峰的随机生产模拟方法获得随机生产模拟数据。
与现有技术相比，本发明具有以如下有益效果：
1)本发明在随机生产模拟中考虑了核电参与调峰的影响，可以平抑负荷曲线的波动，降低峰谷差，进一步提高区域电网的调峰能力。
2)本发明在随机生产模拟中分为年调型水电机组、季调型水电机组、日调型水电机组考虑了水电参与调峰的影响，根据水电机组的调节能力的不同合理安排各类型水电机组参与调峰，可以充分消纳水电，提高发电厂年发电计划的经济性。
附图说明
图1为本发明的流程示意图；
图2为持续负荷曲线示意图；
图3为单个水电机组带峰荷情况示意图；
图4为单个水电机组带负荷位置的确定示意图；
图5为H电网的总生产成本示意图；
图6为H电网的LOLP指标示意图；
图7为H电网的EENS指标示意图；
图8为H电网的核电运行成本示意图；
图9为H电网的火电运行成本示意图；
图10为H电网的水电运行成本示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
本发明提供一种考虑多类型电源参与调峰的随机生产模拟方法，该方法在传统的电力系统随机生产模拟基础上进一步考虑核电、日调型、季调型、年调型水电机组参与调峰的情况，主要技术要点包括：
1)获取全年负荷曲线、全年风电、光伏出力预测曲线、火电机组煤耗系数、水电机组给定电量以及所有类型机组强迫停运率、发电成本等基础数据；
2)设计考虑核电与日调型、季调型、年调型水电机组参与调峰的随机生产模拟流程，形成持续负荷曲线以及负荷频率曲线进而用等效电量函数法对各类型机组进行随机生产模拟计算；
3)求得多类型发电机组参与调峰的情形下，各类型机组年发电量、生产成本、电量不足期望值、电力不足概率等可靠性指标，分析多类型发电机组参与调峰对随机生产模拟结果的影响。
上述随机生产模拟方法中，对核电、火电机组及水电机组的运行情况进行了设计，具体原理如下。
在核电站实际运行中，压水堆一般采用控制棒和硼浓度作为控制量来实现对反应堆负荷跟踪的控制。其中控制棒控制快速的反应性变化，如机端负荷变化；硼浓度控制慢速的反应变化，主要用于补偿燃耗、堆芯毒物的反应性变化。核电机组在可调节区域运行时，功率调节速率一般约为±0.2％Pn/min-±0.3％Pn/min，最高不能超过±5％Pn/min的极限，且能够进行功率变化范围在±10％Pn以内的阶跃变化，但机组不能进行连续的功率阶跃变化和线性变化，考虑到调峰的深度和速度限制，压水堆核电机组可采用“12-3-6-3”的日负荷跟踪模式，即以额定功率运行12h，在3h内完成降功率，以30％Pn-70％Pn运行6h(压水堆核电机组的负荷跟踪阈值为30％Pn-100％Pn。留有70％Pn的裕度)，再在3h内完成升功率继续以12h额定功率运行。另外，考虑到核电机组进行频繁的功率调节，会加剧反应堆压力容器的辐照脆化、冷却剂循环系统及蒸汽供应系统的某些关键金属部件的金属疲劳，可能降低机组设备的寿命，核电机组不应全年均采用“12-3-6-3”的日负荷跟踪模式参与调峰。
火电机组根据其经济性进行安排。
水电机组按调峰能力由强到弱依次分为：年调型水电机组、季调型水电机组、日调型水电机组。与火电机组相比，水电机组具有以下特点：
1)水电机组不消耗燃料，其发电费用几乎与发电量无关；
2)水电机组的发电量和发电出力是由水文条件所决定的，在随机生产模拟中其发电量应该看作给定量水电机组属于给定电量的发电机组，在进行生产模拟时其发电量作为给定条件。
3)水电机组开停灵活，增减负荷迅速，适合担任负荷曲线的峰荷部分。
为了充分发挥水电机组的效益，在安排水电运行时应遵循以下两条原则：
1)充分利用水能资源，尽量避免弃水。
2)利用水电机组来带峰荷部分，这样可以取代煤耗量大的火电机组，从而可以取得较好的经济效益和环境效益，因此当系统中有水电机组时应尽可能用水电就来承担峰荷，以到达降低煤耗的效果。
如图1所示，上述随机生产模拟方法包括数据输入与预处理、随机生产模拟计算、可靠性指标与系统生产成本计算三个部分，步骤具体描述如下：
步骤S101，获取全年负荷曲线和全年风电光伏出力预测曲线，从所述全年负荷曲线中剔除风电、光伏出力。
步骤S102，安排核电机组跟踪负荷曲线。
安排核电参与调峰的方法是：在全年时序负荷曲线中峰谷差较大的典型日，优先安排核电以上述“12-3-6-3”模式参与调峰，全年其余时段核电以额定功率带基荷。
步骤S103，根据上述得出核电全年出力时序曲线，将此曲线从负荷曲线中减去，得到净负荷曲线，然后形成持续负荷曲线如图2所示，然后根据全年负荷水平，将持续负荷曲线分为基荷、腰荷、峰荷三段。
t＝F(x)为研究周期T内系统的负荷持续曲线，其横坐标表示系统负荷，纵坐标表示负荷的持续时间，PLmax表示系统最大负荷，(x,t)表示系统负荷大于或等于x的持续时间为t。
步骤S104，选取发电机组容量最大公因子作为步长，形成初始等效电量函数。
具体方法为：在持续负荷曲线上，将横轴x按Δx分段，定义一个离散的电量函数：
式中：k＝＜x/Δx＞+1；＜x/Δx＞表示不大于x/Δx的整数；E(k)等于从x到x+Δx这一段负荷曲线对应的负荷电量。
为了考虑发电机组的停运因素，在等效电量函数中需要用发电机组随机停运的影响来修正电量函数。设原始持续负荷曲线为F(0)(x)，对应的电量函数为E(0)(k)，在安排完第i-1台发电机组运行以后得到等效持续负荷曲线F(i-1)(x)，对应的等效电量函数为E(i-1)(k)。设第i台发电机组的装机容量为Ci，强迫停运率为qi，则安排第i台发电机组运行后的等效电量函数为：
E(i)(k)＝piE(i-1)(k)+qiE(i-1)(k-mi)
式中：pi＝1-qi；mi＝Ci/Δx。
步骤S105，根据经济性对火电机组进行安排，优先安排所有火电机组按最小出力段带基荷。
步骤S106，在基荷段优先调峰能力最弱的水电机组出力，即优先安排日调型水电机组带基荷。
单台水电机组担任峰荷的情况如图3所示。图中曲线cg是由原负荷曲线向左平移相当于水电机组容量CH而来。阴影部分的面积应等于水电机组的给定量EA。在这种情况下其余机组应担任的负荷为oacgfh所围成的部分。在距a点CH(水电机组容量)的b点作垂线be，则图形acg的面积和图形bde的面积相等。也就是说其余机组所担任的负荷可以看成由Oafh和bde两部分组成。这相当于水电机组承担了图中abef部分的负荷。
因此将单个水电机组的生产模拟中的处理原则归结为：
在等效负荷曲线下寻找相当于水电机组容量CH的一段，其间面积恰好等于水电机组的给定电量EA。即，水电机组在生产模拟中应满足以下条件：
式中PHL为水电机组担任的最大负荷功率，CH为水电机组容量，EL为负荷电量，EA为水电机组给定电量。
图4表示在生产模拟中确定水电机组运行位置的过程。首先在等效负荷曲线下作出水电机组的特征矩形abb’a’，其底为CH，高为水电机组在模拟周期内的利用小时数TH。当把这个矩形向右移动，使负荷曲线在相应区间内的面积等于该矩形面积时，找到了水电机组运行位置。
上述水电机组的特征矩形向左移功的过程，实际上是顺次安排水电机组运行的过程，每安排一台水电机组特征矩形将向右移动一段与该火电机组容量相应的距离。由于这种移动是不连续的。很难保证上述条件恰好成立。因此在算法上，将水电机组容量CH和水电机组的给定电量EA作为安排水电机组时的上边界条件进行计算。对于多台水电机组的情况，设系统中有NH台水电机组。把它们的特征矩形按其高度(利用小时数)大小从左到右排列，形成水电机组特征矩形序列图。当把序列图从左向右移动，在等效持续负荷曲线的某一区间满足以下条件时：
前n个水电机组即可合并为一个等效水电机组，并带相应位置的负荷。剩余的NH-n台水电机组的矩形序列图应继续向右移动，并在边界条件的区间合并为另一等效水电机组，带该区间的负荷。
步骤S107，判断是否有剩余基荷，若有剩余基荷，则执行步骤S108，否则进入步骤S109，进行腰荷段负荷分配。
步骤S108，依次安排季调型机组、年调型机组带基荷。
步骤S109，优先安排所有火电机组带腰荷。
步骤S110，依次安排日调型、季调型水电机组带腰荷。
考虑到与火电机组相比，水电机组带峰荷具有更高经济效益，在腰荷段优先安排火电机组带腰荷，然后按照调节能力由弱到强的顺序依次安排日调型、季调型水电机组带腰荷。
步骤S111，判断是否有剩余腰荷，若有剩余腰荷则执行步骤S112，否则进入步骤S113进行峰荷段负荷分配。
步骤S112，安排年调型机组带腰荷。
步骤S113，在峰荷段，优先安排年调型水电机组带峰荷。
步骤S114，安排季调型水电机组带峰荷。
步骤S115，安排日调型水电机组带峰荷。
步骤S116，判断是否有剩余峰荷，若还有剩余峰荷则执行步骤S117，否则进入步骤S118。
步骤S117，安排火电机组带峰荷。
步骤S118，计算区域电网生产成本和可靠性指标。
(1)发电量与可靠性指标
第i台发电机组的发电量应根据等效电量函数E(i-1)(k)来进行计算：
式中：ki-1＝xi-1/Δx；ki＝(xi-1+Ci)/Δx。假设系统中共有n台发电机组，总容量为Cs，当全部发电机组卷积运算结束以后，等效电量函数为E(n)(k)。
1)系统电量不足期望值EENS：
2)电力不足概率LOLP：
3)电力不足期望值EDNS(Expected Demand Not Supplied)：
表示所研究的一段时间内系统电力不足的期望值：
EDNS＝EENS/T
上式中：EENS表示电量不足期望值，T表示研究周期。
(2)生产成本分析
系统总的生产成本包括燃料费用、火电机组启动费用、停机费用。启动费用包括机组启动的燃料费用、人工费用以及频繁启停造成的机组磨损及寿命折损。汽轮机启动费用还包括机组启动的运行和维护费用。在不同的启动条件下，启动费用略有不同，但可以近似看作是固定的。在某一时间段内，机组的停机次数和机组的启动次数是一样的，因此，停机费用一般也是固定的，可以综合到机组的启动费用里，即：
Ctotal＝Cfuel+Ctsu+Csd
式中：Cfuel为燃料费用；Ctsu为火电机组的汽轮机启动费用；Csd为停机费用。
式中：cfuel,i为第i台机组单位发电量的燃料费用；ctsu,i为第i台机组单位发电量的汽轮机启动费用；csd,i为第i台机组单位发电量的停机费用；为第i台机组的期望开机次数；为第i台发电机组的发电量，T为研究周期，pi等于1减去强迫停运率的值。
实施例
本实施例将考虑多类型电源参与调峰的电力系统随机生产模拟方法应用于某水平年某实际H区域电网，计算了核电机组、各类型水电机组参与调峰时的可靠性指标及生产成本指标。该年度H区域电网总装机容量372598MW，其中燃煤机组装机容量227751MW，占比61.3％，水电机组装机容量8062MW，占比2.16％；核电装机容量35724MW，占比9.5％，最高负荷284800MW。部分计算数据设置如下：煤价取为560元/t，气价为4000元/t，燃煤机组环境成本为101.246元/MWh，燃气机组环境成本为19.285元/MWh；水电平均运行成本为60元/MWh，日调型、季调型、年调型水电机组的基荷电量、腰荷电量、峰荷电量设置如表1；核电机组，以CPR1000核电机组为基准，平均运行成本为222.4元/MWh，日峰谷差率高于20％时核电机组参与调峰，调峰深度为60％。同时设置一组对比算例，其核电全年以额定功率带基荷，同时未考虑类型各水电机组调峰能力差异，计算可靠性指标及生产成本指标。
表1各类型水电机组的基荷、腰荷、峰荷段电量
考虑多类型机组参与调峰的某水平年H区域电网随机生产模拟结果如表2所示，H电网的可靠性指标及各项成本如图5-图10所示。
表2考虑多类型机组参与调峰的某水平年H区域电网随机生产模拟结果
表3为现有一些发电可靠性标准(数据来自清华大学学术专著：《电力系统可靠性分析》，郭永基著)，可靠性方面，对比表3与表2可以看出，以电力不足期望值EDNS作为指标，不考虑核电及各类型水电机组调峰性能差异时，H电网的EDNS为1.06，属于较高水平，考虑多类型机组调峰后的EDNS为0.22，可靠性水平得到显著提升。一方面是由于系统装机容量充足；另一方面，核电与各类型水电机组参与调峰后有效减轻系统调峰压力，平抑负荷曲线波动，进一步提高了系统的可靠性水平。
表3现有发电可靠性标准
经济性方面，根据表2可以看出，考虑多类型电源参与调峰后，系统总生产成本降低，经济性得到显著提升，其中火电运行成本降低81.5亿元，核电运行成本减低32.5亿元，水电成本增加2.3亿元。火电成本大幅降低是因为峰荷优先安排核电与水电机组承担，基荷与腰荷优先安排火电机组承担，减少了火电机组以较高经济代价承担的峰荷；核电成本降低是因为核电在峰谷差率较大的典型日按“12-3-6-3”模式参与调峰，减少了全年满功率运行时长，导致运行成本降低；水电成本增加是因为与传统随机生产模拟安排完火电机组后再安排水电机组带峰荷的方法相比，充分考虑各类型水电机组调节能力不同后，安排调节能力差的机组优先带基荷、腰荷，调节能力强的机组优先带峰荷，可以优先消纳更多水电，增强了系统水电消纳能力，因而水电成本上升。综合来看，由于水电成本远低于核电、火电运行成本，系统总生产成本大幅降低，经济效益显著。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。
Claims (10)
1.一种考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，包括以下步骤：
1)获取全年负荷曲线；
2)从所述全年负荷曲线中剔除风电、光伏和核电全年出力，获得持续负荷曲线，并根据全年负荷水平将所述持续负荷曲线划分为基荷、腰荷、峰荷三个阶段；
3)根据所述持续负荷曲线建立等效电量函数；
4)根据基荷、腰荷和峰荷的等效电量，实现火电机组和水电机组在不同阶段的出力，获得电力系统随机生产模拟数据。
2.根据权利要求1所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，所述持续负荷曲线的获取步骤包括：
201)获取全年风电和光伏出力预测曲线，将风电与光伏出力从所述全年负荷曲线中剔除；
202)考虑核电参与调峰，设计核电在峰谷差大于设定值的典型日以“12-3-6-3”调峰模式出力，其余时段以额定功率带基荷，获得核电全年出力时序曲线；
203)基于步骤201)，再将所述核电全年出力时序曲线从全年负荷曲线中剔除，获得持续负荷曲线。
3.根据权利要求1所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，所述等效电量函数的建立过程包括：
301)在持续负荷曲线上，将横轴x按Δx分段，将所述持续负荷曲线转换为离散的初始电量函数：
式中，k＝＜x/Δx＞+1，＜x/Δx＞表示不大于x/Δx的整数，E(k)等于从x到x+Δx这一段负荷曲线对应的负荷电量，F(x)表示持续负荷曲线；
302)考虑发电机组随机停运的影响，以初始电量函数为基础进行修正，获得最终的等效电量函数：
E(i)(k)＝piE(i-1)(k)+qiE(i-1)(k-mi)
式中，E(i)(k)为运行第i台发电机组运行后的等效电量函数，qi为强迫停运率，pi＝1-qi，mi＝Ci/Δx，Ci为第i台发电机组的装机容量。
4.根据权利要求1所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，所述实现火电机组和水电机组在不同阶段的出力具体包括以下步骤：
401)根据经济性对火电机组排序，优先运行所有火电机组按最小出力段带基荷；
402)根据调节能力对水电机组排序，优先运行所有日调型机组带基荷；
403)判断基荷是否分配完成，若是，则执行步骤404)，若否，则依次运行季调型机组、年调型机组带基荷；
404)优先运行所有火电机组带腰荷，然后依次运行日调型、季调型机组带腰荷；
405)判断腰荷是否分配完成，若是，则执行步骤406)，若否，则运行年调型机组带腰荷；
406)依次运行年调型、季调型、日调型机组带峰荷；
407)判断峰荷是否分配完成，若是，则结束，若否，则运行火电机组带峰荷。
5.根据权利要求1所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，所述水电机组在生产模拟中满足以下条件：
式中，PHL为水电机组担任的最大负荷功率，CH为水电机组容量，EL为负荷电量，EA为水电机组给定电量。
6.根据权利要求4所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，单台所述水电机组带负荷的位置通过以下方式确定：
在所述持续负荷曲线下构建一台水电机组的特征矩形，该特征矩形的宽度为水电机组容量，高度为水电机组在模拟周期内的利用小时数，特征矩形底部与特征矩形的横轴重合，使该特征矩形在持续负荷曲线下移动，当该特征矩形的面积等于持续负荷曲线在相应区间内的面积时，将所述相应区间作为该台水电机组带负荷的位置。
7.根据权利要求1所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，该方法还包括步骤：
5)根据电力系统随机生产模拟数据计算发电量、生产成本和可靠性指标。
8.根据权利要求7所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，所述可靠性指标包括系统电量不足期望值、电力不足概率和电力不足期望值。
9.根据权利要求7所述的考虑多类型电源参与调峰的随机生产模拟方法，其特征在于，所述生产成本包括燃料费用、火电机组启动费用和停机费用。
10.一种电力系统规划与运行方法，其特征在于，该方法包括：
采用如权利要求1所述的考虑多类型电源参与调峰的随机生产模拟方法获得随机生产模拟数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811457102.0A CN109347152B (zh) | 2018-11-30 | 2018-11-30 | 考虑多类型电源参与调峰的随机生产模拟方法及应用 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811457102.0A CN109347152B (zh) | 2018-11-30 | 2018-11-30 | 考虑多类型电源参与调峰的随机生产模拟方法及应用 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109347152A true CN109347152A (zh) | 2019-02-15 |
| CN109347152B CN109347152B (zh) | 2022-01-18 |
Family
ID=65319151
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811457102.0A Active CN109347152B (zh) | 2018-11-30 | 2018-11-30 | 考虑多类型电源参与调峰的随机生产模拟方法及应用 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109347152B (zh) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110854933A (zh) * | 2019-11-26 | 2020-02-28 | 三峡大学 | 利用灵活性资源的月度机组组合优化方法 |
| CN112398176A (zh) * | 2020-11-05 | 2021-02-23 | 贵州电网有限责任公司 | 一种考虑燃煤机组启停调峰的水火风互济系统日前优化调度方法 |
| CN113541204A (zh) * | 2021-08-10 | 2021-10-22 | 上海电力大学 | 一种燃煤发电机组全负荷调峰的自适应调节方法及系统 |
| CN114220560A (zh) * | 2021-12-13 | 2022-03-22 | 中国核动力研究设计院 | 一种压水堆核电厂硼控调峰方法及系统 |
| CN115400864A (zh) * | 2021-05-28 | 2022-11-29 | 北京京能电力股份有限公司 | 一种自动判断调峰火电机组磨煤机自动启停时机的方法 |
| US20220385064A1 (en) * | 2019-11-27 | 2022-12-01 | Robert F. Cruickshank, III | An optimized load shaping system, method & apparatus for optimizing production and consumption of energy |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015031581A1 (en) * | 2013-08-28 | 2015-03-05 | Robert Bosch Gmbh | System and method for energy asset sizing and optimal dispatch |
| CN104617591A (zh) * | 2015-01-28 | 2015-05-13 | 西安交通大学 | 基于多场景新能源发电模拟的日运行方式安排和调峰方法 |
| CN107104467A (zh) * | 2017-06-05 | 2017-08-29 | 武汉大学 | 一种考虑核电调峰及其安全约束的机组组合优化方法 |
| CN107644116A (zh) * | 2017-08-02 | 2018-01-30 | 广东电网有限责任公司肇庆供电局 | 一种适应于间歇性能源接入的随机生产模拟算法 |
| CN107910883A (zh) * | 2017-12-01 | 2018-04-13 | 中国南方电网有限责任公司电网技术研究中心 | 基于抽水蓄能电站修正时序负荷曲线的随机生产模拟方法 |
| CN108876026A (zh) * | 2018-06-08 | 2018-11-23 | 国网山东省电力公司青岛供电公司 | 兼顾特高压电网及区外来电的电源优化配置方法 |
-
2018
- 2018-11-30 CN CN201811457102.0A patent/CN109347152B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015031581A1 (en) * | 2013-08-28 | 2015-03-05 | Robert Bosch Gmbh | System and method for energy asset sizing and optimal dispatch |
| CN104617591A (zh) * | 2015-01-28 | 2015-05-13 | 西安交通大学 | 基于多场景新能源发电模拟的日运行方式安排和调峰方法 |
| CN107104467A (zh) * | 2017-06-05 | 2017-08-29 | 武汉大学 | 一种考虑核电调峰及其安全约束的机组组合优化方法 |
| CN107644116A (zh) * | 2017-08-02 | 2018-01-30 | 广东电网有限责任公司肇庆供电局 | 一种适应于间歇性能源接入的随机生产模拟算法 |
| CN107910883A (zh) * | 2017-12-01 | 2018-04-13 | 中国南方电网有限责任公司电网技术研究中心 | 基于抽水蓄能电站修正时序负荷曲线的随机生产模拟方法 |
| CN108876026A (zh) * | 2018-06-08 | 2018-11-23 | 国网山东省电力公司青岛供电公司 | 兼顾特高压电网及区外来电的电源优化配置方法 |
Non-Patent Citations (1)
| Title |
|---|
| 江伟民: "抽水蓄能电站经济效应计算方法研究", 《CNKI优秀硕士学位论文全文库》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110854933A (zh) * | 2019-11-26 | 2020-02-28 | 三峡大学 | 利用灵活性资源的月度机组组合优化方法 |
| CN110854933B (zh) * | 2019-11-26 | 2023-05-02 | 三峡大学 | 利用灵活性资源的月度机组组合优化方法 |
| US20220385064A1 (en) * | 2019-11-27 | 2022-12-01 | Robert F. Cruickshank, III | An optimized load shaping system, method & apparatus for optimizing production and consumption of energy |
| US11735919B2 (en) * | 2019-11-27 | 2023-08-22 | Robert F. Cruickshank, III | Optimized load shaping system, method and apparatus for optimizing production and consumption of energy |
| EP4066367A4 (en) * | 2019-11-27 | 2024-01-03 | Cruickshank, Robert F., III | OPTIMIZED LOAD SHAPING SYSTEM, METHOD AND APPARATUS FOR OPTIMIZING ENERGY PRODUCTION AND CONSUMPTION |
| CN112398176A (zh) * | 2020-11-05 | 2021-02-23 | 贵州电网有限责任公司 | 一种考虑燃煤机组启停调峰的水火风互济系统日前优化调度方法 |
| CN115400864A (zh) * | 2021-05-28 | 2022-11-29 | 北京京能电力股份有限公司 | 一种自动判断调峰火电机组磨煤机自动启停时机的方法 |
| CN113541204A (zh) * | 2021-08-10 | 2021-10-22 | 上海电力大学 | 一种燃煤发电机组全负荷调峰的自适应调节方法及系统 |
| CN114220560A (zh) * | 2021-12-13 | 2022-03-22 | 中国核动力研究设计院 | 一种压水堆核电厂硼控调峰方法及系统 |
| CN114220560B (zh) * | 2021-12-13 | 2024-05-07 | 中国核动力研究设计院 | 一种压水堆核电厂硼控调峰方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109347152B (zh) | 2022-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Energy management and optimization of vehicle-to-grid systems for wind power integration | |
| CN109347152A (zh) | 考虑多类型电源参与调峰的随机生产模拟方法及应用 | |
| CN107276127B (zh) | 考虑联络线电力交易计划的多区域互联电力系统的风电消纳优化方法 | |
| CN106099993B (zh) | 一种适应新能源大规模接入的电源规划方法 | |
| CN109767078B (zh) | 一种基于混合整数规划的多类型电源检修安排方法 | |
| CN109586284A (zh) | 考虑弃能约束的送端电力系统随机生产模拟方法及应用 | |
| CN109103926A (zh) | 基于多辐照特性年气象场景的光伏发电接纳能力计算方法 | |
| CN105207259B (zh) | 基于能量管理的并网状态下的微电网系统调度方法 | |
| CN113095791B (zh) | 一种综合能源系统运行方法及系统 | |
| CN106485352A (zh) | 一种多类电源发电计划组合决策方法 | |
| CN111614121A (zh) | 考虑需求响应的含电动汽车的多能源园区日前经济调度方法 | |
| CN112464477A (zh) | 计及需求响应的多能耦合综合能源运行仿真方法 | |
| CN110061524B (zh) | 一种基于深度神经网络的分布式电源虚拟电厂有功调度等值聚合方法及系统 | |
| CN108876000A (zh) | 一种光伏充电站光、储、变压器容量协调优化配置方法 | |
| CN111210079B (zh) | 一种分布式能源虚拟电厂运行优化方法及系统 | |
| CN114707289B (zh) | 基于机会约束的电热耦合综合能源系统多目标优化方法 | |
| CN108039737A (zh) | 一种源网荷协调运行模拟系统 | |
| CN110661301B (zh) | 一种水光蓄多能互补发电系统的容量配置优化方法 | |
| CN114595868A (zh) | 一种综合能源系统源网荷储协同规划方法及系统 | |
| CN115577929A (zh) | 基于多场景分析的乡村综合能源系统随机优化调度方法 | |
| CN108667077A (zh) | 一种风蓄联合系统优化调度方法 | |
| CN112134271A (zh) | 一种多能互补系统装机容量优化方法及系统 | |
| CN112510702B (zh) | 一种可再生能源发电制气的综合能源站规划方法及系统 | |
| CN108767855A (zh) | 一种时序持续混合的电力系统随机生产模拟方法 | |
| CN117114281A (zh) | 一种灵活资源多阶段规划方案的确定方法 |
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 |