基于Elman神經(jīng)網(wǎng)絡(luò)的阿拉善荒漠嚙齒動物群落組成物種數(shù)量預(yù)測研究

盧志宏，武曉東*，郭利彪，付和平，滿都呼，岳闖，柴享賢，包達(dá)爾罕，楊素文，葉麗娜，李燕妮

1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)生態(tài)環(huán)境學(xué)院，內(nèi)蒙古 呼和浩特 010018；2. 草業(yè)與草地資源教育部重點實驗室，內(nèi)蒙古 呼和浩特 010018；3. 中國農(nóng)業(yè)科學(xué)研究院草原科學(xué)研究所，內(nèi)蒙古 呼和浩特 010010

基于Elman神經(jīng)網(wǎng)絡(luò)的阿拉善荒漠嚙齒動物群落組成物種數(shù)量預(yù)測研究

盧志宏1,2，武曉東1,2*，郭利彪3，付和平1,2，滿都呼1,2，岳闖1,2，柴享賢1,2，包達(dá)爾罕1,2，楊素文1,2，葉麗娜1,2，李燕妮1,2

1. 內(nèi)蒙古農(nóng)業(yè)大學(xué)生態(tài)環(huán)境學(xué)院，內(nèi)蒙古 呼和浩特 010018；2. 草業(yè)與草地資源教育部重點實驗室，內(nèi)蒙古 呼和浩特 010018；3. 中國農(nóng)業(yè)科學(xué)研究院草原科學(xué)研究所，內(nèi)蒙古 呼和浩特 010010

群落的數(shù)量變動及預(yù)測是生態(tài)學(xué)研究的重要內(nèi)容，將神經(jīng)網(wǎng)絡(luò)技術(shù)應(yīng)用到嚙齒動物群落數(shù)量預(yù)測中是一種新嘗試。Elman神經(jīng)網(wǎng)絡(luò)通過在前饋網(wǎng)絡(luò)中增加延時算子，實現(xiàn)了記憶功能，能夠?qū)X動物組成物種數(shù)量進行動態(tài)模擬和預(yù)測。以騰格里沙漠東緣荒漠為試驗區(qū)，以嚙齒動物數(shù)量為研究對象，采用標(biāo)志重捕法，逐月監(jiān)測2006─2014年每年的4─10月捕獲量，建立Elman神經(jīng)網(wǎng)絡(luò)預(yù)測模型，利用2006─2013年的捕獲量建立訓(xùn)練網(wǎng)絡(luò)，以2014年的數(shù)據(jù)進行驗證與測試，比較3種數(shù)據(jù)處理方法建立預(yù)測模型后的平均誤差和擬合度，確立最優(yōu)模型，預(yù)測阿拉善荒漠嚙齒動物組成物種數(shù)量動態(tài)。結(jié)果表明：（1）未經(jīng)歸一化處理預(yù)測結(jié)果的平均誤差mse為5.30，最小誤差1.52%，擬合度為0.80；（2）經(jīng)[0, 1]歸一化處理的預(yù)測結(jié)果平均誤差mse為4.51，最小誤差1.54%，擬合度為0.82；（3）經(jīng)[-1, 1]歸一化處理預(yù)測結(jié)果的平均誤差mse為5.03，最小誤差1.63%，擬合度為0.69；（4）3種歸一化處理后Elman神經(jīng)網(wǎng)絡(luò)模型差異不顯著。通過平均誤差和擬合度的比較，文章認(rèn)為采用[0, 1]歸一化建立的Elman神經(jīng)網(wǎng)絡(luò)能較好的預(yù)測荒漠嚙齒動物數(shù)量的變化規(guī)律，應(yīng)用該網(wǎng)絡(luò)可以預(yù)測阿拉善荒漠嚙齒動物組成物種數(shù)量變化趨勢，對指導(dǎo)當(dāng)?shù)厥笄楸O(jiān)控和防治具有重要意義。

Elman神經(jīng)網(wǎng)絡(luò)；阿拉善荒漠；嚙齒動物；標(biāo)志重捕法

LU Zhihong, WU Xiaodong, GUO Libiao, FU Heping, MAN Duhu, YUE Chuang, CHAI Xiangxian, BAO Daerhan, YANG Suwen, YE Lina, LI Yanni. Prediction of the Number of Rodent Community Composition Species Based on Elman Neural Network in Alasan Desert [J]. Ecology and Environmental Sciences, 2015, 24(12): 1976-1982.

種群的數(shù)量動態(tài)主要是研究種群組成數(shù)量在時間和空間上的變動規(guī)律及其調(diào)節(jié)機制，是種群生態(tài)學(xué)研究的核心理論和核心問題，為此生態(tài)學(xué)家提出了許多數(shù)學(xué)模型進行模擬和預(yù)測種群數(shù)量動態(tài)（孫儒泳，2001）。長期以來，種群生態(tài)學(xué)家通過大量研究揭示動物種群的數(shù)量變動，認(rèn)為調(diào)節(jié)機制復(fù)雜多變，提出了許多種群數(shù)量動態(tài)研究的假說和理論（Chitty，1960；Christian et al.，1964；Wynne-Edwards，1962；Krebs et al.，1973；張志彬等，1998；張志彬，1999；Lima et al.，1999）。

前人應(yīng)用的主要預(yù)測方法是建立線性回歸方程，運用不同的預(yù)測指標(biāo)，針對特定的嚙齒動物建立一元或多元回歸模型預(yù)測種群數(shù)量或發(fā)生等級（陳安國等，1981；朱盛侃等，1993；姜運良等，1994；王勇等，1997；董維惠等，1998；王利民等，1998；侯希賢等，2000a，2000b；楊再學(xué)等，2010）。由于建立回歸方程的指標(biāo)、研究對象不同，建立的預(yù)測模型存在較大差異，同時人們越來越發(fā)現(xiàn)用線性方程的關(guān)系反應(yīng)種群數(shù)量的變動預(yù)測，具有較大的局限性。自然界中物種種群的數(shù)量變動與環(huán)境之間的關(guān)系并非是簡單的線性關(guān)系。因此，對于嚙齒動物種群數(shù)量動態(tài)預(yù)測，研究人員運用了指數(shù)法、灰色理論、馬爾科夫鏈進行了模擬預(yù)測（何淼等，1996；王玉志等，1997；馮志勇等，2000；武文華等，2007）；戴應(yīng)貴（2001）研究了四川邛崍縣農(nóng)田生境中嚙齒動物結(jié)合群落密度動態(tài)變化規(guī)律。

