CN115566689A - 一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法 - Google Patents
一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法 Download PDFInfo
- Publication number
- CN115566689A CN115566689A CN202211112654.4A CN202211112654A CN115566689A CN 115566689 A CN115566689 A CN 115566689A CN 202211112654 A CN202211112654 A CN 202211112654A CN 115566689 A CN115566689 A CN 115566689A
- Authority
- CN
- China
- Prior art keywords
- peak
- valley
- time
- load
- user
- 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
- 230000005611 electricity Effects 0.000 title claims abstract description 121
- 238000000034 method Methods 0.000 title claims abstract description 66
- 238000005457 optimization Methods 0.000 title claims abstract description 40
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 37
- 230000004044 response Effects 0.000 claims abstract description 21
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 238000012937 correction Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 16
- 239000013598 vector Substances 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 7
- 239000002245 particle Substances 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000000739 chaotic effect Effects 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 6
- 238000013459 approach Methods 0.000 claims description 5
- 239000000523 sample Substances 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 3
- 238000006529 Prey cleavage reaction Methods 0.000 claims description 2
- 238000007621 cluster analysis Methods 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000005315 distribution function Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 206010036590 Premature baby Diseases 0.000 claims 1
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims 1
- 239000000126 substance Substances 0.000 claims 1
- 238000011410 subtraction method Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 description 6
- 238000009472 formulation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000003631 expected effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
- H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
- H02J3/144—Demand-response operation of the power transmission or distribution network
-
- 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
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0206—Price or cost determination based on market factors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
- H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
-
- 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
- H02J2310/00—The network for supplying or distributing electric power characterised by its spatial reach or by the load
- H02J2310/50—The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads
- H02J2310/56—The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads characterised by the condition upon which the selective controlling is based
- H02J2310/58—The condition being electrical
- H02J2310/60—Limiting power consumption in the network or in one section of the network, e.g. load shedding or peak shaving
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Economics (AREA)
- Accounting & Taxation (AREA)
- Finance (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- Power Engineering (AREA)
- Computational Linguistics (AREA)
- Biomedical Technology (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Biophysics (AREA)
- Game Theory and Decision Science (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Human Resources & Organizations (AREA)
- Primary Health Care (AREA)
- Tourism & Hospitality (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，1)采用模糊半梯形隶属度函数，分别计算出各时点的偏大型和偏小型半梯度隶属度；2)以计算的各时点的偏大型和偏小型半梯度隶属度为统计指标，基于模糊聚类分析方法，得出基础时段划分方案；3)根据基础时段划分方案计算阈值指标，进一步确定偏大型和偏小型半梯形隶属函数的最优阈值，经过修正策略得出最终峰谷时段划分方案；4)建立用户负荷的需求响应模型、建立同时考虑削峰填谷效果和用户满意度的组合模型；5)利用改进天鹰优化算法找到组合模型下的最优峰谷电价。该方法可以使用户综合满意度更高，同时更好的达到削峰填谷的目的，实现供电公司与用户的双赢。
Description
技术领域
本发明属于电力市场领域，具体涉及一种峰谷时段划分和峰谷电价的优化方法。
背景技术
需求侧管理(DSM)是通过价格信号引导电力消费者采取合理的用电结构和方式的一种手段，在一些国家已经取得了一定成果。采用分时电价是DSM的重要途径之一，其中峰谷分时电价是分时电价的最主要组成部分。其最基本思想是体现电能在负荷高峰时作为短缺商品的价值，运用价格杠杆的作用引导终端用户的用电行为，提高电网安全性以及负荷率水平，引导用户根据自身生产方式的可调节性和利益改变用电方式，进而影响系统负荷。峰谷时段的划分和峰谷电价比价关系的确定是制定峰谷分时电价的基础。
根据峰谷时段划分的理论，目前已开展的划分方法主要有以下2种：①基于供电成本变化分析。该方法结合机组运行状况建立发电成本-负荷函数，利用该函数在负荷点的突变特征来划分峰谷时段区间，该方法较为复杂，在实际中不易实施。②基于负荷曲线分布分析。该方法一般采用模糊半梯形隶属度函数法，根据负荷曲线上各点分别处于峰、谷时段的可能性来划分峰谷时段区间，该方法可操作性强，但是根据隶属度函数阈值难以合理界定各时段区间分界时刻点。
峰谷电价(Time of use Price，TOU Price)是一种基于价格的需求侧管理措施，合理的峰谷电价水平可以达到削峰填谷、实现供电公司和用户双赢的目的。但是，如果峰谷电价比太高将导致用户对电价响应过度，峰谷时段产生了较大的漂移，甚至产生峰谷倒置，调峰失败的同时电网经济利益受损。峰谷电价比太低又会使用户响应不足，无法达到峰谷电价制定的预期效果。所以，有效地测量和量化用户对峰谷分时电价的响应是十分必要的。现有的峰谷电价难以充分挖掘需求响应的效益，峰谷时段划分也不太合理，另一方面，制定峰谷分时电价时也没有充分考虑用户对该政策的满意度，寻找经济效益和用户满意度以及其他系统目标之间的均衡点。
发明内容
为解决上述问题，本发明提出一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，利用模糊隶属度函数法和模糊聚类组合方法，根据构造的阈值指标函数λtif和修正方案，解决模糊隶属度函数法难以界定峰、谷区间分解点的问题，同时采用改进的天鹰优化算法优化峰谷电价，建立的目标函数综合考虑了削峰填谷的效果和用户的综合满意度，通过与某地区原峰谷电价方案相比，本发明提出的方法峰谷时段划分更符合负荷曲线特征，优化后的峰谷定价有效减少了最大负荷和负荷峰谷差，改进的天鹰优化算法与传统粒子群算法优化结果进行对比，解决粒子群算法易陷入局部最优，收敛速度慢的缺点。
本发明采用以下技术方案：
一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，包括以下步骤：
步骤一，根据电力系统负荷数据，以1h为单位将一天分为24个时点，组成时间点集合T＝{t1，t2，…，t24}，各时间点对应的负荷值集合为Q＝{q1，q2，…，q24}，利用模糊半梯形隶属函数法，分别采用偏大型半梯度隶属度函数和偏小型半梯度隶属度函数，初步确定ti时刻的峰、谷隶属度xi1和xi2；
步骤二，采用模糊聚类分析方法，以各时刻的峰、谷隶属度xi1和xi2做为统计指标，经过标定以及聚类等步骤，得到峰、平、谷时段的聚类集合Tf、Tp和Tg；
步骤三，运用集合分类思想，基于模糊聚类分析方法得出的基础时段划分结果，计算构造的阈值指标函数λtif，进一步确定偏大型和偏小型半梯形隶属函数的最优阈值λ1和λ2；
步骤四，根据偏大型和偏小型半梯形隶属函数的最优阈值λ1和λ2，将其与ti时刻的峰、谷隶属度xi1和xi2比较得出峰谷时段划分结果，针对隶属度函数阈值难以合理界定各时段区间分界时刻点的问题，参考基础方案和已实施的时段划分方案的特点，经过修正策略得出最终峰谷时段划分方案；
步骤五，根据最终峰谷时段划分方案，建立用户弹性响应矩阵，并建立用户负荷的需求响应模型；
步骤六，构建用户用电方式满意度和用户电费支出满意度模型，经过线性加权得出用户综合满意度模型，构建以系统谷值负荷最大、负荷峰谷差值最小、用户综合满意度最高为目标函数的峰谷分时电价及负荷优化模型；约束条件分别为电价调整前后总用电负荷需求量不变、峰、谷电价边界约束，用户和供电公司收益约束；利用改进天鹰优化算法获得优化后的分时电价和该分时电价下的需求响应负荷曲线。
进一步的，所述步骤一中，ti时刻的峰、谷隶属度xi1和xi2，由如下方法计算得到：
进一步的，在所述步骤二中，采用模糊聚类分析方法包括以下步骤：
⑴以各时点负荷为分类对象，以各时刻的峰、谷隶属度xi1和xi2做为统计指标，此时：xi＝(xi1，xi2)，式中：i＝1，2，…，24
得到特征指标矩阵X为：
⑵将特征指标矩阵X做标准化变换，基于绝对值减数法建立模糊相似矩阵R(r)24×24；其中R(xi1，xi2)＝rij，rij由下式计算：
式中：i，j＝1，2，…，24；m＝1，2；σ为适当选取的参数，使得rij∈[0,1]，
|xim-xjm|表示xi与xj之间的距离；
⑶对相似矩阵R(r)24×24求平方，即R→R2→…→R2k，直到第一次出现Rk×Rk＝Rk，Rk即为代求的传递闭包，记t(R)＝(tij)24×24；
⑷求传递闭包的截矩阵Rδ，Rδ＝(δγij)24×24，δγij由下式计算：
式中：i，j＝1，2，…，24；δ∈[0,1]，令δ从1逐渐减小，按Rδ进行动态聚类，聚类数取3，即可得到峰、平、谷时段的聚类聚合Tf，Tp和Tg。
进一步的，在所述步骤三中，运用集合分类思想，根据模糊聚类分析方法所得的峰、平、谷时段的聚类集合Tf，Tp和Tg，阈值指标函数λtif由如下方法计算得到：
Pfl、Ppu、Ppl和Pgu分别为Tf集合时点对应负荷中最小负荷值、Tp集合时点对应负荷中最大负荷值、Tp集合时点对应负荷中最小负荷值和Tg集合时点对应负荷中最大负荷值，Pi、Pj为Tf集合时点对应负荷中任意两个不同的负荷值，βi、βj为Tp集合时点对应负荷中任意两个不同的负荷值，αi、αj为Tg集合时点对应负荷中任意两个不同的负荷值；阈值指标函数λtif分子表征相邻集合中样本之间的最小距离，分母表征集合内样本间的最大差距，根据集合分类思想，使不同集合中的样本差异明显化，同一集合内的样本差异模糊化，因此，λtif越大说明峰、平、谷时段划分越好。
进一步的，在所述步骤四中，针对隶属度函数阈值难以界定各时段区间分界时刻点问题，作出以下修正方案：
①高峰、平时、低谷3个时段集合内的时间点数量应控制在6～10个；划分方案应符合典型日负荷曲线特征，例如典型的双峰性，上午和下午一般各出现一次用电高峰，低谷通常出现在夜间。具体可通过修正隶属度函数阈值，然后比较阈值指标函数λtif来合理规范峰、平、谷时段集合的大小；
②所划分的时段应具有可行性和实施的便利性，因此高峰、平时、低谷集合中的每个时段不应小于2h；如果出现某个时间点被孤立出来时，要根据模糊聚类分析方法得出的基础方案和孤立点及其相邻点对应负荷值的相对大小进行调整。
进一步的，在所述步骤五中，基于最终峰谷时段划分方案，确定用户弹性响应矩阵，并建立的用户负荷的需求响应模型如下：
用户弹性响应矩阵M如下：
式中，下标f、p、g分别代表峰、平、谷3个时段，对角线上的元素表征自弹性系数，非对角线上的元素表征交叉弹性系数；
用户负荷需求响应后的用电量为各时段原用电量与用电量变化的和，模型如下：
式中，ETOU＝[ef ep eg]T为实施峰谷分时电价后的各时段用电量；E0＝[e0f e0p e0g]T为原各时段用电量；p0f、p0p、p0g分别为实施分时电价前各时段的电价；Δpf、Δpp、Δpg分别为实施分时电价前后各时段的电价变化量；在此基础上将峰、平、谷各时段求的总的用电量的变化量，分别以各时段内每小时原用电量对应的比例系数分摊到各个小时，便可得出实施分时电价后各个小时用电量的变化，进而得出实施后的各时点负荷值。
进一步的，在所述步骤六中，建立的用户用电方式满意度和用户电费支出满意度模型如下：
用电方式满意度μ表示为：
其中，
表示实行峰谷分时电价后各时间点用电量变化的和，式中，STOU,t(Pf,Pp,Pg)为实行分时电价后的t时段的用电负荷，Pf、Pp、Pg分别为峰、平、谷时段的电价，St(Pt)表示未实行峰谷分时电价时时段的用电负荷，它是t时段电价Pt的函数。
μ体现了用户用电时间调整后的舒适度，μ∈[0,1]，即用户各时间点用电量未改变时用电方式满意度最高；
用电支出满意度
表示为：
其中，CTOU为实行峰谷分时电价后用户总的电费支出，C0为实施前用户总的电费支出；
用户综合满意度采用用电方式满意度和电费支出满意度的加权平均数，具体模型为：
其中，ε1+ε2＝1，ε1和ε2根据用户对用电方式和电费支出的重视程度灵活选取，可根据用户类型参考模糊描述的赋值方法选取；
根据步骤六中以系统谷值负荷最大、负荷峰谷差值最小、用户综合满意度最高为目标函数，构建的峰谷分时电价及负荷优化模型表示如下：
其中，ω1+ω2＝1；ω1和ω2表征用户对削峰填谷和综合满意度的权重，minQ和maxQ-minQ分别表示未实行峰谷分时电价前最小负荷和峰谷差值，minQTOU和maxQTOU-minQTOU分别表示实行峰谷分时电价后最小负荷和峰谷差值，二者做比值计算是为将其取值范围设定为1左右，以便与用户综合满意度Ψ在同一水平比较；
改进天鹰优化算法IAO包括以下步骤：
步骤1，算法初始化：
步骤1.1，初始化种群过程中引入混沌映射，嵌入粒子之间的位置信息，使得搜素范围更广；设置种群数量为N，解空间维度为Dim，最大迭代次数为Miter，设置两个计数器C1＝0，C2＝0；初始化N组位置向量，第一组X1＝[x11,x12,…,x1Dim]随机生成，剩下的N-1组由式：Xn+1＝θXn(1-Xn)计算得到，θ为变化因子，可根据θ的变化产生不同的混沌状态；初始化t＝0；
步骤1.2，利用混沌映射初始化过程中，需要满足条件cos(Xi,Xj)≤ε，i，j＝1，2，…，N，i≠j；其中，cos(Xi,Xj)表示位置向量之间的距离，该值越小表示两个粒子之间越远，计算方式如下：
ε∈[0,1]为阈值，可凭经验和统计实验给定。
步骤1.3，将初始化完成的N组位置向量在解空间[Xmin,Xmax]中进行映射；
步骤2，计算种群的适应度值fitness，并对这些位置向量进行排序，初步确定全局最优解；
步骤3，在每次迭代的过程中，天鹰的位置更新描述如下：
步骤3.1，扩大搜索
此时天鹰识别猎物区域，并通过垂直弯腰的高飞选择最佳的狩猎区域；为了增加种群的多样性，解决“早熟”问题，可在位置更新中加入余弦变化因子，数学模型如下：
其中，X1(t+1)为第一种搜索方法生成的t+1次迭代的解；Xbest(t)代表第t次迭代之前获得的最佳解，这反映了猎物的近似位置；rand是介于0和1之间的随机值；XM(t)表示在第t次迭代时当前解的平均值，计算方式如下：
步骤3.2，缩小搜索
在第二种方法中，当天鹰从高空发现猎物区域时，天鹰在目标猎物上方盘旋，准备着陆陆地，然后攻击，这种方法称为短滑翔攻击的等高线飞行；在这里，AO狭窄地探索目标猎物的选定区域，为攻击做准备；该行为的数学模式如下式所示：
X2(t+1)＝Xbest(t)×Levy(Dim)+XR(t)+(y-x)×rand
式中，X2(t+1)是由第二个搜索方法生成的第t+1次迭代的解；Levy(Dim)是Levy飞行分布函数，XR(t)是在第t次迭代时在[1,N]范围内获得的随机解；
式中，s＝0.01，u和v为介于0和1之间的随机数，β＝1.5，σ计算式如下：
x和y表示搜索中的螺旋形状，其中，x＝r×cos(θ)；y＝r×sin(θ)；r＝r1+U×D1；
r1取1到20之间的值，用于固定搜索周期数，U＝0.00565，D1是从1到搜索空间维数Dim的整数，ω＝0.005；
步骤3.3，扩展开发
在第三种方法中，当天鹰准确地指定了猎物区域，并且准备好着陆和攻击时，天鹰垂直下降并进行初步攻击，以试探猎物反应；这种方法称为低空飞行下降攻击；在这里，AO利用目标的选定区域接近猎物并进行攻击，该行为的数学模式如下式所示：
X3(t+1)＝(Xbest(t)-XM(t))×α-rand+((UB-LB)×rand+LB)×δ
式中，X3(t+1)是由第三种搜索方法生成的第t+1次迭代的解；α和δ是固定为较小值(0.1)的开发调整参数，UB和LB是位置向量上下界；
步骤3.4，缩小开发
在第四种方法中，当天鹰接近猎物时，天鹰根据其随机运动在陆地上攻击猎物，这种方法称为“行走并抓住猎物”，最后，AO在最后一个位置攻击猎物，该行为的数学模式如下式所示：
X4(t+1)＝QF(t)×Xbest(t)-(G1×X(t)×rand)-G2×Levy(D)+Rand×G1
式中，X4(t+1)是由第四种搜索方法生成的第t+1次迭代的解；QF(t)表示在第t次迭代时用于平衡搜索策略的质量函数，G1表示在搜索猎物期间用于跟踪猎物的的各种运动，G2表示天鹰跟踪猎物的飞行速率；Levy(D)是莱维飞行函数，X(t)是第t次迭代时的位置。
质量函数计算式为：
G1＝2×rand-1；飞行速率由线性递减修改为：
步骤4，计算出更新后的适应度值，并且更新群体极值Xbest，如果此次迭代Xbest没有进行更新，则令C1＝C1+1；
步骤5，评估算法是否存在早熟现象并修正；
步骤5.1，若C1≥2，则C2＝C2+1，同时令C1＝0；
步骤5.2，若C2≥γ，γ∈[5,10]，则令群体极值Xbest沿其负梯度方向进行更新，计算方式如下：
式中，
为Xbest最近一次更新之前的值；
步骤6，判断是否满足终止条件，若满足则输出最优解，结束程序；否则，重复以上天鹰优化算法流程继续进行寻优迭代处理。
本发明的有益效果：
通过与某地区原峰谷电价方案相比，基于本发明提出的方法峰谷时段划分更符合负荷曲线特征，利用改进的天鹰算法优化后的峰谷定价有效减少了最大负荷和负荷峰谷差，同时用户综合满意度较高。
附图说明
图1为本发明IAO算法优化负荷峰谷时段划分以及峰谷电价方法流程图；
图2为负荷数据与优化后阶梯电价图；
图3为PSO算法优化迭代曲线图；
图4为IAO算法优化迭代曲线图；
图5原负荷曲线与原峰谷电价下用户响应后负荷曲线比较图；
图6原负荷曲线与优化后峰谷电价下用户响应后负荷曲线比较图。
具体实施方式
图1为本发明IAO算法优化负荷峰谷时段划分以及峰谷电价方法流程图。下面以一个具体实例说明本方法的实施流程，该实例以某用户实际负荷数据验证方法有效性。未实行峰谷电价前系统电价为0.49元/(kW·h)，实行峰谷电价后的峰平谷电价分别为0.658元/(kW·h)，0.49元/(kW·h)，0.358元/(kW·h),原峰谷时段划分结果为：24：00～7：00为谷段，8：00-11：00以及18：00～24：00为峰段，7：00～8：00以及11：00～18：00为平段。
首先，采用模糊聚类分析方法，以各时刻的峰、谷隶属度xi1和xi2做为统计指标，经过标定以及聚类等步骤，得到峰、平、谷时段的聚类集合Tf、Tp和Tg。然后运用集合分类思想，基于模糊聚类分析方法得出的基础时段划分结果，计算构造的阈值指标函数λtif，进一步确定偏大型和偏小型半梯形隶属函数的最优阈值λ1和λ2。最后经过修正策略得出最终峰谷时段划分方案如表1所示。
表1时段划分方案
得到最终峰谷时段的划分结果之后，利用电量电价弹性矩阵和建立的用户需求响应模型，以系统谷值负荷最大、负荷峰谷差值最小、用户综合满意度最高为目标函数，利用改进的天鹰优化算法获得优化后的峰谷电价结果如表2所示：
表2峰谷电价制定结果
图2为负荷数据与优化后阶梯电价图。
根据峰谷电价优化结果，进一步计算对模型的合理性及优势进行分析，下表3和表4给出了实施优化后的峰谷电价与原峰谷电价的指标对比数据。
表3峰谷电价优化前后满意度结果对比
表4峰谷电价优化前后削峰填谷效果对比
本实施算例中算法参数取值如下：
PSO算法：粒子群群体大小为N＝100；维数Dim＝2，最大迭代次数Miter＝200。
IAO算法：天鹰群体大小为N＝100，维数Dim＝2，最大迭代次数Miter＝200，
该实例用户为工业用户，由于其为三班制生产，工艺流程调整较为容易，且产品产值电耗较高，因此取用户综合满意度为：
同时由于该用户在拉闸限电表中位置较靠前，故取ω1＝0.7，ω2＝0.3。
采用本发明提出的改进天鹰优化算法(IAO)优化负荷峰谷时段划分以及峰谷电价方法，并与粒子群算法(PSO)优化结果进行对比，通过图3-图6的结果可知，本发明相比传统PSO算法优化，收敛速度更快，寻优精度更高，通过表3和表4结果可知，优化后的分时电价比原峰时电价削峰填谷效果更好，用户满意度更高。
最后说明的是，以上仅对本发明具体实施例进行详细描述说明。但本发明并不限制于以上描述具体实施例。本领域的技术人员对本发明进行的等同修改和替代也都在本发明的范畴之中。因此，在不脱离本发明的精神和范围下所作的均等变换和修改，都涵盖在本发明范围内。
Claims (7)
1.一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，包括以下步骤：
步骤一，根据电力系统负荷数据，以1h为单位将一天分为24个时点，组成时间点集合T＝{t1，t2，...，t24}，各时间点对应的负荷值集合为Q＝{q1，q2，...，q24}，利用模糊半梯形隶属函数法，分别采用偏大型半梯度隶属度函数和偏小型半梯度隶属度函数，初步确定ti时刻的峰、谷隶属度xi1和xi2；
步骤二，采用模糊聚类分析方法，以各时刻的峰、谷隶属度xi1和xi2做为统计指标，经过标定以及聚类步骤，得到峰、平、谷时段的聚类集合Tf、Tp和Tg；
步骤三，运用集合分类思想，基于模糊聚类分析方法得出的基础时段划分结果，计算构造的阈值指标函数λtif，进一步确定偏大型和偏小型半梯形隶属函数的最优阈值λ1和λ2；
步骤四，根据偏大型和偏小型半梯形隶属函数的最优阈值λ1和λ2，将其与ti时刻的峰、谷隶属度xi1和xi2比较得出峰谷时段划分结果，针对隶属度函数阈值难以合理界定各时段区间分界时刻点的问题，参考基础方案和已实施的时段划分方案的特点，经过修正策略得出最终峰谷时段划分方案；
步骤五，根据最终峰谷时段划分方案，建立用户弹性响应矩阵，并建立用户负荷的需求响应模型；
步骤六，构建用户用电方式满意度和用户电费支出满意度模型，经过线性加权得出用户综合满意度模型，构建以系统谷值负荷最大、负荷峰谷差值最小、用户综合满意度最高为目标函数的峰谷分时电价及负荷优化模型；约束条件分别为电价调整前后总用电负荷需求量不变、峰、谷电价边界约束，用户和供电公司收益约束；利用改进天鹰优化算法获得优化后的分时电价和该分时电价下的需求响应负荷曲线。
2.根据权利要求1所述的一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，所述步骤一中，ti时刻的峰、谷隶属度xi1和xi2，由如下方法计算得到：
3.根据权利要求1所述的一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，在所述步骤二中，采用模糊聚类分析方法包括以下步骤：
(1)以各时点负荷为分类对象，以各时刻的峰、谷隶属度xi1和xi2做为统计指标，此时：xi＝(xi1，xi2)，式中：i＝1，2，…，24
得到特征指标矩阵X为：
(2)将特征指标矩阵X做标准化变换，基于绝对值减数法建立模糊相似矩阵R(r)24×24；其中R(xi1，xi2)＝rij，rij由下式计算：
式中：i，j＝1，2，…，24；m＝1，2；σ为适当选取的参数，使得rij∈[0，1]，
|xim-xjm|表示xi与xj之间的距离；
(3)对相似矩阵R(r)24×24求平方，即R→R2→...→R2k，直到第一次出现Rk×Rk＝Rk，Rk即为代求的传递闭包，记t(R)＝(tij)24×24；
(4)求传递闭包的截矩阵Rδ，Rδ＝(δγij)24×24，δγij由下式计算：
式中：i，j＝1，2，…，24；δ∈[0，1]，令δ从1逐渐减小，按Rδ进行动态聚类，聚类数取3，即可得到峰、平、谷时段的聚类聚合Tf，Tp和Tg。
4.根据权利要求1所述的一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，在所述步骤三中，运用集合分类思想，根据模糊聚类分析方法所得的峰、平、谷时段的聚类集合Tf，Tp和Tg，阈值指标函数λtif由如下方法计算得到：
Pfl、Ppu、Ppl和Pgu分别为Tf集合时点对应负荷中最小负荷值、Tp集合时点对应负荷中最大负荷值、Tp集合时点对应负荷中最小负荷值和Tg集合时点对应负荷中最大负荷值，Pi、Pj为Tf集合时点对应负荷中任意两个不同的负荷值，βi、βj为Tp集合时点对应负荷中任意两个不同的负荷值，αi、αj为Tg集合时点对应负荷中任意两个不同的负荷值；阈值指标函数λtif分子表征相邻集合中样本之间的最小距离，分母表征集合内样本间的最大差距，根据集合分类思想，使不同集合中的样本差异明显化，同一集合内的样本差异模糊化，因此，λtif越大说明峰、平、谷时段划分越好。
5.根据权利要求1所述的一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，在所述步骤四中，针对隶属度函数阈值难以界定各时段区间分界时刻点问题，作出以下修正方案：
①高峰、平时、低谷3个时段集合内的时间点数量应控制在6～10个；划分方案应符合典型日负荷曲线特征；具体可通过修正隶属度函数阈值，然后比较阈值指标函数λtif来合理规范峰、平、谷时段集合的大小；
②所划分的时段应具有可行性和实施的便利性，因此高峰、平时、低谷集合中的每个时段不应小于2h；如果出现某个时间点被孤立出来时，要根据模糊聚类分析方法得出的基础方案和孤立点及其相邻点对应负荷值的相对大小进行调整。
6.根据权利要求1所述的一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，在所述步骤五中，基于最终峰谷时段划分方案，确定用户弹性响应矩阵，并建立的用户负荷的需求响应模型如下：
用户弹性响应矩阵M如下：
式中，下标f、p、g分别代表峰、平、谷3个时段，对角线上的元素表征自弹性系数，非对角线上的元素表征交叉弹性系数；
用户负荷需求响应后的用电量为各时段原用电量与用电量变化的和，模型如下：
式中，ETOU＝[ef ep eg]T为实施峰谷分时电价后的各时段用电量；E0＝[e0f e0p e0g]T为原各时段用电量；p0f、p0p、p0g分别为实施分时电价前各时段的电价；Δpf、Δpp、Δpg分别为实施分时电价前后各时段的电价变化量；在此基础上将峰、平、谷各时段求的总的用电量的变化量，分别以各时段内每小时原用电量对应的比例系数分摊到各个小时，得出实施分时电价后各个小时用电量的变化，进而得出实施后的各时点负荷值。
7.根据权利要求1所述的一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法，其特征在于，在所述步骤六中，建立的用户用电方式满意度和用户电费支出满意度模型如下：
用电方式满意度μ表示为：
其中，
表示实行峰谷分时电价后各时间点用电量变化的和，式中，STOU，t(Pf，Pp，Pg)为实行分时电价后的t时段的用电负荷，Pf、Pp、Pg分别为峰、平、谷时段的电价，St(Pt)表示未实行峰谷分时电价时时段的用电负荷，它是t时段电价Pt的函数；
μ体现了用户用电时间调整后的舒适度，μ∈[0，1]，即用户各时间点用电量未改变时用电方式满意度最高；
用电支出满意度
表示为：
其中，CTOU为实行峰谷分时电价后用户总的电费支出，C0为实施前用户总的电费支出；
用户综合满意度采用用电方式满意度和电费支出满意度的加权平均数，具体模型为：
其中，ε1+ε2＝1，ε1和ε2根据用户对用电方式和电费支出的重视程度灵活选取，可根据用户类型参考模糊描述的赋值方法选取；
根据步骤六中以系统谷值负荷最大、负荷峰谷差值最小、用户综合满意度最高为目标函数，构建的峰谷分时电价及负荷优化模型表示如下：
其中，ω1+ω2＝1；ω1和ω2表征用户对削峰填谷和综合满意度的权重，min Q和max Q-min Q分别表示未实行峰谷分时电价前最小负荷和峰谷差值，min QTOU和max QTOU-min QTOU分别表示实行峰谷分时电价后最小负荷和峰谷差值，二者做比值计算是为将其取值范围设定为1左右，以便与用户综合满意度Ψ在同一水平比较；
改进天鹰优化算法IAO包括以下步骤：
步骤1，算法初始化：
步骤1.1，初始化种群过程中引入混沌映射，嵌入粒子之间的位置信息，使得搜素范围更广；设置种群数量为N，解空间维度为Dim，最大迭代次数为Miter，设置两个计数器C1＝0，C2＝0；初始化N组位置向量，第一组X1＝[x11，x12，...，x1Dim]随机生成，剩下的N-1组由式：Xn+1＝θXn(1-Xn)计算得到，θ为变化因子，可根据θ的变化产生不同的混沌状态；初始化t＝0；
步骤1.2，利用混沌映射初始化过程中，需要满足条件cos(Xi，Xj)≤ε，i，j＝1，2，...，N，i≠j；其中，cos(Xi，Xj)表示位置向量之间的距离，该值越小表示两个粒子之间越远，计算方式如下：
ε∈[0，1]为阈值；
步骤1.3，将初始化完成的N组位置向量在解空间[Xmin，Xmax]中进行映射；
步骤2，计算种群的适应度值fitness，并对这些位置向量进行排序，初步确定全局最优解；
步骤3，在每次迭代的过程中，天鹰的位置更新描述如下：
步骤3.1，扩大搜索
此时天鹰识别猎物区域，并通过垂直弯腰的高飞选择最佳的狩猎区域；为了增加种群的多样性，解决“早熟”问题，可在位置更新中加入余弦变化因子，数学模型如下：
其中，X1(t+1)为第一种搜索方法生成的t+1次迭代的解；Xbest(t)代表第t次迭代之前获得的最佳解，这反映了猎物的近似位置；rand是介于0和1之间的随机值；XM(t)表示在第t次迭代时当前解的平均值，计算方式如下：
步骤3.2，缩小搜索
在第二种方法中，当天鹰从高空发现猎物区域时，天鹰在目标猎物上方盘旋，准备着陆陆地，然后攻击，这种方法称为短滑翔攻击的等高线飞行；在这里，AO狭窄地探索目标猎物的选定区域，为攻击做准备；该行为的数学模式如下式所示：
X2(t+1)＝Xbest(t)×Levy(Dim)+XR(t)+(y-x)×rand
式中，X2(t+1)是由第二个搜索方法生成的第t+1次迭代的解；Levy(Dim)是Levy飞行分布函数，XR(t)是在第t次迭代时在[1，N]范围内获得的随机解；
式中，s＝0.01，u和v为介于0和1之间的随机数，β＝1.5，σ计算式如下：
x和y表示搜索中的螺旋形状，其中，x＝r×cos(θ)；y＝r×sin(θ)；r＝r1+U×D1；
r1取1到20之间的值，用于固定搜索周期数，U＝0.00565，D1是从1到搜索空间维数Dim的整数，ω＝0.005；
步骤3.3，扩展开发
在第三种方法中，当天鹰准确地指定了猎物区域，并且准备好着陆和攻击时，天鹰垂直下降并进行初步攻击，以试探猎物反应；这种方法称为低空飞行下降攻击；在这里，AO利用目标的选定区域接近猎物并进行攻击，该行为的数学模式如下式所示：
X3(t+1)＝(Xbest(t)-XM(t))×α-rand+((UB-LB)×rand+LB)×δ
式中，X3(t+1)是由第三种搜索方法生成的第t+1次迭代的解；α和δ是固定为较小值(0.1)的开发调整参数，UB和LB是位置向量上下界；
步骤3.4，缩小开发
在第四种方法中，当天鹰接近猎物时，天鹰根据其随机运动在陆地上攻击猎物，这种方法称为“行走并抓住猎物”，最后，AO在最后一个位置攻击猎物，该行为的数学模式如下式所示：
X4(t+1)＝QF(t)×Xbest(t)-(G1×X(t)×rand)-G2×Levy(D)+Rand×G1
式中，X4(t+1)是由第四种搜索方法生成的第t+1次迭代的解；QF(t)表示在第t次迭代时用于平衡搜索策略的质量函数，G1表示在搜索猎物期间用于跟踪猎物的的各种运动，G2表示天鹰跟踪猎物的飞行速率；Levy(D)是莱维飞行函数，X(t)是第t次迭代时的位置；
质量函数计算式为：
G1＝2×rand-1；飞行速率由线性递减修改为：
步骤4，计算出更新后的适应度值，并且更新群体极值Xbest，如果此次迭代Xbest没有进行更新，则令C1＝C1+1；
步骤5，评估算法是否存在早熟现象并修正；
步骤5.1，若C1≥2，则C2＝C2+1，同时令C1＝0；
步骤5.2，若C2≥γ，γ∈[5，10]，则令群体极值Xbest沿其负梯度方向进行更新，计算方式如下：
式中，
为Xbest最近一次更新之前的值；
步骤6，判断是否满足终止条件，若满足则输出最优解，结束程序；否则，重复以上天鹰优化算法流程继续进行寻优迭代处理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211112654.4A CN115566689A (zh) | 2022-09-14 | 2022-09-14 | 一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211112654.4A CN115566689A (zh) | 2022-09-14 | 2022-09-14 | 一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115566689A true CN115566689A (zh) | 2023-01-03 |
Family
ID=84741745
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211112654.4A Pending CN115566689A (zh) | 2022-09-14 | 2022-09-14 | 一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115566689A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116187099A (zh) * | 2023-04-24 | 2023-05-30 | 山东理工大学 | 一种基于双层迭代的用户侧储能配置方法 |
| CN117060425A (zh) * | 2023-10-12 | 2023-11-14 | 国网浙江省电力有限公司宁波供电公司 | 一种基于强化学习的配网峰谷差自适应控制方法及系统 |
-
2022
- 2022-09-14 CN CN202211112654.4A patent/CN115566689A/zh active Pending
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116187099A (zh) * | 2023-04-24 | 2023-05-30 | 山东理工大学 | 一种基于双层迭代的用户侧储能配置方法 |
| CN117060425A (zh) * | 2023-10-12 | 2023-11-14 | 国网浙江省电力有限公司宁波供电公司 | 一种基于强化学习的配网峰谷差自适应控制方法及系统 |
| CN117060425B (zh) * | 2023-10-12 | 2024-04-09 | 国网浙江省电力有限公司宁波供电公司 | 一种基于强化学习的配网峰谷差自适应控制方法及系统 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108416366B (zh) | 一种基于气象指数的加权ls-svm的电力系统短期负荷预测方法 | |
| Liu et al. | Random forest solar power forecast based on classification optimization | |
| Li et al. | Short-term wind power forecasting based on support vector machine with improved dragonfly algorithm | |
| Mahmoud et al. | An advanced approach for optimal wind power generation prediction intervals by using self-adaptive evolutionary extreme learning machine | |
| CN106600059B (zh) | 基于改进rbf神经网络的智能电网短期负荷预测方法 | |
| CN115566689A (zh) | 一种改进天鹰优化算法优化负荷峰谷时段划分以及峰谷电价方法 | |
| CN112288164B (zh) | 一种计及空间相关性和修正数值天气预报的风功率组合预测方法 | |
| Huang et al. | A clustering based grouping method of nearly zero energy buildings for performance improvements | |
| CN109214095A (zh) | 电动汽车充放电多目标优化调度方法 | |
| CN113300380B (zh) | 一种基于负荷曲线分段的配电网无功优化补偿方法 | |
| CN112186761B (zh) | 一种基于概率分布的风电功率场景生成方法及系统 | |
| CN110380444A (zh) | 一种基于变结构Copula的多场景下分散式风电有序接入电网的容量规划方法 | |
| CN115986728A (zh) | 一种计及不确定性因素的电网态势预测方法与终端 | |
| Wang et al. | Big data analytics for price forecasting in smart grids | |
| CN108694475B (zh) | 基于混合模型的短时间尺度光伏电池发电量预测方法 | |
| CN113688960A (zh) | 基于灰狼优化ghfcm的居民电力数据聚类方法及装置 | |
| CN113659631A (zh) | 考虑时变特性的风光电站群出力描述方法 | |
| CN115829418B (zh) | 适用于负荷管理的电力用户负荷特性画像构建方法和系统 | |
| CN112381315A (zh) | 一种基于pso优化的ls-svm智能台区负荷预测方法及系统 | |
| CN115622056B (zh) | 一种基于线性加权和选解法的储能优化配置方法及系统 | |
| CN110163437A (zh) | 基于DPK-means的日前光伏发电功率预测方法 | |
| CN112508280B (zh) | 一种综合能源系统实时性能的评价和优化方法及系统 | |
| CN116468138A (zh) | 空调负荷预测方法、系统、电子设备及计算机存储介质 | |
| CN115689614A (zh) | 一种新型供需关系下地区峰谷分时电价优化方法 | |
| CN109767037A (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 |