本發(fā)明涉環(huán)境監(jiān)測(cè)領(lǐng)域，尤其涉及一種赤潮生物量預(yù)測(cè)方法及裝置。

背景技術(shù)：

赤潮是一種由于局部海區(qū)的浮游生物突發(fā)性急劇繁殖，并聚集在一起而引起海平面顏色出現(xiàn)異常和發(fā)臭的現(xiàn)象。赤潮的發(fā)生是生物、化學(xué)、水文、氣象等因素綜合影響的結(jié)果。赤潮生物的生長(zhǎng)、繁殖、消亡過(guò)程與環(huán)境因子之間存在著復(fù)雜的非線性關(guān)系，且這些環(huán)境因子具有時(shí)間上的連續(xù)性以及空間上的異質(zhì)性。近年來(lái)近岸海域赤潮頻發(fā)，赤潮的爆發(fā)不僅會(huì)導(dǎo)致水環(huán)境生態(tài)系統(tǒng)失衡，還會(huì)造成水產(chǎn)養(yǎng)殖災(zāi)難。

傳統(tǒng)的赤潮預(yù)測(cè)方法是通過(guò)各個(gè)環(huán)境因子得到當(dāng)前的赤潮狀況，但是當(dāng)赤潮的爆發(fā)時(shí)已經(jīng)造成水環(huán)境生態(tài)系統(tǒng)失衡和水產(chǎn)養(yǎng)殖災(zāi)難，因此無(wú)法提前對(duì)赤潮進(jìn)行預(yù)知。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明實(shí)施例通過(guò)提供一種赤潮生物量預(yù)測(cè)方法及裝置，解決了無(wú)法提前對(duì)赤潮進(jìn)行預(yù)知的技術(shù)問(wèn)題。

第一方面，本發(fā)明實(shí)施例提供了一種赤潮生物量預(yù)測(cè)方法，包括：

獲得M組有效海洋環(huán)境因子，其中，M為所述有效海洋環(huán)境因子的種類(lèi)數(shù)乘以子海域的個(gè)數(shù)，同一組所述有效海洋環(huán)境因子包括多幀在不同時(shí)間點(diǎn)采集的同一種類(lèi)的有效海洋環(huán)境因子；

對(duì)所述M組有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理，以得到對(duì)應(yīng)的M組歸一化后海洋環(huán)境因子；

分別將所述M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至對(duì)應(yīng)的ARIMA模型中進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出M組海洋環(huán)境因子預(yù)測(cè)值，m依次為1至M，其中，一個(gè)ARIMA模型對(duì)應(yīng)一組歸一化后海洋環(huán)境因子；

將所述M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值。

優(yōu)選的，所述獲得M組有效海洋環(huán)境因，包括：

建立原始海洋環(huán)境因子的關(guān)系數(shù)據(jù)庫(kù)，其中，在所述關(guān)系數(shù)據(jù)庫(kù)中保存有多個(gè)歷史監(jiān)測(cè)的原始海洋環(huán)境因子，以及保存有各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)位置信息和監(jiān)測(cè)時(shí)間信息；

根據(jù)各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)時(shí)間信息和監(jiān)測(cè)位置信息，從所述關(guān)系數(shù)據(jù)庫(kù)中提取M組原始海洋環(huán)境因子，其中，所述M組原始海洋環(huán)境因子中的同一組原始海洋環(huán)境因子在同一子海域，且距離當(dāng)前時(shí)刻小于預(yù)設(shè)時(shí)間間隔監(jiān)測(cè)到；

從所述M組原始海洋環(huán)境因子中的每組原始海洋環(huán)境因子中篩選出滿足預(yù)設(shè)條件的因子組成為所述M組有效海洋環(huán)境因子。

優(yōu)選的，在所述分別將所述M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至對(duì)應(yīng)的ARIMA模型中進(jìn)行預(yù)測(cè)之前，所述方法還包括：通過(guò)如下步驟建立所述ARIMA模型：

根據(jù)所述關(guān)系數(shù)據(jù)庫(kù)中各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)位置信息、監(jiān)測(cè)時(shí)間信息確定針對(duì)所述M組原始海洋環(huán)境因子中每組原始海洋環(huán)境因子的ARIMA模型；

基于AIC準(zhǔn)則分別針對(duì)所述每組原始海洋環(huán)境因子的ARIMA模型進(jìn)行參數(shù)估計(jì)；

選取基于所述AIC準(zhǔn)則進(jìn)行參數(shù)估計(jì)得到最小AIC值時(shí)的參數(shù)組合，作為針對(duì)所述每組原始海洋環(huán)境因子的ARIMA模型的模型參數(shù)。

優(yōu)選的，在所述將所述M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)之前，所述方法還包括：通過(guò)如下步驟建立所述BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型：

獲取赤潮生物量測(cè)量值樣本；

基于如下對(duì)數(shù)公式對(duì)所述赤潮生物量測(cè)量值樣本進(jìn)行歸一化處理得到對(duì)應(yīng)的歸一化后赤潮生物量樣本：

