• <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
    av有码第一页| 久久伊人香网站| 国产av精品麻豆| 欧美黄色片欧美黄色片| 国产在线观看jvid| 欧美精品亚洲一区二区| 精品一区二区三区av网在线观看| www.熟女人妻精品国产| 精品第一国产精品| 夜夜躁狠狠躁天天躁| 男女高潮啪啪啪动态图| 欧美另类亚洲清纯唯美| 国产欧美日韩综合在线一区二区| 看片在线看免费视频| 欧美乱妇无乱码| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区精品| 国产精品成人在线| 老鸭窝网址在线观看| 91av网站免费观看| 国产伦人伦偷精品视频| 黑丝袜美女国产一区| 精品一区二区三卡| 亚洲全国av大片| 国产精品综合久久久久久久免费 | 亚洲av日韩精品久久久久久密| 国产又爽黄色视频| 国产一区在线观看成人免费| 涩涩av久久男人的天堂| 男人舔女人的私密视频| 日韩中文字幕欧美一区二区| 色老头精品视频在线观看| 在线国产一区二区在线| 精品一区二区三区av网在线观看| 精品高清国产在线一区| 真人做人爱边吃奶动态| 黄色怎么调成土黄色| 欧美日韩亚洲国产一区二区在线观看| x7x7x7水蜜桃| 国产欧美日韩一区二区精品| 男人的好看免费观看在线视频 | 国产有黄有色有爽视频| 丰满的人妻完整版| 久久久久久免费高清国产稀缺| 日韩国内少妇激情av| 久久人妻av系列| 热re99久久精品国产66热6| 又黄又爽又免费观看的视频| 亚洲国产精品合色在线| 午夜福利影视在线免费观看| 又黄又粗又硬又大视频| 高清欧美精品videossex| 一本大道久久a久久精品| 国产精品久久久久成人av| 91九色精品人成在线观看| 国产成人免费无遮挡视频| 91老司机精品| 9色porny在线观看| 搡老乐熟女国产| 国产人伦9x9x在线观看| 成人三级做爰电影| 首页视频小说图片口味搜索| bbb黄色大片| 亚洲欧美激情综合另类| 日韩一卡2卡3卡4卡2021年| 无人区码免费观看不卡| 别揉我奶头~嗯~啊~动态视频| 麻豆av在线久日| 99久久综合精品五月天人人| 久热这里只有精品99| 精品一区二区三卡| 亚洲视频免费观看视频| 午夜影院日韩av| 精品电影一区二区在线| 亚洲 国产 在线| 国产三级在线视频| 精品一区二区三卡| 国产av一区二区精品久久| a级片在线免费高清观看视频| 久久久水蜜桃国产精品网| 国产成人精品在线电影| 男女下面插进去视频免费观看| tocl精华| 中文字幕另类日韩欧美亚洲嫩草| 精品国产超薄肉色丝袜足j| 乱人伦中国视频| 国产xxxxx性猛交| 国产精品一区二区精品视频观看| av网站免费在线观看视频| 最近最新免费中文字幕在线| 国产伦人伦偷精品视频| 男女之事视频高清在线观看| 国产精品自产拍在线观看55亚洲| av在线天堂中文字幕 | 久久精品国产亚洲av高清一级| 一级,二级,三级黄色视频| 夜夜爽天天搞| 视频在线观看一区二区三区| 一边摸一边抽搐一进一小说| 国产色视频综合| 嫁个100分男人电影在线观看| 麻豆国产av国片精品| 亚洲欧美激情综合另类| 精品欧美一区二区三区在线| 日韩欧美在线二视频| 制服人妻中文乱码| 久久热在线av| 9色porny在线观看| 黄色毛片三级朝国网站| 午夜日韩欧美国产| 国产乱人伦免费视频| 动漫黄色视频在线观看| 18禁观看日本| 亚洲第一av免费看| 国产精品久久久人人做人人爽| 深夜精品福利| 亚洲av日韩精品久久久久久密| 久久人妻福利社区极品人妻图片| 色精品久久人妻99蜜桃| 99热国产这里只有精品6| 精品久久久久久电影网| 97超级碰碰碰精品色视频在线观看| 91麻豆av在线| 免费久久久久久久精品成人欧美视频| 久久久久久久久免费视频了| 脱女人内裤的视频| 99国产精品一区二区三区| 亚洲自偷自拍图片 自拍| 国产视频一区二区在线看| 丝袜在线中文字幕| 新久久久久国产一级毛片| av视频免费观看在线观看| 亚洲熟妇中文字幕五十中出 | 久久人妻av系列| 国产1区2区3区精品| 一本大道久久a久久精品| 国产精品98久久久久久宅男小说| 色哟哟哟哟哟哟| 国产真人三级小视频在线观看| 午夜福利,免费看| 色综合站精品国产| 狠狠狠狠99中文字幕| a在线观看视频网站| 在线观看www视频免费| 狂野欧美激情性xxxx| 亚洲免费av在线视频| 欧美成人免费av一区二区三区| 男女高潮啪啪啪动态图| 免费一级毛片在线播放高清视频 | 久久人人爽av亚洲精品天堂| 丝袜美足系列| 老司机在亚洲福利影院| 99久久久亚洲精品蜜臀av| 亚洲国产毛片av蜜桃av| 一进一出好大好爽视频| 国产成人精品在线电影| 久久人妻熟女aⅴ| 91国产中文字幕| 人人妻人人澡人人看| 成熟少妇高潮喷水视频| 热re99久久国产66热| 国产亚洲精品一区二区www| 在线观看www视频免费| 无遮挡黄片免费观看| 天天影视国产精品| 亚洲av成人av| 欧美成人性av电影在线观看| 后天国语完整版免费观看| 最新在线观看一区二区三区| 欧美久久黑人一区二区| 一本综合久久免费| 丝袜在线中文字幕| 在线观看免费视频网站a站| 亚洲九九香蕉| 香蕉国产在线看| 亚洲一区高清亚洲精品| 久久国产精品人妻蜜桃| 欧美国产精品va在线观看不卡| 免费在线观看影片大全网站| 美女午夜性视频免费| 一二三四社区在线视频社区8| 中文字幕另类日韩欧美亚洲嫩草| 中文欧美无线码| 操美女的视频在线观看| 91老司机精品| 亚洲国产精品合色在线| 亚洲伊人色综图| 欧美老熟妇乱子伦牲交| 欧美日韩亚洲综合一区二区三区_| 十八禁人妻一区二区| aaaaa片日本免费| 亚洲国产中文字幕在线视频| 亚洲国产欧美网| 国产亚洲精品久久久久久毛片| 午夜影院日韩av| 在线免费观看的www视频| 香蕉丝袜av| 欧美色视频一区免费| 成人亚洲精品av一区二区 | 精品乱码久久久久久99久播| 天堂中文最新版在线下载| 日韩有码中文字幕| 欧美中文综合在线视频| 看片在线看免费视频| 黑人巨大精品欧美一区二区mp4| 韩国精品一区二区三区| 国产91精品成人一区二区三区| 99国产综合亚洲精品| 制服诱惑二区| 丰满人妻熟妇乱又伦精品不卡| 欧美 亚洲 国产 日韩一| 97碰自拍视频| 91精品三级在线观看| 一个人观看的视频www高清免费观看 | 亚洲精品一区av在线观看| 国产欧美日韩精品亚洲av| 电影成人av| 精品日产1卡2卡| 好看av亚洲va欧美ⅴa在| 亚洲 国产 在线| 色老头精品视频在线观看| 欧美大码av| 一级黄色大片毛片| 欧美黄色片欧美黄色片| 99久久国产精品久久久| 国产伦一二天堂av在线观看| 欧美在线黄色| 香蕉国产在线看| 亚洲午夜理论影院| 亚洲成人免费av在线播放| а√天堂www在线а√下载| 免费久久久久久久精品成人欧美视频| 又紧又爽又黄一区二区| 男女下面进入的视频免费午夜 | 国产精品野战在线观看 | 老司机福利观看| 亚洲三区欧美一区| 亚洲第一青青草原| www.www免费av| 无人区码免费观看不卡| 国产精品偷伦视频观看了| 人妻久久中文字幕网| 国产精品久久久久成人av| 亚洲国产精品sss在线观看 | 国产精品成人在线| 日日爽夜夜爽网站| 乱人伦中国视频| 黄色视频不卡| 久久精品国产99精品国产亚洲性色 | www.www免费av| 搡老乐熟女国产| 91精品三级在线观看| 久久久久国内视频| 午夜成年电影在线免费观看| 少妇 在线观看| 亚洲第一青青草原| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 国产乱人伦免费视频| 亚洲成人国产一区在线观看| 满18在线观看网站| 久久草成人影院| 久久精品国产综合久久久| 亚洲中文字幕日韩| 久久久久久大精品| 国产激情久久老熟女| 动漫黄色视频在线观看| 国产免费现黄频在线看| 性欧美人与动物交配| 9热在线视频观看99| 久久精品aⅴ一区二区三区四区| 天堂中文最新版在线下载| 欧美大码av| 神马国产精品三级电影在线观看 | 夜夜爽天天搞| 国产精品一区二区免费欧美| www.自偷自拍.com| 国产精品自产拍在线观看55亚洲| 久久久国产一区二区| 欧美成人午夜精品| 国产精品九九99| 久久久精品国产亚洲av高清涩受| 午夜91福利影院| 日本 av在线| 高清黄色对白视频在线免费看| 制服人妻中文乱码| 天天躁夜夜躁狠狠躁躁| 女性生殖器流出的白浆| 国产成人av教育| 91大片在线观看| 成人18禁高潮啪啪吃奶动态图| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品无人区| 久久人妻av系列| 99久久人妻综合| 久久亚洲真实| 午夜成年电影在线免费观看| 亚洲三区欧美一区| 国产一区二区三区综合在线观看| 女警被强在线播放| 欧美日韩国产mv在线观看视频| av天堂在线播放| 日日干狠狠操夜夜爽| 日韩欧美在线二视频| 久久香蕉精品热| 国产午夜精品久久久久久| 日本wwww免费看| 99国产精品一区二区三区| 亚洲激情在线av| 中文欧美无线码| 精品卡一卡二卡四卡免费| 在线观看66精品国产| 999精品在线视频| 亚洲av日韩精品久久久久久密| 日本三级黄在线观看| 正在播放国产对白刺激| 妹子高潮喷水视频| 男女之事视频高清在线观看| 级片在线观看| www.自偷自拍.com| 1024视频免费在线观看| 欧美日韩亚洲国产一区二区在线观看| 热re99久久精品国产66热6| 俄罗斯特黄特色一大片| √禁漫天堂资源中文www| 久久婷婷成人综合色麻豆| 黄色视频,在线免费观看| 中出人妻视频一区二区| 亚洲五月婷婷丁香| 免费在线观看黄色视频的| 99久久综合精品五月天人人| 国产一区二区三区视频了| 国产精品九九99| 欧美黄色淫秽网站| 亚洲欧美激情在线| 成年人免费黄色播放视频| 久久这里只有精品19| 又黄又粗又硬又大视频| 日本 av在线| 男人操女人黄网站| 欧美日韩国产mv在线观看视频| 两个人免费观看高清视频| 99香蕉大伊视频| 亚洲专区字幕在线| 国产一卡二卡三卡精品| 一进一出抽搐动态| 欧美激情久久久久久爽电影 | 中文亚洲av片在线观看爽| 琪琪午夜伦伦电影理论片6080| 在线av久久热| 777久久人妻少妇嫩草av网站| 国产黄a三级三级三级人| 中文字幕色久视频| 欧美精品亚洲一区二区| 午夜福利在线免费观看网站| 自拍欧美九色日韩亚洲蝌蚪91| 成人国产一区最新在线观看| 在线观看免费午夜福利视频| ponron亚洲| 欧美色视频一区免费| 久久精品91蜜桃| 亚洲一区二区三区欧美精品| 侵犯人妻中文字幕一二三四区| 在线看a的网站| 美国免费a级毛片| 亚洲人成77777在线视频| 国产成人免费无遮挡视频| 成年人免费黄色播放视频| 人人妻人人添人人爽欧美一区卜| 亚洲熟妇中文字幕五十中出 | 国产成人精品在线电影| netflix在线观看网站| 国产亚洲欧美在线一区二区| 99国产精品99久久久久| 久久久精品欧美日韩精品| 丁香六月欧美| 久久久久久久午夜电影 | 一夜夜www| 国产熟女xx| 天天添夜夜摸| 亚洲专区字幕在线| 欧美成人性av电影在线观看| 老熟妇乱子伦视频在线观看| 一级,二级,三级黄色视频| 国产精品99久久99久久久不卡| 啪啪无遮挡十八禁网站| 一级黄色大片毛片| 男人的好看免费观看在线视频 | 久久精品亚洲av国产电影网| 欧美精品亚洲一区二区| 亚洲人成77777在线视频| 天堂影院成人在线观看| 欧美黄色片欧美黄色片| 国产99久久九九免费精品| 国产日韩一区二区三区精品不卡| 涩涩av久久男人的天堂| 国产97色在线日韩免费| 久久久久亚洲av毛片大全| 老司机午夜十八禁免费视频| 亚洲欧美日韩无卡精品| 亚洲avbb在线观看| 两个人免费观看高清视频| 色尼玛亚洲综合影院| 丰满的人妻完整版| 淫妇啪啪啪对白视频| 亚洲第一青青草原| 国产99久久九九免费精品| 十八禁网站免费在线| 免费在线观看日本一区| 久久久久久久午夜电影 | 交换朋友夫妻互换小说| 97人妻天天添夜夜摸| a级毛片在线看网站| 一二三四在线观看免费中文在| 国产黄色免费在线视频| 久久久久久免费高清国产稀缺| 亚洲色图 男人天堂 中文字幕| 丰满饥渴人妻一区二区三| 在线观看免费日韩欧美大片| 久久久久久人人人人人| 国产麻豆69| 黑人巨大精品欧美一区二区蜜桃| 亚洲一区高清亚洲精品| av欧美777| 久久伊人香网站| 18禁国产床啪视频网站| 亚洲 欧美一区二区三区| 国产三级黄色录像| 色综合婷婷激情| 欧美日韩亚洲综合一区二区三区_| 高清欧美精品videossex| 高清黄色对白视频在线免费看| 精品福利观看| 国产精品偷伦视频观看了| 亚洲一区二区三区色噜噜 | 国产精品免费视频内射| 欧美乱妇无乱码| 纯流量卡能插随身wifi吗| 日本wwww免费看| 无遮挡黄片免费观看| 久久久国产一区二区| 日韩欧美一区二区三区在线观看| 欧美久久黑人一区二区| 色综合站精品国产| 免费在线观看视频国产中文字幕亚洲| a级毛片在线看网站| 天堂俺去俺来也www色官网| 欧美成狂野欧美在线观看| 女人被躁到高潮嗷嗷叫费观| e午夜精品久久久久久久| 啦啦啦免费观看视频1| 午夜成年电影在线免费观看| 中文字幕人妻丝袜一区二区| 中文字幕高清在线视频| 日本免费a在线| 老司机在亚洲福利影院| 欧美精品一区二区免费开放| 欧美激情极品国产一区二区三区| 久久国产亚洲av麻豆专区| av在线天堂中文字幕 | 亚洲中文av在线| 咕卡用的链子| 国产成人精品久久二区二区91| 淫妇啪啪啪对白视频| 91成人精品电影| av网站在线播放免费| 免费av毛片视频| 国产成人av激情在线播放| 亚洲av日韩精品久久久久久密| 美国免费a级毛片| 交换朋友夫妻互换小说| 国产一区二区三区视频了| 国产成人av教育| 在线视频色国产色| 国产精品永久免费网站| 成年人免费黄色播放视频| 国产熟女xx| 一个人观看的视频www高清免费观看 | 久久热在线av| 久久久久久大精品| 又黄又爽又免费观看的视频| 青草久久国产| 亚洲少妇的诱惑av| 日韩欧美在线二视频| 天堂动漫精品| 香蕉久久夜色| 麻豆av在线久日| 亚洲一区二区三区色噜噜 | 久久精品91蜜桃| 热99re8久久精品国产| 在线视频色国产色| 91大片在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 日本三级黄在线观看| 中文字幕人妻丝袜一区二区| 欧美乱码精品一区二区三区| 交换朋友夫妻互换小说| 老司机靠b影院| 最新在线观看一区二区三区| 他把我摸到了高潮在线观看| 日日摸夜夜添夜夜添小说| 欧美日韩乱码在线| 中文亚洲av片在线观看爽| 国产国语露脸激情在线看| 巨乳人妻的诱惑在线观看| 国产高清国产精品国产三级| 午夜福利免费观看在线| 成年版毛片免费区| 一级作爱视频免费观看| 亚洲第一av免费看| 淫妇啪啪啪对白视频| 国产精品av久久久久免费| 国产乱人伦免费视频| 69精品国产乱码久久久| 亚洲国产精品sss在线观看 | 啦啦啦免费观看视频1| 丰满迷人的少妇在线观看| 午夜福利在线免费观看网站| 国产精品1区2区在线观看.| 中文字幕人妻熟女乱码| 美女午夜性视频免费| 50天的宝宝边吃奶边哭怎么回事| 极品教师在线免费播放| 侵犯人妻中文字幕一二三四区| 一进一出抽搐动态| 精品一品国产午夜福利视频| 亚洲精品av麻豆狂野| 亚洲国产精品一区二区三区在线| 亚洲精品中文字幕在线视频| 可以免费在线观看a视频的电影网站| 欧美午夜高清在线| 亚洲午夜理论影院| 麻豆av在线久日| www.自偷自拍.com| 国内久久婷婷六月综合欲色啪| 国产一区二区激情短视频| 亚洲一区中文字幕在线| 午夜免费激情av| 国产男靠女视频免费网站| 成人18禁在线播放| 国产精品爽爽va在线观看网站 | 欧美一级毛片孕妇| 亚洲久久久国产精品| 又黄又粗又硬又大视频| 欧美成狂野欧美在线观看| 欧美av亚洲av综合av国产av| 国产aⅴ精品一区二区三区波| 色综合站精品国产| 一二三四在线观看免费中文在| 黄色片一级片一级黄色片| avwww免费| 窝窝影院91人妻| 国产成人精品无人区| 午夜视频精品福利| 亚洲男人的天堂狠狠| 日韩三级视频一区二区三区| 色老头精品视频在线观看| 中文字幕人妻熟女乱码| 在线观看免费视频网站a站| 亚洲精品av麻豆狂野| 啪啪无遮挡十八禁网站| 伊人久久大香线蕉亚洲五| 亚洲精品国产一区二区精华液| 国产av一区二区精品久久| 最近最新中文字幕大全电影3 | 变态另类成人亚洲欧美熟女 | 亚洲七黄色美女视频| 在线观看66精品国产| 老司机午夜十八禁免费视频| 亚洲精品一二三| 国产伦一二天堂av在线观看| 两性夫妻黄色片| 久久久国产欧美日韩av| 高潮久久久久久久久久久不卡| 黄频高清免费视频| 成人av一区二区三区在线看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久水蜜桃国产精品网| 欧美精品啪啪一区二区三区| av福利片在线| 国产av一区二区精品久久| 亚洲午夜理论影院| 国产精华一区二区三区| 天天影视国产精品| 黄色毛片三级朝国网站| 亚洲精华国产精华精| 亚洲专区中文字幕在线| 国产片内射在线| 国产有黄有色有爽视频| 999精品在线视频| 亚洲avbb在线观看| 91成年电影在线观看| 亚洲精品成人av观看孕妇| 欧美黄色淫秽网站| 国产一区二区三区综合在线观看| 欧美激情久久久久久爽电影 | 日日爽夜夜爽网站| 在线看a的网站| 天天躁夜夜躁狠狠躁躁| 在线播放国产精品三级| 99精国产麻豆久久婷婷| 亚洲一区二区三区色噜噜 | 99国产极品粉嫩在线观看| 久久午夜综合久久蜜桃| 午夜亚洲福利在线播放| 国产精品综合久久久久久久免费 | 国产成人av激情在线播放| 国产成人精品无人区| 韩国精品一区二区三区| 久久人妻av系列| 满18在线观看网站| 韩国精品一区二区三区| 国产精品一区二区免费欧美| 亚洲一码二码三码区别大吗| 伊人久久大香线蕉亚洲五|