CN109301817B - 一种考虑需求响应的多时间尺度源网荷协调调度方法 - Google Patents
一种考虑需求响应的多时间尺度源网荷协调调度方法 Download PDFInfo
- Publication number
- CN109301817B CN109301817B CN201811128233.4A CN201811128233A CN109301817B CN 109301817 B CN109301817 B CN 109301817B CN 201811128233 A CN201811128233 A CN 201811128233A CN 109301817 B CN109301817 B CN 109301817B
- Authority
- CN
- China
- Prior art keywords
- load
- scheduling
- day
- time
- output
- 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 100
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000008569 process Effects 0.000 claims abstract description 6
- 230000003247 decreasing effect Effects 0.000 claims description 32
- 230000005540 biological transmission Effects 0.000 claims description 19
- 230000009194 climbing Effects 0.000 claims description 17
- 238000010248 power generation Methods 0.000 claims description 15
- 238000005457 optimization Methods 0.000 claims description 14
- 230000000452 restraining effect Effects 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- AGBQKNBQESQNJD-ZETCQYMHSA-N (S)-lipoic acid Chemical compound OC(=O)CCCC[C@H]1CCSS1 AGBQKNBQESQNJD-ZETCQYMHSA-N 0.000 description 4
- 230000002354 daily effect Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000002776 aggregation Effects 0.000 description 3
- 238000004220 aggregation Methods 0.000 description 3
- 239000003795 chemical substances by application Substances 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005096 rolling process Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
Images
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
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Power Engineering (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- Primary Health Care (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Strategic Management (AREA)
- Water Supply & Treatment (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Control Of Eletrric Generators (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种考虑需求响应的多时间尺度源网荷协调调度方法,属于电网调度技术领域。所述方法根据需求响应和调度的时间特性,将整体调度过程分为日前调度、日内调度和实时调度三个阶段,建立网络约束条件下基于多时间尺度需求响应的源网荷经济调度模型,包括日前调度模型、日内调度模型、实时调度模型,并基于调度模型的求解结果进行调度。与现有日前‑实时调度模型进行对比,利用本发明所提模型和方法能够得到更低的调度成本和更低的调度响应量。
Description
技术领域
本发明属于电网调度技术领域,具体涉及一种考虑需求响应的多时间尺度源网荷协调调度方法。
背景技术
随着风电等新能源的不断接入,电网的安全可靠运行面临着挑战。风电接入电网后,其间歇性和不确定性给电力***的调度和运行增加了难度。风电功率预测存在误差,且误差大小与预测时间有关,越接近预测点时误差越小。日前风电预测的误差一般为25%-40%,日内4h风电预测误差为10%-20%,1h风电预测误差则在10%以内。可见对于时间尺度的不同,风电预测的误差也会不一样,因此在调度过程中需要考虑到多个时间尺度。由于日前调度是提前制定的调度计划,受新能源和负荷预测精度以及电力***实际状况变化影响,要通过实时调度来调整日前调度计划的客观技术需求。现有模型与方法中对于需求响应(Demand Response,DR)资源的调度多集中在某一个固定的时间尺度上,忽略了DR资源具有的多时间尺度特性,不能充分调用所有的DR资源。
发明内容
发明目的:针对现有技术的不足,本发明提出一种考虑需求响应的多时间尺度源网荷协调调度方法,能够充分调动DR资源,实现DR和电源侧、电网侧资源的协调优化。
技术方案:为实现以上目的,本发明所采用以下技术方案:
一种考虑需求响应的多时间尺度源网荷协调调度方法,根据需求响应和调度的时间特性,将整体调度过程分为日前调度、日内调度和实时调度三个阶段,建立网络约束条件下基于多时间尺度需求响应的源网荷经济调度模型,并基于调度模型的求解结果制定调度计划,具体包括以下步骤:
(1)在已知单个风电场预测出力、负荷预测曲线和A类负荷响应潜力的情况下,以发电机成本、弃风成本和负荷响应成本总和最小为目标,建立考虑需求响应的日前调度模型,确定日前调度计划中常规发电机组出力、风电场出力和A类负荷的调整量,以第一时间间隔制定日前调度计划;
(2)将日前调度的优化决策变量中的发电机组出力和A类负荷的调整量作为日内调度的已知量,输入风电场的日内预测曲线和B类负荷响应潜力,以总成本最小为目标建立考虑需求响应的日内调度模型,进而确定日内调度计划中发电机组出力、风电场出力和B类负荷调整量,以第二时间间隔制定日内调度计划;
(3)将日内调度的优化决策变量中的发电机组出力和B类负荷调整量作为实时调度的已知量,输入风电场的实时预测曲线和当前时刻C类负荷响应潜力,以总成本最小为目标建立考虑需求响应的实时调度模型,进而确定实时调度计划中发电机组出力、风电场出力和C类负荷调整量,以第三时间间隔制定实时调度计划;
其中,第一时间间隔>第二时间间隔>第三时间间隔;
(4)调度中心向参与DR计划的智能家电下发状态指令,智能家电进行多次响应至满足***C类负荷调整量的需求;C类负荷响应动作之后,更新下一阶段的C类负荷总功率,对下一时间段进行实时调度。
作为优选,所述步骤(1)中日前调度将第一时间间隔划分为若干个时段,计算每个时段的日前优化决策量,日前调度模型具体如下:
日前调度模型的目标函数为:
其中,NT为调度时段个数;NG为发电机个数;CG,i,t和PG,i,t分别为第i台发电机在t时段的发电成本和输出功率;NLA为A类可响应负荷个数;和C-LA,j,t分别为第j个A类负荷在t时段响应增加负荷信号和减少负荷信号的成本;和S-LA,j,t分别为第j个A类负荷在t时段是否响应增加负荷信号和减少负荷信号的状态,表示第j个A类负荷在t时段参与响应增加负荷信号,表示第j个A类负荷在t时段未响应增加负荷信号,S-LA,j,t=1表示第j个A类负荷在t时段参与响应减少负荷信号,S-LA,j,t=0表示第j个A类负荷在t时段未响应减少负荷信号;PLA,j为第j个A类负荷的额定功率;NW为风电场个数;CW,k,t为第k个风电场在t时段的弃风成本;PW,ahead,k,t为第k个风电场在t时段的日前风电预测功率;PW,k,t为第k个风电场在t时段的风电输出功率;
日前调度模型的约束条件为:
①有功功率平衡约束:
其中,NL为除A类负荷以外的负荷个数;Pl,t为第l个负荷的额定功率;
②A类负荷响应平衡约束:
③线路潮流约束:
Pmn,t=Bmn(θn,t-θm,t)
其中,Pmn,t为线路mn在t时段的传输功率;Bmn为线路mn的电纳;θm,t和θn,t分别为t时段线路mn中m节点和n节点的相角;
④发电机出力上下限约束:
PG,i,min≤PG,i,t≤PG,i,max
其中,PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值;
⑤发电机爬坡约束:
-Rd,iΔT1≤PG,i,t-PG,i,t-1≤Ru,iΔT1
其中,PG,i,t和PG,i,t-1分别为第i台发电机在t时段和t-1时段的输出功率;Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT1为每个t时段的时长;
⑥安全约束:
|Pmn,t|≤Pmn,lim
-π≤θt≤π
其中,Pmn,lim为线路mn的传输极限;θt为各节点的相角;
⑦风电出力上下限约束:
0≤PW,k,t≤PW,ahead,k,t
⑧A类负荷响应约束:
其中,DRPahead,u和DRPahead,d分别为A类负荷响应增加负荷信号和减少负荷信号的潜力最大值。
作为优选,所述日前调度模型确定的决策变量包括各发电机组日前调度计划中的发电功率PG,i,t、风电出力PW,k,t和A类负荷的调整量,t时段的调整量为其中常规火电机组的发电功率和A类负荷的调整量用于日内调度模型的输入。
作为优选,所述步骤(2)中日内调度在基于日前调度的第一时间间隔划分的时间段内,以第二时间间隔计算日内调度的优化决策变量,日内调度模型具体如下:
日内调度模型的目标函数为:
其中,为日内调度中第i台发电机出力变化成本;为第i台发电机日内调度输出功率;PG,i为日前调度的第i台发电机输出功率;NLB为B类负荷个数;和C-LB,j分别为第j个B类负荷响应增加负荷信号和减少负荷信号的成本;和S-LB,j分别为第j个B类负荷是否响应增加负荷信号和减少负荷信号的状态,表示第j个B类负荷参与响应增加负荷信号,表示第j个B类负荷未响应增加负荷信号,S-LB,j=1表示第j个B类负荷参与响应减少负荷信号,S-LB,j=0表示第j个B类负荷未响应减少负荷信号;PLB,j为第j个B类负荷的额定功率;CW,k为第k个风电场的弃风成本;PW,Nei,k为基于第一时间间隔划分的时间段内日内风电预测功率;为第k个风电场日内调度输出功率;
日内调度的约束条件为:
①有功功率平衡约束:
其中,NL1为除B类负荷以外的负荷个数;Pl为第l个负荷的额定功率;δBj为第j个未参与响应的B类负荷运行状态,δBj=1表示第j个B类负荷正在工作,δBj=0表示第j个B类负荷停止工作;
②线路潮流约束:
Pmn=Bmn(θn-θm)
其中,Pmn为线路mn的传输功率;Bmn为线路mn的电纳;θm和θn分别为线路mn中m节点和n节点的相角;
③发电机出力上下限约束:
其中,PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值;
④发电机爬坡约束:
其中,Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT2为第二时间间隔;
⑤安全约束:
|Pmn|≤Pmn,lim
-π≤θ≤π
其中,Pmn,lim为线路mn的传输极限;θ为各节点的相角;
⑥风电出力上下限约束:
⑦B类负荷响应约束:
其中,DRPNei,u和DRPNei,d分别为B类负荷响应增加负荷信号和减少负荷信号的潜力最大值。
作为优选,所述步骤(3)中实时调度在第二时间间隔范围内以第三时间间隔计算实时调度的优化决策变量,实时调度模型具体如下:
实时调度模型的目标函数为:
其中,为实时调度中第i台发电机出力变化成本;为第i台发电机实时输出功率;和分别为第j个C类负荷响应增加负荷信号和减少负荷信号的成本;和分别为第j个C类负荷是否响应增加负荷信号和减少负荷信号的状态,表示第j个C类负荷参与响应增加负荷信号,表示第j个C类负荷未响应增加负荷信号,表示第j个C类负荷参与响应减少负荷信号,表示第j个C类负荷未响应减少负荷信号;PLC,j为第j个C类负荷的额定功率;PW,real,k为实时风电预测功率;PrealW,k为第k个风电场实时输出功率;
实时调度的约束条件为:
①有功功率平衡约束:
其中,NL2为除C类负荷以外的负荷个数;Pl为第l个负荷的额定功率;δCj为第j个未参与响应的C类负荷运行状态,δCj=1表示第j个C类负荷正在工作,δCj=0表示第j个C类负荷停止工作;
②线路潮流约束:
Pmn=Bmn(θn-θm)
其中,Pmn为线路mn的传输功率;Bmn为线路mn的电纳;θm和θn分别为线路mn中m节点和n节点的相角;
③发电机出力上下限约束:
其中,PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值;
④发电机爬坡约束:
其中,Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT3为第三时间间隔;
⑤安全约束:
|Pmn|≤Pmn,lim
-π≤θt≤π
其中,Pmn,lim为线路mn的传输极限;θ为各节点的相角;
⑥风电出力上下限约束:
⑦C类负荷响应约束:
其中,
其中,DRPreal,u和DRPreal,d分别为C类负荷响应增加负荷信号和减少负荷信号的潜力最大值;N1为空调个数;为第a个空调的额定功率(kW);为第a个空调在t时刻响应增加负荷信号或者减少负荷信号的DR潜力状态;N2为热水器个数;为第h个热水器的额定功率(kW);为第h个热水器在t时刻响应增加负荷信号或者减少负荷信号的DR潜力状态;N3为电动汽车的个数;为第e个电动汽车的额定功率(kW);为第e个电动汽车在t时刻响应增加负荷信号或者减少负荷信号DR潜力状态;Re(t)为根据历史数据预估的t时刻智能家电反悔度。
作为优选,所述第一时间间隔为24h。
作为优选,所述第二时间间隔为15min。
作为优选,所述第三时间间隔为1min。
有益效果:本发明在高风电渗透率和高预测误差情况下,充分利用智能家电参与响应,利用智能家电的需求响应潜力,分析多时间尺度的需求响应过程中用户行为的时变性和作用,建立网络约束条件下基于多时间尺度需求响应的源网荷经济调度模型,实现日前、日内和实时调度。与现有日前-实时调度模型进行对比,利用本发明所提模型和方法能够得到更低的调度成本和显著更低的调度响应量。
附图说明
图1是本发明考虑需求响应的多时间尺度源网荷协调调度方法的整体技术框架图;
图2是本发明考虑需求响应的多时间尺度源网荷协调调度方法的流程图;
图3是实现算例的***仿真线路图;
图4是风电不同时间尺度下预测功率曲线。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
本发明提供了一种考虑需求响应的多时间尺度源网荷协调调度方法,在一个实施例中,根据需求响应和调度的时间特性,将整体过程分为日前24h调度、日内15min调度和实时1min调度三个阶段,建立了网络约束条件下基于多时间尺度需求响应的源网荷经济调度模型,包括日前调度、日内调度和实时调度模型。但是本领域技术人员应当理解,这里划分的24h、15min和1min的时间尺度仅仅是用于说明本发明,而不是对本发明的限制。
如图1所示是根据该实施例的考虑需求响应的多时间尺度源网荷协调调度方法的整体技术框架图,每个时间尺度的调度分为输入层、调度控制层、代理协调层和本地响应层。输入层将风电场和负荷侧常规负荷(不包括参与DR的负荷)的预测功率作为输入参与各阶段的调度,其中Pwfa和Plfa、Pwfn和Plfn、Pwfr和Plfr分别为日前调度、日内调度、实时调度的风电场和负荷侧常规负荷的预测功率。调度控制层负责制定并实施调度计划,日前调度计划为24h执行一次,分辨率为1h,调度任务包括确定发电机组出力和需要提前一天通知的负荷所在各负荷聚合商的响应量,需要提前一天通知的负荷下文中称为A类负荷;日内调度计划为每15min执行一次,分辨率为15min,调度任务包括确定发电机组出力变化量和提前1h通知的负荷所在各负荷聚合商的响应量,需要提前1h通知的负荷在下文中称为B类负荷;实时调度计划为每1min执行一次,分辨率为1min,调度任务包括发电机组出力变化量和能够实时参与调度的负荷所在各负荷聚合商的响应量,能够实时参与调度的负荷在下文中称为C类负荷。日前调度计划每天24:00制定一次,与此同时每15min滚动制定一次日内调度计划,每1min滚动制定一次实时调度计划。随着时间的推移,日内、实时调度计划对应的时段不断向前推移。也就是说,日前调度是提前24h制定接下来24h的计划,1h为分辨率,即未来24个1h的出力,包括常规发电机组的出力、风电场出力和A类负荷的调整量。而日内调度计划就是每15min制定接下来15min的计划,实时是1min执行一次,日内计划其实是对日前计划的修正,常规发电机组出力会发生改变,B类负荷参与调整,而风电场出力由于时间尺度不同导致预测数据不一样,使得出力也会不一样,实时调度也是如此。PGa和PLa、PGn和PLn、PGr和PLr分别为日前调度、日内调度、实时调度的发电机和负荷的输出结果。代理协调层协调***侧调度信息和负荷侧响应资源,做出针对某一优化目标的最优决策,向参与响应的负荷下发控制信号。在实时调度中,代理协调层的各负荷聚合商向调度控制层上传负荷实时聚合需求响应潜力DLr。本地响应层向各负荷聚合商上传参与响应的负荷用电信息。
如图2所示是考虑需求响应的多时间尺度源网荷协调调度方法的流程图。本发明的方法具体包括以下步骤:
步骤(1)、已知单个风电场预测出力、负荷预测曲线和A类负荷的响应潜力,以发电机成本、弃风成本和负荷响应成本总和最小为目标,建立考虑需求响应的日前调度模型,确定日前调度计划中常规发电机组出力、风电场出力和A类负荷的调整量;日前调度计划24h执行一次,分辨率为1h,即包括24个时间段。
所述日前调度模型为:
1)日前调度计划模型的目标函数:
式中:NT为调度时段个数;NG为发电机个数;CG,i,t和PG,i,t分别为第i台发电机在t时段的发电成本和输出功率;NLA为A类可响应负荷个数;和C-LA,j,t分别为第j个A类负荷在t时段响应增加负荷信号和减少负荷信号的成本;和S-LA,j,t分别为第j个A类负荷在t时段是否响应增加负荷信号和减少负荷信号的状态,表示第j个A类负荷在t时段参与响应增加负荷信号,表示第j个A类负荷在t时段未响应增加负荷信号,S-LA,j,t=1表示第j个A类负荷在t时段参与响应减少负荷信号,S-LA,j,t=0表示第j个A类负荷在t时段未响应减少负荷信号;PLA,j为第j个A类负荷的额定功率;NW为风电场个数;CW,k,t为第k个风电场在t时段的弃风成本;PW,ahead,k,t为第k个风电场在t时段的日前风电预测功率;PW,k,t为第k个风电场在t时段的风电输出功率。
2)约束条件:
①有功功率平衡约束:
式中:NL为除A类负荷以外的负荷个数;Pl,t为第l个负荷的额定功率。
②A类负荷响应平衡约束:
③线路潮流约束:
Pmn,t=Bmn(θn,t-θm,t)
式中:Pmn,t为线路mn在t时段的传输功率;Bmn为线路mn的电纳;θm,t和θn,t分别为t时段线路mn中m节点和n节点的相角。
④发电机出力上下限约束:
PG,i,min≤PG,i,t≤PG,i,max
式中:PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值。
⑤发电机爬坡约束:
-Rd,iΔT1≤PG,i,t-PG,i,t-1≤Ru,iΔT1
式中:PG,i,t和PG,i,t-1分别为第i台发电机在t时段和t-1时段的输出功率;Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT1为t时段时长,即1h。
⑥安全约束:
|Pmn,t|≤Pmn,lim
-π≤θt≤π
式中:Pmn,lim为线路mn的传输极限;θt为各节点的相角。
⑦风电出力上下限约束:
0≤PW,k,t≤PW,ahead,k,t
⑧A类负荷响应约束:
式中:DRPahead,u和DRPahead,d分别为A类负荷响应增加负荷信号和减少负荷信号的潜力最大值。
上述1)和2)决定的日前调度计划模型确定的决策变量包括各发电机组日前调度计划中的发电功率、风电出力和A类负荷的调整量将常规火电机组的发电功率和A类负荷的调整量带入日内调度模型并作为基准值求解后续的优化模型。
步骤(2)、将日前调度的24个时间点优化决策变量作为日内调度的已知量,输入风电场的日内预测曲线和B类负荷响应潜力,以总成本最小建立考虑需求响应的日内调度模型,进而确定发电机出力变化量、风电场出力和B类负荷调整量。
所述日内调度模型为:
1)日前调度计划模型的目标函数:
式中:为日内调度中第i台发电机出力变化成本;为第i台发电机日内调度输出功率;PG,i为日前调度的第i台发电机输出功率;NLB为B类负荷个数;和C-LB,j分别为第j个B类负荷响应增加负荷信号和减少负荷信号的成本;和S-LB,j分别为第j个B类负荷是否响应增加负荷信号和减少负荷信号的状态,表示第j个B类负荷参与响应增加负荷信号,表示第j个B类负荷未响应增加负荷信号,S-LB,j=1表示第j个B类负荷参与响应减少负荷信号,S-LB,j=0表示第j个B类负荷未响应减少负荷信号;PLB,j为第j个B类负荷的额定功率;CW,k为第k个风电场的弃风成本;PW,Nei,k为日内1h风电预测功率;为第k个风电场日内调度输出功率。
2)约束条件
①有功功率平衡约束:
式中:NL1为除B类负荷以外的负荷个数;Pl为第l个负荷的额定功率;δBj为第j个未参与响应的B类负荷运行状态,δBj=1表示第j个B类负荷正在工作,δBj=0表示第j个B类负荷停止工作。
②线路潮流约束:
Pmn=Bmn(θn-θm)
式中:Pmn为线路mn的传输功率;Bmn为线路mn的电纳;θm和θn分别为线路mn中m节点和n节点的相角。
③发电机出力上下限约束:
式中:PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值。
④发电机爬坡约束:
式中:Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT2为第二时间间隔,即15min。
⑤安全约束:
|Pmn|≤Pmn,lim
-π≤θ≤π
式中:Pmn,lim为线路mn的传输极限;θ为各节点的相角。
⑥风电出力上下限约束:
⑦B类负荷响应约束:
式中:DRPNei,u和DRPNei,d分别为B类负荷响应增加负荷信号和减少负荷信号的潜力最大值。
步骤(3)、日内调度的优化决策变量作为实时调度的已知量,输入风电场的实时预测曲线和当前时刻C类负荷响应潜力,以总成本最小建立考虑需求响应的实时调度模型,进而确定发电机出力变化量、风电场出力和C类负荷调整量。在本发明的模型中,实时参与调度的C类负荷主要包括智能家电,可根据C类负荷调整量的大小实时动作。考虑到用户存在反悔的可能性,也就是参与DR计划的智能家电未能及时动作,多次响应至满足***负荷调整量的需求。多次响应的意思是调度策略的多次实施,先对某智能家电下发指令,看其是否响应状态反馈,如果不能响应则对其他的智能家电继续下发指令,直至满足***负荷调整量的需求。C类负荷响应动作之后,更新下一阶段的C类负荷总功率,对下一时间段进行实时调度。
所述实时调度模型为:
1)实时调度计划模型的目标函数:
式中:为实时调度中第i台发电机出力变化成本;为第i台发电机实时输出功率;和分别为第j个C类负荷响应增加负荷信号和减少负荷信号的成本;和分别为第j个C类负荷是否响应增加负荷信号和减少负荷信号的状态,表示第j个C类负荷参与响应增加负荷信号,表示第j个C类负荷未响应增加负荷信号,表示第j个C类负荷参与响应减少负荷信号,表示第j个C类负荷未响应减少负荷信号;PLC,j为第j个C类负荷的额定功率;PW,real,k为实时风电预测功率;PrealW,k为第k个风电场实时输出功率。
2)约束条件
①有功功率平衡约束:
式中:NL2为除C类负荷以外的负荷个数;Pl为第l个负荷的额定功率;δCj为第j个未参与响应的C类负荷运行状态,δCj=1表示第j个C类负荷正在工作,δCj=0表示第j个C类负荷停止工作。
②线路潮流约束:
Pmn=Bmn(θn-θm)
式中:Pmn,t为线路mn的传输功率;Bmn为线路mn的电纳;θm和θn分别为线路mn中m节点和n节点的相角。
③发电机出力上下限约束:
④发电机爬坡约束:
⑤安全约束:
|Pmn|≤Pmn,lim
-π≤θt≤π
式中:Pmn,lim为线路mn的传输极限;θ为各节点的相角。
⑥风电出力上下限约束:
⑦C类负荷响应约束:
其中:
其中,DRPreal,u和DRPreal,d分别为C类负荷响应增加负荷信号和减少负荷信号的潜力最大值;N1为空调个数;为第a个空调的额定功率(kW);为第a个空调在t时刻响应增加负荷信号或者减少负荷信号的DR潜力状态;N2为热水器个数;为第h个热水器的额定功率(kW);为第h个热水器在t时刻响应增加负荷信号或者减少负荷信号的DR潜力状态;N3为电动汽车的个数;为第e个电动汽车的额定功率(kW);为第e个电动汽车在t时刻响应增加负荷信号或者减少负荷信号DR潜力状态;Re(t)为根据历史数据预估的t时刻智能家电反悔度。
图3所示为根据一个实现算例的仿真线路图,该***包括12台发电机和17个负荷节点,其中风电场位于节点19。设置10个负荷聚合商,每个聚合商均包含A类、B类和C类负荷,负荷聚合商所在节点分别为节点3、4、5、6、8、9、10、14、19和20。各发电机所在节点位置和发电成本如表1所示,弃风成本为21$/MWh,A类、B类和C类负荷的响应成本分别为9.87$/MWh、12$/MWh和14$/MWh。
表1
假设***中可调用的A类和B类负荷的容量不大于总负荷的5%,***中可调用的C类负荷容量随着时间会发生改变,其额定功率总和不大于总负荷的25%。总负荷各时间段的功率如表2所示。为阐述本发明模型影响因素的应用效果,***中各线路传输功率极限减少至所允许容量的50%。模型采用MATLAB软件进行求解。日前24h、日内15min和实时***风电出力预测曲线如图4所示。
表2
为了验证本发明提出的日前-日内-实时多级调度模型的有效性,设置2个仿真条件:条件一为考虑日前-日内-实时调度配合的调度策略;条件二为考虑日前-实时调度配合的策略。相关仿真结果如表3所示。
表3
由表3的仿真结果可见,两种调度测量下最终都能实现风电消纳率为1。在条件一所描述的调度策略下,虽然包括了日前、日内和实时三个阶段的成本,但是总调度成本仍然比条件二调度策略下的成本低,并且负荷响应量显著降降低。
Claims (9)
1.一种考虑需求响应的多时间尺度源网荷协调调度方法,其特征在于,所述方法根据需求响应和调度的时间特性,将整体调度过程分为日前调度、日内调度和实时调度三个阶段,建立网络约束条件下基于多时间尺度需求响应的源网荷经济调度模型,并基于调度模型的求解结果制定调度计划,具体包括以下步骤:
(1)在已知单个风电场预测出力、负荷预测曲线和A类负荷响应潜力的情况下,以发电机成本、弃风成本和负荷响应成本总和最小为目标,建立考虑需求响应的日前调度模型,确定日前调度计划中常规发电机组出力、风电场出力和A类负荷的调整量,以第一时间间隔制定日前调度计划;
(2)将日前调度的优化决策变量中的发电机组出力和A类负荷的调整量作为日内调度的已知量,输入风电场的日内预测曲线和B类负荷响应潜力,以总成本最小为目标建立考虑需求响应的日内调度模型,进而确定日内调度计划中发电机组出力、风电场出力和B类负荷调整量,以第二时间间隔制定日内调度计划;
(3)将日内调度的优化决策变量中的发电机组出力和B类负荷调整量作为实时调度的已知量,输入风电场的实时预测曲线和当前时刻C类负荷响应潜力,以总成本最小为目标建立考虑需求响应的实时调度模型,进而确定实时调度计划中发电机组出力、风电场出力和C类负荷调整量,以第三时间间隔制定实时调度计划;
其中,第一时间间隔>第二时间间隔>第三时间间隔;
(4)调度中心向参与DR计划的智能家电下发状态指令,智能家电进行多次响应至满足***C类负荷调整量的需求;C类负荷响应动作之后,更新下一阶段的C类负荷总功率,对下一时间段进行实时调度;
其中,所述步骤(1)中日前调度将第一时间间隔划分为若干个时段,计算每个时段的日前优化决策量,日前调度模型具体如下:
日前调度模型的目标函数为:
其中,NT为调度时段个数;NG为发电机个数;CG,i,t和PG,i,t分别为第i台发电机在t时段的发电成本和输出功率;NLA为A类可响应负荷个数;和分别为第j个A类负荷在t时段响应增加负荷信号和减少负荷信号的成本;和分别为第j个A类负荷在t时段是否响应增加负荷信号和减少负荷信号的状态,表示第j个A类负荷在t时段参与响应增加负荷信号,表示第j个A类负荷在t时段未响应增加负荷信号,表示第j个A类负荷在t时段参与响应减少负荷信号,表示第j个A类负荷在t时段未响应减少负荷信号;PLA,j为第j个A类负荷的额定功率;NW为风电场个数;CW,k,t为第k个风电场在t时段的弃风成本;PW,ahead,k,t为第k个风电场在t时段的日前风电预测功率;PW,k,t为第k个风电场在t时段的风电输出功率;
日前调度模型的约束条件为:
①有功功率平衡约束:
其中,NL为除A类负荷以外的负荷个数;Pl,t为第l个负荷的额定功率;
②A类负荷响应平衡约束:
③线路潮流约束:
Pmn,t=Bmn(θn,t-θm,t)
其中,Pmn,t为线路mn在t时段的传输功率;Bmn为线路mn的电纳;θm,t和θn,t分别为t时段线路mn中m节点和n节点的相角;
④发电机出力上下限约束:
PG,i,min≤PG,i,t≤PG,i,max
其中,PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值;
⑤发电机爬坡约束:
-Rd,iΔT1≤PG,i,t-PG,i,t-1≤Ru,iΔT1
其中,PG,i,t和PG,i,t-1分别为第i台发电机在t时段和t-1时段的输出功率;Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT1为每个t时段的时长;
⑥安全约束:
|Pmn,t|≤Pmn,lim
-π≤θt≤π
其中,Pmn,lim为线路mn的传输极限;θt为各节点的相角;
⑦风电出力上下限约束:
0≤PW,k,t≤PW,ahead,k,t
⑧A类负荷响应约束:
其中,DRPahead,u和DRPahead,d分别为A类负荷响应增加负荷信号和减少负荷信号的潜力最大值。
3.根据权利要求1所述的考虑需求响应的多时间尺度源网荷协调调度方法,其特征在于,所述步骤(2)中日内调度在基于日前调度的第一时间间隔划分的时间段内,以第二时间间隔计算日内调度的优化决策变量,日内调度模型具体如下:
日内调度模型的目标函数为:
其中,为日内调度中第i台发电机出力变化成本;为第i台发电机日内调度输出功率;PG,i为日前调度的第i台发电机输出功率;NLB为B类负荷个数;和分别为第j个B类负荷响应增加负荷信号和减少负荷信号的成本;和分别为第j个B类负荷是否响应增加负荷信号和减少负荷信号的状态,表示第j个B类负荷参与响应增加负荷信号,表示第j个B类负荷未响应增加负荷信号,表示第j个B类负荷参与响应减少负荷信号,表示第j个B类负荷未响应减少负荷信号;PLB,j为第j个B类负荷的额定功率;CW,k为第k个风电场的弃风成本;PW,Nei,k为基于第一时间间隔划分的时间段内日内风电预测功率;为第k个风电场日内调度输出功率;
日内调度的约束条件为:
①有功功率平衡约束:
其中,NL1为除B类负荷以外的负荷个数;Pl为第l个负荷的额定功率;δBj为第j个未参与响应的B类负荷运行状态,δBj=1表示第j个B类负荷正在工作,δBj=0表示第j个B类负荷停止工作;
②线路潮流约束:
Pmn=Bmn(θn-θm)
其中,Pmn为线路mn的传输功率;Bmn为线路mn的电纳;θm和θn分别为线路mn中m节点和n节点的相角;
③发电机出力上下限约束:
其中,PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值;
④发电机爬坡约束:
其中,Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT2为第二时间间隔;
⑤安全约束:
|Pmn|≤Pmn,lim
-π≤θ≤π
其中,Pmn,lim为线路mn的传输极限;θ为各节点的相角;
⑥风电出力上下限约束:
⑦B类负荷响应约束:
其中,DRPNei,u和DRPNei,d分别为B类负荷响应增加负荷信号和减少负荷信号的潜力最大值。
5.根据权利要求1所述的考虑需求响应的多时间尺度源网荷协调调度方法,其特征在于,所述步骤(3)中实时调度在第二时间间隔范围内以第三时间间隔计算实时调度的优化决策变量,实时调度模型具体如下:
实时调度模型的目标函数为:
其中,NLC为C类负荷个数;为实时调度中第i台发电机出力变化成本;为第i台发电机实时输出功率;和分别为第j个C类负荷响应增加负荷信号和减少负荷信号的成本;和分别为第j个C类负荷是否响应增加负荷信号和减少负荷信号的状态,表示第j个C类负荷参与响应增加负荷信号,表示第j个C类负荷未响应增加负荷信号,表示第j个C类负荷参与响应减少负荷信号,表示第j个C类负荷未响应减少负荷信号;PLC,j为第j个C类负荷的额定功率;PW,real,k为实时风电预测功率;PrealW,k为第k个风电场实时输出功率;
实时调度的约束条件为:
①有功功率平衡约束:
其中,NL2为除C类负荷以外的负荷个数;Pl为第l个负荷的额定功率;δCj为第j个未参与响应的C类负荷运行状态,δCj=1表示第j个C类负荷正在工作,δCj=0表示第j个C类负荷停止工作;
②线路潮流约束:
Pmn=Bmn(θn-θm)
其中,Pmn为线路mn的传输功率;Bmn为线路mn的电纳;θm和θn分别为线路mn中m节点和n节点的相角;
③发电机出力上下限约束:
其中,PG,i,min和PG,i,max分别为第i台发电机的出力最小值和最大值;
④发电机爬坡约束:
其中,Rd,i和Ru,i分别为第i台发电机单位时间内可以降低和升高的发电出力,即爬坡率;ΔT3为第三时间间隔;
⑤安全约束:
|Pmn|≤Pmn,lim
-π≤θt≤π
其中,Pmn,lim为线路mn的传输极限;θt为各节点的相角;
⑥风电出力上下限约束:
⑦C类负荷响应约束:
其中,
7.根据权利要求1所述的考虑需求响应的多时间尺度源网荷协调调度方法,其特征在于,所述第一时间间隔为24h。
8.根据权利要求1所述的考虑需求响应的多时间尺度源网荷协调调度方法,其特征在于,所述第二时间间隔为15min。
9.根据权利要求1所述的考虑需求响应的多时间尺度源网荷协调调度方法,其特征在于,所述第三时间间隔为1min。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811128233.4A CN109301817B (zh) | 2018-09-27 | 2018-09-27 | 一种考虑需求响应的多时间尺度源网荷协调调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811128233.4A CN109301817B (zh) | 2018-09-27 | 2018-09-27 | 一种考虑需求响应的多时间尺度源网荷协调调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109301817A CN109301817A (zh) | 2019-02-01 |
CN109301817B true CN109301817B (zh) | 2020-09-18 |
Family
ID=65164526
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811128233.4A Active CN109301817B (zh) | 2018-09-27 | 2018-09-27 | 一种考虑需求响应的多时间尺度源网荷协调调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109301817B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112332460B (zh) * | 2020-10-30 | 2024-06-04 | 重庆大学 | 一种考虑能流特性差异的电-气互联***异步调度方法 |
CN113610339A (zh) * | 2021-06-24 | 2021-11-05 | 华北电力大学 | 一种基于多时间尺度协调的需求响应*** |
CN113420728B (zh) * | 2021-08-23 | 2021-11-09 | 国网江苏省电力有限公司营销服务中心 | 融合多时间尺度信息的非侵入式空调负荷辨识方法及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194514A (zh) * | 2017-05-27 | 2017-09-22 | 重庆大学 | 一种针对风电功率预测误差的需求响应多时间尺度调度方法 |
CN107732977A (zh) * | 2017-09-21 | 2018-02-23 | 东南大学 | 一种基于需求响应的agc实时调度方法 |
CN108288854A (zh) * | 2018-01-04 | 2018-07-17 | 国网江苏省电力有限公司电力科学研究院 | 一种源网荷协调控制方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108416536A (zh) * | 2018-04-10 | 2018-08-17 | 国网江苏省电力有限公司电力科学研究院 | 一种消纳新能源的需求响应资源多时间尺度滚动调度方法 |
-
2018
- 2018-09-27 CN CN201811128233.4A patent/CN109301817B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107194514A (zh) * | 2017-05-27 | 2017-09-22 | 重庆大学 | 一种针对风电功率预测误差的需求响应多时间尺度调度方法 |
CN107732977A (zh) * | 2017-09-21 | 2018-02-23 | 东南大学 | 一种基于需求响应的agc实时调度方法 |
CN108288854A (zh) * | 2018-01-04 | 2018-07-17 | 国网江苏省电力有限公司电力科学研究院 | 一种源网荷协调控制方法及*** |
Non-Patent Citations (1)
Title |
---|
考虑大规模风电接入并计及多时间尺度需求响应资源协调优化的滚动调度模型;包宇庆 等;《中国电机工程学报》;20160905;第4589-4599页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109301817A (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hui et al. | A transactive energy framework for inverter-based HVAC loads in a real-time local electricity market considering distributed energy resources | |
CN109301817B (zh) | 一种考虑需求响应的多时间尺度源网荷协调调度方法 | |
CN104134994B (zh) | 一种利用agc进行稳定断面潮流越限的校正控制方法 | |
CN107069776B (zh) | 一种平滑微网联络线功率的储能前瞻分布式控制方法 | |
CN107732977B (zh) | 一种基于需求响应的agc实时调度方法 | |
Hui et al. | Coordination control of distributed generators and load resources for frequency restoration in isolated urban microgrids | |
CN104915737A (zh) | 基于风险控制含大规模风电的电力***协调经济调度方法 | |
Saez-de-Ibarra et al. | Management strategy for market participation of photovoltaic power plants including storage systems | |
CN108092324A (zh) | 一种风电参与调峰调频的agc控制***和控制方法 | |
CN113725913B (zh) | 一种适应源网荷储多类型资源协同的自动功率控制方法 | |
CN106160000B (zh) | 一种分调-省调两级协调的风电增发控制方法 | |
CN105006846A (zh) | 一种风电场场站级有功功率优化方法 | |
CN107346893B (zh) | 一种风电场有功功率控制方法及*** | |
CN106602613B (zh) | 一种省地两级调度断面协调控制方法 | |
CN103151802A (zh) | 多时间尺度的主动配电网dg协调控制***及方法 | |
CN103259274A (zh) | 多区域互联电网特高压联络线功率控制方法及控制器 | |
CN105743089A (zh) | 基于柔性负荷的需求侧负荷调控方法及*** | |
Zhao et al. | Model predictive based frequency control of power system incorporating air-conditioning loads with communication delay | |
CN104239966B (zh) | 一种主动配电网基于用电成本差异化的运营方法 | |
CN106300417A (zh) | 基于模型预测控制的风电场群无功电压优化控制方法 | |
CN106340893A (zh) | 一种基于虚拟发电厂的分布式光伏发电功率控制方法 | |
Lu et al. | Two-stage robust scheduling and real-time load control of community microgrid with multiple uncertainties | |
CN104182808A (zh) | 一种基于等比例限电的新能源场站发电计划制定方法 | |
CN105490268A (zh) | 一种交直流互联电网的负荷跟踪方法及*** | |
CN104578174A (zh) | 一种基于指令的主站agc机组控制方法 |
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 |