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

特高壓井口節流閥節流性能及結構分析

2025-08-07 00:00:00周天明何強王俊玉戴啟平韓傳軍
石油機械 2025年7期
關鍵詞:節流閥沖蝕楔形

中圖分類號:TE931 文獻標識碼:A DOI:10.12473/CPM.202407066

Zhou Tianming,He Qiang,Wang Junyu,et al.Throtling performance and structural analysis of UHP wellhead throttle valve[J].China PetroleumMachinery,2025,53(7):75-82.

Throttling Performance and Structural AnalysisofUHP WellheadThrottleValve

Zhou Tianming1.2He Qiang1Wang Junyu3Dai Qiping1.2Han Chuanjun' (1.MOEKeyLaboratoryforOilandGasEquipmentatSouthwestPetroleumUniversity;2.CNPCBaoj OilfieldMachineryCo., Ltd.;3.CCDCDownhole Service Company)

Abstract:The throtle valve,,as an ultra-high pressure (UHP)wellhead device,is susceptible to erosion failure.This paper proposes thedesigns of spool and anti-sprint protection sleeve for double-wedge throttle valve,and analyzes theflow field characteristics of nine spools of the double-wedge throtle valve at diferent openings.The results show thatthespools ofthedouble-wedgethrotle valvearebetter-performed in erosionresistance than thatof the single-wedge throtle valve.With the increaseof spool stroke,the velocityof drillingfluid passing throughthe wedge-shapedsurfaceofthespol increases,andthedrillingfluid hasa strong erosion to the spooland the wall surfacebelow thespol.Among the9types of spoolsof the double-wedge throttle valve,the convex symmetrical spool has the best throtling and anti-erosion performance,and its linear regression coeffcientof pressure drop curve is 99. 6% . The cavity size has a significant impact on the performance of the throtle valve:expansion of the cavity size by 30 mm leads to a decrease of maximum velocity by 3.87% on average.After the spool anti-sprint protection sleeve is installed,the maximum fluid velocity decreasesby an average of 30. 35% ,extending the service life of the spool.The research results provide areference forthe design and application of UHP wellhead throtle valves.

Keywords: UHP wellhead; wedge throtle valve;double-wedge surface; erosion and wear; design optimization

0 引言

節流閥是石油和天然氣鉆采過程中使用的關鍵裝置之一,其主要作用是控制井口流體的流量和壓力,以確保鉆井或油氣生產過程的安全。隨著油氣開發向深層發展和井口條件的多樣化,對節流閥的性能和可靠性提出了越來越高的要求。節流閥根據閥芯結構的不同可分為楔形、圓柱、針形和孔板節流閥,其中楔形節流閥是應用較廣泛的一類節流閥。在工作過程中流體沖蝕磨損楔形閥芯,可導致節流閥失效。為了提升楔形節流閥的可靠性,延長其使用壽命,需要對楔形節流閥閥芯結構進行優化設計。

