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

大斜度井投撈操作動力學模型及投撈器參數優化設計

2016-07-21 07:52:15宋顯民黃重光劉士英范曉東中國石油冀東油田公司中國石油玉門油田公司
石油鉆采工藝 2016年2期
關鍵詞:優化模型

宋顯民黃重光劉士英范曉東.中國石油冀東油田公司;.中國石油玉門油田公司

大斜度井投撈操作動力學模型及投撈器參數優化設計

宋顯民1黃重光1劉士英2范曉東2
1.中國石油冀東油田公司;2.中國石油玉門油田公司

大斜度井投撈是氣舉采油檢閥和水井分注調配的關鍵技術。目前大斜度井投撈研究僅限于工藝和工具設計,缺乏相關的理論模型指導工作。因此,以大斜度井氣舉投撈操作為研究對象,利用動力學原理及不考慮鋼絲拖曳和考慮鋼絲拖曳兩種情況下投送器下沖速度公式,建立了能量形式的大斜度井投撈操作動力學模型,分析了下沖剩余能與投撈段井斜、投送器直徑、投送器長度的關系;并以投撈技術適應的最大井斜為目標建立了投送器參數優化設計模型,進行了相關設計計算,為目前大斜度井投撈工藝和相關工具研究提供了理論基礎,對今后投撈技術的改進與提高提供參考。

大斜度井;投撈;下沖速度;能量關系;動力學模型;參數優化

對于大斜度井氣舉閥和大斜度井注水堵塞器的井下投撈技術,目前國內外工藝及工具研究較多:如宋顯民、劉清友等研究了大斜度井鋼絲投撈式氣舉采油技術及主體工具[1];林書國、劉洪海、劉德基、張澤華等人分別在國外引進技術基礎上研究形成了可投撈氣舉采油工具及投撈技術[2-8];王巨鵬、李日寧等在最初偏心分注井投撈技術基礎上分別進行了斜井和大斜度井投撈測配技術研究[9-11]。投撈技術中投撈運動過程的基礎理論研究比較少,雖然有鞏志蓉、唐高峰等人研究了液力投撈(無鋼絲)測試調配時儀器下井運動規律[12-13],王旭升等人研究了注水井管柱內測調儀下行阻力計算[14],但是目前國內外仍沒有建立起鋼絲投撈操作的動力學模型,和起投撈器參數優化設計模型。為此,以大斜度氣舉投撈操作為對象進行了相關研究。

1 投送器大斜度井段下沖速度分析Undershoot velocity of fishing tool in highly-deviated well

目前大斜度氣舉井氣舉閥和注水井配水堵塞器的井下投撈過程相同,以氣舉閥井下投送為例進行說明。首先在油管內下放投送器至預定工作筒下部;然后上提投送器,上提高度隨成功投送所需下沖距離而增加,此時投送器的導向塊徑向釋放;然后地面絞車下放鋼絲,投送器開始向下沖擊,當進入到工作筒內時,受工作筒導向筆尖驅動而沿軸向旋轉,并使氣舉閥對準工作筒偏孔,緊接著氣舉閥快速進入到工作筒偏孔內并被鎖定;最后地面絞車上提鋼絲,投送器與氣舉閥之間的剪釘被剪斷,投送器被鋼絲上提至井口[1]。以投送器在工作筒上部開始向下沖擊至氣舉閥完全鎖定于工作筒偏孔內的過程為研究對象,進行動力學模型研究。

1.1不考慮鋼絲拖曳的投送器在斜井段下沖速度

Undershoot velocity of fishing tool in deviated well section without consideration to wireline towing

投送器下沖速度為[12-14]

其中

式中,v為不考慮鋼絲拖曳時的下沖速度,m/s;v*11為不考慮鋼絲拖曳時的臨界下沖速度,m/s;t為下沖時間,s;D為投送器所在處油管內徑,m;d為投送器外徑,m;l為投送器長度,m ;α為井斜角度,°;k為浮力系數;μ為摩擦系數,取0.15。

1.2考慮鋼絲拖曳的投送器斜井段下沖速度

Undershoot velocity of fishing tool in deviated well section with consideration to wireline towing

