廣花盆地地下水三氮時(shí)空分布特征及影響因素分析

龐園，李志威，張明珠

廣州市水務(wù)科學(xué)研究所，廣東 廣州 510220

廣花盆地地下水資源豐富，年均天然地下水資源量為 2.6×108m3，年均地下水允許開(kāi)采量為1.97×108m3。2011年，廣花盆地被廣東省劃定為廣州市的地下水應(yīng)急水源區(qū)。作為廣州市唯一的地下水應(yīng)急水源區(qū)，廣花盆地的地下水質(zhì)量直接關(guān)系著廣州市的供水安全。據(jù)調(diào)查，目前廣州市地下水的主要污染類型為氮污染。地下水環(huán)境中，氮主要以離子態(tài)的硝態(tài)氮、亞硝態(tài)氮及氨氮等3種形式存在，簡(jiǎn)稱三氮（於嘉聞等，2016）。地下水中的氮主要來(lái)源于人類活動(dòng)，包括化學(xué)肥料、農(nóng)家肥、生活污水和生活垃圾以及工業(yè)污染等（許可等，2011）。農(nóng)藥化肥的大量使用、生活污水的隨意排放以及工業(yè)污染等，都為氮元素在地下水中的累積提供了條件（於嘉聞等，2016），這也是國(guó)內(nèi)外地下水三氮普遍超標(biāo)的主要原因。關(guān)于地下水三氮的研究較多，主要有地下水硝酸鹽的轉(zhuǎn)化過(guò)程與機(jī)理，地下水硝酸鹽污染的來(lái)源以及避免污染的方法和措施，地下水硝酸鹽氮、亞硝酸鹽氮的確定方法，如何減少或除去地下水中的硝酸鹽，以及一些地下水硝酸鹽污染與人體健康的研究等（陳建耀等，2006）。近幾年，國(guó)外的研究主要集中在三氮的處理技術(shù)研究（Glenn et al.，2017；Vance et al.，2015；Huang et al.，2014；Nissim et al.，2014）和三氮的來(lái)源及轉(zhuǎn)化機(jī)理研究（Couturier et al.，2017；Nilsson et al.，2017；Ki et al.，2015；Moran et al.，2014）兩個(gè)方面。中國(guó)馮海波等（2017）、於嘉聞等（2016）、張婷等（2014）、萬(wàn)長(zhǎng)園等（2014）各自對(duì)地下水的三氮時(shí)空分布特征以及影響因素進(jìn)行了研究，結(jié)果表明地下水中三氮的含量受多種因素影響，包括氣象條件、農(nóng)業(yè)活動(dòng)、工業(yè)污水、氧化還原條件、地表水水源、包氣帶巖性、地下水徑流條件、潛水埋深、土地利用類型和生活污染等。本文在前人研究的基礎(chǔ)上，基于廣花盆地80個(gè)地下水監(jiān)測(cè)井2012—2016年雨季和旱季共 10次采樣數(shù)據(jù)，利用ArcGIS的數(shù)據(jù)管理、空間分析和制圖工具分雨季和旱季分析了研究區(qū)松散巖類孔隙水、基巖裂隙水和巖溶水中含氮污染物的濃度分布和時(shí)空分布特征，并探討了降雨、地下水埋藏深度和人類活動(dòng)與含氮污染物之間的關(guān)系，包括利用相關(guān)系數(shù)定量分析了降雨量與含氮污染物之間的相關(guān)性，采用灰色關(guān)聯(lián)法計(jì)算了耕地面積、化肥施用量、降雨量與含氮污染物之間的關(guān)聯(lián)度，以及從垂向和水平兩個(gè)方向探討了含氮污染物含量變化的原因。研究結(jié)果可為廣州市制定地下水保護(hù)與利用規(guī)劃提供理論基礎(chǔ)和科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)域概況

廣花盆地位于珠江三角洲東北部邊緣，地理坐標(biāo)為北緯 23°00′～23°30′，東經(jīng) 112°57′～113°26′，面積約858 km2（龐園等，2017）。廣花盆地地處南亞熱帶，雨量充沛，多年平均降雨量1818.7 mm，多年平均蒸發(fā)量1611.7 mm，大氣降雨是地下水補(bǔ)給的主要來(lái)源。區(qū)內(nèi)江、河、湖、庫(kù)、渠等地表水體較發(fā)達(dá)，主要河流有流溪河、巴江河、新街水、天馬河等（龐園等，2017），地表水入滲也是地下水的重要來(lái)源。

