一種用於風資源評估的年風速生成方法
2023-12-02 00:06:46
一種用於風資源評估的年風速生成方法
【專利摘要】本發明涉及風力資源發電技術,特別是一種用於風資源評估的年風速生成方法,其步驟如下:首先根據當地的年月均風速、威布爾係數、自相關係數、晝夜模式強度、最大風速時刻生成一個含8760個數據點的服從均值為0、方差為1的正態分布的自相關序列;再根據每月的月均風速,以及晝夜模式強度係數、最大風速時刻,生成年8760個小時的風速初始序列;然後將風速初始序列進行概率變換後與自相關序列相加得到序列,序列在經過概率變化最後擬合得到一年8760小時的風速序列;本發明利用當地各種有效數據進行年風速擬合的方法,最終得到有效的評估數據以適應電風場的建立。
【專利說明】一種用於風資源評估的年風速生成方法
【技術領域】
[0001]本發明涉及風力資源發電技術,特別是一種用於風資源評估的年風速生成方法。【背景技術】
[0002]隨著風力發電的進一步發展,我國風電場建設逐步呈現速度快、規模大的特點,而風電場的建設需要風電場所在地測風數據處理形成的資料和長期站(氣象站、海洋站)的測風資料,伴隨著風電場的快速擴建以及新建,導致對風電場所在地的風能資源的評估的風資源數據匱乏,因此,在當前能夠獲取的有限的風資源統計參數的情況下,進行風電場所在地的年風速的生成,對風資源進行評估就顯得尤為必要。
[0003]現有風資源評估方法一般通過氣象觀測記錄的代表年或特徵年的風資源情況,進行數據統計分析,採用修正、擬合的方法得到年均或月均風速,雖然能夠反映當地的風能情況,但是年均或月均風速對於評估風電場年發電量情況意義不大。因此,在現有的氣象觀測記錄的數據基礎上,提供一種能夠生成一年8760小時或更短時間間隔的年風速,將有助於風電場風力機選型(不同風力機的風功曲線不一致,導致同樣年風速情況下,年風力發電量就不一致),提供所建風電場的預期發電量。
【發明內容】
[0004]本發明的目的在於提供一種用於風資源評估的年風速生成方法,根據年月均風速、威布爾係數、自相關係數、晝夜模式強度、最大風速時刻進行年風速擬合的方法,可以得到有效的評估數據以適應電風場的建立。
[0005]為實現上述發明目的,本發明的技術方案如下:
一種用於風資源評估的年風速生成方法,其特徵在於包括以下步驟:
第一步,根據當地氣象統計數據得到一年12個月的月均風速,以及威布爾係數、自相關係數、晝夜模式強度係數、最大風速時刻;
第二步,生成一個含8760個數據點的服從均值為O、方差為I的正態分布的自相關序列
Z ;
第三步,根據每月的月均風速,以及晝夜模式強度係數、最大風速時刻,生成年8760個小時的風速初始序列I/
第四步,將第三步中得到的風速初始序列I/進行概率變換,得到與第二步的自相關序列Z服從相同的正態分布序列的Cf ;
第五步,將第四步得到的概率變換後的序列(6/與第二步產生的自相關數列Z相加,得到新的符合正態分布的序列,;
第六步,採用第四步的概率變換方法,將第五步產生的序列IT做概率變換,得到服從以第一步所述的威布爾係數作為參數的威布爾分布序列Hr,由此便得到了擬合所得一年8760小時的風速序列IT。[0006]第二步中所述自相關序列Z的生成模型為:
【權利要求】
1.一種用於風資源評估的年風速生成方法,其特徵在於包括以下步驟: 第一步,根據當地氣象統計數據得到一年12個月的月均風速,以及威布爾係數、自相關係數、晝夜模式強度係數、最大風速時刻; 第二步,生成一個含8760個數據點的服從均值為O、方差為I的正態分布的自相關序列Z ; 第三步,根據每月的月均風速,以及晝夜模式強度係數、最大風速時刻,生成年8760個小時的風速初始序列U 第四步,將第三步中得到的風速初始序列σ進行概率變換得到序列I/,序列I/與第二步中得到的自相關序列Z服從相同的正態分布; 第五步,將第四步得到的概率變換後的序列?/與第二步產生的自相關數列Z相加,得到新的符合正態分布的序列紀; 第六步,採用第四步的概率變換方法,將第五步產生的序列IT做概率變換,得到服從以第一步所述的威布爾係數作為參數的威布爾分布序列Wr,由此便得到了擬合所得一年8760小時的風速序列JT
2.根據權利要求1所述的一種用於風資源評估的年風速生成方法,其特徵在於:第二步中所述自相關序列Z的生成模型為:
3.根據權利要求1所述的一種用於風資源評估的年風速生成方法,其特徵在於:第三步中所述風速初始序列V中第況月中任意一天第i小時的風速生成模型如下:
4.根據權利要求1所述的一種用於風資源評估的年風速生成方法,其特徵在於: 第四步中所述概率變換的方法如下:將風速初始序列I/記為序列A,A=IA1, -,A8760I,期望分布序列的累計分布函數為P(B),求期望分布序列B: ①求序列A的累計概率分布函數P(A):根據序列A的最大值maxA,劃分為N個區間,得到[O, maxA/N)、[maxA/N, maxA*2/N)、…、[maxA*( N-l) /N, maxA];分別統計每個區間段的數據個數,除以8760便得到這個區間段的數據概率,數據在第i (I < i < N)個數據段的累計概率分布為第I到(1-Ι)個數據段的概率之和,由此得到P(A); ②根據已知序列A及P(A)、期望分布序列的累計分布函數P(B),求序列B:對於Aj(l ( i ( 8760),根據P(A),得到Aj對應的累計概率P (Aj),從期望分布序列的累計分布函數P(B)中,找到累計概率為P(Ap的B值,記為由此反覆便可得到期望分布的序列B=IB1,...,B_},所得的序列B即為序列U。
【文檔編號】G06F19/00GK103886185SQ201410077792
【公開日】2014年6月25日 申請日期:2014年3月5日 優先權日:2014年3月5日
【發明者】舒軍, 田軍, 劉徵宇, 劉靜波, 吳建東, 唐健 申請人:中國東方電氣集團有限公司