以投送器與鋼絲整體作為研究對象,在下沖過程中沿油管軸向受6個力作用。(1)投送器和鋼絲的重力沿井筒軸向分力,為下沖動力;(2)投送器與鋼絲受垂直向上的浮力,為下沖阻力;(3)沿井筒軸向上受摩擦力,為下沖阻力;(4)井內鋼絲上端因井口密封和絞車拖曳而受到拉力,為下沖阻力,實測值為40 N;(5)鋼絲在充滿液體的油管內下行時,受到黏滯力,為下沖阻力。氣舉管柱洗井后下部原油由于密度差而向上運移,根據現場實際情況,原油占油管容積的5%、水占油管容積的95%。根據實際情況,鋼絲直徑為3.2 mm、鋼絲下深1 000 m時,黏滯力與速度的關系為fsw=24.93v;當以油為介質時,設定油的動力黏度系數為水的80倍,則黏滯力與速度的關系fsw=1994v;則鋼絲受到油管內流體黏滯力與速度關系為[15]

fsw=0.05×1 994v+0.95×24.93v=123.4v (4)

式中,fsw為鋼絲受到油管內流體黏滯力,N。

(6)投送器在充滿液體的油管內下行時受到壓差阻力,為下沖阻力[ 12-14]

式中,f為投送器壓差阻力,N;M為投送器質量,kg。1.2.1 投送器斜井段下沖速度分析 將投送器與鋼絲整體為研究對象,建立斜井段鋼絲投送過程的動力學關系式為

其中

積分后得

其中

式中,m為鋼絲質量,kg;v2為考慮鋼絲拖曳時的下沖速度,m/s。

1.2.2下沖長度與時間關系分析

將式(10)代入式(12)積分后得

式中,S為下沖長度,m。

1.2.3極限情況分析 隨著下沖時間延續,式(10)的極限情況為

式中,v2*為考慮鋼絲拖曳時下沖臨界速度,m/s。

可見投送器下沖過程是由初始速度向極限速度無限趨近的過程,并最終以極限速度進行勻速運動。設定實際速度達到極限速度0.95倍時,便進入到勻速運動階段了,則下沖達到極限速度的時間為

式中,t*為鋼絲下沖到達極限速度的時間,s。

由以上公式計算可知,用長度1.3 m、直徑40 mm的投送器在投撈深度1 000 m、初速度0向下沖擊:當井斜從0°到70°時逐漸增加時,投送器從初始速度0到達極限速度的時間也逐漸增加,從1.49 s增加至1.89 s;投送器從初始速度0到達極限速度的下沖距離逐漸減小,由4.22 m減至1.15 m。

1.3實際投撈過程中下沖速度公式選擇

Selection of undershoot velocity formulas in actual fishing operations

不考慮鋼絲拖曳的投送器下沖速度公式適用于投送器下沖速度低于其上部鋼絲運行速度的情況,此時鋼絲對投送器不產生拖曳作用,投送器實際下沖速度即是投送器在充滿液體的斜井油管中的自由落體速度的情況;考慮鋼絲拖曳的投送器下沖速度公式適用于以投送器為研究對象的下沖速度高于其上部的鋼絲運行速度,此時實際下沖速度等于有鋼絲拖曳情況下的投送器下沖速度。令不考慮鋼絲拖曳的下沖速度與考慮鋼絲拖曳的下沖速度相等,可得2種下沖速度公式選擇的臨界條件

實際投送器直徑變化范圍在0.03 m到0.046 m之間、投送器長度變化范圍在1 m到10 m之間,此時單獨投送器下沖時的速度高于其上部鋼絲的運行速度,實際速度公式采用考慮鋼絲拖曳的投送器下沖速度公式,所以基于考慮鋼絲拖曳的投送器下沖速度公式建立了投送過程的動力學模型和參數優化設計模型。

1.4下沖速度公式驗證

Verification of undershoot velocity formulas

將工作筒與油管連接,下入到實驗井1 050 m處,斜度33°,完成井口;然后將長度1.80 m、直徑0.04 m的模擬投送器用鋼絲攜帶下入油管中,下入到1 000 m處,分別測試不同下沖長度時的下沖速度;之后在另外2口斜度更大的井上重新下入管柱,也在1 000 m處測試不同下沖長度時下沖速度。測試結果見表1,可以看出以實測速度為基準,計算速度相對于實測速度的誤差基本在2%以內,不同井斜情況下實際速度與計算速度吻合程度較好,說明模型符合實際情況。

