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

    頁巖氣藏滲流機理及壓力動態(tài)分析

    2015-11-30 11:52:58杜殿發(fā)王妍妍付金剛孫召勃中國石油大學(xué)華東石油工程學(xué)院山東青島266580
    計算物理 2015年1期
    關(guān)鍵詞:氣藏機理滲流

    杜殿發(fā),王妍妍,付金剛,孫召勃,喬 妮(中國石油大學(xué)(華東)石油工程學(xué)院,山東青島 266580)

    文章編號:1001?246X(2015)01?0051?07

    頁巖氣藏滲流機理及壓力動態(tài)分析

    杜殿發(fā),王妍妍,付金剛,孫召勃,喬 妮
    (中國石油大學(xué)(華東)石油工程學(xué)院,山東青島 266580)

    針對頁巖氣藏的特點,分析國內(nèi)外對解吸、擴散等數(shù)學(xué)描述方法,建立三種機理作用下的滲流數(shù)學(xué)模型:①考慮有機基質(zhì)表面解吸,解吸氣和游離氣共同竄流的數(shù)學(xué)模型;②考慮有機基質(zhì)解吸和擴散的數(shù)學(xué)模型;③考慮有機基質(zhì)中的氣體解吸、擴散,以及無機基質(zhì)中氣體竄流綜合作用的數(shù)學(xué)模型.運用點源函數(shù)方法,得到模型的基本點源解,對點源解作積分變換得到垂直壓裂井、壓裂水平井兩種開采方式下的地層壓力解,繪制出三種機理影響下的壓力動態(tài)曲線,并從滲流機理上分析曲線特征以及不同曲線存在差異的原因.

    頁巖氣藏;滲流機理;滲流模型;試井分析

    0 引言

    目前對頁巖氣藏滲流機理的研究已越來越深入,已有的滲流數(shù)學(xué)模型對頁巖氣藏進行壓力動態(tài)分析研究,得到了比較符合頁巖氣藏滲流規(guī)律的結(jié)果[1].但是由于頁巖氣藏滲流介質(zhì)十分復(fù)雜,存在有機基質(zhì)、無機基質(zhì)、天然微裂縫、大的裂縫等多種孔隙介質(zhì),導(dǎo)致滲流機理研究上存在著很大的爭議.一般認為吸附氣在氣藏中占很大比例,但以美國福特沃斯盆地Barnett頁巖氣藏為代表,解吸氣量遠少于游離氣量,因此氣體進入裂縫系統(tǒng)的方式也需要依具體氣藏實際情況而定.另外,即便是吸附氣量較大,解吸作用在使孔隙內(nèi)氣體濃度增大的同時也必然使壓力增大,因此濃度差引發(fā)的擴散和壓力差導(dǎo)致的滲流在頁巖氣流動過程中起的作用也沒有定論.針對這些問題,本文考慮不同的滲流機理,建立與之對應(yīng)的數(shù)學(xué)模型,并在此基礎(chǔ)上進行壓力動態(tài)分析.

    1 滲流機理及數(shù)學(xué)模型

    為簡化求解,對不同滲流機理下模型作如下規(guī)定:①解吸作用符合Langmuir等溫方程,解吸作用僅發(fā)生在有機基質(zhì)表面;②擴散作用符合Fick擬穩(wěn)態(tài)擴散定律;③氣體流動為單相流,滲流過程等溫;④儲層各向同性,外邊界為無限大;⑤忽略毛管力和重力影響.

    1.1 考慮基質(zhì)解吸、竄流機理的數(shù)學(xué)模型

    該模型適用于氣藏開采后期,吸附氣含量較少的氣藏.模型認為頁巖氣藏屬于基質(zhì)-裂縫雙重介質(zhì)油藏,滲流分為以下三個階段:

    1)氣井開井生產(chǎn),地層壓力下降,氣體在有機質(zhì)表面向內(nèi)解吸,進入基質(zhì)孔隙內(nèi);

    2)基質(zhì)孔隙內(nèi)氣體量的增加導(dǎo)致壓力升高,解吸的氣體與游離氣體混合在一起,在壓力差的作用下從基質(zhì)孔隙竄流進入裂縫系統(tǒng);

    3)裂縫中的氣體在流體勢的作用下遵循達西定律流向井筒.

    由于解吸的氣體會進入基質(zhì)孔隙,與基質(zhì)孔隙內(nèi)部的游離氣混合,其作用相當(dāng)于增加了基質(zhì)系統(tǒng)的儲容能力.因此可以通過增加基質(zhì)彈性儲容系數(shù)的方式定量描繪解吸的氣體量.此處認為解吸的氣體量使得基質(zhì)系統(tǒng)等溫壓縮系數(shù)增加[2],此時滲流微分方程為-

    為了降低方程的非線性,許多參數(shù)采用了地層平均壓力下的值,選擇以下公式[3]

    引入?yún)?shù)

    對方程進行化簡,并進行Laplace變換,最終得到:

    1.2 考慮基質(zhì)解吸、基質(zhì)擴散機理的數(shù)學(xué)模型

    該模型適用于有機質(zhì)孔隙均比較發(fā)育、剛投入開發(fā)且吸附氣含量比較多的氣藏.認為頁巖氣藏屬于基質(zhì)-裂縫雙重介質(zhì)油藏,將頁巖氣滲流過程描述為以下3個階段[4]:

    1)氣井開始生產(chǎn),地層壓力下降,頁巖氣在基質(zhì)表面發(fā)生解吸,解吸的氣體進入基質(zhì)孔隙;

    2)基質(zhì)孔隙內(nèi)的氣體濃度大于裂縫系統(tǒng),氣體在濃度差的作用下擴散進入裂縫系統(tǒng);

    3)裂縫中的氣體在流體勢的作用下遵循達西定律流向井筒.考慮了解吸和擴散作用的滲流數(shù)學(xué)模型為[5]

    定義一組無因次參數(shù)

    模型無因次化后進行Laplace變換,可得

    1.3 考慮竄流、擴散綜合作用的數(shù)學(xué)模型

    此模型認為頁巖氣藏的基質(zhì)孔隙可以細分為有機和無機孔隙兩種類型,屬于三重介質(zhì)[6],擴散、竄流分開考慮,是一般化的模型,適用范圍較廣.滲流過程可分為

    1)氣井開始生產(chǎn)后,地層壓力下降,氣體在有機質(zhì)表面向內(nèi)解吸,進入有機質(zhì)孔隙;

    2)解吸的氣體使得有機質(zhì)孔隙內(nèi)氣體濃度增大,在濃度差的作用下氣體擴散進入裂縫系統(tǒng);

    3)無機基質(zhì)孔隙內(nèi)存在的游離氣在壓力差作用下竄流進入裂縫;

    4)裂縫中的氣體在流體勢的作用下遵循達西規(guī)律流向井筒.引入擬壓力后,該滲流機理下對應(yīng)的滲流方程為

    進行無因次化時,需要將上一個模型中的部分無因次參數(shù)重新定義

    模型化簡為

    進行Laplace變換得

    2 點源解的求取

    通過引入各自不同的參數(shù)對上述三種模型化簡,最終模型對應(yīng)的基本形式一致.根據(jù)點源函數(shù)的基本理論[7-9]可直接寫出三種情況下點源產(chǎn)生的壓力響應(yīng)公式.

    針對不同的滲流機理只需將f(s)進行替換即可.

    3 垂直裂縫井模型

    3.1 模型假設(shè)

    模型的假設(shè)條件主要包括:①油藏上下為不滲透邊界,外邊界考慮成無限大情況;②裂縫兩翼對稱,沿井軸均勻分布,并且完全壓開氣層,裂縫半長為yf,縫高等于氣層厚度h,裂縫為無限導(dǎo)流能力裂縫,中心在(xw,yw,zw)處.其余假設(shè)均符合點源解推導(dǎo)時的規(guī)定.

    3.2 垂直裂縫井壓力響應(yīng)

    若裂縫為均勻流量模型,可將點源壓力響應(yīng)沿縫高和縫長方向積分,得到垂直裂縫井壓力響應(yīng),并將流量進行無因次化,最終得到(lref=y(tǒng)f):

    式(13)是Laplace空間中垂直裂縫井在地層中任意一點產(chǎn)生的壓力響應(yīng),通過數(shù)值反演[10]可得到井底壓力響應(yīng).根據(jù)無限導(dǎo)流能力裂縫模型與均勻流量模型的關(guān)系,上式中yD=0.732處的壓力即可看作是無限導(dǎo)流能力裂縫模型的井底壓力.

    4 壓裂水平井模型

    4.1 模型假設(shè)

    模型的示意圖如圖1,假設(shè)條件主要包括:①油藏上下為不滲透邊界,外邊界無限大;②一口水平井(長度為L)被不可變形的橫向裂縫(條數(shù)n,且為單數(shù))垂直貫穿,裂縫兩翼對稱,沿井軸均勻分布,并且完全壓開氣層,裂縫半長為yf,縫高等于氣層厚度h,裂縫為無限導(dǎo)流能力裂縫.其余假設(shè)均符合點源解推導(dǎo)時的規(guī)定.

    圖1 壓裂水平井及坐標示意圖Fig.1 Fractured horizontal well and coordinate system

    4.2 壓裂水平井壓力響應(yīng)

    以水平無限大油藏為例,對于均勻流量模型,第j條裂縫的流體流動對第i條裂縫上的某一點(xDi,yDi)產(chǎn)生的壓力響應(yīng)為在每條裂縫上另yD=0.732,可近似得到無限導(dǎo)流能力裂縫井底不穩(wěn)定壓力響應(yīng)[11].

    上式可寫出n+1個方程,需要求解n+1個未知數(shù),即n個qDi和1個.對于其他類邊界的油藏,可用類似的方法得到井底壓力解[12].如需考慮井筒存儲效應(yīng)和表皮系數(shù),則可根據(jù)如下公式[13]

    5 壓力動態(tài)分析

    5.1 垂直壓裂井壓力動態(tài)

    根據(jù)第3部分編程可繪制出頁巖氣藏中垂直壓裂井的壓力動態(tài)響應(yīng),如圖2.

    分析圖2發(fā)現(xiàn)井筒存儲效應(yīng)過后,出現(xiàn)裂縫線性流,繼而出現(xiàn)兩個凹陷段,最終達到邊界反應(yīng)階段.三種機理下壓力動態(tài)曲線的區(qū)別體現(xiàn)在竄流階段.前兩種機理與雙重介質(zhì)油藏壓力動態(tài)曲線類似,第三種滲流機理下曲線出現(xiàn)兩個“凹子”,表現(xiàn)出三重介質(zhì)特點[14].

    分析三種模型的滲流機理,初期階段采出的均是裂縫中的游離氣,竄流和擴散現(xiàn)象還未顯現(xiàn),三種模型的壓力動態(tài)曲線重合;隨著裂縫中游離氣的采出,裂縫和基質(zhì)之間出現(xiàn)壓力差和濃度差.機理1認為解吸的氣體與原基質(zhì)內(nèi)氣體混合,壓力差作用下竄流進入天然裂縫,補充裂縫中的氣體,使得壓力導(dǎo)數(shù)曲線最先出現(xiàn)“凹子”,機理3也考慮到了竄流作用,但是由于它將基質(zhì)細分成了無機基質(zhì)、有機基質(zhì)兩種類型,僅無機基質(zhì)中的氣體發(fā)生竄流作用,因此與機理1相比,“凹子”更淺、更窄,而且出現(xiàn)的時間略晚;接下來擴散作用出現(xiàn),機理2對應(yīng)的曲線出現(xiàn)“凹子”,與機理2的曲線相比,機理3曲線的第二個“凹子”較小,出現(xiàn)時間也略晚,原因與前一階段類似,也是由于基質(zhì)考慮方式不同引起的.

    5.2 壓裂水平井壓力動態(tài)

    壓裂水平井在頁巖氣藏開采中占有十分重要的地位,對其壓力動態(tài)做詳細分析.由于第三種模型將解吸、吸附、擴散全部考慮在內(nèi),以第三種滲流機理為基礎(chǔ),分析頁巖氣藏壓裂水平井的壓力動態(tài),如圖3.分析發(fā)現(xiàn)壓力動態(tài)曲線分為7個階段:①井筒存儲階段;②過渡階段;③裂縫之間地層線性流階段;④裂縫之間地層擬徑向流動階段(該階段只有在裂縫間距取得較大時出現(xiàn),為一縱坐標為(1/2n)的水平線,此處裂縫條數(shù)n=5);⑤擴散、串流階段(圖中兩種作用出現(xiàn)的時機相近,“凹子”重合.);⑥地層整體線性流動階段;⑦地層整體擬徑向流動階段.

    圖2 無限大邊界頁巖氣藏垂直壓裂井壓力動態(tài)曲線Fig.2 Typical curves for vertical fractured well in shale reservoir with infinite boundary

    圖3 無限大頁巖氣藏中壓裂水平井壓力動態(tài)曲線Fig.3 Typical curves for fractured horizontalwell in shale reservoir with infinite boundary

    由于裂縫間距等取值的不同可能會使得“凹子”出現(xiàn)的時間發(fā)生變化,即可能出現(xiàn)在裂縫徑向流之后、地層線性流之后或其它時期.如圖4,降低了裂縫間距,與圖3不同,“凹子”出現(xiàn)在地層線性流階段之后.分析壓裂水平井的壓力動態(tài),三種機理的壓力曲線在除井筒存儲以及邊界反映階段以外的其它階段均有所不同,其原因與垂直壓裂井模型類似,是由于基質(zhì)類型不同導(dǎo)致了滲流機理出現(xiàn)差異,使得壓力導(dǎo)數(shù)曲線出現(xiàn)多次交叉.但是不同機理作用下的壓力曲線的趨勢大體相似.

    對于垂直壓裂井和壓裂水平井,擴散系數(shù)、解吸系數(shù)、竄流系數(shù)的取值不同可能會改變兩個“凹子”的大小以及出現(xiàn)的先后順序.以壓裂水平井為例,通過改變解吸系數(shù)的值,得到圖5所示結(jié)果.

    圖4 不同滲流機理下壓裂水平井的壓力動態(tài)曲線Fig.4 Typical curves for fractured horizontalwellwithdifferent seepagemechanisms

    圖5 解吸系數(shù)對壓力導(dǎo)數(shù)曲線影響Fig.5 Effect of desorption coefficient on pressure derivative curves

    分析機理1和機理2所對應(yīng)曲線中“凹子”出現(xiàn)的時機,發(fā)現(xiàn)系數(shù)變化時,機理1的曲線不發(fā)生變化,α=1,機理1的模型先出現(xiàn)“凹子”,此時先發(fā)生竄流作用,且竄流作用較弱,“凹子”較小;α=20時,機理2對應(yīng)的模型先出現(xiàn)“凹子”,發(fā)生擴散作用,且擴散作用較強,“凹子”較大.對于機理3,發(fā)現(xiàn)當(dāng)解吸系數(shù)取某一值時(α=20),兩個“凹子”重合,且兩個作用疊合使“凹子”最大.

    6 結(jié)論

    1)頁巖氣藏存在著解吸、擴散、滲流等多種形式的氣體的流動,復(fù)雜的滲流機理造成了數(shù)學(xué)描述方法的差異,在分析不同的流動方式對滲流方程的影響的基礎(chǔ)上建立了三種滲流數(shù)學(xué)模型.

    2)三種數(shù)學(xué)模型的不同點主要體現(xiàn)在竄流項上,將其融合在f(s)表達式中,模型可化簡成同一種形式,利用常規(guī)油氣藏的求解方法可寫出頁巖氣藏的地層壓力解.

    3)根據(jù)源函數(shù)的基本理論寫出了無限大油藏的點源解,通過積分得到了垂直壓裂井和壓裂水平井的井底壓力表達式,編程繪制了壓力動態(tài)曲線,著重分析了壓裂水平井的壓力動態(tài),分析得出三種模型壓力動態(tài)曲線的主要區(qū)別出現(xiàn)在竄流階段.

    7 參數(shù)說明

    ψ:擬壓力,Pa;φ:裂縫系統(tǒng)孔隙度;k:裂縫系統(tǒng)和基質(zhì)系統(tǒng)孔隙度,m2;Cg:氣體等溫壓縮系數(shù),Pa-1;ρg:氣體密度,kg·m-3;VL:Langmuir體積,sm3·m-3;pL:Langmuir壓力,Pa;pgi:初始壓力,Pa;:修正的綜合壓縮系數(shù),Pa-1;t:時間,s;r:徑向距離,m;a:權(quán)重系數(shù),可取0.1;λ′:竄流因子;α0:形狀因子;Lref:參考長度,m;ω:彈性儲容比;pg:氣藏壓力,Pa;μ:氣體粘度,Pa·s;Z:氣體壓縮因子;T:溫度;V:濃度,sm3/m3;D:擴散系數(shù),m2·s-1;Bg:氣體體積系數(shù);q:流量,m3·s-1;α:解吸系數(shù);s:Laplace變量;:解吸時間,s;n:裂縫條數(shù);下標說明:f:裂縫系統(tǒng);m:基質(zhì)系統(tǒng);sc:標準狀況;w:井底;D:無因次參數(shù);i:初始時刻.

    [1] Duan Yonggang,WeiMingqiang,Li Jianqiu,et al.Shale gas seepagemechanism and fractured wells'production evaluation [J].Journal of Chongqing University,2011,34(4):62-66.

    [2] Chao Gao,John LeeW.Modelingmultilayer gas reservoirs including sorption effects[C].SPE Eastern Regional Meeting,8 -10 November,1994,Charleston,West Vinginin:Paper SPE29173,1994.

    [3] Bumb A C,Mckee C R.Gas well testing in the presence of desorption for coalbed methane and Devonian shale[J]. SPE15227,1998,3(11):179-185.

    [4] Duan Yonggang,Li Jianqiu.Transient pressure analysis of infinite conductivity fractured wells for shale gas[J]Natural Gas Industry,2010,30(3):26-29.

    [5] Ren Junjie,Guo Ping,Wang Delong,etal.Productivitymodelof fractured horizontalwells in shale gas reservoirs and analysis of influential factors[J].Journal of Northeast Petroleum University,2012,36(6):76-81.

    [6] Cheng Yuanfang,Dong Bingxiang,Shi Xian,et al.Seepagemechanism of a triple?porosiy/dual permeability model for shale gas reservoirs[J].Natural Gas Industry,2012,32(9):44-47.

    [7] Gringarten A C,Ramey H J.The use of source and green function in solving unsteady state problems in the reservoirs[J]. SPE3818,1973,13(5):285-296.

    [8] Ozkane,Raghavanr.Some new solutions to solve problems in well test analysis:Part 1?Analytical considerations[J].SPE Formation Evaluation,1991,6(3):359-368.

    [9] Ozkane,Raghavanr.Some new solutions to solve problems in well testanalysis:Part2?Computational considerations[J].SPE Formation Evaluation,1991,6(3):369-378.

    [10] Tong Dengke,Chen Qinlei.Some annotation about Stehfest numerical inversion[J].Acta Petrolei Sinica,2001,22(6):91 -92.

    [11] Lian Peiqing,Cheng Linsong,Cao Renyi,et al.A coupling model of low permeability reservoir and fractured horizontal wellbore in nonsteady state[J].Chinese JComput Phys,2007,27(2):203-210.

    [12] Chen Linsong,Pi Jian,Lian Peiqing,et al.A computationalmethod for productivity of horizontal well in naturally fractured reservoirs[J].Chinese JComput Phys,2011,28(2):230-236.

    [13] Brown M,Ozkan E,et al.Practical solutions for pressure transient responses of fractured horizontal wells in un?conventional reservoirs[J].SPE125043,2011,14(6):663-676.

    [14] Zhao Dongmei,Yao Jun,Wang Zisheng.The pressure?transient characteristics of a stress?sensitive triple?medium reservoir [J].Chinese JComput Phys,2005,22(5):444-448.

    Shale Gas Seepage M echanism and Transient Pressure Analysis

    DU Dianfa,WANG Yanyan,F(xiàn)U Jingang,SUN Zhaobo,QIAO Ni (College ofPetroleum Engineering in China University ofPetroleum,Shandong Qingdao 266580,China)

    According to characteristic of shale gas reservoir,mathematical description methods of diffusion and desorption are analyzed.Three seepage models controlled by different kinds of seepage mechanisms are built:① a mathematical model which considers diffusion and interporosity of desorbed gas and free gas,②amathematicalmodelwhich considers desorption and diffusion,③amathematical model which considers desorption and diffusion in organic porosity and interporosity in inorganic porosity.With modified Lord Kelvin point source function pressure solutions are botained.With corresponding transformation for point source functions,formation pressure solutions for vertical fractured well and fractured horizontal well are obtained and typical curves are shown.At last,characteristics of different curves and threemechanisms'influences on pressure are analyzed.

    shale reservoir;seepagemechanism;mathematicalmodel;well test analysis

    TE33

    A

    2013-12-30;

    2014-06-13

    中國石油科技創(chuàng)新基金(2011D?5006?0210)資助項目

    杜殿發(fā)(1972-),男,博士,教授,主要從事油氣田開發(fā)等研究及教學(xué)工作,E?mail:dudf@upc.edu.cn

    Received date: 2013-12-30;Revised date: 2014-06-13

    猜你喜歡
    氣藏機理滲流
    隔熱纖維材料的隔熱機理及其應(yīng)用
    煤層氣吸附-解吸機理再認識
    中國煤層氣(2019年2期)2019-08-27 00:59:30
    霧霾機理之問
    致密氣藏壓裂傾斜縫壓力動態(tài)分析
    DNTF-CMDB推進劑的燃燒機理
    塔里木油田超深超高壓氣藏的成功改造
    斷塊油氣田(2014年6期)2014-03-11 15:34:03
    LG地區(qū)復(fù)雜礁灘氣藏描述技術(shù)
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    尾礦壩滲流計算及排滲設(shè)計
    金屬礦山(2013年6期)2013-03-11 16:54:05
    人妻一区二区av| 亚洲av一区综合| 亚洲色图av天堂| 欧美一级a爱片免费观看看| 狂野欧美白嫩少妇大欣赏| 激情五月婷婷亚洲| av播播在线观看一区| 国产毛片a区久久久久| 青春草视频在线免费观看| 日韩欧美精品免费久久| 日本欧美国产在线视频| 成人鲁丝片一二三区免费| 久久人人爽人人爽人人片va| 夜夜看夜夜爽夜夜摸| 国产伦精品一区二区三区四那| 简卡轻食公司| 有码 亚洲区| 中文欧美无线码| 久久久亚洲精品成人影院| 亚洲在久久综合| 久久99热6这里只有精品| 99九九线精品视频在线观看视频| 精品人妻视频免费看| 国产精品一区二区在线观看99 | 如何舔出高潮| 日日摸夜夜添夜夜爱| 日本wwww免费看| 成人毛片a级毛片在线播放| 内射极品少妇av片p| 丝瓜视频免费看黄片| 婷婷色综合www| 日韩不卡一区二区三区视频在线| 精品少妇黑人巨大在线播放| 免费观看精品视频网站| 成人二区视频| 亚洲精品影视一区二区三区av| 国产伦精品一区二区三区视频9| 亚洲国产高清在线一区二区三| 成人综合一区亚洲| 亚洲精品国产av蜜桃| 成人性生交大片免费视频hd| 激情五月婷婷亚洲| 亚洲精品,欧美精品| 精品国产一区二区三区久久久樱花 | 亚洲最大成人av| 国产 一区精品| 欧美精品一区二区大全| 亚洲国产日韩欧美精品在线观看| 亚洲国产成人一精品久久久| 亚洲欧美中文字幕日韩二区| 搡老乐熟女国产| 色综合亚洲欧美另类图片| 亚洲激情五月婷婷啪啪| 久久久久久久久久久免费av| 欧美潮喷喷水| 久久草成人影院| 午夜福利视频精品| 搡老妇女老女人老熟妇| 中文精品一卡2卡3卡4更新| av播播在线观看一区| 亚洲欧美精品自产自拍| av国产免费在线观看| 91狼人影院| 成人午夜精彩视频在线观看| 成年av动漫网址| 成人二区视频| 51国产日韩欧美| 国产美女午夜福利| 欧美zozozo另类| 精品国产三级普通话版| 免费看a级黄色片| 男人舔奶头视频| 波野结衣二区三区在线| 成年女人在线观看亚洲视频 | 一本一本综合久久| 女人十人毛片免费观看3o分钟| 国产一区二区亚洲精品在线观看| 在线观看av片永久免费下载| 人妻制服诱惑在线中文字幕| 亚洲图色成人| 最近的中文字幕免费完整| 青青草视频在线视频观看| 亚洲怡红院男人天堂| 人妻夜夜爽99麻豆av| 大又大粗又爽又黄少妇毛片口| 午夜福利视频1000在线观看| 久久久精品94久久精品| 免费电影在线观看免费观看| 亚洲精品亚洲一区二区| 十八禁网站网址无遮挡 | 麻豆成人av视频| 久久久久久久国产电影| 国产淫片久久久久久久久| 国产精品三级大全| 精华霜和精华液先用哪个| 亚洲成人精品中文字幕电影| 一个人看的www免费观看视频| 人人妻人人看人人澡| 久久99热6这里只有精品| 99久久中文字幕三级久久日本| 久久精品国产亚洲av涩爱| 一个人观看的视频www高清免费观看| 夫妻午夜视频| 亚洲av男天堂| 自拍偷自拍亚洲精品老妇| 国产av码专区亚洲av| 成人一区二区视频在线观看| 成人亚洲精品av一区二区| 午夜爱爱视频在线播放| 一区二区三区乱码不卡18| 日本爱情动作片www.在线观看| 亚洲av成人av| 2018国产大陆天天弄谢| 少妇丰满av| 搡女人真爽免费视频火全软件| 精品久久久久久久久亚洲| 亚洲精品日韩av片在线观看| 精品人妻视频免费看| 97人妻精品一区二区三区麻豆| 五月伊人婷婷丁香| 国产精品av视频在线免费观看| 日本wwww免费看| 日韩欧美三级三区| 国产一区亚洲一区在线观看| 国产精品一及| 又粗又硬又长又爽又黄的视频| 国产亚洲91精品色在线| 国产一区亚洲一区在线观看| 亚州av有码| 激情五月婷婷亚洲| 有码 亚洲区| 久久久精品94久久精品| 青春草亚洲视频在线观看| 国产一区亚洲一区在线观看| 九九久久精品国产亚洲av麻豆| 91aial.com中文字幕在线观看| 高清在线视频一区二区三区| 91精品国产九色| 国产淫片久久久久久久久| 精品人妻熟女av久视频| 久久久色成人| av免费观看日本| 只有这里有精品99| 91久久精品电影网| 精品国产一区二区三区久久久樱花 | 国产视频首页在线观看| 久久久久久久久中文| 好男人视频免费观看在线| 综合色丁香网| 亚洲怡红院男人天堂| 美女大奶头视频| 午夜福利网站1000一区二区三区| 午夜福利网站1000一区二区三区| av.在线天堂| 少妇熟女aⅴ在线视频| 黄色欧美视频在线观看| 亚洲av在线观看美女高潮| 亚洲欧美一区二区三区国产| 欧美性猛交╳xxx乱大交人| 99久久中文字幕三级久久日本| 亚洲最大成人中文| 免费看不卡的av| 国产精品蜜桃在线观看| 久久久久久久国产电影| 国产单亲对白刺激| 国产亚洲一区二区精品| 男人狂女人下面高潮的视频| 亚洲欧美一区二区三区国产| 天天躁日日操中文字幕| 91精品伊人久久大香线蕉| 观看美女的网站| 丰满乱子伦码专区| 国产av不卡久久| 国产成人91sexporn| 国产极品天堂在线| 国产一级毛片七仙女欲春2| 熟妇人妻久久中文字幕3abv| 中文字幕av在线有码专区| 久久这里有精品视频免费| 日本黄色片子视频| 最近手机中文字幕大全| 日本爱情动作片www.在线观看| 亚洲怡红院男人天堂| 久久久久久九九精品二区国产| 午夜福利网站1000一区二区三区| 日本免费在线观看一区| 不卡视频在线观看欧美| 亚洲在线自拍视频| 午夜精品国产一区二区电影 | 97精品久久久久久久久久精品| 嫩草影院新地址| 亚洲自偷自拍三级| 亚洲成人中文字幕在线播放| 日韩av在线免费看完整版不卡| 亚洲av二区三区四区| 亚洲国产日韩欧美精品在线观看| 一级爰片在线观看| 婷婷色av中文字幕| 性色avwww在线观看| 色综合站精品国产| av免费在线看不卡| 婷婷色综合大香蕉| 国产熟女欧美一区二区| 熟妇人妻久久中文字幕3abv| 内射极品少妇av片p| 小蜜桃在线观看免费完整版高清| 日本一本二区三区精品| 99热全是精品| 国语对白做爰xxxⅹ性视频网站| 亚洲精品久久午夜乱码| 国产男人的电影天堂91| 91av网一区二区| 汤姆久久久久久久影院中文字幕 | 国产欧美日韩精品一区二区| 麻豆精品久久久久久蜜桃| 国产亚洲精品久久久com| 日日摸夜夜添夜夜爱| 激情 狠狠 欧美| 最近视频中文字幕2019在线8| 免费电影在线观看免费观看| 亚洲成人精品中文字幕电影| 人人妻人人澡人人爽人人夜夜 | 精品久久久久久久久av| 日韩精品青青久久久久久| 色5月婷婷丁香| 少妇的逼水好多| 国产免费又黄又爽又色| 欧美 日韩 精品 国产| 成人av在线播放网站| 亚洲婷婷狠狠爱综合网| 一区二区三区四区激情视频| 中国国产av一级| 国产成人免费观看mmmm| 亚洲美女视频黄频| 国产免费一级a男人的天堂| 老师上课跳d突然被开到最大视频| 女人被狂操c到高潮| 国产一级毛片在线| 色视频www国产| 啦啦啦啦在线视频资源| 99热6这里只有精品| 国内精品宾馆在线| 在线免费十八禁| 在线免费十八禁| 能在线免费观看的黄片| 国产成人a∨麻豆精品| 亚洲婷婷狠狠爱综合网| 免费观看无遮挡的男女| 国产成人福利小说| 一级毛片 在线播放| 亚洲av免费在线观看| 99久久中文字幕三级久久日本| 美女黄网站色视频| 在线观看一区二区三区| 日韩不卡一区二区三区视频在线| 男人舔奶头视频| 国产真实伦视频高清在线观看| 国产亚洲5aaaaa淫片| 国产极品天堂在线| eeuss影院久久| 日本wwww免费看| 久久久久精品久久久久真实原创| 精华霜和精华液先用哪个| 日韩精品有码人妻一区| eeuss影院久久| 黄片wwwwww| 草草在线视频免费看| 少妇人妻一区二区三区视频| 国产亚洲5aaaaa淫片| 日本色播在线视频| 欧美日韩精品成人综合77777| 亚洲国产精品专区欧美| 久久久久久久久中文| 成人国产麻豆网| 国产欧美日韩精品一区二区| 成人欧美大片| 午夜免费观看性视频| 大又大粗又爽又黄少妇毛片口| 国产淫语在线视频| 午夜免费男女啪啪视频观看| 午夜视频国产福利| 99re6热这里在线精品视频| 日韩欧美一区视频在线观看 | 欧美一级a爱片免费观看看| 精品一区二区三卡| 中文字幕人妻熟人妻熟丝袜美| av在线老鸭窝| 80岁老熟妇乱子伦牲交| 中国美白少妇内射xxxbb| 看免费成人av毛片| 天天一区二区日本电影三级| 国产精品1区2区在线观看.| 成年女人看的毛片在线观看| 寂寞人妻少妇视频99o| 日韩精品青青久久久久久| 久久精品人妻少妇| 国产黄频视频在线观看| 亚洲成人精品中文字幕电影| 黄片wwwwww| 国产 亚洲一区二区三区 | 国产一区二区三区av在线| 欧美日韩精品成人综合77777| 日韩制服骚丝袜av| 一边亲一边摸免费视频| 18禁在线播放成人免费| 久久97久久精品| 男插女下体视频免费在线播放| 精品人妻熟女av久视频| 水蜜桃什么品种好| 中文乱码字字幕精品一区二区三区 | 一级毛片久久久久久久久女| 69人妻影院| 午夜福利在线观看吧| 亚洲国产色片| 久久精品国产亚洲网站| 亚洲精品一区蜜桃| 熟女人妻精品中文字幕| 麻豆成人午夜福利视频| 亚洲四区av| 最近的中文字幕免费完整| 欧美 日韩 精品 国产| 真实男女啪啪啪动态图| 欧美97在线视频| 欧美潮喷喷水| 国产av在哪里看| av国产免费在线观看| 亚洲av男天堂| 男女视频在线观看网站免费| 2021少妇久久久久久久久久久| 国产v大片淫在线免费观看| 免费av观看视频| 欧美日韩综合久久久久久| 在线观看av片永久免费下载| 26uuu在线亚洲综合色| 日韩欧美三级三区| 国产男人的电影天堂91| 国产黄色小视频在线观看| 欧美性感艳星| 校园人妻丝袜中文字幕| 99久久精品国产国产毛片| 中文字幕制服av| 欧美日本视频| 国产免费福利视频在线观看| 免费观看a级毛片全部| 久久热精品热| 草草在线视频免费看| 丰满乱子伦码专区| 久久韩国三级中文字幕| 久久精品夜色国产| 丝瓜视频免费看黄片| 啦啦啦韩国在线观看视频| 欧美+日韩+精品| 午夜福利网站1000一区二区三区| 国产不卡一卡二| 日韩精品有码人妻一区| 国产精品日韩av在线免费观看| 一个人观看的视频www高清免费观看| 青青草视频在线视频观看| 国产精品嫩草影院av在线观看| 免费看日本二区| 中文字幕人妻熟人妻熟丝袜美| 免费大片黄手机在线观看| 日产精品乱码卡一卡2卡三| 免费不卡的大黄色大毛片视频在线观看 | 亚洲图色成人| 日韩视频在线欧美| 国产成人午夜福利电影在线观看| 成人亚洲欧美一区二区av| eeuss影院久久| 亚洲18禁久久av| 欧美zozozo另类| 亚洲欧美一区二区三区国产| 高清在线视频一区二区三区| 十八禁网站网址无遮挡 | 亚洲国产精品专区欧美| 久久久久久久久久黄片| 日本免费a在线| 国产老妇伦熟女老妇高清| 99久国产av精品| 波多野结衣巨乳人妻| 精品不卡国产一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 非洲黑人性xxxx精品又粗又长| 高清在线视频一区二区三区| 久久久久久久久久成人| 国产av码专区亚洲av| 日产精品乱码卡一卡2卡三| 高清av免费在线| 91久久精品国产一区二区三区| 欧美成人一区二区免费高清观看| 特大巨黑吊av在线直播| 80岁老熟妇乱子伦牲交| 国产白丝娇喘喷水9色精品| 国产精品一区二区在线观看99 | av线在线观看网站| 国产成人精品一,二区| 免费少妇av软件| 伊人久久精品亚洲午夜| 国产精品久久久久久久电影| 91aial.com中文字幕在线观看| 精品午夜福利在线看| 99九九线精品视频在线观看视频| 成人性生交大片免费视频hd| 久久人人爽人人片av| 天堂影院成人在线观看| 国产午夜精品一二区理论片| 国产男女超爽视频在线观看| 亚洲真实伦在线观看| 精品一区二区免费观看| 视频中文字幕在线观看| 国产精品美女特级片免费视频播放器| 国语对白做爰xxxⅹ性视频网站| 色视频www国产| 国产色婷婷99| 亚洲人成网站在线观看播放| 日韩一区二区视频免费看| 国产精品不卡视频一区二区| 国产精品1区2区在线观看.| 日产精品乱码卡一卡2卡三| 丰满人妻一区二区三区视频av| 97在线视频观看| 午夜福利在线观看吧| 一级av片app| 国产精品不卡视频一区二区| 亚洲欧美一区二区三区黑人 | 波多野结衣巨乳人妻| 亚洲精品亚洲一区二区| 美女内射精品一级片tv| 亚洲av成人精品一区久久| av播播在线观看一区| 国产av国产精品国产| 一级黄片播放器| 晚上一个人看的免费电影| 成人一区二区视频在线观看| 夫妻性生交免费视频一级片| 亚洲国产精品sss在线观看| av一本久久久久| 嫩草影院精品99| 国产精品美女特级片免费视频播放器| 联通29元200g的流量卡| 国产白丝娇喘喷水9色精品| 亚洲欧美成人综合另类久久久| 免费看美女性在线毛片视频| 国内精品宾馆在线| 国产 一区 欧美 日韩| 日韩伦理黄色片| 亚洲精品久久久久久婷婷小说| 日韩三级伦理在线观看| 能在线免费看毛片的网站| 亚洲四区av| 国产一区有黄有色的免费视频 | 亚洲精品乱码久久久v下载方式| 中文在线观看免费www的网站| 黄色一级大片看看| 亚洲精品日韩av片在线观看| 99久久精品热视频| 99九九线精品视频在线观看视频| 在线免费十八禁| 中文在线观看免费www的网站| 肉色欧美久久久久久久蜜桃 | 99九九线精品视频在线观看视频| 五月天丁香电影| 别揉我奶头 嗯啊视频| 亚洲欧美日韩东京热| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件| 午夜福利成人在线免费观看| 久久久久国产网址| 日韩制服骚丝袜av| 美女cb高潮喷水在线观看| 伊人久久国产一区二区| 欧美日韩一区二区视频在线观看视频在线 | 色吧在线观看| 男女视频在线观看网站免费| 蜜桃久久精品国产亚洲av| 女人久久www免费人成看片| 3wmmmm亚洲av在线观看| 欧美激情在线99| 天堂影院成人在线观看| 亚洲成人av在线免费| 色综合亚洲欧美另类图片| 一区二区三区乱码不卡18| 国产成人精品一,二区| 亚洲av日韩在线播放| 99久久精品热视频| 精品久久久精品久久久| 免费少妇av软件| 亚洲欧美一区二区三区黑人 | 秋霞伦理黄片| 国产成人免费观看mmmm| 18禁动态无遮挡网站| 亚洲精品视频女| 国产老妇伦熟女老妇高清| av福利片在线观看| 深爱激情五月婷婷| 国产伦精品一区二区三区四那| 精品一区二区免费观看| 国产精品久久视频播放| 99热这里只有精品一区| 精品少妇黑人巨大在线播放| 哪个播放器可以免费观看大片| 精品久久久久久久末码| 国产伦一二天堂av在线观看| 日韩强制内射视频| 天堂av国产一区二区熟女人妻| 国产成人福利小说| 亚洲在线观看片| 国产真实伦视频高清在线观看| 99热全是精品| 亚洲欧美成人精品一区二区| 国产亚洲5aaaaa淫片| 一个人观看的视频www高清免费观看| a级毛色黄片| 欧美xxⅹ黑人| 国产精品三级大全| 久久久久久久大尺度免费视频| av.在线天堂| 国产亚洲91精品色在线| 99热这里只有是精品50| 亚洲欧美精品专区久久| 亚洲四区av| 久久鲁丝午夜福利片| 国产av不卡久久| 亚洲国产成人一精品久久久| ponron亚洲| 青春草亚洲视频在线观看| 久久久a久久爽久久v久久| 麻豆成人av视频| 久久人人爽人人爽人人片va| 联通29元200g的流量卡| 一级毛片电影观看| 精品久久久久久久人妻蜜臀av| 日韩欧美精品v在线| 天堂俺去俺来也www色官网 | 视频中文字幕在线观看| 特级一级黄色大片| 女人十人毛片免费观看3o分钟| 精品少妇黑人巨大在线播放| 午夜亚洲福利在线播放| 99re6热这里在线精品视频| 美女xxoo啪啪120秒动态图| 在线天堂最新版资源| 18禁裸乳无遮挡免费网站照片| 夜夜看夜夜爽夜夜摸| 老司机影院成人| 精品人妻偷拍中文字幕| 国产精品国产三级国产专区5o| freevideosex欧美| 2021少妇久久久久久久久久久| 中文天堂在线官网| 青春草国产在线视频| 亚洲熟女精品中文字幕| 99热这里只有精品一区| 国产成人91sexporn| 99久久人妻综合| 免费av毛片视频| 女的被弄到高潮叫床怎么办| 哪个播放器可以免费观看大片| 国产91av在线免费观看| 99热网站在线观看| 日韩欧美国产在线观看| 观看美女的网站| 99久久九九国产精品国产免费| 欧美成人一区二区免费高清观看| av.在线天堂| 成人性生交大片免费视频hd| 午夜福利在线在线| 干丝袜人妻中文字幕| 国产欧美另类精品又又久久亚洲欧美| 国产淫语在线视频| 国产精品爽爽va在线观看网站| 欧美+日韩+精品| 亚洲国产精品成人综合色| 最近中文字幕2019免费版| 午夜福利视频1000在线观看| 国产成人精品久久久久久| 久99久视频精品免费| 午夜福利高清视频| 99热6这里只有精品| 亚洲一级一片aⅴ在线观看| 夜夜爽夜夜爽视频| 久久久久久国产a免费观看| 91狼人影院| 久久精品久久精品一区二区三区| 国产精品1区2区在线观看.| 一级毛片我不卡| 国产精品久久久久久久电影| 国产精品嫩草影院av在线观看| 少妇丰满av| 国产乱人偷精品视频| 亚州av有码| 久久精品熟女亚洲av麻豆精品 | 啦啦啦中文免费视频观看日本| 亚洲最大成人中文| 国模一区二区三区四区视频| 欧美日韩亚洲高清精品| 国产免费一级a男人的天堂| 简卡轻食公司| 亚洲乱码一区二区免费版| 麻豆av噜噜一区二区三区| 神马国产精品三级电影在线观看| 国产精品嫩草影院av在线观看| 久久精品人妻少妇| 啦啦啦啦在线视频资源| 国产一区二区在线观看日韩| 啦啦啦啦在线视频资源| 国产三级在线视频| 国产视频内射| 少妇人妻一区二区三区视频| 国产片特级美女逼逼视频| 亚洲精品自拍成人| 天堂影院成人在线观看| 一级a做视频免费观看| 国产精品1区2区在线观看.| 麻豆成人午夜福利视频|