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

    二次壓力梯度三孔滲流模型及非線性滲流特征

    2020-08-13 03:00:52周琳瑯甘笑非淦文杰
    巖性油氣藏 2020年5期
    關(guān)鍵詞:壓力梯度溶洞滲流

    余 燕,周琳瑯,甘笑非,胡 燕,淦文杰,鄧 莊

    (中國(guó)石油西南油氣田公司川中油氣礦,四川遂寧 629000)

    0 引言

    碳酸鹽巖儲(chǔ)集空間類型多樣、孔隙結(jié)構(gòu)復(fù)雜,儲(chǔ)層裂縫、孔洞發(fā)育,具有很強(qiáng)的非均質(zhì)性[1]。鄭松青等[2]針對(duì)碳酸鹽巖縫洞型儲(chǔ)層呈離散介質(zhì)的特征,提出縫洞網(wǎng)絡(luò)模型,并在此基礎(chǔ)上建立了油水兩相流數(shù)學(xué)模型,簡(jiǎn)化了儲(chǔ)層描述參數(shù),計(jì)算成本低,可以進(jìn)行大規(guī)模的工程應(yīng)用;賈永祿等[3]基于基質(zhì)、微裂縫、溶蝕孔洞和大裂縫組成的四重介質(zhì),建立了縫洞型四重介質(zhì)的碳酸鹽巖試井模型;熊鈺等[4]針對(duì)大尺度的碳酸鹽巖縫洞型油氣藏,考慮縫洞的多尺度特征,建立了井打在大尺度裂縫上的裂縫線性流試井解釋模型;程飛[5]針對(duì)縫洞型碳酸鹽巖儲(chǔ)層的非均質(zhì)性和各向異性強(qiáng)的特征,建立了碳酸鹽巖油藏儲(chǔ)層類型識(shí)別模板,并劃分了3 種儲(chǔ)層類型;姜瑞忠等[6]針對(duì)低滲油藏存在的壓敏效應(yīng)和啟動(dòng)壓力梯度現(xiàn)象,建立了雙重介質(zhì)低滲油藏的斜井?dāng)?shù)學(xué)模型,并繪制了斜井試井曲線;杜鑫等[7]針對(duì)縫洞型油氣藏,建立了井筒與溶洞相連通情況下的流動(dòng)方程并給出了井筒和溶洞之間的壓力差表達(dá)式;劉建峰等[8]建立了孔洞和裂縫組成滲流空間的典型雙重介質(zhì)模型;王蓓等[9]針對(duì)磨溪龍王廟組縫洞型碳酸鹽巖儲(chǔ)層,建立了多尺度離散裂縫模型,明確了該氣藏高、低滲分布區(qū)域;姜瑞忠等[10]建立了縫洞型碳酸鹽巖油藏水平井分形非線性滲流模型,并結(jié)合測(cè)井?dāng)?shù)據(jù)將新模型應(yīng)用于礦場(chǎng)滲流參數(shù)的解釋,但這些有關(guān)碳酸鹽巖儲(chǔ)層的滲流模型均未考慮二次壓力梯度非線性項(xiàng)的影響。

    聶仁仕等[11-13]建立了考慮井儲(chǔ)和表皮效應(yīng)的不穩(wěn)定滲流數(shù)學(xué)模型,同時(shí),該模型考慮了二次壓力梯度非線性項(xiàng)的影響;Nie 等[14]針對(duì)底水油藏和厚油層油藏,建立了不同邊界條件下考慮二次壓力梯度非線性項(xiàng)的球面滲流方程;張林等[15]建立了不同邊界條件下考慮二次壓力梯度非線性項(xiàng)的三區(qū)復(fù)合滲流方程;張強(qiáng)等[16]建立了考慮啟動(dòng)壓力梯度和二次梯度非線性項(xiàng)影響的滲流模型;Wang 等[17]針對(duì)多層復(fù)合儲(chǔ)層,建立了考慮二次壓力梯度非線性項(xiàng)的多區(qū)復(fù)合非線性模型;Guo 等[18]針對(duì)地下地層中氣井生產(chǎn),建立了考慮二次壓力梯度非線性項(xiàng)的直井模型;王美楠等[19]建立了考慮啟動(dòng)壓力梯度、動(dòng)邊界、二次梯度非線性項(xiàng)的低滲透變形介質(zhì)油藏滲流模型;Lu 等[20]建立了考慮二次壓力梯度非線性項(xiàng)的多孔介質(zhì)水平井滲流模型。目前,針對(duì)碳酸鹽巖縫洞型三孔介質(zhì)儲(chǔ)層,尚無有關(guān)二次壓力梯度三孔滲流模型及非線性滲流特征的研究。

    流體在地下儲(chǔ)層中的滲流具有二次壓力梯度的非線性滲流特性[18-20],縫洞型碳酸鹽巖儲(chǔ)層也不例外,采用常規(guī)的三孔介質(zhì)線性滲流模型無法精準(zhǔn)描述流體在縫洞型碳酸鹽巖儲(chǔ)層中的滲流特征。筆者首次建立二次壓力梯度三孔滲流模型,研究縫洞型碳酸鹽巖儲(chǔ)層受二次壓力梯度控制的非線性滲流模型的缺失問題,以期實(shí)現(xiàn)對(duì)縫洞型碳酸鹽巖儲(chǔ)層滲流特征的精準(zhǔn)描述。

    1 滲流模型

    1.1 物理模型

    平面徑向流物理模型假設(shè)如下:①裂縫孔洞型三孔介質(zhì)水平等厚地層中心一口井以定產(chǎn)量q生產(chǎn);②地層厚度為h,地層滲透率為常數(shù)且與壓力無關(guān),原始地層壓力為pi;③井半徑為rw,井底壓力為pwf;④井到地層外邊界的距離為re,外邊界可以是無窮大地層、封閉地層或恒壓邊界;⑤在彈性驅(qū)動(dòng)方式下,黏度為μ的單相液體滲流符合等溫達(dá)西定律;⑥考慮天然裂縫與井筒溝通,溶洞系統(tǒng)和基質(zhì)系統(tǒng)中流體均向天然裂縫系統(tǒng)發(fā)生擬穩(wěn)態(tài)竄流,因溶洞滲透率大于基質(zhì)滲透率,故溶洞系統(tǒng)中流體優(yōu)先向天然裂縫系統(tǒng)竄流;⑦考慮井筒儲(chǔ)集效應(yīng)和表皮效應(yīng)的影響。

    1.2 數(shù)學(xué)模型及其求解

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

    一個(gè)完整的數(shù)學(xué)模型由滲流控制微分方程、初始條件、內(nèi)外邊界條件構(gòu)成。對(duì)縫洞型碳酸鹽巖儲(chǔ)層,其滲流控制微分方程包括裂縫系統(tǒng)方程、基質(zhì)系統(tǒng)方程和溶洞系統(tǒng)方程。

    (1)滲流控制微分方程

    裂縫系統(tǒng)方程:

    式中:pfD為無因次裂縫系統(tǒng)壓力;pmD為無因次基質(zhì)系統(tǒng)壓力;pvD為無因次溶洞系統(tǒng)壓力;rD為無因次徑向距離;β為無因次二次壓力梯度項(xiàng)系數(shù);λmf為基質(zhì)系統(tǒng)向裂縫系統(tǒng)的竄流系數(shù);λvf為溶洞系統(tǒng)向裂縫系統(tǒng)的竄流系數(shù);ωf為裂縫系統(tǒng)的彈性儲(chǔ)容比;S為井的表皮系數(shù);tD為無因次時(shí)間。下標(biāo)f,m,v,D 分別代表裂縫(fracture)、基質(zhì)(matrix)、溶洞(vug)和無因次(Dimensionless)。

    基質(zhì)系統(tǒng)方程:

    式中:ωm為基質(zhì)系統(tǒng)的彈性儲(chǔ)容比。

    溶洞系統(tǒng)方程:

    式中:ωv為溶洞系統(tǒng)的彈性儲(chǔ)容比。

    (2)初始條件

    (3)內(nèi)邊界條件

    式中:CD為無因次井筒儲(chǔ)集系數(shù)。

    (4)外邊界條件

    無限大外邊界

    恒壓邊界

    式中:reD為無因次外邊界距離。

    封閉邊界

    上述數(shù)學(xué)模型的無因次參數(shù)定義如下:無因次壓力

    式中:pf為裂縫系統(tǒng)壓力,MPa;pm為基質(zhì)系統(tǒng)壓力,MPa;pv為溶洞系統(tǒng)壓力,MPa;pi為初始油藏壓力,MPa;q為流量,m3/d;μ為黏度,mPa·s;B為地層流體的體積系數(shù),無因次;kf為裂縫滲透率,D;h為儲(chǔ)層厚度,m。

    基于有效井徑定義的無因次距離

    式中:r為徑向距離,m;rw為井半徑,m。

    無因次時(shí)間

    式中:t為時(shí)間,h;φf為裂縫孔隙度;φv為溶洞孔隙度;φm為基質(zhì)孔隙度;Cft為裂縫彈性綜合壓縮系數(shù),MPa-1;Cvt為溶洞彈性綜合壓縮系數(shù),MPa-1;Cmt為基質(zhì)彈性綜合壓縮系數(shù),MPa-1。

    無因次井筒儲(chǔ)集系數(shù)

    式中:Cs為井的彈性儲(chǔ)容系數(shù),m3/MPa。

    彈性儲(chǔ)容比

    無因次二次壓力梯度項(xiàng)系數(shù)

    式中:Cρ為地層流體的壓縮系數(shù),MPa-1。

    1.2.2 積分變換求解模型

    引入基于tD的拉氏變換

    式中:L[]為L(zhǎng)aplace 算子為變量pD對(duì)應(yīng)的La‐place 空間變量;u為L(zhǎng)aplace 變量。

    對(duì)式(1)—(3)經(jīng)基于tD的Laplace 變換后,得

    當(dāng)ωf=1 時(shí),非線性模型則簡(jiǎn)化為常規(guī)單一介模型,式(22)可以簡(jiǎn)化為質(zhì)

    換句話講,在相同的油井生產(chǎn)和地層邊界條件下,Laplace 空間下的三孔介質(zhì)模型和均質(zhì)模型的差異僅僅表現(xiàn)在f(u)的函數(shù)表達(dá)式不同。因此,如果首先在Laplace 空間中得到均質(zhì)模型的解,就可以直接通過改變f(u)來得到Laplace 空間下三孔介質(zhì)模型的解。

    先令ωf=1,則考慮單一介質(zhì)情形,滲流控制方程中就不存在竄流項(xiàng),則式(1)變?yōu)?/p>

    作如下變量代換,將式(24)線性化

    式中:x為無因次裂縫壓力的中間代換變量。

    再經(jīng)基于tD的Laplace 變換后,可得到Laplace空間中的數(shù)學(xué)模型

    1.2.3 數(shù)學(xué)模型的解

    用式(22)替換式(27),就可得到經(jīng)變量代換線性化后的裂縫孔洞型三孔介質(zhì)滲流偏微分控制方程

    (1)無窮大外邊界模型的解

    故無窮大外邊界條件模型在Laplace 空間的通解為

    井底rD=1,當(dāng)p=pwf時(shí),pD=pwD,則x=xw,故無窮大地層模型井底動(dòng)態(tài)壓力在Laplace 空間為

    (2)恒壓外邊界模型的解

    將式(34)代入恒壓外邊界條件式(30),有

    聯(lián)解式(35),(36),(39),將其寫為矩陣形式,利用克萊姆法則或高斯迭代消元法,即可解出系數(shù)A和B,從而得出外邊界條件為恒壓邊界時(shí)井底壓力動(dòng)態(tài)在Laplace 空間的解。

    (3)封閉外邊界模型的解

    將式(34)代入封閉外邊界條件式(31),有

    聯(lián)解式(35),(36),(40),將其寫為矩陣形式,利用克萊姆法則或高斯迭代消元法,即可解出系數(shù)A和B,從而得出外邊界條件為封閉邊界時(shí)井底壓力動(dòng)態(tài)在Laplace 空間的解。

    對(duì)上述三類外邊界條件下Laplace 空間的解,利用Stehfest 數(shù)值算法進(jìn)行數(shù)值反演,得到實(shí)空間的解xw,再根據(jù)壓力pD與x的變換關(guān)系式,可求出裂縫孔洞型三孔介質(zhì)非線性滲流模型實(shí)空間井底壓力動(dòng)態(tài)響應(yīng)的數(shù)值解,進(jìn)而編程繪制pwD~tD/CD與p'wDtD/CD~tD/CD的無因次雙對(duì)數(shù)滲流特征曲線。

    2 非線性滲流特征分析

    2.1 非線性滲流過程分析

    圖1 為受非線性項(xiàng)系數(shù)影響的無窮大裂縫孔洞型三孔介質(zhì)地層擬穩(wěn)態(tài)竄流非線性滲流模型的典型曲線。該曲線主要受無因次二次梯度非線性項(xiàng)系數(shù)β控制。在其他參數(shù)一定時(shí)(CD=10,S=10,ωf=0.001,ωm=0.989,ωv=0.01,λmf=1×10-10,λvf=1×10-8),通過分別設(shè)定β值為0,0.01 和0.05,來分析非線性滲流特征曲線與線性滲流特征曲線的差異。當(dāng)β=0 時(shí),該曲線則蛻化為線性模型的滲流典型曲線(圖1 曲線①)。從圖1 可以看出,非線性滲流特征曲線向下偏離線性滲流特征曲線,隨β增大,偏離量越大。圖1 主要有如下7 個(gè)流動(dòng)階段:

    第Ⅰ階段:純井筒儲(chǔ)集階段。此流動(dòng)階段的曲線特征不受非線性二次壓力梯度項(xiàng)的影響。

    第Ⅱ階段:井筒儲(chǔ)集和表皮效應(yīng)階段。在此階段,非線性二次壓力梯度項(xiàng)的影響開始出現(xiàn),非線性滲流特征曲線也開始偏離線性滲流特征曲線。

    第Ⅲ階段:早期裂縫徑向流階段。裂縫系統(tǒng)中的流體徑向流入井筒,而溶洞系統(tǒng)和基質(zhì)系統(tǒng)中的流體還未開始流動(dòng)。與第Ⅱ階段相比,線性模型和非線性模型典型曲線的差異更為明顯。

    圖1 受非線性項(xiàng)系數(shù)(β)影響的三孔介質(zhì)擬穩(wěn)態(tài)竄流模型典型曲線Fig.1 Type curves of pressure transients controlled by varying nonlinear coefficient(β)

    第Ⅳ階段:溶洞系統(tǒng)向裂縫系統(tǒng)竄流階段。壓力導(dǎo)數(shù)曲線呈“V 形”,反映了溶洞系統(tǒng)向裂縫系統(tǒng)的竄流。線性模型和非線性模型典型曲線的差異進(jìn)一步增加。

    第Ⅴ階段:裂縫系統(tǒng)和溶洞系統(tǒng)整體徑向流階段。非線性二次壓力梯度項(xiàng)的影響越來越明顯,其壓力導(dǎo)數(shù)曲線不符合“0.5 線”定律,而是位于“0.5線”下方。

    第Ⅵ階段:基質(zhì)系統(tǒng)向裂縫系統(tǒng)竄流階段。壓力導(dǎo)數(shù)曲線呈“V 形”,表明基質(zhì)系統(tǒng)開始向裂縫系統(tǒng)竄流。線性模型和非線性模型特征曲線的差異越來越大。

    第Ⅶ階段:裂縫系統(tǒng)、溶洞系統(tǒng)和基質(zhì)系統(tǒng)整體徑向流階段。非線性二次壓力梯度項(xiàng)的影響越來越明顯,其壓力導(dǎo)數(shù)曲線不符合“0.5 線”定律,而是位于“0.5 線”下方,隨著時(shí)間的推移,壓力導(dǎo)數(shù)曲線逐漸偏離“0.5 線”,如圖1 中曲線②和曲線③所示。隨著時(shí)間的推移,非線性模型的壓力導(dǎo)數(shù)曲線逐漸偏離線性模型的壓力導(dǎo)數(shù)曲線。

    非線性三孔介質(zhì)曲線的外邊界反映特征與常規(guī)線性模型曲線的外邊界反映特征類似,此處不再贅述。

    該曲線除了受β控制,同時(shí)還受別的因素影響。圖2 為受溶洞向裂縫竄流的竄流系數(shù)(λvf)的非線性滲流模型典型曲線。在其他參數(shù)一定時(shí)(CD=10,S=1,ωf=0.001,ωm=0.989,ωv=0.01,β=0.05),分別設(shè)定λvf值為1×10-7,1×10-8和1×10-9。隨λvf增加,溶洞向裂縫竄流時(shí)間推遲,在圖2 中表現(xiàn)為壓力導(dǎo)數(shù)曲線的“V 形”向左移動(dòng)。

    圖2 受溶洞向裂縫竄流的竄流系數(shù)(λvf)影響的典型曲線Fig.2 Type curves of pressure transients controlled by varying inter-porosity flow factor(λvf)

    圖3 為受基質(zhì)向裂縫竄流的竄流系數(shù)(λmf)的非線性滲流模型典型曲線。在其他參數(shù)一定時(shí)(CD=10,S=1,ωf=0.001,ωm=0.989,ωv=0.01,β=0.05),分別設(shè)定λmf值為1×10-9,1×10-10和1×10-11。溶洞向裂縫竄流時(shí)間推遲,在圖3 中表現(xiàn)為壓力導(dǎo)數(shù)曲線的“V 形”向右移動(dòng)。

    圖3 受基質(zhì)向裂縫竄流的竄流系數(shù)(λmf)影響的典型曲線Fig.3 Type curves of pressure transients controlled by varying inter-porosity flow factor(λmf)

    圖4 為受裂縫系統(tǒng)彈性儲(chǔ)容比(ωf)、基質(zhì)系統(tǒng)彈性儲(chǔ)容比(ωm)影響的非線性滲流模型典型曲線。其他參數(shù)一定時(shí)(CD=10,S=1,ωv=0.1,λvf=1×10-10,λmf=1×10-8,β=0.05),分別設(shè)定ωf值為0.001,0.01 和0.08,相應(yīng)的ωm值為0.889,0.89 和0.82。隨ωf的增加,在圖4 中表現(xiàn)為壓力導(dǎo)數(shù)曲線的“V 形”變淺、變窄。

    圖5 為受溶洞系統(tǒng)彈性儲(chǔ)容比(ωv)、基質(zhì)系統(tǒng)彈性儲(chǔ)容比(ωm)影響的非線性滲流模型典型曲線。其他參數(shù)一定時(shí)(CD=10,S=1,ωf=0.001,λvf=1×10-10,λmf=1×10-8,β=0.05),分別設(shè)定ωv值為0.01,0.01 和0.5,相應(yīng)的ωm值為0.989,0.889 和0.499。隨ωv的增加,在圖5 中表現(xiàn)為壓力導(dǎo)數(shù)曲線的第1個(gè)“V 形”變得更深、更寬,第2 個(gè)“V 形”變得更窄、更淺。

    圖4 受裂縫系統(tǒng)彈性儲(chǔ)容比(ωf)影響的典型曲線Fig.4 Type curves of pressure transients controlled by fluid capacitance coefficient(ωf)

    圖5 受溶洞系統(tǒng)彈性儲(chǔ)容比(ωv)影響的典型曲線Fig.5 Type curves of pressure transients controlled by fluid capacitance coefficient(ωv)

    2.2 非線性滲流模型典型曲線的偏差分析

    為了定量描述非線性項(xiàng)系數(shù)的影響,定義絕對(duì)偏差(DV)與相對(duì)偏差(RDV):

    其他參數(shù)一定時(shí)(CD=10,S=1,ωf=0.001,ωm=0.989,ωv=0.01,λmf=1×10-10,λvf=1×10-8),分別設(shè)定β值為0.01 和0.05,分析非線性項(xiàng)系數(shù)對(duì)滲流典型曲線的影響(表1—2)。

    當(dāng)β=0.01 時(shí),由表1 中數(shù)據(jù)知:當(dāng)tD/CD=103時(shí),無因次壓力曲線的絕對(duì)偏差和相對(duì)偏差分別為0.579 8 和8.83%,無因次壓力導(dǎo)數(shù)曲線的絕對(duì)偏差和相對(duì)偏差分別為0.085 4 和17.37%;當(dāng)tD/CD=107時(shí),無因次壓力曲線的絕對(duì)偏差和相對(duì)偏差分別為1.4 和13.43%,無因次壓力導(dǎo)數(shù)曲線的絕對(duì)偏差和相對(duì)偏差分別為0.089 7 和25.40%。

    表1 非線性與線性滲流典型曲線的定量偏差數(shù)據(jù)(β=0.01)Table 1 Theoretical offset of type curves between the linear and nonlinear models(β=0.01)

    表2 非線性與線性滲流典型曲線的定量偏差數(shù)據(jù)(β=0.05)Table 2 Theoretical offset of type curves between the linear and nonlinear models(β=0.05)

    當(dāng)β=0.05,由表2 中數(shù)據(jù)知:當(dāng)tD/CD=103時(shí),無因次壓力曲線的絕對(duì)偏差和相對(duì)偏差分別為2.135 0 和32.52%,無因次壓力導(dǎo)數(shù)曲線的絕對(duì)偏差和相對(duì)偏差分別為0.272 4 和55.40%;當(dāng)tD/CD=107時(shí),無因次壓力曲線的絕對(duì)偏差和相對(duì)偏差分別為4.527 2 和43.43%,無因次壓力導(dǎo)數(shù)曲線的絕對(duì)偏差和相對(duì)偏差分別為0.241 1 和68.28%。

    絕對(duì)偏差和相對(duì)偏差均隨著無因次時(shí)間的增加而增加。在同一時(shí)刻,壓力導(dǎo)數(shù)的相對(duì)偏差值大于壓力的相對(duì)偏差值,例如,當(dāng)tD/CD=103,β=0.01時(shí),壓力導(dǎo)數(shù)的相對(duì)偏差值就(17.37%)大于壓力的相對(duì)偏差(8.83%)。絕對(duì)偏差和相對(duì)偏差值均隨β的增加而增加,例如,當(dāng)tD/CD=103時(shí),β=0.01的壓力相對(duì)偏差值為8.83%,而β=0.05 的壓力相對(duì)偏差值為32.52%。

    3 結(jié)論

    (1)針對(duì)裂縫孔洞型三孔介質(zhì)儲(chǔ)層建立的二次壓力梯度非線性滲流模型比不考慮二次壓力梯度的常規(guī)線性滲流模型能更精準(zhǔn)地描述儲(chǔ)層的滲流特征。

    (2)隨油井生產(chǎn)時(shí)間的增加,裂縫孔洞型三孔介質(zhì)儲(chǔ)層非線性滲流特征曲線逐漸偏離常規(guī)線性滲流特征曲線,且低于常規(guī)線性滲流特征曲線。

    (3)非線性與線性滲流特征曲線之間的偏離量可用“絕對(duì)偏差”和“相對(duì)偏差”等2 個(gè)參數(shù)來定量描述,其偏差隨著二次壓力梯度非線性系數(shù)的增加而增加。當(dāng)非線性系數(shù)為0 時(shí),非線性三孔滲流模型可蛻化為常規(guī)線性三孔滲流模型。

    猜你喜歡
    壓力梯度溶洞滲流
    別有洞天
    出發(fā)吧,去溶洞
    神秘的溶洞
    幼兒100(2017年31期)2017-11-27 02:37:45
    壓力梯度在油田開發(fā)中的應(yīng)用探討
    疊加原理不能求解含啟動(dòng)壓力梯度滲流方程
    致密砂巖啟動(dòng)壓力梯度數(shù)值的影響因素
    斷塊油氣田(2014年5期)2014-03-11 15:33:45
    簡(jiǎn)述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計(jì)算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    神秘的溶洞
    尾礦壩滲流計(jì)算及排滲設(shè)計(jì)
    金屬礦山(2013年6期)2013-03-11 16:54:05
    一级毛片黄色毛片免费观看视频| 中国三级夫妇交换| 国产av国产精品国产| 九色亚洲精品在线播放| 日韩一区二区视频免费看| 女性生殖器流出的白浆| 少妇人妻 视频| 欧美xxxx性猛交bbbb| 午夜精品国产一区二区电影| 日韩,欧美,国产一区二区三区| 亚洲精品视频女| 亚洲综合色网址| 人人澡人人妻人| 亚洲人与动物交配视频| 国产不卡av网站在线观看| 你懂的网址亚洲精品在线观看| 26uuu在线亚洲综合色| av一本久久久久| www日本在线高清视频| 日日啪夜夜爽| 一级毛片黄色毛片免费观看视频| 日韩大片免费观看网站| 2021少妇久久久久久久久久久| 老熟女久久久| 久久久国产精品麻豆| 制服诱惑二区| 欧美日韩视频高清一区二区三区二| 下体分泌物呈黄色| 亚洲欧美一区二区三区国产| 亚洲精品久久成人aⅴ小说| 日韩一区二区视频免费看| 国产精品一区二区在线观看99| 日本色播在线视频| 久久97久久精品| 亚洲av综合色区一区| 国产不卡av网站在线观看| 亚洲精品视频女| 亚洲av福利一区| 国产男人的电影天堂91| 精品亚洲成a人片在线观看| 国产成人欧美| 色婷婷久久久亚洲欧美| 国产成人免费无遮挡视频| 亚洲 欧美一区二区三区| xxxhd国产人妻xxx| 午夜福利影视在线免费观看| 亚洲欧美日韩另类电影网站| 99热6这里只有精品| 永久免费av网站大全| 亚洲色图 男人天堂 中文字幕 | 亚洲成人av在线免费| 中文字幕亚洲精品专区| 黄网站色视频无遮挡免费观看| 亚洲国产精品成人久久小说| 精品久久久精品久久久| 日本欧美国产在线视频| 一级黄片播放器| 免费日韩欧美在线观看| 亚洲欧美清纯卡通| 最近中文字幕高清免费大全6| 国产日韩一区二区三区精品不卡| 久久影院123| 久久久久久久精品精品| 成人漫画全彩无遮挡| 午夜免费鲁丝| 在线观看美女被高潮喷水网站| 黄片播放在线免费| 久久精品国产综合久久久 | 亚洲精华国产精华液的使用体验| 亚洲欧美日韩另类电影网站| 国产亚洲精品久久久com| 内地一区二区视频在线| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 99热这里只有是精品在线观看| 日韩在线高清观看一区二区三区| 久久99精品国语久久久| 亚洲成人手机| 亚洲图色成人| 亚洲经典国产精华液单| xxx大片免费视频| 国产亚洲精品久久久com| 国产xxxxx性猛交| 高清黄色对白视频在线免费看| 在线观看三级黄色| 国产精品一区www在线观看| 精品视频人人做人人爽| 在线观看免费日韩欧美大片| 在线观看免费高清a一片| 18禁观看日本| 久久女婷五月综合色啪小说| 国产一区亚洲一区在线观看| 在线免费观看不下载黄p国产| av视频免费观看在线观看| 一级片免费观看大全| 免费av中文字幕在线| 97精品久久久久久久久久精品| 人妻 亚洲 视频| 日韩av不卡免费在线播放| 激情五月婷婷亚洲| 99热全是精品| 伦理电影大哥的女人| 少妇人妻 视频| 亚洲欧美清纯卡通| 亚洲精品久久成人aⅴ小说| 午夜日本视频在线| 日日爽夜夜爽网站| 男女无遮挡免费网站观看| 一级,二级,三级黄色视频| 久久ye,这里只有精品| 韩国高清视频一区二区三区| 欧美最新免费一区二区三区| 精品一区二区免费观看| 热99国产精品久久久久久7| 久久狼人影院| 国产精品一国产av| 国语对白做爰xxxⅹ性视频网站| 国产高清三级在线| 欧美精品一区二区免费开放| 午夜91福利影院| 丝袜喷水一区| 2022亚洲国产成人精品| 亚洲婷婷狠狠爱综合网| 午夜免费观看性视频| av在线观看视频网站免费| 综合色丁香网| 我要看黄色一级片免费的| 亚洲三级黄色毛片| 亚洲一区二区三区欧美精品| 少妇猛男粗大的猛烈进出视频| 日日啪夜夜爽| 日韩一本色道免费dvd| 校园人妻丝袜中文字幕| 久久久久久久大尺度免费视频| 成人亚洲精品一区在线观看| 国产激情久久老熟女| 久久久精品区二区三区| 最近中文字幕2019免费版| 精品午夜福利在线看| a 毛片基地| a级片在线免费高清观看视频| 中文欧美无线码| 日日爽夜夜爽网站| 日韩一区二区视频免费看| 亚洲美女黄色视频免费看| 啦啦啦视频在线资源免费观看| 在线观看人妻少妇| 国产精品一国产av| 亚洲精品日本国产第一区| 宅男免费午夜| 日韩精品有码人妻一区| 2021少妇久久久久久久久久久| 久久毛片免费看一区二区三区| 一区在线观看完整版| 美女内射精品一级片tv| 久热这里只有精品99| 一区二区三区精品91| 免费观看性生交大片5| 精品亚洲成a人片在线观看| 久久久久久久久久久久大奶| 啦啦啦视频在线资源免费观看| 国产一区有黄有色的免费视频| 久久久国产精品麻豆| 国产免费一级a男人的天堂| www日本在线高清视频| 2021少妇久久久久久久久久久| 国产 精品1| 国产精品蜜桃在线观看| 日韩电影二区| 久久国产精品大桥未久av| 在线观看免费日韩欧美大片| 国产欧美另类精品又又久久亚洲欧美| 最近中文字幕高清免费大全6| 久久精品夜色国产| 高清毛片免费看| 久久久欧美国产精品| 日本免费在线观看一区| 亚洲综合色惰| 亚洲综合精品二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 丝袜脚勾引网站| 午夜福利,免费看| 久久久久久伊人网av| 欧美xxxx性猛交bbbb| 少妇高潮的动态图| 国产亚洲午夜精品一区二区久久| 新久久久久国产一级毛片| 日韩av不卡免费在线播放| 啦啦啦视频在线资源免费观看| 日本欧美视频一区| 最新中文字幕久久久久| 亚洲精品成人av观看孕妇| 国产一区二区激情短视频 | 水蜜桃什么品种好| 国产有黄有色有爽视频| 人人妻人人添人人爽欧美一区卜| 日韩视频在线欧美| 一区在线观看完整版| 国内精品宾馆在线| 岛国毛片在线播放| 伦理电影免费视频| 97精品久久久久久久久久精品| 国产又爽黄色视频| 亚洲国产日韩一区二区| 国产成人精品无人区| 制服人妻中文乱码| 国产日韩欧美在线精品| 日本wwww免费看| 人妻少妇偷人精品九色| 亚洲av电影在线进入| 男女边吃奶边做爰视频| 乱码一卡2卡4卡精品| 少妇高潮的动态图| 国产一级毛片在线| 国产成人欧美| h视频一区二区三区| 亚洲久久久国产精品| 国产成人精品福利久久| 国产永久视频网站| 少妇的逼好多水| 在现免费观看毛片| 男人添女人高潮全过程视频| 纵有疾风起免费观看全集完整版| 亚洲内射少妇av| 成年人午夜在线观看视频| 欧美精品高潮呻吟av久久| 免费观看av网站的网址| 中国三级夫妇交换| 亚洲av免费高清在线观看| 久久久久久久亚洲中文字幕| 最近中文字幕2019免费版| 日本vs欧美在线观看视频| 欧美xxxx性猛交bbbb| 国产高清不卡午夜福利| 校园人妻丝袜中文字幕| 蜜桃在线观看..| 色婷婷av一区二区三区视频| 久热久热在线精品观看| 一区二区三区精品91| 亚洲四区av| 欧美激情极品国产一区二区三区 | 一级爰片在线观看| 精品久久久久久电影网| 国产无遮挡羞羞视频在线观看| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久av不卡| 国产成人午夜福利电影在线观看| 午夜激情久久久久久久| 免费黄频网站在线观看国产| 久久精品久久久久久噜噜老黄| 国产免费一级a男人的天堂| 欧美人与性动交α欧美精品济南到 | 国产乱来视频区| 99国产综合亚洲精品| 色视频在线一区二区三区| 日韩电影二区| 国产一区亚洲一区在线观看| 91成人精品电影| 国产不卡av网站在线观看| 成人国产麻豆网| 日韩一区二区三区影片| 中文天堂在线官网| 五月天丁香电影| 精品少妇黑人巨大在线播放| 欧美精品av麻豆av| 国产淫语在线视频| 日本vs欧美在线观看视频| 国产精品久久久久久精品古装| 在线看a的网站| 久久国产亚洲av麻豆专区| 亚洲人成网站在线观看播放| 精品久久国产蜜桃| 国产日韩欧美在线精品| 18禁裸乳无遮挡动漫免费视频| 超碰97精品在线观看| 大话2 男鬼变身卡| 又黄又粗又硬又大视频| 99视频精品全部免费 在线| 人人澡人人妻人| 人妻 亚洲 视频| 国产精品国产三级国产av玫瑰| 欧美亚洲日本最大视频资源| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久成人av| 久久久久国产网址| 国产综合精华液| 男人添女人高潮全过程视频| 久久久久久伊人网av| 女人被躁到高潮嗷嗷叫费观| 国产精品国产三级专区第一集| 久久国内精品自在自线图片| 久久久久久久久久久免费av| 国产免费视频播放在线视频| 日本黄色日本黄色录像| 久久热在线av| 欧美变态另类bdsm刘玥| 黑人猛操日本美女一级片| 91成人精品电影| 成人国产麻豆网| 999精品在线视频| 免费观看在线日韩| 波多野结衣一区麻豆| 捣出白浆h1v1| 一级片免费观看大全| 国产男人的电影天堂91| 午夜老司机福利剧场| 欧美激情国产日韩精品一区| 黄片播放在线免费| 热99久久久久精品小说推荐| 亚洲欧洲精品一区二区精品久久久 | av在线观看视频网站免费| 国产老妇伦熟女老妇高清| 亚洲成色77777| 91在线精品国自产拍蜜月| 日韩成人伦理影院| 欧美性感艳星| 熟女av电影| 亚洲内射少妇av| 欧美日韩av久久| 国产av一区二区精品久久| 亚洲成人av在线免费| 美女国产视频在线观看| 夜夜骑夜夜射夜夜干| 久久久欧美国产精品| 少妇被粗大猛烈的视频| 国产精品久久久久成人av| 亚洲精华国产精华液的使用体验| 久久久久精品久久久久真实原创| 美国免费a级毛片| 国产av一区二区精品久久| 久久久久精品久久久久真实原创| 国产精品嫩草影院av在线观看| 亚洲性久久影院| 国产一区二区在线观看av| 熟女人妻精品中文字幕| 岛国毛片在线播放| 国产免费视频播放在线视频| 国产免费现黄频在线看| 毛片一级片免费看久久久久| 免费高清在线观看日韩| 精品久久久精品久久久| 免费观看在线日韩| 色婷婷久久久亚洲欧美| 中文天堂在线官网| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 国产 精品1| 人成视频在线观看免费观看| av女优亚洲男人天堂| 制服丝袜香蕉在线| 丝袜脚勾引网站| 亚洲精品aⅴ在线观看| 我要看黄色一级片免费的| 美女脱内裤让男人舔精品视频| 精品国产国语对白av| 亚洲精品美女久久久久99蜜臀 | 色94色欧美一区二区| 日日摸夜夜添夜夜爱| 看免费成人av毛片| 国产1区2区3区精品| 亚洲精品美女久久av网站| 久久精品久久精品一区二区三区| 日韩,欧美,国产一区二区三区| 日本av手机在线免费观看| 午夜福利,免费看| 国产精品女同一区二区软件| 97超碰精品成人国产| 2018国产大陆天天弄谢| 考比视频在线观看| 国产高清国产精品国产三级| 色5月婷婷丁香| 久久影院123| 亚洲国产精品999| 国产精品一区二区在线不卡| 亚洲国产精品999| 人人澡人人妻人| 久久精品国产亚洲av涩爱| 久久青草综合色| 亚洲在久久综合| 免费观看性生交大片5| 久久久久国产网址| 亚洲成人手机| 久久人妻熟女aⅴ| 婷婷色综合www| 亚洲欧美一区二区三区黑人 | 精品一品国产午夜福利视频| 免费av中文字幕在线| 另类亚洲欧美激情| 两个人看的免费小视频| 久久久久精品人妻al黑| 热re99久久国产66热| 色网站视频免费| 我要看黄色一级片免费的| 国产成人精品久久久久久| av不卡在线播放| 天堂中文最新版在线下载| 久久久久网色| 一级,二级,三级黄色视频| 国产日韩欧美视频二区| 欧美人与性动交α欧美精品济南到 | 中文乱码字字幕精品一区二区三区| 国产免费一级a男人的天堂| 欧美日韩av久久| 亚洲av国产av综合av卡| 国产69精品久久久久777片| 国产精品成人在线| 内地一区二区视频在线| 爱豆传媒免费全集在线观看| 免费人成在线观看视频色| 欧美另类一区| 久久av网站| 国产成人av激情在线播放| 精品少妇黑人巨大在线播放| 亚洲高清免费不卡视频| 国产免费福利视频在线观看| 亚洲欧美日韩另类电影网站| 免费观看a级毛片全部| 麻豆精品久久久久久蜜桃| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 我的女老师完整版在线观看| 日韩视频在线欧美| 成人18禁高潮啪啪吃奶动态图| 尾随美女入室| 国产精品久久久久久av不卡| 制服丝袜香蕉在线| 麻豆乱淫一区二区| 哪个播放器可以免费观看大片| 日韩制服丝袜自拍偷拍| 免费不卡的大黄色大毛片视频在线观看| 一个人免费看片子| www.熟女人妻精品国产 | 18禁在线无遮挡免费观看视频| 性高湖久久久久久久久免费观看| 一本久久精品| 成人毛片a级毛片在线播放| 国产黄色视频一区二区在线观看| 午夜av观看不卡| 久久精品国产亚洲av天美| www.色视频.com| 夜夜爽夜夜爽视频| 国产视频首页在线观看| 国产在线免费精品| 一本久久精品| 1024视频免费在线观看| 在线观看www视频免费| 亚洲精品,欧美精品| 丝袜在线中文字幕| 免费在线观看黄色视频的| 老女人水多毛片| 2018国产大陆天天弄谢| 亚洲成av片中文字幕在线观看 | 日本欧美国产在线视频| av国产久精品久网站免费入址| 少妇高潮的动态图| 国内精品宾馆在线| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久久久免| 日本vs欧美在线观看视频| 青春草视频在线免费观看| 亚洲国产精品一区三区| 波野结衣二区三区在线| a 毛片基地| 制服诱惑二区| 老女人水多毛片| 黑人猛操日本美女一级片| 在线精品无人区一区二区三| 午夜福利视频精品| 天堂中文最新版在线下载| 一级黄片播放器| av福利片在线| av在线老鸭窝| 国产视频首页在线观看| 久久久久久久久久久免费av| 久久久久网色| 成人黄色视频免费在线看| 老熟女久久久| 国产精品欧美亚洲77777| 国产探花极品一区二区| 中文字幕精品免费在线观看视频 | 成年av动漫网址| 精品午夜福利在线看| 亚洲,欧美精品.| 亚洲高清免费不卡视频| 母亲3免费完整高清在线观看 | 一级爰片在线观看| 亚洲情色 制服丝袜| 久久久久久久久久久免费av| 亚洲欧美中文字幕日韩二区| 日本av免费视频播放| av在线app专区| 欧美少妇被猛烈插入视频| 欧美激情极品国产一区二区三区 | 国产精品成人在线| 乱码一卡2卡4卡精品| 精品一区在线观看国产| 看十八女毛片水多多多| 中文乱码字字幕精品一区二区三区| av播播在线观看一区| av黄色大香蕉| 国产高清国产精品国产三级| 欧美最新免费一区二区三区| 久久精品夜色国产| 草草在线视频免费看| 亚洲精品久久成人aⅴ小说| 久久青草综合色| 成年美女黄网站色视频大全免费| 一边摸一边做爽爽视频免费| 搡女人真爽免费视频火全软件| 久久婷婷青草| 久久国产亚洲av麻豆专区| 男女下面插进去视频免费观看 | 七月丁香在线播放| 国产精品不卡视频一区二区| 咕卡用的链子| 午夜日本视频在线| 国产乱人偷精品视频| 国产片内射在线| 国产精品一区二区在线观看99| 黄色视频在线播放观看不卡| 又大又黄又爽视频免费| 久久狼人影院| 国产成人免费无遮挡视频| 人成视频在线观看免费观看| 18禁观看日本| 国国产精品蜜臀av免费| 久久女婷五月综合色啪小说| 少妇人妻精品综合一区二区| 少妇人妻久久综合中文| 熟女av电影| 美国免费a级毛片| 只有这里有精品99| 一个人免费看片子| 日韩成人av中文字幕在线观看| av.在线天堂| 亚洲av综合色区一区| av又黄又爽大尺度在线免费看| 99国产精品免费福利视频| 亚洲欧美清纯卡通| 如日韩欧美国产精品一区二区三区| 国产日韩欧美亚洲二区| 成人影院久久| 中文字幕另类日韩欧美亚洲嫩草| 麻豆精品久久久久久蜜桃| 亚洲av男天堂| a级片在线免费高清观看视频| 国语对白做爰xxxⅹ性视频网站| 超色免费av| 极品少妇高潮喷水抽搐| 黑人欧美特级aaaaaa片| 制服丝袜香蕉在线| 久久久久久久大尺度免费视频| 国产精品免费大片| av有码第一页| 亚洲欧美日韩另类电影网站| 一级毛片 在线播放| 18禁国产床啪视频网站| 国产亚洲精品久久久com| 免费观看a级毛片全部| 美女中出高潮动态图| 在线观看免费高清a一片| 久久久亚洲精品成人影院| 如何舔出高潮| 日韩av免费高清视频| 国产在线免费精品| 亚洲精品第二区| 最近最新中文字幕免费大全7| 在线观看一区二区三区激情| 在线观看国产h片| 精品熟女少妇av免费看| 蜜桃国产av成人99| 一级a做视频免费观看| 免费高清在线观看日韩| 久久精品久久久久久久性| 亚洲综合色网址| videosex国产| 黄网站色视频无遮挡免费观看| 黑人高潮一二区| 国产成人精品福利久久| 少妇 在线观看| 一本—道久久a久久精品蜜桃钙片| 精品少妇黑人巨大在线播放| 人妻少妇偷人精品九色| 日日撸夜夜添| 国产av一区二区精品久久| 看免费av毛片| 久久女婷五月综合色啪小说| 亚洲精品日本国产第一区| 国产精品一区二区在线观看99| 我要看黄色一级片免费的| 久久免费观看电影| 三级国产精品片| 观看av在线不卡| 黑人巨大精品欧美一区二区蜜桃 | 99久久人妻综合| 七月丁香在线播放| 午夜视频国产福利| 色视频在线一区二区三区| 久久久久久久久久人人人人人人| 精品一区二区三区视频在线| 国产黄频视频在线观看| 中文字幕制服av| 七月丁香在线播放| 午夜视频国产福利| 亚洲色图综合在线观看| 免费女性裸体啪啪无遮挡网站| 97精品久久久久久久久久精品| 欧美人与性动交α欧美精品济南到 | 久久人人爽人人爽人人片va| 久久av网站| 下体分泌物呈黄色| av免费在线看不卡| 国产一区二区三区综合在线观看 | 久久精品国产自在天天线| 国产av一区二区精品久久|