• <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ì)損失估算
    狠狠狠狠99中文字幕| 欧美变态另类bdsm刘玥| 在线a可以看的网站| 国产免费一级a男人的天堂| 看非洲黑人一级黄片| 亚洲欧洲日产国产| 亚洲国产色片| 国产又黄又爽又无遮挡在线| 欧美最新免费一区二区三区| 99久久精品热视频| 精品久久久久久久末码| 在线播放国产精品三级| 欧美日韩综合久久久久久| 亚洲乱码一区二区免费版| 免费电影在线观看免费观看| 国产在视频线精品| 99久久成人亚洲精品观看| 国产精品国产高清国产av| 国产亚洲5aaaaa淫片| 国产极品精品免费视频能看的| 国产亚洲精品av在线| 高清在线视频一区二区三区 | 国产精品野战在线观看| 亚洲国产高清在线一区二区三| 男人舔女人下体高潮全视频| 国产高清有码在线观看视频| 久久精品久久久久久久性| 欧美变态另类bdsm刘玥| 欧美日本亚洲视频在线播放| 青春草亚洲视频在线观看| 亚洲国产色片| kizo精华| 国产 一区精品| 亚洲综合精品二区| 久久久久网色| 国产淫语在线视频| 国语自产精品视频在线第100页| 一区二区三区高清视频在线| 久久久久久久午夜电影| 韩国高清视频一区二区三区| 少妇熟女aⅴ在线视频| 免费看a级黄色片| 日韩强制内射视频| 婷婷色av中文字幕| 一个人看视频在线观看www免费| 欧美精品一区二区大全| 成年女人永久免费观看视频| 少妇熟女欧美另类| 国产亚洲精品av在线| 日本五十路高清| 听说在线观看完整版免费高清| 三级毛片av免费| 国产伦精品一区二区三区四那| 人妻夜夜爽99麻豆av| 精品午夜福利在线看| 国产亚洲av嫩草精品影院| 午夜久久久久精精品| 精品国产一区二区三区久久久樱花 | 欧美一区二区国产精品久久精品| 麻豆成人午夜福利视频| 日韩欧美在线乱码| 免费看a级黄色片| 18禁在线播放成人免费| 午夜福利成人在线免费观看| 九九在线视频观看精品| 一边亲一边摸免费视频| 青春草国产在线视频| 久久久欧美国产精品| 国产乱人视频| 观看美女的网站| 免费看日本二区| 国产成人freesex在线| av免费观看日本| 一区二区三区四区激情视频| 欧美极品一区二区三区四区| 中文亚洲av片在线观看爽| av免费在线看不卡| 精品欧美国产一区二区三| 国产色婷婷99| 蜜臀久久99精品久久宅男| 国产精品精品国产色婷婷| 国产亚洲一区二区精品| 九草在线视频观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 精品国产三级普通话版| 日韩欧美精品免费久久| 国产精品野战在线观看| 少妇被粗大猛烈的视频| 2021少妇久久久久久久久久久| 一级爰片在线观看| 国产精品综合久久久久久久免费| 天美传媒精品一区二区| 一卡2卡三卡四卡精品乱码亚洲| 大香蕉久久网| 久久久久久伊人网av| 2021少妇久久久久久久久久久| 大又大粗又爽又黄少妇毛片口| 成人毛片a级毛片在线播放| 国产一区二区在线av高清观看| 欧美xxxx性猛交bbbb| 免费av观看视频| 中文欧美无线码| 狠狠狠狠99中文字幕| 亚洲图色成人| 欧美一区二区精品小视频在线| 九九热线精品视视频播放| a级毛色黄片| 伦精品一区二区三区| 欧美日韩精品成人综合77777| 色播亚洲综合网| 国产在线一区二区三区精 | 亚洲av熟女| 日本三级黄在线观看| 欧美区成人在线视频| 亚洲精品色激情综合| 亚洲国产精品成人久久小说| 波多野结衣巨乳人妻| 婷婷色麻豆天堂久久 | 亚洲欧美日韩高清专用| 欧美激情在线99| 最近的中文字幕免费完整| 高清视频免费观看一区二区 | 狠狠狠狠99中文字幕| 亚洲国产精品国产精品| 中文字幕av成人在线电影| 我要看日韩黄色一级片| 一卡2卡三卡四卡精品乱码亚洲| 在线a可以看的网站| 久久人人爽人人爽人人片va| 国产精品久久久久久久久免| 美女脱内裤让男人舔精品视频| 精品午夜福利在线看| 热99re8久久精品国产| 亚洲性久久影院| 在线观看一区二区三区| 午夜爱爱视频在线播放| h日本视频在线播放| 国产精品,欧美在线| 久久久久网色| 99久久九九国产精品国产免费| 七月丁香在线播放| 欧美丝袜亚洲另类| 99视频精品全部免费 在线| 变态另类丝袜制服| 成人亚洲精品av一区二区| 插逼视频在线观看| 中国国产av一级| 久热久热在线精品观看| 色网站视频免费| 中文精品一卡2卡3卡4更新| 三级毛片av免费| 成人高潮视频无遮挡免费网站| 国产私拍福利视频在线观看| 久久久久性生活片| 亚洲国产高清在线一区二区三| 国产精品日韩av在线免费观看| 亚洲一区高清亚洲精品| 欧美区成人在线视频| 男女边吃奶边做爰视频| 老司机福利观看| 九九在线视频观看精品| 偷拍熟女少妇极品色| 国产av一区在线观看免费| 国产高潮美女av| 欧美一区二区亚洲| 国产亚洲av片在线观看秒播厂 | 成人亚洲精品av一区二区| 国产老妇女一区| 男女边吃奶边做爰视频| 亚洲国产精品国产精品| 汤姆久久久久久久影院中文字幕 | 毛片女人毛片| 日日摸夜夜添夜夜添av毛片| 日韩欧美精品v在线| 高清av免费在线| 超碰97精品在线观看| 麻豆一二三区av精品| 精品人妻偷拍中文字幕| 69人妻影院| av在线老鸭窝| 亚洲自拍偷在线| 成人无遮挡网站| 天堂av国产一区二区熟女人妻| 深夜a级毛片| 国模一区二区三区四区视频| 午夜视频国产福利| 在线免费观看的www视频| 99热6这里只有精品| 精品酒店卫生间| 高清av免费在线| 嘟嘟电影网在线观看| 午夜福利在线观看免费完整高清在| 国产精品爽爽va在线观看网站| 国产91av在线免费观看| 日本五十路高清| 一边摸一边抽搐一进一小说| 欧美+日韩+精品| av免费在线看不卡| 女的被弄到高潮叫床怎么办| 日韩精品青青久久久久久| 欧美zozozo另类| 中文字幕久久专区| 草草在线视频免费看| 国产人妻一区二区三区在| 久久鲁丝午夜福利片| 黑人高潮一二区| 久久精品久久久久久久性| 麻豆成人av视频| 午夜福利在线观看免费完整高清在| 日韩大片免费观看网站 | 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 极品教师在线视频| 婷婷色麻豆天堂久久 | 男女啪啪激烈高潮av片| 久久久久九九精品影院| 天天一区二区日本电影三级| 国产亚洲午夜精品一区二区久久 | 久久人人爽人人片av| 成人亚洲欧美一区二区av| 99在线人妻在线中文字幕| 久久久久九九精品影院| 一卡2卡三卡四卡精品乱码亚洲| 伦理电影大哥的女人| 欧美又色又爽又黄视频| 久久久久精品久久久久真实原创| 丝袜美腿在线中文| 久久99热6这里只有精品| 欧美一级a爱片免费观看看| 你懂的网址亚洲精品在线观看 | 国产91av在线免费观看| www.av在线官网国产| 欧美一区二区精品小视频在线| 亚洲精品日韩av片在线观看| 欧美性猛交黑人性爽| 久久久久九九精品影院| 观看免费一级毛片| 夫妻性生交免费视频一级片| 久久久成人免费电影| 久99久视频精品免费| 一个人看的www免费观看视频| 国产精品日韩av在线免费观看| 色尼玛亚洲综合影院| 自拍偷自拍亚洲精品老妇| 日韩高清综合在线| 国产精品一区二区三区四区免费观看| 91久久精品国产一区二区三区| videos熟女内射| 国产一区二区三区av在线| 在现免费观看毛片| av在线亚洲专区| 精品一区二区三区人妻视频| 色吧在线观看| 观看美女的网站| 最近2019中文字幕mv第一页| 亚洲丝袜综合中文字幕| 欧美日韩综合久久久久久| 久久精品国产自在天天线| 噜噜噜噜噜久久久久久91| 国产精品美女特级片免费视频播放器| 蜜桃久久精品国产亚洲av| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 欧美一区二区亚洲| 18禁动态无遮挡网站| 不卡视频在线观看欧美| 1000部很黄的大片| 国产视频首页在线观看| 联通29元200g的流量卡| 亚洲精品乱久久久久久| videos熟女内射| www日本黄色视频网| 精品人妻熟女av久视频| av在线蜜桃| 性插视频无遮挡在线免费观看| 小蜜桃在线观看免费完整版高清| 国产男人的电影天堂91| 国产成人a区在线观看| www日本黄色视频网| 亚洲四区av| 在线播放无遮挡| 村上凉子中文字幕在线| 国产一区二区亚洲精品在线观看| 韩国av在线不卡| 国产精品不卡视频一区二区| av免费观看日本| 少妇猛男粗大的猛烈进出视频 | 青青草视频在线视频观看| 99久久人妻综合| 男人和女人高潮做爰伦理| 亚洲欧美精品综合久久99| 小说图片视频综合网站| 97在线视频观看| 久久综合国产亚洲精品| 欧美极品一区二区三区四区| 中文亚洲av片在线观看爽| a级毛片免费高清观看在线播放| 韩国高清视频一区二区三区| 大香蕉久久网| 久久精品影院6| 一卡2卡三卡四卡精品乱码亚洲| 国产极品精品免费视频能看的| 成人特级av手机在线观看| 三级国产精品片| h日本视频在线播放| 色吧在线观看| 春色校园在线视频观看| 中文资源天堂在线| 国产成年人精品一区二区| 人妻系列 视频| 麻豆乱淫一区二区| 色视频www国产| 国产综合懂色| 国产高清不卡午夜福利| 99久久中文字幕三级久久日本| 一区二区三区免费毛片| 日韩欧美国产在线观看| 亚洲成人av在线免费| 久久人人爽人人爽人人片va| 国产片特级美女逼逼视频| 免费一级毛片在线播放高清视频| 亚洲精品亚洲一区二区| 国产精品女同一区二区软件| a级毛色黄片| 久久久久国产网址| 99热网站在线观看| 午夜爱爱视频在线播放| 麻豆成人av视频| 欧美一区二区国产精品久久精品| av专区在线播放| 亚洲精品色激情综合| 国产av一区在线观看免费| 精品国产三级普通话版| 久久久久久久久久久丰满| 久久久久精品久久久久真实原创| 国产91av在线免费观看| 久久这里只有精品中国| 国产乱人偷精品视频| 久久久久精品久久久久真实原创| 高清毛片免费看| 色噜噜av男人的天堂激情| 国产v大片淫在线免费观看| 亚洲成人精品中文字幕电影| 久久久久性生活片| 亚洲自拍偷在线| 久久99热6这里只有精品| 日韩精品青青久久久久久| 亚洲欧美一区二区三区国产| 久久热精品热| 麻豆乱淫一区二区| 欧美日韩在线观看h| 青春草国产在线视频| 久久久午夜欧美精品| 狂野欧美白嫩少妇大欣赏| 在线观看av片永久免费下载| 国产中年淑女户外野战色| 国产精品综合久久久久久久免费| 亚洲人与动物交配视频| 午夜福利在线在线| 久久久国产成人精品二区| 中文欧美无线码| 欧美3d第一页| 五月伊人婷婷丁香| 成年av动漫网址| 亚洲欧美成人综合另类久久久 | 国产 一区精品| 国产单亲对白刺激| 最新中文字幕久久久久| 午夜福利成人在线免费观看| 色噜噜av男人的天堂激情| 精品久久久久久成人av| 亚洲av.av天堂| 久久人妻av系列| 亚洲av免费在线观看| 国产高清三级在线| 国产国拍精品亚洲av在线观看| 丝袜喷水一区| 日本猛色少妇xxxxx猛交久久| 一级毛片aaaaaa免费看小| 少妇的逼好多水| 美女被艹到高潮喷水动态| 国产黄a三级三级三级人| 可以在线观看毛片的网站| 九色成人免费人妻av| 大又大粗又爽又黄少妇毛片口| 高清av免费在线| 国产伦精品一区二区三区四那| 久久精品综合一区二区三区| 欧美日本视频| 日日干狠狠操夜夜爽| 大香蕉久久网| 亚洲高清免费不卡视频| 国产午夜福利久久久久久| 日本欧美国产在线视频| 国产乱人视频| 免费看av在线观看网站| 久久精品久久久久久噜噜老黄 | 日本猛色少妇xxxxx猛交久久| 免费一级毛片在线播放高清视频| av在线天堂中文字幕| 内射极品少妇av片p| av在线播放精品| 欧美不卡视频在线免费观看| 成人无遮挡网站| 成年女人看的毛片在线观看| 亚洲国产色片| 亚洲精品亚洲一区二区| 亚洲精品456在线播放app| eeuss影院久久| 蜜桃久久精品国产亚洲av| 99热网站在线观看| 联通29元200g的流量卡| 欧美性猛交╳xxx乱大交人| 国产探花极品一区二区| 欧美成人a在线观看| 高清午夜精品一区二区三区| 中文乱码字字幕精品一区二区三区 | 人妻夜夜爽99麻豆av| 在线免费十八禁| 亚洲av中文av极速乱| 超碰97精品在线观看| 男人舔女人下体高潮全视频| 十八禁国产超污无遮挡网站| 国产成人精品婷婷| 天堂av国产一区二区熟女人妻| 乱人视频在线观看| 伦理电影大哥的女人| 日韩高清综合在线| 午夜福利高清视频| 国产69精品久久久久777片| 少妇人妻一区二区三区视频| 亚洲国产精品国产精品| 高清av免费在线| kizo精华| 卡戴珊不雅视频在线播放| 国产熟女欧美一区二区| 亚洲成人久久爱视频| 久久99蜜桃精品久久| 久久久久久久国产电影| 亚洲精品色激情综合| 精品久久久久久久人妻蜜臀av| 成人漫画全彩无遮挡| 久久久久国产网址| 美女cb高潮喷水在线观看| 熟女电影av网| 色播亚洲综合网| 两性午夜刺激爽爽歪歪视频在线观看| 有码 亚洲区| 国产一区有黄有色的免费视频 | 国产极品精品免费视频能看的| 美女内射精品一级片tv| 精品少妇黑人巨大在线播放 | 在线播放无遮挡| 精品久久国产蜜桃| 久久久a久久爽久久v久久| 只有这里有精品99| 中国国产av一级| 亚洲不卡免费看| 久久99蜜桃精品久久| 国产淫片久久久久久久久| 亚洲第一区二区三区不卡| 18禁在线无遮挡免费观看视频| 伊人久久精品亚洲午夜| 一级二级三级毛片免费看| 国产欧美日韩精品一区二区| 在线免费观看的www视频| 国产乱人偷精品视频| 日韩精品有码人妻一区| 国产单亲对白刺激| 2021天堂中文幕一二区在线观| 最近的中文字幕免费完整| 22中文网久久字幕| 成人性生交大片免费视频hd| 永久网站在线| 一级毛片我不卡| a级毛片免费高清观看在线播放| 国产精品永久免费网站| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类| 婷婷色麻豆天堂久久 | 亚洲av二区三区四区| 91精品伊人久久大香线蕉| 欧美激情久久久久久爽电影| 国产av不卡久久| 久久综合国产亚洲精品| 国产欧美日韩精品一区二区| 亚洲自偷自拍三级| 亚洲天堂国产精品一区在线| 国产不卡一卡二| 韩国高清视频一区二区三区| 国产精品福利在线免费观看| 人妻少妇偷人精品九色| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 三级男女做爰猛烈吃奶摸视频| 日日干狠狠操夜夜爽| 久久久欧美国产精品| 国产视频首页在线观看| 国产又色又爽无遮挡免| 狂野欧美激情性xxxx在线观看| 麻豆一二三区av精品| 国产亚洲午夜精品一区二区久久 | 国国产精品蜜臀av免费| 久久精品熟女亚洲av麻豆精品 | 日韩欧美精品v在线| 久久久国产成人精品二区| 欧美日韩在线观看h| 国产单亲对白刺激| 97人妻精品一区二区三区麻豆| 啦啦啦啦在线视频资源| 黄片无遮挡物在线观看| 免费看日本二区| 看非洲黑人一级黄片| 成人亚洲精品av一区二区| 看非洲黑人一级黄片| 欧美成人免费av一区二区三区| 免费av毛片视频| 亚洲av中文字字幕乱码综合| 国产成人a∨麻豆精品| 一区二区三区高清视频在线| 久久久久久久国产电影| 亚洲国产欧美人成| 色综合色国产| 亚洲欧美中文字幕日韩二区| eeuss影院久久| 大又大粗又爽又黄少妇毛片口| 精品免费久久久久久久清纯| 麻豆久久精品国产亚洲av| 极品教师在线视频| 天堂√8在线中文| 亚洲欧美日韩东京热| 久久国内精品自在自线图片| 亚洲欧美日韩无卡精品| 国产色婷婷99| 69av精品久久久久久| 男的添女的下面高潮视频| 中文字幕亚洲精品专区| 国产精品三级大全| 三级经典国产精品| 非洲黑人性xxxx精品又粗又长| 热99在线观看视频| 国产午夜精品一二区理论片| 丰满少妇做爰视频| 村上凉子中文字幕在线| 亚洲在线观看片| 我要搜黄色片| 日本黄大片高清| 91精品国产九色| 丰满人妻一区二区三区视频av| 欧美+日韩+精品| 2021少妇久久久久久久久久久| 亚洲精品日韩av片在线观看| 国产探花在线观看一区二区| 色播亚洲综合网| 亚洲精品乱码久久久久久按摩| 久久精品国产亚洲av天美| 狂野欧美激情性xxxx在线观看| 最近中文字幕2019免费版| 青春草国产在线视频| 亚洲精品乱久久久久久| 日韩,欧美,国产一区二区三区 | 91av网一区二区| 国产精品麻豆人妻色哟哟久久 | 18禁动态无遮挡网站| 麻豆国产97在线/欧美| h日本视频在线播放| 免费在线观看成人毛片| 免费电影在线观看免费观看| 久久精品国产亚洲av涩爱| 男人和女人高潮做爰伦理| 免费观看a级毛片全部| 联通29元200g的流量卡| 欧美成人a在线观看| 蜜臀久久99精品久久宅男| 日本一本二区三区精品| 一区二区三区四区激情视频| 日韩av不卡免费在线播放| 国产在视频线在精品| 偷拍熟女少妇极品色| 天堂中文最新版在线下载 | 色播亚洲综合网| 欧美又色又爽又黄视频| 亚洲人成网站高清观看| 免费观看性生交大片5| 国产乱人视频| 国产亚洲av片在线观看秒播厂 | 午夜福利在线观看吧| 久久99蜜桃精品久久| 亚洲三级黄色毛片| 亚洲高清免费不卡视频| 亚洲自拍偷在线| 美女被艹到高潮喷水动态| av国产久精品久网站免费入址| 欧美日韩在线观看h| 可以在线观看毛片的网站| 国模一区二区三区四区视频| 深夜a级毛片| 天堂av国产一区二区熟女人妻| 最后的刺客免费高清国语| 少妇人妻一区二区三区视频| 免费一级毛片在线播放高清视频| 18禁裸乳无遮挡免费网站照片| 精品人妻视频免费看| a级一级毛片免费在线观看| 97热精品久久久久久| 哪个播放器可以免费观看大片| 精品久久久久久久人妻蜜臀av| 亚洲欧美清纯卡通| 男人舔女人下体高潮全视频| 国产精品不卡视频一区二区| 久久久久网色| 精品一区二区免费观看| 欧美人与善性xxx| 亚洲欧美一区二区三区国产| 亚洲精品乱久久久久久|