王世忠 張利紅 王曉慶 袁建明 黃程璐

摘要? 農(nóng)村面源污染問(wèn)題是當(dāng)前研究的熱點(diǎn)和難點(diǎn)問(wèn)題。以鲹魚(yú)河流域會(huì)東段為例，構(gòu)建了農(nóng)村面源污染強(qiáng)度與生態(tài)非敏感性強(qiáng)度的空間耦合協(xié)調(diào)度模型，并對(duì)各斑塊的農(nóng)村面源污染強(qiáng)度、生態(tài)非敏感性強(qiáng)度和空間耦合協(xié)調(diào)度進(jìn)行了計(jì)算和空間分析。結(jié)果表明：鲹魚(yú)河流域會(huì)東段雖然有60.55%的土地處于優(yōu)質(zhì)耦合協(xié)調(diào)區(qū)和良好耦合協(xié)調(diào)區(qū)，但仍然有11.96%的土地處于極度失調(diào)區(qū)，農(nóng)村面源污染帶來(lái)的生態(tài)環(huán)境問(wèn)題依然非常嚴(yán)峻；耦合協(xié)調(diào)度的熱點(diǎn)區(qū)和冷點(diǎn)區(qū)的空間分布非常明顯，99%置信區(qū)間的熱點(diǎn)區(qū)面積占56.08%，99%置信區(qū)間的冷點(diǎn)區(qū)面積占21.21%；養(yǎng)殖業(yè)污染是鲹魚(yú)河流域會(huì)東段農(nóng)村主要的面源污染源，是影響空間耦合協(xié)調(diào)性的主要因素。建議各鄉(xiāng)鎮(zhèn)科學(xué)地劃定畜牧業(yè)的禁養(yǎng)限養(yǎng)區(qū)，合理確定養(yǎng)殖規(guī)模，推行種養(yǎng)結(jié)合和生態(tài)養(yǎng)殖模式；科學(xué)布設(shè)生態(tài)田埂、生態(tài)溝渠和生態(tài)池塘等系統(tǒng)，以防止農(nóng)業(yè)尾水和農(nóng)村生活污水未經(jīng)任何處理直接排入河道；逐步降低和減少農(nóng)村經(jīng)濟(jì)社會(huì)活動(dòng)對(duì)生態(tài)環(huán)境尤其是生態(tài)高敏感區(qū)的農(nóng)村生態(tài)環(huán)境的影響和干擾。

關(guān)鍵詞? 面源污染強(qiáng)度；生態(tài)非敏感性強(qiáng)度；空間耦合協(xié)調(diào)模型；耦合協(xié)調(diào)度

中圖分類號(hào)? X71? 文獻(xiàn)標(biāo)識(shí)碼? A? 文章編號(hào)? 0517-6611（2024）04-0062-09

doi：10.3969/j.issn.0517-6611.2024.04.013

開(kāi)放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

Spatial Coupling Coordination Analysis of Rural Non.point Source Pollution Intensity and Ecological Non.Sensitivity Intensity—Taking Huidong Section of Trevally River Basin as an Example

WANG Shi.zong1， ZHANG Li.hong2， WANG Xiao.qing2? et al

（1. School of Public Administration， Zhejiang University of Finance and Economy， Hangzhou， Zhejiang 310018; 2. East China Institute of Survey， Design and Research， Hangzhou，Zhejiang 311122）

Abstract? Rural non.point source pollution is a hot and difficult problem in current research. Taking the Huidong section of the Trevally River basin as an example， this paper constructs a spatial coupling coordination model of rural non.point source pollution intensity and ecological non sensitivity intensity， and the rural non.point source pollution intensity， ecological non sensitivity intensity and spatial coupling coordination degree of each patch were calculated and analyzed. The results show that： Although 60.55% of the land in Huidong section of Trevally River basin is in high.quality coupling coordination area and good coupling coordination area， 11.96% of the land is still in extreme imbalance area， and the ecological and environmental problems caused by rural non.point source pollution are still very serious; The spatial distribution of hot spots and cold spots of coupling coordination degree is very obvious， the hot spot area of 99% confidence interval accounts for 56.08%， and the cold spot area of 99% confidence interval accounts for 21.21%; Aquaculture pollution is the main non.point source pollution source in rural areas in Huidong section of Trevally River basin， and it is the main factor affecting the spatial coupling coordination. It suggests that each villages and towns scientifically delimit the no breeding and limited breeding areas of animal husbandry， reasonably determine the breeding scale， and promote the combination of planting and breeding and ecological breeding mode;Ecological ridge， ecological ditch and ecological pond systems shall be scientifically arranged to prevent agricultural tail water and rural domestic sewage from being directly discharged into the river without any treatment;The impact and interference of rural economic and social activities on the ecological environment， especially the rural ecological environment in ecological highly sensitive areas should be gradually reduced.

Key words? Non.point source pollution intensity;Ecological non sensitivity intensity;Spatial coupling coordination model;Coupling coordination degree

作者簡(jiǎn)介? 王世忠（1975—），男，浙江衢州人，副教授，博士，從事土地資源管理研究。

收稿日期? 2023-02-02

