CN114399115A - 以多目标均衡为导向的热力控制平台用户响应优化方法 - Google Patents
以多目标均衡为导向的热力控制平台用户响应优化方法 Download PDFInfo
- Publication number
- CN114399115A CN114399115A CN202210049074.9A CN202210049074A CN114399115A CN 114399115 A CN114399115 A CN 114399115A CN 202210049074 A CN202210049074 A CN 202210049074A CN 114399115 A CN114399115 A CN 114399115A
- Authority
- CN
- China
- Prior art keywords
- thermal
- response
- load
- user
- thermal load
- 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
- 230000004044 response Effects 0.000 title claims abstract description 172
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000005457 optimization Methods 0.000 title claims abstract description 28
- 239000002245 particle Substances 0.000 claims description 39
- 238000005338 heat storage Methods 0.000 claims description 24
- 230000005611 electricity Effects 0.000 claims description 18
- 230000008901 benefit Effects 0.000 claims description 14
- 230000014509 gene expression Effects 0.000 claims description 9
- 238000010438 heat treatment Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000003828 downregulation Effects 0.000 claims description 3
- 230000003827 upregulation Effects 0.000 claims description 3
- 230000006870 function Effects 0.000 description 12
- 230000009286 beneficial effect Effects 0.000 description 9
- 230000003044 adaptive effect Effects 0.000 description 8
- 238000012546 transfer Methods 0.000 description 6
- 230000033228 biological regulation Effects 0.000 description 5
- 230000000295 complement effect Effects 0.000 description 5
- 239000003795 chemical substances by application Substances 0.000 description 4
- 238000010248 power generation Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 239000000178 monomer Substances 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000006116 polymerization reaction Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 238000009423 ventilation 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/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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or 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)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种以多目标均衡为导向的热力控制平台用户响应优化方法，属于热力控制平台用户响应技术领域，本发明通过整合分布式热力负荷资源方式使同一热力用户管控多个不同的热力负荷，提升热力负荷响应潜力，同时，通过评估不同热力负荷的用户舒适度、用户效益指标，采用费效比函数实现了多优化目标的均衡量化，考虑同一热力用户的不同热力负荷之间协同优化，提升地区不平衡电量消纳和热力负荷跨时段转移能力。
Description
技术领域
本发明属于热力控制平台用户响应技术领域，尤其涉及一种以多目标均衡为导向的热力控制平台用户响应优化方法。
背景技术
新能源发电过程中的随机性和波动性将导致电力潮流不确定性增加、发电侧灵活性降低等问题，如何提高系统灵活性，减小发电侧高比例新能源并网带来的负面影响，是建立新型电力系统亟待解决的问题。在此背景下，依靠发电侧资源维持电力电量平衡变得愈发困难，具有热储能性质的热力负荷是一种优质需求侧响应资源，在响应过程中表现出可调容量大、响应速率快、调节时间灵活、持续时间长的特点，通过负荷中断、负荷转移以及负荷调节的方式参与削峰填谷、新能源消纳等需求侧响应。考虑到部分热力负荷单体容量小、分散分布的特点，其响应潜力会随着聚合程度增长，因此通过热力控制平台整合热力资源是提升负荷响应潜力的关键，而建立合适的热力控制平台负荷调控策略能够在提高用户效益的基础上保证用户用热满意度，有效缓解电网压力。
现有研究通过市场机制或调度控制的方式，在保障电力系统的安全可靠运行的前提下提升热力负荷消纳不平衡电量的能力，例如热力负荷群的需求侧响应，热力负荷群的集中控制，热力负荷群的分布式控制等方法。然而热力负荷作为需求侧响应资源的重要组成部分，不同的热力负荷具有不同用热特性和用电特性，其中用热特性差异主要体现在热力负荷的供热对象不同，例如工商业用热的差异；用电特性差异主要体现在热力负荷种类上，例如蓄热式电锅炉和电热泵在响应特性上的差别较大。
由于参与热力控制平台的一个用户主体可能管理着多个不同种类的热力负荷，并且不同热力负荷具有各自的响应特征和用热特性，因此用户主体参与热力控制平台响应时应考虑其控制的居民用户、工业用户、商业用户等多种用热特性的热力负荷协同响应，同时考虑其控制的蓄热式电锅炉和电热泵等不同用电特性的热力负荷互补调控。并且，传统电力调度模式下的蓄热式电锅炉和电热泵等热力负荷通常以集中供热站的形式直接参与用电调度，较少考虑不同热力负荷的协同响应与互补调控，且无法利用一些容量较小的热力负荷资源如暖通空调。
发明内容
针对现有技术中的上述不足，本发明提供的一种以多目标均衡为导向的热力控制平台用户响应优化方法，解决了热力平台用户在制定响应计划时难以最大化挖掘自身响应潜力的问题，同时提高了热力平台用户的响应意愿，以及解决了用户对热力平台指令响应积极性不高的问题。
为了达到以上目的，本发明采用的技术方案为：
本方案提供一种以多目标均衡为导向的热力控制平台用户响应优化方法，包括以下步骤：
S1、由热力控制平台通过最大化利用自身的响应潜力，得到热力控制平台的总负荷指令，并向下辖的集中供热站、负荷聚合商以及虚拟电厂的热力用户公示所述总负荷指令；
S2、由所述热力用户根据所述总负荷指令，向热力控制平台提交响应信息；
S3、判断是否由热力控制平台对热力用户的用电和用热进行优化，若是，则进入步骤S4，否则，结束流程；
S4、由热力用户向热力控制平台提交次日的用电计划和用热计划；
S5、构建热力控制平台负荷响应的约束条件以及目标函数，并根据所述响应信息、用电计划以及用热计划，利用粒子群算法计算得到热力用户的最优响应计划；
S6、利用约束条件对所述最优响应计划进行校核，并根据校核结果由热力控制平台将所述最优响应计划发送至热力用户；
S7、由热力用户确认并提交所述最优响应计划至热力控制平台，完成以多目标均衡为导向的热力控制平台用户响应优化。
本发明的有益效果是：本发明提出的热力控制平台用户响应优化方法主要针对现货市场模式下，通过整合分布式热力负荷资源方式使同一热力用户管控多个不同的热力负荷，提升热力负荷响应潜力，同时，该方法通过评估不同热力负荷的用户舒适度、用户效益指标，采用费效比函数实现了多优化目标的均衡量化，考虑同一热力用户的不同热力负荷之间协同优化，提升地区不平衡电量消纳和热力负荷跨时段转移能力。
再进一步地，所述步骤S5中约束条件包括：热力控制平台热力负荷功率约束、第一热力负荷响应潜力约束、第二热力负荷响应潜力约束以及热力负荷的用电平衡约束。
上述进一步方案的有益效果是：本发明通过设计约束条件，保证热力负荷在响应优化后，得到的响应计划具有较高的可行性，且用户可以有效的执行该响应计划。
再进一步地，所述热力控制平台热力负荷功率约束包括：热力负荷最大功率约束Pmax、基准用电曲线的功率约束P0(t)以及最小负荷约束Pmin；
所述热力负荷最大功率约束Pmax和最小负荷约束Pmin的表达式如下：
Pmin≤P(t)≤Pmax或P(t)＝0
其中，P(t)表示t时段热力负荷参与平台响应后的响应曲线；
所述基准用电曲线的功率约束P0(t)的表达式如下：
0.9*min(P(t-24*n))≤P0(t)≤1.1*max(P(t-24*n)),n＝1,2,3,4,5
其中，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束，n表示历史天数。
上述进一步方案的有益效果是：本发明通过设置以上的约束条件保证热力用户申报的基准用电曲线不会与申报日前5日用户不参与平台响应的正常的用电负荷偏差太多，防止用户篡改申报的基准曲线。
再进一步地，所述第一热力负荷响应潜力约束条件的表达式如下：
其中，v-和v+分别表示t时段热力负荷参与平台响应的上调响应速率和下调响应速率，C-(t)和C+(t)分别表示t时段热力负荷参与平台响应的下调容量和上调容量，D-和D+分别表示热力负荷参与平台响应的最小持续时段和最大持续时段，Δt表示热力负荷实际参与响应的时间段，P(t)表示t时段热力负荷参与平台响应后的响应曲线，P(t-1)表示t-1时段热力负荷参与平台响应后的响应曲线，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束。
上述进一步方案的有益效果是：考虑到实际中热力负荷存在最大响应速率限制、最大调节容量限制以及最大响应持续时间限制，通过上述约束条件，确保了热力负荷在响应优化后得到的响应计划具有执行可行性。
再进一步地，所述第二热力负荷响应潜力约束条件的表达式如下：
其中，T-和T+分别表示t时段热力负荷温度T(t)变化的上下限，T(t)表示t时段热力负荷的温度，TC-和TC+表示储热装置的最小储热量和最大储热量，TC(t)表示t时段蓄热罐的蓄热量。
上述进一步方案的有益效果是：响应优化会改变热力负荷的用电计划，上述储热量约束保证了带有蓄热装置热力负荷的用热计划不会改变，用户的用热感知和满意度也不会改变；温度约束确保了不带有蓄热装置热力负荷虽然改变了负荷原本的用热计划，但不会降低用户的用热感知和满意度。
再进一步地，所述热力负荷的用电平衡约束的表达式如下：
ΔP(t)＝P(t)-P0(t)+C(t)
其中，C(t)表示由于是热力负荷无法完全响应指令时出现的松弛变量，ΔP(t)表示热力负荷实际响应曲线与基准用电曲线的差值，是负荷的净响应曲线,P(t)表示t时段热力负荷参与平台响应后的响应曲线，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束。
上述进一步方案的有益效果是：考虑到大部分热力负荷在执行负荷响应计划时会存在一定的偏差，无法完全响应指令，因此本约束通过松弛变量的方式，保证等式约束一定成立，且该变量可作为热力负荷执行响应计划的考核指标。
再进一步地，所述步骤S5中构建热力控制平台负荷响应的目标函数，其具体为：利用费效比函数构建热力负荷用电效益和用户满意度的热力控制平台负荷响应的目标函数。
上述进一步方案的有益效果是：用费效比函数实现了多优化目标的均衡量化，考虑同一热力用户的不同热力负荷之间协同优化，提升地区不平衡电量消纳和热力负荷跨时段转移能力。
再进一步地，所述热力负荷用电效益的表达式如下：
其中，i表示热力用户管控下的第i个热力负荷，I表示热力用户管控下的热力负荷数量，
表示t时段热力负荷i参与平台响应前用电量，
表示t时段热力负荷i参与平台响应后用电量，
表示t时段热力负荷i参与平台响应的平台激励，E(t)表示t时段热力用户管控下各热力负荷的总效益。
上述进一步方案的有益效果是：本发明通过热力负荷用电效益和用户满意度均衡优化的方法实现了热力控制平台负荷用户的有效整合，充分调动了热力负荷资源的响应潜力，通过不同热力负荷之间协同优化的方法实现了同一热力用户的不同热力负荷协同响应与互补调控，有效降低了地区弃风量，提高了热力负荷效益，一定程度的改善电力系统负荷峰谷差异。
再进一步地，所述用户满意度的表达式如下：
其中，i表示热力用户管控下的第i个热力负荷，I表示热力用户管控下的热力负荷数量，α表示热力负荷类型，当热力负荷具有蓄热罐时，α取0，反之α取1，Ti(t)表示热力用户管控下的热力负荷i在t时段的温度，Tmax表示用户感知的最高温度，Tmin表示用户感知的最低温度，Hi,0(t)表示热力负荷i不参与响应的基准用热，Hi(t)表示热力负荷i参与响应后的实际用热，S(t)表示t时段热力用户的满意度，其取值为热力用户管控下全部热力负荷满意度的平均值。
上述进一步方案的有益效果是：本发明通过热力负荷用电效益和用户满意度均衡优化的方法实现了热力控制平台负荷用户的有效整合，充分调动了热力负荷资源的响应潜力，通过不同热力负荷之间协同优化的方法实现了同一热力用户的不同热力负荷协同响应与互补调控，有效降低了地区弃风量，提高了热力负荷效益，一定程度的改善电力系统负荷峰谷差异。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
实施例
如图1所示，本发明提供了一种以多目标均衡为导向的热力控制平台用户响应优化方法，其实现方法如下：
S1、由热力控制平台通过最大化利用自身的响应潜力，得到热力控制平台的总负荷指令，并向下辖的集中供热站、负荷聚合商以及虚拟电厂的热力用户公示所述总负荷指令；
本实施例中，热力控制平台通过参与电力调度、现货市场、辅助服务等方式最大化利用自身的响应潜力，并得到热力控制平台的总负荷指令，同时向下辖的集中供热站、负荷聚合商以及虚拟电厂三类热力用户公示该指令，通过完善热力控制平台用户响应优化策略以充分调动热力负荷响应潜力，有效提升热力用户主体的不平衡电量消纳能力和负荷跨时段转移能力。
S2、由所述热力用户根据所述总负荷指令，向热力控制平台提交响应信息；
本实施例中，热力控制平台管控下的三类热力用户需要向平台补全、修改并提交相关信息，包括响应速率、可调容量、持续时间、储热量、温度阈值等。用户未提交的部分信息由平台根据用户的历史信息、数据补全。
S3、判断是否由热力控制平台对热力用户的用电和用热进行优化，若是，则进入步骤S4，否则，结束流程；
本实施例中，参与热力控制平台的各类热力用户根据自身情况和热力控制平台指令，通过自己的决策人员制定并提交响应计划，或委托热力控制平台对其热力负荷用电和用热进行优化，代为制定并提交响应计划。
S4、由热力用户向热力控制平台提交次日的用电计划和用热计划；
本实施例中，用户选择委托热力控制平台制定响应计划，则需要向平台提交次日的用电计划和用热计划。
S5、构建热力控制平台负荷响应的约束条件以及目标函数，并根据所述响应信息、用电计划以及用热计划，利用粒子群算法计算得到热力用户的最优响应计划；
本实施例中，约束条件包括：热力控制平台热力负荷功率约束、第一热力负荷响应潜力约束、第二热力负荷响应潜力约束以及热力负荷的用电平衡约束。
本实施例中，热力负荷功率约束包含对热力负荷功率的三种基本电气约束：热力负荷最大功率约束Pmax、基准用电曲线的功率约束P0(t)以及最小负荷约束Pmin，其中，热力负荷最大功率约束Pmax和最小负荷约束Pmin如下所示：
其中，P(t)表示t时段热力负荷参与平台响应后的响应曲线，Pmax是热力负荷的最大功率，Pmin表示热力负荷的最小功率。
通常情况下，对于可削减热力负荷，其功率应保证在0和最大值之间连续可调，但实际中受到最小负荷限制，热力负荷应在最小负荷和最大负荷之间连续可调，当热力负荷具有可中断特性时，热力负荷可以通过负荷中断取到0，因此上式可以改写为：
Pmin≤P(t)≤Pmax或P(t)＝0 (2)
建立历史功率约束P0(t)如下，保证热力用户申报的基准用电曲线不会与申报日前5日不参与平台响应的正常的用电负荷偏差太多。主要目的是防止用户通过虚报基准曲线来增大与响应曲线的差值。
0.9*min(P(t-24*n))≤P0(t)≤1.1*max(P(t-24*n)),n＝1,2,3,4,5 (3)
其中，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束，n表示历史天数。
本实施例中，第一热力负荷响应潜力约束包含对热力负荷响应潜力的三种基本电气约束：响应速率v、可调容量C、持续时间D。
其中，v-和v+分别表示t时段热力负荷参与平台响应的上调响应速率和下调响应速率，C-(t)和C+(t)分别表示t时段热力负荷参与平台响应的下调容量和上调容量，D-和D+分别表示热力负荷参与平台响应的最小持续时段和最大持续时段，Δt表示热力负荷实际参与响应的时间段，P(t)表示t时段热力负荷参与平台响应后的响应曲线，P(t-1)表示t-1时段热力负荷参与平台响应后的响应曲线，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束。
本实施例中，热力负荷响应潜力约束，包含对热力负荷响应潜力的两种基本热力约束：温度阈值T、蓄热罐容量TC。
其中，T-和T+分别表示t时段热力负荷温度T(t)变化的上下限，T(t)表示t时段热力负荷的温度，TC-和TC+表示储热装置的最小储热量和最大储热量，TC(t)表示t时段蓄热罐的蓄热量。
本实施例中，第二热力负荷的用电平衡约束：
ΔP(t)＝P(t)-P0(t)+C(t) (6)
其中，C(t)表示由于是热力负荷无法完全响应指令时出现的松弛变量，ΔP(t)表示热力负荷实际响应曲线与基准用电曲线的差值，P(t)表示t时段热力负荷参与平台响应后的响应曲线。
本实施例中，采用费效比函数构造热力负荷用电效益和用户满意度的综合目标函数，用电效益S(t)和用户满意度E(t)可以分别由公式(8)和(9)确定。其中，费效比函数法构造响应优化的目标函数如下所示：
其中，Si(t)表示t时段用户满意度，Ei(t)表示t时段热力负荷的用户效益。
本实施例中，热力负荷用户满意度计算公式，当热力负荷类型为不配置储热装置的温控负荷时，热力负荷用户的满意度与用户温度变化相关；当热力负荷配置储热装置时，热力负荷用户的满意度与实际用热量和计划用热量的差值相关。
其中，i表示热力用户管控下的第i个热力负荷，I表示热力用户管控下的热力负荷数量，α表示热力负荷类型，当热力负荷具有蓄热罐时，α取0，反之α取1，Ti(t)表示热力用户管控下的热力负荷i在t时段的温度，Tmax表示用户感知的最高温度，Tmin表示用户感知的最低温度，Hi,0(t)表示热力负荷i不参与响应的基准用热，Hi(t)表示热力负荷i参与响应后的实际用热，S(t)表示t时段热力用户的满意度，其取值为热力用户管控下全部热力负荷满意度的平均值。
本实施例中，热力负荷效益指标计算公式如下：
其中，i表示热力用户管控下的第i个热力负荷，I表示热力用户管控下的热力负荷数量，
表示t时段热力负荷i参与平台响应前用电量，
表示t时段热力负荷i参与平台响应后用电量，
表示t时段热力负荷i参与平台响应的平台激励，E(t)表示t时段热力用户管控下各热力负荷的总效益。
本实施例中，粒子群算法属于启发式算法的一种，该方法根据初始设定的随机解通过不断进行迭代寻找最优解，在热力控制平台用户响应优化的过程中，每个粒子可以向其他粒子传递信息，让各个粒子知道自己的位置，通过这样的信息传递来判断每个粒子的位置是不是最优解。在一次迭代中，粒子通过追踪两个“极值”来更新自己，第一个是粒子本身寻找到的最优解，称为个体极值pBest；另一个极值是整个种群目前找到的极值，称为全局极值gBest。
首先，粒子具有两个基本属性：粒子位置和粒子速度，其中，粒子的位置表示粒子的移动方向，粒子的速度表示粒子移动的快慢。在n维空间中粒子的位置表示为空间矢量Xi，粒子的速度表示为速度矢量Vi，表示形式分别如公式(10)、(11)所示：
Xi＝(x1,x2,…xn) (10)
Vi＝(v1,v2,…vn) (11)
每一个粒子所在的位置不同其由目标函数所决定的适应值也不同，储存截止到最新一次迭代时的最好位置pBest和当前位置作为本粒子的飞行经验，整个粒子群中的最好位置gBest作为同伴经验，在下一次迭代时由飞行经验和同伴经验来决定下一步运动，在每一次迭代过程中，粒子通过追踪两个极值(pBest,gBest)对位置进行更新，速度更新方法的标准形式如公式(12)和公式(13)所示：
Vi＝Vi-1+c1×rand×(pBesti-Xi)+c2×rand×(gBesti-Xi) (12)
Xi＝Xi+Vi (13)
其中，Vi表示粒子在本次迭代中的速度，其最大值为Vmax，若迭代时出现Vi≥Vmax则取Vi＝Vmax，rand表示介于(0,1)之间的随机数，Xi表示粒子在本次迭代中的位置，c1和c2均表示学习因子，通常取c1＝c2＝2。
公式(12)中的第一项为上次迭代的速度大小和方向，称为记忆项；第二项是从当前位置指向自身最好位置pBest的矢量，用于表示粒子自身经验对于速度的修正程度，成为自身认知项；第三项是从当前位置指向种群最好位置gBest的矢量，用于表示粒子之间的信息交流对于粒子速度的修正程度，称为群体认知项。
粒子群算法的具体流程如下：
(1)设定种群规模为n，初始化粒子群的随机位置和速度；
(2)计算每个粒子的适应度；
(3)对于每一个粒子，判断其适应值和经过最好位置pBest进行比较，如果较好，则将其作为当前最好位置pBest；
(4)对于每一个粒子，判断其适应值和经过最好位置gBest进行比较，如果较好，则将其作为当前最好位置gBest；
(5)根据速度和位置更新公式对粒子的速度和位置进行更新；
(6)判断迭代结果是否满足设定的收敛条件，若满足则输出最后结果，得到热力用户的最优响应计划，若不满足则转到第(2)步再次进行迭代。
S6、利用约束条件对所述最优响应计划进行校核，并根据校核结果由热力控制平台将所述最优响应计划发送至热力用户；
S7、由热力用户确认并提交所述最优响应计划至热力控制平台，完成以多目标均衡为导向的热力控制平台用户响应优化。
本实施例中，热力控制平台利用约束条件对最优响应计划进行校核，确保热力用户可以正常执行优化后的响应计划，热力控制平台将该响应计划发还给热力用户，由热力用户确认并提交该响应计划到热力控制平台。
本发明通过以上设计，整合分布式热力负荷资源方式使同一热力用户管控多个不同的热力负荷，提升热力负荷响应潜力，同时，该方法通过评估不同热力负荷的用户舒适度、用户效益指标，采用费效比函数实现了多优化目标的均衡量化，考虑同一热力用户的不同热力负荷之间协同优化，提升地区不平衡电量消纳和热力负荷跨时段转移能力。
Claims (9)
1.一种以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，包括以下步骤：
S1、由热力控制平台通过最大化利用自身的响应潜力，得到热力控制平台的总负荷指令，并向下辖的集中供热站、负荷聚合商以及虚拟电厂的热力用户公示所述总负荷指令；
S2、由所述热力用户根据所述总负荷指令，向热力控制平台提交响应信息；
S3、判断是否由热力控制平台对热力用户的用电和用热进行优化，若是，则进入步骤S4，否则，结束流程；
S4、由热力用户向热力控制平台提交次日的用电计划和用热计划；
S5、构建热力控制平台负荷响应的约束条件以及目标函数，并根据所述响应信息、用电计划以及用热计划，利用粒子群算法计算得到热力用户的最优响应计划；
S6、利用约束条件对所述最优响应计划进行校核，并根据校核结果由热力控制平台将所述最优响应计划发送至热力用户；
S7、由热力用户确认并提交所述最优响应计划至热力控制平台，完成以多目标均衡为导向的热力控制平台用户响应优化。
2.根据权利要求1所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述步骤S5中约束条件包括：热力控制平台热力负荷功率约束、第一热力负荷响应潜力约束、第二热力负荷响应潜力约束以及热力负荷的用电平衡约束。
3.根据权利要求2所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述热力控制平台热力负荷功率约束包括：热力负荷最大功率约束Pmax、基准用电曲线的功率约束P0(t)以及最小负荷约束Pmin；
所述热力负荷最大功率约束Pmax和最小负荷约束Pmin的表达式如下：
Pmin≤P(t)≤Pmax或P(t)＝0
其中，P(t)表示t时段热力负荷参与平台响应后的响应曲线；
所述基准用电曲线的功率约束P0(t)的表达式如下：
0.9*min(P(t-24*n))≤P0(t)≤1.1*max(P(t-24*n)),n＝1,2,3,4,5
其中，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束，n表示历史天数。
4.根据权利要求2所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述第一热力负荷响应潜力约束的表达式如下：
其中，v-和v+分别表示t时段热力负荷参与平台响应的上调响应速率和下调响应速率，C-(t)和C+(t)分别表示t时段热力负荷参与平台响应的下调容量和上调容量，D-和D+分别表示热力负荷参与平台响应的最小持续时段和最大持续时段，Δt表示热力负荷实际参与响应的时间段，P(t)表示t时段热力负荷参与平台响应后的响应曲线，P(t-1)表示t-1时段热力负荷参与平台响应后的响应曲线，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束。
5.根据权利要求2所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述第二热力负荷响应潜力约束条件的表达式如下：
其中，T-和T+分别表示t时段热力负荷温度T(t)变化的上下限，T(t)表示t时段热力负荷的温度，TC-和TC+表示储热装置的最小储热量和最大储热量，TC(t)表示t时段蓄热罐的蓄热量。
6.根据权利要求2所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述热力负荷的用电平衡约束的表达式如下：
ΔP(t)＝P(t)-P0(t)+C(t)
其中，C(t)表示由于是热力负荷无法完全响应指令时出现的松弛变量，ΔP(t)表示热力负荷实际响应曲线与基准用电曲线的差值，是负荷的净响应曲线,P(t)表示t时段热力负荷参与平台响应后的响应曲线，P0(t)表示t时段热力负荷参与平台响应前的基准用电曲线的功率约束。
7.根据权利要求1所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述步骤S5中构建热力控制平台负荷响应的目标函数，其具体为：利用费效比函数构建热力负荷用电效益和用户满意度的热力控制平台负荷响应的目标函数。
8.根据权利要求7所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述热力负荷用电效益的表达式如下：
其中，i表示热力用户管控下的第i个热力负荷，I表示热力用户管控下的热力负荷数量，
表示t时段热力负荷i参与平台响应前用电量，
表示t时段热力负荷i参与平台响应后用电量，
表示t时段热力负荷i参与平台响应的平台激励，E(t)表示t时段热力用户管控下各热力负荷的总效益。
9.根据权利要求7所述的以多目标均衡为导向的热力控制平台用户响应优化方法，其特征在于，所述用户满意度的表达式如下：
其中，i表示热力用户管控下的第i个热力负荷，I表示热力用户管控下的热力负荷数量，α表示热力负荷类型，当热力负荷具有蓄热罐时，α取0，反之α取1，Ti(t)表示热力用户管控下的热力负荷i在t时段的温度，Tmax表示用户感知的最高温度，Tmin表示用户感知的最低温度，Hi,0(t)表示热力负荷i不参与响应的基准用热，Hi(t)表示热力负荷i参与响应后的实际用热，S(t)表示t时段热力用户的满意度，其取值为热力用户管控下全部热力负荷满意度的平均值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210049074.9A CN114399115A (zh) | 2022-01-17 | 2022-01-17 | 以多目标均衡为导向的热力控制平台用户响应优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210049074.9A CN114399115A (zh) | 2022-01-17 | 2022-01-17 | 以多目标均衡为导向的热力控制平台用户响应优化方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114399115A true CN114399115A (zh) | 2022-04-26 |
Family
ID=81231211
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210049074.9A Pending CN114399115A (zh) | 2022-01-17 | 2022-01-17 | 以多目标均衡为导向的热力控制平台用户响应优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114399115A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114612021A (zh) * | 2022-05-12 | 2022-06-10 | 四川大学 | 计及多粒度属性的热力负荷协同调控方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108599144A (zh) * | 2018-03-22 | 2018-09-28 | 国网天津市电力公司 | 一种提高电力系统新能源利用率且最小化发电成本的方法 |
| US20200209901A1 (en) * | 2018-01-04 | 2020-07-02 | State Grid Jiangsu Electric Power Co., Ltd | User side load response method based on adjustment and control on temperature of load clusters |
| CN113283709A (zh) * | 2021-04-29 | 2021-08-20 | 国网浙江省电力有限公司湖州供电公司 | 一种面向农村电网的分层分布式负荷协同调度方法 |
| CN113809755A (zh) * | 2021-08-19 | 2021-12-17 | 中国华电集团有限公司福建分公司 | 一种基于需求响应的智慧楼宇节能优化控制方法 |
-
2022
- 2022-01-17 CN CN202210049074.9A patent/CN114399115A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200209901A1 (en) * | 2018-01-04 | 2020-07-02 | State Grid Jiangsu Electric Power Co., Ltd | User side load response method based on adjustment and control on temperature of load clusters |
| CN108599144A (zh) * | 2018-03-22 | 2018-09-28 | 国网天津市电力公司 | 一种提高电力系统新能源利用率且最小化发电成本的方法 |
| CN113283709A (zh) * | 2021-04-29 | 2021-08-20 | 国网浙江省电力有限公司湖州供电公司 | 一种面向农村电网的分层分布式负荷协同调度方法 |
| CN113809755A (zh) * | 2021-08-19 | 2021-12-17 | 中国华电集团有限公司福建分公司 | 一种基于需求响应的智慧楼宇节能优化控制方法 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114612021A (zh) * | 2022-05-12 | 2022-06-10 | 四川大学 | 计及多粒度属性的热力负荷协同调控方法 |
| CN114612021B (zh) * | 2022-05-12 | 2022-09-13 | 四川大学 | 计及多粒度属性的热力负荷协同调控方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112529256B (zh) | 考虑多重不确定性的分布式电源集群日前调度方法和系统 | |
| CN110044020B (zh) | 计及空调用户舒适度的需求侧响应方法 | |
| CN108875992A (zh) | 一种考虑需求响应的虚拟电厂日前优化调度方法 | |
| CN108494012B (zh) | 一种计及电转气技术的区域综合能源系统在线优化方法 | |
| CN113991655A (zh) | 定频空调负荷聚合需求响应潜力评估方法、装置及介质 | |
| CN107706932A (zh) | 一种基于动态自适应模糊逻辑控制器的能量调度优化方法 | |
| CN112348283A (zh) | 蓄热式电采暖虚拟电厂日前可调度潜力评估方法及装置 | |
| CN112734116A (zh) | 一种含分布式储能的主动配电网优化调度方法 | |
| CN114399115A (zh) | 以多目标均衡为导向的热力控制平台用户响应优化方法 | |
| Chen et al. | Day-ahead scheduling of large numbers of thermostatically controlled loads based on equivalent energy storage model | |
| CN114239372A (zh) | 一种考虑机组组合的多目标机组检修双层优化方法和系统 | |
| CN117559526A (zh) | 一种基于光储充一体化充电站的拟路由器能量调控方法 | |
| CN114825371B (zh) | 基于调节前后节点电压约束的聚合温控负荷多层调控方法 | |
| CN115549111A (zh) | 一种面向微电网的温控负荷集群控制方法、系统及介质 | |
| CN114612021B (zh) | 计及多粒度属性的热力负荷协同调控方法 | |
| CN115147245A (zh) | 一种工业负荷参与调峰辅助服务的虚拟电厂优化调度方法 | |
| CN115438906A (zh) | 一种多虚拟电厂点对点交易方法、电子设备及存储介质 | |
| CN114123175B (zh) | 考虑广义需求侧资源联盟特性的配电网分布式控制方法 | |
| Xing et al. | A rolling optimization method of reserve capacity considering wind power frequency control | |
| CN115526504A (zh) | 泵站供水系统节能调度方法、系统、电子设备及存储介质 | |
| CN112862282B (zh) | 一种基于ddqn算法的综合能源系统源-荷协同调度优化方法 | |
| CN115051374A (zh) | 电采暖设备参与电力调峰的控制方法、装置和存储介质 | |
| CN108321851A (zh) | 基于分层递阶控制的风电集群并网有功功率优化调度方法 | |
| CN113078629A (zh) | 一种用于集群温控负荷聚合商功率调控的聚合体功率分配模型及分布式一致性控制方法 | |
| CN117833372B (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 |