表1 不同井斜和下沖長度條件下模型計算速度與現場測試速度對比Table 1 Comparison between simulated velocity and measured velocity

2 大斜度井投送系統動力學模型Kinetic model of fishing in highly-deviated well

2.1模型建立

Model building

投送器在下沖投送過程中經歷3個關鍵點,一是投送器被上提一定高度并開始下沖的起始位置,設為A位置;二是工作筒的主體頂部位置,設為B位置;三是工作筒底端位置,設為C位置。投送器由工作筒下部10 m上提到A位置(工作筒上部20 m處)之后,投送器底端由A位置開始向下沖擊,此下沖過程經歷2個階段:在A-B段,投送器在內徑62 mm油管內向下沖擊;在B-C段,投送器在內徑46 mm的工作筒主通道內向下沖擊。

(1)投送器在A-B階段所獲得的速度和能量。在A位置處,投送器下沖速度為0 m/s,下沖動能為0 m/s。由于A-B階段距離為20 m,即下沖距離為20 m。計算所得的臨界下沖長度小于20 m,因此投送器在B點獲得的下沖速度等于臨界速度,則投送器在B點所獲得能量為

式中,vb為B位置處投撈器的下沖速度,m/s;Eb為B位置處投撈器的下沖動能,J。

(2)投送器在C位置處所獲得的速度和能量。由式(5)得

則C點所具有的動能為

式中,vc為C位置處投送器下沖速度,m/s;Ec為C位置處投送器下沖動能,J。

(3)大斜度井投送系統動力學模型建立。按k取值0.87,B-C段長度0.2 m,D取0.046 m,獲得實際情況下的C位置處動能。氣舉閥投入偏孔的條件為:投送器在向下沖擊至C位置時的剩余能要大于氣舉閥投入偏孔需要的能耗。實測氣舉閥入偏孔所需能耗為10 J,大斜度井投送系統動力學模型為

2.2大斜度井能量關系式的應用分析

Application of energy relation for highlydeviated well

投送器直徑對下沖剩余能的影響:按井斜0°、長度1.8 m、下沖長度20 m計算,結果如圖1所示。

圖1 投送器直徑與下沖剩余能之間關系曲線Fig.1 Diameter of fishing tool vs. residual undershoot energy

隨著直徑從0.024 m增加到0.044 m,油管內下沖獲得能量快速增加;而工作筒壓差阻力消耗能初期緩慢增加,當直徑0.036 mm后快速增加;而有效下沖剩余能先由51.7 J逐漸上升到最大值108.9 J,之后逐漸下降為6.3 J。

投送器長度對下沖剩余能的影響:按井斜0°、直徑0.04 m、下沖長度20 m計算,結果如圖2所示。

圖2 投送器長度與下沖剩余能之間關系曲線Fig. 2 Length of fishing tool vs. residual undershoot energy

隨著投送器長度從1 m上升到4 m時,投送器在油管內下沖獲得能量由75 J快速增加到593 J;而工作筒壓差阻力消耗能由47 J緩慢增加到130 J;有效下沖剩余能由27.7 J逐漸增加到462.7 J。

投送段井斜對下沖剩余能的影響:按鋼絲長度1 000m、鋼絲直徑3.2 mm、投送器長度1.8 m、投送器直徑0.04 m、下沖長度20 m計算,結果如圖3所示。隨著投送段斜度逐漸增加,油管內下沖獲得能量由166.9 J快速下降到1.9 J;而工作筒壓差阻力消耗能由70.9 J下降到0.8 J;有效下沖剩余能由96 J緩慢下降到1.1 J。當投入偏孔所需能量為8 J,有效下沖剩余能曲線與投入偏孔所需能量曲線交于62°。說明當井斜小于62°時,投送器獲得的有效下沖剩余能大于氣舉閥投入偏孔所需能量,滿足投送操作能量條件;當井斜大于62°時,投送器獲得的有效下沖剩余能小于氣舉閥投入偏孔所需能量,不滿足投送操作的能量條件。

圖3 投送段井斜與下沖剩余能之間關系曲線Fig.3 Deviation of fishing section vs. residual undershoot energy

3 投送器參數優化模型建立及求解Building and solution of parameter optimization model of fishing tool

