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

    油藻混合水體遙感反射比光譜分析

    2021-11-17 09:13:22駱蔚健黃妙芬楊俊杰
    海洋科學(xué) 2021年10期

    駱蔚健, 黃妙芬, 楊俊杰, 王 欣

    油藻混合水體遙感反射比光譜分析

    駱蔚健1, 黃妙芬2, 楊俊杰3, 王 欣2

    (1.廣東海洋大學(xué)海洋與氣象學(xué)院, 廣東 湛江 524088; 2.廣東海洋大學(xué)數(shù)學(xué)與計(jì)算機(jī)學(xué)院, 廣東 湛江 524088; 3.廣東海洋大學(xué)電子與信息工程學(xué)院, 廣東 湛江 524088)

    輻射傳輸模式HydroLight是研究水體輻射傳輸特性的有效工具, 同時(shí)也是進(jìn)行石油類水體輻射傳輸特性的有效模型。本文基于2018年8月在遼寧大連港海域?qū)崪y的石油類污染水體的表觀及固有光學(xué)量數(shù)據(jù), 通過設(shè)置不同濃度的配比模擬試驗(yàn), 利用HydroLight對油類物質(zhì)和藻類物質(zhì)的單一組分和兩者混合水體的遙感反射比(remote sensing reflectance,rs)光譜進(jìn)行了模擬。模擬結(jié)果表明: (1)在僅含油類物質(zhì)單一組分的水體中,rs隨著油濃度oil的變化分為兩個(gè)特征波段: 400~480 nm和480~700 nm。在400~480 nm波段范圍內(nèi)rs隨oil的增加而減小, 在480~700 nm隨oil的增大而增大; (2)在僅含藻類物質(zhì)的單一成分水體中, 葉綠素的光譜特性需要其濃度chl達(dá)到一定值才能表現(xiàn)出來, 在低chl時(shí)的rs光譜特性反映為純水的光譜特性; (3)在油藻混合水體中, 隨oil的增加會(huì)增大葉綠素的rs, 但不會(huì)明顯改變?nèi)~綠素的rs隨波長的變化趨勢, 這說明油藻混合水體的光譜形狀主要受葉綠素的影響, 油類物質(zhì)的存在只改變光譜的量值。利用HydroLight對石油類污染水體的rs展開研究, 有助于豐富水色遙感基礎(chǔ)研究, 對完善各類水體生物光學(xué)模型研究具有重要的理論意義。

    遙感反射比; HydroLight; 油類物質(zhì); 葉綠素

    由于水體石油類污染問題會(huì)給近海生態(tài)環(huán)境和漁業(yè)生產(chǎn)帶來巨大災(zāi)難性的后果, 因而對水體(河流、湖泊、海洋)油污染檢測是目前水環(huán)境監(jiān)測的熱點(diǎn)問題。遙感監(jiān)測技術(shù)因其大面積、快速、動(dòng)態(tài)和低成本的優(yōu)勢, 成為檢測石油類污染的重要手段之一。在水色遙感領(lǐng)域, 遙感反射比(remote sensing reflectance,rs)屬于表觀光學(xué)量, 它包含一定深度內(nèi)水層的光學(xué)信息, 因而成為反演水體組分濃度的基本物理量。目前以機(jī)理分析為基礎(chǔ), 反演Ⅱ類水體各水色因子的生物-光學(xué)遙感模型也較為成熟, 該模型的核心思想是基于rs與水體組分之間的吸收特性和散射特性的相互關(guān)系建立的。因而要提高相應(yīng)模型的反演精度, 需確定各組分對rs的貢獻(xiàn)。

    在石油類污染水體中, 利用遙感技術(shù)來提取水體石油類污染信息的相關(guān)研究已經(jīng)取得了長足的進(jìn)展, 主要集中在水體石油類物質(zhì)固有光學(xué)特性(inhe-rent optical properties, IOPs)、表觀光學(xué)特征(apparent optical properties, AOPs)和熒光特性等方面[1-5], 并分別建立了基于這些特性反演水體石油類含量的遙感反演模型[4, 6]。由于懸浮物質(zhì)會(huì)吸附油類物質(zhì), 黃妙芬等[7]對油和石英砂混合水體各組分對后向散射系數(shù)的貢獻(xiàn)分離算法展開了研究, 提出了分離算法, 為油類物質(zhì)和顆粒物后向散射系數(shù)分離提供技術(shù)支持。掌握油類物質(zhì)在含油水體中對rs的影響, 有助于提高遙感反演油類物質(zhì)濃度模型的精度, 對于構(gòu)建精度較高的反演含油水體組分的生物光學(xué)模型和推進(jìn)含油水體的相關(guān)研究具有重要作用。如今關(guān)于水色三要素(葉綠素、懸浮泥沙和黃色物質(zhì))的研究均取得了豐碩的成果, 對于葉綠素光譜, GITELSON[8]的研究表明含藻類水體rs最顯著的光譜特征是在680~715 nm出現(xiàn)反射峰, 其位置從680 nm開始, 在700 nm附近峰值達(dá)到最高, 峰值高度與葉綠素濃度密切相關(guān)。但是, 目前關(guān)于石油類物質(zhì)與水體組分混合后對rs的影響和貢獻(xiàn)的研究還鮮見報(bào)道。

    輻射傳輸模式HydroLight是國際公認(rèn)的研究水體中輻射傳輸過程的有效工具, 廣泛應(yīng)用于國際水色遙感領(lǐng)域, 側(cè)重于rs和水下光場模擬以及對遙感反演算法評估等[9-11]。已有的分析表明, HydroLight是進(jìn)行石油類污染水體輻射傳輸特性的研究的有效模型, 但要將HydroLight應(yīng)用于石油類污染水體的輻射傳輸特性的研究, 需要知道水體組分的IOPs隨波長和深度的參數(shù)化模型, 如吸收系數(shù), 散射系數(shù)和散射相函數(shù)[12], 以及外界環(huán)境參數(shù)(風(fēng)速、底部介質(zhì)等)。目前對Ⅱ類水體的水色因子(葉綠素、懸浮物、黃色物質(zhì)等)的分布、吸收系數(shù)和后向散射系數(shù)的研究比較成熟, 并形成了眾多的參數(shù)化模型[13-14], 為石油類物質(zhì)與各水色因子混合后對rs的影響和貢獻(xiàn)的模擬奠定了基礎(chǔ)。

    本文基于在遼寧大連港現(xiàn)場測定的油含量、吸收系數(shù)和后向散射系數(shù)等實(shí)測水體數(shù)據(jù), 結(jié)合HydroLight提供的部分參數(shù), 進(jìn)行油類物質(zhì)和葉綠素單一組分和混合組分的rs光譜模擬, 從理論上探究含油水體對rs光譜變化的影響, 為提高遙感反演石油類濃度模型的精度以及開展石油類污染水體水下光場變化特征研究提供參考依據(jù)。

    1 方法

    1.1 試驗(yàn)區(qū)與觀察時(shí)間

    本文選擇大連港作為研究區(qū)域, 在該海域3個(gè)站點(diǎn)(A、B、C)進(jìn)行了觀測, 站點(diǎn)分布如圖1所示。A站點(diǎn)曾發(fā)生石油輸油管爆炸事件, 有原油的殘留。B站點(diǎn)原油泊位東側(cè), 受到一定程度石油類物質(zhì)污染[15]。C站點(diǎn)受油類污染較小。測量時(shí)間為2018年8月25日—27日, 每天7: 00開始、17: 00結(jié)束, 整點(diǎn)觀測, 每隔1 h觀測1次, 時(shí)間序號分別為1~11, 樣品編號采用站點(diǎn)編號加觀測時(shí)間編號的方式, 例如A1為A站點(diǎn)在7: 00進(jìn)行測管的樣品, B6為B站點(diǎn)在12: 00進(jìn)行觀測的樣品。

    圖1 站點(diǎn)分布圖

    1.2 現(xiàn)場數(shù)據(jù)采集

    1.2.1 油濃度測量

    油濃度的測量采用美國特納的TD-500D熒光測油儀。該設(shè)備的測定原理與分子熒光光度法(《SL366- 2006水質(zhì)石油類的測定分子熒光光度法》[16]相同, 首先使用正己烷將樣品中的油類物質(zhì)萃取出來, 其次使用特定波長的紫外光照射萃取后的油類物質(zhì), 使其釋放出熒光, 最后根據(jù)熒光強(qiáng)度和油濃度在一定范圍內(nèi)接近線性關(guān)系計(jì)算出油濃度。

    1.2.2 遙感反射比rs的測量

    遙感反射比rs的測量運(yùn)用的是水面以上測量法[17], 參考板為反射率約為30%的灰板, 采用的測量儀器為美國ASD公司生產(chǎn)的可見光-近紅外地物光譜儀(ASD-FieldSoec 3 350 nm~2 500 nm)。

    1.2.3 油類物質(zhì)吸收系數(shù)測量

    黃妙芬等[18]將油質(zhì)量濃度小于0.1 mg/L的樣品作為無石油類污染水體, 對它們的黃色物質(zhì)吸收系數(shù)求算術(shù)平均值, 作為本底值, 然后用油質(zhì)量濃度大于1.0 mg/L水樣的黃色物質(zhì)吸收系數(shù)減去本底值, 作為油類物質(zhì)的吸收系數(shù)。作者采用該處理方式, 測量儀器選擇日本日立UV-3900紫外-可見光分光光度計(jì), 測量時(shí)波長設(shè)置為400~700 nm, 在樣品的制備、測量和分析過程都遵循NASA水色要素測量規(guī)范。分別測量無油水樣和有油水樣的吸收系數(shù), 然后將無油水樣作為本底值, 用有油水樣的吸收系數(shù)減去無油水樣的吸收系數(shù), 從而得到油類物質(zhì)的吸收系數(shù)。

    1.2.4 水體后向散射系數(shù)測量

    后向散射系數(shù)的測量采用美國HOBI labs的HydroScat-6(HS-6)儀器, 測量方法按照《海洋光學(xué)調(diào)查技術(shù)規(guī)程》[19]進(jìn)行。為降低海水溫度和氣泡對儀器測量結(jié)果的影響, 將儀器浸沒水中后, 靜置2 min再布放, 然后將HS-6測量的數(shù)據(jù)取表層(0~0.5 m)的平均值作為后向散射系數(shù)。

    1.2.5 油類物質(zhì)的比吸收系數(shù)和比散射系數(shù)的計(jì)算

    油類物質(zhì)吸收系數(shù)可以通過比吸收系數(shù)與油濃度oil計(jì)算得到[20], 關(guān)系如下:

    =*×oil, (1)

    式中,為吸收系數(shù)(m–1),為比吸收系數(shù)(Specific absorption coefficient)(L·mg–1·m–1),oil為油質(zhì)量濃度(mg/L)。根據(jù)公式(1), 可以計(jì)算出*。

    宋慶君等[21]提出的石油類單位后向散射系數(shù)(mass-specific backscattering coefficient)與石油類物質(zhì)濃度之間的擬合關(guān)系, 見公式(2)。根據(jù)公式(2)求出石油類物質(zhì)特定波長上的比散射系數(shù)(Specific scattering coefficient), 接著對其進(jìn)行線性插值, 插值范圍為400~700 nm。

    表1 不同波段石油質(zhì)量比后向散射系數(shù)[21]

    1.3 HydroLight模擬參數(shù)設(shè)置

    本文選擇輻射傳輸模式HydroLight 5.3所提供的“A User Defined Model”的功能, 該功能可實(shí)現(xiàn)對單一水體組分和混合水體的rs進(jìn)行模擬。在本研究中水體組分主要包括純水、油類物質(zhì)和葉綠素。在進(jìn)行模擬之前需要輸入水體相關(guān)組分的濃度、吸收系數(shù)、散射系數(shù)和相函數(shù), 以及外界條件等參數(shù)。

    圖2和表2分別為2018年8月25日~27日在大連港實(shí)測的oil隨深度的變化和3個(gè)站點(diǎn)oil的統(tǒng)計(jì)信息, 分析可見,oil的范圍在0.20~7.80 mg/L, 考慮到自然界多數(shù)量級的變化是按照指數(shù)規(guī)律, 故在配比時(shí),oil的取值設(shè)置為0.1、0.2、0.3、0.5、1.0、1.5、2.0、3.0、5.0、10.0。在HydroLight中, 基于比吸收系數(shù)(比散射系數(shù))為常數(shù)的假設(shè), 利用已知的比吸收系數(shù)(比散射系數(shù))和水體組分濃度來求解出相應(yīng)組分的吸收系數(shù)(散射系數(shù))是獲取組分吸收系數(shù)(散射系數(shù))的方法之一[22]。對于油類物質(zhì)的比吸收系數(shù)和比散射系數(shù)的設(shè)置, 由于不同站點(diǎn)在不同時(shí)刻上兩者變化較大, 如表3所示。因此, 作者以A10站點(diǎn)測量的最大比吸收系數(shù)和比散射系數(shù)作為油物質(zhì)的比吸收系數(shù)和比散射系數(shù), 即0.019 1和0.147 3 (L·mg–1·m–1)。而油類物質(zhì)相函數(shù)則選擇HydroLight提供的典型顆粒相函數(shù), 該函數(shù)適合大部分輻射傳輸運(yùn)算, 是MOBLEY[23]根據(jù)PETZOLD[24]數(shù)據(jù), 提出一種取b/=0.018 3通過代入Fournier and Forand (FF)體散射函數(shù)求解得到隨波長和深度變化的顆粒體散射函數(shù)。

    對于葉綠素的IOPs設(shè)置參考文獻(xiàn)[25], 比吸收系數(shù)和比散射系數(shù)見圖3, 葉綠素濃度chl設(shè)置同oil。另外模擬深度設(shè)為無限深, 具體情況見表4。

    圖2 不同站點(diǎn)油類物質(zhì)濃度隨深度的變化曲線

    表2 不同站點(diǎn)油濃度觀測值的統(tǒng)計(jì)(單位:mg·L–1)

    表3 各站點(diǎn)油物質(zhì)的比吸收系數(shù)與比后向散射系數(shù)統(tǒng)計(jì)(單位:L·mg–1·m–1)

    續(xù)表

    圖3 葉綠素組分比吸收系數(shù)a*和比散射系數(shù)b*

    表4 HydroLight模擬油藻混合水體遙感反射比參數(shù)輸入

    *輻射傳輸模式HydroLight中自帶的數(shù)據(jù)

    2 結(jié)果與討論

    2.1 純油水體的Rrs模擬

    圖4為純油水體模擬的rs。由圖4可知, 純油水體的光譜特征大致可以分為400~480 nm、480~ 700 nm兩個(gè)特征段。在400~480 nm段,rs的大小隨波長λ的增加而減少, 也隨油類物質(zhì)濃度oil的增加而減少, 且rs的差異性逐漸減少, 而rs隨后向散射系數(shù)的增大而增大, 隨吸收系數(shù)的增大而減少, 由于前期oil較低, 主要受純水的影響, 隨著oil的增加, 油類物質(zhì)的吸收特性逐漸體現(xiàn)。而480~700 nm段,rs隨λ的增加而減少, 但隨oil增大而增大, 由于純水在該波段里具有強(qiáng)吸收, 故在此波段范圍, 油的散射特性隨oil的增加逐漸體現(xiàn)。

    2.2 葉綠素水體的Rrs模擬

    葉綠素水體的rs模擬由于葉綠素水體rs模擬結(jié)果在量級上存在差異, 為了更好地描述葉綠素的光譜圖, 分為低葉綠素濃度chl(≤1.0 mg/m3)和高chl(>1.0 mg/m3)的rs光譜, 如圖5a和b所示。圖5a表明低chl的rs光譜特性反映為純水的光譜特性, 葉綠素的光譜特性需要一定濃度才能表現(xiàn)出來。由圖5b可知, 在高chl下, 由于受葉綠素在紫藍(lán)波段的強(qiáng)吸收影響, 使得光譜在440 nm和675 nm波長附近存在谷值; 在510~610 nm波長范圍存在反射峰, 并隨著濃度升高峰值逐漸增大, 這是由葉綠素的弱吸收所引起的; 在675 nm波長之后, 由葉綠素?zé)晒馓匦运碌诙瓷浞宀惶黠@, 這是由于葉綠素濃度設(shè)置較低所致。

    圖4 純油水體Rrs的模擬

    圖5 葉綠素水體Rrs的模擬

    2.3 油藻混合水體Rrs模擬

    為了弄清葉綠素和油類物質(zhì)兩者在混合水體中對rs的影響, 分別模擬了混合水體rs隨葉綠素和油類物質(zhì)濃度的變化曲線。

    圖6為oil一定和低chl情況下,rs隨chl的變化的曲線圖。圖6a和6b表明, 隨chl的增大, 即使oil較大時(shí),rs光譜形狀特性逐漸反映為葉綠素的光譜特性。由于低chl下, 水體光譜特性并不體現(xiàn)為葉綠素的光譜特性, 葉綠素的光譜特性要在一定chl才能表現(xiàn)出來, 故以下分析均在chl>1.0 mg/m3下進(jìn)行分析。

    圖7為特定oil下, 油藻混合水體的rs(實(shí)線)與純藻水體rs(虛線)的比較圖, 結(jié)果表明, 油類物質(zhì)的存在, 只增大了葉綠素的rs量值, 并不改變其隨波長的變化。另外, 油類物質(zhì)濃度oil對葉綠素rs光譜的影響主要體現(xiàn)在葉綠素的峰值區(qū)(500~600 nm)與熒光區(qū)(670~700 nm), 純藻水體rs隨chl的增大而增大, 而混合水體rs的量值隨chl的增大而下降。這是由于油類物質(zhì)的“阻礙”作用, 削弱了誘發(fā)葉綠素發(fā)生反射和熒光的能量, 使得葉綠素的反射峰值與熒光峰值隨chl的增大而減少, 但不影響葉綠素的光譜形狀。

    圖6 油濃度固定, 在低葉綠素下, 水體Rrs隨葉綠素濃度的變化

    圖7 油濃度一定下, 油藻混合水體Rrs光譜與純藻水體的比較

    圖8為chl分別為5.0 mg/m3和10.0 mg/m3情況下rs隨oil的變化圖。圖8a和8b的結(jié)果表明: 隨oil增大,rs逐漸增大但隨波長λ的變化趨勢沒有明顯改變, 即chl對水體rs光譜形狀的影響是強(qiáng)于oil的影響。綜上, 油類物質(zhì)的存在, 改變了油藻混合水體的光譜量值, 但對光譜形狀影響甚小, 即油藻混合水體的光譜形狀主要受葉綠素的影響。

    圖8 葉綠素濃度一定下, 油藻混合水體的Rrs隨油濃度的變化

    2.4 模擬光譜與實(shí)測光譜的對比分析

    對A、B和C站點(diǎn)相應(yīng)實(shí)測的rs光譜進(jìn)行線性插值和平滑(移動(dòng)平均法)處理后, 結(jié)果如圖9所示。分析圖9可知,rs光譜普遍存在2個(gè)峰值, 前鋒大于后峰, 并且前峰值區(qū)位于570 nm附近, 后峰峰值區(qū)位于680 nm附近。A站點(diǎn)的rs在400~570 nm波長范圍內(nèi)隨波長的增大而增大, 而B和C站點(diǎn)在該范圍內(nèi)除了隨波長增大rs增大的趨勢外, 還存在光譜值在波長400~440 nm先減少在波長440~570 nm再增大的趨勢, 如B2、C4的光譜曲線。

    本文模擬的rs隨波長的變化趨勢類似于B, C站點(diǎn)的rs隨波長的變化趨勢, 在光譜形狀上具有一定的相似性, 但與實(shí)測光譜的差異主要體現(xiàn)在: (1) 模擬光譜第二峰不明顯。由上述的分析可知, 該峰主要受葉綠素?zé)晒獾挠绊? 可能是由于葉綠素組分濃度設(shè)置較小所致; (2) 光譜量值上的不同。實(shí)測光譜測量的rs是主要由葉綠素, 黃色物質(zhì)和懸浮泥沙的貢獻(xiàn), 在含油水體中還包括油物質(zhì)的貢獻(xiàn)。而本模擬僅考慮油與葉綠素兩組分, 并且外部條件理想的情況下的混合光譜。實(shí)測光譜的測量環(huán)境較模擬要復(fù)雜, 太陽的方位與云的覆蓋、風(fēng)速、潮流和生物降解等變化會(huì)使得水下光場發(fā)生動(dòng)態(tài)改變。另外, RUD?等[27]認(rèn)為輻射傳輸理論能很好地評價(jià)油物質(zhì)對rs的貢獻(xiàn)和模擬的準(zhǔn)確性與散射相函數(shù)的選擇有強(qiáng)相關(guān)性。

    2.5 油藻光譜特征對水色遙感算法的影響分析

    在水色遙感中, 遙感反射比rs, 定義為w與d的比值, 其表達(dá)式為,

    式中, Ed為太陽在水表面的下行輻照度(W·m–2·nm–1), Lw為離水輻射(water-leaving radiance)(W·m–2·nm–1·sr–1), 是指經(jīng)水/氣(water-air)界面反射和透射后的輻射能。水體中的各個(gè)重要光學(xué)成分濃度發(fā)生變化時(shí), 必然引起水體光學(xué)性質(zhì)的變化, 主要表現(xiàn)為水體的吸收和散射信號的變化, 綜合起來看, 即是Rrs的變化。通過衛(wèi)星傳感器接受信號的變化, 可以反演得到水體中的各個(gè)光學(xué)重要成分, 即水體中的浮游植物、懸浮泥沙、黃色物質(zhì)以及其他物質(zhì)的含量。

    另外, 在水色遙感中,rs是將表觀光學(xué)量與固有光學(xué)量連接起來的一個(gè)橋梁, 見式(4), 可將rs與水體的總吸收系數(shù)和總的后向散射系數(shù)b聯(lián)系起來, 從而得到反演水體石油類污染的生物-光學(xué)模型,

    其中,為經(jīng)驗(yàn)參數(shù)約為0.32~0.33,和都是太陽天頂角的函數(shù),為水氣界面的投射系數(shù),是水的折射指數(shù)。

    在含油水體中, 水體的主要組分主要包括石油類污染、懸浮物和葉綠素, 遙感得到的rs是這些要素共同貢獻(xiàn)的結(jié)果。我們的研究目標(biāo)是尋找石油類污染的遙感探測特征。在這種情況下, 顯然懸浮物和葉綠素對光譜曲線的影響成為一種噪聲, 所以我們首先要了解清楚僅含油水體的波譜特征以及疊加了油物質(zhì)后, 葉綠素和懸浮物對水體波譜曲線的影響, 然后才能確定石油類污染的水體光譜特征。由于篇幅所限, 本文僅討論了油藻混合的情況。根據(jù)“2.3油藻混合水體rs模擬”小節(jié)模擬結(jié)果, 發(fā)現(xiàn)在油藻混合水體中, 水體的光譜形狀主要受葉綠素的影響, 油類物質(zhì)的存在只改變光譜的量值。這表明藻類作為一種顆粒物, 油類物質(zhì)會(huì)吸附在其上, 共同影響遙感反射比。同時(shí), 由于油類物質(zhì)的存在只改變了混合水體的rs的量值, 這表明在利用半分析模型反演水體組分濃度時(shí), 直接把油類物質(zhì)作為一個(gè)水色要素進(jìn)行反演是可行的。

    3 結(jié)語

    本文基于大連港實(shí)測數(shù)據(jù)分離出油類物質(zhì)的濃度和固有光學(xué)參數(shù), 通過輻射傳輸模式HydroLight模擬出油類物質(zhì), 葉綠素和混合組分的水體rs光譜進(jìn)行理論分析。模擬結(jié)果表明: (1) 在純油水體中,rs隨波長增加而下降, 隨油濃度的變化有兩個(gè)特征波段: 在400~480 nm波長范圍中, 隨油濃度的增加而減少, 在480~700 nm處, 反而隨油濃度的增加而增加; (2) 純?nèi)~綠素水體中, 葉綠素的光譜特征需要在一定葉綠素濃度下才顯現(xiàn), 在低葉綠素濃度下水體光譜反映為純水的光譜特征; (3) 在油和藻混合水體中, 油類物質(zhì)的存在使得葉綠素的rs的增強(qiáng), 但不改變r(jià)s隨波長的變化趨勢, 這說明混合水體光譜形狀主要受葉綠素的影響, 油類物質(zhì)對光譜形狀的影響較小。

    本文所得結(jié)論是在理想情況和固定相關(guān)參數(shù)下進(jìn)行的, 模擬的光譜形狀與實(shí)測光譜形狀具有一定的相似性, 但也有所差異, 主要體現(xiàn)在: (1) 模擬僅涉及葉綠素與油物質(zhì), 但實(shí)測光譜水體成分和外部條件復(fù)雜; (2) 風(fēng)速與潮流對油類物質(zhì)固有光學(xué)的測量有一定的影響; (3) 對于油類物質(zhì)的比吸收系數(shù)和比散射系數(shù)的計(jì)算使用的是經(jīng)驗(yàn)算法并采用固定值, 其適用性以及精確度需要更多數(shù)據(jù)來驗(yàn)證。本文對光譜曲線模擬的結(jié)果在一定程度上解釋了油類物質(zhì)對水遙感反射比的影響, 但仍需要更多油類物質(zhì)的固有光學(xué)特性參數(shù)數(shù)據(jù)進(jìn)行模擬驗(yàn)證, 這是下一步研究的重點(diǎn)。

    [1] 邢旭峰, 王忠林, 黃妙芬, 等.石油污染水體油濃度及后向散射系數(shù)垂向變化特性分析[J].廣東海洋大學(xué)學(xué)報(bào), 2020, 40(4): 75-81.

    XING Xufeng, WANG Zhonglin, HUANG Miaofen, et al.Vertical variation characteristics of petroleum concentration and backscattering coefficient in petroleum-polluted water body[J].Journal of Guangdong Ocean University, 2020, 40(4): 75-81.

    [2] 黃妙芬, 莊煬, 劉楊, 等.大連港海域石油類物質(zhì)比吸收系數(shù)特征分析[J].海洋湖沼通報(bào), 2020(3): 141- 149.

    HUANG Miaofen, ZHUANG Yang, LIU Yang, et al.Analysis on Characteristics of specific absorption coefficient of petroleum-polluted water in Dalian port[J].Transactions of Oceanology and Limnology, 2020(3): 141-149.

    [3] HUANG M F, XING X F, SONG Q J, et al.A new algorithm of retrieving a petroleum substances absorption coefficient in sea water based on a remote sensing image[J].Acta Oceanologica Sinica, 2016, 35(11): 97-104.

    [4] 黃妙芬, 宋慶君, 陳利搏, 等.基于歸一化遙感反射比反演水中石油含量模式研究[J].海洋技術(shù)學(xué)報(bào), 2015, 34(1): 1-9.

    HUANG Miaofen, SONG Qingjun, CHEN Libo, et al.Research on an inverse model of petroleum content in water bodies based on the normalized remote sensing reflectance[J].Journal of Ocean Technology, 2015, 34(1): 1-9.

    [5] 黃妙芬, 宋慶君, 邢旭峰, 等.石油污染水體熒光圖譜特征分析[J].光譜學(xué)與光譜分析, 2014, 34(9): 2466- 2471.

    HUANG Miaofen, SONG Qingjun, XING Xufeng, et al.Analysis of fluorescence spectrum of petroleum-polluted water[J].Spectroscopy and Spectral Analysis, 2014, 34(9): 2466-2471.

    [6] HUANG M F, SONG Q J, XING X F, et al.A biooptical model of retrieving petroleum concentration in seawater[J].Acta Oceanologica Sinica, 2014, 33(5): 81-85.

    [7] 黃妙芬, 邢旭峰, 宋慶君, 等.油砂混合水體后向散射系數(shù)光譜貢獻(xiàn)分離算法I: 實(shí)驗(yàn)理論[J].光譜學(xué)與光譜分析, 2017, 37(1): 205-211.

    HUANG Miaofen, XING Xufeng, SONG Qingjun, et al.New Algorithms to separate the contribution of petroleum substances and suspended particulate matter on the scattering coefficient spectrum from mixed water[J].Spectroscopy and Spectral Analysis, 2017, 37(1): 205- 211.

    [8] GITELSON A.The peak near 700 nm on radiance spectra of algae and water: Relationships of its magnitude and position with chlorophyll concentration[J].International Journal of Remote Sensing, 1992, 13(14): 3367-3373.

    [9] MA C Y, XU S, WANG H S, et al.A real-time photo- realistic rendering algorithm of ocean color based on bio-optical model[J].Journal of Ocean University of China, 2016, 15(6): 996-1006.

    [10] SUNDARABALAN V B, SHANMUGAM P, AHN Y H.Modeling the underwater light field fluctuations in coastal oceanic waters: Validation with experimental data[J].Ocean Science Journal, 2016, 51(1): 67-86.

    [11] PAHLEVAN N, SCHOTT J R, FRANZ B A, et al.Landsat 8 remote sensing reflectance (rs) products: Evaluations, intercomparisons, and enhancements[J].Remote Sens Environ, 2017, 190: 289-301.

    [12] 黃妙芬, 駱蔚健, 劉楊, 等.應(yīng)用Hydrolight模式研究含油水體輻射傳輸特性的關(guān)鍵技術(shù)[J].海洋技術(shù)學(xué)報(bào), 2019, 38(3): 7-14.

    HUANG Miaofen, LUO Weijian, LIU Yang, et al.The key technology for studying the radiation transmission characteristics of petroleum-polluted water using the Hydrolight mode[J].Journal of Ocean Technology, 2019, 38(3): 7-14.

    [13] CHARLES L G.Partitioning spectral absorption in case 2 waters: discrimination of dissolved and particulate components[J].Applied Optics, 2002, 21(41): 4220-4233.

    [14] 宋慶君, 唐軍武.黃海、東海海區(qū)水體散射特性研究[J].海洋學(xué)報(bào)(中文版), 2006, 28(4): 56-63.

    SONG Qingjun, TANG Junwu.The study on the scatteringproperties in the Huanghai sea and east China sea[J].Acta Oceanologica Sinica, 2006, 28(4): 56-63.

    [15] 黃妙芬, 劉楊, 邢旭峰, 等.大連港海域油含量時(shí)空變化特征分析[J].大連海洋大學(xué)學(xué)報(bào), 2020, 35(2): 273-279.

    HUANG Miaofen, LIU Yang, XING Xufeng, et al.Analysis of temporal and spatial variation of petroleum content in Dalian port[J].Journal of Dalian Ocean University, 2020, 35(2): 273-279.

    [16] SL366-2006.水質(zhì)石油類的測定分子熒光光度法[S].

    SL366-2006.Water quality – Determination of petroleum oil – moleclar fluorescence spectrometric method[S].

    [17] 唐軍武, 田國良, 汪小勇, 等.水體光譜測量與分析Ⅰ: 水面以上測量法[J].遙感學(xué)報(bào), 2004, 8(1): 37-44.

    TANG Junwu, TIAN Guoliang, WANG Xiaoyong, et al.The methods of water spectra measurement and analysisⅠ: Above-Water method[J].Journal of Remote Sensing, 2004, 8(1): 37-44.

    [18] 黃妙芬, 唐軍武, 宋慶君.石油類污染水體吸收特性分析[J].遙感學(xué)報(bào), 2010, 14(1): 131-147.

    HUANG Miaofen, TANG Junwu, SONG Qingjun.Analysis of petroleum-polluted water absorption spectral properties[J].Journal of Remote Sensing, 2010, 14(1): 131-147.

    [19] 國家海洋局908專項(xiàng)辦公室.海洋光學(xué)調(diào)查技術(shù)規(guī)程[M].北京: 海洋出版社, 2006: 8-11.

    908 special office of State Oceanic Administration.The technical regulations for oceanographic optical survey[M].Beijing: China Ocean Press, 2006: 8-11.

    [20] BRICAUD A, BABIN M, MOREL A, et al.Variability in the chlorophyll-specific absorption coefficients of natural phytoplankton: Analysis and parameterization[J].Journal of Geophysical Research, 1995, 100(C7): 13321- 13332.

    [21] 宋慶君, 黃妙芬, 唐軍武, 等.水中石油類含量對后向散射系數(shù)光譜的影響[J].光譜學(xué)與光譜分析, 2010, 30(9): 2438-2442.

    SONG Qingjun, HUANG Miaofen, TANG Junwu, et al.Influence of petroleum concentration in water on spectral backscattering coefficient[J].Spectroscopy and Spectral Analysis, 2010, 30(9): 2438-2442.

    [22] 黃佳慧, 段洪濤, 馬榮華, 等.基于Hydrolight的太湖水體區(qū)域化固有光學(xué)模型參數(shù)率定與驗(yàn)證[J].中國科學(xué)院大學(xué)學(xué)報(bào), 2014, 31(5): 613-625.

    HUANG Jiahui, DUAN Hongtao, MA Ronghua, et al.Study on regional parameters of radiative transfer simulation based on Hydrolight in Taihu Lake[J].Journal of University of Chinese Academy of Sciences, 2014, 31(5): 613-625.

    [23] MOBLEY C, GENTILI B, GORDON H, et al.Comparison of numerical models for computing underwater light fields[J].Applied Optics, 1993, 32: 7484-7504.

    [24] PETZOLD T J.Volume Scattering Functions for Selected Ocean Waters[R].San Diego: Scripps Institute of Oceanography, 1972: 72-78.

    [25] 李云梅, 王橋, 黃家柱, 等.太湖水體光學(xué)特性及水色遙感[M].北京: 科學(xué)出版社, 2010: 129.

    LI Yunmei, WANG Qiao, HUANG Jiazhu, et al.Optical properties and Remote sensing of water color of Taihu Lake[M].Beijing: Science Press, 2010: 129.

    [26] POPE R M, FRY E S.Absorption spectrum (380–700 nm) of pure water.II.Integrating cavity measurements[J].Applied Optics, 1997, 36(33): 8710-8723.

    [27] RUD? K, DARECKI M, TOCZEK H.Modelling the influence of oil content on optical properties of sea-water in the Baltic Sea [J].Journal of the European Optical Society-Rapid Publications, 2013, 8: 13063.

    Remote sensing reflectance spectral analysis of oil and algae mixed water

    LUO Wei-jian1, HUANG Miao-fen2, YANG Jun-jie3, WANG Xin2

    (1.College of Oceanography and Meteorology, Guangdong Ocean University, Zhanjiang 524088, Guangdong, China; 2.Faculty of Mathematics and Computer Science, Guangdong Ocean University, Zhanjiang 524088, Guangdong, China; 3.College of Electronic and Information Engineering, Guangdong Ocean University, Zhanjiang 524088, Guangdong, China)

    remote sensing reflectance; HydroLight; oil; chlorophyll

    The radiation transfer model HydroLight is an effective tool to study the radiation transfer characteristics of water body, and it is also an effective model to study the radiation transfer characteristics of petroleum water body.Based on the apparent and inherent optical data of oil-contaminated water bodies measured in the sea area of Dalian Port, Liaoning, in August 2018, we employ the radiation transfer model HydroLight to simulate the remote sensing reflectance (rs)of the single component and mixed water of oil and algae by the method of different concentrations of proportioning.The results show that: (1) In a water body containing only oil,rscan be divided into two characteristic bands: 400–480 nm and 480–700 nm: in the range of 400–480 nm, it decreases with the increasing oil concentrationoil, while in the range of 480–700 nm, it increases; (2) In single-component water containing only algae, the spectral characteristics of chlorophyll require its concentrationchlto reach a certain value to show up; (3) In mixed water body including oil and algae, the increase ofoilwill increasersof chlorophyll, but it will not alter the trend of thersof chlorophyll with wavelength.This demonstrates that chlorophyll has a major influence on the spectral shape of oil algae mixed water, while the presence of oil only modifies the spectral value.The use of HydroLight in researching thersof oil-contaminated water contributes to ocean color remote sensing research.It has substantial theoretical significance for improving the bio-optical model of various water bodies.

    [National Natural Science Foundation of China, No.41806186; Guangdong Ocean University’s “Innovation and Strong School Project” Independent Innovation Capability Enhancement Project in 2017, No.GDOU2017052501; Guangdong Ocean University Scientific Research Startup Funding Project, No.E16187]

    Oct.30, 2020

    P733.3

    A

    1000-3096(2021)10-0049-12

    10.11759/hykx20201030002

    2020-10-30

    2021-02-12

    國家自然科學(xué)基金項(xiàng)目(41771384); 廣東海洋大學(xué)2017年“創(chuàng)新強(qiáng)校工程”自主創(chuàng)新能力提升項(xiàng)目(GDOU2017052501); 廣東海洋大學(xué)科研啟動(dòng)經(jīng)費(fèi)資助項(xiàng)目(E16187)

    駱蔚健(1995—)男, 廣東惠州人, 碩士研究生, 主要從事海洋水色遙感研究, E-mail: 1136953630@qq.com; 黃妙芬(1963—), 通信作者, E-mail: hmf808@163.com

    (本文編輯: 譚雪靜)

    亚洲av成人不卡在线观看播放网| av片东京热男人的天堂| 国产在线视频一区二区| 久热这里只有精品99| 97人妻天天添夜夜摸| 午夜久久久在线观看| 母亲3免费完整高清在线观看| 99re6热这里在线精品视频| 精品亚洲乱码少妇综合久久| 亚洲国产成人一精品久久久| 99国产精品99久久久久| 欧美日韩亚洲国产一区二区在线观看 | 满18在线观看网站| 亚洲精品久久成人aⅴ小说| av免费在线观看网站| 国产成人免费无遮挡视频| 一本色道久久久久久精品综合| 桃红色精品国产亚洲av| 亚洲性夜色夜夜综合| 国产亚洲午夜精品一区二区久久| 纵有疾风起免费观看全集完整版| 欧美午夜高清在线| 99国产极品粉嫩在线观看| 亚洲专区国产一区二区| xxxhd国产人妻xxx| 男人操女人黄网站| 国内毛片毛片毛片毛片毛片| 国产精品久久久人人做人人爽| 激情视频va一区二区三区| 亚洲自偷自拍图片 自拍| 黑人欧美特级aaaaaa片| 首页视频小说图片口味搜索| 国产精品成人在线| 黄频高清免费视频| 精品亚洲成a人片在线观看| 中文亚洲av片在线观看爽 | 国产在线精品亚洲第一网站| aaaaa片日本免费| 视频在线观看一区二区三区| 国产亚洲欧美在线一区二区| 日本欧美视频一区| 狠狠精品人妻久久久久久综合| 一区二区三区乱码不卡18| 制服人妻中文乱码| 三级毛片av免费| 大香蕉久久成人网| 99国产综合亚洲精品| 久久久久精品人妻al黑| 国产高清激情床上av| 免费看a级黄色片| 侵犯人妻中文字幕一二三四区| 一级毛片精品| 色婷婷av一区二区三区视频| 午夜精品久久久久久毛片777| 麻豆av在线久日| 大陆偷拍与自拍| 欧美日本中文国产一区发布| 久热这里只有精品99| 国产成人系列免费观看| 免费在线观看影片大全网站| 亚洲精品成人av观看孕妇| 高清欧美精品videossex| 99国产精品99久久久久| 高潮久久久久久久久久久不卡| 免费久久久久久久精品成人欧美视频| 国产成人影院久久av| 少妇裸体淫交视频免费看高清 | 欧美日韩亚洲高清精品| 老熟妇乱子伦视频在线观看| 一区二区三区国产精品乱码| 国产精品久久电影中文字幕 | 亚洲精品粉嫩美女一区| 国产高清激情床上av| 午夜福利视频在线观看免费| 精品国产亚洲在线| 超碰成人久久| 久久青草综合色| √禁漫天堂资源中文www| 日韩一卡2卡3卡4卡2021年| tocl精华| 国产精品 国内视频| 中文字幕另类日韩欧美亚洲嫩草| 美女高潮喷水抽搐中文字幕| 久久久国产欧美日韩av| 蜜桃国产av成人99| 在线观看免费视频日本深夜| 日韩免费高清中文字幕av| 看免费av毛片| 国产高清国产精品国产三级| 久久久久视频综合| 69精品国产乱码久久久| 1024视频免费在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文av在线| 免费在线观看完整版高清| 亚洲成av片中文字幕在线观看| 在线天堂中文资源库| 国产淫语在线视频| 亚洲午夜理论影院| 成人三级做爰电影| 免费不卡黄色视频| 日韩欧美免费精品| 他把我摸到了高潮在线观看 | 黑人猛操日本美女一级片| 最黄视频免费看| 亚洲精品中文字幕一二三四区 | 精品国产一区二区久久| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲熟女精品中文字幕| 免费高清在线观看日韩| 岛国毛片在线播放| 久久中文看片网| 搡老熟女国产l中国老女人| 久久久久视频综合| 啪啪无遮挡十八禁网站| 少妇被粗大的猛进出69影院| 国产精品电影一区二区三区 | 亚洲色图 男人天堂 中文字幕| 免费观看人在逋| 1024香蕉在线观看| 男人舔女人的私密视频| 欧美日韩av久久| 人人澡人人妻人| 一级毛片女人18水好多| 国产片内射在线| 一区二区av电影网| 麻豆av在线久日| 国产免费福利视频在线观看| 老汉色av国产亚洲站长工具| 国产不卡一卡二| 精品亚洲乱码少妇综合久久| 丰满饥渴人妻一区二区三| 欧美激情久久久久久爽电影 | 国产精品欧美亚洲77777| 国产成人免费观看mmmm| 国产熟女午夜一区二区三区| 国产一区二区在线观看av| 亚洲欧洲日产国产| 欧美黄色淫秽网站| 国产伦人伦偷精品视频| 亚洲精品粉嫩美女一区| 国产91精品成人一区二区三区 | 欧美日韩福利视频一区二区| 国产免费福利视频在线观看| 午夜福利,免费看| 日韩视频在线欧美| 涩涩av久久男人的天堂| 999久久久国产精品视频| 国产在线一区二区三区精| 国产亚洲av高清不卡| 成年人午夜在线观看视频| 精品卡一卡二卡四卡免费| 亚洲欧美日韩高清在线视频 | 亚洲第一av免费看| 亚洲精品在线观看二区| 亚洲国产欧美日韩在线播放| 9热在线视频观看99| 亚洲成人手机| 黄色毛片三级朝国网站| av一本久久久久| 国产麻豆69| 国产av又大| 天堂动漫精品| 少妇精品久久久久久久| 成年版毛片免费区| 亚洲中文日韩欧美视频| 777米奇影视久久| 国产色视频综合| 精品国产乱码久久久久久男人| 老汉色av国产亚洲站长工具| 欧美 日韩 精品 国产| 又紧又爽又黄一区二区| 一本综合久久免费| 免费在线观看日本一区| 国产一区二区在线观看av| 久久精品91无色码中文字幕| 亚洲成av片中文字幕在线观看| 久久婷婷成人综合色麻豆| 国产精品久久久久久精品电影小说| av天堂久久9| 十分钟在线观看高清视频www| 热99re8久久精品国产| 久久久国产成人免费| 久久午夜亚洲精品久久| 别揉我奶头~嗯~啊~动态视频| 中文字幕人妻丝袜制服| 美女高潮到喷水免费观看| 日韩欧美三级三区| 免费在线观看完整版高清| 久久中文看片网| 女人高潮潮喷娇喘18禁视频| 一区二区三区激情视频| 精品一区二区三区视频在线观看免费 | 丁香六月欧美| 我的亚洲天堂| 精品第一国产精品| 淫妇啪啪啪对白视频| 精品一区二区三区av网在线观看 | 女性被躁到高潮视频| 欧美性长视频在线观看| 无限看片的www在线观看| av不卡在线播放| 美女国产高潮福利片在线看| 99国产精品一区二区三区| 女性被躁到高潮视频| 9热在线视频观看99| 国产精品99久久99久久久不卡| 日韩免费高清中文字幕av| 精品久久久精品久久久| 天天躁日日躁夜夜躁夜夜| 亚洲午夜理论影院| 高清欧美精品videossex| 国产精品熟女久久久久浪| 我要看黄色一级片免费的| 一本大道久久a久久精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一区中文字幕在线| av免费在线观看网站| 亚洲成人免费电影在线观看| 国产亚洲一区二区精品| 悠悠久久av| 亚洲国产成人一精品久久久| 亚洲少妇的诱惑av| 天天影视国产精品| 成年动漫av网址| 久久天堂一区二区三区四区| 高清视频免费观看一区二区| 亚洲男人天堂网一区| 日本黄色日本黄色录像| 国产午夜精品久久久久久| 视频区欧美日本亚洲| 少妇裸体淫交视频免费看高清 | 老汉色∧v一级毛片| 免费在线观看影片大全网站| 飞空精品影院首页| 18禁观看日本| 久久精品国产a三级三级三级| 一个人免费看片子| 在线观看人妻少妇| 不卡一级毛片| 日本五十路高清| 少妇精品久久久久久久| 国产高清激情床上av| 久久毛片免费看一区二区三区| 国产精品免费一区二区三区在线 | 国产精品98久久久久久宅男小说| 日韩中文字幕欧美一区二区| 亚洲精品国产区一区二| 999久久久国产精品视频| 极品少妇高潮喷水抽搐| avwww免费| 在线永久观看黄色视频| 亚洲久久久国产精品| 五月开心婷婷网| 在线看a的网站| 91av网站免费观看| 亚洲国产精品一区二区三区在线| 交换朋友夫妻互换小说| 国产成人影院久久av| av福利片在线| 王馨瑶露胸无遮挡在线观看| 免费黄频网站在线观看国产| 一本一本久久a久久精品综合妖精| 美女高潮喷水抽搐中文字幕| 欧美激情极品国产一区二区三区| 丝袜美腿诱惑在线| 欧美在线一区亚洲| 欧美成狂野欧美在线观看| 一本色道久久久久久精品综合| 9热在线视频观看99| 免费在线观看影片大全网站| 涩涩av久久男人的天堂| 一区福利在线观看| 性高湖久久久久久久久免费观看| 母亲3免费完整高清在线观看| 女警被强在线播放| 在线播放国产精品三级| 欧美日韩国产mv在线观看视频| 在线观看66精品国产| 日韩成人在线观看一区二区三区| 亚洲精品美女久久久久99蜜臀| 国产成人精品久久二区二区91| 日韩有码中文字幕| 久久久久久亚洲精品国产蜜桃av| 视频区图区小说| 极品人妻少妇av视频| 日韩大片免费观看网站| 9191精品国产免费久久| 一本大道久久a久久精品| 精品国产超薄肉色丝袜足j| 免费少妇av软件| 国产免费现黄频在线看| 青草久久国产| 久久久国产一区二区| netflix在线观看网站| 999精品在线视频| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜制服| 视频区图区小说| 色视频在线一区二区三区| 国产精品久久久久久精品电影小说| 久久精品国产a三级三级三级| 深夜精品福利| 夜夜爽天天搞| 国产一区二区在线观看av| 精品国产乱子伦一区二区三区| 欧美日韩av久久| 国产在线观看jvid| 女性生殖器流出的白浆| 制服诱惑二区| 亚洲avbb在线观看| 又黄又粗又硬又大视频| 国产欧美日韩一区二区三区在线| 热99久久久久精品小说推荐| 91av网站免费观看| 国产成人av教育| 50天的宝宝边吃奶边哭怎么回事| 黄色成人免费大全| 国产又色又爽无遮挡免费看| 亚洲精品粉嫩美女一区| 中亚洲国语对白在线视频| 老司机福利观看| 久久 成人 亚洲| 亚洲精品一卡2卡三卡4卡5卡| 亚洲专区中文字幕在线| 亚洲精品一二三| 性色av乱码一区二区三区2| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品欧美亚洲77777| 精品久久久精品久久久| 老司机靠b影院| 久热这里只有精品99| 亚洲精华国产精华精| 久久ye,这里只有精品| 成年人午夜在线观看视频| 欧美av亚洲av综合av国产av| 桃红色精品国产亚洲av| 99九九在线精品视频| 亚洲国产毛片av蜜桃av| 亚洲人成电影免费在线| a级毛片黄视频| 日韩三级视频一区二区三区| 亚洲人成电影免费在线| 99国产精品一区二区三区| 亚洲五月色婷婷综合| 亚洲九九香蕉| 正在播放国产对白刺激| 久久精品亚洲熟妇少妇任你| 国产精品影院久久| 亚洲精品美女久久久久99蜜臀| 男女边摸边吃奶| 男女下面插进去视频免费观看| 欧美久久黑人一区二区| 亚洲成人免费av在线播放| 肉色欧美久久久久久久蜜桃| 久久精品国产99精品国产亚洲性色 | 午夜福利免费观看在线| 一进一出抽搐动态| 在线十欧美十亚洲十日本专区| 少妇精品久久久久久久| 亚洲人成电影观看| av电影中文网址| 欧美激情 高清一区二区三区| 69av精品久久久久久 | 搡老岳熟女国产| 亚洲人成电影观看| 一区二区三区国产精品乱码| 亚洲国产毛片av蜜桃av| 亚洲av第一区精品v没综合| 午夜福利在线观看吧| 久久精品成人免费网站| 国产成人欧美在线观看 | 一区二区三区激情视频| 国产在线视频一区二区| 中亚洲国语对白在线视频| 亚洲午夜理论影院| 欧美激情 高清一区二区三区| 精品久久久久久电影网| 色婷婷久久久亚洲欧美| 天堂动漫精品| 人妻久久中文字幕网| 国产精品亚洲一级av第二区| 69av精品久久久久久 | 国产亚洲欧美在线一区二区| 欧美久久黑人一区二区| 国产精品久久久av美女十八| 极品教师在线免费播放| 国产区一区二久久| 国产成人欧美| cao死你这个sao货| 一区二区三区乱码不卡18| 黄色视频不卡| 制服人妻中文乱码| 久久久久久久国产电影| 捣出白浆h1v1| 波多野结衣一区麻豆| 99精品在免费线老司机午夜| 少妇粗大呻吟视频| 黄片小视频在线播放| 欧美乱妇无乱码| 欧美成人午夜精品| 在线永久观看黄色视频| 麻豆av在线久日| svipshipincom国产片| 免费日韩欧美在线观看| 在线观看一区二区三区激情| 99精国产麻豆久久婷婷| 十分钟在线观看高清视频www| 嫩草影视91久久| 亚洲人成77777在线视频| 午夜老司机福利片| 亚洲国产欧美网| 亚洲精品在线观看二区| av有码第一页| 一区二区av电影网| 天天躁日日躁夜夜躁夜夜| 日本vs欧美在线观看视频| 真人做人爱边吃奶动态| 人人妻,人人澡人人爽秒播| 香蕉丝袜av| 欧美国产精品一级二级三级| 在线观看免费高清a一片| 国产日韩欧美亚洲二区| 18在线观看网站| 午夜91福利影院| 变态另类成人亚洲欧美熟女 | 国产野战对白在线观看| 满18在线观看网站| 成人影院久久| 国产成人欧美| 色94色欧美一区二区| 久久久久久亚洲精品国产蜜桃av| 91九色精品人成在线观看| 十八禁高潮呻吟视频| 狠狠婷婷综合久久久久久88av| 国产激情久久老熟女| 999久久久国产精品视频| 免费高清在线观看日韩| 欧美一级毛片孕妇| 又紧又爽又黄一区二区| www.999成人在线观看| 久久久国产成人免费| 在线观看舔阴道视频| 日韩 欧美 亚洲 中文字幕| 18禁黄网站禁片午夜丰满| 女人爽到高潮嗷嗷叫在线视频| 久久久国产欧美日韩av| 久久精品人人爽人人爽视色| 日韩中文字幕视频在线看片| 色婷婷久久久亚洲欧美| 久久久欧美国产精品| 十八禁人妻一区二区| 一本大道久久a久久精品| 女人久久www免费人成看片| 夜夜夜夜夜久久久久| 亚洲欧美一区二区三区久久| 亚洲精品成人av观看孕妇| 国产精品二区激情视频| 在线观看免费午夜福利视频| 色播在线永久视频| 深夜精品福利| 极品人妻少妇av视频| 下体分泌物呈黄色| 久久av网站| 国产精品秋霞免费鲁丝片| 女性被躁到高潮视频| 人人妻人人爽人人添夜夜欢视频| 视频区图区小说| 精品国产超薄肉色丝袜足j| 俄罗斯特黄特色一大片| av一本久久久久| 女人被躁到高潮嗷嗷叫费观| 亚洲九九香蕉| 成人18禁高潮啪啪吃奶动态图| 自线自在国产av| 欧美精品一区二区大全| 看免费av毛片| 日本精品一区二区三区蜜桃| 十分钟在线观看高清视频www| a级毛片黄视频| 一二三四社区在线视频社区8| 日本撒尿小便嘘嘘汇集6| 一级黄色大片毛片| 大香蕉久久网| 国产激情久久老熟女| 国产欧美日韩一区二区三区在线| 精品熟女少妇八av免费久了| 高清欧美精品videossex| 国产精品免费一区二区三区在线 | 国产精品二区激情视频| 国产成人一区二区三区免费视频网站| 欧美人与性动交α欧美软件| 亚洲精品美女久久av网站| 国产人伦9x9x在线观看| 久久人人97超碰香蕉20202| 母亲3免费完整高清在线观看| 日日摸夜夜添夜夜添小说| 成年动漫av网址| 欧美日韩福利视频一区二区| 757午夜福利合集在线观看| 人人妻人人爽人人添夜夜欢视频| 美女扒开内裤让男人捅视频| www.自偷自拍.com| 免费看十八禁软件| 亚洲欧美一区二区三区久久| 宅男免费午夜| av片东京热男人的天堂| 建设人人有责人人尽责人人享有的| 亚洲性夜色夜夜综合| 成人18禁高潮啪啪吃奶动态图| 免费黄频网站在线观看国产| 一本色道久久久久久精品综合| 老汉色∧v一级毛片| 肉色欧美久久久久久久蜜桃| 人人澡人人妻人| 国产又爽黄色视频| 国产精品美女特级片免费视频播放器 | 乱人伦中国视频| 一级片免费观看大全| 久久 成人 亚洲| 人妻 亚洲 视频| 国产成人av教育| 老司机在亚洲福利影院| 国产高清激情床上av| 操出白浆在线播放| 国产精品免费视频内射| 国产av精品麻豆| 日本av手机在线免费观看| 国产xxxxx性猛交| 99热国产这里只有精品6| 9热在线视频观看99| 久久中文字幕人妻熟女| 香蕉久久夜色| 亚洲中文av在线| 夫妻午夜视频| 国产男靠女视频免费网站| 两性夫妻黄色片| 黄色成人免费大全| 91av网站免费观看| 一级毛片精品| kizo精华| 中文字幕色久视频| 一区二区日韩欧美中文字幕| 日韩大片免费观看网站| 激情在线观看视频在线高清 | 亚洲男人天堂网一区| 亚洲第一欧美日韩一区二区三区 | 在线天堂中文资源库| 久久久久久久久免费视频了| 国产91精品成人一区二区三区 | 精品亚洲乱码少妇综合久久| 青草久久国产| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 精品少妇内射三级| 一二三四社区在线视频社区8| 男女高潮啪啪啪动态图| 久久久精品免费免费高清| 老司机午夜十八禁免费视频| av电影中文网址| 日本黄色日本黄色录像| 久久ye,这里只有精品| 巨乳人妻的诱惑在线观看| 三级毛片av免费| 岛国在线观看网站| 日本精品一区二区三区蜜桃| 精品午夜福利视频在线观看一区 | 久久久国产欧美日韩av| 国产欧美日韩综合在线一区二区| 亚洲三区欧美一区| 这个男人来自地球电影免费观看| 欧美国产精品va在线观看不卡| 黄色视频,在线免费观看| 国精品久久久久久国模美| 日韩欧美三级三区| svipshipincom国产片| 久久久精品免费免费高清| 青青草视频在线视频观看| 亚洲视频免费观看视频| 午夜久久久在线观看| 日韩视频在线欧美| 18禁裸乳无遮挡动漫免费视频| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| 丝袜在线中文字幕| 黄片大片在线免费观看| 久热爱精品视频在线9| 国产精品偷伦视频观看了| 国产亚洲一区二区精品| 久久久久国内视频| 国产精品久久久久久精品电影小说| 一夜夜www| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产欧美日韩av| 色婷婷av一区二区三区视频| 精品亚洲成a人片在线观看| 性少妇av在线| 国产日韩欧美亚洲二区| 日本撒尿小便嘘嘘汇集6| 女人高潮潮喷娇喘18禁视频| 无人区码免费观看不卡 | 欧美成狂野欧美在线观看| 丝瓜视频免费看黄片| 日本一区二区免费在线视频| 久久人人97超碰香蕉20202| 18禁黄网站禁片午夜丰满| 亚洲国产成人一精品久久久| 免费观看a级毛片全部| 久久久久久免费高清国产稀缺| 一级,二级,三级黄色视频| 在线十欧美十亚洲十日本专区| 欧美精品一区二区大全| 一本久久精品| 亚洲一码二码三码区别大吗| 成人永久免费在线观看视频 |