• <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)
    亚洲五月色婷婷综合| 国产欧美日韩综合在线一区二区| 国产精品一区二区三区四区免费观看| 精品99又大又爽又粗少妇毛片| 午夜老司机福利剧场| 97超视频在线观看视频| 国产精品久久久久久久电影| 日韩欧美一区视频在线观看| 人妻少妇偷人精品九色| 一区二区三区四区激情视频| 多毛熟女@视频| 欧美最新免费一区二区三区| 建设人人有责人人尽责人人享有的| 好男人视频免费观看在线| 亚洲欧美成人精品一区二区| 男女无遮挡免费网站观看| 久热这里只有精品99| 久久亚洲国产成人精品v| 日韩在线高清观看一区二区三区| 精品一品国产午夜福利视频| 又黄又爽又刺激的免费视频.| videos熟女内射| 亚洲精品一区蜜桃| 国产深夜福利视频在线观看| 国产精品久久久久久精品古装| 最近2019中文字幕mv第一页| 少妇熟女欧美另类| 亚洲五月色婷婷综合| 一边摸一边做爽爽视频免费| 亚洲,一卡二卡三卡| 国产欧美亚洲国产| 如日韩欧美国产精品一区二区三区 | 欧美3d第一页| 国产日韩一区二区三区精品不卡 | a 毛片基地| 亚洲国产色片| 欧美精品亚洲一区二区| 亚洲人与动物交配视频| 久久精品久久久久久久性| 亚洲美女搞黄在线观看| 十八禁高潮呻吟视频| 又黄又爽又刺激的免费视频.| 精品久久久噜噜| 久久国产精品男人的天堂亚洲 | 春色校园在线视频观看| 久久久久久久精品精品| 精品久久久久久久久亚洲| 亚洲精品日韩av片在线观看| 久久久久精品久久久久真实原创| 欧美bdsm另类| 免费播放大片免费观看视频在线观看| 一区二区av电影网| 美女主播在线视频| 亚洲精品日韩av片在线观看| 日本色播在线视频| 久久精品国产自在天天线| 草草在线视频免费看| 久久精品久久久久久噜噜老黄| 国产一区二区在线观看av| 亚洲av.av天堂| 97超视频在线观看视频| 久久人人爽人人爽人人片va| 日本猛色少妇xxxxx猛交久久| 国国产精品蜜臀av免费| 亚洲国产精品一区三区| 黄片播放在线免费| 久久99热这里只频精品6学生| 亚洲成人手机| 91精品国产国语对白视频| 夜夜骑夜夜射夜夜干| 亚洲国产精品成人久久小说| 2021少妇久久久久久久久久久| 春色校园在线视频观看| 最近的中文字幕免费完整| av女优亚洲男人天堂| 亚洲精品日本国产第一区| 日韩 亚洲 欧美在线| 久久人人爽人人片av| 在线观看www视频免费| 国产一级毛片在线| 内地一区二区视频在线| 在线亚洲精品国产二区图片欧美 | 国产精品无大码| 日本色播在线视频| 日本vs欧美在线观看视频| 国产男人的电影天堂91| 免费久久久久久久精品成人欧美视频 | 妹子高潮喷水视频| 性色av一级| 免费黄网站久久成人精品| 国产深夜福利视频在线观看| 丝袜喷水一区| 欧美一级a爱片免费观看看| 国产淫语在线视频| 热99久久久久精品小说推荐| 亚洲国产精品999| 亚洲国产日韩一区二区| 欧美日韩综合久久久久久| 日本免费在线观看一区| av网站免费在线观看视频| 国产精品成人在线| 国产欧美日韩综合在线一区二区| 日韩三级伦理在线观看| 午夜久久久在线观看| 国产亚洲最大av| 亚洲国产精品一区三区| 日韩伦理黄色片| av播播在线观看一区| av免费观看日本| 久久精品国产自在天天线| 美女主播在线视频| 亚洲综合色网址| 久久综合国产亚洲精品| 免费观看性生交大片5| 国产精品久久久久久av不卡| 午夜福利影视在线免费观看| 丝袜脚勾引网站| 妹子高潮喷水视频| 精品亚洲成国产av| 一区二区日韩欧美中文字幕 | 韩国高清视频一区二区三区| 久久精品国产鲁丝片午夜精品| 久久久久久久亚洲中文字幕| 热99久久久久精品小说推荐| 色94色欧美一区二区| 亚洲av日韩在线播放| 哪个播放器可以免费观看大片| 日日撸夜夜添| 中文字幕最新亚洲高清| 视频中文字幕在线观看| 国产亚洲精品久久久com| 久久精品夜色国产| 久久精品国产a三级三级三级| 精品人妻在线不人妻| 国产精品久久久久久久电影| 日本与韩国留学比较| 国产一区二区在线观看av| 日韩强制内射视频| 性色avwww在线观看| 国产69精品久久久久777片| 亚洲图色成人| 一级毛片aaaaaa免费看小| 亚洲无线观看免费| 国产精品久久久久久久电影| 久久99蜜桃精品久久| 亚洲欧美成人综合另类久久久| 2018国产大陆天天弄谢| 黑人欧美特级aaaaaa片| av电影中文网址| 免费黄色在线免费观看| 国产深夜福利视频在线观看| tube8黄色片| 欧美老熟妇乱子伦牲交| 人体艺术视频欧美日本| 日产精品乱码卡一卡2卡三| av免费在线看不卡| 成人漫画全彩无遮挡| 日韩免费高清中文字幕av| 亚洲国产欧美日韩在线播放| 日韩精品免费视频一区二区三区 | 成年av动漫网址| 久久婷婷青草| 夫妻性生交免费视频一级片| 男的添女的下面高潮视频| 精品国产一区二区久久| 久久久久人妻精品一区果冻| 一本大道久久a久久精品| 国产爽快片一区二区三区| 人妻 亚洲 视频| 国产免费现黄频在线看| 91午夜精品亚洲一区二区三区| 最近最新中文字幕免费大全7| 国产高清国产精品国产三级| 国语对白做爰xxxⅹ性视频网站| 免费看av在线观看网站| 免费少妇av软件| 国产高清不卡午夜福利| 欧美日韩视频精品一区| xxx大片免费视频| 亚洲精品av麻豆狂野| 午夜激情av网站| 97在线人人人人妻| 欧美激情国产日韩精品一区| 我的老师免费观看完整版| 国产极品粉嫩免费观看在线 | 亚洲国产色片| 久久久久久久国产电影| 久久亚洲国产成人精品v| 精品熟女少妇av免费看| 国产熟女午夜一区二区三区 | 一个人免费看片子| 十八禁网站网址无遮挡| 成人毛片a级毛片在线播放| 日韩制服骚丝袜av| 在线观看美女被高潮喷水网站| 亚洲精品第二区| 久久久a久久爽久久v久久| 国产欧美亚洲国产| 这个男人来自地球电影免费观看 | 狂野欧美激情性xxxx在线观看| 国模一区二区三区四区视频| 丰满饥渴人妻一区二区三| 黄色欧美视频在线观看| 国产成人精品在线电影| 日本色播在线视频| 激情五月婷婷亚洲| 插阴视频在线观看视频| 国产爽快片一区二区三区| 久久久久久久久久成人| 天天操日日干夜夜撸| 99热网站在线观看| 国产高清有码在线观看视频| 久久精品夜色国产| av黄色大香蕉| 欧美最新免费一区二区三区| 免费大片黄手机在线观看| 日产精品乱码卡一卡2卡三| 国产国拍精品亚洲av在线观看| 欧美日韩国产mv在线观看视频| 亚洲av综合色区一区| 亚洲精品一二三| 国产永久视频网站| 国产成人91sexporn| 国产综合精华液| 亚洲精品第二区| 久久久精品免费免费高清| 欧美日韩综合久久久久久| 亚洲美女视频黄频| 91在线精品国自产拍蜜月| 一边摸一边做爽爽视频免费| 天天影视国产精品| 老司机影院毛片| 欧美日韩视频精品一区| 黑人猛操日本美女一级片| 美女cb高潮喷水在线观看| 新久久久久国产一级毛片| 人人妻人人澡人人爽人人夜夜| 蜜臀久久99精品久久宅男| 亚洲伊人久久精品综合| 免费观看性生交大片5| 国产亚洲一区二区精品| 国产午夜精品一二区理论片| 久久99一区二区三区| 精品久久蜜臀av无| 日韩欧美精品免费久久| 婷婷成人精品国产| 亚洲无线观看免费| 亚洲精品乱码久久久久久按摩| 中国三级夫妇交换| 中文字幕av电影在线播放| 91精品国产九色| 伦精品一区二区三区| 蜜臀久久99精品久久宅男| 97在线视频观看| 热99国产精品久久久久久7| 老女人水多毛片| 免费av不卡在线播放| 一区在线观看完整版| 亚洲色图综合在线观看| 亚洲av综合色区一区| 日韩在线高清观看一区二区三区| 91久久精品电影网| 看非洲黑人一级黄片| 欧美精品国产亚洲| xxxhd国产人妻xxx| 91在线精品国自产拍蜜月| 国产精品久久久久久精品古装| 午夜福利,免费看| 亚洲精品自拍成人| 在线观看国产h片| 国产精品一二三区在线看| 久久亚洲国产成人精品v| 亚洲精品自拍成人| 多毛熟女@视频| 亚洲五月色婷婷综合| kizo精华| 久久久久久久久久久久大奶| 免费久久久久久久精品成人欧美视频 | 精品亚洲成国产av| 色网站视频免费| 亚洲综合色惰| 一级毛片我不卡| 乱码一卡2卡4卡精品| 美女脱内裤让男人舔精品视频| 亚洲欧美成人综合另类久久久| 国模一区二区三区四区视频| 亚洲国产精品国产精品| 欧美日韩一区二区视频在线观看视频在线| 日韩不卡一区二区三区视频在线| 久久韩国三级中文字幕| 丝瓜视频免费看黄片| 国产成人精品在线电影| 我的女老师完整版在线观看| 一区二区三区乱码不卡18| 伊人久久国产一区二区| 久久精品人人爽人人爽视色| 少妇人妻久久综合中文| 色婷婷久久久亚洲欧美| 国产在线一区二区三区精| 成人亚洲欧美一区二区av| 九九久久精品国产亚洲av麻豆| 亚洲色图综合在线观看| 精品99又大又爽又粗少妇毛片| videosex国产| 日韩大片免费观看网站| 日韩亚洲欧美综合| 国产高清三级在线| 久久久国产精品麻豆| 日韩精品有码人妻一区| 卡戴珊不雅视频在线播放| 国产日韩欧美亚洲二区| 插逼视频在线观看| 一本色道久久久久久精品综合| 亚洲精品国产av成人精品| 天天影视国产精品| 亚洲精品自拍成人| 成人18禁高潮啪啪吃奶动态图 | 秋霞伦理黄片| 亚洲精品乱码久久久久久按摩| 亚洲欧洲精品一区二区精品久久久 | av.在线天堂| 久久久久视频综合| 下体分泌物呈黄色| 国产伦理片在线播放av一区| 国产精品无大码| 国产免费一区二区三区四区乱码| 日本91视频免费播放| 三级国产精品欧美在线观看| 一边亲一边摸免费视频| 久久久久国产精品人妻一区二区| 欧美+日韩+精品| 少妇人妻久久综合中文| 熟女人妻精品中文字幕| 久久婷婷青草| 69精品国产乱码久久久| 男男h啪啪无遮挡| 国产成人精品一,二区| 美女大奶头黄色视频| 多毛熟女@视频| 欧美3d第一页| 午夜精品国产一区二区电影| 高清午夜精品一区二区三区| 成人国产av品久久久| 一个人看视频在线观看www免费| 日韩不卡一区二区三区视频在线| 夜夜骑夜夜射夜夜干| 亚洲三级黄色毛片| 视频在线观看一区二区三区| 91在线精品国自产拍蜜月| 涩涩av久久男人的天堂| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 内地一区二区视频在线| 中文字幕人妻丝袜制服| 各种免费的搞黄视频| 97精品久久久久久久久久精品| 国产探花极品一区二区| 亚洲av男天堂| 九色亚洲精品在线播放| 五月开心婷婷网| 日韩三级伦理在线观看| 精品一区二区三区视频在线| 少妇被粗大的猛进出69影院 | 免费高清在线观看日韩| xxx大片免费视频| 蜜桃国产av成人99| 狠狠婷婷综合久久久久久88av| 97在线视频观看| 狠狠精品人妻久久久久久综合| 欧美日韩综合久久久久久| a级毛片免费高清观看在线播放| 久热久热在线精品观看| 久久久久久久久久人人人人人人| 色哟哟·www| 飞空精品影院首页| 国产白丝娇喘喷水9色精品| 中文精品一卡2卡3卡4更新| 韩国高清视频一区二区三区| 国产免费现黄频在线看| 日韩三级伦理在线观看| 亚洲成人一二三区av| 人体艺术视频欧美日本| 男女国产视频网站| 亚洲av不卡在线观看| 人体艺术视频欧美日本| av一本久久久久| 熟妇人妻不卡中文字幕| 三级国产精品欧美在线观看| 国产成人av激情在线播放 | 一级毛片黄色毛片免费观看视频| 国国产精品蜜臀av免费| 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 国产高清有码在线观看视频| 亚洲国产毛片av蜜桃av| 亚洲精品自拍成人| 国产成人a∨麻豆精品| 久久精品熟女亚洲av麻豆精品| 日本午夜av视频| 亚洲国产精品成人久久小说| 欧美日本中文国产一区发布| 成人国产av品久久久| 午夜福利网站1000一区二区三区| 激情五月婷婷亚洲| av线在线观看网站| 人成视频在线观看免费观看| 校园人妻丝袜中文字幕| 高清在线视频一区二区三区| 满18在线观看网站| 欧美性感艳星| 又黄又爽又刺激的免费视频.| 这个男人来自地球电影免费观看 | 国产老妇伦熟女老妇高清| 麻豆成人av视频| 女人精品久久久久毛片| av女优亚洲男人天堂| 久久久久久久久久久免费av| 国产成人精品婷婷| 狠狠精品人妻久久久久久综合| 在线看a的网站| 美女福利国产在线| 午夜视频国产福利| 最黄视频免费看| 午夜激情福利司机影院| av有码第一页| 久久人人爽人人爽人人片va| 亚洲国产欧美日韩在线播放| 国产免费福利视频在线观看| 久久久久久久精品精品| 一区二区三区精品91| 国产精品女同一区二区软件| 亚洲图色成人| 国产黄色免费在线视频| 最后的刺客免费高清国语| 97超碰精品成人国产| 亚洲av成人精品一二三区| 久久精品夜色国产| 伦理电影大哥的女人| 免费久久久久久久精品成人欧美视频 | 亚洲av国产av综合av卡| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 麻豆乱淫一区二区| av有码第一页| 亚洲精品国产色婷婷电影| 伊人久久国产一区二区| 国产精品成人在线| 亚洲国产精品成人久久小说| 成人毛片60女人毛片免费| 少妇猛男粗大的猛烈进出视频| 午夜福利,免费看| √禁漫天堂资源中文www| 高清在线视频一区二区三区| 男女啪啪激烈高潮av片| 欧美日韩视频高清一区二区三区二| 国产视频内射| av视频免费观看在线观看| 国产有黄有色有爽视频| 亚洲av中文av极速乱| 美女主播在线视频| 爱豆传媒免费全集在线观看| 欧美性感艳星| 国产精品熟女久久久久浪| 少妇的逼水好多| 久久精品国产自在天天线| 精品一区二区三区视频在线| 国产白丝娇喘喷水9色精品| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧洲国产日韩| 国产视频内射| 欧美精品一区二区大全| 丰满少妇做爰视频| 欧美+日韩+精品| 亚洲欧洲精品一区二区精品久久久 | 欧美少妇被猛烈插入视频| 91精品国产九色| 国产黄色免费在线视频| 欧美 亚洲 国产 日韩一| 黄色一级大片看看| 国国产精品蜜臀av免费| 国产一区有黄有色的免费视频| 精品久久久久久电影网| 国产无遮挡羞羞视频在线观看| 天天操日日干夜夜撸| 亚洲婷婷狠狠爱综合网| 另类亚洲欧美激情| 高清av免费在线| 老女人水多毛片| 国产欧美日韩综合在线一区二区| 国产精品一二三区在线看| 欧美精品人与动牲交sv欧美| 美女cb高潮喷水在线观看| av有码第一页| 国产成人av激情在线播放 | 午夜福利网站1000一区二区三区| 国模一区二区三区四区视频| 免费大片黄手机在线观看| 夫妻午夜视频| 久久久国产欧美日韩av| 日本wwww免费看| 18禁动态无遮挡网站| 欧美精品一区二区免费开放| 日韩 亚洲 欧美在线| 纵有疾风起免费观看全集完整版| 一本久久精品| 夫妻性生交免费视频一级片| 91久久精品国产一区二区成人| 91精品三级在线观看| 欧美另类一区| 秋霞伦理黄片| 有码 亚洲区| 亚洲三级黄色毛片| 午夜激情久久久久久久| 亚洲少妇的诱惑av| 亚洲国产精品专区欧美| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 亚洲精品中文字幕在线视频| 国产午夜精品久久久久久一区二区三区| 狂野欧美激情性bbbbbb| 免费观看在线日韩| av不卡在线播放| 插逼视频在线观看| 欧美日韩国产mv在线观看视频| 看免费成人av毛片| 日本爱情动作片www.在线观看| 成人午夜精彩视频在线观看| 大片电影免费在线观看免费| 天天影视国产精品| 亚州av有码| 少妇精品久久久久久久| 国产一区二区在线观看av| 国产深夜福利视频在线观看| 看非洲黑人一级黄片| 久热久热在线精品观看| 久久热精品热| 最近的中文字幕免费完整| 伊人亚洲综合成人网| 国产精品99久久99久久久不卡 | 欧美日韩视频精品一区| 嘟嘟电影网在线观看| 日本色播在线视频| 国产69精品久久久久777片| 如日韩欧美国产精品一区二区三区 | 熟女av电影| 国产国拍精品亚洲av在线观看| 亚洲国产毛片av蜜桃av| 亚洲一级一片aⅴ在线观看| 免费黄频网站在线观看国产| 亚洲国产av影院在线观看| 亚洲av在线观看美女高潮| 少妇丰满av| 内地一区二区视频在线| 亚洲美女视频黄频| 国产在线一区二区三区精| 国产一区二区在线观看av| 99久久人妻综合| 日韩不卡一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 99热国产这里只有精品6| 蜜桃久久精品国产亚洲av| 免费播放大片免费观看视频在线观看| 熟妇人妻不卡中文字幕| 少妇人妻久久综合中文| 免费黄频网站在线观看国产| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 精品久久久久久久久av| 久久婷婷青草| 丰满乱子伦码专区| 午夜福利,免费看| 天美传媒精品一区二区| 亚洲在久久综合| 久久 成人 亚洲| .国产精品久久| 插阴视频在线观看视频| 精品熟女少妇av免费看| 简卡轻食公司| 一级,二级,三级黄色视频| 97精品久久久久久久久久精品| 亚洲激情五月婷婷啪啪| 久久精品久久久久久久性| 成人黄色视频免费在线看| 在线亚洲精品国产二区图片欧美 | 亚洲五月色婷婷综合| 亚洲美女黄色视频免费看| 飞空精品影院首页| 国产精品秋霞免费鲁丝片| 一级爰片在线观看| 日本猛色少妇xxxxx猛交久久| 国产午夜精品一二区理论片| 欧美性感艳星| 女性生殖器流出的白浆| 午夜91福利影院| 亚洲欧美一区二区三区黑人 | 国产男女超爽视频在线观看| 丰满少妇做爰视频| 欧美xxⅹ黑人| 免费黄色在线免费观看| 国产永久视频网站| 亚洲欧美一区二区三区国产| 大话2 男鬼变身卡| 水蜜桃什么品种好| 街头女战士在线观看网站| 纯流量卡能插随身wifi吗| 九九久久精品国产亚洲av麻豆| 欧美精品高潮呻吟av久久| 黑人猛操日本美女一级片| 午夜福利,免费看| 卡戴珊不雅视频在线播放| 99精国产麻豆久久婷婷| 一本大道久久a久久精品| av网站免费在线观看视频| 国产精品99久久99久久久不卡 | 精品久久久久久电影网|