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

滑溜水用減阻劑室內性能測試與現場摩阻預測

2017-08-28 21:03:34孟磊周福建劉曉瑞楊釗石華強長慶油田分公司油氣工藝研究院低滲透油氣田勘探開發國家工程實驗室西安7006中國石油大學北京非常規天然氣研究院北京099
鉆井液與完井液 2017年3期
關鍵詞:實驗

孟磊, 周福建, 劉曉瑞, 楊釗, 石華強(.長慶油田分公司油氣工藝研究院低滲透油氣田勘探開發國家工程實驗室,西安 7006;.中國石油大學(北京)非常規天然氣研究院,北京 099)

滑溜水用減阻劑室內性能測試與現場摩阻預測

孟磊1, 周福建2, 劉曉瑞1, 楊釗2, 石華強1
(1.長慶油田分公司油氣工藝研究院低滲透油氣田勘探開發國家工程實驗室,西安 710016;2.中國石油大學(北京)非常規天然氣研究院,北京 102299)

孟磊, 周福建, 劉曉瑞, 等.滑溜水用減阻劑室內性能測試與現場摩阻預測[J].鉆井液與完井液, 2017, 34(3):105-110.

MENG Lei,ZHOU Fujian,LIU Xiaorui,et al.Friction reducers for slick water: performance evaluation and friction prediction[J].Drilling Fluid & Completion Fluid,2017,34(3):105-110.

經驗預估滑溜水壓裂現場的減阻效果具有差異較大的不足,而以往的減阻劑減阻性能評價實驗多側重于相關產品的對比優選,對其測試結果如何應用于現場實際的說明并不多見。針對以上問題,搭建了室內環路摩阻測試系統,依據Prandtl-Karman定律,通過清水湍流率定出測試用的3條直管管徑分別為10.46,7.59和5.86 mm。對不同管徑,在一定泵送排量梯度下對5種不同濃度的DR-800減阻劑溶液的摩阻壓降進行了測定。對于0.07%和0.10%濃度的減阻劑溶液,將測試壓降換算成摩阻系數,與考慮黏度項計算出的雷諾數相對應,可以很好地擬合于Virk漸近線,證明了DR-800減阻劑的優良減阻性能。在定量測試減阻性能后,采用摩阻放大法中的阻力速度法,將室內實驗測得的降阻比與阻力速度進行回歸得出表達式,確定出各阻力速度所對應的具體減阻率大小。運用該式,通過迭代計算,可以得到每種工況下相對確定的減阻率,比起經驗預測,精度明顯提高。因此,上述的室內研究方法對減阻劑的現場應用具有指導意義。

滑溜水;減阻劑;阻力速度;摩阻預測Virk漸近線

美國的“頁巖氣革命”揭開了非常規油氣開發的面紗,自1997年,Mitchell能源公司首次將滑溜水應用在Barnett頁巖的壓裂作業中,滑溜水壓裂液成為非常規儲層開發中最價廉而且質優的改造液[1-2]。滑溜水中的水含量可達99%以上,其他添加劑盡管含量較低卻發揮著重要的作用,特別是減阻劑的加入,可有效地降低壓裂液流動時的摩阻系數,從而降低了壓裂施工的摩阻損耗[3]。雖然由此帶動了減阻劑的研發與評價實驗[4-11],但是目前很少有研究針對于實驗室測試結果與現場施工減阻效果的對應進行展開。實驗室與現場條件總是存在差別,最明顯的便是影響減阻的“管徑效應”[12]。筆者選用了減阻劑DR-800,將實驗室摩阻測試結果與Virk最大減阻漸近線進行對比,確定了其優良的減阻性能[13]。結合石油行業壓裂液紊流摩阻放大方法中的阻力速度法,將其降阻比與阻力速度變化進行了對應[14]。運用阻力速度法對現場壓裂設計參數的模擬計算結果精度優于以往的經驗預測,實現了室內減阻性能測試對現場應用的指導。

1 減阻劑室內性能評價

1.1 實驗器材

