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

    基于波形統(tǒng)計(jì)分析的頻率振蕩相關(guān)聯(lián)機(jī)組快速識(shí)別技術(shù)

    2021-03-06 09:12:38朱余啟陳汝昌戈本星沈力
    廣東電力 2021年2期
    關(guān)鍵詞:皮爾森標(biāo)準(zhǔn)差出力

    朱余啟, 陳汝昌,戈本星,沈力

    (1. 云南電網(wǎng)有限責(zé)任公司云南電力調(diào)度控制中心, 云南 昆明 650041;2. 南京南瑞繼保電氣有限公司,江蘇 南京 211102)

    云南電網(wǎng)與南方主網(wǎng)異步運(yùn)行之后[1],由于云南電網(wǎng)水電裝機(jī)容量大且占比高[2]、省內(nèi)負(fù)荷水平低、直流外送比例高的運(yùn)行特征,網(wǎng)內(nèi)先后出現(xiàn)數(shù)次頻率振蕩情況[3]。當(dāng)存在一次調(diào)頻設(shè)置不合理[4-5]、二次調(diào)頻參數(shù)設(shè)置不合理[6]、切機(jī)切負(fù)荷產(chǎn)生錯(cuò)誤動(dòng)作[7]、直流頻率限制器動(dòng)作[8]等情況時(shí),均有可能發(fā)生頻率振蕩。

    電力系統(tǒng)振蕩的診斷與緊急控制的核心問題是振蕩源識(shí)別,目前有2種主流方法:一是基于離線數(shù)據(jù)進(jìn)行仿真與分析;二是基于相量測量裝置(phasor measurement unit,PMU)毫秒級(jí)實(shí)測數(shù)據(jù)進(jìn)行實(shí)時(shí)數(shù)據(jù)處理與分析。前者對模型參數(shù)的準(zhǔn)確性要求較高,且對設(shè)備異常、模型與參數(shù)變化所引起的振蕩往往無能為力[9]。后者需要基于廣域測量系統(tǒng)實(shí)時(shí)數(shù)據(jù)進(jìn)行大量的處理與分析,通過分析各種電氣量(包括電壓、電流、出力、功角等)來診斷振蕩。本文只討論后者實(shí)時(shí)分析方法,并基于此提出解決方案,即采用秒級(jí)數(shù)據(jù)進(jìn)行波形統(tǒng)計(jì)分析的方法。文獻(xiàn)[10]使用振蕩特征迭代函數(shù)的方法證明了秒級(jí)曲線也可以用作識(shí)別電網(wǎng)頻率振蕩。

    實(shí)時(shí)分析方法所采取的方案有電氣量的暫態(tài)分量判別法和穩(wěn)態(tài)分量判別法,大部分文獻(xiàn)采用的是暫態(tài)分量判別法。文獻(xiàn)[11]利用電壓行波暫態(tài)量分析定位振蕩源頭;在電力系統(tǒng)發(fā)生機(jī)電振蕩時(shí),發(fā)電機(jī)轉(zhuǎn)子的動(dòng)能、勢能與網(wǎng)絡(luò)中的勢能不斷轉(zhuǎn)化,文獻(xiàn)[12]利用發(fā)電機(jī)轉(zhuǎn)子的動(dòng)能、勢能與網(wǎng)絡(luò)中勢能的暫態(tài)量分析定位振蕩源頭;文獻(xiàn)[13]利用振蕩時(shí)各個(gè)機(jī)組的轉(zhuǎn)子相位、電壓相位的不同分布情況判別振蕩源頭;文獻(xiàn)[14]根據(jù)失步中心同側(cè)和異側(cè)的站點(diǎn)電壓頻率傳播方向不同,在大區(qū)域互聯(lián)電網(wǎng)判斷失步振蕩中心。利用暫態(tài)分量進(jìn)行頻率振蕩判別的方法,優(yōu)點(diǎn)是物理模型準(zhǔn)確,能夠從理論的角度確定分析振蕩源頭,缺點(diǎn)是對PMU數(shù)據(jù)的可靠性依賴較高,算法復(fù)雜,難以工程化。

    利用電氣量的穩(wěn)態(tài)分量判別系統(tǒng)振蕩的應(yīng)用情況較少。文獻(xiàn)[15]提出了使用發(fā)電機(jī)組功率譜密度來判別振蕩的方案,但該方案局限于定位發(fā)電機(jī)機(jī)械功率擾動(dòng)源。使用穩(wěn)態(tài)分量判別法的優(yōu)點(diǎn)是算法簡單,易工程化應(yīng)用,缺點(diǎn)是物理意義不清晰,目前僅在高壓線路的振蕩源頭定位有明確的物理含義及原理推導(dǎo),對于其他頻率振蕩情況并不能從理論上分析定位振蕩源頭。

    本文并不需要嚴(yán)格對振蕩源頭進(jìn)行定位,而只需要判斷出頻率振蕩并識(shí)別影響頻率振蕩的關(guān)聯(lián)機(jī)組,因此采用穩(wěn)態(tài)分量判別方法。同時(shí),通過分析波形數(shù)獲取波形特征進(jìn)行振蕩機(jī)組判別的研究較少。文獻(xiàn)[16]通過波形相關(guān)性分析了汽輪機(jī)閥門開度與機(jī)組功率振蕩的關(guān)系;文獻(xiàn)[17]通過分析電壓波形與衰減正弦波的相關(guān)性判斷是否發(fā)生低頻振蕩;文獻(xiàn)[18]通過對低頻含噪等幅振蕩信號(hào)波形進(jìn)行自相關(guān)功率譜分析來判別是否振蕩。所以,本文提出通過分析頻率和有功功率的穩(wěn)態(tài)振蕩波形來檢測頻率振蕩,并判別相關(guān)聯(lián)的振蕩機(jī)組。該方法計(jì)算量小,結(jié)果直觀,易于工程化實(shí)現(xiàn)。

    文獻(xiàn)[19]通過云南異步聯(lián)網(wǎng)擾動(dòng)實(shí)驗(yàn)實(shí)際數(shù)據(jù)說明,水電機(jī)組由于負(fù)阻尼效應(yīng)而助增了頻率振蕩,火電機(jī)組由于正阻尼抑制了頻率振蕩,并且通過波形可以看出助增振蕩的有功功率穩(wěn)態(tài)波形基本與頻率穩(wěn)態(tài)波形同相位,而抑制頻率振蕩的有功功率穩(wěn)態(tài)波形基本與頻率反相位;同時(shí)通過BPA對水錘效應(yīng)的水電機(jī)組調(diào)速器進(jìn)行建模,成功復(fù)現(xiàn)了水電機(jī)組因?yàn)樨?fù)阻尼助增頻率擾動(dòng)的現(xiàn)象,仿真出的波形與實(shí)際實(shí)驗(yàn)的波形近似。據(jù)此原理,本文的頻率振蕩判別與振蕩機(jī)組快速識(shí)別技術(shù)的分析步驟如下:首先使用Prony算法對頻率波形進(jìn)行分析,判斷是否發(fā)生振蕩[20-21],對出力波形計(jì)算標(biāo)準(zhǔn)差,當(dāng)標(biāo)準(zhǔn)差滿足一定大小時(shí)才認(rèn)為其影響頻率振蕩;然后將頻率和出力2個(gè)不同量綱的數(shù)據(jù)標(biāo)準(zhǔn)化為均值為0、方差為1的數(shù)據(jù);計(jì)算頻率與出力的皮爾森系數(shù),當(dāng)皮爾森系數(shù)絕對值滿足一定大小時(shí)才認(rèn)為其與頻率振蕩相關(guān)聯(lián),并且通過皮爾森系數(shù)的正負(fù)號(hào)確定該出力是助增還是抑制頻率振蕩;最后以標(biāo)準(zhǔn)差和皮爾森系數(shù)之積作為綜合相關(guān)系數(shù),對關(guān)聯(lián)振蕩機(jī)組助增或抑制頻率振蕩的能力進(jìn)行排序。

    1 頻率振蕩辨識(shí)

    Prony算法是電力系統(tǒng)中普遍使用的辨識(shí)電網(wǎng)低頻振蕩的方法,本文采用Prony算法對頻率振蕩進(jìn)行識(shí)別。同時(shí),使用調(diào)度主站系統(tǒng)采集的電網(wǎng)頻率秒級(jí)數(shù)據(jù),以便實(shí)現(xiàn)快速采集、快速計(jì)算、快速分析。

    1.1 Prony算法原理

    Prony算法用一組衰減正弦函數(shù)的線性組合來擬和等間距采樣數(shù)據(jù),從中分析出信號(hào)的幅值、頻率、相位、阻尼因子等信息。

    y=∑Ae-λtcos (2πft+θ).

    (1)

    式中:y為待分析的信號(hào);t為時(shí)間;A為信號(hào)幅值;λ為衰減因子;f為頻率;θ為相位。

    當(dāng)輸入一串等距離散序列x(0),x(1),…,x(N-1),N為離散序列的點(diǎn)數(shù),需要用P階的衰減正弦函數(shù)來表達(dá)這一串序列,即

    (2)

    式中:bi=Aejθi,zi=eλi+j2πfi。即已知x(0),x(1),…,x(N-1),求解bi與zi。

    Prony算法的求解流程如下:

    a)構(gòu)造二階樣本矩陣。定義二階樣本矩陣函數(shù)

    i,j=0,1,…,P.

    (3)

    (4)

    為了使樣本矩陣全面反映系統(tǒng)信息,維數(shù)L取其最大值N/2。

    b)確定計(jì)算階數(shù)。在計(jì)算前需要通過一定的規(guī)則確定模型階數(shù)P,一般采用奇異值分解(singular value decomposition,SVD)算法。SVD算法的實(shí)質(zhì)是求解一個(gè)最優(yōu)化問題,利用Frobineous范數(shù)定階方法確定樣本矩陣的有效秩序P。

    c)求解回歸參數(shù)。求解式(5)的回歸參數(shù)a1,a2,…,aP,其中εP為自定義的最小誤差能量:

    (5)

    d)求解特征方程。特征方程的形式為1+a1z-1+…+aPz-P=0,求解得到特征根zi(i=1,2,…,P),即前文的zi=e(λi+j2πfi),也稱為Prony極點(diǎn)。

    (6)

    得到整個(gè)采樣區(qū)間的表達(dá)式

    最終利用bi=Aejθi,zi=eλi+j2πfi,求解θi、fi。

    1.2 基于Prony算法的頻率振蕩辨識(shí)方法

    綜合異步聯(lián)網(wǎng)的云南電網(wǎng)多年運(yùn)行情況,給定的頻率振蕩定值設(shè)置如下:頻率振蕩幅值定值為0.1 Hz,當(dāng)振蕩分量幅值大于0.1 Hz時(shí)認(rèn)為發(fā)生了振蕩;衰減因子絕對值定值為0.1,衰減因子絕對值過大說明振蕩衰減迅速,不納入考慮;當(dāng)振蕩分量的幅值大于0.1 Hz并且衰減因子絕對值小于0.1時(shí),判定該振蕩分量有效。

    2 關(guān)聯(lián)振蕩機(jī)組識(shí)別及排序方法

    2.1 計(jì)算出力標(biāo)準(zhǔn)差

    在進(jìn)行數(shù)據(jù)標(biāo)準(zhǔn)化之前,首先計(jì)算出力曲線的標(biāo)準(zhǔn)差,以判斷出力波動(dòng)的程度。當(dāng)標(biāo)準(zhǔn)差過小時(shí)認(rèn)為對振蕩沒有影響,當(dāng)標(biāo)準(zhǔn)差滿足一定條件時(shí)才認(rèn)為能夠影響頻率振蕩。后續(xù)將采用標(biāo)準(zhǔn)差與皮爾森系數(shù)之積作為綜合相關(guān)系數(shù),作為排序的標(biāo)準(zhǔn)。

    2.2 數(shù)據(jù)標(biāo)準(zhǔn)化

    進(jìn)行相關(guān)性分析之前,需要對電網(wǎng)頻率數(shù)據(jù)與機(jī)組出力數(shù)據(jù)進(jìn)行預(yù)處理,將不同量綱的頻率和出力數(shù)據(jù)標(biāo)準(zhǔn)化為均值為0、方差為1的數(shù)據(jù)。

    (7)

    式中:Xi為第i個(gè)原始數(shù)據(jù),z(Xi)為其標(biāo)準(zhǔn)化后的值;μ為均值;σ為標(biāo)準(zhǔn)差。

    2.3 計(jì)算頻率與出力皮爾森系數(shù)

    2.3.1 皮爾森系數(shù)

    皮爾森相關(guān)系數(shù)

    (8)

    式中:cov(X,Y)為向量X與Y的協(xié)方差;σX為向量X的標(biāo)準(zhǔn)差;σY為向量Y的標(biāo)準(zhǔn)差。皮爾森相關(guān)系數(shù)(Pearson correlation coefficient)也稱皮爾森積矩相關(guān)系數(shù)(Pearson product-moment correlation coefficient),是一種線性相關(guān)系數(shù)。皮爾森相關(guān)系數(shù)是用來反映2個(gè)變量線性相關(guān)程度的統(tǒng)計(jì)量。rc描述的是2個(gè)變量間線性相關(guān)強(qiáng)弱的程度,rc的絕對值越大表明相關(guān)性越強(qiáng)。

    簡單的相關(guān)系數(shù)的分類:rc值為0.8~1.0表示極強(qiáng)相關(guān);rc值為0.6~0.8表示強(qiáng)相關(guān);rc值為0.4~0.6表示中等程度相關(guān);rc值為0.2~0.4表示弱相關(guān);rc值為0.0~0.2表示極弱相關(guān)或無相關(guān)。

    rc的取值在-1與+1之間:若rc>0,表明2個(gè)變量是正相關(guān),即一個(gè)變量的值越大,另一個(gè)變量的值也會(huì)越大;若rc<0,表明2個(gè)變量是負(fù)相關(guān),即一個(gè)變量的值越大,另一個(gè)變量的值反而會(huì)越小。

    2.3.2 仿真計(jì)算結(jié)果

    雖然皮爾森系數(shù)rc描述的是2個(gè)變量間線性相關(guān)強(qiáng)弱的程度,但當(dāng)變量為正弦波時(shí),亦可用皮爾森系數(shù)描述相關(guān)性。仿真結(jié)果如圖1、圖2、圖3所示,其中幅值為歸一化數(shù)值。

    圖1 幅值不同的正弦波相關(guān)系數(shù)Fig.1 Correlation coefficients of two sine waves with different amplitudes

    圖2 頻率不同的正弦波相關(guān)系數(shù)Fig.2 Correlation coefficients of two sine waves with different frequency

    圖3 相位不同的正弦波相關(guān)系數(shù)Fig.3 Correlation coefficients of two sine waves with different phases

    圖2中,2個(gè)正弦波頻率分別為1/2π和1.5/2π時(shí),皮爾森系數(shù)計(jì)算結(jié)果為-0.243 5,即弱負(fù)相關(guān)。

    圖3中:2個(gè)正弦波相位差為π/3時(shí),皮爾森系數(shù)計(jì)算結(jié)果為0.505 4,即正相關(guān);2個(gè)正弦波相位差為2π/3時(shí),皮爾森系數(shù)計(jì)算結(jié)果為-0.433 9,即負(fù)相關(guān)。

    通過比較頻率波形與機(jī)組出力波形的相關(guān)系數(shù)可以判斷出頻率波形與機(jī)組出力波形的相位差,從而判斷出機(jī)組出力是助增頻率振蕩還是抑制頻率振蕩,以及助增頻率振蕩或是抑制頻率振蕩關(guān)聯(lián)性的大小,從而對關(guān)聯(lián)機(jī)組進(jìn)行排序,提供給調(diào)度員分析決策。

    2.4 關(guān)聯(lián)振蕩機(jī)組排序

    本文采用皮爾森系數(shù)與標(biāo)準(zhǔn)差之積作為綜合相關(guān)系數(shù),以綜合相關(guān)系數(shù)的大小對機(jī)組出力進(jìn)行排序。根據(jù)運(yùn)行經(jīng)驗(yàn),需要對皮爾森系數(shù)和標(biāo)準(zhǔn)差設(shè)立計(jì)算門檻。皮爾森系數(shù)絕對值必須大于0.4,即出力與頻率的相關(guān)性必須滿足一定要求;標(biāo)準(zhǔn)差必須大于2 MW,即認(rèn)為出力波動(dòng)到一定情況才認(rèn)為對電網(wǎng)有影響。滿足這兩者定值才參與排序,否則認(rèn)為和頻率振蕩無關(guān)。

    當(dāng)皮爾森系數(shù)和標(biāo)準(zhǔn)差滿足門檻時(shí):綜合相關(guān)系數(shù)大于0,認(rèn)為助增頻率振蕩;綜合相關(guān)系數(shù)小于0,認(rèn)為抑制頻率振蕩;綜合相關(guān)系數(shù)絕對值越大,影響頻率振蕩越強(qiáng);綜合相關(guān)系數(shù)絕對值越小,助增頻率振蕩越弱。

    3 實(shí)際算例分析

    3.1 頻率振蕩辨識(shí)結(jié)果

    對2020年7月某日發(fā)生的頻率振蕩數(shù)據(jù)進(jìn)行分析,選取5 min(即樣本點(diǎn)數(shù)為300)的振蕩數(shù)據(jù),頻率數(shù)據(jù)如圖4所示。對實(shí)際頻率進(jìn)行Prony分析,結(jié)果見表1。

    表1 頻率Prony分析Tab.1 Prony analysis of frequency

    根據(jù)振蕩幅值定值0.1以及衰減因子絕對值定值0.1,即振蕩分量的幅值大于0.1并且衰減因子絕對值小于0.1時(shí)判定該振蕩分量有效。這里認(rèn)為模式1的振蕩分量有效,即認(rèn)為發(fā)生了頻率振蕩。

    3.2 關(guān)聯(lián)振蕩機(jī)組識(shí)別結(jié)果

    分別以糯扎渡電廠、小灣電廠、普西橋電廠、滇東電廠為例進(jìn)行分析。

    糯扎渡電廠出力與頻率曲線如圖5所示,計(jì)算結(jié)果為:皮爾森系數(shù)0.69、標(biāo)準(zhǔn)差15.43、綜合相關(guān)系數(shù)10.65。從皮爾森系數(shù)可以看出糯扎渡電廠出力助增了頻率振蕩,從出力標(biāo)準(zhǔn)差看出電廠發(fā)生了一定程度的振蕩,所以認(rèn)為綜合相關(guān)系數(shù)有效,該電廠為關(guān)聯(lián)振蕩的電廠,并且助增頻率振蕩。

    小灣電廠出力與頻率曲線如圖6所示,計(jì)算結(jié)果為:皮爾森系數(shù)0.12、標(biāo)準(zhǔn)差49.94、綜合相關(guān)系數(shù)6.17。從皮爾森系數(shù)可以看出小灣電廠出力與頻率振蕩無關(guān),從出力標(biāo)準(zhǔn)差看出電廠出力變化較大,所以認(rèn)為綜合相關(guān)系數(shù)無效,該電廠與振蕩無關(guān)。

    普西橋電廠出力與頻率曲線如圖7所示,計(jì)算結(jié)果為:皮爾森系數(shù)0.50、標(biāo)準(zhǔn)差4.00、綜合相關(guān)系數(shù)2.00。從皮爾森系數(shù)可以看出普西橋電廠出力助增了頻率振蕩,從出力標(biāo)準(zhǔn)差看出電廠發(fā)生了一定程度的振蕩,所以認(rèn)為綜合相關(guān)系數(shù)有效,該電廠為關(guān)聯(lián)振蕩的電廠,并且助增頻率振蕩。

    圖4 頻率曲線Fig.4 Frequency curve

    圖5 糯扎渡電廠出力與頻率曲線Fig.5 Output and frequency curves of Nuozhadu power plant

    圖6 小灣電廠出力與頻率曲線Fig.6 Output and frequency curves of Xiaowan power plant

    圖7 普西橋電廠出力與頻率曲線Fig.7 Output and frequency curves of Puxiqiao power plant

    滇東電廠出力與頻率曲線如圖8所示。計(jì)算結(jié)果如下:皮爾森系數(shù)-0.85、標(biāo)準(zhǔn)差3.85、綜合相關(guān)系數(shù)-3.27。從皮爾森系數(shù)可以看出滇東電廠出力抑制頻率振蕩,從出力標(biāo)準(zhǔn)差看出電廠出力變化較大,所以認(rèn)為綜合相關(guān)系數(shù)有效,該電廠為關(guān)聯(lián)振蕩的電廠,并且抑制頻率振蕩。

    圖8 滇東電廠出力與頻率曲線Fig.8 Output and frequency curves of Diandong power plant

    4 結(jié)束語

    云南電網(wǎng)與南方電網(wǎng)主網(wǎng)異步運(yùn)行之后,云南電網(wǎng)頻率穩(wěn)定問題逐漸突出。當(dāng)電網(wǎng)頻率發(fā)生振蕩時(shí),需要監(jiān)控系統(tǒng)能快速響應(yīng),在判別頻率振蕩的同時(shí)快速定位振蕩源頭。本文提出:首先,使用Prony算法對頻率秒級(jí)曲線進(jìn)行分析,判斷頻率是否發(fā)生振蕩;其次,對出力波形計(jì)算標(biāo)準(zhǔn)差,當(dāng)標(biāo)準(zhǔn)差滿足一定大小才認(rèn)為其影響頻率振蕩;然后,對出力數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,將頻率和出力2個(gè)不同量綱的數(shù)據(jù)標(biāo)準(zhǔn)化為均值為0、方差為1的數(shù)據(jù);接著,計(jì)算頻率與出力的皮爾森系數(shù),當(dāng)皮爾森系數(shù)絕對值滿足一定大小才認(rèn)為其和頻率振蕩相關(guān)聯(lián),并且通過皮爾森系數(shù)的正負(fù)確定該有功出力是助增還是抑制頻率振蕩;最后,以標(biāo)準(zhǔn)差和皮爾森系數(shù)之積作為綜合相關(guān)系數(shù),對關(guān)聯(lián)振蕩機(jī)組助增或抑制頻率振蕩的能力進(jìn)行排序。實(shí)際算例分析結(jié)果表明,本文提出的算法能夠正確識(shí)別系統(tǒng)頻率振蕩,并且能夠快速定位關(guān)聯(lián)振蕩機(jī)組。

    猜你喜歡
    皮爾森標(biāo)準(zhǔn)差出力
    用Pro-Kin Line平衡反饋訓(xùn)練儀對早期帕金森病患者進(jìn)行治療對其動(dòng)態(tài)平衡功能的影響
    風(fēng)電場有功出力的EEMD特性分析
    要爭做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    風(fēng)電場群出力的匯聚效應(yīng)分析
    電測與儀表(2015年5期)2015-04-09 11:31:12
    對于平均差與標(biāo)準(zhǔn)差的數(shù)學(xué)關(guān)系和應(yīng)用價(jià)值比較研究
    數(shù)字翹楚皮爾森:忍過100多次整形的女軍人
    有夢的青春不易“殘”
    醫(yī)學(xué)科技論文中有效數(shù)字的確定
    談數(shù)據(jù)的變化對方差、標(biāo)準(zhǔn)差的影響
    日本91视频免费播放| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品自拍成人| 国产亚洲欧美在线一区二区| 久久久欧美国产精品| 亚洲少妇的诱惑av| 欧美成人午夜精品| 下体分泌物呈黄色| 亚洲成人手机| 高潮久久久久久久久久久不卡| 美女扒开内裤让男人捅视频| 可以免费在线观看a视频的电影网站| 一级毛片 在线播放| 少妇人妻 视频| 国产精品国产三级专区第一集| 久久久久久久久久久久大奶| 男女免费视频国产| 亚洲国产精品一区二区三区在线| 99国产综合亚洲精品| 美女高潮到喷水免费观看| 亚洲av美国av| 久久av网站| 在线观看免费视频网站a站| 丁香六月天网| 国产在线一区二区三区精| 丝瓜视频免费看黄片| 黄色 视频免费看| 久久久久久久国产电影| 一区二区三区激情视频| 麻豆乱淫一区二区| 黄色视频在线播放观看不卡| 日本欧美国产在线视频| 岛国毛片在线播放| 亚洲,欧美,日韩| 国产xxxxx性猛交| 韩国精品一区二区三区| 黑人欧美特级aaaaaa片| e午夜精品久久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 精品国产乱码久久久久久小说| 99精品久久久久人妻精品| 麻豆国产av国片精品| 国产熟女午夜一区二区三区| 久久精品亚洲av国产电影网| 国产黄频视频在线观看| 免费人妻精品一区二区三区视频| 久久九九热精品免费| 18在线观看网站| 日韩大码丰满熟妇| 亚洲人成电影观看| 久久九九热精品免费| 后天国语完整版免费观看| 欧美精品一区二区免费开放| 亚洲精品日韩在线中文字幕| 婷婷色麻豆天堂久久| 亚洲五月婷婷丁香| 欧美精品一区二区免费开放| 男的添女的下面高潮视频| 色精品久久人妻99蜜桃| 婷婷色综合大香蕉| 一区二区三区精品91| 婷婷色麻豆天堂久久| 嫁个100分男人电影在线观看 | 在线观看国产h片| 国产一区二区 视频在线| 亚洲国产精品一区三区| 两个人看的免费小视频| 青青草视频在线视频观看| 久久99热这里只频精品6学生| 伦理电影免费视频| 久久国产精品人妻蜜桃| 黄色a级毛片大全视频| 国产日韩一区二区三区精品不卡| 欧美日韩视频高清一区二区三区二| 亚洲色图综合在线观看| 欧美久久黑人一区二区| 日本av手机在线免费观看| 午夜福利视频精品| netflix在线观看网站| 美国免费a级毛片| av国产久精品久网站免费入址| 一本久久精品| 久久精品亚洲熟妇少妇任你| 老司机影院成人| 一本综合久久免费| 欧美国产精品va在线观看不卡| 久久久久网色| 精品国产国语对白av| 丰满少妇做爰视频| www日本在线高清视频| 国产精品熟女久久久久浪| 亚洲欧美一区二区三区国产| 一级黄色大片毛片| 黄色视频在线播放观看不卡| 国产精品国产三级国产专区5o| 国产日韩欧美视频二区| 国产视频一区二区在线看| 中文乱码字字幕精品一区二区三区| 一本综合久久免费| 亚洲精品国产av成人精品| 精品福利观看| 曰老女人黄片| 狂野欧美激情性bbbbbb| 国产成人欧美| av在线老鸭窝| 两个人免费观看高清视频| 蜜桃在线观看..| 日韩,欧美,国产一区二区三区| 丝袜美腿诱惑在线| 久久精品成人免费网站| av福利片在线| 99热网站在线观看| 日韩中文字幕视频在线看片| 欧美中文综合在线视频| 两个人看的免费小视频| 91精品伊人久久大香线蕉| 超碰成人久久| 另类精品久久| 亚洲av在线观看美女高潮| 亚洲国产毛片av蜜桃av| 80岁老熟妇乱子伦牲交| 日韩,欧美,国产一区二区三区| 人体艺术视频欧美日本| 日韩av不卡免费在线播放| 嫁个100分男人电影在线观看 | 欧美国产精品va在线观看不卡| 久热爱精品视频在线9| 在线亚洲精品国产二区图片欧美| 亚洲国产毛片av蜜桃av| 国产又色又爽无遮挡免| 亚洲,欧美,日韩| 老司机影院毛片| 日韩一本色道免费dvd| 女警被强在线播放| 777米奇影视久久| 久久人人爽av亚洲精品天堂| 91九色精品人成在线观看| 97精品久久久久久久久久精品| 老汉色∧v一级毛片| 国产成人系列免费观看| 午夜福利一区二区在线看| 久久女婷五月综合色啪小说| 亚洲av欧美aⅴ国产| 欧美日韩福利视频一区二区| 日韩av免费高清视频| 亚洲欧美日韩高清在线视频 | 国产成人精品无人区| 99精品久久久久人妻精品| 久久久久视频综合| 日日夜夜操网爽| 国产精品99久久99久久久不卡| 中文字幕制服av| 久久精品熟女亚洲av麻豆精品| 国产在线视频一区二区| 日韩av不卡免费在线播放| 日本wwww免费看| 日韩av免费高清视频| 91精品伊人久久大香线蕉| 亚洲精品成人av观看孕妇| 精品人妻1区二区| 熟女av电影| 最近最新中文字幕大全免费视频 | 在现免费观看毛片| 国产欧美日韩一区二区三区在线| 国产精品久久久人人做人人爽| 亚洲精品美女久久av网站| 国产精品久久久久久人妻精品电影 | 国产在线免费精品| 欧美日韩亚洲综合一区二区三区_| 亚洲人成电影观看| 精品欧美一区二区三区在线| 国精品久久久久久国模美| 老司机在亚洲福利影院| 亚洲专区中文字幕在线| h视频一区二区三区| 亚洲中文av在线| 美女大奶头黄色视频| 免费不卡黄色视频| 亚洲情色 制服丝袜| 婷婷色综合www| 啦啦啦啦在线视频资源| 日韩电影二区| 亚洲av日韩精品久久久久久密 | 精品久久蜜臀av无| 美女福利国产在线| 99热全是精品| 在线观看免费视频网站a站| 丝袜美腿诱惑在线| 亚洲一卡2卡3卡4卡5卡精品中文| 高清av免费在线| 日日夜夜操网爽| 精品少妇一区二区三区视频日本电影| 午夜两性在线视频| 精品免费久久久久久久清纯 | 亚洲三区欧美一区| av天堂在线播放| 欧美老熟妇乱子伦牲交| 你懂的网址亚洲精品在线观看| 丝袜人妻中文字幕| 亚洲综合色网址| 免费一级毛片在线播放高清视频 | 一本久久精品| 一区二区三区激情视频| 欧美精品高潮呻吟av久久| 高清视频免费观看一区二区| 国产女主播在线喷水免费视频网站| 午夜老司机福利片| 国产高清videossex| 欧美国产精品va在线观看不卡| 自线自在国产av| 中文字幕人妻丝袜一区二区| 亚洲欧美精品自产自拍| 99久久精品国产亚洲精品| 青春草亚洲视频在线观看| 国产老妇伦熟女老妇高清| 日本五十路高清| 国产精品.久久久| 午夜av观看不卡| 亚洲精品国产色婷婷电影| 丝袜脚勾引网站| 国产免费又黄又爽又色| 午夜福利免费观看在线| 青草久久国产| 美女扒开内裤让男人捅视频| 少妇人妻 视频| 亚洲激情五月婷婷啪啪| 成人国产av品久久久| 纯流量卡能插随身wifi吗| 久久久精品免费免费高清| 大话2 男鬼变身卡| 国产97色在线日韩免费| 日日夜夜操网爽| 丝袜在线中文字幕| 尾随美女入室| 麻豆av在线久日| 午夜激情久久久久久久| 五月天丁香电影| 国产日韩欧美亚洲二区| 九草在线视频观看| 婷婷色综合大香蕉| 黄片播放在线免费| 亚洲精品第二区| 丁香六月天网| 色94色欧美一区二区| 亚洲天堂av无毛| 国产人伦9x9x在线观看| 日韩 亚洲 欧美在线| 国产在线观看jvid| 国产成人精品久久久久久| 操美女的视频在线观看| 精品高清国产在线一区| 纯流量卡能插随身wifi吗| 成人免费观看视频高清| 欧美日韩福利视频一区二区| 亚洲情色 制服丝袜| 国产精品一区二区精品视频观看| 精品视频人人做人人爽| 黄频高清免费视频| 少妇人妻 视频| 99国产精品99久久久久| 伊人久久大香线蕉亚洲五| 国产在线视频一区二区| 视频在线观看一区二区三区| 高清av免费在线| 在线 av 中文字幕| 国产日韩欧美在线精品| 久久国产精品大桥未久av| 欧美人与性动交α欧美精品济南到| 国产视频首页在线观看| 国产精品久久久久久精品电影小说| 少妇精品久久久久久久| 一级毛片黄色毛片免费观看视频| 日本wwww免费看| 久久 成人 亚洲| 亚洲精品国产av蜜桃| 国产在线免费精品| 国产精品一国产av| 亚洲国产精品一区三区| 999精品在线视频| 国产一区二区三区综合在线观看| 亚洲色图 男人天堂 中文字幕| 青青草视频在线视频观看| 性色av一级| 亚洲成av片中文字幕在线观看| 一区在线观看完整版| 啦啦啦中文免费视频观看日本| 久久人人97超碰香蕉20202| 夫妻性生交免费视频一级片| 久久久精品国产亚洲av高清涩受| 青青草视频在线视频观看| 中文字幕亚洲精品专区| 女性生殖器流出的白浆| 亚洲精品久久久久久婷婷小说| 亚洲欧洲精品一区二区精品久久久| 日韩一区二区三区影片| 欧美激情 高清一区二区三区| 日韩大码丰满熟妇| 黄色a级毛片大全视频| 亚洲精品国产av蜜桃| 视频区欧美日本亚洲| 女警被强在线播放| 亚洲成人手机| 精品人妻1区二区| 久久中文字幕一级| 成年女人毛片免费观看观看9 | 久久综合国产亚洲精品| 婷婷色av中文字幕| 午夜福利免费观看在线| 视频区欧美日本亚洲| 亚洲欧洲国产日韩| 丰满少妇做爰视频| 亚洲第一av免费看| 国产成人精品无人区| 日本猛色少妇xxxxx猛交久久| 久久人人爽人人片av| 国产欧美日韩一区二区三 | 午夜免费男女啪啪视频观看| 亚洲精品日韩在线中文字幕| 亚洲av电影在线观看一区二区三区| 狂野欧美激情性bbbbbb| 十八禁网站网址无遮挡| 久久精品久久久久久噜噜老黄| av电影中文网址| 晚上一个人看的免费电影| 亚洲av电影在线进入| 久久久久国产一级毛片高清牌| 久久久久网色| 国产精品香港三级国产av潘金莲 | 一本综合久久免费| 国产深夜福利视频在线观看| 视频区欧美日本亚洲| 国产成人av教育| 真人做人爱边吃奶动态| 91成人精品电影| 久久精品国产亚洲av高清一级| 国产精品三级大全| 久久久久久免费高清国产稀缺| 亚洲av日韩精品久久久久久密 | av一本久久久久| 亚洲精品中文字幕在线视频| 黄色视频不卡| 亚洲人成电影免费在线| 中文字幕精品免费在线观看视频| av天堂在线播放| 中文字幕人妻丝袜制服| 啦啦啦在线观看免费高清www| 欧美日韩av久久| 亚洲成人国产一区在线观看 | 九色亚洲精品在线播放| 免费在线观看完整版高清| av福利片在线| 免费一级毛片在线播放高清视频 | 天天影视国产精品| 大型av网站在线播放| 欧美精品av麻豆av| 超碰97精品在线观看| 黄色一级大片看看| 一本色道久久久久久精品综合| 欧美日韩精品网址| 久久亚洲国产成人精品v| 国产熟女午夜一区二区三区| 国产成人a∨麻豆精品| 亚洲成人免费av在线播放| 亚洲国产精品一区二区三区在线| 少妇的丰满在线观看| 秋霞在线观看毛片| 亚洲精品一二三| 日本一区二区免费在线视频| 不卡av一区二区三区| 亚洲男人天堂网一区| 亚洲伊人久久精品综合| 男人爽女人下面视频在线观看| 高清视频免费观看一区二区| 女人精品久久久久毛片| 亚洲精品一卡2卡三卡4卡5卡 | 美女国产高潮福利片在线看| h视频一区二区三区| 99香蕉大伊视频| 2021少妇久久久久久久久久久| 午夜福利,免费看| 99热网站在线观看| 天堂俺去俺来也www色官网| 大香蕉久久网| 亚洲综合色网址| 欧美 日韩 精品 国产| 男女免费视频国产| 高潮久久久久久久久久久不卡| 亚洲少妇的诱惑av| 亚洲av国产av综合av卡| 丝袜美足系列| 久久精品aⅴ一区二区三区四区| 99精品久久久久人妻精品| 女性被躁到高潮视频| 亚洲,欧美,日韩| 亚洲精品乱久久久久久| 激情视频va一区二区三区| 男女之事视频高清在线观看 | 久久人人爽人人片av| 成人免费观看视频高清| 国产精品三级大全| 韩国精品一区二区三区| 婷婷色麻豆天堂久久| 一区二区三区乱码不卡18| 免费在线观看完整版高清| 久久久久视频综合| 超碰97精品在线观看| 亚洲七黄色美女视频| 国产成人av激情在线播放| 日本欧美视频一区| 欧美精品一区二区免费开放| 老司机深夜福利视频在线观看 | 日本黄色日本黄色录像| 自拍欧美九色日韩亚洲蝌蚪91| 妹子高潮喷水视频| 日本a在线网址| 日日夜夜操网爽| 久热这里只有精品99| 十分钟在线观看高清视频www| 亚洲av综合色区一区| 免费女性裸体啪啪无遮挡网站| 国产一区二区 视频在线| 精品一区二区三卡| 国产极品粉嫩免费观看在线| 亚洲国产精品成人久久小说| 亚洲av日韩精品久久久久久密 | 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 久久人人爽av亚洲精品天堂| 久久国产精品男人的天堂亚洲| 国产熟女午夜一区二区三区| 精品人妻在线不人妻| 性少妇av在线| 另类精品久久| cao死你这个sao货| 亚洲欧美一区二区三区久久| 深夜精品福利| 欧美日韩亚洲国产一区二区在线观看 | 亚洲成av片中文字幕在线观看| 久久人妻熟女aⅴ| 亚洲色图 男人天堂 中文字幕| 老司机午夜十八禁免费视频| 亚洲精品乱久久久久久| 啦啦啦在线免费观看视频4| 国产视频一区二区在线看| 人妻一区二区av| 首页视频小说图片口味搜索 | 日韩一本色道免费dvd| 欧美乱码精品一区二区三区| 欧美激情 高清一区二区三区| 国产精品亚洲av一区麻豆| 午夜视频精品福利| 久久精品成人免费网站| 一级毛片我不卡| 亚洲精品美女久久久久99蜜臀 | 天天影视国产精品| 久久精品久久久久久噜噜老黄| 麻豆乱淫一区二区| 亚洲av日韩在线播放| 丝袜美足系列| 十八禁网站网址无遮挡| 多毛熟女@视频| 免费黄频网站在线观看国产| 咕卡用的链子| 亚洲伊人色综图| 一区在线观看完整版| 侵犯人妻中文字幕一二三四区| 少妇裸体淫交视频免费看高清 | 日韩一区二区三区影片| 精品免费久久久久久久清纯 | 激情视频va一区二区三区| 大型av网站在线播放| 欧美乱码精品一区二区三区| 国产在线免费精品| 黑人欧美特级aaaaaa片| 欧美变态另类bdsm刘玥| 大码成人一级视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲色图 男人天堂 中文字幕| 国产精品99久久99久久久不卡| 国产有黄有色有爽视频| av国产精品久久久久影院| 国产一卡二卡三卡精品| 亚洲综合色网址| 精品亚洲成国产av| 国产极品粉嫩免费观看在线| 黄色 视频免费看| 国产精品国产三级专区第一集| 下体分泌物呈黄色| 最近中文字幕2019免费版| 在线观看免费高清a一片| 国产精品二区激情视频| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 国产一区亚洲一区在线观看| 97在线人人人人妻| 成人黄色视频免费在线看| 天堂俺去俺来也www色官网| 男女下面插进去视频免费观看| 精品久久久精品久久久| 亚洲av日韩精品久久久久久密 | 黄色 视频免费看| 欧美精品一区二区免费开放| 老司机亚洲免费影院| 日韩熟女老妇一区二区性免费视频| 欧美大码av| 男女无遮挡免费网站观看| 国产97色在线日韩免费| 99久久人妻综合| 免费观看a级毛片全部| 日韩一区二区三区影片| 一级片免费观看大全| 男人爽女人下面视频在线观看| 亚洲av日韩在线播放| 不卡av一区二区三区| 一级毛片黄色毛片免费观看视频| 丁香六月欧美| 91精品三级在线观看| 久久午夜综合久久蜜桃| 亚洲人成网站在线观看播放| 午夜久久久在线观看| 国产一区有黄有色的免费视频| 少妇人妻久久综合中文| 韩国精品一区二区三区| bbb黄色大片| 18禁观看日本| 国产国语露脸激情在线看| 亚洲专区中文字幕在线| 中文字幕人妻熟女乱码| 亚洲国产欧美日韩在线播放| av福利片在线| 飞空精品影院首页| 欧美人与性动交α欧美精品济南到| 国产无遮挡羞羞视频在线观看| 青草久久国产| 一区二区三区激情视频| 中文欧美无线码| 亚洲精品在线美女| 久久九九热精品免费| 亚洲中文字幕日韩| 免费看十八禁软件| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 久久久久国产精品人妻一区二区| 男女边摸边吃奶| 一级a爱视频在线免费观看| a级毛片在线看网站| 97人妻天天添夜夜摸| 亚洲黑人精品在线| a 毛片基地| 免费看十八禁软件| 97人妻天天添夜夜摸| 1024视频免费在线观看| 久久国产精品影院| 这个男人来自地球电影免费观看| 欧美成人午夜精品| 1024视频免费在线观看| 国产精品国产三级国产专区5o| 免费人妻精品一区二区三区视频| 久久久精品区二区三区| 亚洲精品日韩在线中文字幕| 麻豆乱淫一区二区| 日本av免费视频播放| 每晚都被弄得嗷嗷叫到高潮| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 精品久久久久久电影网| 精品人妻一区二区三区麻豆| 成年女人毛片免费观看观看9 | 午夜精品国产一区二区电影| 国产国语露脸激情在线看| 亚洲,欧美,日韩| 久久人人爽av亚洲精品天堂| 亚洲久久久国产精品| 高潮久久久久久久久久久不卡| 精品人妻熟女毛片av久久网站| 一级毛片女人18水好多 | 欧美大码av| 宅男免费午夜| 国产成人av激情在线播放| 国产在线视频一区二区| 91精品国产国语对白视频| 久久久国产精品麻豆| 女人爽到高潮嗷嗷叫在线视频| 中文字幕人妻丝袜制服| 午夜激情av网站| 婷婷色麻豆天堂久久| 亚洲精品中文字幕在线视频| 久久精品久久久久久噜噜老黄| av天堂在线播放| 日韩中文字幕视频在线看片| 在线观看免费视频网站a站| 国产人伦9x9x在线观看| 亚洲成av片中文字幕在线观看| 日日摸夜夜添夜夜爱| 90打野战视频偷拍视频| 夜夜骑夜夜射夜夜干| 永久免费av网站大全| 一级黄色大片毛片| 国产一区二区三区av在线| www.自偷自拍.com| 精品欧美一区二区三区在线| 少妇被粗大的猛进出69影院| 看免费av毛片| 成在线人永久免费视频| 亚洲欧美中文字幕日韩二区| 一区二区三区精品91| 欧美成人精品欧美一级黄| 精品一区二区三区四区五区乱码 | 欧美激情高清一区二区三区| 曰老女人黄片| 大香蕉久久成人网| 女人被躁到高潮嗷嗷叫费观| 操美女的视频在线观看| 一区二区日韩欧美中文字幕| 天堂8中文在线网|