一种时序风速模拟方法技术

本发明专利技术为一种时序风速模拟方法，该方法主要通过将历史实际风速分段求相邻两个时间序列变化量的累积分布函数来体现风速的时序性，即相邻时间序列的变化情况，可以准确地表现出与实际风速相似的分布状况。根据历史风速数据中出现次数最少的n个风速点，作为分割点，整个历史风速的变化量都被统计，避免了以往风速建模时风速分区间选取对模拟精度影响较大的缺点。

【技术实现步骤摘要】
一种时序风速模拟方法
本专利技术属于风速建模
，具体涉及一种时序风速模拟方法，根据历史风速数据的时序风速变化状态，建立新的模拟风速模型。
技术介绍
近年来，随着风电行业的不断发展，风力发电并网问题也日益严峻。其中一个主要问题是地区的风速建模与实际风速的分布有差距。目前，风速模拟大多采用威布尔分布点抽样的方法，尽管可以完整地表征实际风速的概率密度，但却缺少对风速时序性的描述。这使得在风力发电并网规划的结果存在一定的偏差。风速的不确定性，是风速建模的核心问题，天气不同、季节不同、温度变化等都对风速的变化产生影响，所以统计以往的风速数据，根据以往的时序风速分析风速变化规律，进而模拟出新的时序风速。基于此，文献《蒋平，霍雨翀，张龙，罗建裕，李海峰.基于改进一阶马尔可夫链的风速时间序列模型[J].电力系统自动化，2014，38(19)：22-27》用一阶马尔科夫链计算风速时序变化规律，根据不同的季节、温度、时间计算不同的转换概率矩阵，进而模拟出更为符合实际的风速，但马尔科夫链方法受风速区间选择多少的影响较大，设置不同的区间即设置不同的场景，求取场景之间的转换概率来确定下一时间序列的风速，这就导致相邻数个时间序列的风速可能一模一样，与实际风速的随时变化不符，所以模拟的风速变化与实际的风速变化依然存在差距。
技术实现思路
本专利技术的目的在于，提供一种全新的可有效解决上述技术问题的时序风速模拟方法。该方法主要通过将历史实际风速分段求相邻两个时间序列变化量的累积分布函数来体现风速的时序性，即相邻时间序列的变化情况，该方法可以准确地表现出与实际风速相似的分布状况。为达到本专利技术目的，采用如下技术方案：一种时序风速模拟方法，其特征在于该方法包括以下步骤：步骤1、获取同地区一段时间的历史风速数据；步骤2、分析这一段时间的历史风速数据得出一年的历史风速分布分割点a1、a2、a3…an，并确定历史风速的最大风速Vmax；步骤3、将历史风速数据根据风速分割点将(0，Vmax)的风速范围划分为n+1个风速区间，分别为(0，a1)、(a1，a2)、(a2，a3)…(an，Vmax)；提取每个风速区间的风速数据，计算每个风速区间中各个风速与下一个时间序列的风速变化量ΔV，从而获取n+1组风速变化量数据，将n+1组风速变化量对应用X1、X2、X3…Xn+1表示；步骤4、分别求取n+1组风速变化量数据即X1、X2、X3…Xn+1的累积分布函数，累积分布函数用式(2)表示；Fi(Δv)＝P(ΔV≤Δv)(i∈(1，n+1))(2)其中，Δv为自变量，ΔV为小于Δv的风速变化量，P为ΔV落在(0，Δv)区间的概率；Fi(Δv)为Xi的累积分布分数；步骤5、求取n+1组风速变化量数据累积分布函数的相应反函数；步骤6、计算整个风速历史数据的平均值作为初始风速V1，判断V1所在n+1组风速区间为(aj-1，aj)，j∈(1，n+1)，随机生成0到1之间的数b1，带入集合Xj的累积分布函数的反函数中，得到一个风速变化量ΔV1，将其加到V1上得到风速V2，再判断V2所在风速区间的位置(ak-1，ak)，k∈(1，n+1)，随机生成0到1之间的数b2，带入集合Xk的累积分布函数的反函数，得到一个风速变化量ΔV2，将其加到V2上得到风速V3，如此迭加N次，N为模拟的风速的点数，从而生成一段时间的模拟风速。步骤1中，所述历史风速数据的时间尺度为一年，一年的历史数据点时间间隔为10min。步骤2中，风速历史风速分布分割点的选取，以该地区历史风速数据中出现概率较低的风速点作为风速分布分割点，a1、a2、a3…an的数值依次增大。风速a1、a2、a3…an风速点是根据地区数年的历史风速分析得到的，n取决于不同地区实际情况。n为大于2的整数。步骤5中，求取累积分布函数的相应反函数，即得到反函数曲线，横坐标为各小于变化量ΔV的分布概率，范围为0到1，纵坐标为变化量ΔV。与现有技术相比，本专利技术有益效果在于：(1)本方法与以往威布尔等随机抽样方法相比，模拟了风速的时序性，符合实际风速在时间序列上变化的特点。(2)本方法通过统计相邻时序风速的变化量，建立累积分布函数，进而模拟出了该地区与选取的历史实际风速数据相似分布情况的一段风速序列。(3)本方法根据历史风速数据中出现次数最少的n个风速点，作为分割点，但分割点的选取对于模拟风速的准确度影响较小，在于整个历史风速的变化量都被统计，即分割点n取大于2的值即可，所以避免了以往风速建模时风速分区间选取对模拟精度影响较大的缺点。附图说明图1为某地区历史的5000个实际风速点曲线。图2为根据历史风速基于马尔科夫链模拟出来的5000个风速点曲线。图3为根据历史风速基于本专利技术方法模拟出来的5000个风速点曲线。图4为本专利技术方法流程图。具体实施方式下面给出本专利技术的具体实施例。具体实施例仅用于进一步详细说明本专利技术，不限制本申请的保护范围。本专利技术一种风速时序模拟方法(参见图4)包括如下步骤：步骤1、获取同地区一年的历史风速数据M；步骤2、找出历史风速数据中出现次数最少的n个风速点，确定这几个风速点为风速分布分割点，a1、a2、a3…an的数值依次增大；步骤3、将一年的历史风速数据根据风速分割点划分为n+1个风速区间，分别为(0，a1)、(a1，a2)、(a2，a3)…(an，Vmax)，其中Vmax是历史风速的最大风速，提取每个风速区间的风速数据，计算每个区间中各个风速与下一个时间序列的风速变化量ΔV，从而获取n+1组风速变化量数据，将n+1组风速变化量用X1、X2、X3…Xn+1表示。如t时刻风速为Va，t+1时刻风速为Vb，且Va∈(ai-1，ai)，计算ΔVXi如式(1)所示，ΔVXi∈Xi。ΔVXi＝Vb-Va(1)其中X1是区间(0，a1)中每个风速点与下一个时间序列之间变化量的集合，因为Va∈(ai-1，ai)，所以记为ΔVXi(ΔVXi∈Xi)。步骤4、通过MATLAB软件进行计算，用ecdf函数分别求取n+1组风速变化量数据即X1、X2、X3…Xn+1的累积分布函数，累积分布函数如式(2)；Fi(Δv)＝P(ΔV≤Δv)(i∈(1，n+1))(2)其中Δv为自变量，ΔV为小于Δv的风速变化量。P为ΔV落在(0，Δv)区间的概率。步骤5、通过createFit函数求取n+1组风速变化量数据累积分布函数的反函数，如式(3)；Fi-1(P)＝Δv(i∈(1，n+1))(3)步骤6、计算整个风速历史数据的平均值作为初始风速V1，判断V1所在n+1组风速区间为(aj-1，aj)，j∈(1，n+1)，随机生成0到1之间的数b1，带入集合Xj的累积分布函数的反函数中，得到一个风速变化量ΔV1，将其加到V1上得到风速V2，再判断V2所在风速区间的位置(ak-1，ak)，本文档来自技高网...

