以地理信息系統(tǒng)為基礎(chǔ)的農(nóng)業(yè)地質(zhì)環(huán)境研究

胡品

摘要：根據(jù)地理信息系統(tǒng)（GIS）制圖方法對四川省內(nèi)江市雙橋鄉(xiāng)農(nóng)業(yè)地質(zhì)環(huán)境進(jìn)行分析。結(jié)果表明，雙橋鄉(xiāng)土壤Cr、Cu、Hg、Pb、As元素含量較低，大部分符合地質(zhì)環(huán)境一級標(biāo)準(zhǔn)；Zn、Ni、Cd元素含量較高，但也符合地質(zhì)環(huán)境質(zhì)量二級標(biāo)準(zhǔn)。通過疊加法將8種重金屬元素綜合考慮得到雙橋鄉(xiāng)絕大部分土壤區(qū)域?yàn)棰蝾愅寥溃急壤?7%左右。對土壤肥力分析可知，雙橋鄉(xiāng)土壤全鉀含量富足，有機(jī)質(zhì)碳與全氮含量水平一般，土壤全磷含量整體水平偏低，在種植水稻過程中需多施加含磷化肥。

關(guān)鍵詞：地理信息系統(tǒng)；農(nóng)業(yè)地質(zhì)環(huán)境；重金屬元素

中圖分類號：X141 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2017）07-1245-05

DOI：10.14088/j.cnki.issn0439-8114.2017.07.012

Agricultural Geological Environment Monitoring Based on

Geographic Information System

HU Pin

（Faculty of Geomatics， East China Institute of Technology， Nanchang 330000， China）

Abstract： The agricultural geological environment of Shuangqiao township， Neijiang city， Sichuan province was analyzed by the geographic information system（GIS） mapping method. The results showed that the content of Cr， Cu， Hg， Pb， As in study area was low， and the most soil belonged to the first standards of geological enviromental quality； the contents of Zn， Ni and Cd was relative high， but kept within the secondary standards. The results of comprehensive assessment on soil heavy metal by the superposition of 8 heavy metals showed that the most soil area of Shuangqiao township belonged to grade II soil， the proportion was about 97%. The analysis results of soil fertility indicated that soil total potassium content was rich， organic carbon and total nitrogen content was normal level， and total phosphorus content was low， so phosphorus fertilizer should be applied more during rice cultivation.

Key words： geographic information system； agricultural geological environment； heavy metal elements

隨著計算機(jī)軟件技術(shù)發(fā)展，地理信息系統(tǒng)（GIS）已經(jīng)越來越多地被運(yùn)用到農(nóng)業(yè)地質(zhì)數(shù)據(jù)處理之中。GIS在進(jìn)行環(huán)境治理與監(jiān)測、制定應(yīng)急預(yù)案以及防止自然災(zāi)害方面發(fā)揮了日益重要的作用。農(nóng)業(yè)地質(zhì)數(shù)據(jù)分析主要包括定量分析與定性分析。定性分析是運(yùn)用歸納演繹法對數(shù)據(jù)進(jìn)行抽象概括與總結(jié)，對獲取的材料進(jìn)行思維加工，進(jìn)而去粗取精、去偽存真、由表及里揭示數(shù)據(jù)內(nèi)在規(guī)律[1]。定量分析則是建立數(shù)學(xué)模型，通過法則計算分析對象的各個指標(biāo)，定量分析是將研究對象內(nèi)部關(guān)系以及外部環(huán)境以數(shù)字或公式表達(dá)出來，對研究對象得到一個精確的理性結(jié)果。在農(nóng)業(yè)地質(zhì)環(huán)境數(shù)據(jù)分析中，各重金屬元素間存在相互關(guān)聯(lián)，難以得到所需要的各種標(biāo)準(zhǔn)數(shù)據(jù)，導(dǎo)致結(jié)果誤差[2，3]。本研究將定量分析與定性分析結(jié)合起來，根據(jù)規(guī)范標(biāo)準(zhǔn)將地質(zhì)環(huán)境分為3類，采用地理信息系統(tǒng)（GIS）制圖方法對四川省內(nèi)江市雙橋鄉(xiāng)農(nóng)業(yè)地質(zhì)環(huán)境進(jìn)行分析，首先結(jié)合GIS原理給出了制圖基本方法，然后根據(jù)樣本數(shù)據(jù)繪制得到二維直觀評價圖，旨在為農(nóng)業(yè)地質(zhì)環(huán)境評價提供參考。

1 材料與方法

1.1 區(qū)域概況

