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

    基于支持向量回歸的油藏瞬態(tài)壓力擬合

    2020-12-05 02:41:40許恩華查文舒李道倫
    關(guān)鍵詞:試井油藏向量

    許恩華, 查文舒, 李道倫, 陳 剛

    (合肥工業(yè)大學(xué) 數(shù)學(xué)學(xué)院,安徽 合肥 230601)

    可靠的地質(zhì)模型是油藏開發(fā)、管理的基礎(chǔ)之一,也是可靠預(yù)測油井的保證。油藏建模無論是對工作人員知識和經(jīng)驗(yàn)的要求,還是對資料的要求都很高;且建模工程不確定因素多、情況復(fù)雜,使得地質(zhì)模型與油藏真實(shí)情況之間存在巨大差異。因此,需要使用已有資料對地質(zhì)模型進(jìn)行修正,尤其是包括了基于生產(chǎn)響應(yīng)的地質(zhì)信息的動態(tài)資料,這對于全面、準(zhǔn)確評價(jià)地質(zhì)模型非常重要。評價(jià)地質(zhì)模型與真實(shí)油藏之間的差異是對比動態(tài)數(shù)據(jù)與地質(zhì)模型計(jì)算的結(jié)果,即歷史擬合。動態(tài)歷史擬合方法是目前普遍采取的方式。

    油藏工程中的歷史擬合是一個(gè)典型的反問題[1]。與利用已知模型直接推斷結(jié)果的正問題不同,反問題是利用已知結(jié)果反向推斷模型參數(shù)再進(jìn)行預(yù)測。由于不同參數(shù)值會得到同樣的好的擬合結(jié)果,歷史擬合的過程中會出現(xiàn)不適定性。

    獲得多個(gè)可靠模型且準(zhǔn)確預(yù)測油藏情況是油藏工程和油藏?cái)?shù)值模擬技術(shù)不可缺少的一部分,自動歷史擬合作為其中的關(guān)鍵環(huán)節(jié),其問題復(fù)雜、工程量和計(jì)算量巨大,一直是該領(lǐng)域內(nèi)的研究熱點(diǎn)[2-4]。

    早期的傳統(tǒng)歷史擬合是人工調(diào)節(jié)參數(shù),即油藏工程師憑借自身的經(jīng)驗(yàn)手動不斷地調(diào)整1個(gè)或者幾個(gè)參數(shù)值,并對產(chǎn)生的模擬計(jì)算值與實(shí)際觀測值之間的誤差進(jìn)行分析,從而達(dá)到最小誤差,又被稱為人工歷史擬合。人工歷史擬合不僅效率低,而且難以保證結(jié)果的準(zhǔn)確性,并且無法滿足油田開采發(fā)展提出的海量參數(shù)的需求。

    在油田需求增加、計(jì)算機(jī)和數(shù)學(xué)工具發(fā)展的背景下,自動擬合技術(shù)得到迅速發(fā)展[5-6]。自動歷史擬合是利用計(jì)算機(jī)、優(yōu)化算法和擬合算法自動地調(diào)整油藏工程模型參數(shù),使得模擬模型與真實(shí)油藏模型之間的差異在可接受的范圍內(nèi)。從最初基于線性回歸的多項(xiàng)式擬合和優(yōu)化算法結(jié)合的擬合技術(shù)[7]發(fā)展到基于貝葉斯方法[8]、響應(yīng)面方法[9]、卡爾曼濾波方法[10]的自動歷史擬合方法,自動歷史擬合技術(shù)利用計(jì)算機(jī)解決了傳統(tǒng)歷史擬合方法中無法解決的問題。橫向比較自動歷史擬合中采用的各種方法,在收斂速度、精確性、非線性擬合、多參數(shù)擬合方面各有優(yōu)劣。

    響應(yīng)面方法由數(shù)學(xué)方法和統(tǒng)計(jì)方法結(jié)合而成,可以對多個(gè)變量影響下生產(chǎn)響應(yīng)的問題進(jìn)行建模和分析,并得到體現(xiàn)生產(chǎn)響應(yīng)與變量之間關(guān)系的顯示表達(dá)式。20世紀(jì)90年代,響應(yīng)面法被引入油藏工程,通過構(gòu)造響應(yīng)面模型并定義反映模擬計(jì)算值與實(shí)測值之間偏差的目標(biāo)函數(shù),響應(yīng)面方法顯著提升了地質(zhì)模型的準(zhǔn)確性。

    貝葉斯方法起源于英國學(xué)者貝葉斯提出的貝葉斯公式和由其發(fā)展出的歸納推理方法,基本思想是利用貝葉斯公式將已知的概率密度參數(shù)和先驗(yàn)概率轉(zhuǎn)換成后驗(yàn)概率,據(jù)此進(jìn)行決策分類。在地質(zhì)模擬和歷史擬合的應(yīng)用過程中,貝葉斯方法在模型推斷、在線預(yù)測方面有顯著的優(yōu)勢[8]。相比于傳統(tǒng)統(tǒng)計(jì)推斷方法,貝葉斯方法能利用計(jì)算機(jī)處理更大的數(shù)據(jù)集。先驗(yàn)信息的收集是貝葉斯方法中的重要部分。

    卡爾曼濾波方法是一種基于統(tǒng)計(jì)和蒙特卡洛方法生成集合的平均梯度類方法,其中最廣泛采用的是集合卡爾曼濾波方法。卡爾曼濾波方法的引入為描述油藏情況及動態(tài)預(yù)測中評價(jià)不確定性提供了一種新的方法。卡爾曼濾波方法具有計(jì)算效率高、魯棒性強(qiáng)的特點(diǎn),在吸收動態(tài)數(shù)據(jù)方面能做到持續(xù)、及時(shí)、迅速,有效減少計(jì)算時(shí)間,非常適合大規(guī)模非線性實(shí)時(shí)動態(tài)油藏模型更新的歷史擬合,在計(jì)算機(jī)多核進(jìn)行并行計(jì)算的條件下,能相互獨(dú)立計(jì)算多個(gè)油藏模型并進(jìn)行生產(chǎn)預(yù)測。模型更新的過程中,完全獨(dú)立于所使用的模型方程,不需要對不同的油藏模擬器進(jìn)行復(fù)雜的敏感計(jì)算,無需關(guān)注過程,在程序的開發(fā)、調(diào)試、維護(hù)方面有非常大的優(yōu)勢。

    支持向量機(jī)(support vector machines,SVM)是一種基于結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則(risk minimization inductive principle)尋找集合中最優(yōu)元素的機(jī)器學(xué)習(xí)算法。結(jié)構(gòu)風(fēng)險(xiǎn)最小化涉及到同時(shí)最小化經(jīng)驗(yàn)風(fēng)險(xiǎn)和VC維數(shù)。SVM主要有支持向量分類機(jī)(support vector classifier,SVC)和支持向量回歸機(jī)(support vector regression,SVR) 2種。SVR在1997年被提出,與傳統(tǒng)統(tǒng)計(jì)學(xué)中依據(jù)于所有訓(xùn)練數(shù)據(jù)與預(yù)測數(shù)據(jù)之間最小偏離的回歸方法不同,支持向量回歸依據(jù)的是所有數(shù)據(jù)最小誤差界,因此支持向量回歸只需要訓(xùn)練數(shù)據(jù)中的一部分子集即可,計(jì)算效率有了巨大提升的同時(shí)還保證了準(zhǔn)確性。經(jīng)過多年的發(fā)展,支持向量回歸在字體和物體識別、回歸和時(shí)間序列預(yù)測等方面取得了顯著的效果。

    可靠的地質(zhì)模型是解決石油開采工程問題、進(jìn)行產(chǎn)量預(yù)測的關(guān)鍵。試井分析是以測量井的井底壓力、溫度、流量等數(shù)據(jù)及其隨時(shí)間的變化并結(jié)合擬合方法分析地層和井筒參數(shù)的方法,其為油藏工程開發(fā)提供準(zhǔn)確的地層數(shù)據(jù),為油藏動態(tài)特性和各種生產(chǎn)方式下的生產(chǎn)變化趨勢提供可靠依據(jù)。經(jīng)過多年的發(fā)展,該方法已經(jīng)成為重要的油藏描述工具。

    數(shù)值試井是近些年發(fā)展起來的一項(xiàng)試井解釋技術(shù)[11-14],不同于基于滲流方程解析解的解析試井,數(shù)值試井是通過大量數(shù)值模擬來精確描述物理過程,應(yīng)用場景更廣。作為一項(xiàng)新技術(shù),數(shù)值試井存在計(jì)算參數(shù)多、計(jì)算時(shí)間長、需要手動調(diào)節(jié)不確定參數(shù)等困難。解釋一口井花費(fèi)的時(shí)間跨度在幾周到幾個(gè)月之間,需花費(fèi)大量的時(shí)間和人工成本,而自動調(diào)節(jié)不確定參數(shù)的方法可以大大減少人員工作強(qiáng)度。

    現(xiàn)有的自動反演方法適用范圍有限。響應(yīng)面方法在多變量情況下的處理效果不好;集合卡爾曼濾波方法的計(jì)算量大、計(jì)算時(shí)間長。為此,本文基于支持向量回歸,提出支持向量回歸的自動反求參數(shù)的試井分析方法和利用多組數(shù)據(jù)衡量擬合效果的整體誤差,并對試井?dāng)?shù)據(jù)進(jìn)行模擬、驗(yàn)證。

    1 理論背景

    1.1 試井分析評價(jià)的目標(biāo)函數(shù)

    在單井試井分析中,需要通過“試湊”一組合適的“地層參數(shù)和井筒參數(shù)”使得滲流方程的計(jì)算壓力逼近實(shí)測壓力,這顯然是一個(gè)反問題。一般采用目標(biāo)函數(shù)OF反映隨時(shí)間變化的計(jì)算壓力、壓力導(dǎo)數(shù)與井底實(shí)測壓力之間的偏差,通常采用一個(gè)簡單的平方和,即

    其中,p為壓力;f為壓力降落;d為壓力導(dǎo)數(shù);i、j、k分別為壓力、壓力降落、壓力導(dǎo)數(shù)的時(shí)間點(diǎn)序號;l、m、n分別為實(shí)測數(shù)據(jù)時(shí)間點(diǎn)的總數(shù)據(jù);obs、sim代表實(shí)測數(shù)據(jù)值和模擬值;x為不確定參數(shù);ω為權(quán)重。

    在油氣工程和試井分析的發(fā)展過程中,逐漸發(fā)現(xiàn)僅僅靠1組觀測值是不能衡量模型誤差,多組觀測值能更加有效衡量模型誤差,因此提出基于所有測量值的誤差函數(shù),即

    其中,α為觀測值組數(shù);τ為觀測值序號。

    1.2 SVM和支持向量回歸

    假設(shè)訓(xùn)練數(shù)據(jù)集是{(x1,y1),…,(xl,yl)}?A×R,A代表數(shù)據(jù)點(diǎn)的模式。ε-SVR的目的是尋找f(x),使得盡可能多的數(shù)據(jù)點(diǎn)(xi,yi)包含在ε邊界內(nèi)。線性函數(shù)f(x)定義為:

    f(x)=〈ω,x〉+b,ω∈A,b∈R

    (1)

    服從

    (2)

    (3)

    (4)

    (5)

    (6)

    將(4)~(6)式代入對偶問題(3)式可得:

    服從

    (5)式可寫為:

    (7)

    并且有:

    (8)

    這就是支持向量擴(kuò)張,即ω可被xi線性表示,因此求f(x)時(shí)不需計(jì)算ω。由KKT條件求解得到b的求解式為:

    b=yi-〈ω,xi〉-ε,αi∈(0,C)

    (9)

    (10)

    通過同樣的推導(dǎo)方法可以得到引入核函數(shù)的ω表達(dá)式為:

    (11)

    分割函數(shù)f(x)的表達(dá)式為:

    f(x)=sgn(〈ω,x〉+b)=

    (12)

    常用的核函數(shù)如下:

    (1) 多項(xiàng)式核函數(shù),即

    k(μi,μj)=(μi·μjd)。

    (2) 非勻次多項(xiàng)式核函數(shù),即

    k(μi,μj)=(μi·μj+1)d。

    (3) 高斯核函數(shù),即

    k(μi,μj)=exp(-γ‖μi-μj‖2),γ>0,

    (4) 雙曲正切核函數(shù),即

    k(μi,μj)=tanh(kμi·μj+c),

    有時(shí)k>0,c<0。

    1.3 BFGS算法

    BFGS算法是一種梯度類優(yōu)化算法,用于搜索函數(shù)的極值點(diǎn)。

    由牛頓法和擬牛頓條件可知:

    (13)

    yk=Hk+1sk

    (14)

    其中,Hk為海森矩陣;gk為函數(shù)f(x)在xk處的梯度;sk=xk+1-xk;yk=gk+1-gk。

    BFGS算法的基本思想是構(gòu)造一個(gè)矩陣直接逼近海森矩陣,即Bk≈Hk。設(shè)迭代格式為:

    Bk+1=Bk+ΔBk,k=0,1,2,…

    (15)

    其中,B0常取為單位矩陣。通過待定系數(shù)法可以求得Bk的迭代公式為:

    (16)

    下面給出BFGS算法的一個(gè)完整描述。

    (1) 給定初值x0和精度閾值ε,并令B0=I,k∶=0。

    (3) 利用線性搜索方法得到步長αk,令sk=αkpk,xk+1∶=xk+sk。

    (4) 若‖gk+1‖<ε,則算法結(jié)束。

    (5) 計(jì)算yk=gk+1-gk。

    (7) 令k∶=k+1,轉(zhuǎn)至步驟(2)。

    (17)

    進(jìn)一步展開可得到:

    (18)

    1.4 方法與流程

    (1) 確定試井?dāng)M合所用不確定參數(shù)及取值范圍。

    (2) 利用抽樣算法抽取若干組試算算例。

    (3) 將每組算例放入試井軟件中進(jìn)行計(jì)算,獲得壓力、壓力變化、壓力導(dǎo)數(shù)等數(shù)據(jù)。

    (4) 使用基于不同核函數(shù)的支持向量回歸模型對數(shù)據(jù)進(jìn)行擬合,得到回歸表達(dá)式和整體誤差。

    (5) 通過整體誤差判斷最優(yōu)的支持向量回歸模型,并作為試井分析所用模型。

    (6) 定義試井分析評價(jià)的目標(biāo)函數(shù),并利用BFGS算法對其進(jìn)行優(yōu)化,求得最優(yōu)的不確定參數(shù)值。

    (7) 試井軟件計(jì)算(6)式中所求的參數(shù)值組,求得壓力值等數(shù)據(jù),并與真實(shí)壓力值對比;若滿足誤差,則計(jì)算結(jié)束,否則返回步驟(1)。

    2 模型描述與案例研究

    五點(diǎn)井網(wǎng)兩相模型如圖1所示。

    圖1 油藏形狀及井分布

    油藏大小為600 m×400 m×10 m,厚度為10 m,孔隙度為0.2。圖1中,井1為生產(chǎn)井;井2、井3、井4、井5為注入井;井1、井2、井3周邊存在徑向復(fù)合區(qū)域,所有井都開井80 d,然后關(guān)井2 d。生產(chǎn)井生產(chǎn)流量為40 m3/d,注入井注入流量均為10 m3/d。

    采用PEBI網(wǎng)格劃分技術(shù),網(wǎng)格如圖2所示。其中,實(shí)驗(yàn)平臺是Interl?CoreTMi7-7 700 K CPU@4.20 GHz,實(shí)驗(yàn)環(huán)境是64位Windows7操作系統(tǒng)。

    圖2 油藏網(wǎng)格劃分

    本文方法選定的9個(gè)不確定性參數(shù)為:地層滲透率K,范圍為(100 mD,1 000 mD);5個(gè)區(qū)域滲透率K1、K2、K3、K4、K5,范圍為(100 mD,1 000 mD);生產(chǎn)井和注入井表皮因子S1、S2、S3,范圍為(-1,3)。本文的研究在注采平衡的前提下進(jìn)行,擬合對象是壓力恢復(fù)段生產(chǎn)井壓力曲線、壓力變化和壓力導(dǎo)數(shù)雙對數(shù)曲線。

    利用抽樣算法隨機(jī)抽取2 000組試算算例,所有算例都滿足滲流方程。將所有算例放入數(shù)值試井模擬器中計(jì)算壓力值等數(shù)據(jù)。2 000組試算算例中的前500組被假設(shè)為觀測值,其余1 500組用來構(gòu)造支持向量回歸模型并定義整體誤差,不同模型的整體誤差見表1所列。

    表1 模型整體誤差

    比較整體誤差后,選取三次支持向量回歸模型(Cubic)作為試井分析所使用的擬合模型。隨機(jī)選取初始點(diǎn),使用BFGS算法進(jìn)行優(yōu)化,得到不確定參數(shù)的估算值,誤差判斷標(biāo)準(zhǔn)是使目標(biāo)函數(shù)OF盡可能小。將2 000組數(shù)據(jù)的第1組作為真實(shí)值,由此得到估算值和真實(shí)值所對應(yīng)的壓力恢復(fù)段的井底壓力的對比、壓力變化及壓力導(dǎo)數(shù)的對比,如圖3所示。

    擬合不同核函數(shù)支持向量回歸模型的耗時(shí)見表2所列。由表2可知,使用BFGS算法進(jìn)行優(yōu)化耗時(shí)為90.488 132 s。可見支持向量回歸模型可以準(zhǔn)確確定試井分析所應(yīng)采用的擬合模型,并能準(zhǔn)確反演正確的地層參數(shù)及井筒參數(shù)。

    圖3 基于支持向量回歸方法的自動擬合結(jié)果對比

    表2 擬合支持向量回歸模型耗時(shí)

    3 結(jié) 論

    本文提出了基于支持向量回歸和BFGS優(yōu)化算法的試井自動反演方法,主要結(jié)論如下:

    (1) 基于實(shí)測的井底壓力、壓力降落、壓力導(dǎo)數(shù),支持向量回歸擬合方法和BFGS優(yōu)化方法能有效地進(jìn)行地層參數(shù)反演。

    (2) 通過對比不同核函數(shù)下的支持向量回歸模型誤差發(fā)現(xiàn),三次支持向量回歸模型最適宜估計(jì)試井分析模型。

    (3) 本文提出的試井自動反演方法能大大提高現(xiàn)場人員的工作效率,具有很好的應(yīng)用前景。

    另外,本文方法還有待于深入研究,例如方法的穩(wěn)定性以及抽樣方法、優(yōu)化方法等對反演結(jié)果的影響。

    猜你喜歡
    試井油藏向量
    向量的分解
    海上邊水氣藏利用試井資料確定水侵狀況研究
    海洋石油(2021年3期)2021-11-05 07:42:26
    聚焦“向量與三角”創(chuàng)新題
    低滲油藏注采結(jié)構(gòu)調(diào)整與挖潛技術(shù)
    云南化工(2020年11期)2021-01-14 00:51:02
    基于模糊數(shù)學(xué)的油藏干層識別研究
    反褶積試井技術(shù)在計(jì)算低滲儲層滲透率中的應(yīng)用
    向量垂直在解析幾何中的應(yīng)用
    向量五種“變身” 玩轉(zhuǎn)圓錐曲線
    注CO2混相驅(qū)候選油藏篩選評價(jià)新方法
    致密油藏多級壓裂水平井?dāng)?shù)值模擬及應(yīng)用
    国产熟女午夜一区二区三区| 欧美成人性av电影在线观看| 国产av一区二区精品久久| 在线观看66精品国产| 欧美久久黑人一区二区| 精品久久久久久久久久久久久| 一级毛片女人18水好多| 亚洲欧美日韩东京热| 亚洲熟女毛片儿| 在线观看舔阴道视频| 韩国av一区二区三区四区| 亚洲国产精品合色在线| 久久精品国产综合久久久| 白带黄色成豆腐渣| 午夜福利在线在线| 一区二区三区高清视频在线| 久久人人精品亚洲av| 两个人免费观看高清视频| 在线免费观看的www视频| 亚洲av电影在线进入| 免费电影在线观看免费观看| 叶爱在线成人免费视频播放| 又大又爽又粗| 亚洲欧美精品综合久久99| 最近最新中文字幕大全电影3| 国产激情久久老熟女| 看免费av毛片| 黄色成人免费大全| 2021天堂中文幕一二区在线观| 国产高清激情床上av| 久久久久久国产a免费观看| 国产高清有码在线观看视频 | 久久午夜亚洲精品久久| 1024视频免费在线观看| 国产日本99.免费观看| 色尼玛亚洲综合影院| 黄色毛片三级朝国网站| 国产精品野战在线观看| 精品一区二区三区视频在线观看免费| 五月伊人婷婷丁香| 亚洲精品粉嫩美女一区| а√天堂www在线а√下载| 99国产精品一区二区蜜桃av| 无遮挡黄片免费观看| 亚洲黑人精品在线| 国产精品免费一区二区三区在线| 在线十欧美十亚洲十日本专区| 国产成+人综合+亚洲专区| 91九色精品人成在线观看| 制服人妻中文乱码| 特大巨黑吊av在线直播| or卡值多少钱| 男女视频在线观看网站免费 | 脱女人内裤的视频| 免费在线观看视频国产中文字幕亚洲| 身体一侧抽搐| 亚洲国产看品久久| 久久精品影院6| 母亲3免费完整高清在线观看| 黄色片一级片一级黄色片| 国产视频内射| 精品第一国产精品| 久久伊人香网站| x7x7x7水蜜桃| 成年人黄色毛片网站| 十八禁人妻一区二区| x7x7x7水蜜桃| 真人一进一出gif抽搐免费| 免费看美女性在线毛片视频| a级毛片a级免费在线| 午夜福利免费观看在线| 国产一区二区在线观看日韩 | 一级作爱视频免费观看| 啪啪无遮挡十八禁网站| 一进一出好大好爽视频| ponron亚洲| 亚洲午夜精品一区,二区,三区| 脱女人内裤的视频| 亚洲欧洲精品一区二区精品久久久| 三级毛片av免费| 国产亚洲av嫩草精品影院| 久久草成人影院| 久久久精品大字幕| 香蕉久久夜色| 伊人久久大香线蕉亚洲五| 国产精品久久久久久精品电影| 国产精品香港三级国产av潘金莲| 久久人妻福利社区极品人妻图片| 久99久视频精品免费| 成熟少妇高潮喷水视频| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品中文字幕在线视频| www.999成人在线观看| 国产高清有码在线观看视频 | 亚洲一区高清亚洲精品| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产高清在线一区二区三| 日韩欧美精品v在线| 免费高清视频大片| 在线观看免费日韩欧美大片| 男女那种视频在线观看| av福利片在线观看| 国产不卡一卡二| 亚洲无线在线观看| 在线观看免费日韩欧美大片| 观看免费一级毛片| 日韩大尺度精品在线看网址| 欧美3d第一页| 亚洲av第一区精品v没综合| 中文字幕人妻丝袜一区二区| 欧美黑人巨大hd| 久久这里只有精品19| 亚洲第一欧美日韩一区二区三区| 国产伦一二天堂av在线观看| 久久久久久人人人人人| 国产欧美日韩精品亚洲av| 老鸭窝网址在线观看| 国产精品亚洲av一区麻豆| 可以免费在线观看a视频的电影网站| 最近最新中文字幕大全电影3| 成人国产综合亚洲| 国产av在哪里看| 变态另类丝袜制服| 日韩大码丰满熟妇| 国产精品 欧美亚洲| 91九色精品人成在线观看| 嫩草影视91久久| 久久久久久久久久黄片| 午夜福利视频1000在线观看| xxx96com| 国产午夜福利久久久久久| 国产精品一及| a级毛片在线看网站| 中出人妻视频一区二区| 一区二区三区国产精品乱码| 黄色成人免费大全| 一进一出好大好爽视频| 俺也久久电影网| 校园春色视频在线观看| 国产免费男女视频| 村上凉子中文字幕在线| 黄色视频,在线免费观看| 99在线视频只有这里精品首页| 黄色丝袜av网址大全| 亚洲av五月六月丁香网| 国产精品香港三级国产av潘金莲| 99久久无色码亚洲精品果冻| 午夜福利成人在线免费观看| 国产精品免费视频内射| 久久这里只有精品中国| 18禁裸乳无遮挡免费网站照片| 成熟少妇高潮喷水视频| 成人国产综合亚洲| 成人三级黄色视频| 久久精品影院6| 美女午夜性视频免费| 无遮挡黄片免费观看| 琪琪午夜伦伦电影理论片6080| 后天国语完整版免费观看| 级片在线观看| 一二三四在线观看免费中文在| 99国产极品粉嫩在线观看| АⅤ资源中文在线天堂| 欧美日韩福利视频一区二区| 亚洲av中文字字幕乱码综合| 少妇粗大呻吟视频| 国产一区二区激情短视频| 后天国语完整版免费观看| 天堂影院成人在线观看| 色哟哟哟哟哟哟| 国产私拍福利视频在线观看| 成人手机av| 中文在线观看免费www的网站 | 亚洲av日韩精品久久久久久密| 久久性视频一级片| 男男h啪啪无遮挡| 男人舔奶头视频| 免费高清视频大片| 在线观看免费午夜福利视频| 亚洲中文日韩欧美视频| 欧美+亚洲+日韩+国产| 麻豆国产97在线/欧美 | 成人国产综合亚洲| 两个人免费观看高清视频| 久久香蕉激情| 成年版毛片免费区| 九九热线精品视视频播放| 精品欧美一区二区三区在线| 动漫黄色视频在线观看| 亚洲国产欧洲综合997久久,| 高清毛片免费观看视频网站| 欧美乱色亚洲激情| 国产蜜桃级精品一区二区三区| 听说在线观看完整版免费高清| 精品一区二区三区四区五区乱码| 一级毛片精品| 日本黄大片高清| 欧美av亚洲av综合av国产av| 久久中文字幕人妻熟女| 制服诱惑二区| 很黄的视频免费| 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看 | 久久香蕉精品热| 又粗又爽又猛毛片免费看| 老熟妇仑乱视频hdxx| 五月伊人婷婷丁香| 亚洲激情在线av| 日本精品一区二区三区蜜桃| 久久中文字幕人妻熟女| 欧美成人一区二区免费高清观看 | 久久精品国产亚洲av香蕉五月| 亚洲国产精品合色在线| 亚洲人成网站在线播放欧美日韩| 午夜激情av网站| 亚洲人与动物交配视频| 真人一进一出gif抽搐免费| 狂野欧美激情性xxxx| 法律面前人人平等表现在哪些方面| 人妻久久中文字幕网| 最近最新中文字幕大全电影3| 搞女人的毛片| 亚洲欧美激情综合另类| 床上黄色一级片| 美女高潮喷水抽搐中文字幕| 俺也久久电影网| 日本三级黄在线观看| 亚洲一区中文字幕在线| 好男人电影高清在线观看| 久久久久久国产a免费观看| 真人一进一出gif抽搐免费| 在线十欧美十亚洲十日本专区| 长腿黑丝高跟| 国产aⅴ精品一区二区三区波| 99国产极品粉嫩在线观看| 18美女黄网站色大片免费观看| 波多野结衣巨乳人妻| 国产精品久久久人人做人人爽| 三级毛片av免费| 在线观看免费视频日本深夜| 日本免费a在线| 亚洲在线自拍视频| 欧美黑人精品巨大| 亚洲成a人片在线一区二区| 日本成人三级电影网站| 国产高清有码在线观看视频 | 日韩精品中文字幕看吧| 久久中文看片网| 国产69精品久久久久777片 | 午夜福利免费观看在线| 日韩大尺度精品在线看网址| 色综合亚洲欧美另类图片| 精品无人区乱码1区二区| 亚洲专区字幕在线| 国产精品久久久av美女十八| 亚洲成人久久爱视频| 成人一区二区视频在线观看| 国产三级中文精品| 天堂√8在线中文| 女人被狂操c到高潮| 亚洲电影在线观看av| 岛国在线免费视频观看| 香蕉久久夜色| 天堂动漫精品| 高清在线国产一区| 国产一区在线观看成人免费| 亚洲成a人片在线一区二区| 在线观看www视频免费| 51午夜福利影视在线观看| 午夜精品一区二区三区免费看| 国内揄拍国产精品人妻在线| 日本撒尿小便嘘嘘汇集6| 日本黄大片高清| 国产伦人伦偷精品视频| 欧美高清成人免费视频www| 五月玫瑰六月丁香| 高清在线国产一区| 国产在线观看jvid| 精品一区二区三区四区五区乱码| 日本免费一区二区三区高清不卡| 欧美成狂野欧美在线观看| 日韩精品中文字幕看吧| 国产亚洲精品久久久久5区| 精品国内亚洲2022精品成人| 亚洲国产欧洲综合997久久,| 制服人妻中文乱码| 他把我摸到了高潮在线观看| 色综合欧美亚洲国产小说| 国产高清视频在线播放一区| 亚洲欧美日韩高清专用| 精品少妇一区二区三区视频日本电影| cao死你这个sao货| 国产伦在线观看视频一区| 哪里可以看免费的av片| 亚洲第一电影网av| 两个人免费观看高清视频| 欧美黄色淫秽网站| 男女午夜视频在线观看| 非洲黑人性xxxx精品又粗又长| 国产亚洲av高清不卡| 欧美zozozo另类| 久久久久久久久中文| 免费高清视频大片| 亚洲人成网站高清观看| 香蕉av资源在线| 女人被狂操c到高潮| 国产午夜福利久久久久久| 琪琪午夜伦伦电影理论片6080| 午夜福利高清视频| 日韩欧美 国产精品| 国产亚洲欧美在线一区二区| 亚洲av成人av| 日韩av在线大香蕉| 国产人伦9x9x在线观看| 国产精品国产高清国产av| 国产av一区二区精品久久| 一本一本综合久久| 91在线观看av| 不卡av一区二区三区| 国产1区2区3区精品| 动漫黄色视频在线观看| 99精品欧美一区二区三区四区| 18禁裸乳无遮挡免费网站照片| 中文在线观看免费www的网站 | 久久天躁狠狠躁夜夜2o2o| 91在线观看av| 九九热线精品视视频播放| 啦啦啦韩国在线观看视频| 18禁国产床啪视频网站| 欧美久久黑人一区二区| 日韩欧美三级三区| 久久人妻av系列| 久久精品aⅴ一区二区三区四区| 中出人妻视频一区二区| 亚洲国产精品sss在线观看| 久久这里只有精品19| 国产高清视频在线播放一区| 亚洲成人久久爱视频| 国产一区二区三区视频了| 黄片大片在线免费观看| 又紧又爽又黄一区二区| 午夜福利在线观看吧| 亚洲va日本ⅴa欧美va伊人久久| 香蕉丝袜av| 麻豆一二三区av精品| 天堂√8在线中文| 久久热在线av| 国产精品免费一区二区三区在线| 久久久久久久久免费视频了| 日本精品一区二区三区蜜桃| 免费在线观看日本一区| 亚洲精品美女久久久久99蜜臀| 精品一区二区三区av网在线观看| 美女黄网站色视频| 欧美一级a爱片免费观看看 | 欧美日本视频| 亚洲狠狠婷婷综合久久图片| 9191精品国产免费久久| 国产精品av久久久久免费| 看片在线看免费视频| 长腿黑丝高跟| 丁香六月欧美| 19禁男女啪啪无遮挡网站| 精品第一国产精品| 亚洲成人国产一区在线观看| 欧美一级a爱片免费观看看 | 亚洲国产精品合色在线| 日韩av在线大香蕉| 超碰成人久久| 18禁国产床啪视频网站| 国产69精品久久久久777片 | 国产精品久久久久久精品电影| 精品午夜福利视频在线观看一区| 亚洲美女视频黄频| 黄色丝袜av网址大全| 欧美黑人精品巨大| 久久午夜综合久久蜜桃| 深夜精品福利| 看免费av毛片| 久久久水蜜桃国产精品网| 精品国产乱子伦一区二区三区| 一二三四社区在线视频社区8| 欧美久久黑人一区二区| 床上黄色一级片| 成人一区二区视频在线观看| a级毛片a级免费在线| 精品高清国产在线一区| 精品国产乱码久久久久久男人| 亚洲av片天天在线观看| 免费看美女性在线毛片视频| 亚洲自偷自拍图片 自拍| 亚洲av日韩精品久久久久久密| 精品国内亚洲2022精品成人| 天堂√8在线中文| 啦啦啦观看免费观看视频高清| www.999成人在线观看| 亚洲 欧美一区二区三区| 哪里可以看免费的av片| 亚洲午夜理论影院| 叶爱在线成人免费视频播放| 亚洲一区中文字幕在线| 麻豆av在线久日| 午夜两性在线视频| 男插女下体视频免费在线播放| 999久久久精品免费观看国产| 两性夫妻黄色片| 日本黄色视频三级网站网址| 久9热在线精品视频| 国产高清视频在线播放一区| 亚洲成a人片在线一区二区| 欧美日本亚洲视频在线播放| 婷婷丁香在线五月| 国产av一区二区精品久久| 制服人妻中文乱码| 天堂√8在线中文| 亚洲中文日韩欧美视频| 露出奶头的视频| 欧美日本视频| 亚洲aⅴ乱码一区二区在线播放 | 国产精品美女特级片免费视频播放器 | 老司机深夜福利视频在线观看| 香蕉丝袜av| 天天一区二区日本电影三级| 久久久久国内视频| 国产午夜福利久久久久久| 成人av在线播放网站| 男人舔女人的私密视频| 欧美人与性动交α欧美精品济南到| 激情在线观看视频在线高清| 动漫黄色视频在线观看| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看 | 久久久国产成人免费| 成人国产综合亚洲| 成在线人永久免费视频| 香蕉久久夜色| 亚洲黑人精品在线| 999久久久精品免费观看国产| 国产熟女午夜一区二区三区| videosex国产| 在线观看美女被高潮喷水网站 | 在线观看美女被高潮喷水网站 | 亚洲熟女毛片儿| 午夜精品在线福利| 午夜福利在线在线| 欧美中文日本在线观看视频| 色综合亚洲欧美另类图片| 黄色a级毛片大全视频| 久久香蕉国产精品| 欧美 亚洲 国产 日韩一| 免费人成视频x8x8入口观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲五月婷婷丁香| 12—13女人毛片做爰片一| 免费在线观看成人毛片| www日本黄色视频网| 美女扒开内裤让男人捅视频| 精品乱码久久久久久99久播| 国产激情偷乱视频一区二区| 国产伦人伦偷精品视频| 一边摸一边做爽爽视频免费| 在线观看日韩欧美| 久久人人精品亚洲av| 午夜福利成人在线免费观看| 亚洲,欧美精品.| 国产激情欧美一区二区| 听说在线观看完整版免费高清| 日韩精品中文字幕看吧| 久久久久九九精品影院| 此物有八面人人有两片| 成人三级做爰电影| 免费在线观看完整版高清| 国产成人av激情在线播放| videosex国产| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久久末码| 久久精品91蜜桃| 精品熟女少妇八av免费久了| 午夜福利在线在线| 久久伊人香网站| 国产av不卡久久| 9191精品国产免费久久| 亚洲国产欧美一区二区综合| 精品一区二区三区视频在线观看免费| 97超级碰碰碰精品色视频在线观看| 国产一区在线观看成人免费| xxx96com| 欧美日韩亚洲综合一区二区三区_| 韩国av一区二区三区四区| 亚洲国产欧美一区二区综合| 日本黄大片高清| 亚洲免费av在线视频| 国产精品久久久av美女十八| 两个人看的免费小视频| 极品教师在线免费播放| 午夜激情av网站| 久久午夜亚洲精品久久| 人人妻人人澡欧美一区二区| 亚洲欧洲精品一区二区精品久久久| www.999成人在线观看| 在线观看舔阴道视频| 久久天躁狠狠躁夜夜2o2o| 悠悠久久av| 国产日本99.免费观看| 在线视频色国产色| 一本精品99久久精品77| 国产精品爽爽va在线观看网站| www日本在线高清视频| 毛片女人毛片| 午夜影院日韩av| 午夜久久久久精精品| 国产黄a三级三级三级人| 无限看片的www在线观看| 18美女黄网站色大片免费观看| 亚洲专区中文字幕在线| 法律面前人人平等表现在哪些方面| 美女高潮喷水抽搐中文字幕| 欧美乱色亚洲激情| 国产亚洲av嫩草精品影院| 久久久久亚洲av毛片大全| 美女扒开内裤让男人捅视频| 哪里可以看免费的av片| 久久久久国产精品人妻aⅴ院| 日韩高清综合在线| 免费看a级黄色片| 丰满人妻熟妇乱又伦精品不卡| 欧美一区二区国产精品久久精品 | 久久这里只有精品中国| 最近最新免费中文字幕在线| 亚洲一码二码三码区别大吗| www国产在线视频色| 国产又色又爽无遮挡免费看| 国产单亲对白刺激| 村上凉子中文字幕在线| 欧美丝袜亚洲另类 | avwww免费| 午夜精品一区二区三区免费看| 久久中文字幕一级| 熟女电影av网| 亚洲人成伊人成综合网2020| 国产私拍福利视频在线观看| 久久精品亚洲精品国产色婷小说| 俺也久久电影网| 久久人人精品亚洲av| 亚洲免费av在线视频| 色在线成人网| 欧美中文日本在线观看视频| ponron亚洲| 美女 人体艺术 gogo| 成人av一区二区三区在线看| 一级作爱视频免费观看| 欧美+亚洲+日韩+国产| 欧美中文综合在线视频| videosex国产| 亚洲精品一卡2卡三卡4卡5卡| 久久久精品国产亚洲av高清涩受| 久久精品aⅴ一区二区三区四区| 窝窝影院91人妻| 母亲3免费完整高清在线观看| 伊人久久大香线蕉亚洲五| 精品久久久久久成人av| 久久香蕉国产精品| 国产黄片美女视频| 老司机午夜十八禁免费视频| 久久精品国产清高在天天线| 亚洲av片天天在线观看| 日本一二三区视频观看| 亚洲真实伦在线观看| 亚洲天堂国产精品一区在线| 黄色视频,在线免费观看| 国产视频一区二区在线看| 欧美丝袜亚洲另类 | 午夜精品在线福利| 巨乳人妻的诱惑在线观看| 亚洲精品av麻豆狂野| 国产精品久久久久久人妻精品电影| 欧美+亚洲+日韩+国产| 午夜久久久久精精品| 国产伦人伦偷精品视频| 久久久久久人人人人人| 色尼玛亚洲综合影院| 欧美性猛交╳xxx乱大交人| 中文亚洲av片在线观看爽| 可以免费在线观看a视频的电影网站| 一二三四社区在线视频社区8| 又黄又爽又免费观看的视频| 非洲黑人性xxxx精品又粗又长| 久久天堂一区二区三区四区| 免费在线观看视频国产中文字幕亚洲| 精品熟女少妇八av免费久了| 久久婷婷人人爽人人干人人爱| 久久久久免费精品人妻一区二区| 一本一本综合久久| 一a级毛片在线观看| 亚洲国产欧美网| 久久九九热精品免费| 啦啦啦免费观看视频1| 中文字幕av在线有码专区| 午夜免费激情av| 少妇人妻一区二区三区视频| 日日爽夜夜爽网站| 国产精品 国内视频| 嫩草影视91久久| 欧美大码av| 免费看美女性在线毛片视频| 一进一出抽搐gif免费好疼| 精品久久久久久久久久久久久| 久久中文字幕一级| 国产精品综合久久久久久久免费| 国产久久久一区二区三区| 亚洲性夜色夜夜综合| 999久久久精品免费观看国产| 午夜福利成人在线免费观看| 后天国语完整版免费观看| 国内揄拍国产精品人妻在线|