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

    基于ERT技術(shù)的含水合物沉積物可視化探測(cè)模擬實(shí)驗(yàn)

    2021-12-30 03:29:10趙金環(huán)劉昌嶺鄒長(zhǎng)春陳強(qiáng)孟慶國(guó)劉洋卜慶濤
    關(guān)鍵詞:層析成像塊狀反應(yīng)釜

    趙金環(huán),劉昌嶺,鄒長(zhǎng)春,陳強(qiáng),孟慶國(guó),劉洋,卜慶濤

    1. 自然資源部天然氣水合物重點(diǎn)實(shí)驗(yàn)室,中國(guó)地質(zhì)調(diào)查局青島海洋地質(zhì)研究所,青島 266237

    2. 青島海洋科學(xué)與技術(shù)試點(diǎn)國(guó)家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室,青島 266237

    3. 中國(guó)地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院,北京 100083

    天然氣水合物被視為21世紀(jì)可供開發(fā)的新能源。天然氣水合物儲(chǔ)層識(shí)別和飽和度估算是水合物資源勘探開發(fā)的重要任務(wù)之一[1-2]。近幾十年來(lái),世界范圍內(nèi)開展了大量的天然氣水合物的實(shí)驗(yàn)?zāi)M研究和現(xiàn)場(chǎng)勘探工作,致力于通過(guò)了解含天然氣水合物沉積物的物性參數(shù)演化識(shí)別水合物并估算飽和度[3-12]。基于電學(xué)響應(yīng)特征的儲(chǔ)層定性與定量評(píng)價(jià)在天然氣水合物的勘探開發(fā)中發(fā)揮著重要作用,是天然氣水合物資源評(píng)價(jià)的重要手段之一。

    隨著越來(lái)越多的原位水合物樣品被獲取,天然氣水合物在沉積物內(nèi)的賦存形態(tài)已基本明確,分為孔隙型(分散狀)和裂隙型(塊狀、脈狀和層狀等)。天然氣水合物賦存形態(tài)多樣性使得儲(chǔ)層各向異性特征顯著,現(xiàn)場(chǎng)勘探獲取的水合物飽和度相近時(shí),電阻率響應(yīng)差異十分明顯[13-16],極大地影響現(xiàn)場(chǎng)基于測(cè)井資料評(píng)估儲(chǔ)層中水合物飽和度的準(zhǔn)確性。物理模擬實(shí)驗(yàn)結(jié)果同樣表明,即使在相同平均水合物飽和度條件下,采用電測(cè)法獲得的沉積物體系電阻率也不盡相同[17-19]。目前,水合物物理模擬實(shí)驗(yàn)大部分采用的是傳統(tǒng)的電阻率點(diǎn)測(cè)方法,獲取的是測(cè)量電極周圍沉積物的綜合電學(xué)響應(yīng)。然而,含水合物沉積物的電阻率響應(yīng)不僅受沉積物固有性質(zhì)、孔隙流體和孔隙尺度上水合物的分布模式的影響,而且還受水合物宏觀賦存形態(tài)的控制[20]。傳統(tǒng)的電阻率點(diǎn)測(cè)不能反映水合物的賦存形態(tài),一定程度上限制了水合物的有效識(shí)別和飽和度的準(zhǔn)確估算。

    電阻率層析成像技術(shù)(Eleactrical Resistivity Tomography,ERT)通過(guò)獲取被測(cè)物體周圍的電壓信號(hào),建立被測(cè)物體的電阻率分布圖像,可用于直觀地呈現(xiàn)水合物的分布情況。成像算法是電阻率層析成像技術(shù)的核心。對(duì)于非均質(zhì)性和各向異性較強(qiáng)的水合物儲(chǔ)層,成像算法是研究的重點(diǎn),以期獲得高分辨率和高準(zhǔn)確度的水合物電阻率分布圖像。目前,電阻率層析成像技術(shù)已在孔隙型水合物樣品上進(jìn)行應(yīng)用[21-26],而裂隙型水合物的電阻率成像還有待進(jìn)一步的研究。

    本文利用水合物電阻率層析成像實(shí)驗(yàn)裝置,模擬分散狀和塊狀水合物的生成過(guò)程;基于物理模擬實(shí)驗(yàn)測(cè)試結(jié)果,建立了含水合物沉積物電阻率層析成像反演方法,實(shí)現(xiàn)分散狀和塊狀水合物的可視化;基于含不同賦存類型水合物沉積物的電阻率圖像,揭示水合物生長(zhǎng)規(guī)律,研究水合物賦存形態(tài)對(duì)電阻率響應(yīng)特征的影響,以期為天然氣水合物賦存類型識(shí)別和飽和度估算提供支持。

    1 實(shí)驗(yàn)

    本研究采用的是青島海洋地質(zhì)研究所自主研制的水合物電阻率層析成像實(shí)驗(yàn)裝置。實(shí)驗(yàn)裝置可實(shí)時(shí)采集原位水合物生成和分解過(guò)程中電阻率層析成像數(shù)據(jù)。根據(jù)電阻率層析成像算法,重建不同時(shí)刻沉積物內(nèi)水合物賦存形態(tài)電阻率圖像。

    1.1 實(shí)驗(yàn)裝置

    水合物電阻率層析成像實(shí)驗(yàn)裝置包括水合物生成/分解和數(shù)據(jù)采集兩個(gè)模塊(圖1)。水合物生成/分解模塊主要包括高壓反應(yīng)釜、恒溫水浴槽和內(nèi)筒,是水合物生成和分解的場(chǎng)所。高壓反應(yīng)釜和恒溫水浴槽用于模擬水合物生成的高壓低溫環(huán)境。內(nèi)筒為一直徑10 cm的圓柱狀絕緣筒,用于承裝沉積物,筒壁和筒底部均布置透氣不透水通道;數(shù)據(jù)采集模塊主要包括測(cè)量電極、電僑儀和控制器,采集水合物生成/分解過(guò)程中的電學(xué)參數(shù)。測(cè)量電極布設(shè)在內(nèi)筒的同一高度,徑向均布16個(gè)測(cè)點(diǎn)、間距22.5°;電僑儀與控制器的相互配合,完成電學(xué)參數(shù)的采集。水合物生成過(guò)程中電學(xué)數(shù)據(jù)的采集方式為四電極法,同一高度的16個(gè)電極中,選相鄰的兩個(gè)電極作為供電電極,其余電極相鄰的兩兩一對(duì)作為測(cè)量電極,直至16個(gè)電極皆完成作為供電電極任務(wù),每一輪次的采集可獲取208個(gè)有效數(shù)據(jù)。

    圖1 水合物電阻率層析成像實(shí)驗(yàn)裝置Fig.1 Experiment device of hydrate electrical resistivity tomography

    1.2 實(shí)驗(yàn)方法及條件

    (1)分散狀水合物

    將充分飽水(濃度為3.5%的NaCl溶液)的石英砂(>500 μm)填入內(nèi)筒后放入高壓反應(yīng)釜,密封反應(yīng)釜;向密封反應(yīng)釜內(nèi)充入甲烷氣體(純度99.99%)至實(shí)驗(yàn)壓力;穩(wěn)定后,開啟水浴循環(huán)進(jìn)行系統(tǒng)降溫,開始水合物生成實(shí)驗(yàn)。實(shí)驗(yàn)過(guò)程中,實(shí)時(shí)采集含分散狀水合物沉積物體系的電阻率層析成像數(shù)據(jù)。

    (2)塊狀水合物

    為提高實(shí)驗(yàn)效率,實(shí)驗(yàn)將已經(jīng)制備好的圓柱狀水合物裝入沉積物內(nèi),以模擬塊狀的裂隙型水合物。參照已有的裂隙充填型水合物合成實(shí)驗(yàn)[27],將去離子水冰粉裝入圓柱狀模具后,放入液氮中冷卻;然后將其放入水合物樣品制備反應(yīng)釜,密閉后充入甲烷氣(純度99.99%)至水合物生成所需的壓力;將反應(yīng)釜放入循環(huán)水浴控溫裝置中進(jìn)行降溫,生成圓柱狀水合物。水合物開始生成后,釜內(nèi)的壓力會(huì)不斷地降低。水合物的生成是一個(gè)從成核到不斷生長(zhǎng)的過(guò)程,十分耗時(shí)。在壓力下降過(guò)程中,也就是冰粉還未完全轉(zhuǎn)換成水合物時(shí),開釜后取出的圓柱狀水合物樣品為冰粉與水合物的混合物(圖2)。

    圖2 圓柱狀水合物樣品(模擬塊狀水合物)Fig.2 Photos of cylindrical hydrate sample which simulates massive hydrate

    為了盡量避免圓柱狀水合物與沉積物接觸導(dǎo)致水合物分解,將充分飽水(濃度為3.5%的NaCl溶液)的石英砂(>500 μm)放入恒溫箱進(jìn)行預(yù)冷,并提前開啟水浴循環(huán)進(jìn)行高壓反應(yīng)釜系統(tǒng)的降溫。待溫度穩(wěn)定后,取出壓力還處在下降過(guò)程中的圓柱狀水合物,放入液氮中冷卻,盡量降低水合物的分解。取出部分預(yù)冷后的石英砂填入內(nèi)筒,將圓柱狀水合物埋入石英砂中,并將內(nèi)筒整體放入反應(yīng)釜內(nèi),快速密封反應(yīng)釜,并迅速向反應(yīng)釜內(nèi)充入甲烷氣體至實(shí)驗(yàn)壓力,盡量避免水合物的分解。由于裝樣過(guò)程和充入的室溫環(huán)境下的甲烷氣體導(dǎo)致沉積物體系溫度升高,高于水合物穩(wěn)定存在的溫度范圍,水合物將存在分解現(xiàn)象。維持水浴循環(huán)的低溫對(duì)沉積物體系進(jìn)行降溫。隨著溫度的降低,水合物將繼續(xù)生成。實(shí)驗(yàn)過(guò)程中,實(shí)時(shí)采集含塊狀水合物沉積物體系的電阻率層析成像數(shù)據(jù)。

    2 層析成像反演算法

    本文電阻率層析成像圖像重建算法采用引入了正則化矩陣和參數(shù)獲取穩(wěn)定解、收斂速度快的高斯牛頓迭代方法求解目標(biāo)函數(shù)。電阻率層析成像屬于非線性的不適定反問(wèn)題,反問(wèn)題等價(jià)于非線性最優(yōu)化問(wèn)題,其最小化目標(biāo)函數(shù)為:

    其中:G為加權(quán)靈敏度矩陣;X為模型參數(shù);d為測(cè)量數(shù)據(jù);λ為正則化參數(shù);C為光滑度矩陣。

    第i+1次迭代得到的新模型為第i次模型X與修正量?X之和,則目標(biāo)函數(shù)可寫為:

    對(duì)目標(biāo)函數(shù)進(jìn)行泰勒展開后,忽略高階項(xiàng),并令其偏導(dǎo)為零,則有:

    其中,J為雅可比矩陣。精度時(shí),迭代終止,Xi即所求模型。

    當(dāng)新模型的解與測(cè)量數(shù)據(jù)的均方根誤差滿足

    3 實(shí)驗(yàn)結(jié)果

    水合物物理模擬實(shí)驗(yàn)開展了多組重復(fù)試驗(yàn),驗(yàn)證實(shí)驗(yàn)裝置和層析成像反演算法的可靠性,本文選取其中的一組數(shù)據(jù)進(jìn)行分析。需要說(shuō)明的是,壓力的變化對(duì)電阻率的影響可忽略不計(jì)[26],而溫度的變化對(duì)孔隙水中導(dǎo)電離子的移動(dòng)影響較大,不可忽略。在層析成像反演過(guò)程中,已將常溫條件下的電阻率作為參考值對(duì)測(cè)量數(shù)據(jù)進(jìn)行了溫度校正。因此,分散狀與塊狀水合物生成的起始溫壓條件的不同不影響水合物物賦存形態(tài)差異引起的電學(xué)響應(yīng)差異的研究。

    對(duì)比分散狀水合物和塊狀水合物的電阻率圖像可以發(fā)現(xiàn),分散狀水合物的高值電阻率分布存在明顯的零散性,塊狀水合物的高值電阻率分布則是相對(duì)比較聚集,且隨著水合物的不斷生成,分散狀水合物的零散性和塊狀水合物的聚集性變得更加明顯。利用二維電阻率圖像上這一特性可初步判斷水合物的賦存形態(tài)。

    3.1 分散狀水合物

    實(shí)驗(yàn)過(guò)程中壓力維持穩(wěn)定10 h后,進(jìn)行了一次補(bǔ)壓和升溫,使沉積物體系內(nèi)的氣水重新分布,促進(jìn)水合物的生成(圖3)。反應(yīng)過(guò)程中隨著水合物生成,體系壓力下降明顯,經(jīng)補(bǔ)壓升溫后,水合物繼續(xù)生成。初始狀態(tài)下,各層的平均電阻率在1.55 Ω·m左右,屬正常海底沉積物電阻率分布范圍內(nèi)(圖4)。分散狀水合物生成過(guò)程中內(nèi)筒內(nèi)不同高度沉積物介質(zhì)的電阻率分布情況如圖5所示。隨著時(shí)間的推移,孔隙水不斷固化成水合物,每層不同位置的電阻率和平均電阻率皆逐漸升高,局部電阻率升高明顯,尤其是上部層位。分析認(rèn)為,水合物的聚集特征受氣源的供給和輸導(dǎo)通道共同控制[28]。Chen等研究認(rèn)為沉積物中水合物首先在氣水接觸面生長(zhǎng),天然氣聚集區(qū)域?qū)⑿纬纱罅康乃衔颷29]??紤]實(shí)驗(yàn)裝置內(nèi)筒上部為開口式,上部氣水接觸充分,下部層位依賴內(nèi)筒側(cè)壁的氣孔進(jìn)氣,氣水接觸相對(duì)較弱。從電阻率圖像上可以觀測(cè)到上部水合物最早生成,生成速度相對(duì)較快,且隨著時(shí)間的變化,大量水合物聚集。

    圖3 分散狀和塊狀水合物生成過(guò)程中溫度和壓力的變化Fig.3 Temperature and pressure during formation of dispersed hydrate and massive hydrate

    圖4 分散狀和塊狀水合物生成過(guò)程每層平均電阻率的變化Fig.4 Average resistivity of each layer during formation of dispersed hydrate and massive hydrate

    從分散狀水合物生成過(guò)程中電阻率圖像的變化可以看出,水合物圍繞水合物優(yōu)先生成區(qū)域不斷擴(kuò)張生長(zhǎng);層與層之間的電阻率和層內(nèi)的電阻率不同指示了水合物的生長(zhǎng)具有明顯的非均質(zhì)性。國(guó)內(nèi)外學(xué)者利用電阻率層析成像技術(shù)在研究水合物生成/分解規(guī)律以及物性演化規(guī)律時(shí)也曾發(fā)現(xiàn)了水合物的上述生長(zhǎng)行為[22-23]。

    隨著時(shí)間的變化,分散狀水合物部分區(qū)域(圖5中橢圓形區(qū)域)電阻率明顯降低,尤其是t3時(shí)刻。分析認(rèn)為是由于水合物生成過(guò)程中的排鹽效應(yīng)導(dǎo)致孔隙流體電阻率降低[30]。水合物生成過(guò)程中孔隙流體被水合物圈閉,被圈閉的這部分流體對(duì)含水合物沉積物的導(dǎo)電貢獻(xiàn)相當(dāng)于水合物[31]。也就是說(shuō),當(dāng)?shù)V化度升高的孔隙流體被水合物圈閉后,同樣會(huì)呈現(xiàn)出高電阻率值。圖5中的長(zhǎng)方形區(qū)域,電阻率出現(xiàn)了升高(t1)-降低(t2)升高(t3)的特征。圖5中t3時(shí)刻方形區(qū)域周圍并沒(méi)有形成高電阻率值的圈閉現(xiàn)象,認(rèn)為圖5中t3時(shí)刻長(zhǎng)方形區(qū)域電阻率的升高是水合物的影響。因此,不同時(shí)刻電阻率圖像的對(duì)比可以明確受水合物影響引起的電阻率的變化,對(duì)水合物的有效識(shí)別十分有利。

    圖5 分散狀水合物生成過(guò)程中不同時(shí)刻電阻率圖像Fig.5 The changes of resistivity distribution of dispersed hydrate with time

    3.2 塊狀水合物

    受實(shí)驗(yàn)裝置硬件條件的影響,塊狀水合物僅成功獲取2層測(cè)量數(shù)據(jù)。塊狀水合物生成過(guò)程中內(nèi)筒內(nèi)沉積物介質(zhì)不同高度的電阻率分布情況如圖6所示。加壓穩(wěn)定后,塊狀水合物的大小呈現(xiàn)準(zhǔn)確,位置稍微有些偏差。實(shí)驗(yàn)發(fā)現(xiàn),隨著時(shí)間的變化圓柱體水合物電阻率不斷增加。非圓柱體水合物區(qū)域處的電阻率值在溫壓條件達(dá)到后的1.5 h(t1時(shí)刻),電阻率值即達(dá)到20 Ω·m,而在分散狀水合物生成過(guò)程歷經(jīng)幾十個(gè)小時(shí)后,電阻率僅有5~7 Ω·m(圖3,圖4)。Lei和Seol通過(guò)對(duì)熱刺激法和降壓法開采的實(shí)驗(yàn)?zāi)M,認(rèn)為外界壓力會(huì)迫使水合物分解產(chǎn)生的氣體滯留在沉積物孔隙內(nèi)[32]。圓柱體水合物放入釜內(nèi)后,受溫度影響水合物仍在分解,分解產(chǎn)生的游離氣被滯留在沉積物內(nèi),驅(qū)替了部分孔隙水。游離氣為不導(dǎo)電介質(zhì),造成其他區(qū)域的電阻率值過(guò)大,進(jìn)而導(dǎo)致圓柱體水合物的位置發(fā)生偏移。同時(shí),孔隙內(nèi)的游離氣與孔隙水接觸充分,促進(jìn)了水合物再次生成,導(dǎo)致非圓柱體水合物區(qū)域處的電阻率值升高較快。

    圖6 塊狀水合物生成過(guò)程中不同時(shí)刻電阻率圖像Fig.6 The changes of resistivity distribution of massive hydrate with time

    4 討論

    在相同沉積物樣品體積條件下,分散狀水合物在壓力下降0.5 MP的情況下,平均電阻率升高了1.22 Ω·m;塊狀水合物在壓力下降約0.12 MP的情況下,平均電阻率升高了14.37 Ω·m(圖3,圖4)。對(duì)比分散狀水合物和塊狀水合物,塊狀水合物電阻率升高較快,而從壓力下降情況可知塊狀水合物生成量又較少。Liu等人通過(guò)數(shù)值模擬研究了水合物形態(tài)對(duì)含天然氣水合物沉積物電阻率特性的影響,發(fā)現(xiàn)低水合物飽和度范圍內(nèi)裂隙型水合物電阻率同樣存在急劇增加的現(xiàn)象[33]。

    電阻率測(cè)量是經(jīng)數(shù)值解法求出測(cè)量電極周圍介質(zhì)的電場(chǎng)分布,然后求取電阻率的計(jì)算公式。測(cè)量電極周圍存在不同賦存形態(tài)的水合物時(shí),電場(chǎng)的分布不同,求取的電阻率也將不同。分散狀水合物是指水合物以顆粒形式分布于粒間孔隙內(nèi),在含量較少的情況下,對(duì)導(dǎo)電路徑(孔隙水)的堵塞作用不明顯;當(dāng)大量的水合物生成后,沉積物內(nèi)的導(dǎo)電路徑將被逐漸堵塞。塊狀水合物模擬實(shí)驗(yàn)的沉積物中,圓柱體水合物為冰粉、孔隙水、水合物和氣的混合。圓柱體水合物是用去離子水制成的冰粉生成的,壓實(shí)后的孔隙空間也很小,可被視為體積相對(duì)較大的不導(dǎo)電介質(zhì),即使水合物生成量相對(duì)較低,對(duì)導(dǎo)電路徑的阻礙作用也會(huì)比較明顯,對(duì)獲取的電學(xué)響應(yīng)的貢獻(xiàn)突出,電阻率升高的就會(huì)較快。

    不同賦存形態(tài)的水合物飽和度估算需要采用不同的電阻率巖石物理模型[12,33-34]。對(duì)于取心難度較大的水合物儲(chǔ)層,水合物的賦存形態(tài)很難獲得。電阻率層析成像技術(shù)在孔隙尺度的微觀探測(cè)與現(xiàn)場(chǎng)尺度的宏觀探測(cè)中具有廣闊的前景[25]。劉洋等人已研發(fā)了一種新的電阻率層析成像陣列,應(yīng)用于含水合物沉積物的電阻率成像模擬實(shí)驗(yàn)裝置[26]。結(jié)合本文研究,電阻率層析成像技術(shù)可作為識(shí)別水合物賦存形態(tài)的一種手段應(yīng)用于現(xiàn)場(chǎng)水合物儲(chǔ)層,為選取電阻率巖石物理模型提供依據(jù),提高基于現(xiàn)場(chǎng)測(cè)井資料評(píng)估儲(chǔ)層中水合物飽和度的準(zhǔn)確性。

    5 結(jié)論

    (1)利用電阻率層析成像結(jié)果可初步判斷沉積物中水合物的賦存形態(tài)。分散狀水合物的高值電阻率分布具有零散性,塊狀水合物的高值電阻率分布比較聚集,且隨著水合物的不斷生成,電阻率高值的零散性和聚集性將越加明顯。

    (2)分散狀水合物不同時(shí)刻電阻率圖像的對(duì)比,可以明確由水合物生成所引起的電阻率的變化。被滯留在沉積物內(nèi)的游離氣對(duì)沉積物的電阻率影響較大,影響了塊狀水合物電阻率層析成像的位置,多物理參數(shù)的聯(lián)合探測(cè)將有助于區(qū)分游離氣和水合物對(duì)沉積物電阻率響應(yīng)的影響。

    (3)水合物的賦存形態(tài)是影響電阻率響應(yīng)特征的重要因素,分散狀水合物和塊狀水合物的電阻率響應(yīng)特征差異明顯;相對(duì)于分散狀水合物,塊狀水合物作為大體積的不導(dǎo)電介質(zhì),對(duì)導(dǎo)電路徑的堵塞作用明顯,導(dǎo)致塊狀水合物電阻率隨著水合物的生成升高較快。

    猜你喜歡
    層析成像塊狀反應(yīng)釜
    基于大數(shù)據(jù)量的初至層析成像算法優(yōu)化
    貴溪冶煉廠臥式反應(yīng)釜自動(dòng)控制方法的研究
    基于快速行進(jìn)法地震層析成像研究
    改進(jìn)PID在反應(yīng)釜溫度控制系統(tǒng)中的應(yīng)用研究
    對(duì)于反應(yīng)釜機(jī)械密封失效與改造的分析
    Ghosts in the shell: identif i cation of microglia in the human central nervous system by P2Y12 receptor
    厚層塊狀特低滲礫巖油藏水平井壓裂參數(shù)優(yōu)化
    某反應(yīng)釜進(jìn)料加熱器管板的應(yīng)力和疲勞分析
    LTE上行塊狀導(dǎo)頻的信道估計(jì)研究
    基于分布式無(wú)線網(wǎng)絡(luò)的無(wú)線電層析成像方法與實(shí)驗(yàn)研究
    亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av第一区精品v没综合| 国产精品98久久久久久宅男小说| 免费观看人在逋| 欧美精品啪啪一区二区三区| 99久久国产精品久久久| 制服人妻中文乱码| 咕卡用的链子| av中文乱码字幕在线| 国产精品二区激情视频| 欧美日韩成人在线一区二区| 国内久久婷婷六月综合欲色啪| 在线国产一区二区在线| 久久香蕉精品热| 国产精品久久久av美女十八| 国产有黄有色有爽视频| 亚洲成a人片在线一区二区| 亚洲情色 制服丝袜| 一边摸一边抽搐一进一出视频| 热re99久久精品国产66热6| 国产在线观看jvid| 色94色欧美一区二区| 国产黄色免费在线视频| 正在播放国产对白刺激| 一级a爱视频在线免费观看| 亚洲精品在线观看二区| 亚洲熟妇熟女久久| 一区在线观看完整版| 欧美国产精品一级二级三级| av免费在线观看网站| tube8黄色片| 亚洲精品国产精品久久久不卡| 麻豆av在线久日| 99热网站在线观看| 女人爽到高潮嗷嗷叫在线视频| 丝袜在线中文字幕| 国产一区二区三区综合在线观看| 视频在线观看一区二区三区| 国产精品免费大片| 亚洲精品久久午夜乱码| 美女视频免费永久观看网站| 国产精品二区激情视频| 亚洲人成电影观看| 日韩欧美在线二视频 | 激情在线观看视频在线高清 | 男人的好看免费观看在线视频 | 丰满人妻熟妇乱又伦精品不卡| 国产成+人综合+亚洲专区| 999久久久国产精品视频| 天堂俺去俺来也www色官网| 在线看a的网站| 一区二区三区激情视频| 黑人操中国人逼视频| 欧美乱码精品一区二区三区| 手机成人av网站| 久久久久久免费高清国产稀缺| 日韩成人在线观看一区二区三区| av在线播放免费不卡| 亚洲av熟女| 精品午夜福利视频在线观看一区| 下体分泌物呈黄色| 99热只有精品国产| 欧美乱妇无乱码| 精品久久久久久久久久免费视频 | 色老头精品视频在线观看| 久久狼人影院| 一区二区三区精品91| 久久热在线av| 亚洲 国产 在线| 视频区欧美日本亚洲| 热99re8久久精品国产| 99国产精品99久久久久| 午夜福利视频在线观看免费| 国产av精品麻豆| 免费不卡黄色视频| 欧美精品啪啪一区二区三区| 国产精品国产av在线观看| 1024视频免费在线观看| 超碰成人久久| 成人精品一区二区免费| 中文字幕人妻丝袜制服| 欧美色视频一区免费| 一级片'在线观看视频| 亚洲欧美日韩高清在线视频| 国产成人精品久久二区二区免费| 老司机亚洲免费影院| 变态另类成人亚洲欧美熟女 | 国产精品一区二区免费欧美| 王馨瑶露胸无遮挡在线观看| 亚洲国产欧美日韩在线播放| 一区二区日韩欧美中文字幕| 精品国产一区二区三区久久久樱花| 亚洲五月婷婷丁香| 老司机福利观看| 黄色怎么调成土黄色| 久久久久国内视频| 麻豆国产av国片精品| 精品午夜福利视频在线观看一区| 99精品在免费线老司机午夜| 亚洲成人手机| 97人妻天天添夜夜摸| 亚洲精品在线美女| 精品电影一区二区在线| 免费少妇av软件| 91字幕亚洲| 国产精品一区二区在线观看99| 亚洲成a人片在线一区二区| av电影中文网址| 无遮挡黄片免费观看| 久久人妻熟女aⅴ| 久久精品国产99精品国产亚洲性色 | 激情视频va一区二区三区| 老汉色av国产亚洲站长工具| 欧美黄色淫秽网站| 亚洲精品一卡2卡三卡4卡5卡| 在线免费观看的www视频| 两性夫妻黄色片| 丁香六月欧美| 免费看a级黄色片| 性色av乱码一区二区三区2| 欧美成人免费av一区二区三区 | 久久人妻福利社区极品人妻图片| 国产精品免费大片| 三上悠亚av全集在线观看| 成在线人永久免费视频| av一本久久久久| 999精品在线视频| 人人妻人人澡人人爽人人夜夜| 最新在线观看一区二区三区| 女同久久另类99精品国产91| 欧美日本中文国产一区发布| 老司机午夜十八禁免费视频| 热99国产精品久久久久久7| 一区福利在线观看| 日本一区二区免费在线视频| 亚洲色图av天堂| 欧美激情 高清一区二区三区| 操出白浆在线播放| 下体分泌物呈黄色| 亚洲av日韩精品久久久久久密| 国产97色在线日韩免费| 不卡av一区二区三区| 80岁老熟妇乱子伦牲交| 国产国语露脸激情在线看| 色播在线永久视频| 久久精品国产亚洲av香蕉五月 | 精品一区二区三区视频在线观看免费 | 黄网站色视频无遮挡免费观看| 天堂√8在线中文| 午夜福利影视在线免费观看| 亚洲视频免费观看视频| 很黄的视频免费| 欧美不卡视频在线免费观看 | 激情在线观看视频在线高清 | 美国免费a级毛片| 村上凉子中文字幕在线| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 国产一区有黄有色的免费视频| 电影成人av| 人妻丰满熟妇av一区二区三区 | 操出白浆在线播放| 一二三四在线观看免费中文在| 免费黄频网站在线观看国产| 老熟妇仑乱视频hdxx| 中国美女看黄片| 在线av久久热| www.精华液| 欧美国产精品va在线观看不卡| 757午夜福利合集在线观看| 久久久久久人人人人人| 99热只有精品国产| 免费看十八禁软件| 日韩一卡2卡3卡4卡2021年| 精品国内亚洲2022精品成人 | 国产伦人伦偷精品视频| 亚洲欧美一区二区三区久久| 国产精品98久久久久久宅男小说| 午夜免费成人在线视频| 免费看a级黄色片| 国产日韩欧美亚洲二区| 大香蕉久久网| 亚洲精品自拍成人| 大码成人一级视频| 咕卡用的链子| 丁香欧美五月| 中亚洲国语对白在线视频| 黄色丝袜av网址大全| 国产免费男女视频| 男女床上黄色一级片免费看| 操美女的视频在线观看| 妹子高潮喷水视频| 亚洲av熟女| 免费观看精品视频网站| 国产一区有黄有色的免费视频| 99riav亚洲国产免费| 久久精品国产综合久久久| 久久亚洲精品不卡| 夜夜夜夜夜久久久久| 亚洲欧美精品综合一区二区三区| 国产精品九九99| 成熟少妇高潮喷水视频| 乱人伦中国视频| 久久天躁狠狠躁夜夜2o2o| 国产激情欧美一区二区| 久久天堂一区二区三区四区| 久久这里只有精品19| 欧美精品亚洲一区二区| 91av网站免费观看| xxx96com| 国产主播在线观看一区二区| 午夜老司机福利片| 日本黄色视频三级网站网址 | 国产1区2区3区精品| 国产又爽黄色视频| 日韩有码中文字幕| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人爽人人夜夜| 男女高潮啪啪啪动态图| 精品亚洲成a人片在线观看| 色精品久久人妻99蜜桃| 韩国av一区二区三区四区| 久久天堂一区二区三区四区| 久久这里只有精品19| 国产精品一区二区精品视频观看| 亚洲九九香蕉| 亚洲 欧美一区二区三区| 日韩 欧美 亚洲 中文字幕| 男人操女人黄网站| 亚洲欧美日韩另类电影网站| 99国产综合亚洲精品| 国产精品一区二区免费欧美| 99re在线观看精品视频| 午夜老司机福利片| 999精品在线视频| 亚洲aⅴ乱码一区二区在线播放 | 激情视频va一区二区三区| 亚洲片人在线观看| 国产亚洲欧美精品永久| 制服诱惑二区| 国产男女内射视频| 看片在线看免费视频| 18禁观看日本| 久久久国产精品麻豆| 国产成人av教育| 男男h啪啪无遮挡| 老熟妇仑乱视频hdxx| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品乱久久久久久| 亚洲av成人不卡在线观看播放网| 在线观看午夜福利视频| 亚洲午夜精品一区,二区,三区| 久久精品人人爽人人爽视色| 亚洲精品久久午夜乱码| 精品乱码久久久久久99久播| 日韩欧美在线二视频 | 高清在线国产一区| 国产淫语在线视频| 久久久国产成人精品二区 | 亚洲精品粉嫩美女一区| 国产在线精品亚洲第一网站| 精品人妻在线不人妻| 很黄的视频免费| 超碰97精品在线观看| 大型av网站在线播放| 亚洲中文日韩欧美视频| 国产三级黄色录像| 黄色a级毛片大全视频| 午夜福利欧美成人| 精品国产一区二区久久| 久久久久久人人人人人| 两性夫妻黄色片| 水蜜桃什么品种好| 久久99一区二区三区| 一区二区三区国产精品乱码| 捣出白浆h1v1| 9色porny在线观看| 欧美激情 高清一区二区三区| av网站在线播放免费| 日韩制服丝袜自拍偷拍| 男女之事视频高清在线观看| 91精品国产国语对白视频| 精品国产美女av久久久久小说| 国产精品亚洲av一区麻豆| 国产精品综合久久久久久久免费 | 18禁国产床啪视频网站| 熟女少妇亚洲综合色aaa.| 亚洲 国产 在线| 亚洲五月婷婷丁香| 黄网站色视频无遮挡免费观看| 99国产精品一区二区蜜桃av | 性色av乱码一区二区三区2| 国产人伦9x9x在线观看| 曰老女人黄片| 夜夜躁狠狠躁天天躁| 亚洲中文日韩欧美视频| 亚洲av成人一区二区三| 18禁国产床啪视频网站| 老司机午夜福利在线观看视频| 美女视频免费永久观看网站| 制服诱惑二区| 建设人人有责人人尽责人人享有的| 亚洲精品中文字幕一二三四区| 中文字幕av电影在线播放| 午夜福利,免费看| 法律面前人人平等表现在哪些方面| 午夜福利视频在线观看免费| 亚洲一区高清亚洲精品| 国产欧美日韩一区二区三区在线| 黄片大片在线免费观看| 久久香蕉精品热| 一级片免费观看大全| 欧美 亚洲 国产 日韩一| videosex国产| 最近最新中文字幕大全电影3 | 国产黄色免费在线视频| 黑人巨大精品欧美一区二区mp4| 国产亚洲精品久久久久5区| 真人做人爱边吃奶动态| 美女高潮到喷水免费观看| 精品久久久久久久久久免费视频 | 捣出白浆h1v1| 我的亚洲天堂| videos熟女内射| www.999成人在线观看| 动漫黄色视频在线观看| av视频免费观看在线观看| 亚洲欧美一区二区三区黑人| 久久久国产成人精品二区 | 熟女少妇亚洲综合色aaa.| 精品国产乱子伦一区二区三区| 日韩中文字幕欧美一区二区| 午夜福利,免费看| 操出白浆在线播放| 欧美国产精品va在线观看不卡| 午夜福利乱码中文字幕| 在线观看免费日韩欧美大片| 免费观看精品视频网站| 热99久久久久精品小说推荐| 色尼玛亚洲综合影院| 久久久国产欧美日韩av| 精品无人区乱码1区二区| 亚洲全国av大片| 国产av又大| 日韩中文字幕欧美一区二区| 如日韩欧美国产精品一区二区三区| 一级a爱视频在线免费观看| 亚洲av日韩在线播放| 一本大道久久a久久精品| 久久国产精品大桥未久av| 看片在线看免费视频| 国产一区二区三区综合在线观看| 免费观看精品视频网站| 美女高潮到喷水免费观看| 精品少妇久久久久久888优播| 99热只有精品国产| 人人澡人人妻人| 91麻豆精品激情在线观看国产 | 亚洲欧美精品综合一区二区三区| 最近最新免费中文字幕在线| 中文亚洲av片在线观看爽 | 精品电影一区二区在线| 18禁黄网站禁片午夜丰满| 美女视频免费永久观看网站| 国产精品美女特级片免费视频播放器 | svipshipincom国产片| 国产又爽黄色视频| 亚洲avbb在线观看| 一进一出好大好爽视频| 丰满迷人的少妇在线观看| a级片在线免费高清观看视频| 午夜成年电影在线免费观看| 在线永久观看黄色视频| 国产高清国产精品国产三级| 国产一卡二卡三卡精品| 黑人欧美特级aaaaaa片| www.999成人在线观看| 亚洲精品久久午夜乱码| 岛国毛片在线播放| 亚洲精品av麻豆狂野| 国产欧美亚洲国产| 久久国产亚洲av麻豆专区| videosex国产| 搡老熟女国产l中国老女人| 色94色欧美一区二区| 亚洲国产欧美网| 露出奶头的视频| 王馨瑶露胸无遮挡在线观看| 丝袜美腿诱惑在线| 亚洲一区中文字幕在线| 久久久久久久精品吃奶| 水蜜桃什么品种好| 国产精品成人在线| 国产精品久久久久久精品古装| 制服诱惑二区| 999久久久国产精品视频| 日本wwww免费看| 国产在线一区二区三区精| 夜夜夜夜夜久久久久| 国产野战对白在线观看| 色综合欧美亚洲国产小说| 欧美av亚洲av综合av国产av| 日韩欧美在线二视频 | 欧美国产精品一级二级三级| 在线观看66精品国产| 色播在线永久视频| 久热这里只有精品99| 在线观看日韩欧美| 在线观看免费视频日本深夜| 欧美日韩乱码在线| 欧美乱码精品一区二区三区| 黄片小视频在线播放| svipshipincom国产片| 久久国产精品男人的天堂亚洲| 制服诱惑二区| 老司机深夜福利视频在线观看| 热99久久久久精品小说推荐| 黄色片一级片一级黄色片| 免费日韩欧美在线观看| 欧美日韩中文字幕国产精品一区二区三区 | svipshipincom国产片| 伊人久久大香线蕉亚洲五| 性少妇av在线| 国产成人系列免费观看| 黄色丝袜av网址大全| 国产一卡二卡三卡精品| 亚洲精品av麻豆狂野| 成人亚洲精品一区在线观看| 一本大道久久a久久精品| cao死你这个sao货| 侵犯人妻中文字幕一二三四区| 18禁裸乳无遮挡免费网站照片 | 欧美最黄视频在线播放免费 | 国产亚洲精品久久久久5区| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久久久久久大奶| 亚洲成人免费电影在线观看| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 一进一出抽搐动态| 亚洲欧美色中文字幕在线| 成熟少妇高潮喷水视频| 国产成人系列免费观看| 久久久久久亚洲精品国产蜜桃av| 久久精品国产a三级三级三级| 国产欧美日韩综合在线一区二区| 大香蕉久久网| 国产精品欧美亚洲77777| 国产在线一区二区三区精| 嫩草影视91久久| 老司机影院毛片| 久久人妻av系列| 国产不卡一卡二| 99热网站在线观看| 亚洲少妇的诱惑av| 操出白浆在线播放| 亚洲专区国产一区二区| 国产三级黄色录像| 黑人欧美特级aaaaaa片| videos熟女内射| 成年人黄色毛片网站| 757午夜福利合集在线观看| 国产精品av久久久久免费| 黄色怎么调成土黄色| 国产主播在线观看一区二区| 性少妇av在线| 亚洲精品在线观看二区| 成人黄色视频免费在线看| www日本在线高清视频| 亚洲中文av在线| 看免费av毛片| 80岁老熟妇乱子伦牲交| 亚洲综合色网址| 激情视频va一区二区三区| 日本a在线网址| 日日夜夜操网爽| 91麻豆精品激情在线观看国产 | 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又粗又硬又大视频| 亚洲熟妇中文字幕五十中出 | 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| 精品久久久久久电影网| 精品国产亚洲在线| 在线看a的网站| 免费在线观看视频国产中文字幕亚洲| 成人亚洲精品一区在线观看| 99国产精品99久久久久| 成人免费观看视频高清| 国产99久久九九免费精品| 欧美午夜高清在线| 99久久精品国产亚洲精品| 纯流量卡能插随身wifi吗| 午夜福利在线观看吧| 少妇猛男粗大的猛烈进出视频| av在线播放免费不卡| 久99久视频精品免费| 麻豆乱淫一区二区| 啪啪无遮挡十八禁网站| 岛国在线观看网站| 日本a在线网址| 久久精品国产综合久久久| 午夜福利,免费看| 女人被狂操c到高潮| 在线观看免费视频日本深夜| 天堂中文最新版在线下载| 99久久综合精品五月天人人| 久久人人97超碰香蕉20202| 久久香蕉激情| 美女午夜性视频免费| 久久热在线av| 精品国产乱码久久久久久男人| 亚洲精品一卡2卡三卡4卡5卡| 午夜精品久久久久久毛片777| 国产在线一区二区三区精| 国产精品免费大片| 美女高潮到喷水免费观看| 国产精品乱码一区二三区的特点 | 狠狠狠狠99中文字幕| 一区二区三区激情视频| 99久久精品国产亚洲精品| 热99久久久久精品小说推荐| 我的亚洲天堂| 亚洲精品粉嫩美女一区| 啦啦啦在线免费观看视频4| 国产高清激情床上av| 欧美另类亚洲清纯唯美| 在线观看一区二区三区激情| 国产成人一区二区三区免费视频网站| 一区二区三区激情视频| 亚洲五月婷婷丁香| 伊人久久大香线蕉亚洲五| 欧美日韩亚洲国产一区二区在线观看 | 最近最新中文字幕大全免费视频| 动漫黄色视频在线观看| 欧美丝袜亚洲另类 | 国产精品久久久av美女十八| 亚洲情色 制服丝袜| 俄罗斯特黄特色一大片| 精品国产亚洲在线| 18禁裸乳无遮挡动漫免费视频| 人妻一区二区av| 好看av亚洲va欧美ⅴa在| 在线观看免费日韩欧美大片| 91成人精品电影| 亚洲情色 制服丝袜| 久久人妻福利社区极品人妻图片| 亚洲熟女毛片儿| 女警被强在线播放| 久久精品国产亚洲av香蕉五月 | 国产成人精品久久二区二区91| 身体一侧抽搐| 熟女少妇亚洲综合色aaa.| 成年人免费黄色播放视频| 下体分泌物呈黄色| 免费在线观看黄色视频的| 在线观看免费视频网站a站| 中文字幕av电影在线播放| 电影成人av| 成人国产一区最新在线观看| 亚洲精品国产一区二区精华液| 99在线人妻在线中文字幕 | 男人操女人黄网站| 国产xxxxx性猛交| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久5区| 999久久久精品免费观看国产| 777米奇影视久久| 一二三四社区在线视频社区8| 欧美人与性动交α欧美软件| 欧美成狂野欧美在线观看| 免费黄频网站在线观看国产| 免费观看精品视频网站| 大型av网站在线播放| 黑人巨大精品欧美一区二区蜜桃| 制服人妻中文乱码| 免费观看人在逋| 欧美日韩亚洲高清精品| av网站在线播放免费| 成人亚洲精品一区在线观看| 亚洲中文字幕日韩| 国产单亲对白刺激| 亚洲精品国产精品久久久不卡| 变态另类成人亚洲欧美熟女 | 午夜精品久久久久久毛片777| 一a级毛片在线观看| 真人做人爱边吃奶动态| 亚洲人成电影观看| 在线视频色国产色| 中文字幕av电影在线播放| 视频在线观看一区二区三区| videos熟女内射| 精品国内亚洲2022精品成人 | 免费日韩欧美在线观看| 一级作爱视频免费观看| 国产精品国产av在线观看| 国产精品久久久久成人av| 国产熟女午夜一区二区三区| 80岁老熟妇乱子伦牲交| 在线观看舔阴道视频| tocl精华| 在线观看免费视频网站a站| 亚洲欧美一区二区三区久久| 波多野结衣一区麻豆| 亚洲片人在线观看| 香蕉久久夜色| 操美女的视频在线观看| 丰满的人妻完整版| 在线观看舔阴道视频| 亚洲色图 男人天堂 中文字幕| 精品视频人人做人人爽| 中文字幕高清在线视频| 亚洲五月婷婷丁香| 久久狼人影院| 亚洲成人手机|