CN106485362A - 一种基于高维风电预测误差模型及降维技术的发电调度方法 - Google Patents

一种基于高维风电预测误差模型及降维技术的发电调度方法 Download PDF

Info

Publication number
CN106485362A
CN106485362A CN201610907866.XA CN201610907866A CN106485362A CN 106485362 A CN106485362 A CN 106485362A CN 201610907866 A CN201610907866 A CN 201610907866A CN 106485362 A CN106485362 A CN 106485362A
Authority
CN
China
Prior art keywords
energy turbine
wind energy
turbine set
function
probability density
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
Application number
CN201610907866.XA
Other languages
English (en)
Other versions
CN106485362B (zh
Inventor
卫鹏
刘建坤
周前
汪成根
徐青山
黄煜
陈静
陈哲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Southeast University
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Southeast University
Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Southeast University, Electric Power Research Institute of State Grid Jiangsu Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610907866.XA priority Critical patent/CN106485362B/zh
Publication of CN106485362A publication Critical patent/CN106485362A/zh
Application granted granted Critical
Publication of CN106485362B publication Critical patent/CN106485362B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • YGENERAL 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS 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/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Supply And Distribution Of Alternating Current (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种基于高维风电预测误差模型及降维技术的发电调度方法,包括:获取多个风电场一年内每小时的历史出力数据及相应的点预测数据;利用混合偏态模型对每个单一风电场的实际出力和预测出力的累积分布函数进行建模;利用每个风电场各自的CDF将实际出力和预测值转换为0‑1区间分布的数据点;通过匹配前一步骤中得到的所有数据点,找到最优的Copula函数并进行参数估计;建立多个风电场预测误差的高维条件概率模型,并通过边缘变换得到降维后的边缘条件概率模型;根据风电场预测误差的边缘条件概率模型,计算发电机组的日前调度计划及旋转备用容量。本发明相较于常用的高斯分布、β分布有着更高的精度,且可以考虑多个风电场之间相关性的影响。

Description

一种基于高维风电预测误差模型及降维技术的发电调度方法
技术领域
本发明涉及电力系统的新能源发电领域,具体涉及一种基于高维风电预测误差模型及降维技术的发电调度方法。
背景技术
随着能源消耗的日益增加,能源供给也持续紧张。以风电为代表的可再生能源得到了很大的发展。但是风电出力具有较强随机波动性、较差的功率调节能力等特点,受气象环境、风场布局等影响性较大,因此风电集群规模化接入对给电力系统的安全运行与调度规划带来巨大挑战。此外目前风电建设规划存在着“重发、轻供、不管用”,风电发展超前于相应区域电网规划,两者规划发展不协调,这些使得我国风电存在着严重的弃风现象,深究导致风电并网困难产生弃风现象的主要原因是风电出力具有随机波动性。因此为解决大规模风电并网带来的一系列问题,风电功率预测技术应运而生。虽然风电功率具有极强的波动性和随机性,但是随着时间尺度的变小,风电功率的规律性也逐渐增强,使得超短期风电功率预测成为可能。超短期风电功率预测将有助于电力系统调度部门进一步了解即将入网的风电功率,为小时级别的发电运行调度提供依据。
风电功率预测误差分布主要有正态分布、β分布、柯西分布以及拉普拉斯分布等几种传统的经验分布模型。但是β分布其求解过程非常复杂,某些预测功率区间上会出现误差概率密度无穷大的异常;柯西分布和拉普拉斯分布都是对称性分布,形状不够灵活,适用范围较固定。针对短期风电功率预测误差分布的多样性,尚未有一种分布模型,能够适应更多的误差分布特性,既能描述对称分布和单峰分布,又能描述不对称分布和多峰分布。
发明内容
发明目的:本发明旨在提供一种有更高的精度,且考虑多个风电场之间相关性的影响的基于高维风电预测误差模型及降维技术的发电调度方法。
技术方案:一种基于高维风电预测误差模型及降维技术的发电调度方法,包括如下步骤:
S1:通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型;
S2:根据步骤S1得到的风电场预测误差的边缘条件概率密度函数模型,计算发电机组的日前调度计划及旋转备用容量。
进一步的,所述步骤S1通过风电场在预设时间内的数据建立多个风电场预测误差的条件概率密度函数模型,并边缘变换得到降维后的边缘条件概率密度函数模型具体为:
S1-1:获取多个风电场在预设时间内的历史出力数据及相应的点预测数据;
S1-2:利用混合偏态模型对步骤S1-1中获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模;
S1-3:利用每个风电场各自的累积分布函数将实际出力和预测值转换为0-1区间分布的数据点;
S1-4:匹配步骤S1-3中得到的所有数据点,并找到最优的Copula函数并进行参数估计;
S1-5:利用步骤S4中找到的最优Copula函数建立多个风电场预测误差的条件概率密度函数模型,并通过边缘变换得到降维后的边缘条件概率。
更进一步的,所述步骤S1-1中获取多个风电场一年内每小时的历史出力数据及相应的点预测数据具体为:将获取的多个风电场实际出力的历史数据记为pi,t,其中,i∈[1,N]表示风电场的编号,t∈[1,T]表示历史数据的时间范围;在一年内按小时为最小单位进行记录;与实际出力相对应的点预测值记为其由各风电场的NWP数据及神经网络方法生成。
进一步的,所述步骤S1-2中利用混合偏态模型对步骤S1-1中获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模具体为:利用的混合偏态模型包括混合偏态分布,混合偏态分布由正态分布线性叠加得到;若一个随机变量X服从位置参数为μ∈R、尺度参数为σ2∈(0,∞)和偏度参数为λ∈R的偏正态分布,则它的概率密度函数是:
式中分别为标准正态分布的密度函数和分布函数,分别表示为:
简记为X~SN(μ,σ2,λ);当μ=0,σ2=1时,X为标准偏正态随机变量,记为X~SN(λ);
混合偏态分布是多个偏态分布的线性组合;具有K个成份的基于偏态分布的混合分布模型的概率密度函数为:
式中αk为不同偏态分布在混合偏态分布中的权重系数,αk>0且α1+…+αk+…+αK=1。
进一步的,所述步骤S1-3具体为:将各风电场实际出力pi,t及相应的点预测值通过各自的累积分布函数 转换成[0,1]区间内分布的数据点ui,t
其中,FX(x)~u(0,1),FY(y)~u(0,1)。
进一步的,所述步骤S1-4中找到最优的Copula函数采用基于经验Copula函数和理论Copula函数的最短欧式距离法;具体为:
设由步骤S3得到的N个风电场在[0,1]区间的数据向量U=(u1,u2,…,uN),其中ui(i=1,2,…,N)的概率密度函数fi(ui)及累积分布函数Fi(ui)已知,其经验Copula函数可由下式求得:
式中,为示性函数,括号中的条件满足时,I=1,反之则I=0;{(u1,j,u2,j,…,uN,j),j=1,2,…,T}是容量为T的时间序列;为顺序统计量,1≤i1,…,iN≤T;
根据式(5)得到经验Copula函数后,由下式计算经验Copula函数与理论Copula函数之间的欧氏距离:
其中为相应的Copula分布函数;
Copula函数与理论Copula函数之间的欧式距离表征了Copula函数与随机变量实际分布的拟合程度;欧式距离越小,说明该Copula模型拟合程度越好,反之则越差。
更进一步的,所述步骤S1-5中建立多个风电场预测误差的条件概率密度函数模型,并通过边缘变换得到降维后的边缘条件概率密度函数模型具体为:风力发电的实际出力与它的预测值之间的联合分布以Copula函数来模拟,设x为实际出力值,y为出力预测值,则它们的联合概率密度函数fXY可以写为:
式中fX(x)和fY(y)分别为x和y的边缘概率密度函数,c为Copula函数;当给定点预测值则实际值的条件概率密度函数fX|Y表示为:
设预测误差为e=x-y,则预测误差的条件概率密度函数fE|Y为:
e的取值范围为
将式(9)推广到模拟多个风电场的情况;假设是风电场1到N的点预测值,e1,e2,…eN表示预测误差,e1,e2,…eN的联合条件概率密度函数在给定点预测值为时,
式中表示所有风电场的点预测值的联合概率密度函数,表示它们的实际出力以及预测值的联合概率密度函数;cN和c2N是分别与相对应的Copula密度函数;是每个风电场预测出力的边缘分布函数;
通过边缘变换得到高维风电场预测误差模型降维后的每个单一风电场的边缘条件概率密度函数
式中ci,N是Copula函数。
更进一步的,所述步骤S1-6计算发电机组的日前调度计划及旋转备用容量具体为:
通过边缘变换得到高维风电场预测误差模型降维后的每个单一风电场的边缘条件概率密度函数模型,对单一风电场的边缘条件概率密度函数求均值记为旋转备用容量Z,发电机组的日前调度计划W,
W=Z+W0
其中W0为除风电场外机组的发电计划容量。
本发明还提供了一种基于高维风电预测误差模型及降维技术的发电调度系统,包括:
用于通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型的单元;
用于根据风电场预测误差的边缘条件概率密度函数模型计算发电机组日前调度计划及旋转备用容量的单元。
进一步的,所述用于通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型的单元包括;
用于获取多个风电场在预设时间内的历史出力数据及相应点预测数据的单元;
用于利用混合偏态模型对获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模的单元;
用于将每个风电场各自的累积分布函数将实际出力和预测值转换为0-1区间分布的数据点的单元;
用于匹配所有数据点并找到最优的Copula函数进行参数估计的单元;
用于利用最优Copula函数建立多个风电场预测误差的条件概率密度函数模型并通过边缘变换得到降维后的边缘条件概率的单元。
有益效果:相对于现有技术,本发明采用通用分布来拟合不同风电预测区间上实际风电功率的概率密度函数PDF和累积分布函数CDF,能较好地适用任意时间尺度和幅值的风电预测误差建模,其模型拟合精度高于常规的高斯分布和β分布,且本发明不仅考虑了每个风电场预测误差与点预测值间的关系,而且考虑了多个风电场之间相关性的影响。除此之外,本发明还采用了基于混合偏态分布的预测误差分布模型,该模型是偏态分布与混合分布的结合,相对于传统误差分布模型,更加灵活可变。正态分布考虑偏态后,增加了偏态系数,可以解决风电出力条件预测误差分布中的多峰现象和重尾现象。
附图说明
图1为本发明的流程示意图;
图2为本发明提供的双峰分布的混合偏态分布图;
图3为本发明提供的风电场1的实际值-预测值散点图(CDF);
图4为本发明提供的所有风电场实际值-预测值的Gaussia copula相关系数图;
图5为本发明提供的点预测值分别为4.99MW、9.81MW时风电场1与风电场2实际出力值的联合条件分布估计图;
图6为本发明提供的点预测值分别为5.30MW、14.12MW时风电场1的条件分布估计图。
具体实施方式
下面将结合附图,对本发明的实施案例进行详细的描述;
建在同一地区的多个风电场由于地理位置相对集中,所处的气象条件和环境也大概相同,所以这些风电场间的风速具有一定的相关性,而与风速密切相关的风电的输出功率也具有一定的相关性。本发明基于考虑到多风电场出力之间的相关性,提出一种基于高维风电预测误差模型及降维技术的发电调度方法。
如图1所示,本发明所述的一种基于高维风电预测误差模型及降维技术的发电调度方法,包括如下步骤:
1):获取多个风电场一年内每小时的历史出力数据及相应的点预测数据;具体为:将获取的江苏多个风电场实际出力的历史数据记为pi,t,这里i∈[1,N]表示风电场的编号,t∈[1,T]表示历史数据的时间范围,在一年内按小时为最小单位进行记录,共8760个数据点。与实际出力相对应的点预测值记为其由各风电场的NWP数据(numerical weatherprediection,数值天气预报)及神经网络方法生成。
2):利用混合偏态模型对每个单一风电场的实际出力和预测出力的累积分布函数(Cumulative Distribution Function,CDF)进行建模;具体包括:混合偏态模型中的混合偏态分布是由正态分布线性叠加得到的。若一个随机变量X服从位置参数为μ∈R、尺度参数为σ2∈(0,∞)和偏度参数为λ∈R的偏正态分布,则它的概率密度函数是:
式中分别为标准正态分布的密度函数和分布函数,分别表示为:
可简记为X~SN(μ,σ2,λ)。当μ=0,σ2=1时,X为标准偏正态随机变量,记为X~SN(λ)。
混合偏态模型中的混合偏态分布是多个偏态分布的线性组合。具有K个成份的基于偏态分布的混合分布模型的概率密度函数为:
式中αk为不同偏态分布在混合偏态分布中的权重系数,αk>0且α1+…+αk+…+αK=1。
混合偏态分布可以描述非常丰富的分布情况,能够更好地描述短期风电功率误差分布的有偏性、重尾性和多峰性。双峰分布的混合偏态分布如图2所示。对概率密度函数f在(-∞,+∞)区间内积分,得到单一风电场实际出力和预测值的累积分布函数(CDF)FX,FY
3):利用每个风电场各自的累积分布函数CDF将实际出力和预测值转换为0-1区间分布的数据点;具体为:将各风电场实际出力pi,t及相应的点预测值通过各自的累积分布函数 转换成[0,1]区间内分布的数据点ui,t
其中,FX(x)~u(0,1),FY(y)~u(0,1)。图3为标准化转换后的风电场1的实际值-预测值散点图。
4):通过匹配步骤3)中得到的所有数据点,找到最优的Copula函数并进行参数估计;其中采用基于经验Copula函数(empirical Copula,EMC)和理论Copula函数(theoretical Copula Function,THC)的最短欧式距离法实现最优Copula函数(把多维随机变量的联合分布用一维边际分布连接起来的函数)的选取。测试表明,在Gumbel copula、Gaussian copula、t copula中,Gumbel copula有着最优的性能。考虑到在多变量系统中Gaussian copula比Gumbel copula更加灵活,本发明采用Gaussian copula(又称正态Copula,是多元正态分布的连接函数)进行研究。具体为:
设由步骤3)得到的N个风电场在[0,1]区间的数据向量U=(u1,u2,…,uN),其中ui(i=1,2,…,N)的概率密度函数fi(ui)及累积分布函数Fi(ui)已知,其EMC可由下式求得:
式中,为示性函数,括号中的条件满足时,I=1,反之则I=0;{(u1,j,u2,j,…,uN,j),j=1,2,…,T}是容量为T的时间序列;为顺序统计量,1≤i1,…,iN≤T。
根据式(5)得到EMC后,可由下式计算EMC与THC之间的欧氏距离:
其中为相应的Copula分布函数。
EMC与THC之间的欧式距离表征了Copula函数与随机变量实际分布的拟合程度。欧式距离越小,说明该Copula模型拟合程度越好,反之则越差。因此选取最优的Copula函数作为描述各风电场实际出力与预测值之间的连接函数。
5):利用步骤(4)中找到的最优Copula函数建立多个风电场预测误差的条件概率密度函数模型,并通过边缘变换得到降维后的边缘条件概率密度函数模型;
由于风速的随机性,风力发电的出力预测值与实际值之间并没有一个确定性的关系,因此它们看作一对具有相关性的随机变量,即其中一个变量的取值会对另一个变量的概率分布产生影响。风力发电的实际出力与它的预测值之间的联合分布可以Copula函数来模拟,设x为实际出力值,y为出力预测值,则它们的联合PDF(概率密度函数)fXY可以写为:
式中fX(x)和fY(y)分别为x和y的边缘PDF,c为Copula函数。当给定点预测值则实际值的CPDF(条件概率密度函数)fX|Y可以表示为:
设预测误差为e=x-y,则预测误差的条件概率密度函数fE|Y为:
e的取值范围为
式(9)提供了一种利用点预测值将风力发电的条件预测误差量化的方法。它表明预测的分布包括两个方面:作为基础函数的实际出力的概率分布,以及作为可变乘子的Copula密度函数。因此,对条件预测误差的计算可以转变为对这两个函数的计算。
式(9)可以被推广到模拟多个风电场的情况。假设是风电场1到N的点预测值,e1,e2,…eN表示预测误差,e1,e2,…eN的联合在给定点预测值为时,如式(10)所示,
式中表示所有风电场的点预测值的联合PDF,表示它们的实际出力以及预测值的联合PDF。cN和c2N是分别与相对应的Copula密度函数。是每个风电场预测出力的CDF(边缘分布函数)。
通过边缘变换可以得到高维风电场预测误差模型降维后的每个单一风电场的边缘
式中ci,N是Copula函数,它只跟单一风电场实际出力与所有风电场出力预测值之间的相关性有关。
式(10)和式(11)不仅考虑了每个风电场预测误差与点预测值间的关系,而且考虑了所有风电场之间的相关性。式(11)中的各个预测误差的边缘CPDF之间也有联系。
本实施例中:建立多个风电场预测误差的边缘条件概率密度函数模型,并通过边缘变换得到降维后的边缘条件概率密度函数模型。图4展示了相关系数的确立结果,它清楚的显示了不同风电场的出力间、不同风电场预测值间的空间相关性,以及二者之间的联系。图5计算了当点预测值分别为4.99MW(总量的24.95%)与9.81MW(总量的49.3%)时,风电场1和2的关联条件预测误差。上图清楚地表明两个风电场条件预测误差之间的关联性。这种相关性理解为:由于相邻风电场总是处在相同的天气与环境中。而这些,预测工具是不能测出的。预测误差具有关联性。所以,当风电场发生一个较大的预测误差,相邻的其他风电场也很有可能发生误差。图6通过将估计的CPDF与风电场1的实际分布相比较判断模型的性能。左侧的直方图展示了点预测值为5.30MW下估计与实际的CPDF数据,而右侧展示了点预测值为14.12MW下的情况。这里通过低维copula函数和高维copula函数对预测值误差的条件概率密度函数模型进行了估计。
6):根据风电场预测误差的边缘条件概率密度函数模型,计算发电机组的日前调度计划及旋转备用容量;
利用风电场预测误差的边缘条件概率密度函数模型来计算日前调度计划,它能考虑每种风电出力场景并分配相应的备用容量。为了清楚地显示之前建立的条件预测误差模型的价值,这里只考虑风电的不确定性。风电场预测误差的边缘条件概率密度函数模型包含即时和观望决策的变量及约束,风电出力的日前点预测值可以形成即时决策的约束条件。而风电实际出力的场景,通过对建立的条件预测误差模型进行采样,形成对应的观望决策约束条件。采样过程包含三个步骤:A.根据点预测值,通过式(11)建立降维后风电预测误差的条件概率密度函数模型;B.由条件概率密度函数模型生成一天内风电实际出力的大规模场景;C.通过K均值聚类技术减少场景数量以提高计算效率。具体为:
通过边缘变换得到高维风电场预测误差模型降维后的每个单一风电场的边缘条件概率密度函数模型,对单一风电场的边缘条件概率密度函数求均值记为旋转备用容量Z,发电机组的日前调度计划W,
W=Z+W0
其中W0为除风电场外机组的发电计划容量。
本发明所述的一种基于高维风电预测误差模型及降维技术的发电调度系统,包括:
用于通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型的单元;
用于根据风电场预测误差的边缘条件概率密度函数模型计算发电机组日前调度计划及旋转备用容量的单元。
其中,用于通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型的单元还包括:
用于获取多个风电场在预设时间内的历史出力数据及相应点预测数据的单元;
用于利用混合偏态模型对获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模的单元;
用于将每个风电场各自的累积分布函数将实际出力和预测值转换为0-1区间分布的数据点的单元;
用于匹配所有数据点并找到最优的Copula函数进行参数估计的单元;
用于利用最优Copula函数建立多个风电场预测误差的条件概率密度函数模型并通过边缘变换得到降维后的边缘条件概率的单元。
对应上述方法实施例的描述,该发电调度系统的各组成单元能够逐步实现本发明的基于高维风电预测误差模型及降维技术的发电调度方法,达到相同的技术效果,对该系统的工作过程不再赘述。本发明先用混合偏态分布构建单个风电场的分布模型,混合偏态分布相对于正态分布考虑偏态后,增加了偏态系数,用来描述不对称分布,可以解决短期风电出力条件预测误差分布中的多峰现象和重尾现象。之后利用Copula函数对风电场实际出力与预测值及多个风电场之间进行相关性建模,建立多个风电场的高维预测误差模型并通过边缘变换进行降维,降维后的风电预测误差模型相较于常用的高斯分布、β分布有着更高的精度,不仅考虑了每个风电场预测误差与点预测值间的关系,而且可以考虑多个风电场之间相关性的影响。利用Copula理论可以将对变量联合分布的建模转变为对变量边缘分布和相依结构分别建模,极大地方便了对非独立随机变量的研究。本发明对于风电预测误差高维与低维建模和多个风电场之间的相关性的考虑具有较好的准确性和有效性。

