999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

通球過程皮碗球運動機理研究綜述

2023-03-13 06:51:40王學成張雷雷龔治海張方曉王國振陳俊丞趙磊
遼寧化工 2023年1期
關鍵詞:研究

王學成,張雷雷,龔治海,張方曉,王國振,陳俊丞,趙磊*

通球過程皮碗球運動機理研究綜述

王學成1,張雷雷1,龔治海1,張方曉1,王國振1,陳俊丞2,趙磊2*

(1. 東莞市盛源石油化工有限公司,廣東 東莞 523000; 2. 遼寧石油化工大學,遼寧 撫順 113001)

通球過程利用氮氣推動皮碗球,將管道內(nèi)的汽油或柴油輸送至儲罐,是連接碼頭與廠區(qū)完成油品運輸?shù)挠行侄巍Fね肭蚴峭ㄇ蜻^程中的重要設備,在管道清理操作中也可作為清管器使用。 皮碗球的速度決定著通球過程的效率和安全性,受到后方氣體的推力,管壁的摩擦阻力和前方液體的阻力等因素的共同影響,皮碗球在管道中的運動是一個復雜的運動過程。本文對近10年通球過程皮碗球運動機理研究進行綜述,理論研究部分給出了氣體速度,皮碗球速度和阻力的計算公式,實驗研究部分介紹了旁通率對通球效率的影響,模擬研究部分介紹了通球過程中管道內(nèi)的流場情況。研究結果為相關企業(yè)設計效率高,經(jīng)濟性好的通球清管過程提供參考。同時發(fā)現(xiàn),利用數(shù)值模擬方法解決工業(yè)生產(chǎn)中的復雜問題具有廣闊前景。

通球過程;皮碗球;油氣輸送;旁路清管;清管器

運輸汽油與柴油的船舶需要在碼頭或港口完成卸油操作,卸下的油品通過輸油管道輸送至廠區(qū)儲罐中進行儲存。在油品被注入管道之后,向管道中放置皮碗球設備,皮碗球放置完成后,利用輸氣軟管向管道內(nèi)輸入氮氣,在氮氣的推動下,皮碗球向前移動,進而推動油品沿管道向儲罐方向移動,這一過程通常被稱作通球過程。除此之外,皮碗球設備也可作為清管器使用,當管道內(nèi)存在殘存液體時,可以利用皮碗球?qū)⑵渫瞥龉艿溃辉诠艿莱笗r,也可利用皮碗球推動清洗液,將垢層溶解在清洗液中一起推出管道。皮碗球在油品輸送與管道清管過程中均具有重要作用。受到后方氣體的推力,管壁的摩擦阻力和前方液體的阻力等因素的共同影響,皮碗球在管道中的運動是一個復雜的運動過程。在通球過程中,通入氣體的速度和皮碗球的速度決定著通球過程的效率和安全性,速度過低,油品的輸送時間變長;速度過高,管道內(nèi)的壓力以及流體的沖擊力會對管道與儲罐造成破壞。皮碗球與管壁之間的摩擦阻力又會使皮碗球發(fā)生磨損。為了解決這些可能出現(xiàn)的問題,有必要對通球過程和清管過程中,推動氣體速度與皮碗球速度之間的關系,皮碗球受到的阻力以及影響通球和清管效率的因素進行分析,為相關企業(yè)設計效率高,經(jīng)濟性好的通球清管過程。通過調(diào)研近十年通球過程的文獻發(fā)現(xiàn),研究方法主要有三種,分別是理論研究,實驗研究與模擬研究,成果多集中于近五年,且目前并沒有關于這些研究的系統(tǒng)性綜述。基于此,本文對上述三種研究方法的相關成果進行介紹,包括但不限于氣體速度,皮碗球速度和阻力的計算公式,旁通率對通球效率的影響,通球過程中管道內(nèi)的流場情況等。同時發(fā)現(xiàn),利用數(shù)值模擬方法解決工業(yè)生產(chǎn)中的復雜問題具有廣闊前景。

1 理論研究

