• <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)用
    国产精品成人在线| 大片免费播放器 马上看| 亚洲国产av影院在线观看| 狠狠婷婷综合久久久久久88av| 久久精品熟女亚洲av麻豆精品| 另类精品久久| 日韩免费高清中文字幕av| 亚洲精品国产一区二区精华液| 日本午夜av视频| 母亲3免费完整高清在线观看 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲精品一二三| 午夜免费观看性视频| 日韩制服丝袜自拍偷拍| 国产老妇伦熟女老妇高清| 亚洲精品av麻豆狂野| 有码 亚洲区| 中文字幕色久视频| 五月伊人婷婷丁香| 婷婷色麻豆天堂久久| 国产免费福利视频在线观看| 成年动漫av网址| 精品亚洲成a人片在线观看| 精品国产超薄肉色丝袜足j| 91精品国产国语对白视频| 国产激情久久老熟女| 性少妇av在线| 视频区图区小说| 只有这里有精品99| 久久人妻熟女aⅴ| av有码第一页| 国产探花极品一区二区| 午夜免费男女啪啪视频观看| 美女午夜性视频免费| 久久精品国产a三级三级三级| 亚洲精品在线美女| 97人妻天天添夜夜摸| 日韩欧美精品免费久久| 国精品久久久久久国模美| 久久久久人妻精品一区果冻| 国产一区亚洲一区在线观看| 五月伊人婷婷丁香| 超色免费av| 亚洲av欧美aⅴ国产| 肉色欧美久久久久久久蜜桃| 欧美最新免费一区二区三区| 欧美激情 高清一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产精品麻豆人妻色哟哟久久| 熟妇人妻不卡中文字幕| 妹子高潮喷水视频| 少妇 在线观看| 亚洲四区av| 国产精品久久久久成人av| 热re99久久国产66热| 伦理电影免费视频| 亚洲av福利一区| 日韩制服骚丝袜av| 亚洲av成人精品一二三区| 街头女战士在线观看网站| 一级,二级,三级黄色视频| 国产男女超爽视频在线观看| 夫妻性生交免费视频一级片| 欧美日韩成人在线一区二区| 精品国产露脸久久av麻豆| 久久久久国产网址| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 18+在线观看网站| 一级爰片在线观看| 欧美日韩亚洲高清精品| 丰满饥渴人妻一区二区三| 九色亚洲精品在线播放| 国产成人aa在线观看| 国产精品久久久久久精品电影小说| 亚洲精品久久午夜乱码| 亚洲内射少妇av| 日韩中文字幕视频在线看片| 七月丁香在线播放| 欧美日韩一级在线毛片| 制服诱惑二区| 免费高清在线观看视频在线观看| 亚洲一区二区三区欧美精品| 国产成人精品久久久久久| 叶爱在线成人免费视频播放| 国产精品一区二区在线观看99| 中文字幕人妻丝袜一区二区 | 99久久人妻综合| 日韩制服丝袜自拍偷拍| 中文字幕制服av| 如何舔出高潮| av线在线观看网站| 久久精品国产a三级三级三级| 亚洲精品一区蜜桃| 18禁观看日本| 久久精品国产鲁丝片午夜精品| 在线亚洲精品国产二区图片欧美| 亚洲国产成人一精品久久久| 高清在线视频一区二区三区| 精品亚洲成国产av| av在线app专区| 在线观看一区二区三区激情| 老女人水多毛片| 波多野结衣av一区二区av| 新久久久久国产一级毛片| 国产xxxxx性猛交| 亚洲精品国产av蜜桃| 婷婷成人精品国产| 亚洲美女搞黄在线观看| 精品亚洲成a人片在线观看| 国产色婷婷99| 国产色婷婷99| 国产精品国产av在线观看| 尾随美女入室| 精品少妇黑人巨大在线播放| 久久鲁丝午夜福利片| 女人被躁到高潮嗷嗷叫费观| 免费女性裸体啪啪无遮挡网站| 99国产综合亚洲精品| 久久久精品免费免费高清| 国产xxxxx性猛交| 久久久国产欧美日韩av| 18禁动态无遮挡网站| 18禁动态无遮挡网站| 精品99又大又爽又粗少妇毛片| 欧美成人精品欧美一级黄| 国产精品久久久久久久久免| 自线自在国产av| 成年av动漫网址| 国产亚洲精品第一综合不卡| 精品亚洲成国产av| 男女无遮挡免费网站观看| av视频免费观看在线观看| 欧美激情 高清一区二区三区| 中文字幕人妻丝袜一区二区 | 大香蕉久久网| 最近中文字幕2019免费版| 午夜福利在线观看免费完整高清在| 天堂中文最新版在线下载| 国产熟女欧美一区二区| 黄频高清免费视频| 日本-黄色视频高清免费观看| 纵有疾风起免费观看全集完整版| 久久99一区二区三区| 电影成人av| 亚洲av中文av极速乱| 汤姆久久久久久久影院中文字幕| 国产黄频视频在线观看| 在线观看免费日韩欧美大片| 如日韩欧美国产精品一区二区三区| 少妇被粗大猛烈的视频| 精品一区二区免费观看| 大码成人一级视频| 伦理电影免费视频| 国产日韩一区二区三区精品不卡| 亚洲色图综合在线观看| 中国国产av一级| 国产成人精品婷婷| 久久精品国产自在天天线| 亚洲美女视频黄频| av福利片在线| 成人国语在线视频| 有码 亚洲区| 综合色丁香网| 国语对白做爰xxxⅹ性视频网站| 婷婷色综合www| 最近2019中文字幕mv第一页| 国产成人欧美| 精品亚洲成国产av| 国产成人午夜福利电影在线观看| 不卡视频在线观看欧美| 色网站视频免费| 老鸭窝网址在线观看| 成人黄色视频免费在线看| 亚洲av.av天堂| 青春草亚洲视频在线观看| 精品卡一卡二卡四卡免费| 下体分泌物呈黄色| 黄网站色视频无遮挡免费观看| 九九爱精品视频在线观看| 久久精品熟女亚洲av麻豆精品| 91精品三级在线观看| 99国产综合亚洲精品| 欧美精品人与动牲交sv欧美| 亚洲美女黄色视频免费看| 99热网站在线观看| 亚洲欧洲精品一区二区精品久久久 | 久久人人爽人人片av| 免费黄网站久久成人精品| 久久毛片免费看一区二区三区| 伦理电影大哥的女人| 最新的欧美精品一区二区| 精品久久久久久电影网| 91久久精品国产一区二区三区| 在线观看免费高清a一片| 在现免费观看毛片| 麻豆乱淫一区二区| www.自偷自拍.com| h视频一区二区三区| 久久久久国产网址| 制服丝袜香蕉在线| 国产熟女午夜一区二区三区| 一区二区av电影网| 人人澡人人妻人| 久久韩国三级中文字幕| 春色校园在线视频观看| 亚洲五月色婷婷综合| 男人舔女人的私密视频| 亚洲欧美色中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 国产精品久久久久成人av| 丰满饥渴人妻一区二区三| 久久精品国产亚洲av天美| 欧美激情 高清一区二区三区| 亚洲成人av在线免费| 欧美亚洲 丝袜 人妻 在线| 一本色道久久久久久精品综合| 91精品伊人久久大香线蕉| av不卡在线播放| 另类精品久久| 国产精品香港三级国产av潘金莲 | 免费av中文字幕在线| 午夜福利视频精品| 丰满饥渴人妻一区二区三| 国产又爽黄色视频| 国产一区二区激情短视频 | 国产成人精品久久二区二区91 | 久久精品国产亚洲av涩爱| 精品国产露脸久久av麻豆| 久久精品国产自在天天线| 日本爱情动作片www.在线观看| 在线免费观看不下载黄p国产| 狠狠精品人妻久久久久久综合| 国产精品香港三级国产av潘金莲 | av国产久精品久网站免费入址| 国产精品国产三级国产专区5o| 亚洲国产看品久久| 大陆偷拍与自拍| 两性夫妻黄色片| 国产高清国产精品国产三级| 欧美日韩国产mv在线观看视频| 免费看av在线观看网站| 国产av一区二区精品久久| av.在线天堂| 久久这里只有精品19| 青青草视频在线视频观看| 丝袜脚勾引网站| 在线观看人妻少妇| 日韩,欧美,国产一区二区三区| 亚洲av欧美aⅴ国产| xxx大片免费视频| 高清不卡的av网站| 丝袜脚勾引网站| 国产xxxxx性猛交| 在线天堂最新版资源| 黄片无遮挡物在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲,一卡二卡三卡| 天天躁日日躁夜夜躁夜夜| 99九九在线精品视频| 亚洲精品第二区| 久久精品国产鲁丝片午夜精品| 老汉色av国产亚洲站长工具| 熟女电影av网| 一区在线观看完整版| 成人国语在线视频| 一二三四中文在线观看免费高清| 亚洲av免费高清在线观看| 中文字幕亚洲精品专区| www.av在线官网国产| 黄片无遮挡物在线观看| 国产成人91sexporn| 美女xxoo啪啪120秒动态图| 超碰成人久久| 精品国产超薄肉色丝袜足j| 亚洲av在线观看美女高潮| 赤兔流量卡办理| 精品酒店卫生间| 欧美精品人与动牲交sv欧美| 亚洲欧美成人综合另类久久久| 搡老乐熟女国产| 国产老妇伦熟女老妇高清| 成年av动漫网址| 午夜福利一区二区在线看| 在线天堂最新版资源| 午夜av观看不卡| 美女国产视频在线观看| 丰满乱子伦码专区| 婷婷色综合www| 成人二区视频| av在线老鸭窝| 99国产综合亚洲精品| 电影成人av| 亚洲综合色网址| 亚洲国产色片| 国产免费又黄又爽又色| 国产av精品麻豆| 亚洲精品,欧美精品| 一边亲一边摸免费视频| 好男人视频免费观看在线| 麻豆乱淫一区二区| 精品一品国产午夜福利视频| 国产精品一区二区在线观看99| 黄片无遮挡物在线观看| 可以免费在线观看a视频的电影网站 | 国产精品久久久久久久久免| 免费黄频网站在线观看国产| 国产成人精品在线电影| 国产黄色免费在线视频| 国产精品人妻久久久影院| 久久精品国产综合久久久| av不卡在线播放| 三上悠亚av全集在线观看| 久久久久久久亚洲中文字幕| 午夜福利在线免费观看网站| 丝袜人妻中文字幕| av网站在线播放免费| 嫩草影院入口| 成人二区视频| 国产色婷婷99| 你懂的网址亚洲精品在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲精品国产一区二区精华液| 最近手机中文字幕大全| 在线观看www视频免费| 精品一品国产午夜福利视频| 日本-黄色视频高清免费观看| 欧美精品国产亚洲| 精品国产国语对白av| 晚上一个人看的免费电影| 高清视频免费观看一区二区| av在线观看视频网站免费| 成年女人毛片免费观看观看9 | 飞空精品影院首页| 国产色婷婷99| 亚洲人成77777在线视频| 国产熟女午夜一区二区三区| 综合色丁香网| 人妻系列 视频| 国产淫语在线视频| 女人久久www免费人成看片| 看免费成人av毛片| 人人妻人人爽人人添夜夜欢视频| 亚洲一区二区三区欧美精品| 又粗又硬又长又爽又黄的视频| 咕卡用的链子| 在线亚洲精品国产二区图片欧美| 两个人免费观看高清视频| 国产白丝娇喘喷水9色精品| 亚洲精品久久久久久婷婷小说| 男女边吃奶边做爰视频| 卡戴珊不雅视频在线播放| 在线观看www视频免费| 日韩欧美一区视频在线观看| 亚洲国产精品一区二区三区在线| 中文精品一卡2卡3卡4更新| 亚洲四区av| 久久久a久久爽久久v久久| 午夜福利影视在线免费观看| 韩国精品一区二区三区| 亚洲精品美女久久久久99蜜臀 | 日韩成人av中文字幕在线观看| 色播在线永久视频| 日韩中文字幕欧美一区二区 | 欧美老熟妇乱子伦牲交| 日韩一本色道免费dvd| 亚洲精品国产色婷婷电影| 亚洲av免费高清在线观看| 欧美成人精品欧美一级黄| 免费在线观看完整版高清| 视频在线观看一区二区三区| 美女xxoo啪啪120秒动态图| 午夜福利视频精品| 狂野欧美激情性bbbbbb| 免费少妇av软件| 久久久国产精品麻豆| 婷婷色综合大香蕉| 一区二区三区四区激情视频| 中文字幕制服av| 免费观看av网站的网址| 精品国产乱码久久久久久小说| 天堂8中文在线网| 老汉色∧v一级毛片| 人妻一区二区av| 如日韩欧美国产精品一区二区三区| 国产综合精华液| 少妇熟女欧美另类| 亚洲,欧美,日韩| 男人操女人黄网站| 青春草亚洲视频在线观看| 国产无遮挡羞羞视频在线观看| 人人妻人人爽人人添夜夜欢视频| 国产精品无大码| www.熟女人妻精品国产| 一级黄片播放器| 18+在线观看网站| 最新的欧美精品一区二区| 中国国产av一级| 妹子高潮喷水视频| 精品少妇黑人巨大在线播放| 中文欧美无线码| 黄色怎么调成土黄色| 国产精品久久久久久精品古装| 精品人妻熟女毛片av久久网站| 国产国语露脸激情在线看| 精品国产国语对白av| 久久婷婷青草| 欧美bdsm另类| 久久青草综合色| 久久精品久久久久久久性| 国产一区二区激情短视频 | 亚洲国产精品成人久久小说| 老女人水多毛片| 男女啪啪激烈高潮av片| 国产精品一区二区在线不卡| 婷婷色av中文字幕| 熟女电影av网| 少妇人妻久久综合中文| 婷婷色综合www| 宅男免费午夜| 香蕉国产在线看| 色播在线永久视频| 一区福利在线观看| a级毛片在线看网站| 男女国产视频网站| 高清欧美精品videossex| 国产亚洲精品第一综合不卡| 人成视频在线观看免费观看| 午夜激情久久久久久久| av女优亚洲男人天堂| 人妻 亚洲 视频| 91国产中文字幕| xxxhd国产人妻xxx| 国产亚洲av片在线观看秒播厂| 亚洲欧美精品自产自拍| 搡老乐熟女国产| 三上悠亚av全集在线观看| 国产在线免费精品| 日韩一卡2卡3卡4卡2021年| 人妻少妇偷人精品九色| 国产毛片在线视频| 国产精品一国产av| 在现免费观看毛片| 国产精品不卡视频一区二区| 久久久久精品久久久久真实原创| 在线 av 中文字幕| 精品人妻在线不人妻| 老司机影院毛片| 搡老乐熟女国产| 精品99又大又爽又粗少妇毛片| 国产有黄有色有爽视频| 在线观看国产h片| 韩国高清视频一区二区三区| 美女脱内裤让男人舔精品视频| 国产精品三级大全| 女的被弄到高潮叫床怎么办| 国产av精品麻豆| 国产高清国产精品国产三级| 亚洲国产成人一精品久久久| 80岁老熟妇乱子伦牲交| 2022亚洲国产成人精品| 免费观看性生交大片5| 老熟女久久久| 成人毛片60女人毛片免费| 好男人视频免费观看在线| 一区福利在线观看| 少妇被粗大猛烈的视频| 欧美精品亚洲一区二区| 91精品三级在线观看| 大码成人一级视频| 亚洲精品国产av蜜桃| 寂寞人妻少妇视频99o| 久久久久人妻精品一区果冻| 亚洲欧美色中文字幕在线| 大香蕉久久成人网| 久久精品国产综合久久久| 久久av网站| 久久久欧美国产精品| 男女边摸边吃奶| 青青草视频在线视频观看| 国产日韩欧美亚洲二区| 最新的欧美精品一区二区| 亚洲av电影在线观看一区二区三区| 精品国产超薄肉色丝袜足j| 午夜免费观看性视频| 超色免费av| 日日摸夜夜添夜夜爱| 久久久a久久爽久久v久久| 日韩av在线免费看完整版不卡| av国产精品久久久久影院| 天天躁狠狠躁夜夜躁狠狠躁| 视频在线观看一区二区三区| 卡戴珊不雅视频在线播放| 亚洲第一av免费看| 街头女战士在线观看网站| 亚洲欧美一区二区三区久久| 国产毛片在线视频| 国产极品粉嫩免费观看在线| 国产黄色视频一区二区在线观看| 国产精品国产av在线观看| 波多野结衣一区麻豆| 婷婷色综合www| 日韩一区二区视频免费看| 亚洲色图 男人天堂 中文字幕| 久久国产亚洲av麻豆专区| 在线观看人妻少妇| 国产精品偷伦视频观看了| av卡一久久| 少妇的逼水好多| 七月丁香在线播放| 欧美成人午夜精品| 国产一区二区三区综合在线观看| 国产欧美日韩综合在线一区二区| 丝袜脚勾引网站| 看十八女毛片水多多多| 日韩三级伦理在线观看| 精品国产超薄肉色丝袜足j| 肉色欧美久久久久久久蜜桃| 亚洲精品av麻豆狂野| 中文字幕人妻丝袜一区二区 | 大陆偷拍与自拍| 国产精品三级大全| 乱人伦中国视频| 国产精品 国内视频| 欧美精品av麻豆av| 亚洲图色成人| 日韩av不卡免费在线播放| 午夜福利视频精品| 国产亚洲精品第一综合不卡| 99精国产麻豆久久婷婷| 亚洲av成人精品一二三区| 国产成人精品婷婷| 黄色怎么调成土黄色| 性色avwww在线观看| 一边摸一边做爽爽视频免费| 亚洲美女搞黄在线观看| 精品人妻在线不人妻| 看非洲黑人一级黄片| 国产精品久久久久久精品电影小说| 老汉色av国产亚洲站长工具| av网站在线播放免费| 性色av一级| 亚洲一码二码三码区别大吗| 国产一区亚洲一区在线观看| 欧美成人午夜精品| 人体艺术视频欧美日本| 多毛熟女@视频| 国产深夜福利视频在线观看| 观看av在线不卡| 国产精品熟女久久久久浪| 国产熟女欧美一区二区| xxxhd国产人妻xxx| 1024香蕉在线观看| 午夜激情久久久久久久| 国产野战对白在线观看| 只有这里有精品99| 欧美日韩av久久| 免费观看av网站的网址| 天堂8中文在线网| 午夜福利一区二区在线看| 高清av免费在线| 99久久人妻综合| 香蕉国产在线看| 亚洲国产欧美日韩在线播放| 老鸭窝网址在线观看| 国产精品久久久久久av不卡| 久久久精品免费免费高清| 色婷婷av一区二区三区视频| 最近2019中文字幕mv第一页| 亚洲精品美女久久av网站| 久久影院123| 波多野结衣av一区二区av| 亚洲av成人精品一二三区| 欧美亚洲日本最大视频资源| 国产一区二区三区av在线| 777久久人妻少妇嫩草av网站| 久久女婷五月综合色啪小说| 婷婷色av中文字幕| 亚洲综合色惰| 亚洲国产欧美日韩在线播放| 国产 一区精品| 亚洲国产成人一精品久久久| 日韩中字成人| 国产日韩欧美视频二区| 中文欧美无线码| 丝袜在线中文字幕| 国产黄频视频在线观看| 久久精品久久久久久久性| 一区二区三区精品91| 麻豆乱淫一区二区| 亚洲av中文av极速乱| 天堂俺去俺来也www色官网| 一级片免费观看大全| 精品久久久精品久久久| 99精国产麻豆久久婷婷| 我要看黄色一级片免费的| 成人毛片a级毛片在线播放| 下体分泌物呈黄色| 超色免费av| 久久久久久人人人人人| 亚洲国产av新网站| 老女人水多毛片| 高清av免费在线| 国产精品不卡视频一区二区| 国产日韩一区二区三区精品不卡| 亚洲国产精品999| 卡戴珊不雅视频在线播放| 99精国产麻豆久久婷婷| 国产在视频线精品| 成人毛片a级毛片在线播放| 国产精品国产三级国产专区5o| 如日韩欧美国产精品一区二区三区| 亚洲精品国产av蜜桃| 亚洲国产av新网站| 午夜福利视频在线观看免费| 一级毛片 在线播放|