農村面源污染強度與生態非敏感性強度的空間耦合協調性分析

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

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

關鍵詞? 面源污染強度；生態非敏感性強度；空間耦合協調模型；耦合協調度

中圖分類號? X71? 文獻標識碼? A? 文章編號? 0517-6611（2024）04-0062-09

doi：10.3969/j.issn.0517-6611.2024.04.013

開放科學（資源服務）標識碼（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

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

收稿日期? 2023-02-02

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

在前人研究的基礎上，以鲹魚河流域會東段為例，通過構建空間耦合協調度模型，進行了鲹魚河流域會東段農村面源污染強度與生態非敏感性強度的空間耦合協調度研究，以期為農村面源污染提供新的研究思路和研究方法。鲹魚河流域會東段，作為一個具有區域特色的局部地區，屬于云貴高原邊緣地帶，地處金沙江的上游，山高地陡，“一山有四季”，生態環境比較惡劣，農村面源污染問題對當地自然環境和社會經濟環境的影響比較明顯。尤其是烏東德電站并網發電后，金沙江的流速和流態均發生了變化，改變了水環境的物理和化學條件，降低了污染物在水體中的稀釋、降解擴散和轉化等過程。因此，很有必要對金沙江的上游鲹魚河流域會東段的農村面源污染問題進行系統研究，以期進一步加強會東縣的農村面源污染治理，優化金沙江流域的空間規劃布局。

1? 研究方法、研究區域與數據來源

1.1? 研究方法

1.1.1? 農村污染物排放量估計。

（1）農村生活污水污染排放量估計。農村生活污水污染排放量計算公式如下：

鄉村生活污水污染排放量=鄉村人口總數×農村生活污水排放系數×污水平均含量×入河系數（1）

由于數據缺失，可以參照同樣是西南山區的重慶市調查結果［9］。根據抽樣調查結果顯示，重慶市農村人口人均每天生活污水排放量為0.67 L/d，參考重慶市環境監測中心的監測結果，COD、BOD5、TN、TP分別取292.69、138.33、44.14、4.49 mg/L，入河系數取0.30。

（2）農村生活垃圾污染排放量估計。農村生活垃圾排放量計算公式如下：

鄉村生活垃圾污染排放量=鄉村人口總數×農村生活垃圾排放系數×垃圾滲漏液平均含量×入河系數（2）

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

（3）種植業化肥的污染物排放量估計。

根據會東縣的農村統計年鑒可知，會東縣的農用化肥主要是氮肥、磷肥、鉀肥和復合肥4種。化肥的農業面源污染產生量根據輸出系數法進行估算，計算公式如下：

化肥污染物產生量=化肥用量×化肥產污系數（3）

化肥潛在污染量=化肥污染物產生量×（1-化肥利用率）=化肥用量×化肥產污系數×（1-化肥利用率）（4）

化肥的面源污染物負荷量=化肥污染物產生量×化肥面源污染排放系數（5）

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

（4）種植業秸稈的污染物排放量估計。

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

秸稈污染物產生量=作物產量×秸稈與作物產量比×秸稈產污系數（6）

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

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

式中：秸稈與作物產量比可以參照《非常規飼料資源的開發與利用》研究組（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產污系數如表1所示［10］。秸稈綜合利用率可以通過實地調研得到，秸稈資源利用主要以秸稈能源化、飼料化、秸稈還田等為主，根據《四川省秸稈綜合利用規劃（2016—2020）》，到2020年，秸稈還田量占秸稈資源量的44.20%。秸稈作為肥料還田時的COD、TN、TP排放系數（入河系數）分別取20%、10%、5%［11］。

（5）養殖業的污染物排放量分析。

養殖業的農業面源污染產生量根據輸出系數法進行估算，具體的計算公式如下所示：

畜禽養殖產污量=畜禽飼養量×產污系數 =畜禽糞便量×各污染物含量（9）

畜禽糞便量=畜禽飼養量×飼養周期×糞尿排泄系數（10）

畜禽養殖面源污染負荷量=畜禽養殖產污量×畜禽養殖面源污染排放系數（11）

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

