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

    環(huán)境衛(wèi)星遙感反演在水源地水質(zhì)應(yīng)急監(jiān)測(cè)預(yù)警中的精度探討及修正方法研究

    2020-03-11 03:35:58
    關(guān)鍵詞:水源地反射率回歸方程

    巴 諾

    (遼寧省本溪水文局,遼寧 本溪 117000)

    1 概述

    水源地水質(zhì)環(huán)境狀況是當(dāng)前社會(huì)關(guān)注的熱點(diǎn)和焦點(diǎn)。水環(huán)境監(jiān)測(cè)是保障水源地水環(huán)境狀況的重要措施。近些年來,水源地水質(zhì)狀況頻發(fā),特別是富營(yíng)養(yǎng)化影響下的水源地“水華”現(xiàn)成呈現(xiàn)多發(fā)、頻發(fā)的態(tài)勢(shì),急需要對(duì)水源地環(huán)境狀況進(jìn)行有效的應(yīng)急監(jiān)測(cè)和預(yù)警。目前,對(duì)于水源地的應(yīng)急監(jiān)測(cè)主要通過采集水樣,通過及時(shí)測(cè)定水樣數(shù)據(jù)來對(duì)水質(zhì)狀況進(jìn)行監(jiān)測(cè),這種方式的優(yōu)點(diǎn)在于測(cè)定數(shù)據(jù)精度較高,缺點(diǎn)在于很難對(duì)水源地水質(zhì)狀況進(jìn)行全覆蓋監(jiān)測(cè),存在以點(diǎn)代面的局限,此外很難實(shí)現(xiàn)水質(zhì)預(yù)警,因而很難實(shí)現(xiàn)對(duì)水源地水環(huán)境狀況的提前保護(hù)。當(dāng)前,環(huán)境衛(wèi)星遙感數(shù)據(jù)逐步在水環(huán)境監(jiān)測(cè)中得到應(yīng)用,但是由于環(huán)境衛(wèi)星遙感數(shù)據(jù)的精度和時(shí)效性,其中水源地水質(zhì)應(yīng)急監(jiān)測(cè)預(yù)警中應(yīng)用還較少,為此本文以遼寧東部某水源地為研究實(shí)例,結(jié)合7期環(huán)境衛(wèi)星遙感數(shù)據(jù),分析其在水源地應(yīng)急監(jiān)測(cè)的預(yù)警精度,并通過建立回歸方程對(duì)反演指標(biāo)進(jìn)行修正。研究成果對(duì)于水源地應(yīng)急遙感監(jiān)測(cè)、預(yù)警具有重要的參考價(jià)值。

    2 環(huán)境衛(wèi)星遙感影像研究方法

    2.1 數(shù)據(jù)獲取方法

    為保障水環(huán)境應(yīng)急監(jiān)測(cè)對(duì)數(shù)據(jù)的時(shí)效性和精度的需求,當(dāng)前應(yīng)用較為廣泛的環(huán)境衛(wèi)星遙感數(shù)據(jù)主要有8個(gè),具體見表1。這些衛(wèi)星遙感數(shù)據(jù)為水質(zhì)應(yīng)急監(jiān)測(cè)提供了較好的數(shù)據(jù)支撐。本文選擇HJ-1A/1B環(huán)境衛(wèi)星數(shù)據(jù)作為分析數(shù)據(jù),該數(shù)據(jù)的下載平臺(tái)見參考文獻(xiàn)[16],該環(huán)境衛(wèi)星包含藍(lán)、綠、紅和近紅外共4個(gè)波段,空間分辨率可達(dá)30m。本文還應(yīng)用Landsat8 OLI遙感影像作為基本影響數(shù)據(jù),對(duì)環(huán)境衛(wèi)星遙感數(shù)據(jù)進(jìn)行幾何修正處理。為保證遙感影響數(shù)據(jù)的同步性,選擇環(huán)境衛(wèi)星遙感過境時(shí)間為2天的數(shù)據(jù)作為研究數(shù)據(jù)。

    表1 常用的環(huán)境衛(wèi)星遙感數(shù)據(jù)

    2.2 數(shù)據(jù)預(yù)處理方法

    (1)輻射定標(biāo)

    由于不同衛(wèi)星傳感器對(duì)反演物體波段的響應(yīng)效果不同,使得相同輸入項(xiàng)后得到不同的輸出值,使得遙感圖像上出現(xiàn)條紋現(xiàn)象,因此需要對(duì)這種影響進(jìn)行輻射定標(biāo)的處理,輻射定標(biāo)系數(shù)的計(jì)算方程為:

    L=DN/A+L0

    (1)

    式中,L0—偏移量,mm;A—定標(biāo)的增益系數(shù)。

    各遙感數(shù)據(jù)輻射定標(biāo)系數(shù)見表2。

    表2 常用環(huán)境遙感衛(wèi)星不同參量的輻射定標(biāo)系數(shù)

    (2)幾何校正

    由于遙感器的狀態(tài)以及運(yùn)行速度會(huì)使得遙感影像數(shù)據(jù)發(fā)生不同程度的形變,使得遙感影像數(shù)據(jù)出現(xiàn)不同程度的幾何變形,幾何校正是遙感影像數(shù)據(jù)預(yù)處理的重要環(huán)節(jié)。本文應(yīng)用Landsat8 OLI作為基準(zhǔn)影像數(shù)據(jù)進(jìn)行幾何校正。

    (3)大氣校正

    由于大氣散射所引起的地表反射率影像使得遙感圖像中出現(xiàn)相應(yīng)的誤差,為消除這類誤差,需要對(duì)獲取的遙感影響數(shù)據(jù)進(jìn)行大氣校正。大氣校正需要輸入的主要參數(shù)包括圖像中心經(jīng)緯度、傳感器飛行高度、圖像區(qū)域平均海拔、圖像像素大小、傳感器種類、溶膠反演方法、氣溶膠模型、成像時(shí)間等。

    2.3 水質(zhì)反演原理

    通過計(jì)算環(huán)境衛(wèi)星遙感數(shù)據(jù)中各波段總的輻射亮度值來計(jì)算各水質(zhì)指標(biāo)的反演值,計(jì)算方程為:

    Lt=LP+LS+LV+Lb

    (2)

    式中,LP—水體表面下太陽(yáng)光和天空輻射;Ls—水體表層輻射;Lv—達(dá)到水體表層的太陽(yáng)輻射以及水體中有機(jī)物和無機(jī)物與水體相互作用的輻射;Lb—水體底部反射的太陽(yáng)光輻射。

    由于水體中不同物質(zhì)成分受到水體吸收和散射的影響,不同水體有機(jī)或者無機(jī)物質(zhì)在一定波段內(nèi)的反射率不同,結(jié)合該反射率是推求水源地水質(zhì)各指標(biāo)的重要基礎(chǔ)。

    2.4 水質(zhì)綜合預(yù)警原理

    結(jié)合環(huán)境衛(wèi)星遙感影像數(shù)據(jù),對(duì)各水質(zhì)指標(biāo)進(jìn)行綜合反演,計(jì)算不同指標(biāo)的綜合指數(shù),再綜合各水質(zhì)指標(biāo)數(shù)進(jìn)行加權(quán),得到綜合預(yù)警指標(biāo)值,綜合預(yù)警指數(shù)方程為:

    (3)

    其中各指標(biāo)權(quán)重計(jì)算方程為:

    (4)

    式中,TLI(∑)—綜合預(yù)警值;Wj—各綜合指標(biāo)的權(quán)重值;TLI(j)—各反演水質(zhì)指標(biāo)修正值;m—評(píng)價(jià)指標(biāo)個(gè)數(shù);rij—各指標(biāo)相關(guān)系數(shù)。

    3 研究實(shí)例

    3.1 水源地概況及水樣監(jiān)測(cè)點(diǎn)分布

    以遼寧東部某水源地作為研究實(shí)例,該水源地為城市居民生活飲用水重要的水源地,每年的城市生活飲用水供水量為1.3億m3,水源地控制的面積為2795km2,庫(kù)區(qū)降水主要來源于夏季,庫(kù)區(qū)以上多年平均降水量為700~900mm,多年平均來水量為11.3×108m3,有效調(diào)節(jié)庫(kù)容為2.35億m3。近些年來,水源地加大水環(huán)境保護(hù),水環(huán)境狀況得到逐步改善,為有效監(jiān)測(cè)水源地水質(zhì)狀況,在水源地庫(kù)區(qū)內(nèi)布設(shè)6個(gè)水樣監(jiān)測(cè)點(diǎn)。

    3.2 數(shù)據(jù)預(yù)處理

    結(jié)合前面敘述的遙感影像獲取數(shù)據(jù),采用輻射定標(biāo)系數(shù)方法、幾何校正、大氣校正方法對(duì)遙感影像數(shù)據(jù)進(jìn)行預(yù)處理,其中選用的遙感影像數(shù)據(jù)輻射定標(biāo)系數(shù)見表3,大氣校正所輸入的參數(shù)見表4。

    表3 選用的遙感影像數(shù)據(jù)的輻射定標(biāo)系數(shù)

    表4 遙感影像數(shù)據(jù)大氣校正輸入的參數(shù)

    3.3 波段組合的選擇

    水質(zhì)指標(biāo)反演精度主要取決于選擇的波段,波段敏感程度對(duì)水質(zhì)監(jiān)測(cè)十分重要,HJ1A-CCD衛(wèi)星遙感數(shù)據(jù)共包含4個(gè)波段,結(jié)合Pearson相關(guān)分析法對(duì)各波段及波段組合及相對(duì)應(yīng)的水質(zhì)指標(biāo)的相關(guān)性進(jìn)行統(tǒng)計(jì)分析,分析結(jié)果見表5。

    表5 各波段及波段組合與相對(duì)應(yīng)的水質(zhì)指標(biāo)的相關(guān)性檢驗(yàn)結(jié)果

    注:**表示為在0.01水平(雙側(cè))上顯著相關(guān);*表示為在0.05水平(雙側(cè))上顯著相關(guān)。

    從相關(guān)性統(tǒng)計(jì)分析結(jié)果可看出,組合波段與水質(zhì)指標(biāo)的相關(guān)度均高于單一波段下的相關(guān)系數(shù)。因此一般選用與水質(zhì)指標(biāo)相關(guān)度較高的波段組合作為反演模型。從表中可看出,葉綠素的濃度與波段組合B2-B1的相關(guān)度最高,相關(guān)度為0.616,且通過雙側(cè)0.01水平的相關(guān)顯著性檢驗(yàn)。而與透明度相關(guān)性最高的組合波段為(B1+B2)/(B3+B4),其相關(guān)系數(shù)為0.413,也通過雙側(cè)0.01水平的相關(guān)顯著性檢驗(yàn)。懸浮物與B2-B1的波段組合相關(guān)系數(shù)最大,相關(guān)度為0.635,通過雙側(cè)0.05水平的相關(guān)顯著性檢驗(yàn)。對(duì)于富營(yíng)養(yǎng)指數(shù)而言,其相關(guān)度最高的波段組合為(B1+B2)/(B3+B4),相關(guān)系數(shù)為-0.325,相關(guān)性通過雙側(cè)0.01水平的相關(guān)顯著性檢驗(yàn)。從波段分析結(jié)果還可看出,水體污染物指標(biāo)濃度越大,其波段反射率越大,通過波段組合優(yōu)先分析,波段反射率在460~750nm適合于葉綠素、懸浮物、透明度以及富營(yíng)養(yǎng)化指數(shù)的反演計(jì)算。

    3.4 波段反演回歸方程的建立

    在波段組合相關(guān)性分析的基礎(chǔ)上,確定了各反演水質(zhì)指標(biāo)的最優(yōu)波段,結(jié)合所選取的環(huán)境衛(wèi)星遙感數(shù)據(jù),提取各最優(yōu)波段的反射率數(shù)據(jù),通過反演方程計(jì)算各組合波段的反演值,作為自變量x,各水質(zhì)反演指標(biāo)作為應(yīng)變量y,建立各水質(zhì)指標(biāo)的回歸方程,見表6,回歸方程建立的依據(jù)如圖1所示。

    表6 各水質(zhì)指標(biāo)反演回歸方程的建立結(jié)果

    以葉綠素波段組合下反射率值為其回歸方程的自變量,以其水質(zhì)指標(biāo)測(cè)定值為其應(yīng)變量,建立葉綠素的反演回歸方程,從其建立的回歸方程可看出,較好的擬合了中、下部區(qū)域的散點(diǎn),相關(guān)系數(shù)為0.6167,相關(guān)度較高,可以應(yīng)用該回歸方程來建立遙感模型,對(duì)水源地葉綠素指標(biāo)進(jìn)行反演。建立懸浮物最優(yōu)波段組合下反射率下的回歸方程,從圖1中可看出,其建立的回歸方程較好的擬合了中下部散點(diǎn)數(shù)據(jù),相關(guān)系數(shù)達(dá)到0.5226,也可用來進(jìn)行水源地懸浮物的遙感反演計(jì)算。從透明度的反演值可看出,以透明度組合波段下反射率和透明度測(cè)定值建立的回歸方程的系數(shù)可達(dá)到0.5139,可用來進(jìn)行透明度的反演計(jì)算。從富營(yíng)養(yǎng)化指數(shù)的反演回歸方程可看出,其相關(guān)系數(shù)達(dá)到0.4965,相關(guān)系數(shù)在各指標(biāo)中最低,主要是因?yàn)楦粻I(yíng)養(yǎng)化指數(shù)為綜合指數(shù),其波段組合反射率影響因素較多,因此其建立的回歸方程相關(guān)系數(shù)較低,但從圖1中可看出,其回歸方程的相關(guān)線較好的擬合了其散點(diǎn),可用來進(jìn)行富營(yíng)養(yǎng)指數(shù)的反演計(jì)算。

    圖1 各水質(zhì)指標(biāo)的波段組合值與反演值的回歸方程

    3.5 反演修正后的精度對(duì)比分析

    結(jié)合建立的反演回歸方程,構(gòu)建水源地各水質(zhì)指標(biāo)遙感監(jiān)測(cè)模型,對(duì)不同指標(biāo)進(jìn)行反演計(jì)算,并結(jié)合各采樣點(diǎn)實(shí)測(cè)的水質(zhì)監(jiān)測(cè)值,對(duì)比分析修正前后的反演計(jì)算誤差,分析結(jié)果見表7—10。

    表7 葉綠素指標(biāo)修正前后的反演誤差對(duì)比結(jié)果

    表8 懸浮物指標(biāo)修正前后的反演誤差對(duì)比結(jié)果

    表10 富營(yíng)養(yǎng)化指數(shù)修正前后的反演誤差對(duì)比結(jié)果 單位:%

    從表7—10中可看出,修正前后水源地葉綠素反演誤差得到較為明顯的改善,修正前期反演誤差在40.41%~43.03%之間,而通過建立的回歸方程進(jìn)行修正后,其葉綠素的反演誤差在21.09%~23.88%之間,相比于修正前,誤差降低22.7%,而從懸浮物的反演誤差對(duì)比結(jié)果可看出,修正前其誤差在±40%以上,而采用建立的回歸方程進(jìn)行修正后,其懸浮物的反演誤差在±25%以內(nèi),反演精度得到有效提高。從透明度修正前后的反演誤差可看出,修正前其反演誤差最高值為48.15%,而同一個(gè)采樣點(diǎn),其誤差降低到22.22%,可見通過建立回歸方程對(duì)其反演進(jìn)行修正,可有效降低水源地水質(zhì)指標(biāo)的反演誤差。

    3.6 各水質(zhì)指標(biāo)空間反演結(jié)果

    在反演修正的基礎(chǔ)上,結(jié)合環(huán)境衛(wèi)星遙感數(shù)據(jù)對(duì)不同指標(biāo)同一期數(shù)據(jù)進(jìn)行空間反演,反演結(jié)果如圖2所示。

    從圖2中可看出,各水質(zhì)指標(biāo)反演空間值總體呈現(xiàn)從上游~下游逐步遞減的變化,污染物濃度較高的區(qū)域主要位于其上游段,這主要是因?yàn)槭艿缴嫌蝸硭挠绊?,上游降雨徑流中的攜帶的污染物進(jìn)行水源地,使得其上游污染物濃度偏高,而下游受到生態(tài)前置庫(kù)的影響,下游水體污染物濃度得到控制。從葉綠素濃度空間分布可看出,45.0~60.0ug/L的區(qū)域比例占總面積比例為8.5%,均位于水源地上游,而0.5~1.0ug/L區(qū)域占總面積比例為45.2%,位于水庫(kù)的下游主要水源區(qū)。從懸浮物的空間分析結(jié)果可看出,懸浮物濃度在200~230mg/l的比例為15.7%,也主要位于水源地上游段。從透明度分析結(jié)果可看出,透明度和懸浮物的空間反演分布較為一致,但中等指標(biāo)面積比例較懸浮物比例有所減少。從富營(yíng)養(yǎng)化指數(shù)可看出,由于其實(shí)葉綠素、透明度、懸浮物的綜合指數(shù)構(gòu)建,因此其反演空間變化是其他幾個(gè)指標(biāo)的綜合反演,因此空間分布具有一定的相似性,從分布值可看出,富營(yíng)養(yǎng)指數(shù)在60%~70%之間的比例為14.1%。

    3.7 水質(zhì)綜合預(yù)警研究

    由于當(dāng)前環(huán)境衛(wèi)星遙感數(shù)據(jù)的時(shí)效性得到顯著提高,可以結(jié)合衛(wèi)星遙感數(shù)據(jù)對(duì)水源地進(jìn)行多指標(biāo)的綜合預(yù)警,本文構(gòu)建了綜合預(yù)警的指標(biāo)等級(jí),見表11,并結(jié)合實(shí)例對(duì)水源地進(jìn)行綜合預(yù)警,結(jié)果見表12。

    圖2 各水質(zhì)指標(biāo)空間反演結(jié)果

    表11 水源地水質(zhì)綜合分級(jí)預(yù)警值

    表12 水源地水質(zhì)綜合預(yù)警結(jié)果

    從表11中可看出,主要分為5個(gè)等級(jí),不同等級(jí)對(duì)應(yīng)其相應(yīng)的預(yù)警指數(shù),通過環(huán)境衛(wèi)星遙感數(shù)據(jù),可提取逐月不同水質(zhì)指標(biāo)對(duì)應(yīng)的組合波段的反射率,通過構(gòu)建的回歸方程對(duì)各指標(biāo)進(jìn)行反演修正,得到不同指標(biāo)的反演濃度,通過設(shè)定權(quán)重,結(jié)合綜合預(yù)警指數(shù)計(jì)算方法,對(duì)其綜合預(yù)警指數(shù)進(jìn)行計(jì)算,實(shí)現(xiàn)不同程度的預(yù)警,表12為研究水源地綜合預(yù)警的一個(gè)實(shí)例,結(jié)合水源地2017年9月到2018年10月的環(huán)境衛(wèi)星遙感數(shù)據(jù),通過組合波段反射率的選取,計(jì)算其預(yù)警指數(shù),并結(jié)合表11確定不同月份的預(yù)警等級(jí),該水質(zhì)預(yù)警可有效實(shí)現(xiàn)水源地水環(huán)境的生態(tài)保護(hù)。

    4 結(jié)語(yǔ)

    (1)水體污染物指標(biāo)濃度越大,其波段反射率越大,通過波段組合優(yōu)先分析,波段反射率在460~750nm適合于葉綠素、懸浮物、透明度以及富營(yíng)養(yǎng)化指數(shù)的反演計(jì)算。

    (2)結(jié)合環(huán)境衛(wèi)星遙感數(shù)據(jù)可實(shí)現(xiàn)區(qū)域全覆蓋應(yīng)急反演監(jiān)測(cè),解決傳統(tǒng)監(jiān)測(cè)方法以點(diǎn)代面的局限,可在水源地水生態(tài)監(jiān)測(cè)中得到推廣和應(yīng)用。

    (3)基于衛(wèi)星反演數(shù)據(jù)建立的綜合預(yù)警方案,可實(shí)現(xiàn)水源地動(dòng)態(tài)預(yù)警,具體應(yīng)用中可結(jié)合水源地實(shí)際情況對(duì)預(yù)警等級(jí)及指標(biāo)權(quán)重進(jìn)行調(diào)整。

    (4)本文只針對(duì)水源地主要4個(gè)影響指標(biāo)反演精度進(jìn)行探討,其他指標(biāo)比如重金屬、水溫還需要以后進(jìn)行深入研究。

    猜你喜歡
    水源地反射率回歸方程
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    具有顏色恒常性的光譜反射率重建
    采用直線回歸方程預(yù)測(cè)桑癭蚊防治適期
    線性回歸方程的求解與應(yīng)用
    線性回歸方程要點(diǎn)導(dǎo)學(xué)
    走進(jìn)回歸分析,讓回歸方程不再是你高考的絆腳石
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    基于因果分析圖法的飲用水源地保護(hù)探討
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟(jì)損失估算
    日本av手机在线免费观看| 一区二区三区精品91| 黑人欧美特级aaaaaa片| 欧美亚洲日本最大视频资源| 男女床上黄色一级片免费看| 女警被强在线播放| 美女国产高潮福利片在线看| 久久精品人人爽人人爽视色| 我的亚洲天堂| 无遮挡黄片免费观看| 人人妻人人澡人人看| 亚洲精品久久成人aⅴ小说| 欧美激情久久久久久爽电影 | 午夜激情av网站| 波多野结衣一区麻豆| 男女无遮挡免费网站观看| 97人妻天天添夜夜摸| 国产成人啪精品午夜网站| 狠狠狠狠99中文字幕| 精品国产一区二区三区四区第35| 欧美日韩精品网址| 亚洲国产中文字幕在线视频| 一区二区三区精品91| 大香蕉久久网| 伊人久久大香线蕉亚洲五| 国产精品国产高清国产av | 高清黄色对白视频在线免费看| 乱人伦中国视频| 999精品在线视频| 老熟妇乱子伦视频在线观看| 亚洲精品国产一区二区精华液| 国产成人精品无人区| 久久狼人影院| 色婷婷久久久亚洲欧美| 亚洲国产欧美在线一区| 狠狠婷婷综合久久久久久88av| 精品视频人人做人人爽| 欧美日韩av久久| 两个人看的免费小视频| 黄色片一级片一级黄色片| 精品久久久久久久毛片微露脸| 亚洲成人免费av在线播放| 一级毛片精品| 午夜精品久久久久久毛片777| 性高湖久久久久久久久免费观看| 99久久人妻综合| 久久精品亚洲熟妇少妇任你| 久久久久久久久免费视频了| 精品国产一区二区久久| 欧美激情极品国产一区二区三区| 免费看a级黄色片| 午夜福利视频在线观看免费| 成人精品一区二区免费| 亚洲人成伊人成综合网2020| 最近最新中文字幕大全免费视频| 日韩欧美三级三区| 国产日韩欧美亚洲二区| 欧美日韩精品网址| 国产aⅴ精品一区二区三区波| 免费观看a级毛片全部| 成在线人永久免费视频| 国产在线观看jvid| 国产精品1区2区在线观看. | 蜜桃在线观看..| 久久人人97超碰香蕉20202| 欧美国产精品va在线观看不卡| 午夜成年电影在线免费观看| 中国美女看黄片| 精品乱码久久久久久99久播| 激情视频va一区二区三区| 国产av国产精品国产| 91老司机精品| 成人手机av| 香蕉丝袜av| 亚洲男人天堂网一区| 乱人伦中国视频| 建设人人有责人人尽责人人享有的| 丰满少妇做爰视频| 伦理电影免费视频| 人人妻人人澡人人看| 一级片免费观看大全| 19禁男女啪啪无遮挡网站| 少妇精品久久久久久久| 手机成人av网站| 国产一区二区 视频在线| 国产色视频综合| 满18在线观看网站| 亚洲av美国av| 精品亚洲成a人片在线观看| 9热在线视频观看99| 久久 成人 亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 国产福利在线免费观看视频| 国产精品二区激情视频| 婷婷成人精品国产| 精品高清国产在线一区| 久久午夜综合久久蜜桃| 国产极品粉嫩免费观看在线| 国产精品av久久久久免费| 久久av网站| 国产成人系列免费观看| 成人特级黄色片久久久久久久 | 免费在线观看日本一区| 国产单亲对白刺激| 美女主播在线视频| 视频在线观看一区二区三区| 国产精品一区二区在线观看99| 精品一品国产午夜福利视频| 午夜福利视频在线观看免费| 亚洲专区中文字幕在线| 欧美日韩精品网址| 欧美 亚洲 国产 日韩一| 在线观看66精品国产| 日本黄色视频三级网站网址 | av网站在线播放免费| 精品国产乱码久久久久久小说| 亚洲精品自拍成人| 一个人免费在线观看的高清视频| 人妻一区二区av| 久久精品熟女亚洲av麻豆精品| 欧美亚洲日本最大视频资源| 9191精品国产免费久久| 国产亚洲欧美在线一区二区| 亚洲五月色婷婷综合| 亚洲av日韩精品久久久久久密| 18禁国产床啪视频网站| 五月天丁香电影| 国产成人免费无遮挡视频| 欧美精品av麻豆av| 亚洲精品乱久久久久久| 国产国语露脸激情在线看| 高清毛片免费观看视频网站 | 日本wwww免费看| 免费女性裸体啪啪无遮挡网站| 亚洲中文av在线| 亚洲精品国产精品久久久不卡| 精品人妻在线不人妻| tube8黄色片| 好男人电影高清在线观看| 中文亚洲av片在线观看爽 | 2018国产大陆天天弄谢| av片东京热男人的天堂| 两人在一起打扑克的视频| 国产国语露脸激情在线看| 别揉我奶头~嗯~啊~动态视频| 桃红色精品国产亚洲av| 日韩大码丰满熟妇| 夜夜爽天天搞| 欧美精品一区二区大全| 欧美精品一区二区大全| 久久人妻熟女aⅴ| 国内毛片毛片毛片毛片毛片| 男女床上黄色一级片免费看| 久久中文看片网| 欧美激情 高清一区二区三区| 日韩 欧美 亚洲 中文字幕| 女人高潮潮喷娇喘18禁视频| 成人国产一区最新在线观看| 悠悠久久av| 亚洲专区中文字幕在线| 大型黄色视频在线免费观看| 欧美+亚洲+日韩+国产| 天天躁夜夜躁狠狠躁躁| 午夜免费成人在线视频| av线在线观看网站| 天天躁夜夜躁狠狠躁躁| 免费在线观看影片大全网站| 悠悠久久av| 不卡av一区二区三区| av国产精品久久久久影院| 黄色a级毛片大全视频| 69精品国产乱码久久久| 国产在线视频一区二区| 色老头精品视频在线观看| 熟女少妇亚洲综合色aaa.| 国产精品 国内视频| 99国产极品粉嫩在线观看| 日本av手机在线免费观看| 久久久久久久久久久久大奶| 一区福利在线观看| netflix在线观看网站| 欧美乱妇无乱码| 色在线成人网| 久久久精品94久久精品| 久久久久久久久久久久大奶| 天天躁狠狠躁夜夜躁狠狠躁| 91精品国产国语对白视频| 精品熟女少妇八av免费久了| 精品乱码久久久久久99久播| 国产日韩欧美视频二区| √禁漫天堂资源中文www| av又黄又爽大尺度在线免费看| 久久香蕉激情| 亚洲精品国产一区二区精华液| 最新美女视频免费是黄的| 正在播放国产对白刺激| 自拍欧美九色日韩亚洲蝌蚪91| 成人国产av品久久久| 免费不卡黄色视频| 最黄视频免费看| 别揉我奶头~嗯~啊~动态视频| 美女主播在线视频| 女警被强在线播放| 法律面前人人平等表现在哪些方面| 精品少妇久久久久久888优播| 日韩中文字幕欧美一区二区| 亚洲三区欧美一区| 黄色视频在线播放观看不卡| 黄色成人免费大全| 亚洲三区欧美一区| 不卡一级毛片| 三上悠亚av全集在线观看| 好男人电影高清在线观看| 亚洲少妇的诱惑av| 国产高清视频在线播放一区| 亚洲国产欧美日韩在线播放| 欧美日韩黄片免| 18禁观看日本| 新久久久久国产一级毛片| 国产深夜福利视频在线观看| 国产精品影院久久| 欧美 日韩 精品 国产| 婷婷成人精品国产| 18禁观看日本| 涩涩av久久男人的天堂| 午夜福利视频在线观看免费| 久热这里只有精品99| 国产av一区二区精品久久| avwww免费| 亚洲第一青青草原| 国产在线免费精品| 久久精品国产亚洲av高清一级| 一级片免费观看大全| 欧美日韩视频精品一区| 日韩有码中文字幕| 日韩大片免费观看网站| 亚洲国产中文字幕在线视频| 精品国产乱子伦一区二区三区| 色综合婷婷激情| 老熟妇乱子伦视频在线观看| 精品国产亚洲在线| 韩国精品一区二区三区| 午夜福利欧美成人| a级片在线免费高清观看视频| 中文字幕人妻丝袜一区二区| 亚洲精品国产精品久久久不卡| 亚洲国产精品一区二区三区在线| 在线观看免费日韩欧美大片| 亚洲国产av新网站| 久久久久国产一级毛片高清牌| 国产成人啪精品午夜网站| 国产精品电影一区二区三区 | av福利片在线| 亚洲人成伊人成综合网2020| 成人三级做爰电影| 99热网站在线观看| 99热国产这里只有精品6| 精品一区二区三区四区五区乱码| 精品欧美一区二区三区在线| 久久国产精品人妻蜜桃| 亚洲欧美一区二区三区久久| 曰老女人黄片| 手机成人av网站| 精品国产亚洲在线| 丝袜喷水一区| 久热这里只有精品99| 国产欧美日韩精品亚洲av| 国产一区二区三区在线臀色熟女 | 成人国产av品久久久| 亚洲伊人色综图| 久久中文字幕一级| 亚洲人成电影免费在线| 大码成人一级视频| 黄色片一级片一级黄色片| 成人三级做爰电影| 啦啦啦中文免费视频观看日本| 成人精品一区二区免费| 五月天丁香电影| 妹子高潮喷水视频| 丝瓜视频免费看黄片| avwww免费| 黄色视频不卡| 久久久久久久精品吃奶| 波多野结衣一区麻豆| 久久久久久亚洲精品国产蜜桃av| 一级片'在线观看视频| 精品第一国产精品| 9191精品国产免费久久| 国产精品久久久久久精品电影小说| 婷婷成人精品国产| tube8黄色片| 大陆偷拍与自拍| 大香蕉久久成人网| 精品午夜福利视频在线观看一区 | 久久人人97超碰香蕉20202| 亚洲av日韩精品久久久久久密| 日韩中文字幕视频在线看片| 美女国产高潮福利片在线看| 1024香蕉在线观看| 最近最新免费中文字幕在线| 欧美日韩一级在线毛片| 高清av免费在线| 男男h啪啪无遮挡| 国产精品亚洲一级av第二区| 女人高潮潮喷娇喘18禁视频| 欧美黄色片欧美黄色片| 黄频高清免费视频| 青草久久国产| 欧美中文综合在线视频| 色综合婷婷激情| 日韩一区二区三区影片| 国产在线观看jvid| 亚洲精品粉嫩美女一区| 纵有疾风起免费观看全集完整版| 成人黄色视频免费在线看| 人妻久久中文字幕网| 午夜日韩欧美国产| 建设人人有责人人尽责人人享有的| 新久久久久国产一级毛片| 男女床上黄色一级片免费看| 欧美变态另类bdsm刘玥| 日韩成人在线观看一区二区三区| 好男人电影高清在线观看| 精品一区二区三区视频在线观看免费 | 精品人妻在线不人妻| 亚洲伊人色综图| 亚洲 国产 在线| 久久久久久久国产电影| 老司机靠b影院| 妹子高潮喷水视频| 日韩一区二区三区影片| 亚洲国产毛片av蜜桃av| 岛国毛片在线播放| 成人国语在线视频| 久久精品亚洲精品国产色婷小说| 日韩免费av在线播放| 久久久欧美国产精品| 久久天躁狠狠躁夜夜2o2o| 少妇的丰满在线观看| 丰满少妇做爰视频| 亚洲全国av大片| 久久久水蜜桃国产精品网| 免费人妻精品一区二区三区视频| 久久午夜综合久久蜜桃| 久久国产精品影院| 在线永久观看黄色视频| 51午夜福利影视在线观看| 国产成人免费观看mmmm| 王馨瑶露胸无遮挡在线观看| 久久久久久久大尺度免费视频| 久久天堂一区二区三区四区| 欧美精品人与动牲交sv欧美| 天堂8中文在线网| 男女午夜视频在线观看| 亚洲五月婷婷丁香| 国产精品一区二区免费欧美| 亚洲中文字幕日韩| 91麻豆精品激情在线观看国产 | 国产精品电影一区二区三区 | 丝袜美足系列| 97人妻天天添夜夜摸| 香蕉久久夜色| 丝袜美腿诱惑在线| 国产成人精品久久二区二区免费| 91九色精品人成在线观看| 久久精品亚洲熟妇少妇任你| 男女高潮啪啪啪动态图| 久久99热这里只频精品6学生| 久久婷婷成人综合色麻豆| 欧美日本中文国产一区发布| 69精品国产乱码久久久| 窝窝影院91人妻| 午夜视频精品福利| 十八禁网站网址无遮挡| 男人操女人黄网站| 国产成人av激情在线播放| 国产无遮挡羞羞视频在线观看| 另类亚洲欧美激情| 免费人妻精品一区二区三区视频| 亚洲情色 制服丝袜| 亚洲中文字幕日韩| 亚洲熟妇熟女久久| 女警被强在线播放| 成人精品一区二区免费| 午夜视频精品福利| 深夜精品福利| 久久精品亚洲av国产电影网| 最近最新中文字幕大全电影3 | 色播在线永久视频| 久久久久久久国产电影| av有码第一页| 中亚洲国语对白在线视频| 欧美人与性动交α欧美精品济南到| 操美女的视频在线观看| 国产aⅴ精品一区二区三区波| 在线观看免费视频日本深夜| av视频免费观看在线观看| 1024视频免费在线观看| 免费久久久久久久精品成人欧美视频| 女性被躁到高潮视频| 美女午夜性视频免费| 亚洲熟女精品中文字幕| 国产成人免费无遮挡视频| 国产精品麻豆人妻色哟哟久久| 国产成人精品久久二区二区91| 18禁裸乳无遮挡动漫免费视频| 丰满少妇做爰视频| 国产一区有黄有色的免费视频| 香蕉丝袜av| 国产成人啪精品午夜网站| 老熟妇乱子伦视频在线观看| 丁香欧美五月| 欧美精品人与动牲交sv欧美| 最新美女视频免费是黄的| 欧美亚洲日本最大视频资源| 亚洲熟妇熟女久久| 热99re8久久精品国产| 亚洲av片天天在线观看| 国产在线观看jvid| 久久99一区二区三区| av一本久久久久| 麻豆成人av在线观看| 欧美老熟妇乱子伦牲交| 丝袜人妻中文字幕| 国产日韩一区二区三区精品不卡| 动漫黄色视频在线观看| 黄片小视频在线播放| 91国产中文字幕| av国产精品久久久久影院| 久久久久久人人人人人| 99国产综合亚洲精品| 99热网站在线观看| 99国产精品一区二区三区| 在线观看免费视频网站a站| 国产成人影院久久av| 国产aⅴ精品一区二区三区波| 婷婷成人精品国产| 亚洲av日韩精品久久久久久密| 亚洲伊人久久精品综合| 女人爽到高潮嗷嗷叫在线视频| 久久99热这里只频精品6学生| 一区二区三区乱码不卡18| av天堂在线播放| 黑丝袜美女国产一区| 大码成人一级视频| 久久国产精品影院| 99精品欧美一区二区三区四区| 国产欧美日韩综合在线一区二区| 丝袜美足系列| av欧美777| 久久性视频一级片| 99国产精品99久久久久| 在线观看人妻少妇| 欧美在线一区亚洲| 亚洲欧洲精品一区二区精品久久久| 久久久精品区二区三区| 国产一区二区激情短视频| 免费不卡黄色视频| 欧美日韩精品网址| 深夜精品福利| 国产97色在线日韩免费| 热re99久久精品国产66热6| 午夜福利免费观看在线| 精品午夜福利视频在线观看一区 | 免费黄频网站在线观看国产| 中国美女看黄片| 色播在线永久视频| av视频免费观看在线观看| 精品熟女少妇八av免费久了| 超碰成人久久| 国产精品二区激情视频| 久久久久久亚洲精品国产蜜桃av| 在线观看66精品国产| 午夜福利免费观看在线| 一级毛片电影观看| 乱人伦中国视频| 日韩人妻精品一区2区三区| 亚洲色图综合在线观看| 在线观看一区二区三区激情| 久久九九热精品免费| netflix在线观看网站| 在线观看免费视频网站a站| 最黄视频免费看| 国产成人免费无遮挡视频| 欧美日韩黄片免| 男女之事视频高清在线观看| 我的亚洲天堂| 一进一出好大好爽视频| 成人三级做爰电影| 欧美精品高潮呻吟av久久| 在线av久久热| 成人特级黄色片久久久久久久 | 一边摸一边抽搐一进一小说 | 精品第一国产精品| 免费黄频网站在线观看国产| 欧美乱妇无乱码| 在线十欧美十亚洲十日本专区| 99在线人妻在线中文字幕 | 狂野欧美激情性xxxx| 每晚都被弄得嗷嗷叫到高潮| 亚洲av第一区精品v没综合| 国产av国产精品国产| 日韩欧美一区二区三区在线观看 | av有码第一页| 国产精品一区二区免费欧美| 天天躁狠狠躁夜夜躁狠狠躁| 男人操女人黄网站| 国产精品久久久久久精品电影小说| 国产精品.久久久| 中文字幕最新亚洲高清| 精品人妻1区二区| 黄片大片在线免费观看| 老司机深夜福利视频在线观看| 99re6热这里在线精品视频| 亚洲免费av在线视频| 妹子高潮喷水视频| 欧美在线一区亚洲| 一区二区三区国产精品乱码| 国产精品国产av在线观看| 国产国语露脸激情在线看| 亚洲欧美日韩另类电影网站| 免费高清在线观看日韩| 久久午夜综合久久蜜桃| 搡老岳熟女国产| 国产欧美日韩综合在线一区二区| 激情在线观看视频在线高清 | 欧美精品高潮呻吟av久久| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美亚洲日本最大视频资源| 久久免费观看电影| 香蕉久久夜色| 久久天堂一区二区三区四区| 三级毛片av免费| 国产精品免费视频内射| 久久国产精品大桥未久av| 久久午夜综合久久蜜桃| 国产精品久久久av美女十八| 美女主播在线视频| 一级a爱视频在线免费观看| 黄色视频不卡| 中文字幕最新亚洲高清| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 欧美精品亚洲一区二区| 俄罗斯特黄特色一大片| 淫妇啪啪啪对白视频| av欧美777| 首页视频小说图片口味搜索| 日韩免费av在线播放| 嫁个100分男人电影在线观看| 九色亚洲精品在线播放| 日韩欧美免费精品| 成人精品一区二区免费| 人人妻人人爽人人添夜夜欢视频| 欧美中文综合在线视频| 一个人免费看片子| 国产国语露脸激情在线看| 久久久久久免费高清国产稀缺| 久久久久久久久久久久大奶| 久久精品亚洲精品国产色婷小说| 女人久久www免费人成看片| 天堂8中文在线网| 人妻 亚洲 视频| 国产成人欧美在线观看 | 一级a爱视频在线免费观看| 狠狠婷婷综合久久久久久88av| 亚洲av片天天在线观看| 国产三级黄色录像| 一级毛片女人18水好多| 中文字幕人妻丝袜一区二区| 欧美大码av| 热99久久久久精品小说推荐| 免费久久久久久久精品成人欧美视频| 天天躁狠狠躁夜夜躁狠狠躁| 操美女的视频在线观看| 欧美乱码精品一区二区三区| 亚洲av日韩精品久久久久久密| 50天的宝宝边吃奶边哭怎么回事| 欧美激情 高清一区二区三区| 成人18禁在线播放| 亚洲精品久久成人aⅴ小说| 日韩中文字幕视频在线看片| 丝袜喷水一区| 一本色道久久久久久精品综合| 日韩欧美国产一区二区入口| 免费在线观看黄色视频的| 丝袜在线中文字幕| 亚洲天堂av无毛| 999精品在线视频| 汤姆久久久久久久影院中文字幕| 亚洲专区国产一区二区| 亚洲国产欧美一区二区综合| 国产精品久久久久久精品电影小说| 男女下面插进去视频免费观看| 老司机亚洲免费影院| 亚洲精品中文字幕一二三四区 | 国产在线免费精品| 国精品久久久久久国模美| 在线天堂中文资源库| 脱女人内裤的视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲伊人色综图| 桃花免费在线播放| 免费黄频网站在线观看国产| 18禁裸乳无遮挡动漫免费视频| 亚洲美女黄片视频| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久国产一级毛片高清牌| 巨乳人妻的诱惑在线观看| 高清av免费在线| 久热这里只有精品99| 中文字幕最新亚洲高清| 日韩欧美三级三区| 怎么达到女性高潮|