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

    計(jì)及分布式電源的改進(jìn)配網(wǎng)潮流計(jì)算方法

    2016-08-16 02:30:51焦彥軍蔣晨陽(yáng)
    關(guān)鍵詞:雅克牛頓潮流

    梁 宵,焦彥軍,蔣晨陽(yáng)

    (華北電力大學(xué) 電氣與電子工程學(xué)院, 河北 保定 071003)

    ?

    計(jì)及分布式電源的改進(jìn)配網(wǎng)潮流計(jì)算方法

    梁宵,焦彥軍,蔣晨陽(yáng)

    (華北電力大學(xué) 電氣與電子工程學(xué)院, 河北 保定 071003)

    針對(duì)多節(jié)點(diǎn)形式分布式電源并入配電網(wǎng)后導(dǎo)致傳統(tǒng)潮流計(jì)算方法并不適用的問(wèn)題,對(duì)牛頓潮流計(jì)算方法進(jìn)行改進(jìn),提出每次迭代時(shí)結(jié)合梯度下降法與牛頓法對(duì)雅克比矩陣中各元素進(jìn)行優(yōu)化。改進(jìn)算法解決了配網(wǎng)中由于分布式電源節(jié)點(diǎn)類型豐富導(dǎo)致傳統(tǒng)算法計(jì)算困難的問(wèn)題,同時(shí),也解決了由于配網(wǎng)線路阻抗比值較大導(dǎo)致雅克比矩陣收斂困難的問(wèn)題。根據(jù)改進(jìn)算法提出以零阻抗連接編號(hào)方案為基礎(chǔ)的潮流計(jì)算流程。通過(guò)數(shù)學(xué)理論分析與IEEE33節(jié)點(diǎn)仿真驗(yàn)證,改進(jìn)算法在收斂性、收斂速度等方面較傳統(tǒng)方法均有所提高。經(jīng)由驗(yàn)證結(jié)果表明,改進(jìn)算法適用范圍較廣,對(duì)于含分布式電源的配網(wǎng)快速潮流分析與實(shí)時(shí)調(diào)度計(jì)算提供有效的參考。

    分布式電源;配網(wǎng)潮流計(jì)算;雅克比矩陣;改進(jìn)算法;零阻抗連接編號(hào)

    0 引 言

    受環(huán)境問(wèn)題、未來(lái)能源問(wèn)題和電力體制改革等綜合因素的影響,分布式電源(Distributed Generation, DG)近年成為研究熱點(diǎn)。分布式電源是指在配電系統(tǒng)靠近用戶側(cè)引入的功率為數(shù)千瓦至50 MW的小型模塊式、與環(huán)境兼容的獨(dú)立電源,具有投資少、占地小、建設(shè)周期短、節(jié)能、環(huán)保等特點(diǎn)[1-2]。分布式電源并網(wǎng)后,改變?cè)潆娋W(wǎng)單端電源輻射型網(wǎng)絡(luò)為多電源網(wǎng)絡(luò),從而改變潮流方向[3-4]。因此,亟需在保證收斂性前提下,得到計(jì)算速度快的新潮流計(jì)算方法,滿足含分布式電源配電網(wǎng)中分析和調(diào)度的需要。

    目前常見(jiàn)的分布式發(fā)電并網(wǎng)系統(tǒng)有微型燃?xì)廨啓C(jī)并網(wǎng)系統(tǒng)、風(fēng)力發(fā)電系統(tǒng)、燃料電池與光伏并網(wǎng)發(fā)電系統(tǒng)[5-8]。對(duì)于微型燃?xì)廨啓C(jī)并網(wǎng)系統(tǒng)根據(jù)并網(wǎng)PWM變流器類型分為電流型與電壓型,對(duì)于電流型并網(wǎng)系統(tǒng)其節(jié)點(diǎn)類型為PI節(jié)點(diǎn);而對(duì)于電壓型并網(wǎng)系統(tǒng),其節(jié)點(diǎn)類型由運(yùn)行模式?jīng)Q定,聯(lián)網(wǎng)模式下為PV節(jié)點(diǎn),孤島模式下則為VQ節(jié)點(diǎn)。風(fēng)力發(fā)電系統(tǒng)以異步電機(jī)直接并網(wǎng)發(fā)電系統(tǒng)和雙饋發(fā)電系統(tǒng)較為常見(jiàn),其節(jié)點(diǎn)形式與其運(yùn)行方式及控制策略有關(guān)。異步電機(jī)節(jié)點(diǎn)發(fā)出的有功功率確定,無(wú)功功率則與機(jī)端電壓有關(guān),稱該類電源節(jié)點(diǎn)為電壓靜特性節(jié)點(diǎn),表示為P-Q(V);雙饋發(fā)電系統(tǒng)一般采用恒功率因數(shù)控制方式,因此節(jié)點(diǎn)類型為PQ節(jié)點(diǎn)。對(duì)于燃料電池的輸出控制與常規(guī)發(fā)電機(jī)控制原理近似,故在潮流計(jì)算中,燃料電池發(fā)電站并網(wǎng)節(jié)點(diǎn)可作為PV節(jié)點(diǎn)處理。光伏并網(wǎng)發(fā)電系統(tǒng)的節(jié)點(diǎn)類型同樣受并網(wǎng)逆變器類型決定,對(duì)于電流型逆變器其節(jié)點(diǎn)類型為PI節(jié)點(diǎn);而電壓型節(jié)點(diǎn)類型為PV節(jié)點(diǎn)。

    傳統(tǒng)配電網(wǎng)的潮流算法有前推回代法、回路阻抗法、牛頓法和Zbus高斯法,具體算法參考文獻(xiàn)[9-12]。傳統(tǒng)潮流計(jì)算方法適用節(jié)點(diǎn)類型有限,對(duì)于分布式電源的節(jié)點(diǎn)類型并不適用。文獻(xiàn)[13]對(duì)前推回代法進(jìn)行改進(jìn),以節(jié)點(diǎn)電抗矩陣形式將PV節(jié)點(diǎn)轉(zhuǎn)化為PQ節(jié)點(diǎn),對(duì)于無(wú)功初值的確定,提出無(wú)功分?jǐn)傇?,提高了?jì)算速度。文獻(xiàn)[14]從分布式電源的數(shù)學(xué)模型出發(fā),著重分析風(fēng)機(jī)并網(wǎng)后,對(duì)雅克比矩陣的影響,引入?yún)?shù)修正雅克比矩陣。該方法可以有效解決風(fēng)機(jī)并網(wǎng)后傳統(tǒng)牛頓法的問(wèn)題,但是,由于并未考慮其他分布式電源并網(wǎng)形式,局限性較強(qiáng)。文獻(xiàn)[15]以回路分析法為基礎(chǔ),利用配電網(wǎng)相角相差不大的特點(diǎn)將回路方程化簡(jiǎn)為線性代數(shù)方程。該方法在較為簡(jiǎn)單的配電網(wǎng)中適用性較強(qiáng),但隨著網(wǎng)絡(luò)拓?fù)鋸?fù)雜度的不斷提升收斂速度明顯降低。

    多節(jié)點(diǎn)類型的分布式電源經(jīng)數(shù)次迭代后均可轉(zhuǎn)化成PQ、PV節(jié)點(diǎn)[16],目前,對(duì)于無(wú)分布式電源并網(wǎng)的配網(wǎng)廣泛應(yīng)用的前推回代潮流計(jì)算方法主要處理PQ節(jié)點(diǎn),不適用于含PV節(jié)點(diǎn)的潮流計(jì)算中,同時(shí),該方法的收斂速度隨著網(wǎng)絡(luò)規(guī)模的擴(kuò)大而降低。牛頓法潮流計(jì)算在處理PQ、PV節(jié)點(diǎn)理論上較為成熟,且迭代次數(shù)于網(wǎng)絡(luò)規(guī)模無(wú)關(guān)。影響牛頓法在配電網(wǎng)潮流計(jì)算的主要問(wèn)題是在配電網(wǎng)中電阻電抗比值較大使雅克比矩陣收斂困難,較易出現(xiàn)“病態(tài)”潮流。

    針對(duì)上述問(wèn)題,本文提出每次迭代時(shí)梯度下降法與牛頓法相結(jié)合對(duì)雅克比矩陣元素進(jìn)行調(diào)整,解決牛頓法在配網(wǎng)潮流計(jì)算中收斂困難的問(wèn)題。通過(guò)數(shù)學(xué)理論分析與IEEE33節(jié)點(diǎn)實(shí)例仿真對(duì)比改進(jìn)算法與傳統(tǒng)算法的收斂性,收斂速度等問(wèn)題,驗(yàn)證改進(jìn)算法適用于多分布式電源并網(wǎng)的配網(wǎng)潮流計(jì)算中。

    1 改進(jìn)潮流算法數(shù)學(xué)模型

    1.1數(shù)學(xué)模型

    牛頓法潮流計(jì)算的基本思路[17]是將潮流方程F(X)=0用泰勒級(jí)數(shù)展開(kāi),在初始值接近于真值的前提下,略去二階項(xiàng)以及更高次項(xiàng),將非線性方程逐次線性化,反復(fù)形成并求解修正方程。直角坐標(biāo)下修正方程為

    (1)

    式中:X和ΔX分別表示n個(gè)狀態(tài)變量和其修正量組成的n維列向量;F(X)是由n個(gè)函數(shù)組成的n維列向量;J為n階雅克比矩陣,第i行第j列元素Jij=?fi/?xj,角標(biāo)k表示第k次迭代。

    由上式可知,雅克比矩陣是潮流計(jì)算的關(guān)鍵,因此,對(duì)雅克比矩陣的修正計(jì)算是解決牛頓法收斂性問(wèn)題的根源。在初值與真值相差較大的條件下,將非線性方程逐次線性化的過(guò)程中,牛拉法迭代收斂困難,為避免該問(wèn)題的發(fā)生,利用非單調(diào)搜索準(zhǔn)則,現(xiàn)將式(1)改進(jìn)為

    (2)

    式中:μ為迭代參數(shù);I為與雅克比矩陣J同階單位陣。

    該方法主旨是將梯度下降法與牛頓法的算法優(yōu)缺點(diǎn)互補(bǔ),從而達(dá)到高效結(jié)合。從式(2)中可以看出,若迭代參數(shù)μ(k)接近于1,此時(shí)為梯度下降法,若迭代參數(shù)μ(k)=0,則轉(zhuǎn)變?yōu)榕nD法。梯度下降法在迭代過(guò)程中下降速度較快,對(duì)于收斂速度的改善有很大幫助,但在最優(yōu)值附近繼續(xù)迭代時(shí),由于下降梯度趨于零,導(dǎo)致收斂速度減緩;而牛頓法在最優(yōu)解附近依然能夠產(chǎn)生一個(gè)理想的搜索方向,彌補(bǔ)了梯度下降法的不足。通過(guò)將上述方法的有效結(jié)合,克服了雅克比矩陣J幾乎奇異或壞條件時(shí)牛頓法所帶來(lái)的收斂困難,提高了收斂速度。

    初始μ(0)賦值較大,則可向著殘余平方和急劇下降的方向進(jìn)行,在迭代的過(guò)程中,如果F(X(k))遞增則說(shuō)明迭代過(guò)程是發(fā)散的,此時(shí)μ(k)應(yīng)增加,使F(X(k))向著遞減的方向進(jìn)行。隨著迭代次數(shù)的不斷增加,越趨近于最優(yōu)解,μ(k)值應(yīng)逐漸減小,有效避免提前進(jìn)入局部最優(yōu)解。

    對(duì)于每次迭代μ(k)值的選取與實(shí)際減少速度Ared和預(yù)測(cè)減少速度Pred有關(guān)。定義:

    實(shí)際減少速度表達(dá)式為

    (3)

    預(yù)測(cè)減少速度表達(dá)式為

    (4)

    兩者的比值為

    (5)

    當(dāng)比值ρ(k)接近于1時(shí),說(shuō)明線性模型比較準(zhǔn)確,當(dāng)比值ρ(k)接近于0或負(fù)值時(shí),說(shuō)明迭代次數(shù)過(guò)多,以至線性化出現(xiàn)無(wú)效。故:

    (6)

    1.2收斂性證明

    設(shè)潮流方程F(X)=0解收斂于X*,雅克比矩陣的奇異值分解結(jié)果表示如下:

    (7)

    式中:

    迭代過(guò)程中,中間迭代解X(k)所對(duì)應(yīng)的雅克比矩陣奇異值分解為

    (8)

    由J(k)的奇異值分解結(jié)果可知當(dāng)前的迭代值為

    (9)

    (10)

    (11)

    (12)

    根據(jù)式(9)~(12)可得

    (13)

    (14)

    (15)

    因此,X(k)二次收斂于方程組的解,即該算法具有二階收斂性,優(yōu)于只具備一階收斂性的牛頓法。

    2 含分布式電源的改進(jìn)算法的實(shí)現(xiàn)

    2.1節(jié)點(diǎn)編號(hào)方案

    對(duì)于多節(jié)點(diǎn)類型分布式電源并網(wǎng)的配網(wǎng)節(jié)點(diǎn)編號(hào)采用斷點(diǎn)零阻抗連接的編號(hào)方案。在編號(hào)時(shí),首先忽略所有分布式電源,從始端節(jié)點(diǎn)依次向末端節(jié)點(diǎn)編號(hào)1,2,…,N,分布式電源從N+1開(kāi)始編號(hào),看作獨(dú)立的節(jié)點(diǎn),通過(guò)一條具有零阻抗的線路與配電網(wǎng)相連。如圖1所示,原配電網(wǎng)標(biāo)號(hào)為0-11,分布式電源通過(guò)一條零阻抗線路在節(jié)點(diǎn)3出并網(wǎng),編號(hào)為12。

    圖1 簡(jiǎn)單含分布式電源配電網(wǎng)絡(luò)圖Fig.1 Simple distribution network with DG

    該方法保證了原配電網(wǎng)絡(luò)編號(hào)的完整性,分布式電源并網(wǎng)并不影響原網(wǎng)絡(luò)編號(hào),并網(wǎng)點(diǎn)的隨機(jī)選擇性增強(qiáng),對(duì)于程序修改較為便捷。

    2.2改進(jìn)算法實(shí)現(xiàn)步驟

    為確保潮流計(jì)算的可行性,作出以下假設(shè):

    (1)為保證整個(gè)配網(wǎng)是嚴(yán)格吸收型網(wǎng)絡(luò),限制分布式電源并網(wǎng)總?cè)萘?,即?fù)荷總?cè)萘看笥诜植际诫娫纯側(cè)萘俊?/p>

    (2)分布在不同饋線上的負(fù)荷為恒功率負(fù)荷。

    規(guī)定流入節(jié)點(diǎn)功率為正,流出節(jié)點(diǎn)功率為負(fù)。多節(jié)點(diǎn)類型分布式電源并網(wǎng)的配網(wǎng)潮流計(jì)算步驟如下:

    ①節(jié)點(diǎn)編號(hào)

    對(duì)于多節(jié)點(diǎn)類型分布式電源并網(wǎng)的配網(wǎng)節(jié)點(diǎn)編號(hào)按照2.1中闡述的斷點(diǎn)零阻抗連接的編號(hào)方案進(jìn)行。

    ②輸入原始數(shù)據(jù)和信息

    (a)輸入支路阻抗數(shù)據(jù);

    (b)確定分布式電源節(jié)點(diǎn)類型;

    (c)輸入負(fù)荷注入的有功和無(wú)功;(d)輸入平衡節(jié)點(diǎn)的電壓

    根據(jù)引言部分對(duì)常見(jiàn)并網(wǎng)分布式單元的節(jié)點(diǎn)類型的介紹,控制方式與運(yùn)行模式?jīng)Q定了常見(jiàn)分布式電源的節(jié)點(diǎn)類型,主要為PQ節(jié)點(diǎn)、PV節(jié)點(diǎn)、PI節(jié)點(diǎn)以及P-Q(V)節(jié)點(diǎn)4種類型。各類型節(jié)點(diǎn)已知的基本信息以及在潮流計(jì)算前的初值設(shè)定如表1所示。

    表1 分布式電源各類型節(jié)點(diǎn)基本信息及初值設(shè)定

    (3)設(shè)定迭代參數(shù)μ0初值以及計(jì)算精度ε

    (5)若max(ΔP(k),ΔQ(k),ΔV(k)2)≤ε,停止計(jì)算;否則,計(jì)算:

    (16)

    (6)根據(jù)式(3)~(5)計(jì)算Ared(k)、Pred(k)、ρ(k)

    (7)根據(jù)式(6)更新μk,返回第(5)步。

    3 算例驗(yàn)證

    為驗(yàn)證算法的有效性與收斂性,在VC++6.0環(huán)境下進(jìn)行編譯仿真,設(shè)計(jì)如下2組仿真實(shí)驗(yàn)。

    實(shí)驗(yàn)1:在不含分布式電源的配網(wǎng)系統(tǒng)中,分別采用本文改進(jìn)算法與其他傳統(tǒng)潮流計(jì)算對(duì)同一網(wǎng)絡(luò)進(jìn)行潮流計(jì)算,并從收斂速度與收斂精度等方面進(jìn)行比較。

    實(shí)驗(yàn)2:驗(yàn)證算法在含分布式潮流計(jì)算中的實(shí)用性。對(duì)給定節(jié)點(diǎn)安裝分布式電源并采用本文算法進(jìn)行潮流計(jì)算,驗(yàn)證算法在單一節(jié)點(diǎn)并網(wǎng)與多節(jié)點(diǎn)并網(wǎng)中的實(shí)用性,并分析不同節(jié)點(diǎn)類型的分布式電源并網(wǎng)后對(duì)配網(wǎng)的影響。

    3.1系統(tǒng)基本資料

    本文以IEEE33節(jié)點(diǎn)[18]為算法驗(yàn)證系統(tǒng),為便于計(jì)算,采用標(biāo)幺值,設(shè)定基準(zhǔn)值VΒ=12.66kV,SΒ=10 000kVA。系統(tǒng)拓?fù)鋱D如圖2所示。

    圖2 IEEE 33節(jié)點(diǎn)網(wǎng)絡(luò)拓?fù)鋱DFig.2 IEEE 33-node test system

    3.2計(jì)算結(jié)果與分析

    (一)實(shí)驗(yàn)1計(jì)算結(jié)果與分析

    由于牛頓法在潮流計(jì)算中不適用于電阻電抗比值較大的網(wǎng)絡(luò)中,該網(wǎng)絡(luò)所形成的雅克比矩陣不易收斂。故與配網(wǎng)常用的前推回代法計(jì)算結(jié)果相比較,從而驗(yàn)證本文提出算法在收斂性與收斂速度的優(yōu)越性。實(shí)驗(yàn)1計(jì)算結(jié)果如圖2所示。

    表2 各節(jié)點(diǎn)電壓實(shí)際值與計(jì)算值

    兩種算法結(jié)果誤差統(tǒng)計(jì)見(jiàn)表3。通過(guò)與參考值[19]相比較的誤差結(jié)果可知,改進(jìn)算法與實(shí)際值更為接近,誤差率低,更適用于配電網(wǎng)的潮流計(jì)算中。

    表3 兩種算法結(jié)果誤差統(tǒng)計(jì)

    兩種算法計(jì)算時(shí)間見(jiàn)表4。由表中數(shù)據(jù)可知,本文算法大大縮短了計(jì)算時(shí)間,提高了計(jì)算速度。在大規(guī)模配電網(wǎng)潮流計(jì)算時(shí),本文所提出的算法對(duì)于快速潮流分析與實(shí)時(shí)調(diào)度計(jì)算更加適合。

    表4 兩種算法計(jì)算時(shí)間比較

    (二)實(shí)驗(yàn)2計(jì)算結(jié)果與分析

    為驗(yàn)證算法在含分布式潮流計(jì)算中的實(shí)用性,設(shè)定分布式電源并網(wǎng)點(diǎn)與參數(shù)類型如表5所示,由于P-Q(V)型分布式電源根據(jù)給定函數(shù)關(guān)系在求出無(wú)功功率后即可視為PQ型DG,不具備特殊性,故在此不再專門設(shè)置P-Q(V)型節(jié)點(diǎn)。

    表5 分布式電源基本資料

    (1)單一分布式電源并網(wǎng)

    令DG1、DG2、DG3分別在相應(yīng)并網(wǎng)點(diǎn)并入電網(wǎng),通過(guò)改進(jìn)算法進(jìn)行計(jì)算,單一分布式電源并網(wǎng)后迭代次數(shù)與迭代時(shí)間如表6所示。

    表6 單一分布式電源并網(wǎng)潮流計(jì)算迭代次數(shù)與計(jì)算時(shí)間

    由表6數(shù)據(jù)可知,不同類型分布式電源并網(wǎng)后潮流計(jì)算迭代次數(shù)與并網(wǎng)前迭代次數(shù)相差不大,由此說(shuō)明改進(jìn)算法在潮流計(jì)算過(guò)程中的迭代次數(shù)與網(wǎng)絡(luò)規(guī)模相關(guān)性較弱。同時(shí),計(jì)算時(shí)間仍較傳統(tǒng)算法時(shí)間短,說(shuō)明改進(jìn)算法在單一分布式電源并網(wǎng)潮流計(jì)算中能夠提高收斂速度,大大縮短計(jì)算時(shí)間。

    通過(guò)潮流計(jì)算可知不同情況下各節(jié)點(diǎn)電壓分布情況如圖3所示。

    圖3 DG并網(wǎng)后與并網(wǎng)前各節(jié)點(diǎn)電壓Fig.3 The voltage before and after DGs connected

    由圖3中電壓分布情況可知,PI、PQ型DG并網(wǎng)提高了并網(wǎng)點(diǎn)的電壓,同時(shí)也提高了所在支路的電壓,其容量越大提高幅度越明顯,但對(duì)于相鄰線路電壓的改善并不明顯;PV型DG為保證所在節(jié)點(diǎn)的電壓穩(wěn)定,同樣對(duì)相鄰節(jié)點(diǎn)電壓有提升的作用,同時(shí)對(duì)所在支路電壓的提升也有貢獻(xiàn)。

    (2)多種分布式電源并網(wǎng)

    根據(jù)表5中分布式電源基本資料將DG1、 DG2與DG3分別于節(jié)點(diǎn)8、18、31同時(shí)并網(wǎng),通過(guò)改進(jìn)算法所得各節(jié)點(diǎn)電壓分布情況如圖4所示。

    圖4 DG1-3并網(wǎng)后與并網(wǎng)前各節(jié)點(diǎn)電壓Fig.4 The voltage before and after DG1-3 connected

    潮流計(jì)算中迭代次數(shù)為7次,計(jì)算時(shí)間為0.086 s。說(shuō)明隨著分布式電源并網(wǎng)數(shù)目的不斷增多,類型不斷豐富,其迭代次數(shù)與計(jì)算時(shí)間也會(huì)隨之相應(yīng)提高,但仍較傳統(tǒng)算法計(jì)算速度快。由圖可知,隨著配網(wǎng)中多節(jié)點(diǎn)類型的分布式電源廣泛接入,對(duì)各節(jié)點(diǎn)電壓均有較為明顯的改善。

    4 結(jié) 論

    由于DG類型眾多,在潮流計(jì)算中如何正確考慮其節(jié)點(diǎn)類型是問(wèn)題的關(guān)鍵。本文以配電網(wǎng)常見(jiàn)并網(wǎng)分布式電源系統(tǒng)為背景,分析傳統(tǒng)配網(wǎng)潮流計(jì)算方法不足之處,考慮配電網(wǎng)電阻電抗比值較大的實(shí)際問(wèn)題,對(duì)牛頓法進(jìn)行改進(jìn),提出在每次迭代時(shí)對(duì)雅克比矩陣元素進(jìn)行修正,在保證收斂性的前提下提高收斂速度,通過(guò)理論推導(dǎo)與仿真驗(yàn)證得到以下結(jié)論:

    (1)改進(jìn)算法在對(duì)配電網(wǎng)潮流計(jì)算中從準(zhǔn)確性、收斂性以及計(jì)算時(shí)間上均優(yōu)于傳統(tǒng)算法。同時(shí),該方法有效解決了配網(wǎng)中電阻電抗比值較大導(dǎo)致修正方程收斂困難的問(wèn)題。

    (2)本文提出的算法可有效地對(duì)多節(jié)點(diǎn)類型分布式電源并網(wǎng)的配網(wǎng)進(jìn)行潮流計(jì)算。不僅適用于單節(jié)點(diǎn)并網(wǎng),同樣適用于多種分布式電源并網(wǎng)。該算法兼具牛頓法與梯度法的特點(diǎn),如迭代次數(shù)與網(wǎng)絡(luò)規(guī)模無(wú)關(guān),計(jì)算速度快等等。對(duì)于單一分布式電源并網(wǎng)的潮流計(jì)算,其DG類型對(duì)于計(jì)算的迭代次數(shù)與迭代時(shí)間影響不大,與傳統(tǒng)潮流計(jì)算方法相較有大幅提高。然而隨著分布式電源并網(wǎng)數(shù)目的不斷增多,類型不斷豐富,其迭代次數(shù)與計(jì)算時(shí)間也會(huì)隨之相應(yīng)提高,但仍較傳統(tǒng)算法計(jì)算速度快。

    [1] 谷彩連,冷曉華. 分布式電源并網(wǎng)對(duì)配電系統(tǒng)的影響研究[J]. 電器開(kāi)關(guān), 2014, (1):51-58.

    [2] 楊旭英,段建東,楊文宇,等. 含分布式發(fā)電的配電網(wǎng)潮流計(jì)算[J]. 電網(wǎng)技術(shù), 2009, 33(18):139-143.

    [3] 楊雄,衛(wèi)志農(nóng),孫國(guó)強(qiáng),等. 含分布式電源的配電網(wǎng)三相解耦潮流計(jì)算方法[J]. 電力自動(dòng)化設(shè)備, 2014, 34(3):99-107.

    [4] 何昌奇,邵龍,高楊,等. 分布式電源并網(wǎng)后對(duì)電壓分布影響的研究[J]. 中國(guó)電力, 2014, 47(5): 88-91.

    [5] 彭彬,劉寧,吳迪. 配電網(wǎng)潮流計(jì)算中的分布式電源建模[J]. 電力系統(tǒng)及其自動(dòng)化學(xué)報(bào), 2011, 23(2):152-156.

    [6] 巨云濤,吳文傳,張伯明,等. 分布式電源三相穩(wěn)態(tài)模型 [J]. 中國(guó)電機(jī)工程學(xué)報(bào), 2014, 34(10):1509-1518.

    [7] FEIJOO Andres E,CIDRAS Jose. Modeling of wind farms in the load flow analysis [J]. IEEE Trans on Power Systems, 2000, 15(1):110-115.

    [8] 徐青山,劉中澤,陳瀟鵬,等. 含分布式電源的配電網(wǎng)三相時(shí)域仿真[J]. 電力自動(dòng)化設(shè)備, 2015, 35(8):22-30.

    [9] 顧晨,樂(lè)秀璠,張曉明. 基于改進(jìn)前推回代法的弱環(huán)配電網(wǎng)三相潮流計(jì)算[J]. 電 力 系 統(tǒng) 保 護(hù) 與 控 制, 2010, 38(19):160-164.

    [10] 劉耀年,豈小梅,李國(guó)鵬,等. 基于回路阻抗法的配電網(wǎng)潮流計(jì)算[J]. 繼電器, 2004, 32(8):8-10.

    [11] 索南加樂(lè),李懷強(qiáng),羅云照,等. 一種新的配網(wǎng)潮流常Jacobian牛頓算法[J]. 西安交通大學(xué)學(xué)報(bào), 2002, 36(12):1222-1226.

    [12] 孫健,江道灼. 基于牛頓法的配電網(wǎng)絡(luò)Zbus潮流計(jì)算方法[J]. 電網(wǎng)技術(shù), 2004, 28(15):40-44.

    [13] 閆麗梅,謝明霞,徐建軍,等. 含分布式電源的配電網(wǎng)潮流改進(jìn)算法[J]. 電力系統(tǒng)保護(hù)與控制, 2013, 41(5):17-22.

    [14] 劉波,張焰,楊娜.計(jì)及分布式電源的配電網(wǎng)潮流和斷線故障分析算法[J].上海交通大學(xué)學(xué)報(bào), 2008, 42(12):2028-2036.

    [15] 牛煥娜,井天軍,李漢成,等.基于回路分析的含分布式電源配電網(wǎng)簡(jiǎn)化潮流計(jì)算[J]. 電網(wǎng)技術(shù),2013, 37(4):1033-1038.

    [16] 彭彬,劉寧,吳迪. 配電網(wǎng)潮流計(jì)算中的分布式電源建模[J].電力系統(tǒng)及其自動(dòng)化學(xué)報(bào), 2011, 23(2):152-156.

    [17] 李庚銀. 電力系統(tǒng)分析基礎(chǔ)[M]. 北京:機(jī)械工業(yè)出版社, 2011.

    [18] BARAN M E, WU F F.Network reconfiguration in distribution systems for loss reduction and load balancing [J].IEEE Transactions on Power Delivery, 1989, 4(2):1401-1406.

    [19] 王守相,王成山. 現(xiàn)代配電系統(tǒng)分析(第二版)[M]. 北京:高等教育出版社, 2014.

    Improved Power Flow Calculation Method for Distribution Network with DGs

    LIANG Xiao, JIAO Yanjun, JIANG Chenyang

    (School of Electrical and Electronic Engineering, North China Electric Power University, Baoding 071003, China)

    The algorithm is improved by adjusting the Newton Method, which solves the problem of the traditional power flow calculation method for distribution network with multi-node Distributed Generation (DG). The elements in the Jacobian matrix are optimized by Steepest Descend Method and Newton Method at each iteration. The improved method solves the difficulty of the rich types of the DG nodes and the convergence difficulties of the high ratio of resistance and reactance in the Jacobian matrix. The new power flow calculation method uses the zero-impedance connection numbering program. Theoretical analysis and the IEEE33 node stimulation show that the improved algorithm can do a better job in terms of convergence and convergence rate than other traditional methods. According to the results, the new algorithm has a wider range of applications and is more suitable for fast flow analysis of power flow for distribution network with DGs and real-time scheduling computing.

    Distributed Generation; power flow calculation for distribution network; Jacobian matrix; improved algorithm; zero-impedance connection numbering program

    10.3969/j.ISSN.1007-2691.2016.04.10

    2015-11-09.

    梁宵(1991-),女,碩士研究生,研究方向?yàn)殡娏ο到y(tǒng)保護(hù)及其控制;焦彥軍(1963-),男,教授,研究方向?yàn)殡娏ο到y(tǒng)保護(hù)及其控制;蔣晨陽(yáng)(1992-),男,碩士研究生,研究方向?yàn)殡娏ο到y(tǒng)分析、運(yùn)行與控制及電力系統(tǒng)繼電保護(hù)。

    TM744

    A

    1007-2691(2016)04-0059-07

    猜你喜歡
    雅克牛頓潮流
    讀書(shū)的快樂(lè)
    讀書(shū)的快樂(lè)
    曾擔(dān)任過(guò)12年國(guó)際奧委會(huì)主席的雅克·羅格逝世,享年79歲
    牛頓忘食
    風(fēng)中的牛頓
    潮流
    足球周刊(2016年14期)2016-11-02 11:47:59
    潮流
    足球周刊(2016年15期)2016-11-02 11:44:02
    潮流
    足球周刊(2016年10期)2016-10-08 18:50:29
    失信的牛頓
    勇于探索的牛頓
    中文字幕精品免费在线观看视频 | av在线播放精品| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲网站| 中文在线观看免费www的网站| 欧美激情极品国产一区二区三区 | av网站免费在线观看视频| 亚洲精品日韩av片在线观看| 欧美国产精品一级二级三级 | 夜夜看夜夜爽夜夜摸| 亚洲欧洲国产日韩| 91精品国产国语对白视频| 精品人妻熟女av久视频| 久久国产精品男人的天堂亚洲 | 在线精品无人区一区二区三| 国产 一区精品| 日韩强制内射视频| 久久99热6这里只有精品| 一个人免费看片子| 嘟嘟电影网在线观看| av.在线天堂| 国产日韩欧美亚洲二区| 成人午夜精彩视频在线观看| 三级经典国产精品| 日产精品乱码卡一卡2卡三| 亚洲一级一片aⅴ在线观看| 亚洲精品久久午夜乱码| 国产免费一区二区三区四区乱码| 久久久久国产精品人妻一区二区| 国产一级毛片在线| 精品一品国产午夜福利视频| 精品一区二区三卡| 亚洲精品一二三| 国产白丝娇喘喷水9色精品| freevideosex欧美| 国产一区二区在线观看日韩| 久久热精品热| 这个男人来自地球电影免费观看 | 欧美3d第一页| 国产精品伦人一区二区| 一级av片app| 蜜桃久久精品国产亚洲av| 国产精品蜜桃在线观看| 18禁动态无遮挡网站| 热re99久久精品国产66热6| 曰老女人黄片| 狂野欧美白嫩少妇大欣赏| 亚洲精品第二区| 伊人久久国产一区二区| 精品少妇黑人巨大在线播放| 熟女电影av网| 国产精品国产三级国产av玫瑰| 99九九在线精品视频 | 亚洲欧洲日产国产| 麻豆精品久久久久久蜜桃| 一个人看视频在线观看www免费| 老司机影院毛片| 视频区图区小说| 啦啦啦视频在线资源免费观看| 久热久热在线精品观看| 国产成人freesex在线| 街头女战士在线观看网站| 欧美国产精品一级二级三级 | 国产极品天堂在线| 久久久久人妻精品一区果冻| 免费在线观看成人毛片| 尾随美女入室| 人人妻人人澡人人看| 免费观看无遮挡的男女| 中文在线观看免费www的网站| 在线天堂最新版资源| 亚洲综合色惰| 简卡轻食公司| 国产淫语在线视频| 9色porny在线观看| 久久久久精品性色| 午夜免费男女啪啪视频观看| 亚洲精品视频女| 狂野欧美激情性xxxx在线观看| 国产精品99久久99久久久不卡 | 免费观看av网站的网址| 午夜激情久久久久久久| 99热全是精品| a级片在线免费高清观看视频| 国产精品国产三级专区第一集| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品乱码久久久v下载方式| 夫妻性生交免费视频一级片| 亚洲一级一片aⅴ在线观看| 国产精品久久久久成人av| 男人添女人高潮全过程视频| 一区二区三区精品91| 久久人人爽av亚洲精品天堂| 99久久人妻综合| 插逼视频在线观看| √禁漫天堂资源中文www| 久久精品夜色国产| 国产毛片在线视频| 九九久久精品国产亚洲av麻豆| videossex国产| 中文字幕人妻熟人妻熟丝袜美| 精品国产露脸久久av麻豆| 国产男人的电影天堂91| 免费看av在线观看网站| 熟女电影av网| 中国国产av一级| 亚洲精品中文字幕在线视频 | 亚洲av电影在线观看一区二区三区| a级毛色黄片| 十八禁网站网址无遮挡 | 久久久久久久亚洲中文字幕| freevideosex欧美| 天美传媒精品一区二区| 91成人精品电影| 两个人的视频大全免费| 国产综合精华液| 亚洲国产成人一精品久久久| 亚洲av二区三区四区| 日韩大片免费观看网站| 人妻系列 视频| 十八禁网站网址无遮挡 | 精品一品国产午夜福利视频| 美女主播在线视频| 亚洲欧美清纯卡通| 熟女av电影| av免费观看日本| 卡戴珊不雅视频在线播放| 久久久久久久久久久丰满| 你懂的网址亚洲精品在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 哪个播放器可以免费观看大片| 国产成人a∨麻豆精品| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品无人区| 中国国产av一级| 内地一区二区视频在线| 在线观看av片永久免费下载| 卡戴珊不雅视频在线播放| videos熟女内射| 我的老师免费观看完整版| 精品一区在线观看国产| 久热久热在线精品观看| 一级毛片黄色毛片免费观看视频| 9色porny在线观看| 日日啪夜夜爽| 欧美亚洲 丝袜 人妻 在线| 一区二区三区四区激情视频| 精品一区二区三卡| 在线观看免费日韩欧美大片 | 欧美成人午夜免费资源| 亚洲精品自拍成人| 国产高清有码在线观看视频| videos熟女内射| 久久久亚洲精品成人影院| 午夜免费鲁丝| 亚洲第一av免费看| 91aial.com中文字幕在线观看| 久久韩国三级中文字幕| 国产av一区二区精品久久| 日本午夜av视频| 国产极品天堂在线| 国产综合精华液| 高清在线视频一区二区三区| 日产精品乱码卡一卡2卡三| 美女福利国产在线| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| 精品久久久精品久久久| 97超碰精品成人国产| 99热国产这里只有精品6| 高清欧美精品videossex| 国产精品伦人一区二区| 纵有疾风起免费观看全集完整版| 免费观看无遮挡的男女| 久久久欧美国产精品| 精品国产乱码久久久久久小说| 最后的刺客免费高清国语| 80岁老熟妇乱子伦牲交| 一区二区av电影网| 久久精品久久久久久久性| 国产黄片美女视频| 99九九线精品视频在线观看视频| 久久国产精品男人的天堂亚洲 | 日韩一本色道免费dvd| 久久青草综合色| 日本91视频免费播放| 伊人亚洲综合成人网| 嫩草影院新地址| 在线观看一区二区三区激情| av福利片在线| 最近2019中文字幕mv第一页| 性高湖久久久久久久久免费观看| 成人影院久久| 国产一区有黄有色的免费视频| 汤姆久久久久久久影院中文字幕| 一本大道久久a久久精品| 另类亚洲欧美激情| 国产精品国产三级国产av玫瑰| 麻豆精品久久久久久蜜桃| 成人亚洲欧美一区二区av| 国产精品三级大全| 高清毛片免费看| 天天躁夜夜躁狠狠久久av| 精品少妇久久久久久888优播| 老女人水多毛片| 欧美精品人与动牲交sv欧美| av国产久精品久网站免费入址| 久久久精品免费免费高清| 国产真实伦视频高清在线观看| 人妻人人澡人人爽人人| 十八禁网站网址无遮挡 | 欧美激情国产日韩精品一区| 久久人人爽人人爽人人片va| 久久精品国产亚洲av天美| 国内少妇人妻偷人精品xxx网站| 久久精品久久精品一区二区三区| 久热这里只有精品99| 精华霜和精华液先用哪个| 在现免费观看毛片| 黄色日韩在线| 日本av手机在线免费观看| 国产色爽女视频免费观看| 免费人妻精品一区二区三区视频| 老司机亚洲免费影院| 国产欧美亚洲国产| 成人18禁高潮啪啪吃奶动态图 | 免费看不卡的av| 久久久久久伊人网av| 韩国高清视频一区二区三区| 男的添女的下面高潮视频| 在线观看免费高清a一片| 国产午夜精品一二区理论片| 国产精品人妻久久久影院| 国产精品伦人一区二区| 亚洲精品久久久久久婷婷小说| 久久午夜福利片| 国产精品不卡视频一区二区| 高清av免费在线| 三上悠亚av全集在线观看 | 一级av片app| 在线看a的网站| 高清不卡的av网站| 久久毛片免费看一区二区三区| 久久久久久久久大av| 妹子高潮喷水视频| 亚洲精品自拍成人| 简卡轻食公司| 99视频精品全部免费 在线| 国产日韩一区二区三区精品不卡 | 日韩在线高清观看一区二区三区| 国产成人精品婷婷| 九九爱精品视频在线观看| 婷婷色麻豆天堂久久| 春色校园在线视频观看| 伦理电影大哥的女人| 成人影院久久| 高清在线视频一区二区三区| 人人妻人人澡人人看| 麻豆精品久久久久久蜜桃| 寂寞人妻少妇视频99o| 日韩成人av中文字幕在线观看| 国产色婷婷99| 久久久久久伊人网av| 嘟嘟电影网在线观看| 22中文网久久字幕| av不卡在线播放| 看免费成人av毛片| 一级二级三级毛片免费看| 久久久久久久亚洲中文字幕| 一级毛片 在线播放| 欧美日韩亚洲高清精品| 如日韩欧美国产精品一区二区三区 | 国产 一区精品| 老司机影院成人| 亚洲av二区三区四区| 国产综合精华液| 午夜免费观看性视频| 国产成人91sexporn| 国产精品伦人一区二区| 街头女战士在线观看网站| 精品人妻熟女av久视频| 边亲边吃奶的免费视频| 高清欧美精品videossex| 国产精品.久久久| 久久午夜综合久久蜜桃| 午夜视频国产福利| 国产一区亚洲一区在线观看| 极品少妇高潮喷水抽搐| www.色视频.com| 如何舔出高潮| 欧美日韩在线观看h| 下体分泌物呈黄色| 大话2 男鬼变身卡| 国产亚洲精品久久久com| 三级经典国产精品| 国精品久久久久久国模美| 亚洲欧美成人精品一区二区| 黑人猛操日本美女一级片| 在线观看一区二区三区激情| 黑人高潮一二区| 蜜桃久久精品国产亚洲av| 国产精品无大码| 国产成人aa在线观看| 制服丝袜香蕉在线| 亚洲欧美清纯卡通| 香蕉精品网在线| 男的添女的下面高潮视频| 亚洲精品国产av成人精品| 不卡视频在线观看欧美| 天堂8中文在线网| 午夜免费男女啪啪视频观看| 国产视频首页在线观看| 成人午夜精彩视频在线观看| 黑人猛操日本美女一级片| 精品一区二区三卡| 一边亲一边摸免费视频| 日韩精品免费视频一区二区三区 | av.在线天堂| 毛片一级片免费看久久久久| 伦理电影大哥的女人| 久久99热这里只频精品6学生| 国产精品国产av在线观看| 国产一区亚洲一区在线观看| 日本欧美视频一区| 在线观看三级黄色| 丝袜喷水一区| 亚洲美女视频黄频| 亚洲美女搞黄在线观看| av福利片在线| 九九爱精品视频在线观看| 69精品国产乱码久久久| 中文字幕亚洲精品专区| 久久久久久久久大av| 晚上一个人看的免费电影| 久久99精品国语久久久| 国产精品国产三级国产专区5o| 天天躁夜夜躁狠狠久久av| 国产淫片久久久久久久久| 91久久精品国产一区二区三区| 18禁裸乳无遮挡动漫免费视频| 麻豆成人午夜福利视频| 男的添女的下面高潮视频| 成年av动漫网址| 午夜精品国产一区二区电影| 日本黄色片子视频| 我要看黄色一级片免费的| 日韩免费高清中文字幕av| 99久国产av精品国产电影| 在线观看www视频免费| a级毛色黄片| 九九爱精品视频在线观看| 久久久久久久久久久丰满| a级片在线免费高清观看视频| 亚洲av二区三区四区| 一级,二级,三级黄色视频| 久久精品久久久久久噜噜老黄| 少妇的逼好多水| 卡戴珊不雅视频在线播放| 国产片特级美女逼逼视频| 欧美日韩亚洲高清精品| 秋霞伦理黄片| 人人妻人人爽人人添夜夜欢视频 | 国产精品偷伦视频观看了| 日韩大片免费观看网站| 国产高清国产精品国产三级| 国语对白做爰xxxⅹ性视频网站| 久久人人爽人人片av| 在线观看人妻少妇| 在线精品无人区一区二区三| 亚洲国产精品一区二区三区在线| 观看av在线不卡| 免费久久久久久久精品成人欧美视频 | 91成人精品电影| 内射极品少妇av片p| 中文天堂在线官网| 丰满人妻一区二区三区视频av| 91久久精品国产一区二区三区| 高清欧美精品videossex| 中文字幕免费在线视频6| 国产av一区二区精品久久| 熟女电影av网| 看非洲黑人一级黄片| 中国国产av一级| 亚洲欧美中文字幕日韩二区| 99久久人妻综合| 日本黄色日本黄色录像| a级毛片在线看网站| 啦啦啦在线观看免费高清www| 男男h啪啪无遮挡| 亚洲人与动物交配视频| 亚洲国产毛片av蜜桃av| 看免费成人av毛片| 免费人妻精品一区二区三区视频| 又粗又硬又长又爽又黄的视频| 亚洲高清免费不卡视频| 蜜桃久久精品国产亚洲av| 女性生殖器流出的白浆| 最新中文字幕久久久久| 久久国产亚洲av麻豆专区| 99热这里只有精品一区| 国产精品一区二区三区四区免费观看| 婷婷色综合大香蕉| 亚洲,一卡二卡三卡| 国产91av在线免费观看| 成人免费观看视频高清| 一区二区av电影网| 午夜91福利影院| 丁香六月天网| 免费观看在线日韩| 亚洲国产欧美在线一区| av女优亚洲男人天堂| 两个人免费观看高清视频 | 男人爽女人下面视频在线观看| 日韩视频在线欧美| 男女啪啪激烈高潮av片| 日韩强制内射视频| 久久99热6这里只有精品| 日韩,欧美,国产一区二区三区| 久久精品国产自在天天线| 精品少妇内射三级| videos熟女内射| 91久久精品电影网| 免费av中文字幕在线| 国产 精品1| 尾随美女入室| 欧美日韩精品成人综合77777| 国产成人免费无遮挡视频| 丝瓜视频免费看黄片| 精品国产一区二区三区久久久樱花| 国产欧美另类精品又又久久亚洲欧美| 在现免费观看毛片| 亚洲精华国产精华液的使用体验| 在线观看免费日韩欧美大片 | 国产精品女同一区二区软件| 国产精品伦人一区二区| 看十八女毛片水多多多| 一本—道久久a久久精品蜜桃钙片| 老司机影院成人| 麻豆成人av视频| 亚洲伊人久久精品综合| 简卡轻食公司| 亚洲欧美日韩东京热| 菩萨蛮人人尽说江南好唐韦庄| 一本一本综合久久| kizo精华| 免费播放大片免费观看视频在线观看| 精品久久久精品久久久| 久久婷婷青草| 欧美日韩视频精品一区| 日本黄色片子视频| 99国产精品免费福利视频| 国产伦精品一区二区三区视频9| 中文乱码字字幕精品一区二区三区| 在线观看三级黄色| 2021少妇久久久久久久久久久| 亚洲人成网站在线观看播放| 婷婷色综合大香蕉| 美女cb高潮喷水在线观看| 人妻 亚洲 视频| 亚洲熟女精品中文字幕| 亚洲综合精品二区| 午夜视频国产福利| 亚洲天堂av无毛| 一级毛片我不卡| 国产又色又爽无遮挡免| 最黄视频免费看| 色婷婷久久久亚洲欧美| 麻豆成人av视频| 我的女老师完整版在线观看| 国产精品一区www在线观看| 欧美日本中文国产一区发布| 97超视频在线观看视频| 免费黄色在线免费观看| 九草在线视频观看| 各种免费的搞黄视频| 免费在线观看成人毛片| 波野结衣二区三区在线| 久久久精品94久久精品| 丝瓜视频免费看黄片| 91午夜精品亚洲一区二区三区| 国产黄片视频在线免费观看| 在线观看一区二区三区激情| 久久99精品国语久久久| 精品久久久精品久久久| 尾随美女入室| 亚州av有码| 香蕉精品网在线| 亚洲欧美一区二区三区黑人 | 高清不卡的av网站| 成人综合一区亚洲| 一区二区三区精品91| 国模一区二区三区四区视频| 日韩电影二区| 欧美日韩视频精品一区| 在现免费观看毛片| www.av在线官网国产| 一级片'在线观看视频| 国产一级毛片在线| 日本黄大片高清| 免费av中文字幕在线| a 毛片基地| 国产欧美亚洲国产| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人手机| 亚洲人成网站在线观看播放| 少妇人妻久久综合中文| 啦啦啦啦在线视频资源| 免费观看av网站的网址| 91在线精品国自产拍蜜月| 99久久精品一区二区三区| 视频区图区小说| 国产片特级美女逼逼视频| 亚洲成色77777| 亚洲av电影在线观看一区二区三区| av在线app专区| 欧美日韩在线观看h| 2021少妇久久久久久久久久久| 免费大片黄手机在线观看| √禁漫天堂资源中文www| 精品亚洲成a人片在线观看| 国产精品无大码| 亚洲国产色片| 在线天堂最新版资源| 国产黄片视频在线免费观看| 80岁老熟妇乱子伦牲交| 日本av免费视频播放| 欧美xxxx性猛交bbbb| 插逼视频在线观看| 免费大片18禁| 赤兔流量卡办理| 国产男人的电影天堂91| 免费黄网站久久成人精品| 街头女战士在线观看网站| av不卡在线播放| 国产精品.久久久| 午夜免费鲁丝| 97超碰精品成人国产| 亚洲精品色激情综合| 女的被弄到高潮叫床怎么办| av.在线天堂| 午夜福利视频精品| 高清视频免费观看一区二区| 九色成人免费人妻av| 又爽又黄a免费视频| 久久国产乱子免费精品| 麻豆乱淫一区二区| 久久久久久久久久久久大奶| 国语对白做爰xxxⅹ性视频网站| 国产成人精品福利久久| 亚洲av电影在线观看一区二区三区| 在线观看免费视频网站a站| 国产精品人妻久久久影院| 亚洲图色成人| 美女主播在线视频| 99热6这里只有精品| 美女大奶头黄色视频| 在线 av 中文字幕| 99视频精品全部免费 在线| 亚洲成人手机| 看十八女毛片水多多多| 免费观看性生交大片5| 欧美高清成人免费视频www| 国产亚洲91精品色在线| 又黄又爽又刺激的免费视频.| 亚洲av成人精品一二三区| 美女视频免费永久观看网站| 日产精品乱码卡一卡2卡三| 97精品久久久久久久久久精品| 亚洲av日韩在线播放| 99久久精品国产国产毛片| 国产熟女欧美一区二区| 美女福利国产在线| 色5月婷婷丁香| 99九九线精品视频在线观看视频| a级一级毛片免费在线观看| 在线观看三级黄色| 国产免费又黄又爽又色| av又黄又爽大尺度在线免费看| 精品人妻一区二区三区麻豆| av免费观看日本| 三上悠亚av全集在线观看 | 人妻系列 视频| 最新的欧美精品一区二区| 麻豆乱淫一区二区| 免费人妻精品一区二区三区视频| 99热全是精品| 日日啪夜夜撸| 丰满饥渴人妻一区二区三| 亚洲欧美精品专区久久| 免费看av在线观看网站| 一区在线观看完整版| 晚上一个人看的免费电影| 亚洲精品日韩av片在线观看| 久热久热在线精品观看| 欧美日韩视频精品一区| 晚上一个人看的免费电影| 国产一区二区在线观看日韩| 国产精品久久久久久久久免| 亚洲熟女精品中文字幕| 欧美日韩亚洲高清精品| 最新中文字幕久久久久| 麻豆精品久久久久久蜜桃| 久久99蜜桃精品久久| 色吧在线观看| 在线观看国产h片| 日本av免费视频播放| 99热网站在线观看| 亚洲av福利一区| 少妇熟女欧美另类| 精品国产一区二区三区久久久樱花| 97在线视频观看| 亚洲国产成人一精品久久久| 少妇人妻一区二区三区视频| 精品久久久精品久久久| 少妇被粗大的猛进出69影院 |