規模化畜禽養殖雖然在產生源頭上體現為點源污染，但由于其處于廣大的農村區域，就環境影響上實際體現為面源污染。另外，由于不同的養殖方式導致處理糞便的方式、配套措施、飼養環境等會有所不同，參考已有的研究［11］，可得不同養殖方式下的不同畜禽種類的排放系數，具體如表3～4所示。鲹魚河流域會東段的養殖業采用規模養殖場、養殖專業戶和放養的比例，可以由現場實地調研獲得。

1.1.2? 農村面源污染強度與生態非敏感性強度的空間耦合協調度模型。

1.1.2.1? 空間耦合協調作用機理。

隨著農村經濟社會的發展，人類活動的頻次和范圍的加大，必然給農村帶來生活和農業生產的面源污染問題。農村的生活面源污染主要是生活污水污染和生活垃圾污染，農業生產污染主要是種植業的化肥污染、秸稈污染以及養殖業的污染。綠水青山，就是金山銀山。在發展農村經濟的同時，必須考慮農村生態環境的敏感性，時刻保護農村生態環境。農村生態敏感地區主要是水源涵養功能區、水土保持功能區、生物多樣性保護區以及生態脆弱區。因此，人類活動的時候，要盡可能地減少或者回避在生態敏感區的活動，盡可能地在生態非敏感區進行生活和生產活動。農村經濟與農村生態相互耦合協調的過程就是農村生活生產與農村生態保護兩者相互影響、相互制約的過程。

1.1.2.2? 空間耦合協調度模型［15］。

假設用P代表農村面源污染強度模塊的m個指標，則可以表示為：P=（p1，p2，…，pm）；用S代表生態非敏感性強度模塊的n個指標，則可以表示為：S=（s1，s2，…，sn），對應指標標準化后的2個模塊標準化值可以用P′=（p′1，p′2，…，p′m），S′=（s′1，s′2，…，s′n）來表示。那么，可用函數g（P）=mjrjp′j表示農村面源污染強度模塊的總體評價，用函數f（S）=niuis′i表示生態非敏感性強度的總體評價。其中，rj表示農村面源污染強度模塊的第j個指標的權重，ui表示生態非敏感性強度模塊的第i個指標的權重。考慮數據間不同量綱和數量級對測算結果的影響，需要對數據進行標準化處理。

對于2個系統的耦合度C可以表示為：

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

對于2個系統的發展度模型為：

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

對于2個系統的耦合協調度模型為：

D=C×T（14）

1.1.2.3? 指標體系的選取。

按照科學性、全局性、對比性、可獲得性、精確性等原則來選取指標。農村面源污染強度采用：單位面積NH3-N污染量（kg/km2）、單位面積COD污染量（kg/km2）、單位面積TN污染量（kg/km2）、單位面積TP污染量（kg/km2）、單位面積BOD5污染量（kg/km2）5個指標。農村生態非敏感性強度采用：水源涵養功能重要性（極重要、重要、一般）、水土保持功能重要性（極重要、重要、一般）、生物多樣性保護重要性（極重要、重要、一般）、生態脆弱性（極脆弱、脆弱、一般）4個指標，其中：重要性是“極重要”的取值1，重要性是“重要”的取值2，重要性是“一般”的取值3，生態脆弱性是“極脆弱”的取值1，生態脆弱性是“脆弱”的取值2，生態脆弱性是“一般”的取值3。

1.1.2.4? 指標權重的確定。

通過熵權法來確定指標權重，其基本思路是根據指標變異性的大小來確定客觀權重。一般來說，若某個指標的信息熵Ej越小，表明指標值變異程度越大，提供的信息量越多，在綜合評價中所能起到的作用也越大，其權重也就越大。相反，某個指標的信息熵越大，表明指標值變異程度越小，提供的信息量也越小，在綜合評價中所起到的作用也越小，其權重也就越小。熵權法賦權步驟如下所示：

（1）數據標準化。

