• <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ì)比研究
    三種不銹鋼材料抗固體顆粒沖蝕性能研究
    捣出白浆h1v1| 男女啪啪激烈高潮av片| 欧美bdsm另类| 亚洲,一卡二卡三卡| 日韩 亚洲 欧美在线| 黑人巨大精品欧美一区二区蜜桃| 国产精品不卡视频一区二区| av在线观看视频网站免费| 一级片'在线观看视频| 亚洲国产最新在线播放| 久久久久久人人人人人| 亚洲在久久综合| 男女边摸边吃奶| 成年女人毛片免费观看观看9 | 精品少妇一区二区三区视频日本电影 | 午夜免费观看性视频| 亚洲美女搞黄在线观看| 九色亚洲精品在线播放| 亚洲欧美中文字幕日韩二区| 99久久精品国产国产毛片| 久久热在线av| 国产精品欧美亚洲77777| 搡女人真爽免费视频火全软件| 各种免费的搞黄视频| 一级片'在线观看视频| 天堂8中文在线网| 一二三四在线观看免费中文在| tube8黄色片| 国产探花极品一区二区| 国产熟女午夜一区二区三区| 久久午夜综合久久蜜桃| 亚洲av综合色区一区| kizo精华| 亚洲一区二区三区欧美精品| 国产亚洲午夜精品一区二区久久| 亚洲精品成人av观看孕妇| 久久精品熟女亚洲av麻豆精品| 熟女电影av网| 又粗又硬又长又爽又黄的视频| 人妻人人澡人人爽人人| 美女福利国产在线| 乱人伦中国视频| 女性生殖器流出的白浆| 国产欧美日韩综合在线一区二区| 91精品伊人久久大香线蕉| 久久久久久久大尺度免费视频| 精品人妻熟女毛片av久久网站| 欧美亚洲 丝袜 人妻 在线| 男女边吃奶边做爰视频| 亚洲av免费高清在线观看| av在线老鸭窝| 成人毛片60女人毛片免费| 精品久久久久久电影网| 99久久综合免费| 久久久久久久国产电影| 人体艺术视频欧美日本| 亚洲内射少妇av| 中文乱码字字幕精品一区二区三区| 老女人水多毛片| 久久国产精品大桥未久av| 国产极品粉嫩免费观看在线| 国产免费一区二区三区四区乱码| 狠狠精品人妻久久久久久综合| 国产在线视频一区二区| 成人国产麻豆网| 国精品久久久久久国模美| 亚洲国产欧美日韩在线播放| 在线天堂中文资源库| 免费观看a级毛片全部| 亚洲av电影在线观看一区二区三区| 一二三四在线观看免费中文在| 精品国产一区二区三区四区第35| 日韩在线高清观看一区二区三区| 女性生殖器流出的白浆| 亚洲成国产人片在线观看| 一本久久精品| 成年女人在线观看亚洲视频| 亚洲av日韩在线播放| 可以免费在线观看a视频的电影网站 | 午夜激情av网站| 日韩电影二区| 久久精品亚洲av国产电影网| 国产1区2区3区精品| 女人被躁到高潮嗷嗷叫费观| 老司机影院毛片| 亚洲国产av影院在线观看| 国产成人a∨麻豆精品| 精品99又大又爽又粗少妇毛片| 亚洲中文av在线| 成人国产麻豆网| 精品卡一卡二卡四卡免费| 久久精品国产鲁丝片午夜精品| 波野结衣二区三区在线| 爱豆传媒免费全集在线观看| 我要看黄色一级片免费的| 国产xxxxx性猛交| 国产成人精品无人区| www.精华液| 国产精品99久久99久久久不卡 | 交换朋友夫妻互换小说| 亚洲av电影在线观看一区二区三区| 久久精品久久久久久久性| 国产精品二区激情视频| 日韩制服骚丝袜av| 国产片特级美女逼逼视频| 日韩精品免费视频一区二区三区| 美女大奶头黄色视频| 97在线视频观看| 少妇人妻精品综合一区二区| 久久鲁丝午夜福利片| 26uuu在线亚洲综合色| 亚洲成人手机| 成年动漫av网址| 精品人妻偷拍中文字幕| 久久精品久久精品一区二区三区| 纯流量卡能插随身wifi吗| 免费高清在线观看日韩| 2018国产大陆天天弄谢| 亚洲av综合色区一区| 国产在线免费精品| 成人国语在线视频| 日韩一区二区三区影片| 国产成人精品在线电影| 欧美另类一区| 免费黄网站久久成人精品| 99精国产麻豆久久婷婷| 亚洲精品国产色婷婷电影| 麻豆乱淫一区二区| 亚洲精品美女久久久久99蜜臀 | 亚洲精品久久久久久婷婷小说| 观看美女的网站| 如何舔出高潮| 久久99一区二区三区| 国产片内射在线| 欧美亚洲日本最大视频资源| 大香蕉久久网| 一区二区三区精品91| 日韩中文字幕视频在线看片| 国产一区亚洲一区在线观看| 欧美精品av麻豆av| av卡一久久| 麻豆乱淫一区二区| 午夜影院在线不卡| 尾随美女入室| 人成视频在线观看免费观看| 日日撸夜夜添| 亚洲av电影在线观看一区二区三区| 一本大道久久a久久精品| 热re99久久精品国产66热6| 男女无遮挡免费网站观看| 国产成人精品在线电影| 街头女战士在线观看网站| 一级片'在线观看视频| 免费黄色在线免费观看| 国产免费现黄频在线看| 欧美另类一区| 久久精品国产亚洲av天美| 免费观看在线日韩| 青青草视频在线视频观看| 久久久久久人妻| 狠狠婷婷综合久久久久久88av| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| av有码第一页| 欧美国产精品一级二级三级| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| 久久综合国产亚洲精品| videos熟女内射| a级毛片在线看网站| 久久久久久久精品精品| 日韩电影二区| 搡女人真爽免费视频火全软件| av国产精品久久久久影院| 九九爱精品视频在线观看| 在线看a的网站| 亚洲欧洲精品一区二区精品久久久 | av女优亚洲男人天堂| 宅男免费午夜| 18禁裸乳无遮挡动漫免费视频| 国产黄色视频一区二区在线观看| 久久精品人人爽人人爽视色| 2018国产大陆天天弄谢| 国产一级毛片在线| 制服丝袜香蕉在线| 免费不卡的大黄色大毛片视频在线观看| 久久热在线av| 国产亚洲精品第一综合不卡| 亚洲精品日韩在线中文字幕| 国语对白做爰xxxⅹ性视频网站| 日韩精品免费视频一区二区三区| 搡女人真爽免费视频火全软件| av女优亚洲男人天堂| 亚洲精品美女久久久久99蜜臀 | 婷婷色综合大香蕉| 亚洲国产色片| 久久女婷五月综合色啪小说| 两个人看的免费小视频| 男人添女人高潮全过程视频| 视频区图区小说| 不卡av一区二区三区| 亚洲精华国产精华液的使用体验| 国产麻豆69| 一二三四在线观看免费中文在| 国产成人精品福利久久| 久久久久人妻精品一区果冻| av在线老鸭窝| 国产国语露脸激情在线看| 久久精品国产综合久久久| 老女人水多毛片| 最新的欧美精品一区二区| 嫩草影院入口| 精品少妇黑人巨大在线播放| 免费av中文字幕在线| 伊人久久大香线蕉亚洲五| 中文天堂在线官网| av电影中文网址| 免费看不卡的av| 最新中文字幕久久久久| 黄网站色视频无遮挡免费观看| 久久国产精品大桥未久av| 久久ye,这里只有精品| 十八禁高潮呻吟视频| 婷婷色麻豆天堂久久| 王馨瑶露胸无遮挡在线观看| 国产一区二区三区综合在线观看| 欧美中文综合在线视频| 国产亚洲av片在线观看秒播厂| 视频在线观看一区二区三区| 日韩欧美精品免费久久| 日本黄色日本黄色录像| 国产男人的电影天堂91| 极品少妇高潮喷水抽搐| 久久精品久久久久久噜噜老黄| 97人妻天天添夜夜摸| 久久女婷五月综合色啪小说| 久久青草综合色| 久久国产精品男人的天堂亚洲| 国产精品一国产av| 男人爽女人下面视频在线观看| 美女高潮到喷水免费观看| 久久久国产一区二区| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 精品卡一卡二卡四卡免费| 亚洲精品一二三| 日韩大片免费观看网站| 九九爱精品视频在线观看| 婷婷色麻豆天堂久久| 亚洲精品国产av成人精品| 久久久久精品人妻al黑| 丰满迷人的少妇在线观看| 在线观看人妻少妇| 久久久精品国产亚洲av高清涩受| 美女视频免费永久观看网站| 亚洲av欧美aⅴ国产| 国产成人免费无遮挡视频| 久久精品熟女亚洲av麻豆精品| 成人国产麻豆网| 国产乱来视频区| 狠狠婷婷综合久久久久久88av| 人人妻人人添人人爽欧美一区卜| kizo精华| 丰满乱子伦码专区| 精品国产一区二区三区四区第35| 成人国产av品久久久| 亚洲欧美日韩另类电影网站| 午夜91福利影院| 成人漫画全彩无遮挡| av不卡在线播放| 天天躁夜夜躁狠狠躁躁| 人人妻人人澡人人爽人人夜夜| 亚洲伊人久久精品综合| 一二三四中文在线观看免费高清| 国产成人欧美| 九草在线视频观看| 欧美日韩精品成人综合77777| 国产成人av激情在线播放| 欧美精品亚洲一区二区| 日韩电影二区| 电影成人av| 国产精品女同一区二区软件| 亚洲国产欧美日韩在线播放| 黄色配什么色好看| 国产国语露脸激情在线看| 美国免费a级毛片| 亚洲国产毛片av蜜桃av| 欧美xxⅹ黑人| 少妇人妻久久综合中文| 黄频高清免费视频| 美女国产视频在线观看| 又粗又硬又长又爽又黄的视频| 国精品久久久久久国模美| 亚洲精品成人av观看孕妇| 18禁国产床啪视频网站| 亚洲成人av在线免费| 久久精品久久久久久噜噜老黄| 青春草亚洲视频在线观看| 午夜日韩欧美国产| 一个人免费看片子| 叶爱在线成人免费视频播放| 午夜激情久久久久久久| 午夜激情av网站| 国产伦理片在线播放av一区| av福利片在线| 咕卡用的链子| 午夜福利乱码中文字幕| 午夜日本视频在线| 亚洲欧美一区二区三区国产| av网站免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 人人妻人人澡人人爽人人夜夜| 黄色怎么调成土黄色| 国产精品国产av在线观看| 国产男人的电影天堂91| 亚洲三区欧美一区| 两个人免费观看高清视频| 国产精品蜜桃在线观看| 国产av码专区亚洲av| 久久这里只有精品19| 性色avwww在线观看| 国产日韩欧美亚洲二区| 不卡视频在线观看欧美| 亚洲精品视频女| 欧美国产精品一级二级三级| 最近2019中文字幕mv第一页| 国产精品不卡视频一区二区| 午夜日本视频在线| 十分钟在线观看高清视频www| 国产av国产精品国产| 国产精品蜜桃在线观看| 99re6热这里在线精品视频| 免费观看a级毛片全部| videossex国产| 日韩电影二区| 亚洲av国产av综合av卡| 最近中文字幕高清免费大全6| 亚洲美女搞黄在线观看| 亚洲av日韩在线播放| 啦啦啦在线观看免费高清www| 欧美国产精品va在线观看不卡| 飞空精品影院首页| 国产野战对白在线观看| 大陆偷拍与自拍| 亚洲综合色惰| 一区二区三区四区激情视频| 黄色视频在线播放观看不卡| 国产精品香港三级国产av潘金莲 | 精品少妇久久久久久888优播| 丝瓜视频免费看黄片| 纵有疾风起免费观看全集完整版| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三区在线| 九草在线视频观看| 国产极品粉嫩免费观看在线| 亚洲欧洲精品一区二区精品久久久 | 免费少妇av软件| 男女免费视频国产| 欧美精品国产亚洲| 亚洲欧美精品综合一区二区三区 | 18在线观看网站| 校园人妻丝袜中文字幕| 999久久久国产精品视频| 成年美女黄网站色视频大全免费| 99久久中文字幕三级久久日本| av国产久精品久网站免费入址| 精品99又大又爽又粗少妇毛片| 欧美xxⅹ黑人| 极品人妻少妇av视频| 国产日韩一区二区三区精品不卡| 国产无遮挡羞羞视频在线观看| 亚洲男人天堂网一区| 国产精品免费大片| 国产午夜精品一二区理论片| 成人18禁高潮啪啪吃奶动态图| 18+在线观看网站| 亚洲av成人精品一二三区| 下体分泌物呈黄色| 在线精品无人区一区二区三| 只有这里有精品99| av国产精品久久久久影院| 在线观看www视频免费| 国产成人精品福利久久| 亚洲国产av新网站| 曰老女人黄片| 成年动漫av网址| 在线观看免费日韩欧美大片| 国产亚洲午夜精品一区二区久久| 性色av一级| 国产白丝娇喘喷水9色精品| 日韩一卡2卡3卡4卡2021年| 久久青草综合色| 久久久精品区二区三区| 看免费成人av毛片| 日韩大片免费观看网站| a级毛片黄视频| 中文字幕人妻丝袜一区二区 | 国产极品天堂在线| 久久久久久久精品精品| 亚洲国产成人一精品久久久| 日本免费在线观看一区| 久久久精品94久久精品| 女性被躁到高潮视频| h视频一区二区三区| 精品人妻偷拍中文字幕| 人人澡人人妻人| 韩国精品一区二区三区| 人成视频在线观看免费观看| 精品一区在线观看国产| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 观看美女的网站| 国产xxxxx性猛交| 视频区图区小说| 韩国av在线不卡| 人人妻人人澡人人爽人人夜夜| 国产片内射在线| 免费黄频网站在线观看国产| 亚洲精品久久久久久婷婷小说| 午夜老司机福利剧场| 午夜影院在线不卡| 天天躁夜夜躁狠狠久久av| 国产男女内射视频| 婷婷色av中文字幕| 男女下面插进去视频免费观看| 国产精品欧美亚洲77777| 国产日韩欧美视频二区| 色94色欧美一区二区| 啦啦啦在线观看免费高清www| 欧美亚洲日本最大视频资源| 国产av码专区亚洲av| 久久99热这里只频精品6学生| 美女中出高潮动态图| 国产麻豆69| 婷婷色av中文字幕| 交换朋友夫妻互换小说| 午夜91福利影院| 欧美激情 高清一区二区三区| 日本wwww免费看| 亚洲一区中文字幕在线| 精品少妇黑人巨大在线播放| 日韩大片免费观看网站| 一边摸一边做爽爽视频免费| 日韩,欧美,国产一区二区三区| 国产一区亚洲一区在线观看| 天堂俺去俺来也www色官网| 最近最新中文字幕大全免费视频 | 欧美精品高潮呻吟av久久| 欧美精品高潮呻吟av久久| 日本午夜av视频| 欧美人与性动交α欧美软件| 久久精品国产亚洲av涩爱| 亚洲国产欧美在线一区| 欧美成人午夜免费资源| 亚洲一区二区三区欧美精品| 国产精品久久久久久av不卡| 久久精品aⅴ一区二区三区四区 | 一本—道久久a久久精品蜜桃钙片| 我要看黄色一级片免费的| 91精品三级在线观看| 激情五月婷婷亚洲| 欧美成人午夜免费资源| 在线观看免费日韩欧美大片| 一区二区三区乱码不卡18| 不卡视频在线观看欧美| 欧美中文综合在线视频| 18禁动态无遮挡网站| 在线天堂最新版资源| 青青草视频在线视频观看| www.自偷自拍.com| 母亲3免费完整高清在线观看 | 少妇人妻精品综合一区二区| 大陆偷拍与自拍| xxxhd国产人妻xxx| 精品一品国产午夜福利视频| av免费在线看不卡| 在线观看www视频免费| 亚洲欧美精品自产自拍| 免费观看性生交大片5| 巨乳人妻的诱惑在线观看| 又大又黄又爽视频免费| 99热全是精品| 边亲边吃奶的免费视频| 校园人妻丝袜中文字幕| 国产精品香港三级国产av潘金莲 | 久热久热在线精品观看| 欧美xxⅹ黑人| 免费观看性生交大片5| 久久这里有精品视频免费| 国产成人免费观看mmmm| 女人被躁到高潮嗷嗷叫费观| 亚洲美女黄色视频免费看| 精品少妇一区二区三区视频日本电影 | 黄片播放在线免费| 日本欧美视频一区| 欧美日韩一级在线毛片| 一级黄片播放器| 制服丝袜香蕉在线| 卡戴珊不雅视频在线播放| 欧美+日韩+精品| 啦啦啦在线观看免费高清www| 黑人巨大精品欧美一区二区蜜桃| 观看av在线不卡| 丰满饥渴人妻一区二区三| 1024视频免费在线观看| 狂野欧美激情性bbbbbb| 成人毛片a级毛片在线播放| 亚洲国产最新在线播放| 久热这里只有精品99| 中文欧美无线码| 亚洲四区av| 亚洲经典国产精华液单| 黑人欧美特级aaaaaa片| 一二三四在线观看免费中文在| 午夜福利一区二区在线看| 久久精品国产自在天天线| 久久国产亚洲av麻豆专区| 午夜福利乱码中文字幕| 国产精品亚洲av一区麻豆 | 精品少妇一区二区三区视频日本电影 | 国语对白做爰xxxⅹ性视频网站| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产一区二区 视频在线| 日韩电影二区| 成年女人毛片免费观看观看9 | 精品少妇内射三级| 日本黄色日本黄色录像| 人人妻人人添人人爽欧美一区卜| 伦精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 欧美变态另类bdsm刘玥| 中国三级夫妇交换| 18禁动态无遮挡网站| 麻豆av在线久日| 国产成人aa在线观看| 国产白丝娇喘喷水9色精品| 日本av手机在线免费观看| 日韩精品免费视频一区二区三区| 亚洲精品乱久久久久久| 久久精品国产a三级三级三级| 欧美日韩精品成人综合77777| 婷婷成人精品国产| 青草久久国产| 性色av一级| 久久99精品国语久久久| 久久精品夜色国产| 精品福利永久在线观看| 王馨瑶露胸无遮挡在线观看| 国产在线免费精品| 制服人妻中文乱码| 精品一品国产午夜福利视频| 可以免费在线观看a视频的电影网站 | av网站在线播放免费| 精品亚洲成国产av| 搡女人真爽免费视频火全软件| 韩国精品一区二区三区| 中文字幕色久视频| 精品人妻熟女毛片av久久网站| 免费日韩欧美在线观看| 美女主播在线视频| 亚洲精品久久成人aⅴ小说| 中文字幕av电影在线播放| 色网站视频免费| 亚洲四区av| 一级,二级,三级黄色视频| 青春草国产在线视频| 日韩av在线免费看完整版不卡| 麻豆精品久久久久久蜜桃| 亚洲av中文av极速乱| 亚洲欧美日韩另类电影网站| 视频在线观看一区二区三区| 97人妻天天添夜夜摸| 人成视频在线观看免费观看| 国产 精品1| 欧美精品人与动牲交sv欧美| 热99国产精品久久久久久7| 亚洲一码二码三码区别大吗| 在现免费观看毛片| 精品福利永久在线观看| 免费黄色在线免费观看| freevideosex欧美| 欧美亚洲日本最大视频资源| 99热国产这里只有精品6| 日韩欧美一区视频在线观看| 亚洲国产欧美网| 免费女性裸体啪啪无遮挡网站| 午夜激情av网站| 国产一区二区三区av在线| 丝瓜视频免费看黄片| 欧美人与性动交α欧美精品济南到 | 波多野结衣一区麻豆| 热99国产精品久久久久久7| 国产欧美亚洲国产| 国产精品 欧美亚洲| 韩国精品一区二区三区| 国产黄频视频在线观看| 人妻少妇偷人精品九色| 老鸭窝网址在线观看| 丝袜美足系列| 十八禁网站网址无遮挡| 国产精品 欧美亚洲| 国产xxxxx性猛交| 国产人伦9x9x在线观看 | 男女边吃奶边做爰视频| 国产精品二区激情视频| 王馨瑶露胸无遮挡在线观看| 亚洲av日韩在线播放| 欧美成人午夜精品| 免费在线观看黄色视频的| 天天躁狠狠躁夜夜躁狠狠躁| 各种免费的搞黄视频| av在线app专区| 久久这里只有精品19| 久久久久久久大尺度免费视频| 中文天堂在线官网|