廣花盆地地勢(shì)大致為北部高、南部低，地貌類型有丘陵、臺(tái)地和平原等，以平原為主，臺(tái)地和丘陵次之（龐園等，2017）。基巖山區(qū)及丘陵臺(tái)地的地下水主要接受大氣降雨的入滲補(bǔ)給，山塘、水庫(kù)及渠道等對(duì)其入滲補(bǔ)給的規(guī)模較小。平原區(qū)地下水的補(bǔ)給來(lái)源較為豐富，除主要接受大氣降雨的入滲補(bǔ)給外，還接受地表河流及灌溉的入滲、渠道滲流、基巖山區(qū)及丘陵臺(tái)地地下水的側(cè)向補(bǔ)給等。

廣花盆地地下水主要有松散巖類孔隙水、基巖裂隙水和巖溶水三大類，松散巖類孔隙水廣泛分布于廣花盆地平原、山間谷地及山前地帶；基巖裂隙水廣泛分布于低山、丘陵及臺(tái)地；巖溶水呈“北東—南西向”條帶狀展布于廣花盆地中，以覆蓋型巖溶水分布為主（龐園等，2017）。總體上盆地地下水從北向南流。松散巖類孔隙水受到壟狀丘陵的限制，其總體上自北部、東部及西北部匯入廣花盆地后轉(zhuǎn)向南徑流，流向大致與巖溶水相同；古河道地下水流向受古河道展布制約，徑流以沿古河道滲流為主。外圍丘陵山區(qū)一般地勢(shì)較高，地形起伏較大，地下水水力坡度大，加上溝谷切割較深和巖石節(jié)理裂隙發(fā)育，故淺部基巖裂隙水獲得補(bǔ)給后經(jīng)過(guò)短暫的徑流，便以泉或滲流形式排入附近的溪流；深層地下水則通過(guò)斷層和裂隙向廣花盆地匯流。覆蓋型巖溶水分布于背、向斜構(gòu)造軸部，兩側(cè)受砂頁(yè)巖限制，故其沿灰?guī)r條帶自北東、北北東向運(yùn)移；東部良田—肖崗一帶山前平原地下水自東向西匯入廣花盆地后轉(zhuǎn)向南徑流。

廣花盆地地下水水位動(dòng)態(tài)變化主要受大氣降雨影響，地下水水位變化受開(kāi)發(fā)利用的干擾較少，水位變化主要隨降雨量的增減而升降，并具明顯的季節(jié)性（龐園等，2017）。

1.2 數(shù)據(jù)來(lái)源與檢測(cè)方法

本研究共布設(shè)采樣點(diǎn)80個(gè)（圖1），其中松散巖類孔隙水采樣點(diǎn)28個(gè)，井深范圍0～10 m；基巖裂隙水采樣點(diǎn)20個(gè)，井深范圍34～95 m；巖溶水采樣點(diǎn)32個(gè)，井深范圍23～100 m。松散巖類孔隙水采樣點(diǎn)均利用居民生活用井，基巖裂隙水采樣點(diǎn)和巖溶水采樣點(diǎn)均利用廣州市水務(wù)局建造的地下水監(jiān)測(cè)專用井。所有采樣點(diǎn)的水位均隨降雨增減而升降，具有明顯的季節(jié)性特征，均屬淺層地下水。采樣時(shí)間為2012—2016年，每年均于雨季（4—9月）和旱季（10月—次年3月）各采樣1次。受不可抗力影響，未能采集全部800組水樣，實(shí)際共采集水樣775組，缺失水樣對(duì)本文結(jié)論影響不大。為使地下水水樣更具代表性，每次采樣前均用潛水泵抽干井內(nèi)的水并直至水位恢復(fù)原狀或?qū)τ克枯^大的井進(jìn)行長(zhǎng)時(shí)間抽水使井內(nèi)的地下水得到充分循環(huán)更新。所有水樣均委托廣州市二次供水技術(shù)咨詢服務(wù)中心檢測(cè)，檢測(cè)方法均根據(jù)《地下水質(zhì)量標(biāo)準(zhǔn)》

圖1 研究區(qū)土地利用類型及地下水監(jiān)測(cè)井位置Fig. 1 Location of groundwater monitoring wells and land use condition in Guanghua Basin

（GB/T14848-93）按《生活飲用水標(biāo)準(zhǔn)檢驗(yàn)法》（GB5750—2006）執(zhí)行。

1.3 評(píng)價(jià)方法

