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

    稠油熱采四通管沖蝕特性影響因素?cái)?shù)值模擬研究

    2022-12-30 02:10:42朱麗云王森王國(guó)濤石佳瑞王振波劉岑凡
    表面技術(shù) 2022年11期
    關(guān)鍵詞:四通沖蝕肩部

    朱麗云,王森,王國(guó)濤,石佳瑞,王振波,劉岑凡

    稠油熱采四通管沖蝕特性影響因素?cái)?shù)值模擬研究

    朱麗云1,王森1,王國(guó)濤2,石佳瑞1,王振波1,劉岑凡3

    (1.中國(guó)石油大學(xué)(華東),山東 青島 266580;2.青島港(集團(tuán))有限公司, 山東 青島 266011;3.中國(guó)特種設(shè)備檢測(cè)研究院 特種設(shè)備安全與節(jié)能?chē)?guó)家市場(chǎng)監(jiān)管重點(diǎn)實(shí)驗(yàn)室,北京 100029)

    對(duì)稠油熱采過(guò)程中的四通管進(jìn)行氣固兩相流沖蝕研究,基于不同影響因素分析其沖蝕變化規(guī)律,并構(gòu)建四通管最大沖蝕速率的預(yù)測(cè)模型。基于氣固兩相流理論,構(gòu)建CFD-DPM-EPM數(shù)值模型,研究不同流速、顆粒質(zhì)量流量以及顆粒粒徑對(duì)四通管沖蝕的影響,預(yù)測(cè)四通管最大沖蝕速率發(fā)生位置及數(shù)值大小,并建立了關(guān)聯(lián)流速、顆粒質(zhì)量流量和顆粒粒徑的四通管沖蝕速率數(shù)學(xué)模型。在四通管的肩部位置即豎直管與水平管交匯處和豎直管的封閉端頂部?jī)商幋嬖跊_蝕,沖蝕模擬結(jié)果與實(shí)際失效四通管相吻合。隨著四通管入口流速的增加,其最大沖蝕速率呈指數(shù)形式增長(zhǎng);隨著顆粒質(zhì)量流量的增加,其最大沖蝕速率近似呈線(xiàn)性增長(zhǎng);隨著四通管內(nèi)顆粒粒徑增大,其最大沖蝕速率先減小后增大,存在沖蝕速率最小的臨界粒徑。構(gòu)建的四通管沖蝕數(shù)學(xué)模型擬合值同模擬值對(duì)比,吻合度很高。四通管出口段肩部位置沖蝕速率高于封閉端頂部位置,肩部為沖蝕磨損高危區(qū)。因此在實(shí)際應(yīng)用過(guò)程中要重點(diǎn)關(guān)注肩部的沖蝕磨損程度并及時(shí)進(jìn)行防護(hù)處理,通過(guò)局部加厚或添加擾流內(nèi)構(gòu)件來(lái)減輕肩部沖蝕;同時(shí)要適當(dāng)降低流速并盡可能減少顆粒夾帶,以增加四通管使用壽命。

    稠油熱采;四通管;沖蝕磨損;氣固兩相流;數(shù)值模擬

    近年來(lái),在稠油開(kāi)采過(guò)程中,熱采技術(shù)應(yīng)用日趨廣泛[1-3]。在稠油熱采過(guò)程的注汽階段,管內(nèi)氣體攜砂高速流動(dòng),對(duì)管道內(nèi)壁產(chǎn)生一定沖擊,其中對(duì)于四通管內(nèi)壁的沖擊尤為明顯,在連續(xù)工作一段時(shí)間后往往會(huì)引起四通管內(nèi)壁損傷甚至穿孔,很大程度上縮短了四通管道的服役時(shí)間。四通管以其操作簡(jiǎn)單、穩(wěn)定性高等優(yōu)點(diǎn)在熱采過(guò)程中得到廣泛使用[4],而氣體攜砂流動(dòng)會(huì)導(dǎo)致四通管失效泄漏,嚴(yán)重影響工業(yè)生產(chǎn)進(jìn)程,威脅生產(chǎn)人員安全[5-6]。然而,目前對(duì)稠油開(kāi)采過(guò)程中四通管道沖蝕失效問(wèn)題的研究較少,缺乏不同參數(shù)對(duì)于四通管沖蝕失效的影響研究,許多結(jié)論往往來(lái)源于工程實(shí)際,對(duì)四通管的沖蝕失效機(jī)理還缺乏相關(guān)的研究和分析。

    針對(duì)沖蝕磨損問(wèn)題各國(guó)學(xué)者已開(kāi)展了大量的研究工作,Sedrez等[7]通過(guò)CFD數(shù)值模擬對(duì)彎管進(jìn)行了研究,分析了不同流動(dòng)狀態(tài)下彎管的沖蝕變化情況。易先中等[8]采用DPM模型并引入重復(fù)性分析法,研究了不同參數(shù)對(duì)壓裂彎管的沖蝕影響規(guī)律。Ko-sinska等[9]采用數(shù)值模擬和理論分析相結(jié)合的方法,著重研究了顆粒粒徑對(duì)彎管的沖蝕影響規(guī)律。Farok-hipour等[10]考慮了粒子轉(zhuǎn)動(dòng)與粒子間的作用力,通過(guò)CFD-DEM的方法研究了在不同幾何結(jié)構(gòu)下不同流速顆粒的流動(dòng)情況。王彥驊等[11]利用數(shù)值仿真的方法,研究了氣液兩相流攜帶顆粒的情況下對(duì)π型管磨損的影響情況。曹學(xué)文等[12]通過(guò)模擬方法,分析了不同參數(shù)下彎管的失效損傷情況??梢钥吹?,目前國(guó)內(nèi)外學(xué)者對(duì)于沖蝕磨損的研究多是對(duì)于彎管的研究[13-15],對(duì)于其他復(fù)雜管道的沖蝕研究較少,而對(duì)于稠油開(kāi)采過(guò)程中常用的四通管的研究更為欠缺。作為稠油開(kāi)采過(guò)程中的關(guān)鍵管道,四通管沖蝕失效會(huì)影響到稠油開(kāi)采過(guò)程的連續(xù)性,甚至造成高溫蒸汽泄漏,引發(fā)嚴(yán)重的安全事故,因此對(duì)四通管沖蝕規(guī)律的研究顯得尤為重要。

    本文以稠油開(kāi)采過(guò)程中發(fā)生真實(shí)失效的四通管為研究對(duì)象,基于氣固兩相流理論,研究了不同流速、質(zhì)量流量以及顆粒粒徑對(duì)四通管沖蝕特性的影響,建立了關(guān)聯(lián)流速、質(zhì)量流量和顆粒粒徑的可準(zhǔn)確預(yù)測(cè)四通管沖蝕特性的數(shù)學(xué)模型,以期對(duì)四通管在使用過(guò)程中可能出現(xiàn)沖蝕磨損的位置以及沖蝕速率進(jìn)行預(yù)測(cè),為四通管的防護(hù)提供參考。

    1 計(jì)算模型

    1.1 控制方程

    連續(xù)相控制方程如式(1)—(2)所示[16]。

    式中:g表示氣相密度,kg/m3;g表示流速,m/s;表示壓力,Pa;表示氣體動(dòng)力黏度,N·s/m2;gg表示重力項(xiàng),N;F表示粒子作用于連續(xù)相的附加源項(xiàng)。

    1.2 湍流模型

    式中:G表示平均速度梯度引起的湍流動(dòng)能;b表示浮力引起的湍流動(dòng)能;M表示可壓縮湍流中波動(dòng)膨脹對(duì)總體耗散率的作用;表示湍動(dòng)能,;表示湍流耗散率,m2/s3。

    1.3 顆粒軌跡模型

    砂礫為離散相,采用DPM方法模擬,離散相的運(yùn)動(dòng)方程為[18]:

    式中:s表示顆粒速度,m/s;g表示流動(dòng)速度,m/s;s表示顆粒粒徑,m;s表示顆粒密度,kg/m3;g表示流體密度,kg/m3;(s?g)/s表示單位質(zhì)量顆粒受到的曳力,N;s表示顆粒受到的其他作用力;s表示相對(duì)雷諾數(shù);D表示曳力系數(shù)。

    1.4 沖蝕磨損模型

    四通管沖蝕速率定義為單位時(shí)間內(nèi)顆粒對(duì)四通管內(nèi)壁面造成的質(zhì)量減薄量。本文使用王國(guó)濤等[19]在Ahlert半經(jīng)驗(yàn)?zāi)P蚚20]基礎(chǔ)上通過(guò)沖蝕實(shí)驗(yàn)得到的沖蝕速率預(yù)測(cè)模型:

    式中:為沖蝕速率,kg/(m2·s);為模型常數(shù),2.67×10?7;s為顆粒形狀系數(shù);為速度指數(shù),1.74;()為沖擊角度函數(shù),表達(dá)式見(jiàn)(7)。

    1.5 幾何模型及參數(shù)

    以實(shí)際工程應(yīng)用四通管建立模型,四通管尺寸結(jié)構(gòu)如圖1所示(所標(biāo)注尺寸為流體在四通管內(nèi)部流動(dòng)區(qū)域,除入口和出口外的兩端為封閉端)。數(shù)值計(jì)算過(guò)程中,網(wǎng)格的劃分尤為重要,為保證計(jì)算效率和精度,網(wǎng)格劃分在ICEM中完成,網(wǎng)格采用六面體網(wǎng)格,并對(duì)邊界層網(wǎng)格進(jìn)行細(xì)化,網(wǎng)格如圖2所示。

    圖1四通管幾何結(jié)構(gòu)圖

    圖2 四通管網(wǎng)格劃分模型

    1.6 網(wǎng)格無(wú)關(guān)性驗(yàn)證

    由計(jì)算流體力學(xué)知識(shí)可知,網(wǎng)格劃分?jǐn)?shù)量會(huì)對(duì)流體域的流場(chǎng)結(jié)構(gòu)(如流速、壓力等)產(chǎn)生一定影響[21]。而流速的變化會(huì)導(dǎo)致流體對(duì)于顆粒的攜帶效果發(fā)生變化,進(jìn)而改變顆粒動(dòng)能,最終影響管壁沖蝕速率的大小。因此在進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證時(shí),選用四通管的最大沖蝕速率作為判斷依據(jù)。如圖3所示,分別選用6種不同網(wǎng)格數(shù)的計(jì)算模型進(jìn)行沖蝕模擬。當(dāng)網(wǎng)格數(shù)量較小時(shí),管壁最大沖蝕速率隨著網(wǎng)格數(shù)量的增加而增加;而當(dāng)網(wǎng)格數(shù)量達(dá)到30萬(wàn)后,隨著網(wǎng)格數(shù)量繼續(xù)增加,最大沖蝕率幾乎不再變化。因此,為了綜合考慮計(jì)算精度以及計(jì)算成本,采用303 584的網(wǎng)格數(shù)進(jìn)行后續(xù)研究。

    圖3 網(wǎng)格數(shù)量對(duì)最大沖蝕速率的影響

    2 四通管沖蝕機(jī)理分析

    以實(shí)際稠油熱采過(guò)程中注汽階段的實(shí)際操作參數(shù)為依據(jù)[4],氣速為10 m/s,真實(shí)砂粒密度為2 000 kg/m3,質(zhì)量流量取0.3 kg/(m2·s),砂粒粒徑取0.1 mm,顆粒球形度取0.8。四通管實(shí)際失效位置與模擬位置對(duì)比如圖4所示,從沖蝕云圖可以看到,四通管出口段肩部和封閉端頂部位置兩處存在沖蝕缺陷[19],其中四通管肩部位置沖蝕最為嚴(yán)重,這與真實(shí)生產(chǎn)過(guò)程中四通管的失效位置相吻合。

    首先通過(guò)分析四通管的壓力場(chǎng)和速度場(chǎng)來(lái)探究四通管的流場(chǎng)分布[22]。圖5為四通管的中心截面速度分布。圖6為四通管的中心截面壓力分布。由圖5可以看出,不同流速下四通管內(nèi)各位置的速度分布基本保持不變。在四通管入口一側(cè),管道內(nèi)不同位置流速基本保持相同;然而在經(jīng)過(guò)四通管肩部時(shí),由于流體流向發(fā)生改變,由豎直流入轉(zhuǎn)為水平流出,使得肩部的速度變化顯著,在入口管和出口管交匯的彎頭處流速達(dá)到最大。隨后高速流體逐漸往水平出口管的頂部轉(zhuǎn)移,造成出口管的徑向速度分布不均,且流速越大,出口管的徑向速度梯度越大。由圖6可以看出,不同流速下四通管內(nèi)各位置的壓力分布基本一致。在入口端附近時(shí),壓力分布較為均勻;而在四通管肩部壓力達(dá)到最小值,隨后流體從水平出口管流出,壓力逐漸趨于穩(wěn)定。隨著流速的增加,進(jìn)出口壓降也隨之變大。此外,在四通管的兩封閉端,流體的速度存在最小值,壓力存在最大值。

    進(jìn)一步分析顆粒的運(yùn)動(dòng)軌跡來(lái)探明顆粒對(duì)四通管內(nèi)壁碰撞造成的沖蝕磨損情況[23]。顆粒沖蝕軌跡見(jiàn)圖7,沖蝕顆粒分布見(jiàn)圖8。由圖可知,四通管沖蝕嚴(yán)重區(qū)域主要集中在肩部[24-25]。這是由于在入口直管段,砂粒在氣流推動(dòng)作用下獲得動(dòng)能向上運(yùn)動(dòng),其運(yùn)動(dòng)方向與管道軸線(xiàn)方向平行,沿管壁法向的速度分量較小,因此入口直管段基本不受沖蝕影響;而在豎直管與水平管交匯處,氣固兩相流向發(fā)生突變,砂粒在其離心力作用下被甩到壁面[19],以高流速?zèng)_擊四通管肩部,引起肩部沖蝕磨損,在長(zhǎng)時(shí)間使用后極易造成肩部位置穿孔引發(fā)泄漏。

    圖4 四通實(shí)際失效位置與模擬位置對(duì)比

    圖5 四通管中心截面速度場(chǎng)云圖

    圖6 四通管中心截面壓力場(chǎng)云圖

    圖7 顆粒沖蝕軌跡圖

    圖8 沖蝕顆粒分布圖

    3 四通管沖蝕模擬結(jié)果與分析

    為了探究流速、顆粒質(zhì)量流量以及顆粒粒徑對(duì)四通管沖蝕的影響,依據(jù)文獻(xiàn)[4]以及文獻(xiàn)[19]中的實(shí)際工況設(shè)置了如表1的變量范圍,分別對(duì)不同條件下四通管氣固兩相流沖蝕進(jìn)行模擬研究,具體的模擬結(jié)果及分析如下。

    表1 變量范圍

    Tab.1 Scope of variables

    3.1 流速對(duì)沖蝕結(jié)果的影響

    本節(jié)對(duì)不同流速(10、15、20、25、30 m/s)條件下四通管的沖蝕變化規(guī)律進(jìn)行研究,其中顆粒粒徑取0.15 mm,顆粒質(zhì)量流量分別為0.001、0.005、0.01 kg/s。圖9給出了質(zhì)量流量為0.005 kg/s時(shí)不同流速下四通管的沖蝕速率云圖,圖10給出了不同流速下四通管最大沖蝕速率的變化趨勢(shì)。

    圖9 不同流速下沖蝕速率云圖

    圖10 流速對(duì)最大沖蝕速率的影響

    由圖9沖蝕云圖可以看到,四通管的沖蝕高危區(qū)主要在豎直管與水平管交匯的肩部區(qū)域以及豎直管的封閉末端,肩部沖蝕更為嚴(yán)重。隨著流速的增大,肩部以及豎直管頂部的沖蝕程度均顯著增加,同時(shí)出口管頂部的沖蝕區(qū)域也有所增加。當(dāng)入口流速增加時(shí),流體對(duì)砂礫的攜帶效果增強(qiáng),導(dǎo)致砂礫自身動(dòng)能增加,對(duì)四通管內(nèi)壁的沖擊碰撞以及切削力度變大,造成了肩部以及豎直管封閉端的沖蝕區(qū)域和沖蝕速率的增大[26]。同時(shí)從圖5中可以看到,流速的增加會(huì)導(dǎo)致出口管出現(xiàn)較大徑向速度梯度,且速度越大,徑向梯度越顯著,流體和顆粒對(duì)出口管頂部的沖擊越大,造成了出口管頂部沖蝕程度的加深。從圖10可以看到,隨流速增加,四通管的最大沖蝕速率近似呈指數(shù)增長(zhǎng),且顆粒質(zhì)量流量越大,沖蝕速率增長(zhǎng)趨勢(shì)越顯著。顆粒質(zhì)量流量的增加導(dǎo)致四通管內(nèi)顆粒濃度增加,在高顆粒濃度和高流速共同作用下,四通管的沖蝕程度顯著增大。

    3.2 質(zhì)量流量對(duì)沖蝕結(jié)果的影響

    本節(jié)對(duì)不同顆粒質(zhì)量流量(0.001、0.003、0.005、0.008、0.01 kg/s)條件下四通管的沖蝕變化規(guī)律進(jìn)行研究,其中顆粒粒徑取0.15 mm,流速分別為10、20、30 m/s。圖11給出了流速為20 m/s時(shí)不同質(zhì)量流量下四通管的沖蝕云圖。圖12為不同顆粒質(zhì)量流量條件下四通管最大沖蝕速率的變化規(guī)律。

    圖11 不同質(zhì)量流量下沖蝕速率云圖

    圖12 質(zhì)量流量對(duì)最大沖蝕速率的影響

    由圖11中四通管沖蝕云圖可知,隨著顆粒質(zhì)量流量的增加,四通管各位置的沖蝕速率均有不同程度的增大。當(dāng)相同流速條件下,顆粒質(zhì)量流量的增大使得四通管內(nèi)顆粒濃度增大即顆粒數(shù)量增加,更多的顆粒沖擊、碰撞以及切削四通管內(nèi)壁面,使得相同時(shí)間四通管內(nèi)壁面受到的顆粒沖擊頻次增加,導(dǎo)致四通管最大沖蝕速率的增大。從圖12可以看到,隨四通管顆粒質(zhì)量流量增大,四通管最大沖蝕速率近似為線(xiàn)性增大。

    3.3 顆粒粒徑對(duì)沖蝕結(jié)果的影響

    本節(jié)對(duì)不同顆粒粒徑(0.05、0.1、0.15、0.20、0.25 mm)條件下的四通管沖蝕速率變化情況進(jìn)行研究,流速為20 m/s,質(zhì)量流量分別取0.001、0.005、0.01 kg/s。圖13是質(zhì)量流量為0.005 kg/s時(shí)不同顆粒粒徑下四通管的沖蝕速率云圖。圖14為不同顆粒粒徑條件下四通管最大沖蝕速率的變化規(guī)律。

    圖13 不同顆粒粒徑下沖蝕速率云圖

    由圖13可知,隨著顆粒粒徑增大,四通管肩部的沖蝕速率出現(xiàn)先減小再上升的變化趨勢(shì)。當(dāng)粒徑較?。?.05 mm)時(shí),砂礫對(duì)流體的跟隨性較強(qiáng),在氣流攜帶作用下造成了四通管肩部的沖蝕區(qū)域,而豎直管頂部幾乎不受顆粒沖蝕的影響。隨著顆粒粒徑的增大(0.05~0.15 mm),四通管豎直封閉末端的沖蝕程度有所增加,但其肩部沖蝕速率下降。當(dāng)顆粒粒徑為0.15 mm時(shí),四通管沖蝕速率最小。而當(dāng)顆粒粒徑繼續(xù)增加時(shí),四通管肩部的沖蝕速率又逐漸增大。圖13的沖蝕云圖表明顆粒粒徑對(duì)四通管最大沖蝕速率的影響存在臨界值(0.15 mm)[27]。結(jié)合圖14可以看到,當(dāng)四通管內(nèi)顆粒粒徑小于臨界值時(shí),最大沖蝕速率隨顆粒粒徑增加而逐漸降低;當(dāng)顆粒粒徑大于臨界值時(shí),最大沖蝕速率隨顆粒粒徑增加而逐漸增加(顆粒質(zhì)量流量較小時(shí)臨界粒徑為0.1 mm)。當(dāng)顆粒粒徑小于臨界粒徑時(shí),顆粒粒徑的增加使得管道內(nèi)的顆粒數(shù)量減少,使得在相同時(shí)間內(nèi)四通管受到的顆粒沖擊頻次減少,造成了肩部的沖蝕速率降低。隨著顆粒粒徑超過(guò)臨界粒徑,在質(zhì)量流量不變的情況下,盡管其顆粒數(shù)量變少,但由于單個(gè)顆粒的動(dòng)能顯著增加,使得四通管內(nèi)壁受到的沖擊動(dòng)能隨之增加,因此四通管最大沖蝕速率又逐漸變大。此外可以看到,隨著顆粒質(zhì)量流量的增大,四通管的沖蝕速率顯著增加。

    圖14 顆粒粒徑對(duì)最大沖蝕速率的影響

    4 四通管最大沖蝕速率預(yù)測(cè)模型

    結(jié)合氣固兩相流四通管在實(shí)際工況下的沖蝕模擬結(jié)果,通過(guò)回歸分析方法對(duì)模擬數(shù)據(jù)進(jìn)行擬合處理,提出了基于流速、質(zhì)量流量和顆粒粒徑這3個(gè)主要影響因素的最大沖蝕速率預(yù)測(cè)模型,如式(8)所示。

    式中:表示沖蝕速率,kg/(m2·s);表示流體流速,m/s;表示顆粒質(zhì)量流量,kg/s;p表示顆粒粒徑,mm;為系數(shù)項(xiàng);、、分別表示速度指數(shù)、質(zhì)量流量指數(shù)和粒徑指數(shù)。數(shù)值模擬結(jié)果如表2。根據(jù)表2沖蝕速率數(shù)值模擬結(jié)果,利用回歸分析法擬合,最終獲得四通管最大沖蝕速率預(yù)測(cè)模型如式:

    擬合結(jié)果與模擬結(jié)果對(duì)比如圖15所示,擬合值與模擬值誤差均在20%以下,二者吻合較好。

    從得到的四通管最大沖蝕速率預(yù)測(cè)模型中可以看到,最大沖蝕速率與流速的二次冪、顆粒質(zhì)量流量的一次冪成正比。因此,在實(shí)際工程應(yīng)用中,應(yīng)適當(dāng)降低流速并盡可能減少注汽過(guò)程中的顆粒夾帶,避免由于流速過(guò)大或顆粒含量過(guò)高引起四通管沖蝕程度增加。所建立的四通管最大沖蝕速率預(yù)測(cè)模型可對(duì)未知工況下四通管的實(shí)際沖蝕程度進(jìn)行預(yù)測(cè),并預(yù)先對(duì)四通管進(jìn)行防護(hù)處理,如局部加厚、添加擾流內(nèi)構(gòu)件等。此外,根據(jù)預(yù)測(cè)模型以及允許的壁厚損失可對(duì)注汽過(guò)程的操作參數(shù)進(jìn)行調(diào)節(jié),避免操作參數(shù)設(shè)置不當(dāng)引起四通管沖蝕程度增加。

    表2 沖蝕速率模擬結(jié)果

    Tab.2 Simulation results of erosion rate

    圖15 不同條件下模擬值與擬合值對(duì)比

    5 結(jié)論

    1)四通管沖蝕主要發(fā)生在四通管的肩部以及豎直管封閉末端兩處位置,其中肩部沖蝕最為嚴(yán)重,在工程應(yīng)用過(guò)程中要加強(qiáng)對(duì)四通管肩部的厚度監(jiān)測(cè)。此外,可通過(guò)對(duì)肩部進(jìn)行加厚處理,或在四通管內(nèi)部添加擾流內(nèi)構(gòu)件來(lái)減輕肩部的沖蝕磨損。

    2)數(shù)值模擬結(jié)果顯示,隨流速的增加,四通管最大沖蝕速率呈指數(shù)增大的趨勢(shì);隨質(zhì)量流量增加,顆粒數(shù)量增加,更多的顆粒沖擊四通管內(nèi)壁,使得四通管最大沖蝕速率呈線(xiàn)性增大的趨勢(shì);隨顆粒粒徑的增加,四通管的沖蝕速率出現(xiàn)先減小后增大的變化趨勢(shì)。在實(shí)際工況中,要適當(dāng)降低流速并盡可能減少顆粒夾帶,以增加四通管使用壽命。

    3)四通管最大沖蝕率預(yù)測(cè)模型是3個(gè)主要影響因素(流速、質(zhì)量流量、顆粒粒徑)的函數(shù)。該模型對(duì)工程應(yīng)用中四通管沖蝕速率的預(yù)測(cè)具有重要作用,即對(duì)于未知工況下四通管沖蝕程度的預(yù)測(cè),并根據(jù)預(yù)測(cè)結(jié)果及時(shí)對(duì)四通管進(jìn)行防護(hù)處理;同時(shí)也可以基于該模型和允許的壁厚損失對(duì)稠油熱采注汽階段的操作參數(shù)進(jìn)行合理調(diào)節(jié),避免操作參數(shù)設(shè)置不當(dāng)引起的四通管沖蝕程度增加。

    [1] 梁志艷, 王磊磊, 唐照星. 塔河油田螺桿泵配合水溶性降黏劑稠油開(kāi)采技術(shù)[J]. 油氣藏評(píng)價(jià)與開(kāi)發(fā), 2020, 10(2): 111-115.

    LIANG Zhi-yan, WANG Lei-lei, TANG Zhao-xing. Tech-nology of Heavy Oil Recovered by Screw Pump Com-bined with Water Soluble Viscosity in Tahe Oilfield[J]. Re-servoir Evaluation and Development, 2020, 10(2): 111-115.

    [2] 王小兵, 王多琦, 李森, 等. 空心抽油桿熱水循環(huán)工藝用于稠油開(kāi)采的影響因素研究[J]. 石油工程建設(shè), 2019, 45(6): 1-6.

    WANG Xiao-bing, WANG Duo-qi, LI Sen, et al. Study on Influence Factors of Using Hot Water Circulation Process in Hollow Sucker Rod for Heavy Oil Production[J]. Pe-troleum Engineering Construction, 2019, 45(6): 1-6.

    [3] 單學(xué)軍, 張士誠(chéng), 王文雄, 等. 稠油開(kāi)采中井筒溫度影響因素分析[J]. 石油勘探與開(kāi)發(fā), 2004, 31(3): 136-139.

    SHAN Xue-jun, ZHANG Shi-cheng, WANG Wen-xiong, et al. The Factors Influencing the Wellbore Temperature in Heavy Oil Production[J]. Petroleum Exploration and Development, 2004, 31(3): 136-139.

    [4] 邱福壽, 王國(guó)濤, 彭輝, 等. 含凹坑缺陷稠油熱采井口用四通管的沖蝕數(shù)值模擬研究[J]. 壓力容器, 2021, 38(2): 28-36, 54.

    QIU Fu-shou, WANG Guo-tao, PENG Hui, et al. The Erosion Numerical Simulation Analysis of the Four-Way Pipe with Pit Defect Used in the Heavy Oil Thermal Production Wellhead[J]. Pressure Vessel Technology, 2021, 38(2): 28-36, 54.

    [5] 胡炳濤. 基于流固耦合的彎管沖蝕失效數(shù)值模擬和優(yōu)化設(shè)計(jì)[D]. 徐州: 中國(guó)礦業(yè)大學(xué), 2019.

    HU Bing-tao. Numerical Simulation and Optimization De-sign of Erosion Failure of Bent-Pipe Based on Fluid- Solid Coupling[D]. Xuzhou: China University of Mining and Technology, 2019.

    [6] 廖華林, 董林, 牛繼磊, 等. 礫石充填條件下篩管堵塞與沖蝕特性試驗(yàn)[J]. 中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 2019, 43(3): 90-97.

    LIAO Hua-lin, DONG Lin, NIU Ji-lei, et al. An Experimental Study on Plugging and Erosion Failures of Sand Screen in Grave-Packing Conditions[J]. Journal of China University of Petroleum (Edition of Natural Science), 2019, 43(3): 90-97.

    [7] SEDREZ T A, SHIRAZI S A, RAJKUMAR Y R, et al. Experiments and CFD Simulations of Erosion of a 90° Elbow in Liquid-Dominated Liquid-Solid and Dispersed- Bubble-Solid Flows[J]. Wear, 2019, 426-427: 570-580.

    [8] 易先中, 彭灼, 周元華, 等. 高壓壓裂液對(duì)JY-50壓裂彎管沖蝕行為影響的數(shù)值模擬[J]. 表面技術(shù), 2019, 48(2): 144-151.

    YI Xian-zhong, PENG Zhuo, ZHOU Yuan-hua, et al. Numerical Simulation for Erosion Behavior of High- Pressure Fracturing Fluids on JY-50 Fracturing Bend Pipe [J]. Surface Technology, 2019, 48(2): 144-151.

    [9] KOSINSKA A, BALAKIN B V, KOSINSKI P. Theore-tical Analysis of Erosion in Elbows Due to Flows with Nano- and Micro-Size Particles[J]. Powder Technology, 2020, 364: 484-493.

    [10] FAROKHIPOUR A, MANSOORI Z, RASTEH A, et al. Study of Erosion Prediction of Turbulent Gas-Solid Flow in Plugged Tees via CFD-DEM[J]. Powder Technology, 2019, 352: 136-150.

    [11] 王彥驊, 吳玉國(guó), 張紹川, 等. π型管的沖刷腐蝕數(shù)值模擬[J]. 表面技術(shù), 2020, 49(12): 259-266.

    WANG Yan-hua, WU Yu-guo, ZHANG Shao-chuan, et al. Numerical Simulation of Erosion Prediction in Π-Shaped Tube[J]. Surface Technology, 2020, 49(12): 259-266.

    [12] 曹學(xué)文, 胥錕, 彭文山. 彎管液固兩相流沖蝕失效模擬分析[J]. 表面技術(shù), 2016, 45(8): 124-131.

    CAO Xue-wen, XU Kun, PENG Wen-shan. Simulation and Analysis of Liquid-Solid Two-Phase Flow Erosion Failure in Pipe Bends[J]. Surface Technology, 2016, 45(8): 124-131.

    [13] 余禮. 水力噴砂射孔時(shí)滾筒卷繞段連續(xù)油管的沖蝕研究[D]. 荊州: 長(zhǎng)江大學(xué), 2018.

    YU Li. Study on Erosion of Coiled Roller Section in Hydro-blasting Perforation[D]. Jingzhou: Yangtze University, 2018.

    [14] 張孟昀, 馬貴陽(yáng), 李存磊, 等. 彎管與盲通管沖蝕磨損對(duì)比分析研究[J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù), 2017, 13(3): 76-81.

    ZHANG Meng-yun, MA Gui-yang, LI Cun-lei, et al. Com-parative Analysis on Erosion Wear of Elbow Pipe and Blind Tube[J]. Journal of Safety Science and Technology, 2017, 13(3): 76-81.

    [15] 王森, 朱麗云, 王振波, 等. 彎頭間連接管長(zhǎng)對(duì)π形管沖蝕影響的數(shù)值模擬[J]. 油氣儲(chǔ)運(yùn), 2021, 40(11): 1285-1292.

    WANG Sen, ZHU Li-yun, WANG Zhen-bo, et al. Nume-rical Simulation on Effect of Connection Pipe Length bet-ween Elbows on Erosion of Π-Shaped Pipeline[J]. Oil & Gas Storage and Transportation, 2021, 40(11): 1285-1292.

    [16] ZHANG Hao, LI Gong, AN Xi-zhong, et al. Numerical Study on the Erosion Process of the Low Temperature Economizer Using Computational Fluid Dynamics-Discrete Particle Method[J]. Wear, 2020, 450-451: 203269.

    [17] 侍吉清. SAGD井口在線(xiàn)檢測(cè)技術(shù)與剩余壽命評(píng)價(jià)方法研究[D]. 上海: 華東理工大學(xué), 2016.

    SHI Ji-qing. Research on the Online Testing Technology and the Residual Life Assessment Method for the Stream Assisted Gravity Drainage (SAGD) Wellhead[D]. Shanghai: East China University of Science and Technology, 2016.

    [18] 錢(qián)東良. 番禺35-2海底輸氣管道沖蝕規(guī)律研究[D]. 成都: 西南石油大學(xué), 2015.

    QIAN Dong-liang. Study on Erosion Law of Panyu 35-2 Submarine Gas Pipeline[D]. Chengdu: Southwest Petroleum University, 2015.

    [19] 王國(guó)濤, 朱麗云, 劉岑凡, 等. 基于試驗(yàn)和CFD模擬的稠油熱采井口四通管沖蝕規(guī)律分析[J]. 表面技術(shù), 2021, 50(8): 247-256.

    WANG Guo-tao, ZHU Li-yun, LIU Cen-fan, et al. Analysis of Erosion Law of Four-Way Pipe in Heavy Oil Thermal Production Wellhead Based on Experiment and CFD Simulation[J]. Surface Technology, 2021, 50(8): 247-256.

    [20] ZHOU Jia-wei, LIU Yu, LIU Song-yong, et al. Effects of Particle Shape and Swirling Intensity on Elbow Erosion in Dilute-Phase Pneumatic Conveying[J]. Wear, 2017, 380- 381: 66-77.

    [21] 陶文銓. 數(shù)值傳熱學(xué)[M]. 2版. 西安: 西安交通大學(xué)出版社, 2001.

    TAO Wen-quan. Numerical heat transfer[M]. 2nd ed. Xi'an: Xi'an Jiaotong University Press, 2001.

    [22] 彭文山. 含固體顆粒多相流彎管沖蝕機(jī)理研究[D]. 東營(yíng): 中國(guó)石油大學(xué)(華東), 2017.

    PENG Wen-shan. Study on the Solid Particle Erosion Mechanism of Pipe Bend for Multiphase Flow[D]. Dong-ying: China University of Petroleum (Huadong), 2017.

    [23] 李井洋. 大口徑輸流管道沖蝕失效影響因素研究[D]. 大慶: 東北石油大學(xué), 2019.

    LI Jing-yang. Study on Influencing Factors of Erosion Failure of Large Diameter Pipeline[D]. Daqing: Northeast Petroleum University, 2019.

    [24] 鐘功祥, 胥卜軒, 謝銳, 等. 壓裂液對(duì)四通管沖蝕磨損的仿真分析[J]. 機(jī)械科學(xué)與技術(shù), 2021, 40(4): 534-541.

    ZHONG Gong-xiang, XU Bo-xuan, XIE Rui, et al. Simu-lation Analysis of Erosion Wear Induced by Hydraulic Fracturing Fluid in Cross-Pipe[J]. Mechanical Science and Technology for Aerospace Engineering, 2021, 40(4): 534-541.

    [25] 戴義明, 丁智超, 王堯, 等. 140 MPa井口壓裂四通管道沖蝕分析[J]. 裝備制造技術(shù), 2021(4): 89-91, 96.

    DAI Yi-ming, DING Zhi-chao, WANG Yao, et al. Erosion Analysis of 140 MPa Wellhead Fracturing Cross Pipe[J]. Equipment Manufacturing Technology, 2021(4): 89-91, 96.

    [26] 胡炳濤, 朱榮濤, 李超永, 等. 彎管沖蝕失效模擬研究及影響因素分析[J]. 常州大學(xué)學(xué)報(bào)(自然科學(xué)版), 2019, 31(2): 27-34.

    HU Bing-tao, ZHU Rong-tao, LI Chao-yong, et al. Simu-lation Study on Erosion Failure of Elbow and Analysis of Influencing Factors[J]. Journal of Changzhou University (Natural Science Edition), 2019, 31(2): 27-34.

    [27] 王軍成. 集輸管線(xiàn)典型管件沖蝕磨損數(shù)值模擬[D]. 撫順: 遼寧石油化工大學(xué), 2019.

    WANG Jun-cheng. Numerical Simulation of Erosion Wear of Typical Pipes in Gathering Pipeline[D]. Fushun: Liaonign Shihua University, 2019.

    Numerical Simulation of Influence Factors on Erosion Characteristics of Four-way Pipe in Heavy Oil Thermal Recovery

    1,1,2,1,1,3

    (1. China University of Petroleum (East China), Shandong Qingdao 266580, China; 2. Qingdao Port (Group) Co. Ltd., Shandong Qingdao 266011, China; 3. Key Laboratory of Special Equipment Safety and Energy-saving for State Market Regulation, China Special Equipment Inspection and Research Institute (CSEI), Beijing 100029, China)

    The gas-solid two-phase flow erosion of the four-way pipe in the process of heavy oil thermal recovery was studied, the erosion variation law was analyzed based on different influencing factors, and the prediction model of the maximum erosion rate of the four-way pipe was built. Based on the gas-solid two-phase flow theory, a numerical model of CFD- DPM-EPM was established to research the influence of different flow velocity, mass flow rate, and particle size on the four-way pipe’s erosion.

    The location and magnitude of maximum erosion rate of four-way pipe are predicted. There are erosion defects in the shoulder position of the four-way pipe, namely the intersection of vertical pipe and horizontal pipe, and the top of the closed end of the vertical pipe. The simulation results are consistent with the actual failure of the four-way pipe. Flow field structures and particles trajectories were analyzed to verify the erosion distribution on the four-way pipe. Due to the change of the pipeline structure at the intersection, the flow direction of gas-solid phase changes significantly, and the particles impact the shoulder position of the four-way pipe with high velocity under the action of centrifugal force, eventually forming a serious erosion area of the shoulder. The maximum erosion rate of the four-way pipe increases exponentially with the increase of inlet velocity. Furthermore, a continuous erosion zone was observed at the top of the outlet pipe when the velocity was higher.With the increase of velocity, the carrying effect of fluid on particles is enhanced, resulting in the increase of the kinetic energy of particles. The collision on the inner wall of the four-way pipe and the cutting force become larger. Meanwhile, the increase of velocity will lead to a large radial velocity gradient in the outlet pipe, and the greater the velocity, the more significant the radial gradient, the greater the collision of fluid and particles on the top of the outlet pipe, resulting in the erosion of the top of the outlet pipe. The maximum erosion rate increases linearly with the increase of particle mass flow rate. When the velocity remains constant, the number of particles increases with the increase of particle mass flow rate. And then the number of particle collision per unit time of four-way pipe wall increases, and finally leads to the increase of erosion rate. With the increase of particle size, the maximum erosion rate decreases first and then increases, and there is a critical particle size with the minimum erosion rate. The study shows that the erosion rate at the shoulder of the outlet section of the four-way pipe is higher than the top of the closed end, and the shoulder is a higher risk area of erosion wear. Hence, it is necessary to pay attention to the erosion wear degree of the shoulder and carry out timely protection treatment, and reduce the erosion of the shoulder by local thickening or adding turbulence internal components in the actual application. At the same time, the velocity should be appropriately reduced and the entrainment of particles should be minimized to increase the service life of the four-way pipe.

    A mathematical model for predicting erosion rate of the four-way pipe was established. The parameters of velocity, particle mass flow rate and particle size were considered in this model. The fitted value of the mathematical model of four-way pipe was compared with the simulated value. The results show that the fitting value is in good agreement with the simulated value. The prediction model of erosion rate can be used to predict erosion degree of the four-way pipe and adjust operation parameters in the process of steam injection in heavy oil thermal recovery.

    heavy oil thermal recovery; four-way pipe; erosion wear; gas solid two-phase flow; numerical simulation

    TH117.1

    A

    1001-3660(2022)11-0244-09

    10.16490/j.cnki.issn.1001-3660.2022.11.022

    2021–10–12;

    2022–03–07

    2021-10-12;

    2022-03-07

    朱麗云(1987—),女,博士,副教授,主要研究方向?yàn)槭突ぴO(shè)備。

    ZHU Li-yun (1987-), Female, Doctor, Associate professor, Research focus: petrochemical equipment.

    朱麗云, 王森, 王國(guó)濤, 等.稠油熱采四通管沖蝕特性影響因素?cái)?shù)值模擬研究[J]. 表面技術(shù), 2022, 51(11): 244-252.

    ZHU Li-yun, WANG Sen, WANG Guo-tao, et al. Numerical Simulation of Influence Factors on Erosion Characteristics of Four-way Pipe in Heavy Oil Thermal Recovery[J]. Surface Technology, 2022, 51(11): 244-252.

    責(zé)任編輯:萬(wàn)長(zhǎng)清

    猜你喜歡
    四通沖蝕肩部
    140MPa井口壓裂四通管道沖蝕分析
    四個(gè)小竅門(mén)輕松緩解肩部疲勞
    幸福家庭(2018年2期)2018-02-07 00:23:59
    四個(gè)小竅門(mén)輕松緩解肩部疲勞
    幸福家庭(2017年2期)2017-04-18 10:13:18
    輸氣管道砂沖蝕的模擬實(shí)驗(yàn)
    環(huán)氧樹(shù)脂及其復(fù)合材料的固體顆粒沖蝕磨損
    空調(diào)四通換向閥的檢測(cè)與更換
    兒童乒乓球啟蒙訓(xùn)練中肩部損傷與預(yù)防
    肩部巨大包塊同時(shí)伴毛鞘癌及鱗癌1例
    采用三換熱器和四通閥的兩種車(chē)用熱泵系統(tǒng)的對(duì)比研究
    三種不銹鋼材料抗固體顆粒沖蝕性能研究
    高清午夜精品一区二区三区| 18禁在线播放成人免费| 国产黄色视频一区二区在线观看| 精品久久久噜噜| 免费大片黄手机在线观看| 蜜臀久久99精品久久宅男| 久久久久视频综合| 偷拍熟女少妇极品色| 亚洲精品国产成人久久av| 美女国产视频在线观看| 日产精品乱码卡一卡2卡三| 国产欧美日韩精品一区二区| 老熟女久久久| 国产成人freesex在线| 纵有疾风起免费观看全集完整版| 久久久国产一区二区| 一本久久精品| av在线app专区| 国产成人精品福利久久| 观看av在线不卡| 99九九线精品视频在线观看视频| 欧美3d第一页| 久久鲁丝午夜福利片| 国产亚洲5aaaaa淫片| 日本猛色少妇xxxxx猛交久久| 欧美 亚洲 国产 日韩一| 中文字幕亚洲精品专区| 亚洲av二区三区四区| 精华霜和精华液先用哪个| av天堂久久9| 狂野欧美白嫩少妇大欣赏| 亚洲国产精品国产精品| 久久久久久久国产电影| 久久99热6这里只有精品| 高清毛片免费看| 亚洲第一av免费看| 午夜激情福利司机影院| 一区二区三区乱码不卡18| 国产男人的电影天堂91| 午夜激情久久久久久久| 亚洲欧洲精品一区二区精品久久久 | 亚洲久久久国产精品| 亚洲国产精品国产精品| tube8黄色片| 在线观看免费视频网站a站| a级片在线免费高清观看视频| 欧美丝袜亚洲另类| 中文字幕精品免费在线观看视频 | 成人免费观看视频高清| 精品一品国产午夜福利视频| 不卡视频在线观看欧美| 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| 免费看日本二区| 下体分泌物呈黄色| 男人和女人高潮做爰伦理| 亚洲婷婷狠狠爱综合网| 91久久精品国产一区二区三区| 亚洲av成人精品一二三区| 亚洲欧洲国产日韩| 亚洲精品日韩av片在线观看| 中文字幕人妻丝袜制服| 熟女电影av网| 久久亚洲国产成人精品v| 男人添女人高潮全过程视频| 夜夜骑夜夜射夜夜干| 极品人妻少妇av视频| 免费人妻精品一区二区三区视频| 免费观看a级毛片全部| 午夜日本视频在线| 99视频精品全部免费 在线| 三级经典国产精品| 久久久久久久久久成人| 黄色欧美视频在线观看| 丝瓜视频免费看黄片| a 毛片基地| 啦啦啦视频在线资源免费观看| 2018国产大陆天天弄谢| 在线观看免费日韩欧美大片 | 免费看光身美女| 中文资源天堂在线| 国产毛片在线视频| 看十八女毛片水多多多| 国产男女超爽视频在线观看| 免费观看在线日韩| 亚洲经典国产精华液单| 国产黄片美女视频| 国产欧美亚洲国产| 黄色配什么色好看| 国产成人午夜福利电影在线观看| 人人妻人人澡人人爽人人夜夜| 大话2 男鬼变身卡| 亚洲电影在线观看av| 久久久久久久久久人人人人人人| 91久久精品电影网| 亚洲精品自拍成人| 2022亚洲国产成人精品| 免费高清在线观看视频在线观看| 久久精品国产自在天天线| 中文天堂在线官网| 欧美日韩一区二区视频在线观看视频在线| 国产精品.久久久| 国产精品久久久久久久电影| 熟女电影av网| 你懂的网址亚洲精品在线观看| 亚洲欧美成人精品一区二区| 国产毛片在线视频| 麻豆成人av视频| 欧美+日韩+精品| 日本免费在线观看一区| 亚洲精品一二三| 纯流量卡能插随身wifi吗| 久久国产亚洲av麻豆专区| 日韩三级伦理在线观看| 久久人人爽人人爽人人片va| 欧美另类一区| 久久久久久久精品精品| 一二三四中文在线观看免费高清| 一本大道久久a久久精品| 国产精品一区www在线观看| 国内精品宾馆在线| 黄色视频在线播放观看不卡| 国语对白做爰xxxⅹ性视频网站| 国内精品宾馆在线| 在线观看www视频免费| 高清在线视频一区二区三区| 精品酒店卫生间| 永久免费av网站大全| 自拍偷自拍亚洲精品老妇| av播播在线观看一区| 99久久精品热视频| 午夜av观看不卡| 日本与韩国留学比较| 国产一区二区三区av在线| 爱豆传媒免费全集在线观看| 欧美性感艳星| 高清欧美精品videossex| 高清av免费在线| 精品久久国产蜜桃| 一级片'在线观看视频| 国产黄色免费在线视频| 日韩欧美 国产精品| 国产熟女欧美一区二区| 91久久精品电影网| 永久网站在线| 人体艺术视频欧美日本| 午夜福利在线观看免费完整高清在| 伊人亚洲综合成人网| av免费观看日本| 久久久国产一区二区| 18禁在线无遮挡免费观看视频| 有码 亚洲区| 免费久久久久久久精品成人欧美视频 | 精品久久久久久电影网| 嫩草影院新地址| 少妇熟女欧美另类| 一个人免费看片子| av天堂中文字幕网| 26uuu在线亚洲综合色| 一级a做视频免费观看| 亚洲综合精品二区| 麻豆精品久久久久久蜜桃| 黄色欧美视频在线观看| 欧美激情极品国产一区二区三区 | 色婷婷久久久亚洲欧美| 亚洲精品国产av蜜桃| 久久亚洲国产成人精品v| 男女免费视频国产| 国产精品国产三级专区第一集| 日本vs欧美在线观看视频 | 国产免费又黄又爽又色| 欧美老熟妇乱子伦牲交| 天天操日日干夜夜撸| 最黄视频免费看| 深夜a级毛片| 亚洲精品色激情综合| 日本av手机在线免费观看| 性色av一级| 99久久精品一区二区三区| 国产精品不卡视频一区二区| 三级国产精品欧美在线观看| 亚洲国产精品专区欧美| 成人影院久久| a 毛片基地| 亚洲国产毛片av蜜桃av| 黄色日韩在线| 久久精品国产自在天天线| 精品久久久久久久久av| 人妻人人澡人人爽人人| 18禁动态无遮挡网站| 欧美97在线视频| 免费大片黄手机在线观看| 寂寞人妻少妇视频99o| 精华霜和精华液先用哪个| 精品一区在线观看国产| 亚洲国产最新在线播放| 欧美日韩亚洲高清精品| 亚洲av成人精品一区久久| 国产成人freesex在线| 国产精品一区二区性色av| 国产成人免费观看mmmm| 一区二区三区免费毛片| 婷婷色av中文字幕| 少妇被粗大的猛进出69影院 | 只有这里有精品99| 建设人人有责人人尽责人人享有的| 亚洲av在线观看美女高潮| 午夜免费观看性视频| 肉色欧美久久久久久久蜜桃| 夫妻午夜视频| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 我的女老师完整版在线观看| 欧美丝袜亚洲另类| 精品酒店卫生间| 亚洲自偷自拍三级| 国产精品伦人一区二区| 美女脱内裤让男人舔精品视频| 免费不卡的大黄色大毛片视频在线观看| 下体分泌物呈黄色| 大香蕉97超碰在线| 久久久久人妻精品一区果冻| 亚洲,欧美,日韩| 曰老女人黄片| 一级a做视频免费观看| 九九久久精品国产亚洲av麻豆| 少妇人妻一区二区三区视频| freevideosex欧美| 天天躁夜夜躁狠狠久久av| 亚洲电影在线观看av| 日本黄大片高清| 精品人妻熟女毛片av久久网站| 国产视频内射| 成人亚洲欧美一区二区av| 一边亲一边摸免费视频| 亚洲美女搞黄在线观看| 汤姆久久久久久久影院中文字幕| 最近手机中文字幕大全| 免费看光身美女| 日韩三级伦理在线观看| 亚洲美女搞黄在线观看| 国产精品无大码| 男男h啪啪无遮挡| 一级毛片aaaaaa免费看小| 欧美日韩视频精品一区| 人妻少妇偷人精品九色| 久久久久久久久久人人人人人人| 久久综合国产亚洲精品| 新久久久久国产一级毛片| 男人舔奶头视频| 91精品伊人久久大香线蕉| 日韩精品免费视频一区二区三区 | 久久久精品免费免费高清| 国产一区二区在线观看日韩| 97精品久久久久久久久久精品| 国产精品.久久久| 亚洲内射少妇av| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 免费看不卡的av| 精品一区二区免费观看| 亚洲成人一二三区av| 寂寞人妻少妇视频99o| 色视频在线一区二区三区| 中文字幕久久专区| 国产一区二区在线观看av| 亚洲精品第二区| 老司机影院毛片| 久久ye,这里只有精品| 丁香六月天网| 又大又黄又爽视频免费| 久久精品久久久久久久性| 永久网站在线| 18禁在线无遮挡免费观看视频| 国产精品国产av在线观看| 极品人妻少妇av视频| 麻豆精品久久久久久蜜桃| 插逼视频在线观看| 人人澡人人妻人| 日本欧美视频一区| 国产一区二区在线观看av| 国产亚洲一区二区精品| 91久久精品电影网| 日本欧美国产在线视频| 欧美区成人在线视频| 亚洲情色 制服丝袜| 久久久久久久亚洲中文字幕| 亚洲美女视频黄频| 婷婷色av中文字幕| 亚洲一区二区三区欧美精品| 国产精品久久久久久精品电影小说| 女人久久www免费人成看片| 久久久精品94久久精品| 欧美精品高潮呻吟av久久| av免费观看日本| 欧美精品亚洲一区二区| 永久网站在线| 最近手机中文字幕大全| 国产高清不卡午夜福利| 美女福利国产在线| 伦精品一区二区三区| 国产亚洲午夜精品一区二区久久| 少妇精品久久久久久久| 国产一区二区在线观看av| 热re99久久精品国产66热6| 寂寞人妻少妇视频99o| 99久久综合免费| 亚洲性久久影院| 国产亚洲91精品色在线| 国产一区有黄有色的免费视频| 26uuu在线亚洲综合色| 午夜激情福利司机影院| 中文天堂在线官网| 中文乱码字字幕精品一区二区三区| 国产精品.久久久| 欧美性感艳星| 久久毛片免费看一区二区三区| 亚洲一级一片aⅴ在线观看| 如何舔出高潮| av播播在线观看一区| 亚洲国产色片| 我的老师免费观看完整版| 六月丁香七月| 大香蕉久久网| 久久99热这里只频精品6学生| 久热这里只有精品99| 色吧在线观看| 国产在线一区二区三区精| 婷婷色综合大香蕉| 色视频www国产| 91午夜精品亚洲一区二区三区| 夜夜骑夜夜射夜夜干| 内射极品少妇av片p| 欧美国产精品一级二级三级 | 99九九线精品视频在线观看视频| kizo精华| 91久久精品国产一区二区三区| 欧美精品高潮呻吟av久久| 国产精品秋霞免费鲁丝片| 晚上一个人看的免费电影| 男人舔奶头视频| 欧美日韩av久久| 亚洲国产精品一区二区三区在线| 一区二区三区免费毛片| 成人免费观看视频高清| 黄色怎么调成土黄色| 国产精品国产av在线观看| 99久久精品热视频| 国产黄色视频一区二区在线观看| 中文字幕人妻丝袜制服| 国产高清有码在线观看视频| 精品久久久噜噜| 精品久久久久久久久亚洲| 99热这里只有是精品在线观看| 交换朋友夫妻互换小说| 亚洲精品日本国产第一区| 99精国产麻豆久久婷婷| videossex国产| 亚洲中文av在线| av国产精品久久久久影院| 日日摸夜夜添夜夜添av毛片| 噜噜噜噜噜久久久久久91| 国产淫语在线视频| 一边亲一边摸免费视频| 午夜免费观看性视频| 大片电影免费在线观看免费| 中文字幕制服av| 黄片无遮挡物在线观看| 久久精品久久久久久久性| 国产极品天堂在线| 久久精品国产亚洲网站| 久久99蜜桃精品久久| 亚洲欧洲精品一区二区精品久久久 | 免费播放大片免费观看视频在线观看| 国产成人精品婷婷| 国产免费一区二区三区四区乱码| 一级av片app| 成年女人在线观看亚洲视频| 黄色日韩在线| 欧美日韩在线观看h| 夫妻午夜视频| 国产精品久久久久成人av| 性色av一级| 亚洲欧美成人精品一区二区| 成人漫画全彩无遮挡| 免费人成在线观看视频色| 大码成人一级视频| 国产色婷婷99| 欧美另类一区| 国产精品人妻久久久久久| 亚洲精品国产色婷婷电影| 另类亚洲欧美激情| av天堂久久9| 亚洲欧美成人精品一区二区| 综合色丁香网| 五月开心婷婷网| 91在线精品国自产拍蜜月| 国产在线视频一区二区| 男男h啪啪无遮挡| 亚洲中文av在线| 毛片一级片免费看久久久久| 美女脱内裤让男人舔精品视频| 能在线免费看毛片的网站| 国产男女超爽视频在线观看| 精品久久久精品久久久| 夜夜看夜夜爽夜夜摸| 美女cb高潮喷水在线观看| 国产探花极品一区二区| 日本黄大片高清| 深夜a级毛片| 一区二区av电影网| 日日爽夜夜爽网站| 自线自在国产av| 欧美三级亚洲精品| 亚洲综合色惰| 天美传媒精品一区二区| 丝瓜视频免费看黄片| 国产乱人偷精品视频| 成人亚洲精品一区在线观看| 久久久久精品久久久久真实原创| 午夜福利影视在线免费观看| 涩涩av久久男人的天堂| 99国产精品免费福利视频| 黑人高潮一二区| 精品一区二区三区视频在线| 日本黄色日本黄色录像| 少妇的逼好多水| 精品亚洲乱码少妇综合久久| 少妇人妻一区二区三区视频| 久久精品夜色国产| 夜夜看夜夜爽夜夜摸| 91精品伊人久久大香线蕉| 国产成人freesex在线| 成人18禁高潮啪啪吃奶动态图 | 我的老师免费观看完整版| 中文字幕av电影在线播放| 大又大粗又爽又黄少妇毛片口| av播播在线观看一区| 男人和女人高潮做爰伦理| 精品酒店卫生间| 免费观看性生交大片5| 国产欧美亚洲国产| 色视频在线一区二区三区| a级片在线免费高清观看视频| 五月玫瑰六月丁香| 国产欧美日韩一区二区三区在线 | av免费观看日本| 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 免费大片18禁| 久久婷婷青草| 美女国产视频在线观看| 特大巨黑吊av在线直播| 日产精品乱码卡一卡2卡三| 亚洲欧美一区二区三区国产| 永久网站在线| 亚洲欧美清纯卡通| 黄色日韩在线| 久久6这里有精品| 国产伦理片在线播放av一区| 国产精品免费大片| 免费黄频网站在线观看国产| 国产欧美日韩综合在线一区二区 | 五月开心婷婷网| 最近中文字幕高清免费大全6| 亚洲无线观看免费| 男男h啪啪无遮挡| 亚洲国产精品一区二区三区在线| 青春草国产在线视频| 五月伊人婷婷丁香| 在线观看一区二区三区激情| 亚洲精品国产色婷婷电影| 国产欧美亚洲国产| 色婷婷久久久亚洲欧美| 天堂俺去俺来也www色官网| 哪个播放器可以免费观看大片| 精品亚洲乱码少妇综合久久| 欧美国产精品一级二级三级 | 婷婷色综合www| 亚洲精品自拍成人| 自线自在国产av| 久久人人爽人人爽人人片va| 99国产精品免费福利视频| videossex国产| 国产永久视频网站| 香蕉精品网在线| 亚洲综合色惰| 欧美高清成人免费视频www| 久久 成人 亚洲| 久久99一区二区三区| 日韩一区二区三区影片| www.色视频.com| h视频一区二区三区| 久久综合国产亚洲精品| 在线天堂最新版资源| 秋霞在线观看毛片| 亚洲伊人久久精品综合| 久久这里有精品视频免费| 视频中文字幕在线观看| 精品国产露脸久久av麻豆| 国产精品99久久99久久久不卡 | 国产免费视频播放在线视频| 边亲边吃奶的免费视频| 久久久久久久大尺度免费视频| 国产一区亚洲一区在线观看| 在线播放无遮挡| 久久午夜福利片| 亚洲av欧美aⅴ国产| 18禁动态无遮挡网站| 欧美精品国产亚洲| 亚洲国产日韩一区二区| 国产片特级美女逼逼视频| 街头女战士在线观看网站| 久久精品国产鲁丝片午夜精品| 2018国产大陆天天弄谢| 精品酒店卫生间| 久久久久精品性色| 国产 一区精品| 精品亚洲成国产av| 久久人人爽av亚洲精品天堂| 欧美日韩一区二区视频在线观看视频在线| 午夜激情福利司机影院| 国产成人91sexporn| 80岁老熟妇乱子伦牲交| 卡戴珊不雅视频在线播放| av网站免费在线观看视频| 91精品国产九色| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品久久久久成人av| 在线观看免费日韩欧美大片 | 少妇被粗大的猛进出69影院 | 亚洲国产精品专区欧美| 久久精品久久久久久久性| 亚洲精品日韩av片在线观看| 丁香六月天网| 激情五月婷婷亚洲| 久久女婷五月综合色啪小说| 亚洲精品成人av观看孕妇| 精品国产一区二区久久| 人人妻人人看人人澡| 99九九在线精品视频 | 免费黄网站久久成人精品| 一级片'在线观看视频| 九草在线视频观看| 国产片特级美女逼逼视频| 少妇丰满av| 久久婷婷青草| 大片免费播放器 马上看| 久久久亚洲精品成人影院| av线在线观看网站| av一本久久久久| 国产精品一区二区在线观看99| a 毛片基地| 亚洲怡红院男人天堂| 亚洲欧美成人精品一区二区| 亚洲欧美清纯卡通| 久久国产乱子免费精品| 黄色欧美视频在线观看| 一二三四中文在线观看免费高清| 亚洲激情五月婷婷啪啪| 国产亚洲av片在线观看秒播厂| 我要看日韩黄色一级片| 欧美区成人在线视频| 免费看av在线观看网站| 亚洲自偷自拍三级| av有码第一页| 少妇熟女欧美另类| 免费不卡的大黄色大毛片视频在线观看| 国产探花极品一区二区| 七月丁香在线播放| 国产69精品久久久久777片| 在线观看av片永久免费下载| 在线观看www视频免费| 久久国产乱子免费精品| 亚洲国产欧美日韩在线播放 | 精品少妇久久久久久888优播| 久久久久久久久久久免费av| 午夜老司机福利剧场| 寂寞人妻少妇视频99o| 久久精品国产亚洲网站| 一区在线观看完整版| tube8黄色片| 久久久久人妻精品一区果冻| 亚洲成色77777| 国产成人aa在线观看| av免费观看日本| 成人美女网站在线观看视频| 成人漫画全彩无遮挡| 伦理电影大哥的女人| 亚洲欧美清纯卡通| 久久热精品热| 超碰97精品在线观看| av免费在线看不卡| 日韩在线高清观看一区二区三区| 嫩草影院入口| a级毛片在线看网站| 如何舔出高潮| 久久久久视频综合| 久久久欧美国产精品| 久久 成人 亚洲| 国产女主播在线喷水免费视频网站| 熟女人妻精品中文字幕| 日日啪夜夜撸| 亚洲国产成人一精品久久久| 一级毛片 在线播放| kizo精华| 一区二区三区乱码不卡18| 亚洲精品自拍成人| 亚洲精品乱码久久久v下载方式| 97超视频在线观看视频| 久久国产精品男人的天堂亚洲 | 熟女电影av网| 久久精品夜色国产| 久久精品国产自在天天线|