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

    基于改進(jìn)量子遺傳算法的油田井位及數(shù)量優(yōu)化

    2019-03-05 10:23:14,,
    計(jì)算機(jī)測量與控制 2019年2期
    關(guān)鍵詞:井位油藏遺傳算法

    ,,

    (1.長江大學(xué) 石油工程學(xué)院,武漢 430100;2.湖北工業(yè)大學(xué) 電氣與電子工程學(xué)院,武漢 430068)

    0 引言

    為了提高油田開發(fā)效果以及采收率,許多技術(shù),如注水驅(qū)油、二氧化碳驅(qū)油以及化學(xué)驅(qū)油等被陸續(xù)開發(fā)出來,這些技術(shù)都有利于于從油田儲層中提取更多的石油,在開采技術(shù)的提高創(chuàng)新之外,適當(dāng)?shù)木慌c合理的井?dāng)?shù)也是盡可能地提取石油并提高采收率的重要途徑,因此,為了應(yīng)對不斷增長的世界能源需求,布井優(yōu)化一直是行業(yè)內(nèi)眾多研究的焦點(diǎn)。

    井位的優(yōu)劣受油藏地質(zhì)條件、流體特性以及開發(fā)時機(jī)等一系列復(fù)雜因素的影響,傳統(tǒng)的井位選擇方式是由人為制定若干并網(wǎng)方案,再采用數(shù)值模擬的方法選出一個較優(yōu)的方案,該方案在理論上十分簡便,但十分依賴方案制定者的經(jīng)驗(yàn),且從方案的制定到評價(jià)過程都十分耗時。布井的另一個問題是確定合理的井?dāng)?shù),作為產(chǎn)油的直接渠道,油井的多少直接決定了油田的經(jīng)濟(jì)效益,油井太少,會使油田開發(fā)速度低下,無法滿足生產(chǎn)任務(wù);井?dāng)?shù)太多,油田的建設(shè)與維護(hù)成本也隨之增加,且大密度的布井會造成油井間的相互干擾,降低油田經(jīng)濟(jì)效益。所以布井的優(yōu)化問題往往是非線性且多模態(tài)的[1]。針對這種復(fù)雜的優(yōu)化問題,建立合適的數(shù)學(xué)模型通過智能優(yōu)化算法來進(jìn)行求解是一種主流的處理方法。

    目前,研究者們采用了許多各具特色的優(yōu)化算法來解決布井優(yōu)化問題。2010年,Bukhamsin[2]等,提出了一種多段井井位的優(yōu)化方法,并采用了連續(xù)遺傳算法(CGA);2011年,Morales[3]等在利用遺傳算法對井網(wǎng)優(yōu)化的同時,引入了風(fēng)險(xiǎn)約束,使井位的分布更加合理;2013年,安瑤[6]等人運(yùn)用了一種改進(jìn)的粒子群算法對井位進(jìn)行優(yōu)化,并與遺傳算法做了比較,證明了在速度方面PSO優(yōu)于GA;同年,Lyons[4]等,將卡爾曼濾波與遺傳算法相結(jié)合,在考慮了時間不確定性的情況下對井位進(jìn)行了優(yōu)化2015年,為了進(jìn)一步提高優(yōu)化效率,Sampaio[5]等提出一種混合算法,該算法將GA與共軛梯度算法(CC)相結(jié)合,算法參數(shù)設(shè)計(jì)合適的條件下,該算法具有良好的全局與局部最優(yōu)的性能。從上述研究進(jìn)程可以看出,在加速井位優(yōu)化進(jìn)程方面,結(jié)合多種機(jī)制的混合算法極具競爭力。

    量子遺傳算法(QGA)作為量子計(jì)算與進(jìn)化算法結(jié)合的產(chǎn)物,量子理論中量子態(tài)相干、疊加與糾纏等特性大大 了傳統(tǒng)算法中欠缺的并行計(jì)算能力,其計(jì)算速度與收斂效果均優(yōu)于傳統(tǒng)遺傳算法,目前已逐步替代傳統(tǒng)遺傳算法應(yīng)用于許多工程領(lǐng)域。因此本文在利用Matlab建立油藏?cái)?shù)值模擬軟件建立模型,將井?dāng)?shù)和井位作為變量,以油田凈現(xiàn)值為目標(biāo)函數(shù)結(jié)合QGA對井位進(jìn)行優(yōu)化研究。

    1 模型建立

    1.1 地質(zhì)模型

    為了精細(xì)描述油藏資源,本文利用Matlab軟件在20×20網(wǎng)格上建立了一個非均質(zhì)油藏模擬模型,該模型為一層,且厚度為50 m,X,Y方向上的網(wǎng)絡(luò)步長為100 m。該油藏在開采時主要依靠其中所儲存的油、水的彈性驅(qū)動,在沒有外界能量供給的條件下將原油從地層驅(qū)至井底。本文在該模型的基礎(chǔ)上以油田累計(jì)凈現(xiàn)值最大化為目標(biāo),通過算法運(yùn)算選取合適的布井位置及數(shù)量。油藏模擬模型如圖1所示。

    圖1中油藏滲透率的變化范圍為100~400 md,其分布場中各節(jié)點(diǎn)的滲透率值由Kriging(克里金)插值法生成。

    圖1 油藏地址模型

    1.2 數(shù)學(xué)模型

    在地質(zhì)模型的基礎(chǔ)上,布井優(yōu)化的數(shù)學(xué)模型將以井位及井?dāng)?shù)為變量并結(jié)合合理的目標(biāo)函數(shù)共同構(gòu)成。

    本文的目標(biāo)函數(shù)模型根據(jù)油田累計(jì)凈現(xiàn)值建立。凈現(xiàn)值(NPV)指在項(xiàng)目建設(shè)以及一定的生成年限內(nèi)的現(xiàn)金流量根據(jù)折現(xiàn)率、通貨膨脹率折算后得到的具體數(shù)值。本文中建立的NPV表達(dá)式為:

    (1)

    式(1)中,T為生產(chǎn)年限;δ為通貨膨脹率;η為基準(zhǔn)折現(xiàn)率;Cin(t)、Cout(t)為第t年的總收入及支出(元),二者相減則為油田在第t年的凈現(xiàn)金流量。其中:

    Cin=PoQ(φ,k,x,y,n)

    (2)

    式(2)中,Po為原油價(jià)格;Q原油產(chǎn)量,根據(jù)油田實(shí)際開發(fā)過程中的影響因素,原油產(chǎn)量的決策變量主要由孔隙度φ、滲透率k、井位坐標(biāo)x,y以及井?dāng)?shù)n構(gòu)成。

    Cout=nPm(Pinit,Pop,Pre)

    (3)

    式(3)中,Pm為總成本,其中包括建設(shè)成本Pinit、操作成本Pop以及維護(hù)成本Pre,為簡化問題的數(shù)學(xué)描述,本文以單個井為單位計(jì)算各類成本。

    通過上述分析,本文設(shè)立的目標(biāo)函數(shù)最終形式為:

    (4)

    同時為了提高油田的安全性,消除油田生產(chǎn)工作時井間的相互影響,在布井時應(yīng)當(dāng)保持一定的安全距離。故本文將井距限制作為井位優(yōu)化的約束條件,以井泄油半徑作為井距限制,則有:

    (5)

    其中:dij為任意兩井i、j之間的距離,Rj、Ri分別為井i、j的泄油半徑。

    2 改進(jìn)的量子遺傳算法

    2.1 量子遺傳算法簡介

    QGA作為量子計(jì)算與遺傳算法相融合的產(chǎn)物,是一種結(jié)合了量子進(jìn)化算法的多樣性和遺傳算法的全局尋優(yōu)能力的混合啟發(fā)式算法。

    在QGA中,染色體不再由傳統(tǒng)的二進(jìn)制數(shù)或浮點(diǎn)數(shù)來表示,而是采用了量子比特的方式來表示,其特點(diǎn)在于該方法能同時表達(dá)任意疊加態(tài),即單個量子位并不表示一個確定的數(shù)值,而是可以同時表示0和1以及他們之間任意的疊加狀態(tài)。量子比特的染色體結(jié)構(gòu)可以表示為:

    (6)

    量子旋轉(zhuǎn)門作為QGA中的個體進(jìn)化策略,可根據(jù)具體問題進(jìn)行選擇,目前學(xué)者們針對不同的問題設(shè)計(jì)了許多量子選擇門,其基本工作原理如下:

    (7)

    式中,αi和βi表示第t代染色體上第i量子比特通過旋轉(zhuǎn)門更新前后的概率幅;θi為旋轉(zhuǎn)角,其值根據(jù)算法策略事先設(shè)定,在算法進(jìn)行時,將根據(jù)當(dāng)前解與當(dāng)前最優(yōu)解之間相對應(yīng)的部分基因的關(guān)系以查表的形式來選擇合適的值替換。

    與傳統(tǒng)GA相比,量子比特的概率幅表示使得QGA更具并行性,在增加種群多樣性的同時加快了算法的求解速度;量子門的種群更新方式也能在保持種群多樣性的同時避免選擇壓力問題,提高了算法的搜索效率,使算法擁有更好的收斂特性。

    2.2 改進(jìn)的量子遺傳算法

    為了使量子遺傳算法更好地解決所研究的問題,本文在標(biāo)準(zhǔn)量子遺傳算法的基礎(chǔ)上做出了改進(jìn)。

    量子遺傳算法是進(jìn)化算法與量子計(jì)算的結(jié)合,它對于初始數(shù)據(jù)的生成、數(shù)據(jù)的處理等方面偏向于進(jìn)化算法的處理方式,而進(jìn)化算法對于數(shù)據(jù)的處理方面具有一定的隨機(jī)性,就本文數(shù)學(xué)模型的約束條件而言,對于隨機(jī)生成的多個坐標(biāo),很難制定一條簡便規(guī)則來保持它們兩兩之間的距離在一個定值之外。所以,為了簡化算法,加快算法的計(jì)算速度,本文采取的改進(jìn)方案受到蟻群算法的啟發(fā),將禁忌表策略與量子遺傳算法相結(jié)合,使算法更好的處理井位的優(yōu)化問題。

    禁忌表是禁忌搜索的一個重要應(yīng)用領(lǐng)域,其主要思想是在算法進(jìn)程中,標(biāo)記已搜索的局部里被約束條件所限制的對象,并在進(jìn)一步的迭代搜索中避開這些對象。這種策略在本文中將嵌入到井位優(yōu)化算法的初始種群生成以及種群進(jìn)化的過程里,其方法如下步驟如下:

    1)算法根據(jù)預(yù)設(shè)井?dāng)?shù)n,依次隨機(jī)生成n個井位坐標(biāo)。當(dāng)?shù)谝粋€井位坐標(biāo)生成時,將位于該坐標(biāo)井距限制內(nèi)的其余坐標(biāo)全部納入禁忌表Tabu;

    2)當(dāng)生成下一個井位坐標(biāo)時,算法將自動排除掉禁忌表Tabu中的點(diǎn),在剩余坐標(biāo)點(diǎn)中繼續(xù)隨機(jī)生成。當(dāng)?shù)诙€井位坐標(biāo)產(chǎn)生后,根據(jù)步驟1)的方法,繼續(xù)擴(kuò)充禁忌表Tabu;

    3)重復(fù)步驟2)直到井位數(shù)為n。

    以上步驟是禁忌表策略在初始種群生成時的應(yīng)用,在種群進(jìn)化時,該步驟同樣適用,只不過由于算法需要根據(jù)適應(yīng)度函數(shù)(目標(biāo)函數(shù))朝著求解的最優(yōu)的方向進(jìn)行,故步驟2)中的坐標(biāo)點(diǎn)不再是隨機(jī)生成,而是改為在進(jìn)化更新后的坐標(biāo)向量中,以第一坐標(biāo)點(diǎn)為基準(zhǔn)建立禁忌表,其余不符合約束條件的坐標(biāo)在禁忌表以外的坐標(biāo)點(diǎn)中選取最近的坐標(biāo)進(jìn)行替換。

    本文利用改進(jìn)后的量子遺傳算法優(yōu)化井位及數(shù)量的流程如圖2所示。

    圖2 改進(jìn)后算法流程圖

    1)算法中的參數(shù)設(shè)置包括油藏參數(shù)、算法參數(shù)設(shè)置。油藏參數(shù)主要包括油藏?cái)?shù)據(jù)以及預(yù)設(shè)油井個數(shù),油藏?cái)?shù)據(jù)主要在目標(biāo)函數(shù)的計(jì)算中使用,通過計(jì)算得數(shù)來判斷井位坐標(biāo)的優(yōu)劣,預(yù)設(shè)油井個數(shù)則決定了算法生成初始坐標(biāo)的個數(shù);算法參數(shù)則主要決定了算法迭代次數(shù)與運(yùn)算速度。

    2)針對本問題是以NPV最大為目的進(jìn)行優(yōu)化,故算法中的全局最優(yōu)解即為特定井?dāng)?shù)下一組NPV最大的井位坐標(biāo),算法每經(jīng)過一次迭代都會通過NPV對當(dāng)代種群進(jìn)行評價(jià),并選取其中NPV最大的解與歷史全局最優(yōu)解進(jìn)行比較,選擇其中較大的值作為當(dāng)前全局最優(yōu)解,不斷循環(huán),直到算法停止迭代時輸出該解。

    3 仿真實(shí)驗(yàn)與結(jié)果分析

    3.1 實(shí)驗(yàn)參數(shù)設(shè)置

    本文在算法運(yùn)行時涉及到的計(jì)算參數(shù)有以下幾種:

    生產(chǎn)年限T:T= 5 (年);

    孔隙度φ:φ=[12.0156+2.2613logk(x,y)]/100

    壓縮系數(shù)Cp:Cp=φ×10-3+10-6(MPa-1);

    油層厚度h:h= 50(米);

    油層粘度μ:μ= 5.0(mPa·s);

    生產(chǎn)壓差ΔP:ΔP=10(MPa);

    泄油半徑R:R=100(m)

    主要經(jīng)濟(jì)參數(shù)有:建設(shè)成本Pinit100萬元;操作成本Pop200萬元;維護(hù)成本Pre50萬元;原油價(jià)格Po500元/m3;折現(xiàn)率η=0.1;通貨膨脹率δ=0.04。

    主要算法參數(shù)有:種群大小100;進(jìn)化代數(shù)100。

    3.2 實(shí)驗(yàn)結(jié)果及分析

    本文在分別在擬定的生產(chǎn)井?dāng)?shù)n={1 2 3…10}的條件下,對模型進(jìn)行了多次運(yùn)算,分別計(jì)算出了不同的擬定井?dāng)?shù)下,5年時間后的NPV值以及最佳井位坐標(biāo)。其中不同井?dāng)?shù)與其對應(yīng)的最優(yōu)NPV值如圖2所示。根據(jù)圖2曲線可以看出,針對該油藏模型通過算法優(yōu)化后的合理井?dāng)?shù)為5口。當(dāng)井?dāng)?shù)小于5時,NPV隨著井?dāng)?shù)的增加而增加;當(dāng)井?dāng)?shù)超過5口時,NPV隨著井?dāng)?shù)的增加而減小,這是因?yàn)殡S著井?dāng)?shù)的增加,雖然油田的累積產(chǎn)量會增加但其增長幅度小于成本的增長幅度,因此會導(dǎo)致NPV的減小。

    部分井?dāng)?shù)最優(yōu)坐標(biāo)如圖3所示。可以看出最優(yōu)的布井位置與油田滲透率的分布有密切聯(lián)系,滲透率大的坐標(biāo)區(qū)域內(nèi)井?dāng)?shù)較多,而滲透率小的坐標(biāo)區(qū)域內(nèi)井?dāng)?shù)很少。通過圖4(a)、(b)、(c)、(d)的對比可以發(fā)現(xiàn),在優(yōu)化布井時滲透率高的坐標(biāo)點(diǎn)總是被優(yōu)先考慮,例如,在坐標(biāo)(11,16)附近的區(qū)域內(nèi)點(diǎn),無論總井?dāng)?shù)多少均會被選中。同時也可以發(fā)現(xiàn),當(dāng)井?dāng)?shù)較多時優(yōu)化后的井位較為分散,這是由于本文所研究的油藏模型的非均質(zhì)程度較大,滲透率分布不均勻。對于該種地質(zhì)區(qū)塊,算法優(yōu)化布井比常規(guī)布井更具優(yōu)勢,下面將具體說明。

    圖3 井?dāng)?shù)-NPV關(guān)系曲線

    圖4 部分井?dāng)?shù)對應(yīng)的最優(yōu)井位

    為了進(jìn)一步說明算法優(yōu)化布井的優(yōu)勢,本文將傳統(tǒng)布井方式與優(yōu)化后的布井方式進(jìn)行了對比,以井?dāng)?shù)為5、6為例,傳統(tǒng)布井方式如圖4所示。

    圖5 傳統(tǒng)布井方法

    通過優(yōu)化前后的NPV來說明優(yōu)化效果,NVP的計(jì)算年限同樣為5年,在不同預(yù)設(shè)井?dāng)?shù)下,傳統(tǒng)布井方式的NPV與優(yōu)化后的NPV如表1所示。

    由表1可以看出,使用優(yōu)化布井在五年內(nèi)所獲得的收益在任何井?dāng)?shù)上均優(yōu)于常規(guī)布井,尤其在最佳布井?dāng)?shù)(井?dāng)?shù)為5)時,增幅達(dá)到最大。同時,從常規(guī)布井的NPV值隨井?dāng)?shù)變化可以看出,當(dāng)井?dāng)?shù)達(dá)到6口時,繼續(xù)開井反而會使NPV值降低,這從說明了另一方面說明了無論采取任何布井方式,合理規(guī)劃井?dāng)?shù)都是一定年限內(nèi)取得最佳收益的重要途徑。

    表1 常規(guī)布井與優(yōu)化布井NPV比較

    4 結(jié)束語

    本文研究將改進(jìn)量子遺傳算法應(yīng)用到油田井位數(shù)量及位置優(yōu)化中,其中目標(biāo)函數(shù)為凈現(xiàn)值。并通過數(shù)值模擬的方法在Matlab上建立仿真模型進(jìn)行實(shí)驗(yàn),通過優(yōu)化實(shí)驗(yàn)可知,基于量子遺傳算法的井位及數(shù)量優(yōu)化能在知曉油藏模型的前提下確定油田最優(yōu)井位及數(shù)量,不僅克服了常規(guī)布井對經(jīng)驗(yàn)的依賴性,還通過優(yōu)化NPV前后的NPV值體現(xiàn)了優(yōu)化方案良好的經(jīng)濟(jì)性。

    由于本次研究是通過模擬油藏模型來驗(yàn)證優(yōu)化方案的可行性,為了簡化問題的數(shù)學(xué)描述,本次研究的約束條件僅為井間距離,但在實(shí)際應(yīng)用中,許多因素例如:地質(zhì)條件、驅(qū)動方式、流體特殊特性等等都會對油田的NPV值造成影響。因此進(jìn)一步結(jié)合實(shí)際情況,提升優(yōu)化方案的泛化性是本次研究的主要后繼工作,需要繼續(xù)深入研究。

    猜你喜歡
    井位油藏遺傳算法
    基于改進(jìn)Nelder-Mead 算法的井位優(yōu)選方法
    基于Microstation二次開發(fā)在油田井位展繪的應(yīng)用
    河北省行唐縣貧水山區(qū)物探找水研究
    渤海油田側(cè)鉆井槽口優(yōu)選設(shè)計(jì)方法
    海洋石油(2021年1期)2021-03-30 13:28:06
    低滲油藏注采結(jié)構(gòu)調(diào)整與挖潛技術(shù)
    云南化工(2020年11期)2021-01-14 00:51:02
    基于模糊數(shù)學(xué)的油藏干層識別研究
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    基于遺傳算法和LS-SVM的財(cái)務(wù)危機(jī)預(yù)測
    基于改進(jìn)的遺傳算法的模糊聚類算法
    一区二区三区国产精品乱码| 人成视频在线观看免费观看| 久久久久久久午夜电影| 一区二区三区激情视频| 精品久久久久久久久久久久久| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久人人做人人爽| 757午夜福利合集在线观看| 精品午夜福利视频在线观看一区| 99国产精品一区二区蜜桃av| 午夜a级毛片| 在线观看免费视频日本深夜| 日韩欧美国产在线观看| 这个男人来自地球电影免费观看| 黄片小视频在线播放| 悠悠久久av| 香蕉久久夜色| 特大巨黑吊av在线直播| 一级片免费观看大全| 亚洲成人国产一区在线观看| 欧美三级亚洲精品| 看免费av毛片| 亚洲成a人片在线一区二区| 欧美日韩乱码在线| avwww免费| 中文亚洲av片在线观看爽| 啦啦啦观看免费观看视频高清| 法律面前人人平等表现在哪些方面| 日本在线视频免费播放| 亚洲aⅴ乱码一区二区在线播放 | 国产欧美日韩精品亚洲av| 午夜福利高清视频| 国产亚洲精品av在线| 精品日产1卡2卡| 欧美成人一区二区免费高清观看 | 午夜福利欧美成人| 亚洲熟女毛片儿| 久久中文看片网| 成人一区二区视频在线观看| 高清在线国产一区| 亚洲 欧美一区二区三区| 中文字幕久久专区| 中文资源天堂在线| 老司机深夜福利视频在线观看| 精品不卡国产一区二区三区| 激情在线观看视频在线高清| 色哟哟哟哟哟哟| 人人妻人人澡欧美一区二区| 亚洲av美国av| 久久亚洲真实| 精品一区二区三区四区五区乱码| 老司机午夜福利在线观看视频| 日日爽夜夜爽网站| 搡老熟女国产l中国老女人| 亚洲电影在线观看av| 精品久久久久久,| 国产午夜精品久久久久久| 成人午夜高清在线视频| 视频区欧美日本亚洲| 午夜福利在线观看吧| 女警被强在线播放| 亚洲av电影不卡..在线观看| 国内精品久久久久精免费| 欧美另类亚洲清纯唯美| 丁香六月欧美| 看黄色毛片网站| 午夜激情av网站| 亚洲欧美日韩无卡精品| 欧美一区二区国产精品久久精品 | 69av精品久久久久久| 老鸭窝网址在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美日韩一级在线毛片| 非洲黑人性xxxx精品又粗又长| 成人国语在线视频| 亚洲国产中文字幕在线视频| 欧美久久黑人一区二区| 黄色a级毛片大全视频| 99国产精品99久久久久| 欧美乱码精品一区二区三区| 国产精品一区二区三区四区免费观看 | 亚洲精品国产一区二区精华液| 欧美一级a爱片免费观看看 | 亚洲精品国产精品久久久不卡| 99在线视频只有这里精品首页| 黄色a级毛片大全视频| 99久久久亚洲精品蜜臀av| 欧美黑人欧美精品刺激| 精品熟女少妇八av免费久了| 丁香六月欧美| 男人舔奶头视频| 欧美性猛交黑人性爽| 国内精品久久久久久久电影| 91在线观看av| 国产蜜桃级精品一区二区三区| 日韩大码丰满熟妇| 一级毛片精品| 亚洲乱码一区二区免费版| 我要搜黄色片| 伦理电影免费视频| 日韩欧美 国产精品| 亚洲av日韩精品久久久久久密| 亚洲欧美激情综合另类| 国产午夜精品论理片| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品成人综合色| 精品日产1卡2卡| 国产乱人伦免费视频| 亚洲成人国产一区在线观看| 老司机午夜十八禁免费视频| 老熟妇乱子伦视频在线观看| 亚洲熟妇中文字幕五十中出| 欧美色欧美亚洲另类二区| 一本大道久久a久久精品| 午夜免费成人在线视频| 啪啪无遮挡十八禁网站| 黑人欧美特级aaaaaa片| 人人妻,人人澡人人爽秒播| 制服丝袜大香蕉在线| 蜜桃久久精品国产亚洲av| 亚洲色图av天堂| 一进一出好大好爽视频| 一本一本综合久久| 亚洲一区二区三区不卡视频| 久久亚洲精品不卡| 欧美乱妇无乱码| 99精品在免费线老司机午夜| 亚洲国产精品sss在线观看| 亚洲黑人精品在线| 免费在线观看日本一区| 日韩欧美精品v在线| 最近最新中文字幕大全免费视频| 91老司机精品| 最好的美女福利视频网| 亚洲精品中文字幕一二三四区| 精品午夜福利视频在线观看一区| 成人特级黄色片久久久久久久| 欧美黑人欧美精品刺激| 成年女人毛片免费观看观看9| 香蕉久久夜色| 久久久久久久久久黄片| 日韩有码中文字幕| 日本熟妇午夜| 亚洲国产欧洲综合997久久,| 亚洲成人久久爱视频| 给我免费播放毛片高清在线观看| 亚洲七黄色美女视频| 欧美成人一区二区免费高清观看 | 特大巨黑吊av在线直播| 老司机深夜福利视频在线观看| 国产三级中文精品| 亚洲精品色激情综合| 久久久国产精品麻豆| 国模一区二区三区四区视频 | x7x7x7水蜜桃| 久久伊人香网站| 亚洲专区字幕在线| 国产高清视频在线播放一区| 亚洲片人在线观看| 欧美黑人欧美精品刺激| 亚洲黑人精品在线| aaaaa片日本免费| 最新美女视频免费是黄的| 欧美激情久久久久久爽电影| 亚洲成人久久性| 国产免费av片在线观看野外av| 老汉色av国产亚洲站长工具| 69av精品久久久久久| 亚洲精品久久国产高清桃花| 午夜激情福利司机影院| 天堂影院成人在线观看| 欧美+亚洲+日韩+国产| 两个人视频免费观看高清| 亚洲天堂国产精品一区在线| 欧美日本亚洲视频在线播放| 99久久无色码亚洲精品果冻| 啦啦啦观看免费观看视频高清| 欧美日本亚洲视频在线播放| 欧美成狂野欧美在线观看| svipshipincom国产片| 国产亚洲精品久久久久久毛片| 哪里可以看免费的av片| 一卡2卡三卡四卡精品乱码亚洲| 日本黄大片高清| 亚洲性夜色夜夜综合| www.熟女人妻精品国产| 精品久久久久久,| 这个男人来自地球电影免费观看| 一本综合久久免费| 亚洲男人天堂网一区| 欧美色欧美亚洲另类二区| 精品福利观看| videosex国产| 精品一区二区三区av网在线观看| 精品久久久久久,| 日韩高清综合在线| 亚洲成av人片在线播放无| 欧美黑人精品巨大| 久99久视频精品免费| 色老头精品视频在线观看| АⅤ资源中文在线天堂| 亚洲男人的天堂狠狠| 亚洲电影在线观看av| 久久人妻福利社区极品人妻图片| 中文字幕久久专区| 国产精品久久久人人做人人爽| 国产成人av教育| 国产乱人伦免费视频| 国产亚洲欧美在线一区二区| 午夜久久久久精精品| 婷婷精品国产亚洲av在线| 免费搜索国产男女视频| 成人永久免费在线观看视频| 免费高清视频大片| 可以免费在线观看a视频的电影网站| 精品一区二区三区视频在线观看免费| 国产成人精品无人区| 国产高清videossex| www.自偷自拍.com| 国产亚洲av嫩草精品影院| 熟妇人妻久久中文字幕3abv| 久久久国产精品麻豆| 亚洲第一电影网av| 99re在线观看精品视频| 国产视频内射| 日韩欧美 国产精品| 日韩av在线大香蕉| 黑人操中国人逼视频| 丰满人妻一区二区三区视频av | 丝袜美腿诱惑在线| 精品国内亚洲2022精品成人| av欧美777| 成人av在线播放网站| 国产黄片美女视频| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 波多野结衣高清作品| 亚洲第一欧美日韩一区二区三区| 成人av一区二区三区在线看| 免费搜索国产男女视频| 超碰成人久久| 国产午夜精品久久久久久| 法律面前人人平等表现在哪些方面| www.自偷自拍.com| 两性夫妻黄色片| 日本熟妇午夜| a级毛片a级免费在线| 成人18禁高潮啪啪吃奶动态图| 国产野战对白在线观看| 久99久视频精品免费| 久热爱精品视频在线9| 亚洲人成77777在线视频| 久久久久精品国产欧美久久久| av天堂在线播放| 午夜精品在线福利| 亚洲熟女毛片儿| 国产高清有码在线观看视频 | 亚洲国产欧洲综合997久久,| 亚洲欧洲精品一区二区精品久久久| 午夜免费成人在线视频| 黄色片一级片一级黄色片| 亚洲成人免费电影在线观看| 亚洲黑人精品在线| 欧美日韩瑟瑟在线播放| 俄罗斯特黄特色一大片| 日韩欧美在线乱码| 88av欧美| 黄色丝袜av网址大全| 欧美av亚洲av综合av国产av| 欧美色欧美亚洲另类二区| 999精品在线视频| a级毛片在线看网站| 91在线观看av| 国产精品爽爽va在线观看网站| 国产又黄又爽又无遮挡在线| 久久亚洲精品不卡| 欧美日韩亚洲国产一区二区在线观看| www.自偷自拍.com| 在线观看日韩欧美| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲av嫩草精品影院| 欧美日韩瑟瑟在线播放| 日本免费a在线| 国产av麻豆久久久久久久| 久久精品国产清高在天天线| 亚洲av电影不卡..在线观看| 亚洲自拍偷在线| 真人一进一出gif抽搐免费| 亚洲av五月六月丁香网| 听说在线观看完整版免费高清| 最近最新中文字幕大全电影3| 久久久国产欧美日韩av| 亚洲国产精品sss在线观看| 免费看十八禁软件| 一本大道久久a久久精品| 19禁男女啪啪无遮挡网站| 99热这里只有精品一区 | avwww免费| 亚洲成a人片在线一区二区| 99久久久亚洲精品蜜臀av| 欧美成人一区二区免费高清观看 | 午夜福利18| 可以在线观看的亚洲视频| 欧美成人午夜精品| 午夜免费激情av| 成人国产综合亚洲| av在线天堂中文字幕| 高潮久久久久久久久久久不卡| 一区福利在线观看| 成人亚洲精品av一区二区| 久久欧美精品欧美久久欧美| 女警被强在线播放| 日韩欧美 国产精品| 日韩大码丰满熟妇| 小说图片视频综合网站| 色精品久久人妻99蜜桃| 国产精品国产高清国产av| 中国国产av一级| 国内精品久久久久精免费| 日韩欧美一区二区三区在线观看| 成年av动漫网址| 国产 一区精品| 久久精品夜色国产| 大型黄色视频在线免费观看| 久久九九热精品免费| 国产精品1区2区在线观看.| 一级毛片久久久久久久久女| 在线播放无遮挡| 一边亲一边摸免费视频| 久久精品国产99精品国产亚洲性色| 变态另类成人亚洲欧美熟女| 哪个播放器可以免费观看大片| 网址你懂的国产日韩在线| 欧美性猛交╳xxx乱大交人| 国产v大片淫在线免费观看| 亚洲性久久影院| 美女黄网站色视频| 蜜桃亚洲精品一区二区三区| 夜夜看夜夜爽夜夜摸| 有码 亚洲区| 精品日产1卡2卡| 日本免费一区二区三区高清不卡| 中文字幕熟女人妻在线| 亚洲天堂国产精品一区在线| 三级国产精品欧美在线观看| 99久久无色码亚洲精品果冻| av天堂中文字幕网| 男插女下体视频免费在线播放| 伦精品一区二区三区| 一级毛片电影观看 | 亚洲久久久久久中文字幕| 亚洲精品日韩av片在线观看| 99久国产av精品| 深夜a级毛片| 性欧美人与动物交配| 日本与韩国留学比较| 美女 人体艺术 gogo| 亚洲成人久久性| 国产片特级美女逼逼视频| av在线蜜桃| 亚洲最大成人手机在线| 1000部很黄的大片| 美女 人体艺术 gogo| 国产三级中文精品| 亚洲一区高清亚洲精品| 国产精品福利在线免费观看| 国产精品久久久久久久久免| 村上凉子中文字幕在线| 成熟少妇高潮喷水视频| 久久亚洲精品不卡| 老女人水多毛片| 日韩欧美三级三区| 日韩大尺度精品在线看网址| 欧美一区二区国产精品久久精品| 久久亚洲国产成人精品v| 亚洲最大成人av| 久久亚洲国产成人精品v| 国产精品女同一区二区软件| 亚洲无线观看免费| 狠狠狠狠99中文字幕| 国产成人freesex在线| 一个人免费在线观看电影| 国产精品99久久久久久久久| 亚洲国产精品成人综合色| 亚洲av二区三区四区| 卡戴珊不雅视频在线播放| 亚洲经典国产精华液单| 一边摸一边抽搐一进一小说| 三级国产精品欧美在线观看| 三级毛片av免费| 欧美又色又爽又黄视频| 美女高潮的动态| 免费观看a级毛片全部| 久久草成人影院| 97人妻精品一区二区三区麻豆| 欧美成人免费av一区二区三区| 男人舔女人下体高潮全视频| 午夜久久久久精精品| 黑人高潮一二区| 欧美最黄视频在线播放免费| 中国国产av一级| 欧美+日韩+精品| av在线播放精品| 久久久久久久久久久免费av| 日本一二三区视频观看| 日韩欧美 国产精品| 三级毛片av免费| 国产三级在线视频| 成年av动漫网址| 久久久久久久久大av| 99热这里只有是精品50| 国产一区二区激情短视频| 97超碰精品成人国产| 中文字幕熟女人妻在线| www.av在线官网国产| 亚洲精品影视一区二区三区av| 美女脱内裤让男人舔精品视频 | 内射极品少妇av片p| 亚洲欧美日韩卡通动漫| 国内精品一区二区在线观看| 18禁在线无遮挡免费观看视频| 夜夜看夜夜爽夜夜摸| 99在线视频只有这里精品首页| 国产精品一区二区在线观看99 | 99在线人妻在线中文字幕| 国产午夜精品一二区理论片| 观看美女的网站| 国产伦在线观看视频一区| 激情 狠狠 欧美| 免费看光身美女| 美女内射精品一级片tv| 国产午夜福利久久久久久| 亚洲成人久久爱视频| 日韩精品青青久久久久久| 成人二区视频| 久久久久久伊人网av| 国产色婷婷99| 国产一区二区三区在线臀色熟女| 国产毛片a区久久久久| 亚洲精品久久国产高清桃花| 国产精品综合久久久久久久免费| 啦啦啦啦在线视频资源| 久久久成人免费电影| 亚洲在线自拍视频| 黄色视频,在线免费观看| 欧美激情久久久久久爽电影| 99久久精品国产国产毛片| 国产乱人视频| 男人狂女人下面高潮的视频| 亚洲av男天堂| 午夜视频国产福利| av视频在线观看入口| 亚洲美女视频黄频| 直男gayav资源| 国产男人的电影天堂91| 国产成人午夜福利电影在线观看| 国产黄片视频在线免费观看| .国产精品久久| 日韩亚洲欧美综合| 中文欧美无线码| 久久精品国产亚洲av天美| 爱豆传媒免费全集在线观看| 97超碰精品成人国产| 日日摸夜夜添夜夜爱| 免费观看在线日韩| 少妇人妻精品综合一区二区 | 午夜激情欧美在线| 在线观看午夜福利视频| 国产精品综合久久久久久久免费| 久久草成人影院| 免费av毛片视频| av专区在线播放| 精品免费久久久久久久清纯| 亚洲国产高清在线一区二区三| 欧美xxxx性猛交bbbb| av国产免费在线观看| 国产亚洲av片在线观看秒播厂 | 日韩av不卡免费在线播放| 国产精品伦人一区二区| 国产蜜桃级精品一区二区三区| 99热精品在线国产| 国产视频首页在线观看| 亚洲欧美日韩东京热| 99riav亚洲国产免费| 日韩人妻高清精品专区| 欧美成人精品欧美一级黄| 亚洲av.av天堂| 青春草视频在线免费观看| 亚洲av一区综合| 欧美区成人在线视频| 久久久久网色| 亚洲成人中文字幕在线播放| 亚洲在线自拍视频| 日本免费a在线| 国产熟女欧美一区二区| 国产老妇伦熟女老妇高清| 蜜臀久久99精品久久宅男| 内地一区二区视频在线| 国产一区二区激情短视频| 少妇熟女aⅴ在线视频| 给我免费播放毛片高清在线观看| 中文欧美无线码| 午夜久久久久精精品| 国产精品爽爽va在线观看网站| 欧美不卡视频在线免费观看| 1024手机看黄色片| 国产一区二区三区av在线 | 久久久久久国产a免费观看| 成人高潮视频无遮挡免费网站| 日本一二三区视频观看| 亚洲乱码一区二区免费版| 午夜久久久久精精品| 淫秽高清视频在线观看| av天堂中文字幕网| 少妇熟女欧美另类| 欧美日韩精品成人综合77777| 亚洲在线观看片| 综合色丁香网| 国产一级毛片七仙女欲春2| 麻豆av噜噜一区二区三区| 国产精品,欧美在线| 人人妻人人看人人澡| 嫩草影院入口| 神马国产精品三级电影在线观看| 精品久久久久久久久久久久久| 国产一区二区激情短视频| 婷婷亚洲欧美| 简卡轻食公司| 热99在线观看视频| 亚洲人成网站在线观看播放| 午夜福利在线观看吧| 波多野结衣高清无吗| 国产乱人偷精品视频| 亚洲最大成人av| 在线观看免费视频日本深夜| 国产亚洲欧美98| 亚洲av男天堂| 成人永久免费在线观看视频| 亚洲四区av| 最好的美女福利视频网| 99热这里只有是精品50| а√天堂www在线а√下载| 亚洲一区二区三区色噜噜| 亚洲欧美精品综合久久99| 国产av不卡久久| 国产精品乱码一区二三区的特点| 国产片特级美女逼逼视频| 国产不卡一卡二| 国产激情偷乱视频一区二区| 18禁在线无遮挡免费观看视频| 神马国产精品三级电影在线观看| 中文字幕人妻熟人妻熟丝袜美| 美女xxoo啪啪120秒动态图| 亚洲国产高清在线一区二区三| 干丝袜人妻中文字幕| 久久久久久久久久黄片| 天天一区二区日本电影三级| 国产精品一区www在线观看| 午夜a级毛片| 看黄色毛片网站| 性色avwww在线观看| 中文欧美无线码| 午夜精品一区二区三区免费看| 少妇裸体淫交视频免费看高清| 身体一侧抽搐| 内地一区二区视频在线| 久久精品综合一区二区三区| 国产老妇女一区| 久久久国产成人精品二区| 国产男人的电影天堂91| 免费无遮挡裸体视频| 国产中年淑女户外野战色| 国产精品,欧美在线| 伊人久久精品亚洲午夜| 99在线视频只有这里精品首页| 国内揄拍国产精品人妻在线| 成人午夜精彩视频在线观看| 性色avwww在线观看| 精品人妻视频免费看| 亚洲丝袜综合中文字幕| 一区二区三区免费毛片| 中文欧美无线码| 啦啦啦韩国在线观看视频| 卡戴珊不雅视频在线播放| 国产极品天堂在线| 天堂网av新在线| 夜夜爽天天搞| 午夜福利在线在线| 欧美一级a爱片免费观看看| 色播亚洲综合网| 亚洲av中文av极速乱| 国产高清三级在线| 欧美在线一区亚洲| 亚洲精品自拍成人| 国产精品麻豆人妻色哟哟久久 | 免费观看人在逋| 欧美一区二区精品小视频在线| 亚洲国产精品国产精品| 免费大片18禁| 国产爱豆传媒在线观看| 国产大屁股一区二区在线视频| а√天堂www在线а√下载| 又粗又爽又猛毛片免费看| 午夜久久久久精精品| 狠狠狠狠99中文字幕| 亚洲精品日韩在线中文字幕 | 免费大片18禁| 国内少妇人妻偷人精品xxx网站| 高清在线视频一区二区三区 | 亚洲欧美日韩高清专用| 九九久久精品国产亚洲av麻豆| 欧美激情国产日韩精品一区| 国产成人aa在线观看| 免费不卡的大黄色大毛片视频在线观看 | 久久精品久久久久久久性|