面源污染，即非點(diǎn)源污染，是相對(duì)于排污點(diǎn)集中、排污途徑明確的點(diǎn)源污染而言的區(qū)域環(huán)境污染問(wèn)題。根據(jù)美國(guó)清潔水法修正案的定義，所謂面源污染是指“進(jìn)入地表及地下水體的，并以廣域、分散和微量的形式存在的一種污染物”。相對(duì)于點(diǎn)源污染，面源污染具有分布隨機(jī)、污染源復(fù)雜、對(duì)土壤肥力破壞性大、對(duì)水源造成污染后防治較難等特點(diǎn)［1］。農(nóng)業(yè)面源污染已經(jīng)成為水體污染的重要來(lái)源之一［2］。從農(nóng)業(yè)面源污染研究方法上看，SWAT（soil and water assessment tool）、AGNPS（agricultural non-point source）和HSPF（hydrological simulation program-fortran）是目前應(yīng)用較多且較為成熟的面源污染模擬模型［3］。從農(nóng)業(yè)面源污染治理研究上，BMPs（Best Manangement Practices）框架體系［4］，即在污染物進(jìn)入水體前，通過(guò)各種經(jīng)濟(jì)高效、工藝簡(jiǎn)單、滿足生態(tài)環(huán)境要求、適應(yīng)污染特性的措施使其得到有效的控制，是目前最常見(jiàn)的治理研究框架體系。近年來(lái)，我國(guó)學(xué)者不斷提出適合我國(guó)面源污染實(shí)際情況的理論和方法體系，并基于SWAT模型、模糊二層多目標(biāo)規(guī)劃模型（EC-IFBLMOP模型）、BMPs技術(shù)、4R技術(shù)體系，分析了我國(guó)面源污染產(chǎn)生的根本原因并進(jìn)行了控制預(yù)防方案的模擬研究［5-8］。但如何治理我國(guó)的面源污染問(wèn)題，減輕各大流域水源污染富營(yíng)養(yǎng)嚴(yán)重的現(xiàn)象，仍是目前研究爭(zhēng)論和探討的熱點(diǎn)。

在前人研究的基礎(chǔ)上，以鲹魚(yú)河流域會(huì)東段為例，通過(guò)構(gòu)建空間耦合協(xié)調(diào)度模型，進(jìn)行了鲹魚(yú)河流域會(huì)東段農(nóng)村面源污染強(qiáng)度與生態(tài)非敏感性強(qiáng)度的空間耦合協(xié)調(diào)度研究，以期為農(nóng)村面源污染提供新的研究思路和研究方法。鲹魚(yú)河流域會(huì)東段，作為一個(gè)具有區(qū)域特色的局部地區(qū)，屬于云貴高原邊緣地帶，地處金沙江的上游，山高地陡，“一山有四季”，生態(tài)環(huán)境比較惡劣，農(nóng)村面源污染問(wèn)題對(duì)當(dāng)?shù)刈匀画h(huán)境和社會(huì)經(jīng)濟(jì)環(huán)境的影響比較明顯。尤其是烏東德電站并網(wǎng)發(fā)電后，金沙江的流速和流態(tài)均發(fā)生了變化，改變了水環(huán)境的物理和化學(xué)條件，降低了污染物在水體中的稀釋、降解擴(kuò)散和轉(zhuǎn)化等過(guò)程。因此，很有必要對(duì)金沙江的上游鲹魚(yú)河流域會(huì)東段的農(nóng)村面源污染問(wèn)題進(jìn)行系統(tǒng)研究，以期進(jìn)一步加強(qiáng)會(huì)東縣的農(nóng)村面源污染治理，優(yōu)化金沙江流域的空間規(guī)劃布局。

1? 研究方法、研究區(qū)域與數(shù)據(jù)來(lái)源

1.1? 研究方法

1.1.1? 農(nóng)村污染物排放量估計(jì)。

（1）農(nóng)村生活污水污染排放量估計(jì)。農(nóng)村生活污水污染排放量計(jì)算公式如下：

鄉(xiāng)村生活污水污染排放量=鄉(xiāng)村人口總數(shù)×農(nóng)村生活污水排放系數(shù)×污水平均含量×入河系數(shù)（1）

由于數(shù)據(jù)缺失，可以參照同樣是西南山區(qū)的重慶市調(diào)查結(jié)果［9］。根據(jù)抽樣調(diào)查結(jié)果顯示，重慶市農(nóng)村人口人均每天生活污水排放量為0.67 L/d，參考重慶市環(huán)境監(jiān)測(cè)中心的監(jiān)測(cè)結(jié)果，COD、BOD5、TN、TP分別取292.69、138.33、44.14、4.49 mg/L，入河系數(shù)取0.30。

（2）農(nóng)村生活垃圾污染排放量估計(jì)。農(nóng)村生活垃圾排放量計(jì)算公式如下：

鄉(xiāng)村生活垃圾污染排放量=鄉(xiāng)村人口總數(shù)×農(nóng)村生活垃圾排放系數(shù)×垃圾滲漏液平均含量×入河系數(shù)（2）

由于數(shù)據(jù)缺失，可以參照同樣是西南山區(qū)的重慶市調(diào)查結(jié)果［9］。根據(jù)抽樣調(diào)查結(jié)果顯示，重慶市農(nóng)村人口人均生活垃圾排放量為0.67 kg/d，其COD、BOD5、TN、TP參考垃圾滲濾液，分別取50.00、5.00、1.00、0.20 mg/kg，入河系數(shù)取0.20。

（3）種植業(yè)化肥的污染物排放量估計(jì)。

根據(jù)會(huì)東縣的農(nóng)村統(tǒng)計(jì)年鑒可知，會(huì)東縣的農(nóng)用化肥主要是氮肥、磷肥、鉀肥和復(fù)合肥4種。化肥的農(nóng)業(yè)面源污染產(chǎn)生量根據(jù)輸出系數(shù)法進(jìn)行估算，計(jì)算公式如下：

化肥污染物產(chǎn)生量=化肥用量×化肥產(chǎn)污系數(shù)（3）

化肥潛在污染量=化肥污染物產(chǎn)生量×（1-化肥利用率）=化肥用量×化肥產(chǎn)污系數(shù)×（1-化肥利用率）（4）

