• <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| 国产黄a三级三级三级人| 最近最新中文字幕免费大全7| 亚洲人成网站高清观看| 欧美一级a爱片免费观看看| 日日啪夜夜爽| 一级爰片在线观看| 少妇的逼水好多| 国产伦理片在线播放av一区| h日本视频在线播放| 成人毛片60女人毛片免费| 免费高清在线观看视频在线观看| 亚洲欧洲国产日韩| 国产欧美另类精品又又久久亚洲欧美| 噜噜噜噜噜久久久久久91| 少妇裸体淫交视频免费看高清| 男女视频在线观看网站免费| 国产老妇伦熟女老妇高清| 中文字幕人妻熟人妻熟丝袜美| 欧美不卡视频在线免费观看| 亚洲熟女精品中文字幕| 毛片一级片免费看久久久久| 十八禁国产超污无遮挡网站| 日韩欧美 国产精品| 精品国产三级普通话版| 亚洲国产最新在线播放| 午夜精品一区二区三区免费看| 精品不卡国产一区二区三区| 九九久久精品国产亚洲av麻豆| av在线播放精品| 亚洲精品中文字幕在线视频 | 久久国产乱子免费精品| 亚洲国产精品成人综合色| 亚洲精品中文字幕在线视频 | 一级毛片电影观看| 久久久久性生活片| 97超碰精品成人国产| 男女边摸边吃奶| 亚洲国产欧美人成| 99热全是精品| 国产成人freesex在线| 日本与韩国留学比较| 亚洲最大成人av| 亚洲精品第二区| 波多野结衣巨乳人妻| 日韩人妻高清精品专区| 国产真实伦视频高清在线观看| 五月玫瑰六月丁香| 午夜激情欧美在线| 日本熟妇午夜| 精品一区二区免费观看| 亚洲最大成人手机在线| 亚洲av中文字字幕乱码综合| 国产精品福利在线免费观看| 日韩成人av中文字幕在线观看| 美女高潮的动态| 亚洲无线观看免费| 少妇丰满av| 性插视频无遮挡在线免费观看| 免费黄色在线免费观看| 天堂影院成人在线观看| 午夜爱爱视频在线播放| 麻豆乱淫一区二区| 在现免费观看毛片| 极品教师在线视频| 日日啪夜夜爽| 97超碰精品成人国产| 夫妻性生交免费视频一级片| 国产精品综合久久久久久久免费| 3wmmmm亚洲av在线观看| 午夜激情久久久久久久| 国产黄片视频在线免费观看| 男人和女人高潮做爰伦理| 久久久久久久大尺度免费视频| 中文字幕久久专区| 99久久九九国产精品国产免费| 国产免费一级a男人的天堂| eeuss影院久久| 亚洲av成人av| 丝袜美腿在线中文| 秋霞在线观看毛片| 亚洲人成网站高清观看| 精品少妇黑人巨大在线播放| 99re6热这里在线精品视频| 成人鲁丝片一二三区免费| 最近最新中文字幕大全电影3| 一级爰片在线观看| 人人妻人人看人人澡| 久久久久久久久久人人人人人人| 久久热精品热| 3wmmmm亚洲av在线观看| 色综合站精品国产| 性色avwww在线观看| 国产一级毛片在线| 亚洲在线自拍视频| 亚洲av二区三区四区| 日产精品乱码卡一卡2卡三| 国产精品国产三级国产av玫瑰| 国产欧美另类精品又又久久亚洲欧美| 精品久久久久久电影网| 综合色av麻豆| 男女视频在线观看网站免费| 99热这里只有是精品在线观看| 精品一区二区三卡| 免费av不卡在线播放| 亚洲欧美一区二区三区黑人 | 久久久久久久久久黄片| 久久综合国产亚洲精品| 美女cb高潮喷水在线观看| 亚洲精品456在线播放app| 日韩电影二区| 精品国产一区二区三区久久久樱花 | 1000部很黄的大片| 日日撸夜夜添| 国产av国产精品国产| 亚洲成人一二三区av| 岛国毛片在线播放| 2018国产大陆天天弄谢| 国产精品嫩草影院av在线观看| 亚洲人成网站在线播| 啦啦啦韩国在线观看视频| 成人亚洲欧美一区二区av| 亚洲欧美日韩无卡精品| 偷拍熟女少妇极品色| 建设人人有责人人尽责人人享有的 | 国产有黄有色有爽视频| 亚洲国产av新网站| 天天躁日日操中文字幕| 国产久久久一区二区三区| 亚洲,欧美,日韩| av在线观看视频网站免费| 男女国产视频网站| 久久久久精品性色| 国产日韩欧美在线精品| 一个人看视频在线观看www免费| 人妻夜夜爽99麻豆av| 最新中文字幕久久久久| 深爱激情五月婷婷| 最近的中文字幕免费完整| 国产一区有黄有色的免费视频 | 色综合亚洲欧美另类图片| 国产高清有码在线观看视频| 午夜爱爱视频在线播放| 十八禁国产超污无遮挡网站| 黄色日韩在线| 精品久久久久久久久亚洲| 18禁在线播放成人免费| 成年免费大片在线观看| 成人午夜精彩视频在线观看| 精品一区二区三区人妻视频| 免费观看无遮挡的男女| 久久久精品免费免费高清| 亚洲自拍偷在线| 毛片一级片免费看久久久久| 国产亚洲精品av在线| 日韩伦理黄色片| 午夜福利视频1000在线观看| 日韩成人伦理影院| 亚洲av免费在线观看| 少妇人妻精品综合一区二区| 久久久久久久久久久免费av| 国产一区二区亚洲精品在线观看| 精品99又大又爽又粗少妇毛片| 亚洲国产高清在线一区二区三| 丰满乱子伦码专区| 99热这里只有精品一区| 亚洲精品日韩在线中文字幕| 高清欧美精品videossex| 免费av观看视频| 亚洲欧洲国产日韩| 青春草国产在线视频| 亚洲经典国产精华液单| 深爱激情五月婷婷| 联通29元200g的流量卡| 搡女人真爽免费视频火全软件| 禁无遮挡网站| 亚洲欧美日韩东京热| 亚洲精品乱久久久久久| 亚洲国产成人一精品久久久| 欧美精品一区二区大全| 最近中文字幕高清免费大全6| 国产高清三级在线| 国产午夜精品久久久久久一区二区三区| 亚洲国产最新在线播放| 十八禁网站网址无遮挡 | 日韩在线高清观看一区二区三区| 亚洲精品乱久久久久久| 韩国高清视频一区二区三区| 麻豆久久精品国产亚洲av| 国产综合精华液| 国产色爽女视频免费观看| 成人毛片a级毛片在线播放| 天美传媒精品一区二区| 秋霞伦理黄片| 国产成人91sexporn| 非洲黑人性xxxx精品又粗又长| 伦理电影大哥的女人| 国产国拍精品亚洲av在线观看| 欧美 日韩 精品 国产| 国内精品美女久久久久久| 日本免费a在线| 国产综合懂色| 欧美日韩精品成人综合77777| 2021少妇久久久久久久久久久| 特级一级黄色大片| 91久久精品国产一区二区三区| 别揉我奶头 嗯啊视频| 亚洲精品自拍成人| 欧美成人一区二区免费高清观看| 特级一级黄色大片| 一级毛片电影观看| 国产成人精品一,二区| 深夜a级毛片| 在线 av 中文字幕| 国产男女超爽视频在线观看| 嫩草影院新地址| 国产淫语在线视频| 一级毛片 在线播放| 国产乱来视频区| 国产熟女欧美一区二区| 69人妻影院| 欧美潮喷喷水| 亚洲伊人久久精品综合| 久久久久久久大尺度免费视频| 男的添女的下面高潮视频| 日韩av在线大香蕉| 亚洲av中文字字幕乱码综合| xxx大片免费视频| 中文字幕亚洲精品专区| 国产成人精品久久久久久| 欧美性猛交╳xxx乱大交人| 国产精品麻豆人妻色哟哟久久 | 久久人人爽人人片av| 国产伦在线观看视频一区| 亚洲精品乱码久久久久久按摩| 久久久国产一区二区| 国产在视频线精品| 亚洲在线观看片| 永久免费av网站大全| 淫秽高清视频在线观看| 久久99热6这里只有精品| 99re6热这里在线精品视频| 少妇被粗大猛烈的视频| 大片免费播放器 马上看| 亚洲天堂国产精品一区在线| 成年女人看的毛片在线观看| 精品午夜福利在线看| 人妻少妇偷人精品九色| 日韩成人伦理影院| 最新中文字幕久久久久| 18禁裸乳无遮挡免费网站照片| 国产精品麻豆人妻色哟哟久久 | 99热这里只有是精品50| 人人妻人人澡欧美一区二区| 国产一区亚洲一区在线观看| 亚洲精品日韩av片在线观看| 欧美另类一区| 久99久视频精品免费| 80岁老熟妇乱子伦牲交| 国产永久视频网站| 国产高清国产精品国产三级 | 美女主播在线视频| 在线观看免费高清a一片| 乱码一卡2卡4卡精品| 纵有疾风起免费观看全集完整版 | 韩国高清视频一区二区三区| 免费大片黄手机在线观看| 久久99热这里只有精品18| 春色校园在线视频观看| 国产伦精品一区二区三区四那| 国产精品福利在线免费观看| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 麻豆国产97在线/欧美| 精品久久久久久久久久久久久| 久久国产乱子免费精品| 干丝袜人妻中文字幕| 精品人妻偷拍中文字幕| 久久久久网色| 国产精品一区二区性色av| 国产精品久久视频播放| 三级经典国产精品| 日日摸夜夜添夜夜爱| 成人美女网站在线观看视频| 神马国产精品三级电影在线观看| 国产高清国产精品国产三级 | 国产精品无大码| 日本av手机在线免费观看| 亚洲久久久久久中文字幕| 全区人妻精品视频| 色综合站精品国产| 女的被弄到高潮叫床怎么办| 老司机影院毛片| 国产国拍精品亚洲av在线观看| 非洲黑人性xxxx精品又粗又长| 人妻制服诱惑在线中文字幕| 亚洲伊人久久精品综合| 国产人妻一区二区三区在| 不卡视频在线观看欧美| 99久久中文字幕三级久久日本| 简卡轻食公司| 国产毛片a区久久久久| 欧美性感艳星| 18+在线观看网站| 精品久久久久久久久av| 欧美潮喷喷水| 国国产精品蜜臀av免费| 成人美女网站在线观看视频| 五月玫瑰六月丁香| 成人综合一区亚洲| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美 国产精品| 日本一二三区视频观看| 国产午夜福利久久久久久| 欧美成人一区二区免费高清观看| 国产高清不卡午夜福利| 日韩欧美精品免费久久| 不卡视频在线观看欧美| 亚洲av免费在线观看| 色吧在线观看| 国产一区亚洲一区在线观看| 男人爽女人下面视频在线观看| 国产精品麻豆人妻色哟哟久久 | 日本猛色少妇xxxxx猛交久久| 最近2019中文字幕mv第一页| 亚洲最大成人中文| 亚州av有码| 成人鲁丝片一二三区免费| 九九在线视频观看精品| a级毛片免费高清观看在线播放| 蜜桃久久精品国产亚洲av| 两个人的视频大全免费| 男女国产视频网站| 最近2019中文字幕mv第一页| 亚洲性久久影院| 国产成人91sexporn| 国产精品嫩草影院av在线观看| 尾随美女入室| 国产精品久久久久久精品电影| 狂野欧美白嫩少妇大欣赏| 99久久九九国产精品国产免费| 精品久久久久久久久av| 99久国产av精品| 国产伦一二天堂av在线观看| av免费观看日本| 欧美xxxx黑人xx丫x性爽| 国产色婷婷99| 免费高清在线观看视频在线观看| 午夜福利在线观看吧| 欧美变态另类bdsm刘玥| 男人舔女人下体高潮全视频| 简卡轻食公司| 国产黄片美女视频| 男女边吃奶边做爰视频| 国产精品久久久久久久久免| 亚洲成色77777| 天堂√8在线中文| 简卡轻食公司| 国产毛片a区久久久久| 一级a做视频免费观看| 中文字幕亚洲精品专区| ponron亚洲| 亚洲国产精品成人综合色| 国产激情偷乱视频一区二区| 亚洲成人中文字幕在线播放| 日日摸夜夜添夜夜爱| 成人一区二区视频在线观看| 婷婷色综合大香蕉| 色播亚洲综合网| 精品人妻一区二区三区麻豆| 国产男女超爽视频在线观看| 午夜激情久久久久久久| 搡老乐熟女国产| 免费黄频网站在线观看国产| 亚洲精品自拍成人| 大片免费播放器 马上看| 午夜老司机福利剧场| 干丝袜人妻中文字幕| 亚洲色图av天堂| 国产精品人妻久久久影院| 99热6这里只有精品| 国产av码专区亚洲av| 搡老妇女老女人老熟妇| 免费少妇av软件| 婷婷六月久久综合丁香| 久久久久性生活片| 精品熟女少妇av免费看| 69人妻影院| 赤兔流量卡办理| 亚洲伊人久久精品综合| 日本黄色片子视频| 永久网站在线| 91在线精品国自产拍蜜月| 国产日韩欧美在线精品| 男女下面进入的视频免费午夜| 简卡轻食公司| 人人妻人人澡人人爽人人夜夜 | 一级毛片我不卡| 欧美+日韩+精品| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 夫妻午夜视频| 伦精品一区二区三区| 一区二区三区四区激情视频| 精品不卡国产一区二区三区| 在线免费十八禁| 在线观看av片永久免费下载| 免费看av在线观看网站| 久久这里只有精品中国| 蜜桃久久精品国产亚洲av| 国产精品av视频在线免费观看| 中文在线观看免费www的网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 联通29元200g的流量卡| 免费电影在线观看免费观看| 午夜福利成人在线免费观看| 国产精品久久久久久久电影| 亚洲av二区三区四区| 色吧在线观看| 极品教师在线视频| 亚洲电影在线观看av| 国产成人a∨麻豆精品| 久久亚洲国产成人精品v| 久久久国产一区二区| 国产大屁股一区二区在线视频| 美女大奶头视频| 国产精品国产三级国产专区5o| 一区二区三区高清视频在线| 2021少妇久久久久久久久久久| 大香蕉久久网| 青青草视频在线视频观看| 国产伦精品一区二区三区四那| 美女黄网站色视频| 99热这里只有精品一区| 亚洲熟女精品中文字幕| 国产乱来视频区| 亚洲精品亚洲一区二区| 日韩欧美精品免费久久| 久久鲁丝午夜福利片| 丝袜美腿在线中文| 欧美激情久久久久久爽电影| 舔av片在线| av线在线观看网站| 婷婷色麻豆天堂久久| 国内精品宾馆在线| 成人亚洲欧美一区二区av| 男插女下体视频免费在线播放| 午夜免费男女啪啪视频观看| 久久久久久久久久黄片| 国产色婷婷99| 国产精品女同一区二区软件| av天堂中文字幕网| 久99久视频精品免费| 亚洲成人精品中文字幕电影| 伊人久久国产一区二区| 2018国产大陆天天弄谢| 熟妇人妻久久中文字幕3abv| 能在线免费看毛片的网站| 美女脱内裤让男人舔精品视频| 亚洲精品影视一区二区三区av| 日本熟妇午夜| 噜噜噜噜噜久久久久久91| 精品一区二区三卡| 国产 亚洲一区二区三区 | 国产免费视频播放在线视频 | 男女边吃奶边做爰视频| 日本一本二区三区精品| 国产精品久久久久久精品电影| 午夜精品国产一区二区电影 | 日本免费在线观看一区| 欧美性猛交╳xxx乱大交人| 国产精品久久视频播放| 国产在线男女| 天堂√8在线中文| 丝袜喷水一区| 国产伦精品一区二区三区视频9| 一级av片app| 亚洲在久久综合| 男的添女的下面高潮视频| 禁无遮挡网站| 99久久九九国产精品国产免费| 亚洲aⅴ乱码一区二区在线播放| 女人久久www免费人成看片| av卡一久久| 美女脱内裤让男人舔精品视频| 午夜福利视频1000在线观看| 午夜福利网站1000一区二区三区| 少妇猛男粗大的猛烈进出视频 | 欧美成人午夜免费资源| 亚洲av成人精品一二三区| 久久99热6这里只有精品| 精品人妻偷拍中文字幕| 国产老妇伦熟女老妇高清| 白带黄色成豆腐渣| 观看美女的网站| 亚洲无线观看免费| 我要看日韩黄色一级片| 欧美变态另类bdsm刘玥| 色播亚洲综合网| 成人毛片60女人毛片免费| 国产单亲对白刺激| av又黄又爽大尺度在线免费看| 色5月婷婷丁香| 能在线免费看毛片的网站| 男女边摸边吃奶| 免费黄网站久久成人精品| 久久久国产一区二区| 看黄色毛片网站| 嫩草影院精品99| 欧美人与善性xxx| 99九九线精品视频在线观看视频| 精品亚洲乱码少妇综合久久| av播播在线观看一区| 99久久精品国产国产毛片| 天天一区二区日本电影三级| 日韩成人av中文字幕在线观看| 成年女人在线观看亚洲视频 | 久久精品夜夜夜夜夜久久蜜豆| 免费黄网站久久成人精品| 美女黄网站色视频| 日韩欧美一区视频在线观看 | 丝袜喷水一区| 精品久久国产蜜桃| 国产伦理片在线播放av一区| 免费观看无遮挡的男女| 99视频精品全部免费 在线| 亚洲国产精品成人久久小说| 亚洲乱码一区二区免费版| 成年av动漫网址| 国产欧美另类精品又又久久亚洲欧美| 国产免费又黄又爽又色| 男女边吃奶边做爰视频| 在线播放无遮挡| 视频中文字幕在线观看| 亚洲欧美日韩东京热| 国产男女超爽视频在线观看| 青青草视频在线视频观看| 欧美激情国产日韩精品一区| 欧美性感艳星| 国产午夜精品论理片| 欧美性感艳星| 少妇熟女aⅴ在线视频| av播播在线观看一区| 最近的中文字幕免费完整| 91av网一区二区| 中文在线观看免费www的网站| 有码 亚洲区| 特级一级黄色大片| 国产精品久久视频播放| 91久久精品电影网| 国产精品一二三区在线看| 亚洲av成人精品一二三区| 国产色婷婷99| 久久综合国产亚洲精品| 欧美最新免费一区二区三区| 亚洲欧美中文字幕日韩二区| 亚洲欧美精品专区久久| 别揉我奶头 嗯啊视频| 九色成人免费人妻av| 国产精品99久久久久久久久| 97热精品久久久久久| 中文字幕av成人在线电影| 晚上一个人看的免费电影| 日韩欧美三级三区| 免费看不卡的av| 国产精品久久久久久av不卡| 国产高潮美女av| 国产69精品久久久久777片| 免费av毛片视频| 国产欧美日韩精品一区二区| 嫩草影院新地址| 精品熟女少妇av免费看| 三级国产精品片| 久久这里只有精品中国| 亚洲丝袜综合中文字幕| 国产高清有码在线观看视频| 欧美xxⅹ黑人| 麻豆av噜噜一区二区三区| 日日撸夜夜添| 免费看不卡的av| 精品久久久久久久末码| 男女啪啪激烈高潮av片| 熟妇人妻不卡中文字幕| 亚洲内射少妇av| 黄色配什么色好看| 亚洲精品日本国产第一区| 亚洲精品日韩av片在线观看| 国产一区有黄有色的免费视频 | 最近的中文字幕免费完整| 男人舔女人下体高潮全视频| 大又大粗又爽又黄少妇毛片口| 搡老妇女老女人老熟妇| 免费观看av网站的网址| 日本黄色片子视频| 99久久中文字幕三级久久日本| 亚洲精华国产精华液的使用体验| 美女内射精品一级片tv| 日韩精品青青久久久久久| av线在线观看网站| 成人特级av手机在线观看| 亚洲av二区三区四区| 欧美xxxx性猛交bbbb| 国产伦精品一区二区三区视频9| 成人午夜精彩视频在线观看| 欧美成人a在线观看| 亚洲精品视频女| 久久精品国产鲁丝片午夜精品| 91久久精品国产一区二区三区| .国产精品久久| 国产中年淑女户外野战色| 国产精品av视频在线免费观看|