一類帶白噪聲和Lévy跳躍的營養(yǎng)-浮游植物模型的分析

譚 楊， 郭子君， 楊 林， 王海揚

(1. 銅仁職業(yè)技術(shù)學(xué)院， 貴州 銅仁 554300；2. 華南農(nóng)業(yè)大學(xué) 理學(xué)院應(yīng)用數(shù)學(xué)研究所， 廣東 廣州 510275)

0 引 言

浮游生物生長在湖泊、 河流及海洋中， 其中的浮游植物是指在水中以浮游方式生活的微小植物， 通常浮游植物就是指浮游藻類. 浮游植物直接或間接地是魚類等水生動物的食物來源， 在人類的水產(chǎn)經(jīng)濟活動中起著重要的作用[1]. 隨著世界經(jīng)濟的發(fā)展， 人類活動產(chǎn)生的大量廢棄物流入水體， 導(dǎo)致水體的富營養(yǎng)化， 這樣的結(jié)果導(dǎo)致了藻華現(xiàn)象的出現(xiàn)， 即浮游植物在短時間內(nèi)爆發(fā)式的增長. 浮游植物自身巨大的耗氧量和有機物以及天氣的變化都可能導(dǎo)致大量藻類死亡， 從而導(dǎo)致水質(zhì)迅速惡化到變黑甚至發(fā)臭， 影響了魚類資源的自然循環(huán)， 也會造成巨大的經(jīng)濟損失[2-3].

近些年來， 利用數(shù)學(xué)模型的方法研究浮游生物群落的發(fā)展規(guī)律成為生物數(shù)學(xué)的一個重要課題， 已經(jīng)建立了一些浮游生物系統(tǒng)的數(shù)學(xué)模型[4-6]. 例如， 文獻[4] 研究了幾類營養(yǎng)-浮游植物模型， 分析了產(chǎn)毒浮游植物(Toxin-Producing Phytoplankton, TPP)釋放的有毒化學(xué)物質(zhì)季節(jié)性反復(fù)出現(xiàn)的水華現(xiàn)象的動力學(xué)性質(zhì)； 文獻[5] 研究了一類具有非線性互抑制項和反饋控制的非自治化感作用浮游植物模型， 得到了系統(tǒng)持續(xù)存在及滅絕的一些充分條件； 在文獻[6]中， 作者提出了一類具有年齡結(jié)構(gòu)和非局部效應(yīng)的水柱中浮游植物成熟、 生長和空間分布的反應(yīng)擴散模型， 其結(jié)果表明浮游植物的死亡率和成熟時間都會影響模型的動力學(xué)特性. 本文研究的模型是基于文獻[4]中的一個確定性模型(1).

(1)

