馮建綱， 阿布都熱西提·阿布都外力

(新疆大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院, 新疆 烏魯木齊 830046)

干旱區(qū)天然植被作為生物生態(tài)系統(tǒng)中的重要部分，是生態(tài)系統(tǒng)中物質(zhì)循環(huán)和能量流動(dòng)的樞紐，作為生態(tài)環(huán)境治理的一項(xiàng)重要措施，是退化生態(tài)系統(tǒng)恢復(fù)的關(guān)鍵所在[1]，在抑制沙漠化過程和保護(hù)生物多樣性等方面有著重要的生態(tài)意義.水作為干旱區(qū)最重要的生態(tài)環(huán)境因子，不僅是干旱區(qū)綠洲生態(tài)系統(tǒng)構(gòu)成、發(fā)展和穩(wěn)定的基礎(chǔ)和依據(jù)，而且決定著干旱區(qū)綠洲化過程與荒漠化過程這兩類極具對立與沖突性的生態(tài)環(huán)境演化過程[2].水資源是干旱區(qū)綠洲系統(tǒng)形成、存在和發(fā)展的基礎(chǔ)[3]，干旱區(qū)植被的需水來源主要來自于降水，在干旱區(qū)，植物對水的依存性更大，由于降雨量稀少,水資源的有限性，加之日照強(qiáng)烈，缺水使分布于綠洲周邊地區(qū)的耐旱植被退化嚴(yán)重，植被分散而且稀疏，地表覆蓋率很低甚至大面積死亡、消失[4]，造成環(huán)境進(jìn)一步惡化，特別是在干旱、半干旱地區(qū)往往是災(zāi)難性的[5].基于以上的情況，考慮水分脅迫條件下植被的生長模型，在捕食模型的基礎(chǔ)上，根據(jù)數(shù)學(xué)物理原理,給出了如下模型[6-10]：

(1)

式中：a,b,c,n,r,β皆為正常數(shù)；x(t)表示t時(shí)刻區(qū)域內(nèi)地下水資源(土壤含水量)；y(t)表示t時(shí)刻區(qū)域內(nèi)植被密度；r表示地下水的增量(主要是降雨量)；c表示因土壤滲透而導(dǎo)致的水的損失率；n表示單位密度植被蒸散量；β表示因缺水而導(dǎo)致的植被死亡率；bxy表示植被對水的吸收量；axy表示植被的生長量.

1 預(yù)備知識(shí)

為了研究模型的穩(wěn)定性態(tài)，我們需要用到如下知識(shí)點(diǎn).

定理3[13]對應(yīng)于矩陣A的特征根的不同情況，關(guān)于微分方程組的零解的穩(wěn)定性態(tài)有以下結(jié)論：

(1)若特征根均具有負(fù)實(shí)部，則方程組的零解是漸進(jìn)穩(wěn)定的.

(2)若有正實(shí)部的特征根，則方程組的零解是不穩(wěn)定的.

2 主要結(jié)果

定理4 系統(tǒng)(1)在第一象限內(nèi)無閉軌.

證明取Dulac函數(shù)B(x,y)=xkym，其中k,m為待定常數(shù).則

krxk-1ym-(k+1)cxkym-(k+1)bxkym+1+

amxk+1ym-1-(m+1)βxkym

取k=m=-1，則D=-rx-2y-1-ay-2.因?yàn)閞>0，a>0,在第一象限內(nèi)x>0,y>0，故D<0.由Dulac準(zhǔn)則，系統(tǒng)(1)在第一象限內(nèi)無閉軌.

定理5 說明植被的退化是不會(huì)周期發(fā)生的.

下面我們考慮系統(tǒng)(1)的平衡點(diǎn)的穩(wěn)定性情況，考慮如下方程

(2)

(3)

其一次近似系統(tǒng)上的Jacobi矩陣為

(4)

其一次近似系統(tǒng)上的Jacobi矩陣為

λ2+(c+by0)λ+ar-βc=0

定理9 說明當(dāng)降雨量充足時(shí)，植被密度取決于水的供給數(shù)量.

3 結(jié)束語

水是干旱區(qū)各種資源中最為重要的資源之一，是植被生存極其重要的因子.通過對模型的定性分析可知：當(dāng)干旱區(qū)的降雨過少，土壤滲透和蒸騰作用過強(qiáng)時(shí)，植被的生長受到抑制，將導(dǎo)致植被大量死亡，植被稀疏，出現(xiàn)地表裸露.干旱區(qū)植被密度取決于水的數(shù)量，有水成綠洲，無水成荒漠.因此，增加干旱區(qū)水的供給，如對干旱區(qū)采取人工灌溉等輸水工程，植被密度將會(huì)增加，這也是防止植被退化和保護(hù)綠洲規(guī)模的主要手段之一.

[1] 郝博，粟曉玲，馬孝義，等.干旱區(qū)植被生態(tài)需水的研究進(jìn)展[J].水資源與工程學(xué)報(bào)，2009,20(4):1-5.

[2] 陳亞寧，李衛(wèi)紅，徐海量，等.塔里木河下游地下水位對植被的影響[J].地理學(xué)報(bào)，2003,58(4):542-549.

[3] 魏偉，石培基,周俊菊,等.近 20 年和田綠洲水資源變化及其驅(qū)動(dòng)力分析[J].干旱區(qū)資源與環(huán)境，2013,27(2):156-161.

[4] 阿布都熱合曼·哈力克.基于生態(tài)環(huán)境保護(hù)的且末綠洲生態(tài)需水量研究[J].干旱區(qū)資源與環(huán)境，2012,26(7):20-25.

[5] 何炎紅.烏蘭布和沙漠植被與水資源相互影響的研究[D].呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué),2006.

[6] 陳蘭蓀.數(shù)學(xué)生態(tài)學(xué)模型與研究方法[M].北京:科學(xué)出版社,1988.

[7] 王樹禾.微分方程模型與混沌[M].合肥:中國科技大學(xué)出版社,1999.

[8] 李婷,劉萍,史峻平,等.一類基于一階傳遞偏微分方程的水污染模型[J].哈爾濱師范大學(xué)自然科學(xué)學(xué)報(bào)，2009，25(1)：4-6.

[9]陳軍杰,潘國衛(wèi). 一個(gè)具暫時(shí)免疫且總?cè)丝跀?shù)可變的傳染病模型[J].生物數(shù)學(xué)報(bào),2003,18(4):401-405.

[10]Li M Y,Graef J R, Wang L C,etal. Global dynamic of SEIR model with varying total population size [J]. Mathematical Biosciences, 1999, 160:191-195.

[11] 馬知恩,周義倉.常微分方程定性與穩(wěn)定性方法[M].北京:科學(xué)出版社,2001.

[12] 葉宗澤,楊萬祿.常微分方程組與運(yùn)動(dòng)穩(wěn)定性理論[M].天津:天津大學(xué)出版社,1985.

[13] 蔡燧林,盛驟.常微分方程與穩(wěn)定性理論[M].北京:高等教育出版社,1986.