• <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ā)時段
    欧美丝袜亚洲另类| 秋霞伦理黄片| 亚洲欧洲日产国产| 香蕉国产在线看| 22中文网久久字幕| 国产av精品麻豆| 免费高清在线观看日韩| 久久国产精品大桥未久av| av有码第一页| 黄色配什么色好看| 蜜桃在线观看..| 亚洲精华国产精华液的使用体验| 精品人妻在线不人妻| 亚洲少妇的诱惑av| 中文字幕亚洲精品专区| 中文字幕制服av| 欧美日韩亚洲高清精品| 99热6这里只有精品| 午夜激情av网站| 亚洲欧美清纯卡通| 久久精品久久精品一区二区三区| 青青草视频在线视频观看| 国产精品欧美亚洲77777| 中文精品一卡2卡3卡4更新| 七月丁香在线播放| 五月玫瑰六月丁香| 国产精品免费大片| 一边亲一边摸免费视频| 午夜久久久在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲国产av新网站| 国产一区二区在线观看av| 自拍欧美九色日韩亚洲蝌蚪91| 大码成人一级视频| 丰满乱子伦码专区| 人人澡人人妻人| 久久 成人 亚洲| 午夜福利影视在线免费观看| 精品国产国语对白av| 久久精品国产综合久久久 | 赤兔流量卡办理| 男女无遮挡免费网站观看| 久久这里有精品视频免费| 午夜免费男女啪啪视频观看| 亚洲精华国产精华液的使用体验| 亚洲精品美女久久久久99蜜臀 | 国产av一区二区精品久久| 考比视频在线观看| 久久久久精品人妻al黑| 亚洲成色77777| 寂寞人妻少妇视频99o| 两性夫妻黄色片 | 两个人看的免费小视频| 免费av中文字幕在线| 久久人人爽av亚洲精品天堂| 色视频在线一区二区三区| 毛片一级片免费看久久久久| 亚洲精品国产av成人精品| 久久久久网色| 妹子高潮喷水视频| 18禁动态无遮挡网站| 国产欧美日韩一区二区三区在线| 韩国精品一区二区三区 | 国产精品秋霞免费鲁丝片| 成人国语在线视频| √禁漫天堂资源中文www| 天堂8中文在线网| 啦啦啦中文免费视频观看日本| 亚洲精品av麻豆狂野| av电影中文网址| 国产成人av激情在线播放| 国产成人免费无遮挡视频| 精品少妇内射三级| 人妻人人澡人人爽人人| 久久久久网色| a级毛色黄片| 少妇 在线观看| 如何舔出高潮| 亚洲少妇的诱惑av| 成人综合一区亚洲| 熟女av电影| 亚洲av电影在线进入| 国产av码专区亚洲av| 香蕉丝袜av| 成人黄色视频免费在线看| 制服丝袜香蕉在线| 伊人亚洲综合成人网| 亚洲av国产av综合av卡| 久久精品国产自在天天线| 少妇的丰满在线观看| 多毛熟女@视频| 国产综合精华液| 热re99久久国产66热| 精品国产国语对白av| 午夜激情av网站| 黄色 视频免费看| 高清在线视频一区二区三区| 精品一区二区三卡| 久久婷婷青草| 亚洲性久久影院| 岛国毛片在线播放| 国产免费福利视频在线观看| 午夜av观看不卡| www.色视频.com| 久久久久人妻精品一区果冻| 精品少妇内射三级| 国产免费一区二区三区四区乱码| 国产亚洲精品久久久com| 麻豆精品久久久久久蜜桃| 国产高清三级在线| av视频免费观看在线观看| 国产探花极品一区二区| 夫妻性生交免费视频一级片| 99精国产麻豆久久婷婷| 视频在线观看一区二区三区| 69精品国产乱码久久久| 99久久中文字幕三级久久日本| 黑人高潮一二区| 狂野欧美激情性xxxx在线观看| av国产精品久久久久影院| av国产久精品久网站免费入址| 免费日韩欧美在线观看| 麻豆乱淫一区二区| 久久久久久久国产电影| 夜夜爽夜夜爽视频| 一区二区日韩欧美中文字幕 | 熟女电影av网| 精品人妻偷拍中文字幕| 亚洲av国产av综合av卡| 日本午夜av视频| 精品一区二区三区四区五区乱码 | 久久免费观看电影| 日本av免费视频播放| 精品国产乱码久久久久久小说| 日韩大片免费观看网站| av一本久久久久| 久久97久久精品| 99久久综合免费| 丝袜脚勾引网站| 久久久久久人人人人人| 欧美97在线视频| 国产一区二区三区综合在线观看 | 午夜福利视频精品| 美女国产高潮福利片在线看| 日韩一区二区三区影片| 亚洲av在线观看美女高潮| 成年动漫av网址| 国产探花极品一区二区| 你懂的网址亚洲精品在线观看| 国产精品国产三级专区第一集| 国产成人91sexporn| 色哟哟·www| 卡戴珊不雅视频在线播放| 下体分泌物呈黄色| 久热久热在线精品观看| 久久狼人影院| av有码第一页| 纵有疾风起免费观看全集完整版| 免费播放大片免费观看视频在线观看| 国产亚洲午夜精品一区二区久久| 国产1区2区3区精品| 亚洲精品第二区| 亚洲第一av免费看| 精品酒店卫生间| 丝袜人妻中文字幕| 精品一区二区免费观看| 久久久精品94久久精品| 男男h啪啪无遮挡| 久久久久国产精品人妻一区二区| 自线自在国产av| a级片在线免费高清观看视频| 肉色欧美久久久久久久蜜桃| 欧美精品一区二区免费开放| 日日撸夜夜添| 亚洲情色 制服丝袜| 亚洲五月色婷婷综合| 国产69精品久久久久777片| 日本91视频免费播放| 国产伦理片在线播放av一区| 日本-黄色视频高清免费观看| 欧美精品一区二区免费开放| 久久人人爽人人爽人人片va| 精品久久久精品久久久| 啦啦啦啦在线视频资源| 久久国产精品男人的天堂亚洲 | 丰满乱子伦码专区| 美国免费a级毛片| 国产精品久久久久久精品古装| 亚洲精品久久成人aⅴ小说| 久久久久久久久久人人人人人人| 欧美精品一区二区免费开放| av免费在线看不卡| 80岁老熟妇乱子伦牲交| 制服丝袜香蕉在线| 18禁裸乳无遮挡动漫免费视频| 中文字幕av电影在线播放| 99热这里只有是精品在线观看| 国产爽快片一区二区三区| 国产精品一区二区在线不卡| 亚洲欧美一区二区三区国产| 亚洲伊人色综图| 18在线观看网站| 男女边吃奶边做爰视频| 咕卡用的链子| 久久久国产欧美日韩av| videosex国产| 日韩中字成人| 午夜免费鲁丝| 久久精品夜色国产| 天天影视国产精品| 亚洲欧美一区二区三区黑人 | 日本色播在线视频| av天堂久久9| 亚洲精品日本国产第一区| 伊人久久国产一区二区| 精品酒店卫生间| 日韩伦理黄色片| xxxhd国产人妻xxx| 三级国产精品片| 在线观看免费视频网站a站| 国产福利在线免费观看视频| 亚洲欧美成人综合另类久久久| 国产综合精华液| 国产日韩欧美视频二区| 80岁老熟妇乱子伦牲交| 日本-黄色视频高清免费观看| 99久久中文字幕三级久久日本| 王馨瑶露胸无遮挡在线观看| 中文字幕制服av| 黄色一级大片看看| 国产熟女午夜一区二区三区| 视频中文字幕在线观看| 在线免费观看不下载黄p国产| 男人添女人高潮全过程视频| 永久网站在线| 久久久久视频综合| 夜夜骑夜夜射夜夜干| 亚洲少妇的诱惑av| 国产精品久久久久久久电影| 女人精品久久久久毛片| 又黄又粗又硬又大视频| 涩涩av久久男人的天堂| 91精品三级在线观看| 日韩av不卡免费在线播放| 黑人高潮一二区| 午夜福利影视在线免费观看| 丁香六月天网| 国产av码专区亚洲av| 精品久久久久久电影网| 午夜福利影视在线免费观看| av在线老鸭窝| 天天影视国产精品| 日韩av不卡免费在线播放| 视频中文字幕在线观看| 高清毛片免费看| 欧美xxⅹ黑人| 搡老乐熟女国产| 飞空精品影院首页| 精品人妻熟女毛片av久久网站| 日韩中字成人| 色视频在线一区二区三区| 国产成人欧美| 久久精品国产亚洲av涩爱| 亚洲美女视频黄频| 国语对白做爰xxxⅹ性视频网站| 精品一区二区三卡| 日韩av不卡免费在线播放| 中文字幕最新亚洲高清| 久久99热6这里只有精品| 伦精品一区二区三区| xxxhd国产人妻xxx| 久久ye,这里只有精品| 黄网站色视频无遮挡免费观看| 日韩,欧美,国产一区二区三区| 国产成人91sexporn| 丝袜脚勾引网站| 日韩在线高清观看一区二区三区| 日韩伦理黄色片| 免费观看无遮挡的男女| 久久精品国产鲁丝片午夜精品| 欧美xxⅹ黑人| 日本午夜av视频| 天堂中文最新版在线下载| 日韩中文字幕视频在线看片| 亚洲精品国产av成人精品| 欧美97在线视频| a级毛片在线看网站| av天堂久久9| 欧美亚洲日本最大视频资源| 精品人妻熟女毛片av久久网站| 男人操女人黄网站| 啦啦啦视频在线资源免费观看| 国产乱人偷精品视频| 99视频精品全部免费 在线| av免费在线看不卡| 咕卡用的链子| 少妇精品久久久久久久| 成人毛片a级毛片在线播放| 国产片内射在线| 丰满饥渴人妻一区二区三| 欧美日本中文国产一区发布| 日韩伦理黄色片| 少妇猛男粗大的猛烈进出视频| 久久精品国产鲁丝片午夜精品| 寂寞人妻少妇视频99o| 国产精品一区二区在线不卡| 妹子高潮喷水视频| 90打野战视频偷拍视频| www.色视频.com| 国产亚洲精品久久久com| 久热这里只有精品99| 国产精品国产三级专区第一集| 国产日韩一区二区三区精品不卡| 国产精品久久久久久精品古装| 国产在线免费精品| 国产成人91sexporn| 久久久久网色| 日韩成人伦理影院| 国产一区二区三区av在线| 在线观看免费高清a一片| 乱码一卡2卡4卡精品| 亚洲欧美一区二区三区黑人 | 2022亚洲国产成人精品| 久久人人爽av亚洲精品天堂| 18禁观看日本| 午夜福利在线观看免费完整高清在| 国产精品欧美亚洲77777| 超碰97精品在线观看| 夫妻性生交免费视频一级片| 我的女老师完整版在线观看| 亚洲精品国产av成人精品| 国产xxxxx性猛交| 男女边摸边吃奶| 91在线精品国自产拍蜜月| 乱人伦中国视频| 国产福利在线免费观看视频| 黄片播放在线免费| 高清毛片免费看| 国产成人精品无人区| 久久午夜综合久久蜜桃| 精品99又大又爽又粗少妇毛片| 久久人妻熟女aⅴ| 久久影院123| 久久久久久人妻| 满18在线观看网站| √禁漫天堂资源中文www| 久久99热6这里只有精品| 色5月婷婷丁香| 少妇人妻 视频| a级毛片黄视频| 精品少妇黑人巨大在线播放| av免费在线看不卡| 亚洲欧美一区二区三区国产| 在线天堂中文资源库| 久久久久视频综合| 少妇熟女欧美另类| 久久久久国产网址| av片东京热男人的天堂| 五月天丁香电影| 免费大片黄手机在线观看| 国产一区亚洲一区在线观看| 一边摸一边做爽爽视频免费| videosex国产| 夫妻午夜视频| 在线亚洲精品国产二区图片欧美| 免费观看在线日韩| 亚洲欧洲精品一区二区精品久久久 | av电影中文网址| 97超碰精品成人国产| 亚洲国产精品专区欧美| 欧美xxxx性猛交bbbb| av有码第一页| 少妇的逼好多水| 女人被躁到高潮嗷嗷叫费观| 成年美女黄网站色视频大全免费| 国国产精品蜜臀av免费| 久久久久久人人人人人| 亚洲欧美一区二区三区黑人 | 制服人妻中文乱码| 亚洲国产欧美日韩在线播放| 国产精品久久久久久久电影| 丰满乱子伦码专区| 日日啪夜夜爽| 两性夫妻黄色片 | 一级爰片在线观看| 亚洲美女黄色视频免费看| 最近手机中文字幕大全| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产又色又爽无遮挡免| 日本午夜av视频| 国产片特级美女逼逼视频| 久久久久精品人妻al黑| 亚洲成国产人片在线观看| 国产亚洲一区二区精品| 九九在线视频观看精品| 搡老乐熟女国产| 极品少妇高潮喷水抽搐| 2022亚洲国产成人精品| 亚洲综合精品二区| 另类亚洲欧美激情| 五月伊人婷婷丁香| 黄色视频在线播放观看不卡| 女人久久www免费人成看片| 一边亲一边摸免费视频| 欧美成人精品欧美一级黄| 交换朋友夫妻互换小说| 日韩成人av中文字幕在线观看| 久久久国产欧美日韩av| 男女高潮啪啪啪动态图| 三上悠亚av全集在线观看| 精品人妻偷拍中文字幕| 在线观看免费高清a一片| 午夜免费观看性视频| 精品一区二区三区四区五区乱码 | 国产精品嫩草影院av在线观看| 99国产精品免费福利视频| 精品午夜福利在线看| 三上悠亚av全集在线观看| 亚洲精品视频女| 国产精品久久久久成人av| 日本黄色日本黄色录像| 在线观看三级黄色| 国产深夜福利视频在线观看| 欧美日韩成人在线一区二区| 亚洲综合精品二区| 亚洲精品美女久久久久99蜜臀 | 一级毛片我不卡| 欧美日本中文国产一区发布| 国产亚洲最大av| tube8黄色片| 国产一区二区在线观看日韩| 一二三四在线观看免费中文在 | 国产精品人妻久久久影院| 成年av动漫网址| 久久99蜜桃精品久久| 五月天丁香电影| 免费观看av网站的网址| 一区二区三区四区激情视频| 午夜免费观看性视频| 最近最新中文字幕免费大全7| 美女视频免费永久观看网站| 亚洲av欧美aⅴ国产| 最后的刺客免费高清国语| 国产成人一区二区在线| 亚洲精品色激情综合| 亚洲图色成人| 美女主播在线视频| 侵犯人妻中文字幕一二三四区| 最近中文字幕2019免费版| 亚洲,一卡二卡三卡| 亚洲人成网站在线观看播放| 精品久久蜜臀av无| 国产乱人偷精品视频| 看非洲黑人一级黄片| 妹子高潮喷水视频| 男女边摸边吃奶| 日本欧美国产在线视频| 成人免费观看视频高清| 9色porny在线观看| 久久久精品免费免费高清| 国产麻豆69| 飞空精品影院首页| 只有这里有精品99| 久久久久国产精品人妻一区二区| 亚洲欧美精品自产自拍| 精品少妇内射三级| 九草在线视频观看| 欧美性感艳星| 精品一品国产午夜福利视频| 中国三级夫妇交换| 老熟女久久久| 久久精品国产自在天天线| 热re99久久精品国产66热6| 国产永久视频网站| 亚洲国产精品成人久久小说| 一级毛片黄色毛片免费观看视频| 18禁在线无遮挡免费观看视频| 美女国产视频在线观看| 中文字幕av电影在线播放| 色5月婷婷丁香| 国产精品欧美亚洲77777| 久久狼人影院| av又黄又爽大尺度在线免费看| 高清在线视频一区二区三区| 草草在线视频免费看| 久久久久久久大尺度免费视频| 两性夫妻黄色片 | 欧美日本中文国产一区发布| 精品第一国产精品| 国产av一区二区精品久久| 夫妻午夜视频| 亚洲丝袜综合中文字幕| 国产精品免费大片| 久久狼人影院| 有码 亚洲区| 韩国精品一区二区三区 | 国产成人免费无遮挡视频| 黄色 视频免费看| 国产男女超爽视频在线观看| 日韩伦理黄色片| 看免费av毛片| 狠狠精品人妻久久久久久综合| 丝袜喷水一区| 国产白丝娇喘喷水9色精品| 成人国产麻豆网| 欧美国产精品va在线观看不卡| 欧美激情极品国产一区二区三区 | 亚洲国产精品999| 久久热在线av| 亚洲欧美成人精品一区二区| 国产精品三级大全| 久久婷婷青草| 免费av中文字幕在线| 午夜免费男女啪啪视频观看| 国产精品欧美亚洲77777| 精品一品国产午夜福利视频| 国产激情久久老熟女| 日韩电影二区| 一边亲一边摸免费视频| 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 免费人成在线观看视频色| 中文精品一卡2卡3卡4更新| 国产高清三级在线| 日韩中文字幕视频在线看片| 久久久久久久亚洲中文字幕| 久久久久久伊人网av| 亚洲少妇的诱惑av| 亚洲国产精品专区欧美| 新久久久久国产一级毛片| 久久久久久伊人网av| 伊人久久国产一区二区| 国产成人免费观看mmmm| 精品亚洲乱码少妇综合久久| 少妇人妻久久综合中文| 丁香六月天网| xxx大片免费视频| 免费人成在线观看视频色| 国产精品免费大片| 国产淫语在线视频| 欧美另类一区| 国产精品人妻久久久久久| 亚洲av电影在线观看一区二区三区| 18禁国产床啪视频网站| 青春草亚洲视频在线观看| 中国美白少妇内射xxxbb| 黑人欧美特级aaaaaa片| 午夜福利,免费看| 欧美日韩视频高清一区二区三区二| 免费在线观看完整版高清| 午夜免费鲁丝| 国产一区二区激情短视频 | 精品酒店卫生间| 黄色怎么调成土黄色| 日韩,欧美,国产一区二区三区| 亚洲精品av麻豆狂野| 成人免费观看视频高清| av卡一久久| 美女脱内裤让男人舔精品视频| av电影中文网址| 狠狠婷婷综合久久久久久88av| 国产亚洲精品久久久com| 永久免费av网站大全| 少妇 在线观看| 久久99蜜桃精品久久| 黄色 视频免费看| 久久久久精品久久久久真实原创| 高清毛片免费看| 日本wwww免费看| 男女午夜视频在线观看 | 少妇精品久久久久久久| 一本—道久久a久久精品蜜桃钙片| 欧美精品高潮呻吟av久久| 自拍欧美九色日韩亚洲蝌蚪91| 日韩精品免费视频一区二区三区 | 国产黄色视频一区二区在线观看| 91午夜精品亚洲一区二区三区| 国产成人aa在线观看| 只有这里有精品99| 久久久久久人人人人人| 亚洲综合色惰| 巨乳人妻的诱惑在线观看| 男女高潮啪啪啪动态图| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄片无遮挡物在线观看| 女人被躁到高潮嗷嗷叫费观| 极品人妻少妇av视频| 国产精品麻豆人妻色哟哟久久| 中国三级夫妇交换| 青春草视频在线免费观看| 精品一区在线观看国产| 亚洲精品美女久久久久99蜜臀 | 三级国产精品片| 日韩中字成人| 丁香六月天网| 久久人人爽人人爽人人片va| 亚洲精华国产精华液的使用体验| av国产精品久久久久影院| 黑丝袜美女国产一区| 下体分泌物呈黄色| 成人手机av| 日韩制服丝袜自拍偷拍| 99久久中文字幕三级久久日本| 精品视频人人做人人爽| 成人漫画全彩无遮挡| 国产白丝娇喘喷水9色精品| 亚洲欧美清纯卡通| 永久免费av网站大全| 青春草国产在线视频| 国产一区二区在线观看日韩| 久久久a久久爽久久v久久| 午夜久久久在线观看| 国产成人一区二区在线| 少妇猛男粗大的猛烈进出视频| 一级毛片黄色毛片免费观看视频|