以《地下水質(zhì)量標(biāo)準(zhǔn)》（GB/T14848-93）的地下水質(zhì)量分類標(biāo)準(zhǔn)作為評(píng)價(jià)依據(jù)，對(duì)氨氮、硝酸鹽氮和亞硝酸鹽氮進(jìn)行分類，以Ⅲ類標(biāo)準(zhǔn)限值（適用于集中式生活飲用水水源及工、農(nóng)業(yè)用水）作為超標(biāo)標(biāo)準(zhǔn)，氨氮、硝酸鹽氮、亞硝酸鹽氮的超標(biāo)標(biāo)準(zhǔn)分別為0.2、20、0.02 mg·L-1。由于部分指標(biāo)濃度低于最低檢測(cè)限值，在計(jì)算時(shí)取限值的一半，氨氮、硝酸鹽氮、亞硝酸鹽氮的最低檢測(cè)限值分別為0.02、0.15、0.001 mg·L-1。

2 結(jié)果

2.1 三氮污染水平分析

根據(jù)表1中三氮的含量范圍和平均質(zhì)量濃度可知，近5年氨氮的極大值和平均值分別是Ⅲ類標(biāo)準(zhǔn)限值的329倍和4.7倍，硝酸鹽氮的極大值和平均值分別是Ⅲ類標(biāo)準(zhǔn)限值的5.5倍和0.2倍，亞硝酸鹽氮的極大值和平均值分別是Ⅲ類標(biāo)準(zhǔn)限值的137.5倍和2.9倍，這表明廣花盆地地下水受到了氨氮、硝酸鹽氮和亞硝酸鹽氮不同程度的污染。根據(jù)表1中三氮的超標(biāo)率可以判斷污染的嚴(yán)重程度，整體而言，廣花盆地地下水中氨氮污染最嚴(yán)重，其次為亞硝酸鹽氮，硝酸鹽氮的污染相對(duì)最輕。

表1 廣花盆地地下水三氮質(zhì)量濃度統(tǒng)計(jì)表Table 1 Statistics of Three-Nitrogen mass concentrations in Guanghua basin groundwater

為進(jìn)一步探清廣花盆地不同地下水中含氮污染物的污染程度，分別統(tǒng)計(jì)了松散巖類孔隙水、基巖裂隙水和巖溶水中的氮含量，結(jié)果見(jiàn)表 2。根據(jù)統(tǒng)計(jì)結(jié)果，結(jié)合Ⅲ類標(biāo)準(zhǔn)限值可知，松散巖類孔隙水、基巖裂隙水和巖溶水均受到了含氮污染物不同程度的污染。由超標(biāo)率可知，3種地下水中含氮污染物的超標(biāo)率按由大至小排列依次為氨氮、亞硝酸鹽氮和硝酸鹽氮。由質(zhì)量濃度的均值可知，3種地下水均以硝酸鹽氮為主，其次為氨氮，亞硝酸鹽氮最小。總體而言，松散巖類孔隙水中的氮含量最高，污染最嚴(yán)重，其次為基巖裂隙水，巖溶水相對(duì)最輕。

2.2 三氮時(shí)間變化特征