目前，研究嚙齒動物種群數(shù)量預(yù)測絕大多數(shù)只限于在種群水平上的研究，很少在群落水平上研究物種的數(shù)量動態(tài)及其預(yù)測，而未來種群和群落生態(tài)學(xué)發(fā)展的趨勢是更加關(guān)注種間關(guān)系、群落動態(tài)及其維持機制，同時研究種群維持機制一直是種群生態(tài)學(xué)討論的核心問題，也是生態(tài)學(xué)家長期面臨的挑戰(zhàn)。因此，應(yīng)用神經(jīng)網(wǎng)絡(luò)技術(shù)研究嚙齒動物群落組成物種的數(shù)量動態(tài)機制具有重要的理論意義和實際意義。

人工神經(jīng)網(wǎng)絡(luò)（Artificial Neural Network，ANN），簡稱神經(jīng)網(wǎng)絡(luò)，是由大量簡單的信息單元（即神經(jīng)元）通過廣泛地相互連接組成的復(fù)雜網(wǎng)絡(luò)，其建立是以仿生物學(xué)對人腦組織結(jié)構(gòu)進行了全面系統(tǒng)研究為基礎(chǔ)，可以模擬大腦組成結(jié)構(gòu)及思維活動。人工神經(jīng)網(wǎng)絡(luò)是分析輸入變量和輸出變量兩者之間復(fù)雜的非線性全局變化映射關(guān)系的數(shù)據(jù)處理方法。從理解模式和分析模式兩者關(guān)系的本質(zhì)上來看，由輸入神經(jīng)元和輸出神經(jīng)元構(gòu)成的神經(jīng)網(wǎng)絡(luò)屬于非線性動力系統(tǒng)，通過對已知事件進行學(xué)習(xí)理解分析及適應(yīng)，實現(xiàn)非線性分析，相比較于傳統(tǒng)的數(shù)據(jù)分析方法，神經(jīng)網(wǎng)絡(luò)具有以下優(yōu)勢：（1）處理邏輯關(guān)系模糊的數(shù)據(jù)；（2）特征表現(xiàn)模糊的數(shù)據(jù)；（3）非線性特征的數(shù)據(jù)；（4）存在噪聲的隨機數(shù)據(jù)。最突出的特點是任何能用傳統(tǒng)的統(tǒng)計方法解決的問題，用神經(jīng)網(wǎng)絡(luò)能處理得更好、更準(zhǔn)確、更科學(xué)（祖元剛，2004）。目前神經(jīng)網(wǎng)絡(luò)已被廣泛應(yīng)用于工業(yè)控制、農(nóng)業(yè)產(chǎn)量預(yù)測、水華預(yù)測、氣象預(yù)測等諸多領(lǐng)域的擬合和預(yù)測（徐興梅，2014；殷高方等，2012；彭昱忠，2015）。

生態(tài)學(xué)數(shù)據(jù)多具有隨機性特征，由于生物的進化和環(huán)境的變化，采用確定的數(shù)學(xué)方程預(yù)測和分析生態(tài)系統(tǒng)行為是十分困難的，人工神經(jīng)網(wǎng)絡(luò)在繼承了統(tǒng)計數(shù)據(jù)處理的優(yōu)點之外，表現(xiàn)出了數(shù)據(jù)處理的靈活性和分析結(jié)構(gòu)較強的規(guī)律性，可以將不同種類、不同性質(zhì)的變量聯(lián)系在一起，數(shù)字化系統(tǒng)信息，并將客觀存在系統(tǒng)的特征映射到數(shù)字化系統(tǒng)中，由此建立的模型能夠模擬和分析生態(tài)系統(tǒng)的本質(zhì)狀態(tài)。Elman神經(jīng)網(wǎng)絡(luò)是在BP神經(jīng)網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)增加延時算子達(dá)到輸出神經(jīng)元與下一時刻輸入隱含層神經(jīng)元的同步性，從而實現(xiàn)了記憶過去動態(tài)的功能，適合于時間序列的預(yù)測研究（陳明，2013）。相關(guān)人員對存在時間序列變化規(guī)律的事件進行了預(yù)測研究，訓(xùn)練曲線光滑，收斂速度快，預(yù)測精度高（吳萍等，2007；王中凱等，2013；崔東文，2013；邵月紅等，2014）。

對環(huán)境造成壓力甚至產(chǎn)生危害的往往是多種嚙齒動物共同作用產(chǎn)生的結(jié)果，嚙齒動物包括不同科屬的物種，生活型亦不同，對環(huán)境的利用及產(chǎn)生的影響亦不同，因此應(yīng)用Elman神經(jīng)網(wǎng)絡(luò)模擬預(yù)測嚙齒動物群落組成物種數(shù)量變動規(guī)律，可以更好地研究群落變動機理。本研究嘗試應(yīng)用Elman神經(jīng)網(wǎng)絡(luò)模型對阿拉善荒漠嚙齒動物群落組成物種數(shù)量進行模擬預(yù)測研究，為評估鼠害危害程度、鼠疫流行強度及制定防治方案提供科學(xué)的理論指導(dǎo)。

1　研究地概況與方法

1.1研究地概況