將各個指標的數據進行標準化處理。假設給定了k個指標x1，x2，…，xk，其中：xi={x1，x2，…，xn}。假設對各個指標數據標準化的值為Y1，Y2，…，Yk，那么：

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

（2）求各指標的信息熵。

根據信息論中信息熵的定義，一組數據的信息熵：

Ej=-1lnnni=1pijlnpij（15）

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

（3）確定各指標權重。

根據信息熵的計算公式，計算出各個指標的信息熵為：E1，E2，…，Ek。通過信息熵計算各指標的權重為：

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

1.1.3? 耦合協調度的熱點分析。

空間熱點分析，即局部聚類檢驗（Getis-Ord Gi*），可以用來分析耦合協調度在局部地區是否存在明顯的高值集聚特征和低值集聚特征，即各個土地斑塊單元的耦合協調度的空間集聚模式，其具體原理如下：通過計算某個土地斑塊要素及其給定距離范圍內相鄰土地斑塊要素的局部總和，將計算結果與該區域范圍內所有要素的總和進行比較，用于分析該土地斑塊要素的屬性值在局部空間水平上的集聚程度，其結果可以體現出土地斑塊規模的空間分異情況，即在局部地區是否存在“冷點”區和“熱點”區［16-17］。其數學模型如下：

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

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

1.2? 研究區域及數據來源

鲹魚河流域位于會東縣的中西部地區，而會東縣隸屬涼山彝族自治州，位于四川省涼山彝族自治州南端，102°13′00″～103°3′15″E，26°12′～26°55′N，地處川南滇北交匯之處、云貴高原的邊緣地帶，西鄰會理縣，北接寧南縣，縣境東、南面隔金沙江與云南省巧家縣、昆明市東川區、祿勸彝族苗族自治縣相望。鲹魚河流域會東段的土地數據來自會東縣最新的“三調”數據，其他數據則來自歷年的會東縣鄉鎮統計年鑒及會東縣相關政府部門所提供的資料。鲹魚河流域會東段的具體區位如圖1所示。

2? 結果與分析

2.1? 鲹魚河流域農村面源污染及污染強度分析

2.1.1? 農村生活污染排放量。

農村生活污染排放量包括農村生活污水污染排放量和農村生活垃圾污染排放量，根據估算結果，可得鲹魚河流域會東段各鄉鎮鄉村生活污染排放量的匯總估計值，具體如表5所示。

2.1.2? 農業的面源污染物排放量。

2.1.2.1? 種植業的污染物排放量分析。

（1）化肥。

根據會東縣的農村統計年鑒可知，會東縣的農用化肥主要是氮肥、磷肥、鉀肥和復合肥4種。根據計算可得鲹魚河流域會東段各種化肥的面源污染情況，具體如表6所示。

（2）秸稈。

計算得到鲹魚河流域會東段各鄉鎮秸稈面源污染情況的匯總結果，具體如表7所示。

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

鲹魚河流域會東段種植業的面源污染主要包括化肥的使用導致的面源污染以及種植過程中產生的秸稈返田利用產生的面源污染。根據估算結果，進行匯總，可得到鲹魚河流域會東段各鄉鎮的種植業的面源污染情況，具體如表8～10所示。

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

2.1.2.2? 養殖業的污染物排放量分析。

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

2.1.2.3? 農業的污染物排放量分析。

鲹魚河流域會東段的農業面源污染，主要包括種植業的面源污染和養殖業的面源污染，根據估算進行匯總，可得各鄉鎮的農業面源污染情況，具體如表12～14所示。

由表14可知，鲹魚河流域會東段的農業面源污染主要由養殖業的面源污染產生，養殖業的面源污染量占所有農業面源污染量的70%以上，而有些鄉鎮的養殖業面源污染量甚至占所有農業面源污染量的90%以上。

2.1.3? 農村的面源污染物排放量。

鲹魚河流域會東段的農村面源污染排放量包括農村生活污染排放量和農業面源污染排放量，把農村的生活污染排放量和農業的面源污染排放量的計算單位進行統一核算和匯總，可得到農村的面源污染的總排放量，具體如表15～17所示。由表15可知，農村面源污染絕大部分由農業面源污染產生。

