基于BP 神經(jīng)網(wǎng)絡(luò)對(duì)某電鍍廠土壤重金屬預(yù)測(cè)及人體健康風(fēng)險(xiǎn)評(píng)價(jià)

黃煜韜, 施維林,* , 紀(jì)娟, 孫守鈞, 劉陽(yáng), 吳錦平

1. 蘇州科技大學(xué)環(huán)境科學(xué)與工程學(xué)院,蘇州 215000

2. 蘇州市太倉(cāng)生態(tài)環(huán)境局土壤環(huán)境科,蘇州 215400

3. 蘇州科技城科技創(chuàng)新服務(wù)中心,蘇州 215400

4. 蘇州逸凡特環(huán)境修復(fù)有限公司,蘇州 215010

5. 昆山市巴城鎮(zhèn)環(huán)境保護(hù)辦,蘇州 215311

隨著國(guó)民經(jīng)濟(jì)的不斷提升,高能耗、高污染的電鍍企業(yè)逐漸關(guān)停,遺留了大量污染地塊,這些地塊可能會(huì)造成環(huán)境污染[1-2]。 關(guān)?；虬徇w的電鍍企業(yè)由于其污染具有長(zhǎng)期性、累積性和潛伏性,將會(huì)對(duì)人體健康產(chǎn)生巨大威脅[3]。 因此,迫切需要對(duì)這類電鍍遺留地塊進(jìn)行健康風(fēng)險(xiǎn)評(píng)價(jià)。

目前,針對(duì)電鍍企業(yè)環(huán)境健康風(fēng)險(xiǎn)的相關(guān)研究主要集中在電鍍排污對(duì)周邊環(huán)境的影響評(píng)價(jià)上,而對(duì)遺留的電鍍企業(yè)場(chǎng)地本身的土壤環(huán)境健康風(fēng)險(xiǎn)研究相對(duì)比較薄弱[4-8]。 美國(guó)環(huán)境保護(hù)局(United States Environmental Protection Agency, US EPA)最早提出“brown plot”的健康風(fēng)險(xiǎn)評(píng)估理念,其RBCA 模型已經(jīng)運(yùn)行多年。 并且相繼制定了技術(shù)導(dǎo)則、風(fēng)險(xiǎn)評(píng)估指南和相關(guān)法律法規(guī),并進(jìn)行典型污染場(chǎng)地風(fēng)險(xiǎn)評(píng)價(jià)和后續(xù)治理工作[9-11]。 我國(guó)的土壤保護(hù)、風(fēng)險(xiǎn)評(píng)估和修復(fù)治理仍處于發(fā)展初期。 目前國(guó)內(nèi)研究多側(cè)重于利用簡(jiǎn)單公式來(lái)計(jì)算污染程度,已有學(xué)者提出用簡(jiǎn)單公式計(jì)算健康風(fēng)險(xiǎn)的不足之處,Gu 等[12]提出簡(jiǎn)單公式不足以評(píng)估污染土壤對(duì)人類的潛在影響,因?yàn)橐环N重金屬在環(huán)境中的毒性和流動(dòng)性強(qiáng)烈依賴于重金屬的空間分布。 Mashal 等[13]認(rèn)為公式法無(wú)法全面有效說(shuō)明重金屬對(duì)人體的危害,需要綜合考慮其在空間的流動(dòng)性。 簡(jiǎn)單公式缺乏從區(qū)域縱向空間角度進(jìn)行的場(chǎng)地重金屬污染程度和人體健康風(fēng)險(xiǎn)評(píng)價(jià)[14-15]。 然而通過(guò)污染預(yù)測(cè)的方法能夠直觀反映縱向空間的污染程度[16]。 因此,本文將江蘇省某電鍍企業(yè)作為研究對(duì)象,通過(guò)人工神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)土壤重金屬含量與人體健康風(fēng)險(xiǎn)評(píng)價(jià)相結(jié)合的方法,分析重金屬在該地塊土壤中污染狀況和健康風(fēng)險(xiǎn)的空間特征,為地塊污染修復(fù)提供理論支持。

1 材料與方法(Materials and methods)

1.1 研究區(qū)概況