研究區(qū)（E104°10′～E105°30′，N37°24′～N38°25′）位于內(nèi)蒙古阿拉善左旗南部荒漠，是典型的溫性荒漠草原，植物種類少，密度低，蓋度一般在1%～20%之間，多為旱生、超旱生及鹽生植物。冬季嚴(yán)寒，夏季酷熱，全年氣候干燥，晝夜溫差大。2003─2014年最低氣溫-36 ℃，最高氣溫42 ℃；降水量45～288 mm，主要集中在7─8月份。年平均氣溫9.0 ℃，無霜期156 d；年潛在蒸發(fā)量3000～4700 mm。土壤貧瘠，質(zhì)地松散，淋溶作用微弱。主要嚙齒動物有子午沙鼠（Meriones meridianus）、三趾跳鼠（Dipus sagitta）、五趾跳鼠（Allacraga sibirica）、長爪沙鼠（Meriones unguiculatus）、小毛足鼠（Phodopus roborovskii）、 阿 拉 善 黃 鼠 （ Spermophilus alaschanicus）等。

1.2研究方法

在阿拉善左旗孿井灘建立固定的標(biāo)志重捕樣地，樣地面積為2 hm2，以方格式布籠，籠距15 m，每個網(wǎng)格點上布設(shè) 1個活捕籠（42 cm×17 cm×13 cm），共計112個活捕籠，以新鮮花生米為誘餌。于2006─2014年每年的4─10月，每月月初連捕4 d。記錄捕獲個體種名、性別、繁殖狀況、體重及捕獲位置，統(tǒng)計捕獲量。以連續(xù)4 d的捕獲個體數(shù)量（不計重捕個體）表示種群數(shù)量。2006─2014年4─10月捕獲比例見表1。

表1　2006─2014年4─10月捕獲比例Table 1 The catch proportion of April to October in 2006─2014 %

本文采用[0, 1]、[-1, 1]歸一化方法對原始數(shù)據(jù)進行處理，對比研究不同歸一化神經(jīng)網(wǎng)絡(luò)模型對荒漠嚙齒動物組成物種數(shù)量變化規(guī)律。歸一化是指將不同量綱和數(shù)量級的數(shù)據(jù)轉(zhuǎn)變成具有可比性的相同量綱和數(shù)量級的數(shù)據(jù)，反歸一化是歸一化的逆過程。歸一化和反歸一化公式如下。

其中公式（1）、（2）表示將原始數(shù)據(jù)歸一化到[0, 1]區(qū)間及反歸一化，xi表示第i個數(shù)據(jù)，xmax表示所有輸入數(shù)據(jù)的最大值，xmin表示所有輸入數(shù)據(jù)的最小值， y1i表示經(jīng)[0, 1]區(qū)間歸一化所得到的相應(yīng)的值， x1i′表示 Elman神經(jīng)網(wǎng)絡(luò)模擬得到的預(yù)測值， y1i′表示經(jīng)[0, 1]區(qū)間反歸一化后對應(yīng)的預(yù)測值；公式（3）、（4）表示將原始數(shù)據(jù)歸一化到[-1, 1]區(qū)間及反歸一化，同xi、xmax、xmin同上，其中常數(shù)項（-1）的作用是將歸一化后的值限定在[-1, 1]之間，

y2i表示表示經(jīng)[-1, 1]區(qū)間歸一化所得到的相應(yīng)的值， x2i′表示Elman神經(jīng)網(wǎng)絡(luò)模擬得到的預(yù)測值， y2i′表示經(jīng)[-1, 1]區(qū)間反歸一化后對應(yīng)的預(yù)測值。

1.3Elman神經(jīng)網(wǎng)絡(luò)原理和方法

1.3.1Elman神經(jīng)網(wǎng)絡(luò)網(wǎng)絡(luò)原理

Elman神經(jīng)網(wǎng)絡(luò)是典型的反饋神經(jīng)網(wǎng)絡(luò)，通過在隱含層中增加承接作用的延時算子，使上一時刻輸出神經(jīng)元與下一時刻輸入神經(jīng)元同時作為隱含層輸入單元，實現(xiàn)了記憶功能，使網(wǎng)絡(luò)具有適應(yīng)時變能力，能夠反應(yīng)系統(tǒng)動態(tài)過程的特征（高瑋等，2001；王小川等，2013）。Elman神經(jīng)網(wǎng)絡(luò)由Elman于 1990年針對語音處理問題提出來的，該神經(jīng)網(wǎng)絡(luò)通過延時算子將隱含層輸出神經(jīng)元延遲到與下一時刻輸入神經(jīng)元信息同時作為輸入信號，實現(xiàn)了對歷史數(shù)據(jù)進行反饋，增加了Elman網(wǎng)絡(luò)處理輸入與輸出數(shù)據(jù)的能力，從而達(dá)到了動態(tài)建模的目的（劉天舒，2001）。

反饋型神經(jīng)網(wǎng)絡(luò)，又稱遞歸網(wǎng)絡(luò)，注重全局穩(wěn)定性方面的研究，是反饋動力學(xué)系統(tǒng)，反向傳播的實現(xiàn)比前向傳播更具強大的計算能力，主要作用是優(yōu)化和預(yù)測，Elman網(wǎng)絡(luò)是一種重要的預(yù)測模型。

Elman型回歸神經(jīng)元網(wǎng)絡(luò)一般分為4層：輸入層、中間層（隱含層）、承接層和輸出層（圖 1）。輸入層主要是傳輸信號，輸出層通過加權(quán)計算實現(xiàn)信息輸出，隱含層實現(xiàn)前一時刻和下一時刻信號的傳遞，有線性和非線性函數(shù)；承接層實現(xiàn)了上一時刻的輸出信號與下一時刻輸入信號同時輸入，可以認(rèn)為這一步是一個延時算子，可實現(xiàn)記憶的功能。以圖1（崔東文，2013；王小川，2013）為例，Elman網(wǎng)絡(luò)的非線性狀態(tài)空間表達(dá)式（崔東文，2013；王小川，2013）為：

圖1　Elman網(wǎng)絡(luò)的非線性狀態(tài)空間表達(dá)式Fig. 1 Nonlinear state space expression of the Elman NeuralNetwork

