李功克 吳 捷 姚國興

（1.華南理工大學(xué)電力學(xué)院，廣州 510640；2.廣東省綠色能源技術(shù)重點(diǎn)實(shí)驗(yàn)室，廣州 510640）

1 引言

當(dāng)今世界，能源已成為人們生活中最基本的動力。但隨著常規(guī)能源的日益枯竭，人們正在尋找綠色環(huán)?？稍偕履茉磥頋M足人們?nèi)找姘l(fā)展的需要。風(fēng)力發(fā)電作為新能源的典型代表，在各種新能源中具有很大的潛力。作為風(fēng)力發(fā)電系統(tǒng)關(guān)鍵技術(shù)的風(fēng)力發(fā)電機(jī)，無刷雙饋電機(jī)的誕生促進(jìn)了風(fēng)力發(fā)電技術(shù)的發(fā)展[1]。

風(fēng)力發(fā)電分為離網(wǎng)風(fēng)力發(fā)電和并網(wǎng)風(fēng)力發(fā)電。離網(wǎng)風(fēng)力發(fā)電發(fā)電是指區(qū)域獨(dú)立風(fēng)力發(fā)電和分戶獨(dú)立風(fēng)力發(fā)電，一般應(yīng)用于遠(yuǎn)離電網(wǎng)的農(nóng)村邊疆地區(qū)，作為獨(dú)立的供電系統(tǒng)。而并網(wǎng)風(fēng)力發(fā)電是將風(fēng)力發(fā)電機(jī)發(fā)出來的電經(jīng)過變換器變換成符合電網(wǎng)要求的電能輸送回電網(wǎng)。因此，并網(wǎng)型風(fēng)力發(fā)電系統(tǒng)具有較廣的發(fā)展前景[2]。本文首先介紹了應(yīng)用于并網(wǎng)風(fēng)力發(fā)電技術(shù)中無刷雙饋電機(jī)的數(shù)學(xué)模型和仿真模型，在此基礎(chǔ)上對并網(wǎng)型無刷雙饋風(fēng)力發(fā)電機(jī)組進(jìn)行了仿真模型的建立，最后在 Matlab/Simulink中對并網(wǎng)型無刷雙饋風(fēng)力發(fā)電機(jī)組分別對啟動過程、風(fēng)速變化、漿距角變化和無刷雙饋電機(jī)參數(shù)變化進(jìn)行了動態(tài)仿真

2 無刷雙饋電機(jī)的數(shù)學(xué)模型

無刷雙饋電機(jī)由于采用比較特殊的定子、轉(zhuǎn)子結(jié)構(gòu)，電機(jī)內(nèi)部的電磁關(guān)系比一般的電機(jī)要復(fù)雜許多，需要針對它的特殊的電磁關(guān)系建立新的數(shù)學(xué)模型。目前無刷雙饋電機(jī)的數(shù)學(xué)模型有網(wǎng)絡(luò)數(shù)學(xué)模型、同步坐標(biāo)數(shù)學(xué)模型、轉(zhuǎn)子速d-q軸模型。

無刷雙饋電機(jī)的網(wǎng)絡(luò)數(shù)學(xué)模型是以定轉(zhuǎn)子的各線圈組的電流作為自變量，線圈組兩端的電壓作為輸入量，用網(wǎng)絡(luò)分析方法建立起來的一種模型。由于在無刷雙饋電機(jī)中存在兩種不同極對數(shù)的定子磁場，在建立數(shù)學(xué)模型時(shí)要采用機(jī)械角進(jìn)行計(jì)算。網(wǎng)絡(luò)數(shù)學(xué)模型能夠較好地反映無刷雙饋電機(jī)中電流的分布，但模型的建立和計(jì)算比較復(fù)雜，不利于分解解耦對無刷雙饋電機(jī)進(jìn)行控制。

同步坐標(biāo)數(shù)學(xué)模型是將功率繞組和控制繞組分別建立在各自的同步坐標(biāo)下，使各輸入輸出變量由時(shí)變量轉(zhuǎn)換為常量，有利于對無刷雙饋電機(jī)的控制。但在風(fēng)力發(fā)電系統(tǒng)中，無刷雙饋發(fā)電機(jī)的控制繞組要與變頻器相連，變換為雙軸模型的控制繞組物理變量不易與變頻器接口直接耦合，給仿真實(shí)驗(yàn)帶來了困難。

轉(zhuǎn)子速d-q軸模型是在同步坐標(biāo)數(shù)學(xué)模型的基礎(chǔ)上建立起來的，將兩相垂直靜止坐標(biāo)系（D，Q）變換成兩相同步旋轉(zhuǎn)坐標(biāo)系（d，q）。無刷雙饋電機(jī)的d-q軸變換有利于控制繞組變量與變頻器銜接和對其進(jìn)行控制。本文的研究就是在這種數(shù)學(xué)模型的基礎(chǔ)上。

無刷雙饋電機(jī)中，轉(zhuǎn)子速坐標(biāo)系與功率繞組同步坐標(biāo)系、控制繞組同步坐標(biāo)系的關(guān)系如圖1所示[3]。

