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

    基于改進模態(tài)分析法的柔性多端交直流混聯(lián)系統(tǒng)靜態(tài)電壓穩(wěn)定性評估

    2021-09-16 10:14:38郭小穎唐俊杰謝開貴
    電工技術(shù)學報 2021年17期
    關(guān)鍵詞:交直流換流站直流

    郭小穎 唐俊杰 舒 銅 彭 穗 謝開貴

    (1. 輸配電裝備及系統(tǒng)安全與新技術(shù)國家重點實驗室(重慶大學) 重慶 400044 2. 湖北省電力勘測設(shè)計院有限公司 武漢 430040 3. 廣東電網(wǎng)有限責任公司電網(wǎng)規(guī)劃研究中心 廣州 510699)

    0 引言

    電壓源型高壓直流輸電(Voltage Source Converter-High Voltage Direct Current, VSC-HVDC)技術(shù)相比于傳統(tǒng)的電流源型高壓直流輸電(Line Current Converter-High Voltage Direct Current, LCCHVDC)技術(shù)有著控制靈活、無換相失敗等優(yōu)點,它可以實現(xiàn)向無源網(wǎng)絡(luò)供電及多落點受電的功能,且功率定向可控。VSC換流器靈活的控制方式更有利于形成多端交直流混聯(lián)系統(tǒng)。當多端直流電網(wǎng)通過VSC換流站接入傳統(tǒng)的交流電網(wǎng)后,會對交流電網(wǎng)的電壓穩(wěn)定性造成巨大的影響,且直流電網(wǎng)本身也存在電壓穩(wěn)定性的問題。因此,研究AC/VSC-MTDC混聯(lián)系統(tǒng)的電壓穩(wěn)定性問題成為當前的焦點之一,其對保證交直流混聯(lián)系統(tǒng)大電網(wǎng)安全可靠的運行具有極其重要的意義[1-6]。

    文獻[7-9]采用靈敏度法對電壓穩(wěn)定性問題進行了分析研究。目前靈敏度法主要應(yīng)用于純交流系統(tǒng)的電壓穩(wěn)定性問題。該方法利用輸入變量與目標變量之間的偏微分關(guān)系求得靈敏度指標,計算簡單且物理概念明確;但該方法的參考價值嚴重依賴所選變量的偏微分靈敏度公式。在實際應(yīng)用中,難以運用變量之間的偏微分靈敏度指標全面表征系統(tǒng)的運行情況,所以靈敏度法存在一定的局限性。

    文獻[10-13]采用連續(xù)潮流法(Continuation Power Flow, CPF)對交流電網(wǎng)進行電壓穩(wěn)定性分析。該方法主要用于計算純交流系統(tǒng)當前運行點至電壓穩(wěn)定臨界點之間的最大負荷裕度。文獻[14]將CPF應(yīng)用于VSC-MTDC交直流混聯(lián)系統(tǒng)的電壓穩(wěn)定性研究中,但該方法主要通過模擬負荷增長求解電壓穩(wěn)定裕度,無法兼顧其他故障場景,如直流電網(wǎng)線路故障;也不能精確定位系統(tǒng)中影響電壓穩(wěn)定的關(guān)鍵節(jié)點或線路,且計算速度較慢。

    文獻[15-17]中將短路比指標應(yīng)用于交直流系統(tǒng)中,但大部分短路比指標的表達式中,均有部分變量取值于系統(tǒng)額定運行工況。而實際運行中,系統(tǒng)參數(shù)是實時變化的(如直流閉鎖后的緊急功率恢復過程),因此,這些指標不能及時、準確地反映系統(tǒng)的變化,導致其常應(yīng)用于交直流混聯(lián)系統(tǒng)的規(guī)劃和設(shè)計階段,在指導實際電力系統(tǒng)安全穩(wěn)定運行時準確性不足。

    文獻[18]應(yīng)用模態(tài)分析法對電壓穩(wěn)定性進行分析。模態(tài)分析法能夠判斷系統(tǒng)在某一運行狀態(tài)下的電壓穩(wěn)定情況,得到系統(tǒng)當前運行點距電壓不穩(wěn)定點之間的距離,并找出影響系統(tǒng)電壓穩(wěn)定的關(guān)鍵節(jié)點和線路。因此,模態(tài)分析法是靜態(tài)電壓穩(wěn)定性分析方法中常用的一種。但是,傳統(tǒng)的模態(tài)分析法是基于降階Q-V雅可比矩陣JR計算的,忽略了系統(tǒng)有功功率對電壓穩(wěn)定的影響[18];文獻[19]雖然對傳統(tǒng)的模態(tài)分析法進行了改進,但也只在系統(tǒng)PQ節(jié)點達到無功功率限制時考慮有功功率對電壓穩(wěn)定的影響;文獻[20]對模態(tài)分析法進一步改進,基于交流系統(tǒng)全雅可比矩陣進行分析,提出一種新的狀態(tài)-模態(tài)指標,能夠量化表征系統(tǒng)狀態(tài)對臨界模態(tài)的貢獻,但該文的分析對象是純交流系統(tǒng)的電壓穩(wěn)定性;文獻[21]將模態(tài)分析法應(yīng)用于交直流系統(tǒng),但針對的是基于LCC(line current converter)技術(shù)的交直流混聯(lián)系統(tǒng),且僅考慮了一條直流線路,直流比例低且未能成網(wǎng)。

    因此,本文旨在對傳統(tǒng)模態(tài)分析法進行改進,并提出將交直流混聯(lián)系統(tǒng)進行解耦分析。解耦之后,交流電網(wǎng)基于其全雅可比矩陣而非降階Q-V矩陣進行計算與模態(tài)分析,所以考慮了有功功率對電壓穩(wěn)定性的影響;而直流電網(wǎng)由于僅傳輸有功功率,改進后的模態(tài)分析法能夠推廣應(yīng)用到直流系統(tǒng)的電壓穩(wěn)定性評估過程中。因此,本文將改進的模態(tài)分析法推廣應(yīng)用到AC/VSC-MTDC混聯(lián)系統(tǒng)的靜態(tài)電壓穩(wěn)定性評估領(lǐng)域。并且,本文將狀態(tài)-模態(tài)指標推廣應(yīng)用到交直流系統(tǒng)中,對影響電壓穩(wěn)定性的關(guān)鍵節(jié)點和線路進行準確辨識,能夠分析交流、直流系統(tǒng)負荷增長及線路故障對交直流混聯(lián)系統(tǒng)電壓穩(wěn)定性的影響。此外,通過改變換流站的控制功率和換流站在交流電網(wǎng)的連接位置,能夠量化分析換流站、交流電網(wǎng)及直流電網(wǎng)之間的相互影響,從而進一步量化換流站作為連接交流電網(wǎng)與直流電網(wǎng)的“橋梁”所發(fā)揮的關(guān)鍵作用,進而改善整個交直流混聯(lián)系統(tǒng)的電壓穩(wěn)定性。

    1 VSC換流站模型建立

    當前許多研究在對VSC換流站建模時,沒有考慮換流站的損耗或是直接省略了濾波器元件[22-23],這不符合實際,會對潮流計算結(jié)果產(chǎn)生一定的影響。因此,本文采用精確計及VSC換流站損耗的模型。VSC可采用兩電平、三電平中點鉗位或模塊化多電平換流器結(jié)構(gòu)。換流站直流側(cè)采用雙極接線,連接到直流母線。交流電網(wǎng)通過換流變壓器、濾波器和換相電抗器連接在VSC換流器上。VSC換流站的結(jié)構(gòu)如圖1所示。

    圖1 VSC換流站的結(jié)構(gòu)Fig.1 Diagram for the structure of VSC station

    圖1 中,換流變壓器的等值阻抗為Zt=Rtf+jXtf,相應(yīng)的導納1/Zt=Gtf+jBtf;電抗器阻抗為Zc=Rcf+jXcf,相應(yīng)的導納1/Zc=Gcf+jBcf;濾波器的電納為jBf;Us∠δs表示交流節(jié)點電壓;Uf∠δf表示濾波器節(jié)點的電壓;Uc∠δc為VSC換流站交流側(cè)的電壓。

    以圖1中的功率方向為正方向,則VSC換流站注入到交流電網(wǎng)的有功和無功功率則可分別表示為

    VSC換流站交流側(cè)的有功和無功功率可以分別表示為

    本文采用精確計及換流站損耗的模型,具體的損耗為流過換流站電流Ic的二次函數(shù),即

    式中,a、b、c為損耗系數(shù),通常從實際工程中獲取,本文采用的換流站及直流電網(wǎng)參數(shù)可參見文獻[24]。

    2 交直流系統(tǒng)雅可比矩陣建立

    電力系統(tǒng)的潮流方程可分為以節(jié)點注入功率、電流和電壓為網(wǎng)絡(luò)注入量的三種類型。根據(jù)節(jié)點注入功率,可以建立一組非線性潮流方程。為了計算方便,本文將交直流混聯(lián)系統(tǒng)的潮流模型進行統(tǒng)一表達,主要分為三個部分:交流電網(wǎng)、換流站和直流電網(wǎng)。進而采用牛頓拉夫遜法進行迭代求解,并取最后一次迭代收斂的雅可比矩陣用于后續(xù)的模態(tài)分析。

    交流電網(wǎng)、換流站以及直流電網(wǎng)的功率失配方程如下[24]。

    交流電網(wǎng)功率失配方程為

    式中,XAC為交流電網(wǎng)的狀態(tài)變量向量,XAC=[δTUT]T,δ為交流電網(wǎng)所有非平衡節(jié)點的電壓相角,U為交流電網(wǎng)所有PQ節(jié)點電壓幅值;ΔP為所有PQ節(jié)點和PV節(jié)點有功不平衡量;ΔQ為所有PQ節(jié)點的無功不平衡量。

    VSC換流站功率失配方程為

    式中,XVSC為換流站的狀態(tài)變量向量,本文VSC換流站采用單點電壓控制,故XVSC=[δTfUTfδTcUTc]T;ΔPf和ΔQf為濾波器節(jié)點的功率不平衡量;ΔPs和ΔQs為VSC換流器注入到交流電網(wǎng)的功率不平衡量。

    直流電網(wǎng)的失配方程為

    式中,XDC為直流電網(wǎng)的狀態(tài)變量向量為所有直流節(jié)點的電壓幅值向量;ΔPdc為直流節(jié)點的有功不平衡量。

    由上述可知,式(6)~式(8)組成了交直流系統(tǒng)潮流計算的非線性方程組。其中所有待求電壓幅值的初值設(shè)為1.0(pu),相位初值設(shè)為0°。采用牛頓拉夫遜法,對非線性方程組式(6)~式(8)進行迭代求解,忽略泰勒展開式的高階項,即可得到相應(yīng)的修正方程式為

    為了方便后續(xù)推導計算,將式(8)簡記為F(X)=JΔX,其中J為交直流系統(tǒng)的雅可比矩陣。

    3 交直流系統(tǒng)靜態(tài)電壓性穩(wěn)定分析

    交流電網(wǎng)的失配方程中不包括直流電網(wǎng)的變量,同樣地,直流電網(wǎng)的失配方程也不包括交流電網(wǎng)的變量[24]。而交流、直流電網(wǎng)的失配方程同時都包含了換流站的變量。由式(9)中雅可比矩陣J的元素也可以看出,?FAC/?XDC=0,?FDC/?XAC=0,這說明交流電網(wǎng)中控制變量的改變不會直接引起直流電網(wǎng)變量的改變;反之亦然。二者之間的相互作用是通過換流站傳遞的。換流站通常分為整流站和逆變站,整流站從交流電網(wǎng)吸收有功功率,注入到直流電網(wǎng);逆變站則從直流電網(wǎng)吸收有功功率,注入到交流電網(wǎng)。這兩種類型的換流站對交流電網(wǎng)和直流電網(wǎng)的電壓穩(wěn)定性影響是完全相反的。因此,將傳統(tǒng)的模態(tài)分析法直接套用在整個交直流系統(tǒng)是不可取的,必須對傳統(tǒng)的模態(tài)分析法進行改進。因此,本文提出對交直流混聯(lián)系統(tǒng)進行解耦。下面將證明對交直流系統(tǒng)進行解耦分析的可行性。

    在傳統(tǒng)的模態(tài)分析法中,若系統(tǒng)的降階Q-V雅可比矩陣奇異,則判斷系統(tǒng)電壓失穩(wěn)[18]。同樣地,在交直流混聯(lián)系統(tǒng)中,若交直流系統(tǒng)的雅可比矩陣J奇異,則整個交直流混聯(lián)系統(tǒng)電壓失穩(wěn)。但是,由于不同類型的換流站對交流電網(wǎng)和直流電網(wǎng)的電壓穩(wěn)定性影響是完全相反的,因此,直接對交直流系統(tǒng)的雅可比矩陣J進行模態(tài)分析是不可取的。本文對交流分塊矩陣JAC和直流分塊矩陣JDC分別進行模態(tài)分析。接下來將證明JAC和JDC其中之一出現(xiàn)奇異,則J也會奇異[25]。

    根據(jù)式(9),對矩陣J求逆,可以推導出交直流混聯(lián)系統(tǒng)的狀態(tài)變量的修正量表達式為

    其中

    由式(10)可以得出一些重要結(jié)論:

    (1)交直流系統(tǒng)雅可比矩陣J的逆矩陣在很大程度上同時取決于交流分塊矩陣JAC和直流分塊矩陣JDC。

    (2)JAC和JDC是分別代表交直流系統(tǒng)在交流電網(wǎng)和直流電網(wǎng)的矩陣,如果這兩個矩陣中有一個矩陣奇異,則不可能求出J的逆,故J也會奇異。

    上述結(jié)果說明JAC和JDC均包含了足夠的信息來表征交流、直流系統(tǒng)的電壓穩(wěn)定性;而整體交直流混聯(lián)系統(tǒng)中,一旦交流電網(wǎng)或直流電網(wǎng)任何一方失穩(wěn),都會導致整個交直流系統(tǒng)電壓失穩(wěn),從而在理論上保證了將交直流系統(tǒng)解耦分析的可行性。

    3.1 交流電網(wǎng)電壓穩(wěn)定性分析

    將交直流系統(tǒng)中的交流電網(wǎng)修正方程展開為

    與傳統(tǒng)的模態(tài)分析法不同,本文基于整個交流分塊矩陣JAC對交流電網(wǎng)進行靜態(tài)電壓穩(wěn)定分析。因此,能夠考慮有功功率對電壓穩(wěn)定性的影響。

    對交流電網(wǎng)矩陣JAC進行模態(tài)分析,JAC為

    則特征值矩陣可進一步表示為

    相應(yīng)地,交流電網(wǎng)的最小特征值可以表示為

    式中,ηmin、ξmin分別為λmin對應(yīng)的左右特征向量。與傳統(tǒng)模態(tài)分析法中的參與因子指標不同,本文采用一種新的分析視角,通過對JAC進行分解,將交流電網(wǎng)的最小特征值λmin表示為系統(tǒng)狀態(tài)的加權(quán)和,從而產(chǎn)生了狀態(tài)-模態(tài)參與因子(State-in-Mode Participation Factor, SIMPF)和狀態(tài)-模態(tài)靈敏度(State-in-Mode Sensitivity, SIMS)的定義[20]。SIMPF可以量化系統(tǒng)狀態(tài)對臨界模態(tài)的直接貢獻,而SIMS則用于量化系統(tǒng)狀態(tài)對臨界模態(tài)的控制靈敏度。相應(yīng)地,將λmin重新改寫為[20]

    λmin由三個部分組成,Cnet、Cinj、Cline越小,λmin就越小,系統(tǒng)則越趨于不穩(wěn)定,現(xiàn)對其進行具體分析。

    Cnet表示電力網(wǎng)絡(luò)耦合強度對最小特征值的貢獻,且Cnet恒為正[20],表明網(wǎng)絡(luò)耦合對電壓穩(wěn)定性有積極影響,這與我們的實際經(jīng)驗一致[20]。因此,臨界模態(tài)的形成主要取決于Cinj和Cline。Cinj量化了節(jié)點注入功率對最小特征值的貢獻大?。欢鳦line則體現(xiàn)了線路傳輸功率對最小特征值的貢獻。

    將λmin對應(yīng)的左右特征向量ηmin、ξmin的有功、無功部分分別記為ηmp、ηmq、ξmp、ξmq,i∈vpq和j∈vpv分別表示交流節(jié)點i、j的類型,則Cinj可具體展開為[20]

    同理,可將Cline具體展開為[18]

    根據(jù)式(16),定義節(jié)點i的有功功率和無功功率的SIMPF指標為

    式中,SPi和SQi分別為相應(yīng)的SIMS指標,將其分別展開為

    一般而言,可認為Pij≈Pji,因此,根據(jù)式(17),可定義表征線路Lij之間傳輸功率的SIMPF指標為

    Sij為線路Lij的SIMS指標,將其展開為

    SIMPF和SIMS分別衡量系統(tǒng)狀態(tài)對臨界模態(tài)的直接貢獻大小和控制靈敏度。因此,具有負SIMPF的系統(tǒng)狀態(tài)不利于電壓穩(wěn)定,具有正SIMPF的系統(tǒng)狀態(tài)有利于電壓穩(wěn)定。而SIMS指標的符號則指明了調(diào)整相應(yīng)系統(tǒng)狀態(tài)以改善靜態(tài)電壓穩(wěn)定性的方向。且SIMS的絕對值越大,同樣的狀態(tài)調(diào)整量可產(chǎn)生更好的效果。

    相比傳統(tǒng)模態(tài)分析法中參與因子(Mode-in-State Participation Factor, MISPF)指標[20],本文所采用的狀態(tài)-模態(tài)指標具有以下優(yōu)點:

    (1)MISPF是一種模態(tài)-狀態(tài)指標,只能表征在臨界模態(tài)下,系統(tǒng)節(jié)點的無功功率對該節(jié)點電壓幅值的影響。而狀態(tài)-模態(tài)指標將臨界模態(tài)表示為系統(tǒng)狀態(tài)的加權(quán)和,能夠具體量化系統(tǒng)狀態(tài)對臨界模態(tài)形成的直接影響。

    (2)MISPF指標基于純交流系統(tǒng)中降階Q-V雅可比矩陣JR進行分析,假定各節(jié)點的有功功率不變,只考慮節(jié)點的無功功率對電壓穩(wěn)定的影響。這一近似處理在距系統(tǒng)臨界運行點較遠時是成立的。但隨著負荷加重,系統(tǒng)逼近臨界運行點時,節(jié)點的無功功率逐漸達到其上限,此時,節(jié)點有功功率的變化無法忽略,假設(shè)ΔP=0是不符合實際的。因此,MISPF指標應(yīng)用在臨界運行點時存在一定局限。而狀態(tài)-模態(tài)指標基于整個交流分塊雅可比矩陣JAC分析,能夠同時考慮節(jié)點注入有功、無功功率對電壓穩(wěn)定性的影響,在實際應(yīng)用中更為全面、合理。

    (3)狀態(tài)-模態(tài)指標能夠計及有功功率對電壓穩(wěn)定性的影響,而直流部分只傳輸有功功率,因此,狀態(tài)-模態(tài)指標能夠推廣應(yīng)用到直流電網(wǎng),并且能夠具體量化分析直流電網(wǎng)與交流電網(wǎng)之間的相互影響。

    3.2 MISPF與SIMS的區(qū)別與聯(lián)系

    無論是運用本文方法對整個交流分塊雅可比矩陣JAC進行模態(tài)分析,還是運用傳統(tǒng)的模態(tài)分析法對純交流系統(tǒng)的降階Q-V雅可比矩陣JR進行模態(tài)分析,當系統(tǒng)處于臨界模態(tài)時,系統(tǒng)的最小特征值均為0。因此有JACξmin=0和ηTminJAC=0,其中ηmin和ξmin分別為JAC的臨界特征值對應(yīng)的左右特征向量;同理,也有JRξq=0和ηTqJR=0,ηq和ξq分別為JR的臨界特征值對應(yīng)的左右特征向量。因此,JR在臨界模態(tài)時的左右特征向量ηq和ξq實際上是JAC的臨界特征向量ηmin和ξmin的無功部分[20]。根據(jù)文獻[18],節(jié)點i在臨界模態(tài)時的MISPF指標為ηqiξqi,該指標與式(21)中的SQi成正比。所以,MISPF在數(shù)值上等于SQi/2。但傳統(tǒng)的MISPF指標側(cè)重于臨界模態(tài)對系統(tǒng)狀態(tài)的影響,該指標只能表征在臨界模態(tài)下,系統(tǒng)節(jié)點的無功功率對該節(jié)點電壓幅值的影響;而本文采用的SIMS指標從狀態(tài)-模態(tài)的角度出發(fā),將系統(tǒng)狀態(tài)表示為臨界模態(tài)的加權(quán)和,因此能夠量化系統(tǒng)狀態(tài)對臨界模態(tài)形成的直接貢獻。

    3.3 直流電網(wǎng)電壓穩(wěn)定性分析

    狀態(tài)-模態(tài)指標SIMPF和SIMS能夠計及有功功率對系統(tǒng)電壓穩(wěn)定的影響,而直流電網(wǎng)只傳輸有功功率,因此,本文將SIMPF和SIMS指標推廣應(yīng)用于直流電網(wǎng)的電壓穩(wěn)定性評估。

    設(shè)直流部分共有ndc個直流節(jié)點,與交流部分類似,直流部分的修正方程式為

    同理,將直流雅可比矩陣的最小特征值λcri用直流電網(wǎng)系統(tǒng)狀態(tài)表示為

    式中,ηcri和ξcri分別為直流電網(wǎng)最小特征值λcri對應(yīng)的左、右特征向量。

    同樣地,將直流部分的雅可比矩陣JDC進行分解,JDC可表示為[26]

    不同于交流部分的雅可比矩陣,直流雅可比矩陣JDC中元素只有?Pdc/?Udc,且直流導納矩陣Y中只有實部G而沒有虛部B,因此JDC分解后只含有兩部分,具體為

    其中

    實際中,電力系統(tǒng)多為微損系統(tǒng),因此JDC的不對稱程度較低又由于恒為正,因此恒為正。由此可知,λcri的變小甚至變負主要由決定,現(xiàn)將其展開為

    式中,i為直流的節(jié)點,i∈ndc。

    因此,直流節(jié)點i的SIMPF指標可以定義為

    相應(yīng)的控制靈敏度指標SIMS則為

    SIMPF指標反映了節(jié)點注入有功功率對直流系統(tǒng)最小特征值的直接貢獻,而SIMS則定位了影響直流系統(tǒng)電壓穩(wěn)定的關(guān)鍵節(jié)點。

    4 算例

    本文采用經(jīng)修改的IEEE 57節(jié)點交直流混聯(lián)系統(tǒng)[24],并在文獻[24]算例的第2、3、4、7個直流節(jié)點分別增加了0.1(pu)的直流負荷,同時將2號換流站的連接位置從交流12節(jié)點更換至13節(jié)點。如圖2所示,換流站分別通過1、13、16、17、15交流節(jié)點接入交流系統(tǒng)。

    圖2 改進IEEE 57節(jié)點7端VSC-MTDC系統(tǒng)示意圖Fig.2 Diagram for modified IEEE 57-bus system with 7-terminal VSC-MTDC

    換流站的控制方式為單點電壓控制,控制參數(shù)參見文獻[24]。潮流計算中,交直流混聯(lián)系統(tǒng)的基準功率設(shè)為100MV·A,直流電網(wǎng)的基準電壓設(shè)為200kV,交流電網(wǎng)基準電壓設(shè)為220kV,收斂的精度設(shè)為10-6。為了增強實驗效果,交流電網(wǎng)的原始負荷被提高到1.4倍,作為本算例的基本負荷水平。

    4.1 最小特征值指標的有效性驗證

    增加交流電網(wǎng)的負荷直至潮流計算結(jié)果不收斂,分別繪制交流、直流電網(wǎng)最小特征值指標隨著交流電網(wǎng)負荷增加的變化趨勢;而直流電網(wǎng)部分,由于直流電網(wǎng)初始的最小特征值指標λcri較大,即電壓穩(wěn)定裕度較高,因此,為了明確量化λcri的變化情況,直流負荷倍數(shù)需提高至很高的水平,才能使潮流計算結(jié)果不收斂。分別繪制交流、直流電網(wǎng)最小特征值指標隨直流電網(wǎng)負荷增加的變化趨勢,其結(jié)果如圖3所示。

    由圖3a可知,隨著交流電網(wǎng)負荷增加,交流電網(wǎng)最小特征值λmin逐漸減小,當系統(tǒng)潮流計算結(jié)果臨界收斂時,λmin接近于0,系統(tǒng)極不穩(wěn)定,而直流電網(wǎng)的最小特征值變化卻不大。由圖3b可知,隨著直流電網(wǎng)負荷增加,直流電網(wǎng)最小特征值λcri逐漸減小,當系統(tǒng)潮流計算結(jié)果臨界收斂時,λcri接近于0,系統(tǒng)的電壓趨于崩潰,而交流電網(wǎng)最小特征值幾乎不變。

    圖3 交直流電網(wǎng)最小特征值隨交直流負荷增長變化曲線Fig.3 AC-DC minimum eigenvalue change curves with increase of AC-DC load

    上述實驗結(jié)果表明,將交、直流最小特征值作為表征交直流系統(tǒng)失穩(wěn)的指標是可行的;同時也表明,在換流站控制功率不變的前提下,交流、直流系統(tǒng)無法傳遞對彼此之間的影響,二者之間的影響是間接的。這一結(jié)果也進一步驗證了本文將交直流系統(tǒng)進行解耦分析的可行性。

    4.2 交流電網(wǎng)電壓穩(wěn)定性分析

    4.2.1 交流電網(wǎng)關(guān)鍵節(jié)點識別

    Cinj指標能夠準確定位交流電網(wǎng)中對電壓穩(wěn)定性影響較大的關(guān)鍵節(jié)點。表征交流節(jié)點的有功、無功功率的SIMPF和SIMS指標見表1和表2。SIMPF指標表征系統(tǒng)狀態(tài)對最小特征值的直接貢獻:若SIMPF指標為負,則說明該系統(tǒng)狀態(tài)不利于電壓穩(wěn)定。而SIMS指標的符號則指明了調(diào)整相應(yīng)系統(tǒng)狀態(tài)以增強靜態(tài)電壓穩(wěn)定性的方向,且SIMS的絕對值越大,意味著同樣的狀態(tài)調(diào)整量可產(chǎn)生更好的效果。

    表1 SIMS指標排名靠前的節(jié)點注入有功率Tab.1 Node active power injection with high SIMS ranking

    表1 、表2分別列出了SIMS絕對值取值最大的前5個系統(tǒng)狀態(tài)。表中SIMPF指標為負,說明相應(yīng)的系統(tǒng)狀態(tài)不利于電壓穩(wěn)定,而SIMS指標為正,說明減小相應(yīng)的系統(tǒng)狀態(tài)能夠改善系統(tǒng)的電壓穩(wěn)定性,即需要對上述的負荷節(jié)點采取相應(yīng)的補償措施。在以上節(jié)點分別通過有功儲能裝備額外增加該節(jié)點的有功出力10MW、通過無功補償增加無功出力10Mvar。采取有功、無功補償措施后,交流電網(wǎng)的最小特征值λmin變化如圖4所示。

    表2 SIMS指標排名靠前的節(jié)點注入無功功率Tab.2 Node reactive power injection with high SIMS ranking

    圖4 交流電網(wǎng)最小特征值對比圖Fig.4 Comparison for the minimum eigenvalue of AC system

    分別在前五個關(guān)鍵節(jié)點采取有功、無功補償措施后,λmin從0.163 8增大到0.2左右,且SIMS指標越大,采取同樣的補償量,λmin增大程度越大,交流系統(tǒng)電壓穩(wěn)定性改善越好。作為對比,在表1中所涉及節(jié)點外的任一位置如44節(jié)點處增加10MW有功出力,λmin變?yōu)?.164 3;同理,在表2指標所涉及節(jié)點外的任一位置如19節(jié)點處增加10Mvar無功補償,λmin變?yōu)?.168 7。

    由圖4可知,在表1、表2所示關(guān)鍵節(jié)點之外的其他位置采取同樣的補償之后,λmin雖有增大,但較圖4中對關(guān)鍵節(jié)點采取補償而言,λmin增大幅度微弱。實驗結(jié)果表明,通過SIMS指標識別系統(tǒng)的關(guān)鍵節(jié)點,進而采取相應(yīng)的補償方案是有效的。

    4.2.2 交流電網(wǎng)關(guān)鍵線路識別

    Cline指標能夠準確定位交流電網(wǎng)中對其電壓穩(wěn)定性影響較大的關(guān)鍵線路,表3列出了SIMS絕對值取值最大的前5個線路狀態(tài)。為了驗證該指標識別關(guān)鍵線路的準確性,分別斷開上述關(guān)鍵的交流線路,交流電網(wǎng)的最小特征值λmin變化見表4。

    表3 SIMS指標排名高的交流線路Tab.3 AC lines with high SIMS ranking

    表4 交流線路開斷后的最小特征值Tab.4 Minimum eigenvalue after tripping off AC lines

    若線路L34,32故障斷開,此時潮流計算的結(jié)果不收斂,系統(tǒng)電壓失穩(wěn),說明線路L34,32對維持系統(tǒng)電壓穩(wěn)定具有重要意義;而線路L24,25為雙回線,若兩條線都發(fā)生故障斷開,此時潮流計算結(jié)果也不收斂,即系統(tǒng)失穩(wěn);若只有一回線路發(fā)生故障,λmin的下降程度也比其他線路大。

    實驗結(jié)果表明,Cline指標能夠?qū)涣麟娋W(wǎng)中的關(guān)鍵線路進行識別,且線路SIMS指標的絕對值越大,線路發(fā)生故障后,λmin下降程度越大,說明線路對維持交流電網(wǎng)的電壓穩(wěn)定性越重要。

    4.3 直流電網(wǎng)電壓穩(wěn)定分析

    4.3.1 直流負荷對直流電壓穩(wěn)定的影響

    同樣地,將直流節(jié)點的SIMS指標按其絕對值從大到小排列,結(jié)果見表5。

    表5 直流電網(wǎng)節(jié)點注入有功SIMS指標Tab.5 SIMS of node active power injection in DC Grid

    表5 中SIMPF指標為負,說明直流節(jié)點的有功負荷不利于電壓穩(wěn)定,而SIMS指標為正,則說明減小節(jié)點有功負荷有利于電壓穩(wěn)定,即要對相應(yīng)節(jié)點采取補償措施。在以上節(jié)點位置分別額外增加有功出力10 MW,采取補償措施后的直流系統(tǒng)最小特征值λcri變化趨勢如圖5所示。

    圖5 采取補償措施后直流最小特征值對比圖Fig.5 Comparison for DC minimum eigenvalue after compensation measure

    由圖5可知,分別在直流節(jié)點2、4、3、7處額外增加有功出力,λcri變大,直流電網(wǎng)電壓穩(wěn)定性得到改善,且SIMS指標越大,采取同樣的有功補償量,λcri的提高幅度越大。

    4.3.2 直流線路對直流電壓穩(wěn)定的影響

    由于直流雅可比矩陣中元素只有?Pdc/?Udc,且直流導納矩陣Y只有實部G,所以直流系統(tǒng)的狀態(tài)-模態(tài)指標中并沒有識別關(guān)鍵直流線路的指標。因此,本文采取模擬直流線路故障的方式,對直流線路的關(guān)鍵程度進行判斷。

    該算例中共有7條直流線路,斷開線路L17、L16和L56會使得潮流計算結(jié)果不收斂,結(jié)合圖2分析可知,這三條線路中的任一條線路故障都會使交流與直流系統(tǒng)解列,因此這三條線路對維持交直流系統(tǒng)的電壓穩(wěn)定性至關(guān)重要。而分別斷開剩余4條直流線路后,原本流過該線路的有功功率,就會被其他直流線路分擔,而導致這些線路承擔的傳輸功率比之前增加了,因此直流電網(wǎng)電壓穩(wěn)定性會下降。繪制直流電網(wǎng)最小特征值λcri隨直流線路故障的變化如圖6所示。當線路L37和L47發(fā)生故障時,λcri下降程度要比線路L24和L32大,這是由于線路L37和L47傳輸?shù)挠泄β时萀24和L32大,因此線路L37和L47對維持直流系統(tǒng)電壓穩(wěn)定性有著更為重要的作用??傊?,通過模擬直流線路故障的思路,可以對直流系統(tǒng)的關(guān)鍵線路進行識別。

    圖6 直流最小特征值變化趨勢Fig.6 Trend for DC minimum eigenvalue change

    4.4 換流站對交、直流電壓穩(wěn)定性的影響

    4.4.1 換流站控制功率對交、直流電壓穩(wěn)定性的影響

    由第3節(jié)可知,交直流混聯(lián)系統(tǒng)中,換流站與交流、直流系統(tǒng)直接相連,換流站的控制功率的增大(或減?。苯佑绊懡涣鳌⒅绷飨到y(tǒng)之間的電壓穩(wěn)定性。分別繪制直流、交流系統(tǒng)最小特征值隨換流站控制功率增加的趨勢如圖7、圖8所示。

    圖7 直流最小特征值隨換流站控制功率變化曲線Fig.7 DC minimum eigenvalue change curves with VSC control power

    圖8 交流最小特征值隨換流站控制功率變化曲線Fig.8 AC minimum eigenvalue change curves with VSC control power

    由圖7可知,隨著換流站控制功率的增加,直流系統(tǒng)的最小特征值λcri逐漸減小,當換流站的有功功率增加到原來的12.6倍時,λcri變負,此時直流系統(tǒng)中部分節(jié)點的電壓幅值已降至很低的水平,從而最終導致直流電網(wǎng)的電壓失穩(wěn)。

    如圖8所示,隨著換流站控制功率的增加,交流系統(tǒng)的最小特征值λmin先增大后減小。這是由于換流站控制功率增大后,逆變站注入到與其直接相連的交流節(jié)點的有功功率也會增大,會導致與這幾個交流節(jié)點直接相連的交流線路傳輸?shù)挠泄β实姆较虬l(fā)生改變,從而改變了交流線路的SIMPF指標。而SIMPF反映了線路傳輸功率對λmin的直接影響。因此,圖8中同時繪制了所有線路的SIMPF指標之和隨換流站控制功率的變化趨勢圖。由圖8可知,λmin的變化趨勢與SIMPF指標的變化趨勢是一致的,這說明換流站控制功率的變化通過改變交流系統(tǒng)線路之間的潮流分布,從而影響到交流系統(tǒng)的電壓穩(wěn)定性。盡管λmin先增大后減小,但最終都比VSC控制功率未改變時大。因此,交流系統(tǒng)的電壓穩(wěn)定性仍然得到改善。由圖7、圖8可知,換流站控制功率的改變,對交流、直流系統(tǒng)的電壓穩(wěn)定性的影響是完全相反的,這也進一步驗證了將交直流系統(tǒng)進行解耦分析的正確性。

    4.4.2 直流電網(wǎng)對交流電壓穩(wěn)定性的影響

    直流電網(wǎng)與交流電網(wǎng)通過換流站相連,二者之間的影響是間接的。直流電網(wǎng)通過直接影響VSC換流站,再由換流站影響交流節(jié)點的注入功率,從而影響交流電網(wǎng)的電壓穩(wěn)定性。因此,通過改變換流站在交流電網(wǎng)的連接位置,可以分析直流系統(tǒng)對交流系統(tǒng)的電壓穩(wěn)定性影響。 換流站與交流電網(wǎng)最初的連接位置如圖2所示,與換流站相連接的交流節(jié)點的有功SIMS與SIMPF指標見表6。

    表6 與換流站連接的交流節(jié)點有功指標Tab.6 Index for active power of AC bus connecting to the VSC

    對交流系統(tǒng)而言,逆變站相當于有功電源,逆變站接入交流系統(tǒng),可以改善交流系統(tǒng)的電壓穩(wěn)定性。根據(jù)表1、表2中的交流節(jié)點指標,將逆變站的位置分別更改為SIMS指標為前四的節(jié)點處,SIMS從大到小分別為交流節(jié)點31、33、32、30。繪制交流電網(wǎng)的最小特征值指標如圖9所示。

    由圖9可知,根據(jù)SIMS指標改變換流站在交流電網(wǎng)的連接位置,交流電網(wǎng)的最小特征值指標λmin明顯提高,均提高至0.3以上,交流電網(wǎng)的電壓穩(wěn)定性上升。作為對比,將換流站2的落點更換到除表1外的任意位置,如節(jié)點44,λmin從原來的0.163 8增加到0.194 6,提升幅度有限。

    圖9 改變VSC位置后交流最小特征值變化圖Fig.9 Diagram for minimum eigenvalue after relocation of VSC

    該結(jié)果表明,換流站在交流電網(wǎng)的連接位置不同,對交流電網(wǎng)的電壓穩(wěn)定性影響也不同。根據(jù)SIMS指標,改變換流站在交流電網(wǎng)的連接位置,不會對直流電網(wǎng)的電壓穩(wěn)定性造成影響,但卻能夠?qū)⒅绷麟娋W(wǎng)對交流電網(wǎng)的積極影響發(fā)揮到最大,從而改善整個交直流混聯(lián)電網(wǎng)的電壓穩(wěn)定性。這一應(yīng)用也可以為換流站在交流電網(wǎng)的落點選址提供依據(jù)。

    4.4.3 交流電網(wǎng)對直流電壓穩(wěn)定性的影響

    同樣地,交流系統(tǒng)對直流系統(tǒng)的影響也是間接的。在實際運行中,交流電網(wǎng)內(nèi)部功率傳輸?shù)淖兓?,會對與它相連的換流站節(jié)點處的功率提出不同的要求。也就是說,根據(jù)交流電網(wǎng)的不同運行情況,會對換流站的功率進行調(diào)整,而換流站的控制功率變化會進一步改變直流電網(wǎng)傳輸?shù)挠泄β?,進而引起直流系電網(wǎng)的電壓穩(wěn)定性變化。因此,通過改變換流站的功率分配,可以量化交流電網(wǎng)對直流電網(wǎng)的電壓穩(wěn)定性的影響。

    由圖2可知,該算例共有五個換流站,其中1號站為整流站,其余四個為逆變站,整流站與逆變站之間滿足有功功率平衡:即整流站向直流電網(wǎng)注入的有功功率等于四個逆變站從直流電網(wǎng)吸收的有功功率之和,總計200MW。在1號站總功率不變的情況下,可以根據(jù)交流電網(wǎng)的運行情況,將四個逆變站的控制功率進行重新分配。又因為與逆變站相連的四個直流節(jié)點的SIMPF指標全為負,因此,逆變站對直流電網(wǎng)的電壓穩(wěn)定性影響是不利的,且直流節(jié)點SIMS指標越大,該不利因素的影響越大。

    因此,在1號站總功率不變的情況下,連接節(jié)點2、3、4、5的逆變站的控制功率應(yīng)該與這幾個節(jié)點的SIMS指標大小成負相關(guān)進行分配;直流節(jié)點2、3、4的SIMS指標已在表4中展示,節(jié)點5的SIMS指標約為0。因此,對換流站的控制功率重新分配后,其功率大小應(yīng)滿足VSC2≤VSC4≤VSC3≤VSC5。改變逆變站的控制功率后,直流電網(wǎng)的最小特征值指標λcri見表7。

    表7 換流站控制功率對λcri的影響Tab.7 Influence of VSC control power on λcri

    根據(jù)與換流站直接相連的直流節(jié)點的SIMS指標,改變換流站的控制功率后,直流電網(wǎng)最小特征值λcri變大,直流電網(wǎng)的電壓穩(wěn)定性得到改善。且VSC2的功率越小,同時VSC5處的功率越大,λcri增大幅度越明顯,即直流電網(wǎng)的電壓穩(wěn)定性越好。

    改變換流站的控制功率后,直流電網(wǎng)中,直流節(jié)點電壓幅值增加的平均數(shù)量級為10-3;作為對比,交流電網(wǎng)中,交流節(jié)點電壓幅值變化的平均值的數(shù)量級為10-4,可認為交流電網(wǎng)的電壓穩(wěn)定性基本不變。因此,該實驗結(jié)果表明,改變換流站的功率分配,不會影響交流電網(wǎng)的電壓穩(wěn)定性,卻能夠改善直流電網(wǎng)的電壓穩(wěn)定水平,從而提升整個交直流混聯(lián)系統(tǒng)的電壓穩(wěn)定性。

    5 結(jié)論

    本文對傳統(tǒng)的模態(tài)分析法改進,將交直流混聯(lián)系統(tǒng)進行解耦分析,并將模態(tài)分析法應(yīng)用到VSCMTDC交直流混聯(lián)系統(tǒng)的靜態(tài)電壓穩(wěn)定性分析中。同時,將狀態(tài)-模態(tài)指標引入并推廣到交直流系統(tǒng)中,能夠計及有功功率對電壓穩(wěn)定性的影響。所提方法與指標不僅可以量化交、直流負荷變化對電壓穩(wěn)定性的影響,還能夠分析交、直流線路故障情況。

    此外,本文通過改變換流站的控制功率及換流站在交流電網(wǎng)的連接位置,詳細探討了交流電網(wǎng)、VSC換流站及直流電網(wǎng)之間的直接、間接影響。該方法為VSC-MTDC交直流混聯(lián)系統(tǒng)的安全穩(wěn)定運行提供了有價值的參考信息,運行人員可以對交直流系統(tǒng)的關(guān)鍵節(jié)點、線路采取相應(yīng)的補償措施或者調(diào)節(jié)換流站的控制功率,從而改善整個交直流混聯(lián)系統(tǒng)的電壓穩(wěn)定性。

    猜你喜歡
    交直流換流站直流
    基于直流載波通信的LAMOST控制系統(tǒng)設(shè)計
    交直流混合微電網(wǎng)多時間尺度協(xié)同控制
    能源工程(2021年1期)2021-04-13 02:06:12
    直流輸電換流站閥冷卻水系統(tǒng)故障實例統(tǒng)計分析
    一款高效的30V直流開關(guān)電源設(shè)計
    換流站電阻冷卻面板鼓包分析
    湖北電力(2016年8期)2016-05-09 09:03:08
    換流站閥廳避雷器停電例行試驗研究
    湖北電力(2016年8期)2016-05-09 09:03:06
    基于兩次采樣的交直流電功率迭代算法
    電測與儀表(2016年1期)2016-04-12 00:35:08
    非隔離型光伏并網(wǎng)逆變器直流注入抑制方法
    配網(wǎng)20kV變電站交直流一體化電源的應(yīng)用與研究
    基于MATLAB的輕型直流輸電系統(tǒng)的仿真
    亚洲人成网站高清观看| 国产精品乱码一区二三区的特点| 99久国产av精品| 性欧美人与动物交配| 亚洲va在线va天堂va国产| 日本免费一区二区三区高清不卡| 日韩,欧美,国产一区二区三区 | 亚洲一级一片aⅴ在线观看| 精华霜和精华液先用哪个| 看十八女毛片水多多多| 赤兔流量卡办理| 久久人妻av系列| 丝袜喷水一区| 日本撒尿小便嘘嘘汇集6| 亚洲欧美清纯卡通| 日韩欧美免费精品| 久久人人爽人人爽人人片va| 国产精品av视频在线免费观看| 三级国产精品欧美在线观看| av.在线天堂| 亚洲精品在线观看二区| 亚洲精品在线观看二区| 不卡视频在线观看欧美| 网址你懂的国产日韩在线| 你懂的网址亚洲精品在线观看 | 三级男女做爰猛烈吃奶摸视频| 亚洲国产精品国产精品| 精品午夜福利视频在线观看一区| 精品久久久久久久久久免费视频| 欧美激情久久久久久爽电影| 大型黄色视频在线免费观看| 久久久成人免费电影| 俄罗斯特黄特色一大片| 男人的好看免费观看在线视频| 麻豆av噜噜一区二区三区| 久久久国产成人精品二区| 亚洲第一电影网av| 大又大粗又爽又黄少妇毛片口| 在线国产一区二区在线| 国产精品电影一区二区三区| 精品国内亚洲2022精品成人| 午夜激情欧美在线| 我要看日韩黄色一级片| 国产黄色小视频在线观看| 色播亚洲综合网| 国产黄色视频一区二区在线观看 | 深夜精品福利| 偷拍熟女少妇极品色| 少妇裸体淫交视频免费看高清| 亚洲国产色片| 欧美成人免费av一区二区三区| 真实男女啪啪啪动态图| 欧美最黄视频在线播放免费| 啦啦啦韩国在线观看视频| 精品久久久久久久人妻蜜臀av| 99精品在免费线老司机午夜| 久久精品夜夜夜夜夜久久蜜豆| 亚洲综合色惰| 91久久精品电影网| 1000部很黄的大片| 成人特级黄色片久久久久久久| 欧美潮喷喷水| 日本在线视频免费播放| 国产91av在线免费观看| 啦啦啦啦在线视频资源| av女优亚洲男人天堂| 国模一区二区三区四区视频| 久久久久九九精品影院| 亚洲成a人片在线一区二区| 国产aⅴ精品一区二区三区波| 国产免费男女视频| 国产熟女欧美一区二区| 午夜福利18| 国产黄a三级三级三级人| 一级毛片电影观看 | 啦啦啦韩国在线观看视频| 人妻制服诱惑在线中文字幕| 永久网站在线| av在线观看视频网站免费| 久久亚洲精品不卡| 久久久久久伊人网av| 国产亚洲精品久久久久久毛片| 精品一区二区三区视频在线观看免费| 久久人人精品亚洲av| 精品一区二区免费观看| 少妇猛男粗大的猛烈进出视频 | 性插视频无遮挡在线免费观看| 最近在线观看免费完整版| 嫩草影院入口| 婷婷精品国产亚洲av| 黄色配什么色好看| 91久久精品电影网| 国产高清视频在线观看网站| 欧美成人a在线观看| 欧美成人一区二区免费高清观看| 欧美色视频一区免费| 又爽又黄无遮挡网站| 国产aⅴ精品一区二区三区波| 熟女人妻精品中文字幕| 少妇裸体淫交视频免费看高清| 一个人观看的视频www高清免费观看| 国产欧美日韩精品亚洲av| 99久久无色码亚洲精品果冻| 国产中年淑女户外野战色| 久久久久精品国产欧美久久久| 可以在线观看毛片的网站| 少妇高潮的动态图| 国产av在哪里看| 国产一区亚洲一区在线观看| 欧美一区二区国产精品久久精品| 日日摸夜夜添夜夜爱| 热99re8久久精品国产| 极品教师在线视频| 伦精品一区二区三区| 久久久久久九九精品二区国产| 人人妻,人人澡人人爽秒播| 日本 av在线| 美女黄网站色视频| 国产一区二区三区在线臀色熟女| 久久午夜福利片| 久久久久九九精品影院| 亚洲精品456在线播放app| 99久国产av精品国产电影| 午夜福利在线在线| 美女大奶头视频| a级毛片a级免费在线| av视频在线观看入口| 男人和女人高潮做爰伦理| 极品教师在线视频| 在线国产一区二区在线| 三级男女做爰猛烈吃奶摸视频| 91av网一区二区| or卡值多少钱| 最近手机中文字幕大全| 成人性生交大片免费视频hd| 能在线免费观看的黄片| 91久久精品国产一区二区三区| 欧美激情在线99| 卡戴珊不雅视频在线播放| 99热这里只有精品一区| 亚洲国产精品合色在线| 97人妻精品一区二区三区麻豆| 一级毛片久久久久久久久女| 性欧美人与动物交配| 亚洲精华国产精华液的使用体验 | 亚洲性久久影院| 我的老师免费观看完整版| 免费电影在线观看免费观看| 国产成人影院久久av| 色av中文字幕| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区免费观看 | av视频在线观看入口| 久久精品人妻少妇| 热99在线观看视频| 日韩av在线大香蕉| a级毛色黄片| 中文在线观看免费www的网站| 少妇熟女aⅴ在线视频| 日日啪夜夜撸| 国产精品一及| or卡值多少钱| 国产麻豆成人av免费视频| 免费人成在线观看视频色| 欧美日韩一区二区视频在线观看视频在线 | 欧美一区二区国产精品久久精品| 听说在线观看完整版免费高清| 夜夜爽天天搞| 麻豆精品久久久久久蜜桃| 可以在线观看毛片的网站| 欧美区成人在线视频| 此物有八面人人有两片| 久久精品国产自在天天线| 国产成人freesex在线 | 欧美又色又爽又黄视频| 日韩三级伦理在线观看| 国产精品美女特级片免费视频播放器| 国产精品电影一区二区三区| 你懂的网址亚洲精品在线观看 | 亚洲av一区综合| 日本黄大片高清| av中文乱码字幕在线| 淫秽高清视频在线观看| 可以在线观看的亚洲视频| 日本在线视频免费播放| 欧美最新免费一区二区三区| 成人国产麻豆网| 欧美日韩乱码在线| 能在线免费观看的黄片| 看片在线看免费视频| 99国产极品粉嫩在线观看| 久久精品影院6| 午夜视频国产福利| 高清毛片免费看| 中国美女看黄片| 亚洲国产精品国产精品| 成人永久免费在线观看视频| 午夜福利高清视频| 国产精品一及| 日日摸夜夜添夜夜爱| 国产中年淑女户外野战色| 韩国av在线不卡| 九九热线精品视视频播放| 国产精品久久久久久av不卡| 亚洲精品456在线播放app| 日韩欧美在线乱码| 国产精品久久久久久久电影| 老熟妇乱子伦视频在线观看| 12—13女人毛片做爰片一| 国产高清视频在线播放一区| 十八禁网站免费在线| 男人和女人高潮做爰伦理| av黄色大香蕉| av天堂中文字幕网| 97超视频在线观看视频| 国产精品久久久久久av不卡| 欧洲精品卡2卡3卡4卡5卡区| 日韩,欧美,国产一区二区三区 | 精品一区二区三区视频在线| 亚洲国产精品久久男人天堂| 精品久久久久久成人av| 女人十人毛片免费观看3o分钟| 午夜老司机福利剧场| 欧美极品一区二区三区四区| 99久久精品一区二区三区| 日本黄大片高清| 91麻豆精品激情在线观看国产| 日日干狠狠操夜夜爽| 亚洲av熟女| 国产精品免费一区二区三区在线| 69av精品久久久久久| 久久精品国产亚洲网站| 18+在线观看网站| 又粗又爽又猛毛片免费看| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 成人鲁丝片一二三区免费| 在线免费观看不下载黄p国产| 欧美zozozo另类| 成年免费大片在线观看| av在线亚洲专区| 免费看a级黄色片| 九九热线精品视视频播放| 国产成人aa在线观看| 国产成人a区在线观看| 天天一区二区日本电影三级| 久久精品国产99精品国产亚洲性色| 日韩 亚洲 欧美在线| 成人永久免费在线观看视频| 日韩欧美一区二区三区在线观看| 99久久精品国产国产毛片| a级一级毛片免费在线观看| 欧美日韩在线观看h| 听说在线观看完整版免费高清| 九九爱精品视频在线观看| 午夜福利在线观看免费完整高清在 | 午夜久久久久精精品| 黄色一级大片看看| 国产一区二区三区av在线 | 九九热线精品视视频播放| 日韩亚洲欧美综合| 国产在线精品亚洲第一网站| 99久久九九国产精品国产免费| 国产精品亚洲一级av第二区| 精品人妻偷拍中文字幕| 九九热线精品视视频播放| 成人无遮挡网站| 99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| 欧美日本视频| 99久久久亚洲精品蜜臀av| 国产高清视频在线观看网站| 91在线精品国自产拍蜜月| 99热这里只有是精品在线观看| 日日摸夜夜添夜夜添av毛片| 男人舔奶头视频| 我要看日韩黄色一级片| 国产久久久一区二区三区| 在线观看免费视频日本深夜| 美女黄网站色视频| 亚洲成人久久爱视频| 男女那种视频在线观看| 秋霞在线观看毛片| 成年av动漫网址| 亚洲第一区二区三区不卡| 波多野结衣高清作品| 国内精品宾馆在线| 国产亚洲精品综合一区在线观看| 久久亚洲精品不卡| 国产麻豆成人av免费视频| 精品福利观看| 99热只有精品国产| 亚洲国产色片| 午夜视频国产福利| 国产午夜精品久久久久久一区二区三区 | 国产三级在线视频| 欧美一区二区亚洲| 老司机午夜福利在线观看视频| av在线蜜桃| 深爱激情五月婷婷| 内地一区二区视频在线| 九九在线视频观看精品| 亚洲精品色激情综合| 少妇的逼好多水| 最好的美女福利视频网| 夜夜爽天天搞| 国内精品久久久久精免费| 午夜福利高清视频| 精品不卡国产一区二区三区| 亚洲欧美成人综合另类久久久 | 成人精品一区二区免费| 三级男女做爰猛烈吃奶摸视频| 国产国拍精品亚洲av在线观看| 国产欧美日韩一区二区精品| 亚洲三级黄色毛片| 中文字幕av成人在线电影| 激情 狠狠 欧美| 一卡2卡三卡四卡精品乱码亚洲| 你懂的网址亚洲精品在线观看 | 久久久a久久爽久久v久久| 精品熟女少妇av免费看| 亚洲自偷自拍三级| 少妇熟女欧美另类| 国产毛片a区久久久久| av在线蜜桃| 欧美激情久久久久久爽电影| 久久精品久久久久久噜噜老黄 | 午夜精品一区二区三区免费看| 亚洲中文字幕一区二区三区有码在线看| 精品久久国产蜜桃| 国产精品一区www在线观看| 日韩成人伦理影院| 深夜a级毛片| 我要看日韩黄色一级片| 亚洲精品色激情综合| 99久久无色码亚洲精品果冻| 十八禁网站免费在线| 能在线免费观看的黄片| 两性午夜刺激爽爽歪歪视频在线观看| 深夜a级毛片| 一区二区三区四区激情视频 | 午夜爱爱视频在线播放| 18+在线观看网站| 精品福利观看| 看黄色毛片网站| 成人av一区二区三区在线看| 亚洲最大成人av| 免费黄网站久久成人精品| 亚洲成人av在线免费| 一级毛片电影观看 | 亚洲欧美日韩东京热| 成人亚洲欧美一区二区av| 国产成人a区在线观看| 99久久久亚洲精品蜜臀av| 成人漫画全彩无遮挡| 欧美色欧美亚洲另类二区| 国内精品一区二区在线观看| 亚洲精品色激情综合| 女人被狂操c到高潮| 精品不卡国产一区二区三区| 久久精品影院6| 午夜福利成人在线免费观看| 少妇熟女aⅴ在线视频| 精品日产1卡2卡| 国产午夜福利久久久久久| 亚洲无线在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品1区2区在线观看.| 亚洲精品国产成人久久av| 中文字幕免费在线视频6| 两个人的视频大全免费| 性插视频无遮挡在线免费观看| 麻豆一二三区av精品| 最好的美女福利视频网| 夜夜爽天天搞| 中文字幕av在线有码专区| 日韩三级伦理在线观看| 亚洲欧美精品综合久久99| 嫩草影视91久久| av视频在线观看入口| 国产又黄又爽又无遮挡在线| 黄色日韩在线| 一区福利在线观看| 最近手机中文字幕大全| 看黄色毛片网站| 少妇丰满av| 日本欧美国产在线视频| 成人永久免费在线观看视频| 丰满人妻一区二区三区视频av| aaaaa片日本免费| 99视频精品全部免费 在线| АⅤ资源中文在线天堂| 干丝袜人妻中文字幕| 色哟哟哟哟哟哟| 国产高清有码在线观看视频| 精品国内亚洲2022精品成人| 国产 一区精品| 久久久国产成人免费| 婷婷亚洲欧美| 欧美精品国产亚洲| 免费看日本二区| 国内精品美女久久久久久| 国产精品爽爽va在线观看网站| 99热这里只有是精品在线观看| 精品福利观看| 一区福利在线观看| 欧美成人一区二区免费高清观看| 国产欧美日韩精品一区二区| 99久久中文字幕三级久久日本| 夜夜夜夜夜久久久久| 国产女主播在线喷水免费视频网站 | 2021天堂中文幕一二区在线观| 午夜免费激情av| 国产探花极品一区二区| 我的女老师完整版在线观看| 亚洲欧美精品综合久久99| 精品一区二区三区人妻视频| 嫩草影院新地址| 亚洲人成网站高清观看| 午夜福利18| 欧美色视频一区免费| 久久亚洲国产成人精品v| 国产精品人妻久久久影院| 欧美中文日本在线观看视频| 亚洲欧美精品自产自拍| 狠狠狠狠99中文字幕| 免费一级毛片在线播放高清视频| 国产成人91sexporn| 亚洲精品456在线播放app| 国产三级中文精品| 国产淫片久久久久久久久| 亚洲最大成人中文| 日日啪夜夜撸| 国产三级中文精品| 久久天躁狠狠躁夜夜2o2o| 国产精品人妻久久久久久| 亚洲av免费在线观看| 国产伦精品一区二区三区四那| 亚洲丝袜综合中文字幕| 美女黄网站色视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人精品中文字幕电影| 麻豆成人午夜福利视频| 国产精品久久久久久久电影| 欧美激情在线99| 国产伦在线观看视频一区| 成人特级av手机在线观看| 99热这里只有是精品50| 亚洲av.av天堂| 亚洲国产精品久久男人天堂| 婷婷色综合大香蕉| 欧美丝袜亚洲另类| 一区二区三区高清视频在线| 乱系列少妇在线播放| 蜜臀久久99精品久久宅男| 老司机福利观看| 久久久久久久久大av| 99精品在免费线老司机午夜| 午夜亚洲福利在线播放| .国产精品久久| 国产片特级美女逼逼视频| 男人舔奶头视频| 亚洲欧美精品综合久久99| 高清日韩中文字幕在线| 赤兔流量卡办理| 亚洲av一区综合| 亚洲色图av天堂| 一本一本综合久久| 亚洲真实伦在线观看| 精品久久久久久成人av| 国产av不卡久久| 亚洲性久久影院| 18禁在线无遮挡免费观看视频 | 校园人妻丝袜中文字幕| 亚洲欧美日韩东京热| 在线观看午夜福利视频| 日韩人妻高清精品专区| 欧美中文日本在线观看视频| 国产精品免费一区二区三区在线| 亚洲av中文字字幕乱码综合| 99国产精品一区二区蜜桃av| 我的老师免费观看完整版| 成年免费大片在线观看| 日韩制服骚丝袜av| 免费一级毛片在线播放高清视频| 亚洲av熟女| 国产极品精品免费视频能看的| aaaaa片日本免费| 国产精品av视频在线免费观看| 亚洲欧美日韩高清在线视频| 超碰av人人做人人爽久久| 一级黄色大片毛片| 中文字幕精品亚洲无线码一区| 最新中文字幕久久久久| 可以在线观看毛片的网站| 精品免费久久久久久久清纯| 永久网站在线| 欧美中文日本在线观看视频| 日韩三级伦理在线观看| 伦精品一区二区三区| 最后的刺客免费高清国语| 22中文网久久字幕| 亚洲人成网站在线播| АⅤ资源中文在线天堂| 看免费成人av毛片| 日本一二三区视频观看| 日本五十路高清| 国产av一区在线观看免费| 久久久精品94久久精品| 亚洲国产日韩欧美精品在线观看| 超碰av人人做人人爽久久| 大香蕉久久网| 国产精品久久视频播放| 在线a可以看的网站| 男人的好看免费观看在线视频| 精品国内亚洲2022精品成人| 亚洲精品在线观看二区| 免费看日本二区| 亚洲欧美成人综合另类久久久 | 国产欧美日韩精品亚洲av| 真实男女啪啪啪动态图| 在线免费观看的www视频| 人妻久久中文字幕网| 此物有八面人人有两片| 国产美女午夜福利| 日本免费一区二区三区高清不卡| 午夜爱爱视频在线播放| 国内精品美女久久久久久| 久久九九热精品免费| 亚洲一区二区三区色噜噜| 亚洲欧美日韩高清专用| 国产高清视频在线播放一区| 一进一出抽搐动态| 人人妻,人人澡人人爽秒播| 久久99热6这里只有精品| 国产真实伦视频高清在线观看| 老司机午夜福利在线观看视频| 亚洲精品成人久久久久久| 午夜福利成人在线免费观看| 精华霜和精华液先用哪个| 特大巨黑吊av在线直播| 国产一区二区在线观看日韩| 成熟少妇高潮喷水视频| 无遮挡黄片免费观看| 如何舔出高潮| 男女做爰动态图高潮gif福利片| 老师上课跳d突然被开到最大视频| 国产高清视频在线观看网站| 99久久精品一区二区三区| 我要看日韩黄色一级片| 日韩,欧美,国产一区二区三区 | 亚洲欧美日韩无卡精品| 亚洲久久久久久中文字幕| 草草在线视频免费看| 久久久久久久亚洲中文字幕| 免费电影在线观看免费观看| 亚洲av免费在线观看| 免费看日本二区| 亚洲乱码一区二区免费版| 最新在线观看一区二区三区| 99国产极品粉嫩在线观看| 久久草成人影院| 亚洲av.av天堂| 少妇的逼水好多| АⅤ资源中文在线天堂| 亚洲欧美成人精品一区二区| 免费人成视频x8x8入口观看| 国产一级毛片七仙女欲春2| 日本撒尿小便嘘嘘汇集6| 五月伊人婷婷丁香| 国产在线精品亚洲第一网站| 欧美性猛交黑人性爽| 欧美绝顶高潮抽搐喷水| 毛片一级片免费看久久久久| 国语自产精品视频在线第100页| 亚洲欧美日韩高清专用| 波多野结衣高清无吗| 午夜福利成人在线免费观看| 在线免费观看的www视频| 看黄色毛片网站| 91在线精品国自产拍蜜月| av天堂中文字幕网| 亚洲最大成人手机在线| 日韩人妻高清精品专区| 51国产日韩欧美| 国产精品福利在线免费观看| 内地一区二区视频在线| 色吧在线观看| 亚洲丝袜综合中文字幕| 成人高潮视频无遮挡免费网站| 欧美日韩在线观看h| 精品日产1卡2卡| 丰满人妻一区二区三区视频av| 国产欧美日韩精品一区二区| 日本a在线网址| 成年女人看的毛片在线观看| 精品无人区乱码1区二区| 日本 av在线| 精品国内亚洲2022精品成人| 超碰av人人做人人爽久久| 欧美区成人在线视频| 少妇人妻精品综合一区二区 | 久久久a久久爽久久v久久| 长腿黑丝高跟| 国产激情偷乱视频一区二区| 我要搜黄色片| 午夜福利成人在线免费观看| 九九久久精品国产亚洲av麻豆| 久久久久久久亚洲中文字幕| 中国国产av一级| 99久久九九国产精品国产免费| 国产男人的电影天堂91| 一个人观看的视频www高清免费观看| 人人妻人人澡人人爽人人夜夜 | 嫩草影院精品99|