式中：y、x、u、xc分別表示m維輸出節(jié)點向量、n維中間層節(jié)點單元向量、r維輸入向量和n維反饋狀態(tài)向量；w3、w2、w1分別表示中間層到輸出層、輸入層到中間層、承接層到中間層的連接權(quán)值；g為輸出神經(jīng)元的傳遞函數(shù)，是中間層輸出的線性組合；f為中間層神經(jīng)元的傳遞函數(shù)，常采用S函數(shù)。

Elman神經(jīng)網(wǎng)絡(luò)也采用BP算法進行權(quán)值修正，學(xué)習(xí)指標(biāo)函數(shù)采用誤差平方和函數(shù)。表達(dá)式為：

式中：y^k(k)為目標(biāo)輸入向量。

Elman神經(jīng)網(wǎng)絡(luò)通過前饋神經(jīng)網(wǎng)絡(luò)中加入延時算子，從而加強了記憶功能，通過誤差反向傳播，逐層調(diào)整連接權(quán)值和閾值，最終確定最優(yōu)的權(quán)值和閾值，并將其儲存，進而對真實數(shù)據(jù)進行模擬和預(yù)測，從而解決實際問題。

1.3.2Elman神經(jīng)網(wǎng)絡(luò)的設(shè)計

在MATLAB神經(jīng)網(wǎng)絡(luò)工具箱中函數(shù)initnw初始化函數(shù)及其參數(shù)值確定網(wǎng)絡(luò)權(quán)值及初始值，以Nguyen-Widrore網(wǎng)絡(luò)初始化方法進行，權(quán)值和閾值的方向傳播函數(shù)為 leargdm。Elman網(wǎng)絡(luò)訓(xùn)練采用train函數(shù)，在梯度下降BP算法的基礎(chǔ)上，與目標(biāo)序列進行比較并產(chǎn)生一序列誤差進行反向傳播，以確定每一個權(quán)值和閾值的誤差梯度，進行反向權(quán)值修正，從而提高網(wǎng)絡(luò)性能。訓(xùn)練結(jié)束后，采用 sim函數(shù)來仿真網(wǎng)絡(luò)輸出，用plot函數(shù)繪制誤差曲線，用 postreg函數(shù)對訓(xùn)練結(jié)果進行分析，用誤差性能函數(shù)mse檢驗誤差大小。

1.3.3樣本設(shè)計

本研究以嚙齒動物群落物種數(shù)量作為研究對象，以歷史觀測數(shù)據(jù)（表 1）作為研究樣本，是按時間序列排列的數(shù)字序列，它們之間具有某種統(tǒng)計意義上的關(guān)系，該樣本是隨時間變化而變化的。高瑋等（2001）指出在時序分析方面采用遞歸網(wǎng)絡(luò)更合理。因此引入對時間序列和歷史數(shù)據(jù)更為敏感的遞歸型Elman神經(jīng)網(wǎng)絡(luò)動態(tài)預(yù)測模型，利用連續(xù)N（N≥1）個的數(shù)據(jù)來預(yù)測未來 M個時刻的值（王小川等，2013；劉天舒，2011），從而對嚙齒動群落物種群數(shù)量進行分析。

嚙齒動物種群數(shù)量是基于時間序列對復(fù)雜環(huán)境因子的綜合表現(xiàn)，種群數(shù)量動態(tài)變化受到歷史種群數(shù)量的制約，因此選擇2006─2013年這8年的數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的訓(xùn)練樣本，前3年的種群數(shù)量作為輸入向量，第4年的種群數(shù)量作為目標(biāo)向量，可以得到6組訓(xùn)練數(shù)據(jù)，2014年的數(shù)據(jù)作為網(wǎng)絡(luò)的測試樣本，用于驗證預(yù)測的準(zhǔn)確性。

2　結(jié)果與分析

2.1未歸一化Elman神經(jīng)網(wǎng)絡(luò)

由表2可知，4─9月份的預(yù)測誤差均小于25%， 8月份的預(yù)測誤差僅為1.52，誤差百分比為5.24%，10月預(yù)測值和實測值的差值較大，預(yù)測趨勢和實際變化趨勢基本吻合。由圖2可知，預(yù)測值變化趨勢接近實測值，經(jīng)回歸分析可知，擬合度為 0.80（P=0.0065）。2014年平均誤差mse為5.30。

表2　2014年未歸一化Elman神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Table 2 Predict result of 2014’ not normalized Elman ANN

圖2　2014年未歸一化Elman神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Fig. 2 Predict result of 2014’ not normalized Elman ANN

2.2[0, 1]區(qū)間歸一化Elman神經(jīng)網(wǎng)絡(luò)

由表3可知，4─9月份的預(yù)測誤差均小于16%，8月份的預(yù)測誤差僅為1.54，誤差百分比為5.31%，10月預(yù)測值和實測值的差值較大，預(yù)測趨勢和實際變化趨勢基本吻合。由圖3可知，預(yù)測值變化趨勢接近實測值，經(jīng)回歸分析可知，擬合度為 0.82（P=0.0047）。2014年平均誤差mse為4.51。

表3　2014年[0, 1]區(qū)間歸一化Elman神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Table 3 Predict result of 2014’ [0, 1] normalized Elman ANN

圖3　2014年[0, 1]歸一化Elman神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Fig. 3 Predict result of 2014’ [0, 1] normalized Elman ANN

表4　2014年[-1, 1]區(qū)間歸一化Elman神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Table 4 Predict result of 2014’ [-1, 1] normalized Elman ANN

2.3[-1, 1]區(qū)間歸一化結(jié)果Elman神經(jīng)網(wǎng)絡(luò)

由表4可知，4─9月份的預(yù)測誤差均小于20%，7月份的預(yù)測誤差為1.83，誤差百分比為5.72%，10月預(yù)測值和實測值的差值較大，預(yù)測值變化趨勢和實測值基本吻合。由圖4可知，預(yù)測值變化趨勢接近實測值，經(jīng)回歸分析可知，擬合度為 0.685（P=0.0216）。2014年平均誤差mse為5.03。

圖4　[-1, 1]歸一化Elman神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果Fig. 4 Predict result of 2014’ [-1, 1] normalized Elman ANN

2.43種處理方式分析