根據(jù)80個(gè)采樣點(diǎn)2012—2016年雨季和旱季共10次采樣數(shù)據(jù)，分別對(duì)氨氮、硝酸鹽氮、亞硝酸鹽氮以及總氮的質(zhì)量濃度和超標(biāo)率隨時(shí)間的變化關(guān)系進(jìn)行分析，結(jié)果見(jiàn)圖2和圖3。由圖2可知，當(dāng)擬合方程的一次項(xiàng)系數(shù)為正時(shí)，表明質(zhì)量濃度隨時(shí)間呈上升趨勢(shì)，反之，則呈下降趨勢(shì)。雨季線性擬合方程的一次項(xiàng)系數(shù)均為負(fù)，表明三氮以及總氮的質(zhì)量濃度在雨季均隨時(shí)間呈下降趨勢(shì)。近5年雨季氨氮、硝酸鹽氮、亞硝酸鹽氮以及總氮的質(zhì)量濃度分別下降了 0.82、0.15、0.009 和 0.977 mg·L-1，年均下降幅度分別為0.16、0.03、0.002和0.195 mg·L-1。旱季線性擬合方程的一次項(xiàng)系數(shù)既有正也有負(fù)，表明各含氮污染物隨時(shí)間的變化趨勢(shì)有所不同。根據(jù)旱季三氮以及總氮與時(shí)間的線性擬合方程可知，氨氮、硝酸鹽氮和總氮的質(zhì)量濃度隨時(shí)間呈下降趨勢(shì)，亞硝酸鹽氮的質(zhì)量濃度隨時(shí)間呈上升趨勢(shì)。近5年旱季氨氮、硝酸鹽氮和總氮的質(zhì)量濃度分別下降了1.36、0.82和2.163 mg·L-1，年均下降幅度分別為0.27、0.16和0.433 mg·L-1；亞硝酸鹽氮的質(zhì)量濃度上升了0.012 mg·L-1，年均上升幅度為0.002 mg·L-1。對(duì)比三氮以及總氮的質(zhì)量濃度在雨季和旱季的變化速度可知，氨氮、硝酸鹽氮和總氮的質(zhì)量濃度在旱季的下降速度分別是雨季的1.7倍、5.3倍和2.2倍，亞硝酸鹽氮的質(zhì)量濃度在旱季的上升速度和雨季的下降速度相同。同時(shí)由圖2可以看出，三氮以及總氮的質(zhì)量濃度隨雨季和旱季的交替而增減，對(duì)比三氮以及總氮在雨季和旱季的質(zhì)量濃度均值可知，旱季氨氮、硝酸鹽氮、亞硝酸鹽氮以及總氮的質(zhì)量濃度分別是雨季的1.4倍、1.1倍、1.1倍和1.1倍，表明旱季三氮以及總氮的質(zhì)量濃度高于雨季。

表2 不同種類地下水中三氮質(zhì)量濃度統(tǒng)計(jì)表Table 2 Statistics of Three-Nitrogen mass concentrations in the different type of groundwater

圖2 2012—2016年地下水中三氮和總氮的質(zhì)量濃度Fig. 2 The mass concentrations of Three-Nitrogen and TN in groundwater from 2012 to 2016

圖3 2012—2016年地下水中三氮的超標(biāo)率Fig. 3 The over standard rate of Three-Nitrogen in groundwater from 2012 to 2016

同理，由圖3可知，雨季氨氮和亞硝酸鹽氮的超標(biāo)率隨時(shí)間呈下降趨勢(shì)，年均下降幅度分別為8%和4%，硝酸鹽氮的超標(biāo)率隨時(shí)間呈上升趨勢(shì)，年均上升幅度為0.4%。旱季氨氮和硝酸鹽氮的超標(biāo)率隨時(shí)間呈下降趨勢(shì)，年均下降幅度分別為 4%和1%，亞硝酸鹽氮的超標(biāo)率隨時(shí)間呈上升趨勢(shì)，年均上升幅度為 1%。同時(shí)，三氮的超標(biāo)率隨雨季和旱季的交替而增減，對(duì)比三氮在雨季和旱季的超標(biāo)率可知，雨季氨氮、亞硝酸鹽氮和硝酸鹽氮的超標(biāo)率分別是旱季的1.2倍、1.2倍和0.7倍，表明氨氮和亞硝酸鹽氮在雨季的超標(biāo)率高于旱季，而硝酸鹽氮的超標(biāo)率則在旱季較高。

2.3 三氮空間變化特征

圖4 氮的質(zhì)量濃度分布圖Fig. 4 Distribution diagram of mass concentration of nitrogen