化肥的面源污染物負(fù)荷量=化肥污染物產(chǎn)生量×化肥面源污染排放系數(shù)（5）

式中：年鑒中化肥施用量采用折純量，根據(jù)化學(xué)組成成分分析，氮肥、磷肥和復(fù)合肥（氮磷鉀含量相同）的總氮（TN）產(chǎn)污系數(shù)分別為1、0和0.33，相應(yīng)的總磷（TP）產(chǎn)污系數(shù)分別為0、0.44和0.15［10］。鲹魚(yú)河流域會(huì)東段的化肥利用率采用2017年全國(guó)化肥的平均利用率36.50%［11］。根據(jù)馬國(guó)霞等［12］給出的全國(guó)種植業(yè)排放的NH3-N量（氨氮含量指標(biāo)）約占TN排放量的8.3%來(lái)估算化肥產(chǎn)生的NH3-N量。由于數(shù)據(jù)缺失，可以參照同樣是西南山區(qū)的重慶市調(diào)查結(jié)果。根據(jù)重慶市農(nóng)業(yè)環(huán)境監(jiān)測(cè)站調(diào)查研究，氮、磷肥的入河系數(shù)分別為0.100 7和0.059 9［9］。另外，以氮肥、磷肥的折純量分別計(jì)算的TN、TP的排放量，不考慮化肥的化學(xué)需氧量（COD）、生化需氧量（BOD5）的排放量。

（4）種植業(yè)秸稈的污染物排放量估計(jì)。

根據(jù)會(huì)東縣的農(nóng)村統(tǒng)計(jì)年鑒可知，會(huì)東縣的主要農(nóng)作物包括：谷物（稻谷、玉米、小麥、大麥、燕麥、蕎子等）、豆類（大豆、胡豆、豌豆等）、薯類（洋芋、紅笤等）、油料作物（花生、油菜籽、芝麻、向日葵等）、甘蔗、煙葉、蔬菜、瓜果等。秸稈的農(nóng)業(yè)面源污染產(chǎn)生量的計(jì)算公式如下所示：

秸稈污染物產(chǎn)生量=作物產(chǎn)量×秸稈與作物產(chǎn)量比×秸稈產(chǎn)污系數(shù)（6）

秸稈潛在污染量=秸稈污染物產(chǎn)生量×（1-秸稈綜合利用率）（7）

秸稈的面源污染物負(fù)荷量=秸稈潛在污染量×秸稈面源污染排放系數(shù)（8）

式中：秸稈與作物產(chǎn)量比可以參照《非常規(guī)飼料資源的開(kāi)發(fā)與利用》研究組（1996）的研究成果［13］，具體為：

水稻0.97，小麥1.03，玉米1.37，高粱1.44，谷子1.51，其他雜糧1.60，大豆1.71，薯類0.61，

花生1.52，油菜3.00，芝麻0.64，向日葵0.60，棉花3.00，麻類1.70，甘蔗0.25，蔬菜、瓜果0.10。

各類作物秸稈的COD、TN、TP產(chǎn)污系數(shù)如表1所示［10］。秸稈綜合利用率可以通過(guò)實(shí)地調(diào)研得到，秸稈資源利用主要以秸稈能源化、飼料化、秸稈還田等為主，根據(jù)《四川省秸稈綜合利用規(guī)劃（2016—2020）》，到2020年，秸稈還田量占秸稈資源量的44.20%。秸稈作為肥料還田時(shí)的COD、TN、TP排放系數(shù)（入河系數(shù)）分別取20%、10%、5%［11］。

（5）養(yǎng)殖業(yè)的污染物排放量分析。

養(yǎng)殖業(yè)的農(nóng)業(yè)面源污染產(chǎn)生量根據(jù)輸出系數(shù)法進(jìn)行估算，具體的計(jì)算公式如下所示：

畜禽養(yǎng)殖產(chǎn)污量=畜禽飼養(yǎng)量×產(chǎn)污系數(shù) =畜禽糞便量×各污染物含量（9）

畜禽糞便量=畜禽飼養(yǎng)量×飼養(yǎng)周期×糞尿排泄系數(shù)（10）

畜禽養(yǎng)殖面源污染負(fù)荷量=畜禽養(yǎng)殖產(chǎn)污量×畜禽養(yǎng)殖面源污染排放系數(shù)（11）

式中：豬、奶牛、肉牛、雞的糞尿排泄系數(shù)分別為3.54、31.39、20.42、0.18 kg/d；羊、馬、騾的糞尿排泄系數(shù)分別為0.87、5.90、5.00 kg/d［14］。根據(jù)已有的研究文獻(xiàn)［11］，可得不同畜禽種類的產(chǎn)污系數(shù)，其中：豬的平均飼養(yǎng)周期取157 d，其他畜禽的飼養(yǎng)周期均以年計(jì)算，具體如表2所示。

規(guī)?；笄蒺B(yǎng)殖雖然在產(chǎn)生源頭上體現(xiàn)為點(diǎn)源污染，但由于其處于廣大的農(nóng)村區(qū)域，就環(huán)境影響上實(shí)際體現(xiàn)為面源污染。另外，由于不同的養(yǎng)殖方式導(dǎo)致處理糞便的方式、配套措施、飼養(yǎng)環(huán)境等會(huì)有所不同，參考已有的研究［11］，可得不同養(yǎng)殖方式下的不同畜禽種類的排放系數(shù)，具體如表3～4所示。鲹魚(yú)河流域會(huì)東段的養(yǎng)殖業(yè)采用規(guī)模養(yǎng)殖場(chǎng)、養(yǎng)殖專業(yè)戶和放養(yǎng)的比例，可以由現(xiàn)場(chǎng)實(shí)地調(diào)研獲得。

