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

    考慮電網(wǎng)特性的動(dòng)態(tài)無(wú)功優(yōu)化配置方法

    2016-05-23 07:02:08黃小慶阮馳騁鄒佳芯曹一家關(guān)維德
    電力自動(dòng)化設(shè)備 2016年9期
    關(guān)鍵詞:暫態(tài)排序容量

    黃小慶,阮馳騁,鄒佳芯,曹一家,關(guān)維德

    (1.湖南大學(xué) 電氣與信息工程學(xué)院,湖南 長(zhǎng)沙 410082;2.國(guó)網(wǎng)青海省電力公司電力科學(xué)研究院,青海 西寧 810008)

    0 引言

    電能需求的不斷增加使電力系統(tǒng)的電壓穩(wěn)定問(wèn)題備受關(guān)注。靜止同步補(bǔ)償器STATCOM(STATic synchronous COMpensator)和靜止無(wú)功補(bǔ)償器SVC(Static Var Compensator)等動(dòng)態(tài)無(wú)功補(bǔ)償裝置能在系統(tǒng)發(fā)生故障后快速提供動(dòng)態(tài)電壓支撐,提高系統(tǒng)的暫態(tài)電壓穩(wěn)定水平。

    為充分發(fā)揮動(dòng)態(tài)無(wú)功補(bǔ)償裝置的作用,必須選擇合適的安裝點(diǎn)和容量。目前,對(duì)安裝點(diǎn)的研究多采用靜態(tài)電壓穩(wěn)定的分析方法,如以提高功率裕度和電壓穩(wěn)定裕度[1-2]為目標(biāo)的方法、基于 U-Q 曲線[3]的方法、先導(dǎo)節(jié)點(diǎn)法[4]和主導(dǎo)節(jié)點(diǎn)法[5]等,或基于向量場(chǎng)正規(guī)形理論,根據(jù)非線性參與因子[6]確定動(dòng)態(tài)無(wú)功補(bǔ)償裝置的安裝點(diǎn),以提高系統(tǒng)的靜態(tài)電壓穩(wěn)定性。上述方法選擇安裝點(diǎn),既沒(méi)有考慮感應(yīng)電動(dòng)機(jī)等快速響應(yīng)設(shè)備的動(dòng)態(tài)特性,也沒(méi)有考慮全局補(bǔ)償效果。文獻(xiàn)[7-8]基于軌跡靈敏度指標(biāo)TSI(Trajectory Sensitivity Index),文獻(xiàn)[9]以最大限度延長(zhǎng)電網(wǎng)關(guān)鍵故障下的極限切除時(shí)間為目標(biāo),進(jìn)行動(dòng)態(tài)無(wú)功優(yōu)化配置。上述方法完全依靠時(shí)域仿真計(jì)算,仿真場(chǎng)景的選擇對(duì)其結(jié)果有一定的影響,而且配置方案只考慮了安裝效果,并未考慮經(jīng)濟(jì)性問(wèn)題。

    相應(yīng)地,動(dòng)態(tài)無(wú)功優(yōu)化配置方法得到研究,如基于阻抗模裕度指標(biāo)IMMI(Impedance Modulus Margin Index)[10]、采用靜態(tài)負(fù)荷裕度分析法和動(dòng)態(tài)仿真分析[11]等方法,綜合考慮了多種故障模式下的無(wú)功補(bǔ)償裝置的整體規(guī)劃。特別地,兼顧配置效果和經(jīng)濟(jì)性的無(wú)功配置方法得到重視,如文獻(xiàn)[12]提出了基于暫態(tài)電壓嚴(yán)重性指標(biāo)TVSI(Transient Voltage Severity Index)和多目標(biāo)分解進(jìn)化算法 MOEA/D(Multi-Objective Evolutionary Algorithm based on Decomposition)的無(wú)功優(yōu)化。文獻(xiàn)[13]提出了一種可同時(shí)考慮靜態(tài)與動(dòng)態(tài)無(wú)功補(bǔ)償?shù)臒o(wú)功補(bǔ)償優(yōu)化模型,在此基礎(chǔ)上通過(guò)動(dòng)態(tài)仿真和靜態(tài)分析相結(jié)合的方法迭代求解最優(yōu)方法。

    但是,目前的動(dòng)態(tài)無(wú)功優(yōu)化配置方法大多僅從某一特定的角度進(jìn)行分析。分析角度不同,造成配置方案差別較大??紤]到各種方法之間并不對(duì)立,如果能利用不同方法的優(yōu)勢(shì),動(dòng)態(tài)地對(duì)無(wú)功優(yōu)化配置問(wèn)題進(jìn)行全面的分析,將減少計(jì)算量并提高配置裝置效率。

    基于此,本文提出了一種考慮電網(wǎng)特性的多準(zhǔn)則動(dòng)態(tài)無(wú)功優(yōu)化配置方法。準(zhǔn)則包括:節(jié)點(diǎn)電壓穩(wěn)定的重要度(下文簡(jiǎn)稱重要度)、節(jié)點(diǎn)電壓穩(wěn)定的薄弱度(下文簡(jiǎn)稱薄弱度)以及改進(jìn)TSI。其中,重要度指標(biāo)可識(shí)別電網(wǎng)拓?fù)浣Y(jié)構(gòu)和功率特性對(duì)節(jié)點(diǎn)電壓穩(wěn)定的影響,明確需重點(diǎn)進(jìn)行電壓監(jiān)測(cè)和保護(hù)的節(jié)點(diǎn);薄弱度指標(biāo)可識(shí)別自身暫態(tài)電壓較不穩(wěn)定的節(jié)點(diǎn);改進(jìn)TSI可描述無(wú)功注入量對(duì)全系統(tǒng)暫態(tài)電壓水平的實(shí)際提升效果。

    各準(zhǔn)則存在一定的內(nèi)在聯(lián)系,但各有側(cè)重。本文利用不同準(zhǔn)則的優(yōu)勢(shì),多角度、全面地篩選對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定有重要影響且自身易出現(xiàn)暫態(tài)電壓失穩(wěn)的節(jié)點(diǎn)作為候選安裝點(diǎn);選擇能兼顧全系統(tǒng)暫態(tài)電壓穩(wěn)定的節(jié)點(diǎn)作為最終安裝點(diǎn);并結(jié)合暫態(tài)電壓穩(wěn)定標(biāo)準(zhǔn)選取最優(yōu)配置容量。

    1 動(dòng)態(tài)無(wú)功優(yōu)化配置的3個(gè)準(zhǔn)則

    1.1 節(jié)點(diǎn)電壓穩(wěn)定的重要度指標(biāo)

    重要度指標(biāo)是指考慮電網(wǎng)拓?fù)浣Y(jié)構(gòu)影響時(shí),維持各節(jié)點(diǎn)電壓穩(wěn)定的重要程度。電力系統(tǒng)是由大量的節(jié)點(diǎn)和節(jié)點(diǎn)之間的連線(邊)組成的復(fù)雜網(wǎng)絡(luò)。節(jié)點(diǎn)作為網(wǎng)絡(luò)的核心元素,其連接方式影響著網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)、信息能量流通的效率和成本。而電力系統(tǒng)規(guī)模龐大,無(wú)法監(jiān)視和調(diào)整每一個(gè)節(jié)點(diǎn)的電壓,通常只能關(guān)注一些重要節(jié)點(diǎn)。

    由于電壓等級(jí)較高的節(jié)點(diǎn)在電網(wǎng)中起著輸送和分配大容量功率的作用,這些節(jié)點(diǎn)一旦發(fā)生電壓失穩(wěn)故障,對(duì)電網(wǎng)的沖擊將很大。以中樞節(jié)點(diǎn)為例,一旦其因電壓失穩(wěn)退出運(yùn)行后,可能造成網(wǎng)絡(luò)結(jié)構(gòu)的大破壞,進(jìn)而引起潮流轉(zhuǎn)移或連鎖性事故擾動(dòng),最終導(dǎo)致系統(tǒng)電壓崩潰。需重點(diǎn)監(jiān)測(cè)和保護(hù)這些節(jié)點(diǎn)電壓。

    本文利用重要度指標(biāo)對(duì)維持各節(jié)點(diǎn)的電壓穩(wěn)定的重要程度進(jìn)行量化評(píng)估。考慮到網(wǎng)絡(luò)結(jié)構(gòu)和節(jié)點(diǎn)注入功率均對(duì)節(jié)點(diǎn)電壓穩(wěn)定重要度指標(biāo)存在影響,本文將節(jié)點(diǎn)的重要度定義為:

    其中,Ii為節(jié)點(diǎn)i的電壓穩(wěn)定重要度;Di為節(jié)點(diǎn)i的度,表示直接與節(jié)點(diǎn)i相連的邊的數(shù)量;η1、η2為權(quán)重系數(shù),且為節(jié)點(diǎn) i的功率特性,計(jì)算公式見(jiàn)式(2)。

    其中,Si為節(jié)點(diǎn)i的注入功率;Sbase為選定的系統(tǒng)基準(zhǔn)功率。反映了節(jié)點(diǎn)i在系統(tǒng)中傳輸和分配功率的重要程度,其值越大,說(shuō)明此節(jié)點(diǎn)在系統(tǒng)中傳輸和分配的功率越多,一旦出現(xiàn)故障,影響就越大。

    式(1)中,η1Di反映節(jié)點(diǎn)的拓?fù)浣Y(jié)構(gòu)特性,它不受潮流分布的影響,其值越大,說(shuō)明該節(jié)點(diǎn)退出運(yùn)行后越有可能引起更多節(jié)點(diǎn)或線路的故障,系統(tǒng)維持拓?fù)浣Y(jié)構(gòu)完整和正常運(yùn)行的能力越低,在系統(tǒng)中越重要隨著電網(wǎng)潮流分布的變化而變化,反映了節(jié)點(diǎn)的功率特性,其值越大,說(shuō)明該節(jié)點(diǎn)越重要。

    根據(jù)節(jié)點(diǎn)的拓?fù)浣Y(jié)構(gòu)特性與功率特性對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定性的相對(duì)影響,選取合適的η1和η2的值,由式(1)即可得到各節(jié)點(diǎn)的電壓穩(wěn)定的重要度Ii。

    1.2 基于偏差量的電壓薄弱度指標(biāo)

    靈敏度法具有物理概念簡(jiǎn)單明確、判據(jù)嚴(yán)格準(zhǔn)確的優(yōu)點(diǎn)[14]。傳統(tǒng)的靈敏度法判定系統(tǒng)電壓穩(wěn)定的判據(jù)如下[15]:

    (1)dUL/dQL<0,當(dāng)負(fù)荷節(jié)點(diǎn)無(wú)功需求 QL減少(或增加)時(shí),該節(jié)點(diǎn)電壓UL隨之上升(或下降),此時(shí)系統(tǒng)電壓是穩(wěn)定的;

    (2)dQL/dQg>0,當(dāng)負(fù)荷節(jié)點(diǎn)的無(wú)功需求 QL增加(或減少)引起發(fā)電機(jī)無(wú)功輸出Qg的增加(或減少)時(shí),系統(tǒng)電壓是穩(wěn)定的。

    當(dāng)時(shí)間足夠短時(shí),判據(jù)(1)和(2)可寫為差分形式 ΔUL/ΔQL和 ΔQL/ΔQg來(lái)判斷電壓穩(wěn)定。 負(fù)荷節(jié)點(diǎn)的有功和無(wú)功隨負(fù)荷特性變化,但在短時(shí)間內(nèi)變化不明顯。因此,當(dāng)負(fù)荷節(jié)點(diǎn)的負(fù)荷水平變化較小時(shí),即使系統(tǒng)遭受大擾動(dòng)或運(yùn)行方式發(fā)生突變,負(fù)荷節(jié)點(diǎn)的無(wú)功的差分量變化也較小。相對(duì)應(yīng)地,負(fù)荷節(jié)點(diǎn)的電壓幅值卻隨系統(tǒng)擾動(dòng)變化明顯,且引起系統(tǒng)電壓失穩(wěn)的主要原因是系統(tǒng)無(wú)法維持無(wú)功功率的動(dòng)態(tài)平衡和缺乏合適的電壓支持[16]。因此,利用發(fā)電機(jī)節(jié)點(diǎn)的無(wú)功功率以及負(fù)荷節(jié)點(diǎn)的電壓幅值判斷系統(tǒng)的電壓穩(wěn)定性更為合適。

    定義基于偏差量的節(jié)點(diǎn)電壓的薄弱度指標(biāo)為:

    其中,ΩL為負(fù)荷節(jié)點(diǎn)的集合;ΩG為發(fā)電機(jī)節(jié)點(diǎn)的集合;ΔUi為負(fù)荷節(jié)點(diǎn)i的電壓幅值變化量;ΔQj為發(fā)電機(jī)節(jié)點(diǎn)j的無(wú)功功率變化量。

    因此,僅需相鄰時(shí)刻電壓和無(wú)功功率的變化量即可計(jì)算式(3)定義的節(jié)點(diǎn)電壓薄弱度指標(biāo)。由判據(jù)(1)和(2)可知,當(dāng)系統(tǒng)電壓穩(wěn)定時(shí),負(fù)荷節(jié)點(diǎn)無(wú)功需求減少(或增加),其電壓隨之上升(或下降),同時(shí)引起發(fā)電機(jī)無(wú)功輸出的減少(或增加)。即電壓穩(wěn)定的判據(jù)為:

    Uwi越靠近0,表明無(wú)功功率對(duì)節(jié)點(diǎn)電壓變化調(diào)節(jié)能力越弱,即節(jié)點(diǎn)的電壓越薄弱。據(jù)此,依據(jù)各節(jié)點(diǎn)的電壓薄弱度排序,可以辨識(shí)電壓薄弱節(jié)點(diǎn)。為減少計(jì)算誤差并避免僅憑一次判斷可能出現(xiàn)的錯(cuò)判,以一段時(shí)間內(nèi)的指標(biāo)均值判斷節(jié)點(diǎn)電壓的薄弱度。將式(3)改寫為:

    其中為單位時(shí)間內(nèi)節(jié)點(diǎn)電壓薄弱度指標(biāo)的均值;Uwi(λ,Δti)為 Δti內(nèi)的指標(biāo)值,Δti為采樣周期,λ 為Δti內(nèi)與該指標(biāo)計(jì)算相關(guān)的參數(shù),包括各發(fā)電機(jī)節(jié)點(diǎn)的無(wú)功功率和負(fù)荷節(jié)點(diǎn)的電壓;T為采樣總時(shí)間;m=T/Δti,為時(shí)間T內(nèi)的采樣周期數(shù)。

    1.3 改進(jìn)TSI

    軌跡靈敏度可從側(cè)面反映節(jié)點(diǎn)電壓隨無(wú)功功率實(shí)時(shí)變化的程度,適用于暫態(tài)穩(wěn)定的分析與控制中。本文采用TSI確定對(duì)全系統(tǒng)暫態(tài)電壓水平提升最明顯的安裝點(diǎn)。多機(jī)電力系統(tǒng)的動(dòng)態(tài)行為通??捎萌缦挛⒎执鷶?shù)方程表示:

    其中,f為發(fā)電機(jī)及其勵(lì)磁系統(tǒng)、負(fù)荷和動(dòng)態(tài)無(wú)功補(bǔ)償裝置等的動(dòng)態(tài)特性;g為系統(tǒng)的潮流模型;x和y分別為系統(tǒng)的狀態(tài)變量(如發(fā)電機(jī)功角和角速度)和代數(shù)變量(如節(jié)點(diǎn)的電壓和相角);λ為系統(tǒng)參數(shù),本文表示節(jié)點(diǎn)注入的無(wú)功功率。式(6)對(duì)應(yīng)的初值為:

    若將式(6)的解記為 x(λ,t)和 y(λ,t),對(duì)其按泰勒級(jí)數(shù)展開(kāi)并忽略高階項(xiàng)可得:

    其中,xλ(t)和 yλ(t)分別為變量 x、y 對(duì)參數(shù) λ 的軌跡靈敏度。當(dāng)Δλ足夠小時(shí),可用式(9)近似計(jì)算軌跡靈敏度 xλ(t)和 yλ(t)。

    文獻(xiàn)[7]將 TSI定義為:

    其中,i、j(i,j=1,2,…,n)為節(jié)點(diǎn)編號(hào),n 為系統(tǒng)的節(jié)點(diǎn)總數(shù);tk(k=1,2,…,Nk)為第 k 個(gè)采樣時(shí)刻,Nk為采樣時(shí)間點(diǎn)總數(shù);Uj為節(jié)點(diǎn)j的電壓;Qi為節(jié)點(diǎn)i注入系統(tǒng)的無(wú)功功率;(?Uj/?Qi)t=tk為 t=tk時(shí) Uj關(guān)于Qi的軌跡靈敏度值。

    文獻(xiàn)[7]采用式(9)的方法近似計(jì)算軌跡靈敏度,即從故障發(fā)生時(shí)刻向節(jié)點(diǎn)i直接注入恒定的無(wú)功功率Q,取有、無(wú)注入Q時(shí)電壓Uj在t=tk時(shí)的差值除以 Q 得到(?Uj/?Qi)t=tk。 文獻(xiàn)[8]在考慮動(dòng)態(tài)無(wú)功補(bǔ)償裝置實(shí)際特性的基礎(chǔ)上,取有、無(wú)在節(jié)點(diǎn)i安裝額定容量為Q的動(dòng)態(tài)無(wú)功補(bǔ)償裝置下,電壓Uj在t=tk時(shí)的差值除以 Q 得到(?Uj/?Qi)t=tk,結(jié)果更準(zhǔn)確,但計(jì)算時(shí)無(wú)功功率的增量仍為恒定值,與真實(shí)的軌跡跟蹤過(guò)程存在誤差,無(wú)法準(zhǔn)確體現(xiàn)實(shí)際電壓支撐效果。

    事實(shí)上,動(dòng)態(tài)無(wú)功補(bǔ)償裝置注入系統(tǒng)的無(wú)功功率隨安裝點(diǎn)電壓的變化而變化,并不能維持恒定。因此本文采用動(dòng)態(tài)無(wú)功補(bǔ)償裝置的實(shí)時(shí)、實(shí)際無(wú)功注入計(jì)算TSI,該無(wú)功功率可由Q-t曲線得出。安裝前不存在補(bǔ)償裝置注入的無(wú)功功率,即Q-t曲線面積為0。安裝后,Q-t曲線的面積即可反映補(bǔ)償裝置實(shí)際注入系統(tǒng)的無(wú)功功率。根據(jù)定積分的定義,將TSI的求解轉(zhuǎn)化為對(duì)各節(jié)點(diǎn)U-t曲線和安裝點(diǎn)Q-t曲線在相同時(shí)間內(nèi)多個(gè)采樣時(shí)間點(diǎn)數(shù)據(jù)的處理。本文基于文獻(xiàn)[7]和文獻(xiàn)[8]對(duì) TSI改進(jìn)如下:

    其中,Uj(tk,Qi0)為故障后 t=tk(k=1,2,…,Nk)時(shí)節(jié)點(diǎn)j的電壓,此時(shí)節(jié)點(diǎn)i注入系統(tǒng)的無(wú)功功率為Qi0;ΔQi(tk)為節(jié)點(diǎn)i安裝動(dòng)態(tài)無(wú)功補(bǔ)償裝置后在t=tk時(shí)注入系統(tǒng)的無(wú)功功率增量,隨時(shí)間不斷變化。

    本文所做改進(jìn)使TSI的計(jì)算方法更符合其物理意義,理論上可更準(zhǔn)確地反映各候選點(diǎn)安裝動(dòng)態(tài)無(wú)功補(bǔ)償裝置對(duì)全系統(tǒng)暫態(tài)電壓水平的提升效果,計(jì)算結(jié)果更準(zhǔn)確。

    2 多準(zhǔn)則動(dòng)態(tài)無(wú)功優(yōu)化配置方法

    由于電網(wǎng)的本質(zhì)特性對(duì)其性能有深刻的影響,電網(wǎng)的拓?fù)浣Y(jié)構(gòu)和功率特性決定了系統(tǒng)中存在關(guān)系電壓穩(wěn)定性的重要節(jié)點(diǎn)。這些節(jié)點(diǎn)因電壓穩(wěn)定問(wèn)題退出運(yùn)行后將造成系統(tǒng)結(jié)構(gòu)的大破壞,進(jìn)而引起潮流轉(zhuǎn)移或連鎖性事故,使系統(tǒng)維持正常運(yùn)行的能力大幅降低,造成更大規(guī)模的電壓穩(wěn)定問(wèn)題。盡管維持這些節(jié)點(diǎn)的電壓穩(wěn)定較其他節(jié)點(diǎn)更為重要,但其實(shí)際運(yùn)行電壓可能處于相對(duì)穩(wěn)定的狀態(tài)。另一方面,電壓薄弱度指標(biāo)只能反映該節(jié)點(diǎn)的運(yùn)行狀態(tài),不能反映在該節(jié)點(diǎn)安裝動(dòng)態(tài)無(wú)功償裝置對(duì)全系統(tǒng)暫態(tài)電壓穩(wěn)定水平的提升效果。因此,可采用多準(zhǔn)則的方法兼顧不同影響因素,選擇最優(yōu)安裝點(diǎn)。

    極限切除時(shí)間表示切除故障的極限允許時(shí)間,可反映系統(tǒng)的暫態(tài)電壓穩(wěn)定水平。不同動(dòng)態(tài)無(wú)功配置方案對(duì)于關(guān)鍵故障的極限切除時(shí)間會(huì)有不同程度的提高[17]。因此,本文通過(guò)比較不同配置方式對(duì)多個(gè)關(guān)鍵故障極限切除時(shí)間的提升效果,確定最佳配置位置,并以滿足暫態(tài)電壓穩(wěn)定標(biāo)準(zhǔn)的最小容量作為兼顧技術(shù)性能與經(jīng)濟(jì)性的最優(yōu)容量,進(jìn)而確定最優(yōu)配置方案。具體步驟如下。

    a.由式(1)計(jì)算各節(jié)點(diǎn)電壓穩(wěn)定重要度指標(biāo)并排序,指標(biāo)值較大者對(duì)維持系統(tǒng)電壓穩(wěn)定更為重要。

    b.通過(guò)時(shí)域掃描確定威脅系統(tǒng)暫態(tài)電壓穩(wěn)定的關(guān)鍵故障,并對(duì)關(guān)鍵故障進(jìn)行時(shí)域仿真。由式(3)和式(5)計(jì)算各節(jié)點(diǎn)在不同時(shí)段T內(nèi)的電壓薄弱度指標(biāo)并排序,統(tǒng)計(jì)各節(jié)點(diǎn)電壓薄弱程度在各時(shí)段排名靠前的次數(shù),次數(shù)較多的節(jié)點(diǎn)電壓較薄弱。

    c.篩選出在電壓穩(wěn)定重要度指標(biāo)排序靠前的節(jié)點(diǎn)和在電壓薄弱度指標(biāo)排序中靠前的節(jié)點(diǎn)作為候選安裝點(diǎn)集。

    d.依次在各候選安裝點(diǎn)安裝動(dòng)態(tài)無(wú)功補(bǔ)償裝置,再次對(duì)關(guān)鍵故障進(jìn)行時(shí)域仿真。由式(11)計(jì)算各候選安裝點(diǎn)的改進(jìn)TSI并排序,指標(biāo)值最大者為最優(yōu)安裝點(diǎn)。

    e.考慮實(shí)際安裝容量的限制,選取一定的動(dòng)態(tài)無(wú)功容量,按平均和步驟d的排序等方式(改進(jìn)TSI越大,分配的容量越多),集中或分散地給各安裝點(diǎn)分配動(dòng)態(tài)無(wú)功補(bǔ)償容量。比較不同配置方式對(duì)多個(gè)關(guān)鍵故障后極限切除時(shí)間的延長(zhǎng)值,確定對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定性提升效果最好的配置方式。

    f.針對(duì)步驟e得出的動(dòng)態(tài)無(wú)功補(bǔ)償?shù)奈恢煤腿萘浚孟到y(tǒng)暫態(tài)電壓穩(wěn)定標(biāo)準(zhǔn)對(duì)安裝效果進(jìn)行校核。選取滿足要求所需的最小容量作為最優(yōu)容量,從而確定最優(yōu)配置方案。

    3 算例分析

    本文以IEEE 39節(jié)點(diǎn)系統(tǒng)和STATCOM為例,說(shuō)明上述動(dòng)態(tài)無(wú)功優(yōu)化配置方法實(shí)施的具體步驟。

    3.1 節(jié)點(diǎn)電壓穩(wěn)定的重要度指標(biāo)的計(jì)算與分析

    IEEE 39節(jié)點(diǎn)系統(tǒng)的拓?fù)浣Y(jié)構(gòu)如圖1所示。取η1=η2=0.5,Sbase=100 MV·A,節(jié)點(diǎn)的視在功率參考新英格蘭10機(jī)39節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)數(shù)據(jù),計(jì)算IEEE 39節(jié)點(diǎn)系統(tǒng)各節(jié)點(diǎn)的電壓穩(wěn)定的重要度。重要度由大到小的順序排序,如表1所示。

    圖1 IEEE 39節(jié)點(diǎn)系統(tǒng)拓?fù)浣Y(jié)構(gòu)圖Fig.1 Topology of IEEE 39-bus system

    表1 節(jié)點(diǎn)電壓穩(wěn)定重要度排序Table 1 Nodes reordered according to voltage stability importance index

    由表1可知:

    (1)節(jié)點(diǎn)38、39的電壓穩(wěn)定重要度都很大,但由于它們是PV節(jié)點(diǎn),電壓幅值恒定,故不考慮其作為候選安裝節(jié)點(diǎn);

    (2)節(jié)點(diǎn)8的電壓穩(wěn)定重要度在PQ節(jié)點(diǎn)中是最大的,說(shuō)明它的電壓一旦失穩(wěn),可能引起連鎖性事故擾動(dòng),甚至導(dǎo)致整個(gè)系統(tǒng)崩潰,需要重點(diǎn)關(guān)注;

    (3)節(jié)點(diǎn)度數(shù)和電壓穩(wěn)定重要度的排序并不完全一致,因?yàn)殡妷悍€(wěn)定重要度指標(biāo)不僅與電網(wǎng)的拓?fù)浣Y(jié)構(gòu)有關(guān),還與其功率特性有關(guān)。

    3.2 節(jié)點(diǎn)電壓穩(wěn)定的薄弱度指標(biāo)的計(jì)算與分析

    對(duì)系統(tǒng)進(jìn)行時(shí)域掃描后發(fā)現(xiàn)線路16-17的三相永久性故障是可能引起全系統(tǒng)暫態(tài)電壓失穩(wěn)的關(guān)鍵故障。因此本算例的時(shí)域仿真均在t=1s時(shí)線路16-17發(fā)生三相永久性故障,故障持續(xù)時(shí)間為0.1 s,采樣時(shí)間間隔為0.01 s,各節(jié)點(diǎn)負(fù)荷水平保持不變的條件下進(jìn)行。

    采用30%感應(yīng)電動(dòng)機(jī)負(fù)荷+70%恒阻抗負(fù)荷進(jìn)行時(shí)域仿真。 取式(5)中的 T=0.1s,在 1.0~2.5s 內(nèi)選取負(fù)荷節(jié)點(diǎn)的電壓和發(fā)電機(jī)節(jié)點(diǎn)的無(wú)功功率各150個(gè)采樣時(shí)間點(diǎn)的值,由式(3)和式(5)計(jì)算各節(jié)點(diǎn)的電壓薄弱度指標(biāo)。統(tǒng)計(jì)各負(fù)荷節(jié)點(diǎn)在各時(shí)段指標(biāo)值排序前10名的次數(shù),并按電壓薄弱度由高到低的順序排序,如表2所示。

    表2 故障下各節(jié)點(diǎn)電壓薄弱度排序Table 2 Nodes reordered according to voltage stability weakness index for specified contingency

    由表2可知:

    (1)節(jié)點(diǎn)13、14、15的電壓薄弱度指標(biāo)排名前10的次數(shù)較多,說(shuō)明當(dāng)系統(tǒng)出現(xiàn)擾動(dòng)時(shí),這些節(jié)點(diǎn)相對(duì)其他節(jié)點(diǎn)而言更容易出現(xiàn)暫態(tài)電壓穩(wěn)定問(wèn)題;

    (2)綜合表1和表2的排序結(jié)果可知,表1中部分電壓穩(wěn)定重要度未達(dá)到平均值的PQ節(jié)點(diǎn)(節(jié)點(diǎn) 1、2、9、10、17、18、22),在表 2 中排名也靠后,反映了電網(wǎng)本質(zhì)特性與其性能的聯(lián)系;

    (3)把表1排序前10的PQ節(jié)點(diǎn)和表2排序前10的節(jié)點(diǎn)作為候選安裝點(diǎn)集,計(jì)算其改進(jìn)TSI。

    3.3 改進(jìn)TSI的計(jì)算與分析

    采用30%感應(yīng)電動(dòng)機(jī)負(fù)荷+70%恒阻抗負(fù)荷進(jìn)行時(shí)域仿真,計(jì)算候選安裝點(diǎn)的改進(jìn)TSI。首先驗(yàn)證本文對(duì)TSI所做改進(jìn)的有效性。文獻(xiàn)[8]已驗(yàn)證采用其改進(jìn)TSI確定安裝點(diǎn)比文獻(xiàn)[7]中傳統(tǒng)的TSI更準(zhǔn)確。因此本文在STATCOM容量同為100 Mvar和仿真條件相同的情況下,針對(duì)IEEE 39節(jié)點(diǎn)算例對(duì)文獻(xiàn)[8]和本文的TSI排序結(jié)果進(jìn)行對(duì)比。本文改進(jìn)TSI排名前 10 的節(jié)點(diǎn)依次為 19、20、16、22、23、21、15、14、13、11;文獻(xiàn)[8]的 TSI排名前 10 的節(jié)點(diǎn)依次為 19、16、22、21、23、15、20、14、13、10。

    文獻(xiàn)[8]與本文的排序結(jié)果相比,2種方法得出的最優(yōu)安裝點(diǎn)均是節(jié)點(diǎn)19,說(shuō)明節(jié)點(diǎn)19確為對(duì)系統(tǒng)暫態(tài)電壓水平提升效果最好的安裝點(diǎn)。兩者排序結(jié)果并不完全相同。當(dāng)系統(tǒng)負(fù)荷為50%感應(yīng)電動(dòng)機(jī)負(fù)荷+50%恒阻抗負(fù)荷時(shí),在2種方法的TSI排序靠前的節(jié)點(diǎn)安裝100 Mvar的STATCOM,考察在不同節(jié)點(diǎn)安裝動(dòng)態(tài)無(wú)功補(bǔ)償裝置后,電壓最薄弱的節(jié)點(diǎn)15在故障后的電壓恢復(fù)情況如圖2所示(U為標(biāo)幺值)。

    圖2 節(jié)點(diǎn)15故障下的電壓-時(shí)間曲線Fig.2 Voltage-time curve for specified contingency of node 15

    由圖2可知,分別在節(jié)點(diǎn)19和在文獻(xiàn)[8]中排名第7的節(jié)點(diǎn)20安裝STATCOM均能使節(jié)點(diǎn)15的電壓最終恢復(fù)穩(wěn)定,安裝效果優(yōu)于節(jié)點(diǎn)16、21、23。圖2所示結(jié)果與本文的排序完全相符,證明本文的改進(jìn)TSI能更真實(shí)地反映實(shí)際軌跡跟蹤,計(jì)算結(jié)果更準(zhǔn)確。

    為確定最優(yōu)安裝點(diǎn),利用各節(jié)點(diǎn)的U-t曲線和STATCOM注入系統(tǒng)的實(shí)際無(wú)功功率的Q-t曲線在1.1~2.6s內(nèi) 150 個(gè)采樣時(shí)間點(diǎn)的值,由式(11)計(jì)算候選安裝點(diǎn)的改進(jìn)TSI,按由高到低的順序排序于表3。

    表3 改進(jìn)TSI排序Table 3 Nodes reordered according to improved TSI

    為研究動(dòng)態(tài)無(wú)功補(bǔ)償裝置的容量對(duì)改進(jìn)TSI排序結(jié)果的影響,本文分別在各候選安裝點(diǎn)安裝容量為 50 Mvar、100 Mvar和 150 Mvar的 STATCOM,計(jì)算各節(jié)點(diǎn)的改進(jìn)TSI,結(jié)果如圖3所示。由圖3可知,不同安裝容量下各節(jié)點(diǎn)的排序總體上是一致的,這是由于單個(gè)動(dòng)態(tài)無(wú)功補(bǔ)償裝置的容量相對(duì)于實(shí)際系統(tǒng)的裝機(jī)容量是很小的,安裝容量的選擇一般不會(huì)影響各節(jié)點(diǎn)改進(jìn)TSI的排序結(jié)果。

    圖3 不同STATCOM容量下各節(jié)點(diǎn)的改進(jìn)TSIFig.3 Improved TSI of different nodes for different STATCOM capacities

    3.4 動(dòng)態(tài)無(wú)功最優(yōu)配置方案的確定

    設(shè)單臺(tái)STATCOM的最大容量為100 Mvar。為進(jìn)一步研究不同配置方式對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定性的提升效果,由表3的排序結(jié)果制定了6種不同的配置方案,如表4所示。其中,第7種配置方案為后文中求出的最優(yōu)配置方案。各方案負(fù)荷均為50%感應(yīng)電動(dòng)機(jī)負(fù)荷+50%恒阻抗負(fù)荷,各節(jié)點(diǎn)所配置的容量比例與表3中節(jié)點(diǎn)的排序?qū)?yīng)。各方案在線路發(fā)生三相永久性故障后極限切除時(shí)間的延長(zhǎng)值如表5所示。

    表4 7種配置方案Table 4 Seven configuration schemes

    表5 不同方案的動(dòng)態(tài)無(wú)功配置效果Table 5 Effect of dynamic var configuration for different schemes

    由表5可知,僅在單個(gè)節(jié)點(diǎn)安裝STATCOM對(duì)系統(tǒng)的暫態(tài)電壓穩(wěn)定性的提升效果有限,應(yīng)分散安裝在多個(gè)節(jié)點(diǎn)。增加STATCOM的安裝容量有助于提升系統(tǒng)的暫態(tài)電壓穩(wěn)定性。比較方案2、5可知,在總安裝容量相同的情況下,選取少數(shù)幾個(gè)改進(jìn)TSI排名靠前的候選點(diǎn)安裝STATCOM的效果優(yōu)于在更多節(jié)點(diǎn)分散安裝,且設(shè)備的投資和管理費(fèi)用較少。比較方案3、4或方案5、6可知,根據(jù)改進(jìn)TSI的排序,按一定比例分配各點(diǎn)的安裝容量比平均分配對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定性的提升更明顯。綜上,在節(jié)點(diǎn)19、16、15、14按其改進(jìn)TSI排序分配STATCOM的安裝容量的配置方式效果最佳。按照三相永久性故障0.1 s切除,故障清除后1 s時(shí)負(fù)荷母線電壓恢復(fù)到0.75 p.u.以上的暫態(tài)電壓穩(wěn)定標(biāo)準(zhǔn)[15]校核,得出最優(yōu)配置方案7:節(jié)點(diǎn)19為48Mvar,節(jié)點(diǎn)16為36Mvar,節(jié)點(diǎn) 15為 24 Mvar,節(jié)點(diǎn) 14為 12 Mvar。

    3.5 配置方案經(jīng)濟(jì)性分析

    STATCOM較其他無(wú)功補(bǔ)償裝置有著優(yōu)越的性能,但同時(shí)其高昂的成本也是限制其發(fā)展的重要因素。因此,關(guān)于配置方案的投資經(jīng)濟(jì)性考慮也就成為了一個(gè)重要的方面。選取7種配置方案進(jìn)行經(jīng)濟(jì)性對(duì)比,對(duì)比結(jié)果如表6所示。其中,總的經(jīng)濟(jì)成本分為安裝成本和設(shè)備購(gòu)買成本[18],設(shè)備購(gòu)買成本約為 32萬(wàn)元 /Mvar[19],安裝成本每個(gè)節(jié)點(diǎn) 30萬(wàn)元。

    表6 不同方案的經(jīng)濟(jì)性比較Table 6 Comparison of costs among different schemes

    由表6可知,7種方案中成本投入最少的為方案1,其總成本為3230萬(wàn)元,但是由表5可知,方案1的總延長(zhǎng)值僅為0.07 s,其暫態(tài)電壓穩(wěn)定性的提升效果有限;而方案7的總成本為3960萬(wàn)元,略高于方案1,遠(yuǎn)低于其他幾種配置方案,但它對(duì)暫態(tài)電壓穩(wěn)定性的提升約為0.12 s,比方案1要高得多,故方案7比方案1更優(yōu)。方案2—6的成本均遠(yuǎn)大于方案7,但是對(duì)暫態(tài)電壓穩(wěn)定性的提升效果并不理想,其資金投入和提升效果不成比例,故相比之下方案7較優(yōu)。

    4 結(jié)語(yǔ)

    本文從維持整個(gè)系統(tǒng)電壓穩(wěn)定的角度,從對(duì)維持系統(tǒng)電壓穩(wěn)定重要且自身電壓穩(wěn)定相對(duì)薄弱的節(jié)點(diǎn)中,選擇動(dòng)態(tài)無(wú)功優(yōu)化裝置的配置候選點(diǎn)。利用改進(jìn)的TSI排序,描述候選節(jié)點(diǎn)對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定性提升效果,進(jìn)一步篩選出安裝節(jié)點(diǎn)。結(jié)合考慮電力系統(tǒng)暫態(tài)電壓穩(wěn)定標(biāo)準(zhǔn)和經(jīng)濟(jì)性約束,得出補(bǔ)償裝置的最佳配置位置和具體配置容量。

    本文的多準(zhǔn)則配置方法,既考慮了系統(tǒng)拓?fù)浣Y(jié)構(gòu)和功率特性對(duì)無(wú)功補(bǔ)償效果的影響,又考慮了節(jié)點(diǎn)電壓本身的薄弱程度,兼顧了單位無(wú)功補(bǔ)償對(duì)系統(tǒng)暫態(tài)電壓穩(wěn)定的調(diào)節(jié)效率。該方法同時(shí)考慮了系統(tǒng)的靜態(tài)物理特性和節(jié)點(diǎn)、系統(tǒng)的動(dòng)態(tài)運(yùn)行狀態(tài),顧及注入動(dòng)態(tài)無(wú)功對(duì)全系統(tǒng)暫態(tài)電壓穩(wěn)定的影響,以及電力系統(tǒng)相關(guān)標(biāo)準(zhǔn)和經(jīng)濟(jì)性的約束,為動(dòng)態(tài)無(wú)功配置裝置的選點(diǎn)和定容,提供了有別于“基于單一準(zhǔn)則選點(diǎn)”的新思路。

    參考文獻(xiàn):

    [1]劉傳銓,張焰.電力系統(tǒng)無(wú)功補(bǔ)償點(diǎn)及其補(bǔ)償容量的確定[J].電網(wǎng)技術(shù),2007,31(12):78-81.LIU Chuanquan,ZHANG Yan.Confirmation of reactive power compensation node and its optimal compensation capacity [J].Power System Technology,2007,31(12):78-81.

    [2]POURBEIK P,KOESSLER R,QUAINTANCE N,et al.Performing comprehensive voltage stability studies for the determination of optimal location,size and type of reactive compensation[C]∥Proceedings of IEEE Power Systems Conference and Exposition.Atlanta,GA,USA:IEEE,2006:118.

    [3]CIGERE TF 38-01-03.Planning against voltage collapse [R].Paris,F(xiàn)rance:CIGRE,1987.

    [4]LAKKARAJU T,F(xiàn)ELIACHI A.Selection of pilot buses for VAR support considering N-1 contingency criteria[C]∥Proceedings of IEEE Power Systems Conference and Exposition.Atlanta,GA,USA:IEEE,2006:1513-1517.

    [5]劉君華,江秀臣,方鴿飛.一種確定無(wú)功源最佳配置地點(diǎn)與數(shù)量的新方法[J].電網(wǎng)技術(shù),2007,31(16):72-76.LIU Junhua,JIANG Xiuchen,F(xiàn)ANG Gefei.A new method to decide optimal allocation sites and number of reactive sources[J].Power System Technology,2007,31(16):72-76.

    [6]張靖,程時(shí)杰,文勁宇,等.通過(guò)選擇SVC安裝地點(diǎn)提高靜態(tài)電壓穩(wěn)定性的新方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2007,27(34):7-11.ZHANG Jing,CHENG Shijie,WEN Jinyu,et al.A novel steadystate voltage stability enhancementmethod based on SVC allocation[J].Proceedings of the CSEE,2007,27(34):7-11.

    [7]SAPKOTA B,VITTAL V.Dynamicvarplanningin a large power system using trajectory sensitivity[J].IEEE Transactions on Power Systems,2010,25(1):461-469.

    [8]黃弘揚(yáng),楊汾艷,徐政,等.基于改進(jìn)軌跡靈敏度指標(biāo)的動(dòng)態(tài)無(wú)功優(yōu)化配置方法[J].電網(wǎng)技術(shù),2012,36(2):88-94.HUANG Hongyang,YANG Fenyan,XU Zheng,et al.A dynamic var configuration method based on improved trajectory sensitivity index[J].Power System Technology,2012,36(2):88-94.

    [9]袁志昌,劉文華,宋強(qiáng).基于暫態(tài)電壓穩(wěn)定指標(biāo)的動(dòng)態(tài)無(wú)功優(yōu)化配置方法[J].電力系統(tǒng)自動(dòng)化,2009,33(14):17-21.YUAN Zhichang,LIU Wenhua,SONG Qiang.Optimal allocation method of dynamic var compensation based on transient voltage stability index[J].Automation of Electric Power Systems,2009,33(14):17-21.

    [10]李帥虎,曹一家,劉光曄.基于阻抗模裕度指標(biāo)的動(dòng)態(tài)無(wú)功補(bǔ)償裝置優(yōu)化配置方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,34(22):3791-3798.LI Shuaihu,CAO Yijia,LIU Guangye.Optimal allocation method of dynamic var compensator based on the impedance modulus margin index[J].Proceedings of the CSEE,2014,34(22):3791-3798.

    [11]ZHANG Yongjun,CHEN Chao,LI Yong,et al.Dynamic voltage support planning for receiving end power systems based on evaluation of state separating and transferring risks[J].Electric Power Systems Research,2010,80(12):1520-1527.

    [12]XU Yan,ZHAO Yangdong,MENG Ke,etal.Multi-objective dynamic varplanning againstshort-term voltage instability using a decomposition-based evolutionary algorithm [J].IEEE Transactions on Power Systems,2014,29(6):2813-2822.

    [13]張勇軍,韓東昆,劉巍.受端電網(wǎng)靜態(tài)/動(dòng)態(tài)無(wú)功補(bǔ)償綜合優(yōu)化[J].電力自動(dòng)化設(shè)備,2009,29(3):32-35.ZHANG Yongjun,HAN Dongkun,LIU Wei.Comprehensive optimization of static/dynamic reactive power compensation in the receiving end[J].Electric Power Automation Equipment,2009,29(3):32-35.

    [14]蘇永春,程時(shí)杰,文勁宇,等.電力系統(tǒng)電壓穩(wěn)定性及其研究現(xiàn)狀(一)[J].電力自動(dòng)化設(shè)備,2006,26(6):97-101.SU Yongchun,CHENG Shijie,WEN Jinyu,etal.Research status of voltage stability in power system[J].Electric Power Automation Equipment,2006,26(6):97-101.

    [15]袁駿,段獻(xiàn)忠,何仰贊,等.電力系統(tǒng)電壓穩(wěn)定靈敏度分析方法綜述[J].電網(wǎng)技術(shù),1997,21(9):7-10.YUAN Jun,DUAN Xianzhong,HE Yangzan,et al.Summarization of the sensitivity analysis method of voltage stability in power system[J].Power System Technology,1997,21(9):7-10.

    [16]段俊東,郭志忠.一種可在線確定電壓穩(wěn)定運(yùn)行范圍的方法[J].中國(guó)電機(jī)工程學(xué)報(bào),2006,26(4):113-120.DUAN Jundong,GUO Zhizhong.A new methodforon-line determination of the capacity curves of voltage stability [J].Proceedings of the CSEE,2006,26(4):113-120.

    [17]李海琛,劉明波,林舜江.考慮暫態(tài)電壓安全的STATCOM安裝地點(diǎn)選擇和容量?jī)?yōu)化[J].電力系統(tǒng)保護(hù)與控制,2011,39(5):69-76.LI Haichen,LIU Mingbo,LIN Shunjiang.Selection of STATCOM’s installation site and capacity optimization considering transient voltage security[J].Power System Protection and Control,2011,39(5):69-76.

    [18]張勇軍,任震,鐘紅梅,等.基于災(zāi)變遺傳算法的無(wú)功規(guī)劃優(yōu)化[J].電力系統(tǒng)自動(dòng)化,2002,26(23):29-32.ZHANG Yongjun,REN Zhen,ZHONG Hongmei,et al.Optimal reactive power planning based on Catastrophic genetic algorithm[J].Automation of Electric Power Systems,2002,26(23):29-32.

    [19]TIWARI A,AJJARAPU V.Optimal allocation of dynamic VAR supportusing mixed integerdynamicoptimization [J].IEEE Transactions on Power Systems,2011,26(1):305-314.

    猜你喜歡
    暫態(tài)排序容量
    排序不等式
    300Mvar空冷隱極同步調(diào)相機(jī)暫態(tài)特性仿真分析
    恐怖排序
    節(jié)日排序
    電力系統(tǒng)全網(wǎng)一體化暫態(tài)仿真接口技術(shù)
    電子制作(2018年14期)2018-08-21 01:38:28
    刻舟求劍
    兒童繪本(2018年5期)2018-04-12 16:45:32
    除氧器暫態(tài)計(jì)算研究
    SnO2納米片容量異常行為的新解釋
    2015年上半年我國(guó)風(fēng)電新增并網(wǎng)容量916萬(wàn)千瓦
    風(fēng)能(2015年8期)2015-02-27 10:15:12
    2015年一季度我國(guó)風(fēng)電新增并網(wǎng)容量470萬(wàn)千瓦
    風(fēng)能(2015年5期)2015-02-27 10:14:46
    在线观看午夜福利视频| 青春草视频在线免费观看| 久久午夜亚洲精品久久| 国产精品三级大全| 日日啪夜夜撸| 久久久久久久久久成人| 晚上一个人看的免费电影| 国产黄色视频一区二区在线观看 | 悠悠久久av| kizo精华| 老司机福利观看| 亚洲精品影视一区二区三区av| 国产私拍福利视频在线观看| 一进一出抽搐gif免费好疼| 天天一区二区日本电影三级| 国产精品野战在线观看| 色视频www国产| 久久久久久国产a免费观看| 黄片wwwwww| 国产高潮美女av| 我的老师免费观看完整版| 欧美日本视频| 国内精品久久久久精免费| 欧美色欧美亚洲另类二区| 亚洲精品国产av成人精品| 欧美日韩乱码在线| 色综合亚洲欧美另类图片| 变态另类成人亚洲欧美熟女| 99久久精品一区二区三区| 美女高潮的动态| 国产淫片久久久久久久久| 岛国在线免费视频观看| 有码 亚洲区| av国产免费在线观看| 女同久久另类99精品国产91| 日韩成人伦理影院| 成人性生交大片免费视频hd| 色播亚洲综合网| 欧美精品国产亚洲| 国产av在哪里看| 亚洲精华国产精华液的使用体验 | 日本黄大片高清| 国内精品一区二区在线观看| 麻豆一二三区av精品| 欧美成人免费av一区二区三区| 给我免费播放毛片高清在线观看| 日日摸夜夜添夜夜爱| 亚洲精品粉嫩美女一区| 成人国产麻豆网| 欧美精品一区二区大全| 三级经典国产精品| 人人妻人人看人人澡| 亚洲欧美日韩高清在线视频| 一卡2卡三卡四卡精品乱码亚洲| 国产亚洲精品久久久com| 一个人免费在线观看电影| 欧美精品国产亚洲| 美女国产视频在线观看| 久久人人精品亚洲av| 人妻夜夜爽99麻豆av| 亚洲性久久影院| 九草在线视频观看| 日韩欧美国产在线观看| 日韩一区二区三区影片| 九九爱精品视频在线观看| 国产精品日韩av在线免费观看| av免费在线看不卡| 男人舔奶头视频| 国产精品乱码一区二三区的特点| 国产精品永久免费网站| 午夜精品在线福利| 禁无遮挡网站| 亚洲婷婷狠狠爱综合网| 免费看美女性在线毛片视频| 有码 亚洲区| 精品午夜福利在线看| 男女下面进入的视频免费午夜| av在线观看视频网站免费| 成年女人看的毛片在线观看| 黄色一级大片看看| 尾随美女入室| 国产精品电影一区二区三区| 免费观看精品视频网站| 一本精品99久久精品77| 亚洲av不卡在线观看| 美女 人体艺术 gogo| 久久久色成人| 免费观看的影片在线观看| 国产精品日韩av在线免费观看| 国产精品av视频在线免费观看| 久久久午夜欧美精品| 变态另类成人亚洲欧美熟女| 成人亚洲欧美一区二区av| 亚洲欧美成人精品一区二区| 国产激情偷乱视频一区二区| 天堂av国产一区二区熟女人妻| 日韩欧美国产在线观看| 日韩 亚洲 欧美在线| 国产片特级美女逼逼视频| 日韩欧美精品免费久久| 亚洲国产精品成人综合色| 亚洲最大成人av| 国产精品人妻久久久影院| 一边亲一边摸免费视频| av免费观看日本| 欧美一区二区精品小视频在线| 在线天堂最新版资源| 美女脱内裤让男人舔精品视频 | 亚洲最大成人中文| 国产伦精品一区二区三区四那| 亚洲美女搞黄在线观看| 秋霞在线观看毛片| 国产麻豆成人av免费视频| 国产老妇伦熟女老妇高清| 欧美成人一区二区免费高清观看| 国产精品女同一区二区软件| 亚洲欧洲日产国产| 欧美精品国产亚洲| 亚洲国产欧洲综合997久久,| 久久久久久久久大av| 国产精品日韩av在线免费观看| 日韩欧美一区二区三区在线观看| 国产真实伦视频高清在线观看| 色哟哟哟哟哟哟| 日本成人三级电影网站| 三级男女做爰猛烈吃奶摸视频| 两性午夜刺激爽爽歪歪视频在线观看| 97超碰精品成人国产| 啦啦啦观看免费观看视频高清| 国产高清三级在线| 美女高潮的动态| 国产视频内射| 精品无人区乱码1区二区| 一本久久中文字幕| 国内精品宾馆在线| 国产麻豆成人av免费视频| 日韩亚洲欧美综合| 国产黄片视频在线免费观看| av视频在线观看入口| 我要看日韩黄色一级片| 99久久无色码亚洲精品果冻| 亚洲av一区综合| 久久精品综合一区二区三区| 亚洲av.av天堂| 久久婷婷人人爽人人干人人爱| 欧美bdsm另类| 26uuu在线亚洲综合色| 欧美日本亚洲视频在线播放| av在线蜜桃| 全区人妻精品视频| 人妻夜夜爽99麻豆av| 天天躁夜夜躁狠狠久久av| 欧美三级亚洲精品| 只有这里有精品99| 不卡一级毛片| 男的添女的下面高潮视频| 2022亚洲国产成人精品| 国产乱人视频| 国产精品无大码| 99久久精品一区二区三区| 别揉我奶头 嗯啊视频| 日韩大尺度精品在线看网址| 又爽又黄a免费视频| 成熟少妇高潮喷水视频| 一本久久精品| 国内精品宾馆在线| 久久久久网色| 久久精品夜色国产| 国产大屁股一区二区在线视频| 国产大屁股一区二区在线视频| 国产色爽女视频免费观看| 国产在视频线在精品| 亚洲在线自拍视频| 国产黄片美女视频| 神马国产精品三级电影在线观看| 搡女人真爽免费视频火全软件| 成人特级av手机在线观看| 亚洲最大成人手机在线| 老司机福利观看| 国产成人午夜福利电影在线观看| 特级一级黄色大片| 欧美又色又爽又黄视频| 精品久久久久久久久av| 亚洲欧美日韩卡通动漫| 亚洲欧美日韩高清专用| 性欧美人与动物交配| 国产大屁股一区二区在线视频| 成人欧美大片| 老熟妇乱子伦视频在线观看| 大香蕉久久网| 全区人妻精品视频| 日本免费一区二区三区高清不卡| 日韩一本色道免费dvd| 亚洲精品自拍成人| 亚洲乱码一区二区免费版| 国产不卡一卡二| h日本视频在线播放| 少妇高潮的动态图| 色5月婷婷丁香| 中文亚洲av片在线观看爽| 亚洲成人精品中文字幕电影| 成人国产麻豆网| 国产一区二区激情短视频| 嫩草影院精品99| 一个人观看的视频www高清免费观看| 国产麻豆成人av免费视频| 成人美女网站在线观看视频| 天天一区二区日本电影三级| 中文资源天堂在线| 国产极品精品免费视频能看的| 人人妻人人澡欧美一区二区| 99久久精品热视频| 国产成年人精品一区二区| 日韩av在线大香蕉| 久久久久性生活片| 男女那种视频在线观看| h日本视频在线播放| 人人妻人人澡欧美一区二区| 一边亲一边摸免费视频| 少妇猛男粗大的猛烈进出视频 | 亚洲精品国产av成人精品| 日本撒尿小便嘘嘘汇集6| 最近最新中文字幕大全电影3| 久久久欧美国产精品| 伦理电影大哥的女人| 悠悠久久av| 亚洲真实伦在线观看| 九九久久精品国产亚洲av麻豆| 爱豆传媒免费全集在线观看| 亚洲精品国产成人久久av| 日本免费a在线| 亚洲va在线va天堂va国产| 男人狂女人下面高潮的视频| 亚洲图色成人| 麻豆成人av视频| 在线a可以看的网站| 中文字幕av成人在线电影| 欧美日韩在线观看h| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 波多野结衣巨乳人妻| 欧美日本视频| 三级经典国产精品| 男的添女的下面高潮视频| 欧美日韩在线观看h| 亚洲精品色激情综合| 九色成人免费人妻av| 国产黄片视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美 国产精品| 狂野欧美激情性xxxx在线观看| 人妻久久中文字幕网| 国内精品美女久久久久久| 欧美一区二区亚洲| 校园春色视频在线观看| 1000部很黄的大片| 综合色丁香网| 午夜久久久久精精品| 成人毛片60女人毛片免费| 国产一区二区激情短视频| 成熟少妇高潮喷水视频| 亚洲一区二区三区色噜噜| 欧美性猛交╳xxx乱大交人| 成人亚洲精品av一区二区| 国产激情偷乱视频一区二区| 日本与韩国留学比较| av视频在线观看入口| 欧美一级a爱片免费观看看| 欧美三级亚洲精品| 麻豆国产97在线/欧美| 免费在线观看成人毛片| 亚洲国产精品成人综合色| 男插女下体视频免费在线播放| 精品日产1卡2卡| 亚洲中文字幕日韩| 国产一区二区在线av高清观看| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频| ponron亚洲| 男人舔奶头视频| 嫩草影院精品99| 秋霞在线观看毛片| 伦精品一区二区三区| 热99在线观看视频| 日韩 亚洲 欧美在线| 免费看光身美女| 黄片无遮挡物在线观看| 一个人免费在线观看电影| 久久久久久久久久久丰满| 午夜免费男女啪啪视频观看| 国产av在哪里看| 少妇裸体淫交视频免费看高清| 婷婷色av中文字幕| 天堂av国产一区二区熟女人妻| 悠悠久久av| 国国产精品蜜臀av免费| 99久国产av精品国产电影| 看片在线看免费视频| 男的添女的下面高潮视频| 乱人视频在线观看| 人妻久久中文字幕网| 寂寞人妻少妇视频99o| 超碰av人人做人人爽久久| 国产精品永久免费网站| 成人毛片60女人毛片免费| 淫秽高清视频在线观看| 欧美成人精品欧美一级黄| 婷婷精品国产亚洲av| 精品国产三级普通话版| 免费观看的影片在线观看| 亚洲天堂国产精品一区在线| 欧美不卡视频在线免费观看| 搡女人真爽免费视频火全软件| 亚洲精品影视一区二区三区av| 成人鲁丝片一二三区免费| 亚洲精华国产精华液的使用体验 | 一边亲一边摸免费视频| 亚洲精品国产av成人精品| 亚洲国产精品国产精品| 国产精品99久久久久久久久| 日韩人妻高清精品专区| 国产男人的电影天堂91| 亚洲婷婷狠狠爱综合网| 亚洲av不卡在线观看| 久久久久免费精品人妻一区二区| 一区二区三区高清视频在线| 国产精品久久久久久精品电影| 高清毛片免费观看视频网站| 又爽又黄a免费视频| ponron亚洲| 亚洲国产欧美人成| 春色校园在线视频观看| 99久久成人亚洲精品观看| 日本熟妇午夜| 国产成人91sexporn| 亚洲国产精品成人综合色| 女人被狂操c到高潮| 99久久精品一区二区三区| 久久热精品热| 国产成年人精品一区二区| 最后的刺客免费高清国语| 在线播放国产精品三级| 精品熟女少妇av免费看| 免费大片18禁| 天堂影院成人在线观看| 亚洲精品久久久久久婷婷小说 | 国产精品一区二区性色av| 国产成年人精品一区二区| 亚洲最大成人手机在线| 啦啦啦观看免费观看视频高清| 日韩一本色道免费dvd| 听说在线观看完整版免费高清| 人妻系列 视频| 久久久精品大字幕| 在线观看免费视频日本深夜| 少妇被粗大猛烈的视频| 免费大片18禁| 性欧美人与动物交配| 国产精品不卡视频一区二区| 久久精品人妻少妇| 综合色av麻豆| 51国产日韩欧美| 国产亚洲精品av在线| 国产女主播在线喷水免费视频网站 | 在线观看免费视频日本深夜| 波多野结衣高清无吗| 日韩成人av中文字幕在线观看| 欧美成人免费av一区二区三区| 18+在线观看网站| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 美女高潮的动态| 国产白丝娇喘喷水9色精品| 爱豆传媒免费全集在线观看| 熟女人妻精品中文字幕| 国产伦理片在线播放av一区 | 免费av毛片视频| 51国产日韩欧美| 看免费成人av毛片| 日韩中字成人| 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 国产 一区精品| 国产真实乱freesex| 国产探花在线观看一区二区| 日本与韩国留学比较| 好男人在线观看高清免费视频| 国产精品99久久久久久久久| 18禁黄网站禁片免费观看直播| 天天躁日日操中文字幕| 波多野结衣巨乳人妻| 只有这里有精品99| 亚洲最大成人手机在线| 看十八女毛片水多多多| 91久久精品国产一区二区三区| 午夜久久久久精精品| 亚洲精品亚洲一区二区| 欧美xxxx黑人xx丫x性爽| 久久久a久久爽久久v久久| 天天躁夜夜躁狠狠久久av| 综合色丁香网| 欧美激情在线99| 亚洲国产精品国产精品| 五月伊人婷婷丁香| 亚洲久久久久久中文字幕| 日韩欧美一区二区三区在线观看| 亚洲精品456在线播放app| 美女大奶头视频| 精品日产1卡2卡| 哪个播放器可以免费观看大片| 久久草成人影院| 男女下面进入的视频免费午夜| 午夜精品国产一区二区电影 | 色综合站精品国产| 久久亚洲国产成人精品v| 一级二级三级毛片免费看| 搞女人的毛片| 一个人观看的视频www高清免费观看| av福利片在线观看| 亚洲四区av| 一进一出抽搐gif免费好疼| 搡老妇女老女人老熟妇| 亚洲国产精品国产精品| 看黄色毛片网站| 国产在视频线在精品| 国产亚洲5aaaaa淫片| 麻豆国产av国片精品| 美女脱内裤让男人舔精品视频 | 国产精品久久视频播放| 联通29元200g的流量卡| 深夜a级毛片| 人妻少妇偷人精品九色| 国语自产精品视频在线第100页| a级毛片免费高清观看在线播放| 久久精品国产亚洲av天美| 日本与韩国留学比较| 亚洲av免费高清在线观看| 国产伦理片在线播放av一区 | 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 波多野结衣高清作品| 日本三级黄在线观看| 亚洲精品久久久久久婷婷小说 | 国产高清视频在线观看网站| 精品日产1卡2卡| 村上凉子中文字幕在线| 麻豆成人午夜福利视频| 日本黄大片高清| 校园人妻丝袜中文字幕| 啦啦啦啦在线视频资源| 中文字幕精品亚洲无线码一区| 大香蕉久久网| 九色成人免费人妻av| 蜜桃久久精品国产亚洲av| 亚洲精品久久久久久婷婷小说 | 精品一区二区免费观看| 久久久国产成人免费| 波野结衣二区三区在线| av在线播放精品| 国产免费男女视频| 91久久精品国产一区二区三区| 99热这里只有是精品50| 国产亚洲精品久久久久久毛片| 国产一区二区三区av在线 | 国产精品一区二区在线观看99 | 国产亚洲av嫩草精品影院| 亚洲人成网站在线观看播放| 99在线人妻在线中文字幕| 老师上课跳d突然被开到最大视频| 高清日韩中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 2021天堂中文幕一二区在线观| 久久精品夜色国产| 12—13女人毛片做爰片一| 亚洲va在线va天堂va国产| 国产黄片美女视频| 久久久精品欧美日韩精品| 精品久久久噜噜| 久久精品国产自在天天线| 久久中文看片网| 亚洲aⅴ乱码一区二区在线播放| 亚洲,欧美,日韩| 国产私拍福利视频在线观看| 国产v大片淫在线免费观看| 永久网站在线| 18禁在线无遮挡免费观看视频| 伦精品一区二区三区| 乱人视频在线观看| 亚洲成av人片在线播放无| 天堂√8在线中文| 亚洲人成网站在线播| 日日啪夜夜撸| 国产精品一区www在线观看| 熟女电影av网| 夜夜爽天天搞| 亚洲国产日韩欧美精品在线观看| 久久人人爽人人爽人人片va| 一本精品99久久精品77| 3wmmmm亚洲av在线观看| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99 | 嫩草影院入口| 三级经典国产精品| 激情 狠狠 欧美| a级一级毛片免费在线观看| 免费av毛片视频| 午夜精品国产一区二区电影 | 麻豆成人av视频| av在线天堂中文字幕| 亚洲欧美成人精品一区二区| 97在线视频观看| 国产精品野战在线观看| 亚洲一级一片aⅴ在线观看| 国产在视频线在精品| 婷婷色综合大香蕉| 非洲黑人性xxxx精品又粗又长| 欧美三级亚洲精品| 美女内射精品一级片tv| 午夜视频国产福利| 国产精品免费一区二区三区在线| 成人鲁丝片一二三区免费| 国产免费一级a男人的天堂| av卡一久久| 看黄色毛片网站| 99久久中文字幕三级久久日本| 亚洲成a人片在线一区二区| 国产亚洲精品久久久久久毛片| 午夜亚洲福利在线播放| 午夜精品国产一区二区电影 | 国产精品一区www在线观看| 白带黄色成豆腐渣| or卡值多少钱| 国产日韩欧美在线精品| www日本黄色视频网| 久久午夜亚洲精品久久| 老熟妇乱子伦视频在线观看| 国产黄片视频在线免费观看| 精品少妇黑人巨大在线播放 | 精品久久国产蜜桃| 尤物成人国产欧美一区二区三区| 国产一区二区三区av在线 | 国产老妇女一区| 淫秽高清视频在线观看| 日韩欧美国产在线观看| 婷婷六月久久综合丁香| 色尼玛亚洲综合影院| av免费观看日本| 国产真实乱freesex| 在线观看av片永久免费下载| 午夜免费激情av| 亚洲色图av天堂| 最近中文字幕高清免费大全6| 尾随美女入室| 亚洲精品色激情综合| 如何舔出高潮| 自拍偷自拍亚洲精品老妇| 乱系列少妇在线播放| 亚洲中文字幕日韩| 午夜老司机福利剧场| 成年版毛片免费区| 91在线精品国自产拍蜜月| 欧美高清成人免费视频www| 淫秽高清视频在线观看| 亚洲丝袜综合中文字幕| 欧美精品国产亚洲| 老女人水多毛片| 日日撸夜夜添| 精品人妻一区二区三区麻豆| 亚洲美女视频黄频| 又爽又黄无遮挡网站| 2022亚洲国产成人精品| 久久精品国产亚洲网站| 日韩一本色道免费dvd| 简卡轻食公司| 亚洲中文字幕日韩| 乱系列少妇在线播放| 日本-黄色视频高清免费观看| 国模一区二区三区四区视频| 婷婷精品国产亚洲av| 亚洲aⅴ乱码一区二区在线播放| 少妇高潮的动态图| 卡戴珊不雅视频在线播放| 欧美高清成人免费视频www| 国产精品国产三级国产av玫瑰| 久久久久久久久久久丰满| 少妇人妻一区二区三区视频| 91精品国产九色| 久久精品国产亚洲av香蕉五月| 国产精品永久免费网站| 男人舔女人下体高潮全视频| 一级av片app| 一区福利在线观看| 国产精品永久免费网站| 国内揄拍国产精品人妻在线| 成熟少妇高潮喷水视频| 在线观看一区二区三区| 欧美成人免费av一区二区三区| 国产 一区精品| 亚洲美女搞黄在线观看| 亚洲成人中文字幕在线播放| 国产精品蜜桃在线观看 | 少妇人妻一区二区三区视频| 国产一级毛片在线| 国产高清激情床上av| av在线老鸭窝| a级毛片a级免费在线| 国产探花极品一区二区| 毛片一级片免费看久久久久| 天堂中文最新版在线下载 | 一进一出抽搐动态| 日本与韩国留学比较| 中文字幕制服av| av免费在线看不卡| 黄色欧美视频在线观看| 免费看美女性在线毛片视频| 国产黄片美女视频|