為分析廣花盆地地下水中氮含量的空間分布特征，利用 ArcGIS的數(shù)據(jù)管理、空間分析和制圖功能分別制作了總氮、氨氮、硝酸鹽氮和亞硝酸鹽氮的質(zhì)量濃度分布圖，結(jié)果見(jiàn)圖 4。由于每個(gè)采樣點(diǎn)的地下水受污染的原因較多，且不同時(shí)間的污染嚴(yán)重程度也具有較大的隨機(jī)性，為使分析結(jié)果更加可靠，均采用質(zhì)量濃度的多年平均值進(jìn)行分析。根據(jù)統(tǒng)計(jì)結(jié)果，總氮、氨氮、硝酸鹽氮和亞硝酸鹽氮在城鄉(xiāng)居民區(qū)中的質(zhì)量濃度均值分別為 6.712、1.90、4.74和0.075 mg·L-1，分別是耕地的1.4倍、4.1倍、1.1倍和2.7倍，是林地的1.3倍、3.0倍、1.1倍和1.0倍。結(jié)合圖4以及土地利用類型可知，城鄉(xiāng)居民區(qū)中三氮和總氮的質(zhì)量濃度普遍高于其他土地利用類型區(qū)，表明城鄉(xiāng)居民區(qū)地下水中的氮污染較耕地和林地嚴(yán)重。由表2可知，氨氮在松散巖類孔隙水、基巖裂隙水和巖溶水中的質(zhì)量濃度均值分別為0.729、0.985和1.112 mg·L-1，硝酸鹽氮在松散巖類孔隙水、基巖裂隙水和巖溶水中的質(zhì)量濃度均值分別為9、2.003和1.824 mg·L-1，亞硝酸鹽氮在松散巖類孔隙水、基巖裂隙水和巖溶水中的質(zhì)量濃度分別為0.047、0.047和0.05 mg·L-1。對(duì)比3種地下水中三氮的質(zhì)量濃度大小可知，氨氮和亞硝酸鹽氮的質(zhì)量濃度在3種地下水中按從大至小排列依次是巖溶水、基巖裂隙水和松散巖類孔隙水，硝酸鹽氮的質(zhì)量濃度在3種地下水中按從大至小排列依次是松散巖類孔隙水、基巖裂隙水和巖溶水。根據(jù)廣花盆地地下水的埋藏特點(diǎn)，按由淺至深排列依次是松散巖類孔隙水、基巖裂隙水、巖溶水。這表明氨氮和亞硝酸鹽的含量隨地下水埋藏深度增加而增大，硝酸鹽氮的含量隨地下水埋藏深度增加而減小。

3 討論

3.1 降雨與氮含量的關(guān)系

研究表明（韓宇平等，2017；蔣然等，2016；陳建平等，2015；張婷等，2014），地下水中氮的濃度會(huì)受到降雨的影響。根據(jù)廣州市水資源公報(bào)統(tǒng)計(jì)了廣州市 2012—2016年雨季和旱季的降雨量，由圖5可知，雨季和旱季的降雨量均隨時(shí)間呈增大趨勢(shì)。與前文三氮的濃度隨時(shí)間的變化特征進(jìn)行對(duì)比可知，雨季總氮以及三氮的濃度均隨降雨量增大而減小；旱季總氮、氨氮和硝酸鹽氮的濃度均隨降雨量增大而減小，亞硝酸鹽氮的濃度則隨降雨量增大而增大。與前文三氮的超標(biāo)率隨時(shí)間的變化特征進(jìn)行對(duì)比可知，雨季氨氮和亞硝酸鹽的超標(biāo)率都隨降雨量增大而減小，硝酸鹽氮的超標(biāo)率則隨降雨量增大而增大；旱季氨氮和硝酸鹽氮的超標(biāo)率隨降雨量增大而減小，亞硝酸鹽氮的超標(biāo)率則隨降雨量增大而增大。

圖5 2012—2016年廣州市降雨量Fig. 5 Rainfall in Guangzhou from 2012 to 2016

為定量分析研究區(qū)降雨量與地下水中氮含量以及污染程度的關(guān)系，對(duì)降雨量與氮含量以及超標(biāo)率進(jìn)行了線性擬合，結(jié)果見(jiàn)圖6～圖9。由圖6和圖7可知，雨季氨氮、亞硝酸鹽氮和總氮濃度均與降雨量呈負(fù)相關(guān)，僅硝酸鹽氮濃度與降雨呈正相關(guān)；旱季氨氮和總氮濃度與降雨量呈負(fù)相關(guān)，硝酸鹽氮和亞硝酸鹽氮濃度與降雨量呈正相關(guān)。由圖8和圖9可知，雨季氨氮和亞硝酸鹽氮的超標(biāo)率與降雨量呈負(fù)相關(guān)，硝酸鹽氮的超標(biāo)率與降雨量呈正相關(guān)；旱季三氮的超標(biāo)率均與降雨量呈正相關(guān)。總體而言，研究區(qū)降雨量的增大會(huì)促使地下水中的氮含量減少，氮污染程度降低。

圖6 2012—2016年雨季降雨量與氮質(zhì)量濃度的線性擬合圖Fig. 6 Linear fitting chart of rainfall and nitrogen mass concentration in the rainy seasons from 2012 to 2016

圖7 2012—2016年旱季降雨量與氮質(zhì)量濃度的線性擬合圖Fig. 7 Linear fitting chart of rainfall and nitrogen mass concentration in the dry seasons from 2012 to 2016

3.2 地下水的埋藏深度與氮含量的關(guān)系

