• <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é)分析
    大码成人一级视频| 成人免费观看视频高清| 国产精品偷伦视频观看了| 免费久久久久久久精品成人欧美视频| 高清欧美精品videossex| 丝瓜视频免费看黄片| 亚洲欧美成人精品一区二区| 在线观看免费日韩欧美大片| 色视频在线一区二区三区| 亚洲欧美精品综合一区二区三区 | 免费看不卡的av| 亚洲av日韩在线播放| 精品一区二区三区四区五区乱码 | 欧美日韩视频精品一区| 制服人妻中文乱码| a级毛片黄视频| 狠狠精品人妻久久久久久综合| 国产伦理片在线播放av一区| 亚洲精品,欧美精品| 丰满乱子伦码专区| 97人妻天天添夜夜摸| 巨乳人妻的诱惑在线观看| 国产精品av久久久久免费| 久久久久精品久久久久真实原创| 男男h啪啪无遮挡| 香蕉精品网在线| 欧美+日韩+精品| 日产精品乱码卡一卡2卡三| 多毛熟女@视频| 亚洲精品日本国产第一区| 丁香六月天网| 一级毛片黄色毛片免费观看视频| 亚洲国产av影院在线观看| 中文字幕av电影在线播放| 国产成人欧美| 精品酒店卫生间| 国产熟女午夜一区二区三区| 亚洲成人手机| 精品亚洲成a人片在线观看| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 性色avwww在线观看| 高清在线视频一区二区三区| 午夜福利在线观看免费完整高清在| 久久久久久伊人网av| 亚洲国产欧美网| 精品午夜福利在线看| 校园人妻丝袜中文字幕| 久久久久久人人人人人| 熟女av电影| 汤姆久久久久久久影院中文字幕| 黄色怎么调成土黄色| 免费看av在线观看网站| 日韩熟女老妇一区二区性免费视频| 最黄视频免费看| 国产高清国产精品国产三级| 飞空精品影院首页| 成年人免费黄色播放视频| 午夜av观看不卡| 人人澡人人妻人| 性色av一级| 人妻一区二区av| 黄片无遮挡物在线观看| 亚洲色图综合在线观看| av在线app专区| 大陆偷拍与自拍| 国产97色在线日韩免费| 9热在线视频观看99| 在线观看免费视频网站a站| 欧美精品一区二区大全| 精品第一国产精品| 在线亚洲精品国产二区图片欧美| 国产精品av久久久久免费| 中国三级夫妇交换| 国产成人aa在线观看| 午夜福利一区二区在线看| 男女啪啪激烈高潮av片| 久久久久久久精品精品| 国产乱人偷精品视频| 久久国产精品大桥未久av| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美网| 丁香六月天网| 亚洲精品视频女| 国产人伦9x9x在线观看 | 久久精品久久久久久噜噜老黄| 日日摸夜夜添夜夜爱| 久久久久久人妻| 波野结衣二区三区在线| 国产精品免费大片| 婷婷色av中文字幕| 成年人免费黄色播放视频| 人人妻人人澡人人看| 国产免费一区二区三区四区乱码| 久久婷婷青草| 日韩人妻精品一区2区三区| 伦理电影免费视频| 欧美人与性动交α欧美软件| 九色亚洲精品在线播放| 97人妻天天添夜夜摸| 丝袜在线中文字幕| av在线app专区| 18+在线观看网站| 成人亚洲欧美一区二区av| 久久久久久久久久久久大奶| 99热国产这里只有精品6| av网站在线播放免费| 黄色毛片三级朝国网站| 欧美日韩精品网址| 黑人巨大精品欧美一区二区蜜桃| 热99久久久久精品小说推荐| 成人亚洲精品一区在线观看| 老鸭窝网址在线观看| 26uuu在线亚洲综合色| 国产精品 欧美亚洲| 嫩草影院入口| 日韩欧美一区视频在线观看| 欧美人与性动交α欧美精品济南到 | 国产亚洲最大av| 韩国高清视频一区二区三区| 最近中文字幕2019免费版| av又黄又爽大尺度在线免费看| 久久久久精品性色| 国产毛片在线视频| 一级毛片 在线播放| 国产老妇伦熟女老妇高清| 国产av一区二区精品久久| 久久99蜜桃精品久久| 看免费av毛片| 精品一区二区三卡| 一区在线观看完整版| 亚洲国产欧美网| 国产又爽黄色视频| 国产黄频视频在线观看| 一本久久精品| 亚洲,欧美精品.| 精品国产乱码久久久久久小说| 亚洲国产最新在线播放| 热99久久久久精品小说推荐| 性少妇av在线| 下体分泌物呈黄色| 伦理电影大哥的女人| 各种免费的搞黄视频| 狂野欧美激情性bbbbbb| 国产一区二区三区综合在线观看| 777米奇影视久久| 亚洲五月色婷婷综合| 在线天堂中文资源库| 亚洲av男天堂| 天堂俺去俺来也www色官网| 伦理电影大哥的女人| 精品久久蜜臀av无| 人妻系列 视频| 久久鲁丝午夜福利片| 美女国产视频在线观看| 欧美成人午夜精品| 叶爱在线成人免费视频播放| 777米奇影视久久| 国产精品一二三区在线看| 在线天堂中文资源库| 中国三级夫妇交换| 亚洲成人一二三区av| 亚洲av电影在线进入| 一边亲一边摸免费视频| 亚洲精品中文字幕在线视频| 99精国产麻豆久久婷婷| 精品一区在线观看国产| 国产在线一区二区三区精| 午夜福利网站1000一区二区三区| 黄色毛片三级朝国网站| 久久毛片免费看一区二区三区| 精品国产乱码久久久久久小说| 欧美 亚洲 国产 日韩一| 一区福利在线观看| 亚洲国产av影院在线观看| 欧美日韩av久久| 久久鲁丝午夜福利片| 国产精品人妻久久久影院| 日韩一区二区三区影片| 一级毛片黄色毛片免费观看视频| 看免费av毛片| 亚洲人成网站在线观看播放| 香蕉国产在线看| 黄色怎么调成土黄色| 国产精品嫩草影院av在线观看| 少妇猛男粗大的猛烈进出视频| 免费久久久久久久精品成人欧美视频| 国产熟女午夜一区二区三区| 99久久中文字幕三级久久日本| 97在线视频观看| 国产精品久久久久久av不卡| 欧美精品亚洲一区二区| 欧美激情 高清一区二区三区| 成人黄色视频免费在线看| 国产精品秋霞免费鲁丝片| 中文字幕av电影在线播放| 一级黄片播放器| 久久精品国产鲁丝片午夜精品| 精品国产乱码久久久久久男人| 观看av在线不卡| 卡戴珊不雅视频在线播放| 亚洲精品成人av观看孕妇| 99热全是精品| 香蕉精品网在线| 国产一区二区三区综合在线观看| 久久99蜜桃精品久久| 一区二区三区四区激情视频| 久久久久久人人人人人| 午夜福利,免费看| 一级,二级,三级黄色视频| 中文精品一卡2卡3卡4更新| 国产成人a∨麻豆精品| 伊人久久大香线蕉亚洲五| 五月开心婷婷网| 亚洲综合色惰| 国产免费一区二区三区四区乱码| 韩国精品一区二区三区| 最近中文字幕2019免费版| 国产成人精品婷婷| 国产老妇伦熟女老妇高清| 久久精品熟女亚洲av麻豆精品| 欧美xxⅹ黑人| 婷婷色综合大香蕉| 在线观看一区二区三区激情| 人妻 亚洲 视频| 黄片无遮挡物在线观看| 中国国产av一级| 侵犯人妻中文字幕一二三四区| 免费黄网站久久成人精品| 久久韩国三级中文字幕| 建设人人有责人人尽责人人享有的| 国产又色又爽无遮挡免| 国产精品偷伦视频观看了| 我的亚洲天堂| 下体分泌物呈黄色| 人人澡人人妻人| 狂野欧美激情性bbbbbb| 国产av一区二区精品久久| 街头女战士在线观看网站| www.精华液| 一级a爱视频在线免费观看| 亚洲精品美女久久久久99蜜臀 | 久久久久久人人人人人| 伦理电影大哥的女人| 又大又黄又爽视频免费| 国产1区2区3区精品| 在线亚洲精品国产二区图片欧美| 欧美最新免费一区二区三区| 国产熟女午夜一区二区三区| 久久99热这里只频精品6学生| 国产片特级美女逼逼视频| 丁香六月天网| 久久国产精品大桥未久av| 99热国产这里只有精品6| 欧美少妇被猛烈插入视频| 在线亚洲精品国产二区图片欧美| 波多野结衣av一区二区av| 边亲边吃奶的免费视频| 久久久国产一区二区| 中文字幕精品免费在线观看视频| 18禁国产床啪视频网站| 亚洲激情五月婷婷啪啪| 制服人妻中文乱码| a 毛片基地| tube8黄色片| 国产极品天堂在线| 丰满迷人的少妇在线观看| 天堂俺去俺来也www色官网| 一区福利在线观看| 欧美日韩成人在线一区二区| 9191精品国产免费久久| 乱人伦中国视频| 寂寞人妻少妇视频99o| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美一区二区三区久久| 天天躁夜夜躁狠狠躁躁| 欧美人与善性xxx| 哪个播放器可以免费观看大片| 久久97久久精品| 欧美av亚洲av综合av国产av | 国产免费福利视频在线观看| 97在线视频观看| 久热久热在线精品观看| 亚洲人成77777在线视频| 一边亲一边摸免费视频| 国产日韩欧美视频二区| 精品酒店卫生间| 最近最新中文字幕免费大全7| 一区二区av电影网| 99香蕉大伊视频| 夜夜骑夜夜射夜夜干| 赤兔流量卡办理| 欧美bdsm另类| 亚洲欧美清纯卡通| 午夜福利,免费看| 日韩三级伦理在线观看| 性色av一级| 久久久亚洲精品成人影院| 亚洲欧洲日产国产| 国产乱来视频区| 国产精品.久久久| 久久久精品94久久精品| 国产日韩欧美亚洲二区| 国产精品国产三级专区第一集| 欧美日韩国产mv在线观看视频| 人妻 亚洲 视频| 人人妻人人澡人人看| 国产精品免费视频内射| 久久久国产精品麻豆| 国产精品亚洲av一区麻豆 | 国产精品免费视频内射| 国产探花极品一区二区| 一级毛片 在线播放| 韩国高清视频一区二区三区| av有码第一页| 亚洲欧洲精品一区二区精品久久久 | 26uuu在线亚洲综合色| 最近2019中文字幕mv第一页| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 十八禁高潮呻吟视频| 欧美亚洲 丝袜 人妻 在线| 亚洲国产欧美网| 亚洲精华国产精华液的使用体验| 亚洲伊人久久精品综合| 一级a爱视频在线免费观看| 精品国产露脸久久av麻豆| 亚洲av电影在线观看一区二区三区| 国产成人欧美| 久久精品人人爽人人爽视色| av天堂久久9| 免费黄频网站在线观看国产| 欧美日韩成人在线一区二区| 欧美 亚洲 国产 日韩一| 亚洲中文av在线| 欧美精品av麻豆av| 女人高潮潮喷娇喘18禁视频| 91国产中文字幕| av不卡在线播放| 久久久久久久久久久免费av| 91久久精品国产一区二区三区| 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 亚洲国产精品成人久久小说| 日韩熟女老妇一区二区性免费视频| 大码成人一级视频| 免费黄网站久久成人精品| 男女边摸边吃奶| 午夜91福利影院| kizo精华| 国产一区二区在线观看av| 一级黄片播放器| 热99久久久久精品小说推荐| 国产麻豆69| 国产成人aa在线观看| 性高湖久久久久久久久免费观看| 可以免费在线观看a视频的电影网站 | 久久久久国产一级毛片高清牌| 少妇人妻精品综合一区二区| 伦理电影大哥的女人| 有码 亚洲区| 国产爽快片一区二区三区| 纵有疾风起免费观看全集完整版| 日韩成人av中文字幕在线观看| 波多野结衣av一区二区av| 麻豆av在线久日| 欧美日韩国产mv在线观看视频| 久久久精品免费免费高清| 日韩中字成人| 亚洲欧美中文字幕日韩二区| 久久鲁丝午夜福利片| 国产男人的电影天堂91| 国产精品蜜桃在线观看| 亚洲精品aⅴ在线观看| 日日爽夜夜爽网站| 久久这里有精品视频免费| 久久久久精品久久久久真实原创| 欧美日韩精品网址| 你懂的网址亚洲精品在线观看| 人妻人人澡人人爽人人| 在线看a的网站| 亚洲国产色片| 最新的欧美精品一区二区| 韩国av在线不卡| 欧美亚洲 丝袜 人妻 在线| 一区二区三区激情视频| 日本午夜av视频| 在线观看www视频免费| av.在线天堂| 校园人妻丝袜中文字幕| 成年动漫av网址| 婷婷色综合大香蕉| 久久人人爽人人片av| 人人妻人人澡人人爽人人夜夜| 国产亚洲午夜精品一区二区久久| 涩涩av久久男人的天堂| 国产一区二区三区av在线| 老司机亚洲免费影院| 99re6热这里在线精品视频| 国产激情久久老熟女| 精品亚洲乱码少妇综合久久| 一级,二级,三级黄色视频| 精品国产国语对白av| 亚洲在久久综合| av天堂久久9| 午夜福利视频在线观看免费| 91国产中文字幕| 9191精品国产免费久久| 亚洲精华国产精华液的使用体验| 日韩三级伦理在线观看| 亚洲国产欧美日韩在线播放| 免费观看av网站的网址| 色播在线永久视频| 精品人妻熟女毛片av久久网站| 色哟哟·www| 久久久a久久爽久久v久久| 一区二区三区精品91| 亚洲欧美一区二区三区久久| 久久av网站| 国产成人免费观看mmmm| 久热久热在线精品观看| 交换朋友夫妻互换小说| videossex国产| 不卡视频在线观看欧美| 色哟哟·www| 一个人免费看片子| 26uuu在线亚洲综合色| 激情视频va一区二区三区| 在线看a的网站| 久久精品久久精品一区二区三区| 精品一品国产午夜福利视频| 美女国产视频在线观看| 亚洲国产欧美日韩在线播放| av网站免费在线观看视频| 肉色欧美久久久久久久蜜桃| 国产综合精华液| 午夜免费观看性视频| 色吧在线观看| 久久久久久人人人人人| 久久99精品国语久久久| 最近的中文字幕免费完整| 只有这里有精品99| 欧美成人午夜免费资源| 欧美精品国产亚洲| 亚洲第一青青草原| 中文字幕av电影在线播放| 最近中文字幕2019免费版| www日本在线高清视频| 一级片免费观看大全| 国产乱来视频区| 欧美bdsm另类| 国产成人精品一,二区| 免费av中文字幕在线| 深夜精品福利| 色播在线永久视频| 永久网站在线| 久久久久网色| 成人手机av| 久久久久久久国产电影| 少妇人妻久久综合中文| 美女中出高潮动态图| 搡老乐熟女国产| 欧美日本中文国产一区发布| 婷婷色综合www| 美女脱内裤让男人舔精品视频| 黄色毛片三级朝国网站| 国产免费福利视频在线观看| 日韩中文字幕欧美一区二区 | 日本-黄色视频高清免费观看| 国产精品99久久99久久久不卡 | 在线观看国产h片| 老汉色av国产亚洲站长工具| 美国免费a级毛片| 天天躁夜夜躁狠狠久久av| 老汉色∧v一级毛片| 亚洲综合色惰| 国产精品国产av在线观看| 国精品久久久久久国模美| 中文字幕人妻丝袜制服| 亚洲三区欧美一区| 国产一区二区在线观看av| 国产精品香港三级国产av潘金莲 | 欧美中文综合在线视频| 国产一区二区三区综合在线观看| 日韩免费高清中文字幕av| 男女边摸边吃奶| 纯流量卡能插随身wifi吗| 亚洲成国产人片在线观看| 免费看不卡的av| 少妇猛男粗大的猛烈进出视频| 寂寞人妻少妇视频99o| 国产亚洲一区二区精品| 亚洲av.av天堂| 国产午夜精品一二区理论片| 美女福利国产在线| 人成视频在线观看免费观看| 国产无遮挡羞羞视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 在线观看免费视频网站a站| 午夜免费鲁丝| 日本免费在线观看一区| 黄网站色视频无遮挡免费观看| 精品国产一区二区久久| 成人午夜精彩视频在线观看| 国产免费视频播放在线视频| 男的添女的下面高潮视频| 欧美 日韩 精品 国产| 成年女人毛片免费观看观看9 | 国产精品久久久久久久久免| 精品久久蜜臀av无| 深夜精品福利| 欧美人与性动交α欧美软件| 黄片无遮挡物在线观看| tube8黄色片| 黑人欧美特级aaaaaa片| 日本免费在线观看一区| 热re99久久精品国产66热6| 七月丁香在线播放| 欧美激情 高清一区二区三区| 国产免费一区二区三区四区乱码| 亚洲精品视频女| 久久精品国产亚洲av天美| 高清av免费在线| 免费在线观看完整版高清| 国语对白做爰xxxⅹ性视频网站| freevideosex欧美| 高清在线视频一区二区三区| 国产高清国产精品国产三级| 精品国产一区二区三区久久久樱花| 亚洲五月色婷婷综合| 亚洲国产av影院在线观看| 国产午夜精品一二区理论片| videossex国产| 最新的欧美精品一区二区| 十八禁高潮呻吟视频| 成人影院久久| 国产有黄有色有爽视频| 性色av一级| 成年女人毛片免费观看观看9 | 久久久久久人妻| √禁漫天堂资源中文www| 成人免费观看视频高清| 最近的中文字幕免费完整| 日本wwww免费看| 国产av一区二区精品久久| 午夜日韩欧美国产| 久久久亚洲精品成人影院| 久久精品国产亚洲av天美| 亚洲欧洲精品一区二区精品久久久 | 国产又色又爽无遮挡免| 咕卡用的链子| 免费人妻精品一区二区三区视频| 免费黄频网站在线观看国产| 我的亚洲天堂| 成人毛片60女人毛片免费| 中文字幕制服av| 国产亚洲午夜精品一区二区久久| 亚洲经典国产精华液单| 亚洲欧美日韩另类电影网站| 性少妇av在线| 亚洲,欧美,日韩| www日本在线高清视频| 看非洲黑人一级黄片| 秋霞伦理黄片| 亚洲欧美一区二区三区久久| 一区二区三区精品91| 一级爰片在线观看| 国产av国产精品国产| 久久免费观看电影| 涩涩av久久男人的天堂| 综合色丁香网| 欧美精品人与动牲交sv欧美| 精品亚洲乱码少妇综合久久| 久久久久久久久免费视频了| 亚洲欧美精品自产自拍| 18禁观看日本| av天堂久久9| 成年美女黄网站色视频大全免费| 亚洲国产成人一精品久久久| www.熟女人妻精品国产| 香蕉丝袜av| 国产成人欧美| 欧美国产精品一级二级三级| 久久久久久久久久人人人人人人| 久久久精品区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 春色校园在线视频观看| 国产熟女午夜一区二区三区| 久热这里只有精品99| 一级毛片我不卡| 波多野结衣一区麻豆| 成年女人在线观看亚洲视频| av国产久精品久网站免费入址| 欧美精品一区二区免费开放| 日韩精品免费视频一区二区三区| 精品卡一卡二卡四卡免费| 黄片无遮挡物在线观看| 久久精品熟女亚洲av麻豆精品| 国产免费一区二区三区四区乱码| 欧美精品人与动牲交sv欧美| 免费看av在线观看网站| 极品人妻少妇av视频| 欧美97在线视频| 丰满迷人的少妇在线观看| 九草在线视频观看| 欧美日韩国产mv在线观看视频| 美女大奶头黄色视频| 天堂俺去俺来也www色官网| 成人国产麻豆网| 国产淫语在线视频| 日本av手机在线免费观看| 香蕉丝袜av| 欧美成人午夜精品| 精品国产一区二区久久| 日韩精品免费视频一区二区三区|