• <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ì)比研究
    三種不銹鋼材料抗固體顆粒沖蝕性能研究
    伦精品一区二区三区| 性色avwww在线观看| 国产极品粉嫩免费观看在线 | 一本一本综合久久| 午夜久久久在线观看| 婷婷色麻豆天堂久久| 久久亚洲国产成人精品v| 亚洲精品国产av蜜桃| 一本久久精品| 99九九线精品视频在线观看视频| 免费人妻精品一区二区三区视频| 久久精品国产亚洲av涩爱| 欧美成人精品欧美一级黄| 最后的刺客免费高清国语| 国产亚洲午夜精品一区二区久久| 国产精品福利在线免费观看| 国产高清不卡午夜福利| 久久久精品免费免费高清| 久热久热在线精品观看| 性色avwww在线观看| 伦理电影免费视频| 亚洲欧洲国产日韩| 寂寞人妻少妇视频99o| 人妻一区二区av| 精品亚洲乱码少妇综合久久| 熟女人妻精品中文字幕| 高清不卡的av网站| 国产精品免费大片| av播播在线观看一区| .国产精品久久| 欧美变态另类bdsm刘玥| 中文字幕免费在线视频6| 日韩免费高清中文字幕av| √禁漫天堂资源中文www| 欧美日韩视频精品一区| 肉色欧美久久久久久久蜜桃| 蜜桃在线观看..| 99热这里只有精品一区| 日本欧美国产在线视频| 久久久精品94久久精品| 最新中文字幕久久久久| 亚洲国产精品一区二区三区在线| 欧美+日韩+精品| 在线观看人妻少妇| 色视频www国产| 春色校园在线视频观看| 青春草亚洲视频在线观看| 亚洲国产精品专区欧美| 最新的欧美精品一区二区| 大陆偷拍与自拍| 18+在线观看网站| 涩涩av久久男人的天堂| 亚洲精品成人av观看孕妇| 午夜福利视频精品| 免费av中文字幕在线| 国产精品秋霞免费鲁丝片| 男人舔奶头视频| 熟女电影av网| 99re6热这里在线精品视频| 少妇精品久久久久久久| 色视频www国产| 国产精品嫩草影院av在线观看| 91在线精品国自产拍蜜月| 久久久国产一区二区| 成人特级av手机在线观看| 久久午夜综合久久蜜桃| 欧美国产精品一级二级三级 | 午夜福利影视在线免费观看| 成人18禁高潮啪啪吃奶动态图 | 女性生殖器流出的白浆| 人人妻人人看人人澡| 精品亚洲乱码少妇综合久久| 三级国产精品片| 亚洲四区av| 91久久精品国产一区二区成人| 成人18禁高潮啪啪吃奶动态图 | 免费少妇av软件| 久久久久国产精品人妻一区二区| 王馨瑶露胸无遮挡在线观看| 成人特级av手机在线观看| 七月丁香在线播放| 国内揄拍国产精品人妻在线| 成人毛片a级毛片在线播放| 成年人免费黄色播放视频 | 波野结衣二区三区在线| 色5月婷婷丁香| 国产精品麻豆人妻色哟哟久久| 精品99又大又爽又粗少妇毛片| 亚洲精品第二区| 乱系列少妇在线播放| 国产在视频线精品| 亚洲色图综合在线观看| 视频中文字幕在线观看| 欧美另类一区| 一本色道久久久久久精品综合| 日本色播在线视频| 中文天堂在线官网| 熟女电影av网| 女人久久www免费人成看片| 伦精品一区二区三区| 国产精品人妻久久久久久| 亚洲第一av免费看| 精品一品国产午夜福利视频| av视频免费观看在线观看| 日本91视频免费播放| 亚洲人成网站在线播| 国产无遮挡羞羞视频在线观看| 国产黄色视频一区二区在线观看| 成人无遮挡网站| 国产成人免费无遮挡视频| 一级毛片久久久久久久久女| 狂野欧美激情性xxxx在线观看| 国产精品三级大全| 97超视频在线观看视频| 久久久久精品久久久久真实原创| 欧美精品一区二区免费开放| 午夜福利网站1000一区二区三区| 成人亚洲欧美一区二区av| 亚洲av二区三区四区| 中国美白少妇内射xxxbb| 国产av一区二区精品久久| 日本欧美视频一区| av线在线观看网站| 中文字幕制服av| 一级av片app| 亚洲精品中文字幕在线视频 | 国产精品麻豆人妻色哟哟久久| 国产在线一区二区三区精| 熟妇人妻不卡中文字幕| 纯流量卡能插随身wifi吗| 超碰97精品在线观看| 熟女人妻精品中文字幕| 中文乱码字字幕精品一区二区三区| 午夜福利影视在线免费观看| 久久精品久久久久久久性| 精品人妻熟女av久视频| 曰老女人黄片| 国产69精品久久久久777片| 亚洲图色成人| 日韩,欧美,国产一区二区三区| 久久久亚洲精品成人影院| 欧美日韩视频精品一区| 菩萨蛮人人尽说江南好唐韦庄| 精品国产一区二区久久| 国产 一区精品| 亚洲欧美精品专区久久| 国产午夜精品久久久久久一区二区三区| 天天躁夜夜躁狠狠久久av| 纵有疾风起免费观看全集完整版| 女人精品久久久久毛片| 欧美日韩视频高清一区二区三区二| 亚洲综合精品二区| 亚洲欧美中文字幕日韩二区| 日韩亚洲欧美综合| 久久午夜福利片| 能在线免费看毛片的网站| 亚洲va在线va天堂va国产| 乱系列少妇在线播放| 秋霞在线观看毛片| 久久人人爽av亚洲精品天堂| 午夜免费鲁丝| 久久99热这里只频精品6学生| 久久综合国产亚洲精品| 日韩成人av中文字幕在线观看| 久久久久久久精品精品| 国产黄片视频在线免费观看| av国产精品久久久久影院| 免费观看无遮挡的男女| 国产精品久久久久久av不卡| 久久国产亚洲av麻豆专区| 国产精品久久久久久av不卡| 亚洲综合精品二区| 国产精品熟女久久久久浪| 我要看黄色一级片免费的| 男人爽女人下面视频在线观看| 麻豆成人午夜福利视频| 韩国av在线不卡| 性色av一级| 国产精品伦人一区二区| 亚洲av在线观看美女高潮| 中文在线观看免费www的网站| 国产精品麻豆人妻色哟哟久久| 亚洲丝袜综合中文字幕| 麻豆成人av视频| 黄色一级大片看看| 黄色一级大片看看| 国产熟女欧美一区二区| 久久热精品热| 亚洲精品乱久久久久久| 国国产精品蜜臀av免费| 丰满迷人的少妇在线观看| 日日摸夜夜添夜夜爱| 五月天丁香电影| 三级国产精品片| 国产91av在线免费观看| tube8黄色片| 少妇的逼水好多| www.av在线官网国产| 在线天堂最新版资源| 国产极品天堂在线| 日韩伦理黄色片| av又黄又爽大尺度在线免费看| 伦理电影免费视频| 中文字幕av电影在线播放| 欧美 亚洲 国产 日韩一| 麻豆精品久久久久久蜜桃| 久久精品国产a三级三级三级| 午夜视频国产福利| 女性生殖器流出的白浆| 天天躁夜夜躁狠狠久久av| 日本色播在线视频| 一级片'在线观看视频| 久久久久人妻精品一区果冻| 一二三四中文在线观看免费高清| 99久久精品国产国产毛片| 国产精品一二三区在线看| 国产日韩欧美亚洲二区| 人人妻人人看人人澡| .国产精品久久| 国产免费一级a男人的天堂| 亚洲精品视频女| 国产成人免费观看mmmm| 夫妻午夜视频| 大片免费播放器 马上看| 在线 av 中文字幕| 在线观看www视频免费| 99热这里只有精品一区| 青春草视频在线免费观看| 国产一区二区三区综合在线观看 | 亚洲国产精品一区三区| 国精品久久久久久国模美| 亚洲国产精品成人久久小说| 午夜影院在线不卡| 夫妻性生交免费视频一级片| av免费在线看不卡| 久久青草综合色| 熟女人妻精品中文字幕| 丝瓜视频免费看黄片| 国产成人午夜福利电影在线观看| 大陆偷拍与自拍| 中文欧美无线码| 亚洲av福利一区| 亚洲成人手机| 免费在线观看成人毛片| 亚洲第一av免费看| 女性生殖器流出的白浆| 久久久午夜欧美精品| 精品久久久久久电影网| 国产精品三级大全| 一本大道久久a久久精品| 国产乱来视频区| 亚洲国产精品国产精品| 两个人免费观看高清视频 | 亚洲av福利一区| 久久青草综合色| av女优亚洲男人天堂| 久久久久久久大尺度免费视频| 日韩精品免费视频一区二区三区 | 亚洲精品日韩av片在线观看| av又黄又爽大尺度在线免费看| 国产黄片视频在线免费观看| 亚洲四区av| 欧美一级a爱片免费观看看| 欧美日韩综合久久久久久| 91精品国产九色| 少妇精品久久久久久久| 久久久久久久久久成人| 欧美xxⅹ黑人| 少妇被粗大猛烈的视频| 久久久国产一区二区| 久久人人爽av亚洲精品天堂| 在现免费观看毛片| 午夜精品国产一区二区电影| 日韩电影二区| 国产男女超爽视频在线观看| 国内少妇人妻偷人精品xxx网站| 国产伦精品一区二区三区视频9| 激情五月婷婷亚洲| 日韩欧美精品免费久久| 亚洲av综合色区一区| 国产淫语在线视频| 伦理电影大哥的女人| 成年美女黄网站色视频大全免费 | 亚洲欧美成人精品一区二区| 黄色一级大片看看| 色5月婷婷丁香| 日韩熟女老妇一区二区性免费视频| 国产毛片在线视频| 中文天堂在线官网| 亚洲一级一片aⅴ在线观看| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女av久视频| 高清在线视频一区二区三区| 精品少妇黑人巨大在线播放| 美女大奶头黄色视频| 大片电影免费在线观看免费| 久久国产亚洲av麻豆专区| 色视频www国产| 人妻 亚洲 视频| 制服丝袜香蕉在线| 亚洲欧美成人综合另类久久久| 成人影院久久| 国产美女午夜福利| 国产伦精品一区二区三区四那| 日韩欧美 国产精品| 女人精品久久久久毛片| 一二三四中文在线观看免费高清| 天堂俺去俺来也www色官网| av专区在线播放| 一级,二级,三级黄色视频| 精品午夜福利在线看| 一区二区三区免费毛片| 91精品国产九色| 尾随美女入室| av播播在线观看一区| 久久午夜综合久久蜜桃| 亚洲欧美中文字幕日韩二区| 日本欧美视频一区| 高清午夜精品一区二区三区| 视频区图区小说| 亚洲欧美日韩卡通动漫| 国产精品久久久久久精品古装| 国产成人午夜福利电影在线观看| 欧美精品一区二区大全| 黑人高潮一二区| 美女内射精品一级片tv| 中文字幕亚洲精品专区| 国产伦理片在线播放av一区| 一区二区三区精品91| 久久影院123| 国产精品伦人一区二区| 亚洲欧美日韩东京热| 看非洲黑人一级黄片| 18禁在线无遮挡免费观看视频| 中文字幕免费在线视频6| 欧美日本中文国产一区发布| 深夜a级毛片| 黄色怎么调成土黄色| 色视频www国产| kizo精华| 男人舔奶头视频| 在线观看美女被高潮喷水网站| 欧美成人午夜免费资源| av天堂中文字幕网| 女人精品久久久久毛片| 99九九在线精品视频 | 国产精品国产三级专区第一集| 久久久久人妻精品一区果冻| 建设人人有责人人尽责人人享有的| 精品视频人人做人人爽| 91久久精品电影网| 啦啦啦在线观看免费高清www| 亚洲国产色片| 久久精品熟女亚洲av麻豆精品| 天天操日日干夜夜撸| 2018国产大陆天天弄谢| 国产淫片久久久久久久久| 久久狼人影院| 一级,二级,三级黄色视频| 搡老乐熟女国产| 建设人人有责人人尽责人人享有的| 欧美人与善性xxx| 少妇丰满av| 亚洲精品国产av成人精品| 涩涩av久久男人的天堂| 全区人妻精品视频| 欧美精品国产亚洲| 下体分泌物呈黄色| 精品视频人人做人人爽| a级毛片在线看网站| 成人影院久久| 天堂俺去俺来也www色官网| 国产伦在线观看视频一区| av播播在线观看一区| 色视频在线一区二区三区| 欧美 日韩 精品 国产| 久久久久久久久久久免费av| 搡女人真爽免费视频火全软件| 嫩草影院入口| 乱码一卡2卡4卡精品| 亚洲av日韩在线播放| 欧美另类一区| 久久99热6这里只有精品| 麻豆成人av视频| 日本色播在线视频| 黄色怎么调成土黄色| 男人狂女人下面高潮的视频| 久久女婷五月综合色啪小说| 久久影院123| 国产精品久久久久久精品古装| 亚洲av成人精品一二三区| 一级毛片电影观看| 在线免费观看不下载黄p国产| 18+在线观看网站| av播播在线观看一区| 五月天丁香电影| 夜夜爽夜夜爽视频| 天天躁夜夜躁狠狠久久av| 成人影院久久| 超碰97精品在线观看| 中文字幕人妻丝袜制服| 成人18禁高潮啪啪吃奶动态图 | 黄色日韩在线| 国产精品国产三级专区第一集| 欧美性感艳星| 91精品国产九色| 亚洲av二区三区四区| 桃花免费在线播放| 能在线免费看毛片的网站| 又爽又黄a免费视频| 色94色欧美一区二区| 日韩制服骚丝袜av| 热99国产精品久久久久久7| 赤兔流量卡办理| 91aial.com中文字幕在线观看| 久久综合国产亚洲精品| 免费看日本二区| 伊人久久精品亚洲午夜| 一区在线观看完整版| 校园人妻丝袜中文字幕| 国产成人精品一,二区| 自拍偷自拍亚洲精品老妇| 国产综合精华液| 国产深夜福利视频在线观看| 最近中文字幕2019免费版| 亚洲综合色惰| 人妻少妇偷人精品九色| 亚洲性久久影院| 国产精品秋霞免费鲁丝片| 2021少妇久久久久久久久久久| 两个人免费观看高清视频 | 亚洲欧美中文字幕日韩二区| 蜜桃在线观看..| 欧美bdsm另类| 日韩免费高清中文字幕av| av专区在线播放| 国产精品成人在线| 美女内射精品一级片tv| 黄色欧美视频在线观看| av在线观看视频网站免费| 男男h啪啪无遮挡| 国产精品麻豆人妻色哟哟久久| 男女边摸边吃奶| 最近中文字幕高清免费大全6| 国产在线一区二区三区精| 欧美精品一区二区大全| 观看美女的网站| av在线app专区| 男女边摸边吃奶| 少妇猛男粗大的猛烈进出视频| 99热这里只有是精品50| 亚洲精品亚洲一区二区| 91久久精品国产一区二区成人| 日韩欧美一区视频在线观看 | 日韩亚洲欧美综合| 欧美激情国产日韩精品一区| 成人国产av品久久久| 最近中文字幕高清免费大全6| 91成人精品电影| 亚洲真实伦在线观看| 午夜福利影视在线免费观看| 在线播放无遮挡| 国产精品国产三级国产av玫瑰| 国产精品国产三级国产专区5o| 大陆偷拍与自拍| 天堂8中文在线网| 涩涩av久久男人的天堂| 欧美老熟妇乱子伦牲交| 最黄视频免费看| 亚洲欧美中文字幕日韩二区| av有码第一页| 一级毛片电影观看| 又爽又黄a免费视频| 国产精品福利在线免费观看| 国产亚洲精品久久久com| 国产一区二区三区综合在线观看 | 中文字幕免费在线视频6| 一级,二级,三级黄色视频| 丰满少妇做爰视频| 日韩成人伦理影院| 日韩三级伦理在线观看| 国产乱来视频区| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 黄色欧美视频在线观看| 国产一区二区在线观看日韩| 欧美日韩av久久| 欧美区成人在线视频| 精品少妇久久久久久888优播| 亚洲丝袜综合中文字幕| 日韩制服骚丝袜av| 国产69精品久久久久777片| 人妻夜夜爽99麻豆av| 国产av国产精品国产| 少妇被粗大的猛进出69影院 | 日韩免费高清中文字幕av| 国产成人freesex在线| 黑人猛操日本美女一级片| 97在线人人人人妻| 观看av在线不卡| 国产视频内射| 我的老师免费观看完整版| 久久国产亚洲av麻豆专区| 高清视频免费观看一区二区| 2018国产大陆天天弄谢| 乱人伦中国视频| 国产午夜精品久久久久久一区二区三区| 久久国内精品自在自线图片| 亚洲国产日韩一区二区| 9色porny在线观看| 日韩成人伦理影院| 亚洲欧洲精品一区二区精品久久久 | 人人妻人人澡人人看| 亚洲国产色片| av天堂中文字幕网| 日本av手机在线免费观看| 亚洲人成网站在线观看播放| 黑人高潮一二区| 黑人巨大精品欧美一区二区蜜桃 | 欧美精品亚洲一区二区| 国产av码专区亚洲av| av在线播放精品| 日日摸夜夜添夜夜添av毛片| 2018国产大陆天天弄谢| 丁香六月天网| 视频区图区小说| 男人狂女人下面高潮的视频| 亚洲精品视频女| 亚洲欧美日韩另类电影网站| 午夜福利网站1000一区二区三区| 少妇猛男粗大的猛烈进出视频| 青春草视频在线免费观看| 2022亚洲国产成人精品| 亚洲四区av| 99视频精品全部免费 在线| 啦啦啦视频在线资源免费观看| 亚洲无线观看免费| 久久午夜综合久久蜜桃| h视频一区二区三区| 在线观看三级黄色| 黄色欧美视频在线观看| 2018国产大陆天天弄谢| 亚洲第一av免费看| 18禁裸乳无遮挡动漫免费视频| 伊人久久国产一区二区| 人人妻人人澡人人看| 如何舔出高潮| 哪个播放器可以免费观看大片| 国产中年淑女户外野战色| 亚洲av欧美aⅴ国产| 日韩亚洲欧美综合| 肉色欧美久久久久久久蜜桃| 嫩草影院入口| 看免费成人av毛片| 视频中文字幕在线观看| 在线播放无遮挡| 成人特级av手机在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲综合色惰| 亚洲图色成人| 最黄视频免费看| 国产精品久久久久久精品古装| 免费高清在线观看视频在线观看| 老司机亚洲免费影院| 国产av精品麻豆| 久久 成人 亚洲| 看免费成人av毛片| 校园人妻丝袜中文字幕| 日日摸夜夜添夜夜添av毛片| a级毛片免费高清观看在线播放| 精品一区二区三卡| 日韩中文字幕视频在线看片| videos熟女内射| 久久久久久伊人网av| av线在线观看网站| 人人妻人人澡人人看| 中文字幕免费在线视频6| 观看av在线不卡| 欧美日韩亚洲高清精品| 亚洲熟女精品中文字幕| 在线观看www视频免费| av在线观看视频网站免费| 日韩一区二区三区影片| 亚洲av国产av综合av卡| 国产美女午夜福利| 男女啪啪激烈高潮av片| 最近2019中文字幕mv第一页| 国产永久视频网站| 欧美xxxx性猛交bbbb| 天天躁夜夜躁狠狠久久av| 国产 精品1| 日韩av在线免费看完整版不卡| 久久精品夜色国产| 婷婷色麻豆天堂久久| 久久人人爽人人片av| 久久久午夜欧美精品| 青春草国产在线视频| 只有这里有精品99| 国产亚洲一区二区精品| 日本与韩国留学比较| 成人二区视频| 69精品国产乱码久久久| 久久久久久久国产电影| 亚洲怡红院男人天堂| av黄色大香蕉| 中国三级夫妇交换| 日本欧美视频一区| 国产亚洲午夜精品一区二区久久| 80岁老熟妇乱子伦牲交| 久久ye,这里只有精品| 日日摸夜夜添夜夜添av毛片| 亚洲情色 制服丝袜| 国产成人freesex在线| 日日爽夜夜爽网站| 黄色日韩在线|