• <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é)分析
    男女那种视频在线观看| 美女 人体艺术 gogo| 天堂动漫精品| 久久99热这里只有精品18| 午夜免费成人在线视频| 久久久久国产精品人妻aⅴ院| 久久人妻av系列| 国产成人av教育| 午夜亚洲福利在线播放| 欧美性猛交黑人性爽| 午夜福利18| 嫩草影视91久久| 国产黄片美女视频| 日本熟妇午夜| 日日摸夜夜添夜夜添小说| 黄色片一级片一级黄色片| 国产亚洲欧美在线一区二区| 日韩欧美国产一区二区入口| 亚洲第一电影网av| 亚洲人与动物交配视频| 亚洲精品久久国产高清桃花| 亚洲精品色激情综合| 日韩成人在线观看一区二区三区| 韩国av一区二区三区四区| 制服人妻中文乱码| 精品无人区乱码1区二区| 亚洲精品粉嫩美女一区| 日韩精品青青久久久久久| 中文资源天堂在线| 最新在线观看一区二区三区| 91成年电影在线观看| 久久人人精品亚洲av| 在线十欧美十亚洲十日本专区| 成人18禁在线播放| 男插女下体视频免费在线播放| 老汉色av国产亚洲站长工具| 午夜福利欧美成人| 国产精品国产高清国产av| 窝窝影院91人妻| 免费搜索国产男女视频| 一进一出好大好爽视频| 一区二区三区激情视频| 亚洲成人精品中文字幕电影| 亚洲乱码一区二区免费版| 免费在线观看视频国产中文字幕亚洲| 少妇的丰满在线观看| 久久天躁狠狠躁夜夜2o2o| 久久国产乱子伦精品免费另类| 亚洲av中文字字幕乱码综合| 51午夜福利影视在线观看| 亚洲国产精品sss在线观看| 麻豆成人av在线观看| 午夜日韩欧美国产| 怎么达到女性高潮| 久久久久久久久免费视频了| 少妇的丰满在线观看| 中国美女看黄片| 91成年电影在线观看| 淫秽高清视频在线观看| 亚洲自拍偷在线| 日本免费一区二区三区高清不卡| 动漫黄色视频在线观看| 在线十欧美十亚洲十日本专区| 欧美一级a爱片免费观看看 | 1024视频免费在线观看| 久久久国产精品麻豆| 亚洲男人天堂网一区| 人妻夜夜爽99麻豆av| 欧美日韩国产亚洲二区| 看黄色毛片网站| 日本撒尿小便嘘嘘汇集6| 色精品久久人妻99蜜桃| 亚洲av电影在线进入| 午夜福利免费观看在线| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲熟女毛片儿| 亚洲中文日韩欧美视频| 99久久综合精品五月天人人| a级毛片在线看网站| 久久久久久久久久黄片| 亚洲av成人不卡在线观看播放网| 看黄色毛片网站| 色综合欧美亚洲国产小说| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av熟女| 成在线人永久免费视频| or卡值多少钱| 色综合站精品国产| 丁香欧美五月| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三区视频在线观看免费| 一边摸一边做爽爽视频免费| 国产视频一区二区在线看| 国产成人av激情在线播放| 一进一出抽搐动态| 狠狠狠狠99中文字幕| 欧美性猛交╳xxx乱大交人| 亚洲人成网站在线播放欧美日韩| 国产成人精品久久二区二区免费| 一进一出抽搐动态| 毛片女人毛片| 又大又爽又粗| 亚洲,欧美精品.| 他把我摸到了高潮在线观看| 日本 av在线| 国产蜜桃级精品一区二区三区| 精品乱码久久久久久99久播| 日韩精品青青久久久久久| 老司机在亚洲福利影院| 天堂av国产一区二区熟女人妻 | 国产激情久久老熟女| 久久婷婷成人综合色麻豆| 黄色女人牲交| 好男人电影高清在线观看| 日本成人三级电影网站| 色综合婷婷激情| 伦理电影免费视频| 老司机午夜十八禁免费视频| 久久性视频一级片| 免费搜索国产男女视频| 欧美性猛交╳xxx乱大交人| 哪里可以看免费的av片| 1024香蕉在线观看| 18禁黄网站禁片免费观看直播| 久久精品人妻少妇| 999久久久精品免费观看国产| 日韩免费av在线播放| 此物有八面人人有两片| 人人妻人人澡欧美一区二区| 国内精品一区二区在线观看| 国产精品亚洲美女久久久| 岛国在线观看网站| 亚洲成人久久性| 手机成人av网站| 欧美精品亚洲一区二区| 麻豆久久精品国产亚洲av| 中文资源天堂在线| 亚洲av成人精品一区久久| 亚洲av成人精品一区久久| 欧美一级毛片孕妇| 国产久久久一区二区三区| www.熟女人妻精品国产| 国产69精品久久久久777片 | 性色av乱码一区二区三区2| 欧美性猛交黑人性爽| 99在线人妻在线中文字幕| videosex国产| 国产成人精品无人区| 禁无遮挡网站| 精品一区二区三区四区五区乱码| 亚洲av电影不卡..在线观看| tocl精华| 麻豆av在线久日| 国产黄色小视频在线观看| 亚洲精品粉嫩美女一区| 精品国产亚洲在线| 久久久久久九九精品二区国产 | 亚洲片人在线观看| 国产野战对白在线观看| 国内少妇人妻偷人精品xxx网站 | 特级一级黄色大片| 9191精品国产免费久久| 国产99久久九九免费精品| 少妇熟女aⅴ在线视频| www.自偷自拍.com| 少妇粗大呻吟视频| 99国产极品粉嫩在线观看| 女人高潮潮喷娇喘18禁视频| 久久精品国产亚洲av香蕉五月| 女人高潮潮喷娇喘18禁视频| 欧美av亚洲av综合av国产av| 久久这里只有精品中国| 男女床上黄色一级片免费看| or卡值多少钱| 99国产精品一区二区蜜桃av| 亚洲av成人av| 91麻豆精品激情在线观看国产| 久久婷婷人人爽人人干人人爱| 99re在线观看精品视频| 免费看a级黄色片| 久久欧美精品欧美久久欧美| 99久久久亚洲精品蜜臀av| 亚洲电影在线观看av| 久热爱精品视频在线9| 欧美日韩精品网址| 国产精品一区二区三区四区免费观看 | 日韩欧美国产一区二区入口| 午夜激情av网站| 两个人的视频大全免费| 波多野结衣高清无吗| 国产精品香港三级国产av潘金莲| 99热只有精品国产| 久久久久久国产a免费观看| 国产黄a三级三级三级人| 亚洲激情在线av| 久久久久久久久久黄片| 男女视频在线观看网站免费 | 看黄色毛片网站| 一二三四在线观看免费中文在| 午夜免费成人在线视频| 19禁男女啪啪无遮挡网站| 真人一进一出gif抽搐免费| 欧美另类亚洲清纯唯美| 熟妇人妻久久中文字幕3abv| 国产亚洲精品av在线| 中文资源天堂在线| 国产午夜精品论理片| 亚洲一区二区三区不卡视频| 亚洲av片天天在线观看| 国产亚洲精品一区二区www| 成人高潮视频无遮挡免费网站| 成人国产综合亚洲| 亚洲欧美精品综合一区二区三区| 激情在线观看视频在线高清| 欧美三级亚洲精品| 神马国产精品三级电影在线观看 | 亚洲黑人精品在线| 国产亚洲av嫩草精品影院| 欧美极品一区二区三区四区| 国产1区2区3区精品| 久久久久久免费高清国产稀缺| 男女做爰动态图高潮gif福利片| 久久天躁狠狠躁夜夜2o2o| 国产精品久久电影中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 一边摸一边做爽爽视频免费| 女人高潮潮喷娇喘18禁视频| 久久精品aⅴ一区二区三区四区| 国产主播在线观看一区二区| 免费高清视频大片| 日日干狠狠操夜夜爽| 叶爱在线成人免费视频播放| 九色成人免费人妻av| 性色av乱码一区二区三区2| 国产精华一区二区三区| 18禁观看日本| 日韩精品青青久久久久久| 亚洲中文av在线| www.999成人在线观看| 黄频高清免费视频| 亚洲人成77777在线视频| 高清毛片免费观看视频网站| 大型av网站在线播放| 18禁裸乳无遮挡免费网站照片| 国产精品 欧美亚洲| 亚洲av五月六月丁香网| 五月伊人婷婷丁香| 久久精品国产亚洲av高清一级| 99久久综合精品五月天人人| 免费在线观看视频国产中文字幕亚洲| 午夜影院日韩av| 淫秽高清视频在线观看| 午夜激情福利司机影院| 中文在线观看免费www的网站 | 午夜a级毛片| 日韩高清综合在线| 国产精品久久电影中文字幕| 不卡av一区二区三区| 日韩免费av在线播放| 亚洲成人国产一区在线观看| 这个男人来自地球电影免费观看| 又爽又黄无遮挡网站| 真人做人爱边吃奶动态| 亚洲国产高清在线一区二区三| 97碰自拍视频| АⅤ资源中文在线天堂| 桃红色精品国产亚洲av| www.熟女人妻精品国产| 午夜激情av网站| 国产精品久久久久久精品电影| 听说在线观看完整版免费高清| 国产成人av教育| 亚洲国产欧美网| 日韩欧美国产一区二区入口| 嫩草影院精品99| 午夜精品久久久久久毛片777| 国产精品永久免费网站| 搡老熟女国产l中国老女人| 麻豆久久精品国产亚洲av| 99精品欧美一区二区三区四区| 成人手机av| 精品国产乱子伦一区二区三区| 久久精品亚洲精品国产色婷小说| 欧美三级亚洲精品| 欧美午夜高清在线| 一个人观看的视频www高清免费观看 | 嫩草影视91久久| 91成年电影在线观看| 亚洲中文av在线| 一本一本综合久久| 色综合欧美亚洲国产小说| 91麻豆av在线| 欧美3d第一页| 91老司机精品| 欧美精品啪啪一区二区三区| 人妻夜夜爽99麻豆av| 日本 av在线| 欧美一区二区国产精品久久精品 | 国产成人aa在线观看| 正在播放国产对白刺激| 国产亚洲精品久久久久久毛片| 美女免费视频网站| 久久中文字幕人妻熟女| 亚洲欧美日韩高清专用| 亚洲欧洲精品一区二区精品久久久| 欧美激情久久久久久爽电影| 美女午夜性视频免费| 欧美性猛交╳xxx乱大交人| 欧美日韩一级在线毛片| 国产人伦9x9x在线观看| 免费在线观看影片大全网站| 久久中文看片网| 欧美黑人巨大hd| 日韩欧美在线二视频| 麻豆成人av在线观看| 岛国视频午夜一区免费看| 国产熟女午夜一区二区三区| 波多野结衣巨乳人妻| 亚洲精品在线观看二区| 床上黄色一级片| 最近在线观看免费完整版| 国产精品 欧美亚洲| 国产精品免费视频内射| 国产主播在线观看一区二区| 欧美zozozo另类| 色综合欧美亚洲国产小说| 久热爱精品视频在线9| 国产av不卡久久| 在线看三级毛片| 欧美高清成人免费视频www| 日韩欧美国产一区二区入口| 亚洲熟妇中文字幕五十中出| 亚洲性夜色夜夜综合| 国产亚洲av嫩草精品影院| 欧美一级毛片孕妇| 少妇熟女aⅴ在线视频| 在线视频色国产色| 久久中文字幕人妻熟女| 欧美3d第一页| 国产精品免费视频内射| a级毛片a级免费在线| 最新在线观看一区二区三区| 中文资源天堂在线| 99国产极品粉嫩在线观看| 午夜福利免费观看在线| 国产成人欧美在线观看| 国内精品一区二区在线观看| 岛国视频午夜一区免费看| 亚洲av第一区精品v没综合| 精品第一国产精品| 亚洲全国av大片| 欧美在线一区亚洲| 亚洲人成网站高清观看| 变态另类成人亚洲欧美熟女| 无遮挡黄片免费观看| 男女之事视频高清在线观看| 亚洲第一欧美日韩一区二区三区| 色尼玛亚洲综合影院| 国产精品一区二区免费欧美| 91在线观看av| 国产精品国产高清国产av| 99国产极品粉嫩在线观看| 麻豆国产97在线/欧美 | 成年人黄色毛片网站| 99热这里只有精品一区 | 一级片免费观看大全| 婷婷丁香在线五月| 久久久久久久久中文| 亚洲精品在线美女| 男人舔女人下体高潮全视频| 露出奶头的视频| 国产成人av教育| 看免费av毛片| 欧美成人午夜精品| 国产精品香港三级国产av潘金莲| 18禁美女被吸乳视频| x7x7x7水蜜桃| 成人精品一区二区免费| 每晚都被弄得嗷嗷叫到高潮| 香蕉av资源在线| 日韩大码丰满熟妇| 精品国产乱子伦一区二区三区| 熟女电影av网| 欧美黑人欧美精品刺激| 青草久久国产| 亚洲在线自拍视频| 国产午夜福利久久久久久| 99精品欧美一区二区三区四区| 青草久久国产| 国产成人av激情在线播放| 日本黄色视频三级网站网址| 一二三四社区在线视频社区8| 高清毛片免费观看视频网站| 午夜成年电影在线免费观看| 亚洲一区中文字幕在线| www日本黄色视频网| 在线观看美女被高潮喷水网站 | 国产精品免费一区二区三区在线| 欧美在线一区亚洲| 午夜成年电影在线免费观看| 国产精品98久久久久久宅男小说| 亚洲片人在线观看| 日本 欧美在线| svipshipincom国产片| 在线观看66精品国产| 欧美黑人精品巨大| 欧美黑人巨大hd| 欧美国产日韩亚洲一区| 麻豆成人av在线观看| 午夜免费观看网址| 黄片大片在线免费观看| xxx96com| 成人av一区二区三区在线看| 大型av网站在线播放| 国产97色在线日韩免费| 一本大道久久a久久精品| 在线视频色国产色| 欧美成人性av电影在线观看| 色av中文字幕| 婷婷精品国产亚洲av| 中文字幕人妻丝袜一区二区| 给我免费播放毛片高清在线观看| 精品久久久久久久久久免费视频| 国产99白浆流出| 美女高潮喷水抽搐中文字幕| 成年人黄色毛片网站| 国产探花在线观看一区二区| 男人舔女人的私密视频| 香蕉丝袜av| 特大巨黑吊av在线直播| 亚洲色图av天堂| 色综合亚洲欧美另类图片| 俺也久久电影网| 夜夜躁狠狠躁天天躁| 日韩免费av在线播放| 亚洲七黄色美女视频| 久9热在线精品视频| 少妇的丰满在线观看| 久久精品国产99精品国产亚洲性色| 国产成人啪精品午夜网站| 日韩中文字幕欧美一区二区| 国产黄a三级三级三级人| 国产成人av教育| 韩国av一区二区三区四区| 精华霜和精华液先用哪个| 波多野结衣巨乳人妻| 日日爽夜夜爽网站| 97碰自拍视频| 极品教师在线免费播放| 欧美成人免费av一区二区三区| 亚洲国产欧美网| 少妇人妻一区二区三区视频| 亚洲七黄色美女视频| 国产高清videossex| 男人舔奶头视频| 99久久无色码亚洲精品果冻| 国产在线精品亚洲第一网站| 国内少妇人妻偷人精品xxx网站 | 天天躁狠狠躁夜夜躁狠狠躁| 久久精品91无色码中文字幕| 久久 成人 亚洲| 亚洲国产看品久久| 别揉我奶头~嗯~啊~动态视频| 国产免费男女视频| 国产精品,欧美在线| 97超级碰碰碰精品色视频在线观看| 国产精品亚洲一级av第二区| 露出奶头的视频| 在线播放国产精品三级| 神马国产精品三级电影在线观看 | 欧美高清成人免费视频www| 99久久精品热视频| 一级作爱视频免费观看| 99国产精品一区二区三区| 国产欧美日韩精品亚洲av| 蜜桃久久精品国产亚洲av| ponron亚洲| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品色激情综合| 精品午夜福利视频在线观看一区| 不卡一级毛片| 中亚洲国语对白在线视频| 午夜日韩欧美国产| 国产欧美日韩一区二区三| 国产区一区二久久| 夜夜看夜夜爽夜夜摸| 少妇裸体淫交视频免费看高清 | 老熟妇乱子伦视频在线观看| 90打野战视频偷拍视频| 一本精品99久久精品77| 国产男靠女视频免费网站| 中国美女看黄片| 在线观看舔阴道视频| 国产精品一及| 首页视频小说图片口味搜索| 欧美日韩黄片免| 国产高清有码在线观看视频 | 变态另类成人亚洲欧美熟女| 黄色女人牲交| 18禁黄网站禁片午夜丰满| 久久香蕉精品热| 国产精品精品国产色婷婷| 免费在线观看日本一区| 国产一区二区在线观看日韩 | 日韩欧美在线乱码| 美女高潮喷水抽搐中文字幕| 成年女人毛片免费观看观看9| 国内少妇人妻偷人精品xxx网站 | 亚洲aⅴ乱码一区二区在线播放 | 床上黄色一级片| 黄色丝袜av网址大全| 午夜免费成人在线视频| 麻豆国产av国片精品| 一二三四在线观看免费中文在| 欧美性猛交╳xxx乱大交人| 啦啦啦观看免费观看视频高清| 一级a爱片免费观看的视频| 狂野欧美白嫩少妇大欣赏| 三级国产精品欧美在线观看 | 一级片免费观看大全| 国产午夜精品论理片| 久久精品综合一区二区三区| 观看免费一级毛片| 国内揄拍国产精品人妻在线| 久久久国产成人精品二区| 午夜免费成人在线视频| 一级作爱视频免费观看| 在线观看www视频免费| 免费搜索国产男女视频| 91九色精品人成在线观看| 久久久水蜜桃国产精品网| www.熟女人妻精品国产| avwww免费| 成人国语在线视频| 波多野结衣高清无吗| 午夜福利视频1000在线观看| 成人永久免费在线观看视频| 亚洲av日韩精品久久久久久密| 中文在线观看免费www的网站 | 日韩欧美免费精品| 日本一二三区视频观看| 亚洲av成人av| 久久久久久大精品| 欧美成人一区二区免费高清观看 | 国内精品久久久久精免费| 国产aⅴ精品一区二区三区波| 国产精品免费视频内射| 国产成人精品无人区| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 色尼玛亚洲综合影院| 国产不卡一卡二| 女人被狂操c到高潮| 男女午夜视频在线观看| 亚洲成av人片免费观看| 在线a可以看的网站| 国产高清有码在线观看视频 | 亚洲成人久久性| 国产亚洲av嫩草精品影院| 国产视频一区二区在线看| 欧美zozozo另类| 欧美绝顶高潮抽搐喷水| 欧美日韩乱码在线| 级片在线观看| 国产亚洲精品第一综合不卡| e午夜精品久久久久久久| 精品电影一区二区在线| 好男人在线观看高清免费视频| 中文亚洲av片在线观看爽| 亚洲国产欧美一区二区综合| 日本熟妇午夜| 国产黄色小视频在线观看| 日韩国内少妇激情av| 在线观看美女被高潮喷水网站 | 亚洲人成伊人成综合网2020| 午夜视频精品福利| 国产高清有码在线观看视频 | 国产男靠女视频免费网站| 美女 人体艺术 gogo| 精品电影一区二区在线| 国产成人av教育| 国产爱豆传媒在线观看 | 欧美精品亚洲一区二区| 母亲3免费完整高清在线观看| 国产亚洲精品一区二区www| 久久精品综合一区二区三区| 国产精品电影一区二区三区| 日韩欧美三级三区| 成人18禁在线播放| 草草在线视频免费看| 日本在线视频免费播放| 国产私拍福利视频在线观看| 蜜桃久久精品国产亚洲av| 午夜福利欧美成人| 成人三级做爰电影| 精华霜和精华液先用哪个| 久久精品国产99精品国产亚洲性色| 亚洲国产精品999在线| 久久久精品国产亚洲av高清涩受| 精品电影一区二区在线| 亚洲人成伊人成综合网2020| 午夜福利在线在线| 成人国语在线视频| 国产三级在线视频| 久久久久性生活片| 亚洲 国产 在线| xxxwww97欧美| 国产在线观看jvid| 婷婷丁香在线五月| 一区福利在线观看| www.精华液| 日韩精品免费视频一区二区三区| 成熟少妇高潮喷水视频| 欧美黑人精品巨大| 国产久久久一区二区三区|