2.1.4? 農村的面源污染物排放強度。

以各鄉鎮的面積為分母，以各鄉鎮的面源污染物數量作為分子，得到的值作為各鄉鎮農村的面源污染強度。根據信息熵和權重的計算公式可知各指標的信息熵和權重如表18所示。根據權重可計算得到各鄉鎮的綜合污染強度，具體如表19所示。

2.2? 鲹魚河流域農村生態非敏感性強度分析

根據會東縣相關政府部門提供的農村生態敏感地區（水源涵養功能區、水土保持功能區、生物多樣性保護區以及生態脆弱區）的空間分布圖，根據信息熵和權重的計算公式計算可知各指標的信息熵和權重如表20所示。根據各指標的權重，可得到各斑塊的農村生態非敏感性強度的數值，采用自然斷點法，分為5大類，每類的具體空間布局如圖2所示。

2.3? 農村面源污染強度和生態非敏感性強度的空間耦合協調度分析

根據耦合協調度模型可計算出鲹魚河流域會東段農村各斑塊的耦合協調度數值，采用逯進等［18］對耦合協調度的等級分類標準，可得出各斑塊的耦合協調度分類等級（表21），其具體的空間布局如圖3所示。

由圖3可知，鲹魚河流域會東段農村面源污染強度和生態非敏感性強度的空間耦合協調度分類等級主要有八大類，

其中：優質耦合協調區面積為322.285 7 km2，占24.63%，主要

分布于鲹魚河鎮和嘎吉鎮；良好耦合協調區面積為470.111 4 km2，占35.92%；中級耦合協調區面積為151.298 9 km2，占11.56%；初級耦合協調區面積為102.411 1 km2，占7.82%；勉強耦合協調區面積為83.688 0 km2，占6.39%；瀕臨失調區面積為22.212 8 km2，占1.70%；輕度失調區面積為0.201 2 km2，占0.02%；極度失調區面積為156.565 3 km2，占11.96%，主要分布于野租鄉。

2.4? 耦合協調度的空間熱點分析

鲹魚河流域會東段農村面源污染強度與生態非敏感性強度的耦合協調度的空間熱點分析的結果如圖4所示。由圖4可知，耦合協調度的熱點區和冷點區的空間分布非常明顯，其中99%置信區間的熱點區面積為733.931 1 km2，占56.08%，達到總面積的一半以上，主要分布在鲹魚河鎮和嘎吉鎮，99%置信區間的冷點區面積為277.596 9 km2，占21.21%，達到總面積的1/5以上，主要分布在野租鄉和拉馬鄉。

3? 結論與展望

3.1? 結論

（1）鲹魚河流域會東段雖然有60.55%的土地處于農村面源污染強度和生態非敏感性強度的優質耦合協調區和良好耦合協調區，但仍然有11.96%的土地處于極度失調區，再加上鲹魚河流域會東段地處云貴高原邊緣地帶，山高地陡，山體滑坡、泥石流、水土流失等地質災害頻繁，石漠化問題突出，農村面源污染帶來的生態環境問題依然非常嚴峻。

（2）鲹魚河流域會東段農村面源污染強度與生態非敏感性強度的耦合協調度的熱點區和冷點區的空間分布非常明顯，其中：99%置信區間的熱點區面積占56.08%，占總面積的一半以上，99%置信區間的冷點區面積占21.21%，占總面積的1/5以上。

（3）養殖業污染是鲹魚河流域會東段農村主要的面源污染源，是影響空間耦合協調性的主要因素。各鄉鎮需要根據當地資源環境承載能力，結合“三生”空間規劃，科學地劃定畜牧業的禁養限養區，合理確定養殖規模，以地定養，以養促種，加快鲹魚河流域會東段尤其是野租鄉、拉馬鄉的畜牧業生產方式的轉變，推行種養結合和生態養殖模式。