【技术保护点】
1.一种时序风速模拟方法，其特征在于该方法包括以下步骤：/n步骤1、获取同地区一段时间的历史风速数据；/n步骤2、分析这一段时间的历史风速数据得出一年的历史风速分布分割点a

【技术特征摘要】
1.一种时序风速模拟方法，其特征在于该方法包括以下步骤：
步骤1、获取同地区一段时间的历史风速数据；
步骤2、分析这一段时间的历史风速数据得出一年的历史风速分布分割点a1、a2、a3…an，并确定历史风速的最大风速Vmax；
步骤3、将历史风速数据根据风速分割点将(0，Vmax)的风速范围划分为n+1个风速区间，分别为(0,a1)、(a1,a2)、(a2,a3)…(an,Vmax)；提取每个风速区间的风速数据，计算每个风速区间中各个风速与下一个时间序列的风速变化量ΔV，从而获取n+1组风速变化量数据，将n+1组风速变化量对应用X1、X2、X3…Xn+1表示；
步骤4、分别求取n+1组风速变化量数据即X1、X2、X3…Xn+1的累积分布函数，累积分布函数用式(2)表示；
Fi(Δv)＝P(ΔV≤Δv)(i∈(1,n+1))(2)
其中，Δv为自变量，ΔV为小于Δv的风速变化量，P为ΔV落在(0，Δv)区间的概率；Fi(Δv)为Xi的累积分布分数；
步骤5、求取n+1组风速变化量数据累积分布函数的相...

【专利技术属性】
技术研发人员：张家安，刘东，王琨玥，
申请(专利权)人：河北工业大学，
类型：发明
国别省市：天津;12