CN106056251A - 一种电‑热耦合多能流系统的优化调度方法 - Google Patents
一种电‑热耦合多能流系统的优化调度方法 Download PDFInfo
- Publication number
- CN106056251A CN106056251A CN201610408959.8A CN201610408959A CN106056251A CN 106056251 A CN106056251 A CN 106056251A CN 201610408959 A CN201610408959 A CN 201610408959A CN 106056251 A CN106056251 A CN 106056251A
- Authority
- CN
- China
- Prior art keywords
- electric
- streaming system
- thermal
- thermal coupling
- multipotency streaming
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000005457 optimization Methods 0.000 title claims abstract description 13
- 230000008878 coupling Effects 0.000 claims description 115
- 238000010168 coupling process Methods 0.000 claims description 115
- 238000005859 coupling reaction Methods 0.000 claims description 115
- 238000005485 electric heating Methods 0.000 claims description 21
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000010977 unit operation Methods 0.000 claims description 8
- 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
- 230000009194 climbing 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
- 238000004458 analytical method Methods 0.000 abstract description 4
- 238000010438 heat treatment Methods 0.000 abstract 1
- 230000005611 electricity Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 239000002699 waste material Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 230000004899 motility Effects 0.000 description 1
- 230000001737 promoting effect 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Remote Monitoring And Control Of Power-Distribution Networks (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种电‐热耦合多能流系统的优化调度方法,属于含多种能源形式的电网运行和控制技术领域。本方法考虑电‑热系统的相互影响,实现了电‑热耦合多能流系统的优化调度。相比独立地对供电、供热系统进行优化调度分析,不仅能得到更优的调度方案(总运行成本或网络损耗更小等),还提高了调度的灵活性。该方法可以应用于电‑热耦合多能流系统的调度计划制定,有利于提高电‑热耦合多能流系统的用能效率,减少运行成本。
Description
技术领域
本发明涉及一种电‐热耦合多能流系统的优化调度方法,属于含多种能源形式的电网运行和控制技术领域。
背景技术
能源综合利用是提高综合能源利用效率、促进可再生能源消纳的重要途径,通过打破原来电、热、冷、气、交通等能流子系统相对割裂的状态,实现多类型能源开放互联,构建多能流系统。多能流是指多种类型的能量流,表示电、热、冷、气、交通等能量流的相互耦合、转换和传输。多能流系统相比传统相互割裂的能源系统,其带来的效益包括:1)通过多类型能源的梯级开发利用和智能管理,可以降低能源消耗和浪费,提高综合能源利用效率,并有助于减少总的用能成本;2)利用不同能源的特性差异和互补、转换,有助于提高消纳间歇式可再生能源的能力;3)通过多能源的转供、互补和协调控制,有助于提高供能的可靠性,并为电网的运行提供更多可调控资源;4)通过多能流系统的协同规划和建设,可以减少基础设施的重复建设和浪费,提高资产利用率。
多能流系统一方面具有可观的效益,另一方面也使原本复杂的能源系统更加复杂。多能流系统由多个能流子系统组成,这些能流子系统之间相互作用和影响,使得多能流系统复杂度显著增加,体现出许多新的特性,传统各个能流单独分析的方法已经难以适应新的要求,亟需发展出新的多能流分析方法。在我国,越来越多的热电联产机组、热泵、电锅炉等耦合元件客观上增强了电-热之间的互联,促进了电-热耦合多能流系统的发展,也对电-热耦合多能流系统的运行和控制技术提出了新的要求。
多能源系统优化调度是指当系统的结构参数和负荷情况都已给定时,调节可利用的控制变量(如电网中发电机的输出功率、热网中泵的扬程等)来找到能满足所有运行约束条件的,并使系统的某一性能指标(如总运行成本或网络损耗)达到最优值下的潮流分布。目前这方面的研究主要集中在单个独立的系统,为了使得电-热耦合多能流系统的运行成本最低,需要研究电-热耦合多能流系统优化调度方法。
发明内容
本发明的目的是提出一种电‐热耦合多能流系统的优化调度方法,以弥补现有领域研 究的空白,建立电-热耦合多能流系统优化调度模型,实现电-热耦合多能流系统的优化调度。
本发明提出的电‐热耦合多能流系统的优化调度方法,包括以下步骤:
(1)建立一个电-热耦合多能流系统优化调度的目标函数:
其中,pb为电-热耦合多能流系统中第b台电‐热联供机组的有功功率,qb为电-热耦合多能流系统中第b台电‐热联供机组的热功率,N为电-热耦合多能流系统中电‐热联供机组的总台数,F(pb,qb)为电-热耦合多能流系统中第b台电‐热联供机组的运行成本,px为电-热耦合多能流系统中第x台火电机组的有功功率,NTU为电-热耦合多能流系统中火电机组的总台数,FTU(px)为电-热耦合多能流系统中第x台火电机组的运行成本;
(2)设定电-热耦合多能流系统中电网与热网稳态安全运行的等式约束条件,包括:
(2-1)电-热耦合多能流系统中的电网潮流方程如下:
其中,Pi为电网中节点i的注入有功功率,Qi为电网中节点i的注入无功功率,θi、θj分别为节点i、节点j的电压相角,Ui和Uj分别为节点i和节点j的电压幅值,Gij为电网节点导纳矩阵Y第i行、第j列元素的实部,Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部,电网节点导纳矩阵Y从电-热耦合多能流系统的能量管理系统中获取;
(2-2)电‐热耦合多能流系统中热网的管道压力损失方程如下:
ΔHl=Slml|ml|,
其中,ΔHl为热网中第l条管道的压力损失,Sl为第l条管道的阻力特性系数,Sl取值范围为10Pa/(kg/s)2≤Sl≤500Pa/(kg/s)2,ml为第l条管道的流量;
(2-3)电‐热耦合多能流系统中热网的循环泵水力特性方程如下:
HP=H0-Spm2,
其中,HP为循环泵扬程,H0为循环泵静扬程,Sp为循环泵阻力系数,H0和Sp由循 环泵的出厂说明书获取,m为流过循环泵的流量;
(2-4)电‐热耦合多能流系统中热网管道热量损失方程如下:
其中,Te,l为热网中第l条管道的末端温度,Th,l为第l条管道的首端温度,Ta,l为第l条管道所在的环境温度,ml为第l条管道的流量,Ll为第l条管道的长度,Cp为水的比热容,比热容的取值为4182焦耳/(千克·摄氏度),λ为管道单位长度的传热系数,λ从电-热耦合多能流系统的能量管理系统中获取;
(2-5)电‐热耦合多能流系统的热网中多管道汇合点的温度方程:
其中,为流出多管道汇合点的流量,为流入多管道汇合点的流量,Tout为流出多管道汇合点的水的温度,Tin为流入多管道汇合点的水的温度,QJ是多管道汇合点的热功率;
(2-6)通过电‐热联供机组耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,p为电-热联供机组的有功功率,q为电-热联供机组的热功率,Pk为电-热联供机组运行可行域近似多边形的第k个顶点的横坐标,Qk为电-热联供机组运行可行域近似多边形的第k个顶点的纵坐标,αk为组合系数,0≤αk≤1,NK为电-热联供机组的运行可行域近似多边形的顶点个数,电-热联供机组运行可行域近似多边形从电-热联供机组的出厂说明书中获取;
(2-7)通过循环泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,PP为循环泵消耗的有功功率,g为重力加速度,ηP为循环泵效率,ηP的取值范围为0~1,mP为流过循环泵的流量,HP为循环泵的扬程;
(2-8)通过热泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
Php=ChpQhp
其中,Qhp为电‐热耦合多能流系统中热泵发出的热功率,Php为热泵消耗的电功率,Chp为热泵的产热效率,Chp从热泵的出厂说明书中获得;
(3)设定电-热耦合多能流系统中电网与热网稳态安全运行的不等式约束条件,包括:
(3-1)电-热耦合多能流系统的电网中第i个节点的电压幅值Ui在设定的电网安全运行电压的上、下限值U i、之间运行,U i为第i个节点额定电压的0.95倍,为第i个节点额定电压的1.05倍:
(3-2)电-热耦合多能流系统的电网中第l条线路的传输容量小于或等于设定的电网安全运行传输容量的最大值
(3-3)电-热耦合多能流系统的电网中电‐热联供机组或有功功率的爬坡约束:
其中,和分别为第b台电-热联供机组有功功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,pb,t和pb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-4)电-热耦合多能流系统的电网中非燃气机组有功功率的爬坡约束:
其中,和分别为第x台火电机组有功功率的向上和向下爬坡速率,和从火电机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,px,t和px,t-1分别为第x台火电机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-5)电-热耦合多能流系统的电网中第b台电-热联供机组的有功功率pb在设定的电网安全运行第b台电-热联供机组有功功率的上、下限值 p b之间:
(3-6)电-热耦合多能流系统的电网中第x台火电机组的有功功率px在设定的电网安全运行第x台火电机组有功功率的上、下限值 p x之间:
(3-7)电-热耦合多能流系统的热网中第l条管道的流量ml小于或等于热网安全运行流量的上限值
(3-8)电-热耦合多能流系统的热网中换热站回水温度T在设定的热网安全运行回水温度的上、下限值 T之间:
(4)采用内点法,将步骤(1)中的方程作为目标函数,将上述步骤(2)和步骤(3)的所有方程作为约束条件,求解得到电-热耦合多能流系统中每台电‐热联供机组的有功功率和热功率,作为电‐热耦合多能流系统的优化调度方案。
本发明提出的电-热耦合多能流系统优化调度方法,其特点和效果是:本方法考虑电-热系统的相互影响,实现了电-热耦合多能流系统的优化调度。相比独立地对供电、供热系统进行优化调度分析,不仅能得到更优的调度方案(总运行成本更低),还提高了调度的灵活性。该方法可以应用于电-热耦合多能流系统的调度计划制定,有利于提高电-热耦合多能流系统的用能效率,减少运行成本。
具体实施方式
本发明提出的电‐热耦合多能流系统的优化调度方法,包括以下步骤:
(1)建立一个电-热耦合多能流系统优化调度的目标函数:
其中,pb为电-热耦合多能流系统中第b台电‐热联供机组的有功功率,qb为电-热耦合多能流系统中第b台电‐热联供机组的热功率,N为电-热耦合多能流系统中电‐热联供机组的总台数,F(pb,qb)为电-热耦合多能流系统中第b台电‐热联供机组的运行成本,px为电-热耦合多能流系统中第x台火电机组的有功功率,NTU为电-热耦合多能流系统中火电机组的总台数,FTU(px)为电-热耦合多能流系统中第x台火电机组的运行成本;
(2)设定电-热耦合多能流系统中电网与热网稳态安全运行的等式约束条件,包括:
(2-1)电-热耦合多能流系统中的电网潮流方程如下:
其中,Pi为电网中节点i的注入有功功率,Qi为电网中节点i的注入无功功率,θi、θj分别为节点i、节点j的电压相角,Ui和Uj分别为节点i和节点j的电压幅值,Gij为电网节点导纳矩阵Y第i行、第j列元素的实部,Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部,电网节点导纳矩阵Y从电-热耦合多能流系统的能量管理系统中获取;
(2-2)电‐热耦合多能流系统中热网的管道压力损失方程如下:
ΔHl=Slml|ml|,
其中,ΔHl为热网中第l条管道的压力损失,Sl为第l条管道的阻力特性系数,Sl取值范围为10Pa/(kg/s)2≤Sl≤500Pa/(kg/s)2,ml为第l条管道的流量;
(2-3)电‐热耦合多能流系统中热网的循环泵水力特性方程如下:
HP=H0-Spm2,
其中,HP为循环泵扬程,H0为循环泵静扬程,Sp为循环泵阻力系数,H0和Sp由循环泵的出厂说明书获取,m为流过循环泵的流量;
(2-4)电‐热耦合多能流系统中热网管道热量损失方程如下:
其中,Te,l为热网中第l条管道的末端温度,Th,l为第l条管道的首端温度,Ta,l为第l条管道所在的环境温度,ml为第l条管道的流量,Ll为第l条管道的长度,Cp为水的比热容,比热容的取值为4182焦耳/(千克·摄氏度),λ为管道单位长度的传热系数,λ从电-热耦合多能流系统的能量管理系统中获取;
(2-5)电‐热耦合多能流系统的热网中多管道汇合点的温度方程:
其中,为流出多管道汇合点的流量,为流入多管道汇合点的流量,Tout为流出多管道汇合点的水的温度,Tin为流入多管道汇合点的水的温度,QJ是多管道汇合点的热 功率;
(2-6)通过电‐热联供机组耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,p为电-热联供机组的有功功率,q为电-热联供机组的热功率,Pk为电-热联供机组运行可行域近似多边形的第k个顶点的横坐标,Qk为电-热联供机组运行可行域近似多边形的第k个顶点的纵坐标,αk为组合系数,0≤αk≤1,NK为电-热联供机组的运行可行域近似多边形的顶点个数,电-热联供机组运行可行域近似多边形从电-热联供机组的出厂说明书中获取;
(2-7)通过循环泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,PP为循环泵消耗的有功功率,g为重力加速度,ηP为循环泵效率,ηP的取值范围为0~1,mP为流过循环泵的流量,HP为循环泵的扬程;
(2-8)通过热泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
Php=ChpQhp
其中,Qhp为电‐热耦合多能流系统中热泵发出的热功率,Php为热泵消耗的电功率,Chp为热泵的产热效率,Chp从热泵的出厂说明书中获得;
(3)设定电-热耦合多能流系统中电网与热网稳态安全运行的不等式约束条件,包括:
(3-1)电-热耦合多能流系统的电网中第i个节点的电压幅值Ui在设定的电网安全运行电压的上、下限值U i、之间运行,U i为第i个节点额定电压的0.95倍,为第i个节点额定电压的1.05倍:
(3-2)电-热耦合多能流系统的电网中第l条线路的传输容量小于或等于设定的电网安全运行传输容量的最大值
(3-3)电-热耦合多能流系统的电网中电‐热联供机组或有功功率的爬坡约束:
其中,和分别为第b台电-热联供机组有功功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,pb,t和pb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-4)电-热耦合多能流系统的电网中非燃气机组有功功率的爬坡约束:
其中,和分别为第x台火电机组有功功率的向上和向下爬坡速率,和从火电机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,px,t和px,t-1分别为第x台火电机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-5)电-热耦合多能流系统的电网中第b台电-热联供机组的有功功率pb在设定的电网安全运行第b台电-热联供机组有功功率的上、下限值 p b之间:
(3-6)电-热耦合多能流系统的电网中第x台火电机组的有功功率px在设定的电网安全运行第x台火电机组有功功率的上、下限值 p x之间:
(3-7)电-热耦合多能流系统的热网中第l条管道的流量ml小于或等于热网安全运行流量的上限值
(3-8)电-热耦合多能流系统的热网中换热站回水温度T在设定的热网安全运行回水温度的上、下限值 T之间:
(4)采用内点法,将步骤(1)中的方程作为目标函数,将上述步骤(2)和步骤(3)的所有方程作为约束条件,求解得到电-热耦合多能流系统中每台电‐热联供机组的有功功率和热功率,作为电‐热耦合多能流系统的优化调度方案。
本发明方法中使用的内点法(Interior Point Method)求解方程是一种求解线性规划或非线性凸优化问题的算法,是一种公知技术。
Claims (1)
1.一种电‐热耦合多能流系统的优化调度方法,其特征在于,该方法包括以下步骤:
(1)建立一个电-热耦合多能流系统优化调度的目标函数:
其中,pb为电-热耦合多能流系统中第b台电‐热联供机组的有功功率,qb为电-热耦合多能流系统中第b台电‐热联供机组的热功率,N为电-热耦合多能流系统中电‐热联供机组的总台数,F(pb,qb)为电-热耦合多能流系统中第b台电‐热联供机组的运行成本,px为电-热耦合多能流系统中第x台火电机组的有功功率,NTU为电-热耦合多能流系统中火电机组的总台数,FTU(px)为电-热耦合多能流系统中第x台火电机组的运行成本;
(2)设定电-热耦合多能流系统中电网与热网稳态安全运行的等式约束条件,包括:
(2-1)电-热耦合多能流系统中的电网潮流方程如下:
其中,Pi为电网中节点i的注入有功功率,Qi为电网中节点i的注入无功功率,θi、θj分别为节点i、节点j的电压相角,Ui和Uj分别为节点i和节点j的电压幅值,Gij为电网节点导纳矩阵Y第i行、第j列元素的实部,Bij为电网节点导纳矩阵Y第i行、第j列元素的虚部,电网节点导纳矩阵Y从电-热耦合多能流系统的能量管理系统中获取;
(2-2)电‐热耦合多能流系统中热网的管道压力损失方程如下:
ΔHl=Slml|ml|,
其中,ΔHl为热网中第l条管道的压力损失,Sl为第l条管道的阻力特性系数,Sl取值范围为10Pa/(kg/s)2≤Sl≤500Pa/(kg/s)2,ml为第l条管道的流量;
(2-3)电‐热耦合多能流系统中热网的循环泵水力特性方程如下:
HP=H0-Spm2,
其中,HP为循环泵扬程,H0为循环泵静扬程,Sp为循环泵阻力系数,H0和Sp由循环泵的出厂说明书获取,m为流过循环泵的流量;
(2-4)电‐热耦合多能流系统中热网管道热量损失方程如下:
其中,Te,l为热网中第l条管道的末端温度,Th,l为第l条管道的首端温度,Ta,l为第l条管道所在的环境温度,ml为第l条管道的流量,Ll为第l条管道的长度,Cp为水的比热容,比热容的取值为4182焦耳/(千克·摄氏度),λ为管道单位长度的传热系数,λ从电-热耦合多能流系统的能量管理系统中获取;
(2-5)电‐热耦合多能流系统的热网中多管道汇合点的温度方程:
其中,为流出多管道汇合点的流量,为流入多管道汇合点的流量,Tout为流出多管道汇合点的水的温度,Tin为流入多管道汇合点的水的温度,QJ是多管道汇合点的热功率;
(2-6)通过电‐热联供机组耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,p为电-热联供机组的有功功率,q为电-热联供机组的热功率,Pk为电-热联供机组运行可行域近似多边形的第k个顶点的横坐标,Qk为电-热联供机组运行可行域近似多边形的第k个顶点的纵坐标,αk为组合系数,0≤αk≤1,NK为电-热联供机组的运行可行域近似多边形的顶点个数,电-热联供机组运行可行域近似多边形从电-热联供机组的出厂说明书中获取;
(2-7)通过循环泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
其中,PP为循环泵消耗的有功功率,g为重力加速度,ηP为循环泵效率,ηP的取值范围为0~1,mP为流过循环泵的流量,HP为循环泵的扬程;
(2-8)通过热泵耦合的电‐热耦合多能流系统中电网与热网之间的耦合方程:
Php=ChpQhp
其中,Qhp为电‐热耦合多能流系统中热泵发出的热功率,Php为热泵消耗的电功率,Chp为热泵的产热效率,Chp从热泵的出厂说明书中获得;
(3)设定电-热耦合多能流系统中电网与热网稳态安全运行的不等式约束条件,包括:
(3-1)电-热耦合多能流系统的电网中第i个节点的电压幅值Ui在设定的电网安全运行电压的上、下限值U i、之间运行,U i为第i个节点额定电压的0.95倍,为第i个节点额定电压的1.05倍:
(3-2)电-热耦合多能流系统的电网中第l条线路的传输容量小于或等于设定的电网安全运行传输容量的最大值
(3-3)电-热耦合多能流系统的电网中电‐热联供机组或有功功率的爬坡约束:
其中,和分别为第b台电-热联供机组有功功率的向上和向下爬坡速率,和从电-热联供机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,pb,t和pb,t-1分别为第b台电-热联供机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-4)电-热耦合多能流系统的电网中非燃气机组有功功率的爬坡约束:
其中,和分别为第x台火电机组有功功率的向上和向下爬坡速率,和从火电机组的出厂说明书中获得,Δt为相邻两个调度时段的时间间隔,px,t和px,t-1分别为第x台火电机组在第t个调度时段和第t-1个调度时段的有功功率;
(3-5)电-热耦合多能流系统的电网中第b台电-热联供机组的有功功率pb在设定的电网安全运行第b台电-热联供机组有功功率的上、下限值 p b之间:
(3-6)电-热耦合多能流系统的电网中第x台火电机组的有功功率px在设定的电网安全运行第x台火电机组有功功率的上、下限值 p x之间:
(3-7)电-热耦合多能流系统的热网中第l条管道的流量ml小于或等于热网安全运行流量的上限值
(3-8)电-热耦合多能流系统的热网中换热站回水温度T在设定的热网安全运行回水温度的上、下限值 T之间:
(4)采用内点法,将步骤(1)中的方程作为目标函数,将上述步骤(2)和步骤(3)的所有方程作为约束条件,求解得到电-热耦合多能流系统中每台电‐热联供机组的有功功率和热功率,作为电‐热耦合多能流系统的优化调度方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610408959.8A CN106056251B (zh) | 2016-06-12 | 2016-06-12 | 一种电-热耦合多能流系统的优化调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610408959.8A CN106056251B (zh) | 2016-06-12 | 2016-06-12 | 一种电-热耦合多能流系统的优化调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106056251A true CN106056251A (zh) | 2016-10-26 |
CN106056251B CN106056251B (zh) | 2019-06-18 |
Family
ID=57170865
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610408959.8A Active CN106056251B (zh) | 2016-06-12 | 2016-06-12 | 一种电-热耦合多能流系统的优化调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106056251B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563674A (zh) * | 2017-10-09 | 2018-01-09 | 清华大学 | 一种考虑管道动态特性的电‑热耦合系统状态估计方法 |
CN107732983A (zh) * | 2017-10-22 | 2018-02-23 | 清华大学 | 一种电‑热耦合多能流系统的日内滚动调度方法 |
CN107800158A (zh) * | 2017-10-17 | 2018-03-13 | 清华大学 | 一种兼顾经济与能效的电‑热耦合多能流系统优化调度方法 |
CN108491973A (zh) * | 2018-03-21 | 2018-09-04 | 浙江大学 | 一种考虑需求侧响应的电-热联合系统优化调度方法 |
CN109472413A (zh) * | 2018-11-14 | 2019-03-15 | 南方电网科学研究院有限责任公司 | 考虑热管网传输特性的园区综合能源系统优化调度方法 |
WO2019075871A1 (zh) * | 2017-10-16 | 2019-04-25 | 清华大学 | 一种基于双侧等效模型的热网稳态运行状态估计方法 |
CN109711601A (zh) * | 2018-11-28 | 2019-05-03 | 国网浙江省电力有限公司电力科学研究院 | 电-气-热综合能源系统分布式优化调度方法和装置 |
CN110991845A (zh) * | 2019-11-25 | 2020-04-10 | 国网节能服务有限公司 | 一种电-热耦合系统分布式协同调度方法 |
US20210141405A1 (en) * | 2018-07-20 | 2021-05-13 | Tsinghua University | Method, apparatus, and storage medium for controlling heating system |
CN113048547A (zh) * | 2019-12-27 | 2021-06-29 | 中国电力科学研究院有限公司 | 一种综合能源供热系统的功率分配方法及装置 |
WO2021159893A1 (zh) * | 2020-02-13 | 2021-08-19 | 清华大学 | 一种基于供热相量模型的电-热多能流系统优化调度方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239981A (zh) * | 2014-10-11 | 2014-12-24 | 中冶南方工程技术有限公司 | 一种基于多目标优化算法的能源动态平衡与优化调度方法 |
CN105046369A (zh) * | 2015-08-13 | 2015-11-11 | 河海大学 | 一种基于能源中心的电气混联系统建模和优化调度方法 |
CN105576710A (zh) * | 2016-02-18 | 2016-05-11 | 东南大学 | 一种综合能源系统中的分布式电源的配置方法 |
-
2016
- 2016-06-12 CN CN201610408959.8A patent/CN106056251B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104239981A (zh) * | 2014-10-11 | 2014-12-24 | 中冶南方工程技术有限公司 | 一种基于多目标优化算法的能源动态平衡与优化调度方法 |
CN105046369A (zh) * | 2015-08-13 | 2015-11-11 | 河海大学 | 一种基于能源中心的电气混联系统建模和优化调度方法 |
CN105576710A (zh) * | 2016-02-18 | 2016-05-11 | 东南大学 | 一种综合能源系统中的分布式电源的配置方法 |
Non-Patent Citations (2)
Title |
---|
孙秋野 等: "能源互联网动态协调优化控制体系构建", 《中国电机工程学报》 * |
顾泽鹏 等: "考虑热网约束的电热能源集成系统运行优化及其风电消纳效益分析", 《中国电机工程学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563674A (zh) * | 2017-10-09 | 2018-01-09 | 清华大学 | 一种考虑管道动态特性的电‑热耦合系统状态估计方法 |
CN107563674B (zh) * | 2017-10-09 | 2020-07-10 | 清华大学 | 一种考虑管道动态特性的电-热耦合系统状态估计方法 |
WO2019075871A1 (zh) * | 2017-10-16 | 2019-04-25 | 清华大学 | 一种基于双侧等效模型的热网稳态运行状态估计方法 |
CN107800158B (zh) * | 2017-10-17 | 2019-09-27 | 清华大学 | 一种兼顾经济与能效的电-热耦合多能流系统优化调度方法 |
CN107800158A (zh) * | 2017-10-17 | 2018-03-13 | 清华大学 | 一种兼顾经济与能效的电‑热耦合多能流系统优化调度方法 |
CN107732983A (zh) * | 2017-10-22 | 2018-02-23 | 清华大学 | 一种电‑热耦合多能流系统的日内滚动调度方法 |
US11306923B2 (en) | 2017-10-22 | 2022-04-19 | Tsinghua University | Intra-day rolling scheduling method for integrated heat and electricity system |
CN108491973A (zh) * | 2018-03-21 | 2018-09-04 | 浙江大学 | 一种考虑需求侧响应的电-热联合系统优化调度方法 |
US20210141405A1 (en) * | 2018-07-20 | 2021-05-13 | Tsinghua University | Method, apparatus, and storage medium for controlling heating system |
CN109472413A (zh) * | 2018-11-14 | 2019-03-15 | 南方电网科学研究院有限责任公司 | 考虑热管网传输特性的园区综合能源系统优化调度方法 |
CN109711601A (zh) * | 2018-11-28 | 2019-05-03 | 国网浙江省电力有限公司电力科学研究院 | 电-气-热综合能源系统分布式优化调度方法和装置 |
CN110991845A (zh) * | 2019-11-25 | 2020-04-10 | 国网节能服务有限公司 | 一种电-热耦合系统分布式协同调度方法 |
CN110991845B (zh) * | 2019-11-25 | 2023-06-23 | 国网节能服务有限公司 | 一种电-热耦合系统分布式协同调度方法 |
CN113048547A (zh) * | 2019-12-27 | 2021-06-29 | 中国电力科学研究院有限公司 | 一种综合能源供热系统的功率分配方法及装置 |
WO2021159893A1 (zh) * | 2020-02-13 | 2021-08-19 | 清华大学 | 一种基于供热相量模型的电-热多能流系统优化调度方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106056251B (zh) | 2019-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106056251A (zh) | 一种电‑热耦合多能流系统的优化调度方法 | |
CN107800158B (zh) | 一种兼顾经济与能效的电-热耦合多能流系统优化调度方法 | |
CN107732983A (zh) | 一种电‑热耦合多能流系统的日内滚动调度方法 | |
CN109978625B (zh) | 一种计及电热气网络的综合能源系统多目标运行优化方法 | |
CN106096777B (zh) | 一种电-气耦合多能流系统的优化调度方法 | |
CN106339772B (zh) | 基于供热管网储热效益的热-电联合优化调度方法 | |
CN106056478B (zh) | 一种电-热耦合系统中热网的区间潮流计算方法 | |
CN107609680B (zh) | 一种多热源环状集中供热管网水力工况优化调度方法 | |
CN106067677B (zh) | 一种电-热耦合多能流系统静态安全分析方法 | |
CN110110913A (zh) | 大型园区综合能源系统能源站优化配置方法 | |
CN111324849B (zh) | 一种考虑热网特性的电热联合系统优化调度方法 | |
CN105490312B (zh) | 一种电力系统多源无功优化控制方法 | |
CN104779611A (zh) | 基于集中式和分布式双层优化策略的微电网经济调度方法 | |
CN103728881A (zh) | 一种多楼宇冷热电联供系统的优化运行方法 | |
CN103514374B (zh) | 电力系统在线滚动调度中不可行传输断面约束的辨识方法 | |
CN105427011A (zh) | 一种温控负荷聚合系统稳定性分析方法 | |
CN107666155A (zh) | 基于Markov模型的多能互补系统随机稳定性分析方法 | |
CN109255489A (zh) | 一种基于半不变量法的电-热互联综合能源系统概率能量流计算方法 | |
CN106169750A (zh) | 一种基于二阶锥松弛的主动配电网最大供电能力计算方法 | |
CN109255194A (zh) | 一种含p2h的电-热互联综合能源系统光伏消纳分析方法 | |
CN104242355A (zh) | 考虑最小弃风的风电场接入电网的位置和容量的控制方法 | |
CN110365062A (zh) | 一种基于Markov模型的多能系统协调控制方法 | |
CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
CN109256766A (zh) | 一种能源互联网中电热耦合系统 | |
CN107332290A (zh) | 一种基于直流线路的区域负荷转供方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20210127 Address after: 102208 817, unit 1, 8th floor, building 1, yard 30, Yuzhi East Road, Changping District, Beijing Patentee after: Beijing Qingyuan Intelligent Technology Co.,Ltd. Address before: 100084 No. 1 Tsinghua Yuan, Beijing, Haidian District Patentee before: TSINGHUA University |