文 | 馬友林

本文主要描述了IEC 61400-12-1 Annex D 和Annex E中不確定度（Uncertainty）相關(guān)部分的處理計算，主要包含以下三個部分：

（1）不確定度的基本概念。

（2）IEC 61400-12-1不確定度介紹。

（3）IEC 61400-12-1風(fēng)速測試不確定度。

本文主要遵循以下標(biāo)準(zhǔn)的相關(guān)內(nèi)容（下稱“標(biāo)準(zhǔn)”）：IEC 61400-12-1:2017, Wind turbine-Part 12-1: Power performance measurements of electricity producing wind turbines。

不確定度的基本概念

測量不確定度意指根據(jù)所用到的信息，表征賦予被測量值分散性的非負(fù)參數(shù)。測試不確定度是基于被測量的概率分布而估計的離散程度，用標(biāo)準(zhǔn)差表示的不確定度又被稱為標(biāo)準(zhǔn)不確定度。在某些情況下需要對不確定度進(jìn)行更可靠的估計，即對標(biāo)準(zhǔn)不確定度乘以一個包含因子k得到擴(kuò)展不確定度。

不確定度分為A、B兩類，它們的本質(zhì)是一樣的，其區(qū)別在于獲取被測量概率分布的方式不同；不能簡單地用“系統(tǒng)誤差”和“隨機(jī)誤差”進(jìn)行區(qū)分。

A類不確定度：其概率分布是直接的，通過多次測量的方式得到的。

B類不確定度：其概率分布是間接的，通過數(shù)學(xué)模型、假定等方式得到的。

將各不確定度分量進(jìn)行合并，可以得到合成不確定度。

測量模型為：

敏感系數(shù)為：

如果各分量不相關(guān)：

式中：

U(xi)為不確定度分量。

如果分量相關(guān)：

[注：不確定度的概念和計算主要參考以下標(biāo)準(zhǔn)：JJF1059.1-2012、JCGM 100: 2018 (GUM1995 with minor corrections)]。

風(fēng)速不確定度

一、概述

在測試中需要使用傳感器采集風(fēng)速相關(guān)數(shù)據(jù)，針對不同傳感器，其不確定度的計算也有所差異。一般在風(fēng)電測試中，主要采用安裝在測風(fēng)塔上的風(fēng)杯式風(fēng)速計和遙感測風(fēng)設(shè)備進(jìn)行風(fēng)信息的采集和測量。

關(guān)于各個不確定度分量，標(biāo)準(zhǔn)中已經(jīng)有相關(guān)規(guī)定。表1為IEC61400-12-1，2017版中附錄D所列的不確定度分量。

由表1可知，IEC61400-12-1功率測試中考慮的不確定度可以分為圖1所示的幾類。

圖1 風(fēng)速不確定度的分類

表1 數(shù)據(jù)篩選原則

Measured parameter Uncertainty component Uncertainty category Terrain (with site calibration)Anemometer calibration before test Post calibration/In-situ calibration Anemometer classification Mounting effects Standard mounting Alternative mounting Side mounted Lightning finial Data acquisition system (see below)Change of correction (adjacent wind direction bins)Removal of wind direction sensor between site calibration and power curve measurement Seasonal variation Statistical variability in site calibration BBBBBBBBBBBBA Method Air density correction Wind conditions - missing shear information Wind conditions - missing veer information Wind conditions - missing upflow information Wind conditions - missing turbulence information Seasonal effects Turbulence normalisation (or lack of turbulence normalisation)Cold climate measurements BBBBBBBB

表2 風(fēng)杯式風(fēng)速計所需計算的不確定度

表3 遙感測風(fēng)設(shè)備所需計算的不確定度

表4 風(fēng)速測量不確定度分量

風(fēng)杯式風(fēng)速計和遙感測風(fēng)設(shè)備的不確定度在大類別上并沒有太大的差異，但是由于設(shè)備不同，其計算方法有較大的區(qū)別。因此，在計算不確定度時，有必要將風(fēng)杯式風(fēng)速計和遙感測風(fēng)設(shè)備的不確定度分別計算。根據(jù)標(biāo)準(zhǔn)，風(fēng)杯式風(fēng)速計和遙感測風(fēng)設(shè)備的不確定度分量見表2和表3。

二、風(fēng)杯式風(fēng)速計不確定度的計算

1.概述

參照標(biāo)準(zhǔn)，風(fēng)速測量的不確定度分量見表4。

風(fēng)速測試的合成不確定度按照以下公式進(jìn)行計算：

式中：

uVS,precal,i—風(fēng)速計標(biāo)定所產(chǎn)生的不確定度；