M. Lesani等[1]推導和求解了皮碗球通過液體管道時的二維和三維動力學方程,采用動量方程和能量方程研究了流場對皮碗球運行軌跡的影響。M. Mirshamsi和M. Rafeeyan[2]利用動量方程和狀態(tài)方程計算氣體沿管道流動時的密度、速度和壓力等參數(shù),導出了皮碗球運動的微分方程。同時,他們[3]在二維燃氣管道長清管器的動態(tài)分析中將清管器分成若干個小部分,依照牛頓第二定律對所有小部分分別列出方程,然后將所有方程相加,最終得到清管器總的運動動力學方程。在龍格庫塔法的基礎上,提出了求解運動微分方程的算法。

上述研究給出了皮碗球的總運動方程,為后續(xù)氣體與皮碗球之間關系的研究打下了基礎。在此之后,鄢華士和余恒東[4]對清管器的運行狀態(tài)進行了分析:清管器在運行過程中,速度呈現(xiàn)加速、減速的交替變化。清管器運行的最大阻力出現(xiàn)在初始階段,即清管器從初始速度為零加速至勻速這一過程,此時沿程摩阻最大,清管器所需的氣體驅(qū)動力最大,對應的氣體流量最大。因此,以清管器初始運行階段作為研究對象,計算出驅(qū)動清管器所需的氣體流量,即可作氣源管道設計的依據(jù)。清管器所需的驅(qū)動力由(1)-(4)聯(lián)立計算可得,驅(qū)動清管器所需的氣體質(zhì)量流量則由式(5)計算:

式中:F—清管器受到的氣體驅(qū)動力,N;

F—清管器受到的液體阻力(包括管道摩阻及液位高差引起的阻力),N;

F—清管器與管壁的摩擦力,N;

—清管器的質(zhì)量,kg;

—清管器的重量,N;

Δ—清管器加速時間,s;

—清管器勻速運行速度,m·s-1。

式中:λ—管道的摩阻系數(shù);

ρ—液體化工品的密度,kg·m-3;

L—通球管線的長度(包括管線長度及管件當量摩阻長度),m;

—重力加速度,m·s-2;

H—收球端與發(fā)球端之間的液位差壓,Pa。

式中:—清管器皮碗與管內(nèi)壁的摩擦系數(shù);

—管道內(nèi)徑,m;

—清管器皮碗數(shù)量;

1—清管器中心軸直徑,m;

—清管器的過盈率;

—清管器材料的彈性模量;

—清管器材料的泊桑系數(shù)。

式中:P—清管器初始受到的氣體最大驅(qū)動壓力,Pa。

式中:Q —氣體質(zhì)量流量,kg·s-1;

ρ0—氣體標準狀態(tài)下的密度,kg·m-3;

P0—氣體標準壓力,取 98 000 Pa;

0—氣體標準溫度,K;

T—氣體實際溫度,K。

王鋒等研究發(fā)現(xiàn)[5]皮碗清管器的磨損量與滑動距離、法向載荷、皮碗過盈量、材料的屈服極限等有關。通球管道內(nèi)的輸送介質(zhì)不同,所使用的公式也存在差異。Chen等[6]通過對天然氣管道通球?qū)嶒炈脭?shù)據(jù)的整理分析,提出了表征皮碗球速度與后方推進氣體速度之間關系的公式。Zhang等[7]研究了皮碗球在丘陵天然氣管道中運行過程的動力學模型。根據(jù)響應面法(RSM)的模擬結果,提出了一個預測水平輸氣管道皮碗球穩(wěn)態(tài)速度的公式

在明確了皮碗球與推動氣體之間的關系后,就可以針對管道的實際情況進行皮碗球的選型。白港生等[8]推導出了清管器通過凹陷變形直管段、三通、厚壁及閥門、彎頭時需要滿足的長度、骨架鋼體直徑、皮碗密封間距等基本設計參數(shù),并給出具體的指導公式和經(jīng)驗值。何畏等[9]通過分析清管器運行過程中碟形皮碗相較于直板皮碗的優(yōu)勢,優(yōu)選碟形皮碗作為清管器的皮碗類型;對碟形皮碗進行力學分析并采用有限元軟件 Abaqus 分析得出,過盈量和皮碗唇部厚度與接觸長度呈正相關,皮碗唇部長度與接觸長度呈負相關;接觸應力隨著過盈量和唇部厚度增加而逐漸增大,隨著唇部長度增加而逐漸減小。

2 實驗研究