其中，y′為歸一化后赤潮生物量的值，y為赤潮生物量測(cè)量值，y_max為所述赤潮生物量測(cè)量值樣本的最大值，y_min為所述赤潮生物量測(cè)量值樣本的最小值；

基于所述歸一化后赤潮生物量樣本確定BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)以及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)；

基于確定的BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)建立針對(duì)所述M組海洋環(huán)境因子預(yù)測(cè)值的BP神經(jīng)網(wǎng)絡(luò)的最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)；

將所述最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)經(jīng)過(guò)學(xué)習(xí)過(guò)程得到的最優(yōu)值作為所述BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

優(yōu)選的，在所述預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值之后，所述方法還包括：

將所述赤潮生物量預(yù)測(cè)值與對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值進(jìn)行對(duì)比；

輸出所述赤潮生物量預(yù)測(cè)值與所述對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值之間的均方根誤差、平均絕對(duì)誤差、平均絕對(duì)百分誤差、擬合度中的至少一種對(duì)比結(jié)果。

優(yōu)選的，所述有效海洋環(huán)境因子的種類(lèi)包括：酸堿度、溶解氧、水溫、溶解氧飽和度、葉綠素-a、磷酸鹽、氨氮、亞硝酸鹽氮、硝酸鹽氮、鹽度、化學(xué)耗氧量和硅酸鹽。

第二方面，本發(fā)明實(shí)施例提供了一種赤潮生物量預(yù)測(cè)裝置，包括：

獲得單元，用于獲得M組有效海洋環(huán)境因子，其中，M為所述有效海洋環(huán)境因子的種類(lèi)數(shù)乘以子海域的個(gè)數(shù)，同一組所述有效海洋環(huán)境因子包括多幀在不同時(shí)間點(diǎn)采集的同一種類(lèi)的有效海洋環(huán)境因子；

預(yù)處理單元，用于對(duì)所述M組有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理，以得到對(duì)應(yīng)的M組歸一化后海洋環(huán)境因子；

環(huán)境因子預(yù)測(cè)單元，用于分別將所述M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至對(duì)應(yīng)的ARIMA模型中進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出M組海洋環(huán)境因子預(yù)測(cè)值，m依次為1至M，其中，一個(gè)ARIMA模型對(duì)應(yīng)一組歸一化后海洋環(huán)境因子；

赤潮生物量預(yù)測(cè)單元，用于將所述M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值。

優(yōu)選的，所述獲得單元包括：

數(shù)據(jù)庫(kù)建立子單元，用于建立原始海洋環(huán)境因子的關(guān)系數(shù)據(jù)庫(kù)，其中，在所述關(guān)系數(shù)據(jù)庫(kù)中保存有多個(gè)歷史監(jiān)測(cè)的原始海洋環(huán)境因子，以及保存有各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)位置信息和監(jiān)測(cè)時(shí)間信息；

提取子單元，用于根據(jù)各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)時(shí)間信息和監(jiān)測(cè)位置信息，從所述關(guān)系數(shù)據(jù)庫(kù)中提取M組原始海洋環(huán)境因子，其中，所述M組原始海洋環(huán)境因子中的同一組原始海洋環(huán)境因子在同一子海域，且距離當(dāng)前時(shí)刻小于預(yù)設(shè)時(shí)間間隔監(jiān)測(cè)到；

篩選子單元，用于從所述M組原始海洋環(huán)境因子中的每組原始海洋環(huán)境因子中篩選出滿足預(yù)設(shè)條件的因子組成為所述M組有效海洋環(huán)境因子。

優(yōu)選的，所述赤潮生物量預(yù)測(cè)裝置還包括：第一建模單元，用于通過(guò)如下步驟建立所述ARIMA模型：

用于根據(jù)所述關(guān)系數(shù)據(jù)庫(kù)中各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)位置信息、監(jiān)測(cè)時(shí)間信息確定針對(duì)所述M組原始海洋環(huán)境因子中每組原始海洋環(huán)境因子的ARIMA模型；

基于AIC準(zhǔn)則分別針對(duì)所述每組原始海洋環(huán)境因子的ARIMA模型進(jìn)行參數(shù)估計(jì)；

選取基于所述AIC準(zhǔn)則進(jìn)行參數(shù)估計(jì)得到最小AIC值時(shí)的參數(shù)組合，作為針對(duì)所述每組原始海洋環(huán)境因子的ARIMA模型的模型參數(shù)。

優(yōu)選的，所述赤潮生物量預(yù)測(cè)裝置還包括：第二建模單元，用于通過(guò)如下步驟建立所述BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型：

獲取赤潮生物量測(cè)量值樣本；

基于如下對(duì)數(shù)公式對(duì)所述赤潮生物量測(cè)量值樣本進(jìn)行歸一化處理得到對(duì)應(yīng)的歸一化后赤潮生物量樣本：

