• <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)差的影響
    91av网站免费观看| 国产不卡一卡二| 欧美日韩亚洲国产一区二区在线观看 | 国产伦人伦偷精品视频| 亚洲欧美日韩高清在线视频| 久久天躁狠狠躁夜夜2o2o| 亚洲精品日韩av片在线观看 | 国产精品久久视频播放| 日韩欧美国产一区二区入口| xxx96com| 亚洲精品久久国产高清桃花| 97超视频在线观看视频| 精品国产亚洲在线| 国产亚洲精品久久久com| 无限看片的www在线观看| 岛国在线免费视频观看| 日本在线视频免费播放| 亚洲国产日韩欧美精品在线观看 | 日本成人三级电影网站| 国产精品98久久久久久宅男小说| 国产乱人伦免费视频| 天天一区二区日本电影三级| 一本综合久久免费| 看片在线看免费视频| 怎么达到女性高潮| 国产探花在线观看一区二区| 国产aⅴ精品一区二区三区波| 国产精品久久电影中文字幕| 国产日本99.免费观看| 女警被强在线播放| 午夜精品在线福利| 久久久久国内视频| 成人午夜高清在线视频| 欧美激情在线99| 久久香蕉国产精品| 美女黄网站色视频| 成人欧美大片| 女警被强在线播放| 操出白浆在线播放| 中文字幕精品亚洲无线码一区| 99在线视频只有这里精品首页| tocl精华| 国产真实乱freesex| 欧美色视频一区免费| 免费观看的影片在线观看| 亚洲aⅴ乱码一区二区在线播放| 精品免费久久久久久久清纯| 18禁美女被吸乳视频| 亚洲性夜色夜夜综合| 久久精品影院6| 国产精品,欧美在线| 久久久久久久久久黄片| 国产91精品成人一区二区三区| 国产精品av视频在线免费观看| 国产真人三级小视频在线观看| 波多野结衣高清无吗| 成人国产一区最新在线观看| 久久久久久人人人人人| 两个人的视频大全免费| 国产极品精品免费视频能看的| 国产午夜福利久久久久久| 国语自产精品视频在线第100页| 少妇熟女aⅴ在线视频| 中文字幕久久专区| 色尼玛亚洲综合影院| 内地一区二区视频在线| av黄色大香蕉| 午夜久久久久精精品| 欧美一级a爱片免费观看看| 欧美一级a爱片免费观看看| 香蕉久久夜色| 午夜免费男女啪啪视频观看 | 在线看三级毛片| 在线观看免费午夜福利视频| 欧美日韩一级在线毛片| 国产中年淑女户外野战色| 老司机午夜十八禁免费视频| 欧美黄色淫秽网站| 欧美国产日韩亚洲一区| 天堂网av新在线| 国产探花在线观看一区二区| 欧美日本视频| 欧美激情久久久久久爽电影| 又紧又爽又黄一区二区| 身体一侧抽搐| 波多野结衣高清无吗| 国产真实伦视频高清在线观看 | 欧美日韩一级在线毛片| 亚洲aⅴ乱码一区二区在线播放| 嫁个100分男人电影在线观看| 欧美成人性av电影在线观看| 69av精品久久久久久| a在线观看视频网站| 亚洲成人久久爱视频| 两人在一起打扑克的视频| 99热这里只有精品一区| 高清在线国产一区| 啦啦啦观看免费观看视频高清| netflix在线观看网站| 国语自产精品视频在线第100页| 亚洲电影在线观看av| 两个人的视频大全免费| 欧美av亚洲av综合av国产av| 中文字幕人成人乱码亚洲影| 18禁国产床啪视频网站| 变态另类丝袜制服| 国产黄色小视频在线观看| www.www免费av| 国产淫片久久久久久久久 | 国产伦人伦偷精品视频| 母亲3免费完整高清在线观看| 久9热在线精品视频| 国产视频内射| 日韩欧美在线二视频| 免费人成在线观看视频色| 国产精品亚洲美女久久久| 白带黄色成豆腐渣| 久久久久久人人人人人| 亚洲欧美日韩无卡精品| 在线免费观看的www视频| 亚洲性夜色夜夜综合| 嫩草影院入口| 成年女人毛片免费观看观看9| 精品久久久久久久久久免费视频| 天堂av国产一区二区熟女人妻| 在线十欧美十亚洲十日本专区| 国产精品一区二区三区四区免费观看 | 亚洲成av人片免费观看| 嫩草影院入口| 美女 人体艺术 gogo| 欧美国产日韩亚洲一区| 国产精品电影一区二区三区| 一个人看的www免费观看视频| 成人特级黄色片久久久久久久| av中文乱码字幕在线| 欧美黄色片欧美黄色片| 国产精品女同一区二区软件 | 久久久久国产精品人妻aⅴ院| 日本黄色视频三级网站网址| 国产乱人视频| 亚洲成人中文字幕在线播放| 欧美黑人欧美精品刺激| 69人妻影院| 亚洲人成电影免费在线| 亚洲人成网站在线播| 亚洲在线自拍视频| 特级一级黄色大片| 不卡一级毛片| 久久国产乱子伦精品免费另类| 精品国产美女av久久久久小说| 国产精品爽爽va在线观看网站| 18美女黄网站色大片免费观看| 亚洲欧美激情综合另类| 久久国产精品人妻蜜桃| 91在线观看av| 九九在线视频观看精品| 99国产极品粉嫩在线观看| 丰满人妻一区二区三区视频av | 精品熟女少妇八av免费久了| 狂野欧美激情性xxxx| 精品国内亚洲2022精品成人| 久久久久国产精品人妻aⅴ院| 夜夜看夜夜爽夜夜摸| 国产午夜精品久久久久久一区二区三区 | 色av中文字幕| 两个人看的免费小视频| 白带黄色成豆腐渣| 亚洲 国产 在线| 中文字幕人成人乱码亚洲影| 亚洲av免费在线观看| av视频在线观看入口| 亚洲国产色片| 最新美女视频免费是黄的| 又黄又粗又硬又大视频| 久久精品91蜜桃| 精华霜和精华液先用哪个| 久久久久久久久久黄片| 两个人视频免费观看高清| 无人区码免费观看不卡| 日韩欧美国产在线观看| 成人高潮视频无遮挡免费网站| ponron亚洲| 国产精品久久久人人做人人爽| 亚洲av电影在线进入| 18禁裸乳无遮挡免费网站照片| 中出人妻视频一区二区| 成人欧美大片| 九九在线视频观看精品| 搡老岳熟女国产| 国产一区二区激情短视频| 久久6这里有精品| 国产亚洲精品久久久com| 麻豆国产av国片精品| 国产日本99.免费观看| 国产高清videossex| 九色国产91popny在线| 日韩国内少妇激情av| 每晚都被弄得嗷嗷叫到高潮| 欧美在线黄色| 天堂√8在线中文| 成人午夜高清在线视频| 国产美女午夜福利| 免费看日本二区| 看片在线看免费视频| 亚洲av五月六月丁香网| 国产野战对白在线观看| 久久精品影院6| 老鸭窝网址在线观看| 亚洲第一欧美日韩一区二区三区| 真人做人爱边吃奶动态| 亚洲专区国产一区二区| 亚洲精品一区av在线观看| 成年女人永久免费观看视频| 亚洲精品影视一区二区三区av| 特大巨黑吊av在线直播| 免费一级毛片在线播放高清视频| 日韩中文字幕欧美一区二区| 国产一区在线观看成人免费| 午夜免费观看网址| 91在线精品国自产拍蜜月 | 国产精品99久久99久久久不卡| 国产精品久久久久久亚洲av鲁大| 99精品欧美一区二区三区四区| 国产精品 国内视频| 精品日产1卡2卡| 亚洲成av人片在线播放无| 精品一区二区三区人妻视频| 亚洲激情在线av| 18禁美女被吸乳视频| x7x7x7水蜜桃| 脱女人内裤的视频| 亚洲av二区三区四区| 日韩欧美精品v在线| 90打野战视频偷拍视频| 别揉我奶头~嗯~啊~动态视频| 最近最新中文字幕大全电影3| 又黄又爽又免费观看的视频| 日韩欧美精品v在线| 中文字幕熟女人妻在线| 亚洲国产色片| 怎么达到女性高潮| 不卡一级毛片| 午夜免费观看网址| 男女视频在线观看网站免费| 黄片大片在线免费观看| 欧美日韩福利视频一区二区| 99精品在免费线老司机午夜| 全区人妻精品视频| 热99在线观看视频| 久久6这里有精品| 国产成人欧美在线观看| 国产亚洲精品综合一区在线观看| 国产午夜福利久久久久久| www.色视频.com| 岛国在线免费视频观看| 美女高潮的动态| 小蜜桃在线观看免费完整版高清| 欧美日韩黄片免| 国产私拍福利视频在线观看| 日韩av在线大香蕉| 精品日产1卡2卡| 天堂网av新在线| 久久香蕉精品热| 色综合亚洲欧美另类图片| 欧美激情在线99| 首页视频小说图片口味搜索| 波多野结衣巨乳人妻| 亚洲国产高清在线一区二区三| 亚洲欧美日韩无卡精品| 我要搜黄色片| 国产单亲对白刺激| 成人高潮视频无遮挡免费网站| www日本黄色视频网| 丝袜美腿在线中文| 熟女少妇亚洲综合色aaa.| 岛国在线观看网站| 美女大奶头视频| 久久久久久久精品吃奶| 精品日产1卡2卡| 伊人久久精品亚洲午夜| 男女午夜视频在线观看| 听说在线观看完整版免费高清| 成年免费大片在线观看| 内射极品少妇av片p| 欧美在线黄色| 亚洲七黄色美女视频| 中出人妻视频一区二区| 国产探花极品一区二区| 91av网一区二区| 成年女人永久免费观看视频| 免费在线观看成人毛片| 欧美极品一区二区三区四区| 国产探花在线观看一区二区| 不卡一级毛片| 国产精品一及| 欧美xxxx黑人xx丫x性爽| 欧美色欧美亚洲另类二区| 9191精品国产免费久久| 国产三级在线视频| 久久久久久大精品| 中文字幕精品亚洲无线码一区| 无限看片的www在线观看| 啦啦啦免费观看视频1| 女人十人毛片免费观看3o分钟| 99热这里只有精品一区| 少妇高潮的动态图| 一区福利在线观看| 午夜免费激情av| 国产精品电影一区二区三区| av在线蜜桃| 国产精品三级大全| 久久久久性生活片| 午夜免费男女啪啪视频观看 | 五月伊人婷婷丁香| 99久久无色码亚洲精品果冻| 欧美日韩亚洲国产一区二区在线观看| 日本熟妇午夜| 91九色精品人成在线观看| 国产伦在线观看视频一区| 国内精品美女久久久久久| 亚洲无线在线观看| 国产中年淑女户外野战色| 亚洲无线观看免费| 五月伊人婷婷丁香| 国产精品 国内视频| 精品久久久久久久久久免费视频| 国内精品美女久久久久久| 女生性感内裤真人,穿戴方法视频| 热99re8久久精品国产| 一本综合久久免费| 国产真实伦视频高清在线观看 | 欧美+日韩+精品| 国产真人三级小视频在线观看| 国产精品久久久久久人妻精品电影| www.熟女人妻精品国产| 夜夜看夜夜爽夜夜摸| 亚洲精品粉嫩美女一区| 国产中年淑女户外野战色| 少妇高潮的动态图| 日韩欧美一区二区三区在线观看| 别揉我奶头~嗯~啊~动态视频| 熟妇人妻久久中文字幕3abv| 日本 av在线| 黄片大片在线免费观看| 免费一级毛片在线播放高清视频| 精品久久久久久,| 午夜福利在线观看吧| 国产又黄又爽又无遮挡在线| 动漫黄色视频在线观看| 女生性感内裤真人,穿戴方法视频| av女优亚洲男人天堂| 18+在线观看网站| 男人舔女人下体高潮全视频| 久久九九热精品免费| 午夜免费男女啪啪视频观看 | 一进一出好大好爽视频| 日韩亚洲欧美综合| 免费人成在线观看视频色| 人人妻人人澡欧美一区二区| 日韩欧美精品免费久久 | 国产亚洲精品一区二区www| avwww免费| av女优亚洲男人天堂| 又黄又粗又硬又大视频| 日本熟妇午夜| 美女大奶头视频| 一本一本综合久久| 午夜免费男女啪啪视频观看 | 少妇的逼好多水| 国产69精品久久久久777片| 老司机在亚洲福利影院| 亚洲熟妇中文字幕五十中出| www国产在线视频色| 久久香蕉精品热| av在线天堂中文字幕| 黄色丝袜av网址大全| 一区二区三区高清视频在线| 极品教师在线免费播放| 亚洲国产中文字幕在线视频| 国产99白浆流出| 亚洲狠狠婷婷综合久久图片| 国产精品永久免费网站| 国产乱人视频| 久久久久久久精品吃奶| 日日夜夜操网爽| 亚洲精品成人久久久久久| 成人无遮挡网站| 亚洲一区二区三区不卡视频| 国产蜜桃级精品一区二区三区| 黄色片一级片一级黄色片| 成人国产综合亚洲| 国产精品一区二区三区四区久久| 桃色一区二区三区在线观看| 黄色成人免费大全| 国产精品日韩av在线免费观看| 午夜福利成人在线免费观看| 日本黄大片高清| 欧美日韩乱码在线| 怎么达到女性高潮| 久久国产乱子伦精品免费另类| 母亲3免费完整高清在线观看| 免费在线观看日本一区| 神马国产精品三级电影在线观看| 国模一区二区三区四区视频| 日韩欧美在线乱码| 久久亚洲精品不卡| 别揉我奶头~嗯~啊~动态视频| 白带黄色成豆腐渣| 久久久久久久午夜电影| 亚洲久久久久久中文字幕| 国产精品综合久久久久久久免费| 亚洲av一区综合| 亚洲人成网站高清观看| 国产真实伦视频高清在线观看 | 国产伦精品一区二区三区四那| 少妇人妻一区二区三区视频| 最后的刺客免费高清国语| 9191精品国产免费久久| 亚洲,欧美精品.| 夜夜夜夜夜久久久久| 日本黄色视频三级网站网址| 我要搜黄色片| 夜夜躁狠狠躁天天躁| 精华霜和精华液先用哪个| 成人av在线播放网站| 成年版毛片免费区| 久久久国产精品麻豆| 午夜久久久久精精品| 午夜免费成人在线视频| 俺也久久电影网| 内射极品少妇av片p| 长腿黑丝高跟| 国产在视频线在精品| 人人妻,人人澡人人爽秒播| 午夜老司机福利剧场| 久久精品亚洲精品国产色婷小说| 99精品欧美一区二区三区四区| 欧美极品一区二区三区四区| 97碰自拍视频| 18禁裸乳无遮挡免费网站照片| 内射极品少妇av片p| 女生性感内裤真人,穿戴方法视频| www.熟女人妻精品国产| 成人三级黄色视频| 欧美色视频一区免费| 夜夜爽天天搞| 在线看三级毛片| 动漫黄色视频在线观看| 欧美成人性av电影在线观看| 国产伦精品一区二区三区四那| 免费av毛片视频| 免费在线观看成人毛片| 动漫黄色视频在线观看| 老鸭窝网址在线观看| 亚洲一区二区三区不卡视频| 少妇的逼好多水| 久久国产精品人妻蜜桃| 亚洲精品日韩av片在线观看 | 国产探花在线观看一区二区| 国产主播在线观看一区二区| 亚洲一区高清亚洲精品| 日韩欧美 国产精品| 一区二区三区免费毛片| 老司机深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 国产高清videossex| 97人妻精品一区二区三区麻豆| 一进一出抽搐动态| 欧美最新免费一区二区三区 | 一区福利在线观看| 午夜激情欧美在线| 夜夜看夜夜爽夜夜摸| 天天躁日日操中文字幕| 国产精品亚洲美女久久久| 亚洲内射少妇av| 最新美女视频免费是黄的| 久久99热这里只有精品18| 婷婷丁香在线五月| 亚洲真实伦在线观看| 亚洲成人久久爱视频| 国产99白浆流出| 嫩草影院入口| 少妇熟女aⅴ在线视频| 国产成人福利小说| 国产aⅴ精品一区二区三区波| 五月伊人婷婷丁香| 欧美3d第一页| 免费看a级黄色片| 欧美日韩福利视频一区二区| 国产老妇女一区| 亚洲精品影视一区二区三区av| 国产探花极品一区二区| 嫩草影院精品99| 国产激情欧美一区二区| 99精品欧美一区二区三区四区| 国内久久婷婷六月综合欲色啪| 99久久九九国产精品国产免费| a级毛片a级免费在线| 国内精品久久久久精免费| 97超视频在线观看视频| 久久久久精品国产欧美久久久| 欧美bdsm另类| 国产v大片淫在线免费观看| 精品久久久久久,| 欧美丝袜亚洲另类 | 国产69精品久久久久777片| 99视频精品全部免费 在线| 精品久久久久久久毛片微露脸| 国产久久久一区二区三区| 我的老师免费观看完整版| 国产av不卡久久| 中国美女看黄片| 观看免费一级毛片| 小说图片视频综合网站| 男女做爰动态图高潮gif福利片| 99国产精品一区二区蜜桃av| 99热这里只有精品一区| 日本五十路高清| 麻豆国产av国片精品| 亚洲成人免费电影在线观看| 色精品久久人妻99蜜桃| 亚洲av二区三区四区| 99久久无色码亚洲精品果冻| 欧美精品啪啪一区二区三区| 九九热线精品视视频播放| 小蜜桃在线观看免费完整版高清| 麻豆国产97在线/欧美| 免费在线观看日本一区| 天堂动漫精品| 国产色爽女视频免费观看| 亚洲激情在线av| 国产在视频线在精品| 香蕉av资源在线| 日本一本二区三区精品| 深夜精品福利| 精品久久久久久久久久免费视频| 美女黄网站色视频| 中文在线观看免费www的网站| 国产私拍福利视频在线观看| 久久伊人香网站| 国产欧美日韩一区二区三| 一本精品99久久精品77| 久久精品国产清高在天天线| 中文字幕人妻熟人妻熟丝袜美 | 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 欧美色欧美亚洲另类二区| 女人高潮潮喷娇喘18禁视频| 一级黄片播放器| 午夜免费成人在线视频| 美女高潮的动态| 少妇人妻一区二区三区视频| 国产单亲对白刺激| 18禁黄网站禁片免费观看直播| 成人亚洲精品av一区二区| 3wmmmm亚洲av在线观看| 亚洲久久久久久中文字幕| 亚洲真实伦在线观看| 国产视频内射| 日韩成人在线观看一区二区三区| 两人在一起打扑克的视频| 村上凉子中文字幕在线| 欧美一级毛片孕妇| 最新中文字幕久久久久| 内地一区二区视频在线| 久久这里只有精品中国| 制服人妻中文乱码| 国产成年人精品一区二区| 脱女人内裤的视频| 久久久久久国产a免费观看| 国产激情偷乱视频一区二区| 极品教师在线免费播放| 少妇熟女aⅴ在线视频| 亚洲黑人精品在线| 国产视频一区二区在线看| 男人和女人高潮做爰伦理| 国内精品美女久久久久久| 日韩有码中文字幕| av天堂在线播放| 色综合欧美亚洲国产小说| 国产高清激情床上av| 99久久成人亚洲精品观看| 精品乱码久久久久久99久播| 一级毛片女人18水好多| 亚洲熟妇熟女久久| 国产亚洲精品综合一区在线观看| 午夜亚洲福利在线播放| 欧美日韩瑟瑟在线播放| 午夜福利在线观看吧| 精品国产超薄肉色丝袜足j| 国产精品美女特级片免费视频播放器| 婷婷精品国产亚洲av在线| 国产v大片淫在线免费观看| 国产一区在线观看成人免费| svipshipincom国产片| 国产私拍福利视频在线观看| 久久国产精品人妻蜜桃| 九色成人免费人妻av| 人人妻人人看人人澡| 午夜亚洲福利在线播放| 欧美最新免费一区二区三区 | 哪里可以看免费的av片| 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| 国产精品一区二区三区四区免费观看 | 午夜福利在线观看免费完整高清在 | 小说图片视频综合网站| 久久九九热精品免费| 在线观看av片永久免费下载| 精品一区二区三区视频在线 | 久久天躁狠狠躁夜夜2o2o| 国产老妇女一区| 国产淫片久久久久久久久 | 中文在线观看免费www的网站| 久久99热这里只有精品18|