由表5可知，3種歸一化處理建立的神經(jīng)網(wǎng)絡(luò)模型的預(yù)測值、誤差、誤差百分比差異不顯著；通過比較神經(jīng)網(wǎng)絡(luò)性能評價指標(biāo)（擬合度 R2和平均誤差mse）發(fā)現(xiàn)，，，均以[0, 1]歸一化建立的模型最小，因此[0, 1]歸一化建立的模型較為理想。

表5　不同歸一化預(yù)測模型Table 5 Predict model of different normalized

3　討論

3.1種群數(shù)量變動規(guī)律

種群數(shù)量變動規(guī)律的研究始終是種群生態(tài)學(xué)研究的核心內(nèi)容，有關(guān)嚙齒動物種群數(shù)量變動規(guī)律，國內(nèi)外取得了較好的研究結(jié)果（房繼明，1989，1994；楊春文，1996；陳安國等，1981；董維惠等，1993，1998；姜運良等，1994；侯希賢等，2000a，2000b；楊再學(xué)等，2013；Ostfeld et al.，2000；Adamík et al.，2008；Letnic et al.，2011；Meserve et al.，2011；Leanne et al.，2014；Georgy，2014）。小型獸類種群數(shù)量一般存在季節(jié)消長規(guī)律，長期定點研究是最基本的研究方法（武文華等，2007；丁平等，1989；楊澤龍等，2005）。楊再學(xué)等（2015）經(jīng)過30多年監(jiān)測，發(fā)現(xiàn)黑線姬鼠種群數(shù)量的消長變化規(guī)律發(fā)生了改變，由前峰型變?yōu)榱撕蠓逍停毁F州地區(qū)黑線姬鼠在1年內(nèi)出現(xiàn)次數(shù)不等、數(shù)量高峰出現(xiàn)時間不同的現(xiàn)象（楊再學(xué)等，2015），黑腹絨鼠種群數(shù)量全年出現(xiàn)兩個數(shù)量高峰，不同月份變異較大（楊再學(xué)等，2013）。戴應(yīng)貴（2001）發(fā)現(xiàn)四川邛崍縣平原農(nóng)田嚙齒動物結(jié)合群落密度在84個月中有8個月有明顯的峰值，表現(xiàn)出不規(guī)律的交替升降動態(tài)變化，大足鼠種群與結(jié)合種群密度具有顯著相關(guān)性，與結(jié)合種群密度具有相似的動態(tài)變化規(guī)律。降雨促使食物資源增加，導(dǎo)致以莖葉和種子為食物的嚙齒動物數(shù)量增加，不同種類存在不同的延遲效應(yīng)（Letnic et al.，2011；Ostfeld et al.，2000；Georgy，2014）。干旱到來時，數(shù)量下降很快（Leanne et al.，2014）。Adamík et al.（2008）發(fā)現(xiàn)溫度對鹿鼠數(shù)量有很大的影響。不同生命跨度的物種在降雨后幾周到一年出現(xiàn)數(shù)量高峰（Letnic et al.，2011；Meserve et al.，2011）。影響嚙齒動物群落數(shù)量變化的因素包括氣象因子、食物資源（植物）、干擾等外部因子和嚙齒動物本身生命跨度、繁殖、躲避天敵、競爭等內(nèi)部因子共同發(fā)揮作用。在群落水平上研究嚙齒動物組成物種數(shù)量變化規(guī)律，這種規(guī)律呈現(xiàn)出復(fù)雜性和可變性，增加了預(yù)測的難度。造成鼠害的鼠種往往不是一種嚙齒動物，而是幾種嚙齒動物的結(jié)合種群密度，實際生產(chǎn)中對草地（農(nóng)田、林地）造成損失的往往是嚙齒動物群落。本研究旨在根據(jù)歷史嚙齒動物群落數(shù)量進行學(xué)習(xí)預(yù)測。如上所述，以往的研究大部分針對1種嚙齒動物從不同角度進行研究，很少一部分研究針對嚙齒動物群落數(shù)量進行研究。嚙齒動物群落數(shù)量受到氣候、植物、食物、天敵以及種間、種內(nèi)關(guān)系綜合作用達(dá)到暫時平衡狀態(tài)，這種平衡狀態(tài)隨時間變化而表現(xiàn)出復(fù)雜的變化規(guī)律。

經(jīng)過歸一化處理的數(shù)據(jù)能夠促進神經(jīng)網(wǎng)絡(luò)的收斂速度和預(yù)測精度，歸一化方法因研究對象不同而存在差異，本研究表明經(jīng)過[0, 1]歸一化能夠提高嚙齒動物物種組成數(shù)量的預(yù)測精度，更加準(zhǔn)確模擬群落數(shù)量動態(tài)規(guī)律；對未來群落數(shù)量進行了預(yù)測，取得了較好的結(jié)果（R2>0.8）。

3.2種群數(shù)量預(yù)測方法

嚙齒動物種群數(shù)量預(yù)測主要是對等級預(yù)測和種群數(shù)量（密度）進行預(yù)測。武文華等（2007）應(yīng)用馬爾科夫鏈各階轉(zhuǎn)移概率矩陣，馮志勇等（2000）利用平滑指數(shù)法，姜運良等（1994）通過逐步回歸分析，準(zhǔn)確預(yù)測不同鼠種捕獲率等級和發(fā)生等級；王玉志等（1997）運用灰色系統(tǒng)，以種群數(shù)量為主行為序列，其它因子為相關(guān)因子序列，建立了動態(tài)模型。侯希賢等（2000a；2000b）應(yīng)用捕獲率、繁殖指數(shù)建立了線性回歸模型分別對小毛足鼠和子午沙鼠進行了短、中、長期預(yù)測。不同科研工作者運用不同的方法對嚙齒動物數(shù)量進行了預(yù)測，取得了良好的效果。本文運用捕獲量建立了Elman神經(jīng)網(wǎng)絡(luò)，實現(xiàn)了連續(xù)預(yù)測，并取得了較好的效果，最小誤差僅為5.13%。

