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

    基于混合粒子群算法的煤層氣井位優(yōu)化方法

    2020-11-09 03:07:32劉明明田中政
    巖性油氣藏 2020年6期
    關(guān)鍵詞:井位模擬退火現(xiàn)值

    劉明明,王 全,馬 收,田中政,叢 顏

    (1.華美孚泰油氣增產(chǎn)技術(shù)服務(wù)有限責(zé)任公司,北京 100101;2.石化盈科信息技術(shù)有限責(zé)任公司,北京 100007)

    0 引言

    煤層氣作為非常規(guī)油氣的重要組成部分,具有重要的研究?jī)r(jià)值[1-3]。煤層氣開(kāi)發(fā)是集多種技術(shù)于一體的復(fù)雜系統(tǒng)工程,而井位選擇是煤層氣勘探開(kāi)發(fā)的基礎(chǔ),涉及到地質(zhì)構(gòu)造條件、煤層發(fā)育特征、煤層氣物性、施工工程條件等一系列因素,是國(guó)內(nèi)外研究的熱點(diǎn)[4-5]。國(guó)內(nèi)研究人員多依據(jù)油藏地質(zhì)情況制定若干井網(wǎng)方案,通過(guò)油藏?cái)?shù)值模擬進(jìn)行井位優(yōu)化研究,此方法對(duì)于研究人員的經(jīng)驗(yàn)要求較高,但油藏地質(zhì)由于不確定性因素過(guò)多,此方法不易得到最優(yōu)解[6-9]。

    國(guó)外研究人員主要通過(guò)優(yōu)化算法來(lái)篩選最優(yōu)解,首先選定一個(gè)目標(biāo)函數(shù)(凈現(xiàn)值或采收率等),然后以算法為手段,調(diào)用數(shù)值模擬軟件優(yōu)化方案,迭代循環(huán)得到最佳方案。目前,廣泛應(yīng)用于井位優(yōu)化的無(wú)梯度優(yōu)化算法多是一些全局優(yōu)化算法,比如遺傳算法(GA)、粒子群算法(PSO)、模擬退火算法(SA)、蟻群算法(ACO)等。其中,模擬退火算法、粒子群算法等在井位優(yōu)化中的應(yīng)用越來(lái)越廣泛。1995 年Beckner 等[10]將模擬退火算法引入到井位優(yōu)化研究中;2002 年Norrena 等[11]采用模擬退火算法并考慮經(jīng)濟(jì)約束進(jìn)行井位優(yōu)化研究;2010 年Onwunalu等[12]將粒子群算法引入到井位優(yōu)化研究中;2011 年Onwunalu 等[13]采用粒子群算法進(jìn)行規(guī)則井網(wǎng)優(yōu)化研究。國(guó)內(nèi)姜瑞忠等[14-15]分別于2014 年和2018 年提出用遺傳算法和新型遺傳算法進(jìn)行井位優(yōu)化研究。

    本次煤層氣井位優(yōu)化研究以沁水盆地沁端區(qū)塊煤層氣開(kāi)發(fā)為例,基于粒子群算法的劣勢(shì)和模擬退火算法的優(yōu)勢(shì),提出一種混合粒子群算法,其以?xún)衄F(xiàn)值為目標(biāo)函數(shù),單井控制面積和井位為變量,采用優(yōu)化理論結(jié)合油藏?cái)?shù)值模擬的方法,優(yōu)選出凈現(xiàn)值最大的單井控制面積和井位,并將優(yōu)化結(jié)果與常規(guī)矩形布井方式進(jìn)行對(duì)比,以期為煤層氣開(kāi)發(fā)井位優(yōu)化提供新方法。

    1 混合粒子群算法

    粒子群算法是由Kennedy 等[16]于1995 年提出的一種智能優(yōu)化算法。它類(lèi)似于遺傳算法,也是從隨機(jī)解出發(fā),通過(guò)適應(yīng)度評(píng)價(jià)解的優(yōu)劣,迭代尋找最優(yōu)解。其優(yōu)點(diǎn)是概念比較清晰,過(guò)程簡(jiǎn)單,涉及的計(jì)算參數(shù)比較少,所需粒子數(shù)不多,比較容易實(shí)現(xiàn),但是它也有缺陷,比如容易早熟,收斂速度慢,容易陷入局部最優(yōu)等。對(duì)粒子群算法的改進(jìn)大致歸為2 類(lèi):一類(lèi)是在基本的算法基礎(chǔ)上通過(guò)增加某些策略改進(jìn)其基本算法的性能;另一類(lèi)則是將粒子群算法與其他的優(yōu)化算法混合而形成一種混合算法?;旌线\(yùn)用多種智能算法可以取長(zhǎng)補(bǔ)短,加快收斂速度的同時(shí)有效地避免早熟問(wèn)題。

    模擬退火算法的思想最早是由Metropolis等[17]于1953 年提出來(lái)的,1983 年Kirkpatrick 等[18]將其應(yīng)用于組合優(yōu)化。模擬退火算法是基于蒙特卡羅迭代求解策略的一種隨機(jī)尋優(yōu)算法,其出發(fā)點(diǎn)是基于固體物質(zhì)的退火過(guò)程與組合優(yōu)化問(wèn)題之間的相似性,它具有跳出局部最優(yōu)和搜索精度高的優(yōu)點(diǎn)。

    基于此,筆者提出一種混合粒子群算法,該算法結(jié)合粒子群算法和模擬退火算法,將模擬退火思想引入粒子群算法中,增強(qiáng)了粒子群算法跳出局部最優(yōu)的能力[19-21],并通過(guò)matlab 編程來(lái)實(shí)現(xiàn)?;旌狭W尤核惴ò凑漳M退火接受準(zhǔn)則允許目標(biāo)函數(shù)依概率在有限范圍內(nèi)變壞[22]。

    本文將混合粒子群算法運(yùn)用到井位優(yōu)化過(guò)程的主要思路為:

    2 地質(zhì)模型建立

    2.1 區(qū)塊簡(jiǎn)介

    沁端區(qū)塊位于沁水盆地南部,隸屬于山西省沁水縣,區(qū)內(nèi)由老到新依次發(fā)育奧陶系、石炭系、二疊系、三疊系和第四系。本區(qū)構(gòu)造形態(tài)總體為一走向北北東、傾向北西西的單斜構(gòu)造,在此基礎(chǔ)上發(fā)育了一系列近南北—北北東向的寬緩褶曲,形成區(qū)內(nèi)地層的波狀起伏,巖層傾角一般不超過(guò)15°,個(gè)別地段受構(gòu)造影響地層傾角變化大,斷層不發(fā)育。

    區(qū)塊內(nèi)主要含煤地層為上石炭統(tǒng)太原組(C3t)和下二疊統(tǒng)山西組(P1s)。含煤17 層,煤層總厚度為14.67 m,其中主要可采煤層2 層,分別為3# 煤層和15#煤層,3#煤層厚度為6.05~6.60 m,平均為6.24 m;15#煤層厚度為3.8~4.5 m,平均為4.14 m。開(kāi)發(fā)方式多采用直井壓裂合采,單井穩(wěn)定日產(chǎn)氣為1 000~2 000 m3。

    2.2 模型建立

    基于Petrel 地質(zhì)建模軟件,利用生產(chǎn)井測(cè)井資料,建立區(qū)塊的地質(zhì)模型(圖1),縱向上分為3 個(gè)小層(包括隔層),x,y方向上的網(wǎng)格步長(zhǎng)為20 m,模型的總結(jié)點(diǎn)數(shù)為241 029(圖2)。根據(jù)測(cè)井解釋的成果,使用克里金插值方法建立區(qū)塊的煤層氣含量、孔隙度、滲透率和凈毛比等屬性模型。

    圖1 沁端區(qū)塊滲透率分布及試驗(yàn)區(qū)位置Fig.1 Permeability distribution and pilot site location in Qinduan block

    本區(qū)塊煤層氣藏整體面積大,若選取整個(gè)氣藏進(jìn)行井位優(yōu)化的數(shù)模工作,計(jì)算量大,時(shí)間長(zhǎng)。因此,依據(jù)各層的物性分布和儲(chǔ)量豐度特征,選取1.2 km×1.2 km 的試驗(yàn)區(qū)進(jìn)行井位優(yōu)化數(shù)值模擬工作。模型網(wǎng)格劃分為61×61×3,x,y方向上的網(wǎng)格步長(zhǎng)均為20 m,所選模型原始地質(zhì)儲(chǔ)量為3.52 億m3。生產(chǎn)制度為直井壓裂,定壓生產(chǎn)(0.5 MPa),總共生產(chǎn)5 000 d。

    3 基于混合粒子群算法的井位優(yōu)化

    3.1 井位優(yōu)化目標(biāo)函數(shù)(凈現(xiàn)值)

    煤層氣開(kāi)發(fā)的井網(wǎng)部署主要考慮3 個(gè)方面的因素:地質(zhì)因素、開(kāi)發(fā)因素與經(jīng)濟(jì)效益。對(duì)于滲透率、煤層氣含量、孔隙度等地質(zhì)因素以及井間干擾、開(kāi)發(fā)層系、采氣速度等開(kāi)發(fā)因素,可以通過(guò)油藏?cái)?shù)值模擬進(jìn)行研究,而經(jīng)濟(jì)效益決定該煤層氣區(qū)塊能否進(jìn)行商業(yè)開(kāi)發(fā)或者采取增產(chǎn)措施,這就要求合理的井位不僅須要滿(mǎn)足采氣速度的要求,更要在經(jīng)濟(jì)上可行。因此,本次井位優(yōu)化研究選用凈現(xiàn)值為目標(biāo)函數(shù)[23]。

    凈現(xiàn)值是指在一個(gè)項(xiàng)目的整個(gè)建設(shè)和生產(chǎn)服務(wù)年限內(nèi)各時(shí)間段的凈現(xiàn)金流量按照設(shè)定的折現(xiàn)率折成現(xiàn)值后求和所得到的值。凈現(xiàn)值越大,項(xiàng)目收益越大。其表達(dá)式為

    式中:NPV為凈現(xiàn)值,元;t0為項(xiàng)目的開(kāi)始年份;tf為項(xiàng)目的結(jié)束年份;Cin(t)為t時(shí)刻的現(xiàn)金流入量,主要是銷(xiāo)售收入,元;Cout(t)為t時(shí)刻的現(xiàn)金流出量,主要是投資、成本及銷(xiāo)售稅金,元;r為折現(xiàn)率,指將未來(lái)有限期預(yù)期收益折算成現(xiàn)值的比率;Qg(t)為第t年的累計(jì)產(chǎn)氣量,m3;Qw(t)為第t年的累計(jì)產(chǎn)水量,m3;Gas Price為煤層氣價(jià)格,元/m3;Water OPEX為污水處理價(jià)格,元/m3;OPEX為操作成本,主要包括人員工資、水電費(fèi)、運(yùn)輸費(fèi)用、維護(hù)費(fèi)、設(shè)備更新費(fèi)等,元;TAX為稅費(fèi),元;CAPEX為基建費(fèi)用,主要包括鉆井、完井、固井、射孔、壓裂、酸化等措施費(fèi)用以及所需物資的費(fèi)用,元。

    油田開(kāi)發(fā)中采用凈現(xiàn)值方法計(jì)算利潤(rùn)時(shí),須要將開(kāi)發(fā)時(shí)間劃分為若干個(gè)時(shí)間間隔,因此,凈現(xiàn)值的目標(biāo)函數(shù)變?yōu)殡x散型

    計(jì)算凈現(xiàn)值用的參數(shù)如表1 所列。

    表1 經(jīng)濟(jì)參數(shù)Table 1 Economic parameters for net present value computation

    3.2 變量及約束條件

    井位坐標(biāo)X(I,J)為井位優(yōu)化研究的變量,而且在地質(zhì)模型中,井位位于網(wǎng)格內(nèi)。假定井?dāng)?shù)為Nw,則問(wèn)題的維數(shù)為D=2Nw。井位坐標(biāo)X如下:

    由于網(wǎng)格在x和y方向是有限個(gè)的,因此,約束條件為:

    (1)1

    (2)1

    3.3 井位優(yōu)化流程

    基于建立的地質(zhì)模型及目標(biāo)函數(shù),將混合粒子群應(yīng)用到井位優(yōu)化中?;诨旌狭W尤核惴ǖ木粌?yōu)化流程如圖3 所示。

    圖3 基于混合粒子群算法的井位優(yōu)化流程圖Fig.3 Well location optimization flowchart based on hybrid particle swarm optimization algorithm

    井位優(yōu)化過(guò)程中混合粒子群算法中慣性權(quán)重ω、最大迭代次數(shù)kmax、學(xué)習(xí)因子c1/c2、粒子速度V、允許變壞的范圍Fval和退溫系數(shù)α等基本參數(shù)取值如表2 所列。

    表2 PSOSA 基本參數(shù)Table 2 Basic parameters for hybrid particle swarm optimization algorithm

    4 井位優(yōu)化方法應(yīng)用

    4.1 可行性研究

    在本次可行性研究中,通過(guò)研究一口生產(chǎn)井的最佳井位來(lái)證明本方法的優(yōu)勢(shì)。針對(duì)本模型,共有61×61=3 721 個(gè)位置,不同位置的含氣量與滲透率不同。因此,通過(guò)3 721 次窮舉計(jì)算可以獲得最佳井位W1(50,52)(參見(jiàn)圖2)與全區(qū)的凈現(xiàn)值分布(圖4)。

    圖4 全區(qū)凈現(xiàn)值分布Fig.4 Net present value distribution

    對(duì)比圖2(a)與圖2(b)可以發(fā)現(xiàn),與含氣量分布相比,滲透率分布對(duì)煤層氣開(kāi)發(fā)的效果影響更大。分析其原因?yàn)闈B透率決定井的采氣速度,而含氣量決定井的生產(chǎn)時(shí)間,只有當(dāng)采氣速度和生產(chǎn)時(shí)間達(dá)到一個(gè)最優(yōu)組合時(shí)才能使單井效益最大化。

    分析圖4 可知,凈現(xiàn)值面有若干個(gè)極大值點(diǎn),這是由于儲(chǔ)層的非均質(zhì)性引起的,最優(yōu)井位為W1(50,52)。圖5 為利用混合粒子群算法進(jìn)行井位優(yōu)化研究的凈現(xiàn)值圖,分析發(fā)現(xiàn)混合粒子群算法可以快速收斂(Ns k=270次)到最優(yōu)井位W1(50,52),獲得最大凈現(xiàn)值。

    圖5 混合粒子群算法井位優(yōu)化方法凈現(xiàn)值Fig.5 Net present value of hybrid particle swarm optimization algorithm

    可行性研究結(jié)果表明,基于混合粒子群算法的井位優(yōu)化算法較窮舉法計(jì)算量大幅度降低,可以準(zhǔn)確確定最優(yōu)井位,用于實(shí)際井位優(yōu)化。

    4.2 單井控制面積對(duì)開(kāi)發(fā)效果的影響

    單井控制面積是一口井所控制的開(kāi)發(fā)面積的大?。?4]。煤層氣主要采用矩形、菱形井網(wǎng)等,單井控制面積為0.09~0.64 km2。單井控制面積涉及氣田開(kāi)發(fā)指標(biāo)和經(jīng)濟(jì)效益的評(píng)價(jià),是煤層氣氣田開(kāi)發(fā)的重要參數(shù),其大小與井型和井間距大小有關(guān),它取決于儲(chǔ)層的性質(zhì)以及煤層氣開(kāi)發(fā)的規(guī)模。

    針對(duì)試驗(yàn)區(qū),利用混合粒子群井位優(yōu)化算法,用數(shù)值模擬技術(shù)研究不同的單井控制面積(0.06,0.08,0.10,0.20,0.30,0.40,0.50,0.80,1.44 km2)對(duì)開(kāi)發(fā)效果的影響。圖6 為不同單井控制面積下的最優(yōu)井位分布圖,分析發(fā)現(xiàn)不同單井控制面積下生產(chǎn)井?dāng)?shù)不同,最優(yōu)井位會(huì)發(fā)生變化,且位于滲透率高的位置,無(wú)明顯井網(wǎng)形式。

    圖6 不同單井控制面積下的最優(yōu)井位Fig.6 Optimum well location distribution under different single well control area

    圖7 為不同單井控制面積下的生產(chǎn)指標(biāo)曲線(xiàn),單井控制面積越大,區(qū)塊的凈現(xiàn)值越小,不存在明顯的拐點(diǎn);單井控制面積對(duì)氣井單井凈現(xiàn)值影響較大,存在明顯的拐點(diǎn)。當(dāng)單井控制面積大于0.2 km2時(shí),單井凈現(xiàn)值隨其增大而減小,原因是煤層氣井通過(guò)排水降壓達(dá)到解吸開(kāi)采,單井控制面積越大,導(dǎo)致井間干擾越小,壓力降低幅度越小,不利于煤層氣降壓開(kāi)采。因此,為獲得最大的區(qū)塊開(kāi)發(fā)收益及單井效益最大化,選擇最佳單井控制面積為0.2 km2。

    4.3 最優(yōu)單井控制面積下的井位優(yōu)選

    采用單井控制面積為0.2 km2進(jìn)行沁水盆地煤層氣田井位優(yōu)化設(shè)計(jì),當(dāng)種群數(shù)量設(shè)置為30、迭代次數(shù)設(shè)置為100 時(shí),混合粒子群算法得到的凈現(xiàn)值為34 115 452 元,大于常規(guī)矩形井網(wǎng)的凈現(xiàn)值30 312 502元,凈現(xiàn)值增加幅度為12.55%。

    圖7 生產(chǎn)指標(biāo)與單井控制面積的關(guān)系Fig.7 Relationship between single well control area and net present value

    從圖8 還可以看出混合粒子群算法尋優(yōu)能力較強(qiáng),能夠快速定位最佳井位。

    圖9 為混合粒子群算法和常規(guī)矩形布井的井位分布圖。從圖9 可以看出,最優(yōu)井位與含氣量及滲透率分布密切相關(guān),與含氣量分布相比,滲透率分布對(duì)煤層氣井位選取的影響較大。最優(yōu)井位多分布在滲透率高的部位,而滲透率低的部位少;含氣量大的部位布井較少,原因是此處滲透率低,無(wú)法滿(mǎn)足一定的采氣速度。最優(yōu)井位是含氣量與滲透率的最優(yōu)組合,此處生產(chǎn)井具有豐富的物質(zhì)基礎(chǔ)與合理的采氣速度。

    圖8 不同井位優(yōu)化方法凈現(xiàn)值對(duì)比Fig.8 Comparison of net present value between hybrid particle swarm optimization algorithm and conventional well pattern

    對(duì)比圖9(a)與圖9(b)還可以發(fā)現(xiàn),混合粒子群算法優(yōu)選出的井位與常規(guī)矩形井網(wǎng)類(lèi)似,分布比較均勻,這與本區(qū)塊非均質(zhì)程度較小有關(guān)。

    圖9 區(qū)塊屬性與最優(yōu)井位分布Fig.9 Geologic property and the optimal well location distribution

    5 結(jié)論

    (1)井位優(yōu)化方法可行性研究驗(yàn)證了基于混合粒子群算法的井位優(yōu)化方法較窮舉法有較大優(yōu)勢(shì),該方法能夠快速收斂到最優(yōu)井位,減少了計(jì)算量。

    (2)基于混合粒子群算法的井位優(yōu)化方法克服了常規(guī)井網(wǎng)部署的經(jīng)驗(yàn)依賴(lài)性,能夠確定儲(chǔ)層非均質(zhì)條件下的最優(yōu)井位,最大限度地提高了煤層氣田的經(jīng)濟(jì)效益。

    (3)選取不同的單井控制面積(0.06~1.44 km2)進(jìn)行沁水盆地煤層氣井位優(yōu)化研究,得到最優(yōu)單井控制面積為0.2 km2,對(duì)于最優(yōu)單井控制面積,混合粒子群算法得到的凈現(xiàn)值比常規(guī)矩形井網(wǎng)凈現(xiàn)值增加12.55%。

    (4)最優(yōu)井位分布與含氣量、滲透率均密切相關(guān),其中滲透率影響尤為重要,最優(yōu)井位是含氣量與滲透率的最優(yōu)組合,既要有豐富的物質(zhì)基礎(chǔ),又要有合理的采氣速度。

    猜你喜歡
    井位模擬退火現(xiàn)值
    基于改進(jìn)Nelder-Mead 算法的井位優(yōu)選方法
    基于Microstation二次開(kāi)發(fā)在油田井位展繪的應(yīng)用
    河北省行唐縣貧水山區(qū)物探找水研究
    渤海油田側(cè)鉆井槽口優(yōu)選設(shè)計(jì)方法
    海洋石油(2021年1期)2021-03-30 13:28:06
    模擬退火遺傳算法在機(jī)械臂路徑規(guī)劃中的應(yīng)用
    基于模糊自適應(yīng)模擬退火遺傳算法的配電網(wǎng)故障定位
    資金時(shí)間價(jià)值中的系數(shù)關(guān)系探析
    資金時(shí)間價(jià)值基礎(chǔ)運(yùn)算解讀
    SOA結(jié)合模擬退火算法優(yōu)化電容器配置研究
    凈現(xiàn)值法對(duì)比煤層氣與常規(guī)天然氣經(jīng)濟(jì)效益
    精品亚洲成国产av| 国产精品国产三级专区第一集| 成年动漫av网址| 国产精品熟女久久久久浪| 校园人妻丝袜中文字幕| 日本猛色少妇xxxxx猛交久久| 日韩伦理黄色片| av网站免费在线观看视频| 丝袜脚勾引网站| 在线观看国产h片| 男女国产视频网站| 成人国语在线视频| 777久久人妻少妇嫩草av网站| 1024香蕉在线观看| 国产老妇伦熟女老妇高清| 欧美日韩成人在线一区二区| 成年动漫av网址| 国产片特级美女逼逼视频| 两个人看的免费小视频| 免费观看a级毛片全部| 天天添夜夜摸| 久久久久久久大尺度免费视频| 一本—道久久a久久精品蜜桃钙片| 国产成人精品久久二区二区91 | 久久女婷五月综合色啪小说| 国产乱人偷精品视频| 香蕉国产在线看| 成人午夜精彩视频在线观看| 亚洲精品中文字幕在线视频| 久久久久人妻精品一区果冻| 制服丝袜香蕉在线| 人妻 亚洲 视频| 国产片特级美女逼逼视频| 综合色丁香网| 又大又黄又爽视频免费| 精品国产乱码久久久久久男人| 亚洲视频免费观看视频| 亚洲av在线观看美女高潮| 色播在线永久视频| 午夜老司机福利片| 国产精品一国产av| 日韩视频在线欧美| 亚洲精品成人av观看孕妇| 一区二区三区四区激情视频| 午夜老司机福利片| 美女大奶头黄色视频| 夫妻性生交免费视频一级片| 国产精品二区激情视频| 亚洲精品成人av观看孕妇| 久久韩国三级中文字幕| av福利片在线| 一级,二级,三级黄色视频| 午夜福利在线免费观看网站| 看十八女毛片水多多多| 热99久久久久精品小说推荐| 中文字幕最新亚洲高清| 亚洲欧洲日产国产| 91老司机精品| av卡一久久| 亚洲欧美精品综合一区二区三区| 高清在线视频一区二区三区| 国产伦理片在线播放av一区| 亚洲成色77777| 欧美亚洲 丝袜 人妻 在线| 日韩大片免费观看网站| 久久ye,这里只有精品| 精品国产国语对白av| 女性被躁到高潮视频| 悠悠久久av| 日韩电影二区| 性高湖久久久久久久久免费观看| 天天添夜夜摸| 纵有疾风起免费观看全集完整版| 国产极品天堂在线| 欧美人与性动交α欧美精品济南到| 别揉我奶头~嗯~啊~动态视频 | 热re99久久国产66热| 亚洲精品第二区| 十八禁高潮呻吟视频| 亚洲,欧美精品.| 亚洲av日韩精品久久久久久密 | 大陆偷拍与自拍| 久久人妻熟女aⅴ| 如日韩欧美国产精品一区二区三区| 国产成人免费观看mmmm| 亚洲欧美成人精品一区二区| 欧美日韩亚洲高清精品| 国产精品久久久久久久久免| videos熟女内射| 高清视频免费观看一区二区| 亚洲av日韩在线播放| 国产黄频视频在线观看| 日本vs欧美在线观看视频| 又粗又硬又长又爽又黄的视频| 欧美日韩亚洲国产一区二区在线观看 | 午夜免费男女啪啪视频观看| 女人高潮潮喷娇喘18禁视频| 日韩制服丝袜自拍偷拍| 亚洲伊人色综图| 亚洲国产成人一精品久久久| 亚洲精品国产av蜜桃| 亚洲国产精品国产精品| 午夜福利乱码中文字幕| 国产精品二区激情视频| 欧美日韩亚洲综合一区二区三区_| 午夜福利免费观看在线| 最近最新中文字幕大全免费视频 | 51午夜福利影视在线观看| 国产黄色免费在线视频| 国产精品亚洲av一区麻豆 | 晚上一个人看的免费电影| 女人被躁到高潮嗷嗷叫费观| 巨乳人妻的诱惑在线观看| 亚洲男人天堂网一区| 精品少妇黑人巨大在线播放| 国产av精品麻豆| 午夜福利在线免费观看网站| 国产男女超爽视频在线观看| 一本一本久久a久久精品综合妖精| 久久鲁丝午夜福利片| 亚洲七黄色美女视频| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇人妻 视频| 国产亚洲欧美精品永久| 久久午夜综合久久蜜桃| 9191精品国产免费久久| 97人妻天天添夜夜摸| 国产爽快片一区二区三区| avwww免费| 美女高潮到喷水免费观看| 亚洲成av片中文字幕在线观看| 男男h啪啪无遮挡| 日韩欧美精品免费久久| 亚洲精品在线美女| 免费黄频网站在线观看国产| 国产精品国产av在线观看| 夫妻性生交免费视频一级片| 男的添女的下面高潮视频| 天天操日日干夜夜撸| 欧美xxⅹ黑人| 国产伦人伦偷精品视频| 亚洲精品日韩在线中文字幕| 在线观看免费日韩欧美大片| 观看av在线不卡| 国产高清国产精品国产三级| 亚洲av电影在线观看一区二区三区| 色视频在线一区二区三区| 精品国产国语对白av| 少妇 在线观看| 欧美日韩亚洲高清精品| 亚洲国产毛片av蜜桃av| 1024香蕉在线观看| av.在线天堂| 国产一区二区在线观看av| 又粗又硬又长又爽又黄的视频| 在现免费观看毛片| 91精品三级在线观看| 在线观看免费午夜福利视频| 欧美 亚洲 国产 日韩一| 男人添女人高潮全过程视频| 国产一区二区三区av在线| 卡戴珊不雅视频在线播放| 免费黄频网站在线观看国产| 大香蕉久久成人网| 亚洲欧美一区二区三区久久| 国产又色又爽无遮挡免| 精品国产乱码久久久久久小说| 午夜免费鲁丝| 在线观看人妻少妇| 十八禁人妻一区二区| 久久青草综合色| 熟妇人妻不卡中文字幕| 一级毛片黄色毛片免费观看视频| 欧美少妇被猛烈插入视频| 啦啦啦啦在线视频资源| 黑丝袜美女国产一区| 爱豆传媒免费全集在线观看| 美女视频免费永久观看网站| 制服丝袜香蕉在线| xxxhd国产人妻xxx| 天天躁狠狠躁夜夜躁狠狠躁| 在线观看免费日韩欧美大片| 国产精品久久久久成人av| 国产淫语在线视频| 久久久亚洲精品成人影院| 日韩av免费高清视频| 亚洲欧美激情在线| 最近中文字幕2019免费版| xxx大片免费视频| 日韩一区二区三区影片| 亚洲av日韩精品久久久久久密 | 国产精品无大码| 丁香六月欧美| 欧美黑人精品巨大| 精品一区二区免费观看| 亚洲精品成人av观看孕妇| 在线天堂中文资源库| 免费不卡黄色视频| www日本在线高清视频| 热99久久久久精品小说推荐| 最近最新中文字幕免费大全7| 在线观看免费高清a一片| 国产亚洲精品第一综合不卡| 免费高清在线观看视频在线观看| 99九九在线精品视频| 嫩草影视91久久| 制服丝袜香蕉在线| 国精品久久久久久国模美| 天天躁狠狠躁夜夜躁狠狠躁| 婷婷色综合大香蕉| 成年人午夜在线观看视频| 午夜免费男女啪啪视频观看| 一本色道久久久久久精品综合| 国产精品免费视频内射| 国产av码专区亚洲av| 精品一区在线观看国产| 在线观看人妻少妇| 国产黄频视频在线观看| 午夜91福利影院| 在线观看免费视频网站a站| 中文欧美无线码| 午夜免费鲁丝| 老汉色∧v一级毛片| 亚洲图色成人| 久久影院123| 亚洲国产最新在线播放| 亚洲国产欧美日韩在线播放| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av高清一级| 一本—道久久a久久精品蜜桃钙片| 人妻人人澡人人爽人人| 国产深夜福利视频在线观看| 精品亚洲成国产av| 看十八女毛片水多多多| 观看av在线不卡| 中文字幕高清在线视频| 国产日韩欧美亚洲二区| 欧美最新免费一区二区三区| 日韩精品有码人妻一区| 亚洲av电影在线观看一区二区三区| 国产在视频线精品| 韩国高清视频一区二区三区| 午夜福利一区二区在线看| 丝袜美腿诱惑在线| 久久久久精品久久久久真实原创| 美女福利国产在线| 午夜精品国产一区二区电影| 黄片小视频在线播放| 国产精品成人在线| 亚洲av电影在线观看一区二区三区| 婷婷色综合www| 18禁国产床啪视频网站| 免费高清在线观看视频在线观看| 在线观看免费高清a一片| 性色av一级| 欧美激情极品国产一区二区三区| 欧美精品一区二区免费开放| 看免费av毛片| 久久人人97超碰香蕉20202| 最近中文字幕2019免费版| 丝袜人妻中文字幕| 在线观看免费高清a一片| 日韩 欧美 亚洲 中文字幕| 国产成人精品在线电影| 欧美日韩成人在线一区二区| www.精华液| 久热爱精品视频在线9| 两个人免费观看高清视频| 午夜福利网站1000一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 国产精品蜜桃在线观看| 1024香蕉在线观看| 欧美精品一区二区免费开放| 天堂8中文在线网| 黄网站色视频无遮挡免费观看| 久久婷婷青草| 亚洲一区二区三区欧美精品| 久久热在线av| 蜜桃国产av成人99| 亚洲av日韩精品久久久久久密 | 亚洲国产欧美一区二区综合| 男女边吃奶边做爰视频| 超碰97精品在线观看| 亚洲久久久国产精品| 国产国语露脸激情在线看| 在线观看国产h片| 欧美另类一区| 777久久人妻少妇嫩草av网站| 国产探花极品一区二区| 精品久久蜜臀av无| 青春草国产在线视频| 久久久久久久久久久免费av| 一边亲一边摸免费视频| 免费黄色在线免费观看| 亚洲国产日韩一区二区| 亚洲一级一片aⅴ在线观看| 曰老女人黄片| 久久精品熟女亚洲av麻豆精品| 又大又爽又粗| 久久韩国三级中文字幕| 另类亚洲欧美激情| 日韩制服骚丝袜av| 97人妻天天添夜夜摸| 高清视频免费观看一区二区| 国产一区二区在线观看av| 少妇人妻精品综合一区二区| 亚洲中文av在线| 国产一区二区在线观看av| 一本大道久久a久久精品| 人妻 亚洲 视频| 一区在线观看完整版| 亚洲国产精品一区三区| 国产精品久久久久久精品古装| 久久鲁丝午夜福利片| 成人黄色视频免费在线看| 91精品国产国语对白视频| 精品亚洲成国产av| 香蕉丝袜av| 少妇人妻久久综合中文| 天天躁日日躁夜夜躁夜夜| 波多野结衣av一区二区av| 人体艺术视频欧美日本| 久久这里只有精品19| 欧美人与性动交α欧美软件| h视频一区二区三区| 大香蕉久久成人网| 久久韩国三级中文字幕| 新久久久久国产一级毛片| 国产伦人伦偷精品视频| 色综合欧美亚洲国产小说| av卡一久久| 美女扒开内裤让男人捅视频| a 毛片基地| 亚洲精品美女久久av网站| 亚洲欧美成人综合另类久久久| 中文乱码字字幕精品一区二区三区| 久久国产精品男人的天堂亚洲| 国产av码专区亚洲av| 人人妻人人添人人爽欧美一区卜| 久久人人97超碰香蕉20202| 亚洲欧美精品自产自拍| 欧美国产精品一级二级三级| 啦啦啦啦在线视频资源| 另类亚洲欧美激情| 哪个播放器可以免费观看大片| 人妻人人澡人人爽人人| 日韩欧美一区视频在线观看| 欧美精品一区二区免费开放| 黑人欧美特级aaaaaa片| 久久久久久人妻| 国产熟女欧美一区二区| 国产日韩一区二区三区精品不卡| 女人被躁到高潮嗷嗷叫费观| 欧美精品人与动牲交sv欧美| 亚洲精品av麻豆狂野| 久久久国产欧美日韩av| 一级,二级,三级黄色视频| 久久影院123| 曰老女人黄片| 人妻 亚洲 视频| 精品一区二区三区四区五区乱码 | 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品日本国产第一区| 国产一区二区三区综合在线观看| 日韩一卡2卡3卡4卡2021年| 久久人人爽人人片av| 国产精品.久久久| 熟女少妇亚洲综合色aaa.| 悠悠久久av| 国产女主播在线喷水免费视频网站| 亚洲情色 制服丝袜| av线在线观看网站| 国产成人精品无人区| 我的亚洲天堂| 免费黄频网站在线观看国产| 亚洲熟女精品中文字幕| 丰满少妇做爰视频| 免费看不卡的av| 国产伦人伦偷精品视频| 9色porny在线观看| 精品人妻一区二区三区麻豆| 男的添女的下面高潮视频| 啦啦啦中文免费视频观看日本| av在线app专区| 婷婷色麻豆天堂久久| 国产成人免费无遮挡视频| 国产成人系列免费观看| 天天操日日干夜夜撸| 丁香六月欧美| 国产精品麻豆人妻色哟哟久久| 亚洲人成电影观看| 中文字幕av电影在线播放| 国产日韩欧美在线精品| 大香蕉久久成人网| 中文精品一卡2卡3卡4更新| 国产乱来视频区| 欧美激情高清一区二区三区 | 人妻 亚洲 视频| 日日撸夜夜添| 另类精品久久| 高清黄色对白视频在线免费看| 日本爱情动作片www.在线观看| 亚洲国产毛片av蜜桃av| 亚洲av成人精品一二三区| 男人爽女人下面视频在线观看| 国产欧美亚洲国产| 性少妇av在线| 日韩不卡一区二区三区视频在线| 亚洲成人手机| 一本一本久久a久久精品综合妖精| 精品久久蜜臀av无| 亚洲精品国产av蜜桃| 欧美日韩亚洲综合一区二区三区_| 一级片'在线观看视频| 侵犯人妻中文字幕一二三四区| 成人免费观看视频高清| 男女边吃奶边做爰视频| 晚上一个人看的免费电影| 黄频高清免费视频| 精品第一国产精品| 三上悠亚av全集在线观看| 欧美日韩一区二区视频在线观看视频在线| √禁漫天堂资源中文www| 亚洲国产欧美日韩在线播放| 大陆偷拍与自拍| 欧美少妇被猛烈插入视频| 久久久久久人妻| 亚洲七黄色美女视频| 久久久久人妻精品一区果冻| 丝袜在线中文字幕| 人成视频在线观看免费观看| 久久热在线av| videos熟女内射| 亚洲精品自拍成人| 久久青草综合色| 亚洲第一av免费看| 欧美亚洲 丝袜 人妻 在线| 丁香六月欧美| 国产免费又黄又爽又色| 日本av免费视频播放| 国产 精品1| 性高湖久久久久久久久免费观看| 国产亚洲午夜精品一区二区久久| 欧美亚洲日本最大视频资源| 天堂中文最新版在线下载| 天天添夜夜摸| 男的添女的下面高潮视频| 久久久久久久久免费视频了| 晚上一个人看的免费电影| 这个男人来自地球电影免费观看 | 在线观看一区二区三区激情| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频| 国产成人午夜福利电影在线观看| 看非洲黑人一级黄片| 欧美精品亚洲一区二区| 国产高清不卡午夜福利| 色视频在线一区二区三区| 极品少妇高潮喷水抽搐| 中文天堂在线官网| 色播在线永久视频| 久久狼人影院| 美女主播在线视频| 黄色视频在线播放观看不卡| 亚洲av电影在线进入| 久久久久精品国产欧美久久久 | 只有这里有精品99| 中文字幕精品免费在线观看视频| 岛国毛片在线播放| 亚洲国产av影院在线观看| 国产精品一二三区在线看| 亚洲精品aⅴ在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲男人天堂网一区| 国产成人欧美在线观看 | 国产一区二区 视频在线| 亚洲国产毛片av蜜桃av| 亚洲熟女精品中文字幕| 国产 精品1| 久久久久久久久久久免费av| 欧美精品av麻豆av| 国产免费又黄又爽又色| 日本wwww免费看| 麻豆精品久久久久久蜜桃| 亚洲av日韩在线播放| 国精品久久久久久国模美| 精品人妻在线不人妻| 亚洲精品乱久久久久久| 男人舔女人的私密视频| 久久性视频一级片| 国产精品一区二区精品视频观看| 免费高清在线观看日韩| 一区二区三区四区激情视频| 国产亚洲av片在线观看秒播厂| 日本午夜av视频| 国产精品免费视频内射| 国产精品久久久久久人妻精品电影 | 亚洲欧洲精品一区二区精品久久久 | 巨乳人妻的诱惑在线观看| 日日撸夜夜添| 国产精品蜜桃在线观看| 国产精品无大码| 韩国高清视频一区二区三区| 亚洲av欧美aⅴ国产| 高清欧美精品videossex| 99精国产麻豆久久婷婷| 久久99热这里只频精品6学生| 亚洲男人天堂网一区| 中文精品一卡2卡3卡4更新| 十八禁高潮呻吟视频| 日韩制服骚丝袜av| 18禁观看日本| 99香蕉大伊视频| 妹子高潮喷水视频| 又大又爽又粗| 少妇 在线观看| 亚洲国产欧美一区二区综合| 不卡av一区二区三区| 90打野战视频偷拍视频| 嫩草影院入口| 女人久久www免费人成看片| av在线播放精品| 亚洲国产精品999| 自拍欧美九色日韩亚洲蝌蚪91| 熟妇人妻不卡中文字幕| 最黄视频免费看| 亚洲av综合色区一区| 欧美精品一区二区免费开放| 狂野欧美激情性bbbbbb| 一级,二级,三级黄色视频| 国产精品熟女久久久久浪| 成年美女黄网站色视频大全免费| 啦啦啦视频在线资源免费观看| 久久精品国产亚洲av涩爱| netflix在线观看网站| 岛国毛片在线播放| 亚洲 欧美一区二区三区| 免费在线观看完整版高清| 国产精品三级大全| 在线观看免费视频网站a站| 精品午夜福利在线看| 精品一区二区三卡| 一本色道久久久久久精品综合| tube8黄色片| 桃花免费在线播放| 天天躁夜夜躁狠狠躁躁| 男人操女人黄网站| 一本大道久久a久久精品| 国产成人精品久久久久久| 夜夜骑夜夜射夜夜干| 免费女性裸体啪啪无遮挡网站| 一区二区三区乱码不卡18| 免费女性裸体啪啪无遮挡网站| 夜夜骑夜夜射夜夜干| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕高清在线视频| 国产欧美亚洲国产| 校园人妻丝袜中文字幕| 中文字幕人妻丝袜一区二区 | 18禁国产床啪视频网站| 久久久久人妻精品一区果冻| 亚洲久久久国产精品| 精品第一国产精品| 亚洲精华国产精华液的使用体验| 国产视频首页在线观看| 亚洲在久久综合| 国精品久久久久久国模美| 免费高清在线观看日韩| 毛片一级片免费看久久久久| 色94色欧美一区二区| 一区二区三区四区激情视频| 九草在线视频观看| 高清不卡的av网站| 亚洲人成77777在线视频| 天天躁日日躁夜夜躁夜夜| 国产精品一区二区在线观看99| 无限看片的www在线观看| 中文字幕人妻丝袜制服| 好男人视频免费观看在线| 国产精品.久久久| 人妻人人澡人人爽人人| 天堂中文最新版在线下载| 欧美日韩精品网址| 1024视频免费在线观看| 精品国产超薄肉色丝袜足j| 亚洲人成电影观看| 国产一级毛片在线| 大香蕉久久成人网| 狂野欧美激情性bbbbbb| 最近中文字幕2019免费版| 精品第一国产精品| 日本欧美视频一区| 日本vs欧美在线观看视频| 国产精品久久久久久精品古装| 亚洲精品视频女| 五月天丁香电影| 国产av码专区亚洲av| 高清不卡的av网站| 中文字幕人妻丝袜一区二区 | 亚洲精品中文字幕在线视频| 在线天堂最新版资源| 嫩草影院入口| 不卡视频在线观看欧美| 亚洲欧洲精品一区二区精品久久久 | 国产免费一区二区三区四区乱码| 在线精品无人区一区二区三| 亚洲av电影在线观看一区二区三区| 日本欧美视频一区| 日韩中文字幕视频在线看片| 九九爱精品视频在线观看| 免费在线观看黄色视频的| 热re99久久精品国产66热6| 午夜日本视频在线| 女性被躁到高潮视频| 久久精品国产亚洲av高清一级|