某電鍍企業(yè)建于1968年,2013年停產(chǎn),占地面積約15 349.6 m2,主要從事自行車配件、五金拉手等生產(chǎn)。 目前廠區(qū)內(nèi)仍留存著各種電鍍生產(chǎn)工具和原輔材料存儲(chǔ)罐,該廠南端污水處理池已停止運(yùn)行,而廢水處理設(shè)備仍在原地。

1.2 土壤樣品的采集與分析

明確研究區(qū)域內(nèi)的功能區(qū)劃分后,通過(guò)專業(yè)判斷布點(diǎn)法在研究區(qū)域內(nèi)布設(shè)33個(gè)采樣點(diǎn),采樣深度為 1.5、3、4.5 和 6 m,采樣點(diǎn)的布設(shè)如圖 1 所示[17]。各采樣點(diǎn)利用GPS 工具進(jìn)行定位,利用Geoprobe 7822 DT 無(wú)擾動(dòng)直推采樣鉆機(jī)分層采樣,去除動(dòng)植物殘?bào)w、石子等異物后將各層土樣混合均勻,根據(jù)四分法取樣裝入樣品袋編號(hào)備用[18]。 在自然條件下風(fēng)干、磨碎、過(guò)篩后對(duì)樣品進(jìn)行檢測(cè)。

檢測(cè)重金屬時(shí)所用試劑均為優(yōu)級(jí)純,使用空白和國(guó)家標(biāo)準(zhǔn)物質(zhì)(GBW07405)進(jìn)行分析質(zhì)量控制,重金屬As 和Hg 依據(jù)《土壤質(zhì)量 總汞、總砷、總鉛的測(cè)定 原子熒光法》(GB/T 22105.2—2008)[19]檢測(cè),檢測(cè)限為 0.008 mg·kg-1和 0.001 mg·kg-1,重金屬 Cd、Cr、Cu、Ni、Pb 和 Zn 采用普通酸分解法-電感耦合等離子體發(fā)射光譜法檢測(cè),前處理參考《土壤環(huán)境監(jiān)測(cè)技術(shù)規(guī)范》(HJ/T 166—2004)[20],分析參考《電感耦合等離子體原子發(fā)射光譜法》(US EPA 6010C-2007),檢測(cè)限分別為 0.042、0.097、0.034、0.192、0.321和 0.008 mg·kg-1。

1.3 重金屬污染評(píng)價(jià)

地質(zhì)累積指數(shù)法是德國(guó)科學(xué)家Muller 在1979年提出的利用某一種重金屬含量與地球化學(xué)背景值的比較來(lái)確定土壤重金屬污染水平的定量指標(biāo)[21-22],表達(dá)式如下:

式中:Cn為重金屬的測(cè)定濃度(mg·kg-1);Bn為重金屬相應(yīng)的背景參考值(mg·kg-1);1.5 為中和背景值變化的常數(shù)。Igeo與污染水平的對(duì)應(yīng)關(guān)系如下:無(wú)污染(Igeo≤0)、輕度污染(0＜Igeo≤1)、中度污染(1＜Igeo≤2)、偏重度污染(2＜Igeo≤3)、重度污染(3＜Igeo≤4)、嚴(yán)重污染(4＜Igeo≤5)、極嚴(yán)重污染(Igeo＞5)[23]。

1.4 人體健康風(fēng)險(xiǎn)評(píng)價(jià)

根據(jù)US EPA 推薦的土壤重金屬暴露風(fēng)險(xiǎn)評(píng)估方法,將暴露途徑分為攝入、吸入、皮膚接觸。 將成人作為研究對(duì)象用于評(píng)估和分析。 潛在有毒金屬的日平均劑量通過(guò)對(duì)成人的3 條路徑進(jìn)行計(jì)算,公式如下:

式中:ADDing、ADDinh和 ADDdermal分別為每日平均通過(guò)土壤攝入、吸入和皮膚接觸的暴露量(mg·kg-1·d-1),公式中土壤暴露參數(shù)及毒性數(shù)據(jù)如表1 所示。

表1 土壤暴露參數(shù)及毒性數(shù)據(jù)Table 1 Soil exposure parameters and toxicity data

