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

    淺海天然氣管道泄漏擴(kuò)散過(guò)程模擬研究

    2022-06-11 02:38:42李延豪侯磊于巧燕柴沖肖開(kāi)喜
    關(guān)鍵詞:氣團(tuán)水流氣泡

    李延豪,侯磊,于巧燕,柴沖,肖開(kāi)喜

    淺海天然氣管道泄漏擴(kuò)散過(guò)程模擬研究

    李延豪1,2,侯磊1,2,于巧燕1,2,柴沖1,2,肖開(kāi)喜1,2

    (1.中國(guó)石油大學(xué)(北京) 機(jī)械與儲(chǔ)運(yùn)工程學(xué)院,北京 102249; 2.中國(guó)石油大學(xué)(北京) 油氣管道輸送安全國(guó)家工程實(shí)驗(yàn)室,北京,102249)

    為評(píng)價(jià)淺海海底天然氣管道泄漏事故后果,根據(jù)計(jì)算流體力學(xué)與多相流動(dòng)理論,針對(duì)國(guó)外某天然氣管道海峽穿越段,建立淺海海底管道泄漏擴(kuò)散過(guò)程的計(jì)算模型。將泄漏孔徑、泄漏速率、水流速度3個(gè)主要影響因素作為條件變量,模擬不同情況下的氣液兩相運(yùn)動(dòng)過(guò)程。結(jié)果表明,水下氣體擴(kuò)散分為三個(gè)階段,即泄漏口上方形成氣團(tuán)、氣團(tuán)呈蘑菇狀上升、氣團(tuán)由大氣泡分裂為小氣泡;泄漏孔徑和泄漏速率對(duì)水下氣體擴(kuò)散到水面的時(shí)間具有顯著影響,泄漏孔徑與泄漏速率越大,氣體泄漏量越大;氣體泄漏量越大,水下氣團(tuán)體積越大,到達(dá)水面的時(shí)間越短;水流速度顯著影響氣體的擴(kuò)散軌跡,水流速度越大,氣體運(yùn)動(dòng)軌跡與海底的夾角越小,沿海流方向擴(kuò)散的距離越遠(yuǎn)。研究結(jié)果可為水下天然氣管道泄漏事故應(yīng)急處理提供一定的科學(xué)指導(dǎo)。

    海底管道; 氣體泄漏; 多相流動(dòng); 計(jì)算流體力學(xué)

    海底管道具有輸量大、效率高的優(yōu)點(diǎn),是海上油氣最主要的輸送方式。不同于陸上埋地管道,海底管道每時(shí)每刻都遭受著海流的沖刷及海底地震、臺(tái)風(fēng)等自然條件的威脅,加之海水腐蝕、第三方破壞等因素,海底管道的安全性和完整性面臨巨大的考驗(yàn)[1]。海底天然氣管道一旦發(fā)生泄漏,其擴(kuò)散至海面的可燃?xì)怏w有可能導(dǎo)致爆炸、火災(zāi)等重大事故,嚴(yán)重威脅海上工作人員和生產(chǎn)設(shè)施的安全[2]。因此,掌握海底管道天然氣泄漏擴(kuò)散規(guī)律,確定氣液兩相在不同環(huán)境條件下的運(yùn)動(dòng)分布,對(duì)海底天然氣管道的風(fēng)險(xiǎn)防范管控及事故應(yīng)急處理具有重要意義。

    目前,關(guān)于水下氣體泄漏擴(kuò)散的研究多基于氣泡羽流理論[3?5]。國(guó)內(nèi)學(xué)者多數(shù)以氣泡羽流理論模型為參考,運(yùn)用有限元軟件進(jìn)行仿真,并輔以一定的實(shí)驗(yàn)驗(yàn)證,對(duì)已有的研究成果進(jìn)行重現(xiàn)及進(jìn)一步完善。景海泳等[6]采用VOF多相流模型和DPM離散相模型耦合的方式進(jìn)行數(shù)值模擬,得到了氣泡粒子的擴(kuò)散半徑和軸向位移,但模擬氣體的泄漏量較小,沒(méi)有考慮到水流速度對(duì)水中氣泡運(yùn)動(dòng)的影響,也未考慮水下氣體的可壓縮性。文闖等[7]采用VOF多相流模型進(jìn)行數(shù)值模擬,考慮了水下泄漏氣體的可壓縮性,通過(guò)實(shí)驗(yàn)驗(yàn)證了模擬結(jié)果的準(zhǔn)確性和可靠性,但是仍以靜水條件作為模擬設(shè)置,不符合海底管道實(shí)際工作環(huán)境的特點(diǎn)。竇梓元[8]在數(shù)值模擬的基礎(chǔ)上,通過(guò)無(wú)量綱分析得出泄漏氣體的上浮時(shí)間、水平遷移距離與管道壓力、水流速度、泄漏孔徑的無(wú)量綱關(guān)系。李新宏等[9]對(duì)水下氣泡的運(yùn)動(dòng)軌跡進(jìn)行跟蹤求解,并以水下氣體擴(kuò)散結(jié)果作為依據(jù),分析了天然氣在大氣中的擴(kuò)散。王少雄等[10]通過(guò)建立三維CFD模型和積分?jǐn)?shù)學(xué)模型,研究了水下輸氣管道泄漏擴(kuò)散特性,并對(duì)兩種模型預(yù)測(cè)結(jié)果的準(zhǔn)確性進(jìn)行了對(duì)比評(píng)估。綜上所述,國(guó)內(nèi)已有研究對(duì)不同條件下天然氣管道泄漏氣體運(yùn)動(dòng)分布規(guī)律的多因素分析仍不完善。因此,針對(duì)國(guó)外某天然氣管道海峽穿越段,本研究以泄漏口孔徑、泄漏速率和水流速度為變量,運(yùn)用Fluent軟件對(duì)水下泄漏擴(kuò)散氣體的時(shí)空分布進(jìn)行有限元數(shù)值模擬,分析不同條件下淺海海底管道泄漏氣體的擴(kuò)散規(guī)律。

    1 海底天然氣管道泄漏擴(kuò)散理論模型

    1.1 基本控制方程

    流體流動(dòng)遵循質(zhì)量、動(dòng)量和能量守恒定律,對(duì)這3個(gè)守恒方程進(jìn)行整合,其表達(dá)式如式(1)所示[11]。

    1.2 湍流模型

    選用標(biāo)準(zhǔn)?湍流模型對(duì)氣液兩相運(yùn)動(dòng)過(guò)程進(jìn)行求解,其輸運(yùn)方程為[12]:

    1.3 VOF模型

    VOF多相流模型在歐拉坐標(biāo)系下,通過(guò)求解基本控制方程及每個(gè)網(wǎng)格內(nèi)的流體體積分?jǐn)?shù)來(lái)進(jìn)行相間流體的表面跟蹤[13]。

    每個(gè)單元的流體體積分?jǐn)?shù)F[14]:

    式中,為相,取1或2;F=0為單元內(nèi)無(wú)該相流體;0<F<1為單元內(nèi)既存在該相流體又存在其他相流體;F=1為單元內(nèi)充滿該相流體。

    F滿足如下輸運(yùn)方程:

    式中,、分別為流體在和方向的流動(dòng)速度,m/s。

    2 淺海天然氣管道泄漏擴(kuò)散仿真模型

    2.1 網(wǎng)格模型

    以國(guó)外某天然氣管道海峽穿越段為研究對(duì)象,該輸氣管道管徑為1 016 mm,輸送壓力為9.8 MPa。為了簡(jiǎn)化數(shù)值模型,取實(shí)際工程穿越段的最大水深22 m作為管道數(shù)值模型的水深,忽略海底地形因素的影響。受水體運(yùn)動(dòng)的影響,泄漏氣體在沿海流方向具有水平最大擴(kuò)散距離,因此選取沿海流方向的垂直切面為研究平面。根據(jù)預(yù)先模擬結(jié)果,在22 m水深條件下泄漏氣體上浮到水面的過(guò)程中,沿海流方向的最大擴(kuò)散距離不超過(guò)30 m,最終選擇長(zhǎng)40 m、寬22 m的矩形建立二維幾何模型,泄漏孔位于模型底部距離左側(cè)10 m處(見(jiàn)圖1)。采用結(jié)構(gòu)化網(wǎng)格劃分方法,相比于非結(jié)構(gòu)化網(wǎng)格,整個(gè)數(shù)值模型的網(wǎng)格數(shù)量更少,更有利于計(jì)算結(jié)果的收斂[15]??紤]到泄漏口處氣體與水體的相互摻混屬于劇烈、復(fù)雜的湍流運(yùn)動(dòng),對(duì)泄漏口及氣泡擴(kuò)散區(qū)域的網(wǎng)格進(jìn)一步細(xì)化加密以防止計(jì)算結(jié)果發(fā)散。

    圖1 二維網(wǎng)格模型

    2.2 邊界條件設(shè)置

    計(jì)算域底部的氣體泄漏口設(shè)置為速度入口,底部其余邊界設(shè)置為無(wú)滑移的壁面。頂部邊界作為水的自由表面,設(shè)置為自由出流邊界。在靜水條件時(shí),計(jì)算域左右兩側(cè)邊界均設(shè)置為自由出流;考慮水流速度對(duì)氣體擴(kuò)散的影響,計(jì)算域左側(cè)采用速度入口,通過(guò)UDF函數(shù)定義水流速度隨深度的變化規(guī)律,計(jì)算域右側(cè)仍采用自由出流邊界。

    2.3 條件變量設(shè)置

    以國(guó)外某天然氣管道海峽穿越段重現(xiàn)期1年的海流速度作為實(shí)際水流速度(見(jiàn)表1),擬合出水流速度隨水深變化的公式(見(jiàn)式(7)):

    式中,v為距海底高度處的海流流速,m/s;為距海底的高度,m。

    以靜水條件和實(shí)際水面流速的2倍(2.72 m/s)設(shè)置對(duì)照組。2倍水面流速情況下不同深度的水流速度按式(8)求得[20]:

    式中,v為海面流速,m/s;為水深,m。

    表1 不同深度處的海流速度

    根據(jù)不同的泄漏孔徑、泄漏速率和水流速度,設(shè)置多組模擬工況進(jìn)行計(jì)算,詳細(xì)條件變量設(shè)置如表2所示。

    表2 條件變量設(shè)置

    2.4 求解器設(shè)置

    選用非穩(wěn)態(tài)壓力基求解器,采用雙精度模式提高計(jì)算結(jié)果的精確性。操作條件中勾選重力選項(xiàng),考慮重力對(duì)流體運(yùn)動(dòng)的影響,選取VOF多相流模型和標(biāo)準(zhǔn)湍流模型作為基本計(jì)算模型,運(yùn)用VOF模型中的隱式體積力公式,提高迭代計(jì)算的穩(wěn)定性[21]。選擇PISO算法作為求解器的壓力?速度耦合方式,相比其他類(lèi)型的非穩(wěn)態(tài)求解算法,PISO的殘差收斂更穩(wěn)定,求解速度更快[22?24]。

    3 數(shù)值模擬結(jié)果及分析

    3.1 水下擴(kuò)散過(guò)程分析

    對(duì)國(guó)外某天然氣管道海峽穿越段的實(shí)際情況進(jìn)行探討,即泄漏孔徑60 mm,泄漏速率260 m/s,水流速度為實(shí)際海流速度。

    不同時(shí)刻氣液兩相體積分?jǐn)?shù)云圖見(jiàn)圖2。由圖2可以看出,在泄漏初始階段,天然氣在管道內(nèi)壓的作用下從泄漏口噴射進(jìn)海水中,高速氣體因海水的阻礙在泄漏口上方急劇減慢,被海水包裹形成氣團(tuán)。氣團(tuán)中心部分天然氣上升較快,兩側(cè)因回流作用上升較慢,導(dǎo)致氣團(tuán)呈現(xiàn)蘑菇狀(見(jiàn)圖2(a)),并在海流的作用下開(kāi)始往右偏移。隨著泄漏的進(jìn)行,后續(xù)的泄漏氣體不斷地對(duì)上方的氣團(tuán)進(jìn)行動(dòng)能和體積的補(bǔ)充,在浮力的作用和下方氣體的推動(dòng)下,氣團(tuán)繼續(xù)向上運(yùn)動(dòng),體積也緩慢變大,強(qiáng)烈的湍流作用將蘑菇狀氣團(tuán)逐步分裂為多個(gè)塊狀氣團(tuán)(見(jiàn)圖2(b))。海水壓力的不斷降低使氣團(tuán)不再容易聚集,在氣液兩相之間不斷地摩擦和摻混下,塊狀氣團(tuán)分裂得越來(lái)越小,兩相間的相互滲透越來(lái)越強(qiáng)。氣團(tuán)的持續(xù)向上運(yùn)動(dòng)引起兩側(cè)海水的回流,使兩側(cè)的外部氣團(tuán)在運(yùn)動(dòng)過(guò)程中伴隨著一定程度的旋轉(zhuǎn),氣體的擴(kuò)散范圍越來(lái)越大(見(jiàn)圖3)。此外,越往上運(yùn)移海流的速度越來(lái)越大,海流沖刷作用的增強(qiáng)使大氣泡進(jìn)一步分散成許多小氣泡,并進(jìn)行大范圍擴(kuò)散(見(jiàn)圖2(c))。泄漏發(fā)生5.5 s后,氣體上浮至水面,此時(shí)的水平最大擴(kuò)散距離約為20 m(見(jiàn)圖2(d))。經(jīng)對(duì)比可知,整個(gè)擴(kuò)散過(guò)程的氣液兩相分布與氣泡羽流模型基本一致。

    圖2 不同時(shí)刻氣液兩相體積分?jǐn)?shù)云圖

    圖3 4.6 s時(shí)泄漏氣體左右兩側(cè)水流流線

    3.2 泄漏孔徑的影響

    選取氣體泄漏速率為260 m/s、水流速度為實(shí)際海流速度,分別對(duì)10、30、60、80、100 mm孔徑下的氣體擴(kuò)散過(guò)程進(jìn)行模擬,結(jié)果見(jiàn)圖4。不同泄漏孔徑下氣泡上升高度隨時(shí)間變化的曲線見(jiàn)圖5。

    圖4 泄漏孔徑不同時(shí)氣體擴(kuò)散至海面的氣液兩相體積分?jǐn)?shù)云圖

    圖5 不同泄漏孔徑下氣泡上升高度隨時(shí)間的變化

    由圖4、5可知,泄漏孔徑為10、30 mm時(shí),泄漏口處氣體以少量小氣團(tuán)和較多小氣泡的形式存在,小氣團(tuán)在向上擴(kuò)散的過(guò)程中又分裂為許多小氣泡。因泄漏量較少,下部氣團(tuán)對(duì)上方氣團(tuán)的補(bǔ)充不足,導(dǎo)致氣團(tuán)的體積小,從而受到的浮力小,上移速度較慢,最終分別于7.1、6.0 s時(shí)擴(kuò)散至海面;當(dāng)泄漏孔徑增大至60、80、100 mm時(shí),由于泄漏量變大,氣體主要以大氣泡的形式存在,上浮過(guò)程中大氣泡與海水摻混逐步分裂出大量的小氣泡。因氣團(tuán)體積較大,引起的回流作用較為強(qiáng)烈,泄漏氣體的左右兩側(cè)氣泡出現(xiàn)了明顯的旋轉(zhuǎn)現(xiàn)象,上部氣團(tuán)不斷受到下部氣團(tuán)動(dòng)能和體積的補(bǔ)充,導(dǎo)致氣體上移速度較快,分別于5.5、4.9、4.1 s時(shí)擴(kuò)散至海面。由圖5還可以看出,隨著泄漏孔徑的增大,氣泡上升得越快,到達(dá)水面的時(shí)間越短。

    3.3 泄漏速率的影響

    以實(shí)際海流速度建立流場(chǎng),取60 mm泄漏孔徑,經(jīng)計(jì)算該管道天然氣的泄漏速率約為260 m/s,并以其0.25、0.50、0.75、1.25倍,即65、130、195、325 m/s設(shè)置對(duì)照組,對(duì)5種泄漏速率下的氣體擴(kuò)散過(guò)程進(jìn)行研究。不同泄漏速率下氣體擴(kuò)散至海面時(shí)氣液兩相體積分?jǐn)?shù)云圖見(jiàn)圖6,氣泡上升高度隨時(shí)間的變化曲線見(jiàn)圖7。

    圖6 不同泄漏速率下氣體擴(kuò)散至海面時(shí)氣液兩相體積分?jǐn)?shù)云圖

    由圖6、7可知,泄漏發(fā)生1.0 s時(shí),相同孔徑下氣體泄漏速率越大,泄漏口上方的初始?xì)鈭F(tuán)體積越大,5種泄漏速率下的氣團(tuán)直徑分別為2.3、3.0、4.0、4.8、5.0 m。在水下擴(kuò)散的過(guò)程中,泄漏速率為65、130 m/s的氣體主要以小氣團(tuán)的形式存在,而泄漏速率為195、260、325 m/s的氣體均以塊狀氣團(tuán)和大氣泡的形式存在。泄漏速率越大,氣團(tuán)初始動(dòng)能越大,泄漏過(guò)程中上部氣團(tuán)受到下部氣團(tuán)動(dòng)能的補(bǔ)充越大;泄漏速率越大,氣體泄漏量越大,導(dǎo)致水中氣團(tuán)體積越大,受到浮力的作用也越大。同時(shí),動(dòng)能和體積兩方面的影響使氣團(tuán)引起的回流作用更為強(qiáng)烈,泄漏氣體的左右兩側(cè)出現(xiàn)了明顯的旋轉(zhuǎn)現(xiàn)象。隨著泄漏速率的增大,氣泡上升得越快,到達(dá)水面的時(shí)間越短,5種泄漏速率下的氣團(tuán)分別于7.0、5.8、5.7、5.5、4.5 s上升至海面。

    圖7 泄漏速率不同時(shí)氣泡上升高度隨時(shí)間的變化

    3.4 水流速度的影響

    由于受不同季節(jié)海面風(fēng)向、洋流流向的變化以及海底地形的起伏等諸多因素的影響,實(shí)際的海洋流場(chǎng)較為復(fù)雜,只能對(duì)其簡(jiǎn)化進(jìn)行模擬。以60 mm泄漏孔徑和260 m/s泄漏速率為基準(zhǔn),選取靜水條件、實(shí)際海流速度以及2倍水面流速進(jìn)行計(jì)算,對(duì)得到的海流速度場(chǎng)進(jìn)行模擬。圖8為不同水流速度下氣體擴(kuò)散至海面時(shí)氣液兩相體積分?jǐn)?shù)云圖。

    圖8 不同水流速度下氣體擴(kuò)散至海面時(shí)氣液兩相體積分?jǐn)?shù)云圖

    由圖8可知,由于靜水條件下的氣體在水中只受浮力的作用,水平方向的擴(kuò)散分布較為對(duì)稱(chēng),最終以5.2 s上升至海平面,此時(shí)的水平最大擴(kuò)散距離約為22 m,單側(cè)氣體距泄漏口的最大水平距離約為10 m。在水流速度的影響下,氣體不僅受到浮力的作用,還受到左側(cè)水流的作用。由于實(shí)際工況下的水流速度較小,氣體所受浮力的影響大于左側(cè)水流的影響,氣體主要向上部空間擴(kuò)散,擴(kuò)散軌跡與海底的夾角較大。上部氣體最終以5.5 s上升至海平面,此時(shí)的水平最大擴(kuò)散距離為20 m,最右側(cè)氣體水平運(yùn)動(dòng)距離約為16 m。在水面流速為2倍的海流流場(chǎng)下,左側(cè)水流對(duì)氣體的推動(dòng)作用較強(qiáng),主要向斜上方進(jìn)行擴(kuò)散,擴(kuò)散軌跡與海底的夾角較小。氣體最終以5.2 s上升至海平面,此時(shí)的水平最大擴(kuò)散距離約為23 m,最右側(cè)氣體距泄漏口的水平最大距離約為22 m。根據(jù)3種水流速度下的模擬結(jié)果發(fā)現(xiàn),水流速度越大,氣體受水流的推動(dòng)作用越明顯,氣體運(yùn)動(dòng)軌跡與海底的夾角越小,沿海流方向的水平運(yùn)移距離越大,而氣體上升到水面的時(shí)間與水流速度并無(wú)明顯關(guān)系。

    4 結(jié) 論

    針對(duì)國(guó)外某天然氣管道海峽穿越段,根據(jù)計(jì)算流體力學(xué)與多相流動(dòng)理論,建立淺海海底管道泄漏擴(kuò)散過(guò)程的計(jì)算模型。以泄漏孔徑、泄漏速率和水流速度為條件變量,運(yùn)用Fluent軟件進(jìn)行有限元數(shù)值模擬,得出以下結(jié)論:

    (1)分析不同時(shí)刻的氣液兩相體積分?jǐn)?shù)云圖,發(fā)現(xiàn)水下氣體擴(kuò)散分為三個(gè)階段,即泄漏口上方形成氣團(tuán)、氣團(tuán)呈蘑菇狀上升、氣團(tuán)由大氣泡分裂為小氣泡,模擬結(jié)果與氣泡羽流模型的描述基本一致。

    (2)分析各個(gè)對(duì)照組的模擬結(jié)果,發(fā)現(xiàn)泄漏孔徑和泄漏速率對(duì)水下氣體擴(kuò)散過(guò)程具有顯著影響。泄漏孔徑與泄漏速率越大,氣體泄漏量越大,上浮過(guò)程中的氣團(tuán)體積越大;氣團(tuán)體積越大,受浮力的作用越強(qiáng),向上運(yùn)動(dòng)得越快,泄漏孔徑為100 mm時(shí)氣體上浮至水面僅需4.1 s。

    (3)研究不同水流速度下的水下泄漏氣體擴(kuò)散軌跡,發(fā)現(xiàn)水流速度對(duì)氣泡運(yùn)動(dòng)的水平最大運(yùn)移距離起決定性作用。水流速度越大,氣體受水流的推動(dòng)作用越明顯,氣體運(yùn)動(dòng)軌跡與海底的夾角越小,沿海流方向的水平運(yùn)移距離越大。2倍水面流速條件下擴(kuò)散至海面的最右側(cè)氣體距泄漏口的水平最大距離可達(dá)22 m,擴(kuò)散距離最遠(yuǎn)。

    [1] 趙冬巖,余建星.海底管道完整性管理研究[J].海洋技術(shù)學(xué)報(bào),2008,27(4):71?74.

    Zhao D Y,Yu J X.Study on submarine pipeline integrity management[J].Ocean Technology,2008,27(4):71?74.

    [2] 田輝.天然氣管道泄漏爆炸動(dòng)力效應(yīng)研究及危險(xiǎn)區(qū)域的劃分[D].蕪湖:安徽理工大學(xué),2012.

    [3] Trevor J M.Bubble plumes in stratified environments[J].Journal of Fluid Mechanics,1978,85(4):655?672.

    [4] Milgram J H.Mean flow in round bubble plumes[J].Journal of Fluid Mechanics,1983,133(1):345.

    [5] Riedl M J,F(xiàn)annel?p T K.Bubble plumes and their interaction with the water surface[J].Applied Ocean Research,2000,22(2):119?128.

    [6] 景海泳,余建星,杜尊峰,等.海底管道水下氣體擴(kuò)散FLUENT仿真分析[J].海洋技術(shù),2012,31(3):82?85.

    Jing H Y,Yu J X,Du Z F,et al.Analysis of underwater gas diffusion from piping based on FLUENT[J].Ocean Technology,2012,31(3):82?85.

    [7] 文闖,延斌,王憲全,等.海底管線天然氣泄漏過(guò)程數(shù)值模擬[J].常州大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,27(2):72?77.

    Wen C, Yan B,Wang X Q,et al.Numerical simulation of natural gas leakage in subsea pipelines[J].Journal of Changzhou University(Natural Science Edition),2015,27(2):72?77.

    [8] 竇梓元.水體中天然氣管道泄漏擴(kuò)散規(guī)律及后果分析[D].成都:西南石油大學(xué),2015.

    [9] 李新宏,陳國(guó)明,朱紅衛(wèi),等.海底輸氣管道泄漏天然氣擴(kuò)散風(fēng)險(xiǎn)研究[J].石油科學(xué)通報(bào),2016,1(3):390?400.

    Li X H,Chen G M,Zhu H W,et al.Research into the risk of natural gas spread from submarine natural gas pipeline leakage[J].Petroleum Science Bulletin,2016,1(3):390?400.

    [10] 王少雄,李玉星,劉翠偉,等.水下輸氣管道泄漏擴(kuò)散特性模擬研究[J].化工學(xué)報(bào),2020,71(4):496?509.

    Wang S X,Li Y X,Liu C W,et al.Numerical simulation on leakage and diffusion characteristics of underwater gas pipeline[J].化工學(xué)報(bào),2020,71(4):496?509.

    [11] 王瑞金,張凱,王剛.Fluent技術(shù)基礎(chǔ)與應(yīng)用實(shí)例[M].北京:清華大學(xué)出版社,2007.

    [12] 熊莉芳,林源,李世武.湍流模型及其在FLUENT軟件中的應(yīng)用[J].工業(yè)加熱,2007,36(4):13?15.

    Xiong L F,Lin Y,Li S W.turbulent model and its application to the FLUENT[J].Industrial Heating,2007,36(4):13?15.

    [13] 張健,方杰,范波芹.VOF方法理論與應(yīng)用綜述[J].水利水電科技進(jìn)展,2005,25(2):67?70.

    Zhang J,F(xiàn)ang J,F(xiàn)an B Q.Advances in research of VOF method[J].Advances in Science and Technology of Water Resources,2005,25(2):67?70.

    [14] 韓金珂.水下輸氣管道泄漏擴(kuò)散特性研究[D].青島:中國(guó)石油大學(xué)(華東),2018.

    [15] 王明強(qiáng),朱永梅,劉文欣.有限元網(wǎng)格劃分方法應(yīng)用研究[J].機(jī)械設(shè)計(jì)與制造,2004(1):22?24.

    Wang M Q,Zhu Y M,Liu W X.The research on finite element mesh generation method[J].Machinery Design & Manufacture,2004(1):22?24.

    [16] Sun Y,Cao X,Liang F,et al.Investigation on underwater gas leakage and dispersion behaviors based on coupled Eulerian?Lagrangian CFD model?ScienceDirect[J].Process Safety and Environmental Protection,2020,136:268?279.

    [17] Zhang Y,Zhu J,Peng Y,et al.Experimental research of flow rate and diffusion behavior of nature gas leakage underwater[J].Journal of Loss Prevention in the Process Industries,2020,65:104?119.

    [18] 霍春勇,董玉華,余大濤,等.長(zhǎng)輸管線氣體泄漏率的計(jì)算方法研究[J].石油學(xué)報(bào),2004,25(1):101?105.

    Huo C Y,Dong Y H,Yu D T,et al.Estimation of accidental gas release flow rate in long transmission pipelines[J].Acta Petrolei Sinica,2004,25(1):101?105.

    [19] 王大慶,霍春勇,高惠臨.長(zhǎng)輸管線氣體泄漏率簡(jiǎn)化計(jì)算方法[J].天然氣工業(yè),2008,28(1):116?118.

    Wang D Q,Huo C Y,Gao H L.Simplified method for calculation long?distance pipeline leakage rate[J].Natural Gas Industry,2008,28(1):116?118.

    [20] 徐興平.海洋石油工程概論[M].東營(yíng):中國(guó)石油大學(xué)出版社,2007.

    [21] 朱約鈞.FLUENT 15.0流場(chǎng)分析實(shí)戰(zhàn)指南[M].北京:人民郵電出版社,2015.

    [22] 王彤,谷傳綱,楊波,等.非定常流動(dòng)計(jì)算的PISO算法[J].水動(dòng)力學(xué)研究與進(jìn)展,2003,18(2):233?239.

    Wang T,Gu C G,Yang B,et al.PISO algorithm for unsteady flow field[J].Chinese Journal of Hydrodynamics,2003,18(2):233?239.

    [23] 李遠(yuǎn)朋,范潮海,張茹.油田集輸管道體系區(qū)域風(fēng)險(xiǎn)評(píng)價(jià)方法[J].東北石油大學(xué)學(xué)報(bào),2019,43(6):118?124.

    Li Y P,F(xiàn)an C H,Zhang R,et al.Systematic regional risk assessment method for oil field gathering and transportation pipeline[J].Journal of Northeast Petroleum University,2019,43(6):118?124.

    [24] 何國(guó)璽,梁永圖,李巖松,等.埋地液烴管道泄漏擴(kuò)散及泄漏量測(cè)算研究進(jìn)展[J].油氣儲(chǔ)運(yùn),2017,36(1):8?20.

    He G X,Liang Y T,Li Y S,et al.Research progress of leakage diffusion and leakage volume calculation for buried liquid hydrocarbon pipeline[J].Oil & Gas Storage and Transportation,2017,36(1):8?20.

    Simulation Study of Leakage and Diffusion of Shallow Subsea Natural Gas Pipeline

    Li Yanhao1,2, Hou Lei1,2, Yu Qiaoyan1,2, Chai Chong1,2, Xiao Kaixi1,2

    (1. College of Mechanical and Transportation Engineering, China University of Petroleum(Beijing),Beijing 102249,China;2. National Engineering Laboratory for Pipeline Safety, China University of Petroleum(Beijing),Beijing 102249,China)

    To evaluate the consequences of leakage accidents of the shallow subsea natural gas pipeline, this paper builds a calculation model of the shallow subsea pipeline leakage and diffusion process in light of computational fluid dynamics and multiphase flow theories for the strait?crossing section of a natural gas pipeline abroad. Three main factors, leakage aperture, leakage rate, and water flow velocity, are selected as condition variables to simulate the motion process of the gas?liquid two?phase flow under different conditions. The results show that underwater gas diffusion can be divided into three stages: the formation of gas masses above the leakage hole, the mushroom?like rise of gas masses, and the splitting of gas masses from large bubbles into small ones. Leakage aperture and leakage rate have a significant effect on the time when the underwater gas diffuses to the water surface. A larger leakage aperture and a higher leakage rate lead to a larger gas leakage amount, which further results in larger volumes of underwater air masses and ultimately a shorter time for them to reach the water surface. Water flow velocity has a significant effect on gas diffusion trajectory. As water flow velocity increases, the angle between the gas trajectory and the seabed decreases, and the diffusion distance along the ocean current direction becomes longer. This study can provide scientific guidance for emergency treatment of underwater natural gas pipeline leakage accidents.

    Submarine pipeline; Gas leakage; Multiphase flow; Computational fluid dynamics

    TE88

    A

    10.3969/j.issn.1006?396X.2022.02.012

    1006?396X(2022)02?0074?07

    2021?03?08

    2021?09?18

    中國(guó)石油天然氣集團(tuán)有限公司——中國(guó)石油大學(xué)(北京)戰(zhàn)略合作科技專(zhuān)項(xiàng)(ZLZX2020?05)。

    李延豪(1996?),男,碩士研究生,從事管道完整性技術(shù)方面研究;E?mail:cup_lyh@163.com。

    侯磊(1966?),男,博士,教授,博士生導(dǎo)師,從事油氣管道輸送技術(shù)方面研究;E?mail:houleicup@126.com。

    (編輯 王戩麗)

    猜你喜歡
    氣團(tuán)水流氣泡
    檸檬氣泡水
    欣漾(2024年2期)2024-04-27 15:19:49
    哪股水流噴得更遠(yuǎn)
    能俘獲光的水流
    輸水管線啟動(dòng)填充過(guò)程含滯留氣團(tuán)瞬變流數(shù)值模擬
    SIAU詩(shī)杭便攜式氣泡水杯
    新潮電子(2021年7期)2021-08-14 15:53:12
    浮法玻璃氣泡的預(yù)防和控制對(duì)策
    基于后向軌跡模式的合肥市臭氧來(lái)源分析
    我只知身在水中,不覺(jué)水流
    文苑(2020年6期)2020-06-22 08:41:56
    氣團(tuán)來(lái)源對(duì)沿海城市PM2.5中二次水溶性無(wú)機(jī)離子形成特征的影響研究
    HZ330CL橫折缺陷原因分析與改進(jìn)
    欧美性猛交╳xxx乱大交人| 免费观看精品视频网站| 欧美日韩中文字幕国产精品一区二区三区| 欧美精品啪啪一区二区三区| 长腿黑丝高跟| 一边摸一边抽搐一进一小说| 人妻夜夜爽99麻豆av| 国产亚洲欧美98| 国产乱人视频| 蜜桃久久精品国产亚洲av| 精品久久久久久久久av| 欧美最黄视频在线播放免费| 免费大片18禁| 日本精品一区二区三区蜜桃| 亚洲人与动物交配视频| 久久久久亚洲av毛片大全| 非洲黑人性xxxx精品又粗又长| 午夜a级毛片| 亚洲av成人精品一区久久| 一个人观看的视频www高清免费观看| 99国产极品粉嫩在线观看| 别揉我奶头 嗯啊视频| 欧美bdsm另类| 很黄的视频免费| 亚洲国产精品sss在线观看| 日韩大尺度精品在线看网址| 观看美女的网站| 欧美三级亚洲精品| 色视频www国产| 91av网一区二区| 久久精品国产99精品国产亚洲性色| 久久久久久久久久黄片| 久久久久久久久久成人| 免费观看的影片在线观看| 日本五十路高清| av黄色大香蕉| 一区福利在线观看| 日本三级黄在线观看| 变态另类丝袜制服| 美女被艹到高潮喷水动态| a级毛片免费高清观看在线播放| 成年女人毛片免费观看观看9| 亚洲内射少妇av| 久久人妻av系列| 丰满人妻一区二区三区视频av| 尤物成人国产欧美一区二区三区| 一二三四社区在线视频社区8| 日日摸夜夜添夜夜添小说| 国产不卡一卡二| 国产精品电影一区二区三区| 18禁在线播放成人免费| 亚洲色图av天堂| 可以在线观看的亚洲视频| 91麻豆精品激情在线观看国产| a级一级毛片免费在线观看| 深夜精品福利| 激情在线观看视频在线高清| 51国产日韩欧美| 欧美黑人欧美精品刺激| 中文字幕高清在线视频| 亚洲五月婷婷丁香| 男女那种视频在线观看| 男女视频在线观看网站免费| 午夜福利18| 91久久精品电影网| 高清毛片免费观看视频网站| av国产免费在线观看| 成人精品一区二区免费| 男插女下体视频免费在线播放| 成人av一区二区三区在线看| xxxwww97欧美| 国产一区二区在线av高清观看| 国产免费一级a男人的天堂| 蜜桃久久精品国产亚洲av| 色综合亚洲欧美另类图片| 黄色一级大片看看| 99在线人妻在线中文字幕| 少妇的逼水好多| 性色avwww在线观看| 中文字幕人成人乱码亚洲影| 日本撒尿小便嘘嘘汇集6| 精品午夜福利在线看| 又爽又黄a免费视频| 亚洲 国产 在线| aaaaa片日本免费| 欧美黄色片欧美黄色片| 又黄又爽又免费观看的视频| 欧美成狂野欧美在线观看| 国产精品免费一区二区三区在线| 国产精品亚洲美女久久久| 非洲黑人性xxxx精品又粗又长| 美女黄网站色视频| 国产激情偷乱视频一区二区| 日韩亚洲欧美综合| 亚洲熟妇熟女久久| 全区人妻精品视频| 一本精品99久久精品77| 日日干狠狠操夜夜爽| 亚洲va日本ⅴa欧美va伊人久久| 亚洲欧美精品综合久久99| 欧美激情在线99| 日韩欧美 国产精品| 亚洲av中文字字幕乱码综合| 精品人妻熟女av久视频| 少妇的逼好多水| 每晚都被弄得嗷嗷叫到高潮| 日本三级黄在线观看| 亚洲欧美清纯卡通| 搡女人真爽免费视频火全软件 | 国产探花极品一区二区| 午夜激情欧美在线| 麻豆国产97在线/欧美| 大型黄色视频在线免费观看| 国产在视频线在精品| 亚洲人成电影免费在线| 国产成年人精品一区二区| 网址你懂的国产日韩在线| 久久九九热精品免费| 在线观看av片永久免费下载| 国产麻豆成人av免费视频| 听说在线观看完整版免费高清| 久久欧美精品欧美久久欧美| 成人精品一区二区免费| 51国产日韩欧美| 欧美精品国产亚洲| 国产黄片美女视频| 国产 一区 欧美 日韩| a在线观看视频网站| 国内揄拍国产精品人妻在线| 国产乱人伦免费视频| 91麻豆av在线| 国产亚洲欧美在线一区二区| 91字幕亚洲| 99久久99久久久精品蜜桃| 俺也久久电影网| 五月玫瑰六月丁香| 亚洲av成人不卡在线观看播放网| 在线国产一区二区在线| 在线看三级毛片| 51午夜福利影视在线观看| 精品国产三级普通话版| 全区人妻精品视频| 精品福利观看| 午夜日韩欧美国产| 亚洲在线自拍视频| 亚洲avbb在线观看| 热99re8久久精品国产| 久久久久九九精品影院| 变态另类丝袜制服| 老熟妇乱子伦视频在线观看| 亚洲av不卡在线观看| 日韩欧美精品免费久久 | 日本五十路高清| 内射极品少妇av片p| 丰满人妻熟妇乱又伦精品不卡| 1024手机看黄色片| 黄色视频,在线免费观看| 婷婷亚洲欧美| 欧美丝袜亚洲另类 | 国产三级中文精品| 91狼人影院| 久久久精品大字幕| 久久久久性生活片| 搡老妇女老女人老熟妇| 一进一出好大好爽视频| av中文乱码字幕在线| 白带黄色成豆腐渣| 99久久精品一区二区三区| 成人亚洲精品av一区二区| 日韩欧美精品免费久久 | 亚洲综合色惰| 亚洲av二区三区四区| 成人一区二区视频在线观看| 午夜福利成人在线免费观看| 成年人黄色毛片网站| 免费电影在线观看免费观看| 久久久久免费精品人妻一区二区| 亚洲人成网站在线播| 欧美成狂野欧美在线观看| 少妇的逼水好多| 亚洲精品粉嫩美女一区| 97碰自拍视频| 国产伦一二天堂av在线观看| 老司机福利观看| 亚洲欧美清纯卡通| 偷拍熟女少妇极品色| av欧美777| 国产亚洲精品久久久久久毛片| 99久久精品一区二区三区| 午夜福利成人在线免费观看| 好男人在线观看高清免费视频| 国产一区二区激情短视频| 亚洲国产精品sss在线观看| 小说图片视频综合网站| 欧美另类亚洲清纯唯美| 亚洲av日韩精品久久久久久密| 亚洲真实伦在线观看| 一个人看视频在线观看www免费| 午夜福利成人在线免费观看| 狠狠狠狠99中文字幕| 热99在线观看视频| h日本视频在线播放| 欧美国产日韩亚洲一区| 国产真实乱freesex| 色噜噜av男人的天堂激情| 欧美激情久久久久久爽电影| 夜夜夜夜夜久久久久| or卡值多少钱| 国内揄拍国产精品人妻在线| 亚洲精品一卡2卡三卡4卡5卡| 在线天堂最新版资源| 亚洲不卡免费看| 国产成人aa在线观看| 国产欧美日韩一区二区精品| 18禁黄网站禁片午夜丰满| 精品一区二区三区视频在线| 欧美色视频一区免费| 久久久国产成人免费| 国产成人aa在线观看| 麻豆成人午夜福利视频| 美女cb高潮喷水在线观看| 国产成人影院久久av| 久久性视频一级片| 久久亚洲真实| а√天堂www在线а√下载| 精品免费久久久久久久清纯| 在线观看美女被高潮喷水网站 | 亚洲精品日韩av片在线观看| 精品久久久久久,| 色噜噜av男人的天堂激情| 午夜福利欧美成人| 99久久精品热视频| 亚洲激情在线av| www.色视频.com| 精品人妻偷拍中文字幕| 99久久久亚洲精品蜜臀av| 国产高潮美女av| 国产三级中文精品| 欧美日本亚洲视频在线播放| 亚洲国产精品sss在线观看| 97人妻精品一区二区三区麻豆| 免费在线观看影片大全网站| 国产亚洲精品av在线| 日韩欧美三级三区| 午夜免费激情av| 99久久精品一区二区三区| 深夜精品福利| 午夜免费成人在线视频| 亚洲国产精品sss在线观看| 欧美bdsm另类| 亚洲自偷自拍三级| 国产精品自产拍在线观看55亚洲| 国产一区二区三区在线臀色熟女| av在线天堂中文字幕| 欧美在线一区亚洲| 99riav亚洲国产免费| 国产真实乱freesex| 久久久久性生活片| 日本黄色片子视频| 日韩有码中文字幕| 欧美黄色淫秽网站| 中出人妻视频一区二区| 日本三级黄在线观看| 高清在线国产一区| 国产色婷婷99| 在线播放无遮挡| 亚洲国产精品合色在线| 婷婷精品国产亚洲av| 日本黄大片高清| 成人特级黄色片久久久久久久| 黄色一级大片看看| 国产一级毛片七仙女欲春2| 人妻夜夜爽99麻豆av| 成人国产一区最新在线观看| 亚洲精品一区av在线观看| 久久精品久久久久久噜噜老黄 | 精品乱码久久久久久99久播| 欧美黄色片欧美黄色片| 亚洲精品色激情综合| 乱人视频在线观看| 欧美成人性av电影在线观看| 久久亚洲真实| 久久九九热精品免费| 三级毛片av免费| 婷婷丁香在线五月| 色尼玛亚洲综合影院| 亚洲在线观看片| a级一级毛片免费在线观看| 乱码一卡2卡4卡精品| 校园春色视频在线观看| 中文在线观看免费www的网站| 别揉我奶头 嗯啊视频| 国产69精品久久久久777片| 久99久视频精品免费| 午夜精品久久久久久毛片777| 亚洲欧美激情综合另类| 中文字幕免费在线视频6| 精品久久久久久久人妻蜜臀av| 免费av观看视频| 波野结衣二区三区在线| 亚洲专区中文字幕在线| 国产高清激情床上av| 露出奶头的视频| 日本 欧美在线| 无人区码免费观看不卡| 欧美黄色淫秽网站| 亚洲精品乱码久久久v下载方式| 中文字幕av在线有码专区| 看十八女毛片水多多多| 又粗又爽又猛毛片免费看| 一a级毛片在线观看| 亚洲无线在线观看| 桃红色精品国产亚洲av| 国产美女午夜福利| 亚洲天堂国产精品一区在线| 亚洲一区二区三区不卡视频| 国产午夜福利久久久久久| 丁香欧美五月| 丰满乱子伦码专区| 亚洲五月婷婷丁香| 成年免费大片在线观看| 天堂影院成人在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲成人免费电影在线观看| 日韩免费av在线播放| 久久精品久久久久久噜噜老黄 | 性插视频无遮挡在线免费观看| 久久久国产成人精品二区| 成人精品一区二区免费| 97碰自拍视频| 麻豆国产97在线/欧美| 在线观看美女被高潮喷水网站 | 精品久久久久久成人av| 国产亚洲精品av在线| 国产精品久久久久久亚洲av鲁大| 国产成年人精品一区二区| 午夜老司机福利剧场| 欧美日韩亚洲国产一区二区在线观看| 国产av麻豆久久久久久久| 国产一级毛片七仙女欲春2| 久久99热这里只有精品18| 性插视频无遮挡在线免费观看| 特大巨黑吊av在线直播| 欧美日本视频| 国产乱人视频| 欧美一区二区国产精品久久精品| 好男人在线观看高清免费视频| 一本精品99久久精品77| 网址你懂的国产日韩在线| eeuss影院久久| 亚洲七黄色美女视频| av女优亚洲男人天堂| 亚洲成a人片在线一区二区| 精品午夜福利视频在线观看一区| 日韩欧美在线二视频| 久久伊人香网站| 成年免费大片在线观看| 如何舔出高潮| 极品教师在线视频| 国产成人啪精品午夜网站| 免费一级毛片在线播放高清视频| 成人国产综合亚洲| 欧美一区二区精品小视频在线| 欧美日本亚洲视频在线播放| 国产精品,欧美在线| 欧美精品国产亚洲| 国产精品爽爽va在线观看网站| 狠狠狠狠99中文字幕| 桃红色精品国产亚洲av| 熟女人妻精品中文字幕| 欧美黄色淫秽网站| 99精品久久久久人妻精品| 亚洲综合色惰| 少妇被粗大猛烈的视频| 久久中文看片网| 国产成年人精品一区二区| 变态另类丝袜制服| 亚洲av熟女| 99久久无色码亚洲精品果冻| 日韩有码中文字幕| 免费看日本二区| 欧美高清性xxxxhd video| 国产视频内射| 欧美另类亚洲清纯唯美| 免费观看精品视频网站| 国产精品久久电影中文字幕| 99久久成人亚洲精品观看| 成人欧美大片| 久久国产精品影院| 亚洲一区二区三区色噜噜| 18禁黄网站禁片免费观看直播| 精品乱码久久久久久99久播| 免费看美女性在线毛片视频| 色播亚洲综合网| 一级作爱视频免费观看| 老司机午夜福利在线观看视频| 国产色婷婷99| 激情在线观看视频在线高清| 1000部很黄的大片| av国产免费在线观看| 最新中文字幕久久久久| 久久久久久久久大av| 亚洲精品粉嫩美女一区| 亚洲精品久久国产高清桃花| 亚洲成a人片在线一区二区| 夜夜躁狠狠躁天天躁| 亚洲最大成人中文| 亚洲成人久久性| 国产男靠女视频免费网站| 99国产极品粉嫩在线观看| 久久精品国产自在天天线| 久久久久久久久久黄片| 亚洲七黄色美女视频| 久久午夜亚洲精品久久| 精品日产1卡2卡| a级一级毛片免费在线观看| 国产乱人伦免费视频| 欧美+亚洲+日韩+国产| 欧美午夜高清在线| 亚洲av熟女| 中文字幕精品亚洲无线码一区| 色尼玛亚洲综合影院| 色5月婷婷丁香| 精品一区二区三区视频在线观看免费| 美女免费视频网站| 一个人免费在线观看的高清视频| 精品久久久久久久末码| 国产视频内射| 啦啦啦观看免费观看视频高清| av天堂中文字幕网| 日韩人妻高清精品专区| 亚洲熟妇熟女久久| 日韩高清综合在线| 蜜桃亚洲精品一区二区三区| 99精品在免费线老司机午夜| 午夜视频国产福利| 欧美区成人在线视频| 久久欧美精品欧美久久欧美| 99久久精品热视频| 最近视频中文字幕2019在线8| 欧美成人性av电影在线观看| 国产真实乱freesex| 国内久久婷婷六月综合欲色啪| 我的老师免费观看完整版| 深夜精品福利| 天美传媒精品一区二区| 99riav亚洲国产免费| 国产69精品久久久久777片| 国产精品野战在线观看| 国产老妇女一区| 精品久久久久久久久亚洲 | 久久国产精品影院| 人妻夜夜爽99麻豆av| 3wmmmm亚洲av在线观看| 床上黄色一级片| 亚洲欧美清纯卡通| 国内久久婷婷六月综合欲色啪| 欧美性猛交╳xxx乱大交人| or卡值多少钱| 精品一区二区三区视频在线| 级片在线观看| av福利片在线观看| 国内揄拍国产精品人妻在线| 欧美日韩乱码在线| 国产视频内射| a级一级毛片免费在线观看| 不卡一级毛片| 日韩欧美三级三区| 91久久精品国产一区二区成人| 久久99热6这里只有精品| 午夜视频国产福利| 国产极品精品免费视频能看的| 国产蜜桃级精品一区二区三区| 亚洲七黄色美女视频| 无人区码免费观看不卡| 99热精品在线国产| 国内精品美女久久久久久| 色综合亚洲欧美另类图片| 在线观看美女被高潮喷水网站 | 俺也久久电影网| 日韩有码中文字幕| 波多野结衣巨乳人妻| 欧美日韩综合久久久久久 | 搡女人真爽免费视频火全软件 | 国产午夜精品久久久久久一区二区三区 | 男女视频在线观看网站免费| 村上凉子中文字幕在线| 亚洲人成网站在线播放欧美日韩| 午夜久久久久精精品| 婷婷丁香在线五月| 午夜久久久久精精品| 亚洲avbb在线观看| 91九色精品人成在线观看| 国产av在哪里看| 久久久久久大精品| 国产成人av教育| 久久久精品大字幕| 一级黄色大片毛片| 亚洲av成人不卡在线观看播放网| 动漫黄色视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚洲精品粉嫩美女一区| 淫秽高清视频在线观看| 在线看三级毛片| 一本精品99久久精品77| 国产午夜精品论理片| 成人国产一区最新在线观看| 在线免费观看的www视频| 成人午夜高清在线视频| 色av中文字幕| 国产伦一二天堂av在线观看| www.色视频.com| 天美传媒精品一区二区| 精品人妻一区二区三区麻豆 | 久久草成人影院| 91狼人影院| 国产蜜桃级精品一区二区三区| 一本久久中文字幕| 精品国内亚洲2022精品成人| 国产伦精品一区二区三区视频9| 亚洲黑人精品在线| 国产精品永久免费网站| 欧美日韩综合久久久久久 | 简卡轻食公司| 啦啦啦韩国在线观看视频| 两个人的视频大全免费| 亚洲av熟女| 日韩欧美国产在线观看| 国产午夜福利久久久久久| 亚洲av美国av| 老司机午夜福利在线观看视频| 看十八女毛片水多多多| 一进一出好大好爽视频| 国产精华一区二区三区| 国产探花极品一区二区| 午夜福利欧美成人| 久久精品91蜜桃| 日韩成人在线观看一区二区三区| 国产黄色小视频在线观看| 亚洲欧美激情综合另类| 伊人久久精品亚洲午夜| 变态另类丝袜制服| 国产麻豆成人av免费视频| 国产精品99久久久久久久久| 少妇人妻一区二区三区视频| 少妇裸体淫交视频免费看高清| 欧美中文日本在线观看视频| 欧美+日韩+精品| 久久精品久久久久久噜噜老黄 | 老司机深夜福利视频在线观看| 欧美一区二区精品小视频在线| 免费高清视频大片| 午夜精品在线福利| 久久久久久久亚洲中文字幕 | 日韩精品中文字幕看吧| 色吧在线观看| 欧美色视频一区免费| 韩国av一区二区三区四区| 精品久久久久久久久亚洲 | 少妇丰满av| 成年女人毛片免费观看观看9| 欧美+亚洲+日韩+国产| 成人av在线播放网站| 国产一区二区三区视频了| 高清毛片免费观看视频网站| 成年版毛片免费区| 色播亚洲综合网| 亚洲一区二区三区不卡视频| 久久久精品欧美日韩精品| 天堂av国产一区二区熟女人妻| 一二三四社区在线视频社区8| 尤物成人国产欧美一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 特大巨黑吊av在线直播| 无遮挡黄片免费观看| 国产视频内射| 久久久久久久久久成人| 99国产综合亚洲精品| 成人亚洲精品av一区二区| 亚洲电影在线观看av| 欧美性猛交黑人性爽| 少妇人妻精品综合一区二区 | 蜜桃亚洲精品一区二区三区| www.色视频.com| 中文字幕精品亚洲无线码一区| 一个人看视频在线观看www免费| 我的女老师完整版在线观看| 熟女人妻精品中文字幕| av天堂中文字幕网| 国产精品爽爽va在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品在线美女| 亚洲av日韩精品久久久久久密| 精品久久久久久久久久免费视频| 国产av在哪里看| 精品久久久久久久末码| 国产在线男女| 日韩欧美一区二区三区在线观看| 精品一区二区三区视频在线| 日本熟妇午夜| 色吧在线观看| 中出人妻视频一区二区| 国产高清视频在线播放一区| 性色av乱码一区二区三区2| 天堂网av新在线| 色噜噜av男人的天堂激情| 国产精品美女特级片免费视频播放器| 成人高潮视频无遮挡免费网站| 国产av不卡久久| 国产精品久久久久久久久免 | 99久久成人亚洲精品观看| 久久久久久久久中文| 久久久久久久久大av| 国产淫片久久久久久久久 | 免费看日本二区|