實驗材料為DR-800,一種陰離子型聚丙烯酰胺減阻劑,溶劑為自來水。實驗時該減阻劑體積分數的加量分別為0.01%、0.03%、0.05%、0.07%和0.10%。實驗儀器主要為一套包含3根3 m長可切換測試管段的環路摩阻測試裝置,見圖1。為消除兩端變徑對流動的影響,取其中間2.5 m段開孔設置測壓點。供選壓差計量程分別為0.5 和5 MPa,依據測試壓差大小和實驗精度可以進行調整切換。實驗動力源為一臺變頻控制的螺桿泵,泵的最大排量為2.5 m3/h。實驗循環管路中的彎管接頭通徑均大于19 mm,連接構成循環的膠皮軟管選用19 mm標準軟管,盡量消除測試段之外的剪切破壞。在循環測試管段的出口端下游,設置了質量流量計。通過數字流量信號的反饋,可以按照需求控制變頻器,調節輸出確定的流體排量,為后期獲取同一流速下的摩阻壓降提供了便利。

圖1 減阻率環路摩阻測試系統結構示意圖

1.2 管徑的率定

由于管道直徑對實驗結果的影響顯著,實驗首先用清水做了管徑的率定。清水的摩阻系數λ在水力光滑區符合普蘭德-卡門(Prandtl-Karman)定律[15]:

式中,λ為達西摩阻系數;Re為雷諾數。

通過控制不同的排量,在25±3 ℃下流經實驗管道。摩阻系數及雷諾數按照下式計算:

式中,△P為實驗管段壓降,Pa;d為管道內徑,m;l為實驗段長度,m;ρ為流體密度,kg/m3;v為斷面平均流速,m/s;Q為排量,m3/s;μ為流體動力黏度,mPas。

圖2 清水湍流實驗結果

由圖2可知,實驗數據點很好地分布在率定基準線附近,說明實驗管路具有很高的精確度,能夠為減阻實驗可靠性提供有力證據。

1.3 減阻性能測試實驗

首先配制50 L待測液體,然后利用摩阻測試儀分次先后在粗、中、細3條管徑條件下連接好環路系統。在已有清水基準數據的基礎上,對不同濃度的減阻劑溶液按照相同的排量梯度調整螺桿泵工況,得到相應的直管段壓力降。實驗的減阻率可用如下公式計算:

式中,DR為減阻劑的減阻率,%;ΔP0為同一流速下未加減阻劑時的摩阻壓降,Pa;ΔPDR為同一流速下加入減阻劑時的摩阻壓降,Pa。

以0.07%濃度中管變排量實驗為例,減阻率隨流速的變化見圖3。

圖3 中管0.07%濃度液樣減阻率隨流速的變化

由圖3可知,DR-800減阻劑減阻性能良好,在管路內平均流速為2 m/s時就已經達到了55%以上的減阻效果,當流速大于10 m/s后,減阻率穩定在70%以上,且尚未出現剪切稀釋帶來的明顯下降。

1.4 減阻性能評價與分析

高分子聚合物稀溶液的減阻效應存在著極限減阻現象,根據Virk等人的研究,高聚物的湍流減阻以Virk線為漸近線[16-17]。Virk漸近線的表達式為:

為評價減阻劑DR-800的減阻性能,引入該線作為比較基準和依據。首先,在每次配液及實驗后分別及時取樣,利用烏氏黏度計測定各液體的黏度,結果見表1。由表1可知,以新配液體標定出的黏度值為準,結合公式(3)可以計算出考慮黏度后的減阻液體的Re。然后利用公式(2)將0.07%、0.10%濃度液體在粗、中、細3根管路中變排量下測試的所有壓降實驗數據換算成確定的λ值,并與雷諾數一一對應。最后把所得數據點與Virk漸近線以及Prandtl-Karman定律基準線繪制在一起,得到圖4。由圖4可知,所有的實驗點都分布在Virk漸近線的附近,說明了減阻率隨流速變化沒有下降,濃度為0.07%及0.10%的液樣,即便更換管路,DR-800減阻劑均貼近了湍流減阻中的最大減阻。早在1987年,劉鶴年等人的研究中曾指出,國產的聚丙烯酰胺減阻劑(FLA-801)在濃度大于50 mg/L時,進入湍流區即達到最大減阻,最大減阻流動的摩阻曲線與Virk漸近線符合。這也支持了DR-800減阻劑具有優良的減阻性能。同時該結果還說明,由Virk漸近線與Prandtl-Karman定律基準線可以反算水力光滑區不同管徑、不同流速下使用DR-800的減阻率大小。由圖4還可以看出,部分離散點位置較Virk漸近線偏低,說明對于滑溜水這樣的假塑性流體,其黏度標定初值與流動中的表觀黏度存在差異,為了得到與λ對應的準確Re值,減阻流體流動中的表觀黏度需要在下一步的工作中確定。