其中，y′為歸一化后赤潮生物量的值，y為赤潮生物量測(cè)量值，y_max為所述赤潮生物量測(cè)量值樣本的最大值，y_min為所述赤潮生物量測(cè)量值樣本的最小值；

基于所述歸一化后赤潮生物量樣本確定BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)以及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)；

基于確定的BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)建立針對(duì)所述M組海洋環(huán)境因子預(yù)測(cè)值的BP神經(jīng)網(wǎng)絡(luò)的最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)；

將所述最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)經(jīng)過(guò)學(xué)習(xí)過(guò)程得到的最優(yōu)值作為所述BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

優(yōu)選的，所述赤潮生物量預(yù)測(cè)裝置還包括：對(duì)比單元，用于將所述赤潮生物量預(yù)測(cè)值與對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值進(jìn)行對(duì)比；結(jié)果輸出單元，用于輸出所述赤潮生物量預(yù)測(cè)值與所述對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值之間的均方根誤差、平均絕對(duì)誤差、平均絕對(duì)百分誤差、擬合度中的至少一種對(duì)比結(jié)果。

優(yōu)選的，所述有效海洋環(huán)境因子的種類(lèi)包括：酸堿度、溶解氧、水溫、溶解氧飽和度、葉綠素-a、磷酸鹽、氨氮、亞硝酸鹽氮、硝酸鹽氮、鹽度、化學(xué)耗氧量和硅酸鹽。

通過(guò)上述本發(fā)明實(shí)施例提供的一個(gè)或多個(gè)技術(shù)方案，至少實(shí)現(xiàn)了如下技術(shù)效果或優(yōu)點(diǎn)：

通過(guò)獲得M組有效海洋環(huán)境因子；對(duì)M組有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理，以得到對(duì)應(yīng)的M組歸一化后海洋環(huán)境因子；分別將M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至ARIMA模型中進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出M組海洋環(huán)境因子預(yù)測(cè)值；將M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值。其中，M為所述有效海洋環(huán)境因子的種類(lèi)數(shù)乘以子海域的個(gè)數(shù)，一組有效海洋環(huán)境因子包括多幀在不同時(shí)間點(diǎn)采集的同一種類(lèi)的有效海洋環(huán)境因子；從而充分考慮了影響赤潮的海洋環(huán)境因子的橫向時(shí)序相關(guān)性及縱向空間異質(zhì)性，能夠針對(duì)不同子海域的不同因子建立不同的ARIMA模型來(lái)描述其時(shí)間上的連續(xù)性及空間異質(zhì)性，從而預(yù)測(cè)出未來(lái)的各海洋環(huán)境因子，再此基礎(chǔ)上利用BP神經(jīng)網(wǎng)絡(luò)來(lái)表達(dá)各海洋環(huán)境因子與赤潮之間的復(fù)雜關(guān)系，進(jìn)而實(shí)現(xiàn)赤潮生物量的短期準(zhǔn)確預(yù)測(cè)，以預(yù)報(bào)未來(lái)幾天是否發(fā)生赤潮，以便進(jìn)行防備。

附圖說(shuō)明

為了更清楚地說(shuō)明本發(fā)明實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本發(fā)明的實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)的前提下，還可以根據(jù)提供的附圖獲得其他的附圖。

圖1為本發(fā)明實(shí)施例中赤潮生物量預(yù)測(cè)方法的流程示意圖；

圖2為本發(fā)明實(shí)施例中赤潮生物量預(yù)測(cè)方法的功能單元圖。

具體實(shí)施方式

為使本發(fā)明實(shí)施例的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚，下面將結(jié)合本發(fā)明實(shí)施例中的附圖，對(duì)本發(fā)明實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例是本發(fā)明一部分實(shí)施例，而不是全部的實(shí)施例?；诒景l(fā)明中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有作出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本發(fā)明保護(hù)的范圍。

參考圖1所示，本發(fā)明實(shí)施例提供了一種赤潮生物量預(yù)測(cè)方法，包括如下步驟：

S101、獲得M組有效海洋環(huán)境因子，其中，M為有效海洋環(huán)境因子的種類(lèi)數(shù)乘以子海域的個(gè)數(shù)，同一組有效海洋環(huán)境因子包括多幀在不同時(shí)間點(diǎn)采集的同一種類(lèi)的有效海洋環(huán)境因子。

在一具體實(shí)施例中，有效海洋環(huán)境因子的種類(lèi)包括：酸堿度、溶解氧、水溫、溶解氧飽和度、葉綠素-a、磷酸鹽、氨氮、亞硝酸鹽氮、硝酸鹽氮、鹽度、化學(xué)耗氧量和硅酸鹽這12類(lèi)。需要進(jìn)行赤潮生物量預(yù)測(cè)的海域劃分為多個(gè)子海域，每個(gè)子海域?qū)?yīng)一個(gè)經(jīng)緯度范圍，比如，劃分為6個(gè)子海域，則獲得12*6組有效海洋環(huán)境因子。每個(gè)子海域?qū)?yīng)有12組有效海洋環(huán)境因子：一組酸堿度因子、一組溶解氧因子、一組溶解氧因子、一組水溫因子、一組溶解氧飽和度因子、一組葉綠素-a因子、一組磷酸鹽因子、一組氨氮因子、一組亞硝酸鹽氮因子、一組鹽度因子、一組化學(xué)耗氧量因子、一組硅酸鹽因子。

