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

    儲(chǔ)層巖石的三維逾滲模型研究

    2020-07-01 08:17:24黃紅亮王旭東鄭智冬
    關(guān)鍵詞:滲流集群介質(zhì)

    蘇 嶄,李 瑋,黃紅亮,王旭東,鄭智冬

    (1.東北石油大學(xué)石油工程學(xué)院,黑龍江大慶163318;2.渤海鉆探工程有限公司第五鉆井分公司,河北河間062465;3.中海油研究總院有限責(zé)任公司海洋石油高效開發(fā)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京100028;4.渤海鉆探第三鉆井工程公司,天津300280)

    儲(chǔ)層中流體的流通狀況是制定油田開發(fā)方案的重要依據(jù),它對(duì)提高采收率具有重要作用。裂縫性地層中的裂縫是流體流動(dòng)的主要通道,研究裂縫連通情況對(duì)儲(chǔ)層滲流特性的影響,一直是油田開采領(lǐng)域的一個(gè)難點(diǎn)。隨著時(shí)代的發(fā)展,計(jì)算機(jī)技術(shù)在油田中的應(yīng)用也愈加廣泛。許多在現(xiàn)實(shí)中難以實(shí)現(xiàn)的實(shí)驗(yàn)、難以驗(yàn)證的猜想,都可以用計(jì)算機(jī)進(jìn)行模擬。計(jì)算機(jī)模擬已成為除實(shí)驗(yàn)外解決問題最為有效的手段。

    隨著介質(zhì)中的孔隙被隨機(jī)堵塞,孔隙度減小,滲透率變差。當(dāng)孔隙度下降到某一臨界值Pc時(shí),介質(zhì)由可滲透轉(zhuǎn)變?yōu)椴豢蓾B透,這種現(xiàn)象被定義為逾滲。逾滲模型是一種描述流體在無(wú)序介質(zhì)中隨機(jī)流動(dòng)的數(shù)學(xué)模型。彭自正等[1]利用計(jì)算機(jī)建立了巖石的二維逾滲模型。馮增朝等[2]詳細(xì)研究了二維條件下帶裂縫介質(zhì)的逾滲機(jī)制,建立了二維逾滲規(guī)律研究模型。朱大勇等[3]基于二維逾滲模型,用相同尺寸的均勻網(wǎng)格劃分巖石,并使每個(gè)連接節(jié)點(diǎn)的鍵具有相同的連通概率,研究逾滲模型中滲透概率遞推矩陣??抵厩诘萚4]基于顯微CT技術(shù)研究了三維孔隙結(jié)構(gòu)的逾滲規(guī)律。以上研究運(yùn)用不同方法建立了逾滲模型,從不同角度研究逾滲規(guī)律?,F(xiàn)實(shí)中的逾滲發(fā)生在三維物體中,二維模型模擬結(jié)果不夠直觀、精確。目前三維逾滲模型的建模方式較少,對(duì)逾滲規(guī)律的認(rèn)識(shí)也不夠深入。

    本文運(yùn)用MATLAB建立三維逾滲模型,將儲(chǔ)層巖石理想地看作由均勻網(wǎng)格劃分的規(guī)則多孔介質(zhì)。與傳統(tǒng)的研究方法不同,它從多孔介質(zhì)的微觀角度出發(fā),將其視為鍵和節(jié)點(diǎn)組成的網(wǎng)絡(luò)體。此模型能夠模擬裂縫連通團(tuán)的分布情況,得出均勻多孔介質(zhì)的逾滲規(guī)律;描述巖石的形狀、大小對(duì)其內(nèi)部流體流動(dòng)的影響,對(duì)高效開發(fā)油藏具有重要的意義。

    1 三維逾滲模型的建立

    1.1 構(gòu)建模型的思路

    逾滲理論研究的是多孔介質(zhì)由于內(nèi)部連接程度的改變所引起的滲透狀態(tài)的改變,即多孔介質(zhì)由不滲透到滲透的過程。裂縫被視為流體流通的通道,裂縫之間的相互連接使其“叢集”最終達(dá)到“貫通”[5]。建立三維逾滲模型,首先將所研究的儲(chǔ)層巖石理想化為一個(gè)規(guī)則的長(zhǎng)方體,然后用均勻網(wǎng)格劃分長(zhǎng)方體。將儲(chǔ)層中的裂縫視為逾滲網(wǎng)格中的節(jié)點(diǎn),裂縫之間的連接情況用逾滲網(wǎng)格中的鍵來(lái)表示。每個(gè)鍵都有一定的“連通概率(P)”,被鍵連接到一起的節(jié)點(diǎn)構(gòu)成一個(gè)簇,即逾滲集團(tuán)。

    由逾滲的定義可知,流體在多孔介質(zhì)中流動(dòng)時(shí),若介質(zhì)被隨機(jī)堵塞,一旦被堵塞的比例達(dá)到某個(gè)值,多孔介質(zhì)就會(huì)由可滲透轉(zhuǎn)變?yōu)椴豢蓾B透。也就是說,存在著一個(gè)臨界值,使模型是否可滲透的狀態(tài)發(fā)生改變,這個(gè)臨界值被稱為逾滲閾值。

    1.2 構(gòu)建模型的方法

    逾滲主要分為鍵逾滲和點(diǎn)逾滲。鍵逾滲是指所有的節(jié)點(diǎn)都是通的,每條鍵以概率P連通;點(diǎn)逾滲是指所有的鍵都是連通的,節(jié)點(diǎn)有可能通有可能不通,通過的點(diǎn)叫作被占據(jù)的點(diǎn)。本文中的模型以鍵逾滲理論為基礎(chǔ)。

    三維點(diǎn)陣的節(jié)點(diǎn)之間的鍵有兩種狀態(tài):連通或不連通,每個(gè)鍵連通的概率是P,不連通的概率是(1-P),每個(gè)鍵的連通狀況與相鄰的鍵無(wú)關(guān)。在三維點(diǎn)陣上,每個(gè)節(jié)點(diǎn)周圍有六個(gè)節(jié)點(diǎn)與之相鄰,假設(shè)每個(gè)節(jié)點(diǎn)代表相同的裂縫,可以認(rèn)為每個(gè)簇的大小僅由內(nèi)部的節(jié)點(diǎn)數(shù)量確定,這樣就建立了三維逾滲模型[6]。當(dāng)連通概率達(dá)到某個(gè)值時(shí)(即逾滲閾值),三維點(diǎn)陣上就會(huì)形成使模型可滲透的簇,即滲流集群。

    假設(shè)模型中鍵的連通概率(P)取值0.01~0.90,模型的長(zhǎng)、寬、高分別用n1、n2、n3表示。設(shè)計(jì)五種不同尺寸的模型N1-N5進(jìn)行模擬,如圖1所示。其中,N1、N2、N3的尺寸等比例增大,N2、N4、N5的節(jié)點(diǎn)數(shù)相近,幾何形態(tài)分別為長(zhǎng)方體、正方體和扁長(zhǎng)方體。

    圖1 五種不同尺寸模型Fig.1 Five different size models

    隨著連通概率(P)的改變,N1、N2、N3模型的簇的分布情況如圖2所示,其中藍(lán)點(diǎn)代表節(jié)點(diǎn),黑色、綠色、紅色的線分別代表節(jié)點(diǎn)之間x、y、z三個(gè)方向的鍵。模型能較直觀地反映多孔介質(zhì)的內(nèi)部連接情況,還能得到不同連通概率下,每個(gè)模型中各個(gè)尺寸的簇的數(shù)量和最大簇。

    圖2 不同尺寸模型的簇分布情況Fig.2 Cluster distribution of different size models

    由不同尺寸的簇的數(shù)量可以計(jì)算出孔隙率和平均簇,計(jì)算孔隙率的公式為:

    式中,Q 為孔隙率;s為簇的大??;Smax為最大簇;n(s)為簇大小為s的簇的個(gè)數(shù);N為模型總節(jié)點(diǎn)數(shù)。

    計(jì)算平均簇的公式為:

    式中,Sav為平均簇。

    相互連接的節(jié)點(diǎn)形成的集團(tuán)稱為簇,一個(gè)簇至少包含兩個(gè)節(jié)點(diǎn),節(jié)點(diǎn)數(shù)最多的簇為最大簇[7]。逾滲概率為最大簇節(jié)點(diǎn)數(shù)與總節(jié)點(diǎn)數(shù)的比值,它可以作為判斷多孔介質(zhì)滲透性優(yōu)劣的依據(jù),是逾滲問題的一個(gè)重要概念。

    計(jì)算逾滲概率的公式為:

    式中,Yg為逾滲概率;Nmax為最大簇節(jié)點(diǎn)數(shù)。

    2 模型模擬結(jié)果的分析

    2.1 連通概率對(duì)總簇?cái)?shù)的影響

    模型可以模擬得出多孔介質(zhì)中的總簇?cái)?shù)。在不同尺寸的模型中,改變連通概率進(jìn)行模擬,其對(duì)應(yīng)的總簇?cái)?shù)如表1所示,兩者的關(guān)系如圖3所示。

    表1 總簇?cái)?shù)Table 1 Total number of clusters

    圖3 不同尺寸模型總簇?cái)?shù)與連通概率的關(guān)系Fig.3 Relationship between total cluster number and P of different size models

    由圖3可知,當(dāng)連通概率較小時(shí),相互連接的節(jié)點(diǎn)少,聚集程度低,難以形成簇。隨著連通概率增加、模型尺寸增大,總簇?cái)?shù)增加得越快。此時(shí)節(jié)點(diǎn)之間快速連接形成簇,雖然已經(jīng)形成的簇有可能會(huì)連接成一個(gè)新簇,使總簇?cái)?shù)減少,但總體來(lái)說,總簇?cái)?shù)還是呈增加趨勢(shì)。此時(shí)連通概率繼續(xù)增加,已經(jīng)形成的簇相互連接,總簇?cái)?shù)的增長(zhǎng)變緩[7]。當(dāng)連通概率達(dá)到一定值時(shí),總簇?cái)?shù)開始隨連通概率的增加而減少,這是因?yàn)榇刂g連接的趨勢(shì)初始階段大于單個(gè)節(jié)點(diǎn)之間連接的趨勢(shì)。

    N1、N2、N3模型尺寸不同,其總簇?cái)?shù)與連通概率的關(guān)系曲線有較大差別,但曲線峰值對(duì)應(yīng)的連通概率相同,均為0.17。N2、N4、N5的曲線基本重合,說明總簇?cái)?shù)只與模型的尺寸有關(guān),與模型的形狀無(wú)關(guān)。

    2.2 孔隙率與逾滲概率的關(guān)系

    孔隙率可以用式(1)計(jì)算,當(dāng)可滲透的多孔介質(zhì)被隨機(jī)堵塞時(shí),多孔介質(zhì)的孔隙率下降,多孔介質(zhì)由可滲透轉(zhuǎn)變?yōu)椴豢蓾B透[8]。用孔隙率可以粗略評(píng)價(jià)介質(zhì)的滲透性。不同連通概率下,N1-N5模型的孔隙率如圖4所示。

    圖4 不同尺寸模型孔隙率與連通概率的關(guān)系Fig.4 Relationship between porosity and P of different size models

    對(duì)不同尺寸的模型,在不同連通概率下進(jìn)行模擬,得到最大簇,運(yùn)用式(3)計(jì)算逾滲概率。同一連通概率下,孔隙率和逾滲概率的關(guān)系如圖5所示。

    圖5 孔隙率和逾滲概率的關(guān)系Fig.5 Relationship between porosity and percolation probability

    由圖5可知,連通概率小于0.20時(shí),孔隙率和逾滲概率之間的關(guān)系沒有規(guī)律性。逾滲概率的大小反映了多孔介質(zhì)內(nèi)部最大簇的分布程度,可避免因過多統(tǒng)計(jì)對(duì)滲透性無(wú)貢獻(xiàn)的“小級(jí)別死端及封閉裂縫”而造成的對(duì)多孔介質(zhì)滲透性優(yōu)劣的誤判[8]。逾滲概率可以對(duì)多孔介質(zhì)的滲透性進(jìn)行判斷,比孔隙率更為準(zhǔn)確。

    2.3 連通概率對(duì)逾滲概率的影響

    不同連通概率下,N1-N5模型的逾滲概率如圖6所示。由圖6可知,不同尺寸、不同形狀的模型,其逾滲概率隨連通概率的變化呈現(xiàn)相同的趨勢(shì)。當(dāng)連通概率較小時(shí),逾滲概率增加比較緩慢,對(duì)N1、N5模型來(lái)說,當(dāng)連通概率超過第一臨界值0.30時(shí),逾滲概率突然大幅度增加,當(dāng)連通概率超過第二臨界值0.70時(shí),逾滲概率再次開始緩慢增加并最終趨近于1。

    圖6 不同尺寸模型逾滲概率與連通概率的關(guān)系Fig.6 Relationship between percolation probability and P of different size models

    對(duì)N2、N3、N4模型來(lái)說,當(dāng)連通概率超過第一臨界值0.25時(shí),逾滲概率突然大幅度增加,當(dāng)連通概率超過第二臨界值0.60時(shí),逾滲概率再次開始緩慢增加并最終趨近于1。

    當(dāng)連通概率小于第一臨界值時(shí),單個(gè)節(jié)點(diǎn)之間連接的趨勢(shì)較小,逾滲概率的增長(zhǎng)率較?。划?dāng)連通概率在第一、第二臨界值之間時(shí),簇與簇之間迅速連接形成更大的簇;當(dāng)連通概率大于第二臨界值時(shí),節(jié)點(diǎn)與簇的連接狀態(tài)比較穩(wěn)定,逾滲概率增加較為緩慢,最終趨向于 1[9‐10]。

    2.4 連通概率對(duì)滲流集群的影響

    在實(shí)際應(yīng)用中,用逾滲概率判斷多孔介質(zhì)是否可滲透存在一定誤差。本模型提出了滲流集群的概念,滲流集群是可以使多孔介質(zhì)滲透的簇,它可以準(zhǔn)確判斷介質(zhì)是否可滲透,精準(zhǔn)捕捉逾滲閾值。本模型中,逾滲閾值指的是當(dāng)多孔介質(zhì)由不滲透轉(zhuǎn)變?yōu)闈B透時(shí)對(duì)應(yīng)的連通概率,即滲流集群≥1時(shí)對(duì)應(yīng)的連通概率。不同連通概率下,N1-N5模型的滲流集群如圖7所示。滲流集群等于0時(shí),模型不可滲透;滲流集群在0~1時(shí),模型有滲透的可能,且滲透的可能性隨著滲流集群數(shù)值的增大而增大;滲流集群等于或大于1時(shí),模型一定可滲透。

    圖7 不同尺寸模型滲流集群數(shù)與連通概率的關(guān)系Fig.7 Relationship between the number of seepage clusters and the P of different size models

    由圖7可以看出,N1-N5模型的滲流集群開始大于0時(shí)對(duì)應(yīng)的連通概率分別是0.05、0.10、0.17、0.10、0.07,即五個(gè)模型各自的逾滲閾值。對(duì)于N1、N2、N3來(lái)說,模型尺寸越大,逾滲閾值就越大;對(duì)于N2、N4、N5來(lái)說,N5的幾何形狀較為特殊,有一個(gè)較小的面,因此其逾滲閾值小于N2、N4。

    五條曲線的最大值都不相同,N3最大,N1最小,N2、N4、N5相近可看作一類??芍P统叽缭酱螅?jié)點(diǎn)越多,可形成滲流集群的上限越高,反映到圖像中就是曲線的最大值越大。

    隨著連通概率的增加,不管是何種尺寸、何種形狀的模型,滲流集群的數(shù)量最終都會(huì)趨向于1。因?yàn)楫?dāng)連通概率達(dá)到一定值時(shí),模型內(nèi)的節(jié)點(diǎn)連接狀況趨于穩(wěn)定,最終會(huì)形成一個(gè)滲流集群。

    在圖線波動(dòng)階段,尺寸越小的模型波動(dòng)得越劇烈。因?yàn)槌叽缧〉哪P凸?jié)點(diǎn)少,由連通概率所導(dǎo)致的連接不確定性比較明顯[11‐12]。

    2.5 平均簇與連通概率之間的關(guān)系

    對(duì)N1-N5模型進(jìn)行模擬,可以得出各個(gè)尺寸的簇的個(gè)數(shù),將得到的數(shù)據(jù)代入式(2),可計(jì)算平均簇,繪制圖8。由圖8可以看出,連通概率小于0.5時(shí),平均簇比較小且增長(zhǎng)緩慢;連通概率在0.5~0.8時(shí),平均簇的增長(zhǎng)率突然變大,這是因?yàn)榇嘏c簇之間迅速連接,使簇的尺寸變大,簇的個(gè)數(shù)變少;當(dāng)連通概率大于0.8時(shí),節(jié)點(diǎn)之間的連接狀況趨于穩(wěn)定,平均簇大小的增長(zhǎng)率有所緩和。

    2.6 對(duì)錯(cuò)誤逾滲規(guī)律的否定

    對(duì)N1-N5模型在相同的連通概率下多次模擬,可以明確一些容易誤判的逾滲規(guī)律。

    2.6.1 最大簇與平均簇的關(guān)系 在同一模型、同一連通概率下多次模擬,得到圖9所示結(jié)果。

    圖8 模型平均簇大小與連通概率的關(guān)系Fig.8 Relationship between model average cluster size and P

    圖9 最大簇與平均簇間的關(guān)系Fig.9 Relationship between the largest cluster and the average cluster size

    由圖9可以看出,平均簇和最大簇的大小隨著連通概率的增加總體呈增加趨勢(shì)。在同一模型,同一連通概率下多次模擬,可以得知最大簇和平均簇之間的關(guān)系沒有明顯的規(guī)律性。

    2.6.2 最大簇與總簇?cái)?shù)的關(guān)系 在同一模型、同一連通概率下多次模擬,得到的結(jié)果如圖10所示。由圖10可以看出,最大簇的大小隨著連通概率的增加總體呈增加趨勢(shì),總簇?cái)?shù)隨著連通概率的增加總體呈減少趨勢(shì)。在同一模型,同一連通概率下多次模擬,可以得知最大簇和總簇?cái)?shù)之間沒有明顯的規(guī)律性。

    圖10 最大簇與總簇?cái)?shù)的關(guān)系Fig.10 Relationship between the largest cluster and the total number of clusters

    3 結(jié) 論

    (1)本文所建立的三維逾滲模型能夠模擬多孔介質(zhì)中由于連通概率變化所引起的長(zhǎng)程聯(lián)結(jié)性突變。它能模擬裂縫連通團(tuán)的分布情況,得出均勻多孔介質(zhì)的逾滲規(guī)律、逾滲閾值,并判斷其是否可滲透,模型模擬的結(jié)果與其它模型模擬結(jié)果及實(shí)驗(yàn)結(jié)果基本一致并有了一定的優(yōu)化。

    (2)在連通概率從0開始增加的初期,隨著連通概率增加,單個(gè)節(jié)點(diǎn)之間相互連接形成簇,總簇?cái)?shù)增加。當(dāng)連通概率增加到0.17后,總簇?cái)?shù)開始減少??偞?cái)?shù)的變化規(guī)律只與介質(zhì)的尺寸有關(guān),與介質(zhì)的形狀無(wú)關(guān)。

    (3)孔隙率、逾滲概率都可以作為判斷多孔介質(zhì)滲透性的依據(jù),逾滲概率可避免過多統(tǒng)計(jì)對(duì)滲透性無(wú)貢獻(xiàn)的“小級(jí)別死端及封閉裂縫”,因此更為準(zhǔn)確。當(dāng)連通概率較小時(shí),逾滲概率隨著連通概率的增加而緩慢增加,當(dāng)連通概率超過第一臨界值時(shí),逾滲概率的增加幅度突然增大,當(dāng)連通概率超過第二臨界值時(shí),逾滲概率再次開始緩慢增加且最終趨近于1。

    (4)滲流集群可以準(zhǔn)確地判定多孔介質(zhì)是否可滲透,精準(zhǔn)捕捉逾滲閾值。N1-N5模型的滲流集群開始大于0時(shí)對(duì)應(yīng)的連通概率分別是0.05、0.10、0.17、0.10、0.07,即五個(gè)模型各自的逾滲閾值。逾滲閾值的大小與介質(zhì)的尺寸和形狀都有關(guān),尺寸越大,逾滲閾值就越大;當(dāng)尺寸相同時(shí),若介質(zhì)有一個(gè)較小的面,逾滲閾值就會(huì)較小。

    (5)最大簇與平均簇、最大簇與總簇?cái)?shù)的關(guān)系之間沒有明顯的規(guī)律性。

    猜你喜歡
    滲流集群介質(zhì)
    信息交流介質(zhì)的演化與選擇偏好
    淬火冷卻介質(zhì)在航空工業(yè)的應(yīng)用
    海上小型無(wú)人機(jī)集群的反制裝備需求與應(yīng)對(duì)之策研究
    一種無(wú)人機(jī)集群發(fā)射回收裝置的控制系統(tǒng)設(shè)計(jì)
    電子制作(2018年11期)2018-08-04 03:25:40
    Python與Spark集群在收費(fèi)數(shù)據(jù)分析中的應(yīng)用
    勤快又呆萌的集群機(jī)器人
    簡(jiǎn)述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計(jì)算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    考慮中間介質(zhì)換熱的廠際熱聯(lián)合
    尾礦壩滲流計(jì)算及排滲設(shè)計(jì)
    金屬礦山(2013年6期)2013-03-11 16:54:05
    最近2019中文字幕mv第一页| 国产欧美日韩精品一区二区| 婷婷色av中文字幕| 女人十人毛片免费观看3o分钟| 亚洲精品影视一区二区三区av| 卡戴珊不雅视频在线播放| 综合色丁香网| 男人的好看免费观看在线视频| 国产黄色视频一区二区在线观看 | 欧美日本视频| 人人妻人人看人人澡| 欧美日韩国产亚洲二区| 国产精品国产高清国产av| 精品国产一区二区三区久久久樱花 | 国产真实伦视频高清在线观看| 精品人妻偷拍中文字幕| 自拍偷自拍亚洲精品老妇| 婷婷色麻豆天堂久久 | 国国产精品蜜臀av免费| 欧美一区二区国产精品久久精品| 亚洲av福利一区| 一级毛片电影观看 | 欧美成人a在线观看| 精品久久久久久久末码| 国产精华一区二区三区| 久久综合国产亚洲精品| 国内揄拍国产精品人妻在线| 精品人妻偷拍中文字幕| 成人亚洲欧美一区二区av| 看片在线看免费视频| 亚洲国产成人一精品久久久| 午夜爱爱视频在线播放| av视频在线观看入口| 日本猛色少妇xxxxx猛交久久| 最新中文字幕久久久久| 在线免费观看的www视频| 午夜福利在线观看免费完整高清在| 日本熟妇午夜| 亚洲在线观看片| 人妻系列 视频| 18禁在线无遮挡免费观看视频| 丰满乱子伦码专区| 精品国产三级普通话版| 一级黄色大片毛片| 久久久精品大字幕| 大又大粗又爽又黄少妇毛片口| 国产精华一区二区三区| 国内揄拍国产精品人妻在线| 精华霜和精华液先用哪个| 午夜福利在线观看免费完整高清在| 亚洲电影在线观看av| 亚洲精品久久久久久婷婷小说 | 九色成人免费人妻av| 精品久久国产蜜桃| 青青草视频在线视频观看| 亚洲婷婷狠狠爱综合网| 99久久精品国产国产毛片| 超碰97精品在线观看| 国产又黄又爽又无遮挡在线| 波多野结衣巨乳人妻| 日韩在线高清观看一区二区三区| 岛国毛片在线播放| 不卡视频在线观看欧美| 免费黄网站久久成人精品| 久久99热6这里只有精品| 午夜免费激情av| 亚洲欧美日韩卡通动漫| 中文字幕熟女人妻在线| 男人舔女人下体高潮全视频| 中文天堂在线官网| 女人久久www免费人成看片 | 少妇被粗大猛烈的视频| 人人妻人人澡欧美一区二区| 国产精品久久久久久久电影| 特级一级黄色大片| 亚洲av免费高清在线观看| 精品免费久久久久久久清纯| 亚洲内射少妇av| 亚洲美女视频黄频| 久久欧美精品欧美久久欧美| 欧美精品一区二区大全| 色网站视频免费| 看黄色毛片网站| 亚洲在线观看片| 级片在线观看| 又爽又黄无遮挡网站| 亚洲欧美日韩高清专用| 色播亚洲综合网| 激情 狠狠 欧美| 精品久久久久久久久av| 亚洲精品乱久久久久久| 久久这里只有精品中国| av在线蜜桃| 国产老妇女一区| 亚洲色图av天堂| 亚洲国产精品合色在线| 国产伦一二天堂av在线观看| 国产免费又黄又爽又色| 亚洲av电影不卡..在线观看| 少妇熟女欧美另类| 国产精品一区www在线观看| 中文乱码字字幕精品一区二区三区 | 久久精品熟女亚洲av麻豆精品 | 日日摸夜夜添夜夜爱| 成人午夜高清在线视频| 一个人免费在线观看电影| 综合色av麻豆| 在线观看美女被高潮喷水网站| 国产老妇伦熟女老妇高清| 91精品国产九色| 啦啦啦观看免费观看视频高清| 插逼视频在线观看| 草草在线视频免费看| 淫秽高清视频在线观看| 人妻少妇偷人精品九色| 99久久人妻综合| 看黄色毛片网站| 日韩国内少妇激情av| 中文字幕av在线有码专区| 亚洲国产精品合色在线| 久久精品91蜜桃| 亚洲激情五月婷婷啪啪| 丰满人妻一区二区三区视频av| av在线天堂中文字幕| 国产精品美女特级片免费视频播放器| 国产乱来视频区| 51国产日韩欧美| 简卡轻食公司| 简卡轻食公司| 狂野欧美白嫩少妇大欣赏| 久久久久网色| 一本一本综合久久| 亚洲国产最新在线播放| 最后的刺客免费高清国语| 亚洲国产最新在线播放| 亚洲va在线va天堂va国产| 亚洲欧美清纯卡通| 欧美zozozo另类| 亚洲av熟女| 一二三四中文在线观看免费高清| 亚洲aⅴ乱码一区二区在线播放| 我要看日韩黄色一级片| 丰满乱子伦码专区| 亚洲欧洲日产国产| 成年女人永久免费观看视频| 亚洲国产欧美在线一区| 国产男人的电影天堂91| 有码 亚洲区| 丰满人妻一区二区三区视频av| 国产精品国产三级国产av玫瑰| 亚洲精品乱久久久久久| 国产成年人精品一区二区| 国产精品精品国产色婷婷| 最近中文字幕高清免费大全6| 国产亚洲最大av| av在线天堂中文字幕| 美女xxoo啪啪120秒动态图| 日本av手机在线免费观看| 欧美色视频一区免费| 精品久久久久久久久久久久久| 亚洲在线观看片| 久久久欧美国产精品| 精品久久久久久久人妻蜜臀av| 少妇的逼水好多| 国产在线男女| 日日摸夜夜添夜夜添av毛片| 一个人看视频在线观看www免费| 国产白丝娇喘喷水9色精品| 人人妻人人澡欧美一区二区| av播播在线观看一区| 国产av在哪里看| 国产真实伦视频高清在线观看| 亚洲欧洲国产日韩| 久久草成人影院| 级片在线观看| 成人高潮视频无遮挡免费网站| 欧美又色又爽又黄视频| 婷婷六月久久综合丁香| av视频在线观看入口| 久99久视频精品免费| 国产色婷婷99| 美女脱内裤让男人舔精品视频| 亚洲欧美日韩无卡精品| 亚洲av成人av| 日本猛色少妇xxxxx猛交久久| 国产不卡一卡二| 3wmmmm亚洲av在线观看| 久久精品国产自在天天线| 亚洲美女视频黄频| 变态另类丝袜制服| 中国国产av一级| 欧美一区二区亚洲| 午夜福利视频1000在线观看| 日韩亚洲欧美综合| 亚洲国产精品成人久久小说| 99热这里只有是精品在线观看| 色视频www国产| 久久精品国产亚洲av天美| 国产亚洲精品av在线| 日韩av不卡免费在线播放| 亚洲自偷自拍三级| 午夜免费激情av| 久久久精品大字幕| 午夜福利高清视频| 亚洲av福利一区| 最近2019中文字幕mv第一页| 欧美成人精品欧美一级黄| 国产综合懂色| 国产精品电影一区二区三区| 国产高清不卡午夜福利| 亚洲av二区三区四区| 1000部很黄的大片| 哪个播放器可以免费观看大片| 久久久久久久久大av| 一级毛片电影观看 | 乱人视频在线观看| 热99re8久久精品国产| 成人av在线播放网站| 99在线视频只有这里精品首页| 亚洲人与动物交配视频| 老司机福利观看| 黄色欧美视频在线观看| 国产又色又爽无遮挡免| 特级一级黄色大片| 精品不卡国产一区二区三区| 一区二区三区高清视频在线| 汤姆久久久久久久影院中文字幕 | 色综合色国产| 在线天堂最新版资源| 精品一区二区免费观看| 嫩草影院入口| 直男gayav资源| 亚洲国产精品合色在线| 寂寞人妻少妇视频99o| 99久久精品热视频| 中国国产av一级| 国产老妇伦熟女老妇高清| 欧美高清性xxxxhd video| 久久久精品大字幕| 熟妇人妻久久中文字幕3abv| 久久精品91蜜桃| 亚洲美女视频黄频| 一个人观看的视频www高清免费观看| 高清毛片免费看| 精品酒店卫生间| 久久婷婷人人爽人人干人人爱| 亚洲精品国产av成人精品| 高清毛片免费看| 高清视频免费观看一区二区 | 日韩,欧美,国产一区二区三区 | 99久国产av精品国产电影| 亚洲国产精品成人久久小说| 久久久久性生活片| 中文字幕免费在线视频6| 一级毛片aaaaaa免费看小| 精品一区二区三区人妻视频| 日韩人妻高清精品专区| 人人妻人人澡人人爽人人夜夜 | 久久久久性生活片| 直男gayav资源| 91狼人影院| 欧美精品一区二区大全| a级毛片免费高清观看在线播放| 一级黄片播放器| 亚洲最大成人中文| 七月丁香在线播放| 久久久久网色| kizo精华| 国产一区二区在线观看日韩| 久久韩国三级中文字幕| 亚洲一区高清亚洲精品| 91精品一卡2卡3卡4卡| 国产高潮美女av| 我要搜黄色片| 欧美高清性xxxxhd video| 亚洲成人av在线免费| 亚洲精品成人久久久久久| 日日摸夜夜添夜夜爱| 亚洲人成网站高清观看| 午夜日本视频在线| 色播亚洲综合网| 联通29元200g的流量卡| 国产精品福利在线免费观看| 免费大片18禁| 极品教师在线视频| 免费av不卡在线播放| 欧美成人a在线观看| 久久精品国产99精品国产亚洲性色| 一级毛片aaaaaa免费看小| 国产成人精品久久久久久| 在线免费观看的www视频| 久久精品久久精品一区二区三区| 中文字幕免费在线视频6| 高清午夜精品一区二区三区| 精品久久国产蜜桃| 2021少妇久久久久久久久久久| 听说在线观看完整版免费高清| 99久久无色码亚洲精品果冻| 青春草视频在线免费观看| 看非洲黑人一级黄片| 简卡轻食公司| 黄色日韩在线| 91久久精品国产一区二区三区| 亚洲欧美精品综合久久99| 午夜福利网站1000一区二区三区| 国产高清有码在线观看视频| 亚洲在线自拍视频| 国产精品无大码| 国产亚洲精品久久久com| av在线亚洲专区| 中文欧美无线码| 久久久久网色| 一级黄色大片毛片| 中国国产av一级| 中文在线观看免费www的网站| av在线观看视频网站免费| 日本爱情动作片www.在线观看| 欧美一级a爱片免费观看看| 成人国产麻豆网| 非洲黑人性xxxx精品又粗又长| 欧美日韩综合久久久久久| 床上黄色一级片| 直男gayav资源| 国产精品无大码| 乱系列少妇在线播放| 在线观看66精品国产| 我要搜黄色片| 联通29元200g的流量卡| 日本色播在线视频| 最近手机中文字幕大全| 中国美白少妇内射xxxbb| 色综合亚洲欧美另类图片| 麻豆乱淫一区二区| 亚洲第一区二区三区不卡| 午夜亚洲福利在线播放| 国产 一区 欧美 日韩| 国产精品99久久久久久久久| 亚洲精品成人久久久久久| 在现免费观看毛片| 欧美激情国产日韩精品一区| 少妇被粗大猛烈的视频| 免费黄色在线免费观看| 亚洲三级黄色毛片| 一级爰片在线观看| 亚洲欧美精品专区久久| 国内精品宾馆在线| 久久鲁丝午夜福利片| 亚州av有码| 嘟嘟电影网在线观看| 午夜福利在线观看免费完整高清在| 我要搜黄色片| 大又大粗又爽又黄少妇毛片口| 26uuu在线亚洲综合色| 久99久视频精品免费| 亚洲最大成人中文| 少妇的逼好多水| 1024手机看黄色片| 老女人水多毛片| h日本视频在线播放| 日本猛色少妇xxxxx猛交久久| 天美传媒精品一区二区| 色尼玛亚洲综合影院| 国产精品福利在线免费观看| 男人的好看免费观看在线视频| 国产伦理片在线播放av一区| 五月玫瑰六月丁香| 久久久久久久国产电影| 22中文网久久字幕| 18禁裸乳无遮挡免费网站照片| 国产午夜精品一二区理论片| 国产精品av视频在线免费观看| 国产老妇伦熟女老妇高清| 欧美高清成人免费视频www| 国产一区二区在线观看日韩| 日韩一区二区三区影片| 国产精品一区二区性色av| 97超视频在线观看视频| 亚洲精品456在线播放app| av.在线天堂| 69av精品久久久久久| 国产在线男女| 亚洲国产精品专区欧美| 亚洲一区高清亚洲精品| 国产精品乱码一区二三区的特点| 九九久久精品国产亚洲av麻豆| 人妻系列 视频| 成人高潮视频无遮挡免费网站| 亚洲一级一片aⅴ在线观看| 春色校园在线视频观看| 又黄又爽又刺激的免费视频.| 国产一级毛片七仙女欲春2| 免费看美女性在线毛片视频| 特大巨黑吊av在线直播| 国产午夜精品论理片| 哪个播放器可以免费观看大片| 亚洲国产精品成人综合色| 亚洲国产精品国产精品| 夜夜看夜夜爽夜夜摸| 国产精品av视频在线免费观看| 亚洲欧美日韩卡通动漫| 小蜜桃在线观看免费完整版高清| 日本五十路高清| 午夜精品在线福利| 亚洲,欧美,日韩| 国产精品一区二区三区四区久久| 国产精品蜜桃在线观看| 老司机影院毛片| 偷拍熟女少妇极品色| 欧美一级a爱片免费观看看| 亚洲国产精品sss在线观看| 国产精华一区二区三区| 中文字幕av在线有码专区| www日本黄色视频网| 日日摸夜夜添夜夜爱| 精品久久久噜噜| 亚洲中文字幕日韩| 亚洲自偷自拍三级| 亚洲精品国产av成人精品| 亚洲三级黄色毛片| 汤姆久久久久久久影院中文字幕 | av女优亚洲男人天堂| kizo精华| 亚洲伊人久久精品综合 | 99久久成人亚洲精品观看| 成年女人看的毛片在线观看| 97超视频在线观看视频| 日韩 亚洲 欧美在线| 欧美人与善性xxx| 少妇丰满av| 熟女人妻精品中文字幕| 不卡视频在线观看欧美| 中文精品一卡2卡3卡4更新| 国产精品久久久久久久电影| 日韩制服骚丝袜av| 欧美色视频一区免费| 国产午夜精品论理片| 搡女人真爽免费视频火全软件| 日韩亚洲欧美综合| 日韩成人伦理影院| 最近最新中文字幕免费大全7| a级毛色黄片| 精品不卡国产一区二区三区| 91久久精品国产一区二区三区| 直男gayav资源| 亚洲精品久久久久久婷婷小说 | 亚洲不卡免费看| 91久久精品电影网| 国产高清国产精品国产三级 | 亚洲国产精品国产精品| 久久精品人妻少妇| 99国产精品一区二区蜜桃av| 久久精品熟女亚洲av麻豆精品 | 男女那种视频在线观看| 人人妻人人澡人人爽人人夜夜 | 日日啪夜夜撸| 国产高清有码在线观看视频| 久久久a久久爽久久v久久| 在线观看av片永久免费下载| 日韩成人av中文字幕在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久热精品热| 日本爱情动作片www.在线观看| 精品人妻一区二区三区麻豆| 22中文网久久字幕| 久久精品久久精品一区二区三区| av在线老鸭窝| 精品99又大又爽又粗少妇毛片| 91精品伊人久久大香线蕉| 免费人成在线观看视频色| 五月伊人婷婷丁香| 国产老妇伦熟女老妇高清| 国产成人freesex在线| 亚洲18禁久久av| 日本黄色片子视频| 亚洲av二区三区四区| 九九热线精品视视频播放| av在线老鸭窝| 亚洲av男天堂| 中文字幕免费在线视频6| 亚洲欧美精品专区久久| 午夜激情福利司机影院| 精品国内亚洲2022精品成人| 久久久久久久久大av| 欧美极品一区二区三区四区| 久久久午夜欧美精品| 欧美成人一区二区免费高清观看| 秋霞伦理黄片| 五月伊人婷婷丁香| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 日本爱情动作片www.在线观看| 99热这里只有是精品50| 免费观看a级毛片全部| 99久久精品国产国产毛片| 特级一级黄色大片| 国产国拍精品亚洲av在线观看| 国产单亲对白刺激| 三级国产精品欧美在线观看| 99久久无色码亚洲精品果冻| 国产伦理片在线播放av一区| 少妇的逼好多水| 国产色爽女视频免费观看| 精品国产三级普通话版| 成年免费大片在线观看| 纵有疾风起免费观看全集完整版 | 黄色日韩在线| 欧美性猛交╳xxx乱大交人| 国产成人a区在线观看| av免费观看日本| 能在线免费观看的黄片| 免费av不卡在线播放| 观看美女的网站| 国产伦精品一区二区三区四那| 亚洲美女视频黄频| 国产精品久久久久久久久免| 成人毛片60女人毛片免费| 精品久久久久久久末码| 日韩在线高清观看一区二区三区| 国产私拍福利视频在线观看| 老司机福利观看| 2022亚洲国产成人精品| 日本一二三区视频观看| 黄色配什么色好看| 国产精品精品国产色婷婷| 午夜精品一区二区三区免费看| 精品国产三级普通话版| 麻豆av噜噜一区二区三区| 麻豆久久精品国产亚洲av| 精品久久久久久久末码| 日韩国内少妇激情av| 成人漫画全彩无遮挡| 久久久久免费精品人妻一区二区| 国产一区亚洲一区在线观看| 国产在线一区二区三区精 | 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 嫩草影院入口| 在线天堂最新版资源| 日韩亚洲欧美综合| 国产高清三级在线| 成人二区视频| 美女国产视频在线观看| 亚洲丝袜综合中文字幕| 视频中文字幕在线观看| 99久久九九国产精品国产免费| 欧美激情在线99| 久久久久久久久久久免费av| 日本一二三区视频观看| 日韩av在线免费看完整版不卡| 久久精品久久久久久久性| 久久精品国产自在天天线| 国产综合懂色| 亚洲av.av天堂| 国产一区有黄有色的免费视频 | 熟妇人妻久久中文字幕3abv| 国产高清有码在线观看视频| 亚洲色图av天堂| 简卡轻食公司| 搞女人的毛片| 中文字幕熟女人妻在线| 99热精品在线国产| 亚洲熟妇中文字幕五十中出| 啦啦啦观看免费观看视频高清| 国产精品久久久久久精品电影| 欧美极品一区二区三区四区| 精品人妻视频免费看| 免费黄色在线免费观看| 国产爱豆传媒在线观看| 亚洲四区av| 天堂中文最新版在线下载 | 最近2019中文字幕mv第一页| 日本黄大片高清| 毛片女人毛片| 国产一级毛片七仙女欲春2| 美女cb高潮喷水在线观看| 插逼视频在线观看| 亚洲高清免费不卡视频| 国产黄片美女视频| 91精品伊人久久大香线蕉| 色综合色国产| 亚洲aⅴ乱码一区二区在线播放| 国产精品乱码一区二三区的特点| 亚洲av一区综合| 国产伦在线观看视频一区| 偷拍熟女少妇极品色| av在线亚洲专区| 偷拍熟女少妇极品色| 亚洲性久久影院| 青青草视频在线视频观看| 久久久久久久久久成人| 亚洲欧美日韩卡通动漫| 伦理电影大哥的女人| 最近中文字幕高清免费大全6| 一个人看视频在线观看www免费| 亚洲av一区综合| 精品少妇黑人巨大在线播放 | 国产成人a区在线观看| 亚洲18禁久久av| 精品久久久噜噜| 久久久精品94久久精品| 精品久久久久久久人妻蜜臀av| 男女国产视频网站| 国产伦一二天堂av在线观看| 亚洲国产精品成人综合色| 亚洲av二区三区四区| 国产精华一区二区三区| 男女啪啪激烈高潮av片| 桃色一区二区三区在线观看| 久久国产乱子免费精品| 亚洲国产成人一精品久久久| 亚洲美女视频黄频| 日韩欧美国产在线观看| 国产精品国产三级国产av玫瑰| 99热这里只有是精品50| 变态另类丝袜制服|