• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    考慮不變槳風(fēng)速范圍的風(fēng)電機(jī)組有功功率控制

    2023-02-27 07:38:44陳載宇殷明慧
    電力系統(tǒng)自動(dòng)化 2023年3期
    關(guān)鍵詞:距角變槳風(fēng)輪

    顧 偉,陳載宇,殷明慧,李 群,李 強(qiáng),鄒 云

    (1. 南京理工大學(xué)自動(dòng)化學(xué)院,江蘇省南京市 210094;2. 國(guó)網(wǎng)江蘇省電力有限公司電力科學(xué)研究院,江蘇省南京市 211103)

    0 引言

    風(fēng)電大規(guī)模、高滲透率接入電網(wǎng)將逐步形成新能源為主體的新型電力系統(tǒng)[1-2]。在此背景下,迫切需要變速-變槳風(fēng)電機(jī)組從常規(guī)采用的最大功率點(diǎn)跟蹤(maximum power point tracking,MPPT)控制轉(zhuǎn)變?yōu)楦欕娋W(wǎng)功率指令的有功功率控制(active power control,APC)[3-4],以參與電網(wǎng)的自動(dòng)發(fā)電控制,支撐系統(tǒng)的有功功率平衡。

    現(xiàn)有APC 方法主要有基于閉環(huán)轉(zhuǎn)速控制[3,5-7]和基于預(yù)設(shè)功率給定[8-12]兩種實(shí)現(xiàn)方式,文獻(xiàn)[13]分別將兩者定義為風(fēng)輪主動(dòng)變速和被動(dòng)變速方法。主動(dòng)變速方法在閉環(huán)反饋控制框架下令風(fēng)電機(jī)組跟蹤設(shè)定穩(wěn)定平衡點(diǎn),從而快速運(yùn)行至平衡點(diǎn)處以響應(yīng)功率指令;被動(dòng)變速方法可以利用大慣量風(fēng)輪慣性響應(yīng)來(lái)應(yīng)對(duì)風(fēng)速波動(dòng),具有功率指令響應(yīng)精確、傳動(dòng)鏈載荷小和變槳?jiǎng)幼髁可俚忍攸c(diǎn),更適用于湍流風(fēng)速場(chǎng)景[13],因此也是本文的研究對(duì)象。

    對(duì)于被動(dòng)變速APC 方法,現(xiàn)有研究的關(guān)注點(diǎn)之一是如何利用風(fēng)輪被動(dòng)變速緩沖風(fēng)速波動(dòng),以減輕變槳機(jī)構(gòu)動(dòng)作疲勞。轉(zhuǎn)速調(diào)節(jié)(rotor speed control,RSC)方法[8]最早提出了優(yōu)先利用風(fēng)輪慣性響應(yīng)平抑因風(fēng)速波動(dòng)帶來(lái)的風(fēng)電機(jī)組氣動(dòng)(輸入)功率與電磁(輸出)功率的不平衡。但因僅利用了零度槳距角風(fēng)輪的被動(dòng)變速,RSC 方法會(huì)在風(fēng)速升高或電網(wǎng)功率指令降低時(shí)退化為只依賴變槳調(diào)節(jié)的恒轉(zhuǎn)速控制。為此,文獻(xiàn)[9]提出了集成變速-變槳APC(下文簡(jiǎn)稱IAPC),通過(guò)利用任意槳距角風(fēng)輪被動(dòng)變速應(yīng)對(duì)風(fēng)速波動(dòng),有效減少了變槳?jiǎng)幼?。類似?文獻(xiàn)[14]提出的分段槳距控制也具有利用任意槳距角風(fēng)輪被動(dòng)變速的效果。

    總結(jié)現(xiàn)有被動(dòng)變速APC 方法,槳距角調(diào)節(jié)均發(fā)生在風(fēng)輪轉(zhuǎn)速達(dá)到轉(zhuǎn)速邊界時(shí),均為限轉(zhuǎn)速控制。這使得槳距角設(shè)定與風(fēng)速、風(fēng)電機(jī)組動(dòng)態(tài)和反饋控制器有關(guān),具有很強(qiáng)的隨機(jī)性和不確定性??梢?jiàn),現(xiàn)有研究中的槳距角設(shè)定只關(guān)注限轉(zhuǎn)速調(diào)節(jié),忽視了其對(duì)風(fēng)輪被動(dòng)變速緩沖風(fēng)速波動(dòng)的影響。然而,本文發(fā)現(xiàn)槳距角會(huì)明顯改變被動(dòng)變速風(fēng)輪的動(dòng)態(tài)特性,體現(xiàn)在不同槳距角風(fēng)輪通過(guò)被動(dòng)變速能夠獨(dú)立應(yīng)對(duì)的風(fēng)速波動(dòng)范圍(即不變槳風(fēng)速范圍[15])存在明顯差異。因此,隨機(jī)設(shè)定的槳距角勢(shì)必會(huì)影響風(fēng)輪被動(dòng)變速對(duì)風(fēng)速波動(dòng)的緩沖效果,即風(fēng)電機(jī)組更容易達(dá)到轉(zhuǎn)速邊界而出現(xiàn)風(fēng)輪超速或電磁功率跌落,削弱了APC 的控制性能。

    為此,本文基于槳距角對(duì)不變槳風(fēng)速范圍的影響進(jìn)行分析,以保證風(fēng)電機(jī)組的穩(wěn)定和盡可能利用風(fēng)輪被動(dòng)變速緩沖風(fēng)速波動(dòng)為原則,討論了槳距角與風(fēng)速變化范圍之間的適配關(guān)系。在此基礎(chǔ)上,提出了考慮不變槳風(fēng)速范圍的風(fēng)電機(jī)組APC 方法。該方法通過(guò)增加適配槳距角的設(shè)定環(huán)節(jié),使得被動(dòng)變速風(fēng)輪的不變槳風(fēng)速范圍動(dòng)態(tài)匹配風(fēng)速波動(dòng)范圍,降低了風(fēng)輪達(dá)到轉(zhuǎn)速邊界的頻次,進(jìn)而緩解了風(fēng)輪超速和電磁功率跌落。最后,基于風(fēng)電機(jī)組動(dòng)模實(shí)驗(yàn)平臺(tái),實(shí)驗(yàn)驗(yàn)證了本文方法的有效性。

    1 風(fēng)電機(jī)組模型與APC

    本章主要介紹變速-變槳風(fēng)電機(jī)組的數(shù)學(xué)模型、APC 的基本原理以及基于被動(dòng)變速運(yùn)行模式的風(fēng)電機(jī)組APC。

    1.1 風(fēng)電機(jī)組數(shù)學(xué)模型

    完整的變速-變槳永磁同步風(fēng)力發(fā)電系統(tǒng)包括風(fēng)輪、傳動(dòng)鏈、發(fā)電機(jī)與變流器等部分,如附錄A 圖A1 所示??紤]到電磁動(dòng)態(tài)的響應(yīng)速度要遠(yuǎn)快于機(jī)械動(dòng)態(tài)的響應(yīng)速度[16],本文主要關(guān)注機(jī)械動(dòng)態(tài),忽略電磁動(dòng)態(tài)過(guò)程,認(rèn)為現(xiàn)有電機(jī)控制可以快速且精準(zhǔn)地響應(yīng)電磁轉(zhuǎn)矩參考值。

    由水平軸風(fēng)電機(jī)組的氣動(dòng)理論可知,風(fēng)輪捕獲的氣動(dòng)功率Pm為:

    式中:ρ為空氣密度;R為風(fēng)輪半徑;v為風(fēng)速;CP為風(fēng)能利用系數(shù),是葉尖速比λ與槳距角β的函數(shù),當(dāng)槳距角為零時(shí),CP在最佳葉尖速比λopt處取得最大值CP,max;ω為風(fēng)輪轉(zhuǎn)速。

    風(fēng)電機(jī)組的機(jī)電動(dòng)態(tài)是指含風(fēng)輪、齒輪箱和發(fā)電機(jī)轉(zhuǎn)子的傳動(dòng)軸轉(zhuǎn)速變化過(guò)程。若齒輪箱變速比為ng,機(jī)電動(dòng)態(tài)可表示為如式(3)所示的雙質(zhì)量塊模型[16]。

    式中:Jr和Jg分別為風(fēng)輪與發(fā)電機(jī)轉(zhuǎn)子的轉(zhuǎn)動(dòng)慣量;ωg為發(fā)電機(jī)轉(zhuǎn)速;Tls和Ths分別為低速側(cè)和高速側(cè)轉(zhuǎn)矩;Dr和Dg分別為風(fēng)輪與發(fā)電機(jī)的外阻尼系數(shù);Tm和Tg分別為氣動(dòng)轉(zhuǎn)矩與電磁轉(zhuǎn)矩。

    將低速軸假定為理想的剛性軸,則風(fēng)電機(jī)組機(jī)電動(dòng)態(tài)可簡(jiǎn)化為如式(4)所示的單質(zhì)量塊集總模型[17]。

    1.2 基于被動(dòng)變速的APC

    在風(fēng)電機(jī)組最大出力、最高轉(zhuǎn)速等約束條件[18]下,風(fēng)電場(chǎng)控制將電網(wǎng)有功功率需求分配給場(chǎng)內(nèi)的風(fēng)電機(jī)組。風(fēng)電機(jī)組APC 的控制目標(biāo)是讓風(fēng)電機(jī)組按功率指令Pcmd輸出電磁功率,同時(shí)維持風(fēng)電機(jī)組機(jī)電動(dòng)態(tài)穩(wěn)定。在穩(wěn)態(tài)視角下,風(fēng)電機(jī)組在風(fēng)速v下存在穩(wěn)定平衡點(diǎn)(ω0,β0)(ω0和β0分別為穩(wěn)定平衡點(diǎn)處的風(fēng)輪轉(zhuǎn)速和槳距角),使得捕獲氣動(dòng)功率、輸出電磁功率與電網(wǎng)功率指令三者相等,即

    那么風(fēng)電機(jī)組運(yùn)行在該穩(wěn)定平衡點(diǎn),便實(shí)現(xiàn)了APC的控制目標(biāo)。基于被動(dòng)變速的IAPC 見(jiàn)圖1,一般包括3 個(gè)環(huán)節(jié):變速區(qū)間內(nèi)的風(fēng)輪被動(dòng)變速、限轉(zhuǎn)速的變槳調(diào)節(jié)和低轉(zhuǎn)速段的風(fēng)電機(jī)組穩(wěn)定。

    圖1 基于被動(dòng)變速的IAPCFig.1 IAPC based on passive speed variation

    1)變速區(qū)間內(nèi)的風(fēng)輪被動(dòng)變速。當(dāng)風(fēng)輪轉(zhuǎn)速ω處于變速區(qū)間[ωopt,ωN]內(nèi)時(shí),風(fēng)電機(jī)組槳距角保持不變,電磁功率直接響應(yīng)電網(wǎng)功率指令,發(fā)送至變流器的電磁轉(zhuǎn)矩指令Tg,ref按式(6)給定。

    式中:Te,ref為等效在低速側(cè)的電磁轉(zhuǎn)矩指令。

    此時(shí),風(fēng)輪轉(zhuǎn)速不依靠電磁功率主動(dòng)調(diào)節(jié),而隨傳動(dòng)鏈承受不平衡功率(即Pm-Pe)自由變化。被動(dòng)變速能夠自然利用風(fēng)輪動(dòng)能緩沖風(fēng)電機(jī)組輸入/輸出功率的短時(shí)不平衡,無(wú)須時(shí)刻調(diào)節(jié)電磁功率或槳距角。

    2)限轉(zhuǎn)速的變槳調(diào)節(jié)。當(dāng)風(fēng)輪轉(zhuǎn)速超出變速區(qū)間,即ω>ωN或ω<ωopt時(shí),通過(guò)槳距角調(diào)節(jié)將風(fēng)輪轉(zhuǎn)速限制在變速區(qū)間內(nèi)。具體如圖1 所示,分別在轉(zhuǎn)速上邊界ωN上調(diào)槳距角和轉(zhuǎn)速下邊界ωopt下調(diào)槳距角,并通常根據(jù)實(shí)際轉(zhuǎn)速相對(duì)于轉(zhuǎn)速邊界的偏差,經(jīng)由比例-積分(PI)環(huán)節(jié)確定槳距角指令[8-9,19]。

    3)低轉(zhuǎn)速段的風(fēng)電機(jī)組穩(wěn)定。為避免低風(fēng)速下風(fēng)電機(jī)組失穩(wěn),當(dāng)風(fēng)速降低使風(fēng)輪轉(zhuǎn)速下降至ωopt時(shí),風(fēng)電機(jī)組運(yùn)行切換至MPPT 模式,電磁轉(zhuǎn)矩指令按式(7)給定。

    值得注意的是,雖然是有限轉(zhuǎn)速的變槳調(diào)節(jié),但由于風(fēng)輪的大慣量和變槳機(jī)構(gòu)的動(dòng)作延遲,風(fēng)電機(jī)組在達(dá)到變速區(qū)間邊界時(shí)均容易出現(xiàn)不同程度的轉(zhuǎn)速超調(diào),即ω>ωN或ω<ωopt。由式(7)可知,后者會(huì)造成電磁功率跌落。

    由上述分析可知,現(xiàn)有被動(dòng)變速APC 中的槳距角調(diào)節(jié)大都在風(fēng)輪轉(zhuǎn)速達(dá)到邊界時(shí)啟動(dòng),并基于轉(zhuǎn)速反饋控制律實(shí)現(xiàn)。槳距角設(shè)定與風(fēng)速、風(fēng)電機(jī)組機(jī)電動(dòng)態(tài)和控制器密切相關(guān),具有很強(qiáng)的隨機(jī)性和不確定性。這種只關(guān)注限轉(zhuǎn)速控制的變槳調(diào)節(jié)未考慮到槳距角對(duì)風(fēng)輪被動(dòng)變速緩沖風(fēng)速波動(dòng)效果的影響,忽視了針對(duì)風(fēng)輪被動(dòng)變速過(guò)程的槳距角優(yōu)化設(shè)定。

    2 被動(dòng)變速風(fēng)輪槳距角與風(fēng)速波動(dòng)范圍的適配分析

    針對(duì)上述問(wèn)題,本章首先分析風(fēng)電機(jī)組不變槳風(fēng)速范圍的變化規(guī)律,其次通過(guò)討論實(shí)際風(fēng)速波動(dòng)范圍與不變槳風(fēng)速范圍之間的關(guān)系,確定與風(fēng)速波動(dòng)范圍相匹配的槳距角。

    2.1 不變槳風(fēng)速范圍及其隨槳距角的變化分析

    風(fēng)電機(jī)組被動(dòng)變速利用風(fēng)輪動(dòng)能緩沖輸入/輸出功率之間的不平衡。其中,槳距角是影響氣動(dòng)功率的主要因素,不同槳距角下的風(fēng)輪被動(dòng)變速特性存在差異?;诓蛔儤L(fēng)速范圍[15],本節(jié)首先分析槳距角對(duì)風(fēng)輪被動(dòng)變速的影響。

    不變槳風(fēng)速范圍表征了APC 方法利用風(fēng)輪變速減少變槳?jiǎng)幼鞯男Ч?duì)于功率指令為Pcmd和變速區(qū)間設(shè)為[ωL,ωU]的APC 方法,當(dāng)槳距角為β時(shí),存在一個(gè)最大的風(fēng)速波動(dòng)范圍Unpv,即

    可見(jiàn),只要風(fēng)速波動(dòng)不超出Unpv(β),其可由槳距角為β的風(fēng)輪在變速區(qū)間[ωL,ωU]內(nèi)被動(dòng)變速完全應(yīng)對(duì),無(wú)須啟動(dòng)槳距角調(diào)節(jié),否則,若風(fēng)速超出Unpv(β)(v>vU(β)或v

    由式(9)和式(10)可以看出,Unpv與風(fēng)電機(jī)組參數(shù)(R、CP、ωN)、環(huán)境參數(shù)(ρ)和功率指令相關(guān),但是風(fēng)電機(jī)組參數(shù)和環(huán)境參數(shù)不可調(diào),Unpv主要受功率指令影響。以額定容量為1.5 MW 的風(fēng)電機(jī)組[20]為例,當(dāng)應(yīng)用IAPC 方法且功率指令設(shè)定為300 kW時(shí),不同槳距角下Unpv的上邊界vU和下邊界vL分別如附錄B 圖B2 中虛線和實(shí)線所示。在每一個(gè)槳距角下,vU和vL共同確定了風(fēng)電機(jī)組在該槳距角下僅利用風(fēng)輪被動(dòng)變速(而不啟動(dòng)變槳調(diào)節(jié))可應(yīng)對(duì)的風(fēng)速變化范圍??梢园l(fā)現(xiàn),槳距角變化對(duì)Unpv的影響主要體現(xiàn)在風(fēng)速水平和波動(dòng)范圍兩個(gè)方面。

    1)Unpv內(nèi)的風(fēng)速水平隨槳距角的增大呈現(xiàn)升高趨勢(shì)。當(dāng)槳距角增大時(shí),風(fēng)輪捕獲氣動(dòng)功率的能力下降(即槳距角對(duì)應(yīng)的CP-λ曲線整體下移)。為實(shí)現(xiàn)給定的功率指令目標(biāo),風(fēng)電機(jī)組則需要在更高的風(fēng)速下捕獲足夠的氣動(dòng)功率。這意味著槳距角越大,被動(dòng)變速風(fēng)輪可以應(yīng)對(duì)幅值越高的風(fēng)速波動(dòng)。

    2)Unpv隨槳距角增大呈現(xiàn)擴(kuò)大趨勢(shì)。當(dāng)槳距角增大時(shí),風(fēng)輪捕獲氣動(dòng)功率對(duì)轉(zhuǎn)速變化的敏感度升高,使得相同風(fēng)速擾動(dòng)下的風(fēng)輪轉(zhuǎn)速變化量減小。因此,槳距角越大,相同變速區(qū)間的風(fēng)輪被動(dòng)變速可以緩沖更大范圍的風(fēng)速波動(dòng)。

    3)功率指令越高,相同槳距角下的Unpv上邊界和下邊界均增大,則Unpv內(nèi)的風(fēng)速幅值越大。同時(shí),由于功率指令增大使得風(fēng)電機(jī)組變速區(qū)間縮小,使得Unpv也有所縮小。

    綜上,槳距角會(huì)顯著影響被動(dòng)變速風(fēng)輪的動(dòng)態(tài)特性,主要體現(xiàn)在不同槳距角對(duì)應(yīng)的不變槳風(fēng)速范圍存在明顯差異。因此,面對(duì)實(shí)際湍流風(fēng)速的波動(dòng)范圍變化(即平均風(fēng)速和湍流強(qiáng)度變化),有必要設(shè)定合適的槳距角,從而利用該槳距角風(fēng)輪的被動(dòng)變速來(lái)更有效地緩沖風(fēng)速波動(dòng)。

    2.2 Unp v 與風(fēng)速變化范圍的關(guān)系分析

    在風(fēng)電機(jī)組運(yùn)行過(guò)程中,風(fēng)速波動(dòng)范圍隨機(jī)變化,且未必處于當(dāng)前槳距角對(duì)應(yīng)的不變槳風(fēng)速范圍內(nèi)。兩者之間的關(guān)系直接影響到APC 的控制性能。附錄B 圖B3 所示為一段由Bladed 軟件[21]生成的湍流風(fēng)速序列,可以看出風(fēng)速特征在分鐘級(jí)時(shí)間尺度上就會(huì)存在明顯的差異[22]。

    風(fēng)電機(jī)組在某一槳距角下的Unpv是有限的,很難全面覆蓋長(zhǎng)時(shí)段湍流風(fēng)速波動(dòng)范圍。但將時(shí)間尺度縮短到分鐘級(jí),風(fēng)速波動(dòng)范圍Uv=[vmin,vmax]與Unpv存在可比性。為使風(fēng)輪被動(dòng)變速更有效地緩沖湍流風(fēng)速,緩解風(fēng)電機(jī)組超速和電磁功率跌落,有必要比較分析Unpv與Uv之間的關(guān)系,為合理設(shè)定匹配分鐘級(jí)風(fēng)速波動(dòng)的槳距角奠定基礎(chǔ)。

    圖2 所示為波動(dòng)范圍相同、幅值水平不同的Uv,并將其與Unpv(β)比較。由圖可知,隨著槳距角設(shè)置的不同,Unpv(β)與Uv之間至少存在6 種關(guān)系。

    圖2 Uv 與U np v 之間的6 種關(guān)系Fig.2 Six kinds of relationship between Uv and U np v

    進(jìn)一步地,為了闡述每種關(guān)系意味著的被動(dòng)變速風(fēng)輪的動(dòng)態(tài)過(guò)程及其達(dá)到變速區(qū)間邊界(ωU或ωL)的情況,在風(fēng)輪轉(zhuǎn)速-風(fēng)速平面上,對(duì)應(yīng)固定β(即β1至β6)的風(fēng)電機(jī)組穩(wěn)定平衡點(diǎn)的集合(黑色實(shí)線所示)見(jiàn)附錄B 圖B4(a)至(f)。根據(jù)圖B4 可知,變速區(qū)間[ωL,ωU]確定的風(fēng)速范圍即為Unpv(β)。共用風(fēng)速縱軸的風(fēng)速序列確定了Uv,其中的最高風(fēng)速vmax(或最低風(fēng)速vmin)對(duì)應(yīng)的風(fēng)輪轉(zhuǎn)速若落在變速區(qū)間內(nèi),則是可行穩(wěn)定平衡點(diǎn),否則,便是不可行穩(wěn)定平衡點(diǎn),且意味著風(fēng)電機(jī)組有可能運(yùn)行至變速區(qū)間邊界而啟動(dòng)變槳調(diào)節(jié)。通過(guò)比較Unpv(β)與Uv,它們的6 種關(guān)系描述如下:

    1)關(guān)系1:Unpv(β1)∩Uv=?,且vU(β1)

    2)關(guān) 系 2:vmin

    3)關(guān) 系 3:Unpv(β3)∩Uv=[vmin,vU(β3)],即vL(β3)

    4)關(guān) 系4:Unpv(β4)∩Uv=[vL(β4),vmax],即vmin

    5)關(guān) 系 5:vL(β5)

    6)關(guān)系6:Unpv(β6)∩Uv=?,且vmax

    由上述不變槳風(fēng)速范圍與風(fēng)速變化范圍之間關(guān)系的討論可知,槳距角會(huì)明顯改變被動(dòng)變速風(fēng)輪觸及轉(zhuǎn)速邊界的情況,進(jìn)而影響到風(fēng)電機(jī)組APC 的控制性能。兩類范圍的關(guān)系討論為適配湍流風(fēng)速的槳距角確定奠定了基礎(chǔ)。

    2.3 適配風(fēng)速波動(dòng)范圍的槳距角確定

    基于2.2 節(jié)Unpv(β)與Uv的關(guān)系分析,本文以充分利用風(fēng)輪被動(dòng)變速,盡量避免風(fēng)電機(jī)組超出變速區(qū)間為目標(biāo),得到與風(fēng)速波動(dòng)相適配的槳距角選取原則如下:

    1)首先,應(yīng)滿足關(guān)系5(Uv?Unpv(β)),這樣僅利用風(fēng)輪被動(dòng)變速便可緩沖風(fēng)速波動(dòng),而無(wú)須啟動(dòng)轉(zhuǎn)速邊界處的變槳?jiǎng)幼鳌?/p>

    2)若不存在滿足關(guān)系5 的槳距角,則應(yīng)讓Unpv(β)與Uv的交集最大,且風(fēng)輪轉(zhuǎn)速不觸及ωL。前者可以最大程度地利用被動(dòng)變速,而減少風(fēng)電機(jī)組達(dá)到轉(zhuǎn)速邊界的頻次;后者是為了避免風(fēng)電機(jī)組因電磁功率跌落而不能實(shí)現(xiàn)APC 的控制功能。

    綜合上述原則,對(duì)于給定Uv,本節(jié)利用Unpv(β)與Uv的關(guān)系描述,給出適配槳距角β?確定策略:

    3 考慮Unp v 的APC 方法

    依據(jù)第2 章適配槳距角的確定原則,本章提出考慮Unpv的風(fēng)電機(jī)組APC 方法。該方法通過(guò)主動(dòng)設(shè)定適配槳距角β?,使得被動(dòng)變速風(fēng)輪的Unpv動(dòng)態(tài)匹配湍流風(fēng)速的實(shí)際變化范圍,進(jìn)而提升APC 性能。

    改進(jìn)APC 方法的控制原理框圖如圖3 所示,在IAPC 方法[9]的基礎(chǔ)上增加了基于β?搜索的槳距角設(shè)定環(huán)節(jié)。除了原有的變速區(qū)間內(nèi)風(fēng)輪被動(dòng)變速、限轉(zhuǎn)速的變槳調(diào)節(jié)和低轉(zhuǎn)速段風(fēng)電機(jī)組穩(wěn)定3 個(gè)環(huán)節(jié)之外,新增了基于β?搜索的槳距角設(shè)定環(huán)節(jié)。該槳距角設(shè)定環(huán)節(jié)具體如下。

    圖3 所提APC 方法的原理圖Fig.3 Principle diagram of proposed APC method

    1)Uv估計(jì)模塊。 應(yīng)用基于長(zhǎng)短期記憶(LSTM)網(wǎng)絡(luò)模型的風(fēng)速預(yù)測(cè)方法[24-25](參數(shù)設(shè)定見(jiàn)附錄C 表C1),預(yù)測(cè)得到下一周期的平均風(fēng)速vˉfc,并結(jié)合風(fēng)場(chǎng)的湍流等級(jí)特征估計(jì)風(fēng)速變化范圍Uv。基于該風(fēng)速預(yù)測(cè)方法,本文對(duì)分鐘級(jí)時(shí)間尺度平均風(fēng)速預(yù)測(cè)進(jìn)行仿真。仿真結(jié)果表明,預(yù)測(cè)精度較高,預(yù)測(cè)誤差在10%以內(nèi)(如附錄C 圖C1 所示)。鑒于分鐘級(jí)的Uv與Unpv具有可比性,并對(duì)比不同周期設(shè)定時(shí)的風(fēng)電機(jī)組APC 的控制性能,本文方法將Uv的周期經(jīng)驗(yàn)性地設(shè)定為1 min。

    2)β?搜 索 模 塊。根 據(jù)Uv估 計(jì) 值,并 基 于2.3 節(jié)所述原則和策略,搜索緊鄰下一周期的適配槳距角β?。由于是有限范圍內(nèi)的單變量尋優(yōu),本文在槳距角運(yùn)行范圍[βmin,βmax]內(nèi)采用遍歷方法,搜索滿足式(11)的β?,遍歷搜索步長(zhǎng)為0.1°。特別地,如果存在多個(gè)滿足式(11)的槳距角,則選擇與當(dāng)前槳距角偏差最小的值作為β?。

    需要注意的是,雖然通過(guò)估計(jì)風(fēng)速范圍尋優(yōu)得到的適配槳距角與實(shí)際風(fēng)速對(duì)應(yīng)的適配槳距角會(huì)存在偏差,但是由于大慣量風(fēng)輪具有慢動(dòng)態(tài)特性,可以應(yīng)對(duì)短時(shí)的風(fēng)速大幅變化,對(duì)APC 的控制性能影響較小。

    改進(jìn)方法對(duì)于Uv估計(jì)、β?搜索和設(shè)定是周期性的,即在每個(gè)周期開(kāi)始前,根據(jù)預(yù)測(cè)得到的平均風(fēng)速和風(fēng)場(chǎng)內(nèi)湍流等級(jí)特征信息,估計(jì)出這一周期的Uv,并據(jù)此搜索這一周期的β?,進(jìn)而將槳距角設(shè)定為該β?。需要注意的是,盡管主動(dòng)設(shè)定了適配槳距角,但風(fēng)電機(jī)組依然存在達(dá)到變速區(qū)間邊界的可能。如果發(fā)生,則基于PI 的變槳調(diào)節(jié)會(huì)啟動(dòng),將風(fēng)輪轉(zhuǎn)速限制在變速區(qū)間內(nèi)。

    綜上,由于保留了變速區(qū)間內(nèi)的風(fēng)輪被動(dòng)變速和限轉(zhuǎn)速的變槳調(diào)節(jié)兩個(gè)環(huán)節(jié),本文方法繼承了IAPC 方法的改進(jìn)機(jī)理:分別利用被動(dòng)變速和PI 變槳調(diào)節(jié)應(yīng)對(duì)湍流風(fēng)速中的高頻率、小幅值和低頻率、大幅值的波動(dòng)分量。在此基礎(chǔ)上,新增的槳距角設(shè)定環(huán)節(jié)不僅克服了被動(dòng)變速風(fēng)輪槳距角設(shè)定的隨機(jī)性和不確定性,而且Unpv(β?)與Uv的適配性使得大慣量風(fēng)輪能更好地緩沖高頻率、小幅值的風(fēng)速波動(dòng)分量,減小因觸及轉(zhuǎn)速邊界而引發(fā)的風(fēng)輪超速或電磁功率跌落。

    4 實(shí)驗(yàn)驗(yàn)證與分析

    本章應(yīng)用風(fēng)電機(jī)組動(dòng)模實(shí)驗(yàn)平臺(tái),實(shí)驗(yàn)驗(yàn)證本文提出方法的有效性。首先,介紹風(fēng)電機(jī)組動(dòng)模實(shí)驗(yàn)平臺(tái);其次,介紹APC 性能的評(píng)價(jià)指標(biāo);最后,針對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行比較分析。

    4.1 風(fēng)電機(jī)組動(dòng)模實(shí)驗(yàn)平臺(tái)

    風(fēng) 電 機(jī) 組 動(dòng) 模 實(shí) 驗(yàn) 平 臺(tái)[26]如 附 錄D 圖D1 所示,包括風(fēng)輪模擬器、發(fā)電變流部分和APC 控制器。該實(shí)驗(yàn)平臺(tái)能夠模擬美國(guó)國(guó)家可再生能源實(shí)驗(yàn)室CART 3 風(fēng)電機(jī)組[27](參數(shù)見(jiàn)附錄D 表D1)的氣動(dòng)特性和機(jī)械動(dòng)態(tài)。

    1)風(fēng)輪模擬器由驅(qū)動(dòng)變流器、18.5 kW 三相異步電機(jī)、機(jī)械飛輪和可編程邏輯控制器(PLC)組成。通過(guò)應(yīng)用功率比例縮放、氣動(dòng)-彈性耦合仿真、慣量補(bǔ)償和阻尼補(bǔ)償算法,該實(shí)驗(yàn)平臺(tái)可基于小容量、低轉(zhuǎn)動(dòng)慣量機(jī)組模擬大容量、大轉(zhuǎn)動(dòng)慣量風(fēng)輪的氣動(dòng)特性和慢動(dòng)態(tài)特性[26,28-30]。此外,對(duì)變槳執(zhí)行機(jī)構(gòu)的動(dòng)作特性是通過(guò)一階慣性環(huán)節(jié)和變槳速率限制來(lái)模擬的。

    2)發(fā)電變流部分與實(shí)際風(fēng)電機(jī)組基本相同,包括15 kW 永磁同步發(fā)電機(jī)和并網(wǎng)變流器(含機(jī)側(cè)變流器和網(wǎng)側(cè)變流器)。機(jī)側(cè)變流器接收來(lái)自APC 控制器的電磁轉(zhuǎn)矩指令Tg,ref,控制發(fā)電機(jī)的電磁轉(zhuǎn)矩Tg。由于采用了功率縮放技術(shù)(縮放系數(shù)為50),模擬風(fēng)電機(jī)組的輸出功率是實(shí)驗(yàn)電機(jī)輸出電磁功率的50 倍。

    3)APC 控制器。在PLC 中編程實(shí)現(xiàn)APC 控制策略,包括RSC 方法、IAPC 方法和本文方法??刂破髟诿恳粋€(gè)控制周期計(jì)算電磁轉(zhuǎn)矩指令Tg,ref和槳距角指令βref,并分別發(fā)送至機(jī)側(cè)變流器和風(fēng)輪模擬程序。

    4.2 APC 控制性能的指標(biāo)

    本節(jié)從風(fēng)輪轉(zhuǎn)速調(diào)節(jié)性能、變槳機(jī)構(gòu)疲勞載荷和功率指令響應(yīng)性能3 個(gè)方面比較分析不同APC方法,具體指標(biāo)介紹如下。

    1)為了評(píng)價(jià)風(fēng)輪轉(zhuǎn)速性能,選取風(fēng)輪超速(即轉(zhuǎn)速大于ωN)累積時(shí)間tU和最大轉(zhuǎn)速偏差率ηmaxU來(lái)描述風(fēng)輪超速的嚴(yán)重程度。其中,ωmaxU為統(tǒng)計(jì)時(shí)段內(nèi)的風(fēng)輪最大轉(zhuǎn)速。類似地,根據(jù)實(shí)驗(yàn)數(shù)據(jù)統(tǒng)計(jì)得到轉(zhuǎn)速低于ωL的累積時(shí)間tL,用于描述風(fēng)輪轉(zhuǎn)速越過(guò)ωL而導(dǎo)致電磁功率跌落的嚴(yán)重程度。

    2)選取變槳?jiǎng)幼髁喀う聅um[5,9]來(lái)反映變槳執(zhí)行機(jī)構(gòu)疲勞載荷。其中,β(i)為第i個(gè)采樣時(shí)刻的槳距角,N為統(tǒng)計(jì)時(shí)段內(nèi)采樣次數(shù)。

    3)考慮到APC 的控制目標(biāo),選取電磁功率相對(duì)電網(wǎng)功率指令的偏差統(tǒng)計(jì)σPe指標(biāo)[31]來(lái)綜合反映風(fēng)電機(jī)組能否快速、準(zhǔn)確地跟蹤功率指令。其中,Pe(i)為第i個(gè)采樣時(shí)刻的輸出電磁功率。需要指出的是,電磁功率跌落僅因風(fēng)輪持續(xù)減速越過(guò)ωL而導(dǎo)致,因此tL與σPe具有正相關(guān)性。

    4.3 實(shí)驗(yàn)結(jié)果分析

    在600 s 湍流風(fēng)速(平均風(fēng)速10 m/s,湍流強(qiáng)度C 級(jí))下,針對(duì)不同功率指令情況進(jìn)行實(shí)驗(yàn)驗(yàn)證。圖4 所示為功率指令為150 kW 時(shí)風(fēng)電機(jī)組在RSC 方法、IAPC 方法和本文方法下的轉(zhuǎn)速、槳距角和電磁功率軌跡,表1 所示為相應(yīng)控制性能評(píng)價(jià)指標(biāo)。其他功率指令下的風(fēng)電機(jī)組轉(zhuǎn)速、槳距角和電磁功率軌跡如附錄D 所示。

    圖4 不同APC 方法的比較Fig.4 Comparison of different APC methods

    表1 控制性能評(píng)價(jià)指標(biāo)Table 1 Evaluation indices of control performance

    1)風(fēng)輪轉(zhuǎn)速調(diào)節(jié)性能

    考慮到風(fēng)輪轉(zhuǎn)速調(diào)節(jié)性能的改善是其他APC控制性能提升的根本原因,本文首先針對(duì)該類性能指標(biāo),比較分析3 種APC 方法。

    對(duì)于RSC 方法,因?yàn)閮H運(yùn)用了零度槳距角風(fēng)輪被動(dòng)變速,風(fēng)速增大使得風(fēng)輪很容易加速到ωN,進(jìn)而頻繁觸發(fā)ωN處的限轉(zhuǎn)速變槳調(diào)節(jié)。相應(yīng)地,轉(zhuǎn)速軌跡表現(xiàn)為持續(xù)在ωN附近波動(dòng),并反復(fù)出現(xiàn)風(fēng)輪超速,tU高達(dá)204.92 s。而且,由于變槳機(jī)構(gòu)動(dòng)作緩慢,且電磁功率因響應(yīng)電網(wǎng)指令而無(wú)法參與限轉(zhuǎn)速調(diào)節(jié),風(fēng)輪超速明顯(ηmaxU=11.20%)。此外,同樣因?yàn)榱愣葮嘟怯兄讷@得強(qiáng)勁的氣動(dòng)功率,零度槳距角風(fēng)輪不易減速到ωopt,tL為0。

    由于IAPC 方法應(yīng)用了任意(零度和非零度)槳距角風(fēng)輪被動(dòng)變速,轉(zhuǎn)速不再是RSC 方法下的在ωN附近持續(xù)波動(dòng),而是在變速區(qū)間[ωopt,ωN]大范圍變化,tU大幅減小到76.88 s。但是,限轉(zhuǎn)速變槳調(diào)節(jié)只關(guān)注轉(zhuǎn)速控制效果,通過(guò)其設(shè)定的槳距角不僅不考慮后續(xù)風(fēng)輪被動(dòng)變速對(duì)風(fēng)速波動(dòng)的緩沖效果,而且具有很強(qiáng)的隨機(jī)性和不確定性。例如,觀察圖4 風(fēng)輪轉(zhuǎn)速和槳距角實(shí)驗(yàn)曲線可以看出,t=7 min 內(nèi)風(fēng)輪轉(zhuǎn)速低于ωopt與相鄰前次ωN處限轉(zhuǎn)速變槳調(diào)節(jié)(紫色虛線圓框所示)引發(fā)的槳距角上調(diào)幅度過(guò)大有關(guān)。相較于RSC 方法,IAPC 方法即便更充分利用了風(fēng)輪被動(dòng)變速,但隨機(jī)且與實(shí)際風(fēng)速無(wú)關(guān)的槳距角設(shè)定依然會(huì)使風(fēng)電機(jī)組容易觸及轉(zhuǎn)速邊界,tU與tL之和達(dá)到了120.44 s。

    對(duì)于本文方法,由于在IAPC 方法的基礎(chǔ)上增加了考慮Unpv的槳距角設(shè)定環(huán)節(jié),根據(jù)風(fēng)速信息動(dòng)態(tài)更新的適配槳距角提升了風(fēng)輪被動(dòng)變速對(duì)湍流風(fēng)速的緩沖效果,進(jìn)而有效降低了風(fēng)輪達(dá)到轉(zhuǎn)速邊界的頻次(tU和tL大幅降低至20.64 s 和0 s),緩解了風(fēng)電機(jī)組超速(ηmaxU=4.36%)并避免了電磁功率跌落。

    2)變槳?jiǎng)幼髁?/p>

    因?yàn)橄揶D(zhuǎn)速變槳調(diào)節(jié)是引發(fā)變槳機(jī)構(gòu)動(dòng)作的主要原因,所以風(fēng)輪轉(zhuǎn)速調(diào)節(jié)性能(即風(fēng)輪被動(dòng)變速對(duì)湍流風(fēng)速的緩沖效果)越好,則風(fēng)電機(jī)組越少觸及轉(zhuǎn)速邊界(tU和tL越小),進(jìn)而Δβsum越低。由表1 可見(jiàn),在功率指令為150 kW(0.25 p.u.)時(shí),僅利用零度槳距角風(fēng)輪被動(dòng)變速的RSC 方法,Δβsum高達(dá)173.71°;IAPC 方法應(yīng)用了非零度槳距角風(fēng)輪被動(dòng)變速,但槳距角設(shè)定未適配風(fēng)速變化范圍,Δβsum緩解為89.10°;在此基礎(chǔ)上,本文方法增加適配槳距角的設(shè)定環(huán)節(jié),使得被動(dòng)變速風(fēng)輪能更好地緩沖高頻率、小幅值的風(fēng)速波動(dòng)分量,Δβsum相應(yīng)降低到42.24°,其中還包括了主動(dòng)設(shè)定適配槳距角導(dǎo)致的25.42°。

    3)功率指令響應(yīng)性能

    當(dāng)風(fēng)電機(jī)組在ωopt處出現(xiàn)轉(zhuǎn)速超調(diào)(即ω<ωopt)時(shí),為維持風(fēng)電機(jī)組轉(zhuǎn)速穩(wěn)定,風(fēng)電機(jī)組運(yùn)行切換到MPPT 模式,導(dǎo)致電磁功率出現(xiàn)跌落。風(fēng)輪轉(zhuǎn)速低于ωopt的情況是影響風(fēng)電機(jī)組輸出電磁功率對(duì)功率指令響應(yīng)性能的主要因素。tL越大,則功率指令響應(yīng)性能越差,σPe越大。由表1 可見(jiàn),相較于IAPC 方法,由于確定適配槳距角遵循了風(fēng)輪轉(zhuǎn)速不觸及ωopt的原則,本文方法將σPe降低為0,改善了功率指令響應(yīng)性能。

    雖然本文方法需要估計(jì)風(fēng)速波動(dòng)范圍和遍歷搜索適配槳距角,但由于風(fēng)速波動(dòng)范圍估計(jì)耗時(shí)較少,且適配槳距角的確定是有限范圍內(nèi)的單變量尋優(yōu),計(jì)算負(fù)擔(dān)不重。經(jīng)測(cè)試,單次適配槳距角確定所需平均時(shí)間為18.8 ms,而且風(fēng)電機(jī)組變槳執(zhí)行機(jī)構(gòu)變槳速率較快,為5~10(°)/s[32]。可見(jiàn),本文方法完全滿足適配槳距角的更新、設(shè)定周期(分鐘級(jí)),不會(huì)影響APC 的工程應(yīng)用和性能。

    當(dāng)功率指令為180 kW(0.3 p.u.)時(shí),相較于RSC 和IAPC 方法,本文方法在風(fēng)輪轉(zhuǎn)速調(diào)節(jié)性能、變槳?jiǎng)幼髁糠矫娴目刂菩阅芏加休^為明顯的改善,見(jiàn)附錄D 圖D2 和表D2。與功率指令為150 kW 的場(chǎng)景相似,此處不再贅述APC 的控制性能提升的具體原因。由于t=7 min 內(nèi)的風(fēng)速幅值較低且持續(xù)時(shí)間較長(zhǎng),3 種APC 方法下的風(fēng)電機(jī)組運(yùn)行均切換到MPPT 模式,通過(guò)降低電磁功率避免風(fēng)電機(jī)組持續(xù)減速,以此保證風(fēng)電機(jī)組轉(zhuǎn)速穩(wěn)定。

    進(jìn)一步地,在附錄D 圖D3 所示的風(fēng)速波動(dòng)下,本文方法通過(guò)設(shè)定與風(fēng)速波動(dòng)范圍相適配的風(fēng)輪槳距角,可以完全利用轉(zhuǎn)速區(qū)間內(nèi)的風(fēng)輪被動(dòng)變速來(lái)應(yīng)對(duì)風(fēng)速波動(dòng),風(fēng)輪轉(zhuǎn)速未達(dá)到轉(zhuǎn)速邊界而產(chǎn)生邊界處的恒轉(zhuǎn)速變槳。而RSC 和IAPC 方法下風(fēng)輪轉(zhuǎn)速則多次觸發(fā)轉(zhuǎn)速邊界,導(dǎo)致風(fēng)輪轉(zhuǎn)速超速和電磁功率跌落。

    5 結(jié)語(yǔ)

    對(duì)于被動(dòng)變速APC,變速區(qū)間內(nèi)風(fēng)輪變速對(duì)于湍流風(fēng)速的緩沖效果成為提升控制性能的關(guān)鍵所在。雖然現(xiàn)有APC 方法在變槳控制環(huán)節(jié)采用先進(jìn)控制算法,實(shí)現(xiàn)了對(duì)風(fēng)輪轉(zhuǎn)速超速的抑制,但是其僅關(guān)注限轉(zhuǎn)速控制效果,忽視了槳距角優(yōu)化設(shè)定對(duì)風(fēng)輪被動(dòng)變速過(guò)程的影響。這使得槳距角運(yùn)行值具有很強(qiáng)的隨機(jī)性,并導(dǎo)致風(fēng)輪容易觸及轉(zhuǎn)速邊界而引發(fā)風(fēng)輪超速和電磁功率跌落問(wèn)題。為此,本文分析了不變槳風(fēng)速范圍隨槳距角的變化關(guān)系,建立了槳距角與風(fēng)速變化范圍之間的適配關(guān)系。在此基礎(chǔ)上,本文提出了考慮不變槳風(fēng)速范圍的APC 方法。該方法通過(guò)動(dòng)態(tài)設(shè)定與風(fēng)速波動(dòng)范圍相適配的槳距角,降低了風(fēng)輪達(dá)到轉(zhuǎn)速邊界的頻次,提升了風(fēng)電機(jī)組APC 在風(fēng)輪轉(zhuǎn)速調(diào)節(jié)、變槳機(jī)構(gòu)疲勞載荷和功率指令響應(yīng)3 個(gè)方面的性能。

    需要指出的是,對(duì)于本文方法,當(dāng)風(fēng)輪轉(zhuǎn)速達(dá)到變速區(qū)間邊界時(shí),同樣可以采樣先進(jìn)控制算法來(lái)進(jìn)一步改善風(fēng)輪轉(zhuǎn)速超速程度。

    本文在撰寫(xiě)過(guò)程中獲得國(guó)家自然科學(xué)基金項(xiàng)目(61773214)資助,特此感謝!

    附錄見(jiàn)本刊網(wǎng)絡(luò)版(http://www.aeps-info.com/aeps/ch/index.aspx),掃英文摘要后二維碼可以閱讀網(wǎng)絡(luò)全文。

    猜你喜歡
    距角變槳風(fēng)輪
    基于數(shù)據(jù)驅(qū)動(dòng)的風(fēng)電機(jī)組最優(yōu)槳距角 辨識(shí)方法
    考慮風(fēng)速差異的風(fēng)電場(chǎng)減載方案與一次調(diào)頻策略
    風(fēng)力發(fā)電機(jī)組最優(yōu)槳距角自尋優(yōu)控制算法設(shè)計(jì)
    葉片數(shù)目對(duì)風(fēng)輪位移和應(yīng)力的影響
    從五臟相關(guān)理論淺析祛風(fēng)退翳法在風(fēng)輪疾病的應(yīng)用
    兆瓦級(jí)風(fēng)電機(jī)組變槳距系統(tǒng)設(shè)計(jì)
    變速風(fēng)力發(fā)電機(jī)組變槳距建模與仿真
    基于最小化多變量的獨(dú)立變槳距控制研究
    變槳距風(fēng)力機(jī)在全風(fēng)速段內(nèi)的槳距角控制策略仿真研究
    兆瓦級(jí)風(fēng)電機(jī)組變槳系統(tǒng)的設(shè)計(jì)與應(yīng)用
    少妇人妻一区二区三区视频| 最近最新中文字幕大全电影3| 欧美一区二区国产精品久久精品| 97碰自拍视频| .国产精品久久| 看黄色毛片网站| 国产精品电影一区二区三区| 级片在线观看| 国产精品美女特级片免费视频播放器| 午夜免费激情av| 一个人免费在线观看电影| 97在线视频观看| 国产精品爽爽va在线观看网站| 色综合站精品国产| 国产亚洲91精品色在线| 深夜a级毛片| 欧美日韩在线观看h| 国产成人精品久久久久久| 国产黄片美女视频| 婷婷精品国产亚洲av| 看十八女毛片水多多多| 三级经典国产精品| 变态另类成人亚洲欧美熟女| 日日摸夜夜添夜夜添av毛片| 国产 一区 欧美 日韩| 精品人妻一区二区三区麻豆 | 色吧在线观看| 一区福利在线观看| 亚洲美女黄片视频| 精品久久久久久久久亚洲| 18禁黄网站禁片免费观看直播| 国产视频一区二区在线看| 在线观看免费视频日本深夜| 成人亚洲欧美一区二区av| 在线观看免费视频日本深夜| 免费观看在线日韩| 波野结衣二区三区在线| 日日摸夜夜添夜夜添av毛片| 少妇的逼水好多| 又爽又黄a免费视频| 欧美极品一区二区三区四区| 免费看av在线观看网站| 久久精品夜色国产| 免费大片18禁| 日韩强制内射视频| 亚洲国产欧洲综合997久久,| 久久亚洲精品不卡| 美女大奶头视频| 日产精品乱码卡一卡2卡三| 九九在线视频观看精品| 日韩国内少妇激情av| 精品久久久久久成人av| 亚洲成a人片在线一区二区| 麻豆成人午夜福利视频| 国产黄色小视频在线观看| 国产一区二区亚洲精品在线观看| avwww免费| 国产亚洲精品久久久久久毛片| 欧美激情国产日韩精品一区| 日韩成人伦理影院| 亚洲欧美日韩东京热| 91在线观看av| 亚洲中文字幕一区二区三区有码在线看| av天堂在线播放| 亚洲av.av天堂| 女生性感内裤真人,穿戴方法视频| 91在线精品国自产拍蜜月| 精品欧美国产一区二区三| 精品午夜福利视频在线观看一区| 日韩在线高清观看一区二区三区| 中文资源天堂在线| 亚洲国产高清在线一区二区三| 日本与韩国留学比较| 色吧在线观看| 免费人成视频x8x8入口观看| 精品欧美国产一区二区三| 国产av一区在线观看免费| 一进一出抽搐动态| 久久久国产成人精品二区| 国产成人一区二区在线| 最好的美女福利视频网| 亚州av有码| 日韩成人av中文字幕在线观看 | 国产成人影院久久av| 免费在线观看成人毛片| 成人精品一区二区免费| 女同久久另类99精品国产91| 黄片wwwwww| 色吧在线观看| 一边摸一边抽搐一进一小说| 精品久久久久久成人av| 国产精品人妻久久久久久| 中文亚洲av片在线观看爽| 日本色播在线视频| 听说在线观看完整版免费高清| 国模一区二区三区四区视频| 国产精品久久电影中文字幕| 成人欧美大片| 亚洲中文字幕日韩| 亚洲精品日韩av片在线观看| 国产男人的电影天堂91| 精品久久国产蜜桃| 久久久久九九精品影院| av女优亚洲男人天堂| 国产精品久久久久久久久免| 国产成人一区二区在线| 国产精品福利在线免费观看| 精品久久久久久久久av| 97热精品久久久久久| 在线看三级毛片| 一本精品99久久精品77| 国产v大片淫在线免费观看| 日韩欧美精品免费久久| 搞女人的毛片| 久久久久久久亚洲中文字幕| av在线播放精品| 亚洲美女视频黄频| 亚洲av中文av极速乱| 久久热精品热| 亚洲欧美清纯卡通| 我的老师免费观看完整版| 午夜视频国产福利| 国产乱人偷精品视频| 最近2019中文字幕mv第一页| 欧美国产日韩亚洲一区| 久久久久国产精品人妻aⅴ院| 一进一出抽搐gif免费好疼| 99久久精品热视频| 欧美性感艳星| 一本一本综合久久| 男女那种视频在线观看| 大又大粗又爽又黄少妇毛片口| 成年版毛片免费区| 午夜福利18| 亚洲精品国产成人久久av| 久久久久九九精品影院| 又爽又黄无遮挡网站| 午夜日韩欧美国产| 人人妻人人澡欧美一区二区| 我要看日韩黄色一级片| 麻豆国产av国片精品| 国产精品伦人一区二区| 乱人视频在线观看| 校园人妻丝袜中文字幕| 国产蜜桃级精品一区二区三区| 久久国产乱子免费精品| 久久这里只有精品中国| 精品一区二区免费观看| 亚洲18禁久久av| 亚洲精品粉嫩美女一区| 精品久久久久久久久av| 男女之事视频高清在线观看| 日本精品一区二区三区蜜桃| 色哟哟·www| 欧美性猛交╳xxx乱大交人| 一本久久中文字幕| 好男人在线观看高清免费视频| 精品欧美国产一区二区三| 国产精华一区二区三区| 黄片wwwwww| 天天躁夜夜躁狠狠久久av| 国产 一区精品| 日日撸夜夜添| 女人十人毛片免费观看3o分钟| 男插女下体视频免费在线播放| 国产69精品久久久久777片| 级片在线观看| 美女被艹到高潮喷水动态| 国产白丝娇喘喷水9色精品| 亚洲专区国产一区二区| 中文字幕av成人在线电影| 精品国产三级普通话版| 亚洲18禁久久av| 女人十人毛片免费观看3o分钟| 国产美女午夜福利| 免费黄网站久久成人精品| 成人欧美大片| 久久久久久久久久成人| 在线观看免费视频日本深夜| 免费看a级黄色片| 国内揄拍国产精品人妻在线| 免费av毛片视频| 亚洲熟妇熟女久久| 色av中文字幕| 美女 人体艺术 gogo| 免费在线观看成人毛片| 成年女人毛片免费观看观看9| 丝袜喷水一区| 欧美最黄视频在线播放免费| 真实男女啪啪啪动态图| 联通29元200g的流量卡| 欧美激情国产日韩精品一区| 99视频精品全部免费 在线| 日韩制服骚丝袜av| 久久久久九九精品影院| 日韩欧美国产在线观看| 欧美区成人在线视频| 国产精品av视频在线免费观看| 日日摸夜夜添夜夜爱| 老司机午夜福利在线观看视频| 丝袜喷水一区| 国产伦一二天堂av在线观看| 99久久精品一区二区三区| 赤兔流量卡办理| 赤兔流量卡办理| 三级男女做爰猛烈吃奶摸视频| 国产乱人视频| a级毛片a级免费在线| 嫩草影院精品99| 国产国拍精品亚洲av在线观看| 一个人免费在线观看电影| 久久国内精品自在自线图片| 亚洲一区高清亚洲精品| 久久久久精品国产欧美久久久| 少妇的逼水好多| 亚洲色图av天堂| 久久精品人妻少妇| 在线观看一区二区三区| 我要看日韩黄色一级片| 搡老妇女老女人老熟妇| 乱人视频在线观看| 我要看日韩黄色一级片| 狂野欧美激情性xxxx在线观看| 波多野结衣高清作品| 久久久久久久久大av| 99久久精品热视频| 日日撸夜夜添| 高清日韩中文字幕在线| 色噜噜av男人的天堂激情| 久久国产乱子免费精品| 精品久久久久久久久亚洲| 美女免费视频网站| 人人妻人人澡人人爽人人夜夜 | 欧美激情国产日韩精品一区| 亚洲无线观看免费| 99视频精品全部免费 在线| 热99在线观看视频| 午夜亚洲福利在线播放| 久久精品国产亚洲网站| 少妇的逼水好多| 99久久久亚洲精品蜜臀av| 精品乱码久久久久久99久播| 亚洲五月天丁香| 中出人妻视频一区二区| 中国美白少妇内射xxxbb| 亚洲精品在线观看二区| 女生性感内裤真人,穿戴方法视频| 一本久久中文字幕| 级片在线观看| 精品不卡国产一区二区三区| 久久亚洲精品不卡| 卡戴珊不雅视频在线播放| 亚洲一区高清亚洲精品| 亚洲无线观看免费| 亚洲精品一区av在线观看| 亚洲18禁久久av| 日本欧美国产在线视频| 国产在视频线在精品| 可以在线观看毛片的网站| 亚洲成人久久性| 十八禁网站免费在线| 日韩欧美精品免费久久| 精品人妻一区二区三区麻豆 | 欧美+日韩+精品| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 亚洲va在线va天堂va国产| 99久国产av精品| 18+在线观看网站| 亚洲中文字幕日韩| 在线免费观看的www视频| 99九九线精品视频在线观看视频| 在现免费观看毛片| 女人被狂操c到高潮| 国产精品野战在线观看| 欧美激情久久久久久爽电影| 久久久久免费精品人妻一区二区| 黄色配什么色好看| 欧美又色又爽又黄视频| 亚洲av美国av| 亚洲精品影视一区二区三区av| 最新在线观看一区二区三区| 丝袜喷水一区| 不卡视频在线观看欧美| 久久精品人妻少妇| 自拍偷自拍亚洲精品老妇| 熟女电影av网| av中文乱码字幕在线| 亚洲精品粉嫩美女一区| 亚洲四区av| 直男gayav资源| 久久鲁丝午夜福利片| 亚洲成人久久性| 一级a爱片免费观看的视频| 日韩强制内射视频| 99久国产av精品国产电影| 啦啦啦韩国在线观看视频| 免费看a级黄色片| 精品久久久久久久末码| 亚洲三级黄色毛片| 日韩精品青青久久久久久| 国产一区二区激情短视频| 国产片特级美女逼逼视频| 国内精品久久久久精免费| 寂寞人妻少妇视频99o| 日日撸夜夜添| 日韩欧美精品免费久久| 久久久久久久久久久丰满| 精品福利观看| 成人国产麻豆网| 欧美zozozo另类| 亚洲av电影不卡..在线观看| 欧美激情国产日韩精品一区| 国产三级在线视频| 99九九线精品视频在线观看视频| 免费在线观看成人毛片| 国产精品野战在线观看| 国产高清视频在线播放一区| 哪里可以看免费的av片| 亚洲欧美中文字幕日韩二区| 啦啦啦观看免费观看视频高清| 身体一侧抽搐| 黑人高潮一二区| 亚洲av免费高清在线观看| av在线老鸭窝| 欧美三级亚洲精品| 国产高清视频在线观看网站| 国产精品嫩草影院av在线观看| 亚洲欧美日韩高清在线视频| 午夜福利在线观看免费完整高清在 | 男女之事视频高清在线观看| 啦啦啦啦在线视频资源| 国产精品亚洲美女久久久| 美女大奶头视频| av视频在线观看入口| 国产精品久久久久久久久免| 看免费成人av毛片| 少妇被粗大猛烈的视频| 国产成人a∨麻豆精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲电影在线观看av| 中国国产av一级| 欧美成人精品欧美一级黄| 一级黄色大片毛片| 国产成人福利小说| 国产视频内射| 国产精品1区2区在线观看.| 久久中文看片网| 欧美xxxx性猛交bbbb| 亚洲av五月六月丁香网| 丰满人妻一区二区三区视频av| 亚洲av免费高清在线观看| 可以在线观看毛片的网站| 波多野结衣高清无吗| 午夜亚洲福利在线播放| 亚洲人成网站在线播| 欧美最新免费一区二区三区| 精品乱码久久久久久99久播| 男人和女人高潮做爰伦理| 中文字幕久久专区| 在线观看av片永久免费下载| 九九久久精品国产亚洲av麻豆| 少妇裸体淫交视频免费看高清| 亚洲成a人片在线一区二区| 乱码一卡2卡4卡精品| 国产午夜福利久久久久久| АⅤ资源中文在线天堂| 国产成人aa在线观看| 伦精品一区二区三区| 男女下面进入的视频免费午夜| 欧美成人a在线观看| 别揉我奶头 嗯啊视频| 亚洲成人av在线免费| 十八禁网站免费在线| 非洲黑人性xxxx精品又粗又长| 亚洲无线观看免费| 男女下面进入的视频免费午夜| 亚洲精品国产av成人精品 | 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av中文字字幕乱码综合| 午夜福利在线在线| 亚洲欧美日韩卡通动漫| 国产高清三级在线| 亚洲欧美成人综合另类久久久 | 香蕉av资源在线| 色视频www国产| 精品不卡国产一区二区三区| 噜噜噜噜噜久久久久久91| 欧美国产日韩亚洲一区| av在线蜜桃| 欧美激情在线99| 精品人妻一区二区三区麻豆 | 香蕉av资源在线| 国产人妻一区二区三区在| 啦啦啦观看免费观看视频高清| 日韩av不卡免费在线播放| 日韩欧美免费精品| 国内久久婷婷六月综合欲色啪| 精品一区二区三区人妻视频| 黄色视频,在线免费观看| 色尼玛亚洲综合影院| 亚洲欧美精品自产自拍| 欧美性感艳星| av黄色大香蕉| 91午夜精品亚洲一区二区三区| 最近中文字幕高清免费大全6| 丝袜喷水一区| 欧美又色又爽又黄视频| 亚洲欧美日韩卡通动漫| 天堂网av新在线| 亚洲美女黄片视频| 亚洲欧美清纯卡通| 97超碰精品成人国产| 搡女人真爽免费视频火全软件 | 啦啦啦观看免费观看视频高清| 日产精品乱码卡一卡2卡三| 性插视频无遮挡在线免费观看| 久久精品国产亚洲av天美| 天天躁夜夜躁狠狠久久av| 简卡轻食公司| 日本黄大片高清| 一个人免费在线观看电影| 人人妻人人看人人澡| 国产高清视频在线观看网站| 欧美成人精品欧美一级黄| 12—13女人毛片做爰片一| 俺也久久电影网| 女同久久另类99精品国产91| 亚洲欧美日韩卡通动漫| 精品少妇黑人巨大在线播放 | 在线天堂最新版资源| 女人十人毛片免费观看3o分钟| 高清毛片免费看| 国内少妇人妻偷人精品xxx网站| 亚洲av熟女| 日韩欧美在线乱码| 噜噜噜噜噜久久久久久91| 青春草视频在线免费观看| 国产精品国产三级国产av玫瑰| 晚上一个人看的免费电影| ponron亚洲| 久久6这里有精品| 91在线观看av| 亚洲国产欧洲综合997久久,| 精品久久久久久久久久久久久| 精品人妻视频免费看| 99在线视频只有这里精品首页| 99久久九九国产精品国产免费| 中文字幕熟女人妻在线| 国产精品一区二区免费欧美| 精品国产三级普通话版| 亚洲五月天丁香| 黄色配什么色好看| 欧美一区二区国产精品久久精品| 国产aⅴ精品一区二区三区波| av中文乱码字幕在线| 中文字幕人妻熟人妻熟丝袜美| 蜜臀久久99精品久久宅男| 简卡轻食公司| 免费在线观看影片大全网站| 久久午夜亚洲精品久久| 婷婷精品国产亚洲av在线| 麻豆av噜噜一区二区三区| 午夜a级毛片| 亚洲av成人精品一区久久| 久久久久久大精品| 伦理电影大哥的女人| 久久国内精品自在自线图片| 99国产极品粉嫩在线观看| 午夜精品国产一区二区电影 | 国产黄a三级三级三级人| 免费av观看视频| 91在线精品国自产拍蜜月| 国产亚洲精品久久久久久毛片| 久99久视频精品免费| 国产 一区精品| 中文字幕人妻熟人妻熟丝袜美| 国产精品三级大全| 国内精品一区二区在线观看| 午夜精品在线福利| 亚洲成人精品中文字幕电影| 国内久久婷婷六月综合欲色啪| АⅤ资源中文在线天堂| 国产免费男女视频| 春色校园在线视频观看| 国产男靠女视频免费网站| 别揉我奶头 嗯啊视频| 变态另类丝袜制服| 欧美人与善性xxx| 亚洲精品久久国产高清桃花| 嫩草影视91久久| 久久草成人影院| 久久久久久伊人网av| 波多野结衣巨乳人妻| 99国产精品一区二区蜜桃av| 免费在线观看影片大全网站| 国产v大片淫在线免费观看| 国产精品国产三级国产av玫瑰| 人人妻人人看人人澡| 少妇猛男粗大的猛烈进出视频 | 国产一区二区在线av高清观看| 国产探花极品一区二区| 午夜免费激情av| av天堂在线播放| 亚洲第一电影网av| 露出奶头的视频| 国语自产精品视频在线第100页| 日本欧美国产在线视频| 少妇人妻一区二区三区视频| 乱码一卡2卡4卡精品| 露出奶头的视频| 黄色日韩在线| 非洲黑人性xxxx精品又粗又长| 精品国产三级普通话版| 真人做人爱边吃奶动态| 嫩草影院入口| 国产乱人视频| 欧美国产日韩亚洲一区| 91久久精品国产一区二区三区| 网址你懂的国产日韩在线| 精品一区二区三区视频在线| eeuss影院久久| 淫妇啪啪啪对白视频| 99九九线精品视频在线观看视频| 亚洲人成网站在线播放欧美日韩| 桃色一区二区三区在线观看| 国产精品免费一区二区三区在线| 亚洲性久久影院| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 欧美一级a爱片免费观看看| 国产精品女同一区二区软件| av.在线天堂| 99热6这里只有精品| 日韩欧美一区二区三区在线观看| 97超碰精品成人国产| 97热精品久久久久久| 蜜臀久久99精品久久宅男| 午夜久久久久精精品| 欧美性感艳星| 午夜福利在线观看免费完整高清在 | 特大巨黑吊av在线直播| 国产国拍精品亚洲av在线观看| 久久久久国产精品人妻aⅴ院| 一级毛片我不卡| 1000部很黄的大片| 一级黄片播放器| 老女人水多毛片| 蜜臀久久99精品久久宅男| 久久精品国产99精品国产亚洲性色| 97人妻精品一区二区三区麻豆| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 欧美高清性xxxxhd video| 亚洲av.av天堂| 国产伦精品一区二区三区四那| 色综合站精品国产| 欧美激情国产日韩精品一区| 男女视频在线观看网站免费| 成人特级黄色片久久久久久久| 1024手机看黄色片| 男女啪啪激烈高潮av片| 欧美中文日本在线观看视频| 日韩精品有码人妻一区| 成人美女网站在线观看视频| 精品久久国产蜜桃| 搡老熟女国产l中国老女人| 女人十人毛片免费观看3o分钟| 欧美日韩综合久久久久久| 1024手机看黄色片| 最好的美女福利视频网| 国产v大片淫在线免费观看| 国产精品久久久久久精品电影| 精品久久久久久久久久免费视频| 国产黄a三级三级三级人| 国产伦在线观看视频一区| 一个人观看的视频www高清免费观看| 成年女人看的毛片在线观看| 国产真实乱freesex| 亚洲乱码一区二区免费版| 蜜桃久久精品国产亚洲av| 香蕉av资源在线| 老司机午夜福利在线观看视频| 亚洲自拍偷在线| 啦啦啦啦在线视频资源| 热99在线观看视频| 久久久精品大字幕| 最近视频中文字幕2019在线8| 91狼人影院| 国产黄色小视频在线观看| 九九在线视频观看精品| 国产人妻一区二区三区在| 国产精品福利在线免费观看| 老司机午夜福利在线观看视频| 丰满的人妻完整版| 国内精品久久久久精免费| 黄色视频,在线免费观看| 少妇熟女aⅴ在线视频| 97热精品久久久久久| 一进一出抽搐gif免费好疼| 亚洲一区高清亚洲精品| 久久99热这里只有精品18| 黄色欧美视频在线观看| 九色成人免费人妻av| 插逼视频在线观看| 亚洲欧美日韩高清专用| 91午夜精品亚洲一区二区三区| 欧美成人免费av一区二区三区| 久久久久久伊人网av| 看片在线看免费视频| 欧美日韩乱码在线| 亚洲专区国产一区二区| 欧美xxxx性猛交bbbb| 国产高清激情床上av| 可以在线观看毛片的网站|