地力夏提·艾木熱拉，丁建麗*，穆艾塔爾·賽地,3，米熱古力·艾尼瓦爾，鄒 杰

(1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院智慧城市與環(huán)境建模自治區(qū)普通高校重點(diǎn)實(shí)驗(yàn)室，新疆 烏魯木齊 830046；2.綠洲生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室，新疆 烏魯木齊 830046；3.新疆水利水電科學(xué)研究院，新疆 烏魯木齊 830049)

【研究意義】土壤鹽漬化是指在特定氣候、地質(zhì)及土壤質(zhì)地和人為用水灌溉不當(dāng)?shù)雀鞣N因素的綜合作用下，所引起的土地質(zhì)量退化過程，一般在干旱、半干旱地區(qū)出現(xiàn)[1]。近幾年，人口不斷增加，氣候和土地利用的變化頻繁，使土壤及其結(jié)構(gòu)性退化正加速危害土壤質(zhì)量和健康[2]。土壤鹽漬化作為土地荒漠化形式之一，主要分布于我國的西部干旱和半干旱地區(qū)，是西部農(nóng)業(yè)發(fā)展及綠洲生態(tài)環(huán)境穩(wěn)定主要障礙之一[3-5]。

新疆鹽漬土面積廣泛，受到鹽漬化災(zāi)害比較嚴(yán)重，大量的鹽漬土造成土壤生產(chǎn)力下降、水土失衡，不僅影響著新疆在我國糧食和棉花生產(chǎn)基地的戰(zhàn)略地位，同時(shí)對(duì)干旱區(qū)域大尺度生態(tài)環(huán)境安全和綠洲經(jīng)濟(jì)可持續(xù)發(fā)展造成阻礙。因此，土壤鹽漬化災(zāi)害評(píng)價(jià)和預(yù)測(cè)研究對(duì)研究區(qū)合理利用土地資源、提高糧食和棉花安全生產(chǎn)具有重要的科學(xué)意義?！厩叭搜芯窟M(jìn)展】隨著人工智能系統(tǒng)的不斷發(fā)展，許多國內(nèi)外學(xué)者在利用人工智能技術(shù)針對(duì)數(shù)據(jù)預(yù)測(cè)方面進(jìn)行了大量的研究。其中，Vijendra等[6]建立了一套農(nóng)業(yè)干旱預(yù)警模型，對(duì)印度拉賈斯坦地區(qū)的干旱農(nóng)業(yè)耕地進(jìn)行了初步研究。Bijanzadeh等[7]通過模糊隸屬函數(shù)的土壤與環(huán)境之間的關(guān)系來預(yù)測(cè)土壤肥力，然后采用模糊層次分析法(FAHP)對(duì)小麥的發(fā)育階段及其與土壤鹽分的關(guān)系進(jìn)行了相關(guān)研究。買買提·沙吾提等[8]利用BP神經(jīng)網(wǎng)絡(luò)的3種算法，選取7個(gè)影響因子作為輸入變量，將土壤脫鹽量作為輸出變量，建立了BP神經(jīng)網(wǎng)絡(luò)土壤鹽漬化預(yù)測(cè)模型，并分別對(duì)3種BP網(wǎng)絡(luò)算法進(jìn)行了對(duì)比分析研究。謝姆斯葉·艾尼瓦爾等[9]建立了預(yù)測(cè)表層土壤鹽分的BP神經(jīng)網(wǎng)絡(luò)模型和RBF神經(jīng)網(wǎng)絡(luò)模型，發(fā)現(xiàn)BP網(wǎng)絡(luò)的預(yù)測(cè)誤差略低于RBF網(wǎng)絡(luò)。丁建麗等[10]采用BP神經(jīng)網(wǎng)絡(luò)結(jié)合Adaboost算法，利用遙感技術(shù)結(jié)合地下水埋深、鹽分指數(shù)、海拔、歸一化干旱指數(shù)4個(gè)影響因子對(duì)研究區(qū)的土壤鹽漬化災(zāi)害程度進(jìn)行了預(yù)測(cè)評(píng)價(jià)?！颈狙芯壳腥朦c(diǎn)】近幾年把人工智能技術(shù)的數(shù)據(jù)預(yù)測(cè)思想應(yīng)用到農(nóng)業(yè)、生態(tài)環(huán)境、氣候等領(lǐng)域取得較大的成功。模糊神經(jīng)網(wǎng)絡(luò)(Fuzzy Neural Network，F(xiàn)NN)作為這幾年發(fā)展較快的一種新興網(wǎng)絡(luò)模型，具有收斂速度快等特點(diǎn)，逐漸受到數(shù)據(jù)分析工作者的青睞。【擬解決的關(guān)鍵問題】本文以渭干河-庫車河三角洲綠洲(以下簡(jiǎn)稱渭庫綠洲)為研究區(qū)，根據(jù)研究區(qū)實(shí)際情況及野外考察所獲得的數(shù)據(jù)為基礎(chǔ)，并參考前人的研究成果，針對(duì)土壤鹽分嘗試構(gòu)建模糊邏輯與神經(jīng)網(wǎng)絡(luò)有機(jī)結(jié)合的表層土壤(0～10 cm)鹽分預(yù)測(cè)模型。