非致癌風(fēng)險(xiǎn)危險(xiǎn)商(HQ)用于評(píng)估單一重金屬對(duì)人體的非致癌風(fēng)險(xiǎn)[24]。 HQ 的總和即為HI,可用于評(píng)估大量重金屬的總體潛在非致癌風(fēng)險(xiǎn)。 致癌風(fēng)險(xiǎn)(CR)是個(gè)體由于暴露于致癌危害而發(fā)展為癌癥的概率[25-26]。 HQ、HI 和 CR 使用以下公式確定:

式中:RfDi是參考劑量,SFi是癌癥斜率因子(kg·d·mg-1)。 以下水平表示非致癌風(fēng)險(xiǎn):HI＜1 或 HQ＜1表示無(wú)不良健康影響:HI＞1 或HQ＞1 表示有不良健康影響。 以下水平表示致癌風(fēng)險(xiǎn):CR＜10-6,無(wú)致癌風(fēng)險(xiǎn);10-6＜CR＜10-4,可接受水平的致癌風(fēng)險(xiǎn);當(dāng)CR＞10-4,致癌風(fēng)險(xiǎn)較高[27]。

1.5 BP 神經(jīng)網(wǎng)絡(luò)

人工神經(jīng)網(wǎng)絡(luò)是由大量簡(jiǎn)單的處理單元(神經(jīng)元)組成的高度網(wǎng)絡(luò)化的人工網(wǎng)絡(luò)[28]。 人工神經(jīng)網(wǎng)絡(luò)可以模擬人體大腦神經(jīng)系統(tǒng)的功能和構(gòu)造,確定輸入輸出之間非線性映射關(guān)系的全局性分析方法[29]。 人工神經(jīng)網(wǎng)絡(luò)在環(huán)境領(lǐng)域主要應(yīng)用于環(huán)境質(zhì)量變化以及綜合污染評(píng)價(jià)[30]。

BP 神經(jīng)網(wǎng)絡(luò)構(gòu)建中,分別選取采樣點(diǎn)位的經(jīng)度、緯度、深度和pH 值作為模型的輸入層。 重金屬預(yù)測(cè)值作為模型的輸出值。 隱含層的確定較為復(fù)雜,通過(guò)定義輸入向量為X=(x1,x2,x3,...,xm)T,隱藏層輸出向量為S=(s1,s2,s3,...,sn)T,期望輸出向量Y=(y1,y2,y3,...,yj)T,則第n個(gè)神經(jīng)元輸出推算為:

式中:vn為每個(gè)節(jié)點(diǎn)的輸出值;f1為激活函數(shù);ω 為各節(jié)點(diǎn)m與n之間的權(quán)值;An為節(jié)點(diǎn)的閾值。

輸出層第j個(gè)神經(jīng)元輸出為:

式中:yj為每個(gè)節(jié)點(diǎn)的輸出值;f2為激活函數(shù);ω 為各節(jié)點(diǎn)n與j之間的權(quán)值;Bn為節(jié)點(diǎn)的閾值。

在BP 神經(jīng)網(wǎng)絡(luò)算法中,輸入層的結(jié)果作為輸入變量輸出到隱藏層。 在隱藏層中,數(shù)據(jù)通過(guò)算法進(jìn)行處理,隱藏層的輸出進(jìn)入輸出層。 當(dāng)輸出層的結(jié)果與期望值不同時(shí),誤差會(huì)在神經(jīng)網(wǎng)絡(luò)中反向傳播,權(quán)值會(huì)通過(guò)原路徑逐層修改,最終達(dá)到期望值[31]。 通過(guò)上述公式采用試湊法和經(jīng)驗(yàn)公式法相結(jié)合的方法,對(duì)不同數(shù)量的隱藏層和神經(jīng)元個(gè)數(shù)的模擬精度進(jìn)行比較,最終確定隱藏層層數(shù)及神經(jīng)元個(gè)數(shù)。

2 結(jié)果與討論(Results and discussion)

2.1 土壤重金屬統(tǒng)計(jì)分析

