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

    非均質(zhì)低滲透油藏非線性流動數(shù)值模擬

    2014-06-09 12:33:45黃朝琴
    計算物理 2014年5期
    關鍵詞:壓力梯度算例滲流

    姚 軍, 黃 濤, 黃朝琴

    (中國石油大學(華東)石油工程學院,山東青島 266580)

    非均質(zhì)低滲透油藏非線性流動數(shù)值模擬

    姚 軍, 黃 濤, 黃朝琴

    (中國石油大學(華東)石油工程學院,山東青島 266580)

    提出一種適用于非均質(zhì)低滲透油藏非線性流動數(shù)值模擬的方法.通過引入單元邊界壓力和速度兩組物理量,在單個網(wǎng)格單元中構造數(shù)值計算格式,可精確地獲取每個網(wǎng)格單元的壓力梯度.以單相非線性滲流問題為出發(fā)點,詳細闡述該方法的基本原理;通過數(shù)值算例驗證方法的正確性和適用性.

    非均質(zhì);低滲透油藏;非線性滲流;油藏數(shù)值模擬

    0 引言

    隨著我國油氣田勘探開發(fā)的不斷深入,低滲透油藏已成為我國增儲上產(chǎn)的主要領域[1].大量低滲透滲流實驗和生產(chǎn)實踐表明[2-4]:油氣流動存在最小啟動壓力梯度,其滲流曲線偏離了經(jīng)典的達西定律,具有顯著的非線性特征.許多專家學者對其流動機理進行了大量研究,提出了多種非線性滲流模式,包括擬啟動壓力梯度模型、三階段模型、非線性連續(xù)性模型等[5-7],并進行了相應的數(shù)值模擬研究[8-11].現(xiàn)有的低滲透油藏非線性流動數(shù)值模擬,大都是采用傳統(tǒng)的有限差分法或控制體積有限差分.此類方法均基于相鄰單元中心點的物理量值來構造數(shù)值計算格式,求取驅替壓力梯度,進而結合滲流曲線來進行判斷和求解.然而,對于非均質(zhì)低滲透油藏,不同區(qū)域甚至不同網(wǎng)格單元的滲流模式不盡相同,數(shù)值模擬中必須獲取各網(wǎng)格單元的壓力梯度,并選擇相應的滲流模式,此時上述傳統(tǒng)的數(shù)值計算方法已無法滿足要求.

    對此,在控制體積有限差分法的基礎上,通過引入界面中心點壓力,提出了一種新的數(shù)值算法.該算法基于單個網(wǎng)格單元來構造數(shù)值計算格式,能夠精確地獲取每一網(wǎng)格單元的壓力梯度,單元間則通過流量連續(xù)條件聯(lián)系起來.并以此為基礎,建立了一套非均質(zhì)低滲透油藏非線性流動數(shù)值模擬方法.

    1 滲流數(shù)學模型

    由于受啟動壓力梯度的影響,低滲透油藏流體流動存在明顯的動邊界[12-13].實驗結果表明[14]:低滲透油藏滲流模式一般為如圖1所示的滲流曲線(實線OA-AE),A點表示孔隙介質(zhì)中最大孔道中的流體流動時的啟動壓力梯度,B點表示最小孔道中流體流動時所需的啟動壓力梯度,曲線AD上一點Di的平均啟動壓力梯度為Ci,C為直線DE與橫坐標軸的交點,即擬啟動壓力梯度,D為曲線轉直線的折點.

    之前很多研究采用簡化的擬啟動壓力梯度滲流模式,即用直線OC、CE來描述滲流過程.擬啟動壓力梯度模型在計算時對于在低壓力梯度時阻力較小的大孔道中的流動估計偏低,因而綜合經(jīng)濟技術指標會偏低.考慮某一封閉低滲透油藏區(qū)域Ω中的不可壓縮單相液體的等溫滲流過程,為簡單起見,本文僅研究各向同性問題,并假設巖石微可壓縮,對于擬啟動壓力梯度其數(shù)學模型如下:

    圖1 典型非達西滲流曲線Fig.1 Non-Darcy flow curve

    其中,υ為滲流速度,K為滲透率張量,μ為流體粘度,p為油藏壓力(或流動勢),λ為啟動壓力梯度,φ為巖石孔隙度,q表示源匯項.

    對于實際非線性滲流曲線(圖1實線部分),文獻[15-17]建議采用動態(tài)滲透率的概念,即滲透率和啟動壓力梯度是隨著壓力梯度的增大而增大的,當壓力梯度大于某一臨界值時滲透率和啟動壓力梯度不再發(fā)生變化,如圖2所示.此時,可求出每個驅替壓力梯度對應的滲透率及啟動壓力梯度,具體如下:①由巖心做出低滲介質(zhì)滲流曲線圖;②由驅替壓力梯度,找到對應于滲流曲線上的滲流點;③對該滲流點沿滲流曲線方向做切線,該切線的斜率即為滲透率,而切線與驅替壓力梯度坐標軸的交點即為該壓力梯度時孔道的啟動壓力梯度.顯然,對于動態(tài)滲透率模型,每一個壓力梯度值對應一個擬啟動壓力梯度模型,因此,本文將重點闡述擬啟動壓力梯度滲流數(shù)學模型及其數(shù)值求解方法.對于圖2實線所示的非線性滲流模式,其滲流數(shù)學模型僅需改變方程(1),其中,在曲線AD段滲透率是驅替壓力梯度函數(shù),具體如下

    圖2 滲透介質(zhì)的滲透率隨壓力梯度的變化Fig.2 Mutative permeability vs.pressuregradient

    2 數(shù)值計算格式構造

    假設油藏區(qū)域Ω∈R2被一組互不重疊的多邊形網(wǎng)格Ω={Ωi}所剖分,如圖3所示.

    圖3 網(wǎng)格單元分析示意圖Fig.3 Schematic of grid analysis

    由方程(1)易知邊界面上的法向滲流速度vf與壓力梯度Δp、啟動壓力梯度λ之間有

    然,方程(15)中的第一行對應于方程(1)Darcy定律,第二行對應于方程(2)質(zhì)量守恒定律,第三行則表征了單元邊界面上的法向速度連續(xù)性條件.上述方程的系數(shù)矩陣表達式:

    其中,Ne為網(wǎng)格單元總數(shù);Ii=Em.從方程(16)可知:方程(15)中的各系數(shù)矩陣僅與網(wǎng)格單元的幾何信息和油藏參數(shù)有關,而對網(wǎng)格的幾何形狀沒有特殊要求,求解方便,原則上適用于任何復雜網(wǎng)格.并且方程(14)中啟動壓力項只與每個網(wǎng)格單元有關,對于非均質(zhì)油藏只需改變每個單元啟動壓力梯度,所以具有很好的適應性.

    3 數(shù)值算例

    通過簡單的二維單相流算例對低滲透油藏數(shù)值模擬進行分析.考慮如圖4(a)所示二維1/4五點井網(wǎng)模型,油藏長100 m,寬100 m,在(0,0)點有一注入井,在(100,100)點有一采出井,注入和采出的速度均為17.28 m3·d-1.采用三角網(wǎng)格剖分,節(jié)點數(shù)為190,單元數(shù)為330,如圖4(b)所示.

    圖4 油藏幾何信息圖Fig.4 Reservoir geometry information

    3.1 算例1

    擬啟動壓力梯度模型:設滲透率Kx=Ky=0.1 md(10-3μm2),Kxy=Kyx=0,初始油藏壓力Pe=Pf=10 MPa,擬啟動壓力梯度λ=0.1 MPa·m-1,巖石壓縮系數(shù)Cf=8.64×10-4MPa-1,流體粘度μ=1 mPa·s,t=200 d(天).

    從圖5可以看出,在擬啟動壓力梯度模型滲流過程中,存在明顯的動邊界,當驅替壓力小于啟動壓力時,油藏中的流體是不流動的,油藏壓力保持初始值不變.當壓力傳播到整個油藏時,會形成優(yōu)勢通道,流體流速較快.

    圖5 算例1壓力及速度分布Fig.5 Distributions of pressure and velocity in Example 1

    3.2 算例2

    非線性滲流模型:該模型與上述擬啟動壓力梯度模型基本相同,唯一區(qū)別是,滲透率和啟動壓力梯度隨驅替壓力梯度增大而增大,當達到某一臨界值時保持不變,其對應關系見表1.圖6為非線性滲流模型不同時刻壓力及速度分布圖.

    圖6 算例2壓力及速度分布Fig.6 Distributions of pressure and velocity in Example 2

    如圖7所示,對比兩個模型不同時間步t=80、160、200 d注采井點連線之間壓力分布,可看出擬啟動壓力梯度模型在同一時間步中,壓力傳播范圍較非線性滲流模型小,且靠近源點處相對于初始地層壓力壓差大;而在非線性滲流模型中,壓力傳播較快,油藏較早達到穩(wěn)態(tài).這是因為擬啟動壓力梯度模型是一種簡化的滲流模式,模型假設了一個啟動壓力梯度(圖1中的C點),當驅替壓力梯度小于該啟動壓力梯度時,地層中的流體是不流動的.而實驗結果表明實際滲流模式應為圖1實線所示,當驅替壓力梯度大于某值(圖1中的A點)時,流體已經(jīng)開始流動.所以說,對比2類模型,由于流體流動所需要的啟動壓力大,擬啟動壓力梯度模型在計算時對于在低壓力梯度時阻力較小的大孔道中的流動估計偏低,因而綜合經(jīng)濟技術指標會偏低,非線性滲流模型近似圖1實線部分所示滲流模式,計算時更符合實際.

    圖7 不同時刻兩種模型井間連線壓力值Fig.7 Pressures between two model's wells

    3.3 算例3

    非線性滲流模式算例中均采用了相同的滲透率、啟動壓力梯度隨驅替壓力梯度變化曲線,對于非均質(zhì)油藏,由于計算過程中針對每個網(wǎng)格單元均求出了驅替壓力梯度,所以只要給定區(qū)塊甚至網(wǎng)格單元的變化曲線,都可以求出網(wǎng)格單元相應的滲透率以及啟動壓力梯度,這也是該方法的優(yōu)勢所在,為了簡單起見,均采用同一變化曲線,見表1.

    考慮與算例2中非線性滲流相同的模型,唯一區(qū)別是其滲透率分布不同,在此算例中滲透率Kx、Ky通過MatLab的隨機函數(shù)randn生成,Kxy=Kyx=0,如圖8(a)所示,圖8(b)、8(c)、8(d)為不同時刻壓力及速度分布.

    低滲透油藏非均質(zhì)性極強,網(wǎng)格之間滲透率差別較大,如圖8(a)所示.從數(shù)值模擬結果看出,流體主要沿著滲透率較大的網(wǎng)格流動、流速較大,壓力傳播快,符合實際,并且本文提出的方法適用于非均質(zhì)性油藏.

    圖8 算例3壓力及速度分布Fig.8 Distributions of pressure and velocity in Example 3

    4 結論

    1)提出一種適用于非均質(zhì)低滲透油藏流動數(shù)值模擬方法,以單相滲流問題為出發(fā)點,詳細闡述了方法的基本原理,并通過數(shù)值算例驗證了方法的正確性.

    2)新方法僅基于單個網(wǎng)格單元的信息來構造數(shù)值計算格式,理論上適用于任何復雜網(wǎng)格,可精確地獲取每個網(wǎng)格單元的壓力梯度,適用于非均質(zhì)低滲透油藏數(shù)值模擬,尤其是不同區(qū)域或網(wǎng)格單元具有完全不同滲流模式的情形.

    3)本文僅對非均質(zhì)低滲透油藏的單相滲流問題予以研究,下一步將該方法推廣至兩相、多相滲流以及三維問題.

    [1]Wang Guangfu,Liao Rongfeng,Li Jianglong,et al.The development situation and future of low permeability oil reservoirs of SINOPEC[J].Petroleum Geology and Recovery Efficiency,2007,14(3):84-89.

    [2]Miller R J,Low P F.Threshold gradient for water flow in clay systems[J].Soil Science Society of America Proceedings,963,27(6):605-609.

    [3]Lv Chengyuan,Wang Jian,Sun Zhigang,et al.An experimental study on starting pressure gradient of fluids flow in lowpermeability sandstone porous media[J].Petroleum Exploration and Development,2002,29(2):86-89.

    [4]Wan Quanlin,Tan Hai,Lv Dongliang,et al.An experimental study on threshold pressure gradient in low permeability reservoir[J].Petroleum Geology and Recovery Efficiency,2011,18(1):97-100.

    [5]Huang Yanzhang.Nonlinear percolation feature in low permeability reservoir[J].Special Oil&Gas Reservoirs,1997,4(1):9-14.

    [6]Yao Yuedong,Ge Jiali.Study on non Darcy flow pattern in low permeability oil reservoir[J].Xinjiang Petroleum Geology,2000,21(3):213-215.

    [7]Xiao Zengli,Pu Chunsheng,Qin Wenlong,et al.Nonlinear percolation feature of low permeability reservoirs and its effect [J].Oil Drilling&Production Technology,2007,29(3):105-107.

    [8]Guo Yongcun,Lu Detan,Ma Linxiao,et al.Numerical simulation of fluid flow in low permeability reservoir using finite difference method[J].Journal of Hydrodynamics,2004,19(3):288-293.

    [9]Han Hongbao,Cheng Linsong,Zhang Minglu,et al.Physical simulation and numerical simulation of ultra-low permeability reservoir in consideration of starting pressure gradient[J].Journal of University of Petroleum,China(Natural Science Edition),2004,28(6):49-53.

    [10]Xue Guoqing,Li Ming,Wang Daocheng,et al.Numerical simulation of low permeability and heterogeneous reservoir[J]. Fault-Block Oil and Gas Field,2007,14(4):43-44.

    [11]Yang Zhengmin,Yu Rongcheng,Su Zhixin,et al.Numerical simulation of the nonlinear flow in ultra-low permeability reservoirs[J].Petroleum Exploration and Development,2010,7(1):94-98.

    [12]Guo Yongcun,Zeng Qinghong,Wang Zhongxun.Numerical simulation of low permeability flow with moving-boundary using meshless methods[J].Engineering Mechanics,2006,23(11):188-192.

    [13]Liu Wenchao,Yao Jun,Wang Jianzhong.Front tracking for moving boundary of non-Darcy seepage flows in low-permeability porous media[J].Chinese Journal of Computational Physics,2012,29(6):823-827.

    [14]Yan Qinglai,He Qiuxuan,Yu Ligang.A laboratory study on percolation characteristics of single phase flow in low-permeability reserviors[J].Journal of Xi'an Petroleum Institute,1990,5(2):1-6.

    [15]Lei Qun,Xiong Wei,Yuan Jiangru,et al.Behavior of flow through low-permeability reservoirs[A].SPE 113144,2008.

    [16]Yao Jun,Liu Shun.Well test interpretation model based on mutative permeability effects for low-permeability reservoir[J]. Acta Petrolei Sinica,2009,30(3):430-433.

    [17]Yin Zhilin,Sun Wenjing,Yao Jun.Numerical simulation of the 3D oil-water phase dynamic permeability for low-permeability reservoirs[J].Acta Petrolei Sinica,2011,32(1):117-121.

    Numerical Simulation of Nonlinear Flow in Heterogeneous and Low-permeability Reservoirs

    YAO Jun,HUANG Tao,HUANG Zhaoqin
    (School of Petroleum Engineering,China University of Petroleum,Qingdao 266580,China)

    A numerical method is developed to model nonlinear flow in heterogeneous and low-permeability reservoirs.In order to obtain pressure gradient of each gridcells precisely,we introduced two sets of physical quantities,cell face pressures and surface efflux for each gridcell.Principle of the method is elaborated based on single-phase fluid flow problem.Several numerical examples are presented to demonstrate efficiency and applicability of the scheme.

    heterogeneity;low-permeability reservoirs;nonlinear porous flow;numerical simulation

    date:2013-10-15;Revised date:2014-02-19

    TE348;TE319;O241.82

    A

    2013-10-15;

    2014-02-19

    國家重點基礎研究發(fā)展計劃(2011CB201004),國家自然科學基金(51234007),國家重大科技專項(2011ZX05014-005-003HZ)和中央高?;究蒲袠I(yè)務費專項資金(13CX05007A,14CX02042A,13CX05017A,14CX06091A)資助項目

    姚軍(1964-),男,山東臨沂,博士生導師,主要從事油氣田開發(fā)工程教研工作,E-mail:RCOGFR_UPC@126.com

    1001-246X(2014)05-0552-07

    猜你喜歡
    壓力梯度算例滲流
    壓力梯度在油田開發(fā)中的應用探討
    基于振蕩能量的低頻振蕩分析與振蕩源定位(二)振蕩源定位方法與算例
    互補問題算例分析
    疊加原理不能求解含啟動壓力梯度滲流方程
    基于CYMDIST的配電網(wǎng)運行優(yōu)化技術及算例分析
    燃煤PM10湍流聚并GDE方程算法及算例分析
    致密砂巖啟動壓力梯度數(shù)值的影響因素
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關于渠道滲流計算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    尾礦壩滲流計算及排滲設計
    金屬礦山(2013年6期)2013-03-11 16:54:05
    人人妻人人澡人人看| 99国产精品一区二区蜜桃av| 亚洲男人天堂网一区| 国产又色又爽无遮挡免费看| 免费看a级黄色片| 999久久久国产精品视频| 精品国产超薄肉色丝袜足j| www.www免费av| 国产黄色小视频在线观看| 亚洲av电影在线进入| 欧美性猛交黑人性爽| 久久久久久国产a免费观看| 亚洲国产欧洲综合997久久, | 嫁个100分男人电影在线观看| 亚洲五月天丁香| 人成视频在线观看免费观看| 久久精品国产亚洲av高清一级| 国产亚洲精品久久久久5区| 亚洲av电影在线进入| 亚洲av电影在线进入| 90打野战视频偷拍视频| 亚洲中文av在线| 久久久久久久午夜电影| 国产一区二区激情短视频| 老汉色av国产亚洲站长工具| 亚洲av五月六月丁香网| 欧美国产精品va在线观看不卡| 欧美精品啪啪一区二区三区| 婷婷亚洲欧美| 国产成人精品久久二区二区免费| 亚洲av片天天在线观看| 国产伦一二天堂av在线观看| 搞女人的毛片| 成年版毛片免费区| 婷婷亚洲欧美| 国产亚洲精品久久久久久毛片| 女同久久另类99精品国产91| 男人舔女人的私密视频| 午夜a级毛片| 久久精品影院6| 欧美黑人欧美精品刺激| 啦啦啦 在线观看视频| 国产伦人伦偷精品视频| 一本久久中文字幕| 桃色一区二区三区在线观看| 国产片内射在线| 中文字幕av电影在线播放| 99国产精品一区二区三区| 日本成人三级电影网站| 亚洲av成人不卡在线观看播放网| 中文字幕人成人乱码亚洲影| 亚洲av成人av| 性色av乱码一区二区三区2| 嫁个100分男人电影在线观看| 97超级碰碰碰精品色视频在线观看| 首页视频小说图片口味搜索| 亚洲欧美激情综合另类| 中文字幕最新亚洲高清| 国产亚洲精品av在线| 一级黄色大片毛片| 日本免费一区二区三区高清不卡| 国产真人三级小视频在线观看| 啦啦啦 在线观看视频| 99久久无色码亚洲精品果冻| av福利片在线| 99热这里只有精品一区 | 99热只有精品国产| 久久中文字幕一级| 成年人黄色毛片网站| 国产真实乱freesex| 欧美日韩亚洲国产一区二区在线观看| 窝窝影院91人妻| 国产精品国产高清国产av| 日韩欧美 国产精品| 日韩精品免费视频一区二区三区| 日韩三级视频一区二区三区| 99久久无色码亚洲精品果冻| 久久这里只有精品19| 亚洲欧洲精品一区二区精品久久久| 久久久久久久久中文| 黑人巨大精品欧美一区二区mp4| 欧美又色又爽又黄视频| 久久天堂一区二区三区四区| 欧美日韩乱码在线| 亚洲色图 男人天堂 中文字幕| 精品午夜福利视频在线观看一区| 久久久久国产一级毛片高清牌| 夜夜躁狠狠躁天天躁| 天天一区二区日本电影三级| 每晚都被弄得嗷嗷叫到高潮| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲国产精品sss在线观看| 国产精品一区二区三区四区久久 | 一本精品99久久精品77| 香蕉久久夜色| 国产一区二区三区视频了| 欧美在线黄色| or卡值多少钱| 青草久久国产| 国产久久久一区二区三区| 国产成人av激情在线播放| www日本在线高清视频| 999久久久国产精品视频| 国产欧美日韩一区二区三| 香蕉丝袜av| 麻豆成人午夜福利视频| √禁漫天堂资源中文www| 亚洲av片天天在线观看| 色哟哟哟哟哟哟| 亚洲精品久久成人aⅴ小说| 黑丝袜美女国产一区| 天天躁狠狠躁夜夜躁狠狠躁| 日韩欧美 国产精品| 免费av毛片视频| 久久久精品欧美日韩精品| 禁无遮挡网站| 日本在线视频免费播放| 色播在线永久视频| 国产99白浆流出| 精品国产一区二区三区四区第35| 色综合站精品国产| 在线观看免费日韩欧美大片| 午夜福利成人在线免费观看| 成人欧美大片| 中文亚洲av片在线观看爽| www日本黄色视频网| 成在线人永久免费视频| 18禁裸乳无遮挡免费网站照片 | 一区二区三区国产精品乱码| 日本黄色视频三级网站网址| 男人舔女人下体高潮全视频| 露出奶头的视频| 后天国语完整版免费观看| 级片在线观看| 婷婷亚洲欧美| 成人三级黄色视频| 欧美日韩中文字幕国产精品一区二区三区| 丝袜人妻中文字幕| 欧美成人免费av一区二区三区| 色综合婷婷激情| 丝袜在线中文字幕| 一进一出好大好爽视频| 欧美在线一区亚洲| 亚洲免费av在线视频| 久久人妻av系列| 亚洲av五月六月丁香网| 精品久久久久久,| 欧美日韩中文字幕国产精品一区二区三区| 久久性视频一级片| √禁漫天堂资源中文www| 日本熟妇午夜| av中文乱码字幕在线| videosex国产| 国产欧美日韩一区二区精品| 91大片在线观看| 麻豆一二三区av精品| 观看免费一级毛片| 免费观看精品视频网站| 亚洲中文日韩欧美视频| 动漫黄色视频在线观看| 一区二区日韩欧美中文字幕| av在线天堂中文字幕| 国产极品粉嫩免费观看在线| 日本一本二区三区精品| 女同久久另类99精品国产91| 久久精品国产亚洲av香蕉五月| 亚洲欧美一区二区三区黑人| av片东京热男人的天堂| 久久国产亚洲av麻豆专区| 麻豆久久精品国产亚洲av| 1024香蕉在线观看| 欧美黑人巨大hd| 99热这里只有精品一区 | 久久精品国产亚洲av高清一级| 免费在线观看黄色视频的| 黑人欧美特级aaaaaa片| 一边摸一边做爽爽视频免费| 久久精品国产亚洲av香蕉五月| 99国产精品一区二区蜜桃av| 99久久综合精品五月天人人| 美国免费a级毛片| 国产精品久久久久久亚洲av鲁大| АⅤ资源中文在线天堂| 青草久久国产| 亚洲在线自拍视频| 亚洲电影在线观看av| 国产精品久久久久久精品电影 | 国产黄片美女视频| 啦啦啦免费观看视频1| 俄罗斯特黄特色一大片| 久久这里只有精品19| 国内久久婷婷六月综合欲色啪| 欧美乱色亚洲激情| 人人妻人人澡欧美一区二区| av在线天堂中文字幕| 夜夜夜夜夜久久久久| 亚洲一区二区三区不卡视频| 亚洲精品在线观看二区| 成人亚洲精品一区在线观看| 国产极品粉嫩免费观看在线| 一级毛片高清免费大全| 精品电影一区二区在线| 国产亚洲欧美在线一区二区| 日本一区二区免费在线视频| 亚洲av中文字字幕乱码综合 | 国产成人一区二区三区免费视频网站| 看片在线看免费视频| 999久久久国产精品视频| 国产av不卡久久| 俺也久久电影网| 91大片在线观看| 动漫黄色视频在线观看| 中文字幕人妻丝袜一区二区| 亚洲成av人片免费观看| 精品久久久久久久久久免费视频| 精品国产亚洲在线| 十八禁网站免费在线| 99精品久久久久人妻精品| 啦啦啦观看免费观看视频高清| 亚洲精品中文字幕在线视频| 亚洲国产欧洲综合997久久, | 亚洲专区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 热99re8久久精品国产| 热re99久久国产66热| 亚洲真实伦在线观看| 成在线人永久免费视频| 桃色一区二区三区在线观看| 黑人操中国人逼视频| av中文乱码字幕在线| 制服丝袜大香蕉在线| 狂野欧美激情性xxxx| 天天一区二区日本电影三级| av在线播放免费不卡| 国产精品 欧美亚洲| 国产精品1区2区在线观看.| 91麻豆精品激情在线观看国产| 十八禁人妻一区二区| 两个人看的免费小视频| 丝袜在线中文字幕| 久久久久久久午夜电影| 国产97色在线日韩免费| 日韩国内少妇激情av| 国产真人三级小视频在线观看| 啪啪无遮挡十八禁网站| 高清在线国产一区| 精品一区二区三区视频在线观看免费| 亚洲成av片中文字幕在线观看| 国产欧美日韩一区二区精品| 男男h啪啪无遮挡| 91成人精品电影| 18禁美女被吸乳视频| 日本撒尿小便嘘嘘汇集6| 一本综合久久免费| 亚洲第一av免费看| 久久亚洲真实| 欧美大码av| 国产高清视频在线播放一区| 欧美性猛交╳xxx乱大交人| 免费看日本二区| 亚洲经典国产精华液单| 99热只有精品国产| 最近2019中文字幕mv第一页| 亚洲欧美成人精品一区二区| 亚洲成人中文字幕在线播放| 男女视频在线观看网站免费| av在线老鸭窝| 国产高清视频在线播放一区| 乱人视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人av在线免费| 久久人人爽人人爽人人片va| 熟女人妻精品中文字幕| 亚洲欧美成人精品一区二区| 精品久久久久久久久亚洲| 欧美成人a在线观看| 亚洲第一电影网av| 亚洲成人久久性| 中文字幕熟女人妻在线| 91精品国产九色| 婷婷精品国产亚洲av| 亚洲精品色激情综合| 国产麻豆成人av免费视频| 搞女人的毛片| 久久精品国产清高在天天线| 国产精品嫩草影院av在线观看| 免费av观看视频| 日韩中字成人| 亚洲精品国产av成人精品 | 九九热线精品视视频播放| 欧美3d第一页| 婷婷精品国产亚洲av| 久久中文看片网| 午夜爱爱视频在线播放| 熟妇人妻久久中文字幕3abv| 亚洲综合色惰| 精品人妻熟女av久视频| 日本熟妇午夜| 深夜精品福利| 少妇熟女欧美另类| 国产毛片a区久久久久| 国产精品,欧美在线| 简卡轻食公司| 欧美+日韩+精品| 男人舔女人下体高潮全视频| 午夜精品在线福利| 亚洲欧美成人精品一区二区| 久久这里只有精品中国| 成人av在线播放网站| 国产精品99久久久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 国产午夜福利久久久久久| 97在线视频观看| 精品日产1卡2卡| 日本黄大片高清| 免费大片18禁| 老师上课跳d突然被开到最大视频| 老熟妇仑乱视频hdxx| 国产高清有码在线观看视频| 精品一区二区三区av网在线观看| 人人妻,人人澡人人爽秒播| 免费人成视频x8x8入口观看| 国产色婷婷99| 亚洲精品日韩av片在线观看| 国产精华一区二区三区| 亚洲无线观看免费| 九九久久精品国产亚洲av麻豆| 亚洲国产色片| 日韩,欧美,国产一区二区三区 | 在线a可以看的网站| 精品午夜福利在线看| 色视频www国产| 免费av不卡在线播放| 久久久久久国产a免费观看| 精品99又大又爽又粗少妇毛片| 日韩欧美 国产精品| 女人十人毛片免费观看3o分钟| 亚洲人成网站高清观看| 亚洲自拍偷在线| 久久九九热精品免费| 成人永久免费在线观看视频| 亚洲av.av天堂| 秋霞在线观看毛片| 欧美区成人在线视频| 听说在线观看完整版免费高清| 国内久久婷婷六月综合欲色啪| 欧美+亚洲+日韩+国产| 夜夜爽天天搞| 日韩一区二区视频免费看| 一区福利在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产欧美人成| 女的被弄到高潮叫床怎么办| 免费在线观看影片大全网站| 校园人妻丝袜中文字幕| 日韩大尺度精品在线看网址| 18禁黄网站禁片免费观看直播| 一本久久中文字幕| 久久精品国产亚洲av涩爱 | 我要搜黄色片| 成人美女网站在线观看视频| www日本黄色视频网| 亚洲欧美中文字幕日韩二区| 国产毛片a区久久久久| 久久久久久国产a免费观看| 国内久久婷婷六月综合欲色啪| 欧美又色又爽又黄视频| 美女cb高潮喷水在线观看| 精品少妇黑人巨大在线播放 | 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 亚洲人成网站在线播| 寂寞人妻少妇视频99o| 精品免费久久久久久久清纯| 少妇的逼好多水| 亚洲久久久久久中文字幕| 又粗又爽又猛毛片免费看| 美女内射精品一级片tv| 国产亚洲精品av在线| 国产精品人妻久久久影院| 色av中文字幕| 午夜福利在线在线| 亚洲欧美清纯卡通| 两性午夜刺激爽爽歪歪视频在线观看| 九色成人免费人妻av| 99久久成人亚洲精品观看| av中文乱码字幕在线| 最近在线观看免费完整版| 欧美不卡视频在线免费观看| 乱人视频在线观看| 精品久久久噜噜| 亚洲美女视频黄频| 桃色一区二区三区在线观看| 1000部很黄的大片| 久久午夜福利片| 伦精品一区二区三区| 十八禁网站免费在线| 亚洲久久久久久中文字幕| 大型黄色视频在线免费观看| 亚洲精品久久国产高清桃花| 91狼人影院| 欧美精品国产亚洲| 国内久久婷婷六月综合欲色啪| 欧美日韩在线观看h| 可以在线观看的亚洲视频| 成人亚洲欧美一区二区av| 免费观看人在逋| 成人一区二区视频在线观看| 久久婷婷人人爽人人干人人爱| 麻豆av噜噜一区二区三区| 国产高清视频在线播放一区| 性色avwww在线观看| 欧美激情久久久久久爽电影| 露出奶头的视频| 亚洲三级黄色毛片| 九九久久精品国产亚洲av麻豆| 性欧美人与动物交配| 午夜福利高清视频| 日本一二三区视频观看| 三级经典国产精品| 校园人妻丝袜中文字幕| 18+在线观看网站| av免费在线看不卡| 午夜福利成人在线免费观看| 免费看a级黄色片| 国产成人aa在线观看| 国产极品精品免费视频能看的| 色播亚洲综合网| 国产精品99久久久久久久久| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 高清日韩中文字幕在线| 特大巨黑吊av在线直播| 久久人人精品亚洲av| 美女 人体艺术 gogo| 两个人的视频大全免费| 午夜老司机福利剧场| 欧美激情久久久久久爽电影| 免费看美女性在线毛片视频| 少妇的逼水好多| 99在线人妻在线中文字幕| 久久精品国产鲁丝片午夜精品| 精品久久久久久久久亚洲| 国产一区二区亚洲精品在线观看| 乱码一卡2卡4卡精品| 男女做爰动态图高潮gif福利片| 中出人妻视频一区二区| 免费看a级黄色片| 久99久视频精品免费| 一级黄色大片毛片| 99久久无色码亚洲精品果冻| 国产精品乱码一区二三区的特点| 亚洲性夜色夜夜综合| 成人精品一区二区免费| 免费高清视频大片| 真实男女啪啪啪动态图| 丰满的人妻完整版| 欧美高清性xxxxhd video| 欧美日本亚洲视频在线播放| 最近2019中文字幕mv第一页| 12—13女人毛片做爰片一| 亚洲四区av| 欧美三级亚洲精品| 亚洲精品久久国产高清桃花| 亚洲经典国产精华液单| 色播亚洲综合网| 内射极品少妇av片p| 国产欧美日韩精品亚洲av| 久久久久久久亚洲中文字幕| 国产私拍福利视频在线观看| 欧美高清性xxxxhd video| 久久精品国产自在天天线| 国产v大片淫在线免费观看| 欧美日韩在线观看h| 欧美色欧美亚洲另类二区| 亚洲精品在线观看二区| 中文字幕精品亚洲无线码一区| 成人毛片a级毛片在线播放| 老司机影院成人| 亚洲一区二区三区色噜噜| 一本一本综合久久| 寂寞人妻少妇视频99o| 亚洲四区av| 寂寞人妻少妇视频99o| 最近手机中文字幕大全| 三级男女做爰猛烈吃奶摸视频| 亚洲五月天丁香| 免费av不卡在线播放| 久久天躁狠狠躁夜夜2o2o| 白带黄色成豆腐渣| 女人十人毛片免费观看3o分钟| 伊人久久精品亚洲午夜| 综合色丁香网| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| 免费看光身美女| av视频在线观看入口| 日韩精品中文字幕看吧| 国产高清有码在线观看视频| 久久久久国产精品人妻aⅴ院| 日日干狠狠操夜夜爽| 日韩国内少妇激情av| а√天堂www在线а√下载| 人妻制服诱惑在线中文字幕| 两个人视频免费观看高清| 欧美激情久久久久久爽电影| 日韩av在线大香蕉| 亚洲色图av天堂| 国产精品久久视频播放| 女同久久另类99精品国产91| 久久久久九九精品影院| 国产熟女欧美一区二区| 国产成人影院久久av| 一级黄色大片毛片| 国模一区二区三区四区视频| 全区人妻精品视频| 99久国产av精品国产电影| 精品一区二区三区人妻视频| 成人性生交大片免费视频hd| 欧美+亚洲+日韩+国产| 国产精品一二三区在线看| 自拍偷自拍亚洲精品老妇| 成人二区视频| 男人舔女人下体高潮全视频| 亚洲图色成人| 国产高清激情床上av| 午夜久久久久精精品| 日本精品一区二区三区蜜桃| 精品国内亚洲2022精品成人| 97在线视频观看| 日日摸夜夜添夜夜添av毛片| 免费大片18禁| 婷婷亚洲欧美| 性色avwww在线观看| 寂寞人妻少妇视频99o| 乱系列少妇在线播放| 波多野结衣高清无吗| 麻豆精品久久久久久蜜桃| 亚洲精品日韩av片在线观看| 亚洲人成网站在线播| 91麻豆精品激情在线观看国产| 最近最新中文字幕大全电影3| 一进一出好大好爽视频| 亚洲第一电影网av| 色吧在线观看| 欧美潮喷喷水| 乱人视频在线观看| 亚洲无线观看免费| 日韩av在线大香蕉| 亚洲中文字幕日韩| 校园春色视频在线观看| 久久婷婷人人爽人人干人人爱| 国产色爽女视频免费观看| 国产成人福利小说| 在线播放无遮挡| 亚洲欧美日韩卡通动漫| 在线a可以看的网站| 欧美色视频一区免费| 日韩精品青青久久久久久| 麻豆av噜噜一区二区三区| 一级a爱片免费观看的视频| 亚洲久久久久久中文字幕| 免费电影在线观看免费观看| 国产色婷婷99| a级毛片免费高清观看在线播放| 久久久久久久久久成人| 日日干狠狠操夜夜爽| 熟女电影av网| 亚洲最大成人手机在线| 午夜福利高清视频| 欧美色欧美亚洲另类二区| 国产真实伦视频高清在线观看| 高清毛片免费看| 亚洲18禁久久av| 国产亚洲精品av在线| 91久久精品国产一区二区成人| 国产亚洲精品久久久久久毛片| 午夜爱爱视频在线播放| 亚洲自拍偷在线| 日本在线视频免费播放| 亚洲一级一片aⅴ在线观看| 亚洲在线观看片| 舔av片在线| 欧美激情久久久久久爽电影| 午夜免费激情av| 男女啪啪激烈高潮av片| 欧美3d第一页| 国产成人91sexporn| 99精品在免费线老司机午夜| 免费看美女性在线毛片视频| 欧美最新免费一区二区三区| 午夜精品在线福利| 日产精品乱码卡一卡2卡三| 99久国产av精品| 美女黄网站色视频| 亚洲av中文字字幕乱码综合| 最近中文字幕高清免费大全6| 99视频精品全部免费 在线| 亚洲国产精品成人久久小说 | 日本与韩国留学比较| 日本黄色片子视频| 国产成人91sexporn| 欧美最黄视频在线播放免费| 综合色丁香网| h日本视频在线播放| 美女cb高潮喷水在线观看| 午夜福利在线在线| 亚洲精品乱码久久久v下载方式| 天堂√8在线中文| 欧美日本视频| 91久久精品国产一区二区成人| 高清午夜精品一区二区三区 | 国产91av在线免费观看| 日韩三级伦理在线观看|