張湑行 劉錦源

摘 要：研究表明海洋溫度正在逐年上升，許多海洋生物會(huì)受到溫度影響而進(jìn)行大規(guī)模的遷移，對(duì)于依靠捕撈為生的漁民和捕撈公司這將是嚴(yán)重的打擊。本文主要針對(duì)魚(yú)類(lèi)種群的遷移進(jìn)行研究。首先利用時(shí)間回歸（AR）模型進(jìn)行海水溫度的預(yù)測(cè)。為了研究魚(yú)群位置的長(zhǎng)期變化趨勢(shì)并消除周期性洄游的影響，將計(jì)算的時(shí)間間隔取為一年，結(jié)合魚(yú)類(lèi)選擇適宜環(huán)境的生物學(xué)特點(diǎn)，建立了規(guī)劃模型來(lái)模擬每年魚(yú)群受溫度影響遷移的過(guò)程。最終收集了鯡魚(yú)、鯖魚(yú)的所需參數(shù)，按照以上模型對(duì)兩種魚(yú)群的遷移進(jìn)行了預(yù)測(cè)。

引言

海洋溫度不僅影響海洋生物物種的穩(wěn)定性，而且與人類(lèi)社會(huì)和經(jīng)濟(jì)密切相關(guān)。 隨著溫度的升高，海洋魚(yú)類(lèi)和無(wú)脊椎動(dòng)物可能遷移到高緯度和更深的水域。 在過(guò)去的40年里，海洋溫度的升高影響了全球漁業(yè)的發(fā)展，魚(yú)類(lèi)棲息地的不斷遷移給漁業(yè)公司帶來(lái)了更多的不確定性。為了穩(wěn)定這些公司的生計(jì)，盡量減少全球變暖對(duì)漁業(yè)的影響，研究魚(yú)類(lèi)遷移的規(guī)律有著巨大的價(jià)值。

模型的建立與求解

1.基于AR模型海水溫度預(yù)測(cè)

本文從國(guó)際海洋物理科學(xué)協(xié)會(huì)（ICES）數(shù)據(jù)庫(kù)獲取了1997年-2018年北大西洋海面溫度數(shù)據(jù)。根據(jù)“蘇格蘭海洋地圖集：國(guó)家海洋計(jì)劃”提供的海域參考圖，蘇格蘭海域的范圍為-15 ～ 2 ° 和 55 ～ 62 °。 為了避免魚(yú)類(lèi)繁殖的特殊階段，提取1997年至2018年蘇格蘭（9月，10月，11月）的平均海面溫度進(jìn)行研究。 采用三次樣條法對(duì)每年數(shù)據(jù)進(jìn)行插值，補(bǔ)充選定經(jīng)緯度范圍的溫度數(shù)據(jù)?？梢杂^察到到海水溫度隨經(jīng)度、緯度、時(shí)間變化的大致規(guī)律：當(dāng)經(jīng)度處于-15°到-5°時(shí)，隨經(jīng)度的改變溫度有明顯上升趨勢(shì)，而經(jīng)度處于-5°到0°時(shí)溫度隨經(jīng)度改變有一定的下降趨勢(shì)。各年份高溫區(qū)域（這里指10°C以上）主要分布在緯度55°到58°范圍，且高溫區(qū)域逐漸向高經(jīng)度、低緯度方向移動(dòng)。

結(jié)合以上分析我們將所研究海域按緯度分為三個(gè)部分：D1（55°到58°）、D2（58°到60°）、D3（60°到62°）。將三個(gè)部分海水溫度取平均值，得到三個(gè)部分平均溫度的時(shí)間序列。分別建立自回歸AR模型，并利用數(shù)據(jù)計(jì)算出模型參數(shù)的估計(jì)值，最終模型為（Ti表示對(duì)應(yīng)區(qū)域第i年溫度）：

將三個(gè)部分分別按照模型對(duì)未來(lái)的海水溫度進(jìn)行預(yù)測(cè)，僅展示2031，2041，2051，2061這四年的溫度預(yù)測(cè)結(jié)果：

2.基于規(guī)劃模型魚(yú)群遷移預(yù)測(cè)：

魚(yú)類(lèi)通過(guò)洄游的方式更換生活水域，以滿(mǎn)足各時(shí)期對(duì)外界環(huán)境的需求，進(jìn)而完成繁殖、索餌等生命活動(dòng)并尋找更加適宜的棲息地 。洄游是魚(yú)類(lèi)對(duì)自然環(huán)境變化的長(zhǎng)期適應(yīng)行為 ，能夠使種群獲得更好的繁殖與生存條件。魚(yú)類(lèi)洄游可以分為主動(dòng)洄游和被動(dòng)洄游，其中主動(dòng)洄游是指魚(yú)類(lèi)在進(jìn)行繁殖前種群遷移過(guò)程，而被動(dòng)洄游是指繁殖結(jié)束后回到原來(lái)?xiàng)⒌氐倪w移過(guò)程。鯡魚(yú)、鯖魚(yú)在一年之內(nèi)由于繁殖或季節(jié)變化離開(kāi)棲息地而最終又回到原棲息地，將類(lèi)似的過(guò)程稱(chēng)為周期性洄游;而某年之中外界環(huán)境相比于上一年整體發(fā)生改變而導(dǎo)致魚(yú)群棲息地的遷移稱(chēng)為非周期性遷移，如全球海洋溫度升高而導(dǎo)致的魚(yú)群棲息地的遷移就屬于非周期性遷移。