四川省內(nèi)江市雙橋鄉(xiāng)轄區(qū)地質(zhì)結(jié)構(gòu)簡單，為第四系河流松散堆積層。地勢東北高西南低，海拔330～473 m，地貌以中、低丘陵為主，多為階梯狀與饅頭狀，谷寬100～200 m，平地少，坡地多。區(qū)域內(nèi)富含鈣、鎂鹽質(zhì)物質(zhì)，土壤pH中性，通透性好，適宜多種經(jīng)濟(jì)作物生長，但由于土壤坡度較大導(dǎo)致水土流失嚴(yán)重，肥力降低，土層變薄。由于發(fā)展水稻基礎(chǔ)條件較好，故目前多以種植冬水稻為主，年種植一季水稻。本研究農(nóng)業(yè)地質(zhì)環(huán)境以水稻為基礎(chǔ)。參照地質(zhì)環(huán)境質(zhì)量評價二級標(biāo)準(zhǔn)（GB 15618-95），將水田地質(zhì)環(huán)境質(zhì)量分為3個等級，如表1所示。

1.2 數(shù)據(jù)采集

采集雙橋鄉(xiāng)區(qū)域內(nèi)294組水稻根系土樣品，測定各樣品根系土元素含量。

1.3 基于GIS的制圖方法

以水稻根系元素含量為指導(dǎo)，結(jié)合數(shù)據(jù)的空間統(tǒng)計方法及繪圖矢量軟件展開農(nóng)業(yè)地質(zhì)環(huán)境分析?？臻g統(tǒng)計分析核心為認(rèn)識地理位置與數(shù)據(jù)之間的空間依賴性，或自相關(guān)程度，涉及空間權(quán)重構(gòu)建、空間自相關(guān)性與檢驗(yàn)相關(guān)性分析。空間自相關(guān)反映一個區(qū)域單元某一屬性與鄰近區(qū)域單元同一屬性值的相關(guān)程度，也反映某一屬性觀測數(shù)據(jù)間的潛在依賴性，是一種有效的分析技術(shù)[4]。自相關(guān)分為正相關(guān)與負(fù)相關(guān)，正相關(guān)指位置緊密的觀測單元傾向于相似，負(fù)相關(guān)指的是位置緊密的觀測單元傾向于相反，本研究運(yùn)用數(shù)據(jù)空間統(tǒng)計方法的正相關(guān)性。

由于全局統(tǒng)計不易發(fā)現(xiàn)位于不同區(qū)域的空間關(guān)聯(lián)模式，必須運(yùn)用局部自相關(guān)系數(shù)，通常采用局部Moran指數(shù)計算局部空間相關(guān)系數(shù)[5]。局部相關(guān)指數(shù)LISA是能夠滿足以下兩種要求的統(tǒng)計量：①每個區(qū)域單元LISA是描述周圍顯著相似值的空間集聚程度指標(biāo)；②所有觀測值LISA之和與全局空間聯(lián)系指標(biāo)成正比。一個LISA對應(yīng)的局部統(tǒng)計關(guān)系表達(dá)式為：

■Li=？姿？撰 （1）

式中，Li 為空間相關(guān)局部指標(biāo)，？撰為空間相關(guān)全局指標(biāo)，？姿為比例因子。

局部Moran指數(shù)Ii是局部LISA的一個特例，具備公式（2）中的含義。

Ii=■ （2）

式中，Xi為區(qū)域i的觀測值，Xj為區(qū)域j的觀測值，■為樣本均值，Wij為空間權(quán)重矩陣。檢驗(yàn)統(tǒng)計量可對有意義的局部空間關(guān)聯(lián)進(jìn)行顯著性檢驗(yàn)，局部Moran指數(shù)檢驗(yàn)公式為：

Z（Ii）=■ （3）

式中，Z（Ii）為檢驗(yàn)值，E（Ii）為樣本數(shù)據(jù)期望值，VAR（Ii）為樣本方差。

E（Ii）=W/（n-1） （4）

VAR（Ii）=■-E（Ii）2 （5）

其中，Wi=■Wij，Wi（2）=■Wij2，

b2=n■。

通過局部Moran指數(shù)檢驗(yàn)公式做出散點(diǎn)圖研究局部空間不穩(wěn)定性，散點(diǎn)圖用于數(shù)據(jù)的二維可視化，形成初步可觀察示圖，然后結(jié)合MapGIS軟件以及模糊綜合評價進(jìn)行下一步的圖像數(shù)據(jù)處理[6-8]。

2 結(jié)果與分析

2.1 土壤各重金屬元素含量

