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

    基于Runge-Kutta法的輸氣管道泄漏定位技術(shù)研究

    2020-11-17 05:47:04吳波楊旭東張小芳司宗慶牟燕梁昌晶
    油氣田地面工程 2020年11期
    關(guān)鍵詞:管輸定位精度起點(diǎn)

    吳波 楊旭東 張小芳 司宗慶 牟燕 梁昌晶

    1中國(guó)石油華北油田公司第三采油廠

    2中國(guó)石油華北油田公司第一采油廠

    3中國(guó)石油華北油田公司二連分公司

    輸氣管道的泄漏定位技術(shù)是管道完整性管理的重要組成部分,管道防腐層的缺失、陰極保護(hù)的失效、土壤介質(zhì)對(duì)管道外壁的腐蝕,以及管輸介質(zhì)對(duì)管道內(nèi)壁產(chǎn)生的腐蝕,均會(huì)引起管道不同程度的泄漏[1]。據(jù)統(tǒng)計(jì),1995—2005 年,中國(guó)石油管道分公司每年大約產(chǎn)生12~15 起重大輸氣管道泄漏事故,其中2004 年因腐蝕引起的管道泄漏事故就多達(dá)20余起[2-3]。目前,已有多種針對(duì)管道泄漏的檢測(cè)與定位技術(shù),但每種檢測(cè)方法的原理和手段差別較大,同時(shí)由于管輸天然氣組分的多樣性,敷設(shè)方式的差異性(如埋地、直敷、架空、穿跨越等),以及泄漏方式的多元性(如滲透、穿孔、斷裂等),使得各種檢測(cè)方法在定位精度和檢測(cè)時(shí)間上都存在一定的局限性[4-5]。

    韓玲娟等[6]采用分布式光纖傳感器對(duì)水下輸氣管道的泄漏進(jìn)行了有效定位,采用小波神經(jīng)網(wǎng)絡(luò)對(duì)信號(hào)進(jìn)行了有效提取,用聚類分析判斷管道泄漏的嚴(yán)重程度,并采用最小二乘法擬合信號(hào)參數(shù)進(jìn)行定位,該方法與傳統(tǒng)定位方法相比誤報(bào)率下降了6.6%;焦敬品等[7]通過對(duì)20種不同的泄漏信號(hào)進(jìn)行識(shí)別,運(yùn)用BP 神經(jīng)網(wǎng)絡(luò)進(jìn)行泄漏點(diǎn)的檢測(cè)和定位,整體識(shí)別率≥92.5%;田野[8]設(shè)計(jì)了一套基于次生波原理的泄漏監(jiān)測(cè)系統(tǒng),結(jié)合了比例-積分-微分(PID)整定調(diào)節(jié)和自適應(yīng)算法,對(duì)輸氣管道的泄漏定位誤差≤20 m;李鳳等[9]將負(fù)壓波和聲波技術(shù)進(jìn)行了耦合,通過前處理軟件(ICEM)對(duì)管道進(jìn)行了網(wǎng)格劃分,應(yīng)用計(jì)算流體力學(xué)(CFD)軟件進(jìn)行了流場(chǎng)變化分析,有效解決了單一信號(hào)噪聲大、魯棒性差的缺點(diǎn)。綜上所述,管道的定位檢測(cè)需要在首末點(diǎn)設(shè)置大量的傳感器,并且采用多種軟件進(jìn)行模型構(gòu)建和數(shù)據(jù)模擬,工作量大,一般操作人員較難掌握。

    我國(guó)大部分輸氣管道由于管道沿程高差的區(qū)別和地理跨度的差異,屬于非等溫輸氣管道,隨著計(jì)算機(jī)科學(xué)和數(shù)值分析理論的不斷發(fā)展,通過現(xiàn)有傳感器測(cè)量的溫度、壓力和流量數(shù)據(jù),采用改進(jìn)的變步長(zhǎng)龍格-庫(kù)塔(Runge-Kutta)法對(duì)非等溫輸氣管道進(jìn)行檢測(cè)和定位,同時(shí)通過實(shí)例進(jìn)行仿真計(jì)算,并與傳統(tǒng)的負(fù)壓波檢測(cè)結(jié)果進(jìn)行對(duì)比,以期可以快速、準(zhǔn)確地對(duì)泄漏點(diǎn)進(jìn)行定位。

    1 模型建立

    1.1 輸氣管道模型

    將輸氣管道中介質(zhì)的流動(dòng)視為穩(wěn)態(tài)流動(dòng),而穩(wěn)態(tài)流動(dòng)與時(shí)間變量無關(guān),根據(jù)流體力學(xué)和熱力學(xué)相關(guān)公式建立沿管輸方向的穩(wěn)態(tài)流動(dòng)方程[10-11]。其中連續(xù)性方程、動(dòng)量方程、能量方程和狀態(tài)方程分別為

    式中:ρ為管輸氣體密度,kg/m3;ω為氣體流速,m/s;x為管道長(zhǎng)度,m;p為管輸氣體壓力,Pa;g為重力加速度,取9.8 m/s2;θ為管道傾角,rand;λ為水力摩阻系數(shù);d為管道內(nèi)徑,m;h為管輸氣體的焓,J/kg;s為管道的位置高度,m;Q為單位質(zhì)量氣體向外界釋放的能量,J/kg;T為管輸氣體溫度,K。

    在已知?dú)怏w組分及溫度、壓力的前提下,可以由實(shí)際氣體狀態(tài)方程對(duì)氣體密度和壓縮系數(shù)進(jìn)行求解,比較實(shí)用的狀態(tài)方程有RK、SRK、BWRS 和PR方程,在此采用BWRS方程進(jìn)行計(jì)算[12]。

    1.2 Runge-Kutta法

    該算法是一種在工程上被廣泛應(yīng)用的高精度單步算法,主要原理是通過對(duì)常微分方程組進(jìn)行離散化處理,求解該方程的初值;隨后選擇不同的步長(zhǎng),在之間計(jì)算多個(gè)點(diǎn)的斜率值,將其加權(quán)平均作為平均斜率。為了便于計(jì)算,通常采用較為經(jīng)典的四階Runge-Kutta 進(jìn)行計(jì)算[13],通過與泰勒級(jí)數(shù)展開式進(jìn)行對(duì)比后求解未知數(shù),簡(jiǎn)化后的公式如下

    式中:xi、yi分別為自變量和因變量;yi+1為下一個(gè)步長(zhǎng)的因變量;h為步長(zhǎng)系數(shù);K1、K2、K3、K4為待定系數(shù),表示輸出變量的一階倒數(shù)。

    將管道分成n段,每段管道長(zhǎng)度定為Δx(代表每段管道的微元)。Runge-Kutta法的關(guān)鍵是對(duì)于Δx的合理選擇,單從每一個(gè)截?cái)鄟砜?,步長(zhǎng)越小,截?cái)嗾`差越小,但隨著步長(zhǎng)的不斷減少,計(jì)算時(shí)間也在不斷延長(zhǎng),一定范圍內(nèi)的求解步數(shù)也會(huì)增加。不同的迭代計(jì)算會(huì)導(dǎo)致舍入誤差的大量累計(jì),在計(jì)算效率和精度上有很大缺陷,因此首先設(shè)置較大的計(jì)算步長(zhǎng),在整個(gè)管段內(nèi)進(jìn)行泄漏位置的初選,隨后逐步縮小泄漏定位的搜索范圍,直到滿足定位精度的要求,一般取0.5%~1%的初始步長(zhǎng)進(jìn)行運(yùn)算。

    在經(jīng)典四階Runge-Kutta 算法的基礎(chǔ)上采用精細(xì)Runge-Kutta 法計(jì)算,該方法融合了精細(xì)積分法和Runge-Kutta 法的優(yōu)點(diǎn),可以將特解的數(shù)值代入Runge-Kutta公式中積分,并利用Runge-Kutta法對(duì)數(shù)值積分點(diǎn)的未知參數(shù)進(jìn)行預(yù)估,得到更精細(xì)的算法。根據(jù)現(xiàn)場(chǎng)傳感器采集到的管道起點(diǎn)溫度Ti、壓力pi和流量mi,計(jì)算從起點(diǎn)到末點(diǎn)管道各截面的溫度、壓力和流量數(shù)值。首先計(jì)算第一個(gè)微元Δx1的斜率k1,從第一個(gè)積分點(diǎn)出發(fā),前進(jìn)半個(gè)步長(zhǎng),得到第二個(gè)積分點(diǎn)的預(yù)測(cè)斜率k2;然后以k2為基礎(chǔ),再次從第一個(gè)積分點(diǎn)出發(fā),前進(jìn)半個(gè)步長(zhǎng),重新計(jì)算第二個(gè)積分點(diǎn)的預(yù)測(cè)斜率k′2;將兩次計(jì)算的k2和k′2進(jìn)行加權(quán)平均,作為第二個(gè)積分點(diǎn)的近似斜率;以此類推,直到計(jì)算出最后一個(gè)積分點(diǎn)也就是管道末點(diǎn)的溫度、壓力和流量數(shù)值。其中以第一個(gè)微元管道Δx1為例,運(yùn)用Runge-Kutta法求解公式(1)~(4)來計(jì)算管道截面的參數(shù),公式如下

    以此類推,得到管道末點(diǎn)的溫度Ti+1、壓力pi+1和流量mi+1的計(jì)算公式如下

    從管道起點(diǎn)開始,運(yùn)用Matlab 對(duì)公式(7)進(jìn)行反復(fù)的迭代計(jì)算,即可完成從管道起點(diǎn)到末點(diǎn)不同截面的溫度、壓力和流量的計(jì)算;從末點(diǎn)到起點(diǎn)的參數(shù)計(jì)算也類似,只是將公式(2)~(3)的摩擦阻力項(xiàng)和熱損失項(xiàng)進(jìn)行調(diào)整,具體見式(8)~(9),這種模型非常適合非等溫輸氣管道的模擬和仿真。

    2 泄漏定位算法

    通過現(xiàn)場(chǎng)傳感器采集到的起點(diǎn)溫度Ti、壓力pi和流量mi,依據(jù)Runge-Kutta法計(jì)算末點(diǎn)管道的溫度、壓力和流量,同樣通過采集到的末點(diǎn)溫度Tn、壓力pn和流量mn,反算起點(diǎn)管道的溫度、壓力和流量,得到圖1。

    圖1 泄漏定位原理Fig.1 Leak location principle

    假設(shè)在某一時(shí)刻t,由管道末點(diǎn)計(jì)算到管道起點(diǎn)的參數(shù)為Tnt、pnt、mnt,將判斷泄漏事故發(fā)生的壓力閾值定義為εt,則判斷泄漏的邊界條件為

    當(dāng)公式(10)成立,表示管內(nèi)無泄漏點(diǎn);當(dāng)存在泄漏點(diǎn)時(shí),管道內(nèi)各點(diǎn)的相關(guān)參數(shù)會(huì)偏離計(jì)算結(jié)果,式(10)不再成立。當(dāng)pnt<pn時(shí),管道內(nèi)存在泄漏點(diǎn)。在管道的運(yùn)行過程中,難免會(huì)出現(xiàn)關(guān)停壓縮機(jī)或停輸再啟動(dòng)等工況,此時(shí)的監(jiān)測(cè)數(shù)據(jù)會(huì)有較大波動(dòng),容易產(chǎn)生誤報(bào),因此不考慮改變操作引起的擾動(dòng)在時(shí)間內(nèi)的報(bào)警,t1為泄漏點(diǎn)傳播到起點(diǎn)的時(shí)間,t2為泄漏點(diǎn)傳播到末點(diǎn)的時(shí)間,t3為擾動(dòng)的熵變傳播到末點(diǎn)的時(shí)間,n為劃分管道的段數(shù),Δt為報(bào)警閾值。

    圖1 中L1代表未發(fā)生泄漏時(shí)的壓降,L2代表泄漏發(fā)生后從起點(diǎn)計(jì)算到末點(diǎn)的壓降,L3代表泄漏發(fā)生后從末點(diǎn)計(jì)算到起點(diǎn)的壓降。當(dāng)發(fā)生泄漏時(shí),泄漏點(diǎn)處氣體的密度和壓縮系數(shù)會(huì)發(fā)生變化,壓力也會(huì)發(fā)生瞬變。對(duì)于L2,泄漏點(diǎn)之前的仿真結(jié)果無限接近L1,泄漏點(diǎn)之后的仿真結(jié)果偏離L1,且距離x0越遠(yuǎn),偏離程度越嚴(yán)重;對(duì)于L3,泄漏點(diǎn)之后的仿真結(jié)果接近L1,泄漏點(diǎn)之前的仿真結(jié)果偏離L1,L2、L3的交點(diǎn)為x0即為泄漏點(diǎn)。泄漏定位流程圖見圖2。

    3 實(shí)例計(jì)算

    華北油田采油二廠從別古莊到霸一聯(lián)的輸氣管道全長(zhǎng)2.7 km,管道規(guī)格為Φ213 mm×5.5 mm,管道輸量200 m3/min,實(shí)際運(yùn)行壓力0.4~0.8 MPa,管輸介質(zhì)成分CH4占98%,另外還有少量的C2H6、C3H8等。利用管線全面清管期間,在一些特定位置(充分考慮安全、環(huán)保的要求)人工制造幾處泄漏點(diǎn),通過球閥控制泄漏過程。為了避免放氣時(shí)間與信號(hào)采集不同步,先進(jìn)行閥門放氣,一段時(shí)間后關(guān)閉閥門,設(shè)置首站壓力0.65 MPa。圖3 和圖4 分別為泄漏點(diǎn)位于40.5 m時(shí),發(fā)生泄漏后管道起點(diǎn)和末點(diǎn)的壓力隨時(shí)間變化的曲線。

    圖2 泄漏定位流程圖Fig.2 Flow chart of leak location

    圖3 起點(diǎn)壓力隨時(shí)間變化曲線Fig.3 Change curve of start point pressure with time

    圖4 末點(diǎn)壓力隨時(shí)間變化曲線Fig.4 Change curve of final point pressure change with time

    由圖3 和圖4 可知,在15 s 處管道發(fā)生泄漏,15~25 s 的時(shí)間段內(nèi)首末點(diǎn)壓力急劇下降,而后趨于平穩(wěn),可見整個(gè)泄漏過程經(jīng)歷了穩(wěn)態(tài)—快瞬態(tài)—慢瞬態(tài)—穩(wěn)態(tài)四個(gè)過程,流體從恒定流狀態(tài)向瞬變流狀態(tài)轉(zhuǎn)變,穩(wěn)定之后的壓力值低于泄漏前的壓力值。為了更好地驗(yàn)證改進(jìn)后的Runge-Kutta 法對(duì)泄漏定位預(yù)測(cè)的準(zhǔn)確性,控制閥門開度,以不同的泄漏率(即泄漏孔徑與管徑的比值)對(duì)模型進(jìn)行驗(yàn)證,同時(shí)與傳統(tǒng)的負(fù)壓超聲波定位精度進(jìn)行對(duì)比,結(jié)果如表1所示。

    表1 不同泄漏率下泄漏定位結(jié)果對(duì)比Tab.1 Comparison of leakage location results under different leakage rates

    由表1可知,在泄漏位置一定的情況下,隨著泄漏孔徑和泄漏率的增加,Runge-Kutta 法和負(fù)壓波定位的定位誤差都在逐漸變小,通過結(jié)果對(duì)比,Runge-Kutta 法的定位誤差為0.4~2.5 m,小于傳統(tǒng)的負(fù)壓波定位誤差(5~12 m)。

    地形地貌及管輸介質(zhì)成分的不同決定了實(shí)際工況下的輸氣管道大多數(shù)為非等溫管道,為考察首末段溫差對(duì)定位精度的影響,在泄漏位置和泄漏率一定的前提下,通過提高首站和末站的溫度,在泄漏率5%的條件下用Runge-Kutta法和負(fù)壓波對(duì)泄漏進(jìn)行定位,結(jié)果如表2所示。

    由表2可知,隨著首末點(diǎn)溫差的不斷增加,兩種方法的定位誤差也隨之增大,但Runge-Kutta 法在定位精度上明顯優(yōu)于負(fù)壓波定位,因此在溫差較大的情況下,Runge-Kutta 法可用于非等溫輸氣管道的泄漏檢測(cè)和定位。

    管道輸量的大小對(duì)于不同泄漏孔徑也有一定影響??刂剖渍玖髁繛?50~300 m3/min,在泄漏率5%條件下,首站壓力設(shè)置為0.5~0.7 MPa,對(duì)326 m處的泄漏點(diǎn)進(jìn)行定位,結(jié)果如表3所示。

    表2 首末點(diǎn)溫差變化對(duì)定位精度的影響Tab.2 Influence of temperature difference between start and final point on positioning accuracy

    表3 不同管道輸量對(duì)定位精度的影響Tab.3 Influence of different pipeline volumes on positioning accuracy

    由表3可知,隨著管道輸量的不斷增加,定位誤差整體減小。這是因?yàn)檩斄吭酱螅艿纼?nèi)的氣體流速越大,在泄漏孔處的流出速度越大,壓力波動(dòng)越明顯。

    綜上所述,對(duì)比不同情況下的定位誤差,對(duì)定位精度的影響因素排序由大至小是管道輸量、泄漏率、首末點(diǎn)溫差。

    4 結(jié)束語(yǔ)

    (1)通過采集到的首末點(diǎn)溫度、壓力和流量數(shù)據(jù),應(yīng)用改進(jìn)后的Runge-Kutta 法可以計(jì)算管道各點(diǎn)的壓力分布,進(jìn)而對(duì)泄漏點(diǎn)進(jìn)行有效定位,該方法原理簡(jiǎn)單,對(duì)現(xiàn)場(chǎng)操作人員的專業(yè)要求不高,易于實(shí)現(xiàn),與傳統(tǒng)的負(fù)壓波相比具有一定的優(yōu)越性。

    (2)泄漏定位技術(shù)涉及多學(xué)科的融合和交叉,不應(yīng)只注重現(xiàn)場(chǎng)傳感器的數(shù)據(jù),隨著機(jī)器學(xué)習(xí)、數(shù)值分析的發(fā)展,今后應(yīng)加強(qiáng)各學(xué)科之間的融合,提高定位精度。

    (2)該方法只適用于干氣管道,對(duì)于含液較多的管道由于氣液流速不一致、氣液滑脫等現(xiàn)象并不適用,同時(shí)對(duì)于同一時(shí)刻有多處泄漏點(diǎn)的數(shù)據(jù)擾動(dòng)考慮不足,今后應(yīng)深入研究,加強(qiáng)模型的完善以提高定位精度。

    猜你喜歡
    管輸定位精度起點(diǎn)
    北斗定位精度可達(dá)兩三米
    軍事文摘(2023年4期)2023-04-05 13:57:35
    儀長(zhǎng)管輸原油加工高等級(jí)道路石油瀝青試驗(yàn)研究
    石油瀝青(2022年1期)2022-04-19 13:10:28
    席海宏:管輸+營(yíng)銷“雙冠王”
    GPS定位精度研究
    組合導(dǎo)航的AGV定位精度的改善
    弄清楚“起點(diǎn)”前面有多少
    起點(diǎn)
    我的“新”起點(diǎn)
    管輸天然氣泄露模型定量分析應(yīng)用
    新年的起點(diǎn)
    中文字幕av在线有码专区| 国产精品一区二区三区四区久久| 美女内射精品一级片tv| 午夜精品国产一区二区电影 | 国产欧美日韩精品一区二区| 亚洲人成网站在线播放欧美日韩| 免费电影在线观看免费观看| 久久久久九九精品影院| 夜夜夜夜夜久久久久| 成人一区二区视频在线观看| 身体一侧抽搐| 少妇的逼好多水| 午夜免费激情av| 国内久久婷婷六月综合欲色啪| 欧美日韩国产亚洲二区| 日本免费一区二区三区高清不卡| 成年女人永久免费观看视频| 男人和女人高潮做爰伦理| 亚洲天堂国产精品一区在线| 成人亚洲欧美一区二区av| 深夜a级毛片| 俄罗斯特黄特色一大片| 麻豆乱淫一区二区| 国产伦在线观看视频一区| 又粗又爽又猛毛片免费看| 亚洲熟妇熟女久久| 99热全是精品| 国产伦在线观看视频一区| 日产精品乱码卡一卡2卡三| 久久久久久国产a免费观看| 亚洲欧美日韩高清专用| 日韩一区二区视频免费看| av在线蜜桃| 国产成年人精品一区二区| 高清毛片免费观看视频网站| 蜜桃久久精品国产亚洲av| 黄色配什么色好看| 又粗又爽又猛毛片免费看| 婷婷六月久久综合丁香| 香蕉av资源在线| 天天一区二区日本电影三级| 国产精品一区www在线观看| 好男人在线观看高清免费视频| 亚洲av电影不卡..在线观看| 人人妻人人看人人澡| 色综合色国产| 在线a可以看的网站| 波多野结衣高清无吗| 男女做爰动态图高潮gif福利片| 亚洲精品国产成人久久av| 亚洲五月天丁香| 又爽又黄无遮挡网站| 国产色婷婷99| 亚洲欧美成人精品一区二区| 亚洲成人av在线免费| 国产精品一区二区性色av| 国产男靠女视频免费网站| 国产午夜精品久久久久久一区二区三区 | 超碰av人人做人人爽久久| 九九热线精品视视频播放| 日韩欧美免费精品| 一级毛片我不卡| 成人二区视频| 精品国内亚洲2022精品成人| 白带黄色成豆腐渣| av.在线天堂| 少妇丰满av| 免费高清视频大片| 搡老妇女老女人老熟妇| 亚洲国产色片| 如何舔出高潮| 直男gayav资源| 中文字幕av在线有码专区| 黄色日韩在线| 国产高潮美女av| 日韩欧美 国产精品| 国产一区亚洲一区在线观看| 成人漫画全彩无遮挡| 春色校园在线视频观看| 久久久国产成人免费| 久久精品国产99精品国产亚洲性色| 97超级碰碰碰精品色视频在线观看| 欧美潮喷喷水| 波多野结衣巨乳人妻| 久久久久国产精品人妻aⅴ院| 九九在线视频观看精品| 午夜福利18| 国产精品人妻久久久影院| 免费大片18禁| 免费搜索国产男女视频| 国产亚洲av嫩草精品影院| 热99re8久久精品国产| 中文在线观看免费www的网站| 国内精品一区二区在线观看| 久久精品91蜜桃| 欧美又色又爽又黄视频| 最新中文字幕久久久久| 国产伦精品一区二区三区视频9| 精品乱码久久久久久99久播| 俄罗斯特黄特色一大片| 久久久久九九精品影院| 亚洲人成网站在线播| 精品国内亚洲2022精品成人| 特级一级黄色大片| 亚洲中文字幕一区二区三区有码在线看| av中文乱码字幕在线| 亚洲精品粉嫩美女一区| 在线观看美女被高潮喷水网站| 午夜老司机福利剧场| aaaaa片日本免费| 亚洲av成人精品一区久久| 亚洲丝袜综合中文字幕| av福利片在线观看| 自拍偷自拍亚洲精品老妇| 日韩欧美三级三区| 亚洲av第一区精品v没综合| 一本一本综合久久| 在线播放国产精品三级| 美女被艹到高潮喷水动态| 三级经典国产精品| 日韩精品青青久久久久久| 婷婷色综合大香蕉| 九九热线精品视视频播放| 禁无遮挡网站| 蜜桃久久精品国产亚洲av| 成人美女网站在线观看视频| 国产伦在线观看视频一区| 一区二区三区高清视频在线| 久久精品夜夜夜夜夜久久蜜豆| 国产男人的电影天堂91| 毛片女人毛片| 久久久国产成人精品二区| av在线老鸭窝| .国产精品久久| 三级男女做爰猛烈吃奶摸视频| 亚洲精品一区av在线观看| 国产私拍福利视频在线观看| 听说在线观看完整版免费高清| 精品无人区乱码1区二区| 国产在线精品亚洲第一网站| 日本熟妇午夜| 色哟哟哟哟哟哟| 免费av毛片视频| 69人妻影院| 国产亚洲欧美98| 国产欧美日韩精品一区二区| 国产熟女欧美一区二区| 波多野结衣高清无吗| 免费电影在线观看免费观看| 丝袜美腿在线中文| 99热这里只有精品一区| 欧美日本视频| 午夜福利在线在线| 99久久中文字幕三级久久日本| av卡一久久| 成年免费大片在线观看| 国产精品一区二区免费欧美| 国产伦在线观看视频一区| 国产午夜福利久久久久久| 国产 一区精品| 国产男人的电影天堂91| 国产精品久久久久久av不卡| aaaaa片日本免费| 99精品在免费线老司机午夜| 夜夜看夜夜爽夜夜摸| 国产精品久久久久久精品电影| 国产黄色小视频在线观看| 人妻少妇偷人精品九色| 99精品在免费线老司机午夜| 全区人妻精品视频| 人人妻人人澡欧美一区二区| АⅤ资源中文在线天堂| 久久精品人妻少妇| 3wmmmm亚洲av在线观看| 免费大片18禁| 欧美成人一区二区免费高清观看| 亚洲七黄色美女视频| 欧美xxxx性猛交bbbb| 一级黄片播放器| 精品人妻一区二区三区麻豆 | 亚洲成人久久性| 一a级毛片在线观看| 欧美性猛交黑人性爽| 99久国产av精品国产电影| 国产91av在线免费观看| 欧美+日韩+精品| 不卡一级毛片| 免费av毛片视频| 一区二区三区四区激情视频 | 99热这里只有精品一区| 91狼人影院| 亚洲精品日韩在线中文字幕 | 久久久久久国产a免费观看| 日韩 亚洲 欧美在线| 男女那种视频在线观看| 看免费成人av毛片| av天堂在线播放| 久久久精品欧美日韩精品| 亚洲国产色片| 亚洲三级黄色毛片| 国产高潮美女av| 欧美精品国产亚洲| 日韩高清综合在线| 成人亚洲精品av一区二区| 国产免费男女视频| 精品国产三级普通话版| 男插女下体视频免费在线播放| 中文字幕免费在线视频6| 日日摸夜夜添夜夜添小说| 你懂的网址亚洲精品在线观看 | 久久久a久久爽久久v久久| 精品一区二区三区人妻视频| 国产精品日韩av在线免费观看| 国产成人aa在线观看| 精品久久久久久久人妻蜜臀av| 成人av一区二区三区在线看| 久久综合国产亚洲精品| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站 | www日本黄色视频网| 激情 狠狠 欧美| 91精品国产九色| 国产精品亚洲美女久久久| 在线观看av片永久免费下载| 精品乱码久久久久久99久播| 日韩亚洲欧美综合| 直男gayav资源| 在线观看免费视频日本深夜| 亚洲av.av天堂| 男人的好看免费观看在线视频| 日韩欧美国产在线观看| 女同久久另类99精品国产91| 日本与韩国留学比较| 观看美女的网站| 亚洲精品456在线播放app| 成人高潮视频无遮挡免费网站| 国产精品福利在线免费观看| 欧美性猛交╳xxx乱大交人| 国内精品宾馆在线| 午夜福利成人在线免费观看| 熟妇人妻久久中文字幕3abv| 自拍偷自拍亚洲精品老妇| 少妇裸体淫交视频免费看高清| 性欧美人与动物交配| 国产精品国产高清国产av| 少妇人妻一区二区三区视频| 亚洲av熟女| 成人鲁丝片一二三区免费| 欧美激情国产日韩精品一区| 亚洲成人中文字幕在线播放| 日日摸夜夜添夜夜爱| 亚洲第一电影网av| 日韩精品中文字幕看吧| 美女免费视频网站| 岛国在线免费视频观看| 热99re8久久精品国产| 可以在线观看毛片的网站| 国产精品一区二区性色av| 欧美精品国产亚洲| 可以在线观看的亚洲视频| 午夜视频国产福利| 在线观看av片永久免费下载| 中国国产av一级| 中出人妻视频一区二区| 色5月婷婷丁香| eeuss影院久久| 美女内射精品一级片tv| 丝袜喷水一区| ponron亚洲| 日韩强制内射视频| 精品午夜福利视频在线观看一区| 在线观看66精品国产| 美女黄网站色视频| 1024手机看黄色片| 久久久久久久久中文| 91av网一区二区| 久久韩国三级中文字幕| 国产精品一区二区免费欧美| 12—13女人毛片做爰片一| 国产精品三级大全| 在线免费观看不下载黄p国产| 91在线观看av| 免费看a级黄色片| 国产高清有码在线观看视频| 精品一区二区三区av网在线观看| 日本黄大片高清| 床上黄色一级片| 亚洲最大成人手机在线| 久久久久久久久中文| 国产成人91sexporn| 男插女下体视频免费在线播放| 老师上课跳d突然被开到最大视频| 婷婷精品国产亚洲av| 亚州av有码| 在现免费观看毛片| 亚洲av中文字字幕乱码综合| 国国产精品蜜臀av免费| 色播亚洲综合网| 久久久久精品国产欧美久久久| 亚洲中文字幕日韩| 乱人视频在线观看| 校园春色视频在线观看| 高清午夜精品一区二区三区 | 亚洲国产欧洲综合997久久,| 男人和女人高潮做爰伦理| 91在线观看av| 男女下面进入的视频免费午夜| 成人无遮挡网站| 人妻丰满熟妇av一区二区三区| 国产三级在线视频| 我要搜黄色片| 中文在线观看免费www的网站| 欧美+亚洲+日韩+国产| 69av精品久久久久久| 亚洲人成网站高清观看| 久久精品人妻少妇| 欧美在线一区亚洲| 午夜精品一区二区三区免费看| 悠悠久久av| 日本免费a在线| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区 | 国产一区亚洲一区在线观看| 1024手机看黄色片| 午夜视频国产福利| 一区福利在线观看| 久久人人爽人人片av| 欧美激情久久久久久爽电影| av免费在线看不卡| 久久精品夜色国产| 不卡视频在线观看欧美| 99久久中文字幕三级久久日本| 最近手机中文字幕大全| 免费一级毛片在线播放高清视频| 国产一区二区三区av在线 | 成人无遮挡网站| 免费观看精品视频网站| 亚洲精品在线观看二区| 久久天躁狠狠躁夜夜2o2o| 可以在线观看毛片的网站| 精华霜和精华液先用哪个| 最近手机中文字幕大全| 内射极品少妇av片p| 99热只有精品国产| 99热这里只有是精品在线观看| av福利片在线观看| 91麻豆精品激情在线观看国产| 亚洲欧美中文字幕日韩二区| 内地一区二区视频在线| 在线观看免费视频日本深夜| 国产精品日韩av在线免费观看| 色视频www国产| 三级男女做爰猛烈吃奶摸视频| 国产精品女同一区二区软件| 日本一本二区三区精品| 国产午夜精品久久久久久一区二区三区 | 国产精华一区二区三区| 午夜免费男女啪啪视频观看 | 色尼玛亚洲综合影院| 久久精品91蜜桃| 舔av片在线| 99热这里只有是精品50| av天堂在线播放| 亚洲精品久久国产高清桃花| 淫妇啪啪啪对白视频| 美女免费视频网站| 久久久国产成人免费| 国产精品国产高清国产av| 久久久久久国产a免费观看| 久久草成人影院| 三级国产精品欧美在线观看| 一区二区三区四区激情视频 | 亚洲av中文字字幕乱码综合| av免费在线看不卡| 色视频www国产| 亚洲欧美精品自产自拍| 亚洲性夜色夜夜综合| 直男gayav资源| av在线观看视频网站免费| 久久久色成人| 国产伦一二天堂av在线观看| 成人三级黄色视频| 一本一本综合久久| 亚洲色图av天堂| 乱人视频在线观看| 中文在线观看免费www的网站| 村上凉子中文字幕在线| 毛片女人毛片| 一个人看的www免费观看视频| 久久久久久久午夜电影| 欧美成人精品欧美一级黄| 一区二区三区高清视频在线| 亚洲国产精品成人综合色| 成人美女网站在线观看视频| 精华霜和精华液先用哪个| 一进一出抽搐动态| 欧美日韩综合久久久久久| 日韩强制内射视频| 久久久久精品国产欧美久久久| 此物有八面人人有两片| 毛片一级片免费看久久久久| 免费一级毛片在线播放高清视频| 热99在线观看视频| 国产黄色视频一区二区在线观看 | 午夜福利视频1000在线观看| 又黄又爽又免费观看的视频| 久久久久久大精品| 99热这里只有精品一区| 久久欧美精品欧美久久欧美| 毛片一级片免费看久久久久| 在线天堂最新版资源| 国内精品宾馆在线| 欧美+亚洲+日韩+国产| 亚洲成a人片在线一区二区| 美女大奶头视频| 精品99又大又爽又粗少妇毛片| 男女下面进入的视频免费午夜| 在线播放国产精品三级| 国产欧美日韩精品亚洲av| 看片在线看免费视频| 亚洲性夜色夜夜综合| 国产探花极品一区二区| 中国美女看黄片| 全区人妻精品视频| 亚洲成人久久爱视频| 亚洲真实伦在线观看| 亚洲精品456在线播放app| 男女下面进入的视频免费午夜| av黄色大香蕉| 国产一区二区在线av高清观看| 亚洲成av人片在线播放无| 久久综合国产亚洲精品| 国产黄片美女视频| 老司机影院成人| 亚洲精品日韩在线中文字幕 | 久久精品国产亚洲av涩爱 | 国产av一区在线观看免费| 日韩欧美三级三区| 久久久久久久久中文| 国产黄色小视频在线观看| 国产黄片美女视频| 精品欧美国产一区二区三| 国产三级中文精品| 欧美3d第一页| 看非洲黑人一级黄片| 午夜久久久久精精品| 国产亚洲91精品色在线| 尾随美女入室| 色播亚洲综合网| av在线蜜桃| 国产黄色视频一区二区在线观看 | 久久久久久九九精品二区国产| a级毛色黄片| 久久精品久久久久久噜噜老黄 | 精品99又大又爽又粗少妇毛片| 久久久久久伊人网av| 亚洲在线观看片| 综合色av麻豆| 成人高潮视频无遮挡免费网站| 久久久久国内视频| 午夜福利高清视频| 成人二区视频| 免费电影在线观看免费观看| 亚洲中文日韩欧美视频| 欧美zozozo另类| 22中文网久久字幕| 99热这里只有是精品在线观看| 在线观看一区二区三区| 欧美一级a爱片免费观看看| 噜噜噜噜噜久久久久久91| 少妇熟女aⅴ在线视频| 国产成人影院久久av| 日韩人妻高清精品专区| 特大巨黑吊av在线直播| 亚洲av美国av| 尾随美女入室| 97热精品久久久久久| 国产精品综合久久久久久久免费| 欧美丝袜亚洲另类| 久久婷婷人人爽人人干人人爱| 国产aⅴ精品一区二区三区波| 欧美一区二区国产精品久久精品| 黄色一级大片看看| 男人和女人高潮做爰伦理| av天堂在线播放| 亚洲美女搞黄在线观看 | 特级一级黄色大片| 一本久久中文字幕| av中文乱码字幕在线| 亚洲精品乱码久久久v下载方式| 亚州av有码| 欧美日本视频| 久久久久国内视频| 神马国产精品三级电影在线观看| 国产91av在线免费观看| 哪里可以看免费的av片| 尾随美女入室| 国产成人影院久久av| 日韩精品青青久久久久久| 中文字幕免费在线视频6| 国产高清视频在线观看网站| 国产免费男女视频| 99久久精品国产国产毛片| av在线观看视频网站免费| 久久久久久久久久久丰满| АⅤ资源中文在线天堂| 久久久久国产网址| 麻豆成人午夜福利视频| 成人三级黄色视频| 久久久精品大字幕| 精品免费久久久久久久清纯| 亚洲国产高清在线一区二区三| 麻豆精品久久久久久蜜桃| 菩萨蛮人人尽说江南好唐韦庄 | 少妇人妻精品综合一区二区 | 男女边吃奶边做爰视频| 婷婷亚洲欧美| 亚洲成人久久爱视频| 久久精品久久久久久噜噜老黄 | 久久人妻av系列| 最近视频中文字幕2019在线8| h日本视频在线播放| 熟女人妻精品中文字幕| 热99在线观看视频| 中文字幕av在线有码专区| 亚洲精品影视一区二区三区av| 天天躁夜夜躁狠狠久久av| 欧美另类亚洲清纯唯美| 国产一区二区在线av高清观看| 日韩精品中文字幕看吧| 亚洲欧美中文字幕日韩二区| 日韩精品有码人妻一区| 真人做人爱边吃奶动态| 久久天躁狠狠躁夜夜2o2o| 男人舔奶头视频| 高清毛片免费看| 有码 亚洲区| 老熟妇乱子伦视频在线观看| 日韩 亚洲 欧美在线| 国产精品av视频在线免费观看| 日日摸夜夜添夜夜添小说| 色播亚洲综合网| 日韩人妻高清精品专区| 久久久久久久久久久丰满| 淫秽高清视频在线观看| 国产视频一区二区在线看| 中文字幕精品亚洲无线码一区| 亚洲美女视频黄频| 免费av观看视频| 性色avwww在线观看| 2021天堂中文幕一二区在线观| 观看免费一级毛片| 舔av片在线| av.在线天堂| 国产精品久久久久久久久免| 国产av不卡久久| 久久久久国产精品人妻aⅴ院| 简卡轻食公司| 男女做爰动态图高潮gif福利片| 人人妻人人看人人澡| 最后的刺客免费高清国语| 村上凉子中文字幕在线| 黄色欧美视频在线观看| av在线播放精品| 午夜福利高清视频| 国产乱人偷精品视频| 欧美性感艳星| 精品一区二区三区人妻视频| 中文字幕免费在线视频6| 欧美性猛交黑人性爽| 国产精品亚洲一级av第二区| 女人被狂操c到高潮| 长腿黑丝高跟| 一个人看的www免费观看视频| 韩国av在线不卡| 亚洲av免费在线观看| 人妻制服诱惑在线中文字幕| 春色校园在线视频观看| 欧美+亚洲+日韩+国产| 中出人妻视频一区二区| 精品人妻视频免费看| 亚洲熟妇熟女久久| 欧美另类亚洲清纯唯美| 久久中文看片网| 美女被艹到高潮喷水动态| 欧美性猛交黑人性爽| 成人高潮视频无遮挡免费网站| 性色avwww在线观看| 国产精品久久久久久久电影| 精品乱码久久久久久99久播| 久久久欧美国产精品| 乱人视频在线观看| 一个人看的www免费观看视频| 少妇猛男粗大的猛烈进出视频 | 舔av片在线| 三级男女做爰猛烈吃奶摸视频| 狂野欧美白嫩少妇大欣赏| 寂寞人妻少妇视频99o| 欧美xxxx性猛交bbbb| www日本黄色视频网| 日本三级黄在线观看| www.色视频.com| 人人妻人人看人人澡| 亚洲精品影视一区二区三区av| 国产白丝娇喘喷水9色精品| 免费观看的影片在线观看| 神马国产精品三级电影在线观看| 欧美成人精品欧美一级黄| 高清日韩中文字幕在线| 午夜久久久久精精品| 日韩强制内射视频| 精品不卡国产一区二区三区| 国产真实乱freesex| 国产精品久久久久久亚洲av鲁大| 午夜精品国产一区二区电影 |