1 材料與方法

1.1 研究區(qū)概況

渭干河-庫車河流域三角洲綠洲(圖1)位于塔里木盆地北緣，天山南麓，塔克拉瑪干沙漠以北，地勢(shì)北高南低，自西北向東南傾斜。渭庫綠洲氣候?qū)儆诖箨懶耘瘻貛Ц珊禋夂?，多年平均氣?0.5～11.4 ℃，極端最高氣溫40.8 ℃，極端最低氣溫-27.8 ℃，日照2888.7 h，全年無霜期209 d，多年平均蒸發(fā)量超過2000 mm，多年平均降水量55.45 mm，降水少，蒸發(fā)強(qiáng)烈，干燥度指數(shù)為44.3，屬于典型的干旱與極端干旱地區(qū)。該地區(qū)植被覆蓋稀少，主要以檉柳、鹽節(jié)木、鹽穗木、花花柴等為主，地下水位高、礦化度大，土層構(gòu)成物顆粒細(xì)，透水性差，綠洲及其外圍鹽類沉積規(guī)模大，鹽漬化普遍存在[11-13]。農(nóng)牧業(yè)是該區(qū)的主要經(jīng)濟(jì)來源，隨著人類不合理灌溉的逐漸擴(kuò)大，地下水位增高，鹽漬化程度逐漸加劇，對(duì)當(dāng)?shù)氐霓r(nóng)業(yè)生產(chǎn)及生態(tài)系統(tǒng)造成嚴(yán)重的危害(表1)。

1.2 數(shù)據(jù)采集與處理

通過遙感影像上的光譜差異選取具有一定代表性的鹽漬地，并記錄其坐標(biāo)，使用GPS精確導(dǎo)航定位，分別從綠洲內(nèi)部和綠洲荒漠交錯(cuò)帶選擇39個(gè)采樣點(diǎn)。土壤采集需滿足室內(nèi)試驗(yàn)要求，室內(nèi)試驗(yàn)主要測(cè)定土壤pH值、總?cè)芙夤腆w(TDS)、電導(dǎo)率以及含鹽量等指標(biāo)，在樣點(diǎn)區(qū)采集土壤樣品為理化指標(biāo)分析和研究鹽漬地建模提供科學(xué)有效的數(shù)據(jù)，采樣時(shí)間為2015年7月。此外，地下水埋深和地下水礦化度數(shù)據(jù)是渭干河流域管理局提供的38個(gè)觀測(cè)井的年平均數(shù)據(jù)。高程和坡度數(shù)據(jù)是從Aster 90 m空間分辨率的DEM數(shù)據(jù)中獲取。

1.3 研究方法

1.3.1 灰色關(guān)聯(lián)分析法 灰色關(guān)聯(lián)度分析法(Gray Relational Analysis，GRA)是由鄧聚龍教授1982年首先提出的，其基本思想是根據(jù)序列曲線幾何形狀的相似度來判斷其關(guān)聯(lián)程度是否緊密的一種多因素統(tǒng)計(jì)分析法。它可以定量的描述事物或因素之間的大小、方向和次序等變化情況的相對(duì)性。若事物或因素的變化形勢(shì)大致相同，可以認(rèn)為它們之間的關(guān)聯(lián)程度較大。目前，灰色關(guān)聯(lián)度分析法在農(nóng)業(yè)、環(huán)保、教育、地理、水文以及氣象等領(lǐng)域的應(yīng)用十分廣泛，尤其是通過已知信息研究和預(yù)測(cè)未知領(lǐng)域達(dá)到了解整個(gè)系統(tǒng)的目的。利用關(guān)聯(lián)度分析方法先計(jì)算第i個(gè)評(píng)價(jià)對(duì)象的第k個(gè)最優(yōu)指標(biāo)的關(guān)聯(lián)系數(shù)，即：

圖1 研究區(qū)Fig.1 The study area

(1)