影響嚙齒動物種群數(shù)量的內(nèi)部因素和外部因素眾多，包括氣象、植物、干擾以及不同物種繁殖能力、壽命、對環(huán)境的適應(yīng)能力，但究竟是內(nèi)部因素還是外部因素起主要作用，以往的研究得出不同的結(jié)論，實際上是內(nèi)部因素和外部因素共同作用影響了種群和群落的數(shù)量變化，因此使預(yù)測變得十分困難。Letnic et al.（2011）認(rèn)為預(yù)測的困難主要是在什么時間、什么地點出現(xiàn)數(shù)量爆發(fā)的把握上。侯希賢等（2000a）認(rèn)為選用諸如捕獲率、繁殖指數(shù)等這些指標(biāo)均未體現(xiàn)外部因素的具體數(shù)值，但均受到了外部因素的影響，這種影響受種群繁殖、生存策略等的影響，選用這些指標(biāo)使模型簡單，便于推廣。

從分析對象來看，神經(jīng)網(wǎng)絡(luò)的優(yōu)勢在于分析背景知識模糊的、相互影響關(guān)系不清楚的，即存在頻繁干擾的大通量信息；從預(yù)測結(jié)果來看，神經(jīng)網(wǎng)絡(luò)的特長在于搜索相對精度較低而比較滿意的結(jié)果，放棄高精度但難以實現(xiàn)的結(jié)果，從而大大增加了滿意結(jié)果的效率以及解決實際問題的能力。組成神經(jīng)網(wǎng)絡(luò)的基本單元及神經(jīng)元通過從臨近的大量其它基本單元接受信號及同步傳遞信號的相互作用機制來實現(xiàn)其信息處理功能。相互獨立的神經(jīng)元通過物理聯(lián)系實現(xiàn)存儲功能，連接權(quán)值的動態(tài)變化實現(xiàn)了網(wǎng)絡(luò)的學(xué)習(xí)和識別，若輸入與輸出間映射關(guān)系為F，則神經(jīng)網(wǎng)絡(luò)通過線性轉(zhuǎn)換實現(xiàn) F的簡單線性函數(shù)復(fù)合的最近逼近，實現(xiàn)樣本的輸入與輸出。本文正是根據(jù)這樣的思路，通過比較不同歸一化方法建立Elman神經(jīng)網(wǎng)絡(luò)，預(yù)測嚙齒動物數(shù)量，過程簡單，結(jié)果較為滿意。本文為進一步研究群落格局變化、物種的消失與更替奠定了一定基礎(chǔ)。Vellend（2010）指出生態(tài)學(xué)家面臨著清晰闡述群落格局和過程的連續(xù)動態(tài)過程，這種變化過程是下一步研究的重點內(nèi)容和方向。

4　結(jié)論

Elman神經(jīng)網(wǎng)絡(luò)模型憑借著其自身輸入與輸出高度的非線性映射的優(yōu)勢，為基于時間序列的嚙齒動物種群數(shù)量的預(yù)測提出了新方法，相對于傳統(tǒng)的預(yù)測方法，Elman神經(jīng)網(wǎng)絡(luò)模型基于歷史數(shù)據(jù)，反應(yīng)復(fù)雜環(huán)境因子對種群數(shù)量的綜合作用，能夠較好地預(yù)測未來的發(fā)展趨勢。

本文通過對比不同歸一化 Elman神經(jīng)網(wǎng)絡(luò)模型，結(jié)果表明[0, 1]歸一化處理建立的神經(jīng)網(wǎng)絡(luò)在嚙齒動物種群數(shù)量方面的預(yù)測，取得了較好的擬合效果（R2=0.82，P=0.0065），平均誤差較?。╩se=4.51），4─9月份的預(yù)測準(zhǔn)確率超過了80%，可以指導(dǎo)實際應(yīng)用。

致謝：本文使用了近十年的野外數(shù)據(jù)，飽含著內(nèi)蒙古阿拉善盟草原工作站工作人員、歷屆師兄師姐師弟師妹們的辛苦工作，在此對你們表示衷心的感謝！

ADAMíK P, KRáL M. 2008. Climate- and resource-driven long-term changes in dormice populations negatively affect hole-nesting songbirds [J]. Journal of Zoology, 275(3): 209-215.

CHITTY D. 1960. Population processes in the vole and their reference to general theory [J]. Canadian Journal of Zoology, 38(1): 99-113.

CHRISTIAN J J, DAVIS D E. 1964. Endocrines,behavior and population [J]. Science, 146(3651): 1550-1560.

GEORGY SHENBROT. 2014. Population and community dynamics and habitat selection of rodents in complex desert landscapes [J]. Mammalia, 78(1): 1-10.

KREBS, C J, GAINES M S, KELLER B L, et al. 1973. Population cycles in small rodents [J]. Science, 179: 35-41.

LEANNE M HEISLER, CHRISTOPHER M SOMERS, RAY G POULIN. 2014. Rodent populations on the northern Great Plains respond to weather variation at a landscape scale [J]. Journal of Mammalogy, 95(1): 82-90.

LETNIC M, STORY P, STORY G, et al. 2011. Resource pulses, switching trophic control, and the dynamics of small mammal assemblages in arid Australia [J]. Journal of Mammalogy, 92(6): 1210-1222.

LIMA M, KEYMER J, JAKSIC F. 1999. El Nino-Southern Oscillation-Drived rainfall variability and delayed density dependence cause rodent outbreaks in Western South America: Linking demography and population dynamics [J]. The America Naturalist, 153(5): 476-491.

MESERVE P L, KELT D A, PREVITALI M A, et al. 2011. Global climate change and small mammal populations in north-central [J]. Journal of Mammalogy, 92(6): 1223-1235.

OSTFELD R S, KEESING F. 2000. Pulsed resources and community dynamics of consumers in terrestrial ecosystems [J]. Trends in Ecology and Evolution, 15(6): 232-237.

VELLEND M. 2010. Conceptual synthesis in community ecology [J]. The Quarterly Review of Biology, 85(5): 183-206.