Chen等[10]采用水平,豎直和傾斜管道相結合的實驗模型,以皮碗球的旁通率為自變量進行多組實驗,對通球過程中管內(nèi)壓力波動、皮碗球速度以及管道末端排出液體流量進行監(jiān)測。實驗結果表明,當管內(nèi)液體體積由2L提高到8L時,旁路管道壓力峰值增加14.4%,無旁路管道壓力峰值增加34%,使用帶旁路的皮碗球裝置可以顯著降低通球過程中管內(nèi)的壓力波動,同時對管內(nèi)液體加載量的變化具有更好的適應性。皮碗球兩端的壓差與旁路分數(shù)無關,與皮碗球運動過程中的沿程阻力有關。皮碗球在管道中的平均速度隨旁路分數(shù)的的增大線性減小。在此之后,Chen等[11]在只包含水平管道的實驗裝置上繼續(xù)使用不同旁通率的皮碗球進行實驗,旁通率從0%提高到2%,液體流出速率降低了75%。

此外,相關企業(yè)結合已有線路對海上油氣田通球過程進行了有效的實驗研究[12-15]。

3 模擬研究

左施施等[16]利用fluent對清管器模型內(nèi)部流場進行仿真,獲得不同形狀的旁通閥在清管器工作時的壓力云圖和湍流動能圖。田宏軍等[17]利用fluent進行數(shù)值模擬。結果表明: 在旁通狀態(tài)下,高速流體通過旁路,對管道內(nèi)壁產(chǎn)生強剪切作用,流過清管器后,在下游形成渦流區(qū)。壓力場分為上游高壓區(qū)、清管器中壓區(qū)、頭部低壓區(qū)和下游穩(wěn)定區(qū); 旁路上和軸線上最大流體速度隨間隙量的增加而減小,流場壓降及所受的驅(qū)動力隨間隙量的增加而急劇減小。Asghar Talbizadeh 和Mohammad Mehdi Keshtkar[18]對通球過程中皮碗球周圍的全湍流流動進行建模,對影響皮碗球兩端壓降的參數(shù)進行了研究,確定了這些基本參數(shù)在通球過程中的最優(yōu)狀態(tài)。

Li等[19]提出了一種基于移動網(wǎng)格法的皮碗球模型,將皮碗球與流體的接觸面視為離散網(wǎng)格節(jié)點,采用精細的運動網(wǎng)格對皮碗球周圍的流場進行了研究,模擬著重考慮了氣體的可壓縮性和通球過程中皮碗球受到的摩擦力。他們[20]為了快速準確地描述皮碗球在通球過程中的運動狀態(tài),首先建立了一個適用于瞬態(tài)與穩(wěn)態(tài)條件下的分析模型,然后建立了一個包含管道進口氣體質(zhì)量流量、皮碗球摩擦力和旁路參數(shù)等變量的模型,用以計算通球過程所需時間。在考慮氣體可壓縮性和摩擦剖面變化的基礎上,在模擬程序中添加摩擦段,研究了皮碗球在通球過程中可能出現(xiàn)的卡死和滑移現(xiàn)象。在設計的摩擦段內(nèi)正摩擦偏差是造成突然減速然后突然停止的主要原因,而當皮碗球進入低摩擦區(qū)域時,氣體可壓縮性則成為皮碗球速度波動或卡死的主要原因。皮碗球與管壁之間的配合關系、皮碗球表面的的摩擦系數(shù)、管材特性(內(nèi)徑、內(nèi)壁粗糙度)和皮碗球速度對摩擦力均有影響。

Liu等[21]采用多相流瞬態(tài)模擬軟件對中國普光氣田P102集氣站通球清管過程進行了研究。左施施等人[22]以可控速清管器為例,建立了清管器的流體區(qū)域泄流流場模型,通過仿真得到了旁通閥不同開度下清管器的流場壓力分布,得到了旁通閥相對開度與壓差的關系曲線(圖1)。隨著旁通閥開度增大,清管器前后壓差先急劇下降再趨于平緩,對附近流場的湍動能影響逐漸變小,清管器在管道中的運行速度逐漸趨于平穩(wěn)。為得到一個合理的開度區(qū)間來平穩(wěn)地調(diào)節(jié)清管器運行速度,擬合了開度與壓差的方程式(6),根據(jù)擬合函數(shù)可對閥門不同開度時的壓差進行估算。