圖1 旋轉(zhuǎn)坐標(biāo)系

圖中上標(biāo)“pe”、“ce”、“r”分別代表功率繞組同步坐標(biāo)系、控制繞組同步坐標(biāo)系、轉(zhuǎn)子速坐標(biāo)系，下標(biāo)“p”、“c”、“r”分別表示功率繞組，控制繞組和轉(zhuǎn)子繞組，下標(biāo)“d”，“q”分別表示d軸和q軸。在本文中，取功率繞組極對數(shù)pp=3，控制繞組極數(shù)pc=1。d-q軸轉(zhuǎn)子速坐標(biāo)系中無刷雙饋電機(jī)的電壓方程為

電磁轉(zhuǎn)矩方程為

式中，rp，Lsp，rc，Lsc，rr，Lr分別是功率繞組和控制繞組的電阻和自感。Mp，Mc是分別是功率繞組與轉(zhuǎn)子，控制繞組與轉(zhuǎn)子的互感[4]。在M atlab/Simulink中建立的無刷雙饋電機(jī)仿真模型如圖2所示。

圖2 基于M atlab/Sim ulink的BDFM仿真模型

其中Subsystem中的無刷雙饋電機(jī)數(shù)學(xué)模型如圖3所示。

圖3 BDFM仿真的數(shù)學(xué)模型

3 并網(wǎng)型無刷雙饋風(fēng)力發(fā)電機(jī)組的建模

并網(wǎng)型無刷雙饋風(fēng)力發(fā)電機(jī)組框圖如圖 4所示，風(fēng)能驅(qū)動風(fēng)力機(jī)轉(zhuǎn)動，風(fēng)力機(jī)將捕獲的風(fēng)能通過傳動機(jī)構(gòu)傳遞給無刷雙饋風(fēng)力發(fā)電機(jī)，無刷雙饋電機(jī)通過矩陣式變換器將發(fā)出來的電能反饋回電

圖4 并網(wǎng)型無刷雙饋風(fēng)力發(fā)電機(jī)組框圖

網(wǎng)。其中描述風(fēng)力機(jī)的數(shù)學(xué)模型為[5]其中，Jr為風(fēng)力機(jī)的轉(zhuǎn)動慣量，ωr為風(fēng)力機(jī)轉(zhuǎn)動角速度，Kr為阻尼系數(shù)，Tα為風(fēng)力機(jī)吸收的動力轉(zhuǎn)矩，Tα的大小可表示為

其中，Pm為風(fēng)力機(jī)所吸收的功率， CP(λ,β)為風(fēng)力機(jī)的功率系數(shù)，λ為葉尖速比且 λ=ωrR / v，β為葉片的漿距角，ρ為空氣密度，R為風(fēng)力機(jī)半徑，v為風(fēng)速。由于 CP(λ,β)與風(fēng)力機(jī)結(jié)構(gòu)相關(guān)，本文采用式（5）描述的特性曲線，該曲線具有一定的普遍性。

無刷雙饋發(fā)電機(jī)在風(fēng)力發(fā)電系統(tǒng)中，可根據(jù)風(fēng)力資源狀況工作于同步、亞同步、超同步三種狀態(tài)。在亞同步時(shí)，控制繞組需從電網(wǎng)吸收能量；而在超同步時(shí)，功率繞組和控制繞組二者都向電網(wǎng)輸送電能，要求變頻器能夠?qū)崿F(xiàn)能量的雙向流動。因此在這里采用矩陣式變換器作為雙向功率變換器。

圖5 矩陣式變換器仿真模型

仿真模型如圖5所示，矩陣模型有3條輸入3條輸出，在每條輸入輸出的交點(diǎn)上有1個(gè)雙向開關(guān)，總共有9個(gè)雙向開關(guān)。矩陣變換器是一次變換器，無需儲能單元，通過控制開關(guān)管的開通與關(guān)斷，即改變輸入項(xiàng)與輸出項(xiàng)的連接方式，實(shí)現(xiàn)變頻輸出[6]。

4 并網(wǎng)型無刷雙饋風(fēng)力發(fā)電機(jī)組的仿真及仿真結(jié)果分析

風(fēng)力系統(tǒng)無刷雙饋發(fā)電機(jī)在雙饋運(yùn)行時(shí)具有同步特性，當(dāng)外界擾動使電機(jī)與同步轉(zhuǎn)速有偏差時(shí)電機(jī)會產(chǎn)生一個(gè)整步轉(zhuǎn)矩將電機(jī)重新牽入同步運(yùn)行狀態(tài)。無刷雙饋電機(jī)的轉(zhuǎn)速是由功率繞組和控制繞組的電源頻率共同決定[7]。功率繞組與電網(wǎng)直接連接，其電壓源為50Hz，380V。由式（7）得控制繞組的電壓源為10Hz，350V。

其中，λopt為最優(yōu)葉尖速比，ng為發(fā)電機(jī)轉(zhuǎn)速，。下面對無刷雙饋風(fēng)力發(fā)電機(jī)組不同條件下的動態(tài)響應(yīng)進(jìn)行了仿真實(shí)驗(yàn)風(fēng)速為8m/s時(shí)，風(fēng)力發(fā)電機(jī)組在零初始狀態(tài)下啟動時(shí)的動態(tài)特性如圖6所示。