為了獲得M組有效海洋環(huán)境因子，則在執(zhí)行S101之前，還包括如下步驟：

首先，建立原始海洋環(huán)境因子的關(guān)系數(shù)據(jù)庫(kù)，其中，在關(guān)系數(shù)據(jù)庫(kù)中保存有多個(gè)歷史監(jiān)測(cè)的原始海洋環(huán)境因子，以及保存有各個(gè)原始海洋環(huán)境因子的監(jiān)測(cè)位置信息和監(jiān)測(cè)時(shí)間信息。

具體的，在關(guān)系數(shù)據(jù)庫(kù)中保存的原始海洋環(huán)境因子是設(shè)置在不同子海域的監(jiān)測(cè)設(shè)備在不同時(shí)間采集得到。從而形成了原始海洋環(huán)境因子的大數(shù)據(jù)。具體的，監(jiān)測(cè)位置信息為用經(jīng)緯度坐標(biāo)表示的監(jiān)測(cè)設(shè)備所在的位置，監(jiān)測(cè)時(shí)間信息為采集原始海洋環(huán)境因子的時(shí)間點(diǎn)。

接著，根據(jù)各個(gè)原始海洋環(huán)境因子的監(jiān)測(cè)時(shí)間信息和監(jiān)測(cè)位置信息，從關(guān)系數(shù)據(jù)庫(kù)中提取M組原始海洋環(huán)境因子，其中，M組原始海洋環(huán)境因子中的同一組原始海洋環(huán)境因子在同一子海域，且距離當(dāng)前時(shí)刻小于預(yù)設(shè)時(shí)間間隔檢測(cè)到。比如，針對(duì)需要獲得36組有效海洋環(huán)境因子，每組有效海洋環(huán)境因子需要20幀，以各種原始海洋環(huán)境因子每天監(jiān)測(cè)一幀為例，從關(guān)系數(shù)據(jù)庫(kù)中各個(gè)子海域中的各種海洋環(huán)境因子都提取當(dāng)天～倒退19天的這20幀。

再接著，從M組原始海洋環(huán)境因子中的每組原始海洋環(huán)境因子中篩選出滿足預(yù)設(shè)條件的因子組成為M組有效海洋環(huán)境因子。

具體的，通過(guò)對(duì)M組原始海洋環(huán)境因子中滿足預(yù)設(shè)條件的因子進(jìn)行數(shù)字編碼來(lái)篩除監(jiān)測(cè)錯(cuò)誤或監(jiān)測(cè)空白的原始海洋環(huán)境因子，以保留得到M組有效海洋環(huán)境因子。

比如，由于海水是堿性的，針對(duì)酸堿度的預(yù)設(shè)條件為PH值大于或等于7，對(duì)酸堿度大于7的進(jìn)行數(shù)字編碼，對(duì)酸堿度小于7的不進(jìn)行數(shù)字編碼，從而篩除錯(cuò)誤的酸堿度。同理的，本領(lǐng)域技術(shù)人員可以設(shè)置其他類(lèi)原始海洋環(huán)境因子的預(yù)設(shè)條件，為了說(shuō)明書(shū)的簡(jiǎn)潔，本文不再贅述。

S102、對(duì)M組有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理，以得到對(duì)應(yīng)的M組歸一化后海洋環(huán)境因子。

具體的，通過(guò)線性方法對(duì)M組有效海洋環(huán)境因子中的每組有效海洋環(huán)境因子的每個(gè)有效海洋環(huán)境因子都進(jìn)行歸一化預(yù)處理，具體的，對(duì)每個(gè)有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理的公式如下：

其中，x′為歸一化后海洋環(huán)境因子，x為有效海洋環(huán)境因子，x_min為對(duì)應(yīng)的一類(lèi)海洋環(huán)境因子的最小值，x_max為對(duì)應(yīng)的一類(lèi)海洋環(huán)境因子的最大值。比如，x為有效的化學(xué)耗氧量，則此時(shí)，x_min為監(jiān)測(cè)到化學(xué)耗氧量的最小值，x_max為監(jiān)測(cè)的化學(xué)耗氧量的最大值。