Claims (10)

1.一种基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,包括如下步骤:
S1:通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型;
S2:根据步骤S1得到的风电场预测误差的边缘条件概率密度函数模型,计算发电机组的日前调度计划及旋转备用容量。
2.根据权利要求1所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型具体为:
S1-1:获取多个风电场在预设时间内的历史出力数据及相应的点预测数据;
S1-2:利用混合偏态模型对步骤S1-1中获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模;
S1-3:利用每个风电场各自的累积分布函数将实际出力和预测值转换为0-1区间分布的数据点;
S1-4:匹配步骤S1-3中得到的所有数据点,并找到最优的Copula函数并进行参数估计;
S1-5:利用步骤S1-4中找到的最优Copula函数建立多个风电场预测误差的条件概率密度函数模型,并通过边缘变换得到降维后的边缘条件概率。
3.根据权利要求2所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1-1中获取多个风电场一年内每小时的历史出力数据及相应的点预测数据具体为:将获取的多个风电场实际出力的历史数据记为pi,t,其中,i∈[1,N]表示风电场的编号,t∈[1,T]表示历史数据的时间范围;在一年内按小时为最小单位进行记录;与实际出力相对应的点预测值记为其由各风电场的NWP数据及神经网络方法生成。
4.根据权利要求2所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1-2中利用混合偏态模型对步骤S1-1中获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模具体为:利用的混合偏态模型包括混合偏态分布,混合偏态分布由正态分布线性叠加得到;若一个随机变量X服从位置参数为μ∈R、尺度参数为σ2∈(0,∞)和偏度参数为λ∈R的偏正态分布,则它的概率密度函数是:
f S N ( x ; μ , σ 2 , λ ) = 2 σ φ ( x - μ σ ) Φ ( λ x - μ σ ) - - - ( 1 )
式中分别为标准正态分布的密度函数和分布函数,分别表示为:
φ ( x ) = 1 2 π exp ( - x 2 2 ) Φ ( x ) = 1 2 π ∫ - ∞ x exp ( - t 2 2 ) d t - - - ( 2 )
简记为X~SN(μ,σ2,λ);当μ=0,σ2=1时,X为标准偏正态随机变量,记为X~SN(λ);
混合偏态分布是多个偏态分布的线性组合;具有K个成份的基于偏态分布的混合分布模型的概率密度函数为:
f M S N ( x ; μ 1 , ... , μ K , σ 1 2 , ... , σ K 2 , λ 1 , ... , λ K ) = Σ k = 1 K α k f S N ( x ; μ k , σ k 2 , λ k ) - - - ( 3 )
式中αk为不同偏态分布在混合偏态分布中的权重系数,αk>0且α1+…+αk+…+αK=1。
5.根据权利要求2所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1-3具体为:将各风电场实际出力pi,t及相应的点预测值通过各自的累积分布函数转换成[0,1]区间内分布的数据点ui,t
u i , t = F X i ( p i , t ) , u ^ i , t = F Y i ( p i , t ) i ∈ [ 1 , N ] , t ∈ [ 1 , T ] - - - ( 4 )
其中,FX(x)~u(0,1),FY(y)~u(0,1)。
6.根据权利要求1所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1-4中找到最优的Copula函数采用基于经验Copula函数和理论Copula函数的最短欧式距离法;具体为:
设由步骤S3得到的N个风电场在[0,1]区间的数据向量U=(u1,u2,…,uN),其中ui(i=1,2,…,N)的概率密度函数fi(ui)及累积分布函数Fi(ui)已知,其经验Copula函数可由下式求得:
C e ( i 1 n , i 2 n , ... , i N n ) = 1 n Σ j = 1 T I ( u 1 , j ≤ u 1 i 1 , u 2 , j ≤ u 2 i 2 , ... , u N , j ≤ u 1 i N ) - - - ( 5 )
式中,为示性函数,括号中的条件满足时,I=1,反之则I=0;{(u1,j,u2,j,…,uN,j),j=1,2,…,T}是容量为T的时间序列;为顺序统计量,1≤i1,…,iN≤T;
根据式(5)得到经验Copula函数后,由下式计算经验Copula函数与理论Copula函数之间的欧氏距离:
d ( C , C e ) = { Σ i 2 = 1 n ... Σ i N = 1 n [ C ( i 1 n , ... , i N n ) - C e ( i 1 n , ... , i N n ) ] 2 } - - - ( 6 )
其中为相应的Copula分布函数;
Copula函数与理论Copula函数之间的欧式距离表征了Copula函数与随机变量实际分布的拟合程度;欧式距离越小,说明该Copula模型拟合程度越好,反之则越差。
7.根据权利要求2所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1-5中建立多个风电场预测误差的条件概率密度函数模型,并通过边缘变换得到降维后的边缘条件概率密度函数模型具体为:风力发电的实际出力与它的预测值之间的联合分布以Copula函数来模拟,设x为实际出力值,y为出力预测值,则它们的联合概率密度函数fXY可以写为:
f X Y ( x , y ) = F X Y ( x , y ) ∂ x ∂ y = c ( F X ( x ) , F Y ( y ) ) f X ( x ) f Y ( y ) - - - ( 7 )
式中fX(x)和fY(y)分别为x和y的边缘概率密度函数,c为Copula函数;当给定点预测值则实际值的条件概率密度函数fX|Y表示为:
f X | Y ( x | y = p ^ ) = f X Y ( x , p ^ ) f Y ( p ^ ) = c ( F X ( x ) , F Y ( p ^ ) ) f X ( x ) - - - ( 8 )
设预测误差为e=x-y,则预测误差的条件概率密度函数fE|Y为:
f E | Y ( e | y = p ^ ) = f X Y ( x , p ^ ) f Y ( p ^ ) = c ( F X ( e + p ^ ) , F Y ( p ^ ) ) f X ( e + p ^ ) - - - ( 9 )
e的取值范围为
将式(9)推广到模拟多个风电场的情况;假设是风电场1到N的点预测值,e1,e2,…eN表示预测误差,e1,e2,…eN的联合条件概率密度函数在给定点预测值为时,
f E 1 , ... , E N | Y 1 , ... , Y N ( e 1 , ... e N | y 1 = p ^ 1 , ... , y N = p ^ N ) = f X 1 , ... , X N , Y 1 , ... , Y N ( e 1 + P 1 , ... e N + P N , P 1 , ... , P N ) f Y 1 , ... , Y N ( P 1 , ... , P N ) = c 2 N ( F X 1 ( e 1 + P 1 ) , F X 2 ( e N + P N ) , F Y 1 ( P 1 ) , ... , F Y N ( P N ) ) c N ( F N 1 ( P 1 ) , ... , F Y N ( P N ) ) Π i = 1 N f X i ( e i + P i ) - - - ( 10 )
式中表示所有风电场的点预测值的联合概率密度函数,表示它们的实际出力以及预测值的联合概率密度函数;cN和c2N是分别与相对应的Copula密度函数;是每个风电场预测出力的边缘分布函数;
通过边缘变换得到高维风电场预测误差模型降维后的每个单一风电场的边缘条件概率密度函数
f E i | Y 1 , ... , Y N ( e i | y 1 = P 1 , ... , y N = P N ) = c i , N ( F X i ( e i + P i ) , F Y 1 ( P 1 ) , ... , F Y N ( P N ) ) c N ( F Y 1 ( P 1 ) , ... , F Y N ( P N ) ) f X i ( e i + P i ) , i ∈ [ 1 , N ] - - - ( 11 )
式中ci,N是Copula函数。
8.根据权利要求2所述的基于高维风电预测误差模型及降维技术的发电调度方法,其特征在于,所述步骤S1-6计算发电机组的日前调度计划及旋转备用容量具体为:
通过边缘变换得到高维风电场预测误差模型降维后的每个单一风电场的边缘条件概率密度函数模型,对单一风电场的边缘条件概率密度函数求均值记为旋转备用容量Z,发电机组的日前调度计划W,
Z = Σ i = 1 N Ω E i | Y 1 , ... , Y N
W=Z+W0
其中W0为除风电场外机组的发电计划容量。
9.一种基于高维风电预测误差模型及降维技术的发电调度系统,其特征在于,包括:
用于通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型的单元;
用于根据风电场预测误差的边缘条件概率密度函数模型计算发电机组日前调度计划及旋转备用容量的单元。
10.根据权利要求9所述的基于高维风电预测误差模型及降维技术的发电调度系统,其特征在于,所述用于通过风电场在预设时间内的数据建立多个风电场预测误差的边缘条件概率密度函数模型的单元包括;
用于获取多个风电场在预设时间内的历史出力数据及相应点预测数据的单元;
用于利用混合偏态模型对获取的每个单一风电场的实际出力和预测出力的累积分布函数进行建模的单元;
用于将每个风电场各自的累积分布函数将实际出力和预测值转换为0-1区间分布的数据点的单元;
用于匹配所有数据点并找到最优的Copula函数进行参数估计的单元;
用于利用最优Copula函数建立多个风电场预测误差的条件概率密度函数模型并通过边缘变换得到降维后的边缘条件概率的单元。
CN201610907866.XA 2016-10-18 2016-10-18 一种基于高维风电预测误差模型及降维技术的发电调度方法 Active CN106485362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610907866.XA CN106485362B (zh) 2016-10-18 2016-10-18 一种基于高维风电预测误差模型及降维技术的发电调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610907866.XA CN106485362B (zh) 2016-10-18 2016-10-18 一种基于高维风电预测误差模型及降维技术的发电调度方法

Publications (2)

Publication Number Publication Date
CN106485362A true CN106485362A (zh) 2017-03-08
CN106485362B CN106485362B (zh) 2019-10-18

Family

ID=58270143

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610907866.XA Active CN106485362B (zh) 2016-10-18 2016-10-18 一种基于高维风电预测误差模型及降维技术的发电调度方法

Country Status (1)

Country Link
CN (1) CN106485362B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292439A (zh) * 2017-06-23 2017-10-24 广东工业大学 一种基于时变混合Copula函数的短期风速预测的方法与装置
CN108985515A (zh) * 2018-07-24 2018-12-11 国网河南省电力公司电力科学研究院 一种基于独立循环神经网络的新能源出力预测方法及系统
CN109755959A (zh) * 2018-12-11 2019-05-14 清华大学 基于风/光出力柯西分布的火电机组动态实时调度方法
CN109816149A (zh) * 2018-12-28 2019-05-28 中国电力科学研究院有限公司 一种风电场随机出力场景生成方法及装置
CN110378504A (zh) * 2019-04-12 2019-10-25 东南大学 一种基于高维Copula技术的光伏发电爬坡事件概率预测方法
CN111274542A (zh) * 2020-01-10 2020-06-12 国网江苏省电力有限公司扬州供电分公司 基于藤copula和混合偏移正态分布的置信容量评估方法及装置
CN111563623A (zh) * 2020-04-30 2020-08-21 国网山东省电力公司威海供电公司 一种风电系统规划的典型场景提取方法及系统
CN111598612A (zh) * 2020-04-28 2020-08-28 西安理工大学 一种分时电价制定方法
CN112232886A (zh) * 2020-10-30 2021-01-15 南方电网能源发展研究院有限责任公司 一种电价概率预测方法、系统、计算机设备和存储介质
CN112380714A (zh) * 2020-11-19 2021-02-19 国网宁夏电力有限公司电力科学研究院 基于统计规律的新能源出力概率区间建模方法
CN112417768A (zh) * 2020-12-07 2021-02-26 天津大学 一种基于藤结构Pair-Copula的风电相关性条件采样方法
CN112581312A (zh) * 2020-12-21 2021-03-30 国网陕西省电力公司电力科学研究院 风电功率预测误差分布分析方法、装置、计算机设备和可读存储介质
CN112685915A (zh) * 2021-01-18 2021-04-20 重庆大学 一种风电出力条件概率分布建模方法
CN112711864A (zh) * 2021-01-18 2021-04-27 国网浙江省电力有限公司电力科学研究院 一种配网电缆质量指标相关性模型构建及数据拓展方法
CN113937769A (zh) * 2021-11-05 2022-01-14 国网甘肃省电力公司 考虑风能出力不确定性的机组组合确定方法及装置
CN112711864B (zh) * 2021-01-18 2024-06-07 国网浙江省电力有限公司电力科学研究院 一种配网电缆质量指标相关性模型构建及数据拓展方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140172329A1 (en) * 2012-12-17 2014-06-19 Junshan Zhang System and method for wind generation forecasting
CN105846425A (zh) * 2016-04-08 2016-08-10 江苏省电力试验研究院有限公司 一种基于通用风电预测误差模型的经济调度方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140172329A1 (en) * 2012-12-17 2014-06-19 Junshan Zhang System and method for wind generation forecasting
CN105846425A (zh) * 2016-04-08 2016-08-10 江苏省电力试验研究院有限公司 一种基于通用风电预测误差模型的经济调度方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘燕华等: "短期风电功率预测误差的混合偏态分布模型", 《中国电机工程学报》 *
叶希等: "基于集群风电虚拟机组的源端系统分层协调日前计划", 《电力系统自动化》 *
李玉敦等: "基于Copula函数的多维时序风速相依模型及其在可靠性评估中的应用", 《电网技术》 *
赵俊华等: "计及电动汽车和风电出力不确定性的随机经济调度", 《电力系统自动化》 *

Cited By (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107292439A (zh) * 2017-06-23 2017-10-24 广东工业大学 一种基于时变混合Copula函数的短期风速预测的方法与装置
CN108985515A (zh) * 2018-07-24 2018-12-11 国网河南省电力公司电力科学研究院 一种基于独立循环神经网络的新能源出力预测方法及系统
CN108985515B (zh) * 2018-07-24 2021-11-26 国网河南省电力公司电力科学研究院 一种基于独立循环神经网络的新能源出力预测方法及系统
CN109755959A (zh) * 2018-12-11 2019-05-14 清华大学 基于风/光出力柯西分布的火电机组动态实时调度方法
CN109816149A (zh) * 2018-12-28 2019-05-28 中国电力科学研究院有限公司 一种风电场随机出力场景生成方法及装置
CN110378504A (zh) * 2019-04-12 2019-10-25 东南大学 一种基于高维Copula技术的光伏发电爬坡事件概率预测方法
CN111274542A (zh) * 2020-01-10 2020-06-12 国网江苏省电力有限公司扬州供电分公司 基于藤copula和混合偏移正态分布的置信容量评估方法及装置
CN111598612A (zh) * 2020-04-28 2020-08-28 西安理工大学 一种分时电价制定方法
CN111598612B (zh) * 2020-04-28 2023-04-18 西安理工大学 一种分时电价制定方法
CN111563623A (zh) * 2020-04-30 2020-08-21 国网山东省电力公司威海供电公司 一种风电系统规划的典型场景提取方法及系统
CN111563623B (zh) * 2020-04-30 2022-05-10 国网山东省电力公司威海供电公司 一种风电系统规划的典型场景提取方法及系统
CN112232886A (zh) * 2020-10-30 2021-01-15 南方电网能源发展研究院有限责任公司 一种电价概率预测方法、系统、计算机设备和存储介质
CN112380714A (zh) * 2020-11-19 2021-02-19 国网宁夏电力有限公司电力科学研究院 基于统计规律的新能源出力概率区间建模方法
CN112417768A (zh) * 2020-12-07 2021-02-26 天津大学 一种基于藤结构Pair-Copula的风电相关性条件采样方法
CN112417768B (zh) * 2020-12-07 2022-04-26 天津大学 一种基于藤结构Pair-Copula的风电相关性条件采样方法
CN112581312A (zh) * 2020-12-21 2021-03-30 国网陕西省电力公司电力科学研究院 风电功率预测误差分布分析方法、装置、计算机设备和可读存储介质
CN112581312B (zh) * 2020-12-21 2024-03-08 国网陕西省电力公司电力科学研究院 风电功率预测误差分布分析方法、装置、计算机设备和可读存储介质
CN112711864A (zh) * 2021-01-18 2021-04-27 国网浙江省电力有限公司电力科学研究院 一种配网电缆质量指标相关性模型构建及数据拓展方法
CN112685915A (zh) * 2021-01-18 2021-04-20 重庆大学 一种风电出力条件概率分布建模方法
CN112685915B (zh) * 2021-01-18 2023-06-30 重庆大学 一种风电出力条件概率分布建模方法
CN112711864B (zh) * 2021-01-18 2024-06-07 国网浙江省电力有限公司电力科学研究院 一种配网电缆质量指标相关性模型构建及数据拓展方法
CN113937769A (zh) * 2021-11-05 2022-01-14 国网甘肃省电力公司 考虑风能出力不确定性的机组组合确定方法及装置

Also Published As

Publication number Publication date
CN106485362B (zh) 2019-10-18

Similar Documents

Publication Publication Date Title
CN106485362A (zh) 一种基于高维风电预测误差模型及降维技术的发电调度方法
Jacobs et al. SOCRATES: A system for scheduling hydroelectric generation under uncertainty
Shi et al. Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features
Zhang et al. Modeling conditional forecast error for wind power in generation scheduling
CN104485665B (zh) 计及风速预测误差时空相关性的动态概率潮流计算方法
Kuo et al. A mechanistic semi-empirical wake interaction model for wind farm layout optimization
CN103713336B (zh) 基于gis分区的水电站流域面雨量气象预报方法
CN104463371B (zh) 基于风电功率变化量的Markov链建模、预测方法
CN107862466A (zh) 考虑系统双侧随机性的源荷跨时空互补效益评价方法
CN105243502A (zh) 一种基于径流区间预测的水电站调度风险评估方法及系统
CN103887813B (zh) 基于风功率预测不确定度的风电系统运行的控制方法
CN107145707A (zh) 一种计及光伏出力不确定性和全寿命周期成本的配电网变压器规划方法
CN106655190A (zh) 一种求解风电场概率最优潮流的方法
CN105512760A (zh) 基于神经网络的光伏电站发电量的计算方法及系统
Deng et al. A new wind speed scenario generation method based on spatiotemporal dependency structure
Milligan et al. Assessment of simulated wind data requirements for wind integration studies
CN106611243A (zh) 一种基于garch模型的风速预测残差修正方法
CN105634018A (zh) 一种基于随机响应面法和内点法的含风电随机最优潮流求解方法
CN115600722A (zh) 企业电碳因子长期预测方法
CN105071388B (zh) 一种基于极大似然估计的配电网状态估计方法
CN105048451B (zh) 一种基于新能源发电量区间预测的区间潮流计算方法
CN106372440A (zh) 一种并行计算的配电网自适应抗差状态估计方法及装置
CN103489037B (zh) 一种可发电风资源的预测方法
CN105741192A (zh) 一种风电场风电机机舱短期风速组合预报方法
CN114336607B (zh) 计及源荷双侧不确定性的电网经济运行域生成方法及系统

Legal Events

Date Code Title Description
C06 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