1.1.2? 農(nóng)村面源污染強(qiáng)度與生態(tài)非敏感性強(qiáng)度的空間耦合協(xié)調(diào)度模型。

1.1.2.1? 空間耦合協(xié)調(diào)作用機(jī)理。

隨著農(nóng)村經(jīng)濟(jì)社會(huì)的發(fā)展，人類活動(dòng)的頻次和范圍的加大，必然給農(nóng)村帶來(lái)生活和農(nóng)業(yè)生產(chǎn)的面源污染問(wèn)題。農(nóng)村的生活面源污染主要是生活污水污染和生活垃圾污染，農(nóng)業(yè)生產(chǎn)污染主要是種植業(yè)的化肥污染、秸稈污染以及養(yǎng)殖業(yè)的污染。綠水青山，就是金山銀山。在發(fā)展農(nóng)村經(jīng)濟(jì)的同時(shí)，必須考慮農(nóng)村生態(tài)環(huán)境的敏感性，時(shí)刻保護(hù)農(nóng)村生態(tài)環(huán)境。農(nóng)村生態(tài)敏感地區(qū)主要是水源涵養(yǎng)功能區(qū)、水土保持功能區(qū)、生物多樣性保護(hù)區(qū)以及生態(tài)脆弱區(qū)。因此，人類活動(dòng)的時(shí)候，要盡可能地減少或者回避在生態(tài)敏感區(qū)的活動(dòng)，盡可能地在生態(tài)非敏感區(qū)進(jìn)行生活和生產(chǎn)活動(dòng)。農(nóng)村經(jīng)濟(jì)與農(nóng)村生態(tài)相互耦合協(xié)調(diào)的過(guò)程就是農(nóng)村生活生產(chǎn)與農(nóng)村生態(tài)保護(hù)兩者相互影響、相互制約的過(guò)程。

1.1.2.2? 空間耦合協(xié)調(diào)度模型［15］。

假設(shè)用P代表農(nóng)村面源污染強(qiáng)度模塊的m個(gè)指標(biāo)，則可以表示為：P=（p1，p2，…，pm）；用S代表生態(tài)非敏感性強(qiáng)度模塊的n個(gè)指標(biāo)，則可以表示為：S=（s1，s2，…，sn），對(duì)應(yīng)指標(biāo)標(biāo)準(zhǔn)化后的2個(gè)模塊標(biāo)準(zhǔn)化值可以用P′=（p′1，p′2，…，p′m），S′=（s′1，s′2，…，s′n）來(lái)表示。那么，可用函數(shù)g（P）=mjrjp′j表示農(nóng)村面源污染強(qiáng)度模塊的總體評(píng)價(jià)，用函數(shù)f（S）=niuis′i表示生態(tài)非敏感性強(qiáng)度的總體評(píng)價(jià)。其中，rj表示農(nóng)村面源污染強(qiáng)度模塊的第j個(gè)指標(biāo)的權(quán)重，ui表示生態(tài)非敏感性強(qiáng)度模塊的第i個(gè)指標(biāo)的權(quán)重?？紤]數(shù)據(jù)間不同量綱和數(shù)量級(jí)對(duì)測(cè)算結(jié)果的影響，需要對(duì)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。

對(duì)于2個(gè)系統(tǒng)的耦合度C可以表示為：

C=2f（S）×g（P）f（S）+g（P）（12）

對(duì)于2個(gè)系統(tǒng)的發(fā)展度模型為：

T=αf（S）+βg（P）=0.5f（S）+0.5g（P）（13）

對(duì)于2個(gè)系統(tǒng)的耦合協(xié)調(diào)度模型為：

D=C×T（14）

1.1.2.3? 指標(biāo)體系的選取。

按照科學(xué)性、全局性、對(duì)比性、可獲得性、精確性等原則來(lái)選取指標(biāo)。農(nóng)村面源污染強(qiáng)度采用：?jiǎn)挝幻娣eNH3-N污染量（kg/km2）、單位面積COD污染量（kg/km2）、單位面積TN污染量（kg/km2）、單位面積TP污染量（kg/km2）、單位面積BOD5污染量（kg/km2）5個(gè)指標(biāo)。農(nóng)村生態(tài)非敏感性強(qiáng)度采用：水源涵養(yǎng)功能重要性（極重要、重要、一般）、水土保持功能重要性（極重要、重要、一般）、生物多樣性保護(hù)重要性（極重要、重要、一般）、生態(tài)脆弱性（極脆弱、脆弱、一般）4個(gè)指標(biāo)，其中：重要性是“極重要”的取值1，重要性是“重要”的取值2，重要性是“一般”的取值3，生態(tài)脆弱性是“極脆弱”的取值1，生態(tài)脆弱性是“脆弱”的取值2，生態(tài)脆弱性是“一般”的取值3。

1.1.2.4? 指標(biāo)權(quán)重的確定。

通過(guò)熵權(quán)法來(lái)確定指標(biāo)權(quán)重，其基本思路是根據(jù)指標(biāo)變異性的大小來(lái)確定客觀權(quán)重。一般來(lái)說(shuō)，若某個(gè)指標(biāo)的信息熵Ej越小，表明指標(biāo)值變異程度越大，提供的信息量越多，在綜合評(píng)價(jià)中所能起到的作用也越大，其權(quán)重也就越大。相反，某個(gè)指標(biāo)的信息熵越大，表明指標(biāo)值變異程度越小，提供的信息量也越小，在綜合評(píng)價(jià)中所起到的作用也越小，其權(quán)重也就越小。熵權(quán)法賦權(quán)步驟如下所示：

（1）數(shù)據(jù)標(biāo)準(zhǔn)化。