1.3.2 T-S模糊神經(jīng)網(wǎng)絡(luò)模型 模糊神經(jīng)網(wǎng)絡(luò)(Fuzzy Neural Network，F(xiàn)NN)技術(shù)是模糊邏輯和神經(jīng)網(wǎng)絡(luò)學(xué)習(xí)逼近的特點(diǎn)相結(jié)合的，能自動(dòng)更新和不斷修正模糊集隸屬度函數(shù)的一種處理模糊信息的神經(jīng)網(wǎng)絡(luò)系統(tǒng)，它匯聚了模糊系統(tǒng)和神經(jīng)網(wǎng)絡(luò)的優(yōu)勢(shì)，克服了它們各方面的缺點(diǎn)[15-17]。近幾年來在數(shù)據(jù)預(yù)測(cè)方面被廣泛應(yīng)用，在生態(tài)環(huán)境方面主要在災(zāi)害預(yù)測(cè)、危險(xiǎn)度評(píng)價(jià)、地下水位以及水質(zhì)預(yù)測(cè)等方面應(yīng)用較多。

T-S模型是1985年由2名日本科學(xué)家Takagi和Sugeno提出的一種模糊推理模型[18]。T-S模糊模型利用“if-then”規(guī)則方式來表示，在規(guī)則方式為Ri的情況下，模糊推理方式用下式來描述：

(2)

為模糊規(guī)則獲得的輸出，輸入部分if是模糊的，輸出部分then是確定的，說明模糊推理可以描述輸出為輸入的線性結(jié)合。

T-S模糊神經(jīng)網(wǎng)絡(luò)模型的結(jié)構(gòu)可以分為4層，分別是輸入層、模糊化層、模糊規(guī)則計(jì)算層以及輸出層(圖2)。輸入層和輸入變量x1接連，輸入變量和節(jié)點(diǎn)數(shù)的維數(shù)一樣。模糊化層是利用隸屬度函數(shù)對(duì)輸入層進(jìn)行模糊化求得模糊隸屬度μ。模糊規(guī)則計(jì)算是采用連乘算子來求得ω，其計(jì)算公式為：

(3)

根據(jù)模糊規(guī)則計(jì)算結(jié)果，利用公式(4)求得模糊神經(jīng)網(wǎng)絡(luò)模型的輸出值yi：

(4)

最后進(jìn)行誤差計(jì)算，其公式為：

(5)

其中，yd是預(yù)測(cè)輸出值，yc是實(shí)測(cè)值，e為預(yù)測(cè)值和實(shí)測(cè)值的誤差。

2 結(jié)果與分析

2.1 土壤鹽分預(yù)測(cè)指標(biāo)的獲取與分析

2.1.1 地下水埋深和地下水礦化度 地下水埋深和地下水礦化度與土壤鹽分有密切的關(guān)系，由于地下水位從高到低的回降過程中，土壤因蒸發(fā)而積鹽，如果地下水位回降過慢，土壤積鹽就越多，這就導(dǎo)致土壤鹽堿化的發(fā)生[19]。利用ArcGIS軟件的地統(tǒng)計(jì)分析模塊，將2015年渭庫綠洲38個(gè)觀測(cè)井作為地下水埋深及地下水礦化度數(shù)據(jù)。獲取的地下水埋深和地下水礦化度數(shù)據(jù)基本上服從正太分布，符合半方差函數(shù)的計(jì)算需求。從地下水埋深和地下水礦化度的空間變異參數(shù)(表1)看出，相關(guān)系數(shù)分別為0.843和0.825，說明選擇球面模型和高斯模型進(jìn)行普通克拉格插值法，可以較好的反映地下水埋深和地下水礦化度的空間分布狀況(圖3)。

圖2 T-S模糊神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.2 Fuzzy Neural Network Structure

表1 地下水埋深和地下水礦化度的空間變異參數(shù)

2.1.2 高程和坡度分析 土壤鹽漬化的形成在一定程度上取決于地形起伏和地貌類型，會(huì)直接影響地下水與地表水的徑流運(yùn)動(dòng)，導(dǎo)致鹽分隨水流動(dòng)進(jìn)而出現(xiàn)土壤鹽分的再分布。坡度和地形起伏對(duì)于輕、中度鹽漬化的分布影響較小，對(duì)重度鹽漬化的分布影響較大[20]。利用ArcGIS軟件的三維分析工具分別提取各采樣點(diǎn)的高程和坡度信息。由圖4可見，研究區(qū)東部綠洲荒漠交錯(cuò)帶海拔在945～980 m內(nèi)重度鹽漬化分布較廣泛；坡度大于1.3°時(shí)，基本上沒有重度鹽漬化分布。

