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

    電力網(wǎng)絡(luò)的諧振穩(wěn)定性分析方法研究

    2017-11-27 07:07:00徐政王世佳邢法財(cái)肖晃慶
    電力建設(shè) 2017年11期
    關(guān)鍵詞:系統(tǒng)

    徐政,王世佳,邢法財(cái),肖晃慶

    (浙江大學(xué)電氣工程學(xué)院,杭州市 310027)

    電力網(wǎng)絡(luò)的諧振穩(wěn)定性分析方法研究

    徐政,王世佳,邢法財(cái),肖晃慶

    (浙江大學(xué)電氣工程學(xué)院,杭州市 310027)

    隨著電力系統(tǒng)電力電子化程度的不斷加深,近年來出現(xiàn)了多起機(jī)理不明的新的振蕩現(xiàn)象。提出了電力網(wǎng)絡(luò)諧振穩(wěn)定性的概念,試圖將上述機(jī)理不明的振蕩現(xiàn)象納入到電力網(wǎng)絡(luò)的諧振不穩(wěn)定范疇,從而基于線性網(wǎng)絡(luò)理論對(duì)眾多復(fù)雜振蕩現(xiàn)象進(jìn)行分析,為借助數(shù)學(xué)上完全成熟的線性系統(tǒng)理論解決電力系統(tǒng)實(shí)際問題提供一條途徑。圍繞判斷諧振穩(wěn)定性的分析方法展開,通過引入電力網(wǎng)絡(luò)的s域節(jié)點(diǎn)導(dǎo)納矩陣,將電力網(wǎng)絡(luò)的諧振穩(wěn)定性問題歸結(jié)為判斷s域節(jié)點(diǎn)導(dǎo)納矩陣行列式的零點(diǎn)在復(fù)平面上的分布問題。首先,理論證明了s域節(jié)點(diǎn)導(dǎo)納矩陣行列式的零點(diǎn)就是系統(tǒng)的特征值。其次,給出了求解s域節(jié)點(diǎn)導(dǎo)納矩陣行列式零點(diǎn)的實(shí)部-虛部交叉迭代法。接著,推導(dǎo)了特定諧振模式下的節(jié)點(diǎn)電壓振型和參與因子矩陣,這2個(gè)指標(biāo)可用來定位特定諧振模式發(fā)生的位置。最后,通過算例展示了所提方法在分析電力網(wǎng)絡(luò)諧振穩(wěn)定性方面的有效性。

    電力網(wǎng)絡(luò);狀態(tài)空間模型;s域節(jié)點(diǎn)導(dǎo)納矩陣;諧振模式;諧振穩(wěn)定性;節(jié)點(diǎn)電壓振型;參與因子矩陣

    0 引 言

    振蕩是電力系統(tǒng)運(yùn)行過程中的一種常見現(xiàn)象[1-2]。一般來說,電力系統(tǒng)振蕩可以分為3種類型:(1)發(fā)電機(jī)軸系的振蕩(扭振);(2)發(fā)電機(jī)轉(zhuǎn)子之間的振蕩;(3)電力網(wǎng)絡(luò)內(nèi)部的固有諧振。隨著大量電力電子裝置接入電力系統(tǒng),電力網(wǎng)絡(luò)的固有諧振問題變得更加嚴(yán)重和復(fù)雜。由于某些電力電子裝置在一定的頻段內(nèi)存在負(fù)電阻效應(yīng),使得原先諧振穩(wěn)定的電力網(wǎng)絡(luò)有可能變得諧振不穩(wěn)定,從而造成嚴(yán)重的后果。例如,2011年至今,河北沽源地區(qū)的雙饋風(fēng)電場(chǎng)與其電容串補(bǔ)裝置相互作用,引發(fā)了上百次頻率為 3~10 Hz的次同步振蕩,造成變壓器異常振動(dòng)和大量風(fēng)機(jī)脫網(wǎng)[3]。福建廈門柔性直流輸電工程在直流側(cè)出現(xiàn)過次同步振蕩問題,該振蕩現(xiàn)象為非衰減和非等幅振蕩,振蕩頻率為25 Hz左右[4]。新疆哈密“7·1”直驅(qū)風(fēng)機(jī)引起的火電機(jī)組次同步振蕩事件[5],可以認(rèn)為是由電力網(wǎng)絡(luò)諧振與發(fā)電機(jī)軸系扭振相互作用而導(dǎo)致的機(jī)網(wǎng)復(fù)合共振問題。

    當(dāng)電力網(wǎng)絡(luò)發(fā)生諧振時(shí),很容易引起系統(tǒng)的過電壓、過電流,甚至?xí)c發(fā)電機(jī)軸系扭振相互作用而導(dǎo)致機(jī)網(wǎng)復(fù)合共振,威脅電力系統(tǒng)的安全穩(wěn)定運(yùn)行。因此,亟需對(duì)電力網(wǎng)絡(luò)的固有諧振結(jié)構(gòu)和諧振穩(wěn)定性進(jìn)行透徹的研究。

    本文深入研究基于s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)的電力網(wǎng)絡(luò)諧振穩(wěn)定性分析方法。首先,在理論上論證s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)行列式為0的根就是系統(tǒng)的特征值,并用一個(gè)簡(jiǎn)單算例進(jìn)行了驗(yàn)證。其次,提出一種實(shí)部-虛部交叉迭代的方法,用以確定s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)行列式為0的根的位置;這種求解方法與采用Newton-Raphson法迭代求解det[Y(s)]零點(diǎn)的方法相比,可以大幅度減少計(jì)算量。接著,推導(dǎo)特定諧振模式下的節(jié)點(diǎn)電壓振型和參與因子矩陣,這2個(gè)指標(biāo)可用來定位特定諧振模式發(fā)生的位置。最后,通過IEEE次同步諧振第一標(biāo)準(zhǔn)測(cè)試系統(tǒng)和IEEE 39節(jié)點(diǎn)測(cè)試系統(tǒng),展示對(duì)實(shí)際電網(wǎng)進(jìn)行諧振穩(wěn)定性分析的結(jié)果。

    1 理論基礎(chǔ)

    1.1 諧振穩(wěn)定性的定義

    當(dāng)考慮輸電線路等分布參數(shù)元件時(shí),因描述分布參數(shù)元件特性的方程是偏微分方程,整個(gè)電力網(wǎng)絡(luò)已不能用線性定常系統(tǒng)的標(biāo)準(zhǔn)狀態(tài)空間模型來描述。另外,若進(jìn)一步考慮元件參數(shù)隨頻率而變化的特性,那么即使對(duì)于由集總參數(shù)元件構(gòu)成的電力網(wǎng)絡(luò),也無法用線性定常系統(tǒng)的標(biāo)準(zhǔn)狀態(tài)空間模型來描述。因此,當(dāng)電力網(wǎng)絡(luò)包含分布參數(shù)元件和頻變參數(shù)元件時(shí),所謂的電力網(wǎng)絡(luò)諧振穩(wěn)定性該如何定義就成為一個(gè)重要問題。

    1999年,文獻(xiàn)[7]首先提出了采用s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)來分析復(fù)雜電力網(wǎng)絡(luò)小信號(hào)穩(wěn)定性的概念和方法,這里所謂的“復(fù)雜電力網(wǎng)絡(luò)小信號(hào)穩(wěn)定性”本質(zhì)上與本文所稱的“電力網(wǎng)絡(luò)諧振穩(wěn)定性”相一致。因此,對(duì)于包含分布參數(shù)元件和頻變參數(shù)元件的復(fù)雜電力網(wǎng)絡(luò),諧振穩(wěn)定性的概念將基于s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)來定義。2001年,文獻(xiàn)[8]對(duì)采用s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)進(jìn)行小信號(hào)穩(wěn)定性分析的方法做了進(jìn)一步的發(fā)展和完善。其要點(diǎn)如下:(1)所謂的“s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)”,國內(nèi)也稱為“運(yùn)算導(dǎo)納矩陣”,是早已存在的概念;例如,電容C的運(yùn)算導(dǎo)納是sC,電感L的運(yùn)算導(dǎo)納是1/(sL);簡(jiǎn)單地說,將交流穩(wěn)態(tài)分析時(shí)元件導(dǎo)納模型中的j用s來替換就構(gòu)成了對(duì)應(yīng)元件的運(yùn)算導(dǎo)納,對(duì)于分布參數(shù)的輸電線路,也有類似的結(jié)果[8];在得到各元件的運(yùn)算導(dǎo)納模型后,構(gòu)建運(yùn)算導(dǎo)納矩陣的步驟與交流穩(wěn)態(tài)分析時(shí)構(gòu)建節(jié)點(diǎn)導(dǎo)納矩陣的步驟完全一致。(2)對(duì)于包含分布參數(shù)元件和頻變參數(shù)元件的電力網(wǎng)絡(luò),構(gòu)建其s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)并不存在特殊困難;因而對(duì)于一般性的電力網(wǎng)絡(luò),基于s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)分析其諧振穩(wěn)定性具有普遍的適用性。(3)令Y(s)的行列式為det[Y(s)],那么det[Y(s)]=0的根(以下統(tǒng)稱為det[Y(s)]的零點(diǎn))就是該電力網(wǎng)絡(luò)的諧振模式,det[Y(s)]的所有零點(diǎn)就是該電力網(wǎng)絡(luò)的所有諧振模式,如果det[Y(s)]的所有零點(diǎn)都位于復(fù)平面的左半平面,那么該電力網(wǎng)絡(luò)就是諧振穩(wěn)定的。容易證明,包含分布參數(shù)元件的電力網(wǎng)絡(luò),其諧振模式有無限個(gè)。

    由上面的介紹可知,判斷電力網(wǎng)絡(luò)諧振穩(wěn)定性與判斷det[Y(s)]的所有零點(diǎn)是否都位于復(fù)平面的左半平面等價(jià)。因此,判斷電力網(wǎng)絡(luò)諧振穩(wěn)定性的最直接方法就是求出det[Y(s)]的所有零點(diǎn),或者求出det[Y(s)]在指定頻段內(nèi)的所有零點(diǎn)。更進(jìn)一步,可以對(duì)電力網(wǎng)絡(luò)的諧振結(jié)構(gòu)進(jìn)行分析,電力網(wǎng)絡(luò)的諧振結(jié)構(gòu)包含4方面的信息:(1)電力網(wǎng)絡(luò)在分析的頻段內(nèi)存在哪些固有諧振模式;(2)各諧振模式的頻率;(3)各諧振模式的阻尼;(4)各諧振模式的振型。文獻(xiàn)[8]提出了采用Newton-Raphson法迭代求解det[Y(s)]零點(diǎn)的方法。該方法需要同時(shí)建立2個(gè)s域節(jié)點(diǎn)導(dǎo)納矩陣,一個(gè)是Y(s),一個(gè)是Y(s)關(guān)于s的導(dǎo)數(shù)矩陣dY(s)/ds,計(jì)算量比較大。

    1.2 系統(tǒng)特征值的定義

    電力網(wǎng)絡(luò)的動(dòng)態(tài)特性可以通過KVL電壓方程、KCL電流方程和儲(chǔ)能元件動(dòng)態(tài)特性方程進(jìn)行描述;其中前2類方程為代數(shù)方程,后1類方程為一階微分方程。當(dāng)電力網(wǎng)絡(luò)由集總參數(shù)元件構(gòu)成,且元件參數(shù)不隨頻率變化時(shí),根據(jù)上述3類方程可以列寫出單輸入單輸出系統(tǒng)如下的狀態(tài)空間模型[9-10]:

    (1)

    式中:{A,T}共同構(gòu)成狀態(tài)空間的系統(tǒng)矩陣;x為狀態(tài)變量向量,包括儲(chǔ)能元件狀態(tài)(電容電壓、電感電流等);輸出y為第k個(gè)節(jié)點(diǎn)的電壓;輸入u為第j個(gè)節(jié)點(diǎn)的注入電流;b是一維常數(shù)列向量;c是一維常數(shù)行向量。

    對(duì)式(1)進(jìn)行拉氏變換,有

    (2)

    因此s域下轉(zhuǎn)移阻抗zkj(s)的表達(dá)式為

    (3)

    式中:分母det(sT-A)表示矩陣sT-A的行列式,“*”表示相應(yīng)矩陣的伴隨矩陣。根據(jù)線性代數(shù)的基本理論可知,det(sT-A)是系統(tǒng)的特征多項(xiàng)式,det(sT-A)=0的根是系統(tǒng)的特征值。系統(tǒng)特征值對(duì)應(yīng)于電力網(wǎng)絡(luò)的所有諧振模式;當(dāng)特征值都位于復(fù)平面的左半平面時(shí),該電力網(wǎng)絡(luò)就是諧振穩(wěn)定的。由于狀態(tài)空間模型需要在集總參數(shù)元件的條件下列寫,此時(shí)諧振模式為有限個(gè)。

    1.3 基于s域節(jié)點(diǎn)導(dǎo)納矩陣求取特征值

    除狀態(tài)空間模型外,可以采用s域節(jié)點(diǎn)導(dǎo)納矩陣分析系統(tǒng)特征值[7-8]。s域下系統(tǒng)節(jié)點(diǎn)電壓方程為

    Y(s)Vnode(s)=Inode(s)

    (4)

    式中:Y(s)為s域下的節(jié)點(diǎn)導(dǎo)納矩陣;Vnode(s)和Inode(s)分別為s域下的節(jié)點(diǎn)電壓向量和節(jié)點(diǎn)注入電流向量。同樣考慮單輸入單輸出系統(tǒng),設(shè)在節(jié)點(diǎn)j上的注入電流為ij(s),在節(jié)點(diǎn)k上的輸出電壓為vk(s)。則方程(4)可以寫成如下形式:

    (5)

    這里b′是一維列向量,其第j個(gè)元素為1,其余元素為0;而c′是一維行向量,其第k個(gè)元素為1,其余元素為0。因此s域下轉(zhuǎn)移阻抗zkj(s)的表達(dá)式為

    (6)

    式中:det[Y(s)]為Y(s)的行列式;Y(s)*為Y(s)的伴隨矩陣。

    由于Y(s)的元素不一定都是s的多項(xiàng)式,例如對(duì)于電感元件L,其運(yùn)算導(dǎo)納的表達(dá)式是1/(sL),因此通過對(duì)det[Y(s)]外乘s的冪函數(shù)sm(m為正整數(shù)),一定能夠?qū)et[Y(s)]化成s的多項(xiàng)式。對(duì)比式(3)和式(6),det(sT-A)與smdet[Y(s)]都是系統(tǒng)的特征多項(xiàng)式。這樣,在s≠0的條件下,det(sT-A)=0與det[Y(s)]=0的根是完全一致的。因此,det[Y(s)]=0的非零根一定是系統(tǒng)的特征值,從而論證了采用狀態(tài)空間法和s域節(jié)點(diǎn)導(dǎo)納矩陣法在求取系統(tǒng)特征值上具有一致性。

    值得指出的是,在構(gòu)建s域節(jié)點(diǎn)導(dǎo)納矩陣Y(s)時(shí),并不要求網(wǎng)絡(luò)元件使用集總參數(shù)模型。因此對(duì)于一般性的電力網(wǎng)絡(luò)(包括分布參數(shù)元件和頻變參數(shù)元件),相對(duì)于狀態(tài)空間法,基于s域節(jié)點(diǎn)導(dǎo)納矩陣分析諧振穩(wěn)定性具有更為普遍的適用性。容易證明,包含分布參數(shù)元件的電力網(wǎng)絡(luò),其諧振模式有無限個(gè)。

    1.4 實(shí)例驗(yàn)證

    圖1 RLC串聯(lián)電路Fig.1 RLC series circuit

    圖1所示RLC電路的狀態(tài)空間方程為

    (7)

    進(jìn)而可以求出系統(tǒng)的特征方程為

    (8)

    若采用s域節(jié)點(diǎn)導(dǎo)納矩陣法,則該電路的s域節(jié)點(diǎn)導(dǎo)納矩陣為:

    (9)

    (10)

    顯然,在s≠0的條件下,式(8)和式(10)具有相同的零點(diǎn)。這說明采用狀態(tài)空間法或s域節(jié)點(diǎn)導(dǎo)納矩陣法進(jìn)行電力網(wǎng)絡(luò)的諧振穩(wěn)定性分析,所得結(jié)果是完全一致的。

    2 基于Y(s)的特征值求解算法

    以上分析可知,電力網(wǎng)絡(luò)特征值求取問題可轉(zhuǎn)化為求解s域節(jié)點(diǎn)導(dǎo)納矩陣行列式等于0的根,即求解如下方程的根:

    det[Y(s)]=0

    (11)

    由于s是復(fù)平面上的變量,直接求解上述方程是一件非常困難的事情?,F(xiàn)有文獻(xiàn)主要采用Newton-Raphson法迭代求解,需要求取Y(s)關(guān)于s的導(dǎo)數(shù)矩陣[8]。為了能夠快速準(zhǔn)確地求出det(Y(s))的零點(diǎn),本文提出實(shí)部-虛部交叉迭代的求解方法。

    2.1 算法原理

    當(dāng)Y(s)行列式不等于0時(shí),其行列式的倒數(shù)為一個(gè)有限值;當(dāng)Y(s)行列式趨近于0時(shí),其行列式的倒數(shù)趨向于無窮大。故定義hm(s)為:

    (12)

    hm(s)在某個(gè)特征值附近的示意圖如圖2所示。需要注意的是,該曲面的頂點(diǎn)在復(fù)平面(σ-jω平面)上的投影就是相應(yīng)的特征值。

    當(dāng)s的實(shí)部固定為σi,虛部在一定范圍內(nèi)變化時(shí),如圖2所示,相當(dāng)于平面σ=σi與曲面hm(s)相交的曲線;找出該曲線極大值所對(duì)應(yīng)的虛部坐標(biāo),記為jωi。接下來,將s的虛部固定為jωi,而實(shí)部在一定范圍內(nèi)變化,如圖3所示,相當(dāng)于平面jω=jωi與曲面hm(s)的交線;找出這條交線的極大值點(diǎn)對(duì)應(yīng)的實(shí)部坐標(biāo),記為σi+1。

    重復(fù)上述過程,直到所得的實(shí)部和虛部不再改變,或者達(dá)到所設(shè)定的精度要求為止。最后求出的實(shí)部和虛部的組合就是s域下系統(tǒng)的特征值。在電力網(wǎng)絡(luò)諧振穩(wěn)定性分析時(shí),系統(tǒng)特征值又稱為諧振模式。因此一旦求出系統(tǒng)特征值,則對(duì)應(yīng)此諧振模式的阻尼和頻率就是已知的。

    圖2 平面σ=σi與曲面相交情況Fig.2 Intersection of σ=σi and tapered surface

    圖3 平面jω=jωi與曲面相交情況Fig.3 Intersection of jω=jωi and tapered surface

    2.2 算法流程圖

    上述的求解思路只是針對(duì)1個(gè)特征值,對(duì)含有多個(gè)特征值的系統(tǒng),其求解過程類似。電力網(wǎng)絡(luò)特征值的詳細(xì)求解流程如圖4所示。

    3 諧振模式的節(jié)點(diǎn)電壓振型和參與因子矩陣

    借鑒文獻(xiàn)[11]的做法,可以定義諧振模式sk的節(jié)點(diǎn)電壓振型和參與因子矩陣。設(shè)電力網(wǎng)絡(luò)的第k個(gè)諧振模式為sk,則det[Y(sk)]=0。由矩陣?yán)碚撝仃囆辛惺降闹档扔诰仃囁刑卣髦档某朔e,因此,復(fù)常數(shù)矩陣Y(sk)必有一個(gè)零特征值=0。設(shè)Y(sk)為n階

    圖4 電力網(wǎng)絡(luò)諧振模式的求解流程圖Fig.4 Solving process of electric network resonance modes

    Y(sk)=RΛR-1=RΛL

    (13)

    式中:Λ=diag(λ1,λ2,…,λn)是對(duì)角元素為特征值的對(duì)角矩陣,不妨設(shè)1=0;R=[R1,R2,…,Rn]是Y(sk)的右特征向量矩陣;L=R-1=T是Y(sk)的左特征向量矩陣。且由于Y(sk)為對(duì)稱矩陣,因此又有R-1=RT。根據(jù)s域節(jié)點(diǎn)導(dǎo)納方程式(4),有

    Y(sk)Vnode(sk)=Inode(sk)

    (14)

    Umode=LVnode(sk)

    (15)

    Jmode=LInode(sk)

    (16)

    則式(14)可以變換為

    (17)

    (18)

    (19)

    因此,我們定義R1為對(duì)應(yīng)諧振模式sk的節(jié)點(diǎn)電壓振型,表示在諧振模式sk下電網(wǎng)中各節(jié)點(diǎn)電壓的相對(duì)大小和相位。而根據(jù)式(16),有

    J1=L1Inode(sk)

    (20)

    因此,根據(jù)式(18)~(20)有

    (21)

    (22)

    因此定義矩陣

    (23)

    為參與因子矩陣,P的元素pij表示在諧振模式sk下電網(wǎng)中節(jié)點(diǎn)j的注入電流對(duì)節(jié)點(diǎn)i電壓的相對(duì)作用大小。

    4 算例及結(jié)果分析

    4.1 IEEE次同步諧振第一標(biāo)準(zhǔn)測(cè)試系統(tǒng)分析

    IEEE次同步諧振第一標(biāo)準(zhǔn)測(cè)試系統(tǒng)結(jié)構(gòu)如圖5所示,發(fā)電機(jī)參數(shù)見文獻(xiàn)[12]。發(fā)電機(jī)的阻抗頻率特性采用測(cè)試信號(hào)法[13]計(jì)算得出,如圖6所示。

    圖5 IEEE 次同步振蕩第一標(biāo)準(zhǔn)測(cè)試系統(tǒng)接線圖Fig.5 Single line diagram of IEEE subsynchronousresonance (SSR) first benchmark model

    對(duì)該測(cè)試系統(tǒng)在5~115 Hz頻率范圍進(jìn)行固有諧振結(jié)構(gòu)分析,發(fā)現(xiàn)存在1個(gè)次同步諧振模式,頻率fssr=39.3 Hz,衰減因子σssr=-0.454 51/s。由于衰減因子為負(fù),說明該測(cè)試系統(tǒng)是諧振不穩(wěn)定的。該諧振模式的節(jié)點(diǎn)電壓振型如圖7所示,參與因子矩陣如式(24)所示。由節(jié)點(diǎn)電壓振型圖和參與因子列表可見,模式SSR的主要參與節(jié)點(diǎn)是3號(hào)節(jié)點(diǎn),諧振類型是全局性的。

    圖6 發(fā)電機(jī)的阻抗頻率特性Fig.6 Impedance-frequency characteristicsof generator

    (24)

    圖7 模式SSR的節(jié)點(diǎn)電壓振型圖Fig.7 Nodal voltage mode shape of mode SSR

    4.2 IEEE 39節(jié)點(diǎn)測(cè)試系統(tǒng)分析

    IEEE 39節(jié)點(diǎn)測(cè)試系統(tǒng)[14]如圖8所示。對(duì)該測(cè)試系統(tǒng)在0~1 500 Hz頻率范圍進(jìn)行固有諧振結(jié)構(gòu)分析,結(jié)果如表1所示。表1顯示該系統(tǒng)在0~1 500 Hz頻率范圍存在19個(gè)諧振模式,所有諧振模式的衰減因子均大于0,說明系統(tǒng)是諧振穩(wěn)定的。下面對(duì)諧振模式1(131.6 Hz)進(jìn)行詳細(xì)分析。模式1 (131.6 Hz) 的節(jié)點(diǎn)電壓振型圖如圖9所示,數(shù)值大于0.1的參與因子矩陣元素只有1個(gè),為p28,28=0.104 1。

    由節(jié)點(diǎn)電壓振型圖和參與因子矩陣可見,模式1的主要參與節(jié)點(diǎn)是28號(hào)和26號(hào)、29號(hào)節(jié)點(diǎn),諧振類型為局部諧振,主要是28號(hào)節(jié)點(diǎn)的無功負(fù)荷和所連線路的對(duì)地電容之間的諧振。

    圖8 IEEE 39節(jié)點(diǎn)測(cè)試系統(tǒng)Fig.8 IEEE 39-bus test system

    4.3 加入風(fēng)電場(chǎng)的IEEE 39節(jié)點(diǎn)測(cè)試系統(tǒng)分析

    在IEEE 39節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)上增加1個(gè)風(fēng)電場(chǎng)節(jié)點(diǎn),考慮風(fēng)電場(chǎng)通過遠(yuǎn)距離輸電線路加裝50%串補(bǔ)接入到IEEE 39節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)的25號(hào)節(jié)點(diǎn)上,如圖10所示。

    由雙饋風(fēng)力發(fā)電機(jī)聚合而成的風(fēng)電場(chǎng)額定容量 1 500 MW,風(fēng)機(jī)阻抗模型采用文獻(xiàn)[15]給出的模型,聚合后風(fēng)電場(chǎng)在0~100 Hz范圍內(nèi)的阻抗頻率特性如圖11所示,顯示風(fēng)電場(chǎng)在25~45 Hz頻率范圍內(nèi)呈現(xiàn)負(fù)電阻-電感特性(相頻特性大于90°)。

    表1IEEE39節(jié)點(diǎn)測(cè)試系統(tǒng)的諧振模式
    Table1ResonancemodesofIEEE39-bustestsystem

    圖9 模式1(131.6 Hz)的節(jié)點(diǎn)電壓振型圖Fig.9 Nodal voltage mode shape of mode 1 (131.6 Hz)

    圖10 加入到IEEE 39節(jié)點(diǎn)標(biāo)準(zhǔn)測(cè)試系統(tǒng)上的風(fēng)電場(chǎng)Fig.10 Wind farm connected to IEEE 39-bus test system

    對(duì)該系統(tǒng)在0~1 500 Hz頻率范圍進(jìn)行固有諧振結(jié)構(gòu)分析,結(jié)果如表2所示。表2與表1對(duì)比顯示,雙饋風(fēng)機(jī)通過串補(bǔ)系統(tǒng)接入原系統(tǒng)25號(hào)節(jié)點(diǎn)后,部分諧振模式發(fā)生了變化,如新增了29.8 Hz的諧振模式,而且29.8 Hz的諧振模式衰減因子為負(fù)(發(fā)散),說明雙饋風(fēng)機(jī)通過串補(bǔ)送出可以引起不穩(wěn)定的次同步諧振。進(jìn)一步對(duì)不穩(wěn)定的29.8 Hz諧振模式進(jìn)行詳細(xì)分析,得到模式1(29.8 Hz)的節(jié)點(diǎn)電壓振

    圖11 風(fēng)機(jī)阻抗頻率特性Fig.10 Impedance-frequency characteristic ofwind generator

    圖12 模式1(29.8 Hz)的節(jié)點(diǎn)電壓振型圖Fig.12 Nodal voltage mode shape of mode 1 (29.8 Hz)

    5 結(jié) 論

    本文通過構(gòu)建電力網(wǎng)絡(luò)的s域節(jié)點(diǎn)導(dǎo)納矩陣,提出了電力網(wǎng)絡(luò)諧振穩(wěn)定性的概念,并提出了相應(yīng)的穩(wěn)定判據(jù)——s域節(jié)點(diǎn)導(dǎo)納矩陣行列式的零點(diǎn)就是系統(tǒng)的特征值,若零點(diǎn)均位于左半平面,則電力網(wǎng)絡(luò)是穩(wěn)定的。這一穩(wěn)定判據(jù)對(duì)含頻變?cè)头植紖?shù)元件的電力網(wǎng)絡(luò)具有很好的適用性。在求解零點(diǎn)的過程中,本文提出了一種實(shí)部-虛部交叉迭代法,這種求解方法與采用Newton-Raphson迭代求解的方法相比,可以大幅度減少計(jì)算量。本文還闡述了在固定諧振模式下電力網(wǎng)絡(luò)參與因子矩陣和節(jié)點(diǎn)電壓振型的意義。通過算例分析可以發(fā)現(xiàn):

    (1)元件的負(fù)電阻效應(yīng)是導(dǎo)致電力網(wǎng)絡(luò)諧振不穩(wěn)定的主要原因,同步發(fā)電機(jī)在次同步低頻段(低于基頻)由于異步發(fā)電機(jī)效應(yīng)會(huì)表現(xiàn)出負(fù)電阻;

    (2)雙饋風(fēng)機(jī)在次同步頻段也存在著負(fù)電阻效應(yīng),它的接入可能會(huì)引起電力網(wǎng)絡(luò)在次同步頻段的諧振不穩(wěn)定;

    (3)電力網(wǎng)絡(luò)的諧振模式通常表現(xiàn)為2種形式:一種是節(jié)點(diǎn)對(duì)地的諧振模式,主要特征為節(jié)點(diǎn)電壓的振型接近于同一方向;另一種是節(jié)點(diǎn)對(duì)節(jié)點(diǎn)的諧振模式,主要特征為節(jié)點(diǎn)電壓的振型會(huì)存在2個(gè)明顯相反的方向。

    [1] Kundur Prabha, Balu Neal J, Lauby Mark G. Power system stability and control[M]. McGraw-hill New York, 1994.

    [2] 徐政. 交直流電力系統(tǒng)動(dòng)態(tài)行為分析[M]. 北京: 機(jī)械工業(yè)出版社, 2005.

    [3] 呂敬, 董鵬, 施剛,等. 大型雙饋風(fēng)電場(chǎng)經(jīng)MMC-HVDC并網(wǎng)的次同步振蕩及其抑制[J]. 中國電機(jī)工程學(xué)報(bào), 2015, 35(19): 4852-4860.

    Lü Jing, DONG Peng, SHI Gang,et al. Subsynchronous Oscillation and Its Mitigation of MMC-Based HVDC With Large Doubly-Fed Induction Generator-Based Wind Farm Integration[J]. Proceedings of the CSEE, 2015, 35(19): 4852-4860.

    [4] 李云豐, 湯廣福, 賀之淵,等. MMC型直流輸電系統(tǒng)阻尼控制策略研究[J]. 中國電機(jī)工程學(xué)報(bào), 2016,36(20): 5492-5503.

    LI Yunfeng, TANG Guangfu, HE Zhiyuan,et al. Damping Control Strategy Research for MMC Based HVDC System[J]. Proceedings of the CSEE, 2016, 36(20): 5492-5503.

    [5] 李明節(jié), 于釗, 許濤,等. 新能源并網(wǎng)系統(tǒng)引發(fā)的復(fù)雜振蕩問題及其對(duì)策研究[J]. 電網(wǎng)技術(shù), 2017, 41(4): 1035-1042.

    LI Mingjie, YU Zhao, XU Tao,et al. Study of Complex Oscillation Caused by Renewable Energy Integration and Its Solution[J]. Power System Technology, 2017, 41(4): 1035-1042.

    [6] 王錫凡, 方萬良, 杜正春. 現(xiàn)代電力系統(tǒng)分析[M]. 北京:科學(xué)出版社, 2003.

    [7] Semlyen A I. s-domain methodology for assessing the small signal stability of complex systems in nonsinusoidal steady state[J]. IEEE Transactions on Power Systems, 1999, 14(1): 132-137.

    [8] Gomes S, Martins N, Portela C. Modal analysis applied to s-domain models of AC networks[C]//IEEE Power and Energy Society.Proceedings of 2001 IEEE Power Engineering Society Winter Meeting,2001.1.28—2001.2.1, Columbus, Ohio, USA: 1305-1310.

    [9] Lima L T G, Martins N, Carneiro S. Augmented state-space formulation for the study of electric networks including distributed-parameter transmission line models [C]// IPST commission.Proceedings of 1999 International conference on power systems transients,1999.6.20-24, Budapest, Hungary:87-92.

    [10] Varricchio S L, Martins N, Lima L T G. A Newton-Raphson method based on eigenvalue sensitivities to improve harmonic voltage performance[J]. IEEE Transactions on Power Delivery, 2003, 18(1): 334-342.

    [11] Xu W, Huang Z, Cui Y, et al. Harmonic resonance mode analysis[J]. IEEE Transactions on Power Delivery, 2005, 20(2): 1182-1190.

    [12] IEEE subsynchronous resonance task force. First benchmark model for computer simulation of subsynchronous resonance[J]. IEEE Trans on Power Apparatus and Systems, 1977,96(5):1565-1572.

    [13] Xu Zheng, Feng Zhouyan. A novel unified approach for analyzing small-signal stability of power systems[C]//IEEE Power and Energy Society.Proceedings of 2000 IEEE Power Engineering Society Winter Meeting,2000.1.23-27, Singapore:963-967.

    [14] Ramos R A, Hiskens I, Canizares C, et al. Benchmark systems for small-signal stability analysis and control[R]. IEEE Power and Energy Society Technical Report PES-TR18, 2015.

    [15] Miao Z. Impedance-Model-Based SSR Analysis for Type 3 Wind Generator and Series-Compensated Network[J]. IEEE Transactions on Energy Conversion, 2012, 27(4): 984-991.

    2017-09-10

    徐政(1962),男,博士,教授,通信作者,主要研究方向:大規(guī)模交直流電力系統(tǒng)分析、直流輸電與柔性交流輸電、風(fēng)力發(fā)電技術(shù)與風(fēng)電場(chǎng)并網(wǎng)技術(shù);

    王世佳(1991),男,博士研究生,主要研究方向:交直流電力系統(tǒng)次同步振蕩;

    邢法財(cái)(1993),男,博士研究生,主要研究方向:新能源并網(wǎng)的交直流系統(tǒng)穩(wěn)定性研究;

    肖晃慶(1990),男,博士研究生,主要研究方向:直流輸電與柔性交流輸電。

    (編輯 魏希輝)

    QualitativeAnalysisMethodofElectricNetworkResonanceStability

    XU Zheng, WANG Shijia, XING Facai, XIAO Huangqing

    (College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China)

    With the increasing utilization of power electronic equipment in power systems, in recent years,a number of new oscillations with unknown mechanisms have emerged. This paper puts forward the concept of the electric network resonance stability, and tries to classify the unclear reason oscillations mentioned above into the electric network resonance instability category. Thus, many complex power system oscillations can be analyzed by the linear network theory, which provides an approach to solve the actual power system problems by the mathematically mature linear system theory. The objective of this paper is to establish a method for analyzing the electric network resonance stability. By introducing the s-domain nodal admittance matrix of the electric network, this paper transforms the discrimination of the electric network resonance stability into the distribution problem in complex plane of zero point of the determinant of the s-domain nodal admittance matrix. Firstly, it is proved that the zero points of the determinant of the s-domain nodal admittance matrix are actually the eigenvalues of the system. Secondly, we use the cross iteration method of the real part and the imaginary part of the zero point to solve the zero points of the determinant of the s-domain nodal admittance matrix. Thirdly, we derive the nodal voltage mode shape and the participation factor matrix corresponding to a particular resonance mode, which can be used to locate the resonant region of this particular resonance mode in the network. Finally, we illustrate the effectiveness of the proposed method for analyzing the resonance stability of electric networks by several studied cases.

    electric network; state space model; s-domain nodal admittance matrix; resonance mode; resonance stability;nodal voltage mode shape; participation factor matrix

    國家高技術(shù)研究發(fā)展計(jì)劃項(xiàng)目(863計(jì)劃)(2011AA05A119);國家電網(wǎng)公司科技項(xiàng)目(柔性輸電網(wǎng)規(guī)劃評(píng)估方法及應(yīng)用關(guān)鍵技術(shù)研究)

    Project supported by The National High Technology Research and Development of China (863 Program)( 2011AA05A119)

    TM711

    A

    1000-7229(2017)11-0001-08

    10.3969/j.issn.1000-7229.2017.11.001

    猜你喜歡
    系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    WJ-700無人機(jī)系統(tǒng)
    ZC系列無人機(jī)遙感系統(tǒng)
    基于PowerPC+FPGA顯示系統(tǒng)
    基于UG的發(fā)射箱自動(dòng)化虛擬裝配系統(tǒng)開發(fā)
    半沸制皂系統(tǒng)(下)
    FAO系統(tǒng)特有功能分析及互聯(lián)互通探討
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    一德系統(tǒng) 德行天下
    PLC在多段調(diào)速系統(tǒng)中的應(yīng)用
    亚洲欧洲日产国产| 中文字幕久久专区| 成人午夜高清在线视频| 91aial.com中文字幕在线观看| 伦理电影大哥的女人| 美女被艹到高潮喷水动态| 免费黄网站久久成人精品| 嘟嘟电影网在线观看| 精品一区二区三区视频在线| 欧美在线一区亚洲| 哪里可以看免费的av片| 一进一出抽搐gif免费好疼| 久久久久国产网址| 欧美潮喷喷水| 国产乱人偷精品视频| 亚洲欧美清纯卡通| 校园人妻丝袜中文字幕| 国产精品一区二区三区四区久久| 欧美日韩在线观看h| 国产亚洲5aaaaa淫片| 精品人妻偷拍中文字幕| 亚洲第一区二区三区不卡| 成年版毛片免费区| 在线播放无遮挡| 天堂网av新在线| 高清毛片免费看| 亚洲精品乱码久久久v下载方式| 久久精品久久久久久久性| 日本欧美国产在线视频| 婷婷精品国产亚洲av| 桃色一区二区三区在线观看| 国产中年淑女户外野战色| 国产亚洲av嫩草精品影院| 午夜老司机福利剧场| 国产成人影院久久av| 国产精华一区二区三区| 一区二区三区四区激情视频 | 舔av片在线| 国产精品一区二区三区四区久久| 丰满的人妻完整版| 草草在线视频免费看| 亚洲精品亚洲一区二区| 你懂的网址亚洲精品在线观看 | 国产精品不卡视频一区二区| 亚洲在线自拍视频| 高清在线视频一区二区三区 | 搡老妇女老女人老熟妇| 亚洲一级一片aⅴ在线观看| 国产不卡一卡二| 在线观看av片永久免费下载| 91久久精品电影网| 在线天堂最新版资源| 日韩国内少妇激情av| 天堂网av新在线| 久久久久久久久久久丰满| 成年免费大片在线观看| 丝袜喷水一区| 欧美成人a在线观看| 午夜爱爱视频在线播放| 成人无遮挡网站| av女优亚洲男人天堂| 亚洲国产色片| 校园春色视频在线观看| 一本精品99久久精品77| 大又大粗又爽又黄少妇毛片口| 中文字幕免费在线视频6| 非洲黑人性xxxx精品又粗又长| 成人无遮挡网站| 亚洲va在线va天堂va国产| 色视频www国产| a级毛片免费高清观看在线播放| 一区福利在线观看| 国产爱豆传媒在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲久久久久久中文字幕| 99热全是精品| 麻豆国产av国片精品| 欧美一区二区亚洲| 亚洲熟妇中文字幕五十中出| 亚洲一级一片aⅴ在线观看| 18禁黄网站禁片免费观看直播| 精品一区二区三区人妻视频| 欧美+日韩+精品| av专区在线播放| 国产精品永久免费网站| 干丝袜人妻中文字幕| 国产成人aa在线观看| 亚洲精品久久国产高清桃花| 欧美性感艳星| 久久久久久国产a免费观看| 99热只有精品国产| 久久久久久久久大av| 日韩三级伦理在线观看| 日本爱情动作片www.在线观看| 日韩高清综合在线| 99国产精品一区二区蜜桃av| 18禁在线播放成人免费| 淫秽高清视频在线观看| 久久久久久国产a免费观看| 亚洲av成人av| 色吧在线观看| 只有这里有精品99| 欧美三级亚洲精品| 在线播放无遮挡| 看黄色毛片网站| 国产探花极品一区二区| 国产精品不卡视频一区二区| 黄色配什么色好看| 免费av不卡在线播放| 色综合色国产| 天堂中文最新版在线下载 | 亚洲色图av天堂| 亚洲av中文av极速乱| 少妇丰满av| 爱豆传媒免费全集在线观看| АⅤ资源中文在线天堂| av在线天堂中文字幕| 插阴视频在线观看视频| 欧美色视频一区免费| 国产精品嫩草影院av在线观看| 99热这里只有是精品在线观看| 免费看av在线观看网站| 国产亚洲欧美98| 熟女人妻精品中文字幕| 国产69精品久久久久777片| 欧美高清性xxxxhd video| 久久久久久久久久久丰满| 亚洲第一电影网av| 十八禁国产超污无遮挡网站| 又黄又爽又刺激的免费视频.| 尤物成人国产欧美一区二区三区| 国产 一区 欧美 日韩| 免费观看精品视频网站| 日本撒尿小便嘘嘘汇集6| 日日干狠狠操夜夜爽| 99在线视频只有这里精品首页| 久久精品久久久久久久性| 丝袜喷水一区| 大又大粗又爽又黄少妇毛片口| 欧美高清成人免费视频www| 国产亚洲精品av在线| 91狼人影院| 中国美白少妇内射xxxbb| 少妇熟女aⅴ在线视频| 欧美激情久久久久久爽电影| 亚洲av中文av极速乱| 国产精品免费一区二区三区在线| 听说在线观看完整版免费高清| 深爱激情五月婷婷| 最近视频中文字幕2019在线8| 亚洲欧美清纯卡通| 亚洲精品乱码久久久v下载方式| 日韩欧美在线乱码| 免费无遮挡裸体视频| 在线a可以看的网站| 国产精品国产高清国产av| 婷婷亚洲欧美| 亚洲国产高清在线一区二区三| av在线蜜桃| 日韩人妻高清精品专区| 国产三级中文精品| 久久国产乱子免费精品| 成人性生交大片免费视频hd| 一本久久精品| 久99久视频精品免费| 99热全是精品| 国产精品伦人一区二区| 九九爱精品视频在线观看| 亚洲天堂国产精品一区在线| 亚洲在线自拍视频| 看片在线看免费视频| 欧美日韩国产亚洲二区| 黄色一级大片看看| 国语自产精品视频在线第100页| 色综合色国产| 插逼视频在线观看| 国产女主播在线喷水免费视频网站 | 久久久久久久久久久免费av| 欧美日韩综合久久久久久| 欧美日韩综合久久久久久| 免费观看在线日韩| 欧美一区二区精品小视频在线| 丰满的人妻完整版| 18禁在线无遮挡免费观看视频| 内地一区二区视频在线| 亚洲经典国产精华液单| 色吧在线观看| av女优亚洲男人天堂| 中文字幕av在线有码专区| 99国产极品粉嫩在线观看| 日韩成人av中文字幕在线观看| 成人鲁丝片一二三区免费| 少妇高潮的动态图| 亚洲av男天堂| 亚洲av免费高清在线观看| 免费观看精品视频网站| 神马国产精品三级电影在线观看| 中国国产av一级| 免费看日本二区| 欧美不卡视频在线免费观看| www.色视频.com| 国产片特级美女逼逼视频| 欧美另类亚洲清纯唯美| 黄色一级大片看看| 国产成年人精品一区二区| 好男人在线观看高清免费视频| 18+在线观看网站| 哪个播放器可以免费观看大片| 亚洲婷婷狠狠爱综合网| 老女人水多毛片| 色哟哟哟哟哟哟| 国产一区二区激情短视频| 99在线视频只有这里精品首页| 欧美最黄视频在线播放免费| 91午夜精品亚洲一区二区三区| 长腿黑丝高跟| 欧美色视频一区免费| 国产单亲对白刺激| 成人美女网站在线观看视频| 午夜免费男女啪啪视频观看| 人妻制服诱惑在线中文字幕| 国产一级毛片七仙女欲春2| 自拍偷自拍亚洲精品老妇| 一区福利在线观看| 99在线人妻在线中文字幕| 欧美在线一区亚洲| 99在线视频只有这里精品首页| 亚洲一区高清亚洲精品| 精品久久久久久久久久久久久| 国产精品野战在线观看| 大香蕉久久网| 欧美三级亚洲精品| 亚洲真实伦在线观看| 婷婷六月久久综合丁香| 免费大片18禁| 能在线免费看毛片的网站| 热99re8久久精品国产| 亚洲人成网站在线播| 成人美女网站在线观看视频| 成人永久免费在线观看视频| 成年女人看的毛片在线观看| 亚洲精品日韩在线中文字幕 | 国产精品一区www在线观看| 99精品在免费线老司机午夜| 精品少妇黑人巨大在线播放 | 国产精品乱码一区二三区的特点| 在线观看66精品国产| 男人舔奶头视频| 久久人人精品亚洲av| 久久中文看片网| 噜噜噜噜噜久久久久久91| 美女内射精品一级片tv| 热99re8久久精品国产| av黄色大香蕉| 欧美日韩国产亚洲二区| 日本黄色视频三级网站网址| 国产三级在线视频| 丰满的人妻完整版| 久久综合国产亚洲精品| 久久国内精品自在自线图片| 一级av片app| 久久久久久久亚洲中文字幕| 伊人久久精品亚洲午夜| 寂寞人妻少妇视频99o| 哪个播放器可以免费观看大片| 中文精品一卡2卡3卡4更新| 日本爱情动作片www.在线观看| 国产淫片久久久久久久久| 久久99热6这里只有精品| av黄色大香蕉| 搞女人的毛片| 中文亚洲av片在线观看爽| 色吧在线观看| 精品无人区乱码1区二区| 亚洲国产精品合色在线| 久久中文看片网| 秋霞在线观看毛片| 欧美最新免费一区二区三区| 波多野结衣高清作品| 久久精品国产自在天天线| 国产精品免费一区二区三区在线| 午夜爱爱视频在线播放| 不卡视频在线观看欧美| 亚洲欧洲国产日韩| 老司机影院成人| 色5月婷婷丁香| 国产精品av视频在线免费观看| 久久久国产成人免费| 亚洲欧美日韩东京热| 久久久久久久久大av| 免费在线观看成人毛片| 亚洲色图av天堂| 国产三级在线视频| 长腿黑丝高跟| 国产精华一区二区三区| 国产乱人偷精品视频| 欧美不卡视频在线免费观看| 国产午夜精品久久久久久一区二区三区| 久久人人爽人人爽人人片va| 中国美女看黄片| 精品99又大又爽又粗少妇毛片| 少妇熟女aⅴ在线视频| 村上凉子中文字幕在线| 日韩强制内射视频| 成人av在线播放网站| 精品久久久久久久久久免费视频| 波野结衣二区三区在线| 国产精品人妻久久久影院| av女优亚洲男人天堂| 中文字幕人妻熟人妻熟丝袜美| 婷婷亚洲欧美| 村上凉子中文字幕在线| 亚洲最大成人av| 老司机影院成人| 99热精品在线国产| 麻豆久久精品国产亚洲av| 我的女老师完整版在线观看| 国产精品一二三区在线看| 成人高潮视频无遮挡免费网站| 丰满的人妻完整版| 乱系列少妇在线播放| 啦啦啦韩国在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 色吧在线观看| 欧美成人一区二区免费高清观看| 一卡2卡三卡四卡精品乱码亚洲| 免费看光身美女| 黄片无遮挡物在线观看| 免费不卡的大黄色大毛片视频在线观看 | 日本免费一区二区三区高清不卡| 深夜精品福利| 卡戴珊不雅视频在线播放| 夜夜夜夜夜久久久久| 国产伦理片在线播放av一区 | 哪里可以看免费的av片| 日韩欧美三级三区| 欧美性猛交╳xxx乱大交人| 在线免费观看不下载黄p国产| av专区在线播放| 在线观看av片永久免费下载| 天美传媒精品一区二区| 中文字幕久久专区| 国产午夜精品久久久久久一区二区三区| av卡一久久| 国产精品电影一区二区三区| 日本与韩国留学比较| 免费人成视频x8x8入口观看| 别揉我奶头 嗯啊视频| 久久久久久九九精品二区国产| 淫秽高清视频在线观看| 免费在线观看成人毛片| 伊人久久精品亚洲午夜| 久久久久久久久大av| 麻豆国产av国片精品| 欧美日韩国产亚洲二区| 丰满的人妻完整版| 久久久国产成人免费| 在线播放无遮挡| 国产不卡一卡二| 美女cb高潮喷水在线观看| 成人欧美大片| 日本五十路高清| 岛国在线免费视频观看| 在线播放国产精品三级| 欧美日韩一区二区视频在线观看视频在线 | 一卡2卡三卡四卡精品乱码亚洲| h日本视频在线播放| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲成人中文字幕在线播放| 狠狠狠狠99中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 精品熟女少妇av免费看| 亚洲内射少妇av| 精品无人区乱码1区二区| 看免费成人av毛片| 国产精品,欧美在线| 麻豆成人午夜福利视频| 国产午夜精品论理片| 亚洲欧美成人综合另类久久久 | 69人妻影院| 午夜福利成人在线免费观看| 一级黄色大片毛片| 国产黄a三级三级三级人| 校园春色视频在线观看| 大型黄色视频在线免费观看| 深爱激情五月婷婷| 日韩强制内射视频| 日韩欧美国产在线观看| 亚洲精品日韩av片在线观看| 99久久精品国产国产毛片| 天天一区二区日本电影三级| 我要看日韩黄色一级片| 国产av麻豆久久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产精品乱码一区二三区的特点| 国产精品爽爽va在线观看网站| 亚洲在线自拍视频| 九九久久精品国产亚洲av麻豆| 国产人妻一区二区三区在| 看黄色毛片网站| 久久6这里有精品| 日本黄色片子视频| 国产精品.久久久| 精品午夜福利在线看| 久久人人爽人人片av| 国产白丝娇喘喷水9色精品| 久久久久久国产a免费观看| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| 美女cb高潮喷水在线观看| 黄色一级大片看看| 少妇裸体淫交视频免费看高清| 看黄色毛片网站| 九色成人免费人妻av| 丝袜喷水一区| 久久亚洲精品不卡| 青春草国产在线视频 | 亚洲乱码一区二区免费版| 高清毛片免费观看视频网站| 国产精品麻豆人妻色哟哟久久 | 久久久久网色| 波多野结衣巨乳人妻| 少妇人妻一区二区三区视频| 欧美在线一区亚洲| 老女人水多毛片| 又粗又爽又猛毛片免费看| 深爱激情五月婷婷| 蜜桃久久精品国产亚洲av| 麻豆久久精品国产亚洲av| 一区二区三区四区激情视频 | 国产蜜桃级精品一区二区三区| 国产真实乱freesex| 看免费成人av毛片| 亚洲av成人av| 日韩一区二区三区影片| 国产毛片a区久久久久| 国产成年人精品一区二区| 欧美xxxx性猛交bbbb| 久久热精品热| 十八禁国产超污无遮挡网站| 日日啪夜夜撸| 毛片一级片免费看久久久久| 欧美潮喷喷水| 最新中文字幕久久久久| 美女cb高潮喷水在线观看| 日日摸夜夜添夜夜添av毛片| 亚洲av.av天堂| 日本撒尿小便嘘嘘汇集6| 全区人妻精品视频| av在线天堂中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美 国产精品| 亚洲精品影视一区二区三区av| av专区在线播放| 亚洲国产精品sss在线观看| 久久午夜亚洲精品久久| 日韩国内少妇激情av| 99热只有精品国产| 久久久国产成人精品二区| 黄色日韩在线| 日韩强制内射视频| 夜夜夜夜夜久久久久| 三级国产精品欧美在线观看| 少妇裸体淫交视频免费看高清| 亚洲人成网站高清观看| 日韩欧美精品v在线| 能在线免费观看的黄片| 亚洲色图av天堂| 色综合站精品国产| 亚洲三级黄色毛片| 欧美bdsm另类| 干丝袜人妻中文字幕| 99久久久亚洲精品蜜臀av| 狠狠狠狠99中文字幕| 成人欧美大片| 在线免费观看的www视频| 亚洲精品456在线播放app| 国产成人精品久久久久久| 国产黄片视频在线免费观看| 天堂av国产一区二区熟女人妻| 女同久久另类99精品国产91| 亚洲欧美日韩卡通动漫| 赤兔流量卡办理| 亚洲av成人av| 中文字幕免费在线视频6| 久久久午夜欧美精品| 男插女下体视频免费在线播放| 亚洲无线观看免费| 亚洲国产欧洲综合997久久,| 国产91av在线免费观看| 神马国产精品三级电影在线观看| 男插女下体视频免费在线播放| 伊人久久精品亚洲午夜| 久久久精品94久久精品| 国产又黄又爽又无遮挡在线| 3wmmmm亚洲av在线观看| 99热全是精品| 欧美+亚洲+日韩+国产| 亚洲四区av| 国产亚洲精品久久久久久毛片| 青青草视频在线视频观看| 成人综合一区亚洲| av又黄又爽大尺度在线免费看 | 日韩欧美精品v在线| 夜夜夜夜夜久久久久| 国产极品天堂在线| 狂野欧美白嫩少妇大欣赏| 国内精品美女久久久久久| 亚洲五月天丁香| 一个人看视频在线观看www免费| 床上黄色一级片| 国产在视频线在精品| 亚洲国产精品sss在线观看| 国产中年淑女户外野战色| 午夜福利在线观看免费完整高清在 | 好男人视频免费观看在线| 国产午夜精品论理片| 床上黄色一级片| 午夜久久久久精精品| 精华霜和精华液先用哪个| 高清在线视频一区二区三区 | 国产午夜精品一二区理论片| 亚洲一级一片aⅴ在线观看| 嘟嘟电影网在线观看| 亚洲精品日韩在线中文字幕 | 白带黄色成豆腐渣| 欧美性猛交╳xxx乱大交人| 精品国内亚洲2022精品成人| 久久久成人免费电影| 99久国产av精品| 日日啪夜夜撸| 老司机影院成人| 欧美日韩综合久久久久久| 男插女下体视频免费在线播放| 免费观看a级毛片全部| 午夜老司机福利剧场| 欧美不卡视频在线免费观看| 嫩草影院精品99| 免费人成在线观看视频色| 男女下面进入的视频免费午夜| 99热这里只有是精品在线观看| 我的老师免费观看完整版| 给我免费播放毛片高清在线观看| 成人高潮视频无遮挡免费网站| 色哟哟·www| 黄色欧美视频在线观看| 国产精品,欧美在线| 蜜桃久久精品国产亚洲av| 3wmmmm亚洲av在线观看| 日本黄色片子视频| 国内少妇人妻偷人精品xxx网站| 99久久精品国产国产毛片| av天堂中文字幕网| 日日摸夜夜添夜夜爱| 欧美日韩国产亚洲二区| 日韩强制内射视频| 99久久精品一区二区三区| 国产精品一区二区三区四区久久| 中文字幕人妻熟人妻熟丝袜美| 99热网站在线观看| 亚洲不卡免费看| 国产精品久久久久久亚洲av鲁大| 国产精品美女特级片免费视频播放器| av.在线天堂| 18禁在线无遮挡免费观看视频| 日韩,欧美,国产一区二区三区 | www日本黄色视频网| 亚洲一区高清亚洲精品| av在线观看视频网站免费| 99热6这里只有精品| 日韩av不卡免费在线播放| 在线免费十八禁| 亚洲成人久久性| 麻豆国产av国片精品| 亚洲一区二区三区色噜噜| 免费观看的影片在线观看| 亚洲人成网站在线播放欧美日韩| 18禁在线无遮挡免费观看视频| 日本成人三级电影网站| 亚洲av免费高清在线观看| 日本-黄色视频高清免费观看| 国产老妇伦熟女老妇高清| 成人毛片60女人毛片免费| 久久久久网色| 久久久久久九九精品二区国产| 一进一出抽搐动态| 2022亚洲国产成人精品| 国产69精品久久久久777片| 亚洲自拍偷在线| 成人二区视频| 精品久久久久久久久亚洲| 精品99又大又爽又粗少妇毛片| 一区福利在线观看| 欧美高清性xxxxhd video| 我的女老师完整版在线观看| 最近最新中文字幕大全电影3| 在线免费观看的www视频| 成人美女网站在线观看视频| 亚洲18禁久久av| 欧美潮喷喷水| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲av涩爱 | 赤兔流量卡办理| 偷拍熟女少妇极品色| 伊人久久精品亚洲午夜| 丝袜美腿在线中文| 亚洲人成网站在线播放欧美日韩| 日韩亚洲欧美综合| 国产精品久久久久久亚洲av鲁大| 干丝袜人妻中文字幕| 看非洲黑人一级黄片| 人妻系列 视频| 日韩一本色道免费dvd| 亚洲成人久久爱视频| 国产精品福利在线免费观看|