表1 各濃度液體實驗前后的黏度值匯總

圖4 不同濃度減阻劑溶液湍流實驗結果

2 利用阻力速度法實現現場摩阻預測

室內實驗表明,在水力光滑區DR-800的λ與Re有很好的相關性,借助Virk漸近線也可以量化減阻劑在室內的減阻性能。但是從實驗室到現場,由于管徑接近10倍的差異,Re值也相差了一個數量級,現場的流動早已超越了水力光滑區。陳鵬飛等人的研究也指出,剪切速率的相似原則不能完全有效地評價減阻性能[18]。因此,減阻的相似模擬需要引入其他的參數。

在石油工業壓裂液紊流摩阻的放大方法中,Lord將放大方法分為3大類:①博溫(Bowen)方法;②阻力系數與雷諾數法;③阻力速度法。但前2種方法都只適用于某些條件,而阻力速度法可以無條件地滿足所有不同的實驗情況[19-21]。通過模擬計算,現場與室內的阻力速度也在0~0.5 m/s的同一范圍之間,因此采用阻力速度法來實現實驗室管路到施工實際的預測。

2.1 阻力速度定義

Whitsitt發現降阻比公式(σ=ΔPDR/ΔP0)在不同的直徑下可以與阻力速度相關聯。溶劑的阻力速度定義為:

式中,v*

s為阻力速度,m/s;Rws為管壁上的剪應力,kg/m2;ρ為密度,kg/m3。

2.2 阻力速度與降阻比在本實驗的應用

結合阻力速度與降阻比的定義,對室內實驗濃度為0.07%的減阻劑溶液的實驗數據進行再處理。將采集到的壓差、管徑、減阻率代入公式(6)、(7)計算相應的阻力速度,可以得到粗、中、細3條管路減阻率隨阻力速度的變化,見圖5。由圖5可知,對于每條管路都可以得到減阻率與阻力速度之間的一條平滑曲線,隨著管徑的變粗,曲線的相對位置逐漸往圖版的左上方移動。表現出更好的減阻性能,即對于不同管徑管路在同一阻力速度下,管徑越粗減阻越好。同時,隨著管徑的變細,該管路覆蓋的阻力速度范圍變大。

圖5 3條管路減阻率隨阻力速度的變化

將圖5中的數據點做相應轉換,可以得出變形后的對應散點圖,結果見圖6。

將該圖6的數據進行回歸,可以得出降阻比與阻力速度相應的數學表達式。

依據該公式計算的阻力速度與降阻比以及減阻率的對應表見表2。

2.3 現場摩阻預測

某油田一口氣井的壓裂設計中可以采集到如下參數:壓裂方式為套管壓裂,套管尺寸為139.7 mm,井深為3 000 m;滑溜水壓裂液密度為1 000 kg/m3。清水的沿程摩阻在10~32 MPa之間。滑溜水的減阻率預計在60%~80%。

根據以上參數,假定運用滑溜水后實現了65%的減阻,則對應3 000 m管段的沿程摩阻為3.5~11.2 MPa。將這些數據代入公式(6),可以計算出對應的阻力速度范圍為0.201 8~0.361 m/s。由該阻力速度范圍比照,可得在不同的工況條件下,使用0.07%濃度的減阻劑溶液預測得到的減阻率為69.2%~71.5%,利用所得的減阻率迭代計算,可以分別得到各工況下更為精確的阻力速度值以及減阻率值。對比之前60%~80%的預計精度明顯提高,對比65%綜合減阻的假設,該計算結果的差別小于3%。由此可得,利用室內性能測試得到的降阻比與阻力速度的對應關系,可以更為合理的預測滑溜水壓裂現場的沿程摩阻大小,為施工參數的制定提供參考。

表2 降阻比、減阻率與阻力速度的對應關系

3 結論

1.經過管徑率定后的清水,實驗數據可以很好地分布在Prandtl-Karman基準線附近,說明該系統具有很高的精確性。

2.濃度為0.07%的DR-800溶液在實驗管路內流速大于10 m/s后,減阻率在70%以上。濃度為0.10%的DR-800溶液動力黏度小于3 mPas。黏度的實驗處理結果可以很好地擬合于Virk漸近線,證明了相對較高濃度的DR-800具備極限減阻的優良性能。

3.通過模擬計算,現場與室內的阻力速度值分布在0~0.5 m/s的同一范圍之間,比Re、剪切速率的相似性更好;通過對3條管路實驗數據的回歸,得到了濃度為0.07%的DR-800溶液降阻比與阻力速度的表達式。