研究區(qū)土壤重金屬元素含量如表2 所示。 土壤中 Hg、Cd、Cr、Cu、Ni、Pb 和 Zn 這 7 種重金屬均值都超過(guò)江蘇省土壤背景值,超標(biāo)率分別為38.6%、63.6%、72.7%、81.1%、77.3%、72.0% 和 71.2%。 As元素均值低于江蘇省土壤背景值,但仍有部分點(diǎn)位超標(biāo),超標(biāo)率為22.7%。 研究區(qū)土壤中 Cr 含量遠(yuǎn)遠(yuǎn)超過(guò)《土壤環(huán)境質(zhì)量 建設(shè)用地土壤污染風(fēng)險(xiǎn)管控標(biāo)準(zhǔn)》(GB 36600—2018)[32]中“第二類用地”Cr的篩選值,超標(biāo)率為100%;Ni 和Pb 分別在一車間和廢水處理區(qū)超標(biāo)。 表明該研究區(qū)可能已受到了這8 種重金屬的污染。

表2 研究區(qū)域土壤重金屬含量Table 2 Heavy metal content in soil of study area

變異系數(shù)(CV)可反映重金屬的離散性以及人類活動(dòng)對(duì)重金屬的影響,值越大,表明人類活動(dòng)干擾就越強(qiáng)[33]。 8 種重金屬變異系數(shù)依次為 Pb＞Ni＞Cr＞Zn＞Cu＞Hg＞Cd＞As＞36%,根據(jù)變異系數(shù)的劃分,CV＜15%屬于輕度變異,15%＜CV＜36%屬于中等變異,CV＞36%屬于高度變異[34]。 因此8 種重金屬均屬于高度變異,說(shuō)明8 種重金屬樣本含量的離散性較大,且受到人類影響較大。

2.2 重金屬污染評(píng)價(jià)

對(duì)不同深度的重金屬地質(zhì)累積指數(shù)進(jìn)行可視化繪圖,由圖 2 可知,As 在 1.5 m 處Igeo＜0,屬于無(wú)污染,而在6.0 m 處1＜Igeo≤2,屬于中度污染,可能的原因是污染物已經(jīng)下滲到下層土壤而長(zhǎng)時(shí)間未進(jìn)行生產(chǎn)作業(yè)對(duì)于表層土壤的污染有所減輕。 Cr 在1.5 m 處屬于偏重度污染,污染較為嚴(yán)重,在6.0 m 處污染主要集中在11 號(hào)、14 號(hào)和26 號(hào)點(diǎn)位。 Ni 在1.5 m 處 2＜Igeo≤3,屬于偏重度污染,而在 6.0 m 處僅在26 號(hào)點(diǎn)位屬于偏重度污染。 Hg、Cd、Cu、Pb 和 Zn均處于無(wú)污染與輕度污染之間,其中污染較為嚴(yán)重的點(diǎn)位仍是11 號(hào)、14 號(hào)和26 號(hào)點(diǎn)位,其功能區(qū)為廢水處理區(qū)、一車間。 超標(biāo)的污染因子與電鍍企業(yè)的電鍍品種一致,主要是電鍍過(guò)程中對(duì)含重金屬?gòu)U水和灰渣的處理不當(dāng)導(dǎo)致滲漏滲入土壤造成土壤污染。

圖2 土壤重金屬地質(zhì)累積指數(shù)法可視圖(Igeo)注:(a) As,(b) Hg,(c) Cd,(d) Cr,(e) Cu,(f) Ni,(g) Pb,(h) Zn;上部為1.5 m 處重金屬地質(zhì)累積指數(shù),下部為6.0 m 處重金屬地質(zhì)累積指數(shù)。Fig.2 A visibility graph of the soil heavy metal geological accumulation index method (Igeo)Note: (a) As, (b) Hg, (c) Cd, (d) Cr, (e) Cu, (f) Ni, (g) Pb, (h) Zn; the upper part shows the geoaccumulation index of heavy metals at 1.5 m, and the lower part shows the geoaccumulation index of heavy metals at 6.0 m.

2.3 人體健康風(fēng)險(xiǎn)評(píng)價(jià)

