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

    基于WorldView-Ⅱ高分?jǐn)?shù)據(jù)的鹽湖礦化度含量定量反演研究
    ——以尕斯庫勒鹽湖為例

    2021-02-02 09:37:16王俊虎陸冬華周覓武鼎郝偉林
    鈾礦地質(zhì) 2021年1期
    關(guān)鍵詞:礦化度鹽湖反射率

    王俊虎,陸冬華,周覓,武鼎,郝偉林

    (1.核工業(yè)北京地質(zhì)研究院 遙感信息與圖像分析技術(shù)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,北京 100029;2.中國(guó)核工業(yè)集團(tuán)有限公司,北京 100822)

    我國(guó)鹽湖以化學(xué)類型全、分布數(shù)量多、礦產(chǎn)資源豐富和富含稀散元素而聞名于世,且尤以青藏高原、柴達(dá)木盆地、二連盆地等地區(qū)鹽湖數(shù)量最多。我國(guó)鹽湖湖水不僅礦化度高,而且蘊(yùn)含著豐富的鉀、鈉、鎂、硼、鋰、銣、銫及鈾、釷等礦產(chǎn)資源,開發(fā)利用價(jià)值極大[1-2]。據(jù)白朝軍等學(xué)者研究,鹽湖不同水域的礦化度高低不同,只有高礦化度的鹽湖水域才能成礦[3]。張彭熹等通過開展柴達(dá)木盆地鹽湖水礦化度與成礦元素含量的相關(guān)性研究,得出了湖水中鈉、鉀、鎂、硼、鋰等成礦元素含量存在著隨湖水礦化度增高而增加的趨勢(shì)[4]。因此,定量識(shí)別鹽湖水體礦化度對(duì)湖水含礦性評(píng)價(jià)和開發(fā)具有重要的指導(dǎo)意義。傳統(tǒng)獲取鹽湖水礦化度數(shù)據(jù)的方法主要為大量的人工取樣分析并計(jì)算而得,但鹽湖分布區(qū)鹽沼密布、通行困難,自然環(huán)境十分惡劣。因此,傳統(tǒng)的水樣采集和礦化度分析周期長(zhǎng)、難度大、成本高,不利于鹽湖水域礦化度的快速識(shí)別和評(píng)價(jià)。

    遙感技術(shù)具有探測(cè)范圍廣、時(shí)效性強(qiáng)、不受自然條件限制、性價(jià)比高等優(yōu)勢(shì),已在鹽湖成礦環(huán)境識(shí)別和礦產(chǎn)資源調(diào)查評(píng)價(jià)領(lǐng)域發(fā)揮了重要作用[5-8]。如田淑芳等(2005),張大林等(2007)分別基于TM 遙感數(shù)據(jù)定量反演了西藏扎布耶鹽湖總鹽含量和氧化硼含量[6-7],取得了一定研究效果,但限于TM 數(shù)據(jù)空間分辨率偏低,定量反演精度受到一定的限制。王俊虎等(2015)基于SPOT 5 遙感數(shù)據(jù)開展了尕斯庫勒含鈾鹽湖礦化度的定性估測(cè),識(shí)別出了鹽湖水域礦化度變異部位[8],但限于SPOT 5 數(shù)據(jù)波段少,未開展定量反演。WorldView-II 高分辨率遙感衛(wèi)星自2009 年成功發(fā)射以來,因其在可見—近紅外譜帶具有8個(gè)波段且空間分辨率達(dá)1.84 m 的優(yōu)勢(shì),已在水體信息提取等領(lǐng)域開展了廣泛應(yīng)用并取得了良好效果[9-12],但未見應(yīng)用于鹽湖礦化度信息提取領(lǐng)域的相關(guān)研究報(bào)導(dǎo)。

    基于上述背景,筆者基于WorldView-II高分辨率衛(wèi)星數(shù)據(jù),開展了尕斯庫勒鹽湖礦化度含量定量反演研究,評(píng)價(jià)了WorldView-II數(shù)據(jù)的礦化度定量反演能力,旨在為鹽湖湖表鹵水高礦化水域定量識(shí)別提供新的技術(shù)支撐。

    1 研究區(qū)概況和數(shù)據(jù)獲取

    1.1 研究區(qū)概況

    尕斯庫勒鹽湖位于青海省西北部海西州,地理坐標(biāo)為90°40′E~91°10′E,37°57′N~38°11′N,為柴達(dá)木盆地發(fā)展到晚更新世后,湖水變淺后,變?yōu)榇笮〔坏?、散布各地、互不相連的多個(gè)咸水湖和鹽湖之一,類型為硫酸鎂亞型鹽湖。鹽湖湖盆為封閉型內(nèi)流盆地,呈北西向展布,其北東、北西、南西分別由油砂山、阿爾金山和祁漫塔格山環(huán)繞,東南面被東柴山丘陵相隔。阿爾金山和祁漫塔格山主峰海拔均在4 000 m 以上,終年積雪不化,為鹽湖主要的淡水補(bǔ)給來源。鹽湖區(qū)平均海拔高度為2 835 m,為該區(qū)主要匯水洼地。湖表鹵水面積約為103 m2,水深一般位于0.5~0.8 m 之間,鹽湖中南部最深,可達(dá)1.3 m。該湖東部水域淡水補(bǔ)給少,蒸發(fā)量大,平均礦化度約為305 g/L。湖的西部和南部受多條淡水河流補(bǔ)給,礦化度有所降低[13]。

    尕斯庫勒鹽湖是固液相并存的特大型石鹽鹽湖礦床,湖表鹵水和晶間鹵水十分豐富,鹵水中富含鉀、鋰等礦產(chǎn)資源。核工業(yè)北京地質(zhì)研究院鹽湖項(xiàng)目團(tuán)隊(duì)于2010—2012 年,對(duì)該鹽湖開展了系統(tǒng)的含鈾性調(diào)查評(píng)價(jià)。經(jīng)室內(nèi)化學(xué)分析,湖表水平均鈾含量為212 mg/L,屬于潛在的鹽湖型鈾資源[14]。為形象展示尕斯庫勒鹽湖-蝕源區(qū)的地形、地貌及補(bǔ)給-徑流-排泄特征,本文基于中等分辨率ETM+741彩色合成影像,疊加DEM 高程數(shù)據(jù)制作出了鹽湖-蝕源區(qū)補(bǔ)-徑-排三維影像解譯圖(圖1)。

    圖1 基于ETM+741 和DEM 高程數(shù)據(jù)的尕斯庫勒鹽湖補(bǔ)給-徑流-排泄三維解譯圖Fig.1 Interpreted 3D map of the recharge-runoff-discharge zones in Gasikule salt lake based on ETM+741 and DEM

    1.2 數(shù)據(jù)獲取

    1.2.1遙感數(shù)據(jù)獲取

    本次研究獲取了覆蓋尕斯庫勒鹽湖的WorldView-II 多光譜和全色波段數(shù)據(jù),數(shù)據(jù)采集時(shí)間為2011 年8 月23 日,光譜波段范圍為400~1 040 nm,多光譜波段空間分辨率為1.84 m,全色波段空間分辨率為0.46 m(表1)[10]。該數(shù)據(jù)信噪比高,質(zhì)量良好,且均已經(jīng)過輻射定標(biāo)和幾何校正,可滿足本文研究的精度需求。為了提高礦化度定量反演精度,基于ENVI 4.8 軟件平臺(tái)的FLAASH 模塊對(duì)數(shù)據(jù)進(jìn)行了大氣校正,并進(jìn)行歸一化處理,獲得了WorldView-II 數(shù)據(jù)的反射率。

    1.2.2湖水礦化度實(shí)測(cè)數(shù)據(jù)獲取

    為了獲取鹽湖水域湖表鹵水礦化度實(shí)測(cè)值,項(xiàng)目團(tuán)隊(duì)于2010 年劃船進(jìn)入尕斯庫勒湖區(qū)水域,按照預(yù)先設(shè)計(jì)的采樣點(diǎn)位,均勻采集了40 個(gè)湖表水樣,并對(duì)水樣進(jìn)行室內(nèi)陰陽離子分析和礦化度計(jì)算,獲得了40 個(gè)點(diǎn)位湖表水的礦化度值。選取跨越鹽湖水域20 個(gè)水樣的礦化度數(shù)值用于反演建模,剩余20 個(gè)水樣的礦化度數(shù)值用于模型精度評(píng)價(jià)。

    2 鹽湖水域礦化度定量反演

    2.1 鹽湖水際線提取

    鹽湖水際線,即鹽灘和湖水的分界線,可以呈現(xiàn)鹽湖水域分布范圍及其變化狀況。在水體遙感中,因水體對(duì)近紅外波段具有很強(qiáng)的吸收特性,使得陸、水影像特征截然不同。因此,河流、湖泊和海洋的水際線可以用遙感數(shù)據(jù)的近紅外波段來準(zhǔn)確提?。?5]??紤]到WorldView-II 數(shù)據(jù)有兩個(gè)近紅外波段,本文選用了對(duì)水體吸收效果最好的WorldView-II 第8 波段,采用計(jì)算機(jī)自動(dòng)分類結(jié)合人機(jī)對(duì)話解譯,提取出了尕斯庫勒鹽湖的水際線。利用提取出的水際線對(duì)WorldView-II 其余各波段數(shù)據(jù)進(jìn)行掩模處理,獲得了尕斯庫勒鹽湖WorldView-II 各波段的水域分布圖。通過提取鹽湖水際線,消除了湖區(qū)以外信息的干擾,為鹽湖水體礦化度含量定量反演奠定了良好基礎(chǔ)。

    2.2 不同礦化度鹽湖水體光譜信息提取和診斷特征分析

    筆者對(duì)尕斯庫勒鹽湖不同礦化度的鹽湖水進(jìn)行了室內(nèi)光譜測(cè)試和診斷特征分析,得出了鹽湖湖水的光譜反射率與湖水礦化度的高低呈明顯的正相關(guān)性[8]。因此,本文直接從尕斯庫勒鹽湖WorldView-II 彩色影像中依據(jù)反射率值的高低,選取一條反射率值由高到低的剖面,并在剖面上均勻的選取了8 個(gè)點(diǎn)(圖2a),從每個(gè)波段影像中讀取這8 個(gè)點(diǎn)反射率值,形成了WorldView-II 全波段光譜曲線圖(圖2b)。1~8 號(hào)點(diǎn)位湖水的光譜曲線形態(tài)基本相同,在第2 波段 (藍(lán)光波段)處為明顯的吸收特征,在第3 波段(綠光)和第4波段(黃光)處為明顯反射特征,從第5 波段(紅光)~第8 波段(近紅外)影像反射率快速降低。此外,不同礦化度湖水點(diǎn)位在WorldView-II 第1 波段、第3 波段和第4 波段處反射率值最高,在WorldView-II 第7 波段和第8 波段反射率值最低。

    表1 WorldView-II 數(shù)據(jù)波段參數(shù)表Table 1 Band parameters of WorldView-II data

    2.3 鹽湖水體礦化度信息遙感指數(shù)構(gòu)建

    通過研究不同礦化度湖水在WorldView-II多波譜影像的光譜曲線特征得出:在可見光光譜區(qū)間,以WorldView-II 第1 波段、第3波段和第4 波段反射率最高,在近紅外光譜區(qū)間以WorldView-II 第8 波段反射率最低。用這種反射率最高與最低譜段的比值便可以反映不同區(qū)域湖水的礦化度信息。因此,本文建立了基于WorldView-II 數(shù)據(jù)的鹽湖水體礦化度信息指數(shù)SIWVII (即WorldView-II 第3 波段反射率/WorldView-II 第8 波段反射率)。為了充分體現(xiàn)WorldView-II 不同波段對(duì)湖水礦化度信息的綜合反映,將該礦化度指數(shù)進(jìn)行拓展,得到了WorldView-II 綜合礦化度信息指數(shù)T-SIWVII,即WorldView-II 第3 波段反射率/WorldView-II 第8 波段反射率、WorldView-II 第4 波段反射率/WorldView-II第8 波段反射率、WorldView-II 第1 波段反射率/WorldView-II 第 8 波段反射率。WorldView-II 多波段綜合礦化度指數(shù)相較單一波段礦化度指數(shù)可反映更為豐富的礦化度信息,為后續(xù)湖水礦化度信息識(shí)別和提取奠定了堅(jiān)實(shí)基礎(chǔ)。

    2.4 鹽湖水體礦化度信息分離

    WorldView-II 遙感數(shù)據(jù)獲取的鹽湖水體信息包含了深度信息、礦化度信息以及可能存在的懸浮物信息等,尤以湖水深度信息對(duì)WorldView-II 數(shù)據(jù)光譜反射率影響較大。因此,為了剝離出WorldView-II 綜合礦化度指數(shù)所包含的深度及其他信息,本文基于主成分分析法對(duì)WorldView-II 綜合礦化度指數(shù)進(jìn)行了信息分離,獲得了PC1、PC2 和PC3 3 個(gè)主成分影像(圖3)。因PC3 影像主要為水體雜質(zhì)信息和數(shù)據(jù)自身噪聲,本文不做闡述。

    圖2 WV-II 彩色影像不同礦化度湖水點(diǎn)位選取圖(a)、選取點(diǎn)位光譜曲線圖(b)Fig.2 The sampling location of different salinity brine in WV-II color image (a)and spectral curve of samples (b)

    圖3 WorldView-Ⅱ綜合礦化度信息分離PC1 主分量圖(a)和PC2 主分量圖(b)Fig.3 Extracted salinity information in PC1(a)and PC2 (b)of WV-2 PCA image

    PC2 影像呈現(xiàn)出鹽湖湖岸淺水區(qū)色調(diào)淺、亮度高,深水區(qū)色調(diào)深、亮度低的特性,主要反映了湖水水深變化以及水底物質(zhì)反射所引起的圖像色調(diào)和亮度值變化信息。PC1 影像集中了WorldView-II 綜合礦化度指數(shù)的主要信息,影像清晰且亮度值變化有明顯的層次性。鹽湖東部和東北部亮度值總體偏高,反映了湖水因淡水補(bǔ)給量少、蒸發(fā)量大而形成的高礦化度,湖水西部和南部亮度值總體偏低,反映了湖水因淡水補(bǔ)給量大而形成的低礦化度。為了進(jìn)一步論證PC1 影像亮度值非湖水深度信息的反映,本文引入實(shí)測(cè)的尕斯庫勒鹽湖水深等深線圖進(jìn)行對(duì)比(圖4)[14]。由圖4 可見,鹽湖湖岸的水深值大致相同,約為40 cm 左右。越靠近鹽湖中央水深越深,最深約為120 cm 左右。因此,如果PC1 影像東部和北部湖區(qū)亮度值高代表了鹽湖深度,則依據(jù)鹽湖水深等深線分布,PC1 影像西部和南部湖區(qū)也應(yīng)同東部和北部湖區(qū)一樣高亮顯示。但是,PC1 影像西部和南部湖區(qū)則主要呈現(xiàn)低亮度值。因此,本文分離的PC1 影像亮度值的高低非鹽湖深度信息而為湖水的礦化度信息。

    圖4 尕斯庫勒鹽湖人工實(shí)測(cè)等深線圖[14]Fig.4 Contour map of manually measured water depth in Gasikule salt lake

    2.5 鹽湖礦化度含量定量反演

    以1.2.2 節(jié)所述的20 個(gè)實(shí)測(cè)鹽湖水礦化度樣本為建模數(shù)據(jù),利用最小二乘法對(duì)自變量20 個(gè)實(shí)測(cè)鹽湖水礦化度點(diǎn)位的PC1 影像值與其對(duì)應(yīng)的實(shí)測(cè)礦化度值進(jìn)行線性回歸擬合,得到了一次直線函數(shù)、對(duì)數(shù)曲線函數(shù)、指數(shù)曲線函數(shù)擬合模型(表2)。由表2 中可見,一次直線函數(shù)決定系數(shù)(R2)最大,擬合效果最好,故選其作為礦化度定量反演模型。基于ENVI 4.8 軟件進(jìn)行一次直線反演函數(shù)運(yùn)算,得到了鹽湖礦化度定量反演灰度影像圖,再對(duì)其進(jìn)行彩色密度分割,得到了尕斯庫勒鹽湖礦化度定量反演結(jié)果圖(圖5)。

    2.6 鹽湖礦化度含量定量反演模型精度評(píng)價(jià)與誤差分析

    為評(píng)價(jià)上節(jié)建立的尕斯庫勒鹽湖礦化度含量定量反演模型的反演精度,選取了1.2.2節(jié)所述的20 個(gè)實(shí)測(cè)尕斯庫勒鹽湖水礦化度值與模型反演值進(jìn)行統(tǒng)計(jì)分析(表3)。據(jù)表3 可以得出,礦化度反演值與實(shí)測(cè)值相比較,絕對(duì)誤差平均值為9.33 g/L,反演精度平均值達(dá)到96.61%。由圖5 可見,絕對(duì)誤差大于10 g/L的水樣點(diǎn)分別是3、5、8、13、14、15、18和20 號(hào)點(diǎn),除14 號(hào)點(diǎn)位于鹽湖中央深水區(qū)外,其余水樣點(diǎn)位均分布在靠近湖岸的淺水區(qū)。而湖岸淺水區(qū)受到外部淡水不定量的動(dòng)態(tài)補(bǔ)給,湖水礦化度的波動(dòng)變化也比較大。此外,水樣采集間距較大、取樣部位不完全為湖表水、WorldView-II 遙感影像與取樣單元空間分辨率的不匹配、室內(nèi)化學(xué)分析測(cè)試誤差等眾多因素,也是造成礦化度定量反演產(chǎn)生誤差的主要原因。

    表2 鹽湖礦化度定量反演回歸擬合函數(shù)Table 2 Regressing fitting function of quantitatively inversed salinity

    圖5 基于WorldView-II 數(shù)據(jù)的尕斯庫勒鹽湖礦化度定量反演圖Fig.5 Quantitatively inversed salinity map of Gasikule salt lake based on WV-2 data

    表3 尕斯庫勒鹽湖礦化度定量反演模型精度評(píng)價(jià)表Table 3 Accuracy evaluation on the salinity quantitative inversion model for Gasikule salt lake

    3 結(jié)論與討論

    本文基于WorldView-II 高分辨率衛(wèi)星遙感數(shù)據(jù)開展了鹽湖礦化度含量定量反演研究,并以尕斯庫勒鹽湖為例進(jìn)行了礦化度定量反演和精度驗(yàn)證,得出結(jié)論如下:

    1)基于WorldView-II 第1、3、4、8 波段建立了可表征鹽湖礦化度信息的WorldView-II 綜合礦化度指數(shù),利用主成分分析法實(shí)現(xiàn)了對(duì)WorldView-II 綜合礦化度指數(shù)中鹽湖礦化度和深度信息的分離,獲取的PC1主成分影像亮度值即可表征湖水的礦化度高低。

    2)一次直線函數(shù)對(duì)基于WorldView-II 數(shù)據(jù)構(gòu)建鹽湖礦化度定量反演模型的效果最好,經(jīng)精度評(píng)價(jià),模型總體反演精度達(dá)到96.61%,反演成功。

    3)影響鹽湖礦化度定量反演模型精度的因素主要為湖岸淺水區(qū)淡水補(bǔ)給、水樣采集間距較大、取樣部位不完全為湖表水、WorldView-II 遙感影像與取樣單元空間分辨率的不匹配、室內(nèi)化學(xué)分析測(cè)試誤差等。

    總之,由于受到各種因素影響,基于WorldView-II 遙感數(shù)據(jù)定量反演的鹽湖礦化度值與實(shí)際值有一定的誤差,但反演的礦化度變化趨勢(shì)及空間展布規(guī)律與實(shí)際值基本一致,完全可以應(yīng)用于鹽湖礦化度含量的定量識(shí)別,發(fā)揮遙感技術(shù)宏觀、快速、不受地域條件限制的綜合優(yōu)勢(shì),降低鹽湖區(qū)湖水取樣與分析的勞力和財(cái)力成本,為湖水礦化度分布規(guī)律分析和成礦預(yù)測(cè)提供重要的數(shù)據(jù)支撐。

    猜你喜歡
    礦化度鹽湖反射率
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    天空之境——新疆柴窩堡鹽湖
    具有顏色恒常性的光譜反射率重建
    山西運(yùn)城:冬日鹽湖色彩斑斕
    離子組成及礦化度對(duì)低礦化度水驅(qū)采收率的影響
    巖性油氣藏(2019年5期)2019-09-11 03:54:44
    如冰如雪,貌美鹽湖
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    西昆侖新發(fā)現(xiàn)鹽湖型鹵水硼鋰礦床
    利用綜合測(cè)井資料計(jì)算地下水等效NaCl溶液礦化度方法
    一区二区三区精品91| 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产| 久久久久久伊人网av| 街头女战士在线观看网站| 22中文网久久字幕| 欧美人与善性xxx| 国产精品不卡视频一区二区| 在线观看一区二区三区激情| 国产精品伦人一区二区| 国产美女午夜福利| 人妻少妇偷人精品九色| 精品一区二区三区视频在线| 男人和女人高潮做爰伦理| 精品少妇内射三级| 搡老乐熟女国产| 3wmmmm亚洲av在线观看| 夜夜看夜夜爽夜夜摸| 岛国毛片在线播放| 日韩,欧美,国产一区二区三区| 又黄又爽又刺激的免费视频.| 精品99又大又爽又粗少妇毛片| 国产精品一二三区在线看| 狂野欧美激情性bbbbbb| 国产美女午夜福利| 国产日韩欧美视频二区| 亚洲欧美一区二区三区国产| videossex国产| 最近手机中文字幕大全| 一区二区三区免费毛片| 亚洲欧美清纯卡通| 免费观看的影片在线观看| 亚洲va在线va天堂va国产| 91精品国产国语对白视频| 美女福利国产在线| 大码成人一级视频| 极品人妻少妇av视频| 全区人妻精品视频| 欧美xxⅹ黑人| 两个人的视频大全免费| 久久99精品国语久久久| 国产伦精品一区二区三区视频9| 国产成人免费观看mmmm| 少妇裸体淫交视频免费看高清| 亚洲欧美中文字幕日韩二区| 黑人巨大精品欧美一区二区蜜桃 | 日韩三级伦理在线观看| 日韩大片免费观看网站| 伦理电影免费视频| av女优亚洲男人天堂| 乱人伦中国视频| 亚洲久久久国产精品| 亚洲av国产av综合av卡| 久久影院123| 免费久久久久久久精品成人欧美视频 | 亚洲av二区三区四区| 欧美区成人在线视频| 大陆偷拍与自拍| 亚洲精品日韩av片在线观看| 国产精品国产三级国产av玫瑰| 最新中文字幕久久久久| 少妇裸体淫交视频免费看高清| 国产日韩欧美视频二区| 少妇被粗大的猛进出69影院 | 日韩中文字幕视频在线看片| 午夜激情久久久久久久| 亚洲国产最新在线播放| 亚洲精品第二区| 国产91av在线免费观看| 久久精品国产a三级三级三级| 国产精品不卡视频一区二区| 亚洲电影在线观看av| 日本wwww免费看| 久久久精品94久久精品| 久久ye,这里只有精品| 一级毛片电影观看| .国产精品久久| 成人亚洲精品一区在线观看| 中文字幕人妻丝袜制服| 男女免费视频国产| 在线看a的网站| 美女大奶头黄色视频| a级毛片在线看网站| 亚洲激情五月婷婷啪啪| 国产精品无大码| 亚洲熟女精品中文字幕| 少妇丰满av| 亚洲欧洲国产日韩| 亚洲怡红院男人天堂| av又黄又爽大尺度在线免费看| 成人特级av手机在线观看| 久久精品久久精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产精品嫩草影院av在线观看| 少妇 在线观看| 天堂俺去俺来也www色官网| 欧美成人精品欧美一级黄| 新久久久久国产一级毛片| 国产伦精品一区二区三区四那| 午夜福利,免费看| 国产黄片美女视频| 边亲边吃奶的免费视频| 久久精品久久久久久噜噜老黄| 少妇人妻久久综合中文| 久久久久视频综合| 国产精品一区二区在线不卡| 在线观看免费高清a一片| 18禁在线无遮挡免费观看视频| 欧美日韩av久久| 男人狂女人下面高潮的视频| 九九爱精品视频在线观看| tube8黄色片| 人妻 亚洲 视频| 黄色日韩在线| 美女主播在线视频| 我要看日韩黄色一级片| 六月丁香七月| 中文字幕制服av| 欧美一级a爱片免费观看看| 国产亚洲一区二区精品| 美女cb高潮喷水在线观看| 最近手机中文字幕大全| 国产乱来视频区| 超碰97精品在线观看| 国产色婷婷99| 亚洲欧美精品自产自拍| 亚洲国产精品专区欧美| 久久韩国三级中文字幕| 女的被弄到高潮叫床怎么办| 国产成人freesex在线| 久久99蜜桃精品久久| 日本wwww免费看| 爱豆传媒免费全集在线观看| 成人毛片60女人毛片免费| 国产一区二区三区av在线| 成人国产av品久久久| 在现免费观看毛片| 热re99久久精品国产66热6| 九九爱精品视频在线观看| 男人添女人高潮全过程视频| 亚洲欧美成人精品一区二区| 黄色日韩在线| 91aial.com中文字幕在线观看| 国产一区二区在线观看av| 激情五月婷婷亚洲| 六月丁香七月| 国产 精品1| 国产爽快片一区二区三区| 国产伦精品一区二区三区四那| av天堂中文字幕网| 亚洲激情五月婷婷啪啪| 亚洲国产毛片av蜜桃av| 亚洲图色成人| 成人午夜精彩视频在线观看| 久久久久久久久久久免费av| 男的添女的下面高潮视频| 日韩av免费高清视频| 久久国产亚洲av麻豆专区| 汤姆久久久久久久影院中文字幕| av福利片在线| 欧美精品一区二区免费开放| 欧美+日韩+精品| 亚洲av二区三区四区| 精品人妻偷拍中文字幕| 精品久久久久久电影网| 久久影院123| 亚洲四区av| 欧美 日韩 精品 国产| 一级毛片aaaaaa免费看小| 人妻制服诱惑在线中文字幕| 大陆偷拍与自拍| 日本av免费视频播放| 中国三级夫妇交换| 我要看日韩黄色一级片| 亚洲av免费高清在线观看| 91精品国产九色| 欧美变态另类bdsm刘玥| 亚洲经典国产精华液单| 亚洲精品aⅴ在线观看| 最新中文字幕久久久久| 欧美少妇被猛烈插入视频| 日日摸夜夜添夜夜爱| 最近的中文字幕免费完整| 乱码一卡2卡4卡精品| 一个人看视频在线观看www免费| 日韩av不卡免费在线播放| 亚洲国产av新网站| 视频中文字幕在线观看| 欧美xxxx性猛交bbbb| 久久精品国产a三级三级三级| 欧美变态另类bdsm刘玥| 国产免费福利视频在线观看| 亚洲精品日韩av片在线观看| 噜噜噜噜噜久久久久久91| 2021少妇久久久久久久久久久| 久久久久久久大尺度免费视频| 成年女人在线观看亚洲视频| 99国产精品免费福利视频| 高清午夜精品一区二区三区| 成人毛片a级毛片在线播放| 久久99精品国语久久久| 自拍偷自拍亚洲精品老妇| 午夜av观看不卡| 丰满人妻一区二区三区视频av| 一级二级三级毛片免费看| 黄色怎么调成土黄色| 丰满少妇做爰视频| 男人爽女人下面视频在线观看| 青春草国产在线视频| 丝袜在线中文字幕| 黄色欧美视频在线观看| 国产探花极品一区二区| 99精国产麻豆久久婷婷| 一区二区三区精品91| 国产一区亚洲一区在线观看| 免费观看无遮挡的男女| 国产黄色视频一区二区在线观看| 久久精品国产鲁丝片午夜精品| 一级爰片在线观看| 91精品一卡2卡3卡4卡| 在线观看av片永久免费下载| 国产精品一区二区性色av| 精品久久久久久久久亚洲| 有码 亚洲区| 美女内射精品一级片tv| 麻豆精品久久久久久蜜桃| 久久国产精品男人的天堂亚洲 | 黑人高潮一二区| 一区二区三区乱码不卡18| 一级av片app| 91精品一卡2卡3卡4卡| 亚洲av欧美aⅴ国产| 色哟哟·www| 国产美女午夜福利| 在线观看人妻少妇| 亚洲成人av在线免费| 国产成人免费观看mmmm| 最近手机中文字幕大全| 熟妇人妻不卡中文字幕| 国产精品久久久久久久电影| 自拍欧美九色日韩亚洲蝌蚪91 | 国产又色又爽无遮挡免| 久久99一区二区三区| 日本免费在线观看一区| 国产一区二区三区av在线| 免费久久久久久久精品成人欧美视频 | 久久久久久伊人网av| 亚洲精品久久午夜乱码| 久热久热在线精品观看| 亚洲欧美成人精品一区二区| 国产女主播在线喷水免费视频网站| 亚洲美女黄色视频免费看| 成人国产av品久久久| 男女免费视频国产| 国产一区二区在线观看av| 91久久精品国产一区二区三区| 建设人人有责人人尽责人人享有的| 欧美成人午夜免费资源| 国产伦精品一区二区三区四那| 亚洲av二区三区四区| 亚洲av中文av极速乱| 黑人猛操日本美女一级片| 99热6这里只有精品| 久久鲁丝午夜福利片| 精品99又大又爽又粗少妇毛片| av黄色大香蕉| 精品酒店卫生间| 永久免费av网站大全| 久热这里只有精品99| 欧美激情国产日韩精品一区| 国产免费一级a男人的天堂| 久久久久精品久久久久真实原创| 亚洲欧洲精品一区二区精品久久久 | 国产一区亚洲一区在线观看| 日日摸夜夜添夜夜爱| 国产日韩欧美视频二区| 国产在线一区二区三区精| 久久精品国产亚洲网站| 人人妻人人看人人澡| 中文天堂在线官网| 亚洲伊人久久精品综合| 色视频在线一区二区三区| 国模一区二区三区四区视频| 天堂俺去俺来也www色官网| 在线观看免费视频网站a站| 夜夜看夜夜爽夜夜摸| 26uuu在线亚洲综合色| 一本色道久久久久久精品综合| 人人澡人人妻人| av女优亚洲男人天堂| 欧美精品高潮呻吟av久久| 亚洲av日韩在线播放| 亚洲国产av新网站| 毛片一级片免费看久久久久| 三级国产精品片| 久久狼人影院| 最近中文字幕2019免费版| 欧美亚洲 丝袜 人妻 在线| 精品一区二区免费观看| 久久人人爽av亚洲精品天堂| 曰老女人黄片| 久久久a久久爽久久v久久| 免费av不卡在线播放| 久久精品国产自在天天线| 久久午夜综合久久蜜桃| 男人和女人高潮做爰伦理| 极品教师在线视频| 尾随美女入室| 美女大奶头黄色视频| 亚洲综合色惰| 亚洲国产色片| 国产色爽女视频免费观看| 91aial.com中文字幕在线观看| 免费观看的影片在线观看| 秋霞伦理黄片| 99热国产这里只有精品6| 国产精品99久久99久久久不卡 | 曰老女人黄片| 日本91视频免费播放| 国产精品偷伦视频观看了| 观看美女的网站| 自线自在国产av| freevideosex欧美| 欧美日韩国产mv在线观看视频| 日韩视频在线欧美| 成年美女黄网站色视频大全免费 | 日日爽夜夜爽网站| 久热这里只有精品99| 欧美日韩精品成人综合77777| 久久久久精品性色| 国产男女超爽视频在线观看| 黄色怎么调成土黄色| 成年女人在线观看亚洲视频| 蜜臀久久99精品久久宅男| 我要看黄色一级片免费的| 国产精品人妻久久久久久| 亚洲精品乱久久久久久| 麻豆乱淫一区二区| 免费看光身美女| 国产精品免费大片| 啦啦啦啦在线视频资源| 永久网站在线| av卡一久久| 亚洲性久久影院| 久久人人爽人人爽人人片va| 欧美成人精品欧美一级黄| 永久免费av网站大全| 91精品伊人久久大香线蕉| 成年女人在线观看亚洲视频| 国产伦精品一区二区三区视频9| 少妇 在线观看| 最近中文字幕2019免费版| 亚洲精品国产成人久久av| 国产一级毛片在线| 一级毛片黄色毛片免费观看视频| 狂野欧美激情性xxxx在线观看| av在线app专区| 精品久久久久久久久亚洲| 欧美日韩av久久| 麻豆成人av视频| 自拍偷自拍亚洲精品老妇| 高清欧美精品videossex| 国产伦理片在线播放av一区| 日韩三级伦理在线观看| 亚洲国产欧美日韩在线播放 | 中文字幕av电影在线播放| 少妇的逼水好多| 久久久久久久久大av| 热re99久久精品国产66热6| 亚洲成人av在线免费| av在线播放精品| 免费少妇av软件| 亚洲精品自拍成人| 校园人妻丝袜中文字幕| 精品一品国产午夜福利视频| 久久久久久伊人网av| 嘟嘟电影网在线观看| 免费观看在线日韩| 免费av中文字幕在线| 国产精品一二三区在线看| 亚洲经典国产精华液单| 一本—道久久a久久精品蜜桃钙片| 曰老女人黄片| 国产精品不卡视频一区二区| 欧美性感艳星| 狂野欧美激情性xxxx在线观看| 亚洲成人一二三区av| 国产精品.久久久| 日日撸夜夜添| 99热6这里只有精品| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产精品一区三区| 成年人午夜在线观看视频| 欧美日韩av久久| 亚洲在久久综合| 97超视频在线观看视频| 街头女战士在线观看网站| 久久综合国产亚洲精品| 九九爱精品视频在线观看| 国产成人精品久久久久久| 夫妻午夜视频| 久久久久久久大尺度免费视频| 国产成人精品久久久久久| 噜噜噜噜噜久久久久久91| 亚洲美女搞黄在线观看| 高清黄色对白视频在线免费看 | 丝瓜视频免费看黄片| 在线看a的网站| 9色porny在线观看| 香蕉精品网在线| 永久免费av网站大全| 在线观看美女被高潮喷水网站| 97在线人人人人妻| 亚洲精品,欧美精品| av国产久精品久网站免费入址| 久热这里只有精品99| 亚洲美女视频黄频| 国产精品99久久99久久久不卡 | 毛片一级片免费看久久久久| 国产日韩欧美在线精品| 国内少妇人妻偷人精品xxx网站| 天美传媒精品一区二区| 少妇精品久久久久久久| 久久99一区二区三区| 天堂8中文在线网| 乱人伦中国视频| 欧美精品一区二区大全| 18禁在线播放成人免费| 亚洲人成网站在线播| 大话2 男鬼变身卡| 在线观看美女被高潮喷水网站| 亚洲精品乱码久久久v下载方式| 高清午夜精品一区二区三区| av天堂久久9| 丝瓜视频免费看黄片| 日韩人妻高清精品专区| 亚洲精品久久午夜乱码| 午夜影院在线不卡| 亚洲av中文av极速乱| av卡一久久| 成年女人在线观看亚洲视频| 日韩精品有码人妻一区| 久久久久网色| 日本黄色日本黄色录像| 一本色道久久久久久精品综合| 国产中年淑女户外野战色| 岛国毛片在线播放| 国产亚洲欧美精品永久| videossex国产| 18+在线观看网站| 少妇人妻 视频| 天堂俺去俺来也www色官网| 久久国产精品大桥未久av | 男女无遮挡免费网站观看| 国产真实伦视频高清在线观看| 亚洲va在线va天堂va国产| 国产精品人妻久久久久久| 国产精品秋霞免费鲁丝片| 国产美女午夜福利| 女人久久www免费人成看片| 久久人人爽人人爽人人片va| 一级av片app| 亚洲欧美一区二区三区黑人 | 欧美精品人与动牲交sv欧美| 久久久久久久久久久丰满| 大陆偷拍与自拍| 青春草亚洲视频在线观看| 极品少妇高潮喷水抽搐| 日本与韩国留学比较| 国产老妇伦熟女老妇高清| 五月开心婷婷网| 狂野欧美激情性xxxx在线观看| 精品亚洲成国产av| 欧美bdsm另类| 国产69精品久久久久777片| 天堂中文最新版在线下载| 久久 成人 亚洲| 日韩成人伦理影院| 97超视频在线观看视频| 国产69精品久久久久777片| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 国产成人精品福利久久| 欧美日韩视频高清一区二区三区二| 久久狼人影院| 久久久久久伊人网av| 久久久精品94久久精品| 丝瓜视频免费看黄片| 亚洲精品,欧美精品| 简卡轻食公司| 免费观看在线日韩| 精品少妇内射三级| 国产探花极品一区二区| 国产亚洲最大av| 欧美高清成人免费视频www| 狂野欧美激情性xxxx在线观看| 少妇的逼水好多| 肉色欧美久久久久久久蜜桃| 国产一区二区在线观看av| 国产一级毛片在线| 制服丝袜香蕉在线| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 国产av精品麻豆| 成人18禁高潮啪啪吃奶动态图 | 观看美女的网站| 蜜桃在线观看..| 国产男女内射视频| 99热这里只有是精品50| 国产熟女午夜一区二区三区 | 日韩 亚洲 欧美在线| 高清午夜精品一区二区三区| 啦啦啦视频在线资源免费观看| 久久av网站| 伦理电影免费视频| 一级毛片aaaaaa免费看小| 美女国产视频在线观看| 女的被弄到高潮叫床怎么办| 少妇人妻一区二区三区视频| 麻豆精品久久久久久蜜桃| 22中文网久久字幕| 国产片特级美女逼逼视频| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 日日摸夜夜添夜夜添av毛片| 2018国产大陆天天弄谢| 99热这里只有是精品50| 丰满人妻一区二区三区视频av| 亚洲精品一区蜜桃| 最新中文字幕久久久久| 亚洲av综合色区一区| 免费高清在线观看视频在线观看| 国产精品嫩草影院av在线观看| 九九在线视频观看精品| 久久午夜福利片| 久久国产精品大桥未久av | 久久av网站| 狂野欧美激情性xxxx在线观看| 赤兔流量卡办理| 色哟哟·www| 男女无遮挡免费网站观看| 精品午夜福利在线看| 国产一区有黄有色的免费视频| 女人精品久久久久毛片| 久久久久久人妻| 观看av在线不卡| 亚洲精品中文字幕在线视频 | 美女xxoo啪啪120秒动态图| 国产又色又爽无遮挡免| 欧美国产精品一级二级三级 | 亚洲国产欧美在线一区| 91成人精品电影| 日韩三级伦理在线观看| 亚洲av成人精品一区久久| 亚洲一级一片aⅴ在线观看| 乱码一卡2卡4卡精品| 美女国产视频在线观看| 久久精品国产亚洲av天美| 日韩不卡一区二区三区视频在线| 狠狠精品人妻久久久久久综合| av福利片在线| 99热这里只有精品一区| av福利片在线| 晚上一个人看的免费电影| 哪个播放器可以免费观看大片| 国产精品久久久久久久久免| 深夜a级毛片| 久久久久久久大尺度免费视频| 少妇被粗大猛烈的视频| 蜜桃久久精品国产亚洲av| 成人国产麻豆网| 蜜桃久久精品国产亚洲av| 一区二区av电影网| 亚洲精品一区蜜桃| 亚洲欧美成人综合另类久久久| 国产精品伦人一区二区| 精品久久久久久电影网| 日韩电影二区| 久久午夜综合久久蜜桃| 久久久久久久久久成人| 天美传媒精品一区二区| 伊人久久精品亚洲午夜| 热re99久久国产66热| 大陆偷拍与自拍| 国产老妇伦熟女老妇高清| 亚洲真实伦在线观看| 纯流量卡能插随身wifi吗| 日本vs欧美在线观看视频 | 国产成人91sexporn| 久久久久久久亚洲中文字幕| 免费av中文字幕在线| 熟女av电影| 免费大片18禁| 一个人看视频在线观看www免费| 好男人视频免费观看在线| 国产免费一级a男人的天堂| 性色av一级| 菩萨蛮人人尽说江南好唐韦庄| 乱人伦中国视频| 美女国产视频在线观看| 午夜视频国产福利| 免费观看无遮挡的男女| 制服丝袜香蕉在线| 日本黄色片子视频| 亚洲国产精品专区欧美| 午夜免费观看性视频| 久久精品久久精品一区二区三区| 精品国产一区二区三区久久久樱花| av女优亚洲男人天堂| 3wmmmm亚洲av在线观看| 精品亚洲成国产av| 亚洲国产精品一区三区| 人人妻人人澡人人爽人人夜夜| 黄色视频在线播放观看不卡| 国产极品粉嫩免费观看在线 |