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

    溫排水自動提取方法研究

    2023-11-29 01:34:16朱秀芳
    中國環(huán)境科學(xué) 2023年11期
    關(guān)鍵詞:核電站紅外水體

    朱秀芳,李 原

    溫排水自動提取方法研究

    朱秀芳1,2*,李 原3,4

    (1.北京師范大學(xué),遙感科學(xué)國家重點實驗室,北京 100875;2.北京師范大學(xué)地理科學(xué)學(xué)部,遙感科學(xué)與工程研究院,北京 100875;3.內(nèi)蒙古自治區(qū)河流與湖泊生態(tài)重點實驗室,內(nèi)蒙古 呼和浩特 010021;4.內(nèi)蒙古大學(xué)生態(tài)與環(huán)境學(xué)院,內(nèi)蒙古 呼和浩特 010021)

    為解決現(xiàn)有的遙感溫排水檢測的算法自動化程度低、時效性差的問題,提出了一種基于單時相熱紅外波段的溫排水自動化提取方法(IForest-SVM),并以福建省福清核電站為例,使用2018年10月至2019年1月內(nèi)的4期Landsat-8遙感影像對所提方法進行了測試與驗證.該方法首先利用孤立森林進行溫排水和正常水體樣本的自動提取,接著通過溫排水和排水口的空間鄰接關(guān)系純化溫排水樣本,進而利用支持向量機監(jiān)督分類提取溫排水像元,最后通過溫排水和排水口的空間鄰接關(guān)系剔除誤判像元得到最終的溫排水空間分布.測試結(jié)果顯示本研究所提出的方法在2018年11月15日、2018年12月1日和2019年1月18日三期圖像上檢測得到的溫排水的用戶精度和生產(chǎn)者精度的平均值分別為89.69%、94.97%和90.04%.在2018年10月30日無溫排水的影像上沒有出現(xiàn)誤檢,有效的避免了“虛警”的出現(xiàn).IForest-SVM只需要輸入遙感影像的熱紅外波段,無需額外代入其他參數(shù),具有可移植性好、普適性強、自動化程度高的優(yōu)點,在溫排水的實時監(jiān)測和快速發(fā)現(xiàn)業(yè)務(wù)中有很好的應(yīng)用前景.

    孤立森林;異常檢測;支持向量機;鄰接關(guān)系;溫排水

    隨著能源需求的增加和核能技術(shù)的發(fā)展,中國在沿海地區(qū)興建了許多的核電站.然而,核電站的發(fā)電效率要低于火電站,僅有30%~35%的核能轉(zhuǎn)化為了電能[1],剩余的大部分能量以熱能的形式隨冷卻水排出,使得周圍水體溫度迅速升高,從而影響了水生動植物的生長和繁殖[2].藻類對溫度變化尤其敏感,海水溫度的升高可能會導(dǎo)致藻類的異常爆發(fā).在風(fēng)和潮汐的驅(qū)使下,藻類會侵入核電站的水冷系統(tǒng),影響冷卻水的排放,阻礙核電站的安全運行[3-4].因此,研究核電附近水域溫排水的分布具有重要意義[5].

    溫排水監(jiān)測的方法主要包括數(shù)值模擬法、大面積實地測量法和遙感監(jiān)測法.數(shù)值模擬法可以獲得大范圍、高時空分辨率的二、三維溫度場分布[6],但其模擬結(jié)果的精度依賴于研究者所使用的模型[7].實地測量法能夠獲取最真實的數(shù)據(jù),但其只能提供關(guān)于溫度的分散的點數(shù)據(jù),不能提供完整的空間分布[8],且耗費的人力物力較多,部分研究中將其用于遙感算法的定標(biāo)與驗證工作[9].遙感技術(shù)憑借同步性好、觀測范圍大、可重復(fù)觀測和成本相對較低等優(yōu)勢成為溫排水監(jiān)測的重要方法[10-14].目前,基于遙感手段進行溫排水監(jiān)測也廣泛應(yīng)用在我國的諸多核電站附近的水域中,如大亞灣核電站、田灣核電站等[15-17].

    國際上早期的遙感溫排水監(jiān)測使用的是AVHRR數(shù)據(jù),后期隨著國內(nèi)核電站的興建也有相關(guān)的應(yīng)用[18-19],但由于其空間分辨率較低,無法準(zhǔn)確獲取溫排水的細節(jié)分布.21世紀以后隨著Landsat、MODIS、環(huán)境一號、高分五號等系列衛(wèi)星數(shù)據(jù)的發(fā)布,這個問題得到了一定程度上的解決,國內(nèi)外涌現(xiàn)出了一批相關(guān)的應(yīng)用成果[20-29].之后,也有部分學(xué)者通過X波段雷達數(shù)據(jù)進行溫排水的檢測[30],但滿足具有高時空分辨率的數(shù)據(jù)通常是地基雷達,其造價昂貴難以進行大規(guī)模應(yīng)用.載人航空遙感靈活性及空間分辨率都優(yōu)于衛(wèi)星遙感,但其對起降場地有特殊要求,備航手續(xù)繁瑣且成本昂貴,難以業(yè)務(wù)推廣[31].近些年隨著無人機技術(shù)的發(fā)展,利用無人機攜帶熱紅外傳感器進行核電廠溫排水監(jiān)測取得了部分成果[32],但其硬件設(shè)備與圖像處理技術(shù)仍在發(fā)展過程中,其飛行范圍以及飛行條件也受到一定的制約.總體來看,基于衛(wèi)星遙感數(shù)據(jù)的溫排水提取仍然是實際應(yīng)用中使用最廣泛的.

    目前基于遙感手段開展的溫排水監(jiān)測主要是使用含有熱紅外波段的遙感數(shù)據(jù)進行水體表面溫度的反演,再基于統(tǒng)計學(xué)的方法進行溫排水的判別,整個過程費時費力,很難做到快速實時的溫排水監(jiān)測.在溫度反演中常用的算法包括單通道法[33]、分裂窗法[34]和多通道法[35].單通道算法基于地表溫度輻射傳導(dǎo)方程推導(dǎo)而來,需要大氣透過率、地表比輻射率和大氣平均作用溫度這3個參數(shù),這些參數(shù)需要同步實測獲取;分裂窗法與多通道法雖然需要的實測數(shù)據(jù)較少但其模型復(fù)雜程度要遠高于單通道法,且對熱紅外波段數(shù)量有一定的要求.

    孤立森林是一種適用于異常檢測的機器學(xué)習(xí)算法[36].溫排水是排入自然水體中的溫度高于自然水體的水.受溫排水影響的水體像元可以看作是水體像元中的異常像元,其在熱紅外波段的值明顯區(qū)別于正常的像元,且溫排水的像元與核電站排放口在空間位置上存在鄰接關(guān)系.為此,本文提出了一種綜合使用孤立森林、空間鄰接關(guān)系和支持向量機的溫排水自動化監(jiān)測方法,助力溫排水實時監(jiān)測業(yè)務(wù)的開展,其主要創(chuàng)新在于無需進行水體表面溫度的反演,能夠?qū)崿F(xiàn)快速、自動化的溫排水監(jiān)測.

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

    1.1 研究區(qū)概況

    圖1 研究區(qū)概況

    a、b分別為2018年12月1日的Landsat-8真彩色合成圖和熱紅外波段圖像

    本研究選擇福建省福清市三山鎮(zhèn)南段的福清核電站作為研究目標(biāo),將其周圍的海域設(shè)置為研究區(qū)(圖1).福清核電站是中國第一座采用三代核電技術(shù)的核電站,占地面積約18km2,共有6個核電機組,總裝機容量為6600MW,是東南沿海的重要能源基地之一,也是中國核電行業(yè)的重要標(biāo)志性工程[37].福清核電站三面環(huán)海,東北部位置與陸地連接,周邊的道路發(fā)達、交通便利.其氣候為典型亞熱帶海洋性季風(fēng)氣候,區(qū)域內(nèi)海潮以半日為周期,每潮汐日內(nèi)高低潮各有兩次,漲潮時其流向偏西北,落潮時其流向偏東南,海域周圍多年平均水溫在22.3℃左右,夏季和冬季的海水溫度差異在10-15℃左右.前人研究發(fā)現(xiàn)該區(qū)域受溫排水影響的范圍大概在6km以內(nèi),受溫排水影響水溫升高一般不超過4℃,不同季節(jié)和不同潮汐條件下,溫升情況有所差異,冬季溫升分布范圍小于同潮態(tài)溫升分布范圍,漲潮時刻溫升范圍小于落潮時刻溫深分布范圍[37].

    1.2 研究數(shù)據(jù)及其預(yù)處理

    通過Google Earth Engine獲取經(jīng)過大氣校正與幾何校正后的Landsat-8衛(wèi)星數(shù)據(jù),獲取時間如表1所示,Landsat-8數(shù)據(jù)的波段信息如表2所示.

    表1 所用數(shù)據(jù)列表

    表2 Landsat-8各波段信息

    圖2 實驗數(shù)據(jù)

    首先基于2019年1月18日的Landsat-8數(shù)據(jù)對常年不變的海岸邊界進行矢量化,接著對4幅圖像分別計算自動水體提取指數(shù)(AWEI, Automated Water Extraction Index),具體公式如下所示:

    式中:2為綠波段,4為近紅外波段,5與7分別為兩個短波紅外波段.

    設(shè)置AWEI大于0的像元為水體像元[38],提取海岸邊界內(nèi)的水體像元得到4幅圖像上的監(jiān)測水體范圍,以該范圍為掩膜對4個時相的Landsat-8熱紅外波段進行裁剪,獲得最終用于監(jiān)測的影像,圖2是以波段10為例進行的展示.觀察比較四個時相的數(shù)據(jù),目視解譯和手工數(shù)字化得到各圖像上的溫排水范圍,作為算法驗證的真值.

    2 研究方法與技術(shù)路線

    2.1 孤立森林算法原理介紹

    孤立森林(IForest)是在第八屆IEEE數(shù)據(jù)挖掘國際會議上首次被提出的用于離群點檢測的算法[39].在孤立森林中離群點是那些分布稀疏且離密度高的群體較遠的樣本.孤立森林是通過構(gòu)建個孤立的二叉樹得到的,其基本流程是:1)從樣本空間中,隨機選取個樣本,構(gòu)成子空間,并將這些樣本放入樹的根節(jié)點;2)隨機選取一個樣本空間中的某一個特征,并選取該特征數(shù)據(jù)集中介于最大值和最小值中的某個數(shù)值,以為節(jié)點將特征向量中小于的數(shù)據(jù)放在當(dāng)前節(jié)點的左側(cè),大于的數(shù)據(jù)放在當(dāng)前節(jié)點的右側(cè),由此形成左右兩個新的特征子空間;3)分別在左右兩個新的特征子空間中,不斷重復(fù)上一個步驟,構(gòu)造出新的子節(jié)點,直到數(shù)據(jù)本身只包含一個樣本或者樹的高度已經(jīng)達到了預(yù)設(shè)的高度時,孤立樹建立完成.循環(huán)建立新的孤立樹,當(dāng)構(gòu)建完棵孤立二叉樹后,就形成了一片孤立森林.計算待測數(shù)據(jù)樣本到達每顆孤立樹的葉子節(jié)點的平均路徑長度,進而計算異常程度值,具體計算見公式(2)和(3)[40],路徑長度越短,對應(yīng)的異常程度值越高,越先被孤立出來.

    給定一個包含個樣本的數(shù)據(jù)集,樹的平均路徑長度為:

    式中:()為調(diào)和數(shù),該值可以被估計為ln()+ 0.5772156649.()為給定樣本數(shù)時路徑長度的平均值,用于標(biāo)準(zhǔn)化樣本的路徑長度().

    樣本的異常評分(,)定義為:

    式中:(())為樣本在一批孤立樹中找到自己所處的位置所花費的平均路徑長度.

    基于Python實現(xiàn)孤立森林的過程中,在計算異常評分時,根據(jù)數(shù)據(jù)點的局部密度來估算這個點的異常程度.密度越低的點會得到更高的異常評分.因為負數(shù)更能直觀的表示低密度,所以scikit-learn庫在原有公式的基礎(chǔ)上進行了調(diào)整,調(diào)整后的計算方法如公式4所示,調(diào)整后值越小對應(yīng)異常程度越高.

    2.2 溫排水自動監(jiān)測方法的基本處理流程

    受溫排水影響的水體的溫度會明顯高于同時刻周邊沒有受溫排水影響的水體,且溫排水與核電站排放口在空間位置上存在鄰接關(guān)系.為此,本文提出了一種綜合使用孤立森林、空間鄰接關(guān)系和支持向量機的溫排水自動化監(jiān)測方法,簡稱為IForest- SVM,其基本流程包括四個步驟:1)使用孤立森林計算監(jiān)測水域中各水體像元的熱紅外波段的異常評分;2)提取異常得分小于均值減2倍標(biāo)準(zhǔn)差的像元作為異常水體樣本,保留其中與核電站排水口鄰接的樣本圖斑作為溫排水樣本;3)提取異常評分為正值的像元,取中位數(shù)至最大值區(qū)間內(nèi)的像元作為正常的水體樣本;4)將兩類樣本代入支持向量機進行監(jiān)督分類,得到初步的溫排水圖斑,保留其中與核電站排水口鄰接的圖斑從而得到最終的溫排水空間分布.

    2.3 技術(shù)路線與實驗設(shè)計

    為了測試本文所提出的方法在實時監(jiān)測溫排水中的有效性,設(shè)計了如圖3所示的技術(shù)路線.具體步驟如下:1)通過目視解譯標(biāo)定核電站及其進出水口的位置;2)利用本文所提出的方法從不同時相的(2018年11月15日、2018年12月1日和2019年1月18日)Landsat-8的熱紅外波段衛(wèi)星數(shù)據(jù)中提取溫排水的空間分布.3)通過目視解譯,確定每期影像上溫排水的區(qū)域,以此為真值樣本,利用公式5和6計算溫排水提取的用戶精度與生產(chǎn)者精度.

    式中:UA為用戶精度;PA為生產(chǎn)者精度,R為正確分類的像元個數(shù);W為正常水體錯分為溫排水的像元個數(shù);A為實際是溫排水而錯分為正常水體的像元個數(shù).

    為了探討鄰接關(guān)系判斷對所提出方法(簡稱方法Ⅰ)的作用,本研究還進行了三組對照實驗:1)在利用孤立森林提取出異常水體像元的基礎(chǔ)上,提取與核電站排水口具有鄰接關(guān)系的像元作為溫排水的訓(xùn)練樣本,將樣本帶入到支持向量機的監(jiān)督分類器中,得到溫排水提取結(jié)果,即孤立森林+鄰接關(guān)系判斷+SVM,簡稱方法Ⅱ;2)將經(jīng)過孤立森林提取的樣本直接帶入支持向量機的監(jiān)督分類器中,得到溫排水提取結(jié)果,即孤立森林+無鄰接關(guān)系判斷+SVM,簡稱方法Ⅲ;3)在對比實驗2的基礎(chǔ)上,在監(jiān)督分類完成之后,提取與核電站排水口具有鄰接關(guān)系的圖斑作為最終結(jié)果,即孤立森林+無鄰接關(guān)系判斷+ SVM+鄰接關(guān)系判斷,簡稱方法Ⅳ.

    圖3 技術(shù)路線

    3 研究結(jié)果

    不同方法識別的溫排水的空間分布如圖4所示.根據(jù)目視解譯的結(jié)果,使用ArcGIS計算溫排水區(qū)域的面積得到2018年11月15日、12月1日和2019年1月18日溫排水總面積分別為0.76 km2、8.30km2和2.25km2.對不同方法監(jiān)測出的溫排水進行驗證,結(jié)果如表3所示.綜合圖4和表2可以得出:

    無論哪種方法三個時期提取的溫排水的生產(chǎn)者精度均明顯高于用戶精度.用戶精度是正確分類的溫排水像元個數(shù)與分類得到的所有溫排水像元個數(shù)的比例,用戶精度越高說明錯分誤差越低.生產(chǎn)者精度是正確分類的溫排水像元個數(shù)與實際溫排水像元個數(shù)的比例,生產(chǎn)者精度越高說明漏分誤差越低.因此,總體來說上述幾種方法的漏分誤差均低于錯分誤差.

    在四種方法中,本研究所提出的方法在三個時期檢測出的溫排水的精度都是最高的,2018年11月15日、2018年12月1日和2019年1月18日三期圖像上檢測得到的溫排水的用戶精度和生產(chǎn)者精度的平均值分別為89.69%、94.97%和90.04%.

    圖4 不同方法提取的溫排水

    Ⅰ、Ⅱ、Ⅲ、Ⅳ分別指孤立森林+鄰接關(guān)系判斷+SVM+鄰接關(guān)系判斷的方法(即本研究所提方法)、孤立森林+鄰接關(guān)系判斷+SVM的方法、孤立森林+無鄰接關(guān)系判斷+SVM的方法、孤立森林+無鄰接關(guān)系判斷+SVM+鄰接關(guān)系判斷的方法

    表3 精度評價表(%)

    注:Ⅰ、Ⅱ、Ⅲ、Ⅳ的含義同圖4.

    對比添加和未添加鄰接關(guān)系的方法,可以看出總體來說添加了兩次鄰接關(guān)系后可以大大提高用戶精度、減少錯分誤差.例如2018年11月15日,加入兩次鄰接關(guān)系的方法Ⅰ用戶精度達到了99.82%,比僅在孤立森林后添加鄰接關(guān)系的方法Ⅱ提升了44.85%.兩次鄰接關(guān)系的使用中,第一次鄰接關(guān)系的重要性要高于第二次.在制備樣本時考慮空間鄰接關(guān)系的方法(方法I和方法Ⅱ)要好于不考慮的方法,且在僅有少量溫排水爆發(fā)的情況下(2018年11月15日和2019年1月18日)表現(xiàn)的尤為明顯.錯誤的樣本選擇甚至可以導(dǎo)致無法檢測出溫排水,例如2018年11月15日和2019年1月18方法Ⅳ得到結(jié)果.

    4 討論

    4.1 本文所提方法的優(yōu)勢

    傳統(tǒng)的溫排水監(jiān)測方法需要在進行輻射定標(biāo)、大氣校正等一系列預(yù)處理的前提下,進行水體表面溫度的反演,再進行統(tǒng)計學(xué)上的分析,篩選出異常位置.整套流程涉及方法諸多,過程較為繁瑣,尤其是對于溫度的反演可能還需要輸入遙感數(shù)據(jù)之外的其他因子(如地表比輻射率、大氣在熱紅外的透過率、大氣向下輻射亮度、大氣向上輻射亮度等).本研究提出的IForest-SVM方法無需進行溫度反演,相較于傳統(tǒng)方法更能做到實時監(jiān)測.

    IForest-SVM通過孤立森林計算監(jiān)測區(qū)水體像元的異常評分,結(jié)合均值標(biāo)準(zhǔn)差法確定閾值來提取溫排水樣本,不同區(qū)域或者同一區(qū)域不同時間計算得到的溫排水閾值均不同,也即所提取的閾值是動態(tài)和自適應(yīng)的,保證了方法的可移植性.此外,自動化的訓(xùn)練樣本提取大大減少了樣本采集上的人力、物力和時間投入.

    要想實現(xiàn)溫排水的實時監(jiān)測,需要對監(jiān)測的水域進行持續(xù)跟蹤,所用算法在無溫排水發(fā)生時不能產(chǎn)生虛警,在有溫排水發(fā)生時要盡量減少錯報.以往研究往往利用溫排水發(fā)生時的影像測試算法的有效性,這類算法往往是“馬后炮”式的,不適用于實時監(jiān)測.為了進一步驗證IForest-SVM方法在溫排水實時監(jiān)測中的有效性,本文不僅對有溫排水發(fā)生的圖像進行了溫排水的提取實驗(圖4),還對沒有溫排水發(fā)生的2018年10月30日的Landsat-8的熱紅外波進行了相同的處理,觀察是否有“虛警”出現(xiàn),結(jié)果如圖5所示,在未爆發(fā)溫排水時IForest-SVM不會產(chǎn)出虛警.結(jié)合圖5可知IForest-SVM在沒有溫排水發(fā)生、有少量溫排水發(fā)生和大量溫排水發(fā)生時都是適用的.

    圖5 未爆發(fā)溫排水時的檢測結(jié)果

    4.2 影響本文所提方法的精度和效率的因素

    首先,在基于孤立森林制備樣本時,標(biāo)準(zhǔn)差倍數(shù)的選取對樣本的制備具有一定的影響.圖6展示了溫排水樣本提取結(jié)果的生產(chǎn)者精度和用戶精度隨著標(biāo)準(zhǔn)差倍數(shù)增加(設(shè)定標(biāo)準(zhǔn)差倍數(shù)在1到2之間以0.1為步長變化)而變化的情況.總體來說,隨著標(biāo)準(zhǔn)差倍數(shù)的增加生產(chǎn)者精度變化并不大,最低也在0.98以上,而用戶精度則表現(xiàn)出了明顯的提升.這說明孤立森林加鄰接關(guān)系所得到的樣本基本能夠包含所有的核電站溫排水像元,漏選的樣本很少;但當(dāng)標(biāo)準(zhǔn)差倍數(shù)過小時錯提進來的溫排水樣本多,隨著標(biāo)準(zhǔn)差倍數(shù)的提高,用戶精度提高,錯提的溫排水樣本變少.此外,在不同時相的數(shù)據(jù)中用戶精度隨標(biāo)準(zhǔn)差倍數(shù)增加而提高的速率不一樣,但整體來說接近2倍標(biāo)準(zhǔn)差時三個時相下用戶精度提升的速率均明顯減小了.

    其次,監(jiān)督分類算法的選擇會對最終結(jié)果造成影響.本文選擇了常用的SVM方法與自動化樣本提取進行組合,所提方法的核心是自動化的樣本提取和鄰接關(guān)系的判定,SVM監(jiān)督分類器也可以替換為其他的分類器.不同分類器對樣本的要求不同,相同樣本得到的分類精度會存在差異.未來可以嘗試和其他監(jiān)督分類器進行組合,優(yōu)選精度高且穩(wěn)定的分類器.

    再者,本研究所用的Landsat數(shù)據(jù)重訪周期較長,不利于實時監(jiān)測,未來可以考慮架設(shè)熱紅外攝像頭或者使用更高時空分辨率的數(shù)據(jù)(如哨兵)進行監(jiān)測.對于時間分辨率高而空間分辨率偏低的數(shù)據(jù)(如MODIS)也可以結(jié)合降尺度算法,在提高其空間分辨率的基礎(chǔ)上結(jié)合本文所提出的算法進行溫排水的監(jiān)測.

    圖6 基于不同閾值選取的樣本進行溫排水提取的精度

    最后,熱帶海域、溫帶海域、寒帶海域的溫排水與周圍海水的溫差是不同的,反應(yīng)在熱紅外波段上的溫排水和正常水體的像元值差異大小不同.本文的算法僅僅在福清核電站進行了測試,福清核電站位于亞熱帶海域,本文算法的適用性和精度在其他溫差條件不一樣的海域有待于進一步驗證.

    5 結(jié)論

    5.1 本文提出了一種基于單時相熱紅外波段的溫排水自動化提取方法,并以福清核電站不同時期的溫排水檢測為例進行了算法的精度測試.測試結(jié)果顯示2018年11月15、12月1日和2019年1月18日的溫排水檢測的生產(chǎn)者精度分別為79.56%、95.17%和98.46%,用戶精度分別為99.82%、94.77%和81.62%.

    5.2 加入溫排水像元和排水口鄰接關(guān)系的判定,大大提高了溫排水識別的精度.在樣本提取階段加入鄰接關(guān)系的判斷尤為重要,樣本提取階段不加入鄰接關(guān)系的分類方法(方法Ⅲ和Ⅳ)的精度甚至低至0%.

    5.3 在無溫排水的時相(2018年10月30日),本文所提方法未造成誤檢.

    [1] Liu M, Yin X, Xu Q, et al. Monitoring of fine-scale warm drain-off water from nuclear power stations in the daya bay based on landsat 8data [J]. Remote Sensing, 2020,12(4):627-648.

    [2] Verones F, Hanafiah M M., Pfister S, et al. Characterization factors for thermal pollution in freshwater aquatic environments [J]. Environmental Science & Technology, 2010,44(24):9364-9369.

    [3] 於 凡,張永興.濱海核電站溫排水對海洋生態(tài)系統(tǒng)影響的研究 [J]. 輻射防護通訊, 2008,(1):1-7. Yu F, ZhangY X. The review onthe effects of thermal effluent from nuclear plants on the marine ecosystem [J]. Radiation Protection Newsletter, 2008,(1):1-7.

    [4] 賀立燕,宋秀賢,於 凡,等.潛在影響防城港核電冷源系統(tǒng)的藻類暴發(fā)特點及其監(jiān)測防控技術(shù) [J]. 海洋與湖沼, 2019,50(3):700-706. He L Y, Song X X, Yu F, et al. potential risk and prevention of phytoplankton outbreak to water-cooling system in nuclear power plant in Fangchenggang, Guangxi [J]. Oceanologia Et Limnologia Sinica, 2019,50(3):700-706.

    [5] 劉永葉,陳曉秋.核電廠溫排水熱影響研究的建議[J]. 輻射防護通訊, 2011,31(6):20-23. Liu Y Y, Chen X Q. Suggestion of the study on thermal impact of thermal discharge from NPPs [J]. Radiation Protection Newsletter, 2011,31(6):20-23.

    [6] 張 琨,張愛玲,覃春麗,等.CORMIX與二維數(shù)模、物模在濱海核電廠溫排水近區(qū)模擬中的比對研究 [J]. 海洋科學(xué)進展, 2017,35(1): 117-123. Zhang K, Zhang A L, Qin C L, et al. Comparison between CORMIX and traditional models on the simulation of coastal nuclear power plant water discharge [J]. Advances in Marine Science, 2017,35(1):117-123.

    [7] 趙 強,曹 維,邢 健,等.基于現(xiàn)場觀測數(shù)據(jù)的三門核電站冬季溫排水時空特征分析 [J]. 海洋環(huán)境科學(xué), 2022,41(6):847-856. Zhao Q, Cao W, Xing J, et al. Spatio-temporal characteristics of thermal discharge from Sanmen Nuclear Power Plant in winter based on field observations [J]. Marine Environmental Science, 2022,41(6): 847-856.

    [8] Muthulakshmi A, Natesan U, Ferrer V A,et al. A novel technique to monitor thermal discharges using thermal infrared imaging [J]. Environmental Science: Processes & Impacts, 2013,15(9):1729-1734.

    [9] 何 園.遙感技術(shù)在核電站溫排水影響范圍監(jiān)測中的應(yīng)用 [J]. 科技創(chuàng)新與應(yīng)用, 2016,(18):71. He Y. Application of remote sensing technology in monitoring the scope of influence of warm drainage in nuclear power plants [J]. Technological Innovation and Application, 2016,(18):71.

    [10] 姜 晟,李俊龍,李旭文,等.核電站溫排水遙感監(jiān)測方法研究——以田灣核電站為例 [J]. 中國環(huán)境監(jiān)測, 2013,29(6):212-216. Jiang S, Li J L, Li X W, et al. A research on the remote sensing monitoring method with the thermal discharge of Tian-wan nuclear power station [J]. Environmental Monitoring in China, 2013,29(6): 212-216.

    [11] 姚沛林.濱海核電廠溫排水的監(jiān)測與實踐 [J]. 紅外, 2013,34(11): 43-48. Yao P L. Monitoring of warm water discharged from coastal nuclear power plant [J]. Infraied, 2013,34(11):43-48.

    [12] 賀佳惠,梁春利,李名松.核電站近岸溫度場航空熱紅外遙感測量數(shù)據(jù)處理研究 [J]. 國土資源遙感, 2010,(3):51-53. He J H, Liang C L, Li M S. Temperature field airborne thermal remote sensing surveyof the alongshore nuclear power station [J]. Remote Sensing for Land & Resources, 2010,(3):51-53.

    [13] 張曉峰,徐京萍,張曼祺,等.Landsat-8熱紅外數(shù)據(jù)監(jiān)測田灣核電站溫排水分布 [J]. 海洋科學(xué)進展, 2019,37(3):518-525. Zhang X F, Xu J P, Zhang M Q, et al. Distribution of thermal discharge from Tianwan nuclear power plant: analysis of thermal infrared data from landsat-8 [J]. Advances in Marine Science, 2019,37(3):518-525.

    [14] 王 祥,蘇 岫,張 浩,等.不同空間分辨率遙感數(shù)據(jù)在核電站溫排水監(jiān)測中的應(yīng)用研究 [J]. 海洋環(huán)境科學(xué), 2020,39(4):646-651. Wang X, Su X, Zhang H, et al. Remote sensing based application research of nuclear power plant thermal plume monitoring with different spatial resolution imagery [J]. Marine Environmental Science, 2020,39(4):646-651.

    [15] 張 彩,朱 利,賈 祥,等.不同空間分辨率熱紅外數(shù)據(jù)在近海核電廠溫排水監(jiān)測一致性研究 [J]. 遙感信息, 2015,30(2):71-76. Zhang C, Zhu L, Jia X, et al. Consistency of different spatial resolution thermal infrared data in monitoring thermal discharge of nuclear power plant in coastal area [J]. Remote Sensing Information, 2015, 30(2):71-76.

    [16] 梁珊珊,張 兵,李俊生,等.環(huán)境一號衛(wèi)星熱紅外數(shù)據(jù)監(jiān)測核電站溫排水分布——以大亞灣為例 [J]. 遙感信息, 2012,(2):41-46. Liang S S, Zhang B, Li J S,et al. Distribution of therm-water pollution of nuclear power plant using the thermal infrared band of HJ-IRS data-taking Daya Bay as an example [J]. Remote Sensing Information, 2012, (2):41-46.

    [17] 周 穎,鞏彩蘭,匡定波,等.基于環(huán)境減災(zāi)衛(wèi)星熱紅外波段數(shù)據(jù)研究核電廠溫排水分布 [J]. 紅外與毫米波學(xué)報, 2012,31(6):544-549. Zhou Y, Gong C L, Kuang D B, et al. Distribution of thermal discharge from a power plant:Analysis of thermal infrared data from the environmental mitigation satellite [J]. J. Infrared Millim. Waves, 2012,31(6):544-549.

    [18] Tang D, Kester D R, Wang Z,et al. AVHRR satellite remote sensing and shipboard measurements of the thermal plume from the Daya Bay, nuclear power station, China [J]. Remote Sensing of Environment, 2003,84(4):506-515.

    [19] McClain E P, Pichel W G, Walton C C. Comparative performance of AVHRR-based multichannel sea surface temperatures [J]. Journal of Geophysical Research: Oceans, 1985,90(C6):11587-11601.

    [20] Chen C, Shi P, Mao Q. Application of remote sensing techniques for monitoring the thermal pollution of cooling-water discharge from nuclear power plant [J]. Journal of Environmental Science and Health, Part A, 2003,38(8):1659-1668.

    [21] Mustard J, Carney M, Sen A. The use of satellite data to quantify thermal effluent impacts [J]. Estuarine, Coastal and Shelf Science, 1999,49(4):509-524.

    [22] Gibbons D, Wukelic G, Leighton J, et al. Application of landsat thematic mapper data for coastal thermal plume analysis at Diablo Canyon [J]. Photogrammetric Engineering and Remote Sensing;(USA), 1989,55:

    [23] 朱 利,趙利民,王 橋,等.核電站溫排水分布衛(wèi)星遙感監(jiān)測及驗證 [J]. 光譜學(xué)與光譜分析, 2014,34(11):3079-3084. Zhu L, Zhao L M,Wang Q, et al. Monitoring the thermal plume from coastal nuclear power plant using satellite remote sensing data: modeling and validation [J]. Spectroscopy and Spectral Analysis, 2014,34(11):3079-3084.

    [24] 許 靜,朱 利,姜 建,等.基于HJ-1B與TM熱紅外數(shù)據(jù)的大亞灣核電基地溫排水遙感監(jiān)測 [J]. 中國環(huán)境科學(xué), 2014,34(5):1181-1186.Xu J, Niu L, Jiang J, et al. Monitoring thermal discharge in Daya Bay plant based on thermal infrared band of HJ-1B and TM remote sensing data [J]. China Environmental Science, 2014,34(5):1181-1186.

    [25] Benoit J, Fox-Kemper B. Contextualizing thermal effluent impacts in narragansett bay using Landsat-derived surface temperature [J]. Frontiers in Marine Science, 2021,8:705204.

    [26] Liu Y, Huang S, Yi P,et al. The thermal discharge retrieve and analysis of Hongyanhe district with Landsat-8 and GF-5 Image [C]// Proceedings of the IOP Conference Series: Earth and Environmental Science, F, 2020.IOP Publishing.

    [27] Wu C, Wang Q, Yang Z,et al. Monitoring heated water pollution of the DaYaWan nuclear power plant using TM images [J]. International Journal of Remote Sensing, 2007,28(5):885-890.

    [28] Nie P, Zhu H, Xu H,et al. Monitoring of Tianwan nuclear power plant thermal pollution based on remotely sensed landsat data [C]// Proceedings of the IGARSS 2020-2020 IEEE International Geoscience and Remote Sensing Symposium, F, 2020. IEEE.

    [29] Ma P, Dai X, Guo Z,et al. Detection of thermal pollution from power plants on China's eastern coast using remote sensing data [J]. Stochastic Environmental Research and Risk Assessment, 2017,(31): 1957-1975.

    [30] Bau H Y, Wu L C, Chuang L Z H. Detection of thermal effluent discharge using sea surface nautical radar images [J]. IEEE Geoscience and Remote Sensing Letters, 2012,10(3):612-616.

    [31] 王 祥,王新新,蘇 岫,等.無人機平臺航空遙感監(jiān)測核電站溫排水——以遼寧省紅沿河核電站為例 [J]. 自然資源遙感, 2018,30(4): 182-186. Wang X, Wang X X, Su X, et al. Thermal discharge monitoring of nuclear power plant with aerial remote sensing technology using a UAV platform: Take Hongyanhe Nuclear Power Plant, Liaoning Province, as example [J]. Remote Sensing for Land & Resources, 2018,30(4):182-186.

    [32] 孫芹芹,張加晉,姬厚德,等.基于衛(wèi)星和無人機的后石電廠溫排水分布研究 [J]. 應(yīng)用海洋學(xué)學(xué)報, 2020,39(2):261-265. Sun Q Q, Zhang J J, Ji H D, et al. Distribution of thermal water around Houshi Power Plant based on landsat 8 and Uva study [J]. Journal of Applied Oceanography, 2020,39(2):261-265.

    [33] Jiménez-Mu?oz J C, Sobrino J A. A generalized single‐channel method for retrieving land surface temperature from remote sensing data [J]. Journal of geophysical research: Atmospheres, 2003,108(D22):

    [34] Wan Z, Dozier J. A generalized split-window algorithm for retrieving land-surface temperature from space [J]. IEEE Transactions on Geoscience and Remote Sensing, 1996,34(4):892-905.

    [35] Wan Z M, Li Z L. A physics-based algorithm for retrieving land- surface emissivity and temperature from EOS/MODIS data [J]. IEEE Transactions on Geoscience and Remote Sensing, 1997,35(4):980-996.

    [36] Liu F T, Ting K M, Zhou Z H. Isolation-based anomaly detection [J]. ACM Transactions on Knowledge Discovery from Data, 2012:6(1): 1-39.

    [37] 董雙發(fā),范 曉,石海崗,等.基于Landsat8和無人機的福清核電溫排水分布研究 [J]. 自然資源遙感, 2022,34(3):112-120. Dong S F, Fan X, Shi H G, et al. Study on distribution of thermal discharge in Fuqing nuclear power plant based on landsat8 and UAV [J]. Remote Sensing for Natural Resources, 2022,34(3):112-120.

    [38] Feyisa G L, Meilby H, Fensholt R,et al. Automated water extraction index: a new technique for surface water mapping using Landsat imagery [J]. Remote Sensing of Environment, 2014,140:23-35.

    [39] Liu F T, Ting K M, Zhou Z H. Isolation forest [C]//Proceedings of the 2008 Eighth Ieee International Conference on Data Mining, F, 2008. IEEE.

    [40] Liu F T, Ting K M, Zhou Z H. Isolation-based anomaly detection [J]. ACM Transactions on Knowledge Discovery from Data (TKDD), 2012,6(1):1-39.

    Research on automatic extraction method of thermal discharge.

    ZHU Xiu-fang1,2*, LI Yuan3,4

    (1.State Key Laboratory of Remote Sensing Science, Beijing Normal University, Beijing 100875, China;2.Institute of Remote Sensing Science and Engineering, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China;3.The Inner Mongolia Key Laboratory of River and Lake Ecology, Hohhot 010021, China;4.School of Ecology and Environment, Inner Mongolia University, Hohhot 010021, China)., 2023,43(11):6115~6122

    To solve the problems of low automation and poor timeliness of existing algorithms for remote sensing thermal discharge detection, this study proposes an automated extraction method for thermal discharge based on the single-temporal thermal infrared band, called IForest-SVM. The Fuqing Nuclear Power Plant in Fujian Province was taken as an example to test and verify the proposed method using four Landsat-8 remote sensing images from October 2018 to January 2019. IForest-SVM first uses isolation forest to automatically extract samples of thermal discharge and normal water bodies, then purifies the thermal discharge samples by their spatial adjacency relationship with discharge outlets. After that, a support vector machine is used to extract thermal discharge pixels, and misclassified pixels are removed based on their spatial adjacency relationship with discharge outlets to obtain the final spatial distribution of thermal discharge. The test results show that the average values of user accuracy and producer accuracy of thermal discharge detected by the proposed method on three images on November 15, 2018, December 1, 2018, and January 18, 2019 were 89.69%, 94.97%, and 90.04%, respectively. No false alarms occurred on the image without thermal discharge on October 30, 2018, effectively avoiding false detections. IForest-SVM only requires input of the thermal infrared band of remote sensing images without additional parameters, and has the advantages of good portability, strong applicability, and high automation. It has a good application potential in the real-time monitoring and rapid detection of thermal discharge.

    isolation forest;anomaly detection;support vector machine;adjacency relationship;thermal discharge

    X87

    A

    1000-6923(2023)11-6115-08

    朱秀芳(1982-),女,浙江天合人,教授,博士,研究方向為遙感應(yīng)用.發(fā)表論文160余篇.zhuxiufang@bnu.edu.cn.

    朱秀芳,李 原.溫排水自動提取方法研究 [J]. 中國環(huán)境科學(xué), 2023,43(11):6115-6122.

    Zhu X F, Li Y.Research on automatic extraction method of thermal discharge [J]. China Environmental Science, 2023,43(11):6115-6122.

    2023-03-21

    國家自然科學(xué)基金資助項目(41292583)

    * 責(zé)任作者, 教授, zhuxiufang@bnu.edu.cn

    猜你喜歡
    核電站紅外水體
    如果離我不遠的核電站出現(xiàn)泄漏該怎么辦
    飛碟探索(2022年2期)2022-05-31 01:28:03
    網(wǎng)紅外賣
    農(nóng)村黑臭水體治理和污水處理淺探
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    閃亮的中國紅外『芯』
    金橋(2021年4期)2021-05-21 08:19:20
    TS系列紅外傳感器在嵌入式控制系統(tǒng)中的應(yīng)用
    電子制作(2019年7期)2019-04-25 13:17:14
    核電站護衛(wèi)隊
    核電站的“神經(jīng)中樞”
    中國核電(2017年2期)2017-08-11 08:01:12
    廣元:治理黑臭水體 再還水清岸美
    海上核電站提速
    能源(2016年3期)2016-12-01 05:11:05
    麻豆一二三区av精品| 美女高潮喷水抽搐中文字幕| 国产国拍精品亚洲av在线观看| 亚洲 欧美 日韩 在线 免费| 亚洲国产高清在线一区二区三| 色5月婷婷丁香| 欧美高清成人免费视频www| 欧美潮喷喷水| 9191精品国产免费久久| 久久久久国产精品人妻aⅴ院| 在线十欧美十亚洲十日本专区| 在线观看66精品国产| 深夜a级毛片| 舔av片在线| 国产人妻一区二区三区在| 老司机午夜福利在线观看视频| 在线播放无遮挡| 999久久久精品免费观看国产| 丰满的人妻完整版| 高清毛片免费观看视频网站| 日本一本二区三区精品| 成人特级黄色片久久久久久久| 精品国产三级普通话版| 宅男免费午夜| 免费电影在线观看免费观看| 3wmmmm亚洲av在线观看| 久久国产乱子免费精品| 亚洲精品在线观看二区| www日本黄色视频网| 精品乱码久久久久久99久播| 麻豆国产97在线/欧美| 给我免费播放毛片高清在线观看| 国产成人a区在线观看| 精品福利观看| 男人舔女人下体高潮全视频| 久久天躁狠狠躁夜夜2o2o| 国产单亲对白刺激| 伊人久久精品亚洲午夜| 亚洲精品一区av在线观看| 国产三级中文精品| 久久热精品热| 91av网一区二区| 国产精品免费一区二区三区在线| 久久久久久久午夜电影| 日韩欧美免费精品| 麻豆av噜噜一区二区三区| 免费av观看视频| 日韩人妻高清精品专区| 99国产综合亚洲精品| 老司机午夜十八禁免费视频| 欧美乱色亚洲激情| 久久国产乱子免费精品| 99久久无色码亚洲精品果冻| 日韩 亚洲 欧美在线| www.熟女人妻精品国产| 黄色一级大片看看| 嫩草影视91久久| 综合色av麻豆| 免费看日本二区| 国产伦在线观看视频一区| 久久久久亚洲av毛片大全| 亚洲无线在线观看| 99久国产av精品| 在线观看舔阴道视频| 91狼人影院| a级毛片免费高清观看在线播放| 色噜噜av男人的天堂激情| 少妇被粗大猛烈的视频| 一二三四社区在线视频社区8| 国内毛片毛片毛片毛片毛片| 亚洲五月天丁香| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品456在线播放app | 窝窝影院91人妻| 欧美一区二区精品小视频在线| 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 一级黄色大片毛片| 97超视频在线观看视频| 午夜激情福利司机影院| 亚洲欧美日韩东京热| 丁香欧美五月| 嫁个100分男人电影在线观看| 日本一本二区三区精品| 91久久精品电影网| 日韩欧美免费精品| 天堂影院成人在线观看| 大型黄色视频在线免费观看| 淫秽高清视频在线观看| 日本黄色片子视频| 国产老妇女一区| 99热只有精品国产| 亚洲人成伊人成综合网2020| 国产伦精品一区二区三区四那| 啦啦啦韩国在线观看视频| 一级毛片久久久久久久久女| 欧美一区二区国产精品久久精品| 久久久久久国产a免费观看| 日日干狠狠操夜夜爽| 美女 人体艺术 gogo| 国产aⅴ精品一区二区三区波| 国产免费男女视频| .国产精品久久| 午夜激情福利司机影院| 一进一出好大好爽视频| av国产免费在线观看| 久久久久久大精品| 夜夜爽天天搞| 亚洲精华国产精华精| 久久久久九九精品影院| 国产精品伦人一区二区| 窝窝影院91人妻| 三级男女做爰猛烈吃奶摸视频| 亚洲专区中文字幕在线| 天美传媒精品一区二区| av在线老鸭窝| 人人妻人人澡欧美一区二区| 两个人的视频大全免费| 露出奶头的视频| 国产成人欧美在线观看| 亚洲 欧美 日韩 在线 免费| 国产成人福利小说| 欧美成人免费av一区二区三区| 国产白丝娇喘喷水9色精品| 91麻豆精品激情在线观看国产| 久久久久九九精品影院| 欧美日韩国产亚洲二区| 国产高清激情床上av| 嫁个100分男人电影在线观看| 欧美潮喷喷水| 一区福利在线观看| 午夜精品久久久久久毛片777| 午夜日韩欧美国产| 日本 av在线| 久久天躁狠狠躁夜夜2o2o| 淫妇啪啪啪对白视频| 欧美一区二区精品小视频在线| 国产精品精品国产色婷婷| 精品午夜福利视频在线观看一区| 国产精品一区二区性色av| 亚洲av电影在线进入| 美女高潮的动态| 激情在线观看视频在线高清| 国产精品免费一区二区三区在线| 国产中年淑女户外野战色| 午夜日韩欧美国产| 一进一出好大好爽视频| 国产成人aa在线观看| 99国产精品一区二区三区| www.色视频.com| 久久久色成人| 内射极品少妇av片p| 精品久久久久久久久亚洲 | 久9热在线精品视频| 在线国产一区二区在线| 他把我摸到了高潮在线观看| 热99re8久久精品国产| 99国产精品一区二区三区| 国产欧美日韩一区二区三| 在现免费观看毛片| 国产午夜精品论理片| 久久精品久久久久久噜噜老黄 | 又爽又黄a免费视频| 人妻夜夜爽99麻豆av| 美女xxoo啪啪120秒动态图 | a级毛片a级免费在线| 男人狂女人下面高潮的视频| 亚洲 欧美 日韩 在线 免费| 亚洲av第一区精品v没综合| 天堂影院成人在线观看| 久久国产乱子伦精品免费另类| 中亚洲国语对白在线视频| 人妻丰满熟妇av一区二区三区| 国产精品影院久久| 性色avwww在线观看| 波多野结衣巨乳人妻| 日韩中字成人| 国产精品美女特级片免费视频播放器| 精品福利观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲第一电影网av| 国产爱豆传媒在线观看| 久久午夜福利片| 看片在线看免费视频| 97人妻精品一区二区三区麻豆| 午夜老司机福利剧场| 丁香六月欧美| 久久久久久久久大av| 日韩有码中文字幕| 99视频精品全部免费 在线| 国产一区二区激情短视频| 91麻豆av在线| 亚洲av成人不卡在线观看播放网| 最后的刺客免费高清国语| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久精品电影| 简卡轻食公司| 一级黄色大片毛片| 嫩草影院新地址| 最近最新中文字幕大全电影3| av天堂在线播放| 天天躁日日操中文字幕| 国产亚洲精品av在线| 一进一出好大好爽视频| 99久国产av精品| 麻豆久久精品国产亚洲av| 最近中文字幕高清免费大全6 | 在线观看舔阴道视频| 欧美在线黄色| 国内毛片毛片毛片毛片毛片| 夜夜夜夜夜久久久久| 国产精品人妻久久久久久| 国产亚洲精品久久久久久毛片| 亚洲成人精品中文字幕电影| 亚洲七黄色美女视频| 女人被狂操c到高潮| 欧美精品啪啪一区二区三区| 少妇人妻一区二区三区视频| 免费看美女性在线毛片视频| 无人区码免费观看不卡| 又爽又黄a免费视频| 婷婷精品国产亚洲av| 99久久精品热视频| 欧美xxxx黑人xx丫x性爽| 直男gayav资源| 亚洲中文字幕日韩| 免费av毛片视频| 亚洲专区中文字幕在线| 亚洲欧美日韩无卡精品| 国产黄a三级三级三级人| 精品久久久久久久末码| 国产av一区在线观看免费| av福利片在线观看| 国产成年人精品一区二区| 国产亚洲精品久久久com| 亚洲欧美日韩高清在线视频| 黄色一级大片看看| 精品久久久久久久久久免费视频| 免费无遮挡裸体视频| 亚洲熟妇熟女久久| 首页视频小说图片口味搜索| 在线看三级毛片| 黄色一级大片看看| 在线观看一区二区三区| 一区二区三区免费毛片| 日日夜夜操网爽| 国产午夜精品久久久久久一区二区三区 | bbb黄色大片| 久久久久久久亚洲中文字幕 | 真人一进一出gif抽搐免费| 嫩草影视91久久| 一进一出好大好爽视频| 亚洲人成网站在线播| 欧美+亚洲+日韩+国产| 色精品久久人妻99蜜桃| 亚洲无线观看免费| 一个人免费在线观看的高清视频| 97超级碰碰碰精品色视频在线观看| 成人美女网站在线观看视频| 国产精品久久久久久久久免 | 日本熟妇午夜| av中文乱码字幕在线| 51午夜福利影视在线观看| 中文字幕免费在线视频6| 欧美不卡视频在线免费观看| 精品乱码久久久久久99久播| 精品一区二区三区视频在线| 在线观看免费视频日本深夜| 天天躁日日操中文字幕| 黄色配什么色好看| 我要搜黄色片| 久久久久免费精品人妻一区二区| 亚洲熟妇熟女久久| 高潮久久久久久久久久久不卡| 国产主播在线观看一区二区| 亚洲av不卡在线观看| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 免费无遮挡裸体视频| 国内毛片毛片毛片毛片毛片| 亚洲国产日韩欧美精品在线观看| АⅤ资源中文在线天堂| 亚洲av成人精品一区久久| 无遮挡黄片免费观看| 久久天躁狠狠躁夜夜2o2o| 中文字幕久久专区| 2021天堂中文幕一二区在线观| 国内毛片毛片毛片毛片毛片| 午夜福利成人在线免费观看| 精品不卡国产一区二区三区| 男女床上黄色一级片免费看| 熟女人妻精品中文字幕| 国产午夜精品论理片| 婷婷色综合大香蕉| 日韩免费av在线播放| 日韩欧美三级三区| 亚洲欧美日韩高清专用| 色综合婷婷激情| 日本 av在线| 99精品久久久久人妻精品| 国内精品久久久久精免费| 色在线成人网| 色播亚洲综合网| 亚洲欧美激情综合另类| 韩国av一区二区三区四区| 婷婷色综合大香蕉| 免费观看的影片在线观看| 人人妻人人看人人澡| 最新中文字幕久久久久| 免费观看的影片在线观看| 国产亚洲精品综合一区在线观看| 亚洲18禁久久av| 亚洲第一电影网av| 色综合亚洲欧美另类图片| 亚洲av不卡在线观看| 亚洲自拍偷在线| 亚洲人成伊人成综合网2020| 久久久久亚洲av毛片大全| 可以在线观看的亚洲视频| 成人特级黄色片久久久久久久| 欧美成人a在线观看| 欧美绝顶高潮抽搐喷水| 国产一区二区激情短视频| 成年女人毛片免费观看观看9| 国产在视频线在精品| 在线播放国产精品三级| 日韩精品中文字幕看吧| 91麻豆精品激情在线观看国产| 国产精华一区二区三区| 午夜日韩欧美国产| 国产伦精品一区二区三区视频9| 久久性视频一级片| 男女那种视频在线观看| 欧美+亚洲+日韩+国产| 男女之事视频高清在线观看| 老女人水多毛片| 免费看a级黄色片| 国产精品久久视频播放| 色综合婷婷激情| 国产私拍福利视频在线观看| 亚洲av不卡在线观看| 日韩欧美精品v在线| 免费在线观看亚洲国产| 亚洲真实伦在线观看| 狂野欧美白嫩少妇大欣赏| 一区二区三区高清视频在线| 成人性生交大片免费视频hd| 成年女人看的毛片在线观看| 搡老妇女老女人老熟妇| 亚洲在线观看片| 日韩精品青青久久久久久| 天美传媒精品一区二区| 一本精品99久久精品77| 欧美+亚洲+日韩+国产| 99国产精品一区二区蜜桃av| 亚洲三级黄色毛片| 亚洲精品成人久久久久久| 午夜精品在线福利| 国产精品国产高清国产av| 久久精品影院6| av女优亚洲男人天堂| 在线a可以看的网站| 国产精品自产拍在线观看55亚洲| 国产伦在线观看视频一区| 亚洲,欧美精品.| 毛片女人毛片| 欧美日韩乱码在线| 欧美成人免费av一区二区三区| 色哟哟哟哟哟哟| 成人特级av手机在线观看| 亚洲人成网站在线播| ponron亚洲| 最近最新免费中文字幕在线| 看片在线看免费视频| 国产免费av片在线观看野外av| 欧美黄色淫秽网站| 欧美bdsm另类| 尤物成人国产欧美一区二区三区| 亚州av有码| 午夜精品一区二区三区免费看| www.999成人在线观看| 中文字幕久久专区| 亚洲无线在线观看| 欧美乱色亚洲激情| 岛国在线免费视频观看| 免费av不卡在线播放| 日韩欧美精品免费久久 | 国产高清有码在线观看视频| 亚洲经典国产精华液单 | 日韩欧美 国产精品| 简卡轻食公司| 日本五十路高清| 麻豆av噜噜一区二区三区| 最好的美女福利视频网| 成年女人毛片免费观看观看9| 97人妻精品一区二区三区麻豆| 色综合欧美亚洲国产小说| 五月伊人婷婷丁香| 国产av一区在线观看免费| 我的女老师完整版在线观看| 国产探花在线观看一区二区| 一个人免费在线观看的高清视频| 内射极品少妇av片p| 岛国在线免费视频观看| 日本成人三级电影网站| 久久国产精品影院| 一边摸一边抽搐一进一小说| 亚洲第一区二区三区不卡| 男人和女人高潮做爰伦理| 欧美日韩黄片免| 精品熟女少妇八av免费久了| 最好的美女福利视频网| 中出人妻视频一区二区| 九九久久精品国产亚洲av麻豆| www.熟女人妻精品国产| 免费高清视频大片| 久久久久久久亚洲中文字幕 | 欧美成人一区二区免费高清观看| 国产精品美女特级片免费视频播放器| 又黄又爽又免费观看的视频| 琪琪午夜伦伦电影理论片6080| 久久99热6这里只有精品| 九九在线视频观看精品| 两个人视频免费观看高清| 在线国产一区二区在线| 五月玫瑰六月丁香| 丰满人妻熟妇乱又伦精品不卡| 国产午夜精品久久久久久一区二区三区 | 高潮久久久久久久久久久不卡| 国产一区二区在线av高清观看| 人人妻人人看人人澡| 亚洲av一区综合| 白带黄色成豆腐渣| 国产精华一区二区三区| 国产三级中文精品| 亚洲av不卡在线观看| 亚洲综合色惰| 日日摸夜夜添夜夜添小说| 噜噜噜噜噜久久久久久91| 欧美激情国产日韩精品一区| 亚洲精品影视一区二区三区av| 亚洲成人免费电影在线观看| а√天堂www在线а√下载| 精品久久久久久久人妻蜜臀av| 美女xxoo啪啪120秒动态图 | 国产真实乱freesex| 婷婷精品国产亚洲av在线| 极品教师在线视频| 丰满人妻一区二区三区视频av| 成人一区二区视频在线观看| 亚洲av一区综合| 黄色日韩在线| 别揉我奶头~嗯~啊~动态视频| 午夜免费成人在线视频| 亚洲男人的天堂狠狠| 久久久久久久久久黄片| 亚洲成av人片在线播放无| 岛国在线免费视频观看| 精品一区二区三区视频在线观看免费| 国产中年淑女户外野战色| 亚洲成人中文字幕在线播放| 欧美高清性xxxxhd video| 午夜免费成人在线视频| 国产亚洲av嫩草精品影院| 国产精品亚洲一级av第二区| 欧美成人性av电影在线观看| 在线观看舔阴道视频| 亚洲国产欧美人成| 中文字幕久久专区| 亚州av有码| 国产综合懂色| 久久这里只有精品中国| 蜜桃亚洲精品一区二区三区| av中文乱码字幕在线| 精品无人区乱码1区二区| 精品人妻1区二区| 一级作爱视频免费观看| 国产免费一级a男人的天堂| 在线免费观看不下载黄p国产 | 精品国产三级普通话版| 真实男女啪啪啪动态图| 丰满乱子伦码专区| 看十八女毛片水多多多| 成人av在线播放网站| 日韩成人在线观看一区二区三区| 乱码一卡2卡4卡精品| 欧美日韩中文字幕国产精品一区二区三区| 欧美日韩瑟瑟在线播放| 在现免费观看毛片| 深夜a级毛片| 国产aⅴ精品一区二区三区波| 欧美一区二区国产精品久久精品| 国产精品一区二区三区四区免费观看 | 亚洲av二区三区四区| 夜夜爽天天搞| 搞女人的毛片| 亚洲av一区综合| 成人永久免费在线观看视频| 乱人视频在线观看| 91在线观看av| a级毛片a级免费在线| 人人妻人人看人人澡| 午夜亚洲福利在线播放| 久久精品国产亚洲av涩爱 | 日本在线视频免费播放| 波多野结衣高清无吗| 久久午夜福利片| 搡老熟女国产l中国老女人| 午夜福利免费观看在线| 日日夜夜操网爽| 床上黄色一级片| 免费一级毛片在线播放高清视频| 99在线人妻在线中文字幕| 日韩大尺度精品在线看网址| 97超视频在线观看视频| 日日摸夜夜添夜夜添av毛片 | 国产久久久一区二区三区| 哪里可以看免费的av片| 天天一区二区日本电影三级| 好男人在线观看高清免费视频| 欧美日本亚洲视频在线播放| 国产又黄又爽又无遮挡在线| 深爱激情五月婷婷| 亚洲电影在线观看av| 亚洲欧美激情综合另类| 91麻豆精品激情在线观看国产| 人人妻人人看人人澡| 最近最新免费中文字幕在线| 97碰自拍视频| 高清在线国产一区| 中文字幕熟女人妻在线| 国产老妇女一区| 免费在线观看成人毛片| 最新中文字幕久久久久| 亚洲av中文字字幕乱码综合| 性插视频无遮挡在线免费观看| 精品不卡国产一区二区三区| 高清毛片免费观看视频网站| 国产极品精品免费视频能看的| 国产精品一及| 两个人视频免费观看高清| 免费看光身美女| 露出奶头的视频| 人妻丰满熟妇av一区二区三区| 无遮挡黄片免费观看| 久久久久性生活片| 此物有八面人人有两片| 18禁黄网站禁片午夜丰满| 波多野结衣高清无吗| 成人av在线播放网站| 精品人妻偷拍中文字幕| 免费观看的影片在线观看| 中文亚洲av片在线观看爽| 久久天躁狠狠躁夜夜2o2o| 免费无遮挡裸体视频| 国产精品一区二区三区四区久久| 中文资源天堂在线| 最近最新免费中文字幕在线| 1000部很黄的大片| 精品无人区乱码1区二区| 婷婷精品国产亚洲av在线| 国产成人欧美在线观看| 婷婷六月久久综合丁香| 成年女人永久免费观看视频| 精品欧美国产一区二区三| 国产精品一区二区性色av| 亚洲av免费高清在线观看| 欧美不卡视频在线免费观看| 丝袜美腿在线中文| 午夜两性在线视频| 成人av在线播放网站| 自拍偷自拍亚洲精品老妇| 免费在线观看亚洲国产| 亚洲自拍偷在线| 欧美黄色淫秽网站| 久久国产乱子伦精品免费另类| 一级av片app| 男女做爰动态图高潮gif福利片| av在线老鸭窝| 少妇的逼好多水| 搡老熟女国产l中国老女人| 欧美激情久久久久久爽电影| 国产一区二区三区在线臀色熟女| 亚洲,欧美,日韩| 欧美成人a在线观看| 一进一出抽搐gif免费好疼| 国内精品久久久久久久电影| 色噜噜av男人的天堂激情| 最近最新中文字幕大全电影3| 天堂√8在线中文| 国产午夜精品论理片| 亚洲精华国产精华精| 国产精品一区二区性色av| 免费在线观看亚洲国产| 亚洲片人在线观看| 九色国产91popny在线| 国产一区二区在线av高清观看| 91九色精品人成在线观看| 亚洲国产精品成人综合色| 亚洲人成网站在线播放欧美日韩| 一级作爱视频免费观看| 亚洲自拍偷在线| 亚洲精品乱码久久久v下载方式| 男女下面进入的视频免费午夜| 男女之事视频高清在线观看| netflix在线观看网站| 亚洲七黄色美女视频| 少妇熟女aⅴ在线视频| 国产伦人伦偷精品视频| 老熟妇乱子伦视频在线观看| 国产欧美日韩一区二区精品| 亚洲自拍偷在线| 亚洲在线自拍视频| 一夜夜www| 女同久久另类99精品国产91|