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

    橢圓歧型三通沖蝕磨損數(shù)值模擬

    2022-09-27 12:43:30莫麗馮滿陳行郭振興
    表面技術(shù) 2022年9期
    關(guān)鍵詞:管匯沖蝕三通

    莫麗,馮滿,陳行,郭振興

    橢圓歧型三通沖蝕磨損數(shù)值模擬

    莫麗1,2,馮滿1,陳行1,郭振興1

    (1.西南石油大學(xué) 機(jī)電工程學(xué)院,成都 610500;2.石油天然氣裝備技術(shù)四川省科技資源共享服務(wù)平臺,成都 610500)

    高壓管匯;沖蝕磨損;橢圓截面;歧型三通管;數(shù)值模擬

    近年來,隨著全球?qū)δ茉葱枨蟮脑黾?,各類非常?guī)油氣資源的開發(fā)越來越受到重視[1]。國家能源結(jié)構(gòu)調(diào)整改革的實(shí)施,使得以頁巖氣為代表的非常規(guī)天然氣的開發(fā)逐步推進(jìn),其年產(chǎn)量已占據(jù)中國石化天然氣總產(chǎn)量的25%。頁巖氣開采過程中,由于低滲致密儲層而導(dǎo)致的采氣效率下降的問題,可采用分段壓裂技術(shù)對滲流環(huán)境和井底流動條件進(jìn)行改善,提升油氣井的產(chǎn)能[2]。在加砂壓裂過程中,高壓管匯作為壓裂系統(tǒng)的重要組成部分,其工況較為惡劣:壓裂液和所含化學(xué)物質(zhì)引發(fā)的化學(xué)腐蝕和電化學(xué)腐蝕,固相顆粒沖擊導(dǎo)致的管道沖蝕磨損,超高壓強(qiáng)下的交變載荷等。以上因素導(dǎo)致高壓管匯極易發(fā)生失效,在壓裂作業(yè)時(shí)出現(xiàn)壓裂液泄露、管匯爆裂等諸多危險(xiǎn)情況,對井場工作人員的人生安全和壓裂作業(yè)的順利進(jìn)行構(gòu)成威脅。圖1為2008年到2015年高壓管匯事故統(tǒng)計(jì),期間高壓管匯爆裂事故高達(dá)36起,其中,活動彎頭爆裂事故和三通管爆裂事故共計(jì)24起,占比67%[3]。而引發(fā)此類事故的主要原因是:固相顆粒在流動過程中對管壁的沖蝕磨損,導(dǎo)致管壁變薄,甚至穿孔,最終導(dǎo)致管道無法承受工作壓力,引發(fā)爆裂事故。因此,探究固相顆粒對管道的沖蝕磨損機(jī)理,進(jìn)而提出能降低管道的沖蝕磨損率的可行性方案,有較高的研究價(jià)值。

    季楚凌等[4]通過耐磨生物體表的表面形態(tài),建立了加厚彎管物理模型,并運(yùn)用沖蝕理論進(jìn)行數(shù)值模擬,分析了顆粒運(yùn)動參數(shù)與耐磨性的關(guān)系。張繼信等[5-6]利用CFD軟件建立了水力壓裂的固液兩相流模型,通過數(shù)值模擬得到了彎管沖蝕磨損最嚴(yán)重區(qū)域,并根據(jù)失效樣品分析了管道沖蝕的宏觀、微觀形貌特征。姜心等[7]討論了在不同沖蝕角度下,40CrNi2Mo合金鋼材料的沖蝕磨損機(jī)制,并采用電鏡分析了材料沖蝕磨損后的形貌特征。金雪梅等[3]對加砂壓裂過程中高壓管匯的失效原因進(jìn)行分析,提出了相關(guān)建議。成芳等[8]基于液–固兩相流理論和Fluent沖蝕模型,建立了彎管沖蝕磨損模型,并引入斯托克斯數(shù)探究了彎管沖蝕磨損區(qū)域的變化情況。劉洪斌等[9]基于數(shù)值模擬計(jì)算研究了三通管沖蝕發(fā)生的機(jī)理及位置,探討了高壓管匯的沖蝕特點(diǎn)。彭方現(xiàn)等[10-11]采用標(biāo)準(zhǔn)–epsilon模型和SMPLE算法,對T型三通管道內(nèi)部流場進(jìn)行了仿真分析,并對氣固兩相流作用下的T型管和歧型管的沖蝕泄漏問題進(jìn)行了模擬分析,為減緩管道沖蝕磨損和快速定位、封堵破損管道等提供參考依據(jù)。胡開松等[12]基于液固兩相流和沖蝕磨損理論,研究了魚尾管不同夾角下,顆粒直徑、顆粒形狀、質(zhì)量流量、流體速度等4種因素對魚尾管沖蝕的影響,總結(jié)了高壓魚尾管匯內(nèi)壁的沖蝕磨損規(guī)律。祝效華等[13]研究了含砂壓裂液對高壓管匯三通沖蝕磨損的影響規(guī)律及主要影響因素,并預(yù)測了三通管在使用過程中易發(fā)生沖蝕磨損的部位。Zhang等[14]基于CFD系統(tǒng)地研究了壓力管匯中常見的T型管和歧型管,探討了壓力管匯不同部位和工藝過程環(huán)境中的流體場分布、沖蝕機(jī)理和沖蝕磨損的規(guī)律。

    圖1 2008—2015年高壓管匯事故統(tǒng)計(jì)

    眾多研究者對三通管和彎管的沖蝕理論以及沖蝕特性做了大量研究,并總結(jié)出了普遍性的規(guī)律,極大推動了管道沖蝕方面的研究進(jìn)展。但這些研究對于解決降低管道沖蝕率的問題,所提出的方案都是基于改變流體的沖蝕參數(shù)和改變管道角度(歧型管夾角和彎管角度)的,尚沒有做關(guān)于改變管道截面形狀的研究,對管道沖蝕規(guī)律的研究仍有所欠缺。基于此,本文通過歧型三通管沖蝕機(jī)理的研究,提出一種橢圓截面的歧型三通管。本文基于固液兩相流建立了沖蝕模型,系統(tǒng)地研究了橢圓歧型三通管的橢圓形狀變化對最大沖蝕率的影響,并研究比較不同工況下,橢圓歧型三通管和普通歧型三通管最大沖蝕率的情況。通過分析橢圓歧型三通管的沖蝕規(guī)律,為橢圓歧型三通管的設(shè)計(jì)提供理論依據(jù)和新的方向。

    1 歧型三通管結(jié)構(gòu)模型

    為改善滲流環(huán)境和井底流動條件,提升油氣井的產(chǎn)能,需要對油氣儲層進(jìn)行壓裂作業(yè)。壓裂時(shí)需要提供極強(qiáng)的注入壓力和較大的流量,因此,壓裂作業(yè)通常由多臺壓裂車并聯(lián),將壓裂液泵出至高壓管匯匯集,再經(jīng)由高壓管線接入到井口。高壓管匯主要包含Y型三通管、歧型三通管和U型彎管,其結(jié)構(gòu)如圖2所示。圖2a為高壓管匯實(shí)物。圖2b為高壓管匯模型,其中接口1~6為壓裂車接口,接口7為高壓管線接口[9]。

    歧型三通管是整個高壓管匯的關(guān)鍵部件,有著匯集支管來流、引導(dǎo)流體流向的作用,其承受了壓裂液和所含化學(xué)物質(zhì)引發(fā)的化學(xué)腐蝕和電化學(xué)腐蝕、固相顆粒沖擊導(dǎo)致的管道沖蝕磨損以及超高壓強(qiáng)下的交變載荷等惡劣工況,極易發(fā)生失效。因此,本文旨在不改變歧型三通管支管流量的情況下,通過改變截面形狀以降低其沖蝕率。本研究所采用的歧型三通管結(jié)構(gòu)尺寸模型如圖3所示。

    圖2 高壓管匯結(jié)構(gòu)

    圖3 歧型三通管及實(shí)驗(yàn)?zāi)P统叽缃Y(jié)構(gòu)

    2 沖蝕磨損理論

    2.1 控制方程

    1)連續(xù)性方程

    流體流動滿足質(zhì)量守恒定律,有:

    2)動量守恒方程

    流體微元滿足牛頓第二定律[15],有:

    3)能量守恒方程

    2.3 沖蝕磨損計(jì)算模型

    根據(jù)管道內(nèi)固相顆粒的碰撞規(guī)律[18],與管壁的沖蝕磨損模型如下:

    2.4 壁面碰撞恢復(fù)方程

    根據(jù)本文所研究的對象,選擇的恢復(fù)系數(shù)方程[19]如下:

    2.5 沖擊角函數(shù)

    對沖擊角函數(shù)采用分段線性方式定義,定義數(shù)據(jù)見表1。

    表1 沖擊角函數(shù)定義

    Tab.1 Definition of impact angle function

    3 數(shù)值仿真計(jì)算

    3.1 CFD模型

    本文提出了一種橢圓截面三通管,模型如圖4a所示,橢圓管末端經(jīng)過渡管段可與正常圓形管連接。為便于分析,本文僅選取橢圓管段與主管段進(jìn)行仿真分析。分析模型如圖4b和圖4c所示,改變支管截面–的形狀,利用fluent進(jìn)行數(shù)值模擬仿真。為使橢圓歧型三通管與普通圓截面歧型三通管的最大沖蝕率形成直觀的對比,需保證支管流速一致的情況下,流量相同,故截面–的面積基本不隨截面形狀的變化而變化。截面形狀變化如圖4c所示。定義截面橢圓長半軸為,短半軸為,=/。越小,橢圓度越大,=1時(shí)為普通圓截面歧型三通管。為研究橢圓支管對管道沖蝕率的影響,本文對各橢圓截面具體參數(shù)設(shè)定如表2所示。

    3.2 數(shù)值模擬參數(shù)設(shè)定

    本文以壓裂現(xiàn)場工況條件為參考,設(shè)定仿真邊界條件如表3所示。

    支管(Inlet 1)主要與壓裂車相連,流動環(huán)境變化不大,根據(jù)壓裂現(xiàn)場壓裂液指標(biāo),設(shè)定流速為10.19 m/s,黏度、質(zhì)量流量、粒子直徑與主管(Inlet 2)設(shè)定保持一致。主管相關(guān)參數(shù)根據(jù)不同壓裂工況,參數(shù)設(shè)定如表4所示。

    圖4 歧型三通管橢圓截面與i值

    表2 截面參數(shù)

    Tab.2 Section parameters

    表3 邊界條件

    Tab.3 Boundary condition

    表4 仿真變量

    Tab.4 Simulation variables

    3.3 網(wǎng)格劃分及網(wǎng)格無關(guān)性驗(yàn)證

    對數(shù)值計(jì)算的結(jié)構(gòu)模型采用ICEM進(jìn)行六面體網(wǎng)格劃分,并設(shè)定邊界層,網(wǎng)格質(zhì)量良好,網(wǎng)格劃分如圖5所示。

    網(wǎng)格數(shù)量過少會影響計(jì)算結(jié)果的準(zhǔn)確性和穩(wěn)定性,網(wǎng)格數(shù)量過多并不會對計(jì)算結(jié)果有正面影響,還會增加計(jì)算中心的工作量,過多占用資源,造成不必要的浪費(fèi)[20]。本文計(jì)算模型數(shù)量較多,耗時(shí)較長,因此,對網(wǎng)格數(shù)量有必要進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,以選取合適的網(wǎng)格數(shù)量。網(wǎng)格無關(guān)性驗(yàn)證結(jié)果如圖6所示,當(dāng)網(wǎng)格數(shù)量達(dá)到38.5萬及以上時(shí),對計(jì)算結(jié)果基本沒有影響,且計(jì)算結(jié)果穩(wěn)定。計(jì)算迭代步長大于800步時(shí),已經(jīng)收斂,計(jì)算結(jié)果穩(wěn)定。故模型網(wǎng)格劃分?jǐn)?shù)量應(yīng)選擇40萬左右,迭代步數(shù)設(shè)定為800步。

    圖5 歧型三通管網(wǎng)格劃分

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

    4 數(shù)值仿真結(jié)果

    4.1 壓力云圖與速度云圖

    歧型三通管數(shù)值模擬結(jié)果的壓力云圖和速度云圖如圖7所示。從流場壓力云圖可以看出,主管端壓強(qiáng)為85 MPa,支管端壓強(qiáng)為85.2 MPa,主管與支管的交匯處(靠近出口端的相貫線)附近,壓力為84.9 MPa,出口(Outlet)壓強(qiáng)分層現(xiàn)象明顯。從流場速度云圖可以看出,支管與主管相交匯時(shí),流體流速增大,流速最大處位于靠近出口端的相貫線處,流速最大區(qū)域之后出現(xiàn)明顯的相對低速區(qū)域,并在出口端產(chǎn)生流速分層現(xiàn)象。

    由于支管內(nèi)流體流向主管時(shí),具有較大的慣性,對主管的流體有向遠(yuǎn)離支管端擠壓的作用,故在靠近出口端的相貫線附近形成了明顯的低壓區(qū);當(dāng)主管與支管流體交匯后,由于流體慣性和管道幾何結(jié)構(gòu),將會在靠近出口端的相貫線附近形成回旋的渦流,持續(xù)沖擊管道,故形成流體相對低速區(qū)。

    圖7 壓力云圖和速度云圖

    4.2 夾角與i值變化對最大沖蝕率的影響

    歧型三通管不同夾角下不同橢圓截面(值)的最大沖蝕率如圖8所示。相同橢圓截面,夾角為60°時(shí),管道最大沖蝕率最大;夾角為50°和70°時(shí),管道最大沖蝕率比較接近;夾角為40°時(shí),管道最大沖蝕率最小;最大沖蝕率比最小沖蝕率增大了2.18倍。試驗(yàn)選取的所有夾角,均滿足:當(dāng)由0.76增大至1過程中,管道最大沖蝕率先緩慢減小,后迅速增大,最小值在=0.85附近;同一夾角不同橢圓截面時(shí),最大沖蝕率比最小沖蝕率增大了3.66倍。

    其沖蝕機(jī)理是:管道內(nèi)流體在交匯后速度增加,并在主管與支管后相貫線附近達(dá)到最大值,且主管與支管存在一定角度,流體在交匯后,固相顆粒間相互碰撞產(chǎn)生大量不規(guī)則運(yùn)動,并高速撞擊管道內(nèi)壁。改變主管與支管間夾角,當(dāng)夾角由40°增加至60°過程中,當(dāng)角度較小時(shí),支管內(nèi)流體能較為順利匯入主管中,當(dāng)角度逐漸增大時(shí),支管與主管內(nèi)的流體碰撞程度加重,引起相貫線及其附近的沖蝕加劇;當(dāng)夾角由60°增加至80°過程中,支管與主管內(nèi)流體交匯時(shí)碰撞程度進(jìn)一步增加,并引發(fā)較強(qiáng)湍流,同時(shí)支管內(nèi)流體由于慣性流向主管底部,因而相貫線附近的沖蝕情況得到緩解。隨著支管截面橢圓度逐漸增加,支管與主管間的相貫線曲率逐漸減小,使主管與支管的流體交匯區(qū)域逐漸轉(zhuǎn)移到主管上部,且交匯區(qū)域在流速方向拉長,固相顆粒間的碰撞劇烈程度有所下降,對管道相貫線附近的撞擊減小。但在相貫線后部區(qū)域,會形成較強(qiáng)渦流,部分固相顆粒在此處向上運(yùn)動,對主管道上臂形成撞擊。

    圖8 夾角與i值對沖蝕率的影響

    為便于比較,后文統(tǒng)一選取在沖蝕率最大的60°夾角時(shí),研究不同橢圓截面的歧型三通管的沖蝕磨損情況。

    4.3 60°夾角時(shí)不同i值的沖蝕率

    夾角為60°時(shí),不同橢圓截面的歧型三通管沖蝕結(jié)果如圖9所示。由數(shù)值仿真結(jié)果可知,從1減小至0.85過程中,最大沖蝕率逐漸減??;從0.85減小至0.76過程中,最大沖蝕率逐漸增大,沖蝕率最小值出現(xiàn)在=0.85時(shí)。從1減小至0.76過程中,相貫線處沖蝕率逐漸減小,主管中部的沖蝕率逐漸增加,且沖蝕區(qū)域逐漸由主管中部轉(zhuǎn)移到主管上部,沖蝕最嚴(yán)重的區(qū)域由相貫線處轉(zhuǎn)變?yōu)橄嘭灳€附近的回旋渦流區(qū)(低流速區(qū)),但渦流區(qū)最大沖蝕率仍低于相貫線處。

    4.4 i為0.85和1的沖蝕率對比

    圖10為歧型三通管在=0.85和=1時(shí),不同固相顆粒質(zhì)量流率和主管流量對管道最大沖蝕率的影響。由圖10可知,=0.85和=1時(shí),管道最大沖蝕率均隨主管流量增大而增大,隨固相顆粒質(zhì)量流率減小而減??;=1時(shí),在主管不同流量下最大沖蝕率均明顯大于=0.85時(shí)的,最高倍數(shù)達(dá)1.9倍;=1時(shí),在不同質(zhì)量流率下最大沖蝕率均明顯大于=0.85時(shí)的,最高倍數(shù)達(dá)1.93倍,且=1時(shí),最大沖蝕率隨質(zhì)量流率變化的速度遠(yuǎn)高于=0.85。

    圖9 沖蝕區(qū)域隨i值的變化情況

    圖10 流量、質(zhì)量流率對管道最大沖蝕率的影響

    圖11為歧型三通管在=0.85和=1時(shí),不同固相顆粒直徑和流體黏度對管道最大沖蝕率的影響。由圖11可知,=0.85和=1時(shí),管道最大沖蝕率均隨固相顆粒直徑增大而減小,隨流體黏度增大而先增大后減小,在流體黏度為15 mpa·s時(shí)達(dá)到最大值;=1時(shí),在不同固相顆粒直徑下最大沖蝕率均明顯大于=0.85時(shí)的,最高倍數(shù)達(dá)2.74倍;=1時(shí),在不同流體黏度下最大沖蝕率均明顯大于=0.85時(shí)的,最高倍數(shù)達(dá)2.36倍。

    圖11 粒子直徑、流體黏度對管道最大沖蝕率的影響

    5 結(jié)論

    1)橢圓歧型三通管不同值時(shí)均在夾角為60°時(shí)最大沖蝕率最大;不同夾角時(shí),值由0.76增大到1過程中,最大沖蝕率均先減小后增大,且在=0.85附近有最小值,即=0.85時(shí)降低沖蝕的表現(xiàn)最為優(yōu)異。

    2)橢圓歧型三通管值由1減小到0.76過程中,相貫線處的沖蝕率逐漸減小,主管中部的沖蝕率逐漸增加,且沖蝕區(qū)域逐漸由主管中部轉(zhuǎn)移到主管上部;沖蝕最嚴(yán)重的區(qū)域由相貫線處轉(zhuǎn)變到回旋渦流區(qū)(低流速區(qū)),但渦流區(qū)的最大沖蝕率仍低于相貫線處。

    3)橢圓歧型三通管=0.85時(shí),最大沖蝕率在不同工況下的變化趨勢與=1時(shí)基本一致,滿足普通圓截面歧型三通管的沖蝕磨損規(guī)律。

    4)在不同工況下,橢圓歧型三通管=0.85相較于=1,降低沖蝕率的效果十分顯著。

    [1] 步玉環(huán), 常智楊, 邵子璇, 等. 適用于稠油熱采井熱膨脹水泥的外摻料優(yōu)選[J]. 鉆井液與完井液, 2016, 33(4): 87-91.

    BU Yu-huan, CHANG Zhi-yang, SHAO Zi-xuan, et al. Optimization of Admixture for Thermal Expansion Ce-ment Used for Cementing Heavy Oil Thermal Reco-very Well[J]. Drilling Fluid & Completion Fluid, 2016, 33(4): 87-91.

    [2] 徐玉龍, 程智遠(yuǎn), 劉志斌, 等. 水平井分段壓裂投球滑套承壓接觸分析[J]. 石油機(jī)械, 2012, 40(8): 93-96.

    XU Yu-long, CHENG Zhi-yuan, LIU Zhi-bin, et al. Pres-sure Bearing and Contact Analysis of Ball Injection Sli-ding Sleeve in Horizontal Well Staged Fracturing[J]. China Petroleum Machinery, 2012, 40(8): 93-96.

    [3] 金雪梅, 張祥來, 廖浩, 等. 加砂壓裂過程中高壓管匯失效爆裂分析[J]. 安全, 2017, 38(1): 17-18.

    JIN Xue-mei, ZHANG Xiang-lai, LIAO Hao, et al. Failure and Burst Analysis of High Pressure Manifold during Sand Fracturing[J]. Safety, 2017, 38(1): 17-18.

    [4] 季楚凌, 李長俊, 馬樹鋒, 等. 彎管仿生耐磨方法數(shù)值模擬[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2016, 47(10): 3582-3589.

    JI Chu-ling, LI Chang-jun, MA Shu-feng, et al. Numerical Simulation of Bionic Method for Improving Property of Elbow Erosion Resistance[J]. Journal of Central South Uni-versity (Science and Technology), 2016, 47(10): 3582-3589.

    [5] 張繼信, 樊建春, 汪彤, 等. 壓裂液對高壓彎管沖蝕作用的數(shù)值分析[J]. 潤滑與密封, 2013, 38(4): 27-31.

    ZHANG Ji-xin, FAN Jian-chun, WANG Tong, et al. Re-search on Erosion Wear of Fracturing Fluid on High Pre-ssure Elbow[J]. Lubrication Engineering, 2013, 38(4): 27-31.

    [6] 張繼信, 康健, 樊建春, 等. 高壓彎管沖蝕失效分析及數(shù)值模擬[J]. 潤滑與密封, 2016, 41(10): 40-46.

    ZHANG Ji-xin, KANG Jian, FAN Jian-chun, et al. Fail-ure Analysis and Numerical Simulation of High Pressure Elbow Erosion Wear[J]. Lubrication Engineering, 2016, 41(10): 40-46.

    [7] 姜心, 張來斌, 樊建春, 等. 沖蝕角度對40CrNi2Mo材料沖蝕磨損性能的影響[J]. 潤滑與密封, 2012, 37(6): 24-26, 45.

    JIANG Xin, ZHANG Lai-bin, FAN Jian-chun, et al. Eff-ect of Impact Angle on Erosive Resistance of 40CrNi2Mo[J]. Lubrication Engineering, 2012, 37(6): 24-26, 45.

    [8] 成芳, 易先中, 彭灼, 等. 頁巖氣壓裂彎管中液固兩相流沖蝕磨損的數(shù)值模擬[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2019, 15(11): 29-35.

    CHENG Fang, YI Xian-zhong, PENG Zhuo, et al. Num-erical Simulation on Erosion Wear in Shale Gas Frac-turing Bend Pipe with Liquid-Solid Two-Phase Flow[J]. Journal of Safety Science and Technology, 2019, 15(11): 29-35.

    [9] 劉洪斌, 牟浩. 加砂壓裂中固體支撐劑對壓裂管道的沖蝕磨損研究[J]. 中國安全生產(chǎn)科學(xué)技術(shù), 2018, 14(1): 87-94.

    LIU Hong-bin, MU Hao. Study on Erosion Wear of Frac-turing Pipe Caused by Solid Proppant in Sand Frac-tu-ring[J]. Journal of Safety Science and Technology, 2018, 14(1): 87-94.

    [10] 彭方現(xiàn). 三通管道沖蝕泄漏流場數(shù)值模擬分析研究[D]. 太原: 中北大學(xué), 2020.

    PENG Fang-xian. Numerical Simulation Analysis of Ero-sion Leakage Flow Field of Tee Pipeline[D]. Taiyuan: North University of China, 2020.

    [11] 彭方現(xiàn), 閆宏偉, 李亞杰, 等. T型管道的沖蝕磨損數(shù)值模擬分析[J]. 當(dāng)代化工, 2020, 49(3): 733-736, 752.

    PENG Fang-xian, YAN Hong-wei, LI Ya-jie, et al. Nume-r-ical Simulation Analysis of Erosion Wear of T-Junction[J]. Contemporary Chemical Industry, 2020, 49(3): 733-736, 752.

    [12] 胡開松, 周思柱, 華劍, 等. 高壓魚尾管匯沖蝕規(guī)律研究[J]. 中國科技論文, 2020, 15(10): 1177-1181.

    HU Kai-song, ZHOU Si-zhu, HUA Jian, et al. Analysis on Erosion Laws of High-Pressure Fishtail Manifold[J]. China Sciencepaper, 2020, 15(10): 1177-1181.

    [13] 祝效華, 張覃, 張洋銘, 等. 高壓管匯三通沖蝕磨損特性研究[J]. 表面技術(shù), 2021, 50(7): 258-265.

    ZHU Xiao-hua, ZHANG Qin, ZHANG Yang-ming, et al. Study on Erosion Wear Characteristics of High Pressure Manifold Tee[J]. Surface Technology,2021, 50(7): 258- 265.

    [14] ZHANG J, ZHANG Hao, LIU Yang, et al. Erosion Wear Characteristics of Tee Tubes with Gas-Solid Two-Phase Flow[J]. Journal of Pressure Vessel Technology, 2021, 143(6): 256.

    [15] ZHANG Yong-li, MCLAURY B S, SHIRAZI S A. Improvements of Particle Near-Wall Velocity and Erosion Predictions Using a Commercial CFD Code[J]. Journal of Fluids Engineering, 2009, 131(3): 456.

    [16] 邱華, 牟俊玲, 徐燃霞, 等. 基于ANSYS CFX的網(wǎng)絡(luò)噴嘴內(nèi)部流場計(jì)算[J]. 紡織學(xué)報(bào), 2009, 30(2): 104-108.

    QIU Hua, MOU Jun-ling, XU Ran-xia, et al. Computation of Flow Patterns Inside Interlacer Jet Based on ANSYS CFX[J]. Journal of Textile Research, 2009, 30(2): 104- 108.

    [17] 姬長春, 張開源, 王玉棟, 等. 熔噴三維氣流場的數(shù)值計(jì)算與分析[J]. 紡織學(xué)報(bào), 2019, 40(8): 175-180.

    JI Chang-chun, ZHANG Kai-yuan, WANG Yu-dong, et al.Numerical Calculation and Analysis of Three-Dimensional Flow Field in Melt-Blown Process[J]. Journal of Textile Research, 2019, 40(8): 175-180.

    [18] 毛靖儒, 柳成文, 相曉偉. 彎管內(nèi)二次流對固粒磨損壁面的影響[J]. 西安交通大學(xué)學(xué)報(bào), 2004, 38(7): 746-749.

    MAO Jing-ru, LIU Cheng-wen, XIANG Xiao-wei. Effect of Secondary Flow on Erosion from Solid Particles in 90° Curved Duct of Quadrate Section[J]. Journal of Xi'an Jiaotong University, 2004, 38(7): 746-749.

    [19] FORDER A, THEW M, HARRISON D. A Numerical Investigation of Solid Particle Erosion Experienced within Oilfield Control Valves[J]. Wear, 1998, 216(2): 184-193.

    [20] 楊周瑾. 集砂桶液位升高對旋流式除砂器氣相分離效果影響的研究[D]. 西安: 西安石油大學(xué), 2019.

    YANG Zhou-jin. Study on the Separation Efficiency of Cyclone Separator while Sand Level Rising in Sand Col-lecting Barrel[D]. Xi'an: Xi'an Shiyou University, 2019.

    Numerical Simulation of Erosion Wear of Oval Disproportionate Tee

    1,2,1,1,1

    (1. School of Mechanical and Electrical Engineering, Southwest Petroleum University, Chengdu 610500, China; 2. Sichuan Province Science and Technology Resource Sharing Service Platform for Petroleum and Natural Gas Equipment Technology, Chengdu 610500, China)

    To solve the problem of high failure rate of high-pressure tee manifolds in the petrochemical industry due to erosion and wear, reduce uncontrollable costs in petrochemical production, and increase the service life of disproportionate tee manifolds and the safety factor of the manifold area. Based on the erosion theory of solid-liquid two-phase flow, a numerical simulation model of erosion and wear of the elliptical cross-section disproportionate tee is established to predict the geometric location of erosion and wear of the elliptical cross-section disproportionate tee during operation. Study the influence of the elliptical shape and azimuth angle of the elliptical cross-section disproportionate tee on the maximum erosion rate of the pipeline. Compare the elliptical cross-section when=0.85 and the ordinary circular cross-section when=1. The two manifold three-way pipes are different influence of mass flow rate, particle diameter, liquid flow rate, and fluid viscosity on the maximum erosion rate of the pipeline. The experimental results found that when the included angle is 60°, the maximum erosion rate of the pipeline is the largest. When the included angle is 50° and 70°, the maximum erosion rate of the pipeline is relatively close. When the included angle is 40°, the maximum erosion rate of the pipeline is the smallest; The included angles are all satisfied: whenincreases from 0.76 to 1, the maximum erosion rate of the pipeline decreases slowly at first and then increases rapidly. There is a minimum near=0.85, and the maximum is 3.66 times larger than the minimum; In the process of=1 decreasing to=0.76, the erosion rate at the intersecting line gradually decreases, the erosion rate in the middle of the main pipe gradually increases, and the erosion area gradually shifts from the middle of the main pipe to the upper part of the main pipe, with the most severe erosion The area changes from the intersecting line to the swirling vortex area near the intersecting line (low velocity area); when=1, under different solid particle mass flow rates and main flow rates, the maximum erosion rate is significantly greater than=0.85, The increase times were 1.9 times and 1.93 times respectively; when=1, under different solid particle diameters and fluid viscosities, the maximum erosion rate was significantly greater than=0.85, and the increase times were 2.74 times and 2.36 times, respectively. The following conclusions are drawn when thevalue of the elliptical cross-section disproportionate tee increases from 0.76 to 1, the maximum erosion rate first decreases and then increases, and there is a minimum near=0.85; when=0.85, the maximum erosion The changing trend of the erosion rate under different working conditions is basically the same as that when=1, but the reduction effect of the erosion rate is very significant.

    high pressure manifolds;erosion wear;elliptic section;disproportionate tee;numerical simulation

    2021-07-12;Revised:2021-12-30

    MO Li (1968-), Female, Master, Professor, Research focus: fluid machinery and oil and gas equipment.

    莫麗,馮滿,陳行,等. 橢圓歧型三通沖蝕磨損數(shù)值模擬[J]. 表面技術(shù), 2022, 51(9): 151-159.

    TH117

    A

    1001-3660(2022)09-0151-09

    10.16490/j.cnki.issn.1001-3660.2022.09.000

    2021–07–12;

    2021–12–30

    四川省省院省??萍己献餮邪l(fā)項(xiàng)目(21SYSX0054)

    Fund:Sichuan Provincial Academy Provincial School Science and Technology Cooperation Research and Development Project (21SYSX0054)

    莫麗(1968—),女,碩士,教授,主要研究方向?yàn)榱黧w機(jī)械與油氣裝備。

    MO Li, FENG Man, CHEN Hang, et al. Numerical Simulation of Erosion Wear of Oval Disproportionate Tee[J]. Surface Technology, 2022, 51(9): 151-159.

    責(zé)任編輯:萬長清

    猜你喜歡
    管匯沖蝕三通
    140MPa井口壓裂四通管道沖蝕分析
    基于模糊綜合評價(jià)的水下管匯結(jié)構(gòu)可靠性分析*
    隆德縣三通古碑考釋
    西夏研究(2017年4期)2017-08-17 08:53:10
    輸氣管道砂沖蝕的模擬實(shí)驗(yàn)
    拋錨撞擊水下管匯的數(shù)值模擬研究
    海洋工程(2016年4期)2016-10-12 03:21:26
    “三通兩平臺”推動下的我國在線教育出版發(fā)展
    新聞傳播(2016年4期)2016-07-18 10:59:22
    讀寫結(jié)合教學(xué)中“三合三通”
    環(huán)氧樹脂及其復(fù)合材料的固體顆粒沖蝕磨損
    番禺35-1/35-2氣田水下管匯嵌入式在線安裝方案設(shè)計(jì)及建造關(guān)鍵技術(shù)
    海洋石油管匯安裝方法及安裝過程力學(xué)分析
    欧美极品一区二区三区四区| 欧美97在线视频| 国产伦精品一区二区三区四那| 麻豆成人av视频| 秋霞伦理黄片| 最近中文字幕2019免费版| 午夜福利在线在线| 亚洲18禁久久av| 日本免费a在线| 一个人观看的视频www高清免费观看| 午夜爱爱视频在线播放| 国产高清三级在线| 欧美激情国产日韩精品一区| 亚洲精品亚洲一区二区| 亚洲熟妇中文字幕五十中出| 久久99热这里只频精品6学生| 精品欧美国产一区二区三| 国产一级毛片在线| 欧美变态另类bdsm刘玥| 91久久精品电影网| 亚洲av.av天堂| 亚洲av电影在线观看一区二区三区 | 丰满少妇做爰视频| 午夜福利成人在线免费观看| 国产亚洲精品av在线| 九九在线视频观看精品| 国产午夜福利久久久久久| 国产 一区精品| 九九久久精品国产亚洲av麻豆| 男女国产视频网站| av女优亚洲男人天堂| 精品久久久久久成人av| 看黄色毛片网站| 欧美日韩一区二区视频在线观看视频在线 | 婷婷色av中文字幕| 成人一区二区视频在线观看| 亚洲最大成人中文| 天堂中文最新版在线下载 | 国产午夜精品论理片| 国产老妇女一区| 亚洲欧美日韩无卡精品| 我要看日韩黄色一级片| 在线 av 中文字幕| 乱码一卡2卡4卡精品| 国产伦一二天堂av在线观看| 午夜福利高清视频| 久久久国产一区二区| 午夜福利高清视频| 日韩av免费高清视频| 日韩av在线大香蕉| 国产精品久久久久久久久免| 国产亚洲精品av在线| 永久网站在线| 久久久久精品久久久久真实原创| 久久久久久国产a免费观看| 午夜福利在线观看吧| 亚洲精品国产av成人精品| 插阴视频在线观看视频| 综合色丁香网| 国产精品女同一区二区软件| 2022亚洲国产成人精品| 91午夜精品亚洲一区二区三区| 男人舔奶头视频| 亚洲精品,欧美精品| 欧美日韩在线观看h| 久久国产乱子免费精品| eeuss影院久久| 亚洲精品成人久久久久久| av国产免费在线观看| 国产一区二区三区av在线| 免费观看a级毛片全部| 男人舔女人下体高潮全视频| 精品一区二区免费观看| 日韩精品有码人妻一区| 国产成人免费观看mmmm| 亚洲国产精品专区欧美| 亚洲av男天堂| 国产精品人妻久久久影院| 2018国产大陆天天弄谢| 床上黄色一级片| 听说在线观看完整版免费高清| 一二三四中文在线观看免费高清| 嫩草影院入口| 午夜福利视频1000在线观看| 欧美性猛交╳xxx乱大交人| 欧美zozozo另类| 嫩草影院入口| 久久久久网色| 国产精品国产三级国产专区5o| 内地一区二区视频在线| 国产一级毛片在线| 精品久久久久久久久久久久久| 国产激情偷乱视频一区二区| 十八禁网站网址无遮挡 | 男女下面进入的视频免费午夜| 国产精品人妻久久久影院| 搡老妇女老女人老熟妇| 亚洲精品日本国产第一区| 久久精品夜色国产| 九九爱精品视频在线观看| 久久这里有精品视频免费| 欧美不卡视频在线免费观看| 一级毛片aaaaaa免费看小| 久久久久久久久久久丰满| 91精品一卡2卡3卡4卡| 黄色日韩在线| 人体艺术视频欧美日本| 国产精品国产三级专区第一集| 国产成人免费观看mmmm| 精品一区在线观看国产| 日韩成人av中文字幕在线观看| 天天躁日日操中文字幕| 最近2019中文字幕mv第一页| 直男gayav资源| 精品久久久久久久人妻蜜臀av| 三级男女做爰猛烈吃奶摸视频| 街头女战士在线观看网站| 国产精品一区二区三区四区久久| 综合色av麻豆| 亚洲高清免费不卡视频| 狂野欧美白嫩少妇大欣赏| 内地一区二区视频在线| 少妇的逼水好多| 国产高清不卡午夜福利| 日韩av不卡免费在线播放| 国内揄拍国产精品人妻在线| 国产欧美另类精品又又久久亚洲欧美| 天美传媒精品一区二区| 国产亚洲av嫩草精品影院| 日本黄色片子视频| 男女那种视频在线观看| 免费看光身美女| 国产成人精品久久久久久| 亚洲国产精品sss在线观看| 午夜激情久久久久久久| 色播亚洲综合网| 国产美女午夜福利| 欧美高清性xxxxhd video| 中文字幕av在线有码专区| 在线a可以看的网站| 亚洲真实伦在线观看| 五月玫瑰六月丁香| 亚洲国产精品成人综合色| 午夜激情欧美在线| 男女那种视频在线观看| www.av在线官网国产| 男女啪啪激烈高潮av片| 99久久中文字幕三级久久日本| 搡老妇女老女人老熟妇| 熟女人妻精品中文字幕| 亚洲av在线观看美女高潮| 亚洲国产日韩欧美精品在线观看| 国产亚洲最大av| av卡一久久| 干丝袜人妻中文字幕| 日韩视频在线欧美| 国产白丝娇喘喷水9色精品| 神马国产精品三级电影在线观看| 热99在线观看视频| 嫩草影院精品99| 国产淫片久久久久久久久| 久久精品国产亚洲av天美| 日韩在线高清观看一区二区三区| 天堂中文最新版在线下载 | av在线蜜桃| www.色视频.com| 免费大片18禁| 亚洲国产色片| 日本一本二区三区精品| 岛国毛片在线播放| 日日啪夜夜撸| 亚洲内射少妇av| av国产久精品久网站免费入址| 亚洲成人一二三区av| 小蜜桃在线观看免费完整版高清| 人体艺术视频欧美日本| 日本wwww免费看| 成人午夜高清在线视频| 亚洲精品第二区| 天堂网av新在线| 99视频精品全部免费 在线| 欧美日韩精品成人综合77777| 欧美性感艳星| 国产麻豆成人av免费视频| 亚洲人成网站高清观看| 又黄又爽又刺激的免费视频.| 久久午夜福利片| 建设人人有责人人尽责人人享有的 | 丝袜喷水一区| 熟妇人妻久久中文字幕3abv| 免费观看av网站的网址| 精华霜和精华液先用哪个| 青春草亚洲视频在线观看| 亚洲av免费高清在线观看| 精品久久久久久久久久久久久| 国产白丝娇喘喷水9色精品| 亚洲成人精品中文字幕电影| 色综合站精品国产| 久久久色成人| 尤物成人国产欧美一区二区三区| 97超碰精品成人国产| 亚洲第一区二区三区不卡| 亚洲欧洲国产日韩| 最近中文字幕高清免费大全6| 99热全是精品| 国产在线男女| 日本黄大片高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久这里只有精品中国| 听说在线观看完整版免费高清| 老司机影院成人| videos熟女内射| 国产色爽女视频免费观看| 伦精品一区二区三区| 久久精品人妻少妇| 一级片'在线观看视频| 亚洲综合精品二区| 丝瓜视频免费看黄片| 日本熟妇午夜| 久久99蜜桃精品久久| 欧美另类一区| 国产美女午夜福利| 99久久精品一区二区三区| 国产精品一区二区性色av| 国产精品人妻久久久影院| 婷婷色综合www| av在线老鸭窝| 亚洲第一区二区三区不卡| 美女黄网站色视频| 一级毛片我不卡| 极品少妇高潮喷水抽搐| h日本视频在线播放| 久久精品国产自在天天线| 深爱激情五月婷婷| 免费在线观看成人毛片| 日韩欧美一区视频在线观看 | 亚州av有码| 在线观看免费高清a一片| 韩国高清视频一区二区三区| 国产免费又黄又爽又色| 欧美区成人在线视频| 日韩不卡一区二区三区视频在线| 一个人免费在线观看电影| 国产黄色视频一区二区在线观看| 精品久久久噜噜| 国产免费一级a男人的天堂| 精品久久久久久电影网| 国产av国产精品国产| 亚洲精品456在线播放app| 麻豆成人av视频| 国产一区二区三区av在线| 综合色丁香网| 91精品国产九色| 大香蕉久久网| 日韩精品有码人妻一区| 午夜福利在线观看吧| 天堂影院成人在线观看| 国产老妇女一区| 亚洲精品乱久久久久久| 日本一二三区视频观看| 国产精品麻豆人妻色哟哟久久 | 亚洲最大成人手机在线| 99热这里只有是精品50| 国产成人精品久久久久久| 九色成人免费人妻av| 人妻少妇偷人精品九色| 天堂影院成人在线观看| 又粗又硬又长又爽又黄的视频| 亚洲aⅴ乱码一区二区在线播放| 久久久久久久国产电影| 亚洲经典国产精华液单| 搡女人真爽免费视频火全软件| 美女内射精品一级片tv| 直男gayav资源| 熟女电影av网| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品日本国产第一区| 亚洲丝袜综合中文字幕| 男人舔奶头视频| 欧美一区二区亚洲| 亚洲成人一二三区av| www.av在线官网国产| 搡女人真爽免费视频火全软件| 国模一区二区三区四区视频| 十八禁国产超污无遮挡网站| 两个人的视频大全免费| 男女国产视频网站| 国产精品一二三区在线看| 熟妇人妻久久中文字幕3abv| 3wmmmm亚洲av在线观看| 97热精品久久久久久| 免费少妇av软件| 午夜日本视频在线| 国产综合精华液| 中文字幕免费在线视频6| 插逼视频在线观看| 亚洲熟妇中文字幕五十中出| 欧美xxxx性猛交bbbb| 亚洲无线观看免费| eeuss影院久久| 久久久久九九精品影院| 夜夜爽夜夜爽视频| 亚洲激情五月婷婷啪啪| 日本av手机在线免费观看| 国产91av在线免费观看| 国产精品久久久久久久电影| 一级片'在线观看视频| 日韩一区二区三区影片| 欧美激情久久久久久爽电影| 国产精品久久久久久精品电影| 免费观看的影片在线观看| 啦啦啦韩国在线观看视频| 男人舔奶头视频| 欧美zozozo另类| 久久久a久久爽久久v久久| 国产爱豆传媒在线观看| 偷拍熟女少妇极品色| 亚洲久久久久久中文字幕| 最近中文字幕2019免费版| 久久久国产一区二区| 国产亚洲一区二区精品| 色综合色国产| 国产 亚洲一区二区三区 | 亚洲av一区综合| 狠狠精品人妻久久久久久综合| 国产高清不卡午夜福利| 精品久久久久久久久av| 色视频www国产| 久久久久久久久久久免费av| 国产精品人妻久久久影院| 韩国高清视频一区二区三区| 亚洲经典国产精华液单| 欧美精品国产亚洲| 色视频www国产| 在线a可以看的网站| 国产亚洲午夜精品一区二区久久 | 午夜福利成人在线免费观看| 九九在线视频观看精品| 国产成人免费观看mmmm| 国产精品美女特级片免费视频播放器| 久久亚洲国产成人精品v| 菩萨蛮人人尽说江南好唐韦庄| 精品一区二区三区视频在线| 欧美日韩视频高清一区二区三区二| 午夜福利成人在线免费观看| 黄色一级大片看看| 一级毛片 在线播放| 免费在线观看成人毛片| 黄色一级大片看看| 我的女老师完整版在线观看| 一级毛片久久久久久久久女| 色吧在线观看| 免费看美女性在线毛片视频| 嘟嘟电影网在线观看| 欧美日韩综合久久久久久| 草草在线视频免费看| 少妇猛男粗大的猛烈进出视频 | 亚洲精品乱久久久久久| 乱系列少妇在线播放| www.av在线官网国产| 精品久久国产蜜桃| 精品99又大又爽又粗少妇毛片| 日本黄大片高清| 综合色丁香网| 日日啪夜夜爽| 中文天堂在线官网| av又黄又爽大尺度在线免费看| 99热6这里只有精品| 精品酒店卫生间| 精品一区在线观看国产| 精品久久久久久成人av| 大话2 男鬼变身卡| 校园人妻丝袜中文字幕| 99热这里只有是精品50| 欧美日韩在线观看h| 美女cb高潮喷水在线观看| a级一级毛片免费在线观看| 最新中文字幕久久久久| 水蜜桃什么品种好| av天堂中文字幕网| 狂野欧美激情性xxxx在线观看| 99久久人妻综合| 神马国产精品三级电影在线观看| 国模一区二区三区四区视频| 久久亚洲国产成人精品v| 丰满人妻一区二区三区视频av| 在线免费观看不下载黄p国产| 国产精品.久久久| 夫妻午夜视频| 伊人久久国产一区二区| 蜜臀久久99精品久久宅男| 美女内射精品一级片tv| 看黄色毛片网站| 最近最新中文字幕免费大全7| 超碰97精品在线观看| 国产免费一级a男人的天堂| 久久97久久精品| 午夜免费观看性视频| 日韩不卡一区二区三区视频在线| 少妇熟女欧美另类| 亚洲国产高清在线一区二区三| 嫩草影院新地址| 国产淫语在线视频| 九九久久精品国产亚洲av麻豆| 亚洲欧美一区二区三区国产| 成人亚洲精品一区在线观看 | 欧美变态另类bdsm刘玥| 草草在线视频免费看| 午夜福利在线观看吧| 精品久久久噜噜| 日本午夜av视频| 97精品久久久久久久久久精品| 国产精品一区二区三区四区免费观看| 国产精品人妻久久久影院| 蜜桃亚洲精品一区二区三区| av在线蜜桃| 在线观看免费高清a一片| 波多野结衣巨乳人妻| 日韩制服骚丝袜av| 国产亚洲一区二区精品| 国产成人精品久久久久久| 国精品久久久久久国模美| 中文字幕久久专区| 黄片wwwwww| av女优亚洲男人天堂| h日本视频在线播放| or卡值多少钱| 日本一本二区三区精品| 3wmmmm亚洲av在线观看| 狂野欧美白嫩少妇大欣赏| 插逼视频在线观看| 久久精品久久精品一区二区三区| 亚洲丝袜综合中文字幕| 亚洲国产av新网站| 国产亚洲一区二区精品| 一区二区三区高清视频在线| 91精品伊人久久大香线蕉| 91狼人影院| 久久人人爽人人爽人人片va| 色尼玛亚洲综合影院| 自拍偷自拍亚洲精品老妇| 免费观看的影片在线观看| 欧美激情久久久久久爽电影| 男女国产视频网站| 国产成人精品一,二区| 在线观看一区二区三区| 97精品久久久久久久久久精品| 肉色欧美久久久久久久蜜桃 | 一区二区三区免费毛片| 亚洲性久久影院| 亚洲熟妇中文字幕五十中出| 国产久久久一区二区三区| 亚洲自拍偷在线| 成人av在线播放网站| 国产熟女欧美一区二区| 校园人妻丝袜中文字幕| 中文资源天堂在线| 日本猛色少妇xxxxx猛交久久| 成年av动漫网址| 国产精品福利在线免费观看| 乱码一卡2卡4卡精品| 日韩一本色道免费dvd| 简卡轻食公司| 国内精品美女久久久久久| 一个人免费在线观看电影| 97超碰精品成人国产| 美女黄网站色视频| av在线老鸭窝| 黄片wwwwww| 国产精品伦人一区二区| 国产永久视频网站| 精品欧美国产一区二区三| 国产av码专区亚洲av| 亚洲精品影视一区二区三区av| av在线亚洲专区| 日韩强制内射视频| 欧美97在线视频| 国产成人精品婷婷| 丝袜喷水一区| 高清日韩中文字幕在线| 久久精品国产鲁丝片午夜精品| 2022亚洲国产成人精品| 亚洲四区av| 熟女电影av网| 能在线免费看毛片的网站| 一个人看的www免费观看视频| 免费观看无遮挡的男女| 国产探花极品一区二区| 精品久久国产蜜桃| 国产免费又黄又爽又色| 亚洲熟妇中文字幕五十中出| 亚洲三级黄色毛片| 国产成人福利小说| 国产探花在线观看一区二区| 国产精品一区二区三区四区久久| 久久久a久久爽久久v久久| 美女国产视频在线观看| 高清毛片免费看| 日本wwww免费看| 内地一区二区视频在线| 久久这里只有精品中国| 高清毛片免费看| 免费电影在线观看免费观看| 午夜精品一区二区三区免费看| 精品久久国产蜜桃| 国产成人午夜福利电影在线观看| 久久久久久国产a免费观看| 亚洲av中文字字幕乱码综合| 久久久精品94久久精品| 欧美高清成人免费视频www| 综合色av麻豆| 日韩成人伦理影院| 在线观看美女被高潮喷水网站| 特级一级黄色大片| 亚洲av电影在线观看一区二区三区 | 非洲黑人性xxxx精品又粗又长| 汤姆久久久久久久影院中文字幕 | 亚洲色图av天堂| 国模一区二区三区四区视频| 熟女电影av网| 久久久久网色| 91午夜精品亚洲一区二区三区| 免费黄色在线免费观看| 久久热精品热| 在线 av 中文字幕| ponron亚洲| 亚洲人与动物交配视频| 汤姆久久久久久久影院中文字幕 | 久久鲁丝午夜福利片| 日韩av在线大香蕉| 黄色日韩在线| 婷婷色综合www| 日本黄大片高清| 久99久视频精品免费| 亚洲精品亚洲一区二区| 国产亚洲一区二区精品| 亚洲欧美日韩无卡精品| 干丝袜人妻中文字幕| 亚洲成人中文字幕在线播放| 亚洲av电影不卡..在线观看| 久久这里只有精品中国| 午夜久久久久精精品| 亚洲精品乱久久久久久| 精品一区二区三卡| 免费电影在线观看免费观看| 亚洲在线自拍视频| 亚洲内射少妇av| 波多野结衣巨乳人妻| 国产精品精品国产色婷婷| 欧美+日韩+精品| 国产黄色小视频在线观看| 黄色一级大片看看| 欧美xxxx黑人xx丫x性爽| 肉色欧美久久久久久久蜜桃 | 国产精品国产三级国产专区5o| 一区二区三区高清视频在线| 久久人人爽人人片av| 少妇人妻一区二区三区视频| 国产精品久久久久久久久免| 热99在线观看视频| 99久久人妻综合| 又大又黄又爽视频免费| 亚洲天堂国产精品一区在线| 日韩欧美三级三区| 欧美激情在线99| 成人一区二区视频在线观看| 九色成人免费人妻av| 亚洲精品一二三| 日日干狠狠操夜夜爽| 欧美激情国产日韩精品一区| 色播亚洲综合网| 欧美区成人在线视频| 国精品久久久久久国模美| 国产又色又爽无遮挡免| 超碰av人人做人人爽久久| 人妻制服诱惑在线中文字幕| 国产精品久久久久久av不卡| av天堂中文字幕网| 六月丁香七月| 亚洲最大成人手机在线| 美女脱内裤让男人舔精品视频| 精品欧美国产一区二区三| 午夜福利高清视频| 亚洲精品日本国产第一区| 亚洲精品亚洲一区二区| 成年免费大片在线观看| 国产精品国产三级专区第一集| 人妻系列 视频| av免费观看日本| 日韩亚洲欧美综合| 国产极品天堂在线| 亚洲精品中文字幕在线视频 | 国产一区二区三区综合在线观看 | 舔av片在线| kizo精华| 卡戴珊不雅视频在线播放| 欧美最新免费一区二区三区| 好男人视频免费观看在线| 人妻夜夜爽99麻豆av| 秋霞在线观看毛片| 亚洲av国产av综合av卡| 又爽又黄无遮挡网站| 1000部很黄的大片| 免费看av在线观看网站| 男女下面进入的视频免费午夜| 亚洲不卡免费看| 777米奇影视久久| 成年女人在线观看亚洲视频 | 亚洲18禁久久av| 日日摸夜夜添夜夜添av毛片| 最近中文字幕高清免费大全6| 99re6热这里在线精品视频| 亚洲精品影视一区二区三区av| 亚洲成人av在线免费| av女优亚洲男人天堂|