將目前氣舉投撈技術所適應的最大井斜作為優化設計所追求的目標,投送器直徑、長度作為需要優化的變量,建立投送器參數優化模型。

3.1投送器參數優化模型目標函數

Obective function for parameter optimization model of fishing tool

依據大斜度井投送系統動力學模型,氣舉閥投入偏孔最低條件的動能形式為

式中,M*為投送系統下沖當量質量,kg。

將鋼絲直徑3.2 mm、長度1 000 m時,投送器優化設計的目標函數為

3.2投送器參數優化模型約束條件

Constraints for parameter optimization model of fishing tool

投送器結構包括上下接頭、投撈手臂、凸輪、內槽、閥底支撐裝置、氣舉閥投入口袋等部分,為滿足工藝要求,工程上所需的最小長度為1.2 m。由于井口防噴器、井眼軌跡的限制,投送器限料最大長度為2 m。投送器的投撈手臂及氣舉閥收縮時需要完全容納于內槽中,工程上最小直徑為40 mm。投送器需要從工作筒主通道順利通過且無卡阻現象,且工作筒主通道內徑為46 mm,故投送器最大直徑為46 mm。因此投送器約束條件為

1.2≤l≤2 0.040≤d≤0.046 (24)

3.3優化結果

Results of optimization

應用優化設計原理,采用matlab編程,最終優化得到在可行域內的最優值為:投送器最優長度2 m、投送器最優直徑40 mm時、最優目標函數值64.7°。圖4為圖示法的大斜度井投送器參數優化結果,在投送器長度和直徑約束條件所形成的矩形可行域內,通過可行域的最大斜度線是64.7°線,與可行域交于此時點(40 mm,2 m)。

圖4 大斜度井投送器參數優化結果Fig.4 The delivery tool parameter optimization results in graphic method

4 結論Conclusions

(1)研究了投送器在斜井段上有鋼絲拖曳和無鋼絲拖曳的下沖速度公式,并提出了2種公式適用條件,基于投送器和工作筒的可能外徑參數,有鋼絲拖曳的投送器下沖速度公式符合實際情況。

(2)基于有鋼絲拖曳的下沖速度公式,建立了能量形式的大斜度井投送系統動力學模型,分析了下沖剩余能與投送段井斜、投送器直徑、投送器長度的關系,為目前大斜度井投撈工藝和相關工具研究提供了理論基礎。

(3)以目前投撈技術適應的最大井斜為目標函數,建立了大斜度井投送器參數優化設計模型,并進行了優化計算。投送器最優長度2 m,投送器最優直徑40 mm,最大井斜64.7°,為投撈工具和工作筒尺寸參數設計提供依據。

References:

[1]宋顯民,劉清友,孫成林,肖國華,胡慧莉,王金忠,靳彥

卿.大斜度井鋼絲投撈式氣舉采油技術[J].石油學報,2014,35(5):979-986.

SONG Xianmin, LIU Qingyou, SUN Chenglin, XIAO Guohua, HU Huili, WANG Jinzhong, JIN Yanqing. Gaslift production technology of wire casting and fishing in highly deviated wells[J]. Acta Petrolei Sinica, 2014, 35 (5), 979-986.

[2]林書國,李桂英,陳偉良.國產化可投撈氣舉采油工具的研制[J].石油鉆采工藝,1989,11(6):85-90.

LIN Shuguo, LI Guiying, CHEN Weiliang. Development of domestic production tools for gas lift [J]. Oil Drilling & Production Technology, 1989, 11(6): 85-90.

[3]袁業啟,陶利萍,何漢坤,常衛群,范海玲.影響氣舉投撈工藝的因素及對策[J].鉆采工藝,2004,27(3):99-100.

YUAN Yeqi, TAO Liping, HE Hankun, CHANG Weiqun,FAN Hailing. The affecting factors and countermeasures of gas lift and fishing [J]. Drilling Production Technology, 2004, 27 (3): 99-100.

[4]劉洪海,樊文利,魏瑞金,張俊法.氣舉投撈工藝技術及其應用[J].西部探礦工程,2001,13(6):63-64.

LIU Honghai, FAN Wenli, WEI Ruijin, ZHANG Junfa. Gas-lift salvage technique and application [J].West-China Exploration Engineering, 2001, 13(6): 63-64.

[5]中國石油天然氣股份有限公司.氣舉注水雙用工作筒:中國,201220176949.3[P].2012-12-12.

