CN107910866A - 一种考虑需求侧响应不确定性的电力系统日前优化调度方法 - Google Patents

一种考虑需求侧响应不确定性的电力系统日前优化调度方法 Download PDF

Info

Publication number
CN107910866A
CN107910866A CN201711218112.4A CN201711218112A CN107910866A CN 107910866 A CN107910866 A CN 107910866A CN 201711218112 A CN201711218112 A CN 201711218112A CN 107910866 A CN107910866 A CN 107910866A
Authority
CN
China
Prior art keywords
demand
price
curve
uncertainty
node
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
Application number
CN201711218112.4A
Other languages
English (en)
Other versions
CN107910866B (zh
Inventor
黎华盛
乔静
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jingdezhen Power Supply Branch Jiangxi Electric Power Co Ltd
Original Assignee
Jingdezhen Power Supply Branch Jiangxi Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jingdezhen Power Supply Branch Jiangxi Electric Power Co Ltd filed Critical Jingdezhen Power Supply Branch Jiangxi Electric Power Co Ltd
Priority to CN201711218112.4A priority Critical patent/CN107910866B/zh
Publication of CN107910866A publication Critical patent/CN107910866A/zh
Application granted granted Critical
Publication of CN107910866B publication Critical patent/CN107910866B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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
    • 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/06315—Needs-based resource requirements planning or analysis
    • 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
    • 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
    • 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
    • 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]
    • 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • 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)
  • Entrepreneurship & Innovation (AREA)
  • Tourism & Hospitality (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Physics & Mathematics (AREA)
  • Educational Administration (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Power Engineering (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明属于电网调度领域，具体是一种考虑需求侧响应不确定性的电力系统日前优化调度方法，用来求解含新能源的电力系统动态经济调度问题。首先基于价格弹性需求曲线，以社会福利最大为优化目标，同时考虑需求侧响应不确定性和系统的实际运行约束，构建了需求侧间接控制的电力系统优化调度模型；然后，针对模型求解复杂度高的问题，提出了一种适用于不确定性环境的序优化算法对模型进行求解。基于IEEE30节点算例的仿真结果表明，本发明提出的调度方法，由于充分考虑了源荷侧调度资源的不确定性，有效降低了电力系统的运行成本，提升其运行精细化水平。

Description

一种考虑需求侧响应不确定性的电力系统日前优化调度方法
技术领域
本发明一种考虑需求侧响应不确定性的电力系统日前优化调度方法，涉及电力系统调度领域。
背景技术
风电及光伏是无污染、绿色的可再生能源，其分布广泛、能量密度高，适合大规模开发，因此，风力和光伏发电技术已受到了世界各国的高度重视。新能源的大规模并网虽然能在一定程度上缓解环境压力和能源危机，但由于其大都具有随机性、波动性和出力不可控的特点，因而其大规模接入会对电网的安全稳定运行带来一定的影响。随着智能电网建设的稳步推进，系统内的需求侧响应手段逐步增加，从而为平抑系统内新能源出力波动性提供了一种高效、廉价的调度资源。因此，充分考虑需求侧的调度资源，研究基于源荷互动的电力系统日前调度模式已成为当前学术界的研究热点。
发明内容
针对现有方法的不足，本发明提出了一种考虑需求侧响应不确定性的日前优化调度方法。首先基于价格弹性需求曲线，以社会福利最大为优化目标，同时考虑需求侧响应不确定性和系统的实际运行约束，构建了需求侧间接控制的电力系统优化调度模型；然后，针对模型求解复杂度高的问题，提出了一种适用于不确定性环境的序优化算法对模型进行求解。基于IEEE30节点算例的仿真结果表明，本发明提出的调度方法，由于充分考虑了源荷侧调度资源的不确定性，有效降低了电力系统的运行成本，提升其运行精细化水平。
本发明采取的技术方案为：
一种考虑需求侧响应不确定性的电力系统日前优化调度方法，包括以下步骤：
步骤1：源荷侧不确定性建模
步骤1.1：需求侧响应的不确定性建模
本发明针对基于电价的需求侧响应，根据电价与电能需求量之间的变化关系，考虑价格弹性需求曲线不确定性，如图1所示，在此基础上建立需求侧响应模型。
价格弹性需求曲线的数学表达式为：
或者
式中：表示用于描述价格弹性需求曲线的不确定性的偏差。按照近似价格弹性需求的步骤，曲线作为分段函数，对于价格弹性需求曲线中的每个电价相应的电能需求量允许在范围内变化，其中表示基准电价参考值，是的偏差，是的上限。
需求响应曲线的不确定性集可以描述如下：
式中：是的偏差，是的上限。T表示时间段集，B表示总线，K表示引入的近似价格性需求曲线。
步骤1.2：风电出力的不确定性建模
本发明利用基数不确定性集来调整风电出力模型的保守性。引入整数πb作为基数预算，以限制风力输出远离其在总线b处的预测值的时间段的数量。在每个总线b处，当风力输出达到其上限、下限或预测值以及风力输出不处于其预测值的周期的总数时，发生最坏情况的风力输出情形应该不超过预算值πb。
不确定性集合可以描述如下：
式中：代表在时间段t内总线b的风力输出的预测值，Wt b+Wt b-分别代表高于和低于允许的最大偏差值，和是二进制变量。如果风力输出将达到其上限，如果风力输出将达到其下限，如果两者都为0，就会实现预测值。
步骤2：考虑不确定性需求侧响应的日前调度建模
步骤2.1：确定性需求侧响应的日前优化调度建模
本发明首先建立确定性需求侧响应的日前调度模型。由于在日前调度优化模型考虑了电价因素，因此目标函数为社会福利最大化。此外在该模型中，风电出力被假设为确定值，价格弹性需求曲线也是确定的。模拟需求曲线和供给曲线，如图2所示。
目标函数描述如下：
式中：T表示时间段的集合，B表示节点的集合，Gb表示节点b处的发电机组，表示发电机i在节点b的启动成本，为发电机i在节点b的停机成本，为发电机i在节点b处在时间段t内产生的电量，为节点b在时间段t内的实际用电需求，为节点b在时间段t内的价格弹性需求曲线的积分，fi b为发电机i在节点b处的燃料成本函数。
为保证电力系统能够安全可靠地运行，决策变量还需满足以下常规约束条件：
(1)机组技术约束条件
式中:表示发电机i在节点b处的最小正常运行时间，表示发电机i在节点b处的最小停机时间，表示发电机i在节点b的最小发电量，为发电机i在节点b的最大发电量,为二进制变量，用于指示在时间段t中发电机i是否在节点b上，为二进制变量，用于指示发电机i是否在时间段t内在节点b上启动，为二进制变量，用于指示发电机i在时间段t是否在节点b处关闭。
(2)系统约束条件
系统的约束条件包括：系统的功率平衡约束和线路的传输容量约束。其数学表达如下式所示：
式中:Ω为连接两个节点的传输线路，为节点b在时间段t内的实际用电需求，为连接节点i和节点b的传输线的传输容量。
(3)考虑需求响应后新添加的约束条件。
式中:为节点b在时间段t内的非弹性需求部分，为节点b在时间段t内的最大需求。
步骤2.2：线性化机组燃料成本函数
本发明将机组的燃料成本表示为二次函数，使用N段线性函数来近似燃料成本函数：
式中:和是第j段函线的截距和斜率，是辅助变量。
步骤2.3：线性化消费者需求响应模型
本发明将消费者剩余中的弹性用电需求响应模型转化为价格弹性需求曲线模型，如果价格弹性是恒定的，可以将价格弹性需求曲线表示为：
式中：是节点b在时间段t内的电价，是t时间段中节点b处所给定的价格弹性值，并且使可以由给定参考点决定的参数值。
除此之外，本发明所提出的解决方法可以应用到其他弹性需求的建模中，对于一些负荷聚合商或负荷代理商，价格弹性需求曲线本身是一个分段函数，如图3所示，对于一般的价格弹性需求曲线，可以应用分段函数来近似该价格弹性需求曲线，因此可以将近似为：
其中:是分段函数的第k段，是k段处的相应价格，是在k段处处为需求引入的辅助变量，K是所有的段长集合。
由于本发明是将最大化，所以是严格地随着k的增加而减少，因此可以得到：
当存在某个s0使得成立时，可以证明是价格弹性需求曲线的近似积分，即是合理的。
步骤2.4：不确定性需求侧响应的日前调度建模
在步骤2.1中，假设价格弹性需求曲线是确定的。当作出日前调度决策时，必须允许价格弹性需求曲线在一定范围内变化，实际的价格弹性需求曲线是不确定的。为了调整保守性，引入参数以限制偏差的总量，即可以通过改变的值来调整所提出的方法的保守性。值越小，需求响应曲线的不确定性也就越小。
需求响应曲线的不确定性集可以描述如下：
目标函数为：
约束条件如下：
在上述公式中,(23)中的目标函数是使得社会福利最大化(没有常数部分)。约束(24)、(25)分别表示最小的开机时间和最小的停机时间限制。约束条件(26)和(27)计算机组的启动和关闭状态变量。约束(28)强制执行每台发电机组的功率输出上限和下限。约束(29)和(30)实施每台机组的斜率限制。约束(31)确保负载平衡，并要求电源满足需求。约束(32)是传输线容量限制。最后，由于弹性需求的贡献，约束(33)强制实际需求的下限和上限。
步骤3：模型求解
本发明提出了一种序优化的求解思想，基于较高的概率寻求足够好的解集来对模型进行求解。
步骤3.1：按照机会均等原则，从总的决策空间Θ里挑选出N个决策，并构成有限集合ΘN。
步骤3.2：利用某种选择规则，通过快速粗糙评估，从N个决策中确定OPC类型。利用选择规则选出一子决策集S，其元素的个数为s。子决策集S将包含k个足够好的解(以至少α的概率)，相关的数学公式如下所示：
P{|G∩S|>k}≥α (35)
式中：G表示足够好的解的集合，解集中元素的个数为g，表示将各个决策经精确评估后的前g个最好的解；||表示集合元素的个数。
步骤3.3：精确评估子决策集中的每一个元素，从s个元素中选择出一个最好的解。
通过上述步骤，完成考虑需求侧响应不确定性的电力系统日前优化调度。
本发明一种考虑需求侧响应不确定性的电力系统日前优化调度方法，技术效果如下：
1)、本发明基于价格弹性需求曲线，以社会福利最大为优化目标，同时考虑需求侧响应的不确定性和系统的实际运行约束，构建计及需求侧间接控制的电力系统优化调度模型，利用一种不确定性的序优化求解算法对模型进行求解，降低了模型求解复杂度。
2)、本发明提出的调度方法，由于充分考虑了源荷侧调度资源的不确定性，可以有效提升系统的运行效益，提升其精细化水平。
附图说明
图1为本发明中价格弹性需求曲线不确定性描述图。
图2为本发明中需求曲线与供给曲线模拟图。
图3为本发明中分段函数近似价格弹性需求曲线图。
具体实施方式
一种考虑需求侧响应不确定性的电力系统日前优化调度方法，包括以下步骤：
步骤1：源荷侧不确定性建模
步骤1.1：需求侧响应的不确定性建模
本发明针对基于电价的需求侧响应，根据电价与电能需求量之间的变化关系，考虑价格弹性需求曲线不确定性，如图1所示，在此基础上建立需求侧响应模型。
价格弹性需求曲线的数学表达式为：
或者
式中：表示用于描述价格弹性需求曲线的不确定性的偏差。按照近似价格弹性需求的步骤，曲线作为分段函数，对于价格弹性需求曲线中的每个电价相应的电能需求量允许在范围内变化，其中表示基准电价参考值，是的偏差，是的上限。
需求响应曲线的不确定性集可以描述如下：
式中：是的偏差，是的上限。T表示时间段集，B表示总线，K表示引入的近似价格性需求曲线。
步骤1.2：风电出力的不确定性建模
本发明利用基数不确定性集来调整风电出力模型的保守性。引入整数πb作为基数预算，以限制风力输出远离其在总线b处的预测值的时间段的数量。在每个总线b处，当风力输出达到其上限、下限或预测值以及风力输出不处于其预测值的周期的总数时，发生最坏情况的风力输出情形应该不超过预算值πb。
不确定性集合可以描述如下：
式中：代表在时间段t内总线b的风力输出的预测值，Wt b+Wt b-分别代表高于和低于允许的最大偏差值，和是二进制变量。如果风力输出将达到其上限，如果风力输出将达到其下限，如果两者都为0，就会实现预测值。
步骤2：考虑不确定性需求侧响应的日前调度建模
步骤2.1：确定性需求侧响应的日前优化调度建模
本发明首先建立确定性需求侧响应的日前调度模型。由于在日前调度优化模型考虑了电价因素，因此目标函数为社会福利最大化。此外在该模型中，风电出力被假设为确定值，价格弹性需求曲线也是确定的。模拟需求曲线和供给曲线，如图2所示。
目标函数描述如下：
式中：T表示时间段的集合，B表示节点的集合，Gb表示节点b处的发电机组，表示发电机i在节点b的启动成本，为发电机i在节点b的停机成本，为发电机i在节点b处在时间段t内产生的电量，为节点b在时间段t内的实际用电需求，为节点b在时间段t内的价格弹性需求曲线的积分，fi b为发电机i在节点b处的燃料成本函数。
为保证电力系统能够安全可靠地运行，决策变量还需满足以下常规约束条件：
(1)机组技术约束条件
式中:表示发电机i在节点b处的最小正常运行时间，表示发电机i在节点b处的最小停机时间，表示发电机i在节点b的最小发电量，为发电机i在节点b的最大发电量,为二进制变量，用于指示在时间段t中发电机i是否在节点b上，为二进制变量，用于指示发电机i是否在时间段t内在节点b上启动，为二进制变量，用于指示发电机i在时间段t是否在节点b处关闭。
(2)系统约束条件
系统的约束条件包括：系统的功率平衡约束和线路的传输容量约束。其数学表达如下式所示：
式中:Ω为连接两个节点的传输线路，为节点b在时间段t内的实际用电需求，为连接节点i和节点b的传输线的传输容量。
(3)考虑需求响应后新添加的约束条件。
式中:为节点b在时间段t内的非弹性需求部分，为节点b在时间段t内的最大需求。
步骤2.2：线性化机组燃料成本函数
本发明将机组的燃料成本表示为二次函数，使用N段线性函数来近似燃料成本函数：
式中:和是第j段函线的截距和斜率，是辅助变量。
步骤2.3：线性化消费者需求响应模型
本发明将消费者剩余中的弹性用电需求响应模型转化为价格弹性需求曲线模型，如果价格弹性是恒定的，可以将价格弹性需求曲线表示为：
式中：是节点b在时间段t内的电价，是t时间段中节点b处所给定的价格弹性值，并且使可以由给定参考点决定的参数值。
除此之外，本发明所提出的解决方法可以应用到其他弹性需求的建模中，对于一些负荷聚合商或负荷代理商，价格弹性需求曲线本身是一个分段函数，如图3所示，对于一般的价格弹性需求曲线，可以应用分段函数来近似该价格弹性需求曲线，因此可以将近似为：
其中:是分段函数的第k段，是k段处的相应价格，是在k段处处为需求引入的辅助变量，K是所有的段长集合。
由于本发明是将最大化，所以是严格地随着k的增加而减少，因此可以得到：
当存在某个s0使得成立时，可以证明是价格弹性需求曲线的近似积分，即是合理的。
步骤2.4：不确定性需求侧响应的日前调度建模
在步骤2.1中，假设价格弹性需求曲线是确定的。当作出日前调度决策时，必须允许价格弹性需求曲线在一定范围内变化，实际的价格弹性需求曲线是不确定的。为了调整保守性，引入参数以限制偏差的总量，即可以通过改变的值来调整所提出的方法的保守性。值越小，需求响应曲线的不确定性也就越小。
需求响应曲线的不确定性集可以描述如下：
目标函数为：
约束条件如下：
在上述公式中,(23)中的目标函数是使得社会福利最大化(没有常数部分)。约束(24)、(25)分别表示最小的开机时间和最小的停机时间限制。约束条件(26)和(27)计算机组的启动和关闭状态变量。约束(28)强制执行每台发电机组的功率输出上限和下限。约束(29)和(30)实施每台机组的斜率限制。约束(31)确保负载平衡，并要求电源满足需求。约束(32)是传输线容量限制。最后，由于弹性需求的贡献，约束(33)强制实际需求的下限和上限。
步骤3：模型求解
本发明提出了一种序优化的求解思想，基于较高的概率寻求足够好的解集来对模型进行求解。
步骤3.1：按照机会均等原则，从总的决策空间Θ里挑选出N个决策，并构成有限集合ΘN。
步骤3.2：利用某种选择规则，通过快速粗糙评估，从N个决策中确定OPC类型。利用选择规则选出一子决策集S，其元素的个数为s。子决策集S将包含k个足够好的解(以至少α的概率)，相关的数学公式如下所示：
P{|G∩S|>k}≥α (35)
式中：G表示足够好的解的集合，解集中元素的个数为g，表示将各个决策经精确评估后的前g个最好的解；||表示集合元素的个数。
步骤3.3：精确评估子决策集中的每一个元素，从s个元素中选择出一个最好的解。
步骤4：确立运行模式
为对比分析本发明所建调度模型的有效性与正确性，确立以下两种运行模式：
模式1：不同弹性值下考虑需求侧响应不确定性的电力系统日前优化调度。
模式2：不同偏差值下考虑需求侧响应不确定性的电力系统日前优化调度。
本发明以修改的IEEE-30节点系统来验证所提模型的正确性。该系统包含30个节点、6台发电机，调度周期为24小时，使用四段分段线性函数近似机组燃料成本函数，基准电价参考值定为25$/MWh。所有的实验使用CPLEX 12.1，相关计算均在英特尔酷睿i3-3240处理器3.40GHz，8GB内存的Intel Quad Core 2.40GHz的计算机上实现。
为了比较不同需求响应场景下系统的运行效益，本发明计算了不同弹性值下(α＝0的时候，既为不考虑需求侧响应)，基于本发明中的方法得出运行方案的总社会福利，对比结果如表1所示：
表1不同需求响应方案对比
由表1可知，不同需求响应场景下的总生产成本、单位负荷成本和总的社会福利均有不同。总的生产成本和总的社会福利随着弹性值的增加而上升，同时单位负荷成本随着弹性值的增加呈现下降趋势。
可见，引入需求响应有助于降低系统的单位负荷成本，虽然会使得总的生产成本有所增加，但总的社会福利却是增大的。这是由于，将需求侧响应纳入到调度体系之中，有利于平抑源荷侧的波动性因素，提升其运行效益。
为了对比不同的需求响应方案，本发明假设需求弹性曲线是确定的，并将设置为5，但在实际生活中，由于多种原因，实际的价格弹性需求曲线实质上是不确定的。将弹性值固定为α＝-1，在此情景下本发明测试价格弹性需求曲线的三个不同偏差值，和为了比较不同偏差值下的性能，利用单位负荷成本和社会福利值来进行对比，其结果如下表2所示：
表2需求响应不确定性对比
由表2可知，当有不同偏差值时，且当系统有更多的不确定性需求时，即偏差值增大时，社会福利值和单位负荷成本均减小，其中总生产成本也呈现减少的趋势。分析可知，当需求响应曲线具有更多的不确定时，也就是说偏差值增大时，需求曲线移动以实现较小的需求平衡，这就对应于较小的总生产成本和单位负荷成本。
本发明按照优选实施例进行了说明，但上述实施例不以任何形式限定本发明，凡采用等同替换或等效变换的形式所获得的技术方案，均落在本发明技术方案的范围内。

Claims (1)

1.一种考虑需求侧响应不确定性的电力系统日前优化调度方法，其特征在于，包括以下步骤：
步骤1：源荷侧不确定性建模
步骤1.1：需求侧响应的不确定性建模
针对基于电价的需求侧响应，根据电价与电能需求量之间的变化关系，考虑价格弹性需求曲线不确定性，在此基础上建立需求侧响应模型，
价格弹性需求曲线的数学表达式为：
<mrow> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>=</mo> <msubsup> <mi>A</mi> <mi>t</mi> <mi>b</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>&amp;alpha;</mi> <mi>t</mi> <mi>b</mi> </msubsup> </msup> <mo>+</mo> <msubsup> <mi>&amp;epsiv;</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
或者
<mrow> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>=</mo> <msubsup> <mi>A</mi> <mi>t</mi> <mi>b</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;epsiv;</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>&amp;alpha;</mi> <mi>t</mi> <mi>b</mi> </msubsup> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中：表示用于描述价格弹性需求曲线的不确定性的偏差，按照近似价格弹性需求的步骤，曲线作为分段函数，对于价格弹性需求曲线中的每个电价相应的电能需求量允许在范围内变化，其中表示基准电价参考值，是的偏差，是的上限；
需求响应曲线的不确定性集可以描述如下：
<mrow> <mi>&amp;Pi;</mi> <mo>=</mo> <mo>{</mo> <mi>&amp;epsiv;</mi> <mo>:</mo> <mo>-</mo> <msubsup> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>&amp;epsiv;</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msubsup> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>k</mi> <mo>&amp;Element;</mo> <mi>K</mi> <mo>}</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中：是的偏差，是的上限。T表示时间段集，B表示总线，K表示引入的近似价格性需求曲线；
步骤1.2：风电出力的不确定性建模
利用基数不确定性集来调整风电出力模型的保守性，引入整数πb作为基数预算，以限制风力输出远离其在总线b处的预测值的时间段的数量，在每个总线b处，当风力输出达到其上限、下限或预测值以及风力输出不处于其预测值的周期的总数时，发生最坏情况的风力输出情形应该不超过预算值πb；
不确定性集合可以描述如下：
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>W</mi> <mo>:</mo> <mo>=</mo> <mo>{</mo> <mi>w</mi> <mo>&amp;Element;</mo> <msup> <mi>R</mi> <mrow> <mo>|</mo> <mi>B</mi> <mo>|</mo> <mo>&amp;times;</mo> <mo>|</mo> <mi>T</mi> <mo>|</mo> </mrow> </msup> <mo>:</mo> <msubsup> <mi>w</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>*</mo> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>z</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>+</mo> </mrow> </msubsup> <msubsup> <mi>W</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>+</mo> </mrow> </msubsup> <mo>-</mo> <msubsup> <mi>z</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>-</mo> </mrow> </msubsup> <msubsup> <mi>W</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>-</mo> </mrow> </msubsup> <mo>,</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> </mrow> </munder> <mrow> <mo>(</mo> <msubsup> <mi>z</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>+</mo> </mrow> </msubsup> <mo>+</mo> <msubsup> <mi>z</mi> <mi>t</mi> <mrow> <mi>b</mi> <mo>-</mo> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <msup> <mi>&amp;pi;</mi> <mi>b</mi> </msup> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>}</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中：代表在时间段t内总线b的风力输出的预测值，分别代表高于和低于允许的最大偏差值，和是二进制变量；如果风力输出将达到其上限，如果风力输出将达到其下限，如果两者都为0，就会实现预测值；
步骤2：考虑不确定性需求侧响应的日前调度建模
步骤2.1：确定性需求侧响应的日前优化调度建模
首先建立确定性需求侧响应的日前调度模型，由于在日前调度优化模型考虑了电价因素，因此目标函数为社会福利最大化，此外在该模型中，风电出力被假设为确定值，价格弹性需求曲线也是确定的，模拟需求曲线和供给曲线，
目标函数描述如下：
<mrow> <mi>max</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> </mrow> </munder> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <msubsup> <mi>r</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <mo>-</mo> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> </mrow> </munder> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <munder> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> </mrow> </munder> <mrow> <mo>(</mo> <mrow> <msubsup> <mi>f</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>SU</mi> <mi>i</mi> <mi>b</mi> </msubsup> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>SD</mi> <mi>i</mi> <mi>b</mi> </msubsup> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> </mrow> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
式中：T表示时间段的集合，B表示节点的集合，Gb表示节点b处的发电机组，表示发电机i在节点b的启动成本，为发电机i在节点b的停机成本，为发电机i在节点b处在时间段t内产生的电量，为节点b在时间段t内的实际用电需求，为节点b在时间段t内的价格弹性需求曲线的积分，fi b为发电机i在节点b处的燃料成本函数；
步骤2.2：线性化机组燃料成本函数
将机组的燃料成本表示为二次函数，使用N段线性函数来近似燃料成本函数：
<mrow> <msubsup> <mi>&amp;phi;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;GreaterEqual;</mo> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mrow> <mi>j</mi> <mi>b</mi> </mrow> </msubsup> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>&amp;beta;</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mrow> <mi>j</mi> <mi>b</mi> </mrow> </msubsup> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>j</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mi>L</mi> <mo>,</mo> <mi>N</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式中:和是第j段函线的截距和斜率，是辅助变量；
步骤2.3：线性化消费者需求响应模型
将消费者剩余中的弹性用电需求响应模型转化为价格弹性需求曲线模型，如果价格弹性是恒定的，可以将价格弹性需求曲线表示为：
<mrow> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>=</mo> <msubsup> <mi>A</mi> <mi>t</mi> <mi>b</mi> </msubsup> <msup> <mrow> <mo>(</mo> <msubsup> <mi>p</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>&amp;alpha;</mi> <mi>t</mi> <mi>b</mi> </msubsup> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式中：是节点b在时间段t内的电价，是t时间段中节点b处所给定的价格弹性值，并且使可以由给定参考点决定的参数值，
对于一般的价格弹性需求曲线，应用分段函数来近似该价格弹性需求曲线，因此将近似为：
<mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>K</mi> </mrow> </munder> <msubsup> <mi>p</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <msubsup> <mi>h</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>&amp;Element;</mo> <mi>K</mi> </mrow> </munder> <msubsup> <mi>h</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mn>0</mn> <mo>&amp;le;</mo> <msubsup> <mi>h</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>l</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>k</mi> <mo>&amp;Element;</mo> <mi>K</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中:是分段函数的第k段，是k段处的相应价格，是在k段处处为需求引入的辅助变量，K是所有的段长集合；
将最大化，是严格地随着k的增加而减少，得到：
<mrow> <msubsup> <mi>h</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>s</mi> </mrow> </msubsup> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <msubsup> <mi>l</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>s</mi> </mrow> </msubsup> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <mi>s</mi> <mo>&lt;</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>s</mi> <mn>0</mn> </msub> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <msubsup> <mi>l</mi> <mi>t</mi> <mrow> <mi>b</mi> <mi>k</mi> </mrow> </msubsup> </mrow> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <mi>s</mi> <mo>=</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>i</mi> <mi>f</mi> <mi> </mi> <mi>s</mi> <mo>&gt;</mo> <msub> <mi>s</mi> <mn>0</mn> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
当存在某个s0使得成立时，可以证明是价格弹性需求曲线的近似积分，即是合理的；
步骤2.4：不确定性需求侧响应的日前调度建模
在步骤2.1中，假设价格弹性需求曲线是确定的，当作出日前调度决策时，必须允许价格弹性需求曲线在一定范围内变化，实际的价格弹性需求曲线是不确定的；为调整保守性，引入参数以限制偏差的总量，即通过改变的值来调整所提出的方法的保守性，值越小，需求响应曲线的不确定性也就越小；
需求响应曲线的不确定性集可以描述如下：
目标函数为：
<mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>&amp;lsqb;</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> </mrow> </munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <msubsup> <mi>r</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>(</mo> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> <mo>-</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> </mrow> </munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> </mrow> </munder> <mrow> <mo>(</mo> <msubsup> <mi>f</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>(</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>)</mo> <mo>+</mo> <msubsup> <mi>SU</mi> <mi>i</mi> <mi>b</mi> </msubsup> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>SD</mi> <mi>i</mi> <mi>b</mi> </msubsup> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
约束条件如下：
<mrow> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <mn>0</mn> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>k</mi> <mo>:</mo> <mn>1</mn> <mo>&amp;le;</mo> <mi>k</mi> <mo>-</mo> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <msubsup> <mi>MU</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>k</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <mn>1</mn> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>k</mi> <mo>:</mo> <mn>1</mn> <mo>&amp;le;</mo> <mi>k</mi> <mo>-</mo> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <msubsup> <mi>MD</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <mn>0</mn> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <mn>0</mn> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>MI</mi> <mi>i</mi> <mi>b</mi> </msubsup> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>MA</mi> <mi>i</mi> <mi>b</mi> </msubsup> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <msubsup> <mi>RU</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>MA</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <msubsup> <mi>RD</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <msubsup> <mi>MA</mi> <mi>i</mi> <mi>b</mi> </msubsup> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <mrow> <mo>(</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> </mrow> </munder> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>w</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mo>-</mo> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>&amp;le;</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> </mrow> </munder> <msubsup> <mi>F</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>b</mi> </msubsup> <mrow> <mo>(</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>q</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> </mrow> </munder> <msubsup> <mi>x</mi> <mrow> <mi>q</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>+</mo> <msubsup> <mi>w</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>-</mo> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>)</mo> </mrow> <mo>&amp;le;</mo> <msub> <mi>U</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>,</mo> <mo>&amp;ForAll;</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo> </mrow> <mo>&amp;Element;</mo> <mi>&amp;Omega;</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>D</mi> <mrow> <mi>t</mi> <mi>b</mi> </mrow> <mn>0</mn> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>&amp;le;</mo> <msubsup> <mi>D</mi> <mrow> <mi>t</mi> <mi>b</mi> </mrow> <mi>M</mi> </msubsup> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>T</mi> <mo>&amp;Element;</mo> <mi>t</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msubsup> <mi>y</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>,</mo> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>,</mo> <msubsup> <mi>v</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>&amp;Element;</mo> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>}</mo> <mo>;</mo> <msubsup> <mi>x</mi> <mrow> <mi>i</mi> <mi>t</mi> </mrow> <mi>b</mi> </msubsup> <mo>,</mo> <msubsup> <mi>d</mi> <mi>t</mi> <mi>b</mi> </msubsup> <mo>&amp;GreaterEqual;</mo> <mn>0</mn> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>t</mi> <mo>&amp;Element;</mo> <mi>T</mi> <mo>;</mo> <mo>&amp;ForAll;</mo> <mi>b</mi> <mo>&amp;Element;</mo> <mi>B</mi> <mo>,</mo> <mo>&amp;ForAll;</mo> <mi>i</mi> <mo>&amp;Element;</mo> <msub> <mi>G</mi> <mi>b</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
步骤3：模型求解
步骤3.1：按照机会均等原则，从总的决策空间Θ里挑选出N个决策，并构成有限集合ΘN；
步骤3.2：利用某种选择规则，通过快速粗糙评估，从N个决策中确定OPC类型。利用选择规则选出一子决策集S，其元素的个数为s，子决策集S将包含k个足够好的解(以至少α的概率)，相关的数学公式如下所示：
P{|G∩S|>k}≥α (25)
式中：G表示足够好的解的集合，解集中元素的个数为g，表示将各个决策经精确评估后的前g个最好的解；||表示集合元素的个数；
步骤3.3：精确评估子决策集中的每一个元素，从s个元素中选择出一个最好的解。
CN201711218112.4A 2017-11-28 2017-11-28 一种考虑需求侧响应不确定性的电力系统日前优化调度方法 Active CN107910866B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711218112.4A CN107910866B (zh) 2017-11-28 2017-11-28 一种考虑需求侧响应不确定性的电力系统日前优化调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711218112.4A CN107910866B (zh) 2017-11-28 2017-11-28 一种考虑需求侧响应不确定性的电力系统日前优化调度方法

Publications (2)

Publication Number Publication Date
CN107910866A true CN107910866A (zh) 2018-04-13
CN107910866B CN107910866B (zh) 2021-09-03

Family

ID=61848987

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711218112.4A Active CN107910866B (zh) 2017-11-28 2017-11-28 一种考虑需求侧响应不确定性的电力系统日前优化调度方法

Country Status (1)

Country Link
CN (1) CN107910866B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108667012A (zh) * 2018-05-21 2018-10-16 国网山东省电力公司电力科学研究院 基于多场景的区域能源互联网源荷储双阶段协调优化方法
CN114662756A (zh) * 2022-03-23 2022-06-24 福州大学 考虑源-荷不确定性的热电微能网双层仿射优化调度方法
CN117852737A (zh) * 2024-01-17 2024-04-09 中国电子科技集团公司第十五研究所 一种用于精确制导弹药战时实时调配的方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014078336A1 (en) * 2012-11-14 2014-05-22 Autogrid, Inc. Identifying operability failure in dr assets
CN106779180A (zh) * 2016-11-29 2017-05-31 国网陕西省电力公司电力科学研究院 基于曲线分段线性化的电力系统旋转备用优化方法
CN107017630A (zh) * 2017-05-18 2017-08-04 南京师范大学 一种多时间尺度的电力需求侧响应方法
CN107276096A (zh) * 2017-06-28 2017-10-20 国网江苏省电力公司电力科学研究院 一种光伏与空调负荷协调优化的配电网电压控制方法
CN107330546A (zh) * 2017-06-14 2017-11-07 武汉大学 一种考虑风电出力和需求响应不确定性的优化调度方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2014078336A1 (en) * 2012-11-14 2014-05-22 Autogrid, Inc. Identifying operability failure in dr assets
CN106779180A (zh) * 2016-11-29 2017-05-31 国网陕西省电力公司电力科学研究院 基于曲线分段线性化的电力系统旋转备用优化方法
CN107017630A (zh) * 2017-05-18 2017-08-04 南京师范大学 一种多时间尺度的电力需求侧响应方法
CN107330546A (zh) * 2017-06-14 2017-11-07 武汉大学 一种考虑风电出力和需求响应不确定性的优化调度方法
CN107276096A (zh) * 2017-06-28 2017-10-20 国网江苏省电力公司电力科学研究院 一种光伏与空调负荷协调优化的配电网电压控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
罗纯坚等: "需求响应不确定性对日前优化调度的影响分析", 《电力系统自动化》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108667012A (zh) * 2018-05-21 2018-10-16 国网山东省电力公司电力科学研究院 基于多场景的区域能源互联网源荷储双阶段协调优化方法
CN114662756A (zh) * 2022-03-23 2022-06-24 福州大学 考虑源-荷不确定性的热电微能网双层仿射优化调度方法
CN117852737A (zh) * 2024-01-17 2024-04-09 中国电子科技集团公司第十五研究所 一种用于精确制导弹药战时实时调配的方法及装置

Also Published As

Publication number Publication date
CN107910866B (zh) 2021-09-03

Similar Documents

Publication Publication Date Title
Luna et al. Mixed-integer-linear-programming-based energy management system for hybrid PV-wind-battery microgrids: Modeling, design, and experimental verification
Karami et al. An optimal dispatch algorithm for managing residential distributed energy resources
Soliman et al. Supervisory energy management of a hybrid battery/PV/tidal/wind sources integrated in DC-microgrid energy storage system
Ma et al. Optimal allocation of hybrid energy storage systems for smoothing photovoltaic power fluctuations considering the active power curtailment of photovoltaic
Wang et al. Adaptive control for energy storage systems in households with photovoltaic modules
Hong et al. Multiscenario underfrequency load shedding in a microgrid consisting of intermittent renewables
Pourmousavi et al. Real-time energy management of a stand-alone hybrid wind-microturbine energy system using particle swarm optimization
Chen et al. Jump and shift method for multi-objective optimization
Colson et al. Improving sustainability of hybrid energy systems part II: Managing multiple objectives with a multiagent system
Mengi et al. A new energy management technique for PV/wind/grid renewable energy system
KR20210100699A (ko) 하이브리드 발전소
Cheng et al. Adaptive robust method for dynamic economic emission dispatch incorporating renewable energy and energy storage
CN107910866B (zh) 一种考虑需求侧响应不确定性的电力系统日前优化调度方法
CN110350512A (zh) 一种智能园区新能源发电站调度优化方法及系统
Tang et al. Optimal operation method for microgrid with wind/PV/diesel generator/battery and desalination
CN108288132B (zh) 一种基于源荷互动电力系统调度的建模方法
Shahzad et al. Model predictive control strategies in microgrids: A concise revisit
Li et al. Flexible scheduling of microgrid with uncertainties considering expectation and robustness
Belkhier et al. Novel design and adaptive coordinated energy management of hybrid fuel‐cells/tidal/wind/PV array energy systems with battery storage for microgrids
Sarasua et al. Hybrid frequency control strategies based on hydro‐power, wind, and energy storage systems: Application to 100% renewable scenarios
Ma et al. Optimal economic operation of islanded microgrid by using a modified pso algorithm
Prakash et al. Intelligent Energy Management for Distributed Power Plants and Battery Storage
CN109711605A (zh) 一种多类型电源并网定容方法和装置
Li et al. Optimization scheduling model based on source-load-energy storage coordination in power systems
Sechilariu et al. A conceptual framework for full optimal operation of a grid-connected DC microgrid

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