利用污染場(chǎng)地健康風(fēng)險(xiǎn)模型,計(jì)算研究區(qū)內(nèi)致癌和非致癌風(fēng)險(xiǎn)。 利用Origin 繪制重金屬在不同深度下3 條暴露途徑的致癌風(fēng)險(xiǎn)和非致癌風(fēng)險(xiǎn)。 選取了US EPA 綜合風(fēng)險(xiǎn)信息系統(tǒng)(Integrated Risk and Information System, IRIS)所提供的非致癌物與致癌物作為健康風(fēng)險(xiǎn)的評(píng)估依據(jù),重金屬 As、Cd、Cr 和Ni 屬于致癌物,致癌風(fēng)險(xiǎn)如圖3 所示。 重金屬Cr在各深度致癌風(fēng)險(xiǎn)均＞10-4,在1.5 m 處致癌風(fēng)險(xiǎn)最高達(dá)到1.7×10-3,隨著采樣點(diǎn)深度的增加,Cr 的致癌風(fēng)險(xiǎn)略有降低,但仍超過(guò)可接受的致癌風(fēng)險(xiǎn)10-6。重金屬As 和Ni 在各深度致癌風(fēng)險(xiǎn)亦均＞10-6。

圖3 不同深度下各重金屬致癌風(fēng)險(xiǎn)圖注:(a)1.5 m;(b)3.0 m;(c)4.5 m;(d)6.0 m。Fig.3 Carcinogenic risk maps of various heavy metals at different depthsNote: (a)1.5 m; (b)3.0 m; (c)4.5 m; (d)6.0 m.

重金屬 As、Hg、Cd、Cr、Cu、Ni、Pb 和 Zn 的非致癌風(fēng)險(xiǎn)如圖4 所示。 在各個(gè)深度下,僅重金屬Cr 的非致癌風(fēng)險(xiǎn)危險(xiǎn)商HQ＞1,說(shuō)明研究區(qū)內(nèi)Cr 已經(jīng)超過(guò)可接受的非致癌風(fēng)險(xiǎn)。 各重金屬HQ 依次為Cr＞Pb＞Ni＞As＞Cu＞Zn＞Hg＞Cd,基本與地質(zhì)累積指數(shù)的結(jié)論相符。 主要是因?yàn)殡婂冃袠I(yè)在退鍍過(guò)程中會(huì)產(chǎn)生大量 Cd、As、Cu 和 Cr 粉塵,拋光時(shí)也會(huì)產(chǎn)生重金屬粉塵導(dǎo)致污染,結(jié)合研究區(qū)內(nèi)工藝,在銅件酸洗、鋁件鈍化、鋁件電化學(xué)拋光等流程中,經(jīng)常使用鉻酸鹽,由于處理不當(dāng),而導(dǎo)致廢水中含有大量Cr,滲透進(jìn)入土壤造成污染。

圖4 不同深度下各重金屬非致癌風(fēng)險(xiǎn)圖注:(a)1.5 m;(b)3.0 m;(c)4.5 m;(d)6.0 m。Fig.4 Non-carcinogenic risk maps of heavy metals at different depthsNote: (a)1.5 m; (b)3.0 m; (c)4.5 m; (d)6.0 m.

2.4 BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)

在本研究區(qū)中,根據(jù)數(shù)據(jù)的可得性和合理性,設(shè)置了4個(gè)輸入層節(jié)點(diǎn),分別為采樣點(diǎn)位的經(jīng)度、緯度、深度和pH 值,將重金屬濃度預(yù)測(cè)值作為輸出層進(jìn)而計(jì)算深層致癌風(fēng)險(xiǎn)和非致癌風(fēng)險(xiǎn)。 最終確定為4-20-16-8-6-1 的BP 神經(jīng)網(wǎng)絡(luò),各層之間以 Purelin函數(shù)連接,然后進(jìn)行線性擬合,從而得到重金屬實(shí)際值與預(yù)測(cè)值之間的相關(guān)性,其拓?fù)淠Ｐ徒Y(jié)構(gòu)如圖5所示。 選取超過(guò)可接受致癌風(fēng)險(xiǎn)的 Ni、Cr、As 以及超過(guò)可接受非致癌風(fēng)險(xiǎn)的Cr 作為預(yù)測(cè)對(duì)象,預(yù)測(cè)模型擬合回歸評(píng)價(jià)如圖6 所示。 為驗(yàn)證回歸模型的精度與準(zhǔn)確度,引入評(píng)價(jià)指標(biāo)決定系數(shù)(R2)、均方根誤差(RMSE)、平均絕對(duì)誤差(MAE)和平均相對(duì)誤差(MRE)。 精度評(píng)價(jià)指標(biāo)如表3 所示,R2是衡量回歸模型擬合程度的指標(biāo),從預(yù)測(cè)模型擬合回歸中可以看出,R2均接近于1,說(shuō)明回歸方程的擬合效果較好。 RMSE、MAE 和 MRE 較小,MRE 不超過(guò) 30%即可說(shuō)明預(yù)測(cè)模型精度較高,該模型對(duì)深層土壤中重金屬Ni、Cr 和As 的預(yù)測(cè)有很好的適用性。