廣花盆地地下水的埋藏深度按由淺至深排列依次是松散巖類孔隙水、基巖裂隙水和巖溶水，前面已分析了3種地下水中三氮的含量和超標(biāo)率的變化特征，由此可以判斷不同埋藏深度地下水中氮含量的變化特征。總體而言，廣花盆地地下水中的氮含量隨地下水埋藏深度的增加而減少，其中氨氮和亞硝酸鹽的含量隨地下水埋藏深度增加而增大，硝酸鹽氮的含量隨地下水埋藏深度增加而減小。三氮的變化規(guī)律表現(xiàn)不一致，主要是由于隨著地下水埋藏深度的增加，三氮之間可通過(guò)硝化作用和反硝化作用相互轉(zhuǎn)化。當(dāng)?shù)叵滤h(huán)境呈氧化條件時(shí)，有利于硝化過(guò)程，氨氮在好氧生物的硝化作用下生成硝酸鹽氮和亞硝酸鹽氮。當(dāng)?shù)叵滤h(huán)境呈還原條件時(shí)，有利于反硝化過(guò)程，此時(shí)硝酸鹽氮被還原，含量減少。而在還原條件下，氨氮和亞硝酸鹽氮的性質(zhì)比較穩(wěn)定。隨著地下水埋藏深度的增加，地下水的徑流條件以及巖性結(jié)構(gòu)發(fā)生變化，導(dǎo)致地下水的氧化還原環(huán)境發(fā)生變化，間接引起地下水中氮含量的變化。

3.3 人類活動(dòng)與氮含量的關(guān)系

由圖1可知，研究區(qū)土地利用類型主要分為城鄉(xiāng)居民用地、耕地、林地和水域4類，耕地和城鄉(xiāng)居民用地面積較大，林地和水域面積較小。對(duì)比圖10可知，研究區(qū)的林地主要分布在山丘區(qū)，受人類活動(dòng)影響較小，該區(qū)域地下水中的氮含量基本不受人類活動(dòng)影響。因此，影響研究區(qū)地下水中氮含量的人類活動(dòng)主要有農(nóng)業(yè)污染、生活和工業(yè)廢污水等。

3.3.1 農(nóng)業(yè)污染

圖8 2012—2016年雨季降雨量與超標(biāo)率的線性擬合圖Fig. 8 Linear fitting chart of rainfall and over standard rate in the rainy seasons from 2012 to 2016

農(nóng)業(yè)上長(zhǎng)期施用高量氮肥是造成地下水氮污染的重要原因之一（王佳音等，2013）。化學(xué)氮肥和有機(jī)肥均含有大量的氮化物，進(jìn)入土壤后最終以氨氮形式存在，其中一部分被農(nóng)作物吸收；一部分經(jīng)硝化作用或反硝化作用等轉(zhuǎn)化為硝酸鹽氮和亞硝酸鹽氮隨入滲水進(jìn)入含水層；另一部分則吸附在土壤中，在降雨或灌溉條件下又可轉(zhuǎn)化為硝酸鹽氮或亞硝酸鹽氮而持續(xù)進(jìn)入地下水（劉曉晨等，2008；王佳音等，2013）。研究表明，氮肥的利用率較低，只有 20%～30%被農(nóng)作物吸收利用，近 40%通過(guò)灌溉水及降雨的淋溶作用滲透至含水層中（於嘉聞等，2016）。根據(jù)廣州市統(tǒng)計(jì)局?jǐn)?shù)據(jù)，2012—2016年廣州市花都區(qū)和白云區(qū)耕地面積逐年減小，與三氮的時(shí)間變化特征一致，而化肥施用量卻逐年增大，即單位面積的化肥施用量逐年增大，與三氮的時(shí)間變化特征相反（圖 11）。耕地面積減少將會(huì)使研究區(qū)農(nóng)業(yè)污染的面積減少，單位面積的化肥施用量增大將會(huì)使相應(yīng)區(qū)域地下水的氮含量增加。同時(shí)，農(nóng)業(yè)污染進(jìn)入地下水還受降雨影響。為定量分析農(nóng)業(yè)污染和降雨與耕地區(qū)域地下水中氮含量的相關(guān)性，運(yùn)用灰色關(guān)聯(lián)分析法（李雯等，2009）分別計(jì)算了研究區(qū)降雨、耕地面積以及化肥施用量與耕地區(qū)域地下水中氮含量的關(guān)聯(lián)度，計(jì)算得到降雨、耕地面積、化肥施用量的r值分別為0.51、0.84、0.51，表明耕地面積與氮含量之間的關(guān)聯(lián)度最大。因此，在化肥施用量逐年升高的情況下，受耕地面積的減少以及降雨量的增大，耕地區(qū)域地下水中的氮含量仍呈下降趨勢(shì)。

