CN110490421A - 一种基于模糊c-均值聚类的微电网中小水电站容量配置方法 - Google Patents
一种基于模糊c-均值聚类的微电网中小水电站容量配置方法 Download PDFInfo
- Publication number
- CN110490421A CN110490421A CN201910661007.0A CN201910661007A CN110490421A CN 110490421 A CN110490421 A CN 110490421A CN 201910661007 A CN201910661007 A CN 201910661007A CN 110490421 A CN110490421 A CN 110490421A
- Authority
- CN
- China
- Prior art keywords
- matrix
- fuzzy
- micro
- minimum
- capacitance sensor
- 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 34
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 238000013480 data collection Methods 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims abstract description 3
- 239000011159 matrix material Substances 0.000 claims description 72
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 23
- 238000007621 cluster analysis Methods 0.000 claims description 6
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 4
- 238000010248 power generation Methods 0.000 abstract description 19
- 230000005611 electricity Effects 0.000 description 17
- 241000196324 Embryophyta Species 0.000 description 10
- 230000001932 seasonal effect Effects 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 3
- 240000002853 Nelumbo nucifera Species 0.000 description 2
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- -1 Electricity Substances 0.000 description 1
- 241000452036 Hydrometra stagnorum Species 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000153 supplemental effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- 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
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Economics (AREA)
- Data Mining & Analysis (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Tourism & Hospitality (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Primary Health Care (AREA)
- Probability & Statistics with Applications (AREA)
- Power Engineering (AREA)
- Educational Administration (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于模糊C‑均值聚类的微电网中小水电站容量配置方法,包括下述步骤:数据集的构建;从多年历史和实时数据集中获取微电网中水电站每年1月‑12月流域降雨量、来水量、径流量、入库流量和出库流量的数据,通过处理、计算和分析,构建水电站1月‑12月包含流域降雨量、来水量、径流量、入库流量和出库流量五个特征值的数据集X={X1,X2,...,Xn};本发明可计算小水电微电网中小水电站的容量配置,这种小水电微电网电源容量配置方法,同时反映了多年1‑12月入库流量变化的随机性,为小水电微电网电源容量配置、发电出力预测、运行调度提供理论指导,为分布式新能源发电及智能电网调度运行提供必要的技术支撑。
Description
技术领域
本发明涉及电力系统及其自动化技术领域,具体涉及一种基于模糊C-均值聚类的微电网中小水电站容量配置方法。
背景技术
微电网是在其内部以一定方式集成分布式源(小水电、小风电、光伏发电)-荷(水、电、气、冷、热负荷)的一种电网形式。微电网以380V、10kV、35kV等电压等级与主电网联接,在一般正常运行条件下与主电网并网运行,在重负荷时微电网从主电网吸收功率,在轻负荷时可能向主电网注入功率;在主电网局部故障情况下或在邻近微电网故障情况下可以孤网运行,在确保电能质量的前提下由微电网内部分布式电源向负荷提供电力电量,实现无故障微电网正常的供电状态,减小停电时间,提高供电可靠性。
微电网构建和运行的目的是可持续高效利用/消纳微电网内部分布式电源电力电量,最小化与主电网交换的电力电量。
基于分布式小水电的微电网,是一种以小水电供电为主要形式的微电网。小水电微电网中,水电站大多数为径流式,水坝一般没有蓄水功能,水库没有蓄水和调水能力,小水电站水能利用完全取决于水库来水量,小水电机组发电状态与出力规模也完全取决于水库来水量。在这种情况下,小水电站要实现高效利用水能发电,就必须做到来多少水发多少电。而小水电站水库来水量具有随机性,在不同水文周期来水量完全不同,在丰水期来水量很大,在枯水期来水量很小。因此,小水电站流域河流流量往往表现为最小流量、最大流量、平均流量、多年平均流量、计算平均流量、加权平均流量、数学平均流量等等表式形式。采用不同流量的表式形式,小水电站会获得不同装机容量水平。不同装机容量水平,在不同水文周期小水电站发电功率和发电量也往往不同,最优导致小水电站水能利用率、发电设备利用率、发电设备年最大利用小时数也不同。
小水电-风微电网是一种集成小水电和小风电两种分布式电源、一定容量负荷并以一定方式将分布式电源和负荷联接起来的微电网。小水电-风微电网中,不仅小水电站水库来水量、水库流量、发电流量等具有不确定性和随机性,而且小风电场风速具有不确定性和随机性。当风速小于切入风速或大于切出风速,风电机组均无输出功率;当风速大于切入风速而小于额定风速时,风电机组输出功率小于额定功率;当风速大于额定风速而小于切出风速时,风电机组输出额定功率。在一年的不同季节风速完全不同,在一天的不同时段风速完全不同,风速具有随机性、波动性和间歇性。因此,小风电场风速往往也表现为最小风速、最大风速、平均风速、多年平均风速、计算平均风速、加权平均风速、数学平均风速等等表式形式。采用不同风速的表式形式,小风电场会获得不同装机容量水平。不同装机容量水平,在不同季节小风电场发电功率和发电量也往往不同,最优导致小风电站风能利用率、发电设备利用率、发电设备年最大利用小时数也不同。
小水电-风-光微电网是一种集成小水电、小风电、光伏发电三种分布式电源、一定容量负荷并以一定方式将分布式电源和负荷联接起来的微电网。小水电-风-光微电网中,不仅小水电站水库来水量、水库流量、发电流量等具有不确定性和随机性,小风电场风速具有不确定性和随机性,而且日照强度也具有不确定性和随机性。日照强度越大,光伏发电系统输出功率越大。在一年的不同季节日照强度完全不同,在一天的不同时段日照强度完全不同,日照强度具有随机性、波动性和间歇性。因此,光伏发电站日照强度往往也表现为最小日照强度、最大日照强度、平均日照强度、多年平均日照强度、计算平均日照强度、加权平均日照强度、数学平均日照强度等等表式形式。采用不同日照强度的表式形式,光伏发电站会获得不同装机容量水平。不同装机容量水平,在不同季节光伏发电站发电功率和发电量也往往不同,最优导致光伏发电站风能利用率、发电设备利用率、发电设备年最大利用小时数也不同。
不同负荷水平及分布式电源容量规模在微电网中集成,使微电网结构形态、潮流特性发生了改变。由于小水电、小风电、光伏发电等多种分布式电源接入,会因接入电源容量规模的不同而采用多种不同等级的电压。由于用电的随机性,负荷功率在不同时空尺度上总会发生变化,具有明显的时段性;同时,风电、光伏发电等分布式电源出力具有间歇性、随机性和时段性,小水电机组出力具有季节性。因此,微电网负荷功率和电源功率平衡关系往往很难能够保持,当负荷功率大于电源功率时微电网需要从主电网获得补充功率,而负荷功率小于电源功率时微电网剩余功率需要注入主电网,形成了随机的双向潮流特性。随机的双向潮流特性对于微电网内部节点电压质量,在分布式电源出力大而负荷轻时就会造成微电网内部局部区域节点电压偏高,而在分布式电源出力小而负荷重时就会造成微电网内部局部区域节点电压偏低。因此,微电网内部节点电压的限制条件和要求对微电网中分布式电源容量配置、运行方式和电压控制策略均有影响和制约,微电网中分布式电源容量配置、运行方式和电压控制策略需要考虑微电网内部节点电压的限制条件和要求。微电网在不同电压等级配电网的节点接入,会因微电网从配电网吸收或向配电网注入功率的大小不同而导致配电网节点电压发生偏高或偏低的变化,微电网中分布式电源容量配置、运行方式和电压控制策略需要考虑配电网节点电压的限制条件和要求。
微电网分布式电源系统是一种同时具有关系复杂并交互作用的随机和模糊不确定性事件或参量的系统。在各种不确定性随机和模糊事件或参量的影响下,微电网分布式电源发电功率和发电量变得更具随机特性和模糊特性,这些特性极大影响微电网分布式电源容量配置。以往微电网分布式电源系统发电功率和发电量通常采用确定性的计算方法,有些也采用概率分析的不确定性计算方法。确定性计算的方法通常是在假设小水电站来水量和流量、区域内日照强度和风速都确定的情况下计算微电网分布式电源系统发电功率、发电量及装机容量,也没有考虑微电网和配电网调压要求和柔性控制方式等因素的影响,计算结果是唯一性和确定性的,往往不能反应微电网分布式电源系统发电功率、发电量及装机容量的实际情况。而概率分析的计算方法通常是在只假设小水电站来水量和流量、区域内日照强度和风速等单一因素为不确定性因素的情况下计算微电网分布式电源系统发电功率、发电量及装机容量,计算结果是具有一定置信水平的概率值。实际上,微电网分布式电源系统发电功率、发电量及装机容量受多种不确定性因素影响。而且,这些影响因素通常都有随机不确定性或模糊不确定性,或者他是具有随机和模糊不确定性,往往以随机和模糊不确定性事件或参量而存在。可见,微电网分布式电源系统发电功率、发电量及装机容量计算的现有技术都没有全面考虑影响因素的不确定性和随机性,计算方法适用性、实用性和应用性也难以得到满足。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种基于模糊C-均值聚类的微电网中小水电站容量配置方法,该方法针对由小水电、小风电和光伏发电组成的微电网系统,并考虑小水电站入库流量、风电场风速、光伏发电站日照强度的不确定性和随机性,还考虑日照时间、日照阴影、日照偏角等的不确定性和随机性,提出基于模糊C-均值聚类的微电网中小水电站容量配置方法。
本发明的目的通过下述技术方案实现:
一种基于模糊C-均值聚类的微电网中小水电站容量配置方法,包括下述步骤:
S1,数据集的构建;从多年历史和实时数据集中获取微电网中水电站每年1月-12月流域降雨量、来水量、径流量、入库流量和出库流量的数据,通过处理、计算和分析,构建水电站1月-12月包含流域降雨量、来水量、径流量、入库流量和出库流量五个特征值的数据集X={X1,X2,...,Xn},月数据集特征向量为:
Xi={xi1,xi2,xi3,xi4,xi5}={YJYL,LLSL,QJLL,QI,QO}
其中YMJYL、LMLSL、QMJLL、QMI、QMO分别为降雨量、来水量、径流量、入库流量和出库流量,其单位分别为:毫米、立方米、立方米、立方米/秒、立方米/秒;数据样本的数量n=12,数据集的数量m=5,i=1,2,...,n;
S2,模糊聚类中心矩阵的构建;采用模糊聚类分析的方法,选取c=3,构建数据集的模糊聚类中心矩阵:C={C1,C2,C3},其中Ci={Ci1,Ci2,...,Cim}(i=1,2,...,c),C1、C2、C3分别为最小量、平均量、最大量模糊聚类中心矩阵;
S3,模糊聚类中心矩阵的初始化;按照微电网中小水电站容量配置原则,根据流域降雨量、来水量、径流量、入库流量和出库流量的数据变化特征,设定最小量、平均量、最大量模糊聚类中心矩阵的初始值;假设流域降雨量、来水量、径流量、入库流量和出库流量的多年平均值为YMJYL、LMLSL、QMJLL、QMI、QMO,那么设定最小量、平均量、最大量模糊聚类中心矩阵的初始值分别为:
其中
kMJYL1=kMLSL1=kMJLL1=kMI1=kMO1=0.2
kMJYL2=kMLSL2=kMJLL2=kMI2=kMO2=1.0
kMJYL3=kMLSL3=kMJLL3=kMI3=kMO3=1.5
模糊聚类中心矩阵的初始值为
S4,最优模糊聚类矩阵和最优模糊聚类中心矩阵的迭代计算;采用模糊聚类分析的方法,通过迭代计算方式,确定模糊C-均值聚类矩阵的最优值,确定最小量、平均量、最大量模糊聚类中心矩阵的最优值;
模糊C-均值聚类矩阵的最优值为:
最小量、平均量、最大量模糊聚类中心矩阵的最优值为:
其中YCMJYL1、YCMJYL2、YCMJYL3分别为小水电站流域降雨量最小量、平均量、最大量模糊聚类中心矩阵的最优值,LCMLSL1、LCMLSL2、LCMLSL3分别为小水电站来水量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMJLL1、QCMJLL2、QCMJLL3分别为小水电站径流量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMI1、QCMI2、QCMI3分别为小水电站入库流量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMO1、QCMO2、QCMO3分别为小水电站出库流量最小量、平均量、最大量模糊聚类中心矩阵的最优值;
S5,入库流量平均值计算;计算微电网中小水电站入库流量平均值:
QI=max(u21,u22,u23,u24,u25,u26,u27,u28,u29,u210,u211,u212)QCMI2
其中max()为取最大值函数;
S6,装机容量计算;计算小水电站小水电发电机组装机容量:
PS=0.0098HQI
其中H为小水电站水头。
优选地,所述S4中最优模糊聚类矩阵和最优模糊聚类中心矩阵的迭代计算,其具体步骤为:
S4.1,设定模糊系数其中n为数据集数据样本的数量;设定迭代计算结束判定误差值ε=0.001;
S4.2,设置降雨量对小水电站最小量、平均量、最大量的影响系数分别为kYc1、kYc2、kYc3,来水量对小水电站最小量、平均量、最大量的影响系数分别为kLc1、kLc2、kLc3,径流量对小水电站最小量、平均量、最大量的影响系数分别为kJc1、kJc2、kJc3,入库流量对小水电站最小量、平均量、最大量的影响系数分别为kIc1、kIc2、kIc3,出库流量对小水电站最小量、平均量、最大量的影响系数分别为kOc1、kOc2、kOc3;构建流域降雨量、来水量、径流量、入库流量、出库流量对最小量、平均量、最大量的影响系数矩阵:
其中j=1,2,...,n;
S4.3,利用随机函数在[0,1]随机产生随机数,以初始化模糊聚类矩阵U(0);
S4.4,设定迭代次数t=1;
S4.5,在第t次迭代中计算模糊聚类矩阵U(t)的元素值,元素值迭代更新公式为:
其中d()为距离函数;
S4.6,在第t次迭代中计算模糊聚类中心矩阵C(t)的元素值,最小量、平均量、最大量模糊聚类中心矩阵元素值迭代更新公式分别为:
其中i=1,2,...,c,j=1,2,...,m;
S4.7,若则迭代计算结束;否则使t=t+1并转到S4.5。
本发明与现有技术相比具有以下的有益效果:
本发明可计算小水电微电网中小水电站的容量配置,这种小水电微电网电源容量配置方法,同时反映了多年1-12月入库流量变化的随机性,为小水电微电网电源容量配置、发电出力预测、运行调度提供理论指导,为分布式新能源发电及智能电网调度运行提供必要的技术支撑。
附图说明
图1为本发明的流程框图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
本发明针对小水电站入库流量、风电场风速、光伏发电站日照强度的不确定性,采用概率分析的方法,计算小水电站入库流量、风电场风速、光伏发电站日照强度按照正态分布规律变化的概率和平均值,计算小水电站、风电场、光伏发电站装机容量初值,根据微电网内部及微电网接入配电网节点电压变化的限制,按照小水电站、风电场、光伏发电站装机容量初值的占比缩减其装机容量的大小,其期望目标是实现新能源利用率、发电效率、发电设备利用率、发电设备年利用小时数最大化,增加在日、月、年、多年等运行周期内的发电量。
小水电站装机容量由小水电站入库流量决定,风电场装机容量由风电场风速决定,光伏发电站装机容量由日照强度决定。由对新能源利用率、发电量、发电设备利用率和年利用小时数决定,由微电网、配电网输电能力、调压要求、网损控制等决定,在不同时间和空间上日照强度、日照时间、日照阴影、日照偏角度及其不确定性和随机性对光伏发电系统装机容量有极大的影响。水、风、光的时空特征相互独立,但是集成在微电网中其容量配置相互影响、相互制约,微电网和配电网输电能力、调压要求、网损控制等使微电网水风光电源容量配置变得更为复杂。
图1中的步骤1描述小水电站数据集构建的过程和方法。从多年历史和实时数据集中获取微电网中水电站每年1月-12月流域降雨量、来水量、径流量、入库流量、出库流量等的数据,通过进行处理、计算和分析,构建水电站1月-12月包含流域降雨量、来水量、径流量、入库流量、出库流量5个特征值的数据集X={X1,X2,...,Xn},月数据集特征向量为:
Xi={xi1,xi2,xi3,xi4,xi5}={YJYL,LLSL,QJLL,QI,QO}
其中YMJYL、LMLSL、QMJLL、QMI、QMO分别为降雨量、来水量、径流量、入库流量、出库流量,其单位分别为:毫米、立方米、立方米、立方米/秒、立方米/秒。数据样本的数量n=12,数据集的数量m=5,i=1,2,...,n。
图1中的步骤2描述模糊聚类中心矩阵的构建的过程和方法。采用模糊聚类分析的方法,选取c=3,构建数据集的模糊聚类中心矩阵:C={C1,C2,C3},其中Ci={Ci1,Ci2,...,Cim}(i=1,2,...,c),C1、C2、C3分别为最小量、平均量、最大量模糊聚类中心矩阵。
图1中的步骤3描述模糊聚类中心矩阵初始化的过程和方法。按照微电网中小水电站容量配置原则,根据流域降雨量、来水量、径流量、入库流量、出库流量等的数据变化特征,设定最小量、平均量、最大量模糊聚类中心矩阵的初始值。假设流域降雨量、来水量、径流量、入库流量、出库流量的多年平均值为YMJYL、LMLSL、QMJLL、QMI、QMO,那么设定最小量、平均量、最大量模糊聚类中心矩阵的初始值分别为:
其中
kMJYL1=kMLSL1=kMJLL1=kMI1=kMO1=0.2
kMJYL2=kMLSL2=kMJLL2=kMI2=kMO2=1.0
kMJYL3=kMLSL3=kMJLL3=kMI3=kMO3=1.5
模糊聚类中心矩阵的初始值为
图1中的步骤4描述最优模糊聚类矩阵和最优模糊聚类中心矩阵的迭代计算的过程和方法。采用模糊聚类分析的方法,通过迭代计算方式,确定C-均值模糊聚类矩阵的最优值,确定最小量、平均量、最大量模糊聚类中心矩阵的最优值。
C-均值模糊聚类矩阵的最优值为:
最小量、平均量、最大量模糊聚类中心矩阵的最优值为:
其中YCMJYL1、YCMJYL2、YCMJYL3分别为小水电站流域降雨量最小量、平均量、最大量模糊聚类中心矩阵的最优值,LCMLSL1、LCMLSL2、LCMLSL3分别为小水电站来水量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMJLL1、QCMJLL2、QCMJLL3分别为小水电站径流量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMI1、QCMI2、QCMI3分别为小水电站入库流量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMO1、QCMO2、QCMO3分别为小水电站出库流量最小量、平均量、最大量模糊聚类中心矩阵的最优值。
其中迭代计算包括如下步骤:
(1)设定模糊系数其中n为数据集数据样本的数量。设定迭代计算结束判定误差值ε=0.001。
(2)设置降雨量对小水电站最小量、平均量、最大量的影响系数分别为kYc1、kYc2、kYc3,来水量对小水电站最小量、平均量、最大量的影响系数分别为kLc1、kLc2、kLc3,径流量对小水电站最小量、平均量、最大量的影响系数分别为kJc1、kJc2、kJc3,出库流量对小水电站最小量、平均量、最大量的影响系数分别为kOc1、kOc2、kOc3。构建流域降雨量、来水量、径流量、入库流量、出库流量对最小量、平均量、最大量的影响系数矩阵:
其中j=1,2,...,n。
(3)利用随机函数在[0,1]随机产生随机数,以初始化模糊聚类矩阵U(0)。
(4)设定迭代次数t=1。
(5)在第t次迭代中计算模糊聚类矩阵U(t)的元素值,元素值迭代更新公式为:
其中d()为距离函数。
(6)在第t次迭代中计算模糊聚类中心矩阵C(t)的元素值,最小量、平均量、最大量模糊聚类中心矩阵元素值迭代更新公式分别为:
其中i=1,2,...,c,j=1,2,...,m。
(7)若则迭代计算结束;否则使t=t+1并转到步骤(5)。
图1中的步骤5描述微电网中小水电站入库流量平均值计算的过程和方法。计算微电网中小水电站入库流量平均值:
QI=max(u21,u22,u23,u24,u25,u26,u27,u28,u29,u210,u211,u212)QCMI2
其中max()为取最大值函数。
图1中的步骤6描述小水电站小水电发电机组装机容量计算的过程和方法。计算小水电站小水电发电机组装机容量:
PS=0.0098HQI
其中H为小水电站水头。
上述为本发明较佳的实施方式,但本发明的实施方式并不受上述内容的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (2)
1.一种基于模糊C-均值聚类的微电网中小水电站容量配置方法,其特征在于,包括下述步骤:
S1,数据集的构建;从多年历史和实时数据集中获取微电网中水电站每年1月-12月流域降雨量、来水量、径流量、入库流量和出库流量的数据,通过处理、计算和分析,构建水电站1月-12月包含流域降雨量、来水量、径流量、入库流量和出库流量五个特征值的数据集X={X1,X2,...,Xn},月数据集特征向量为:
Xi={xi1,xi2,xi3,xi4,xi5}={YJYL,LLSL,QJLL,QI,QO}
其中YMJYL、LMLSL、QMJLL、QMI、QMO分别为降雨量、来水量、径流量、入库流量和出库流量,其单位分别为:毫米、立方米、立方米、立方米/秒、立方米/秒;数据样本的数量n=12,数据集的数量m=5,i=1,2,...,n;
S2,模糊聚类中心矩阵的构建;采用模糊聚类分析的方法,选取c=3,构建数据集的模糊聚类中心矩阵:C={C1,C2,C3},其中Ci={Ci1,Ci2,...,Cim}(i=1,2,...,c),C1、C2、C3分别为最小量、平均量、最大量模糊聚类中心矩阵;
S3,模糊聚类中心矩阵的初始化;按照微电网中小水电站容量配置原则,根据流域降雨量、来水量、径流量、入库流量和出库流量的数据变化特征,设定最小量、平均量、最大量模糊聚类中心矩阵的初始值;假设流域降雨量、来水量、径流量、入库流量和出库流量的多年平均值为YMJYL、LMLSL、QMJLL、QMI、QMO,那么设定最小量、平均量、最大量模糊聚类中心矩阵的初始值分别为:
其中
kMJYL1=kMLSL1=kMJLL1=kMI1=kMO1=0.2
kMJYL2=kMLSL2=kMJLL2=kMI2=kMO2=1.0
kMJYL3=kMLSL3=kMJLL3=kMI3=kMO3=1.5
模糊聚类中心矩阵的初始值为
S4,最优模糊聚类矩阵和最优模糊聚类中心矩阵的迭代计算;采用模糊聚类分析的方法,通过迭代计算方式,确定模糊C-均值聚类矩阵的最优值,确定最小量、平均量、最大量模糊聚类中心矩阵的最优值;
模糊C-均值聚类矩阵的最优值为:
最小量、平均量、最大量模糊聚类中心矩阵的最优值为:
其中YCMJYL1、YCMJYL2、YCMJYL3分别为小水电站流域降雨量最小量、平均量、最大量模糊聚类中心矩阵的最优值,LCMLSL1、LCMLSL2、LCMLSL3分别为小水电站来水量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMJLL1、QCMJLL2、QCMJLL3分别为小水电站径流量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMI1、QCMI2、QCMI3分别为小水电站入库流量最小量、平均量、最大量模糊聚类中心矩阵的最优值,QCMO1、QCMO2、QCMO3分别为小水电站出库流量最小量、平均量、最大量模糊聚类中心矩阵的最优值;
S5,入库流量平均值计算;计算微电网中小水电站入库流量平均值:
QI=max(u21,u22,u23,u24,u25,u26,u27,u28,u29,u210,u211,u212)QCMI2
其中max()为取最大值函数;
S6,装机容量计算;计算小水电站小水电发电机组装机容量:
PS=0.0098HQI
其中H为小水电站水头。
2.根据权利要求1所述的基于模糊C-均值聚类的微电网中小水电站容量配置方法,其特征在于,所述S4中最优模糊聚类矩阵和最优模糊聚类中心矩阵的迭代计算,其具体步骤为:
S4.1,设定模糊系数其中n为数据集数据样本的数量;设定迭代计算结束判定误差值ε=0.001;
S4.2,设置降雨量对小水电站最小量、平均量、最大量的影响系数分别为kYc1、kYc2、kYc3,来水量对小水电站最小量、平均量、最大量的影响系数分别为kLc1、kLc2、kLc3,径流量对小水电站最小量、平均量、最大量的影响系数分别为kJc1、kJc2、kJc3,入库流量对小水电站最小量、平均量、最大量的影响系数分别为kIc1、kIc2、kIc3,出库流量对小水电站最小量、平均量、最大量的影响系数分别为kOc1、kOc2、kOc3;构建流域降雨量、来水量、径流量、入库流量、出库流量对最小量、平均量、最大量的影响系数矩阵:
其中j=1,2,...,n;
S4.3,利用随机函数在[0,1]随机产生随机数,以初始化模糊聚类矩阵U(0);
S4.4,设定迭代次数t=1;
S4.5,在第t次迭代中计算模糊聚类矩阵U(t)的元素值,元素值迭代更新公式为:
其中d()为距离函数;
S4.6,在第t次迭代中计算模糊聚类中心矩阵C(t)的元素值,最小量、平均量、最大量模糊聚类中心矩阵元素值迭代更新公式分别为:
其中i=1,2,...,c,j=1,2,...,m;
S4.7,若则迭代计算结束;否则使t=t+1并转到S4.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910661007.0A CN110490421B (zh) | 2019-07-22 | 2019-07-22 | 一种基于模糊c-均值聚类的微电网中小水电站容量配置方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910661007.0A CN110490421B (zh) | 2019-07-22 | 2019-07-22 | 一种基于模糊c-均值聚类的微电网中小水电站容量配置方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110490421A true CN110490421A (zh) | 2019-11-22 |
CN110490421B CN110490421B (zh) | 2022-11-22 |
Family
ID=68547667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910661007.0A Active CN110490421B (zh) | 2019-07-22 | 2019-07-22 | 一种基于模糊c-均值聚类的微电网中小水电站容量配置方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110490421B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111934358A (zh) * | 2020-06-28 | 2020-11-13 | 广东电网有限责任公司 | 一种微电网中光伏发电系统输出功率平均值的计算方法 |
CN112862215A (zh) * | 2021-03-10 | 2021-05-28 | 广东电网有限责任公司 | 一种微电网能源需求预测方法、系统及设备 |
CN114118812A (zh) * | 2021-11-29 | 2022-03-01 | 湖南五凌电力科技有限公司 | 基于改进模糊均值聚类的水电站能效分析评价方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080292467A1 (en) * | 2005-10-13 | 2008-11-27 | Sway As | Direct-Drive Generator/Motor for a Windmill/Hydropower Plant/Vessel Where the Generator/Morot is Configured as a Hollow Profile and a Method to Assemble Such a Windmill/Hydropower Plant |
CN104392284A (zh) * | 2014-11-28 | 2015-03-04 | 大连理工大学 | 一种基于情境分析的大中小水电短期可消纳电量最大协调优化调度方法 |
US20150302332A1 (en) * | 2013-05-03 | 2015-10-22 | Gansu Electric Power Corporation Wind Power Technology Center | Maintenance schedule optimization method for electric power system including large-scale wind power |
CN109636140A (zh) * | 2018-11-27 | 2019-04-16 | 广东电网有限责任公司韶关供电局 | 考虑径流量的微电网中小水电站弃水电量计算方法 |
CN109902915A (zh) * | 2019-01-11 | 2019-06-18 | 国网浙江省电力有限公司 | 一种基于模糊c均值模型的电水气用能行为分析方法 |
-
2019
- 2019-07-22 CN CN201910661007.0A patent/CN110490421B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080292467A1 (en) * | 2005-10-13 | 2008-11-27 | Sway As | Direct-Drive Generator/Motor for a Windmill/Hydropower Plant/Vessel Where the Generator/Morot is Configured as a Hollow Profile and a Method to Assemble Such a Windmill/Hydropower Plant |
US20150302332A1 (en) * | 2013-05-03 | 2015-10-22 | Gansu Electric Power Corporation Wind Power Technology Center | Maintenance schedule optimization method for electric power system including large-scale wind power |
CN104392284A (zh) * | 2014-11-28 | 2015-03-04 | 大连理工大学 | 一种基于情境分析的大中小水电短期可消纳电量最大协调优化调度方法 |
CN109636140A (zh) * | 2018-11-27 | 2019-04-16 | 广东电网有限责任公司韶关供电局 | 考虑径流量的微电网中小水电站弃水电量计算方法 |
CN109902915A (zh) * | 2019-01-11 | 2019-06-18 | 国网浙江省电力有限公司 | 一种基于模糊c均值模型的电水气用能行为分析方法 |
Non-Patent Citations (2)
Title |
---|
吴杰康等: "基于数据包络分析的梯级水电站发电效率优化模型", 《中国电机工程学报》 * |
贺志尧等: "基于Elman神经网络的锦屏一级水电站年平均径流量集合预报研究", 《水电能源科学》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111934358A (zh) * | 2020-06-28 | 2020-11-13 | 广东电网有限责任公司 | 一种微电网中光伏发电系统输出功率平均值的计算方法 |
CN111934358B (zh) * | 2020-06-28 | 2023-05-26 | 广东电网有限责任公司 | 一种微电网中光伏发电系统输出功率平均值的计算方法 |
CN112862215A (zh) * | 2021-03-10 | 2021-05-28 | 广东电网有限责任公司 | 一种微电网能源需求预测方法、系统及设备 |
CN114118812A (zh) * | 2021-11-29 | 2022-03-01 | 湖南五凌电力科技有限公司 | 基于改进模糊均值聚类的水电站能效分析评价方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN110490421B (zh) | 2022-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Eltamaly et al. | Optimal sizing and designing of hybrid renewable energy systems in smart grid applications | |
Zhang et al. | Optimal photovoltaic capacity of large-scale hydro-photovoltaic complementary systems considering electricity delivery demand and reservoir characteristics | |
Lin et al. | Prefeasibility study of a distributed photovoltaic system with pumped hydro storage for residential buildings | |
CN109636140A (zh) | 考虑径流量的微电网中小水电站弃水电量计算方法 | |
CN103218757B (zh) | 一种确定光伏发电容量可信度的方法 | |
Wang et al. | Coordinated operation of conventional hydropower plants as hybrid pumped storage hydropower with wind and photovoltaic plants | |
CN109687506A (zh) | 微电网中小水电站发电量预测方法 | |
CN110490421A (zh) | 一种基于模糊c-均值聚类的微电网中小水电站容量配置方法 | |
Chuang et al. | Analyzing major renewable energy sources and power stability in Taiwan by 2030 | |
Thalamttathu Thankappan et al. | Pico‐hydel hybrid power generation system with an open well energy storage | |
CN108808730A (zh) | 考虑光伏时空分布规律的配电网系统负荷备用容量计算方法及系统 | |
CN108281959A (zh) | 一种高比例可再生能源电力系统的骨干网架优化方法 | |
CN110460090A (zh) | 一种水风光微电网电源容量配置方法 | |
CN110707754B (zh) | 微电网中水风光电源容量配置的优化方法 | |
Yang et al. | The potential for photovoltaic-powered pumped-hydro systems to reduce emissions, costs, and energy insecurity in rural China | |
CN110707682B (zh) | 基于模糊c-均值聚类的微电网中水风光电源容量配置方法 | |
Chisale | Design of stand-alone solar-wind-hydro based hybrid power system: Case of rural village in Malawi | |
Gaddam | Optimal unit commitment using swarm intelligence for secure operation of solar energy integrated smart grid | |
Gauché | Spatial-temporal model to evaluate the system potential of concentrating solar power towers in South Africa | |
CN106557867A (zh) | 适用中长时间尺度电网分析的光伏发电概率模型建模方法 | |
Stucchi et al. | Future hydropower production under the framework of NextGenerationEU: The case of Santa Giustina reservoir in Italian Alps | |
Stocks et al. | Pumped Hydro Energy Storage to support a 50-100% renewable electricity grid | |
CN110889592B (zh) | 一种小水电微电网装机容量配置方法 | |
CN110991798B (zh) | 一种小水电微电网发电设备利用率计算方法 | |
CN110570078B (zh) | 一种基于可能性c-均值聚类的微电网中小水电站发电量计算方法 |
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 |