• <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免费视频播放| 亚洲中文av在线| 人妻人人澡人人爽人人| 久久久国产一区二区| 啦啦啦中文免费视频观看日本| 亚洲午夜精品一区,二区,三区| 亚洲精品一区蜜桃| 天天添夜夜摸| 亚洲美女黄色视频免费看| 多毛熟女@视频| av在线播放精品| 男女国产视频网站| 操出白浆在线播放| 国产成人精品无人区| 日本91视频免费播放| 日本欧美视频一区| 亚洲精品久久成人aⅴ小说| 免费少妇av软件| 日韩欧美一区视频在线观看| 成年动漫av网址| 亚洲精品一区蜜桃| 国产成人免费观看mmmm| 国产欧美日韩一区二区精品| 欧美日韩亚洲综合一区二区三区_| 日韩欧美免费精品| 在线观看免费日韩欧美大片| 香蕉国产在线看| 丝袜人妻中文字幕| bbb黄色大片| kizo精华| 新久久久久国产一级毛片| 天天影视国产精品| 国产欧美亚洲国产| 国内毛片毛片毛片毛片毛片| 叶爱在线成人免费视频播放| 国内毛片毛片毛片毛片毛片| 日韩视频在线欧美| 999久久久国产精品视频| 亚洲欧美日韩另类电影网站| 欧美日韩精品网址| 正在播放国产对白刺激| 99香蕉大伊视频| 久久狼人影院| 性高湖久久久久久久久免费观看| av在线老鸭窝| 俄罗斯特黄特色一大片| 老司机福利观看| 亚洲自偷自拍图片 自拍| 精品亚洲乱码少妇综合久久| 高清黄色对白视频在线免费看| 中文字幕色久视频| 久热爱精品视频在线9| 精品视频人人做人人爽| 超色免费av| 中国国产av一级| 亚洲精品成人av观看孕妇| 欧美变态另类bdsm刘玥| 国产精品成人在线| 国产成人免费无遮挡视频| 女人被躁到高潮嗷嗷叫费观| 老司机影院毛片| 久久这里只有精品19| 美女福利国产在线| 亚洲国产毛片av蜜桃av| 老鸭窝网址在线观看| 美女扒开内裤让男人捅视频| 日韩一区二区三区影片| 99精品欧美一区二区三区四区| 免费人妻精品一区二区三区视频| 天堂8中文在线网| 亚洲中文字幕日韩| 丁香六月天网| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 久久天堂一区二区三区四区| 国产成人系列免费观看| 国产精品香港三级国产av潘金莲| 韩国高清视频一区二区三区| 夜夜骑夜夜射夜夜干| 亚洲国产欧美网| xxxhd国产人妻xxx| 国精品久久久久久国模美| 午夜91福利影院| 欧美激情久久久久久爽电影 | 精品免费久久久久久久清纯 | 精品熟女少妇八av免费久了| 国产成人欧美在线观看 | 亚洲美女黄色视频免费看| 精品一区二区三区av网在线观看 | 国产日韩欧美在线精品| 中国美女看黄片| 久久影院123| 一级黄色大片毛片| 女性生殖器流出的白浆| 成年女人毛片免费观看观看9 | 国产97色在线日韩免费| 亚洲国产精品成人久久小说| 欧美变态另类bdsm刘玥| 一进一出抽搐动态| 99久久99久久久精品蜜桃| 乱人伦中国视频| 叶爱在线成人免费视频播放| 婷婷丁香在线五月| 一二三四在线观看免费中文在| 热99久久久久精品小说推荐| 99精品久久久久人妻精品| 国产成人精品久久二区二区免费| 国产又爽黄色视频| 1024视频免费在线观看| 亚洲国产精品一区三区| 乱人伦中国视频| 国产亚洲精品久久久久5区| 91麻豆av在线| 妹子高潮喷水视频| 国产在线免费精品| 亚洲综合色网址| 国产精品99久久99久久久不卡| 91麻豆精品激情在线观看国产 | 99久久人妻综合| 免费观看av网站的网址| 成人影院久久| 丝瓜视频免费看黄片| 久久中文字幕一级| 国产亚洲午夜精品一区二区久久| 欧美 日韩 精品 国产| 成年人免费黄色播放视频| 亚洲伊人久久精品综合| 免费在线观看视频国产中文字幕亚洲 | 国产三级黄色录像| 亚洲综合色网址| 日韩一区二区三区影片| 俄罗斯特黄特色一大片| 制服人妻中文乱码| 丝瓜视频免费看黄片| 天堂中文最新版在线下载| 一级片'在线观看视频| 国产精品免费视频内射| 视频区欧美日本亚洲| 亚洲成国产人片在线观看| 9热在线视频观看99| 精品乱码久久久久久99久播| 69精品国产乱码久久久| 曰老女人黄片| 成年人免费黄色播放视频| 两人在一起打扑克的视频| 老司机影院成人| 国产精品国产av在线观看| 别揉我奶头~嗯~啊~动态视频 | 国产成人精品无人区| 国产日韩一区二区三区精品不卡| 国产野战对白在线观看| 亚洲一码二码三码区别大吗| 久久久久久久久免费视频了| 亚洲成人免费av在线播放| 交换朋友夫妻互换小说| 一级毛片精品| 国产主播在线观看一区二区| av福利片在线| 日韩 亚洲 欧美在线| www.999成人在线观看| 三级毛片av免费| 免费在线观看视频国产中文字幕亚洲 | 两人在一起打扑克的视频| 亚洲一区二区三区欧美精品| 青春草亚洲视频在线观看| 欧美+亚洲+日韩+国产| 国产伦人伦偷精品视频| 黄色怎么调成土黄色| 极品少妇高潮喷水抽搐| 国产黄频视频在线观看| 久久精品国产综合久久久| 国产精品香港三级国产av潘金莲| 国产成人一区二区三区免费视频网站| 久久国产精品影院| 中文字幕人妻丝袜制服| 久久ye,这里只有精品| 嫩草影视91久久| 国产免费福利视频在线观看| 中文字幕av电影在线播放| 欧美av亚洲av综合av国产av| 免费高清在线观看视频在线观看| 亚洲欧美一区二区三区久久| 久久av网站| 免费不卡黄色视频| 超碰97精品在线观看| 最近中文字幕2019免费版| 免费av中文字幕在线| 老熟妇仑乱视频hdxx| 亚洲欧美一区二区三区黑人| 久久国产亚洲av麻豆专区| 日韩电影二区| 操美女的视频在线观看| 不卡一级毛片| 午夜福利在线观看吧| 啦啦啦在线免费观看视频4| 国产精品麻豆人妻色哟哟久久| 久久人妻熟女aⅴ| 亚洲avbb在线观看| 午夜福利影视在线免费观看| 真人做人爱边吃奶动态| 国产亚洲精品一区二区www | 一本—道久久a久久精品蜜桃钙片| 各种免费的搞黄视频| 啦啦啦视频在线资源免费观看| 九色亚洲精品在线播放| 国产日韩欧美亚洲二区| 精品一品国产午夜福利视频| 美女高潮到喷水免费观看| 如日韩欧美国产精品一区二区三区| 亚洲欧美成人综合另类久久久| av国产精品久久久久影院| 亚洲欧美精品综合一区二区三区| 建设人人有责人人尽责人人享有的| 无遮挡黄片免费观看| 久久精品人人爽人人爽视色| 久久久国产成人免费| 啦啦啦中文免费视频观看日本| 亚洲国产毛片av蜜桃av| 青草久久国产| 久久久久网色| 午夜福利在线免费观看网站| 亚洲精品国产色婷婷电影| 十八禁高潮呻吟视频| 最新的欧美精品一区二区| 91九色精品人成在线观看| 热99久久久久精品小说推荐| 久久久久国产一级毛片高清牌| 亚洲熟女毛片儿| 成年人免费黄色播放视频| 大码成人一级视频| 999久久久精品免费观看国产| 男女之事视频高清在线观看| 国产成人精品久久二区二区免费| 精品人妻熟女毛片av久久网站| 精品欧美一区二区三区在线| 亚洲精品中文字幕在线视频| 又大又爽又粗| 国产老妇伦熟女老妇高清| 日本黄色日本黄色录像| 中文字幕最新亚洲高清| 丝袜喷水一区| 亚洲天堂av无毛| 久久久久网色| 黄频高清免费视频| 精品少妇内射三级| 十八禁网站免费在线| 日韩中文字幕视频在线看片| av线在线观看网站| 久久精品亚洲av国产电影网| 国产精品一区二区精品视频观看| 亚洲性夜色夜夜综合| 91老司机精品| 久久久欧美国产精品| 亚洲免费av在线视频| 婷婷成人精品国产| 免费高清在线观看视频在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 色视频在线一区二区三区| 国产淫语在线视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲三区欧美一区| 母亲3免费完整高清在线观看| 大片免费播放器 马上看| 亚洲精品一二三| 可以免费在线观看a视频的电影网站| 久久久久久久大尺度免费视频| 亚洲一区二区三区欧美精品| 一边摸一边抽搐一进一出视频| 欧美日韩中文字幕国产精品一区二区三区 | 一本色道久久久久久精品综合| 亚洲国产欧美日韩在线播放| 欧美精品一区二区免费开放| 国产亚洲av高清不卡| 欧美xxⅹ黑人| 国产精品久久久久成人av| 女性生殖器流出的白浆| 黑人欧美特级aaaaaa片| 欧美97在线视频| 精品久久久久久电影网| 成人国产一区最新在线观看| 久9热在线精品视频| 亚洲欧美精品综合一区二区三区| 欧美精品一区二区大全| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久久久久免费视频了| 国产91精品成人一区二区三区 | 欧美精品亚洲一区二区| 国产成人精品无人区| 黑人操中国人逼视频| 99国产精品免费福利视频| 一本色道久久久久久精品综合| 亚洲国产欧美网| 国产免费av片在线观看野外av| 欧美日韩福利视频一区二区| 999久久久精品免费观看国产| 性高湖久久久久久久久免费观看| 免费不卡黄色视频| 欧美激情久久久久久爽电影 | 视频在线观看一区二区三区| 美女福利国产在线| av有码第一页| 看免费av毛片| 久久久久精品人妻al黑| 国产精品免费大片| 悠悠久久av| 麻豆av在线久日| 少妇猛男粗大的猛烈进出视频| 日本91视频免费播放| 精品人妻1区二区| 国产熟女午夜一区二区三区| 黑人巨大精品欧美一区二区mp4| 亚洲av日韩在线播放| 美国免费a级毛片| 五月开心婷婷网| 高潮久久久久久久久久久不卡| 国产av一区二区精品久久| 亚洲国产欧美在线一区| 国产成人系列免费观看| 欧美一级毛片孕妇| 欧美av亚洲av综合av国产av| 大码成人一级视频| 国产xxxxx性猛交| 日韩大片免费观看网站| 久久久国产欧美日韩av| 国产成人欧美在线观看 | 精品国产超薄肉色丝袜足j| 欧美精品人与动牲交sv欧美| 亚洲精品一卡2卡三卡4卡5卡 | 99九九在线精品视频| 美女高潮到喷水免费观看| 国产成人av激情在线播放| 另类亚洲欧美激情| 精品一区在线观看国产| h视频一区二区三区| 国产男人的电影天堂91| 日韩精品免费视频一区二区三区| 欧美中文综合在线视频| 亚洲精品久久成人aⅴ小说| 18禁裸乳无遮挡动漫免费视频| 亚洲伊人色综图| 免费观看a级毛片全部| 久久精品亚洲熟妇少妇任你| 成年人免费黄色播放视频| 各种免费的搞黄视频| a级毛片黄视频| 一个人免费看片子| 亚洲性夜色夜夜综合| 成人国产av品久久久| 91精品国产国语对白视频| 丝袜喷水一区| 两性夫妻黄色片| 又黄又粗又硬又大视频| 蜜桃在线观看..| 色婷婷av一区二区三区视频| 亚洲精品成人av观看孕妇| 成在线人永久免费视频| 伦理电影免费视频| 午夜免费成人在线视频| 这个男人来自地球电影免费观看| 亚洲成人手机| 在线永久观看黄色视频| 一区二区av电影网| 亚洲熟女毛片儿| 精品国产一区二区三区久久久樱花| 人人妻人人添人人爽欧美一区卜| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩高清在线视频 | 美女大奶头黄色视频| 男人爽女人下面视频在线观看| 性色av乱码一区二区三区2| 99久久99久久久精品蜜桃| 一个人免费在线观看的高清视频 | kizo精华| 欧美久久黑人一区二区| 99re6热这里在线精品视频| 国产高清国产精品国产三级| 亚洲国产精品999| 免费久久久久久久精品成人欧美视频| 精品国产一区二区三区四区第35| 国产精品欧美亚洲77777| 考比视频在线观看| 久久热在线av| 久久久水蜜桃国产精品网| 亚洲精品久久午夜乱码| 女人久久www免费人成看片| av国产精品久久久久影院| 黄色视频不卡| 中文字幕高清在线视频| 中文字幕另类日韩欧美亚洲嫩草| 男人爽女人下面视频在线观看| 交换朋友夫妻互换小说| 一边摸一边做爽爽视频免费| 欧美精品一区二区免费开放| 狠狠精品人妻久久久久久综合| 久9热在线精品视频| 中文字幕制服av| 免费在线观看黄色视频的| 丝瓜视频免费看黄片| 99久久精品国产亚洲精品| 精品人妻一区二区三区麻豆| 久久久久视频综合| 国产主播在线观看一区二区| 国产人伦9x9x在线观看| 日本黄色日本黄色录像| 99国产精品99久久久久| 久久人妻福利社区极品人妻图片| 欧美激情久久久久久爽电影 | 成人18禁高潮啪啪吃奶动态图| 少妇被粗大的猛进出69影院| 9色porny在线观看| 国产欧美亚洲国产| 女人爽到高潮嗷嗷叫在线视频| 中文字幕另类日韩欧美亚洲嫩草| 男人爽女人下面视频在线观看| 最近最新免费中文字幕在线| 欧美日韩黄片免| 亚洲一卡2卡3卡4卡5卡精品中文| av天堂在线播放| 国产精品亚洲av一区麻豆| 精品人妻一区二区三区麻豆| 久久九九热精品免费| 黄色视频在线播放观看不卡| 免费不卡黄色视频| 乱人伦中国视频| 成人国产av品久久久| 午夜福利在线免费观看网站| 热99国产精品久久久久久7| 久久久久久久久久久久大奶| 欧美日韩黄片免| 18禁黄网站禁片午夜丰满| 亚洲全国av大片| 三级毛片av免费| 一个人免费看片子| 国产片内射在线| 午夜福利,免费看| 久久 成人 亚洲| 欧美日韩亚洲综合一区二区三区_| 在线十欧美十亚洲十日本专区| 交换朋友夫妻互换小说| 精品国产超薄肉色丝袜足j| 99热国产这里只有精品6| 日本五十路高清| 亚洲综合色网址| 国产不卡av网站在线观看| 日本av手机在线免费观看| 90打野战视频偷拍视频| 岛国在线观看网站| 欧美精品一区二区免费开放| 国产精品影院久久| 国产激情久久老熟女| 亚洲国产精品999| 真人做人爱边吃奶动态| 国产免费福利视频在线观看| 午夜两性在线视频| 国产精品秋霞免费鲁丝片| 国产亚洲精品第一综合不卡| 免费av中文字幕在线| 精品少妇黑人巨大在线播放| 精品卡一卡二卡四卡免费| 黄频高清免费视频| 国产日韩欧美视频二区| 国产欧美日韩综合在线一区二区| 国产日韩欧美在线精品| 丰满迷人的少妇在线观看| 五月开心婷婷网| 国产精品久久久久成人av| 麻豆av在线久日| 欧美日韩亚洲高清精品| 欧美黄色淫秽网站| 一区福利在线观看| 人人妻人人添人人爽欧美一区卜| 丝袜美腿诱惑在线| 99热全是精品| 国产精品秋霞免费鲁丝片| 日韩视频在线欧美| 黄色a级毛片大全视频| 国产一区二区三区在线臀色熟女 | 人成视频在线观看免费观看| 老司机深夜福利视频在线观看 | 国产一级毛片在线| 99精品久久久久人妻精品| 婷婷成人精品国产| 国产又色又爽无遮挡免| 丁香六月欧美| 纵有疾风起免费观看全集完整版| 啦啦啦啦在线视频资源| 久久久久久亚洲精品国产蜜桃av| 两个人免费观看高清视频| 亚洲色图综合在线观看| 成年美女黄网站色视频大全免费| 亚洲欧美激情在线| 性少妇av在线| 国产欧美日韩一区二区三 | 久久亚洲国产成人精品v| 男男h啪啪无遮挡| 欧美在线一区亚洲| 欧美亚洲日本最大视频资源| 十分钟在线观看高清视频www| 国产免费av片在线观看野外av| 日日爽夜夜爽网站| 永久免费av网站大全| 久久久精品区二区三区| 十分钟在线观看高清视频www| 他把我摸到了高潮在线观看 | 黑人巨大精品欧美一区二区蜜桃| www.精华液| 国产欧美日韩精品亚洲av| 一本久久精品| 国产精品成人在线| 欧美av亚洲av综合av国产av| 青青草视频在线视频观看| 黄色视频,在线免费观看| 侵犯人妻中文字幕一二三四区| 99香蕉大伊视频| 亚洲欧美色中文字幕在线| 国产高清国产精品国产三级| 天天添夜夜摸| 欧美成人午夜精品| 极品人妻少妇av视频| av视频免费观看在线观看| 在线观看免费高清a一片| 国产精品久久久久久人妻精品电影 | 久久人人97超碰香蕉20202| 欧美另类亚洲清纯唯美| 大香蕉久久网| 岛国在线观看网站| 国产av又大| 久久久久久人人人人人| 国产国语露脸激情在线看| 天天躁夜夜躁狠狠躁躁| 成年人黄色毛片网站| 在线观看www视频免费| 汤姆久久久久久久影院中文字幕| 一区二区三区激情视频| 亚洲伊人久久精品综合| 亚洲成人国产一区在线观看| 亚洲专区字幕在线| 久久精品亚洲av国产电影网| 亚洲国产毛片av蜜桃av| 久久99热这里只频精品6学生| 性色av乱码一区二区三区2| 亚洲精品国产区一区二| 欧美一级毛片孕妇| 岛国毛片在线播放| 丝袜脚勾引网站| 视频区图区小说| 在线永久观看黄色视频| 91九色精品人成在线观看| 免费人妻精品一区二区三区视频| 一个人免费看片子| 看免费av毛片| 美女大奶头黄色视频| 韩国精品一区二区三区| 国产高清国产精品国产三级| 日韩欧美免费精品| 999精品在线视频| 一级片'在线观看视频| 黄色视频不卡| 国产精品1区2区在线观看. | 最近中文字幕2019免费版| 激情视频va一区二区三区| 十分钟在线观看高清视频www| 欧美黑人欧美精品刺激| 不卡av一区二区三区| 99精品久久久久人妻精品| 一区二区日韩欧美中文字幕| 国产免费一区二区三区四区乱码| 18在线观看网站| tube8黄色片| 国产99久久九九免费精品| 老鸭窝网址在线观看| 亚洲精品一区蜜桃| 9色porny在线观看| 高清黄色对白视频在线免费看| 18禁观看日本| 搡老乐熟女国产| 精品国产乱码久久久久久男人| 老汉色∧v一级毛片| 亚洲精品国产一区二区精华液| 色94色欧美一区二区| 亚洲熟女毛片儿| 欧美一级毛片孕妇| 国产精品.久久久| 精品少妇一区二区三区视频日本电影| 成人国语在线视频| 日本黄色日本黄色录像| 亚洲欧美精品综合一区二区三区| 女人精品久久久久毛片| 男女床上黄色一级片免费看| 亚洲国产精品成人久久小说| 日韩制服骚丝袜av| 日本精品一区二区三区蜜桃| 久久久久网色| 极品人妻少妇av视频| 国产精品二区激情视频| 久久久久网色| 亚洲伊人久久精品综合| 精品少妇黑人巨大在线播放| 九色亚洲精品在线播放| 亚洲免费av在线视频| 日韩三级视频一区二区三区| 一级黄色大片毛片| 老司机午夜福利在线观看视频 | 久久久久久久久久久久大奶| 极品少妇高潮喷水抽搐| 脱女人内裤的视频| 这个男人来自地球电影免费观看| 亚洲欧美色中文字幕在线| 99国产精品99久久久久| 国产黄频视频在线观看| 一区福利在线观看| 国产亚洲精品一区二区www | 亚洲,欧美精品.| 咕卡用的链子| 汤姆久久久久久久影院中文字幕|