4.根據現場壓裂設計所采集的參數,運用上述表達式,通過迭代試算的方法可以得到每種工況下相對確定的減阻率,結果比單純的經驗估計更為合理。因此,可以運用該方法對現場施工進行指導。

[1]杜凱,黃鳳興,伊卓,等.頁巖氣滑溜水壓裂用降阻劑研究與應用進展[J].中國科學:化學,2014,44(11):1696-1704. DU Kai,HUANG Fengxing,YI Zhuo,et al.Recent advances on friction reducer for slick-water fracturing of shale gas reservoirs[J].Science China:Chemistry,2014,44(11):1696-1704.

[2]MAYERHOFER M J,RICHARDSON M F,WALKER JR R N,et al.Proppants we don't need no proppants[C]// SPE annual technical conference and exhibition.Society of Petroleum Engineers,1997.

[3]蔣官澄,許偉星,李穎穎,等.國外減阻水壓裂液技術及其研究進展[J].特種油氣藏,2013,20(1):1-6. JIANG Guancheng,XU Weixing,LI Yingying,et al. Technology of fracturing fluids with friction-reducing water overseas and its research progress[J].Special Oil and Gas Reservoir,2013,20(1):1-6.

[4]魏娟明,劉建坤,杜凱,等.反相乳液型減阻劑及滑溜水體系的研發與應用[J].石油鉆探技術,2015,43(1):27-32. WEI Juanming,LIU Jiankun,DU Kai,et al.The development and application of inverse emulsified friction reducer and slickwater system[J].Petroleum Drilling Techniques,2015,43(1):27-32.

[5]王娟娟,劉通義,趙眾從,等.減阻劑乳液的合成與性能評價[J].應用化工,2014,43(2):308-310. WANG Juanjuan,LIU Tongyi,ZHAO Zhongcong,et al. Synthesis and performance evaluation of dragreduction agent emulsion[J].Applied Chemical Industry,2014,43(2):308-310.

[6]肖博,張士誠,雷鑫,等.頁巖氣藏清水壓裂減阻劑優選與性能評價[J]. 油氣地質與采收率,2014, 21(2):102-105. XIAO Bo,ZHANG Shicheng,LEI Xin,et al. Optimization and evaluation of drag reducer of slickwater fracturing in shale gas reservoir[J].Petroleum Geology and Recovery Efficiency,2014,21(2):102-105.

[7]陳鵬飛,劉友權,鄧素芬,等.頁巖氣體積壓裂滑溜水的研究及應用[J].石油與天然氣化工,2013,42(3):270-273. CHEN Pengfei,LIU Youquan,DENG Sufen, et al. Research and application of slick water for shale vblume fracturing[J].Chemical Engineering of Oil and Gas,2013,42(3):270-273.

[8]陳朝霞.高分子減阻劑性能評價與表征[D].烏魯木齊:新疆大學,2007. CHEN Zhaoxia.The evaluation and token of the macromolecule drag reduction agent[D].Urumchi:Xinjiang University,2007.

[9]吳雷澤,薛苗,魯紅升.礁石壩區塊頁巖氣壓裂用滑溜水體系研究[J].鉆采工藝,2015,38(6):84-86. WU Leize,XUE Miao,LU Hongsheng.Research of slick water for shale gas fracturing in jiaoshiba block[J]. Drilling and Production Technology, 2015, 38(6): 84-86.

[10]龍學莉,葉智,程曉亮,等.減阻劑體系的室內研究及應用[J].當代化工,2015,44(11):2563-2565. LONG Xueli,YE Zhi,CHENG Xiaoliang,et al. Research and application of the drag reducing agent system[J].Contemporary Chemical Industry,2015,44(11):2563-2565.

[11]彭飛,方波,盧擁軍,等.四元兩性疏水聚合物溶液流變和減阻性能[J].鉆井液與完井液,2016,33(1):92-96. PENG Fei,FANG Bo,LU Yongjun,et al.Rheology and drag reduction of hydrophobic amphoteric quadripolymer solution[J].Drilling Fluid & Completion Fluid,2016,33(1):92-96.

[12]ZHOU M J,SUN H,QU Q,et al.An effective model of pipe friction prediction from laboratory characterization to field applications for friction reducers[C]//SPE Annual Technical Conference and Exhibition. Society of Petroleum Engineers,2011.