S103、分別將M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至對(duì)應(yīng)的ARIMA模型(全稱(chēng)為差分自回歸移動(dòng)平均模型(Autoregressive Integrated Moving Average Model,簡(jiǎn)記為ARIMA)中進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出M組海洋環(huán)境因子預(yù)測(cè)值，m依次為1至M，其中，一個(gè)ARIMA模型對(duì)應(yīng)一組歸一化后海洋環(huán)境因子。

在具體實(shí)施過(guò)程中，需要在執(zhí)行S103之前，通過(guò)如下步驟建立分別針對(duì)M組歸一化后海洋環(huán)境因子的ARIMA模型，從而建立M個(gè)ARIMA模型：

首先，根據(jù)關(guān)系數(shù)據(jù)庫(kù)中各個(gè)原始海洋環(huán)境因子的監(jiān)測(cè)位置信息、監(jiān)測(cè)時(shí)間信息確定針對(duì)M組原始海洋環(huán)境因子中每組原始海洋環(huán)境因子的ARIMA模型。從而確定M個(gè)ARIMA模型。

具體的，確定的ARIMA模型為ARMA(自回歸移動(dòng)平均模型)模型，AR(自回歸)模型、MA(移動(dòng)平均)模型中的一種。具體來(lái)講，根據(jù)各個(gè)原始海洋環(huán)境因子的監(jiān)測(cè)位置信息，將各個(gè)原始海洋環(huán)境因子對(duì)應(yīng)至所在的子海域，再根據(jù)各個(gè)原始海洋環(huán)境因子的監(jiān)測(cè)時(shí)間信息，將對(duì)應(yīng)至所在子海域的12類(lèi)原始海洋環(huán)境因子的監(jiān)測(cè)測(cè)時(shí)間信息形成12個(gè)時(shí)間序列。從而形成“子海域個(gè)數(shù)*12”個(gè)時(shí)間序列。對(duì)形成的各時(shí)間序列進(jìn)行單位根檢驗(yàn)，以判斷各個(gè)時(shí)間序列是否為平穩(wěn)序列。平穩(wěn)序列為：均值和方差在時(shí)間過(guò)程上是常數(shù)，并且在任何兩時(shí)期的協(xié)方差值僅依賴(lài)該兩時(shí)期間的距離或滯后，否則為非平穩(wěn)序列。若進(jìn)行單位根檢驗(yàn)時(shí)存在單位根，則該時(shí)間序列為非平穩(wěn)序列，非平穩(wěn)序列需要進(jìn)行差分形成平穩(wěn)序列，差分次數(shù)記為d。

具體來(lái)講，對(duì)平穩(wěn)序列或差分后的平穩(wěn)序列通過(guò)相關(guān)性分析確定每類(lèi)原始海洋環(huán)境因子適用的ARIMA模型：具體來(lái)講，是對(duì)平穩(wěn)序列或差分后的平穩(wěn)序列進(jìn)行相關(guān)性分析，繪出自相關(guān)圖和偏相關(guān)圖，基于自相關(guān)圖和偏相關(guān)圖確定偏相關(guān)函數(shù)和自相關(guān)函數(shù)是截尾的還是拖尾的，若平穩(wěn)序列的偏相關(guān)函數(shù)是截尾的、自相關(guān)數(shù)是拖尾的，則該平穩(wěn)序列適用AR模型；若平穩(wěn)序列的偏相關(guān)函數(shù)拖尾、自相關(guān)函數(shù)截尾，則適用MR模型；若平穩(wěn)序列的偏相關(guān)函數(shù)、自相關(guān)函數(shù)均是拖尾的，則適用ARMA模型。其中，截尾是指時(shí)間序列的自相關(guān)函數(shù)(ACF)或偏自相關(guān)函數(shù)(PACF)在某階后均為0的性質(zhì)；拖尾是ACF或PACF并不在某階后均為0的性質(zhì)。

接著，基于AIC準(zhǔn)則分別針對(duì)每組原始海洋環(huán)境因子的ARIMA模型進(jìn)行參數(shù)估計(jì)，選取基于AIC準(zhǔn)則進(jìn)行參數(shù)估計(jì)得到最小AIC值時(shí)的參數(shù)組合，作為針對(duì)每組原始海洋環(huán)境因子的ARIMA模型的模型參數(shù)ARIMA(p,d,q)。

具體來(lái)講，根據(jù)AIC準(zhǔn)則確定ARMA模型的自回歸項(xiàng)數(shù)p與移動(dòng)平均項(xiàng)數(shù)q。根據(jù)AIC準(zhǔn)則確定AR模型的自回歸項(xiàng)p；根據(jù)AIC準(zhǔn)則確定MA模型的平均移動(dòng)項(xiàng)數(shù)q。AIC準(zhǔn)則的公式為AIC＝2K-2ln(F)，其中，K是項(xiàng)數(shù)，F(xiàn)是似然函數(shù)。

最后，設(shè)t,t-1,…t-p時(shí)刻的時(shí)間序列值分別為y_t,y_t-1,…,y_t-p，在t,t-1,…t-q時(shí)刻的誤差記作n_t,n_t-1,…n_t-q，n_t,n_t-1,…n_t-q是相互獨(dú)立且符合高斯分布的白噪聲序列。a₁,a₂,…,a_p是自回歸(AR)系數(shù)，b₁,b₂,…,b_q是滑動(dòng)平均(MA)系數(shù)，則ARMA模型的表達(dá)式如下所示：

y_t＝a₁y_t-1+a₂y_t-2+…+a_py_t-p+n_t+b₁n_t-1+b₂n_t-2+b_qn_t-q。

S104、將M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP(Back Propagation反向傳播)神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值。

在具體實(shí)施過(guò)程中，預(yù)設(shè)時(shí)期設(shè)置為小于7天，則預(yù)測(cè)出7天內(nèi)的赤潮生物量預(yù)測(cè)值，從而避免長(zhǎng)期預(yù)測(cè)的值不準(zhǔn)確。

進(jìn)一步的，在執(zhí)行S104之前，通過(guò)如下步驟建立針對(duì)海洋環(huán)境因子預(yù)測(cè)值進(jìn)行預(yù)測(cè)的BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型：

步驟一、獲取赤潮生物量測(cè)量值樣本；

具體的，赤潮生物量測(cè)量值樣本可以為設(shè)置在各個(gè)子海域中的監(jiān)測(cè)設(shè)備監(jiān)測(cè)得到的。

步驟二、基于如下對(duì)數(shù)公式對(duì)采集的赤潮生物量測(cè)量值樣本進(jìn)行歸一化處理得到對(duì)應(yīng)的歸一化后赤潮生物量樣本：其中，y′為歸一化后赤潮生物量的值，y為赤潮生物量測(cè)量值，y_max為赤潮生物量測(cè)量值樣本的最大值，y_min為赤潮生物量測(cè)量值樣本的最小值。

具體的，將赤潮生物量測(cè)量值樣本的每個(gè)赤潮生物量測(cè)量值分別經(jīng)過(guò)對(duì)數(shù)公式得到每個(gè)歸一化后赤潮生物量的值，全部歸一化后赤潮生物量構(gòu)成了以形成用于訓(xùn)練BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)的歸一化后赤潮生物量樣本。