uVS,postcal,i—采集器現(xiàn)場/后標(biāo)定的不確定度；

uVS,class,i—風(fēng)速計分級所產(chǎn)生的不確定度；

uVS,mnt,i—風(fēng)速計安裝所產(chǎn)生的不確定度；

uVS,lgt,i—避雷針?biāo)a(chǎn)生的不確定度；

uVS,mnt,i—風(fēng)速計安裝所產(chǎn)生的不確定度；

udVS,i—采集系統(tǒng)的不確定度。

2.標(biāo)定不確定度

在風(fēng)速測試的合成不確定度計算公式中，uVS,precal,i為風(fēng)速計標(biāo)定所帶來的不確定度，如圖2所示。

通過風(fēng)速計的校準(zhǔn)證書，可以查看校準(zhǔn)實驗室資質(zhì)、校準(zhǔn)日期和風(fēng)速計信息，如圖3所示。

圖2 風(fēng)速計標(biāo)定不確定度

圖3 風(fēng)速計校準(zhǔn)證書1

通過風(fēng)速計的校準(zhǔn)證書，還可以查看執(zhí)行標(biāo)準(zhǔn)、支撐桿的直徑，如圖4所示。

同時，風(fēng)速計標(biāo)定不確定度也可以從標(biāo)定報告中獲取。在這里，標(biāo)準(zhǔn)不確定度uVS,precal,i可以取成(0.05m/s)/k=0.025m/s，如圖5所示。

3.風(fēng)速計運(yùn)行特性（分級）不確定度

在風(fēng)速測試的合成不確定度計算公式中，uVS,class,i為風(fēng)速計運(yùn)行特性（分級）不確定度。

IEC61400-12-1附錄I和附錄J中列出了風(fēng)速計分級測試的相關(guān)要求。溫度、入流角、湍流度、空氣密度等在測試周期內(nèi)都會發(fā)生變化，因此，風(fēng)速計將在不同的環(huán)境條件下運(yùn)行，并會受到環(huán)境變化的影響。而風(fēng)速計分級的目的是得到風(fēng)速計在不同環(huán)境條件下的運(yùn)行特性，并在數(shù)據(jù)處理中考慮運(yùn)行特性帶來的不確定度，如圖6與圖7所示。

根據(jù)不同的使用環(huán)境應(yīng)用對應(yīng)的分級系數(shù)k計算不確定度。

例如，10m/s的風(fēng)速，使用A1.8的分級結(jié)果，對應(yīng)的運(yùn)行不確定度為：

4.風(fēng)速計安裝不確定度

在風(fēng)速測試的合成不確定度計算公式中，uVS,mnt,i為風(fēng)速計安裝所產(chǎn)生的不確定度，如圖8所示。

風(fēng)速計安裝不確定度可參照標(biāo)準(zhǔn)中的默認(rèn)值進(jìn)行計算，如表5所示。

5.避雷針不確定度

在風(fēng)速測試的合成不確定度計算公式中，uVS,lgt,i為避雷針?biāo)a(chǎn)生的不確定度，如圖9所示。按照標(biāo)準(zhǔn)要求，避雷針不確定度的默認(rèn)值為 0.1% ～ 0.2%。

6.采集系統(tǒng)不確定度

在風(fēng)速測試的合成不確定度計算公式中，udVS,i為采集系統(tǒng)的不確定度，如圖10、圖11所示。

7.后標(biāo)定/現(xiàn)場標(biāo)定不確定度

圖4 風(fēng)速計校準(zhǔn)證書2

圖5 風(fēng)速計校準(zhǔn)證書3

圖6 溫度變化對測試結(jié)果的影響

圖7 分級測試結(jié)果

圖8 風(fēng)速計安裝不確定度

圖9 避雷針不確定度

圖10 采集系統(tǒng)不確定度

圖11 采集器使用說明書

表5 風(fēng)速計安裝的不確定度

在風(fēng)速測試的合成不確定度計算公式中，uVS,postcal,i為采集器后標(biāo)定/現(xiàn)場的不確定度，如圖12所示。

利用以下公式進(jìn)行計算：

通過風(fēng)速比對可以得到現(xiàn)場標(biāo)定風(fēng)速比對關(guān)系，如圖13所示。

如果進(jìn)行了后標(biāo)定：不確定度可以估計為4～12m/s范圍內(nèi)兩次標(biāo)定差異的最大值，上限為0.2m/s。

如果進(jìn)行了現(xiàn)場標(biāo)定：不確定度可以估計為4～12m/s范圍內(nèi)最大的δ值，上限為0.2m/s，如圖14所示。

