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

    砂泥巖儲(chǔ)層孔隙度預(yù)測(cè)的改進(jìn)模型

    2017-12-06 18:11:28鄧社根滕新保華桂錢王美均張宏兵
    關(guān)鍵詞:方法模型

    鄧社根, 滕新保, 華桂錢, 王美均, 張宏兵

    (1.江蘇省工程勘測(cè)研究院有限責(zé)任公司, 江蘇 揚(yáng)州225000; 2.河海大學(xué) 地球科學(xué)與工程學(xué)院, 南京210098)

    砂泥巖儲(chǔ)層孔隙度預(yù)測(cè)的改進(jìn)模型

    鄧社根1, 滕新保1, 華桂錢1, 王美均1, 張宏兵2

    (1.江蘇省工程勘測(cè)研究院有限責(zé)任公司, 江蘇 揚(yáng)州225000; 2.河海大學(xué) 地球科學(xué)與工程學(xué)院, 南京210098)

    針對(duì)支持向量回歸機(jī)(SVR)預(yù)測(cè)孔隙度精度不高的問(wèn)題,將線性擬合方法與SVR方法相結(jié)合,在模型建立過(guò)程中融入巖性信息,并為孔隙度預(yù)測(cè)值設(shè)定約束,建立一種新的儲(chǔ)層孔隙度預(yù)測(cè)模型。建模過(guò)程中,首先利用密度測(cè)井對(duì)孔隙度進(jìn)行線性擬合,然后在模型的輸入屬性中加入視巖性特征屬性,以實(shí)際孔隙度與孔隙度擬合值差值作為預(yù)測(cè)值建立SVR預(yù)測(cè)模型,最后利用密度測(cè)井與孔隙度的圖版關(guān)系對(duì)孔隙度預(yù)測(cè)值進(jìn)行約束。預(yù)測(cè)結(jié)果表明,采用綜合方法建立的孔隙度預(yù)測(cè)模型預(yù)測(cè)精度明顯高于傳統(tǒng)線性方法和SVR方法。該研究為儲(chǔ)層孔隙度預(yù)測(cè)提供了一種新方法。

    孔隙度; 預(yù)測(cè)模型; 支持向量回歸機(jī); 密度測(cè)井; 巖性特征屬性; 圖版約束

    0 引 言

    油氣勘探中,儲(chǔ)層孔隙度是表征儲(chǔ)層中孔隙體積的參數(shù)之一[1-2],同時(shí)也是建立各種地質(zhì)模型的必不可少的參數(shù)[3],因此,提高儲(chǔ)層孔隙度預(yù)測(cè)精度是很有必要的。影響儲(chǔ)層孔隙度的因素主要有顆粒粒徑的分布[4]、顆粒形狀、顆粒的排列方式、巖性類別、膠結(jié)物的類型及膠結(jié)程度等[5]。

    測(cè)井資料是預(yù)測(cè)儲(chǔ)層孔隙度的理想資料,常用的測(cè)井方法有自然伽馬、中子測(cè)井、聲波測(cè)井、密度測(cè)井等,不同的測(cè)井方法[6]都能一定程度上反映儲(chǔ)層孔隙度的變化。當(dāng)儲(chǔ)層孔隙度較大時(shí),孔隙度的變化對(duì)儲(chǔ)層巖石密度有很大影響,這時(shí)密度測(cè)井曲線與儲(chǔ)層孔隙度具有很好的線性相關(guān)性;然而隨著儲(chǔ)層孔隙度逐漸減小,孔隙度的變化對(duì)儲(chǔ)層密度影響力隨之降低,密度與孔隙度之間逐漸由線性相關(guān)轉(zhuǎn)變?yōu)榉蔷€性相關(guān)??紫抖纫彩怯绊懣v波速度的重要因素,同理縱波波速與孔隙度也存在上述關(guān)系。

    在孔隙度預(yù)測(cè)中,傳統(tǒng)線性擬合方法[7]是在測(cè)井屬性(如密度測(cè)井)與孔隙度之間建立線性擬合關(guān)系,根據(jù)求解的擬合關(guān)系式來(lái)計(jì)算孔隙度的大小。這種方法具有一定的可靠性(泛化性能),但其預(yù)測(cè)精度不高,這是由于測(cè)井屬性與儲(chǔ)層孔隙度之間并不是絕對(duì)線性相關(guān)的。

    近年來(lái),方中于等[8-11]提出使用非線性算法支持向量回歸機(jī)(SVR)對(duì)儲(chǔ)層參數(shù)進(jìn)行預(yù)測(cè),雖然預(yù)測(cè)精度有了很大的提高,但仍不能滿足實(shí)際需要。原因有三點(diǎn):其一,SVR沒(méi)有充分利用原始樣本空間中與孔隙度線性相關(guān)性好的測(cè)井屬性;其二,研究發(fā)現(xiàn),在影響孔隙度的眾多因素中,巖性對(duì)孔隙度影響最大,隨著儲(chǔ)層巖性不同,其孔隙類型往往不同,孔隙結(jié)構(gòu)也存在較大差異[12-13],而前人在孔隙度預(yù)測(cè)過(guò)程中沒(méi)有考慮儲(chǔ)層巖性的信息,直接將不同巖性的樣本混合在一起構(gòu)建一個(gè)預(yù)測(cè)模型;其三,SVR對(duì)輸入屬性進(jìn)行非線性變換,不僅無(wú)法知道確切的孔隙度擬合關(guān)系式,而且使得映射在高維空間的特征屬性沒(méi)有明確的物理含義,因此,模型的外推性能無(wú)法得到可靠的保證。

    為了充分利用實(shí)際資料中的已有信息,提高孔隙度的預(yù)測(cè)精度,保證模型的外推性能,筆者提出將傳統(tǒng)線性擬合方法與非線性算法SVR相結(jié)合,并在模型建立過(guò)程中加入巖性信息,在模型預(yù)測(cè)過(guò)程中設(shè)定孔隙度預(yù)測(cè)值與密度之間的圖版約束,由此建立新的儲(chǔ)層孔隙度預(yù)測(cè)模型,并以某油田三口井的測(cè)井?dāng)?shù)據(jù)展開(kāi)實(shí)例研究,力求找到精度高且方便快捷的孔隙度預(yù)測(cè)方法。

    1 孔隙度建模方法

    1.1支持向量回歸機(jī)方法原理

    支持向量機(jī) (Support Vector Machine, SVM)是針對(duì)經(jīng)典的二分類問(wèn)題提出的一種方法,它不僅可以應(yīng)用在分類領(lǐng)域,也可以應(yīng)用到回歸估計(jì)領(lǐng)域——支持向量回歸機(jī)(Support Vector Regression, SVR)。SVR模型中的樣本點(diǎn)不分類別,其輸出值是一組實(shí)數(shù)值,所尋求的最優(yōu)超平面也不能使兩類樣本點(diǎn)之間的分類間隔最大,而是使所有樣本點(diǎn)到最優(yōu)超平面的距離總和最小。算法的目標(biāo)與SVM相同,也是使兩條邊界線之間的距離盡可能的大,間隔內(nèi)容納更多的點(diǎn),這時(shí)最優(yōu)回歸超平面的求解等價(jià)于求解兩條邊界線的最大間隔。具體原理如下。

    作為訓(xùn)練樣本集

    T={(x1,y1),(x2,y2),…,(xi,yi),…,(xl,yl)},

    i=1,2,…,l,

    式中:l——樣本數(shù);

    xi——輸入向量,xin;

    yi——輸出值,yi。

    通常數(shù)據(jù)在原始樣本空間中線性函數(shù)回歸效果不佳,應(yīng)用映射函數(shù)Φ(·)對(duì)樣本進(jìn)行非線性變換,使其映射達(dá)到所需線性回歸精度的高維特征空間。為了簡(jiǎn)化運(yùn)算,SVR用核函數(shù)替換樣本高維特征空間的內(nèi)積,將原本應(yīng)在特征空間直接進(jìn)行的計(jì)算替換成在原始樣本空間中進(jìn)行,很大程度減少了程序的計(jì)算時(shí)間,使得該算法可以應(yīng)用到各個(gè)程序領(lǐng)域。目前,常用的核函數(shù)主要有以下三種:

    (1)多項(xiàng)式型

    K(xi·xj)=(p〈xi·xj〉+1)q。

    (1)

    (2)RBF型

    K(xi·xj)=exp{-γ‖xi-xj‖2}。

    (2)

    (3)S型

    K(xi·xj)=tanh(υ〈xi·xj〉+d),υ>0,d>0 。

    (3)

    式(1)~(3)中,p、q為多項(xiàng)式參數(shù),γ為徑向核參數(shù),υ是一個(gè)標(biāo)量,d為位移參數(shù)。

    非線性回歸函數(shù)可表示為

    f(x)=〈w,Φ(x)〉+b,

    (4)

    其中,Φ(x)是x映射到高維特征空間的映射函數(shù);b為函數(shù)的截距,b;w為函數(shù)的斜率,wn。根據(jù)凸二次規(guī)劃問(wèn)題可得到式(1)中的回歸函數(shù):

    (5)

    引入Lagrange函數(shù),并將問(wèn)題求解轉(zhuǎn)化為對(duì)偶形式,解得:

    (6)

    (7)

    (8)

    1.2SVR孔隙度預(yù)測(cè)模型的改進(jìn)

    1.2.1 測(cè)井屬性

    以國(guó)內(nèi)某油田測(cè)區(qū)三口井(井1、井2、井3)的測(cè)井?dāng)?shù)據(jù)為例,所使用的測(cè)井屬性有中子、自然伽馬、縱波速度、橫波速度、密度五種屬性。儲(chǔ)層巖性共有泥質(zhì)粉砂巖、粉砂巖、砂巖三種,其樣本數(shù)量如表1所示。各測(cè)井屬性與孔隙度的相關(guān)性圖版如圖1所示,其中φs為實(shí)際孔隙度。

    表1 巖性樣本數(shù)量

    從圖1可以看出,各測(cè)井屬性與孔隙度之間都有一定的線性相關(guān)性,除密度屬性與孔隙度的線性相關(guān)性相對(duì)明顯(圖1c),其他屬性與孔隙度的線性相關(guān)性效果都一般。當(dāng)孔隙度較大的時(shí)候,密度測(cè)井與孔隙度的相關(guān)性較好。

    a 自然伽馬

    b 中子測(cè)井

    c 密度測(cè)井

    d 縱波速度

    e 橫波速度

    1.2.2 改進(jìn)方法

    針對(duì)SVR砂泥巖儲(chǔ)層孔隙度預(yù)測(cè)模型,以下三個(gè)方面提出改進(jìn)措施:

    (1)線性方法與SVR結(jié)合

    傳統(tǒng)線性擬合方法利用密度測(cè)井與孔隙度的擬合關(guān)系建立預(yù)測(cè)模型,其外推性能穩(wěn)定;SVR方法建立孔隙度預(yù)測(cè)模型,預(yù)測(cè)精度較高。結(jié)合兩種方法的優(yōu)點(diǎn),提出將線性擬合方法和SVR結(jié)合,構(gòu)建一種新的孔隙度預(yù)測(cè)模型。首先建立密度測(cè)井與孔隙度的傳統(tǒng)線性擬合關(guān)系,用擬合關(guān)系式計(jì)算孔隙度的擬合值;然后,求出孔隙度剩余值(孔隙度剩余值=真實(shí)孔隙度-孔隙度擬合值);最后,利用SVR建立孔隙度剩余值的預(yù)測(cè)模型。

    (2)將巖性信息融入擬合模型

    巖性對(duì)孔隙度預(yù)測(cè)的影響較大[11-12],因此,文中在利用SVR建立孔隙度剩余值預(yù)測(cè)模型的過(guò)程中,引入巖性信息,以此對(duì)預(yù)測(cè)模型進(jìn)行改進(jìn)。但巖性是對(duì)巖芯使用文字表述的一種定性分類,只有通過(guò)對(duì)巖性資料的量化,得到一種與樣本一一對(duì)應(yīng)的數(shù)值,才能將巖性信息加入到預(yù)測(cè)模型中。文中將不同巖性設(shè)置為不同數(shù)值,如將砂巖類設(shè)置為數(shù)值1,粉砂巖設(shè)置為數(shù)值2,泥質(zhì)粉砂巖設(shè)置為數(shù)值3,所得到的屬性列命名為巖性特征屬性。實(shí)際過(guò)程中,鉆井的巖心資料很少,為了能夠預(yù)測(cè)整個(gè)井段的孔隙度值,提出首先在測(cè)井資料和少量巖性資料的基礎(chǔ)上,應(yīng)用支持向量機(jī)算法建立巖性識(shí)別模型,然后量化模型對(duì)測(cè)井巖性的分類結(jié)果,所得到的屬性列稱為視巖性特征屬性。

    (3)孔隙度預(yù)測(cè)值的條件約束

    在密度屬性與實(shí)際孔隙度所形成的圖版(圖1c)上,隨著孔隙度的減小,孔隙度對(duì)密度變化的影響力逐漸減弱,二者逐漸由線性相關(guān)轉(zhuǎn)變?yōu)榉蔷€性相關(guān),整體上樣本點(diǎn)基本包裹于兩條直線(H1、H2)之間,故孔隙度預(yù)測(cè)值與密度屬性在圖版上的投影點(diǎn)應(yīng)當(dāng)也在這兩條直線之間,可以使用斜交的兩條直線來(lái)約束孔隙度預(yù)測(cè)值。

    2 預(yù)測(cè)模型與結(jié)果分析

    2.1模型的建立

    建立預(yù)測(cè)模型的基本步驟如下:

    (1)樣本數(shù)據(jù)的選取及歸一化

    首先將樣本數(shù)據(jù)進(jìn)行歸一化處理,從中隨機(jī)抽取一半樣本數(shù)據(jù)作為訓(xùn)練樣本,剩余作為預(yù)測(cè)樣本。

    (2)視巖性特征屬性的求解

    使用RBF核函數(shù)構(gòu)建支持向量機(jī)巖性識(shí)別(SVC)模型。模型的構(gòu)建取決于兩個(gè)模型參數(shù):徑向核參數(shù)γ和懲罰參數(shù)C。運(yùn)用遺傳算法同步優(yōu)選這兩個(gè)參數(shù),優(yōu)選結(jié)果為γ=8.344 7,C=8.002 5。

    巖性分類準(zhǔn)確率如表2所示,量化分類結(jié)果,并進(jìn)行歸一化處理,即得到視巖性特征屬性。

    表2 巖性分類準(zhǔn)確率

    (3)孔隙度的線性擬合

    建立密度屬性與孔隙度之間線性擬合關(guān)系式,即

    φ=-0.645ρ+1.67,

    (9)

    其中,φ為孔隙度擬合值,ρ為密度屬性值。使用擬合關(guān)系式計(jì)算孔隙度擬合值,實(shí)際孔隙度與孔隙度擬合值之差即為孔隙度剩余值。

    (4)孔隙度剩余值預(yù)測(cè)模型的建立

    分別將各測(cè)井屬性(方法一)、各測(cè)井屬性與巖性特征屬性(方法二)、各測(cè)井屬性與視巖性特征屬性(方法三)作為輸入屬性,將孔隙度剩余值作為輸出值,建立孔隙度剩余值預(yù)測(cè)模型。文中利用孔隙度預(yù)測(cè)中最常用的RBF核函數(shù)構(gòu)建SVR預(yù)測(cè)模型,模型的構(gòu)建需確定懲罰參數(shù)C、不敏感參數(shù)ε和徑向核參數(shù)γ三個(gè)模型參數(shù)。運(yùn)用遺傳算法同步優(yōu)選三參數(shù),優(yōu)選結(jié)果為C=50、ε=0.011、γ=0.55。通過(guò)構(gòu)建的預(yù)測(cè)模型計(jì)算出孔隙度剩余值的預(yù)測(cè)值。

    (5)孔隙度預(yù)測(cè)值的條件約束

    將孔隙度剩余值的預(yù)測(cè)值與孔隙度擬合值求和,計(jì)算出孔隙度的預(yù)測(cè)值。在孔隙度與巖性屬性的圖版(圖1c)中,兩條邊界線的方程分別為:

    H1φ+1.11ρ-2.66=0,

    (10)

    H2φ+0.43ρ-1.24=0。

    (11)

    若樣本的孔隙度預(yù)測(cè)值與巖性屬性的投影點(diǎn)在兩條邊界線之外,將孔隙度預(yù)測(cè)值修正為該點(diǎn)密度測(cè)井屬性所對(duì)應(yīng)的邊界直線上的孔隙度值。

    2.2預(yù)測(cè)結(jié)果與分析

    傳統(tǒng)線性擬合、SVR及改進(jìn)方法(方法一、方法二、方法三+條件約束)的孔隙度預(yù)測(cè)精度如表3所示。

    表3 不同方法對(duì)儲(chǔ)層孔隙度的預(yù)測(cè)精度

    從表3可以看出,方法二的預(yù)測(cè)精度遠(yuǎn)高于方法一,說(shuō)明在孔隙度預(yù)測(cè)過(guò)程中加入巖性信息是有效的。方法三的預(yù)測(cè)精度雖略低于方法二,但依然遠(yuǎn)高于方法一,這是由于視巖性特征屬性是使用測(cè)井?dāng)?shù)據(jù)分類得到的結(jié)果,與巖性特征屬性還存在一定的誤差,因此,它只能近似地代表巖性特征屬性。這一方面說(shuō)明使用視巖性特征屬性代替巖性特證屬性是可行的,另一方面更加肯定了孔隙度預(yù)測(cè)模型中加入巖性信息的正確性。方法三+條件約束的孔隙度預(yù)測(cè)精度相對(duì)于方法三有了很大的提高,說(shuō)明文中提出的條件約束方法是有效的。

    非線性算法SVR對(duì)孔隙度的預(yù)測(cè)精度遠(yuǎn)高于線性擬合方法,這說(shuō)明SVR方法具有一定的優(yōu)越性;改進(jìn)方法(方法三+條件約束)對(duì)孔隙度的預(yù)測(cè)精度遠(yuǎn)高于傳統(tǒng)方法,證明文中提出的改進(jìn)方法是有效的。

    3 結(jié)束語(yǔ)

    (1)將傳統(tǒng)方法與SVR相結(jié)合,提出一種新的孔隙度預(yù)測(cè)模型。線性擬合方法、SVR方法及改進(jìn)方法對(duì)孔隙度的預(yù)測(cè)結(jié)果表明,文中改進(jìn)方法所建立的模型明顯提高了孔隙度的預(yù)測(cè)精度。

    (2)該模型將巖性信息融入到輸入屬性中,對(duì)孔隙度預(yù)測(cè)模型有一定的優(yōu)化效果。文中提出視巖性特征屬性的概念,解決了實(shí)際巖性數(shù)據(jù)少的問(wèn)題。

    (3)在孔隙度預(yù)測(cè)過(guò)程中,利用密度測(cè)井與孔隙度的圖版關(guān)系對(duì)孔隙度預(yù)測(cè)值進(jìn)行約束,使預(yù)測(cè)精度明顯提高。

    [1] 張晉言, 孫建孟. 利用測(cè)井資料評(píng)價(jià)泥頁(yè)巖油氣“五性”指標(biāo)[J]. 測(cè)井技術(shù), 2012, 36(2): 146-153.

    [2] 陳歡慶, 曹 晨, 梁淑賢, 等. 儲(chǔ)層孔隙結(jié)構(gòu)研究進(jìn)展[J]. 天然氣地球科學(xué), 2013, 24(2): 227-237.

    [3] 鄧少貴, 盧春利. 裂縫-孔隙雙孔介質(zhì)的雙側(cè)向測(cè)井解釋方法初探[J]. 測(cè)井技術(shù), 2010, 34(6): 537-541.

    [4] 黃思靜, 郎咸國(guó), 蘭葉芳, 等. 儲(chǔ)層孔隙度-滲透率關(guān)系曲線中的截止孔隙度與儲(chǔ)層質(zhì)量[J]. 成都理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 38(6): 593-602.

    [5] 史基安, 王 琪. 影響碎屑巖天然氣儲(chǔ)層物性的主要控制因素[J]. 沉積學(xué)報(bào), 1995, 13(2): 128-139.

    [6] 范銘濤, 沈全意, 吳 輝, 等. 復(fù)雜巖性裂縫-孔隙型儲(chǔ)層孔隙度計(jì)算方法研究[J]. 天然氣工業(yè), 2005, 25(5): 29-30.

    [7] 王鐵利, 施為兵. 利用密度測(cè)井計(jì)算孔隙度[J]. 煤炭技術(shù), 2010, 29(2): 151.

    [8] 方中于, 萬(wàn) 歡, 李 勇. W地區(qū)儲(chǔ)層孔隙度地震預(yù)測(cè)技術(shù)及應(yīng)用研究[J]. 礦物巖石, 2012, 32(3): 110-115.

    [9] 朱永才, 薛坤林. 基于支持向量機(jī)的儲(chǔ)層參數(shù)反演[J]. 茂名學(xué)院學(xué)報(bào), 2012, 22(1): 44-46.

    [10] 張彥周, 王春香. 支持向量機(jī)的測(cè)井?dāng)?shù)據(jù)預(yù)測(cè)儲(chǔ)層滲透率方法[J]. 科技傳播, 2009, 6(5): 83-85.

    [11] 陳 華, 鄧少貴, 范宜仁. 基于LS-SVM的測(cè)井物性參數(shù)的預(yù)測(cè)方法[J]. 計(jì)算機(jī)工程與應(yīng)用, 2007, 43(23): 208-210.

    [12] 賴 錦, 王貴文, 王書(shū)南, 等. 碎屑巖儲(chǔ)層成巖相研究現(xiàn)狀及進(jìn)展[J]. 地球科學(xué)進(jìn)展, 2013, 28(1): 39-50.

    [13] 閆國(guó)亮, 孫建孟, 劉學(xué)峰, 等. 儲(chǔ)層巖石微觀孔隙結(jié)構(gòu)特征及其對(duì)滲透率影響[J]. 測(cè)井技術(shù), 2014, 38(1): 28-32.

    (編校荀海鑫)

    Anovelpredictionmodelforporosityinsandshalereservoir

    DengShegen1,TengXinbao1,HuaGuiqian1,WangMeijun1,ZhangHongbing2

    (1.Jiangsu Province Engineering Investigation & Research Institute, Yangzhou 225000, China; 2.School of Earth Sciences & Engineering, Hohai University, Nanjing 210098, China)

    This paper proposes a novel prediction model as an alternative to Support Vector Regression (SVR) which suffers from an inaccurate porosity prediction. The model combining the linear fitting with the SVR works by adding the lithology information to the model establishment process and thereby constraining porosity prediction. The modeling process begins with a linear fitting of porosity using the density log, followed by the addition of lithology feature attribute into the input attributes and the subsequent development of SVR prediction model using the difference between the actual porosity and porosity fitting value as model predictions; and ends with constraining porosity prediction value using plates relations of density log and porosity. The results demonstrate that the proposed model affords a significantly higher prediction accuracy than the existing linear method and SVR method. This study may provide a new method for reservoir porosity prediction.

    porosity; prediction model; SVR; density logging; lithology feature attribute; plate constraint

    10.3969/j.issn.2095-7262.2017.06.010

    P618.13

    2095-7262(2017)06-0621-05

    A

    2017-05-26

    國(guó)家自然科學(xué)基金項(xiàng)目(41374116)

    鄧社根(1974-),男,江蘇省泰州人,高級(jí)工程師,研究方向:地震勘探、測(cè)井,E-mail:644968553@qq.com。

    猜你喜歡
    方法模型
    一半模型
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    學(xué)習(xí)方法
    可能是方法不對(duì)
    3D打印中的模型分割與打包
    用對(duì)方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    午夜福利在线观看吧| 亚洲精品成人久久久久久| 亚洲精品日本国产第一区| 国产视频首页在线观看| 国产真实伦视频高清在线观看| 黄色一级大片看看| 国产精品99久久久久久久久| 特级一级黄色大片| 免费大片黄手机在线观看| 日韩欧美一区视频在线观看 | 嫩草影院新地址| 国产精品三级大全| 听说在线观看完整版免费高清| 日韩av在线大香蕉| 精品一区二区三区视频在线| 超碰97精品在线观看| 亚洲国产成人一精品久久久| 嘟嘟电影网在线观看| 国产淫语在线视频| 熟妇人妻久久中文字幕3abv| 夜夜看夜夜爽夜夜摸| 午夜日本视频在线| 日韩欧美精品v在线| 亚洲av成人精品一二三区| 一级a做视频免费观看| 观看美女的网站| 一区二区三区免费毛片| 国产精品熟女久久久久浪| 91久久精品国产一区二区成人| 青春草亚洲视频在线观看| 女人十人毛片免费观看3o分钟| 久久久久九九精品影院| 一级黄片播放器| 中文字幕人妻熟人妻熟丝袜美| 精品久久久噜噜| 亚洲无线观看免费| 精品国产一区二区三区久久久樱花 | 一个人观看的视频www高清免费观看| 国产精品av视频在线免费观看| 天堂√8在线中文| 天美传媒精品一区二区| 色尼玛亚洲综合影院| 国产91av在线免费观看| 国产精品伦人一区二区| 国产精品久久久久久久久免| 国产精品伦人一区二区| 日韩欧美国产在线观看| 亚洲最大成人中文| 婷婷色麻豆天堂久久| 在线免费观看的www视频| 天美传媒精品一区二区| 国产精品久久久久久精品电影小说 | 三级男女做爰猛烈吃奶摸视频| 日本与韩国留学比较| 国产成人福利小说| 老司机影院成人| 女的被弄到高潮叫床怎么办| 国产伦精品一区二区三区四那| 亚洲成色77777| 日产精品乱码卡一卡2卡三| 日韩不卡一区二区三区视频在线| 国产不卡一卡二| 97超视频在线观看视频| 国产av不卡久久| or卡值多少钱| 亚洲四区av| 高清毛片免费看| 18+在线观看网站| 高清午夜精品一区二区三区| 亚洲精品成人久久久久久| 乱人视频在线观看| 最后的刺客免费高清国语| 午夜福利成人在线免费观看| 国产精品不卡视频一区二区| 国产永久视频网站| 亚洲自偷自拍三级| 毛片一级片免费看久久久久| 欧美成人精品欧美一级黄| 久久久久免费精品人妻一区二区| 午夜福利网站1000一区二区三区| 亚洲伊人久久精品综合| 日本爱情动作片www.在线观看| 99热6这里只有精品| 男女那种视频在线观看| 国产毛片a区久久久久| av网站免费在线观看视频 | 两个人视频免费观看高清| 99热6这里只有精品| 国产高清三级在线| 免费大片黄手机在线观看| 国产黄频视频在线观看| 美女大奶头视频| 99热网站在线观看| 日本爱情动作片www.在线观看| 秋霞在线观看毛片| 九九爱精品视频在线观看| 又爽又黄a免费视频| 国产成人一区二区在线| 99视频精品全部免费 在线| 亚洲成人精品中文字幕电影| 免费av毛片视频| 久久久成人免费电影| 美女黄网站色视频| 99久久中文字幕三级久久日本| 亚洲怡红院男人天堂| 国产精品久久久久久精品电影| 精品国产三级普通话版| 草草在线视频免费看| 91在线精品国自产拍蜜月| 非洲黑人性xxxx精品又粗又长| 亚洲精品456在线播放app| 免费看av在线观看网站| 国产淫语在线视频| 国产老妇女一区| 亚洲怡红院男人天堂| 亚洲精品乱久久久久久| 最近最新中文字幕免费大全7| 99热网站在线观看| 两个人视频免费观看高清| 熟女人妻精品中文字幕| 亚洲av在线观看美女高潮| 可以在线观看毛片的网站| 精品国产一区二区三区久久久樱花 | av.在线天堂| 国产精品综合久久久久久久免费| 神马国产精品三级电影在线观看| 日韩一本色道免费dvd| 久久精品夜夜夜夜夜久久蜜豆| 99久久精品热视频| 国产成人免费观看mmmm| 久久久久久久久久黄片| 一级片'在线观看视频| 欧美精品国产亚洲| 欧美精品国产亚洲| av一本久久久久| 欧美日韩视频高清一区二区三区二| 青青草视频在线视频观看| 久久午夜福利片| 禁无遮挡网站| 在现免费观看毛片| 大香蕉97超碰在线| 婷婷色综合www| 亚洲av一区综合| 午夜激情久久久久久久| 久久久精品欧美日韩精品| 视频中文字幕在线观看| 成人午夜高清在线视频| 久久久精品94久久精品| 丰满乱子伦码专区| 最近视频中文字幕2019在线8| 听说在线观看完整版免费高清| 欧美日韩视频高清一区二区三区二| 国产人妻一区二区三区在| 蜜桃久久精品国产亚洲av| 免费大片黄手机在线观看| 黄片无遮挡物在线观看| 亚洲av一区综合| 性插视频无遮挡在线免费观看| 最后的刺客免费高清国语| 亚洲四区av| 99热全是精品| 国产v大片淫在线免费观看| 国产在线一区二区三区精| 亚洲成色77777| 亚洲欧美成人综合另类久久久| 啦啦啦中文免费视频观看日本| 久久久久性生活片| 亚洲在线自拍视频| 高清午夜精品一区二区三区| 国产精品一及| 国产免费视频播放在线视频 | 在线观看av片永久免费下载| 亚洲欧洲国产日韩| av卡一久久| videos熟女内射| 欧美日韩综合久久久久久| 99热网站在线观看| 国内精品宾馆在线| 日韩视频在线欧美| 男女边摸边吃奶| 中国国产av一级| 亚洲一区高清亚洲精品| 国产一区二区三区综合在线观看 | 秋霞在线观看毛片| 国产在线男女| freevideosex欧美| 免费大片18禁| 晚上一个人看的免费电影| 久久久久久久久久成人| 深夜a级毛片| 欧美成人a在线观看| 亚洲不卡免费看| 亚洲无线观看免费| 婷婷色麻豆天堂久久| 久久久久久久午夜电影| 久久国内精品自在自线图片| videos熟女内射| 国产精品1区2区在线观看.| 国内揄拍国产精品人妻在线| 免费观看性生交大片5| 91精品国产九色| 看黄色毛片网站| 国产精品综合久久久久久久免费| 国产不卡一卡二| a级毛片免费高清观看在线播放| 啦啦啦中文免费视频观看日本| 国产淫片久久久久久久久| 免费看日本二区| 99热这里只有是精品50| 久久久久久久久久成人| 午夜福利在线观看免费完整高清在| 欧美日韩在线观看h| 精品午夜福利在线看| 深爱激情五月婷婷| 婷婷色综合www| 国产精品久久久久久久久免| 精品久久久噜噜| 成人亚洲欧美一区二区av| 男人舔奶头视频| 亚洲av男天堂| 一级二级三级毛片免费看| 日韩亚洲欧美综合| 婷婷色av中文字幕| 亚洲天堂国产精品一区在线| 亚洲欧美成人精品一区二区| 国产精品一区二区在线观看99 | 国产亚洲91精品色在线| 亚洲aⅴ乱码一区二区在线播放| 国产黄色视频一区二区在线观看| 天堂网av新在线| 一区二区三区高清视频在线| 精品久久久久久久久av| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 国产精品久久久久久久电影| 男女下面进入的视频免费午夜| 国产精品日韩av在线免费观看| 欧美成人a在线观看| 午夜激情欧美在线| 高清在线视频一区二区三区| 精品久久久久久久人妻蜜臀av| av福利片在线观看| 日本与韩国留学比较| 日韩三级伦理在线观看| 国产综合精华液| 国产精品一区www在线观看| 极品少妇高潮喷水抽搐| 欧美日本视频| 国产高清不卡午夜福利| 国产三级在线视频| 亚洲精品色激情综合| 免费观看无遮挡的男女| 秋霞在线观看毛片| 亚洲人成网站高清观看| 十八禁网站网址无遮挡 | 中文天堂在线官网| 亚洲国产精品专区欧美| 春色校园在线视频观看| 亚洲在线自拍视频| 成人亚洲欧美一区二区av| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 日本欧美国产在线视频| a级一级毛片免费在线观看| 日本色播在线视频| 91久久精品国产一区二区成人| 真实男女啪啪啪动态图| 中文欧美无线码| 高清欧美精品videossex| 国产成人午夜福利电影在线观看| 成人国产麻豆网| 国产精品熟女久久久久浪| 亚洲精华国产精华液的使用体验| 亚洲熟女精品中文字幕| 精品久久久精品久久久| 国产精品人妻久久久久久| 欧美bdsm另类| 精品欧美国产一区二区三| 国产午夜精品论理片| 日本wwww免费看| 在线 av 中文字幕| 小蜜桃在线观看免费完整版高清| 久久精品国产自在天天线| 丰满少妇做爰视频| 亚洲三级黄色毛片| 国产亚洲5aaaaa淫片| 搡老妇女老女人老熟妇| videos熟女内射| 亚洲精品久久久久久婷婷小说| 国产免费视频播放在线视频 | 日韩不卡一区二区三区视频在线| 好男人在线观看高清免费视频| 亚洲人成网站在线观看播放| 欧美日韩在线观看h| 国产精品久久久久久av不卡| 老师上课跳d突然被开到最大视频| 精品熟女少妇av免费看| 久久久成人免费电影| 国产人妻一区二区三区在| 午夜免费观看性视频| 男女下面进入的视频免费午夜| 好男人在线观看高清免费视频| 欧美人与善性xxx| 亚洲欧美成人综合另类久久久| 在线免费十八禁| kizo精华| 热99在线观看视频| 麻豆成人午夜福利视频| 干丝袜人妻中文字幕| 日韩,欧美,国产一区二区三区| www.色视频.com| 熟女人妻精品中文字幕| 黄色日韩在线| 噜噜噜噜噜久久久久久91| 亚洲不卡免费看| 日韩av不卡免费在线播放| 久久午夜福利片| 午夜视频国产福利| 免费av毛片视频| 少妇熟女欧美另类| h日本视频在线播放| 午夜激情久久久久久久| 天堂俺去俺来也www色官网 | 自拍偷自拍亚洲精品老妇| 丰满少妇做爰视频| 国产伦精品一区二区三区四那| 纵有疾风起免费观看全集完整版 | 久久这里只有精品中国| 春色校园在线视频观看| 国产精品人妻久久久久久| 99热网站在线观看| 99热这里只有是精品在线观看| 国产白丝娇喘喷水9色精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产黄a三级三级三级人| 免费看av在线观看网站| 寂寞人妻少妇视频99o| 国产一区二区亚洲精品在线观看| 国产精品久久久久久久久免| 国产熟女欧美一区二区| 搡老乐熟女国产| 色尼玛亚洲综合影院| 国产乱人视频| 一级毛片我不卡| 亚洲性久久影院| 我的女老师完整版在线观看| 简卡轻食公司| 免费无遮挡裸体视频| 亚洲国产欧美在线一区| 国产精品伦人一区二区| 免费无遮挡裸体视频| 精品久久国产蜜桃| 中文字幕制服av| 精品国产一区二区三区久久久樱花 | 国产一区二区亚洲精品在线观看| 国产亚洲一区二区精品| av福利片在线观看| 中国国产av一级| 熟妇人妻不卡中文字幕| 国产永久视频网站| 男女视频在线观看网站免费| 日韩电影二区| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| 91狼人影院| 亚洲怡红院男人天堂| 久久久久精品性色| 国产白丝娇喘喷水9色精品| 淫秽高清视频在线观看| 亚洲精品456在线播放app| 一级毛片黄色毛片免费观看视频| 亚洲av一区综合| av在线老鸭窝| 中文欧美无线码| 春色校园在线视频观看| 国产精品人妻久久久久久| 国产不卡一卡二| 成年人午夜在线观看视频 | 日韩av在线免费看完整版不卡| 亚洲aⅴ乱码一区二区在线播放| 国产综合懂色| av在线观看视频网站免费| 在线观看免费高清a一片| 青春草视频在线免费观看| 日韩中字成人| 亚洲av免费高清在线观看| 九草在线视频观看| 成人特级av手机在线观看| 在线观看人妻少妇| 精华霜和精华液先用哪个| 久久久久网色| 黄色日韩在线| 亚洲av免费在线观看| 亚洲欧美一区二区三区国产| 国产伦一二天堂av在线观看| 色综合亚洲欧美另类图片| 2021少妇久久久久久久久久久| 少妇的逼水好多| 久久国产乱子免费精品| 成年人午夜在线观看视频 | 国产欧美另类精品又又久久亚洲欧美| 国产精品蜜桃在线观看| 亚洲人与动物交配视频| 色综合色国产| 床上黄色一级片| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 日韩 亚洲 欧美在线| 亚洲电影在线观看av| 一个人免费在线观看电影| 国产午夜精品一二区理论片| 亚州av有码| 男人狂女人下面高潮的视频| 高清午夜精品一区二区三区| 一夜夜www| 国产一区有黄有色的免费视频 | 一级a做视频免费观看| 99热这里只有精品一区| 精品人妻偷拍中文字幕| 好男人视频免费观看在线| 国产伦在线观看视频一区| 七月丁香在线播放| 午夜精品国产一区二区电影 | 日韩,欧美,国产一区二区三区| av又黄又爽大尺度在线免费看| 国产精品蜜桃在线观看| 亚洲精品自拍成人| 日韩欧美精品免费久久| 久久精品国产亚洲av涩爱| 国产精品一区二区三区四区久久| 久久久久久久大尺度免费视频| 成人毛片60女人毛片免费| 天天躁日日操中文字幕| 成人欧美大片| 99热网站在线观看| 国产单亲对白刺激| 国产免费一级a男人的天堂| 色播亚洲综合网| kizo精华| 国产一区亚洲一区在线观看| 别揉我奶头 嗯啊视频| 免费高清在线观看视频在线观看| 婷婷色综合www| 国产午夜精品论理片| 国产精品三级大全| 99久久人妻综合| 亚洲自拍偷在线| av女优亚洲男人天堂| 成人亚洲精品一区在线观看 | 久久草成人影院| 国产极品天堂在线| 国产人妻一区二区三区在| 大话2 男鬼变身卡| 亚洲精品日本国产第一区| 麻豆久久精品国产亚洲av| 亚洲国产精品成人久久小说| 国产高潮美女av| 一本一本综合久久| 在线天堂最新版资源| av国产久精品久网站免费入址| 免费无遮挡裸体视频| 日日干狠狠操夜夜爽| 免费看光身美女| 观看美女的网站| 欧美日本视频| 一级毛片电影观看| 秋霞伦理黄片| 九色成人免费人妻av| 国产一区有黄有色的免费视频 | 色综合亚洲欧美另类图片| 亚洲一区高清亚洲精品| 2022亚洲国产成人精品| 97精品久久久久久久久久精品| 两个人的视频大全免费| 国产在视频线精品| 成年版毛片免费区| 亚洲aⅴ乱码一区二区在线播放| 亚洲,欧美,日韩| 色综合站精品国产| 天堂网av新在线| 七月丁香在线播放| 高清毛片免费看| 日本与韩国留学比较| 建设人人有责人人尽责人人享有的 | 亚洲国产最新在线播放| 中文天堂在线官网| 晚上一个人看的免费电影| 亚洲一级一片aⅴ在线观看| 看黄色毛片网站| 女人十人毛片免费观看3o分钟| 欧美三级亚洲精品| 欧美丝袜亚洲另类| 国产伦精品一区二区三区视频9| 国产色婷婷99| 嫩草影院新地址| av.在线天堂| 韩国av在线不卡| 国产伦精品一区二区三区视频9| 日韩电影二区| 日日干狠狠操夜夜爽| 亚洲电影在线观看av| 国产精品国产三级国产专区5o| 又大又黄又爽视频免费| 成年女人在线观看亚洲视频 | 中文字幕亚洲精品专区| 国产成人精品婷婷| 亚洲在久久综合| 禁无遮挡网站| 亚洲欧美精品专区久久| 天堂影院成人在线观看| 我要看日韩黄色一级片| 十八禁网站网址无遮挡 | 天堂中文最新版在线下载 | 伊人久久国产一区二区| a级毛片免费高清观看在线播放| 色综合色国产| 国产成人aa在线观看| 男女下面进入的视频免费午夜| 欧美97在线视频| 最近中文字幕高清免费大全6| 天堂中文最新版在线下载 | 乱码一卡2卡4卡精品| 国产精品爽爽va在线观看网站| 久久久久久久久久久免费av| 狂野欧美白嫩少妇大欣赏| 亚洲美女视频黄频| 久久久精品免费免费高清| 在线观看免费高清a一片| 22中文网久久字幕| 熟妇人妻不卡中文字幕| 色综合站精品国产| 夫妻性生交免费视频一级片| 在线免费观看不下载黄p国产| 亚洲美女视频黄频| 人人妻人人看人人澡| 久久午夜福利片| 精品人妻偷拍中文字幕| 亚洲图色成人| 我要看日韩黄色一级片| 在线天堂最新版资源| 青春草国产在线视频| 美女黄网站色视频| 99re6热这里在线精品视频| 亚洲四区av| 日本色播在线视频| 成人一区二区视频在线观看| 一区二区三区四区激情视频| 午夜福利网站1000一区二区三区| 亚洲精品国产av蜜桃| 亚洲国产欧美人成| av天堂中文字幕网| 麻豆精品久久久久久蜜桃| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本熟妇午夜| 亚洲国产精品成人久久小说| 久久久久精品久久久久真实原创| 成人午夜精彩视频在线观看| 国产黄片视频在线免费观看| 欧美精品一区二区大全| 九九爱精品视频在线观看| 色综合站精品国产| 亚洲av二区三区四区| 99热6这里只有精品| 免费电影在线观看免费观看| av在线亚洲专区| 日本一二三区视频观看| 精华霜和精华液先用哪个| 美女被艹到高潮喷水动态| 六月丁香七月| av在线观看视频网站免费| 一级爰片在线观看| 女人十人毛片免费观看3o分钟| 国产午夜福利久久久久久| 国产精品美女特级片免费视频播放器| av在线观看视频网站免费| 亚洲av电影不卡..在线观看| 高清日韩中文字幕在线| 亚洲丝袜综合中文字幕| 一个人免费在线观看电影| 久久久久久国产a免费观看| 白带黄色成豆腐渣| 三级毛片av免费| 亚洲av成人精品一二三区| 免费看不卡的av| 一区二区三区乱码不卡18| 男女边摸边吃奶| 久久99热6这里只有精品| 男女下面进入的视频免费午夜| 国产在线男女| 国产在视频线精品| 精品人妻偷拍中文字幕| 秋霞在线观看毛片| av免费观看日本| 国产在线男女| 最近手机中文字幕大全| 国产黄色免费在线视频| 精品99又大又爽又粗少妇毛片| 美女国产视频在线观看| 大又大粗又爽又黄少妇毛片口| 天美传媒精品一区二区| 国产综合懂色| 国产三级在线视频| 丝瓜视频免费看黄片| 中文资源天堂在线| 久久久久九九精品影院| 自拍偷自拍亚洲精品老妇| 欧美成人午夜免费资源| 国产成人福利小说| 国产伦在线观看视频一区| 欧美zozozo另类| 成人毛片60女人毛片免费| 特大巨黑吊av在线直播| 亚洲18禁久久av| 内射极品少妇av片p|