步驟三、基于歸一化后赤潮生物量樣本確定需要的BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)的隱含層數(shù)以及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)。

具體的，將每一個(gè)歸一化后赤潮生物量分別帶入的BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)進(jìn)行實(shí)驗(yàn)，基于實(shí)驗(yàn)結(jié)果好壞確定所需要的隱含層數(shù)以及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)。

步驟四、基于確定的BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)建立BP神經(jīng)網(wǎng)絡(luò)的最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)；

步驟五、將最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)經(jīng)過(guò)學(xué)習(xí)過(guò)程得到最優(yōu)值作為BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

具體的，學(xué)習(xí)過(guò)程如下：

第一步、隨機(jī)初始化權(quán)重θ非零值，并通過(guò)前向傳播算法計(jì)算出損失函數(shù)J(θ)，并計(jì)算初始化權(quán)重θ的偏導(dǎo)D，計(jì)算公式如下：

第二步、給定學(xué)習(xí)過(guò)程的最大學(xué)習(xí)次數(shù)、設(shè)定誤差函數(shù)、預(yù)設(shè)精度值，給定輸入海洋環(huán)境因子樣本集：給定赤潮生物量輸出樣本集：

第三步、隨機(jī)從輸入海洋環(huán)境因子樣本集選取一個(gè)輸入樣本，對(duì)應(yīng)的，從赤潮生物量輸出樣本集中選取一個(gè)輸出樣本，使用反向傳播的梯度下降算法進(jìn)行計(jì)算得到全局誤差，判斷全局誤差是否達(dá)到預(yù)設(shè)精度值，或?qū)W習(xí)次數(shù)是否大于最大學(xué)習(xí)次數(shù)，當(dāng)全局誤差達(dá)到預(yù)設(shè)精度值或大于最大學(xué)習(xí)次數(shù)，則結(jié)束梯度下降算法，以得到使損失函數(shù)J(θ)最小的權(quán)重值θ為最優(yōu)值作為BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型，否則選取下一個(gè)輸入樣本及對(duì)應(yīng)的輸出樣本返回進(jìn)行下一次學(xué)習(xí)。得到的使損失函數(shù)J(θ)最小的權(quán)重值θ為最優(yōu)值。

在預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值之后，將赤潮生物量預(yù)測(cè)值與對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值進(jìn)行對(duì)比；輸出赤潮生物量預(yù)測(cè)值與對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值之間的均方根誤差，平均絕對(duì)誤差、平均絕對(duì)百分誤差、擬合度中的至少一種對(duì)比結(jié)果。

具體的，均方根誤差的計(jì)算公式如下：其中，RMSE為，y_i為赤潮生物量測(cè)量值，y″_i為赤潮生物量預(yù)測(cè)值，N為赤潮生物量測(cè)量值的個(gè)數(shù)，根據(jù)均方根誤差RMSE能夠很好地反映出赤潮生物量預(yù)測(cè)值的精密度。

