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

    基于頻率安全約束與臨界慣量計(jì)算的分時段限制風(fēng)電出力方法

    2022-08-19 06:34:06李世春田冰杰李惠子黃森焰徐松林
    電力系統(tǒng)保護(hù)與控制 2022年15期
    關(guān)鍵詞:慣量電功率時段

    李世春,田冰杰,李惠子,羅 穎,黃森焰,徐松林

    基于頻率安全約束與臨界慣量計(jì)算的分時段限制風(fēng)電出力方法

    李世春,田冰杰,李惠子,羅 穎,黃森焰,徐松林

    (三峽大學(xué)電氣與新能源學(xué)院,湖北 宜昌 443002)

    高比例風(fēng)電接入導(dǎo)致系統(tǒng)等效慣量降低,增加了電網(wǎng)引發(fā)頻率失穩(wěn)的可能性。為保障日內(nèi)頻率安全穩(wěn)定運(yùn)行,同時避免大面積切除風(fēng)電的不合理操作,提出了基于頻率安全約束與臨界慣量計(jì)算的分時段限制風(fēng)電出力方法。通過引入頻率最低點(diǎn)和最大頻率變化率約束指標(biāo),并考慮日內(nèi)運(yùn)行方式變化,分時段求解電網(wǎng)臨界慣量及對應(yīng)的最大風(fēng)電并網(wǎng)容量。結(jié)合電網(wǎng)實(shí)際慣量判定小于臨界慣量的時段,對這些時段依據(jù)最大風(fēng)電并網(wǎng)容量確定風(fēng)電切除量,從而在日內(nèi)各時段消除頻率安全隱患。最后通過算例系統(tǒng)證明:所提方法能準(zhǔn)確篩查電網(wǎng)日內(nèi)存在頻率安全風(fēng)險的時段,求解的風(fēng)功率分時段切除量可有效消除頻率安全隱患。

    高比例風(fēng)電;頻率安全約束;臨界慣量;最大風(fēng)電并網(wǎng)容量;風(fēng)電切除量

    0 引言

    在國家碳達(dá)峰、碳中和的目標(biāo)引領(lǐng)下,風(fēng)電大規(guī)模、高比例接入電網(wǎng)[1-3],將會導(dǎo)致電力系統(tǒng)等效慣量持續(xù)降低、頻率穩(wěn)定運(yùn)行能力顯著削弱[4-6]。目前,在應(yīng)對大規(guī)模風(fēng)電并網(wǎng)引起的潛在電網(wǎng)安全穩(wěn)定問題時,主要采取限制風(fēng)電出力的方法[7-9],通過對風(fēng)電場直接解列的“一刀切”方式來實(shí)施,往往會造成大量棄風(fēng),犧牲不必要的發(fā)電效益。鑒于此,當(dāng)面臨頻率安全威脅時,制定合理的風(fēng)電切除方案、通過準(zhǔn)確計(jì)算頻率安全約束下的風(fēng)電最大并網(wǎng)容量來確定切除量,對于科學(xué)指導(dǎo)風(fēng)電并網(wǎng)、保障電網(wǎng)頻率安全具有重要意義。

    關(guān)于風(fēng)電最大并網(wǎng)容量的計(jì)算問題,一些學(xué)者從不同角度開展過相關(guān)研究。文獻(xiàn)[10-11]考慮風(fēng)速的隨機(jī)性,應(yīng)用機(jī)會約束規(guī)劃理論,在滿足系統(tǒng)潮流約束下,將風(fēng)電接入容量最大化作為優(yōu)化目標(biāo)來求解風(fēng)電穿透功率極限。文獻(xiàn)[12]在考慮風(fēng)電時空分布特征和功率平衡條件下,通過引入切負(fù)荷量和棄風(fēng)量等風(fēng)險指標(biāo),建立受限風(fēng)險閾值下的多變量、非線性優(yōu)化模型,求解了風(fēng)電準(zhǔn)入容量。上述研究從規(guī)劃或穩(wěn)態(tài)優(yōu)化運(yùn)行角度研究風(fēng)電最大并網(wǎng)容量,未涉及電網(wǎng)安全穩(wěn)定約束因素。隨著風(fēng)電滲透率不斷提高,電網(wǎng)頻率安全穩(wěn)定水平已經(jīng)成為限制風(fēng)電并網(wǎng)容量的主要因素之一。

    近年來,也有學(xué)者基于頻率穩(wěn)定約束研究風(fēng)電最大并網(wǎng)容量計(jì)算問題。文獻(xiàn)[13]基于風(fēng)電參與調(diào)頻的前提條件,考慮穩(wěn)態(tài)頻率偏差和頻率變化率約束求解風(fēng)電穿透功率極限,但該研究方法不能適應(yīng)電網(wǎng)運(yùn)行方式變化下的計(jì)算和應(yīng)用。文獻(xiàn)[14]通過對澳大利亞電網(wǎng)設(shè)置不同的風(fēng)電滲透率,進(jìn)行多次動態(tài)仿真試驗(yàn),選擇滿足頻率變化率和頻率最低點(diǎn)約束的場景,以此求得系統(tǒng)的風(fēng)電穿透功率極限,但該文獻(xiàn)缺乏理論計(jì)算過程。文獻(xiàn)[15]在分析不同風(fēng)電滲透率對電力系統(tǒng)的等效慣量和發(fā)電機(jī)單位調(diào)節(jié)功率影響的基礎(chǔ)上,計(jì)算評估了考慮最大頻率跌落值和穩(wěn)態(tài)頻率偏差穩(wěn)定約束的風(fēng)電并網(wǎng)接納能力。文獻(xiàn)[16]建立了包含常規(guī)電源和可再生能源的電網(wǎng)頻率響應(yīng)綜合模型,通過設(shè)置不同風(fēng)電滲透率,分析擾動后頻率變化率和最大頻率偏移量,直至滿足頻率安全約束條件,最終確定可再生能源可接納容量。上述研究通過計(jì)及頻率穩(wěn)定約束來求解風(fēng)電最大并網(wǎng)容量,其計(jì)算過程和結(jié)果是粗略的。實(shí)際上,在1天24 h的各時段,電網(wǎng)慣量具有時變特征,風(fēng)電等可再生能源出力也具有波動性,因此日內(nèi)各時段的頻率安全裕度存在差異,需要按時段確定最大風(fēng)電并網(wǎng)容量,以精細(xì)化指導(dǎo)電網(wǎng)日常安全穩(wěn)定運(yùn)行。

    綜上所述,現(xiàn)有文獻(xiàn)從規(guī)劃或穩(wěn)態(tài)優(yōu)化運(yùn)行角度研究風(fēng)電最大并網(wǎng)容量,一部分文獻(xiàn)雖然也考慮頻率穩(wěn)定約束求解風(fēng)電最大并網(wǎng)容量,但存在理論計(jì)算過程缺乏和計(jì)算太粗略的問題,無法科學(xué)細(xì)致地指導(dǎo)電網(wǎng)日常頻率安全穩(wěn)定運(yùn)行的具體操作。鑒于此,本文提出基于頻率安全約束與臨界慣量計(jì)算的分時段限制風(fēng)電出力方法,以下詳細(xì)論述其理論計(jì)算過程。

    1 基于頻率安全約束的風(fēng)電最大并網(wǎng)容量計(jì)算

    1.1 頻率安全約束指標(biāo)

    1.2 計(jì)算風(fēng)電最大并網(wǎng)容量的基本原理

    若將電力系統(tǒng)等效為一個整體,其等值轉(zhuǎn)子運(yùn)動方程[17]可表示為

    根據(jù)式(2)可求得最大頻率偏差和頻率變化率。

    若電網(wǎng)包含臺同步機(jī)組和個風(fēng)電場,日內(nèi)時段的等效慣量[13]為

    2 臨界慣量和電網(wǎng)實(shí)際慣量計(jì)算

    在式(7)中,若要求解WFmax,須先計(jì)算得到電網(wǎng)臨界慣量和低于臨界慣量的時段。因此,電網(wǎng)分時段實(shí)際慣量和臨界慣量計(jì)算是本文研究的核心問題,以下詳細(xì)論述其過程。

    2.1 電網(wǎng)慣量計(jì)算

    由式(5)可知,計(jì)算次日電網(wǎng)分時段實(shí)際慣量,需知次日同步機(jī)組啟停狀況和風(fēng)電并網(wǎng)容量。從調(diào)度中心獲取次日發(fā)電計(jì)劃,可獲得機(jī)組組合及其對應(yīng)的慣性時間常數(shù)。短期風(fēng)電功率預(yù)測可以預(yù)測次日零時起3天內(nèi)的風(fēng)電輸出功率[18],預(yù)測點(diǎn)時間分辨率為15 min。短期風(fēng)電功率負(fù)荷預(yù)測技術(shù)在電網(wǎng)調(diào)度應(yīng)用方面比較成熟,且可用超短期風(fēng)電功率預(yù)測值、實(shí)時風(fēng)電功率滾動修正風(fēng)電功率出力曲線。本文默認(rèn)已得到風(fēng)電功率出力曲線。

    2.2 基于頻率安全指標(biāo)求解臨界慣量

    2.2.1頻率最低點(diǎn)約束的臨界慣量

    可采用平均系統(tǒng)頻率(Average System Frequency, ASF)響應(yīng)模型[19-21]來描述包含多臺同步機(jī)組的電力系統(tǒng)頻率響應(yīng)特性,并可將各機(jī)組調(diào)速系統(tǒng)模型采用一階慣性環(huán)節(jié)來近似反映其響應(yīng)過程[22],如圖1所示。

    圖1 多機(jī)系統(tǒng)的頻率響應(yīng)閉環(huán)傳遞函數(shù)模型

    圖1中同步機(jī)組一次調(diào)頻響應(yīng)的機(jī)械功率增量可表示為

    由式(9)可見,系統(tǒng)頻率偏差由初始功率缺額和系統(tǒng)慣量決定,并與時間t呈線性比例關(guān)系,在該時間段內(nèi)、易于求解。另一方面,頻率跌落最低點(diǎn)通常在功率缺額擾動后3~5 s到達(dá),此刻一次調(diào)頻響應(yīng)即開始動作。因此,可將頻率跌落最低點(diǎn)到達(dá)之前、經(jīng)式(9)計(jì)算的頻率偏差作為計(jì)算頻率跌落最低點(diǎn)時刻的同步機(jī)組一次調(diào)頻機(jī)械功率增量的輸入,這實(shí)際上是將圖1的閉環(huán)傳遞函數(shù)模型打開為如圖2所示開環(huán)模型,可簡化求解。

    根據(jù)圖2所示的傳遞函數(shù)模型,將式(9)的時域表達(dá)式代入式(8),并進(jìn)行拉式反變換后可得

    在研究電網(wǎng)頻率降到最低點(diǎn)的過程中,可對同步機(jī)組的機(jī)械功率增量進(jìn)行線性化處理,用常系數(shù)一次函數(shù)近似表達(dá)其機(jī)械功率增量。對于典型系統(tǒng),C可取頻率最低點(diǎn)時刻對應(yīng)的值[22]。

    設(shè)系統(tǒng)頻率最低點(diǎn)時刻為min,則同步機(jī)組對應(yīng)的機(jī)械功率增量為

    根據(jù)上述分析并結(jié)合圖2可知,受擾后整個系統(tǒng)的動態(tài)頻率方程為

    式中,表示同步機(jī)組的臺數(shù)。

    對式(12)兩端積分并整理后可得

    將式(14)代入式(13),計(jì)算頻率最低點(diǎn)對應(yīng)的最大角頻率偏差為

    再轉(zhuǎn)換為最大系統(tǒng)頻率偏差的有名值,如式(16)所示。

    2.2.2頻率變化率約束的臨界慣量

    系統(tǒng)最大頻率變化率出現(xiàn)在擾動發(fā)生后瞬間,同步機(jī)組機(jī)械功率增量和等值阻尼系數(shù)都可忽略,根據(jù)式(3),最大頻率變化率[23]可等效為

    由前述分析可知,電網(wǎng)臨界慣量應(yīng)取式(18)和式(20)中較大值。

    3 風(fēng)電最大并網(wǎng)容量計(jì)算

    將電力系統(tǒng)各時段電網(wǎng)實(shí)際慣量H與其對應(yīng)臨界慣量進(jìn)行比較。如果H<max,表明此時段的電網(wǎng)慣量不滿足臨界慣量要求,需要切除部分風(fēng)電出力,以消除頻率安全隱患;如果H≥max,則電網(wǎng)慣量滿足臨界慣量要求,電網(wǎng)無頻率安全隱患,基于文中前述理論分析過程,風(fēng)電最大并網(wǎng)容量計(jì)算流程圖如附圖A1。

    系統(tǒng)大功率缺額分為負(fù)荷突增和同步機(jī)組脫網(wǎng)兩類,以下分別探討兩者對應(yīng)的風(fēng)電最大并網(wǎng)容量。

    3.1 負(fù)荷突增下的風(fēng)電最大并網(wǎng)容量計(jì)算

    根據(jù)上述公式計(jì)算各時段的電網(wǎng)慣量H和電網(wǎng)臨界慣量1max。對于電網(wǎng)慣量低于臨界慣量時段,令切除風(fēng)電后的電網(wǎng)慣量恰等于此時的臨界慣量,根據(jù)式(7)可求得風(fēng)電最大并網(wǎng)容量為

    3.2 同步機(jī)組脫網(wǎng)下的風(fēng)電最大并網(wǎng)容量計(jì)算

    同步機(jī)組脫網(wǎng)會造成電網(wǎng)損失掉一部分慣量,從而呈現(xiàn)和負(fù)荷突增擾動不同的頻率響應(yīng)過程。因此,此類型故障下臨界慣量應(yīng)有所提高,保證電網(wǎng)在損失一部分慣量后仍滿足動態(tài)頻率安全約束。計(jì)算得到各時段的電網(wǎng)慣量H和同步機(jī)組脫網(wǎng)故障下的電網(wǎng)臨界慣量2max。對于電網(wǎng)慣量低于臨界慣量時段,令切除風(fēng)電后的電網(wǎng)慣量恰等于臨界慣量,根據(jù)式(7)可求得風(fēng)電最大并網(wǎng)容量為

    4 算例分析與驗(yàn)證

    采用如圖4所示IEEE-39算例系統(tǒng),在Matlab/ Simulink下驗(yàn)證計(jì)算方法和結(jié)果的正確性、有效性。同步機(jī)組參數(shù)(發(fā)電機(jī)慣性時間常數(shù)、功頻特性系數(shù)、調(diào)速器響應(yīng)時間常數(shù))如附表A1所示,同步機(jī)組各時段發(fā)電計(jì)劃如附錄表A2所示。

    圖4 算例系統(tǒng)

    圖5 各時段風(fēng)電功率值

    以下分別針對負(fù)荷突增和同步機(jī)組脫網(wǎng)兩種典型功率缺額狀況,求解對應(yīng)的臨界電網(wǎng)慣量和風(fēng)電最大并網(wǎng)容量,并進(jìn)行仿真計(jì)算與分析。

    4.1 負(fù)荷突增時的計(jì)算分析

    考慮較嚴(yán)重的實(shí)際情況,設(shè)置突增負(fù)荷占總負(fù)荷的15%。根據(jù)附表A2的同步機(jī)組各時段發(fā)電計(jì)劃和風(fēng)電功率預(yù)測結(jié)果,利用式(5)計(jì)算得到算例系統(tǒng)1天24 h的電網(wǎng)慣量,并通過式(21)求解電網(wǎng)臨界慣量,如圖6所示。

    圖6 負(fù)荷突增時的電網(wǎng)慣量與臨界慣量

    由圖6可知,電網(wǎng)慣量和臨界慣量均具有時變特征,使得一天中某些時段電網(wǎng)慣量小于臨界慣量,如01:00—02:00、18:00—22:00時段。因此,需對這些時段限制風(fēng)電功率,以消除頻率安全隱患。

    可利用式(22)求得上述時段對應(yīng)的風(fēng)電最大并網(wǎng)容量及其切除量,具體結(jié)果如附圖A2和附表A3所示。在不同時段需切除對應(yīng)的風(fēng)電功率,例如在19:00—21:00時段需切除較大風(fēng)電功率,來滿足電網(wǎng)臨界慣量要求。

    為驗(yàn)證上述風(fēng)電最大并網(wǎng)容量及其切除量的正確性和準(zhǔn)確性,分別在01:00—01:15和19:00—19:15時段150 s處設(shè)置突增15%額定負(fù)荷,兩個時段風(fēng)速分別為12 m/s和10 m/s。對比仿真切除/不切除風(fēng)電的系統(tǒng)頻率受擾曲線,具體如圖7—圖10所示。

    在01:00—01:15時段發(fā)生負(fù)荷突增功率缺額時,由圖7、圖8可知:1) 當(dāng)不對風(fēng)電出力限制時(不切除),頻率跌落最低點(diǎn)達(dá)到48.86 Hz,低于最低頻率跌落安全值49 Hz。通過限制風(fēng)電出力,切除408 MW風(fēng)電功率,可使頻率跌落最低點(diǎn)恰好滿足49 Hz的頻率安全值要求;2) 無論是否限制風(fēng)電功率,系統(tǒng)頻率變化率均可滿足-0.5 Hz/s安全值要求,這說明在該狀況下,通過式(21)計(jì)算臨界慣量和利用式(22)計(jì)算最大風(fēng)電并網(wǎng)容量時,系統(tǒng)頻率跌落最低點(diǎn)約束指標(biāo)是主導(dǎo)因素。

    圖7 01:00—01:15時段負(fù)荷突增的頻率受擾曲線

    圖8 01:00—01:15時段負(fù)荷突增的頻率變化率曲線

    圖9 19:00—19:15時段負(fù)荷突增的頻率受擾曲線

    圖10 19:00—19:15時段負(fù)荷突增的頻率變化率曲線

    在19:00—19:15時段發(fā)生負(fù)荷突增功率缺額時,由圖9、圖10可知:1) 當(dāng)不對風(fēng)電出力限制時(不切除),頻率變化率最大值達(dá)到-0.625 Hz/s,超過了頻率變化率安全值-0.5 Hz/s。通過限制風(fēng)電功率,切除489 MW風(fēng)電功率,可使頻率變化率最大值滿足-0.5 Hz/s安全值要求。2) 另一方面,當(dāng)不限制風(fēng)電功率時,系統(tǒng)頻率跌落最低點(diǎn)低于最低頻率跌落安全值49 Hz。當(dāng)限制功率時(切除489 MW),系統(tǒng)頻率最低點(diǎn)顯著高于49 Hz。這說明在該狀況下,通過式(21)計(jì)算臨界慣量和利用式(22)計(jì)算最大風(fēng)電并網(wǎng)容量時,系統(tǒng)頻率變化率約束指標(biāo)是主導(dǎo)因素,為了保證該約束指標(biāo),切除了489 MW風(fēng)電功率,而實(shí)際上該切除量足以保證系統(tǒng)頻率跌落最低點(diǎn)安全約束。

    上述仿真計(jì)算與分析表明,在突增負(fù)荷時,本文對于風(fēng)電最大并網(wǎng)容量及其切除量的計(jì)算結(jié)果能非常準(zhǔn)確地控制在最大頻率跌落最低點(diǎn)約束指標(biāo)值或最大頻率變化率約束指標(biāo)值范圍內(nèi),證明了本文計(jì)算方法的正確性和計(jì)算結(jié)果的準(zhǔn)確性。

    4.2 同步機(jī)組脫網(wǎng)時的計(jì)算分析

    在此設(shè)置算例系統(tǒng)中機(jī)組G9脫網(wǎng),由于同步機(jī)組脫網(wǎng)會影響到系統(tǒng)慣量,因此電網(wǎng)慣量和臨界慣量計(jì)算結(jié)果與負(fù)荷突增的情況存在一定差異。此時,同樣根據(jù)附表A2的同步機(jī)組各時段發(fā)電計(jì)劃和風(fēng)電功率預(yù)測結(jié)果,利用式(5)計(jì)算得到算例系統(tǒng)1天24 h的電網(wǎng)慣量,并通過式(21)求解臨界慣量,如圖11所示。

    圖11 同步機(jī)組脫網(wǎng)時的電網(wǎng)慣量與臨界慣量

    由圖11可知,電網(wǎng)慣量小于臨界慣量的時段主要在01:00—05:00。同樣地,需對這些時段限制風(fēng)電功率,以消除頻率安全隱患。

    可利用式(23)求得上述時段對應(yīng)的風(fēng)電最大并網(wǎng)容量及其切除量,具體結(jié)果如附圖A3和附表A4所示。在01:00—05:00時段切除一定量的風(fēng)電功率,即可滿足臨界慣量的要求。

    為驗(yàn)證上述風(fēng)電最大并網(wǎng)容量及其切除量的正確性和準(zhǔn)確性,考慮算例系統(tǒng)中較大容量機(jī)組脫網(wǎng)情況,分別在01:00—01:15和02:00—02:15時段150 s處設(shè)置G9脫網(wǎng),兩個時段風(fēng)速分別為12 m/s和10.8 m/s。對比仿真切除/不切除風(fēng)電的系統(tǒng)頻率受擾曲線,具體如圖12—圖15所示。

    圖12 01:00—01:15時段機(jī)組脫網(wǎng)的頻率受擾曲線

    圖13 01:00—01:15時段機(jī)組脫網(wǎng)的頻率變化率曲線

    圖14 02:00—02:15時段機(jī)組脫網(wǎng)的頻率受擾曲線

    圖15 02:00—02:15時段機(jī)組脫網(wǎng)的頻率變化率曲線

    在01:00—01:15時段發(fā)生機(jī)組G9脫網(wǎng)時,由圖12、圖13可知:1) 當(dāng)不對風(fēng)電出力限制時(不切除),頻率變化率最大值達(dá)到-0.58 Hz/s,超過了頻率變化率安全值-0.5 Hz/s。通過限制風(fēng)電出力,切除363 MW風(fēng)電功率,可使頻率變化率最大值滿足-0.5 Hz/s安全值要求。2) 另一方面,當(dāng)不限制風(fēng)電功率時,系統(tǒng)頻率跌落最低點(diǎn)超過最低頻率跌落安全值49 Hz。當(dāng)限制出力時(切除363 MW),系統(tǒng)頻率最低點(diǎn)顯著高于49 Hz。這說明在該狀況下,通過式(21)計(jì)算臨界慣量和利用式(23)計(jì)算最大風(fēng)電并網(wǎng)容量時,系統(tǒng)頻率變化率約束指標(biāo)是主導(dǎo)因素,為了保證該約束指標(biāo),切除了363 MW風(fēng)電功率,而實(shí)際上該切除量足以保證系統(tǒng)頻率跌落最低點(diǎn)安全約束。

    在02:00—02:15時段發(fā)生機(jī)組G9脫網(wǎng)時,由圖14、圖15可知:1) 當(dāng)不對風(fēng)電功率限制時(不切除),頻率跌落最低點(diǎn)達(dá)到48.93 Hz,低于最低頻率跌落安全值49 Hz。通過限制風(fēng)電出力,切除335 MW風(fēng)電功率,可使頻率跌落最低點(diǎn)恰好滿足49 Hz安全值要求。2) 無論是否限制風(fēng)電功率,系統(tǒng)頻率變化率均可滿足-0.5 Hz/s的安全值要求,這說明在該狀況下,通過式(21)計(jì)算臨界慣量和利用式(23)計(jì)算最大風(fēng)電并網(wǎng)容量時,系統(tǒng)頻率跌落最低點(diǎn)約束指標(biāo)是主導(dǎo)因素。

    上述仿真計(jì)算與分析表明,在同步機(jī)組脫網(wǎng)時,本文對于風(fēng)電最大并網(wǎng)容量及其切除量的計(jì)算結(jié)果能準(zhǔn)確地控制在最大頻率跌落最低點(diǎn)約束指標(biāo)值或最大頻率變化率約束指標(biāo)值范圍內(nèi),證明了本文計(jì)算方法的正確性和計(jì)算結(jié)果的準(zhǔn)確性。

    綜上可知,在負(fù)荷突增和機(jī)組脫網(wǎng)兩種典型功率缺額狀況下,切除一定量的風(fēng)電后,可保障頻率變化率和頻率最低點(diǎn)均在安全值內(nèi),證明了本文所提風(fēng)電最大并網(wǎng)容量及其切除量計(jì)算方法的正確性和計(jì)算結(jié)果的準(zhǔn)確性。

    4.3 風(fēng)電功率預(yù)測誤差對本文風(fēng)電切除方案影響分析

    對風(fēng)電出力實(shí)際值和經(jīng)超短期風(fēng)電功率滾動修正后平均絕對誤差為5%的預(yù)測值[26-29]進(jìn)行對比,如圖16所示。

    分別對風(fēng)電功率實(shí)際值和預(yù)測值對應(yīng)的風(fēng)電切除量進(jìn)行計(jì)算,可得圖17。

    由圖17可知,圖中風(fēng)電切除量最小誤差僅為0.8%,最大誤差為9.2%,其主要原因在于風(fēng)電功率預(yù)測誤差為11.1%。次日風(fēng)電切除量平均絕對誤差僅為3.8%,整體誤差較小。由此可得超短期風(fēng)電功率預(yù)測滾動修正下的風(fēng)電切除方案精確度較高。在工程實(shí)際操作運(yùn)行時,可通過超短期風(fēng)電功率預(yù)測、實(shí)時風(fēng)電功率對次日風(fēng)電功率預(yù)測值進(jìn)行滾動修正求解,從而提高風(fēng)電功率切除的精確性。

    圖16 風(fēng)電功率預(yù)測值和實(shí)際值的對比

    圖17 風(fēng)電功率預(yù)測值和實(shí)際值對應(yīng)的風(fēng)電切除量

    4.4 風(fēng)電最大并網(wǎng)容量及切除量的確定

    上述部分驗(yàn)證了兩種功率缺額場景下風(fēng)電最大并網(wǎng)容量及其切除量計(jì)算方法的正確性和計(jì)算結(jié)果的準(zhǔn)確性,并針對風(fēng)電功率預(yù)測誤差下的風(fēng)電切除方案進(jìn)行相應(yīng)修正處理。在實(shí)際取值時還須根據(jù)式(24)取兩種場景下、各時段的風(fēng)電并網(wǎng)容量較小值和切除量較大值,據(jù)此得出如圖18所示結(jié)果。

    圖18 一天中風(fēng)電最大并網(wǎng)容量及切除量

    5 結(jié)論

    本文應(yīng)用頻率安全約束和臨界慣量,對風(fēng)電場最大并網(wǎng)容量及其切除量進(jìn)行計(jì)算,得出如下結(jié)論:

    1) 根據(jù)頻率跌落最低點(diǎn)和最大頻率變化率約束指標(biāo)分時段求得的臨界慣量,可快速準(zhǔn)確地判斷電網(wǎng)存在頻率安全隱患的時段。

    2) 電網(wǎng)慣量和臨界慣量均具有時變特征,對于電網(wǎng)實(shí)際慣量小于臨界慣量時段,應(yīng)用本文計(jì)算方法可對風(fēng)電場并網(wǎng)容量進(jìn)行定量切除,從而保證電網(wǎng)故障下的頻率安全。

    3) 隨著風(fēng)電虛擬慣性控制和儲能技術(shù)的進(jìn)一步發(fā)展,未來風(fēng)電場將會具備一定的慣性響應(yīng)能力。在動態(tài)頻率安全約束下的風(fēng)電并網(wǎng)能力將會發(fā)生改變,這將是下一步研究的內(nèi)容。

    附表A1 發(fā)電機(jī)組參數(shù)

    Attached Table A1 Generator set parameters

    發(fā)電機(jī)H/sKT/s G16.716.678 G25.214.286 G34.816.677 G45.216.676 G54.514.286 G65.416.678 G75.416.677 G84.216.677 G95.416.676 G104.214.286

    附圖A1 次日分時段風(fēng)電最大并網(wǎng)功率計(jì)算流程圖

    Attached Fig. A1 Calculation flow chart of the maximum grid-connected capacity of wind power in different periods of the next day

    附表A2 各時段發(fā)電計(jì)劃

    Attached Table A2 Power generation plan for each time period

    時間/hG1并網(wǎng)功率/MWG2并網(wǎng)功率/MWG3并網(wǎng)功率/MWG4并網(wǎng)功率/MWG5并網(wǎng)功率/MW 1800520650640510 2900320450440510 3900320550440510 4900320750440510 5900320750440510 6900320750440510 7900520750440510 8900520750740510 9900520750740510 101300720750740510 111300720750740510 121000720750440510 131000720750440510 141000520750440510 151000520750440510 161000520750440510 171000520950440510 181000520950440510 191300920950440510 201300920950440510 21900920950140510 22900920950140510 23900920950140510 24900920950140510 時間/hG6并網(wǎng)功率/MWG7并網(wǎng)功率/MWG8并網(wǎng)功率/MWG9并網(wǎng)功率/MWG10并網(wǎng)功率/MW 14505405401030250 24505405401030250 34505405401030350 44505405401030350 54505405401030350 64505405401030350 74505405401030350 84505405401030350 94505405401030750 104505405401030750 114505409401030750 124505409401030750 134505409401030750 144505409401030350 154505409401030350 164505409401030350 174505409401030350 184505409401030350 194505409401030350 204505409401030350 214505409401030350 224505409401030350 234505406401030350 244505406401030150

    附圖A2 負(fù)荷突增下的風(fēng)電最大并網(wǎng)功率和風(fēng)電切除量

    Attached Fig. A2 Maximum grid-connected capacity and removal of wind power under sudden load increase

    附A3 負(fù)荷突增下的風(fēng)電最大并網(wǎng)功率和風(fēng)電切除量

    Attached Table A3 Maximum grid-connected capacity and removal of wind power under sudden load increase

    時間/h風(fēng)電最大并網(wǎng)功率/MW風(fēng)電切除量/MW時段/h風(fēng)電最大并網(wǎng)功率/MW風(fēng)電切除量/MW 1940.67408.0013505.310.00 21149.4877.0014465.910.00 31025.640.0015491.870.00 4968.820.0016570.130.00 51034.000.0017657.090.00 6929.240.0018748.4161.00 7959.560.0019599.10489.00 8806.880.0020924.50451.00 9600.6923.0021738.74471.00 10435.43175.0022929.27167.00 11592.820.00231026.370.00 12527.280.0024799.610.00

    附圖A3 同步機(jī)組脫網(wǎng)下風(fēng)電最大并網(wǎng)功率和風(fēng)電切除量

    Attached Fig. A3 Maximum grid-connected capacity and cut-off amount of wind power when synchronized generators are off-grid

    附表A4 同步機(jī)組脫網(wǎng)下風(fēng)電最大并網(wǎng)功率和風(fēng)電切除量

    Attached Table A4 Maximum grid-connected capacity and cut-off amount of wind power when synchronized generators are off-grid

    時間/h風(fēng)電最大并網(wǎng)功率/MW風(fēng)電切除量/MW時段/h風(fēng)電最大并網(wǎng)功率/MW風(fēng)電切除量/MW 1985.67363.0013505.310.00 2891.48335.0014465.910.00 3780.64245.0015491.870.00 4811.82157.0016570.130.00 5870.00164.0017657.090.00 6929.240.0018809.410.00 7959.560.00191088.100.00 8806.880.00201375.500.00 9623.690.00211209.740.00 10610.430.00221096.270.00 11592.820.00231026.370.00 12527.280.0024799.610.00

    [1] 趙福林, 俞嘯玲, 杜詩嘉, 等. 計(jì)及需求響應(yīng)的含大規(guī)模風(fēng)電并網(wǎng)下電力系統(tǒng)靈活性評估[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(1): 42-51.

    ZHAO Fulin, YU Xiaoling, DU Shijia, et al. Assessment on flexibility of a power grid with large-scale wind farm integration considering demand response[J]. Power System Protection and Control, 2021, 49(1): 42-51.

    [2] 李德鑫, 宗崇林, 黃大為, 等. 大規(guī)模風(fēng)電匯聚外送的區(qū)域送端電力系統(tǒng)兩級優(yōu)化調(diào)度[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(11): 35-44.

    LI Dexin, ZONG Chonglin, HUANG Dawei, et al. Two-level optimal dispatching of regional transmission end power system for large-scale wind power convergence and external transmission[J]. Power System Protection and Control, 2021, 49(11): 35-44.

    [3] 孫偉卿, 羅靜, 張婕. 高比例風(fēng)電接入的電力系統(tǒng)儲能容量配置及影響因素分析[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(15): 9-18.

    SUN Weiqing, LUO Jing, ZHANG Jie. Energy storage capacity allocation and influence factor analysis of a power system with a high proportion of wind power[J]. Power System Protection and Control, 2021, 49(15): 9-18.

    [4] 孫華東, 許濤, 郭強(qiáng), 等. 英國“8?9”大停電事故分析及對中國電網(wǎng)的啟示[J]. 中國電機(jī)工程學(xué)報(bào), 2019, 39(21): 6183-6192.

    SUN Huadong, XU Tao, GUO Qiang, et al. Analysis on blackout in Great Britain power grid on August 9th, 2019 and its enlightenment to power grid in China[J]. Proceedings of the CSEE, 2019, 39(21): 6183-6192.

    [5] ERDIWANSYAH, MAHIDIN, HUSIN H, et al. A critical review of the integration of renewable energy sources with various technologies[J]. Protection and Control of Modern Power Systems, 2021, 6(1): 37-54.

    [6] 劉英琪, 謝敏, 韋薇, 等. 高比例風(fēng)電接入的電力系統(tǒng)靈活性評估與優(yōu)化[J]. 電力建設(shè), 2019, 40(9): 1-10.

    LIU Yingqi, XIE Min, WEI Wei, et al. Assessment and optimization for power system flexibility with high proportion of wind power[J]. Electric Power Construction, 2019, 40(9): 1-10.

    [7] 陳赟, 陳得治, 馬世英, 等. 風(fēng)光火打捆交直流外送系統(tǒng)的高頻切機(jī)方案研究[J]. 電網(wǎng)技術(shù), 2016, 40(1): 186-192.

    CHEN Yun, CHEN Dezhi, MA Shiying, et al. Studies on high-frequency generator tripping strategy for wind-photovoltaic-thermal-bundled power transmitted by AC/DC system[J]. Power System Technology, 2016, 40(1): 186-192.

    [8] 丁劍, 邱躍豐, 孫華東, 等. 大規(guī)模風(fēng)電接入下風(fēng)電機(jī)組切機(jī)措施研究[J]. 中國電機(jī)工程學(xué)報(bào), 2011, 31(19): 25-36.

    DING Jian, QIU Yuefeng, SUN Huadong, et al. Consideration of wind generator tripping under large-scale wind power integration[J]. Proceedings of the CSEE, 2011, 31(19): 25-36.

    [9] 岳涵, 邵廣惠, 孫銘澤, 等. 含大規(guī)模新能源集中外送型電網(wǎng)高頻第三道防線配置方法[J]. 電網(wǎng)技術(shù), 2019, 43(5): 1766-1771.

    YUE Han, SHAO Guanghui, SUN Mingze, et al. A configuration method of high frequency third defensive line for power-integrated-exporting grid with large-scale renewable energy[J]. Power System Technology, 2019, 43(5): 1766-1771.

    [10]王銳, 顧偉, 孫蓉, 等. 基于概率最優(yōu)潮流的風(fēng)電接入能力分析[J]. 電網(wǎng)技術(shù), 2011, 35(12): 214-220.

    WANG Rui, GU Wei, SUN Rong, et al. Analysis on wind power penetration limit based on probabilistically optimal power flow[J]. Power System Technology, 2011, 35(12): 214-220.

    [11]湯雪松, 殷明慧, 鄒云. 考慮風(fēng)速相關(guān)性的風(fēng)電穿透功率極限的改進(jìn)計(jì)算[J]. 電網(wǎng)技術(shù), 2015, 39(2): 420-425.

    TANG Xuesong, YIN Minghui, ZOU Yun. An improved method to calculate wind power penetration limit considering wind speed correlation[J]. Power System Technology, 2015, 39(2): 420-425.

    [12]曲正偉, 王京波, 王云靜, 等. 考慮運(yùn)行風(fēng)險約束的風(fēng)電場群準(zhǔn)入容量分析[J]. 電網(wǎng)技術(shù), 2014, 38(7): 1861-1866.

    QU Zhengwei, WANG Jingbo, WANG Yunjing, et al. Analysis on acceptable capacity of wind farm group considering operation risk constraints[J]. Power System Technology, 2014, 38(7): 1861-1866.

    [13] 李世春, 唐紅艷, 鄧長虹, 等. 考慮頻率約束及風(fēng)電機(jī)組調(diào)頻的風(fēng)電穿透功率極限計(jì)算[J]. 電力系統(tǒng)自動化, 2019, 43(4): 33-43.

    LI Shichun, TANG Hongyan, DENG Changhong, et al. Calculation of wind power penetration limit involving frequency constraints and frequency regulation of wind turbines[J]. Automation of Electric Power Systems, 2019, 43(4): 33-43.

    [14] AHMADYAR A, RIAZ S, VERBIC G, et al. A framework for assessing renewable integration limits with respect to frequency performance[J]. IEEE Transactions on Power Systems, 2018, 33(4): 4444-4453.

    [15] 謝巖, 馬世英, 唐曉駿. 頻率穩(wěn)定約束條件下風(fēng)電并網(wǎng)能力的確定方法[J]. 電力系統(tǒng)及其自動化學(xué)報(bào), 2021, 33(8): 130-137.

    XIE Yan, MA Shiying, TANG Xiaojun. Determination of wind power grid connection capability under frequency stability constraints[J]. Proceedings of the CSU-EPSA, 2021, 33(8): 130-137.

    [16] 郭小龍, 畢天姝, 劉方蕾, 等. 風(fēng)、光高滲透率電網(wǎng)中考慮頻率穩(wěn)定的可再生能源承載力研究[J]. 可再生能源, 2020, 38(1): 84-90.

    GUO Xiaolong, BI Tianshu, LIU Fanglei, et al. Estimating maximum penetration level of renewable energy based on frequency stability constrains in networks with high penetration wind and photovoltaic energy[J]. Renewable Energy, 2020, 38(1): 84-90.

    [17] 李世春, 夏智雄, 程緒長, 等. 基于類噪聲擾動的電網(wǎng)慣量常態(tài)化連續(xù)估計(jì)方法[J]. 中國電機(jī)工程學(xué)報(bào), 2020, 40(14): 4430-4439, 4723.

    LI Shichun, XIA Zhixiong, CHENG Xuchang, et al. Continuous estimation method of power system inertia based on ambient disturbance[J]. Proceedings of the CSEE, 2020, 40(14): 4430-4439, 4723.

    [18] LI Zhi, YE Lin, ZHAO Yongning, et al. Short-term wind power prediction based on extreme learning machine with error correction[J]. Protection and Control of Modern Power Systems, 2016, 1(1): 9-16.

    [19] HUANG H, JU P, JIN Y, et al. Generic system frequency response model for power grids with different generations[J]. IEEE Access, 2020, 8: 14314-14321.

    [20] JU Ping, ZHENG Yi, JIN Yuqing, et al. Analytic assessment of the power system frequency security[J]. IET Generation, Transmission & Distribution, 2021, 15(15): 2215-2225.

    [21] 吳雪蓮, 李兆偉, 劉福鎖, 等. 大功率擾動下計(jì)及系統(tǒng)頻率分布特性的緊急控制策略研究[J]. 電力系統(tǒng)保護(hù)與控制, 2021, 49(3): 104-114.

    WU Xuelian, LI Zhaowei, LIU Fusuo, et al. Analysis of the emergency control strategy of system frequency considering system frequency distribution characteristics under large power disturbance[J]. Power System Protection and Control, 2021, 49(3): 104-114.

    [22] EGIDO I, FERNANDEZ-BERNAL F, CENTENO P, et al. Maximum frequency deviation calculation in small isolated power systems[J]. IEEE Transactions on Power Systems, 2009, 24(4): 1731-1738.

    [23] 王霞, 應(yīng)黎明, 盧少平. 考慮動態(tài)頻率約束的一次調(diào)頻和二次調(diào)頻聯(lián)合優(yōu)化模型[J]. 電網(wǎng)技術(shù), 2020, 44(8): 2858-2867.

    WANG Xia, YING Liming, LU Shaoping. Joint optimization model for primary and secondary frequency regulation considering dynamic frequency constraint[J]. Power System Technology, 2020, 44(8): 2858-2867.

    [24]葉婧, 林濤, 張磊, 等. 考慮動態(tài)頻率約束的含高滲透率光伏電源的孤立電網(wǎng)機(jī)組組合[J]. 電工技術(shù)學(xué)報(bào), 2017, 32(13): 194-202.

    YE Jing, LIN Tao, ZHANG Lei, et al. Isolated grid unit commitment with dynamic frequency constraint considering photovoltaic power plants participating in frequency regulation[J]. Transactions of China Electrotechnical Society, 2017, 32(13): 194-202.

    [25] 樊艷芳, 鐘顯, 常喜強(qiáng), 等. 頻率的時空分布對低頻減載的影響研究[J]. 電力系統(tǒng)保護(hù)與控制, 2015, 43(1): 55-60.

    FAN Yanfang, ZHONG Xian, CHANG Xiqiang, et al. Research of the UFLS under the frequency characteristics of temporal and spatial distribution[J]. Power System Protection and Control, 2015, 43(1): 55-60.

    [26] 楊茂, 張羅賓. 基于數(shù)據(jù)驅(qū)動的超短期風(fēng)電功率預(yù)測綜述[J]. 電力系統(tǒng)保護(hù)與控制, 2019, 47(13): 171-186.

    YANG Mao, ZHANG Luobin. Review on ultra-short term wind power forecasting based on data-driven approach[J]. Power System Protection and Control, 2019, 47(13): 171-186.

    [27] 趙征, 汪向碩. 基于CEEMD和改進(jìn)時間序列模型的超短期風(fēng)功率多步預(yù)測[J]. 太陽能學(xué)報(bào), 2020, 41(7): 352-358.

    ZHAO Zheng, WANG Xiangshuo. Ultra-short-term multi-step wind power prediction based on CEEMD and improved time series model[J]. Acta Energiae Solaris Sinica, 2020, 41(7): 352-358.

    [28] LIU J, SHI Q, HAN R, et al. A hybrid GA–PSO–CNN model for ultra-short-term wind power forecasting[J]. Energies, 2021, 14(20): 1-22.

    [28] WU Q, GUAN F, Lü C, et al. Ultra-short-term multi-step wind power forecasting based on CNN-LSTM[J]. IET Renewable Power Generation, 2021, 15(5): 1019-1029.

    Method for limiting wind power output in time periods based on frequency safety constraints and a critical inertia calculation

    LI Shichun, TIAN Bingjie, LI Huizi, LUO Ying, HUANG Senyan, XU Songlin

    (College of Electrical Engineering and New Energy, China Three Gorges University, Yichang 443002, China)

    The high proportion of wind power access to power systems leads to a reduction of the equivalent inertia of the system. This increases the possibility of frequency instability caused by the power grid. In order to ensure the safe and stable operation of the frequency during the day, and at the same time avoid the unreasonable step of removing wind power in a large area, a method of limiting wind power output based on frequency safety constraints and a critical inertia calculation is proposed. By introducing a frequency nadir and maximum rate of change of frequency constraint index, and considering the changes in the operation mode within the day, the critical inertia of the grid and the corresponding maximum wind power grid-connected capacity are calculated in time periods. A time period less than the critical inertia is determined based on the actual inertia of the grid, to determine the amount of wind power removal based on the maximum wind power grid-connected capacity. This is to eliminate hidden frequency safety hazards at various times of the day. Finally, a calculation example proves that the method proposed in this paper can accurately screen the time period of frequency safety risk in the grid day, and the calculated amount of wind power cut in time can effectively eliminate the hidden frequency safety hazard.

    high proportion of wind power; frequency safety constraints; critical inertia; maximum wind power grid-connected capacity; wind power removal

    10.19783/j.cnki.pspc.211322

    2021-09-29;

    2021-12-23

    李世春(1984—),男,博士,副教授,研究方向?yàn)楹履茉措娏ο到y(tǒng)運(yùn)行與控制;E-mail: lschunu_023@126.com

    田冰杰(1997—),男,碩士研究生,研究方向?yàn)楹履茉措娏ο到y(tǒng)運(yùn)行與控制。E-mail: 1354993460@qq.com

    國家自然科學(xué)基金項(xiàng)目資助(51907104)

    This work is supported by the National Natural Science Foundation of China (No. 51907104).

    (編輯 周金梅)

    猜你喜歡
    慣量電功率時段
    基于PCC-CNN-GRU的短期風(fēng)電功率預(yù)測
    并網(wǎng)模式下虛擬同步發(fā)電機(jī)的虛擬慣量控制策略
    輕松上手電功率
    你會計(jì)算電功率嗎
    解讀電功率
    四個養(yǎng)生黃金時段,你抓住了嗎
    一種基于模擬慣量偏差的電慣量控制算法
    低階可約慣量任意符號模式矩陣的刻畫
    三階不可約零-非零模式中的幾乎慣量任意模式
    傍晚是交通事故高發(fā)時段
    国产精品久久久久久久电影| 一级毛片电影观看 | 色综合站精品国产| 亚洲一区二区三区色噜噜| 级片在线观看| 六月丁香七月| 婷婷六月久久综合丁香| 看十八女毛片水多多多| 特大巨黑吊av在线直播| 久久久久久国产a免费观看| 一区二区三区高清视频在线| 亚洲av中文字字幕乱码综合| 日韩 亚洲 欧美在线| 国产一区二区三区在线臀色熟女| 22中文网久久字幕| 啦啦啦观看免费观看视频高清| 精品一区二区三区视频在线| 亚洲自偷自拍三级| 日本黄大片高清| 欧美潮喷喷水| 热99re8久久精品国产| 欧美性感艳星| 91aial.com中文字幕在线观看| 不卡视频在线观看欧美| 中文字幕免费在线视频6| 能在线免费看毛片的网站| 国产精品免费一区二区三区在线| 国产午夜福利久久久久久| 国内精品一区二区在线观看| 天堂av国产一区二区熟女人妻| 一区二区三区免费毛片| 99久久精品热视频| 人体艺术视频欧美日本| 亚洲欧美日韩高清在线视频| 国产精品一区www在线观看| 欧美性猛交黑人性爽| 亚洲自拍偷在线| 亚洲国产精品合色在线| 毛片女人毛片| 一级二级三级毛片免费看| 熟女电影av网| 国产探花极品一区二区| 欧美激情国产日韩精品一区| 黄色一级大片看看| av天堂中文字幕网| 日本黄色视频三级网站网址| 欧美三级亚洲精品| 插逼视频在线观看| 成人毛片60女人毛片免费| 国产成人freesex在线| 一个人看的www免费观看视频| 人体艺术视频欧美日本| 九九在线视频观看精品| 波野结衣二区三区在线| 日本成人三级电影网站| 国产午夜福利久久久久久| 不卡一级毛片| 日韩av不卡免费在线播放| 日日啪夜夜撸| www日本黄色视频网| 亚洲欧美成人综合另类久久久 | 人人妻人人澡人人爽人人夜夜 | 欧美一区二区精品小视频在线| 精品一区二区三区人妻视频| .国产精品久久| 亚洲欧美日韩东京热| 看黄色毛片网站| 日韩视频在线欧美| 亚洲欧美日韩高清在线视频| 欧美最黄视频在线播放免费| 国产综合懂色| 精品久久久久久久久久免费视频| 久久久a久久爽久久v久久| 国产精品野战在线观看| 美女被艹到高潮喷水动态| 亚洲在线观看片| 国产精品电影一区二区三区| 99久久精品一区二区三区| 精品国内亚洲2022精品成人| 国产又黄又爽又无遮挡在线| 亚洲国产精品久久男人天堂| 国产久久久一区二区三区| 又爽又黄a免费视频| 又爽又黄无遮挡网站| 亚洲一区高清亚洲精品| 成年av动漫网址| 亚洲第一电影网av| 久久精品国产亚洲av天美| 嘟嘟电影网在线观看| 免费看日本二区| 国产乱人偷精品视频| 欧美潮喷喷水| 白带黄色成豆腐渣| 午夜福利高清视频| 大型黄色视频在线免费观看| 亚洲精品乱码久久久久久按摩| 麻豆久久精品国产亚洲av| 久久热精品热| 在线观看av片永久免费下载| 蜜臀久久99精品久久宅男| 亚洲精品日韩av片在线观看| 97人妻精品一区二区三区麻豆| 国产三级在线视频| 悠悠久久av| 此物有八面人人有两片| 国产av在哪里看| 久久久午夜欧美精品| 亚洲精品日韩av片在线观看| 久久欧美精品欧美久久欧美| 久久99精品国语久久久| 在线观看一区二区三区| 丰满的人妻完整版| 国产精品三级大全| a级毛片a级免费在线| 黄色视频,在线免费观看| 两个人视频免费观看高清| 亚洲五月天丁香| 久久久久网色| 精品无人区乱码1区二区| 麻豆av噜噜一区二区三区| 99久国产av精品| 男女那种视频在线观看| 亚洲人成网站高清观看| 看黄色毛片网站| 亚洲在久久综合| 国产精品乱码一区二三区的特点| 婷婷精品国产亚洲av| 身体一侧抽搐| 国产av在哪里看| 日日啪夜夜撸| 日日啪夜夜撸| 黄片无遮挡物在线观看| 久久婷婷人人爽人人干人人爱| 精品午夜福利在线看| 国产精品爽爽va在线观看网站| 免费人成视频x8x8入口观看| 乱码一卡2卡4卡精品| 国产乱人偷精品视频| 日本一本二区三区精品| av在线亚洲专区| 久久久久久久久中文| 伊人久久精品亚洲午夜| av福利片在线观看| 男的添女的下面高潮视频| 国产一区二区三区av在线 | 亚洲欧美日韩高清专用| 非洲黑人性xxxx精品又粗又长| 桃色一区二区三区在线观看| 亚洲av电影不卡..在线观看| 国产亚洲av片在线观看秒播厂 | 波野结衣二区三区在线| 波多野结衣巨乳人妻| 国产伦理片在线播放av一区 | 婷婷色综合大香蕉| 黑人高潮一二区| 日韩强制内射视频| 国产黄色小视频在线观看| 免费人成视频x8x8入口观看| 亚洲欧美清纯卡通| 99久久精品国产国产毛片| 能在线免费观看的黄片| 一边摸一边抽搐一进一小说| 亚洲成人久久爱视频| 老司机福利观看| 日韩视频在线欧美| 黄色配什么色好看| 日本-黄色视频高清免费观看| 夜夜看夜夜爽夜夜摸| 男人舔奶头视频| 国产精品无大码| 日韩欧美精品v在线| 精品不卡国产一区二区三区| 国产精品.久久久| av福利片在线观看| 伊人久久精品亚洲午夜| 亚洲精华国产精华液的使用体验 | 久久亚洲精品不卡| АⅤ资源中文在线天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 国产精品电影一区二区三区| 99九九线精品视频在线观看视频| 变态另类成人亚洲欧美熟女| 国产综合懂色| av视频在线观看入口| 欧美日韩乱码在线| 最近2019中文字幕mv第一页| 如何舔出高潮| 久久久a久久爽久久v久久| 国产免费男女视频| 亚洲一区二区三区色噜噜| 成人三级黄色视频| 色噜噜av男人的天堂激情| 亚洲国产欧美在线一区| 97超视频在线观看视频| 最近视频中文字幕2019在线8| 午夜免费激情av| 成年免费大片在线观看| 久久久a久久爽久久v久久| 亚洲,欧美,日韩| 不卡一级毛片| av免费观看日本| 欧美一区二区国产精品久久精品| 嫩草影院新地址| 在线免费十八禁| 亚洲图色成人| av在线亚洲专区| 国产精品一区二区性色av| 久久国产乱子免费精品| 欧美一区二区国产精品久久精品| 一区二区三区免费毛片| 国产精品一二三区在线看| 国产在线精品亚洲第一网站| 麻豆一二三区av精品| 国产老妇女一区| 欧美又色又爽又黄视频| 国产伦精品一区二区三区视频9| 国产精品日韩av在线免费观看| 免费人成视频x8x8入口观看| 久久精品国产亚洲网站| 亚洲国产精品国产精品| 亚洲国产日韩欧美精品在线观看| 亚洲中文字幕日韩| www日本黄色视频网| 乱人视频在线观看| 久久国内精品自在自线图片| 成年av动漫网址| 欧美日韩国产亚洲二区| 亚洲成人久久性| 久久这里只有精品中国| 免费大片18禁| 成人性生交大片免费视频hd| 婷婷色av中文字幕| 联通29元200g的流量卡| 在线a可以看的网站| 男人的好看免费观看在线视频| 日本欧美国产在线视频| 国产又黄又爽又无遮挡在线| 黄色日韩在线| 国产蜜桃级精品一区二区三区| 在线免费观看的www视频| 国产成人午夜福利电影在线观看| 人妻夜夜爽99麻豆av| 久久人人爽人人爽人人片va| 一本久久精品| 一级毛片aaaaaa免费看小| 国产精品久久久久久亚洲av鲁大| 免费大片18禁| 国产精品蜜桃在线观看 | 久久草成人影院| 一级毛片电影观看 | 看黄色毛片网站| 日本三级黄在线观看| 国产成人精品一,二区 | 亚洲av男天堂| 亚洲精品粉嫩美女一区| 欧美成人a在线观看| 毛片女人毛片| 99热6这里只有精品| 如何舔出高潮| 中文在线观看免费www的网站| 赤兔流量卡办理| 内射极品少妇av片p| 白带黄色成豆腐渣| 日本-黄色视频高清免费观看| 日韩一本色道免费dvd| 中文欧美无线码| 99久久中文字幕三级久久日本| 国产国拍精品亚洲av在线观看| 九九热线精品视视频播放| 91久久精品国产一区二区成人| 免费观看人在逋| 久久久久久国产a免费观看| 丰满乱子伦码专区| 日韩三级伦理在线观看| 3wmmmm亚洲av在线观看| 日韩一本色道免费dvd| 国产女主播在线喷水免费视频网站 | 一本精品99久久精品77| 国产欧美日韩精品一区二区| 欧美高清成人免费视频www| 国产精品1区2区在线观看.| 精品熟女少妇av免费看| av天堂在线播放| 久久久成人免费电影| 国产精品一区二区三区四区久久| 国产一区二区三区在线臀色熟女| 两个人视频免费观看高清| 老司机影院成人| 日本撒尿小便嘘嘘汇集6| 特大巨黑吊av在线直播| 精品久久久久久久久久免费视频| 亚洲在久久综合| 夜夜看夜夜爽夜夜摸| 村上凉子中文字幕在线| 丰满的人妻完整版| 久久久精品94久久精品| av卡一久久| 久99久视频精品免费| 亚洲久久久久久中文字幕| 啦啦啦观看免费观看视频高清| av在线亚洲专区| 日本-黄色视频高清免费观看| 人人妻人人澡欧美一区二区| 菩萨蛮人人尽说江南好唐韦庄 | 久久久欧美国产精品| 九草在线视频观看| 男的添女的下面高潮视频| 少妇被粗大猛烈的视频| 丰满人妻一区二区三区视频av| 亚洲无线在线观看| 99国产精品一区二区蜜桃av| 亚洲av免费高清在线观看| 婷婷色av中文字幕| 美女被艹到高潮喷水动态| 日韩,欧美,国产一区二区三区 | 国产探花在线观看一区二区| 国产成人一区二区在线| 成人毛片60女人毛片免费| 哪个播放器可以免费观看大片| 日韩欧美三级三区| 久久综合国产亚洲精品| 麻豆国产av国片精品| 亚洲在线观看片| 成人永久免费在线观看视频| 中文亚洲av片在线观看爽| 成人国产麻豆网| 亚洲最大成人av| 国产成人a∨麻豆精品| 国产91av在线免费观看| 国内精品美女久久久久久| 中国国产av一级| 欧美高清性xxxxhd video| 国产熟女欧美一区二区| 日韩强制内射视频| 99在线视频只有这里精品首页| 久久久久久久久久久丰满| 亚洲,欧美,日韩| 成人午夜精彩视频在线观看| 小蜜桃在线观看免费完整版高清| 高清午夜精品一区二区三区 | 国产中年淑女户外野战色| 国产乱人偷精品视频| 久久精品国产亚洲av涩爱 | 一边摸一边抽搐一进一小说| 午夜精品国产一区二区电影 | 中国美女看黄片| 老司机福利观看| 亚洲,欧美,日韩| 日韩一本色道免费dvd| 又粗又硬又长又爽又黄的视频 | 久久精品夜色国产| 乱码一卡2卡4卡精品| 黄色配什么色好看| 日韩欧美精品v在线| 嫩草影院精品99| 国产视频内射| 少妇的逼水好多| 午夜福利在线观看免费完整高清在 | 国产成人精品婷婷| 一级毛片我不卡| 日本在线视频免费播放| 国产精品乱码一区二三区的特点| 亚洲欧美精品专区久久| 亚洲国产欧洲综合997久久,| 天天躁日日操中文字幕| 亚洲成a人片在线一区二区| 成人一区二区视频在线观看| 亚洲av二区三区四区| 国内揄拍国产精品人妻在线| 亚洲av一区综合| 成人性生交大片免费视频hd| 午夜久久久久精精品| 九草在线视频观看| 91aial.com中文字幕在线观看| 亚洲精品粉嫩美女一区| 噜噜噜噜噜久久久久久91| 卡戴珊不雅视频在线播放| 久久久久性生活片| 老女人水多毛片| 国产午夜精品论理片| 伊人久久精品亚洲午夜| 男的添女的下面高潮视频| 午夜福利在线在线| 哪里可以看免费的av片| 一本久久精品| 天堂影院成人在线观看| 成年版毛片免费区| av免费观看日本| 在线观看美女被高潮喷水网站| 又爽又黄无遮挡网站| 97超碰精品成人国产| 日韩欧美 国产精品| 97超视频在线观看视频| 一边摸一边抽搐一进一小说| 欧美高清性xxxxhd video| 免费看光身美女| 26uuu在线亚洲综合色| 偷拍熟女少妇极品色| 在线免费十八禁| 午夜精品一区二区三区免费看| 精品不卡国产一区二区三区| 国产蜜桃级精品一区二区三区| 性插视频无遮挡在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 九草在线视频观看| 26uuu在线亚洲综合色| av视频在线观看入口| 如何舔出高潮| www.av在线官网国产| 草草在线视频免费看| 18禁在线无遮挡免费观看视频| 亚洲不卡免费看| 久久久久久久久久久免费av| 最新中文字幕久久久久| av在线播放精品| 性欧美人与动物交配| 国产精品乱码一区二三区的特点| 看非洲黑人一级黄片| 18禁在线无遮挡免费观看视频| 国产综合懂色| 国产高潮美女av| 国产精品三级大全| 亚洲成人久久爱视频| 特级一级黄色大片| 男女下面进入的视频免费午夜| 成人无遮挡网站| 日本一本二区三区精品| 国产精品综合久久久久久久免费| 男女边吃奶边做爰视频| 在线免费十八禁| 在线播放无遮挡| 欧美日韩精品成人综合77777| 国产高清有码在线观看视频| 又爽又黄a免费视频| 悠悠久久av| 深爱激情五月婷婷| 成年av动漫网址| 熟女电影av网| 亚洲熟妇中文字幕五十中出| 中文字幕精品亚洲无线码一区| 国产成人freesex在线| 免费电影在线观看免费观看| 精品欧美国产一区二区三| 婷婷亚洲欧美| 少妇的逼好多水| 99热全是精品| 白带黄色成豆腐渣| 在线免费观看的www视频| 秋霞在线观看毛片| 少妇高潮的动态图| 国产成人91sexporn| 国产亚洲av嫩草精品影院| 久久这里有精品视频免费| 最好的美女福利视频网| 久久6这里有精品| 亚洲av成人精品一区久久| 最好的美女福利视频网| 69av精品久久久久久| 最近最新中文字幕大全电影3| 日韩av不卡免费在线播放| 最近中文字幕高清免费大全6| 成人av在线播放网站| 久久久久免费精品人妻一区二区| 精品国产三级普通话版| 午夜精品在线福利| 亚洲欧美成人精品一区二区| 最新中文字幕久久久久| 人妻少妇偷人精品九色| 草草在线视频免费看| 欧美日韩乱码在线| 久久久久久久久久久丰满| 全区人妻精品视频| 91久久精品国产一区二区成人| 天堂网av新在线| 3wmmmm亚洲av在线观看| 国产亚洲av嫩草精品影院| 特级一级黄色大片| 黑人高潮一二区| 日本黄色片子视频| 亚洲中文字幕一区二区三区有码在线看| 免费大片18禁| av在线老鸭窝| 国产三级在线视频| 最近手机中文字幕大全| 国产伦精品一区二区三区四那| 婷婷色综合大香蕉| 在现免费观看毛片| 国产精品国产高清国产av| 国产男人的电影天堂91| 女人被狂操c到高潮| 国产精品久久久久久久电影| 综合色av麻豆| 日本三级黄在线观看| 日韩一本色道免费dvd| 少妇猛男粗大的猛烈进出视频 | 欧美高清性xxxxhd video| 三级国产精品欧美在线观看| 嫩草影院精品99| 直男gayav资源| 99久久久亚洲精品蜜臀av| 精品免费久久久久久久清纯| 成人一区二区视频在线观看| 男女下面进入的视频免费午夜| 亚洲国产精品成人久久小说 | 日韩强制内射视频| 中文在线观看免费www的网站| 精品人妻视频免费看| 美女cb高潮喷水在线观看| 国产成人一区二区在线| 亚洲熟妇中文字幕五十中出| 特大巨黑吊av在线直播| 又粗又爽又猛毛片免费看| 男女视频在线观看网站免费| 国产精品人妻久久久久久| 日日撸夜夜添| 国产精品久久久久久久久免| 不卡一级毛片| 狠狠狠狠99中文字幕| 免费观看精品视频网站| 久久精品夜夜夜夜夜久久蜜豆| 中文字幕精品亚洲无线码一区| 国产高清不卡午夜福利| 嫩草影院精品99| 一本精品99久久精品77| 亚洲自拍偷在线| 乱系列少妇在线播放| 久久久久久久久久久免费av| 激情 狠狠 欧美| 欧美一级a爱片免费观看看| 欧美日韩精品成人综合77777| 能在线免费观看的黄片| av天堂中文字幕网| 日本免费a在线| 午夜精品一区二区三区免费看| 91av网一区二区| 国产精品美女特级片免费视频播放器| 波野结衣二区三区在线| 我要搜黄色片| 18+在线观看网站| 免费不卡的大黄色大毛片视频在线观看 | 欧美zozozo另类| 国产伦精品一区二区三区四那| 精品一区二区免费观看| 美女高潮的动态| 两个人的视频大全免费| 春色校园在线视频观看| 搞女人的毛片| 91久久精品国产一区二区三区| 国产在线精品亚洲第一网站| 久久精品国产自在天天线| 日本三级黄在线观看| 亚洲欧美精品专区久久| 久久精品国产亚洲网站| 日日啪夜夜撸| 婷婷色av中文字幕| 亚洲欧美日韩无卡精品| 亚洲精品456在线播放app| 精华霜和精华液先用哪个| 嘟嘟电影网在线观看| 亚洲精品久久久久久婷婷小说 | 99九九线精品视频在线观看视频| 可以在线观看毛片的网站| 国产国拍精品亚洲av在线观看| 亚洲在线观看片| 国产精品国产三级国产av玫瑰| 别揉我奶头 嗯啊视频| avwww免费| 极品教师在线视频| 国产真实伦视频高清在线观看| 在线免费观看不下载黄p国产| 日本免费一区二区三区高清不卡| 国产淫片久久久久久久久| av在线蜜桃| 国产精品久久久久久精品电影小说 | 国产精品av视频在线免费观看| 久久草成人影院| 99国产极品粉嫩在线观看| 亚洲国产精品成人久久小说 | 69人妻影院| 国产精品不卡视频一区二区| 国产一区二区在线av高清观看| 日韩欧美国产在线观看| 国产高清激情床上av| 亚洲欧美精品综合久久99| 国产av在哪里看| 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 午夜精品在线福利| 夜夜看夜夜爽夜夜摸| 尤物成人国产欧美一区二区三区| 欧美日韩在线观看h| 少妇熟女欧美另类| 亚洲精品自拍成人| 麻豆精品久久久久久蜜桃| 亚洲成人中文字幕在线播放| 一卡2卡三卡四卡精品乱码亚洲| 久久精品夜夜夜夜夜久久蜜豆| 亚洲图色成人| 国产一区亚洲一区在线观看| 黄色欧美视频在线观看| 亚洲精品色激情综合| 午夜久久久久精精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产私拍福利视频在线观看| 尤物成人国产欧美一区二区三区| 欧美bdsm另类| 中文字幕免费在线视频6| avwww免费| 精品99又大又爽又粗少妇毛片| 18禁黄网站禁片免费观看直播| 中出人妻视频一区二区| 少妇人妻一区二区三区视频| 国产精品蜜桃在线观看 |