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

    基于機器學習的中西太平洋黃鰭金槍魚漁場預報模型

    2022-11-13 08:03:30周為峰唐峰華石永闖
    農(nóng)業(yè)工程學報 2022年15期
    關(guān)鍵詞:黃鰭漁場金槍魚

    張 聰,周為峰,唐峰華,石永闖,樊 偉

    基于機器學習的中西太平洋黃鰭金槍魚漁場預報模型

    張 聰1,2,周為峰1※,唐峰華1,石永闖1,樊 偉1

    (1. 中國水產(chǎn)科學研究院東海水產(chǎn)研究所,上海 200090;2. 中國農(nóng)業(yè)科學院研究生院,北京 100081)

    為提供準確的中西太平洋黃鰭金槍魚漁場預報信息,該研究利用2008-2019年中國水產(chǎn)集團43艘遠洋延繩釣漁船在中西太平洋海域(0°~30°S;110°E~170°W)作業(yè)的漁業(yè)數(shù)據(jù),通過方差膨脹因子篩選、歸一化處理,選取時空因子、海洋環(huán)境因子及大尺度氣候數(shù)據(jù)等共35種特征因子,構(gòu)建了一種隨機森林和極端梯度提升決策樹相結(jié)合的XGBRF模型,并利用五折交叉驗證法確定最佳參數(shù),選擇邏輯回歸、分類與回歸樹、最近鄰、自適應增強、梯度提升決策樹、極端梯度提升決策樹和隨機森林等模型作為對照,建立8種黃鰭金槍魚漁場預測模型并進行模型間的比較分析。結(jié)果表明,XGBRF模型對中西太平洋黃鰭金槍魚漁場的預測性能比其他模型更好,其準確率、漁場召回率、漁場F1得分、非漁場查準率和曲線下面積值A(chǔ)UC均最高,分別為75.39%、87.36%、82.64%、66.32%和79.48%,且模型的受試者工作特征曲線ROC更靠近左上角;海表溫度是影響中西太平洋黃鰭金槍魚漁場分布最重要的環(huán)境因子,其他因子依次是300 m水層溫度、50 m水層鹽度、葉綠素a濃度、南方濤動指數(shù)以及表層鹽度因子,時空因子和其余大尺度氣候因子的影響程度較低;基于XGBRF預報模型得到的漁場預測結(jié)果與實際作業(yè)范圍總體一致。XGBRF集成模型對中西太平洋海域黃鰭金槍魚的漁場預報具有較好的效果,可為漁場預報提供參考。

    機器學習;模型;中西太平洋;黃鰭金槍魚;漁場預報

    0 引 言

    黃鰭金槍魚()屬鱸形目、金槍魚屬,廣泛分布于熱帶和亞熱帶水域,中西太平洋是金槍魚漁業(yè)重要的作業(yè)區(qū)域[1]。延繩釣被認為是最合理的金槍魚捕撈方法之一。近年來,由于過度捕撈、漁業(yè)管理不當?shù)葐栴},漁場時空分布的穩(wěn)定性降低,加之金槍魚高度洄游的特性,金槍魚中心漁場范圍分散且分布廣泛[2]。精確的漁場預報模型能夠提高對潛在漁區(qū)和非漁區(qū)的識別能力,從而有效縮短尋漁時間,減少尋漁成本,提高漁獲產(chǎn)量,對金槍魚漁業(yè)的發(fā)展具有重要意義。

    魚類的生長發(fā)育與周圍的海洋環(huán)境密切相關(guān),利用不同的海洋環(huán)境因子建立漁場預報模型是當前漁場預報的常用手段[3]。貝葉斯(Bayesian)是一種常見的漁場預報模型建模方法,周為峰等[4]基于貝葉斯分類器,利用環(huán)境因子及不同分類策略構(gòu)建8種南海外海黃鰭金槍魚預報模型,最高準確率達75%?;貧w模型也很常見,如廣義相加模型(Generalized Additive Model,GAM)[5-6],Perez等[7]依據(jù)GAM模型,利用單位捕撈努力量指數(shù)預測菲律賓周圍海域的黃鰭金槍魚分布。Zagaglia等[8]基于GAM統(tǒng)計模型分析了巴西東北部黃鰭金槍魚延繩釣數(shù)據(jù)與海表溫度等環(huán)境變量之間的關(guān)系,證明了漁獲量與環(huán)境數(shù)據(jù)之間的非線性關(guān)系。隨著計算機技術(shù)的發(fā)展,人工神經(jīng)網(wǎng)絡(luò)(Artificial neural network, ANN)逐漸應用到黃鰭金槍魚漁場預報模型中,鄭志輝[9]基于人工神經(jīng)網(wǎng)絡(luò)模型有效預測了中西太平洋馬紹爾海域黃鰭金槍魚的漁場分布。棲息地指數(shù)模型(Habitat Suitability Index,HSI)近年來受到較多關(guān)注,趙海龍等[10]以產(chǎn)量來表征資源量,結(jié)合海表溫度、海面高度按季度建立東太平洋黃鰭金槍魚棲息地適應性指數(shù),預報的準確率達到66%以上。盡管上述漁場預報模型的結(jié)構(gòu)簡單便利,但這些模型大多是單一的預報模型且只使用海洋表面環(huán)境數(shù)據(jù)。由于海洋漁業(yè)信息具有模糊性、復雜性、不精確性的特點,單一模型容易過擬合而泛化性能差,無法滿足漁場預報應用對于精度的要求。且上述模型大多采用具體的海洋環(huán)境因子,對于與海洋因素密切相關(guān)的各種大尺度氣候數(shù)據(jù)的研究較少。

    針對數(shù)據(jù)復雜和模型單一的問題,為提供準確的中西太平洋黃鰭金槍魚漁場預報信息,本文提出隨機森林[11](Random Forest,RF)和極端梯度提升決策樹[12](eXtreme Gradient Boosting,XGBoost)相結(jié)合的XGBRF(XGBoost with RF)混合集成模型,根據(jù)2008-2019年中西太平洋黃鰭金槍魚的漁業(yè)數(shù)據(jù),結(jié)合31種特征因子和4種大尺度氣候因子,通過方差膨脹因子(Variance Inflation Factor,VIF)篩選、歸一化處理,建立了8種黃鰭金槍魚漁場預測模型。通過模型檢驗和對比,篩選預報精度較高、預測性能較好的模型,分析環(huán)境因子影響漁場分布的重要程度,以提高大數(shù)據(jù)和人工智能方法在黃鰭金槍漁業(yè)中的應用,為今后漁場預報的模型構(gòu)建提供參考。

    1 材料與方法

    1.1 數(shù)據(jù)來源

    以中西太平洋黃鰭金槍魚延繩釣漁船的作業(yè)范圍(110°E~170°W、0°~30°S)為研究海域,使用的漁業(yè)生產(chǎn)數(shù)據(jù)來自2008-2019年中國水產(chǎn)集團43艘遠洋延繩釣漁船漁撈日志,包括船名、作業(yè)日期(年/月/日)、作業(yè)位置(經(jīng)度、緯度)、漁獲信息(魚種、產(chǎn)量、尾數(shù)和下鉤數(shù)等)。

    葉綠素a濃度數(shù)據(jù)來源于美國國家海洋和大氣管理局數(shù)據(jù)庫(https://oceancolor.gsfc.nasa.gov/)。海面高度異常(Sea Level Abnormal,SLA)數(shù)據(jù)來源于AVISO(Archiving, Validation and Interpretation of Satellite Oceanographic data)數(shù)據(jù)庫(https://www.aviso.oceanobs.com)。渦動能(Eddy Kinetic Energy,EKE)、0~500 m水層的垂直溫度和鹽度數(shù)據(jù)來源于哥白尼海洋環(huán)境監(jiān)測服務(wù)中心網(wǎng)站(http://marine.copernicus.eu)。環(huán)境數(shù)據(jù)的時間分辨率為月,SLA、EKE和0~500 m水層垂直溫度和鹽度數(shù)據(jù)的空間分辨率為0.25°×0.25°,Chla的空間分辨率為4 km。使用python將環(huán)境數(shù)據(jù)的空間分辨率統(tǒng)一為0.5°×0.5°的空間格網(wǎng)并與漁獲量數(shù)據(jù)進行匹配。

    大尺度氣候數(shù)據(jù)均為月度數(shù)據(jù),南方濤動指數(shù)(Southern Oscillation Index,SOI)和北極濤動指數(shù)(Arctic Oscillation Index,AOI)來源于美國海洋與大氣局氣候預測中心;太平洋年代際振蕩指數(shù)(Pacific Decadal Oscillation Index,PDOI)來源于華盛頓大學(http://research.jisao.washington.edu /pdo),北太平洋環(huán)流振蕩指數(shù)(North Pacific Gyre Oscillation Index,NPGOI)來源于http://www.o3d.org /npgo /。

    1.2 數(shù)據(jù)處理

    1.2.1 CPUE計算

    將漁區(qū)網(wǎng)格以0.5°×0.5°進行劃分,按月對其作業(yè)位置、尾數(shù)和放鉤數(shù)進行統(tǒng)計,計算各漁區(qū)網(wǎng)格內(nèi)的黃鰭金槍魚單位捕撈努力量漁獲量[13](Catch Per Unit Effort,CPUE),公式為

    1.2.2 衍生環(huán)境因子計算

    本文選取的衍生環(huán)境因子包括海表溫度距平(SSTdt)、海表溫度梯度(SSTG)、葉綠素距平(Chladt)和葉綠素梯度(ChlaG)。SSTdt,j、Chladt,j、SSTG,j和ChlaG,j計算公式如下:

    1.2.3 CPUE與各環(huán)境因子的相關(guān)性分析

    采用Python的seaborn包計算Pearson相關(guān)系數(shù)并進行顯著性檢驗,結(jié)果如表1所示。假設(shè)顯著性水平為0.05(即超過95%的概率有相關(guān)關(guān)系),本文中CPUE與各指標均呈現(xiàn)顯著相關(guān)關(guān)系(<0.05)。

    1.2.4 數(shù)據(jù)標準化與共線性分析

    由于漁業(yè)數(shù)據(jù)和海洋環(huán)境數(shù)據(jù)具有不同的量綱和數(shù)量級,直接用原始指標值進行分析,會突出數(shù)值較高的指標在綜合分析中的作用,為了保證結(jié)果的可靠性,需要對原始數(shù)據(jù)進行標準化處理。本文采用數(shù)據(jù)歸一化方法將各個特征映射到[0,1]的范圍,計算公式為

    式中、、max、min分別為歸一化后的值、實際值、最大值、最小值。

    對于環(huán)境因子之間可能存在共線性的問題,本文采用方差膨脹因子VIF進行篩選,VIF的計算公式如下:

    式中是第個變量在所有其他變量上回歸時的復相關(guān)系數(shù)。各環(huán)境變量之間的Person相關(guān)系數(shù)矩陣圖如圖1所示。

    表1 CPUE與變量間的相關(guān)系數(shù)和P值

    注:Lat表示緯度;Lon表示經(jīng)度;Chla表示葉綠素a濃度;Chladt表示葉綠素距平;SSTdt表示海表溫度距平;SSTG表示海表溫度梯度;ChlaG表示葉綠素梯度;SLA表示海面高度異常;EKE表示渦動能;S、T(=0,50,100,…,500)分別表示對應深度處的鹽度和溫度;PDOI為太平洋年代際振蕩指數(shù);SOI為南方濤動指數(shù);AOI為北極濤動指數(shù);NPGOI為北太平洋環(huán)流振蕩指數(shù)。

    Note: Lat indicates latitude; Lon indicates longitude; Chla denotes chlorophyll a concentration; Chladt indicates chlorophyll anomaly; SSTdt indicates sea surface temperature anomaly; SSTG indicates sea surface temperature gradient; ChlaG indicates chlorophyll gradient; SLA indicates abnormal sea level; EKE denotes eddy kinetic energy; Sand T(= 0,50, 100, ..., 500) respectively represent the salinity and temperature at the corresponding depth; PDOI is the Pacific interdecadal oscillation index; SOI is the Southern Oscillation Index; AOI is the Arctic oscillation index; NPGOI is the North Pacific circulation oscillation index.

    圖1 各變量間Pearson相關(guān)系數(shù)

    1.3 XGBRF漁場預報模型的建立

    隨機森林算法結(jié)合自助采樣和隨機子空間思想對決策樹進行集成,可降低使用單一決策樹時可能出現(xiàn)的較大誤差和過擬合現(xiàn)象,從而得到更加穩(wěn)定和準確的預測效果。XGBoost算法作為一種可擴展的樹提升系統(tǒng),通過正則化項、二階泰勒展開等策略控制模型的復雜度并提高運行效率,能有效應對高維度的復雜數(shù)據(jù)。

    由于漁業(yè)數(shù)據(jù)數(shù)量龐大、特征維度高,考慮到隨機森林并行集成、不易過擬合的優(yōu)勢及XGBoost算法在分類問題中的優(yōu)秀性能,本文提出XGBRF混合集成模型,將boosting與bagging思想相結(jié)合,選擇XGBoost算法作為預報模型的基本框架來訓練隨機森林。本文中XGBRF模型的算法流程如圖2所示,參數(shù)設(shè)置如表2所示。

    注:L為損失函數(shù),P′為迭代次數(shù),T為樹的數(shù)量,θ為一個隨機向量,b(x)為基學習器。訓練集D由M個樣本組成,D=(xi, yi),i=1,2,3,…,M。G為一階導數(shù)和,H為二階導數(shù)和。

    表2 XGBRF集成模型中的參數(shù)設(shè)置

    2008-2019年中西太平洋黃鰭金槍魚的延繩釣數(shù)據(jù)共18 013條。本研究采用五折交叉方法隨機劃分為訓練數(shù)據(jù)(80%數(shù)據(jù))和測試數(shù)據(jù)(20%數(shù)據(jù))。對數(shù)據(jù)集進行預處理,將單位捕撈努力量按照三分位數(shù)劃分為0和1,分別對應非漁場和漁場。為了驗證模型的有效性,以邏輯回歸(Logistic Regression,LR)、分類與回歸(Classification and regression tree,CART)、近鄰(K Nearest Neighbors,KNN)、自適應增強(Adaptive boosting,Adaboost)、梯度提升樹(Gradient Boosting Decision Tree,GBDT)、XGBoost和RF共7種模型作為對照,對照模型均在相同的實驗環(huán)境下使用同一數(shù)據(jù)集。為了消除由于數(shù)據(jù)劃分造成模型訓練的性能誤差,本文中所使用的模型均采用網(wǎng)格搜索和兩次五折交叉驗證訓練法以得到最優(yōu)參數(shù)。各模型的參數(shù)設(shè)置如表3所示。

    表3 模型參數(shù)設(shè)置

    1.4 模型精度評價

    模型精度是判別預報結(jié)果與真實漁場信息符合程度的標準。本文使用準確率(Accuracy)、查準率(Precision)、召回率(Recall)、F1得分(F1-score)、受試者工作特征曲線(Receiver Operating Characteristic,ROC)及曲線下的面積(Area under Curve,AUC)值綜合評價模型的性能。

    將模型中少數(shù)類樣本記為負例(非漁場),多數(shù)類樣本記為正例(漁場),則各個評價指標的計算公式為

    式中TP、TN、FP、FN分別代表預測值和真實值均為漁場、預測值和真實值均為非漁場、預測值為漁場而真實值為非漁場以及預測值為非漁場而真實值為漁場。ROC曲線由真正率和假正率分別作為縱坐標和橫坐標繪制得到,AUC值是ROC曲線與橫坐標圍成的面積,值域在 0~1之間。

    2 結(jié)果與分析

    2.1 基于VIF的特征選擇

    基于方差膨脹因子VIF的特征選擇過程中,剔除環(huán)境因子中VIF較高的值后重新計算剩余變量的VIF值,重復操作直至剩余變量的VIF值均小于10。VIF選擇后的變量如表4所示。方差膨脹因子篩選后的變量共19種,除了SLA、S0和T150,其余環(huán)境因子的VIF值均小于5,可作為模型輸入變量。

    表4 特征選擇后的環(huán)境變量VIF值

    2.2 預測結(jié)果及模型對比分析

    將20%的測試數(shù)據(jù)分別代入LR、CART、KNN、Adaboost、GBDT、XGBoost、RF和XGBRF模型,分別計算各個模型的Accuracy、Precision、Recall、F1-score、AUC值,結(jié)果如表5所示。

    表5 各模型預測結(jié)果對比

    由表5可知,XGBRF模型在準確率和AUC的值均最高,分別為75.39%和79.48%,證明模型的預測效果較好。對于漁場和非漁場的預測結(jié)果,8種模型的漁場召回率、查準率和F1-score相對于非漁場均較高,可以看出所有模型對于漁場的識別均要好于非漁場。RF模型對于非漁場的識別效果較好,其非漁場F1-score(60.27%)最高,XGBRF模型表現(xiàn)次之,其非漁場查準率(66.32%)最高。對于遠洋漁業(yè)而言,漁船的作業(yè)圍繞中心漁場進行以提高漁獲產(chǎn)量和捕撈效率,因此中心漁場的預報精度具有重要意義,XGBRF模型的漁場查準率略低于RF模型,但其召回率(87.36%)和F1-score(82.64%)均最高,證明模型對于真實漁場的識別精度高,XGBRF模型的性能總體上優(yōu)于其他模型。

    各個模型的ROC曲線如圖3所示,由圖3可知,XGBRF模型較其他模型的ROC曲線更靠近左上角,能夠較好地包含其他模型,證明模型的分類和預測效果較好。

    圖3 不同模型的ROC曲線

    2.3 模型因子的重要性分析

    借助sklearn中的feature_importances_方法,得到入模因子對中西太平洋黃鰭金槍魚漁場預報的重要程度圖(圖4)。由圖4可知,研究結(jié)果顯示海表溫度是影響黃鰭金槍魚漁場布的最主要因素,為7.573%,300 m水層溫度、50 m水層鹽度、葉綠素a濃度(Chla)、南方濤動指數(shù)(SOI)和表層鹽度(S0)對黃鰭金槍魚的影響也較大,分別為7.369%、6.765%、6.533%、6.284%和6.212%。除SOI以外,其余大尺度氣候因子的影響相對較低。時空因子的相對重要性較環(huán)境因子低。

    圖4 模型輸入因子的相對重要性排序

    2.4 模型驗證

    模型利用2008-2019年的真實數(shù)據(jù)進行訓練和預測。將20%測試集的實際數(shù)據(jù)與模型預測的數(shù)據(jù)疊加,結(jié)果如圖5所示。

    由圖5b可知,測試的漁場主要分布在10°~20°S,155°~175°E海域,這一海域模型的識別程度較高,但也存在一定的誤判。在27°S海域附近也有少量漁場分布,但該海域的識別程度較低。非漁場的實際與預測結(jié)果如圖5a所示,非漁場的分布與漁場的分布范圍相似,在15°S以下的海域具有較好的識別度,但非漁場的總體識別率低于漁場的識別率。總體而言,預測的漁場位置與實際漁場位置存在少量的偏差,相較于其他模型,預測結(jié)果的精確度與可信度高。

    注:審圖號GS(2016)1665。

    3 討 論

    3.1 環(huán)境因子選擇的合理性

    海洋環(huán)境因子是影響黃鰭金槍魚活動和空間分布的重要外在因素,現(xiàn)有的研究表明,黃鰭金槍魚漁場的形成極易受到海洋環(huán)境因子的影響[5,15-16]。葉綠素a濃度通過海洋食物鏈原理影響漁場分布;溫度是影響魚類活動最基本的海洋環(huán)境要素之一[17];海平面高度異常作為與平均海平面的差值[18],反映了包括海流、冷暖水團、鹽度、海洋鋒面等海洋動力特征,常被應用于海洋漁場分析;渦動能是通過影響環(huán)流、海水溫度以及葉綠素a的分布,從而影響黃鰭金槍魚的漁場分布;研究表明,黃鰭金槍魚的垂直移動現(xiàn)象明顯[19-20],其主要的垂直移動范圍從幾十米至幾百米,且不同水層的溫度會影響黃鰭金槍魚漁場的分布;南方濤動伴隨季節(jié)上升流的延遲可能與某些海域浮游植物的繁盛有關(guān)[21];太平洋年代際振蕩被證明與SST具有一定的關(guān)聯(lián)[22];北太平洋環(huán)流振蕩、北極濤動等其他大尺度氣候也有助于全面的了解海洋環(huán)境的變化[23]。因此,本研究選取的環(huán)境因子包括葉綠素a濃度、海面高度異常、渦動能、0~500 m水層的垂直溫度和鹽度數(shù)據(jù),為了增大樣本間差異,提高模型對中心漁場的識別能力,計算了SST和Chla的衍生變量(SSTdt、SSTG、Chladt和ChlaG)。

    影響漁業(yè)的海洋環(huán)境因子關(guān)系錯綜復雜,多個特征變量之間存在較強的相關(guān)性(圖1),研究表明相關(guān)系數(shù)大于0.85[24],方差膨脹因子大于10[25],變量之間可能存在共線性,造成重要變量的預測結(jié)果不顯著[26]。本研究綜合考慮各個海洋環(huán)境因子及共線性因素對模型精度和效率的影響,最終確定的入模變量具有一定的科學性和合理性。

    3.2 模型因子的相對重要性

    海表溫度通過影響黃鰭金槍魚的生長、覓食、洄游遷徙等直接影響漁場的分布和變動,研究表明絕大多數(shù)黃鰭金槍魚可至10 m以內(nèi)的淺層活動[20],結(jié)果表明海表溫度是最重要的環(huán)境因子。300 m水層溫度對黃鰭金槍魚的影響同樣重要,這可能與黃鰭金槍魚的垂直活動且主要集中在300 m深度水層有關(guān)[27],除此之外S50、S0、S150、T150等對黃鰭金槍魚的影響也間接證明了其垂直游動的特性。葉綠素a濃度作為海洋初級生產(chǎn)力的指標,其對黃鰭金槍魚漁場分布的影響已被多人研究證實,且Chla的梯度分布(ChlaG)與鋒面的形成有關(guān),也與中心漁場的分布存在聯(lián)系[28]。大尺度氣候數(shù)據(jù)中,南方濤動指數(shù)SOI對黃鰭金槍魚的影響相對較大,這與Torres-Faurrieta[29]的研究結(jié)果一致,NPGOI、PDOI、AOI的影響相對較小,這可能與大尺度氣候影響的滯后效應有關(guān)[30]。

    3.3 漁場與非漁場的劃分

    三分位數(shù)是統(tǒng)計學研究中最常用的方法之一。在進行漁場分類的研究中,對于“中心漁場”的界定并沒有一個統(tǒng)一的標準,已有的研究大多采用三分位數(shù)對漁場進行劃分[11,31-32],部分學者利用三分位數(shù)將漁區(qū)劃分為高、中、低3個漁區(qū),基于預報模型進行分別預測。但是需要考慮到的是CPUE作為連續(xù)變量,處于三分位數(shù)邊界區(qū)域的值往往大小相近,會給模型的預測造成一定的困難,且各種復雜環(huán)境要素如氣候異常、海洋大尺度事件、漁類洄游等的影響進一步導致中產(chǎn)漁區(qū)的預測準確率較低[11,32-33],造成分類損失。本研究結(jié)合真實漁業(yè)數(shù)據(jù),將大于整年單位捕撈努力量第一三分位數(shù)T1的區(qū)域定義為漁場,其余部分作為非漁場,減小由于臨界區(qū)域造成的模型性能損失,實現(xiàn)簡單且具有一定的實際意義。

    3.4 模型性能對比

    本文中XGBRF模型得到的預報準確率為75.39%,漁場召回率和AUC值分別為87.36%和79.48%,預報效果優(yōu)于棲息地指數(shù)模型[10](準確率66%)、GAM模型[6](解釋率37.2%)、stacking算法[34](準確率68.72%)、貝葉斯分類器[4](準確率75%),且實驗結(jié)果顯示XGBRF模型相比于RF和XGBoost模型在各個指標上的表現(xiàn)均有一定程度的提升,這是因為XGBRF模型結(jié)合了boosting與bagging思想的優(yōu)點,XGBoost作為基本框架能夠針對錯誤樣本進行學習,顯著提高學習效率,同時隨機采樣和隨機特征選擇的加入增大了模型的靈活性,使其具有更好的泛化性能。LR和CART屬于簡單模型,簡單直觀、解釋性強,但LR容易欠擬合,CART容易過擬合,兩種模型的預測精度均不高。KNN通過計算不同數(shù)據(jù)之間的歐式距離來進行分類,準確度高,對異常點不敏感,在本研究中的表現(xiàn)較好。Boosting方法能夠針對樣本進行學習進而提高學習效率,Adaboost利用錯誤分類的樣本來識別問題;GBDT通過計算負梯度來擬合殘差,該類模型能充分學習每個弱分類器,具有很高的精度,結(jié)果顯示兩者在本研究中的表現(xiàn)互有優(yōu)劣。XGBoost算法在GBDT的基礎(chǔ)上進行了優(yōu)化,模型的效率更高、泛化性能更強、應用范圍更廣,研究結(jié)果也證明其效果好于Adaboost和GBDT。RF作為bagging的升級,不僅具有Bagging并行集成的優(yōu)勢,同時其基學習器的決策邊界往往更加靈活,有效解決了單一決策樹的性能瓶頸,具有較高的泛化能力。實驗中RF模型的漁場和非漁場召回率均較高,間接證明了模型的泛化能力。

    3.5 模型精度評價

    模型的評價指標是直接判斷一個模型好壞的關(guān)鍵,漁場預報發(fā)展至今,各種預報模型層出不窮,需要對模型進行合理合適的評估。準確率是現(xiàn)有漁場預報模型采用較多的評估依據(jù)[31-32,34],準確率表示模型預測正確的正例、負例占總樣本的比例,易于計算且復雜度低。但準確率是一種偏向全局化的衡量標準,信息量少且當模型對于非漁場(少數(shù)類)的識別率較高而對于漁場的識別率較低時也可能帶來準確率的偏高,這種模型往往不是我們想要的。查準率表示正例的分類準確率,表明被預測為正例的樣本中真實的正例的比例,查準率越高,表明漏掉的漁場就越少。召回率表示被預測正確的正例占所有真實的正例比例,召回率越高,代表漁場的識別率就越高。在漁業(yè)生產(chǎn)中,如果想盡可能覆蓋某片海域內(nèi)的所有漁場可以選擇查準率較高的模型;如果想盡量避免無獲而返的情況可以選擇召回率較高的模型。實際情況下,需要在二者中尋找平衡,F(xiàn)1得分作為查準率和查全率的綜合考慮解決了這一問題。F1得分被認為在某些方面表現(xiàn)優(yōu)于準確率[35],F(xiàn)1得分越高,表明漁場的識別率高且誤報的漁場較少。受試者工作特征曲線ROC曲線被認為衡量模型性能的有用指標[36],ROC曲線越靠近左上角表明模型的預測能力越強,AUC值大于0.75的模型被認為是“有用”的[37]。本文綜合考慮各個模型的特點,選擇準確率、查準率、召回率、F1-score、ROC曲線和AUC值作為模型的綜合評判指標,確保了模型結(jié)果的可靠性。

    4 結(jié) 論

    本研究利用中西太平洋黃鰭金槍魚的漁業(yè)生產(chǎn)數(shù)據(jù),結(jié)合漁場時空因子、海洋環(huán)境因子及大尺度氣候數(shù)據(jù)作為訓練集構(gòu)建了邏輯回歸(LR)、分類與回歸分類與回歸(CART)、近鄰(KNN)、自適應增強(Adaboost)、梯度提升樹(GBDT)、極端梯度提升決策樹(XGBoost)、隨機森林(RF)和XGBRF共8種黃鰭金槍魚漁場預報模型,并利用真實的漁業(yè)數(shù)據(jù)進行了模型驗證,研究結(jié)果如下:

    1)通過多個指標結(jié)果的綜合分析顯示,XGBRF混合集成模型能在一定程度上提高XGBoost和RF模型的性能,其準確率和AUC的值均最高,分別為75.39%和79.48%,且XGBRF模型的ROC曲線能較好包含其他模型,證明模型的性能表現(xiàn)最好,優(yōu)于傳統(tǒng)的預報模型。

    2)模型因子的重要性分析結(jié)果表明,海表溫度是影響黃鰭金槍魚漁場分布最重要的因素,為7.573%。接下來的是300 m水層的溫度,其重要性達7.369%。除此之外,50 m水層鹽度、南方濤動指數(shù)、葉綠素a濃度和表層鹽度對黃鰭金槍魚的影響也較大。除SOI以外,其余大尺度氣候因子的影響相對較低。

    3)20%測試集實際漁場和非漁場的疊加圖顯示,XGBRF模型預報得到的漁場位置與實際的漁場位置僅存在少量的偏差,預測結(jié)果的精確度和可信度較高,是一種有效的黃鰭金槍魚漁場預報方法。

    本研究構(gòu)建了一種隨機森林和極端梯度提升決策樹結(jié)合的XGBRF模型作為中西太平洋黃鰭金槍魚的漁場預報方法,通過與多種模型對比和實際漁業(yè)數(shù)據(jù)的檢驗,證明了該模型的可行性。在今后的研究中,可以進一步探究其他環(huán)境因子、模型參數(shù)優(yōu)化、漁區(qū)劃分標準、海洋模式數(shù)據(jù)等因素對黃鰭金槍魚漁場的影響,進一步提高模型的精度和實用性。

    [1] 張春玲,蔣逸,王冰洋,等. 中西太平洋黃鰭金槍魚隨附魚群區(qū)域溫度垂直結(jié)構(gòu)的構(gòu)建與初步分析[J]. 上海海洋大學學報,2022,31(1):233-241.

    ZhangChunling,JiangYi,Wang Bingyang, et al. Construction and analysis of vertical temperature structure for FAD yellowfin tuna in the central and western Pacific Ocean[J]. Journal of Shanghai Ocean University, 2022, 31(1): 233-241. (in Chinese with English abstract)

    [2] Ms A, Gmp B, Cr A, et al. Trade-offs for the southern longline fishery in achieving a candidate South Pacific albacore target reference point-ScienceDirect[J]. Marine Policy, 2019, 100(C): 66-75.

    [3] 朱浩朋,伍玉梅,唐峰華,等. 采用卷積神經(jīng)網(wǎng)絡(luò)構(gòu)建西北太平洋柔魚漁場預報模型[J]. 農(nóng)業(yè)工程學報,2020,36(24):153-160.

    Zhu Haopeng, Wu Yumei, Tang Fenghua, et al. Construction of fishing ground forecast model of Ommastrephes bartramii using convolutional neural network in the Northwest Pacific[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(24): 153-160. (in Chinese with English abstract)

    [4] 周為峰,黎安舟,紀世建,等. 基于貝葉斯分類器的南海黃鰭金槍魚漁場預報模型[J]. 海洋湖沼通報,2018(1):116-122.

    Zhou Weifeng, Li Anzhou, Ji Shijian, et al. Forecasting model for yellowfin tuna () fishing ground in the south China sea based on bayes classifier[J]. Transactions of Oceanology and Limnology, 2018(1): 116-122. (in Chinese with English abstract)

    [5] 宋利明,沈智賓,周建坤,等. 庫克群島海域海洋環(huán)境因子對黃鰭金槍魚漁獲率的影響[J]. 上海海洋大學學報,2016,25(3):454-464.

    Song Liming, Shen Zhibin, Zhou Jiankun, et al. Effects of environmental variables on catch rates of yellowfin tuna () in waters near Cook Islands[J]. Journal of Shanghai Ocean University, 2016, 25(3): 454-464. (in Chinese with English abstract)

    [6] 徐國強,朱文斌,張洪亮,等. 基于GAM模型分析印度洋大眼金槍魚和黃鰭金槍魚漁場分布與不同環(huán)境因子關(guān)系[J]. 海洋學報,2018,40(12):68-80.

    Xu Guoqiang, Zhu Wenbing, Zhang Hongliang, et al. Relationship between fishing grounds ofandwith environmental factors in the Indian Ocean based on generalized additive model[J]. Haiyang Xuebao, 2018, 40(12): 68-80. (in Chinese with English abstract)

    [7] Perez G J P, Leonardo E M. Predicting the Distribution of Yellowfin Tuna in Philippine Waters[C]. Washington DC: American Geophysical Union, 2015.

    [8] Zagaglia C R, Lorenzzetti J A, Stech José L. Remote sensing data and longline catches of yellowfin tuna () in the equatorial Atlantic[J]. Remote Sensing of Environment, 2004, 93(1/2): 267-281.

    [9] 鄭志輝. 基于貝葉斯方法、分位數(shù)回歸和人工神經(jīng)網(wǎng)絡(luò)模型的黃鰭金槍魚漁情預報模型比較[D]. 上海:上海海洋大學,2017.

    Zheng Zhihui. A Comparison of Yellowfin Tuna Fishing Condition Forecast Models Based on Bayesian Method, Quantile Regression and Artificial Neural Network Model[D]. Shanghai: Shanghai Ocean University, 2017. (in Chinese with English abstract)

    [10] 趙海龍,陳新軍,方學燕. 基于棲息地指數(shù)的東太平洋黃鰭金槍魚漁場預報[J]. 生態(tài)學報,2016,36(3):778-785.

    Zhao Hailong, Chen Xinjun, Fang Xueyan. Forecasting fishing ground of yellowfin tuna in the eastern Pacific Ocean based on the habitat suitability index[J]. Acta Ecologica Sinica, 2016, 36(3): 778-785. (in Chinese with English abstract)

    [11] 陳雪忠,樊偉,崔雪森,等. 基于隨機森林的印度洋長鰭金槍魚漁場預報[J]. 海洋學報, 2013,35(1):158-164.

    Chen Xuezhong, Fan Wei, Cui Xuesen, et al. Fishing ground forecasting of Thunnus alalung in Indian Ocean based on random forest[J]. Haiyang Xuebao, 2013, 35(1): 158-164. (in Chinese with English abstract)

    [12] Breiman L. Random forests[J]. Machine Learning, 2001, 45(1): 5-32.

    [13] Feng Y, Chen X, Gao F, et al. Impacts of changing scale on Getis-Ord Gi* hotspots of CPUE: A case study of the neon flying squid () in the northwest Pacific Ocean[J]. Acta Oceanologica Sinica, 2018, 37(5): 67-76.

    [14] Ismail A I, Morrison E C, Burt B A, et al. Natural history of periodontal disease in adults: findings from the Tecumseh Periodontal Disease study, 1959-87[J]. Journal of Dental Research, 1990, 69(2): 430-435.

    [15] 王少琴,許柳雄,朱國平,等. 中西太平洋金槍魚圍網(wǎng)的黃鰭金槍魚CPUE時空分布及其與環(huán)境因子的關(guān)系[J]. 大連海洋大學學報,2014,29(3):303-308.

    Wang Shaoqin, Xu Liuxiong, Zhu Guoping, et al. Spatial-temporal profiles of CPUE and relations to environmental factors for yellowfin tuna () from purse-seine fishery in Western and Central Pacific Ocean[J]. Journal of Dalian Ocean University, 2014, 29(3): 303-308. (in Chinese with English abstract)

    [16] 劉勇,陳新軍. 中西太平洋金槍魚圍網(wǎng)黃鰭金槍魚產(chǎn)量的時空分布及與表溫的關(guān)系[J]. 海洋漁業(yè),2007(4):296-301.

    Liu Yong, Chen Xinjun. Spatio-temporal distribution ofand its relationship with sea surface temperature in the tuna purse seine fishery of the Central and Western Pacific[J]. Marine Fisheries, 2007(4): 296-301. (in Chinese with English abstract)

    [17] 崔雪森,樊偉,張晶. 太平洋黃鰭金槍魚延繩釣漁獲分布及漁場水溫淺析[J]. 海洋通報,,2005(5):54-59.

    Cui Xuesen, Fan Wei, Zhang Jing. Primary study on the distribution of yellowfin tuna () longline fishery catch and the fishing-ground environment of the Pacific Ocean[J]. Marine Science Bulletin, 2005(5): 54-59. (in Chinese with English abstract)

    [18] 宋婷婷,樊偉,伍玉梅. 衛(wèi)星遙感海面高度數(shù)據(jù)在漁場分析中的應用綜述[J]. 海洋通報,,2013,32(4):474-480.

    Song Tingting, Fan Wei, Wu Yumei. Review on the fishery analysis with satellite remote sensing sea surface height data[J]. Marine Science Bulletin, 2013, 32(4): 474-480. (in Chinese with English abstract)

    [19] Cayre P. Behaviour of Yellowfin Tuna () and Skipjack Tuna () around fish aggregating devices (FADs) in the Comoros Islands as determined by ultrasonic tagging[J]. Aquatic Living Resources, 1991, 4(1): 1-12.

    [20] 張衡,戴陽,楊勝龍,等. 基于分離式衛(wèi)星標志信息的金槍魚垂直移動特性[J]. 農(nóng)業(yè)工程學報,2014,30(20):196-203.

    Zhang Heng, Dai Yang, Yang Shenglong, et al. Vertical movement characteristics of tuna () in Pacific Ocean determined using pop-up satellite archival tags[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(20): 196-203. (in Chinese with English abstract)

    [21] 王靚. 西北太平洋主要小型中上層魚類資源量對大尺度氣候-海洋環(huán)境變化的響應研究[D]. 上海:上海海洋大學,2021.

    Wang Jing. Study on the Response of Major Small Pelagic Fish Resource in the Northwest Pacific to Large Scale Climate-Marine Environmental Changes[D]. Shanghai: Shanghai Ocean University, 2021. (in Chinese with English abstract)

    [22] Tian Y J, Akamine T, Suda M. Variations in the abundance of Pacific saury () from the northwestern Pacific in relation to oceanic-climate changes[J]. Fisheries Research, 2003, 60(2/3): 439-454.

    [23] Vimont D J. The Contribution of the interannual ENSO cycle to the spatial pattern of decadal ENSO-like variability[J]. Journal of Climate, 2010, 18(12): 2080-2092.

    [24] Elith J, Graham C H, Anderson R P. et al. Novel methods improve prediction of species' distributions from occurrence data[J]. Ecography (Copenhagen), 2006, 29(2): 129-151.

    [25] Dormann C F, Elith J, Bacher S, et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance[J]. Ecography, 2013, 36(1): 27-46.

    [26] Ohlemuller R, Anderson B J, Araujo M B, et al. The coincidence of climatic and species rarity: High risk to small-range species from climate change[J]. Biology letters, 2008, 4(5): 568-572.

    [27] 楊勝龍,張忭忭,張衡,等. 黃鰭金槍魚垂直移動及水層分布研究進展[J]. 水產(chǎn)科學,2019,38(1):119-126.

    Yang Shenglong, Zhang Bianbian, Zhang Heng, et al. A review:vertical swimming and distribution of yellowfin tuna (thunnus albacares)[J]. Fisheries Science, 2019, 38(1): 119-126. (in Chinese with English abstract)

    [28] 毛志華,朱乾坤,龔芳. 衛(wèi)星遙感北太平洋漁場葉綠素a濃度[J]. 水產(chǎn)學報,2005(2):270-274.

    Mao Zhihua, Zhu Qiankun, Gong Fang. Satellite remote sensing of chlorophyll a concentration in the north Pacific Fishery[J]. Journal of Fisheries of China, 2005(2): 270-274. (in Chinese with English abstract)

    [29] Torres-Faurrieta L K, Dreyfus-León M J, Rivas D. Recruitment forecasting of yellowfin tuna in the eastern Pacific Ocean with artificial neuronal networks[J]. Ecological Informatics, 2016, 36: 106-113.

    [30] Báez J C, Czerwinski I A, Ramos M L. Climatic oscillations effect on the yellowfin tuna () Spanish captures in the Indian Ocean[J]. Fisheries Oceanography, 2020, 29(6): 572-583.

    [31] 侯娟,周為峰,樊偉,等. 基于集成學習的南太平洋長鰭金槍魚漁場預報模型研究[J]. 南方水產(chǎn)科學,2020,16(5):42-50.

    Hou Juan, Zhou Weifeng, Fan Wei, et al. Research on fishing grounds forecasting models of albacore tuna based on ensemble learning in South Pacific[J]. South China Fisheries Science, 2020, 16(5): 42-50. (in Chinese with English abstract)

    [32] 宮鵬,王德興,袁紅春,等. 基于 LightGBM 的南太平洋長鰭金槍魚漁場預報模型研究[J]. 水產(chǎn)科學, 2021,40(5):762-767.

    Gong Peng, Wang Dexing, Yuan Hongchun, et al. Fishing ground forecast model of albacore tuna based on LightGBM in the South Pacific Ocean[J]. Fisheries Science, 2021, 40(5): 762-767. (in Chinese with English abstract)

    [33] 袁紅春,陳冠奇,張?zhí)祢?,? 基于全卷積網(wǎng)絡(luò)的南太平洋長鰭金槍魚漁場預報模型[J]. 江蘇農(nóng)業(yè)學報,2020,36(2):423-429.

    Yuan Hongchun, Chen Guanqi, Zhang Tianjiao, et al. Fishing ground forecast model of albacore tuna based on fully convolutional networks in the South Pacific[J]. Jiangsu Journal of Agricultural Sciences, 2020, 36(2): 423-429. (in Chinese with English abstract)

    [34] 宋利明,任士雨,張敏,等. 基于集成學習的大西洋熱帶海域黃鰭金槍魚漁情預報[J]. 中國水產(chǎn)科學,2021,28(8):1069-1078.

    Song Liming, Ren Shiyu, Zhang Min, et al. Fishing ground forecasting models for yellowfin tuna () in the tropical waters of the Atlantic Ocean based on ensemble learning[J]. Journal of Fishery Sciences of China, 2021, 28(8): 1069-1078. (in Chinese with English abstract)

    [35] Hossin M, Sulaiman M N. A review on evaluation metrics for data classification evaluations[J]. International Journal of Data Mining & Knowledge Management Process, 2015, 5(2): 1-11.

    [36] Guo H, Li Y, Shang J, et al. Learning from class-imbalanced data: Review of methods and applications[J]. Expert Systems with Applications, 2017, 73: 220-239.

    [37] Swets J A. Measuring the accuracy of diagnostic systems[J]. Science, 1988, 240(4857):1285-1293.

    Forecasting models for yellowfin tuna fishing ground in the central and western Pacific based on machine learning

    Zhang Cong1,2, Zhou Weifeng1※, Tang Fenghua1, Shi Yongchuang1, Fan Wei1

    (1.200090; 2.100081)

    An accurate forecast can be greatly contributed to the yellowfin tuna fishing ground in the western and Central Pacific. However, a large amount of fishery data, and high feature dimension have posed a great over-fitting on the various classification in recent years. The random forest parallel integration can be expected to achieve the excellent performance of the extreme gradient boosting decision tree algorithm. In this study, a hybrid integration model was proposed to combine the Xgboost with Random Forest (XGBRF) with the random forest and extreme gradient lifting decision tree. The fishery production data was also collected from the operation data of 43 distant-water longline fishing vessels of China Aquatic Group in the western and Central Pacific (0°-30°S; 110°E-170°W) from 2008 to 2019, including catch information, such as amount, job date, as well as the job latitude and longitude. A comparison was performed on the fishery data, including the concentration of chlorophyll, eddy kinetic energy, sea surface height anomalies, temperature and salinity of the 0-500 m mixed water layer. A total of 36 variable combinations were used as the original data set, including the Southern Oscillation Index (SOI), the Arctic Oscillation Index (AOI), the Pacific Decadal Oscillation Index (PDOI), and North Pacific Gyre Oscillation Index (NPGOI). The original data set was divided into the training set and test set after the screening and normalization of the variance expansion factor, accounting for 80% and 20%, respectively. The training set was used to train eight models, including classification and regression, logistic regression, k-nearest neighbor, adaptive boosting, gradient boosting decision tree, xgboost, random forest, and XGBRF. The five-fold cross-validation was used for each model to determine the optimal parameters. Finally, the model was verified to superimpose the actual fishing ground of the test set. The experimental results showed that: 1) There was a significant correlation between the catch per unit fishing effort and various variable factors. There was also a great decrease in the degree of collinearity between the variables that were filtered by variance inflation factor. 2) The XGBRF hybrid ensemble model also significantly improved the performance of XGBoost and RF models. Specifically, the highest accuracy rate and Area Under Curve (AUC) were 75.39%, and 79.48%, respectively. The Receiver Operator Characteristic (ROC) curve of the XGBRF model was closer to the upper left, indicating the best performance of the forecasting model than before. 3) The sea surface temperature was the most important factor to dominate the distribution of yellowfin tuna fishing ground, accounting for 7.573%. The temperature of the 300 m water layer was equally important for the yellowfin tuna, which was 7.369%. In addition, the greater impact was also found in the salinity of the 50-meter water layer, the SOI, the concentration of chlorophyll, and the surface salinity. There was a relatively low influence of other large-scale climatic factors, except for the SOI. 4) There was only a small deviation between the fishing ground predicted by the XGBRF model and the actual fishing ground, indicating the high accuracy and reliability of the prediction. Overall, the XGBRF ensemble model performed the best on the fishing ground forecast of yellowfin tuna in the western and Central Pacific. The finding can also provide a strong reference for the fishing ground forecast.

    machine learning; models; western and central Pacific;yellowfin tuna; fishing ground forecast

    10.11975/j.issn.1002-6819.2022.15.036

    S931.3

    A

    1002-6819(2022)-15-0330-09

    張聰,周為峰,唐峰華,等. 基于機器學習的中西太平洋黃鰭金槍魚漁場預報模型[J]. 農(nóng)業(yè)工程學報,2022,38(15):330-338. doi:10.11975/j.issn.1002-6819.2022.15.036 http://www.tcsae.org

    Zhang Cong, Zhou Weifeng, Tang Fenghua, et al. Forecasting models for yellowfin tuna fishing ground in the central and western Pacific based on machine learning[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(15): 330-338. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2022.15.036 http://www.tcsae.org

    2022-04-29

    2022-07-19

    國家重點研發(fā)計劃項目(2019YFD0901405);農(nóng)業(yè)農(nóng)村部南海漁業(yè)資源開發(fā)利用重點實驗室開放基金項目(LOF 2022-05);中央級公益性科研院所基本科研業(yè)務(wù)費項目(2019T09);核電廠冷源安全保障相關(guān)技術(shù)開發(fā)項目(21FW018)

    張聰,研究方向為機器學習與海洋漁場預報。Emai:zhangcong20190404@163.com

    周為峰,博士,副研究員,研究方向為漁業(yè)遙感與漁場海洋學。Email:zhwfzhwf@163.com

    猜你喜歡
    黃鰭漁場金槍魚
    黃鰭鯛選育新品系“中珠1 號”養(yǎng)殖試驗
    黃鰭金槍魚深水網(wǎng)箱養(yǎng)殖技術(shù)初探
    不同規(guī)格黃鰭鯛的形態(tài)性狀與體質(zhì)量的相關(guān)性研究
    金槍魚淹死事件
    可以加速的金槍魚機器人
    軍事文摘(2021年22期)2022-01-18 06:22:00
    《關(guān)于加強遠洋魷釣漁船作業(yè)管理的通知》公布
    金槍魚與海豚
    黃鰭鯛苗種淡化技術(shù)研究
    埃及超大規(guī)模養(yǎng)魚場
    連續(xù)三年保持50%增長,入駐通威、新希望養(yǎng)殖基地,這家漁機企業(yè)將打造未來無人漁場
    国产亚洲欧美精品永久| 欧美日韩av久久| 色老头精品视频在线观看| 亚洲精品国产区一区二| 91麻豆精品激情在线观看国产 | 久久精品国产亚洲av香蕉五月 | 最新的欧美精品一区二区| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 国产老妇伦熟女老妇高清| 好男人电影高清在线观看| 欧美午夜高清在线| 色精品久久人妻99蜜桃| √禁漫天堂资源中文www| 精品国产国语对白av| 超碰97精品在线观看| 欧美黑人精品巨大| 精品国产乱码久久久久久男人| 69精品国产乱码久久久| 免费在线观看日本一区| 老汉色∧v一级毛片| 男女无遮挡免费网站观看| 成在线人永久免费视频| 国产国语露脸激情在线看| 熟女少妇亚洲综合色aaa.| 丰满少妇做爰视频| 麻豆av在线久日| 亚洲一卡2卡3卡4卡5卡精品中文| 久久99热这里只频精品6学生| 黄色成人免费大全| 国产精品麻豆人妻色哟哟久久| 国产三级黄色录像| 久久精品亚洲av国产电影网| svipshipincom国产片| 亚洲自偷自拍图片 自拍| 中文字幕制服av| a级片在线免费高清观看视频| 天堂中文最新版在线下载| 热99re8久久精品国产| 9色porny在线观看| 亚洲中文字幕日韩| 高清视频免费观看一区二区| 青草久久国产| 丰满人妻熟妇乱又伦精品不卡| 在线十欧美十亚洲十日本专区| 日韩欧美国产一区二区入口| 亚洲九九香蕉| 久久久久视频综合| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久久国内视频| 搡老乐熟女国产| 亚洲av片天天在线观看| 午夜免费成人在线视频| 俄罗斯特黄特色一大片| 亚洲精品乱久久久久久| 一本久久精品| 国产成人系列免费观看| 久久中文字幕一级| 老司机午夜福利在线观看视频 | 国产精品二区激情视频| 久久精品成人免费网站| √禁漫天堂资源中文www| 在线 av 中文字幕| 久热这里只有精品99| 精品免费久久久久久久清纯 | 女人爽到高潮嗷嗷叫在线视频| 亚洲全国av大片| 国产在线视频一区二区| 女警被强在线播放| 免费一级毛片在线播放高清视频 | 日本一区二区免费在线视频| 一级毛片精品| 亚洲伊人色综图| 人人妻人人澡人人爽人人夜夜| 美女主播在线视频| 99国产精品99久久久久| 午夜视频精品福利| 国产成人精品在线电影| 十分钟在线观看高清视频www| 王馨瑶露胸无遮挡在线观看| 丁香六月天网| 国产成人精品久久二区二区91| 十分钟在线观看高清视频www| 91九色精品人成在线观看| 国产精品影院久久| 99国产综合亚洲精品| 三级毛片av免费| 国产免费视频播放在线视频| 男女午夜视频在线观看| 亚洲国产精品一区二区三区在线| 十分钟在线观看高清视频www| 一区二区三区乱码不卡18| 久久精品国产a三级三级三级| 国产亚洲欧美在线一区二区| 首页视频小说图片口味搜索| 老司机在亚洲福利影院| 一本一本久久a久久精品综合妖精| 亚洲av国产av综合av卡| 脱女人内裤的视频| 新久久久久国产一级毛片| 人人妻人人澡人人爽人人夜夜| 三级毛片av免费| 色综合婷婷激情| 99国产极品粉嫩在线观看| 国产一区二区三区在线臀色熟女 | 99riav亚洲国产免费| 国产精品一区二区在线不卡| 少妇粗大呻吟视频| 精品久久久久久久毛片微露脸| 操出白浆在线播放| 亚洲自偷自拍图片 自拍| 丁香欧美五月| a在线观看视频网站| 搡老乐熟女国产| 精品国产超薄肉色丝袜足j| 高清视频免费观看一区二区| 国产精品秋霞免费鲁丝片| 少妇猛男粗大的猛烈进出视频| 日韩免费av在线播放| 亚洲精品国产精品久久久不卡| 69精品国产乱码久久久| 超碰成人久久| 日韩欧美一区二区三区在线观看 | 国产欧美日韩综合在线一区二区| 精品卡一卡二卡四卡免费| 亚洲av成人不卡在线观看播放网| 亚洲精品中文字幕一二三四区 | www.999成人在线观看| 在线观看免费日韩欧美大片| 9色porny在线观看| 久久久久久久国产电影| 日本a在线网址| 亚洲avbb在线观看| 亚洲精品在线美女| 黑丝袜美女国产一区| av片东京热男人的天堂| 国产精品1区2区在线观看. | 在线天堂中文资源库| 国产精品美女特级片免费视频播放器 | av天堂在线播放| 少妇裸体淫交视频免费看高清 | 女人爽到高潮嗷嗷叫在线视频| 好男人电影高清在线观看| 99re在线观看精品视频| 啪啪无遮挡十八禁网站| 夫妻午夜视频| 日本黄色视频三级网站网址 | 成人手机av| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区三区视频了| 亚洲精品美女久久av网站| 波多野结衣一区麻豆| 丁香六月欧美| 国产精品久久久av美女十八| 久久精品亚洲av国产电影网| 亚洲午夜理论影院| 韩国精品一区二区三区| 天堂动漫精品| a级毛片在线看网站| 在线观看免费高清a一片| 麻豆乱淫一区二区| 精品国产一区二区三区四区第35| 国产在线视频一区二区| 免费观看a级毛片全部| 黑人欧美特级aaaaaa片| 水蜜桃什么品种好| 大型黄色视频在线免费观看| 一进一出好大好爽视频| 色94色欧美一区二区| 免费人妻精品一区二区三区视频| 午夜精品久久久久久毛片777| 国产成人免费观看mmmm| 久久国产精品人妻蜜桃| 中文字幕精品免费在线观看视频| 亚洲精品中文字幕在线视频| 国产不卡一卡二| 90打野战视频偷拍视频| 久久精品成人免费网站| 日韩中文字幕视频在线看片| 国产免费视频播放在线视频| 建设人人有责人人尽责人人享有的| 啦啦啦免费观看视频1| 视频区欧美日本亚洲| 国产国语露脸激情在线看| 99九九在线精品视频| 精品亚洲成国产av| 国产精品自产拍在线观看55亚洲 | 亚洲免费av在线视频| 老熟妇乱子伦视频在线观看| 性高湖久久久久久久久免费观看| 多毛熟女@视频| 国产亚洲一区二区精品| 日韩视频一区二区在线观看| 在线播放国产精品三级| 天天躁夜夜躁狠狠躁躁| av超薄肉色丝袜交足视频| 我的亚洲天堂| 丁香欧美五月| 操出白浆在线播放| 又紧又爽又黄一区二区| 欧美日本中文国产一区发布| 亚洲综合色网址| 在线观看www视频免费| 国产免费现黄频在线看| 久久影院123| 亚洲国产成人一精品久久久| 最近最新中文字幕大全免费视频| 国产99久久九九免费精品| 国产亚洲精品第一综合不卡| 男人操女人黄网站| 色老头精品视频在线观看| 亚洲久久久国产精品| 亚洲免费av在线视频| 黄色成人免费大全| 亚洲性夜色夜夜综合| 国产男女内射视频| 国产在线一区二区三区精| 交换朋友夫妻互换小说| 99九九在线精品视频| 国产精品久久久人人做人人爽| 成年人黄色毛片网站| av电影中文网址| 国产免费av片在线观看野外av| 69精品国产乱码久久久| 国产精品免费视频内射| 黑丝袜美女国产一区| 高清在线国产一区| 日韩成人在线观看一区二区三区| 这个男人来自地球电影免费观看| 夜夜骑夜夜射夜夜干| 午夜福利乱码中文字幕| 精品国产一区二区久久| 一区二区三区激情视频| 大型av网站在线播放| 天堂8中文在线网| 新久久久久国产一级毛片| 这个男人来自地球电影免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 91九色精品人成在线观看| av视频免费观看在线观看| 精品国产乱码久久久久久小说| 成年人午夜在线观看视频| 国产一区二区三区综合在线观看| 男女高潮啪啪啪动态图| 波多野结衣av一区二区av| 又大又爽又粗| 超碰成人久久| 午夜福利免费观看在线| 亚洲情色 制服丝袜| 极品少妇高潮喷水抽搐| 麻豆av在线久日| 亚洲午夜理论影院| 国产在线视频一区二区| 中文字幕最新亚洲高清| 99re在线观看精品视频| 免费不卡黄色视频| 欧美日韩亚洲高清精品| 免费看十八禁软件| 欧美国产精品va在线观看不卡| 日本精品一区二区三区蜜桃| 亚洲第一欧美日韩一区二区三区 | 国产成人精品久久二区二区91| 啦啦啦在线免费观看视频4| 美女午夜性视频免费| 免费观看a级毛片全部| 国产精品国产高清国产av | 亚洲色图 男人天堂 中文字幕| 50天的宝宝边吃奶边哭怎么回事| 大型黄色视频在线免费观看| 久久精品亚洲av国产电影网| 精品国产国语对白av| 国产精品影院久久| 99香蕉大伊视频| 国产精品久久久av美女十八| av超薄肉色丝袜交足视频| 大片免费播放器 马上看| 国产欧美日韩一区二区三| 欧美日本中文国产一区发布| 精品久久蜜臀av无| 日韩免费av在线播放| 一级,二级,三级黄色视频| 搡老熟女国产l中国老女人| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| 伊人久久大香线蕉亚洲五| 中亚洲国语对白在线视频| 黑人操中国人逼视频| 一本一本久久a久久精品综合妖精| 精品欧美一区二区三区在线| 乱人伦中国视频| 最近最新免费中文字幕在线| 91九色精品人成在线观看| 亚洲成av片中文字幕在线观看| 精品欧美一区二区三区在线| 精品亚洲成a人片在线观看| 亚洲色图 男人天堂 中文字幕| 国产精品久久久久久精品电影小说| 亚洲av日韩精品久久久久久密| 国产人伦9x9x在线观看| 日韩欧美一区视频在线观看| 国产精品秋霞免费鲁丝片| 久久久精品区二区三区| 免费看十八禁软件| 国产精品99久久99久久久不卡| 俄罗斯特黄特色一大片| 丝袜美足系列| 19禁男女啪啪无遮挡网站| 伦理电影免费视频| 男女无遮挡免费网站观看| 亚洲伊人久久精品综合| 欧美 日韩 精品 国产| 母亲3免费完整高清在线观看| 午夜91福利影院| 久久久久精品人妻al黑| 一级毛片电影观看| 久久亚洲精品不卡| 国产单亲对白刺激| 国产高清视频在线播放一区| 男人操女人黄网站| 999久久久精品免费观看国产| 交换朋友夫妻互换小说| www日本在线高清视频| 久久久国产精品麻豆| 国产在线精品亚洲第一网站| 久久久精品免费免费高清| 中文字幕色久视频| 精品人妻1区二区| tocl精华| 不卡一级毛片| 国产色视频综合| 曰老女人黄片| 国产成人精品无人区| 热re99久久国产66热| 人妻一区二区av| 国产激情久久老熟女| 免费黄频网站在线观看国产| 国产一区二区在线观看av| 女人被躁到高潮嗷嗷叫费观| 久久性视频一级片| 久久久精品区二区三区| 在线观看免费日韩欧美大片| 天堂俺去俺来也www色官网| 老鸭窝网址在线观看| 久久精品亚洲熟妇少妇任你| 亚洲av第一区精品v没综合| 精品国产国语对白av| 欧美激情极品国产一区二区三区| 91老司机精品| 91精品三级在线观看| 久久影院123| 一区二区日韩欧美中文字幕| 18禁裸乳无遮挡动漫免费视频| 999久久久国产精品视频| 欧美乱妇无乱码| 亚洲欧美激情在线| 一边摸一边抽搐一进一小说 | 啦啦啦免费观看视频1| 两人在一起打扑克的视频| 亚洲专区中文字幕在线| 人人妻人人澡人人看| 超碰97精品在线观看| 日韩大码丰满熟妇| 亚洲欧美激情在线| 大香蕉久久成人网| 久久久久久免费高清国产稀缺| 免费看十八禁软件| 国产日韩欧美视频二区| 桃红色精品国产亚洲av| 欧美成人午夜精品| 国产精品偷伦视频观看了| 一级黄色大片毛片| 亚洲精华国产精华精| 欧美日韩av久久| 69av精品久久久久久 | 亚洲人成电影观看| 激情视频va一区二区三区| 又大又爽又粗| 激情视频va一区二区三区| 国产成人影院久久av| 97人妻天天添夜夜摸| 交换朋友夫妻互换小说| 亚洲专区国产一区二区| 国产成人影院久久av| 在线观看免费视频日本深夜| 国产黄色免费在线视频| 一区二区三区激情视频| tocl精华| 日韩免费av在线播放| 精品国产乱码久久久久久小说| 免费在线观看影片大全网站| 国产精品美女特级片免费视频播放器 | 老司机影院毛片| 一级毛片精品| 国产成人精品在线电影| 天天躁日日躁夜夜躁夜夜| 国产一区二区三区在线臀色熟女 | 精品第一国产精品| 国产精品国产av在线观看| 丰满人妻熟妇乱又伦精品不卡| 久久久精品免费免费高清| 夫妻午夜视频| 欧美激情 高清一区二区三区| 搡老乐熟女国产| 18禁国产床啪视频网站| 50天的宝宝边吃奶边哭怎么回事| 黄色视频不卡| 王馨瑶露胸无遮挡在线观看| 亚洲精品久久成人aⅴ小说| 熟女少妇亚洲综合色aaa.| 日韩精品免费视频一区二区三区| 丰满人妻熟妇乱又伦精品不卡| av视频免费观看在线观看| 黄色丝袜av网址大全| kizo精华| 看免费av毛片| 91成年电影在线观看| 成人精品一区二区免费| 老司机影院毛片| 女人爽到高潮嗷嗷叫在线视频| 精品国产乱码久久久久久男人| 一本大道久久a久久精品| 亚洲七黄色美女视频| 在线观看舔阴道视频| 精品国产一区二区三区久久久樱花| 丝袜美足系列| 久久人妻熟女aⅴ| 欧美国产精品va在线观看不卡| 亚洲七黄色美女视频| 亚洲自偷自拍图片 自拍| 狠狠狠狠99中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲九九香蕉| 曰老女人黄片| 一区二区三区精品91| 欧美亚洲日本最大视频资源| 一区二区av电影网| 国产一区有黄有色的免费视频| 亚洲第一欧美日韩一区二区三区 | 国产精品免费大片| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 久久精品亚洲精品国产色婷小说| videosex国产| 亚洲精品国产色婷婷电影| 久久久国产欧美日韩av| 操出白浆在线播放| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 夜夜骑夜夜射夜夜干| 成人亚洲精品一区在线观看| 精品国产乱子伦一区二区三区| 欧美黑人精品巨大| 十分钟在线观看高清视频www| 国产高清视频在线播放一区| 亚洲五月婷婷丁香| 久久国产精品影院| 高清黄色对白视频在线免费看| 一区二区三区国产精品乱码| 欧美变态另类bdsm刘玥| 欧美成人免费av一区二区三区 | 中文字幕av电影在线播放| 免费看a级黄色片| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲一级av第二区| 国产精品 国内视频| 欧美性长视频在线观看| www.999成人在线观看| 可以免费在线观看a视频的电影网站| 少妇粗大呻吟视频| 国产色视频综合| 亚洲欧美日韩另类电影网站| 国内毛片毛片毛片毛片毛片| 亚洲精品久久午夜乱码| 国产亚洲精品第一综合不卡| 美女主播在线视频| 黄色成人免费大全| 国产精品.久久久| 黑人欧美特级aaaaaa片| 人妻 亚洲 视频| 美女扒开内裤让男人捅视频| 久久人人97超碰香蕉20202| 18在线观看网站| 国产精品熟女久久久久浪| 啦啦啦 在线观看视频| 黑人猛操日本美女一级片| 免费少妇av软件| 国产国语露脸激情在线看| 国产亚洲一区二区精品| 亚洲国产欧美一区二区综合| 夜夜骑夜夜射夜夜干| 色94色欧美一区二区| 91老司机精品| 欧美午夜高清在线| 久久久国产一区二区| 午夜免费成人在线视频| 香蕉国产在线看| 在线观看免费视频日本深夜| 久久久久久久精品吃奶| 日日摸夜夜添夜夜添小说| 深夜精品福利| 亚洲欧洲日产国产| 色在线成人网| 黄色片一级片一级黄色片| 飞空精品影院首页| 欧美乱妇无乱码| 精品国产超薄肉色丝袜足j| 黄片大片在线免费观看| 亚洲欧美精品综合一区二区三区| 久9热在线精品视频| 日本五十路高清| 人妻久久中文字幕网| 国产精品国产av在线观看| xxxhd国产人妻xxx| 免费女性裸体啪啪无遮挡网站| 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 精品亚洲成a人片在线观看| svipshipincom国产片| 日韩三级视频一区二区三区| 不卡av一区二区三区| 一区二区三区精品91| 乱人伦中国视频| 国产激情久久老熟女| 精品人妻熟女毛片av久久网站| 黑人操中国人逼视频| 国产精品 国内视频| 国产高清视频在线播放一区| 777米奇影视久久| 欧美黑人欧美精品刺激| 国产在视频线精品| 久久精品91无色码中文字幕| 欧美+亚洲+日韩+国产| 中文字幕色久视频| 欧美黑人欧美精品刺激| 亚洲专区字幕在线| 一区二区三区精品91| 亚洲少妇的诱惑av| 欧美中文综合在线视频| www.999成人在线观看| 又黄又粗又硬又大视频| 久久性视频一级片| 无遮挡黄片免费观看| 久久人妻av系列| 成在线人永久免费视频| 国产视频一区二区在线看| 免费在线观看黄色视频的| 久久狼人影院| 亚洲色图综合在线观看| 久久久久久久精品吃奶| 国产高清视频在线播放一区| 久久久国产一区二区| 欧美日本中文国产一区发布| www日本在线高清视频| 欧美日韩视频精品一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久影院123| 国产精品av久久久久免费| 麻豆成人av在线观看| 精品久久久久久电影网| 91麻豆精品激情在线观看国产 | 少妇的丰满在线观看| 国产精品二区激情视频| 热99国产精品久久久久久7| 亚洲五月婷婷丁香| tube8黄色片| 两个人看的免费小视频| 免费女性裸体啪啪无遮挡网站| 久久久久久免费高清国产稀缺| 久久精品国产99精品国产亚洲性色 | 俄罗斯特黄特色一大片| 亚洲三区欧美一区| 国产区一区二久久| 欧美黑人欧美精品刺激| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣一区麻豆| 精品国产一区二区三区四区第35| 日日爽夜夜爽网站| a级毛片在线看网站| 国产精品影院久久| 国产人伦9x9x在线观看| 国内毛片毛片毛片毛片毛片| 国产无遮挡羞羞视频在线观看| 狂野欧美激情性xxxx| 欧美黑人欧美精品刺激| 热99国产精品久久久久久7| 别揉我奶头~嗯~啊~动态视频| 精品视频人人做人人爽| 精品少妇内射三级| 亚洲精品自拍成人| 极品少妇高潮喷水抽搐| 免费av中文字幕在线| 亚洲av电影在线进入| 黄网站色视频无遮挡免费观看| 国产单亲对白刺激| 国产欧美亚洲国产| 777久久人妻少妇嫩草av网站| 菩萨蛮人人尽说江南好唐韦庄| 欧美激情 高清一区二区三区| 欧美 亚洲 国产 日韩一| 99香蕉大伊视频| 欧美变态另类bdsm刘玥| 一本久久精品| 久9热在线精品视频| 一区福利在线观看| 一区二区三区国产精品乱码| 国产成人影院久久av| 亚洲欧美色中文字幕在线| 日日摸夜夜添夜夜添小说| 亚洲av国产av综合av卡| 天天躁狠狠躁夜夜躁狠狠躁| 国产av国产精品国产| 母亲3免费完整高清在线观看| www.自偷自拍.com| 91大片在线观看| 免费观看av网站的网址|