將各個(gè)指標(biāo)的數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理。假設(shè)給定了k個(gè)指標(biāo)x1，x2，…，xk，其中：xi={x1，x2，…，xn}。假設(shè)對(duì)各個(gè)指標(biāo)數(shù)據(jù)標(biāo)準(zhǔn)化的值為Y1，Y2，…，Yk，那么：

Yij=xij-min（xi）max（xi）-min（xi）（14）

（2）求各指標(biāo)的信息熵。

根據(jù)信息論中信息熵的定義，一組數(shù)據(jù)的信息熵：

Ej=-1lnnni=1pijlnpij（15）

式中：pij=Yijni=1Yij，如果pij=0，則定義：limpij→0pijlnpij=0。

（3）確定各指標(biāo)權(quán)重。

根據(jù)信息熵的計(jì)算公式，計(jì)算出各個(gè)指標(biāo)的信息熵為：E1，E2，…，Ek。通過(guò)信息熵計(jì)算各指標(biāo)的權(quán)重為：

wi=1-Eik-Ei（i=1，2，…，k）（16）

1.1.3? 耦合協(xié)調(diào)度的熱點(diǎn)分析。

空間熱點(diǎn)分析，即局部聚類檢驗(yàn)（Getis-Ord Gi*），可以用來(lái)分析耦合協(xié)調(diào)度在局部地區(qū)是否存在明顯的高值集聚特征和低值集聚特征，即各個(gè)土地斑塊單元的耦合協(xié)調(diào)度的空間集聚模式，其具體原理如下：通過(guò)計(jì)算某個(gè)土地斑塊要素及其給定距離范圍內(nèi)相鄰?fù)恋匕邏K要素的局部總和，將計(jì)算結(jié)果與該區(qū)域范圍內(nèi)所有要素的總和進(jìn)行比較，用于分析該土地斑塊要素的屬性值在局部空間水平上的集聚程度，其結(jié)果可以體現(xiàn)出土地斑塊規(guī)模的空間分異情況，即在局部地區(qū)是否存在“冷點(diǎn)”區(qū)和“熱點(diǎn)”區(qū)［16-17］。其數(shù)學(xué)模型如下：

G*i（d）=nj=1wij（d）xj/nj=1xj（17）

式中：G*i（d）是指待測(cè)算土地斑塊要素i的G*值；n是土地斑塊總數(shù)；wij（d）為距離d范圍內(nèi)的空間權(quán)重；xj是第j個(gè)土地斑塊屬性值，若該土地斑塊單元（j）與待測(cè)算土地斑塊單元（i）之間的距離小于臨界距離d，則其空間權(quán)重矩陣wij（d）為1，否則為0。

1.2? 研究區(qū)域及數(shù)據(jù)來(lái)源

鲹魚(yú)河流域位于會(huì)東縣的中西部地區(qū)，而會(huì)東縣隸屬?zèng)錾揭妥遄灾沃?，位于四川省涼山彝族自治州南端?02°13′00″～103°3′15″E，26°12′～26°55′N，地處川南滇北交匯之處、云貴高原的邊緣地帶，西鄰會(huì)理縣，北接寧南縣，縣境東、南面隔金沙江與云南省巧家縣、昆明市東川區(qū)、祿勸彝族苗族自治縣相望。鲹魚(yú)河流域會(huì)東段的土地?cái)?shù)據(jù)來(lái)自會(huì)東縣最新的“三調(diào)”數(shù)據(jù)，其他數(shù)據(jù)則來(lái)自歷年的會(huì)東縣鄉(xiāng)鎮(zhèn)統(tǒng)計(jì)年鑒及會(huì)東縣相關(guān)政府部門所提供的資料。鲹魚(yú)河流域會(huì)東段的具體區(qū)位如圖1所示。

2? 結(jié)果與分析

2.1? 鲹魚(yú)河流域農(nóng)村面源污染及污染強(qiáng)度分析

2.1.1? 農(nóng)村生活污染排放量。

農(nóng)村生活污染排放量包括農(nóng)村生活污水污染排放量和農(nóng)村生活垃圾污染排放量，根據(jù)估算結(jié)果，可得鲹魚(yú)河流域會(huì)東段各鄉(xiāng)鎮(zhèn)鄉(xiāng)村生活污染排放量的匯總估計(jì)值，具體如表5所示。

2.1.2? 農(nóng)業(yè)的面源污染物排放量。

2.1.2.1? 種植業(yè)的污染物排放量分析。

（1）化肥。

根據(jù)會(huì)東縣的農(nóng)村統(tǒng)計(jì)年鑒可知，會(huì)東縣的農(nóng)用化肥主要是氮肥、磷肥、鉀肥和復(fù)合肥4種。根據(jù)計(jì)算可得鲹魚(yú)河流域會(huì)東段各種化肥的面源污染情況，具體如表6所示。

（2）秸稈。

計(jì)算得到鲹魚(yú)河流域會(huì)東段各鄉(xiāng)鎮(zhèn)秸稈面源污染情況的匯總結(jié)果，具體如表7所示。

（3）種植業(yè)的面源污染情況估算。

鲹魚(yú)河流域會(huì)東段種植業(yè)的面源污染主要包括化肥的使用導(dǎo)致的面源污染以及種植過(guò)程中產(chǎn)生的秸稈返田利用產(chǎn)生的面源污染。根據(jù)估算結(jié)果，進(jìn)行匯總，可得到鲹魚(yú)河流域會(huì)東段各鄉(xiāng)鎮(zhèn)的種植業(yè)的面源污染情況，具體如表8～10所示。