（4）科學布設生態田埂、生態溝渠和生態池塘等系統，以防止農業尾水和農村生活污水未經任何處理直接排入河道。可以在河道兩側或者湖庫周邊地區，對原有的田埂進行升級處理，建設生態田埂；結合農田灌排系統，設計和建設生態溝渠，減少農業和農村生活污水對河流的破壞；利用現有的池塘，進行浮葉植物的養殖，實現對農業污染物和農村生活污水的吸收。

（5）逐步降低和減少農村經濟社會活動對生態環境尤其是生態高敏感區的農村生態環境的影響和干擾。根據會東縣“三調”數據可知，鲹魚河流域會東段生態保護極重要區內有耕地7.113 3 km2，園地0.285 3 km2，草地35.724 8 km2，農村宅基地0.226 1 km2，而生態保護極重要區就是生態高度敏感區，隨著農村經濟社會的發展，農村的生產、生活活動所產生的面源污染，極大地影響了當地農村的生態環境，尤其是嚴重影響和干擾了生態高度敏感區的農村生態環境安全。

3.2? 研究展望

在前人研究的基礎上，通過構建耦合協調度模型，對鲹魚河流域會東段農村面源污染強度和生態非敏感性強度的空間耦合協調度進行了研究。農村的面源污染問題研究是個很復雜的系統工程，其中所涉及的大量參數、系數都需要通過實地調研和長期的野外數據監測和積累獲得，由于受到研究條件限制，該研究很多參數和系數都是參照前人的研究成果，必然對研究的精確度有一定程度的影響，這是需要進一步研究和拓展的方向。

參考文獻

［1］李秀芬，朱金兆，顧曉君，等.農業面源污染現狀與防治進展［J］.中國人口·資源與環境，2010，20（4）：81-84.

［2］閔繼勝，孔祥智.我國農業面源污染問題的研究進展［J］.華中農業大學學報（社會科學版），2016（2）：59-66.

［3］袁一斌，毛萍，昝曉輝，等.農業非點源污染SWAT模型研究態勢及研究前沿知識圖譜［J］.應用與環境生物學報，2018，24（5）：1050-1057.

［4］倉恒瑾，許煉峰，李志安，等.農業非點源污染控制中的最佳管理措施及其發展趨勢［J］.生態科學，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］耿潤哲，梁璇靜，殷培紅，等.面源污染最佳管理措施多目標協同優化配置研究進展［J］.生態學報，2019，39（8）：2667-2675.

［8］朱金格，張曉姣，劉鑫，等.生態溝-濕地系統對農田排水氮磷的去除效應［J］.農業環境科學學報，2019，38（2）：405-411.

［9］張廣納，邵景安，王金亮，等.三峽庫區重慶段農村面源污染時空格局演變特征［J］.自然資源學報，2015，30（7）：1197-1209.

［10］賴斯蕓.非點源污染調查評估方法及其應用研究［D］.北京：清華大學，2004.

［11］陶園，王少麗，管孝艷，等.青海省農業面源污染源特征分析［J］.農業工程學報，2019，35（10）：164-172.

［12］馬國霞，於方，曹東，等.中國農業面源污染物排放量計算及中長期預測［J］.環境科學學報，2012，32（2）：489-497.

［13］《非常規飼料資源的開發與利用》研究組.非常規飼料資源的開發與利用［M］.北京：中國農業出版社，1996：14-18.

［14］中國農業科學院，環境保護部南京環境科學研究所.第一次全國污染源普查畜禽養殖業源產排污系數手冊［R］.2009.

［15］姜磊，柏玲，吳玉鳴.中國省域經濟、資源與環境協調分析：兼論三系統耦合公式及其擴展形式［J］.自然資源學報，2017，32（5）：788-799.

［16］陳宗峰，李裕瑞，劉彥隨.黃土丘陵溝壑區鄉村聚落分布格局特征與類型［J］.農業工程學報，2017，33（14）：266-274.

［17］王曼曼，吳秀芹，吳斌，等.鹽池北部風沙區鄉村聚落空間格局演變分析［J］.農業工程學報，2016，32（8）：260-271.

［18］逯進，周惠民.中國省域人力資本與經濟增長耦合關系的實證分析［J］.數量經濟技術經濟研究，2013，30（9）：3-19，36.