參照地質(zhì)環(huán)境質(zhì)量標(biāo)準(zhǔn)（GB 15618-1995）以及樣本數(shù)據(jù)對測試區(qū)土壤元素含量進(jìn)行分析，首先根據(jù)國標(biāo)（GB 15618-1995）將土壤分為Ⅰ類、Ⅱ類、Ⅲ類。Ⅰ類土壤重金屬元素含量值接近背景值未受后期污染，是最接近原始自然環(huán)境的土壤[9]；Ⅱ類土壤重金屬存在積累狀況，但不對植物與環(huán)境造成污染傷害；Ⅲ類土壤重金屬積累已超過臨界值，此類土壤種植的作物危害人體健康。根據(jù)表層土壤樣品數(shù)據(jù)，得到區(qū)域內(nèi)As、Hg元素含量所受污染情況最小，幾乎全部為Ⅰ類土，其余6種重金屬分別受到不同程度的污染，基于GIS重金屬評價獲取各元素分級面積及所占比例，具體見表2。

由表2可知，土壤As元素含量范圍為2.2～8.3 mg/kg，平均含量5.0 mg/kg，全區(qū)含量低，為Ⅰ類土壤區(qū)；Cr、Cu、Hg、Pb含量較低，Ⅰ類土壤面積所占比例均在94%以上；土壤Cd元素含量較高，均為Ⅱ類；土壤Zn和Ni元素Ⅱ類土壤所占比例較大。其中，Zn含量范圍在64.5～123.6 mg/kg，平均為91.0 mg/kg，Ⅰ類土壤區(qū)面積為783.9 hm2，Ⅱ類土壤區(qū)面積為438.9 hm2，分別占總面積的59.0%與41.0%，具體分布情況如圖1所示。土壤元素Ni含量較高，范圍在23.6～52.6 mg/kg，均值為38.2 mg/kg，觀測區(qū)域內(nèi)主要為Ⅱ類土壤，占總面積的73.8%，存在少量Ⅲ類土壤，面積2.4 hm2，占總面積的0.2%，Ⅲ類土壤區(qū)域集中分布于馬龍村附近，其余區(qū)域?yàn)棰耦愅寥溃伎偯娣e的26.0%，具體分布情況如圖2所示。從總體上來說，雙橋鄉(xiāng)As、Cr、Cu、Hg、Pb 5種金屬元素含量較低，大部分符合地質(zhì)環(huán)境質(zhì)量一級標(biāo)準(zhǔn)，Zn、Ni、Cd含量相對較高，但也符合地質(zhì)環(huán)境質(zhì)量二級標(biāo)準(zhǔn)。

2.2 土壤重金屬綜合評價

綜合評價是將As、Cr、Cu、Hg、Pb、Zn、Ni、Cd 8種元素含量進(jìn)行標(biāo)準(zhǔn)化處理后采取數(shù)據(jù)空間疊加方式編制重金屬綜合評價空間信息圖像[10]。疊合原則為量變引起質(zhì)變，8種元素分別單獨(dú)考慮時為Ⅰ類土壤，當(dāng)8種元素疊合考慮時，區(qū)域可能由Ⅰ類土壤變換為Ⅱ類土壤甚至Ⅲ類土壤?？紤]樣本數(shù)據(jù)Cd元素的異常分布特征，在進(jìn)行疊合考慮綜合評價時分為剔除與不剔除Cd元素兩種情況，通過MapGIS軟件進(jìn)行數(shù)據(jù)處理得到土壤重金屬評價情況，見表3、圖3、圖4。

綜合評價結(jié)果顯示，在剔除Cd元素情況下， Ⅰ類土壤區(qū)域所占比例僅為0.2%，絕大部分土壤區(qū)域?yàn)棰蝾愅寥溃急壤哌_(dá)98.1%，存在少量Ⅲ類土壤區(qū)域，所占比例為1.3%，Ⅲ類土壤區(qū)域主要分布于馬龍村附近（圖3）；在加入Cd元素情況下，Ⅰ類土壤區(qū)域完全消失，Ⅲ類土壤區(qū)域所占比例增大為3.5%，其余區(qū)域?yàn)棰蝾愅寥?，所占比例?6.5%（圖4）。

2.3 土壤肥力評價

土壤肥力是土壤生物性質(zhì)的綜合反映，土壤中有機(jī)質(zhì)碳、全氮、全磷、全鉀含量是評價土壤肥力的重要性指標(biāo)[11]，根據(jù)四川省項(xiàng)目“金土地工程”中農(nóng)業(yè)地質(zhì)調(diào)查成果確定分級評價標(biāo)準(zhǔn)，具體見表4。

