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

    基于BP神經網絡方法的近岸數(shù)值海溫預報釋用技術*

    2016-01-15 11:22:18匡曉迪王兆毅張苗茵何恩業(yè)鄧小花
    海洋與湖沼 2016年6期
    關鍵詞:海溫臺站數(shù)值

    匡曉迪 王兆毅 張苗茵 何恩業(yè) 鄧小花

    (國家海洋環(huán)境預報中心 北京 100086)

    隨著數(shù)值預報和同化技術的快速發(fā)展, 從全球到近岸的不同水平分辨率、不同動力過程、不同預報精度的海洋數(shù)值預報模型不斷被建立, 并具備一定的精度(李冬等, 2004; 張權等, 2011; 李云等, 2011;方長芳等, 2013; 魏澤勛等, 2015; 王輝, 2016)。然而,區(qū)域海洋數(shù)值預報受地形精度、物理過程描述及近岸觀測數(shù)據(jù)等因素的限制, 在近岸和淺水區(qū)預報表現(xiàn)受到制約, 對定點海溫預報準確度不高。對于有觀測資料的站點, 定點數(shù)值預報與觀測的結合大多只通過人工經驗判定, 往往不能同時滿足穩(wěn)定的預報表現(xiàn)和客觀計算兩方面的需求。因此對近岸海溫預報來說, 數(shù)值預報的結果不能直接作為預報產品進行發(fā)布。

    預報釋用技術是數(shù)值預報和預報產品發(fā)布之間的橋梁。該技術可以充分挖掘數(shù)值預報產品中的有用信息, 進一步提高其精度, 實現(xiàn)預報產品的覆蓋率、客觀性和準確性的同步提升, 突出預報產品的客觀解釋應用特征, 使其具有優(yōu)良的服務體驗、更好的適用性和更高的準確性。這些釋用技術基本上可以歸納為兩大類: 一類是以人的經驗為主的方法(嚴明良,2004; 堵盤軍等, 2013); 另一類是客觀定量方法, 如動力釋用方法(苗世光等, 2009; 胡波等, 2014)、統(tǒng)計釋用方法(程戴暉等, 2001; 張宏秋等, 2004; 胡春梅等, 2006; 李江萍等, 2008; 曾瑾瑜等, 2015)等。

    人工神經網絡方法是統(tǒng)計釋用方法中一種非線性擬合方法。該方法是人腦及其活動的一個理論化的數(shù)學模型, 是一個大規(guī)模的非線性自適應系統(tǒng)(蔣宗禮, 2002)。近些年, 人工神經網絡在海洋預報中得到廣泛運用, 水位(Filippo et al, 2012)、鹽度(張文孝等,2006)、海浪(齊義泉等, 2005; 王華等, 2007;Deshmukh et al, 2016)、風暴潮(Tissot et al, 2002) 、赤潮(黃良民等, 1994; 吳京洪等, 2001; 楊建強等,2003; 何恩業(yè)等, 2008)、潮汐(王建華等, 2007; 歐素英等, 2008)、海冰(盧海, 2007)等預報及氣候預測(張韌等, 2000; 王業(yè)宏等, 2003; 周林等, 2009; Liu et al,2016)和垂直結構計算(劉輝等, 2008; Chaudhari et al,2008; Zhao et al, 2015)等領域已有初步應用。Vilibi?等(2016)采用自組織映射網絡, 基于地波雷達和氣象數(shù)值預報建立了實時海流預報系統(tǒng), 通過與主流區(qū)域海洋模式 ROMS比較, 認為神經網絡方法建立的模型與基于海洋動力方程建立的模型精度相當甚至略優(yōu), 給基于人工智能的近岸海洋要素定量預報開辟了新的思路。對于近岸定點海溫的定量預報, 由于近岸海洋是強非線性復雜系統(tǒng), 地域和基礎數(shù)據(jù)依賴性明顯, 使得數(shù)值預報具有較大的不確定性, 因此預報一般采用經驗方法或者分型統(tǒng)計方法(張宏秋等,2004), 采用神經網絡方法尚無前例?;谌斯ど窠浘W絡的誤差后傳算法(back-propagation neural network,簡稱BP)是人工神經網絡中應用最為廣泛的一種多層前饋網絡的學習算法, 它可以通過神經網絡的自學習功能, 確定神經元之間的耦合權值, 從而使網絡整體具有近似函數(shù)的功能, 非常適用于非線性系統(tǒng)的建模研究(張宏秋等, 2004)。在理論上, BP模型能高精度逼近任意復雜系統(tǒng), 可有效用于常規(guī)方法不易處理的關系辨識和復雜預測, 并為建立合理、可靠和準確的預報模型提供依據(jù)(Vilibi? et al, 2016)。對于確定的海洋要素, 通過不同的因子場得到多個典型因子,進而通過BP神經網絡實現(xiàn)典型因子場與海洋要素之間的非線性建模, 在海洋預報中也有成功提高預報精度的實例(齊義泉等, 2005; 何恩業(yè)等, 2008; Tisso et al, 2014; Deshmukh et al, 2016)。本文以提高數(shù)值預報的近岸精細化預報精度為目的, 將近岸臺站的定點表層海溫作為預報對象, 開展基于 BP神經網絡的預報釋用技術研究, 旨在解決海溫數(shù)值預報產品在近岸、單點預報誤差大的問題, 有針對性地提高數(shù)值海溫的預報精度和質量。第一節(jié)介紹預報釋用建模使用的數(shù)值預報及觀測數(shù)據(jù), 簡述了BP神經網絡方法的基本結構、算法和比測實驗方案; 第二節(jié)基于各組調訓實驗的誤差分析, 研究了固定影響因子和一般影響因子的作用和權重, 分析了釋用模型對預報精度的提升效果, 開展預報檢驗, 并與人工經驗預報誤差進行了比較, 篩選出了最佳預報模型。

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

    1.1 研究數(shù)據(jù)

    1.1.1 氣象預報數(shù)據(jù) 數(shù)值氣象資料為國家海洋環(huán)境預報中心的 WRF風場預報網格資料, 調訓模型數(shù)據(jù)時間范圍為2014年1月1日至12月31日, 預報實驗數(shù)據(jù)時間長度為2015年1月1日至4月30日。每天 8、14、20時 3時次(北京時)進行平均, 共 365個樣本, 網格精度為 1/12°×1/12°, 利用雙線性插值,將數(shù)據(jù)插值到臺站所在位置。具體氣象要素選取了海面 2米比濕(Q2)、海面氣壓(SLP)、凈熱通量(SWDOWN)、海面2米氣溫(T2), 緯向10米風速(U10)、經向10米風速(V10)、10米風速大小(VEL)7個要素。1.1.2 臺站觀測數(shù)據(jù) 臺站觀測的主要觀測要素涵蓋潮汐、表層水溫、表層鹽度、海浪、風向、風速、氣壓、氣溫、相對濕度、能見度和降水量等。臺站依據(jù)《海濱觀測規(guī)范》、《地面氣象觀測規(guī)范》和《海洋自動化觀測通用技術要求》等觀測工作執(zhí)行標準, 開展各觀測項目和要素的數(shù)據(jù)采集處理、傳輸, 其觀測資料能在一定程度上反映出觀測海區(qū)環(huán)境的基本特征和變化規(guī)律。本文使用舟山站的常規(guī)臺站觀測資料中的海溫和氣溫數(shù)據(jù), 分別作為預報因子(調訓和預報)和檢驗數(shù)據(jù)使用。調訓使用的觀測資料時間長度為2014年1月1日至12月31日, 預報實驗使用的觀測資料時間長度為2015年1月1日至4月30日。每天8、14、20時3個時次(北京時)進行平均, 共365個樣本, 觀測要素和具體影響要素選取了觀測氣溫(Tao)和表層海溫(Tso)。

    1.1.3 海溫預報數(shù)據(jù) 數(shù)值海溫采用國家海洋環(huán)境預報中心開發(fā)的全球業(yè)務化海洋學系統(tǒng)(Chinese Global Ocean Forecasting System, CGOFS), 調訓模型數(shù)據(jù)時間長度為2014年1月1日至12月31日, 預報實驗數(shù)據(jù)時間長度為2015年1月1日至4月30日。釋用采用每天20時(北京時)的數(shù)值模擬海溫場, 共365個樣本, 檢驗采用每天8、14、20時3個時次(北京時)進行平均, 網格精度為 1°/30×1°/30, 選取離站點最近的模式水點代表臺站的海溫數(shù)值模擬結果, 并計算了位置誤差。

    1.1.4 經驗預報 采用2014年全年時效為24小時的日平均海溫經驗預報值的時間序列和臺站觀測海溫對比, 并采用與釋用預報相同方法進行誤差分析,與釋用預報的結果進行對比。

    1.2 人工神經網絡方法

    BP網絡是一種非線性映射人工神經網絡, 具有很強的學習和記憶能力, 一般由一個輸入層、若干隱含層和一個輸出層組成。它是一種有導師的學習算法,通常是應用最速下降法來調整各層權值, 并最終達到預先給定的誤差為止。本文采用單隱層BP網絡建立非線性回歸模型, 隱層神經元采用S型對數(shù)激活函數(shù)(單極性 Sigmoid函數(shù)), 輸出層神經元采用線性激活函數(shù), 為了有利于使網絡聯(lián)接權加快脫離誤差曲面的平坦區(qū), 引進動量項。

    1.2.1 BP網絡結構模型

    單隱層BP網絡結構如圖1所示, 共分3層:

    第1層為變量輸入層, xj, j=1, 2, …n0為輸入變量,n0為輸入變量的個數(shù)。若設x0為第一隱層中激活函數(shù)的域值, 則輸入向量總共為 n0+1維。x0一般取為-1,將其增廣到輸入量中, 作為一個分量, 則有

    第2層為隱層, 設有n1個神經元, 則其輸出向量為g=(g0, g1, g2, …,1ng ), 其中g0為第一隱層中激活函數(shù)的域值, 一般取為-1。

    第3層為輸出層, 設有m個神經元, 網絡的輸出向量為y=(y0,y1,y2,…, ym) 。

    圖1 單隱層BP網絡結構Fig.1 Structure of single-hidden-layer BP network

    隱層神經元采用 S型對數(shù)激活函數(shù)(單極性Sigmoid函數(shù)), 輸出層神經元采用線性激活函數(shù)。S型對數(shù)激活函數(shù)表達式為:

    隱層輸出向量為:

    輸出層輸出向量為:

    1.2.2 BP網絡權值的動態(tài)學習算法 設一共有α個樣本對, 第 p個樣本對(pptx, )輸入模式后輸出方差為:

    α個樣本對經正向傳播運算后網絡總誤差為:

    輸入層到隱層的聯(lián)接權修正量為:

    設n為迭代次數(shù), 應用梯度法可以得到網絡各層聯(lián)接權的迭代關系為:

    為了有利于使網絡聯(lián)接權加快脫離誤差曲面的平坦區(qū), 引進動量項mcΔwuv, 則以上方程組變?yōu)?

    上面各式中, η表示學習率, 其值通常在0.01— 1.0之間。學習率η選得太小, 會導致網絡參數(shù)(耦合權值)修改量過小, 收斂緩慢; 選得太大, 雖然可以加快學習速度, 但可能致使權值修改量在穩(wěn)定點附近持續(xù)振蕩, 難以收斂, mc為動量因子, 一般取0.9左右。動量項的作用在于記憶前一時刻聯(lián)接權值的變化方向(即變化量的值), 增加動量項, 利用其“慣性效應”來抑制可能產生的振蕩, 起到平滑作用, 這樣可以采用較大的學習率η, 以提高學習速度。本文調訓模型目標函數(shù)誤差為0.005, 學習率0.2, 動量因子取0.9。

    1.2.3 集合樣本的構造 由于研究使用的各數(shù)據(jù)單位不一致, 對 BP建模效果存在較大影響, 在影響因子的樣本序列輸入神經網絡前, 需要對集合樣本進行去量綱化、歸一化和類間交叉處理; 本文使用的去量綱化和歸一化計算方法為:

    其中X為集合樣本, Xmin為樣本最小值, Xmax為樣本最大值, σX為樣本方差, x為去量綱化和歸一化后的集合樣本。同時, 如果給定的樣本要素隨著時間有整體變化趨勢, 不利于模型的泛化, 因此樣本的組織要注意將不同類別的樣本交叉輸入, 對樣本的類間交叉預處理, 本文采用隨機數(shù)方法對集合樣本進行打亂。

    調訓樣本為2014年1月1日至12月31日共365個樣本, 預報樣本為2015年1月1日至4月30共120個樣本。其中由于觀測數(shù)據(jù)具有滯后性, 在構建觀測海溫和氣溫樣本時, 以第T-1日的觀測數(shù)據(jù)和第T日的數(shù)值預報構建樣本集合。

    在BP神經網絡計算結束輸出結果后, 進行反算恢復樣本量綱、量級和原有序列。

    1.2.4 計算步驟 單隱層BP網絡反傳學習算法的實現(xiàn)分為5步:

    (1) 初始化。設置網絡各聯(lián)接權值(包括閾值), 置所有可調參數(shù)(權和閾值)為均勻分布的較小隨機數(shù),使每個神經元的凈輸入值較小, 以保證工作在激活函數(shù)變化較靈敏的區(qū)域。

    (2) 信息正向傳遞過程的計算。輸入學習樣本對(xP, tP), 對每個輸入樣本計算各神經元的輸出值。

    (3) 誤差反向傳播。從輸出層至輸入層, 反向依次計算各神經元的等效誤差, 然后返回步驟(2)對其他學習樣本對進行正向傳播計算和誤差反向傳播,一直到所有α個學習樣本對都進行類似的運算為止。

    (4) 調整各層的聯(lián)接權值。按照權值調整公式,修改各層的聯(lián)接權值。

    (5) 返回步驟(2), 根據(jù)新的聯(lián)接權值, 進行正向計算。若對每一個學習樣本對(xP, tP),和輸出層的每一個神經元l均滿足精度要求, 即:

    1.3 研究方案

    本文選用舟山站為例進行預報因子的敏感性試驗、調訓和預報模型建模。首先將影響因子分為固定因子(數(shù)值和觀測海溫, 共2個)和一般因子(數(shù)值氣象7個因子和觀測氣溫1個因子, 共8個)兩類, 設計4組敏感性試驗: 控制實驗組 NONE(無固定影響因子)、數(shù)值海溫組TM(數(shù)值海溫為固定影響因子)、觀測海溫組 TO(觀測海溫為固定影響因子)、數(shù)值及觀測海溫組TMTO(兩者均作為固定的影響因子)進行模型調訓。每組調訓實驗中, 將一般影響因子按 5—8個因子數(shù)進行排列組合后, 加上各組的固定影響因子, 共設計了 93組影響因子, 按照調訓模型的誤差(如無特別說明, 指均方根誤差, 下同)分析進行排序, 并比較不同實驗組相同一般影響因子實驗之間的誤差。

    (1) 組1: 控制實驗組NONE(93組):

    實驗數(shù)據(jù): 第T日數(shù)值氣象+第T-1日的臺站氣溫觀測

    氣象變量: Q2, SLP, SWDOWN, T2, U10, V10,VEL(比濕, 海面氣壓, 凈熱通量, 氣溫, 緯向、經向風速, 風速)

    臺站觀測: Tao(氣溫, 下同)

    (2) 組2: 數(shù)值海溫組TM(93組)

    實驗數(shù)據(jù): 第T日數(shù)值氣象+第T-1日的臺站氣溫觀測+海溫數(shù)值預報

    氣象變量: Q2, SLP, SWDOWN, T2, U10, V10, VEL;

    臺站觀測: Tao;

    模型數(shù)據(jù): Tm(數(shù)值海溫)

    (3) 組3: 觀測海溫組TO(93組)

    實驗數(shù)據(jù): 第T日數(shù)值氣象+第T-1日的臺站氣溫觀測+第T-1日臺站海溫觀測

    氣象變量: Q2, SLP, SWDOWN, T2, U10, V10, VEL;

    臺站觀測: Tao, Tso(觀測海溫, 下同);

    (4) 組4: 數(shù)值及觀測海溫組TMTO(93組)

    實驗數(shù)據(jù): 第T日數(shù)值氣象+海溫+第T-1日的臺站氣溫+第T-1日臺站海溫觀測;

    氣象變量: Q2, SLP, SWDOWN, T2, U10, V10, VEL;

    臺站觀測: Tao, Tso;

    模型數(shù)據(jù):Tm。

    通過分析4種釋用方案各93組配置實驗的誤差,計算了372個實驗的絕對誤差、平均誤差、平均絕對誤差和均方根誤差。對4個釋用方案的相同配置實驗進行平行比較, 對同一釋用方案下的不同配置實驗也進行誤差對比, 分析數(shù)值海溫、觀測海溫在建模中的影響和作用, 比較其他影響因子的權重, 并篩選誤差較小、穩(wěn)定性較好的調訓模型進行24小時海溫預報實驗, 同時以釋用結果與 24小時的經驗預報結果的對比結果作為參考, 選定預報釋用試行方案。

    2 調訓實驗

    2.1 固定影響因子的作用

    首先, 對4個釋用方案進行整體誤差比較, 其中NONE組誤差較大且分散, 按照NONE、TM、TO、TMTO的順序, 誤差依次減小且逐步穩(wěn)定, NONE為1.1°C, 且有顯著高誤差值出現(xiàn), 有1組實驗最大誤差高達28°C以上, 其他3個方案均無此現(xiàn)象發(fā)生; TM誤差低于NONE, 為0.9°C, 另外2組誤差接近且都低于TM; TMTO誤差略低于TO, 且整體誤差均穩(wěn)定的處在較低水平, TO誤差為0.8°C, TMTO誤差最小, 為0.7°C(表 1)。

    通過誤差分布可以看出(圖 2), 數(shù)值海溫和觀測海溫分別作為固定影響因子, 使得實驗誤差從 1.1°C分別降至 0.9°C和 0.8°C, 對降低模型誤差起到了十分重要的作用; 同時, 觀測海溫作為固定影響因子,起到了穩(wěn)定模型誤差的作用, 不同配置實驗之間誤差的差別明顯縮小。因此, 數(shù)值海溫和觀測海溫作為固定影響因子, 可以使釋用模型表現(xiàn)優(yōu)良且穩(wěn)定, 因此, TMTO方案應作為首選的釋用方案。

    圖2 釋用方案均方差分布Fig.2 Mean square error distribution of the four interpretation plans

    表1 釋用方案實驗誤差比較(單位: °C)Tab.1 Comparison in the errors among the four plans(°C)

    2.2 一般影響因子權重分析

    對 TMTO93組配置實驗的誤差進行排序, 對誤差最低的10個實驗組進行影響因子分析, 記錄前10位的實驗中, 對每個一般影響因子出現(xiàn)的頻次計數(shù),可以看出 8個一般影響因子在釋用模型中的地位和權重(見表2)??傮w上, 海面氣溫(T2)、觀測氣溫(Tao)和海表面氣壓(SLP)的出現(xiàn)頻次較高; TMTO方案中,凈熱通量(SWDOWN)、海面風速(VEL)和海面氣壓(SLP)權重較大。

    表2 影響因子在誤差最小的前十組實驗中出現(xiàn)頻次Tab.2 Occurency number of affectors in the best ten combinations

    2.3 預報釋用與數(shù)值預報

    數(shù)值預報本身存在較多的影響因素, 因此單點的海溫預報誤差一般不能保證, 且模式在近岸的臺站附近預報能力十分有限。如圖 3a所示, 數(shù)值預報的海溫能夠較好地表現(xiàn)海溫的季節(jié)變化, 但預報誤差較大, 均方差約為 2.2°C; 釋用模型對海溫的改善是明顯的, TM方案將預報誤差大幅減低至1.2°C, 且整體變化曲線與觀測值十分接近; TMTO方案則進一步降低了誤差, 提升了在海溫突變時的表現(xiàn), 兩條曲線的變化趨勢和量值都十分接近, 將圖3c和a相比可以看出, 釋用技術大大改善了數(shù)值海溫的預報表現(xiàn),提升了數(shù)值預報產品在指定區(qū)域或站點的預報準確性, 拓展了數(shù)值海溫產品的應用前景。

    2.4 預報釋用與經驗預報

    對同一時段的經驗預報進行誤差分析(圖 4, 剔除了兩次海溫異常峰值)發(fā)現(xiàn), 經驗預報總體誤差較小; 但如果天氣導致了短期內海溫顯著變化, 經驗預報 1—2天的預報表現(xiàn)會較差, 經驗預報的總體方差約為 0.6°C, 最大誤差也可達 4°C以上。而釋用模型可以較好地體現(xiàn)海溫的突出變化(圖3c), 總體誤差約為 0.7°C, 且誤差水平比較穩(wěn)定, 很少有誤差極大值出現(xiàn), 因此經驗預報與數(shù)值釋用水平基本持平。

    圖3 觀測、數(shù)值和釋用海溫比較Fig.3 Comparison in the results of observation, numerical modeling, and interpretation

    圖4 經驗預報與觀測值的比較Fig.4 Comparison between empirical forecast and observation on temperature

    3 預報實驗

    鑒于TMTO的93組實驗的模型調訓誤差都相差較小(前 10位誤差如表 3, 影響因子作為配置實驗名稱), 進行預報實驗時, 對每個模型都進行預報實驗及誤差分析, 檢驗模型的預報表現(xiàn)。基于 2015年1—5月的風場、數(shù)值海溫和臺站觀測數(shù)據(jù)進行了為期五個月的釋用試預報, 調訓誤差和預報誤差比較如圖5所示。

    從圖上可以看出, 各個模型的調訓誤差相對穩(wěn)定, 預報誤差比調訓誤差略高, 最小預報誤差為0.88°C, 雖然模型之間的調訓誤差差別不大(表 3),但從模型的預報表現(xiàn)來看, 預報誤差水平差別較為明顯。

    結合調訓和預報實驗結果來看, 可以選擇預報誤差最小的Q2SWDOWNT2V10Tao模型作為舟山站的海溫釋用模型, 該模型因子包含 TMTO中權重最大的短波輻射和海表氣溫, 兩者皆與海面感熱交換關系密切; 同時, 模型包含權重較大的觀測氣溫和比濕,比濕能夠影響海表蒸發(fā)過程, 從而影響海面的潛熱交換; 另外, 考慮到舟山站的地理位置, 緯向風(南風或北風)除了影響海表失熱過程, 還可以通過平流輸運帶來暖濕或干冷的空氣, 間接影響表層海溫,因此該影響因子的配置對舟山站的表層海溫均合理有效。

    表3 TMTO組前10位模型配置及誤差(單位: °C)Tab 3 Training errors of the best 10 combinations in plan TMTO

    圖5 配置實驗的調訓誤差(紅)和預報誤差(藍)Fig.5 Training errors (red) and forecast errors (blue)

    4 結論和展望

    本文以海溫為突破口, 以舟山站為示范, 開展了定點數(shù)值海溫預報釋用技術研究, 在數(shù)值預報和預報產品之間建立橋梁, 客觀地提高定點海溫預報精度, 形成了如下結論:

    (1) 基于 BP人工神經網絡的預報釋用技術,可以大幅降低近岸數(shù)值預報的誤差, 且不依賴主觀經驗;

    (2) 數(shù)值海溫和觀測海溫具有降低和穩(wěn)定調訓誤差的作用, 是近岸數(shù)值海溫預報釋用不可忽略的重要影響因子;

    (3) 通過學習算法, BP神經網絡預報誤差略高于調訓誤差, 但仍保持在合理水平, 與經驗預報的預報能力持平;

    (4) 舟山站表層海溫釋用的影響因子應包含短波輻射通量、氣溫、比濕和緯向風, 通過這些因子建立的釋用模型的預報誤差約為0.88°C。

    通過本文的研究可以看出, 基于人工神經網絡的預報釋用技術確實可以充分挖掘數(shù)值系統(tǒng)和觀測網絡中的重要信息, 同步提高預報產品的客觀性和準確性, 突出預報產品的客觀解釋應用特征, 提供具有優(yōu)良服務體驗、更好適用性和更高準確性的定點海溫預報產品, 具有十分廣闊的拓展空間和應用前景;同時, 結合海洋基本原理和物理背景的預報因子的優(yōu)選和模型調訓, 以及單點預報模型在區(qū)域預報產品中的適用性等, 仍將是這一工作繼續(xù)深化提升的難點和重點。

    王 華, 姚圣康, 珣龔茂等, 2007. 東海區(qū)域災害性海浪長期預測方法研究. 海洋通報, 26(5): 35—42

    王業(yè)宏, 金 龍, 2003. 基于自然正交展開的神經網絡長期預報模型. 自然災害學報, 12(2): 127—132

    王建華, 于紅兵, 宋運法等, 2007. 人工神經網絡在潮汐數(shù)值預報中的應用. 海洋預報, 24(2): 47—51

    王 輝, 2016. 全球海洋環(huán)境數(shù)值預報關鍵技術系統(tǒng)集成研究及應用. 中國科技成果, 17(3): 18—19

    方長芳, 張 翔, 尹建平等, 2013. 21世紀初海洋預報系統(tǒng)發(fā)展現(xiàn)狀和趨勢. 海洋預報, 30(4): 93—102

    盧 海, 2007. 基于人工神經網絡的冰情預測研究. 天津: 天津大學碩士學位論文

    劉 輝, 王 靜, 2008. 人工神經網絡方法估算海洋上混合層深度的初步研究. 熱帶海洋學報, 27(3): 9—13

    齊義泉, 張志旭, 李志偉等, 2005. 人工神經網絡在海浪數(shù)值預報中的應用. 水科學進展, 16(1): 32—35

    嚴明良, 2004. 數(shù)值產品釋用方法在預報業(yè)務系統(tǒng)中的應用.推進氣象科技創(chuàng)新加快氣象事業(yè)發(fā)展——中國氣象學2004年年會討論文集. 北京: 中國氣象協(xié)會

    李 云, 萬莉穎, 劉欽政, 2011. 衛(wèi)星高度計資料在三維海溫和鹽度數(shù)值預報中的應用. 海洋學報, 33(6): 64—72

    李 冬, 劉克修, 韓桂軍, 2004. 變分伴隨數(shù)據(jù)同化方法在斷面海溫數(shù)值計算中的應用研究. 海洋通報, 23(4): 1—8

    李江萍, 王式功, 2008. 統(tǒng)計降尺度法在數(shù)值預報產品釋用中的應用. 氣象, 34(6): 41—45

    楊建強, 羅先香, 丁德文等, 2003. 赤潮預測的人工神經網絡方法初步研究. 海洋科學進展, 21(3): 318—324

    吳京洪, 楊秀環(huán), 唐寶英等, 2001. 大亞灣澳頭增養(yǎng)殖區(qū)赤潮與環(huán)境的關系研究 I. 浮游植物總生物量與環(huán)境因子的關系. 中山大學學報(自然科學版), 40(3): 37—40

    何恩業(yè), 李 海, 任湘湘等, 2008. BP人工神經網絡在渤海灣葉綠素預測中的應用. 海洋預報, 25(2): 1—10

    張 權, 王 燕, 王關鎖等, 2011. 高分辨率浪-潮-流耦合數(shù)值預報系統(tǒng)的業(yè)務化運行與檢驗分析. 海洋科學進展,29(2): 129—136

    張 韌, 蔣國榮, 2000. 赤道東太平洋海溫及El Nino/La Nina的反演及預測. 解放軍理工大學學報(自然科學版), 1(5):7—12

    張文孝, 高國棟, 慕光宇等, 2006. 基于BP神經網絡的海水鹽度建模研究. 海洋技術, 25(4):39—41

    張宏秋, 汪麗新, 張繼贏, 2004. 夏季溫度預報綜合統(tǒng)計釋用技術. 遼寧氣象, (3): 18—19

    苗世光, 孫桂平, 馬 艷等, 2009. 青島奧帆賽高分辨率數(shù)值模式系統(tǒng)研制與應用. 應用氣象學報, 20(3): 370—379

    歐素英, 楊清書, 2008. 人工神經網絡模型在航道、港口潮水位預報中的應用. 水利水運工程學報, (2): 67—70

    周 林, 楊成蔭, 王漢杰等, 2009. 基于 CCA-BP-BPNN釋用模型的太平洋 SST預報. 解放軍理工大學學報(自然科學版), 10(4): 391—396

    胡 波, 杜惠良, 郝世峰等, 2014. 一種統(tǒng)計技術結合動力釋用的沿海海霧預報方法. 海洋預報, 31(5): 82—86

    胡春梅, 余 暉, 陳佩燕, 2006. 西北太平洋熱帶氣旋強度統(tǒng)計釋用預報方法研究. 氣象, 32(8): 64—69

    堵盤軍, 王麗琳, 管琴樂等, 2013. 人機交互式臺風浪業(yè)務預報系統(tǒng)的研制及應用. 中國海洋大學學報(自然科學版),43(10): 16—24

    黃良民, 錢宏林, 李錦蓉, 1994. 大鵬灣赤潮多發(fā)區(qū)的葉綠素α分布與環(huán)境關系初探. 海洋與湖沼, 25(2): 197—205

    蔣宗禮, 2002. 人工神經網絡導論. 北京: 高等教育出版社, 10—13

    程戴暉, 楊美川, 2001. 上海浦西地區(qū)霧持續(xù)時間的統(tǒng)計釋用預報. 氣象, 27(7): 16—19

    曾瑾瑜, 劉愛鳴, 高 珊等, 2015. 福建省沿海冬半年東北大風的數(shù)值預報釋用方法研究. 海洋預報, 32(5): 61—68

    魏澤勛, 徐騰飛, 王永剛等, 2015. 南海及周邊海域風浪流耦合同化精細化數(shù)值預報與信息服務系統(tǒng)簡介. 海洋技術學報, 34(3): 86—90

    Chaudhari S, Balasubramanian R, Gangopadhyay A, 2008.Upwelling detection in avhrr sea surface temperature (SST)images using neural-network framework. In: IEEE International Geoscience and Remote Sensing Symposium,2008. Boston, MA, USA: IEEE, IV-926—IV-929

    Deshmukh A N, Deo M C, Bhaskaran P K et al, 2016.Neural-network-based data assimilation to improve numerical ocean wave forecast. IEEE Journal of Oceanic Engineering, http://dx.doi.org/10.1109/JOE.2016.2521222 Filippo A, Torres A R Jr, Kjerfve B et al, 2012. Application of artificial neural network (ANN) to improve forecasting of sea level. Ocean & Coastal Management, 55: 101—110

    Liu S Y, Xu L Q, Li D L, 2016. Multi-scale prediction of water temperature using empirical mode decomposition with back-propagation neural networks. Computers & Electrical Engineering, 49: 1—8

    Tissot P E, Cox D T, Michaud P, 2002. Neural network forecasting of storm surges along the gulf of Mexico. In:Edge B L, Hemsley J M eds. Ocean Wave Measurement and Analysis. San Francisco, California, United States: ASCE,1535—1544

    Vilibi? I, ?epi? J, Mihanovi? H et al, 2016. Self-organizing maps-based ocean currents forecasting system. Scientific Reports, 6: 22924

    Zhao N, Han Z, 2015. A Simulation model of seawater vertical temperature by using back-propagation neural network.Polish Maritime Research, 22(S1): 82—88

    猜你喜歡
    海溫臺站數(shù)值
    用固定數(shù)值計算
    中國科學院野外臺站檔案工作回顧
    氣象基層臺站建設
    西藏科技(2021年12期)2022-01-17 08:46:38
    數(shù)值大小比較“招招鮮”
    基于深度學習的海溫觀測數(shù)據(jù)質量控制應用研究
    海洋通報(2021年3期)2021-08-14 02:20:48
    基于Argo、XBT數(shù)據(jù)的蘇拉威西海溫鹽特征分析
    海洋通報(2021年1期)2021-07-23 01:55:22
    南印度洋偶極型海溫與中國西南地區(qū)初秋降水的關系
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    基層臺站綜合觀測業(yè)務管理之我見
    西藏科技(2015年6期)2015-09-26 12:12:13
    波浪輸運通量剩余量對上層海溫的影響分析
    黄色欧美视频在线观看| 日本与韩国留学比较| 直男gayav资源| www.www免费av| 国产精品女同一区二区软件 | 成人性生交大片免费视频hd| 少妇丰满av| av天堂中文字幕网| 看免费成人av毛片| 久久久国产成人免费| 亚洲av成人av| 免费电影在线观看免费观看| 精品人妻视频免费看| 国产午夜精品论理片| 中文字幕精品亚洲无线码一区| 免费看美女性在线毛片视频| 亚洲国产欧美人成| АⅤ资源中文在线天堂| netflix在线观看网站| 午夜精品一区二区三区免费看| 美女高潮的动态| 欧美另类亚洲清纯唯美| 精品日产1卡2卡| av视频在线观看入口| 色噜噜av男人的天堂激情| 成人国产综合亚洲| 亚洲午夜理论影院| 男人狂女人下面高潮的视频| 中出人妻视频一区二区| 午夜激情欧美在线| 欧美成人性av电影在线观看| 天美传媒精品一区二区| 亚洲av日韩精品久久久久久密| 综合色av麻豆| 亚洲内射少妇av| 国产精品免费一区二区三区在线| 最好的美女福利视频网| 亚洲第一电影网av| 日韩 亚洲 欧美在线| 干丝袜人妻中文字幕| 久久久国产成人免费| 男人狂女人下面高潮的视频| 婷婷丁香在线五月| 如何舔出高潮| 精品国产三级普通话版| 国产蜜桃级精品一区二区三区| 国产高清激情床上av| 一进一出好大好爽视频| 成人国产一区最新在线观看| 少妇丰满av| 波多野结衣高清作品| 欧美精品啪啪一区二区三区| 久久国产精品人妻蜜桃| 欧美日韩亚洲国产一区二区在线观看| 久久久久久伊人网av| 国产精品福利在线免费观看| 99热这里只有精品一区| av中文乱码字幕在线| 搡女人真爽免费视频火全软件 | 草草在线视频免费看| 精品国产三级普通话版| 国产男靠女视频免费网站| 欧美日韩国产亚洲二区| 特级一级黄色大片| 少妇的逼水好多| 99久久精品一区二区三区| 成人高潮视频无遮挡免费网站| 黄色配什么色好看| 国产视频一区二区在线看| 久久久久久久久中文| 啦啦啦观看免费观看视频高清| 欧美一级a爱片免费观看看| 久久久色成人| 麻豆精品久久久久久蜜桃| 桃红色精品国产亚洲av| 久久九九热精品免费| 少妇人妻精品综合一区二区 | 免费大片18禁| 国产精品久久久久久亚洲av鲁大| 欧美精品啪啪一区二区三区| 国产一级毛片七仙女欲春2| 国产精品亚洲一级av第二区| 久久久久久九九精品二区国产| 国产一区二区三区视频了| 一进一出抽搐gif免费好疼| 中文字幕久久专区| 国产色爽女视频免费观看| 亚洲成av人片在线播放无| 午夜激情福利司机影院| 男人狂女人下面高潮的视频| 十八禁国产超污无遮挡网站| 日本熟妇午夜| 午夜福利在线观看免费完整高清在 | 免费观看精品视频网站| 黄色欧美视频在线观看| 小蜜桃在线观看免费完整版高清| 韩国av一区二区三区四区| 人妻久久中文字幕网| 老司机午夜福利在线观看视频| 免费看av在线观看网站| 国产av一区在线观看免费| 久久国产精品人妻蜜桃| 如何舔出高潮| 淫妇啪啪啪对白视频| 免费观看精品视频网站| 欧美绝顶高潮抽搐喷水| 91久久精品电影网| 国内揄拍国产精品人妻在线| 99视频精品全部免费 在线| 日韩人妻高清精品专区| 国产三级在线视频| 国产主播在线观看一区二区| 国产三级中文精品| 国产一区二区三区视频了| 亚洲精品一区av在线观看| 日本色播在线视频| 国产男靠女视频免费网站| 色综合亚洲欧美另类图片| 老熟妇仑乱视频hdxx| 男女啪啪激烈高潮av片| 午夜福利18| 亚洲成人精品中文字幕电影| 国产高清激情床上av| 久久亚洲精品不卡| 亚洲精品粉嫩美女一区| 搡老妇女老女人老熟妇| 黄色视频,在线免费观看| 亚洲国产精品合色在线| 九九爱精品视频在线观看| 真实男女啪啪啪动态图| 亚洲久久久久久中文字幕| 日韩大尺度精品在线看网址| 午夜日韩欧美国产| 尾随美女入室| 午夜免费男女啪啪视频观看 | 九九爱精品视频在线观看| 美女被艹到高潮喷水动态| 伦精品一区二区三区| 国产在线男女| 小蜜桃在线观看免费完整版高清| 亚洲在线观看片| 国产精品一区二区三区四区免费观看 | 欧美成人一区二区免费高清观看| 99在线视频只有这里精品首页| 国产精华一区二区三区| 少妇人妻精品综合一区二区 | 亚洲四区av| 校园春色视频在线观看| 99热这里只有精品一区| 少妇的逼水好多| 少妇高潮的动态图| 1000部很黄的大片| 97人妻精品一区二区三区麻豆| 亚洲国产精品sss在线观看| 亚洲真实伦在线观看| 村上凉子中文字幕在线| 亚洲一级一片aⅴ在线观看| 变态另类丝袜制服| 俄罗斯特黄特色一大片| 国产亚洲av嫩草精品影院| 好男人在线观看高清免费视频| 无人区码免费观看不卡| 久久久久久久亚洲中文字幕| 99久久成人亚洲精品观看| 亚洲美女搞黄在线观看 | 嫩草影院入口| 最后的刺客免费高清国语| 狂野欧美白嫩少妇大欣赏| 啪啪无遮挡十八禁网站| 99久久无色码亚洲精品果冻| 亚洲aⅴ乱码一区二区在线播放| 日韩亚洲欧美综合| 成熟少妇高潮喷水视频| 中文资源天堂在线| 欧美黑人欧美精品刺激| 欧美日韩亚洲国产一区二区在线观看| 91午夜精品亚洲一区二区三区 | 国产在线精品亚洲第一网站| 国产精品人妻久久久久久| 久久久久久久久久成人| 麻豆成人av在线观看| 一进一出抽搐动态| 18禁黄网站禁片午夜丰满| 精品无人区乱码1区二区| 中出人妻视频一区二区| 啦啦啦韩国在线观看视频| 国产乱人伦免费视频| 日本精品一区二区三区蜜桃| 亚洲av电影不卡..在线观看| 嫩草影视91久久| 最近中文字幕高清免费大全6 | 伊人久久精品亚洲午夜| 免费在线观看影片大全网站| 国产真实伦视频高清在线观看 | 亚洲人与动物交配视频| 乱码一卡2卡4卡精品| 国产乱人伦免费视频| 床上黄色一级片| 91精品国产九色| 少妇熟女aⅴ在线视频| 欧美一区二区精品小视频在线| 国产一区二区在线观看日韩| 免费不卡的大黄色大毛片视频在线观看 | 国内揄拍国产精品人妻在线| 亚州av有码| 美女 人体艺术 gogo| 伦理电影大哥的女人| АⅤ资源中文在线天堂| 黄色丝袜av网址大全| 亚洲精品亚洲一区二区| 观看美女的网站| 精品人妻熟女av久视频| 精品日产1卡2卡| 99久久精品一区二区三区| 日韩,欧美,国产一区二区三区 | 最近中文字幕高清免费大全6 | 久久婷婷人人爽人人干人人爱| 99国产极品粉嫩在线观看| 老熟妇仑乱视频hdxx| 丰满人妻一区二区三区视频av| 丝袜美腿在线中文| 日本在线视频免费播放| 日韩精品中文字幕看吧| 久久亚洲精品不卡| 亚洲精品国产成人久久av| 偷拍熟女少妇极品色| 欧美一区二区亚洲| 99久久精品热视频| 亚洲精品在线观看二区| 国内精品美女久久久久久| 亚洲最大成人中文| 亚洲,欧美,日韩| 小蜜桃在线观看免费完整版高清| 联通29元200g的流量卡| 免费av不卡在线播放| www日本黄色视频网| 国内毛片毛片毛片毛片毛片| 能在线免费观看的黄片| 久久99热6这里只有精品| 精品久久久噜噜| 日日摸夜夜添夜夜添av毛片 | 我的女老师完整版在线观看| 国国产精品蜜臀av免费| 久久午夜福利片| 一级av片app| 联通29元200g的流量卡| 夜夜看夜夜爽夜夜摸| 十八禁网站免费在线| 91麻豆av在线| 欧美日本亚洲视频在线播放| 一进一出抽搐动态| 日本黄大片高清| 听说在线观看完整版免费高清| 91精品国产九色| 91狼人影院| 人人妻,人人澡人人爽秒播| 免费高清视频大片| 日韩高清综合在线| 如何舔出高潮| 欧美激情久久久久久爽电影| 97碰自拍视频| 天堂√8在线中文| 午夜老司机福利剧场| 国产精品久久电影中文字幕| 两个人的视频大全免费| 亚洲av第一区精品v没综合| 九九爱精品视频在线观看| 老熟妇乱子伦视频在线观看| 日本色播在线视频| 深夜精品福利| 久久久久久久精品吃奶| 日本黄色片子视频| 精品久久久久久久久久免费视频| 一级黄色大片毛片| 一个人免费在线观看电影| 99精品在免费线老司机午夜| 久久中文看片网| 久久久国产成人精品二区| 国产黄片美女视频| 日韩欧美精品v在线| 超碰av人人做人人爽久久| 午夜激情欧美在线| 99热这里只有是精品50| 白带黄色成豆腐渣| 久久这里只有精品中国| 国产黄色小视频在线观看| 国产成人av教育| 黄色丝袜av网址大全| 成人特级黄色片久久久久久久| 舔av片在线| 国产精品一区二区免费欧美| 香蕉av资源在线| 最近最新免费中文字幕在线| 韩国av在线不卡| 中文字幕免费在线视频6| 午夜影院日韩av| 亚洲av中文av极速乱 | 国产精品人妻久久久影院| 嫩草影院精品99| 亚洲精品粉嫩美女一区| 亚洲七黄色美女视频| 18禁裸乳无遮挡免费网站照片| 成人二区视频| 日本熟妇午夜| 波多野结衣高清作品| 国内精品宾馆在线| 日韩欧美在线二视频| 亚洲中文日韩欧美视频| 国产高清不卡午夜福利| 午夜福利欧美成人| av在线天堂中文字幕| 99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 午夜免费成人在线视频| 国产精品一区二区三区四区久久| 日本一二三区视频观看| 最近最新中文字幕大全电影3| 国产熟女欧美一区二区| 免费av毛片视频| 春色校园在线视频观看| 观看免费一级毛片| 国产一级毛片七仙女欲春2| 国产一区二区在线av高清观看| 精品午夜福利视频在线观看一区| 欧美又色又爽又黄视频| 精品人妻熟女av久视频| 色哟哟哟哟哟哟| 日本免费a在线| 麻豆精品久久久久久蜜桃| 亚洲性久久影院| 国产午夜精品久久久久久一区二区三区 | 老司机午夜福利在线观看视频| 搞女人的毛片| 国产精品久久久久久精品电影| 午夜老司机福利剧场| 日韩欧美在线二视频| 99久久九九国产精品国产免费| 99久久中文字幕三级久久日本| 久99久视频精品免费| 网址你懂的国产日韩在线| 国产黄色小视频在线观看| 色5月婷婷丁香| 亚洲狠狠婷婷综合久久图片| 国产在视频线在精品| 精品一区二区三区人妻视频| 色综合站精品国产| 乱系列少妇在线播放| 可以在线观看毛片的网站| 又紧又爽又黄一区二区| 国产精品国产三级国产av玫瑰| 少妇猛男粗大的猛烈进出视频 | 国产色婷婷99| 麻豆成人av在线观看| 九九热线精品视视频播放| 日韩 亚洲 欧美在线| 久久久久久久久久成人| 精品久久久久久久久亚洲 | 国产老妇女一区| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| av天堂中文字幕网| 美女黄网站色视频| 能在线免费观看的黄片| 免费一级毛片在线播放高清视频| 一区二区三区高清视频在线| 极品教师在线免费播放| 精品午夜福利视频在线观看一区| 免费观看在线日韩| a级毛片a级免费在线| 日本爱情动作片www.在线观看 | 97超级碰碰碰精品色视频在线观看| 老熟妇仑乱视频hdxx| 午夜精品久久久久久毛片777| 精品一区二区三区视频在线观看免费| 国产高清有码在线观看视频| 国产精品野战在线观看| 亚洲自拍偷在线| 国产亚洲精品久久久com| 欧美日本视频| 成人三级黄色视频| 午夜亚洲福利在线播放| 级片在线观看| 欧美又色又爽又黄视频| 少妇裸体淫交视频免费看高清| 人妻久久中文字幕网| av.在线天堂| 最近视频中文字幕2019在线8| 亚洲久久久久久中文字幕| 久久国产乱子免费精品| 国产视频内射| 在线观看av片永久免费下载| 日韩高清综合在线| 国产精品人妻久久久久久| 男女做爰动态图高潮gif福利片| 亚洲欧美清纯卡通| 亚洲国产日韩欧美精品在线观看| 日日撸夜夜添| 成人二区视频| 国产老妇女一区| 十八禁网站免费在线| 一进一出抽搐动态| 亚洲最大成人中文| 久久久久久久久久黄片| 天天躁日日操中文字幕| 一本精品99久久精品77| 又粗又爽又猛毛片免费看| 国产伦人伦偷精品视频| 亚洲人成网站高清观看| 精品国产三级普通话版| 日韩亚洲欧美综合| 国产精品乱码一区二三区的特点| www.色视频.com| 99热6这里只有精品| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 成人美女网站在线观看视频| 少妇丰满av| 一进一出抽搐动态| av天堂中文字幕网| 一个人免费在线观看电影| 啦啦啦韩国在线观看视频| 国内揄拍国产精品人妻在线| 国产成人aa在线观看| 天天一区二区日本电影三级| 婷婷六月久久综合丁香| 非洲黑人性xxxx精品又粗又长| 中文亚洲av片在线观看爽| 色综合站精品国产| 久久精品综合一区二区三区| 亚洲av电影不卡..在线观看| 日韩精品青青久久久久久| 高清毛片免费观看视频网站| 午夜日韩欧美国产| 久久婷婷人人爽人人干人人爱| 成年人黄色毛片网站| 免费无遮挡裸体视频| 国产成人福利小说| 日韩欧美在线二视频| 亚洲真实伦在线观看| 亚洲黑人精品在线| 波多野结衣巨乳人妻| 九九久久精品国产亚洲av麻豆| 色视频www国产| 人妻夜夜爽99麻豆av| 久久亚洲真实| 久久久久久久久久成人| 国产一区二区激情短视频| 淫妇啪啪啪对白视频| 日韩中字成人| 一进一出抽搐动态| 免费无遮挡裸体视频| 国产精品女同一区二区软件 | 我的老师免费观看完整版| 亚洲男人的天堂狠狠| 99久久九九国产精品国产免费| 国产一区二区在线av高清观看| 久久精品国产亚洲av天美| 日韩欧美国产一区二区入口| 给我免费播放毛片高清在线观看| 欧美人与善性xxx| 波多野结衣高清作品| 亚洲精品乱码久久久v下载方式| 欧美日本亚洲视频在线播放| 亚洲国产色片| 91狼人影院| 精品免费久久久久久久清纯| 中文资源天堂在线| 亚洲精品粉嫩美女一区| 国内久久婷婷六月综合欲色啪| 美女大奶头视频| 日韩av在线大香蕉| 久久久久免费精品人妻一区二区| 97超视频在线观看视频| 亚洲国产精品久久男人天堂| 免费看a级黄色片| 色精品久久人妻99蜜桃| netflix在线观看网站| 国产精品一区二区三区四区久久| 十八禁国产超污无遮挡网站| 97人妻精品一区二区三区麻豆| 精品无人区乱码1区二区| 91久久精品电影网| 免费搜索国产男女视频| 欧美潮喷喷水| 久久精品国产清高在天天线| 亚洲欧美精品综合久久99| 村上凉子中文字幕在线| 99热6这里只有精品| 国产精品久久久久久久电影| а√天堂www在线а√下载| 日日夜夜操网爽| 日韩欧美一区二区三区在线观看| 夜夜夜夜夜久久久久| 精品一区二区三区人妻视频| a在线观看视频网站| 国产高清不卡午夜福利| 国产欧美日韩一区二区精品| 亚洲国产精品sss在线观看| 中文在线观看免费www的网站| 女生性感内裤真人,穿戴方法视频| 免费看光身美女| 欧美xxxx黑人xx丫x性爽| 在线观看免费视频日本深夜| 99国产精品一区二区蜜桃av| 欧美日韩国产亚洲二区| 麻豆国产97在线/欧美| 级片在线观看| 免费av毛片视频| 亚洲国产精品久久男人天堂| 日本五十路高清| 亚洲四区av| 国语自产精品视频在线第100页| 久久6这里有精品| 联通29元200g的流量卡| 日韩欧美一区二区三区在线观看| 一级av片app| 日韩精品青青久久久久久| 久99久视频精品免费| 亚洲一区二区三区色噜噜| 成人无遮挡网站| 免费在线观看影片大全网站| 精品日产1卡2卡| 久久这里只有精品中国| 在线播放国产精品三级| 一进一出好大好爽视频| 成人二区视频| 国产乱人视频| 99久久精品国产国产毛片| 色综合色国产| 欧美在线一区亚洲| 国产亚洲精品综合一区在线观看| 亚洲国产欧美人成| 免费观看精品视频网站| 久久精品国产亚洲av香蕉五月| 高清在线国产一区| 99热网站在线观看| av在线老鸭窝| 人人妻,人人澡人人爽秒播| 日本三级黄在线观看| 亚洲av日韩精品久久久久久密| 午夜视频国产福利| 18禁在线播放成人免费| 日韩,欧美,国产一区二区三区 | 草草在线视频免费看| 亚洲美女搞黄在线观看 | 国产免费一级a男人的天堂| 免费搜索国产男女视频| 成人高潮视频无遮挡免费网站| 最近中文字幕高清免费大全6 | 欧美三级亚洲精品| 国产欧美日韩一区二区精品| 亚洲 国产 在线| 国产视频内射| 亚洲精品在线观看二区| 国产真实乱freesex| 动漫黄色视频在线观看| 欧美黑人巨大hd| av在线亚洲专区| 亚洲av熟女| 亚洲av二区三区四区| 色精品久久人妻99蜜桃| av天堂中文字幕网| 俺也久久电影网| 无遮挡黄片免费观看| 亚洲人成网站在线播| 在线免费观看的www视频| 国内精品久久久久精免费| АⅤ资源中文在线天堂| 国产av麻豆久久久久久久| 亚洲精品色激情综合| 国产淫片久久久久久久久| netflix在线观看网站| 久久久久久久久大av| av在线老鸭窝| 日本精品一区二区三区蜜桃| 啦啦啦韩国在线观看视频| 久久久久久大精品| 亚洲最大成人中文| 久久6这里有精品| 久久久久国产精品人妻aⅴ院| 久久亚洲真实| 在线观看av片永久免费下载| 国产在视频线在精品| 日本免费a在线| 一区二区三区高清视频在线| 亚洲第一电影网av| 观看美女的网站| 国产爱豆传媒在线观看| 尤物成人国产欧美一区二区三区| 国产精品精品国产色婷婷| 国产精品98久久久久久宅男小说| 亚洲最大成人手机在线| 国内少妇人妻偷人精品xxx网站| 99热6这里只有精品| 免费在线观看影片大全网站| 国产精品综合久久久久久久免费| 亚洲人成伊人成综合网2020| 精华霜和精华液先用哪个| 国产精品永久免费网站| 九九久久精品国产亚洲av麻豆| 国产在视频线在精品| 亚洲内射少妇av| 亚洲中文日韩欧美视频| 少妇裸体淫交视频免费看高清| 欧美激情久久久久久爽电影| 一区二区三区免费毛片| 成人特级黄色片久久久久久久| 可以在线观看的亚洲视频| 亚洲成人精品中文字幕电影| 中文字幕免费在线视频6| 国产精品,欧美在线| 日韩欧美在线乱码| 中亚洲国语对白在线视频| 精品久久久久久,| 看片在线看免费视频|