CN108376262A - 一种风电出力典型特性的分析模型构建方法 - Google Patents
一种风电出力典型特性的分析模型构建方法 Download PDFInfo
- Publication number
- CN108376262A CN108376262A CN201810155688.9A CN201810155688A CN108376262A CN 108376262 A CN108376262 A CN 108376262A CN 201810155688 A CN201810155688 A CN 201810155688A CN 108376262 A CN108376262 A CN 108376262A
- Authority
- CN
- China
- Prior art keywords
- wind power
- power output
- data
- wind
- typical characteristics
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 37
- 238000010276 construction Methods 0.000 title claims abstract description 35
- 230000005611 electricity Effects 0.000 claims abstract description 48
- 238000000034 method Methods 0.000 claims abstract description 20
- 239000006185 dispersion Substances 0.000 claims abstract description 16
- 238000006116 polymerization reaction Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 7
- 238000009434 installation Methods 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 238000012512 characterization method Methods 0.000 abstract description 6
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 abstract description 2
- 230000003044 adaptive effect Effects 0.000 abstract 1
- 238000011156 evaluation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 238000005192 partition Methods 0.000 description 4
- 241001269238 Data Species 0.000 description 3
- 230000001965 increasing effect Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000033228 biological regulation Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 239000002023 wood Substances 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/23211—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with adaptive number of clusters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H02J3/386—
-
- 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
- H02J3/48—Controlling the sharing of the in-phase component
-
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- Water Supply & Treatment (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (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)
- Wind Motors (AREA)
Abstract
本发明公开一种风电出力典型特性的分析模型构建方法,获取原始数据并对数据预处理,采用表征多时空尺度的风电波动性、同时率评价指标对数据进行归一化处理,分段聚合降低维度;基于改进的自适应模糊聚类算法,计算自适应函数α(c),确定最佳聚类数c,对风电出力特性分类;引入变异离散度系数βi,剔除风电出力畸变数据;提出分区加权中位值概率估计法,辨识并提取风电出力日典型特性曲线/区间带。本发明对多空间尺度的大规模风电出力特性,月、季度或年多时间尺度的风电出力特性进行分析计算,提取风电出力日典型特性曲线/区间带,使电力运行人员能更好地把握风电出力规律,为电力系统实时调度、调度计划安排、风电场规划及系统备用配置提供技术参考。
Description
技术领域
本发明涉及电力系统控制技术领域,尤其涉及一种风电场出力日典型特性曲线/区间带分析模型构建方法。
背景技术
随着风力发电技术的发展,风力发电是目前世界上发展最快的可再生能源,装机容量在电网中所占比逐日俱增,风电出力具有间歇性和随机性,与负荷需求不同步、反调峰的特点,风电大规模并网对电网产生的影响也越来越明显,例如对电网安全稳定、风电送出、调频调峰、电能质量、备用安排、运行单位众多、协调困难等问题的影响,不仅影响到电网的安全运行,也限制了电网接纳风电的能力,导致弃风量进一步增大,造成能源浪费。
由于风电出力的间歇性不确定性,风电出力特性与负荷需求不一样,采用常规出力指标(日最大值/最小值、日平均值、日峰谷差、日峰谷差率等)难以直观得到风电出力特性曲线的大概轮廓。而风电出力特性的评价指标体系是风电并网研究的基础,目前风电出力特性的研究大多为定性分析,缺乏系统性的风电出力特性指标体系,从而导致难以深入掌握大规模风电出力的规律和典型特性,缺少对风电场出力典型特性模型构建方法,导致电力工作者无法更好地把握风电出力规律,电力系统实时调度、调度计划安排、风电场规划及系统备用配置缺失依据标准。
发明内容
本发明的目的在于提供一种风电出力典型特性的分析模型构建方法,综合考虑了风电出力的不确定性,波动性及随机性等不同表现形式,更益于识别风电出力日典型特性,保留了风电出力波动过程的完整性与连续性,协助电力工作者更好地把握风电出力规律,为电力系统实时调度、调度计划安排、风电场规划及系统备用配置提供技术参考。
为实现上述目的,本发明的技术方案为:一种风电出力典型特性的分析模型构建方法,所述构建方法包括以下步骤:
1)获取风电出力的样本数据,计算风电出力特性评价指标及空间特性指标并对其进行归一化处理;
2)依据电价表的峰平谷时段,对风电出力的样本数据分段聚合降维处理;
3)对分段后的数据集,初始化模糊聚类算法的参数,获得初始聚类中心和模糊矩阵;
4)根据聚类数的自适应函数α(c),迭代计算,找出最大的α(c),获得最佳模糊聚类数c;
5)根据最佳模糊聚类数,获得样本数据的模糊离散度及变异系数倒数,并计算出样本数据的变异离散度系数βi;
6)根据变异离散度系数βi,设选出当|βi|>ε时的样本,得到畸变样本数据i,删除畸变样本数据;
7)对删除畸变样本数据后的风电出力数据,采用分区加权中位值概率估计法,获取各采样点的典型特性点,获得风电出力日典型特性曲线;
8)计算j采样点的风电出力值pw,j与典型特性点值pc,j的差值,记为差值pw-c,j,再计算j采样点风电出力日典型特性区间带的上界pup,j和下界pdown,j,获得j采样点的典型特性区间带,然后获得风电出力日典型特性区间带。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤1)中,特性评价指标为波动性指标,波动性指标获取的方法为,获取风电出力变化率:风电出力变化即相隔某一时段的两个点出力之差占风电总装机容量的百分比:
其中:P(t+T)、P(t)为t+T、t时刻的风电出力;Ptotal为风电的总装机容量;
一段时间内风电的波动率为:
其中,n表示采样点的个数。风电出力波动性按时间尺度划分,有秒级、分钟级、小时级波动;按空间尺度划分有单台风电机组、风电场/群、区域波动等。风电波动率反映了风功率变化的不平稳程度,数值越大,则其平稳性越差,数值越小,则其平稳性越好,当E(t)为0时风功率最为平稳,此时风电出力大小保持不变。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤1)中,空间特性指标为同时率,同时率为统计时段内多个风电场有功功率和总额定容量之比,用于确定区域电网的风电接纳能力,同时率的确定公式为:
其中:Pimax为第i个风电场或第i日风电出力最大值;PiN为风电并网装机容量;
风电出力同时率反映各风电场间或同一风电场各天间出力是否同步增加或减小,同时率越高,风电系统对电网的影响越大。
所述步骤1)中,归一化处理方式为样本数据原值减去均值除以标准差,表示为
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤2)中,风电出力的样本数据分段聚合降维处理方式为采用低维度数据近似表示高维度数据的基本形态,对于一个长度为n的风电数据,将其划分为若干段长度为m的风电数据,通过将长度为n的风电数据转变为长度为m的风电数据,实现数据降维。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤3)中,模糊聚类算法是将给定的有限数据对象集合X={x1,x2,...,xn},xi∈Rs,数据集合中每个样本是s维向量,把数据对象聚合成c类,各聚类中心向量矩阵为V={v1,v2,...,vc},vj∈Rs,引入一个模糊划分隶属度矩阵(uij)c×n,采用模糊划分表示各样本属于各聚类中心(vij)s×c的程度,经过迭代计算使目标函数最小,隶属度U(k)和聚类中心矩阵V(k+1)确定公式分别为:
其中,q∈{1,+∞}为模糊权重因子,控制分类隶属度矩阵的模糊度,q越大分类的模糊度越高。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤4)中,自适应函数α(c)确定公式为:
对分段后的数据集,初始化模糊聚类算法的参数,c=2,求出初始聚类中心和模糊矩阵。计算自适应函数α(c),通过返回c,如果c<cmax,则c=c+1返回所述步骤3,若干次迭代计算,找出最大的α(c),与之对应有聚类数即为最佳模糊聚类数c。应用模糊聚类算法对给定的数据集进行聚类分析时,关键是样本集的最佳聚类数c的选取和模糊权重因子q的确定。
基于“聚类的目的是将数据分类,并尽量使类间距离d1尽可能的大,而类内数据点距离d2尽可能的小”的基本思想,引入聚类数的自适应函数α(c),根据迭代计算结果,找到最大的α(c),返回最佳聚类数c,且当样本数为n时则最佳聚类数的选择范围为模糊权重因子值的选取直接决定分类结果的模糊性,数据集的划分尽可能分明,且聚类结果不能太模糊,即划分模糊度不能太高,优选的,取模糊权重因子值q=2。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤5)中,模糊离散度系数βi确定公式为:
βi=QF(ai)QP(ai)
其中,QF(ai)为变异系数倒数,QP(ai)为模糊离散度,分别表示为:
其中:ui,j是第j类中第i个元素的隶属度,di,j是第j类中第i个元素到该聚类中心的距离,σ为模糊距离的标准差,Di,k为样本i与其它样本k的模糊距离:Di,k=|(ui,j)2(di,j)2-(uk,j)2(dk,j)2|。为了检测不良数据,剔除异常数据,变异离散度系数βi较大的样本数据说明是与数据特性具有弱相关性或者无相关的数据,将此数据样本删除。与聚类中心的远近用模糊离散度QP(ai)表征,不稳定性用变异系数倒数QF(ai)表征。畸变数据存在会导致模糊距离和模糊离散度增大,在各类中,所有异常点都是到聚类中心最远和相对不稳定的点,从而根据计算的βi值可以检测出畸变样本数据。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤7)中,采用分区加权中位值的概率统计法,计算出采样点的典型特性点,分区加权中位值的概率统计法表达式为:
其中,N为区间数;Median(pw,i,max,pw,i,min)为第i区间采样点的中位值;pw,i,max为i区间采样点的最大值;pw,i,min为i区间采样点的最小值;pri,j为采样点分布在i区间内的概率。
如上所述的一种风电出力典型特性的分析模型构建方法,所述步骤8)中,j采样点差值pw-c,j为:
pw-c,j=pw,j-pc,j
j采样点风电出力日典型特性区间带上界pup,j和下界pdown,j为:
其中:N+、N-分别为差值大于零、小于零的区间数,且存在N++N-=N;分别为i区间差值大于零、小于零的中位值;分别为差值在i区间大于零、小于零的概率,依次计算各采样点的区间带,获得风电出力日典型特性区间带。
本发明具有如下优点:与现有技术相比,综合考虑了风电出力的不确定性,波动性及随机性等不同表现形式,更益于识别风电出力日典型特性,保留了风电出力波动过程的完整性与连续性;提取风电出力日典型特性区间带来近似表征一定时间段的风电场/群出力特性,与单根风电出力日典型特性曲线对比,可以更好地表征多时空尺度的风电出力日典型特性,能够协助电力工作者更好地把握风电出力规律,为电力系统实时调度、调度计划安排、风电场规划及系统备用配置提供技术参考。
附图说明
图1为风电出力典型特性的分析模型构建方法流程示意图;
图2为本发明实施例中负荷峰段时风电出力的聚类分布及各天的变异离散度系数示意图;
图3为本发明实施例中负荷谷段时风电出力的聚类分布及各天的变异离散度系数示意图;
图4为本发明实施例中负荷平段时风电出力的聚类分布及各天的变异离散度系数示意图;
图5为本发明实施例中峰段(19:30-00:30)风电出力分类结果示意图;
图6为本发明实施例中1月份风电出力日典型特性曲线图;
图7为本发明实施例中风电出力日典型特性区间带示意图;
图8为本发明实施例中风电出力曲线比较图。
具体实施方式
以下实施例用于说明本发明,但不用来限制本发明的范围。
一种风电出力典型特性的分析模型构建方法,所述构建方法包括以下步骤:
1)获取风电出力的样本数据,计算风电出力特性评价指标及空间特性指标并对其进行归一化处理;
2)依据电价表的峰平谷时段,对风电出力的样本数据分段聚合降维处理;
3)对分段后的数据集,初始化模糊聚类算法的参数,获得初始聚类中心和模糊矩阵;
4)根据聚类数的自适应函数α(c),迭代计算,找出最大的α(c),获得最佳模糊聚类数c;
5)根据最佳模糊聚类数,获得样本数据的模糊离散度及变异系数倒数,并计算出样本数据的变异离散度系数βi;
6)根据变异离散度系数βi,设选出当|βi|>ε时的样本,得到畸变样本数据i,删除畸变样本数据;
7)对删除畸变样本数据后的风电出力数据,采用分区加权中位值概率估计法,获取各采样点的典型特性点,获得风电出力日典型特性曲线;
8)计算j采样点的风电出力值pw,j与典型特性点值pc,j的差值,记为差值pw-c,j,再计算j采样点风电出力日典型特性区间带的上界pup,j和下界pdown,j,获得j采样点的典型特性区间带,然后获得风电出力日典型特性区间带。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤1)中,特性评价指标为波动性指标,波动性指标获取的方法为,获取风电出力变化率:风电出力变化即相隔某一时段的两个点出力之差占风电总装机容量的百分比:
其中:P(t+T)、P(t)为t+T、t时刻的风电出力;Ptotal为风电的总装机容量;
一段时间内风电的波动率为:
其中,n表示采样点的个数。风电出力波动性按时间尺度划分,有秒级、分钟级、小时级波动;按空间尺度划分有单台风电机组、风电场/群、区域波动等。风电波动率反映了风功率变化的不平稳程度,数值越大,则其平稳性越差,数值越小,则其平稳性越好,当E(t)为0时风功率最为平稳,此时风电出力大小保持不变。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤1)中,空间特性指标为同时率,同时率为统计时段内多个风电场有功功率和总额定容量之比,用于确定区域电网的风电接纳能力,同时率的确定公式为:
其中:Pimax为第i个风电场或第i日风电出力最大值;PiN为风电并网装机容量;风电出力同时率反映各风电场间或同一风电场各天间出力是否同步增加或减小,同时率越高,风电系统对电网的影响越大。
所述步骤1)中,归一化处理方式为样本数据原值减去均值除以标准差,表示为
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤2)中,风电出力的样本数据分段聚合降维处理方式为采用低维度数据近似表示高维度数据的基本形态,对于一个长度为n的风电数据,将其划分为若干段长度为m的风电数据,通过将长度为n的风电数据转变为长度为m的风电数据,实现数据降维。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤3)中,模糊聚类算法是将给定的有限数据对象集合X={x1,x2,...,xn},xi∈Rs,数据集合中每个样本是s维向量,把数据对象聚合成c类,各聚类中心向量矩阵为V={v1,v2,...,vc},vj∈Rs,引入一个模糊划分隶属度矩阵(uij)c×n,采用模糊划分表示各样本属于各聚类中心(vij)s×c的程度,经过迭代计算使目标函数最小,隶属度U(k)和聚类中心矩阵V(k+1)确定公式分别为:
其中,q∈{1,+∞}为模糊权重因子,控制分类隶属度矩阵的模糊度,q越大分类的模糊度越高。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤4)中,自适应函数α(c)确定公式为:
对分段后的数据集,初始化模糊聚类算法的参数,c=2,求出初始聚类中心和模糊矩阵。计算自适应函数α(c),通过返回c,如果c<cmax,则c=c+1返回所述步骤3,若干次迭代计算,找出最大的α(c),与之对应有聚类数即为最佳模糊聚类数c。应用模糊聚类算法对给定的数据集进行聚类分析时,关键是样本集的最佳聚类数c的选取和模糊权重因子q的确定。
基于“聚类的目的是将数据分类,并尽量使类间距离d1尽可能的大,而类内数据点距离d2尽可能的小”的基本思想,引入聚类数的自适应函数α(c),根据迭代计算结果,找到最大的α(c),返回最佳聚类数c,且当样本数为n时则最佳聚类数的选择范围为模糊权重因子值的选取直接决定分类结果的模糊性,数据集的划分尽可能分明,且聚类结果不能太模糊,即划分模糊度不能太高,本实施例中模糊权重因子q=2。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤5)中,模糊离散度系数βi确定公式为:
βi=QF(ai)QP(ai)
其中,QF(ai)为变异系数倒数,QP(ai)为模糊离散度,分别表示为:
其中:ui,j是第j类中第i个元素的隶属度,di,j是第j类中第i个元素到该聚类中心的距离,σ为模糊距离的标准差,Di,k为样本i与其它样本k的模糊距离:Di,k=|(ui,j)2(di,j)2-(uk,j)2(dk,j)2|。为了检测不良数据,剔除异常数据,变异离散度系数βi较大的样本数据说明是与数据特性具有弱相关性或者无相关的数据,将此数据样本删除。与聚类中心的远近用模糊离散度QP(ai)表征,不稳定性用变异系数倒数QF(ai)表征。畸变数据存在会导致模糊距离和模糊离散度增大,在各类中,所有异常点都是到聚类中心最远和相对不稳定的点,从而根据计算的βi值可以检测出畸变样本数据。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤7)中,采用分区加权中位值的概率统计法,计算出采样点的典型特性点,分区加权中位值的概率统计法表达式为:
其中,N为区间数;Median(pw,i,max,pw,i,min)为第i区间采样点的中位值;pw,i,max为i区间采样点的最大值;pw,i,min为i区间采样点的最小值;pri,j为采样点分布在i区间内的概率。
风电出力典型特性的分析模型构建方法的一个实施例中,所述步骤8)中,j采样点差值pw-c,j为:
pw-c,j=pw,j-pc,j
j采样点风电出力日典型特性区间带上界pup,j和下界pdown,j为:
其中:N+、N-分别为差值大于零、小于零的区间数,且存在N++N-=N;分别为i区间差值大于零、小于零的中位值;分别为差值在i区间大于零、小于零的概率,依次计算各采样点的区间带,获得风电出力日典型特性区间带。
以下结合新疆地区2015年冬季1月份的风电出力数据通过本发明方法构建风电出力典型特性的分析模型为例对本发明进一步说明。
新疆地区冬季时间长,冬季供热负荷大、供热机组多,调峰灵活性差,冬季新疆风电出力又较大,导致新疆电网消纳不了全部风电,需要弃风处理。本实施例针对新疆乌鲁木齐电网2015年冬季1月份的风电出力数据,以5分钟为一个采样点,每天采样点数为288个,共计8928个数据样本,应用MATLAB软件编写程序计算分析风电出力典型特性。
风电出力随机强,规律性不强,不好对风电出力区间划分,但负荷波动性具有较强的峰平谷规律性,根据国家电网公司新疆电力公司公布的乌鲁木齐最新电价表中的峰平谷时段,如表1所示,将风电出力区间按五个时间段划分,分别称为峰段、谷段及平段。
表1乌鲁木齐负荷峰平谷时段表
依据2015年1月份31天达坂城地区风电出力数据,按风电的同时率及波动率对风电出力进行归一化处理,再进行分段聚合降维近似,对各段的聚类数c初值选为2,模糊权重因子值q取2,计算自适应函数α(c),返回各时段的最优聚类数:峰段最优聚类数为5、谷段最优聚类数为4、平段最优聚类数为2和5。风电出力在各峰段、谷段及平段的聚类分布及变异离散度系数如图2、3、4所示,依据变异离散度系数可确定畸变数据。
根据计算的变异离散度系数βi如图2、3、4右侧,图中的直线为ε,当的数据为畸变数据,将其剔除。从各天的变异离散度系数可知:风电出力在峰段10:00-13:00畸变数据为第2、8、13、14、23五天,在峰段19:30-00:30畸变数据为7、8两天;风电出力在谷段00:35-8:30畸变数据为8、10、11、21、24、27六天;风电出力在平段08:35-09:55畸变数据为15号一天,在平段13:05-19:25畸变数据为7、9、12、19、24五天。
风电出力在峰段19:30-00:30的最优聚类数为5,畸变数据为7、8两天,将其剔除,将此时间段内29天风电出力分类,结果如图5所示(其余各时段的分类结果略)。可见按本方法分类的各类别风电出力形态具有较明显的特征。
计算出剔除畸变数据后各采样点的典型特性点,得到风电出力日典型特性曲线如图6所示。基于风电出力日典型特性曲线,计算每天的各采样点风电出力与风电出力日典型特性点的差值,计算各采样点风电出力日典型特性区间带的上界pup,j和下界pdown,j,获得表征月风电出力日典型特性区间带如图7所示。将1月份31天风电出力及表征该月风电出力日典型特性区间带的对比如图8所示,可见本发明的方法计算分析典型特性区间带可以较好地表征风电出力特性。
虽然,上文中已经用一般性说明及具体实施例对本发明作了详尽的描述,但在本发明基础上,可以对之作一些修改或改进,这对本领域技术人员而言是显而易见的。因此,在不偏离本发明精神的基础上所做的这些修改或改进,均属于本发明要求保护的范围。
Claims (9)
1.一种风电出力典型特性的分析模型构建方法,其特征在于,所述构建方法包括以下步骤:
1)获取风电出力的样本数据,计算风电出力特性评价指标及空间特性指标并对其进行归一化处理;
2)依据电价表的峰平谷时段,对风电出力的样本数据分段聚合降维处理;
3)对分段后的数据集,初始化模糊聚类算法的参数,获得初始聚类中心和模糊矩阵;
4)根据聚类数的自适应函数α(c),迭代计算,找出最大的α(c),获得最佳模糊聚类数c;
5)根据最佳模糊聚类数,获得样本数据的模糊离散度及变异系数倒数,并计算出样本数据的变异离散度系数βi;
6)根据变异离散度系数βi,设选出当|βi|>ε时的样本,得到畸变样本数据i,删除畸变样本数据;
7)对删除畸变样本数据后的风电出力数据,采用分区加权中位值概率估计法,获取各采样点的典型特性点,获得风电出力日典型特性曲线;
8)计算j采样点的风电出力值pw,j与典型特性点值pc,j的差值,记为差值pw-c,j,再计算j采样点风电出力日典型特性区间带的上界pup,j和下界pdown,j,获得j采样点的典型特性区间带,然后获得风电出力日典型特性区间带。
2.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤1)中,特性评价指标为波动性指标,波动性指标获取的方法为,获取风电出力变化率:风电出力变化即相隔某一时段的两个点出力之差占风电总装机容量的百分比:
其中:P(t+T)、P(t)为t+T、t时刻的风电出力;Ptotal为风电的总装机容量;
一段时间内风电的波动率为:
其中,n表示采样点的个数。
3.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤1)中,空间特性指标为同时率,同时率为统计时段内多个风电场有功功率和总额定容量之比,用于确定区域电网的风电接纳能力,同时率的确定公式为:
其中:Pimax为第i个风电场或第i日风电出力最大值;PiN为风电并网装机容量;
所述步骤1)中,归一化处理方式为样本数据原值减去均值除以标准差,表示为
4.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤2)中,风电出力的样本数据分段聚合降维处理方式为采用低维度数据近似表示高维度数据的基本形态,对于一个长度为n的风电数据,将其划分为若干段长度为m的风电数据,通过将长度为n的风电数据转变为长度为m的风电数据,实现数据降维。
5.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤3)中,模糊聚类算法是将给定的有限数据对象集合X={x1,x2,...,xn},xi∈Rs,数据集合中每个样本是s维向量,把数据对象聚合成c类,各聚类中心向量矩阵为V={v1,v2,...,vc},vj∈Rs,引入一个模糊划分隶属度矩阵(uij)c×n,采用模糊划分表示各样本属于各聚类中心(vij)s×c的程度,经过迭代计算使目标函数最小,隶属度U(k)和聚类中心矩阵V(k+1)确定公式分别为:
其中,q∈{1,+∞}为模糊权重因子,控制分类隶属度矩阵的模糊度,q越大分类的模糊度越高。
6.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤4)中,自适应函数α(c)确定公式为:
计算自适应函数α(c),通过返回c,如果c<cmax,则c=c+1返回所述步骤3,若干次迭代计算,找出最大的α(c),与之对应有聚类数即为最佳模糊聚类数c。
7.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤5)中,模糊离散度系数βi确定公式为:
βi=QF(ai)QP(ai)
其中,QF(ai)为变异系数倒数,QP(ai)为模糊离散度,分别表示为:
其中:ui,j是第j类中第i个元素的隶属度,di,j是第j类中第i个元素到该聚类中心的距离,σ为模糊距离的标准差,Di,k为样本i与其它样本k的模糊距离:Di,k=|(ui,j)2(di,j)2-(uk,j)2(dk,j)2|。
8.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤7)中,采用分区加权中位值的概率统计法,计算出采样点的典型特性点,分区加权中位值的概率统计法表达式为:
其中,N为区间数;Median(pw,i,max,pw,i,min)为第i区间采样点的中位值;pw,i,max为i区间采样点的最大值;pw,i,min为i区间采样点的最小值;pri,j为采样点分布在i区间内的概率。
9.根据权利要求1所述的一种风电出力典型特性的分析模型构建方法,其特征在于,所述步骤8)中,j采样点差值pw-c,j为:
pw-c,j=pw,j-pc,j
j采样点风电出力日典型特性区间带上界pup,j和下界pdown,j为:
其中:N+、N-分别为差值大于零、小于零的区间数,且存在N++N-=N;分别为i区间差值大于零、小于零的中位值;分别为差值在i区间大于零、小于零的概率,依次计算各采样点的区间带,获得风电出力日典型特性区间带。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810155688.9A CN108376262B (zh) | 2018-02-23 | 2018-02-23 | 一种风电出力典型特性的分析模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810155688.9A CN108376262B (zh) | 2018-02-23 | 2018-02-23 | 一种风电出力典型特性的分析模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108376262A true CN108376262A (zh) | 2018-08-07 |
CN108376262B CN108376262B (zh) | 2021-08-10 |
Family
ID=63017862
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810155688.9A Expired - Fee Related CN108376262B (zh) | 2018-02-23 | 2018-02-23 | 一种风电出力典型特性的分析模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108376262B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828375A (zh) * | 2018-08-17 | 2018-11-16 | 国网江苏省电力有限公司电力科学研究院 | 基于双回路智能电能表的空调用电监控系统及监控方法 |
CN109524993A (zh) * | 2018-12-19 | 2019-03-26 | 中国农业大学 | 用于中长期优化调度的风电光伏典型周出力场景生成方法 |
CN110111003A (zh) * | 2019-05-08 | 2019-08-09 | 东南大学 | 一种基于改进fcm聚类算法的新能源典型场景构建方法 |
CN111667135A (zh) * | 2020-03-25 | 2020-09-15 | 国网天津市电力公司 | 一种基于典型特征提取的负荷结构解析方法 |
CN111861278A (zh) * | 2020-08-04 | 2020-10-30 | 长沙理工大学 | 一种电力系统的用电负荷峰谷时段划分方法及系统 |
CN112052981A (zh) * | 2020-07-30 | 2020-12-08 | 南方电网科学研究院有限责任公司 | 基于簇类的风电出力波动区间预测方法、装置及存储介质 |
CN112733316A (zh) * | 2019-10-29 | 2021-04-30 | 新奥数能科技有限公司 | 一种能源系统建模的优化方法及装置 |
CN112766733A (zh) * | 2021-01-21 | 2021-05-07 | 山东大学 | 利用改进的K-means算法加速优化调度算法收敛的方法及系统 |
CN113792778A (zh) * | 2021-09-07 | 2021-12-14 | 国网经济技术研究院有限公司 | 一种风电场分区及其出力相关系数计算方法和系统 |
CN117318126A (zh) * | 2023-11-30 | 2023-12-29 | 山东华天电气有限公司 | 一种微网电能质量优化与储能一体化装置及控制方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258142A (zh) * | 2013-06-04 | 2013-08-21 | 哈尔滨工业大学 | 基于天气模式的风电场短期风速预测方法 |
CN104036073A (zh) * | 2014-05-23 | 2014-09-10 | 电子科技大学 | 适于有功功率特性分析的双馈风电场动态等值建模方法 |
CN104574202A (zh) * | 2014-12-09 | 2015-04-29 | 国网青海省电力公司经济技术研究院 | 一种风电场出力特性分析方法 |
CN105005708A (zh) * | 2015-08-13 | 2015-10-28 | 山东大学 | 一种基于ap聚类算法的广义负荷特性聚类方法 |
CN105069236A (zh) * | 2015-08-13 | 2015-11-18 | 山东大学 | 考虑风电场节点空间相关性的广义负荷联合概率建模方法 |
CN105160151A (zh) * | 2015-07-31 | 2015-12-16 | 国家电网公司 | 一种基于模糊K-means算法的电动汽车应急供电方法 |
CN105244922A (zh) * | 2015-11-13 | 2016-01-13 | 东北电力大学 | 限风情况下双馈型风电场群有功出力优化分配方法 |
CN105703396A (zh) * | 2016-03-21 | 2016-06-22 | 广西电网有限责任公司电力科学研究院 | 一种双馈风力发电场多机表征等值建模方法 |
CN105787588A (zh) * | 2016-02-26 | 2016-07-20 | 南京瑞泽启阳信息科技有限公司 | 一种用于提升新能源消纳能力的动态峰谷分时电价方法 |
CN105956720A (zh) * | 2016-06-07 | 2016-09-21 | 中南大学 | 一种基于t-s模糊模型的风电功率超短期预测方法 |
CN106611260A (zh) * | 2015-10-26 | 2017-05-03 | 国家电网公司 | 小区配电网工程造价指标体系评估模型的建立 |
CN106845863A (zh) * | 2017-02-23 | 2017-06-13 | 沈阳工业大学 | 一种分布式风力发电出力与供暖负荷同步指数预测方法 |
CN107679687A (zh) * | 2017-08-29 | 2018-02-09 | 国网江苏省电力公司经济技术研究院 | 一种光伏出力建模方法以及发电系统可靠性评估方法 |
-
2018
- 2018-02-23 CN CN201810155688.9A patent/CN108376262B/zh not_active Expired - Fee Related
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103258142A (zh) * | 2013-06-04 | 2013-08-21 | 哈尔滨工业大学 | 基于天气模式的风电场短期风速预测方法 |
CN104036073A (zh) * | 2014-05-23 | 2014-09-10 | 电子科技大学 | 适于有功功率特性分析的双馈风电场动态等值建模方法 |
CN104574202A (zh) * | 2014-12-09 | 2015-04-29 | 国网青海省电力公司经济技术研究院 | 一种风电场出力特性分析方法 |
CN105160151A (zh) * | 2015-07-31 | 2015-12-16 | 国家电网公司 | 一种基于模糊K-means算法的电动汽车应急供电方法 |
CN105005708A (zh) * | 2015-08-13 | 2015-10-28 | 山东大学 | 一种基于ap聚类算法的广义负荷特性聚类方法 |
CN105069236A (zh) * | 2015-08-13 | 2015-11-18 | 山东大学 | 考虑风电场节点空间相关性的广义负荷联合概率建模方法 |
CN106611260A (zh) * | 2015-10-26 | 2017-05-03 | 国家电网公司 | 小区配电网工程造价指标体系评估模型的建立 |
CN105244922A (zh) * | 2015-11-13 | 2016-01-13 | 东北电力大学 | 限风情况下双馈型风电场群有功出力优化分配方法 |
CN105787588A (zh) * | 2016-02-26 | 2016-07-20 | 南京瑞泽启阳信息科技有限公司 | 一种用于提升新能源消纳能力的动态峰谷分时电价方法 |
CN105703396A (zh) * | 2016-03-21 | 2016-06-22 | 广西电网有限责任公司电力科学研究院 | 一种双馈风力发电场多机表征等值建模方法 |
CN105956720A (zh) * | 2016-06-07 | 2016-09-21 | 中南大学 | 一种基于t-s模糊模型的风电功率超短期预测方法 |
CN106845863A (zh) * | 2017-02-23 | 2017-06-13 | 沈阳工业大学 | 一种分布式风力发电出力与供暖负荷同步指数预测方法 |
CN107679687A (zh) * | 2017-08-29 | 2018-02-09 | 国网江苏省电力公司经济技术研究院 | 一种光伏出力建模方法以及发电系统可靠性评估方法 |
Non-Patent Citations (8)
Title |
---|
MIIN-SHEN YANG等: "A Feature-Reduction Fuzzy Clustering Algorithm Based on Feature-Weighted Entropy", 《IEEE TRANSACTIONS ON FUZZY SYSTEMS》 * |
NEHA BHARILL等: "Enhanced cluster validity index for the evaluation of optimal number of clusters for Fuzzy C-Means algorithm", 《2014 IEEE INTERNATIONAL CONFERENCE ON FUZZY SYSTEMS (FUZZ-IEEE)》 * |
SIMON FONG等: "Improving classification accuracy using Fuzzy Clustering Coefficients of Variations (FCCV) feature selection algorithm", 《2014 IEEE 15TH INTERNATIONAL SYMPOSIUM ON COMPUTATIONAL INTELLIGENCE AND INFORMATICS (CINTI)》 * |
YANLONG WANG等: "Application of clustering technique to electricity customer classification for load forecasting", 《2015 IEEE INTERNATIONAL CONFERENCE ON INFORMATION AND AUTOMATION》 * |
刘汉强等: "局部搜索自适应核模糊聚类方法", 《计算机工程与科学》 * |
吴琴: "风电场多机等值建模及其应用研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
曲直等: "风电功率变化的一致性和互补性量化评估", 《电网技术》 * |
曾博等: "改进自适应模糊C均值算法在负荷特性分类的应用", 《电力系统》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828375A (zh) * | 2018-08-17 | 2018-11-16 | 国网江苏省电力有限公司电力科学研究院 | 基于双回路智能电能表的空调用电监控系统及监控方法 |
CN108828375B (zh) * | 2018-08-17 | 2020-08-14 | 国网江苏省电力有限公司电力科学研究院 | 基于双回路智能电能表的空调用电监控系统及监控方法 |
CN109524993A (zh) * | 2018-12-19 | 2019-03-26 | 中国农业大学 | 用于中长期优化调度的风电光伏典型周出力场景生成方法 |
CN109524993B (zh) * | 2018-12-19 | 2020-06-02 | 中国农业大学 | 用于中长期优化调度的风电光伏典型周出力场景生成方法 |
CN110111003A (zh) * | 2019-05-08 | 2019-08-09 | 东南大学 | 一种基于改进fcm聚类算法的新能源典型场景构建方法 |
CN112733316A (zh) * | 2019-10-29 | 2021-04-30 | 新奥数能科技有限公司 | 一种能源系统建模的优化方法及装置 |
CN112733316B (zh) * | 2019-10-29 | 2024-04-19 | 新奥数能科技有限公司 | 一种能源系统建模的优化方法及装置 |
CN111667135A (zh) * | 2020-03-25 | 2020-09-15 | 国网天津市电力公司 | 一种基于典型特征提取的负荷结构解析方法 |
CN112052981A (zh) * | 2020-07-30 | 2020-12-08 | 南方电网科学研究院有限责任公司 | 基于簇类的风电出力波动区间预测方法、装置及存储介质 |
CN112052981B (zh) * | 2020-07-30 | 2023-08-22 | 南方电网科学研究院有限责任公司 | 基于簇类的风电出力波动区间预测方法、装置及存储介质 |
CN111861278A (zh) * | 2020-08-04 | 2020-10-30 | 长沙理工大学 | 一种电力系统的用电负荷峰谷时段划分方法及系统 |
CN111861278B (zh) * | 2020-08-04 | 2023-06-09 | 国家电网有限公司西北分部 | 一种电力系统的用电负荷峰谷时段划分方法及系统 |
CN112766733A (zh) * | 2021-01-21 | 2021-05-07 | 山东大学 | 利用改进的K-means算法加速优化调度算法收敛的方法及系统 |
CN113792778A (zh) * | 2021-09-07 | 2021-12-14 | 国网经济技术研究院有限公司 | 一种风电场分区及其出力相关系数计算方法和系统 |
CN113792778B (zh) * | 2021-09-07 | 2023-10-10 | 国网经济技术研究院有限公司 | 一种风电场分区及其出力相关系数计算方法和系统 |
CN117318126A (zh) * | 2023-11-30 | 2023-12-29 | 山东华天电气有限公司 | 一种微网电能质量优化与储能一体化装置及控制方法 |
CN117318126B (zh) * | 2023-11-30 | 2024-02-27 | 山东华天电气有限公司 | 一种微网电能质量优化与储能一体化装置及控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108376262B (zh) | 2021-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108376262A (zh) | 一种风电出力典型特性的分析模型构建方法 | |
CN110428168A (zh) | 一种计及储能的多能源配电系统协调调度综合评价方法 | |
CN107133652A (zh) | 基于K‑means聚类算法的用电客户价值评估方法及系统 | |
CN108805213B (zh) | 计及小波熵降维的电力负荷曲线双层谱聚类方法 | |
CN103473621A (zh) | 风电场短期功率预测方法 | |
CN103489046A (zh) | 风电场短期功率预测方法 | |
CN112186761B (zh) | 一种基于概率分布的风电功率场景生成方法及系统 | |
CN107679687A (zh) | 一种光伏出力建模方法以及发电系统可靠性评估方法 | |
CN115173465A (zh) | 一种基于Copula理论的风光水火储一体化耦合机理分析方法 | |
CN108765194A (zh) | 一种有效的居民用电行为分析系统 | |
CN107730399B (zh) | 基于风力发电特性曲线的理论线损评估方法 | |
CN110380444A (zh) | 一种基于变结构Copula的多场景下分散式风电有序接入电网的容量规划方法 | |
CN111428745A (zh) | 一种基于聚类分析的低压用户用电特征提取方法 | |
CN104915727B (zh) | 多维度同构异质bp神经网络光功率超短期预测方法 | |
CN108898273A (zh) | 一种基于形态分析的用户侧负荷特征聚类评价方法 | |
CN112365164B (zh) | 基于改进密度峰值快速搜索聚类算法的中大型能源用户用能特性画像方法 | |
CN109586309B (zh) | 基于大数据自由熵理论和场景匹配的配电网无功优化方法 | |
CN117613993A (zh) | 一种公路微电网风光氢蓄储系统的日内优化调度方法 | |
CN115859775B (zh) | 一种面向光储虚拟电厂划分的典型日场景提取方法及装置 | |
CN107832928B (zh) | 基于风力发电特性曲线的等效综合线损评价方法 | |
CN109149644B (zh) | 一种基于大数据分析的光储一体在线策略匹配和协同优化方法 | |
Fan et al. | Prediction and Analysis of Power User Energy Consumption Based on Demand Side Management | |
CN114186631A (zh) | 一种基于非侵入式智能终端的负荷识别方法 | |
CN113505943A (zh) | 一种电网短期负荷的预测方法、系统、设备和介质 | |
CN112949938A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210810 |
|
CF01 | Termination of patent right due to non-payment of annual fee |