具體的，平均絕對(duì)誤差的計(jì)算公式為：其中，MAE為平均絕對(duì)誤差，y_i為赤潮生物量測(cè)量值，y″_i為赤潮生物量預(yù)測(cè)值，N為赤潮生物量測(cè)量值的個(gè)數(shù)，根據(jù)平均絕對(duì)誤差MAE能夠很好地反映出赤潮生物量預(yù)測(cè)值的實(shí)際誤差情況。

具體的，平均絕對(duì)百分誤差的計(jì)算公式為：其中，MAPE為平均絕對(duì)百分誤差，y_i為赤潮生物量測(cè)量值，y″_i為赤潮生物量預(yù)測(cè)值，N為赤潮生物量測(cè)量值的個(gè)數(shù)，根據(jù)平均絕對(duì)百分誤差MAPE能夠很好地反映出赤潮生物量預(yù)測(cè)值的準(zhǔn)確性。

具體的，擬合度的計(jì)算公式為：其中，R為平均絕對(duì)百分誤差，y_i為赤潮生物量測(cè)量值，y″_i為赤潮生物量預(yù)測(cè)值，N為赤潮生物量測(cè)量值的個(gè)數(shù)。當(dāng)擬合度的數(shù)值越接近于1，則表明赤潮生物量測(cè)量值與赤潮生物量預(yù)測(cè)值之間的吻合度越強(qiáng)，預(yù)測(cè)效果越佳；反之，則說(shuō)明赤潮生物量測(cè)量值與赤潮生物量預(yù)測(cè)值之間的擬合效果差，預(yù)測(cè)準(zhǔn)確性較差。

基于同一發(fā)明構(gòu)思，本發(fā)明實(shí)施例提供了一種赤潮生物量預(yù)測(cè)裝置，參考圖2所示，包括：獲得單元201，用于獲得M組有效海洋環(huán)境因子，其中，M為所述有效海洋環(huán)境因子的種類(lèi)數(shù)乘以子海域的個(gè)數(shù)，同一組所述有效海洋環(huán)境因子包括多幀在不同時(shí)間點(diǎn)采集的同一種類(lèi)的有效海洋環(huán)境因子；預(yù)處理單元202，用于對(duì)所述M組有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理，以得到對(duì)應(yīng)的M組歸一化后海洋環(huán)境因子；環(huán)境因子預(yù)測(cè)單元203，用于分別將所述M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至對(duì)應(yīng)的ARIMA模型中進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出M組海洋環(huán)境因子預(yù)測(cè)值，m依次為1至M，其中，一個(gè)ARIMA模型對(duì)應(yīng)一組歸一化后海洋環(huán)境因子；赤潮生物量預(yù)測(cè)單元204，用于將所述M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值。

優(yōu)選的，所述獲得單元201包括：

數(shù)據(jù)庫(kù)建立子單元，用于建立原始海洋環(huán)境因子的關(guān)系數(shù)據(jù)庫(kù)，其中，在所述關(guān)系數(shù)據(jù)庫(kù)中保存有多個(gè)歷史監(jiān)測(cè)的原始海洋環(huán)境因子，以及保存有各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)位置信息和監(jiān)測(cè)時(shí)間信息；

提取子單元，用于根據(jù)各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)時(shí)間信息和監(jiān)測(cè)位置信息，從所述關(guān)系數(shù)據(jù)庫(kù)中提取M組原始海洋環(huán)境因子，其中，所述M組原始海洋環(huán)境因子中的同一組原始海洋環(huán)境因子在同一子海域，且距離當(dāng)前時(shí)刻小于預(yù)設(shè)時(shí)間間隔監(jiān)測(cè)到；

篩選子單元，用于從所述M組原始海洋環(huán)境因子中的每組原始海洋環(huán)境因子中篩選出滿足預(yù)設(shè)條件的因子組成為所述M組有效海洋環(huán)境因子。

優(yōu)選的，所述赤潮生物量預(yù)測(cè)裝置還包括：第一建模單元，用于通過(guò)如下步驟建立所述ARIMA模型：

用于根據(jù)所述關(guān)系數(shù)據(jù)庫(kù)中各個(gè)所述原始海洋環(huán)境因子的監(jiān)測(cè)位置信息、監(jiān)測(cè)時(shí)間信息確定針對(duì)所述M組原始海洋環(huán)境因子中每組原始海洋環(huán)境因子的ARIMA模型；

基于AIC準(zhǔn)則分別針對(duì)所述每組原始海洋環(huán)境因子的ARIMA模型進(jìn)行參數(shù)估計(jì)；

選取基于所述AIC準(zhǔn)則進(jìn)行參數(shù)估計(jì)得到最小AIC值時(shí)的參數(shù)組合，作為針對(duì)所述每組原始海洋環(huán)境因子的ARIMA模型的模型參數(shù)。

優(yōu)選的，所述赤潮生物量預(yù)測(cè)裝置還包括：第二建模單元，用于通過(guò)如下步驟建立所述BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型：

