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

    基于數(shù)字巖心技術(shù)研究低滲砂巖滲流特征

    2017-08-30 18:16:02王春生孫啟冀
    物探化探計算技術(shù) 2017年4期
    關(guān)鍵詞:巖心滲流巖石

    王春生, 劉 洋, 孫啟冀, 張 凱, 梁 超

    (東北石油大學(xué) 石油工程學(xué)院, 大慶 163318)

    基于數(shù)字巖心技術(shù)研究低滲砂巖滲流特征

    王春生, 劉 洋, 孫啟冀, 張 凱, 梁 超

    (東北石油大學(xué) 石油工程學(xué)院, 大慶 163318)

    以傳統(tǒng)巖石物理實(shí)驗(yàn)研究低孔、低滲、高復(fù)雜性孔隙結(jié)構(gòu)儲層的宏觀巖石物理性質(zhì),不僅成本高、周期長、誤差大且難以定量研究儲層微觀因素對巖石物理屬性的影像規(guī)律。通過Micro-CT掃描與高級成像處理技術(shù),構(gòu)建與真實(shí)巖心相等價的精確代表性孔隙空間模型,利用Avizo-XLab求解器數(shù)值模擬計算以N-S為控制方程的孔隙流體流動及以達(dá)西定律計算巖石滲流參數(shù)。通過計算結(jié)果的可視化處理,以三維流線展示孔隙空間流體流動,以高度梯度圖展示孔隙壓力分布場,使模擬結(jié)果清晰直觀。研究表明,基于數(shù)字平臺仿真孔隙微流動及測量巖石滲流參數(shù)具有高可信度。該研究方式彌補(bǔ)了傳統(tǒng)巖石物理實(shí)驗(yàn)的不足,為數(shù)字化研究巖石物理屬性提供借鑒。

    巖石物理實(shí)驗(yàn); 微觀因素; Micro-CT掃描; 孔隙微流動; 巖石滲流參數(shù)

    0 前言

    我國的低滲油藏分布廣泛、資源豐富、開發(fā)潛力巨大[1-2]。低滲儲層由于低孔、低滲、高孔隙結(jié)構(gòu)復(fù)雜性等特點(diǎn),給巖石物理實(shí)驗(yàn)帶來巨大挑戰(zhàn)[3]。為了提高低滲油藏原油采收率,急需從微觀層面上對巖石物理屬性進(jìn)行研究[4]。傳統(tǒng)巖石物理實(shí)驗(yàn)(如:滲流實(shí)驗(yàn)[5],壓汞實(shí)驗(yàn)[6]及核磁共振實(shí)驗(yàn)[7]等),無法定量描述巖石微觀性質(zhì),據(jù)此建立的很多理論模型均是以巖石骨架和孔隙為基礎(chǔ)的粗略模型,不能精準(zhǔn)描述孔隙空間分布,且傳統(tǒng)巖石物理實(shí)驗(yàn)具有測試周期長,數(shù)據(jù)結(jié)果誤差大、巖心重復(fù)實(shí)驗(yàn)次數(shù)低等弊端[8]。隨著高分辨率Micro-CT掃描技術(shù)及數(shù)字圖像處理技術(shù)的高速發(fā)展,基于數(shù)字巖心及模擬算法的數(shù)字巖石物理技術(shù),不僅能三維可視化巖石微觀結(jié)構(gòu)且能定量表征孔隙結(jié)構(gòu)參數(shù),更能快速、精準(zhǔn)獲取孔隙度、滲透率、巖電參數(shù)及毛管壓力等巖石物理參數(shù)[9]。同時,基于構(gòu)建的數(shù)字巖心,不僅能夠重復(fù)利用,更能實(shí)現(xiàn)同一巖心孔隙空間不同類型物理場的模擬與耦合,而且能使影響巖石物性的多種微觀因素變得可控,并易于定量考察,從而搭建起了聯(lián)系微觀與宏觀的橋梁,有助于更好地了解儲層特征、提高儲層評價準(zhǔn)確性及開采效率[10]。

    目前,通過高級成像技術(shù)構(gòu)建代表性巖心孔隙空間模型并以此數(shù)字平臺計算研究巖石物理屬性的方法,被國內(nèi)、外眾多學(xué)者所推崇[11-12]。為了估算滲透率,需先對孔隙流體流動進(jìn)行仿真模擬。此類研究主要有2種:①立足于宏觀角度,視孔隙內(nèi)流體為連續(xù)流體,用一組非線性偏微分方程,即反映流體質(zhì)、能、動量守恒的Navier-Stokes方程描述其宏觀運(yùn)動[13-14];②立足于微觀角度,視孔隙內(nèi)流體為大量無規(guī)則運(yùn)動的粒子,視流體流動為粒子間碰撞,經(jīng)過統(tǒng)計平均得到孔隙內(nèi)流體的宏觀運(yùn)動[15-16]。筆者通過巖樣Micro-CT掃描、孔隙空間三維重構(gòu)、數(shù)字滲流模擬,得到孔隙空間流體流動與巖石宏觀滲流特性等結(jié)果,為開展數(shù)字滲流模擬研究提供借鑒。

    1 數(shù)字巖心模型構(gòu)建

    采用Micro-CT構(gòu)建巖心孔隙空間結(jié)構(gòu)不僅準(zhǔn)確快速,而且無損巖樣。本次研究對象為國內(nèi)某油田的一塊天然低滲砂巖,首先對其進(jìn)行全尺寸(25 mm)CT掃描,對獲取的掃描切片,以截面孔隙度為依據(jù),整體考察巖心均質(zhì)性。由于該塊巖心面孔隙度沿軸向波動極小,故其均質(zhì)程度高。為了進(jìn)一步研究孔隙尺度滲流特征,在該塊巖心上隨機(jī)鉆取直徑為2 mm的巖樣,其基本物性參數(shù)見表1。通過德國GE公司進(jìn)口的柜式Micro-CT儀器以1 μm分辨率對鉆取的2 mm巖樣進(jìn)行掃描,用一組包含1 008張的灰度圖像以216個灰度等級可視化展示巖石內(nèi)部不同密度差異的微觀結(jié)構(gòu)。圖1為巖樣掃描切片的灰度直方圖與巖石不同密度材質(zhì)的對應(yīng)關(guān)系。

    表1 常規(guī)實(shí)驗(yàn)測取巖樣參數(shù)Tab.1 Using conventional experiment sample parameters

    圖1 不同密度材質(zhì)與灰度值的對應(yīng)關(guān)系Fig.1 The corresponding relation of different density material and grey value

    根據(jù)Micro-CT掃描獲取的巖心橫截面切片,運(yùn)用Avizo軟件主要進(jìn)行以下四步處理:

    1)灰度圖像濾波。采用中值濾波和高斯濾波去除切片圖像中系統(tǒng)噪點(diǎn)干擾,使巖石骨架和孔隙空間的過度明顯,更易清楚分辨。濾波后的灰度圖像中部分孤立孔隙及巖石骨架像素也相應(yīng)被剔除。

    2)灰度圖像二值化。根據(jù)巖樣灰度統(tǒng)計直方圖,選取位于兩峰之間的波谷對應(yīng)的灰度值,作為分割孔隙空間與巖石骨架的閾值;閾值分割后的二值圖像中數(shù)字1標(biāo)記孔隙空間,數(shù)字0標(biāo)記巖石骨架。通過二值圖像便可構(gòu)建數(shù)字巖心,最后通過驗(yàn)證數(shù)字巖心與真實(shí)巖石有效孔隙度的吻合度來確定閾值分割的精準(zhǔn)性。

    圖2 孔隙空間及球棒模型Fig.2 The model of pore space and sphere-stick(a)孔隙空間;(b)球棒模型

    圖3 巖心孔隙半徑與喉道長度頻率分布Fig.3 The frequency distribution of core pore radius and the channel length(a)巖石孔隙半徑;(b)喉道長度頻率分布

    3)表征單元構(gòu)建。由于計算機(jī)處理性能限制,要求選取合理尺寸的表征單元,既能使此尺寸大小的數(shù)字巖心代表巖石宏觀物理性質(zhì),又能適應(yīng)當(dāng)前計算平臺硬件水平。根據(jù)前人研究結(jié)果,本研究中選取數(shù)字巖心的表征單元體積為500×500×600體素。渲染后的三維數(shù)字巖心如圖2所示,該模型直接用于滲流模擬與絕對滲透率計算。

    4)球棒模型構(gòu)建。為了后續(xù)工作能夠直觀研究孔隙空間拓?fù)浣Y(jié)構(gòu)與流體空間分布間的聯(lián)系,在直接應(yīng)用三維數(shù)字巖心仿真微觀孔隙流動前,通過簡化孔隙為球體,簡化喉道為棍棒,以最大球算法提取出與表征單元體相等價的孔隙網(wǎng)絡(luò)模型(圖2)并對該研究模型的858個孔隙及1 279個喉道統(tǒng)計分析,定量描述了滲流發(fā)生的物理空間的詳細(xì)特征,即孔隙半徑及喉道長度分布特征(圖3)。

    2 巖石物性參數(shù)計算

    2.1 低滲透砂巖儲層流動控制方程

    以傳統(tǒng)計算流體力學(xué)(CFD)數(shù)值模擬計算低滲儲層孔隙模型中微觀流動的方法,目前發(fā)展成熟應(yīng)用廣泛?;谶B續(xù)假設(shè)的傳統(tǒng)CFD法的描述方程為Navier-Stokes方程:

    (1)

    (2)

    (3)

    式中:ρ表示流體密度;t表示時間;u表示流速;μ表示流體動力粘度;▽表示拉普拉斯算子;e表示內(nèi)能;σ表示應(yīng)力張量;q表示熱通量。本次模擬流體基于如下條件假設(shè):孔隙流動為流速較小且不隨時間變化的層流;流體為不可壓縮的牛頓流體。

    孔隙空間內(nèi)流體流動數(shù)值后,通過在進(jìn)口或出口邊界對流速積分獲取經(jīng)過孔隙空間的流體流量,結(jié)合達(dá)西公式,從而求得宏觀滲流參數(shù)。達(dá)西尺度滲流描述方程為式(4)。

    (4)

    式中:Q表示通過數(shù)巖總流量,m3·s-1;S表示滲流橫截面積,m2;k表示絕對滲透率,m2;μ表示流體動力粘度,Pa·s;ΔP表示壓差,Pa;L表示巖樣長度,m。

    2.2 應(yīng)用Avizo-XLab求解滲流特征

    本研究應(yīng)用Avizo-XLab求解器求解上述控制方程,該求解器直接基于構(gòu)成數(shù)字巖心三維圖像的最小單元——體積元素(體素)進(jìn)行求解,不需三維圖像逆向成實(shí)體,因此更不用劃分網(wǎng)格等其他求解方法中的前處理過程?,F(xiàn)以二維圖像闡述其求解原理,如圖4所示。本次求解過程在500×500×600個邊長為1 μm的立方體體素上進(jìn)行,邊界條件設(shè)置如表2所示。

    表2 邊界條件參數(shù)Tab.2 Boundary condition parameters

    2.3 模擬結(jié)果

    通過對模擬結(jié)果的后處理,得到沿Z軸正向及負(fù)向壓力場分布和流線圖。結(jié)果顯示正反滲流方向求得絕對滲透率相同,均為0.316 667×10-3μm2,與常規(guī)水測實(shí)驗(yàn)結(jié)果一致。圖5為Z+方向滲流仿真模擬圖。

    圖4 Avizo-XLab求解器求解原理圖Fig.4 The principle diagram of the Avizo-XLab solver

    通過滲流模擬結(jié)果圖6(a)可知,巖心內(nèi)部存在主要滲流通道,且主流通道趨勢與大孔徑孔隙分布吻合。結(jié)合孔隙網(wǎng)絡(luò)模型與流線分布圖6(a)可知,滲流流速最大處位于與較大孔隙相鄰的喉道中。

    圖5 Z+方向滲流數(shù)值模擬Fig.5 The numerical simulation of seepage follow Z+ direction

    圖6 孔隙流動流線分布與壓力分布圖Fig.6 Pore flow streamline distribution and pressure profile(a)、(b)孔隙流動流線分布圖;(c)壓力分布圖

    3 結(jié)論

    本次研究基于Micro-CT掃描,利用閾值分割算法三維重構(gòu)復(fù)雜孔隙空間幾何模型,應(yīng)用Avizo-XLab求解器計算以N-S為控制方程的孔隙流體流動及以達(dá)西定律計算巖石滲流參數(shù)。通過數(shù)字模擬計算出的絕對滲透率與傳統(tǒng)滲流實(shí)驗(yàn)結(jié)果數(shù)據(jù)比對結(jié)果的一致性,驗(yàn)證了以數(shù)字平臺測量巖石物理屬性參數(shù)可靠信。

    以三維流線展示孔隙空間流體流動,結(jié)合流域壓力梯度場,使模擬結(jié)果清晰直觀,此優(yōu)勢使數(shù)字滲流實(shí)驗(yàn)便利于傳統(tǒng)巖石物理實(shí)驗(yàn)。

    基于數(shù)字巖心研究巖石物理屬性最大優(yōu)勢是,一旦與真實(shí)巖心等價的孔隙空間及巖石骨架構(gòu)建完畢,結(jié)合不同的模擬軟件,便可實(shí)現(xiàn)同一巖心孔隙空間不同類型物理場的模擬與耦合。例如,基于孔隙空間模型,我們可以進(jìn)一步研究孔隙空間內(nèi)水動力擴(kuò)散、化學(xué)藥品的輸運(yùn)及其與孔隙壁面的反應(yīng)、多相流體流動等。

    4 展望

    數(shù)字巖心技術(shù)是當(dāng)下研究微觀孔喉性質(zhì)對宏觀巖石物理屬性影響規(guī)律的有效手段,深入研究對精細(xì)油藏描述及孔隙尺度提高油氣采收率具有重要意義。結(jié)合本次對低滲砂巖滲流特征的研究工作,筆者認(rèn)為數(shù)字巖心技術(shù)未來的研究應(yīng)著重于以下三點(diǎn):

    1)結(jié)合多種技術(shù)與算法多尺度,構(gòu)建非均質(zhì)復(fù)雜儲層的三維數(shù)字巖心模型。

    2)完善能夠同時考慮流體與巖石壁面、微觀基質(zhì)、粘土礦物等物化反應(yīng)更加精確的滲流模型。

    3)改進(jìn)算法或計算程序,在現(xiàn)有計算機(jī)硬件平臺上,擴(kuò)大計算模型尺度的同時提升計算速度,使數(shù)字巖石物理技術(shù)得到更廣泛的應(yīng)用。

    [1] 潘凌, 方全堂, 段永剛. 低滲油藏非均質(zhì)性對采收率的影響因素研究[J]. 西南石油大學(xué)學(xué)報(自然科學(xué)版), 2012(03):111-115. PAN L, FANG Q T, DUAN Y G. Mechanism Research of the Impact of Heterogeneity on the Low Permeability Reservoir Recovery[J]. Journal of Southwest Petroleum University(Science & Technology Edition),2012(03):111-115. (In Chinese)

    [2] 李濤. 特低滲油藏水驅(qū)油規(guī)律研究[D]. 青島:中國石油大學(xué), 2008. LI T. The Study of Water-oil Displacement Law in Extra-low Permeability Reservoir[D]. Qingdao:China University of Petroleum (EastChina), 2008. (In Chinese)

    [3] 郭平, 苑志旺, 易敏, 等. 低滲低壓油藏真實(shí)巖心薄片微觀水驅(qū)試驗(yàn)研究[J]. 石油天然氣學(xué)報, 2009(04):100-105. GUO P, YUAN Z W, YI M, et al. Experimental Study of Microscopic Water Displacement in Low-permeability and Low-pressure Reservoirs by Using Sandstone Slice[J]. Journal of Oil and Gas Technology, 2009(04):100-105. (In Chinese)

    [4] 孫致學(xué), 楚鵬, 張青茹, 等. 微觀驅(qū)油實(shí)驗(yàn)數(shù)值模擬方法[J]. 大慶石油地質(zhì)與開發(fā), 2015(03):124-128. SUN Z X, CHU P, ZHANG Q R , et al. Numerical simulating method of the microscopic oil-flooding experiment[J]. Petroleum Geology and Oilfield Development in Daqing, 2015(03): 124-128. (In Chinese)

    [5] 何文祥, 楊億前, 馬超亞. 特低滲透率儲層水驅(qū)油規(guī)律實(shí)驗(yàn)研究[J]. 巖性油氣藏, 2010(04):109-111. HE W X, YANG Y Q, MA C X. Experimental study on water flooding in ultra-low permeability reservoirs[J]. lithologic reservoirs, 2010(04):109 -111. (In Chinese)

    [6] 高輝, 解偉, 楊建鵬, 等. 基于恒速壓汞技術(shù)的特低—超低滲砂巖儲層微觀孔喉特征[J]. 石油實(shí)驗(yàn)地質(zhì), 2011(02):206-211. GAO H, XIE W, YANG J P, et al. Pore throat characteristics of extra-ultra low permeability sandstone reservoir based on constant-rate mercury penetration technique[J].Petroleum geology & experiment, 2011(02):206-211. (In Chinese)

    [7] 李海波, 郭和坤, 劉強(qiáng), 等. 致密油儲層水驅(qū)油核磁共振實(shí)驗(yàn)研究[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2014(12):4370-4376. LI H B, GUO H K, LIU Q, et al. NMR experimental study of water displacing oil of tight oil reservoir[J], Journal of Central South University (Science and Technology) 2014(12):4370-4376. (In Chinese)

    [8] 孔強(qiáng)夫, 周燦燦, 李潮流, 等. 數(shù)字巖心電性數(shù)值模擬方法及其發(fā)展方向[J]. 中國石油勘探, 2015(01):69-77. KONG Q F, ZHOU C C, LI C L, et al. Numerical simulation method of digital core electrical property and its development orientations[J]. China petroleum exploration, 2015(01):69-77. (In Chinese)

    [9] 徐湖山. 基于CT技術(shù)的數(shù)字巖心重構(gòu)及其應(yīng)用研究[D]. 大連:大連理工大學(xué), 2014. XU H S. Study on reconstruction technology and application of digital core based on the X-ray computed tomography (CT)[D]. Dalian:Dalian University of Technology, 2014. (In Chinese)

    [10]閆國亮. 基于數(shù)字巖心儲層滲透率模型研究[D].青島: 中國石油大學(xué)(華東), 2013. YAN G L. Research of Permeability Models of Reservoirs Based on Digital Cores[D]. Qingdao:University of Petroleum (EastChina), 2013. (In Chinese)

    [11]王冬欣. 基于Micro-CT圖像的數(shù)字巖心孔隙級網(wǎng)絡(luò)建模研究[D].長春:吉林大學(xué), 2015. WANG D X. The Research of Digital Core Network Extraction Based on Micro-CT Images[D].Changchuen:JILIN University, 2015. (In Chinese)

    [12]ANDR·H, COMBARET N, DVORKIN J, et al. Digital rock physics benchmarks—part II: Computing effective properties[J]. Computers & Geosciences, 2013,50(1):33-43.

    [13]朱洪林. 低滲砂巖儲層孔隙結(jié)構(gòu)表征及應(yīng)用研究[D]. 成都:西南石油大學(xué), 2014. ZHU H L. Pore Structure characterization of low permeability sandstone reservoir and application[D].Chengdu:Southwest Petroleum University, 2014. (In Chinese)

    [14]ZHANG S, KLIMENTIDIS R, CO E M U. Porosity and permeability analysis on nanoscale FIB-SEM 3D imaging of shale rock: Society of Core Analysts, 2011[C]. The International Symposium of the Society of Core Analysts held in Austin, Texas, USA 18-21 September, 2011. SCA2011-30

    [15]張磊, 姚軍, 孫海, 等. 基于數(shù)字巖心技術(shù)的氣體解析/擴(kuò)散格子Boltzmann模擬[J]. 石油學(xué)報, 2015(03):361-365. ZHANG L, YAO J, SUN H, et al. Lattice Boltzmann simulation of gas desorption and diffusion based on digital core technology[J].Acta petrolei sinica, 2015(03):361-365. (In Chinese)

    [16]EDO S, VENTUROLI, MADDALENA. Lattice-Boltzmann studies of fluid flow in porous media with realistic rock geometries[J]. Computers & Mathematics with Applications, 2010,59(7):2305-2314.

    Study on seepage characteristics of low permeability sandstone core using digital core technology

    WANG Chunsheng, LIU Yang, SUN Qiji, ZHANG Kai, LIANG Chao

    (Northeast Petroleum University, Institute of Petroleum Engineering,Daqing 163318, China)

    When we adopt the traditional rock physics experiments to study the macroscopic petrophysical properties of low porosity, low permeability and high complexity pore structure reservoirs, there will be such disadvantages as high cost, long cycle, large error of image, as well as the difficulty on how to quantitatively study the microscopic factors in reservoir influence rock physical properties. Through Micro-CT scanning and advanced imaging processing technology, an accurate representative pore space model equivalent to true core is constructed. In addition, we use numerical simulation of Avizo-XLab solver to calculate the flow of pore fluid with N-S as control equation and adopt Darcy's to calculate the parameter of rock seepage. Through the visualization of the results, the pore-space fluid flow is displayed in a three-dimensional streamline, and the pore pressure distribution field is displayed with a high gradient, all these make the simulation result clear and intuitive. The research shows that the numerical simulation based on the digital platform has a high credibility in the micro flow of pore fluid and the measurement of rock seepage parameters. This paper makes up the deficiency of the traditional rock physics experiment, and providing a new method for digital research on rock physical properties.

    traditional rock physics experiments; microscopic factors; Micro-CT scanning; the micro flow of pore fluid; rock seepage parameters

    2016-11-30 改回日期:2017-02-03

    黑龍江省創(chuàng)新訓(xùn)練項(xiàng)目(201610220047);東北石油大學(xué)優(yōu)秀科研人才培育基金(SCXHB201602)

    王春生(1977-),男,副教授,主要從事復(fù)雜流體流動數(shù)值模擬研究工作,E-mail:wangchunsheng@nepu.edu.cn。

    1001-1749(2017)04-0573-06

    TE135+.4

    A

    10.3969/j.issn.1001-1749.2017.04.19

    猜你喜歡
    巖心滲流巖石
    第五章 巖石小專家
    Osteotomized folded scapular tip free flap for complex midfacial reconstruction
    3深源巖石
    一種叫做煤炭的巖石
    海藻與巖石之間
    一種頁巖巖心資料的保存方法
    化工管理(2017年23期)2017-09-11 14:14:22
    Acellular allogeneic nerve grafting combined with bone marrow mesenchymal stem cell transplantation for the repair of long-segment sciatic nerve defects: biomechanics and validation of mathematical models
    長巖心注CO2氣水交替驅(qū)試驗(yàn)?zāi)M研究
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    一个人看视频在线观看www免费| 欧美人与善性xxx| 国产成人a区在线观看| 亚洲精华国产精华液的使用体验 | 国产成人freesex在线 | 亚洲图色成人| 亚洲精品一卡2卡三卡4卡5卡| 精品人妻视频免费看| 亚洲av成人av| 成人性生交大片免费视频hd| 国产激情偷乱视频一区二区| 亚洲精品亚洲一区二区| 有码 亚洲区| 久久九九热精品免费| 国产av一区在线观看免费| av专区在线播放| 成人高潮视频无遮挡免费网站| 精品久久久久久久人妻蜜臀av| 免费在线观看成人毛片| 成人性生交大片免费视频hd| 久久久国产成人精品二区| 欧美极品一区二区三区四区| 久久久久精品国产欧美久久久| 国产高清三级在线| 亚洲中文日韩欧美视频| 日韩国内少妇激情av| 亚洲av五月六月丁香网| 在线观看66精品国产| 三级经典国产精品| 精品人妻视频免费看| 亚洲欧美日韩无卡精品| 成人毛片a级毛片在线播放| 插逼视频在线观看| av福利片在线观看| 国产精品久久久久久久电影| 免费观看在线日韩| 一区二区三区高清视频在线| 国产爱豆传媒在线观看| 91午夜精品亚洲一区二区三区| 国产精品久久久久久亚洲av鲁大| 亚洲欧美日韩无卡精品| 婷婷六月久久综合丁香| 亚洲欧美成人综合另类久久久 | 久久99热6这里只有精品| 久久草成人影院| 女的被弄到高潮叫床怎么办| 亚洲人成网站在线播| 蜜桃久久精品国产亚洲av| 成人亚洲精品av一区二区| 精品99又大又爽又粗少妇毛片| 夜夜夜夜夜久久久久| 成人鲁丝片一二三区免费| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成网站在线播| 国内精品久久久久精免费| 老师上课跳d突然被开到最大视频| 黄色欧美视频在线观看| a级一级毛片免费在线观看| 少妇的逼水好多| 午夜视频国产福利| 真人做人爱边吃奶动态| 日日摸夜夜添夜夜添av毛片| 亚洲人成网站在线观看播放| 欧美日本视频| 亚洲精品一卡2卡三卡4卡5卡| 黄色日韩在线| 国产乱人偷精品视频| 欧美最黄视频在线播放免费| 国产白丝娇喘喷水9色精品| 国产精品亚洲美女久久久| 亚洲欧美成人精品一区二区| 欧美不卡视频在线免费观看| 久久久久国产网址| 国产在线男女| 成人特级黄色片久久久久久久| 日产精品乱码卡一卡2卡三| 久久中文看片网| 嫩草影视91久久| 老女人水多毛片| 日本一本二区三区精品| 国产乱人偷精品视频| 别揉我奶头 嗯啊视频| 久久久久久久久久久丰满| 欧美高清性xxxxhd video| 欧美极品一区二区三区四区| 直男gayav资源| 亚洲丝袜综合中文字幕| 日韩成人伦理影院| 看免费成人av毛片| 成年版毛片免费区| 夜夜看夜夜爽夜夜摸| 国产爱豆传媒在线观看| 色在线成人网| 久久精品国产亚洲av香蕉五月| 精品一区二区免费观看| 日日撸夜夜添| 91麻豆精品激情在线观看国产| 免费观看人在逋| 日韩欧美精品免费久久| 免费搜索国产男女视频| 日本免费a在线| 午夜福利在线观看吧| 亚洲精品影视一区二区三区av| 中国美女看黄片| 一区二区三区免费毛片| 久久久久国内视频| 性色avwww在线观看| 毛片女人毛片| 国产美女午夜福利| 久久精品夜色国产| 国产亚洲欧美98| 特大巨黑吊av在线直播| 成年av动漫网址| 麻豆乱淫一区二区| 成人鲁丝片一二三区免费| 国产精品伦人一区二区| 99久国产av精品| 国产精品永久免费网站| 国产乱人偷精品视频| 久久久久久久亚洲中文字幕| 国产精品电影一区二区三区| 直男gayav资源| 亚洲在线自拍视频| 国产高清三级在线| 亚洲人成网站高清观看| 大又大粗又爽又黄少妇毛片口| 国产又黄又爽又无遮挡在线| 午夜精品一区二区三区免费看| 精品久久久久久久久久免费视频| 3wmmmm亚洲av在线观看| 亚洲成人av在线免费| 国产精品福利在线免费观看| 精品一区二区免费观看| 国产成人freesex在线 | 国产淫片久久久久久久久| 国产大屁股一区二区在线视频| 成年版毛片免费区| 一a级毛片在线观看| 欧美性猛交黑人性爽| 可以在线观看的亚洲视频| 天堂av国产一区二区熟女人妻| 日韩一本色道免费dvd| 狂野欧美激情性xxxx在线观看| 人妻少妇偷人精品九色| 中文字幕久久专区| 精品一区二区三区视频在线观看免费| 99热只有精品国产| 国产精品久久久久久av不卡| 精品久久久久久久久久久久久| av天堂中文字幕网| 久久精品国产清高在天天线| 亚洲精品久久国产高清桃花| 亚洲天堂国产精品一区在线| 国产午夜福利久久久久久| 婷婷六月久久综合丁香| 18禁在线无遮挡免费观看视频 | 国产精品久久久久久亚洲av鲁大| 国产精品三级大全| 欧美xxxx性猛交bbbb| 久久久a久久爽久久v久久| 国产成人一区二区在线| 成人三级黄色视频| 成人午夜高清在线视频| 精品99又大又爽又粗少妇毛片| 国产毛片a区久久久久| 欧美bdsm另类| 一区二区三区四区激情视频 | 成人美女网站在线观看视频| 欧美一区二区精品小视频在线| 久久久久久国产a免费观看| 99热6这里只有精品| 日韩成人伦理影院| 亚洲精品一区av在线观看| 观看美女的网站| 国产精品久久久久久久久免| 一个人看视频在线观看www免费| 亚洲精品成人久久久久久| 免费人成视频x8x8入口观看| 中文字幕人妻熟人妻熟丝袜美| 女同久久另类99精品国产91| 国产一区亚洲一区在线观看| 男人和女人高潮做爰伦理| 亚洲在线自拍视频| 老熟妇仑乱视频hdxx| 久久精品综合一区二区三区| 在线看三级毛片| 国产午夜精品论理片| 97人妻精品一区二区三区麻豆| 91在线精品国自产拍蜜月| 晚上一个人看的免费电影| 内射极品少妇av片p| 国产三级在线视频| 午夜亚洲福利在线播放| 色视频www国产| 小蜜桃在线观看免费完整版高清| 伦理电影大哥的女人| 成人性生交大片免费视频hd| 午夜精品国产一区二区电影 | 波多野结衣巨乳人妻| 看免费成人av毛片| 夜夜看夜夜爽夜夜摸| 欧美日韩在线观看h| 亚洲美女视频黄频| 久久久精品94久久精品| 日本爱情动作片www.在线观看 | 欧美日韩国产亚洲二区| 欧美最黄视频在线播放免费| 天堂影院成人在线观看| 插阴视频在线观看视频| 国产一区二区激情短视频| 国产精品综合久久久久久久免费| 天堂av国产一区二区熟女人妻| 欧美成人免费av一区二区三区| 非洲黑人性xxxx精品又粗又长| 内地一区二区视频在线| 亚洲色图av天堂| 日本五十路高清| 国产精品一及| 国产精品爽爽va在线观看网站| 久久6这里有精品| 免费搜索国产男女视频| 色播亚洲综合网| 日韩高清综合在线| 日韩 亚洲 欧美在线| 老司机午夜福利在线观看视频| 精品久久久噜噜| 国产精品99久久久久久久久| 欧美区成人在线视频| av天堂在线播放| 久久精品国产亚洲av香蕉五月| 久久久久久国产a免费观看| 搡老妇女老女人老熟妇| 亚洲国产色片| 精品99又大又爽又粗少妇毛片| av黄色大香蕉| 少妇的逼水好多| 精品少妇黑人巨大在线播放 | 天堂动漫精品| 给我免费播放毛片高清在线观看| 国产精品一区二区免费欧美| 熟女电影av网| 深夜a级毛片| 亚洲国产欧洲综合997久久,| 嫩草影院新地址| 麻豆成人午夜福利视频| 女的被弄到高潮叫床怎么办| 免费人成在线观看视频色| 人妻制服诱惑在线中文字幕| 亚洲一区二区三区色噜噜| 亚洲国产精品sss在线观看| 亚洲国产欧洲综合997久久,| 最近中文字幕高清免费大全6| 亚洲国产色片| 亚洲人成网站在线播| 国内精品久久久久精免费| 中文资源天堂在线| 女同久久另类99精品国产91| 1000部很黄的大片| 亚洲精品在线观看二区| 在线观看一区二区三区| 男女边吃奶边做爰视频| 久久6这里有精品| 给我免费播放毛片高清在线观看| 久久草成人影院| 天堂网av新在线| 久久久久国产精品人妻aⅴ院| 国产一级毛片七仙女欲春2| 男人舔奶头视频| avwww免费| 露出奶头的视频| 欧美成人一区二区免费高清观看| 成人性生交大片免费视频hd| 精品午夜福利视频在线观看一区| 男女那种视频在线观看| 91在线观看av| 国产午夜福利久久久久久| 精品人妻视频免费看| 精品久久久久久久久av| 亚洲精品久久国产高清桃花| 韩国av在线不卡| 国内久久婷婷六月综合欲色啪| 亚洲国产日韩欧美精品在线观看| 日韩人妻高清精品专区| 99久久中文字幕三级久久日本| 69av精品久久久久久| 日本撒尿小便嘘嘘汇集6| 成人一区二区视频在线观看| 亚洲国产高清在线一区二区三| 少妇被粗大猛烈的视频| 一级毛片久久久久久久久女| 别揉我奶头~嗯~啊~动态视频| 国产在视频线在精品| 亚洲av五月六月丁香网| 日本与韩国留学比较| 欧美一区二区国产精品久久精品| 久久九九热精品免费| 乱码一卡2卡4卡精品| 国产精品一区二区三区四区免费观看 | 特级一级黄色大片| 1000部很黄的大片| 免费不卡的大黄色大毛片视频在线观看 | 欧美一区二区精品小视频在线| 久久婷婷人人爽人人干人人爱| 日韩欧美免费精品| 欧美色欧美亚洲另类二区| 国产在线男女| 三级国产精品欧美在线观看| av免费在线看不卡| av卡一久久| 老司机福利观看| 天天一区二区日本电影三级| 大香蕉久久网| 日本五十路高清| 国产一区二区三区在线臀色熟女| 在线免费十八禁| 亚洲在线自拍视频| 国产午夜福利久久久久久| 久久久久久久久久黄片| 俺也久久电影网| 99久久成人亚洲精品观看| 一个人看的www免费观看视频| 男人舔女人下体高潮全视频| 丝袜美腿在线中文| 欧美成人精品欧美一级黄| 日本熟妇午夜| 美女免费视频网站| 搡老妇女老女人老熟妇| 国产单亲对白刺激| 国产免费男女视频| 亚洲精品一区av在线观看| 丰满乱子伦码专区| 国产人妻一区二区三区在| 一级毛片aaaaaa免费看小| 国产高潮美女av| 国产欧美日韩精品亚洲av| 久久久久九九精品影院| 免费在线观看成人毛片| 内地一区二区视频在线| 欧美+日韩+精品| 日本五十路高清| 男人舔女人下体高潮全视频| 国国产精品蜜臀av免费| 亚洲av第一区精品v没综合| 久久国内精品自在自线图片| 国产在视频线在精品| 亚洲色图av天堂| 国产不卡一卡二| 搡老妇女老女人老熟妇| 波多野结衣巨乳人妻| 两个人的视频大全免费| 国产黄a三级三级三级人| 久久久国产成人免费| 啦啦啦观看免费观看视频高清| 亚洲精品在线观看二区| av视频在线观看入口| 精品久久久久久久久久免费视频| 欧美+亚洲+日韩+国产| 一级毛片电影观看 | 中文字幕人妻熟人妻熟丝袜美| 变态另类成人亚洲欧美熟女| 免费av毛片视频| 成人无遮挡网站| 国产综合懂色| 国产极品精品免费视频能看的| 天堂动漫精品| 国产美女午夜福利| 赤兔流量卡办理| 免费看a级黄色片| 国产黄色小视频在线观看| 精品日产1卡2卡| 日韩精品中文字幕看吧| 久久精品夜夜夜夜夜久久蜜豆| 国产成人a∨麻豆精品| 在线播放国产精品三级| 日本成人三级电影网站| 日韩欧美三级三区| 亚洲欧美日韩高清专用| 少妇的逼水好多| 村上凉子中文字幕在线| 91在线精品国自产拍蜜月| 亚洲不卡免费看| 欧美最新免费一区二区三区| 日韩欧美三级三区| 1024手机看黄色片| 在线国产一区二区在线| 国产午夜精品久久久久久一区二区三区 | 99在线人妻在线中文字幕| 女的被弄到高潮叫床怎么办| 变态另类成人亚洲欧美熟女| 精品午夜福利在线看| 一级黄片播放器| 熟妇人妻久久中文字幕3abv| 少妇人妻一区二区三区视频| 久久6这里有精品| 又爽又黄a免费视频| 色哟哟哟哟哟哟| 舔av片在线| 成熟少妇高潮喷水视频| 亚洲电影在线观看av| 亚洲熟妇熟女久久| 午夜激情福利司机影院| 亚洲欧美日韩东京热| 久久精品国产亚洲av香蕉五月| 亚洲va在线va天堂va国产| 99热精品在线国产| 久久久午夜欧美精品| 校园人妻丝袜中文字幕| 欧美成人免费av一区二区三区| 老熟妇乱子伦视频在线观看| 99热只有精品国产| 久久精品夜色国产| 搞女人的毛片| 日韩欧美一区二区三区在线观看| 99久久精品国产国产毛片| 天堂av国产一区二区熟女人妻| 国产一区二区三区av在线 | 97在线视频观看| 99久国产av精品| 欧美日韩综合久久久久久| 国产精品一区二区三区四区久久| 一进一出抽搐动态| 欧美+日韩+精品| 国内揄拍国产精品人妻在线| 国产精品嫩草影院av在线观看| 日日摸夜夜添夜夜添小说| 最好的美女福利视频网| 欧美成人一区二区免费高清观看| 一级毛片aaaaaa免费看小| 午夜福利18| 97热精品久久久久久| 国产探花在线观看一区二区| 村上凉子中文字幕在线| 亚洲精品影视一区二区三区av| 亚洲在线观看片| 国产精品久久久久久亚洲av鲁大| 成人午夜高清在线视频| 中文字幕久久专区| 亚洲欧美日韩无卡精品| 成人无遮挡网站| 欧美高清成人免费视频www| 国产一区二区在线观看日韩| 亚洲最大成人手机在线| 亚洲欧美日韩高清在线视频| 久久欧美精品欧美久久欧美| 亚洲aⅴ乱码一区二区在线播放| 大香蕉久久网| 色噜噜av男人的天堂激情| 禁无遮挡网站| 搡老熟女国产l中国老女人| 国产伦精品一区二区三区视频9| 久久九九热精品免费| 国产伦一二天堂av在线观看| 99热网站在线观看| 精品无人区乱码1区二区| 亚洲经典国产精华液单| 欧美一区二区国产精品久久精品| 校园春色视频在线观看| 1000部很黄的大片| 日日干狠狠操夜夜爽| 淫妇啪啪啪对白视频| 国产精品福利在线免费观看| 97热精品久久久久久| 中文字幕av在线有码专区| 一区福利在线观看| 一进一出抽搐动态| 婷婷精品国产亚洲av| 国产av不卡久久| 又爽又黄无遮挡网站| 精品一区二区三区视频在线观看免费| 美女xxoo啪啪120秒动态图| 国产色爽女视频免费观看| 亚洲成人中文字幕在线播放| 男插女下体视频免费在线播放| 日日摸夜夜添夜夜添av毛片| 免费电影在线观看免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲av中文av极速乱| а√天堂www在线а√下载| 99久久成人亚洲精品观看| 久久久久久九九精品二区国产| 一区二区三区四区激情视频 | av黄色大香蕉| 成年女人看的毛片在线观看| 免费看美女性在线毛片视频| 精品久久久久久久久久免费视频| 少妇人妻一区二区三区视频| 欧美色视频一区免费| 少妇熟女aⅴ在线视频| 97在线视频观看| 亚洲国产精品国产精品| 村上凉子中文字幕在线| 日本欧美国产在线视频| 不卡视频在线观看欧美| 色哟哟·www| 亚洲不卡免费看| 免费看a级黄色片| 一本一本综合久久| 国产探花极品一区二区| 亚洲精品粉嫩美女一区| 欧美另类亚洲清纯唯美| 免费av不卡在线播放| 欧美性猛交黑人性爽| 尤物成人国产欧美一区二区三区| 久久久午夜欧美精品| 毛片女人毛片| 日韩,欧美,国产一区二区三区 | 一本精品99久久精品77| 村上凉子中文字幕在线| av女优亚洲男人天堂| 亚洲经典国产精华液单| 18禁在线播放成人免费| 国内精品久久久久精免费| 精品国产三级普通话版| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕一区二区三区有码在线看| 嫩草影院新地址| 成人精品一区二区免费| 99热网站在线观看| 极品教师在线视频| 大型黄色视频在线免费观看| 男女视频在线观看网站免费| 直男gayav资源| 欧美日韩综合久久久久久| 亚洲激情五月婷婷啪啪| 久久人人精品亚洲av| 少妇丰满av| 精品乱码久久久久久99久播| 国国产精品蜜臀av免费| 老司机影院成人| 亚洲专区国产一区二区| 亚洲内射少妇av| 亚洲久久久久久中文字幕| 日本五十路高清| 国产亚洲精品久久久com| 亚洲av五月六月丁香网| 偷拍熟女少妇极品色| 成人一区二区视频在线观看| 国产精品,欧美在线| 国产91av在线免费观看| 亚洲精品日韩在线中文字幕 | 嫩草影院精品99| 亚洲美女搞黄在线观看 | 欧美一区二区国产精品久久精品| 亚洲天堂国产精品一区在线| avwww免费| 成人二区视频| 亚洲av中文av极速乱| 1024手机看黄色片| 插逼视频在线观看| 18禁在线无遮挡免费观看视频 | 久久久精品94久久精品| 中文字幕久久专区| 日韩欧美三级三区| 色尼玛亚洲综合影院| 在线观看一区二区三区| 欧美丝袜亚洲另类| 三级国产精品欧美在线观看| 俺也久久电影网| 亚洲经典国产精华液单| 两个人视频免费观看高清| 国产精品亚洲美女久久久| 成人特级av手机在线观看| 丝袜喷水一区| 日日撸夜夜添| 日本三级黄在线观看| 少妇丰满av| 99热这里只有精品一区| www日本黄色视频网| 国产单亲对白刺激| 此物有八面人人有两片| 日韩 亚洲 欧美在线| 亚洲成人久久性| 精品欧美国产一区二区三| 成年女人看的毛片在线观看| 男人的好看免费观看在线视频| 日本熟妇午夜| 亚洲av熟女| 伦理电影大哥的女人| 卡戴珊不雅视频在线播放| 久久鲁丝午夜福利片| 欧美激情在线99| 日韩,欧美,国产一区二区三区 | 黄色配什么色好看| 日本免费a在线| 亚洲,欧美,日韩| 日本爱情动作片www.在线观看 | 乱人视频在线观看| 搡老熟女国产l中国老女人| 欧美一区二区亚洲| 久久久久久久久久成人| 欧美日韩乱码在线| 两性午夜刺激爽爽歪歪视频在线观看| 综合色丁香网| 亚洲国产精品合色在线| 色播亚洲综合网| 变态另类成人亚洲欧美熟女| 日韩欧美精品免费久久| 久久精品91蜜桃| 午夜爱爱视频在线播放| 日韩精品青青久久久久久| 国产中年淑女户外野战色| 免费看a级黄色片| 国产69精品久久久久777片| 日本色播在线视频| 色av中文字幕| 男女做爰动态图高潮gif福利片| 99久久成人亚洲精品观看| 国产精品,欧美在线| 午夜精品一区二区三区免费看| 久久精品影院6| 国产在线男女| 国产精品精品国产色婷婷| 国产探花极品一区二区| 免费人成视频x8x8入口观看| 国产老妇女一区|