• <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| 尤物成人国产欧美一区二区三区| 久久精品久久精品一区二区三区| 亚洲伊人久久精品综合| 婷婷色麻豆天堂久久| 1000部很黄的大片| 婷婷色麻豆天堂久久| 国产亚洲5aaaaa淫片| 只有这里有精品99| 欧美日韩精品成人综合77777| 色综合亚洲欧美另类图片| 国产av不卡久久| 国产亚洲91精品色在线| 男女视频在线观看网站免费| 久久久久性生活片| 国产一区有黄有色的免费视频 | 国产精品国产三级专区第一集| 麻豆久久精品国产亚洲av| 亚洲精品久久久久久婷婷小说| 婷婷六月久久综合丁香| 久久人人爽人人爽人人片va| 中文精品一卡2卡3卡4更新| 欧美一区二区亚洲| 一个人免费在线观看电影| 国语对白做爰xxxⅹ性视频网站| 免费黄色在线免费观看| 又黄又爽又刺激的免费视频.| 尾随美女入室| 精品人妻偷拍中文字幕| 色吧在线观看| 乱人视频在线观看| 国产精品嫩草影院av在线观看| 欧美丝袜亚洲另类| 免费黄色在线免费观看| 国产91av在线免费观看| 亚洲av中文字字幕乱码综合| 欧美日韩一区二区视频在线观看视频在线 | 日日干狠狠操夜夜爽| 男插女下体视频免费在线播放| 99热这里只有是精品在线观看| 欧美极品一区二区三区四区| 亚洲,欧美,日韩| 亚洲国产高清在线一区二区三| 久久久午夜欧美精品| 97人妻精品一区二区三区麻豆| 亚洲人成网站在线观看播放| 亚洲av男天堂| 亚洲av二区三区四区| 亚洲人成网站高清观看| 亚洲欧美中文字幕日韩二区| 国产精品一及| 99久国产av精品| 久久久久网色| 91精品伊人久久大香线蕉| 国产大屁股一区二区在线视频| 久久久久九九精品影院| 麻豆国产97在线/欧美| 精品久久久久久久久亚洲| 网址你懂的国产日韩在线| 成人午夜精彩视频在线观看| 色网站视频免费| 成人毛片60女人毛片免费| av播播在线观看一区| 免费电影在线观看免费观看| 欧美性感艳星| 国产黄色视频一区二区在线观看| 亚洲精品日本国产第一区| 老司机影院成人| 亚洲成人久久爱视频| 久久精品国产亚洲av天美| 亚洲av中文字字幕乱码综合| 我的老师免费观看完整版| 欧美bdsm另类| 性色avwww在线观看| 极品少妇高潮喷水抽搐| 午夜福利高清视频| 搡老乐熟女国产| 非洲黑人性xxxx精品又粗又长| 中文字幕免费在线视频6| 人妻少妇偷人精品九色| 老女人水多毛片| 日本爱情动作片www.在线观看| 午夜福利视频精品| 久久久久精品性色| 免费在线观看成人毛片| 中文乱码字字幕精品一区二区三区 | 午夜福利成人在线免费观看| 搡老妇女老女人老熟妇| 插逼视频在线观看| 国产一级毛片七仙女欲春2| 日日啪夜夜撸| 2022亚洲国产成人精品| 爱豆传媒免费全集在线观看| 亚洲aⅴ乱码一区二区在线播放| 成人av在线播放网站| 亚洲国产精品国产精品| 欧美激情国产日韩精品一区| 2021天堂中文幕一二区在线观| 久久人人爽人人爽人人片va| 免费观看的影片在线观看| 久热久热在线精品观看| 国产亚洲精品av在线| ponron亚洲| av国产免费在线观看| 色哟哟·www| 久久久久九九精品影院| 日日啪夜夜爽| www.av在线官网国产| 永久免费av网站大全| 少妇丰满av| 成人漫画全彩无遮挡| 国产老妇伦熟女老妇高清| av线在线观看网站| 少妇高潮的动态图| 国产真实伦视频高清在线观看| 男女边吃奶边做爰视频| 亚洲精品国产av成人精品| 极品教师在线视频| 熟妇人妻久久中文字幕3abv| 欧美激情国产日韩精品一区| 日本色播在线视频| 一个人免费在线观看电影| 免费在线观看成人毛片| kizo精华| 人妻夜夜爽99麻豆av| 国产精品久久久久久精品电影| 国产黄色小视频在线观看| 水蜜桃什么品种好| 伊人久久国产一区二区| 99热这里只有是精品50| 在线a可以看的网站| 久久国产乱子免费精品| 校园人妻丝袜中文字幕| 久久精品久久久久久久性| 欧美成人午夜免费资源| 免费看光身美女| 国产高清不卡午夜福利| 97超视频在线观看视频| 日韩欧美精品免费久久| 亚洲精品一区蜜桃| 久久久欧美国产精品| 99re6热这里在线精品视频| 人妻夜夜爽99麻豆av| 一级毛片 在线播放| 精品久久国产蜜桃| 成年人午夜在线观看视频 | 秋霞伦理黄片| 夫妻性生交免费视频一级片| a级一级毛片免费在线观看| 成人一区二区视频在线观看| 99视频精品全部免费 在线| 精华霜和精华液先用哪个| 舔av片在线| av女优亚洲男人天堂| 精品一区二区三卡| 亚洲第一区二区三区不卡| 丰满人妻一区二区三区视频av| 久久97久久精品| 国产精品国产三级国产专区5o| 国产成人a区在线观看| 白带黄色成豆腐渣| 久久久久久久久久黄片| 美女cb高潮喷水在线观看| 亚洲人与动物交配视频| 人妻一区二区av| 在线观看av片永久免费下载| 日韩av免费高清视频| 亚洲最大成人手机在线| 日本一本二区三区精品| 九九在线视频观看精品| 日本猛色少妇xxxxx猛交久久| 成人无遮挡网站| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 草草在线视频免费看| 五月天丁香电影| 成人毛片a级毛片在线播放| 色网站视频免费| 色综合站精品国产| 亚洲av在线观看美女高潮| 午夜福利高清视频| 成人午夜高清在线视频| av在线老鸭窝| 2021天堂中文幕一二区在线观| 午夜福利在线观看吧| 秋霞伦理黄片| 亚洲欧美成人综合另类久久久| 日韩强制内射视频| 亚洲va在线va天堂va国产| 乱人视频在线观看| 视频中文字幕在线观看| 色综合亚洲欧美另类图片| 91久久精品国产一区二区成人| 国产精品伦人一区二区| 伦精品一区二区三区| 日日啪夜夜撸| 搡老乐熟女国产| 久久精品夜色国产| 国产激情偷乱视频一区二区| 久99久视频精品免费| 亚洲丝袜综合中文字幕| 肉色欧美久久久久久久蜜桃 | 又大又黄又爽视频免费| 国产精品一区www在线观看| 精品少妇黑人巨大在线播放| 亚洲欧洲日产国产| 国产乱人偷精品视频| 亚洲av中文av极速乱| 亚洲高清免费不卡视频| 亚洲国产日韩欧美精品在线观看| 免费不卡的大黄色大毛片视频在线观看 | 99re6热这里在线精品视频| 免费电影在线观看免费观看| 18禁动态无遮挡网站| 麻豆久久精品国产亚洲av| 啦啦啦中文免费视频观看日本| 亚洲人成网站高清观看| 亚洲欧美成人综合另类久久久| 国产精品一区www在线观看| 国产精品.久久久| 国产黄a三级三级三级人| 99久国产av精品| 永久免费av网站大全| 国精品久久久久久国模美| 欧美潮喷喷水| 99久久中文字幕三级久久日本| 久久精品久久久久久久性| 全区人妻精品视频| 欧美激情在线99| 联通29元200g的流量卡| 99热网站在线观看| 国产单亲对白刺激| 日韩中字成人| 国产精品嫩草影院av在线观看| 久久久久国产网址| 国产亚洲午夜精品一区二区久久 | 日韩国内少妇激情av| 五月玫瑰六月丁香| 国产乱来视频区| 成人亚洲精品一区在线观看 | 日韩av不卡免费在线播放| 久久午夜福利片| 天堂网av新在线| 亚洲欧美日韩无卡精品| 能在线免费观看的黄片| 欧美最新免费一区二区三区| 精品一区在线观看国产| 精品国内亚洲2022精品成人| 亚洲精品乱久久久久久| 亚洲久久久久久中文字幕| 久久久久久久久久成人| 国内揄拍国产精品人妻在线| 国产单亲对白刺激| 日本免费在线观看一区| 日韩大片免费观看网站| 成人二区视频| 亚洲精品,欧美精品| 免费观看的影片在线观看| 人妻制服诱惑在线中文字幕| 最近视频中文字幕2019在线8| 国产成人a区在线观看| 看非洲黑人一级黄片| 水蜜桃什么品种好| 两个人的视频大全免费| 最近手机中文字幕大全| 国产男女超爽视频在线观看| 久久精品久久精品一区二区三区| 成人一区二区视频在线观看| 亚洲自偷自拍三级| 国产 一区 欧美 日韩| 91av网一区二区| 狂野欧美白嫩少妇大欣赏| 国产精品久久久久久久电影| 国模一区二区三区四区视频| 黄色配什么色好看| 免费观看性生交大片5| a级一级毛片免费在线观看| 一级毛片 在线播放| 嘟嘟电影网在线观看| 久久久亚洲精品成人影院| 性色avwww在线观看| 一个人看视频在线观看www免费| 久久久久久国产a免费观看| 中文字幕亚洲精品专区| 国产一区二区三区av在线| 水蜜桃什么品种好| 91久久精品电影网| 男人爽女人下面视频在线观看| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 美女被艹到高潮喷水动态| 色综合亚洲欧美另类图片| 欧美 日韩 精品 国产| 天天躁夜夜躁狠狠久久av| 亚洲国产最新在线播放| 国产白丝娇喘喷水9色精品| 免费看美女性在线毛片视频| 久久久久久久久中文| 国产成年人精品一区二区| 国产久久久一区二区三区| 亚洲欧美精品专区久久| 日本欧美国产在线视频| 网址你懂的国产日韩在线| 噜噜噜噜噜久久久久久91| 丝袜喷水一区| 九九久久精品国产亚洲av麻豆| 非洲黑人性xxxx精品又粗又长| 日本-黄色视频高清免费观看| 国产精品久久久久久av不卡| 高清视频免费观看一区二区 | 亚洲性久久影院| 777米奇影视久久| 丝瓜视频免费看黄片| 久久99精品国语久久久| 亚洲av免费高清在线观看| 少妇高潮的动态图| 综合色丁香网| 97超碰精品成人国产| 黑人高潮一二区| 如何舔出高潮| 亚洲成人精品中文字幕电影| 狂野欧美白嫩少妇大欣赏| 国产免费福利视频在线观看| 欧美日韩国产mv在线观看视频 | 亚洲精品成人av观看孕妇| av福利片在线观看| 欧美激情在线99| 亚洲欧美成人精品一区二区| 亚洲av福利一区| 内地一区二区视频在线| 深夜a级毛片| 熟女人妻精品中文字幕| 午夜福利网站1000一区二区三区| 国内精品宾馆在线| 欧美三级亚洲精品| 性色avwww在线观看| 国内精品一区二区在线观看| 一级毛片久久久久久久久女| 国产伦理片在线播放av一区| 深爱激情五月婷婷| kizo精华| 偷拍熟女少妇极品色| 久久久欧美国产精品| 日本猛色少妇xxxxx猛交久久| 午夜精品一区二区三区免费看| 色吧在线观看| 免费少妇av软件| 99久久中文字幕三级久久日本| 国产91av在线免费观看| 亚洲美女视频黄频| 直男gayav资源| 免费播放大片免费观看视频在线观看| 亚洲精品乱码久久久v下载方式| 亚洲精品日韩av片在线观看| 蜜桃亚洲精品一区二区三区| 少妇被粗大猛烈的视频| 欧美日韩亚洲高清精品| 欧美日韩精品成人综合77777| 女人被狂操c到高潮| 麻豆国产97在线/欧美| 最近最新中文字幕大全电影3| 国产成人免费观看mmmm| 丝袜喷水一区| 国产成年人精品一区二区| 日日啪夜夜爽| 国产一区亚洲一区在线观看| 可以在线观看毛片的网站| 2018国产大陆天天弄谢| 成人欧美大片| 搡老乐熟女国产| 91精品一卡2卡3卡4卡| 日韩强制内射视频| 麻豆成人av视频| 欧美性感艳星| 久久久久国产网址| 美女内射精品一级片tv| 大陆偷拍与自拍| 亚洲婷婷狠狠爱综合网| 亚洲国产高清在线一区二区三| 国产亚洲av嫩草精品影院| 国产熟女欧美一区二区| 国产老妇伦熟女老妇高清| 国产女主播在线喷水免费视频网站 | 午夜视频国产福利| 我要看日韩黄色一级片| 不卡视频在线观看欧美| 69av精品久久久久久| 不卡视频在线观看欧美| 欧美成人a在线观看| 直男gayav资源| 亚洲经典国产精华液单| 国产一区二区三区综合在线观看 | 狂野欧美激情性xxxx在线观看| 欧美日韩亚洲高清精品| av在线观看视频网站免费| 久久久久精品性色| 插逼视频在线观看| 一级二级三级毛片免费看| 狠狠精品人妻久久久久久综合| 两个人的视频大全免费| 日日干狠狠操夜夜爽| 国产淫语在线视频| 搡女人真爽免费视频火全软件| 最近2019中文字幕mv第一页| 精品不卡国产一区二区三区| 久久综合国产亚洲精品| 国产淫语在线视频| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的 | 国产午夜精品久久久久久一区二区三区| 99久久精品国产国产毛片| 尾随美女入室| 一边亲一边摸免费视频| 亚洲美女搞黄在线观看| av播播在线观看一区| 非洲黑人性xxxx精品又粗又长| 国产亚洲5aaaaa淫片| 久久草成人影院| 亚洲av成人精品一二三区| 国产高清国产精品国产三级 | 成人美女网站在线观看视频| 高清av免费在线| 亚洲av中文字字幕乱码综合| 欧美日韩视频高清一区二区三区二| 少妇熟女aⅴ在线视频| 性插视频无遮挡在线免费观看| 国产有黄有色有爽视频| 人妻少妇偷人精品九色| 免费高清在线观看视频在线观看| 久久草成人影院| 草草在线视频免费看| 国产色爽女视频免费观看| 国产av不卡久久| 18+在线观看网站| 搡老妇女老女人老熟妇| 美女脱内裤让男人舔精品视频| 久久久色成人| 免费看a级黄色片| 一级av片app| 又大又黄又爽视频免费| 亚洲在线自拍视频| 国产成人福利小说| 精品少妇黑人巨大在线播放| 午夜老司机福利剧场| 欧美xxⅹ黑人| 99热全是精品| 国内揄拍国产精品人妻在线| 七月丁香在线播放| 亚洲欧洲日产国产| 亚洲,欧美,日韩| 91精品伊人久久大香线蕉| 最近中文字幕2019免费版| 国产片特级美女逼逼视频| 久久精品人妻少妇| 男人狂女人下面高潮的视频| 草草在线视频免费看| 一级爰片在线观看| 最近视频中文字幕2019在线8| 久久鲁丝午夜福利片| 午夜老司机福利剧场| 日韩欧美 国产精品| 成年免费大片在线观看| 久久久久精品久久久久真实原创| 真实男女啪啪啪动态图| 久久精品综合一区二区三区| 国内精品一区二区在线观看| 2021少妇久久久久久久久久久| 国内少妇人妻偷人精品xxx网站| 国产精品国产三级国产专区5o| 黄片wwwwww| 午夜久久久久精精品| 国产亚洲午夜精品一区二区久久 | av.在线天堂| 亚洲成人久久爱视频| 精品人妻视频免费看| 日日啪夜夜爽| 成人特级av手机在线观看| 久久鲁丝午夜福利片| 午夜免费激情av| ponron亚洲| 国产精品久久久久久久久免| 久久久久久久久久黄片| 亚洲国产日韩欧美精品在线观看| 大片免费播放器 马上看| 成人毛片60女人毛片免费| 日本黄色片子视频| 高清日韩中文字幕在线| 亚洲精华国产精华液的使用体验| 性插视频无遮挡在线免费观看| 国产免费一级a男人的天堂| 三级毛片av免费| av播播在线观看一区| 伊人久久国产一区二区| 男女国产视频网站| 精品酒店卫生间| 精品欧美国产一区二区三| 亚洲乱码一区二区免费版| 国产乱人偷精品视频| 亚洲国产欧美人成| 白带黄色成豆腐渣| 精品人妻熟女av久视频| av免费观看日本| 99热这里只有是精品50| 99九九线精品视频在线观看视频| 精品不卡国产一区二区三区| 熟女人妻精品中文字幕| 日本色播在线视频| 精品久久久久久久久久久久久| 亚洲av成人av| 欧美激情国产日韩精品一区| 国产亚洲av嫩草精品影院| 日韩精品青青久久久久久| 免费黄网站久久成人精品| 日韩av免费高清视频| 久久鲁丝午夜福利片| 午夜福利视频1000在线观看| 日韩av在线免费看完整版不卡| 女的被弄到高潮叫床怎么办| 国产国拍精品亚洲av在线观看| 尾随美女入室| 国产在视频线精品| 男女边摸边吃奶| 国产一级毛片在线| 免费电影在线观看免费观看| 国产午夜精品论理片| 亚洲精品国产av成人精品| 亚洲精品,欧美精品| 午夜久久久久精精品| 我要看日韩黄色一级片| 国产av码专区亚洲av| 麻豆精品久久久久久蜜桃| 亚洲av.av天堂| 亚洲精品一二三| 午夜老司机福利剧场| 国产亚洲5aaaaa淫片| videos熟女内射| 18禁在线播放成人免费| 免费观看av网站的网址| 亚洲精品456在线播放app| 国产男人的电影天堂91| 国产亚洲精品av在线| 午夜福利在线观看吧| 午夜福利在线在线| 直男gayav资源| 亚洲最大成人av| 九九在线视频观看精品| 欧美丝袜亚洲另类| 免费av毛片视频| 亚洲人成网站在线观看播放| 欧美成人一区二区免费高清观看| 亚洲真实伦在线观看| 精品人妻一区二区三区麻豆| 欧美一区二区亚洲| 亚洲欧美精品自产自拍| 亚洲精品中文字幕在线视频 | 欧美成人午夜免费资源| 国产成人福利小说| 国产黄片美女视频| 秋霞在线观看毛片| 搡女人真爽免费视频火全软件| 国产午夜精品一二区理论片| 亚洲婷婷狠狠爱综合网| 免费高清在线观看视频在线观看| 真实男女啪啪啪动态图| 亚洲精品久久午夜乱码| 国产高潮美女av| 亚洲精品乱码久久久v下载方式| 丰满乱子伦码专区| 精品国产露脸久久av麻豆 | 亚洲精品自拍成人| 亚洲av在线观看美女高潮| 国产精品精品国产色婷婷| 亚洲第一区二区三区不卡| 晚上一个人看的免费电影| 一级片'在线观看视频| 久久久久久久久久久免费av| 韩国av在线不卡| 亚洲av免费在线观看| 久久久久国产网址| 51国产日韩欧美| 亚洲伊人久久精品综合| 亚洲性久久影院| 国产片特级美女逼逼视频| 亚洲欧美成人综合另类久久久| 免费av观看视频| 插阴视频在线观看视频| 免费观看a级毛片全部| 99久久精品一区二区三区| 成人无遮挡网站| 亚洲熟女精品中文字幕| 午夜免费观看性视频| 一级毛片我不卡| av国产免费在线观看| 最近最新中文字幕大全电影3| 床上黄色一级片| 18禁动态无遮挡网站| 亚洲性久久影院| 18禁在线播放成人免费| 国产黄片美女视频| 老司机影院成人| 欧美日韩在线观看h| 蜜桃久久精品国产亚洲av| 免费av毛片视频| 又爽又黄a免费视频| av又黄又爽大尺度在线免费看| 啦啦啦韩国在线观看视频| 一二三四中文在线观看免费高清| 日韩av不卡免费在线播放| 免费黄频网站在线观看国产| 日本av手机在线免费观看| 中文天堂在线官网| 中文在线观看免费www的网站| 日本免费a在线| 中文字幕av成人在线电影| 国产高清三级在线|