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

    大口徑天然氣長輸管道狀態(tài)空間模型的適用性分析

    2020-04-01 03:02:20殷雄溫凱謝萍閃向營夏子杰宮敬
    石油科學(xué)通報 2020年1期
    關(guān)鍵詞:段長度管段傳遞函數(shù)

    殷雄,溫凱*,謝萍,閃向營,夏子杰,宮敬*

    1 中國石油大學(xué)(北京)城市油氣輸配技術(shù)北京市重點實驗室,北京 102249 2 中石油管道有限責(zé)任公司西部分公司,烏魯木齊 830013*通信作者, ydgj@cup.edu.cn, kewin1983@126.com

    1 天然氣管道流動計算模型介紹

    管道是長距離運輸天然氣的主要設(shè)備,天然氣管道系統(tǒng)是由多個壓氣站及其調(diào)節(jié)單元組成的高度非線性、強耦合、大時滯分布參數(shù)系統(tǒng)。天然氣管道內(nèi)氣體的流動可以用連續(xù)性方程、運動方程、能量方程、氣體狀態(tài)方程等非線性偏微分方程組來描述。要求出該方程組的解析解是非常困難,工程應(yīng)用中一般只求得方程組的近似解或數(shù)值解。數(shù)值解法主要是通過特征線法、隱式中心法等將偏微分方程轉(zhuǎn)化成差分方程求解管內(nèi)氣體流動狀態(tài)分布和動態(tài)變化過程。這些方法可以直接服務(wù)于天然氣管道的設(shè)計計算、流動分析等各個方面。但是這些方法的計算效率較低,并不適合狀態(tài)估計、實時優(yōu)化之類的天然氣管道在線應(yīng)用,無法滿足管道實時控制的要求。

    1984年,Kralik[1]等提出一種簡化的天然氣管道傳遞函數(shù)模型,探討了管道控制偏微分方程通過隱式中心差分空間離散而引入的誤差,研究結(jié)果可以提高復(fù)雜天然氣管道的計算效率。雖然Kralik僅僅把這種方法用于天然氣管道離線仿真,但是卻為傳遞函數(shù)的在線應(yīng)用提供了新思路。其他學(xué)者提出液體管道傳遞函數(shù)模型并在頻域進行分析,隨后傳遞函數(shù)模型被用于水管道的泄漏檢測。2006年,Reddy[2]等引用Kralik提出的天然氣管道傳遞函數(shù)模型用于分輸量發(fā)生變化的在線狀態(tài)估計和管道破裂的在線泄漏檢測。2010年,Behbahani-Nejad[3]等采用二階傳遞函數(shù)在MATLAB-Simulink中構(gòu)建了用于模擬天然氣管道中流體運動的模型。由于傳遞函數(shù)模型不適用于具有復(fù)雜拓撲結(jié)構(gòu)的天然氣管道,而狀態(tài)空間模型可以很好地解決這一問題,Behbahani-Nejad[4]等進一步通過傳遞函數(shù)得到天然氣管道的狀態(tài)空間模型。2012年,Alamian[5]等推導(dǎo)了不同邊界條件下天然氣管道傳遞函數(shù),并且得到不同邊界條件約束的狀態(tài)空間模型。2015年,K.A.Pambour[6]等針對天然氣管網(wǎng)的仿真控制提出了一種集總參數(shù)瞬態(tài)模型。2015年,R.Whalley[7]等建立了天然氣長輸管道的分布式參數(shù)模型,通過魯棒控制算法提出了一種閉環(huán)優(yōu)化調(diào)控策略。同在2015年,S.E.Mohammed[8]等提出了一種基于網(wǎng)絡(luò)分析與仿真的管道狀態(tài)評價模型。2016年,劉亞龍[9]等提出了一種基于自適應(yīng)漸進約束參量和控制目標(biāo)函數(shù)的控制模型。Wen[10]的研究表明集總參數(shù)狀態(tài)空間模型是對傳遞函數(shù)模型的改進,去掉了控制直通項,使方程的響應(yīng)過程更接近于實際管道。

    本文采用機理建模的方法,通過天然氣管道流動的基本方程推導(dǎo)建立狀態(tài)空間模型。一般仿真模型都是從管道的偏微分方程組通過數(shù)值解法求解,本文提出用狀態(tài)空間模型對管道水力進行建模分析。狀態(tài)空間模型是一種集總參數(shù)模型,在機理建模的基礎(chǔ)上將系統(tǒng)簡化為狀態(tài)量、輸入和輸出3類變量。通過綜合空間變量個數(shù),實現(xiàn)在保證計算精度的同時提高計算效率。狀態(tài)空間模型不僅可用于實時控制,還能夠良好的體現(xiàn)管道本身的非線性、大時滯等特性。對于仿真模型,時間上和空間上的劃分精細程度,通常對仿真結(jié)果有很大影響。同時由于在建模過程中進行了一定的近似和線性替換,所以需要對模型的適用范圍進行分析。為展示模型的實用性,本文選定管徑1016 mm的天然氣管道作為研究對象。

    2 天然氣管道的狀態(tài)空間模型

    狀態(tài)空間模型依據(jù)管內(nèi)氣體動力學(xué)方程進行推導(dǎo),其中有2處簡化假設(shè):(1)輸送過程中氣體溫度場變化忽略不計;(2)忽略動量方程中對流項的影響。以下推導(dǎo)以剛性管道的一元流動為基礎(chǔ),針對控制體的動力學(xué)方程,即連續(xù)性方程、運動方程和能量方程[11],加上氣體狀態(tài)方程構(gòu)成方程組(1)。

    式中:m為氣體的質(zhì)量流量,kg/s;p為氣體的壓力,Pa;A為管道流通截面面積,m2;ρ為氣體的密度,kg/m3;D為管道內(nèi)徑,m;g為重力加速度,m/s2;Δh為管道與水平面的垂直高度,m;L為管道的長度,m;R為氣體常數(shù),kJ/(kg·K);λ為管道水力摩阻系數(shù);Z為氣體壓縮因子;c為氣體的波速,m/s;x為管道位置變量,m;t為時間變量,s。

    天然氣長輸管道中氣體運行大多在阻力平方區(qū),摩阻系數(shù)λ可以通過科爾布魯克公式計算,也可以用由Hofer(Hofer,1973)[12]得到的適用于阻力平方區(qū)的近似顯式公式(2)計算;當(dāng)氣體確定,壓縮因子Z是氣體壓力和溫度的函數(shù),可通過公式(3)計算。

    式中:Re為輸氣管道的雷諾數(shù);r為管壁的絕對當(dāng)量粗糙度,m;pc氣體的臨界壓力,Pa;Tc氣體的臨界溫度,K。

    管段的傳遞函數(shù)則是在方程組(1)的基礎(chǔ)上推導(dǎo)得到,是關(guān)于上游的壓力、流量和下游的壓力、流量的頻域表達式。對管段的動力學(xué)方程線性化并進行拉普拉斯變換獲得頻域上的常微分方程組。假定管道的進口壓力和出口流量被指定為邊界條件,可獲得高階的傳遞函數(shù)表達式。但通過這種方式獲得的傳遞函數(shù)很難求得時域解析解,因此需要采用泰勒展開得到簡化的傳遞函數(shù)表達式。

    式中:

    狀態(tài)空間模型表達式如下

    3 狀態(tài)空間模型的適用性分析

    狀態(tài)空間模型是一個集總參數(shù)模型,即將一整段管道作為一個對象進行建模和分析。由于管道的分布參數(shù)特性,當(dāng)我們將管道簡化為某一點的狀態(tài)集合進行分析時,管道的長度將對模型的精度和適用性產(chǎn)生直接影響。根據(jù)狀態(tài)空間模型表達式(6)、(7),編程實現(xiàn)模型參數(shù)的計算。

    管道基礎(chǔ)參數(shù)設(shè)定為:管道管徑為1016 mm,粗糙度為0.04 mm,管道整體絕熱,進口規(guī)定壓力為10 MPa,出口規(guī)定流量為131 m3/s。管段按30、40、60、80、100、120、140、160、200 km共9種長度情況進行分析,對于不同的管道長度分別代入方程進行參數(shù)計算。

    我們首先根據(jù)上述參數(shù)在SPS中搭建對應(yīng)的仿真模型,進口規(guī)定壓力為10 MPa,出口規(guī)定流量為131 m3/s,模型如圖1。

    在上文中推導(dǎo)的狀態(tài)空間模型的基礎(chǔ)上,根據(jù)管道對應(yīng)狀態(tài)參數(shù)和狀態(tài)方程系數(shù)矩陣,計算得到描述管道系統(tǒng)特性的對應(yīng)狀態(tài)空間模型的參數(shù)。在MATLAB-Simulink中建立管道的可視化模型,包括描述管段特性的狀態(tài)空間模型模塊以及對應(yīng)連接的輸入輸出模塊。將計算得到狀態(tài)空間模型的參數(shù)矩陣輸入到對應(yīng)的狀態(tài)空間模塊中,狀態(tài)空間模型的輸入是管段的進口壓力和出口流量,輸出是進口流量和出口壓力。建立的Simulink模型如圖2所示。

    在分析中,我們對兩種輸入變化展開研究:一種是輸入的進口壓力產(chǎn)生階躍擾動(0.2 MPa)而出口流量不變,而另一種是壓力、流量都產(chǎn)生階躍擾動(壓力0.2 MPa、流量4 m3/s)。在這兩種情況下,分別修改管段為30、40、60、80、100、120、140、160、200 km共9種長度進行分析。

    圖1 SPS模型Fig. 1 SPS model

    圖2 MATLAB-Simulink模型Fig. 2 Matlab-simulink model

    3.1 輸入變化為壓力階躍

    當(dāng)輸入信號的變化只有壓力階躍時,不同管段長度下的仿真對比結(jié)果如下所述。

    從圖3-1和圖3-2中可以明顯看出由于管段長度增加導(dǎo)致的穩(wěn)定時間增加,同時模型與SPS的仿真結(jié)果之間的偏差也在增大。由仿真結(jié)果可知,當(dāng)輸入出現(xiàn)擾動時,輸出的入口流量和出口壓力值也隨之產(chǎn)生變化。對于入口流量,無論管段長度如何,其輸出流量都會在擾動出現(xiàn)的時間點瞬時增加。由于出口流量不變,輸出的流量值也會逐漸趨于穩(wěn)定值,即設(shè)定的流量值。對于出口壓力,在擾動出現(xiàn)后,輸出的壓力也會隨之變化,然后達到新的穩(wěn)態(tài)值。

    對于只有壓力變化的情況,在100 km以上長度的管道,流量經(jīng)歷短期的波動后,最后都會回到原來的穩(wěn)定值。這說明狀態(tài)空間模型在流量的關(guān)系上可以達到與SPS模型相同的結(jié)果。圖3-1(a)和圖3-2(a)可以看到狀態(tài)空間模型結(jié)果會出現(xiàn)振蕩,這是由于管道長度過短時,管道本身的壓力波傳播特性占主要因素,會表現(xiàn)出一定的波動。在實際管道中由于摩阻和氣體傳熱特性,并不會出現(xiàn)較大的波動。因此,狀態(tài)空間模型并不適用分段長度過小的情況。由圖3-1(c)和圖3-2(c)可以看出對于過大的管段分段情形,狀態(tài)空間模型的代數(shù)特性使它表現(xiàn)出與SPS模型在反應(yīng)時間上較大的偏差。圖3-1(b)和圖3-2(b)反映了狀態(tài)空間模型的仿真結(jié)果與SPS的仿真結(jié)果較為符合。其中100 km長管段條件下振蕩完全消失,且在仿真反應(yīng)時間偏差較小的情況下與SPS仿真結(jié)果較為接近。因此,對于管徑為1016的氣體管道,將其切分為100 km左右長度的管段采用狀態(tài)空間模型進行分析比較合理。

    3.2 輸入擾動為壓力和流量同時存在

    當(dāng)輸入的變化既有壓力階躍,又有流量階躍(4 m3/s)時,不同管段長度下的仿真結(jié)果如下所述。

    在壓力和流量同時出現(xiàn)階躍擾動的情況下,我們可以得到與單獨存在壓力擾動情況下相同的結(jié)論。

    由圖3-1(a)、圖3-2(a)、圖4-1(a)和圖4-2(a)中發(fā)現(xiàn)無論輸入是否存在流量階躍擾動,當(dāng)管道長度較小時,模型仿真都會存在振蕩現(xiàn)象。由3-1(b)、圖3-2(b)、圖4-1(b)和圖4-2(b)可以發(fā)現(xiàn)振蕩現(xiàn)象直到管道長度增加到80 km時才消失。理論上,隨著管段的增長,模型的精度會逐漸降低。而在本算例中,由于振蕩的存在導(dǎo)致在管段長度較小時,精度不是很高;當(dāng)振蕩消失時由于長度較大,精度也不是很高。狀態(tài)空間模型是由傳遞函數(shù)模型變形得到的,使用這類模型進行仿真時,當(dāng)管段長度足夠小時均會出現(xiàn)振蕩現(xiàn)象。隨著管段長度的增加,振蕩則會消失。管徑越大,振蕩現(xiàn)象消失時對應(yīng)的管段長度越長。

    圖3-1 只有壓力階躍時,不同管道長度下的輸出流量對比Fig. 3-1 Comparison of output flow with different pipeline length

    圖3-2 只有壓力階躍時,不同管道長度下的輸出壓力對比Fig. 3-2 Comparison of output pressure with different pipeline length

    圖4-1 即有壓力階躍,又有流量階躍時,不同管道長度下的輸出流量對比Fig. 4-1 Comparison of output flow with different pipeline length

    由圖3-1(c)、圖3-2(c)、圖4-1(c)和圖4-2(c)可以發(fā)現(xiàn),管道長度較長時,狀態(tài)空間模型仿真穩(wěn)定時間較長且與SPS仿真結(jié)果會有較大偏差。管段狀態(tài)空間模型在一階條件下其實是把管段作為存儲可壓縮氣體的圓筒看待。當(dāng)管段較短時,模型特性無法得到滿足,則會出現(xiàn)振蕩現(xiàn)象,而當(dāng)管段長度較長時,一階模型并不能良好的描述系統(tǒng),會出現(xiàn)較大的偏差。故而對于管徑為1016的天然氣管道,狀態(tài)空間模型的適用范圍為80~120 km長管段,其中以100 km為最佳。

    通過本算例我們發(fā)現(xiàn)狀態(tài)空間模型可以簡化天然氣管道系統(tǒng)的建模,在不影響仿真精度的同時,減少了建模所需的管段數(shù)量。由于能夠兼顧計算效率和精度需求,該模型適合于管道的動態(tài)控制。

    4 狀態(tài)空間模型魯棒性分析

    在由傳遞函數(shù)得到狀態(tài)空間方程的過程中,我們用穩(wěn)態(tài)線性化得到狀態(tài)空間模型用于系統(tǒng)控制特性的分析。如果模型對于穩(wěn)態(tài)點的選取過于敏感,則在此條件下做出的分析也只是針對特定的工況點而失去實際意義。因此,我們將分析不同工況點得到狀態(tài)空間模型的差異性,從而說明狀態(tài)空間模型是否能準(zhǔn)確描述管道在各個工況點的動態(tài)特性。

    在對連續(xù)性方程和運動方程進行線性化處理時,不同的穩(wěn)態(tài)可以得到不同的傳遞函數(shù)。對于長度為60 km直徑為1016 mm的天然氣管道,在穩(wěn)態(tài)流量 為55 m3/s、70 m3/s、85 m3/s、100 m3/s、115 m3/s、130 m3/s和進口壓力為7.6 MPa、8.0 MPa、8.4 MPa、8.8 MPa條件下采用Bode分析展開結(jié)果。

    線性管段模型的不確定性主要來源于一些參數(shù)值取進出口的平均近似結(jié)果以及參數(shù)由于運行點的改變而發(fā)生變化。圖5-1~5-4對不同輸入輸出組合下的狀態(tài)空間模型的魯棒性進行了分析,流量越小阻尼比越小、諧振峰值的幅值越大,進口壓力越大阻尼比越小、諧振峰值的幅值越大,諧振峰值的大小決定著時間響應(yīng)振蕩的強弱。對于1016管徑管道,建立的狀態(tài)空間模型在不同流量點線性化得到的幅頻響應(yīng)均沒有出現(xiàn)較大的偏移,模型自身的穩(wěn)定特性較好。這說明狀態(tài)空間模型在描述天然氣管道時,對于管道不同工況點具有良好的一致性。

    圖5-1 p1、p2組合傳遞函數(shù)bode分析Fig. 5-1 Bode analysis of p1 and p2 combined transfer function

    為了驗證模型的魯棒性還需要對模型在單壓力變化和單流量變化下的穩(wěn)定性進行分析,其bode圖如圖6-1和圖6-2所示。

    圖5-2 m2、p2組合傳遞函數(shù)bode分析Fig. 5-2 Bode analysis of m2 and p2 combined transfer function

    圖5-3 p1、m1組合傳遞函數(shù)bode分析Fig. 5-3 Bode analysis of p1 and m1 combined transfer function

    圖5-4 m2、m1組合傳遞函數(shù)bode分析Fig. 5-4 Bode analysis of m2 and m1 combined transfer function

    圖6-1 流量在標(biāo)稱狀態(tài)附近變化Fig. 6-1 Flow changes under standard state

    圖6-1和圖6-2為分別選取穩(wěn)態(tài)流量為100 m3/s和進口壓力8.0 MPa為標(biāo)稱狀態(tài)。圖6-1中假定穩(wěn)態(tài)流量不變,穩(wěn)態(tài)進口壓力為7.2 MPa、7.6 MPa、8.0 MPa、8.4 MPa、8.8 MPa;而 在 圖6-2中 假 定穩(wěn)態(tài)進口壓力不變,穩(wěn)態(tài)流量為70 m3/s、85 m3/s、100 m3/s、115 m3/s、130 m3/s。當(dāng)壓力、流量分別發(fā)生改變時,系統(tǒng)的頻率特性并無較大偏移,狀態(tài)空間模型的動態(tài)性能仍然非常接近。因此,我們認為狀態(tài)空間模型在描述適用長度的管段時具有較強的魯棒性。

    圖6-2 壓力在標(biāo)稱狀態(tài)附近變化Fig. 6-2 Pressure changes under standard state

    5 結(jié)論

    狀態(tài)方程作為控制理論中常用的建模方法,能夠?qū)⒐艿老到y(tǒng)簡化為狀態(tài)量、輸入和輸出3個模塊,兼顧計算精度和計算效率。選取管段長度是建模時需要考慮的重要因素,選取管段長度過短,會導(dǎo)致整體模型模塊過多,建模復(fù)雜;而長度過長時則會導(dǎo)致模型失真。本文對不同管段長度和工況點的模型準(zhǔn)確性進行了分析,結(jié)果可以分為3大類:第1類建模管段長度過短,模型振蕩嚴(yán)重;第2類長度適中,仿真結(jié)果良好;第3類長度過長,模型失真。對于只關(guān)注管道進出口壓力、流量的情況下,大口徑天然氣管道可以用狀態(tài)空間方程分段近似分析。該模型可以將管道的分段擴大到100 km左右,這對于提高計算速度和進一步綜合系統(tǒng)整體信息有很大的作用。狀態(tài)空間模型對于工況點鄰域內(nèi)的波動不是非常敏感,具有很好的動態(tài)擴展性。狀態(tài)空間模型能夠描述管道本身的非線性、大時滯等特點,在此基礎(chǔ)上可以采用經(jīng)典控制方法進行控制器的設(shè)計。

    (責(zé)任編輯 王雨墨 編輯 馬桂霞)

    猜你喜歡
    段長度管段傳遞函數(shù)
    高溫氣冷堆核電站蒸汽發(fā)生器可拆管段拆裝系統(tǒng)研究
    管段沿線流量簡化前后水頭和流行時間差異性分析
    基于LabVIEW的火焰?zhèn)鬟f函數(shù)測量系統(tǒng)
    過渡段長度對混合梁橋的受力影響
    電站配管設(shè)計中的旋轉(zhuǎn)角度分析及計算
    紅外成像系統(tǒng)的調(diào)制傳遞函數(shù)測試
    某亞臨界鍋爐水冷壁泄漏失效分析
    機械工程師(2015年8期)2015-04-10 09:52:14
    高強鋼組合K型偏心支撐框架耗能梁段長度研究
    基于傳遞函數(shù)自我優(yōu)化的BP網(wǎng)絡(luò)算法改進
    不同自由段長度預(yù)應(yīng)力錨索力學(xué)特性分析
    俺也久久电影网| 午夜免费成人在线视频| 91狼人影院| 国产精品99久久久久久久久| 免费观看精品视频网站| 在线天堂最新版资源| 国产野战对白在线观看| 国产精品爽爽va在线观看网站| 午夜精品在线福利| 在线播放国产精品三级| 免费无遮挡裸体视频| 波多野结衣高清无吗| 久久人人爽人人爽人人片va | 99热这里只有是精品50| 网址你懂的国产日韩在线| 精品人妻偷拍中文字幕| 久久久成人免费电影| 亚洲av熟女| 欧美丝袜亚洲另类 | 三级国产精品欧美在线观看| 亚洲激情在线av| 免费看美女性在线毛片视频| aaaaa片日本免费| av专区在线播放| 国内精品美女久久久久久| 国产欧美日韩一区二区三| 亚洲aⅴ乱码一区二区在线播放| 亚洲成a人片在线一区二区| 免费人成在线观看视频色| 日本黄色视频三级网站网址| 国产精品,欧美在线| 日韩欧美免费精品| 色视频www国产| 国产精品人妻久久久久久| 波多野结衣高清无吗| 久久精品综合一区二区三区| 国产成人aa在线观看| 欧美一级a爱片免费观看看| 熟女电影av网| 精品乱码久久久久久99久播| 亚洲欧美清纯卡通| 国产色爽女视频免费观看| av天堂中文字幕网| 日韩欧美一区二区三区在线观看| 精品久久久久久久末码| 在线播放国产精品三级| 九九在线视频观看精品| 丰满人妻一区二区三区视频av| 精品久久久久久久末码| 观看免费一级毛片| 成人特级黄色片久久久久久久| 黄色女人牲交| 亚洲精品久久国产高清桃花| 乱码一卡2卡4卡精品| 在线天堂最新版资源| 99热只有精品国产| 性色avwww在线观看| 亚洲第一电影网av| 欧美区成人在线视频| 国产免费av片在线观看野外av| 欧美成人性av电影在线观看| 精品久久久久久成人av| 免费av不卡在线播放| 少妇丰满av| 亚洲av日韩精品久久久久久密| 精品熟女少妇八av免费久了| 午夜视频国产福利| 人人妻人人澡欧美一区二区| 欧美日韩乱码在线| 精品久久久久久成人av| 女同久久另类99精品国产91| 波多野结衣高清作品| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久成人av| 亚洲最大成人中文| 日本 av在线| 欧美成人a在线观看| 天堂影院成人在线观看| 我的老师免费观看完整版| 国产精品一区二区三区四区久久| 别揉我奶头~嗯~啊~动态视频| 免费在线观看日本一区| 国产乱人视频| 舔av片在线| 99久久精品国产亚洲精品| 亚洲成av人片在线播放无| 婷婷精品国产亚洲av在线| 国产精品女同一区二区软件 | 99国产精品一区二区三区| 成人三级黄色视频| 日本a在线网址| 久久精品国产99精品国产亚洲性色| 成人永久免费在线观看视频| 少妇被粗大猛烈的视频| 少妇裸体淫交视频免费看高清| xxxwww97欧美| 99热这里只有精品一区| 国产久久久一区二区三区| 在线观看一区二区三区| 好男人电影高清在线观看| 搡女人真爽免费视频火全软件 | 日韩av在线大香蕉| 国模一区二区三区四区视频| 蜜桃亚洲精品一区二区三区| 国产av一区在线观看免费| 亚洲人与动物交配视频| 精品99又大又爽又粗少妇毛片 | 亚洲av电影不卡..在线观看| 欧美国产日韩亚洲一区| 欧美潮喷喷水| 亚洲,欧美,日韩| 日本成人三级电影网站| 国产精品影院久久| 日本免费a在线| 女人被狂操c到高潮| 白带黄色成豆腐渣| 亚洲综合色惰| 成年女人看的毛片在线观看| 午夜视频国产福利| 网址你懂的国产日韩在线| 99久久久亚洲精品蜜臀av| eeuss影院久久| av在线天堂中文字幕| 国产人妻一区二区三区在| 日韩欧美国产在线观看| 国产熟女xx| 欧美一区二区精品小视频在线| 亚洲18禁久久av| 精品人妻熟女av久视频| 国产免费男女视频| 久久婷婷人人爽人人干人人爱| 十八禁人妻一区二区| 国内精品久久久久精免费| 悠悠久久av| 国产一区二区在线观看日韩| 变态另类丝袜制服| 美女xxoo啪啪120秒动态图 | 亚洲精品一区av在线观看| 51午夜福利影视在线观看| 91麻豆av在线| 成人特级黄色片久久久久久久| 亚洲内射少妇av| x7x7x7水蜜桃| 欧美不卡视频在线免费观看| 亚洲国产欧洲综合997久久,| 欧美3d第一页| av福利片在线观看| 女生性感内裤真人,穿戴方法视频| 午夜精品久久久久久毛片777| 日韩欧美国产在线观看| 99热这里只有是精品50| 伦理电影大哥的女人| 成人永久免费在线观看视频| 国产成+人综合+亚洲专区| 老熟妇乱子伦视频在线观看| 精品不卡国产一区二区三区| 久久久久久九九精品二区国产| 男人舔女人下体高潮全视频| 久久人人爽人人爽人人片va | 亚洲avbb在线观看| 国产高潮美女av| 国产精品女同一区二区软件 | 网址你懂的国产日韩在线| 亚洲av美国av| 国产亚洲欧美98| 少妇熟女aⅴ在线视频| 亚洲片人在线观看| 亚洲国产色片| 成人高潮视频无遮挡免费网站| www.色视频.com| 午夜精品久久久久久毛片777| 亚洲,欧美,日韩| www.色视频.com| 国产成人aa在线观看| 人妻制服诱惑在线中文字幕| 亚洲国产欧洲综合997久久,| 婷婷丁香在线五月| 99久久久亚洲精品蜜臀av| 国产伦人伦偷精品视频| 亚洲人成网站高清观看| 丁香欧美五月| 禁无遮挡网站| 天天一区二区日本电影三级| 变态另类丝袜制服| 精品99又大又爽又粗少妇毛片 | 有码 亚洲区| 欧美黄色淫秽网站| 日本黄色视频三级网站网址| 男女之事视频高清在线观看| 久久久久久久亚洲中文字幕 | 国产久久久一区二区三区| 一区二区三区高清视频在线| 两个人视频免费观看高清| 观看免费一级毛片| 亚洲乱码一区二区免费版| 亚洲精品影视一区二区三区av| 国产精品野战在线观看| av欧美777| 欧美色视频一区免费| 欧美日韩国产亚洲二区| 欧美黄色片欧美黄色片| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久久电影| 免费看美女性在线毛片视频| 怎么达到女性高潮| 国产一区二区在线观看日韩| 三级毛片av免费| 激情在线观看视频在线高清| 国产三级黄色录像| 久久精品国产亚洲av天美| 国产高清视频在线播放一区| 免费搜索国产男女视频| 午夜视频国产福利| 国产91精品成人一区二区三区| 热99re8久久精品国产| 婷婷亚洲欧美| 久久久久久久午夜电影| 亚洲第一电影网av| av天堂在线播放| 18美女黄网站色大片免费观看| 高潮久久久久久久久久久不卡| 婷婷色综合大香蕉| 亚洲成人久久爱视频| h日本视频在线播放| 成年女人永久免费观看视频| 亚洲人成网站在线播| 91久久精品电影网| 免费搜索国产男女视频| 色5月婷婷丁香| 久久中文看片网| 麻豆av噜噜一区二区三区| 日韩中文字幕欧美一区二区| 51午夜福利影视在线观看| 国产探花在线观看一区二区| 精品久久久久久久久亚洲 | 99热精品在线国产| 丰满的人妻完整版| 如何舔出高潮| 97碰自拍视频| 在现免费观看毛片| 一个人看的www免费观看视频| 精品久久国产蜜桃| 天堂影院成人在线观看| 国产一区二区在线观看日韩| 久久99热6这里只有精品| 亚洲成人免费电影在线观看| 女人十人毛片免费观看3o分钟| 国产精品电影一区二区三区| 免费av毛片视频| 国产视频内射| 12—13女人毛片做爰片一| 久久亚洲真实| 99国产精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 国内久久婷婷六月综合欲色啪| 日韩精品青青久久久久久| www.www免费av| 色5月婷婷丁香| 免费av观看视频| 国产三级在线视频| 99精品久久久久人妻精品| 看十八女毛片水多多多| 欧美性猛交╳xxx乱大交人| 又黄又爽又刺激的免费视频.| 久久人人精品亚洲av| www.www免费av| 国产综合懂色| 无遮挡黄片免费观看| 亚洲最大成人手机在线| 免费一级毛片在线播放高清视频| eeuss影院久久| 亚洲人与动物交配视频| 免费在线观看影片大全网站| 国产精华一区二区三区| 国产白丝娇喘喷水9色精品| 国产精品女同一区二区软件 | 精品国产亚洲在线| 欧美成人一区二区免费高清观看| 在线十欧美十亚洲十日本专区| 五月玫瑰六月丁香| 一卡2卡三卡四卡精品乱码亚洲| 91午夜精品亚洲一区二区三区 | 亚洲人成网站高清观看| 97碰自拍视频| 精品人妻熟女av久视频| av在线蜜桃| 欧美黑人巨大hd| 熟妇人妻久久中文字幕3abv| 亚洲最大成人av| 日本一本二区三区精品| 色吧在线观看| 中文在线观看免费www的网站| 亚洲成人精品中文字幕电影| 一进一出抽搐动态| 人妻夜夜爽99麻豆av| 中国美女看黄片| 欧美+日韩+精品| 女生性感内裤真人,穿戴方法视频| av在线老鸭窝| 亚洲综合色惰| 丰满人妻熟妇乱又伦精品不卡| 久久人人精品亚洲av| 亚洲人与动物交配视频| 亚洲专区中文字幕在线| 国产成人aa在线观看| 不卡一级毛片| 性色avwww在线观看| 免费看光身美女| 精品乱码久久久久久99久播| 亚洲人成伊人成综合网2020| 搡老熟女国产l中国老女人| av天堂在线播放| 欧美性猛交黑人性爽| 很黄的视频免费| 成人永久免费在线观看视频| 一区福利在线观看| 国产精品久久久久久人妻精品电影| 亚洲国产精品久久男人天堂| 婷婷精品国产亚洲av| 97人妻精品一区二区三区麻豆| 国产高清有码在线观看视频| 日本在线视频免费播放| 一级黄色大片毛片| 久久精品国产清高在天天线| 亚洲av免费高清在线观看| 国产在线男女| 3wmmmm亚洲av在线观看| 久久天躁狠狠躁夜夜2o2o| 真人做人爱边吃奶动态| 88av欧美| av在线老鸭窝| 男女床上黄色一级片免费看| av在线老鸭窝| 久久中文看片网| 99久久成人亚洲精品观看| 美女免费视频网站| 最新在线观看一区二区三区| 国产真实伦视频高清在线观看 | 欧美+日韩+精品| 国产精品精品国产色婷婷| 看黄色毛片网站| 国产精品三级大全| 丁香六月欧美| 亚洲国产欧美人成| 国产激情偷乱视频一区二区| 精品久久久久久久久av| 啦啦啦韩国在线观看视频| 天堂动漫精品| 国产av麻豆久久久久久久| 美女高潮喷水抽搐中文字幕| 亚洲三级黄色毛片| 国产三级黄色录像| 免费av观看视频| 亚洲av电影在线进入| 成人美女网站在线观看视频| 亚洲内射少妇av| 日日夜夜操网爽| 波多野结衣高清作品| 精品一区二区三区av网在线观看| 深夜精品福利| 美女被艹到高潮喷水动态| 亚洲精品粉嫩美女一区| 天堂av国产一区二区熟女人妻| 久久精品影院6| 看免费av毛片| 国产一区二区在线av高清观看| 国产精品电影一区二区三区| 午夜福利免费观看在线| 免费av观看视频| 我要看日韩黄色一级片| 欧美绝顶高潮抽搐喷水| 免费搜索国产男女视频| 久久香蕉精品热| 97热精品久久久久久| 婷婷精品国产亚洲av| 国产精品久久视频播放| 午夜两性在线视频| av女优亚洲男人天堂| 身体一侧抽搐| 久99久视频精品免费| 亚洲色图av天堂| 久久久久久久久大av| 91久久精品电影网| 欧美性感艳星| www.熟女人妻精品国产| 麻豆久久精品国产亚洲av| 国产一区二区三区在线臀色熟女| 久久精品国产99精品国产亚洲性色| av天堂中文字幕网| 国产一区二区激情短视频| 亚洲不卡免费看| 日本免费一区二区三区高清不卡| 日韩av在线大香蕉| 小蜜桃在线观看免费完整版高清| 免费av毛片视频| 日本黄色片子视频| 欧美日韩黄片免| 亚洲av熟女| 麻豆一二三区av精品| 亚洲av.av天堂| 欧美3d第一页| 身体一侧抽搐| 99国产精品一区二区蜜桃av| 中文资源天堂在线| 九色成人免费人妻av| 级片在线观看| 亚洲熟妇熟女久久| 亚洲中文日韩欧美视频| 搡老妇女老女人老熟妇| av福利片在线观看| 好男人在线观看高清免费视频| ponron亚洲| 免费在线观看成人毛片| 国产真实伦视频高清在线观看 | 免费观看精品视频网站| 天堂√8在线中文| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩黄片免| 亚洲av成人不卡在线观看播放网| 一本久久中文字幕| 内射极品少妇av片p| 免费高清视频大片| 亚洲欧美精品综合久久99| 成人av一区二区三区在线看| 欧美xxxx黑人xx丫x性爽| 国产伦在线观看视频一区| 亚洲无线观看免费| 久久精品国产自在天天线| 给我免费播放毛片高清在线观看| 黄色一级大片看看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品永久免费网站| 97人妻精品一区二区三区麻豆| 国产真实乱freesex| 日本黄色视频三级网站网址| 69av精品久久久久久| 中文资源天堂在线| 九色成人免费人妻av| 午夜精品一区二区三区免费看| 九色成人免费人妻av| 国产精品一区二区免费欧美| 日韩精品青青久久久久久| 国产精品综合久久久久久久免费| 国产精品精品国产色婷婷| 精品久久久久久久久久免费视频| 婷婷六月久久综合丁香| 欧美色欧美亚洲另类二区| 国产精品综合久久久久久久免费| 夜夜爽天天搞| 91九色精品人成在线观看| 国产国拍精品亚洲av在线观看| 一卡2卡三卡四卡精品乱码亚洲| 男人舔奶头视频| 国产淫片久久久久久久久 | 好男人电影高清在线观看| 亚洲五月婷婷丁香| 网址你懂的国产日韩在线| 国产黄色小视频在线观看| 亚洲中文日韩欧美视频| 99在线人妻在线中文字幕| 老女人水多毛片| 搡女人真爽免费视频火全软件 | 久久国产乱子免费精品| 热99re8久久精品国产| 亚洲 国产 在线| 久99久视频精品免费| 精品人妻偷拍中文字幕| 精品午夜福利视频在线观看一区| 亚洲av美国av| 国产精品亚洲一级av第二区| 亚洲欧美日韩高清专用| 精品国产亚洲在线| 一本久久中文字幕| 麻豆国产av国片精品| 看免费av毛片| 国产精品av视频在线免费观看| 国产麻豆成人av免费视频| 欧美区成人在线视频| 嫩草影院精品99| 脱女人内裤的视频| 美女黄网站色视频| 欧美国产日韩亚洲一区| 男人舔女人下体高潮全视频| 免费看日本二区| 日韩欧美精品v在线| 日本撒尿小便嘘嘘汇集6| 久久热精品热| 国产淫片久久久久久久久 | 精品人妻视频免费看| 在线看三级毛片| 精品人妻1区二区| 亚洲欧美激情综合另类| 嫩草影视91久久| 狂野欧美白嫩少妇大欣赏| 九九在线视频观看精品| 亚洲五月天丁香| 国产精品乱码一区二三区的特点| 乱人视频在线观看| 天堂av国产一区二区熟女人妻| 国产三级黄色录像| 看黄色毛片网站| 欧美性感艳星| 久久久久免费精品人妻一区二区| 成人无遮挡网站| 12—13女人毛片做爰片一| 一级av片app| 亚洲av.av天堂| 欧美高清成人免费视频www| 免费搜索国产男女视频| 悠悠久久av| 日本撒尿小便嘘嘘汇集6| 国产蜜桃级精品一区二区三区| 国产av麻豆久久久久久久| 亚洲av熟女| 99久久久亚洲精品蜜臀av| 91麻豆av在线| 99在线人妻在线中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 久久性视频一级片| 日韩欧美精品v在线| 十八禁国产超污无遮挡网站| 亚洲狠狠婷婷综合久久图片| 国产精品98久久久久久宅男小说| 别揉我奶头~嗯~啊~动态视频| www日本黄色视频网| 啦啦啦观看免费观看视频高清| 日韩欧美免费精品| 午夜福利在线观看免费完整高清在 | 欧美国产日韩亚洲一区| 欧美性猛交黑人性爽| 国产精品98久久久久久宅男小说| 少妇人妻精品综合一区二区 | 女人被狂操c到高潮| 91av网一区二区| 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉精品热| 国产高清激情床上av| 美女高潮的动态| 精华霜和精华液先用哪个| 内地一区二区视频在线| 色吧在线观看| 啦啦啦韩国在线观看视频| 国产精品三级大全| 真人做人爱边吃奶动态| 九九热线精品视视频播放| 色5月婷婷丁香| 一级毛片久久久久久久久女| 伊人久久精品亚洲午夜| 国内毛片毛片毛片毛片毛片| 黄色配什么色好看| 日本 欧美在线| 国产精品亚洲一级av第二区| 午夜免费激情av| 999久久久精品免费观看国产| 可以在线观看的亚洲视频| av中文乱码字幕在线| 亚洲第一欧美日韩一区二区三区| www日本黄色视频网| www.熟女人妻精品国产| 亚洲精品色激情综合| 高潮久久久久久久久久久不卡| 亚洲av日韩精品久久久久久密| 国产伦一二天堂av在线观看| 69av精品久久久久久| 欧美区成人在线视频| 久久久国产成人免费| 精品人妻1区二区| 亚洲av电影不卡..在线观看| 欧美午夜高清在线| 人妻丰满熟妇av一区二区三区| 色综合站精品国产| 欧美日本视频| 免费人成视频x8x8入口观看| 999久久久精品免费观看国产| 人妻久久中文字幕网| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩综合久久久久久 | 亚洲欧美激情综合另类| 精品乱码久久久久久99久播| 欧美日韩国产亚洲二区| 69av精品久久久久久| 好看av亚洲va欧美ⅴa在| av中文乱码字幕在线| 极品教师在线视频| 搡女人真爽免费视频火全软件 | 久久久久久九九精品二区国产| 久久久久性生活片| 色在线成人网| 欧美3d第一页| 男女做爰动态图高潮gif福利片| 国产 一区 欧美 日韩| 变态另类丝袜制服| 中文资源天堂在线| 亚洲国产精品成人综合色| 一区二区三区激情视频| 少妇人妻精品综合一区二区 | 网址你懂的国产日韩在线| 亚洲五月天丁香| 亚洲第一区二区三区不卡| 国产精品久久视频播放| 白带黄色成豆腐渣| 真人一进一出gif抽搐免费| av中文乱码字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 一二三四社区在线视频社区8| 亚洲午夜理论影院| 男女视频在线观看网站免费| 动漫黄色视频在线观看| 免费在线观看成人毛片| 久久精品久久久久久噜噜老黄 | 无人区码免费观看不卡| 久久国产乱子免费精品| 亚洲精品乱码久久久v下载方式| 丰满人妻一区二区三区视频av| 真人做人爱边吃奶动态|