CN111027790B - 一种多个风电场出力的时间序列构建方法及*** - Google Patents
一种多个风电场出力的时间序列构建方法及*** Download PDFInfo
- Publication number
- CN111027790B CN111027790B CN201811171127.4A CN201811171127A CN111027790B CN 111027790 B CN111027790 B CN 111027790B CN 201811171127 A CN201811171127 A CN 201811171127A CN 111027790 B CN111027790 B CN 111027790B
- Authority
- CN
- China
- Prior art keywords
- hidden
- probability
- wind power
- state
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000010276 construction Methods 0.000 title claims abstract description 12
- 238000009826 distribution Methods 0.000 claims abstract description 130
- 239000011159 matrix material Substances 0.000 claims abstract description 110
- 230000007704 transition Effects 0.000 claims abstract description 100
- 239000000203 mixture Substances 0.000 claims abstract description 51
- 230000001186 cumulative effect Effects 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims description 54
- 239000013598 vector Substances 0.000 claims description 43
- 238000012549 training Methods 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 7
- 230000005251 gamma ray Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 4
- 238000000342 Monte Carlo simulation Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 10
- 238000004590 computer program Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000003860 storage Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 241000135164 Timea Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- 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/067—Enterprise or organisation modelling
-
- 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/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Wind Motors (AREA)
Abstract
本发明提供了一种多个风电场出力的时间序列构建方法及***,包括:基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列。本发明采用混合高斯隐马尔科夫模型描述多个风电场之间的时空相关性,采用多维混合高斯分布描述多个风电场出力在不同相关性状态下的联合概率分布,通过蒙特卡洛仿真生成具有时空相关性的多风电场出力时间序列,提高了风电场出力时间序列的准确性。
Description
技术领域
本发明涉及新能源发电技术领域,具体涉及一种多个风电场出力的时间序列构建方法及***。
背景技术
随着风电装机容量占比在电力***中的增加,使得风电出力的随机性和间歇性特点对电力***的规划和运行带来的挑战也逐渐增加。目前,电力***的时序生产模拟方法是应对风电出力随机性的有效手段,已广泛应用于电力***中长期规划和运行领域,如电力***风电容量规划、***调峰裕度评估和电力电量平衡计算等。但时序生产模拟方法需要以大量长时间尺度的风电场出力时间序列为输入数据,并且风电场出力时间序列的准确性对时序生产模拟方法得到的计算结果有着重要的影响。
随着风电装机容量的增加,电网中所接入的风电场个数不断增加,而目前已有的时间序列建模方法是对多个风电场出力进行单独的建模,并未考虑不同风电场出力之间的时空相关性,由于距离相近的风电场通常受同一天气过程的影响,其出力之间具有较强的时空相关性。因此,不考虑风电场之间的相关性因素而进行独立建模,势必会影响时间序列建模结果的准确性。实际观测结果表明,风电场之间的时空相关性并非固定不变的,而是具有时变特性,这是因为全年中风电场会受到多种类型天气过程的影响,由于不同类型天气过程的形成及演化机理的不同,对风电场出力的时空相关性的影响也会不同。因此,如何提高风电场出力时间序列的准确性有待解决。
发明内容
为了解决现有技术中所存在已有的时间序列建模方法是对多个风电场出力进行单独的建模,并未考虑不同风电场出力之间的时空相关性、风电场之间时空相关性的时变性,以及在生成多个风电场出力时间序列时,需要基于不同的时空相关性状态下多个风电场出力的联合概率分布问题,本发明提供一种多个风电场出力的时间序列构建方法及***,采用混合高斯隐马尔科夫模型描述多个风电场之间的时空相关性,通过状态转移特性刻画时空相关性的时变特点,采用多维混合高斯分布描述多个风电场出力在不同相关性状态下的联合概率分布,通过蒙特卡洛仿真生成具有时空相关性的多风电场出力时间序列,准确描述风电出力的变化规律和随机特性,以及不同风电场之间的时空相关性,提高生成出力时间序列的准确性,为电力***中长期规划和运行模拟提供所需的基础数据,并且所得到的时间序列场景比现有的时间序列建模方法通过独立建模得到的时间序列更加科学、合理。
本发明提供的技术方案是:一种多个风电场出力的时间序列构建方法,包括:
基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;
基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;
由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列。
优选的,所述基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵,包括:
基于多个风电场的出力数据确定相关性状态;
将所述相关性状态设置为所述混合高斯隐马尔可夫模型中隐状态,并设置所述隐状态的个数;
基于多个风电场的出力数据和所述隐状态的个数,通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态对应的混合高斯分布的权值、均值向量和协方差矩阵;
基于每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵确定对应的多维混合高斯分布;
基于所述状态转移矩阵计算累积状态转移概率;
其中,所述累积状态转移概率矩阵由累积状态转移概率组成。
优选的,所述基于多个风电场的出力数据和所述隐状态的个数通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态对应的高斯分布的权值、均值向量和协方差矩阵,包括:
步骤101、为混合高斯隐马尔科夫模型中的状态转移矩阵,混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量设置初始值;
步骤102、计算训练周期内初始时刻下各隐状态的前向概率,并依次向前计算训练周期内其余所有时刻下各隐状态的前向概率;
步骤103、计算训练周期内终止时刻下各隐状态的后向概率,并依次向后计算训练周期内其余所有时刻下各隐状态的后向概率;
步骤104、基于多个风电场的出力数据构成的观测序列、各隐状态的前向概率和各隐状态的后向概率,计算任意时刻隐状态的条件概率;
步骤105、基于多个风电场的出力数据构成的观测序列、各隐状态的前向概率和各隐状态的后向概率,计算任意时刻隐状态和所述任意时刻的下一时刻隐状态的联合条件概率;
步骤106、基于所述条件概率和联合条件概率更新混合高斯隐马尔科夫模型中的状态转移矩阵,以及混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量;
步骤107、若更新后的混合高斯隐马尔科夫模型中的状态转移矩阵,以及混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量满足预设的收敛条件,则结束循环,按当前参数设置混合高斯隐马尔科夫模型,否则执行步骤102。
优选的,所述初始时刻下各隐状态的前向概率,按下式计算:
α1(j)=πj·bj(O1),j=1,2,...,N
式中:α1(j):初始时刻下隐状态为j的前向概率;πj:风电场的隐状态为θj的概率;bj(·):多个风电场在隐状态为j的多维混合高斯分布;O1:初始时刻下各风电场出力的历史观测序列;N:隐状态的个数;
所述训练周期内其余所有时刻下各隐状态的前向概率,按下式计算:
式中:αt+1(j):t+1时刻下隐状态为j的前向概率;αt(j):t时刻下隐状态为j的前向概率;aij:状态转移矩阵A中第i行第j列t时刻的隐状态为θi,t+1时刻的隐状态转换为θj的概率;Ot+1:t+1时刻下各风电场出力的历史观测序列;
其中,所述初始时刻下风电场的隐状态为θj的概率组成初始时刻隐状态的概率分布向量。
优选的,所述训练周期内终止时刻下各隐状态的后向概率,按下式计算:
βT(j)=1,j=1,2,...,N
式中:βT(j):终止时刻T下隐状态为j的后向概率;N:隐状态的个数;
所述训练周期内其余所有时刻下各隐状态的后向概率,按下式计算:
式中:βt(j):t时刻下隐状态为j的后向概率;aji:状态转移矩阵A中第j行第i列t时刻的隐状态为θj,t+1时刻的隐状态转换为θi的概率;bi(·):多个风电场在隐状态为j的多维混合高斯分布;βt+1(i):t+1时刻下隐状态为i的后向概率;Ot+1:t+1时刻下各风电场出力的历史观测序列。
优选的,所述任意时刻隐状态的条件概率,按下式计算:
式中:γt(j,k):t时刻下隐状态j的第k个高斯分布的条件概率;αt(j):t时刻下隐状态为j的前向概率;βt(j):t时刻下隐状态为j的后向概率;αt(i):t时刻下隐状态为i的前向概率;βt(i):t时刻下隐状态为i的后向概率;隐状态为j的第k个高斯分布的权值;/>隐状态为j的第k个高斯分布的均值向量;/>隐状态为j的第k个高斯分布的协方差矩阵;隐状态为j的第n个高斯分布的权值;/>隐状态为j的第n个高斯分布的均值向量;/>隐状态为j的第n个高斯分布的协方差矩阵;N:隐状态的个数;K:隐状态为j对应的混合高斯分布中的高斯分布个数。
优选的,所述联合条件概率,按下式计算:
式中:ξt(i,j):t时刻下隐状态j的第k个高斯分布的联合条件概率;αt(j):t时刻下隐状态为j的前向概率;aij:状态转移矩阵A中第i行第j列t时刻的隐状态为θi,t+1时刻的隐状态转换为θj的概率;bj(·):多个风电场在隐状态为j的多维混合高斯分布;Ot+1:t+1时刻下各风电场出力的历史观测序列;βt+1(j):t+1时刻下隐状态为j的后向概率;αt(r):t时刻下隐状态为r的前向概率;ars:状态转移矩阵A中第r行第s列t时刻的隐状态为θr,t+1时刻的隐状态转换为θs的概率;bs(·):隐状态为s的观测序列服从的概率分布;βt+1(s):t+1时刻下隐状态为s的后向概率;N:隐状态的个数。
优选的,所述基于所述条件概率和联合条件概率更新混合高斯隐马尔科夫模型中的状态转移矩阵,以及混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量,如下式所示:
式中:πj:风电场的相关状态为θj的概率;γt(j,k):t时刻下隐状态j的第k个高斯分布的条件概率;N:隐状态的个数;
式中:aij:状态转移矩阵A中第i行第j列t时刻的隐状态为θi,t+1时刻的隐状态转换为θj的概率;ξt(i,j):t时刻下隐状态j的第k个高斯分布的联合条件概率;γt(j,n):t时刻下隐状态j的第n个高斯分布的条件概率;K:隐状态为j对应的混合高斯分布中的高斯分布个数;
优选的,所述多维混合高斯分布,按下式计算:
优选的,所述累积状态转移概率按下式计算:
式中:dij:累积状态转移概率矩阵D中第i行第j列的累积状态转移概率;aik:状态转移矩阵A中第i行第k列t时刻的隐状态为θi,t+1时刻的隐状态转换为θk的概率。
优选的,所述基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力,包括:
基于所述累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的相关性状态;
基于每一时刻所述多个风电场的相关性状态和对应的多维混合高斯分布抽样生成所述时刻下多个风电场的出力。
优选的,所述基于所述累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的相关性状态,包括:
基于隐状态的个数随机生成一个正整数作为预设周期中初始时刻多个风电场的相关性状态;
基于初始时刻多个风电场的相关性状态与所述累积状态转移概率矩阵迭代计算获得预设周期中其余时刻多个风电场的相关性状态。
优选的,所述基于初始时刻多个风电场的相关性状态与所述累积状态转移概率矩阵迭代计算获得预设周期中其余时刻多个风电场的相关性状态,包括:
步骤201、将所述初始时刻作为已知时刻;
步骤202、生成区间[0,1]内的随机数;
步骤203、所述随机数与所述累积状态转移概率矩阵中已知时刻对应行的所有元素进行比较,当随机数的大小位于所述已知时刻对应行相邻两列元素之间时,则确定下一时刻多个风电场的相关性状态为所述已知时刻对应行相邻两列中前一列对应的值;
步骤204、当所述下一时刻为预设周期的终止时刻时,则循环结束,输出预设周期内所有时刻对应的多个风电场的相关性状态;否则,将所述下一时刻作为已知时刻执行步骤202。
优选的,所述多个风电场出力的时间序列,如下式所示:
式中:Θ:预设周期T内多个风电场出力的时间序列;pM(tT):T时刻风电场M的出力。
基于同一发明构思本发明还提供了一种多风电场出力的时间序列构建***,包括:
处理模块,用于基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;
出力模块,用于基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;
出力序列模块,用于由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列。
优选的,所述处理模块,包括:
相关性状态子模块,用于基于多个风电场的出力数据确定相关性状态;
设置子模块,用于将所述相关性状态设置为所述混合高斯隐马尔可夫模型中隐状态,并设置所述隐状态的个数;
迭代子模块,用于基于多个风电场的出力数据和所述隐状态的个数通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵;
第一计算子模块,用于基于每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵确定对应的多维混合高斯分布;
第二计算子模块,用于基于所述状态转移矩阵计算累积状态转移概率;
其中,所述累积状态转移概率矩阵由累积状态转移概率组成。
与现有技术相比,本发明的有益效果为:
本发明提供的技术方案,基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列,在生成多个风电场出力时间序列时,考虑了不同的时空相关性状态下多个风电场出力的联合概率分布,提高了风电场出力时间序列的准确性。
本发明提供的技术方案,采用混合高斯隐马尔科夫模型描述多个风电场之间的时空相关性,通过状态转移特性刻画时空相关性的时变特点,采用多维混合高斯分布描述多个风电场出力在不同相关性状态下的联合概率分布,通过蒙特卡洛仿真生成具有时空相关性的多风电场出力时间序列。
本发明提供的技术方案,准确描述风电出力的变化规律和随机特性,以及不同风电场之间的时空相关性,提高生成出力时间序列的准确性,为电力***中长期规划和运行模拟提供所需的基础数据。
本发明提供的技术方案,时间序列场景比现有的时间序列建模方法通过独立建模得到的时间序列更加科学、合理。
附图说明
图1为本发明的一种多个风电场出力的时间序列构建方法流程图;
图2为本发明的混合高斯隐马尔科夫模型结构示意图;
图3为本发明实施例中两个风电场生成的出力时间序列与历史出力序列的概率分布对比示意图;
图4为本发明实施例中两个风电场生成的出力时间序列与历史出力序列的自相关系数结果对比示意图;
图5为本发明实施例中两个风电场生成出力时间序列和历史出力序列的互相关系数结果示意图。
具体实施方式
为了更好地理解本发明,下面结合说明书附图和实例对本发明的内容做进一步的说明。
因为从短时间尺度来看,风电场之间的时空相关性具有时变性,并且在生成多个风电场出力时间序列时,需要基于不同的时空相关性状态下多个风电场出力的联合概率分布。因此本发明提供了一种多个风电场出力的时间序列构建方法,首先,建立以风电场之间的相关性关系为隐状态,以同一时刻各风电场出力为观测的混合高斯隐马尔可夫模型;然后,基于多个风电场历史出力时间序列来训练模型,得到混合高斯隐马尔可夫模型的参数、相关性状态之间的转移关系、以及风电场出力的联合概率分布参数;最后,通过蒙特卡洛仿真的方法生成满足相关性关系时变特性和概率分布特性的多风电场出力时间序列。
实施例1:
图1为本发明的一种多个风电场出力的时间序列构建方法流程图,如图1所示,包括:
步骤S1、基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;
步骤S2、基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;
步骤S3、由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列。
对本发明提供的一种多个风电场出力的时间序列构建方法进一步解释:
步骤S1、基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵,具体包括:
步骤101、考虑时空相关性的混合高斯隐马尔科夫模型
本发明采用混合高斯隐马尔科夫模型描述多个风电场的时空相关性。如图2所示,混合高斯隐马尔科夫模型由隐状态序列和观测序列组成,{Q1,Q2,...,Qt}为各时刻下的状态序列,是不可观测量;{O1,O2,...,Ot}是各时刻的状态对应的观测序列。图2中的横向箭头表示在相邻两个时刻下状态之间的转换,竖向箭头表示从状态向观测量的输出。
混合高斯隐马尔科夫模型假设隐状态和观测量具有马尔可夫性,即任意时刻的隐状态仅依赖于前一时刻的隐状态,观测量仅依赖于当前时刻的隐状态。
混合高斯隐马尔科夫模型的数学模型描述如下:
1)假设隐状态的可能取值为θ1,θ2,...,θN,N为隐状态的数目,则t时刻下隐状态的取值范围为Qt∈{θ1,θ2,...,θN};假设观测量包含M个元素,即Ot为M维列向量。
2)隐状态的概率分布向量为π=[π1,π2,...,πN]T,其中
πj=P(Q1=θj),1≤j≤N
表示t1时刻的状态为θj时的概率。
3)将隐状态用离散的马尔可夫链进行描述,状态转移矩阵A=(aij)N×N,为N×N的矩阵,其中第i行第j的元素为:
aij=P(Qt+1=θj|Qt=θi),1≤i,j≤N
aij表示t时刻状态为θi,t+1时刻状态转换为θj的概率。
4)观测概率B={bj(O),j=1,2,...,N},其中
式中,bj(x)表示隐状态j的观测变量O服从的概率分布,K为状态j的观测变量的高斯分布个数(假设所有的状态的观测变量均包含K个多维高斯分布),为状态j的第k个高斯分布的权值,/>为状态j的第k个高斯分布的均值向量,/>为状态j的第k个高斯分布的协方差矩阵。/> 和表示所有N个状态对应的N·K个高斯分布的权值、均值向量和协方差矩阵的集合。
本发明中采用混合高斯隐马尔科夫模型进行多个风电场相关性的建模,具体过程为:
将多个风电场出力的相关性状态作为混合高斯隐马尔科夫模型的隐状态变量,每个时刻对应一个离散的隐状态变量,设隐状态个数为N,规定隐状态变量的取值为1至N之间的正整数,即θ1=1,θ2=2,...,θN=N;
混合高斯隐马尔科夫模型的状态转移矩阵即刻画了风电场之间时空相关关系的时变特性;
将各风电场归一化出力作为混合高斯隐马尔科夫模型的观测量,即当前时刻下所有风电场出力组成观测向量,不同状态下所有风电场出力的联合概率分布采用多维混合高斯分布进行描述。
此时,所建立的混合高斯隐马尔科夫模型可用多元组λ=(A,c,μ,Σ,π)来进行表示。
步骤102、混合高斯隐马尔科夫模型参数估计
基于时刻1,2,...,T各风电场出力的历史观测序列O={O1,O2,...,OT},对步骤101中所建立的混合高斯隐马尔科夫模型的参数λ=(A,c,μ,Σ,π)进行估计。采用前向-后向算法(Baum Welch算法)估计模型参数,其目标是在已知观测序列O的情况下,使得概率P(O|λ)最大。
Baum Welch算法的基本流程如下:
1)随机初始化所有参数λ=(A,c,μ,Σ,π);
2)计算时刻1的各个隐状态的前向概率:
α1(j)=πj·bj(O1),j=1,2,...,N
递推计算时刻2,...,T的前向概率:
3)终止时刻T的各个隐状态的后向概率:
βT(j)=1,j=1,2,...,N
递推计算时刻T-1,T-2,...,1时刻的后向概率:
4)计算在给定模型参数λ=(A,c,μ,Σ,π)和观测序列O,时刻t状态的条件概率:
5)计算在给定模型参数λ=(A,c,μ,Σ,π)和观测序列O,时刻t和t+1状态的联合条件概率:
6)更新模型参数λ=(A,c,μ,Σ,π):
7)判断更新前后的参数是否收敛,若未收敛,则重复执行步骤2)至7);若收敛,则算法结束。
通过以上的Baum Welch算法,便可得到混合高斯隐马尔科夫模型的参数λ=(A,c,μ,Σ,π)。
步骤103、计算混合高斯隐马尔科夫模型的累积状态转移概率矩阵
混合高斯隐马尔科夫模型的状态转移矩阵即描述了不同状态(时空相关性关系)之间的概率转移特性,采用蒙特卡洛仿真的方式随机生成多个风电场的出力时间序列时,需要基于混合高斯隐马尔科夫模型的累积状态转移概率矩阵。
由混合高斯隐马尔科夫模型的状态转移矩阵A=(aij)N×N,计算其累积状态转移概率矩阵D=(bij)N×(N+1),各元素的计算方式如下:
dij为累积状态转移概率矩阵D中的第i行、第j列元素,表示状态在当前时刻为θi,在下一时刻处于θ1,...,θj的概率总和。
步骤S2、基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力,包括:
步骤201、基于蒙特卡洛仿真的多个风电场出力时间序列建模,即假设所需生成的时间序列对应的时段范围为{t1,t2,...,tT},风电场个数为M,时间序列生成步骤如下:
1)基于M个风电场的历史出力数据,通过Baum Welch算法训练模型,得到混合高斯隐马尔科夫模型参数λ=(A,c,μ,Σ,π);
2)假在初始t1时刻,随机生成1至N之间的正整数q1,将其作为t1时刻M个风电场之间的相关性状态Q1=q1,基于状态Q1对应的混合高斯分布,通过随机抽样生成t1时刻M个风电场的出力[p1(t1),p2(t1),...,pM(t1)]T;
3)假设当前时刻为tl,1≤l≤T,相关性状态Ql=ql,基于[0,1]均匀分布,生成区间[0,1]内的随机数xl,将与累积状态转移概率矩阵B中第ql行的所有元素进行比较,当xl的大小位于该行第y列元素和第y+1列元素之间时,则可确定tl+1时刻相关性状态Ql+1=y;
4)基于状态Ql+1对应的混合高斯分布,通过随机抽样生成tl+1时刻M个风电场的出力[p1(tl+1),p2(tl+1),...,pM(tl+1)]T;
5)判断tl+1时刻是否为最后一个时刻tT,若不是,则重复步骤3)至5),若是,结束循环。
步骤S3、由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列,包括:
在预设时刻t1,t2,...,tT下,M个风电场出力的时间序列:
式中:Θ:预设周期T内多个风电场出力的时间序列;pM(tT):T时刻风电场M的出力。
本实施例中选取某2个风电场对所提出方法进行测试,基于全年8760小时的历史出力数据来训练混合高斯隐马尔科夫模型参数,并生成两个风电场8760h的出力时间序列。
采用概率分布指标、自相关系数指标和互相关系数指标来对所生成的风电场出力时间序列进行评价。其中,概率分布指标描述了全年风电出力在不同出力区间内的概率分布情况;自相关性指标描述了单个风电场出力随时间的变化特性;互相关系数用来描述各风电场之间出力的时空相关性,能够反映出两个风电场出力时间序列的相关性随时间的变化情况。
附图3展示了两个风电场生成的出力时间序列与历史出力序列的概率分布对比图,结果显示两个风电场生成的出力序列与历史出力序列的概率分布比较接近,证明所提出的模型能够较好的反映出风电出力的概率分布特性。
附图4展示了两个风电场生成的出力时间序列与历史出力序列的自相关系数结果,结果显示,所生成的出力时间序列的自相关系数与历史出力序列的变化趋势相一致,说明本专利所提出的建模方法能够很好的刻画风电出力的时间变化特性。
附图5展示了两个风电场所生成出力时间序列和历史出力序列的互相关系数对比情况,结果显示所生成序列的互相关系数与历史序列的互相关系数变化偏差较小,变化趋势也保持一致,说明本专利所提出的建模方法能够很好的描述不同风电场之间的时空相关性。
图3‐图5的结果均验证了本发明所提出方法的有效性。
实施例2
1)读取所有风电场的历史出力数据,并用装机容量进行归一化处理;
2)根据步骤101的建模方法,建立以所有风电场出力为观测量、以风电场之间的相关性关系为状态量的混合高斯隐马尔科夫模型参数;
3)基于风电场出力历史数据,根据步骤102中的Baum Welch算法,估计混合高斯隐马尔科夫模型的参数;
4)基于步骤103中的方法,计算所建立的混合高斯隐马尔科夫模型的累积状态转移概率矩阵;
5)给定所需生成序列的时间分辨率和时长,基于训练得到的混合高斯隐马尔科夫模型,采用步骤201中的蒙特卡洛仿真方法,随机生成所有风电场的出力时间序列。
实施例3:
基于同一种发明构思,本发明还提供了一种多风电场出力的时间序列构建***,包括:
处理模块,用于基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;
出力模块,用于基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;
出力序列模块,用于由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列。
实施例中,所述处理模块,包括:
相关性状态子模块,用于基于多个风电场的出力数据确定相关性状态;
设置子模块,用于将所述相关性状态设置为所述混合高斯隐马尔可夫模型中隐状态,并设置所述隐状态的个数;
迭代子模块,用于基于多个风电场的出力数据和所述隐状态的个数通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵;
第一计算子模块,用于基于每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵确定对应的多维混合高斯分布;
第二计算子模块,用于基于所述状态转移矩阵计算累积状态转移概率;
其中,所述累积状态转移概率矩阵由累积状态转移概率组成。
实施例中,所述出力模块,包括:
相关性状态子模块,用于基于所述累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的相关性状态;
出力子模块,用于基于每一时刻所述多个风电场的相关性状态和对应的多维混合高斯分布抽样生成所述时刻下多个风电场的出力。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上仅为本发明的实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均包含在申请待批的本发明的权利要求范围之内。
Claims (2)
1.一种多个风电场出力的时间序列构建方法,其特征在于,包括:
基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;
基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;
由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列;
所述基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵,包括:
基于多个风电场的出力数据确定相关性状态;
将所述相关性状态设置为所述混合高斯隐马尔科夫模型中隐状态,并设置所述隐状态的个数;
基于多个风电场的出力数据和所述隐状态的个数,通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态对应的混合高斯分布的权值、均值向量和协方差矩阵;
基于每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵确定对应的多维混合高斯分布;
基于所述状态转移矩阵计算累积状态转移概率;
其中,所述累积状态转移概率矩阵由累积状态转移概率组成;
所述基于多个风电场的出力数据和所述隐状态的个数通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态对应的高斯分布的权值、均值向量和协方差矩阵,包括:
步骤101、为混合高斯隐马尔科夫模型中的状态转移矩阵,混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量设置初始值;
步骤102、计算训练周期内初始时刻下各隐状态的前向概率,并依次向前计算训练周期内其余所有时刻下各隐状态的前向概率;
步骤103、计算训练周期内终止时刻下各隐状态的后向概率,并依次向后计算训练周期内其余所有时刻下各隐状态的后向概率;
步骤104、基于多个风电场的出力数据构成的观测序列、各隐状态的前向概率和各隐状态的后向概率,计算任意时刻隐状态的条件概率;
步骤105、基于多个风电场的出力数据构成的观测序列、各隐状态的前向概率和各隐状态的后向概率,计算任意时刻隐状态和所述任意时刻的下一时刻隐状态的联合条件概率;
步骤106、基于所述条件概率和联合条件概率更新混合高斯隐马尔科夫模型中的状态转移矩阵,以及混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量;
步骤107、若更新后的混合高斯隐马尔科夫模型中的状态转移矩阵,以及混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量满足预设的收敛条件,则结束循环,按当前参数设置混合高斯隐马尔科夫模型,否则执行步骤102;
所述初始时刻下各隐状态的前向概率,按下式计算:
α1(j)=πj·bj(O1),j=1,2,...,N
式中:α1(j):初始时刻下隐状态为j的前向概率;πj:风电场的隐状态为θj的概率;bj(·):多个风电场在隐状态为j的多维混合高斯分布;O1:初始时刻下各风电场出力的历史观测序列;N:隐状态的个数;
所述训练周期内其余所有时刻下各隐状态的前向概率,按下式计算:
式中:αt+1(j):t+1时刻下隐状态为j的前向概率;αt(j):t时刻下隐状态为j的前向概率;aij:状态转移矩阵A中第i行第j列t时刻的隐状态为θi,t+1时刻的隐状态转换为θj的概率;Ot+1:t+1时刻下各风电场出力的历史观测序列;
其中,所述初始时刻下风电场的隐状态为θj的概率组成初始时刻隐状态的概率分布向量;
所述训练周期内终止时刻下各隐状态的后向概率,按下式计算:
βT(j)=1,j=1,2,...,N
式中:βT(j):终止时刻T下隐状态为j的后向概率;N:隐状态的个数;
所述训练周期内其余所有时刻下各隐状态的后向概率,按下式计算:
式中:βt(j):t时刻下隐状态为j的后向概率;aji:状态转移矩阵A中第j行第i列t时刻的隐状态为θj,t+1时刻的隐状态转换为θi的概率;bi(·):多个风电场在隐状态为j的多维混合高斯分布;βt+1(i):t+1时刻下隐状态为i的后向概率;Ot+1:t+1时刻下各风电场出力的历史观测序列;
所述任意时刻隐状态的条件概率,按下式计算:
式中:γt(j,k):t时刻下隐状态j的第k个高斯分布的条件概率;αt(j):t时刻下隐状态为j的前向概率;βt(j):t时刻下隐状态为j的后向概率;αt(i):t时刻下隐状态为i的前向概率;βt(i):t时刻下隐状态为i的后向概率;隐状态为j的第k个高斯分布的权值;/>隐状态为j的第k个高斯分布的均值向量;/>隐状态为j的第k个高斯分布的协方差矩阵;/>隐状态为j的第n个高斯分布的权值;/>隐状态为j的第n个高斯分布的均值向量;/>隐状态为j的第n个高斯分布的协方差矩阵;N:隐状态的个数;K:隐状态为j对应的混合高斯分布中的高斯分布个数;
所述联合条件概率,按下式计算:
式中:ξt(i,j):t时刻下隐状态j的第k个高斯分布的联合条件概率;αt(j):t时刻下隐状态为j的前向概率;aij:状态转移矩阵A中第i行第j列t时刻的隐状态为θi,t+1时刻的隐状态转换为θj的概率;bj(·):多个风电场在隐状态为j的多维混合高斯分布;Ot+1:t+1时刻下各风电场出力的历史观测序列;βt+1(j):t+1时刻下隐状态为j的后向概率;αt(r):t时刻下隐状态为r的前向概率;ars:状态转移矩阵A中第r行第s列t时刻的隐状态为θr,t+1时刻的隐状态转换为θs的概率;bs(·):隐状态为s的观测序列服从的概率分布;βt+1(s):t+1时刻下隐状态为s的后向概率;N:隐状态的个数;
所述基于所述条件概率和联合条件概率更新混合高斯隐马尔科夫模型中的状态转移矩阵,以及混合高斯分布的权值、均值向量、协方差矩阵和隐状态的概率分布向量,如下式所示:
式中:πj:风电场的相关状态为θj的概率;γt(j,k):t时刻下隐状态j的第k个高斯分布的条件概率;N:隐状态的个数;
式中:aij:状态转移矩阵A中第i行第j列t时刻的隐状态为θi,t+1时刻的隐状态转换为θj的概率;ξt(i,j):t时刻下隐状态j的第k个高斯分布的联合条件概率;γt(j,n):t时刻下隐状态j的第n个高斯分布的条件概率;K:隐状态为j对应的混合高斯分布中的高斯分布个数;
所述多维混合高斯分布,按下式计算:
所述累积状态转移概率按下式计算:
式中:dij:累积状态转移概率矩阵D中第i行第j列的累积状态转移概率;aik:状态转移矩阵A中第i行第k列t时刻的隐状态为θi,t+1时刻的隐状态转换为θk的概率;
所述基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力,包括:
基于所述累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的相关性状态;
基于每一时刻所述多个风电场的相关性状态和对应的多维混合高斯分布抽样生成所述时刻下多个风电场的出力;
所述基于所述累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的相关性状态,包括:
基于隐状态的个数随机生成一个正整数作为预设周期中初始时刻多个风电场的相关性状态;
基于初始时刻多个风电场的相关性状态与所述累积状态转移概率矩阵迭代计算获得预设周期中其余时刻多个风电场的相关性状态;
所述基于初始时刻多个风电场的相关性状态与所述累积状态转移概率矩阵迭代计算获得预设周期中其余时刻多个风电场的相关性状态,包括:
步骤201、将所述初始时刻作为已知时刻;
步骤202、生成区间[0,1]内的随机数;
步骤203、所述随机数与所述累积状态转移概率矩阵中已知时刻对应行的所有元素进行比较,当随机数的大小位于所述已知时刻对应行相邻两列元素之间时,则确定下一时刻多个风电场的相关性状态为所述已知时刻对应行相邻两列中前一列对应的值;
步骤204、当所述下一时刻为预设周期的终止时刻时,则循环结束,输出预设周期内所有时刻对应的多个风电场的相关性状态;否则,将所述下一时刻作为已知时刻执行步骤202;
所述多个风电场出力的时间序列,如下式所示:
式中:Θ:预设周期T内多个风电场出力的时间序列;pM(tT):T时刻风电场M的出力。
2.一种多个风电场出力的时间序列构建***,用于实现一种如权利要求1所述的一种多个风电场出力的时间序列构建方法,其特征在于,包括:
处理模块,用于基于多个风电场的出力数据和预先构建的混合高斯隐马尔科夫模型得到多维混合高斯分布和累积状态转移概率矩阵;
出力模块,用于基于所述多维混合高斯分布和累积状态转移概率矩阵确定预设周期内每一时刻所述多个风电场的出力;
出力序列模块,用于由所述周期内所有时刻对应的多个风电场的出力构成所述多个风电场出力的时间序列;
所述处理模块,包括:
相关性状态子模块,用于基于多个风电场的出力数据确定相关性状态;
设置子模块,用于将所述相关性状态设置为所述混合高斯隐马尔科夫模型中隐状态,并设置所述隐状态的个数;
迭代子模块,用于基于多个风电场的出力数据和所述隐状态的个数通过前向-后向算法,迭代计算所述混合高斯隐马尔科夫模型的状态转移矩阵,以及每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵;
第一计算子模块,用于基于每一个隐状态的混合高斯分布的权值、均值向量和协方差矩阵确定对应的多维混合高斯分布;
第二计算子模块,用于基于所述状态转移矩阵计算累积状态转移概率;
其中,所述累积状态转移概率矩阵由累积状态转移概率组成。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811171127.4A CN111027790B (zh) | 2018-10-09 | 2018-10-09 | 一种多个风电场出力的时间序列构建方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811171127.4A CN111027790B (zh) | 2018-10-09 | 2018-10-09 | 一种多个风电场出力的时间序列构建方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111027790A CN111027790A (zh) | 2020-04-17 |
CN111027790B true CN111027790B (zh) | 2023-06-09 |
Family
ID=70190297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811171127.4A Active CN111027790B (zh) | 2018-10-09 | 2018-10-09 | 一种多个风电场出力的时间序列构建方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111027790B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326903A (zh) * | 2013-07-05 | 2013-09-25 | 华北电力大学 | 基于隐马尔科夫的Internet网络时延预测方法 |
CN104182914A (zh) * | 2014-09-05 | 2014-12-03 | 国家电网公司 | 一种基于波动特性的风电出力时间序列建模方法 |
CN108335010A (zh) * | 2017-12-18 | 2018-07-27 | 中国电力科学研究院有限公司 | 一种风电出力时间序列建模方法及*** |
CN108520094A (zh) * | 2018-03-07 | 2018-09-11 | 南方电网科学研究院有限责任公司 | 一种风电场出力的时空序列模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101300247B1 (ko) * | 2011-11-11 | 2013-08-26 | 경희대학교 산학협력단 | 마르코프 연쇄 은닉 조건부 랜덤 필드 모델 기반의 패턴 인식 방법 |
-
2018
- 2018-10-09 CN CN201811171127.4A patent/CN111027790B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103326903A (zh) * | 2013-07-05 | 2013-09-25 | 华北电力大学 | 基于隐马尔科夫的Internet网络时延预测方法 |
CN104182914A (zh) * | 2014-09-05 | 2014-12-03 | 国家电网公司 | 一种基于波动特性的风电出力时间序列建模方法 |
CN108335010A (zh) * | 2017-12-18 | 2018-07-27 | 中国电力科学研究院有限公司 | 一种风电出力时间序列建模方法及*** |
CN108520094A (zh) * | 2018-03-07 | 2018-09-11 | 南方电网科学研究院有限责任公司 | 一种风电场出力的时空序列模拟方法 |
Non-Patent Citations (10)
Title |
---|
Deep Gaussian Mixture-Hidden Markov Model for Classification of EEG Signals;Min Wang; Sherif Abdelfattah;《IEEE Transactions on Emerging Topics in Computational Intelligence》;第2卷(第4期);第278 – 287页 * |
Hidden Markov Models for Wind Farm Power Output;Debarati Bhaumik; Daan Crommelin;《IEEE Transactions on Sustainable Energy》;第10卷(第2期);第533 – 539页 * |
Pai Li ; Gaofeng Fan ; Weisheng Wang.Proportion optimization of wind and solar power capacity for regional power network.《2016 IEEE PES Asia-Pacific Power and Energy Engineering Conference (APPEEC)》.2016,第2424-2428页. * |
Pierre Pinson ; Henrik Madsen.Probabilistic forecasting of wind power at the minute time-scale with Markov-switching autoregressive models.《Proceedings of the 10th International Conference on Probablistic Methods Applied to Power Systems》.2009,第1-8页. * |
基于改进马尔科夫链的风电功率时间序列模型;赵宇;肖白;顾兵;王逍;张钰;王成龙;;电力建设;38(第07期);第18-24页 * |
基于新型场景划分与考虑时序相关性的光伏出力时间序列模拟方法;江雪辰;朱俊澎;袁越;王跃峰;黄阮明;;电力建设;39(第10期);第63-70页 * |
基于波动特性的风电出力时间序列建模方法研究;李驰;刘纯;黄越辉;王伟胜;《电网技术》;第39卷(第1期);第208-214页 * |
基于混合高斯隐马尔可夫模型的带式输送机堆煤时刻预测方法;钱建生,李小斌,秦文光,秦海初;《工矿自动化》;第40卷(第11期);第26-30页 * |
风电场短期风速预测若干问题研究;张驰;《中国博士学位论文全文数据库 (信息科技辑)》(第2期);I140-7 * |
马尔科夫理论及其在预测中的应用综述;黄麒元;王致杰;王东伟;杜彬;;技术与市场;22(第09期);第12-13+16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111027790A (zh) | 2020-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111027732B (zh) | 一种多风电场出力场景的生成方法及*** | |
CN109871977A (zh) | 基于离散小波变换和优化最小支持向量机的负荷预测方法 | |
CN103986156A (zh) | 一种考虑风电接入的动态概率潮流计算方法 | |
CN110413878B (zh) | 基于自适应弹性网络的用户-商品偏好的预测装置和方法 | |
CN112464567B (zh) | 基于变分同化框架的智能资料同化方法 | |
CN114897277B (zh) | 一种基于lstm的突发型滑坡位移预测方法 | |
CN112100911A (zh) | 一种基于深度bisltm的太阳辐射预测方法 | |
CN112270439B (zh) | 超短期风电功率预测方法、装置、电子设备及存储介质 | |
CN106529185A (zh) | 一种古建筑位移的组合预测方法及*** | |
CN116680540A (zh) | 一种基于深度学习的风电功率预测方法 | |
CN111985719A (zh) | 一种基于改进的长短期记忆网络的电力负荷预测方法 | |
CN111507530A (zh) | 基于分数阶动量梯度下降的rbf神经网络船舶交通流预测方法 | |
CN111563203A (zh) | 一种基于快速非负隐特征分析的智能家居用户-服务兴趣度个性化预测装置和预测方法 | |
CN112836885B (zh) | 组合负荷预测方法、装置、电子设备和存储介质 | |
CN111027790B (zh) | 一种多个风电场出力的时间序列构建方法及*** | |
CN114168574A (zh) | 一种面向工业负荷的数据缺失处理方法及装置 | |
Del Debbio et al. | Machine learning trivializing maps: A first step towards understanding how flow-based samplers scale up | |
CN112801364A (zh) | 基于miv与改进蝙蝠算法的参考作物腾发量预测方法 | |
Tao et al. | Chaotic time series prediction based on radial basis function network | |
CN111832623A (zh) | 一种基于相空间重构的回声状态网络时间序列预测算法 | |
CN109840308B (zh) | 一种区域风电功率概率预报方法及*** | |
CN117200208A (zh) | 基于多尺度分量特征学习的用户级短期负荷预测方法及*** | |
CN112348247A (zh) | 一种光能发电的预测方法、装置、计算机设备及存储介质 | |
CN117291069A (zh) | 一种基于改进de与注意力机制的lstm污水水质预测方法 | |
CN115759455A (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 |