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

    基于燈光遙感圖像的地震災(zāi)情信息提取

    2018-11-05 10:56:26張小詠朱俊杰許建華
    地震研究 2018年2期
    關(guān)鍵詞:差值燈光災(zāi)害

    張小詠 朱俊杰 許建華

    摘要:燈光遙感圖像變化是強(qiáng)震后燈光變化的直接表現(xiàn)?;趯Φ卣鹎昂鬅艄膺b感圖像特點(diǎn)分析,建立了一種基于燈光破壞的地震災(zāi)害信息提取方法,并在2001年印度7.9級地震、2011年“3·11”日本9.0級地震進(jìn)行了應(yīng)用分析。結(jié)果表明:基于顯著性檢驗(yàn)的燈光變化檢測方法簡潔有效,利用燈光遙感數(shù)據(jù)不僅能夠從宏觀上獲得地震破壞范圍及其程度,而且能夠直接反映震后受災(zāi)人口規(guī)模及其空間分布。

    關(guān)鍵詞:燈光遙感;地震災(zāi)害顯著性檢驗(yàn);DMSP/OLS

    中圖分類號(hào):P315.94 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1000-0666(2018)02-0311-08

    0 引言

    災(zāi)害規(guī)模、受災(zāi)程度,以及災(zāi)害的空間分布等信息是災(zāi)后應(yīng)急指揮和決策的依據(jù),特別是在發(fā)生嚴(yán)重地震災(zāi)害時(shí),通訊與電力中斷往往導(dǎo)致災(zāi)區(qū)與外界失去聯(lián)系,如何快速獲取災(zāi)情信息進(jìn)行緊急救援成為指揮決策的迫切需求。長期以來,地震工作者一直探索著強(qiáng)震巨災(zāi)發(fā)生后如何快速有效地獲取災(zāi)情信息。當(dāng)前,國內(nèi)外利用遙感技術(shù)獲取地震災(zāi)害信息主要集中在對震后建筑物、橋梁、斷層等地表形變破壞以及崩塌、滑坡等地震誘發(fā)次生災(zāi)害的研究(張景發(fā)等,2002;單新建等,2005;郭華東等,2011;王曉青等,2015;張小詠等,2016),但地震所帶來的破壞還包括水、電、氣、通訊、交通等生命線工程的損失(陸鳴等,2008;劉愛文等,2012;李東平等,2017)。一旦發(fā)生強(qiáng)震,災(zāi)區(qū)的電力系統(tǒng)可能受損,這種破壞不僅能反映地震造成的直接破壞,如房屋倒塌、道路中斷等明顯破壞,還能反映表面沒有破壞而實(shí)際造成的經(jīng)濟(jì)損失,如電力中斷導(dǎo)致的生產(chǎn)停止,這種軟破壞對經(jīng)濟(jì)和民生也是一種直接的損失。而燈光遙感圖像可以直接反映這類信息,在地震破壞范圍確定、受災(zāi)人口、經(jīng)濟(jì)損失評估等有著重要的應(yīng)用價(jià)值。

    1976年,美國軍事氣象衛(wèi)星(Defense Meteoro-logical Satellite Program,簡稱DMSP)開始搭載的Operational Lines can System(簡稱OLS)傳感器可以在夜間獲取燈光影像,星下點(diǎn)空間分辨率1km,具有每天覆蓋全球的觀測能力,幅寬達(dá)2960km,能夠探測到城市燈光甚至小規(guī)模居民地、車流等發(fā)出的低強(qiáng)度燈光,并使之明顯區(qū)別于黑暗背景(陳晉等,2003;何春陽等,2006)。由于燈光直接與人類活動(dòng)密切相關(guān),DMSP/OLS已經(jīng)廣泛應(yīng)用于多個(gè)方面,包括燈光與經(jīng)濟(jì)活動(dòng)之間的相關(guān)性(Elvidgeet al,1997),使用燈光數(shù)據(jù)定量估算美國、日本、中國,甚至全球人口(Sutton,1997;Kimberly,2010):燈光與地震災(zāi)害之間的相關(guān)性,如日本和美國合作研制的地震災(zāi)害評估系統(tǒng)(Kohiyama et al,2004),且該評估系統(tǒng)已被應(yīng)用到日本阪神地震評估中(Hayashi et al,2000)。

    由于燈光遙感圖像在地震災(zāi)害評估方面具有重要的應(yīng)用價(jià)值,本文利用DMSP/OLS數(shù)據(jù)對2001年印度7.9級地震和2011年“3·11”日本9.0級大地震后燈光遙感圖像破壞信息進(jìn)行提取,并討論了燈光遙感數(shù)據(jù)用于強(qiáng)震巨災(zāi)后災(zāi)情信息獲取的應(yīng)用潛力和局限性。

    1 研究區(qū)與數(shù)據(jù)選取

    北京時(shí)間2001年1月26日11時(shí)16分,印度古吉拉特邦發(fā)生7.9級地震,此次地震造成大量建筑物倒塌,2.5萬人死亡,16.7萬人受傷。2011年3月11日13時(shí)46分,日本東北部近海發(fā)生9.0級特大地震,并引發(fā)高達(dá)40.5m海嘯,造成15985人死亡,13346人失蹤。這2次大地震發(fā)生后震區(qū)內(nèi)電力、通信設(shè)備受到嚴(yán)重破壞。本文收集這2次地震前后燈光遙感圖像用于分析,具體數(shù)據(jù)見圖1、2。從圖中可以看出,夜間燈光遙感圖像相對于白天的對地觀測衛(wèi)星圖像非常簡單,僅有燈光和背景2種信息,還可以清楚識(shí)別出城市和較大居民點(diǎn),包括大城市和中小城市。對比印度和日本地震前后燈光圖像(圖1、2),不同國家、不同經(jīng)濟(jì)發(fā)展水平、不同能源管理方式的燈光亮度差別很大,日本經(jīng)濟(jì)發(fā)達(dá)、燈光成片且亮度大,而印度燈光相對稀疏。同樣,日本地震造成的燈光破壞也更加明顯。

    2 基于顯著性檢驗(yàn)的震害信息提取方法

    燈光代表了人類活動(dòng),正常情況下同一地區(qū)的燈光相對穩(wěn)定(圖2a,b,c)。地震發(fā)生后,震區(qū)燈光強(qiáng)度明顯減弱,燈光范圍減小或者消失,發(fā)生了顯著變化(圖2d),地震破壞引起燈光變化是震害信息的一個(gè)明顯表征。基于上述現(xiàn)象分析,可以通過檢測地震前后燈光遙感圖像的變化提取震害信息。目前,基于像素的遙感圖像變化檢測常用的方法有差值法和比值法,考慮到燈光遙感圖像空間分辨率較低,且由地震引起的變化明顯,本文選定差值法,該方法的關(guān)鍵在于如何確定差值的閾值,如果圖像差值滿足一定概率分布,基于數(shù)理統(tǒng)計(jì)的顯著性檢驗(yàn)方法是獲得差值分割閾值的經(jīng)典方法(李春干,代華兵,2017;高波等,2012),也避免了人為給定一個(gè)閾值的隨意性。

    為了分析基于顯著性檢驗(yàn)方法進(jìn)行變化檢測的合理性,分別對2個(gè)研究區(qū)的不同時(shí)相圖像進(jìn)行了輻射歸一化,然后計(jì)算了不同時(shí)相圖像的差值,并根據(jù)差值直方圖擬合得到高斯分布曲線(圖3)。從圖3可發(fā)現(xiàn):(1)地震前不同時(shí)相圖像差值完全滿足正態(tài)分布:(2)圖3a和3b中對稱軸左側(cè)的差值分布基本一樣,而圖3b中對稱軸右側(cè)明顯不對稱:(3)日本和印度地震前后圖像差值規(guī)律一樣(圖3b,c),整體上符合正態(tài)分布,但差值在大于0方向的分布明顯多于小于0方向,且在大于0方向有很多數(shù)值超過高斯分布曲線。

    進(jìn)一步分析可知,正常情況下燈光遙感圖像沒有明顯變化,震前不同時(shí)相燈光遙感圖像的差值主要是隨機(jī)噪聲。從統(tǒng)計(jì)意義上來講,遙感圖像每個(gè)像素都可以看作是一個(gè)統(tǒng)計(jì)樣本,因此震前不同時(shí)相圖像差值的直方圖分布應(yīng)該滿足高斯分布,日本地震前圖像差值的統(tǒng)計(jì)結(jié)果也證明了這一點(diǎn)(圖3a)。當(dāng)?shù)卣鹨馃艄庾兓?,地震前后圖像的差值直方圖雖然整體上仍然符合正態(tài)分布規(guī)律,但是在大于0的方向很多差值分布在高斯分布曲線之外,這部分像素是由于地震引起燈光熄滅或者亮度減弱導(dǎo)致地震前后差值很大,從而出現(xiàn)圖3b和3c中的非對稱現(xiàn)象。口本地震的破壞范圍大,對應(yīng)的差值直方圖在高斯擬合曲線之外的像素多;而印度地震破壞相對小,在高斯分布曲線之外的像素相對較少,但表現(xiàn)的規(guī)律和特征都一樣。因此,可以通過地震前后差值直方圖擬合得到高斯概率密度分布函數(shù),進(jìn)而利用該函數(shù)的顯著性檢驗(yàn)的方法來提取地震導(dǎo)致的燈光破壞區(qū)域。

    根據(jù)顯著性檢驗(yàn),當(dāng)概率分布函數(shù)和置信度確定后,可計(jì)算出概率密度函數(shù)的自變量值。本研究中,自變量是地震前后圖像的差值,置信度對應(yīng)的自變量數(shù)值是差值圖像分割的閾值。置信度的設(shè)置通常采用0.9、0.95和0.99,0.99的置信度一般用于醫(yī)學(xué)等要求非常高的領(lǐng)域,在遙感影像的變化檢測中0.9的置信度結(jié)果存在較大的虛警率(Kohiyama et al,2004;高波等,2012)。本文的燈光遙感圖像分辨率較低,且由于多種原因使得圖像存在較大的噪聲,因此置信度設(shè)置為0.95,這樣檢測的結(jié)果不會(huì)導(dǎo)致虛警率或漏檢率過高。由于變化區(qū)域僅僅體現(xiàn)在差值直方圖大于。的方向,在利用顯著性檢驗(yàn)確定閾值時(shí),選擇右側(cè)檢驗(yàn)。

    利用地震前后燈光遙感圖像進(jìn)行震害信息提取流程如圖4所示。由于獲取的DMSP/OLS燈光數(shù)據(jù)是DN值圖像,沒有進(jìn)行地表輻射量的定量化轉(zhuǎn)換和去云和去噪聲等處理,所以先進(jìn)行輻射歸一化消除成像過程中產(chǎn)生的誤差。具體采用統(tǒng)計(jì)量法,該方法假設(shè)同一地區(qū)的遙感圖像具有相同的概率分布,即假設(shè)圖像具有相同的均值和方差,所以輻射歸一化可以通過源圖像和參考圖像概率分布的統(tǒng)計(jì)量進(jìn)行匹配。根據(jù)這一假設(shè)則有:對式(1)進(jìn)行變換,得到統(tǒng)計(jì)量法的歸一化模型為式中:b1為原圖像歸一化后的像素灰度值;μ1,σ1分別為參考圖像的像素均值和標(biāo)準(zhǔn)差;b2,μ2,σ2分別為原圖像的像素灰度值、均值和標(biāo)準(zhǔn)差。

    輻射歸一化后進(jìn)行地震前后圖像的差值計(jì)算,同時(shí)統(tǒng)計(jì)差值圖像的直方圖,并擬合得到高斯分布函數(shù),高斯函數(shù)采用:

    利用高斯擬合的概率分布函數(shù)進(jìn)行顯著性檢驗(yàn),計(jì)算置信度為0.95條件下的差值圖像的分割閾值,最后根據(jù)閾值對差值圖像進(jìn)行分割提取震害區(qū)域。

    3 計(jì)算結(jié)果及討論

    根據(jù)本文的計(jì)算方法和計(jì)算流程,首先對地震前后燈光遙感圖像進(jìn)行輻射歸一化。在統(tǒng)計(jì)歸一化轉(zhuǎn)換參數(shù)時(shí),為了降低云及其他噪聲的影響,在圖像中選擇海洋上空若干不受云和噪聲影響的區(qū)域作為統(tǒng)計(jì)樣本,統(tǒng)計(jì)樣本的均值與方差見表1,再根據(jù)式(2)進(jìn)行輻射歸一化,并計(jì)算歸一化圖像的差值,統(tǒng)計(jì)差值圖像的歸一化直方圖,同時(shí)擬合得到高斯概率分布函數(shù)(圖3),表1給出了擬合的高斯概率函數(shù)的參數(shù)。根據(jù)高斯擬合曲線計(jì)算置信度為0.95條件下的差值圖像的分割閾值(表1),印度地震和日本地震分別對應(yīng)差值為10和22。最后,根據(jù)閾值對差值圖像進(jìn)行閾值分割提取震害,差值大于閾值的像素作為提取的震害區(qū)域(圖5)。

    日本地震提取結(jié)果表明,破壞燈光主要分布在青森縣、秋田縣、巖手縣、宮城縣、茨城縣(圖5),其中宮城縣和茨城縣破壞燈光面積最大,巖手縣和宮城縣破壞燈光比例最大,達(dá)90%以上。圖6a是基于燈光變化的震害破壞范圍和地震動(dòng)分布圖。從圖中可以看出,宮城縣和茨城縣中部地震動(dòng)達(dá)到55~63gal,對應(yīng)位置的燈光也出現(xiàn)大范圍破壞;巖手縣破壞燈光呈現(xiàn)帶狀分布,對應(yīng)地震動(dòng)的高值區(qū),強(qiáng)烈的地震動(dòng)導(dǎo)致大規(guī)模的電力設(shè)施破壞,因此,燈光破壞與地震動(dòng)強(qiáng)度存在較好的一致性。本文提取的巖手縣、宮城縣的破壞燈光面積最大,與震前燈光遙感圖像相比,破壞比例達(dá)到90%以上,這與日本氣象廳發(fā)布的巖手縣和宮城縣受災(zāi)范圍大,全縣停電相一致;茨城縣破壞燈光較多,青森縣、秋田縣提取的破壞燈光相對較少,福島縣、山形縣、新潟縣只是出現(xiàn)少量燈光破壞,與日本氣象廳發(fā)布發(fā)布的局部停電也大致一致。由于本文使用的3月12日DMSP/OLS燈光圖像的空間采樣為1km,燈光遙感分辨率低使得這些縣內(nèi)的一些人口密度小的地區(qū)燈光無法識(shí)別,無法提取燈光破壞。對比圖2a,b,c和圖6b可以看出,正常情況下震前燈光遙感圖像與人口密度空間分布具有非常高的相關(guān)性。對日本全年平均燈光強(qiáng)度與人口密度空間分布進(jìn)行統(tǒng)計(jì)分析(以縣為統(tǒng)計(jì)單元),見圖7。從圖7也可以看出,燈光強(qiáng)度與人口密度具有很好的相關(guān)性,因此,通過燈光破壞區(qū)域提取也能間接的獲取地震受災(zāi)人口的規(guī)模和空間分布,這也是利用燈光遙感圖像提取震災(zāi)害信息的優(yōu)勢。

    從利用印度地震前后的燈光遙感圖像檢測結(jié)果來看,震中附近破壞嚴(yán)重的普杰、巴庫和安杰爾3個(gè)地區(qū)破壞燈光均被提取出來,這與Kohiya-ma等(2004)提取結(jié)果一致。由于沒有更多的震后調(diào)查資料,本文不做詳細(xì)的論證,但從圖3c地震前后的差值圖像的統(tǒng)計(jì)特征來看,和“3·11”日本地震前后差值圖像具有相似性。

    從“3·11”日本地震和印度地震應(yīng)用效果來看,利用燈光遙感圖像提取的破壞燈光能夠從宏觀上反映地震災(zāi)害的總體情況,特別是燈光破壞強(qiáng)弱及空間分布與受災(zāi)人口規(guī)模與空間分布也密切相關(guān),對于震后第一時(shí)間的應(yīng)急救援具有重要意義。但目前使用的DMSP/OLS燈光數(shù)據(jù)還存在一些缺陷,無論是空間分辨率還是輻射分辨率都比較低,不能滿足對燈光分布和強(qiáng)度的精細(xì)化觀測:(1)空間分辨率在星下點(diǎn)為1km,邊緣圖像分辨率2.7km,數(shù)據(jù)量化等級為6bit,DN值0~64,這極大降低了數(shù)據(jù)對燈光的描述能力;(2)由于數(shù)據(jù)不能進(jìn)行絕對輻射定標(biāo)和大氣校正等定量化處理,使得數(shù)據(jù)不能定量反演災(zāi)害程度:(3)燈光數(shù)據(jù)受云的影響大,遙感器波段設(shè)置應(yīng)該考慮去云波段的設(shè)置。我國有關(guān)部門正在進(jìn)行燈光遙感衛(wèi)星的論證,應(yīng)提高空間分辨率、輻射分辨率,同時(shí)提高數(shù)據(jù)定量化處理和去云處理的能力。

    4 結(jié)論

    燈光是人類活動(dòng)的標(biāo)志,強(qiáng)烈地震發(fā)生后的燈光變化反映了電力供應(yīng)與使用方面的變化,本質(zhì)上是災(zāi)害造成的發(fā)電、輸電等設(shè)施和線路的破壞。因此,燈光遙感影像可以作為地震破壞的一種顯性載體。本文基于對燈光遙感影像和震后燈光變化原因分析,提供了一種利用燈光遙感影像進(jìn)行地震災(zāi)害信息提取的方法,并在2001年印度地震、2011年“3·11”日本地震進(jìn)行應(yīng)用試驗(yàn),發(fā)現(xiàn)燈光遙感圖像目標(biāo)單一,除了受云等噪聲影響,主要包括燈光和背景信息,可以作為一種地震災(zāi)情快速獲取的信息來源,利用燈光遙感影像不僅能夠準(zhǔn)確提取地震破壞范圍和強(qiáng)度,而且能夠直接反映了震后受災(zāi)人口規(guī)模與空間分布,對于震后應(yīng)急救援具有重要意義。

    2011年“3·11”日本大地震的地震動(dòng)數(shù)據(jù)由日本消防廳Hasokawa博士提供,審稿專家提出了中肯的修改意見和建議,在此表示衷心感謝!

    參考文獻(xiàn):

    陳晉,卓莉,史培軍,等,2003.基于DMSP/OLS數(shù)據(jù)的中國城市化過程研究——反映區(qū)域城市化水平的燈光指數(shù)的構(gòu)建[J].遙感學(xué)報(bào),7(3):168-175.

    單新建,李建華,馬超.2005.昆侖山口西MS8.1級地震地表破裂帶高分辨率衛(wèi)星影像特征研究[J].地球物理學(xué)報(bào),48(2):321-326.

    高波,吳濤,李亮.2012.基于假設(shè)檢驗(yàn)的遙感影像變化檢測[J].測繪技術(shù)裝備,14(2):38-41.

    郭華東,陳玉,馮強(qiáng).2011.2011年日本東北9.0級大地震災(zāi)區(qū)建筑物與農(nóng)田災(zāi)情遙感評估[J].科學(xué)通報(bào),56(12):895-901.

    何春陽,史培軍,李景剛,等.2016.基于DMSP/OLS夜間燈光數(shù)據(jù)和統(tǒng)日數(shù)據(jù)的中國大陸20世紀(jì)90年代城市化空間過程重建研究[J].科學(xué)通報(bào),51(7):856-661.

    李春于,代華兵.2017.基于統(tǒng)計(jì)檢驗(yàn)的面向?qū)ο蟾叻直媛蔬b感圖像森林變化檢測[J].林業(yè)科學(xué),53(5):74-81.

    李東平,劉倩倩,襲俊.2017.基于海量視頻監(jiān)控的中小地震烈度快速判斷方法嘗試[J].40(2):324-332.

    劉愛文,廈珊,呂紅山.2012.汶川地震極重災(zāi)區(qū)通信基站典型震害及原因分析[J].電信工程技術(shù)與標(biāo)準(zhǔn)化,(12):82-86.

    陸鳴,季鴻晶,溫增平,等.2008.都江堰市移動(dòng)通信系統(tǒng)及其建筑物震害特征[J].北京工業(yè)大學(xué)學(xué)報(bào),34(12):1160-1165.

    王曉青,竇愛霞,王龍,等.2015.2013年四川蘆山7.0級地震烈度遙感評估[J].地球物理學(xué)報(bào),58(1):163-171.

    張景發(fā),謝禮立,陶夏新,2002.建筑物震害遙感圖像的變化檢測與震害評估[J].自然災(zāi)害學(xué)報(bào),11(2):59-64.

    張小詠,李慶亭,高娜.2016.基于中分辨率遙感影像的居民區(qū)震害信息提取[J].地震學(xué)報(bào),38(3):486-495.

    Elvidge C P,Bangh K E,Kihn E A,et al.1997.Relation between satelliteobserved visible-near infrared emissions,population,economic ac-tivity and electric power consumption[J],International Journal ofRemote Sensing,18(6):1373-1379.

    Hayashi H,Hashiteras,Kohiyama M,et al.2000.International coRabora-tion for the early damaged area estimation system using DMSP/OLSnighttime images[C].IEEE 2000 International Geoscience and Re-mote Sensing Symposium,6:2697-2699.

    Kimberly B,Christopher D,Tdottama G,et al.2010.Development of a2009 stable lights product using DMSP/OLS data[J].Proceedings ofthe Asia-Pacific Advanced Network,30:114-130.

    Kohiyama M,Hayashi H,Maki N,et al.2004.Early damaged~estima-tion system using DMSP-OLS night-time imagery J].Internation-al Journal of Remote Sensing,25(11):2013-2036.

    Liu A W,Xias,Lv H S,et al.2012.Typical damage analysis far mobilecommunication base stations in the extremely damage area of Wen-churn earthquake[J].Telecom Engineering Technics and Standardi-zation,12:1-5.

    Sutton P.1997.Modeling population density with nighttime satellite im-agery and GIS[J]·Computers Environment and Urban SysteMS,21(3-4):227-244.

    猜你喜歡
    差值燈光災(zāi)害
    水中燈光秀
    河南鄭州“7·20”特大暴雨災(zāi)害的警示及應(yīng)對
    今晚的燈光亮了天
    心聲歌刊(2021年4期)2021-10-13 08:31:38
    我省汛期常見氣象災(zāi)害及防御
    來一場燈光派對
    差值法巧求剛體轉(zhuǎn)動(dòng)慣量
    推動(dòng)災(zāi)害防治工作實(shí)現(xiàn)新跨越
    枳殼及其炮制品色差值與化學(xué)成分的相關(guān)性
    中成藥(2017年6期)2017-06-13 07:30:35
    燈光閃閃
    基于區(qū)域最大值與平均值差值的動(dòng)態(tài)背光調(diào)整
    99国产精品免费福利视频| 国内精品久久久久久久电影| 免费在线观看日本一区| 午夜激情av网站| 亚洲欧美日韩无卡精品| 亚洲av电影在线进入| 嫁个100分男人电影在线观看| 老司机福利观看| 午夜日韩欧美国产| 成人av一区二区三区在线看| 亚洲成av片中文字幕在线观看| 日韩av在线大香蕉| 国产精品亚洲av一区麻豆| 国产精品野战在线观看| 身体一侧抽搐| 99精品欧美一区二区三区四区| a级毛片在线看网站| 成人三级做爰电影| 欧美 亚洲 国产 日韩一| 此物有八面人人有两片| 国产一区二区在线av高清观看| 亚洲国产精品成人综合色| 首页视频小说图片口味搜索| 欧美乱码精品一区二区三区| 久久人人精品亚洲av| 色综合站精品国产| 国产极品粉嫩免费观看在线| 热re99久久国产66热| 女人被躁到高潮嗷嗷叫费观| 99国产精品免费福利视频| 亚洲成av人片免费观看| 正在播放国产对白刺激| 亚洲专区国产一区二区| 婷婷六月久久综合丁香| 久久亚洲真实| 亚洲欧洲精品一区二区精品久久久| 又黄又爽又免费观看的视频| 亚洲精品久久成人aⅴ小说| 国产伦人伦偷精品视频| 亚洲一区高清亚洲精品| 欧美乱色亚洲激情| 无遮挡黄片免费观看| 天天添夜夜摸| 欧美丝袜亚洲另类 | 亚洲国产精品sss在线观看| 黑丝袜美女国产一区| 国产精品美女特级片免费视频播放器 | 黑人欧美特级aaaaaa片| 日本黄色视频三级网站网址| 欧美日本中文国产一区发布| 国产极品粉嫩免费观看在线| 亚洲专区字幕在线| 日本五十路高清| 国产精品日韩av在线免费观看 | 亚洲自拍偷在线| 俄罗斯特黄特色一大片| 成人国语在线视频| 搞女人的毛片| 国产精品久久久人人做人人爽| 亚洲五月天丁香| 女人被狂操c到高潮| 亚洲国产日韩欧美精品在线观看 | 亚洲中文av在线| 在线观看免费日韩欧美大片| 亚洲情色 制服丝袜| 亚洲中文字幕日韩| 免费高清视频大片| 久久香蕉激情| 久久国产亚洲av麻豆专区| 亚洲全国av大片| 亚洲无线在线观看| 精品一区二区三区四区五区乱码| 日韩国内少妇激情av| 97人妻精品一区二区三区麻豆 | 久久久久九九精品影院| 久久人妻av系列| 亚洲三区欧美一区| 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区精品| 国产精品野战在线观看| 国产男靠女视频免费网站| 国产精品自产拍在线观看55亚洲| x7x7x7水蜜桃| 免费在线观看视频国产中文字幕亚洲| 丰满的人妻完整版| 久久午夜亚洲精品久久| 看片在线看免费视频| 日韩免费av在线播放| 欧美激情极品国产一区二区三区| 亚洲免费av在线视频| 国产精品亚洲美女久久久| 无遮挡黄片免费观看| 搡老妇女老女人老熟妇| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 日本精品一区二区三区蜜桃| 黄色片一级片一级黄色片| 久久久久久久久中文| 午夜精品国产一区二区电影| 满18在线观看网站| 成人18禁在线播放| 亚洲精品中文字幕在线视频| 日韩视频一区二区在线观看| 国产激情久久老熟女| 侵犯人妻中文字幕一二三四区| 午夜福利,免费看| 香蕉国产在线看| 可以在线观看毛片的网站| 男人舔女人的私密视频| 久久人妻熟女aⅴ| 成人三级做爰电影| 成人特级黄色片久久久久久久| 午夜福利,免费看| 欧美日韩福利视频一区二区| 国产精品乱码一区二三区的特点 | 日日夜夜操网爽| 一区二区三区高清视频在线| 久久人妻福利社区极品人妻图片| 国产人伦9x9x在线观看| 欧美亚洲日本最大视频资源| 久久久久久国产a免费观看| 久久久精品国产亚洲av高清涩受| 欧美不卡视频在线免费观看 | 一区福利在线观看| 黄片大片在线免费观看| 成人欧美大片| 国产真人三级小视频在线观看| 母亲3免费完整高清在线观看| 最近最新中文字幕大全电影3 | 亚洲第一av免费看| 18美女黄网站色大片免费观看| 久久人妻av系列| 日本免费一区二区三区高清不卡 | 亚洲五月色婷婷综合| 91麻豆精品激情在线观看国产| 久久九九热精品免费| 看免费av毛片| 日本五十路高清| 人成视频在线观看免费观看| 免费一级毛片在线播放高清视频 | 国内久久婷婷六月综合欲色啪| 亚洲熟女毛片儿| 久久久久久国产a免费观看| 欧美中文日本在线观看视频| 又紧又爽又黄一区二区| 女性生殖器流出的白浆| 国产成人欧美在线观看| 变态另类成人亚洲欧美熟女 | 91九色精品人成在线观看| 美女免费视频网站| 国产成人影院久久av| 99精品在免费线老司机午夜| 国产一区在线观看成人免费| 日韩有码中文字幕| 欧美 亚洲 国产 日韩一| 国产精品久久久人人做人人爽| 99在线视频只有这里精品首页| 精品高清国产在线一区| 日韩精品免费视频一区二区三区| 亚洲人成电影免费在线| 亚洲一区二区三区不卡视频| 欧美黄色片欧美黄色片| 国产亚洲欧美98| 欧美日韩精品网址| 亚洲中文av在线| 亚洲欧美精品综合久久99| 久热爱精品视频在线9| 天堂√8在线中文| e午夜精品久久久久久久| 美女国产高潮福利片在线看| 人人妻人人爽人人添夜夜欢视频| 非洲黑人性xxxx精品又粗又长| 成人精品一区二区免费| 久久人人97超碰香蕉20202| videosex国产| 一级作爱视频免费观看| 最近最新中文字幕大全电影3 | 欧美av亚洲av综合av国产av| 桃红色精品国产亚洲av| 搡老熟女国产l中国老女人| 一区福利在线观看| bbb黄色大片| 99国产精品一区二区三区| 在线免费观看的www视频| 国产国语露脸激情在线看| 不卡av一区二区三区| 国产精品免费一区二区三区在线| 首页视频小说图片口味搜索| 一区二区三区高清视频在线| 色播亚洲综合网| 亚洲欧美日韩高清在线视频| 久久人人精品亚洲av| 99riav亚洲国产免费| 中文字幕高清在线视频| 丰满的人妻完整版| 精品国产超薄肉色丝袜足j| 国产精品香港三级国产av潘金莲| 久久精品国产亚洲av高清一级| 黄片播放在线免费| 亚洲第一av免费看| 别揉我奶头~嗯~啊~动态视频| 日韩国内少妇激情av| 男女床上黄色一级片免费看| 日韩欧美在线二视频| 少妇粗大呻吟视频| 亚洲五月天丁香| 丝袜美足系列| 中文字幕人妻丝袜一区二区| 精品日产1卡2卡| 999久久久精品免费观看国产| 美女高潮到喷水免费观看| 91麻豆av在线| 97碰自拍视频| 真人一进一出gif抽搐免费| 日韩国内少妇激情av| 免费高清视频大片| 多毛熟女@视频| 亚洲专区中文字幕在线| 久久精品成人免费网站| 咕卡用的链子| 啦啦啦 在线观看视频| 国产激情久久老熟女| 亚洲 国产 在线| 久久精品91无色码中文字幕| 欧美国产精品va在线观看不卡| 长腿黑丝高跟| 日韩欧美国产在线观看| 国产野战对白在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 变态另类成人亚洲欧美熟女 | 一二三四在线观看免费中文在| 久久久久久国产a免费观看| 亚洲第一电影网av| 国产野战对白在线观看| 午夜精品久久久久久毛片777| 欧美性长视频在线观看| 成人免费观看视频高清| 91大片在线观看| 波多野结衣av一区二区av| 香蕉国产在线看| 色婷婷久久久亚洲欧美| 成人亚洲精品一区在线观看| 亚洲在线自拍视频| 婷婷六月久久综合丁香| 久久精品国产99精品国产亚洲性色 | 91麻豆精品激情在线观看国产| 女人精品久久久久毛片| 色哟哟哟哟哟哟| 日日干狠狠操夜夜爽| 怎么达到女性高潮| 欧美午夜高清在线| 色综合欧美亚洲国产小说| 黄色 视频免费看| 两个人视频免费观看高清| 国产精品电影一区二区三区| 亚洲第一青青草原| 国产亚洲欧美精品永久| 亚洲少妇的诱惑av| 精品欧美国产一区二区三| 国产乱人伦免费视频| 国产av精品麻豆| 日韩大尺度精品在线看网址 | 日韩大尺度精品在线看网址 | 99国产精品99久久久久| 国产成年人精品一区二区| 欧美一级毛片孕妇| 黑人巨大精品欧美一区二区蜜桃| 99精品久久久久人妻精品| 又紧又爽又黄一区二区| 欧美在线一区亚洲| 中文字幕精品免费在线观看视频| 91成年电影在线观看| 国产精品一区二区在线不卡| av欧美777| 很黄的视频免费| 久久久久久免费高清国产稀缺| 人成视频在线观看免费观看| 欧美国产精品va在线观看不卡| videosex国产| 黄色毛片三级朝国网站| 88av欧美| 两个人看的免费小视频| 免费av毛片视频| 欧美一级毛片孕妇| 91成人精品电影| 欧美色视频一区免费| 一级a爱视频在线免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 国产乱人伦免费视频| 国内毛片毛片毛片毛片毛片| 亚洲av成人一区二区三| 黄色视频不卡| 怎么达到女性高潮| 深夜精品福利| 国产日韩一区二区三区精品不卡| 国产一卡二卡三卡精品| 欧洲精品卡2卡3卡4卡5卡区| 97人妻精品一区二区三区麻豆 | √禁漫天堂资源中文www| 亚洲少妇的诱惑av| 老司机午夜十八禁免费视频| 神马国产精品三级电影在线观看 | 午夜福利在线观看吧| 午夜福利影视在线免费观看| 国产精品av久久久久免费| 可以在线观看毛片的网站| 日韩欧美三级三区| www.熟女人妻精品国产| 精品一品国产午夜福利视频| 国产激情久久老熟女| 亚洲国产精品999在线| 国内精品久久久久久久电影| 一进一出抽搐gif免费好疼| 日韩大码丰满熟妇| 亚洲午夜理论影院| 嫩草影院精品99| 国产成人免费无遮挡视频| 午夜福利,免费看| 国产高清videossex| 美女免费视频网站| 欧美一级毛片孕妇| 色在线成人网| 老司机午夜福利在线观看视频| 欧美+亚洲+日韩+国产| 久久香蕉激情| 亚洲专区中文字幕在线| 黄色片一级片一级黄色片| 老汉色∧v一级毛片| 大香蕉久久成人网| 亚洲人成网站在线播放欧美日韩| a在线观看视频网站| 一二三四在线观看免费中文在| 久久久久久久精品吃奶| 亚洲国产精品合色在线| 亚洲七黄色美女视频| 国产亚洲精品av在线| 久久香蕉国产精品| 可以在线观看毛片的网站| 国产精品二区激情视频| 欧美久久黑人一区二区| 亚洲av美国av| 脱女人内裤的视频| 国产一区二区三区综合在线观看| 中文字幕av电影在线播放| 午夜福利,免费看| 久久影院123| 久久精品国产清高在天天线| 午夜久久久久精精品| 国产高清激情床上av| 婷婷丁香在线五月| 国产三级在线视频| 精品久久蜜臀av无| 午夜a级毛片| 国产激情久久老熟女| 久久狼人影院| www.熟女人妻精品国产| 国产精品免费一区二区三区在线| 黄色a级毛片大全视频| 99热只有精品国产| 国产成人av激情在线播放| 美女高潮喷水抽搐中文字幕| 精品一区二区三区av网在线观看| 日本精品一区二区三区蜜桃| 亚洲精品av麻豆狂野| 老汉色av国产亚洲站长工具| 一级毛片精品| videosex国产| 成人永久免费在线观看视频| 国产一区二区三区综合在线观看| 黄色毛片三级朝国网站| 黑丝袜美女国产一区| 亚洲国产日韩欧美精品在线观看 | 在线av久久热| 精品国产一区二区三区四区第35| 亚洲avbb在线观看| 精品国产美女av久久久久小说| 丰满人妻熟妇乱又伦精品不卡| 少妇被粗大的猛进出69影院| 久久久久久久久免费视频了| 国产黄a三级三级三级人| 一级毛片精品| 久久中文字幕人妻熟女| 色婷婷久久久亚洲欧美| av超薄肉色丝袜交足视频| 黄色成人免费大全| 丝袜人妻中文字幕| 动漫黄色视频在线观看| 欧美人与性动交α欧美精品济南到| 999久久久精品免费观看国产| 男人操女人黄网站| 此物有八面人人有两片| 久久精品国产99精品国产亚洲性色 | 制服丝袜大香蕉在线| 亚洲专区字幕在线| 亚洲国产精品成人综合色| 777久久人妻少妇嫩草av网站| 色精品久久人妻99蜜桃| 久久香蕉精品热| 禁无遮挡网站| 女人被躁到高潮嗷嗷叫费观| 日本vs欧美在线观看视频| 涩涩av久久男人的天堂| 亚洲精品国产一区二区精华液| 婷婷六月久久综合丁香| 淫妇啪啪啪对白视频| 日韩三级视频一区二区三区| 精品国产国语对白av| 成在线人永久免费视频| 亚洲一区中文字幕在线| a在线观看视频网站| 国产一级毛片七仙女欲春2 | 丝袜美足系列| а√天堂www在线а√下载| 亚洲av五月六月丁香网| 99riav亚洲国产免费| 亚洲精品国产区一区二| 亚洲人成伊人成综合网2020| 丁香欧美五月| 在线播放国产精品三级| 长腿黑丝高跟| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品久久久久5区| 日韩国内少妇激情av| 一个人观看的视频www高清免费观看 | 亚洲视频免费观看视频| 我的亚洲天堂| 日韩视频一区二区在线观看| 色老头精品视频在线观看| 亚洲av熟女| 亚洲av五月六月丁香网| 欧美成狂野欧美在线观看| ponron亚洲| 97超级碰碰碰精品色视频在线观看| 国产精品久久久av美女十八| 高清在线国产一区| 十八禁网站免费在线| 精品高清国产在线一区| 亚洲精品中文字幕一二三四区| 九色亚洲精品在线播放| 激情在线观看视频在线高清| 亚洲七黄色美女视频| 一二三四在线观看免费中文在| 午夜福利视频1000在线观看 | 午夜福利成人在线免费观看| 亚洲,欧美精品.| 国产激情欧美一区二区| 久久中文字幕人妻熟女| 国产精品一区二区精品视频观看| 大香蕉久久成人网| 国产精品久久久久久亚洲av鲁大| 久久草成人影院| 久久婷婷人人爽人人干人人爱 | 国产精品精品国产色婷婷| 欧美性长视频在线观看| 午夜免费激情av| 亚洲全国av大片| 国产单亲对白刺激| 国产91精品成人一区二区三区| 精品国内亚洲2022精品成人| 久久久精品欧美日韩精品| 超碰成人久久| 欧美成人一区二区免费高清观看 | 亚洲成国产人片在线观看| 久久精品aⅴ一区二区三区四区| 人人澡人人妻人| 国产成人精品在线电影| 免费无遮挡裸体视频| 一边摸一边抽搐一进一小说| 天天躁夜夜躁狠狠躁躁| 婷婷丁香在线五月| 岛国视频午夜一区免费看| 精品国内亚洲2022精品成人| 日本免费a在线| 免费无遮挡裸体视频| 欧美乱妇无乱码| 日本五十路高清| 国产成人啪精品午夜网站| 在线免费观看的www视频| 亚洲中文日韩欧美视频| 亚洲aⅴ乱码一区二区在线播放 | 女同久久另类99精品国产91| 人人妻人人澡人人看| 久久国产精品男人的天堂亚洲| 在线播放国产精品三级| 视频区欧美日本亚洲| 又黄又爽又免费观看的视频| 午夜福利视频1000在线观看 | 国产精品久久电影中文字幕| 欧美日本中文国产一区发布| 欧美性长视频在线观看| 两性夫妻黄色片| 人人妻人人澡欧美一区二区 | 在线天堂中文资源库| 又黄又粗又硬又大视频| 一区二区三区精品91| 亚洲精品在线观看二区| 国产精品永久免费网站| 亚洲无线在线观看| 国产精品1区2区在线观看.| 在线观看日韩欧美| 亚洲国产精品999在线| 精品第一国产精品| 叶爱在线成人免费视频播放| 在线观看免费视频日本深夜| 欧美人与性动交α欧美精品济南到| 国产成人欧美| 国产欧美日韩综合在线一区二区| 搡老岳熟女国产| 国产高清videossex| 91麻豆精品激情在线观看国产| 国产精品一区二区在线不卡| 成人亚洲精品一区在线观看| av免费在线观看网站| 国产欧美日韩综合在线一区二区| av欧美777| 日日爽夜夜爽网站| 亚洲国产精品成人综合色| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| 免费在线观看完整版高清| 国产熟女xx| 1024视频免费在线观看| 乱人伦中国视频| www.熟女人妻精品国产| 一个人免费在线观看的高清视频| 精品福利观看| 精品人妻在线不人妻| 免费久久久久久久精品成人欧美视频| 国产一区二区在线av高清观看| 亚洲一区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 一a级毛片在线观看| 高潮久久久久久久久久久不卡| 黄频高清免费视频| 国产成人欧美| 亚洲avbb在线观看| 国产精品日韩av在线免费观看 | 国产熟女午夜一区二区三区| 中文字幕人成人乱码亚洲影| 怎么达到女性高潮| 久久久久久久精品吃奶| 亚洲全国av大片| 国产精品 欧美亚洲| 国产一区二区在线av高清观看| 手机成人av网站| av福利片在线| 两个人免费观看高清视频| 丝袜美足系列| 久久久久九九精品影院| 欧美av亚洲av综合av国产av| 怎么达到女性高潮| 人成视频在线观看免费观看| 91字幕亚洲| 国内久久婷婷六月综合欲色啪| 国产精品日韩av在线免费观看 | 人人妻人人澡人人看| 黄色女人牲交| 精品久久久精品久久久| 欧美成狂野欧美在线观看| 丁香欧美五月| 可以在线观看毛片的网站| 老鸭窝网址在线观看| 神马国产精品三级电影在线观看 | 欧美日韩亚洲国产一区二区在线观看| 成人免费观看视频高清| 久热爱精品视频在线9| 操美女的视频在线观看| 精品乱码久久久久久99久播| 欧美乱码精品一区二区三区| 久久午夜亚洲精品久久| 国产精品香港三级国产av潘金莲| 中文字幕高清在线视频| 亚洲人成电影观看| 99国产极品粉嫩在线观看| 亚洲精品国产色婷婷电影| 国内毛片毛片毛片毛片毛片| 国产成人欧美| 欧洲精品卡2卡3卡4卡5卡区| 久久久精品欧美日韩精品| 精品国产亚洲在线| 亚洲第一欧美日韩一区二区三区| 国产一区二区在线av高清观看| 国产精品免费视频内射| 久久久久久亚洲精品国产蜜桃av| 亚洲熟妇中文字幕五十中出| 欧美激情久久久久久爽电影 | 国产主播在线观看一区二区| 亚洲欧美日韩高清在线视频| av网站免费在线观看视频| 亚洲av成人不卡在线观看播放网| 久久人妻av系列| 久久精品国产亚洲av高清一级| 亚洲全国av大片| 成人18禁高潮啪啪吃奶动态图| 久久人妻福利社区极品人妻图片| 12—13女人毛片做爰片一| 国内久久婷婷六月综合欲色啪| 露出奶头的视频| 999久久久精品免费观看国产| 免费久久久久久久精品成人欧美视频| 免费高清在线观看日韩| 国产一区二区三区视频了| 午夜成年电影在线免费观看| 欧美日韩亚洲综合一区二区三区_| 久久婷婷成人综合色麻豆| 国产精品 欧美亚洲| 中文字幕另类日韩欧美亚洲嫩草| 色哟哟哟哟哟哟| 国内久久婷婷六月综合欲色啪| 国产伦人伦偷精品视频| 天天躁夜夜躁狠狠躁躁| 最近最新中文字幕大全电影3 | 久9热在线精品视频| 在线永久观看黄色视频| 国产精品乱码一区二三区的特点 |