CN113312779B - 一种低碳型柔性配电网高满意度动态综合规划方法 - Google Patents
一种低碳型柔性配电网高满意度动态综合规划方法 Download PDFInfo
- Publication number
- CN113312779B CN113312779B CN202110629893.6A CN202110629893A CN113312779B CN 113312779 B CN113312779 B CN 113312779B CN 202110629893 A CN202110629893 A CN 202110629893A CN 113312779 B CN113312779 B CN 113312779B
- Authority
- CN
- China
- Prior art keywords
- power
- stage
- distribution network
- decision
- formula
- 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.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 76
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000004519 manufacturing process Methods 0.000 claims abstract description 9
- 238000012423 maintenance Methods 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 5
- 230000005611 electricity Effects 0.000 claims abstract 2
- 238000010606 normalization Methods 0.000 claims abstract 2
- 230000006870 function Effects 0.000 claims description 23
- 238000005457 optimization Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 15
- 230000008901 benefit Effects 0.000 claims description 11
- 238000007599 discharging Methods 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 230000002776 aggregation Effects 0.000 claims description 5
- 238000004220 aggregation Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 5
- 238000011835 investigation Methods 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 3
- 238000003064 k means clustering Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 229910002056 binary alloy Inorganic materials 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 238000011084 recovery Methods 0.000 claims description 2
- 238000003786 synthesis reaction Methods 0.000 claims description 2
- 238000013459 approach Methods 0.000 claims 1
- 230000008859 change Effects 0.000 claims 1
- 229910052731 fluorine Inorganic materials 0.000 claims 1
- 125000001153 fluoro group Chemical group F* 0.000 claims 1
- 229910052739 hydrogen Inorganic materials 0.000 claims 1
- 239000001257 hydrogen Substances 0.000 claims 1
- 125000004435 hydrogen atom Chemical class [H]* 0.000 claims 1
- 238000011160 research Methods 0.000 claims 1
- 238000007621 cluster analysis Methods 0.000 abstract 1
- 238000011161 development Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000007726 management method Methods 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提出一种低碳型柔性配电网高满意度动态综合规划方法,该方法将多群差分进化优化的多层分数阶泰勒动态网络规划应用到低碳型柔性配电网规划中。低碳型柔性配电网规划考虑设备扩容投资成本、设备运行维护成本、设备能源生产成本、设备能源减载成本和用户需求满意程度。该方法对采集后的数据进行聚类分析和归一化处理,对电力负荷以及居民用电满意度进行精准预测,依据预测数据对低碳型柔性配电网进行优化,确定最优低碳型柔性配电网规划方案,保证低碳型柔性配电网绿色高效安全经济运行。
Description
技术领域
本发明属于电力系统配电网的规划领域,涉及综合规划方法,适用于电力系统与智能电网中配电网的规划。
背景技术
目前,传统的配电系统扩展规划模型已经被用来确定新建电力设备的类型和容量,完成选址选线工作。随着分布式可再生能源渗透率的不断提高,配电网实时运行状态的不确定性逐渐增加,配电网各个主体之间仅仅拥有区域内的有限信息,缺乏必要的信息交流与有效的协调机制,易造成分布式电源及配电线路投资不足、重复投资、用电率较低问题。
另外,在配电系统扩展规划模型引入不确定因素时,现有的方法使用鲁棒性优化方法和随机规划方法来模拟不确定因素。鲁棒性优化方法能够在最坏的负载增长情况下实现配电系统最佳的性能。随机规划方法通过代表性的场景样本或者不确定性变量的概率分布来确定配电系统最佳性能的规划方案。但是上述方法未考虑长期发展过程中的不确定性因素,因此一旦使用上述方法确定规划方案,就无法再对规划方案进行进一步地优化调整。
因此,考虑投资成本和用户满意度,提出一种动态综合规划方法,解决分布式电源及配电线路投资不足、重复投资、用电率较低问题,解决无法动态规划的问题。
发明内容
本发明提出一种低碳型柔性配电网高满意度动态综合规划方法,该方法将多群差分进化优化的多层分数阶泰勒动态网络规划应用到低碳型柔性配电网规划中,完成对电力负荷以及居民用电满意度的精准预测,依据预测数据对低碳型柔性配电网进行优化,确定最优低碳型柔性配电网规划方案,保证低碳型柔性配电网绿色高效安全经济运行;低碳型柔性配电网高满意度动态综合规划方法在使用过程中的步骤为:
步骤(1):构建低碳型柔性配电网优化扩展模型,旨在保证负荷增长状态下的预期总成本达到最小化,同时数字化的用户满意程度达到最大化,保证低碳型柔性配电网规划策略处在最优解;总规划方案的目标函数f表示为:
式中,w为负荷预测误差指标序号,λt,w为t阶段第w个负荷预测误差指标水平,Ωt为阶段集,α为年利率,y表示规划年序号;为总规划方案的数学期望,Ny为规划阶段的总年数; A、B、C、D和F分别为设备扩容投资成本、设备运行维护成本、设备能源生产成本、设备能源减载成本和用户需求满意程度标识;和分别为t阶段的设备扩容投资成本、设备运行维护成本、设备能源生产成本、设备能源减载成本和用户需求满意程度;其中,c表示规划方案序号,t表示当前处于第t个阶段;i表示当前处在第i个设备节点上;m 表示当前处在第m条馈线支路上,k表示已经安装的设备类型序号;
设备扩容投资成本为:
式中,Lm表示第m条支路的建设总长度;TS、TR、FL、PV、WT和EV分别表示变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩标识;ΩTS、ΩTR、ΩFL、ΩPV、ΩWT和ΩEV分别表示变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的候选节点集;KTR、KFL、KPV、KWT和KEV分别表示变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的型号类型集;和分别表示在t阶段下第i个设备节点上已经安装的第k个变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的投资变量;和分别表示在已经安装的设备类型 k中,变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的投资成本系数;τTS、τTR、τFL、τPV、τWT和τEV分别表示变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的设备回收率;
设备运行维护成本为:
设备能源生产成本为:
式中,υc为方案c的持续时间;CH和DIS分别表示电动汽车充电桩的充电过程和放电过程标识;和分别为在方案c中变电站的能源购买成本、光伏电站和风力涡轮机的生产成本;和分别为方案c中电动汽车充电桩的充电成本和放电成本;和分别为方案c在t阶段的变电站、光伏电站和风力涡轮机注入的有功功率;和分别为方案c在t阶段的电动汽车充电桩充电时的有功功率和放电时的有功功率;
设备能源减载成本为:
式中,FH表示负荷标识,ΩFH表示负荷节点集;和分别为方案c在t阶段的变电站和负荷的视在功率削减量,和分别为方案c在t阶段的光伏电站和风力涡轮机的有功功率削减量;和分别为方案c中变电站、负荷、光伏电站和风力涡轮机的削减成本;
用户需求满意程度为:
式中,和分别为方案c中变电站、变压器、光伏电站、风力涡轮机和电动汽车充电桩拟定选址的用户需求满意程度系数;和分别表示第c个规划方案的在t阶段下第i个设备节点上变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的投资变量;
低碳型柔性配电网规划是对现有变电站、变压器以及输电线路进行改造或者新建,对太阳能光伏板、风力发电机和电动汽车充电桩进行选址,对各类设备的约束为:
为了保证新建输电线路不在两个及两个以上的方向同时规划,约束为:
为了保证低碳型柔性配电网有功功率和无功功率的平衡,对各个节点的功率限制为:
式中,和分别为以节点i为线路首端和以节点i为线路末端的馈线集,为方案c 在t阶段流经馈线的有功功率,和分别为方案c在t阶段的电动汽车充电桩充电时和放电时的有功功率,和分别为方案c在t阶段的节点预期和实际需求的有功功率;μc为方案c中的需求水平因子,ψc为方案c中的负载功率因数;为方案c在t阶段流经馈线的无功功率;和分别为方案c在t阶段的变电站、光伏电站和风力涡轮机注入的无功功率,和分别为方案c在t阶段的节点预期和实际消耗的无功功率;
馈线在决策后与决策前的视在功率关系为:
馈线的视在功率与有功功率和无功功率的关系为:
对馈线的功率限制为:
变电站在决策后与决策前的视在功率关系为:
变电站的视在功率与有功功率和无功功率的关系为:
对变电站的功率限制为:
光伏电站在决策后与决策前的有功功率关系为:
对光伏电站的功率限制为:
风力涡轮机在决策后与决策前的有功功率关系为:
对风力涡轮机的功率限制为:
电动汽车充电桩在决策后与决策前的有功功率关系为:
对电动汽车充电桩的最大充电功率和放电功率限制为:
为了保证低碳型柔性配电网在运行过程中网络枢纽点的电压水平以及网络中各个节点的电压水平满足技术要求,对电压的限制为:35kV及以上供电电压正、负偏差的绝对值之和不超过系统额定电压的10%;10kV及以下三相供电电压允许偏差为系统额定电压的±7%;220V 单相供电电压允许偏差为系统额定电压的-10%至+7%;
步骤(2):对目标区域的低碳型柔性配电网进行初步可行性研究,将计划扩建的电源场址制作成调查问卷,以入户访问、召开座谈会以及线上调研的方式,组织开展关于扩建电源选址满意程度的调查工作,统计各电源选址的满意度情况;与此同时,由电力用户智能信息采集终端对电力用户的用电情况,目标区域的用电经济效益以及电源地理分布数据进行实时采集;
步骤(3):在数字化信息平台上传步骤(2)中的采集数据样本,将数据样本进行k-means 聚类技术处理;其中,拟合数据和原始数据的误差平方和定义为:
x具体为[日负荷,季度负荷,年负荷,城乡居民用电负荷,工业用电负荷,需求电源点]T;
各个簇类的集聚程度定义为:
式中,RE表示簇类的集聚程度标识,RRE为簇类聚集程度系数,ΩC为总样本集,n为总样本数,假设某样本点xl已经被聚类到簇AA当中,r(xl)表示样本点xl与簇AA内的其他样本点的平均距离,表示样本点xl与除了簇AA以外的其他簇样本点的平均距离;当簇类的集聚程度最大时,聚类效果最优;当K值小于最佳聚类数值时,增大K值,进而加大各个簇类的集聚程度,使得误差平方和急剧下降;当增大K值至K值取到最佳聚类数时,各个簇类的集聚程度变化不明显,误差平方和也趋近于0,此时,数据样本达到了最佳的聚合效果;在确定最佳聚类数值K之后,获得原始数据特征集;
步骤(4):定义预处理数据的类型和特征;遍历原始数据特征集中的每一个数据,记录数据中的最小值ymin和最大值ymax;将ymin和ymax分别映射为数值0和数值1,对原始数据特征集中的其余数值进行归一化处理为:
式中,ynew为归一化处理后的数值,y为原始数据特征集的数值,ymin和ymax原始数据特征集的数据最小值和数据最大值;
步骤(5):将处理过后的数值区间均匀地离散化,以获取每个特征与类别的联合信息熵以及互信息;其中,联合信息熵定义为:
式中,L(b,d)为第b个特征与类别d的联合信息熵;u为区间数,v为类别数,Nb为特征b的区间数量,Nd为总类别数,MLSH为总样本数,MLSH,u,v表示第b个特征落在第u个区间并且类别为v的样本数;
互信息定义为:
式中,H(b;d)为第b个特征与类别d的互信息,MLSH,u为第b个特征落在第u个区间的样本数,MLSH,v为属于第v个类别的样本数;
最大相关性指标定义为:
式中,G(S,d)为最优特征集S与类别d的最大相关性指标;NS最优特征集S中特征的数量;
根据最大相关最小冗余准则对数据进行分析,筛选出冗余性最低并且相关性最高的优质数据特征集;最大相关最小冗余准则定义为:
步骤(6):将筛选后的优质数据特征集输入至步骤(1)建立的低碳型柔性配电网优化扩展模型中,使用多群差分进化优化的多层分数阶泰勒动态网络规划对该模型进行求解;其中,目标函数f(x)中的NTaylor阶泰勒展开为:
式中,f(x)是关于x=a的NTaylor阶泰勒多项式,c介于x与a之间;其中,NTaylor为整数;
利用分数阶微积分将泰勒级数展开成多项式,此时的函数f(x)能进行任意次的微分,公式为:
由公式(29)和公式(30)重新定义步骤(1)中模型的目标函数:
步骤(7):使用多群差分进化优化方法筛选优质种群,具体为步骤(7.1)至步骤(7.4);
步骤(7.1):随机产生初始种群为:
步骤(7.2):对第g代种群中的第i个个体Xi,j(g)进行差分变异操作,产生相应的变异个体Vi,j(g+1);具体操作为:随机选取第g代种群中除了第i个个体Xi,j(g)以外的两个不同个体,将这两个个体的向量差进行缩放后与待变异个体进行向量合成;公式如下:
式中,r1、r2和r3是三个随机数,且区间为[1,NP],NP表示个体总数;g表示种群代数标识; F为缩放因子;
步骤(7.3):对个体Xi,j(g)和变异个体Vi,j(g+1)实施交叉操作,产生相应的试验个体 Ui,j(g+1)为:
式中,CR为交叉概率,通过概率的方式随机生成新个体;
步骤(7.4):对试验个体Ui,j(g+1)组成的临时种群进行评价,择优选择下一代新种群:
式中,xi,j(g+1)为第g+1代新种群中的个体;
步骤(8):使用动态网络规划求解模型;
状态变量定义为:
决策变量定义为:
状态转移方程定义为:
式中,为设备TS、TR、PV、WT和EV的简写;St+1表示投到t+1阶段的设备资源; p(St+1,n|St,n,At,n,λt,n)表示在t阶段经过决策和负载预测误差分析之后,从状态St,n转移至t+1 阶段的状态空间集St+1,n的概率;表示在t阶段进行决策之后,低碳型柔性配电网拓扑结构由变化到的概率;p(λt+1,n|λt,n)表示在t阶段负载增长的预测误差由λt,n变化到λt+1,n的概率;其中,λt,n服从高斯分布,并且满足条件σ为高斯分布的标准差;
阶段指标函数定义为:
ε表示步长,并且满足ε∈(0,1);
最优指标函数定义为:
步骤(8.1):根据第n次迭代的第t个阶段的设备资源集合St,n,结合模型的目标函数和约束条件确定可行的设备决策集合At,n;
步骤(8.3):经过决策和负载预测误差分析,根据状态转移方程p(St+1,n|St,n,At,n,λt,n)将状态St,n转移至t+1阶段的状态空间集St+1,n中;
步骤(8.5):重复步骤(8.1)至步骤(8.4)直到达到最大迭代次数n+1时,停止迭代;根据最优指标函数输出最优决策值矩阵,由此确定最优低碳型柔性配电网规划方案。
本发明相对于现有技术具有如下的优点及效果:
(1)在能源电力发展趋势和“双碳”目标导向下,建设以新能源为主体的新型电力系统是推动电力清洁低碳发展的必然选择之一。本发明构建低碳型柔性配电网优化扩展模型,支持大规模电动汽车充电桩接入现有电网,支持以风能、太阳能、水能为代表清洁能源机组参与电力系统规划及调控过程,更新优化配电网结构,实现新能源绿色低碳化转型,促进能源电力高质量发展。
(2)本发明在传统电网规划的基础上考虑数字化用户满意度因素,在满足配电网建设需求的同时,实现柔性电力系统中电源、电网、负荷、储能的各个环节与数字化电力用户端的协调互动和协同优化,提高电力设施利用效率,促成配电网基础设施在长期建设过程中成本的合理分摊。
(3)针对负荷增长不确定性和配电网拓扑动态性的配电系统扩展规划问题,本发明提出一种多群差分进化优化的多层分数阶泰勒动态网络规划方法。相比于整数阶微积分而言,分数阶微积分具有历史记忆性,它能够更好地实现任意阶次的自然延伸,更好地表征电力系统配电网络的特点,从而更好地描述电力系统的动力学行为。而电力系统配电网规划问题本质上是大规模且多阶段的优化问题,应用动态网络规划能将此类问题分解为一系列能被轻松解决的单阶段的子问题,并且在优化周期内相继解决这些子问题,达到合理配置电源设备和降低投资风险的作用。
附图说明
图1是本发明方法的低碳型柔性配电网优化扩展流程示意图。
图2是本发明方法的低碳型柔性配电网优化扩展结构图。
具体实施方式
本发明提出的一种低碳型柔性配电网高满意度动态综合规划方法,结合附图详细说明如下:
图1是本发明方法的低碳型柔性配电网优化扩展流程示意图。首先,将数字化信息平台上采集的数据样本进行k-means聚类技术处理,获取原始特征集。然后,定义预处理数据的类型和特征,并且进行归一化处理。接着,将处理过后的数值区间均匀地离散化,根据最大相关最小冗余准则筛选出冗余性最低并且相关性最高的优质数据特征集。最后,将筛选后的优质数据特征集输入至低碳型柔性配电网优化扩展模型中,使用多群差分进化优化的多层分数阶泰勒动态网络规划对该模型进行求解。其中,多层分数阶泰勒动态网络规划的具体操作为:首先,设置模型的初始化参数和最大迭代数,训练负荷增长场景,生成多个负荷增长方案。根据状态空间集,结合目标函数和约束条件确定可行的动作空间。然后,遍历决策空间的每一个决策,求解出当前操作场景中的状态空间的最佳决策。接着,由当前操作场景中的状态空间以及最佳动作公式,得出当前阶段实际负荷的增长情况,确定最优决策。再根据当前阶段的最优决策,考虑预计负荷增长情况和环境不确定性因素,由当前操作场景中的状态转移方程进入下一阶段的状态空间集,并且重复以上操作,得到各个阶段低碳型柔性配电网扩展计划模型的最优决策方案。一直递归迭代直到达到最大迭代次数之后,输出最优决策值矩阵,由此确定最佳的低碳型柔性配电网规划方案。
图2是本发明方法的低碳型柔性配电网优化扩展结构图。首先,根据用户历史用电数据、电源地理分布数据、地区用电经济效益情况和用户满意度数据汇总至数字化信息管理平台,深度挖掘需求侧响应潜力,再将信息整合后纳入低碳型柔性配电网优化扩展模型中。其次,依据配电网实时运行数据,对各类电源和负荷进行能量管理,并将能量管理信息输入到低碳型柔性配电网优化扩展模型中。最后,使用多群差分进化优化的多层分数阶泰勒动态网络规划对模型做出动态规划最优决策,合理布局调峰,完成各类能源的存储与消纳,完成对发电厂及各类分布式发电设备的动态规划。
Claims (1)
1.一种低碳型柔性配电网高满意度动态综合规划方法,其特征在于,该方法将多群差分进化优化的多层分数阶泰勒动态网络规划应用到低碳型柔性配电网规划中,完成对电力负荷以及居民用电满意度的精准预测,依据预测数据对低碳型柔性配电网进行优化,确定最优低碳型柔性配电网规划方案,保证低碳型柔性配电网绿色高效安全经济运行;低碳型柔性配电网高满意度动态综合规划方法在使用过程中的步骤为:
步骤(1):构建低碳型柔性配电网优化扩展模型,旨在保证负荷增长状态下的预期总成本达到最小化,同时数字化的用户满意程度达到最大化,保证低碳型柔性配电网规划策略处在最优解;总规划方案的目标函数f表示为:
式中,w为负荷预测误差指标序号,λt,w为t阶段第w个负荷预测误差指标水平,Ωt为阶段集,α为年利率,y表示规划年序号;为总规划方案的数学期望,Ny为规划阶段的总年数;A、B、C、D和F分别为设备扩容投资成本、设备运行维护成本、设备能源生产成本、设备能源减载成本和用户需求满意程度标识;和分别为t阶段的设备扩容投资成本、设备运行维护成本、设备能源生产成本、设备能源减载成本和用户需求满意程度;其中,c表示规划方案序号,t表示当前处于第t个阶段;i表示当前处在第i个设备节点上;m表示当前处在第m条馈线支路上,k表示已经安装的设备类型序号;
设备扩容投资成本为:
式中,Lm表示第m条支路的建设总长度;TS、TR、FL、PV、WT和EV分别表示变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩标识;ΩTS、ΩTR、ΩFL、ΩPV、ΩWT和ΩEV分别表示变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的候选节点集;KTR、KFL、KPV、KWT和KEV分别表示变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的型号类型集;和分别表示在t阶段下第i个设备节点上已经安装的第k个变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的投资变量;和分别表示在已经安装的设备类型k中,变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的投资成本系数;τTS、τTR、τFL、τPV、τWT和τEV分别表示变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的设备回收率;
设备运行维护成本为:
设备能源生产成本为:
式中,υc为方案c的持续时间;CH和DIS分别表示电动汽车充电桩的充电过程和放电过程标识;和分别为在方案c中变电站的能源购买成本、光伏电站和风力涡轮机的生产成本;和分别为方案c中电动汽车充电桩的充电成本和放电成本;和分别为方案c在t阶段的变电站、光伏电站和风力涡轮机注入的有功功率;和分别为方案c在t阶段的电动汽车充电桩充电时的有功功率和放电时的有功功率;
设备能源减载成本为:
式中,FH表示负荷标识,ΩFH表示负荷节点集;和分别为方案c在t阶段的变电站和负荷的视在功率削减量,和分别为方案c在t阶段的光伏电站和风力涡轮机的有功功率削减量;和分别为方案c中变电站、负荷、光伏电站和风力涡轮机的削减成本;
用户需求满意程度为:
式中,和分别为方案c中变电站、变压器、光伏电站、风力涡轮机和电动汽车充电桩拟定选址的用户需求满意程度系数;和分别表示第c个规划方案的在t阶段下第i个设备节点上变电站、变压器、馈线、光伏电站、风力涡轮机和电动汽车充电桩的投资变量;
低碳型柔性配电网规划是对现有变电站、变压器以及输电线路进行改造或者新建,对太阳能光伏板、风力发电机和电动汽车充电桩进行选址,对各类设备的约束为:
为了保证新建输电线路不在两个及两个以上的方向同时规划,约束为:
为了保证低碳型柔性配电网有功功率和无功功率的平衡,对各个节点的功率限制为:
式中,和分别为以节点i为线路首端和以节点i为线路末端的馈线集,为方案c在t阶段流经馈线的有功功率,和分别为方案c在t阶段的电动汽车充电桩充电时和放电时的有功功率,和分别为方案c在t阶段的节点预期和实际需求的有功功率;μc为方案c中的需求水平因子,ψc为方案c中的负载功率因数;为方案c在t阶段流经馈线的无功功率;和分别为方案c在t阶段的变电站、光伏电站和风力涡轮机注入的无功功率,和分别为方案c在t阶段的节点预期和实际消耗的无功功率;
馈线在决策后与决策前的视在功率关系为:
馈线的视在功率与有功功率和无功功率的关系为:
对馈线的功率限制为:
变电站在决策后与决策前的视在功率关系为:
变电站的视在功率与有功功率和无功功率的关系为:
对变电站的功率限制为:
光伏电站在决策后与决策前的有功功率关系为:
对光伏电站的功率限制为:
风力涡轮机在决策后与决策前的有功功率关系为:
对风力涡轮机的功率限制为:
电动汽车充电桩在决策后与决策前的有功功率关系为:
对电动汽车充电桩的最大充电功率和放电功率限制为:
为了保证低碳型柔性配电网在运行过程中网络枢纽点的电压水平以及网络中各个节点的电压水平满足技术要求,对电压的限制为:35kV及以上供电电压正、负偏差的绝对值之和不超过系统额定电压的10%;10kV及以下三相供电电压允许偏差为系统额定电压的±7%;220V单相供电电压允许偏差为系统额定电压的-10%至+7%;
步骤(2):对目标区域的低碳型柔性配电网进行初步可行性研究,将计划扩建的电源场址制作成调查问卷,以入户访问、召开座谈会以及线上调研的方式,组织开展关于扩建电源选址满意程度的调查工作,统计各电源选址的满意度情况;与此同时,由电力用户智能信息采集终端对电力用户的用电情况,目标区域的用电经济效益以及电源地理分布数据进行实时采集;
步骤(3):在数字化信息平台上传步骤(2)中的采集数据样本,将数据样本进行k-means聚类技术处理;其中,拟合数据和原始数据的误差平方和定义为:
x具体为[日负荷,季度负荷,年负荷,城乡居民用电负荷,工业用电负荷,需求电源点]T;
各个簇类的集聚程度定义为:
式中,RE表示簇类的集聚程度标识,RRE为簇类聚集程度系数,ΩC为总样本集,n为总样本数,假设某样本点xl已经被聚类到簇AA当中,r(xl)表示样本点xl与簇AA内的其他样本点的平均距离,表示样本点xl与除了簇AA以外的其他簇样本点的平均距离;当簇类的集聚程度最大时,聚类效果最优;当K值小于最佳聚类数值时,增大K值,进而加大各个簇类的集聚程度,使得误差平方和急剧下降;当增大K值至K值取到最佳聚类数时,各个簇类的集聚程度变化不明显,误差平方和也趋近于0,此时,数据样本达到了最佳的聚合效果;在确定最佳聚类数值K之后,获得原始数据特征集;
步骤(4):定义预处理数据的类型和特征;遍历原始数据特征集中的每一个数据,记录数据中的最小值ymin和最大值ymax;将ymin和ymax分别映射为数值0和数值1,对原始数据特征集中的其余数值进行归一化处理为:
式中,ynew为归一化处理后的数值,y为原始数据特征集的数值,ymin和ymax原始数据特征集的数据最小值和数据最大值;
步骤(5):将处理过后的数值区间均匀地离散化,以获取每个特征与类别的联合信息熵以及互信息;其中,联合信息熵定义为:
式中,L(b,d)为第b个特征与类别d的联合信息熵;u为区间数,v为类别数,Nb为特征b的区间数量,Nd为总类别数,MLSH为总样本数,MLSH,u,v表示第b个特征落在第u个区间并且类别为v的样本数;
互信息定义为:
式中,H(b;d)为第b个特征与类别d的互信息,MLSH,u为第b个特征落在第u个区间的样本数,MLSH,v为属于第v个类别的样本数;
最大相关性指标定义为:
式中,G(S,d)为最优特征集S与类别d的最大相关性指标;NS最优特征集S中特征的数量;
根据最大相关最小冗余准则对数据进行分析,筛选出冗余性最低并且相关性最高的优质数据特征集;最大相关最小冗余准则定义为:
步骤(6):将筛选后的优质数据特征集输入至步骤(1)建立的低碳型柔性配电网优化扩展模型中,使用多群差分进化优化的多层分数阶泰勒动态网络规划对该模型进行求解;其中,目标函数f(x)中的NTaylor阶泰勒展开为:
式中,f(x)是关于x=a的NTaylor阶泰勒多项式,c介于x与a之间;其中,NTaylor为整数;
利用分数阶微积分将泰勒级数展开成多项式,此时的函数f(x)能进行任意次的微分,公式为:
由公式(29)和公式(30)重新定义步骤(1)中模型的目标函数:
步骤(7):使用多群差分进化优化方法筛选优质种群,具体为步骤(7.1)至步骤(7.4);
步骤(7.1):随机产生初始种群为:
步骤(7.2):对第g代种群中的第i个个体Xi,j(g)进行差分变异操作,产生相应的变异个体Vi,j(g+1);具体操作为:随机选取第g代种群中除了第i个个体Xi,j(g)以外的两个不同个体,将这两个个体的向量差进行缩放后与待变异个体进行向量合成;公式如下:
式中,r1、r2和r3是三个随机数,且区间为[1,NP],NP表示个体总数;g表示种群代数标识;F为缩放因子;
步骤(7.3):对个体Xi,j(g)和变异个体Vi,j(g+1)实施交叉操作,产生相应的试验个体Ui,j(g+1)为:
式中,CR为交叉概率,通过概率的方式随机生成新个体;
步骤(7.4):对试验个体Ui,j(g+1)组成的临时种群进行评价,择优选择下一代新种群:
式中,xi,j(g+1)为第g+1代新种群中的个体;
步骤(8):使用动态网络规划求解模型;
状态变量定义为:
决策变量定义为:
状态转移方程定义为:
式中,为设备TS、TR、PV、WT和EV的简写;St+1表示投到t+1阶段的设备资源;p(St+1,n|St,n,At,n,λt,n)表示在t阶段经过决策和负载预测误差分析之后,从状态St,n转移至t+1阶段的状态空间集St+1,n的概率;表示在t阶段进行决策之后,低碳型柔性配电网拓扑结构由变化到的概率;p(λt+1,n|λt,n)表示在t阶段负载增长的预测误差由λt,n变化到λt+1,n的概率;其中,λt,n服从高斯分布,并且满足条件σ为高斯分布的标准差;
阶段指标函数定义为:
ε表示步长,并且满足ε∈(0,1);
最优指标函数定义为:
步骤(8.1):根据第n次迭代的第t个阶段的设备资源集合St,n,结合模型的目标函数和约束条件确定可行的设备决策集合At,n;
步骤(8.3):经过决策和负载预测误差分析,根据状态转移方程p(St+1,n|St,n,At,n,λt,n)将状态St,n转移至t+1阶段的状态空间集St+1,n中;
步骤(8.5):重复步骤(8.1)至步骤(8.4)直到达到最大迭代次数n+1时,停止迭代;根据最优指标函数输出最优决策值矩阵,由此确定最优低碳型柔性配电网规划方案。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110629893.6A CN113312779B (zh) | 2021-06-07 | 2021-06-07 | 一种低碳型柔性配电网高满意度动态综合规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110629893.6A CN113312779B (zh) | 2021-06-07 | 2021-06-07 | 一种低碳型柔性配电网高满意度动态综合规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113312779A CN113312779A (zh) | 2021-08-27 |
CN113312779B true CN113312779B (zh) | 2022-11-08 |
Family
ID=77377820
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110629893.6A Active CN113312779B (zh) | 2021-06-07 | 2021-06-07 | 一种低碳型柔性配电网高满意度动态综合规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113312779B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116452010B (zh) * | 2023-03-06 | 2024-04-12 | 国网河北省电力有限公司正定县供电分公司 | 针对城区的配电网扩张规划平台与规划方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8909950B1 (en) * | 2010-04-18 | 2014-12-09 | Aptima, Inc. | Systems and methods of power management |
EP3203392A1 (en) * | 2016-02-03 | 2017-08-09 | Technische Universität München | Common rank approximation in distribution grid probabilistic simulation |
CN108074021A (zh) * | 2016-11-10 | 2018-05-25 | 中国电力科学研究院 | 一种配电网风险辨识系统及方法 |
CN106711966B (zh) * | 2016-12-12 | 2019-01-08 | 云南电网有限责任公司玉溪供电局 | 配电网保护定值整定方法及系统 |
US11178625B2 (en) * | 2017-06-06 | 2021-11-16 | Supply, Inc. | Method and system for wireless power delivery |
CN111146821A (zh) * | 2019-12-31 | 2020-05-12 | 国网浙江省电力有限公司嘉兴供电公司 | 一种计及光伏不确定性的dstatcom优化配置方法 |
CN112217202B (zh) * | 2020-09-29 | 2022-06-14 | 东南大学 | 考虑灵活性投资的分布式新能源、储能与配电网规划方法 |
-
2021
- 2021-06-07 CN CN202110629893.6A patent/CN113312779B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN113312779A (zh) | 2021-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109508857B (zh) | 一种主动配电网多阶段规划方法 | |
CN107301470B (zh) | 一种配电网扩展规划与光储选址定容的双层优化方法 | |
CN108306303B (zh) | 一种考虑负荷增长和新能源出力随机的电压稳定评估方法 | |
CN110797919B (zh) | 基于Wasserstein距离和分布鲁棒优化的清洁能源电源规划方法 | |
CN111092451B (zh) | 一种基于配电网网格的光伏消纳预警方法 | |
CN114069687B (zh) | 一种计及逆变器无功调节作用的分布式光伏规划方法 | |
CN115062930A (zh) | 一种配电网源网荷储协调规划方法 | |
CN116523240A (zh) | 一种基于mcmc方法预测的微电网评价反馈调度方法 | |
Du et al. | Deep power forecasting model for building attached photovoltaic system | |
CN108075471B (zh) | 基于随机性电源出力预测的多目标约束优化电网调度策略 | |
CN116822719A (zh) | 一种配电网的多目标规划方法及装置 | |
CN113312779B (zh) | 一种低碳型柔性配电网高满意度动态综合规划方法 | |
CN115833105A (zh) | 一种基于集群划分的配电网规划方法 | |
CN113690930B (zh) | 基于nsga-iii算法的分布式光伏电源中长期选址定容方法 | |
CN111724064A (zh) | 一种基于改进免疫算法的含储能配电网规划方法 | |
CN109829599B (zh) | 基于高比例可再生能源的配电网的集群划分方法及装置 | |
CN108694475B (zh) | 基于混合模型的短时间尺度光伏电池发电量预测方法 | |
CN117748444A (zh) | 一种配电系统的运行模拟方法 | |
CN110867902A (zh) | 基于发电预测的微电网分布式电源去中心优化运行方法 | |
Ma et al. | Distributed control of battery energy storage system in a microgrid | |
CN115395577A (zh) | 基于Copula和K-L散度的综合能源配电网分布鲁棒优化方法 | |
CN114329857A (zh) | 一种基于改进鲸鱼算法的分布式电源规划方法 | |
CN112116131A (zh) | 计及碳排放的综合能源系统多层级优化方法 | |
Ghiassi et al. | On the use of AI as a requirement for improved insolation forecasting accuracy to achieve optimized PV utilization | |
Wu et al. | A distribution network flexible resource capacity configuration method with large renewable energy sources access |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |