CN109978222A - 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 - Google Patents
一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 Download PDFInfo
- Publication number
- CN109978222A CN109978222A CN201910009520.1A CN201910009520A CN109978222A CN 109978222 A CN109978222 A CN 109978222A CN 201910009520 A CN201910009520 A CN 201910009520A CN 109978222 A CN109978222 A CN 109978222A
- Authority
- CN
- China
- Prior art keywords
- probability
- bayesian network
- wind
- climbing
- node
- 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
- 230000009194 climbing Effects 0.000 title claims abstract description 146
- 230000005611 electricity Effects 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000005457 optimization Methods 0.000 claims abstract description 30
- 238000004422 calculation algorithm Methods 0.000 claims description 42
- 238000011156 evaluation Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 18
- 238000010845 search algorithm Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 9
- 239000002245 particle Substances 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 230000006870 function Effects 0.000 claims description 6
- 230000001568 sexual effect Effects 0.000 claims description 6
- 235000015170 shellfish Nutrition 0.000 claims description 6
- 238000009795 derivation Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 241001269238 Data Species 0.000 claims description 2
- 238000013277 forecasting method Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000009826 distribution Methods 0.000 description 6
- 238000012549 training Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000000714 time series forecasting Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 125000002015 acyclic group Chemical group 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000000979 retarding effect Effects 0.000 description 1
- 238000009987 spinning 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/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
- G06F18/24155—Bayesian classification
-
- 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Data Mining & Analysis (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统,根据所掌握的观测样本数据,挖掘风电爬坡事件与风速、风向、温度、气压、湿度等相关气象影响因子间的相依关系,搭建与样本数据拟合程度最高的贝叶斯网络拓扑结构;定量描述爬坡事件与各气象因子间的条件相依关系,估计贝叶斯网络各节点处的条件概率表内各项条件概率的取值,与贝叶斯网络拓扑结构共同组成风电爬坡事件预测的贝叶斯网络模型;由所掌握的预测时刻的数值天气预报信息,推断爬坡事件各状态发生的条件概率;自适应地调整各节点处相应条件概率的取值,从而优化推断出的爬坡事件各状态发生的条件概率结果,实现预测结果可靠性与敏锐性的折中。
Description
技术领域
本发明属于风电功率预测领域,尤其涉及一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统。
背景技术
随着风电在电力系统中渗透率的持续增长,风电出力所固有的随机性、波动性及不确定性对电网安全稳定运行、经济调度和保护控制造成的影响日益严峻。风电场有功出力在短时间内发生的大幅度变化称作风电爬坡事件。我国并网风电具有大规模、高度集中的特点,当风电穿透功率超过一定值时,无预期的风电爬坡事件便会直接导致电力系统发用电不平衡,易造成系统频率越限、恶化电能质量、威胁安全运行,甚至引起切负荷或大面积停电等事故,造成重大经济损失。实现风电爬坡事件的精确量化与准确预警,对辅助调度部门优化常规机组出力,合理配置旋转备用,减缓风/光功率波动对电网的冲击具有重要的理论价值和迫切的现实意义。
在风电爬坡事件预测技术方面,目前使用较多的方法可大致分为间接预测法和直接预测法两类。间接预测法指先借助风电功率预测技术对风电功率时间序列进行预测,再就风电功率时间序列预测结果利用爬坡事件的定义来检测风电爬坡事件的发生,并提取相应爬坡特征量的预测信息,是当前风电爬坡事件预测的主流方法。然而,经典的风电功率预测模型往往有意忽略极端气象样本,并得到相对平滑的功率时序预测曲线以最小化整体预测误差。因此,依赖于功率预测结果的间接法或许会低估爬坡事件的严重程度,或导致爬坡预警的漏报。直接预测法利用历史样本数据挖掘相应爬坡特征量与区域气象信息间的相依关系,建立由气象信息到爬坡特征量的映射,无需进行风电功率预测。该类方法更加直观、预测精度高,但其模型的训练依赖于大量历史观测样本,样本的完备性和数据的准确性直接影响了该类方法的预测性能。此外,爬坡事件的罕见性导致了观测样本的有限性,缺乏足够的训练样本或许会给确定性爬坡事件预测带来难以避免的统计误差。所以,风电爬坡事件观测样本的稀缺性严重制约了现阶段确定性风电爬坡事件预测方法的应用与发展。在此背景下,结合贝叶斯网络理论与非精确概率统计技术,实现风电爬坡事件发生概率区间的可靠估计与爬坡状态的可靠预警,为解决有限样本条件下爬坡事件的有效预测提供了新的思路。
贝叶斯网络是一种表达不确定性知识和进行因果推理的非精确概率图模型,其将图模型与非精确概率理论结合,具有表达更一般不确定性问题和进行不确定性推理的能力。面对电力调度系统获取可靠风电出力预测信息的要求,将贝叶斯网络理论用于风电爬坡事件的预测是非常有利的,尤其在历史爬坡事件观测样本有限的条件下,贝叶斯网络统计技术可以智能挖掘爬坡事件与多气象变量间的隐含的相依关系,并以非精确条件概率的形式量化相依性程度,在获取相关气象证据的条件下实现爬坡事件发生概率的可靠推理,从而提高爬坡事件发生概率的预测精度。
发明内容
为了解决现有风电爬坡事件预测技术受风电功率预测精度的影响较大,对训练样本非完备性和数据量测非准确性预测情景的适应性差的问题,本申请提供了一种基于贝叶斯网络的风电爬坡事件概率预测方法以及一种基于贝叶斯网络的风电爬坡事件概率预测系统。
为实现上述发明目的,本申请具体采用以下技术方案。
一种基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于,所述风电爬坡事件概率预测方法包括以下步骤:
步骤1:根据所掌握的地区电网实测数据,由贝叶斯网络结构学习算法从历史观测样本数据中挖掘风电爬坡事件与风速、风向、温度、气压、湿度等相关气象影响因子间的相依关系,搭建与观测样本数据拟合程度最高的贝叶斯网络拓扑结构,其中所述观测样本数据是指历史观测样本中的各气象影响因子-风电爬坡事件数据对;
步骤2:定量描述爬坡事件与各气象因子间的条件相依关系,从而初步估计步骤1中所搭建的贝叶斯网络拓扑结构中各节点处的条件概率表(CPT)内各项条件概率也即贝叶斯网络的参数的取值,构建风电爬坡事件预测的贝叶斯网络模型;
步骤3:依据步骤2所构建的贝叶斯网络模型,借助贝叶斯网络推断算法,由所掌握的预测时刻的数值天气预报信息,推断爬坡事件各状态发生的条件概率;
步骤4:对预测的条件概率结果的可靠性与敏锐性进行评价,并根据预测者对预测结果的可靠性与敏锐性的风险态度,借助粒子群优化算法自适应地调整步骤2中初步估计的贝叶斯网络参数的取值,从而优化依步骤3所推断出的爬坡事件各状态发生的条件概率,实现预测结果可靠性与敏锐性的折中。
所述基于贝叶斯网络的风电爬坡事件概率预测方法进一步包括以下优选方案。
在步骤1中,首先对地区电网实测数据进行预处理,剔除坏数据、补全缺失数据,对各实测数据实现等频离散划分,由连续量转化为状态量。
在步骤1中,结合最大权重生成树(MWST)算法与贪婪搜索(GS)算法搭建与观测样本数据拟合程度最高的贝叶斯网络拓扑结构。
其中,借助最大权重生成树算法构建与当前具体化样本数据集相匹配的树形结构,获取关于气象变量各要素与风电爬坡事件相依关系的启发性知识,借此初始化贪婪搜索算法,将初始拓扑结构限制在全局最优解的邻域内。
在步骤1中,包括以下内容:
1.1对电网实测数据进行预处理,剔除坏数据、补全缺失数据,并对风电爬坡事件的状态进行定义;
由下式对风电爬坡事件的各状态进行定义:
其中,Pt表示时刻t观测的风电功率,Pt-30min表示前一观测时刻t-30min观测的风电功率,PR为装机容量。
1.2经等频离散化过程分别将风速(S)、风向(D)、温度(T)、湿度(h)、风速变化量(V)以及爬坡事件即功率变化率(H)各划分为3个状态,组成贝叶斯网络的六个节点变量;
1.3借助最大权重生成树算法初步挖掘1.2中离散化所得的六个节点变量之间的相依关系,并搭建与观测样本数据集相匹配的初始树形结构;
1.4通过贪婪搜索算法对1.3得到的初始树形结构进行局部更新,形成一系列候选模型,然后由式(2)所示的BIC评分函数对一系列候选模型与观测样本数据之间的拟合程度进行逐一评价,其中,评价得分最优的候选模型成为“当前模型”;
式(2)中,n为节点个数,对于第i个节点Xi,假设其有ri个状态(x1,x2,...xk,...,xri),其父节点集写作Pa(Xi),共有qi种组合状态;N为观测样本的总数,Nijk则代表观测样本中符合Xi=xk且Pa(Xi)取第j种组合状态时的样本个数;Nij则表示符合Pa(Xi)取第j种组合状态,而Xi任意取值的样本个数,其中,父节点是指:对于网络拓扑结构中的两节点变量A与B,若它们之间的连接线由A指向B,则称A为B的父节点。
1.5在后续的每一次迭代中,均由贪婪搜索算法通过增边、减边、旋转边等操作对“当前模型”局部更新得到一系列候选模型,由式(2)的BIC指标对每一候选模型进行评价;
1.6若本次迭代中最优候选模型的得分优于“当前模型”的得分,则将此最优候选模型作为“当前模型”,并返回步骤1.5继续搜索;否则,停止搜索,且“当前模型”即为由贪婪搜索算法所构建的最优贝叶斯网络拓扑结构。
在步骤2中,采用改进的非精确狄利克雷模型(IDM)来初步估计贝叶斯网络的参数,包括以下步骤:
2.1依据步骤1所构建的与观测样本数据拟合程度最高的贝叶斯网络拓扑结构,确定各节点变量的父节点集合,确定各节点处待估计的条件概率表的内容;
2.2使用改进的非精确狄利克雷模型估计各节点处待估计的条件概率表内各项条件概率的取值,从而完成贝叶斯网络参数的初步估计:
基于改进的非精确狄利克雷模型估计条件概率P(A1|B1)的取值范围的计算公式如下:
其中,上式中的参数含义解释如下:节点变量A共有na个相互独立的取值,节点变量B共有nb个相互独立的取值,通过统计观测样本数据可得,节点变量B取值B1的样本个数为M,节点变量A取值的样本个数分别为节点变量A取值的概率分别为u为引入外生参数。
在步骤3中,按照贝叶斯公式,结合步骤1构建的贝叶斯网络拓扑结构与步骤2估计的贝叶斯网络参数,推导爬坡事件各状态发生的条件概率的计算公式,具体包括以下内容:
3.1根据气象变量的已知状态,写出爬坡事件待估计的条件概率表达式;
3.2经由贝叶斯公式、链式法则化简待估计的条件概率表达式;
3.3步骤2估计的贝叶斯网络参数带入化简后的待估计的条件概率表达式,计算爬坡事件发生的条件概率。
在步骤4中,根据预测者对爬坡事件各状态发生的概率预测结果的可靠性和敏锐性的风险态度,由粒子群优化算法自适应地调整步骤2所估计的贝叶斯网络中各节点处各项条件概率的取值,从而优化步骤3推断出的爬坡事件各状态的发生概率,满足预测者对结果可靠性与敏锐性的要求。
具体包括以下内容:
4.1使用指标score1对预测结果的可靠性进行评价;
预测结果的可靠性由指标score1来评价,表征预测结果对真实概率的覆盖能力,score1得分的计算方式如式(10)所示:
其中,El为由预测时刻的数值天气预报信息确定的各气象变量的所属状态,为预测出的爬坡事件状态Hw发生的概率区间,P*(Hw|El)为该项爬坡状态发生的真实概率;指标score1得分越高,则预测结果越可靠。
4.2使用指标score2对预测结果的敏锐性进行评价;
预测结果的敏锐性由指标score2来评价,表征预测结果的不确定性的程度,score2得分越低,预测结果的敏锐性越高:
4.3预测者自由指定指标score1与指标score2在综合评价指标score中所占的权重wt1与wt2,以体现个人的风险态度,对步骤3推断出的爬坡事件各状态发生概率的预测结果进行综合评价;
预测结果的综合性能由式(12)所示的score1和score2的加权和来评估:
score指标得分越高,则模型的整体预测性能越优异。
4.4使用粒子群优化算法(PSO),以最大化综合评价指标score的得分为目标,迭代优化步骤2中改进的非精确狄利克雷模型引进的超参数u的取值,从而得出与预测者风险态度相匹配的爬坡事件概率预测结果。
在粒子群优化算法的每一步迭代中,首先在超参数u的当前取值下由步骤2中改进的非精确狄利克雷模型估计贝叶斯网络各节点处条件概率表内各项条件概率的取值范围;再由步骤3推导出的爬坡事件条件概率计算公式,在已知相关气象变量状态取值的条件下,推断未知的爬坡事件变量各状态发生的概率区间;接着,由步骤3.4中设计的综合评价指标score对预测结果的可靠性与敏锐性进行综合评价;若此次迭代所得的评价指标score的得分优于上次迭代中score的得分,则将u的当前取值作为最优取值,计算参数u下一步的搜索方向与更新步长,跳转到下一步迭代中继续寻优;若此次迭代所得的评价指标score的得分未能优于上次迭代中score的得分,则停止迭代过程,输出当前u的取值为最优取值,且此次迭代中计算所得的爬坡事件各状态发生的条件概率区间即为优化后的预测结果。
本申请还提供了一种利用前述风电爬坡事件概率预测方法的风电爬坡事件概率预测系统,技术方案如下。
一种利用前述风电爬坡事件概率预测方法的风电爬坡事件概率预测系统,包括贝叶斯网络拓扑结构学习模块、贝叶斯网络参数估计模块、爬坡事件的条件概率推断模块以及模型预测效果评价与自适应参数优化模块,其特征在于:
贝叶斯网络拓扑结构学习模块,其用于构建与所掌握的观测样本数据拟合程度最高的网络拓扑结构,以智能提取并抽象表达爬坡事件与多个气象变量之间的条件相依性关系,所述贝叶斯网络拓扑结构学习模块的输出信息分别上传至贝叶斯网络参数估计模块和爬坡事件的条件概率推断模块;
所述贝叶斯网络参数估计模块读入由贝叶斯网络拓扑结构学习模块传递的变量间父子节点关系,用于由所掌握的观测样本数据估计贝叶斯网络各节点处条件概率表内各项条件概率的取值范围,以量化各变量间的条件相依性程度;所述贝叶斯网络参数估计模块的输出信息将上传至爬坡事件的条件概率推断模块;
爬坡事件的条件概率推断模块,其用于在获取相关气象变量证据后,对爬坡事件各状态发生的条件概率区间进行推导计算;
模型预测效果评价与自适应参数优化模块,其用于对所生成的爬坡事件概率预测结果的可靠性与敏锐性进行综合评价,并依据预测者的风险态度自适应地优化预测结果。
所述风电爬坡事件概率预测系统进一步包括以下优选方案:
所述贝叶斯网络拓扑结构学习模块采集包含历史风电出力以及风速、风向、温度、湿度、风速变化量气象信息在内的观测样本数据,对原始数据中的连续量信息进行离散化处理,经由最大权重生成树算法与贪婪搜索算法,输出与观测样本数据拟合程度最高的贝叶斯网络拓扑结构,以表征挖掘出的变量间的相依性关系。
所述贝叶斯网络拓扑结构学习模块将输出信息上传至贝叶斯网络参数估计模块,提供变量间的父子节点关系,为确定各节点处需估计的条件概率的内容提供参考。
所述贝叶斯网络拓扑结构学习模块还将输出信息上传至爬坡事件的条件概率推断模块,提供未知变量状态取值的概率推断所需的条件独立性信息。
所述贝叶斯网络参数估计模块读入由贝叶斯网络拓扑结构学习模块传递的变量间父子节点关系,并采集观测样本数据用于对待估计的变量间各项条件概率的统计估计,输出网络中各节点处的条件概率表。
所述爬坡事件的条件概率推断模块读入从贝叶斯网络拓扑结构学习模块传递的变量间条件独立性信息,以及从贝叶斯网络参数估计模块传递的网络中各节点处的条件概率表内的数值信息,化简爬坡事件条件概率表达式,输出爬坡事件各状态发生的条件概率预测结果。
该模块需读入从爬坡事件的条件概率推断模块传递的爬坡事件条件概率预测信息,并采集观测样本数据以对爬坡事件各状态发生的真实概率进行统计,计算综合评价指标的得分,并以最大化该评分指标的取值为目标对超参数u进行迭代寻优。在每一步的迭代过程中,该模块输出对当前参数u取值下贝叶斯网络模型预测结果的评价得分,并计算下一步迭代前对参数u的取值所做的更改。若该模块的计算结果并没有满足参数寻优的收敛判据,则该模块输出的关于参数u的更新信息将回传至贝叶斯网络参数估计模块,更新所搭建的贝叶斯网络模型的参数;若收敛判据得以满足,则该模块的输出信息将不再回传,迭代寻优过程结束,超参数u的当前取值即为与预测者风险态度相适应的最优取值。
与现有技术相比,本发明的有益效果是:
(1)经典的风电功率预测方法有意忽略极端样本,并获得相对平滑的功率时序预测曲线以最小化整体预测误差。因此,传统的从预测的风电功率序列中检测爬坡事件发生的预测方法可能会低估爬坡事件发生的概率。相比之下,本发明直接探索了不同气象条件下风电爬坡事件各状态的概率分布,因此可以避免在风电功率预测环节引入的累积误差。
(2)本发明将风电爬坡事件视为多状态随机变量,其状态根据爬坡幅度划分。由非精确概率质量函数表示的预测结果可以根据预测时刻的气象条件有效地估计每个可能的爬坡状态的概率区间,并量化预测中的不确定性,从而为调度运行人员提供更全面的预测信息。
(3)大多数基于场景分类的预测方法均是根据气象条件对观测样本进行划分,以找出易触发爬坡事件的极端气象条件。然而,对应于极端气象的样本量通常很少,或将导致预测结果不可靠。本发明通过提取爬坡事件与相关气象变量之间的相依关系,可以增加用于条件概率预测的有效样本容量(针对某一气象条件进行爬坡概率预测时,并不局限于此气象条件下的观测样本),使得在有限样本条件下也得以提高爬坡事件预测的可靠性。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为本申请预测理论中用于展示贝叶斯网络推理过程的简单三节点网络拓扑结构示意图;
图2为本申请的基于贝叶斯网络的风电爬坡事件概率预测方法流程示意图;
图3为本申请的基于MWST-GS算法的贝叶斯网络拓扑结构学习算法处理流程示意图;
图4为本申请实例验证中由MWST算法搭建的树形结构及由GS算法获取的最优拓扑结构示意图;
图5为本申请实例验证中在8种气象条件下预测的爬坡事件各状态的非精确条件概率结果示意图;
图6是本申请的基于贝叶斯网络的风电爬坡事件概率预测系统的结构示意图。
具体实施方式
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本申请中使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
本申请实施例对宁夏地区的某风电场进行了爬坡事件的概率预测,下面结合附图与实施例对本申请作进一步说明。
如图2所示,基于贝叶斯网络的风电爬坡事件概率预测方法,主要包括如下步骤:
步骤1:根据所掌握的观测样本数据,由贝叶斯网络结构学习算法挖掘风电爬坡事件与风速、风向、温度、气压、湿度等相关气象影响因子间的相依关系,搭建与样本数据拟合程度最高的贝叶斯网络拓扑结构。
按照图3示出的由最大权重生成树算法初始化的贪婪搜索算法搭建贝叶斯网络拓扑结构的处理流程,本实施例在搭建最优贝叶斯网络拓扑结构时的主要实施方式如下:
1.1对宁夏某风电场的量测数据(包括输出功率、风速、风向、温度、湿度)进行预处理,剔除坏数据、补全缺失数据;
1.2经等频离散划分将各变量的取值由连续量转化为状态量;
按照式(1)定义风电爬坡事件的各个状态。此外,5种气象变量,也即风速(S)、风向(D)、温度(T)、湿度(h)、风速变化量(V)均由等频离散化过程划分为3个状态,如表1所示。
表1贝叶斯网络中节点变量的状态
1.3由最大权重生成树算法初步挖掘六变量间的相依关系,搭建与观测样本数据集相匹配的初始树形结构;
在六节点的网络结构空间中搜索最优拓扑结构是一个NP困难的问题,故本申请拟采用最基本的启发式搜索算法—贪婪搜索(GS)算法以搜寻与样本数据拟合程度最高的网络拓扑结构。
贪婪搜索算法中,初始网络结构可以是空图、随即图或利用专家知识构造的先验结构,然而,不合理的初始结构或导致搜索迭代过程复杂化甚至陷入局部最优解。故在本申请中,首先借助最大权重生成树(MWST)算法搭建与具体化样本数据集相匹配的树形结构,获取关于各变量相依关系的启发性知识,借此初始化贪婪搜索算法,将初始拓扑结构限制在全局最优解的邻域内。
最大权重生成树算法是一种基于约束测试的结构学习算法,首先计算两两变量间的联合概率分布,使用互信息函数(式13)度量节点间的相依性程度,将此作为各节点间连接边的权重。
依各边权重从大到小排序,并从类节点出发添加权重最大的一条边,然后按照有向无环图原则逐步添加剩余边中权重最大的边,形成的n-1边生成树即是针对当前历史样本集由最大权重生成树算法构建的最大权重生成树。本申请实施例针对宁夏某风电场的观测样本数据构建的最大权重生成树结构如图4(a)所示。
1.4应用贪婪搜索算法处理1.3得到的最大权重生成树结构,生成一系列候选模型,并逐一评价候选模型与观测样本数据之间的拟合程度;
贪婪搜索算法由一个初始模型出发开始搜索,每一步搜索中首先利用搜索算子(增边、减边、反转边)对当前模型进行局部更新,然后由式(2)所示的BIC评分指标评价所得的一系列候选模型,量化候选模型与观测样本数据之间的拟合程度。其中,评价得分最优的候选模型成为“当前模型”。
1.5在贪婪搜索算法后续的每一次迭代中,均由BIC指标对本次迭代中局部更新生成的候选模型进行评价;
1.6在每一次迭代中,若评价后选出的最优候选模型的得分优于“当前模型”的得分,则将此最优候选模型作为“当前模型”,并返回步骤1.5继续搜索;否则,停止搜索,并输出“当前模型”作为由贪婪搜索算法构建的最优贝叶斯网络拓扑结构。本申请实施例针对宁夏某风电场的观测样本数据构建的最优贝叶斯网络拓扑结构如图4(b)所示。
步骤2:采用改进的非精确狄利克雷模型,初步估计步骤1所搭建的网络拓扑结构中各节点处的相关条件概率。
2.1确定步骤1所搭建的贝叶斯网络拓扑结构中各节点的父节点集合,明确各节点处待估计的条件概率的内容;
根据图4(b)所示的最优拓扑结构,可知,节点V的父节点集合为{节点H},节点S的父节点集合为{节点H与节点V},节点D的父节点集合为{节点H},节点T的父节点集合为{节点H与节点S},节点h的父节点集合为{节点D与节点T},节点H无父节点。
由此,节点V处待估计的条件概率的内容包括:P(Vy|Hw),其中,y=1,2,3;w=1,2,3。
节点S处待估计的条件概率的内容包括:P(Sr|Hw),P(Sr|Vy),其中,y=1,2,3;r=1,2,3;w=1,2,3。
节点D处待估计的条件概率的内容包括:P(Dq|Hw),其中,q=1,2,3;w=1,2,3。
节点T处待估计的条件概率的内容包括:P(Tk|Hw),P(Tk|Sr),其中,r=1,2,3;k=1,2,3;w=1,2,3。
节点h处待估计的条件概率的内容包括:P(hd|Dq),P(hd|Tk),其中,q=1,2,3;k=1,2,3;d=1,2,3。
2.2使用改进的非精确狄利克雷模型估计2.1列出的待估计的条件概率的取值,完成贝叶斯网络参数的初步估计。
在估计2.1中列出的各项条件概率时,通过对观测样本数据的统计估计,便可得出式(4)中参数m1与M的相应取值,随后计算式(4)便得出2.1中各项条件概率的取值范围,也即实现了贝叶斯网络参数的初步估计。
下面以条件概率P(A1|B1)为例,介绍非精确狄利克雷模型对此条件概率进行统计估计的理论,并详细给出本申请对非精确狄利克雷模型所做的改进。
假设节点变量A共有na个相互独立的取值,节点变量B共有nb个相互独立的取值,通过统计观测样本数据可得,节点变量B取值B1的样本个数为M,节点变量A取值的样本个数分别为节点变量A取值的概率分别为则根据非精确狄利克雷模型估计理论,条件概率P(A1|B1)的取值范围可写作:
其中,参数s为非精确狄利克雷模型为解决有限样本估计的非精确性而引入的超参数,其取值大小反应了非精确狄利克雷模型对条件概率估计所引入的不确定性的程度。
然而,本申请在非精确狄利克雷模型的实际应用中发现,对于传统经典模型所设置的常值参数s,当满足父节点取值条件要求的观测样本数量较多时(即M取值较大时),在条件概率估计中所引入的不确定性s/(M+s)衰减过快,使得步骤3中所推断的风电爬坡事件各状态发生的非精确条件概率结果的可靠性难以满足要求。因此,考虑将恒定的参数s替换为满足父节点取值条件要求的观测样本数量M的函数,即s=f(logM),以延缓非精确条件概率估计的不确定性的衰减速度,而这一延缓作用的强度则通过引入外生参数u来控制,即令s=u×logM。参数u的取值将在步骤4中对爬坡事件各状态发生概率的预测结果的评价中由粒子群优化算法自适应地优化,以实现预测结果的可靠性与敏锐性的折中。
基于改进的非精确狄利克雷模型估计条件概率P(A1|B1)的取值范围的计算公式如下:
凭借本申请在改进的非精确狄利克雷模型中所引入的可调参数u,在步骤4中便可依据预测者的风险态度自适应地调整步骤2中估计的各节点处条件概率表内各项条件概率的取值,从而实现在步骤3中对爬坡事件各状态发生概率的预测结果的自适应优化。
步骤3:根据所构建的贝叶斯网络,在获取相关气象证据的条件下,实现爬坡事件各状态发生的非精确条件概率估计。
为详细展示由贝叶斯网络拓扑结构和贝叶斯网络参数推断爬坡事件发生的条件概率的计算方法,下面将借助图1所示的简单三节点网络拓扑结构,依次经由贝叶斯公式、链式法则、以及拓扑结构中隐含的条件独立性对待估计的条件概率表达式进行变形、拆分以及化简,推导根据贝叶斯网络模型,在掌握已知变量取值的证据条件下,推断未知变量取值的算法。
3.1根据已知变量的取值,确定已知变量的所属状态,写出未知变量与已知变量之间待估计的条件概率表达式;
在图1所示的简单三节点网络拓扑结构中,假设A,C为两状态节点,B为三状态节点,也即A={Ai|i=1,2},B={Bk|k=1,2,3},C={Cd|d=1,2}。在已掌握的证据条件{Bk,Cd}下,未知变量A取值A1的条件概率可表示为P(A1|Bk,Cd),即为本示例部分待估计的条件概率。
3.2根据贝叶斯公式将3.1中的条件概率表达式改写:
根据贝叶斯公式,此条件概率表达式可写作:
根据链式法则将3.2中变形后的条件概率表达式拆分;
根据链式法则,概率P(Bk,Cd|Ai)可拆分为:
P(Bk,Cd|Ai)=P(Bk|Ai)·P(Cd|Ai,Bk) (6)
因此,式(5)可改写为:
按照贝叶斯网络拓扑结构中变量间的条件独立性关系将3.3中拆分后的条件概率表达式化简,完成对贝叶斯推断算法的推导;
由图1所示贝叶斯网络拓扑结构可知,节点B的父节点为A,节点C的父节点为A,节点C的取值与节点B的取值之间不存在直接相关关系,即在已知变量A取值状态的条件下,变量B与变量C之间条件独立,故式(7)中有:
P(Cd|Ai,Bk)=P(Cd|Ai) (8)
由此,拆分后的条件概率表达式可进一步化简为:
3.3将步骤2估计的贝叶斯网络参数带入3.4所推导出的待估计条件概率的计算公式,由所掌握的预测时刻的数值天气预报信息,推断爬坡事件各状态发生的条件概率。
式(9)中所含的条件概率子式P(Ai)、P(Bk|Ai)以及P(Cd|Ai)的取值范围均可由步骤2中改进的非精确狄利克雷模型对观测样本数据的统计估计计算得出,即在超参数u的某一取值下,通过式(4)估计各项条件概率子式的取值范围,再由式(9)计算出待估计条件概率P(A1|Bk,Cd)的上限和下限,便可得出由贝叶斯网络推断算法估计的未知变量在已知证据下取值的概率分布信息。
在本申请实施例中,3.1写出爬坡事件与五项气象变量之间待估计的条件概率表达式;
根据表1所确定的六变量的可取状态,预测时刻共存在243种不同的气象证据,表示为El={Vy,Sr,Dq,Tk,hd},y,r,q,k,d∈{1,2,3}。则实施例中待估计的爬坡事件条件概率表达式应写作:P(Hw|El)。本申请所涉及的爬坡事件概率预测方法可同时量化估计中的不确定性,故预测结果为非精确概率的形式,即
3.2根据贝叶斯公式改写3.1中的条件概率表达式:
3.3根据链式法则拆分3.2中变形后的条件概率表达式;
根据链式法则,概率Pim(Vy,Sr,Dq,Tk,hd|Hw)可拆分为:
因此,式(14)可进一步改写为:
3.4根据2.1中列出的六节点之间的父子节点关系,易知六变量间的条件独立性关系,从而待估计的爬坡事件条件概率取值的上下界与P(Hw|El)可简化为:
在式(17)中,先验概率P(Hw)可由状态Hw在观测样本数据中出现的频率来估计,其他所涉及的非精确条件概率均为直接从2.2所构建的条件概率表中读取。
步骤4:评估模型的预测效果,并将模型预测性能与对照模型作比较。依据预测者为预测结果的可靠性与敏锐性赋予的权重自适应地优化2.2中超参数u的取值。
本发明实施例所选取的对照模型为基于中心极限定理的非精确概率估计模型。
中心极限定理(CLT)是估计统计样本数据平均分布的常用方法。假设μ和σ2分别代表样本数据集的均值和方差。当由M所表示的数据集样本量足够大时,根据CLT理论,样本均值可以近似通过正态分布N(μ,σ2/M)来拟合。从而,在设定的置信水平下,可进一步获取均值的置信区间。
表2分析了在243种气象证据条件下基于CLT模型的爬坡事件预测结果和基于贝叶斯网络(BN)模型的爬坡事件预测结果的平均预测性能(CLT模型的置信度选为90%)。在此表中,综合评价指标score中的权重wt1被依次选作0.3,0.5,0.7以体现本申请的预测方法对预测者的不同风险态度的普适性。
表2预测结果分析
风险寻求者或许会指定wt1=0.3,wt2=0.7以反映他们对预测结果敏锐性的关注。在这种情况下,BN模型中所应用的PSO参数优化方法为超参数u的取值寻得的最优结果为u=11.7,BN模型预测的概率区间的平均宽度仅为0.122。超过90%的预测概率区间窄于0.2,只有极少数的区间宽于0.3,这比CLT模型的预测结果要好得多。此外,在此项测试中,78.2%的预测概率区间可以覆盖真实概率,这一比例高于CLT模型的结果,体现了BN模型在预测可靠性方面的优异性能。
相反,风险规避者或许会选择wt1=0.7,wt2=0.3来反映他们对预测结果可靠性的关注。在这种情况下,BN模型中所应用的PSO参数优化方法为超参数u的取值寻得的最优结果为u=18.2,BN模型预测结果所实现的覆盖率接近90%,体现了BN模型在可靠性方面卓越的预测性能。然而,作为达到较高可靠性的代价,此情况下预测区间的平均宽度显着增加,当wt1从0.3增加到0.7时,预测区间的平均宽度增加了50%。
有趣的是,当wt1=wt2=0.5时,BN模型的预测结果与CLT模型的预测结果恰巧实现了几乎相同的平均区间宽度。在这种情况下,BN模型中所应用的PSO参数优化方法为超参数u的取值寻得的最优结果为u=13.5,BN模型预测结果的更高的覆盖率清晰地反映了其在可靠性与敏锐性的综合评价方面相对CLT模型的优势。
总的来说,从表2展示的BN模型和CLT模型的测试结果中可以得出以下结论:
1.在测试中,综合评价指标score中的权重wt1和wt2反映了预测者个人的风险态度。权重wt1越大,所得的预测结果的可靠性越强,同时,作为高可靠性的代价,预测的概率区间将相对较宽。
2.对于不同的风险态度,可以相应地优化预测模型中的超参数u,以使模型获取更好的预测性能。从表2的预测结果分析中可以清楚地观察到,本申请的爬坡事件预测模型预测的概率区间可以根据预测者个人的风险态度进行自适应地调整,体现了该模型的灵活性。
3.从表2的分析中可见,无论预测者对可靠性与敏锐性的风险态度如何,本申请的风电爬坡事件预测模型总是表现出比CLT模型更优异的预测性能。
图5形象展示了当wt1=wt2=0.5时,在8种不同的气象条件下,本申请的BN模型和作为对照的CLT模型对爬坡事件各状态发生的非精确概率预测结果。
在图5中可以观察到,在气象条件E1,E2和E3下,通过CLT模型获得的预测区间过宽,虽覆盖住了真实概率,但过宽的预测区间也模糊了所揭示的统计规律,弱化了预测的有效性。在这三种气象条件下,BN模型在保证覆盖真实概率的基础上,总是可以获得更狭窄的预测区间,体现了本申请的爬坡事件概率预测方法优异的敏锐性。
在气象条件E4和E5下,由于观测样本数量有限,由训练集和验证集所统计的真实概率分布存在很大差异。在这两种预测情景中,CLT模型就样本充裕量方面的应用要求难以满足。在图5中可以观察到,CLT模型估计的概率区间并不能覆盖住真实概率,存在估计偏差。然而,本申请所提出的BN模型在这两种气象条件下均实现了对真实概率的覆盖,体现了该方法在有限样本条件下对揭示真实统计概率的优异性能。
此外,对于气象条件E6和E7,训练集中均无爬坡事件发生,因此,依据CLT模型的预测理论仅可得到单值预测结果P(H1|El)=1且P(H2|El)=P(H3|El)=0,带来无法容忍的预测误差。与此相反,在图5中可以观察到,本申请所提出的BN模型在这两种气象条件下仍可保证良好的预测性能,实现了对统计真实概率的覆盖。
在最不利的预测条件E8下,由于严重缺乏观测样本,训练集和验证集所统计的真实概率之间的偏差十分显著,使得BN模型和CLT模型在可靠性方面的预测性能均被劣化。然而,从图5中可以看出,BN模型预测的概率区间相对验证集统计出的真实概率的偏差更小些,体现了其与CLT模型相比,在这种非常不利的预测情景下也能更准确地揭示潜在的分布规律。
图6是本申请的基于贝叶斯网络的风电爬坡事件概率预测系统的结构示意图。
如图6所示,本申请的一种基于贝叶斯网络的风电爬坡事件概率预测系统,包括:
(1)贝叶斯网络拓扑结构学习模块,其用于构建与所掌握的观测样本数据拟合程度最高的网络拓扑结构,以智能提取并抽象表达爬坡事件与多个气象变量之间的条件相依性关系。
该模块需采集包含历史风电出力以及风速、风向、温度、湿度、风速变化量等气象信息在内的观测样本数据,对原始数据中的连续量信息进行离散化处理,经由最大权重生成树算法与贪婪搜索算法,输出与观测样本数据拟合程度最高的贝叶斯网络拓扑结构以表征挖掘出的变量间的相依性关系。该模块的输出信息将上传至贝叶斯网络参数估计模块,提供变量间的父子节点关系,为确定各节点处需估计的条件概率的内容提供参考;该模块的输出信息还将上传至爬坡事件的条件概率推断模块,提供未知变量状态取值的概率推断所需的条件独立性信息。(2)贝叶斯网络参数估计模块,其用于由所掌握的观测样本数据估计贝叶斯网络中各节点处条件概率表内各项条件概率的取值范围,以量化各变量间的条件相依性程度。
该模块需读入由贝叶斯网络拓扑结构学习模块传递的变量间父子节点关系,并采集观测样本数据用于模块内改进的非精确狄利克雷模型对待估计的变量间各项条件概率的统计估计,输出网络中各节点处的条件概率表。该模块的输出信息将上传至爬坡事件的条件概率推断模块,提供爬坡事件发生的条件概率的计算过程中所需的相关数量信息。
(3)爬坡事件的条件概率推断模块,其用于在获取相关气象变量证据后,对爬坡事件各状态发生的条件概率区间进行推导计算。
该模块需读入从贝叶斯网络拓扑结构学习模块传递的变量间条件独立性信息,以及从贝叶斯网络参数估计模块传递的网络中各节点处的条件概率表内的数值信息,经由贝叶斯公式、链式法则等对待估计的爬坡事件条件概率表达式进行改写、拆分与化简,输出爬坡事件各状态发生的条件概率预测结果。该模块的输出信息将上传至模型预测效果评价与自适应参数优化模块,提供在超参数u的当前取值下,依所搭建的贝叶斯网络模型对爬坡事件发生的概率区间的预测结果。
(4)模型预测效果评价与自适应参数优化模块,其用于对所生成的爬坡事件概率预测结果的可靠性与敏锐性进行综合评价,并依据预测者的风险态度自适应地优化预测结果。
该模块需读入从爬坡事件的条件概率推断模块传递的爬坡事件条件概率预测信息,并采集观测样本数据以对爬坡事件各状态发生的真实概率进行统计,计算综合评价指标的得分,并以最大化该评分指标的取值为目标对超参数u进行迭代选优。在每一步的迭代过程中,该模块输出对当前参数u取值下贝叶斯网络模型预测结果的评价得分,并计算下一步迭代前对参数u的取值所做的更改。若该模块的计算结果并没有达到参数寻优的收敛判据,则该模块输出的关于参数u的更新信息将回传至贝叶斯网络参数估计模块,更新所搭建的贝叶斯网络模型的参数;若收敛判据得以满足,则该模块的输出信息将不再回传,迭代寻优过程结束。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (15)
1.一种基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于,所述风电爬坡事件概率预测方法包括以下步骤:
步骤1:采集地区电网实测数据,由贝叶斯网络结构学习算法从历史观测样本数据中挖掘风电爬坡事件与风速、风向、温度、气压、湿度相关气象影响因子间的相依关系,搭建与观测样本数据拟合程度最高的贝叶斯网络拓扑结构,其中所述观测样本数据是指历史观测样本中的各气象影响因子-风电爬坡事件数据对;
步骤2:定量描述爬坡事件与各气象因子间的条件相依关系,从而初步估计步骤1中所搭建的贝叶斯网络拓扑结构中各节点处的条件概率表(CPT)内各项条件概率也即贝叶斯网络的参数的取值,构建风电爬坡事件预测的贝叶斯网络模型;
步骤3:依据步骤2所构建的贝叶斯网络模型,借助贝叶斯网络推断算法,由所掌握的预测时刻的数值天气预报信息,推断爬坡事件各状态发生的条件概率;
步骤4:对预测的条件概率结果的可靠性与敏锐性进行评价,并根据预测者对预测结果的可靠性与敏锐性的风险态度,借助粒子群优化算法自适应地调整步骤2中初步估计的贝叶斯网络参数的取值,从而优化依步骤3所推断出的爬坡事件各状态发生的条件概率,实现预测结果可靠性与敏锐性。
2.根据权利要求1所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤1中,首先对地区电网实测数据进行预处理,剔除坏数据、补全缺失数据,对各实测数据实现等频离散划分,由连续量转化为状态量。
3.根据权利要求2所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤1中,结合最大权重生成树(MWST)算法与贪婪搜索(GS)算法搭建与观测样本数据拟合程度最高的贝叶斯网络拓扑结构。
4.根据权利要求3所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤1中,借助最大权重生成树算法构建与当前具体化样本数据集相匹配的树形结构,获取关于气象变量各要素与风电爬坡事件相依关系的启发性知识,借此初始化贪婪搜索算法,将初始拓扑结构限制在全局最优解的邻域内。
5.根据权利要求1所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤1中,具体包括以下内容:
1.1对电网实测数据进行预处理,剔除坏数据、补全缺失数据,并对风电爬坡事件的状态进行定义;
由下式对风电爬坡事件的各状态进行定义:
其中,Pt表示时刻t观测的风电功率,Pt-30min表示前一观测时刻t-30min观测的风电功率,PR为装机容量。
1.2经等频离散化过程分别将风速(S)、风向(D)、温度(T)、湿度(h)、风速变化量(V)以及爬坡事件即功率变化率(H)各划分为3个状态,组成贝叶斯网络的六个节点变量;
1.3借助最大权重生成树算法初步挖掘1.2中离散化所得的六个节点变量之间的相依关系,并搭建与观测样本数据集相匹配的初始树形结构;
1.4通过贪婪搜索算法对1.3得到的初始树形结构进行局部更新,形成一系列候选模型,然后由式(2)所示的BIC评分函数对一系列候选模型与观测样本数据之间的拟合程度进行逐一评价,其中,评价得分最优的候选模型成为“当前模型”;
式(2)中,n为节点个数,对于第i个节点Xi,假设其有ri个状态(x1,x2,...xk,...,xri),其父节点集写作Pa(Xi),共有qi种组合状态;N为观测样本的总数,Nijk则代表观测样本中符合Xi=xk且Pa(Xi)取第j种组合状态时的样本个数;Nij则表示符合Pa(Xi)取第j种组合状态,而Xi任意取值的样本个数,其中,父节点是指:对于网络拓扑结构中的两节点变量A与B,若它们之间的连接线由A指向B,则称A为B的父节点;
1.5在后续的每一次迭代中,均由贪婪搜索算法通过增边、减边、旋转边等操作对“当前模型”局部更新得到一系列候选模型,由式(2)的BIC指标对每一候选模型进行评价;
1.6若本次迭代中最优候选模型的得分优于“当前模型”的得分,则将此最优候选模型作为“当前模型”,并返回步骤1.5继续搜索;否则,停止搜索,且“当前模型”即为由贪婪搜索算法所构建的最优贝叶斯网络拓扑结构。
6.根据权利要求1所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤2中,采用改进的非精确狄利克雷模型(IDM)来初步估计贝叶斯网络的参数,包括以下步骤:
2.1依据步骤1所构建的与观测样本数据拟合程度最高的贝叶斯网络拓扑结构,确定各节点变量的父节点集合,确定各节点处待估计的条件概率表的内容;
2.2使用改进的非精确狄利克雷模型估计各节点处待估计的条件概率表内各项条件概率的取值,从而完成贝叶斯网络参数的初步估计:
基于改进的非精确狄利克雷模型估计条件概率P(A1|B1)的取值范围的计算公式如下:
其中,上式中的参数含义解释如下:节点变量A共有na个相互独立的取值,节点变量B共有nb个相互独立的取值,通过统计观测样本数据可得,节点变量B取值B1的样本个数为M,节点变量A取值的样本个数分别为节点变量A取值的概率分别为u为引入外生参数。
7.根据权利要求1所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤3中,按照贝叶斯公式,结合步骤1构建的贝叶斯网络拓扑结构与步骤2估计的贝叶斯网络参数,推导爬坡事件各状态发生的条件概率的计算公式,具体包括以下内容:
3.1根据气象变量的已知状态,写出爬坡事件待估计的条件概率表达式;
3.2经由贝叶斯公式、链式法则化简待估计的条件概率表达式;
3.3步骤2估计的贝叶斯网络参数带入化简后的待估计的条件概率表达式,计算爬坡事件发生的条件概率。
8.根据权利要求1所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
在步骤4中,根据预测者对爬坡事件各状态发生的概率预测结果的可靠性和敏锐性的风险态度,由粒子群优化算法自适应地调整步骤2所估计的贝叶斯网络中各节点处各项条件概率的取值,从而优化步骤3推断出的爬坡事件各状态的发生概率,满足预测者对结果可靠性与敏锐性的要求。
9.根据权利要求8所述的基于贝叶斯网络的风电爬坡事件概率预测方法,其特征在于:
步骤4具体包括以下内容:
4.1使用指标score1对预测结果的可靠性进行评价;
预测结果的可靠性由指标score1来评价,表征预测结果对真实概率的覆盖能力,score1得分的计算方式如式(10)所示:
其中,El为由预测时刻的数值天气预报信息确定的各气象变量的所属状态,为预测出的爬坡事件状态Hw发生的概率区间,P*(Hw|El)为该项爬坡状态发生的真实概率;指标score1得分越高,则预测结果越可靠;
4.2使用指标score2对预测结果的敏锐性进行评价;
预测结果的敏锐性由指标score2来评价,表征预测结果的不确定性的程度,score2得分越低,预测结果的敏锐性越高:
4.3预测者自由指定指标score1与指标score2在综合评价指标score中所占的权重wt1与wt2,以体现个人的风险态度,对步骤3推断出的爬坡事件各状态发生概率的预测结果进行综合评价;
预测结果的综合性能由下式所示的score1和score2的加权和来评估:
score=wt1×score1-wt2×score2
s.t.wt1,wt2>0,wt1+wt2=1
score指标得分越高,则模型的整体预测性能越优异;
4.4使用粒子群优化算法(PSO),以最大化综合评价指标score的得分为目标,迭代优化步骤2中初步估计的贝叶斯网络参数的取值,从而得出与预测者风险态度相匹配的爬坡事件概率预测结果。
10.一种利用前述风电爬坡事件概率预测方法的风电爬坡事件概率预测系统,包括贝叶斯网络拓扑结构学习模块、贝叶斯网络参数估计模块、爬坡事件的条件概率推断模块以及模型预测效果评价与自适应参数优化模块,其特征在于:
贝叶斯网络拓扑结构学习模块,其用于构建与所掌握的观测样本数据拟合程度最高的网络拓扑结构,以智能提取并抽象表达爬坡事件与多个气象变量之间的条件相依性关系,所述贝叶斯网络拓扑结构学习模块的输出信息分别上传至贝叶斯网络参数估计模块和爬坡事件的条件概率推断模块;
所述贝叶斯网络参数估计模块读入由贝叶斯网络拓扑结构学习模块传递的变量间父子节点关系,用于由所掌握的观测样本数据估计贝叶斯网络各节点处条件概率表内各项条件概率的取值范围,以量化各变量间的条件相依性程度;所述贝叶斯网络参数估计模块的输出信息将上传至爬坡事件的条件概率推断模块;
爬坡事件的条件概率推断模块,其用于在获取相关气象变量证据后,对爬坡事件各状态发生的条件概率区间进行推导计算;
模型预测效果评价与自适应参数优化模块,其用于对所生成的爬坡事件概率预测结果的可靠性与敏锐性进行综合评价,并依据预测者的风险态度自适应地优化预测结果。
11.根据权利要求10所述的风电爬坡事件概率预测系统,其特征在于:
所述贝叶斯网络拓扑结构学习模块采集包含历史风电出力以及风速、风向、温度、湿度、风速变化量气象信息在内的观测样本数据,对原始数据中的连续量信息进行离散化处理,经由最大权重生成树算法与贪婪搜索算法,输出与观测样本数据拟合程度最高的贝叶斯网络拓扑结构,以表征挖掘出的变量间的相依性关系。
12.根据权利要求11所述的风电爬坡事件概率预测系统,其特征在于:
所述贝叶斯网络拓扑结构学习模块将输出信息上传至贝叶斯网络参数估计模块,提供变量间的父子节点关系,为确定各节点处需估计的条件概率的内容提供参考。
13.根据权利要求11或12所述的风电爬坡事件概率预测系统,其特征在于:
所述贝叶斯网络拓扑结构学习模块还将输出信息上传至爬坡事件的条件概率推断模块,提供未知变量状态取值的概率推断所需的条件独立性信息。
14.根据权利要求10所述的风电爬坡事件概率预测系统,其特征在于:
所述贝叶斯网络参数估计模块读入由贝叶斯网络拓扑结构学习模块传递的变量间父子节点关系,并采集观测样本数据用于对待估计的变量间各项条件概率的统计估计,输出网络中各节点处的条件概率表。
15.根据权利要求10所述的风电爬坡事件概率预测系统,其特征在于:
所述爬坡事件的条件概率推断模块读入从贝叶斯网络拓扑结构学习模块传递的变量间条件独立性信息,以及从贝叶斯网络参数估计模块传递的网络中各节点处的条件概率表内的数值信息,化简爬坡事件条件概率表达式,输出爬坡事件各状态发生的条件概率预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910009520.1A CN109978222B (zh) | 2019-01-04 | 2019-01-04 | 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910009520.1A CN109978222B (zh) | 2019-01-04 | 2019-01-04 | 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109978222A true CN109978222A (zh) | 2019-07-05 |
CN109978222B CN109978222B (zh) | 2021-04-27 |
Family
ID=67076557
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910009520.1A Active CN109978222B (zh) | 2019-01-04 | 2019-01-04 | 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109978222B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766215A (zh) * | 2019-10-17 | 2020-02-07 | 东北电力大学 | 基于特征自适应选择和wdnn的风电爬坡事件预测方法 |
CN110766069A (zh) * | 2019-10-22 | 2020-02-07 | 中国人民解放军国防科技大学 | 基于最优值迭代的多任务学习模型训练以及预测方法 |
CN110826237A (zh) * | 2019-11-11 | 2020-02-21 | 云南电网有限责任公司电力科学研究院 | 基于贝叶斯信念网络的风电设备可靠性分析方法及装置 |
CN111143573A (zh) * | 2019-12-04 | 2020-05-12 | 华东师范大学 | 基于用户反馈信息进行知识图谱目标节点预测的方法 |
CN111582557A (zh) * | 2020-04-20 | 2020-08-25 | 哈尔滨工业大学 | 基于变差函数的风电爬坡事件多级预警方法 |
CN111612253A (zh) * | 2020-05-22 | 2020-09-01 | 广西电网有限责任公司 | 基于贝叶斯网络的配电网净负荷预测方法及系统 |
CN111736690A (zh) * | 2020-05-25 | 2020-10-02 | 内蒙古工业大学 | 基于贝叶斯网络结构辨识的运动想象脑机接口 |
CN112312443A (zh) * | 2020-10-13 | 2021-02-02 | 西安电子科技大学 | 海量告警数据处理方法、系统、介质、计算机设备及应用 |
CN112464164A (zh) * | 2020-11-25 | 2021-03-09 | 北京航空航天大学 | 一种人因可靠性评估的方法、装置和信息处理设备 |
CN112488399A (zh) * | 2020-12-04 | 2021-03-12 | 国网冀北电力有限公司计量中心 | 电力负荷预测方法及装置 |
CN112598883A (zh) * | 2021-02-07 | 2021-04-02 | 国网江苏省电力有限公司苏州供电分公司 | 基于贝叶斯网络的输电线路落雷概率预警方法及预警系统 |
CN112801350A (zh) * | 2021-01-13 | 2021-05-14 | 中国电建集团华东勘测设计研究院有限公司 | 基于不确定度的深度学习超短时风电功率预测方法 |
CN113191568A (zh) * | 2021-05-21 | 2021-07-30 | 上海市气象灾害防御技术中心(上海市防雷中心) | 基于气象的城市运行管理大数据分析预测方法和系统 |
CN113505448A (zh) * | 2021-06-09 | 2021-10-15 | 上海电力大学 | 一种基于改进贝叶斯网络的风电机组动态可靠性评价方法 |
CN114781731A (zh) * | 2022-04-26 | 2022-07-22 | 成都理工大学 | 基于贝叶斯理论的滑坡运动距离超越概率预测方法和系统 |
CN117196122A (zh) * | 2023-11-02 | 2023-12-08 | 湖南赛能环测科技有限公司 | 一种基于风电功率爬坡时长的风电场调整方法及装置 |
CN117408535A (zh) * | 2023-12-15 | 2024-01-16 | 国网山东省电力公司电力科学研究院 | 一种基于深度学习的新能源爬坡事件滚动预警方法及系统 |
CN117650635A (zh) * | 2024-01-30 | 2024-03-05 | 南京金绿汇成信息科技有限公司 | 基于终端停电事件的配网拓扑分析方法及系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184337A (zh) * | 2011-06-07 | 2011-09-14 | 中国电力科学研究院 | 一种受气象信息影响的新能源发电能力动态组合分析方法 |
US20120046917A1 (en) * | 2010-08-23 | 2012-02-23 | Hsin-Fa Fang | Wind energy forecasting method with extreme wind speed prediction function |
CN103955779A (zh) * | 2014-05-23 | 2014-07-30 | 武汉大学 | 一种风电功率爬坡事件概率场景预测方法 |
US20150154504A1 (en) * | 2012-12-17 | 2015-06-04 | Arizona Board Of Regents On Behalf Of Arizona State University | Support vector machine enhanced models for short-term wind farm generation forecasting |
CN104820868A (zh) * | 2015-04-23 | 2015-08-05 | 北京恒泰实达科技股份有限公司 | 一种基于风电概率进行预测的调度展示方法 |
CN105868853A (zh) * | 2016-03-28 | 2016-08-17 | 山东大学 | 一种短期风电功率组合概率预测方法 |
CN106933778A (zh) * | 2017-01-22 | 2017-07-07 | 中国农业大学 | 一种基于爬坡事件特征识别的风电功率组合预测方法 |
CN107807309A (zh) * | 2017-10-27 | 2018-03-16 | 广东电网有限责任公司中山供电局 | 一种基于大数据的输电线路故障预警方法及系统 |
CN107846034A (zh) * | 2017-11-28 | 2018-03-27 | 甘肃省电力公司风电技术中心 | 一种平抑风电场爬坡率的储能出力控制方法 |
CN108596449A (zh) * | 2018-04-09 | 2018-09-28 | 南京邮电大学 | 一种考虑天气对配电网故障概率影响的配电网可靠性预测方法 |
-
2019
- 2019-01-04 CN CN201910009520.1A patent/CN109978222B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120046917A1 (en) * | 2010-08-23 | 2012-02-23 | Hsin-Fa Fang | Wind energy forecasting method with extreme wind speed prediction function |
CN102184337A (zh) * | 2011-06-07 | 2011-09-14 | 中国电力科学研究院 | 一种受气象信息影响的新能源发电能力动态组合分析方法 |
US20150154504A1 (en) * | 2012-12-17 | 2015-06-04 | Arizona Board Of Regents On Behalf Of Arizona State University | Support vector machine enhanced models for short-term wind farm generation forecasting |
CN103955779A (zh) * | 2014-05-23 | 2014-07-30 | 武汉大学 | 一种风电功率爬坡事件概率场景预测方法 |
CN104820868A (zh) * | 2015-04-23 | 2015-08-05 | 北京恒泰实达科技股份有限公司 | 一种基于风电概率进行预测的调度展示方法 |
CN105868853A (zh) * | 2016-03-28 | 2016-08-17 | 山东大学 | 一种短期风电功率组合概率预测方法 |
CN106933778A (zh) * | 2017-01-22 | 2017-07-07 | 中国农业大学 | 一种基于爬坡事件特征识别的风电功率组合预测方法 |
CN107807309A (zh) * | 2017-10-27 | 2018-03-16 | 广东电网有限责任公司中山供电局 | 一种基于大数据的输电线路故障预警方法及系统 |
CN107846034A (zh) * | 2017-11-28 | 2018-03-27 | 甘肃省电力公司风电技术中心 | 一种平抑风电场爬坡率的储能出力控制方法 |
CN108596449A (zh) * | 2018-04-09 | 2018-09-28 | 南京邮电大学 | 一种考虑天气对配电网故障概率影响的配电网可靠性预测方法 |
Non-Patent Citations (1)
Title |
---|
欧阳庭辉 等: "风电功率爬坡事件预测时间窗选取建模", 《中国电机工程学报》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110766215A (zh) * | 2019-10-17 | 2020-02-07 | 东北电力大学 | 基于特征自适应选择和wdnn的风电爬坡事件预测方法 |
CN110766215B (zh) * | 2019-10-17 | 2022-04-01 | 东北电力大学 | 基于特征自适应选择和wdnn的风电爬坡事件预测方法 |
CN110766069B (zh) * | 2019-10-22 | 2022-03-18 | 中国人民解放军国防科技大学 | 基于最优值迭代的多任务学习模型训练以及预测方法 |
CN110766069A (zh) * | 2019-10-22 | 2020-02-07 | 中国人民解放军国防科技大学 | 基于最优值迭代的多任务学习模型训练以及预测方法 |
CN110826237A (zh) * | 2019-11-11 | 2020-02-21 | 云南电网有限责任公司电力科学研究院 | 基于贝叶斯信念网络的风电设备可靠性分析方法及装置 |
CN110826237B (zh) * | 2019-11-11 | 2024-01-23 | 云南电网有限责任公司电力科学研究院 | 基于贝叶斯信念网络的风电设备可靠性分析方法及装置 |
CN111143573A (zh) * | 2019-12-04 | 2020-05-12 | 华东师范大学 | 基于用户反馈信息进行知识图谱目标节点预测的方法 |
CN111143573B (zh) * | 2019-12-04 | 2023-04-28 | 华东师范大学 | 基于用户反馈信息进行知识图谱目标节点预测的方法 |
CN111582557A (zh) * | 2020-04-20 | 2020-08-25 | 哈尔滨工业大学 | 基于变差函数的风电爬坡事件多级预警方法 |
CN111612253A (zh) * | 2020-05-22 | 2020-09-01 | 广西电网有限责任公司 | 基于贝叶斯网络的配电网净负荷预测方法及系统 |
CN111736690A (zh) * | 2020-05-25 | 2020-10-02 | 内蒙古工业大学 | 基于贝叶斯网络结构辨识的运动想象脑机接口 |
CN111736690B (zh) * | 2020-05-25 | 2023-07-14 | 内蒙古工业大学 | 基于贝叶斯网络结构辨识的运动想象脑机接口 |
CN112312443A (zh) * | 2020-10-13 | 2021-02-02 | 西安电子科技大学 | 海量告警数据处理方法、系统、介质、计算机设备及应用 |
CN112464164A (zh) * | 2020-11-25 | 2021-03-09 | 北京航空航天大学 | 一种人因可靠性评估的方法、装置和信息处理设备 |
CN112488399B (zh) * | 2020-12-04 | 2024-04-23 | 国网冀北电力有限公司计量中心 | 电力负荷预测方法及装置 |
CN112488399A (zh) * | 2020-12-04 | 2021-03-12 | 国网冀北电力有限公司计量中心 | 电力负荷预测方法及装置 |
CN112801350A (zh) * | 2021-01-13 | 2021-05-14 | 中国电建集团华东勘测设计研究院有限公司 | 基于不确定度的深度学习超短时风电功率预测方法 |
CN112598883A (zh) * | 2021-02-07 | 2021-04-02 | 国网江苏省电力有限公司苏州供电分公司 | 基于贝叶斯网络的输电线路落雷概率预警方法及预警系统 |
CN113191568A (zh) * | 2021-05-21 | 2021-07-30 | 上海市气象灾害防御技术中心(上海市防雷中心) | 基于气象的城市运行管理大数据分析预测方法和系统 |
CN113191568B (zh) * | 2021-05-21 | 2024-02-02 | 上海市气象灾害防御技术中心(上海市防雷中心) | 基于气象的城市运行管理大数据分析预测方法和系统 |
CN113505448A (zh) * | 2021-06-09 | 2021-10-15 | 上海电力大学 | 一种基于改进贝叶斯网络的风电机组动态可靠性评价方法 |
CN114781731A (zh) * | 2022-04-26 | 2022-07-22 | 成都理工大学 | 基于贝叶斯理论的滑坡运动距离超越概率预测方法和系统 |
CN114781731B (zh) * | 2022-04-26 | 2023-04-18 | 成都理工大学 | 基于贝叶斯理论的滑坡运动距离超越概率预测方法和系统 |
CN117196122A (zh) * | 2023-11-02 | 2023-12-08 | 湖南赛能环测科技有限公司 | 一种基于风电功率爬坡时长的风电场调整方法及装置 |
CN117196122B (zh) * | 2023-11-02 | 2024-02-06 | 湖南赛能环测科技有限公司 | 一种基于风电功率爬坡时长的风电场调整方法及装置 |
CN117408535A (zh) * | 2023-12-15 | 2024-01-16 | 国网山东省电力公司电力科学研究院 | 一种基于深度学习的新能源爬坡事件滚动预警方法及系统 |
CN117408535B (zh) * | 2023-12-15 | 2024-04-12 | 国网山东省电力公司电力科学研究院 | 一种基于深度学习的新能源爬坡事件滚动预警方法及系统 |
CN117650635B (zh) * | 2024-01-30 | 2024-04-09 | 南京金绿汇成信息科技有限公司 | 基于终端停电事件的配网拓扑分析方法及系统 |
CN117650635A (zh) * | 2024-01-30 | 2024-03-05 | 南京金绿汇成信息科技有限公司 | 基于终端停电事件的配网拓扑分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN109978222B (zh) | 2021-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109978222A (zh) | 一种基于贝叶斯网络的风电爬坡事件概率预测方法及系统 | |
CN109165664A (zh) | 一种基于生成对抗网络的属性缺失数据集补全与预测方法 | |
CN108733966A (zh) | 一种基于决策树群的多维电能表现场状态检验方法 | |
CN110796368B (zh) | 基于贝叶斯网络的社区配电网动态风险评估方法及装置 | |
CN102938708B (zh) | 基于告警传播模式的告警相关性分析系统及其分析方法 | |
CN112633316B (zh) | 一种基于边界估值理论的负荷预测方法和装置 | |
CN107274105A (zh) | 基于线性判别分析的多属性决策树电网稳定裕度评估方法 | |
CN106600115A (zh) | 一种企业信息系统运维智能分析方法 | |
CN112749904A (zh) | 基于深度学习的配电网故障风险预警方法及系统 | |
CN112132430B (zh) | 一种配电主设备分布式状态传感器可靠性评估方法及系统 | |
CN107483487B (zh) | 一种基于topsis的多维网络安全度量方法 | |
CN112464996B (zh) | 基于LSTM-XGBoost的智能电网入侵检测方法 | |
CN103617447B (zh) | 智能变电站的评价系统及评价方法 | |
CN103093097A (zh) | 基于规范切的电力系统脆弱断面识别方法 | |
CN109636066A (zh) | 一种基于模糊时间序列挖掘的风电输出功率预测方法 | |
CN114579407B (zh) | 一种因果关系检验和微服务指标预测报警方法 | |
CN103107902A (zh) | 基于决策树的攻击检测系统 | |
CN109491339B (zh) | 一种基于大数据的变电站设备运行状态预警系统 | |
CN113792754A (zh) | 一种先除异后修复的换流变dga在线监测数据处理方法 | |
CN109711664A (zh) | 一种基于大数据的输变电设备健康评估系统 | |
CN116432123A (zh) | 一种基于cart决策树算法的电能表故障预警方法 | |
CN114357670A (zh) | 一种基于bls和自编码器的配电网用电数据异常预警方法 | |
CN114154766A (zh) | 雷电动态预测下电网的易损性预警方法及系统 | |
CN116842459B (zh) | 一种基于小样本学习的电能计量故障诊断方法及诊断终端 | |
CN106100870A (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 |