假設(shè)遷移過(guò)程中魚(yú)群規(guī)模保持不變且分布均勻，可以認(rèn)為魚(yú)群區(qū)域總面積保持不變，進(jìn)一步假設(shè)魚(yú)群的分布為半徑固定的圓形區(qū)域（記為D），由于水溫隨時(shí)間變化導(dǎo)致區(qū)域中心坐標(biāo)不斷變化，將初始時(shí)刻區(qū)域中點(diǎn)坐標(biāo)記為（x0， y0），取相等的時(shí)間間隔Δt，第i個(gè)時(shí)間點(diǎn)所對(duì)應(yīng)的區(qū)域中心坐標(biāo)記為（xi， yi）。i時(shí)刻魚(yú)群所在區(qū)域D可以表示為（其中：r為圓形區(qū)域的半徑）：

為了預(yù)測(cè)魚(yú)群位置長(zhǎng)期變化趨勢(shì)，將時(shí)間間隔Δt取為一年。魚(yú)群的位置受到遷移的影響，由于所取時(shí)間間隔Δt為一年且周期性洄游的周期也為一年，可以認(rèn)為在時(shí)間間隔為一年的前提下魚(yú)群的位置變化主要由非周期性遷移決定，將周期性洄游的影響忽略不計(jì)。根據(jù)分析，當(dāng)水溫發(fā)生變化時(shí)，在一定區(qū)域內(nèi)，魚(yú)群能夠本能的進(jìn)行遷移并不斷地將棲息地調(diào)整到環(huán)境最適宜的位置。按照這一規(guī)律建立規(guī)劃模型模擬魚(yú)群自主尋找棲息地的過(guò)程。

目標(biāo)函數(shù)：

由于溫度受到位置的影響且隨時(shí)間不斷變化，將溫度定義為三元函數(shù)T（x，y，ti），魚(yú)群會(huì)不斷調(diào)整棲息地使溫度盡可能接近其最適值T0，將魚(yú)群所處區(qū)域溫度與最適值的差距定義為目標(biāo)函數(shù)：

按照對(duì)區(qū)域D的定義對(duì)目標(biāo)函數(shù)進(jìn)行積分變換，得到：

約束條件：

需將魚(yú)群活動(dòng)范圍限制在所研究的海洋區(qū)域（記為Dsea）內(nèi)。且在一年的時(shí)間內(nèi)，魚(yú)群需要花費(fèi)大量時(shí)間完成繁殖、索餌、周期性洄游等基本行為，所以考慮非周期性遷移距離有一定限制（最大距離記為tranmax）：

參考了鯖魚(yú)、鯡魚(yú)在蘇格蘭周?chē)Ｓ虻牟稉莆恢眉安稉屏?，提取了其中近年?lái)秋季的捕撈數(shù)據(jù)，根據(jù)這些信息大致確定了兩類(lèi)魚(yú)群中心的初始位置（xi， yi）及魚(yú)群區(qū)域的近似半徑r。通過(guò)查閱魚(yú)類(lèi)遷移的相關(guān)研究，得知了過(guò)去三十年間全球范圍內(nèi)兩種魚(yú)的魚(yú)群最大移動(dòng)距離為400km，所以近似的將兩種魚(yú)群每年最大移動(dòng)距離r選定為13km。確定了模型中的未知參量如下表：

根據(jù)以上參數(shù)并結(jié)合模型二中對(duì)海水溫度的預(yù)測(cè)結(jié)果我們可以預(yù)測(cè)兩種魚(yú)的遷移情況如下：

結(jié)論

結(jié)合實(shí)際情況不能簡(jiǎn)單的把魚(yú)群區(qū)域當(dāng)作固定的圓，實(shí)際分布很可能隨溫度的改變發(fā)生變化。從這個(gè)角度出發(fā)將對(duì)規(guī)劃模型做出調(diào)整。經(jīng)過(guò)調(diào)整后的模型可以提高預(yù)測(cè)魚(yú)群中心位置變化的準(zhǔn)確性，同時(shí)考慮到魚(yú)群外形的變化，這樣更加接近實(shí)際情況。同時(shí)考慮到，即使魚(yú)群中心處于適宜溫度范圍內(nèi)，中心位置以外的魚(yú)也可能進(jìn)行小規(guī)模遷移導(dǎo)致區(qū)域外形改變。在遇到這種情況時(shí)可以用調(diào)整后的模型解決，保持魚(yú)群中心坐標(biāo)進(jìn)行規(guī)劃，即確定了區(qū)域外形的變化。

參考文獻(xiàn)：

[1]常劍波. 水利水電工程對(duì)魚(yú)類(lèi)的影響及減緩對(duì)策[C]. 中國(guó)水利學(xué)會(huì).中國(guó)水利學(xué)會(huì)2008學(xué)術(shù)年會(huì)論文集（上冊(cè)）.中國(guó)水利學(xué)會(huì)：中國(guó)水利學(xué)會(huì)，2008：696-707.

[2]楊慶，胡鵬，楊澤凡，儲(chǔ)立民，楊靜.草魚(yú)洄游的適宜流速條件與適應(yīng)閾值[J].水生態(tài)學(xué)雜志，2019，40（04）：93-100.