CN114021997A - 计及网损的电热综合能源系统分布式经济调度方法及系统 - Google Patents
计及网损的电热综合能源系统分布式经济调度方法及系统 Download PDFInfo
- Publication number
- CN114021997A CN114021997A CN202111320160.0A CN202111320160A CN114021997A CN 114021997 A CN114021997 A CN 114021997A CN 202111320160 A CN202111320160 A CN 202111320160A CN 114021997 A CN114021997 A CN 114021997A
- Authority
- CN
- China
- Prior art keywords
- electric heating
- energy system
- total
- minimum
- heat
- 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
- 238000005485 electric heating Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000005457 optimization Methods 0.000 claims abstract description 36
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 230000005540 biological transmission Effects 0.000 claims description 56
- 238000010977 unit operation Methods 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 abstract description 8
- 238000005859 coupling reaction Methods 0.000 abstract description 8
- 238000004364 calculation method Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 12
- 238000004891 communication Methods 0.000 description 9
- 239000000126 substance Substances 0.000 description 8
- 238000010438 heat treatment Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 230000008878 coupling Effects 0.000 description 5
- 238000010248 power generation Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000008400 supply water Substances 0.000 description 2
- 241000764238 Isis Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013468 resource allocation Methods 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/06311—Scheduling, planning or task assignment for a person or group
-
- 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
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)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请涉及综合能源系统经济调度技术领域，提供了一种计及网损的电热综合能源系统分布式经济调度方法包括：根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；根据最小值优化模型，利用双层一致性算法求解得到系统运行总成本最小值。在实际应用过程中，系统优先调度增量成本小的机组出力以使系统运行总成本减小；本方法不仅能够很好地求解电热综合能源系统多约束、强耦合优化问题，而且迭代计算均匀分配给各个参与机组，最终具有较快的收敛速度和满意的收敛结果。
Description
技术领域
本申请涉及综合能源系统经济调度技术领域，具体涉及一种计及网损的电热综合能源系统分布式经济调度方法及系统。
背景技术
经济调度作为电力系统运行中技术和经济优化的重要部分，旨在满足发电单元运行约束的前提下，通过优化分配负荷需求并合理安排发电计划使得系统总运行成本最小。电力系统经济调度本质上是资源配置问题，在满足用户侧负荷需求和供电质量的前提下，引导供给侧制定合理的产能方案，减少企业运行成本的同时保障系统安全稳定运行。
经济调度解决方法总体上分为集中式和分布式，集中式算法要求系统控制中心与每一个发电单元进行信息交互，并收集所有需要的信息计算经济调度最优方案，最后通过发布调度指令安排所有发电单元安排出力计划。然而，集中式算法存在以下关键问题：首先系统控制中心需要较高的通信建设成本；其次巨大的计算和通信负担容易造成单点故障；此外集中式算法易受通信故障影响导致无法正常实现经济调度功能。相比于集中式算法，分布式算法要求发电单元获取邻居机组信息局部计算出力计划，分散了计算和通信负担避免单点故障，满足即插即用功能和适应拓扑变化，因此具有更好的鲁棒性和前瞻性。
目前，电热综合能源系统优化调度领域主要集中于系统建模、风电消纳、能源侧与负荷侧的不确定性等研究，未能在能量传输过程中有效考虑传输损耗，忽略了其对系统功率平衡的重要影响，导致以往优化结果由于传输损耗的存在不能很好地满足实际负荷需求，因此在优化调度过程中存在供需不平衡问题；此外，现有电热综合能源系统经济调度文献大多数以集中式方法求解，未能对其分布式求解方法展开进一步的研究。
综上所述，有必要针对电力系统经济调度提供一种新的优化调度方法—一种电热综合能源系统分布式经济调度方法，以解决考虑网络传输损耗下电热综合能源系统经济调度求解困难、计算复杂、多约束含耦合的问题，同时达到协同优化机组多种类型出力和保障系统安全经济运行的目的。
发明内容
为了解决考虑网络传输损耗下电热综合能源系统经济调度求解困难、计算复杂、多约束含耦合的问题，本申请实施例提供一种计及网损的电热综合能源系统分布式经济调度方法及系统。
根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；
根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；
根据最小值优化模型，利用双层一致性算法得到系统运行总成本最小值。
进一步的，所述系统运行总成本最小值、所述供需平衡约束条件和所述运行限制约束条件，具体为：
根据纯发电机组运行总成本、热电联产机组运行总成本和纯产热机组运行总成本得到所述系统运行总成本最小值；
根据系统电负荷需求、纯发电机组电出力、热电联产机组的电出力和系统电传输损耗、系统热负荷需求、纯产热机组热出力、热电联产机组的热出力和系统热传输损耗得到所述供需平衡约束条件；
根据纯发电机组电出力上下限、纯产热机组热出力上下限、热电联产机组热-电可运行域、电网线路传输功率上下限、热网管道供水温度上下限、热网管道传输流量上下限、热网管道传输热量得到所述运行限制约束条件。
进一步的，所述电热综合能源系统经济调度模型为：
其中，Ft、Fp、Fc和Fh分别表示系统运行总成本、纯发电机组运行总成本、热电联产机组运行总成本和纯产热机组运行总成本，fi(Pi)、fj(Pj,Hj)和fk(Hk)分别表示第i个纯发电机组运行成本函数、第j个热电联产机组运行成本函数和第k个纯产热机组运行成本函数。
进一步的，所述最小值优化模型，具体为：
min L＝Ft-λpΔp-λhΔh；
其中，minL为优化模型最小值，λp和λh分别表示电出力等式约束和热出力等式约束的拉格朗日乘子，Δp和Δh分别表示系统电功率偏差和系统热功率偏差。
进一步的，根据最小值优化模型，利用双层一致性算法得到满足供需平衡约束条件与运行限制约束条件的各个机组出力值，根据各个机组出力值得到系统运行总成本最小值。
一种计及网损的电热综合能源系统分布式经济调度系统包括：电热综合能源系统经济调度模型建立模块，根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；
最小值优化模型转化模块，根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；
系统运行总成本最小值求解模块，根据最小值优化模型，利用双层一致性算法得到满足供需平衡约束条件与运行限制约束条件的各个机组出力值，根据各个机组出力值得到系统运行总成本最小值。
由以上技术方案可知，本申请提供的一种计及网损的电热综合能源系统分布式经济调度方法及系统；根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；根据最小值优化模型，利用双层一致性算法得到系统运行总成本最小值。
在实际应用过程中，由于在优化调度求解中，系统优先调度增量成本小的机组出力以使系统运行总成本减小，同时兼顾系统约束条件，使得机组的最优出力与机组的增量成本成负相关；同时调度模型中计及了网络传输损耗和电热耦合约束条件由此得到的优化出力结果不仅能够满足用户实际负荷需求，还能减少企业的产能成本以提高经济效益；最后设计的分布式双一致性算法不仅能够很好地求解电热综合能源系统多约束、强耦合优化问题，而且迭代计算均匀分配给各个参与机组，因此对通信交互要求较低，并有效保护了参与者隐私，最终具有较快的收敛速度和满意的收敛结果。
附图说明
为了更清楚地说明本申请的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例的一种计及网损的电热综合能源系统分布式经济调度方法流程图；
图2为本申请实施例的系统结构示意图；
图3为本申请实施例的双一致性算法流程图；
图4为本申请实施例的机组通信拓扑图；
图5为本申请实施例不考虑网络传输损耗下的电出力一致性变量波形图；
图6为本申请实施例不考虑网络传输损耗下的热出力一致性变量波形图；
图7为本申请实施例不考虑网络传输损耗下的机组出力波形图；
图8为本申请实施例不考虑网络传输损耗下的功率偏差波形图；
图9为本申请实施例考虑网络传输损耗下的电出力一致性变量波形图；
图10为本申请实施例考虑网络传输损耗下的热出力一致性变量波形图；
图11为本申请实施例考虑网络传输损耗下的机组出力波形图；
图12为本申请实施例考虑网络传输损耗下的功率偏差波形图。
具体实施方式
为了解决现有技术中考虑网络传输损耗下电热综合能源系统经济调度求解困难、计算复杂、多约束含耦合的问题，本申请实施例提供一种计及网损的电热综合能源系统分布式经济调度方法及系统。
参见图1，为本申请实施例的一种计及网损的电热综合能源系统分布式经济调度方法流程图，所述一种计及网损的电热综合能源系统分布式经济调度方法包括：根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；根据最小值优化模型，利用双层一致性算法得到系统运行总成本最小值。
进一步的，以所述电热综合能源系统经济调度模型中的所述系统以运行总成本最小值为目标，具体为：
建立电热综合能源系统经济调度模型，包括供需平衡约束条件、运行限制约束条件和系统运行总成本最小值。
建立电热综合能源系统经济调度模型的系统以运行总成本模型：如图2所示本申请实施例的系统结构示意图，系统包含纯发电机组总数为Np编号i＝1,2,3,…,Np，热电联产机组总数为Nc编号j＝1,2,3,…,Nc，以及纯产热机组总数为Nh编号k＝1,2,3,…,Nh；Pi表示第i个纯发电机组的电出力，Pj和Hj分别表示第j个热电联产机组的电出力和热出力，Hk表示第k个纯产热机组的热出力，系统运行总成本，具体为：
其中，Ft、Fp、Fc和Fh分别表示系统运行总成本、纯发电机组运行总成本、热电联产机组运行总成本和纯产热机组运行总成本，fi(Pi)、fj(Pj,Hj)和fk(Hk)分别表示第i个纯发电机组运行成本函数、第j个热电联产机组运行成本函数和第k个纯产热机组运行成本函数，具体为：
其中，αi、βi和γi＞0表示fi(Pi)的拟合参数，αj、βj、γj＞0、δj、θj＞0和εj表示fj(Pj,Hj)的拟合参数，αk、βk和γk＞0表示fk(Hk)的拟合参数。
本申请实施例中，各类机组运行成本函数及出力限制参数取值见表1所示：
表1各类机组运行成本函数及出力限制参数
| 机组 | α | β | γ | δ | θ | ε | P<sub>i</sub><sup>m</sup>/H<sub>k</sub><sup>m</sup>(MW) | P<sub>i</sub><sup>M</sup>/H<sub>k</sub><sup>M</sup>(MW) |
| G<sub>p1</sub> | 25 | 3.0 | 0.020 | -- | -- | -- | 10 | 120 |
| G<sub>p2</sub> | 40 | 3.2 | 0.016 | -- | -- | -- | 25 | 150 |
| G<sub>p3</sub> | 75 | 2.6 | 0.018 | -- | -- | -- | 30 | 200 |
| G<sub>p4</sub> | 100 | 2.4 | 0.012 | -- | -- | -- | 40 | 300 |
| G<sub>c1</sub> | 1250 | 2.2 | 0.032 | 1.2 | 0.032 | 0.008 | -- | -- |
| G<sub>c2</sub> | 680 | 1.2 | 0.048 | 0.4 | 0.044 | 0.021 | -- | -- |
| G<sub>h1</sub> | 650 | 1.6 | 0.036 | -- | -- | -- | 0 | 1695 |
| G<sub>h2</sub> | 520 | 1.2 | 0.024 | -- | -- | -- | 0 | 1250 |
确立电热综合能源系统经济调度模型的供需平衡约束条件，具体为：
其中，Δp和Pd分别表示系统电功率偏差与系统电负荷需求，Pl表示系统电传输损耗，具体为：
其中，Bim、Bij和Bjn表示损耗系数矩阵B中对应的元素。
本申请实施例中，系统电负荷需求取Pd＝700MW，损耗系数矩阵B取值如下：
其中，Δh和Hd分别表示系统热功率偏差与系统热负荷需求，Hl表示系统热传输损耗具体为：
其中，Ng表示热媒流经管道的总段数，la表示热媒流经管道a的长度，ts,b表示热网节点b的供水温度，te,a表示管道a中周围介质的平均温度，Rh表示热媒到周围介质间每千米管道的总热阻。
本申请实施例中，系统热负荷需求取Hd＝380MW，热力网络传输管道参数取值见表2所示：
表2热力网络传输管道参数
| 管道 | L<sub>a</sub>(km) | m<sub>a</sub><sup>m</sup>(m<sup>3</sup>/h) | m<sub>a</sub><sup>M</sup>(m<sup>3</sup>/h) | R<sub>h</sub>(km×℃/kW) | 节点 | t<sub>s,b</sub><sup>m</sup>(℃) | t<sub>s,b</sub><sup>M</sup>(℃) |
| 5-12 | 2.8 | 0 | 3000 | 20 | 5 | 80 | 100 |
| 6-12 | 2.5 | 0 | 3000 | 20 | 6 | 80 | 100 |
| 7-12 | 3.0 | 0 | 3000 | 20 | 7 | 80 | 100 |
| 8-12 | 2.6 | 0 | 3000 | 20 | 8 | 80 | 100 |
确立电热综合能源系统经济调度模型的运行限制约束条件。
Pi m≤Pi≤Pi M；
其中，Pi M和Pi m分别表示第i个纯发电机组电出力的上限和下限。
其中，
和
分别表示第k个纯产热机组热出力的上限和下限。
其中，
和
表示热电联产机组的热-电可运行域参数。
本申请实施例中，热电联产机组热-电可运行域参数取值见表3所示：
表3热电联产机组热-电可运行域参数
| 机组 | 热电可运行域(H<sub>j</sub>,P<sub>j</sub>) |
| G<sub>c1</sub> | A<sub>1</sub>(0,187),B<sub>1</sub>(153,132),C<sub>1</sub>(121,42),D<sub>1</sub>(0,63) |
| G<sub>c2</sub> | A<sub>2</sub>(0,94),B<sub>2</sub>(122,68),C<sub>2</sub>(106,22),D<sub>2</sub>(0,36) |
其中，
和
分别表示第f条电网线路传输功率Pt,f的上限和下限。
本申请实施例中，电力网络传输线路参数取值见表4所示：
表4电力网络传输线路参数
| 线路 | P<sub>t,f</sub><sup>m</sup>(MW) | P<sub>t,f</sub><sup>M</sup>(MW) | 线路 | P<sub>t,f</sub><sup>m</sup>(MW) | P<sub>t,f</sub><sup>M</sup>(MW) |
| 1-11 | 0 | 140 | 2-11 | 0 | 170 |
| 3-11 | 0 | 220 | 4-11 | 0 | 320 |
| 5-11 | 0 | 160 | 6-11 | 0 | 100 |
Hb＝Cma(ts,b-tr,b)；
其中，
和
分别表示ts,b的上限和下限，
和
分别表示热网管道a传输流量ma的上限和下限，Hb表示热网节点b的传输热量，tr,b表示热网节点b的回水温度，C表示热媒的比热容。
本申请实施例中，周围介质的平均温度均取te,a＝0℃，热网节点的回水温度均取tr,b＝45℃，热媒比热容取C＝4.2kJ/(kg×℃)。
根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型。
所述最小值优化模型，具体为：
minL＝Ft-λpΔp-λhΔh；
其中，minL为最小调度值，λp和λh分别表示电出力等式约束和热出力等式约束的拉格朗日乘子，Δp和Δh分别表示系统电功率偏差和系统热功率偏差。
考虑网络传输损耗和不等式约束，根据Pi、Pj、Hj和Hk得到系统运行总成本最小值的最优条件，具体为：
其中，fi p和
分别表示第i个纯发电机组和第j个热电联产机组的电传输损耗惩罚因子，
和
分别表示第j个热电联产机组和第k个纯产热机组的热传输损耗惩罚因子。具体为：
根据最小值优化模型，利用双层一致性算法得到系统运行总成本最小值。如图3为本申请实施例的双一致性算法流程图所示。
输入电热综合能源系统相关参数，包括纯发电机组运行成本函数拟合参数αi、βi和γi，热电联产机组运行成本函数拟合参数αj、βj、γj、δj、θj和εj，纯产热机组运行成本函数拟合参数αk、βk和γk，热网传输管道参数Ng、la和Rh，损耗系数矩阵B，纯发电机组出力上下限约束参数Pi M和Pi m，纯产热机组出力上下限约束参数
和
热电联产机组热-电可运行域约束参数，系统电负荷需求Pd，系统热负荷需求Hd。
设置迭代次数τ＝0,1,2...，当τ＝0时设置各机组出力初值并使其满足如下：
本申请实施例中，各机组出力初值取值如下：
[P1 P2 P3 P4 P5 H5 P6 H6 H7 H8]＝[70 100 150 200 110 100 70 80 90 110]。
测量热网节点的供水温度ts,b[τ]和管道周围介质的平均温度te,a[τ]，并根据系统点功率损耗和系统热功率损耗分别转化系统电传输损耗Pl[τ]和系统热传输损耗Hl[τ]，进一步根据第i个纯发电机组、第j个热电联产机组的电传输损耗惩罚因子、第j个热电联产机组和第k个纯产热机组的热传输损耗惩罚因子分别转化传输损耗惩罚因子fi p[τ]、
和fk h[τ]。
选取各机组增量成本作为一致性变量，具体为：
其中，
和
分别表示第i个纯发电机组和第j个热电联产机组的电出力一致性变量，
和
分别表示第j个热电联产机组和第k个纯产热机组的热出力一致性变量。
选取任意一热电联产机组作为领导者节点，剩余所有机组均作为跟随者节点，进一步更新各节点一致性变量如下：
其中，ωp和ωh分别表示系统电功率偏差和系统热功率偏差的收敛因子，且ωp,ωh∈(0,1)，
分别表示电力子系统、热力子系统梅特罗波利斯加权矩阵Qp、Qh中对应的元素，且矩阵Q由机组通信拓扑图决定。具体为：
其中，dm和dn分别表示节点m和节点n的度，Nm表示节点m的邻居节点集合；
本申请实施例中，选取热电联产机组Gc1作为领导者节点，其余机组作为跟随者节点，功率偏差收敛因子取ωp＝ωh＝0.001，由图4机组通信拓扑图决定的矩阵Qp、Qh如下所示：
根据更新后的节点一致性变量，求解满足约束条件的机组出力，具体描述如下：
其中，Ωp＝{i|Pi＝Pi M∪Pi＝Pi m}表示电出力达到上限或下限的纯发电机组集合。
其中，
表示电出力达到可运行域边界的热电联产机组集合。
其中，
表示热出力达到可运行域边界的热电联产机组集合。
其中，
表示热出力达到上限或下限的纯产热机组集合。
将上述Pl[τ]、Pi[τ+1]、Pj[τ+1]、Hl[τ]、Hj[τ+1]、Hk[τ+1]根据系统电功率偏差和系统热功率偏差分别转化为系统电功率偏差Δp[τ+1]和系统热功率偏差Δh[τ+1]。
判断系统功率偏差是否满足收敛条件，具体为：
μ≥max(|Δp[τ+1]|,|Δh[τ+1]|)；
其中，μ表示收敛判定系数，即不小于功率偏差绝对值的最大值。
若不满足收敛，将τ＝τ+1返回转化传输损耗惩罚因子fi p[τ]、
和fk h[τ]；否则，输出当前迭代时刻的机组出力Pi[τ+1]、Pj[τ+1]、Hj[τ+1]和Hk[τ+1]，并根据系统运行总成本转化为系统运行总成本最小值Ft *。
本申请实施例中，收敛判定系数取值μ＝0.002。
为说明所提求解算法的有效性，本实施例通过以下2个算例进行验证，仿真平台采用Matlab运行实现，算例仿真结果见表6-7所示：
表6机组最优出力(单位：MW)
表7系统最小运行总成本
算例1：不考虑网络传输损耗下分布式经济调度策略的有效性。本算例忽略了网络传输损耗，通过快速迭代，机组电出力、热出力一致性变量收敛为λp＝5.0779和λh＝4.5524，最终在计及网络传输损耗的前提下系统达到了供需平衡，仿真波形如图5为本申请实施例不考虑网络传输损耗下的电出力一致性变量波形图、图6为本申请实施例不考虑网络传输损耗下的热出力一致性变量波形图、图7为本申请实施例不考虑网络传输损耗下的机组出力波形图和图8为本申请实施例不考虑网络传输损耗下的功率偏差波形图。
算例2：考虑网络传输损耗下分布式经济调度策略的有效性。本算例考虑了网络传输损耗，通过快速迭代，机组电出力、热出力一致性变量收敛为λp＝5.2648和λh＝4.5637，网络传输损耗为为Pl＝10.3369MW和Hl＝0.2888MW，最终在计及网络传输损耗的前提下系统达到了供需平衡，仿真波形如图9为本申请实施例考虑网络传输损耗下的电出力一致性变量波形图、图10为本申请实施例考虑网络传输损耗下的热出力一致性变量波形图、图11为本申请实施例考虑网络传输损耗下的机组出力波形图、图12为本申请实施例考虑网络传输损耗下的功率偏差波形图。
在本申请的一种实施例中，一种计及网损的电热综合能源系统分布式经济调度系统，包括：电热综合能源系统经济调度模型建立模块，根据电热综合能源系统的系统供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；最小值优化模型转化模块，根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；系统运行总成本最小值求解模块，根据最小值优化模型，利用双层一致性算法得到满足供需平衡约束条件与运行限制约束条件的各个机组出力值，根据各个机组出力值得到系统运行总成本最小值。
由以上技术方案可知，所述一种计及网损的电热综合能源系统分布式经济调度方法包括：根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；根据最小值优化模型，利用双层一致性算法得到系统运行总成本最小值。
在实际应用过程中，由于在优化调度求解中，系统优先调度增量成本小的机组出力以使系统运行总成本减小，同时兼顾系统约束条件，使得机组的最优出力与机组的增量成本成负相关；同时调度模型中计及了网络传输损耗和电热耦合约束条件由此得到的优化出力结果不仅能够满足用户实际负荷需求，还能减少企业的产能成本以提高经济效益；最后设计的分布式双一致性算法不仅能够很好地求解电热综合能源系统多约束、强耦合优化问题，而且迭代计算均匀分配给各个参与机组，因此对通信交互要求较低，并有效保护了参与者隐私，最终具有较快的收敛速度和满意的收敛结果。
以上结合具体实施方式和范例性实例对本申请进行了详细说明，不过这些说明并不能理解为对本申请的限制。本领域技术人员理解，在不偏离本申请精神和范围的情况下，可以对本申请技术方案及其实施方式进行多种等价替换、修饰或改进，这些均落入本申请的范围内。
Claims (6)
1.一种计及网损的电热综合能源系统分布式经济调度方法，其特征在于，包括：
根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；
根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；
根据最小值优化模型，利用双层一致性算法得到系统运行总成本最小值。
2.根据权利要求1所述的方法，其特征在于，所述系统运行总成本最小值、所述供需平衡约束条件和所述运行限制约束条件，具体为：
根据纯发电机组运行总成本、热电联产机组运行总成本和纯产热机组运行总成本得到所述系统运行总成本最小值；
根据系统电负荷需求、纯发电机组电出力、热电联产机组的电出力和系统电传输损耗、系统热负荷需求、纯产热机组热出力、热电联产机组的热出力和系统热传输损耗得到所述供需平衡约束条件；
根据纯发电机组电出力上下限、纯产热机组热出力上下限、热电联产机组热-电可运行域、电网线路传输功率上下限、热网管道供水温度上下限、热网管道传输流量上下限、热网管道传输热量得到所述运行限制约束条件。
3.根据权利要求1所述的方法，其特征在于，所述电热综合能源系统经济调度模型为：
其中，Ft、Fp、Fc和Fh分别表示系统运行总成本、纯发电机组运行总成本、热电联产机组运行总成本和纯产热机组运行总成本，fi(Pi)、fj(Pj,Hj)和fk(Hk)分别表示第i个纯发电机组运行成本函数、第j个热电联产机组运行成本函数和第k个纯产热机组运行成本函数。
4.根据权利要求1所述的方法，其特征在于，所述最小值优化模型，具体为：
min L＝Ft-λpΔp-λhΔh；
其中，min L为优化模型最小值，λp和λh分别表示电出力等式约束和热出力等式约束的拉格朗日乘子，Δp和Δh分别表示系统电功率偏差和系统热功率偏差。
5.根据权利要求1所述的方法，其特征在于，根据最小值优化模型，利用双层一致性算法得到满足供需平衡约束条件与运行限制约束条件的各个机组出力值，根据各个机组出力值得到系统运行总成本最小值。
6.一种计及网损的电热综合能源系统分布式经济调度系统，其特征在于，包括：
电热综合能源系统经济调度模型建立模块，根据电热综合能源系统的供需平衡约束条件和运行限制约束条件，以电热综合能源系统的系统运行总成本最小值为目标，建立电热综合能源系统经济调度模型；最小值优化模型转化模块，根据所述电热综合能源系统经济调度模型，利用拉格朗日乘子法转化为最小值优化模型；系统运行总成本最小值求解模块，根据最小值优化模型，利用双层一致性算法得到满足供需平衡约束条件与运行限制约束条件的各个机组出力值，根据各个机组出力值得到系统运行总成本最小值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111320160.0A CN114021997A (zh) | 2021-11-09 | 2021-11-09 | 计及网损的电热综合能源系统分布式经济调度方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111320160.0A CN114021997A (zh) | 2021-11-09 | 2021-11-09 | 计及网损的电热综合能源系统分布式经济调度方法及系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114021997A true CN114021997A (zh) | 2022-02-08 |
Family
ID=80062859
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111320160.0A Pending CN114021997A (zh) | 2021-11-09 | 2021-11-09 | 计及网损的电热综合能源系统分布式经济调度方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114021997A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114640138A (zh) * | 2022-05-18 | 2022-06-17 | 浙江大学 | 一种凸包经济运行域的求解方法及装置、电子设备 |
-
2021
- 2021-11-09 CN CN202111320160.0A patent/CN114021997A/zh active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114640138A (zh) * | 2022-05-18 | 2022-06-17 | 浙江大学 | 一种凸包经济运行域的求解方法及装置、电子设备 |
| CN114640138B (zh) * | 2022-05-18 | 2022-09-02 | 浙江大学 | 一种凸包经济运行域的求解方法及装置、电子设备 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110245878B (zh) | 智能楼宇群体的分布式综合能源需求响应协同优化方法 | |
| CN109345045B (zh) | 基于双乘子迭代算法的电热综合能源系统经济调度方法 | |
| Lu et al. | Thermal inertial aggregation model for integrated energy systems | |
| Lin et al. | Decentralized solution for combined heat and power dispatch through benders decomposition | |
| CN110210747B (zh) | 一种电热气互联能源系统灵活性调度方法 | |
| CN106447113B (zh) | 一种基于运行优化模型的多区域综合能源系统运行方法 | |
| CN108921727A (zh) | 考虑热负荷动态特性的区域综合能源系统可靠性评估方法 | |
| CN109190271B (zh) | 一种考虑传输损耗的电热综合能源系统经济优化调度方法 | |
| CN109492325B (zh) | 一种基于扩展能源集线器的多能耦合系统潮流分析方法 | |
| CN109066812B (zh) | 一种基于双层一致性算法的电力系统多目标优化调度方法 | |
| CN109102104B (zh) | 一种基于碳排放流的多能源系统协同低碳规划方法 | |
| CN109711601A (zh) | 电-气-热综合能源系统分布式优化调度方法和装置 | |
| CN108898265A (zh) | 一种综合能源系统一体化规划方法 | |
| CN116341881B (zh) | 计及热网灵活性的电-热系统鲁棒超前调度方法及系统 | |
| CN114077934B (zh) | 一种综合能源微网互联系统及其调度方法 | |
| CN107392366A (zh) | 计及用户行为分析的综合能源系统供需双侧综合优化方法 | |
| CN112380681A (zh) | 基于分布式优化的综合能源系统参数确定方法及相关装置 | |
| CN111476394B (zh) | 一种适用于电热气等多能源系统的鲁棒运行优化方法 | |
| CN115392035A (zh) | 电-气互联综合能源系统优化调度的改进方法 | |
| CN114021997A (zh) | 计及网损的电热综合能源系统分布式经济调度方法及系统 | |
| CN108594658A (zh) | 一种电-气耦合系统概率最大负荷裕度多目标优化模型及其求解方法 | |
| CN113324283B (zh) | 一种计及温度半动态特性的集中供热网络等值方法 | |
| CN111310310B (zh) | 一种用于量调节的热力系统静态潮流快速解耦计算方法 | |
| CN111144642B (zh) | 一种综合能源系统优化方法及系统 | |
| CN112926835A (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 |