• <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é)特性分析
    久久免费观看电影| 丝瓜视频免费看黄片| 午夜激情av网站| 女性被躁到高潮视频| 水蜜桃什么品种好| 午夜免费鲁丝| 黑人猛操日本美女一级片| 亚洲精品色激情综合| 一区在线观看完整版| 街头女战士在线观看网站| 美女大奶头黄色视频| 成人免费观看视频高清| 久久精品国产鲁丝片午夜精品| 大香蕉久久网| 成人综合一区亚洲| 久久久久久人妻| 麻豆精品久久久久久蜜桃| 精品久久久精品久久久| 春色校园在线视频观看| 亚洲高清免费不卡视频| 高清欧美精品videossex| 免费黄网站久久成人精品| 99re6热这里在线精品视频| 汤姆久久久久久久影院中文字幕| 视频中文字幕在线观看| 亚洲av综合色区一区| 久久av网站| 汤姆久久久久久久影院中文字幕| 久久久久久久精品精品| 自线自在国产av| 看免费成人av毛片| 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦视频在线资源免费观看| 国产无遮挡羞羞视频在线观看| 赤兔流量卡办理| 一级,二级,三级黄色视频| 亚洲内射少妇av| 国产日韩一区二区三区精品不卡 | 最近的中文字幕免费完整| 亚洲国产色片| 国产精品久久久久久精品电影小说| 国产在视频线精品| 99九九在线精品视频| 欧美xxxx性猛交bbbb| 一级毛片aaaaaa免费看小| 亚洲欧美精品自产自拍| 色哟哟·www| 亚洲国产色片| 日日摸夜夜添夜夜添av毛片| 五月伊人婷婷丁香| av国产久精品久网站免费入址| 丝瓜视频免费看黄片| 亚洲av二区三区四区| 草草在线视频免费看| 91精品国产国语对白视频| 18禁动态无遮挡网站| 岛国毛片在线播放| 一区二区av电影网| 亚洲精品国产av成人精品| 老女人水多毛片| .国产精品久久| 一区二区av电影网| 久久影院123| 亚洲国产精品一区三区| 精品一区二区三卡| 在线观看免费视频网站a站| 赤兔流量卡办理| 日本免费在线观看一区| 国产黄色免费在线视频| 免费观看a级毛片全部| 观看av在线不卡| 欧美bdsm另类| 免费不卡的大黄色大毛片视频在线观看| 国产精品免费大片| 国国产精品蜜臀av免费| 热re99久久精品国产66热6| 国产精品99久久久久久久久| 亚洲国产av新网站| 五月天丁香电影| 国产色婷婷99| 免费人妻精品一区二区三区视频| 午夜福利视频精品| 91成人精品电影| 一级毛片黄色毛片免费观看视频| a级毛片免费高清观看在线播放| 99热全是精品| 免费看光身美女| 精品久久久久久久久亚洲| 啦啦啦视频在线资源免费观看| 丝瓜视频免费看黄片| 美女脱内裤让男人舔精品视频| 一级毛片 在线播放| 久久国内精品自在自线图片| 亚洲av成人精品一区久久| 欧美最新免费一区二区三区| 人人妻人人澡人人爽人人夜夜| 51国产日韩欧美| 少妇熟女欧美另类| tube8黄色片| 一个人免费看片子| 欧美激情 高清一区二区三区| 全区人妻精品视频| 欧美日韩在线观看h| 久久久久视频综合| 啦啦啦视频在线资源免费观看| 欧美xxxx性猛交bbbb| 国产色婷婷99| 熟女av电影| 久久国产精品大桥未久av| 人人妻人人爽人人添夜夜欢视频| 99精国产麻豆久久婷婷| 日韩人妻高清精品专区| 搡女人真爽免费视频火全软件| 久久精品久久精品一区二区三区| 午夜福利视频在线观看免费| av在线观看视频网站免费| 日本av手机在线免费观看| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 狂野欧美激情性xxxx在线观看| 成人二区视频| 欧美日韩综合久久久久久| 亚洲综合色惰| 免费久久久久久久精品成人欧美视频 | 国产精品麻豆人妻色哟哟久久| 伦精品一区二区三区| 午夜福利在线观看免费完整高清在| 亚洲av国产av综合av卡| 最新中文字幕久久久久| 日本vs欧美在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 2018国产大陆天天弄谢| 免费观看无遮挡的男女| 大又大粗又爽又黄少妇毛片口| 国产黄色免费在线视频| 最后的刺客免费高清国语| 国产高清三级在线| 日本与韩国留学比较| 男男h啪啪无遮挡| 一个人免费看片子| 制服诱惑二区| 亚洲av免费高清在线观看| 美女中出高潮动态图| 2022亚洲国产成人精品| 91精品伊人久久大香线蕉| 丝袜喷水一区| 如日韩欧美国产精品一区二区三区 | 欧美xxxx性猛交bbbb| 99re6热这里在线精品视频| 国产精品秋霞免费鲁丝片| 国产精品嫩草影院av在线观看| 精品视频人人做人人爽| 国产探花极品一区二区| 久久久久网色| 99热6这里只有精品| 极品少妇高潮喷水抽搐| 国产成人午夜福利电影在线观看| 一级毛片黄色毛片免费观看视频| 黄片无遮挡物在线观看| 2018国产大陆天天弄谢| 久久韩国三级中文字幕| 国产精品一区www在线观看| 另类精品久久| 日本91视频免费播放| 国产精品99久久久久久久久| 七月丁香在线播放| 亚洲精品成人av观看孕妇| 亚洲中文av在线| 国产欧美另类精品又又久久亚洲欧美| 韩国高清视频一区二区三区| 三级国产精品欧美在线观看| 日日摸夜夜添夜夜爱| 一级爰片在线观看| 日本午夜av视频| 插阴视频在线观看视频| 精品人妻熟女毛片av久久网站| 国产成人91sexporn| 国产精品久久久久久精品电影小说| 国产精品99久久久久久久久| 有码 亚洲区| 18禁动态无遮挡网站| 亚洲欧美精品自产自拍| 人人澡人人妻人| av在线app专区| 亚洲精品日韩在线中文字幕| 午夜福利影视在线免费观看| 国产一区二区三区综合在线观看 | 日韩大片免费观看网站| 国产免费一区二区三区四区乱码| 久久精品国产a三级三级三级| 人人妻人人澡人人爽人人夜夜| 99久久中文字幕三级久久日本| 成人无遮挡网站| 亚洲精品国产av蜜桃| 亚洲经典国产精华液单| 成人国产麻豆网| 精品久久久精品久久久| 午夜影院在线不卡| 亚洲av欧美aⅴ国产| 精品一区在线观看国产| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放| 日产精品乱码卡一卡2卡三| 搡老乐熟女国产| 青春草国产在线视频| 亚洲,欧美,日韩| 婷婷色麻豆天堂久久| 精品午夜福利在线看| 久久97久久精品| 久久精品久久久久久久性| 一区二区日韩欧美中文字幕 | 国产高清不卡午夜福利| 国产日韩一区二区三区精品不卡 | 国产极品粉嫩免费观看在线 | 插逼视频在线观看| 国产日韩一区二区三区精品不卡 | 日韩精品有码人妻一区| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 成人漫画全彩无遮挡| av在线app专区| 一区二区三区精品91| 一区二区三区乱码不卡18| 日韩成人伦理影院| 国产亚洲午夜精品一区二区久久| 高清视频免费观看一区二区| 一边亲一边摸免费视频| 午夜精品国产一区二区电影| 亚洲激情五月婷婷啪啪| 美女xxoo啪啪120秒动态图| 日日爽夜夜爽网站| 国产免费视频播放在线视频| 在线免费观看不下载黄p国产| 高清在线视频一区二区三区| 亚洲经典国产精华液单| 国产精品无大码| 日本猛色少妇xxxxx猛交久久| 大香蕉97超碰在线| 国产日韩欧美亚洲二区| 肉色欧美久久久久久久蜜桃| 一级片'在线观看视频| 少妇 在线观看| 亚洲久久久国产精品| 亚洲国产av影院在线观看| 夜夜骑夜夜射夜夜干| 国产色婷婷99| 男女免费视频国产| 99久久中文字幕三级久久日本| 成人无遮挡网站| 久久精品国产亚洲网站| 精品一区二区三卡| 成年美女黄网站色视频大全免费 | 国产亚洲一区二区精品| www.av在线官网国产| 免费人妻精品一区二区三区视频| 边亲边吃奶的免费视频| 这个男人来自地球电影免费观看 | 成人影院久久| 在线免费观看不下载黄p国产| 国产探花极品一区二区| 草草在线视频免费看| 一级毛片我不卡| 国产黄色视频一区二区在线观看| 91精品一卡2卡3卡4卡| 久久久久久久久久久免费av| 观看av在线不卡| 老司机影院成人| 亚洲在久久综合| 久久久久久伊人网av| 亚洲国产色片| 久久精品国产亚洲av涩爱| 亚洲熟女精品中文字幕| 免费大片18禁| 亚洲国产成人一精品久久久| 人妻制服诱惑在线中文字幕| 免费高清在线观看视频在线观看| 搡老乐熟女国产| 韩国av在线不卡| 精品人妻熟女毛片av久久网站| 亚洲av成人精品一二三区| 下体分泌物呈黄色| 人人妻人人爽人人添夜夜欢视频| 久久久亚洲精品成人影院| 国产一区二区三区综合在线观看 | 91精品一卡2卡3卡4卡| 91久久精品电影网| 免费久久久久久久精品成人欧美视频 | 亚洲怡红院男人天堂| 超碰97精品在线观看| 久久久亚洲精品成人影院| 国产精品久久久久成人av| 一个人看视频在线观看www免费| 精品久久久久久久久av| 夜夜爽夜夜爽视频| 亚洲av福利一区| 亚洲精品一区蜜桃| 亚洲伊人久久精品综合| 18禁动态无遮挡网站| 在线观看免费日韩欧美大片 | 国产熟女午夜一区二区三区 | 亚洲欧美成人综合另类久久久| 天堂俺去俺来也www色官网| 成人毛片60女人毛片免费| 国语对白做爰xxxⅹ性视频网站| 久久久久国产网址| 91在线精品国自产拍蜜月| 考比视频在线观看| 丰满乱子伦码专区| 免费看不卡的av| 九色成人免费人妻av| 少妇精品久久久久久久| 十分钟在线观看高清视频www| 国产精品久久久久久久久免| 国产熟女午夜一区二区三区 | 国产成人91sexporn| 久久久久久久大尺度免费视频| 亚洲精品亚洲一区二区| 精品人妻熟女av久视频| 久久热精品热| 精品人妻熟女av久视频| 亚洲激情五月婷婷啪啪| 99九九在线精品视频| 国产乱人偷精品视频| 老女人水多毛片| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 26uuu在线亚洲综合色| 国产免费视频播放在线视频| 高清毛片免费看| 日本-黄色视频高清免费观看| 精品少妇内射三级| 久久精品人人爽人人爽视色| 亚洲成人手机| 成人亚洲精品一区在线观看| 亚洲欧美日韩卡通动漫| 日韩熟女老妇一区二区性免费视频| 国产黄片视频在线免费观看| 国产精品一区二区三区四区免费观看| 在线观看免费高清a一片| 亚洲精品美女久久av网站| 亚洲国产精品一区三区| 欧美+日韩+精品| 制服人妻中文乱码| 久久久久久久久久成人| 91在线精品国自产拍蜜月| 99re6热这里在线精品视频| 一级二级三级毛片免费看| kizo精华| 丰满迷人的少妇在线观看| 亚洲精华国产精华液的使用体验| 中国国产av一级| 考比视频在线观看| 伦理电影免费视频| 熟女av电影| 久久影院123| 日韩在线高清观看一区二区三区| 国产高清国产精品国产三级| 午夜福利,免费看| 欧美性感艳星| 日韩,欧美,国产一区二区三区| 久久久久久久精品精品| 国产探花极品一区二区| 在线观看免费高清a一片| 成人免费观看视频高清| 一级毛片我不卡| 国产老妇伦熟女老妇高清| 久久鲁丝午夜福利片| 午夜免费鲁丝| 免费看光身美女| 高清不卡的av网站| 久久久久久久久久人人人人人人| 丰满少妇做爰视频| 寂寞人妻少妇视频99o| 国产精品熟女久久久久浪| 亚洲国产精品一区二区三区在线| 人妻少妇偷人精品九色| 国产黄色免费在线视频| 999精品在线视频| 成人黄色视频免费在线看| 亚洲国产欧美日韩在线播放| 久久久久久久久久久久大奶| 丝袜美足系列| 久久久久久久国产电影| 国产永久视频网站| 美女内射精品一级片tv| 免费看不卡的av| 国产成人91sexporn| 婷婷色av中文字幕| 人妻夜夜爽99麻豆av| 美女福利国产在线| 人妻 亚洲 视频| 美女国产视频在线观看| 在线观看www视频免费| 亚洲成色77777| 欧美性感艳星| 欧美一级a爱片免费观看看| 国产熟女午夜一区二区三区 | 男女边摸边吃奶| 亚洲精品成人av观看孕妇| 只有这里有精品99| 香蕉精品网在线| 视频区图区小说| 日日啪夜夜爽| 亚洲综合色网址| 精品人妻在线不人妻| 亚洲欧洲精品一区二区精品久久久 | 日韩成人伦理影院| 日本午夜av视频| 日韩在线高清观看一区二区三区| 免费人妻精品一区二区三区视频| 亚洲欧美日韩卡通动漫| 国产一级毛片在线| 大陆偷拍与自拍| 色视频在线一区二区三区| 久久婷婷青草| 亚洲美女视频黄频| 51国产日韩欧美| 人人妻人人澡人人爽人人夜夜| 亚洲av电影在线观看一区二区三区| videos熟女内射| 18在线观看网站| 欧美 亚洲 国产 日韩一| 国产欧美亚洲国产| 免费av中文字幕在线| 久久久精品免费免费高清| 青春草亚洲视频在线观看| 日本黄色日本黄色录像| 一区二区三区乱码不卡18| 高清黄色对白视频在线免费看| 亚洲欧美一区二区三区国产| 欧美日韩在线观看h| 免费黄网站久久成人精品| 成人国产av品久久久| 99热这里只有精品一区| 亚洲精品美女久久av网站| 91精品一卡2卡3卡4卡| 国产一区二区在线观看日韩| 国产精品嫩草影院av在线观看| 色5月婷婷丁香| 日日撸夜夜添| 美女国产高潮福利片在线看| 不卡视频在线观看欧美| 精品午夜福利在线看| 久久久久久久大尺度免费视频| 2022亚洲国产成人精品| 国产精品久久久久久精品电影小说| 99久久精品一区二区三区| 日韩大片免费观看网站| 十八禁高潮呻吟视频| 中文欧美无线码| 男女啪啪激烈高潮av片| 免费久久久久久久精品成人欧美视频 | av在线播放精品| 乱码一卡2卡4卡精品| 亚洲美女视频黄频| 国内精品宾馆在线| 精品国产露脸久久av麻豆| 成人手机av| 好男人视频免费观看在线| 久久av网站| 人人澡人人妻人| 国产精品人妻久久久久久| 成人影院久久| 亚洲精华国产精华液的使用体验| 国产欧美另类精品又又久久亚洲欧美| 日本91视频免费播放| 十分钟在线观看高清视频www| 欧美成人午夜免费资源| 日本与韩国留学比较| 国产欧美日韩综合在线一区二区| 亚洲欧美一区二区三区国产| 国产成人精品久久久久久| 男男h啪啪无遮挡| av有码第一页| 亚洲欧美日韩卡通动漫| 亚洲国产精品999| 丰满乱子伦码专区| 亚洲人成网站在线观看播放| 男女边吃奶边做爰视频| h视频一区二区三区| 2018国产大陆天天弄谢| 国产色婷婷99| 国产国拍精品亚洲av在线观看| 少妇精品久久久久久久| 亚洲一级一片aⅴ在线观看| 十八禁高潮呻吟视频| 美女脱内裤让男人舔精品视频| 久久久久久久久大av| 免费久久久久久久精品成人欧美视频 | 日日啪夜夜爽| 汤姆久久久久久久影院中文字幕| 熟女人妻精品中文字幕| 三级国产精品欧美在线观看| a级毛片黄视频| 国产精品 国内视频| 男人爽女人下面视频在线观看| 亚洲av日韩在线播放| 插阴视频在线观看视频| 亚洲性久久影院| 在线观看三级黄色| 视频中文字幕在线观看| 久久久久久久亚洲中文字幕| 免费久久久久久久精品成人欧美视频 | 91久久精品国产一区二区成人| 人人妻人人爽人人添夜夜欢视频| 丰满乱子伦码专区| 少妇的逼水好多| 2018国产大陆天天弄谢| 久久精品熟女亚洲av麻豆精品| 啦啦啦在线观看免费高清www| 99久久人妻综合| 免费黄频网站在线观看国产| 久久99精品国语久久久| 精品久久久精品久久久| 久久久久久久久久久免费av| 亚洲色图综合在线观看| 国产精品99久久久久久久久| 免费观看性生交大片5| 欧美激情极品国产一区二区三区 | 国产精品一区二区在线不卡| 97超碰精品成人国产| 亚洲欧洲精品一区二区精品久久久 | 91精品三级在线观看| 如日韩欧美国产精品一区二区三区 | 欧美日韩精品成人综合77777| 久久精品国产亚洲网站| 亚洲av日韩在线播放| 王馨瑶露胸无遮挡在线观看| 久久精品人人爽人人爽视色| 免费观看无遮挡的男女| 成人免费观看视频高清| 人妻人人澡人人爽人人| 成年美女黄网站色视频大全免费 | 久久人妻熟女aⅴ| 老司机亚洲免费影院| 熟女电影av网| 水蜜桃什么品种好| 精品久久久久久电影网| 国产不卡av网站在线观看| 欧美日韩国产mv在线观看视频| av网站免费在线观看视频| 久久久久久人妻| 黑人猛操日本美女一级片| 人人妻人人澡人人爽人人夜夜| 久久精品夜色国产| 亚洲av.av天堂| 亚洲精品一二三| 亚洲欧美清纯卡通| 亚洲熟女精品中文字幕| 久久精品熟女亚洲av麻豆精品| 国产成人免费无遮挡视频| 中文天堂在线官网| 亚洲图色成人| 狂野欧美激情性xxxx在线观看| 亚洲成人av在线免费| 韩国高清视频一区二区三区| 国产精品人妻久久久影院| 国产在线视频一区二区| a 毛片基地| 精品久久久久久电影网| 人妻一区二区av| 国产精品久久久久成人av| 哪个播放器可以免费观看大片| 精品视频人人做人人爽| 精品少妇内射三级| 久久99精品国语久久久| 精品久久久噜噜| 黄色毛片三级朝国网站| 亚洲av国产av综合av卡| 国产日韩欧美亚洲二区| 午夜免费鲁丝| 美女大奶头黄色视频| 久久精品国产鲁丝片午夜精品| 一级二级三级毛片免费看| 国产免费福利视频在线观看| 午夜精品国产一区二区电影| av一本久久久久| 日本免费在线观看一区| 中文字幕亚洲精品专区| 日韩强制内射视频| 欧美日韩一区二区视频在线观看视频在线| 另类精品久久| 一级片'在线观看视频| 国产成人精品在线电影| 久久精品久久久久久久性| 国产成人精品婷婷| 国产成人精品在线电影| 妹子高潮喷水视频| 午夜久久久在线观看| 欧美日韩综合久久久久久| 婷婷成人精品国产| 欧美少妇被猛烈插入视频| 一级毛片黄色毛片免费观看视频| 欧美 日韩 精品 国产| 激情五月婷婷亚洲| 肉色欧美久久久久久久蜜桃| 国产欧美日韩一区二区三区在线 | 黄色一级大片看看| 国产精品熟女久久久久浪| freevideosex欧美| 丝袜在线中文字幕| 国产伦理片在线播放av一区| 欧美精品国产亚洲| 99久久精品一区二区三区| 欧美另类一区| 日日摸夜夜添夜夜添av毛片| 久久久久久久精品精品| 大陆偷拍与自拍| 国产成人aa在线观看| av免费观看日本| 老女人水多毛片| 日韩强制内射视频| 国产精品.久久久|