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

    基于井間連通性的致密油藏注水井動(dòng)態(tài)裂縫研究

    2022-06-09 09:51:26曾慶橋張亮劉萍王哲王鏈趙良言李治平
    斷塊油氣田 2022年3期
    關(guān)鍵詞:模型

    曾慶橋 ,張亮 ,劉萍 ,王哲 ,王鏈 ,趙良言 ,李治平

    (1.中國石油華北油田分公司勘探開發(fā)研究院,河北 任丘 062552;2.成都北方石油勘探開發(fā)技術(shù)有限公司,四川 成都 610051;3.中國地質(zhì)大學(xué)(北京)能源學(xué)院,北京 100089;4.中國石油華北油田分公司第五采油廠,河北 辛集 052300)

    0 引言

    致密油藏作為重要的非常規(guī)油氣資源,常采用注水開發(fā)的方式進(jìn)行生產(chǎn),因此,研究水驅(qū)注水井動(dòng)態(tài)裂縫對(duì)于合理、高效的配注和科學(xué)開發(fā)致密油藏有著至關(guān)重要的作用[1-4]。注水井動(dòng)態(tài)裂縫的研究主要有流線模型方法、數(shù)值模擬方法、現(xiàn)代試井方法等[5-10],但是,對(duì)于定量分析注水井動(dòng)態(tài)裂縫,這些方法具有理論復(fù)雜、計(jì)算量大和所需測試數(shù)據(jù)多等不足[11-12]。而分析井間連通性對(duì)于研究注水井動(dòng)態(tài)裂縫具有簡潔快速的優(yōu)勢,常用的連通性模型主要包括Spearman相關(guān)分析模型[13]、人工神經(jīng)網(wǎng)絡(luò)(ANN)模型[14-15]、卡爾曼濾波法[16]、多 元 線 性 回 歸 模 型[17-18]、 容 阻 模 型 (CRM)[19-22]等 。Spearman模型計(jì)算過程簡單、易于求解,但相對(duì)理想化;ANN模型收斂速度較慢,操作過程比較復(fù)雜;卡爾曼濾波法和多元線性回歸模型都能夠反演井間關(guān)系和儲(chǔ)層特征,但不能考慮井底壓力變化對(duì)油田生產(chǎn)的影響;容阻模型既可以表征生產(chǎn)井對(duì)注水井信號(hào)響應(yīng)的延遲程度,又可以考慮變井底壓力生產(chǎn)的狀況[23-25]。因此,本研究采用容阻模型分析注水井動(dòng)態(tài)裂縫。

    本文考慮的注水井動(dòng)態(tài)裂縫特征主要有裂縫半長、裂縫高度和裂縫導(dǎo)流系數(shù)。首先,通過抽取不同裂縫參數(shù)樣本,結(jié)合數(shù)值模擬器獲得注采井生產(chǎn)數(shù)據(jù);然后,根據(jù)油藏容阻模型,獲得不同裂縫參數(shù)下的注采井井間連通系數(shù),并建立裂縫參數(shù)與連通系數(shù)的關(guān)系代理模型;最后,將實(shí)際生產(chǎn)數(shù)據(jù)分階段進(jìn)行容阻模型擬合,得到不同階段連通系數(shù),代入代理模型,獲得注水井裂縫動(dòng)態(tài)特征。

    1 油藏容阻模型

    1.1 容阻模型推導(dǎo)

    油藏容阻模型中[26-30],假設(shè)生產(chǎn)井有一個(gè)控制區(qū)域,每口注水井與每口生產(chǎn)井之間通過連通系數(shù)fij(i,j分別為注水井、生產(chǎn)井編號(hào))來表征注水井的注水分配情況,也反映出注采井之間的連通關(guān)系(見圖1)。

    圖1 油藏容阻模型示意

    模型中,假設(shè)油藏系統(tǒng)主要為單相微可壓縮流體,生產(chǎn)井卸油區(qū)內(nèi)流體綜合壓縮系數(shù)可描述為

    式中:ct為卸油區(qū)流體綜合壓縮系數(shù),1/MPa;Vp為卸油區(qū)容積,m3;pˉ為卸油區(qū)平均壓力,MPa;t為生產(chǎn)時(shí)間,d。

    對(duì)于生產(chǎn)井卸油區(qū)的封閉系統(tǒng),系統(tǒng)中隨時(shí)間而積累的流體體積等于流入流體與流出流體體積的代數(shù)和,將式(1)中平均壓力項(xiàng)采用線性產(chǎn)能模型消除,得到:

    式中:I為卸油區(qū)流入液量,m3/d;Q為卸油區(qū)流出液量(生產(chǎn)井產(chǎn)液量),m3/d;J 為生產(chǎn)井產(chǎn)液指數(shù),m3/(d·MPa);pwf為生產(chǎn)井井底流壓,MPa。

    1.2 容阻模型離散化

    為了將油藏開發(fā)各階段實(shí)際生產(chǎn)動(dòng)態(tài)數(shù)據(jù)代入油藏容阻模型進(jìn)行參數(shù)擬合,對(duì)式(4)進(jìn)行離散化處理。假設(shè)注入井Wi注入的流體中有一部分流向生產(chǎn)井Pj,Wi和Pj之間的連通系數(shù)為fij。忽略生產(chǎn)井井底流壓時(shí),生產(chǎn)井Pj在第n個(gè)月末時(shí)的產(chǎn)液量Qjn為

    式中:Qj1為生產(chǎn)井 Pj初期產(chǎn)液量,m3/d;tn,t1,tk分別為第n個(gè)月末、第1個(gè)月末、第k個(gè)時(shí)間步末對(duì)應(yīng)的時(shí)間,d;N為到第n個(gè)月末時(shí)總的時(shí)間步,d;Ni為總的注水井井?dāng)?shù),口;Iik為注水井Wi在第k時(shí)間步的注水量,m3/d。

    式中:Δpwf,jk為生產(chǎn)井Pj在第k個(gè)時(shí)間步的井底壓力壓降,MPa。

    1.3 注采井井間連通系數(shù)求解

    在已知生產(chǎn)井Pj的初始產(chǎn)量Qj1、產(chǎn)液指數(shù)Jj、時(shí)間常數(shù)τj以及各注采井井間連通系數(shù)fij時(shí),結(jié)合各注水井注入量,可根據(jù)式(6)快速計(jì)算各生產(chǎn)井的產(chǎn)液量。本研究中,結(jié)合注采井各階段實(shí)際產(chǎn)液量和井底流壓數(shù)據(jù),獲取注采井井間連通系數(shù)變化情況,以進(jìn)一步表征各階段井間裂縫變化情況。

    2 致密油藏?cái)?shù)值模擬方法

    為獲得水力壓裂后不同裂縫與井間連通系數(shù)的關(guān)系,利用不同裂縫配置下的致密油藏?cái)?shù)值模擬方法,得到不同裂縫下的注采井動(dòng)態(tài)生產(chǎn)數(shù)據(jù)。本研究數(shù)值模擬方法采用對(duì)裂縫表征較為精確的離散裂縫模型,獲得不同裂縫配置下的注采井生產(chǎn)數(shù)據(jù)。

    2.1 考慮水力裂縫的致密油藏?cái)?shù)值模擬方法

    考慮致密儲(chǔ)層非線性滲流特征、強(qiáng)應(yīng)力敏感特征及多組分多相流動(dòng)特征,描述基巖-裂縫系統(tǒng)流體傳導(dǎo)特征,參考文獻(xiàn)[31]以致密油藏滲流模型為基礎(chǔ),建立了致密油基巖-裂縫流動(dòng)耦合數(shù)學(xué)模型。水力裂縫致密油藏?cái)?shù)值模擬質(zhì)量守恒方程中,考慮巖石和流體的可壓縮性,基于控制體積的流動(dòng)方程為

    式中:Ω為積分控制單元體;φ為控制單元體儲(chǔ)層孔隙度;ρ為儲(chǔ)層流體的密度,kg/m3;v為儲(chǔ)層流體的滲流速度,m/s;q為控制單元體內(nèi)儲(chǔ)層流體的體積流量,m3/d;n為控制單元體外圍?Ω法向量;S,σ分別為曲面積分變量和體積分變量。

    考慮致密油藏基質(zhì)啟動(dòng)壓力梯度以及基質(zhì)和裂縫應(yīng)力敏感效應(yīng),得到基質(zhì)和裂縫中流體運(yùn)動(dòng)方程:

    式中:vm,vF分別為基質(zhì)網(wǎng)格和裂縫網(wǎng)格滲流速度,m/s;下標(biāo)m,F(xiàn)分別為基質(zhì)網(wǎng)格和裂縫網(wǎng)格;Km0,KF0分別為基質(zhì)網(wǎng)格和裂縫網(wǎng)格初始滲透率,10-3μm2;αm,αF分別為基質(zhì)和裂縫的應(yīng)力敏感系數(shù),MPa-1;p,p0分別為當(dāng)前壓力和原始?jí)毫Γ琈Pa;▽p為相鄰網(wǎng)格間壓力梯度,MPa/m;▽z為相鄰網(wǎng)格間重力方向高度梯度,m;ξ1,ξ2分別為啟動(dòng)壓力梯度項(xiàng)系數(shù);μ為流體黏度,mPa·s;g為重力加速度,m/s2。

    采用調(diào)和平均方法,計(jì)算基質(zhì)網(wǎng)格與基質(zhì)網(wǎng)格、基質(zhì)網(wǎng)格與裂縫網(wǎng)格、裂縫網(wǎng)格與裂縫網(wǎng)格之間的傳導(dǎo)率。計(jì)算公式分別為

    式中:T為網(wǎng)格間的傳導(dǎo)率;下標(biāo)a,b分別為不同的網(wǎng)格編號(hào)。

    2.2 致密油藏注水井裂縫動(dòng)態(tài)表征及模擬

    低滲致密油藏注水開發(fā)中,注水井動(dòng)態(tài)裂縫是油藏水驅(qū)及剩余油分布的主控因素之一。由于致密油藏物性差,“注不進(jìn)”導(dǎo)致的注入壓力變化往往引起注水井壓裂裂縫的動(dòng)態(tài)變化。注水井動(dòng)態(tài)裂縫受注水壓力和地應(yīng)力共同影響,本文考慮注水井動(dòng)態(tài)裂縫的變化特征主要體現(xiàn)在裂縫半長Lf、裂縫高度H和裂縫導(dǎo)流系數(shù)Fc3個(gè)方面, 用其參數(shù)組合 [Lf, H,F(xiàn)c]來表征某時(shí)刻的裂縫特征,其中,裂縫高度通過裂縫與儲(chǔ)層的縫高比表征,導(dǎo)流系數(shù)通過裂縫滲透率與裂縫寬度的乘積表征。

    圖2(圖中P為生產(chǎn)井,W為注水井)展示了注水井水力裂縫特征分別為[50,1,800]和[100,1,600]時(shí),模型剩余油分布情況。

    圖2 不同水力裂縫特征下的油藏含油飽和度分布

    3 代理模型方法及裂縫參數(shù)預(yù)測

    代理模型方法可以方便快速地建立不同數(shù)據(jù)之間的聯(lián)系,本文引入基于粒子群算法的最小二乘支持向量機(jī)方法(LS-SVM-PSO),建立起不同裂縫配置參數(shù)與注采井井間連通系數(shù)的關(guān)系。LS-SVM-PSO相比較于支持向量機(jī)(SVM)和最小二乘支持向量機(jī)(LSSVM),克服了SVM在大數(shù)據(jù)集中的應(yīng)用受到時(shí)間和內(nèi)存消耗優(yōu)化的限制?;诹W尤簝?yōu)化算法是一種模擬鳥群覓食行為的進(jìn)化算法,即隨機(jī)初始化粒子群,進(jìn)行LS-SVM訓(xùn)練,計(jì)算出每個(gè)粒子的適應(yīng)度值。若滿足終止條件,則輸出最優(yōu)參數(shù)值;否則,計(jì)算每個(gè)粒子和群體的歷史最優(yōu)值,并不斷更新粒子的位置和速度,反復(fù)迭代,直至滿足條件計(jì)算出最優(yōu)值。

    用變量 y=[Lf, H,F(xiàn)c]來表征某個(gè)特定狀態(tài)下的裂縫特征,采用網(wǎng)格搜索抽樣、拉丁立方體抽樣以及蒙特卡洛隨機(jī)抽樣方法獲得M組不同的裂縫參數(shù)組合。將M組裂縫特征變量分別代入數(shù)值模擬器,獲得M組不同裂縫參數(shù)下的注采井生產(chǎn)數(shù)據(jù)Q={Qm,m=1,2,…,M}。 若模擬區(qū)共有Ni口注水井,Np口生產(chǎn)井,則 q=[W1,W2,…,WNi;P1,P2,…,PNp]為Ni+Np維向量。進(jìn)一步將各注采井生產(chǎn)數(shù)據(jù)代入油藏容阻模型式(8)中,獲取不同裂縫參數(shù)配置下的各注采井井間連通系數(shù)。最后,采用LS-SVM-PSO建立不同裂縫參數(shù)對(duì)注采井井間連通系數(shù)的響應(yīng)模型:

    式中:αm為PSO算法優(yōu)化獲得的模型權(quán)重系數(shù);b為支持向量機(jī)模型位移量;K( x,xm)為核函數(shù)。

    4 實(shí)例應(yīng)用

    選取華北油田實(shí)際低滲區(qū)塊某致密油藏反五點(diǎn)井組(見圖 3,圖中 P1,P2,P3,P4為生產(chǎn)井),根據(jù)區(qū)塊的實(shí)際生產(chǎn)動(dòng)態(tài)參數(shù),為注水井設(shè)置不同裂縫參數(shù)進(jìn)行數(shù)值模擬,獲得不同裂縫參數(shù)下的模擬生產(chǎn)數(shù)據(jù)。

    圖3 某致密油藏反五點(diǎn)井組單元模型

    利用網(wǎng)格搜索抽樣原理,建立裂縫半長、裂縫高度和裂縫導(dǎo)流系數(shù)的設(shè)計(jì)樣本數(shù)據(jù):裂縫半長分別取0,20,40,60,80 m;裂縫高度通過縫高比表征,分別取0.2,0.4,0.6,0.8,1.0;裂縫導(dǎo)流系數(shù)分別取 100×10-3,200×10-3,300×10-3,400×10-3,500×10-3μm2·cm。 根據(jù)抽樣共設(shè)計(jì)5×5×5=125套數(shù)值模擬方案,以樣本裂縫半長 60 m、縫高比 0.8、裂縫導(dǎo)流系數(shù) 400×10-3μm2·cm為例,模擬生產(chǎn)2 700 d,數(shù)值模擬得到各生產(chǎn)井的產(chǎn)液量如圖4所示。

    圖4 特定裂縫參數(shù)下數(shù)值模擬各生產(chǎn)井產(chǎn)液量

    采用容阻模型擬合各生產(chǎn)井產(chǎn)液量(見圖5),用前70%生產(chǎn)時(shí)間對(duì)應(yīng)數(shù)據(jù)進(jìn)行模型訓(xùn)練,后30%數(shù)據(jù)進(jìn)行測試。根據(jù)容阻模型預(yù)測數(shù)據(jù)與模擬器真實(shí)產(chǎn)量數(shù)據(jù)擬合得到:P1井訓(xùn)練誤差為4.50%,測試誤差為1.62%;P2井訓(xùn)練誤差為8.09%,測試誤差為12.89%;P3井訓(xùn)練誤差為7.79%,測試誤差為2.82%;P4井訓(xùn)練誤差為5.35%,測試誤差為1.31%。各井平均訓(xùn)練誤差為6.43%,平均測試誤差為4.66%。擬合結(jié)果表明,油藏容阻模型能夠較好地反映油藏注采井間的動(dòng)態(tài)關(guān)系。在[60,0.8,400]條件下,注采井井間連通關(guān)系如圖6所示,其中,注水井至生產(chǎn)井的連通性大小由底邊中點(diǎn)在注水井、頂點(diǎn)位于生產(chǎn)井的等腰三角形表示,三角 形底邊長度與注采井之間的連通系數(shù)成正比。

    根據(jù)LS-SVM-PSO模型和數(shù)值模擬的計(jì)算結(jié)果,對(duì)比驗(yàn)證LS-SVM-PSO模型的正確性。從數(shù)值模擬過程產(chǎn)生的125組數(shù)據(jù)中,選擇100組進(jìn)行學(xué)習(xí),構(gòu)建LS-SVM代理模型,剩下25組對(duì)LS-SVM代理模型進(jìn)行精度驗(yàn)證,分別對(duì)動(dòng)態(tài)裂縫的裂縫長度、縫高比和裂縫導(dǎo)流系數(shù)進(jìn)行測試,其R2分別為0.948 9,0.925 8和0.934 8,驗(yàn)證了LS-SVM代理模型的可信度。

    獲得裂縫參數(shù)與井間連通系數(shù)的代理模型后,通過獲取不同時(shí)間段注采井井間連通系數(shù),進(jìn)一步結(jié)合代理模型,預(yù)測各時(shí)間段注水井壓裂裂縫特征。因?yàn)?口生產(chǎn)井位于模型的邊角區(qū)域,生產(chǎn)井泄壓范圍僅占1/4,根據(jù)流動(dòng)邊界數(shù)據(jù)對(duì)比,進(jìn)行生產(chǎn)井的產(chǎn)量劈分校正。

    針對(duì)各生產(chǎn)井校正后的產(chǎn)量數(shù)據(jù),將生產(chǎn)時(shí)間2 700 d等分為3個(gè)時(shí)間段,各階段分別采用油藏容阻模型擬合(見圖 7),P1,P2,P3,P4井的平均擬合誤差分別為3.46%,3.81%,3.76%,3.81%。各井各階段平均擬合誤差為3.72%,擬合效果較好。擬合出的各階段注采井井間連通關(guān)系如圖8所示。

    圖7 4口井分階段容阻模型擬合結(jié)果

    圖8 各階段注采井井間連通關(guān)系

    將各階段井間連通系數(shù)代入通過樣本數(shù)據(jù)建立的LS-SVM代理模型中,獲得3個(gè)階段注水井裂縫半長分別為26,64,41 m;注水井裂縫縫高比分別為0.34,0.66,0.48;注水井裂縫導(dǎo)流系數(shù)分別為286×10-3,510×10-3,352×10-3μm2·cm。

    將代理模型預(yù)測的裂縫參數(shù)進(jìn)行建模,設(shè)置3個(gè)時(shí)間段不同裂縫參數(shù)的井組數(shù)值模擬方案,數(shù)值模擬結(jié)果與原各井生產(chǎn)數(shù)據(jù)結(jié)果一致,結(jié)果表明,通過代理模型預(yù)測出的注水井壓裂裂縫動(dòng)態(tài)變化特征具有較好的可信性。

    5 結(jié)論

    1)本文建立了基于油藏容阻模型的井間連通系數(shù)反演方法,建立的容阻模型在井組生產(chǎn)數(shù)據(jù)擬合中精度較高,表明擬合參數(shù)能很好反映油藏注采井井間連通關(guān)系。

    2)基于建立的裂縫參數(shù)與連通系數(shù)的關(guān)系代理模型,提出了致密油藏注水井動(dòng)態(tài)裂縫定量分析方法,能快速實(shí)現(xiàn)基于生產(chǎn)動(dòng)態(tài)數(shù)據(jù)的不同階段注水井動(dòng)態(tài)裂縫參數(shù)預(yù)測。

    3)選取華北油田實(shí)際低滲區(qū)塊某致密油藏反五點(diǎn)井組,利用所建立的方法,快速實(shí)現(xiàn)了較為準(zhǔn)確的注水井動(dòng)態(tài)裂縫預(yù)測分析。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    伊人久久国产一区二区| 久久毛片免费看一区二区三区| a级毛片黄视频| 国产精品久久久久久av不卡| 久久午夜综合久久蜜桃| 看非洲黑人一级黄片| 天天操日日干夜夜撸| a级片在线免费高清观看视频| 久久久久久久国产电影| 国产黄色视频一区二区在线观看| 国产高清国产精品国产三级| 草草在线视频免费看| av在线观看视频网站免费| 在线看a的网站| 国产精品一区二区在线观看99| 日日摸夜夜添夜夜添av毛片| 99久久人妻综合| 国产日韩欧美视频二区| 美女内射精品一级片tv| 国国产精品蜜臀av免费| 极品人妻少妇av视频| 黄片无遮挡物在线观看| 欧美精品一区二区大全| .国产精品久久| 色94色欧美一区二区| 尾随美女入室| 午夜视频国产福利| 久久这里有精品视频免费| 成人午夜精彩视频在线观看| 高清视频免费观看一区二区| 在线天堂最新版资源| 欧美+日韩+精品| 日本欧美视频一区| 亚洲国产成人一精品久久久| 国产片特级美女逼逼视频| 在线观看美女被高潮喷水网站| 亚洲精品乱久久久久久| 婷婷色综合大香蕉| 亚洲av欧美aⅴ国产| 日韩,欧美,国产一区二区三区| 美女中出高潮动态图| 高清毛片免费看| 人妻系列 视频| 国产亚洲精品久久久com| 成人18禁高潮啪啪吃奶动态图 | 黄片播放在线免费| 国产亚洲欧美精品永久| 久久久久精品久久久久真实原创| 纵有疾风起免费观看全集完整版| 欧美激情极品国产一区二区三区 | 夜夜爽夜夜爽视频| 美女福利国产在线| 又大又黄又爽视频免费| 色5月婷婷丁香| 久久久久久久大尺度免费视频| 丝袜美足系列| 国产成人免费观看mmmm| 亚洲欧美一区二区三区黑人 | 欧美精品国产亚洲| 亚洲天堂av无毛| 观看av在线不卡| 日本午夜av视频| 亚洲精品自拍成人| 国产精品成人在线| 亚洲精品,欧美精品| 91aial.com中文字幕在线观看| 成人手机av| 制服诱惑二区| 国产熟女午夜一区二区三区 | 简卡轻食公司| 亚洲欧美中文字幕日韩二区| 下体分泌物呈黄色| 午夜激情福利司机影院| 欧美丝袜亚洲另类| 亚洲精品国产av蜜桃| 伊人久久国产一区二区| 国产一区有黄有色的免费视频| 久久久久久久亚洲中文字幕| 亚洲av不卡在线观看| 国产片特级美女逼逼视频| 国产综合精华液| 中文字幕最新亚洲高清| 免费观看性生交大片5| 91精品国产国语对白视频| 国产片内射在线| 成人漫画全彩无遮挡| 亚洲欧美日韩卡通动漫| 美女cb高潮喷水在线观看| 美女国产高潮福利片在线看| 欧美丝袜亚洲另类| 99久久人妻综合| 欧美97在线视频| 99久久精品一区二区三区| 国产精品人妻久久久久久| 少妇高潮的动态图| 午夜福利在线观看免费完整高清在| 成人毛片a级毛片在线播放| 精品久久久久久久久亚洲| 久久人人爽人人爽人人片va| 91在线精品国自产拍蜜月| 狂野欧美激情性xxxx在线观看| 欧美精品国产亚洲| 天美传媒精品一区二区| 99久国产av精品国产电影| 亚洲精品456在线播放app| tube8黄色片| 一区二区日韩欧美中文字幕 | 亚洲四区av| 丰满饥渴人妻一区二区三| 亚洲国产精品国产精品| 国产亚洲av片在线观看秒播厂| 免费黄网站久久成人精品| 亚洲av不卡在线观看| 国产精品蜜桃在线观看| 国产亚洲一区二区精品| 青春草国产在线视频| 韩国av在线不卡| 久久人人爽人人片av| 欧美xxxx性猛交bbbb| 最近2019中文字幕mv第一页| 51国产日韩欧美| 免费黄色在线免费观看| 在线免费观看不下载黄p国产| 美女国产高潮福利片在线看| 亚洲国产精品成人久久小说| 亚洲四区av| 春色校园在线视频观看| 国产在视频线精品| 精品久久久精品久久久| 婷婷色综合www| 最近手机中文字幕大全| 日韩免费高清中文字幕av| 大香蕉97超碰在线| 国产亚洲精品久久久com| av不卡在线播放| 国产在视频线精品| 国产精品嫩草影院av在线观看| av不卡在线播放| 下体分泌物呈黄色| 午夜激情福利司机影院| 国产精品国产三级国产专区5o| 国产精品国产三级国产专区5o| 热re99久久精品国产66热6| 蜜桃国产av成人99| 女的被弄到高潮叫床怎么办| 女的被弄到高潮叫床怎么办| 18禁动态无遮挡网站| 久久精品国产鲁丝片午夜精品| 乱码一卡2卡4卡精品| 国产日韩欧美在线精品| 一级毛片aaaaaa免费看小| 国产午夜精品一二区理论片| 日韩在线高清观看一区二区三区| 亚洲人成网站在线观看播放| 18+在线观看网站| 亚洲人成网站在线观看播放| 午夜福利,免费看| 国产欧美亚洲国产| 免费高清在线观看日韩| 亚洲精品国产av蜜桃| 久久久久久久久久久丰满| 国产乱人偷精品视频| 国产高清三级在线| 国产一区有黄有色的免费视频| 日韩电影二区| 国产精品不卡视频一区二区| 亚洲精品乱码久久久久久按摩| 日韩av免费高清视频| 亚洲欧美日韩卡通动漫| 如何舔出高潮| 成年女人在线观看亚洲视频| 日本-黄色视频高清免费观看| 亚洲欧美精品自产自拍| 国产女主播在线喷水免费视频网站| 午夜日本视频在线| 一区在线观看完整版| 亚洲国产av新网站| 久久久久久伊人网av| 母亲3免费完整高清在线观看 | 母亲3免费完整高清在线观看 | av国产精品久久久久影院| av电影中文网址| 99热这里只有精品一区| 国产色婷婷99| 欧美老熟妇乱子伦牲交| 久久久国产精品麻豆| 三级国产精品欧美在线观看| 日本与韩国留学比较| 观看av在线不卡| 各种免费的搞黄视频| 9色porny在线观看| 毛片一级片免费看久久久久| 性色avwww在线观看| 欧美xxxx性猛交bbbb| 三上悠亚av全集在线观看| 免费黄网站久久成人精品| 日韩电影二区| 少妇高潮的动态图| 三级国产精品欧美在线观看| 亚洲美女黄色视频免费看| 草草在线视频免费看| 一级片'在线观看视频| 日韩成人av中文字幕在线观看| 亚洲av免费高清在线观看| 99热全是精品| 18禁动态无遮挡网站| 精品久久久噜噜| 成年人免费黄色播放视频| 国产有黄有色有爽视频| av在线观看视频网站免费| 国产一区二区在线观看av| 免费久久久久久久精品成人欧美视频 | 国精品久久久久久国模美| 亚洲av免费高清在线观看| 国产无遮挡羞羞视频在线观看| av在线老鸭窝| 欧美激情 高清一区二区三区| 女人久久www免费人成看片| 人妻夜夜爽99麻豆av| 又大又黄又爽视频免费| 国产成人精品婷婷| 最近中文字幕高清免费大全6| 免费观看的影片在线观看| 天堂俺去俺来也www色官网| 中文字幕人妻熟人妻熟丝袜美| 欧美变态另类bdsm刘玥| 女人久久www免费人成看片| 日韩av不卡免费在线播放| 日本wwww免费看| 成人18禁高潮啪啪吃奶动态图 | 涩涩av久久男人的天堂| 亚洲精华国产精华液的使用体验| 人人妻人人澡人人爽人人夜夜| 日韩人妻高清精品专区| 菩萨蛮人人尽说江南好唐韦庄| 激情五月婷婷亚洲| 性色avwww在线观看| 日韩一区二区三区影片| 精品亚洲乱码少妇综合久久| 18+在线观看网站| 成人亚洲欧美一区二区av| freevideosex欧美| 最近最新中文字幕免费大全7| av国产精品久久久久影院| 欧美xxxx性猛交bbbb| 又黄又爽又刺激的免费视频.| 国产免费现黄频在线看| 亚洲美女黄色视频免费看| 高清av免费在线| 久久国内精品自在自线图片| 亚洲不卡免费看| 十八禁网站网址无遮挡| 亚洲精品日韩在线中文字幕| 免费人妻精品一区二区三区视频| 欧美bdsm另类| 日韩大片免费观看网站| 男女高潮啪啪啪动态图| 黄色毛片三级朝国网站| 一区二区三区乱码不卡18| 成年人午夜在线观看视频| 青春草国产在线视频| 熟妇人妻不卡中文字幕| 精品一区二区三卡| 国产毛片在线视频| 卡戴珊不雅视频在线播放| 日韩免费高清中文字幕av| 99热6这里只有精品| 亚洲美女搞黄在线观看| 成人综合一区亚洲| 久久久欧美国产精品| 91精品三级在线观看| 色网站视频免费| 欧美精品人与动牲交sv欧美| 91国产中文字幕| 十分钟在线观看高清视频www| 男女高潮啪啪啪动态图| 少妇猛男粗大的猛烈进出视频| 男人爽女人下面视频在线观看| 天堂俺去俺来也www色官网| 久久久久久人妻| 日韩电影二区| 国产在线免费精品| 午夜福利在线观看免费完整高清在| 欧美日韩成人在线一区二区| 极品人妻少妇av视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品中文字幕在线视频| av线在线观看网站| 日韩精品免费视频一区二区三区 | 亚洲不卡免费看| 人妻 亚洲 视频| 五月开心婷婷网| 久久热精品热| 少妇人妻久久综合中文| 日韩强制内射视频| 一级毛片我不卡| a级毛色黄片| 国产免费又黄又爽又色| 亚洲成人av在线免费| 日韩伦理黄色片| 久久综合国产亚洲精品| 国产精品蜜桃在线观看| 成人免费观看视频高清| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 免费观看的影片在线观看| 乱人伦中国视频| 国产精品.久久久| av福利片在线| xxxhd国产人妻xxx| 日本欧美国产在线视频| av在线老鸭窝| 国产有黄有色有爽视频| 成人国产av品久久久| 日本91视频免费播放| av.在线天堂| 免费久久久久久久精品成人欧美视频 | 成人亚洲欧美一区二区av| 草草在线视频免费看| 丰满少妇做爰视频| 精品亚洲成国产av| 亚洲av欧美aⅴ国产| 欧美少妇被猛烈插入视频| 免费黄色在线免费观看| 中国美白少妇内射xxxbb| 亚洲精品,欧美精品| 人人妻人人添人人爽欧美一区卜| 婷婷成人精品国产| 国产探花极品一区二区| 少妇丰满av| 最后的刺客免费高清国语| 老女人水多毛片| 日本黄大片高清| 有码 亚洲区| 一边摸一边做爽爽视频免费| 亚洲国产精品999| 成人黄色视频免费在线看| 能在线免费看毛片的网站| 精品亚洲成国产av| 人成视频在线观看免费观看| 夫妻午夜视频| 韩国高清视频一区二区三区| 久久久久久久久久久免费av| 午夜激情久久久久久久| 久久人人爽人人片av| 亚洲成色77777| 色哟哟·www| 99热这里只有精品一区| 日韩av不卡免费在线播放| 国产日韩一区二区三区精品不卡 | 欧美激情极品国产一区二区三区 | 人人妻人人澡人人看| 国产黄片视频在线免费观看| 国产精品蜜桃在线观看| 黄色配什么色好看| 又大又黄又爽视频免费| 久久国产亚洲av麻豆专区| 高清毛片免费看| 亚洲成人av在线免费| 最近最新中文字幕免费大全7| 夜夜骑夜夜射夜夜干| 国产精品麻豆人妻色哟哟久久| 你懂的网址亚洲精品在线观看| 成人国产av品久久久| 国产探花极品一区二区| 成人手机av| av在线播放精品| 久久久久久人妻| 色网站视频免费| 久久久久久伊人网av| 国产在视频线精品| 国产色婷婷99| 美女大奶头黄色视频| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| 亚洲欧美一区二区三区黑人 | 国产在线一区二区三区精| 熟女电影av网| 在线观看三级黄色| 国产精品99久久99久久久不卡 | 精品熟女少妇av免费看| 国产成人91sexporn| 日韩免费高清中文字幕av| 天天影视国产精品| 一级,二级,三级黄色视频| 97在线视频观看| 亚洲国产欧美日韩在线播放| 岛国毛片在线播放| av免费在线看不卡| 在现免费观看毛片| 国产午夜精品一二区理论片| 汤姆久久久久久久影院中文字幕| 天美传媒精品一区二区| 内地一区二区视频在线| 男女免费视频国产| 日韩三级伦理在线观看| 亚洲成人手机| 十八禁高潮呻吟视频| 欧美日韩综合久久久久久| www.色视频.com| 性色avwww在线观看| 丁香六月天网| 成年美女黄网站色视频大全免费 | av线在线观看网站| 久久久亚洲精品成人影院| 亚洲一级一片aⅴ在线观看| 午夜福利视频精品| 亚洲精品一区蜜桃| 免费大片黄手机在线观看| 亚洲天堂av无毛| 免费观看无遮挡的男女| 欧美丝袜亚洲另类| 五月开心婷婷网| 自拍欧美九色日韩亚洲蝌蚪91| 丝袜脚勾引网站| 中国美白少妇内射xxxbb| 久久精品国产自在天天线| 亚洲精品国产av成人精品| a级片在线免费高清观看视频| 人人妻人人澡人人看| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 少妇人妻 视频| 国产黄色视频一区二区在线观看| 黄色一级大片看看| 中文精品一卡2卡3卡4更新| 秋霞在线观看毛片| 久久久久久久久久人人人人人人| 国产成人精品婷婷| 欧美日韩亚洲高清精品| 成人毛片60女人毛片免费| 国产高清有码在线观看视频| 精品久久久噜噜| 免费高清在线观看视频在线观看| 精品卡一卡二卡四卡免费| 国产精品国产av在线观看| 亚洲美女视频黄频| 99热6这里只有精品| 大陆偷拍与自拍| 亚洲av电影在线观看一区二区三区| 国产深夜福利视频在线观看| 成人毛片a级毛片在线播放| 久久综合国产亚洲精品| 久久国产精品男人的天堂亚洲 | 成人毛片a级毛片在线播放| 一区二区三区四区激情视频| 国产淫语在线视频| 搡女人真爽免费视频火全软件| 国产黄色视频一区二区在线观看| 波野结衣二区三区在线| 久久午夜综合久久蜜桃| 国产一区亚洲一区在线观看| 人妻人人澡人人爽人人| 久久亚洲国产成人精品v| 亚洲国产欧美在线一区| 久久久久久久久久久丰满| 欧美日韩视频精品一区| 国产精品国产三级国产专区5o| 亚洲成人av在线免费| 国产黄片视频在线免费观看| 午夜免费男女啪啪视频观看| 国产男女内射视频| 国产色爽女视频免费观看| 哪个播放器可以免费观看大片| 国内精品宾馆在线| 国产精品99久久久久久久久| 成人国语在线视频| 国产69精品久久久久777片| 高清午夜精品一区二区三区| 黄片播放在线免费| 国产一区有黄有色的免费视频| 国产精品一国产av| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 国产男人的电影天堂91| 晚上一个人看的免费电影| 亚洲精品aⅴ在线观看| 91久久精品国产一区二区成人| 久热这里只有精品99| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级国产专区5o| 在线天堂最新版资源| 日韩成人av中文字幕在线观看| 成人综合一区亚洲| 视频区图区小说| 一边摸一边做爽爽视频免费| 你懂的网址亚洲精品在线观看| 综合色丁香网| 久久精品熟女亚洲av麻豆精品| 精品国产一区二区久久| 99热这里只有是精品在线观看| 欧美xxⅹ黑人| 在线观看免费高清a一片| 午夜免费男女啪啪视频观看| 七月丁香在线播放| 特大巨黑吊av在线直播| 五月开心婷婷网| 午夜激情久久久久久久| 熟妇人妻不卡中文字幕| 日产精品乱码卡一卡2卡三| 2021少妇久久久久久久久久久| 欧美激情极品国产一区二区三区 | 免费高清在线观看视频在线观看| 能在线免费看毛片的网站| 国产伦理片在线播放av一区| 亚洲综合色网址| 99久久人妻综合| 欧美日韩视频精品一区| 最近手机中文字幕大全| 岛国毛片在线播放| 青春草亚洲视频在线观看| 中文字幕久久专区| 亚洲伊人久久精品综合| 热re99久久精品国产66热6| av播播在线观看一区| 国产亚洲精品久久久com| 不卡视频在线观看欧美| 五月玫瑰六月丁香| 久久久久久久国产电影| 男女免费视频国产| 少妇猛男粗大的猛烈进出视频| 少妇的逼水好多| 少妇人妻久久综合中文| 妹子高潮喷水视频| 边亲边吃奶的免费视频| 大码成人一级视频| 亚洲欧美日韩另类电影网站| 国产精品.久久久| 人人妻人人澡人人看| 亚洲精品日本国产第一区| 女性生殖器流出的白浆| 国产精品一区二区在线不卡| 一个人免费看片子| 五月伊人婷婷丁香| 韩国av在线不卡| 精品人妻偷拍中文字幕| 性色avwww在线观看| 熟女电影av网| 99久久精品国产国产毛片| 亚洲第一区二区三区不卡| 中文字幕av电影在线播放| 亚洲,欧美,日韩| av天堂久久9| 18+在线观看网站| 一边亲一边摸免费视频| 午夜影院在线不卡| 七月丁香在线播放| 国产在线一区二区三区精| 精品久久久久久久久亚洲| 国产69精品久久久久777片| 18在线观看网站| 国产精品嫩草影院av在线观看| 熟女电影av网| 亚洲成人手机| 国产白丝娇喘喷水9色精品| 亚洲色图综合在线观看| 久久久久久久久久成人| 成人18禁高潮啪啪吃奶动态图 | 亚洲中文av在线| 免费观看性生交大片5| 国产伦理片在线播放av一区| 伦理电影免费视频| 一区二区日韩欧美中文字幕 | 国产精品久久久久久久电影| 国产免费一区二区三区四区乱码| 免费大片黄手机在线观看| 日韩成人av中文字幕在线观看| 亚洲,欧美,日韩| 9色porny在线观看| 另类亚洲欧美激情| 国国产精品蜜臀av免费| 青春草视频在线免费观看| 亚洲国产欧美在线一区| tube8黄色片| 亚洲一级一片aⅴ在线观看| 免费黄频网站在线观看国产| 亚洲国产精品专区欧美| 蜜桃久久精品国产亚洲av| 久久精品久久久久久噜噜老黄| 久久久久久久大尺度免费视频| 69精品国产乱码久久久| 午夜精品国产一区二区电影| 最黄视频免费看| 成年av动漫网址| 久久久久久伊人网av| 一个人免费看片子| 亚州av有码| 中文字幕久久专区| 日本猛色少妇xxxxx猛交久久| 久久久精品免费免费高清| 国产视频首页在线观看| 欧美老熟妇乱子伦牲交| 青青草视频在线视频观看| 精品久久久久久久久av| 一二三四中文在线观看免费高清| 99视频精品全部免费 在线| 免费日韩欧美在线观看| 91午夜精品亚洲一区二区三区| 另类精品久久| 亚洲av成人精品一二三区| 久久国产精品男人的天堂亚洲 | 97在线视频观看| 国产在线视频一区二区| 成人无遮挡网站| 亚洲av日韩在线播放| 少妇熟女欧美另类| 最新中文字幕久久久久| 亚洲精品一区蜜桃| 自拍欧美九色日韩亚洲蝌蚪91| 久久精品熟女亚洲av麻豆精品| a级毛片免费高清观看在线播放| kizo精华| 免费观看性生交大片5| 久久精品人人爽人人爽视色| av国产久精品久网站免费入址|