圖12 后標(biāo)定/現(xiàn)場的不確定度

圖13 現(xiàn)場標(biāo)定風(fēng)速比對關(guān)系

圖14 風(fēng)速計in-situ不確定度

8.風(fēng)速測量合成不確定度

按照以上的方法，可以得到風(fēng)速不確定度的各個分量，由各個分量可以合成總的風(fēng)速不確定度，如圖15所示。

至此，風(fēng)速測量不確定度的計算工作已經(jīng)完成。

三、遙感裝置不確定度

1.概述

在標(biāo)準(zhǔn)附錄L中，有與遙感裝置相關(guān)的敘述。附錄L是一個規(guī)范性附錄，也就是附錄中的要求都是必須要遵守的。

附錄L描述了對遠(yuǎn)程風(fēng)速測試裝置（RSD）的操作要求，例如光雷達(dá)（LIDARS）和聲雷達(dá)（SODARS）。

附錄L中規(guī)定了一個流程：確保測量的可追溯性符合國際標(biāo)準(zhǔn)并用于功率曲線評估的相關(guān)不確定性分量。

附錄L包含以下內(nèi)容：

（1）L.1 General；

（2）L.2 Classi fi cation of remote sensing devices；

（3）L.3 Verification of the performance of remote sensing device；

（4）L.4 Evaluation of uncertainty of measurements of remote sensing devices；

（5）L.5 Additional checks；

（6）L.6 Other requirements specific to power curve testing；

（7）L.7 Reporting。

2.遙感裝置不確定度的主要分量

遙感裝置不確定度主要包含圖16所示的幾個部分。

通過以下公式可以求出遙感裝置的不確定度：

圖15 風(fēng)速測量不確定度的計算

圖16 遙感裝置不確定度

式中：

uVR,ver,i—RSD驗證的不確定度；

uVR,isc,i—RSD現(xiàn)場標(biāo)定的不確定度；

uVR,class,i—RSD分級的不確定度；

uVR,mnt,i—RSD安裝的不確定度；

uVR,flow,i—RSD安裝的不確定度；

uVR,mon,i—RSD安裝的不確定度。

3.RSD驗證和分級不確定度

在遙感裝置不確定度計算公式中，uVR,ver,i為RSD驗證的不確定度，uVR,class,i為RSD分級的不確定度，如圖17所示。

若要計算RSD驗證和分級不確定度，必須先了解RSD的驗證和分級測試。參照標(biāo)準(zhǔn)中的要求，遙感裝置（RSD）使用時須遵循圖18所示的流程。

RSD需要經(jīng)過分級和驗證后才可以在測試中使用。通過分級測試，可以得到環(huán)境因素對遙感設(shè)備的影響大小，即敏感性。同時，可以得到獨(dú)立變量對于遙感設(shè)備的精度計算結(jié)果。結(jié)合相關(guān)結(jié)果，得到最終的不確定度，如圖19－圖21。

4.RSD安裝不確定度

在遙感裝置不確定度計算公式中，uVR,mnt,i為RSD安裝的不確定度1，如圖22所示。參照標(biāo)準(zhǔn)，該不確定度的默認(rèn)值為0.1%。

在遙感裝置不確定度計算公式中，uVR,flow,i為RSD安裝的不確定度2，如圖23所示。

參照標(biāo)準(zhǔn)，該不確定度的默認(rèn)值為2%～3%，可按照以下公式進(jìn)行計算：

在遙感裝置不確定度計算公式中，uVR,mon,i為RSD安裝的不確定度3，如圖24所示。參照標(biāo)準(zhǔn)，該不確定度默認(rèn)值為0.5%。

5.RSD現(xiàn)場標(biāo)定不確定度

圖17 RSD驗證和分級不確定度

圖18 RSD使用流程

圖19 精度計算方法

圖20 IEC標(biāo)準(zhǔn)中的案例

圖21 遙感裝置不確定度的計算

圖22 RSD安裝不確定度1

圖23 RSD安裝不確定度2

圖24 RSD安裝不確定度

圖25 RSD現(xiàn)場標(biāo)定不確定度

圖26 RSD合成不確定度計算表

關(guān)于RSD現(xiàn)場標(biāo)定，標(biāo)準(zhǔn)附錄K中有詳細(xì)描述，相關(guān)信息如圖25所示。

6.RSD合成不確定度

RSD合成不確定度計算表如圖26所示。

至此，RSD的不確定度計算完畢。

綜上，實際測試和分析中，需要根據(jù)實際情況，先求取不確定度的分量，再合并成最終的合成不確定度。

攝影：王鵬岳