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

    煤層氣儲層滲透率的預(yù)測方法

    2016-12-09 07:30:52張欣欣李婷婷車向前
    關(guān)鍵詞:煤層氣滲透率交叉

    邊 莉, 張欣欣, 李婷婷, 車向前

    (1.黑龍江科技大學(xué) 電子與信息工程學(xué)院, 哈爾濱 150022; 2.黑龍江科技大學(xué) 電氣與控制工程學(xué)院,

    ?

    煤層氣儲層滲透率的預(yù)測方法

    邊 莉1, 張欣欣2, 李婷婷2, 車向前3

    (1.黑龍江科技大學(xué) 電子與信息工程學(xué)院, 哈爾濱 150022; 2.黑龍江科技大學(xué) 電氣與控制工程學(xué)院,

    哈爾濱 150022; 3.黑龍江科技大學(xué) 計算機與信息工程學(xué)院,哈爾濱 150022)

    煤層氣儲層滲透率參數(shù)的強非線性,使得常規(guī)參數(shù)難以預(yù)測,由此提出一種交叉熵優(yōu)化支持向量機的煤層氣滲透率預(yù)測新方法。利用UCI數(shù)據(jù)庫中的數(shù)據(jù),將該算法預(yù)測結(jié)果與粒子群算法、遺傳算法的預(yù)測結(jié)果進行對比。結(jié)合某礦區(qū)實際地質(zhì)特性,用該算法對該礦區(qū)煤層氣儲層滲透率參數(shù)進行預(yù)測,得出大部分被測試數(shù)據(jù)的相對誤差率達到 1%左右,除了個別數(shù)據(jù)超過20%,總體分布較穩(wěn)定。該礦區(qū)滲透率參數(shù)的預(yù)測值和真實值分布圖大體一致。該算法可行,且優(yōu)于其他算法,可以為煤層氣儲層參數(shù)的預(yù)測提供算法保障。

    煤層氣; 滲透率預(yù)測; 交叉熵算法; 支持向量機

    煤層氣的儲層參數(shù)能夠描述地下巖性特征,它與煤層氣儲量和產(chǎn)能密切相關(guān),因此,在開發(fā)前,對儲層參數(shù)進行準確預(yù)測,能夠指導(dǎo)資源的勘探與開發(fā),可以有效提升開發(fā)效率。當前,主要煤儲層的滲透率預(yù)測方法包括描述技術(shù)、數(shù)理分析等。描述技術(shù)是根據(jù)煤樣找到與滲透率有關(guān)的特征參數(shù)例如充填程度、割理密度、割理壁距等。數(shù)理分析是指通過一些數(shù)學(xué)方法建立對應(yīng)的數(shù)學(xué)模型,從而預(yù)測滲透率[1-2],這種方法在實際應(yīng)用中并沒有進行可行性分析。目前,越來越多的數(shù)學(xué)方法被應(yīng)用到地質(zhì)研究中,例如統(tǒng)計學(xué)、模式識別和人工神經(jīng)網(wǎng)絡(luò)。但是,利用人工神經(jīng)網(wǎng)絡(luò)建立的模型具有結(jié)構(gòu)不易確定、易出現(xiàn)局部最優(yōu)等缺陷[3]。鑒于此,筆者通過將交叉熵與支持向量機相結(jié)合的方法來對煤層氣儲層滲透率進行預(yù)測,并利用UCI數(shù)據(jù)庫中的樣本對該方法進行可行性分析,結(jié)合某實際礦區(qū)的地質(zhì)特性,對實際煤層氣儲層滲透率進行預(yù)測。

    1 支持向量機

    支持向量機(Support vector machine,SVM)形式簡單,尤其在小數(shù)據(jù)問題上具有較強的優(yōu)越性[4],越來越受到研究者的青睞。SVM基本思想如式(1)所示,將非線性問題φ(·)中輸入數(shù)據(jù)T(T∈n)映射到Hilbert空間中,即將非線性問題線性化。在Hilber中,構(gòu)造最佳分類超平面F(φ(T))的基礎(chǔ)上,建立分類決策函數(shù)。其中,最佳分類超平面的參數(shù)是依據(jù)結(jié)構(gòu)風險最小化原則,利用拉格朗日乘子,構(gòu)造與原始分類問題對偶的非線性規(guī)劃問題并求解,得到[5]

    (1)

    當非線性映射滿足Mercer條件時,可采用核函數(shù)來代替輸入向量與支持向量的內(nèi)積運算。核函數(shù)的選取決定最終SVM的分類效果,同時考慮后續(xù)算法的復(fù)雜度,采用徑向基核函數(shù),為

    (2)

    式中:‖T-Y‖2——向量間距;

    g——徑向基寬度。

    如何恰當?shù)倪x擇參數(shù)是算法的關(guān)鍵,這將影響構(gòu)建模型的精確度,為了協(xié)調(diào)SVM的參數(shù),選取交叉熵支持向量機(Cross entropy SVM, CE-SVM)的方法來優(yōu)化SVM中的兩個重要參數(shù),即懲罰因子c與徑向基寬度g[6]。由于,非線性優(yōu)化問題的復(fù)雜度由g決定,所以g的取值不當,會直接影響SVM的適應(yīng)性;c的取值不當,會引起過學(xué)習或欠學(xué)習。

    2 交叉熵算法的原理

    交叉熵(Cross entropy,CE)算法尋優(yōu)迭代的過程主要為[7]:(1) 根據(jù)具體的概率密度函數(shù)產(chǎn)生隨機樣本;(2) 為產(chǎn)生更優(yōu)的樣本,對概率密度函數(shù)進行更新??紤]式(3)所示典型的無約束連續(xù)型優(yōu)化問題。

    G(x*)=γ*=minG(x),

    (3)

    式中:x——有限集;

    G——實值函數(shù),其中γ*是所求問題的最小值,x*是其最優(yōu)解。

    設(shè)在χ上的一個概率密度函數(shù)族為{f(·;v),v∈V}。 對于給定的概率密度f(·;u),u∈V,優(yōu)化問題又可轉(zhuǎn)化為

    (4)

    式中:{I{G(X)≤γ}}——指示函數(shù)的集合;

    X——由f(·;u)產(chǎn)生的隨機樣本;

    Eu——相應(yīng)的期望值。

    當γ逐漸接近γ*時,G(X)≤γ為小概率事件;所以,選取合適的u確保Pu(G(X)≤γ)的值不會太小[8],顯然γ和u應(yīng)恰當?shù)倪x擇。當所求問題的概率很小時,可以通過采用多級別算法,構(gòu)造參數(shù)序列{vt,t≥0}和級別序列{γt,t≥1},然后對vt和γt更新迭代。文中利用CE算法對c和g迭代尋優(yōu),概率密度函數(shù)選取高斯分布N(μ,σ2),SVM交叉驗證概率作為適應(yīng)度函數(shù)[9]。

    3 交叉熵算法優(yōu)化SVM算法步驟

    采用連續(xù)型交叉熵算法,以SVM(c,g)為優(yōu)化目標,SVM的交叉驗證概率為適應(yīng)度函數(shù),具體優(yōu)化步驟如下[10]。

    (5)

    (6)

    (7)

    步驟5 平滑:對于m=1,2,…,M計算

    μL(t)=βμL(t)+(1-β)μL(t-1),

    (8)

    (9)

    4 交叉熵支持向量機預(yù)測模型驗證

    利用UCI中的混凝土抗壓強度數(shù)據(jù)進行驗證,數(shù)據(jù)包含1 030個對象,9個連續(xù)屬性,同煤層氣儲層參數(shù)預(yù)測中,測井數(shù)據(jù)和地震屬性數(shù)據(jù)構(gòu)建的樣本集較類似。在訓(xùn)練過程中,輸出量為抗壓強度,輸入向量為各個屬性。經(jīng)過歸一化處理之后,在數(shù)據(jù)集中隨機選取300個樣本構(gòu)成訓(xùn)練樣本集,測試集為包含歸一化處理后的全部數(shù)據(jù),對混凝土抗壓強度進行回歸分析,評價預(yù)測性能。

    選取軟件Matlab對此問題編程,實現(xiàn)對SVM核函數(shù)的參數(shù)g與c的迭代尋優(yōu),得出懲罰因子c=13.512 4,徑向基寬度g=0.061 760 3,CV=0.888 889,其c與g的變化曲線如圖1a所示。圖1b為c=13.512 4,g=0.061 760 3時的適應(yīng)度變化曲線即交叉驗證概率曲線。從圖1中可以看出c、g收斂曲線和交叉驗證概率曲線均在進化代數(shù)為40次左右趨于穩(wěn)定,該算法具有較快的收斂速度和穩(wěn)定性。圖2a是混凝土抗壓強度的真實值與交叉熵優(yōu)化支持向量機預(yù)測的情況,只有在樣本值為400左右的位置出現(xiàn)了差異,其余樣本預(yù)測效果較好,圖2b是采用粒子群算法PS和遺傳算法GA對混凝土的抗壓強度進行預(yù)測的結(jié)果,和真實結(jié)果對比可以看出,這兩種智能算法與真實值的結(jié)果在400個樣本后,預(yù)測結(jié)果有很大差別。

    a SVM(c,g)

    b 交叉熵算法

    Fig. 1 Convergence curve of SVM and cross entropy algorithm

    a 真實結(jié)果與CE-SVM預(yù)測結(jié)果

    b PS與GA優(yōu)化SVM結(jié)果的對比

    Fig. 2 Comparison of simulation results of cross entropy algorithm and other intelligent algorithms

    表1對訓(xùn)練樣本為300從均方誤差、相關(guān)系數(shù)和訓(xùn)練時間三個方面對三種算法進行了對比。由表1可看出,交叉熵算法的均方誤差要小于粒子群算法和遺傳算法兩種智能算法,相關(guān)系數(shù)要大于粒子群算法和遺傳算法的兩種智能算法,訓(xùn)練時間要比其他兩種智能算法短。

    綜上所述交叉熵算法優(yōu)化的支持向量機確實有良好的預(yù)測能力,在某些性能上,要優(yōu)于其他智能算法。因此,以交叉熵算法作為核函數(shù)優(yōu)化算法,構(gòu)建支持向量機預(yù)測模型是有效可行的。

    表1 訓(xùn)練集樣本為300的性能評價

    Table 1 Performance evaluation of training set samples for 300

    算法均方誤差相關(guān)系數(shù)訓(xùn)練時間/sCE0.015360.690648.406PS0.017580.646482.148GA0.019900.6020114.016

    5 煤層氣滲透率的預(yù)測

    5.1 數(shù)據(jù)的選取及處理

    測井數(shù)據(jù)主要來源于我國山西沁水盆地南部山西組3號煤層,這些數(shù)據(jù)由兩組構(gòu)成:一組是訓(xùn)練樣本集,用來構(gòu)建測試模型;另一組是測試樣本集,用來對模型進行檢測。在模型構(gòu)建前,需要對樣本數(shù)據(jù)進行歸一化處理,主要有以下兩個原因:一是由于采集到的煤儲層滲透率受數(shù)據(jù)的性質(zhì)、采集技術(shù)等眾多因素影響,這些會對支持向量機的建模和預(yù)測的結(jié)果有一定的影響[11]。二是樣本中的數(shù)據(jù)不易變化太大,它不僅會影響模型的建立,而且對支持向量機的識別和收斂會產(chǎn)生一定的影響。為了盡量減少預(yù)測模型的訓(xùn)練時間、防止數(shù)值變化并且使參數(shù)在計算過程中更容易,所以一般在預(yù)測模型構(gòu)建之前,對樣本集統(tǒng)一進行歸一化處理。滲透率部分樣本的數(shù)據(jù),共有100組,表2選取其中的10組數(shù)據(jù),并對其歸一化處理[12]。

    表2 歸一化后的部分樣本

    5.2 煤層氣滲透率預(yù)測結(jié)果

    圖3為適應(yīng)度曲線,從圖3中可以看出當進化代數(shù)為4的時曲線變得很穩(wěn)定,說明了該算法的收斂速度快的優(yōu)點,圖4是原始滲透率和預(yù)測后的對比,從圖中可以看出大體預(yù)測趨勢較好,圖5是相對誤差圖,盡管有幾個點的誤差比較大,但大部分相對誤差都趨近于0,這更說明了交叉熵支持向量機預(yù)測煤層氣滲透率的可靠性。

    由圖5可知,新算法預(yù)測的結(jié)果和滲透率的真實結(jié)果幾乎一致,并且預(yù)測值與實際值的總體趨勢吻合,個別預(yù)測異常點可以忽略。從中得到的均方根誤差為0.032 692 1%,相關(guān)系數(shù)為50.777 71%。

    由圖6可以看到,其相關(guān)系數(shù)相對來說較大,但是均方根誤差較小,由此可以說明支持向量機的預(yù)測效果較好,尤其對于小樣本的訓(xùn)練數(shù)據(jù)模型更為適用。當相對誤差率的分布范圍在 0~45%時,大部分被測試數(shù)據(jù)的相對誤差率相差不大,通常只有 1%左右的波動,除了個別數(shù)據(jù)超過20%,總體分布穩(wěn)定。

    圖3 適應(yīng)度曲線

    圖4 原始數(shù)據(jù)與回歸預(yù)測數(shù)據(jù)的對比

    圖5 相對誤差

    6 結(jié)束語

    交叉熵支持向量機的煤層氣滲透率預(yù)測是一種新方法。將交叉熵支持向量機結(jié)合起來應(yīng)用到滲透

    率預(yù)測中,通過對該算法的驗證、建模和仿真,并與其他智能算法比較分析,得出該算法對滲透率預(yù)測具有更準確的精度,可以證明該算法具有穩(wěn)定性高、收斂速度快、運行時間短等優(yōu)勢。

    將這兩種算法結(jié)合起來應(yīng)用到煤儲層滲透率中可行。今后可以將這種方法應(yīng)用到煤層氣儲層其他參數(shù)的預(yù)測中,可以為煤層氣的勘探與開發(fā)提供很好的幫助。該算法具有較為廣闊的應(yīng)用前景。

    [1] 蘇現(xiàn)波, 馮艷麗, 陳江峰. 煤中裂隙的分類[J]. 煤田地質(zhì)與勘探, 2002(4): 21-24.

    [2] LEVINE J R. Model study of the influence of matrix shrinkage on absolute Permeability of coal bed reservoir[J]. Geological Society Publication, 1996, 109(1):197-212.

    [3] 李志明, 張金珠. 地應(yīng)力與油氣勘探開發(fā)[J]. 北京:石油工業(yè)出版社, 2007(1): 16-27.

    [4] 李艷芳, 程建遠, 王 成. 基于支持向量機的地震屬性優(yōu)選及煤層氣預(yù)測[J]. 煤田地質(zhì)與勘探, 2012, 40(6): 75-78.

    [5] 邊 莉, 邊晨源. 采用交叉熵支持向量機和模糊積分的電網(wǎng)故障診斷[J]. 電機與控制學(xué)報, 2016, 20(2): 112-120.

    [6] 于文彬. 基于偏最小二乘及最小二乘支持向量機算法在儲層參數(shù)預(yù)測上的應(yīng)用[D]. 長春: 吉林大學(xué), 2013.

    [7] PIETER TJERK DE BOER, REUVEN Y RUBINSTEIN. A tutorial on the cross-entropy method[J]. Annals of Operations Research. 2005, 134(2): 19-67.

    [8] 邊 莉, 邊晨源. 基于組合型交叉熵算法的多電源配電網(wǎng)故障定位[J]. 電力科學(xué)與技術(shù)學(xué)報, 2014, 30(03): 86-91.

    [9] 孟召平, 田永東, 李國富. 沁水盆地南部煤儲層滲透性與地應(yīng)力之間關(guān)系和控制機理[J]. 自然科學(xué)進展, 2009, 10(9): 1142-1148.

    [10] 傅雪海, 秦 勇, 姜 波, 等. 煤割理壓縮實驗及滲透率數(shù)值模擬[J]. 煤炭學(xué)報, 2001(6): 967-972.

    [11] 束詩雨. 基于集成學(xué)習的支持向量機預(yù)測優(yōu)化算法及其應(yīng)用[D]. 東華大學(xué), 2015.

    [12] 王 雷. 支持向量機預(yù)測煤儲層滲透性[D]. 北京: 中國石油大學(xué), 2014.

    (編輯 晁曉筠 校對 李德根)

    Prediction method of reservoir permeability of coal bed methane

    BianLi1,ZhangXinxin2,LiTingting2,CheXiangqian3

    (1.School of Electronics & Information Engineering, Heilongjiang University of Science & Technology,Harbin 150022, China; 2.School of Electrical & Control Engineering, Heilongjiang University of Science & Technology, Harbin 150022, China; 3.School of Computer & Information Engineering,Heilongjiang University of Science & Technology, Harbin 150022, China)

    This paper presents a novel method tailored for coal bed gas permeability, based on cross entropy optimization support vector machine to address the notorious prediction of conventional parameters due to the stronger nonlinearity associated with the permeability parameters behind coalbed gas reservoir. The algorithm is validated by using the data in UCI database, and then comparing the prediction results with the those of particle swarm algorithm and genetic one, suggesting its feasibility and consequently distinct advantages over other ones. The prediction of the coalbed gas reservoir permeability in a mining area with the actual geological characteristics, using the proposed algorithm, indicates that the algorithm operates within the relative error of about 1% in the data measurement except more than 20% in some individual cases, with a relatively stable distribution. The general consistency between the predicted value of permeability parameter and the distribution map of true value of the mining area points to the conclusion that the algorithm can provide a safe prediction of the parameters behind coal bed gas reservoir.

    coal bed methane; permeability prediction; cross entropy algorithm; support vector machine

    2016-08-23

    國家自然科學(xué)基金項目(51504085);哈爾濱市科技局項目(2015RQQXJ009)

    邊 莉(1978-),女,河北省涿州人,副教授,博士,研究方向:陣列信號處理、人工智能技術(shù),E-mail:branran@163.com。

    10.3969/j.issn.2095-7262.2016.05.003

    TD713.2

    2095-7262(2016)05-0480-05

    A

    猜你喜歡
    煤層氣滲透率交叉
    2019年《中國煤層氣》征訂單
    中國煤層氣(2019年2期)2019-08-27 00:59:38
    煤層氣吸附-解吸機理再認識
    中國煤層氣(2019年2期)2019-08-27 00:59:30
    “六法”巧解分式方程
    中煤階煤層氣井排采階段劃分及滲透率變化
    不同滲透率巖芯孔徑分布與可動流體研究
    SAGD井微壓裂儲層滲透率變化規(guī)律研究
    連一連
    高滲透率風電并網(wǎng)對電力系統(tǒng)失步振蕩的影響
    DFS-C02V煤層氣井高產(chǎn)因素簡析
    中國煤層氣(2015年6期)2015-08-22 03:25:29
    基于Fast-ICA的Wigner-Ville分布交叉項消除方法
    計算機工程(2015年8期)2015-07-03 12:19:54
    春色校园在线视频观看| 飞空精品影院首页| 中文字幕亚洲精品专区| 91精品国产国语对白视频| 女人精品久久久久毛片| 久久精品国产自在天天线| www.av在线官网国产| 久久亚洲国产成人精品v| tube8黄色片| 韩国精品一区二区三区 | 国产av一区二区精品久久| 色94色欧美一区二区| 亚洲av日韩在线播放| 国产成人91sexporn| 人人妻人人添人人爽欧美一区卜| 飞空精品影院首页| 在线精品无人区一区二区三| 69精品国产乱码久久久| 久久 成人 亚洲| 性色av一级| 国产精品人妻久久久久久| 成人毛片60女人毛片免费| 亚洲一级一片aⅴ在线观看| 制服人妻中文乱码| 最近最新中文字幕免费大全7| 国产一区有黄有色的免费视频| 国产精品一二三区在线看| 国产精品嫩草影院av在线观看| 美女国产高潮福利片在线看| 五月天丁香电影| 在线 av 中文字幕| 国产精品熟女久久久久浪| 建设人人有责人人尽责人人享有的| 欧美日韩亚洲高清精品| 男女午夜视频在线观看 | 国产又色又爽无遮挡免| 久久人人97超碰香蕉20202| 国产不卡av网站在线观看| 国产福利在线免费观看视频| 99re6热这里在线精品视频| 久久久久久久久久久免费av| 男女免费视频国产| 全区人妻精品视频| 国产白丝娇喘喷水9色精品| 欧美人与善性xxx| 中国美白少妇内射xxxbb| 午夜免费男女啪啪视频观看| 日本与韩国留学比较| 欧美精品av麻豆av| 国产成人欧美| 国产探花极品一区二区| 又大又黄又爽视频免费| 99九九在线精品视频| 亚洲av在线观看美女高潮| 日韩熟女老妇一区二区性免费视频| 美女主播在线视频| 成人无遮挡网站| 精品亚洲乱码少妇综合久久| 高清黄色对白视频在线免费看| 亚洲一级一片aⅴ在线观看| 一本大道久久a久久精品| 免费黄频网站在线观看国产| 最近中文字幕2019免费版| 国产福利在线免费观看视频| 狂野欧美激情性xxxx在线观看| 亚洲国产色片| 久久av网站| www.av在线官网国产| freevideosex欧美| 久久人人爽人人爽人人片va| 毛片一级片免费看久久久久| 久久精品国产亚洲av天美| 熟女av电影| 国产欧美日韩一区二区三区在线| 国精品久久久久久国模美| 亚洲精品av麻豆狂野| 国产精品女同一区二区软件| 青春草亚洲视频在线观看| 九草在线视频观看| 亚洲国产精品一区二区三区在线| 蜜臀久久99精品久久宅男| 国产麻豆69| 丝瓜视频免费看黄片| 人妻系列 视频| 王馨瑶露胸无遮挡在线观看| 我要看黄色一级片免费的| 菩萨蛮人人尽说江南好唐韦庄| 熟女电影av网| 99精国产麻豆久久婷婷| 伦理电影大哥的女人| 亚洲精品中文字幕在线视频| www日本在线高清视频| 久久精品久久久久久久性| 欧美亚洲日本最大视频资源| 在线观看三级黄色| 国产精品.久久久| 最近的中文字幕免费完整| 国产69精品久久久久777片| 韩国av在线不卡| 国产爽快片一区二区三区| 国产成人精品一,二区| 天天影视国产精品| 女人被躁到高潮嗷嗷叫费观| 一区二区av电影网| 久久久久久久久久久免费av| 肉色欧美久久久久久久蜜桃| 国产精品国产三级专区第一集| 人人妻人人澡人人看| 亚洲成人手机| 亚洲精品久久成人aⅴ小说| 天天躁夜夜躁狠狠久久av| av一本久久久久| 99热全是精品| 精品酒店卫生间| 国产麻豆69| 久久国内精品自在自线图片| 亚洲精华国产精华液的使用体验| 国产熟女午夜一区二区三区| 久久韩国三级中文字幕| 久久这里有精品视频免费| 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| 亚洲精品一区蜜桃| 日韩人妻精品一区2区三区| 久久ye,这里只有精品| 热re99久久精品国产66热6| 亚洲,欧美精品.| 国产熟女午夜一区二区三区| 国产精品熟女久久久久浪| 久久精品aⅴ一区二区三区四区 | 成人毛片a级毛片在线播放| 男女边吃奶边做爰视频| 午夜精品国产一区二区电影| 午夜福利在线观看免费完整高清在| 国产成人免费观看mmmm| 看免费av毛片| 精品国产一区二区久久| 成人亚洲欧美一区二区av| 波野结衣二区三区在线| 亚洲av日韩在线播放| 丝袜人妻中文字幕| 看非洲黑人一级黄片| 色吧在线观看| 欧美精品国产亚洲| 一级a做视频免费观看| 日韩成人av中文字幕在线观看| 亚洲高清免费不卡视频| 99九九在线精品视频| 99热这里只有是精品在线观看| 婷婷色av中文字幕| 欧美精品国产亚洲| 亚洲 欧美一区二区三区| 精品卡一卡二卡四卡免费| 欧美国产精品va在线观看不卡| 国产 精品1| 亚洲国产成人一精品久久久| 国产亚洲精品第一综合不卡 | 久久精品久久久久久噜噜老黄| 亚洲成国产人片在线观看| 久久影院123| 咕卡用的链子| 亚洲精品第二区| 亚洲精品国产av蜜桃| 日韩成人伦理影院| 22中文网久久字幕| 久久久久视频综合| 日韩精品免费视频一区二区三区 | 欧美性感艳星| 在线观看www视频免费| 高清视频免费观看一区二区| 日韩一本色道免费dvd| 久久ye,这里只有精品| 99热6这里只有精品| 国产午夜精品一二区理论片| 亚洲精品乱码久久久久久按摩| 亚洲精品国产色婷婷电影| 久久人人爽人人片av| 亚洲欧美成人精品一区二区| 亚洲精品,欧美精品| 久久人人爽av亚洲精品天堂| 一级黄片播放器| 国产成人午夜福利电影在线观看| 日本与韩国留学比较| 国产激情久久老熟女| 最近手机中文字幕大全| 超碰97精品在线观看| 青春草视频在线免费观看| 老司机影院毛片| 91国产中文字幕| 日韩伦理黄色片| 久久久久人妻精品一区果冻| 欧美日韩国产mv在线观看视频| 性色av一级| a级片在线免费高清观看视频| 免费在线观看黄色视频的| 亚洲av免费高清在线观看| 尾随美女入室| 久久人妻熟女aⅴ| 一区二区av电影网| 黄色视频在线播放观看不卡| 成人亚洲欧美一区二区av| 久久久精品94久久精品| tube8黄色片| 超色免费av| 咕卡用的链子| 欧美日本中文国产一区发布| 国产亚洲最大av| 免费黄网站久久成人精品| 大话2 男鬼变身卡| 91精品伊人久久大香线蕉| 成年女人在线观看亚洲视频| 亚洲国产精品专区欧美| 大香蕉97超碰在线| 2022亚洲国产成人精品| 丝袜喷水一区| 99久久人妻综合| 久久久久国产网址| 亚洲一码二码三码区别大吗| 精品国产一区二区久久| 亚洲精品av麻豆狂野| 麻豆乱淫一区二区| 少妇被粗大的猛进出69影院 | 国产熟女午夜一区二区三区| av国产精品久久久久影院| 99九九在线精品视频| 久久精品久久精品一区二区三区| 侵犯人妻中文字幕一二三四区| 亚洲在久久综合| 色婷婷av一区二区三区视频| 五月天丁香电影| 久久国产精品男人的天堂亚洲 | 精品国产一区二区久久| 欧美性感艳星| 国产日韩欧美视频二区| 国产免费又黄又爽又色| 母亲3免费完整高清在线观看 | 哪个播放器可以免费观看大片| 欧美精品国产亚洲| 亚洲成av片中文字幕在线观看 | 18+在线观看网站| 制服丝袜香蕉在线| 国产男女内射视频| 天堂俺去俺来也www色官网| 看免费成人av毛片| 亚洲国产精品专区欧美| 精品久久久精品久久久| 涩涩av久久男人的天堂| 插逼视频在线观看| 少妇人妻久久综合中文| 免费日韩欧美在线观看| 建设人人有责人人尽责人人享有的| 国产精品不卡视频一区二区| 91国产中文字幕| 最近中文字幕高清免费大全6| 在线亚洲精品国产二区图片欧美| 亚洲精品一二三| 国产国语露脸激情在线看| 久久久久精品性色| 亚洲精品456在线播放app| 亚洲中文av在线| 菩萨蛮人人尽说江南好唐韦庄| 一二三四中文在线观看免费高清| 蜜桃国产av成人99| 在线天堂中文资源库| 久久人人爽人人爽人人片va| 水蜜桃什么品种好| 亚洲国产看品久久| 肉色欧美久久久久久久蜜桃| 最黄视频免费看| 国产精品嫩草影院av在线观看| 久久久久久久大尺度免费视频| 免费在线观看完整版高清| 亚洲熟女精品中文字幕| 欧美日韩视频高清一区二区三区二| 亚洲欧美清纯卡通| 亚洲av男天堂| 人妻人人澡人人爽人人| 中文字幕最新亚洲高清| 考比视频在线观看| 婷婷色麻豆天堂久久| 成人18禁高潮啪啪吃奶动态图| 亚洲av免费高清在线观看| 熟妇人妻不卡中文字幕| 看十八女毛片水多多多| 久久这里只有精品19| 久久精品国产综合久久久 | 纵有疾风起免费观看全集完整版| 久久久久人妻精品一区果冻| 久久狼人影院| 精品午夜福利在线看| 熟妇人妻不卡中文字幕| 永久免费av网站大全| 国产一区二区在线观看av| 日本av免费视频播放| 日本爱情动作片www.在线观看| 韩国高清视频一区二区三区| 久久鲁丝午夜福利片| 成人黄色视频免费在线看| 午夜免费鲁丝| 久久久久久久国产电影| 亚洲成色77777| 桃花免费在线播放| 国产精品久久久久成人av| 免费人成在线观看视频色| 免费人妻精品一区二区三区视频| 国产黄色视频一区二区在线观看| 午夜福利网站1000一区二区三区| 欧美 亚洲 国产 日韩一| 久久精品人人爽人人爽视色| 国产免费视频播放在线视频| 午夜福利,免费看| 精品久久国产蜜桃| 久久精品国产亚洲av天美| 大话2 男鬼变身卡| 中国三级夫妇交换| 国产精品人妻久久久影院| 国产精品久久久久久av不卡| 欧美精品高潮呻吟av久久| 日韩av在线免费看完整版不卡| 啦啦啦中文免费视频观看日本| av在线老鸭窝| 午夜日本视频在线| 观看美女的网站| 熟妇人妻不卡中文字幕| 国产福利在线免费观看视频| 久久99热这里只频精品6学生| 极品少妇高潮喷水抽搐| 男女啪啪激烈高潮av片| 亚洲,一卡二卡三卡| 精品99又大又爽又粗少妇毛片| 熟女电影av网| 大片电影免费在线观看免费| 亚洲av免费高清在线观看| 99久久精品国产国产毛片| 久久这里有精品视频免费| 精品一区二区免费观看| 欧美精品亚洲一区二区| 久久久国产精品麻豆| 两个人看的免费小视频| 国产又爽黄色视频| 中文字幕人妻熟女乱码| 久久久久国产精品人妻一区二区| 久久女婷五月综合色啪小说| 久久久精品94久久精品| 亚洲图色成人| 天天躁夜夜躁狠狠久久av| 9191精品国产免费久久| 高清视频免费观看一区二区| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 免费观看a级毛片全部| 欧美人与性动交α欧美软件 | 久久精品久久久久久噜噜老黄| 黄色视频在线播放观看不卡| 国产探花极品一区二区| 又大又黄又爽视频免费| 国产日韩欧美视频二区| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 国产色爽女视频免费观看| 午夜日本视频在线| 国产精品久久久久久av不卡| 秋霞伦理黄片| 少妇人妻久久综合中文| 在线亚洲精品国产二区图片欧美| 中文天堂在线官网| 国产免费一区二区三区四区乱码| 欧美激情国产日韩精品一区| 日韩三级伦理在线观看| 在线亚洲精品国产二区图片欧美| 欧美亚洲 丝袜 人妻 在线| 我要看黄色一级片免费的| 欧美人与善性xxx| 少妇的逼好多水| 老熟女久久久| 国产精品偷伦视频观看了| 欧美国产精品va在线观看不卡| 蜜桃在线观看..| 精品久久久久久电影网| 欧美性感艳星| 乱码一卡2卡4卡精品| 亚洲精品美女久久久久99蜜臀 | 精品酒店卫生间| 26uuu在线亚洲综合色| 桃花免费在线播放| 亚洲国产精品成人久久小说| 亚洲,欧美,日韩| 少妇人妻久久综合中文| 亚洲四区av| 国产在线一区二区三区精| 美女大奶头黄色视频| 男女国产视频网站| 夜夜骑夜夜射夜夜干| 日韩,欧美,国产一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 蜜桃国产av成人99| 国产福利在线免费观看视频| 亚洲色图综合在线观看| 国产黄色视频一区二区在线观看| 精品99又大又爽又粗少妇毛片| 亚洲一码二码三码区别大吗| 美女xxoo啪啪120秒动态图| 国产爽快片一区二区三区| 日韩人妻精品一区2区三区| 亚洲婷婷狠狠爱综合网| 亚洲欧美成人精品一区二区| 一级毛片 在线播放| 欧美3d第一页| 亚洲精品aⅴ在线观看| 精品少妇久久久久久888优播| 亚洲精品一二三| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 美女福利国产在线| 免费观看av网站的网址| 黄色视频在线播放观看不卡| 在线观看免费日韩欧美大片| 中文字幕亚洲精品专区| 亚洲av男天堂| 成年女人在线观看亚洲视频| 制服人妻中文乱码| 久久久久精品人妻al黑| 好男人视频免费观看在线| 国产精品久久久久久久电影| 久久久久精品性色| 91国产中文字幕| 女的被弄到高潮叫床怎么办| 日韩制服骚丝袜av| 巨乳人妻的诱惑在线观看| 日本猛色少妇xxxxx猛交久久| av视频免费观看在线观看| 欧美性感艳星| 成人亚洲精品一区在线观看| av不卡在线播放| 亚洲一码二码三码区别大吗| 青春草国产在线视频| 成人黄色视频免费在线看| 内地一区二区视频在线| xxx大片免费视频| 成人午夜精彩视频在线观看| 在线观看国产h片| 欧美激情国产日韩精品一区| 久久午夜福利片| 成人亚洲精品一区在线观看| 在线观看美女被高潮喷水网站| 久久久久久人妻| 久久精品熟女亚洲av麻豆精品| 欧美bdsm另类| 超碰97精品在线观看| 精品人妻偷拍中文字幕| 一本久久精品| 亚洲精品美女久久久久99蜜臀 | 国产男人的电影天堂91| 中文字幕制服av| 男女啪啪激烈高潮av片| 亚洲国产成人一精品久久久| 亚洲国产精品专区欧美| 你懂的网址亚洲精品在线观看| 晚上一个人看的免费电影| 亚洲精品视频女| 免费黄色在线免费观看| 成年美女黄网站色视频大全免费| 国产综合精华液| 国产成人一区二区在线| 久久影院123| 男女无遮挡免费网站观看| 我的女老师完整版在线观看| 中文字幕亚洲精品专区| 欧美xxⅹ黑人| 香蕉国产在线看| 久久精品夜色国产| 肉色欧美久久久久久久蜜桃| 9色porny在线观看| 国产又爽黄色视频| 日本vs欧美在线观看视频| 如何舔出高潮| av播播在线观看一区| 狂野欧美激情性bbbbbb| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 亚洲人成77777在线视频| av国产久精品久网站免费入址| 国产成人精品无人区| 建设人人有责人人尽责人人享有的| 日韩精品有码人妻一区| 男女无遮挡免费网站观看| 两个人看的免费小视频| 成人亚洲欧美一区二区av| www日本在线高清视频| av在线观看视频网站免费| 内地一区二区视频在线| 国产激情久久老熟女| 欧美+日韩+精品| 精品国产国语对白av| 美女大奶头黄色视频| 亚洲国产欧美日韩在线播放| 女人精品久久久久毛片| 欧美变态另类bdsm刘玥| 在线观看免费高清a一片| 制服丝袜香蕉在线| 下体分泌物呈黄色| 人成视频在线观看免费观看| www.熟女人妻精品国产 | 超碰97精品在线观看| 大话2 男鬼变身卡| 欧美日韩av久久| 成年美女黄网站色视频大全免费| 国产综合精华液| 午夜影院在线不卡| 少妇被粗大猛烈的视频| 黄色视频在线播放观看不卡| 亚洲激情五月婷婷啪啪| 春色校园在线视频观看| 日韩,欧美,国产一区二区三区| √禁漫天堂资源中文www| 日日啪夜夜爽| 一区在线观看完整版| 狂野欧美激情性bbbbbb| 国产熟女欧美一区二区| 校园人妻丝袜中文字幕| 久久这里有精品视频免费| 七月丁香在线播放| 欧美日韩亚洲高清精品| 天天操日日干夜夜撸| 精品少妇内射三级| 深夜精品福利| 国国产精品蜜臀av免费| 一个人免费看片子| 国产国语露脸激情在线看| 久久久精品区二区三区| 母亲3免费完整高清在线观看 | 国产精品一国产av| 最近2019中文字幕mv第一页| 国产精品久久久久久久电影| 国产免费视频播放在线视频| 18禁国产床啪视频网站| 欧美国产精品一级二级三级| 久久这里有精品视频免费| 成人国产av品久久久| 亚洲精品日韩在线中文字幕| 美女内射精品一级片tv| 国产免费现黄频在线看| 卡戴珊不雅视频在线播放| 亚洲经典国产精华液单| 国产男女内射视频| 精品久久久精品久久久| av.在线天堂| 欧美国产精品va在线观看不卡| 一级片免费观看大全| 激情视频va一区二区三区| freevideosex欧美| 亚洲精品成人av观看孕妇| 国产黄色免费在线视频| 国产深夜福利视频在线观看| 国产爽快片一区二区三区| 国产一区二区在线观看日韩| 欧美成人午夜免费资源| 菩萨蛮人人尽说江南好唐韦庄| 国产高清三级在线| 男女边吃奶边做爰视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产av精品麻豆| 五月天丁香电影| 亚洲激情五月婷婷啪啪| 999精品在线视频| www.色视频.com| 亚洲国产看品久久| 亚洲国产欧美在线一区| 亚洲人成网站在线观看播放| 视频中文字幕在线观看| 成人免费观看视频高清| 高清欧美精品videossex| 日韩一本色道免费dvd| 蜜桃国产av成人99| 亚洲国产精品专区欧美| 国产淫语在线视频| 9色porny在线观看| 90打野战视频偷拍视频| 亚洲精品日本国产第一区| 在线观看免费日韩欧美大片| 高清毛片免费看| 又粗又硬又长又爽又黄的视频| 亚洲人成网站在线观看播放| 91在线精品国自产拍蜜月| 最近中文字幕2019免费版| 春色校园在线视频观看| 九九在线视频观看精品| 亚洲综合精品二区| 午夜日本视频在线| 中文字幕免费在线视频6| 少妇的丰满在线观看| 97超碰精品成人国产| 美女福利国产在线| 纵有疾风起免费观看全集完整版| 丰满饥渴人妻一区二区三| 成人国产麻豆网| 精品国产一区二区三区四区第35| 久久热在线av| 亚洲人成77777在线视频| 国产精品一区www在线观看| 制服诱惑二区| 丝袜喷水一区| 久久久久人妻精品一区果冻| 久久女婷五月综合色啪小说| 人妻少妇偷人精品九色| 国产一区二区三区av在线| 男女午夜视频在线观看 | 亚洲人成网站在线观看播放| 日日摸夜夜添夜夜爱| 满18在线观看网站| 亚洲av综合色区一区| 十分钟在线观看高清视频www| 午夜免费男女啪啪视频观看| 黄色毛片三级朝国网站| 欧美精品一区二区大全|