圖6 風(fēng)力發(fā)電機(jī)組零初始狀態(tài)下啟動的動態(tài)特性

圖6（a）為無刷雙饋發(fā)電機(jī)的轉(zhuǎn)速從0上升到600rpm，符合穩(wěn)態(tài)運(yùn)行時(shí)轉(zhuǎn)速表達(dá)式圖6（b）為風(fēng)力機(jī)速度從0上升到134r/m in；圖6（c）為捕獲的風(fēng)能。

由于風(fēng)力發(fā)電機(jī)組是非線性系統(tǒng)，為能更直觀表現(xiàn)運(yùn)行條件變化時(shí)無刷雙饋電機(jī)的運(yùn)行特性，采用相對誤差的表達(dá)形式來進(jìn)行仿真。在風(fēng)力發(fā)電機(jī)穩(wěn)態(tài)運(yùn)行的基礎(chǔ)上，風(fēng)速降低1m/s，即風(fēng)速從8m/s降低到7m/s時(shí)，風(fēng)力發(fā)電機(jī)組各變量相對誤差的動態(tài)響應(yīng)如圖7所示。

圖7 風(fēng)速降低1m/s時(shí)風(fēng)力發(fā)電機(jī)組的動態(tài)響應(yīng)曲線

在風(fēng)速不變穩(wěn)態(tài)運(yùn)行的基礎(chǔ)上，槳矩角變化2弧度， 即槳矩角從0弧度增加到2弧度時(shí)，風(fēng)力發(fā)電機(jī)組各變量相對誤差的動態(tài)響應(yīng)如圖8所示。

在風(fēng)力發(fā)電穩(wěn)態(tài)運(yùn)行的基礎(chǔ)上，當(dāng)控制繞組電壓降低5V， 即控制繞組電壓從350V降低到345V時(shí)，風(fēng)力發(fā)電機(jī)組的動態(tài)響應(yīng)見圖9。

圖8 槳矩角變化2弧度時(shí)的風(fēng)力發(fā)電機(jī)組的動態(tài)響應(yīng)曲線

圖9 控制繞組電壓降低5V時(shí)的風(fēng)力發(fā)電機(jī)組的動態(tài)響應(yīng)曲線

從以上仿真圖形中可以得到，經(jīng)過動態(tài)過程后，風(fēng)力發(fā)電機(jī)轉(zhuǎn)速與發(fā)電機(jī)穩(wěn)態(tài)轉(zhuǎn)速均保持為改變前的值，這是因?yàn)闊o刷雙饋風(fēng)力發(fā)電機(jī)工作在雙饋方式下，能夠保持同步特性；當(dāng)風(fēng)速降低、漿距角增加時(shí)，由于輸入能量下降，捕獲的風(fēng)能均減少；當(dāng)控制繞組電壓降低時(shí)，達(dá)到穩(wěn)態(tài)后捕獲的風(fēng)能沒有變化，這是在兩個(gè)穩(wěn)態(tài)之間的過渡過程。

5 結(jié)論

通過對并網(wǎng)型無刷雙饋發(fā)電機(jī)組各個(gè)部分的數(shù)學(xué)分析，建立了并網(wǎng)型無刷雙饋發(fā)電機(jī)的仿真模型，同時(shí)在Matlab/Simulink進(jìn)行了動態(tài)仿真。仿真結(jié)果說明了該仿真模型的正確性，同時(shí)也說明了并網(wǎng)型無刷雙饋發(fā)電機(jī)優(yōu)越的運(yùn)行特性，為以后有關(guān)風(fēng)力發(fā)電的并網(wǎng)運(yùn)行奠定了良好的基礎(chǔ)。

[1]王承煦,張?jiān)?風(fēng)力發(fā)電[M].北京:中國電力出版社,2003.

[2]葉杭冶. 風(fēng)力發(fā)電機(jī)組的控制技術(shù)[M].北京:機(jī)械工業(yè)出版社,2006.

[3]張先勇,吳捷,楊俊華,楊金明. 無刷雙饋風(fēng)力發(fā)電機(jī)滑模功率解耦控制[J].太陽能學(xué),2007,28(1):68-73.

[4]孫國霞,李嘯驄,蔡義明. 大型變速恒頻風(fēng)電系統(tǒng)的建模與仿真[J].電力自動化設(shè)備,2007,27(10)：69-73.

[5]張琦瑋,蔡旭.最大風(fēng)能捕獲風(fēng)力發(fā)電系統(tǒng)及其仿真[J].電機(jī)與控制應(yīng)用,2007,34(5).

[6]黃科元,賀益康,卞松江.矩陣式變換器交流勵(lì)磁的變速恒頻風(fēng)力發(fā)電系統(tǒng)的研究[J].中國電機(jī)工程學(xué)報(bào),2002, 22(11):100-105.

[7]黃守道.無刷雙饋電機(jī)控制方法的研究[D].武漢: 湖南大學(xué), 2005.