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

    鎖相環(huán)失步對(duì)同步發(fā)電機(jī)暫態(tài)功角穩(wěn)定性的影響機(jī)理分析

    2022-11-28 11:03:04李劍澤蘇晨博王潔聰劉崇茹
    電力系統(tǒng)自動(dòng)化 2022年22期
    關(guān)鍵詞:折點(diǎn)功角鎖相環(huán)

    黎 曉,李劍澤,蘇晨博,王潔聰,劉崇茹,鄭 樂

    (1. 新能源電力系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室(華北電力大學(xué)),北京市 102206;2. 中國(guó)電力科學(xué)研究院有限公司,北京市 100192)

    0 引言

    為實(shí)現(xiàn)中國(guó)“雙碳”目標(biāo),據(jù)估計(jì)到2060 年新能源將逐漸成為電力供應(yīng)主體,但傳統(tǒng)的同步電源依然占有一定比例,電網(wǎng)仍通過交流同步機(jī)制實(shí)現(xiàn)功率傳輸[1]。新型電力系統(tǒng)中,多種設(shè)備的多時(shí)間尺度特性給同步發(fā)電機(jī)暫態(tài)功角穩(wěn)定性問題帶來新的挑戰(zhàn)。 目前,電壓源型變流器(voltage source converter,VSC)是一種常見的新能源并網(wǎng)設(shè)備,其中已被廣泛使用的跟網(wǎng)型變流器(grid-following converter)依靠鎖相環(huán)(phase-locked loop,PLL)與電網(wǎng)保持同步[2]。因此,研究變流器的鎖相環(huán)和傳統(tǒng)同步發(fā)電機(jī)互聯(lián)的暫態(tài)穩(wěn)定性具有重要意義。

    在小擾動(dòng)方面,已有較多文獻(xiàn)研究了鎖相環(huán)與同步發(fā)電機(jī)功角穩(wěn)定的關(guān)系,相關(guān)的研究結(jié)論已比較成熟。例如,弱電網(wǎng)情況下鎖相環(huán)與同步發(fā)電機(jī)轉(zhuǎn)子軸承的耦合作用可能引發(fā)次同步振蕩[3-4]。鎖相環(huán)參數(shù)不合理時(shí)可能引起同步發(fā)電機(jī)低頻振蕩[5-8],優(yōu)化鎖相環(huán)參數(shù)可以為同步發(fā)電機(jī)提供正阻尼[9]。

    在大擾動(dòng)方面,基于不同的模型和方法得到的結(jié)論并不完全一致。已有的針對(duì)鎖相環(huán)與同步發(fā)電機(jī)功角的相互影響的研究認(rèn)為:1)當(dāng)鎖相環(huán)穩(wěn)定時(shí),變流器對(duì)電網(wǎng)主要表現(xiàn)出變流器外環(huán)控制的特征,可以近似等值為恒功率或恒阻抗[10]。較多文獻(xiàn)討論了這種情況下新能源電源經(jīng)變流器接入電網(wǎng)的滲透率[10-12]、接入點(diǎn)位置[13-15]及功率控制方式[16-17]等對(duì)同步發(fā)電機(jī)功角穩(wěn)定性的影響。2)大擾動(dòng)下鎖相環(huán)失步時(shí),認(rèn)為鎖相環(huán)動(dòng)態(tài)過程的時(shí)間尺度明顯小于發(fā)電機(jī)的轉(zhuǎn)子時(shí)間尺度,因此,將鎖相環(huán)簡(jiǎn)化為代數(shù)方程來討論同步發(fā)電機(jī)的功角穩(wěn)定性問題[18]。文獻(xiàn)[18]指出,含變流器電網(wǎng)的同步發(fā)電機(jī)暫態(tài)功角失穩(wěn)機(jī)理與傳統(tǒng)電網(wǎng)相比,存在新的失穩(wěn)約束。

    然而,當(dāng)變流器公共連接點(diǎn)(point of common coupling,PCC)處電壓跌落[19-21],或者電網(wǎng)強(qiáng)度較弱[22-23]、故障后PCC 處相位跳變[24]等情況發(fā)生時(shí),鎖相環(huán)可能失去與電網(wǎng)保持同步的能力[18-19]。此時(shí),將換流器簡(jiǎn)化為恒功率源或恒阻抗無法反映出鎖相環(huán)失步情況下對(duì)同步發(fā)電機(jī)的影響。當(dāng)鎖相環(huán)在大擾動(dòng)后的暫態(tài)過程中失去平衡點(diǎn)時(shí),將鎖相環(huán)簡(jiǎn)化為代數(shù)方程的條件不再成立?,F(xiàn)有的研究結(jié)論,即變流器失去平衡點(diǎn)后導(dǎo)致同步發(fā)電機(jī)失去同步的條件值得進(jìn)一步研究。

    本文針對(duì)上述問題,深入研究鎖相環(huán)動(dòng)態(tài)與同步發(fā)電機(jī)的交互影響機(jī)理,重點(diǎn)考慮故障后系統(tǒng)恢復(fù)過程中兩者交互引起的同步穩(wěn)定問題,而非故障直接引起的失穩(wěn)情況。本文的主要貢獻(xiàn)為:

    1)建立了鎖相環(huán)與同步發(fā)電機(jī)功角的動(dòng)力學(xué)方程,分析了鎖相環(huán)與同步發(fā)電機(jī)功角耦合的條件,揭示了大量變流器接入電網(wǎng)時(shí),同步發(fā)電機(jī)功角穩(wěn)定邊界由不穩(wěn)定平衡點(diǎn)約束轉(zhuǎn)為鎖相環(huán)約束的機(jī)理。基于幾何奇異攝動(dòng)理論,從理論層面解釋了鎖相環(huán)失去平衡點(diǎn)時(shí),系統(tǒng)軌跡從同步發(fā)電機(jī)動(dòng)態(tài)主導(dǎo)過渡到鎖相環(huán)動(dòng)態(tài)主導(dǎo)的動(dòng)力學(xué)特征。

    2)通過發(fā)電機(jī)轉(zhuǎn)速的解析,分析了鎖相環(huán)簡(jiǎn)化為代數(shù)方程在功角穩(wěn)定性分析中的局限性,探討了同步發(fā)電機(jī)功角穩(wěn)定受鎖相環(huán)制約的條件,為大規(guī)模新能源饋入電網(wǎng)的暫態(tài)穩(wěn)定域分析提供了理論支撐。

    1 鎖相環(huán)與發(fā)電機(jī)功角動(dòng)態(tài)交互場(chǎng)景描述

    系統(tǒng)中的電力電子接口控制環(huán)節(jié)與同步發(fā)電機(jī)的轉(zhuǎn)子運(yùn)動(dòng)具有明顯的時(shí)間尺度差異。如圖1 所示,內(nèi)環(huán)電流控制的帶寬在100 Hz 及以上,鎖相環(huán)控制的帶寬一般在10~70 Hz,而發(fā)電機(jī)轉(zhuǎn)子的慣性時(shí)間常數(shù)通常在秒級(jí)。因此,按照不同時(shí)間尺度對(duì)系統(tǒng)的動(dòng)態(tài)過程進(jìn)行劃分,能夠簡(jiǎn)化分析過程。圖中:為變流器的電流指令值,在文中用Id和Iq表示;θpll為鎖相環(huán)的輸出角;Z1、Z2、Z3為網(wǎng)絡(luò)阻抗;ZL為PCC 處并聯(lián)阻抗;為大系統(tǒng)母線的電壓相量;Eg和δg分別為發(fā)電機(jī)內(nèi)電勢(shì)幅值和功角。當(dāng)研究鎖相環(huán)的同步穩(wěn)定性時(shí),認(rèn)為同步發(fā)電機(jī)轉(zhuǎn)子的相位不發(fā)生變化;當(dāng)研究同步發(fā)電機(jī)的功角穩(wěn)定性時(shí),假設(shè)鎖相環(huán)在合理的控制參數(shù)下能夠始終保持同步,此時(shí),變流器表現(xiàn)出外環(huán)功率控制或低壓穿越控制的特性。

    圖1 變流器接入交流系統(tǒng)示意圖Fig.1 Schematic diagram of converter connected to AC system

    然而,當(dāng)故障持續(xù)時(shí)間較長(zhǎng)導(dǎo)致發(fā)電機(jī)的功角擺開較大時(shí),振蕩中心區(qū)域的母線電壓大幅下降,進(jìn)一步,可能會(huì)導(dǎo)致其所連接的變流器的鎖相環(huán)因失去平衡點(diǎn)而失去同步。此時(shí),鎖相環(huán)失穩(wěn)主要受發(fā)電機(jī)功角而不是鎖相環(huán)控制參數(shù)的影響。同時(shí),鎖相環(huán)失去穩(wěn)定之后,變流器輸出的功率出現(xiàn)振蕩,也會(huì)影響同步發(fā)電機(jī)功角的恢復(fù)過程。因此,在含電壓源型變流器的電網(wǎng)中,同步發(fā)電機(jī)的功角穩(wěn)定性可能受到鎖相環(huán)的制約[18]。

    鎖相環(huán)對(duì)發(fā)電機(jī)功角穩(wěn)定域的約束條件與傳統(tǒng)電力系統(tǒng)中恒功率性質(zhì)的負(fù)荷對(duì)同步發(fā)電機(jī)的功角穩(wěn)定域的約束條件,具有一定的相似性,因?yàn)閮烧叨际峭ㄟ^影響交流網(wǎng)絡(luò)的功率傳輸能力,間接地改變了發(fā)電機(jī)的功角穩(wěn)定性[25-26]。但由于存在前文所描述的時(shí)間尺度差異,兩者表現(xiàn)出不同的現(xiàn)象。

    簡(jiǎn)而言之,變流器平衡點(diǎn)存在的條件,影響了同步發(fā)電機(jī)功角穩(wěn)定域的邊界條件。這既是鎖相環(huán)動(dòng)態(tài)與發(fā)電機(jī)功角之間多時(shí)間尺度特性的表現(xiàn),也是鎖相環(huán)通過交流網(wǎng)絡(luò)與同步發(fā)電機(jī)耦合的特征。

    2 變流器-同步發(fā)電機(jī)系統(tǒng)的動(dòng)力學(xué)行為

    2.1 鎖相環(huán)與發(fā)電機(jī)功角的雙時(shí)間尺度建模

    發(fā)電機(jī)的簡(jiǎn)化二階方程如式(1)所示。

    式中:ω0=(100π)rad/s 為工頻50 Hz 下的角速度;TJ為標(biāo)幺化的發(fā)電機(jī)慣性時(shí)間常數(shù);Pm、Pg和D分別為發(fā)電機(jī)的機(jī)械功率、電磁功率和阻尼系數(shù);ωg為δg的角頻率標(biāo)幺值,穩(wěn)態(tài)下等于1 p.u.。

    對(duì)于變流器和同步發(fā)電機(jī),本文為突出機(jī)理描述,有以下假設(shè):1)忽略電流內(nèi)環(huán)控制,假設(shè)電流內(nèi)環(huán)控制的帶寬足夠大,使得內(nèi)環(huán)控制與鎖相環(huán)近似解耦[27];2)忽略外環(huán)控制和直流電壓控制,假設(shè)故障及故障后恢復(fù)期間外環(huán)控制達(dá)到限幅或者故障穿越 控 制 啟 動(dòng)[18,20-23,28-30],變 流 器 表 現(xiàn) 出 電 流 源 的 特性;3)同步發(fā)電機(jī)暫態(tài)電抗后的內(nèi)電勢(shì)恒定。

    根據(jù)三相同步坐標(biāo)系的鎖相環(huán)(synchronous reference frame phase-locked loop,SRF-PLL)的 控制框圖(見附錄A 圖A1),建立鎖相環(huán)狀態(tài)方程如下:

    式中:θp為鎖相環(huán)的直接輸出角度θpll與平衡節(jié)點(diǎn)的相角差,是一個(gè)相對(duì)角,假設(shè)平衡節(jié)點(diǎn)相角的轉(zhuǎn)速等于ω0,ω0和θp均 為 有 名 值;xp為 鎖 相 環(huán) 比 例-積 分(PI)控制器中積分的輸出值;Vq為PCC 處電壓相量在q軸上的分量;ωp為角頻率標(biāo)幺值,穩(wěn)態(tài)下等于1 p.u.;kp為 鎖 相 環(huán)PI 環(huán) 節(jié) 的 比 例 系 數(shù);ki為 鎖 相 環(huán)PI 環(huán)節(jié)的積分增益。

    將圖1(a)中的交流電網(wǎng)用一個(gè)含等值大系統(tǒng)的三節(jié)點(diǎn)網(wǎng)絡(luò)代替,本文通過如圖1(b)所示的簡(jiǎn)化系統(tǒng)來進(jìn)行機(jī)理研究。其中,等值大系統(tǒng)的端電壓相角為0°,同步發(fā)電機(jī)和變流器分別代表發(fā)電機(jī)的一個(gè)集合和變流器的一個(gè)集合[11,31]。在風(fēng)機(jī)的穩(wěn)定性研究中,也常常使用單臺(tái)等值風(fēng)機(jī)代替整個(gè)風(fēng)電場(chǎng)以簡(jiǎn)化研究[32]。Pg和Vq可以寫成θp和δg的函數(shù),見附錄B,組合式(1)和式(2)即可得到系統(tǒng)的狀態(tài)方程。假設(shè)鎖相環(huán)的響應(yīng)速度遠(yuǎn)快于發(fā)電機(jī)轉(zhuǎn)子的響應(yīng)速度,系統(tǒng)具有雙時(shí)間尺度特性。通過引入無量綱參數(shù)ε=1/(kpTJ)和坐標(biāo)變換,可以將系統(tǒng)的狀態(tài)方程改寫成如式(3)所示的奇異攝動(dòng)的形式。

    式中:yp=xpkp,為鎖相環(huán)PI 控制器中積分的輸出值xp與kp的乘積;λ=t/TJ為時(shí)間變量。式(3)忽略了頻率變化對(duì)線路阻抗參數(shù)的影響。

    基于奇異攝動(dòng)理論[33],原模型解耦為一個(gè)關(guān)于發(fā)電機(jī)轉(zhuǎn)子動(dòng)態(tài)的慢系統(tǒng)(式(4))與一個(gè)關(guān)于鎖相環(huán)動(dòng)態(tài)的快系統(tǒng)(式(5))。

    針對(duì)圖1(b)所示簡(jiǎn)化網(wǎng)絡(luò)得到式(3)、式(4)和式(5)的推導(dǎo)過程,詳見附錄B。

    注意到快系統(tǒng)(式(5))的平衡點(diǎn)也是慢系統(tǒng)(式(4))中代數(shù)方程Vq=0 的解。因此,在快系統(tǒng)(式(5))穩(wěn)定的情況下,系統(tǒng)方程(式(3))的軌跡會(huì)很快趨近Vq=0 的解,如附錄C 圖C1 所示。此時(shí),慢系統(tǒng)與快系統(tǒng)可以近似解耦,慢系統(tǒng)能夠描述同步發(fā)電機(jī)的暫態(tài)功角穩(wěn)定性,與傳統(tǒng)交流電網(wǎng)的功角穩(wěn)定性沒有本質(zhì)上的區(qū)別。

    2.2 鎖相環(huán)與發(fā)電機(jī)功角耦合條件的分析

    2.2.1 快系統(tǒng)穩(wěn)定的條件

    快系統(tǒng)在擾動(dòng)后能夠恢復(fù)穩(wěn)定需要兩個(gè)條件:1)快系統(tǒng)存在穩(wěn)定的平衡點(diǎn);2)故障切除后軌跡位于快系統(tǒng)的穩(wěn)定域內(nèi)。

    對(duì)于快系統(tǒng)第2 個(gè)穩(wěn)定條件,已有文獻(xiàn)分析鎖相環(huán)的穩(wěn)定域[34]。因?yàn)楣堑淖兓俣蕊@著慢于鎖相環(huán)的調(diào)節(jié)速度,且功角較小時(shí)鎖相環(huán)平衡點(diǎn)距離穩(wěn)定邊界有一定的穩(wěn)定裕度,功角增大造成的鎖相環(huán)收斂域減小難以造成鎖相環(huán)失穩(wěn)。也就是說,在合適的鎖相環(huán)參數(shù)下,如果故障切除后鎖相環(huán)能夠很快恢復(fù)穩(wěn)定,那么鎖相環(huán)可以跟蹤上發(fā)電機(jī)的功角變化。

    忽略線路電導(dǎo)和并聯(lián)電導(dǎo),基于附錄B 中Vq表達(dá)式(B4),Vq=0 等價(jià)于式(6)。

    其中:

    式中:B2和B3分別為阻抗Z2和Z3支路的電納。已將同步發(fā)電機(jī)等值內(nèi)電抗通過電路星形-三角形變換歸算到網(wǎng)絡(luò)支路中;?I=arctan(-Iq/Id),Id和Iq分別為變流器的d、q軸電流;I3為變流器輸出電流的幅值。

    由于aδ(δg)最小值在δg=π 處取得,當(dāng)aδ(π)<I3cos?I時(shí),發(fā)電機(jī)功角增大到π 的過程中就會(huì)使得代數(shù)方程Vq=0 失去平衡點(diǎn)。此時(shí),發(fā)電機(jī)功角增大會(huì)使得快系統(tǒng)的第1 個(gè)穩(wěn)定條件不能滿足。從網(wǎng)絡(luò)結(jié)構(gòu)上表現(xiàn)為,B2與B3Eg越接近,鎖相環(huán)越容易失去平衡點(diǎn);從變流器容量上表現(xiàn)為,變流器輸出電流越大,鎖相環(huán)越容易失去平衡點(diǎn)。

    當(dāng)發(fā)電機(jī)功角增大,代數(shù)方程Vq=0 剛好失去平衡點(diǎn)時(shí),將這一點(diǎn)記作折點(diǎn)F[35],也稱為鞍結(jié)型分岔點(diǎn)和奇異點(diǎn)[25]。此時(shí),子系統(tǒng)(式(4))無法描述交直流電網(wǎng)的同步穩(wěn)定性,鎖相環(huán)與同步發(fā)電機(jī)的動(dòng)態(tài)過程耦合,鎖相環(huán)的動(dòng)態(tài)過程就不能完全被忽略。本文2.3 節(jié)詳細(xì)分析了這種條件下兩者的交互過程。

    2.2.2 慢系統(tǒng)的穩(wěn)定邊界由折點(diǎn)決定的條件

    功角穩(wěn)定邊界的評(píng)估大多基于不穩(wěn)定平衡點(diǎn)處的減速面積、暫態(tài)能量或者其穩(wěn)定流形。當(dāng)同步發(fā)電機(jī)功角在擾動(dòng)后首先經(jīng)過不穩(wěn)定平衡點(diǎn)而不是折點(diǎn)時(shí),基于不穩(wěn)定平衡點(diǎn)的功角穩(wěn)定性評(píng)估方法仍然適用。但是,當(dāng)同步發(fā)電機(jī)功角增大過程中先經(jīng)過折點(diǎn)時(shí),折點(diǎn)將會(huì)影響功角穩(wěn)定域的邊界。

    據(jù)統(tǒng)計(jì),2017年我國(guó)城鎮(zhèn)化率已達(dá)58.52%,農(nóng)業(yè)轉(zhuǎn)移到城鎮(zhèn)的人口日益成為城市建設(shè)的重要參與者,也是培育農(nóng)業(yè)農(nóng)村發(fā)展新動(dòng)能、推動(dòng)農(nóng)業(yè)農(nóng)村現(xiàn)代化的重要探索者。實(shí)施鄉(xiāng)村振興戰(zhàn)略,加快推進(jìn)農(nóng)業(yè)農(nóng)村現(xiàn)代化,不僅能夠推動(dòng)廣大農(nóng)民在共建共享發(fā)展中實(shí)現(xiàn)更多獲得感、幸福感和安全感,而且有利于更好踐行“以人民為中心”的發(fā)展思想。[3]

    本節(jié)基于發(fā)電機(jī)電磁功率Pg的解析表達(dá)式,分析了同步發(fā)電機(jī)不穩(wěn)定平衡點(diǎn)與折點(diǎn)的位置關(guān)系,以及同步發(fā)電機(jī)機(jī)械功率Pm和變流器電流大小I3對(duì)這種位置關(guān)系的影響,揭示了折點(diǎn)影響功角穩(wěn)定邊界的一些運(yùn)行條件。

    基于附錄B 的推導(dǎo)得到消去θp的發(fā)電機(jī)電磁功率Pg的表達(dá)式如式(9)所示。

    其中:

    式中:b1為功率表達(dá)式中正弦分量的系數(shù);B1為阻抗Z1支路的電納;BL為PCC 處的并聯(lián)電納;Pg2(δg)的表達(dá)式詳見附錄B。在δg取值范圍內(nèi),Pg2(δg)變化較小,可以簡(jiǎn)化為常數(shù)。

    因此,鎖相環(huán)保持在其穩(wěn)定平衡點(diǎn)的情況下,Pg關(guān)于δg變化的曲線可以更簡(jiǎn)潔地表示為:

    其中:

    δg取值范圍分以下兩種情況:

    情況1:當(dāng)aδ(π)≥I3cos?I時(shí),δg的取值范圍為實(shí)數(shù)域,這種情況不存在折點(diǎn)。

    情況2:當(dāng)aδ(π)<I3cos?I時(shí),δg的取值范 圍為[-δg,m,δg,m],取值范圍的兩端為折點(diǎn)的位置。δg,m的表達(dá)式如下:

    在δg取值范圍的情況2 下,折點(diǎn)和不穩(wěn)定平衡點(diǎn)的相對(duì)位置受變流器輸出電流I3與同步發(fā)電機(jī)機(jī)械功率Pm的影響。由式(11)可知,Pg-δg曲線接近正弦函數(shù)。Pg=Pm>0 的解中,功角較小的解稱為穩(wěn)定平衡點(diǎn),其他解稱為不穩(wěn)定平衡點(diǎn)。由δg的取值范圍可知,I3越大則δg取值范圍越小。因此,I3增大可以使得在δg的取值范圍之內(nèi)只存在穩(wěn)定平衡點(diǎn),不穩(wěn)定平衡點(diǎn)在功角δg的取值范圍之外。同時(shí),較小的Pm也會(huì)使得不穩(wěn)定平衡點(diǎn)在δg的取值范圍之外。

    以附錄D 表D1 的網(wǎng)絡(luò)參數(shù)1 為例驗(yàn)證上述分析。同步發(fā)電機(jī)等值內(nèi)阻抗為0.05 p.u.,內(nèi)電勢(shì)Eg=1.1 p.u.,電磁功率曲線如圖2 所示。結(jié)果表明,式(11)能夠很好地近似式(9),驗(yàn)證了利用式(11)研究Pg性質(zhì)的正確性。比較圖2 中不同I3條件下的Pg曲線,以及Pg=Pm的解的情況。當(dāng)I3=1 p.u.或3 p.u.時(shí),Pg=4 p.u.有兩個(gè)解:穩(wěn)定平衡點(diǎn)和不穩(wěn)定平衡點(diǎn)。而當(dāng)I3增加到5 p.u.時(shí),Pg=4 p.u.僅有一個(gè)穩(wěn)定平衡點(diǎn)的解,不存在不穩(wěn)定平衡點(diǎn),說明較大的I3使得不穩(wěn)定平衡點(diǎn)消失。比較不同Pm條件下的解。對(duì)于Pg=2 p.u.的解,當(dāng)I3=3 p.u.或5 p.u.時(shí),δg都只有一個(gè)穩(wěn)定平衡點(diǎn)。與Pg=4 p.u.相比,較小的Pm也會(huì)使得不穩(wěn)定平衡點(diǎn)消失。

    圖2 同步發(fā)電機(jī)電磁功率Pg-功角δg曲線Fig.2 Curves of electromagnetic power Pg versus power angle δg of synchronous generators

    因此,較大的變流器輸出電流I3與較小的同步發(fā)電機(jī)機(jī)械功率Pm會(huì)使得慢系統(tǒng)的穩(wěn)定邊界由折點(diǎn)決定而不是不穩(wěn)定平衡點(diǎn)。這也是后文分析的鎖相環(huán)與同步發(fā)電機(jī)動(dòng)態(tài)交互的一個(gè)條件,這個(gè)條件也正是大規(guī)模新能源接入電網(wǎng)的特點(diǎn)。

    2.3 慢動(dòng)態(tài)到快動(dòng)態(tài)過程的定性分析

    針對(duì)前文所描述的鎖相環(huán)與發(fā)電機(jī)功角動(dòng)態(tài)交互的場(chǎng)景,基于幾何奇異攝動(dòng)理論[36],附錄C 中的定理C1 對(duì)鎖相環(huán)與同步發(fā)電機(jī)動(dòng)態(tài)耦合的過程進(jìn)行了理論解釋。

    由該定理可知,如果發(fā)電機(jī)功角δg在逼近折點(diǎn)F之前速度未減為0,那么系統(tǒng)軌跡將越過折點(diǎn)F,且在F點(diǎn)足夠小的領(lǐng)域內(nèi),系統(tǒng)的動(dòng)態(tài)過程將從式(4)所描述的慢動(dòng)態(tài)過程跳躍到式(5)所描述的快動(dòng)態(tài)過程。也即,系統(tǒng)軌跡Sa,ε沿著Vq=0 的穩(wěn)定分支Sa越過折點(diǎn)F之后,鎖相環(huán)輸出角度θp幾乎總是以遠(yuǎn)快于發(fā)電機(jī)功角變化的速度增大,而發(fā)電機(jī)功角δg變化較小。下一章將進(jìn)一步通過分析折點(diǎn)F處的一些性質(zhì),得到鎖相環(huán)對(duì)發(fā)電機(jī)功角穩(wěn)定的影響。

    3 鎖相環(huán)失步對(duì)同步發(fā)電機(jī)動(dòng)態(tài)的影響機(jī)理

    3.1 折點(diǎn)附近同步發(fā)電機(jī)轉(zhuǎn)速的解析分析

    在折點(diǎn)處進(jìn)行泰勒展開,以折點(diǎn)領(lǐng)域內(nèi)的一點(diǎn)作為t=0 時(shí)刻,得到發(fā)電機(jī)轉(zhuǎn)速在折點(diǎn)的一個(gè)很小領(lǐng)域內(nèi)的表達(dá)式如式(14)所示。

    式中:系數(shù)d1等于t=0 時(shí)刻2TJdωg/dt的值;d2等于t=0 時(shí)刻2εTJ(?Pg/?θp)(dθp/dt)的值;等于t=0時(shí)刻的發(fā)電機(jī)轉(zhuǎn)速偏差。推導(dǎo)過程與d1和d2的系數(shù)表達(dá)式見附錄E。

    根據(jù)式(14),軌跡穿過附錄C 圖C2(a)所示的折點(diǎn)F附近區(qū)域之后進(jìn)入折點(diǎn)F的右側(cè)區(qū)域時(shí),有以下兩種情況:

    情況1:當(dāng)折點(diǎn)附近區(qū)域內(nèi)d2>0 時(shí),越過折點(diǎn)后同步發(fā)電機(jī)轉(zhuǎn)速減小。若d1>0,發(fā)電機(jī)轉(zhuǎn)速先增大后減?。蝗鬱1<0,發(fā)電機(jī)轉(zhuǎn)速單調(diào)減小。

    情況2:當(dāng)折點(diǎn)附近區(qū)域內(nèi)d2<0 時(shí),軌跡越過折點(diǎn)后發(fā)電機(jī)轉(zhuǎn)速增大。若d1>0,發(fā)電機(jī)轉(zhuǎn)速單調(diào)增大;若d1<0,發(fā)電機(jī)轉(zhuǎn)速先減小后增大。

    基于2.3 節(jié)的分析,軌跡越過折點(diǎn)后鎖相環(huán)輸出角增加,故dθp/dt>0。折點(diǎn)附近區(qū)域內(nèi)dωg/dt的符號(hào)與?Pg/?θp的符號(hào)確定了d1和d2的負(fù)荷,決定了軌跡越過折點(diǎn)后發(fā)電機(jī)轉(zhuǎn)速的變化趨勢(shì)。

    忽略線路電導(dǎo)和PCC 處并聯(lián)負(fù)荷,得到?Pg/?θp的表達(dá)式(15)。由于函數(shù)的連續(xù)性,?Pg/?θp的符號(hào)在折點(diǎn)的小領(lǐng)域內(nèi)與折點(diǎn)處同號(hào),則可以用折點(diǎn)處?Pg/?θp的符號(hào)來近似表征其領(lǐng)域內(nèi)的性質(zhì)。而且基于式(15),?Pg/?θp=0 得到θp-δg-?I=2kπ,k取整數(shù)。那么在θp-δg平面上,折點(diǎn)右上側(cè)的區(qū)域與折點(diǎn)處的?Pg/?θp保持同號(hào)。因此,可以用折點(diǎn)處的靈敏度?Pg/?θp估計(jì)其右側(cè)軌跡的性質(zhì)。

    當(dāng)折點(diǎn)處?Pg/?θp>0時(shí),折點(diǎn)領(lǐng)域內(nèi)d2>0,那么對(duì)應(yīng)于上述情況1,軌跡越過折點(diǎn)后發(fā)電機(jī)轉(zhuǎn)速減小,有助于發(fā)電機(jī)恢復(fù)功角穩(wěn)定。當(dāng)折點(diǎn)處?Pg/?θp<0時(shí),折點(diǎn)領(lǐng)域內(nèi)d2<0,那么對(duì)應(yīng)于上述第2 種情況,軌跡越過折點(diǎn)后發(fā)電機(jī)轉(zhuǎn)速增加,功角更容易失穩(wěn)。

    兩種情況的示例分別如圖3(b)和(d)所示,其中比較了近似解式(14)與數(shù)值解,也驗(yàn)證了近似解和基于近似解分析的正確性。圖3(b)和(d)的轉(zhuǎn)速曲線分別對(duì)應(yīng)于圖3(a)和(c)中折點(diǎn)右側(cè)領(lǐng)域的綠色加粗實(shí)線軌跡。圖3(a)和(c)中藍(lán)色線表示軌跡,藍(lán)色點(diǎn)表示平衡點(diǎn),紅色點(diǎn)表示折點(diǎn),紫色虛線表示Vq=0 的解。本文后續(xù)圖中的藍(lán)色點(diǎn)和紅色點(diǎn)表示相同的含義。

    圖3 越過折點(diǎn)F 之后的軌跡Fig.3 Trajectory after crossing fold point F

    在圖3(b)的示例中,d2>0,d1<0,根據(jù)式(14),轉(zhuǎn)速為單調(diào)減小。轉(zhuǎn)速減小到小于1.0 p.u.之后功角也減小,故發(fā)電機(jī)功角恢復(fù)穩(wěn)定。在圖3(d)的示例中,折點(diǎn)右側(cè)領(lǐng)域d2<0,d1>0,根據(jù)式(14),轉(zhuǎn)速為單調(diào)增加。在折點(diǎn)左側(cè)領(lǐng)域,d2<0,d1<0,根據(jù)式(14),轉(zhuǎn)速為先減小再增加。因此,軌跡從折點(diǎn)左側(cè)到折點(diǎn)右側(cè)領(lǐng)域時(shí),轉(zhuǎn)速經(jīng)歷先減小再增加的過程。折點(diǎn)左側(cè)領(lǐng)域的轉(zhuǎn)速軌跡未在圖3 中畫出,可以從后文更加完整的轉(zhuǎn)速曲線中得到驗(yàn)證。

    3.2 鎖相環(huán)簡(jiǎn)化為代數(shù)方程的局限性

    根據(jù)3.1 節(jié)對(duì)折點(diǎn)右側(cè)領(lǐng)域內(nèi)同步發(fā)電機(jī)轉(zhuǎn)速的分析,變流器失去平衡點(diǎn)并非單一地惡化同步發(fā)電機(jī)功角穩(wěn)定。簡(jiǎn)單將鎖相環(huán)簡(jiǎn)化為代數(shù)方程,難以解釋如圖3(a)和(b)所示的情況1 的機(jī)理。

    以附錄D 表D1 中的網(wǎng)絡(luò)參數(shù)1 為例,通過對(duì)數(shù)值解物理意義的分析,說明情況1 下軌跡越過折點(diǎn)后較長(zhǎng)一段時(shí)間的物理過程。發(fā)電機(jī)的慣性時(shí)間常數(shù)為3.4 s,發(fā)電機(jī)等值內(nèi)阻抗為0.05 p.u.,等值內(nèi)電勢(shì)為1.1 p.u.,D=10 p.u.,Pm=3 p.u.。鎖相環(huán)參數(shù)kp=80 p.u./s,ki=1 800 p.u./s2,ε=0.003 68,I3=5 p.u.,?3=0。故障切除后,系統(tǒng)軌跡在θp-δg平面上的投影如圖3(a)中藍(lán)色線所示。發(fā)電機(jī)的電磁功率Pg先增大再減少,如圖4(a)所示,其中紅色虛線表示附錄C 圖C2(a)中折點(diǎn)F左邊部分Vq=0 的解Sa,斜線陰影部分表示θp增大過程中所“額外”增加的 減 速 面 積。此 時(shí),折 點(diǎn)F的 坐 標(biāo)(,)=(2.801 76,2.111 98),在折點(diǎn)處?Pg/?θp=1.944 44>0,其靈敏度變化曲線如圖4(b)所示。因此,軌跡越過折點(diǎn)F時(shí)發(fā)電機(jī)的電磁功率Pg也短時(shí)增加,形成了圖4(a)中的陰影區(qū)域。其物理意義可以解釋為鎖相環(huán)輸出角度的快速增加對(duì)發(fā)電機(jī)電磁功率的影響,改變了同步發(fā)電機(jī)功角的減速面積。

    圖4 動(dòng)態(tài)相互作用對(duì)發(fā)電機(jī)功角穩(wěn)定的影響:類型1Fig.4 Influence of dynamic interaction on power angle stability of generators:type 1

    當(dāng)故障更加嚴(yán)重時(shí),鎖相環(huán)輸出角度θp增大過程中所“額外”增加的減速面積也不足以使發(fā)電機(jī)功角減速,發(fā)電機(jī)才可能失去功角穩(wěn)定。在圖4(a)示例的基礎(chǔ)上,增加故障持續(xù)時(shí)間得到圖4(c)所示曲線,其中θp增大過程中對(duì)應(yīng)的發(fā)電機(jī)減速面積為藍(lán)色陰影部分,對(duì)應(yīng)的加速面積為紅色陰影部分。雖然軌跡越過折點(diǎn)F之后,減速面積增大,發(fā)電機(jī)轉(zhuǎn)速ωg迅速降低(對(duì)應(yīng)于圖4(d)中δg=2.8 rad),但ωg并未恢復(fù)到同步轉(zhuǎn)速1 p.u.。隨后電磁功率減小,加速面積進(jìn)一步增大,發(fā)電機(jī)很快失去同步。但是,此時(shí)失穩(wěn)的臨界狀態(tài)卻沒有對(duì)應(yīng)的“不穩(wěn)定平衡點(diǎn)”或者“主導(dǎo)不穩(wěn)定平衡點(diǎn)”。

    3.3 同步發(fā)電機(jī)功角穩(wěn)定受鎖相環(huán)制約的條件

    3.1 節(jié)分析中的情況2 的失穩(wěn)現(xiàn)象體現(xiàn)了鎖相環(huán)對(duì)同步發(fā)電機(jī)功角穩(wěn)定性的制約,因?yàn)殒i相環(huán)動(dòng)態(tài)過程惡化了發(fā)電機(jī)功角穩(wěn)定,鎖相環(huán)失穩(wěn)與同步發(fā)電機(jī)功角失穩(wěn)幾乎同時(shí)發(fā)生。因此,在折點(diǎn)F處?Pg/?θp<0 成為同步發(fā)電機(jī)功角穩(wěn)定受鎖相環(huán)制約的條件。

    將網(wǎng)絡(luò)參數(shù)Z2減小到1/6 p.u.,見附錄D 表D1中的網(wǎng)絡(luò)參數(shù)2。在折點(diǎn)處?Pg/?θp=-0.653 846 <0。系統(tǒng)軌跡在θp-δg平面上的投影如圖3(c)中藍(lán)色線所示,在其他平面上的投影如圖5 所示。圖3(c)中A點(diǎn)到B點(diǎn)之間的軌跡段對(duì)應(yīng)于圖5 中紅色方框里的部分。

    圖5 動(dòng)態(tài)相互作用對(duì)發(fā)電機(jī)功角穩(wěn)定的影響:類型2Fig.5 Influence of dynamic interaction on power angle stability of generators:type 2

    故障期間,發(fā)電機(jī)功角δg增大。故障切除后,軌跡經(jīng)過A點(diǎn)之后功角δg開始減速并回?cái)[,如圖5(a)和(b)所示。但在軌跡經(jīng)過B點(diǎn)之后,功角又開始加速增大,并且發(fā)電機(jī)功角失穩(wěn)。在從A點(diǎn)到B點(diǎn)這一過程中,注意到兩點(diǎn)特征:1)發(fā)電機(jī)功角軌跡在AB段之間有恢復(fù)穩(wěn)定的趨勢(shì),如圖5(b)和(c)所示,ωg<1.0 p.u.滿足發(fā)電機(jī)穩(wěn)定判據(jù);2)功角δg的變化很小,而鎖相環(huán)輸出角度θp變化范圍較大。圖5(c)表明,盡管功角在AB段中回?cái)[,但功率卻仍然往下跌落,與傳統(tǒng)交流系統(tǒng)中發(fā)電機(jī)的功角-功率曲線不同。圖5(d)說明,這個(gè)現(xiàn)象的原因是鎖相環(huán)輸出角度快速增大,且發(fā)電機(jī)的電磁功率對(duì)鎖相環(huán)輸出角度的靈敏度是負(fù)數(shù),導(dǎo)致發(fā)電機(jī)電磁功率快速減小,功角從減速變成加速。 因此,折點(diǎn)處?Pg/?θp<0 時(shí),越過折點(diǎn)F之后,鎖相環(huán)輸出角度增大的動(dòng)態(tài)過程成為功角失穩(wěn)的一個(gè)主要原因。

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

    在PSCAD/EMTDC 上搭建如附錄F 圖F1 所示的電路對(duì)上述分析進(jìn)行仿真驗(yàn)證。逆變側(cè)變流器外環(huán)采用定交流電壓控制和定直流電壓控制,整流側(cè)采用有功功率控制和無功功率控制。假設(shè)故障期間直流回路的卸荷耗能電路維持直流電壓的穩(wěn)定。在故障期間,逆變站PCC 處電壓跌落,觸發(fā)低電壓穿越控制模塊。此時(shí),根據(jù)電網(wǎng)導(dǎo)則要求,為了支撐交流電壓,有功電流指令值等于零,無功電流指令值等于逆變站允許的最大電流-imax。鎖相環(huán)PI環(huán)節(jié)的比例系數(shù)kp和積分增益ki分別為80 p.u./s 和1 800 p.u./s2,帶寬約等于16.16 Hz,阻尼比ξ約等于0.94[28,37]。發(fā)電機(jī)慣性時(shí)間常數(shù)換算到100 MW 的基準(zhǔn)值下為17 s,式(3)中的ε=0.000 74,其他參數(shù)見附錄D。

    4.1 算例1:鎖相環(huán)簡(jiǎn)化為代數(shù)方程的局限性驗(yàn)證

    穩(wěn)態(tài)時(shí),Pm=2.0 p.u.,I3=5.39 p.u.,?3=0.38 rad,網(wǎng)絡(luò)參數(shù)為附錄D 表D1 中網(wǎng)絡(luò)參數(shù)1。該條件屬于圖3 中的情況1。在0.01 s 時(shí)刻,發(fā)電機(jī)出口處母線發(fā)生三相接地短路瞬時(shí)故障。故障持續(xù)時(shí)間0.245 s 時(shí),折點(diǎn)附近的軌跡見附錄F 圖F2 陰影部分。在這段時(shí)間內(nèi),PCC 處電壓跌落嚴(yán)重,鎖相環(huán)因此失去平衡點(diǎn)。角速度ωp快速增大,迅速達(dá)到了上限值1.2 p.u.,導(dǎo)致鎖相環(huán)輸出角θp快速增大,而發(fā)電機(jī)的轉(zhuǎn)速逐漸減小。發(fā)電機(jī)功角δg變化較小,但發(fā)電機(jī)電磁功率快速增大,它與θp快速增大同時(shí)發(fā)生,驗(yàn)證了本文3.2 節(jié)的分析:?Pg/?θp>0 時(shí),θp增大初期提供的發(fā)電機(jī)P-δ平面上的“減速面積”有利于發(fā)電機(jī)保持功角穩(wěn)定。

    當(dāng)故障持續(xù)時(shí)間增加至0.37 s 時(shí),同步發(fā)電機(jī)才出現(xiàn)功角失穩(wěn),見附錄F 圖F3。θp快速增大導(dǎo)致的減速面積增加效應(yīng)仍然不足以使發(fā)電機(jī)功角回?cái)[。

    4.2 算例2:同步發(fā)電機(jī)功角穩(wěn)定受鎖相環(huán)制約條件的驗(yàn)證

    減小變流器與無窮大系統(tǒng)之間的等值線路阻抗至1/6 p.u.。穩(wěn) 態(tài) 時(shí),Pm=2.0 p.u.,I3=5.15 p.u.,?3=0.3 rad。該條件屬于圖3 中的第2 種情況。0.01 s 時(shí)刻發(fā)生三相短路瞬時(shí)故障,故障持續(xù)時(shí)間0.385 s,鎖相環(huán)和同步發(fā)電機(jī)都發(fā)生了失穩(wěn),折點(diǎn)附近軌跡見附錄F 圖F4 陰影部分。但是,系統(tǒng)軌跡越過折點(diǎn)之后表現(xiàn)出了與算例1 不同的特征。同步發(fā)電機(jī)電磁功率Pg在越過折點(diǎn)后繼續(xù)下降,然后才有小幅度的回升。鎖相環(huán)的快速動(dòng)態(tài)沒有立刻給發(fā)電機(jī)的轉(zhuǎn)子運(yùn)動(dòng)起到減速的作用。如果故障持續(xù)時(shí)間從0.385 s 減少到0.384 s,發(fā)電機(jī)和鎖相環(huán)都能保持穩(wěn)定,沒有算例1 中變流器失去平衡點(diǎn)而發(fā)電機(jī)功角恢復(fù)穩(wěn)定的狀態(tài)。

    折點(diǎn)附近?Pg/?θp<0,可以在接近折點(diǎn)F時(shí)給鎖相環(huán)輸出角θp施加-1 rad 的擾動(dòng)來觀察。0.58 s 在鎖相環(huán)輸出施加擾動(dòng)-1 rad,持續(xù)時(shí)間0.02 s。如圖6 中紅色線所示,施加-1 rad 的擾動(dòng)后,發(fā)電機(jī)電磁功率短時(shí)上升,說明?Pg/?θp<0。而且,擾動(dòng)之后發(fā)電機(jī)功角從折點(diǎn)位置開始回?cái)[,θp也恢復(fù)了穩(wěn)定,驗(yàn)證了本文3.3 節(jié)的分析:?Pg/?θp<0 時(shí),變流器失去平衡點(diǎn)對(duì)發(fā)電機(jī)功角穩(wěn)定有惡化的作用。

    圖6 鎖相環(huán)擾動(dòng)對(duì)發(fā)電機(jī)功角穩(wěn)定的影響Fig.6 Influence of PLL disturbance on generator power angle stability

    5 結(jié)語(yǔ)

    本文分析了鎖相環(huán)與同步發(fā)電機(jī)交互的條件,揭示了變流器失去平衡點(diǎn)對(duì)同步發(fā)電機(jī)暫態(tài)功角穩(wěn)定性的影響機(jī)理。本文的結(jié)論如下:

    1)在故障后系統(tǒng)恢復(fù)階段,發(fā)電機(jī)功角擺開過大可能使得變流器失去平衡點(diǎn),導(dǎo)致鎖相環(huán)與同步發(fā)電機(jī)的動(dòng)態(tài)過程耦合。尤其在含有高比例電力電子設(shè)備和較小容量同步發(fā)電機(jī)的系統(tǒng)中,這種耦合更容易主導(dǎo)功角穩(wěn)定域的邊界。因此,研究高比例新能源電力系統(tǒng)的暫態(tài)功角穩(wěn)定需要考慮鎖相環(huán)的影響。

    2)根據(jù)折點(diǎn)附近鎖相環(huán)輸出角度變化對(duì)發(fā)電機(jī)電磁功率的靈敏度符號(hào),鎖相環(huán)與發(fā)電機(jī)的交互作用可分為兩種類型:(1)鎖相環(huán)輸出角的快速變化為發(fā)電機(jī)提供了額外的減速面積,在一定程度上改善了發(fā)電機(jī)的暫態(tài)功角穩(wěn)定性;(2)鎖相環(huán)輸出角的快速變化惡化了發(fā)電機(jī)的暫態(tài)功角穩(wěn)定。鎖相環(huán)對(duì)同步發(fā)電機(jī)功角穩(wěn)定的制約體現(xiàn)在第2 種類型,因此,對(duì)于功角穩(wěn)定性分析需要著重考慮第2 種類型。該結(jié)論對(duì)研究含電壓源型變流器和同步發(fā)電機(jī)的電力系統(tǒng)的功角穩(wěn)定邊界條件、在線穩(wěn)定評(píng)估等具有參考意義。

    3)本文探討了鎖相環(huán)與同步發(fā)電機(jī)之間的非線性現(xiàn)象。由于在解析分析時(shí)截?cái)嗔巳A及以上的高階項(xiàng),然后結(jié)合數(shù)值計(jì)算現(xiàn)象對(duì)解析結(jié)果進(jìn)行了推廣,本文分析結(jié)論適用于鎖相環(huán)失穩(wěn)后的短暫時(shí)間。同時(shí),如果鎖相環(huán)失穩(wěn)時(shí)間較長(zhǎng)也會(huì)觸發(fā)繼電保護(hù)動(dòng)作導(dǎo)致變流器閉鎖。此外,本文的研究基于單個(gè)變流器聚合模型和單個(gè)發(fā)電機(jī)模型。對(duì)于更復(fù)雜的多變流器和多機(jī)系統(tǒng),本文的結(jié)論是否需要修正仍需要進(jìn)一步研究。

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

    猜你喜歡
    折點(diǎn)功角鎖相環(huán)
    虛擬調(diào)速器對(duì)VSG暫態(tài)功角穩(wěn)定影響機(jī)理分析
    含碘水折點(diǎn)加氯過程中消毒副產(chǎn)物的生成影響機(jī)制
    基于改進(jìn) shapelet 挖掘的風(fēng)電并網(wǎng)系統(tǒng)暫態(tài)功角穩(wěn)定評(píng)估
    能源工程(2019年6期)2019-12-02 01:58:20
    鎖相環(huán)HMC832的設(shè)計(jì)與實(shí)現(xiàn)
    新型無鎖相環(huán)DSTATCOM直接電流控制方法
    基于功角測(cè)量和等面積法則的發(fā)電機(jī)暫態(tài)穩(wěn)定在線判別研究
    某車型雙腔比例閥試驗(yàn)研究
    汽車科技(2017年4期)2017-08-08 04:20:50
    CLSI頭孢吡肟折點(diǎn)改變對(duì)腸桿菌科細(xì)菌的藥敏變化分析
    基于PSS/E風(fēng)電并網(wǎng)系統(tǒng)的靜態(tài)功角穩(wěn)定性分析
    新型鎖相環(huán)技術(shù)及仿真分析
    性色avwww在线观看| 一进一出抽搐gif免费好疼| 身体一侧抽搐| 男人舔女人的私密视频| 色在线成人网| 他把我摸到了高潮在线观看| 久久久水蜜桃国产精品网| a级毛片a级免费在线| or卡值多少钱| 久久亚洲真实| 亚洲黑人精品在线| 给我免费播放毛片高清在线观看| 三级毛片av免费| 两人在一起打扑克的视频| 欧美3d第一页| 久久久久国产一级毛片高清牌| 日韩欧美在线二视频| 最新美女视频免费是黄的| 99国产精品99久久久久| 亚洲av五月六月丁香网| 1000部很黄的大片| 欧美又色又爽又黄视频| 精品国产三级普通话版| 欧美黑人欧美精品刺激| 亚洲中文av在线| 国产高清有码在线观看视频| 国产亚洲精品久久久com| 非洲黑人性xxxx精品又粗又长| 校园春色视频在线观看| 亚洲av成人精品一区久久| 中文资源天堂在线| 亚洲成人免费电影在线观看| 久久人人精品亚洲av| 国产精品久久久人人做人人爽| 国产一区二区激情短视频| 日本a在线网址| 中亚洲国语对白在线视频| 狂野欧美白嫩少妇大欣赏| 少妇的丰满在线观看| 久久人妻av系列| 国产激情欧美一区二区| 欧美乱色亚洲激情| 国语自产精品视频在线第100页| 99久久久亚洲精品蜜臀av| 亚洲成人中文字幕在线播放| 波多野结衣高清作品| 天天添夜夜摸| 69av精品久久久久久| 香蕉国产在线看| 国产成人欧美在线观看| 亚洲人成网站高清观看| 免费看光身美女| 成年女人看的毛片在线观看| 日本与韩国留学比较| 免费电影在线观看免费观看| 综合色av麻豆| www.www免费av| 国产精品美女特级片免费视频播放器 | 久久久久久国产a免费观看| 亚洲精品中文字幕一二三四区| 麻豆国产av国片精品| 首页视频小说图片口味搜索| 两个人看的免费小视频| 国产亚洲精品av在线| 亚洲 欧美 日韩 在线 免费| 免费看a级黄色片| 午夜免费观看网址| 中文字幕熟女人妻在线| 久久精品91无色码中文字幕| 久99久视频精品免费| 欧美国产日韩亚洲一区| 午夜日韩欧美国产| 一级黄色大片毛片| 99国产极品粉嫩在线观看| 岛国视频午夜一区免费看| 欧美成人性av电影在线观看| 老熟妇乱子伦视频在线观看| 欧美日韩黄片免| 巨乳人妻的诱惑在线观看| 国产av不卡久久| 国产黄色小视频在线观看| 毛片女人毛片| 国产97色在线日韩免费| 亚洲在线自拍视频| 白带黄色成豆腐渣| 999久久久国产精品视频| avwww免费| 欧美日韩中文字幕国产精品一区二区三区| а√天堂www在线а√下载| 最近在线观看免费完整版| 欧美另类亚洲清纯唯美| 99热这里只有是精品50| 成人国产综合亚洲| 国产三级中文精品| 久久中文看片网| 制服丝袜大香蕉在线| 真人一进一出gif抽搐免费| 亚洲国产中文字幕在线视频| 国产成人精品久久二区二区免费| 久久香蕉精品热| 亚洲成人精品中文字幕电影| 五月玫瑰六月丁香| 亚洲国产欧美一区二区综合| 亚洲av电影不卡..在线观看| 日日摸夜夜添夜夜添小说| 国产乱人伦免费视频| 2021天堂中文幕一二区在线观| 免费电影在线观看免费观看| 亚洲av熟女| 亚洲人与动物交配视频| 亚洲第一电影网av| 欧美乱码精品一区二区三区| 亚洲中文字幕一区二区三区有码在线看 | 偷拍熟女少妇极品色| 给我免费播放毛片高清在线观看| 欧美另类亚洲清纯唯美| 国产伦在线观看视频一区| 亚洲成人久久性| 亚洲片人在线观看| 成在线人永久免费视频| 给我免费播放毛片高清在线观看| 国产精品久久电影中文字幕| 久久精品国产综合久久久| 久久精品91蜜桃| 少妇丰满av| 精品久久久久久久毛片微露脸| 成年人黄色毛片网站| 国产成人啪精品午夜网站| 韩国av一区二区三区四区| 亚洲专区中文字幕在线| 国产91精品成人一区二区三区| 日本免费a在线| 久9热在线精品视频| 国产aⅴ精品一区二区三区波| 激情在线观看视频在线高清| 国产av在哪里看| 日韩三级视频一区二区三区| 国产亚洲av嫩草精品影院| 国产久久久一区二区三区| 日韩国内少妇激情av| 日本免费一区二区三区高清不卡| 18禁国产床啪视频网站| 免费搜索国产男女视频| 美女高潮的动态| 欧美高清成人免费视频www| 人妻夜夜爽99麻豆av| 丰满人妻一区二区三区视频av | 九色成人免费人妻av| 日本一二三区视频观看| 狂野欧美激情性xxxx| 99国产精品99久久久久| 大型黄色视频在线免费观看| 国产精品一区二区三区四区久久| a级毛片在线看网站| 成熟少妇高潮喷水视频| 搡老熟女国产l中国老女人| 久久精品人妻少妇| 国产高清视频在线观看网站| www.www免费av| 成人18禁在线播放| 久久久久国产精品人妻aⅴ院| 亚洲精品一区av在线观看| 国产精品日韩av在线免费观看| 后天国语完整版免费观看| 国产aⅴ精品一区二区三区波| 亚洲专区国产一区二区| 精华霜和精华液先用哪个| 综合色av麻豆| 一进一出抽搐动态| 怎么达到女性高潮| 精品国产亚洲在线| 亚洲精品色激情综合| 在线观看午夜福利视频| 法律面前人人平等表现在哪些方面| 国产精品99久久久久久久久| 亚洲av成人不卡在线观看播放网| 国产精品日韩av在线免费观看| 午夜精品久久久久久毛片777| 国产 一区 欧美 日韩| 亚洲电影在线观看av| 色吧在线观看| 久久性视频一级片| 九色国产91popny在线| 亚洲人成网站在线播放欧美日韩| av国产免费在线观看| 亚洲九九香蕉| 无人区码免费观看不卡| 久久久久久久久免费视频了| 国产精品久久久久久久电影 | 国产午夜精品久久久久久| 欧美日韩一级在线毛片| 啦啦啦韩国在线观看视频| 夜夜看夜夜爽夜夜摸| 日本 av在线| 国产高清有码在线观看视频| 成人特级av手机在线观看| cao死你这个sao货| 欧美黄色淫秽网站| 欧美绝顶高潮抽搐喷水| 亚洲美女黄片视频| 免费看美女性在线毛片视频| 久久久国产成人精品二区| 天天一区二区日本电影三级| 18禁黄网站禁片免费观看直播| 精品国产乱码久久久久久男人| aaaaa片日本免费| 在线国产一区二区在线| 波多野结衣高清无吗| a级毛片a级免费在线| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品av在线| 女人高潮潮喷娇喘18禁视频| a级毛片a级免费在线| e午夜精品久久久久久久| 91九色精品人成在线观看| 久久久国产欧美日韩av| 欧美性猛交╳xxx乱大交人| 免费在线观看成人毛片| 国产亚洲精品综合一区在线观看| 男人舔奶头视频| 日韩欧美免费精品| 亚洲精品456在线播放app | 亚洲欧美日韩无卡精品| 我要搜黄色片| 特大巨黑吊av在线直播| 岛国在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 性色av乱码一区二区三区2| 久久草成人影院| 在线十欧美十亚洲十日本专区| 日韩 欧美 亚洲 中文字幕| 久久亚洲真实| 一本精品99久久精品77| 成人特级黄色片久久久久久久| 丰满人妻一区二区三区视频av | cao死你这个sao货| 成人18禁在线播放| 全区人妻精品视频| 99精品欧美一区二区三区四区| av女优亚洲男人天堂 | 欧美一区二区国产精品久久精品| 99久久国产精品久久久| 久久精品国产综合久久久| 国产精品久久久av美女十八| 成人午夜高清在线视频| www.精华液| 日本成人三级电影网站| 欧美黄色片欧美黄色片| 麻豆国产av国片精品| 一级毛片高清免费大全| 99国产极品粉嫩在线观看| 亚洲欧美一区二区三区黑人| 久久久久国产一级毛片高清牌| 国产黄片美女视频| 国产精品,欧美在线| 国产亚洲av嫩草精品影院| 久久久久久久久免费视频了| 日韩欧美 国产精品| or卡值多少钱| 99久国产av精品| 麻豆成人av在线观看| 亚洲熟妇中文字幕五十中出| 国产精品美女特级片免费视频播放器 | 12—13女人毛片做爰片一| 国产欧美日韩精品亚洲av| 久99久视频精品免费| 在线国产一区二区在线| 亚洲成人中文字幕在线播放| 一个人观看的视频www高清免费观看 | 成人鲁丝片一二三区免费| 欧美不卡视频在线免费观看| 一进一出好大好爽视频| 青草久久国产| 欧美最黄视频在线播放免费| 日韩欧美国产一区二区入口| 精品国产乱码久久久久久男人| 国产成人啪精品午夜网站| 亚洲一区二区三区色噜噜| 欧美日韩中文字幕国产精品一区二区三区| 久久久久性生活片| 亚洲欧美日韩东京热| 午夜a级毛片| 免费看十八禁软件| 脱女人内裤的视频| 成人国产一区最新在线观看| 国产精品女同一区二区软件 | 国产精品电影一区二区三区| 亚洲美女黄片视频| 午夜福利在线观看吧| 高潮久久久久久久久久久不卡| 国产av在哪里看| 日本 欧美在线| 国产综合懂色| 欧美一区二区精品小视频在线| 极品教师在线免费播放| 色av中文字幕| 精品国产超薄肉色丝袜足j| 久久久精品欧美日韩精品| 狂野欧美激情性xxxx| 黑人巨大精品欧美一区二区mp4| 亚洲精品粉嫩美女一区| 久久精品国产99精品国产亚洲性色| 国语自产精品视频在线第100页| 亚洲熟妇中文字幕五十中出| 日本撒尿小便嘘嘘汇集6| 国产乱人伦免费视频| 一二三四社区在线视频社区8| 在线看三级毛片| 国产视频一区二区在线看| 国产精品九九99| 特大巨黑吊av在线直播| 国产精品98久久久久久宅男小说| 亚洲精品色激情综合| xxxwww97欧美| 99视频精品全部免费 在线 | 亚洲中文日韩欧美视频| 国内揄拍国产精品人妻在线| 麻豆国产av国片精品| 成人高潮视频无遮挡免费网站| 国产精品久久久久久精品电影| 18禁裸乳无遮挡免费网站照片| 亚洲中文日韩欧美视频| 久久亚洲精品不卡| 精品欧美国产一区二区三| 一本综合久久免费| 日韩欧美三级三区| av片东京热男人的天堂| 老司机在亚洲福利影院| 美女黄网站色视频| 久久久久性生活片| 三级男女做爰猛烈吃奶摸视频| 哪里可以看免费的av片| 国产69精品久久久久777片 | a在线观看视频网站| 亚洲中文日韩欧美视频| 久久精品夜夜夜夜夜久久蜜豆| 熟女人妻精品中文字幕| 成人国产综合亚洲| 亚洲精品粉嫩美女一区| 欧美日韩亚洲国产一区二区在线观看| 日本黄大片高清| 男女之事视频高清在线观看| 九色国产91popny在线| 久久人人精品亚洲av| 国产精华一区二区三区| 色综合婷婷激情| 在线永久观看黄色视频| 老司机午夜十八禁免费视频| 国产精品1区2区在线观看.| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 午夜视频精品福利| 国产成人精品无人区| 99久久精品一区二区三区| 美女午夜性视频免费| 国产欧美日韩一区二区精品| 少妇熟女aⅴ在线视频| 国产成年人精品一区二区| 久久精品国产99精品国产亚洲性色| 欧美日韩黄片免| 国产成人精品无人区| 一级毛片女人18水好多| 国产真人三级小视频在线观看| 亚洲色图 男人天堂 中文字幕| 欧美一级a爱片免费观看看| 国内精品美女久久久久久| 91av网站免费观看| 国产成人欧美在线观看| 变态另类成人亚洲欧美熟女| 久久久久久久午夜电影| 久久久久精品国产欧美久久久| 亚洲真实伦在线观看| 免费看十八禁软件| 一本综合久久免费| 国产精品一区二区三区四区久久| 精品久久蜜臀av无| 久久久久久人人人人人| 欧美zozozo另类| 国产伦精品一区二区三区四那| 国产高清激情床上av| 黄色片一级片一级黄色片| 国产精品久久久久久久电影 | 欧美乱妇无乱码| 岛国在线观看网站| 亚洲无线在线观看| 久久午夜综合久久蜜桃| av天堂在线播放| 1024手机看黄色片| 国语自产精品视频在线第100页| 中文字幕熟女人妻在线| 亚洲av成人精品一区久久| 少妇丰满av| 欧美色欧美亚洲另类二区| 动漫黄色视频在线观看| 制服丝袜大香蕉在线| 国内精品久久久久精免费| 国产高清激情床上av| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 欧美在线一区亚洲| 国产成人精品久久二区二区免费| 99热只有精品国产| 成人亚洲精品av一区二区| 麻豆久久精品国产亚洲av| 日本 欧美在线| 国产精品99久久久久久久久| 亚洲精品美女久久久久99蜜臀| 校园春色视频在线观看| 日本五十路高清| a在线观看视频网站| 亚洲国产精品合色在线| or卡值多少钱| 看免费av毛片| 精品国产美女av久久久久小说| 国产成人欧美在线观看| 好男人电影高清在线观看| 免费观看的影片在线观看| 久久久久国产精品人妻aⅴ院| 精品久久久久久久毛片微露脸| 国产在线精品亚洲第一网站| www日本黄色视频网| xxxwww97欧美| 欧美一区二区精品小视频在线| 美女cb高潮喷水在线观看 | 色播亚洲综合网| 国产免费男女视频| 哪里可以看免费的av片| 最新美女视频免费是黄的| 国产主播在线观看一区二区| 亚洲精品456在线播放app | 久久中文字幕一级| 日日夜夜操网爽| 五月玫瑰六月丁香| 在线播放国产精品三级| 一级a爱片免费观看的视频| 亚洲人与动物交配视频| 国产99白浆流出| 欧美成人一区二区免费高清观看 | 18禁美女被吸乳视频| 人人妻人人看人人澡| 长腿黑丝高跟| 国产成人av激情在线播放| 91麻豆精品激情在线观看国产| 精品日产1卡2卡| 村上凉子中文字幕在线| 女警被强在线播放| 999精品在线视频| 性色avwww在线观看| 免费看光身美女| 制服丝袜大香蕉在线| 国产综合懂色| 亚洲精品在线观看二区| 操出白浆在线播放| 黄色 视频免费看| 母亲3免费完整高清在线观看| 国产午夜精品论理片| 久久久久久大精品| 国内精品久久久久精免费| 美女黄网站色视频| 黄色丝袜av网址大全| 国产高清激情床上av| 午夜福利成人在线免费观看| 亚洲美女黄片视频| 国产成人av教育| 亚洲成人免费电影在线观看| 久久精品国产清高在天天线| 成人18禁在线播放| 国产一区二区在线观看日韩 | 视频区欧美日本亚洲| 国产欧美日韩精品一区二区| 免费av毛片视频| 91麻豆av在线| 巨乳人妻的诱惑在线观看| 日本黄色片子视频| 成人午夜高清在线视频| 亚洲无线观看免费| 欧美精品啪啪一区二区三区| 一夜夜www| 亚洲一区二区三区不卡视频| 91麻豆精品激情在线观看国产| 又大又爽又粗| 国产毛片a区久久久久| 欧美黑人欧美精品刺激| 最近最新中文字幕大全电影3| 成人性生交大片免费视频hd| 亚洲成人久久性| 两性夫妻黄色片| 亚洲中文日韩欧美视频| 国产精品九九99| 国产精品亚洲av一区麻豆| 精品久久久久久,| 免费看日本二区| 在线视频色国产色| 国产亚洲av嫩草精品影院| 手机成人av网站| 欧美在线一区亚洲| 成人精品一区二区免费| 欧美黄色片欧美黄色片| 国产精品一区二区免费欧美| 午夜久久久久精精品| 免费电影在线观看免费观看| 久久国产精品人妻蜜桃| 亚洲av美国av| www.精华液| 男女做爰动态图高潮gif福利片| 色视频www国产| 999久久久国产精品视频| bbb黄色大片| 成人av一区二区三区在线看| 国产精品香港三级国产av潘金莲| 国产精品久久电影中文字幕| 成人av一区二区三区在线看| 99久久国产精品久久久| 日韩欧美国产一区二区入口| 国产精品日韩av在线免费观看| 99国产精品99久久久久| av在线天堂中文字幕| 国产亚洲欧美98| 最好的美女福利视频网| 亚洲成a人片在线一区二区| 在线a可以看的网站| 最近最新免费中文字幕在线| 免费av不卡在线播放| 男插女下体视频免费在线播放| 国产精品一区二区精品视频观看| 天堂动漫精品| 怎么达到女性高潮| 琪琪午夜伦伦电影理论片6080| 级片在线观看| 日韩免费av在线播放| 在线观看免费午夜福利视频| 全区人妻精品视频| 亚洲色图 男人天堂 中文字幕| 成人精品一区二区免费| 露出奶头的视频| 韩国av一区二区三区四区| 久久久久久大精品| 国产激情久久老熟女| 国产成人系列免费观看| www.精华液| 美女大奶头视频| 精品久久久久久成人av| 国产欧美日韩一区二区三| 99久久成人亚洲精品观看| 久久久久久九九精品二区国产| 热99re8久久精品国产| 人人妻人人澡欧美一区二区| 少妇的丰满在线观看| 少妇丰满av| 免费观看的影片在线观看| 人妻夜夜爽99麻豆av| 国产成人一区二区三区免费视频网站| 国产精品久久久久久亚洲av鲁大| 88av欧美| 麻豆一二三区av精品| 天天添夜夜摸| 中文资源天堂在线| 亚洲七黄色美女视频| 国产成人精品久久二区二区免费| 一个人看视频在线观看www免费 | 男人舔女人下体高潮全视频| 亚洲天堂国产精品一区在线| 哪里可以看免费的av片| 亚洲成av人片在线播放无| 亚洲精品美女久久久久99蜜臀| 真人做人爱边吃奶动态| 最近视频中文字幕2019在线8| 国产人伦9x9x在线观看| 中文字幕人妻丝袜一区二区| 免费观看精品视频网站| 成人特级av手机在线观看| 一区二区三区国产精品乱码| 听说在线观看完整版免费高清| 青草久久国产| 日韩精品青青久久久久久| 99久国产av精品| 午夜免费观看网址| 国内精品久久久久精免费| 露出奶头的视频| 日韩精品中文字幕看吧| 男人的好看免费观看在线视频| 中文字幕av在线有码专区| 亚洲av第一区精品v没综合| 麻豆成人av在线观看| 久久久色成人| 亚洲午夜精品一区,二区,三区| 日韩av在线大香蕉| 小说图片视频综合网站| 国产成人精品久久二区二区91| netflix在线观看网站| 啪啪无遮挡十八禁网站| 亚洲精品在线美女| 黄色片一级片一级黄色片| 亚洲一区二区三区色噜噜| 精品一区二区三区四区五区乱码| 国产成人系列免费观看| 中文字幕最新亚洲高清| 一二三四社区在线视频社区8| 9191精品国产免费久久| 亚洲aⅴ乱码一区二区在线播放| 成年女人看的毛片在线观看| 久久亚洲真实| 悠悠久久av| 久久久久久久午夜电影| 色老头精品视频在线观看| 国产精华一区二区三区| 99国产极品粉嫩在线观看| av片东京热男人的天堂| 久久久成人免费电影| 99在线视频只有这里精品首页| 久久人人精品亚洲av| 天堂影院成人在线观看| 亚洲av日韩精品久久久久久密| 欧美乱码精品一区二区三区| 亚洲精华国产精华精|