獲取赤潮生物量測(cè)量值樣本；

基于如下對(duì)數(shù)公式對(duì)所述赤潮生物量測(cè)量值樣本進(jìn)行歸一化處理得到對(duì)應(yīng)的歸一化后赤潮生物量樣本：

其中，y′為歸一化后赤潮生物量的值，y為赤潮生物量測(cè)量值，y_max為所述赤潮生物量測(cè)量值樣本的最大值，y_min為所述赤潮生物量測(cè)量值樣本的最小值；

基于所述歸一化后赤潮生物量樣本確定BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)以及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)；

基于確定的BP神經(jīng)網(wǎng)絡(luò)的隱含層數(shù)及每一個(gè)隱含層的節(jié)點(diǎn)數(shù)建立針對(duì)所述M組海洋環(huán)境因子預(yù)測(cè)值的BP神經(jīng)網(wǎng)絡(luò)的最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)；

將所述最優(yōu)網(wǎng)絡(luò)結(jié)構(gòu)經(jīng)過(guò)學(xué)習(xí)過(guò)程得到的最優(yōu)值作為所述BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型。

優(yōu)選的，所述赤潮生物量預(yù)測(cè)裝置還包括：對(duì)比單元，用于將所述赤潮生物量預(yù)測(cè)值與對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值進(jìn)行對(duì)比；結(jié)果輸出單元，用于輸出所述赤潮生物量預(yù)測(cè)值與所述對(duì)應(yīng)時(shí)間測(cè)量的赤潮生物量測(cè)量值之間的均方根誤差、平均絕對(duì)誤差、平均絕對(duì)百分誤差、擬合度中的至少一種對(duì)比結(jié)果。

優(yōu)選的，所述有效海洋環(huán)境因子的種類(lèi)包括：酸堿度、溶解氧、水溫、溶解氧飽和度、葉綠素-a、磷酸鹽、氨氮、亞硝酸鹽氮、硝酸鹽氮、鹽度、化學(xué)耗氧量和硅酸鹽。

通過(guò)上述本發(fā)明提供的一個(gè)或多個(gè)實(shí)施例，至少實(shí)現(xiàn)了如下技術(shù)效果或優(yōu)點(diǎn)：

通過(guò)獲得M組有效海洋環(huán)境因子；對(duì)M組有效海洋環(huán)境因子進(jìn)行歸一化預(yù)處理，以得到對(duì)應(yīng)的M組歸一化后海洋環(huán)境因子；分別將M組歸一化后海洋環(huán)境因子中的第m組歸一化后海洋環(huán)境因子輸入至ARIMA模型中進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出M組海洋環(huán)境因子預(yù)測(cè)值；將M組海洋環(huán)境因子預(yù)測(cè)值同時(shí)輸入至BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型進(jìn)行預(yù)測(cè)，以預(yù)測(cè)出在預(yù)設(shè)時(shí)期內(nèi)的赤潮生物量預(yù)測(cè)值。其中，M為有效海洋環(huán)境因子的種類(lèi)數(shù)乘以子海域的個(gè)數(shù)，一組有效海洋環(huán)境因子包括多幀在不同時(shí)間點(diǎn)采集的同一種類(lèi)的有效海洋環(huán)境因子；從而充分考慮了影響赤潮的海洋環(huán)境因子的橫向時(shí)序相關(guān)性及縱向空間異質(zhì)性，能夠針對(duì)不同子海域的不同因子建立不同的ARIMA模型來(lái)描述其時(shí)間上的連續(xù)性及空間異質(zhì)性，從而預(yù)測(cè)出未來(lái)的各海洋環(huán)境因子，再此基礎(chǔ)上利用BP神經(jīng)網(wǎng)絡(luò)來(lái)表達(dá)各海洋環(huán)境因子與赤潮之間的復(fù)雜關(guān)系，進(jìn)而實(shí)現(xiàn)赤潮生物量的短期準(zhǔn)確預(yù)測(cè)，以預(yù)報(bào)未來(lái)幾天是否發(fā)生赤潮，以便進(jìn)行防備。

盡管已描述了本發(fā)明的優(yōu)選實(shí)施例，但本領(lǐng)域內(nèi)的技術(shù)人員一旦得知了基本創(chuàng)造性概念，則可對(duì)這些實(shí)施例作出另外的變更和修改。所以，所附權(quán)利要求意欲解釋為包括優(yōu)選實(shí)施例以及落入本發(fā)明范圍的所有變更和修改。

顯然，本領(lǐng)域的技術(shù)人員可以對(duì)本發(fā)明進(jìn)行各種改動(dòng)和變型而不脫離本發(fā)明的精神和范圍。這樣，倘若本發(fā)明的這些修改和變型屬于本發(fā)明權(quán)利要求及其等同技術(shù)的范圍之內(nèi)，則本發(fā)明也意圖包含這些改動(dòng)和變型在內(nèi)。