由表8～10可知，鲹魚(yú)河流域會(huì)東段的種植業(yè)的面源污染主要有NH3-N、COD、TN和TP 4種，其中由化肥導(dǎo)致的面源污染主要有NH3-N、TN和TP 3種，由秸稈導(dǎo)致的面源污染主要有COD、TN和TP 3種；污染物TP主要由化肥導(dǎo)致的面源污染產(chǎn)生，除鐵柳鎮(zhèn)（72.16%）外，由化肥導(dǎo)致的面源污染物TP量占全部種植業(yè)面源污染物TP量的91.00%以上；污染物TN主要由化肥導(dǎo)致的面源污染產(chǎn)生，除鐵柳鎮(zhèn)（31.03%）外，由化肥導(dǎo)致的面源污染物TN量占全部種植業(yè)面源污染物TN量的67.00%以上。

2.1.2.2? 養(yǎng)殖業(yè)的污染物排放量分析。

根據(jù)會(huì)東縣的農(nóng)村統(tǒng)計(jì)年鑒（2016年）可知，會(huì)東縣的養(yǎng)殖業(yè)主要包括生豬、牛（包括肉牛、役用牛、奶牛）、馬、騾、羊（山羊、綿羊）、家禽（雞、鴨、鵝）、兔子等，其中：豬、肉牛、家禽的飼養(yǎng)量為當(dāng)年的出欄量，役用牛、羊、馬、騾的飼養(yǎng)量為年末存欄量。根據(jù)前文公式和參數(shù)數(shù)據(jù)，可計(jì)算得到鲹魚(yú)河流域會(huì)東段各鄉(xiāng)鎮(zhèn)養(yǎng)殖業(yè)的面源污染情況的匯總結(jié)果，具體如表11所示。

2.1.2.3? 農(nóng)業(yè)的污染物排放量分析。

鲹魚(yú)河流域會(huì)東段的農(nóng)業(yè)面源污染，主要包括種植業(yè)的面源污染和養(yǎng)殖業(yè)的面源污染，根據(jù)估算進(jìn)行匯總，可得各鄉(xiāng)鎮(zhèn)的農(nóng)業(yè)面源污染情況，具體如表12～14所示。

由表14可知，鲹魚(yú)河流域會(huì)東段的農(nóng)業(yè)面源污染主要由養(yǎng)殖業(yè)的面源污染產(chǎn)生，養(yǎng)殖業(yè)的面源污染量占所有農(nóng)業(yè)面源污染量的70%以上，而有些鄉(xiāng)鎮(zhèn)的養(yǎng)殖業(yè)面源污染量甚至占所有農(nóng)業(yè)面源污染量的90%以上。

2.1.3? 農(nóng)村的面源污染物排放量。

鲹魚(yú)河流域會(huì)東段的農(nóng)村面源污染排放量包括農(nóng)村生活污染排放量和農(nóng)業(yè)面源污染排放量，把農(nóng)村的生活污染排放量和農(nóng)業(yè)的面源污染排放量的計(jì)算單位進(jìn)行統(tǒng)一核算和匯總，可得到農(nóng)村的面源污染的總排放量，具體如表15～17所示。由表15可知，農(nóng)村面源污染絕大部分由農(nóng)業(yè)面源污染產(chǎn)生。

2.1.4? 農(nóng)村的面源污染物排放強(qiáng)度。

以各鄉(xiāng)鎮(zhèn)的面積為分母，以各鄉(xiāng)鎮(zhèn)的面源污染物數(shù)量作為分子，得到的值作為各鄉(xiāng)鎮(zhèn)農(nóng)村的面源污染強(qiáng)度。根據(jù)信息熵和權(quán)重的計(jì)算公式可知各指標(biāo)的信息熵和權(quán)重如表18所示。根據(jù)權(quán)重可計(jì)算得到各鄉(xiāng)鎮(zhèn)的綜合污染強(qiáng)度，具體如表19所示。

2.2? 鲹魚(yú)河流域農(nóng)村生態(tài)非敏感性強(qiáng)度分析

根據(jù)會(huì)東縣相關(guān)政府部門提供的農(nóng)村生態(tài)敏感地區(qū)（水源涵養(yǎng)功能區(qū)、水土保持功能區(qū)、生物多樣性保護(hù)區(qū)以及生態(tài)脆弱區(qū)）的空間分布圖，根據(jù)信息熵和權(quán)重的計(jì)算公式計(jì)算可知各指標(biāo)的信息熵和權(quán)重如表20所示。根據(jù)各指標(biāo)的權(quán)重，可得到各斑塊的農(nóng)村生態(tài)非敏感性強(qiáng)度的數(shù)值，采用自然斷點(diǎn)法，分為5大類，每類的具體空間布局如圖2所示。

2.3? 農(nóng)村面源污染強(qiáng)度和生態(tài)非敏感性強(qiáng)度的空間耦合協(xié)調(diào)度分析

根據(jù)耦合協(xié)調(diào)度模型可計(jì)算出鲹魚(yú)河流域會(huì)東段農(nóng)村各斑塊的耦合協(xié)調(diào)度數(shù)值，采用逯進(jìn)等［18］對(duì)耦合協(xié)調(diào)度的等級(jí)分類標(biāo)準(zhǔn)，可得出各斑塊的耦合協(xié)調(diào)度分類等級(jí)（表21），其具體的空間布局如圖3所示。

由圖3可知，鲹魚(yú)河流域會(huì)東段農(nóng)村面源污染強(qiáng)度和生態(tài)非敏感性強(qiáng)度的空間耦合協(xié)調(diào)度分類等級(jí)主要有八大類，

其中：優(yōu)質(zhì)耦合協(xié)調(diào)區(qū)面積為322.285 7 km2，占24.63%，主要

