CN111241459B - 一种考虑自趋优的工业园区分布式综合需求响应方法 - Google Patents
一种考虑自趋优的工业园区分布式综合需求响应方法 Download PDFInfo
- Publication number
- CN111241459B CN111241459B CN202010014848.5A CN202010014848A CN111241459B CN 111241459 B CN111241459 B CN 111241459B CN 202010014848 A CN202010014848 A CN 202010014848A CN 111241459 B CN111241459 B CN 111241459B
- Authority
- CN
- China
- Prior art keywords
- idr
- industrial
- user
- period
- cost
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000004044 response Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000005611 electricity Effects 0.000 claims abstract description 31
- 238000005265 energy consumption Methods 0.000 claims abstract description 18
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000011156 evaluation Methods 0.000 claims abstract description 6
- 239000007789 gas Substances 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 13
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 12
- 238000012423 maintenance Methods 0.000 claims description 7
- 230000008901 benefit Effects 0.000 claims description 6
- 238000004146 energy storage Methods 0.000 claims description 6
- 239000003345 natural gas Substances 0.000 claims description 6
- 238000010248 power generation Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 abstract description 12
- 230000008569 process Effects 0.000 abstract description 5
- 230000003993 interaction Effects 0.000 abstract description 4
- 238000009826 distribution Methods 0.000 abstract description 3
- 238000011160 research Methods 0.000 description 5
- 238000007726 management method Methods 0.000 description 4
- 238000005057 refrigeration Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 101150067055 minC gene Proteins 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 239000002918 waste heat Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0206—Price or cost determination based on market factors
-
- 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
-
- 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
-
- 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
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards
-
- 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
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/16—Energy services, e.g. dispersed generation or demand or load or energy savings aggregation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Accounting & Taxation (AREA)
- Marketing (AREA)
- Data Mining & Analysis (AREA)
- Finance (AREA)
- Mathematical Physics (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Tourism & Hospitality (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Water Supply & Treatment (AREA)
- Pure & Applied Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Public Health (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种考虑自趋优的工业园区分布式综合需求响应方法。工业用户首先对分时电价做出响应，并上报用能计划，若系统总体用能情况不满足安全性评估，则向园区下发调峰需求指令。在设计的IDR机制下，上层园区以IDR总成本最小化为目标，下层工业用户以自身参与IDR后的成本最小为目标，进行双层分布式优化。最终用户对IDR的边际成本达成一致，得到使用户自身和园区IDR总成本均最小的削峰量分配方案，并且在参与IDR的过程中较好地保护了用户的信息私密性本发明规划方法对问题的考虑更为全面合理，解决了现有规划方法中上级电网对响应结果安全性的评估及其与用户的互动不足，对用户的信息私密性保护不够的问题。
Description
技术领域
本发明涉及工业园区综合需求响应领域，尤其涉及一种考虑自趋优的工业园 区分布式综合需求响应方法。
背景技术
随着化石能源的快速枯竭及其对环境、气候的影响日益显著，提高能源利用 效率、充分利用现有的可调控资源而非单纯增加装机容量成为可持续发展的必由 之路。需求响应(demand response，DR)通过整合需求侧的可调控资源与电网进行“双向互动”，是应对日益增长且多样化的电力需求和缓解电力供需不平衡的有效 手段。目前国内外关于区域综合能源系统规划的研究主要集中在系统设备容量规 划和考虑多能耦合的配电网扩展规划，而针对区域综合能源系统中能源站的数量、 位置及供能管线布局规划的研究较少。
目前国内外关于综合需求响应(integrated demand response,IDR)的研究已经获得了一定成果。戚野白等人(戚野白，王丹，贾宏杰，等.基于需求响应的集群温控负荷建模及参与低频减载策略研究[J].中国电机工程学报，2017，37(3)： 751-759.)建立了电热水器集群温控负荷的模型，研究了其在频率响应下的控制策略。Tsui K M等人(Tsui K M，Chan S C.Demand response optimization for smart home scheduling under real-time pricing[J].IEEE Transactions on Smart Grid，2012， 3(4)：1812-1821.)对常见居民用户负荷的需求响应特性进行分类，提出利用家 庭能量管理系统(energymanagement system，EMS)响应实时电价的方法。曾丹等 人(曾丹，姚建国，杨胜春，等.应对风电消纳中基于安全约束的价格型需求响 应优化调度建模[J].中国电机工程学报，2014，34(31)：5571-5578.)研究了具 有固定价格弹性系数的可调控负荷参与平衡风电功率波动的调度过程。上述文献 分析了具有特定价格响应特性和运行约束的负荷参与DR的策略，未考虑不同能源间的互补作用。
综上所述，现有的IDR研究多为集中式的统一优化调度，上级能量管理系 统(energy management system，EMS)需要获取用户的生产、用能等可能涉及私密性的信息，降低了用户参与需求响应的意愿。此外，现有研究中通常考虑用户对 分时电价做出响应，但未考虑上级电网对响应结果安全性的评估及其与用户的互 动。对于工业园区内部用户参与IDR的容量分配和补偿机制，目前也尚未有充 分的研究。
发明内容
针对上述问题，本发明考虑工业园区与上级电网的双向互动，提出一种园区 用户分布式地参与IDR的机制。工业用户首先对分时电价做出响应，考虑生产 过程中的多能协同和热能梯级利用，以用能成本最小化为目标进行自趋优，并上 报用能计划。上级电网对系统总体用能情况进行安全性评估，若不能通过安全校 核，则向园区下发调峰需求指令。进而在设计的IDR机制下，上层园区以IDR 总成本最小化为目标，下层工业用户以自身参与IDR后的成本最小为目标，进行双层分布式优化。最终用户对IDR的边际成本达成一致，得到使用户自身和 园区IDR总成本均最小的削峰量分配方案，并且在参与IDR的过程中较好地保 护了用户的信息私密性。
本发明的目的是通过以下技术方案来实现的：
一种考虑自趋优的工业园区分布式综合需求响应方法，包括如下步骤：
步骤1：在工业园区内部建立以统一价格结算的园区IDR市场；工业用户对 分时电价做出响应，以用能成本最小化为目标进行自趋优，并上报用能计划；上 级电网对工业园区总体用能情况进行安全性评估，若能通过安全校核，则工业用 户按原计划用能，若不能通过安全校核，则上级电网向工业园区下发削峰需求指 令和电网调峰辅助服务价格；
步骤2：工业园区EMS根据上级电网的削峰需求指令和工业用户缴纳的基 本电费，分配各个工业用户应削峰量；
步骤3：各个工业用户根据最新的园区IDR市场出清价格进行自趋优，并申 报期望IDR容量和期望价格，其中初始出清价格为电网调峰辅助服务价格；
步骤4：工业园区EMS计算所有工业用户申报的期望IDR容量的总量，判 断IDR总量和电网下发的削峰需求总量的大小，若IDR总量≤削峰需求总量， 则用户按照所申报的期望IDR容量参与需求响应，园区IDR市场按照上一次更 新的出清价格进行结算；
若IDR总量＞削峰需求总量，工业园区EMS以IDR市场供应方的总成本最 小为目标选择中标用户，并得到更新后的出清价格；重复步骤3～4，直到IDR总 量≤削峰需求总量。
进一步的，步骤1所述的统一价格不应高于电网调峰辅助服务价格。
进一步的，所述的步骤2具体为：
工业园区EMS参考工业用户缴纳的基本电费对应的需量分配各个工业用户 的应削峰量，首先计算各工业用户的可用容量基准值：
其中，t为一天内的时段编号，Pref,i,t为第i个工业用户在时段t的可用容量 基准值，Pi,t为未考虑IDR时工业用户在时段t的计划购电功率，Pd,t为上级电网 在时段t的削峰需求总量，N为工业园区内的工业用户数量，Pb,i为第i个工业用 户的需量；
计算各工业用户的应削峰量：
Pdref,i,t＝Pi,t-Pref,i,t
其中，Pdref,i,t为第i个工业用户在时段t的应削峰量。
进一步的，所述的步骤3具体为：
步骤3.1：将一天分为96个时段，工业用户自趋优的目的是使其参与IDR 后的总成本最小，目标函数如下：
min Ci＝Ceqp,i-Ii
Ceqp,i＝Ce,i+Cgas,i+Ch,i+Com,i+Css,i+Cbsw,i
式中：Ci为第i个工业用户参与IDR后的总成本，Ii为第i个工业用户在园 区IDR市场中获得的收益，pIDR为园区IDR市场出清价格，初始出清价格为电 网调峰辅助服务价格；PDR,i,t为第i个工业用户参与IDR后在时段t的购电需求， Γ为削峰时段集合，Ceqp,i为第i个工业用户的日用能成本，Ce,i、Cgas,i和Ch,i分别 为第i个工业用户的购电、购气和购热成本，Com,i为第i个工业用户的设备运行 维护成本，Css,i为第i个工业用户的燃气轮机和燃气锅炉的启停成本，Cbsw,i为第 i个工业用户的电池储能折旧成本；
为了避免工业用户在需要调峰的一部分时段内用尽所有的可调控资源，导致 其余时段无法提供需求响应，要求工业用户在各调峰时段等比例地进行削峰；
其中，αi为一正数变量；
购电成本的计算公式如下：
其中，ce,t为时段t的电价，T为单位时段的时长；
购气成本的计算公式如下：
其中，cgas为天然气价格，FGT,i,j,t和FGB,i,j,t分别为第i个工业用户的第j台燃气 轮机在时段t的产电功率和天然气消耗速率；
购热成本的计算公式如下：
其中，ch为CHP热蒸汽价格，Hbuy,i,t为第i个工业用户在时段t的购热功率；
设备运行维护成本的计算公式如下：
其中，com,i,j为第i个工业用户的第j台设备单位输出功率的运行维护费用，Pout,i,j,t为第i个工业用户的第j台设备在时段t的输出功率；
启停成本的计算公式如下：
其中，css,i,j为第i个工业用户的第j台设备的开停机费用，ui,j,t为第i个工业 用户的第j台设备在时段t的状态，ui,j,t＝1表示开机，ui,j,t＝0表示停机；
电池储能折旧成本的计算公式如下：
其中，cbsw为蓄电池累计充电1kWh的折旧成本，Pbs,c,i,j,t为第i个工业用户 的第j个电池储能在时段t的充电功率。
步骤3.2：根据步骤3.1得到的min Ci获得PDR,i,t，计算各个工业用户的期望IDR容量为：Pi,t-PDR,i,t；期望价格pi为：
进一步的，所述的步骤4具体为：
工业园区EMS计算所有工业用户申报的期望IDR容量的总量，判断IDR总 量和电网下发的削峰需求总量的大小；
若IDR总量≤削峰需求总量，则用户按照所申报的期望IDR容量参与需求 响应，园区IDR市场按照上一次更新的出清价格进行结算，若出清价格还未进 行迭代更新，则按照电网调峰辅助服务价格进行结算；
若IDR总量＞削峰需求总量，工业园区EMS以IDR市场供应方的总成本最 小为目标选择中标用户，确定用户在时段t的削峰比例αΔ,i,t，其目标函数为：
0≤αΔ,i,t≤αi
其中，C为IDR市场供应方的总成本，Λ为供应方的编号集合，αΔ,i,t为第i 个工业用户在时段t的削峰比例，若αΔ,i,t＞0，则说明该工业用户中标；根据中 标的工业用户在时段t的削峰比例αΔ,i,t，得到各工业用户的中标容量：αΔ,i,t*Pdref,i,t；
将所有中标的工业用户中申报的期望价格最大值作为更新后的出清价格。
本发明具备的有益效果是：
1)与集中式IDR相比，用户提供给园区EMS的信息较少，有利于保护用 户的信息私密性和控制自主权，使得需求响应更容易实际实施；
2)借鉴并结合现有市场机制设计了工业用户参与调峰需求响应的容量与收 益分配机制，从而在现有的市场机制下具备实施的可行性；同时在该机制下用户 可获得比各自独立参与IDR更高的收益，提高了用户参与IDR的积极性；
3)在IDR过程中用户能够以自身成本最小化为目标进行自趋优，可促进用 户充分考虑冷、热、电、气等多种能源在时空上的互补特性，挖掘自身的IDR 潜力；同时本发明分布式需求响应的结果可使园区整体IDR成本也达到最小。
附图说明
图1为工业用户响应调峰需求的流程图；
图2为本实施例中工业用户计划购电功率；
图3为本实施例中用户1设备优化运行结果。
图4-a为各用户的交流电负荷曲线图；
图4-b为各用户的直流电负荷曲线图；
图4-c为各用户的热负荷曲线图；
图4-d为各用户的冷负荷曲线图；
图4-e为各用户的蒸汽驱动负荷曲线图；
图4-f为各用户的光伏预测出力曲线图。
具体实施方式
一种考虑能源网拓扑特性的区域综合能源系统布局规划方法，该方法包括以 下步骤：
步骤一、在工业园区内部建立以统一价格结算的园区IDR市场，所述的统 一价格不应高于电网调峰辅助服务价格。工业用户对分时电价做出响应，以用能 成本最小化为目标进行自趋优，并上报用能计划；上级电网对工业园区总体用能 情况进行安全性评估，若能通过安全校核，则工业用户按原计划用能，若不能通 过安全校核，则上级电网向工业园区下发削峰需求指令和电网调峰辅助服务价格；
步骤二：工业园区EMS根据上级电网的削峰需求指令和工业用户缴纳的基 本电费，分配各个工业用户应削峰量。
首先计算各工业用户的可用容量基准值：
其中，t为一天内的时段编号，Pref,i,t为第i个工业用户在时段t的可用容量基准值，Pi,t为未考虑IDR时工业用户在时段t的计划购电功率，Pd,t为上级电网在时 段t的削峰需求总量，N为工业园区内的工业用户数量，Pb,i为第i个工业用户的 需量；
计算各工业用户的应削峰量：
Pdref,i,t＝Pi,t-Pref,i,t
其中，Pdref,i,t为第i个工业用户在时段t的应削峰量。
步骤三：各个工业用户根据最新的园区IDR市场出清价格进行自趋优，并 申报期望IDR容量和期望价格，其中初始出清价格为电网调峰辅助服务价格；
将一天分为96个时段，工业用户自趋优的目的是使其参与IDR后的总成本 最小，目标函数如下：
min Ci＝Ceqp,i-Ii
Ceqp,i＝Ce,i+Cgas,i+Ch,i+Com,i+Css,i+Cbsw,i
其中：Ci为第i个工业用户参与IDR后的总成本，Ii为第i个工业用户在园区IDR 市场中获得的收益，pIDR为园区IDR市场出清价格，初始出清价格为电网调峰辅 助服务价格；PDR,i,t为第i个工业用户参与IDR后在时段t的购电需求，Γ为削峰 时段集合，Ceqp,i为第i个工业用户的日用能成本，Ce,i、Cgas,i和Ch,i分别为第i个 工业用户的购电、购气和购热成本，Com,i为第i个工业用户的设备运行维护成本， Css,i为第i个工业用户的燃气轮机和燃气锅炉的启停成本，Cbsw,i为第i个工业用 户的电池储能折旧成本；
为了避免工业用户在需要调峰的一部分时段内用尽所有的可调控资源，导致 其余时段无法提供需求响应，要求工业用户在各调峰时段等比例地进行削峰；
其中，αi为一正数变量；
根据得到的minCi获得PDR,it,，计算各个工业用户的期望IDR容量为： Pi,t-PDR,i,t；期望价格pi为：
步骤四：工业园区EMS计算所有工业用户申报的期望IDR容量的总量，判 断IDR总量和电网下发的削峰需求总量的大小；
若IDR总量≤削峰需求总量，则用户按照所申报的期望IDR容量参与需求 响应，园区IDR市场按照上一次更新的出清价格进行结算，若出清价格还未进 行迭代更新，则按照电网调峰辅助服务价格进行结算；
若IDR总量＞削峰需求总量，工业园区EMS以IDR市场供应方的总成本最 小为目标选择中标用户，确定用户在时段t的削峰比例αΔ,i,t，其目标函数为：
0≤αΔ,i,t≤αi
其中，C为IDR市场供应方的总成本，Λ为供应方的编号集合，αΔ,i,t为第i 个工业用户在时段t的削峰比例，若αΔ,i,t＞0，则说明该工业用户中标；根据中 标的工业用户在时段t的削峰比例αΔ,i,t，得到各工业用户的中标容量：αΔ,i,t*Pdref,i,t；
将所有中标的工业用户中申报的期望价格最大值作为更新后的出清价格。
重复步骤三至步骤四，直到IDR总量≤削峰需求总量。
实施例
以某工业园区为例分析所提的分布式IDR方法。在工业园区中有5个用户 参与IDR，各用户工厂综合能源系统中的设备数量如表1所示，各类设备的参数 如表2所示。广州于2019年4月起执行的分时电价如表3所示。
表1各用户的设备
表2设备参数
表3分时电价
各用户的电、热、冷、蒸汽驱动负荷和光伏预测出力如图4-a至4-f所示。 天然气和CHP热蒸汽的价格分别为3.45元/m3和348元/t，折合成单位热值价格 分别为0.349元/kWh和0.465元/kWh。5个用户缴纳的基本电费所对应的需量分 别为：4.5MW、2.1MW、1.6MW、1.8MW和3.8MW。各用户的可中断负荷 参数如表4所示。相邻两次市场出清价格变化上限κ＝0.1元，最小价格单位 ε＝0.001元。
表4各用户可中断负荷参数
1.分时电价响应结果：
根据工业园区IDR机制，工业用户首先对分时电价做出响应，向园区EMS 上报用能计划。各用户利用自趋优模型优化得到的计划购电功率如图2所示。由 图2可以得出，用户对分时电价的响应结果是尽可能在电价谷时段、平时段和峰 时段的购电功率依次减小，从而降低用能成本。同时结合表1，可以得出能源设 备种类、数量较多，多能耦合更为紧密的用户具有更强的IDR能力。以用户1 为例，分析其对分时电价的具体响应情况，该用户经自趋优后的设备运行计划如 图3所示。图3所示的设备优化运行结果具有如下特征：电池储能在电价谷时段 充电，在峰时段放电，在平时段放电或为下一个峰时段补充电量；在峰时段和平时段，燃气轮机的发电成本低于向电网购电成本，处于满发状态；冰蓄冷装置在 谷时段进行蓄冰，在峰时段融冰释放冷量；吸收式制冷机可利用工业生产过程产 生的余热作为热源进行制冷，是工厂中最经济的制冷设备，用户优先使用吸收式制冷机制冷，不足的冷量由冰蓄冷装置和中央空调补充。
根据上述分析，本方法的工业用户自趋优模型可有效利用工厂中的能源设备 实现多能协同互补，充分挖掘用户对分时电价的响应潜力，同时可对热能进行梯 级利用，进一步降低用能成本。
2.调峰需求响应结果：
得到工业用户自趋优的购电计划后，上级电网通过评估系统整体用能情况的 安全性，向园区下发调峰需求：在11:15—12:15的四个时段，调峰需求分别为 4.50MW、3.79MW、2.86MW和2.96MW。同时下发电网的调峰辅助服务价格 为0.540元/kWh。
园区EMS收到电网调峰需求指令后，各用户的应削峰量计算结果如表5所 示。
表5各用户的应削峰量
初始以电网调峰辅助服务价格作为工业园区IDR市场的结算价格。用户利 用考虑IDR的自趋优模型计算最优的削峰量与报价，园区EMS以IDR总成本最 小化为目标确定市场出清结果，进行图1所示的分布式IDR响应过程。
经过6次迭代后计算收敛，得到最终的IDR市场结算价格为0.140元/kWh。 用户1-5申报的削峰比例分别为3.20、0.897、0、0、0，其中用户1在IDR市场中为供应方，用户2、3、4、5为需求方。各个用户最终提供的削峰容量如表6 所示。
表6各用户提供的IDR容量
由上述市场出清结果可知，用户1的IDR能力最强，在相同的结算价格下， 意愿提供的IDR容量最大。用户3、4、5参与IDR成本较高，在最终的市场结 算价格下，自身不提供削峰容量，而是选择从IDR市场购买削峰容量。用户2 提供的削峰容量小于其应削峰量，不足的部分在IDR市场中购买。
5个用户参与IDR后的用能成本分别为95356、73835、66093、45019和59199 元，关于IDR容量边际成本分别为0.131、0.140、0.230、0.200和0.235元/kWh。 由此可看出用户1的边际成本小于结算价格，其作为供应方可从IDR市场获得 收益。用户2的边际成本等于结算价格，此时其削峰成本比从市场购买全部应削 峰量或由自身提供全部应削峰量都更低。用户3、4、5的边际成本此时已经达到 最低值，但仍大于结算价格，因此选择从IDR市场购买削峰容量相比其自身进行削峰的效益更高。
最终各用户在IDR市场获得的收益分别为261、-14、-60、-66和-120元。 上级电网支付给工业园区的IDR补偿为1512元，园区EMS按照各用户可用容 量基准值与其基本电费对应需量的偏差分配从电网的补偿收益，分别为493、230、 175、197和416元。各用户响应分时电价后的初始计划用能成本、参与分布式 IDR后的总成本和各自独立进行需求响应的成本如表7所示。
表7各用户的成本
由表7可知，相比于初始计划用能的成本，各用户按本方法IDR机制参与 削峰需求响应后的总成本更低，提高了用户具有参与IDR的意愿。此外，若园 区用户各自独立地参与IDR，则无法实现用户之间IDR资源的协调，每个用户的成本相比采用本方法机制均有所提高。
通过上述分析可以得出，本方法的工业园区IDR机制能够有效地配置园区 内IDR资源。在削峰需求量一定的约束下，可使得各用户参与IDR的成本达到 最低或使收益达到最大。同时，最终迭代结果使园区整体的IDR成本也达到最 小。
以上列举的仅是本发明的具体实施例。显然，本发明不限于以上实施例，还 可以有许多变形。本领域的普通技术人员能从本发明公开的内容直接导出或联想 到的所有变形，均应认为是本发明的保护范围。
Claims (4)
1.一种考虑自趋优的工业园区分布式综合需求响应方法，其特征在于，包括如下步骤：
步骤1：在工业园区内部建立以统一价格结算的园区IDR市场；工业用户对分时电价做出响应，以用能成本最小化为目标进行自趋优，并上报用能计划；上级电网对工业园区总体用能情况进行安全性评估，若能通过安全校核，则工业用户按原计划用能，若不能通过安全校核，则上级电网向工业园区下发削峰需求指令和电网调峰辅助服务价格；
步骤2：工业园区EMS根据上级电网的削峰需求指令和工业用户缴纳的基本电费，分配各个工业用户应削峰量；
步骤3：各个工业用户根据最新的园区IDR市场出清价格进行自趋优，并申报期望IDR容量和期望价格，其中初始出清价格为电网调峰辅助服务价格；
所述的步骤3具体为：
步骤3.1：将一天分为96个时段，工业用户自趋优的目的是使其参与IDR后的总成本最小，目标函数如下：
min Ci＝Ceqp,i-Ii
Ceqp,i＝Ce,i+Cgas,i+Ch,i+Com,i+Css,i+Cbsw,i
式中：Pref,i,t为第i个工业用户在时段t的可用容量基准值，Ci为第i个工业用户参与IDR后的总成本，Ii为第i个工业用户在园区IDR市场中获得的收益，pIDR为园区IDR市场出清价格，初始出清价格为电网调峰辅助服务价格；PDR,i,t为第i个工业用户参与IDR后在时段t的购电需求，Γ为削峰时段集合，Ceqp,i为第i个工业用户的日用能成本，Ce,i、Cgas,i和Ch,i分别为第i个工业用户的购电、购气和购热成本，Com,i为第i个工业用户的设备运行维护成本，Css,i为第i个工业用户的燃气轮机和燃气锅炉的启停成本，Cbsw,i为第i个工业用户的电池储能折旧成本；
为了避免工业用户在需要调峰的一部分时段内用尽所有的可调控资源，导致其余时段无法提供需求响应，要求工业用户在各调峰时段等比例地进行削峰；
其中，t为一天内的时段编号，Pi,t为未考虑IDR时工业用户在时段t的计划购电功率，αi为一正数变量；
购电成本的计算公式如下：
其中，ce,t为时段t的电价，T为单位时段的时长；
购气成本的计算公式如下：
其中，cgas为天然气价格，FGT,i,j,t和FGB,i,j,t分别为第i个工业用户的第j台燃气轮机在时段t的产电功率和天然气消耗速率；
购热成本的计算公式如下：
其中，ch为CHP热蒸汽价格，Hbuy,i,t为第i个工业用户在时段t的购热功率；
设备运行维护成本的计算公式如下：
其中，com,i,j为第i个工业用户的第j台设备单位输出功率的运行维护费用，Pout,i,j,t为第i个工业用户的第j台设备在时段t的输出功率；
启停成本的计算公式如下：
其中，css,i,j为第i个工业用户的第j台设备的开停机费用，ui,j,t为第i个工业用户的第j台设备在时段t的状态，ui,j,t＝1表示开机，ui,j,t＝0表示停机；
电池储能折旧成本的计算公式如下：
其中，cbsw为蓄电池累计充电1kWh的折旧成本，Pbs,c,i,j,t为第i个工业用户的第j个电池储能在时段t的充电功率；
步骤3.2：根据步骤3.1得到的min Ci获得PDR,i,t，计算各个工业用户的期望IDR容量为：Pi,t-PDR,i,t；期望价格pi为：
步骤4：工业园区EMS计算所有工业用户申报的期望IDR容量的总量，判断IDR总量和电网下发的削峰需求总量的大小，若IDR总量≤削峰需求总量，则用户按照所申报的期望IDR容量参与需求响应，园区IDR市场按照上一次更新的出清价格进行结算；
若IDR总量＞削峰需求总量，工业园区EMS以IDR市场供应方的总成本最小为目标选择中标用户，并得到更新后的出清价格；重复步骤3～4，直到IDR总量≤削峰需求总量。
2.如权利要求1所述的一种考虑自趋优的工业园区分布式综合需求响应方法，其特征在于，步骤1所述的统一价格不应高于电网调峰辅助服务价格。
3.如权利要求1所述的一种考虑自趋优的工业园区分布式综合需求响应方法，其特征在于，所述的步骤2具体为：
工业园区EMS参考工业用户缴纳的基本电费对应的需量分配各个工业用户的应削峰量，首先计算各工业用户的可用容量基准值：
其中，t为一天内的时段编号，Pref,i,t为第i个工业用户在时段t的可用容量基准值，Pi,t为未考虑IDR时工业用户在时段t的计划购电功率，Pd,t为上级电网在时段t的削峰需求总量，N为工业园区内的工业用户数量，Pb,i为第i个工业用户的需量；
计算各工业用户的应削峰量：
Pdref,i,t＝Pi,t-Pref,i,t
其中，Pdref,i,t为第i个工业用户在时段t的应削峰量。
4.如权利要求1所述的一种考虑自趋优的工业园区分布式综合需求响应方法，其特征在于，所述的步骤4具体为：
工业园区EMS计算所有工业用户申报的期望IDR容量的总量，判断IDR总量和电网下发的削峰需求总量的大小；
若IDR总量≤削峰需求总量，则用户按照所申报的期望IDR容量参与需求响应，园区IDR市场按照上一次更新的出清价格进行结算，若出清价格还未进行迭代更新，则按照电网调峰辅助服务价格进行结算；
若IDR总量＞削峰需求总量，工业园区EMS以IDR市场供应方的总成本最小为目标选择中标用户，确定用户在时段t的削峰比例αΔ,i,t，其目标函数为：
0≤αΔ,i,t≤αi
其中，Pdref,i,t为第i个工业用户在时段t的应削峰量，Pd,t为上级电网在时段t的削峰需求总量，C为IDR市场供应方的总成本，Λ为供应方的编号集合，αΔ,i,t为第i个工业用户在时段t的削峰比例，若αΔ,i,t＞0，则说明该工业用户中标；根据中标的工业用户在时段t的削峰比例αΔ,i,t，得到各工业用户的中标容量：αΔ,i,t*Pdref,i,t；
将所有中标的工业用户中申报的期望价格最大值作为更新后的出清价格。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010014848.5A CN111241459B (zh) | 2020-01-07 | 2020-01-07 | 一种考虑自趋优的工业园区分布式综合需求响应方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010014848.5A CN111241459B (zh) | 2020-01-07 | 2020-01-07 | 一种考虑自趋优的工业园区分布式综合需求响应方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111241459A CN111241459A (zh) | 2020-06-05 |
| CN111241459B true CN111241459B (zh) | 2023-07-21 |
Family
ID=70872468
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010014848.5A Active CN111241459B (zh) | 2020-01-07 | 2020-01-07 | 一种考虑自趋优的工业园区分布式综合需求响应方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111241459B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114330843B (zh) * | 2021-12-15 | 2024-04-05 | 国网江苏省电力有限公司常州供电分公司 | 面向虚拟电厂的需求响应调度方法和系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107292416A (zh) * | 2017-05-03 | 2017-10-24 | 国网浙江省电力公司 | 基于回归预测需量的电力用户基本电费决策方法 |
| CN107832979A (zh) * | 2017-12-06 | 2018-03-23 | 浙江大学 | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 |
| CN108808663A (zh) * | 2018-06-12 | 2018-11-13 | 浙江大学 | 一种基于多能互补的工业用户热需求响应方法 |
| CN108898290A (zh) * | 2018-06-13 | 2018-11-27 | 广州供电局有限公司 | 一种面向多参与主体的工业园区综合能源系统互动方法 |
| CN110390443A (zh) * | 2018-04-16 | 2019-10-29 | 中国电力科学研究院有限公司 | 一种考虑需求响应的生产计划调整方法及系统 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9893526B2 (en) * | 2011-03-25 | 2018-02-13 | Green Charge Networks Llc | Networked power management and demand response |
| US10423138B2 (en) * | 2017-03-06 | 2019-09-24 | Con Edison Battery Storage, Llc | Building energy storage system with planning tool |
| US10739742B2 (en) * | 2017-04-27 | 2020-08-11 | Johnson Controls Technology Company | Building energy system with stochastic model predictive control |
-
2020
- 2020-01-07 CN CN202010014848.5A patent/CN111241459B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107292416A (zh) * | 2017-05-03 | 2017-10-24 | 国网浙江省电力公司 | 基于回归预测需量的电力用户基本电费决策方法 |
| CN107832979A (zh) * | 2017-12-06 | 2018-03-23 | 浙江大学 | 一种考虑能量梯级利用的工厂综合能源系统经济优化调度方法 |
| CN110390443A (zh) * | 2018-04-16 | 2019-10-29 | 中国电力科学研究院有限公司 | 一种考虑需求响应的生产计划调整方法及系统 |
| CN108808663A (zh) * | 2018-06-12 | 2018-11-13 | 浙江大学 | 一种基于多能互补的工业用户热需求响应方法 |
| CN108898290A (zh) * | 2018-06-13 | 2018-11-27 | 广州供电局有限公司 | 一种面向多参与主体的工业园区综合能源系统互动方法 |
Non-Patent Citations (3)
| Title |
|---|
| Frederik R.Residential demand response of thermostatically controlled loads using batch reinforcement learning.《IEEE Transactions on Smart Grid》.2017,第8卷(第5期),第2149-2159页. * |
| 姜子卿 ; 刘育权 ; 艾芊 ; 蔡莹 ; 马溪原 ; 郭晓斌 ; .面向多主体的工业园区综合能源系统互动机制.南方电网技术.2018,(第03期),第23-32页. * |
| 王莉.考虑上级电网与工业园区互动的 综合需求响应.《电力建设》.2019,第40卷(第9期),第52-63页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111241459A (zh) | 2020-06-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | Coordinated optimization scheduling operation of integrated energy system considering demand response and carbon trading mechanism | |
| Bao et al. | A multi time-scale and multi energy-type coordinated microgrid scheduling solution—Part I: Model and methodology | |
| Wu et al. | Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers | |
| Ma et al. | Energy flow modeling and optimal operation analysis of the micro energy grid based on energy hub | |
| WO2014034391A1 (ja) | エネルギー管理システム、サーバ、エネルギー管理方法および記憶媒体 | |
| CN111553544B (zh) | 基于一致性算法的工业园区分布式综合需求响应方法 | |
| CN106055773A (zh) | 结合热网模型的多区域综合能源系统配置模型的建立方法 | |
| CN110649594A (zh) | 一种基于多能协同的工业园区综合需求响应调度方法 | |
| Zhang et al. | Bi-layer economic scheduling for integrated energy system based on source-load coordinated carbon reduction | |
| Wang et al. | Research on collaborative operation optimization of multi-energy stations in regional integrated energy system considering joint demand response | |
| CN112182915A (zh) | 一种协同促进风电消纳的优化调度方法及系统 | |
| Ghayour et al. | Optimal scheduling of electrical and thermal resources and appliances in a smart home under uncertainty | |
| Cheng et al. | A review on virtual power plants interactive resource characteristics and scheduling optimization | |
| CN114139837A (zh) | 一种考虑双层碳排放优化分配模型的区域多系统双层分散优化调度方法 | |
| CN111241459B (zh) | 一种考虑自趋优的工业园区分布式综合需求响应方法 | |
| Juyal et al. | Optimized Cost of Energy by a Home Energy Management System Employing Dynamic Power Import Limit Strategy: A Case study Approach | |
| Singabhattu et al. | Distributed energy resources optimization for demand response using MILP | |
| Wen et al. | Thermal and electrical demand response based on robust optimization | |
| CN111125611A (zh) | 面向多场景的冷-热-电微能网群两阶段优化调度方法 | |
| CN112270432B (zh) | 一种考虑多主体利益均衡的综合能源系统的能量管理方法 | |
| Deb et al. | Electric Vehicle Integrated Virtual Power Plants: A Systematic Review | |
| Gao et al. | Vehicle-to-home, vehicle-to-vehicle, and vehicle-to-grid energy systems | |
| Paul et al. | A review on residential area demand response analysis in smart grid era | |
| Wang et al. | Two-stage cooperative operation strategy for home management systems with smart appliances | |
| Barala et al. | Optimal scheduling for residential building based on virtual energy storage system |
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 |