式中,α定義為旁通閥開度,由式(7)計算可得:

式中:S—旁通閥任意狀態(tài)下的泄流面積;

1—旁通閥最大泄流面積。

Liu等[23]采用OLGA和CFD模擬相結合的方法,分析了不同進口質(zhì)量流量和旁通率對提升管底部壓力、皮碗球速度和液相分布形態(tài)等特征參數(shù)的影響。模擬結果表明:立管底部最大壓力隨著皮碗球旁通率的增加先減小后基本保持不變,最小壓力和通球時間增大,壓力變化范圍和通球速度減小; 當旁通率從0%增加到2%時,壓力和皮碗球速度波動下降最顯著,而進口質(zhì)量流量對立管底部壓力的影響不大。過高的旁通率或不足的進口質(zhì)量流量會降低通球操作的效率,并可能不足以提供皮碗球運動所需的壓差。

在長距離氣液管道通球過程中,皮碗球和管道內(nèi)壁之間的摩擦會導致皮碗球皮碗的直徑小于管道的內(nèi)徑,從而產(chǎn)生間隙,間隙的存在對液體的流動有重要影響。為了探究流體在管道中的流動特性,Chen等人[24]基于CFD方法進行了一系列數(shù)值模擬。結果表明:間隙越大,管道內(nèi)殘留液體越多,輸送流體效率越低,此時需要適當增大后方推動氣體速度,進而增大皮碗球速度。

圖1 旁通閥相對開度與壓差的關系[21]

4 結束語

1)利用理論研究得到的公式,可以計算在管道輸送不同介質(zhì)時,皮碗球的速度,后方推動氣體的流量,皮碗球在運動過程中受到的阻力等。通過分析可知,皮碗球速度的主要影響因素是管道中的氣體速度,其次是旁通面積和皮碗球受到的摩擦阻力,而提高皮碗球清管器的旁通率,是提升通球效率的有效手段。

2)通過搭建實驗模型,利用壓力監(jiān)測設備和速度監(jiān)測設備,可以較好的得到通球過程中管道內(nèi)的壓力變化情況和皮碗球的速度,同時使用透明管道可以清晰看到管內(nèi)輸送介質(zhì)的流動情況。實驗結果表明,使用帶旁路的皮碗球裝置可以顯著降低通球過程中管內(nèi)的壓力波動,同時對管內(nèi)液體加載量的變化具有更好的適應性;皮碗球兩端的壓差與旁路分數(shù)無關,與皮碗球運動過程中的沿程阻力有關;皮碗球在管道中的平均速度隨旁路分數(shù)的的增大線性減小。

3)使用數(shù)值模擬軟件對通球過程進行仿真模擬可以得到通球過程中管內(nèi)速度場和壓力場的分布情況,利用流固耦合功能可以分析皮碗球的變形與應力,使用瞬態(tài)模型還可以分析某一時刻的相關結果,仿真模擬所得結果與理論和實驗的結果吻合程度好,這些優(yōu)勢使得模擬方法成為研究通球過程的最佳手段,相關成果在近五年的文獻中得到大量介紹。

[1]LESANI M, RAFEEYAN M, SOHANKAR A. Dynamic Analysis of Small Pig through Two and Three-Dimensional Liquid Pipeline [J]. Journal of Applied Fluid Mechanics, 2012, 5(2): 75-83.

[2]MIRSHAMSI M, RAFEEYAN M. Dynamic Analysis of Pig through Two and Three Dimensional Gas Pipeline [J]. Journal of Applied Fluid Mechanics, 2015, 8(1): 43-54.

[3]MIRSHAMSI M, RAFEEYAN M. Dynamic analysis and simulation of long pig in gas pipeline [J]. Journal of Natural Gas Science and Engineering, 2015, 23:294-303.

[4]鄢華士, 余恒東. 液體化工品管道通球動力氣源設計計算方法研究 [J]. 化學工程與裝備, 2015, (07): 50-52.

[5]王鋒, 聞亞星, 陳金忠. 油氣管道皮碗清管器材料磨損機理研究 [J]. 管道技術與設備, 2019, (01): 40-42.

[6]CHEN J H, HE L M, LUO X M, et al. Characterization of bypass pig velocity in gas pipeline: An experimental and analytical study [J]. Journal of Natural Gas Science and Engineering, 2020, 73.