分布于鲹魚(yú)河鎮(zhèn)和嘎吉鎮(zhèn)；良好耦合協(xié)調(diào)區(qū)面積為470.111 4 km2，占35.92%；中級(jí)耦合協(xié)調(diào)區(qū)面積為151.298 9 km2，占11.56%；初級(jí)耦合協(xié)調(diào)區(qū)面積為102.411 1 km2，占7.82%；勉強(qiáng)耦合協(xié)調(diào)區(qū)面積為83.688 0 km2，占6.39%；瀕臨失調(diào)區(qū)面積為22.212 8 km2，占1.70%；輕度失調(diào)區(qū)面積為0.201 2 km2，占0.02%；極度失調(diào)區(qū)面積為156.565 3 km2，占11.96%，主要分布于野租鄉(xiāng)。

2.4? 耦合協(xié)調(diào)度的空間熱點(diǎn)分析

鲹魚(yú)河流域會(huì)東段農(nóng)村面源污染強(qiáng)度與生態(tài)非敏感性強(qiáng)度的耦合協(xié)調(diào)度的空間熱點(diǎn)分析的結(jié)果如圖4所示。由圖4可知，耦合協(xié)調(diào)度的熱點(diǎn)區(qū)和冷點(diǎn)區(qū)的空間分布非常明顯，其中99%置信區(qū)間的熱點(diǎn)區(qū)面積為733.931 1 km2，占56.08%，達(dá)到總面積的一半以上，主要分布在鲹魚(yú)河鎮(zhèn)和嘎吉鎮(zhèn)，99%置信區(qū)間的冷點(diǎn)區(qū)面積為277.596 9 km2，占21.21%，達(dá)到總面積的1/5以上，主要分布在野租鄉(xiāng)和拉馬鄉(xiāng)。

3? 結(jié)論與展望

3.1? 結(jié)論

（1）鲹魚(yú)河流域會(huì)東段雖然有60.55%的土地處于農(nóng)村面源污染強(qiáng)度和生態(tài)非敏感性強(qiáng)度的優(yōu)質(zhì)耦合協(xié)調(diào)區(qū)和良好耦合協(xié)調(diào)區(qū)，但仍然有11.96%的土地處于極度失調(diào)區(qū)，再加上鲹魚(yú)河流域會(huì)東段地處云貴高原邊緣地帶，山高地陡，山體滑坡、泥石流、水土流失等地質(zhì)災(zāi)害頻繁，石漠化問(wèn)題突出，農(nóng)村面源污染帶來(lái)的生態(tài)環(huán)境問(wèn)題依然非常嚴(yán)峻。

（2）鲹魚(yú)河流域會(huì)東段農(nóng)村面源污染強(qiáng)度與生態(tài)非敏感性強(qiáng)度的耦合協(xié)調(diào)度的熱點(diǎn)區(qū)和冷點(diǎn)區(qū)的空間分布非常明顯，其中：99%置信區(qū)間的熱點(diǎn)區(qū)面積占56.08%，占總面積的一半以上，99%置信區(qū)間的冷點(diǎn)區(qū)面積占21.21%，占總面積的1/5以上。

（3）養(yǎng)殖業(yè)污染是鲹魚(yú)河流域會(huì)東段農(nóng)村主要的面源污染源，是影響空間耦合協(xié)調(diào)性的主要因素。各鄉(xiāng)鎮(zhèn)需要根據(jù)當(dāng)?shù)刭Y源環(huán)境承載能力，結(jié)合“三生”空間規(guī)劃，科學(xué)地劃定畜牧業(yè)的禁養(yǎng)限養(yǎng)區(qū)，合理確定養(yǎng)殖規(guī)模，以地定養(yǎng)，以養(yǎng)促種，加快鲹魚(yú)河流域會(huì)東段尤其是野租鄉(xiāng)、拉馬鄉(xiāng)的畜牧業(yè)生產(chǎn)方式的轉(zhuǎn)變，推行種養(yǎng)結(jié)合和生態(tài)養(yǎng)殖模式。

（4）科學(xué)布設(shè)生態(tài)田埂、生態(tài)溝渠和生態(tài)池塘等系統(tǒng)，以防止農(nóng)業(yè)尾水和農(nóng)村生活污水未經(jīng)任何處理直接排入河道?？梢栽诤拥纼蓚?cè)或者湖庫(kù)周邊地區(qū)，對(duì)原有的田埂進(jìn)行升級(jí)處理，建設(shè)生態(tài)田??；結(jié)合農(nóng)田灌排系統(tǒng)，設(shè)計(jì)和建設(shè)生態(tài)溝渠，減少農(nóng)業(yè)和農(nóng)村生活污水對(duì)河流的破壞；利用現(xiàn)有的池塘，進(jìn)行浮葉植物的養(yǎng)殖，實(shí)現(xiàn)對(duì)農(nóng)業(yè)污染物和農(nóng)村生活污水的吸收。

（5）逐步降低和減少農(nóng)村經(jīng)濟(jì)社會(huì)活動(dòng)對(duì)生態(tài)環(huán)境尤其是生態(tài)高敏感區(qū)的農(nóng)村生態(tài)環(huán)境的影響和干擾。根據(jù)會(huì)東縣“三調(diào)”數(shù)據(jù)可知，鲹魚(yú)河流域會(huì)東段生態(tài)保護(hù)極重要區(qū)內(nèi)有耕地7.113 3 km2，園地0.285 3 km2，草地35.724 8 km2，農(nóng)村宅基地0.226 1 km2，而生態(tài)保護(hù)極重要區(qū)就是生態(tài)高度敏感區(qū)，隨著農(nóng)村經(jīng)濟(jì)社會(huì)的發(fā)展，農(nóng)村的生產(chǎn)、生活活動(dòng)所產(chǎn)生的面源污染，極大地影響了當(dāng)?shù)剞r(nóng)村的生態(tài)環(huán)境，尤其是嚴(yán)重影響和干擾了生態(tài)高度敏感區(qū)的農(nóng)村生態(tài)環(huán)境安全。

