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

    光伏并網(wǎng)發(fā)電系統(tǒng)參數(shù)協(xié)調(diào)優(yōu)化方法研究

    2016-10-21 01:13:38西南交通大學(xué)電氣工程學(xué)院四川成都610031
    西南交通大學(xué)學(xué)報 2016年5期
    關(guān)鍵詞:特征值種群發(fā)電

    (西南交通大學(xué)電氣工程學(xué)院,四川成都610031)

    (西南交通大學(xué)電氣工程學(xué)院,四川成都610031)

    光伏并網(wǎng)發(fā)電系統(tǒng)的參數(shù)設(shè)計對其穩(wěn)定運(yùn)行具有重要意義,為此,提出了一種光伏并網(wǎng)發(fā)電系統(tǒng)參數(shù)的協(xié)調(diào)優(yōu)化方法.該方法通過建立包含光伏電池、并網(wǎng)逆變器、變壓器和電網(wǎng)的光伏發(fā)電系統(tǒng)小信號數(shù)學(xué)模型,利用特征值分析(根軌跡圖)得到各個特征值穩(wěn)定情況下的系統(tǒng)參數(shù)邊界值;在此基礎(chǔ)上,建立了計及小干擾穩(wěn)定性、阻尼比和穩(wěn)定裕度的協(xié)調(diào)優(yōu)化目標(biāo)函數(shù),并采用回溯搜索算法(backtracking search algorithm,BSA)對系統(tǒng)參數(shù)進(jìn)行優(yōu)化;通過仿真測試,對比分析了參數(shù)優(yōu)化前后系統(tǒng)的動態(tài)響應(yīng).研究結(jié)果表明,在系統(tǒng)遭受到小擾動后,經(jīng)過參數(shù)優(yōu)化的光伏并網(wǎng)系統(tǒng)能夠在0.1 s內(nèi)達(dá)到新的平衡點(diǎn)穩(wěn)定運(yùn)行,比未優(yōu)化參數(shù)的系統(tǒng)動態(tài)響應(yīng)提高了50%.

    光伏發(fā)電系統(tǒng);小信號模型;回溯搜索算法;協(xié)調(diào)優(yōu)化目標(biāo)函數(shù);特征值分析法

    分布式發(fā)電系統(tǒng)具有環(huán)境污染小、投資少、效率高和靈活性強(qiáng)等優(yōu)點(diǎn),大量光伏、風(fēng)力和燃料電池等分布式電源并入電網(wǎng)將成為智能電網(wǎng)發(fā)展的必然趨勢[1-5],對其進(jìn)行合理的并網(wǎng)控制和穩(wěn)定性分析至關(guān)重要[6-8].在光伏并網(wǎng)發(fā)電系統(tǒng)運(yùn)行的過程中,太陽光強(qiáng)度和環(huán)境溫度改變、端電壓下降、負(fù)荷參數(shù)變化等因素都會影響系統(tǒng)運(yùn)行的穩(wěn)定性.隨著光伏發(fā)電行業(yè)的迅速發(fā)展,大規(guī)模光伏發(fā)電系統(tǒng)接入電網(wǎng)使得對其穩(wěn)定性的要求越來越高.因此,研究影響光伏發(fā)電系統(tǒng)運(yùn)行穩(wěn)定性的因素和提高其運(yùn)行穩(wěn)定性的改善措施具有重要意義.

    近年來,國內(nèi)外大部分研究工作主要是針對光伏發(fā)電并網(wǎng)系統(tǒng)建模和對接口逆變器控制方法的研究,對從控制角度進(jìn)行穩(wěn)定性分析的研究工作相對較少.文獻(xiàn)[9]研究了光伏發(fā)電系統(tǒng)的等值模型,但該模型沒有全面考慮到系統(tǒng)逆變器模型和并網(wǎng)模型,不利于研究系統(tǒng)運(yùn)行穩(wěn)定性.文獻(xiàn)[10]通過實(shí)驗得到了光伏電池輸出電壓和輸出功率對系統(tǒng)穩(wěn)定性的影響.文獻(xiàn)[11]通過建立光伏發(fā)電系統(tǒng)的小信號模型,利用特征值分析了系統(tǒng)遭受小干擾后的穩(wěn)定性.文獻(xiàn)[12]建立了系統(tǒng)孤島運(yùn)行時的全階小信號模型,得到系統(tǒng)運(yùn)行中存在的3種主導(dǎo)振蕩模式并分析了控制器和網(wǎng)絡(luò)參數(shù)變化對系統(tǒng)小干擾穩(wěn)定性的影響.但文獻(xiàn)[10-12]沒有給出提高系統(tǒng)穩(wěn)定性和控制器參數(shù)的優(yōu)化方法.文獻(xiàn)[13]在小信號模型的基礎(chǔ)上,提出了根據(jù)系統(tǒng)特征值軌跡來定性設(shè)計光伏發(fā)電系統(tǒng)的控制器參數(shù)的方法,但該方法未詳細(xì)考慮各個控制器參數(shù)之間的耦合關(guān)系,僅分析了參數(shù)變化對系統(tǒng)穩(wěn)定性的影響,未能使系統(tǒng)參數(shù)達(dá)到相對全局最優(yōu).

    回溯搜索算法(BSA)是一種基于種群循環(huán)的進(jìn)化算法,受啟發(fā)于自然界生物種群在尋找食物過程中會隨機(jī)前往曾經(jīng)發(fā)現(xiàn)豐富食物的地點(diǎn)搜尋的習(xí)性,在迭代過程中隨機(jī)選擇一個歷史種群來指導(dǎo)當(dāng)前種群的進(jìn)化,具有良好的全局收斂性和收斂速度快的優(yōu)點(diǎn).

    本文在所建立的光伏發(fā)電系統(tǒng)小信號模型基礎(chǔ)上,建立能夠全面衡量系統(tǒng)穩(wěn)定性的協(xié)調(diào)優(yōu)化目標(biāo)函數(shù),通過特征值根軌跡圖確定系統(tǒng)參數(shù)的上下界,采用BSA對系統(tǒng)參數(shù)進(jìn)行尋優(yōu),提高系統(tǒng)穩(wěn)定性.

    1 小信號模型

    1.1 光伏發(fā)電并網(wǎng)系統(tǒng)結(jié)構(gòu)

    典型的光伏發(fā)電并網(wǎng)系統(tǒng)結(jié)構(gòu)如圖1所示,包括光伏陣列、儲能電容器CDC、電壓源型并網(wǎng)逆變器、濾波電感Lf、變比為380 V/10 kV的升壓變壓器、等效線路阻抗Zline和單機(jī)無窮大電力系統(tǒng).

    光伏電池經(jīng)過串并聯(lián)形成光伏陣列,光伏陣列輸出的直流電壓通過儲能電容器接入并網(wǎng)逆變器,逆變器輸出電壓經(jīng)過濾波電感,經(jīng)過升壓變壓器和輸電線路接入電網(wǎng).

    圖1 光伏發(fā)電系統(tǒng)的原理示意圖Fig.1 Schematic diagram of grid-connected photovoltaic system

    1.2 光伏發(fā)電系統(tǒng)模型

    光伏發(fā)電系統(tǒng)模型包括光伏電池模型、逆變器模型和并網(wǎng)模型3部分.

    光伏電池是在工程實(shí)用模型的基礎(chǔ)上,采用多項式擬合的方法得到的光伏電池線性擬合模型[14]:

    式中:I表示光伏電池的輸出電流;t表示在非標(biāo)準(zhǔn)條件下光伏電池的溫度;Isc和Uoc分別為光伏電池的短路電流和開路電壓;Im和Um分別為光伏電池工作在最大功率點(diǎn)時的電流和電壓;S為光照強(qiáng)度;k為溫度系數(shù);tair為環(huán)境溫度;l、m、n為補(bǔ)償系數(shù);ρ、σ為線性擬合參數(shù).

    光伏并網(wǎng)逆變器采用電壓源性逆變器,其控制策略采用直流母線電壓外環(huán)、電流內(nèi)環(huán)[17],如圖2所示.

    根據(jù)圖2得到并網(wǎng)逆變器控制方程如式(2)所示,其在dq軸坐標(biāo)下的數(shù)學(xué)模型如式(3)所示.

    圖2 光伏逆變器控制系統(tǒng)Fig.2 Diagram of PV inverter control system

    式中:ωB為電力系統(tǒng)角頻率基準(zhǔn)值;x1、x2、x3是狀態(tài)變量,單位為s;KP1、KI1為直流母線電壓外環(huán)控制器的比例參數(shù)和積分參數(shù);KP2、KI2、KP3和KI3分別為有功電流內(nèi)環(huán)和無功電流內(nèi)環(huán)控制器的比例參數(shù)和積分參數(shù);上標(biāo)*表示變量的參考值;ud表示dq坐標(biāo)逆變器輸出端的d軸電壓.

    直流母線儲能電容CDC的數(shù)學(xué)模型可根據(jù)功率平衡得到,

    光伏發(fā)電系統(tǒng)與電網(wǎng)之間的連接模型為

    式中:R和X分別為變壓器以及傳輸線路的等效電阻和等效電抗,R=Rline,X=XT+Xline;K為變壓器變比.

    式(1)~(5)完整描述了光伏并網(wǎng)發(fā)電系統(tǒng),將其在穩(wěn)態(tài)值附近進(jìn)行線性化,得到光伏電站發(fā)電系統(tǒng)的狀態(tài)矩陣如式(6)所示,其中,變量括號中的0表示取該標(biāo)量的初值.

    式中:x=[x1x2x3idgiqguDC]T;

    其中:

    2 特征值分析

    根據(jù)式(6)可以計算得到系統(tǒng)的所有特征值,給定系統(tǒng)各參數(shù)值如表1所示.

    表1 系統(tǒng)參數(shù)Tab.1 Parameters of system

    為了得到系統(tǒng)在穩(wěn)定狀態(tài)下,各控制器參數(shù)的取值范圍,畫出系統(tǒng)不同參數(shù)對應(yīng)的特征值根軌跡圖進(jìn)行分析,如圖3所示.

    由圖3可知,系統(tǒng)有6個特征根,其中λ5和λ6是否穩(wěn)定只與參數(shù)KP3和KI3的取值有關(guān),且隨著KP3和KI3的增大,λ5朝穩(wěn)定區(qū)域運(yùn)動,而λ6朝著不穩(wěn)定區(qū)域運(yùn)動.

    圖3 不同參數(shù)對應(yīng)的特征值軌跡Fig.3 Eigenvalue locus under different combinations of controller parameters

    參數(shù)KP1、KI1、KP2和KI2決定了特征根λ1、λ2、λ3和λ4是否穩(wěn)定.KP1的增大會使得λ3不穩(wěn)定,而增大KI1能夠使λ3朝穩(wěn)定方向移動.增大KP1、KI1和KI2均能使λ1,λ2變得穩(wěn)定,而增大KP2會導(dǎo)致λ2朝不穩(wěn)定區(qū)域移動.綜合圖3得到在系統(tǒng)所有特征值穩(wěn)定條件下的KP1、KI1、KP2、KI2、KP3和KI3的取值范圍如表2所示.

    表2 控制器參數(shù)取值范圍Tab.2 Range of control system parameters

    3 控制器參數(shù)優(yōu)化

    3.1 協(xié)調(diào)優(yōu)化目標(biāo)函數(shù)

    傳統(tǒng)的目標(biāo)函數(shù)將特征值實(shí)部小于0的部分相加,尋其最大值.這種方法僅考慮了特征值穩(wěn)定性,因此得到的優(yōu)化參數(shù)不夠全面.

    針對光伏并網(wǎng)發(fā)電系統(tǒng)運(yùn)行特點(diǎn),考慮到在系統(tǒng)穩(wěn)定運(yùn)行過程中的多方面要求,包括系統(tǒng)小干擾穩(wěn)定性、阻尼比和穩(wěn)定裕度,本文對文獻(xiàn)[15]中包含微電網(wǎng)運(yùn)行場景波動性、多變性和多運(yùn)行場景的綜合優(yōu)化目標(biāo)函數(shù)進(jìn)行簡化處理,得到協(xié)調(diào)優(yōu)化目標(biāo)函數(shù)如式(7)所示,

    式中:J1表示系統(tǒng)中實(shí)部為正的特征值所構(gòu)成的目標(biāo)函數(shù),其中:αi為系統(tǒng)中第i個實(shí)部大于0的特征值的實(shí)部,μi為相應(yīng)的權(quán)重系數(shù),可由μi=計算得到,L為大于0的經(jīng)驗參數(shù);

    J2表示系統(tǒng)特征值阻尼比小于給定阻尼比ξ0的特征值所構(gòu)成的目標(biāo)函數(shù),其中:ξ1為第j個小于ξ0的特征值阻尼比,為相應(yīng)的權(quán)重系數(shù),可由式(8)計算得到,其中ξ0=0.75,L′同樣為大于0的經(jīng)驗參數(shù),

    J3表示系統(tǒng)特征值中實(shí)部小于0,大于給定實(shí)部α0的特征值所構(gòu)成的目標(biāo)函數(shù),其中:αk表示第k個特征值實(shí)部小于0且大于α0的特征值實(shí)部,μk為相應(yīng)的權(quán)重系數(shù),可由μk=L″(αk-α0)2計算得到,α0=-50,同樣的L″>0.

    3.2 回溯搜索算法

    BSA是Civicioglu.P于2013年提出的一種基于種群的優(yōu)化算法,具有良好的全局收斂和收斂速度快的優(yōu)點(diǎn)[16].由于光伏發(fā)電并網(wǎng)系統(tǒng)控制參數(shù)之間具有強(qiáng)烈的耦合性,采用BSA對其進(jìn)行尋優(yōu)能夠更好地得到全局最優(yōu)解.

    BSA通過隨機(jī)初始化種群,用交叉變異的方法和目標(biāo)函數(shù)評價種群的適應(yīng)度,兩次選擇粒子,逐步迭代,直至收斂.

    BSA的尋優(yōu)過程分為種群初始化、第Ⅰ次選擇、變異、交叉、越界處理和第Ⅱ次選擇6個步驟.

    (1)種群初始化

    用粒子的維度(1~6)分別表示6個系統(tǒng)參數(shù)KP1、KI1、KP2、KI2、KP3和 KI3.粒子的隨機(jī)初始化表示為

    式中:行數(shù)i表示種群中粒子個數(shù);列數(shù)j表示每個粒子的維度數(shù);dup(j)和dlow(j)分別表示第j個維度的上下界,可通過上一節(jié)的結(jié)論確定.

    (2)第Ⅰ次選擇

    第Ⅰ次選擇用當(dāng)前種群隨機(jī)更新歷史種群,并對更新后的歷史種群中的個體進(jìn)行隨機(jī)排序,其更新方式為

    (3)變異

    通過種群中個體變異來產(chǎn)生實(shí)驗種群的初級種群,其變異方式如式(11),

    式中:Pmutant為變異后種群;F為變異尺度系數(shù)控制個體在各方向變異的幅度.

    式中:Ofit表示目標(biāo)函數(shù).

    3.3 基于BSA的并網(wǎng)發(fā)電系統(tǒng)控制器參數(shù)優(yōu)化本文采用回溯搜索法實(shí)現(xiàn)光伏并網(wǎng)發(fā)電系統(tǒng)控制參數(shù)協(xié)調(diào)優(yōu)化,其流程圖如圖4所示.

    (4)交叉

    BSA算法的交叉步驟是基于兩種交叉策略等概率隨機(jī)調(diào)用的交叉策略.隨機(jī)產(chǎn)生一個規(guī)模為i×j的二進(jìn)制整數(shù)編碼的矩陣Cmap,用于指導(dǎo)試驗種群Ti,j是否受原種群Pi,j支配.當(dāng)Cmap(i,j)=1時,用原種群中Pi,j更新Ti,j,

    (5)越界處理

    考慮到目標(biāo)函數(shù)總是在邊界處取得極值,在判斷粒子越界后,采用的是隨機(jī)初始化該粒子的方法,如式(13)所示,

    (6)第Ⅱ次選擇

    第Ⅱ次選擇的粒子是通過比較變異交叉后的粒子和歷史種群粒子的適應(yīng)度得到的,選擇適應(yīng)度較好的粒子為該次迭代的最優(yōu)粒子,如式(14)所示,

    圖4 BSA算法流程圖Fig.4 Flow chart of backtracking search algorithm

    在BSA算法中,具體設(shè)置如下:粒子維度為6,分別代表控制參數(shù)KP1、KP2、KP3、KI1、KI2和KI3,各個粒子的上下界分別根據(jù)表2得到,種群規(guī)模為50,迭代次數(shù)為100;變異因子F選定為服從正太分布的隨機(jī)數(shù),優(yōu)化目標(biāo)函數(shù)根據(jù)式(7)得到.

    根據(jù)表2所示的控制器參數(shù)取值范圍,采用BSA優(yōu)化算法實(shí)現(xiàn)控制參數(shù)的協(xié)調(diào)優(yōu)化.控制器參數(shù)優(yōu)化結(jié)果如表3所示,優(yōu)化前后系統(tǒng)特征值和阻尼比如表4所示.優(yōu)化結(jié)果表明:在控制參數(shù)優(yōu)化前,系統(tǒng)特征根λ1和λ2的阻尼比較小,且λ3的實(shí)部絕對值較小,容易引起系統(tǒng)的不穩(wěn)定和振蕩;在控制參數(shù)優(yōu)化后,系統(tǒng)特征值均滿足優(yōu)化目標(biāo),系統(tǒng)的穩(wěn)定性較高.

    表3 優(yōu)化前后的控制參數(shù)Tab.3 Control parameters before and after optimization

    表4 優(yōu)化前后系統(tǒng)特征值和阻尼比Tab.4 Eigenvalues and damping ratios before and after optimization

    4 仿真測試

    為了驗證本文所提出方法的有效性,根據(jù)圖1所示的光伏并網(wǎng)發(fā)電系統(tǒng),分別對外界因素影響下系統(tǒng)輸出功率突變情況和電力系統(tǒng)端電壓下降情況,進(jìn)行了優(yōu)化前后的仿真測試及分析.

    (1)功率突變

    光伏電池在運(yùn)行過程中,受到外界因素的影響會導(dǎo)致輸出功率的改變.優(yōu)化前后系統(tǒng)輸出功率的動態(tài)響應(yīng)曲線如圖5所示,系統(tǒng)在1 s時刻由10 kW提高至13 kW.仿真結(jié)果表明:在參數(shù)優(yōu)化前,系統(tǒng)在功率突變情況下,需要大約0.2s時間來重新恢復(fù)到穩(wěn)定運(yùn)行狀態(tài);而在參數(shù)優(yōu)化后,從功率突變時刻達(dá)到新的穩(wěn)態(tài)運(yùn)行點(diǎn)的時間明顯減短,系統(tǒng)具有較好的動態(tài)響應(yīng).

    圖5 功率突變下系統(tǒng)動態(tài)響應(yīng)Fig.5 Dynamic response when power is in step change

    優(yōu)化前后光伏發(fā)電并網(wǎng)系統(tǒng)在輸出功率為10 kW的穩(wěn)態(tài)運(yùn)行條件下,并網(wǎng)側(cè)電流諧波含量如表5所示.以A相電流為例,在進(jìn)行系統(tǒng)參數(shù)優(yōu)化前,其諧波總含量為4.08%,小于IEEE Std1547—2003標(biāo)準(zhǔn)5%的要求,符合并網(wǎng)條件;在進(jìn)行系統(tǒng)參數(shù)優(yōu)化后,A相電流諧波總含量下降為2.62%.數(shù)據(jù)表明:對系統(tǒng)控制參數(shù)的優(yōu)化在提高系統(tǒng)動態(tài)響應(yīng)的同時,還降低了并網(wǎng)電流的諧波含量,提高了系統(tǒng)供電質(zhì)量.

    表5 優(yōu)化前后電流諧波含量Tab.5 Harmonic content of current before and after optimization %

    (2)三相電壓10%跌落

    電力系統(tǒng)三相電壓跌落情況下的仿真結(jié)果如圖6所示.系統(tǒng)在1 s時刻發(fā)生10%電壓跌落,仿真結(jié)果表明:參數(shù)優(yōu)化前,在三相電壓跌落時刻,系統(tǒng)輸出的有功功率和無功功率均有一個持續(xù)時間約為2 s的震蕩狀態(tài),而在系統(tǒng)參數(shù)優(yōu)化后,即使在三相電壓發(fā)生跌落的情況下,輸出的有功和無功功率均比較平滑,系統(tǒng)具有較好的小干擾穩(wěn)定性.

    圖6 三相電壓跌落下系統(tǒng)動態(tài)響應(yīng)Fig.6 Dynamic response under the situation of three-phase voltage step-down

    當(dāng)電力系統(tǒng)側(cè)發(fā)生三相電壓跌落,光伏發(fā)電并網(wǎng)系統(tǒng)恢復(fù)至穩(wěn)態(tài)運(yùn)行后,通過對網(wǎng)側(cè)三相電流的快速傅里葉分析得到優(yōu)化前的電流諧波總含量見表5.在進(jìn)行系統(tǒng)參數(shù)優(yōu)化前,三相諧波總含量分別為6.34%、5.99%、6.02%;參數(shù)優(yōu)化后三相諧波總含量分別為2.65%、2.83%、2.72%.數(shù)據(jù)表明:優(yōu)化前系統(tǒng)在電壓突然跌落后,網(wǎng)側(cè)電流諧波含量已大于5%,不滿足并網(wǎng)條件;而優(yōu)化后的系統(tǒng)在該情況下仍舊符合并網(wǎng)要求.因此,所提出的系統(tǒng)參數(shù)優(yōu)化方法能夠在系統(tǒng)受到小干擾后依舊保持較好的穩(wěn)定運(yùn)行狀態(tài).

    5 結(jié) 論

    本文首先通過建立包含光伏電池、并網(wǎng)逆變器、變壓器和電網(wǎng)的光伏發(fā)電系統(tǒng)小信號數(shù)學(xué)模型,準(zhǔn)確和方便地得到系統(tǒng)的特征根;接著通過特征值分析表明光伏發(fā)電系統(tǒng)參數(shù)在一定程度上影響著系統(tǒng)運(yùn)行;再由根軌跡圖得到各個特征值穩(wěn)定情況下的系統(tǒng)參數(shù)邊界值,并采用基于BSA算法的協(xié)調(diào)優(yōu)化方法對計及小干擾穩(wěn)定性、阻尼比和穩(wěn)定裕度的優(yōu)化目標(biāo)函數(shù)進(jìn)行尋優(yōu).通過比較參數(shù)優(yōu)化前后,系統(tǒng)特征值和阻尼比,以及在功率突變情況和三相電壓跌落條件下的仿真測試結(jié)果,對比表明本文所提出的采用基于BSA優(yōu)化算法,建立協(xié)調(diào)優(yōu)化目標(biāo)函數(shù),用特征值根軌跡圖確定參數(shù)邊界的方法能夠提高系統(tǒng)的動態(tài)響應(yīng)和小干擾穩(wěn)定性.

    [1] 顧偉,吳志,王銳.考慮污染氣體排放的熱電連供微電網(wǎng)多目標(biāo)無功優(yōu)化[J].電力系統(tǒng)自動化,2012,36(14):177-185.

    GUWei,WU Zhi,WANG Rui.Multi-objective optimization of combined heat and power microgrid considering pollutant emission[J]. Automatics of Electric Power Systems,2012,36(14):177-185.

    [2] 李生虎,劉正楷,楊振林.風(fēng)電系統(tǒng)中異步電機(jī)機(jī)電模式研究[J].電力系統(tǒng)保護(hù)與控制,2010,38(21):6-11.

    LI Shenghu, LIU Zhengkai, YANG Zhenlin. Electronmechanical modes of induction generators in wind power system[J].Power System Protection and Control,2010,38(21):6-11.

    [3] 李奇,陳維榮,劉述奎,等.基于虛擬磁鏈的質(zhì)子交換膜燃料電池并網(wǎng)控制系統(tǒng)研究[J].電力系統(tǒng)保護(hù)與控制,2013,41(2):91-96.

    LI Qi,CHEN Weirong,LIU Shukui,et al.Study on proton exchange membrane fuel cell grid-connected control system based on virtual flux-linage[J].Power System Protection and Control,2013,41(2):91-96.

    [4] 韓瑩,陳維榮,李奇,等.分布式電源并網(wǎng)逆變器諧波抑制方法[J].電力系統(tǒng)及其自動化學(xué)報,2014,26(9):1-6.

    HAN Ying,CHEN Weirong,LI Qi,LIU Zhixiang.The method of harmonic suppression on distributed generator grid-inverter[J].Automatics of Electric Power Systems,2014,26(9):1-6.

    [5] 李奇,陳維榮,劉述奎,等.多變量H∞次優(yōu)控制在質(zhì)子交換膜燃料電池壓力控制系統(tǒng)中的應(yīng)用[J].中國電機(jī)工程學(xué)報,2010,30(20):123-128.

    LI Qi,CHEN Weirong,LIU Shukui,et al.Application of multivariable H∞suboptimalcontrolforproton exchange membrane fuel cell pressure control system[J].Proceeding of the CSEE,2010,30(20):123-128.

    [6] LI Qi,CHENn Weirong,LIU Zhixiang,et al.Net power control based on linear matrix inequality for proton exchange membrane fuel cell system[J].IEEE Transactions on Energy Conversion,2014,19(1):1-8.

    [7] 丁磊,潘貞存,王賓.分散電源并網(wǎng)對供電可靠性的影響分析[J].電力系統(tǒng)自動化,2007,31(20):89-93.

    DING Lei,PAN Zhenxiang,WANG Bin.Impact on reliability of distribution networks with dispersed generators[J].Automatics of Electric Power Systems,2007,31(20):89-93.

    [8] 丁磊,潘貞存,叢偉.基于有根樹的分布式發(fā)電孤島搜索[J].中國電機(jī)工程學(xué)報,2008,28(25):62-67.

    DING Lei,PAN Zhencun,CONG Wei.Searching for international islanding strategies of distributed generation based on rooted tree[J].Proceedings of the CSEE,2008,28(25):62-67.

    [9] DU W,WANG H F,DUNN R.Power system smallsignal oscillation stability as affected by large-scale PV penetration[C]∥ SUPERGEN'09 International Conference on Sustainable Power Generation and Supply.Nanjing:IEEE,2009:1-6.

    [10] FIGUERES E,GARCERá S J.Sensitivity study of the dynamics of three-phase photovoltaic inverters with an LCL grid filter[J].IEEE Transactions on Industrial Electronics,2009,56(3):706-717.

    [11] HUANG Hanqi,MAO Chengxiong,LU Jiming,et al. Small-signal modeling and analysis of grid-connected photovoltaic generation systems[J].Proceedings of the CSEE,2012,32(22):7-14.

    [12] 杜威,姜齊榮,洪蘆誠,等.光伏微電網(wǎng)孤島運(yùn)行時多種振蕩模式的小信號建模分析[J].電力系統(tǒng)自動化,2014,38(10):17-23.

    DU Wei,JIANG Qirong,HONG Lucheng,et al. Small signal analysis of multiple oscillation modes in a photovoltaicmicrogrid[J]. Automation ofElectric Power System,2014,38(10):17-23.

    [13] 熊連松,劉小康,卓放,等.光伏發(fā)電系統(tǒng)的小信號建模及其控制器參數(shù)的全局優(yōu)化設(shè)計方法[J].電網(wǎng)技術(shù),2014,38(5):1234-1241.

    XIONG Liansong,LIU Xiaokang,ZHUO Fang,et al. Small-signal modeling of photovoltaic power generation system and global optimal design for its controller parameters[J]. PowerSystem Technology,2014,38(5):1234-1241.

    [14] 蘇建徽,余世杰,趙為,等.硅太陽電池工程用數(shù)學(xué)模型[J].太陽能學(xué)報,2001,22(4):409-412.

    SU Jianhui, YU Shijie, ZHAO Wei, etal. Investigation on engineering analytical model of silicon solar cells[J].Acta Energiae Solaris Sinica,2001,22(4):409-412.

    [15] WANG Chengshan,LIYan,PENGKe,etal. Coordinated optimal design of inverter controllers in a micro-grid with multiple distributed generation units[J].IEEE Transactionson PowerSystems,2013,28(3):2679-2687.

    [16] CIVICIOGLU P. Backtrackingsearch optimization algorithm for numerical optimization problems[J]. Applied Mathematics and Computation, 2013,219(15):8121-8144.

    [17] 張興,張崇.PWM整流器及其控制[M].北京:機(jī)械工業(yè)出版社,2012:1-30.

    [18] LASSETER R H,ETO J H,SCHENKMAN B,et al. CERTS microgrid laboratory testbed[J]. IEEE Transactions on Power Delivery,2011,26(1):325-332.

    [19] GUO Li,F(xiàn)ENG Yibin,LI Xialin,et al.Stability analysis of a DC microgrid with master-slave control structure[C]∥ Energy Conversion Congress and Exposition(ECCE).[S.l.]:IEEE,2014:5682-5689.

    [20] 劉夢璇,王成山,郭力,等.基于多目標(biāo)的獨(dú)立微電網(wǎng)優(yōu)化設(shè)計方法[J].電力系統(tǒng)自動化,2012,36(17):34-39.

    LIU Mengxuan,WANG Chenshan,GUO Li,et al.An optimal design method of multi objective based island microgird[J].Automatics of Electric Power Systems,2010,36(17):34-39.

    [21] 洪博文,郭力,王成山,等.微電網(wǎng)多目標(biāo)動態(tài)優(yōu)化調(diào)度模型與方法[J].電力自動化設(shè)備,2013,33(3):100-107.

    HONG Bowen,GUO Li,WANG Chenshan,et al. Model and method of dynamic multi-objective optimal dispatch for microgird[J].Electric Power Automation Equipment,2013,33(3):100-107.

    光伏并網(wǎng)發(fā)電系統(tǒng)參數(shù)協(xié)調(diào)優(yōu)化方法研究

    李 奇, 楊寒卿, 韓 瑩, 陳維榮

    Method of Parameter Coordination Optimization for Grid-Connected Photovoltaic System

    LI Qi, YANG Hanqing, HAN Ying, CHEN Weirong
    (School of Electrical Engineering,Southwest Jiaotong University,Chengdu 610031,China)

    The parameter design of grid-connected photovoltaic(PV)system is significant for its stable operation.In this work,a method of parameter coordination optimization was proposed for gridconnected PV system,by establishing the small-signal mathematical model of PV system,including PV cell,inverter,transformer and gird.By analyzing eigenvalues(root locus plot)in stable conditions,the boundary values of system parameters are available.Then the objective function for coordination optimization was built involving small-signal stability,damping ratio and stability margin.And the backtracking search algorithm(BSA)was adopted to optimize the system parameters.The simulation test is carried out to compare the dynamic system responses before and after parameter optimization. Simulation test results show that when system suffers small disturbances,the optimized grid-connected PV system can reach a new stable point within 0.1 s,and its dynamic response is improved 50%than that of the original system.

    grid-connected photovoltaic(PV)system;small-signal model;backtracking search algorithm;coordination optimization objective function;eigenvalue analysis

    李奇,楊寒卿,韓瑩,等.光伏并網(wǎng)發(fā)電系統(tǒng)參數(shù)協(xié)調(diào)優(yōu)化方法研究[J].西南交通大學(xué)學(xué)報,2016,51(5):894-901.

    0258-2724(2016)05-0894-08

    10.3969/j.issn.0258-2724.2016.05.011

    TM911.4

    A

    2015-07-22

    國家自然科學(xué)基金資助項目(61473238,51407146);四川省杰出青年基金資助項目(2015JQO016);高等學(xué)校博士學(xué)科點(diǎn)專項科研基金資助項目(20120184120011)

    李奇(1984—),男,副教授,博士生導(dǎo)師,研究方向為分布式發(fā)電并網(wǎng)技術(shù)、電力系統(tǒng)無功優(yōu)化、智能信息處理,E-mail:liqi0800@163.com

    (中文編輯:唐 晴 英文編輯:周 堯)

    猜你喜歡
    特征值種群發(fā)電
    邢氏水蕨成功繁衍并建立種群 等
    山西省發(fā)現(xiàn)刺五加種群分布
    “發(fā)電”
    一類帶強(qiáng)制位勢的p-Laplace特征值問題
    檸檬亦能發(fā)電?
    單圈圖關(guān)聯(lián)矩陣的特征值
    搖晃發(fā)電小圓球
    摩擦發(fā)電
    學(xué)與玩(2017年4期)2017-02-16 07:05:40
    基于商奇異值分解的一類二次特征值反問題
    關(guān)于兩個M-矩陣Hadamard積的特征值的新估計
    日韩中文字幕欧美一区二区 | 人人妻人人添人人爽欧美一区卜| 啦啦啦在线观看免费高清www| 亚洲伊人久久精品综合| 观看av在线不卡| 两性夫妻黄色片| 国产精品久久久久成人av| 97精品久久久久久久久久精品| 国产欧美日韩一区二区三区在线| 18禁观看日本| 亚洲国产精品国产精品| 一个人免费看片子| 亚洲成av片中文字幕在线观看 | 欧美日韩成人在线一区二区| 国产精品无大码| av免费在线看不卡| 一区二区三区精品91| 国产成人精品婷婷| 我的亚洲天堂| 欧美 亚洲 国产 日韩一| 亚洲国产成人一精品久久久| 欧美日韩亚洲高清精品| 亚洲欧美色中文字幕在线| 欧美激情 高清一区二区三区| 青春草视频在线免费观看| 成人免费观看视频高清| 黄色一级大片看看| 欧美亚洲 丝袜 人妻 在线| 看非洲黑人一级黄片| 久久久精品国产亚洲av高清涩受| 久久久久久伊人网av| 久久久精品国产亚洲av高清涩受| 久久久精品国产亚洲av高清涩受| 中文字幕人妻熟女乱码| 十八禁高潮呻吟视频| 国产乱来视频区| 亚洲精品日本国产第一区| 中国国产av一级| 国产精品亚洲av一区麻豆 | 91国产中文字幕| 黄片无遮挡物在线观看| 成人毛片a级毛片在线播放| 亚洲精品日韩在线中文字幕| 国产一区二区三区av在线| 男人操女人黄网站| 99热网站在线观看| 国产精品亚洲av一区麻豆 | 亚洲欧洲日产国产| av在线播放精品| 春色校园在线视频观看| 999精品在线视频| 美女中出高潮动态图| 女人久久www免费人成看片| 亚洲国产精品一区二区三区在线| 亚洲经典国产精华液单| 视频在线观看一区二区三区| 中文字幕色久视频| 中文字幕最新亚洲高清| 中文乱码字字幕精品一区二区三区| 亚洲av免费高清在线观看| 校园人妻丝袜中文字幕| 国语对白做爰xxxⅹ性视频网站| 天天躁日日躁夜夜躁夜夜| 热re99久久国产66热| 国产免费一区二区三区四区乱码| 亚洲婷婷狠狠爱综合网| 亚洲精品国产av成人精品| 中文乱码字字幕精品一区二区三区| 久久精品夜色国产| 欧美日韩亚洲高清精品| 天天影视国产精品| 不卡视频在线观看欧美| 欧美最新免费一区二区三区| 午夜福利在线观看免费完整高清在| 午夜日本视频在线| 色婷婷av一区二区三区视频| 人人澡人人妻人| 99热全是精品| 亚洲精华国产精华液的使用体验| 大话2 男鬼变身卡| 国产精品秋霞免费鲁丝片| 日韩三级伦理在线观看| 国产高清国产精品国产三级| 国产免费一区二区三区四区乱码| 一边摸一边做爽爽视频免费| 午夜福利乱码中文字幕| 午夜福利在线观看免费完整高清在| 91aial.com中文字幕在线观看| 日韩av免费高清视频| 一本大道久久a久久精品| 精品人妻一区二区三区麻豆| 不卡视频在线观看欧美| 国产精品成人在线| 亚洲欧洲日产国产| av女优亚洲男人天堂| 波多野结衣av一区二区av| av有码第一页| 一级片免费观看大全| 国产国语露脸激情在线看| 综合色丁香网| 在线免费观看不下载黄p国产| 久久国产精品男人的天堂亚洲| 人人妻人人澡人人爽人人夜夜| 亚洲精品,欧美精品| 欧美国产精品一级二级三级| 久久久久久久久免费视频了| 最近2019中文字幕mv第一页| 精品少妇一区二区三区视频日本电影 | 99久久人妻综合| 亚洲欧美色中文字幕在线| 18+在线观看网站| 免费大片黄手机在线观看| 日本-黄色视频高清免费观看| 亚洲精品自拍成人| 免费播放大片免费观看视频在线观看| 中文精品一卡2卡3卡4更新| av网站在线播放免费| 最近中文字幕高清免费大全6| 春色校园在线视频观看| 三上悠亚av全集在线观看| 久久久欧美国产精品| 电影成人av| 在线观看国产h片| 国产成人欧美| 午夜福利在线观看免费完整高清在| 夫妻午夜视频| 三上悠亚av全集在线观看| 国产av码专区亚洲av| 777米奇影视久久| 一级a爱视频在线免费观看| 黄片小视频在线播放| 一区在线观看完整版| 热re99久久国产66热| 亚洲欧美色中文字幕在线| 亚洲精品aⅴ在线观看| 成人影院久久| 在线观看人妻少妇| 欧美日韩成人在线一区二区| 性色avwww在线观看| 欧美97在线视频| 视频在线观看一区二区三区| 国产免费一区二区三区四区乱码| 免费看不卡的av| 男人操女人黄网站| 女性被躁到高潮视频| 久久毛片免费看一区二区三区| av免费观看日本| 精品久久久精品久久久| 欧美黄色片欧美黄色片| 日本wwww免费看| 免费在线观看视频国产中文字幕亚洲 | 国产麻豆69| 亚洲一区中文字幕在线| 免费大片黄手机在线观看| 十分钟在线观看高清视频www| 丝袜脚勾引网站| 中文字幕人妻丝袜一区二区 | 午夜91福利影院| 成年女人在线观看亚洲视频| 人人妻人人澡人人看| www日本在线高清视频| 亚洲视频免费观看视频| av片东京热男人的天堂| 一级a爱视频在线免费观看| 看免费成人av毛片| 黑人猛操日本美女一级片| 晚上一个人看的免费电影| 熟女电影av网| 日本爱情动作片www.在线观看| 日日啪夜夜爽| 街头女战士在线观看网站| 久久午夜福利片| 成年女人在线观看亚洲视频| 最新中文字幕久久久久| 精品少妇一区二区三区视频日本电影 | 国产av一区二区精品久久| 国产 一区精品| 人妻人人澡人人爽人人| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 侵犯人妻中文字幕一二三四区| 国产视频首页在线观看| 久久久久久久大尺度免费视频| 天堂中文最新版在线下载| 国产精品欧美亚洲77777| 国产精品久久久久久久久免| 天天躁夜夜躁狠狠久久av| 亚洲精品久久午夜乱码| 久久久国产欧美日韩av| 99香蕉大伊视频| 春色校园在线视频观看| 黄色一级大片看看| 久久久久人妻精品一区果冻| 国产亚洲av片在线观看秒播厂| 亚洲av国产av综合av卡| 一本大道久久a久久精品| 久久亚洲国产成人精品v| 搡女人真爽免费视频火全软件| 国产老妇伦熟女老妇高清| 不卡视频在线观看欧美| 男女高潮啪啪啪动态图| 男人舔女人的私密视频| 18在线观看网站| 亚洲婷婷狠狠爱综合网| 亚洲精品美女久久久久99蜜臀 | 高清黄色对白视频在线免费看| 在线观看免费日韩欧美大片| 精品国产一区二区三区四区第35| 国产有黄有色有爽视频| 日韩大片免费观看网站| 精品久久久精品久久久| 国产一级毛片在线| 99久国产av精品国产电影| 久久精品国产鲁丝片午夜精品| 亚洲欧美一区二区三区黑人 | 另类精品久久| 精品人妻熟女毛片av久久网站| 日韩伦理黄色片| 夜夜骑夜夜射夜夜干| 国产亚洲精品第一综合不卡| 亚洲成色77777| 国产麻豆69| 99久国产av精品国产电影| 人妻 亚洲 视频| av在线观看视频网站免费| 欧美日韩视频高清一区二区三区二| 午夜激情av网站| 五月天丁香电影| 亚洲国产色片| 天天躁夜夜躁狠狠躁躁| 亚洲欧美成人精品一区二区| 母亲3免费完整高清在线观看 | 国产成人a∨麻豆精品| 麻豆精品久久久久久蜜桃| 美女高潮到喷水免费观看| 一级毛片 在线播放| av国产久精品久网站免费入址| 天天影视国产精品| 在线天堂最新版资源| 国产高清国产精品国产三级| 丝瓜视频免费看黄片| 天天操日日干夜夜撸| 9色porny在线观看| 婷婷成人精品国产| 国产av精品麻豆| 欧美国产精品一级二级三级| 妹子高潮喷水视频| 久久久久久久久久久免费av| 国产一区二区三区av在线| 高清不卡的av网站| 美女视频免费永久观看网站| 成人亚洲精品一区在线观看| 国产一区亚洲一区在线观看| 性色av一级| 性高湖久久久久久久久免费观看| 少妇人妻 视频| 国产精品久久久久久av不卡| 午夜精品国产一区二区电影| 肉色欧美久久久久久久蜜桃| 亚洲欧洲日产国产| 另类精品久久| 国产精品.久久久| 最新中文字幕久久久久| 看免费av毛片| 国产深夜福利视频在线观看| 男女免费视频国产| 9色porny在线观看| 丝袜喷水一区| 成人毛片a级毛片在线播放| 国产精品 欧美亚洲| 亚洲欧美精品综合一区二区三区 | 亚洲精品美女久久av网站| 亚洲一区二区三区欧美精品| 成人国语在线视频| 人妻人人澡人人爽人人| 少妇的逼水好多| 婷婷色综合大香蕉| av国产精品久久久久影院| 丰满饥渴人妻一区二区三| 免费在线观看完整版高清| 精品一区二区三区四区五区乱码 | h视频一区二区三区| 亚洲av日韩在线播放| av天堂久久9| 狠狠精品人妻久久久久久综合| 90打野战视频偷拍视频| 少妇的丰满在线观看| 女性被躁到高潮视频| 80岁老熟妇乱子伦牲交| 男人操女人黄网站| 成人18禁高潮啪啪吃奶动态图| 建设人人有责人人尽责人人享有的| 久久久欧美国产精品| 精品国产乱码久久久久久小说| 丰满少妇做爰视频| 久久精品熟女亚洲av麻豆精品| 亚洲精品,欧美精品| 久久国产亚洲av麻豆专区| tube8黄色片| 日韩一本色道免费dvd| 久久青草综合色| 极品少妇高潮喷水抽搐| 人成视频在线观看免费观看| 两个人看的免费小视频| 丝袜在线中文字幕| 叶爱在线成人免费视频播放| 一级片'在线观看视频| 啦啦啦啦在线视频资源| 免费观看av网站的网址| 成人影院久久| 韩国精品一区二区三区| a级毛片在线看网站| 日韩大片免费观看网站| 大陆偷拍与自拍| 成人漫画全彩无遮挡| 黄片无遮挡物在线观看| 国产日韩一区二区三区精品不卡| 我要看黄色一级片免费的| 国产伦理片在线播放av一区| 久久久久视频综合| 蜜桃国产av成人99| 永久网站在线| 男女高潮啪啪啪动态图| 国产在线一区二区三区精| 国产有黄有色有爽视频| 日韩大片免费观看网站| 麻豆乱淫一区二区| 欧美最新免费一区二区三区| 大片电影免费在线观看免费| av在线app专区| 视频区图区小说| 一级a爱视频在线免费观看| 国产免费福利视频在线观看| 欧美国产精品va在线观看不卡| 国产熟女欧美一区二区| 欧美人与性动交α欧美软件| 天堂中文最新版在线下载| 性色avwww在线观看| 亚洲av电影在线进入| 卡戴珊不雅视频在线播放| 久久久久国产网址| 国产成人a∨麻豆精品| 大香蕉久久网| 在线观看免费高清a一片| 捣出白浆h1v1| 亚洲精品久久午夜乱码| 亚洲中文av在线| 黄色怎么调成土黄色| 黄片小视频在线播放| 国产亚洲最大av| 两个人免费观看高清视频| 一级片'在线观看视频| 久久国产精品男人的天堂亚洲| 各种免费的搞黄视频| 丝袜喷水一区| 啦啦啦在线观看免费高清www| 午夜影院在线不卡| 自拍欧美九色日韩亚洲蝌蚪91| 久久久国产欧美日韩av| 啦啦啦在线观看免费高清www| 午夜影院在线不卡| 午夜久久久在线观看| 美女中出高潮动态图| 久久久国产一区二区| 男女无遮挡免费网站观看| 亚洲国产欧美在线一区| 久久婷婷青草| 欧美少妇被猛烈插入视频| h视频一区二区三区| 日本猛色少妇xxxxx猛交久久| h视频一区二区三区| 男女边吃奶边做爰视频| a级片在线免费高清观看视频| 欧美精品av麻豆av| 99re6热这里在线精品视频| 青青草视频在线视频观看| 亚洲国产精品999| 青青草视频在线视频观看| 91精品三级在线观看| 在线观看美女被高潮喷水网站| 高清视频免费观看一区二区| 高清不卡的av网站| 成人亚洲欧美一区二区av| 一本久久精品| 交换朋友夫妻互换小说| 少妇熟女欧美另类| 国产精品无大码| 久热久热在线精品观看| 亚洲欧洲国产日韩| 丝袜喷水一区| 男女边吃奶边做爰视频| 亚洲第一青青草原| 高清不卡的av网站| 美女脱内裤让男人舔精品视频| 亚洲国产最新在线播放| 最近最新中文字幕大全免费视频 | 十分钟在线观看高清视频www| 精品亚洲成国产av| 亚洲综合色网址| 青春草视频在线免费观看| 男女边吃奶边做爰视频| 高清视频免费观看一区二区| 国产1区2区3区精品| 久久久久久久久久人人人人人人| 久久国内精品自在自线图片| 亚洲欧美精品自产自拍| 国产精品免费视频内射| 黄片无遮挡物在线观看| 韩国精品一区二区三区| 日韩中文字幕视频在线看片| 赤兔流量卡办理| 精品第一国产精品| 久久久久国产网址| 精品久久久精品久久久| 欧美精品一区二区大全| 精品少妇一区二区三区视频日本电影 | 黄色配什么色好看| 丰满少妇做爰视频| 交换朋友夫妻互换小说| 性色av一级| 看非洲黑人一级黄片| 九九爱精品视频在线观看| 色婷婷av一区二区三区视频| 国产又爽黄色视频| 国产人伦9x9x在线观看 | 日韩免费高清中文字幕av| 国产成人精品久久久久久| 日韩免费高清中文字幕av| 午夜福利网站1000一区二区三区| 少妇熟女欧美另类| 国产极品粉嫩免费观看在线| 久久久久久久亚洲中文字幕| 亚洲欧洲精品一区二区精品久久久 | 国产乱人偷精品视频| 一本色道久久久久久精品综合| 国产精品香港三级国产av潘金莲 | 国产精品一国产av| 好男人视频免费观看在线| 在线天堂最新版资源| 曰老女人黄片| 青草久久国产| 乱人伦中国视频| 亚洲成国产人片在线观看| 久久婷婷青草| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 亚洲av免费高清在线观看| 男女边摸边吃奶| 另类精品久久| 美女国产视频在线观看| 日韩精品免费视频一区二区三区| 美女xxoo啪啪120秒动态图| 亚洲精品美女久久av网站| 9191精品国产免费久久| 交换朋友夫妻互换小说| 亚洲精品中文字幕在线视频| 国产成人免费观看mmmm| 久久这里有精品视频免费| 9色porny在线观看| 久久久久久久久久人人人人人人| 人妻少妇偷人精品九色| 亚洲美女黄色视频免费看| 18禁观看日本| 午夜激情久久久久久久| 黄片小视频在线播放| 在线观看免费高清a一片| 成年美女黄网站色视频大全免费| 国产精品免费视频内射| 日韩伦理黄色片| 老女人水多毛片| 欧美少妇被猛烈插入视频| 老汉色av国产亚洲站长工具| 天天躁日日躁夜夜躁夜夜| 欧美精品一区二区大全| 如何舔出高潮| 日韩三级伦理在线观看| 国产伦理片在线播放av一区| av国产久精品久网站免费入址| 午夜福利网站1000一区二区三区| 亚洲av综合色区一区| 久久久久久伊人网av| 午夜福利乱码中文字幕| 一区二区三区四区激情视频| 天堂中文最新版在线下载| 国产日韩欧美视频二区| 免费播放大片免费观看视频在线观看| 国产在线一区二区三区精| 国产一区二区在线观看av| 国产在线一区二区三区精| 国产精品一国产av| 啦啦啦啦在线视频资源| 黄片无遮挡物在线观看| av电影中文网址| 大陆偷拍与自拍| 夜夜骑夜夜射夜夜干| 精品国产露脸久久av麻豆| 午夜91福利影院| 黄色配什么色好看| 国产一区二区 视频在线| 免费少妇av软件| 日韩三级伦理在线观看| 少妇精品久久久久久久| 成人国产av品久久久| 日本欧美国产在线视频| 久久久国产精品麻豆| av在线播放精品| 国产精品二区激情视频| 69精品国产乱码久久久| 精品国产乱码久久久久久小说| 久久午夜福利片| 久久久国产精品麻豆| 亚洲成人手机| 国产精品.久久久| 91在线精品国自产拍蜜月| 久久精品国产亚洲av涩爱| 黄片播放在线免费| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 人人妻人人爽人人添夜夜欢视频| 人人妻人人添人人爽欧美一区卜| 亚洲在久久综合| 成人漫画全彩无遮挡| 亚洲久久久国产精品| 校园人妻丝袜中文字幕| 欧美日本中文国产一区发布| 久久97久久精品| 久久精品人人爽人人爽视色| 王馨瑶露胸无遮挡在线观看| 精品国产露脸久久av麻豆| 亚洲国产精品一区三区| 亚洲内射少妇av| 亚洲第一av免费看| 国产亚洲最大av| 熟女av电影| 大片免费播放器 马上看| 午夜福利在线观看免费完整高清在| 亚洲美女黄色视频免费看| 午夜福利影视在线免费观看| 国产有黄有色有爽视频| 久久精品国产综合久久久| 国产毛片在线视频| 亚洲精品久久久久久婷婷小说| 搡女人真爽免费视频火全软件| 又黄又粗又硬又大视频| 久久久久久久久久久免费av| 女性生殖器流出的白浆| 麻豆av在线久日| 国产精品久久久久久精品电影小说| 不卡视频在线观看欧美| 久久午夜综合久久蜜桃| 午夜福利视频在线观看免费| 制服诱惑二区| 国产日韩欧美视频二区| 一级,二级,三级黄色视频| 午夜福利视频精品| 日韩精品免费视频一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 国产精品国产三级国产专区5o| 久久久久久久久久人人人人人人| 国产成人免费无遮挡视频| 不卡视频在线观看欧美| 午夜免费鲁丝| 人人妻人人爽人人添夜夜欢视频| 777米奇影视久久| 一本大道久久a久久精品| 亚洲第一青青草原| 国产激情久久老熟女| 秋霞伦理黄片| 一级毛片 在线播放| 纯流量卡能插随身wifi吗| 午夜av观看不卡| 蜜桃国产av成人99| 久久99一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产色片| 国产精品一区二区在线观看99| 丰满迷人的少妇在线观看| 亚洲,欧美精品.| 精品少妇久久久久久888优播| 伦精品一区二区三区| 少妇的丰满在线观看| 一个人免费看片子| 一区在线观看完整版| 又黄又粗又硬又大视频| 国产av一区二区精品久久| 亚洲天堂av无毛| 亚洲av免费高清在线观看| 黑丝袜美女国产一区| 免费观看无遮挡的男女| 久久久久国产网址| 青青草视频在线视频观看| 欧美变态另类bdsm刘玥| av又黄又爽大尺度在线免费看| 极品少妇高潮喷水抽搐| 成人国产av品久久久| 成人亚洲欧美一区二区av| 免费观看无遮挡的男女| 精品酒店卫生间| 侵犯人妻中文字幕一二三四区| av免费在线看不卡| 99国产精品免费福利视频| 国产精品久久久久久精品电影小说| 亚洲一码二码三码区别大吗| 一本色道久久久久久精品综合| 国产精品一二三区在线看| 男女边摸边吃奶| 国产 精品1| 男女下面插进去视频免费观看| 9191精品国产免费久久| 久久久亚洲精品成人影院| 91成人精品电影| av在线播放精品| 婷婷色麻豆天堂久久| av在线老鸭窝| 最近中文字幕高清免费大全6| 日韩av不卡免费在线播放|