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

    珠江河口潮周期泥沙動(dòng)態(tài)輸移特征遙感反演分析

    2022-05-26 09:31:48潘洪洲何穎清馮佑斌熊龍海吳俊涌
    熱帶地理 2022年5期
    關(guān)鍵詞:懸沙河口區(qū)河口

    潘洪洲,何穎清,馮佑斌,熊龍海,吳俊涌,楊 敏

    (a. 珠江河口動(dòng)力學(xué)及伴生過(guò)程調(diào)控重點(diǎn)實(shí)驗(yàn)室;b. 珠江水利委員會(huì)珠江水利科學(xué)研究院,廣州 510611)

    粵港澳大灣區(qū),核心水域即為珠江河口區(qū)域(以下簡(jiǎn)稱“河口區(qū)”),戰(zhàn)略地位重要。區(qū)內(nèi)河網(wǎng)密布交錯(cuò),形成“三江匯流,八門(mén)出海”的網(wǎng)河特色,河道、河口區(qū)徑流與潮流相互作用,水情復(fù)雜。受人類(lèi)活動(dòng)和海洋潮流共同影響,河口區(qū)水沙動(dòng)力條件、河道演變發(fā)生較大變化,嚴(yán)重威脅粵港澳大灣區(qū)生態(tài)環(huán)境及防洪、供水及航運(yùn)安全。區(qū)內(nèi)未來(lái)人類(lèi)經(jīng)濟(jì)活動(dòng)及相應(yīng)規(guī)劃建設(shè)的節(jié)奏將繼續(xù)加快,河網(wǎng)水沙調(diào)控、咸潮入侵、水質(zhì)安全等水問(wèn)題將更加突出,亟需加強(qiáng)河口區(qū)水沙動(dòng)力環(huán)境監(jiān)測(cè)。水體泥沙是河道演變的物質(zhì)來(lái)源,其淤積過(guò)程直接關(guān)系到河道演變,對(duì)防洪、供水、航道、海岸工程等安全都會(huì)造成直接影響。水體泥沙作為水色要素重要指標(biāo)之一,能夠直接影響水體的光學(xué)特征,如增加水體渾濁度、降低透明度等,進(jìn)而影響太陽(yáng)光在水體中的傳播,影響水生生物的光合作用、初級(jí)生產(chǎn)力,對(duì)沿岸生態(tài)環(huán)境、生物資源、漁業(yè)資源產(chǎn)生重要影響(姜麗君,2018;潘磊劍,2020)。

    由于河口區(qū)懸沙受徑流、潮汐、鹽淡混合及風(fēng)浪作用等多種因素的影響,懸浮泥沙的時(shí)空分布、運(yùn)動(dòng)特性和沉積部位經(jīng)常發(fā)生變化,對(duì)于懸沙的監(jiān)測(cè)具有很強(qiáng)的實(shí)時(shí)性要求,珠江河口區(qū)域無(wú)常規(guī)高頻次遙感衛(wèi)星覆蓋,中、高分遙感衛(wèi)星受重訪周期的限制,以及天氣的影響,使得泥沙質(zhì)量濃度的遙感反演在時(shí)間上不具有連續(xù)性,實(shí)時(shí)性較差,難以實(shí)時(shí)掌握河口泥沙質(zhì)量濃度的時(shí)空分布與動(dòng)態(tài)變化規(guī)律。近年來(lái),常規(guī)泥沙監(jiān)測(cè)手段不斷豐富,如浮標(biāo)、平臺(tái)監(jiān)測(cè)系統(tǒng)或水文監(jiān)測(cè)站點(diǎn),可以對(duì)泥沙質(zhì)量濃度等水質(zhì)參數(shù)進(jìn)行實(shí)時(shí)在線監(jiān)測(cè)。但通過(guò)常規(guī)站點(diǎn)監(jiān)測(cè)獲取的點(diǎn)狀數(shù)據(jù)信息,僅能代表觀測(cè)區(qū)域小范圍內(nèi)的泥沙狀況,且無(wú)法獲取大范圍水體表觀光學(xué)信息。因而導(dǎo)致河口區(qū)潮周期內(nèi)泥沙輸移遙感監(jiān)測(cè)研究較少,無(wú)法掌握河口區(qū)復(fù)雜水情下的泥沙輸移特征。2016年中國(guó)首顆地球同步軌道高分對(duì)地觀測(cè)衛(wèi)星高分四號(hào)(GF-4)投入使用,其具有高時(shí)間和高空間分辨率的觀測(cè)優(yōu)勢(shì),在水色遙感上具有重大應(yīng)用潛力(沈宏,2012;李玨,2016;劉明等,2016;陳曉英等,2018;宋挺等,2018)。其多時(shí)相遙感數(shù)據(jù)采用連續(xù)觀測(cè),能夠較好地反映水域高頻變化的特點(diǎn)(Lamquin et al., 2012;江彬彬等,2015;伍高燕,2015;劉波等,2016;楊雪飛,2016;孫璐等,2017;胡耀躲等,2018)。因此,為研究珠江河口泥沙輸移特征,本文以GF-4衛(wèi)星為數(shù)據(jù)源,結(jié)合野外觀測(cè)實(shí)驗(yàn),構(gòu)建反演模型,開(kāi)展河口區(qū)泥沙質(zhì)量濃度高頻次反演研究。以期探究潮周期內(nèi)珠江河口泥沙動(dòng)態(tài)輸移特征,為河口沉積物輸運(yùn)、泥沙來(lái)源及泥沙數(shù)值模擬研究等提供數(shù)據(jù)支撐。

    1 數(shù)據(jù)與方法

    1.1 研究區(qū)概況

    珠江河口位于廣東省中南部,南臨南海,三角洲河網(wǎng)和河口灣并存,承擔(dān)安全宣泄洪水,保障供水、航運(yùn)、生態(tài)安全等重要任務(wù),是珠江流域防洪、供水、水生態(tài)安全保障體系的重要組成部分。珠江河口是世界上水系結(jié)構(gòu)、動(dòng)力特性最復(fù)雜,人類(lèi)活動(dòng)干擾最顯著的河口之一。主要有以下特點(diǎn):1)河網(wǎng)交錯(cuò),水系結(jié)構(gòu)復(fù)雜,珠江水系西江、北江、東江三江匯流,由八大口門(mén)出海,河網(wǎng)密度高,屬?gòu)?fù)雜河口;2)多級(jí)分流,徑潮動(dòng)力復(fù)雜,珠三角河網(wǎng)的水沙經(jīng)多級(jí)分汊節(jié)點(diǎn)重新分配,潮優(yōu)型與河優(yōu)型河口相互依存;3)灘涂發(fā)育,人類(lèi)活動(dòng)劇烈,灘涂圍墾利用面積大,珠江河口人類(lèi)活動(dòng)對(duì)灘槽的影響已遠(yuǎn)超自然演變;4)通江達(dá)海,航運(yùn)交通發(fā)達(dá),河口區(qū)規(guī)劃分布有多個(gè)主要航道;5)物種豐富,生態(tài)系統(tǒng)多樣。珠江河口極具典型性,水沙輸移環(huán)境復(fù)雜,具有重要的研究?jī)r(jià)值。

    1.2 野外觀測(cè)數(shù)據(jù)獲取

    于2019-12-28和2020-06-20分別開(kāi)展珠江河口枯季和洪季野外定點(diǎn)及走航觀測(cè)實(shí)驗(yàn)。定點(diǎn)觀測(cè)主要通過(guò)珠江河口原型觀測(cè)浮標(biāo)站點(diǎn)自動(dòng)觀測(cè)及船上人工采樣觀測(cè)獲取水體泥沙質(zhì)量濃度,共獲取懸沙質(zhì)量濃度樣點(diǎn)59 個(gè);利用便攜地物光譜儀(ASD)采集水面光譜數(shù)據(jù),測(cè)量方法采用標(biāo)準(zhǔn)水上測(cè)量法(Mueller et al.,2003;唐軍武等,2004),共獲取8條光譜曲線。如圖1 所示,A1~A12 為原型觀測(cè)浮標(biāo)站點(diǎn),A13~A17為枯季船上采樣點(diǎn),B13~B17為洪季船上采樣點(diǎn),C1~C11 為枯季走航水面光譜測(cè)量點(diǎn)。其中定點(diǎn)觀測(cè)時(shí)間(表示北京時(shí)間,下同)為T(mén) 00:00—24:00,頻率為2 h/次,走航測(cè)量為T(mén) 09:00—16:00,頻率約為1 h/次。浮標(biāo)泥沙質(zhì)量濃度由浮標(biāo)標(biāo)定的濁度與泥沙質(zhì)量濃度的轉(zhuǎn)換模型計(jì)算得到,船上采樣通過(guò)實(shí)驗(yàn)室過(guò)濾重量法(American Public Health Association,1992)獲取。觀測(cè)實(shí)驗(yàn)所獲取的T 10:00—16:00 的泥沙質(zhì)量濃度及光譜數(shù)據(jù)見(jiàn)圖2所示。光譜符合內(nèi)陸或河口水體光譜曲線特征,泥沙質(zhì)量濃度梯度分級(jí)明顯。

    圖1 研究區(qū)地理位置及實(shí)測(cè)站點(diǎn)分布Fig.1 The study area and sampling sites'distribution

    圖2 泥沙質(zhì)量濃度野外實(shí)測(cè)數(shù)據(jù)集(a)與水面光譜測(cè)量數(shù)據(jù)(b)Fig.2 The field spectral and sediment concentration data

    1.3 遙感影像數(shù)據(jù)獲取及處理

    通過(guò)向中國(guó)資源衛(wèi)星應(yīng)用中心提交申請(qǐng),進(jìn)行衛(wèi)星編程拍攝,同步觀測(cè)期內(nèi)獲取了2019-12-28和2020-06-20各4景GF-4高頻次衛(wèi)星數(shù)據(jù),具體成像情況見(jiàn)表1。對(duì)獲取的遙感數(shù)據(jù)進(jìn)行幾何校正及大氣校正,使得所獲取的多時(shí)相遙感影像數(shù)據(jù)與野外觀測(cè)站點(diǎn)數(shù)據(jù)在空間位置及輻射水平上得以配準(zhǔn)。大氣校正采用基于影像自身統(tǒng)計(jì)特征的大氣校正法(Lee et al.,2007)。剔除走航實(shí)測(cè)異常光譜數(shù)據(jù)后,獲取2019-12-28 C1#T 10:00、C2#T 12:00、C6# T 14:00三個(gè)星-地同步觀測(cè)點(diǎn)位的實(shí)測(cè)光譜數(shù)據(jù),利用GF4號(hào)光譜響應(yīng)函數(shù)對(duì)實(shí)測(cè)光譜進(jìn)行積分,獲取band2~band5 四個(gè)波段的實(shí)測(cè)反射率值,并通過(guò)對(duì)比實(shí)測(cè)反射率值與GF4大氣校正反射值,評(píng)價(jià)大氣校正的準(zhǔn)確性。

    表1 GF-4高頻次衛(wèi)星數(shù)據(jù)成像情況Table 1 Imaging parameters of GF-4 high frequency satellite data

    1.4 反演模型構(gòu)建及驗(yàn)證

    根據(jù)遙感成像時(shí)間,分洪、枯季挑選對(duì)應(yīng)成像時(shí)刻的同步實(shí)測(cè)懸沙質(zhì)量濃度數(shù)據(jù),以及提取對(duì)應(yīng)實(shí)測(cè)點(diǎn)位空間位置GF4號(hào)遙感反射率數(shù)據(jù),剔除有云及云陰影污染、實(shí)測(cè)水表層懸沙質(zhì)量濃度缺失及數(shù)值異常的點(diǎn)位,獲取2019-12-28 共18 個(gè),2020-06-20共28個(gè)星-地同步觀測(cè)數(shù)據(jù)對(duì)。分別隨機(jī)選取15 和25 個(gè)形成建模數(shù)據(jù)集,用于模型的構(gòu)建,剩余為模型驗(yàn)證數(shù)據(jù)集,用于模型的驗(yàn)證。對(duì)于泥沙質(zhì)量濃度的遙感反演,已有研究構(gòu)建了眾多經(jīng)驗(yàn)?zāi)P?,其中指?shù)模型較為常見(jiàn),而對(duì)于因變量,常選取特征波段的比值,能夠?qū)崿F(xiàn)較好的擬合效果,減少大氣校正誤差的干擾(劉志國(guó)等,2006;Ye et al., 2014;潘洪洲 等,2019)。從2019-12-28 T 12:00成像的GF4數(shù)據(jù)不同渾濁水體反射率(相對(duì)清澈水體<渾濁水體<較為渾濁水體)可知(圖3),隨著泥沙質(zhì)量濃度的增加,紅及近紅外波段反射率抬升比較明顯,藍(lán)波段變化較小,與前人研究結(jié)論一致,可作為懸沙質(zhì)量濃度的特征波段(馬驊,2015;張俊儒,2010)。因此,分別以紅、近紅外與藍(lán)波段的比值作為自變量因子X(jué),構(gòu)建泥沙質(zhì)量濃度(SSC)指數(shù)反演模型。

    圖3 不同泥沙含量水體GF4不同波段反射率光譜差異Fig.3 Spectral differences of GF4 reflectance of water with different sediment content

    1.5 模型驗(yàn)證

    利用驗(yàn)證數(shù)據(jù)集開(kāi)展模型精度驗(yàn)證,誤差評(píng)價(jià)主要采用平均絕對(duì)百分比誤差及均方根誤差,計(jì)算方式為:

    式中:MAPE 為平均絕對(duì)百分比誤差;RMSE 均方根誤差;xir為第i變量;xim為變量的平均值;n為樣本個(gè)數(shù)。

    2 結(jié)果分析

    2.1 GF4影像大氣校正效果

    由圖4、表2可知,GF4影像大氣校正光譜與實(shí)測(cè)光譜較為吻合,較好地反映了水體的光譜特征。C1#、C6#點(diǎn)位藍(lán)、綠波段表現(xiàn)出一定程度的高估,C1#、C2#點(diǎn)位紅波段被略微低估。

    圖4 GF4光譜響應(yīng)函數(shù)及實(shí)測(cè)光譜與大氣校正結(jié)果對(duì)比Fig.4 Spectral response function of GF4 and comparison between measured spectrum and atmospheric correction results

    表2 各點(diǎn)位GF-4影像大氣校正各波段相對(duì)誤差Table 2 MAPE error of atmospheric correction results in each band %

    2.2 懸沙反演模型

    分 別 以GF4 影 像 反 射 率Rrs (band5)/Rrs(band2) [bandX表 示GF4 影 像 第X波 段,Rrs(bandX) 表示對(duì)應(yīng)波段反射率,下同]、Rrs(band4)/Rrs(band2)(紅藍(lán)比值)為自變量構(gòu)建河口區(qū)洪、枯季懸沙質(zhì)量濃度(SSC)反演模型(表3),枯季兩模型擬合系數(shù)接近0.6,差異不大,洪季紅藍(lán)比值模型表現(xiàn)相對(duì)較好,擬合系數(shù)為0.82。利用驗(yàn)證數(shù)據(jù)集對(duì)模型誤差進(jìn)行計(jì)算,枯季模型驗(yàn)證誤差整體小于洪季,最小平均絕對(duì)百分比誤差及均方根誤差分別為23.54%、0.007(見(jiàn)表3),綜合來(lái)看,紅藍(lán)比值模型表現(xiàn)相對(duì)較好,洪、枯季模型整體平均絕對(duì)百分比誤差為26.77%,因此,采用紅藍(lán)比值模型開(kāi)展河口區(qū)泥沙質(zhì)量濃度高頻次反演研究。

    表3 GF-4懸沙質(zhì)量濃度-指數(shù)反演模型Table 3 The exponential function inversion models of Suspended sediment using GF4 image data

    2.3 河口區(qū)潮周期泥沙輸移定量反演

    利用已構(gòu)建的泥沙反演模型,基于2019-12-28和2020-06-20多時(shí)相的GF4數(shù)據(jù),定量提取河口區(qū)潮周期泥沙質(zhì)量濃度(圖5-a)。2019-12-28河口區(qū)整體懸沙質(zhì)量濃度大部分介于0.01~0.056 kg/m3,少數(shù)區(qū)域懸沙質(zhì)量濃度高于0.056 kg/m3,高懸沙區(qū)域主要沿西四口門(mén)呈帶狀分布。外海懸沙質(zhì)量濃度基本都低于0.014 kg/m3。T 10:00—16:00 期間,河口東岸漲潮流較為顯著,外海低懸沙水體沿外伶仃島兩側(cè)上溯,伶仃洋整體懸沙質(zhì)量濃度降低。由圖5-b 可知,2020-06-20 河口整體懸沙質(zhì)量濃度大部分介于0.01~0.28 kg/m3,河 口西岸徑流輸沙較為顯著,各口門(mén)附近區(qū)域懸沙質(zhì)量濃度高于0.34 kg/m3,少部分區(qū)域質(zhì)量濃度高于0.50 kg/m3。而外海懸沙質(zhì)量濃度基本都低于0.02 kg/m3。隨著落潮,口門(mén)懸沙不斷向口外擴(kuò)散輸移,受虎門(mén)下泄潮流的頂托,洪奇門(mén)下泄泥沙擴(kuò)散區(qū)域受限,附近泥沙質(zhì)量濃度東西向呈現(xiàn)較為明顯的梯度變化。

    圖5 河口區(qū)潮周期泥沙質(zhì)量濃度定量提取結(jié)果(a~d.2019-12-28;e~h.2020-06-20)Fig.5 The high frequency remote sensing inversion results of sediment concentration in tidal cycle in estuarine area(a~d.Dec.28th,2019;e~h.June 20th,2020)

    3 討論

    隨著漲、落潮不同的潮汐過(guò)程,河口區(qū)泥沙質(zhì)量濃度不斷發(fā)生變化,特別是伶仃洋區(qū)域,泥沙質(zhì)量濃度表現(xiàn)出明顯的潮汐過(guò)程變化。為了分析該變化,選取伶仃洋水域,提取洪、枯季多時(shí)相泥沙質(zhì)量濃度特征等值線。如圖6-a 所示,從主要潮位站潮位數(shù)據(jù)來(lái)看,2019-12-28主要為枯季漲潮,徑流輸沙較少,外海漲潮水體的沖淡為影響伶仃洋泥沙質(zhì)量濃度空間分布的主要因素,因此,枯季選取低質(zhì)量濃度0.02 kg/m3等值線進(jìn)行分析;如圖6-b 所示,而2020-06-20主要潮汐過(guò)程為洪季落潮,徑流泥沙擴(kuò)散為影響河口區(qū)泥沙質(zhì)量濃度空間分布的主要因素,因此,洪季選取接近徑流輸沙質(zhì)量濃度的0.10 kg/m3等值線進(jìn)行泥沙擴(kuò)散分析。由圖7可知,受漲潮流影響,伶仃洋整體泥沙質(zhì)量濃度降低,其泥沙空間分布主要受沿河口東岸及香港西側(cè)兩股漲潮流影響,其中沿河口東岸漲潮流率先上漲,受沖淡影響,東岸沿岸區(qū)域泥沙質(zhì)量濃度從東南向西北逐漸降低;香港西側(cè)的漲潮于內(nèi)伶仃島兩側(cè)分流形成兩股東北向、西北向漲潮流,流經(jīng)區(qū)域泥沙質(zhì)量濃度沖淡作用顯著,西北向分流直接流入圖7-c 中區(qū)域B,該區(qū)域受人為活動(dòng)影響,海底地形下切顯著(Li et al.,2010;Wu et al.,2016),潮汐動(dòng)力加強(qiáng),河床泥沙源減少,泥沙再懸浮至表層的量大大減少,在漲潮流的影響下形成大片低懸沙質(zhì)量濃度區(qū)域(van Maren et al., 2016;詹偉康 等,2019)。T14:00—16:00內(nèi)伶仃島、澳門(mén)、三灶等河口下游潮位站進(jìn)漲潮階段,上游橫門(mén)、舢板洲延遲落潮(圖7-c),交椅灣西南側(cè)區(qū)域沖淡區(qū)域仍有上移,區(qū)域A附近水域泥沙質(zhì)量濃度仍存在被沖淡現(xiàn)象,區(qū)域B 低質(zhì)量濃度水體沿內(nèi)伶仃島西側(cè)呈西南向開(kāi)始下泄。

    圖6 GF4成像時(shí)刻各潮位站潮位過(guò)程Fig.6 The tidal process of tide stations at GF4 imaging time

    圖7 2019-12-28(a、b、c)和2020-06-20(d、e、f)泥沙質(zhì)量濃度分布、特征等值線及前一時(shí)刻特征等值線Fig.7 The Sediment concentration distribution,isopleth of the current and previous moment at Dec.28th,2019(a,b,c)and June 20th,2020(d,e,f)

    洪季落潮階段,隨著落潮,伶仃洋西岸高懸沙水體下泄。橫門(mén)漲平至初落階段,圖7-d 中區(qū)域A泥沙向東南及偏南向輸移,隨著河口東岸落潮流的增強(qiáng),泥沙輸移受虎門(mén)落潮流頂托,東南向輸移受限,輸移方向向南向偏轉(zhuǎn),圖7-f中區(qū)域B,落潮流與內(nèi)伶仃兩側(cè)分流,沿內(nèi)伶仃兩側(cè)成半環(huán)形下泄,區(qū)域B泥沙輸移受東岸落潮流頂托顯著,隨著泥沙南向下泄推移、東岸落潮流東南向偏移,圖7-f 中區(qū)域C,泥沙向東南有所擴(kuò)散,隨著西岸落潮流不斷向南擴(kuò)散,至T 16:00時(shí),區(qū)域D,泥沙擴(kuò)散至澳門(mén)附近水域。

    根據(jù)影像多時(shí)相反演有效性及連續(xù)性情況,挑選無(wú)云、陰影等影響,時(shí)序較為連續(xù)的6 個(gè)點(diǎn)位,對(duì)比實(shí)測(cè)及遙感多時(shí)相泥沙時(shí)序變化(圖8)。整體上,兩者表現(xiàn)出較為一致的特征,枯季河口上游泥沙質(zhì)量濃度時(shí)序變化基本呈現(xiàn)遞減趨勢(shì),河口下游受漲落潮交替影響,部分點(diǎn)位呈現(xiàn)先減后增的趨勢(shì),如枯季A7 點(diǎn)位;洪季泥沙質(zhì)量濃度最大值出現(xiàn)在落急時(shí)刻(T 14:00)前后,泥沙質(zhì)量濃度變化大致呈先大增后小減的趨勢(shì)。由表4可知,A7點(diǎn)位T 10:00—12:00為漲潮階段,T 12:00—14:00為漲落潮交替階段,T 14:00—16:00轉(zhuǎn)為落潮。A7點(diǎn)位靠近沿岸,上游高質(zhì)量濃度泥沙水體沿澳門(mén)機(jī)場(chǎng)沿岸下泄,漲潮轉(zhuǎn)落潮階段,表層水體易受上游下泄泥沙影響,因此,A7點(diǎn)位泥沙質(zhì)量濃度時(shí)序變化呈現(xiàn)先減后增的趨勢(shì)(圖9)。由圖8可知,部分點(diǎn)位遙感反演值與站點(diǎn)實(shí)測(cè)值之間差異較大,如枯季A1、A2點(diǎn)。從枯季遙感反演及實(shí)測(cè)點(diǎn)位泥沙質(zhì)量濃度整體變化看,枯季潮流攜帶泥沙較少,水體泥沙質(zhì)量濃度低。如圖6、8及表4所示,A1點(diǎn)T 10:00—14:00為漲潮階段,T 14:00—16:00為漲落潮交替階段,隨著漲潮流的上漲,受潮流沖淡作用影響,A1點(diǎn)位泥沙質(zhì)量濃度呈現(xiàn)逐漸遞減趨勢(shì),遙感反演值時(shí)序變化特征與之相符,但與實(shí)測(cè)值時(shí)序趨勢(shì)不一致,主要有以下原因:從空間采樣的差異性看,遙感反演采用的是50 m分辨率影像的面上采樣,而實(shí)測(cè)站點(diǎn)為點(diǎn)狀采樣,采樣尺度存在差異;且浮標(biāo)站點(diǎn)濁度與懸沙質(zhì)量濃度的轉(zhuǎn)換也存在一定的誤差;另外,GF4號(hào)遙感影像大氣校正誤差及反演模型的不穩(wěn)定性,都可能成為實(shí)測(cè)值與反演值不一致性的影響因素。由于缺少測(cè)量時(shí)刻其他描述材料,如水樣采集時(shí)刻船只擾動(dòng)影響等情況,不再對(duì)A2 測(cè)點(diǎn)的差異性作進(jìn)一步分析。

    圖8 部分實(shí)測(cè)站點(diǎn)及對(duì)應(yīng)遙感反演泥沙時(shí)序變化Fig.8 The comparison of sediment concentration daily changes between remote sensing inversion and measured datas

    圖9 2019-12-28 A7實(shí)測(cè)站點(diǎn)T 14:00泥沙分布疊加T 14:00 及T 16:00 0.02 kg/m3等值線Fig.9 The Sediment distribution at 14:00 and 0.02 kg/m3isopleth at T 14:00/16:00 at A7 site on Dec.28th,2019

    表4 A1及A7點(diǎn)位2019-12-28實(shí)測(cè)潮流流向Table 4 The measured tidal current direction at A1 and A7 site on Dec.28th,2019

    4 結(jié)論

    本文利用洪、枯季兩批多時(shí)相GF4 影像數(shù)據(jù),結(jié)合地面同步實(shí)測(cè)站點(diǎn)數(shù)據(jù),通過(guò)構(gòu)建反演模型實(shí)現(xiàn)珠江河口潮周期內(nèi)泥沙質(zhì)量濃度的高頻次反演,并以伶仃洋為例分析成像期間的泥沙動(dòng)態(tài)輸移特征,結(jié)果表明:1)基于GF4 影像紅、藍(lán)波段比值及近紅外、藍(lán)波段比值能很好地指示河口泥沙含量,基于兩類(lèi)比值構(gòu)建的模型表現(xiàn)較優(yōu),其中紅、藍(lán)比值模型相對(duì)較好,整體平均絕對(duì)百分比誤差為26.77%,多時(shí)相遙感反演結(jié)果與實(shí)測(cè)站點(diǎn)的時(shí)序變化特征較為一致;2)枯季漲潮水體對(duì)泥沙分布起明顯的沖淡作用,枯季漲潮階段伶仃洋整體泥沙質(zhì)量濃度有所減小,河口東岸漲潮流較為強(qiáng)勁,由于潮流的沖淡作用,沿伶仃洋東岸及中灘大面積區(qū)域形成低泥沙質(zhì)量濃度區(qū),高泥沙質(zhì)量濃度分布區(qū)域受漲潮擠壓,沿伶仃洋西岸呈帶狀分布;3)洪季河口高泥沙質(zhì)量濃度水體沿伶仃洋西岸呈帶狀分布,落潮階段泥沙最高質(zhì)量濃度出現(xiàn)在落急時(shí)刻左右,隨著落潮向西南或偏南向擴(kuò)散輸移,虎門(mén)下泄徑潮流泥沙質(zhì)量濃度相對(duì)較低,受該潮流下泄影響,泥沙擴(kuò)散受到擠壓,西岸東南向泥沙擴(kuò)散路徑向南推移,南向泥沙輸移擴(kuò)散至澳門(mén)水域附近。

    本研究為復(fù)雜河口泥沙輸移大范圍業(yè)務(wù)化監(jiān)測(cè)提供了新的思路,有助于提升對(duì)珠江河口潮周期內(nèi)泥沙動(dòng)態(tài)輸移變化的認(rèn)識(shí),同時(shí)為珠江河口泥沙輸移數(shù)值模擬與遙感監(jiān)測(cè)提供了交叉驗(yàn)證的新思路,為河口沉積物輸運(yùn)、泥沙來(lái)源及泥沙數(shù)值模擬研究等提供數(shù)據(jù)支撐。隨著遙感技術(shù)的發(fā)展,遙感光譜分辨率及時(shí)間分辨率不斷提高,泥沙潮遙感反演研究由多頻次向高頻次擴(kuò)展,監(jiān)測(cè)維度也將從表層向垂直層面擴(kuò)展,實(shí)分析泥沙質(zhì)量濃度垂直層面分布及準(zhǔn)實(shí)時(shí)輸移,進(jìn)而可以估算水體某過(guò)流斷面的泥沙量,甚至可以估算泥沙通量,可為泥沙動(dòng)態(tài)輸移特征的深入研究提供支撐。

    猜你喜歡
    懸沙河口區(qū)河口
    漢江河口區(qū)四明灘段四明上護(hù)岸坍岸修復(fù)技術(shù)
    近岸懸沙垂線分布多元線性回歸分析
    臺(tái)風(fēng)對(duì)長(zhǎng)江口表層懸沙濃度的影響
    “傳承—信仰—未來(lái)”教育主題的設(shè)計(jì)者
    ——記山東省東營(yíng)市河口區(qū)河安小學(xué)校長(zhǎng)劉波
    他們?yōu)槭裁催x擇河口
    河口,我們的家
    特殊的河口水
    河口
    閩江下游及河口區(qū)表層沉積物多種微量元素特征
    東山灣波浪對(duì)懸沙濃度場(chǎng)影響的數(shù)值模擬研究
    99久久综合免费| 一级黄片播放器| 亚洲av电影在线观看一区二区三区| 国产成人freesex在线| 国产精品一区二区三区四区免费观看| 欧美日韩一区二区视频在线观看视频在线| 日韩免费高清中文字幕av| 久久97久久精品| 亚洲欧洲国产日韩| 在线观看人妻少妇| 少妇裸体淫交视频免费看高清| 亚洲精品日韩av片在线观看| 日韩一本色道免费dvd| 五月开心婷婷网| 亚洲av日韩在线播放| 国产黄片美女视频| 99热这里只有精品一区| a级毛色黄片| 少妇人妻 视频| 国产精品国产三级国产av玫瑰| 亚洲欧洲国产日韩| 国产成人精品婷婷| 久久久午夜欧美精品| 日韩成人伦理影院| 99久久人妻综合| 国产精品国产三级国产专区5o| 天堂俺去俺来也www色官网| 成人二区视频| 99热这里只有是精品50| 亚洲第一av免费看| 午夜影院在线不卡| 国产亚洲精品久久久com| 久久青草综合色| 久久影院123| 91aial.com中文字幕在线观看| 久久女婷五月综合色啪小说| 国产精品久久久久久久电影| 全区人妻精品视频| 人人妻人人看人人澡| 乱码一卡2卡4卡精品| 久久久久久人妻| 99热这里只有精品一区| 插阴视频在线观看视频| 亚洲精品乱久久久久久| 伊人亚洲综合成人网| 九草在线视频观看| 亚洲av二区三区四区| 国产成人精品福利久久| 亚洲精品成人av观看孕妇| 我要看黄色一级片免费的| 女人精品久久久久毛片| 国产熟女欧美一区二区| 欧美三级亚洲精品| 欧美bdsm另类| 国产精品人妻久久久影院| 日韩三级伦理在线观看| 国产精品欧美亚洲77777| 欧美日韩视频高清一区二区三区二| 日本爱情动作片www.在线观看| 国产免费一区二区三区四区乱码| 男人和女人高潮做爰伦理| 国产日韩欧美视频二区| 少妇精品久久久久久久| 成人毛片60女人毛片免费| 精品一区二区免费观看| 亚洲av免费高清在线观看| 91成人精品电影| av天堂中文字幕网| 亚洲精品乱久久久久久| 97精品久久久久久久久久精品| 国产淫语在线视频| 亚洲激情五月婷婷啪啪| 青青草视频在线视频观看| 三级国产精品片| 在线播放无遮挡| 国产极品天堂在线| 欧美日韩视频高清一区二区三区二| 午夜激情福利司机影院| kizo精华| 永久网站在线| 又粗又硬又长又爽又黄的视频| 男人和女人高潮做爰伦理| 黄色欧美视频在线观看| 蜜臀久久99精品久久宅男| av天堂久久9| 秋霞在线观看毛片| 久久久久国产网址| av免费在线看不卡| 熟女人妻精品中文字幕| 麻豆成人av视频| av播播在线观看一区| 99九九线精品视频在线观看视频| 久久久久久久国产电影| 中文字幕人妻丝袜制服| 成人影院久久| 亚洲av不卡在线观看| 人妻制服诱惑在线中文字幕| 桃花免费在线播放| 水蜜桃什么品种好| 国产伦理片在线播放av一区| 又大又黄又爽视频免费| 日本欧美国产在线视频| 国产精品不卡视频一区二区| 亚洲第一区二区三区不卡| av又黄又爽大尺度在线免费看| 欧美 亚洲 国产 日韩一| 一个人免费看片子| 亚洲精品日韩在线中文字幕| 国产av国产精品国产| 亚洲av不卡在线观看| 黄色日韩在线| 日本猛色少妇xxxxx猛交久久| 欧美精品亚洲一区二区| 日本黄色日本黄色录像| 亚洲一级一片aⅴ在线观看| 女性被躁到高潮视频| a级一级毛片免费在线观看| a级毛片在线看网站| 国产视频首页在线观看| 特大巨黑吊av在线直播| 精品熟女少妇av免费看| 国产亚洲最大av| 麻豆精品久久久久久蜜桃| 久久久久久久久久成人| 日日撸夜夜添| 99热6这里只有精品| 色5月婷婷丁香| av免费观看日本| 美女内射精品一级片tv| 99久久精品一区二区三区| 亚洲一区二区三区欧美精品| 国产精品久久久久久av不卡| 中文字幕久久专区| 如日韩欧美国产精品一区二区三区 | 日本黄色日本黄色录像| 97超视频在线观看视频| 国产中年淑女户外野战色| 亚洲av在线观看美女高潮| 久久这里有精品视频免费| 国产精品.久久久| 天天操日日干夜夜撸| 日本午夜av视频| 久久久国产一区二区| 国产高清国产精品国产三级| 男女国产视频网站| 免费不卡的大黄色大毛片视频在线观看| 视频区图区小说| 91久久精品电影网| 黄色日韩在线| 欧美日本中文国产一区发布| 国产在线男女| 日本欧美视频一区| av黄色大香蕉| 亚洲精品国产av蜜桃| 热re99久久精品国产66热6| 黄色视频在线播放观看不卡| 亚洲四区av| 日韩欧美精品免费久久| 国产成人a∨麻豆精品| 中文字幕av电影在线播放| 观看美女的网站| 国产精品伦人一区二区| 99视频精品全部免费 在线| 亚洲国产欧美在线一区| 亚洲人与动物交配视频| 久久久久久久久久久免费av| 色哟哟·www| 免费看av在线观看网站| 亚洲av日韩在线播放| 国产精品一区二区在线不卡| 亚洲av在线观看美女高潮| 国产精品一区二区性色av| 中国国产av一级| 女的被弄到高潮叫床怎么办| 成人毛片60女人毛片免费| 成人18禁高潮啪啪吃奶动态图 | 亚洲精品乱码久久久v下载方式| 中文字幕人妻丝袜制服| 久久青草综合色| 中文在线观看免费www的网站| 久久久久久久久久久久大奶| 国产免费视频播放在线视频| 精品一区二区三卡| 日本午夜av视频| 午夜福利影视在线免费观看| 亚洲欧美成人精品一区二区| 丰满人妻一区二区三区视频av| 久久精品夜色国产| 九九在线视频观看精品| 欧美另类一区| 丁香六月天网| 街头女战士在线观看网站| 欧美日韩国产mv在线观看视频| 日韩电影二区| 免费播放大片免费观看视频在线观看| av女优亚洲男人天堂| 人妻一区二区av| 亚洲怡红院男人天堂| 精品人妻熟女av久视频| 成人亚洲精品一区在线观看| 大片电影免费在线观看免费| 欧美成人午夜免费资源| av福利片在线观看| 亚洲欧洲国产日韩| 亚洲婷婷狠狠爱综合网| 国产欧美日韩一区二区三区在线 | 在线播放无遮挡| 22中文网久久字幕| 97在线视频观看| 乱系列少妇在线播放| 日本欧美视频一区| 纵有疾风起免费观看全集完整版| 男女边吃奶边做爰视频| 三级国产精品片| 精品人妻熟女毛片av久久网站| 亚洲熟女精品中文字幕| 久久国产亚洲av麻豆专区| 亚洲精品国产av蜜桃| 黄色一级大片看看| av免费在线看不卡| 国产精品久久久久久精品电影小说| 少妇被粗大的猛进出69影院 | 亚洲va在线va天堂va国产| 亚洲精品乱码久久久v下载方式| 女人久久www免费人成看片| 日韩中字成人| 亚洲精品久久午夜乱码| 久久精品国产亚洲av天美| 亚洲中文av在线| 日日啪夜夜撸| 中文字幕免费在线视频6| 肉色欧美久久久久久久蜜桃| 少妇猛男粗大的猛烈进出视频| 在线观看一区二区三区激情| 最新中文字幕久久久久| a级毛色黄片| 成人亚洲欧美一区二区av| 噜噜噜噜噜久久久久久91| 精品亚洲乱码少妇综合久久| 国产熟女午夜一区二区三区 | 观看免费一级毛片| 26uuu在线亚洲综合色| 国产精品一区www在线观看| 日日摸夜夜添夜夜添av毛片| 丰满人妻一区二区三区视频av| 欧美3d第一页| 午夜福利在线观看免费完整高清在| 久久人人爽av亚洲精品天堂| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 亚洲精品中文字幕在线视频 | av卡一久久| 国产伦精品一区二区三区四那| 国产一区二区三区av在线| 热99国产精品久久久久久7| 日本黄色日本黄色录像| 丰满乱子伦码专区| 欧美精品高潮呻吟av久久| 成人黄色视频免费在线看| 国产片特级美女逼逼视频| 国产精品.久久久| 国产视频首页在线观看| 色视频www国产| 久久久久久久亚洲中文字幕| 国产在线免费精品| 亚洲国产av新网站| 99九九在线精品视频 | 日韩伦理黄色片| 99久久综合免费| 又黄又爽又刺激的免费视频.| 一级,二级,三级黄色视频| 综合色丁香网| 精品人妻一区二区三区麻豆| 国产av码专区亚洲av| 中文精品一卡2卡3卡4更新| 久久人人爽人人片av| a级片在线免费高清观看视频| 自拍偷自拍亚洲精品老妇| 午夜福利影视在线免费观看| 日韩熟女老妇一区二区性免费视频| 成人影院久久| 秋霞伦理黄片| av又黄又爽大尺度在线免费看| 亚洲精品国产av蜜桃| 日韩一本色道免费dvd| 亚洲va在线va天堂va国产| 亚洲国产日韩一区二区| 亚洲欧洲精品一区二区精品久久久 | 最近2019中文字幕mv第一页| 制服丝袜香蕉在线| 男男h啪啪无遮挡| 欧美人与善性xxx| 欧美另类一区| 两个人免费观看高清视频 | 国产av一区二区精品久久| 精品一区二区免费观看| 黑人巨大精品欧美一区二区蜜桃 | 国产女主播在线喷水免费视频网站| 精品卡一卡二卡四卡免费| 国产伦在线观看视频一区| 少妇的逼好多水| 丝袜喷水一区| 桃花免费在线播放| av线在线观看网站| 又大又黄又爽视频免费| 午夜精品国产一区二区电影| 日本午夜av视频| 人妻系列 视频| 国产精品一区二区性色av| 性色avwww在线观看| 精品人妻熟女毛片av久久网站| 亚洲欧洲精品一区二区精品久久久 | 久久久国产欧美日韩av| 精品国产一区二区久久| 中文字幕制服av| av播播在线观看一区| 欧美xxxx性猛交bbbb| 夜夜看夜夜爽夜夜摸| 卡戴珊不雅视频在线播放| 尾随美女入室| 91久久精品国产一区二区三区| 精品久久久噜噜| 国产精品秋霞免费鲁丝片| 最近的中文字幕免费完整| 国产精品一区二区性色av| 亚洲va在线va天堂va国产| 色94色欧美一区二区| 亚洲色图综合在线观看| 久久久久久久久久成人| 午夜免费鲁丝| 国产探花极品一区二区| 91在线精品国自产拍蜜月| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 男女无遮挡免费网站观看| 色吧在线观看| 国产亚洲av片在线观看秒播厂| 亚洲国产av新网站| 十八禁高潮呻吟视频 | 中文在线观看免费www的网站| 美女主播在线视频| 中文字幕制服av| 在线观看免费日韩欧美大片 | 久久国产精品男人的天堂亚洲 | 成人毛片a级毛片在线播放| 啦啦啦视频在线资源免费观看| 亚洲精品乱码久久久v下载方式| 精品一区二区三区视频在线| 另类亚洲欧美激情| 欧美日韩在线观看h| 国产av国产精品国产| 一个人免费看片子| 91久久精品电影网| 青青草视频在线视频观看| 黑人巨大精品欧美一区二区蜜桃 | 国产av国产精品国产| 三上悠亚av全集在线观看 | 国产精品久久久久成人av| 国产在线男女| 视频区图区小说| 水蜜桃什么品种好| 国产在线免费精品| 久久人人爽人人爽人人片va| 国产精品伦人一区二区| 熟女电影av网| 免费大片18禁| 国产高清三级在线| 亚洲欧美精品专区久久| 最近手机中文字幕大全| 久久久欧美国产精品| 日日撸夜夜添| 超碰97精品在线观看| 你懂的网址亚洲精品在线观看| 一区二区三区乱码不卡18| 久久女婷五月综合色啪小说| 午夜老司机福利剧场| 九九爱精品视频在线观看| 最新中文字幕久久久久| 高清在线视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 中文在线观看免费www的网站| 桃花免费在线播放| 欧美97在线视频| 大码成人一级视频| 男女无遮挡免费网站观看| 免费久久久久久久精品成人欧美视频 | 亚洲欧美成人精品一区二区| 日韩精品有码人妻一区| 老司机影院成人| 精品国产露脸久久av麻豆| 99九九在线精品视频 | 亚洲欧美日韩另类电影网站| 亚洲精品,欧美精品| 日韩一区二区三区影片| 一区在线观看完整版| 三级国产精品片| 一级毛片黄色毛片免费观看视频| 一区二区三区免费毛片| 天堂中文最新版在线下载| 啦啦啦中文免费视频观看日本| 男女无遮挡免费网站观看| 久久人人爽人人爽人人片va| 国产成人午夜福利电影在线观看| 色视频在线一区二区三区| 日韩,欧美,国产一区二区三区| 日本爱情动作片www.在线观看| 欧美激情极品国产一区二区三区 | 青青草视频在线视频观看| 国产探花极品一区二区| 王馨瑶露胸无遮挡在线观看| 午夜福利在线观看免费完整高清在| 女性被躁到高潮视频| 男的添女的下面高潮视频| 亚洲人与动物交配视频| 黄片无遮挡物在线观看| 嘟嘟电影网在线观看| 婷婷色综合大香蕉| 伊人久久国产一区二区| 在线观看三级黄色| 两个人免费观看高清视频 | av天堂中文字幕网| 男人和女人高潮做爰伦理| 久久热精品热| 婷婷色av中文字幕| 日本wwww免费看| 久久久国产欧美日韩av| 男女边吃奶边做爰视频| 欧美激情国产日韩精品一区| 啦啦啦视频在线资源免费观看| 国产男女超爽视频在线观看| 在线 av 中文字幕| 蜜臀久久99精品久久宅男| 免费看日本二区| 又大又黄又爽视频免费| 丁香六月天网| 一级毛片电影观看| 国产成人freesex在线| 国产精品久久久久久av不卡| 午夜福利在线观看免费完整高清在| 久久久久精品性色| 国产精品福利在线免费观看| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄| 免费黄频网站在线观看国产| √禁漫天堂资源中文www| 亚洲国产精品999| 我要看黄色一级片免费的| 亚洲成人av在线免费| 一本一本综合久久| 久久99精品国语久久久| 欧美成人午夜免费资源| 啦啦啦在线观看免费高清www| 日韩欧美一区视频在线观看 | 日韩熟女老妇一区二区性免费视频| 一级毛片我不卡| 午夜福利网站1000一区二区三区| 麻豆乱淫一区二区| 国内揄拍国产精品人妻在线| 久久国产精品男人的天堂亚洲 | 又粗又硬又长又爽又黄的视频| 人人妻人人看人人澡| 亚洲精品久久午夜乱码| 国产有黄有色有爽视频| 久久久久国产精品人妻一区二区| av福利片在线| 成人美女网站在线观看视频| 中国国产av一级| 少妇高潮的动态图| 男人和女人高潮做爰伦理| 婷婷色麻豆天堂久久| 亚洲国产成人一精品久久久| 亚洲精品色激情综合| 国产午夜精品一二区理论片| 久久 成人 亚洲| 熟女电影av网| 狂野欧美白嫩少妇大欣赏| 91久久精品电影网| 色婷婷久久久亚洲欧美| 美女脱内裤让男人舔精品视频| 色视频在线一区二区三区| 久久毛片免费看一区二区三区| 亚洲内射少妇av| 亚洲性久久影院| 午夜福利视频精品| 亚洲精品aⅴ在线观看| 国产69精品久久久久777片| 99久久精品热视频| h视频一区二区三区| av一本久久久久| 黄色一级大片看看| 热99国产精品久久久久久7| 国产极品天堂在线| 高清毛片免费看| 中文精品一卡2卡3卡4更新| 97超碰精品成人国产| 一二三四中文在线观看免费高清| 免费观看av网站的网址| 麻豆乱淫一区二区| 国产伦精品一区二区三区视频9| 国产色婷婷99| 青青草视频在线视频观看| 久久国产精品大桥未久av | 如日韩欧美国产精品一区二区三区 | 日韩一区二区三区影片| 视频区图区小说| 免费人成在线观看视频色| 美女xxoo啪啪120秒动态图| 久久久国产精品麻豆| 在线 av 中文字幕| 深夜a级毛片| tube8黄色片| 久久午夜综合久久蜜桃| 夫妻性生交免费视频一级片| 国产美女午夜福利| 欧美bdsm另类| 国产淫语在线视频| 亚洲精品久久午夜乱码| 成年美女黄网站色视频大全免费 | 色5月婷婷丁香| av国产久精品久网站免费入址| 男的添女的下面高潮视频| 美女视频免费永久观看网站| 一级毛片 在线播放| 少妇的逼水好多| 国产精品久久久久久av不卡| 国内少妇人妻偷人精品xxx网站| 欧美 日韩 精品 国产| 国产精品人妻久久久影院| 色5月婷婷丁香| 午夜老司机福利剧场| 亚洲一区二区三区欧美精品| 久久精品国产亚洲av天美| 国产精品欧美亚洲77777| 女人久久www免费人成看片| 国产老妇伦熟女老妇高清| 一区在线观看完整版| 免费少妇av软件| 天美传媒精品一区二区| 三级国产精品欧美在线观看| 欧美成人精品欧美一级黄| 午夜激情久久久久久久| 人妻制服诱惑在线中文字幕| 亚洲欧美中文字幕日韩二区| 最黄视频免费看| 成人综合一区亚洲| 日日爽夜夜爽网站| 寂寞人妻少妇视频99o| 亚洲精品中文字幕在线视频 | 亚洲图色成人| 免费久久久久久久精品成人欧美视频 | 亚洲国产最新在线播放| 久久久亚洲精品成人影院| 自线自在国产av| 欧美另类一区| 免费av不卡在线播放| 熟女人妻精品中文字幕| 少妇人妻 视频| 国产av一区二区精品久久| 91精品国产国语对白视频| 亚洲精品aⅴ在线观看| 香蕉精品网在线| 一本—道久久a久久精品蜜桃钙片| 尾随美女入室| 我要看黄色一级片免费的| 日本av免费视频播放| 又黄又爽又刺激的免费视频.| 欧美性感艳星| 你懂的网址亚洲精品在线观看| 热re99久久精品国产66热6| 亚洲精品色激情综合| 精品人妻熟女毛片av久久网站| 老司机影院成人| 伦精品一区二区三区| 人妻人人澡人人爽人人| 亚洲美女黄色视频免费看| 国产一区有黄有色的免费视频| 亚洲情色 制服丝袜| 熟女电影av网| 国产精品伦人一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 99久久精品国产国产毛片| 久久国产乱子免费精品| 久久精品熟女亚洲av麻豆精品| 免费观看无遮挡的男女| 大话2 男鬼变身卡| 精品99又大又爽又粗少妇毛片| 国产精品三级大全| 黑丝袜美女国产一区| 热re99久久精品国产66热6| 中文精品一卡2卡3卡4更新| 女性生殖器流出的白浆| 免费黄色在线免费观看| 亚洲av男天堂| 麻豆成人av视频| 精品一区二区三区视频在线| 日韩不卡一区二区三区视频在线| 国产精品人妻久久久久久| 春色校园在线视频观看| 日本午夜av视频| a级一级毛片免费在线观看| 丰满少妇做爰视频| 午夜av观看不卡| 欧美97在线视频| 在线观看免费高清a一片| 汤姆久久久久久久影院中文字幕| 国产高清三级在线| 免费观看无遮挡的男女| 久久久午夜欧美精品| 国产一区有黄有色的免费视频| 日日爽夜夜爽网站| 亚洲av中文av极速乱| 国产精品麻豆人妻色哟哟久久| 日本爱情动作片www.在线观看| 在线观看av片永久免费下载| 中文在线观看免费www的网站|