2.1.3 土地利用狀況 土地利用是多學(xué)科的研究對(duì)象，其主要目的在于了解研究區(qū)土地利用狀況和對(duì)生態(tài)環(huán)境變化的預(yù)測(cè)以及在分析各種驅(qū)動(dòng)因子方面具有重要的研究意義。本文使用中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心的1∶10萬比例尺的土地利用現(xiàn)狀數(shù)據(jù)庫中獲取研究區(qū)土地利用狀況。如圖5所示，研究區(qū)土地利用類型由林地、耕地、水域、草地、居民地和未利用土地6個(gè)一級(jí)類型以及20個(gè)二級(jí)類型組成。野外考察的39個(gè)樣本點(diǎn)大部分分布于鹽堿地和中、低覆蓋草地區(qū)域，另外東北區(qū)域鹽堿地普遍存在。

2.2 各影響因子間的關(guān)聯(lián)度分析

土壤鹽漬化是在多種自然因素綜合作用下，對(duì)水鹽運(yùn)移產(chǎn)生影響的過程。在利用模糊神經(jīng)網(wǎng)絡(luò)模型對(duì)土壤鹽分進(jìn)行預(yù)測(cè)之前，首先要考慮各影響因子和土壤鹽分之間的關(guān)聯(lián)程度。首先對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，并將土壤鹽分設(shè)置為母序列，得出土壤鹽分與各影響因子之間的關(guān)聯(lián)程度。一般來說關(guān)聯(lián)程度在0.5～0.7屬于關(guān)聯(lián)度一般，0.70～0.85屬于關(guān)聯(lián)度較大，0.85～1.00表示關(guān)聯(lián)度大[21]。如表2所示，土壤鹽分與電導(dǎo)率的關(guān)聯(lián)系數(shù)最大，與坡度值的關(guān)聯(lián)系數(shù)最小，選取的7個(gè)影響因子與土壤鹽分的關(guān)聯(lián)度均大于0.5，說明所選擇的影響因子用來預(yù)測(cè)土壤鹽分是比較合適的。

圖3 研究區(qū)地下水埋深和地下水礦化度空間分布Fig.3 The Spatial distribution of the groundwater depth and salinity in the studied area

圖4 研究區(qū)高程 (a)和坡度分析 (b)Fig.4 The elevation (a) and slope analysis (b) in the studied area

圖5 研究區(qū)土地利用類型和土壤樣本分布Fig.5 Land use types of the study area and the distribution of soil samples

2.3 T-S模糊神經(jīng)網(wǎng)絡(luò)模型的建立與分析

本研究以土壤鹽分為研究對(duì)象，多年實(shí)地考察獲得的數(shù)據(jù)為基礎(chǔ)，把前人的研究經(jīng)驗(yàn)作為參考，采用T-S模糊神經(jīng)網(wǎng)絡(luò)模型來預(yù)測(cè)土壤鹽分。數(shù)據(jù)源是2015年7月采集39個(gè)采樣點(diǎn)土壤表層(0～10 cm)的實(shí)驗(yàn)數(shù)據(jù)，選用7個(gè)因素作為輸入樣本，土壤鹽分作為輸出樣本，對(duì)研究區(qū)土壤鹽分進(jìn)行T-S模糊神經(jīng)網(wǎng)絡(luò)模型的土壤含鹽量預(yù)測(cè)實(shí)驗(yàn)。

本文通過matlab軟件編程實(shí)現(xiàn)，以電導(dǎo)率、坡度值、高程、TDS、地下水礦化度、地下水位和pH值等影響土壤鹽漬化的7個(gè)因子作為輸入因子，土壤含鹽量為輸出因子，選擇32個(gè)采樣點(diǎn)的數(shù)據(jù)作為訓(xùn)練樣本，7個(gè)采樣點(diǎn)的數(shù)據(jù)作為檢驗(yàn)樣本，經(jīng)過歸一化處理后進(jìn)行模擬模型。輸入數(shù)據(jù)為7維，輸出數(shù)據(jù)為1維，隱含節(jié)點(diǎn)為25，神經(jīng)網(wǎng)絡(luò)的隱含節(jié)點(diǎn)用來表示隸屬度函數(shù)和模糊規(guī)則，選擇7組系數(shù)p0～p6進(jìn)行初始化模糊神經(jīng)網(wǎng)絡(luò)，模糊隸屬度函數(shù)的中心和寬度可隨機(jī)得到，網(wǎng)絡(luò)學(xué)習(xí)效率為0.05，慣性系數(shù)為0.005，最大迭代次數(shù)為200。

