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

    基于暫態(tài)數(shù)據(jù)戴維南等值參數(shù)辨識的時變廣義短路比計算方法

    2024-03-05 06:53:30孫銳鵬王奕鑫王玉昊吳浩辛煥海但揚清
    能源工程 2024年1期
    關(guān)鍵詞:等值廣義端口

    孫銳鵬,王奕鑫,王玉昊,吳浩,辛煥海,但揚清

    (1.浙江大學(xué)電氣工程學(xué)院,浙江 杭州 310027;2.浙江省電力有限公司經(jīng)濟研究院,浙江 杭州 310008)

    0 引 言

    隨著我國西北風(fēng)電、光伏等新能源基地相繼建成,各區(qū)域電網(wǎng)之間大規(guī)模、長距離輸電形成多饋入直流系統(tǒng)[1]。多條直流饋入同一交流電網(wǎng)一方面消納西北區(qū)域多余的風(fēng)電光伏新能源電力,另一方面緩解了東南區(qū)域負荷需求大的現(xiàn)狀,實現(xiàn)了電網(wǎng)的電力供需平衡,加快了電力系統(tǒng)向“高比例可再生能源”和“高比例電力電子設(shè)備”的新型電力系統(tǒng)轉(zhuǎn)變的速度[2,3]。多條直流饋入同一交流電網(wǎng)加重了多條直流之間的相互耦合形勢,產(chǎn)生交直流系統(tǒng)之間相互影響的問題,是影響電網(wǎng)安全穩(wěn)定水平的因素之一。同時越多直流饋入同一交流電網(wǎng),整個系統(tǒng)存在電壓失穩(wěn)風(fēng)險,連續(xù)換相失敗隱患,甚至導(dǎo)致系統(tǒng)無法安全穩(wěn)定運行。

    目前,單饋入和多饋入直流系統(tǒng)能否安全穩(wěn)定運行一般由受端交流系統(tǒng)電網(wǎng)強度衡量。單饋入直流系統(tǒng)的受端交流電網(wǎng)強度由短路比表征,其定義為交流系統(tǒng)短路容量與直流額定容量之比。為了拓展到多饋入直流系統(tǒng)電網(wǎng)強度量化,國際大電網(wǎng)會議(International Council on Large Electric Systems,CIGRE)提出了多饋入短路比指標(multi infeed short circuit ratio,MISCR)、多饋入有效短路比(multi-infeed effective short circuit ratio,MIESCR)、多饋入相互作用影響因子指標(multi-infeed interaction factor,MIIF)[5]。雖然解決了描述受端交流電網(wǎng)強度評價指標缺失問題,但文獻[6]指出MISCR、MIESCR與MIIF物理意義不明確,邊界條件弱,并根據(jù)電壓穩(wěn)定分岔與系統(tǒng)雅可比矩陣奇異的等價性,建立交直流系統(tǒng)雅可比矩陣,定義了廣義短路比指標(generalized short circuit ratio,gSCR),克服了傳統(tǒng)短路比機理不明確的缺陷。

    gSCR一般通過等值獲得直流饋入節(jié)點的等值導(dǎo)納矩陣計算。等值導(dǎo)納矩陣的計算是電力系統(tǒng)阻抗辨識問題,基于電網(wǎng)的信息來源,該問題的解決方法一般分為兩大類:基于網(wǎng)絡(luò)參數(shù)與拓撲的白箱法和基于機電暫態(tài)仿真端口量測數(shù)據(jù)的黑箱法[7]。

    第一類白箱法是基于已知完整的外網(wǎng)結(jié)構(gòu)、線路元件參數(shù)和潮流信息,通過建立網(wǎng)絡(luò)方程計算等值系統(tǒng)的戴維南等值參數(shù)。如文獻[8]提出基于節(jié)點導(dǎo)納矩陣的戴維南辨識方法,該方法基于全網(wǎng)信息,得到系統(tǒng)節(jié)點導(dǎo)納矩陣,并對導(dǎo)納矩陣進行修正、求逆求解戴維南等值參數(shù),但該方法在電網(wǎng)狀態(tài)改變時,無法準確求解[9-10]。文獻[11]利用仿真分析了基于BPA-SCCP軟件的多端口等值功能,闡明SCCP等值方法,但該方法不適用于系統(tǒng)參數(shù)時變的情況。

    對于白箱法,其自身存在一定的局限性,一方面實際系統(tǒng)內(nèi)部參數(shù)并非可以做到全知;另一方面系統(tǒng)的某些參數(shù)可能在不斷的變化,而又難以精確描述。

    第二類方法黑箱法是基于機電暫態(tài)仿真,獲得等值端口機電仿真數(shù)據(jù),提取機電仿真數(shù)據(jù)之中隱含的系統(tǒng)特征,經(jīng)過數(shù)據(jù)處理等手段反推戴維南等值參數(shù)。最簡單的黑箱法僅利用兩個時間斷面的仿真數(shù)據(jù)作差,可直接解出戴維南等值阻抗,故稱兩點法。此外,文獻[12]提出一種基于時域仿真的戴維南等值參數(shù)辨識方法,其采用補償法,結(jié)合全局網(wǎng)絡(luò)方程計算等值節(jié)點處的開路電壓和短路電流,辨識戴維南等值阻抗。但是該方法在形成全局網(wǎng)絡(luò)方程時,將負荷統(tǒng)一看作恒阻抗負荷,忽略了系統(tǒng)非線性負荷的動態(tài)特性。

    對于黑箱法,其在辨識數(shù)據(jù)上不可避免地會造成戴維南參數(shù)零點漂移問題[13];在方法上,不論是兩點法還是時域仿真并補償?shù)姆椒?,對于負荷都失去了有效性,原因是絕大多數(shù)負荷的阻抗模型是非線性的,節(jié)點負荷阻抗是變化的,其值與該時刻下負荷特性曲線的切線大小相關(guān)[14]。這表明將負荷統(tǒng)一為恒阻抗負荷無法計及非線性負荷對戴維南等值參數(shù)的影響。

    為此,本文在機電暫態(tài)時間尺度下,提出基于戴維南多端口等值模型,采用最小二乘法辨識系統(tǒng)故障后暫態(tài)數(shù)據(jù),得到受端交流電網(wǎng)的時變等值導(dǎo)納矩陣,進而計算gSCR。通過單饋入直流仿真算例和某省級電網(wǎng)算例,本文分析了負荷對戴維南等值參數(shù)的影響,在減少了計算量的同時考慮了系統(tǒng)中各元件的動態(tài)過程,實現(xiàn)了對gSCR的進一步分析,精準量化了受端交流電網(wǎng)的強度,驗證了所提方法的實用性和有效性。

    1 基于戴維南等值的廣義短路比計算

    1.1 多饋入直流系統(tǒng)的戴維南等值模型

    相較于單饋入直流系統(tǒng),多饋入直流系統(tǒng)更為復(fù)雜,多條直流之間相互耦合同時直流系統(tǒng)與交流電網(wǎng)相互影響。因此,對多饋入直流系統(tǒng)在饋入節(jié)點做戴維南等值時,在盡可能簡化其受端交流電網(wǎng)結(jié)構(gòu)的同時突出直流間的耦合關(guān)系。多饋入戴維南等值結(jié)構(gòu)將多饋入直流系統(tǒng)分為多饋入直流(multi-infeed direct current,MIDC)和受端交流電網(wǎng)。其中受端交流電網(wǎng)為研究區(qū)域,等值前多饋入直流系統(tǒng)如圖1所示,等值后多饋入系統(tǒng)戴維南等值結(jié)構(gòu)如圖2所示。

    圖1 多饋入直流系統(tǒng)

    圖2 多饋入系統(tǒng)戴維南等值結(jié)構(gòu)

    圖2中研究區(qū)域等值端口數(shù)量n等于多饋入直流落點,Ui為端口等值節(jié)點的電壓相量,Si為端口等值節(jié)點i注入研究區(qū)域的視在功率,其值等于研究區(qū)域內(nèi)與端口等值節(jié)點i相連線路上傳輸功率之和。在研究區(qū)域每個端口節(jié)點i處,將研究區(qū)域等值為恒定電勢Ei、Paci+JQaci、導(dǎo)納Yii表示虛擬發(fā)電機的電壓、功率和等值端口節(jié)點與系統(tǒng)各發(fā)電機之間的等值聯(lián)絡(luò)導(dǎo)納。在研究區(qū)域端口等值節(jié)點i與端口等值節(jié)點j之間,將研究區(qū)域等值為導(dǎo)納Yij=Gij+jBij,表示直流i和直流j的耦合關(guān)系。

    1.2 廣義短路比計算

    文獻[3]基于以下3個假設(shè)條件,提出多饋入系統(tǒng)廣義短路比指標。

    假設(shè)1:考慮的n條直流系統(tǒng)都是“相似的”,且控制方式都為整流側(cè)定功率控制、逆變側(cè)定熄弧角控制。即除了額定容量外,控制方式、以自身容量為基準的主電路參數(shù)標幺值都相同。

    假設(shè)2:受端交流電網(wǎng)是電感型、拓撲連通的系統(tǒng)。

    假設(shè)3:在額定點處直流系統(tǒng)注入交流系統(tǒng)的無功功率很小,即直流間直流落點間交流聯(lián)絡(luò)功率遠小于其傳輸極限。

    多饋入直流系統(tǒng)廣義短路比定義為拓展導(dǎo)納矩陣的最小特征值,即:

    JB為研究區(qū)域的拓展導(dǎo)納矩陣,具體為:

    式中:Si為n饋入直流系統(tǒng)中第i條直流的額定容量?;诙囵伻氪骶S南等值模型,可計算多饋入等值系統(tǒng)的節(jié)點阻抗矩陣,對節(jié)點阻抗矩陣求逆為多饋入等值系統(tǒng)的節(jié)點導(dǎo)納矩陣,由式(2)可計算gSCR。特別指出,對于單饋入直流系統(tǒng),其廣義短路比為拓展導(dǎo)納矩陣維數(shù)為一的情況[15],并基于式(2)可知,決定廣義短路比的是系統(tǒng)等值阻抗。

    gSCR描述多饋入交直流系統(tǒng)強度和穩(wěn)定性時,須借助臨界短路比(critical generalized short circuit ratio,CgSCR)輔助評判,特別指出一般CgSCR≈2,gSCR略大于2為弱交流系統(tǒng),小于2為極弱交流系統(tǒng)[16,17]。特別指出,廣義短路比計算基于拓展導(dǎo)納矩陣的最小特征值,系統(tǒng)的導(dǎo)納大小不能直接決定廣義短路比的大小,但依據(jù)經(jīng)驗,導(dǎo)納矩陣各元素越小,廣義短路比越小。

    2 基于暫態(tài)數(shù)據(jù)的時變戴維南等值參數(shù)辨識

    2.1 基于暫態(tài)數(shù)據(jù)的多饋入時變戴維南等值參數(shù)辨識方法

    戴維南定理表明:在任意時刻,交流電網(wǎng)可看作由等值電勢通過等值阻抗進行供電的戴維南等值電路。當拓展到多端口時,n饋入交直流系統(tǒng)可看作n個等值電壓源通過n(n+1)/2個等值阻抗進行供電的互聯(lián)網(wǎng)絡(luò)[18]。

    如圖2所示,t時刻注入研究區(qū)域端口節(jié)點多端口戴維南等值方程如下:

    式中:上標t表示第t時刻數(shù)據(jù),]T表示端口流入電流n維列向量,表示端口節(jié)點電壓的n維列向量,Y表示多端口等值系統(tǒng)n×n的節(jié)點導(dǎo)納矩陣,特別指出節(jié)點阻抗矩陣為Zth=Y-1,表示研究區(qū)域的等值電勢,Y'為n×n的對角矩陣表示端口節(jié)點和等值電勢之間的自導(dǎo)納矩陣diag(y11,y22,K,ynn),其中對角元素分別表示等值節(jié)點的對地導(dǎo)納。

    根據(jù)式(3)在t時刻時,取式(3)的第i行,在端口節(jié)點i處有:

    同時在仿真數(shù)據(jù)中,在端口節(jié)點i處注入:

    根據(jù)式(4)和式(5),多饋入戴維南等值參數(shù)最小二乘估計模型為:

    式中:m為在第t時刻的窗口大小,而端口節(jié)點個數(shù)為n時,可用的方程式共nk個。同時,式(3)中待求解等值參數(shù)Ei有n個,節(jié)點導(dǎo)納矩陣Y是對稱陣,有n(n+1)/2個參數(shù),因此待求等值參數(shù)共n(n+3)/2個。隨著多饋入直流系統(tǒng)中直流個數(shù)的增加,待求等值參數(shù)呈平方關(guān)系增加。若使得式(6)有解,則必須nm>n(n+3)/2成立,即m≥[(n+3)/2]。

    將基于式(6)可辨識出在t時刻的戴維南等值導(dǎo)納矩陣與結(jié)合n饋入直流系統(tǒng)中直流的額定容量,由式(2)計算出t時刻的廣義短路比。

    2.2 基于暫態(tài)數(shù)據(jù)的多饋入時變戴維南等值參數(shù)辨識計算步驟

    多端口時變戴維南等值方法具體流程如圖所示;

    基于暫態(tài)數(shù)據(jù)多饋入戴維南等值參數(shù)計算方法的主要步驟如下所述:

    1)對多饋入直流系統(tǒng)暫態(tài)仿真,選取故障恢復(fù)后等值端口的n個節(jié)點的電氣量數(shù)據(jù)(S,U)。

    2)計算端口節(jié)點的節(jié)點電壓以及注入電流,并選取每時刻的窗口大小m。每時刻依據(jù)步驟1)的電氣量數(shù)據(jù),可得到n個端口節(jié)點的電壓相量和n個端口節(jié)點的視在功率。計算在各個端口的電壓電流,將后者除以該節(jié)點的電壓并取共軛便得到流入n個端口節(jié)點處的電流相量,并由端口節(jié)點數(shù)n計算該時刻戴維南值參數(shù)的窗口大小m,將n×m個節(jié)點電壓與注入節(jié)點電流數(shù)據(jù)作為步驟3)的辨識數(shù)據(jù)。

    3)辨識該時刻的多端口戴維南等值參數(shù)。將步驟2)中計算n×m個端口節(jié)點電壓和注入端口節(jié)點電流數(shù)據(jù),建立如式(3)所示的方程,并基于此建立如式(6)所示的多端口戴維南等值參數(shù)的最小二乘模型,獲得戴維南等值參數(shù),形成受端交流電網(wǎng)的等值節(jié)點矩陣。

    4)依據(jù)式(1)計算該時刻的廣義短路比。

    5)將窗口隨時間移動,重復(fù)步驟1至步驟4,直到系統(tǒng)穩(wěn)定結(jié)束,獲得時變廣義短路比。

    3 電力系統(tǒng)負荷對戴維南等值參數(shù)的影響

    3.1 戴維南等值參數(shù)理論值計算

    系統(tǒng)等值阻抗可以根據(jù)網(wǎng)絡(luò)阻抗以及負荷等值阻抗的串并聯(lián)關(guān)系聚合得到,該聚合阻抗即為從等值端口看進去的理論戴維南等值阻抗。因此從外特性等價的角度來看,復(fù)雜系統(tǒng)總可以簡化成圖2中所示的戴維南等值模型[19]。

    戴維南等值參數(shù)辨識的目標就是系統(tǒng)理論等值阻抗Zth,其可以作為驗證辨識戴維南等值阻抗計算準確的依據(jù)。由于網(wǎng)絡(luò)阻抗在暫態(tài)過程中保持恒定,后續(xù)需關(guān)注負荷等值阻抗對于戴維南等值參數(shù)阻抗Zth的影響。

    從物理意義上講,電力系統(tǒng)負荷是消耗功率,但不發(fā)出功率,因此在戴維南理論值計算中,將其按時變的等值阻抗。負荷節(jié)點i處的等值阻抗可結(jié)合t時刻與t?1時刻的負荷節(jié)點電壓與注入電流求出:

    3.2 靜態(tài)負荷模型對于戴維南等值參數(shù)的影響

    靜態(tài)負荷由恒阻抗、恒電流和恒功率負荷組成,其中恒阻抗負荷為線性負荷;恒電流和恒功率負荷為非線性負荷模型。式(8)為靜態(tài)負荷的數(shù)學(xué)模型:

    式中:KPZ、KPI、KPP分別為有功負荷中,恒阻抗負荷、恒電流負荷、恒功率負荷所占的比例,且滿足KPZ+KPI+KPP=1,KQZ、KQI、KQP類似。PL0為負荷有功分量初始值,QL0為負荷無功分量初始值,U0為節(jié)點電壓幅值的初始值;PL為負荷有功分量實際值;QL為負荷無功分量實際值;U為負荷節(jié)點電壓幅值的實際值。其計算值:

    其中ux、uy表示負荷節(jié)點電壓的實部和虛部。已知負荷節(jié)點電壓U,可知負荷功率:

    式中:ix、iy表示注入負荷節(jié)點電流的實部和虛部,可解得電流:

    對式(11)線性化:

    其中消去中間變量可得:

    式中:Gxx為在該時刻線性化后電壓實部對于電流實部的系數(shù),類似的Bxy、?Bxy、Gyy為電壓虛實部對于電流虛實部系數(shù),其值如下:

    引入如下轉(zhuǎn)移矩陣:

    利用上述轉(zhuǎn)移矩陣T對式(13)進行變換,可得:

    其中Ypn各個元素為:

    故靜態(tài)負荷線性化后的戴維南阻抗為:

    將不同的KPZ、KPI、KPP,與KQZ、KQI、KQP帶入式(18)。當且僅當KPZ=KQZ=1,即ZIP負荷為恒阻抗負荷時,負荷戴維南等值阻抗恒定不變,由式(18)可得其值為:

    此外在ZIP負荷系數(shù)取其他值時,由于z12≠0,因此,非線性負荷的戴維南阻抗包含電流向量共軛的耦合項,與各時刻電壓電流相關(guān)。如在電流變化較大的時刻,傳統(tǒng)認為非線性負荷是恒阻抗的做法對于戴維南等值參數(shù)具有較大的誤差。

    4 算例仿真分析

    4.1 單饋入構(gòu)造算例

    基于PSD-BPA暫態(tài)仿真軟件搭建單饋入直流系統(tǒng)模型,以驗證圖3中提出的基于戴維南參數(shù)辨識的時變廣義短路比計算方法的有效性,并分析負荷特性對時變廣義短路比的影響。

    圖3 多饋入直流系統(tǒng)戴維南參數(shù)辨識流程圖

    構(gòu)造單饋入直流系統(tǒng)拓撲結(jié)構(gòu)如圖4所示,在節(jié)點1處對受端交流電網(wǎng)進行戴維南等值。電網(wǎng)中各參數(shù)如下:受端交流電網(wǎng)電壓等級為500kV,基準功率為1000MVA,線路阻抗ZL1=ZL2=0.02+j0.06,負荷有功功率PL、無功功率QL分別為2000MW、1000MVar,其中靜態(tài)負荷功率占比42%,感應(yīng)電機占比58%,其中KPZ、KPI、KPP、KQZ、KQI、KQP系數(shù)分別0.13、0.34、0.53、0.13、0.34、0.53,E節(jié)點接入無窮大電源。在節(jié)點1處設(shè)置三相金屬接地故障,故障持續(xù) 0.1s以故障切除作為辨識數(shù)據(jù)的起始時刻0時刻,每個窗口大小m為2,包含2個時刻的數(shù)據(jù),窗口起始時刻至結(jié)束時刻時間間隔為0.01s,每個窗口計算出的廣義短路比為窗口結(jié)束時刻的廣義短路比。故障前系統(tǒng)戴維南理論阻抗為0.0489+j0.1170,廣義短路比為2.4628,暫態(tài)數(shù)據(jù)中各節(jié)點電壓與節(jié)點注入功率精度為10-5。

    圖4 單饋入直流系統(tǒng)等值算例

    將綜合負荷模型接入節(jié)點2,驗證所提方法的準確性。節(jié)點1電壓及注入功率如圖5所示,故障清除后節(jié)點1的節(jié)點電壓與注入功率在0.3s后基本保持不變。因此,選取節(jié)點1在故障清除后0~0.3s期間的U、S等電氣量作為辨識數(shù)據(jù),采用所提方法辨識系統(tǒng)等值阻抗,辨識與理論等值阻抗對比如圖6所示,時變廣義短路比對比如圖7所示。

    圖5 節(jié)點1電壓及注入功率

    圖6 辨識與理論等值阻抗對比

    圖7 時變廣義短路比對比

    由圖5、圖6和圖7可知:

    1)在圖中,辨識與理論電阻、電抗以及廣義短路比曲線基本重合,這是基于所提方法辨識得到時變電阻、電抗以及廣義短路比,在0.3s范圍內(nèi)均能夠較好擬合戴維南理論電阻、電抗和廣義短路比;

    2)在圖中,廣義短路比的變化趨勢與等值電抗的變化趨勢相反,這與單饋入直流算例中基于式(1)得出廣義短路比與阻抗成反比關(guān)系結(jié)論相吻合;

    3)在0.01~0.05s內(nèi),共有五個窗口,計算得到的每個窗口理論廣義短路比從1.6541減小至0.5576后,迅速增大到7.4124,最后減小到2.9885。這是在期間內(nèi)由于節(jié)點2處的感應(yīng)電機負荷在故障恢復(fù)時,其有功功率從1543MW近似線性增加到1730MW,無功功率從1840MVar近似線性減小至1292MVar,同時負荷節(jié)點電壓從0.9701增加到1.0325后,減小至1.0215。其中最大電壓差出現(xiàn)在0.02~0.03s,為0.0216,使得由式(7)基于0.02s窗口內(nèi)的電壓、功率數(shù)據(jù)計算的負荷理論阻抗大,導(dǎo)致節(jié)點1處等值阻抗大,基于式(2)得到的理論廣義短路比最小,為0.5576;同理整個期間0~0.04s功率線性變化,而0.03~0.04s的電壓差最小僅為0.0011,使得在0.03s處得到的理論廣義短路比最大,為7.4124;

    4)當0.05s以后,系統(tǒng)理論阻抗變化不超過2%,其值與故障前戴維南理論阻抗0.0489+j*0.1170不超過變化2%,理論廣義短路比變化不超過5%。由于感應(yīng)電機暫態(tài)過程緩慢逐漸恢復(fù)至穩(wěn)態(tài),其中注入感應(yīng)電機負荷的功率在1870MW、1220MVar附近,范圍不超過3%,且負荷節(jié)點電壓變化在0.1以內(nèi),負荷等值阻抗變化不超過2%,其值與故障等值前負荷阻抗接近,與恒定的線路阻抗聚合后,計算所得系統(tǒng)廣義短路比接近故障前的廣義短路比。

    為了量化辨識戴維南阻抗、辨識廣義短路比相對于理論值的誤差,定義t時刻阻抗誤差、和為:

    圖8 對比

    為了進一步研究不同負荷對戴維南參數(shù)與時變廣義短路比的影響,分別在單饋入直流系統(tǒng)節(jié)點2處接入恒阻抗、恒電流、恒功率、感應(yīng)電機及綜合負荷,計算五種負荷下系統(tǒng)的等值阻抗以及時變廣義短路比,如圖9所示,計算不同類型負荷的戴維南等值阻抗,如圖10所示,分析不同負荷對單饋入直流系統(tǒng)廣義短路比計算的影響。

    圖9 不同類型負荷的戴維南等值阻抗

    圖10 五種負荷類型下的等值系統(tǒng)時變廣義短路比

    由圖9、圖10可知,節(jié)點2處連接不同負荷類型時,戴維南等值阻抗以及時變廣義短路比是不同的,其中恒阻抗、恒電流、恒功率以及感應(yīng)電機負荷對于戴維南等值阻抗以及廣義短路比影響依次增大,這表明了忽略負荷的動態(tài)特性計算廣義短路比是不準確的。

    對于恒阻抗負荷,其暫態(tài)過程中負荷阻抗不變,由于算例中BPA數(shù)據(jù)精度問題,使得負荷等值阻抗變化略有波動,但不超過1%,因此理論廣義短路比基本不變,故恒阻抗負荷特性對于廣義短路比無影響。

    對于恒電流負荷,相較于恒阻抗負荷,其計算的時變廣義短路比總體更小。

    1) 在0.01~0.27s期間,其時變廣義短路比小于恒阻抗負荷的時變廣義短路比。恒電流負荷在暫態(tài)過程中,電流基本不變,而電壓從1.0214減小到0.8835,然后持續(xù)上升至0.9817,電壓差最小在0.26~0.27s之間,為10-4數(shù)量級,使得在此期間由式(7)計算的恒電流負荷等值阻抗恒大于恒阻抗負荷等值阻抗,并與線路聚合得到的戴維南等值阻抗大于恒阻抗負荷戴維南等值阻抗,導(dǎo)致在恒電流負荷情況下的時變廣義短路比小于恒阻抗情況下的時變廣義短路比;

    2) 在0.27~0.3s期間,其時變廣義短路比增加0.1124,其值為2.6238。這是由于BPA輸出精度為10-5,而在此期間節(jié)點2處節(jié)點電壓差很小,實部與虛部范圍分別為10-4數(shù)量級,注入恒電流負荷的電流差也很小,實虛部范圍分別為10-5數(shù)量級,而該時刻廣義短路比2.6238在考慮帶入BPA誤差計算得到時變廣義短路比范圍內(nèi),其范圍為2.5089 ~2.8667。

    對于恒功率負荷,其對于廣義短路影響復(fù)雜,這是由于在為了滿足負荷為2+j1,電壓實部與虛部變化,導(dǎo)致電流實部虛部中的一個或兩個向相反方向變化,使得基于式(7)的負荷阻抗呈負電阻或負電抗特性,需要結(jié)合電壓電流變化分析負荷阻抗的變化,進而廣義短路比變化復(fù)雜。

    1) 在0.01~0.03s期間,廣義短路比從2.5860增加到2.7042然后減小到2.6571,這是由于該期間負荷節(jié)點電壓實部由1.0150上升至1.0179,而電壓虛部由-0.1508減小到-0.1552,電流實部從1.7838減小到1.7736,電流虛部從-1.2498減小到-1.25334,計算得到負荷阻抗實部為負電阻,其中0.02s與其他時刻相比,其電流變化相較于電壓更大,使得基于式(7)計算的負荷等值阻抗在該期間最小。

    2) 在0.03~0.13s期間,廣義短路比變化從2.6701減小到2.6653,隨后增加到0.14s的2.5400,這是由于負荷節(jié)點電壓快速恢復(fù),實部虛部快速變化,其間電壓實部增加了0.0136,電壓虛部減少了0.0051,進而引起電流相反方向變化,其中電流實部減小了-0.0446,電流虛部增大了0.0173,此時由式(7)計算得到恒功率負荷阻抗與線路阻抗聚合后的系統(tǒng)等值阻抗最大,使得廣義短路比在此階段最小。

    3) 在0.14 ~0.25s期間,廣義短路比變化幅度較小,為7%,這是由于在每個窗口內(nèi),電流實部約線性增加0.0041,虛部約線性減小0.0210,對于電壓實部線性增加0.0034,虛部線性增加0.0021,使得計算戴維南等值阻抗變化不超過10%;

    4) 在0.26~0.30s期間,此時廣義短路比相比0.25s最大增加0.3274,這是由于在0.09~0.25s前后負荷節(jié)點電壓變化為10-4數(shù)量級,電流變化為10-3數(shù)量級,與恒電流負荷類似,由BPA輸出精度造成此期間廣義短路比計算范圍應(yīng)為2.2630~2.9662,該期間的廣義短路比仍在此范圍內(nèi)。

    對于感應(yīng)電機負荷,在圖中感應(yīng)電機與綜合負荷時變廣義短路比兩者變化趨勢一致,但與綜合負荷相比,感應(yīng)電機負荷條件下的廣義短路比最小值降低了0.1855,其最大值降低了1.282,這是由于單饋入直流廣義短路比與戴維南等值阻抗成反比關(guān)系,而感應(yīng)電機的戴維南等值阻抗變化比綜合負荷大,其中戴維南等值最大阻抗約為綜合負荷的1.5倍,最小戴維南等值阻抗約為綜合負荷的1.2倍,導(dǎo)致廣義短路比最值減小。這表明感應(yīng)電機占比越大,時變戴維南等值阻抗最大值越大,計算所得時變短路比最小值越小。

    4.2 實際系統(tǒng)算例

    以某省級實際電網(wǎng)為例,驗證本文方法的實用性。在實際四饋入直流系統(tǒng)中,設(shè)置四個饋入點同時發(fā)生三相金屬接地短路故障,故障持續(xù)時間0.1s,四條直流線路的額定功率分別為8000MW,8000MW,8000MW,8000MW。記故障清除時刻為起始時刻,記錄1s的U、S等電氣量作為辨識數(shù)據(jù)。分別計算負荷為恒阻抗和綜合負荷兩種情況下的實際的饋入直流系統(tǒng)時變廣義短路比,與SCCP等值計算的廣義短路比相比較。其中SCCP等值方法基于網(wǎng)絡(luò)參數(shù)與拓撲,將負荷按照恒阻抗處理求取戴維南等值參數(shù),求取計算廣義短路比。

    圖11為三種情況下實際電網(wǎng)的廣義短路比,其中SCCP等值計算的廣義短路比為2.0370。

    圖11 三種計算條件下實際電網(wǎng)的廣義短路比

    由于實際系統(tǒng)中負荷采用綜合負荷模型,為更好解釋系統(tǒng)的廣義短路比變化,將實際系統(tǒng)各直流端口節(jié)點注入有功功率繪制如圖12。

    圖12 實際系統(tǒng)各直流端口節(jié)點注入有功功率

    由圖11、圖12可知:

    1) 恒阻抗下與SCCP等值方法相比,恒阻抗情況下計算的廣義短路比仍時變,且在0~0.2s廣義短路比尚未穩(wěn)定時大于SCCP計算值,而 0.2s以后廣義短路比穩(wěn)定小于SCCP計算值。一方面是由于系統(tǒng)中除了負荷以外,電力系統(tǒng)中其他元件具有動態(tài)特性,并增強了受端交流系統(tǒng)強度,使得0~0.2s的廣義短路比大于SCCP計算值,并以同步機勵磁系統(tǒng)為例,由于強勵磁維持系統(tǒng)發(fā)電機節(jié)點電壓,使得受端電網(wǎng)強度增強;另一方面,由于在系統(tǒng)穩(wěn)定時,本方法將SCCP計算條件未考慮到的元件內(nèi)部阻抗計算在戴維南等值參數(shù)中,使得戴維南等值阻抗整體偏大,導(dǎo)致在0.2s廣義短路比變?。?/p>

    2) 綜合負荷與SCCP條件下廣義短路比相比,更具有波動性,且兩者相差較大,這表明在時變電力系統(tǒng)中SCCP計算誤差較大,不能正確地反映時變系統(tǒng)的廣義短路比,同時表明本文能夠反映非線性負荷對于廣義短路比的影響;

    3) 綜合負荷與恒阻抗計算條件下廣義短路比,時變廣義短路比相差最大為2.8912。這是由于在0 ~0.2s四條直流功率快速恢復(fù),使得負荷節(jié)點電壓快速增加,引起感應(yīng)電機負荷暫態(tài)特性快速響應(yīng),導(dǎo)致感應(yīng)電機負荷等值阻抗小,使得此期間系統(tǒng)廣義短路比主要受感應(yīng)電機負荷動態(tài)特性而變化。

    5 結(jié)語

    本文基于暫態(tài)數(shù)據(jù),準確地辨識出多饋入系統(tǒng)戴維南等值參數(shù),其中包含負荷特性、網(wǎng)絡(luò)參數(shù)等電網(wǎng)信息,有效的避免將電力系統(tǒng)綜合負荷按照恒阻抗方式處理的誤差,同時反映了系統(tǒng)中廣義短路比的變化,提出了一種計算時變廣義短路比的方法。單饋入直流等值算例與實際四饋入等值系統(tǒng)結(jié)果表明,本文提示的方法具有較高的準確性,相比于SCCP等值方法本方法不僅考慮到不同負荷模型對于戴維南參數(shù)的影響,同時將其他元件的動態(tài)特性反映在內(nèi)具有較好的實用性,在計算廣義短路比更為準確,有助于準確量化受端交流系統(tǒng)電網(wǎng)強度。

    猜你喜歡
    等值廣義端口
    Rn中的廣義逆Bonnesen型不等式
    一種端口故障的解決方案
    異步電動機等值負載研究
    防爆電機(2020年5期)2020-12-14 07:03:50
    從廣義心腎不交論治慢性心力衰竭
    端口阻塞與優(yōu)先級
    有限群的廣義交換度
    電網(wǎng)單點等值下等效諧波參數(shù)計算
    初識電腦端口
    電腦迷(2015年6期)2015-05-30 08:52:42
    生成樹協(xié)議實例探討
    基于戴維南等值模型的靜穩(wěn)極限在線監(jiān)視
    欧美日韩在线观看h| 丝袜喷水一区| 欧美高清成人免费视频www| h日本视频在线播放| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲av香蕉五月| 男女边吃奶边做爰视频| 日韩欧美精品免费久久| 秋霞在线观看毛片| 美女被艹到高潮喷水动态| 国产探花极品一区二区| 亚洲av二区三区四区| 成人亚洲欧美一区二区av| 久久久久久九九精品二区国产| 国产亚洲精品av在线| 国产69精品久久久久777片| 亚洲av中文av极速乱| 两个人的视频大全免费| 亚洲国产日韩欧美精品在线观看| 国产精品久久电影中文字幕| 晚上一个人看的免费电影| 亚洲精品456在线播放app| 中文字幕精品亚洲无线码一区| 男人和女人高潮做爰伦理| 成人无遮挡网站| 午夜老司机福利剧场| av卡一久久| 我要搜黄色片| 亚洲精品乱码久久久v下载方式| 一区二区三区免费毛片| 两个人的视频大全免费| 小蜜桃在线观看免费完整版高清| 中文字幕制服av| 别揉我奶头 嗯啊视频| 欧美成人精品欧美一级黄| 桃色一区二区三区在线观看| 可以在线观看毛片的网站| 99久久九九国产精品国产免费| 小说图片视频综合网站| 小说图片视频综合网站| 国产探花在线观看一区二区| 午夜精品在线福利| 色综合站精品国产| av又黄又爽大尺度在线免费看 | 赤兔流量卡办理| 少妇熟女aⅴ在线视频| 白带黄色成豆腐渣| 午夜精品一区二区三区免费看| 一本久久精品| 国产一区二区在线av高清观看| 久久中文看片网| 波多野结衣高清作品| 有码 亚洲区| av卡一久久| 日韩av在线大香蕉| 三级男女做爰猛烈吃奶摸视频| 狂野欧美白嫩少妇大欣赏| 草草在线视频免费看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲性久久影院| 三级国产精品欧美在线观看| 91午夜精品亚洲一区二区三区| 日韩一本色道免费dvd| 少妇被粗大猛烈的视频| 黑人高潮一二区| 青春草亚洲视频在线观看| 国产精品久久视频播放| 99九九线精品视频在线观看视频| 国产精品av视频在线免费观看| 我的女老师完整版在线观看| 一边亲一边摸免费视频| 亚洲无线在线观看| 高清毛片免费观看视频网站| 嘟嘟电影网在线观看| 内地一区二区视频在线| 乱码一卡2卡4卡精品| 麻豆一二三区av精品| 看免费成人av毛片| 91aial.com中文字幕在线观看| 亚洲最大成人av| 亚洲成人久久爱视频| 国产日韩欧美在线精品| 三级经典国产精品| 亚洲经典国产精华液单| 国产v大片淫在线免费观看| 蜜桃久久精品国产亚洲av| 日韩欧美三级三区| 亚洲av中文字字幕乱码综合| 自拍偷自拍亚洲精品老妇| 国产三级在线视频| 国产单亲对白刺激| 晚上一个人看的免费电影| 亚洲在线自拍视频| 女同久久另类99精品国产91| 亚洲国产精品sss在线观看| 99在线人妻在线中文字幕| 最后的刺客免费高清国语| 久久久国产成人精品二区| 国内精品久久久久精免费| 变态另类成人亚洲欧美熟女| 午夜a级毛片| 综合色丁香网| 国产精品久久久久久精品电影小说 | 三级经典国产精品| 草草在线视频免费看| 一个人免费在线观看电影| 亚洲欧美中文字幕日韩二区| 国产在线精品亚洲第一网站| 九九在线视频观看精品| 亚州av有码| 美女内射精品一级片tv| 最近最新中文字幕大全电影3| 国语自产精品视频在线第100页| 免费观看a级毛片全部| 黄色视频,在线免费观看| 少妇人妻精品综合一区二区 | www日本黄色视频网| av在线天堂中文字幕| 狂野欧美白嫩少妇大欣赏| 免费观看精品视频网站| 日韩av不卡免费在线播放| 国产av不卡久久| 日本欧美国产在线视频| 亚洲在线自拍视频| 欧美+亚洲+日韩+国产| 亚洲精品粉嫩美女一区| 免费在线观看成人毛片| 91狼人影院| av天堂在线播放| 亚洲av成人av| 乱系列少妇在线播放| 一级毛片aaaaaa免费看小| 身体一侧抽搐| 丰满乱子伦码专区| 非洲黑人性xxxx精品又粗又长| 国产精品野战在线观看| 人人妻人人澡欧美一区二区| av在线播放精品| 蜜桃亚洲精品一区二区三区| 亚洲av二区三区四区| avwww免费| 日日干狠狠操夜夜爽| 桃色一区二区三区在线观看| 黄色配什么色好看| 久久精品国产清高在天天线| 乱系列少妇在线播放| 午夜视频国产福利| 99精品在免费线老司机午夜| 狠狠狠狠99中文字幕| 国产大屁股一区二区在线视频| 亚洲精品国产成人久久av| 亚洲精品久久久久久婷婷小说 | 97超视频在线观看视频| 国产伦理片在线播放av一区 | 国产精品99久久久久久久久| 日韩视频在线欧美| 高清午夜精品一区二区三区 | 日韩欧美在线乱码| 3wmmmm亚洲av在线观看| 日韩欧美精品免费久久| 国产精品野战在线观看| eeuss影院久久| 久久韩国三级中文字幕| 亚洲欧洲国产日韩| 高清毛片免费观看视频网站| 成年av动漫网址| 婷婷精品国产亚洲av| 亚洲精品日韩av片在线观看| 日韩人妻高清精品专区| 久久久久久大精品| 一进一出抽搐gif免费好疼| 国产精品久久电影中文字幕| 国产精品永久免费网站| 亚洲av不卡在线观看| 免费在线观看成人毛片| 国产真实乱freesex| 久久亚洲精品不卡| 国产熟女欧美一区二区| 国产伦一二天堂av在线观看| 九九爱精品视频在线观看| 毛片一级片免费看久久久久| 久久精品综合一区二区三区| 亚洲国产精品sss在线观看| 天天躁夜夜躁狠狠久久av| 12—13女人毛片做爰片一| 国产精品.久久久| 久久久久久九九精品二区国产| 久久久久久国产a免费观看| 久久精品国产亚洲网站| 网址你懂的国产日韩在线| 久久久久久伊人网av| 久久久国产成人免费| 美女被艹到高潮喷水动态| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久av不卡| 国产探花在线观看一区二区| 精品人妻视频免费看| 一个人观看的视频www高清免费观看| 日韩欧美国产在线观看| 日日啪夜夜撸| 久久午夜福利片| 国产日韩欧美在线精品| 国产精品无大码| 国产精品久久视频播放| 国产精品不卡视频一区二区| www.色视频.com| 亚洲最大成人手机在线| 搞女人的毛片| 欧美精品国产亚洲| 国产v大片淫在线免费观看| .国产精品久久| 97超视频在线观看视频| 国产精品久久久久久久久免| 色哟哟哟哟哟哟| 亚洲最大成人av| 亚洲久久久久久中文字幕| 在线天堂最新版资源| 久久久久久国产a免费观看| 国产一级毛片在线| 日韩欧美 国产精品| 中国国产av一级| 欧美另类亚洲清纯唯美| 22中文网久久字幕| 亚洲精品久久国产高清桃花| 久久亚洲国产成人精品v| 男人舔女人下体高潮全视频| 能在线免费观看的黄片| 黄色日韩在线| 中文字幕熟女人妻在线| 中文在线观看免费www的网站| 色哟哟哟哟哟哟| 亚洲av一区综合| 国产综合懂色| 精品久久久久久久久亚洲| 给我免费播放毛片高清在线观看| 亚洲在久久综合| 菩萨蛮人人尽说江南好唐韦庄 | 成人无遮挡网站| 成人一区二区视频在线观看| 日韩精品有码人妻一区| 熟女人妻精品中文字幕| 在线观看66精品国产| 麻豆国产97在线/欧美| 成人二区视频| 国产免费一级a男人的天堂| 18禁在线播放成人免费| 天堂网av新在线| 最好的美女福利视频网| 日本欧美国产在线视频| 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 国产大屁股一区二区在线视频| 久久久午夜欧美精品| 成人无遮挡网站| 欧美高清成人免费视频www| 男人舔女人下体高潮全视频| 麻豆精品久久久久久蜜桃| 免费av不卡在线播放| 国产一区二区在线av高清观看| 久久精品国产99精品国产亚洲性色| 在线观看免费视频日本深夜| 综合色丁香网| 看黄色毛片网站| 搡女人真爽免费视频火全软件| 在线免费十八禁| 蜜桃亚洲精品一区二区三区| 高清毛片免费看| 成年免费大片在线观看| 国产麻豆成人av免费视频| 国产成人a∨麻豆精品| 亚洲av一区综合| 1024手机看黄色片| 婷婷色av中文字幕| kizo精华| 欧美+亚洲+日韩+国产| 国产av不卡久久| 久久久久性生活片| 成人毛片60女人毛片免费| 久久韩国三级中文字幕| 天天一区二区日本电影三级| 国产视频首页在线观看| 国产成人a区在线观看| 26uuu在线亚洲综合色| 国产老妇伦熟女老妇高清| 三级经典国产精品| av视频在线观看入口| 日韩一区二区视频免费看| 美女国产视频在线观看| 亚洲国产精品久久男人天堂| 国产精品一区二区在线观看99 | 久久草成人影院| 日韩一区二区视频免费看| 国产成年人精品一区二区| 亚洲精品粉嫩美女一区| 成人性生交大片免费视频hd| 大又大粗又爽又黄少妇毛片口| 欧美日韩一区二区视频在线观看视频在线 | 国产成人a∨麻豆精品| 国产69精品久久久久777片| 成人亚洲精品av一区二区| 亚洲熟妇中文字幕五十中出| 禁无遮挡网站| 一本精品99久久精品77| 亚洲av二区三区四区| 在线观看一区二区三区| 久久6这里有精品| 欧美极品一区二区三区四区| 在线观看66精品国产| 免费人成在线观看视频色| 亚洲欧美成人精品一区二区| 99久久久亚洲精品蜜臀av| 国产国拍精品亚洲av在线观看| 久99久视频精品免费| 天堂av国产一区二区熟女人妻| 99久久九九国产精品国产免费| 国产精品av视频在线免费观看| 听说在线观看完整版免费高清| 国产精品一区www在线观看| 国产精品野战在线观看| 波多野结衣高清无吗| 黄片wwwwww| 丝袜美腿在线中文| 天堂网av新在线| 蜜臀久久99精品久久宅男| 久久人人爽人人爽人人片va| 国产老妇伦熟女老妇高清| 尤物成人国产欧美一区二区三区| 日本熟妇午夜| 国产三级在线视频| 欧美性猛交╳xxx乱大交人| 亚洲最大成人手机在线| 国产探花在线观看一区二区| 天天躁夜夜躁狠狠久久av| 给我免费播放毛片高清在线观看| 夜夜看夜夜爽夜夜摸| 亚洲国产精品成人综合色| 尾随美女入室| 最近视频中文字幕2019在线8| 一级毛片aaaaaa免费看小| 日韩欧美三级三区| 久久久久久久亚洲中文字幕| 日韩制服骚丝袜av| 高清在线视频一区二区三区 | 日韩欧美 国产精品| 禁无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 一个人观看的视频www高清免费观看| 熟妇人妻久久中文字幕3abv| 偷拍熟女少妇极品色| 成人午夜精彩视频在线观看| 九色成人免费人妻av| 丝袜喷水一区| 丰满人妻一区二区三区视频av| 狂野欧美白嫩少妇大欣赏| 色哟哟哟哟哟哟| 久久国内精品自在自线图片| 国产成人精品久久久久久| 最近手机中文字幕大全| 九九在线视频观看精品| 禁无遮挡网站| 国产免费一级a男人的天堂| 久久国内精品自在自线图片| 日本黄大片高清| 18禁裸乳无遮挡免费网站照片| 白带黄色成豆腐渣| 18禁在线无遮挡免费观看视频| 亚洲最大成人中文| а√天堂www在线а√下载| av专区在线播放| 丰满的人妻完整版| 听说在线观看完整版免费高清| 精品午夜福利在线看| 久久久色成人| 欧美+日韩+精品| 日韩三级伦理在线观看| 久久午夜亚洲精品久久| 国产综合懂色| 边亲边吃奶的免费视频| 国产黄片视频在线免费观看| 又爽又黄a免费视频| 国产女主播在线喷水免费视频网站 | 长腿黑丝高跟| 三级毛片av免费| 国产熟女欧美一区二区| 五月玫瑰六月丁香| 麻豆av噜噜一区二区三区| 久久精品91蜜桃| 精品久久久久久久久av| 午夜激情欧美在线| 久久久久免费精品人妻一区二区| 只有这里有精品99| av视频在线观看入口| 国产综合懂色| 日韩欧美精品v在线| 99国产精品一区二区蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看| 国产黄色小视频在线观看| 欧美日本视频| 在线观看美女被高潮喷水网站| 中文资源天堂在线| 欧美日韩乱码在线| 老熟妇乱子伦视频在线观看| 嫩草影院入口| 欧美区成人在线视频| 一本久久中文字幕| 国产女主播在线喷水免费视频网站 | 男的添女的下面高潮视频| 在线免费十八禁| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 国产乱人视频| 少妇熟女欧美另类| 日本色播在线视频| 99久久人妻综合| 久久久久久久久大av| 99久久成人亚洲精品观看| 免费电影在线观看免费观看| 久久久久久久亚洲中文字幕| 成年女人永久免费观看视频| 中国美白少妇内射xxxbb| 青青草视频在线视频观看| 日韩国内少妇激情av| 99国产精品一区二区蜜桃av| 国产老妇女一区| 99久久精品国产国产毛片| 亚洲电影在线观看av| 精品一区二区三区视频在线| 12—13女人毛片做爰片一| 亚洲成人久久爱视频| 神马国产精品三级电影在线观看| 夜夜夜夜夜久久久久| av国产免费在线观看| 中文字幕制服av| 菩萨蛮人人尽说江南好唐韦庄 | 免费一级毛片在线播放高清视频| 听说在线观看完整版免费高清| 国产精品99久久久久久久久| 国产精品久久久久久av不卡| 亚洲第一区二区三区不卡| 在现免费观看毛片| 日韩人妻高清精品专区| 久久热精品热| 人妻系列 视频| 久久人人精品亚洲av| 国产精品蜜桃在线观看 | 99在线视频只有这里精品首页| 国内少妇人妻偷人精品xxx网站| 韩国av在线不卡| 97人妻精品一区二区三区麻豆| 色视频www国产| 12—13女人毛片做爰片一| 少妇裸体淫交视频免费看高清| 精品久久久久久久久av| 美女高潮的动态| 永久网站在线| а√天堂www在线а√下载| 99热这里只有是精品50| 亚洲精品日韩av片在线观看| 久久午夜亚洲精品久久| 午夜爱爱视频在线播放| 亚洲中文字幕一区二区三区有码在线看| 在线免费十八禁| av在线观看视频网站免费| 亚洲精品456在线播放app| 五月伊人婷婷丁香| 高清在线视频一区二区三区 | 男女做爰动态图高潮gif福利片| 尤物成人国产欧美一区二区三区| 少妇的逼好多水| 日韩欧美在线乱码| 91久久精品国产一区二区成人| 插逼视频在线观看| 久久精品夜色国产| 边亲边吃奶的免费视频| 欧美在线一区亚洲| 性色avwww在线观看| 国产在线男女| 国产成人91sexporn| 国产v大片淫在线免费观看| 色哟哟哟哟哟哟| 最近手机中文字幕大全| 精品少妇黑人巨大在线播放 | 色综合亚洲欧美另类图片| 国产黄色视频一区二区在线观看 | 国产女主播在线喷水免费视频网站 | 麻豆成人av视频| 少妇猛男粗大的猛烈进出视频 | 精品一区二区三区视频在线| 村上凉子中文字幕在线| 国产久久久一区二区三区| 免费电影在线观看免费观看| 午夜免费男女啪啪视频观看| 亚洲国产精品成人综合色| 亚洲第一区二区三区不卡| 一级黄片播放器| 色噜噜av男人的天堂激情| 日本一本二区三区精品| 成人特级av手机在线观看| 久久久久国产网址| 亚洲av免费高清在线观看| 日本五十路高清| 欧洲精品卡2卡3卡4卡5卡区| 国产成年人精品一区二区| 91精品一卡2卡3卡4卡| 一区二区三区四区激情视频 | 欧美色欧美亚洲另类二区| 黄片wwwwww| 亚洲在线观看片| 激情 狠狠 欧美| 国产精品一区二区在线观看99 | 变态另类丝袜制服| 久久久久久久久久成人| 97热精品久久久久久| 日韩人妻高清精品专区| 精品99又大又爽又粗少妇毛片| 亚洲国产欧洲综合997久久,| 国产免费男女视频| 久久九九热精品免费| 最好的美女福利视频网| 午夜久久久久精精品| 91久久精品国产一区二区三区| 老熟妇乱子伦视频在线观看| 国产真实乱freesex| 精华霜和精华液先用哪个| 亚洲av成人精品一区久久| av在线天堂中文字幕| 三级男女做爰猛烈吃奶摸视频| 国产黄色视频一区二区在线观看 | 卡戴珊不雅视频在线播放| 国产女主播在线喷水免费视频网站 | av又黄又爽大尺度在线免费看 | 色综合色国产| 成人av在线播放网站| 国产一区二区在线观看日韩| 国产成人a区在线观看| 欧美一区二区亚洲| 日本黄色视频三级网站网址| 欧美3d第一页| 亚洲自偷自拍三级| 日韩国内少妇激情av| 亚洲四区av| 自拍偷自拍亚洲精品老妇| 中文在线观看免费www的网站| 噜噜噜噜噜久久久久久91| 国产白丝娇喘喷水9色精品| 久久久久网色| 免费观看a级毛片全部| 亚洲欧美日韩高清在线视频| 91在线精品国自产拍蜜月| 亚洲国产高清在线一区二区三| 日本黄色片子视频| 国产精品一及| 欧美日韩国产亚洲二区| 久久久久九九精品影院| 国产亚洲5aaaaa淫片| 国产大屁股一区二区在线视频| 日韩强制内射视频| 亚洲人成网站在线播放欧美日韩| 韩国av在线不卡| 我的女老师完整版在线观看| 亚洲精华国产精华液的使用体验 | 成人永久免费在线观看视频| 综合色丁香网| 亚洲av二区三区四区| 97热精品久久久久久| 日本黄色视频三级网站网址| 美女cb高潮喷水在线观看| 午夜福利在线观看免费完整高清在 | 欧美丝袜亚洲另类| 日韩人妻高清精品专区| 日韩一区二区三区影片| 亚洲成人精品中文字幕电影| 性欧美人与动物交配| 精品久久久久久久久亚洲| 色尼玛亚洲综合影院| 中文精品一卡2卡3卡4更新| 久久精品国产亚洲av涩爱 | 日产精品乱码卡一卡2卡三| 亚洲国产精品国产精品| 成人漫画全彩无遮挡| 亚洲第一区二区三区不卡| 干丝袜人妻中文字幕| 国产成人精品婷婷| 国产白丝娇喘喷水9色精品| 可以在线观看的亚洲视频| 精品一区二区免费观看| 亚洲性久久影院| 国产亚洲精品久久久com| 婷婷色综合大香蕉| 午夜视频国产福利| 国产成人影院久久av| 午夜福利在线在线| 国产精品麻豆人妻色哟哟久久 | 色综合站精品国产| 国产精品一区二区在线观看99 | 中国国产av一级| 亚洲第一区二区三区不卡| 亚洲av不卡在线观看| 亚洲无线在线观看| 在线观看免费视频日本深夜| 国产单亲对白刺激| 国产国拍精品亚洲av在线观看| h日本视频在线播放| 夜夜夜夜夜久久久久| 又粗又硬又长又爽又黄的视频 | 99在线人妻在线中文字幕| 国产精品一区二区三区四区免费观看| 给我免费播放毛片高清在线观看| 看非洲黑人一级黄片| 成人午夜精彩视频在线观看| 成年免费大片在线观看| 国产精品一区二区三区四区久久| 老熟妇乱子伦视频在线观看| 欧美色欧美亚洲另类二区| 欧美精品国产亚洲|