[7]ZHANG L, ZHOU J W, HE H G. Modeling and Simulation of Pigging for a Gas Pipeline Using a Bypass Pig [J]. Mathematical Problems in Engineering, 2020, 2020.

[8]白港生, 李衛(wèi)全, 張元, et al. 油氣管道清管器的通過性能設計 [J]. 油氣田地面工程, 2021, 40(04): 71-76.

[9]何畏, 孫磊, 黃飛楊, et al. 清管器聚氨酯皮碗接觸性能分析與結構優(yōu)選 [J]. 工程塑料應用, 2021, 49(03): 118-123+129.

[10] CHEN J H, LUO X M, ZHANG H L, et al. Experimental study on movement characteristics of bypass pig [J]. Journal of Natural Gas Science and Engineering, 2018, 59:212-223.

[11] CHEN J H, HE L M, LUO X M, et al. Bypass pigging technology on amelioration of pigging-induced liquid flow: An experimental and modelling study [J]. Ocean Engineering, 2020, 198.

[12] 宮鵬群. 海上無人平臺油氣海管內(nèi)環(huán)路通球問題及解決方法 [J]. 石化技術, 2021, 28(09): 42-43.

[13] 呂長亮. 某海上油田混輸柔性軟管首次在線通球分析研究 [J]. 石油和化工設備, 2021, 24(09): 20-22.

[14] 馬駿, 陳啟明, 鄭武. 某海上油田清管球收球筒工藝流程優(yōu)化改造 [J]. 天津科技, 2022, 49(05): 18-21.

[15] 王玉偉, 郝敬雷. 海洋石油平臺立管陸地通球施工問題探討 [J]. 遼寧化工, 2018, 47(09): 933-935.

[16] 左施施, 張恕遠, 呂光遠, et al. 旁通閥形狀對清管器管內(nèi)流場的影響研究 [J]. 管道技術與設備, 2019, (04): 59-62.

[17] 田宏軍, 魏云港, 曹宇光, et al. 管道清管器旁通流場三維數(shù)值模擬分析 [J]. 石油機械, 2019, 47(01): 120-129.

[18] TALBIZADEH A, KESHTKAR M M. Numerical and experimental study on a bypass pig motion in oil transmission pipeline: a case study [J]. Journal of Petroleum Exploration and Production Technology, 2020, 10(7): 3007-3023.

[19] LI X W, HE L M, LUO X M, et al. Numerical simulation and experimental study of bypass pigging slippage [J]. Ocean Engineering, 2021, 230.

[20] LI X W, HE L M, LUO X M, et al. Transient pigging dynamics in gas pipelines: Models, experiments, and simulations [J]. Ocean Engineering, 2021, 232.

[21] LIU W, HUI L X, SHANG J F, et al. Study on Pigging Process Analyses and Control Methods for Moisture Gathering Lines: Study in Puguang Gas Field, China [J]. Journal of Pipeline Systems Engineering and Practice, 2021, 12(3): 04021011.

[22] 左施施, 張恕遠, 呂光遠, 等. 旁通閥開度對清管器泄流流場影響研究 [J]. 管道技術與設備, 2021, (01): 33-36.

[23] LIU Y T, ZHU X X, SONG J P, et al. Research on bypass pigging in offshore riser system to mitigate severe slugging [J]. Ocean Engineering, 2022, 246:110606.

[24] CHEN S T, ZHANG Y B, WANG X L, et al. Numerical simulation and experiment of the gas-liquid two-phase flow in the pigging process based on bypass state [J]. Ocean Engineering, 2022, 252:111184.

Research on the Movement Mechanism of the Bowl Ball During the Passing Process

1,1,1,1,1,2,2*

(1. Dongguan Shengyuan Petrochemical Co., Ltd., Dongguan Guangdong 523000, China;2. Liaoning Petrochemical University, Fushun Liaoning 113001, China)