WYNNE-EDWARDS V C. 1962. Animal Dispersion in Relation to Social Behavior [M]. New York: Hafner.

陳安國, 朱盛侃, 李春秋, 等. 1981. 新疆北部農(nóng)業(yè)區(qū)害鼠的研究（六）小家鼠種群數(shù)量消長同氣候的關(guān)系[R]//滅鼠和鼠類生物學(xué)研究報告[M]. 北京: 科學(xué)出版社: 69-93.

陳明. 2013. MATLAB神經(jīng)網(wǎng)絡(luò)原理與實例精解[M]. 北京: 清華大學(xué)出版社: 278-307.

崔東文. 2013. 改進Elman神經(jīng)網(wǎng)絡(luò)在徑流預(yù)測中的應(yīng)用[J]. 水利水運工程學(xué)報, 4(2): 71-77.

戴應(yīng)貴. 2001. 川西平原農(nóng)田嚙齒動物群落動態(tài): 年間變動和季節(jié)變動[J]. 獸類學(xué)報, 21(1): 23-34.

丁平, 鮑毅新, 諸葛陽, 等. 1989. 蕭山圍墾農(nóng)區(qū)小型獸類種群動態(tài)的研究[J]. 獸類學(xué)報, 9(2): 35-42.

董維惠, 侯希賢, 林小泉, 等. 1993. 黑線倉鼠種群數(shù)量動態(tài)預(yù)測研究[J].生態(tài)學(xué)報, 13(4): 300-304.

董維惠, 侯希賢, 周延林, 等. 1998. 黑線毛足鼠種群數(shù)量動態(tài)及預(yù)測研究[J]. 草地學(xué)報, 6(34): 207-211.

房繼明, 孫儒泳. 1989. 布氏田鼠種群數(shù)量的季節(jié)動態(tài)與鼠洞的關(guān)系[J].獸類學(xué)報, 9(3): 202-209.

房繼明, 孫儒泳. 1994. 布氏田鼠數(shù)量和空間分布的年際動態(tài)及周期性初步分析[J]. 動物學(xué)雜志, 29(6): 35-37.

馮志勇, 黃秀清, 顏世祥, 等. 2000. 黃毛鼠種群數(shù)量中、長期測報模型的研究[J]. 廣東農(nóng)業(yè)科學(xué), (5): 47-49.

高瑋, 鄭穎人. 2001. 巖土工程位移預(yù)測神經(jīng)網(wǎng)絡(luò)建模的幾個問題[J].地下空間, 21(5): 369-374.

何淼, 林繼球, 翁文英, 等. 1996. 板齒鼠種群數(shù)量中長期預(yù)測的時間序列模型[J].獸類學(xué)報, 16(4): 297-301.

侯希賢, 董維惠, 楊玉平, 等. 2013. 鄂爾多斯沙地草場小毛足鼠種群數(shù)量動態(tài)分析[J]. 中國媒介生物學(xué)及控制雜志, 14(3): 177-180.

侯希賢, 董維惠, 周延林, 等. 2000a. 鄂爾多斯沙地草場小毛足鼠種群數(shù)量動態(tài)及預(yù)測[J]. 中國媒介生物學(xué)及控制雜志, 11(1): 7-10.

侯希賢, 董維惠, 周延林, 等. 2000b. 子午沙鼠種群數(shù)量動態(tài)及預(yù)測[J].生態(tài)學(xué)報, 20(4): 711-714.

姜運良, 盧浩泉, 李玉春, 等. 1994. 山東陽谷縣黑線倉鼠種群數(shù)量預(yù)測預(yù)報[J].獸類學(xué)報, 14(3): 195-202.

劉冰, 郭海霞. 2014. MATLAB神經(jīng)網(wǎng)絡(luò)超級學(xué)習(xí)手冊[M]. 北京: 人民郵電出版社: 242-252.

劉天舒. 2011. BP神經(jīng)網(wǎng)絡(luò)的改進研究及應(yīng)用[D]. 哈爾濱: 東北農(nóng)業(yè)大學(xué): 20-32.

彭昱忠，王謙，元昌安，等. 2015. 數(shù)據(jù)挖掘技術(shù)在氣象預(yù)報研究中的應(yīng)用[J]. 干旱氣象, 33(1): 19-27.

邵月紅, 林炳章, 葉金印, 等. 2014. 基于 Elman 動態(tài)神經(jīng)網(wǎng)絡(luò)的降雨—徑流模擬研究[J]. 大氣科學(xué)學(xué)報, 37(2): 223-228.

孫儒泳. 2001. 動物生態(tài)學(xué)原理（第三版）[M]. 北京: 北京師范大學(xué)出版社: 122-139.

王利民, 周延林, 董維惠, 等. 1998. 庫布其沙地黑線倉鼠種群數(shù)量動態(tài)的研究[J]. 內(nèi)蒙古大學(xué)學(xué)報, 29(1): 109-113.

王小川, 史峰, 郁磊, 等. 2013. MATLAB神經(jīng)網(wǎng)絡(luò)43個案例分析[M].北京: 北京航空航天大學(xué)出版社: 196-200.

王勇, 陳安國, 郭聰, 等. 1997. 洞庭湖稻區(qū)黑線姬鼠種群數(shù)量預(yù)測[J].獸類學(xué)報, 17(2): 125-130.

王玉志, 盧浩泉, 陳安, 等. 1997. 灰色系統(tǒng)在華北平原旱作區(qū)黑線倉鼠種群數(shù)量預(yù)測預(yù)報中的應(yīng)用[J]. 動物學(xué)報, 43(增刊): 107-112.

王中凱, 梁秀娟, 肖長來, 等. 2013. Elman神經(jīng)網(wǎng)絡(luò)在平原區(qū)降水入滲補給預(yù)測中的應(yīng)用[J]. 節(jié)水灌溉, (7): 42-44.

吳萍, 遲道才, 王殿武, 等. 2007. Elman反饋型神經(jīng)網(wǎng)絡(luò)模型在灌溉用水量預(yù)測中的應(yīng)用[J]. 沈陽農(nóng)業(yè)大學(xué)學(xué)報, 38(4): 586-589.