通過GIS數(shù)據(jù)處理得到各元素分級面積及所占比例，具體見表5，分布情況見圖5～圖8。雙橋鄉(xiāng)土壤有機(jī)質(zhì)碳含量主要在四級水平，占55.8%，五級水平占38.1%，有小部分二級、三級和六級區(qū)域，說明有機(jī)質(zhì)碳含量在整體區(qū)域內(nèi)跨度較大，但分布較為均勻，主要在四級、五級水平；土壤全氮含量主要分布在三級、四級水平，分別占51.3%與42.3%，由此可見土壤全氮含量中等肥力區(qū)占到了93.6%。土壤全磷含量主要分布在四級水平，所占比例為76.4%，五級區(qū)域所占比例為23.3%，說明土壤全磷含量整體水平偏低，在種植水稻過程中可適當(dāng)施加含磷化肥，提高土壤肥力。土壤全鉀含量主要分布在二級水平，所占比例高達(dá)85.3%，說明雙橋鄉(xiāng)土壤鉀含量均衡，且為高水平，利于農(nóng)作物生長。

3 結(jié)論

以四川省內(nèi)江市雙橋鄉(xiāng)為研究對象，對采集的樣本數(shù)據(jù)結(jié)合地理信息系統(tǒng)進(jìn)行分析，得出以下幾點(diǎn)結(jié)論。

1）根據(jù)GIS原理以及樣本數(shù)據(jù)給出了數(shù)據(jù)處理方法，并根據(jù)樣本數(shù)據(jù)得到了各元素二維直觀評價圖，為觀察數(shù)據(jù)規(guī)律做出了貢獻(xiàn)。

2）土壤單項(xiàng)重金屬含量評價結(jié)果表明，Cr、Cu、Hg、Pb、As含量較低，在整個雙橋鄉(xiāng)大部分為Ⅰ類土壤區(qū)；Zn、Ni、Cd元素含量相對較高，但也符合地質(zhì)環(huán)境質(zhì)量二級標(biāo)準(zhǔn)。

3）將8種重金屬元素綜合考慮，進(jìn)行土壤重金屬綜合評價得到雙橋鄉(xiāng)絕大部分土壤區(qū)域?yàn)棰蝾愅寥?，在剔除Cd元素情況下，Ⅱ類土壤所占比例高達(dá)98.1%；在不剔除Cd元素情況下，Ⅱ類土壤所占比例為96.5%。

4）根據(jù)土壤肥力評價得出，雙橋鄉(xiāng)土壤有機(jī)質(zhì)碳與全氮含量水平一般，土壤全鉀含量富足，土壤全磷含量整體水平偏低，全磷含量主要分布在四、五級水平，所占比例高達(dá)99.7%，在種植水稻過程中需多施加含磷化肥。

參考文獻(xiàn)：

[1] 劉應(yīng)平.成都平原土壤地球化學(xué)組成及其生態(tài)響應(yīng)研究[D].成都：成都理工大學(xué)，2012.

[2] 葛 文.山東煙臺地區(qū)土壤地球化學(xué)環(huán)境與優(yōu)質(zhì)蘋果生產(chǎn)的適應(yīng)性評價[D].武漢：中國地質(zhì)大學(xué)，2013.

[3] 范文靜.甘肅黃河石林地質(zhì)景觀可持續(xù)發(fā)展研究[D].北京：中國地質(zhì)大學(xué)（北京），2014.

[4] 魏莉莉.張掖市城市濕地環(huán)境條件評價及保護(hù)研究[D].西安：長安大學(xué)，2014.

[5] 許文超.五大連池地質(zhì)公園生態(tài)農(nóng)業(yè)觀光園規(guī)劃設(shè)計研究[D].哈爾濱：東北農(nóng)業(yè)大學(xué)，2015.

[6] 田言亮，張光輝，王 茜，等.黃淮海平原灌溉農(nóng)業(yè)對地下水依賴程度與保障能力[J].地球?qū)W報，2016，37（3）：257-265.

[7] 張光輝，費(fèi)宇紅，王 茜，等.灌溉農(nóng)業(yè)的地下水保障能力評價方法研究——黃淮海平原為例[J].水利學(xué)報，2016，47（5）：608-615.

[8] 武 超，張兆吉，費(fèi)宇紅，等.天津污灌區(qū)水稻土壤汞形態(tài)特征及其食品安全評估[J].農(nóng)業(yè)工程學(xué)報，2016，32（18）：207-212.

[9] 張光輝，田言亮，王電龍，等.冀中山前農(nóng)業(yè)區(qū)地下水位強(qiáng)降弱升特征與機(jī)制[J].水科學(xué)進(jìn)展，2015，26（2）：227-232.

[10] 童 欣，甘義群，夏 友.武漢市地質(zhì)環(huán)境綜合評價與區(qū)劃[J].水文地質(zhì)工程地質(zhì)，2015，42（3）：149-155.

[11] 盧耀如，張鳳娥，劉 琦，等.建設(shè)生態(tài)文明保障新型城鎮(zhèn)群環(huán)境安全與可持續(xù)發(fā)展[J].地球?qū)W報，2015，36（4）：403-412.