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

    基于遺傳算法優(yōu)化隨機(jī)森林模型的機(jī)械鉆速分類預(yù)測(cè)方法

    2022-02-05 12:48:26張海軍張高峰王國(guó)娜王立輝劉洋任陽峰鄭雙進(jìn)
    科學(xué)技術(shù)與工程 2022年35期
    關(guān)鍵詞:分類優(yōu)化模型

    張海軍, 張高峰, 王國(guó)娜, 王立輝, 劉洋, 任陽峰, 鄭雙進(jìn)*

    (1.中國(guó)石油天然氣股份有限公司大港油田分公司, 天津 300450; 2.中國(guó)石油天然氣股份有限公司西南油氣田分公司, 成都 610000;3.長(zhǎng)江大學(xué)石油工程學(xué)院, 武漢 430100)

    在油氣鉆探領(lǐng)域,機(jī)械鉆速(rate of penetration,ROP)是評(píng)估鉆井效率的重要指標(biāo),直接影響著鉆井周期和鉆井成本。針對(duì)機(jī)械鉆速,中外學(xué)者大多運(yùn)用傳統(tǒng)方法建立模型進(jìn)行預(yù)測(cè),每種模型均有其適用性。Bahari等[1]建立了W.C.Maurer方程,但該方程考慮因素不夠全面;Adebayo等[2]建立了巖石特性、物理性質(zhì)與機(jī)械鉆速之間的經(jīng)驗(yàn)?zāi)P?,但未考慮鉆井參數(shù)對(duì)機(jī)械鉆速的影響;Kumar等[3]研究得到了一套機(jī)械鉆速預(yù)測(cè)模型,但該模型參數(shù)受巖性影響較大,且模型較為復(fù)雜;Hung等[4]研究得到了基于旋轉(zhuǎn)沖擊鉆頭的機(jī)械鉆速預(yù)測(cè)模型,僅在硬度高的地層中較為適用;巨滿成等[5]、杜鐮等[6]基于巖石可鉆性對(duì)機(jī)械鉆速預(yù)測(cè)進(jìn)行了研究,但研究?jī)H引用了少量的鉆井參數(shù),基于大量統(tǒng)計(jì)數(shù)據(jù)并應(yīng)用各種數(shù)學(xué)方法進(jìn)行處理而建立了鉆速方程,由于不具有普遍性,致使其應(yīng)用受到限制;田璐等[7]通過綜合應(yīng)用油基鉆井液、旋轉(zhuǎn)導(dǎo)向地質(zhì)導(dǎo)向、雙凝雙密度水泥漿等技術(shù),并進(jìn)行了鉆頭優(yōu)選和鉆具組合優(yōu)化,有效提高了機(jī)械鉆速。李琪等[8]提出了一種人工智能算法與結(jié)合粒子群的神經(jīng)網(wǎng)絡(luò)算法(particle swarm optimization-back propagation, PSO-BP)神經(jīng)網(wǎng)絡(luò)相結(jié)合的鉆進(jìn)機(jī)械鉆速預(yù)測(cè)模型,并根據(jù)實(shí)際工況對(duì)其進(jìn)行了模型評(píng)價(jià),結(jié)果表明PSO-BP機(jī)械鉆速模型具有良好的預(yù)測(cè)精度。李昌盛[9]利用多元回歸方法改進(jìn)了B-Y(Bourgoyne-Young)鉆速方程,但該方程只適用于鉆井液鉆井;景寧等[10]提出了一種基于層析分析法和神經(jīng)網(wǎng)絡(luò)相組合的智能模型,利用巖石抗壓強(qiáng)度、鉆頭尺寸等影響因素進(jìn)行鉆速預(yù)測(cè);劉勝娃等[11]研究了一種基于人工神經(jīng)網(wǎng)絡(luò)技術(shù)的定向井機(jī)械鉆速預(yù)測(cè)模型,該模型能在數(shù)據(jù)量較充足、數(shù)據(jù)質(zhì)量較高的條件下得出較高預(yù)測(cè)準(zhǔn)確度;石祥超等[12]評(píng)價(jià)了隨機(jī)森林、支持向量機(jī)、梯度提升樹、人工神經(jīng)網(wǎng)絡(luò)4種人工智能算法在四川盆地某區(qū)塊的機(jī)械鉆速預(yù)測(cè)精度,結(jié)果表明隨機(jī)森林算法對(duì)區(qū)塊內(nèi)各單井?dāng)?shù)據(jù)的預(yù)測(cè)精度能達(dá)到90%,對(duì)整個(gè)區(qū)塊數(shù)據(jù)預(yù)測(cè)的準(zhǔn)確度能達(dá)到88%;左迪一[13]通過采用隨機(jī)森林回歸方法對(duì)5種類別的井建立鉆速預(yù)測(cè)模型并進(jìn)行現(xiàn)場(chǎng)應(yīng)用驗(yàn)證,結(jié)果均取得了較好的效果,說明隨機(jī)森林回歸方法具有非常強(qiáng)的針對(duì)性和準(zhǔn)確度。

    調(diào)研分析發(fā)現(xiàn),現(xiàn)有的機(jī)械鉆速預(yù)測(cè)大多依靠經(jīng)驗(yàn),或是在前人已建立的模型基礎(chǔ)上行改進(jìn),或是依靠控制變量等方法進(jìn)行實(shí)驗(yàn),尋找影響機(jī)械鉆速的主要因素,這些預(yù)測(cè)方法難免會(huì)造成較大誤差,難以滿足當(dāng)前鉆井施工需求。近年來,隨著大數(shù)據(jù)與機(jī)器學(xué)習(xí)技術(shù)的興起,基于歷史鉆井?dāng)?shù)據(jù)運(yùn)用機(jī)器學(xué)習(xí)方法進(jìn)行大數(shù)據(jù)分析,開展機(jī)械鉆速預(yù)測(cè)不失為一種好方法?,F(xiàn)通過分析東部某油田機(jī)械鉆速的影響因素,結(jié)合鉆井歷史數(shù)據(jù)建立了機(jī)械鉆速預(yù)測(cè)模型,并利用遺傳算法進(jìn)行優(yōu)化,測(cè)試模型的精度和泛化性,得到了滿足施工設(shè)計(jì)及現(xiàn)場(chǎng)作業(yè)需要的機(jī)械鉆速預(yù)測(cè)方法,有助于指導(dǎo)該區(qū)塊鉆井施工參數(shù)優(yōu)化,提高鉆井施工效益。

    1 機(jī)械鉆速影響因素分析

    現(xiàn)場(chǎng)鉆井施工過程中,鉆壓、轉(zhuǎn)速、鉆井液密度、循環(huán)排量、漏斗黏度、鉆頭尺寸及牙齒磨損等因素都會(huì)影響機(jī)械鉆速。鉆壓通過影響破巖作用來影響機(jī)械鉆速,如果在硬度較高的地層施加小鉆壓,會(huì)直接影響機(jī)械鉆速的提高[14];轉(zhuǎn)速也會(huì)對(duì)鉆頭破巖產(chǎn)生較大影響,但影響略小于鉆壓,當(dāng)鉆遇較軟地層時(shí)提高轉(zhuǎn)速可明顯提高機(jī)械鉆速,但當(dāng)鉆遇較硬地層時(shí),提高轉(zhuǎn)速并不能很好地提高機(jī)械鉆速[14];鉆井液的密度、排量、漏斗黏度共同影響井底巖屑清理程度和輔助破巖程度,進(jìn)而影響機(jī)械鉆速[15];增大鉆頭尺寸會(huì)降低鉆頭破巖效率,進(jìn)而降低機(jī)械鉆速,但減小鉆頭尺寸又會(huì)影響油氣開采效率,所以鉆頭尺寸設(shè)計(jì)需要權(quán)衡鉆頭尺寸與機(jī)械鉆速的關(guān)系來滿足油田開發(fā)需求[16]。另外,鉆頭牙齒磨損、老化、結(jié)構(gòu)設(shè)計(jì)不合理等因素也會(huì)造成機(jī)械鉆速降低[16]。本文中搜集了來自東部某油田的現(xiàn)場(chǎng)鉆井施工數(shù)據(jù),考慮機(jī)械鉆速影響因素包含鉆壓、轉(zhuǎn)速、鉆頭直徑、鉆頭壓降、鉆頭出入井新度、鉆頭磨損程度、鉆井液密度、漏斗黏度及排量,基于以上影響因素建立滿足施工設(shè)計(jì)及現(xiàn)場(chǎng)作業(yè)需要的機(jī)械鉆速分類預(yù)測(cè)方法。

    2 機(jī)械鉆速數(shù)據(jù)處理與分析

    2.1 數(shù)據(jù)預(yù)處理

    本文建模數(shù)據(jù)來自東部某油田的現(xiàn)場(chǎng)鉆井施工數(shù)據(jù),原始數(shù)據(jù)如表1所示,總計(jì)398組,變量名稱包括入井新度、出井新度、牙齒磨損量、鉆壓、轉(zhuǎn)速、鉆井液排量、鉆井液密度、鉆井液漏斗黏度、鉆頭壓降、鉆頭直徑及平均機(jī)械鉆速。其中以前十組變量作為建立模型時(shí)的輸入變量,平均機(jī)械鉆速為模型的輸出變量。

    2.2 機(jī)械鉆速分級(jí)

    為了提高機(jī)械鉆速分類預(yù)測(cè)的針對(duì)性,需要對(duì)機(jī)械鉆速進(jìn)行分級(jí),鉆速分級(jí)使得機(jī)械鉆速數(shù)值大小不受鉆頭尺寸制約,且將建模任務(wù)從回歸問題轉(zhuǎn)化為分類問題,有利于提高模型的精度。根據(jù)現(xiàn)場(chǎng)實(shí)際工況,將不同鉆頭尺寸對(duì)應(yīng)的機(jī)械鉆速進(jìn)行分級(jí)(“低鉆速”“中鉆速”“高鉆速”),分級(jí)規(guī)則如表2所示。

    按照表2所示的分級(jí)規(guī)則對(duì)表1中的機(jī)械鉆速數(shù)據(jù)按不同鉆頭直徑進(jìn)行分級(jí),分級(jí)后的數(shù)據(jù)如表3所示(低鉆速為1;中鉆速為2;高鉆速為3),為消除各特征之間量綱不同造成的數(shù)值規(guī)模差異,需要對(duì)除機(jī)械鉆速外所有特征數(shù)據(jù)進(jìn)行歸一化處理,歸一化的計(jì)算方法為

    (1)

    表1 順南區(qū)塊現(xiàn)場(chǎng)鉆井施工原始數(shù)據(jù)(前十行)

    表2 不同鉆頭尺寸對(duì)應(yīng)的機(jī)械鉆速分級(jí)

    式(1)中:x′i為歸一化后的數(shù)據(jù);xi表示為歸一化前的原始樣本數(shù)據(jù);xmin表示為原始樣本數(shù)據(jù)特征值的最小值;xmax表示為原始樣本數(shù)據(jù)特征值的最大值。歸一化處理后的施工數(shù)據(jù)見表3所示。

    3 機(jī)械鉆速分類預(yù)測(cè)模型建立

    3.1 機(jī)器學(xué)習(xí)算法原理介紹

    隨機(jī)森林是集成學(xué)習(xí)算法之一[17-19],它的基學(xué)習(xí)器是決策樹;該算法的核心思想是采用多顆決策樹的投票機(jī)制,解決分類和預(yù)測(cè)問題。對(duì)于回歸預(yù)測(cè)問題,將多棵樹的回歸結(jié)果進(jìn)行平均得到最終結(jié)果;對(duì)于分類問題,將多棵樹的判斷結(jié)果進(jìn)行投票,基于少數(shù)服從多數(shù)得到最終的分類結(jié)果;該算法有運(yùn)行速度快,精度高等特點(diǎn)。

    K近鄰算法[20-22](K-nearest neighbor,KNN)寓意K個(gè)靠近的“鄰居”,屬于監(jiān)督式學(xué)習(xí)算法;該算法不會(huì)預(yù)先生成一個(gè)分類預(yù)測(cè)模型,而是將建模與預(yù)測(cè)、分類工作同時(shí)進(jìn)行,模型構(gòu)建好后進(jìn)行分類、預(yù)測(cè)結(jié)果也將輸出;該算法可以對(duì)離散型變量進(jìn)行分類,也可以對(duì)連續(xù)型變量進(jìn)行預(yù)測(cè)。

    支持向量機(jī)[23-25](support vector machine,SVM)的核心思想是結(jié)構(gòu)風(fēng)險(xiǎn)最小化原則和統(tǒng)計(jì)學(xué)習(xí)理論VC(Vapnik-Chervonenkis)維理論,該算法屬于監(jiān)督式機(jī)器學(xué)習(xí)算法,可用于離散因變量的分類和連續(xù)因變量的預(yù)測(cè)。

    3.2 建立機(jī)械鉆速分類預(yù)測(cè)模型

    基于預(yù)處理后的鉆井施工數(shù)據(jù),由Python軟件在數(shù)據(jù)中隨機(jī)選出建模所需的訓(xùn)練集和評(píng)估模型精度所需的測(cè)試集,運(yùn)用隨機(jī)森林算法、K近鄰算法、支持向量機(jī)算法分別建立模型,建模過程全程使用Python自編程序,隨機(jī)森林算法、K近鄰算法、支持向量機(jī)算法通過Sklearn庫實(shí)現(xiàn)。

    表3 機(jī)械鉆速分級(jí)預(yù)處理后的數(shù)據(jù)集(前10行)

    以數(shù)據(jù)集中的鉆頭入井新度、鉆頭出井新度、牙齒磨損量、鉆壓、轉(zhuǎn)速、鉆井液排量、鉆井液密度、鉆井液漏斗黏度、鉆頭壓降、鉆頭直徑作為輸入變量,以平均機(jī)械鉆速作為輸出變量,訓(xùn)練數(shù)據(jù)與測(cè)試數(shù)據(jù)的劃分為9∶1,基于隨機(jī)森林算法、K近鄰算法、支持向量機(jī)算法建立默認(rèn)參數(shù)的機(jī)械鉆速分類預(yù)測(cè)模型;3種算法模型默認(rèn)參數(shù)如表4所示,默認(rèn)模型分類預(yù)測(cè)結(jié)果如表5~表7所示。

    由機(jī)械鉆速等級(jí)預(yù)測(cè)結(jié)果可知,隨機(jī)森林算法、K近鄰算法、支持向量機(jī)算法在39組測(cè)試數(shù)據(jù)集的分類準(zhǔn)確率分別為69.2%、59%、71.8%,分類效果最好的是支持向量機(jī)算法,其次是隨機(jī)森林算法、K近鄰算法??紤]隨機(jī)森林算法、K近鄰算法的分類準(zhǔn)確率不高,需要針對(duì)模型參數(shù)進(jìn)行優(yōu)化,以提高模型預(yù)測(cè)精度。

    表4 3種算法對(duì)應(yīng)的模型默認(rèn)參數(shù)

    表5 隨機(jī)森林算法預(yù)測(cè)機(jī)械鉆速等級(jí)結(jié)果混淆矩陣

    表6 K近鄰算法預(yù)測(cè)機(jī)械鉆速等級(jí)結(jié)果混淆矩陣

    表7 支持向量機(jī)算法預(yù)測(cè)機(jī)械鉆速等級(jí)結(jié)果混淆矩陣

    3.3 模型優(yōu)化與預(yù)測(cè)結(jié)果分析

    遺傳算法(genetic algorithm,GA)[26-28]。是基于自然界遺傳機(jī)制和生物進(jìn)化論的一種高效隨機(jī)搜索和優(yōu)化方法,其具有全局優(yōu)化性能,能夠找到機(jī)器學(xué)習(xí)算法參數(shù)最合適的值,使得算法達(dá)到最優(yōu)配置,從而更加準(zhǔn)確地預(yù)測(cè)機(jī)械鉆速分類。遺傳算法優(yōu)化模型參數(shù)流程圖如圖1所示。

    遺傳算法的種群數(shù)量設(shè)定為200,繁衍最大代數(shù)設(shè)置為80,染色體交叉概率為80%,染色體變異概率為染色體長(zhǎng)度的倒數(shù),適應(yīng)值函數(shù)設(shè)定為39組測(cè)試數(shù)據(jù)在各優(yōu)化模型的準(zhǔn)確率?;谶z傳算法優(yōu)化后的3種算法模型參數(shù)如表8所示。

    圖1 遺傳算法優(yōu)化模型參數(shù)流程圖Fig.1 Flow chart of optimizing model parameters by genetic algorithm

    表8 基于遺傳算法優(yōu)化后的3種算法模型參數(shù)

    將優(yōu)化后的模型參數(shù)輸入對(duì)應(yīng)模型,利用39組測(cè)試數(shù)據(jù)驗(yàn)證優(yōu)化后各模型的精度和泛化性,優(yōu)化后的3個(gè)模型預(yù)測(cè)結(jié)果如表9~表11所示。

    表9 優(yōu)化后的隨機(jī)森林算法預(yù)測(cè)機(jī)械鉆速等級(jí)結(jié)果混淆矩陣

    表10 優(yōu)化后的K近鄰算法預(yù)測(cè)機(jī)械鉆速等級(jí)結(jié)果混淆矩陣

    表11 優(yōu)化后的支持向量機(jī)算法預(yù)測(cè)機(jī)械鉆速 等級(jí)結(jié)果混淆矩陣

    機(jī)械鉆速分類預(yù)測(cè)結(jié)果表明,經(jīng)遺傳算法優(yōu)化后3種算法的準(zhǔn)確率均有提升,隨機(jī)森林算法的分類準(zhǔn)確率提高了12.9%,K近鄰算法的分類準(zhǔn)確率提高了12.8%,支持向量機(jī)算法的分類準(zhǔn)確率提高了5.2%,其中隨機(jī)森林算法的分類準(zhǔn)確率最高為82.1%,其次是支持向量機(jī)算法和K近鄰算法,該方法東部某油田機(jī)械鉆速預(yù)測(cè)提供了新思路。

    4 實(shí)例驗(yàn)證

    4.1 某井基本情況

    該井是東部某油田的一口四開井身結(jié)構(gòu)井,設(shè)計(jì)井深6 190 m。一開φ660.4 mm鉆頭鉆至井深494 m,下入φ508 mm套管,鉆井液密度1.10 g/cm3,漏斗黏度8 s,鉆壓40 kN,轉(zhuǎn)速50 r/min,鉆井液循環(huán)排量60 L/s;二開φ444.5 mm鉆頭鉆至井深2 236 m,下入φ339.7 mm套管,鉆井液密度1.25 g/cm3,漏斗黏度55 s,鉆壓220 kN,轉(zhuǎn)速80 r/min,鉆井液循環(huán)排量65 L/s;三開φ311.2 mm鉆頭鉆至井深4 720 m,下入φ244.5 mm套管,鉆井液密度1.40 g/cm3,漏斗黏度60 s,鉆壓80 kN,轉(zhuǎn)速80 r/min,鉆井液循環(huán)排量60 L/s;四開φ215.9 mm鉆頭鉆至井深6 166 m,下入φ139.7 mm套管,鉆井液密度1.20 g/cm3,漏斗黏度45 s,鉆壓40 kN,轉(zhuǎn)速60 r/min,鉆井液循環(huán)排量30 L/s。該井累計(jì)使用鉆頭27只,其中一開使用鉆頭1只,二開使用鉆頭2只,三開使用鉆頭4只,四開使用鉆頭3只,機(jī)械鉆速范圍在1.03~123.67 m/h。

    4.2 結(jié)果驗(yàn)證分析

    選用皮爾遜相關(guān)系數(shù)法進(jìn)行分析評(píng)價(jià),該方法可用于考量?jī)蓚€(gè)變量X和Y之間的相關(guān)度,其值介于-1~1,1表示變量完全正相關(guān),0表示無關(guān),-1表示完全負(fù)相關(guān)。總體相關(guān)系數(shù)定義為

    (2)

    估算樣本的協(xié)方差和標(biāo)準(zhǔn)差的公式為

    (3)

    圖2 機(jī)械鉆速預(yù)測(cè)值與實(shí)際值相關(guān)性對(duì)比圖Fig.2 Correlation diagram between predicted and actual ROP values

    運(yùn)用遺傳算法優(yōu)化后的隨機(jī)森林模型進(jìn)行機(jī)械鉆速分類預(yù)測(cè),預(yù)測(cè)值與實(shí)際值對(duì)比如圖2所示。從圖2可以看出,皮爾遜相關(guān)系數(shù)r=0.987 3,機(jī)械鉆速分類預(yù)測(cè)值與實(shí)際值之間具有很好的相關(guān)性,二者之間誤差很小,表明運(yùn)用遺傳算法優(yōu)化后的隨機(jī)森林模型可用于東部某油田機(jī)械鉆速分類預(yù)測(cè),根據(jù)多因素敏感性分析可進(jìn)一步開展鉆井施工參數(shù)優(yōu)化。

    5 結(jié)論

    (1)分析了東部某油田機(jī)械鉆速的主要影響因素,基于該油田的鉆井歷史數(shù)據(jù)分別運(yùn)用隨機(jī)森林算法、K近鄰算法及支持向量機(jī)算法建立了機(jī)械鉆速分類預(yù)測(cè)基礎(chǔ)模型,測(cè)試得以上3種算法的分類預(yù)測(cè)準(zhǔn)確率分別為69.2%、59%、71.8%,模型預(yù)測(cè)精度有待提高。

    (2)運(yùn)用遺傳算法針對(duì)隨機(jī)森林算法、K近鄰算法及支持向量機(jī)算法建立的機(jī)械鉆速分類預(yù)測(cè)基礎(chǔ)模型進(jìn)行優(yōu)化,優(yōu)化后3種模型的分類預(yù)測(cè)準(zhǔn)確率分別為:82.1%、71.8%、77%,其中經(jīng)遺傳算法優(yōu)化后的隨機(jī)森林模型精確度最高,可用于該油田機(jī)械鉆速分類預(yù)測(cè)及施工參數(shù)優(yōu)化。

    猜你喜歡
    分類優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    分類算一算
    一道優(yōu)化題的幾何解法
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    分類討論求坐標(biāo)
    數(shù)據(jù)分析中的分類討論
    久久午夜综合久久蜜桃| 国产精品一区二区精品视频观看| 久久久国产成人精品二区| 99热6这里只有精品| 成人手机av| 18禁裸乳无遮挡免费网站照片| 久久天躁狠狠躁夜夜2o2o| 国产精品九九99| 波多野结衣高清作品| 91麻豆精品激情在线观看国产| 丰满的人妻完整版| 搡老岳熟女国产| 成人三级黄色视频| 午夜免费成人在线视频| 波多野结衣高清无吗| 999精品在线视频| 我的老师免费观看完整版| 国语自产精品视频在线第100页| 国产黄片美女视频| 色老头精品视频在线观看| 国语自产精品视频在线第100页| 成年版毛片免费区| 搡老熟女国产l中国老女人| 波多野结衣巨乳人妻| 中文字幕av在线有码专区| 日日摸夜夜添夜夜添小说| 精品国内亚洲2022精品成人| 国产成人精品久久二区二区免费| 搡老熟女国产l中国老女人| 久久精品亚洲精品国产色婷小说| 国产精品国产高清国产av| 在线播放国产精品三级| 国产视频内射| 一本大道久久a久久精品| 色综合婷婷激情| 99久久99久久久精品蜜桃| 美女 人体艺术 gogo| 亚洲中文av在线| 亚洲中文av在线| 亚洲人成77777在线视频| 叶爱在线成人免费视频播放| 亚洲狠狠婷婷综合久久图片| 日韩三级视频一区二区三区| 变态另类丝袜制服| 色综合站精品国产| 变态另类丝袜制服| 嫩草影院精品99| 亚洲av第一区精品v没综合| 视频区欧美日本亚洲| 欧美成人性av电影在线观看| 国模一区二区三区四区视频 | 日韩欧美精品v在线| 精华霜和精华液先用哪个| 国产一区在线观看成人免费| 男人舔女人下体高潮全视频| 久久久久久久久久黄片| 国产欧美日韩精品亚洲av| 一个人免费在线观看的高清视频| 亚洲狠狠婷婷综合久久图片| 日韩三级视频一区二区三区| 欧美性猛交黑人性爽| 欧美性猛交黑人性爽| www.999成人在线观看| 少妇裸体淫交视频免费看高清 | 一进一出抽搐动态| 日本一区二区免费在线视频| 99久久精品国产亚洲精品| 校园春色视频在线观看| 最近最新中文字幕大全电影3| 国产精品一区二区免费欧美| 精品高清国产在线一区| 搡老熟女国产l中国老女人| 757午夜福利合集在线观看| 免费在线观看影片大全网站| 亚洲av成人不卡在线观看播放网| 欧美av亚洲av综合av国产av| 久久精品aⅴ一区二区三区四区| 岛国视频午夜一区免费看| 天堂√8在线中文| 嫩草影院精品99| 亚洲五月天丁香| 1024手机看黄色片| 又爽又黄无遮挡网站| 欧美国产日韩亚洲一区| 国产探花在线观看一区二区| 日日爽夜夜爽网站| 亚洲中文日韩欧美视频| 国产亚洲精品av在线| 白带黄色成豆腐渣| 999久久久精品免费观看国产| 亚洲av成人不卡在线观看播放网| 午夜a级毛片| 床上黄色一级片| 亚洲九九香蕉| 最近最新免费中文字幕在线| 国产欧美日韩一区二区精品| 99久久综合精品五月天人人| 国产精品 欧美亚洲| 中文在线观看免费www的网站 | 香蕉av资源在线| 99在线人妻在线中文字幕| 丝袜美腿诱惑在线| 麻豆av在线久日| 精品欧美国产一区二区三| 久久精品国产99精品国产亚洲性色| 国产v大片淫在线免费观看| 精品久久久久久成人av| 午夜久久久久精精品| 成人欧美大片| 国产亚洲精品一区二区www| 熟女电影av网| 日韩精品青青久久久久久| 国产真人三级小视频在线观看| 黄色女人牲交| 黄色视频不卡| 欧美zozozo另类| 禁无遮挡网站| 亚洲人成伊人成综合网2020| 日本在线视频免费播放| 国产探花在线观看一区二区| 超碰成人久久| 亚洲人成网站在线播放欧美日韩| 国产精品影院久久| 国产精品香港三级国产av潘金莲| 亚洲成人国产一区在线观看| 18禁美女被吸乳视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产一区二区精华液| 窝窝影院91人妻| 国产精品一区二区三区四区免费观看 | 亚洲无线在线观看| 夜夜爽天天搞| 99久久综合精品五月天人人| 日日爽夜夜爽网站| 亚洲精品一卡2卡三卡4卡5卡| 精品电影一区二区在线| 草草在线视频免费看| 极品教师在线免费播放| a在线观看视频网站| 久久婷婷成人综合色麻豆| 色在线成人网| 国产精品久久久久久精品电影| 一二三四社区在线视频社区8| 久久人妻福利社区极品人妻图片| 亚洲精华国产精华精| 免费在线观看成人毛片| 欧美中文综合在线视频| 中文字幕久久专区| 99久久99久久久精品蜜桃| 9191精品国产免费久久| 又紧又爽又黄一区二区| av天堂在线播放| 免费电影在线观看免费观看| 久久久久久大精品| av视频在线观看入口| 精品无人区乱码1区二区| 成人18禁高潮啪啪吃奶动态图| 亚洲自偷自拍图片 自拍| 欧美日韩黄片免| 亚洲片人在线观看| 欧美日韩中文字幕国产精品一区二区三区| videosex国产| 又大又爽又粗| 成人国语在线视频| 丁香六月欧美| 国产久久久一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 岛国在线观看网站| 久久欧美精品欧美久久欧美| 丰满人妻一区二区三区视频av | 两个人看的免费小视频| 欧美乱色亚洲激情| 麻豆久久精品国产亚洲av| 99久久精品热视频| 欧美日韩亚洲综合一区二区三区_| 日本一区二区免费在线视频| 久久久久精品国产欧美久久久| 老司机在亚洲福利影院| av国产免费在线观看| 欧美日韩亚洲综合一区二区三区_| bbb黄色大片| 99热这里只有是精品50| 成人国产综合亚洲| 亚洲真实伦在线观看| 免费高清视频大片| 最近最新免费中文字幕在线| 欧美3d第一页| 国产激情欧美一区二区| 精品少妇一区二区三区视频日本电影| 国产野战对白在线观看| 长腿黑丝高跟| 又大又爽又粗| 色av中文字幕| 99久久99久久久精品蜜桃| 少妇粗大呻吟视频| 一区福利在线观看| 搡老妇女老女人老熟妇| 欧美人与性动交α欧美精品济南到| 国产精品 欧美亚洲| 国产伦一二天堂av在线观看| 午夜福利高清视频| 91麻豆精品激情在线观看国产| 男女床上黄色一级片免费看| 人成视频在线观看免费观看| 97超级碰碰碰精品色视频在线观看| 国产1区2区3区精品| 老汉色∧v一级毛片| 在线免费观看的www视频| 欧美在线一区亚洲| 国产精品美女特级片免费视频播放器 | 性色av乱码一区二区三区2| 好男人在线观看高清免费视频| 国产成人啪精品午夜网站| 无限看片的www在线观看| www.自偷自拍.com| 欧美午夜高清在线| 一级片免费观看大全| 久久草成人影院| 好看av亚洲va欧美ⅴa在| 少妇粗大呻吟视频| 国产精品久久久人人做人人爽| 51午夜福利影视在线观看| 一本一本综合久久| 国产av麻豆久久久久久久| 三级国产精品欧美在线观看 | 国产av一区二区精品久久| 亚洲一码二码三码区别大吗| 精品日产1卡2卡| 脱女人内裤的视频| 999精品在线视频| 久久这里只有精品中国| 欧美日韩乱码在线| 在线观看日韩欧美| 国产成人av激情在线播放| 精品不卡国产一区二区三区| 男女下面进入的视频免费午夜| 搡老妇女老女人老熟妇| 久久久久性生活片| 夜夜躁狠狠躁天天躁| 老司机在亚洲福利影院| 小说图片视频综合网站| 成人三级黄色视频| 午夜福利免费观看在线| 亚洲欧美一区二区三区黑人| 久久天躁狠狠躁夜夜2o2o| 亚洲人成电影免费在线| 国产单亲对白刺激| 亚洲国产精品合色在线| 18禁黄网站禁片免费观看直播| 丝袜美腿诱惑在线| 亚洲人成伊人成综合网2020| 真人做人爱边吃奶动态| 12—13女人毛片做爰片一| 日韩欧美三级三区| 桃红色精品国产亚洲av| 亚洲av电影不卡..在线观看| 成年版毛片免费区| 久久这里只有精品中国| 久久亚洲真实| 特大巨黑吊av在线直播| 欧美激情久久久久久爽电影| 午夜a级毛片| 免费人成视频x8x8入口观看| 成人亚洲精品av一区二区| 欧美日韩精品网址| 男女下面进入的视频免费午夜| 欧美zozozo另类| 蜜桃久久精品国产亚洲av| 欧美不卡视频在线免费观看 | 欧美日韩一级在线毛片| 人妻丰满熟妇av一区二区三区| 国产午夜福利久久久久久| 美女黄网站色视频| 成人手机av| 黑人欧美特级aaaaaa片| 久久香蕉精品热| 亚洲成人国产一区在线观看| 欧美成狂野欧美在线观看| 精品第一国产精品| 成年版毛片免费区| 色噜噜av男人的天堂激情| 制服丝袜大香蕉在线| 女人被狂操c到高潮| 亚洲欧美日韩东京热| 亚洲美女黄片视频| 后天国语完整版免费观看| 国产aⅴ精品一区二区三区波| 少妇粗大呻吟视频| 视频区欧美日本亚洲| 国内精品久久久久久久电影| 亚洲精品av麻豆狂野| 日韩精品青青久久久久久| 国产精品日韩av在线免费观看| 18禁黄网站禁片午夜丰满| 亚洲成av人片在线播放无| 男人舔奶头视频| 香蕉丝袜av| 校园春色视频在线观看| 岛国在线观看网站| 亚洲国产高清在线一区二区三| 日本免费一区二区三区高清不卡| 美女大奶头视频| 久久久水蜜桃国产精品网| 亚洲七黄色美女视频| 亚洲五月婷婷丁香| 国产区一区二久久| 欧美国产日韩亚洲一区| 亚洲精品av麻豆狂野| 成人国产综合亚洲| 黄色 视频免费看| 18禁黄网站禁片免费观看直播| 五月玫瑰六月丁香| 成在线人永久免费视频| 国产精品久久久久久人妻精品电影| 色精品久久人妻99蜜桃| 免费人成视频x8x8入口观看| 1024手机看黄色片| 精品久久蜜臀av无| 国产人伦9x9x在线观看| 黄色毛片三级朝国网站| 在线十欧美十亚洲十日本专区| 国产精品久久久久久久电影 | 国内揄拍国产精品人妻在线| 99热只有精品国产| 久久久久久久午夜电影| 日韩大码丰满熟妇| 欧美一级a爱片免费观看看 | 久久香蕉激情| 99久久99久久久精品蜜桃| 国产在线观看jvid| 99riav亚洲国产免费| 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 在线观看午夜福利视频| 日韩大尺度精品在线看网址| 男人舔奶头视频| 欧美成人性av电影在线观看| 欧美久久黑人一区二区| 身体一侧抽搐| 丝袜美腿诱惑在线| 亚洲国产精品久久男人天堂| 久久精品91蜜桃| or卡值多少钱| 欧美一级毛片孕妇| 亚洲色图av天堂| 欧美黑人欧美精品刺激| 狂野欧美白嫩少妇大欣赏| 亚洲自偷自拍图片 自拍| 国产免费av片在线观看野外av| 成年人黄色毛片网站| 99国产精品99久久久久| 搡老岳熟女国产| 欧美乱妇无乱码| 一本综合久久免费| 一a级毛片在线观看| 午夜免费激情av| 黑人操中国人逼视频| 日本 av在线| a在线观看视频网站| 国产探花在线观看一区二区| 亚洲av成人不卡在线观看播放网| 精品国产乱码久久久久久男人| 成人18禁在线播放| 免费在线观看影片大全网站| 黄频高清免费视频| 制服人妻中文乱码| 黄色 视频免费看| 免费看美女性在线毛片视频| 好男人在线观看高清免费视频| 亚洲人成77777在线视频| aaaaa片日本免费| 婷婷亚洲欧美| 美女高潮喷水抽搐中文字幕| 亚洲免费av在线视频| 亚洲国产欧美网| 国产主播在线观看一区二区| 在线观看日韩欧美| 天天一区二区日本电影三级| 久久天堂一区二区三区四区| 亚洲人与动物交配视频| 无限看片的www在线观看| 美女高潮喷水抽搐中文字幕| 在线观看免费午夜福利视频| 日本免费a在线| 国产视频内射| 女同久久另类99精品国产91| 岛国视频午夜一区免费看| 国产亚洲欧美98| 一边摸一边做爽爽视频免费| 19禁男女啪啪无遮挡网站| 亚洲av电影不卡..在线观看| 三级国产精品欧美在线观看 | 欧美日本亚洲视频在线播放| 中国美女看黄片| 欧美日韩国产亚洲二区| 一进一出好大好爽视频| 可以在线观看的亚洲视频| 国产精品1区2区在线观看.| 天堂√8在线中文| 欧美极品一区二区三区四区| 亚洲精品在线美女| 男女之事视频高清在线观看| 一级毛片女人18水好多| 色噜噜av男人的天堂激情| 午夜精品久久久久久毛片777| 久久久久国产精品人妻aⅴ院| 亚洲国产欧洲综合997久久,| 狂野欧美白嫩少妇大欣赏| 欧美日韩精品网址| 亚洲av电影在线进入| 亚洲男人的天堂狠狠| 一区二区三区激情视频| 香蕉丝袜av| 精品国产乱子伦一区二区三区| 国产成人系列免费观看| 国产成人影院久久av| 天堂av国产一区二区熟女人妻 | 久久久久久久午夜电影| 丰满人妻熟妇乱又伦精品不卡| 神马国产精品三级电影在线观看 | 日韩精品中文字幕看吧| 免费高清视频大片| 欧美中文综合在线视频| 精品电影一区二区在线| 亚洲国产欧美一区二区综合| 窝窝影院91人妻| 国产av一区二区精品久久| 久久久国产成人免费| av免费在线观看网站| 国产精品自产拍在线观看55亚洲| 久久午夜综合久久蜜桃| 97人妻精品一区二区三区麻豆| 一边摸一边抽搐一进一小说| 熟妇人妻久久中文字幕3abv| 免费高清视频大片| 欧美一级毛片孕妇| 亚洲全国av大片| 中亚洲国语对白在线视频| 成人欧美大片| 1024香蕉在线观看| 国产精品香港三级国产av潘金莲| 成人手机av| 不卡av一区二区三区| 国产一区二区在线av高清观看| 成熟少妇高潮喷水视频| 国产三级中文精品| 91国产中文字幕| www日本在线高清视频| 久久久久久久午夜电影| 欧美大码av| www日本黄色视频网| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美激情久久久久久爽电影| 午夜福利免费观看在线| 久久热在线av| av在线播放免费不卡| 97超级碰碰碰精品色视频在线观看| 不卡一级毛片| 国产亚洲精品久久久久5区| 亚洲欧美日韩高清在线视频| 一级a爱片免费观看的视频| 久久精品国产亚洲av高清一级| 日韩三级视频一区二区三区| 可以免费在线观看a视频的电影网站| 999久久久精品免费观看国产| 亚洲精品中文字幕一二三四区| 久久久精品大字幕| 国产成人系列免费观看| 成人av在线播放网站| 午夜激情av网站| 麻豆国产97在线/欧美 | 人妻丰满熟妇av一区二区三区| 亚洲国产精品sss在线观看| 久9热在线精品视频| 黄色成人免费大全| 丝袜美腿诱惑在线| 免费观看精品视频网站| 国产野战对白在线观看| 后天国语完整版免费观看| 国产精品99久久99久久久不卡| 亚洲狠狠婷婷综合久久图片| 国产成人av激情在线播放| 在线观看一区二区三区| 午夜激情福利司机影院| 亚洲人成网站在线播放欧美日韩| 性色av乱码一区二区三区2| 中文字幕高清在线视频| 少妇被粗大的猛进出69影院| 三级毛片av免费| 久久久久久免费高清国产稀缺| а√天堂www在线а√下载| 国产1区2区3区精品| 国产精品av视频在线免费观看| e午夜精品久久久久久久| 国产精品综合久久久久久久免费| 最近最新免费中文字幕在线| 特级一级黄色大片| 九九热线精品视视频播放| 两个人看的免费小视频| 一级a爱片免费观看的视频| 免费在线观看亚洲国产| 悠悠久久av| 麻豆av在线久日| 不卡av一区二区三区| 国产一区二区激情短视频| 久久热在线av| 两个人免费观看高清视频| aaaaa片日本免费| 999久久久国产精品视频| 禁无遮挡网站| 亚洲色图av天堂| 男女之事视频高清在线观看| 国产亚洲精品av在线| av有码第一页| 精品国产乱子伦一区二区三区| 午夜免费观看网址| 狂野欧美激情性xxxx| 神马国产精品三级电影在线观看 | 精品第一国产精品| 小说图片视频综合网站| av中文乱码字幕在线| 国产真人三级小视频在线观看| 久久久久久国产a免费观看| 制服丝袜大香蕉在线| 国产欧美日韩一区二区精品| 一卡2卡三卡四卡精品乱码亚洲| 人妻丰满熟妇av一区二区三区| 韩国av一区二区三区四区| 久久久久国产精品人妻aⅴ院| 我的老师免费观看完整版| 国产伦人伦偷精品视频| 国产成人精品久久二区二区免费| 亚洲熟妇中文字幕五十中出| 啦啦啦观看免费观看视频高清| 久久精品国产清高在天天线| 亚洲天堂国产精品一区在线| 99久久精品热视频| svipshipincom国产片| 国产一区在线观看成人免费| 免费看日本二区| 日本五十路高清| 黄色片一级片一级黄色片| 日韩欧美 国产精品| 在线观看日韩欧美| 一区二区三区高清视频在线| 1024手机看黄色片| 熟女电影av网| 国产精品 国内视频| 又黄又爽又免费观看的视频| 亚洲一码二码三码区别大吗| 欧美性长视频在线观看| 国产成人精品久久二区二区91| 欧美激情久久久久久爽电影| 欧美日韩乱码在线| 天天一区二区日本电影三级| 51午夜福利影视在线观看| 国产午夜精品论理片| 搡老妇女老女人老熟妇| 无限看片的www在线观看| 叶爱在线成人免费视频播放| 又黄又粗又硬又大视频| 国产亚洲精品久久久久5区| 91九色精品人成在线观看| 99国产综合亚洲精品| 熟妇人妻久久中文字幕3abv| 性色av乱码一区二区三区2| 不卡av一区二区三区| 精品免费久久久久久久清纯| 又黄又爽又免费观看的视频| 99久久无色码亚洲精品果冻| av福利片在线观看| 久久久久九九精品影院| 亚洲电影在线观看av| 国产av又大| 欧美+亚洲+日韩+国产| 女人被狂操c到高潮| 久久久国产欧美日韩av| 国产视频内射| 国产精品久久久av美女十八| 色av中文字幕| 他把我摸到了高潮在线观看| 两个人看的免费小视频| 国产高清视频在线观看网站| 窝窝影院91人妻| 天天添夜夜摸| 国产午夜精品论理片| 51午夜福利影视在线观看| 亚洲av五月六月丁香网| 欧美成人一区二区免费高清观看 | 欧美成人午夜精品| 嫁个100分男人电影在线观看| 日韩中文字幕欧美一区二区| 亚洲成av人片免费观看| 成人av一区二区三区在线看| 在线观看日韩欧美| 亚洲精品在线观看二区| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 久久伊人香网站| 美女黄网站色视频| 窝窝影院91人妻| 亚洲av电影在线进入| 在线十欧美十亚洲十日本专区| 大型av网站在线播放| 久热爱精品视频在线9| a在线观看视频网站| 久久久久久人人人人人| 欧美色欧美亚洲另类二区| 美女大奶头视频| 国产av一区在线观看免费| 美女黄网站色视频| 日日摸夜夜添夜夜添小说| 成在线人永久免费视频| 日韩精品中文字幕看吧| 嫁个100分男人电影在线观看|