圖9 2012—2016年旱季降雨量與超標(biāo)率的線性擬合圖Fig. 9 Linear fitting chart of rainfall and over standard rate in the dry seasons from 2012 to 2016

3.3.2 廢污水的排放

圖10 研究區(qū)地貌圖Fig. 10 Geomorphic map of the study area

圖11 2012—2016年研究區(qū)耕地面積及化肥施用量Fig. 11 The area of cultivated land and the amount of fertilizer in the study area from 2012 to 2016

研究區(qū)經(jīng)濟(jì)發(fā)達(dá)，區(qū)域內(nèi)人口密度大，每年產(chǎn)生大量廢污水，且大部分廢污水被排入河道等地表水體中。據(jù)廣州市水資源公報(bào)，2016年廣州市廢污水排放量為22.89億噸，其中排入河流、湖泊和水庫(kù)等地表水體的廢污水量為16.48億噸。大量廢污水的排放導(dǎo)致地表水體的水質(zhì)極差，廣州市公布的53段重點(diǎn)整治河涌中，僅5段河涌達(dá)到或優(yōu)于Ⅴ類水體，未達(dá)Ⅴ類水質(zhì)河涌的主要污染指標(biāo)包括氨氮等。由于研究區(qū)地表水與地下水聯(lián)系緊密，大量含氮污染物很快隨河水通過(guò)越流補(bǔ)給地下水，造成地下水的污染。另?yè)?jù)廣州市水資源公報(bào)，廣州市排污管道的漏損率達(dá)到 14%，每年有大量廢污水直接滲漏進(jìn)入地下水，這也是研究區(qū)地下水中含氮污染物的重要來(lái)源。研究區(qū)的城鄉(xiāng)居民用地是產(chǎn)生和排放廢污水的主要區(qū)域，也是研究區(qū)地下水氮污染最嚴(yán)重的區(qū)域。由表 3可知，城鄉(xiāng)居民用地地下水中的氮含量明顯高于其他區(qū)域，其質(zhì)量濃度的多年平均值分別是耕地以及其他區(qū)域地下中的1.4倍和1.3倍。這說(shuō)明廢污水的排放是研究區(qū)最主要的污染源。

表3 2012—2016年不同土地利用類型地下水中的氮含量Table 3 The concentration of Nitrogen in groundwater of different land use land use condition from 2012 to 2016

4 結(jié)論

（1）近5年廣花盆地0～100 m深度地下水中的氮污染較嚴(yán)重，該區(qū)域松散巖類孔隙水、基巖裂隙水和巖溶水中主要的含氮污染物均為氨氮。總體而言，松散巖類孔隙水的水質(zhì)最差，其次為基巖裂隙水，巖溶水相對(duì)最好。

（2）研究區(qū)地下水中的三氮含量與降雨量大小關(guān)系密切，近5年廣花盆地地下水中的氮含量隨降雨量增大而減小。三氮的超標(biāo)率也隨降雨量變化而增減，總體而言，雨季地下水的水質(zhì)優(yōu)于旱季。

（3）研究區(qū)耕地區(qū)域地下水的水質(zhì)受耕地面積、化肥施用量和降雨量等因素的影響，其中與耕地面積關(guān)系最為密切，受耕地面積減少影響，該區(qū)域地下水中的氮含量隨時(shí)間呈下降趨勢(shì)。

（4）研究區(qū)地下水中的含氮污染物主要來(lái)源于廢污水的排放和農(nóng)業(yè)污染，其中廢污水的排放引起的氮污染更嚴(yán)重，導(dǎo)致研究區(qū)城鄉(xiāng)居民區(qū)地下水的水質(zhì)較其他區(qū)域差。

（5）受氧化還原環(huán)境的影響，氨氮和亞硝酸鹽氮的含量均隨地下水埋藏深度增加而增加，硝酸鹽氮和總氮的含量則隨地下水埋藏深度增加而減少。

COUTURIER M, TOMMI-MORIN G, SIROIS M, et al. 2017. Nitrogen transformations along a shallow subterranean estuary [J].Biogeosciences, 14(13): 3321-3336.

GLENN E P, JORDAN F, WAUGH W J. 2017. Phytoremediation of a Nitrogen- Contaminated Desert Soil by Native Shrubs and Microbial Processes [J]. Land Degradation & Development, 28(1): 361-369.