表3 BP 神經(jīng)網(wǎng)絡(luò)訓(xùn)練樣本和檢驗(yàn)樣本精度評(píng)價(jià)指標(biāo)Table 3 BP neural network training sample and test sample accuracy evaluation index

圖5 BP 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)模型Fig.5 BP neural network structure model

圖6 BP 神經(jīng)網(wǎng)絡(luò)污染預(yù)測(cè)線性圖注:(a) Ni;(b) Cr;(c) As。Fig.6 BP neural network pollution prediction linear graphNote: (a) Ni; (b) Cr; (c) As.

通過(guò)上述神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)預(yù)測(cè)超過(guò)可接受水平重金屬 Ni、Cr 和 As 在深層土壤 7.5 m 和 9.0 m 處的人類健康風(fēng)險(xiǎn),以判斷重金屬污染對(duì)深層土壤的影響作用和污染深度。 如圖7 所示,深層土壤中重金屬的致癌風(fēng)險(xiǎn)和非致癌風(fēng)險(xiǎn)明顯下降。 從深層土壤(7.5 m 和9.0 m)來(lái)看,重金屬As 的致癌風(fēng)險(xiǎn)已小于0 ~6.0 m 土壤致癌風(fēng)險(xiǎn)的1/2,但仍＞10-6,超過(guò)可接受水平。 此外重金屬Cr 在深層土壤7.5 m 和9.0 m處的致癌風(fēng)險(xiǎn)亦＞10-6,但與0 ~6.0 m 土壤相比,風(fēng)險(xiǎn)小了一個(gè)數(shù)量級(jí)。 經(jīng)BP 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè),重金屬Ni 在深層土壤(7.5 m 和9.0 m)的致癌風(fēng)險(xiǎn)＜10-6,表明重金屬Ni 對(duì)深層土壤已無(wú)污染。 在深層土壤(7.5 m 和9.0 m),重金屬Cr 的非致癌風(fēng)險(xiǎn)亦＜1,屬于可接受水平,表明Cr 在深層土壤已無(wú)非致癌風(fēng)險(xiǎn)。

圖7 深層土壤污染預(yù)測(cè)圖注:(a) As 致癌風(fēng)險(xiǎn);(b) Cr 致癌風(fēng)險(xiǎn);(c) Ni 致癌風(fēng)險(xiǎn);(d) Cr 非致癌風(fēng)險(xiǎn)。Fig.7 Deep soil pollution slice mapNote: (a) Carcinogenic risk of As; (b) Carcinogenic risk of Cr; (c) Carcinogenic risk of Ni; (d) Non-carcinogenic risks of Cr.

本研究以蘇州市土壤重金屬含量為研究對(duì)象,將BP 神經(jīng)網(wǎng)絡(luò)和人體健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)合,能夠更好地?cái)M合重金屬的空間變異規(guī)律,具有較優(yōu)的RMSE、MAE 和 MRE 值。 通過(guò) BP 神經(jīng)網(wǎng)絡(luò)和人體健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)合的方式可以為場(chǎng)地監(jiān)測(cè)與分析提供科學(xué)依據(jù),管理決策部門能夠通過(guò)模型預(yù)測(cè)結(jié)果對(duì)場(chǎng)地重金屬的空間分布及其污染狀況有全面準(zhǔn)確的了解,對(duì)場(chǎng)地污染狀況做出科學(xué)的判斷,并為污染場(chǎng)地修復(fù)提供理論指導(dǎo)。 本研究?jī)H采用了經(jīng)緯度、

pH 值、深度作為輸入層經(jīng)模型擬合重金屬空間變異規(guī)律,在今后的研究中還可以將其他土壤屬性信息如土壤有機(jī)質(zhì)含量、土壤滲透率以及土壤巖性帶的變化等作為輔助信息輸入,進(jìn)一步提高模型的準(zhǔn)確性。