該模型假設(shè)該項浮游植物本身不產(chǎn)生毒素， 而環(huán)境中的毒素(其他浮游植物產(chǎn)生)對系統(tǒng)中的浮游植物有影響. 其中，x(t)和y(t)分別為營養(yǎng)物質(zhì)和浮游植物在t時刻的數(shù)量或密度.其他常數(shù)的意義如下：a為外部營養(yǎng)物質(zhì)流入的速率；b為營養(yǎng)物質(zhì)的吸收速率；β為接觸率；c(c

在現(xiàn)實環(huán)境中， 種群的發(fā)展不可避免地會受到隨機因素的干擾， 環(huán)境隨機性影響著種群發(fā)展的波動(如持續(xù)存在或滅絕). 因此， 許多研究者在確定性模型中引入隨機白噪聲來描述隨機噪聲的影響[7-10]. 文獻[7]研究了隨機干擾下的具有間接影響的浮游生物群落捕食模型的長期動態(tài)， 證明了在給定任何非負的正初始條件下， 模型的解是非負有界的， 并利用經(jīng)典隨機動力系統(tǒng)理論， 證明了文中的隨機系統(tǒng)具有唯一的隨機吸引子； 文獻[8]利用隨機模型研究了產(chǎn)毒浮游植物對種群持續(xù)性的影響， 證明了在隨機模型中， 無論是輸入營養(yǎng)物濃度較大或較小， 產(chǎn)毒浮游植物都可以終止有害藻華的發(fā)生； 在文獻[9]中， 作者研究了一類基于模型(1)的隨機營養(yǎng)-浮游植物模型， 得到了依平均的滅絕和持續(xù)存在的充分條件， 并研究了遍歷平穩(wěn)分布及周期解的存在性； 文獻[10] 也基于確定性模型(1)提出了具有白噪聲及電報噪聲的隨機模型， 得到了滅絕性及遍歷平穩(wěn)分布的充分條件. 本文的研究不同于文獻[9-10]， 考慮了模型(1)的接觸率β受到的白噪聲干擾：β→β+σB(t)， 其中，σ為白噪聲的干擾強度，B(t)為標(biāo)準(zhǔn)布朗運動， 模型為

(2)

眾所周知， 模型(2)的解是連續(xù)的， 但該模型無法解釋遇到突然或其他嚴(yán)重擾動時的現(xiàn)象， 如熱沖擊、 壓力沖擊和催化劑等因素的突然加入[11]. 為了描述這類環(huán)境噪聲， 目前使用的比較多的方法就是用帶跳躍的隨機微分方程來研究種群模型的動力學(xué)行為[12]. 一般來說， 跳躍時間總是隨機的， 跳躍的等待時間與Lévy跳躍相似， Lévy過程本質(zhì)上是具有平穩(wěn)和獨立增量的隨機過程. 目前， 這類的研究也較多， 例如： 文獻[13]研究了一類具有Lévy跳躍的Beddington-DeAngelis型功能反應(yīng)的捕食-被捕食模型， 得到了滅絕、 依平均不持久及弱持續(xù)存在的充分條件， 并證明了Lévy跳躍的變化能夠影響系統(tǒng)的漸近性質(zhì)； 在文獻[14]中， 作者分析了具有Allee效應(yīng)的隨機種群模型的漸近行為， 并給出了該擾動模型漸近行為的判定條件； 文獻[15]研究了一類具有Lévy跳躍和白噪聲干擾的飽和發(fā)病率和雙重傳染病的隨機SIS流行病模型， 得到了兩種傳染病平均滅絕和持續(xù)的充分條件； 文獻[16]研究了一類基于白噪聲及Lévy跳躍的隨機SIRS模型， 得到了基本滅絕及持久的充分條件， 并證明了大的噪聲強度能夠抑制疾病爆發(fā)的結(jié)果. 本文將Lévy跳躍引入隨機模型(2)， 假設(shè)接觸率系數(shù)β又受到Lévy跳躍的擾動， 得到模型

(3)

下面給出本文的一個必要假設(shè)：

本文證明了對于任意給定的整初始值， 模型(3) 存在唯一的全局正解； 給出了兩個浮游植物滅絕的充分條件； 研究了浮游植物依平均持續(xù)存在的充分條件； 給出了所得結(jié)果的數(shù)值模擬.

1 模型解的存在及其唯一性

為了研究模型(3)的動力學(xué)性質(zhì)， 首先要考慮模型的解是否全局存在， 是否為正. 下面證明模型(3)的解是全局正解.

證明模型(3)的系數(shù)滿足局部利普希茨條件， 則在區(qū)間[0,τe)上存在唯一的局部解(x(t),y(t))， 其中τe為爆炸時刻[14-15].若能證明τ∞=∞ a.s.， 則該局部解是全局存在的. 定義停時

P{τk≤T}≥ε，k≥k1.

(4)

對任意k， 當(dāng)t≥τk時，

d(cx+by)=

[ca-ξ(cx+by)]dt，

其中,ξ=min{u，ε-m}， 因為b>c，ε>m， 則

dV(x,y)=LV(x,y)dt-bσydB(t)+cσxdB(t)-

(5)

其中，

LV(x,y)=

又cx(t)+by(t)≤K， 于是

(6)

由假設(shè)A可知1-bγ(v)y>0， 再由泰勒公式可得

-cln(1-bγ(v)y)-bcγ(v)y=

同理可得

-bln(1+cγ(v)x)+bcγ(v)x=

因此

(7)

對式(5)兩邊從0到ζk∧T積分可得

(8)

對式(8)兩邊取期望可得

EV(x(τk∧T),y(τk∧T))≤EV(x(0),y(0))+

K0E(τk∧T)≤EV(x(0)，y(0))+K0T.

(9)

對任意k≥k1， 令Ωk={τk≤T}， 則由式(4)可知P(Ωk)≥ε， 于是對所有的ω∈Ωk， 使得x(τk，ω)或y(τk，ω)等于k或1/k， 再由ε>0的任意性有

V(x(0)，y(0))+K0T≥

E[IΩk(ω)V(x(τk∧T)，y(τk∧T)]≥

其中，IΩk是Ωk的示性函數(shù).令k→∞， 則有

∞>V(x(0)，y(0))+K0T=∞，

所以， 必然有τ∞=∞ a.s.， 證畢.

2 滅絕性

在研究生物數(shù)學(xué)模型的動力學(xué)性質(zhì)時， 主要關(guān)心的是如何調(diào)節(jié)有害種群的發(fā)展動態(tài)， 使有害種群長期消失. 本節(jié)研究浮游植物滅絕的充分條件.

證明由It公式有

dlny(t)=

(10)

對式(10)兩邊從0到t積分并除以t可得

(11)

由泰勒公式有

因此,

(12)

其中

(13)

因此， 對式(13)兩邊求極限可得

(14)

則

于是

證畢.

注情形1說明只要隨機干擾的強度足夠大， 必定造成浮游植物種群將依指數(shù)趨于滅絕.

3 持久性

在生物系統(tǒng)模型的研究中， 持久性是最重要的性質(zhì)之一. 本節(jié)研究浮游植物依平均持續(xù)存在的充分條件. 為此， 首先引用如下定義：

定義1[20]系統(tǒng)(3)被稱為依平均持續(xù)存在， 若

定理3若假設(shè)A和如下條件

成立， 則系統(tǒng)(3)將依平均持續(xù)存在.

證明由式(11)可知

則由泰勒公式有

-ln(1+cγ(v)x)+cγ(v)x=

因此,

(15)

其中，

又因

則

(16)

式(16)結(jié)合式(15)可得

(17)

則有

(18)

式(18)兩邊取極限可得

(19)

證畢.

4 數(shù)值模擬

采用類似文獻[19]的方法模擬上述得到的結(jié)果.

1) 令a=3，β=0.3，b=1，c=0.9，ε=2，u=1.2，m=0.1，μ=0.6，θ=0.01，σ=0.3，γ(v)=0.15，λ(Z)=1，x(0)=0.5，y(0)=1.8， 則

即定理2的情形1滿足， 則浮游植物種群趨于滅絕 (見圖 1).

圖 1 浮游植物種群趨于滅絕的情形1

2) 令a=3，β=0.5，b=1，c=0.9，ε=2，u=1.2，m=0.1，μ=0.6，θ=0.01，σ=0.15，γ(v)=0.15，λ(Z)=1，x(0)=0.5，y(0)=1.8， 則