HUANG T L, CAO X, ZHANG Q, et al. 2014. Catalytic oxidation of high-concentration ammonia in groundwater by a naturally formed co-oxide filter film [J]. Desalination and Water Treatment, 52(7-9):1615-1623.

KI M G, KOH D C, YOON H, et al. 2015. Temporal variability of nitrate concentration in groundwater affected by intensive agricultural activities in a rural area of Hongseong, South Korea [J]. Environmental Earth Sciences, 74(7): 6147-6161.

MORAN S B, STACHELHAUS S L, KELLY R P, et al. 2014. Submarine Groundwater Discharge as a Source of Dissolved Inorganic Nitrogen and Phosphorus to Coastal Ponds of Southern Rhode Island [J].Estuaries and Coasts, 37(1): 104-118.

NILSSON L, WIDERLUND A. 2017. Tracing nitrogen cycling in mining waters using stable nitrogen isotope analysis [J]. Applied Geochemistry, 84: 41-51.

NISSIM W G, VOICU A, LABRECQUE M. 2014. Willow short-rotation coppice for treatment of polluted groundwater [J]. Ecological Engineering, 62: 102-114.

VANCE D, NEARY L, SUEKER J, et al. 2015. Nitrate Remediation: The Importance of the Advanced Nitrogen Cycle in Remedy Selection [J].Remediation Journal, 26(1): 93-105.

陳建平, 丁際豫, 吳子杰. 2015. 氮肥對(duì)地下水中氮遷轉(zhuǎn)機(jī)理研究[J]. 長(zhǎng)江科學(xué)院院報(bào), 32(2): 24-29.

陳建耀, 王亞, 張洪波, 等. 2006. 地下水硝酸鹽污染研究綜述[J]. 地理科學(xué)進(jìn)展, 25(1): 34-44.

馮海波, 董少剛, 王立新, 等. 2017. 潛水與承壓水三氮污染空間分布及成因分析——以內(nèi)蒙古托克托地區(qū)為例[J]. 中國(guó)農(nóng)村水利水電,(4): 91-96.

韓宇平, 柴奕, 劉中培, 等. 2017. 人民勝利渠灌區(qū)地下水氮含量時(shí)空分布特征研究[J]. 華北水利水電大學(xué)學(xué)報(bào)(自然科學(xué)版), (6): 1-7.

蔣然, 朱小平, 梁志宏, 等. 2016. 桂林毛村地下河水質(zhì)評(píng)價(jià)[J].水資源保護(hù), 32(5): 85-90.

李雯, 趙寶峰, 王杰. 2009. 基于灰色關(guān)聯(lián)度法的涇河徑流量變化因素分析[J]. 水資源與水工程學(xué)報(bào), 20(2): 92-94.

劉景濤, 孫繼朝, 林良俊, 等. 2011. 廣州市地下水環(huán)境三氮污染初探[J]. 中國(guó)地質(zhì), 38(2): 489-494.

劉曉晨, 孫占祥. 2008. 地下水硝態(tài)氮污染現(xiàn)狀及研究進(jìn)展[J]. 遼寧農(nóng)業(yè)科學(xué), (5): 41-45.

龐園, 張明珠, 龐志研, 等. 2017. 廣花盆地地下水水質(zhì)評(píng)價(jià)及變化趨勢(shì)分析[J]. 水電能源科學(xué), 35(11): 40-43, 35.

萬(wàn)長(zhǎng)園, 王明玉, 王慧芳, 等. 2014. 華北平原典型剖面地下水三氮污染時(shí)空分布特征[J]. 地球與環(huán)境, 42(4): 472-479.

王佳音, 張世濤, 王明玉, 等. 2013. 滇池流域大河周邊地下水氮污染的時(shí)空分布特征及影響因素分析[J]. 中國(guó)科學(xué)院研究生院學(xué)報(bào),30(3): 339-346.

許可, 陳鴻漢. 2011. 地下水中三氮污染物遷移轉(zhuǎn)化規(guī)律研究進(jìn)展[J].中國(guó)人口·資源與環(huán)境, 21(S2): 421-424.

於嘉聞, 周金龍, 曾妍妍, 等. 2016. 新疆喀什地區(qū)東部地下水三氮空間分布特征及影響因素[J]. 環(huán)境化學(xué), 35(11): 2402-2410.

張婷, 陳世儉, 傅嬌鳳. 2014. 四湖地區(qū)地下水三氮含量及時(shí)空分布特征分析[J]. 長(zhǎng)江流域資源與環(huán)境, 23(9): 1295-1300.