• <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
    久久九九热精品免费| 三上悠亚av全集在线观看| 日日摸夜夜添夜夜添小说| 国产激情久久老熟女| 大码成人一级视频| 成人国语在线视频| 国产精品美女特级片免费视频播放器 | 一区二区三区精品91| 午夜福利,免费看| 无遮挡黄片免费观看| 午夜日韩欧美国产| 精品少妇久久久久久888优播| 欧美日韩亚洲国产一区二区在线观看 | 午夜精品国产一区二区电影| 国产无遮挡羞羞视频在线观看| 97人妻天天添夜夜摸| 亚洲专区国产一区二区| 另类精品久久| 在线亚洲精品国产二区图片欧美| 日韩一区二区三区影片| 精品一区二区三区视频在线观看免费 | 热re99久久精品国产66热6| 亚洲成人免费电影在线观看| 午夜激情av网站| 成人影院久久| 在线播放国产精品三级| 久久久久视频综合| 18在线观看网站| 丝袜人妻中文字幕| 制服诱惑二区| 欧美成狂野欧美在线观看| 国产av又大| 操出白浆在线播放| 国产成人av教育| 精品一区二区三区四区五区乱码| 日韩一卡2卡3卡4卡2021年| 国产成人欧美在线观看 | 99久久国产精品久久久| 国产成人免费观看mmmm| 亚洲国产av新网站| 老司机午夜十八禁免费视频| 国产色视频综合| 免费在线观看黄色视频的| 一进一出抽搐动态| 精品国产一区二区久久| 视频区图区小说| 两人在一起打扑克的视频| 欧美黄色淫秽网站| 首页视频小说图片口味搜索| 黄色丝袜av网址大全| 黑人猛操日本美女一级片| 香蕉国产在线看| 99精品久久久久人妻精品| 欧美黄色淫秽网站| svipshipincom国产片| 免费av中文字幕在线| 亚洲情色 制服丝袜| av片东京热男人的天堂| 狂野欧美激情性xxxx| 午夜激情久久久久久久| 亚洲精品成人av观看孕妇| 久久热在线av| 高潮久久久久久久久久久不卡| 一级片'在线观看视频| 欧美久久黑人一区二区| 黄色a级毛片大全视频| 国产1区2区3区精品| 99国产精品99久久久久| 国产在视频线精品| 亚洲成人免费电影在线观看| 欧美日韩亚洲高清精品| 纯流量卡能插随身wifi吗| 久久久久久人人人人人| 999精品在线视频| 女人被躁到高潮嗷嗷叫费观| www.自偷自拍.com| 国产野战对白在线观看| 欧美 日韩 精品 国产| 韩国精品一区二区三区| 97在线人人人人妻| 午夜福利欧美成人| 欧美变态另类bdsm刘玥| 男女床上黄色一级片免费看| 五月开心婷婷网| 91成年电影在线观看| 免费高清在线观看日韩| 久久人妻熟女aⅴ| 亚洲黑人精品在线| 日本av免费视频播放| 热99re8久久精品国产| 国产一区二区在线观看av| 18禁裸乳无遮挡动漫免费视频| 国产精品免费大片| 精品第一国产精品| 亚洲中文日韩欧美视频| 欧美日韩国产mv在线观看视频| 精品久久久精品久久久| 在线 av 中文字幕| 国产精品亚洲一级av第二区| 午夜两性在线视频| 激情视频va一区二区三区| 999久久久国产精品视频| 久久精品国产a三级三级三级| 国产野战对白在线观看| 午夜福利欧美成人| 亚洲久久久国产精品| 日韩人妻精品一区2区三区| 操出白浆在线播放| 青青草视频在线视频观看| 一边摸一边抽搐一进一出视频| 9热在线视频观看99| 不卡一级毛片| 极品少妇高潮喷水抽搐| 国产淫语在线视频| 男女床上黄色一级片免费看| 国产有黄有色有爽视频| 超色免费av| 99热网站在线观看| 成人18禁在线播放| 黄色怎么调成土黄色| 亚洲自偷自拍图片 自拍| 国产精品1区2区在线观看. | 18禁裸乳无遮挡动漫免费视频| 国产亚洲午夜精品一区二区久久| av国产精品久久久久影院| 国产人伦9x9x在线观看| 91大片在线观看| 精品少妇内射三级| 色综合婷婷激情| 一本综合久久免费| av免费在线观看网站| 精品一区二区三区视频在线观看免费 | 欧美精品亚洲一区二区| 在线观看66精品国产| www.自偷自拍.com| 亚洲成人免费av在线播放| 人妻 亚洲 视频| 男女床上黄色一级片免费看| 91av网站免费观看| 高清毛片免费观看视频网站 | 无限看片的www在线观看| 九色亚洲精品在线播放| 又紧又爽又黄一区二区| 精品亚洲乱码少妇综合久久| a级片在线免费高清观看视频| 亚洲av日韩精品久久久久久密| 中国美女看黄片| 国产成人免费观看mmmm| 色视频在线一区二区三区| 中文字幕精品免费在线观看视频| 欧美精品亚洲一区二区| 亚洲少妇的诱惑av| 亚洲成人免费av在线播放| 人人妻人人爽人人添夜夜欢视频| 国产伦人伦偷精品视频| 国产黄频视频在线观看| 久久这里只有精品19| www.熟女人妻精品国产| 黑人操中国人逼视频| 九色亚洲精品在线播放| 亚洲成人免费av在线播放| 久久久久久免费高清国产稀缺| 韩国精品一区二区三区| 国产免费av片在线观看野外av| 国产三级黄色录像| 国产成人精品无人区| 在线观看免费日韩欧美大片| 国产成人免费观看mmmm| 中文字幕最新亚洲高清| 日日夜夜操网爽| 精品福利观看| 动漫黄色视频在线观看| av欧美777| 亚洲性夜色夜夜综合| 亚洲国产欧美一区二区综合| 亚洲人成伊人成综合网2020| 伦理电影免费视频| 色婷婷久久久亚洲欧美| 亚洲熟女毛片儿| 黑人巨大精品欧美一区二区蜜桃| 午夜福利一区二区在线看| 久久99一区二区三区| 三级毛片av免费| 久久午夜综合久久蜜桃| 曰老女人黄片| 无遮挡黄片免费观看| 90打野战视频偷拍视频| 成年人黄色毛片网站| 高清在线国产一区| 欧美人与性动交α欧美软件| 狠狠精品人妻久久久久久综合| 在线观看免费视频网站a站| 成人亚洲精品一区在线观看| 欧美在线一区亚洲| 国产精品秋霞免费鲁丝片| 国产亚洲午夜精品一区二区久久| 免费在线观看视频国产中文字幕亚洲| 纯流量卡能插随身wifi吗| 深夜精品福利| 夜夜爽天天搞| 久久性视频一级片| 一级毛片电影观看| 久久久久网色| 午夜免费成人在线视频| 欧美变态另类bdsm刘玥| 精品免费久久久久久久清纯 | 日本av免费视频播放| 在线播放国产精品三级| 又黄又粗又硬又大视频| 不卡一级毛片| avwww免费| 老司机影院毛片| 女警被强在线播放| 在线播放国产精品三级| 色老头精品视频在线观看| 色在线成人网| 叶爱在线成人免费视频播放| 最黄视频免费看| 9色porny在线观看| 欧美人与性动交α欧美软件| 日本五十路高清| 国产一区二区激情短视频| 久久精品亚洲av国产电影网| 熟女少妇亚洲综合色aaa.| 午夜日韩欧美国产| 老汉色∧v一级毛片| 99re6热这里在线精品视频| 这个男人来自地球电影免费观看| 狂野欧美激情性xxxx| 国产1区2区3区精品| 成人av一区二区三区在线看| 国产主播在线观看一区二区| 久久中文字幕一级| 精品一区二区三区四区五区乱码| 亚洲国产中文字幕在线视频| 成人永久免费在线观看视频 | 日韩人妻精品一区2区三区| 亚洲av第一区精品v没综合| 久久久久久人人人人人| 久久热在线av| 国产真人三级小视频在线观看| e午夜精品久久久久久久| 免费高清在线观看日韩| 欧美精品一区二区大全| 久久香蕉激情| svipshipincom国产片| 在线 av 中文字幕| 欧美+亚洲+日韩+国产| 热99久久久久精品小说推荐| 免费观看人在逋| 日韩有码中文字幕| 色综合欧美亚洲国产小说| 女人爽到高潮嗷嗷叫在线视频| 老汉色∧v一级毛片| 美女高潮喷水抽搐中文字幕| 满18在线观看网站| 国产成人一区二区三区免费视频网站| 欧美国产精品一级二级三级| 9热在线视频观看99| 亚洲美女黄片视频| 国产三级黄色录像| 99riav亚洲国产免费| 2018国产大陆天天弄谢| 又黄又粗又硬又大视频| 国产男靠女视频免费网站| 老汉色av国产亚洲站长工具| 久久99一区二区三区| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| av天堂在线播放| av有码第一页| 丝袜人妻中文字幕| 亚洲国产中文字幕在线视频| 欧美精品啪啪一区二区三区| 国产日韩欧美在线精品| 麻豆av在线久日| 国产免费现黄频在线看| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 亚洲精品一卡2卡三卡4卡5卡| 国产真人三级小视频在线观看| 性色av乱码一区二区三区2| 精品国产一区二区久久| 我的亚洲天堂| 精品一区二区三区视频在线观看免费 | 乱人伦中国视频| 精品少妇黑人巨大在线播放| 1024视频免费在线观看| 日韩一卡2卡3卡4卡2021年| 免费在线观看日本一区| 侵犯人妻中文字幕一二三四区| 一区福利在线观看| 中文字幕人妻丝袜制服| 日本撒尿小便嘘嘘汇集6| xxxhd国产人妻xxx| 成人av一区二区三区在线看| 在线天堂中文资源库| 一区二区三区精品91| 欧美日韩亚洲国产一区二区在线观看 | 国产不卡一卡二| 欧美日韩亚洲综合一区二区三区_| 久久九九热精品免费| 日日爽夜夜爽网站| 99国产极品粉嫩在线观看| 天堂中文最新版在线下载| 黄色丝袜av网址大全| 亚洲国产欧美日韩在线播放| 18禁美女被吸乳视频| 久久久精品94久久精品| 男人舔女人的私密视频| 免费女性裸体啪啪无遮挡网站| 欧美大码av| 国产深夜福利视频在线观看| 69av精品久久久久久 | 蜜桃在线观看..| 久久这里只有精品19| 在线观看66精品国产| 国产精品自产拍在线观看55亚洲 | 免费一级毛片在线播放高清视频 | 亚洲av日韩在线播放| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 99精品欧美一区二区三区四区| 亚洲国产精品一区二区三区在线| 一本色道久久久久久精品综合| 757午夜福利合集在线观看| 一区二区日韩欧美中文字幕| a级毛片在线看网站| 久久影院123| 精品国产亚洲在线| 亚洲成av片中文字幕在线观看| 99香蕉大伊视频| av国产精品久久久久影院| 午夜精品国产一区二区电影| 国产精品九九99| 亚洲天堂av无毛| 国产精品一区二区免费欧美| 久久人人97超碰香蕉20202| 丝袜美腿诱惑在线| 欧美一级毛片孕妇| 亚洲三区欧美一区| 99国产精品免费福利视频| 免费人妻精品一区二区三区视频| 精品国内亚洲2022精品成人 | 日韩中文字幕视频在线看片| 国产精品国产高清国产av | 国产91精品成人一区二区三区 | 曰老女人黄片| 91精品三级在线观看| bbb黄色大片| 一本大道久久a久久精品| 久久精品熟女亚洲av麻豆精品| 精品福利永久在线观看| 国产高清videossex| 每晚都被弄得嗷嗷叫到高潮| 在线天堂中文资源库| 人妻久久中文字幕网| 在线观看66精品国产| 国产午夜精品久久久久久| 夜夜夜夜夜久久久久| 岛国在线观看网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精华国产精华精| 国产精品免费大片| 亚洲一区中文字幕在线| 欧美激情 高清一区二区三区| 波多野结衣一区麻豆| 亚洲中文日韩欧美视频| 精品人妻1区二区| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 狠狠精品人妻久久久久久综合| 国产精品av久久久久免费| 色综合欧美亚洲国产小说| 亚洲精品中文字幕在线视频| 免费高清在线观看日韩| 狂野欧美激情性xxxx| 日韩中文字幕欧美一区二区| 一二三四社区在线视频社区8| 99热网站在线观看| 免费在线观看视频国产中文字幕亚洲| 久久亚洲精品不卡| 美女主播在线视频| 国产成+人综合+亚洲专区| av视频免费观看在线观看| 免费在线观看黄色视频的| 我要看黄色一级片免费的| 国产精品国产av在线观看| 最黄视频免费看| 久久精品亚洲熟妇少妇任你| 男女床上黄色一级片免费看| 麻豆国产av国片精品| 一本久久精品| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 天堂动漫精品| 国产成人免费观看mmmm| 激情在线观看视频在线高清 | 一进一出好大好爽视频| 欧美日韩福利视频一区二区| 国产老妇伦熟女老妇高清| 日本黄色日本黄色录像| 久久久久视频综合| 成人特级黄色片久久久久久久 | 久9热在线精品视频| 国产亚洲精品一区二区www | 久久精品aⅴ一区二区三区四区| 欧美乱妇无乱码| 亚洲午夜理论影院| 亚洲国产看品久久| 中亚洲国语对白在线视频| 国产成人av激情在线播放| videos熟女内射| a级毛片在线看网站| 少妇猛男粗大的猛烈进出视频| 国产主播在线观看一区二区| 操出白浆在线播放| 9热在线视频观看99| 岛国毛片在线播放| av又黄又爽大尺度在线免费看| 在线观看66精品国产| 欧美激情极品国产一区二区三区| 成年女人毛片免费观看观看9 | 久久午夜综合久久蜜桃| 99国产极品粉嫩在线观看| 久久久久久久久久久久大奶| 国产免费视频播放在线视频| av欧美777| 亚洲伊人久久精品综合| 一级黄色大片毛片| 亚洲中文av在线| 91麻豆av在线| 每晚都被弄得嗷嗷叫到高潮| 欧美在线一区亚洲| 中文字幕av电影在线播放| 国产在线观看jvid| 日本av免费视频播放| www日本在线高清视频| 人人妻人人爽人人添夜夜欢视频| 可以免费在线观看a视频的电影网站| 日韩免费av在线播放| 国产成人精品久久二区二区91| 国产成人av激情在线播放| 99re在线观看精品视频| 一进一出抽搐动态| 一边摸一边做爽爽视频免费| 国产伦人伦偷精品视频| 欧美久久黑人一区二区| 欧美精品一区二区免费开放| 两个人看的免费小视频| 国产精品麻豆人妻色哟哟久久| 欧美亚洲 丝袜 人妻 在线| 天天躁日日躁夜夜躁夜夜| 欧美黑人欧美精品刺激| 欧美一级毛片孕妇| 狠狠精品人妻久久久久久综合| a级毛片在线看网站| 天天躁日日躁夜夜躁夜夜| www.999成人在线观看| 色婷婷av一区二区三区视频| 黄色 视频免费看| 亚洲精品自拍成人| 色播在线永久视频| 久久精品国产亚洲av香蕉五月 | 久久久精品国产亚洲av高清涩受| 老司机在亚洲福利影院| 极品教师在线免费播放| 怎么达到女性高潮| av视频免费观看在线观看| 首页视频小说图片口味搜索| 精品亚洲乱码少妇综合久久| 国产黄频视频在线观看| 我要看黄色一级片免费的| 国产在线免费精品| av天堂久久9| 精品一区二区三卡| 黄色片一级片一级黄色片| 另类精品久久| 国产成人精品无人区| av电影中文网址| 人妻一区二区av| 高清av免费在线| 日韩欧美国产一区二区入口| 亚洲天堂av无毛| 亚洲精品成人av观看孕妇| 亚洲国产精品一区二区三区在线| 精品少妇久久久久久888优播| 久久天堂一区二区三区四区| 桃花免费在线播放| 成人av一区二区三区在线看| 这个男人来自地球电影免费观看| 不卡av一区二区三区| 黑丝袜美女国产一区| 精品国产乱子伦一区二区三区| 热99久久久久精品小说推荐| 亚洲精品国产一区二区精华液| www.自偷自拍.com| av福利片在线| 动漫黄色视频在线观看| 久久精品91无色码中文字幕| 亚洲午夜精品一区,二区,三区| 久久久国产精品麻豆| 高潮久久久久久久久久久不卡| 人妻 亚洲 视频| 亚洲综合色网址| 丁香六月天网| 亚洲国产av影院在线观看| 妹子高潮喷水视频| 丝袜人妻中文字幕| 九色亚洲精品在线播放| 18禁美女被吸乳视频| 在线观看舔阴道视频| 人人妻人人澡人人爽人人夜夜| 久久 成人 亚洲| 欧美黑人欧美精品刺激| 国产片内射在线| 日本五十路高清| 亚洲人成伊人成综合网2020| 999久久久国产精品视频| 99热国产这里只有精品6| 精品少妇一区二区三区视频日本电影| 中文字幕人妻丝袜制服| 9191精品国产免费久久| 五月开心婷婷网| 免费一级毛片在线播放高清视频 | 高清黄色对白视频在线免费看| 美国免费a级毛片| 成人亚洲精品一区在线观看| 人人澡人人妻人| 亚洲成人国产一区在线观看| 国产精品欧美亚洲77777| 99re在线观看精品视频| 亚洲精品乱久久久久久| 国产精品麻豆人妻色哟哟久久| 夜夜夜夜夜久久久久| 国产精品影院久久| 精品国产乱码久久久久久小说| 成人av一区二区三区在线看| 国产成人系列免费观看| 免费观看a级毛片全部| 婷婷成人精品国产| 日本av免费视频播放| 国产无遮挡羞羞视频在线观看| 久久性视频一级片| 国产精品亚洲一级av第二区| 视频在线观看一区二区三区| 欧美精品亚洲一区二区| 人成视频在线观看免费观看| 国产精品熟女久久久久浪| 国产伦人伦偷精品视频| 757午夜福利合集在线观看| 日本a在线网址| 久久精品亚洲熟妇少妇任你| 丝袜美腿诱惑在线| 我要看黄色一级片免费的| 黄色视频,在线免费观看| 亚洲色图综合在线观看| 亚洲中文av在线| 热re99久久国产66热| 国产欧美亚洲国产| 真人做人爱边吃奶动态| 欧美日本中文国产一区发布| 午夜免费鲁丝| 老汉色av国产亚洲站长工具| 久久99热这里只频精品6学生| 国产欧美日韩精品亚洲av| 午夜免费鲁丝| 在线天堂中文资源库| 久久 成人 亚洲| 午夜福利在线免费观看网站| 精品卡一卡二卡四卡免费| 国产av精品麻豆| 啦啦啦视频在线资源免费观看| 五月开心婷婷网| 99久久99久久久精品蜜桃| 国产av精品麻豆| av不卡在线播放| 久久久久久久精品吃奶| av免费在线观看网站| 桃红色精品国产亚洲av| 精品一区二区三区视频在线观看免费 | 久久久久久久久久久久大奶| 亚洲精品一卡2卡三卡4卡5卡| 久久亚洲真实| 9热在线视频观看99| 久久精品91无色码中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲黑人精品在线| 女同久久另类99精品国产91| 国产熟女午夜一区二区三区| 亚洲精品国产区一区二| 国产精品美女特级片免费视频播放器 | 91麻豆精品激情在线观看国产 | 十八禁网站免费在线| 丁香六月欧美| 亚洲av电影在线进入| 搡老熟女国产l中国老女人| 狠狠婷婷综合久久久久久88av| 两性午夜刺激爽爽歪歪视频在线观看 | 在线看a的网站| 9热在线视频观看99| 久久精品91无色码中文字幕| 女同久久另类99精品国产91| 国产成人影院久久av| 每晚都被弄得嗷嗷叫到高潮| 美女主播在线视频| 国产精品亚洲一级av第二区| 久久人妻熟女aⅴ| 欧美人与性动交α欧美软件| 最近最新中文字幕大全免费视频| 黄色怎么调成土黄色| 久久毛片免费看一区二区三区| 色尼玛亚洲综合影院| 高清视频免费观看一区二区| 看免费av毛片|