[13]王鴻勛,張士誠.水力壓裂設計數值計算方法[M].北京: 石油工業出版社,1998:36-40. WANG Hongxun,ZHANG Shicheng.Numerical calculation method of hydraulic fracturing design[M]. Beijing:Petroleum industry press,1998:36-40.

[14]朱蒙生.管流添加劑減阻的實驗與機理研究[D].哈爾濱:哈爾濱工業大學,2009. ZHU Mengsheng.Experiment and mechanism study on drag-reduction by additives in pipe flow[D].Harbin:Harbin Institute of Technology,2009.

[15]黃趾海.新型滑溜水壓裂液研究[D].成都:西南石油大學,2014. HUANG Zhihai.Study on a new type of slick-water fracturing fluid[D].Chengdu:Southwest Petroleum University,2014.

[16]溫慶志, 高金劍,劉華,等.滑溜水攜砂性能動態實驗[J].石油鉆采工藝,2015,37(2):97-100. WEN Qingzhi,GAO Jinjian,LIU hua,et al.Dynamic experiment on slick-water prop-carrying capacity[J].Oil Drilling & Production Technology,2015,37(2):97-100.

[17]劉通義, 黃趾海, 趙眾從,等.新型滑溜水壓裂液的性能研究[J].鉆井液與完井液,2014,31(1):80-83. LIU Tongyi, HUANG Zhihai, ZHAO Zhongcong, et al. Performance study on a new slickwater fracturing fluid[J]. Drilling Fluid & Completion Fluid,2014,31(1):80-83.

[18]陳鵬飛,唐永帆,劉友權,等.頁巖氣藏滑溜水壓裂用降阻劑性能影響因素研究[J].石油與天然氣化工,2014,43(4):405-408. CHEN Pengfei,TANG Yongfan,LIU Youquan,et al. Influencing facctors of friction reducer in shale slickwater fracturing[J].Chemical Engineering of Oil & Gas,2014,43(4):405-408.

[19]王海燕, 邱曉惠, 翟文.頁巖儲層壓裂減阻劑減阻機理研究[J].鉆井液與完井液, 2015,32(4) :75-77. WANG Haiyan,QIU Xiaohui,ZHAI Wen. Functioning mechanism of drag reducer used in shale reservoir fracturing[J].Drilling Fluid & Completion Fluid,2015,32(4) :75-77.

[20]馮虎, 徐志強. 塔里木油田克深區塊致密砂巖氣藏的儲層改造技術[J].石油鉆采工藝,2014,36(5):93-96. FENG Hu,XU Zhiqiang.Reservoir reconstruction technology for tight sandstone gas pool in Keshen Block of Tarim Oilifeld[J].Oil Drilling & Production Technology,2014,36(5):93-96.

[21]周成香,周玉倉, 李雙明,等.川東南頁巖氣井壓裂降壓技術[J].石油鉆探技術,2014,42(4):42-47. ZHOU Chengxiang,ZHOU Yucang,LI Shuangming,et al.Fracturing pressure reducing technique for shale gas well in the Southeastern Sichuan[J].Petroleum Drilling Techniques,2014,42(4):42-47.

Friction Reducers for Slick Water: Performance Evaluation and Friction Prediction

MENG Lei1, ZHOU Fujian2, LIU Xiaorui1, YANG Zhao2, SHI Huaqiang1
(1. Research Institute of Oil & Gas Technology of PetroChina Changqing Oilfield CompanyNational Engineering Laboratory of Low Permeability Oil & Gas Field Exploration and Development, Xi’an, Shaanxi 710016;2. The Unconventional Natural Gas Institute, China University of Petroleum, Beijing 102299)

Deficiencies exist in empirically predicting the friction reducing efficiency of slick water fracturing fluids in field operation, yet the experimental methods presently in use for evaluating the performance of friction reducers are only focused on the performance comparison and product optimization, while the illustration of how to apply the test results to field operation is rarely seen in literatures. To resolve this problem, a laboratory loop-line friction measurement system has been established. According to the Prandtl-Karman law, the diameters of three straight pipes used in the system were determined to be 10.46 mm, 7.59 mm and 5.86 mm using turbulent flow rate of clear water. For any one of the pipe diameters, the frictional pressure drops of five different concentrations of DR-800 (a friction reducer) water solutions were tested at a certain flow rate. For the 0.07% and 0.10% DR-800 solutions, the test pressure drops were converted into frictional coefficients, and then corresponded with the Reynolds numbers calculated which took into account viscosity. In this way the frictional coefficients can be satisfactorily fitted with Virk asymptote, and the excellent performance of the DR-800 was proven. After quantitatively measuring the friction reducing performance, using friction-velocity method (a submethod of the so-called friction amplification method), the ratios of friction reduction and velocities-under-friction obtained from the experiments can be regressed into a mathematical expression that can be used to work out the rate of friction reduction corresponding to a specific velocity-under-friction. Using this expression, a comparatively affirmative rate of friction reduction under a work conditioncan be determined through iteration. The calculation has better precision than empirical prediction, and the laboratory study method is instructive to the field application of friction reducers.