練章華等對楔形閥芯進行了優化設計,通過閥芯結構引導流體流向,避免了流體對閥芯的直接沖擊,延長了楔形節流閥的使用壽命。王國榮等[3針對鉆井作業中高壓流體對閥芯輪廓沖蝕的問題,研究得出了粒子入射角度與沖蝕情況的關系,有助于楔形閥芯的結構設計。屈俊波等4推導了楔形節流閥閥芯行程與流體壓降的理論關系,為不同節流閥的優化設計提供了壓降控制的依據。曹陽等通過試驗驗證了多種類型節流閥在不同工況下的抗沖蝕性能,并提出了楔形節流閥閥芯形狀和材質等的改進思路。李千登等研究了不同節流閥的沖蝕損傷機制,分析了流體壓力和運動速度等對閥芯沖蝕速率的影響,并利用現場試驗對數值模擬結果進行了驗證。文獻[7為延長節流閥使用壽命,將金剛石應用于節流閥閥座和閥芯的表面強化中,試驗結果表明,強化后的節流閥抗沖蝕能力增強了約 50% 。WANGG.R.等[8為降低節流閥的沖蝕磨損,采用響應面法并結合多目標遺傳算法對節流閥結構進行優化設計,優化后沖蝕速率明顯降低。HUG.等研究了閥門開度、鉆井液流量、固相質量和鉆井液密度對節流閥沖蝕特性的影響,為節流閥結構設計優化提供了理論依據。

國內外學者主要在中低壓工況下分析節流閥的沖蝕情況,缺少對其在特高壓工況下的分析。針對特高壓節流閥,其設計和工作環境與常規節流閥存在顯著差異,特高壓節流閥的流體動能和壓力變化更為劇烈,這使得特高壓節流閥面臨著更為嚴峻的沖蝕失效問題。為解決節流閥在特高壓工作過程中流體沖蝕導致失效的問題,筆者對閥門的流體通道進行了優化設計,通過研究雙楔節流閥楔形面多種形狀結構對閥芯壓力控制能力和抗沖蝕能力的影響,優選出了一種性能較好的雙楔節流閥,該閥能減小瞬時壓差,降低流速變化;并提出了一種閥芯防沖刺保護套,增強了對閥芯的保護效果。優化改進方案降低了閥芯沖蝕失效風險,可為延長節流閥使用壽命提供參考。

1流動控制方程

1.1 流量系數

流量系數[]是描述流體控制設備在指定工況下性能的重要參數,流量系數越大,流體經過節流閥后的壓降就越小。流量系數主要取決于節流閥的內部結構和尺寸。其計算式為:

式中: Cv 為流量系數, m2 ; Q 為體積流量, m3/s ρ 為流體密度, kg/m3 ; ΔP 為閥門壓力損失, Pa 。

1.2 流阻系數

流體通過彎管和截面突變的地方會產生擾動、攪拌和能量損失。因此,流阻系數[是衡量節流閥中流體流過節流閥產生壓力損失的重要參數。其計算式為:

Δp=p1-p2=ξρv22/2

式中: ξ 為流阻系數,無量綱; p1 為進口壓力, Pa p2 為出口壓力, Pa ; v2 為局部阻力下游區的平均流速, m/s 。

1.3 控制方程

在鉆井過程中,節流閥內鉆井液的流動是一種非常復雜的湍流流動[12-13]。這里采用標準 k-ε 方程[14-15]模型作為封閉方程,流體在節流閥內流動所滿足的控制方程如下。

質量守恒方程:

動量守恒方程:

k-ε 湍流方程:

式中: fi 為流體加速度矢量, m/s2 ; p 為壓力, Pa u 為運動黏度, m2/s : 為流體的速度矢量, m/s , μ 為動力黏度, Pa?S . μr 為湍流黏性系數, Pa?s ; ε 為耗散率, m2/s3 ; k 為湍流動能, m2/s2

在標準 k-ε 方程中,相關常數一般取為 Cu= 0.09, C1=1.41~1.45 , C2=1.9~1.92 , σ1=1 ,σ2=1.3 。

2 楔形節流閥模型優化設計

楔形節流閥三維模型與流體計算網格模型如圖1所示。為降低計算量和提高計算可靠性,對內部結構倒角進行了簡化處理,對計算域中流體沖擊閥芯處進行了網格加密處理。

圖1楔形節流閥三維模型及網格模型 Fig.1Three-dimensional model and mesh modelof the wedge throttle valve

2.1 閥芯結構設計

單楔節流閥工作原理如圖2a所示。圖2b為雙楔節流閥工作示意圖。

Fig.2Working principleof thewedge throttle valve

單楔閥芯其楔形面與流體入口方向相對,工作過程中,流體首先沖擊到楔形面上,然后在單楔閥芯下端的結構導流[16作用下向下端流出。雙楔閥芯的對稱面法向與流體入口方向相互垂直,工作過程中,一部分流體首先沖擊閥芯2楔形面間的曲面,另一部分流體繞流后進入腔體內部,最后流體分別從2楔形面的下端流出。

雙楔閥芯的兩側楔形面相互對稱,通過改變閥芯楔形面的形狀(平面、凸面、凹面)和楔形面的方位(對稱、迎流面前銳、迎流面前鈍)形成了9種楔形閥芯結構方案,如圖3所示。

圖3不同結構形式的雙楔閥芯Fig.3Double-wedge spoolswith different structural forms

2.2楔形節流閥模型計算參數及邊界條件

在節流管匯工作過程中,井控壓力為 175MPa 流體主要為鉆井液與甲烷。其中,鉆井液的密度為1080kg/m3 ,黏度為 0.012Pa?s ;甲烷密度為0.6679kg/m3 ,黏度為 1.087×10-5Pa?s 。采用多相流模型,將鉆井液設為主相,次相為甲烷。其余邊界條件具體設置如下:

(1)入口邊界采用壓力入口。根據工況,設總壓為 175MPa ,入口動壓 0.1MPa 。湍流強度和水力直徑分別為 4.7% 和 78.5mm 。

(2)出口邊界條件采用質量流量出口。甲烷的質量流量為 0.0072kg/s ,鉆井液的質量流量為43.2kg/s 。

3 計算結果分析

根據前面所述的邊界條件,針對單楔閥芯與雙楔閥芯(平面對稱)進行了不同開度(15、20、25、30、 35mm )下的數值模擬計算,計算結果如表1所示。從表1可知,在相同開度下,雙楔閥芯的最大流體速度比單楔閥芯的最大速度平均下降21. 03% ,抗沖蝕能力優于傳統的單楔閥芯。

基于上述對單、雙楔閥芯的性能分析方法,進一步分析了9種雙楔閥芯節流閥流動情況,分別進行了不同行程開度(15、20、25、30、 35mm ))下的數值模型計算。其中,凸面對稱雙楔節流閥( 25mm 開度)的yz截面及 xz 截面流體速度云圖如圖4所示,其他開度下流體最大速度如表2所示。

圖4 25mm 行程開度凸面對稱節流閥速度云圖Fig.4Velocity nephograms of the convex symmetricalthrottle valve at 25mm opening
表2不同開度下凸面對稱節流閥最大流體速度

由圖4和表2可知:在特高壓工況下,流體在楔形節流閥入口段的速度場分布均勻;當流體進入閥腔后,過流面積急劇縮小,流體被瞬間加速,流體最大速度出現在閥芯楔形面附近;高速流體沿著閥芯楔形面延伸,流動方向與楔形面平行;當流體運動到閥芯下端導流面處,流體的運動方向沿導流面方向開始發生改變,同時,沿楔形面方向的流體不斷沖擊閥芯下方的管壁;隨著閥芯行程的增加,流體通過楔形面處的速度也呈現增加的趨勢。在特高壓工況下,這些流場特征表明,鉆井液對閥芯及閥芯下方的壁面將產生較強的沖蝕磨損。

閥芯行程與整體壓差、最大速度的關系如圖5、圖6和圖7所示。為判斷閥芯行程與壓差曲線的線性程度,通過線性回歸的決定系數 [17]R2 進行判斷。決定系數值介于0\~1之間, R2 值越高,則分析數據的線性程度越高。決定系數計算式為:

圖5平面雙楔閥芯節流閥不同閥芯行程下壓差與流速 Fig.5Differential pressureand flow velocityof the flat doublewedgethrottlevalveatdifferentspoolstrokes
圖6凸面形雙楔閥芯節流閥不同閥芯行程下壓差與流速 Fig.6Differential pressure and flow velocity of the convex double-wedge throttlevalveatdifferent spool strokes
圖7凹面形雙楔閥芯節流閥不同閥芯行程下壓差與流速 Fig.7Differential pressure and flow velocity of the concave double-wedge throttle valve at different spool strokes
表3不同閥芯與壓差間的決定系數表Table 3 Determinant coefficients for spools

式中: yi 為計算數據; 為計算數據的均值; 為預測值。

表3為不同閥芯與壓差間的決定系數表。由表3可知,在特高壓條件下,壓差與閥芯行程的變化關系越接近線性,說明壓力控制效果越好,這一點對特高壓節流裝置尤為重要。在15、20、25、30、 35mm 行程中,平面鈍迎、凸面對稱和凹面鈍迎閥芯在壓力控制上效果較好。其中,凸面對稱閥芯的壓降曲線的線性回歸決定系數最高為99.6% 。

在最大速度與閥芯行程的關系圖中,最大速度越小,則流體對閥芯所造成的沖蝕影響越小。其中,平面對稱、凸面對稱和凹面對稱閥芯各行程速度均較小,因此能有效降低特高壓流體的沖蝕影響。

結合不同閥芯在各行程下的速度分布情況,凸面對稱閥芯具有較好的壓力控制能力,適合特高壓節流的應用場合,且沖蝕影響的范圍和流體最大速度都較小,因此優選凸面對稱閥芯用于特高壓工況。

4閥腔體尺寸對節流閥性能的影響

當節流閥內部腔體結構尺寸發生改變時,流體對節流閥的沖蝕情況也會改變[18]。閥腔尺寸增大,減弱了流體流入閥腔后直接沖擊閥芯的強度,閥腔內部的流場分布規律也會發生變化。為分析節流閥來流方向腔體尺寸變化對節流閥閥芯沖蝕程度的影響,針對行程與最大速度關系中速度較小的平面對稱、凸面對稱和凹面對稱閥芯,均在20mm 行程開度下開展數值模擬分析。腔體尺寸分別沿 x 軸的正方向和負方向各擴大5、10、 15mm (共擴大10、20、 30mm ),如圖8所示。

圖8節流閥腔體結構尺寸擴大方向Fig.8Expansion direction of the throttlevalve’scavity size

凸面對稱閥芯節流閥擴大 30mm 后的速度云圖如圖9所示。3種閥芯節流閥擴大尺寸與最大流體速度的關系如圖10所示。

圖10中,隨著 x 方向尺寸的擴大,節流閥中最大速度隨之降低,當 x 方向尺寸擴大 30mm 時,最大速度平均下降 3.87% 。在腔體尺寸擴大的同時,閥體結構的強度會被削弱。因此,通過改變腔體尺寸減緩沖蝕效果需確保閥體結構強度滿足使用要求。

圖12凸面閥芯節流閥不同行程下的壓差Fig.12Differential pressure of the convex throttle

5 防沖刺閥芯保護套對節流閥性能的影響

為減弱特高壓流體對節流閥的沖蝕作用,延長特高壓工況下節流閥閥芯的使用壽命,在節流閥閥芯處增加一個防沖刺保護套,如圖11所示。圖11a為無防沖刺保護套的節流閥模型,圖11b為有防沖刺保護套的節流閥模型。通過增加防沖刺保護套,可避免流體進入節流閥后直接沖擊閥芯。

圖11節流閥閥芯局部示意圖Fig.11Partial illustration of the throttle valve spool"

分析最大速度較小的凸面形雙楔閥芯在不同行程開度(15、20、25、30、 35mm )下防沖刺保護套對節流閥的保護效果。不同閥芯行程下,壓差與最大速度如圖12和圖13所示。其中,虛線表示有防沖刺閥芯保護套的計算結果,實線則表示無防沖刺保護套。

計算結果表明,在特高壓工況下,防沖蝕保護套的設計顯著改變了流體運動特性。在有防沖刺閥芯保護套的情況下,當流體通過下閥座時,過流面積再次減小,流體速度再次增加[19]。最大速度產生在閥芯保護套附近,因此可以認為閥芯保護套代替閥芯受到了直接沖蝕。

從圖12的壓差曲線可知,在15、20、25、30、35mm 行程內,3種凸面形閥芯在防沖刺保護套的作用下壓差下降明顯,且有防沖刺保護套的凸面銳迎閥芯整體壓差最小。從圖13可知,在特高壓工況下,安裝防沖刺保護套后,3種凸面閥芯的最大速度相對無防沖刺閥芯保護套時最大速度平均下降 30.35% 。其中,凸面對稱雙楔閥芯的最大流體速度最低。

6結論

為解決單楔節流閥在特高壓工況下易產生沖蝕失效的問題,提出了雙楔形閥芯和閥芯防沖刺保護套方案,通過數值模擬證明了其抗沖蝕性能提升明顯,并根據楔形面平、凸和凹的特點對比分析了9種不同的雙楔閥芯的流場。研究結果表明:(1)綜合壓力控制能力和沖蝕情況,9種雙楔閥芯中凸面對稱閥芯性能最好,其壓降曲線的線性回歸決定系數最高。(2)隨著節流閥腔體 x 方向尺寸的擴大,節流閥中流體最大速度隨之降低;當 x 方向尺寸擴大30mm 時,節流閥中流體最大速度平均下降3.87% ,能夠在一定程度上減弱特高壓條件下的沖蝕磨損。(3)通過應用節流閥閥芯防沖刺保護套,流體最大速度平均下降 30.35% ,緩解了流體對節流閥閥芯的直接沖蝕,延長了節流閥使用壽命。(4)研究結果可為 175MPa 節流閥的設計和應用提供參考,有助于提高節流閥在特高壓工況下的穩定性、可靠性,延長其使用壽命。

參考文獻

[1]周俊澤,練章華,王裕海,等.節流閥閥芯失效分析及其結構改進[J].石油機械,2019,47(12):131-138.ZHOUJZ,LIANZH,WANGYH,etal.Failureanalysis and structure optimization of throttle valve core[J].China Petroleum Machinery,2019,47(12):131-138.

[2] 練章華,劉干,龔建文,等.高壓節流閥流場分析及其結構改進[J].石油機械,2004,32(9):22-24,41.LIANZH,LIUG,GONGJW,etal.Fluidfieldanalysis of high pressure throttlevalve and its structureimprovement[J].China Petroleum Machinery,2004,32(9):22-24,41.

[3]王國榮,錢權,楚飛,等.高壓流體入射角對節流閥材料的沖蝕預測及驗證[J].中國機械工程,2017,28(14):1652-1657,1663.WANGGR,QIANQ,CHUF,etal.Erosion predic-tionandverificationforthrottlevalvematerialsaboutinci-dentangles at high pressure fluids [J].China Mechani-calEngineering,2017,28(14):1652-1657,1663.

[4] 屈俊波,閆鐵,孫曉峰,等.控壓鉆井楔形節流閥節流壓降規律研究[J].化工機械,2018,45(3):328-332.QUJB,YANT,SUNXF,et al.Study on chokepressure drop of wedge - shaped throttle valvesused inmanaged pressure drilling [J]. Chemical Engineeringamp;Machinery,2018,45(3):328-332.

[5]曹陽,方曉慶,張俊藍,等.鉆井用節流閥抗沖蝕性能的實驗評價[J].天然氣工業,2020,40(5):89-93.CAOY,FANG XQ,ZHANGJL,et al. Experimen-tal evaluation on the anti-erosion performance of throttlevalves for drilling[J].Natural Gas Industry,2020,40(5):89-93.

[6]李千登,管志川,樊朝斌,等.高壓液固兩相流節流閥耐沖蝕特性研究[J].石油機械,2022,50(2):89-94,101.LI Q D,GUAN Z C,FAN Z B,et al. Erosion-resistingcharacteristic study of throttle valve in high-pressureliquid-solid Flow [J].China Petroleum Machinery,2022,50(2):89-94,101.

[7]WHEELER D W,WOOD R JK,HARRISON D,etal.Application of diamond to enhance choke valve life inerosive duties[J].Wear,2006,261(10):1087-1094.

[8]WANG GR,CHUF, TAO SY,et al. Optimizationdesign for throttle valve of managed pressure drillingbased on CFD erosion simulation and response surfacemethodology[J].Wear,2015,338-339:114-121.

[9]HUG, ZHANG P,WANG G R,et al. Performancestudy of erosion resistance on throttle valve of managedpressure drilling [J]. Journal of Petroleum Science andEngineering,2017,156:29-40.

[10]王國榮,王德貴,劉清友,等.一種井控節流閥防刺短節的特性分析[J].石油機械,2011,39(11):10-12.WANGGR,WANGDG,LIUQY,et al.Thecharacteristic analysis of well-control throtle valve witherosion resistance sub [J]. China Petroleum Machin-ery,2011,39(11):10-12.

[11]TREVOR G,HARSA M,DAN W,et al. Incompress-ible flow through choke valve:an experimental andcomputational investigation [J]. Journal of Fluids andStructures,2022,113:103669.

[12]汪周華,范琨鵬,趙建飛,等.非常規儲層微納米孔隙介質中流體相態研究進展[J].世界石油工業,2024,31(3):68-77.WANGZH,FANKP,ZHAOJF,etal.Researchprogress of fluid phase behavior in micro -nano porousmedia of unconventional reservoirs[J].World PetroleumIndustry,2024,31(3):68-77.

[13]楊永飛,王金雷,王建忠,等.基于VOF方法的超臨界二氧化碳:水兩相流動孔隙尺度數值模擬[J].天然氣工業,2023,43(3):69-77.YANGYF,WANGJL,WANGJZ,et al.Pore-scalenumerical simulation ofsupercritical CO2 -brine two-phase flow based on VOF method[J].Natural GasIndustry,2023,43(3):69-77.

[14] 賈林,艾志久,王彪,等.節流閥閥芯振動失效分析及結構改進[J].石油礦場機械,2013,42(4):29-32.JIAL,AIZJ,WANGB,etal.Throttle spool vibra-tion failure analysis and structure improvement [J].OilFieldEquipment,2013,42(4):29-32.

[15]HUOVINEN M, KOLEHMAINEN J, KOPONENP,et al.Experimental and numerical study of a chokevalveina turbulentflow[J].FlowMeasurementandInstrumentation,2015,45:151-161.

[16]宋保健,明鑫,孫凱,等.楔形節流閥沖蝕磨損規律研究[J].石油礦場機械,2018,47(4):30-34.SONGBJ,MINGX,SUNK,etal.Studyontheerosion law of wedge-shaped throttle valve[J].OilFieldEquipment,2018,47(4):30-34.

[17] 李瑞月,鐘山,張漓杉,等.多元線性回歸-分光光度法測定電鍍排放廢水中重金屬[J].中國環境監測,2023,39(5):226-235.LIRY,ZHONGS,ZHANGLS,etal.Simultane-ousdetermination of heavy metals in electroplatingwastewater by multiple linear regression spectrophotom-etry[J].Environmental Monitoringin China,2023,39 (5): 226-235.

[18] 常漢卿.高壓氣井節流閥設計及研究[D].青島:中國石油大學(華東),2020.CHANG HQ. Design and research of high pressure gaswell throttlevalve [D].Qingdao:china UniversityofPetroleum (Huadong),2020.

[19]朱軍龍.多通道流場條件下的水下節流閥沖蝕研究[D].北京:中國石油大學(北京),2019.ZHUJL. Erosionresearch ofsubseachokevalvewithmulti-channels flow field [D].Beijing:China Univer-sityofPetroleum(Beijing),2019.

第一

作者簡介:周天明,教授級高級工程師,生于1978年,2001年畢業于西南石油大學機械設計制造及其自動化專業,現從事石油鉆井裝備技術研究工作。地址:(721002)陜西省寶雞市。電話:(0917)3462166。email: zhoutianming@cnpc.com.cn。

通信作者:韓傳軍。email:hanchuanjun@swpu.edu.cn

收稿日期:2024-07-17 修改稿收到日期:2024-12-21(本文編輯劉鋒)

猜你喜歡
節流閥沖蝕楔形
基于CADe_SIMU軟件的機電綜合仿真研究
基于CT平掃的預測模型對肺梗死與胸膜下肺炎的鑒別診斷價值
基于無傳感器無刷直流電動機的節流控制方法
汽車電器(2025年7期)2025-08-10 00:00:00
離心壓縮機葉輪沖蝕失效的可靠性研究
石油機械(2025年7期)2025-08-07 00:00:00
以上海市森蘭楔形綠地一綠帶項目為例探討森林氧吧型公園綠地建設
主站蜘蛛池模板: 国产簧片免费在线播放| 青青操视频在线| 中文字幕人妻无码系列第三区| 亚洲国产精品成人久久综合影院| 国产人成乱码视频免费观看| 久久免费精品琪琪| 免费看av在线网站网址| 国产福利微拍精品一区二区| 在线精品欧美日韩| 国产自在自线午夜精品视频| a网站在线观看| 婷婷午夜天| 国产免费福利网站| 亚洲天堂网在线观看视频| 色噜噜中文网| www.99在线观看| 亚洲中久无码永久在线观看软件| 国产精品永久免费嫩草研究院| 91在线日韩在线播放| 日本欧美一二三区色视频| 国产精品乱偷免费视频| 五月婷婷综合网| 欧美专区日韩专区| 亚洲第一精品福利| 中国国产一级毛片| 亚洲精品无码在线播放网站| 国产av无码日韩av无码网站| 丰满人妻一区二区三区视频| 丁香五月婷婷激情基地| 91小视频在线播放| 婷婷色一区二区三区| 亚洲欧美精品日韩欧美| 亚洲人成影视在线观看| 欧美a在线| 色精品视频| 国产成人精品一区二区三区| 最新国产高清在线| 亚洲一级无毛片无码在线免费视频| 国产毛片久久国产| 国产成人亚洲日韩欧美电影| 日韩精品一区二区三区免费在线观看| 97成人在线视频| 真实国产精品vr专区| 欧美精品亚洲精品日韩专区va| 亚洲最新地址| 日韩AV手机在线观看蜜芽| 亚洲丝袜第一页| 视频在线观看一区二区| 亚洲人成网站色7777| 激情国产精品一区| 国产在线观看高清不卡| 99re热精品视频国产免费| 精品人妻一区无码视频| 中文纯内无码H| 青青青视频蜜桃一区二区| 在线观看亚洲天堂| 婷婷色婷婷| 毛片免费视频| 欧美亚洲国产日韩电影在线| 久久精品国产精品一区二区| 欧美日韩国产成人在线观看| 亚洲精品国产首次亮相| 精品久久久久久中文字幕女| 国产精品漂亮美女在线观看| 伊人激情久久综合中文字幕| 无码在线激情片| 久久久久国产一区二区| 伊人久久精品无码麻豆精品 | 国产尤物视频网址导航| 婷婷亚洲视频| 欧美日韩成人在线观看 | 亚欧乱色视频网站大全| 国产00高中生在线播放| 久久国产拍爱| 在线观看91精品国产剧情免费| 亚洲熟妇AV日韩熟妇在线| 国产精品美乳| 亚洲欧美激情小说另类| 国产欧美亚洲精品第3页在线| 在线播放国产99re| 蜜臀av性久久久久蜜臀aⅴ麻豆 | 被公侵犯人妻少妇一区二区三区|