• <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ì)比研究
    三種不銹鋼材料抗固體顆粒沖蝕性能研究
    一二三四在线观看免费中文在| 老熟妇乱子伦视频在线观看| 级片在线观看| 日本三级黄在线观看| 一级毛片高清免费大全| 亚洲一区二区三区色噜噜 | 亚洲av成人av| 国产免费现黄频在线看| 成人免费观看视频高清| 成人av一区二区三区在线看| 久久午夜亚洲精品久久| 久久九九热精品免费| 一级,二级,三级黄色视频| 亚洲九九香蕉| 极品教师在线免费播放| 大型av网站在线播放| 性色av乱码一区二区三区2| 嫩草影院精品99| 天堂√8在线中文| 婷婷精品国产亚洲av在线| 午夜免费鲁丝| 国产av精品麻豆| 老汉色av国产亚洲站长工具| 亚洲中文av在线| 久久久水蜜桃国产精品网| 99精品欧美一区二区三区四区| 亚洲精品一二三| 精品国产乱码久久久久久男人| 麻豆国产av国片精品| 国产高清videossex| 韩国av一区二区三区四区| 欧美乱妇无乱码| 久热这里只有精品99| 国产亚洲精品第一综合不卡| 国产av精品麻豆| 国产欧美日韩综合在线一区二区| 美女福利国产在线| 97碰自拍视频| 国产精品自产拍在线观看55亚洲| 国产麻豆69| 免费在线观看日本一区| 精品电影一区二区在线| 久久久久久久午夜电影 | 国产成人精品久久二区二区91| 久久精品影院6| 久久天躁狠狠躁夜夜2o2o| 91麻豆av在线| 老司机在亚洲福利影院| 黄色女人牲交| 一区二区三区精品91| 黑人猛操日本美女一级片| 精品国内亚洲2022精品成人| 免费在线观看日本一区| 午夜久久久在线观看| 老司机午夜十八禁免费视频| 亚洲精品国产色婷婷电影| 亚洲一区中文字幕在线| 亚洲一区二区三区色噜噜 | 岛国在线观看网站| 在线十欧美十亚洲十日本专区| 亚洲中文日韩欧美视频| 亚洲五月色婷婷综合| 国产在线观看jvid| 精品福利永久在线观看| 在线播放国产精品三级| 欧美日韩黄片免| 精品乱码久久久久久99久播| 成年人黄色毛片网站| 亚洲精品国产色婷婷电影| av免费在线观看网站| 男女下面进入的视频免费午夜 | 国产亚洲精品综合一区在线观看 | 国产精品野战在线观看 | 国产欧美日韩综合在线一区二区| 中文字幕高清在线视频| 波多野结衣av一区二区av| 免费高清在线观看日韩| 成人国语在线视频| 精品福利观看| 深夜精品福利| 18禁黄网站禁片午夜丰满| 黄色怎么调成土黄色| 动漫黄色视频在线观看| 亚洲熟妇中文字幕五十中出 | av天堂在线播放| 国产成人精品在线电影| 老司机福利观看| 一边摸一边抽搐一进一出视频| 国产av又大| 日韩欧美在线二视频| 美女国产高潮福利片在线看| 日韩欧美一区二区三区在线观看| 亚洲色图 男人天堂 中文字幕| 久久人妻av系列| 国产成人一区二区三区免费视频网站| 一级,二级,三级黄色视频| 黄色片一级片一级黄色片| 一级毛片精品| 国产男靠女视频免费网站| 露出奶头的视频| 在线观看午夜福利视频| 精品国产一区二区久久| 精品久久久久久电影网| 亚洲精品国产一区二区精华液| 国内久久婷婷六月综合欲色啪| e午夜精品久久久久久久| 免费一级毛片在线播放高清视频 | 99re在线观看精品视频| 日韩高清综合在线| av国产精品久久久久影院| 99久久久亚洲精品蜜臀av| 成年人免费黄色播放视频| 午夜精品久久久久久毛片777| 久久青草综合色| ponron亚洲| 国产一区二区在线av高清观看| 操出白浆在线播放| 超碰成人久久| 午夜免费观看网址| 美女 人体艺术 gogo| 一区二区三区国产精品乱码| 精品久久久久久,| 日韩欧美免费精品| 国产亚洲欧美精品永久| 国产精品一区二区三区四区久久 | 久久中文字幕人妻熟女| 日本黄色日本黄色录像| 高清毛片免费观看视频网站 | 啦啦啦免费观看视频1| 日韩av在线大香蕉| 黄色a级毛片大全视频| 一区二区三区激情视频| 久久狼人影院| 亚洲精品久久成人aⅴ小说| 亚洲男人天堂网一区| 男人的好看免费观看在线视频 | 免费高清在线观看日韩| 亚洲成国产人片在线观看| 日韩精品中文字幕看吧| 国产精品二区激情视频| 99久久精品国产亚洲精品| 超碰成人久久| 国产精品一区二区三区四区久久 | 久久香蕉国产精品| 色综合欧美亚洲国产小说| 欧美精品啪啪一区二区三区| 成年人黄色毛片网站| 国产成人影院久久av| 十分钟在线观看高清视频www| 国产亚洲精品一区二区www| 欧美在线黄色| 久久午夜亚洲精品久久| 欧美乱色亚洲激情| 欧美黄色片欧美黄色片| 女人精品久久久久毛片| 亚洲专区国产一区二区| 久久久国产精品麻豆| 国产精品亚洲av一区麻豆| 色老头精品视频在线观看| 亚洲视频免费观看视频| 老司机午夜十八禁免费视频| 777久久人妻少妇嫩草av网站| 亚洲,欧美精品.| 久久久国产成人免费| 久久久久久久久中文| 婷婷精品国产亚洲av在线| 男女之事视频高清在线观看| xxx96com| 一级a爱片免费观看的视频| 一级片免费观看大全| 久久国产精品人妻蜜桃| 亚洲欧美日韩无卡精品| 中文欧美无线码| 色综合婷婷激情| 国产精品久久久av美女十八| 午夜精品在线福利| 欧美国产精品va在线观看不卡| 婷婷六月久久综合丁香| bbb黄色大片| 一本大道久久a久久精品| 精品人妻在线不人妻| 超色免费av| 久久国产亚洲av麻豆专区| 一边摸一边抽搐一进一出视频| 亚洲专区国产一区二区| 国产精品一区二区精品视频观看| xxx96com| 国产精品野战在线观看 | 校园春色视频在线观看| 精品久久久精品久久久| 成人av一区二区三区在线看| 久久国产乱子伦精品免费另类| 国产免费av片在线观看野外av| 中文字幕人妻丝袜一区二区| 久久人妻福利社区极品人妻图片| 午夜老司机福利片| 99久久精品国产亚洲精品| 日本wwww免费看| 高清在线国产一区| 夫妻午夜视频| 久久影院123| 精品一区二区三区视频在线观看免费 | 中文字幕精品免费在线观看视频| 婷婷六月久久综合丁香| 色婷婷av一区二区三区视频| 久久久国产精品麻豆| 90打野战视频偷拍视频| 久久久久久久精品吃奶| 国产97色在线日韩免费| 久久国产精品人妻蜜桃| 在线观看舔阴道视频| 国产一区二区三区视频了| 国产精品自产拍在线观看55亚洲| 两个人看的免费小视频| 国产一区在线观看成人免费| 中文字幕人妻熟女乱码| 色老头精品视频在线观看| 天堂动漫精品| 亚洲精品国产一区二区精华液| 亚洲精品国产区一区二| 啪啪无遮挡十八禁网站| 亚洲av片天天在线观看| 岛国视频午夜一区免费看| 日本vs欧美在线观看视频| 午夜激情av网站| 日韩三级视频一区二区三区| 国产高清videossex| 97碰自拍视频| 成年人黄色毛片网站| 成人永久免费在线观看视频| 长腿黑丝高跟| 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久视频播放| 一区二区三区国产精品乱码| 在线免费观看的www视频| aaaaa片日本免费| 777久久人妻少妇嫩草av网站| 国产三级黄色录像| 91国产中文字幕| 午夜精品久久久久久毛片777| 高清在线国产一区| 国产精品电影一区二区三区| 国产精品99久久99久久久不卡| 欧美日本亚洲视频在线播放| 一二三四社区在线视频社区8| 成人手机av| 天堂影院成人在线观看| 变态另类成人亚洲欧美熟女 | 少妇被粗大的猛进出69影院| 国产av精品麻豆| 两人在一起打扑克的视频| 国产精品免费视频内射| 一夜夜www| 成人18禁在线播放| 悠悠久久av| 欧美激情 高清一区二区三区| 这个男人来自地球电影免费观看| av网站在线播放免费| 成年版毛片免费区| 亚洲一区中文字幕在线| 久久中文看片网| 日本撒尿小便嘘嘘汇集6| 久久香蕉精品热| 国产精品av久久久久免费| 欧美日韩黄片免| 欧美激情 高清一区二区三区| a级毛片黄视频| 欧美日韩精品网址| 欧美一级毛片孕妇| 很黄的视频免费| 咕卡用的链子| 午夜a级毛片| 99国产精品一区二区蜜桃av| 国产亚洲欧美98| 国产熟女xx| 99国产极品粉嫩在线观看| 中文字幕人妻丝袜一区二区| 一级作爱视频免费观看| 亚洲在线自拍视频| 亚洲精品在线美女| 在线免费观看的www视频| 日韩大码丰满熟妇| 高清av免费在线| 日韩国内少妇激情av| 国产精品香港三级国产av潘金莲| 1024视频免费在线观看| 精品久久久久久久久久免费视频 | 免费在线观看影片大全网站| 亚洲熟妇熟女久久| 午夜老司机福利片| 色综合婷婷激情| 亚洲 国产 在线| 在线免费观看的www视频| 桃红色精品国产亚洲av| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品成人av观看孕妇| 国产精品国产av在线观看| bbb黄色大片| 国产成人精品久久二区二区免费| 成人三级黄色视频| 亚洲av熟女| 国产精品一区二区在线不卡| 亚洲欧美日韩无卡精品| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 国产亚洲精品综合一区在线观看 | 99国产极品粉嫩在线观看| 免费在线观看完整版高清| 日本a在线网址| 国产激情久久老熟女| 亚洲五月色婷婷综合| 欧美黑人精品巨大| 高清毛片免费观看视频网站 | 欧美成人午夜精品| 国产不卡一卡二| 久久中文字幕人妻熟女| 亚洲一卡2卡3卡4卡5卡精品中文| 日日摸夜夜添夜夜添小说| 久久天堂一区二区三区四区| 老司机午夜十八禁免费视频| 手机成人av网站| 18禁国产床啪视频网站| 国产色视频综合| 国产成人av激情在线播放| 国产精品一区二区精品视频观看| 男女做爰动态图高潮gif福利片 | 国内久久婷婷六月综合欲色啪| 久久精品91蜜桃| 亚洲中文av在线| 女警被强在线播放| 在线国产一区二区在线| 99国产极品粉嫩在线观看| 国产伦人伦偷精品视频| 久久久国产成人精品二区 | 制服诱惑二区| 黄色成人免费大全| 欧美人与性动交α欧美软件| 变态另类成人亚洲欧美熟女 | 午夜免费成人在线视频| www.999成人在线观看| 亚洲精华国产精华精| 人妻久久中文字幕网| 大型黄色视频在线免费观看| 视频区图区小说| 12—13女人毛片做爰片一| 在线观看免费日韩欧美大片| 中文亚洲av片在线观看爽| 亚洲五月天丁香| 色婷婷久久久亚洲欧美| 亚洲欧美一区二区三区黑人| 日韩国内少妇激情av| 嫩草影院精品99| 看片在线看免费视频| 日本黄色日本黄色录像| 亚洲色图综合在线观看| 国产午夜精品久久久久久| 久久久久久久久免费视频了| 99国产精品免费福利视频| 正在播放国产对白刺激| 国产av一区二区精品久久| 精品久久蜜臀av无| 国产精品久久久av美女十八| 精品久久蜜臀av无| 国产又色又爽无遮挡免费看| 精品久久蜜臀av无| 亚洲成人免费av在线播放| 午夜老司机福利片| 国产av在哪里看| 午夜老司机福利片| 三上悠亚av全集在线观看| tocl精华| 黄色 视频免费看| 午夜免费激情av| 欧美日韩av久久| 欧美日韩亚洲高清精品| 另类亚洲欧美激情| 成人黄色视频免费在线看| 好看av亚洲va欧美ⅴa在| 97碰自拍视频| 精品福利永久在线观看| 久久久国产成人精品二区 | 国产亚洲精品久久久久5区| 精品卡一卡二卡四卡免费| 国产午夜精品久久久久久| 国产精品永久免费网站| 欧美+亚洲+日韩+国产| 欧美激情高清一区二区三区| 国产精品久久电影中文字幕| 精品国产超薄肉色丝袜足j| 久久久久久大精品| 正在播放国产对白刺激| 日本wwww免费看| 久久久久久久久中文| 丁香欧美五月| 丝袜人妻中文字幕| 91老司机精品| 亚洲 欧美 日韩 在线 免费| 叶爱在线成人免费视频播放| 久久热在线av| 在线观看一区二区三区| 午夜成年电影在线免费观看| 精品人妻在线不人妻| 美女 人体艺术 gogo| 欧美成人午夜精品| 波多野结衣一区麻豆| 国产成人欧美| 精品无人区乱码1区二区| 嫩草影视91久久| 热99re8久久精品国产| 日韩精品青青久久久久久| 无人区码免费观看不卡| 国产男靠女视频免费网站| 亚洲一区高清亚洲精品| 国产又爽黄色视频| 欧美中文综合在线视频| 91麻豆精品激情在线观看国产 | 欧美黑人欧美精品刺激| 国产成人免费无遮挡视频| 女生性感内裤真人,穿戴方法视频| 中文字幕av电影在线播放| 国产精品九九99| 美女福利国产在线| 少妇裸体淫交视频免费看高清 | 啦啦啦 在线观看视频| 黄色毛片三级朝国网站| 欧美日韩亚洲综合一区二区三区_| 香蕉久久夜色| 欧美在线一区亚洲| 亚洲精品成人av观看孕妇| 久久精品亚洲av国产电影网| 欧美一区二区精品小视频在线| 无人区码免费观看不卡| 欧美av亚洲av综合av国产av| 性欧美人与动物交配| 亚洲人成伊人成综合网2020| 老司机在亚洲福利影院| 亚洲九九香蕉| 人妻丰满熟妇av一区二区三区| 女人高潮潮喷娇喘18禁视频| 一级片免费观看大全| 国产精品亚洲一级av第二区| www.999成人在线观看| 美女扒开内裤让男人捅视频| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| 精品卡一卡二卡四卡免费| 国产成人福利小说| 亚洲精品一区av在线观看| 久久久久久久精品吃奶| 在线看三级毛片| 小说图片视频综合网站| 欧美性感艳星| 欧美另类亚洲清纯唯美| 国产白丝娇喘喷水9色精品| 亚洲熟妇熟女久久| 性欧美人与动物交配| 18+在线观看网站| 久久精品国产清高在天天线| 直男gayav资源| 美女高潮喷水抽搐中文字幕| 国产激情偷乱视频一区二区| 一个人免费在线观看电影| 又黄又爽又免费观看的视频| 久久午夜福利片| 动漫黄色视频在线观看| 亚洲国产色片| 色在线成人网| 免费在线观看亚洲国产| 亚洲精品亚洲一区二区| 久9热在线精品视频| 小说图片视频综合网站| 免费观看人在逋| 亚洲中文日韩欧美视频| 久久精品国产亚洲av香蕉五月| 亚洲一区高清亚洲精品| 男女做爰动态图高潮gif福利片| av在线老鸭窝| 亚洲最大成人av| 久久精品国产亚洲av涩爱 | 真人做人爱边吃奶动态| 91在线精品国自产拍蜜月| 白带黄色成豆腐渣| 久久性视频一级片| 91久久精品电影网| 成年女人毛片免费观看观看9| 国内揄拍国产精品人妻在线| 人妻丰满熟妇av一区二区三区| 长腿黑丝高跟| 国产成人a区在线观看| 麻豆久久精品国产亚洲av| 一区二区三区四区激情视频 | 精品福利观看| 噜噜噜噜噜久久久久久91| 淫妇啪啪啪对白视频| 在线观看一区二区三区| 精品无人区乱码1区二区| ponron亚洲| 两人在一起打扑克的视频| 少妇丰满av| 精品一区二区三区av网在线观看| 亚洲狠狠婷婷综合久久图片| 国产白丝娇喘喷水9色精品| 久久久色成人| 日韩成人在线观看一区二区三区| 日韩欧美国产在线观看| 淫秽高清视频在线观看| 国语自产精品视频在线第100页| 亚洲av免费高清在线观看| 国产三级中文精品| 亚洲av电影不卡..在线观看| 国产精品一区二区三区四区久久| 亚洲一区二区三区不卡视频| 性色avwww在线观看| 91午夜精品亚洲一区二区三区 | 搞女人的毛片| 18禁裸乳无遮挡免费网站照片| ponron亚洲| 精品午夜福利视频在线观看一区| 久久久久久久久大av| 国产精品野战在线观看| 日日干狠狠操夜夜爽| 动漫黄色视频在线观看| 欧美丝袜亚洲另类 | 欧美日韩乱码在线| 女同久久另类99精品国产91| x7x7x7水蜜桃| 丰满的人妻完整版| 午夜日韩欧美国产| 亚洲精品影视一区二区三区av| 91久久精品电影网| 国产日本99.免费观看| 欧美日韩福利视频一区二区| 国内精品美女久久久久久| 一本久久中文字幕| 99视频精品全部免费 在线| 嫩草影院新地址| 成人永久免费在线观看视频| 黄色视频,在线免费观看| 日日摸夜夜添夜夜添小说| 国产av不卡久久| 中文字幕久久专区| 精品午夜福利在线看| 久久国产精品人妻蜜桃| 免费大片18禁| 又黄又爽又免费观看的视频| 国产精品久久视频播放| 99视频精品全部免费 在线| 欧美成人一区二区免费高清观看| 老熟妇乱子伦视频在线观看| 老鸭窝网址在线观看| 日日摸夜夜添夜夜添av毛片 | 免费大片18禁| 国产午夜福利久久久久久| 欧美中文日本在线观看视频| 88av欧美| 99国产精品一区二区三区| 国产亚洲av嫩草精品影院| 自拍偷自拍亚洲精品老妇| 天堂影院成人在线观看| 俄罗斯特黄特色一大片| 久久精品91蜜桃| 麻豆av噜噜一区二区三区| 97超视频在线观看视频| 搞女人的毛片| 麻豆成人午夜福利视频| 男人舔女人下体高潮全视频| 看片在线看免费视频| 最近最新免费中文字幕在线| 97碰自拍视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 可以在线观看的亚洲视频| 国产男靠女视频免费网站| 欧美成人免费av一区二区三区| 在线观看一区二区三区| 直男gayav资源| 亚洲自偷自拍三级| 精品久久久久久久久久免费视频| 我的女老师完整版在线观看| 亚洲在线自拍视频| 性色av乱码一区二区三区2| 在线天堂最新版资源| 欧美中文日本在线观看视频| 长腿黑丝高跟| 国产黄片美女视频| 99热6这里只有精品| 国产91精品成人一区二区三区| netflix在线观看网站| 亚洲,欧美,日韩| 黄色视频,在线免费观看| 老司机深夜福利视频在线观看| 亚洲av第一区精品v没综合| 在线十欧美十亚洲十日本专区| 女人十人毛片免费观看3o分钟| 少妇高潮的动态图| 欧美3d第一页| 18禁黄网站禁片免费观看直播| 亚洲av不卡在线观看| 婷婷六月久久综合丁香| 九色国产91popny在线| 亚洲综合色惰| 日本一本二区三区精品| 99在线人妻在线中文字幕| 色5月婷婷丁香| 久久亚洲真实| 国内精品美女久久久久久| 狂野欧美白嫩少妇大欣赏| 成人三级黄色视频| 国产久久久一区二区三区| 丰满人妻一区二区三区视频av| 午夜精品一区二区三区免费看| 综合色av麻豆| 国产精品久久久久久久电影| 99久国产av精品| 亚洲无线在线观看|