武文華, 付和平, 武曉東, 等. 2007. 應(yīng)用馬爾可夫鏈模型預(yù)測長爪沙鼠和黑線倉鼠種群數(shù)量[J]. 動物學(xué)雜志, 42(6): 69-78.

徐興梅，曹麗英. 2014. 基于粗糙集和BP神經(jīng)網(wǎng)絡(luò)的糧食產(chǎn)量預(yù)測研究[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報, 45(10): 55-100.

楊再學(xué), 金星, 郭永旺, 等. 2010. 高山姬鼠種群數(shù)量動態(tài)及預(yù)測預(yù)報模型[J]. 生態(tài)學(xué)報, 30(13): 3545-3552.

楊春文，張春美，張廣臣, 等.1996. 長白山地林區(qū)棕背?種群數(shù)量分布及變動的研究[J]. 森林病蟲通訊, (3): 9-11.

楊再學(xué), 金星, 郭永旺, 等. 2015. 貴州省不同地區(qū)黑線姬鼠種群數(shù)量動態(tài)分析[J]. 山地農(nóng)業(yè)生物學(xué)報, 34(1): 13-17.

楊再學(xué), 雷邦海, 金星, 等. 2013. 凱里市黑腹絨鼠種群數(shù)量變動規(guī)律[J].中國農(nóng)學(xué)通報, 29(36): 378-381.

楊澤龍, 傅和平, 武曉東, 等. 2005. 內(nèi)蒙古阿拉善荒漠嚙齒動物群落格局的干擾效應(yīng)[J]. 內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報, 26(2): 28-34.

殷高方, 張玉鈞, 胡麗, 等. 2012. BP神經(jīng)網(wǎng)絡(luò)在水華短期預(yù)測中的應(yīng)用[J]. 北京理工大學(xué)學(xué)報, 32(6): 655-660.

張志彬, 王祖望. 1998. ENSO現(xiàn)象與生物災(zāi)害[J]. 中國科學(xué)院院刊, 13(1): 34-38.

張志彬. 1999. 厄爾尼諾與大尺度、超長期生物災(zāi)害預(yù)警[J]. 中國科學(xué)院院刊, (1): 34-36.

朱盛侃, 陳安國. 1993. 小家鼠生態(tài)特征與預(yù)測[M]. 北京: 科學(xué)出版社: 197-250.

鄒波, 張慧娣, 魏明峰, 等. 2007. 晉西林區(qū)地面害鼠數(shù)量的季節(jié)變化與防治技術(shù)研究[J]. 中國植保導(dǎo)刊, 27(12): 8-11.

祖元剛，趙則海，于景華，等. 2004. 非線性生態(tài)模型[M]. 北京: 科學(xué)出版社: 530-572.

Prediction of the Number of Rodent Community Composition Species Based on Elman Neural Network in Alasan Desert

LU Zhihong1,2, WU Xiaodong1,2*, GUO Libiao3, FU Heping1,2, MAN Duhu1,2, YUE Chuang1,2, CHAI Xiangxian1,2, BAO Daerhan1,2, YANG Suwen1,2, YE Lina1,2, LI Yanni1,2
1. College of Ecology and Environmental Science, Inner Mongolia Agricultural University, Hohhot 010018, China; 2. Key Laboratory of Prataculture Grassland Resources, Ministry of Education, Huhhot 010018, China; 3. Grassland Research Institute, Chinese Academy of Agricultural Sciences, Huhhot 010018, China

The fluctuation and prediction of population is one of important research contents in ecology, however, it is necessary to explore new approaches. In this study, a novel method was used in rodent population prediction by neural network technology. The neural network of Elmam has the function of memory, which can simulate and forecast the quantity of species in rodent by adding delay-units in feedforward networks. Based on the eastern edge of Tengger desert as the study area, with the rodent population as the research object, by use of mark recapture method to continuously looked into 2006─2014 (Apr-Oct) and build Elman neural network forecasting model, data between 2006─2013 were used to build training network, data of 2014 were used for test, the mean error and fitting degree of three processing methods were compared and predicted the number of Alashan desert rodent population dynamics. Results showed that: (1) Average error of the prediction results without normalization mse was 5.30 with Minimum error of 1.52%, fitting degree 0.80. (2) After [0, 1] normalization, Average error of the prediction results mse was 4.51 with Minimum error of 1.54%, fitting degree 0.82. (3) After [-1, 1] normalization, Average error of the prediction results mse was 5.03 with Minimum error of 1.63%, fitting degree 0.69. (4) The difference of Elman neural network model was not significant after three kinds of normalized treatment. While comparison of mean error and fitting degree, the article suggested that by use of [0, 1] normalization to establish Elman neural network can better predicted the animal composition and species dynamic trend. Thus providing theoretical basis for guiding and preventing local rodent infestation

elman neural network; alasan desert; rodent; mark recapture method

10.16258/j.cnki.1674-5906.2015.12.008

S714.6; X171.1

A

1674-5906（2015）12-1976-07

國家自然科學(xué)基金資助項目（30760044；31160096）；公益性行業(yè)科研專項經(jīng)費項目（201203041）

盧志宏（1982年生），男（蒙古族），博士研究生，主要從事嚙齒動物生態(tài)學(xué)與草地保護研究。E-mail: luzhihong305@163.com *通信作者：武曉東（1959年生），男，教授，博士生導(dǎo)師，主要從事嚙齒動物生態(tài)學(xué)與草地保護研究。Email: wuxiaodong_hgb@163.com

2015-11-02

引用格式：盧志宏, 武曉東, 郭利彪, 付和平, 滿都呼, 岳闖, 柴享賢, 包達(dá)爾罕, 楊素文, 葉麗娜, 李燕妮. 基于Elman神經(jīng)網(wǎng)絡(luò)的阿拉善荒漠嚙齒動物群落組成物種數(shù)量預(yù)測研究[J]. 生態(tài)環(huán)境學(xué)報, 2015, 24(12): 1976-1982.