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

    BP和RBF神經(jīng)網(wǎng)絡應用于海表溫鹽短期預測效果對比

    2022-06-06 01:36:56李亞蒙丁軍航孫寶楠官
    海洋科學進展 2022年2期
    關鍵詞:海表鹽度神經(jīng)網(wǎng)絡

    李亞蒙丁軍航孫寶楠官 晟*

    (1.青島大學 自動化學院,山東 青島 266071;2.自然資源部 第一海洋研究所,山東 青島 266061;3.山東省工業(yè)控制技術重點實驗室,山東 青島 266071;4.自然資源部 海洋環(huán)境科學與數(shù)值模擬重點實驗室,山東 青島 266061;5.山東省海洋環(huán)境科學與數(shù)值模擬重點實驗室,山東 青島 266061;6.青島海洋科學與技術試點國家實驗室 區(qū)域海洋動力學與數(shù)值模擬功能實驗室,山東 青島 266237)

    20世紀以來,在全球氣候變暖趨勢下,極端氣候及災害事件頻發(fā),海洋環(huán)境也變得愈加復雜。要準確掌握海洋環(huán)境變化的規(guī)律,進而增強人類認識海洋、開發(fā)海洋的能力,需要具備比以往更準確、更及時和更全面的觀測手段。海洋觀測浮標可利用其搭載的各類傳感器和設備,測量包括溫度和鹽度等海洋及大氣環(huán)境要素,并實時將觀測數(shù)據(jù)回傳岸基臺站。海洋觀測浮標具有長時間連續(xù)、定點觀測能力,搭載傳感器種類、數(shù)量較多,其觀測方式與海洋環(huán)境變化特征相吻合,是海洋立體觀測網(wǎng)絡中重要的節(jié)點形式。對觀測浮標獲取數(shù)據(jù)的有效利用與深度分析,可為全球氣候變化趨勢研究、海洋災害預警和人類海上作業(yè)等領域提供有效的數(shù)據(jù)支撐[1-3]。

    海洋觀測浮標獲取的多參數(shù)數(shù)據(jù)可以幫助人們了解海洋環(huán)境歷經(jīng)的狀態(tài),但對于人類海上活動而言,更需要從已有數(shù)據(jù)中預測未來將發(fā)生的變化,由此產(chǎn)生了人們對海洋環(huán)境要素進行預測的需求。相對于長期預測而言,海洋環(huán)境要素的時間序列短期預測面臨更多偶然因素干擾,實現(xiàn)精確預測難度更大,但短期預測的結果對海上工程作業(yè)在安全保障、方案選擇和成本控制等方面具有很高的應用價值。

    對于絕大多數(shù)由實際觀測資料構成的動力系統(tǒng),往往表現(xiàn)出非線性、非平穩(wěn)性雙重特征。即系統(tǒng)各要素非獨立,存在耦合作用,同時觀測資料的時間序列又表現(xiàn)出具有趨勢或周期性變化特性。對于這些動力系統(tǒng),通常利用統(tǒng)計預測理論和非線性預測理論開展預測,由于其對信號平穩(wěn)性的假定,有可能導致預測水平低下[4]。也有研究[5]通過對數(shù)變換或差分變換將相應的過程平穩(wěn)化,對能夠尋找近似平穩(wěn)信號分量的非平穩(wěn)過程取得了不錯的預測效果。越來越多的研究[6-12]將神經(jīng)網(wǎng)絡預報法應用于氣象和水文等非線性預測領域。常見的神經(jīng)網(wǎng)絡預報法包括反向傳播(Back Propagation,BP)神經(jīng)網(wǎng)絡、徑向基函數(shù)(Radial Basis Function,RBF)神經(jīng)網(wǎng)絡和差分自回歸移動平均模型(Autoregressive Integrated Moving Average model,ARIMA)等。相比之下,海洋動力學預報法需要大型計算機對流體力學微分方程進行長時間積分,所需計算成本和對研究者建模水平的要求都較高。而神經(jīng)網(wǎng)絡方法逼近能力強、收斂速度快,具有信息分布式存儲和并行協(xié)同處理的特點[13],實現(xiàn)相對簡單,在近年海水溫、鹽預報應用研究中更為普遍。

    基于BP神經(jīng)網(wǎng)絡和隨機森林算法的評估模型[14]可以用來精準估算近岸海域的海表鹽度。高國棟等[15]基于RBF神經(jīng)網(wǎng)絡模型建立了海水鹽度的預測方法,并證明了該模型相較于BP神經(jīng)網(wǎng)絡模型收斂速度更快、精度更高。建立基于機器學習的反演模型并利用遙感技術進行大尺度動態(tài)模擬[16],或利用衛(wèi)星輔助數(shù)據(jù)建立BP神經(jīng)網(wǎng)絡預測模型[17]均能較為精確地預測海表鹽度。王穎超等[17]利用BP神經(jīng)網(wǎng)絡方法預測的海表鹽度以及土壤水分和海洋鹽度(Soil Moisture and Ocean Salinity,SMOS)衛(wèi)星提供的3個粗糙度模型給出了鹽度產(chǎn)品與ARGO 浮標實測鹽度的均方根誤差(Root Mean Square Error,RMSE),分別為0.847、2.041、2.028和2.081。基于空間遙感數(shù)據(jù),使用神經(jīng)網(wǎng)絡算法預測[18-19]或者將深度學習神經(jīng)網(wǎng)絡與數(shù)值估值器相結合[20]均能有效預測海表溫度。YANG 等[18]使用長短期記憶(Long Short-Term Memory,LSTM)神經(jīng)網(wǎng)絡預測了渤海1 d的海表溫度數(shù)據(jù),支持向量機、支持向量回歸模型、FC-LSTM 等算法的預測結果與遙感實測值的RMSE為0.076~0.399,預測精度為98.49%~99.58%,證明了所提出的神經(jīng)網(wǎng)絡預測方法可行,但也發(fā)現(xiàn)數(shù)據(jù)集空間分辨率越高,精確預測難度越大。

    以往研究結果證明了多種神經(jīng)網(wǎng)絡預測技術能夠有效應用于海表溫、鹽時間序列數(shù)據(jù)的預測工作。然而在這些研究中,未討論不同訓練時長和預測時長等因素對預測精度的影響,也沒有討論算法的多海域適用性等問題?;谏鲜銮闆r,本文將利用BP和RBF兩種神經(jīng)網(wǎng)絡方法開展研究,以期獲得一種具有穩(wěn)定預報準確率,并具有較廣泛適用性的海表溫、鹽短期預報方法。

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

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

    本研究數(shù)據(jù)主要來自美國國家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)組織實施的OceanSITES觀測系統(tǒng)中PAPA 站(145°W,50°N)、Stratus站(85°W,20°S)、KEO 站(145°E,32°N)、WHOTS站(160°W,25°N)和JKEO 站(146°E,38°N)等多個浮標站位于2009—2012年的監(jiān)測數(shù)據(jù)。這些數(shù)據(jù)包括0~300 m 深度的溫度、0~200 m 深度的鹽度和密度、表層海流、緯向和經(jīng)向海流的聲學多普勒剖面等,以及風速、氣溫、相對濕度和大氣壓等氣象觀測數(shù)據(jù)。本文采用Mat Lab軟件及其機器學習工具箱進行短期預測方法研究。

    1.2 實驗方法

    1.2.1 BP神經(jīng)網(wǎng)絡預測模型的構建

    BP神經(jīng)網(wǎng)絡是一種按照誤差逆向傳播算法訓練的多層前饋神經(jīng)網(wǎng)絡,是目前應用最為廣泛的神經(jīng)網(wǎng)絡。BP網(wǎng)絡按結構分為輸入層、隱藏層和輸出層,如圖1所示。

    圖1 BP神經(jīng)網(wǎng)絡結構Fig.1 Structure diagram of BP neural network

    對于隱藏層,利用激活函數(shù)處理輸入可以得到隱藏層輸出,公式為:

    式中,f(net j)為隱藏層激活函數(shù),其中

    對于輸出層,利用激活函數(shù)處理隱藏層輸出可以得到神經(jīng)元輸出,公式為:

    式中,f(net k)為輸出層激活函數(shù),其中

    BP 神經(jīng)網(wǎng)絡的優(yōu)點是擁有很好的非線性映射能力和靈活的網(wǎng)絡結構。網(wǎng)絡各層的神經(jīng)元個數(shù)可根據(jù)具體情況設定,并且其性能也會隨著結構的差異而有所不同。但是BP神經(jīng)網(wǎng)絡也存在學習速度慢、因容易陷入局部極小值而無法得到全局最優(yōu)值等主要缺陷。建立BP神經(jīng)網(wǎng)絡模型的步驟主要有5步[21]。

    步驟一:讀取研究對象數(shù)據(jù)集,如海表溫度數(shù)據(jù)或海表鹽度數(shù)據(jù),對缺損值和突兀點進行預處理。

    步驟二:劃分訓練集和測試集,對海表溫度和海表鹽度數(shù)據(jù)進行歸一化處理。一般情況下,訓練集的數(shù)據(jù)量會影響B(tài)P神經(jīng)網(wǎng)絡的預測效果。因此,以PAPA 站2009—2012年每年4月海表溫度數(shù)據(jù)為樣本,在預測時長固定為5 d的情況下,對比5 d、10 d、15 d、20 d和25 d不同時長訓練海表溫度數(shù)據(jù)集預測結果的均方誤差,結果如表1所示。由表1可知,選取20 d的觀測海表溫度數(shù)據(jù)作為訓練集,預測結果相對最優(yōu),訓練時長也在可接受范圍內。

    表1 2009—2012年每年4月BP神經(jīng)網(wǎng)絡不同訓練時長對應海表溫度預測結果的均方誤差Table 1 MSE of SST prediction results of different training time by BP neural network in April 2009—2012

    步驟三:確定輸入層節(jié)點和輸出層節(jié)點數(shù)。實驗中將歷史海表溫度和海表鹽度數(shù)據(jù)作為輸入,因此,選取輸入層節(jié)點數(shù)為1,分別選擇下一時刻的海表溫度和海表鹽度作為輸出,因此,實驗中輸出層節(jié)點數(shù)為1。

    步驟四:確定網(wǎng)絡各層之間的傳遞函數(shù)和訓練算法。BP神經(jīng)網(wǎng)絡使用最多的激活函數(shù)為sigmoid函數(shù),其輸入值可取任意值,輸出值范圍為0~1。本文選擇隱藏層激活函數(shù)為tansig,其輸入值為任意值,輸出值為-1~1。選擇輸出層激活函數(shù)為purelin,其輸入和輸出值均為任意值。訓練算法為Levenberg-Marquardt[trainlm],該算法作為系統(tǒng)的默認算法,其優(yōu)勢在于針對中等規(guī)模的BP神經(jīng)網(wǎng)絡,收斂速度最快。

    通過網(wǎng)絡測試,對比隱藏層節(jié)點數(shù)對應的訓練集均方誤差值,可以選擇出最優(yōu)的隱藏層節(jié)點個數(shù),并應用于海洋數(shù)據(jù)預測。同樣以2009—2012年每年4月海表溫度的預測為例,設置每月前20 d海表溫度數(shù)據(jù)為訓練集進行實驗,訓練集均方誤差結果如表2所示。由表2可知,最優(yōu)的隱藏層節(jié)點數(shù)為8。

    表2 BP神經(jīng)網(wǎng)絡隱藏層節(jié)點數(shù)對應海表溫度訓練集的均方誤差Table 2 The MSE of SST training set corresponding to the number of hidden layer nodes of BP neural network

    步驟五:分析BP 神經(jīng)網(wǎng)絡的預測值與實測值誤差。通過誤差分析,研究預測方法的適用性及策略改進。

    1.2.2 RBF神經(jīng)網(wǎng)絡預測模型構建

    徑向基函數(shù)RBF神經(jīng)網(wǎng)絡主要用于解決多變量插值的問題。其取值依賴于空間任意一點x到中心點c的距離,即。最常用的徑向基函數(shù)是高斯核函數(shù)[22],其形式為:

    式中,x c是核函數(shù)中心,σ是函數(shù)的寬度參數(shù)。

    RBF神經(jīng)網(wǎng)絡屬于前向神經(jīng)網(wǎng)絡類型,擁有3層前向網(wǎng)絡結構,包括輸入層、隱藏層和輸出層。輸入層由信號源結點構成;隱藏層神經(jīng)單元的激活函數(shù)是徑向基函數(shù),它是對應中心點徑向對稱且逐漸衰減的非線性函數(shù),BP 神經(jīng)網(wǎng)絡可以有多個隱藏層,但是RBF 只有一個隱藏層;輸出層對輸入模式做出響應。從輸入層到隱藏層的變換是非線性的,而從隱藏層到輸出層的變換是線性的[23]。RBF神經(jīng)網(wǎng)絡結構如圖2所示。

    圖2 RBF神經(jīng)網(wǎng)絡結構Fig.2 Structure diagram of RBF neural network

    RBF神經(jīng)網(wǎng)絡實質上是利用由徑向基函數(shù)構成的隱藏層,將輸入變量直接映射至隱藏空間,通過確定RBF 中心點來控制映射關系。通常徑向基函數(shù)對于空間中心點具備徑向對稱性,當神經(jīng)元輸入距離中心點較遠時,神經(jīng)元的激活程度較低,這屬于RBF 神經(jīng)網(wǎng)絡的局部特性[24]。隱藏層的作用是將輸入向量由低維度n映射成高維度h,使得向量具有線性可分的特性,通常利用高斯核函數(shù)作為激活函數(shù):

    式中,k=1,2,…,n;i=1,2,…,h。

    此時網(wǎng)絡的輸出是隱藏層輸出的線性加權和,權即為可調參數(shù),網(wǎng)絡輸出可表示為:

    式中,k=1,2,…,n,j=1,2,…,m;σ為最小二乘損失函數(shù),其中i=1,2,…,h;d j為理想輸出值。

    RBF神經(jīng)網(wǎng)絡的設計主要需要確認徑向基函數(shù)中心、方差和隱藏層到輸出層的權值三個參數(shù)。對于高斯核函數(shù)的徑向基網(wǎng)絡,首先利用k-means聚類得到h個中心,然后利用(式中,i=1,2,…,h,cmax表示選擇的中心點之間最大距離)計算徑向基函數(shù)方差,最后利用最小二乘法計算權值w,計算式為ω=

    RBF模型建立步驟基本與BP神經(jīng)網(wǎng)絡類似。同樣考慮到訓練集的數(shù)據(jù)量對RBF 神經(jīng)網(wǎng)絡預測效果的影響,以2009—2012年間每年4月海表溫度數(shù)據(jù)預測為例,在預測時長固定5 d情況下,設置不同訓練時長進行對比,預測結果的均方誤差如表3所示。由表3可知,選取20 d的海表溫度作為訓練集,預測結果相對最優(yōu)。

    表3 RBF神經(jīng)網(wǎng)絡訓練時間對應海表溫度預測結果均方誤差Table 3 Mean square error of SST prediction corresponding to different training days by RBF neural network

    1.2.3 評價標準

    本文采取均方誤差(MSE)和平均絕對誤差(Mean Absolute Error,MAE)兩個指標來對檢測結果進行分析。均方誤差是指參數(shù)估計值與參數(shù)真值之差平方的期望值,可以評價數(shù)據(jù)的變化程度,MSE 值越小,說明預測模型描述的實驗數(shù)據(jù)具有更好的精確度;平均絕對誤差能更好地反映預測值誤差的實際情況。2個指標的計算公式為:

    式中,N為模型樣本數(shù),f i為海表溫度或海表鹽度的實測值,y i為海表溫度或海表鹽度的預測值。

    相關系數(shù)是一種研究變量之間線性相關程度的量。相關系數(shù)以變量與各自平均的差值為基礎,按積差的方式進行計算,其表達式為:處于0~1之間越接近1,代表預測值與實測值之間相關程度越高,即預測性能越好;r

    越接近0則代表預測性能越差。

    2 結果與分析

    海洋表面的海水溫度、鹽度等物理參數(shù)受到波浪、海流、風及海氣交換等諸多不穩(wěn)定因素的影響,有時變化劇烈、隨機性較強,利用傳統(tǒng)擬合、客觀分析等方法難以進行準確的短期預測。

    本文首先選擇PAPA 站在2009年的1月、4月、7月和10月的海表溫度和海表鹽度數(shù)據(jù)作為研究數(shù)據(jù)集。1月、4月、7月和10月是常用的北半球冬季、春季、夏季、秋季四個季節(jié)典型環(huán)境特點代表月份。PAPA站浮標搭載傳感器對海表溫度及海表鹽度采樣周期均為1 h。根據(jù)前文所述,預測模型所用海表溫度及海表鹽度實測數(shù)據(jù)的訓練集時間長度為20 d,短期預測時間長度為5 d,共25 d、合計600 h的海表溫度及海表鹽度數(shù)據(jù)變化曲線如圖3和圖4所示。

    從圖3和圖4可見,海表溫度數(shù)據(jù)表現(xiàn)出了明顯的季節(jié)變化特征與日周期的變化特點,在這些趨勢性和固定的周期變化基礎上,疊加了非線性、隨機的擾動。其中1月和10月海表溫度呈明顯下降趨勢,分別從5.9 ℃下降到5.2 ℃、從11.0 ℃下降到9.0 ℃;4月和7月海表溫度呈明顯上升趨勢,分別從7.0 ℃上升到9.0 ℃、從12.0 ℃上升到13.2 ℃。而海表鹽度數(shù)據(jù)的季節(jié)性趨勢相對不明顯,整體變化范圍較小,其中1月和4月大部分時間的海表鹽度實測數(shù)據(jù)均超過32.50,而7月和10月大部分時間的海表鹽度實測數(shù)據(jù)卻均低于32.50。

    圖3 冬、春、夏、秋四個季節(jié)海表溫度隨時間的變化Fig.3 Changes of SST with time in winter,spring,summer and autumn

    圖4 冬、春、夏、秋四個季節(jié)海表鹽度隨時間的變化Fig.4 Changes of SSS with time in winter,spring,summer and autumn

    2.1 預測對比實驗

    分別采用BP和RBF神經(jīng)網(wǎng)絡模型,對PAPA 站4個代表月份的21至25日、共120 h的海表溫度與海表鹽度進行了短期預報。預報結果如圖5和圖6所示。訓練集均為各月前20 d的實測數(shù)據(jù)。BP神經(jīng)網(wǎng)絡隱藏層節(jié)點數(shù)為8,RBF神經(jīng)網(wǎng)絡隱藏層節(jié)點數(shù)為8。

    從圖5和圖6看出,2種神經(jīng)網(wǎng)絡所獲得的預測結果能夠反映實測數(shù)據(jù)的變化趨勢,對局部細節(jié)也有所反應。2種方法的預測效果隨著預測時長增加均出現(xiàn)相對實測數(shù)據(jù)的偏離,其中BP神經(jīng)網(wǎng)絡的預測結果偏離更為顯著。對2種方法預測結果的量化評價,主要從均方誤差(MSE)和平均絕對誤差(MAE)來衡量。

    圖5 4個季節(jié)代表月份21至25日海表溫度預測結果Fig.5 Prediction results of SST from 21 to 25 in representative months of four seasons

    圖6 4個季節(jié)代表月份21至25日海表鹽度預測結果Fig.6 Prediction results of SSS from 21 to 25 in representative months of four seasons

    在溫度短期預報方面(表4),2種方法的MAE均遠小于0.5 ℃。其中誤差超過《海洋調查規(guī)范》[25]規(guī)定的0.02 ℃標準的情況,RBF預測結果有1次,BP預測結果有2次。由于鹽度實測值(表4)較為穩(wěn)定,各項誤差指數(shù)均優(yōu)于溫度預報結果。其中誤差超過《海洋調查規(guī)范》[25]規(guī)定的0.02標準的情況,只有BP預測結果有1次。RBF對海溫和鹽度的預測最小MAE分別為0.99×10-2℃和9.07×10-4,BP對海溫和鹽度的預測最小MAE分別為1.50×10-2℃和0.12×10-2。

    表4 2種預測方法預測海表溫度和海表鹽度誤差分析Table 4 Error analysis of two methods for predicting SST and SSS

    利用式(8)計算RBF和BP兩種方法的預測結果與實測結果的相關系數(shù),結果如表5和表6所示。從相關系數(shù)指標看,RBF預測方法對海表溫度和海表鹽度的預測結果,在所有月份均優(yōu)于BP預測方法的結果。受篇幅所限,本文只展示了2種預測方法的溫度(圖7a)和鹽度(圖7b)單月預測結果的相關性趨勢,其他3個月情況沒有大的差異。在相關性趨勢圖中,RBF預測方法結果分布點也更集中于y=x直線,沒有明顯的偏離點。因此,就整體預測效果而言,RBF神經(jīng)網(wǎng)絡方法更精準。

    圖7 RBF和BP兩種神經(jīng)網(wǎng)絡方法1月海表溫度、4月海表鹽度預測值與實測值的相關性Fig.7 Correlation between the predicted and measured values of monthly SST or SSS by the two methods

    表5 2種預測方法的海表溫度預測值與實測值相關系數(shù)Table 5 Correlation coefficients between predicted and measured SST of the two methods

    表6 2種預測方法的海表鹽度預測值與實測值相關系數(shù)Table 6 Correlation coefficients between predicted and measured SSS of the two methods

    此外,對于不同季節(jié)的海表溫度和海表鹽度數(shù)據(jù),2種神經(jīng)網(wǎng)絡方法預測效果的準確性均達到海表溫度模式預報0.5 ℃的精度,海表溫度和海表鹽度的短期預報精度均接近《海洋調查規(guī)范》[25]規(guī)定的0.02 ℃和0.02標準的要求。預測結果的誤差指標證明,這2種神經(jīng)網(wǎng)絡預測方法具備針對不同季節(jié)的海洋環(huán)境的適應能力,可以預測整年的海洋溫鹽數(shù)據(jù)。

    5 d是典型的短期預報時長,但是在以20 d實測數(shù)據(jù)為訓練集的前提下,對RBF和BP兩種神經(jīng)網(wǎng)絡預測方法而言,對能獲取更高預測準確率的最佳預測周期可開展進一步研究。

    設置預測時長為1 d、2 d、3 d、4 d和5 d,利用2種神經(jīng)網(wǎng)絡模型對PAPA 站2009年各季節(jié)海表溫度和海表鹽度進行預測對比,訓練集仍選最佳訓練時長20 d。預測值各月MAE 的結果如圖8 和圖9所示。由圖8和圖9可知,隨著測試集預測時長的變化,利用RBF神經(jīng)網(wǎng)絡方法預測的海表溫度和海表鹽度的MAE始終低于BP神經(jīng)網(wǎng)絡。由于海表溫度的變化幅度大于海表鹽度,因此隨著預測時長的增加,海表溫度預測的誤差增加較為明顯,除個別情形(如RBF 對7 月海表溫度的預測)外,2 種預測模型的海表溫度預測最佳周期可選為1 d;由于海表鹽度的變化幅度較小,因此預測誤差隨預測時長的變化不太明顯,整體趨勢為隨著預測時長增加,預測誤差也對應增大,最佳預測周期可選2~3 d。在實際應用中,需綜合考慮預測成本與準確性的實際需求確定預測時長。

    圖8 RBF和BP兩種神經(jīng)網(wǎng)絡預測方法在4個代表月份以不同預測時長預測海表溫度結果的平均絕對誤差Fig.8 The MAE of SST predicted by two prediction methods with different prediction days in four representative months

    圖9 RBF和BP兩種神經(jīng)網(wǎng)絡預測方法在4個代表月份以不同預測時長預測海表鹽度結果的平均絕對誤差Fig.9 The MAE of SSS predicted by two prediction methods with different prediction days in four representative months

    2.2 多站點預測驗證

    為了驗證預測模型的適用性和普遍性,選擇PAPA 站以外的其他4個站點的浮標數(shù)據(jù),利用RBF神經(jīng)網(wǎng)絡模型對海表溫鹽數(shù)據(jù)進行預測驗證。這4 個站點分別為Stratus站(85°W,20°S)、KEO 站(145°E,32°N)、WHOTS站(160°W,25°N)和JKEO 站(146°E,38°N)。各站點海表溫度和海表鹽度數(shù)據(jù)采樣周期也為1 h,模型參數(shù)通過上述實驗方法選定,海表溫度預測時長為1 d,海表鹽度預測時間為3 d。其中,利用位于東太平洋距離較近的KEO 站和JKEO 站對2010年的海表溫度和海表鹽度數(shù)據(jù)進行重復性實驗,驗證RBF神經(jīng)網(wǎng)絡模型預測能力的可重復性;利用中太平洋相距較遠的WHOTS站和西太平洋Stratus站2012年的海表溫度和海表鹽度數(shù)據(jù)驗證預測模型的普適性。由于數(shù)據(jù)較多,本文只給出Stratus站(85°W,20°S)各季節(jié)的預測結果(圖10和圖11)。各站點的預測誤差及相關系數(shù)如表7~表10所示。

    表7 RBF神經(jīng)網(wǎng)絡海表溫度預測誤差Table 7 Prediction error of SST by RBF neural network

    表8 RBF神經(jīng)網(wǎng)絡海表鹽度預測誤差Table 8 Prediction error of SSS by RBF neural network

    表9 RBF神經(jīng)網(wǎng)絡預測海表溫度的預測值與實測值相關系數(shù)Table 9 Correlation coefficients between predicted and measured SST predicted by RBF neural network

    表10 RBF神經(jīng)網(wǎng)絡預測海表鹽度的預測值與實測值相關系數(shù)Table 10 Correlation coefficients between predicted and measured SSS predicted by RBF neural network

    圖10 Stratus站4個代表月份1 d海表溫度預測結果Fig.10 1-day SST prediction results of four representative months at Stratus Station

    圖11 Stratus站4個代表月份3 d海表鹽度預測結果Fig.11 3-day SSS prediction results of four representative months at Stratus Station

    由多站點預測實驗結果的誤差分析和相關系數(shù)比對可知,RBF神經(jīng)網(wǎng)絡模型對于海洋表面溫鹽數(shù)據(jù)的超短期預測具備普遍適用性,并能精確擬合出多個海域及其四季海表溫鹽數(shù)據(jù)的變化規(guī)律。這表明,盡管海表溫鹽等數(shù)據(jù)受諸多不穩(wěn)定因素影響而具備隨機性、非線性的特點,但使用RBF神經(jīng)網(wǎng)絡能夠進行有效、高精度的預測。相較于已往研究成果,應用本方法對獨立站位海表溫度數(shù)據(jù)進行1 d預測、對海表鹽度數(shù)據(jù)進行3 d預測,預測結果在均方誤差和平均絕對誤差指標方面,均優(yōu)于已有的海表溫度和海表鹽度遙感數(shù)據(jù)預測研究所取得的結果。

    3 結論

    基于海洋環(huán)境復雜多變,難以做到精準預測的現(xiàn)狀,本文利用BP和RBF神經(jīng)網(wǎng)絡模型對海表溫度和海表鹽度進行預測,在Mat Lab環(huán)境下搭建了RBF和BP兩種神經(jīng)網(wǎng)絡的預測模型,開展了海表溫度和海表鹽度的短期預測和結果比對研究。主要結論如下:

    首先,通過單站點四個季節(jié)代表月份的海表溫度和海表鹽度預測實驗,證明了BP和RBF 神經(jīng)網(wǎng)絡模型均能有效模擬海表溫度和海表鹽度數(shù)據(jù)的短期變化趨勢,對季節(jié)性變化均有較強的適應性,RBF 神經(jīng)網(wǎng)絡擁有更高的預測準確性。

    其次,基于不同預測時間下的預測值絕對平均誤差的對比分析,得出海表溫度和海表鹽度的最佳預測周期分別為1 d和3 d。同時證明了RBF神經(jīng)網(wǎng)絡針對不同預測周期,預測效果始終優(yōu)于BP神經(jīng)網(wǎng)絡。

    最后,本文提出并配置各項訓練參數(shù)的RBF神經(jīng)網(wǎng)絡模型,在多站點的海表溫度和海表鹽度預測實驗中,表現(xiàn)出在多個不同的海域,對海表溫度和海表鹽度普遍適用的精確預測能力。

    后續(xù)將進一步開展年際預測精度變化與全球變化之間關聯(lián)性研究,以及對海洋數(shù)據(jù)預測領域中更多海洋數(shù)據(jù)如輻照度、風向等進行預測。

    猜你喜歡
    海表鹽度神經(jīng)網(wǎng)絡
    基于無人機的海表環(huán)境智能監(jiān)測系統(tǒng)設計與應用
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    融合海表溫度產(chǎn)品在渤黃東海的對比分析及初步驗證
    海洋通報(2020年6期)2020-03-19 02:10:08
    神經(jīng)網(wǎng)絡抑制無線通信干擾探究
    電子制作(2019年19期)2019-11-23 08:42:00
    太陽總輻照度對熱帶中太平洋海表溫度年代際變化的可能影響
    基于神經(jīng)網(wǎng)絡的拉矯機控制模型建立
    重型機械(2016年1期)2016-03-01 03:42:04
    鹽度和pH對細角螺耗氧率和排氨率的影響
    復數(shù)神經(jīng)網(wǎng)絡在基于WiFi的室內LBS應用
    鹽度脅迫對入侵生物福壽螺的急性毒性效應
    適用于高鹽度和致密巖層驅油的表面活性劑
    亚洲一卡2卡3卡4卡5卡精品中文| 国产精品九九99| 一级a爱视频在线免费观看| 国产麻豆69| 欧美成狂野欧美在线观看| 夫妻性生交免费视频一级片| 亚洲精品美女久久av网站| 搡老岳熟女国产| 欧美av亚洲av综合av国产av| 少妇粗大呻吟视频| 黄片播放在线免费| 黄频高清免费视频| 亚洲精品中文字幕在线视频| 2018国产大陆天天弄谢| 午夜福利免费观看在线| 侵犯人妻中文字幕一二三四区| 人人妻人人澡人人看| 蜜桃在线观看..| 2018国产大陆天天弄谢| 国产免费现黄频在线看| 久久亚洲国产成人精品v| 波多野结衣av一区二区av| 三上悠亚av全集在线观看| 婷婷色av中文字幕| 激情视频va一区二区三区| 又紧又爽又黄一区二区| 99国产精品一区二区三区| 精品久久久精品久久久| 99久久综合免费| 色播在线永久视频| 99九九在线精品视频| 国产日韩一区二区三区精品不卡| 国产精品免费大片| 国产精品免费大片| 在线av久久热| 性少妇av在线| 日本色播在线视频| 后天国语完整版免费观看| 亚洲熟女精品中文字幕| 十八禁人妻一区二区| 日韩av不卡免费在线播放| 永久免费av网站大全| 最近最新中文字幕大全免费视频 | 一区二区日韩欧美中文字幕| 青草久久国产| 这个男人来自地球电影免费观看| 在线天堂中文资源库| 免费少妇av软件| 亚洲中文字幕日韩| 男人爽女人下面视频在线观看| 美国免费a级毛片| 国产成人精品久久二区二区免费| 久久影院123| 熟女少妇亚洲综合色aaa.| 国产精品 欧美亚洲| 亚洲精品乱久久久久久| 亚洲国产av新网站| 在线观看免费午夜福利视频| 美女主播在线视频| 69精品国产乱码久久久| 国产一区二区三区av在线| 汤姆久久久久久久影院中文字幕| 欧美亚洲 丝袜 人妻 在线| 婷婷色麻豆天堂久久| 大香蕉久久成人网| 日韩大码丰满熟妇| 一级片免费观看大全| 男女边摸边吃奶| cao死你这个sao货| 国产精品久久久久久人妻精品电影 | 精品福利观看| 美国免费a级毛片| av电影中文网址| 美国免费a级毛片| 人妻人人澡人人爽人人| 国产在视频线精品| 国产在视频线精品| 手机成人av网站| 午夜激情av网站| 777久久人妻少妇嫩草av网站| 午夜激情av网站| 国产亚洲一区二区精品| 久久久精品区二区三区| 国产成人欧美| 久久久久久久国产电影| 91字幕亚洲| 亚洲国产精品国产精品| 自线自在国产av| 嫩草影视91久久| 亚洲熟女毛片儿| 在线观看www视频免费| 国产精品国产三级国产专区5o| 在线观看www视频免费| 好男人电影高清在线观看| 日韩制服丝袜自拍偷拍| 日本五十路高清| 亚洲欧美色中文字幕在线| 日韩精品免费视频一区二区三区| 久久精品aⅴ一区二区三区四区| 午夜免费男女啪啪视频观看| 欧美激情 高清一区二区三区| 亚洲精品美女久久久久99蜜臀 | 国产男女内射视频| 啦啦啦在线观看免费高清www| 一级,二级,三级黄色视频| 久久av网站| 飞空精品影院首页| 日韩免费高清中文字幕av| 亚洲综合色网址| 激情五月婷婷亚洲| 久久精品人人爽人人爽视色| 伊人久久大香线蕉亚洲五| 久久 成人 亚洲| 亚洲综合色网址| 亚洲美女黄色视频免费看| 久久久国产精品麻豆| 午夜免费观看性视频| 精品国产超薄肉色丝袜足j| 精品卡一卡二卡四卡免费| 777米奇影视久久| 一边亲一边摸免费视频| 欧美日韩视频高清一区二区三区二| 99国产精品99久久久久| 高清欧美精品videossex| 一级a爱视频在线免费观看| 老汉色∧v一级毛片| 少妇精品久久久久久久| 亚洲视频免费观看视频| 女人久久www免费人成看片| 少妇裸体淫交视频免费看高清 | 国产熟女欧美一区二区| 国产成人精品久久二区二区91| 亚洲人成电影观看| 国产极品粉嫩免费观看在线| 国产精品亚洲av一区麻豆| www.精华液| 色播在线永久视频| 亚洲 欧美一区二区三区| 美女福利国产在线| 亚洲精品国产av成人精品| 丝袜脚勾引网站| 国产无遮挡羞羞视频在线观看| av线在线观看网站| 久久精品熟女亚洲av麻豆精品| 老汉色av国产亚洲站长工具| 大陆偷拍与自拍| 男的添女的下面高潮视频| 成人手机av| 久久精品熟女亚洲av麻豆精品| 一本一本久久a久久精品综合妖精| 成人亚洲欧美一区二区av| 中国美女看黄片| 亚洲欧美成人综合另类久久久| 成人国语在线视频| 午夜久久久在线观看| 黄色怎么调成土黄色| 国产一区亚洲一区在线观看| 这个男人来自地球电影免费观看| 精品国产一区二区三区久久久樱花| 久久天堂一区二区三区四区| 91字幕亚洲| 亚洲欧美一区二区三区国产| av国产精品久久久久影院| av网站在线播放免费| av欧美777| 天天躁夜夜躁狠狠躁躁| 精品少妇一区二区三区视频日本电影| 美女福利国产在线| 9色porny在线观看| 色婷婷av一区二区三区视频| 久久久亚洲精品成人影院| 欧美激情极品国产一区二区三区| 午夜免费男女啪啪视频观看| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品中文字幕在线视频| 欧美精品啪啪一区二区三区 | 国产黄频视频在线观看| 一本—道久久a久久精品蜜桃钙片| 多毛熟女@视频| 99九九在线精品视频| 麻豆乱淫一区二区| 欧美精品高潮呻吟av久久| 热99久久久久精品小说推荐| 夫妻性生交免费视频一级片| 十八禁人妻一区二区| 在线av久久热| 中文字幕av电影在线播放| 一边摸一边抽搐一进一出视频| 建设人人有责人人尽责人人享有的| 国产成人系列免费观看| 国产欧美日韩精品亚洲av| 肉色欧美久久久久久久蜜桃| 久久99精品国语久久久| 中文字幕色久视频| 人人妻人人添人人爽欧美一区卜| www日本在线高清视频| 丝袜在线中文字幕| 精品视频人人做人人爽| 99精国产麻豆久久婷婷| 国产黄色免费在线视频| 国产成人欧美在线观看 | 中文字幕高清在线视频| 狂野欧美激情性bbbbbb| 女人高潮潮喷娇喘18禁视频| 国产午夜精品一二区理论片| 色网站视频免费| 麻豆国产av国片精品| 欧美中文综合在线视频| 久久国产亚洲av麻豆专区| 免费看av在线观看网站| 精品视频人人做人人爽| 波野结衣二区三区在线| 日韩av不卡免费在线播放| 亚洲伊人久久精品综合| 黄片播放在线免费| 欧美人与性动交α欧美软件| 成人国产一区最新在线观看 | 黄片播放在线免费| 国产黄色视频一区二区在线观看| 亚洲伊人色综图| 午夜福利,免费看| 大话2 男鬼变身卡| 精品一区二区三区四区五区乱码 | 欧美日韩精品网址| 新久久久久国产一级毛片| 久久鲁丝午夜福利片| 18在线观看网站| 婷婷丁香在线五月| 免费一级毛片在线播放高清视频 | 无遮挡黄片免费观看| 午夜福利免费观看在线| 日韩一本色道免费dvd| 久久久久国产一级毛片高清牌| 99久久99久久久精品蜜桃| 国产亚洲一区二区精品| tube8黄色片| 大香蕉久久网| 亚洲激情五月婷婷啪啪| 少妇的丰满在线观看| 99国产精品99久久久久| 日韩中文字幕欧美一区二区 | 人人妻人人澡人人看| 亚洲av电影在线观看一区二区三区| 欧美性长视频在线观看| 97在线人人人人妻| 国产欧美亚洲国产| 婷婷色麻豆天堂久久| 久久久国产一区二区| 亚洲精品av麻豆狂野| 啦啦啦在线免费观看视频4| 一边摸一边做爽爽视频免费| 男女无遮挡免费网站观看| 免费观看人在逋| 久久精品久久久久久噜噜老黄| 国产人伦9x9x在线观看| 亚洲欧美日韩另类电影网站| 美女中出高潮动态图| 日韩熟女老妇一区二区性免费视频| 国产色视频综合| av一本久久久久| 欧美精品亚洲一区二区| 青春草视频在线免费观看| av网站免费在线观看视频| av又黄又爽大尺度在线免费看| 91九色精品人成在线观看| 亚洲专区中文字幕在线| 国产一区二区三区综合在线观看| 久久性视频一级片| 久久久久久久久免费视频了| 亚洲中文字幕日韩| 国产一卡二卡三卡精品| 国产精品二区激情视频| 男女下面插进去视频免费观看| 爱豆传媒免费全集在线观看| 无限看片的www在线观看| 人人妻,人人澡人人爽秒播 | 免费高清在线观看视频在线观看| 亚洲av电影在线观看一区二区三区| 欧美国产精品va在线观看不卡| 亚洲第一av免费看| 亚洲欧美精品自产自拍| 亚洲免费av在线视频| 国产99久久九九免费精品| 晚上一个人看的免费电影| 1024香蕉在线观看| 国产一区二区三区av在线| 国产精品一区二区免费欧美 | 午夜视频精品福利| 大香蕉久久网| 日本vs欧美在线观看视频| 精品国产国语对白av| 国产日韩欧美亚洲二区| 黑丝袜美女国产一区| 老司机午夜十八禁免费视频| 色播在线永久视频| 搡老岳熟女国产| 1024视频免费在线观看| 美女扒开内裤让男人捅视频| 欧美日韩精品网址| 赤兔流量卡办理| 男女边摸边吃奶| 成人影院久久| 日韩中文字幕欧美一区二区 | 成人三级做爰电影| 一边摸一边抽搐一进一出视频| 日韩 欧美 亚洲 中文字幕| 欧美激情高清一区二区三区| 国产精品人妻久久久影院| 五月天丁香电影| 久久亚洲精品不卡| 日韩人妻精品一区2区三区| 国产人伦9x9x在线观看| 97精品久久久久久久久久精品| 欧美黑人精品巨大| 国产成人影院久久av| 免费在线观看日本一区| 中文字幕制服av| 香蕉丝袜av| 午夜久久久在线观看| 秋霞在线观看毛片| 十八禁网站网址无遮挡| 中文字幕人妻熟女乱码| 午夜福利视频精品| 精品一区二区三卡| 欧美日韩成人在线一区二区| 亚洲色图综合在线观看| 中文字幕亚洲精品专区| 亚洲人成77777在线视频| 亚洲成国产人片在线观看| 国产片特级美女逼逼视频| 久久天堂一区二区三区四区| 91国产中文字幕| 啦啦啦在线免费观看视频4| 伊人久久大香线蕉亚洲五| 免费观看av网站的网址| 国产成人啪精品午夜网站| 亚洲熟女毛片儿| 性高湖久久久久久久久免费观看| 色婷婷av一区二区三区视频| 一本久久精品| 天天影视国产精品| 这个男人来自地球电影免费观看| 91精品国产国语对白视频| 国产野战对白在线观看| 高清黄色对白视频在线免费看| 亚洲国产中文字幕在线视频| 亚洲欧美精品综合一区二区三区| 男女床上黄色一级片免费看| 99香蕉大伊视频| 丝袜美足系列| 国产成人系列免费观看| 婷婷色麻豆天堂久久| 国产成人免费无遮挡视频| 女人高潮潮喷娇喘18禁视频| 丝袜美腿诱惑在线| 日本猛色少妇xxxxx猛交久久| 欧美日韩精品网址| 亚洲成人国产一区在线观看 | 另类精品久久| 欧美激情极品国产一区二区三区| 首页视频小说图片口味搜索 | 国产成人精品无人区| 天堂8中文在线网| 久久精品国产a三级三级三级| 精品国产乱码久久久久久男人| 777米奇影视久久| 亚洲精品在线美女| 岛国毛片在线播放| 一本—道久久a久久精品蜜桃钙片| 国产视频一区二区在线看| 久久鲁丝午夜福利片| 女人精品久久久久毛片| 又紧又爽又黄一区二区| 99国产精品免费福利视频| 国产一区亚洲一区在线观看| 午夜两性在线视频| 午夜福利免费观看在线| 亚洲国产欧美在线一区| 日韩一区二区三区影片| 国产精品.久久久| 国产精品人妻久久久影院| 亚洲成人免费电影在线观看 | 啦啦啦在线免费观看视频4| 久久国产精品男人的天堂亚洲| 国产精品偷伦视频观看了| 午夜日韩欧美国产| 精品人妻在线不人妻| 操出白浆在线播放| 国产老妇伦熟女老妇高清| 国产在线视频一区二区| 欧美黄色淫秽网站| 黄网站色视频无遮挡免费观看| 国产高清不卡午夜福利| 久久久久精品国产欧美久久久 | 久久人人爽av亚洲精品天堂| 成人国产一区最新在线观看 | 欧美性长视频在线观看| 亚洲男人天堂网一区| 亚洲av男天堂| 亚洲三区欧美一区| 夫妻性生交免费视频一级片| 久久精品国产a三级三级三级| 99九九在线精品视频| 成人国语在线视频| 美女国产高潮福利片在线看| 久久av网站| 国产极品粉嫩免费观看在线| 日韩制服骚丝袜av| 亚洲午夜精品一区,二区,三区| 国产在线视频一区二区| 国产成人免费观看mmmm| 在线观看www视频免费| 亚洲av日韩在线播放| 中文字幕av电影在线播放| 无遮挡黄片免费观看| 国产欧美日韩综合在线一区二区| 一二三四在线观看免费中文在| 成人亚洲精品一区在线观看| 成人影院久久| 免费少妇av软件| 亚洲欧洲国产日韩| 黄网站色视频无遮挡免费观看| 免费高清在线观看日韩| 人人澡人人妻人| avwww免费| 性少妇av在线| 国产精品久久久久久精品古装| 国产女主播在线喷水免费视频网站| 精品高清国产在线一区| 国产一卡二卡三卡精品| 国产男人的电影天堂91| 午夜精品国产一区二区电影| 欧美日韩成人在线一区二区| 女人爽到高潮嗷嗷叫在线视频| 国产成人一区二区在线| 免费久久久久久久精品成人欧美视频| av线在线观看网站| 如日韩欧美国产精品一区二区三区| 老司机在亚洲福利影院| 欧美日韩亚洲综合一区二区三区_| 久久毛片免费看一区二区三区| 十分钟在线观看高清视频www| 国产欧美日韩一区二区三区在线| 国产黄色免费在线视频| www日本在线高清视频| 久久精品人人爽人人爽视色| 国产成人免费观看mmmm| www日本在线高清视频| 叶爱在线成人免费视频播放| 国产av国产精品国产| 1024香蕉在线观看| 日韩av免费高清视频| 无遮挡黄片免费观看| 老司机午夜十八禁免费视频| 中文字幕亚洲精品专区| 在线观看免费日韩欧美大片| 午夜激情av网站| 日本a在线网址| 精品国产国语对白av| 啦啦啦 在线观看视频| 亚洲av电影在线观看一区二区三区| 美女中出高潮动态图| 一级毛片 在线播放| 丝瓜视频免费看黄片| av福利片在线| 啦啦啦啦在线视频资源| 视频在线观看一区二区三区| 下体分泌物呈黄色| 亚洲欧美一区二区三区黑人| 欧美日韩国产mv在线观看视频| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 国产欧美日韩一区二区三 | 久久精品亚洲av国产电影网| 欧美日韩精品网址| 亚洲欧美一区二区三区黑人| 亚洲精品国产av成人精品| 色94色欧美一区二区| 亚洲色图综合在线观看| 国产精品二区激情视频| 久久影院123| 少妇被粗大的猛进出69影院| 国产激情久久老熟女| 天天影视国产精品| 久久久久久亚洲精品国产蜜桃av| 18禁观看日本| 亚洲五月婷婷丁香| 久久久久国产精品人妻一区二区| 欧美激情高清一区二区三区| av国产精品久久久久影院| 国产淫语在线视频| 精品久久蜜臀av无| 波多野结衣av一区二区av| 国产成人91sexporn| av线在线观看网站| 老司机午夜十八禁免费视频| 亚洲图色成人| 亚洲欧洲国产日韩| 亚洲伊人色综图| 午夜老司机福利片| 亚洲 国产 在线| 亚洲天堂av无毛| 超碰97精品在线观看| 天天操日日干夜夜撸| 国产av国产精品国产| 欧美日本中文国产一区发布| 脱女人内裤的视频| 欧美日本中文国产一区发布| 日日爽夜夜爽网站| 啦啦啦 在线观看视频| 在线观看免费高清a一片| 国产精品九九99| 久久 成人 亚洲| 亚洲中文av在线| 中文字幕高清在线视频| 久久精品人人爽人人爽视色| 99九九在线精品视频| 欧美激情高清一区二区三区| 不卡av一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 国产在视频线精品| 精品国产乱码久久久久久男人| 不卡av一区二区三区| cao死你这个sao货| 精品福利观看| 1024视频免费在线观看| 韩国高清视频一区二区三区| 久久性视频一级片| 欧美日韩黄片免| 国产熟女午夜一区二区三区| 久久国产精品影院| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久精品精品| 欧美黑人精品巨大| 亚洲国产av影院在线观看| 日韩一区二区三区影片| 精品国产一区二区三区久久久樱花| 日韩av免费高清视频| 亚洲美女黄色视频免费看| 夜夜骑夜夜射夜夜干| 精品少妇黑人巨大在线播放| 欧美日韩亚洲国产一区二区在线观看 | 水蜜桃什么品种好| 18在线观看网站| 亚洲精品中文字幕在线视频| 亚洲色图 男人天堂 中文字幕| 久久久久久人人人人人| 国产精品秋霞免费鲁丝片| 亚洲av在线观看美女高潮| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 午夜激情久久久久久久| 久久久欧美国产精品| 欧美久久黑人一区二区| 国产一区二区激情短视频 | 国产黄色视频一区二区在线观看| 麻豆av在线久日| 久久久精品区二区三区| 久久精品久久久久久噜噜老黄| 亚洲精品国产av蜜桃| 久久久久久亚洲精品国产蜜桃av| 日韩电影二区| 国产亚洲午夜精品一区二区久久| 1024香蕉在线观看| 啦啦啦在线免费观看视频4| av有码第一页| 美国免费a级毛片| 久久99一区二区三区| 日韩 亚洲 欧美在线| 国产精品偷伦视频观看了| 晚上一个人看的免费电影| 国产成人av教育| 无遮挡黄片免费观看| svipshipincom国产片| 免费在线观看完整版高清| 国产成人a∨麻豆精品| 欧美日韩亚洲综合一区二区三区_| 亚洲午夜精品一区,二区,三区| 久久99精品国语久久久| 欧美 日韩 精品 国产| 女人精品久久久久毛片| 一区二区日韩欧美中文字幕| 校园人妻丝袜中文字幕| 老司机影院毛片| 亚洲精品日本国产第一区| 免费人妻精品一区二区三区视频| 少妇粗大呻吟视频| 一本色道久久久久久精品综合| 精品一品国产午夜福利视频| 久久国产精品人妻蜜桃| 90打野战视频偷拍视频| 亚洲午夜精品一区,二区,三区| 中文欧美无线码| 久久精品国产a三级三级三级| 一级片'在线观看视频| 宅男免费午夜| 久久国产亚洲av麻豆专区| 国产伦人伦偷精品视频| 亚洲精品第二区| 人人妻人人澡人人看| 午夜影院在线不卡| 久久ye,这里只有精品| 成年人免费黄色播放视频| 美国免费a级毛片| 亚洲一区二区三区欧美精品| 国产爽快片一区二区三区| 狂野欧美激情性xxxx| 国产又色又爽无遮挡免| 国产av国产精品国产| 亚洲欧美一区二区三区黑人| 国产精品一二三区在线看| 亚洲欧洲国产日韩| 精品一品国产午夜福利视频| 久久精品久久精品一区二区三区| 国产av精品麻豆|