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

    加氫空冷器注劑T型管內(nèi)氣-液流動(dòng)特性的數(shù)值模擬

    2017-10-16 11:29:24金浩哲劉文文偶國(guó)富陳小平李鵬軒
    關(guān)鍵詞:型管分率剪切應(yīng)力

    金浩哲, 劉文文, 偶國(guó)富, 陳小平, 李鵬軒

    (1.浙江理工大學(xué) 流動(dòng)腐蝕研究所, 浙江 杭州 310018; 2.杭州市特種設(shè)備檢測(cè)研究院, 浙江 杭州 310003)

    加氫空冷器注劑T型管內(nèi)氣-液流動(dòng)特性的數(shù)值模擬

    金浩哲1, 劉文文1, 偶國(guó)富1, 陳小平1, 李鵬軒2

    (1.浙江理工大學(xué) 流動(dòng)腐蝕研究所, 浙江 杭州 310018; 2.杭州市特種設(shè)備檢測(cè)研究院, 浙江 杭州 310003)

    以加氫反應(yīng)流出物空冷器(簡(jiǎn)稱加氫空冷器)系統(tǒng)前的工藝注劑T型管為研究對(duì)象,分析其工藝過程,基于工藝注劑T型管的幾何特征建立結(jié)構(gòu)離散模型,選擇Mixture模型和Realizablek-ε模型對(duì)其在操作工況下的流動(dòng)特性進(jìn)行計(jì)算流體動(dòng)力學(xué)(CFD)模擬,考察不同注劑T型管結(jié)構(gòu)的管內(nèi)水相分率、湍動(dòng)能、剪切應(yīng)力等流體動(dòng)力學(xué)參數(shù)的分布規(guī)律,揭示氣-液兩相流動(dòng)作用下注劑T型管段的沖蝕減薄現(xiàn)象的形成機(jī)理。研究結(jié)果表明:水相分率及剪切應(yīng)力最大的位置,即橫坐標(biāo)x位于-0.15~-0.26 m之間,管道易發(fā)生穿孔泄漏。超聲波測(cè)厚儀測(cè)試獲得的管道剩余壁厚與數(shù)值預(yù)測(cè)的結(jié)果基本吻合。

    加氫空冷器; 注劑T型管; 多相流; 數(shù)值模擬; 失效分析

    Abstract: In order to reveal the erosion thinning formation mechanism with the action of gas-liquid two-phase in the injection T-tube section before the hydrogenation reaction effluent air cooler (REAC) system, the process injection T-tube was taken as the study object, the relevant technical process was analyzed and the discrete model of structure based on its geometric feature was built. Subsequently the Mixture model and Realizablek-εmodel were chosen in this computational fluid dynamics (CFD) simulation for researching flow characteristics under operation condition, including water phase fraction, turbulent kinetic energy, wall shear stress and other parameters’ distribution.The results indicated that: at the position of maximum value in water phase fraction and shear stress, i.e. betweenx=-0.15 m andx=-0.26 m, the pipe is prone to blockage and leak, which reasonably agrees with the distribution of pipe remaining wall thickness obtained by using ultrasonic thickness gauges.

    Keywords:hydrogenation air cooler; injection T-tube; multiphase flow; numerical simulation; failure analysis

    加氫反應(yīng)流出物空冷器(Reactor Effluent Air Coolers)(簡(jiǎn)稱加氫空冷器)是石油化工、煤化工等流程型工業(yè)的重要裝備系統(tǒng),長(zhǎng)期承受高溫、高壓、臨氫工況,運(yùn)行過程中的風(fēng)險(xiǎn)極高[1-2]。特別是在裝置大型化、原料劣質(zhì)化和工況苛刻化的發(fā)展過程中,加氫空冷器系統(tǒng)普遍出現(xiàn)了因流動(dòng)腐蝕引發(fā)的管束穿孔、泄漏等事故[3-5],不僅嚴(yán)重影響企業(yè)的生產(chǎn)計(jì)劃,而且嚴(yán)重威脅著節(jié)能、環(huán)保與安全。

    原料油富含N和S的化合物,并含有少量的氯化物,經(jīng)加氫反應(yīng)生成NH4Cl和NH4HS,隨著溫度的降低,銨鹽直接由氣相冷凝成固相并迅速堵塞管束[6-7]。金浩哲等[8-9]基于熱力學(xué)推導(dǎo)的方法,建立了銨鹽的結(jié)晶溫度平衡曲線。為避免銨鹽結(jié)晶堵塞空冷管束,美國(guó)石油學(xué)會(huì)(API)給出了推薦性的加氫空冷器系統(tǒng)設(shè)計(jì)、制造、運(yùn)行和檢驗(yàn)標(biāo)準(zhǔn)API 932-A[10],API 932-B[11],并建議在加氫空冷器前設(shè)置工藝注劑(注水)點(diǎn)對(duì)銨鹽進(jìn)行溶解洗滌。由于工藝注劑T型管道內(nèi)的介質(zhì)組成復(fù)雜,氣-液流動(dòng)過程中的流體動(dòng)力學(xué)參數(shù)分布不均,極易引發(fā)注劑T型管道的穿孔泄漏[12]。因此,建立工藝注劑T型管道的氣-液流動(dòng)特性評(píng)價(jià)方法,并進(jìn)行針對(duì)性的設(shè)計(jì)和操作優(yōu)化是學(xué)術(shù)界關(guān)注的焦點(diǎn)之一。當(dāng)前,與工藝注劑T型管道氣-液流動(dòng)特性相關(guān)的研究成果主要包括:Wang等[13]針對(duì)天然氣領(lǐng)域三通管道的腐蝕失效問題,基于流體動(dòng)力學(xué)的方法指出,流動(dòng)沖刷和電偶腐蝕的協(xié)同作用是三通管道失效的危害源;偶國(guó)富等[14]針對(duì)葉片式靜態(tài)混合器內(nèi)的多相流動(dòng)特性進(jìn)行了數(shù)值模擬,建立了流速與湍流強(qiáng)度的對(duì)應(yīng)關(guān)系;Costa等[15]測(cè)試獲得了90° T型三通內(nèi)的壓降,并揭示了三通相連位置的邊緣效應(yīng)對(duì)流動(dòng)特性的影響;Zhu等[16]基于CFD和流固耦合相結(jié)合的方法,通過求解RANS方程獲得了石油輸運(yùn)過程T型三通內(nèi)壓力和流速的分布規(guī)律。

    為了進(jìn)一步揭示注劑T型管道內(nèi)的氣液運(yùn)動(dòng)規(guī)律,筆者以加氫空冷器上游的工藝注劑T型管為研究對(duì)象,根據(jù)實(shí)際工況(管道結(jié)構(gòu)、多相流介質(zhì))進(jìn)行建模,采用流體動(dòng)力學(xué)軟件,分析注劑T型管道內(nèi)流體動(dòng)力學(xué)參數(shù)的分布規(guī)律,研究注劑T型管道的沖蝕失效形成機(jī)理。本研究有望為加氫空冷器系統(tǒng)的優(yōu)化操作、優(yōu)化運(yùn)行提供指導(dǎo)意義,并為工藝注劑T型管段的設(shè)計(jì)優(yōu)化提供理論支持。

    1 加氫空冷器計(jì)算流體力學(xué)模型

    1.1加氫空冷器工藝流程及注劑T型管結(jié)構(gòu)

    某企業(yè)的加氫空冷器系統(tǒng)的工藝流程示意圖見圖1。加熱后的原料油與催化劑在加氫反應(yīng)器反應(yīng)后生成反應(yīng)流出物從反應(yīng)器底部流出,反應(yīng)流出物依次經(jīng)換熱器、空冷器后進(jìn)入高壓分離器,經(jīng)高壓分離器分離出循環(huán)氫、高分油和含硫污水。其中,循環(huán)氫可返回至反應(yīng)器循環(huán)利用,高分油再進(jìn)入低壓分離器繼續(xù)進(jìn)行三相分離。

    圖1 加氫空冷器系統(tǒng)工藝流程示意圖Fig.1 Process flow diagram of hydrogenation REAC system

    根據(jù)API 932-B標(biāo)準(zhǔn)的規(guī)定,在空冷器前的管道上設(shè)置工藝注劑(注水)點(diǎn)(見圖1中的A)。注劑T型管的幾何結(jié)構(gòu)如圖2所示,其主流管徑(d1)和支流管徑(d2)分別為φ170 mm和φ40 mm。定義主流和支流管徑軸線的交點(diǎn)為坐標(biāo)原點(diǎn)(O),水平向右為x軸正方向,豎直向上為y軸正方向,v1為主流管道多相流流速,v2為支流管道注劑流速。為使進(jìn)出口多相流達(dá)到充分發(fā)展?fàn)顟B(tài),主流管道和支流管道的管道長(zhǎng)度分別延長(zhǎng)至24d1和50d2。

    運(yùn)用ICEM對(duì)注劑T型管進(jìn)行網(wǎng)格劃分;考慮到網(wǎng)格質(zhì)量與計(jì)算結(jié)果的準(zhǔn)確性直接相關(guān),建模過程采用分區(qū)劃分網(wǎng)格的方法,并對(duì)模擬計(jì)算過程進(jìn)行了網(wǎng)格無關(guān)性驗(yàn)證,計(jì)算得到各網(wǎng)格密度下注劑T型管主管的進(jìn)出口壓力差,如表1所示。由表1可知,相對(duì)于8×105的網(wǎng)格密度,隨著網(wǎng)格密度的逐漸增加,進(jìn)出口壓力差的相對(duì)變化值均小于1%,認(rèn)為已經(jīng)達(dá)到網(wǎng)格無關(guān)性要求。最終確定工藝注劑T型管的網(wǎng)格數(shù)量為826300個(gè)。

    圖2 注劑T型管的幾何模型示意Fig.2 Geometric modeling of injection T-tube

    GridnumberPressuredrop/PaRelativevariation/%8×105156.510×105156.10.2612×105155.90.38

    1.2湍流控制方程及湍流模型

    加氫反應(yīng)流出物的多相流介質(zhì)含氣、油、水等組分,工藝注劑沖洗水經(jīng)支流管道與主流管道介質(zhì)混合后,構(gòu)成氣-液多相流,故采用Mixture模型描述主流管道與支流管道內(nèi)流體介質(zhì)的氣-液混合過程,流體運(yùn)動(dòng)的連續(xù)性方程和動(dòng)量方程分別表示為:

    (1)

    (2)

    第α相的體積分?jǐn)?shù)方程:

    (3)

    (4)

    考慮主流管道與支流管道內(nèi)的多相流介質(zhì)混合效果及邊緣漩渦效應(yīng),氣液相間存在圓柱射流狀態(tài),故選擇Realizablek-ε湍流模型對(duì)動(dòng)量方程進(jìn)行封閉求解。湍動(dòng)能和湍流耗散輸運(yùn)方程為:

    (5)

    (6)

    式中:k為湍動(dòng)能,ε為耗散率,μi為i方向的速度分量,xi、xj分別為i、j方向上的坐標(biāo)分量,μt為湍流黏度,GK為速度引起的湍動(dòng)能生成項(xiàng),Gb為浮力引起的湍動(dòng)能生成項(xiàng),YM表示可壓縮湍流脈動(dòng)膨脹對(duì)總的耗散率影響;αk、αε分別為湍動(dòng)能及其耗散率的湍流普朗特?cái)?shù),分別取值1.0、1.2;C1ε、C2ε為常數(shù),取值1.44、1.9。C3ε亦為常數(shù),當(dāng)浮力應(yīng)力層中速度方向與重力方向相同時(shí),其值為1;當(dāng)浮力應(yīng)力層中速度方向與重力方向垂直時(shí),其值為0。

    1.3邊界條件及求解格式

    計(jì)算域采用速度作為進(jìn)口參數(shù),氣、油、水為主流管道入口介質(zhì),水相為支流管道入口介質(zhì)。出口邊界為自由出流,按照充分發(fā)展管流條件處理,變量在出口截面的法向梯度上為零,為此模擬過程中將主流管道和支流管道的管路加長(zhǎng),以保證流動(dòng)充分發(fā)展。壁面邊界采用無滑移邊界條件,近壁網(wǎng)格采用標(biāo)準(zhǔn)壁面函數(shù)進(jìn)行處理。求解過程使用有限體積法實(shí)現(xiàn)計(jì)算區(qū)域和控制方程的離散化,體積相分率、湍動(dòng)能、動(dòng)量和湍流耗散率采用一階迎風(fēng)格式進(jìn)行離散,壓力項(xiàng)采用Standard格式,壓力-速度耦合方程的求解采用SIMPLE方法。通過選擇加氫空冷器典型工況和操作參數(shù),運(yùn)用工藝仿真軟件Aspen獲得注劑T型管主管入口的多相流介質(zhì)組成及物性參數(shù),其結(jié)果見表2。

    表2 加氫空冷器注劑T型管道進(jìn)口物性參數(shù)Table 2 Inlet physical characteristics of injection T-tube in the hydrogenation air cooler

    2 注劑T型管數(shù)值模擬結(jié)果與討論

    考慮注劑T型管道內(nèi)多相流混合過程的復(fù)雜性,數(shù)值模擬過程作如下基本假定:1)流場(chǎng)內(nèi)輸運(yùn)的流體為牛頓流體;2)注劑T型管道內(nèi)部溫度恒定,不考慮多相流介質(zhì)的傳質(zhì)傳熱。采用流體動(dòng)力學(xué)軟件進(jìn)行數(shù)值計(jì)算,并對(duì)數(shù)值計(jì)算結(jié)果進(jìn)行網(wǎng)格自適應(yīng)處理,模擬結(jié)果及分析如下。

    2.1注劑T型管內(nèi)流動(dòng)特性的分布規(guī)律

    圖3為三維坐標(biāo)系中z=0 m的平面內(nèi)90° 注劑T型管內(nèi)的水相分率及流線圖。支流管道注水之前,主流管道內(nèi)的介質(zhì)基本被氣相所占據(jù),在主流管道油、氣相的橫向剪切作用下,部分水相被湍流核心區(qū)域流速相對(duì)較高的氣相拖曳輸運(yùn),垂直向下射入的水相區(qū)域出現(xiàn)偏移,其最高位置出現(xiàn)在支流管道左側(cè)區(qū)域。由圖3可以清楚的看出,在x=-0.2 m至x=-0.7 m區(qū)域的主流管道底部流線較為集中,隨著主流管道下游流體的充分發(fā)展,氣相和水相充分混合,靠近出口邊界區(qū)域水相分率逐漸降低至0.05左右。與z=0平面內(nèi)湍動(dòng)能的分布規(guī)律(見圖4)對(duì)比分析可知,在注水點(diǎn)的下游,主流管道內(nèi)的氣-液混合相湍動(dòng)強(qiáng)度顯著增加,其中湍動(dòng)能較高的區(qū)域主要對(duì)應(yīng)于水相分率較高和流線發(fā)生突變的部位。

    圖5為沿x軸負(fù)方向不同截面的T型管內(nèi)水相分率分布圖。由圖5可知,在x=-0.10 m的橫截面上,由于主流管道湍流核心區(qū)域氣相速度的橫向剪切及豎直向下注水速度的聯(lián)合作用,水相分率分布較高的區(qū)域基本上位于主流管道的中心位置,其分布形態(tài)為沿x-y平面對(duì)稱分布;在x=-0.25 m的橫截面上,水相在橫截面上的擴(kuò)散效應(yīng)顯著增強(qiáng),較x=-0.10 m截面的水相分率相比,水相分布區(qū)域擴(kuò)大,且漸漸靠近至管道底部區(qū)域,但其峰值已由0.25逐漸降低至0.15左右;在x為-0.50、-0.75、-1.00 m區(qū)域,水相分率均為沿x-y平面對(duì)稱分布,但隨著氣-液相交界面相間質(zhì)量和動(dòng)量的交換,受相間曳力的影響,油-氣-水多相流的混合效應(yīng)明顯,相間速度梯度逐漸降低,水相分率的峰值與前者相比相對(duì)減小且呈現(xiàn)先向橫截面周向擴(kuò)散然后又向管道底部匯聚的趨勢(shì)。

    圖3 z=0 m截面T型管內(nèi)水相分率與流線分布Fig.3 Water phase fraction distribution andstreamline in section z=0 m

    圖4 z=0 m截面T型管內(nèi)湍動(dòng)能與流線分布Fig.4 Turbulent kinetic energy distribution andstreamline in section z=0 m

    圖5 注劑T型管主流管道不同橫截面水相分率分布Fig.5 Water phase fraction contour in different cross-sectionsof injection T-tube mainstream area

    為對(duì)比呈現(xiàn)注劑T型管內(nèi)的流動(dòng)情況,圖6給出了上述5個(gè)橫截面以及x=-1.25 m處的水相分率等值線圖。由圖6(a)可知,在x=-0.1 m橫截面上,水相集中在面域的核心區(qū)域,原因在于:油氣相與水相間速度梯度相對(duì)較高,多相流混合效應(yīng)不顯著;核心區(qū)域邊緣部分由于水-氣相界面的存在,水相分率由核心區(qū)域的最高值0.25降低至0.09左右;受主流和支流管道結(jié)構(gòu)特征的影響,水相沿z=0平面對(duì)稱分布。由圖6(b)可知,在x=-0.25 m 橫截面上,由于水相從支流管道進(jìn)入時(shí)具有噴射作用,以及在主流管道介質(zhì)的橫向剪切作用,水相分率分布較高的區(qū)域逐漸向橫截面底部遷移,同時(shí)氣液相混合效應(yīng)相應(yīng)提高,水相分率的峰值降低至0.14左右,其中管道底部近壁面水相分率約為0.12,分布區(qū)域約占管道橫截面周長(zhǎng)的1/3。由圖6(c)可知,在x=-0.50 m截面區(qū)域,水相分率已從湍流的核心區(qū)域逐漸向整個(gè)面域內(nèi)擴(kuò)散,峰值進(jìn)一步降低至0.10左右,同時(shí)水相分率沿管道周向的分布區(qū)域有所擴(kuò)大,占管道橫截面周長(zhǎng)的3/4以上。由圖6(d)、(e)、(f)可知,隨著多相流介質(zhì)向管道下游流動(dòng),水相的噴射作用減弱,徑向速度逐步減小,水相分布漸趨均勻,峰值在0.07左右,標(biāo)志著氣液相進(jìn)入充分發(fā)展區(qū)域。

    圖6 注劑T型管主流管道不同橫截面水相分率等值線分布Fig.6 Water phase fraction contour in different cross-sections of injection T-tube mainstream area(a) x=-0.10 m; (b) x=-0.25 m; (c) x=-0.50 m; (d) x=-0.75 ml; (e) x=-1.00 m; (f) x=-1.25 m

    對(duì)比分析圖6(a)~圖6(f)可知,在x為-0.25~-0.50 m區(qū)域內(nèi),主流管道底部的水相分率相對(duì)較高,若管內(nèi)存在結(jié)晶反應(yīng)析出的NH4Cl晶體,銨鹽顆粒易溶解于水并對(duì)管壁造成腐蝕。

    2.290°注劑T型管內(nèi)流動(dòng)特性變化規(guī)律

    為豐富流動(dòng)的研究,進(jìn)而選擇瞬態(tài)求解器對(duì)90° 注劑T型管內(nèi)流動(dòng)特性變化規(guī)律進(jìn)行更進(jìn)一步的研究,設(shè)置時(shí)間步長(zhǎng)為0.002 s,時(shí)間步最大迭代20次,瞬態(tài)研究結(jié)果經(jīng)后處理后,得到流動(dòng)的不同時(shí)刻主流管道下游不同橫截面上的水相分率云圖,如圖7所示,且圖7中各截面位置與圖5各截面位置相同。由圖7可知,4s過后,在截面x=-0.1 m和截面x=-0.25 m處水相分布已經(jīng)穩(wěn)定,與穩(wěn)態(tài)計(jì)算所得的結(jié)果已經(jīng)一致;6 s過后,對(duì)比圖7(d)與圖5可知,整個(gè)注劑T型管的主流管道內(nèi)水相分布已和穩(wěn)態(tài)計(jì)算的水相分布一致。

    圖7 注劑T型管主流管道不同橫截面水相分率發(fā)展過程Fig.7 Development of water phase fraction contour in different cross-sections of injection T-tube mainstream areaTime step/s: (a) 2; (b) 3; (c) 4; (d) 6

    2.3注水的入射角度對(duì)流動(dòng)特性的影響

    圖8為入射角度為90°時(shí),注劑T型管主流管道底部近壁面的水相分率及剪切應(yīng)力分布圖。由圖8 可知,壁面剪切應(yīng)力分布最大的區(qū)域分布在區(qū)域①內(nèi),約為x為-0.15~-0.75 m之間,峰值約0.9 Pa;而沿x軸正方向的管道近壁面區(qū)域,由于主流管道多相流中氣相占比達(dá)98%,壁面剪切應(yīng)力僅為0.2 Pa。主流管道底部近壁面水相分率分布最大的區(qū)域?yàn)閰^(qū)域②,約為x為-0.12~-0.26 m之間;其中近壁面剪切應(yīng)力分布最大的區(qū)域①與水相分率最大的區(qū)域②交集為x為-0.15~-0.26 m之間,在該區(qū)域內(nèi),水相分率高,易使結(jié)晶的銨鹽溶解形成腐蝕溶液對(duì)管壁形成腐蝕,剪切應(yīng)力大則易引發(fā)腐蝕基礎(chǔ)上的沖蝕失效[18]。

    圖8 主流管道下部壁面剪切應(yīng)力和水相分率分布(入射角90°)Fig.8 Wall shear stress and water phase fraction distributionat the bottom of main pipe with merging angle of 90°

    圖9為入射角度為75°時(shí),注劑T型管主流管道底部近壁面水相分率及剪切應(yīng)力分布。由圖9可知,隨著水相注入,受氣、液相混合黏度及相間速度梯度的影響,主流管道近壁面剪切應(yīng)力較注水點(diǎn)前增加了2倍以上,剪切應(yīng)力較大的區(qū)域分布在x為-0.25~-1.9 m之間(見圖9中的區(qū)域③),其峰值約為0.85 Pa;相應(yīng)地,主流管道近壁面的水相分率較注水點(diǎn)前增幅更為顯著,其最大值出現(xiàn)在區(qū)域④,約在x為-0.7~-1.0 m之間,與入射角度為90°的注劑T型管相比,由于水相射流沿主流管道的徑向速度分量減小,水相分率在主流管道下部近壁面的峰值由x=-0.25 m逐漸后移至x=-0.78 m 處。

    圖9 主流管道下部壁面剪切應(yīng)力和水相分率分布(入射角75°)Fig.9 Wall shear stress and water phase fraction distributionat the bottom of main pipe with merging angle of 75°

    圖10為入射角度為60°時(shí),注劑T型管主流管道底部近壁面水相分率及剪切應(yīng)力分布。由圖10可知,在注水點(diǎn)后的區(qū)域⑤,對(duì)應(yīng)于x為-0.35~-0.80 m 之間,出現(xiàn)剪切應(yīng)力峰值(1.18 Pa);水相分率較高的位置位于區(qū)域⑥。其中對(duì)應(yīng)水相分率和剪切應(yīng)力均較大的區(qū)域約為x為-0.35~-0.80 m。對(duì)比圖8可知,隨著入射角度的減小,剪切應(yīng)力峰值出現(xiàn)逐漸向主流管道下游移動(dòng)的趨勢(shì)。

    當(dāng)入射角為45°和30°時(shí),主流管道下部壁面剪切應(yīng)力和水相分率的分布規(guī)律分別如圖11和圖12所示。由圖11可知,水相分率的峰值約為0.28,約位于x=-0.13 m截面,剪切應(yīng)力的峰值約為2.0 Pa,約位于x=-0.24 m截面;當(dāng)入射角進(jìn)一步降低至30°時(shí)(見圖12),主流管道近壁面剪切應(yīng)力的峰值增加至2.2 Pa。對(duì)比圖11、圖12可知,水相分率和剪切應(yīng)力的分布趨勢(shì)、峰值出現(xiàn)的位置基本相似。與高角度入射角相比,低角度入射時(shí)主流管道底部水相分率和剪切應(yīng)力峰值出現(xiàn)的位置基本重合,當(dāng)存在NH4Cl顆粒時(shí),更易出現(xiàn)溶液腐蝕基礎(chǔ)上的沖蝕失效。因此,應(yīng)對(duì)x為-0.13~-0.24 m 之間的主流管道底部區(qū)域重點(diǎn)關(guān)注,定期監(jiān)測(cè)其剩余壁厚,防止突然穿孔泄漏引發(fā)非計(jì)劃停工事故。

    圖10 主流管道下部壁面剪切應(yīng)力和水相分率分布(入射角60°)Fig.10 Wall shear stress and water phase fraction distributionat the bottom of main pipe with merging angle of 60°

    圖11 主流管道下部壁面剪切應(yīng)力和水相分率分布(入射角45°)Fig.11 Wall shear stress and water phase fraction distributionat the bottom of main pipe with merging angle of 45°

    圖12 主流管道下部壁面剪切應(yīng)力和水相分率分布(入射角30°)Fig.12 Wall shear stress and water phase fraction distributionat the bottom of main pipe with merging angle of 30°

    通過分析以上5種不同入射角注劑T型管的流體動(dòng)力學(xué)參數(shù)分布情況,可以幫助選取最優(yōu)的入射角度,優(yōu)化注劑T型管結(jié)構(gòu)的設(shè)計(jì),使工藝設(shè)備能夠長(zhǎng)周期、穩(wěn)定、安全運(yùn)行。

    2.4入射角為90°時(shí)注劑T型管超聲波測(cè)厚驗(yàn)證

    含氯原料油加工過程中,出現(xiàn)工藝注劑T型管穿孔泄漏的現(xiàn)象,采用ECHOMETER1076超聲波測(cè)厚儀(精度為0.01 mm)對(duì)實(shí)際運(yùn)行工況下的入射角為90°時(shí)注劑T型管進(jìn)行剩余壁厚檢測(cè)。測(cè)厚點(diǎn)

    位置劃分如圖13所示,將主流管段自左向右間隔20 mm劃線,沿周向進(jìn)行角度劃分定位,其中165°~195°之間標(biāo)記為主流管道的底部區(qū)域。測(cè)量經(jīng)線和緯線交點(diǎn)的剩余壁厚,檢測(cè)結(jié)果見圖14。

    圖13 泄漏區(qū)域測(cè)厚方案示意Fig.13 Schematic diagram of thickness measurementin leakage area

    圖14 泄漏區(qū)域注劑T型管剩余壁厚分布(入射角為90°)Fig.14 Distribution of remaining wall thickness in leakage area with merging angle of 90°

    超聲波測(cè)厚數(shù)據(jù)表明:穿孔泄漏的位置位于經(jīng)向角度165°~195° 間,位于90° 注劑T型管主流管道的底部,軸向位置位于x=-0.10 m±0.02 m,且穿孔泄漏位置的兩側(cè)均有不同程度的壁厚減薄。由數(shù)值模擬獲得的水相分率及剪切應(yīng)力分布規(guī)律(圖8)對(duì)比分析可知,水相分率及剪切應(yīng)力較大的位置,即x為-0.15~-0.26 m之間的區(qū)域,這一數(shù)值模擬結(jié)果與實(shí)際管道發(fā)生穿孔泄漏的位置基本吻合。

    3 結(jié) 論

    (1)建立了氣-液流動(dòng)的數(shù)理模型,基于數(shù)值模擬獲得了注劑T型管內(nèi)的氣-液流動(dòng)特性,獲得了氣-液相的流體動(dòng)力學(xué)參數(shù)分布規(guī)律,其中水相分率和剪切應(yīng)力大的位置與注劑T型管穿孔泄漏的位置基本吻合,驗(yàn)證了加氫空冷器注劑T型管內(nèi)氣液流動(dòng)特性預(yù)測(cè)方法的正確性。

    (2)工藝注水點(diǎn)下游區(qū)域,在較高氣相速度橫向剪切作用下,在x為-0.15~-0.26 m橫截面,90°注劑T型管主流管道下部近壁面區(qū)域,水相分率較高,若存在NH4Cl結(jié)晶顆粒,易形成腐蝕性水溶液,在近壁面高剪切應(yīng)力的作用下,易對(duì)該區(qū)域的主流管道壁面產(chǎn)生沖蝕失效。

    (3)與90°注劑T型管相比,入射角60°、75°注劑T型管內(nèi)主流管道下部近壁面的最大剪切應(yīng)力峰值均向流體下游移動(dòng);入射角30°、45°注劑T型管內(nèi)的氣液流動(dòng)特性基本相似,水相分率和剪切應(yīng)力最大區(qū)域基本重合,在x為-0.13~-0.24 m間的主流管道下部區(qū)域存在沖蝕失效風(fēng)險(xiǎn)。

    [1] PALMASP. 25 years of development[J].Hydrocarbon Engineering, 2009, 14(6): 53-59.

    [2] HARVEY C, SINGH A. Mitigate failure for reactor effluent air coolers[J].Hydrocarbon Processing, 1999, 10(3): 59-72.

    [3] ZHENG Z J, OU G F, YE H J, et al. Investigation on the deposition failure of a reactor effluent air cooler in hydrocracking unit[J].Engineering Failure Analysis, 2016, 68: 52-63.

    [4] SUN L, ZHU M, OU G F, et al. Corrosion investigation of the inlet section of REAC pipes in the refinery[J].Engineering Failure Analysis, 2016, 66: 468-478.

    [5] 偶國(guó)富, 裘杰, 朱祖超, 等. 異徑管沖蝕失效的流固耦合數(shù)值模擬[J].力學(xué)學(xué)報(bào), 2010, 42(2): 197-204. (OU Guofu, QIU Jie, ZHU Zuchao, et al. Fluid-structure interactional numerical simulation of erosion-corrosion failure of reducer in multiphase flow[J].Chinese Journal of Theoretical and Applied Mechanics, 2010, 42(2): 197-204.)

    [6] JIN H Z, OU G F, WANG Y P, et al. Failure analysis and structure optimization of hydrogenation air-cooler system based on imbalanced degree[C]//ASME Pressure Vessels and Piping Conference, Prague Czech, 2009: 439-448.

    [7] 譚金龍, 夏翔鳴, 胡傳清, 等. 加氫裂化裝置高壓空氣冷卻器的腐蝕失效分析[J].石油化工腐蝕與防護(hù), 2009, 26(2): 52-57. (TAN Jinlong, XIA Xiangming, HU Chuanqing, et al. Corrosion failure analysis of high-pressure air cooler of hydrocracker[J].Corrosion & Protection in Petrochemical Industry, 2009, 26(2): 52-57.)

    [8] 金浩哲, 偶國(guó)富, 王寬心, 等. 加氫處理系統(tǒng)NH4Cl結(jié)晶沉積預(yù)測(cè)及優(yōu)化防控[J].石油學(xué)報(bào)(石油加工), 2014, 30(4): 662-667. (JIN Haozhe, OU Guofu, WANG Kuanxin, et al. Prediction of NH4Cl crystalline deposition and optimized prevention method in hydrotreating unit system[J].Acta Petrolei Sinica (Petroleum Processing Section), 2014, 30(4): 662-667.)

    [9] 金浩哲, 王寬心, 偶國(guó)富, 等. 基于Aspen Plus的加氫反應(yīng)流出物銨鹽結(jié)晶速率計(jì)算模型與分析[J].石油學(xué)報(bào)(石油加工), 2015, 31(6): 1444-1449. (JIN Haozhe, WANG Kuanxin, OU Guofu, et al. Calculation model and analysis of ammonium salt crystallization rate in hydrogenation reaction effluent based on aspen plus[J].Acta Petrolei Sinica (Petroleum Processing Section), 2015, 31(6): 1444-1449.)

    [10] API Publication 932-A. A study of corrosion in hydroprocessing reactor effluent air coolers system[S]. 2002.

    [11] API Publication 932-B. Design, materials, fabrication, operation and inspection guideline for corrosion control in hydroprocessing reactor effluent air cooler (REAC) system[S].Second Edition, 2004.

    [12] JACOBS G E, SHARGAY C A, CABRER A J, et al. Single vs multiple injection points for REAC wash water system interesting process simulation results[C]//California USA, NACE CORROSION, 2016, No.92698.

    [13] WANG Z, LIU B, YANG Y, et al. Experimental and numerical studies on corrosion failure of a three-limb pipe in natural gas field[J].Engineering Failure Analysis, 2016, 62: 21-38.

    [14] 偶國(guó)富, 鄭智劍, 金浩哲. 葉片式靜態(tài)混合器多相流動(dòng)特性的數(shù)值分析[J].高校化學(xué)工程學(xué)報(bào), 2016, 30(1): 40-47. (OU Guofu, ZHENG Zhijian, JIN Haozhe. Numerical analysis on the multiphase flow characteristics of static blade mixer[J].Journal of Chemical Engineering of Chinese Universities, 2016, 30(1): 40-47.)

    [15] COSTA N P, MAIA R, PROENCA M F, et al. Edge effects on the flow characteristics in a 90deg tee junction[J].Journal of Fluids Engineering, 2006, 128(6): 1204-1217.

    [16] ZHU H J, ZHANG W L,FENG G, et al. Fluid-structure interaction computational analysis of flow field, shear stress distribution and deformation of three-limb pipe[J].Engineering Failure Analysis, 2014, 42: 252-262.

    [17] BRACKBILL J U, KOTHE D B, ZEMACH C. A continuum method for modeling surface tension[J].Journal of Computational Physics, 1992, 100(2): 335-354.

    [18] 偶國(guó)富, 許根富, 朱祖超, 等. 彎管沖蝕失效流固耦合機(jī)理及數(shù)值模擬[J].機(jī)械工程學(xué)報(bào), 2009, 45(11): 119-124, 132. (OU Guofu, XU Genfu, ZHU Zuchao, et al. Fluid-structure interaction mechanism and numerical simulation of elbow erosion failure[J].Journal of Mechanical Engineering, 2009, 45(11): 119-124, 132.)

    NumericalSimulationofGas-LiquidPhaseFlowCharacteristicofInjectionT-tubeinHydrogenationReactorEffluentAirCoolerSystems

    JIN Haozhe1, LIU Wenwen1, OU Guofu1, CHEN Xiaoping1, LI Pengxuan2

    (1.TheInstituteofFlow-InducedCorrosion,ZhejiangSci-TechUniversity,Hangzhou310018,China,2.HangzhouSpecialEquipmentInspectionandResearchInstitute,Hangzhou310003,China)

    2016-10-14

    國(guó)家自然科學(xué)基金委員會(huì)-神華集團(tuán)有限公司煤炭聯(lián)合基金(U1361107)、浙江省自然科學(xué)基金(LY17E060008)、浙江省科技計(jì)劃項(xiàng)目(2017C01070)、浙江省公益技術(shù)應(yīng)用研究計(jì)劃項(xiàng)目(2015C31013)和浙江理工大學(xué)研究生創(chuàng)新研究項(xiàng)目(YCX16029)資助

    金浩哲,男,副教授,博士,主要從事多相流沉積與沖蝕的研究;E-mail:haozhe2007@163.com

    1001-8719(2017)05-0910-09

    TQ022.4

    A

    10.3969/j.issn.1001-8719.2017.05.013

    猜你喜歡
    型管分率剪切應(yīng)力
    量率對(duì)應(yīng) 解決問題
    b型管板與筒體溫差應(yīng)力的分析計(jì)算和評(píng)定
    化工管理(2021年7期)2021-05-13 00:46:04
    內(nèi)置加勁環(huán)T型管節(jié)點(diǎn)抗沖擊承載力計(jì)算
    Study on the Letter Words and Letter-word Phrases under Language Contact Theory
    Montgomery T型管置入術(shù)治療聲門下氣管狹窄失敗三例報(bào)道并文獻(xiàn)復(fù)習(xí)
    心瓣瓣膜區(qū)流場(chǎng)中湍流剪切應(yīng)力對(duì)瓣膜損害的研究進(jìn)展
    解分?jǐn)?shù)問題例談
    分?jǐn)?shù)應(yīng)用題常見錯(cuò)例剖析
    剪切應(yīng)力對(duì)聚乳酸結(jié)晶性能的影響
    利用分率巧解題
    少妇丰满av| 俄罗斯特黄特色一大片| 免费观看的影片在线观看| 夜夜夜夜夜久久久久| 国产高清三级在线| 激情在线观看视频在线高清| 亚洲欧美日韩无卡精品| 国产精品精品国产色婷婷| 久久九九热精品免费| 美女高潮的动态| 99热精品在线国产| 三级国产精品欧美在线观看 | 香蕉国产在线看| 色av中文字幕| 亚洲欧美日韩东京热| 久久久久久人人人人人| 美女黄网站色视频| 亚洲国产精品成人综合色| 精品一区二区三区四区五区乱码| 国产激情欧美一区二区| 亚洲电影在线观看av| 男人舔女人下体高潮全视频| 国产高清有码在线观看视频| 老司机午夜福利在线观看视频| 国产激情偷乱视频一区二区| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品国产精品久久久不卡| 最近最新免费中文字幕在线| 婷婷丁香在线五月| 香蕉久久夜色| 国产三级在线视频| 偷拍熟女少妇极品色| 午夜精品久久久久久毛片777| 成人18禁在线播放| 国产精品98久久久久久宅男小说| 国产主播在线观看一区二区| 免费搜索国产男女视频| 亚洲avbb在线观看| 国产精品国产高清国产av| 国产精品久久久久久亚洲av鲁大| 90打野战视频偷拍视频| 国产成人一区二区三区免费视频网站| 啦啦啦韩国在线观看视频| 网址你懂的国产日韩在线| 国产成人系列免费观看| 欧美在线一区亚洲| 久久久久精品国产欧美久久久| 免费无遮挡裸体视频| 欧美在线黄色| 色精品久久人妻99蜜桃| 日韩欧美国产在线观看| 国产一级毛片七仙女欲春2| 18禁国产床啪视频网站| 后天国语完整版免费观看| 国产又黄又爽又无遮挡在线| 国产一级毛片七仙女欲春2| 成人亚洲精品av一区二区| 久久精品亚洲精品国产色婷小说| 欧美在线黄色| 国产精品女同一区二区软件 | 色在线成人网| 国产精品久久久久久人妻精品电影| 精品久久蜜臀av无| 网址你懂的国产日韩在线| 天天躁狠狠躁夜夜躁狠狠躁| 18禁美女被吸乳视频| 亚洲一区高清亚洲精品| 亚洲av五月六月丁香网| 他把我摸到了高潮在线观看| 欧美乱妇无乱码| 久久久久久久久免费视频了| 伊人久久大香线蕉亚洲五| 久久香蕉国产精品| 日本撒尿小便嘘嘘汇集6| 男女那种视频在线观看| 在线国产一区二区在线| 国产男靠女视频免费网站| 欧美黑人欧美精品刺激| 18美女黄网站色大片免费观看| aaaaa片日本免费| 亚洲国产精品久久男人天堂| 757午夜福利合集在线观看| 色哟哟哟哟哟哟| 欧美av亚洲av综合av国产av| 欧美黑人巨大hd| 亚洲人成伊人成综合网2020| 日日夜夜操网爽| 国产精品1区2区在线观看.| 久久欧美精品欧美久久欧美| 国产成年人精品一区二区| 性色av乱码一区二区三区2| 国产成人精品无人区| 夜夜爽天天搞| 国产精品美女特级片免费视频播放器 | 黄色丝袜av网址大全| 日韩欧美在线乱码| 久久久久久人人人人人| 老汉色av国产亚洲站长工具| 变态另类丝袜制服| 两性夫妻黄色片| 99久久综合精品五月天人人| 国产精品亚洲美女久久久| 91av网一区二区| 亚洲无线在线观看| 熟女人妻精品中文字幕| 窝窝影院91人妻| 女同久久另类99精品国产91| 国产蜜桃级精品一区二区三区| 国内毛片毛片毛片毛片毛片| 国产日本99.免费观看| 国产探花在线观看一区二区| 免费无遮挡裸体视频| 国产乱人视频| 91av网站免费观看| 精品日产1卡2卡| 少妇的逼水好多| 欧美三级亚洲精品| 久久久久久久久久黄片| 欧美成人一区二区免费高清观看 | 美女大奶头视频| 哪里可以看免费的av片| 国产欧美日韩一区二区精品| 最近最新中文字幕大全电影3| 少妇熟女aⅴ在线视频| 日韩人妻高清精品专区| 成年女人永久免费观看视频| 麻豆国产av国片精品| 天堂影院成人在线观看| 国产精品99久久99久久久不卡| 日韩欧美精品v在线| 久久久国产精品麻豆| 国产亚洲精品久久久com| 老司机福利观看| 九色国产91popny在线| 午夜视频精品福利| 丰满人妻一区二区三区视频av | 他把我摸到了高潮在线观看| 成人无遮挡网站| 日韩大尺度精品在线看网址| 国产精品综合久久久久久久免费| 亚洲av日韩精品久久久久久密| 男人和女人高潮做爰伦理| 天堂影院成人在线观看| 麻豆久久精品国产亚洲av| 午夜免费观看网址| 最近最新免费中文字幕在线| 最近最新中文字幕大全电影3| 亚洲av五月六月丁香网| 免费av毛片视频| 亚洲精品中文字幕一二三四区| 99热6这里只有精品| 搡老熟女国产l中国老女人| 国产乱人视频| 女生性感内裤真人,穿戴方法视频| 男插女下体视频免费在线播放| www.精华液| 成人国产综合亚洲| 国内揄拍国产精品人妻在线| 日韩av在线大香蕉| 美女扒开内裤让男人捅视频| 又紧又爽又黄一区二区| 神马国产精品三级电影在线观看| 成人午夜高清在线视频| 精品久久久久久成人av| 手机成人av网站| 欧美在线黄色| 国产精品爽爽va在线观看网站| 91在线精品国自产拍蜜月 | 色吧在线观看| 视频区欧美日本亚洲| 视频区欧美日本亚洲| 成人特级av手机在线观看| 90打野战视频偷拍视频| 国产伦人伦偷精品视频| 黄色日韩在线| 99国产精品一区二区三区| 99久久99久久久精品蜜桃| 搞女人的毛片| 网址你懂的国产日韩在线| 日韩中文字幕欧美一区二区| 夜夜看夜夜爽夜夜摸| av天堂在线播放| 综合色av麻豆| xxxwww97欧美| 国语自产精品视频在线第100页| avwww免费| 成人午夜高清在线视频| 欧美另类亚洲清纯唯美| 欧美精品啪啪一区二区三区| 日韩成人在线观看一区二区三区| 叶爱在线成人免费视频播放| 熟女少妇亚洲综合色aaa.| 激情在线观看视频在线高清| 他把我摸到了高潮在线观看| 亚洲精品久久国产高清桃花| 夜夜躁狠狠躁天天躁| 亚洲人成网站在线播放欧美日韩| 女生性感内裤真人,穿戴方法视频| 嫩草影院入口| 国产激情偷乱视频一区二区| 亚洲精品一区av在线观看| 国产精品女同一区二区软件 | 免费无遮挡裸体视频| 欧美在线黄色| a级毛片在线看网站| 美女扒开内裤让男人捅视频| 国产亚洲精品久久久久久毛片| 久久精品国产清高在天天线| 亚洲国产色片| 一级毛片高清免费大全| 午夜福利18| 搡老妇女老女人老熟妇| 级片在线观看| 亚洲国产高清在线一区二区三| 国产精品av久久久久免费| 久久久久性生活片| 久久精品亚洲精品国产色婷小说| 亚洲国产精品999在线| 最近最新中文字幕大全免费视频| 黄色丝袜av网址大全| 亚洲人成网站在线播放欧美日韩| 成年女人永久免费观看视频| 亚洲第一欧美日韩一区二区三区| 精品福利观看| 热99在线观看视频| АⅤ资源中文在线天堂| 亚洲av日韩精品久久久久久密| 一级毛片高清免费大全| 久久香蕉精品热| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av香蕉五月| 俄罗斯特黄特色一大片| 两性夫妻黄色片| 人妻久久中文字幕网| 精品久久久久久久人妻蜜臀av| 在线观看免费视频日本深夜| 久久久久久人人人人人| 国产精品久久久久久精品电影| 狠狠狠狠99中文字幕| 精品熟女少妇八av免费久了| 午夜久久久久精精品| www日本在线高清视频| 一区二区三区高清视频在线| 久久精品aⅴ一区二区三区四区| 亚洲精品色激情综合| 色综合站精品国产| 男女视频在线观看网站免费| av女优亚洲男人天堂 | 黑人巨大精品欧美一区二区mp4| 精品国产美女av久久久久小说| 中国美女看黄片| 日本一本二区三区精品| 日本黄大片高清| 噜噜噜噜噜久久久久久91| 国产单亲对白刺激| 国产乱人伦免费视频| 国产精品免费一区二区三区在线| 老司机福利观看| 男人舔女人的私密视频| 好男人在线观看高清免费视频| 免费无遮挡裸体视频| 色综合站精品国产| 99热这里只有是精品50| 在线永久观看黄色视频| 亚洲国产中文字幕在线视频| 在线视频色国产色| 热99re8久久精品国产| 久久国产精品影院| 啪啪无遮挡十八禁网站| 国产精品亚洲美女久久久| 久久久水蜜桃国产精品网| 色尼玛亚洲综合影院| 精品久久久久久,| 成人av在线播放网站| 欧美日韩综合久久久久久 | 丰满人妻熟妇乱又伦精品不卡| 午夜免费成人在线视频| 国产麻豆成人av免费视频| 又紧又爽又黄一区二区| 日本 欧美在线| 国产精品99久久久久久久久| 国产精品av久久久久免费| 亚洲在线观看片| 国内揄拍国产精品人妻在线| 欧美一级毛片孕妇| 日韩精品青青久久久久久| 久久九九热精品免费| 国产精品99久久99久久久不卡| 国产精品女同一区二区软件 | 少妇裸体淫交视频免费看高清| 国产探花在线观看一区二区| 91麻豆精品激情在线观看国产| www国产在线视频色| 欧美+亚洲+日韩+国产| 国产亚洲精品久久久com| 亚洲色图av天堂| 黄频高清免费视频| 欧美乱码精品一区二区三区| 亚洲国产欧洲综合997久久,| 国产高清三级在线| 久久久成人免费电影| 国产精品亚洲av一区麻豆| 国产精品综合久久久久久久免费| 久9热在线精品视频| 国产毛片a区久久久久| 亚洲中文字幕日韩| 国产欧美日韩一区二区三| 久99久视频精品免费| 亚洲人成网站在线播放欧美日韩| 天天躁日日操中文字幕| 欧美不卡视频在线免费观看| 国产成人精品久久二区二区91| 精品久久久久久久人妻蜜臀av| 欧美黄色片欧美黄色片| 国产高清三级在线| 欧美黄色片欧美黄色片| 久久久久久久久免费视频了| 美女高潮喷水抽搐中文字幕| 最近最新中文字幕大全电影3| 午夜福利欧美成人| 精品午夜福利视频在线观看一区| 亚洲国产精品999在线| 日本撒尿小便嘘嘘汇集6| 久久久久久人人人人人| 麻豆成人午夜福利视频| 一区二区三区激情视频| 欧美成人一区二区免费高清观看 | 这个男人来自地球电影免费观看| av福利片在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产成+人综合+亚洲专区| 精品久久久久久久毛片微露脸| 国产av麻豆久久久久久久| 99精品欧美一区二区三区四区| 国产精品 欧美亚洲| 五月伊人婷婷丁香| 婷婷丁香在线五月| 国产一区在线观看成人免费| 欧美极品一区二区三区四区| 色噜噜av男人的天堂激情| 久久精品国产99精品国产亚洲性色| 日韩欧美国产一区二区入口| 亚洲人成伊人成综合网2020| 免费看a级黄色片| 国产精品美女特级片免费视频播放器 | 欧美日韩乱码在线| 欧美乱码精品一区二区三区| 最近视频中文字幕2019在线8| 一区二区三区国产精品乱码| 亚洲精品国产精品久久久不卡| 午夜精品一区二区三区免费看| 在线免费观看的www视频| 色尼玛亚洲综合影院| 少妇丰满av| 亚洲avbb在线观看| 最新中文字幕久久久久 | 999久久久国产精品视频| 青草久久国产| 悠悠久久av| 麻豆国产av国片精品| 国产一级毛片七仙女欲春2| 禁无遮挡网站| 1000部很黄的大片| 无人区码免费观看不卡| 欧美又色又爽又黄视频| 桃色一区二区三区在线观看| 亚洲成人精品中文字幕电影| 欧美成人免费av一区二区三区| www日本黄色视频网| 91麻豆av在线| 九九久久精品国产亚洲av麻豆 | 久久国产乱子伦精品免费另类| 亚洲人成伊人成综合网2020| 亚洲av成人不卡在线观看播放网| 桃色一区二区三区在线观看| 97人妻精品一区二区三区麻豆| 少妇的丰满在线观看| 成人三级黄色视频| 亚洲国产中文字幕在线视频| 国产精品av视频在线免费观看| 99热这里只有精品一区 | 99精品久久久久人妻精品| 免费在线观看影片大全网站| 99精品久久久久人妻精品| ponron亚洲| 欧美日韩福利视频一区二区| 波多野结衣高清无吗| 国内精品久久久久久久电影| av天堂在线播放| 国产视频内射| 亚洲乱码一区二区免费版| 久久精品国产99精品国产亚洲性色| 91av网一区二区| 亚洲av电影不卡..在线观看| 久久婷婷人人爽人人干人人爱| 欧洲精品卡2卡3卡4卡5卡区| 免费电影在线观看免费观看| 看免费av毛片| 九九热线精品视视频播放| 国产亚洲精品综合一区在线观看| 欧美日韩乱码在线| 欧美精品啪啪一区二区三区| 母亲3免费完整高清在线观看| 亚洲欧美精品综合一区二区三区| 俺也久久电影网| 在线观看免费视频日本深夜| 18禁裸乳无遮挡免费网站照片| 午夜两性在线视频| 麻豆国产av国片精品| 我的老师免费观看完整版| 两性夫妻黄色片| 黄色 视频免费看| 99在线视频只有这里精品首页| 我的老师免费观看完整版| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 午夜精品久久久久久毛片777| 熟女电影av网| 给我免费播放毛片高清在线观看| 亚洲熟女毛片儿| 国产精品 国内视频| 日本熟妇午夜| 亚洲国产精品成人综合色| 国产精品国产高清国产av| 欧美在线一区亚洲| 日韩精品中文字幕看吧| 国产精品久久视频播放| 国内精品美女久久久久久| 精品一区二区三区视频在线 | 久久精品人妻少妇| 在线看三级毛片| 人人妻人人澡欧美一区二区| 搞女人的毛片| 欧美日韩乱码在线| 岛国在线观看网站| 午夜福利在线观看吧| av福利片在线观看| 俺也久久电影网| 老汉色∧v一级毛片| 99久久国产精品久久久| 丁香欧美五月| 国产 一区 欧美 日韩| 99在线人妻在线中文字幕| 免费一级毛片在线播放高清视频| 免费看美女性在线毛片视频| 国产黄片美女视频| 久久久色成人| 欧美又色又爽又黄视频| 亚洲av电影不卡..在线观看| 精品乱码久久久久久99久播| 人人妻人人看人人澡| 国产伦人伦偷精品视频| 麻豆成人av在线观看| 禁无遮挡网站| 久久久久亚洲av毛片大全| 久久草成人影院| 国产蜜桃级精品一区二区三区| 91老司机精品| 天天添夜夜摸| 黄色日韩在线| tocl精华| 欧美日本亚洲视频在线播放| 亚洲国产精品久久男人天堂| 精品乱码久久久久久99久播| 亚洲国产欧美网| 亚洲av第一区精品v没综合| 三级男女做爰猛烈吃奶摸视频| 欧美乱妇无乱码| 999久久久国产精品视频| 国产成人欧美在线观看| 18禁国产床啪视频网站| 看黄色毛片网站| а√天堂www在线а√下载| 日本成人三级电影网站| 久久久色成人| 在线视频色国产色| 午夜久久久久精精品| 亚洲精品粉嫩美女一区| 日韩国内少妇激情av| 性色avwww在线观看| 国内少妇人妻偷人精品xxx网站 | 香蕉av资源在线| 国产伦精品一区二区三区四那| 久久中文字幕一级| 小说图片视频综合网站| 狂野欧美白嫩少妇大欣赏| 少妇人妻一区二区三区视频| 伦理电影免费视频| 欧美日韩中文字幕国产精品一区二区三区| 少妇的逼水好多| 国产成人精品久久二区二区91| 亚洲成人免费电影在线观看| 在线看三级毛片| 国产主播在线观看一区二区| 成人一区二区视频在线观看| 搡老妇女老女人老熟妇| 一区二区三区激情视频| 桃色一区二区三区在线观看| 国产精品久久久久久精品电影| 国产欧美日韩精品一区二区| 国产高清有码在线观看视频| 不卡一级毛片| 少妇的丰满在线观看| 色吧在线观看| 亚洲午夜理论影院| 在线观看午夜福利视频| 国产三级中文精品| 亚洲国产精品久久男人天堂| 亚洲国产欧美人成| 亚洲色图av天堂| 午夜精品久久久久久毛片777| 少妇人妻一区二区三区视频| 日本一本二区三区精品| 少妇裸体淫交视频免费看高清| 国产视频一区二区在线看| 99久久国产精品久久久| 国产主播在线观看一区二区| 真人做人爱边吃奶动态| 免费搜索国产男女视频| ponron亚洲| 成人三级做爰电影| 国产黄a三级三级三级人| 日本熟妇午夜| 亚洲欧美日韩高清专用| www日本黄色视频网| 高潮久久久久久久久久久不卡| 嫩草影院精品99| 网址你懂的国产日韩在线| 在线观看日韩欧美| 五月伊人婷婷丁香| 成人一区二区视频在线观看| 99国产精品一区二区蜜桃av| 麻豆久久精品国产亚洲av| 成人鲁丝片一二三区免费| 免费电影在线观看免费观看| 亚洲无线观看免费| 毛片女人毛片| 国产精品久久久av美女十八| 日本一本二区三区精品| 国产伦人伦偷精品视频| 两个人的视频大全免费| 成熟少妇高潮喷水视频| 国产av不卡久久| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜亚洲精品久久| 美女免费视频网站| 99热精品在线国产| 亚洲黑人精品在线| 嫩草影视91久久| 亚洲精品乱码久久久v下载方式 | 99热6这里只有精品| 男女做爰动态图高潮gif福利片| 国产1区2区3区精品| 51午夜福利影视在线观看| 亚洲电影在线观看av| 国产单亲对白刺激| 亚洲av成人av| 国产精品国产高清国产av| 日韩有码中文字幕| 亚洲美女视频黄频| 国产激情偷乱视频一区二区| 老鸭窝网址在线观看| 91字幕亚洲| www.熟女人妻精品国产| 叶爱在线成人免费视频播放| 欧美三级亚洲精品| 国产v大片淫在线免费观看| 天堂√8在线中文| 免费无遮挡裸体视频| 无人区码免费观看不卡| 国产欧美日韩一区二区三| 一本综合久久免费| 亚洲欧美日韩高清专用| 在线观看一区二区三区| 在线国产一区二区在线| 俺也久久电影网| 精品久久久久久,| 中文字幕久久专区| 嫩草影院入口| 女人高潮潮喷娇喘18禁视频| 亚洲av成人精品一区久久| 成人午夜高清在线视频| 大型黄色视频在线免费观看| av黄色大香蕉| 在线观看美女被高潮喷水网站 | 欧美极品一区二区三区四区| 免费av毛片视频| 人妻丰满熟妇av一区二区三区| 亚洲av成人不卡在线观看播放网| 色综合站精品国产| 香蕉久久夜色| 日日干狠狠操夜夜爽| 亚洲精品在线美女| 熟女人妻精品中文字幕| 91字幕亚洲| aaaaa片日本免费| 2021天堂中文幕一二区在线观| 欧美色欧美亚洲另类二区| 色噜噜av男人的天堂激情| 一本精品99久久精品77| 国产亚洲欧美在线一区二区| 免费在线观看视频国产中文字幕亚洲| 成人性生交大片免费视频hd| 99久久精品热视频| 搡老熟女国产l中国老女人| 国产欧美日韩一区二区精品| 草草在线视频免费看| 一本综合久久免费| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区三| 国产亚洲精品一区二区www| 精品一区二区三区视频在线 | 亚洲熟妇熟女久久| 变态另类成人亚洲欧美熟女| 国产97色在线日韩免费|