3.2? 研究展望

在前人研究的基礎(chǔ)上，通過(guò)構(gòu)建耦合協(xié)調(diào)度模型，對(duì)鲹魚(yú)河流域會(huì)東段農(nóng)村面源污染強(qiáng)度和生態(tài)非敏感性強(qiáng)度的空間耦合協(xié)調(diào)度進(jìn)行了研究。農(nóng)村的面源污染問(wèn)題研究是個(gè)很復(fù)雜的系統(tǒng)工程，其中所涉及的大量參數(shù)、系數(shù)都需要通過(guò)實(shí)地調(diào)研和長(zhǎng)期的野外數(shù)據(jù)監(jiān)測(cè)和積累獲得，由于受到研究條件限制，該研究很多參數(shù)和系數(shù)都是參照前人的研究成果，必然對(duì)研究的精確度有一定程度的影響，這是需要進(jìn)一步研究和拓展的方向。

參考文獻(xiàn)

［1］李秀芬，朱金兆，顧曉君，等.農(nóng)業(yè)面源污染現(xiàn)狀與防治進(jìn)展［J］.中國(guó)人口·資源與環(huán)境，2010，20（4）：81-84.

［2］閔繼勝，孔祥智.我國(guó)農(nóng)業(yè)面源污染問(wèn)題的研究進(jìn)展［J］.華中農(nóng)業(yè)大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2016（2）：59-66.

［3］袁一斌，毛萍，昝曉輝，等.農(nóng)業(yè)非點(diǎn)源污染SWAT模型研究態(tài)勢(shì)及研究前沿知識(shí)圖譜［J］.應(yīng)用與環(huán)境生物學(xué)報(bào)，2018，24（5）：1050-1057.

［4］倉(cāng)恒瑾，許煉峰，李志安，等.農(nóng)業(yè)非點(diǎn)源污染控制中的最佳管理措施及其發(fā)展趨勢(shì)［J］.生態(tài)科學(xué)，2005，24（2）：173-177.

［5］LIU R M，ZHANG P P，WANG X J，et al.Cost.effectiveness and cost.benefit analysis of BMPs in controlling agricultural nonpoint source pollution in China based on the SWAT model［J］.Environmental monitoring and assessment，2014，186（12）：9011-9022.

［6］CAI Y P，RONG Q Q，YANG Z F，et al.An export coefficient based inexact fuzzy bi.level multi.objective programming model for the management of agricultural nonpoint source pollution under uncertainty［J］.Journal of hydrology，2018，557：713-725.

［7］耿潤(rùn)哲，梁璇靜，殷培紅，等.面源污染最佳管理措施多目標(biāo)協(xié)同優(yōu)化配置研究進(jìn)展［J］.生態(tài)學(xué)報(bào)，2019，39（8）：2667-2675.

［8］朱金格，張曉姣，劉鑫，等.生態(tài)溝-濕地系統(tǒng)對(duì)農(nóng)田排水氮磷的去除效應(yīng)［J］.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2019，38（2）：405-411.

［9］張廣納，邵景安，王金亮，等.三峽庫(kù)區(qū)重慶段農(nóng)村面源污染時(shí)空格局演變特征［J］.自然資源學(xué)報(bào)，2015，30（7）：1197-1209.

［10］賴斯蕓.非點(diǎn)源污染調(diào)查評(píng)估方法及其應(yīng)用研究［D］.北京：清華大學(xué)，2004.

［11］陶?qǐng)@，王少麗，管孝艷，等.青海省農(nóng)業(yè)面源污染源特征分析［J］.農(nóng)業(yè)工程學(xué)報(bào)，2019，35（10）：164-172.

［12］馬國(guó)霞，於方，曹東，等.中國(guó)農(nóng)業(yè)面源污染物排放量計(jì)算及中長(zhǎng)期預(yù)測(cè)［J］.環(huán)境科學(xué)學(xué)報(bào)，2012，32（2）：489-497.

［13］《非常規(guī)飼料資源的開(kāi)發(fā)與利用》研究組.非常規(guī)飼料資源的開(kāi)發(fā)與利用［M］.北京：中國(guó)農(nóng)業(yè)出版社，1996：14-18.

［14］中國(guó)農(nóng)業(yè)科學(xué)院，環(huán)境保護(hù)部南京環(huán)境科學(xué)研究所.第一次全國(guó)污染源普查畜禽養(yǎng)殖業(yè)源產(chǎn)排污系數(shù)手冊(cè)［R］.2009.

［15］姜磊，柏玲，吳玉鳴.中國(guó)省域經(jīng)濟(jì)、資源與環(huán)境協(xié)調(diào)分析：兼論三系統(tǒng)耦合公式及其擴(kuò)展形式［J］.自然資源學(xué)報(bào)，2017，32（5）：788-799.

［16］陳宗峰，李裕瑞，劉彥隨.黃土丘陵溝壑區(qū)鄉(xiāng)村聚落分布格局特征與類型［J］.農(nóng)業(yè)工程學(xué)報(bào)，2017，33（14）：266-274.

［17］王曼曼，吳秀芹，吳斌，等.鹽池北部風(fēng)沙區(qū)鄉(xiāng)村聚落空間格局演變分析［J］.農(nóng)業(yè)工程學(xué)報(bào)，2016，32（8）：260-271.

［18］逯進(jìn)，周惠民.中國(guó)省域人力資本與經(jīng)濟(jì)增長(zhǎng)耦合關(guān)系的實(shí)證分析［J］.數(shù)量經(jīng)濟(jì)技術(shù)經(jīng)濟(jì)研究，2013，30（9）：3-19，36.