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

    基于軌跡靈敏度分析的直流調(diào)制參數(shù)優(yōu)化設(shè)計

    2012-11-09 10:06:57王長香房大中薛振宇魏文輝
    關(guān)鍵詞:阻尼靈敏度直流

    王長香,房大中,薛振宇,魏文輝

    (1.天津大學(xué)電氣與自動化工程學(xué)院,天津 300072;2.中國電力科學(xué)研究院,北京 100192)

    基于軌跡靈敏度分析的直流調(diào)制參數(shù)優(yōu)化設(shè)計

    王長香1,房大中1,薛振宇1,魏文輝2

    (1.天津大學(xué)電氣與自動化工程學(xué)院,天津 300072;2.中國電力科學(xué)研究院,北京 100192)

    提出一種基于時域仿真的直流調(diào)制參數(shù)優(yōu)化方法,將直流調(diào)制參數(shù)優(yōu)化問題轉(zhuǎn)化為帶參數(shù)約束的非線性優(yōu)化問題。該方法,首先將約束優(yōu)化問題轉(zhuǎn)化為無約束優(yōu)化問題,通過軌跡靈敏度仿真得到目標(biāo)函數(shù)的梯度信息,應(yīng)用擬牛頓法(DFP)優(yōu)化直流調(diào)制控制器參數(shù)。由于目標(biāo)函數(shù)反映了AC/DC復(fù)雜系統(tǒng)的動力學(xué)特性,該方法可有效地改善系統(tǒng)的阻尼特性,抑制大擾動引起的系統(tǒng)振蕩。在IEEE四機系統(tǒng)上進(jìn)行的特征值分析和時域仿真結(jié)果驗證了該方法的有效性和魯棒性。

    直流調(diào)制;軌跡靈敏度;擬牛頓法;參數(shù)優(yōu)化

    直流輸電功率調(diào)節(jié)快速靈活,不像交流系統(tǒng)存在同步穩(wěn)定性問題,是較理想的大容量遠(yuǎn)距離輸電方式。

    文獻(xiàn)[1]中指出直流系統(tǒng)常規(guī)控制方式雖然可以緩沖一側(cè)交流系統(tǒng)對另一側(cè)交流系統(tǒng)的擾動,但由于換流器像是一個對頻率不敏感的負(fù)荷,不但不會為發(fā)電機提供同步轉(zhuǎn)矩,反而會為發(fā)電機提供負(fù)的阻尼轉(zhuǎn)矩,不利于系統(tǒng)的穩(wěn)定。因此,為了提高交流系統(tǒng)的動態(tài)穩(wěn)定性,常須裝附加阻尼控制來拓展直流系統(tǒng)的可控性。

    為了提高交流系統(tǒng)的機電振蕩阻尼,抑制大擾動引起的系統(tǒng)振蕩,可以在直流系統(tǒng)控制器上附加阻尼控制器,即直流調(diào)制[2,3]。采用直流調(diào)制的一個典型例子是美國太平洋聯(lián)絡(luò)線[4],該工程采用小信號電流調(diào)制來抑制并聯(lián)交流線路上的功率振蕩。文獻(xiàn)[5]采用整流側(cè)的直流功率調(diào)制和逆變側(cè)的熄弧角調(diào)制來阻尼并聯(lián)交流線路上的區(qū)間低頻振蕩,取得了很好的控制效果。文獻(xiàn)[6,7]對直流調(diào)制改善南方電網(wǎng)交直流并聯(lián)系統(tǒng)穩(wěn)定性進(jìn)行了研究,研究表明直流調(diào)制可有效改善交直流并聯(lián)系統(tǒng)的穩(wěn)定性,快速抑制交流系統(tǒng)功率振蕩。但是,文獻(xiàn)[5,6,7]中采用的直流調(diào)制參數(shù)都是來源于工程經(jīng)驗或大量仿真計算確定,因為實際交直流并聯(lián)系統(tǒng)具有結(jié)構(gòu)復(fù)雜、動態(tài)特性很難精確描述的特點。只有合理配置直流調(diào)制參數(shù),才能有效改善系統(tǒng)阻尼,抑制系統(tǒng)功率振蕩,因此大規(guī)模AC/DC復(fù)雜系統(tǒng)的直流調(diào)制參數(shù)整定仍然是一個重要的研究課題。

    長期以來,學(xué)者們對直流調(diào)制控制器參數(shù)優(yōu)化問題進(jìn)行了大量的研究,并提出了多種優(yōu)化方法,這些方法已在實際應(yīng)用中取得了一定的效果[8~12]。近年來,遺傳算法、模擬退化算法、混沌優(yōu)化算法等智能優(yōu)化算法已經(jīng)廣泛應(yīng)用到參數(shù)優(yōu)化中,文獻(xiàn)[9]和[11]提出了基于混沌優(yōu)化算法的直流調(diào)制控制器參數(shù)優(yōu)化算法,并通過中國南方電網(wǎng)和三區(qū)域雙饋入交直流系統(tǒng)進(jìn)行了驗證。然而,智能優(yōu)化算法普遍具有隨機性強、收斂速度慢的缺點。文獻(xiàn)[8]和[10]對南方交直流混合電網(wǎng)直流調(diào)制的參數(shù)協(xié)調(diào)優(yōu)化問題進(jìn)行了研究,也取得了很好的效果。文獻(xiàn)[13]提出了一種基于軌跡靈敏度優(yōu)化電力系統(tǒng)穩(wěn)定器參數(shù)的方法,并通過新英格蘭系統(tǒng)和四機系統(tǒng)進(jìn)行驗證。

    將文獻(xiàn)[13]的思想應(yīng)用到直流調(diào)制控制器參數(shù)優(yōu)化中,提出了一種基于軌跡靈敏度分析的直流調(diào)制控制器參數(shù)優(yōu)化方法。該方法通過時域仿真得到系統(tǒng)的軌跡和軌跡靈敏度,繼而得到目標(biāo)函數(shù)的梯度信息,應(yīng)用擬牛頓法(DFP)優(yōu)化直流調(diào)制控制器參數(shù)。由于該方法目標(biāo)函數(shù)反映了AC/DC復(fù)雜系統(tǒng)的動力學(xué)特性,該方法設(shè)計的直流調(diào)制控制器,不僅可提高小擾動穩(wěn)定性,還可有效抑制大擾動后的系統(tǒng)振蕩。在IEEE四機系統(tǒng)上進(jìn)行的特征值分析和時域仿真結(jié)果驗證了該方法的有效性和魯棒性。

    1 AC/DC系統(tǒng)數(shù)學(xué)模型

    1.1AC/DC系統(tǒng)機電暫態(tài)仿真的基本原理

    為了更準(zhǔn)確地反映直流系統(tǒng)中元件的動態(tài)性能對交流系統(tǒng)穩(wěn)定性的影響,在AC/DC系統(tǒng)機電暫態(tài)仿真中,使用較為精確的直流輸電模型,即準(zhǔn)穩(wěn)態(tài)模型[14]。在準(zhǔn)穩(wěn)態(tài)模型中,換流器方程仍用直流分量有效值與基頻交流分量之間的關(guān)系表示,直流線路用R-L線路模擬。

    在AC/DC系統(tǒng)機電暫態(tài)仿真中,交流系統(tǒng)和直流系統(tǒng)是相互作用的,如圖1所示。對直流線路來說,兩側(cè)的換流器相當(dāng)于兩個電壓源,其數(shù)值決定于觸發(fā)控制角和交流側(cè)換流母線電壓有效值。對交流系統(tǒng)來說,直流系統(tǒng)相當(dāng)于一個變化的負(fù)荷,可用掛在換流站交流母線上的可變功率支路來模擬。

    圖1 AC/DC系統(tǒng)相互作用示意圖

    1.2 直流控制器模型及直流調(diào)制模型

    對于直流控制系統(tǒng),本文采用整流側(cè)定電流控制,逆變側(cè)定熄弧角控制的調(diào)節(jié)運行方式,直流控制器采用文獻(xiàn)[14]中給出的控制器模型,定電流控制器及直流調(diào)制模型如圖2。

    圖2 直流控制器及直流調(diào)制模型

    直流調(diào)制以與直流并聯(lián)交流聯(lián)絡(luò)線傳輸?shù)挠泄β蔖ac作為輸入信號,經(jīng)過測量、隔直和超前滯后環(huán)節(jié),再經(jīng)過限幅環(huán)節(jié)后得到電流調(diào)制信號Imod,Imod和電流控制信號Idref進(jìn)行綜合,經(jīng)過定電流控制器控制整流側(cè)觸發(fā)角α,從而控制直流線路傳輸功率,達(dá)到改善系統(tǒng)穩(wěn)定性的目的。

    圖2中,對于定電流控制器,Id和Idref分別為輸入信號和電流控制信號;Tc3為測量環(huán)節(jié)時間常數(shù);Kc2和Tc2分別為積分環(huán)節(jié)增益系數(shù)和時間常數(shù);Kc1為比例環(huán)節(jié)增益系數(shù);αref為整流側(cè)觸發(fā)角參考值。對于直流調(diào)制控制器,Ks為增益系數(shù);T7和Tw分別為測量環(huán)節(jié)和隔直環(huán)節(jié)時間常數(shù);T1、T2、T3、T4、T5、T6為超前滯后環(huán)節(jié)的時間常數(shù)。下標(biāo)max和min分別表示相應(yīng)變量的上下限值。圖中除角度和時間外,其它所有變量和參數(shù)均為標(biāo)幺值。

    1.3 電力系統(tǒng)暫態(tài)模型

    AC/DC混聯(lián)電力系統(tǒng)的非線性模型由發(fā)電機及其控制器、負(fù)荷和HVDC及其他裝置模型組成,并通過網(wǎng)絡(luò)方程連接,其動態(tài)過程可用如下一般形式的微分-代數(shù)方程組表示:

    (1)

    式中:微分方程組f(·)描述各動態(tài)元件的動態(tài)特性,包括發(fā)電機及其控制系統(tǒng)、直流線路和直流控制器等動態(tài)元件;代數(shù)方程組g(·)包括網(wǎng)絡(luò)方程、各動態(tài)元件與網(wǎng)絡(luò)接口方程、HVDC電壓方程等;向量x表示狀態(tài)變量;向量y表示運行參量;α表示待優(yōu)化參數(shù)向量;x0、y0表示狀態(tài)變量和運行參量的初值,由潮流計算確定。

    2 直流調(diào)制參數(shù)優(yōu)化設(shè)計

    2.1 直流調(diào)制參數(shù)優(yōu)化問題模型

    優(yōu)化直流調(diào)制控制器參數(shù)的目的是改善系統(tǒng)的阻尼特性,抑制系統(tǒng)區(qū)域間的低頻振蕩,使得電力系統(tǒng)以最快的速度恢復(fù)到穩(wěn)定運行點。本文采用系統(tǒng)中發(fā)電機的速度偏差構(gòu)成二次型性能指標(biāo)如下:

    (2)

    式(2)所表示性能指標(biāo)的物理意義如圖3所示,顯而易見,系統(tǒng)發(fā)電機速度偏差平方的積分為圖中陰影部分的面積,面積越小則系統(tǒng)振蕩平息的越快,也就是說系統(tǒng)的阻尼特性越好。因此,最小化性能指標(biāo)(2),系統(tǒng)振蕩的阻尼被最大化,此時對應(yīng)的參數(shù)為直流調(diào)制控制器的最優(yōu)參數(shù)。

    圖3 目標(biāo)函數(shù)的物理意義

    直流調(diào)制參數(shù)優(yōu)化問題模型式(2)是一個有約束非線性優(yōu)化問題,其一般解決方法是先轉(zhuǎn)化為無約束優(yōu)化問題。引入非線性變換[13,15]:

    (3)

    minJ(X)

    (4)

    2.2 軌跡靈敏度分析

    運用軌跡靈敏度方法可定量分析電力系統(tǒng)某些參數(shù)對其動態(tài)性能的影響,在電力系統(tǒng)參數(shù)辨識、電力系統(tǒng)穩(wěn)定分析與控制以及附加阻尼控制器的參數(shù)設(shè)置等領(lǐng)域都得到了廣泛的應(yīng)用。本文僅對軌跡靈敏度作簡單介紹,具體參見文獻(xiàn)[15]。對式(1)求偏導(dǎo)數(shù)可得:

    (5)

    式中:xα(t)和yα(t)分別表示狀態(tài)變量x和運行參量y對直流調(diào)制控制器參數(shù)的靈敏度;fx、fy和fα分別表示f對向量x、y和α的梯度;gx、gy和gα具有相同意義。本文稱式(1)為軌跡系統(tǒng),其伴隨系統(tǒng)式(5)被稱為軌跡靈敏度系統(tǒng)。運用雙向模塊簡化技術(shù)聯(lián)合求解軌跡和軌跡靈敏度系統(tǒng),可大大減少軌跡靈敏度系統(tǒng)的計算量[8]。

    采用擬牛頓法求解式(4)的無約束優(yōu)化問題的前提是要獲得目標(biāo)函數(shù)的梯度信息。通過時域仿真計算得到每一時步的ωi和ω(rj)i,將其帶入式(6),從而可得到目標(biāo)函數(shù)的梯度信息?J/?X=[?J/?β1,?J/?β2,…,?J/?βm]T。目標(biāo)函數(shù)對任意變量βj的導(dǎo)數(shù)表達(dá)式如下:

    (6)

    式中:ω(γj)i表示第i臺發(fā)電機相對同步坐標(biāo)的速度偏差對第j個待優(yōu)化參數(shù)γj的軌跡靈敏度;其它變量的含義,上文已作介紹,這里不再贅述。

    2.3 優(yōu)化算法簡介

    直流調(diào)制控制器參數(shù)優(yōu)化過程:首先給定一組初始參數(shù)及收斂條件,進(jìn)行時域仿真得到系統(tǒng)的軌跡和軌跡靈敏度,利用擬牛頓法優(yōu)化直流調(diào)制參數(shù),直到滿足收斂條件,得到最優(yōu)參數(shù)為止。采用VS2005編寫程序,實現(xiàn)本文所提算法。圖4為程序基本流程圖。

    圖4 直流調(diào)制參數(shù)優(yōu)化流程圖

    本文采用擬牛頓法(DFP)[16]求解直流調(diào)制控制器參數(shù)優(yōu)化問題。擬牛頓法的基本思想是用一階導(dǎo)數(shù)近似得到海森矩陣的逆矩陣,不必計算二階導(dǎo)數(shù),減少了計算量。同時綜合了最速下降法具有較好的整體收斂性和牛頓法收斂速度快的優(yōu)點。DFP算法是第一個被提出的擬牛頓算法,也是無約束最優(yōu)化問題最有效的算法之一。

    DFP算法步驟:

    步驟1給定控制誤差ε,給定直流調(diào)制控制器的初始參數(shù)α0,通過式(3)的非線性變換得到初始點X0,給定初始矩陣H0(通常取單位矩陣E),計算初始梯度g0,令k=0。

    (7)

    步驟2令k步的搜索方向pk。

    pk= -Hkgk

    (8)

    步驟3由精確一維搜索確定步長λk。

    (9)

    步驟4計算Xk+1,由式(3)計算得到αk+1。

    Xk+1=Xk+λkpk

    (10)

    步驟5將αk+1帶入軌跡靈敏度系統(tǒng),計算k+1步的梯度gk+1,若滿足‖gk+1‖≤ε,則最優(yōu)解X*=Xk+1;否則令

    sk=Xk+1-Xk,yk=gk+1-gk

    (11)

    步驟6由DEP修正公式計算Hk+1。令k=k+1,轉(zhuǎn)到步驟2。

    (12)

    3 算例分析

    圖5為IEEE四機系統(tǒng)的結(jié)構(gòu)圖,系統(tǒng)參數(shù)見文獻(xiàn)[1]。算例中:系統(tǒng)功率基準(zhǔn)值為100 MVA;發(fā)電機均采用雙軸暫態(tài)模型,每臺發(fā)電機配有EA模型的勵磁器;負(fù)荷采用恒阻抗負(fù)荷;兩端直流輸電系統(tǒng)采用準(zhǔn)穩(wěn)態(tài)模型,控制系統(tǒng)采用整流側(cè)定電流,逆變側(cè)定熄弧角調(diào)節(jié)運行方式,直流控制器及直流調(diào)制模型見圖2。直流調(diào)制輸出信號上下限Imod max和Imod min分別為0.3Ids(p.u.)和-0.2Ids(p.u.),其中Ids為給定直流線路電流標(biāo)幺值。定電流控制器參數(shù)取值為:Tc3=0.09,Kc1=0.75,Kc2=10,Tc2=1。對于直流調(diào)制控制器[1],Tw一般取經(jīng)驗值,Tw=10,由于本文主要考慮抑制區(qū)域間低頻振蕩,這個數(shù)值可以避免低頻時產(chǎn)生較大的相位超前,而破壞區(qū)域間暫態(tài)穩(wěn)定性;測量環(huán)節(jié)時間常數(shù)T7=0.01;滯后時間常數(shù)T2=T4=T6=0.2;超前時間常數(shù)T1、T3、T5為待優(yōu)化參數(shù),考慮到實際情況,每個超前或滯后環(huán)節(jié)補償量限于最大值60°,其取值范圍為[0.01,5.0];增益Ks為待優(yōu)化參數(shù),其取值范圍為[0.01,2.0]。

    參數(shù)優(yōu)化采用的預(yù)想故障為L7-8在母線7側(cè)發(fā)生三相短路,0.2 s故障,經(jīng)5個周波0.3 s故障消失。應(yīng)用本文方法優(yōu)化后的直流調(diào)制控制器參數(shù)列于表1,其對應(yīng)系統(tǒng)機電振蕩模式特征值列于表2。表2中還列出了無直流調(diào)制控制器系統(tǒng)的機電振蕩模式特征值。由表2可見,沒有安裝直流調(diào)制控制器,系統(tǒng)機電模式最小阻尼比為0.017,一般要求機電模式最小阻尼比不小于5%,故系統(tǒng)阻尼較弱。而安裝本文設(shè)計的直流調(diào)制控制器后,系統(tǒng)機電模式特征值明顯向復(fù)平面左平面大幅移動,提高了系統(tǒng)機電模式阻尼特性,可有效抑制系統(tǒng)的低頻振蕩。

    圖5 IEEE四機算例系統(tǒng)

    表1 直流調(diào)制參數(shù)

    表2 IEEE四機算例系統(tǒng)機電振蕩模式特征根比較

    在上述故障條件下,分別對無直流調(diào)制控制器、直流調(diào)制控制器采用預(yù)選參數(shù)和本文優(yōu)化后參數(shù)的三種情況進(jìn)行時域仿真,仿真結(jié)果如圖6~7。由圖6~7可見,本文優(yōu)化設(shè)計的直流調(diào)制控制器可以有效提高系統(tǒng)的阻尼特性,抑制大擾動后的系統(tǒng)振蕩,使得故障后系統(tǒng)發(fā)電機的功角快速衰減,系統(tǒng)功率振蕩得到快速抑制。

    圖6 1號發(fā)電機相對4號發(fā)電機的功角搖擺曲線

    圖7 交流聯(lián)絡(luò)線L7-8有功功率振蕩曲線

    直流調(diào)制控制器參數(shù)優(yōu)化設(shè)計方法的一個重要要求是魯棒性,優(yōu)化結(jié)果應(yīng)確保直流調(diào)制在系統(tǒng)運行方式或系統(tǒng)結(jié)構(gòu)發(fā)生變化后,仍然可以有效抑制系統(tǒng)振蕩。上面已經(jīng)說明了L7-8線路三相短路故障經(jīng)5個周波后消失,優(yōu)化參數(shù)對抑制系統(tǒng)振蕩的良好效果。直流調(diào)制控制器采用優(yōu)化后參數(shù)分別對以下兩種情況進(jìn)行時域仿真,驗證本文優(yōu)化方法的魯棒性。

    (1)圖5系統(tǒng)中,L8-9的一回線路在母線8側(cè)發(fā)生三相短路,0.2 s時故障,經(jīng)5個周波0.3 s切除一條線路。仿真結(jié)果如圖8,驗證了該優(yōu)化方法對系統(tǒng)結(jié)構(gòu)改變具有一定的魯棒性。

    (2)圖5系統(tǒng)中,改變系統(tǒng)運行方式,交流聯(lián)絡(luò)線L7-8雙回線的輸送功率由201.6 MW增加到298.2 MW,故障取L7-8線路0.2 s在母線7側(cè)三相短路經(jīng)5個周波后故障消失。仿真結(jié)果如圖9,驗證了本文優(yōu)化方法對系統(tǒng)運行方式改變具有一定的魯棒性。

    圖8 交流聯(lián)絡(luò)線L7-8有功功率振蕩曲線

    圖9 交流聯(lián)絡(luò)線L7-8有功功率振蕩曲線

    4 結(jié)語

    本文提出了一種基于時域仿真的直流調(diào)制參數(shù)優(yōu)化方法,優(yōu)化后的直流調(diào)制控制器可以有效提高系統(tǒng)整體的阻尼特性,抑制系統(tǒng)振蕩,使得系統(tǒng)快速恢復(fù)到穩(wěn)定運行狀態(tài)。通過在IEEE四機系統(tǒng)上的特征值分析和時域仿真結(jié)果,驗證了該方法的有效性和魯棒性。將本文提出方法應(yīng)用到實際復(fù)雜電力系統(tǒng)是進(jìn)一步要研究的重點。

    [1]Kundar P.Power System Stability and Control[M].New York:McGraw-Hill,1994.

    [2]Peterson H A,Krause P C.Damping of power swing in a parallel AC and DC system [J].IEEE Trans on Power Apparatus and Systems,1966,85(12):1231-1239.

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

    [4]Cresap R L,Scott D N,Mittelstadt W A.Small-signal modulation of the Pacific HVDC Inter-tie[J].IEEE Trans on Power Apparatus and Systems,1976,95(2):536-541.

    [5]沈梁,陳陳,史慧杰,等(Shen Liang,Chen Chen,Shi Huijie,etal).直流調(diào)制對電網(wǎng)區(qū)間低頻振蕩的抑制作用(Suppression of inter-area low frequency oscillation by HVDC modulation)[J].電力系統(tǒng)及其自動化學(xué)報(Proceedings of the CSU-EPSA),2008,20(4):82-86.

    [6]荊勇,洪潮,楊晉柏,等(Jing Yong,Hong Chao,Yang Jinbai,etal).直流調(diào)制抑制南方電網(wǎng)區(qū)域功率振蕩的研究(Suppression of inter-area power oscillation in southern China power grid by HVDC modulation)[J].電網(wǎng)技術(shù)(Power System Technology),2005,29(20):53-56.

    [7]荊勇,楊晉柏,李柏青,等(Jing Yong,Yang Jinbai,Li Baiqing,etal).直流調(diào)制改善交直流混聯(lián)系統(tǒng)暫態(tài)穩(wěn)定性的研究(Research on improving transient stability of AC/DC hybrid system by HVDC modulation)[J].電網(wǎng)技術(shù)(Power System Technology),2004,28(10):1-4.

    [8]劉紅超,李興源,王路,等(Liu Hongchao,Li Xingyuan,Wang Lu,etal).多饋入直流輸電系統(tǒng)中直流調(diào)制的協(xié)調(diào)優(yōu)化(Coordination and optimization of HVDC modulations in multi-infeed HVDC transmission system)[J].電網(wǎng)技術(shù)(Power System Technology),2004,28(1):5-9.

    [9]鄭希云,李興源,王渝紅(Zheng Xiyun,Li Xingyuan,Wang Yuhong).基于混沌優(yōu)化算法的PSS和直流調(diào)制的協(xié)調(diào)優(yōu)化(Coordination and optimization of PSS and HVDC modulations using chaotic optimization algorithm)[J].電工技術(shù)學(xué)報(Transactions of China Electrotechnical Society),2010,25(5):170-175.

    [10]毛曉明,張堯,管霖,等(Mao Xiaoming,Zhang Yao,Guan Lin,etal).南方交直流混合電網(wǎng)區(qū)域振蕩的協(xié)調(diào)控制策略(Research on coordinated control strategy for inter-area oscillations in China southern AC/DC power grid)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems),2005,29(20):55-59.

    [11]周孝法,陳陳,楊帆,等(Zhou Xiaofa,Chen Chen,Yang Fan,etal).基于自適應(yīng)混沌粒子群優(yōu)化算法的多饋入直流輸電系統(tǒng)優(yōu)化協(xié)調(diào)直流調(diào)制(Optima coordinated HVDC modulation based on adaptive chaos particle swarm optimization algorithm in multi-infeed HVDC transmission system)[J].電工技術(shù)學(xué)報(Transactions of China Electrotechnical Society),2009,24(4):193-201.

    [12]金麗成,劉海峰,徐政(Jin Licheng,Liu Haifeng,Xu Zheng).多饋入直流輸電系統(tǒng)小信號調(diào)制器的協(xié)調(diào)優(yōu)化整定(Coordination and optimization of small signal modulation in multi-infeed HVDC systems)[J].電力系統(tǒng)自動化(Automation of Electric Power Systems),2003,27(16):10-15.

    [13]周保榮,房大中,孫景強(Zhou Baorong,Fang Dazhong,Sun Jingqiang).基于軌跡靈敏度分析的電力系統(tǒng)穩(wěn)定器參數(shù)優(yōu)化設(shè)計(Tuning of PSS parameters using optimization approach based on trajectory sensitivity analysis)[J].電網(wǎng)技術(shù)(Power System Technology),2004,28(19):20-23,84.

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

    [15]Hiskens Ian A,Pai M A.Trajectory sensitivity analysis of hybrid systems [J].IEEE Trans on Circuits and Systems,2000,47(2):204-220.

    [16]解可新,韓健,林友聯(lián).最優(yōu)化方法[M].天津:天津大學(xué)出版社,1997.

    王長香(1986-),女,碩士研究生,研究方向為電力系統(tǒng)穩(wěn)定性。Email:daisywang_415@yahoo.cn

    房大中(1946-),男,教授,博士生導(dǎo)師,主要從事電力系統(tǒng)穩(wěn)定性分析與控制研究。Email:dz_fang@yahoo.com.cn

    薛振宇(1983-),男,博士研究生,研究方向為電力系統(tǒng)安全分析和數(shù)據(jù)庫等。Email:banditxue@163.com

    ParameterOptimizationDesignofDCModulationBasedonTrajectorySensitivityAnalysis

    WANG Chang-xiang1,F(xiàn)ANG Da-zhong1,XUE Zhen-yu1,WEI Wen-hui2

    (1.School of Electrical Engineering and Automation,Tianjin University,Tianjin 300072,China;2.China Electric Power Research Institute,Beijing 100192,China)

    A method for parameter-optimization of DC-modulation controller on the basis of time domain simulation has been proposed in this paper.In this method,the parameter optimization problem of DC-modulation controller is converted to the non-linear optimization problem with parameter constraint,and then the constrained optimization problem has been converted to the unconstrained optimization problem.The gradient information of objective function is obtained by means of the simulation of trajectory sensitivity,and then the parameters of DC-modulation controller are optimized by Quasi-Newton Algorithm(DFP).Because the objective function reflects the power oscillation behavior of AC/DC interconnected power systems,the proposed method can be used to effectively damp the power oscillations of AC/DC system incurred by severe disturbances.The eigenvalue analysis and time domain simulations on the IEEE 4-generator AC/DC system verify the effectiveness and robustness of the proposed method.

    DC modulation;trajectory sensitivity;Quasi-Newton algorithm;parameter optimization

    TM712;TM743

    A

    1003-8930(2012)04-0030-06

    2011-08-22;

    2011-09-16

    國家自然科學(xué)基金資助項目(50777046);國家重點基礎(chǔ)研究發(fā)展計劃(973)資助項目(2009CB219700)

    猜你喜歡
    阻尼靈敏度直流
    基于直流載波通信的LAMOST控制系統(tǒng)設(shè)計
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項的擴散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    導(dǎo)磁環(huán)對LVDT線性度和靈敏度的影響
    一款高效的30V直流開關(guān)電源設(shè)計
    地下水非穩(wěn)定流的靈敏度分析
    具阻尼項的Boussinesq型方程的長時間行為
    穿甲爆破彈引信對薄弱目標(biāo)的靈敏度分析
    非隔離型光伏并網(wǎng)逆變器直流注入抑制方法
    x7x7x7水蜜桃| 亚洲中文日韩欧美视频| 亚洲一区二区三区不卡视频| 欧美亚洲日本最大视频资源| 精品乱码久久久久久99久播| 日本免费a在线| 最近最新免费中文字幕在线| 在线观看日韩欧美| 国产片内射在线| av电影中文网址| 成人av一区二区三区在线看| av视频免费观看在线观看| 日韩欧美免费精品| 日韩一卡2卡3卡4卡2021年| 黑人欧美特级aaaaaa片| 国产精品亚洲一级av第二区| 亚洲精品久久成人aⅴ小说| 国产精品九九99| 国产黄色免费在线视频| √禁漫天堂资源中文www| 亚洲午夜理论影院| 亚洲av成人不卡在线观看播放网| 国产精品九九99| 日韩高清综合在线| 国产人伦9x9x在线观看| 手机成人av网站| 日本五十路高清| 亚洲av成人av| 超色免费av| 一级片'在线观看视频| 男男h啪啪无遮挡| 国产99白浆流出| 91麻豆精品激情在线观看国产 | 人人澡人人妻人| 欧美精品亚洲一区二区| 80岁老熟妇乱子伦牲交| 国产精品成人在线| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| 精品人妻在线不人妻| 人成视频在线观看免费观看| 日韩欧美在线二视频| 欧美激情久久久久久爽电影 | 亚洲av第一区精品v没综合| 两人在一起打扑克的视频| 亚洲国产欧美网| 亚洲伊人色综图| 91成年电影在线观看| 亚洲人成电影免费在线| 久久精品国产综合久久久| 日韩精品青青久久久久久| 国产午夜精品久久久久久| 一级作爱视频免费观看| 国产亚洲精品综合一区在线观看 | 国产亚洲欧美精品永久| 变态另类成人亚洲欧美熟女 | 99riav亚洲国产免费| 久久久国产一区二区| 欧美亚洲日本最大视频资源| 国产黄a三级三级三级人| 一区在线观看完整版| 看免费av毛片| 国产深夜福利视频在线观看| 国产伦人伦偷精品视频| 国产成人啪精品午夜网站| 9热在线视频观看99| 亚洲国产看品久久| 一边摸一边抽搐一进一小说| 亚洲精品久久午夜乱码| 日本五十路高清| 久久中文字幕一级| www.熟女人妻精品国产| 悠悠久久av| 久久久水蜜桃国产精品网| 一级作爱视频免费观看| 又黄又粗又硬又大视频| 午夜福利一区二区在线看| 久久天躁狠狠躁夜夜2o2o| 在线播放国产精品三级| 12—13女人毛片做爰片一| 国产黄a三级三级三级人| 欧美日韩福利视频一区二区| 亚洲专区字幕在线| 丰满的人妻完整版| 青草久久国产| 黑人操中国人逼视频| 最好的美女福利视频网| 一夜夜www| 老司机福利观看| 国产精品98久久久久久宅男小说| 免费高清视频大片| 桃色一区二区三区在线观看| 777久久人妻少妇嫩草av网站| 国产精品久久久人人做人人爽| 在线观看午夜福利视频| 中文字幕av电影在线播放| 伊人久久大香线蕉亚洲五| 90打野战视频偷拍视频| 熟女少妇亚洲综合色aaa.| 黄色成人免费大全| 久久香蕉国产精品| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| 午夜91福利影院| 色在线成人网| 久久精品91无色码中文字幕| 黄色毛片三级朝国网站| 亚洲第一av免费看| 中文欧美无线码| 国产精品电影一区二区三区| 午夜免费成人在线视频| 电影成人av| 成人18禁高潮啪啪吃奶动态图| 亚洲精品国产区一区二| 亚洲av日韩精品久久久久久密| 免费不卡黄色视频| 侵犯人妻中文字幕一二三四区| 88av欧美| 国产精品亚洲av一区麻豆| 美女国产高潮福利片在线看| 天堂动漫精品| 性欧美人与动物交配| 一边摸一边抽搐一进一出视频| 欧美黑人精品巨大| 国产极品粉嫩免费观看在线| 岛国在线观看网站| 99国产精品免费福利视频| 18禁黄网站禁片午夜丰满| 香蕉久久夜色| 亚洲成人国产一区在线观看| 操美女的视频在线观看| 在线观看66精品国产| 午夜福利,免费看| 国产欧美日韩综合在线一区二区| 欧美激情久久久久久爽电影 | 91国产中文字幕| 亚洲 欧美一区二区三区| 日韩欧美国产一区二区入口| 久久精品91无色码中文字幕| 一级毛片高清免费大全| 在线观看免费视频网站a站| 桃红色精品国产亚洲av| 啦啦啦在线免费观看视频4| 亚洲成人免费av在线播放| 母亲3免费完整高清在线观看| 最近最新中文字幕大全电影3 | 一级毛片高清免费大全| 久久久国产欧美日韩av| 桃红色精品国产亚洲av| 久久久久久久久中文| 在线视频色国产色| 男人操女人黄网站| 亚洲av成人一区二区三| 一进一出好大好爽视频| 桃色一区二区三区在线观看| 久久香蕉精品热| 五月开心婷婷网| 成人国语在线视频| 日韩中文字幕欧美一区二区| 美女高潮到喷水免费观看| 国产97色在线日韩免费| 成人特级黄色片久久久久久久| 美女高潮到喷水免费观看| 成年女人毛片免费观看观看9| 亚洲精品一卡2卡三卡4卡5卡| 老汉色av国产亚洲站长工具| 久久久久久亚洲精品国产蜜桃av| 亚洲精品成人av观看孕妇| 国产一区二区三区综合在线观看| 欧美另类亚洲清纯唯美| 欧美+亚洲+日韩+国产| 国产成人精品无人区| 久久人人爽av亚洲精品天堂| 国产深夜福利视频在线观看| 久久久久九九精品影院| 国产成人精品久久二区二区免费| 少妇被粗大的猛进出69影院| 免费不卡黄色视频| 欧美日韩黄片免| 天天添夜夜摸| 91老司机精品| 亚洲第一av免费看| 国产国语露脸激情在线看| 免费日韩欧美在线观看| 亚洲男人天堂网一区| 99久久精品国产亚洲精品| 9热在线视频观看99| 国产一区二区三区综合在线观看| 国产成人啪精品午夜网站| 少妇粗大呻吟视频| 国产精品久久视频播放| 国产成+人综合+亚洲专区| 男女下面进入的视频免费午夜 | 久久香蕉国产精品| 国产人伦9x9x在线观看| 国产亚洲精品第一综合不卡| 欧美色视频一区免费| 国产免费男女视频| 欧美激情极品国产一区二区三区| 久久久久精品国产欧美久久久| 91字幕亚洲| 电影成人av| 91精品国产国语对白视频| 久久午夜亚洲精品久久| 午夜日韩欧美国产| 美女高潮喷水抽搐中文字幕| 成人亚洲精品av一区二区 | 十分钟在线观看高清视频www| 欧美精品亚洲一区二区| 国产亚洲精品久久久久5区| 人人澡人人妻人| 亚洲国产精品sss在线观看 | 久久久久久免费高清国产稀缺| 叶爱在线成人免费视频播放| 亚洲色图综合在线观看| 久久天堂一区二区三区四区| 精品久久久久久电影网| 一级a爱视频在线免费观看| 视频区欧美日本亚洲| √禁漫天堂资源中文www| 自线自在国产av| 亚洲精华国产精华精| 中文字幕高清在线视频| 男女床上黄色一级片免费看| 国产精品永久免费网站| 日本欧美视频一区| 免费在线观看亚洲国产| 不卡一级毛片| 18美女黄网站色大片免费观看| 在线观看免费高清a一片| 亚洲片人在线观看| 性欧美人与动物交配| 麻豆久久精品国产亚洲av | 身体一侧抽搐| 99国产精品一区二区蜜桃av| 精品高清国产在线一区| 一二三四社区在线视频社区8| 亚洲欧美精品综合久久99| 12—13女人毛片做爰片一| 免费搜索国产男女视频| 不卡av一区二区三区| www.999成人在线观看| 高清av免费在线| 久久久久久人人人人人| 国产精品 国内视频| 91在线观看av| 亚洲熟妇熟女久久| 精品第一国产精品| 欧美黑人精品巨大| 色综合欧美亚洲国产小说| 老熟妇仑乱视频hdxx| 交换朋友夫妻互换小说| 老司机深夜福利视频在线观看| 亚洲熟妇熟女久久| 精品久久久久久久久久免费视频 | 夜夜夜夜夜久久久久| 亚洲 欧美 日韩 在线 免费| 中亚洲国语对白在线视频| 午夜91福利影院| 国产精品一区二区三区四区久久 | 日韩高清综合在线| 97人妻天天添夜夜摸| 黑丝袜美女国产一区| 怎么达到女性高潮| 1024香蕉在线观看| 黄色丝袜av网址大全| 欧美av亚洲av综合av国产av| 亚洲成a人片在线一区二区| 国产一区二区三区视频了| 交换朋友夫妻互换小说| 免费日韩欧美在线观看| 国产亚洲精品久久久久5区| 中文字幕av电影在线播放| 亚洲视频免费观看视频| svipshipincom国产片| av中文乱码字幕在线| 波多野结衣一区麻豆| 手机成人av网站| 国产av精品麻豆| 制服人妻中文乱码| 免费av毛片视频| 女人被狂操c到高潮| 日本免费一区二区三区高清不卡 | 亚洲情色 制服丝袜| 国产又爽黄色视频| 首页视频小说图片口味搜索| 国产精品影院久久| 国产麻豆69| 丝袜在线中文字幕| 免费看十八禁软件| 久久国产精品影院| 亚洲精华国产精华精| 欧美性长视频在线观看| 香蕉丝袜av| 黄色毛片三级朝国网站| 日韩欧美三级三区| 两个人看的免费小视频| 18禁国产床啪视频网站| 欧美大码av| 国产精品香港三级国产av潘金莲| 婷婷丁香在线五月| 日本欧美视频一区| 精品日产1卡2卡| 女人高潮潮喷娇喘18禁视频| x7x7x7水蜜桃| 人人妻,人人澡人人爽秒播| 成人黄色视频免费在线看| 久热爱精品视频在线9| 久久精品国产亚洲av高清一级| 亚洲午夜理论影院| 嫩草影院精品99| 亚洲精品av麻豆狂野| 亚洲,欧美精品.| 国产激情久久老熟女| 天天影视国产精品| 91麻豆av在线| 亚洲成人精品中文字幕电影 | 久久久久久免费高清国产稀缺| 亚洲五月天丁香| 日韩国内少妇激情av| 久久中文字幕人妻熟女| 黑人操中国人逼视频| 亚洲中文av在线| 9色porny在线观看| 亚洲久久久国产精品| 国产人伦9x9x在线观看| 国产一卡二卡三卡精品| 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产 | 女生性感内裤真人,穿戴方法视频| 国产日韩一区二区三区精品不卡| 1024香蕉在线观看| 欧美久久黑人一区二区| 他把我摸到了高潮在线观看| 国产aⅴ精品一区二区三区波| 婷婷六月久久综合丁香| 成人18禁在线播放| 久久天堂一区二区三区四区| 久久午夜综合久久蜜桃| 黄片小视频在线播放| svipshipincom国产片| 日韩精品免费视频一区二区三区| 法律面前人人平等表现在哪些方面| 久久国产精品男人的天堂亚洲| 成人国语在线视频| 啦啦啦免费观看视频1| 性少妇av在线| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人免费av一区二区三区| 国产av在哪里看| 免费在线观看亚洲国产| 日韩精品青青久久久久久| 国产精品美女特级片免费视频播放器 | netflix在线观看网站| 女生性感内裤真人,穿戴方法视频| 亚洲精品久久午夜乱码| 女性生殖器流出的白浆| 看黄色毛片网站| 国产一区二区激情短视频| 窝窝影院91人妻| 人人妻人人爽人人添夜夜欢视频| 国产成人av激情在线播放| 欧美黄色片欧美黄色片| 久久青草综合色| 亚洲成av片中文字幕在线观看| 精品第一国产精品| 欧美日本中文国产一区发布| 男男h啪啪无遮挡| videosex国产| 欧美日韩国产mv在线观看视频| 久久久久久久午夜电影 | 久久精品国产亚洲av高清一级| 免费在线观看视频国产中文字幕亚洲| 亚洲av熟女| 岛国在线观看网站| 韩国av一区二区三区四区| 国产精品自产拍在线观看55亚洲| 亚洲一区二区三区不卡视频| 久久99一区二区三区| 久久久久国产精品人妻aⅴ院| 性色av乱码一区二区三区2| 超碰成人久久| av超薄肉色丝袜交足视频| 亚洲欧美精品综合一区二区三区| 18禁裸乳无遮挡免费网站照片 | 91老司机精品| 美女高潮到喷水免费观看| 国产1区2区3区精品| 视频区欧美日本亚洲| 亚洲午夜理论影院| 18禁美女被吸乳视频| 欧美人与性动交α欧美软件| 啪啪无遮挡十八禁网站| 亚洲三区欧美一区| 一级作爱视频免费观看| 午夜免费成人在线视频| 精品午夜福利视频在线观看一区| av在线播放免费不卡| 国产成人啪精品午夜网站| 国产av一区二区精品久久| www.熟女人妻精品国产| 亚洲久久久国产精品| 欧美最黄视频在线播放免费 | 久久热在线av| 天堂√8在线中文| 国产精品电影一区二区三区| 成在线人永久免费视频| 老司机亚洲免费影院| 夜夜爽天天搞| 亚洲熟女毛片儿| 国产精品成人在线| 亚洲全国av大片| 热99re8久久精品国产| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色 | 黄色 视频免费看| 中文字幕人妻熟女乱码| 婷婷精品国产亚洲av在线| 性欧美人与动物交配| 日韩 欧美 亚洲 中文字幕| 色尼玛亚洲综合影院| 一进一出好大好爽视频| 欧美最黄视频在线播放免费 | 高潮久久久久久久久久久不卡| 亚洲熟妇中文字幕五十中出 | 亚洲七黄色美女视频| 一边摸一边抽搐一进一出视频| 中国美女看黄片| 欧美日韩黄片免| 国产国语露脸激情在线看| 欧美日本亚洲视频在线播放| 熟女少妇亚洲综合色aaa.| 久久精品国产清高在天天线| 一级毛片精品| 午夜亚洲福利在线播放| 水蜜桃什么品种好| 中文欧美无线码| 91字幕亚洲| 久久久久久亚洲精品国产蜜桃av| 国产高清激情床上av| 精品第一国产精品| 90打野战视频偷拍视频| 欧美激情久久久久久爽电影 | 两个人看的免费小视频| 亚洲熟妇中文字幕五十中出 | 性色av乱码一区二区三区2| 亚洲av五月六月丁香网| 香蕉国产在线看| 黑人巨大精品欧美一区二区蜜桃| 精品欧美一区二区三区在线| 18禁国产床啪视频网站| 亚洲中文日韩欧美视频| 久久久久久久精品吃奶| 亚洲人成网站在线播放欧美日韩| 99香蕉大伊视频| 久久人人精品亚洲av| 亚洲国产中文字幕在线视频| 亚洲熟妇中文字幕五十中出 | 久久精品91无色码中文字幕| 12—13女人毛片做爰片一| 性少妇av在线| 久久精品国产亚洲av香蕉五月| 国产成人av激情在线播放| 一个人免费在线观看的高清视频| 国产精品国产高清国产av| 嫁个100分男人电影在线观看| 成人av一区二区三区在线看| 亚洲成国产人片在线观看| 黄色女人牲交| 国产成人免费无遮挡视频| 欧美激情高清一区二区三区| 麻豆成人av在线观看| 精品熟女少妇八av免费久了| 免费少妇av软件| 最新美女视频免费是黄的| 琪琪午夜伦伦电影理论片6080| 国产免费男女视频| 亚洲五月婷婷丁香| 国产av在哪里看| 法律面前人人平等表现在哪些方面| 亚洲成av片中文字幕在线观看| 国产亚洲精品综合一区在线观看 | 成人av一区二区三区在线看| 国产成人欧美在线观看| 久久中文看片网| 亚洲中文日韩欧美视频| 啪啪无遮挡十八禁网站| 国产精品综合久久久久久久免费 | 91大片在线观看| 国产免费现黄频在线看| 午夜激情av网站| 国产熟女午夜一区二区三区| 亚洲自拍偷在线| 搡老乐熟女国产| 视频区欧美日本亚洲| 亚洲一区中文字幕在线| 久久香蕉激情| 国产一区在线观看成人免费| 久久人人爽av亚洲精品天堂| 国产有黄有色有爽视频| 夜夜看夜夜爽夜夜摸 | 99久久久亚洲精品蜜臀av| 欧美黄色淫秽网站| 亚洲视频免费观看视频| 亚洲午夜理论影院| 国产黄色免费在线视频| 99国产极品粉嫩在线观看| 午夜成年电影在线免费观看| 国产亚洲欧美在线一区二区| 9热在线视频观看99| 久久香蕉激情| 午夜福利欧美成人| 精品国内亚洲2022精品成人| 满18在线观看网站| 国产成人欧美| 99久久综合精品五月天人人| 天堂动漫精品| www日本在线高清视频| 欧美乱妇无乱码| 精品久久久久久,| 9色porny在线观看| 日本免费a在线| 久久精品91蜜桃| 女人精品久久久久毛片| 精品乱码久久久久久99久播| av电影中文网址| 一区二区三区精品91| 亚洲少妇的诱惑av| av天堂久久9| 亚洲狠狠婷婷综合久久图片| 国产视频一区二区在线看| 性少妇av在线| 欧美激情 高清一区二区三区| 国产熟女xx| 嫁个100分男人电影在线观看| 女人高潮潮喷娇喘18禁视频| 最新美女视频免费是黄的| 叶爱在线成人免费视频播放| 精品国产亚洲在线| 可以免费在线观看a视频的电影网站| 亚洲av美国av| 不卡av一区二区三区| 麻豆国产av国片精品| 精品国产一区二区三区四区第35| 精品国产乱子伦一区二区三区| 亚洲欧美一区二区三区黑人| 麻豆久久精品国产亚洲av | 十八禁网站免费在线| 欧美精品啪啪一区二区三区| 在线观看日韩欧美| 久热这里只有精品99| 丝袜美足系列| 一级,二级,三级黄色视频| 亚洲av成人不卡在线观看播放网| 国产片内射在线| 大陆偷拍与自拍| 精品国产亚洲在线| 精品午夜福利视频在线观看一区| 欧美激情极品国产一区二区三区| 国产精品一区二区精品视频观看| 免费在线观看亚洲国产| 99热只有精品国产| 一本大道久久a久久精品| 后天国语完整版免费观看| 国产不卡一卡二| 亚洲精品av麻豆狂野| 最近最新免费中文字幕在线| 免费看a级黄色片| 久久国产精品人妻蜜桃| 国产精品国产av在线观看| 香蕉久久夜色| 99在线视频只有这里精品首页| 国产成人系列免费观看| 窝窝影院91人妻| 色精品久久人妻99蜜桃| 欧美日本亚洲视频在线播放| 精品国产乱子伦一区二区三区| cao死你这个sao货| 乱人伦中国视频| 亚洲情色 制服丝袜| 亚洲一码二码三码区别大吗| 自线自在国产av| 成人18禁高潮啪啪吃奶动态图| 亚洲成a人片在线一区二区| 国产成年人精品一区二区 | 久久午夜亚洲精品久久| 在线观看一区二区三区| 国产成人精品久久二区二区免费| 桃红色精品国产亚洲av| 国产真人三级小视频在线观看| 91精品国产国语对白视频| 久久久久国内视频| 色在线成人网| 丁香六月欧美| www.www免费av| 国产成人系列免费观看| 久久国产亚洲av麻豆专区| 久久久久九九精品影院| 国产一区在线观看成人免费| 丝袜美足系列| 欧美成狂野欧美在线观看| 日韩 欧美 亚洲 中文字幕| 午夜日韩欧美国产| 丝袜人妻中文字幕| 男女午夜视频在线观看| 亚洲 欧美一区二区三区| 免费看十八禁软件| 亚洲第一av免费看| 色综合婷婷激情| 国产av又大| 国产蜜桃级精品一区二区三区| 亚洲专区字幕在线| 99热国产这里只有精品6| 91麻豆精品激情在线观看国产 | 真人一进一出gif抽搐免费|