PetroChina Company Limited. Gas-lifting & waterinjection barrel: China, 201220176949.3 [P].2012-12-12.

[6]中國石油天然氣股份有限公司.氣舉閥投撈頭:中國,201220050177.9[P].2012-11-14.

PetroChina Company Limited. Fishing head of gas-lifting valve: China, 201220050177.9 [P]. 2012-11-14.

[7]劉德基,溫偉豪,王良,馬帥幫.投撈式氣舉工作筒的研制及應用[C].全國氣舉技術研討會論文集,2007:107-111.

LIU Deji, WEN Weihao, WANG Liang, Ma Shuaibang. Development and application of retrievable gas lift working cylinder [C]. National Gas Lift Technology Conference, 2007: 107-111.

[8]張澤華.氣舉投撈工具串的研制及其應用[J].石油礦場機械,1992,21(1):22-27.

ZHANG Zehua.Development and application of tool string for delivery and fishing in gas lift wells[J].Oil Field Equipment 1992, 21(1): 22-27.

[9]宋顯民,孫成林,陳雷,田晶蒙,彭丙軍,吳訓平,宋利彬.深斜井偏心定量早期分注技術[J].石油鉆采工藝,2007,29(增刊1): 78-81.

SONG Xianmin, SUN Chenglin, CHEN Lei, TIAN Jingmeng, PENG Bingjun, WU Xunping, SONG Libin. Technology of eccentric quantitative early separate injection in deep deviated wells[J]. Oil Drilling & Production Technology, 2007, 29(S1): 78-81.

[10]王巨鵬,史贊絨,宋健.分注井投撈測配技術的研究及應用[J].中國石油和化工,2013,33(10):23-26.

WANG Jupeng, SHI Zanrong, SONG Jian. Research and application of the measurement and matching technology on sub-injection well [J].China Petroleum and Chemical Industry, 2013, 33 (10): 23-26..

[11]李日寧.大斜度偏心分注井測試配套工具的改進與應用[J].石油天然氣學報,2013,35(3):28-31.

LI Rining.High angle eccentric points improvement and application note well testing supporting tools [J].Oil and Gas Technology, 2013, 35(3): 28-31.

[12]金少先,鄭廣漢.井液黏滯力對油管疲勞斷裂的影響[J].大慶石油學院學報,1996,20(1):96-99.

JIN Shaoxian, ZHENG Guanghan. The influence of fluid viscosity on the fatigue fracture of tubing [J].Journal of Daqing Petroleum Institute, 1996, 20(1): 96-99.

[13]鞏志蓉,李欽道.超深天然氣井的測試實踐與分析[J].天然氣工業,1989,9(4):47-50.

GONG Zhirong, LI Qindao. The testing practice and analysis of ultra deep natural gas well[J]. Natural gas industry, 1989, 9(4): 47-50.

[14]王旭升,王春生,王朝麗.注水井管柱內測調儀下行阻力計算[J].石油天然氣學報,2012,34(8):147-149.

WANG Xusheng, WANG Chunsheng, WANG Chaoli. Downstream resistance calculation of measuring instrument in injection well pipe column[J]. Journal of Oil and Gas Technology, 2012, 34(8): 147-149.

(修改稿收到日期 2016-02-14)

〔編輯 李春燕〕

Kinetic model of fishing in highly-deviated well and parameter optimization of fishing tool

SONG Xianmin1, HUANG Chongguang2, LIU Shiying2, FAN Xiaodong2
1. PetroChina Jidong Oilfield Company, Tangshan, Hebei 063004, China; 2. PetroChina Yumen Oilfield Company, Jiuquan, Gansu 735019, China

Fishing in highly-deviated well is critical for valve inspection in gas lift oil production and for separated injection and allocation in water injection wells. The current researches in this aspect only focus on design of processes and tools, but there is no appropriate theoretical model available as a guide. In this paper, an energy-based kinetic model of fishing was established for highlydeviated wells which are developed by gas lift, using the kinetic principle and the formula of tool undershoot velocity with and without consideration to wireline towing. The relationships between the residual undershoot energy and the well deviation in fishing section and diameter/length of the fishing tool were analyzed. Moreover, a parameter optimization model of fishing tool was built on the basis of the maximum well deviation that the fishing technique can adapt to, and relevant calculations were made for the design. The results provide theoretical basis for current researches on fishing in highly-deviated wells and relevant tools, and also contribute to the future improvement of fishing technologies.