從表3可見，檢驗(yàn)樣本的平均相對(duì)誤差為13.092 %，最小誤差為0.875 %，最大相對(duì)誤差為41.733 %。出現(xiàn)誤差波動(dòng)的主要原因在于研究區(qū)鹽漬化程度變化比較大。模糊神經(jīng)網(wǎng)絡(luò)在模擬輕度鹽漬化和重度鹽漬化時(shí)，相對(duì)誤差較小，模擬效果較好；模擬非鹽漬地和中度鹽漬化時(shí)，相對(duì)誤差比較大。從圖8可以看出，基于T-S模型的模糊神經(jīng)網(wǎng)絡(luò)模型計(jì)算出來的預(yù)測(cè)精度較高，相關(guān)系數(shù)達(dá)到0.989，說明預(yù)測(cè)能力良好，可以用來預(yù)測(cè)土壤鹽漬化狀況。

表2 灰色關(guān)聯(lián)度分析結(jié)果

表3 檢驗(yàn)樣本誤差

圖8 實(shí)測(cè)值和預(yù)測(cè)值Fig.8 The measured and predicted values

3 討 論

土壤鹽漬化是多種影響因子共同作用形成的復(fù)雜非線性系統(tǒng)，影響模糊神經(jīng)網(wǎng)絡(luò)模擬土壤鹽分效果的原因主要是由于受實(shí)測(cè)數(shù)據(jù)的限制，只選擇了7個(gè)對(duì)土壤鹽漬化的影響因子，這就影響了土壤鹽分的模擬效果。其次是模糊神經(jīng)網(wǎng)絡(luò)本身的缺欠引起的，以往研究表明，模糊神經(jīng)網(wǎng)絡(luò)的自我調(diào)整能力和適應(yīng)能力較差，對(duì)于一些因素的抵抗性較弱，無法對(duì)系統(tǒng)進(jìn)行實(shí)時(shí)控制；另外，這種方法的實(shí)用性還有待提高，這就對(duì)模糊神經(jīng)網(wǎng)絡(luò)的應(yīng)用性形成了較大的限制。

對(duì)于土壤鹽漬化預(yù)測(cè)以及災(zāi)害評(píng)價(jià)而言，前人研究成果不盡相同。目前土壤鹽漬化預(yù)測(cè)主要是采用不同的神經(jīng)網(wǎng)絡(luò)方法來建立土壤鹽漬化預(yù)測(cè)模型，無論是哪種神經(jīng)網(wǎng)絡(luò)模型，對(duì)樣本訓(xùn)練的要求均較高。網(wǎng)絡(luò)訓(xùn)練是輸入和輸出變量之間的復(fù)雜規(guī)律，土壤鹽漬化預(yù)測(cè)模型效果的好壞主要是取決于輸入變量的數(shù)量和網(wǎng)絡(luò)訓(xùn)練的效果。

4 結(jié) 論

本文從土壤鹽漬化影響因子出發(fā)，結(jié)合土壤鹽漬化自然環(huán)境特征，選取了電導(dǎo)率、海拔高程、地下水埋深、地下水礦化度、TDS、pH值以及坡度。采用T-S模糊神經(jīng)網(wǎng)絡(luò)模型對(duì)干旱區(qū)土壤鹽分進(jìn)行預(yù)測(cè)評(píng)價(jià)，研究結(jié)果表明：

(1)利用野外采樣點(diǎn)所獲得的數(shù)據(jù)，客觀地揭示了地下水埋深、地下水礦化度、海拔高程以及坡度等土壤鹽漬化影響因子的空間分布特征。

(2)利用灰色關(guān)聯(lián)度分析得出，所選取的7個(gè)土壤鹽漬化影響因子可以用來構(gòu)建土壤鹽分預(yù)測(cè)模型。

(3)依據(jù)T-S模糊神經(jīng)網(wǎng)絡(luò)方法構(gòu)建土壤鹽分預(yù)測(cè)模型，將7個(gè)影響因子作為輸入變量，土壤鹽分作為輸出變量進(jìn)行模擬模型。通過模型預(yù)測(cè)結(jié)果得出，T-S模糊神經(jīng)網(wǎng)絡(luò)模型預(yù)測(cè)效果較好，可以用于預(yù)測(cè)土壤鹽漬化狀況。

(4)雖然T-S模糊神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)良好，但是，由于模型自身的原因加上數(shù)據(jù)獲取的限制，影響了土壤鹽分的預(yù)測(cè)模擬效果，需在今后的研究中進(jìn)行進(jìn)一步完善和檢驗(yàn)。