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幾何模型轉換方法初步研究
主站蜘蛛池模板: 日韩在线视频网| 久久精品中文无码资源站| 91视频99| 日韩资源站| a级毛片视频免费观看| 亚洲AⅤ永久无码精品毛片| 日本午夜视频在线观看| 国产在线观看成人91| 久久无码高潮喷水| 久久久久亚洲精品成人网| 日韩人妻无码制服丝袜视频| 欧美一级黄片一区2区| 免费av一区二区三区在线| 美女内射视频WWW网站午夜| 国产杨幂丝袜av在线播放| 91最新精品视频发布页| 青青久久91| 人妻中文久热无码丝袜| 91精品国产福利| 精品91在线| 欧美性精品不卡在线观看| 国产精品视频第一专区| 91亚洲精品国产自在现线| 五月婷婷伊人网| 亚洲精品国产成人7777| AV网站中文| 麻豆精品久久久久久久99蜜桃| 亚洲成网777777国产精品| 国产极品美女在线观看| 国产人人干| 久久久久国色AV免费观看性色| www.亚洲一区| 欧美色视频网站| 激情综合网激情综合| 69视频国产| 国产精品久久自在自线观看| 夜色爽爽影院18禁妓女影院| 高清免费毛片| 在线观看亚洲成人| 国产精品手机在线观看你懂的| 亚洲男人的天堂久久香蕉| 萌白酱国产一区二区| 啪啪永久免费av| 无码日韩视频| 国产尤物在线播放| 亚洲国产精品不卡在线 | 亚洲国产欧美国产综合久久 | 不卡色老大久久综合网| 国产精品页| 国产精品va| 免费国产在线精品一区| 亚洲,国产,日韩,综合一区 | 国产亚洲日韩av在线| 国产高清在线观看| 午夜福利亚洲精品| 亚洲男女在线| 狠狠亚洲婷婷综合色香| 2022国产91精品久久久久久| 久久这里只有精品66| 欧美笫一页| 午夜福利在线观看成人| 青青草原国产av福利网站| 国产男人天堂| 国产成人精品一区二区免费看京| 亚洲美女视频一区| 国产综合无码一区二区色蜜蜜| 亚洲天堂777| 亚洲无码高清一区二区| 国产精鲁鲁网在线视频| 日本AⅤ精品一区二区三区日| 成人国产精品一级毛片天堂| 日韩成人高清无码| a亚洲天堂| 欧美成人综合视频| 亚洲欧洲自拍拍偷午夜色| 天堂网亚洲综合在线| 亚洲中文字幕日产无码2021| 久久久久人妻一区精品| 国禁国产you女视频网站| 亚洲一级毛片| 狂欢视频在线观看不卡| 国产免费久久精品99re丫丫一|