即定理2的情形2滿足， 則浮游植物種群趨于滅絕 (見圖 2).

3) 令a=2，β=0.5，b=1，c=0.9，ε=0.4，u=0.6，m=0.1，μ=6，θ=0.1，σ=0.1，γ(v)=0.05，λ(Z)=1，x(0)=0.5，y(0)=1.8. 則

圖 2 浮游植物種群趨于滅絕的情形2

即定理3的條件滿足， 則浮游植物種群依平均持續(xù)存在(見圖 3).

圖 3 浮游植物種群依平均持續(xù)存在

4) 令a=1，b=1，β=0.8，c=0.6，ε=0.4，u=0.1，m=0.2，μ=0.6，θ=0.02，σ=0.6，γ(v)=0.5，λ(Z)=1，x(0)=0.5，y(0)=1.8. 如圖 4 所示， 確定模型中浮游植物持續(xù)存在而隨機模型中是滅絕的.

圖 4 浮游植物種群在確定性模型中持續(xù)存在而在隨機模型中趨于滅絕的情形

5) 令前一種模擬4)中的白噪聲強度為σ=0.3 而其他值不變， 如圖 5 所示， 確定模型和隨機模型中的浮游植物種群均持續(xù)存在.

圖 5 確定模型和隨機模型中的浮游植物種群均持續(xù)存在

5 結(jié) 論

本文研究了隨機因素干擾下具有Lévy跳躍的營養(yǎng)物質(zhì)-浮游植物模型， 使用了Lévy跳躍過程來模擬突然的環(huán)境擾動對營養(yǎng)物質(zhì)與浮游植物的接觸率的影響. 定理1保證了模型理論解的存在唯一性； 定理2表明， 在一些簡單的假設(shè)下， 大強度的白噪聲及Lévy噪聲有可能造成浮游植物的滅絕. 另外， 定理3得出了營養(yǎng)物質(zhì)與浮游植物種群依平均持續(xù)存在的充分條件.