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

    考慮高斯激勵的多機電力系統(tǒng)的p階隨機穩(wěn)定性分析

    2020-04-10 12:34:00秦子璇盧占會魏軍強
    關(guān)鍵詞:兩區(qū)功角機系統(tǒng)

    秦子璇, 盧占會, 魏軍強

    (華北電力大學 數(shù)理學院,北京 102206)

    0 引 言

    當前全球能源系統(tǒng)面臨的主要問題有能源短缺、環(huán)境惡化以及能源利用效率較低等。以可再生能源逐步替代化石能源,實現(xiàn)能源轉(zhuǎn)型,建設(shè)清潔低碳、安全高效的新一代能源系統(tǒng)是全世界新一輪能源革命的核心目標[1]。與廣泛使用的傳統(tǒng)能源相比,新能源是指基于新技術(shù)開發(fā)和利用的非常規(guī)能源,包括風能、太陽能、生物質(zhì)能、地熱能和海洋能等。近年來,新能源的使用量一直在快速增長,由于新能源發(fā)電的間歇性和隨機波動性,以及新能源接入系統(tǒng)可能出現(xiàn)的低頻和高頻問題等,在將新能源大規(guī)模集成到電網(wǎng)中時,都可能會對電網(wǎng)的穩(wěn)定性產(chǎn)生重大影響[3-5]。例如,風力與光伏發(fā)電,都與風和光的大小和強弱有關(guān),當大規(guī)模風電站集成到電網(wǎng)中時,由于風能的隨機性,風電將隨著風速的變化而隨時變化,這將影響電能的實時平衡,進而影響電網(wǎng)的安全穩(wěn)定運行。

    國內(nèi)外許多學者對電力系統(tǒng)的隨機穩(wěn)定性進行了一些研究,研究成果見諸于不同文獻[6-16]。其中,文獻[6]基于隨機微積分理論,提出了一種新的隨機暫態(tài)穩(wěn)定性評估框架。文獻[7]對擬哈密頓系統(tǒng)利用隨機平均法,對阻尼系數(shù)、激勵強度等對隨機系統(tǒng)的暫態(tài)穩(wěn)定性的影響進行分析。文獻[8]針對單機無窮大系統(tǒng)下的擬哈密頓系統(tǒng)的隨機平均法的一次推廣,同時將單機系統(tǒng)的理論推廣到了多機系統(tǒng),得到了隨機激勵下的多機系統(tǒng)能量函數(shù)的一維擴散過程公式。文獻[9]通過隨機Lyapunov穩(wěn)定性分析,將電力系統(tǒng)瞬態(tài)能量函數(shù)法用于隨機電力系統(tǒng)的穩(wěn)定性。文獻[10]提出了電力系統(tǒng)的均值和均方穩(wěn)定性的概念,1階穩(wěn)定性為均值穩(wěn)定,2階穩(wěn)定性為均方穩(wěn)定,均值穩(wěn)定意味著,一個隨機輸入所激發(fā)的系統(tǒng)響應(yīng)的均值是有界的;均方穩(wěn)定意味著,一個隨機輸入所激發(fā)的系統(tǒng)響應(yīng)的均方差是有界的,并采用Euler-Maruyama(EM)數(shù)值方法得到電力系統(tǒng)功角曲線的軌跡,同時分析了不同激勵強度下功角曲線的不同。文獻[11] 利用Heun數(shù)值方法對單機無窮大系統(tǒng)進行了仿真模擬,分別分析了計算步長、激勵強度以及激勵步長對功角曲線的影響。文獻[12]采用It公式和It等距法推導出系統(tǒng)狀態(tài)方程的解析解,同時提出利用矩穩(wěn)定性的定義來描述隨機穩(wěn)定性,給出了隨機激勵下系統(tǒng)小擾動的穩(wěn)定性判據(jù),并通過雙機無窮大系統(tǒng)的仿真驗證了該方法的有效性。文獻[13-14]研究了區(qū)間穩(wěn)定性,利用Lyapunov函數(shù)和矩陣不等式等數(shù)學工具,分析了異步風力機的小干擾區(qū)間模型,給出了系統(tǒng)的區(qū)間穩(wěn)定性。文獻[15]利用Lyapunov方法和其它相關(guān)理論,證明了一類隨機線性系統(tǒng)的p階穩(wěn)定性定理,其中3階穩(wěn)定性和4階穩(wěn)定性的意義與偏度和峰度有關(guān),峰度值越大,概率密度曲線的尾端越厚,這意味著極值出現(xiàn)的概率更大;相反,如果峰度值很小,則概率分布集中,這意味著極值出現(xiàn)的概率很小,因此隨機系統(tǒng)的穩(wěn)定性更好。但基于隨機微分方程理論研究多機電力系統(tǒng)的p階穩(wěn)定性問題尚未見報道,本文將對多機電力系統(tǒng)在隨機激勵下的p階穩(wěn)定性問題進行研究擴展,擴展思路是先將多機系統(tǒng)進行等效轉(zhuǎn)換,將其等效為單機系統(tǒng),利用已有理論知識對其穩(wěn)定性進行檢驗。

    本文利用隨機微分方程理論,對原多機系統(tǒng)以及由擴展等面積法得到的多機電力系統(tǒng)的等效模型進行線性化處理。選取四機兩區(qū)系統(tǒng)作為實際算例,利用EM數(shù)值方法研究電力系統(tǒng)功角、轉(zhuǎn)速等隨時間變化的趨勢,通過蒙特卡羅法做出等效轉(zhuǎn)換前后電力系統(tǒng)的功角的概率密度函數(shù)曲線,利用其分析等效轉(zhuǎn)換前后系統(tǒng)穩(wěn)定性是否發(fā)生變化。將單機無窮大系統(tǒng)的p階穩(wěn)定性擴展到多機系統(tǒng)的p階穩(wěn)定性。

    1 多機電力系統(tǒng)的隨機建模

    1.1 多機電力系統(tǒng)的隨機模型

    考慮具有高斯型白噪聲的矢量隨機微分方程:

    dX(t)=AX(t)dt+QdB(t),t∈[0,+∞)

    (1)

    式中:X(t)=(X1(t),X2(t),…,Xn(t))T表示n維隨機變量,B(t)=(B1(t),B2(t),…,Bn(t))T表示n維維納過程,dB(t)/dt=W(t),W(t)為高斯過程,A和Q為常數(shù)矩陣。

    (2)

    1.2 多機系統(tǒng)的雙機等效隨機模型

    (3)

    根據(jù)擴展等面積法,將多機系統(tǒng)的失穩(wěn)看作是一臺或是多臺發(fā)電機形成的一個機群相對系統(tǒng)的其余部分失去同步,因此將受擾部分機群稱作S群,其余部分稱為A群,推導出了多機系統(tǒng)的雙機等效模型[7]:

    (4)

    根據(jù)隨機微分方程理論,對(4)式進行線性化處理,可得:

    (5)

    2 多機系統(tǒng)的p階穩(wěn)定性

    2.1 p階穩(wěn)定性基礎(chǔ)知識

    關(guān)于隨機激勵下電力系統(tǒng)的p階穩(wěn)定性,有如下結(jié)論。

    定義1[15].如果隨機微分方程(1)的解過程X(t)滿足以下條件:

    (6)

    則系統(tǒng)(1)具有p階穩(wěn)定性。對于所有的p> 0,其中c> 0,c是常數(shù)。

    在上述定義下,可得出系統(tǒng)(1)的穩(wěn)定性條件。

    引理1[15]. 對于給定的高斯型隨機激勵下的線性常系數(shù)隨機微分系統(tǒng)(1),初始值X(0)=X0是有界的并且與維納過程B(t)獨立。 當矩陣A的所有特征值的實部小于零時,隨機微分系統(tǒng)(1)具有p階穩(wěn)定性。

    2.2 Euler-Maruyama數(shù)值方法

    對于大多數(shù)隨機微分方程來說,其解析解是無法求得的,可以求出解析解的隨機微分方程大多為線性隨機微分方程。一些隨機微分方程雖可求得解析解,但解析解形式復雜,處理起來非常不方便。因此,在實際應(yīng)用中,往往通過數(shù)值積分的方法獲得解過程的軌跡,利用解的軌跡無限逼近精確解即可。

    EM 算法是目前最簡單的求解隨機微分方程的數(shù)值解法。設(shè)隨機微分方程的解過程是X(t),對于某個正整數(shù)N,有Δt=(T-t0)/N,τj=jΔt,j=0,1,2,…,N,Xj=X(τj)。EM數(shù)值方法對應(yīng)的差分迭代格式為

    Xj+1=Xj+f(Xj,τj)Δt+G(Xj,τj)ΔWj

    (7)

    3 案例研究

    3.1 四機兩區(qū)系統(tǒng)

    四機兩區(qū)電力系統(tǒng)的示意圖如圖1所示[19]。

    通過對四機電力系統(tǒng)進行簡化,可得四機兩區(qū)系統(tǒng)的伊藤隨機微分方程為[8]

    圖1 四機兩區(qū)電力系統(tǒng)示意圖Fig.1 Schematic diagram of four-machine two-area power system

    式中:第一臺發(fā)電機的方程為

    (8)

    對(8)式進行線性化處理,可得:

    (9)

    則(9)式的系數(shù)矩陣為

    E1E3B13cos(δ10-δ30)-E1E4B14cos(δ10-δ40))。

    同理,可以寫出四機兩區(qū)系統(tǒng)中其余三臺發(fā)電機對應(yīng)的伊藤隨機微分方程,同時可寫出其系數(shù)矩陣A2、A3、A4,則四機兩區(qū)系統(tǒng)的系數(shù)矩陣為

    3.2 數(shù)值分析

    通過對電力系統(tǒng)四機兩區(qū)模型進行EM數(shù)值分析,觀察發(fā)電機功角、發(fā)電機轉(zhuǎn)速在給定初始值的情況下,隨時間的變化的趨勢,對電力系統(tǒng)的p階穩(wěn)定性進行一系列的分析驗證。四機兩區(qū)系統(tǒng)的各個參數(shù)值分別為[8]:σ1=0.05,σ2=0.06,σ3=0.07,σ4=0.08,M1=13 s,M2=13 s,M3=12.35 s,M4=12.35 s,E1=1.03,E2=1.01,E3=1.03,E4=1,P1=0.998 119 479 653 3,P2=0.319 842 206 470 3,P3=0.241 493 636 334 7,P4=1.076 468 049 788 9,B12=B21=1.49,B13=B31=0.62,B14=B41=0.7,B23=B32=0.7,B24=B42=0.76,B43=B34=1.49,D1=0.5,D2=0.55,D3=0.6,D4=0.65,ωN=2·π·60,系統(tǒng)的穩(wěn)定平衡點為δ1 s-δ4 s=37.2°,δ2s-δ4 s=27.5°,δ3s-δ4 s=10.2°。

    利用已知參數(shù),對等效轉(zhuǎn)換后的多機電力系統(tǒng)進行數(shù)值模擬,得到圖2-圖4的三幅圖,分別是系統(tǒng)功角曲線、轉(zhuǎn)速以及500條系統(tǒng)功角模擬曲線的均值與時間t的關(guān)系圖。

    通過應(yīng)用EM數(shù)值計算方法對系統(tǒng)(5)進行分析,在xt0=(0,0)T的初始條件下,30 s內(nèi)的運算結(jié)果如圖2和圖3所示。

    圖3 ωSA隨時間t的波動圖Fig.3 Movement of ωSA along the time t

    從圖4可以看出,在隨機激勵的影響下,功角曲線在初值δSA=0 rad附近波動。在統(tǒng)計學中常用樣本均值來估計總體均值,選擇500條樣本路徑模擬的均值來估計總體的變動情況,其中路徑1和路徑2是500條樣本模擬中的兩條。單獨的每條樣本模擬曲線隨時間的波動都是比較大的,但是500條路徑的樣本均值則波動較小,基本接近穩(wěn)態(tài)初始值δSA=0 rad。

    圖4 500條樣本路徑模擬的均值Fig.4 Average angle over 500 paths

    利用已有參數(shù),可以分別計算出等效轉(zhuǎn)換前后兩個模型的系數(shù)矩陣A的特征值。其中,原四機兩區(qū)系統(tǒng)的系數(shù)矩陣A的特征值為

    λ1=-0.000 192 307 7+0.882 978 3617j

    λ4=-0.000 192 307 7-0.882 978 3617j

    λ3=-0.000 211 584 6+0.917 693 818 5j

    λ4=-0.000 211 584 6-0.917 693 8185j

    λ5=-0.000 242 949 8+0.923 161 262 2j

    λ6=-0.000 242 949 8-0.923 161 262 2j

    λ7=-0.052 631 587 9

    λ8=-0.037 699 111 2

    根據(jù)引理1的描述,當隨機微分方程系統(tǒng)的系數(shù)矩陣A的特征值實部均小于零時,可以推斷出系統(tǒng)是p階穩(wěn)定的。因此原四機兩區(qū)系統(tǒng)是p階穩(wěn)定的。

    同理可以求出等效轉(zhuǎn)換后的系數(shù)矩陣A的特征值:

    λ1=-0.000 226 824 5+0.820 252 071 9j

    λ2=-0.000 226 824 5-0.820 252 071 9j

    由計算結(jié)果可以看出,系數(shù)矩陣A的特征值實部均小于零。則等效轉(zhuǎn)換前后系統(tǒng)的p階穩(wěn)定性沒有發(fā)生變化。對于多機系統(tǒng)的p階穩(wěn)定性,均可以采用這種方式進行驗證。首先將需要檢驗是否具有p階穩(wěn)定性的多機系統(tǒng),利用擴展等面積法等[20],通過進行等效轉(zhuǎn)換的方式,將多機系統(tǒng)轉(zhuǎn)換成雙機或單機系統(tǒng),其次利用已有定理對等效轉(zhuǎn)換后的系統(tǒng)進行判定,最后根據(jù)等效轉(zhuǎn)換不改變系統(tǒng)穩(wěn)定性這一結(jié)論,可以推斷出原多機系統(tǒng)是否為p階穩(wěn)定。

    圖5是等效轉(zhuǎn)換后t=1 s時多機電力系統(tǒng)的功角δ的概率密度函數(shù)圖像,可以看出δ的概率密度函數(shù)曲線是基本符合正態(tài)分布的趨勢。對于3階、4階穩(wěn)定性的判定,通常也可以使用偏度和峰度這一定義[15]。根據(jù)表1中仿真模擬的結(jié)果可以看出,盡管功角δ在各時刻的偏度和峰度不完全接近0和3,但總體的偏度和峰度的模擬效果還是很接近理想狀態(tài)的。這也驗證了前文對于多機系統(tǒng)穩(wěn)定性的判定。

    圖5 δ的概率密度函數(shù)Fig.5 Probability density function of δ

    表 1 t時刻下δ的偏度和峰度值

    3.3 等效轉(zhuǎn)換前后系統(tǒng)功角的概率密度

    通過利用蒙特卡羅方法,對等效轉(zhuǎn)換前后的兩個系統(tǒng)在不同時刻下的概率密度函數(shù)的大小進行對比。得到如圖6-圖9所示的,4個不同時刻下的基于等效轉(zhuǎn)換后的系統(tǒng)和原多機系統(tǒng)的四幅圖。對于原始多機系統(tǒng)而言,功角δ是第一臺發(fā)電機功角,對于等效轉(zhuǎn)換后的系統(tǒng)而言為ΔδSA。

    圖6 時間t為1 s時δ的概率密度函數(shù)Fig.6 Probability density function of δ when time t is 1 s

    由圖6-圖9可以看出,不同時刻下基于原多機系統(tǒng)和轉(zhuǎn)換后系統(tǒng)的功角δ的概率密度函數(shù)曲線的一致性比較好,由此可以判定出等效轉(zhuǎn)化前后系統(tǒng)的穩(wěn)定性基本不會發(fā)生變化,那么多機電力系統(tǒng)的p階穩(wěn)定性問題,可以通過等效轉(zhuǎn)換的方式轉(zhuǎn)換為單機或雙機系統(tǒng)進行驗證。

    圖7 時間t為2 s時δ的概率密度函數(shù)Fig.7 Probability density function of δ when time t is 2 s

    圖8 時間t為5 s時δ的概率密度函數(shù)Fig.8 Probability density function of δ when time t is 5 s

    圖9 時間t為10 s時δ的概率密度函數(shù)Fig.9 Probability density function of δ when time t is 10 s

    4 結(jié) 論

    本文基于隨機微分方程理論、EM數(shù)值方法以及蒙特卡羅法對多機電力系統(tǒng)的p階穩(wěn)定性進行了分析。在多機電力系統(tǒng)的系數(shù)矩陣形式復雜或參數(shù)不易獲取等情況下,可以將其等效成單機或雙機模型,利用p階穩(wěn)定性定理對等效之后的模型的系數(shù)矩陣進行分析。通過對實際四機兩區(qū)系統(tǒng)的EM數(shù)值模擬的結(jié)果來看,電力系統(tǒng)是p階穩(wěn)定的。根據(jù)蒙特卡羅法,對電力系統(tǒng)的功角的概率密度函數(shù)曲線分析,發(fā)現(xiàn)等效轉(zhuǎn)換前后電力系統(tǒng)的功角穩(wěn)定性不發(fā)生變化。未來的研究工作或?qū)⒓杏诳紤]乘性噪聲下的電力系統(tǒng)的p階穩(wěn)定性。

    猜你喜歡
    兩區(qū)功角機系統(tǒng)
    虛擬調(diào)速器對VSG暫態(tài)功角穩(wěn)定影響機理分析
    基于改進 shapelet 挖掘的風電并網(wǎng)系統(tǒng)暫態(tài)功角穩(wěn)定評估
    能源工程(2019年6期)2019-12-02 01:58:20
    手持式觸摸測試機對閘機系統(tǒng)的維護研究
    基于功角測量和等面積法則的發(fā)電機暫態(tài)穩(wěn)定在線判別研究
    電子測試(2017年12期)2017-12-18 06:35:33
    經(jīng)濟、可靠的自動開關(guān)機系統(tǒng)
    基于PSS/E風電并網(wǎng)系統(tǒng)的靜態(tài)功角穩(wěn)定性分析
    京石高速公路自助發(fā)卡機系統(tǒng)的設(shè)計與應(yīng)用
    后置式自動發(fā)卡機系統(tǒng)應(yīng)用
    欧美日韩视频高清一区二区三区二| 欧美国产精品va在线观看不卡| 麻豆乱淫一区二区| 国产91精品成人一区二区三区 | 捣出白浆h1v1| 国产精品成人在线| www.熟女人妻精品国产| 欧美国产精品va在线观看不卡| av天堂在线播放| 男女床上黄色一级片免费看| 男的添女的下面高潮视频| 两人在一起打扑克的视频| av国产精品久久久久影院| 99热全是精品| a级片在线免费高清观看视频| 亚洲av欧美aⅴ国产| 女人精品久久久久毛片| 人妻人人澡人人爽人人| 久久久国产欧美日韩av| 肉色欧美久久久久久久蜜桃| 精品久久蜜臀av无| www日本在线高清视频| 成年av动漫网址| 精品少妇黑人巨大在线播放| 亚洲色图 男人天堂 中文字幕| 欧美精品av麻豆av| 亚洲熟女毛片儿| 久久精品熟女亚洲av麻豆精品| 欧美日韩亚洲综合一区二区三区_| 国产成人一区二区三区免费视频网站 | 亚洲情色 制服丝袜| 久久国产精品影院| 亚洲成人国产一区在线观看 | 久久久久久人人人人人| 欧美黄色片欧美黄色片| 亚洲色图 男人天堂 中文字幕| av一本久久久久| 日日摸夜夜添夜夜爱| 免费一级毛片在线播放高清视频 | 亚洲精品国产av成人精品| 欧美日韩精品网址| 国产日韩欧美视频二区| 欧美 日韩 精品 国产| 成人国产一区最新在线观看 | 欧美激情高清一区二区三区| 亚洲人成77777在线视频| 肉色欧美久久久久久久蜜桃| 色精品久久人妻99蜜桃| 久9热在线精品视频| 日韩制服骚丝袜av| 一边亲一边摸免费视频| 国产精品一区二区免费欧美 | 欧美日韩精品网址| 每晚都被弄得嗷嗷叫到高潮| 乱人伦中国视频| 精品人妻在线不人妻| 国产在线视频一区二区| 99国产精品一区二区蜜桃av | 国产亚洲av片在线观看秒播厂| 国产精品九九99| 一区在线观看完整版| 久久99热这里只频精品6学生| 国产精品一二三区在线看| 黄片小视频在线播放| 最新的欧美精品一区二区| 一级毛片电影观看| 啦啦啦啦在线视频资源| videosex国产| bbb黄色大片| 人人澡人人妻人| 亚洲精品第二区| 欧美变态另类bdsm刘玥| 电影成人av| 久久久精品94久久精品| 啦啦啦啦在线视频资源| 午夜福利在线免费观看网站| 99热网站在线观看| 久久精品国产综合久久久| 新久久久久国产一级毛片| 黑丝袜美女国产一区| 日本wwww免费看| 热re99久久国产66热| 亚洲精品在线美女| 国产成人精品在线电影| 久久久久国产一级毛片高清牌| 国产精品久久久久成人av| tube8黄色片| 国产成人精品久久久久久| 日韩中文字幕视频在线看片| 一级a爱视频在线免费观看| 91麻豆精品激情在线观看国产 | 久久人人爽av亚洲精品天堂| 韩国精品一区二区三区| 亚洲av电影在线进入| 亚洲av成人精品一二三区| 国产成人精品久久二区二区免费| 在线观看一区二区三区激情| 只有这里有精品99| 欧美性长视频在线观看| 国产主播在线观看一区二区 | 晚上一个人看的免费电影| 亚洲中文av在线| 嫁个100分男人电影在线观看 | 亚洲第一av免费看| 搡老岳熟女国产| 午夜福利免费观看在线| 美女主播在线视频| 啦啦啦中文免费视频观看日本| 欧美 日韩 精品 国产| 亚洲,欧美,日韩| 午夜福利免费观看在线| 国产亚洲精品久久久久5区| 午夜av观看不卡| 两个人看的免费小视频| 欧美成人午夜精品| 免费黄频网站在线观看国产| av在线老鸭窝| 国产精品国产av在线观看| 亚洲中文av在线| 亚洲精品久久久久久婷婷小说| 久久久久久亚洲精品国产蜜桃av| 久久人人爽av亚洲精品天堂| 女性被躁到高潮视频| 亚洲国产成人一精品久久久| 国产午夜精品一二区理论片| 男女之事视频高清在线观看 | av国产精品久久久久影院| 国产黄色视频一区二区在线观看| 国产亚洲精品第一综合不卡| 岛国毛片在线播放| 中文欧美无线码| 在线av久久热| 亚洲欧美激情在线| 国产欧美日韩精品亚洲av| 成人手机av| 国精品久久久久久国模美| 国产日韩欧美亚洲二区| 99国产精品一区二区蜜桃av | 日韩熟女老妇一区二区性免费视频| 色婷婷av一区二区三区视频| 97精品久久久久久久久久精品| 99re6热这里在线精品视频| 亚洲人成电影免费在线| 欧美乱码精品一区二区三区| 国产91精品成人一区二区三区 | 欧美亚洲 丝袜 人妻 在线| 国产成人av激情在线播放| 久久久久久久国产电影| 免费观看人在逋| 欧美少妇被猛烈插入视频| 999精品在线视频| 视频区欧美日本亚洲| 成人国语在线视频| 欧美精品亚洲一区二区| 91麻豆精品激情在线观看国产 | 精品少妇内射三级| 国产一区有黄有色的免费视频| 电影成人av| 看免费av毛片| 国产野战对白在线观看| 国产女主播在线喷水免费视频网站| 在线观看国产h片| 建设人人有责人人尽责人人享有的| 美女午夜性视频免费| 国产精品一区二区在线观看99| 中国美女看黄片| bbb黄色大片| 国产日韩一区二区三区精品不卡| 一级毛片我不卡| 美国免费a级毛片| 又紧又爽又黄一区二区| 亚洲精品国产av蜜桃| 在线观看一区二区三区激情| 日韩熟女老妇一区二区性免费视频| 国产一区二区三区综合在线观看| 老熟女久久久| 国产91精品成人一区二区三区 | 亚洲精品国产av成人精品| 水蜜桃什么品种好| 亚洲av在线观看美女高潮| 18禁国产床啪视频网站| 国产欧美亚洲国产| 国产成人精品在线电影| 男人操女人黄网站| av一本久久久久| 久久久久国产精品人妻一区二区| 午夜福利视频在线观看免费| 成人亚洲欧美一区二区av| 亚洲欧洲国产日韩| 亚洲色图综合在线观看| 亚洲熟女毛片儿| 午夜免费男女啪啪视频观看| 高清黄色对白视频在线免费看| 国产成人欧美在线观看 | 国产91精品成人一区二区三区 | 波多野结衣一区麻豆| 亚洲欧美一区二区三区久久| 免费在线观看完整版高清| 久久久精品国产亚洲av高清涩受| 国产视频一区二区在线看| 五月开心婷婷网| 国产亚洲欧美精品永久| 亚洲精品中文字幕在线视频| 电影成人av| 国产成人a∨麻豆精品| 国产xxxxx性猛交| 日本色播在线视频| 一本久久精品| 国产91精品成人一区二区三区 | 欧美日韩福利视频一区二区| 又紧又爽又黄一区二区| 久久女婷五月综合色啪小说| 色视频在线一区二区三区| 欧美精品av麻豆av| 18禁观看日本| 免费看十八禁软件| 99久久99久久久精品蜜桃| 国产欧美日韩一区二区三 | √禁漫天堂资源中文www| 亚洲专区中文字幕在线| a级毛片在线看网站| 中文字幕av电影在线播放| 欧美97在线视频| 国产一区二区三区av在线| 在现免费观看毛片| 19禁男女啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 国产亚洲精品第一综合不卡| bbb黄色大片| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| 久久热在线av| 美女中出高潮动态图| 日本午夜av视频| 成年人午夜在线观看视频| 久久鲁丝午夜福利片| 少妇的丰满在线观看| 天堂俺去俺来也www色官网| 黄频高清免费视频| 久久精品久久精品一区二区三区| 国产一卡二卡三卡精品| 国产成人免费观看mmmm| 另类精品久久| 免费看十八禁软件| 婷婷色综合www| 捣出白浆h1v1| 国产精品一区二区免费欧美 | 黄片小视频在线播放| a级片在线免费高清观看视频| 乱人伦中国视频| 日韩视频在线欧美| 尾随美女入室| 成人手机av| 国产日韩一区二区三区精品不卡| 日韩 亚洲 欧美在线| 国产精品久久久av美女十八| cao死你这个sao货| 久久久久精品人妻al黑| 亚洲成人国产一区在线观看 | 一级毛片 在线播放| 黄色怎么调成土黄色| 久久久精品区二区三区| 亚洲av日韩在线播放| 一区二区三区四区激情视频| 欧美+亚洲+日韩+国产| 欧美精品高潮呻吟av久久| av不卡在线播放| 亚洲av综合色区一区| 欧美97在线视频| 国产福利在线免费观看视频| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| 在线观看免费日韩欧美大片| 精品福利观看| 日本a在线网址| 99久久综合免费| 国产91精品成人一区二区三区 | 一级片免费观看大全| 天天躁日日躁夜夜躁夜夜| 深夜精品福利| av天堂久久9| 天天躁狠狠躁夜夜躁狠狠躁| 日韩一卡2卡3卡4卡2021年| 欧美亚洲日本最大视频资源| 最近中文字幕2019免费版| 一本色道久久久久久精品综合| 激情五月婷婷亚洲| 欧美 亚洲 国产 日韩一| 91字幕亚洲| 亚洲人成网站在线观看播放| 国产精品一区二区在线观看99| 国产一区二区在线观看av| 亚洲国产欧美网| 国产xxxxx性猛交| 亚洲精品自拍成人| 一级毛片女人18水好多 | 性少妇av在线| 国产成人免费无遮挡视频| 欧美日韩精品网址| 国产亚洲精品第一综合不卡| 久久久久精品人妻al黑| 亚洲精品自拍成人| 国产日韩一区二区三区精品不卡| 晚上一个人看的免费电影| 天天添夜夜摸| 大片电影免费在线观看免费| 国产黄色视频一区二区在线观看| 男女床上黄色一级片免费看| 日本vs欧美在线观看视频| a级毛片黄视频| 午夜福利,免费看| 国产精品秋霞免费鲁丝片| 中文乱码字字幕精品一区二区三区| 免费看不卡的av| 黄色片一级片一级黄色片| 一本大道久久a久久精品| h视频一区二区三区| 成人手机av| 丝袜美足系列| 精品人妻熟女毛片av久久网站| 欧美大码av| 精品一区二区三卡| 久久久久久亚洲精品国产蜜桃av| 女人高潮潮喷娇喘18禁视频| 亚洲免费av在线视频| www.自偷自拍.com| 久9热在线精品视频| 久久久精品区二区三区| 视频在线观看一区二区三区| 秋霞在线观看毛片| 电影成人av| 美女国产高潮福利片在线看| 国产亚洲精品久久久久5区| av国产久精品久网站免费入址| 五月开心婷婷网| 久久青草综合色| 欧美人与善性xxx| 国产片内射在线| 久热这里只有精品99| 国产亚洲欧美精品永久| 国产精品久久久av美女十八| 婷婷丁香在线五月| 在线观看一区二区三区激情| 午夜激情久久久久久久| 精品亚洲成国产av| 日本vs欧美在线观看视频| 日韩制服丝袜自拍偷拍| bbb黄色大片| 在线亚洲精品国产二区图片欧美| 国产色视频综合| 一区二区三区激情视频| 亚洲激情五月婷婷啪啪| 777久久人妻少妇嫩草av网站| 欧美xxⅹ黑人| 国产片特级美女逼逼视频| 亚洲国产精品一区二区三区在线| 在现免费观看毛片| 亚洲av日韩在线播放| 久久精品亚洲熟妇少妇任你| 交换朋友夫妻互换小说| 下体分泌物呈黄色| 亚洲av欧美aⅴ国产| 丝袜喷水一区| 免费观看a级毛片全部| 婷婷成人精品国产| 亚洲熟女精品中文字幕| 蜜桃国产av成人99| 日本色播在线视频| 欧美 日韩 精品 国产| 麻豆国产av国片精品| 国产成人欧美| 电影成人av| 欧美精品亚洲一区二区| 日韩av不卡免费在线播放| 国产黄色视频一区二区在线观看| 最新的欧美精品一区二区| 一级黄片播放器| 天天躁夜夜躁狠狠久久av| 水蜜桃什么品种好| 国产福利在线免费观看视频| 国产精品99久久99久久久不卡| 在线看a的网站| 国产片特级美女逼逼视频| 国产女主播在线喷水免费视频网站| 亚洲一区中文字幕在线| 国产精品久久久久久人妻精品电影 | 制服诱惑二区| 好男人电影高清在线观看| 久久久精品国产亚洲av高清涩受| 日本91视频免费播放| 两人在一起打扑克的视频| 精品亚洲成a人片在线观看| 性少妇av在线| 精品亚洲成a人片在线观看| 国产成人91sexporn| videos熟女内射| 日本色播在线视频| 少妇 在线观看| 国产在线免费精品| 狂野欧美激情性xxxx| 曰老女人黄片| 尾随美女入室| 啦啦啦在线观看免费高清www| 777久久人妻少妇嫩草av网站| 国产在线免费精品| www.熟女人妻精品国产| 天堂俺去俺来也www色官网| 肉色欧美久久久久久久蜜桃| 黑人欧美特级aaaaaa片| 国产在线观看jvid| 亚洲一卡2卡3卡4卡5卡精品中文| 一二三四在线观看免费中文在| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美日韩高清在线视频 | 久久久国产一区二区| 又紧又爽又黄一区二区| 成在线人永久免费视频| 人人澡人人妻人| 男人添女人高潮全过程视频| 免费在线观看黄色视频的| 97精品久久久久久久久久精品| 一级毛片黄色毛片免费观看视频| 亚洲,一卡二卡三卡| 少妇被粗大的猛进出69影院| 在线av久久热| 成年av动漫网址| 一区二区三区精品91| 老司机亚洲免费影院| 午夜福利在线免费观看网站| 91老司机精品| 久热这里只有精品99| 日韩大片免费观看网站| 国产高清视频在线播放一区 | 亚洲熟女精品中文字幕| 少妇 在线观看| 欧美老熟妇乱子伦牲交| 亚洲人成网站在线观看播放| 国产一区二区在线观看av| 777久久人妻少妇嫩草av网站| 欧美日韩综合久久久久久| 9色porny在线观看| av在线app专区| bbb黄色大片| 婷婷成人精品国产| 视频在线观看一区二区三区| 中文字幕精品免费在线观看视频| 国产91精品成人一区二区三区 | 成人手机av| 老司机亚洲免费影院| 久久鲁丝午夜福利片| 婷婷丁香在线五月| 精品人妻1区二区| 精品国产一区二区久久| 欧美在线黄色| 久久综合国产亚洲精品| 国产三级黄色录像| 热99国产精品久久久久久7| 男的添女的下面高潮视频| 韩国高清视频一区二区三区| 精品一区二区三区四区五区乱码 | 欧美日韩福利视频一区二区| 新久久久久国产一级毛片| 国产亚洲欧美精品永久| 男女床上黄色一级片免费看| 午夜91福利影院| 天天操日日干夜夜撸| 久久久久国产精品人妻一区二区| 欧美日韩福利视频一区二区| 亚洲第一av免费看| 亚洲情色 制服丝袜| 国产在线免费精品| 婷婷色综合www| 麻豆av在线久日| 国产亚洲一区二区精品| 久久精品久久久久久噜噜老黄| 亚洲图色成人| 免费久久久久久久精品成人欧美视频| 在线精品无人区一区二区三| 王馨瑶露胸无遮挡在线观看| 一级片免费观看大全| 丝袜在线中文字幕| 爱豆传媒免费全集在线观看| 亚洲第一青青草原| 男男h啪啪无遮挡| 久久99一区二区三区| 亚洲精品成人av观看孕妇| 国产色视频综合| 十八禁高潮呻吟视频| 国产亚洲午夜精品一区二区久久| 亚洲精品美女久久av网站| 一本久久精品| 免费高清在线观看视频在线观看| 国产爽快片一区二区三区| 一本色道久久久久久精品综合| 久久毛片免费看一区二区三区| 一区二区三区四区激情视频| 天堂中文最新版在线下载| 一边摸一边做爽爽视频免费| 亚洲欧洲国产日韩| 国产主播在线观看一区二区 | 精品免费久久久久久久清纯 | 又大又黄又爽视频免费| 午夜福利视频在线观看免费| 久久久国产一区二区| 精品少妇黑人巨大在线播放| 亚洲自偷自拍图片 自拍| 十八禁人妻一区二区| 亚洲专区中文字幕在线| 国产精品免费大片| 国产在线免费精品| 亚洲欧美中文字幕日韩二区| 麻豆国产av国片精品| 欧美精品高潮呻吟av久久| 丰满迷人的少妇在线观看| 国产男女超爽视频在线观看| 水蜜桃什么品种好| 久久精品久久久久久噜噜老黄| 欧美97在线视频| 日本av手机在线免费观看| 久久人妻福利社区极品人妻图片 | 日韩制服骚丝袜av| 国产精品香港三级国产av潘金莲 | 日本a在线网址| 久久精品久久久久久久性| 七月丁香在线播放| 欧美激情 高清一区二区三区| 王馨瑶露胸无遮挡在线观看| 婷婷色麻豆天堂久久| 丝瓜视频免费看黄片| 黑人猛操日本美女一级片| 国产精品九九99| a级毛片黄视频| 午夜福利,免费看| 人妻人人澡人人爽人人| 大香蕉久久成人网| tube8黄色片| 国产成人av激情在线播放| 少妇裸体淫交视频免费看高清 | 亚洲精品国产av蜜桃| 高清av免费在线| 久久亚洲国产成人精品v| 女人爽到高潮嗷嗷叫在线视频| 汤姆久久久久久久影院中文字幕| 亚洲国产精品国产精品| 在线观看国产h片| 最近中文字幕2019免费版| 亚洲成人免费电影在线观看 | 少妇被粗大的猛进出69影院| 少妇猛男粗大的猛烈进出视频| 亚洲一码二码三码区别大吗| 免费在线观看影片大全网站 | 亚洲图色成人| 成人影院久久| 国产熟女欧美一区二区| 啦啦啦啦在线视频资源| 国产男人的电影天堂91| 国产一区二区激情短视频 | 久久久欧美国产精品| 午夜91福利影院| 大香蕉久久成人网| 午夜影院在线不卡| 亚洲色图综合在线观看| 母亲3免费完整高清在线观看| 男人舔女人的私密视频| www.熟女人妻精品国产| 国产精品av久久久久免费| 亚洲精品在线美女| 大码成人一级视频| 精品人妻熟女毛片av久久网站| 中文乱码字字幕精品一区二区三区| 中文字幕人妻熟女乱码| 国产淫语在线视频| 久久青草综合色| 男女下面插进去视频免费观看| 晚上一个人看的免费电影| 国产亚洲午夜精品一区二区久久| 国产一区二区激情短视频 | 黄片小视频在线播放| 2021少妇久久久久久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美黑人精品巨大| 欧美日韩国产mv在线观看视频| 国产精品偷伦视频观看了| 巨乳人妻的诱惑在线观看| 久久 成人 亚洲| 这个男人来自地球电影免费观看| 一个人免费看片子| 黄色a级毛片大全视频| 交换朋友夫妻互换小说| 人人妻,人人澡人人爽秒播 | 成年动漫av网址| 亚洲欧美日韩高清在线视频 | 国产精品九九99| 亚洲欧美日韩另类电影网站| av片东京热男人的天堂| 中文字幕人妻丝袜一区二区| 五月天丁香电影| 日本wwww免费看| 一本一本久久a久久精品综合妖精| 中文乱码字字幕精品一区二区三区| 国产成人av教育| 国产又色又爽无遮挡免| 女人爽到高潮嗷嗷叫在线视频| 免费观看a级毛片全部| 午夜福利一区二区在线看| 丰满少妇做爰视频| 波多野结衣一区麻豆| 国产av一区二区精品久久| 欧美日韩国产mv在线观看视频| 成人黄色视频免费在线看| 国产亚洲av片在线观看秒播厂| 亚洲精品国产一区二区精华液| 国产91精品成人一区二区三区 | 日韩制服骚丝袜av| 中文字幕制服av|