CN110400059B - 一种综合能源系统联络线功率控制方法 - Google Patents
一种综合能源系统联络线功率控制方法 Download PDFInfo
- Publication number
- CN110400059B CN110400059B CN201910605092.9A CN201910605092A CN110400059B CN 110400059 B CN110400059 B CN 110400059B CN 201910605092 A CN201910605092 A CN 201910605092A CN 110400059 B CN110400059 B CN 110400059B
- Authority
- CN
- China
- Prior art keywords
- energy
- power
- control
- time
- ice
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000004519 manufacturing process Methods 0.000 claims abstract description 78
- 238000005457 optimization Methods 0.000 claims abstract description 44
- 230000008878 coupling Effects 0.000 claims abstract description 16
- 238000010168 coupling process Methods 0.000 claims abstract description 16
- 238000005859 coupling reaction Methods 0.000 claims abstract description 16
- 238000012937 correction Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims abstract description 8
- 238000003860 storage Methods 0.000 claims description 105
- 239000007789 gas Substances 0.000 claims description 56
- 230000005611 electricity Effects 0.000 claims description 34
- 238000004146 energy storage Methods 0.000 claims description 31
- 238000002844 melting Methods 0.000 claims description 23
- 230000008018 melting Effects 0.000 claims description 23
- 238000004806 packaging method and process Methods 0.000 claims description 18
- 238000002485 combustion reaction Methods 0.000 claims description 15
- 239000013598 vector Substances 0.000 claims description 15
- 238000005057 refrigeration Methods 0.000 claims description 14
- 230000032683 aging Effects 0.000 claims description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 12
- 230000004044 response Effects 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000005096 rolling process Methods 0.000 claims description 11
- 230000000295 complement effect Effects 0.000 claims description 10
- 238000012423 maintenance Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000009825 accumulation Methods 0.000 claims description 6
- 238000004378 air conditioning Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000001816 cooling Methods 0.000 claims description 6
- 238000001035 drying Methods 0.000 claims description 6
- 238000010438 heat treatment Methods 0.000 claims description 6
- 239000003345 natural gas Substances 0.000 claims description 6
- 230000006870 function Effects 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 239000000446 fuel Substances 0.000 claims description 3
- 239000011265 semifinished product Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 2
- 230000001276 controlling effect Effects 0.000 abstract 1
- 230000001105 regulatory effect Effects 0.000 abstract 1
- 238000007726 management method Methods 0.000 description 20
- 238000009499 grossing Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 238000010248 power generation Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 230000009194 climbing Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- 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/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] 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
-
- 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
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
-
- 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)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Power Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提供了一种综合能源系统联络线功率控制方法,(1)建立能源供应与生产动态耦合的综合能源优化调度模型,计算获取日前能源生产优化计划及联络线功率计划参考值;(2)基于综合能源系统中各能源子系统实时功率平衡方程,构建综合能源系统参与联络线功率控制的状态空间模型;(3)将联络线功率控制的状态空间模型等效转换为二次规划模型,并采用序列二次规划算法对模型求解,获取当前时段向后T个控制时段的可调节设备功率调节量;(4)下发当前时段向后第1个时段的功率调节量,并采样综合能源系统实时状态进行反馈校正。本发明实现复杂时变环境下的综合能源系统联络线功率精确控制,保证综合能源系统联络线功率满足上级电网部门的考核。
Description
技术领域
本发明涉及一种综合能源系统联络线功率控制方法,属于综合能源系统能量管理技术领域。
背景技术
以工业园区为主体的综合能源系统具有区域集中、能源运行设备种类多、供能规律明显、资源互补等特点,同时园区负荷稳定且自动化水平高,与商业、居民用户相比具有更大的响应容量,是能源市场改革背景下参与互动的重要需求响应资源。工业园区综合能源系统参与能源市场调节时,其能量管理中心需要对内部多种能源进行统一的优化调度和管理,并基于实时电价做出快速响应。但在实际运行中,由于光伏、风电等分布式能源的出力和负荷的用电功率预测误差较大,加之其他随机性的扰动因素,导致综合能源系统与上级配电网的联络线实际交互功率与计划值存在较大偏差,而不能实现平稳接入,进而难以满足配电网对园区综合能源系统联络线功率的考核结算需求,带来较大惩罚费用。因此,制定有效的联络线交换功率控制方案,以满足能源市场背景下的功率交互具有重要理论和现实意义。
目前已有相关针对综合能源系统或者微电网的联络线控制方案的研究展开,如专利CN201310234708.9“一种新能源并网联络线功率的平滑控制方法及装置”提出利用储能电池与燃机之间的协调配合,对新能源输出功率中不同频段的波动成分分别补偿,进而保障联络线功率的平滑;专利CN201210056752.0“一种微网联络线功率波动平滑控制方法”中通过实时采集可再生能源功率、不可控型家居设备的负荷、室内温度、温控型家居设备的开关状态及额定功率,结合微网联络线出力的发电计划需求,控制温控型家居设备的聚合负荷来实现联络线功率跟踪;专利CN201611031730.3“一种平抑微网联络线功率波动的空调负荷协调控制方法”提出基于低通滤波器原理,计算空调聚合功率控制目标,然后针对空调个体,基于市场控制,将空调聚合功率控制目标协调分配至各个空调负荷来跟踪联络线功率计划。文献“基于多智能体系统的微电网联络线潮流精确控制”、文献“采用居民温控负荷控制的微网联络线功率波动平滑方法”也针对微电网联络线功率控制分别提出了分布式多智能体控制算法、负荷需求响应队列控制算法等策略,但总体来看这些方案均是滞后的控制调节过程,不能前瞻联络线功率的变化趋势而进行超前调控,会导致出现机组功率响应不及时以及联络线控制误差大等问题。虽然专利CN201710236059.4‘一种平滑微网联络线功率的储能前瞻分布式控制方法’考虑了储能控制的前瞻性而进行联络线功率平滑的超前调节,但其针对的只是纯电网络,而未涉及多能源网络及设备的互补调节能力,也未考虑多个子能源系统共同参与联络线调节的差异化特性,存在一定的局限性。
发明内容
本发明的技术解决问题:提出一种基于模型预测控制的工业园区综合能源系统联络线功率控制方法,确保消除综合能源系统中各类不确定因素干扰的影响,也可有效避免因为机组出力的爬坡率限制或者热、冷设备响应慢而导致的联络线功率控制不及时问题,保障实现综合能源系统更好的参与能源市场互动及辅助服务。
本发明技术解决方案:一种综合能源系统联络线功率控制方法,以工业园区综合能源系统日前经济调度方案为基准参考,结合综合能源系统中各能源子系统实时功率平衡方程,构建工业园区综合能源系统参与联络线功率控制的状态空间模型,配合可再生能源、负荷的超短期功率值,实现对联络线功率变化情况的提前感知;接着将所建立的状态空间模型等效转换为二次规划模型进行快速求解,并通过不断滚动优化和采样实际运行状态的反馈校正,确保消除综合能源系统中各类不确定因素干扰的影响,同时基于超前联络线控制,也可有效避免因为机组出力的爬坡率限制或者热、冷设备响应慢而导致的联络线功率控制不及时问题,保障实现综合能源系统更好的参与能源市场互动及辅助服务。
本发明具体主要包括以下步骤:
(1)以工业园区综合能源系统运行总成本最低为优化目标,同时满足冷、热、电、气多能源子系统的功率平衡约束、能源互补耦合约束,建立能源供应与生产动态耦合的工业园区综合能源日前优化调度模型,计算获取日前能源生产优化计划及联络线功率计划参考值;
(2)以日前计算得到的能源生产优化计划及联络线功率计划参考值为基础,基于工业园区综合能源系统中冷、热、电、气各能源子系统实时功率平衡方程,构建工业园区综合能源系统参与联络线功率控制的状态空间模型;
(3)以工业园区综合能源系统日前优化得到的联络线计划与实际联络线控制功率之间的误差最小为目标,将上述状态空间模型等效转换为二次规划模型,并采用序列二次规划算法对二次规划模型求解,获取当前时段向后个NC控制时段的可调节设备功率调节量,其中NC为单次控制计算时长;
(4)下发当前时段向后第1个时段的功率调节量,并采样工业园区综合能源系统实时状态进行反馈校正,重复上述步骤(3)-步骤(4)进行滚动优化控制,直至达到总控制时长N,其中N为预先设置的联络线功率控制总时长。
所述步骤(1)中,以工业园区综合能源系统运行总成本最低为优化目标,同时满足冷、热、电、气多能源子系统的功率平衡约束、能源互补耦合约束,建立能源供应与生产动态耦合的工业园区综合能源日前优化调度模型,计算获取日前能源生产优化计划及联络线功率计划参考值。
步骤(1-1):根据公式(1)确定工业园区综合能源系统运行成本最小化的目标函数:
其中,Nt表示日前优化周期的时长;Cg,sum表示总购电成本,主要指分布式能源总发电量与各能源子系统总用电量之间由于存在差值造成的向外部电网购电的费用;Ce,sum表示总供能费用,主要包括微燃机的开机费用和运行过程所需的天然气燃料费用、燃气锅炉运行过程所需的天然气燃料费用、冰蓄冷系统的运行维护费用以及空气压缩系统的运行维护费用;Cp,sum表示总生产费用,主要由生产线运行成本、老化机和生产线机组的运行维护费用三部分构成。
步骤(1-2):确定工业园区综合能源运行的约束条件,主要包括能源互补耦合约束、各能源功率平衡约束:
1)能源互补耦合约束描述如公式(2)所示:
其中,Lc,t、Lh,t、Lp,t和La,t分别表示t时刻冷能,热能,电能和气压能的实际供应量,βc、βh、βpl和βa表示冷能,热能,电能和气压能的供应效率,Pc,t、Ph,t、Pp,t和Pa,t表示t时刻的供冷设备,供热设备,生产设备和供气设备的出力功率,αcp、αhp、αap、αpc、αph、αpa分别表示冷-电,热-电,气电,电-冷,电-热和电-气耦合系数,Bc,t、Bh,t、Bp,t和Ba,t表示冷能,热能,电能和气压能的基本负荷量。
2)各能源功率平衡约束如公式(3)所示:
式中:表示冰蓄冷系统的融冰效率,Mis,t表示t时刻冰蓄冷系统的融冰出力功率;ηae为工业空调的冷能供应效率;Pae,t表示t时刻工业空调的出力大小,Mt表示工业园区内的所有冷负荷大小;MT代表微燃机的编号,NMT代表微燃机的数量,PMT,t为第MT个微燃机在t时刻的出力功率,RMT表示第MT个微燃机的额定热电比,Pbr,t是t时刻的锅炉出力功率,Hbase,t表示工业热水和采暖等基本热负荷大小,Hcl表示每条电芯生产线在烘干阶段需要的单位热量大小,Hpl表示每条电池打包生产线在干燥阶段所需的单位热量大小;pl、cl分别代表打包生产线和电芯生产线的编号,Npl、Ncl分别为打包生产线及电芯生产线的数量,upl,t、ucl,t分别表示t时刻打包生产线和电芯生产线的运行状态,为1表示开启,为0表示关闭;Dca,t为压缩空气系统在t时刻的供气量;Dpl表示电池打包生产线在传送电池包半成品时需要的单位用气量;Dcl表示表示每条电芯生产线在传送电芯时所需要的单位用气量;PWT,t表示t时刻风电机组出力功率;PPV,t是t时刻光伏出力功率;Pes,t是t时刻储电设备的运行功率,Pgrid,t表示t时刻工业园区的购电总量,am代表老化机的编号,Nam代表老化机的数量,Pam,t表示第am个老化机在t时刻的运行功率;Pcl,t是t时刻第cl条电芯生产线的运行功率大小,Ppl,t是t时刻第pl条打包生产线的运行功率大小,Pca,t表示t时刻空压机的出力功率,Pis,t表示t时刻冰蓄冷系统制冰功率。
步骤(1-3):依据上述建立的工业园区综合能源运行调度目标函数和约束条件,形成日前优化调度模型,然后调用非线性优化求解算法对日前优化调度模型求解,获取日前能源生产优化计划及联络线功率计划参考值。
所述步骤(2)中,以日前计算得到的能源生产优化计划及联络线功率计划参考值为基础,基于工业园区综合能源系统中冷、热、电、气各能源子系统实时功率平衡方程,构建工业园区综合能源系统参与联络线功率控制的状态空间模型,以如图1所示的典型工业园区综合能源系统为具体实例,该工业园区综合能源系统包括冷、热、电、气四种能源,在能源供应侧,电能主要来源于外部电网和可再生能源,冷能主要依靠工业空调和冰蓄冷系统联合供应,热能靠微燃机和锅炉联合供应,气压能完全由压缩空气系统进行供应。在需求侧,主要包括基本用电管理,温控设备管理和生产设备管理。其中基本用电管理指日常照明和办公用电等,而温控设备管理则依靠冷热能相互协调将室温维持在一定范围内,而对于生产设备的管理,则指满足生产任务和流程的前提下,综合考虑不同类型的生产设备运行需要进行多种能源的协同供应。以该系统为例构建状态空间具体包括以下步骤:
步骤(2-1):选取空压机,空调制冷、制冰机,融冰机,蓄电池储能充放电的功率,储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能储能荷电状态,室内温度以及工业园区综合能源系统与外电网的联络线交换功率构成的向量为状态变量:
x(k)=[Pca(k),Pae(k),Pis(k),Mis(k),Pes(k),Eca(k),Eis(k),SB(k),T(k),Pgrid(k)]T (4)
步骤(2-2):选取空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力增量构成控制变量向量:
u(k)=[ΔPca(k),ΔPae(k),ΔPis(k),Mis(k),ΔPes(k)]T (5)
步骤(2-3):选取风电机组及光伏的超短期预测功率增量和储气罐放气功率变化量、室外温度变化量构成扰动变量向量:
d(k)=[ΔPWT(k),ΔPPV(k),ΔDca(k),ΔTo(K)]T (6)
步骤(2-4):选取储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率构成输出变量向量:
y(k)=[Eca(k),Eis(k),SB(k),Pgrid(k)]T (7)
步骤(2-5):依据综合能源系统中冷、热、电、气各能源子系统实时功率平衡方程,列写如下的状态空间模型方程:
式中,A为状态转移矩阵,B为控制矩阵,R为扰动矩阵,C为输出矩阵;x(k)为状态变量,u(k)为控制变量,d(k)为扰动变量,y(k)为输出变量;Δt为控制时间步长,ηca为空压机的效率,为冰蓄冷系统制冰效率;为冰蓄冷系统融冰效率;σ为储能自放电率;k1、k2为园区厂房热容量模型参数。
所述步骤(3)中,以工业园区综合能源系统日前优化得到的联络线计划与实际联络线控制功率之间的误差最小为目标,将上述状态空间模型等效转换为二次规划模型,并采用序列二次规划算法对二次规划模型求解,获取当前时段向后NC个控制时段的可调节设备功率调节量,其中NC为单次控制计算时长,具体包括以下步骤:
步骤(3-1):通过对步骤(2)所构建的状态空间预测模型的公式(8)反复迭代,直到向前预测p步,得到储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率输出变量在预测时长pΔt内的预估输出值构成的向量Yf:
步骤(3-2):选取当前时刻向前pΔt时段内,储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率等输出变量的日前计划值构成向量Rda为跟踪控制目标:
步骤(3-3):以Yf和Rda之间的差值最小为目标,将对应的实时滚动优化控制转化为二次规划问题:
其中,Werr为权重系数矩阵,Pca,min和Pca,max分别为空气机最小和最大输出功率;Pae,min和Pae,max分别为冰蓄冷系统空调制冷的最小和最大运行功率;Pis,min和Pis,max分别为制冰机最小和最大输出功率;Mis,min和Mis,max分别为融冰机最小和最大输出功率;Vca为储气罐容量体积;Eis,min和Eis,max分别为冰蓄冷系统最小和最大蓄冷量;SB,min和SB,max分别为蓄电池储能最小和最大荷电状态;Ti,min和Ti,max分别为厂房控制温度的最小和最大值。
步骤(3-4):通过对式(11)所示的二次规划问题调用序列二次规划算法求解,求解得到当前时段向后NC个控制时段的空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力调节增量,其中NC为单次控制计算时长。
所述步骤(4)中,下发当前时段向后第1个时段的控制调节增量,并采样工业园区综合能源系统实时运行状态进行反馈校正,重复步骤(3)-步骤(4),直至达到总控制时长N,其中N为预先设置的联络线功率控制总时长。具体步骤描述如下:
步骤(4-1):读取当前时段向后NC个控制时段内所有空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力调节增量构成的优化控制序列;
步骤(4-2):在当前控制时刻仅将当前时刻向后的第1个时段控制指令序列下发;
步骤(4-3):依据当前控制指令序列下发后综合能源系统的实际响应情况,采样综合能源系统实时运行状态数据进行反馈,并将其作为下一个滚动优化控制的基础,更新状态空间模型的状态变量;
步骤(4-4):重复步骤(3)中的滚动优化计算以及步骤(4)的指令下发和反馈校正过程,直至达到总控制时长N为止,其中N为预先设置的联络线功率控制总时长。
本发明与现有技术相比的优点在于:
(1)本发明基于综合能源系统中冷、热、电、气各能源子系统的实时功率平衡方程,构建了综合能源系统参与联络线功率控制的状态空间模型,结合可再生能源、负荷的超短期功率预测,能实现对联络线功率变化情况的提前感知而实施超前控制,相比现有的滞后型潮流断面控制方法,可有效规避机组出力的爬坡率限制或者热、冷设备响应慢而导致的联络线功率控制不及时问题;
(2)本发明通过将综合能源系统联络线功率控制的状态空间模型与模型预测控制方法相结合,采取滚动优化和实际运行状态的不断反馈校正调节,能确保消除综合能源系统中各类不确定因素的干扰影响,相比现有的开环联络线控制方案,提高了综合能源系统联络线功率控制的鲁棒性;
(3)本发明考虑了综合能源系统中冷、热、电等多类型能源子系统和设备参与联络线功率控制的互补调节能力及差异化响应特性,相比现有的纯粹基于电网络及其设备的调节方案,具有更大的调节能力和更好的调节灵活性。
附图说明
图1为综合能源系统基本架构;
图2为本发明的综合能源系统联络线功率控制方法流程图。
具体实施方式
下面结合附图及实施例对本发明进行详细说明。
本发明的一种综合能源系统联络线功率控制方法,主要用于工业园区综合能源系统,主要包括:1)以工业园区综合能源系统运行总成本最低为优化目标,同时满足冷、热、电、气多能源子系统的功率平衡约束、能源互补耦合约束,建立能源供应与生产动态耦合的工业园区综合能源日前优化调度模型,计算获取日前能源生产优化计划及联络线功率计划参考值;2)以日前计算得到的能源生产优化计划及联络线功率计划参考值为基础,基于工业园区综合能源系统中冷、热、电、气各能源子系统实时功率平衡方程,构建工业园区综合能源系统参与联络线功率控制的状态空间模型;3)以工业园区综合能源系统日前优化得到的联络线计划与实际联络线控制功率之间的误差最小为目标,将上述状态空间模型等效转换为二次规划模型,并采用序列二次规划算法对模型求解,获取当前时段向后个NC控制时段的可调节设备功率调节量,其中NC为单次控制计算时长;4)下发当前时段向后第1个时段的功率调节量,并采样工业园区综合能源系统实时状态进行反馈校正,重复上述步骤3)-步骤4)进行滚动优化控制,直至达到总控制时长N,其中N为预先设置的联络线功率控制总时长。本发明提出的综合能源系统联络线功率控制方法,能提前感知联络线功率的变化情况,进而实现超前优化控制,避免因为机组出力的爬坡率限制或者热、冷设备响应慢而导致的联络线功率控制不及时问题,具有较好的应用前景。
具体实施流程如图2所示,包括以下步骤:
步骤(1):预测日前及超短期光伏、风力发电等可再生能源发电功率,预测日前及超短期电负荷、热负荷、气负荷、冷负荷功率,预测电价和天然气价格,设置联络线功率控制总时长N,单次控制时长NC;
步骤(2):根据日前光伏、风力发电预测功率、电负荷、热负荷、气负荷、冷负荷预测功率以及电价和天然气价格预测数据,以如图1所示的典型工业园区综合能源系统为具体实例,该工业园区综合能源系统包括冷、热、电、气四种能源,在能源供应侧,电能主要来源于外部电网和可再生能源,冷能主要依靠工业空调和冰蓄冷系统联合供应,热能靠微燃机和锅炉联合供应,气压能完全由压缩空气系统进行供应。在需求侧,主要包括基本用电管理,温控设备管理和生产设备管理。其中基本用电管理指日常照明和办公用电等,而温控设备管理则依靠冷热能相互协调将室温维持在一定范围内,而对于生产设备的管理,则指满足生产任务和流程的前提下,综合考虑不同类型的生产设备运行需要进行多种能源的协同供应。基于该综合能源系统构建能源供应与生产动态耦合的工业园区综合能源日前优化调度模型(1),计算获取日前能源生产优化计划及联络线功率计划参考值
其中,为工业园区综合能源系统日运行总成本, 为日前经济调度的部分决策变量,分别代指一天内各时段的储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态以及联络线交换功率构成的向量集。Nt表示优化周期的时长;Cg,sum表示总购电成本;Ce,sum表示总供能费用,主要包括微燃机的开机费用和运行过程所需的天然气燃料费用、燃气锅炉运行过程所需的天然气燃料费用、冰蓄冷系统的运行维护费用以及空气压缩系统的运行维护费用;Cp,sum表示总生产费用,主要由生产线运行成本、老化机和生产线机组的运行维护费用三部分构成。Lc,t、Lh,t、Lp,t和La,t分别表示t时刻冷能,热能,电能和气压能的实际供应量,βc、βh、βpl和βa表示冷能,热能,电能和气压能的供应效率,Pc,t、Ph,t、Pp,t和Pa,t表示t时刻的供冷设备,供热设备,生产设备和供气设备的出力功率,αcp、αhp、αap、αpc、αph、αpa分别表示冷-电,热-电,气-电,电-冷,电-热和电-气耦合系数,Bc,t、Bh,t、Bp,t和Ba,t表示冷能,热能,电能和气压能的基本负荷量。表示冰蓄冷系统的融冰效率,Mis,t表示t时刻冰蓄冷系统的融冰出力功率;ηae为工业空调的冷能供应效率;Pae,t表示t时刻工业空调的出力大小,Mt表示工业园区内的所有冷负荷大小;MT代表微燃机的编号,NMT代表微燃机的数量,PMT,t为第MT个微燃机在t时刻的出力功率,RMT表示第MT个微燃机的额定热电比,Pbr,t是t时刻的锅炉出力功率,Hbase,t表示工业热水和采暖等基本热负荷大小,Hcl表示每条电芯生产线在烘干阶段需要的单位热量大小,Hpl表示每条电池打包生产线在干燥阶段所需的单位热量大小;pl、cl分别代表打包生产线和电芯生产线的编号,Npl、Ncl分别为打包生产线及电芯生产线的数量,upl,t表示t时刻第pl条打包生产线的运行状态,为1表示开启,为0表示关闭,ucl,t表示t时刻第cl条电芯生产线的运行状态,为1表示开启,为0表示关闭;Dca,t为压缩空气系统在t时刻的供气量;Dpl表示电池打包生产线在传送电池包半成品时需要的单位用气量;Dcl表示表示每条电芯生产线在传送电芯时所需要的单位用气量;PWT,t表示t时刻风电机组出力功率;PPV,t是t时刻光伏出力功率;Pes,t是t时刻储电设备的运行功率,Pgrid,t表示t时刻工业园区的购电总量,am代表老化机的编号,Nam代表老化机的数量,Pam,t表示第am个老化机在t时刻的运行功率;Pcl,t是t时刻第cl条电芯生产线的运行功率大小,Ppl,t是t时刻第pl条打包生产线的运行功率大小,Pca,t表示t时刻空压机的出力功率,Pis,t表示t时刻冰蓄冷系统制冰功率。
步骤(3):基于工业园区综合能源系统中冷、热、电、气各能源子系统实时功率平衡方程,构建工业园区综合能源系统参与联络线功率控制的状态空间模型,以如图1所示的典型工业园区综合能源系统为具体实例,该工业园区综合能源系统包括冷、热、电、气四种能源,在能源供应侧,电能主要来源于外部电网和可再生能源,冷能主要依靠工业空调和冰蓄冷系统联合供应,热能靠微燃机和锅炉联合供应,气压能完全由压缩空气系统进行供应。在需求侧,主要包括基本用电管理,温控设备管理和生产设备管理。其中基本用电管理指日常照明和办公用电等,而温控设备管理则依靠冷热能相互协调将室温维持在一定范围内,而对于生产设备的管理,则指满足生产任务和流程的前提下,综合考虑不同类型的生产设备运行需要的多种能源的协同供应。构建得到的状态空间模型如式(2)所示;
其中,Pca(k)、Pae(k)、Pis(k)、Mis(k)、Pes(k)、Eca(k)、Eis(k)、SB(k)、T(k)、Pgrid(k)分别为k时段空压机,空调制冷、制冰机,融冰机,蓄电池储能充放电的功率,储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能储能荷电状态,室内温度以及工业园区综合能源系统与外电网的联络线交换功率;ΔPca(k)、ΔPae(k)、ΔPis(k)、ΔMis(k)、ΔPes(k)分别为空压机,空调制冷、制冰机,融冰机和蓄电池储能在k时段的出力增量;ΔPWT(k)、ΔPPV(k)、ΔDca(k)、ΔTo(k)分别为风电机组及光伏的超短期预测功率增量、储气罐放气功率变化量以及室外温度的变化量;Δt为控制时间步长,ηca为空压机的效率,为冰蓄冷系统制冰效率;为冰蓄冷系统融冰效率;σ为储能自放电率;k1、k2为园区厂房的热容量模型参数。
步骤(4):以工业园区综合能源系统日前优化得到的联络线计划与实际联络线控制功率之间的误差最小为目标,将上述状态空间模型等效转换为二次规划模型(3),并结合风电机组及光伏的超短期预测功率增量、储气罐放气功率增量以及室外温度的变化量,采用序列二次规划算法对二次规划模型求解,获取当前时段向后NC个控制时段的空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力调节增量,其中,NC为单次控制计算时长;
其中,J为联络线功率控制的性能指标;Rda为储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率等输出变量在预测时长NPΔt内的日前计划值构成的向量;Yf为储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率输出变量在预测时长NPΔt内的预估输出值构成的向量;Werr为权重系数矩阵,Pca,min和Pca,max分别为空压机最小和最大输出功率;Pae,min和Pae,max分别为冰蓄冷系统空调制冷的最小和最大运行功率;Pis,min和Pis,max分别为制冰机最小和最大输出功率;Mis,min和Mis,max分别为融冰机最小和最大输出功率;Vca为储气罐容量体积;Eis,min和Eis,max分别为冰蓄冷系统最小和最大蓄冷量;SB,min和SB,max分别为蓄电池储能最小和最大荷电状态;Ti,min和Ti,max分别为厂房控制温度的最小和最大值。
步骤(5):下发当前时刻向后的第1个控制时长的空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力调节增量控制序列;
步骤(6):依据当前控制指令序列下发后系统的实际响应情况,采样综合能源系统实时运行状态数据反馈,并将其作为下一个滚动优化控制的基础,更新状态空间模型的状态变量;
步骤(7):重复步骤(3)到步骤(6),直至达到总控制时长N为止,其中N为预先设置的联络线功率控制总时长。
提供以上实施步骤仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。
Claims (3)
1.一种综合能源系统联络线功率控制方法,用于工业园区综合能源系统,其特征在于,包括以下步骤:
(1)以工业园区综合能源系统运行总成本最低为优化目标,同时满足冷、热、电和气多能源子系统的功率平衡约束和能源互补耦合约束,建立能源供应与生产动态耦合的工业园区综合能源日前优化调度模型,计算获取日前能源生产优化计划及联络线功率计划参考值;
(2)以日前计算得到的能源生产优化计划及联络线功率计划参考值,基于工业园区综合能源系统中冷、热、电和气各能源子系统实时功率平衡方程,构建工业园区综合能源系统参与联络线功率控制的状态空间模型;
(3)以工业园区综合能源系统日前优化得到的联络线计划与实际联络线控制功率之间的误差最小为目标,将上述状态空间模型等效转换为二次规划模型,并采用序列二次规划算法对二次规划模型求解,获取当前时段向后个NC控制时段的可调节设备功率调节量,其中NC为单次控制计算时长;
(4)下发当前时段向后第1个时段的可调设备功率调节量,并采样工业园区综合能源系统实时状态进行反馈校正,重复上述步骤(3)-步骤(4)进行滚动优化控制,直至达到总控制时长N,其中N为预先设置的联络线功率控制总时长;
所述步骤(2)具体包括以下步骤:
步骤(2-1):选取空压机,空调制冷或制冰机,融冰机,蓄电池储能充放电的功率,储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态,室内温度以及综合能源系统与外电网的联络线交换功率构成的向量为状态变量:
x(k)=[Pca(k),Pae(k),Pis(k),Mis(k),Pes(k),Eca(k),Eis(k),SB(k),T(k),Pgrid(k)]T (4)
步骤(2-2):选取空压机,空调制冷或制冰机,融冰机和蓄电池储能的出力增量构成控制变量向量:
u(k)=[△Pca(k),△Pae(k),△Pis(k),Mis(k),△Pes(k)]T (5)
步骤(2-3):选取风电机组及光伏的超短期预测功率增量和储气罐放气功率变化量、室外温度变化量构成扰动变量向量:
d(k)=[△PWT(k),△PPV(k),△Dca(k),△To(K)]T (6)
步骤(2-4):选取储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率构成输出变量向量:
y(k)=[Eca(k),Eis(k),SB(k),Pgrid(k)]T (7)
步骤(2-5):依据综合能源系统中冷、热、电、气各能源子系统实时功率平衡方程,列写如下的状态空间模型方程:
式中,A为状态转移矩阵,B为控制矩阵,R为扰动矩阵,C为输出矩阵;x(k)为状态变量,u(k)为控制变量,d(k)为扰动变量,y(k)为输出变量;△t为控制时间步长,ηca为空压机的效率,为冰蓄冷系统制冰效率;为冰蓄冷系统融冰效率;σ为储能自放电率;k1、k2为厂房热容量模型参数;
所述步骤(3)中,以工业园区综合能源系统日前优化得到的联络线计划与实际联络线控制功率之间的误差最小为目标,将上述状态空间模型等效转换为二次规划模型,并采用序列二次规划算法对二次规划模型求解,获取当前时段向后NC个控制时段的可调节设备功率调节量,其中NC为单次控制计算时长,具体包括以下步骤:
步骤(3-1):通过对步骤(2)所构建的状态空间预测模型的公式(8)反复迭代,直到向前预测p步,得到储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率输出变量在预测时长pΔt内的预估输出值构成的向量Yf:
步骤(3-2):选取当前时刻向前pΔt时段内,储气罐储气量,冰蓄冷系统蓄冷量,蓄电池储能荷电状态和联络线交换功率输出变量的日前计划值构成向量Rda为跟踪控制目标:
步骤(3-3):以Yf和Rda之间的差值最小为目标,将对应的实时滚动优化控制转化为二次规划问题:
其中,Werr为权重系数矩阵,Pca,min和Pca,max分别为空气机最小和最大输出功率;Pae,min和Pae,max分别为冰蓄冷系统空调制冷的最小和最大运行功率;Pis,min和Pis,max分别为制冰机最小和最大输出功率;Mis,min和Mis,max分别为融冰机最小和最大输出功率;Vca为储气罐容量体积;Eis,min和Eis,max分别为冰蓄冷系统最小和最大蓄冷量;SB,min和SB,max分别为蓄电池储能最小和最大荷电状态;Ti,min和Ti,max分别为厂房控制温度的最小和最大值;
步骤(3-4):通过对式(11)所示的二次规划问题调用序列二次规划算法求解,求解得到当前时段向后NC个控制时段的空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力调节增量。
2.根据权利要求1所述的综合能源系统联络线功率控制方法,其特征在于:所述步骤(1)具体实现如下:
步骤(1-1):根据以下公式(1)确定工业园区综合能源系统运行成本最小化的目标函数:
其中,Nt表示日前优化周期的时长;Cg,sum表示总购电成本,指分布式能源总发电量与各能源子系统总用电量之间由于存在差值造成的向外部电网购电的费用;Ce,sum表示总供能费用,包括微燃机的开机费用和运行过程所需的天然气燃料费用、燃气锅炉运行过程所需的天然气燃料费用、冰蓄冷系统的运行维护费用以及空气压缩系统的运行维护费用;Cp,sum表示总生产费用,由生产线运行成本、老化机和生产线机组的运行维护费用三部分构成;
步骤(1-2):确定工业园区综合能源运行的约束条件,约束条件包括能源互补耦合约束、各能源功率平衡约束:
1)能源互补耦合约束如公式(2)所示:
其中,Lc,t、Lh,t、Lp,t和La,t分别表示t时刻冷能,热能,电能和气压能的实际供应量,βc、βh、βpl和βa表示冷能,热能,电能和气压能的供应效率,Pc,t、Ph,t、Pp,t和Pa,t表示t时刻的供冷设备,供热设备,生产设备和供气设备的出力功率,αcp、αhp、αap、αpc、αph、αpa分别表示冷-电,热-电,气-电,电-冷,电-热和电-气耦合系数,Bc,t、Bh,t、Bp,t和Ba,t表示冷能,热能,电能和气压能的基本负荷量;
2)各能源功率平衡约束如公式(3)所示:
式中:表示冰蓄冷系统的融冰效率,Mis,t表示t时刻冰蓄冷系统的融冰出力功率;ηae为工业空调的冷能供应效率;Pae,t表示t时刻工业空调的出力大小,Mt表示工业园区内的所有冷负荷大小;MT代表微燃机的编号,NMT代表微燃机的数量,PMT,t为第MT个微燃机在t时刻的出力功率,RMT表示第MT个微燃机的额定热电比,Pbr,t是t时刻的锅炉出力功率,Hbase,t表示工业热水和采暖基本热负荷大小,Hcl表示每条电芯生产线在烘干阶段需要的单位热量大小,Hpl表示每条电池打包生产线在干燥阶段所需的单位热量大小;pl、cl分别代表打包生产线和电芯生产线的编号,Npl、Ncl分别为打包生产线及电芯生产线的数量,upl,t、ucl,t分别表示t时刻打包生产线和电芯生产线的运行状态,为1表示开启,为0表示关闭;Dca,t为压缩空气系统在t时刻的供气量;Dpl表示电池打包生产线在传送电池包半成品时需要的单位用气量;Dcl表示每条电芯生产线在传送电芯时所需要的单位用气量;PWT,t表示t时刻风电机组出力功率;PPV,t是t时刻光伏出力功率;Pes,t是t时刻储电设备的运行功率,Pgrid,t表示t时刻工业园区的购电总量,am代表老化机的编号,Nam代表老化机的数量,Pam,t表示第am个老化机在t时刻的运行功率;Pcl,t是t时刻第cl条电芯生产线的运行功率大小,Ppl,t是t时刻第pl条打包生产线的运行功率大小,Pca,t表示t时刻空压机的出力功率,Pis,t表示t时刻冰蓄冷系统制冰功率;
步骤(1-3):依据上述建立的工业园区综合能源运行调度目标函数和约束条件,形成日前优化调度模型,然后调用非线性优化求解算法对日前优化调度模型求解,获取日前能源生产优化计划及联络线功率计划参考值。
3.根据权利要求1所述的综合能源系统联络线功率控制方法,其特征在于:所述步骤(4)中具体步骤如下:
步骤(4-1):读取当前时段向后NC个控制时段内所有空压机,空调制冷、制冰机,融冰机和蓄电池储能的出力调节增量构成的优化控制序列;
步骤(4-2):在当前控制时刻仅将当前时刻向后的第1个时段控制指令序列下发;
步骤(4-3):依据当前控制指令序列下发后系统的实际响应情况,采样综合能源系统实时运行状态数据进行反馈,并将其作为下一个滚动优化控制的基础,更新状态空间模型的状态变量;
步骤(4-4):重复步骤(3)中的滚动优化计算以及步骤(4)的指令下发和反馈校正过程,直至达到总控制时长N,其中N为预先设置的联络线功率控制总时长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910605092.9A CN110400059B (zh) | 2019-07-05 | 2019-07-05 | 一种综合能源系统联络线功率控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910605092.9A CN110400059B (zh) | 2019-07-05 | 2019-07-05 | 一种综合能源系统联络线功率控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110400059A CN110400059A (zh) | 2019-11-01 |
CN110400059B true CN110400059B (zh) | 2021-11-02 |
Family
ID=68322743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910605092.9A Active CN110400059B (zh) | 2019-07-05 | 2019-07-05 | 一种综合能源系统联络线功率控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110400059B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111313400B (zh) * | 2019-11-11 | 2022-07-12 | 国网吉林省电力有限公司 | 一种基于鲁棒修正的多能虚拟电厂运行参数聚合方法 |
CN112819191A (zh) * | 2019-11-18 | 2021-05-18 | 国网能源研究院有限公司 | 一种面向工程应用的园区综合能源系统运行优化分析方法 |
CN111799806A (zh) * | 2020-06-22 | 2020-10-20 | 长沙理工大学 | 一种基于云边协同和多智能体的降压运行高峰负荷削减方法及装置 |
CN111969623B (zh) * | 2020-08-17 | 2022-03-29 | 中国科学院电工研究所 | 一种电网主导型电压源换流器幅相控制方法 |
CN112366704B (zh) * | 2020-11-13 | 2022-09-27 | 中国科学院电工研究所 | 一种基于激励需求响应的综合能源系统联络线功率控制方法 |
CN113344736B (zh) * | 2021-05-21 | 2023-08-18 | 温州电力设计有限公司 | 一种园区级综合能源系统及其控制方法 |
CN113937788B (zh) * | 2021-10-11 | 2024-02-02 | 国网山东省电力公司营销服务中心(计量中心) | 一种电力负荷侧资源调度方法及系统 |
ES2955464A1 (es) * | 2023-09-13 | 2023-12-01 | Gtd Sist De Informacion Sau | Metodo para la planificacion y optimizacion de almacenes termicos en instalaciones industriales de produccion de agua fria y/o caliente |
CN117973948B (zh) * | 2024-04-01 | 2024-06-21 | 福建昊智检测技术有限公司 | 一种可基于数字式变电站电能质量分析系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101917002A (zh) * | 2010-07-16 | 2010-12-15 | 清华大学 | 基于特高压电网联络线运行状态的电压控制方法 |
CN104268806A (zh) * | 2014-11-03 | 2015-01-07 | 四川慧盈科技有限责任公司 | 一种微电网功率监控系统 |
CN107069776A (zh) * | 2017-04-12 | 2017-08-18 | 东南大学 | 一种平滑微网联络线功率的储能前瞻分布式控制方法 |
CN107958300A (zh) * | 2017-08-21 | 2018-04-24 | 浙江工业大学 | 一种考虑互动响应的多微网互联运行协调调度优化方法 |
CN108494015A (zh) * | 2018-02-09 | 2018-09-04 | 中国科学院电工研究所 | 一种源-荷-储协调互动的综合能源系统设计方法 |
CN108932560A (zh) * | 2018-06-13 | 2018-12-04 | 天津大学 | 基于模型预测控制的园区综合能源系统优化调度方法 |
CN109038589A (zh) * | 2018-08-02 | 2018-12-18 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种多省区电网协调运行生产模拟方法 |
CN109327042A (zh) * | 2018-09-27 | 2019-02-12 | 南京邮电大学 | 一种微电网多能源联合优化调度方法 |
CN109634119A (zh) * | 2018-12-24 | 2019-04-16 | 浙江工业大学 | 一种基于日内滚动优化的能源互联网优化控制方法 |
CN109919399A (zh) * | 2019-04-15 | 2019-06-21 | 中国科学院电工研究所 | 一种综合能源系统日前经济调度方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9507367B2 (en) * | 2012-04-09 | 2016-11-29 | Clemson University | Method and system for dynamic stochastic optimal electric power flow control |
US9625887B2 (en) * | 2012-08-03 | 2017-04-18 | Carnegie Mellon University | Contingency screening in multi-control area electrical power systems using coordinated DC power flow |
-
2019
- 2019-07-05 CN CN201910605092.9A patent/CN110400059B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101917002A (zh) * | 2010-07-16 | 2010-12-15 | 清华大学 | 基于特高压电网联络线运行状态的电压控制方法 |
CN104268806A (zh) * | 2014-11-03 | 2015-01-07 | 四川慧盈科技有限责任公司 | 一种微电网功率监控系统 |
CN107069776A (zh) * | 2017-04-12 | 2017-08-18 | 东南大学 | 一种平滑微网联络线功率的储能前瞻分布式控制方法 |
CN107958300A (zh) * | 2017-08-21 | 2018-04-24 | 浙江工业大学 | 一种考虑互动响应的多微网互联运行协调调度优化方法 |
CN108494015A (zh) * | 2018-02-09 | 2018-09-04 | 中国科学院电工研究所 | 一种源-荷-储协调互动的综合能源系统设计方法 |
CN108932560A (zh) * | 2018-06-13 | 2018-12-04 | 天津大学 | 基于模型预测控制的园区综合能源系统优化调度方法 |
CN109038589A (zh) * | 2018-08-02 | 2018-12-18 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种多省区电网协调运行生产模拟方法 |
CN109327042A (zh) * | 2018-09-27 | 2019-02-12 | 南京邮电大学 | 一种微电网多能源联合优化调度方法 |
CN109634119A (zh) * | 2018-12-24 | 2019-04-16 | 浙江工业大学 | 一种基于日内滚动优化的能源互联网优化控制方法 |
CN109919399A (zh) * | 2019-04-15 | 2019-06-21 | 中国科学院电工研究所 | 一种综合能源系统日前经济调度方法及系统 |
Non-Patent Citations (3)
Title |
---|
基于DMPC的联络线功率波动协调抑制方法;代仕勇等;《广东电力》;20141125;第27卷(第11期);第64-70页 * |
考虑复杂生产约束的电池生产工业园区能源网络与生产管理综合优化运行;马昕等;《电网技术》;20180919;第42卷(第11期);第3566-3574页 * |
考虑空调负荷的微网联络线功率波动;潘浩等;《南方电网技术》;20160820;第10卷(第8期);第56-60页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110400059A (zh) | 2019-11-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110400059B (zh) | 一种综合能源系统联络线功率控制方法 | |
CN109063925B (zh) | 一种计及负荷聚合商的区域综合能源系统优化运行方法 | |
CN109004686B (zh) | 一种考虑冰蓄冷空调多模式的冷热电联供型微网系统 | |
CN111950809B (zh) | 基于主从博弈的综合能源系统分层分区优化运行方法 | |
CN110619110B (zh) | 一种含热泵的综合能源系统协调运行优化方法 | |
CN109752953B (zh) | 一种集成电制冷机的楼宇供能系统模型预测调控方法 | |
CN111400641A (zh) | 一种含蓄热式电采暖的综合能源系统日前优化调度方法 | |
CN111404183B (zh) | 区域综合能源系统多元储能协同配置方法、程序、系统与应用 | |
CN112990523B (zh) | 区域综合能源系统分层优化运行方法 | |
CN110807588B (zh) | 一种多能源耦合综合能源系统的优化调度方法 | |
CN109617142A (zh) | 一种cchp型微电网多时间尺度优化调度方法及系统 | |
CN110474370B (zh) | 一种空调可控负荷、光伏储能系统的协同控制系统及方法 | |
CN110880772A (zh) | 一种基于工业园区负荷聚合的售电公司响应电网控制方法 | |
CN112700094A (zh) | 一种基于mpc及loddlc的综合能源系统多时间尺度优化调度方法 | |
CN105955931A (zh) | 面向高密度分布式光伏消纳的区域能源网络优化调度方法 | |
KR20130074045A (ko) | 마이크로그리드 기반의 수용가 에너지 관리 방법 및 시스템 | |
CN117081143A (zh) | 促进分布式光伏就地消纳的园区综合能源系统协调优化运行方法 | |
CN115912419A (zh) | 一种基于负荷聚合商的储热-储电协同调度方法 | |
Luo et al. | A two-stage energy management strategy for CCHP microgrid considering house characteristics | |
CN105631557B (zh) | 考虑蓄冰空调并具有电冷耦合特征的微电网优化调度方法 | |
CN113673739B (zh) | 分布式综合能源系统的多时空尺度协同优化运行方法 | |
CN116227739B (zh) | 计及能效的含热泵与电热混合储能的建筑微能源网双层优化配置方法 | |
CN113659569A (zh) | 一种电力系统日前优化调度方法及系统 | |
Palma et al. | Application of a demand-response-optimized electrical load profile to a plant supplying an energy micro-community | |
CN115378002A (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 |