CN107732983A - 一种电‑热耦合多能流系统的日内滚动调度方法 - Google Patents
一种电‑热耦合多能流系统的日内滚动调度方法 Download PDFInfo
- Publication number
- CN107732983A CN107732983A CN201710989015.9A CN201710989015A CN107732983A CN 107732983 A CN107732983 A CN 107732983A CN 201710989015 A CN201710989015 A CN 201710989015A CN 107732983 A CN107732983 A CN 107732983A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- electric
- thermal
- streaming system
- 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.)
- Pending
Links
- 230000008878 coupling Effects 0.000 title claims abstract description 102
- 238000010168 coupling process Methods 0.000 title claims abstract description 102
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 102
- 238000005096 rolling process Methods 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 24
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000010977 unit operation Methods 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000005619 thermoelectricity Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000005265 energy consumption Methods 0.000 abstract description 2
- 238000005485 electric heating Methods 0.000 abstract 1
- 238000010438 heat treatment Methods 0.000 abstract 1
- 238000007726 management method Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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
- 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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F24—HEATING; RANGES; VENTILATING
- F24D—DOMESTIC- OR SPACE-HEATING SYSTEMS, e.g. CENTRAL HEATING SYSTEMS; DOMESTIC HOT-WATER SUPPLY SYSTEMS; ELEMENTS OR COMPONENTS THEREFOR
- F24D19/00—Details
- F24D19/10—Arrangement or mounting of control or safety devices
- F24D19/1006—Arrangement or mounting of control or safety devices for water heating systems
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
-
- 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
-
- 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
-
- 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
-
- 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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- 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
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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
- 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
- Y02P80/15—On-site combined power, heat or cool generation or distribution, e.g. combined heat and power [CHP] 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
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Water Supply & Treatment (AREA)
- Chemical & Material Sciences (AREA)
- Public Health (AREA)
- Power Engineering (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Thermal Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- Automation & Control Theory (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明涉及一种电‐热耦合多能流系统的日内滚动调度方法,属于含多种能源形式的电网运行和控制技术领域。本方法考虑电‑热系统的相互影响,实现了电‑热耦合多能流系统的日内滚动调度。本方法可以不断对日前调度计划进行修正,得到更加准确的日内滚动调度方案。相比独立地对供电、供热系统进行优化调度分析,不仅能得到更优的调度方案(总运行成本更低),还提高了调度的灵活性。该方法可以应用于电‑热耦合多能流系统的日内滚动调度计划制定,有利于提高电‑热耦合多能流系统的用能效率,减少运行成本。
Description
技术领域
本发明涉及一种电‐热耦合多能流系统的日内滚动调度方法,属于含多种能源形式的电网运行和控制技术领域。
背景技术
能源综合利用是提高综合能源利用效率、促进可再生能源消纳的重要途径,通过打破原来电、热、冷、气、交通等能流子系统相对割裂的状态,实现多类型能源开放互联,构建多能流系统。多能流是指多种类型的能量流,表示电、热、冷、气、交通等能量流的相互耦合、转换和传输。多能流系统相比传统相互割裂的能源系统,其带来的效益包括:1)通过多类型能源的梯级开发利用和智能管理,可以降低能源消耗和浪费,提高综合能源利用效率,并有助于减少总的用能成本;2)利用不同能源的特性差异和互补、转换,有助于提高消纳间歇式可再生能源的能力;3)通过多能源的转供、互补和协调控制,有助于提高供能的可靠性,并为电网的运行提供更多可调控资源;4)通过多能流系统的协同规划和建设,可以减少基础设施的重复建设和浪费,提高资产利用率。
多能流系统一方面具有可观的效益,另一方面也使原本复杂的能源系统更加复杂。多能流系统由多个能流子系统组成,这些能流子系统之间相互作用和影响,使得多能流系统复杂度显著增加,体现出许多新的特性,传统各个能流单独分析的方法已经难以适应新的要求,亟需发展出新的多能流分析方法。在我国,越来越多的热电联产机组、热泵、电锅炉等耦合元件客观上增强了电-热之间的互联,促进了电-热耦合多能流系统的发展,也对电-热耦合多能流系统的运行和控制技术提出了新的要求。
负荷预测是调度机构每日的重要工作,有效的负荷预测可以保证供能设备出力的合理调度。但是,日前负荷预测相比实际负荷可能存在比较大的偏差,日前制定的调度计划可能无法满足用户的用能需求。日内滚动调度以15分钟为一个步长进行超短期负荷预测,滚动修正供能设备出力的出力计划,可以得到更加准确的日内滚动调度方案。目前日内滚动调度方面的研究主要集中在单个独立的系统,为了使得电-热耦合多能流系统的运行成本最低,,需要研究电-热耦合多能流系统的日内滚动调度方法。
发明内容
本发明的目的是提出一种电‐热耦合多能流系统的日内滚动调度方法,以弥补现有领域研究的空白,建立电-热耦合多能流系统的日内滚动调度模型,实现电-热耦合多能流系统的优化调度。
本发明提出的电‐热耦合多能流系统的日内滚动调度方法,包括以下步骤:
(1)建立一个电-热耦合多能流系统优化调度的目标函数:
其中,pb,t为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的有功功率,qb,t为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的热功率,N为电-热耦合多能流系统中电‐热联供机组的总台数,F(pb,t,qb,t)为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的运行成本,px,t为电-热耦合多能流系统中第x台火电机组在第t个调度时段的有功功率,NTU为电-热耦合多能流系统的电网中火电机组的总台数,FTU(px,t)为电-热耦合多能流系统中电网的第x台火电机组在第t个调度时段的运行成本,Δt是相邻两个调度时刻的时间间隔,Δt的取值为15分钟;
(2)设定电-热耦合多能流系统中电网与热网稳态安全运行的等式约束条件,包括:
(2-1)电-热耦合多能流系统中的电网潮流方程如下:
其中,Pi,t为电-热耦合多能流系统的电网中节点i在第t个调度时段的注入有功功率,Qi,t为电网中节点i在第t个调度时段的注入无功功率,θi,t和θj,t分别为节点i和节点j在第t个调度时段的电压相角,Ui,t和Uj,t分别为节点i和节点j在第t个调度时段的电压幅值,Gij为电网节点导纳矩阵Y中第i行、第j列元素的实部,Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部,电网节点导纳矩阵Y从电-热耦合多能流系统的能量管理系统中获取,n为电网的节点总数;
(2-2)电‐热耦合多能流系统中热网的管道压力损失方程如下:
ΔHl,t=Slml,t|ml,t|,t=1,2,...,16,
其中,ΔHl,t为电‐热耦合多能流系统中热网的第l条管道在第t个调度时段的压力损失,Sl为第l条管道的阻力特性系数,Sl取值范围为10Pa/(kg/s)2≤Sl≤500Pa/(kg/s)2,ml,t为第l条管道在第t个调度时段的流量;
(2-3)电‐热耦合多能流系统中热网的循环泵水力特性方程如下:
其中,HP,t为循环泵在第t个调度时段的扬程,H0为循环泵静扬程,Sp为循环泵阻力系数,H0和Sp由循环泵的出厂说明书获取,mP,t为第t个调度时段流过循环泵的流量;
(2-4)电‐热耦合多能流系统中热网管道热量损失方程如下:
其中,Te,l,t为热网中第l条管道第t个调度时段的末端温度,Th,l,t为第l条管道第t个调度时段的首端温度,Ta,l,t为第l条管道在第t个调度时段的环境温度,ml,t为第l条管道在第t个调度时段的流量,Ll为第l条管道的长度,Cp为水的比热容,比热容的取值为4182焦耳/(千克·摄氏度),λ为管道单位长度的传热系数,λ从电-热耦合多能流系统的能量管理系统中获取,e为自然对数;
(2-5)电‐热耦合多能流系统的热网中多管道汇合点的温度方程:
其中,ml,t为第l条管道在第t个调度时段的流量,Te,l,t为热网中第l条管道第t个调度时段的末端温度,Th,l,t为第l条管道第t个调度时段的首端温度,为流出热网节点n的所有管道编号的集合,为流入热网节点n的所有管道编号的集合,QJ,n,t是第n个热网节点第t个调度时刻的热功率;
(2-6)通过电‐热联供机组耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,pb,t为第b台电-热联供机组在第t个调度时段的有功功率,q为第b台电-热联供机组在第t个调度时段的热功率,为第b台电-热联供机组运行可行域近似多边形的第k个顶点的横坐标,为第b台电-热联供机组运行可行域近似多边形的第k个顶点的纵坐标,为第b台电-热联供机组在第t个调度时段的第k个组合系数, NKb为第b台电-热联供机组的运行可行域近似多边形的顶点个数,电-热联供机组运行可行域近似多边形从电-热联供机组的出厂说明书中获取;
(2-7)通过循环泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,PP,t为循环泵第t个调度时段消耗的有功功率,g为重力加速度,ηP为循环泵效率,ηP的取值范围为0~1,mP,t为第t个调度时段的流过循环泵的流量,HP,t为循环泵第t个调度时段的的扬程;
(2-8)通过热泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
Php,t=ChpQhp,t,t=1,2,...,16
其中,Qhp,t为电‐热耦合多能流系统中第t个调度时段热泵发出的热功率,Php,t为第t个调度时段热泵消耗的电功率,Chp为热泵的产热效率,Chp从热泵的出厂说明书中获得;
(3)设定电-热耦合多能流系统中电网与热网稳态安全运行的不等式约束条件,包括:
(3-1)电-热耦合多能流系统的电网中第i个节点第t个调度时段的电压幅值Ui,t在设定的电网安全运行电压的上限值、下限值U i、之间运行,U i为第i个节点额定电压的0.95倍,为第i个节点额定电压的1.05倍:
(3-2)电-热耦合多能流系统的电网中第l条线路第t个调度时段的传输容量小于或等于设定的电网安全运行传输容量的最大值
(3-3)电-热耦合多能流系统的电网中电‐热联供机组有功功率的爬坡约束:
其中,和分别为第b台电-热联供机组有功功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,pb,t和pb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-4)电-热耦合多能流系统的电网中电‐热联供机组热功率的爬坡约束:
其中,和分别为第b台电-热联供机组热功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,qb,t和qb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-5)电-热耦合多能流系统的电网中火电机组有功功率的爬坡约束:
其中,和分别为第x台火电机组有功功率的向上和向下爬坡速率,和从火电机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,px,t和px,t-1分别为第x台火电机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-6)电-热耦合多能流系统的电网中第b台电-热联供机组第t个调度时段的有功功率pb,t在设定的电网安全运行第b台电-热联供机组有功功率的上、下限值 p b之间:
(3-7)电-热耦合多能流系统的电网中第b台电-热联供机组第t个调度时段的热功率qb,t在设定的电网安全运行第b台电-热联供机组热功率的上限值和下限值q b之间:
(3-8)电-热耦合多能流系统的电网中第x台火电机组第t个调度时段的有功功率px,t在设定的电网安全运行第x台火电机组有功功率的上限值p x和下限值之间:
(3-9)电-热耦合多能流系统的热网中第l条管道第t个调度时段的流量ml小于或等于热网安全运行流量的上限值
(3-10)电-热耦合多能流系统的热网中第i个换热站第t个调度时段回水温度Ti,t在设定的热网安全运行回水温度的上限值T和下限值之间:
(4)采用内点法,将上述步骤(1)中的方程作为目标函数,将上述步骤(2)和步骤(3)的所有方程作为约束条件,求解得到电-热耦合多能流系统中每台电‐热联供机组的有功功率和热功率,每台火电机组的有功功率,每台热泵的热功率,每台循环泵消耗的有功功率,作为电‐热耦合多能流系统的日内滚动调度方案,实现电‐热耦合多能流系统的日内滚动调度。
本发明提出的电-热耦合多能流系统的日内滚动调度方法,其特点和效果是:本方法考虑电-热系统的相互影响,实现了电-热耦合多能流系统的日内滚动调度。本方法可以不断对日前调度计划进行修正,得到更加准确的日内滚动调度方案。相比独立地对供电、供热系统进行优化调度分析,不仅能得到更优的调度方案(总运行成本更低),还提高了调度的灵活性。该方法可以应用于电-热耦合多能流系统的日内滚动调度计划制定,有利于提高电-热耦合多能流系统的用能效率,减少运行成本。
附图说明
图1是本发明方法涉及的电-热耦合多能流系统的结构示意图。
具体实施方式
本发明提出的电‐热耦合多能流系统的日内滚动调度方法,其中涉及的电-热耦合多能流系统的结构示意图如图1所示,该方法包括以下步骤:
(1)建立一个电-热耦合多能流系统优化调度的目标函数:
其中,pb,t为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的有功功率,qb,t为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的热功率,N为电-热耦合多能流系统中电‐热联供机组的总台数,F(pb,t,qb,t)为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的运行成本,px,t为电-热耦合多能流系统中第x台火电机组在第t个调度时段的有功功率,NTU为电-热耦合多能流系统的电网中火电机组的总台数,FTU(px,t)为电-热耦合多能流系统中电网的第x台火电机组在第t个调度时段的运行成本,Δt是相邻两个调度时刻的时间间隔,Δt的取值为15分钟;
(2)设定电-热耦合多能流系统中电网与热网稳态安全运行的等式约束条件,包括:
(2-1)电-热耦合多能流系统中的电网潮流方程如下:
其中,Pi,t为电-热耦合多能流系统的电网中节点i在第t个调度时段的注入有功功率,Qi,t为电网中节点i在第t个调度时段的注入无功功率,θi,t和θj,t分别为节点i和节点j在第t个调度时段的电压相角,Ui,t和Uj,t分别为节点i和节点j在第t个调度时段的电压幅值,Gij为电网节点导纳矩阵Y中第i行、第j列元素的实部,Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部,电网节点导纳矩阵Y从电-热耦合多能流系统的能量管理系统中获取,n为电网的节点总数;
(2-2)电‐热耦合多能流系统中热网的管道压力损失方程如下:
ΔHl,t=Slml,t|ml,t|,t=1,2,...,16,
其中,ΔHl,t为电‐热耦合多能流系统中热网的第l条管道在第t个调度时段的压力损失,Sl为第l条管道的阻力特性系数,Sl取值范围为10Pa/(kg/s)2≤Sl≤500Pa/(kg/s)2,ml,t为第l条管道在第t个调度时段的流量;
(2-3)电‐热耦合多能流系统中热网的循环泵水力特性方程如下:
其中,HP,t为循环泵在第t个调度时段的扬程,H0为循环泵静扬程,Sp为循环泵阻力系数,H0和Sp由循环泵的出厂说明书获取,mP,t为第t个调度时段流过循环泵的流量;
(2-4)电‐热耦合多能流系统中热网管道热量损失方程如下:
其中,Te,l,t为热网中第l条管道第t个调度时段的末端温度,Th,l,t为第l条管道第t个调度时段的首端温度,Ta,l,t为第l条管道在第t个调度时段的环境温度,ml,t为第l条管道在第t个调度时段的流量,Ll为第l条管道的长度,Cp为水的比热容,比热容的取值为4182焦耳/(千克·摄氏度),λ为管道单位长度的传热系数,λ从电-热耦合多能流系统的能量管理系统中获取,e为自然对数;
(2-5)电‐热耦合多能流系统的热网中多管道汇合点的温度方程:
其中,ml,t为第l条管道在第t个调度时段的流量,Te,l,t为热网中第l条管道第t个调度时段的末端温度,Th,l,t为第l条管道第t个调度时段的首端温度,为流出热网节点n的所有管道编号的集合,为流入热网节点n的所有管道编号的集合,QJ,n,t是第n个热网节点第t个调度时刻的热功率;
(2-6)通过电‐热联供机组耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,pb,t为第b台电-热联供机组在第t个调度时段的有功功率,q为第b台电-热联供机组在第t个调度时段的热功率,为第b台电-热联供机组运行可行域近似多边形的第k个顶点的横坐标,为第b台电-热联供机组运行可行域近似多边形的第k个顶点的纵坐标,为第b台电-热联供机组在第t个调度时段的第k个组合系数, NKb为第b台电-热联供机组的运行可行域近似多边形的顶点个数,电-热联供机组运行可行域近似多边形从电-热联供机组的出厂说明书中获取;
(2-7)通过循环泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,PP,t为循环泵第t个调度时段消耗的有功功率,g为重力加速度,ηP为循环泵效率,ηP的取值范围为0~1,mP,t为第t个调度时段的流过循环泵的流量,HP,t为循环泵第t个调度时段的的扬程;
(2-8)通过热泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
Php,t=ChpQhp,t,t=1,2,...,16
其中,Qhp,t为电‐热耦合多能流系统中第t个调度时段热泵发出的热功率,Php,t为第t个调度时段热泵消耗的电功率,Chp为热泵的产热效率,Chp从热泵的出厂说明书中获得;
(3)设定电-热耦合多能流系统中电网与热网稳态安全运行的不等式约束条件,包括:
(3-1)电-热耦合多能流系统的电网中第i个节点第t个调度时段的电压幅值Ui,t在设定的电网安全运行电压的上限值、下限值U i、之间运行,U i为第i个节点额定电压的0.95倍,为第i个节点额定电压的1.05倍:
(3-2)电-热耦合多能流系统的电网中第l条线路第t个调度时段的传输容量小于或等于设定的电网安全运行传输容量的最大值
(3-3)电-热耦合多能流系统的电网中电‐热联供机组有功功率的爬坡约束:
其中,和分别为第b台电-热联供机组有功功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,pb,t和pb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-4)电-热耦合多能流系统的电网中电‐热联供机组热功率的爬坡约束:
其中,和分别为第b台电-热联供机组热功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,qb,t和qb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-5)电-热耦合多能流系统的电网中火电机组有功功率的爬坡约束:
其中,和分别为第x台火电机组有功功率的向上和向下爬坡速率,和从火电机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,px,t和px,t-1分别为第x台火电机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-6)电-热耦合多能流系统的电网中第b台电-热联供机组第t个调度时段的有功功率pb,t在设定的电网安全运行第b台电-热联供机组有功功率的上、下限值 p b之间:
(3-7)电-热耦合多能流系统的电网中第b台电-热联供机组第t个调度时段的热功率qb,t在设定的电网安全运行第b台电-热联供机组热功率的上限值和下限值q b之间:
(3-8)电-热耦合多能流系统的电网中第x台火电机组第t个调度时段的有功功率px,t在设定的电网安全运行第x台火电机组有功功率的上限值p x和下限值之间:
(3-9)电-热耦合多能流系统的热网中第l条管道第t个调度时段的流量ml小于或等于热网安全运行流量的上限值
(3-10)电-热耦合多能流系统的热网中第i个换热站第t个调度时段回水温度Ti,t在设定的热网安全运行回水温度的上限值T和下限值之间:
(4)采用内点法,将上述步骤(1)中的方程作为目标函数,将上述步骤(2)和步骤(3)的所有方程作为约束条件,求解得到电-热耦合多能流系统中每台电‐热联供机组的有功功率和热功率,每台火电机组的有功功率,每台热泵的热功率,每台循环泵消耗的有功功率,作为电‐热耦合多能流系统的日内滚动调度方案,实现电‐热耦合多能流系统的日内滚动调度。
Claims (1)
1.一种电‐热耦合多能流系统的日内滚动调度方法,其特征在于该方法包括以下步骤:
(1)建立一个电-热耦合多能流系统优化调度的目标函数:
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>16</mn>
</munderover>
<mrow>
<mo>(</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>F</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>x</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>T</mi>
<mi>U</mi>
</mrow>
</msub>
</munderover>
<msub>
<mi>F</mi>
<mrow>
<mi>T</mi>
<mi>U</mi>
</mrow>
</msub>
<mo>(</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>,</mo>
</mrow>
其中,pb,t为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的有功功率,qb,t为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的热功率,N为电-热耦合多能流系统中电‐热联供机组的总台数,F(pb,t,qb,t)为电-热耦合多能流系统中第b台电‐热联供机组在第t个调度时段的运行成本,px,t为电-热耦合多能流系统中第x台火电机组在第t个调度时段的有功功率,NTU为电-热耦合多能流系统的电网中火电机组的总台数,FTU(px,t)为电-热耦合多能流系统中电网的第x台火电机组在第t个调度时段的运行成本,Δt是相邻两个调度时刻的时间间隔,Δt的取值为15分钟;
(2)设定电-热耦合多能流系统中电网与热网稳态安全运行的等式约束条件,包括:
(2-1)电-热耦合多能流系统中的电网潮流方程如下:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>16</mn>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>&Element;</mo>
<mi>i</mi>
</mrow>
</munder>
<msub>
<mi>U</mi>
<mi>j</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>16</mn>
</mrow>
其中,Pi,t为电-热耦合多能流系统的电网中节点i在第t个调度时段的注入有功功率,Qi,t为电网中节点i在第t个调度时段的注入无功功率,θi,t和θj,t分别为节点i和节点j在第t个调度时段的电压相角,Ui,t和Uj,t分别为节点i和节点j在第t个调度时段的电压幅值,Gij为电网节点导纳矩阵Y中第i行、第j列元素的实部,Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部,电网节点导纳矩阵Y从电-热耦合多能流系统的能量管理系统中获取,n为电网的节点总数;
(2-2)电‐热耦合多能流系统中热网的管道压力损失方程如下:
ΔHl,t=Slml,t|ml,t|,t=1,2,...,16,
其中,ΔHl,t为电‐热耦合多能流系统中热网的第l条管道在第t个调度时段的压力损失,Sl为第l条管道的阻力特性系数,Sl取值范围为10Pa/(kg/s)2≤Sl≤500Pa/(kg/s)2,ml,t为第l条管道在第t个调度时段的流量;
(2-3)电‐热耦合多能流系统中热网的循环泵水力特性方程如下:
<mrow>
<msub>
<mi>H</mi>
<mrow>
<mi>P</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>p</mi>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>P</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>,</mo>
</mrow>
其中,HP,t为循环泵在第t个调度时段的扬程,H0为循环泵静扬程,Sp为循环泵阻力系数,H0和Sp由循环泵的出厂说明书获取,mP,t为第t个调度时段流过循环泵的流量;
(2-4)电‐热耦合多能流系统中热网管道热量损失方程如下:
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>&lambda;L</mi>
<mi>l</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<msub>
<mi>m</mi>
<mi>l</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
</mrow>
其中,Te,l,t为热网中第l条管道第t个调度时段的末端温度,Th,l,t为第l条管道第t个调度时段的首端温度,Ta,l,t为第l条管道在第t个调度时段的环境温度,ml,t为第l条管道在第t个调度时段的流量,Ll为第l条管道的长度,Cp为水的比热容,比热容的取值为4182焦耳/(千克·摄氏度),λ为管道单位长度的传热系数,λ从电-热耦合多能流系统的能量管理系统中获取,e为自然对数;
(2-5)电‐热耦合多能流系统的热网中多管道汇合点的温度方程:
<mrow>
<mrow>
<mo>(</mo>
<mrow>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>&Element;</mo>
<msubsup>
<mi>S</mi>
<mi>n</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msubsup>
</mrow>
</munder>
<msub>
<mi>m</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<munder>
<mi>&Sigma;</mi>
<mrow>
<mi>l</mi>
<mo>&Element;</mo>
<msubsup>
<mi>S</mi>
<mi>n</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msubsup>
</mrow>
</munder>
<msub>
<mi>m</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mrow>
<mi>e</mi>
<mo>,</mo>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>J</mi>
<mo>,</mo>
<mi>n</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>,</mo>
</mrow>
其中,ml,t为第l条管道在第t个调度时段的流量,Te,l,t为热网中第l条管道第t个调度时段的末端温度,Th,l,t为第l条管道第t个调度时段的首端温度,为流出热网节点n的所有管道编号的集合,为流入热网节点n的所有管道编号的集合,QJ,n,t是第n个热网节点第t个调度时刻的热功率;
(2-6)通过电‐热联供机组耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>NK</mi>
<mi>b</mi>
</msub>
</mrow>
</msubsup>
<msubsup>
<mi>&alpha;</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>k</mi>
</msubsup>
<msubsup>
<mi>P</mi>
<mi>b</mi>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<msub>
<mi>NK</mi>
<mi>b</mi>
</msub>
</mrow>
</msubsup>
<msubsup>
<mi>&alpha;</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mi>k</mi>
</msubsup>
<msubsup>
<mi>Q</mi>
<mi>b</mi>
<mi>k</mi>
</msubsup>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>,</mo>
</mrow>
其中,pb,t为第b台电-热联供机组在第t个调度时段的有功功率,q为第b台电-热联供机组在第t个调度时段的热功率,为第b台电-热联供机组运行可行域近似多边形的第k个顶点的横坐标,为第b台电-热联供机组运行可行域近似多边形的第k个顶点的纵坐标,为第b台电-热联供机组在第t个调度时段的第k个组合系数, NKb为第b台电-热联供机组的运行可行域近似多边形的顶点个数,电-热联供机组运行可行域近似多边形从电-热联供机组的出厂说明书中获取;
(2-7)通过循环泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
<mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>P</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mrow>
<mi>P</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>gH</mi>
<mrow>
<mi>P</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msup>
<mn>10</mn>
<mn>6</mn>
</msup>
<msub>
<mi>&eta;</mi>
<mi>P</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>,</mo>
</mrow>
其中,PP,t为循环泵第t个调度时段消耗的有功功率,g为重力加速度,ηP为循环泵效率,ηP的取值范围为0~1,mP,t为第t个调度时段的流过循环泵的流量,HP,t为循环泵第t个调度时段的的扬程;
(2-8)通过热泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
Php,t=ChpQhp,t,t=1,2,...,16
其中,Qhp,t为电‐热耦合多能流系统中第t个调度时段热泵发出的热功率,Php,t为第t个调度时段热泵消耗的电功率,Chp为热泵的产热效率,Chp从热泵的出厂说明书中获得;
(3)设定电-热耦合多能流系统中电网与热网稳态安全运行的不等式约束条件,包括:
(3-1)电-热耦合多能流系统的电网中第i个节点第t个调度时段的电压幅值Ui,t在设定的电网安全运行电压的上限值、下限值U i、之间运行,U i为第i个节点额定电压的0.95倍,为第i个节点额定电压的1.05倍:
<mrow>
<msub>
<munder>
<mi>U</mi>
<mo>&OverBar;</mo>
</munder>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>U</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(3-2)电-热耦合多能流系统的电网中第l条线路第t个调度时段的传输容量小于或等于设定的电网安全运行传输容量的最大值
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>S</mi>
<mo>&OverBar;</mo>
</mover>
<mi>l</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(3-3)电-热耦合多能流系统的电网中电‐热联供机组有功功率的爬坡约束:
<mrow>
<mo>-</mo>
<msubsup>
<mi>&Delta;p</mi>
<mi>b</mi>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>&Delta;p</mi>
<mi>b</mi>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
其中,和分别为第b台电-热联供机组有功功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,pb,t和pb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-4)电-热耦合多能流系统的电网中电‐热联供机组热功率的爬坡约束:
<mrow>
<mo>-</mo>
<msubsup>
<mi>&Delta;q</mi>
<mi>b</mi>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&le;</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>&Delta;q</mi>
<mi>b</mi>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
其中,和分别为第b台电-热联供机组热功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,qb,t和qb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-5)电-热耦合多能流系统的电网中火电机组有功功率的爬坡约束:
<mrow>
<mo>-</mo>
<msubsup>
<mi>&Delta;p</mi>
<mi>x</mi>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>&le;</mo>
<msubsup>
<mi>&Delta;p</mi>
<mi>x</mi>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
其中,和分别为第x台火电机组有功功率的向上和向下爬坡速率,和从火电机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,px,t和px,t-1分别为第x台火电机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-6)电-热耦合多能流系统的电网中第b台电-热联供机组第t个调度时段的有功功率pb,t在设定的电网安全运行第b台电-热联供机组有功功率的上、下限值 p b之间:
<mrow>
<msub>
<munder>
<mi>p</mi>
<mo>&OverBar;</mo>
</munder>
<mi>b</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>b</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(3-7)电-热耦合多能流系统的电网中第b台电-热联供机组第t个调度时段的热功率qb,t在设定的电网安全运行第b台电-热联供机组热功率的上限值和下限值q b之间:
<mrow>
<msub>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>b</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<munder>
<mi>q</mi>
<mo>&OverBar;</mo>
</munder>
<mi>b</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(3-8)电-热耦合多能流系统的电网中第x台火电机组第t个调度时段的有功功率px,t在设定的电网安全运行第x台火电机组有功功率的上限值p x和下限值之间:
<mrow>
<msub>
<munder>
<mi>p</mi>
<mo>&OverBar;</mo>
</munder>
<mi>x</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>p</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>p</mi>
<mo>&OverBar;</mo>
</mover>
<mi>x</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(3-9)电-热耦合多能流系统的热网中第l条管道第t个调度时段的流量ml小于或等于热网安全运行流量的上限值
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>l</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mover>
<mi>m</mi>
<mo>&OverBar;</mo>
</mover>
<mi>l</mi>
</msub>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(3-10)电-热耦合多能流系统的热网中第i个换热站第t个调度时段回水温度Ti,t在设定的热网安全运行回水温度的上限值T和下限值之间:
<mrow>
<munder>
<mi>T</mi>
<mo>&OverBar;</mo>
</munder>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&le;</mo>
<mover>
<mi>T</mi>
<mo>&OverBar;</mo>
</mover>
<mo>,</mo>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mn>16</mn>
<mo>;</mo>
</mrow>
(4)采用内点法,将上述步骤(1)中的方程作为目标函数,将上述步骤(2)和步骤(3)的所有方程作为约束条件,求解得到电-热耦合多能流系统中每台电‐热联供机组的有功功率和热功率,每台火电机组的有功功率,每台热泵的热功率,每台循环泵消耗的有功功率,作为电‐热耦合多能流系统的日内滚动调度方案,实现电‐热耦合多能流系统的日内滚动调度。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710989015.9A CN107732983A (zh) | 2017-10-22 | 2017-10-22 | 一种电‑热耦合多能流系统的日内滚动调度方法 |
PCT/CN2017/114465 WO2019075872A1 (zh) | 2017-10-22 | 2017-12-04 | 一种电-热耦合多能流系统的日内滚动调度方法 |
US16/842,566 US11306923B2 (en) | 2017-10-22 | 2020-04-07 | Intra-day rolling scheduling method for integrated heat and electricity system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710989015.9A CN107732983A (zh) | 2017-10-22 | 2017-10-22 | 一种电‑热耦合多能流系统的日内滚动调度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107732983A true CN107732983A (zh) | 2018-02-23 |
Family
ID=61213219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710989015.9A Pending CN107732983A (zh) | 2017-10-22 | 2017-10-22 | 一种电‑热耦合多能流系统的日内滚动调度方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11306923B2 (zh) |
CN (1) | CN107732983A (zh) |
WO (1) | WO2019075872A1 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108985524A (zh) * | 2018-08-07 | 2018-12-11 | 南京南瑞继保电气有限公司 | 一种多能互补系统的协调控制方法 |
CN115903549A (zh) * | 2023-01-06 | 2023-04-04 | 国网浙江省电力有限公司金华供电公司 | 基于TwinCAT3的综合能源系统的调度策略筛选方法及装置 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112884191A (zh) * | 2019-11-30 | 2021-06-01 | 国网天津市电力公司电力科学研究院 | 一种基于网源协调的热电日前调度模型及计算方法 |
CN112539449B (zh) * | 2020-10-21 | 2022-02-01 | 西安交通大学 | 一种多能耦合的恒温供水系统及其优化方法 |
CN112906292B (zh) * | 2021-01-26 | 2024-02-23 | 西安热工研究院有限公司 | 热电联产机组厂级热电负荷在线优化分配的方法、系统、设备及存储介质 |
CN114517984B (zh) * | 2022-04-18 | 2023-07-28 | 国网天津市电力公司城南供电分公司 | 一种换热站蓄热式电锅炉设备控制装置及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014241689A (ja) * | 2013-06-12 | 2014-12-25 | 日本電気株式会社 | 負荷分散システム、負荷分散装置、負荷分散方法およびプログラム |
CN105046395A (zh) * | 2015-05-15 | 2015-11-11 | 华南理工大学 | 一种含多类型新能源的电力系统日内滚动计划编制方法 |
CN106056251A (zh) * | 2016-06-12 | 2016-10-26 | 清华大学 | 一种电‑热耦合多能流系统的优化调度方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6021402A (en) * | 1997-06-05 | 2000-02-01 | International Business Machines Corporaiton | Risk management system for electric utilities |
US8126685B2 (en) * | 2006-04-12 | 2012-02-28 | Edsa Micro Corporation | Automatic real-time optimization and intelligent control of electrical power distribution and transmission systems |
US8396572B2 (en) * | 2009-09-11 | 2013-03-12 | Siemens Corporation | System and method for energy plant optimization using mixed integer-linear programming |
US9335748B2 (en) * | 2010-07-09 | 2016-05-10 | Emerson Process Management Power & Water Solutions, Inc. | Energy management system |
CN106022624B (zh) * | 2016-05-27 | 2019-07-26 | 清华大学 | 一种电-热耦合多能流网络状态估计方法 |
CN105958531B (zh) * | 2016-05-27 | 2018-10-19 | 清华大学 | 一种电-气耦合多能流网络状态估计方法 |
CN106845671A (zh) * | 2016-12-12 | 2017-06-13 | 长沙理工大学 | 一种多能流系统多目标最优潮流模型及其求解方法 |
-
2017
- 2017-10-22 CN CN201710989015.9A patent/CN107732983A/zh active Pending
- 2017-12-04 WO PCT/CN2017/114465 patent/WO2019075872A1/zh active Application Filing
-
2020
- 2020-04-07 US US16/842,566 patent/US11306923B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014241689A (ja) * | 2013-06-12 | 2014-12-25 | 日本電気株式会社 | 負荷分散システム、負荷分散装置、負荷分散方法およびプログラム |
CN105046395A (zh) * | 2015-05-15 | 2015-11-11 | 华南理工大学 | 一种含多类型新能源的电力系统日内滚动计划编制方法 |
CN106056251A (zh) * | 2016-06-12 | 2016-10-26 | 清华大学 | 一种电‑热耦合多能流系统的优化调度方法 |
Non-Patent Citations (1)
Title |
---|
卢志刚,等: "考虑储热装置与碳捕集设备的风电消纳低碳经济调度", 《电工技术学报》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108985524A (zh) * | 2018-08-07 | 2018-12-11 | 南京南瑞继保电气有限公司 | 一种多能互补系统的协调控制方法 |
CN108985524B (zh) * | 2018-08-07 | 2020-11-17 | 南京南瑞继保电气有限公司 | 一种多能互补系统的协调控制方法 |
CN115903549A (zh) * | 2023-01-06 | 2023-04-04 | 国网浙江省电力有限公司金华供电公司 | 基于TwinCAT3的综合能源系统的调度策略筛选方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2019075872A1 (zh) | 2019-04-25 |
US11306923B2 (en) | 2022-04-19 |
US20200232654A1 (en) | 2020-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106056251B (zh) | 一种电-热耦合多能流系统的优化调度方法 | |
CN107732983A (zh) | 一种电‑热耦合多能流系统的日内滚动调度方法 | |
CN109978625B (zh) | 一种计及电热气网络的综合能源系统多目标运行优化方法 | |
CN107800158B (zh) | 一种兼顾经济与能效的电-热耦合多能流系统优化调度方法 | |
CN106056478B (zh) | 一种电-热耦合系统中热网的区间潮流计算方法 | |
Geidl et al. | Operational and structural optimization of multi‐carrier energy systems | |
Beigvand et al. | Optimal operation of multicarrier energy systems using time varying acceleration coefficient gravitational search algorithm | |
CN103246263B (zh) | 冷热电联供微网系统的通用优化调度方法 | |
CN107609680B (zh) | 一种多热源环状集中供热管网水力工况优化调度方法 | |
CN110110913A (zh) | 大型园区综合能源系统能源站优化配置方法 | |
CN106067677A (zh) | 一种电‑热耦合多能流系统静态安全分析方法 | |
CN107910871A (zh) | 一种能源互联网多能互补控制方法 | |
CN101882785B (zh) | 一种含分布式电源的智能配电网络潮流分析装置及方法 | |
CN103728881A (zh) | 一种多楼宇冷热电联供系统的优化运行方法 | |
CN109345012A (zh) | 基于综合评价指标的园区能源互联网运行优化方法 | |
CN106651628A (zh) | 基于图论的区域冷热电综合能源优化配置方法及装置 | |
CN106169750A (zh) | 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 | |
CN106056246A (zh) | 一种冷热电多能流微电网考虑运行的储能容量优化方法 | |
CN109255194A (zh) | 一种含p2h的电-热互联综合能源系统光伏消纳分析方法 | |
CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
Tian et al. | Energy network flow model and optimization based on energy hub for big harbor industrial park | |
CN105406517B (zh) | 基于有限时间平均一致性算法的电力系统经济调度方法 | |
CN111340271A (zh) | 一种基于供热相量模型的电-热多能流系统优化调度方法 | |
CN111144642B (zh) | 一种综合能源系统优化方法及系统 | |
CN109446545A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180223 |
|
WD01 | Invention patent application deemed withdrawn after publication |