• <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
    賺錢方法
    亚洲av熟女| 热99在线观看视频| av在线亚洲专区| 久久精品久久久久久噜噜老黄 | 国产不卡一卡二| 精品国产三级普通话版| 美女国产视频在线观看| 精品久久久久久电影网 | 日本黄色片子视频| 欧美另类亚洲清纯唯美| 蜜桃久久精品国产亚洲av| 国产色爽女视频免费观看| 日韩一区二区视频免费看| 熟妇人妻久久中文字幕3abv| 国产精品1区2区在线观看.| 成年女人看的毛片在线观看| 亚洲激情五月婷婷啪啪| 久久精品影院6| 男人的好看免费观看在线视频| 欧美成人午夜免费资源| 小说图片视频综合网站| 婷婷六月久久综合丁香| 中文精品一卡2卡3卡4更新| 一区二区三区四区激情视频| 国产久久久一区二区三区| 国产探花极品一区二区| 免费观看a级毛片全部| 综合色av麻豆| a级毛色黄片| 久久久欧美国产精品| 永久免费av网站大全| 最近最新中文字幕免费大全7| 丝袜美腿在线中文| 欧美成人精品欧美一级黄| 色5月婷婷丁香| 国产极品天堂在线| 观看美女的网站| 青春草国产在线视频| 麻豆一二三区av精品| 中文天堂在线官网| 亚洲精品国产成人久久av| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av天美| 国产免费男女视频| 国内揄拍国产精品人妻在线| 变态另类丝袜制服| 性插视频无遮挡在线免费观看| 国产亚洲av嫩草精品影院| 女人被狂操c到高潮| 国产伦理片在线播放av一区| 少妇被粗大猛烈的视频| 成人无遮挡网站| 亚洲国产精品久久男人天堂| 人体艺术视频欧美日本| av免费在线看不卡| 久久久久网色| 免费人成在线观看视频色| 永久免费av网站大全| 长腿黑丝高跟| 国产精品国产三级国产专区5o | 最近中文字幕高清免费大全6| 成人性生交大片免费视频hd| 国产精品久久电影中文字幕| 麻豆精品久久久久久蜜桃| 日产精品乱码卡一卡2卡三| 三级国产精品欧美在线观看| 黄片无遮挡物在线观看| 午夜日本视频在线| 久久国产乱子免费精品| 爱豆传媒免费全集在线观看| 听说在线观看完整版免费高清| 欧美成人一区二区免费高清观看| 三级毛片av免费| 亚洲精品日韩av片在线观看| 你懂的网址亚洲精品在线观看 | 男女下面进入的视频免费午夜| 免费不卡的大黄色大毛片视频在线观看 | 你懂的网址亚洲精品在线观看 | 国产成人aa在线观看| av线在线观看网站| 久久久久久久久久成人| 国产精品久久久久久久电影| av.在线天堂| 国产老妇女一区| 只有这里有精品99| 亚洲国产欧洲综合997久久,| 国产探花极品一区二区| 成人一区二区视频在线观看| 白带黄色成豆腐渣| 午夜福利在线在线| 国产成人免费观看mmmm| 欧美性猛交╳xxx乱大交人| 性色avwww在线观看| 啦啦啦啦在线视频资源| 99热这里只有是精品50| 精品少妇黑人巨大在线播放 | 99九九线精品视频在线观看视频| 美女cb高潮喷水在线观看| 国产视频内射| 夜夜爽夜夜爽视频| 亚洲成色77777| 久久国产乱子免费精品| 一级毛片电影观看 | 亚洲中文字幕一区二区三区有码在线看| 国产高清视频在线观看网站| 色视频www国产| 精品久久久噜噜| 久久鲁丝午夜福利片| 色网站视频免费| 日本与韩国留学比较| 日韩三级伦理在线观看| 日韩成人av中文字幕在线观看| 国产真实乱freesex| 国产乱人偷精品视频| 日韩成人伦理影院| 日本黄大片高清| 草草在线视频免费看| 国产精品精品国产色婷婷| 亚洲av免费高清在线观看| 国产精品久久电影中文字幕| 国产亚洲精品av在线| 国产成人精品久久久久久| 男女下面进入的视频免费午夜| 乱人视频在线观看| 插阴视频在线观看视频| 一级av片app| 久久这里有精品视频免费| 熟女电影av网| 五月伊人婷婷丁香| 亚洲欧洲日产国产| 日日摸夜夜添夜夜爱| 久久99热6这里只有精品| 乱码一卡2卡4卡精品| 天天一区二区日本电影三级| 午夜福利高清视频| videos熟女内射| 免费人成在线观看视频色| 黄色一级大片看看| 欧美又色又爽又黄视频| 黑人高潮一二区| 亚洲五月天丁香| 三级国产精品片| av专区在线播放| 看十八女毛片水多多多| 只有这里有精品99| 久久久久久国产a免费观看| 18禁裸乳无遮挡免费网站照片| 在线免费观看不下载黄p国产| 最近最新中文字幕大全电影3| 久久精品国产亚洲av涩爱| 国产精品综合久久久久久久免费| 亚洲综合色惰| 国产成人a∨麻豆精品| 又爽又黄无遮挡网站| 久久久欧美国产精品| 91久久精品电影网| 爱豆传媒免费全集在线观看| 精品久久久久久久末码| 91精品国产九色| 国产激情偷乱视频一区二区| 国内精品美女久久久久久| or卡值多少钱| 国产久久久一区二区三区| 国产精品一区二区性色av| 国产亚洲91精品色在线| 国产黄色小视频在线观看| 2022亚洲国产成人精品| 国产不卡一卡二| 在线免费观看的www视频| 我要看日韩黄色一级片| 国产免费视频播放在线视频 | 日本免费a在线| 国产黄色视频一区二区在线观看 | 我的女老师完整版在线观看| 99国产精品一区二区蜜桃av| 99热这里只有是精品在线观看| 日韩欧美在线乱码| 国产乱人视频| 国产国拍精品亚洲av在线观看| 一区二区三区乱码不卡18| 国产在线一区二区三区精 | 亚洲美女视频黄频| 97在线视频观看| 极品教师在线视频| 91午夜精品亚洲一区二区三区| 日本爱情动作片www.在线观看| 国产综合懂色| eeuss影院久久| 日本欧美国产在线视频| 精品少妇黑人巨大在线播放 | 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 亚洲av电影不卡..在线观看| 男人舔奶头视频| 亚洲国产精品成人综合色| 日本一本二区三区精品| www日本黄色视频网| 国产在视频线精品| 成人高潮视频无遮挡免费网站| av又黄又爽大尺度在线免费看 | 男的添女的下面高潮视频| 搡老妇女老女人老熟妇| 中文字幕久久专区| 99久久精品热视频| 美女高潮的动态| 亚洲电影在线观看av| 偷拍熟女少妇极品色| 色尼玛亚洲综合影院| 国产乱人视频| 大香蕉97超碰在线| 色综合色国产| 听说在线观看完整版免费高清| 日韩欧美国产在线观看| 2022亚洲国产成人精品| 国产淫片久久久久久久久| 一个人看视频在线观看www免费| 日韩大片免费观看网站 | 国产精品综合久久久久久久免费| 一区二区三区乱码不卡18| 日韩在线高清观看一区二区三区| 国产成人精品婷婷| 国产在线一区二区三区精 | 99在线视频只有这里精品首页| 精品少妇黑人巨大在线播放 | 欧美人与善性xxx| 五月玫瑰六月丁香| 如何舔出高潮| 日本五十路高清| 六月丁香七月| 成人高潮视频无遮挡免费网站| 少妇的逼好多水| 九九爱精品视频在线观看| 日韩av在线免费看完整版不卡| 嫩草影院入口| 欧美日本视频| 亚洲国产欧美在线一区| 精品久久久久久久久久久久久| 亚洲精品久久久久久婷婷小说 | 网址你懂的国产日韩在线| 丰满乱子伦码专区| 精品久久久噜噜| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 欧美日韩一区二区视频在线观看视频在线 | 色5月婷婷丁香| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品乱码久久久v下载方式| 在线观看av片永久免费下载| 一区二区三区乱码不卡18| 伦精品一区二区三区| 亚洲18禁久久av| 女的被弄到高潮叫床怎么办| 熟女电影av网| 亚洲欧洲日产国产| 免费搜索国产男女视频| 免费看a级黄色片| 视频中文字幕在线观看| 亚洲不卡免费看| 你懂的网址亚洲精品在线观看 | 国产片特级美女逼逼视频| 一级毛片aaaaaa免费看小| 亚洲激情五月婷婷啪啪| 久久99热这里只频精品6学生 | 亚洲精品亚洲一区二区| 日韩 亚洲 欧美在线| 国产亚洲av片在线观看秒播厂 | 波多野结衣巨乳人妻| 国产精品国产三级专区第一集| 一个人免费在线观看电影| 亚洲久久久久久中文字幕| 亚洲中文字幕一区二区三区有码在线看| 国产一区二区亚洲精品在线观看| a级毛色黄片| 国产午夜精品论理片| 国产三级中文精品| 成人特级av手机在线观看| 天天躁日日操中文字幕| 国产69精品久久久久777片| 精品国内亚洲2022精品成人| 欧美区成人在线视频| 久久精品夜夜夜夜夜久久蜜豆| 中国国产av一级| 久久精品国产亚洲网站| 小蜜桃在线观看免费完整版高清| 亚洲精品aⅴ在线观看| 99久久中文字幕三级久久日本| videossex国产| 91精品国产九色| 小说图片视频综合网站| 国产成人精品婷婷| a级毛色黄片| 久久精品国产亚洲av涩爱| 国产成人免费观看mmmm| 熟女人妻精品中文字幕| 少妇高潮的动态图| 国产精品国产三级国产av玫瑰| 日韩亚洲欧美综合| 97超碰精品成人国产| 久久亚洲精品不卡| 91精品伊人久久大香线蕉| 国产欧美日韩精品一区二区| 少妇人妻一区二区三区视频| 一边摸一边抽搐一进一小说| 日韩一区二区视频免费看| 91午夜精品亚洲一区二区三区| 综合色av麻豆| 超碰97精品在线观看| 国产精品乱码一区二三区的特点| 久久久久精品久久久久真实原创| 欧美zozozo另类| 午夜久久久久精精品| 国产精品综合久久久久久久免费| 国产亚洲最大av| 日本黄色片子视频| 亚洲最大成人中文| 少妇被粗大猛烈的视频| 欧美激情在线99| 中文字幕久久专区| 婷婷色麻豆天堂久久 | 国产探花在线观看一区二区| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 国产探花极品一区二区| 久久这里只有精品中国| 精品一区二区三区视频在线| 中文资源天堂在线| av免费在线看不卡| 国产中年淑女户外野战色| 亚洲经典国产精华液单| 久久精品综合一区二区三区| 99热这里只有是精品在线观看| 亚洲色图av天堂| 亚洲18禁久久av| 久久草成人影院| 一个人看的www免费观看视频| 亚洲内射少妇av| 91精品国产九色| 美女脱内裤让男人舔精品视频| 国产一区二区在线av高清观看| 久久久久久久久久久丰满| 日韩一本色道免费dvd| 一二三四中文在线观看免费高清| 久久人妻av系列| 亚洲欧美精品自产自拍| 国产精品日韩av在线免费观看| 99热这里只有是精品在线观看| 亚洲av二区三区四区| 高清av免费在线| 久久草成人影院| av天堂中文字幕网| 久久99热这里只有精品18| 久久午夜福利片| 久久亚洲国产成人精品v| 久久精品夜夜夜夜夜久久蜜豆| 一区二区三区乱码不卡18| 又黄又爽又刺激的免费视频.| 亚洲国产精品sss在线观看| 国产精品久久视频播放| 中国美白少妇内射xxxbb| 少妇丰满av| 欧美激情在线99| 乱码一卡2卡4卡精品| 我的老师免费观看完整版| 最近视频中文字幕2019在线8| 中文字幕制服av| 你懂的网址亚洲精品在线观看 | 亚洲精品456在线播放app| av在线亚洲专区| 亚洲精品日韩av片在线观看| 精品欧美国产一区二区三| 一级毛片aaaaaa免费看小| 哪个播放器可以免费观看大片| 一区二区三区免费毛片| 国产精品野战在线观看| 日韩av不卡免费在线播放| 99国产精品一区二区蜜桃av| 国产高清三级在线| 亚洲av日韩在线播放| 国产视频首页在线观看| 中文字幕制服av| 亚洲第一区二区三区不卡| 亚洲av免费在线观看| 午夜老司机福利剧场| 国产亚洲av片在线观看秒播厂 | 日韩一区二区视频免费看| eeuss影院久久| 好男人在线观看高清免费视频| 小蜜桃在线观看免费完整版高清| 国产伦精品一区二区三区视频9| 国产精品嫩草影院av在线观看| 日本一本二区三区精品| 麻豆久久精品国产亚洲av| 亚洲av熟女| 97超碰精品成人国产| 国产成人精品婷婷| 亚洲av电影不卡..在线观看| 边亲边吃奶的免费视频| 草草在线视频免费看| 午夜爱爱视频在线播放| 欧美最新免费一区二区三区| 国产视频首页在线观看| 成人综合一区亚洲| 国产精品无大码| 国产成人aa在线观看| 2021天堂中文幕一二区在线观| 长腿黑丝高跟| 老司机福利观看| 免费观看在线日韩| 69人妻影院| 91精品一卡2卡3卡4卡| 有码 亚洲区| 亚洲天堂国产精品一区在线| 亚洲精品成人久久久久久| 男人和女人高潮做爰伦理| 成人高潮视频无遮挡免费网站| 国产欧美另类精品又又久久亚洲欧美| 欧美丝袜亚洲另类| 国产黄片美女视频| 高清av免费在线| 精品人妻一区二区三区麻豆| 日韩中字成人| 日韩一区二区三区影片| 成人毛片60女人毛片免费| 久久韩国三级中文字幕| 亚洲国产日韩欧美精品在线观看| 亚洲一级一片aⅴ在线观看| 欧美成人精品欧美一级黄| 大话2 男鬼变身卡| 亚洲五月天丁香| 欧美+日韩+精品| 麻豆国产97在线/欧美| 又爽又黄a免费视频| 五月伊人婷婷丁香| 人妻系列 视频| 国产免费福利视频在线观看| 亚洲精品影视一区二区三区av| 久久久久久久久久成人| 久久精品综合一区二区三区| 真实男女啪啪啪动态图| 免费播放大片免费观看视频在线观看 | 一边摸一边抽搐一进一小说| 欧美激情久久久久久爽电影| 成人漫画全彩无遮挡| 国产伦精品一区二区三区四那| 国产一区亚洲一区在线观看| 成人特级av手机在线观看| 国产精品av视频在线免费观看| 99久国产av精品| 97人妻精品一区二区三区麻豆| 国产精品一二三区在线看| 在线观看av片永久免费下载| 国产淫片久久久久久久久| 三级国产精品欧美在线观看| 好男人在线观看高清免费视频| av在线蜜桃| 亚洲国产精品成人久久小说| 全区人妻精品视频| 一区二区三区免费毛片| 免费无遮挡裸体视频| 91久久精品国产一区二区成人| 精品酒店卫生间| 视频中文字幕在线观看| 国产免费男女视频| 少妇裸体淫交视频免费看高清| 国产伦在线观看视频一区| 国产亚洲最大av| 美女被艹到高潮喷水动态| av又黄又爽大尺度在线免费看 | 欧美变态另类bdsm刘玥| 国产成人精品婷婷| 亚洲人成网站高清观看| 99热网站在线观看| 欧美成人精品欧美一级黄| 人人妻人人澡欧美一区二区| 简卡轻食公司| 午夜福利在线观看免费完整高清在| 国产成人a∨麻豆精品| 日韩,欧美,国产一区二区三区 | 国产午夜精品一二区理论片| 亚洲国产欧洲综合997久久,| 中文字幕人妻熟人妻熟丝袜美| 成人国产麻豆网| 国产91av在线免费观看| 久久精品久久久久久噜噜老黄 | 日韩精品有码人妻一区| 成年av动漫网址| 成人国产麻豆网| 亚洲国产精品国产精品| 七月丁香在线播放| 国产又黄又爽又无遮挡在线| 五月玫瑰六月丁香| 国产淫语在线视频| 久久精品综合一区二区三区| 少妇的逼水好多| 欧美性猛交黑人性爽| 亚洲精华国产精华液的使用体验| 综合色av麻豆| 亚洲无线观看免费| 狠狠狠狠99中文字幕| h日本视频在线播放| 精品久久久久久久末码| 精品不卡国产一区二区三区| 亚洲婷婷狠狠爱综合网| 夫妻性生交免费视频一级片| 能在线免费看毛片的网站| 亚洲中文字幕日韩| 一边摸一边抽搐一进一小说| 女人久久www免费人成看片 | 成年版毛片免费区| 欧美成人一区二区免费高清观看| 老司机影院毛片| 久久韩国三级中文字幕| 成年版毛片免费区| 69人妻影院| 少妇高潮的动态图| 一区二区三区四区激情视频| 精品少妇黑人巨大在线播放 | 内射极品少妇av片p| 天美传媒精品一区二区| 国产精品99久久久久久久久| 亚洲av一区综合| 不卡视频在线观看欧美| 一个人免费在线观看电影| 中文资源天堂在线| 成人漫画全彩无遮挡| 国产午夜精品论理片| 搡老妇女老女人老熟妇| 久久鲁丝午夜福利片| 69av精品久久久久久| 黄色一级大片看看| 日韩av在线大香蕉| 国产精品1区2区在线观看.| 欧美成人一区二区免费高清观看| 国产av码专区亚洲av| 毛片女人毛片| 国产老妇女一区| 插阴视频在线观看视频| 三级男女做爰猛烈吃奶摸视频| av在线播放精品| 精品少妇黑人巨大在线播放 | 欧美人与善性xxx| 免费大片18禁| 少妇熟女aⅴ在线视频| 69人妻影院| 免费av不卡在线播放| 日韩高清综合在线| 久久精品91蜜桃| 老司机影院成人| 天天躁夜夜躁狠狠久久av| 国产中年淑女户外野战色| av在线观看视频网站免费| 亚洲精品456在线播放app| 亚洲av成人精品一二三区| 午夜久久久久精精品| 爱豆传媒免费全集在线观看| 精品国内亚洲2022精品成人| 色综合站精品国产| 国产成人一区二区在线| 最近中文字幕高清免费大全6| 免费无遮挡裸体视频| 午夜老司机福利剧场| 寂寞人妻少妇视频99o| 国产精品美女特级片免费视频播放器| 婷婷色麻豆天堂久久 | 国产精品福利在线免费观看| 男人狂女人下面高潮的视频| 精品久久久久久久久亚洲| 亚洲av不卡在线观看| 亚洲精品成人久久久久久| 欧美另类亚洲清纯唯美| 亚洲美女视频黄频| 人人妻人人澡人人爽人人夜夜 | 视频中文字幕在线观看| 老司机影院毛片| 简卡轻食公司| 久久久午夜欧美精品| 女的被弄到高潮叫床怎么办| 国产精品一及| 亚洲av中文字字幕乱码综合| 在线观看美女被高潮喷水网站| 久久久久久久久久久免费av| 好男人在线观看高清免费视频| 免费看a级黄色片| 99久久人妻综合| 色哟哟·www| 成人高潮视频无遮挡免费网站| 最近手机中文字幕大全| 亚洲精品日韩在线中文字幕| 亚洲av电影不卡..在线观看| 日韩三级伦理在线观看| 18+在线观看网站| 一个人观看的视频www高清免费观看| 夫妻性生交免费视频一级片| 日本一二三区视频观看| 麻豆久久精品国产亚洲av| 51国产日韩欧美| 国产 一区精品| 别揉我奶头 嗯啊视频| 亚洲自偷自拍三级| 看黄色毛片网站| 国产精品女同一区二区软件| 精品少妇黑人巨大在线播放 | 内射极品少妇av片p| 一二三四中文在线观看免费高清| 波野结衣二区三区在线| 大香蕉97超碰在线| 亚洲精品亚洲一区二区| 国产精品伦人一区二区| 免费av毛片视频| 欧美潮喷喷水| 一个人免费在线观看电影| 欧美最新免费一区二区三区| 国产片特级美女逼逼视频| 亚洲成人精品中文字幕电影| 日本三级黄在线观看|