CN106961125B - 风水火电力***动态经济调度的等式约束处理方法 - Google Patents
风水火电力***动态经济调度的等式约束处理方法 Download PDFInfo
- Publication number
- CN106961125B CN106961125B CN201710318716.XA CN201710318716A CN106961125B CN 106961125 B CN106961125 B CN 106961125B CN 201710318716 A CN201710318716 A CN 201710318716A CN 106961125 B CN106961125 B CN 106961125B
- Authority
- CN
- China
- Prior art keywords
- power
- thermal
- generating unit
- adjusted
- hydroelectric
- 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.)
- Expired - Fee Related
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 16
- 230000003828 downregulation Effects 0.000 claims description 30
- 230000003827 upregulation Effects 0.000 claims description 29
- 238000000034 method Methods 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 18
- 230000033228 biological regulation Effects 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 6
- 230000009194 climbing Effects 0.000 claims description 5
- 239000003245 coal Substances 0.000 claims description 5
- 229910052739 hydrogen Inorganic materials 0.000 claims description 5
- 239000001257 hydrogen Substances 0.000 claims description 5
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims description 4
- 238000004134 energy conservation Methods 0.000 claims description 4
- 230000007613 environmental effect Effects 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims 2
- 238000005457 optimization Methods 0.000 abstract description 16
- 239000002245 particle Substances 0.000 description 18
- 230000009467 reduction Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 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
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 150000002431 hydrogen Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000011160 research 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/06316—Sequencing of tasks or work
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
风水火电力***动态经济调度的等式约束处理方法,包括以下步骤:步骤一、建立风水火电力***动态经济调度模型中的等式约束,筛选出优化算法在当前迭代中产生的不满足等式约束的候选解,并计算这些候选解的等式约束违反量,确定***总功率的调整量;步骤二、设定***总功率的上调和下调原则,包括确定各发电机组的出力调整顺序以及***总功率的调整量的分配原则;步骤三、根据步骤二设定的***总功率上调和下调原则分配***总功率的调整量,得到满足等式约束的各发电机组的出力值。本发明解决了现有技术中存在的算法在求解电力***动态经济调度问题时因等式约束无法满足而影响求解效率和收敛性的问题。
Description
技术领域
本发明属于电力***经济调度技术领域,涉及一种风水火电力***动态经济调度的等式约束处理方法。
背景技术
电力***动态经济调度是一个数学优化问题,具有非线性、多约束、高维度、强耦合性等特点,求解非常复杂。含风水火的电力***动态经济调度由于涉及的调度对象较多,且不同能源电力的运行特性差异较大,使得求解难度和计算量剧增。目前已有大量智能优化算法用于电力***动态经济调度的求解,并且新算法的发展也是日新月异,比如遗传算法、粒子群算法、蚁群算法、微分进化算法等。这些算法本质都是从随机生成的一组初始解开始,遵循算法自身的寻优机制,通过一步步迭代更新逼近全局最优解,算法的初始值均是满足等式约束的解,然而经过迭代更新生成的解往往难以满足等式约束,造成大量无效解的产生,严重影响算法的求解效率,如何将无效解调整为可行解成为提高算法求解效率的重要手段,即实时功率平衡约束是含风水火的电力***动态经济调度模型中最难以严格满足的等式约束,也是影响影响算法求解效率和求解质量的关键因素。
而约束条件的处理方式直接影响到优化算法求解效率和调度结果。大部分用于电力***动态经济调度的优化算法对等式约束的处理方式是基于罚函数法,将约束条件以惩罚项的形式融入到目标函数中,将原问题转化为无约束的优化问题。对于严重违反等式约束的候选解,其对应的目标函数值很大,该解就成为无效解,算法将继续搜寻可行的候选解,因此大量无效解的产生严重影响了算法的求解效率,导致无法找到全局最优解。针对这一问题,文献(Sinha N,Chakrabarti R,Chattopadhyay PK.Evolutionary programmingtechniques for economic load dispatch[J].IEEE Transactions on EvolutionaryComputation,2003,7(1):83-94)在利用粒子群算法求解电力***经济调度时,提出一种处理等式约束的策略,即仅利用速度和位置更新公式生成前n-1维变量,剩余最后1维变量的取值则用来平衡等式约束的不满足量。这种方法很容易导致最后1维变量越出其可行域,生成的解仍是无效解。为解决该问题,文献(ParkJB,Lee KS,Shin JR,et al.A particleswarm optimization for economic dispatch with nonsmooth cost functions[J].IEEE Transactions on Power Systems,2005,20(1):34-42)提出当最后1维变量越出其可行域时,则重新返回去调整之前n-1维变量的取值,如此类推,又形成一个周而复始的循环过程。上述研究仅仅是针对含火电机组的电力***静态经济调度所提出的等式约束处理方法,对各机组功率的调整并没有合理有效的依据。因此,对于电力***特别是含风水火电的多源电力***,缺少合理有效且具有普适性的电力***动态经济调度的等式约束处理方法来提高优化算法的求解效率和求解质量。
发明内容
本发明的目的是提供一种风水火电力***动态经济调度的等式约束处理方法,解决了现有技术中存在的算法在求解电力***动态经济调度问题时因等式约束无法满足而影响求解效率和收敛性的问题。
本发明所采用的技术方案是,风水火电力***动态经济调度的等式约束处理方法,包括以下步骤:
步骤一、建立风水火电力***动态经济调度模型中的等式约束,筛选出求解风水火电力***动态经济调度问题的算法在当前迭代中产生的不满足等式约束的候选解,并计算这些候选解的等式约束违反量,判断需要上调还是下调***总功率,并确定***总功率的调整量;
步骤二、设定***总功率的上调和下调原则,包括确定各发电机组的出力调整顺序以及***总功率的调整量的分配原则;
步骤2.1、对风水火电力***中参与调度的各发电机组设定优先级;
步骤2.2、根据优先级确定各发电机组的出力调整顺序;
步骤2.3、根据各发电机组的功率调整限额分配***总功率的调整量;
步骤三、根据步骤二设定的***总功率上调和下调原则分配***总功率的调整量,得到满足等式约束的各发电机组的出力值。
步骤2.1的具体过程为:根据节能环保性和调节灵活性,设定水电机组的优先级高于火电机组的优先级;对每台火电机组根据经济性从高到低的顺序设定优先级,具体为:若***给定火电机组的单耗,则根据单耗由低到高的顺序设定优先级;若***给定火电机组的煤耗曲线,则根据最小比耗量由小到大的顺序设定优先级。
步骤2.2的具体过程为:当某一时段需要下调***总功率时,首先按优先级从低到高的顺序依次下调各火电机组的出力值,直到满足负荷需求,若所有火电机组的出力值调整完仍不能满足负荷需求,开始下调各水电机组的出力值,直到满足负荷需求;当某一时段需要上调***总功率时,首先上调各水电机组的出力值,直到满足负荷需求,若所有水电机组的出力值调整完仍不能满足负荷需求,开始按优先级从高到低的顺序依次上调各火电机组的出力值,直到满足负荷需求。
步骤2.3的具体过程为:对于火电机组,根据各火电机组的功率调整限额按优先级的高低分配***总功率的调整量;对于水电机组,根据各水电机组的功率调整限额按占比分配***总功率的调整量;对于风电场,遵循风电功率全额入网原则。
步骤一的具体步骤为:建立风水火电力***动态经济调度模型中的等式约束,即***实时功率平衡约束:
式中:Nt为火电机组总数;Nh为水电机组总数;Nw为风电场总数;PGit为t时段第i台火电机组的输出功率;PHjt为t时段第j台水电机组的输出功率;PWkt为t时段第k个风电场的所有风电机组出力之和;PDt为t时段的***负荷需求值;
筛选出求解风水火电力***动态经济调度问题的算法在当前迭代中产生的不满足等式约束的候选解,计算这些候选解的等式约束违反量,即***功率不平衡量ΔPDt:
若ΔPDt>0,需要下调***总功率,若ΔPDt<0,则需要上调***总功率,上调和下调时***总功率的调整量均为ΔPDt。
步骤2.3中,需要上调***总功率时,各发电机组的功率调整限额为功率上调限额,需要下调***总功率时,各发电机组的功率调整限额为功率下调限额;
功率上调限额为当前发电功率与最大输出功率的差额,分别计算火电机组的功率上调限额和水电机组的功率上调限额:
式中:和分别是第i台火电机组和第j台水电机组在t时段的功率上调限额;PGmaxi是第i台火电机组的最大输出功率;PHmaxj是第j台水电机组的最大输出功率;PGi(t-1)为t-1时段第i台火电机组的输出功率;rui是第i台火电机组增出力时的速率限值;Δt是t时段的时间长度;ruiΔt是爬坡限额;
功率下调限额为当前发电功率与最小技术出力的差额,分别计算火电机组的功率下调限额和水电机组的功率下调限额:
式中:P Git和P Hjt分别是第i台水电机组和第j台水电机组在t时段的功率下调限额;PGmini和PHminj分别是第i台水电机组和第j台水电机组的最小技术出力;rdi是第i台火电机组减出力时的速率限值;rdiΔt是爬坡限额。
步骤三的具体过程为:
当ΔPDt>0时,判断仅火电机组下调功率还是火电机组和水电机组共同下调功率,若仅需要s台火电机组下调功率即可满足等式约束,则按火电机组的功率下调限额从优先级低的火电机组开始下调功率,调整后满足等式约束的火电机组的出力值为:
若仅火电机组下调功率不足以满足负荷需求,则需要水电机组同时下调功率,依据各水电机组的功率下调限额占所有水电机组的功率下调限额的比例分配出力调整量,第j台水电机组的出力调整量ΔPHjt按下式计算:
式中:ΔPDt′是火电机组下调功率后剩余的负荷不平衡量;
调整后满足等式约束的水电机组的出力值为:
当ΔPDt<0时,判断仅水电机组上调功率还是火电机组和水电机组共同上调功率,若仅需水电机组上调功率,依据各水电机组的功率下调限额占所有水电机组的功率下调限额的比例分配出力调整量,第j台水电机组的出力调整量ΔPHjt按下式计算:
调整后满足等式约束的水电机组的出力值为:
若水电机组上调功率后仍不能满足负荷需求,则火电机组继续上调功率,计算火电机组需要上调的总功率:从优先级最高的火电机组开始依次上调功率,直到满足:
且
第1~s台火电机组就是需要上调功率的火电机组,调整后满足等式约束的火电机组的出力值为:
本发明的有益效果是:风水火电力***动态经济调度的等式约束处理方法,将算法在求解电力***动态经济调度中产生的不满足等式约束的候选解进行调整,使之满足等式约束,有效避免了大量无效解的产生,促进了算法的求解效率,对电力***动态经济调度和静态经济调度均具有普遍适用性:
(1)对不满足等式约束的候选解的调整过程一次性完成,避免了现有方法中周而复始的循环过程,能够提高算法的求解效率;
(2)充分考虑了参与调度的各电力能源的特性,对各发电机组按照节能优先的顺序设定优先级,根据机组的经济性和灵活性,为***总功率的调整量的分配提供了合理依据,克服了现有处理方法无依据可循的缺陷;
(3)依据经济性确定各发电机组调整出力的顺序和出力值,满足等式约束的同时,增加了经济性好的机组的负荷量,有利于降低整个***的煤耗,与电力***动态经济调度节能的优化目标相一致,因此调整的候选解更接近最优值,提高了算法优化结果的质量。
附图说明
图1是风水火电力***动态经济调度的等式约束处理方法的流程图;
图2是基于IEEE 24节点26机测试***改进的风水火电力***动态经济调度丰水期的结果图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
如图1所示,风水火电力***动态经济调度的等式约束处理方法,其特征在于,包括以下步骤:
步骤一、建立风水火电力***动态经济调度模型中的等式约束,筛选出求解风水火电力***动态经济调度问题的算法在当前迭代中产生的不满足等式约束的候选解,并计算这些候选解的等式约束违反量,判断需要上调还是下调***总功率,并确定***总功率的调整量;
步骤二、设定***总功率的上调和下调原则,包括确定各发电机组的出力调整顺序以及***总功率的调整量的分配原则;
步骤2.1、对风水火电力***中参与调度的各发电机组设定优先级;
步骤2.2、根据优先级确定各发电机组的出力调整顺序;
步骤2.3、根据各发电机组的功率调整限额分配***总功率的调整量;
步骤三、根据步骤二设定的***总功率上调和下调原则分配***总功率的调整量,得到满足等式约束的各发电机组的出力值。
步骤2.1的具体过程为:根据节能环保性和调节灵活性,设定水电机组的优先级高于火电机组的优先级;对每台火电机组根据经济性从高到低的顺序设定优先级,假设风水火电力***的动态经济调度周期为T,t时段参与调度的火电机组有Nt台,以***给定火电机组的煤耗曲线为例,计算各火电机组的最小比耗量μmin,最小比耗量μmin的计算公式如下:
式中:μmini为第i台火电机组的最小比耗量/$/(MW·h);ai、bi、ci是第i台火电机组的燃料费用系数;PG0i是与第i台火电机组的最小比耗量点对应的有功出力/MW,PG0i按照下式取值:
根据μmin由小到大的顺序为各火电机组设定优先级,序号为1~Nt。
步骤2.2的具体过程为:当某一时段需要下调***总功率时,首先按优先级从低到高的顺序依次下调各火电机组的出力值,直到满足负荷需求,若所有火电机组的出力值调整完仍不能满足负荷需求,开始下调各水电机组的出力值,直到满足负荷需求;当某一时段需要上调***总功率时,首先上调各水电机组的出力值,直到满足负荷需求,若所有水电机组的出力值调整完仍不能满足负荷需求,开始按优先级从高到低的顺序依次上调各火电机组的出力值,直到满足负荷需求。
步骤2.3的具体过程为:对于火电机组,根据各火电机组的功率调整限额按优先级的高低分配***总功率的调整量;对于水电机组,根据各水电机组的功率调整限额按占比分配***总功率的调整量;对于风电场,遵循风电功率全额入网原则,一般不参与调度,若出现火电机组和水电机组即使按最小技术出力仍高于负荷需求的极限,则按需弃掉一部分风电功率。
步骤一的具体步骤为:
步骤1.1、建立风水火电力***动态经济调度模型中的等式约束,即***实时功率平衡约束:
式中:Nt为火电机组总数;Nh为水电机组总数;Nw为风电场总数;PGit为t时段第i台火电机组的输出功率/MW;PHjt为t时段第j台水电机组的输出功率/MW;PWkt为t时段第k个风电场的所有风电机组出力之和/MW;PDt为t时段的***负荷需求值/MW;
本发明以粒子群算法求解风水火电力***动态经济调度问题为例:
步骤1.2、设置初始时刻t=1;
步骤1.3、更新粒子速度和粒子位置计算当前迭代的候选解,粒子速度更新公式为:粒子位置更新公式分别为:
其中:粒子是i在第k+1次迭代中第d维的速度;ω是惯性权重系数,用于平衡局部搜索能力和全局搜索能力;vk id是粒子i在第k次迭代中第d维的速度;c1、c2称为学习因子;r1、r2是[0,1]区间内的随机数;是粒子i在第k次迭代中第d维的个体极值点;xkid是粒子i在第k次迭代中第d维位置的当前值;是整个种群在第d维的全局极值点;xk+1id为粒子i在第k+1次迭代中第d维位置的当前值;
步骤1.4、筛选出粒子群算法在当前迭代中产生的不满足等式约束的候选解,计算这些候选解的等式约束违反量,即***功率不平衡量ΔPDt:
若ΔPDt>0,需要下调***总功率,若ΔPDt<0,则需要上调***总功率,上调和下调时***总功率的调整量均为ΔPDt。
步骤2.3中,需要上调***总功率时,各发电机组的功率调整限额为功率上调限额,需要下调***总功率时,各发电机组的功率调整限额为功率下调限额;
功率上调限额为当前发电功率与最大输出功率的差额,分别计算粒子每一维变量对应的火电机组的功率上调限额和水电机组的功率上调限额,火电机组的最大输出功率是考虑了爬坡限额后机组可输出的最大功率:
式中:和分别是第i台火电机组和第j台水电机组在t时段的功率上调限额;PGmaxi是第i台火电机组的最大输出功率/MW;PHmaxj是第j台水电机组的最大输出功率/MW;PGi(t-1)为t-1时段第i台火电机组的输出功率/MW;rui是第i台火电机组增出力时的速率限值/MW/min;Δt是t时段的时间长度;ruiΔt是爬坡限额;
功率下调限额为当前发电功率与最小技术出力的差额,分别计算粒子每一维变量对应的火电机组的功率下调限额和水电机组的功率下调限额,火电机组的最小技术出力是考虑了爬坡限额后机组可输出的最小功率:
式中:P Git和P Hjt分别是第i台水电机组和第j台水电机组在t时段的功率下调限额;PGmini和PHminj分别是第i台水电机组和第j台水电机组的最小技术出力/MW;rdi是第i台火电机组减出力时的速率限值/MW/min;rdiΔt是爬坡限额。
步骤三的具体过程为:
当ΔPDt>0时,当前发电机组出力之和大于实时负荷需求,首先火电机组下调功率,设s=0,用于记录下调功率的火电机组数目,递减变量ΔPDt′=ΔPDt,从优先级低的第Nt台火电机组开始,令递减变量其中是第Nt台火电机组的功率下调限额,若ΔPDt′>0且s<Nt,则Nt=Nt-1,s=s+1,重复按公式(4)计算粒子每一维变量对应的火电机组的功率下调限额和水电机组的功率下调限额;
若ΔPDt′≤0且s≤Nt,说明仅需要s台火电机组下调功率即可满足等式约束,按火电机组的功率下调限额从优先级低的火电机组开始下调功率,剩余的Nt-s台火电机组和水电机组保持原出力,调整后满足等式约束的火电机组的出力值为:
由公式(5)能够得出,对于选出的需要下调功率的s台火电机组,其中优先级低的s-1台火电机组的下调功率均取其功率下调限额,第s台火电机组下调剩余的功率盈余,其下调功率小于或者等于其功率下调限额;
若ΔPDt′>0且s=Nt,说明仅靠火电机组下调功率仍有功率盈余,需要水电机组同时下调功率,依据各水电机组的功率下调限额占所有水电机组的功率下调限额的比例分配出力调整量,第j台水电机组的出力调整量ΔPHjt按下式计算:
式中:ΔPDt′是火电机组下调功率后剩余的负荷不平衡量;
调整后满足等式约束的水电机组的出力值为:
当ΔPDt<0时,当前发电机组出力之和小于实时负荷需求,首先水电机组上调功率,若说明仅需水电机组上调功率,依据各水电机组的功率上调限额占所有水电机组的功率上调限额的比例分配出力调整量,第j台水电机组的出力调整量ΔPHjt按下式计算:
调整后满足等式约束的水电机组的出力值为:
若说明水电机组上调功率后仍不能满足负荷需求,需要火电机组上调功率,而每个水电机组的出力调整量均是其功率上调限额,因此调节后满足等式约束的水电机组出力为:
计算火电机组需要上调的总功率:从优先级最高的火电机组开始依次上调功率,直到满足:
且
第1~s台火电机组就是需要上调功率的火电机组,其中优先级靠前的s-1台火电机组的上调功率均为其功率上调限额,第s台火电机组上调剩余的功率盈余,上调功率小于或者等于其功率上调限额,调整后满足等式约束的火电机组的出力值为:
令t=t+1,若t≤T,返回步骤1.2,否则,结束调整,输出结果。
表1采用与未采用等式约束处理方法的26机测试***求解结果对比
图2和表1是基于IEEE 24节点26机测试***改进的含风水火电***动态经济调度丰水期的结果,调度时段为1天24小时,表1给出了采用和未采用本发明的等式约束处理方法得到的统计结果对比,在该算例中,粒子群算法的最大迭代次数取为100,独立运行50次,对***的最小燃煤费用、找到最优解的平均迭代次数和平均运行时间进行了对比,能够发现,采用本发明的等式约束处理方法,能够使算法求得的最优解更优,并且能以更少的迭代次数和运行时间找到最优解,从而提高算法求解效率和精度,该方法具有很好的普适性。
通过上述方式,本发明风水火电力***动态经济调度的等式约束处理方法,事先对参与调度的火电机组和水电机组分别依据经济性和灵活性进行排序,并设计兼顾经济性和快速性的功率上调和下调原则;然后针对优化算法在当前迭代中产生的不满足等式约束的候选解,计算其违反等式约束的量,确定功率上调量或下调量;最后根据设定好的调节原则,将等式约束违反量分配给相关发电机组。在调整过程中遵循节能、灵活的原则,与节能的优化目标相一致,使得调整后的候选解均满足等式约束,比没有调整依据产生的解更接近最优值,保证了算法寻优的顺利进行,提高了算法的求解精度和效率。
Claims (4)
1.风水火电力***动态经济调度的等式约束处理方法,其特征在于,包括以下步骤:
步骤一、建立风水火电力***动态经济调度模型中的等式约束,筛选出求解风水火电力***动态经济调度问题的算法在当前迭代中产生的不满足等式约束的候选解,并计算这些候选解的等式约束违反量,判断需要上调还是下调***总功率,并确定***总功率的调整量,所述步骤一的具体步骤为:
建立风水火电力***动态经济调度模型中的等式约束,即***实时功率平衡约束:
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>W</mi>
<mi>k</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:Nt为火电机组总数;Nh为水电机组总数;Nw为风电场总数;PGit为t时段第i台火电机组的输出功率;PHjt为t时段第j台水电机组的输出功率;PWkt为t时段第k个风电场的所有风电机组出力之和;PDt为t时段的***负荷需求值;
筛选出求解风水火电力***动态经济调度问题的算法在当前迭代中产生的不满足等式约束的候选解,计算这些候选解的等式约束违反量,即***功率不平衡量ΔPDt:
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>w</mi>
</msub>
</munderover>
<msub>
<mi>P</mi>
<mrow>
<mi>W</mi>
<mi>k</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
若ΔPDt>0,需要下调***总功率,若ΔPDt<0,则需要上调***总功率,上调和下调时***总功率的调整量均为ΔPDt;
步骤二、设定***总功率的上调和下调原则,包括确定各发电机组的出力调整顺序以及***总功率的调整量的分配原则;
步骤2.1、对风水火电力***中参与调度的各发电机组设定优先级;
步骤2.2、根据优先级确定各发电机组的出力调整顺序;
步骤2.3、根据各发电机组的功率调整限额分配***总功率的调整量;所述步骤2.3中,需要上调***总功率时,各发电机组的功率调整限额为功率上调限额,需要下调***总功率时,各发电机组的功率调整限额为功率下调限额;
功率上调限额为当前发电功率与最大输出功率的差额,分别计算火电机组的功率上调限额和水电机组的功率上调限额:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>max</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>u</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>max</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:和分别是第i台火电机组和第j台水电机组在t时段的功率上调限额;PGmaxi是第i台火电机组的最大输出功率;PHmaxj是第j台水电机组的最大输出功率;PGi(t-1)为t-1时段第i台火电机组的输出功率;rui是第i台火电机组增出力时的速率限值;Δt是t时段的时间长度;ruiΔt是爬坡限额;
功率下调限额为当前发电功率与最小技术出力的差额,分别计算火电机组的功率下调限额和水电机组的功率下调限额:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>max</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>min</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>r</mi>
<mrow>
<mi>d</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>min</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:P Git和P Hjt分别是第i台水电机组和第j台水电机组在t时段的功率下调限额;PGmini和PHminj分别是第i台水电机组和第j台水电机组的最小技术出力;rdi是第i台火电机组减出力时的速率限值;rdiΔt是爬坡限额;
步骤三、根据所述步骤二设定的***总功率上调和下调原则分配***总功率的调整量,得到满足等式约束的各发电机组的出力值;所述步骤三的具体过程为:
当ΔPDt>0时,判断仅火电机组下调功率还是火电机组和水电机组共同下调功率,若仅需要s台火电机组下调功率即可满足等式约束,则按火电机组的功率下调限额从优先级低的火电机组开始下调功率,调整后满足等式约束的火电机组的出力值为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>2</mn>
</mrow>
</munderover>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mi>s</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
若仅火电机组下调功率不足以满足负荷需求,则需要水电机组同时下调功率,依据各水电机组的功率下调限额占所有水电机组的功率下调限额的比例分配出力调整量,第j台水电机组的出力调整量ΔPHjt按下式计算:
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&prime;</mo>
</msup>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
</munderover>
<msub>
<munder>
<mi>P</mi>
<mo>&OverBar;</mo>
</munder>
<mrow>
<mi>H</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:ΔPDt′是火电机组下调功率后剩余的负荷不平衡量;
调整后满足等式约束的水电机组的出力值为:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
当ΔPDt<0时,判断仅水电机组上调功率还是火电机组和水电机组共同上调功率,若仅需水电机组上调功率,依据各水电机组的功率下调限额占所有水电机组的功率下调限额的比例分配出力调整量,第j台水电机组的出力调整量ΔPHjt按下式计算:
<mrow>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mfrac>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
</munderover>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>H</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
调整后满足等式约束的水电机组的出力值为:
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>H</mi>
<mi>j</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>N</mi>
<mi>h</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
若水电机组上调功率后仍不能满足负荷需求,则火电机组继续上调功率,计算火电机组需要上调的总功率:从优先级最高的火电机组开始依次上调功率,直到满足:
第1~s台火电机组就是需要上调功率的火电机组,调整后满足等式约束的火电机组的出力值为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>s</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;P</mi>
<mrow>
<mi>D</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>s</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mover>
<mi>P</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>G</mi>
<mi>i</mi>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>s</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
2.根据权利要求1所述的风水火电力***动态经济调度的等式约束处理方法,其特征在于,所述步骤2.1的具体过程为:
根据节能环保性和调节灵活性,设定水电机组的优先级高于火电机组的优先级;对每台火电机组根据经济性从高到低的顺序设定优先级,具体为:若***给定火电机组的单耗,则根据单耗由低到高的顺序设定优先级;若***给定火电机组的煤耗曲线,则根据最小比耗量由小到大的顺序设定优先级。
3.根据权利要求1所述的风水火电力***动态经济调度的等式约束处理方法,其特征在于,所述步骤2.2的具体过程为:
当某一时段需要下调***总功率时,首先按优先级从低到高的顺序依次下调各火电机组的出力值,直到满足负荷需求,若所有火电机组的出力值调整完仍不能满足负荷需求,开始下调各水电机组的出力值,直到满足负荷需求;当某一时段需要上调***总功率时,首先上调各水电机组的出力值,直到满足负荷需求,若所有水电机组的出力值调整完仍不能满足负荷需求,开始按优先级从高到低的顺序依次上调各火电机组的出力值,直到满足负荷需求。
4.根据权利要求1所述的风水火电力***动态经济调度的等式约束处理方法,其特征在于,所述步骤2.3的具体过程为:
对于火电机组,根据各火电机组的功率调整限额按优先级的高低分配***总功率的调整量;对于水电机组,根据各水电机组的功率调整限额按占比分配***总功率的调整量;对于风电场,遵循风电功率全额入网原则。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710318716.XA CN106961125B (zh) | 2017-05-08 | 2017-05-08 | 风水火电力***动态经济调度的等式约束处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710318716.XA CN106961125B (zh) | 2017-05-08 | 2017-05-08 | 风水火电力***动态经济调度的等式约束处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106961125A CN106961125A (zh) | 2017-07-18 |
CN106961125B true CN106961125B (zh) | 2018-04-13 |
Family
ID=59482057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710318716.XA Expired - Fee Related CN106961125B (zh) | 2017-05-08 | 2017-05-08 | 风水火电力***动态经济调度的等式约束处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106961125B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108155674B (zh) * | 2018-02-01 | 2020-02-21 | 清华大学 | 考虑不确定分布特性的水火电联合调度方法及*** |
CN108400617B (zh) * | 2018-03-19 | 2021-05-28 | 燕山大学 | 一种经济调度的约束处理方法 |
CN111638693B (zh) * | 2020-06-04 | 2021-04-27 | 杭州和利时自动化有限公司 | 一种多个设备控制方法及*** |
CN114693343A (zh) * | 2022-03-14 | 2022-07-01 | 北京明略软件***有限公司 | 一种广告预算分配方法、装置、设备及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102779322A (zh) * | 2012-06-29 | 2012-11-14 | 西安理工大学 | 基于环境保护和经济效益的火电厂综合调度方法 |
CN103246263A (zh) * | 2013-04-22 | 2013-08-14 | 天津大学 | 冷热电联供微网***的通用优化调度策略 |
CN104036320A (zh) * | 2014-05-30 | 2014-09-10 | 国家电网公司 | 一种基于改进粒子群算法的微网***动态经济调度方法 |
-
2017
- 2017-05-08 CN CN201710318716.XA patent/CN106961125B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102779322A (zh) * | 2012-06-29 | 2012-11-14 | 西安理工大学 | 基于环境保护和经济效益的火电厂综合调度方法 |
CN103246263A (zh) * | 2013-04-22 | 2013-08-14 | 天津大学 | 冷热电联供微网***的通用优化调度策略 |
CN104036320A (zh) * | 2014-05-30 | 2014-09-10 | 国家电网公司 | 一种基于改进粒子群算法的微网***动态经济调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106961125A (zh) | 2017-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107994595B (zh) | 一种削峰填谷控制方法和***以及应用该控制方法的*** | |
CN109599861B (zh) | 考虑本地负荷调峰能力的送端电网电源结构规划方法 | |
CN106961125B (zh) | 风水火电力***动态经济调度的等式约束处理方法 | |
CN109347151B (zh) | 一种新能源参与调峰的送端电网电源结构优化方法 | |
CN109345019B (zh) | 一种基于改进粒子群算法的微电网经济调度优化策略 | |
CN107039977B (zh) | 电力***鲁棒调度不确定集构建方法 | |
CN111555355B (zh) | 水光蓄联合发电的调度策略及其优化方法 | |
CN109390973B (zh) | 一种考虑通道约束的送端电网电源结构优化方法 | |
CN106953316B (zh) | 微电网变时间尺度优化调度方法 | |
CN109167347B (zh) | 基于云自适应粒子群多目标电动汽车充放电优化调度方法 | |
CN108471143A (zh) | 基于正负反馈粒子群算法的微电网多能源调度优化方法 | |
CN109146320B (zh) | 一种考虑配电网安全性的虚拟电厂优化调度方法 | |
CN106877339B (zh) | 一种考虑电动汽车接入配电网后随机模糊潮流的分析方法 | |
CN109447393A (zh) | 一种电力***经济负荷分配的改进的粒子群算法 | |
CN109861302B (zh) | 一种基于主从博弈的能源互联网日前优化控制方法 | |
CN113783224A (zh) | 一种考虑多种分布式能源运行的配电网双层优化规划方法 | |
CN105184426B (zh) | 一种基于随机连续寻优策略的梯级电站调峰方法 | |
Tian et al. | Coordinated planning with predetermined renewable energy generation targets using extended two-stage robust optimization | |
CN108808737A (zh) | 促进可再生分布式电源消纳的主动配电网优化调度方法 | |
CN115117931A (zh) | 考虑电动汽车灵活性和光伏接入的配电网规划方法及*** | |
CN111404204A (zh) | 一种利用改进狼群算法的含风电场电力***优化调度方法 | |
CN114819373A (zh) | 一种基于合作博弈的共享混合储能电站的储能规划方法 | |
CN116090739A (zh) | 一种综合考虑电动汽车有序充电的虚拟电厂优化调度方法 | |
CN117578537A (zh) | 基于碳交易和需求响应的微电网优化调度方法 | |
CN103679284A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180413 Termination date: 20210508 |
|
CF01 | Termination of patent right due to non-payment of annual fee |