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

    黃海綠潮災(zāi)害應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警系統(tǒng)

    2015-12-23 07:43:21吳玲娟曹雅靜靳熙芳
    關(guān)鍵詞:綠潮黃海應(yīng)急

    吳玲娟,高 松,丁 一,曹雅靜,靳熙芳

    (1.國家海洋局北海預(yù)報中心,山東青島 266061;2.山東省海洋生態(tài)環(huán)境與防災(zāi)減災(zāi)重點實驗室,山東青島 266061)

    黃海綠潮災(zāi)害應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警系統(tǒng)

    吳玲娟1,2,高 松1,2,丁 一1,2,曹雅靜1,2,靳熙芳1,2

    (1.國家海洋局北海預(yù)報中心,山東青島 266061;2.山東省海洋生態(tài)環(huán)境與防災(zāi)減災(zāi)重點實驗室,山東青島 266061)

    黃海綠潮災(zāi)害應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警系統(tǒng)是在綠潮應(yīng)急遙感監(jiān)測技術(shù)、應(yīng)急快速漂移預(yù)測等技術(shù)研究的基礎(chǔ)上,基于GIS技術(shù)構(gòu)建了集綠潮遙感信息解譯和提取、多源監(jiān)測數(shù)據(jù)融合、快速漂移預(yù)測和預(yù)警產(chǎn)品制作和發(fā)布綜合業(yè)務(wù)化平臺。在應(yīng)急遙感監(jiān)測數(shù)據(jù)提取方面主要研發(fā)了信息提取準(zhǔn)自動化技術(shù)和多源多時間窗的綠潮監(jiān)測數(shù)據(jù)融合技術(shù);在應(yīng)急預(yù)測技術(shù)方面主要基于精細化大氣和海洋環(huán)境動力場,研發(fā)了綠潮快速應(yīng)急漂移預(yù)測模型,并通過海上和后報實驗提高模擬精度。本系統(tǒng)服務(wù)于國家海洋局和政府防災(zāi)減災(zāi)應(yīng)急部門,為海上運動、水產(chǎn)養(yǎng)殖、濱海旅游、交通運輸、漁業(yè)生產(chǎn)等涉?;顒雍秃Q蠊芾硖峁┝丝茖W(xué)可靠的技術(shù)保障,特別在2008年奧運會帆船/板賽和2012年亞洲沙灘運動會的綠潮監(jiān)測和預(yù)測預(yù)警保障工作中發(fā)揮了重要的作用,具有巨大的社會、經(jīng)濟和生態(tài)效益。

    綠潮;遙感監(jiān)測;預(yù)測預(yù)警;系統(tǒng)

    0 引言

    近7年我國黃海海域綠潮連續(xù)暴發(fā),對2008年青島奧運會帆船/板賽和2012年的海陽亞洲沙灘運動會等諸多海上賽事造成了困難,對水產(chǎn)養(yǎng)殖、濱海旅游、海上交通運輸?shù)认嚓P(guān)產(chǎn)業(yè)的影響十分嚴(yán)重。為了保證賽事、沿海經(jīng)濟和生態(tài)環(huán)境可持續(xù)發(fā)展,我國相關(guān)部門給予了高度重視,啟動相應(yīng)的應(yīng)急預(yù)案,應(yīng)對綠潮災(zāi)害;同時啟動了多項科研和應(yīng)用項目,如2008綠潮專項;2008年“滸苔大規(guī)模暴發(fā)應(yīng)急處置關(guān)鍵技術(shù)研究與應(yīng)用”;2012年國家海洋局公益性項目“黃海綠潮業(yè)務(wù)化預(yù)測預(yù)警關(guān)鍵技術(shù)研究和應(yīng)用”。

    國內(nèi)外學(xué)者主要在大型藻類及滸苔綠潮的遙感監(jiān)測方面進行研究。Gower等[1]第一次用MERIS 300m分辨率數(shù)據(jù),采用MCI(最大葉綠素指數(shù))算法監(jiān)測到墨西哥灣大面積線形馬尾藻類海草,但MERIS衛(wèi)星傳感器掃描寬度窄,限制了遙感監(jiān)測范圍。張娟[2],李三妹[3]和Elijah[4]等根據(jù)滸苔光譜特性建立了其遙感監(jiān)測模型。李穎等[5]利用不同3種星載雷達數(shù)據(jù),判讀雷達圖像上滸苔分布的典型區(qū)域,對滸苔的灰度值和正常海水的灰度值進行分析,得出不同極化方式下滸苔灰度值的主要動態(tài)范圍。蔣興偉等[6]運用一種基于區(qū)域增長面向?qū)ο蟮膱D像尺度分割方法調(diào)整圖像的分割尺度,實現(xiàn)了滸苔信息的快速提取。Gui等[7]利用ENVISAT-ASAR與環(huán)境災(zāi)害小衛(wèi)星(HJ-1A/ B)監(jiān)測2009年黃海滸苔爆發(fā)過程,指出二者平均差異為15%。Hu等[8-9]提出了一種新的海洋浮游藻類指數(shù)(FAI),并運用這個指數(shù)監(jiān)測2008年我國青島海域滸苔,結(jié)果顯示該方法比NDVI(Normalized Difference Vegetation Index)方法更加準(zhǔn)確。Keesing等[10]基于MODIS獲取的灰度值方法,研究了2007~2009年黃海滸苔綠潮的年際變化,分析了其與海岸帶水藻養(yǎng)殖的關(guān)系?;诓ǘ伪戎捣ǖ倪b感藻類指數(shù),Shi等[11]提出的歸一化藻類指數(shù) NDAI。國家海洋局北海預(yù)報中心用RGB圖像和歸一化植被值數(shù)NDVI方法,對MODIS、SAR COSMO和HY-1B等衛(wèi)星遙感數(shù)據(jù)進行提取,并融合航空遙感、船舶、陸岸等其他監(jiān)測結(jié)果,得到綠潮的綜合監(jiān)測

    結(jié)果,為綠潮漂移預(yù)測和溯源提供初始場,同時為研究綠潮分布、發(fā)生發(fā)展和漂移機理等工作提供數(shù)據(jù)[12-14]。國內(nèi)外學(xué)者在綠潮遙感監(jiān)測方面主要集中在綠潮遙感數(shù)據(jù)的提取,但是綠潮遙感監(jiān)測信息自動提取和數(shù)據(jù)融合鮮為報道。

    近年來國內(nèi)外科學(xué)家對綠潮的漂移和生長機制等要素進行多方面研究,并得到了初步的成果,但是對綠潮漂移快速模擬的研究較少。Aurousseau等[15]在法國綠潮災(zāi)害比較嚴(yán)重的Brest灣,建立三維生物地球化學(xué)模式,對該海域的漂浮綠潮的生長和腐?。疲两颠M行模擬。Cugier等[16]建立的浮游植物三維生態(tài)學(xué)模型和三維水動力模型[17],建立適應(yīng)于Brest灣的綠潮特征[18]的三維生物地球化學(xué)模型,同時開展了潮間帶綠潮的預(yù)報研究。Perrot等[19]發(fā)展了預(yù)報潮間帶海藻的簡單方法,并于2007年建立了綠潮在潮間帶生長和漂移的預(yù)報模式。喬方利等[20]利用數(shù)值模式對2010年綠潮漂移進行了模擬研究,發(fā)現(xiàn)風(fēng)場驅(qū)動下的海洋表層流場年際變化是綠潮漂移路徑變異的主要原因。國家海洋環(huán)境預(yù)報中心[21]2008年建立了滸苔漂移路徑預(yù)測系統(tǒng),向有關(guān)部門發(fā)布滸苔漂移軌跡預(yù)報和海洋環(huán)境預(yù)報,協(xié)助前線應(yīng)急指揮中心完成了滸苔的控制和治理。國家海洋局北海預(yù)報中心利用綠潮綜合監(jiān)測數(shù)據(jù),啟動應(yīng)急漂移模式,快速預(yù)測和追溯綠潮的漂移方向和路徑[22-23]。

    目前國內(nèi)外綠潮的監(jiān)測和數(shù)值模擬已取得一定的研究成果,但從2008年黃海綠潮暴發(fā)開始,綠潮應(yīng)急遙感監(jiān)測和預(yù)測關(guān)鍵技術(shù)不成熟,并缺少業(yè)務(wù)化的綜合平臺,因此,北海預(yù)報中心研發(fā)綠潮應(yīng)急監(jiān)測和預(yù)測預(yù)警關(guān)鍵技術(shù),并基于GIS技術(shù)建設(shè)集綠潮遙感信息解譯和提取、多源監(jiān)測數(shù)據(jù)融合、快速漂移預(yù)測和預(yù)警產(chǎn)品制作和發(fā)布于一體的綜合業(yè)務(wù)化,以提升我國綠潮防災(zāi)減災(zāi)能力。

    1 綠潮應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警系統(tǒng)簡介

    綠潮應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警系統(tǒng)(以下簡稱綠潮系統(tǒng))業(yè)務(wù)化運行流程如圖1。綠潮系統(tǒng)可以分為綠潮應(yīng)急遙感監(jiān)測子系統(tǒng)和應(yīng)急漂移預(yù)測子系統(tǒng),并基于兩個子系統(tǒng)建立綜合業(yè)務(wù)化平臺,對服務(wù)對象發(fā)布監(jiān)測和預(yù)警信息產(chǎn)品。

    1.1 綠潮應(yīng)急遙感監(jiān)測子系統(tǒng)

    1.1.1 綠潮信息自動化提取

    綠潮衛(wèi)星遙感數(shù)據(jù)包括光學(xué)遙感數(shù)據(jù)(MODIS-TERRA/AQUA和HJ-1A、B衛(wèi)星)和微波遙感數(shù)據(jù)(COSMO-SAR、Radarsat、Terra-SAR)。光學(xué)衛(wèi)星影像綠潮信息提取包括輻射定標(biāo)、幾何校正、大氣校正、云掩模和綠潮信息提取五個步驟,可獲取產(chǎn)品有NDVI、云、影像圖、綠潮解譯圖和綠潮光學(xué)報告;微波衛(wèi)星影像信息提取包括輻射定標(biāo)、幾何校正和綠潮信息提取三個步驟,可獲得產(chǎn)品有影像圖、解譯圖和綠潮微波報告。系統(tǒng)發(fā)展的綠潮分布輪廓線自動提取算法,經(jīng)過距離分析、邏輯運算、膨脹腐蝕運算,所提取的輪廓線效果達到了人工操作水平(圖2)。根據(jù)綠潮自動提取結(jié)果可計算綠潮覆蓋面積,并自動生成綠潮分布面并計算分布面積,并結(jié)合業(yè)務(wù)需求,可以統(tǒng)計35°以北以及沿海受災(zāi)地市劃定警戒區(qū)覆蓋面積和分布面積。

    1.1.2 多時間窗數(shù)據(jù)源信息融合

    通過綠潮應(yīng)急漂移預(yù)測和溯源技術(shù)[22-23],解決了數(shù)據(jù)源時間不同步的問題。基于不同時間(h1、h2、h3)的應(yīng)急遙感監(jiān)測結(jié)果,利用綠潮應(yīng)急快速漂移(溯源)模型將h1和h3時刻的綠潮數(shù)據(jù),快速預(yù)測或追溯到h2時刻,為遙感信息匹配融合提供不同數(shù)據(jù)源的h2時刻的綠潮模擬位置。對于高、低分辨率光學(xué)遙感監(jiān)測信息,系統(tǒng)研發(fā)融合技術(shù):建立公共區(qū)域內(nèi)低分辨率遙感影像NDVI值與高分辨率遙感影像覆蓋率的統(tǒng)計關(guān)系,修正整景低分辨率遙感影像的綠潮監(jiān)測結(jié)果。基于上述技術(shù),研發(fā)了多源衛(wèi)星影像監(jiān)測信息的融合技術(shù),其主要流程以2014年7月4日6:56 Radarsat和9:46 HJ-1A兩景監(jiān)測數(shù)據(jù)(圖3)的數(shù)據(jù)融合過程為例:1)由于兩景數(shù)據(jù)成像時刻不同,利用綠潮應(yīng)急漂移預(yù)測技術(shù),將6:56 Radarsat數(shù)據(jù)漂移到9:46時刻,解決時間不同的問題;2)判別兩類監(jiān)測解譯結(jié)果是否有重疊區(qū),如果無,則融合結(jié)果為兩種遙感數(shù)據(jù)綠潮探測結(jié)果的直接疊加;3)如果存在疊加區(qū)域,那么判斷疊加區(qū)域是否有云覆蓋,如果有云覆蓋,選用受云覆蓋影響小的影像提取結(jié)果(Radarsat影像);如果無云,那么選擇分辨率較高影像(Radarsat提取結(jié)果)(圖3);4)最終得到9:46時刻的綠潮分布情況。

    圖1 綠潮系統(tǒng)業(yè)務(wù)化運行流程Fig.1 Operational Process of Remote-sensing M onitoring,F(xiàn)orecasting and W arning System on Green Tide in Yellow Sea

    圖2 綠潮二值影像(左圖)和提取結(jié)果(右圖)Fig.2 The Binary Image of Green Tide(left)and the Extraction Resu lt(right)

    1.1.3 綜合分布圖制作

    利用綠潮信息自動化提取技術(shù)和多時間窗數(shù)據(jù)源信息融合技術(shù)提取并融合多源遙感數(shù)據(jù),獲得綠潮的分布面積和覆蓋面積等信息,最終制作綜合分布圖(圖5,紅色線代表綠潮的外緣線),為綠潮應(yīng)急預(yù)測子系統(tǒng)提供初始場、邊界場和驗證場。

    1.2 黃海綠潮快速漂移預(yù)測子系統(tǒng)

    1.2.1 黃海環(huán)境動力業(yè)務(wù)化模式

    圖3 2014年7月4日Radarsat(a)和HJ-1A(b)衛(wèi)星遙感解譯結(jié)果,Radarsat衛(wèi)星遙感數(shù)據(jù)漂移到9:46時刻的模擬結(jié)果(c,藍色點)與同時刻HJ-1A衛(wèi)星數(shù)據(jù)(c,綠色點)疊加結(jié)果(c)、數(shù)據(jù)融合結(jié)果(d)Fig.3 The Interp retation of 6:56 Radarsat(a)and 9:46 HJ-1A(b)on 4 July,2014;Superposition of HJ-1A data(c,green dots)and Drift Sim ulation Result of Radarsat Data to 9:46(c,blue dots)and Fused Result(d)

    黃海環(huán)境動力業(yè)務(wù)化模式主要包括氣象和海洋環(huán)境動力模式。氣象模式采用美國大氣研究中心開發(fā)的WRF(Weather Research and Forecasting Model),采用了多重嵌套技術(shù),建立中國海海域、東中國海海域、北海區(qū)和青島近海模型,并在青島進行加密計算,水平分辨率達3 km,同時采用三維變分方法24 h循環(huán)同化GTS、海洋站和浮標(biāo)等實時觀測數(shù)據(jù)。根據(jù)綠潮所在位置、范圍以及政府部門對應(yīng)急預(yù)測的不同需求,采用多重嵌套技術(shù),建立中國海區(qū)、北海區(qū)和青島近海小區(qū)海流模型。該模型采用WRF風(fēng)場和熱通量場為大氣強迫場,中國海區(qū)采用ROMS模式,計算范圍為99°~148°E,9°S~44.05°N,水平分辨率為0.1°× 0.1°,垂向分25層;北海區(qū)采用ROMS模型,模型范圍為117°E~127°E,32°N~41°N,水平分辨率為1/30°×1/30°,垂向分6層。近海小區(qū)海流預(yù)報采用FVCOM模式,對關(guān)注海域局部加密,分辨率最小可達10 m。同時海洋模型采用改進集合Kalman濾波同化技術(shù)實現(xiàn)對地波雷達海流實時觀測數(shù)據(jù)的同化[24],進一步提高了海洋模式的預(yù)報精度。

    1.2.2 黃海綠潮快速漂移預(yù)測模型

    在不考慮綠潮自身生態(tài)特征的情況下,其在海水中的移動,可以看作是質(zhì)點跟隨海流的物理運動,所以綠潮應(yīng)急漂移預(yù)測,采用拉格朗日粒子追蹤方法?;诶窭嗜沼^點,粒子的漂移速度V?L計算公式為:

    式中:V?W—由風(fēng)力和波浪作用產(chǎn)生的速度分量;V?t—潮流作用產(chǎn)生的速度分量;V?r—潮致余流作用產(chǎn)生的速度分量;V?h—環(huán)流(包括:風(fēng)海流和密度流)作用產(chǎn)生的速度分量;潮流流速分量V?t和V?r,由潮流調(diào)和常數(shù)預(yù)報得到。環(huán)流流速分量

    V?h,綠潮環(huán)境動力業(yè)務(wù)化預(yù)報系統(tǒng)提供。

    基于黃海綠潮環(huán)境動力數(shù)值模擬結(jié)果,利用綠潮斑塊海上漂移實驗修正了海流和風(fēng)力系數(shù)的比值[25],并分析風(fēng)拖曳角度變化對綠潮分布的影響,確定模型參數(shù)。該模式能自動搜索綠潮所在區(qū)域,并考慮圍欄和流網(wǎng)障礙物的阻擋,快速預(yù)測綠潮斑塊漂移路徑[22]。

    1.2.3 綠潮漂移預(yù)測結(jié)果驗證

    由于綠潮(滸苔)的繁殖擴展速度非常之快,相差一天的滸苔衛(wèi)星圖片表現(xiàn)出來的綠潮分布形狀將大為不同,很難辨別出昨天的綠潮在今天是在什么位置。6月底至7月青島至日照近海綠潮密集度比較高,分布特征也比較明顯。我們在青島至日照近海收集到3次時間間隔僅相差幾個小時的衛(wèi)星圖片,可以明顯地辨別出上一時刻在下一時刻的位置,作為實測滸苔漂移位置,與模式計算漂移位置作對比檢驗。三組觀測對比資料分別為2008年7月16日18時16分COSMO-1 SAR和19時4分COSMO-2 SAR,2011年7月11日SAR COSMO 06:02和MODIS 10:59,2014年7月16日RADARSAT-2衛(wèi)星6:53和18:08監(jiān)測數(shù)據(jù)。2014年7月16日RADARSAT-2衛(wèi)星于兩個時刻經(jīng)過該海域,分辨率分別為50m和25m (圖4)。RADARSAT-2成像基本不受云覆蓋的影響,對比兩張影像圖,在6:53的監(jiān)測結(jié)果中選取綠潮形狀變化不大的區(qū)域(藍色方框)提取有效信息作為綠潮漂移模式的初始場,預(yù)測其到18:08綠潮的位置,并與同時刻監(jiān)測結(jié)果作比較(圖5)。受青島近海風(fēng)和流的影響,綠潮向西北向漂移。三組綠潮的漂移軌跡預(yù)報結(jié)果檢驗如表1。綜合考慮綠潮有可能會沉降、不同衛(wèi)星的成像也會出現(xiàn)一些差別等因素,綠潮數(shù)值預(yù)報結(jié)果與監(jiān)測結(jié)果大體一致,因此黃海綠潮快速應(yīng)急模型能夠準(zhǔn)確快速預(yù)測綠潮斑塊的漂移軌跡和方向。

    圖4 2014年6月17日RADARSAT-2衛(wèi)星6:53(a)和18:08(b)綠潮(滸苔)解譯圖Fig.4 The Green Tide In terpretation of RADARSAT-2 at 6:53(a)and 18:08(b),on 17 June,2014

    表1 綠潮漂移軌跡預(yù)報結(jié)果檢驗Tab.1 The Forecasting Results Test of Green Tide

    圖5 2014年6月17日18:08綠潮監(jiān)測(綠點)和預(yù)測(紅點)結(jié)果比對Fig.5 The Observation(green dots)and Prediction(red dots)of Green Tide at 18:08 on June 17,2014

    1.3 黃海綠潮綜合業(yè)務(wù)平臺的建立

    基于綠潮應(yīng)急遙感監(jiān)測子系統(tǒng)和應(yīng)急漂移預(yù)測子系統(tǒng),建立一套集黃海綠潮災(zāi)害應(yīng)急遙感信息解譯、多源數(shù)據(jù)融合、漂移預(yù)測和預(yù)警產(chǎn)品制作和發(fā)布的綜合業(yè)務(wù)平臺。該平臺包含綠潮應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警兩個子平臺;兩者既獨立又相互聯(lián)系,其框架如圖6。綠潮遙感監(jiān)測子平臺的關(guān)鍵技術(shù)為基于K-均值聚類算法的綠潮散點抽稀、業(yè)務(wù)化專題圖制作;應(yīng)急預(yù)測子平臺的關(guān)鍵技術(shù)主要是綠潮數(shù)值模擬與GIS系統(tǒng)的無縫集成、模式自動調(diào)用和并行計算、不同類型的綠潮制圖輸出、業(yè)務(wù)化快報制作和發(fā)布。以2013年6月29日為例,基于綠潮遙感監(jiān)測子平臺準(zhǔn)自動化提取當(dāng)天綠潮衛(wèi)星和航空等遙感數(shù)據(jù)信息得到多源綜合解譯結(jié)果(圖7左圖),并制作綜合圖和綠潮遙感快報,可以進行快報發(fā)布或?qū)⒈O(jiān)測結(jié)果導(dǎo)入綠潮應(yīng)急預(yù)測子系統(tǒng)(圖7右圖),進行綠潮應(yīng)急漂移預(yù)測,根據(jù)需要選擇不同時刻、不同斑點的綠潮預(yù)測結(jié)果,制作綠潮趨勢預(yù)測圖(圖8錯誤!未找到引用源。左圖),并制作主要包括綠潮綜合分布圖、綠潮趨勢預(yù)測圖、綠潮斑點漂移圖、預(yù)警區(qū)域面積統(tǒng)計表、海洋環(huán)境預(yù)報表等信息的預(yù)警信息快報,最后發(fā)布預(yù)警信息(圖8右圖)。

    圖6 綠潮應(yīng)急遙感監(jiān)測(左圖)子平臺和預(yù)測預(yù)警子平臺(右圖)總體框架圖Fig.6 The Frame Diagram onSub-platform of Remote Sensing M onitoring(left figure)and Forecasting and W arning(right figure)on Green Tide.

    圖7 遙感監(jiān)測子平臺中監(jiān)測數(shù)據(jù)的準(zhǔn)自動化提取(左圖)和應(yīng)急預(yù)測子平臺中遙感解譯數(shù)據(jù)導(dǎo)入(右圖)Fig.7 The Quasi-autom atic Extraction ofM onitoring Data on Remote-Sensing M onitoring Sub-p latform(left figure)and Data Im port in the Forecasting and W arning Sub-p latform(righ t figure)

    圖8 綠潮漂移預(yù)測結(jié)果圖(左圖)、綠潮快報制作過程(中圖)和監(jiān)測預(yù)警產(chǎn)品輸出(右圖)Fig.8 The Figure on Forecast Result Exported from Platform(left figure),Express Production(m idd le figure)and Output Products(right figure)for Green Tide

    2 討論

    黃海綠潮應(yīng)急遙感監(jiān)測和預(yù)測預(yù)警系統(tǒng)從2008年開始研發(fā)到2010年底基本建設(shè)完成,并在應(yīng)用中不斷完善升級,大大提高了工作效率,進一步提升了黃海綠潮災(zāi)害應(yīng)急工作預(yù)警能力。本系統(tǒng)服務(wù)于國家海洋局和山東省海洋與漁業(yè)廳、青島市海洋與漁業(yè)局、威海市海洋與漁業(yè)局、煙臺市海洋與漁業(yè)局、日照市海洋與漁業(yè)局等政府防災(zāi)減災(zāi)應(yīng)急部門,為海上運動、水產(chǎn)養(yǎng)殖、濱海旅游、交通運輸、漁業(yè)生產(chǎn)等涉?;顒雍秃Q蠊芾硖峁┛茖W(xué)可靠的技術(shù)保障,特別是在2008年奧帆船賽和2012年亞沙會的綠潮預(yù)測保障工作中,發(fā)揮了重要的作用,并取得巨大的社會、經(jīng)濟和生態(tài)效益。

    海洋中滸苔條帶的寬度通常不大于10 m,而分辨率為25~30 m的SAR圖像中異常條帶的寬度是多個像素點。基于遙感圖像異常對綠潮量的估算可能會偏高1個量級(10倍)左右,所以單獨利用海洋遙感圖像判定綠潮在海面的覆蓋率存在一定問題。所以迫切需要完善海洋綠潮監(jiān)視監(jiān)測系統(tǒng),衛(wèi)星遙感需緊密結(jié)合航空遙感、陸岸巡視等其他監(jiān)測結(jié)果,并實現(xiàn)多源數(shù)據(jù)充分融合,提高系統(tǒng)的監(jiān)測數(shù)據(jù)提取精度。在多源綠潮融合結(jié)果中,不同影像應(yīng)用基于NDVI與滸苔覆蓋率關(guān)系的混合像元分解方法,提取精度基本達到HJ-1 CCD提取精度,其余SAR影像分辨率各不相同,未將結(jié)果統(tǒng)一到HJ-1提取精度,因此為保證融合結(jié)果各星源信息提取精度統(tǒng)一,需對該問題開展進一步的研究。綠潮的發(fā)生發(fā)展過程十分復(fù)雜,在漂移的過程中會有生長、死亡沉降等生態(tài)現(xiàn)象。我們所采用的方法是在不考慮綠潮自身

    生態(tài)特征的情況下,綠潮隨風(fēng)和流等要素漂移的物理過程。特別在6、7月份綠潮暴發(fā)期間,綠潮的生長速度較快,單純的物理模型不能模擬綠潮藻類的生長受溫度、光照、pH、鹽度和營養(yǎng)鹽等諸多環(huán)境因素影響的過程,所以下一步我們將建立綠潮物理—生態(tài)模型進行漂移、分布和生長預(yù)測,從而提高黃海綠潮應(yīng)急監(jiān)測和預(yù)測預(yù)警的能力,為政府相關(guān)部門提供更準(zhǔn)確和豐富的綠潮信息。

    [1] Gower J,King S,Goncalves P.Globalmonitoring of plankton blooms using MERISMCI,29(21),2008,6209-6216.

    [2] 張娟.滸苔遙感監(jiān)測方法研究及軟件實現(xiàn)—以青島奧帆賽場及周邊海域為例[D].成都:電子科技大學(xué),2009.

    [3] 李三妹,李亞君,董海鷹,等.淺析衛(wèi)星遙感在黃海滸苔監(jiān)測中的應(yīng)用[J].應(yīng)用氣象學(xué)報,2010,21(1):76-82.

    [4] Elijah R,Amina R,Mads S T,et al.Spectral definition of the Macro-algae Ulva curvata in the Back-barrier Bays of the eastern shore of Virginia,USA[J].Int.J.Remote Sens.,2012,33(2): 586-603.

    [5] 李穎,梁剛,于水明,等.監(jiān)測滸苔災(zāi)害的微波遙感數(shù)據(jù)選取[J].海洋環(huán)境科學(xué),2011,30(5): 739-742.

    [6] 蔣興偉,鄒亞榮,王華,等.基于SAR快速提取滸苔信息應(yīng)用研究[J].海洋學(xué)報,2009,31(2): 63-68.

    [7] Cui T W,Zhang J,Sun L E,et al.Satellite monitoring of massive green macroalgae bloom (GMB):Imaging ability comparison ofmulti-source data and drifting velocity estimation[J].Int.J. Remote Sens.,2012,33(17):5513-5527.

    [8] Hu C M,He M X.Origin and offshore extent of floating algae in olympic sailing area[J].Eos,American Geophysical Union Transactions,2008,89 (33):302-303.

    [9] Hu C M,Li D Q,Chen C S,et al.On the recurrent Ulvaprolifera blooms in the Yellow Sea and East China Sea[J].Journal of Geophysical Research,2010 (115):C05017.doi:10.1029/2009JC005561.

    [10] Keesing J K,Liu D Y,F(xiàn)earns P,et al.Inter-and Intra-annual patterns of Ulvaprolifera green tides in the Yellow Sea during 2007-2009,their origin and relationship to the expansion of coastal seaweed aquaculture in China[J].Marine Pollution Bulletin,2011,62(6):1169-1182.

    [11] Shi W,Wang M.Green macroalgae blooms in the Yellow Sea during the spring and summer of2008[J] J.Geophys.R.,2009(114):C120010,doi: 10.1029/2009JC005513.

    [12] 高松,黃娟,白濤,等,2008年與2009年黃海綠潮漂移路徑分析[J].海洋科學(xué),38(2):94-98.

    [13] 王寧,曹叢華,黃娟,等.基于遙感監(jiān)測的黃海綠潮漂移路徑及分布面積特征分析[J].防災(zāi)科技學(xué)院學(xué)報,15(4):24-29.

    [14] 黃娟,吳玲娟,高松,等,黃海綠潮分布年際變化分析[J].激光生物學(xué)報,2014.

    [15] AUROUSSEAU P.Les flux d'azote et de phosphore provenantdes bassins versants de la rade de Brest [J].Comparaison avec laBretagne.Océanis,2001,27:137-161.

    [16] CUGIER P,LE HIR P.Three dimensional(3D) ecological modellingof the Bay of Seine(English Channel,F(xiàn)rance)[J].JSea Res,2005,54:104-124.

    [17] CUGIER P.P.LE HIR.Development of a 3D hydrodynamicalmodel for coastal ecosystem modeling,Application to the plume of the Seine River(France) [J].Estuar Coast Shelf Sci,2000,55:673-695.

    [18] Ménesguen,JC Salomon.Eutrophication modelling as a toolfor fightingagainst Ulva coastalmass blooms,B.A.Schrefler and O.C.Zienkiewicz[eds.] [M].Computer modelling in oceanengineering. Balkema,1988,443-450.

    [19] Perrot T,Dion P,Populus J,Guillaumont B,Brault D.A predictiveapproach to the mapping of fucoid beds using Spotmagery and a digital elevation model [C].VIIIth International SeaweedSymposium,June 20-25,2004,Bergen,Norway.

    [20] 喬方利.2008與2010年黃海滸苔漂移輸運特征對比[J].科學(xué)通報,2011,56(18):1470-1476.

    [21] 李燕,李云,劉欽政.滸苔漂移軌跡預(yù)報系統(tǒng)[J].海洋預(yù)報,2011,56(18):1470-1476.

    [22] 黃娟,吳玲娟,高松,等.黃海綠潮應(yīng)急漂移數(shù)值模擬[J],海洋預(yù)報,2011,28(1):25-32.

    [23] 吳玲娟,曹叢華,黃娟,等.黃海綠潮應(yīng)急溯源數(shù)值模擬[J],海洋科學(xué),2011,35(6):44-47.

    [24] Xu Jiangling, Huang Juan, Gao Song, et al.,Assimilation of High Frequency Radar Data into a Shelf Sea Circulation Model[J],Journal of Ocean

    University of China,2014,13(4):572-578.

    [25] 黃娟,徐江玲,高松,等.基于海上試驗對海上漂移物運移軌跡影響因素分析[J].海洋預(yù)報,2014,31(4):97-104.

    Remote-sensing M onitoring,F(xiàn)orecasting and Warning System for Emergency Response on Green Tide in Yellow Sea

    Wu lingJuan1,2,Gao Song1,2,Ding Yi1,2,Cao Yajing1,2,Jin Xifang1,2
    (1.North China Sea Marine Forecast Center ofState Oceanic Administration,Qingdao,Shandong 266061,China;2.Shandong Provincial Key Laboratory ofMarine Ecological Environmentand Disaster Prevention and Mitigation,Qingdao,Shandong 266061,China)

    The system on remote-sensingmonitoring,forecasting and warning for emergency response on green tide in Yellow Sea is developed based on the study of monitoring and forecasting techniques.The system could provide remote-sensingmonitoring,information extraction,data collection and analysis,and broadcast of forecast and warning production.On remote-sensing monitoring,the technique on quasi-automatic information extraction and fusion of multi-source,multi-window data is developed;the fast drifting forecastmodel is developed based on fine atmospheric and oceanic field.The system serves the State Oceanic Adm inistration and related government departments,and provides scientific,reliable and technical support for marine sports,aquaculture,tourism,transportation,fishery,and marinemanagement.The system has played important roles in 2008 Olympic Sailing Games and 2012 Asian Beach Games,and has achieved great social,economic and ecological benefits.

    green tide;remote-sensingmonitoring;forecasting and warning;system

    X

    :B

    :1673-8047(2015)01-0059-09

    2014-10-20

    山東省海洋生態(tài)環(huán)境與防災(zāi)減災(zāi)重點實驗室開放基金(201407);國家海洋局公益性項目(201205010)

    :吳玲娟(1979—),女,博士,高工,主要從事海洋環(huán)境要素預(yù)報和研究和防災(zāi)減災(zāi)工作。

    猜你喜歡
    綠潮黃海應(yīng)急
    多維深入復(fù)盤 促進應(yīng)急搶險
    HY-1衛(wèi)星CZI影像在黃海綠潮監(jiān)測的應(yīng)用
    綠潮監(jiān)測行業(yè)標(biāo)準(zhǔn)通過專家審查
    東方濕地 黃海明珠
    華人時刊(2020年21期)2021-01-14 01:32:28
    完善應(yīng)急指揮機制融嵌應(yīng)急準(zhǔn)備、響應(yīng)、處置全周期
    黃海簡介
    應(yīng)急管理部6個“怎么看”
    勞動保護(2018年5期)2018-06-05 02:12:05
    基于MODIS數(shù)據(jù)的2016年黃海綠潮災(zāi)害動態(tài)監(jiān)測研究
    三角恒等變換去哪兒了
    國際新應(yīng)急標(biāo)準(zhǔn)《核或輻射應(yīng)急的準(zhǔn)備與響應(yīng)》的釋疑
    极品教师在线视频| 亚洲精品一区蜜桃| 六月丁香七月| 高清不卡的av网站| 2018国产大陆天天弄谢| 午夜日本视频在线| av福利片在线观看| 国产午夜精品一二区理论片| 涩涩av久久男人的天堂| 妹子高潮喷水视频| 精品一区二区三区视频在线| 黑人猛操日本美女一级片| 久久久成人免费电影| 日本av免费视频播放| 性高湖久久久久久久久免费观看| 涩涩av久久男人的天堂| 日本黄色片子视频| 欧美一区二区亚洲| xxx大片免费视频| 亚洲一级一片aⅴ在线观看| 亚洲熟女精品中文字幕| av网站免费在线观看视频| 天堂中文最新版在线下载| 精品酒店卫生间| 欧美变态另类bdsm刘玥| 亚洲精品一区蜜桃| 欧美极品一区二区三区四区| 久久久久精品性色| 亚洲va在线va天堂va国产| 一级毛片久久久久久久久女| 五月玫瑰六月丁香| 黄片wwwwww| 亚洲成人中文字幕在线播放| 黄片无遮挡物在线观看| av网站免费在线观看视频| 欧美三级亚洲精品| 国产av码专区亚洲av| 午夜免费观看性视频| 丝袜喷水一区| 韩国高清视频一区二区三区| 九色成人免费人妻av| 国产熟女欧美一区二区| 丰满少妇做爰视频| 你懂的网址亚洲精品在线观看| 国产久久久一区二区三区| 国产爱豆传媒在线观看| 人妻少妇偷人精品九色| 国产乱人视频| 久久精品国产a三级三级三级| 日韩三级伦理在线观看| 精品少妇黑人巨大在线播放| 午夜老司机福利剧场| 观看av在线不卡| av免费观看日本| 午夜福利网站1000一区二区三区| 十八禁网站网址无遮挡 | 精品少妇久久久久久888优播| 丰满乱子伦码专区| 青青草视频在线视频观看| 2018国产大陆天天弄谢| 高清日韩中文字幕在线| 久久久国产一区二区| 伦理电影免费视频| 97超视频在线观看视频| 日韩不卡一区二区三区视频在线| 男女下面进入的视频免费午夜| 在线看a的网站| 亚洲欧美日韩东京热| 日日摸夜夜添夜夜添av毛片| 国产欧美日韩精品一区二区| 又黄又爽又刺激的免费视频.| 婷婷色av中文字幕| 热99国产精品久久久久久7| 国产成人a∨麻豆精品| 成人18禁高潮啪啪吃奶动态图 | 国产成人精品婷婷| 国产精品成人在线| 国产乱人视频| 免费观看av网站的网址| 国产精品一区二区在线观看99| 蜜桃在线观看..| 欧美97在线视频| 老女人水多毛片| 免费黄网站久久成人精品| 高清毛片免费看| videos熟女内射| 直男gayav资源| 中文字幕亚洲精品专区| 亚洲不卡免费看| 美女视频免费永久观看网站| 免费黄色在线免费观看| 51国产日韩欧美| 黄色视频在线播放观看不卡| 在线免费十八禁| 婷婷色综合www| 97热精品久久久久久| 美女内射精品一级片tv| 日日摸夜夜添夜夜爱| 午夜福利视频精品| 国产精品嫩草影院av在线观看| 在现免费观看毛片| av不卡在线播放| 国产精品国产三级国产专区5o| 国产日韩欧美亚洲二区| 国内精品宾馆在线| 在现免费观看毛片| av在线播放精品| 亚洲欧美日韩卡通动漫| 男男h啪啪无遮挡| 久久 成人 亚洲| 毛片一级片免费看久久久久| av国产精品久久久久影院| 一区二区三区免费毛片| 夫妻午夜视频| 九色成人免费人妻av| 黑人高潮一二区| 亚洲国产精品专区欧美| 少妇人妻 视频| 久久久午夜欧美精品| 天堂中文最新版在线下载| 91午夜精品亚洲一区二区三区| 午夜福利网站1000一区二区三区| 国产精品一区二区性色av| 免费av不卡在线播放| 精品少妇黑人巨大在线播放| 国产毛片在线视频| 久久av网站| 欧美人与善性xxx| 中文资源天堂在线| kizo精华| 蜜桃在线观看..| 亚洲人成网站在线观看播放| 国产av一区二区精品久久 | 日日摸夜夜添夜夜添av毛片| 亚洲精品一二三| 国产精品人妻久久久久久| 秋霞伦理黄片| av视频免费观看在线观看| 老师上课跳d突然被开到最大视频| 成年av动漫网址| 中文字幕亚洲精品专区| 中文字幕制服av| 特大巨黑吊av在线直播| 中文字幕免费在线视频6| 老师上课跳d突然被开到最大视频| 国产在线男女| 一级毛片aaaaaa免费看小| 80岁老熟妇乱子伦牲交| 人人妻人人看人人澡| 成人高潮视频无遮挡免费网站| 国产探花极品一区二区| 亚洲经典国产精华液单| 亚洲国产高清在线一区二区三| 一级毛片aaaaaa免费看小| 六月丁香七月| 美女内射精品一级片tv| 美女内射精品一级片tv| av线在线观看网站| 97热精品久久久久久| 国产精品熟女久久久久浪| 秋霞在线观看毛片| 99热这里只有精品一区| 欧美成人午夜免费资源| 成人综合一区亚洲| 麻豆成人av视频| 久久久久久久久久久免费av| 国产欧美亚洲国产| 日韩欧美一区视频在线观看 | 女人十人毛片免费观看3o分钟| 国产亚洲欧美精品永久| 日韩av免费高清视频| 日韩av免费高清视频| 丰满乱子伦码专区| 人人妻人人添人人爽欧美一区卜 | 女性被躁到高潮视频| 熟女人妻精品中文字幕| 久久久久性生活片| 国产精品蜜桃在线观看| 欧美高清性xxxxhd video| 极品少妇高潮喷水抽搐| 夜夜爽夜夜爽视频| 亚洲人成网站在线观看播放| 国产探花极品一区二区| 黑丝袜美女国产一区| 欧美xxⅹ黑人| 22中文网久久字幕| 各种免费的搞黄视频| 一级毛片电影观看| 免费观看无遮挡的男女| 成年人午夜在线观看视频| 美女高潮的动态| 久久久久久久大尺度免费视频| 久久精品国产a三级三级三级| 久久久久性生活片| 人妻系列 视频| 国产爽快片一区二区三区| 亚洲人成网站高清观看| 色5月婷婷丁香| 国产精品一区二区在线观看99| 插逼视频在线观看| 在线观看av片永久免费下载| 亚洲一区二区三区欧美精品| 深夜a级毛片| 一边亲一边摸免费视频| 青春草亚洲视频在线观看| 欧美日韩视频高清一区二区三区二| 视频中文字幕在线观看| 免费观看的影片在线观看| 一区二区三区免费毛片| 嫩草影院入口| 亚洲va在线va天堂va国产| 我要看日韩黄色一级片| 久久影院123| 91aial.com中文字幕在线观看| 一本—道久久a久久精品蜜桃钙片| 深爱激情五月婷婷| 国产黄片视频在线免费观看| 一区二区三区免费毛片| 免费黄色在线免费观看| av线在线观看网站| 美女中出高潮动态图| 国产毛片在线视频| 久久国产精品大桥未久av | 欧美日本视频| 色网站视频免费| 观看av在线不卡| 国产亚洲91精品色在线| 久久久久人妻精品一区果冻| 日本黄色片子视频| 日韩国内少妇激情av| 亚洲av男天堂| 亚洲av免费高清在线观看| 特大巨黑吊av在线直播| 国产精品国产av在线观看| 一级毛片黄色毛片免费观看视频| 国产午夜精品久久久久久一区二区三区| 亚洲av.av天堂| 综合色丁香网| 国产 精品1| 欧美高清性xxxxhd video| 国产av一区二区精品久久 | 黄色日韩在线| 亚洲欧洲国产日韩| 亚洲国产欧美人成| 久久 成人 亚洲| 免费在线观看成人毛片| 蜜桃久久精品国产亚洲av| 内射极品少妇av片p| 在线观看av片永久免费下载| 99久久精品热视频| 国产免费视频播放在线视频| 亚洲欧美日韩另类电影网站 | 蜜桃久久精品国产亚洲av| 亚洲av综合色区一区| 精品一品国产午夜福利视频| 亚洲成色77777| 欧美精品一区二区免费开放| 99九九线精品视频在线观看视频| 亚洲精品第二区| 黄色欧美视频在线观看| 一级爰片在线观看| 尾随美女入室| 亚洲欧美精品自产自拍| 亚洲欧美日韩无卡精品| 一级片'在线观看视频| 在线观看免费高清a一片| 久久精品夜色国产| 老熟女久久久| 少妇的逼水好多| 国产熟女欧美一区二区| 一区在线观看完整版| 亚洲内射少妇av| 亚洲国产最新在线播放| 日本黄色片子视频| 少妇精品久久久久久久| 日韩av不卡免费在线播放| 99久久精品热视频| 色视频在线一区二区三区| 18禁裸乳无遮挡免费网站照片| 久久99热6这里只有精品| tube8黄色片| 国产精品.久久久| 亚洲av欧美aⅴ国产| 少妇人妻久久综合中文| 久久97久久精品| 狂野欧美白嫩少妇大欣赏| 欧美一级a爱片免费观看看| .国产精品久久| 欧美日韩在线观看h| 大片电影免费在线观看免费| 亚洲婷婷狠狠爱综合网| 成年免费大片在线观看| 直男gayav资源| 免费看光身美女| 久久热精品热| 高清日韩中文字幕在线| 亚洲精品久久午夜乱码| 国产美女午夜福利| 国产高清三级在线| 亚洲中文av在线| 2022亚洲国产成人精品| 国产一区亚洲一区在线观看| 日日摸夜夜添夜夜爱| 一级黄片播放器| 黄色欧美视频在线观看| 女人久久www免费人成看片| freevideosex欧美| 精品午夜福利在线看| 国产高清有码在线观看视频| 亚洲欧美日韩另类电影网站 | 亚洲国产精品999| 国产女主播在线喷水免费视频网站| 日韩精品有码人妻一区| 在线观看三级黄色| 精品久久久久久久末码| 婷婷色av中文字幕| 国产 一区 欧美 日韩| 色综合色国产| 亚洲美女黄色视频免费看| 久久久精品94久久精品| 国产av国产精品国产| 日韩一区二区三区影片| 大又大粗又爽又黄少妇毛片口| 日本黄色日本黄色录像| 亚洲av欧美aⅴ国产| 最近最新中文字幕免费大全7| 啦啦啦视频在线资源免费观看| 麻豆精品久久久久久蜜桃| 成年女人在线观看亚洲视频| 色网站视频免费| 久久久久精品久久久久真实原创| 最近中文字幕高清免费大全6| 亚洲成人一二三区av| 80岁老熟妇乱子伦牲交| 国产高潮美女av| 中国三级夫妇交换| 亚洲国产最新在线播放| 国产精品久久久久久久久免| 日日摸夜夜添夜夜添av毛片| 在线观看av片永久免费下载| 国产亚洲91精品色在线| 一本久久精品| 日日啪夜夜爽| 网址你懂的国产日韩在线| 深夜a级毛片| 久久久久人妻精品一区果冻| 久久人人爽人人片av| 日韩一区二区视频免费看| 纵有疾风起免费观看全集完整版| 男女下面进入的视频免费午夜| 日韩成人av中文字幕在线观看| 啦啦啦啦在线视频资源| 国产成人午夜福利电影在线观看| 男的添女的下面高潮视频| 一本色道久久久久久精品综合| 一二三四中文在线观看免费高清| 国产69精品久久久久777片| 熟女av电影| 久久久午夜欧美精品| 一级二级三级毛片免费看| 国产精品熟女久久久久浪| av黄色大香蕉| 五月开心婷婷网| 国产女主播在线喷水免费视频网站| 国产精品人妻久久久影院| 女人十人毛片免费观看3o分钟| 九九爱精品视频在线观看| 国产在线视频一区二区| 精品人妻偷拍中文字幕| 国产成人精品婷婷| 精品久久久精品久久久| 午夜精品国产一区二区电影| 不卡视频在线观看欧美| 观看av在线不卡| 免费观看无遮挡的男女| 欧美高清性xxxxhd video| 亚洲av二区三区四区| 日韩制服骚丝袜av| 午夜福利影视在线免费观看| 亚洲人成网站高清观看| 亚洲av日韩在线播放| 久久人人爽人人爽人人片va| 亚洲欧美中文字幕日韩二区| 午夜视频国产福利| 国产精品福利在线免费观看| 人人妻人人添人人爽欧美一区卜 | 中文字幕免费在线视频6| 免费av不卡在线播放| 中国三级夫妇交换| av女优亚洲男人天堂| 午夜免费男女啪啪视频观看| 国内少妇人妻偷人精品xxx网站| 日韩三级伦理在线观看| 观看美女的网站| 少妇人妻精品综合一区二区| 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 欧美xxxx性猛交bbbb| 亚洲人成网站高清观看| 日韩av在线免费看完整版不卡| 舔av片在线| 亚洲成人一二三区av| 亚洲经典国产精华液单| 日韩中文字幕视频在线看片 | 人人妻人人澡人人爽人人夜夜| 色婷婷av一区二区三区视频| 久久久久网色| 日韩成人伦理影院| 成人黄色视频免费在线看| 91精品伊人久久大香线蕉| 尤物成人国产欧美一区二区三区| 亚洲第一区二区三区不卡| 国产精品三级大全| 免费在线观看成人毛片| 免费看光身美女| 中国三级夫妇交换| 久久久久精品性色| 人人妻人人爽人人添夜夜欢视频 | 青春草视频在线免费观看| av在线观看视频网站免费| 美女xxoo啪啪120秒动态图| 欧美 日韩 精品 国产| av天堂中文字幕网| 国产成人精品婷婷| 亚洲精品第二区| 国产在线视频一区二区| www.色视频.com| av女优亚洲男人天堂| 在线播放无遮挡| 亚洲熟女精品中文字幕| 97热精品久久久久久| 免费人成在线观看视频色| 免费人妻精品一区二区三区视频| 国产高清有码在线观看视频| 精品亚洲成a人片在线观看 | 一本久久精品| 日韩视频在线欧美| 亚洲av国产av综合av卡| 偷拍熟女少妇极品色| 中文字幕制服av| 寂寞人妻少妇视频99o| 天美传媒精品一区二区| 精品亚洲成a人片在线观看 | 日本一二三区视频观看| 久久毛片免费看一区二区三区| 亚洲熟女精品中文字幕| 性色avwww在线观看| 小蜜桃在线观看免费完整版高清| 亚洲成人中文字幕在线播放| 国产精品久久久久成人av| kizo精华| 亚洲av男天堂| 免费黄色在线免费观看| 欧美最新免费一区二区三区| 高清av免费在线| 欧美日韩亚洲高清精品| 天堂中文最新版在线下载| 美女xxoo啪啪120秒动态图| 国产日韩欧美亚洲二区| 亚洲精品自拍成人| 亚洲婷婷狠狠爱综合网| 国产成人aa在线观看| 久久久久国产网址| 精品久久久久久久久av| 久久久久久久大尺度免费视频| 欧美xxⅹ黑人| av国产精品久久久久影院| a级一级毛片免费在线观看| 最新中文字幕久久久久| 老司机影院成人| 看免费成人av毛片| 精品久久国产蜜桃| 亚洲美女搞黄在线观看| 日韩一本色道免费dvd| 一二三四中文在线观看免费高清| av专区在线播放| 国产精品秋霞免费鲁丝片| 久久精品夜色国产| 欧美日韩综合久久久久久| av一本久久久久| 内地一区二区视频在线| 男人舔奶头视频| 午夜老司机福利剧场| 麻豆成人av视频| 大码成人一级视频| 久久精品国产鲁丝片午夜精品| 天天躁夜夜躁狠狠久久av| 国产精品嫩草影院av在线观看| 男人添女人高潮全过程视频| 99热网站在线观看| 国产一级毛片在线| 男人舔奶头视频| av卡一久久| 日韩成人伦理影院| 国产欧美另类精品又又久久亚洲欧美| 蜜桃亚洲精品一区二区三区| 男女免费视频国产| 国产黄片视频在线免费观看| h视频一区二区三区| 少妇人妻 视频| 日本黄色片子视频| 天堂中文最新版在线下载| 视频中文字幕在线观看| 午夜精品国产一区二区电影| 男人舔奶头视频| 国产亚洲精品久久久com| 日韩,欧美,国产一区二区三区| 一级毛片电影观看| 极品少妇高潮喷水抽搐| 美女主播在线视频| 99热全是精品| 成人毛片a级毛片在线播放| 六月丁香七月| 亚洲在久久综合| 国产亚洲最大av| 男女无遮挡免费网站观看| 色吧在线观看| 亚洲av不卡在线观看| 99精国产麻豆久久婷婷| 国精品久久久久久国模美| 蜜桃在线观看..| 爱豆传媒免费全集在线观看| 成年av动漫网址| 日日啪夜夜撸| 免费观看的影片在线观看| 国产欧美亚洲国产| 久久久精品94久久精品| 国产大屁股一区二区在线视频| 亚洲成人中文字幕在线播放| 亚洲人与动物交配视频| 国产精品久久久久久久电影| 久久精品久久精品一区二区三区| 日韩欧美一区视频在线观看 | 亚洲激情五月婷婷啪啪| 性色avwww在线观看| 两个人的视频大全免费| 欧美成人一区二区免费高清观看| 99国产精品免费福利视频| 激情五月婷婷亚洲| 亚洲国产毛片av蜜桃av| 亚洲国产av新网站| 国产精品蜜桃在线观看| 自拍偷自拍亚洲精品老妇| 精品人妻熟女av久视频| 国产熟女欧美一区二区| 美女福利国产在线 | 日韩一本色道免费dvd| 99久国产av精品国产电影| 精品国产露脸久久av麻豆| 美女福利国产在线 | 日本av免费视频播放| 免费在线观看成人毛片| 日韩av在线免费看完整版不卡| 成年人午夜在线观看视频| 我的老师免费观看完整版| 成人无遮挡网站| .国产精品久久| 亚洲欧美清纯卡通| 1000部很黄的大片| 少妇熟女欧美另类| 中文欧美无线码| 欧美3d第一页| 亚洲精品日韩av片在线观看| 91精品国产九色| 亚洲怡红院男人天堂| 中文字幕人妻熟人妻熟丝袜美| 十分钟在线观看高清视频www | 亚洲四区av| 成年美女黄网站色视频大全免费 | 日本黄大片高清| 亚洲精品日本国产第一区| 日本黄大片高清| 少妇人妻精品综合一区二区| 国产精品成人在线| 极品教师在线视频| 亚洲怡红院男人天堂| 97在线人人人人妻| 亚洲一级一片aⅴ在线观看| 成人综合一区亚洲| 国产精品三级大全| a级一级毛片免费在线观看| 丰满迷人的少妇在线观看| 一个人看视频在线观看www免费| 一级毛片久久久久久久久女| 熟女av电影| 色视频在线一区二区三区| 熟女av电影| 在线观看美女被高潮喷水网站| 91狼人影院| 午夜免费鲁丝| 久久久久国产网址| 一区二区三区免费毛片| 亚洲第一区二区三区不卡| 日本欧美国产在线视频| 日韩欧美 国产精品| 亚洲欧美中文字幕日韩二区| 亚洲av国产av综合av卡| 久久人人爽人人片av| 国产黄色免费在线视频| 久久精品夜色国产| 国产黄色免费在线视频| 秋霞在线观看毛片| 肉色欧美久久久久久久蜜桃| 国产免费一区二区三区四区乱码| 十八禁网站网址无遮挡 | av又黄又爽大尺度在线免费看| 国产精品不卡视频一区二区| 国产免费一区二区三区四区乱码| 高清黄色对白视频在线免费看 | 亚洲国产精品国产精品| 丰满人妻一区二区三区视频av| 80岁老熟妇乱子伦牲交| 亚洲欧美日韩卡通动漫| 天堂中文最新版在线下载| 日本黄色片子视频|