highly-deviated well; fishing; undershoot velocity; energy relation; kinetic model; parameter optimization.

SONG Xianmin, HUANG Chongguang, LIU Shiying, FAN Xiaodong. Kinetic model of the delivery operating in highangle wells and tools parameter optimization[J].Oil Drilling & Production Technology, 2016, 38(2): 238-243.

TE355

A

1000 -7393( 2016 ) 02 -0238-06

10.13639/j.odpt.2016.02.021

宋顯民(1971-),西南石油大學石油礦場機械專業在讀博士研究生,現從事采油工藝研究工作,高級工程師。通訊地址:(063004)河北省唐山市路北區光明西里冀東油田鉆采院。電話:13903154635或0315-8768028。E-mail:songxm@petrochina.com.cn 或1925184311@qq.com

引用格式:宋顯民,黃重光,劉士英,范曉東.大斜度井投撈操作動力學模型及投撈器參數優化設計[J].石油鉆采工藝,2016,38(2):238-243.

猜你喜歡
優化模型
一半模型
超限高層建筑結構設計與優化思考
房地產導刊(2022年5期)2022-06-01 06:20:14
民用建筑防煙排煙設計優化探討
關于優化消防安全告知承諾的一些思考
一道優化題的幾何解法
由“形”啟“數”優化運算——以2021年解析幾何高考題為例
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
3D打印中的模型分割與打包
FLUKA幾何模型到CAD幾何模型轉換方法初步研究
主站蜘蛛池模板: 色欲综合久久中文字幕网| 永久成人无码激情视频免费| 国产成人亚洲欧美激情| 真实国产乱子伦视频| 中文字幕在线观看日本| 美女视频黄又黄又免费高清| 久久青青草原亚洲av无码| 亚洲欧美人成电影在线观看| 日本妇乱子伦视频| 91精品国产91久久久久久三级| a级毛片网| 欧美在线三级| 亚洲乱强伦| 黄色三级毛片网站| 婷婷午夜影院| 一级一毛片a级毛片| 日韩高清欧美| 亚洲人网站| 蜜臀AV在线播放| 欧美一级在线看| 精品综合久久久久久97超人| 欧美性爱精品一区二区三区 | v天堂中文在线| 久久精品国产电影| 欧美黑人欧美精品刺激| 亚洲香蕉久久| 无码 在线 在线| 5555国产在线观看| 尤物视频一区| 国产女同自拍视频| 亚洲精品第一页不卡| 青青久久91| 久精品色妇丰满人妻| 国产亚洲精品97AA片在线播放| 在线播放国产99re| 欧美国产日韩在线| 亚洲欧美成人网| 99re热精品视频国产免费| 丰满人妻被猛烈进入无码| 久久99精品久久久久纯品| 毛片基地视频| 国产精品久线在线观看| 国内精品免费| 国产主播喷水| 中国毛片网| 浮力影院国产第一页| 亚洲综合一区国产精品| 99国产精品一区二区| 日韩免费毛片视频| 国产无遮挡猛进猛出免费软件| 午夜啪啪福利| 成人免费一区二区三区| 国产欧美精品一区二区| 毛片a级毛片免费观看免下载| 久久久精品无码一区二区三区| 欧美啪啪一区| 国产三级成人| 草逼视频国产| 国产欧美日韩精品综合在线| 一本一道波多野结衣一区二区| 蝴蝶伊人久久中文娱乐网| 欧美第一页在线| 日本欧美视频在线观看| 特级毛片8级毛片免费观看| 九九九九热精品视频| 欧美成人午夜在线全部免费| 国产精品夜夜嗨视频免费视频| 伊人激情综合| 狼友av永久网站免费观看| 国产精品手机视频| 国产 在线视频无码| 91久久夜色精品国产网站| 色悠久久久| 亚洲国产精品无码AV| 日韩精品少妇无码受不了| 日韩中文精品亚洲第三区| 成人va亚洲va欧美天堂| 国产香蕉97碰碰视频VA碰碰看| 精品一區二區久久久久久久網站| 福利视频久久| 91小视频版在线观看www| 欧美国产菊爆免费观看|