Friction reducers for Slick Water; friction-velocity; Friction prediction Virk asymptote

TE357.12

A

1001-5620(2017)03-0105-06

2017-1-23;HGF=1703C7;編輯 王超)

10.3969/j.issn.1001-5620.2017.03.021

中國石油天然氣集團公司科學研究與技術開發項目“井筒工作液基礎理論關鍵技術研究”(2014A-4212)和低滲透油氣田勘探開發國家工程實驗室開放課題“致密儲層體積壓裂減阻劑減阻機理研究”(ZX20150057)資助。

孟磊,工程師,碩士,1983年生,現在從事油氣田開發壓裂液研究工作。電話 (029)86590661/ 15009285802;E-mail:mlei1_cq@petrochina.com.cn。

猜你喜歡
實驗
我做了一項小實驗
記住“三個字”,寫好小實驗
我做了一項小實驗
我做了一項小實驗
記一次有趣的實驗
有趣的實驗
小主人報(2022年4期)2022-08-09 08:52:06
微型實驗里看“燃燒”
做個怪怪長實驗
NO與NO2相互轉化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
主站蜘蛛池模板: 免费三A级毛片视频| 潮喷在线无码白浆| 亚洲天堂视频在线观看免费| 伊人狠狠丁香婷婷综合色| 精品91自产拍在线| 尤物国产在线| 国产精品无码影视久久久久久久 | 蜜桃视频一区二区| 成人免费网站久久久| 青青草一区二区免费精品| www亚洲精品| m男亚洲一区中文字幕| 亚洲综合在线网| 亚洲欧洲日本在线| 国产一区二区免费播放| 欧美三级自拍| 一本一道波多野结衣av黑人在线| 全部毛片免费看| 日韩欧美在线观看| 成人看片欧美一区二区| 色播五月婷婷| 亚洲男人天堂2020| 亚洲精选高清无码| 欧美不卡二区| 亚洲一级色| 视频一本大道香蕉久在线播放| 青青草原国产| 乱系列中文字幕在线视频| 久久精品丝袜| 91精品国产一区| 最新国产精品第1页| 亚洲第一成年网| 亚洲av无码成人专区| 国产精品久久久久久久久久久久| 国产精品成人观看视频国产 | 伊人天堂网| 一区二区在线视频免费观看| 国产99欧美精品久久精品久久| 色噜噜在线观看| 91免费国产高清观看| 欧美a级完整在线观看| 日韩精品亚洲人旧成在线| 国产美女视频黄a视频全免费网站| 一级毛片免费高清视频| 狠狠色噜噜狠狠狠狠色综合久| 亚洲无码在线午夜电影| 成人午夜亚洲影视在线观看| 国产成人8x视频一区二区| 欧美yw精品日本国产精品| 亚洲综合狠狠| 亚洲成年人片| 波多野结衣视频网站| 中文国产成人久久精品小说| 亚洲啪啪网| 国产a v无码专区亚洲av| 亚洲床戏一区| 色网站免费在线观看| 国产精品手机在线观看你懂的| 69综合网| 在线免费亚洲无码视频| 国产欧美日韩视频怡春院| 91精品啪在线观看国产91九色| 激情五月婷婷综合网| 亚洲区欧美区| 亚洲欧美一区二区三区蜜芽| 国产亚洲欧美另类一区二区| 玖玖免费视频在线观看| 亚洲无线视频| 就去色综合| 精久久久久无码区中文字幕| 久综合日韩| 无码福利日韩神码福利片| 五月婷婷亚洲综合| 欧美97欧美综合色伦图| 日韩欧美91| 日韩一区二区三免费高清| 99精品伊人久久久大香线蕉 | 免费无遮挡AV| 欧美成人亚洲综合精品欧美激情| 亚洲中文字幕在线一区播放| 亚洲中文无码h在线观看 | 国内毛片视频|