• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于高精度波前快速推進法的頁巖氣儲層動用分析

    2014-10-03 16:31:26黃世軍滕柏路程林松
    東北石油大學學報 2014年4期
    關(guān)鍵詞:動用運移氣藏

    黃世軍,滕柏路,程林松,艾 爽,賈 振

    (中國石油大學(北京)石油工程學院,北京 102249)

    基于高精度波前快速推進法的頁巖氣儲層動用分析

    黃世軍,滕柏路,程林松,艾 爽,賈 振

    (中國石油大學(北京)石油工程學院,北京 102249)

    為了快速準確模擬頁巖氣儲層開發(fā)過程中的動用情況,利用高精度波前快速推進法對波前運移程函方程進行求解,通過復合線性流模型對多級壓裂水平井模型中的氣體流態(tài)進行分區(qū),實現(xiàn)飛行時間與真實時間的換算,得到儲層動用隨時間變化示意圖及儲層動用范圍與時間的定量關(guān)系.對比波前快速推進法與高精度波前快速推進法的計算結(jié)果,發(fā)現(xiàn)在相同時間內(nèi),利用高精度波前快速推進法計算時儲層動用范圍更大.在裂縫總長相等的情況下,通過對比不同裂縫形態(tài)下的儲層動用發(fā)現(xiàn),在裂縫間的壓力達到連通前,儲層動用程度與裂縫形態(tài)無關(guān);當裂縫之間壓力連通后,儲層動用與裂縫形態(tài)密切相關(guān).

    頁巖氣;高精度波前快速推進法;儲層動用;波前運移;裂縫形態(tài)

    0 引言

    隨著鉆井及壓裂技術(shù)的進步,頁巖氣藏得以經(jīng)濟高效的開采,頁巖氣也成為北美地區(qū)主要的供氣來源,北美頁巖氣的成功開發(fā)也激起全球開采頁巖氣的高潮[1-2].中國具有豐富的頁巖氣資源,經(jīng)濟價值巨大,加大對頁巖氣的開發(fā)力度可以有效緩解能源緊張的局勢[3-7].在頁巖氣藏生產(chǎn)過程中,儲層動用范圍是描述頁巖氣藏生產(chǎn)能力的一個重要參數(shù),大部分研究人員認為壓裂改造區(qū)(SRV)外側(cè)的動用很少.Fisher M K等利用微地震監(jiān)測技術(shù)研究Barnett頁巖壓裂水平井的裂縫分布情況[8],Mayerhofer M J等闡述從微地震圖數(shù)據(jù)估算SRV的方法,實現(xiàn)利用微地震監(jiān)測技術(shù)研究頁巖氣藏壓裂水平井動用規(guī)律[9].徐兵祥等利用雙孔瞬態(tài)線性流動模型及典型生產(chǎn)曲線,分析Eagle Ford頁巖氣藏產(chǎn)量數(shù)據(jù),計算壓裂區(qū)范圍內(nèi)的游離氣總儲量,從而得到SRV[10].這些研究只是靜態(tài)地利用SRV表征儲層的動用情況,并沒有動態(tài)給出儲層動用范圍隨時間的變化,更沒有考慮到開發(fā)后期改造區(qū)以外的儲層動用.另外,雖然可以采用Eclipse等數(shù)模軟件模擬儲層的動用變化,但在模型網(wǎng)格數(shù)較多的情況下,一般數(shù)模軟件需要花費幾十分鐘甚至幾個小時才能完成模擬運算.Jiang Xie等提出利用波前快速推進法(簡稱FMM)預測儲層動用范圍的方法,即使針對百萬計網(wǎng)格的計算也只需要十幾分鐘時間[11-12].在利用FMM進行計算時,由于插值方法所限,所得結(jié)果的精度還有待提高.

    筆者改變FMM的插值計算方法,利用高精度波前快速推進法(簡稱HAFMM)模擬計算波前運移情況,實現(xiàn)對儲層動用范圍隨時間變化更為準確的預測,并且對比不同裂縫形態(tài)下的儲層動用結(jié)果,對頁巖氣藏的生產(chǎn)開發(fā)具有指導作用.

    1 波前運移程函方程

    超低滲是頁巖氣藏的一個重要特點,與常規(guī)氣藏相比,壓力波不可能在短時間內(nèi)傳到整個儲層,相反,在整個氣藏的開發(fā)過程中,儲層的壓力傳播往往限制在距離井筒不遠的范圍內(nèi).對于氣藏,壓力波傳播的范圍即代表儲層動用的范圍.為了對頁巖氣井產(chǎn)能及氣藏開發(fā)狀況進行準確快速的預測,對波前運移的研究顯得更為重要.文中對于壓力波前緣的預測主要基于調(diào)查半徑,Lee W J提出調(diào)查半徑的概念,并給出調(diào)查半徑r的計算公式[13]:

    式中:K為滲透率;t為時間;φ為孔隙度;μ為氣體黏度;ct為綜合壓縮系數(shù).通常,對于不同的流動模型,壓力波傳播的距離和時間關(guān)系可以寫為,其中η=K/(μφct),β是與流態(tài)有關(guān)的常數(shù),在線性流、平面徑向流及球面徑向流中取值分別為2、4、6[14].波前運移的程函方程可以寫為

    其中,未知量τ(x)為飛行時間,它與波的傳播時間t有關(guān),即

    在將飛行時間換算成真實時間時,利用Van Kruysdijk等提出的多段壓裂水平井的復合線性流模型[15](見圖1),將區(qū)域分為4個流態(tài):(1)裂縫線性流;(2)雙線性流;(3)復合線性流;(4)擬徑向流.在換算過程中β在不同區(qū)域分別選取相應的值.

    基質(zhì)和裂縫的程函方程可以分別寫為

    圖1 復合線性流模型Fig.1 Compound linear flow model

    式中:下標m、f分別代表基質(zhì)與裂縫系統(tǒng),其中ηm=Kami/(μgiφmctmi),ηf=Kf/(μgiφfctfi).

    Kami是在考慮滑脫、擴散情況下基質(zhì)的原始表觀滲透率[16],即

    ctmi是考慮解吸條件下基質(zhì)的原始綜合壓縮系數(shù),即

    ctfi是裂縫的原始綜合壓縮系數(shù),即

    式(6-8)中:r為基質(zhì)孔隙半徑;ρi為氣體原始密度;M為儲層氣體的相對分子質(zhì)量;R為通用氣體常數(shù);T為儲層溫度;pi為儲層原始壓力;α為切向動量供給因數(shù),取為0~1;cmi為基質(zhì)系統(tǒng)原始壓縮系數(shù);cfi為裂縫系統(tǒng)原始壓縮系數(shù);cgi為氣體原始壓縮系數(shù);cdi為基質(zhì)原始解吸壓縮系數(shù).

    2 利用HAFMM求解程函方程

    2.1 計算方法

    HAFMM是Sethian J A提出的可以高效計算程函方程的一種方法[17],利用HAFMM求解計算波前運移的程函方程,可以對儲層動用進行預測.利用HAFMM計算時用到的點見圖2.

    假設 點τi,j未 知,它 周 圍 的 8 個 點τi,j+1、τi,j-1、τi+1,j、τi-1,j、τi,j+2、τi,j-2、τi+2,j、τi-2,j是已知的,Sethian J A 提出計算τi,j方法,則

    由式(2)可得

    圖2 波前快速推進法中的相鄰點Fig.2 Neighborhood points in FMM

    將式(10)代入式(9),可以得到η,即

    式(11)中,差分運算D的計算方法[18]為

    式(16)是一個簡單的一元二次方程,對其求解,選出較大的解作為τi,j值.當τi,j周圍的8個點中,如果有任意點值不為已知,在計算τi,j時可不將其考慮在內(nèi).

    2.2 計算過程

    HAFMM的計算過程見圖3:

    (1)找到所有的已知點標記為初始點,見圖3(a)中的黑點;(2)找到初始點周圍的相鄰點,并利用初始點計算相鄰點的值,標記為備選點,見圖3(b)中的藍點;(3)選出所有備選點中的最小值,標記為已知點,見圖3(c)紅框中的點,并以它為初始點計算相鄰點的值,標記為備選點,見圖3(d);(4)重復步驟(2)和(3),直到所有點被標記為已知點.

    式(12-15)中:Δx、Δy分別代表x、y方向上兩點間的距離.在x方向上,比較、、0的大小,選出最大值,確保壓力波沿x方向一直向外推進,假設為τ,同理也選出y方向的最大值,假設為τ,式(11)可以改寫為

    圖3 HAFMM計算過程Fig.3 Calculation process of HAFMM

    3 結(jié)果分析

    3.1 參數(shù)選取

    綜合文獻[16,19],在計算頁巖氣藏壓力波前緣運移情況時,選用的參數(shù)見表1.

    在計算基質(zhì)表觀滲透率時,考慮基質(zhì)孔隙的吸附作用,在基質(zhì)的孔隙壁上吸附一層分子膜,所以在計算時,需要使用基質(zhì)的有效孔隙半徑re,re=rm-2rg,其中,rm為基質(zhì)孔隙半徑;rg為氣體半徑.

    將re代入到式(6)中,可以得到基質(zhì)的表觀滲透率Kami=5.21×10-6μm2.

    3.2 HAFMM與FMM結(jié)果對比

    與傳統(tǒng)的 FMM 不同,HAFMM 改變 FMM 中的一階近似插值方法(τi+1,j-τi,j)/Δx,替換成二階近似,使得所計算的結(jié)果更加準確.為了比較 HAFMM 與FMM 的計算結(jié)果,選用裂縫模型(見圖4),其中中水平井長度L=1 000 m,裂縫長度Lf=500 m,裂縫間距為200 m.

    表1 儲層參數(shù)Table 1 Reservoir parameters

    圖4 裂縫模型及HAFMM與FMM計算的動用范圍隨時間變化結(jié)果Fig.4 Fracture model and drainage area vs time under FMM and HAFMM

    儲層動用范圍與時間的關(guān)系曲線見圖5(其中箭頭所指處代表裂縫間的壓力達到連通).由圖4及圖5可以看出,利用HAFMM計算時,在相同的時間內(nèi)儲層動用的范圍比利用FMM計算傳播的范圍更大.以氣井生產(chǎn)3 a為例,利用HAFMM計算時儲層的動用范圍為0.823 km2,比利用FMM的計算結(jié)果0.779 km2要多出0.044 km2;當氣井的生產(chǎn)時間達到12 a時,前者的計算結(jié)果為1.402 km2,后者的計算結(jié)果為1.33 km2,兩者相差0.072 km2,且隨著生產(chǎn)時間的增加,兩者的差別越來越大.

    3.3 不同裂縫形態(tài)下的儲層動用

    在頁巖氣藏開發(fā)的過程中,裂縫的形態(tài)對氣井產(chǎn)能具有直接影響.在裂縫總長度相等的條件下,考慮縫間干擾對儲層動用的影響及工程施工的可行性,選取5種裂縫形態(tài)進行對比(見表2).其中裂縫間距為200 m,并且假設氣井生產(chǎn)過程中未受到邊界的影響.

    經(jīng)HAFMM計算的結(jié)果見圖6.圖6的橙色、紫色和粉色區(qū)域分別代表氣井生產(chǎn)30 d、2 a和9 a時儲層的動用范圍.壓力波的傳播范圍與時間的關(guān)系曲線見圖7.

    由圖7可以看出,當氣井生產(chǎn)570 d時,裂縫之間的壓力達到連通,5種裂縫形態(tài)下的儲層動用范圍為0.681 km2,在此之前儲層動用范圍與裂縫形態(tài)無關(guān).這是因為當裂縫之間的壓力達到連通之前,各裂縫之間不存在干擾.當裂縫之間的壓力達到連通后,不同的裂縫形態(tài)下儲層動用具有明顯區(qū)別.在表2(d)時儲層動用的最快,表2(e)的動用速度要低于表2(d)的,但明顯高于其他3種情況,表2(b)和(c)的儲層動用較慢,表2(a)的最慢.這是因為在壓力連通后,裂縫之間正對的區(qū)域形成嚴重的干擾區(qū),阻礙壓力波的運移.由于在表2(a)時各裂縫長度相等,縫間正對面積達到0.4 km2,干擾區(qū)面積最大,所以儲層動用最慢.表2(b)和(c)分別采用紡錘型和啞鈴型裂縫形態(tài),長短裂縫依次排列,縫間正對面積為0.32 km2,一定程度上減少干擾區(qū)的面積,所以它們的動用速度稍快.表2(d)和(e)采用長裂縫與短裂縫交錯排列的辦法,可有效地將縫間正對面積分別降到0.16 km2和0.24 km2,所以表2(d)的動用最快,表2(e)的次之.

    圖5 不同算法的儲層動用與時間關(guān)系曲線Fig.5 Curve of drainage area vs time under different calculation method

    表2 裂縫形態(tài)與長度Table 2 Fracture pattern and length

    圖6 不同裂縫形態(tài)下的氣井生產(chǎn)30 d、2 a及9 a時的儲層動用Fig.6 Drainage area after 30 d,2 years and 9 years under different fracture patterns

    圖7 不同裂縫形態(tài)的儲層動用與時間關(guān)系曲線Fig.7 Curve of drainage area vs time under different fracture patterns

    4 結(jié)論

    (1)利用波前快速推進法對頁巖氣儲層的波前運移進行預測,并且給出儲層動用情況與時間之間的定量關(guān)系.

    (2)對比波前快速推進法與高精度波前快速推進法的結(jié)果,利用高精度波前快速推進法計算時,在相同時間內(nèi),儲層動用范圍更大.

    (3)選取5種不同的裂縫形態(tài)進行對比,在裂縫之間的壓力達到連通前,儲層的動用與裂縫形態(tài)無關(guān);當裂縫之間的壓力達到連通后,不同的裂縫形態(tài)對儲層動用具有明顯的影響.

    (4)對比不同裂縫形態(tài)下的儲層動用情況,縫間正對面積越小,即縫間干擾越小,儲層動用的越快.

    [1] Nobakht M,Clarkson C R.Hybrid forecasting methods for multi-fractured horizontal wells:EUR sensitivities[C].SPE 153220,2012.

    [2] Cipolla C L,Lolon E,Mayerhofer M J.Reservoir modeling and production evaluation in shale-gas reservoirs[C].SPE 13185,2009.

    [3] 孫贊東,賈承造,李相方,等.非常規(guī)油氣勘探與開發(fā)[M].北京:石油工業(yè)出版社,2011:865-1116.

    Sun Zandong,Jia Chengzao,Li Xiangfang,et al.Unconventional oil and gas exploration and development[M].Beijing:Petroleum Industry Press,2011:865-1116.

    [4] 單衍勝,張金川,李曉光,等.遼河盆地東部凸起太原組頁巖氣聚集條件及有力預測[J].大慶石油學院學報,2012,36(1):1-7.

    Shan Yansheng,Zhang Jinchuan,Li Xiaoguang,et al.Shale gas accumulation factors and prediction of favorable area Taiyuan formation in Liaohe eastern uplift[J].Journal of Daqing Petroleum Institute,2012,36(1):1-7.

    [5] 梁超,姜在興,郭嶺,等.陸棚相黑色泥巖發(fā)育特征、沉積演化及頁巖氣勘探前景——以甕安永和剖面牛蹄塘組為例[J].大慶石油學院學報,2011,35(6):13-21.

    Liang Chao,Jiang Zaixing,Guo Ling,et al.Characteristics of black shale,sedimentary evolution and shale gas exploration prospect of shale face taking Weng'an Yonghe profile Niutitang group as an example[J].Journal of Daqing Petroleum Institute,2011,35(6):13-21.

    [6] 久凱,丁文龍,黃文輝,等.渤灣地盆地濟陽拗陷沙三段頁巖氣地質(zhì)條件分析[J].大慶石油學院學報,2012,36(2):65-70.

    Jiu Kai,Ding Wenlong,Huang Wenhui,et al.Analysis of geological condition of the formation of Shahejie shale gas in Jiyang depression in Bohai bay basin[J].Journal of Daqing Petroleum Institute,2012,36(2):65-70.

    [7] 任俊杰,郭平,王德龍,等.頁巖氣藏壓裂水平井產(chǎn)能模型及影響因素[J].東北石油大學學報,2012,36(6):76-81.

    Ren Junjie,Guo Ping,Wang Delong,et al.Productivity model of fractural horizontal wells in shale gas reservoir and analysis of influential factors[J].Journal of Northeast Petroleum University,2012,36(6):76-81.

    [8] Fisher M K,Heinze J R.Optimizing horizontal completion techniques in the Barnett shale using microseismic fracture mapping[C].SPE 2004,26-29.

    [9] Mayerhofer M J,Lolon E P,Warpinski N R,et al.What is stimulated reservoir volume[C].SPE 119890,2008.

    [10] 徐兵祥,李相方.頁巖氣產(chǎn)量數(shù)據(jù)分析方法及產(chǎn)能預測[J].中國石油大學學報:自然科學版,2011,37(3):120-125.

    Xu Bingxiang,Li Xiangfang.Production data analysis and productivity forecast of shale gas reservoir[J]. Journal of China University of Petroleum:Natural Science Edition,2011,37(3):120-125.

    [11] Zhang Yanbin,Yang Changdong,King Michael J,et al.Fast marching methods for complex grids and anisotropic permeabilities:Application to unconventional reservoirs[C].SPE 163637,2013.

    [12] Jiang Xie,Neha Gupta,Michael J King,et al.Depth of investigation and depletion behavior in unconventional reservoirs using fast marching methods[C].SPE 154532,2012.

    [13] Lee W J.Well Testing[M]∥Richardson.SPE Textbook Series.Texas:Society of Petroleum Engineers,1982.

    [14] Kim J U,Data-Gupta A,Brouwer R,et al.Calibration of high-resolution reservoir models using transient pressure data[C].SPE 124834,2009.

    [15] Van Kruysdijk,Dullaert G M.A boundary element solution of the transient pressure response of multiple fractured horizontal wells[C].Cambridge:Paper Presented at The 2nd European Conference of The Mathematicsof Oil Recovery,1989.

    [16] Shabro V,Torres-Verdin C,Sepehrnoori K.Forecasting gas production in organic shale with the combined numerical simulation of gas diffusion in kerogen,Langmuir desorption from kerogen surfaces,and advection in nanopores[C].SPE 159250,2012.

    [17] Sethian J A.Level set methods and fast marching methods[M].Cambridge:Cambridge University Press,1999:86-100.

    [18] Andreas Barentaen J.On the implementation of fast marching methods for 3d lattices[C].Department of Mathematical Modeling,2001.

    [19] 王瑞,張寧生,劉曉娟,等.考慮吸附和擴散的頁巖是滲透率及其與溫度、壓力之關(guān)系[J].西安石油大學學報:自然科學版,2013,28(2):50-53.

    Wang Rui,Zhang Ningsheng,Liu Xiaojuan,et al.Apparent permeability of shale considering the adsorption and diffusion of gas and the effects of temperature and pressure on it[J].Journal of Xi'an Shiyou University:Natural Science Edition,2013,28(2):50-53.

    TE332

    A

    2095-4107(2014)04-0017-06

    2013-10-12;

    陸雅玲

    國家“973”重點基礎研究發(fā)展計劃項目(SQ2012CB027805)

    黃世軍(1974-),男,博士,副教授,主要從事油氣滲流機理方面的研究.

    DOI 10.3969/j.issn.2095-4107.2014.04.003

    猜你喜歡
    動用運移氣藏
    滿語動詞的使動用法
    滿族文學(2023年5期)2023-09-12 03:16:26
    曲流河復合點壩砂體構(gòu)型表征及流體運移機理
    東營凹陷北帶中淺層油氣運移通道組合類型及成藏作用
    基于波前快速推進法的頁巖氣儲層動用預測
    開采過程中上覆急傾斜巖層運移規(guī)律模擬與研究
    煤炭學報(2015年10期)2015-12-21 01:55:49
    川西坳陷孝泉-新場地區(qū)陸相天然氣地球化學及運移特征
    致密氣藏壓裂傾斜縫壓力動態(tài)分析
    汽車起動用蓄電池健康診斷及應急點火儀的研制
    電測與儀表(2015年5期)2015-04-09 11:31:04
    塔里木油田超深超高壓氣藏的成功改造
    斷塊油氣田(2014年6期)2014-03-11 15:34:03
    LG地區(qū)復雜礁灘氣藏描述技術(shù)
    香蕉av资源在线| 久久久久久久久大av| 国产精品久久久久久亚洲av鲁大| 国产男靠女视频免费网站| 18禁黄网站禁片免费观看直播| 宅男免费午夜| 国产一区二区在线av高清观看| 99国产综合亚洲精品| 中文字幕熟女人妻在线| 免费大片18禁| 久久久久国内视频| 国产成人av教育| 久久精品国产亚洲av天美| 日日摸夜夜添夜夜添小说| 欧美3d第一页| 国产高清激情床上av| 十八禁网站免费在线| 变态另类丝袜制服| 国产精品综合久久久久久久免费| 亚洲人成伊人成综合网2020| 免费黄网站久久成人精品 | 中文在线观看免费www的网站| 亚洲中文日韩欧美视频| 成年版毛片免费区| 国产乱人视频| 免费av毛片视频| 一夜夜www| 欧美成人a在线观看| 美女cb高潮喷水在线观看| 嫩草影院入口| 国产真实乱freesex| 久久草成人影院| 一本精品99久久精品77| 欧美日本视频| 精品久久久久久久人妻蜜臀av| 成人国产综合亚洲| 国产精品综合久久久久久久免费| 少妇被粗大猛烈的视频| 哪里可以看免费的av片| 天美传媒精品一区二区| 一边摸一边抽搐一进一小说| 午夜精品久久久久久毛片777| 亚洲第一电影网av| 婷婷丁香在线五月| 久久久精品大字幕| 国产熟女xx| 少妇高潮的动态图| 成人特级黄色片久久久久久久| 一进一出抽搐gif免费好疼| 3wmmmm亚洲av在线观看| 我的女老师完整版在线观看| 色哟哟·www| 观看美女的网站| 亚洲国产色片| 老司机午夜福利在线观看视频| 两人在一起打扑克的视频| 丁香六月欧美| 国产精品国产高清国产av| 国产午夜福利久久久久久| 精品久久久久久,| 亚洲自偷自拍三级| 制服丝袜大香蕉在线| 能在线免费观看的黄片| 婷婷丁香在线五月| 怎么达到女性高潮| 成人国产综合亚洲| 国产高清视频在线播放一区| av欧美777| 久久精品影院6| 在线观看美女被高潮喷水网站 | 三级毛片av免费| 99热这里只有精品一区| 国产欧美日韩精品一区二区| 欧美日韩国产亚洲二区| 亚洲色图av天堂| 身体一侧抽搐| 亚洲自拍偷在线| 88av欧美| 亚洲aⅴ乱码一区二区在线播放| 免费在线观看亚洲国产| 国产亚洲欧美98| 又黄又爽又刺激的免费视频.| 国产伦在线观看视频一区| 国产精品不卡视频一区二区 | 毛片一级片免费看久久久久 | 十八禁网站免费在线| 国产黄色小视频在线观看| 一个人观看的视频www高清免费观看| 国产单亲对白刺激| 国产精品永久免费网站| 看片在线看免费视频| 欧美日本视频| 在线a可以看的网站| 男人舔奶头视频| 亚洲乱码一区二区免费版| 搡老岳熟女国产| 五月玫瑰六月丁香| 亚洲av熟女| 一区二区三区四区激情视频 | 一个人免费在线观看的高清视频| 欧美性猛交╳xxx乱大交人| 欧美在线一区亚洲| 国产精品亚洲av一区麻豆| 欧美另类亚洲清纯唯美| av专区在线播放| 又紧又爽又黄一区二区| 欧美一区二区精品小视频在线| 久久这里只有精品中国| АⅤ资源中文在线天堂| 十八禁人妻一区二区| 免费看日本二区| 国产精品电影一区二区三区| 国产黄片美女视频| 亚洲成人久久性| 国产野战对白在线观看| 99久久精品热视频| 欧美精品国产亚洲| 少妇裸体淫交视频免费看高清| 性色avwww在线观看| 老司机午夜十八禁免费视频| 国产精品一区二区免费欧美| 少妇熟女aⅴ在线视频| 91字幕亚洲| 亚州av有码| 性插视频无遮挡在线免费观看| 欧美色欧美亚洲另类二区| 成人一区二区视频在线观看| 中文亚洲av片在线观看爽| 久久精品国产亚洲av香蕉五月| 天天躁日日操中文字幕| 十八禁国产超污无遮挡网站| 国产av在哪里看| 久久久久久久精品吃奶| 窝窝影院91人妻| 欧美日本亚洲视频在线播放| 色精品久久人妻99蜜桃| 亚洲中文日韩欧美视频| 日韩中字成人| 97热精品久久久久久| 亚洲在线观看片| 精品人妻偷拍中文字幕| 日本免费a在线| 亚洲男人的天堂狠狠| 亚洲电影在线观看av| 亚洲av日韩精品久久久久久密| 最新在线观看一区二区三区| 夜夜看夜夜爽夜夜摸| 久久久久久国产a免费观看| 国产精品亚洲av一区麻豆| 国产精品1区2区在线观看.| 日本精品一区二区三区蜜桃| 少妇人妻精品综合一区二区 | 欧美性猛交黑人性爽| 国产成人欧美在线观看| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站| 色综合欧美亚洲国产小说| 夜夜看夜夜爽夜夜摸| 亚洲av日韩精品久久久久久密| 中出人妻视频一区二区| 成人毛片a级毛片在线播放| 超碰av人人做人人爽久久| 日本三级黄在线观看| 欧美日韩黄片免| 极品教师在线免费播放| 亚洲天堂国产精品一区在线| 我的老师免费观看完整版| 极品教师在线免费播放| 亚洲在线自拍视频| eeuss影院久久| 高清毛片免费观看视频网站| 国产精品亚洲美女久久久| 国产精品av视频在线免费观看| 日韩亚洲欧美综合| 国产亚洲精品av在线| 国产精品av视频在线免费观看| 麻豆国产av国片精品| 久久久久免费精品人妻一区二区| 97超视频在线观看视频| 成年女人毛片免费观看观看9| 一级黄色大片毛片| 淫秽高清视频在线观看| 国产aⅴ精品一区二区三区波| 又粗又爽又猛毛片免费看| 日日干狠狠操夜夜爽| 国产aⅴ精品一区二区三区波| 日韩欧美国产一区二区入口| 精品久久久久久久末码| 欧美+日韩+精品| 色吧在线观看| 婷婷丁香在线五月| 国内少妇人妻偷人精品xxx网站| 国产高清有码在线观看视频| 国产成人a区在线观看| 亚洲精品一区av在线观看| 色精品久久人妻99蜜桃| 国产免费男女视频| 琪琪午夜伦伦电影理论片6080| 色综合欧美亚洲国产小说| 身体一侧抽搐| 久久久久久久精品吃奶| 中文资源天堂在线| 宅男免费午夜| 精品久久久久久成人av| 午夜影院日韩av| 禁无遮挡网站| 精品福利观看| 国产一区二区在线av高清观看| 午夜福利成人在线免费观看| 悠悠久久av| 91在线精品国自产拍蜜月| 淫妇啪啪啪对白视频| 美女大奶头视频| 久久九九热精品免费| 国产精品亚洲美女久久久| 99热只有精品国产| 欧美另类亚洲清纯唯美| 久久人人精品亚洲av| 日本一二三区视频观看| 国产精品国产高清国产av| 亚洲精品粉嫩美女一区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 18禁裸乳无遮挡免费网站照片| 亚洲中文日韩欧美视频| 午夜福利18| 中亚洲国语对白在线视频| 99在线人妻在线中文字幕| 国产真实伦视频高清在线观看 | 亚洲狠狠婷婷综合久久图片| 亚洲av电影不卡..在线观看| 国产精品1区2区在线观看.| 天堂动漫精品| 国产中年淑女户外野战色| 香蕉av资源在线| 午夜亚洲福利在线播放| 色综合欧美亚洲国产小说| 日韩精品中文字幕看吧| 99久久精品热视频| 日本一二三区视频观看| 欧美高清成人免费视频www| 男女下面进入的视频免费午夜| 中国美女看黄片| av专区在线播放| 精品久久久久久久久亚洲 | 91久久精品国产一区二区成人| 狠狠狠狠99中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 变态另类成人亚洲欧美熟女| 精品午夜福利在线看| 欧美黑人巨大hd| 亚洲欧美激情综合另类| 久久久久久久久久成人| 国产av在哪里看| 午夜亚洲福利在线播放| 亚洲不卡免费看| 九九在线视频观看精品| 九九久久精品国产亚洲av麻豆| 岛国在线免费视频观看| 亚洲人成网站在线播| 亚洲成人免费电影在线观看| 国产亚洲欧美在线一区二区| 欧美3d第一页| 变态另类成人亚洲欧美熟女| 亚洲精品色激情综合| 嫁个100分男人电影在线观看| 亚洲人成网站高清观看| 色哟哟哟哟哟哟| 国产黄a三级三级三级人| 亚洲人成伊人成综合网2020| 内射极品少妇av片p| 久久精品夜夜夜夜夜久久蜜豆| 一区二区三区免费毛片| 精品久久久久久久久久免费视频| 两个人的视频大全免费| 丰满乱子伦码专区| 美女高潮的动态| 少妇高潮的动态图| 宅男免费午夜| 精品熟女少妇八av免费久了| 国产伦在线观看视频一区| 久久亚洲精品不卡| 婷婷丁香在线五月| 久久国产乱子免费精品| 成年免费大片在线观看| 亚洲真实伦在线观看| 亚洲欧美日韩高清专用| 久久精品91蜜桃| 男女视频在线观看网站免费| 黄色配什么色好看| 嫩草影院新地址| 国产在线男女| 麻豆国产97在线/欧美| 午夜久久久久精精品| 五月伊人婷婷丁香| 美女高潮的动态| 91字幕亚洲| 久久久色成人| 在线观看免费视频日本深夜| 一个人观看的视频www高清免费观看| 国产精品亚洲美女久久久| 久久久久久久午夜电影| 成人美女网站在线观看视频| 色综合亚洲欧美另类图片| 如何舔出高潮| 五月伊人婷婷丁香| 人妻久久中文字幕网| 人人妻人人澡欧美一区二区| 国产精品嫩草影院av在线观看 | 日本黄大片高清| 性色avwww在线观看| 看十八女毛片水多多多| 男女视频在线观看网站免费| 一个人看的www免费观看视频| 嫩草影院新地址| 欧美日韩国产亚洲二区| 一级毛片久久久久久久久女| 国产伦在线观看视频一区| 成人特级av手机在线观看| 看免费av毛片| 久久精品影院6| 久久久久久久亚洲中文字幕 | 亚洲中文日韩欧美视频| 欧美日韩中文字幕国产精品一区二区三区| 脱女人内裤的视频| www.999成人在线观看| 国产大屁股一区二区在线视频| 国产主播在线观看一区二区| 亚洲不卡免费看| 久久精品国产自在天天线| 男人狂女人下面高潮的视频| 五月伊人婷婷丁香| 精品久久久久久成人av| 久久这里只有精品中国| 久久伊人香网站| www日本黄色视频网| 成人鲁丝片一二三区免费| 男人舔奶头视频| 18禁黄网站禁片午夜丰满| 国产高潮美女av| 亚洲av免费高清在线观看| 成人特级黄色片久久久久久久| 国模一区二区三区四区视频| 熟女电影av网| 亚洲成av人片免费观看| 永久网站在线| 国产精华一区二区三区| 三级国产精品欧美在线观看| 亚洲人成电影免费在线| 97碰自拍视频| 国产亚洲欧美98| 亚洲真实伦在线观看| 午夜激情福利司机影院| 人人妻人人澡欧美一区二区| 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 国产视频内射| 精品人妻视频免费看| 日本黄色片子视频| 亚洲精品成人久久久久久| 蜜桃亚洲精品一区二区三区| 午夜a级毛片| 欧美成人免费av一区二区三区| 91av网一区二区| 国产精品伦人一区二区| 亚洲精品成人久久久久久| 亚洲第一欧美日韩一区二区三区| 国产真实伦视频高清在线观看 | 免费看光身美女| 国产伦人伦偷精品视频| 久久久久国产精品人妻aⅴ院| 岛国在线免费视频观看| 国产精品1区2区在线观看.| 日本 欧美在线| av天堂中文字幕网| 久久性视频一级片| 一进一出抽搐gif免费好疼| 成人性生交大片免费视频hd| 悠悠久久av| 精品熟女少妇八av免费久了| 久久久久亚洲av毛片大全| 国产成人av教育| 国产精品影院久久| 久久草成人影院| 日韩中字成人| 国产av一区在线观看免费| 国产免费av片在线观看野外av| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久精品电影| 久久久久精品国产欧美久久久| 久久99热这里只有精品18| 天堂影院成人在线观看| 一二三四社区在线视频社区8| 国产探花在线观看一区二区| 露出奶头的视频| 亚洲午夜理论影院| 一个人免费在线观看电影| 人妻制服诱惑在线中文字幕| 18美女黄网站色大片免费观看| 国产视频一区二区在线看| 桃红色精品国产亚洲av| ponron亚洲| 国产中年淑女户外野战色| 国产男靠女视频免费网站| 欧美3d第一页| 国产精品av视频在线免费观看| 国产av一区在线观看免费| 性色avwww在线观看| 亚洲国产高清在线一区二区三| 一本一本综合久久| 此物有八面人人有两片| 蜜桃亚洲精品一区二区三区| 99久久精品热视频| 高潮久久久久久久久久久不卡| 国产一区二区在线av高清观看| 小蜜桃在线观看免费完整版高清| 在线观看美女被高潮喷水网站 | 亚洲精品粉嫩美女一区| 在线观看舔阴道视频| 我要搜黄色片| .国产精品久久| 午夜福利视频1000在线观看| 午夜福利高清视频| 日本黄色片子视频| 久久亚洲精品不卡| 天堂网av新在线| 波多野结衣高清作品| 99久久成人亚洲精品观看| 91av网一区二区| 97超视频在线观看视频| 一区二区三区免费毛片| 免费av观看视频| 国产亚洲欧美在线一区二区| 午夜精品一区二区三区免费看| 亚洲精品色激情综合| 欧美xxxx性猛交bbbb| 国语自产精品视频在线第100页| 黄色视频,在线免费观看| 国产69精品久久久久777片| 国产男靠女视频免费网站| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区三区四区免费观看 | 国产欧美日韩精品亚洲av| 国产成人啪精品午夜网站| 2021天堂中文幕一二区在线观| 久久久国产成人免费| 中文亚洲av片在线观看爽| 啪啪无遮挡十八禁网站| 国产伦在线观看视频一区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲av成人精品一区久久| 他把我摸到了高潮在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品乱码一区二三区的特点| 久久国产乱子伦精品免费另类| 一区二区三区免费毛片| 婷婷亚洲欧美| 国产一区二区三区在线臀色熟女| 啦啦啦观看免费观看视频高清| 色噜噜av男人的天堂激情| 麻豆av噜噜一区二区三区| 亚洲真实伦在线观看| 51国产日韩欧美| 精品一区二区三区视频在线观看免费| 亚洲欧美激情综合另类| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 97人妻精品一区二区三区麻豆| 少妇人妻精品综合一区二区 | 每晚都被弄得嗷嗷叫到高潮| 高清毛片免费观看视频网站| 精品福利观看| 亚洲七黄色美女视频| 国产亚洲精品av在线| 有码 亚洲区| 在线播放无遮挡| 亚洲黑人精品在线| 身体一侧抽搐| 两人在一起打扑克的视频| 热99在线观看视频| 成人欧美大片| 一区二区三区四区激情视频 | 桃色一区二区三区在线观看| 日本 av在线| 免费在线观看日本一区| 制服丝袜大香蕉在线| 女同久久另类99精品国产91| 亚洲国产色片| 欧美在线黄色| 国产色爽女视频免费观看| 久久久久久久午夜电影| 少妇高潮的动态图| 一级作爱视频免费观看| av国产免费在线观看| 国产又黄又爽又无遮挡在线| 女生性感内裤真人,穿戴方法视频| 老司机深夜福利视频在线观看| 久久久久国内视频| 国产精品影院久久| 高清在线国产一区| 亚洲无线在线观看| 婷婷亚洲欧美| 精品熟女少妇八av免费久了| 久久精品国产自在天天线| 久久久精品大字幕| 精品久久国产蜜桃| 丁香六月欧美| 丁香欧美五月| 国产成人福利小说| 国产大屁股一区二区在线视频| 色吧在线观看| 午夜久久久久精精品| www.色视频.com| 校园春色视频在线观看| 村上凉子中文字幕在线| 99久国产av精品| 精品人妻视频免费看| 成人无遮挡网站| 国产亚洲精品久久久久久毛片| 欧美成人一区二区免费高清观看| 看十八女毛片水多多多| 99国产精品一区二区蜜桃av| 两个人的视频大全免费| 色av中文字幕| 床上黄色一级片| 757午夜福利合集在线观看| 亚洲欧美日韩高清专用| 日韩中字成人| 国产精品野战在线观看| 免费搜索国产男女视频| 欧美+亚洲+日韩+国产| 精品99又大又爽又粗少妇毛片 | 丝袜美腿在线中文| 日韩免费av在线播放| 国产精品久久电影中文字幕| 色综合站精品国产| 看十八女毛片水多多多| 国产成人av教育| 久久99热这里只有精品18| 午夜免费男女啪啪视频观看 | 亚洲aⅴ乱码一区二区在线播放| 我的老师免费观看完整版| 精品日产1卡2卡| 精品久久久久久久久av| 国语自产精品视频在线第100页| 欧美xxxx性猛交bbbb| 国产蜜桃级精品一区二区三区| 亚洲中文日韩欧美视频| 日韩高清综合在线| 亚洲成人精品中文字幕电影| 国产色婷婷99| 免费人成视频x8x8入口观看| 嫩草影视91久久| 级片在线观看| 少妇裸体淫交视频免费看高清| 国产一级毛片七仙女欲春2| 性欧美人与动物交配| 亚洲av美国av| 久久精品国产亚洲av涩爱 | 亚洲精品日韩av片在线观看| 中文字幕久久专区| 亚洲欧美日韩卡通动漫| 成人国产一区最新在线观看| 亚洲第一电影网av| 偷拍熟女少妇极品色| 精品一区二区三区视频在线观看免费| 99久国产av精品| 国产一区二区激情短视频| 国产乱人伦免费视频| 免费在线观看亚洲国产| 亚洲av美国av| 国产毛片a区久久久久| 日韩有码中文字幕| 亚洲美女视频黄频| 搡老熟女国产l中国老女人| 日本一二三区视频观看| 精品久久久久久久久亚洲 | 国产男靠女视频免费网站| av在线老鸭窝| 成人午夜高清在线视频| 大型黄色视频在线免费观看| 中文在线观看免费www的网站| 99久久九九国产精品国产免费| 熟女人妻精品中文字幕| 男女做爰动态图高潮gif福利片| 在线观看舔阴道视频| 久久久色成人| 婷婷六月久久综合丁香| 中文资源天堂在线| 亚洲欧美日韩高清专用| 婷婷六月久久综合丁香| 国产精品亚洲美女久久久| 亚洲国产色片| 欧美午夜高清在线| 国产免费男女视频| 18禁裸乳无遮挡免费网站照片| 搡老岳熟女国产| 国产在线男女| 国产亚洲精品久久久com| 最新中文字幕久久久久| 韩国av一区二区三区四区| 麻豆国产av国片精品| 日韩欧美一区二区三区在线观看| 欧美乱色亚洲激情| 熟女人妻精品中文字幕| 亚洲,欧美精品.| 日韩欧美在线乱码| 久久久久久久精品吃奶| 精品一区二区三区人妻视频| x7x7x7水蜜桃| 精品久久久久久久久亚洲 | 国产在线精品亚洲第一网站| 女生性感内裤真人,穿戴方法视频| 一本综合久久免费| 色精品久久人妻99蜜桃| 老熟妇仑乱视频hdxx| 国产综合懂色|