The ball passing process uses nitrogen to push the cup ball to deliver the gasoline or diesel in the pipeline to the storage tank.It is an effective means to connect the wharf and the factory to complete the transportation of oil products.The ball is an important equipment in the process of passing the ball, and it can also be used as a pig in the pipeline cleaning operation.The speed of the bowl ball determines the efficiency and safety of the passing process.Affected by the thrust of the gas behind, the frictional resistance of the pipe wall and the resistance of the liquid in the front, the movement of the ball in the pipe is a complex movement process. In this paper, the research on the movement mechanism of the ball during the passing process in the past 10 years was summarized. In the theoretical part, the calculation formulas of gas velocity, cup ball velocity and resistance were given; in the experimental part, the influence of bypass rate on ball passing efficiency was introduced; and in the simulation part, the flow field in the pipeline during ball passing was introduced. The research results can provide some reference for related enterprises to design high efficiency and good economy pigging process. At the same time, it is found that the numerical simulation method has broad prospects for solving complex problems in industrial production.

Ball pass process; Bowl ball; Oil and gas transportation; Bypass pigging; Pig

2022-06-27

王學成(1966-),男,新疆人,研究方向:成品油庫運營管理。

趙磊(1982-),男,講師,博士,儲運動態(tài)系統(tǒng)控制與仿真。

TE83

A

1004-0935(2023)01-0117-04

猜你喜歡
研究
FMS與YBT相關性的實證研究
2020年國內(nèi)翻譯研究述評
遼代千人邑研究述論
視錯覺在平面設計中的應用與研究
科技傳播(2019年22期)2020-01-14 03:06:54
關于遼朝“一國兩制”研究的回顧與思考
EMA伺服控制系統(tǒng)研究
基于聲、光、磁、觸摸多功能控制的研究
電子制作(2018年11期)2018-08-04 03:26:04
新版C-NCAP側面碰撞假人損傷研究
關于反傾銷會計研究的思考
焊接膜層脫落的攻關研究
電子制作(2017年23期)2017-02-02 07:17:19
主站蜘蛛池模板: 久久成人国产精品免费软件| 日本国产一区在线观看| 2048国产精品原创综合在线| 99性视频| 青草视频网站在线观看| 91久久天天躁狠狠躁夜夜| 国产乱子伦视频三区| 久久香蕉国产线看精品| 国产精品美女免费视频大全| 亚洲swag精品自拍一区| 色欲国产一区二区日韩欧美| 好紧太爽了视频免费无码| 色悠久久综合| 免费三A级毛片视频| 精品久久久久无码| 女人18一级毛片免费观看| 国产一区二区三区夜色 | 精品久久综合1区2区3区激情| 欧美成人一区午夜福利在线| 国产黑丝一区| 丁香综合在线| 亚洲黄色视频在线观看一区| 日本欧美中文字幕精品亚洲| 国产乱子伦一区二区=| 波多野结衣视频网站| 国产视频你懂得| 在线人成精品免费视频| 一区二区三区国产| 久久9966精品国产免费| 国产成人无码Av在线播放无广告| 欧美在线国产| 亚洲中文字幕久久精品无码一区| 国产成在线观看免费视频 | 色婷婷狠狠干| 亚洲第一综合天堂另类专| 亚洲无码91视频| 亚洲精品麻豆| 欧美色视频日本| 国产激情国语对白普通话| 精品无码国产自产野外拍在线| 在线免费观看AV| 啪啪国产视频| 亚洲中文字幕av无码区| 日本免费一级视频| 久久永久精品免费视频| 亚洲综合极品香蕉久久网| 精品自窥自偷在线看| 国产丝袜第一页| 四虎永久在线精品影院| 精品国产99久久| 亚洲婷婷在线视频| 欧洲熟妇精品视频| 精品欧美一区二区三区久久久| 欧美人人干| 国产素人在线| 国产精品网址在线观看你懂的| 亚洲无码91视频| 久久亚洲国产视频| 久久九九热视频| 永久毛片在线播| 91高清在线视频| 黄色网站不卡无码| 91小视频版在线观看www| 欧美日韩va| 国产黄在线免费观看| 午夜老司机永久免费看片| 国产成人一级| 亚洲欧美成aⅴ人在线观看| 欧美劲爆第一页| 国产在线一区视频| 久久久久久尹人网香蕉| 波多野结衣一区二区三区AV| 九色免费视频| 久久成人国产精品免费软件 | 亚洲男人的天堂在线| 日本午夜影院| 久久国产香蕉| 国产尤物视频网址导航| 欧美综合中文字幕久久| 国产人成午夜免费看| 国产jizz| 一级毛片免费观看久|