孟祥海,周海川,周海文
(1.淮陰師范學(xué)院經(jīng)濟(jì)與管理學(xué)院,江蘇 淮安 223001;2.中國(guó)林業(yè)科學(xué)院林業(yè)科技信息研究所,北京 100091;3.中國(guó)人民大學(xué)農(nóng)業(yè)與農(nóng)村發(fā)展學(xué)院,北京 100872)
2014年以來(lái),政府密集出臺(tái)了《畜禽規(guī)模養(yǎng)殖污染防治條例》等一系列養(yǎng)殖環(huán)保政策法規(guī),各地紛紛重新劃定“禁養(yǎng)區(qū)”,加大了環(huán)保執(zhí)法力度,畜禽規(guī)模養(yǎng)殖面臨巨大的環(huán)保壓力,同時(shí)新政策更加突出畜禽糞便的資源化利用,以種養(yǎng)結(jié)合、循環(huán)利用為主體,解決糞便還田利用的瓶頸問(wèn)題。江蘇是畜禽養(yǎng)殖大省,2015年全省生豬出欄2 978.32萬(wàn)頭,規(guī)?;捷^高,全省年出欄500頭以上的大中型規(guī)模養(yǎng)殖比例達(dá)64%,居全國(guó)第6,而全省年出欄生豬50頭以下散戶(hù)為47.8萬(wàn)戶(hù),生豬出欄量?jī)H占全省出欄總量的12%,低于全國(guó)平均水平,域內(nèi)河網(wǎng)密布,畜禽養(yǎng)殖污染防治形勢(shì)嚴(yán)峻。2012年9月,江蘇省農(nóng)委與省環(huán)保廳聯(lián)合下發(fā)《關(guān)于進(jìn)一步加強(qiáng)農(nóng)業(yè)源污染減排工作的意見(jiàn)》,將化學(xué)需氧量排放量占到農(nóng)業(yè)源污染96%的畜禽養(yǎng)殖業(yè)污染治理作為農(nóng)業(yè)源減排的“一號(hào)工程”。2016年8月,江蘇省農(nóng)委出臺(tái)《江蘇調(diào)整優(yōu)化畜禽養(yǎng)殖布局推進(jìn)養(yǎng)殖污染防治工作指導(dǎo)意見(jiàn)》,規(guī)劃截至2017年底,完成全省畜禽養(yǎng)殖禁養(yǎng)區(qū)劃定和禁養(yǎng)區(qū)內(nèi)畜禽養(yǎng)殖場(chǎng)(小區(qū))和養(yǎng)殖專(zhuān)業(yè)戶(hù)關(guān)閉或搬遷工作,到2020年,全省穩(wěn)定生豬出欄量在3 000萬(wàn)頭左右,年出欄500頭以上生豬規(guī)模養(yǎng)殖比例達(dá)80%,家禽出欄7.5億只,羊出欄800萬(wàn)頭,奶牛存欄21萬(wàn)頭;規(guī)模場(chǎng)畜禽糞便處理設(shè)施基本全配套,糞便綜合利用率達(dá)95%,形成與資源環(huán)境承載能力相匹配、與種植業(yè)布局相協(xié)調(diào)的生產(chǎn)格局,實(shí)現(xiàn)畜禽養(yǎng)殖污染基本治理和畜禽生產(chǎn)基本穩(wěn)定的雙重目標(biāo)。
畜禽糞便具有資源屬性,種養(yǎng)結(jié)合是實(shí)現(xiàn)畜禽糞便資源化利用和治理養(yǎng)殖污染的重要途徑。在國(guó)家、省域尺度,我國(guó)種養(yǎng)總量失衡與省域間結(jié)構(gòu)差異并存。2009年,我國(guó)單位耕地面積的畜禽氮污染平均負(fù)荷已達(dá)138.13 kg·hm-2,其中四川等6省市高達(dá)202.98 kg·hm-2以上,遠(yuǎn)超出歐盟農(nóng)業(yè)政策規(guī)定的糞肥年施氮量170 kg·hm-2的警戒線(xiàn),其余大部分省區(qū)已接近該警戒值[1];2010年,我國(guó)畜禽養(yǎng)殖環(huán)境總?cè)萘繛?29.56億頭豬當(dāng)量,實(shí)際養(yǎng)殖總量約占環(huán)境容量的1/4,但考慮化肥的使用情況,約有20個(gè)省超過(guò)本省50%環(huán)境容量[2]。在縣域尺度,閻波杰等[3]利用空間自相關(guān)分析方法對(duì)安徽省縣域耕地畜禽養(yǎng)殖廢棄物氮負(fù)荷量的空間格局和空間演變進(jìn)行研究,結(jié)果表明,2004—2012年間,安徽省各縣市耕地畜禽養(yǎng)殖廢棄物氮負(fù)荷具有顯著的空間正相關(guān)性和潛在的空間依賴(lài)性。在農(nóng)場(chǎng)或園區(qū)尺度,研究認(rèn)為種養(yǎng)結(jié)合園區(qū)能夠較好地實(shí)現(xiàn)養(yǎng)殖廢棄物循環(huán)利用和系統(tǒng)能量投入產(chǎn)出平衡,但合理的經(jīng)濟(jì)激勵(lì)政策、科學(xué)的養(yǎng)分管理措施、足夠的農(nóng)田消納面積和適當(dāng)?shù)酿B(yǎng)殖密度是實(shí)現(xiàn)種養(yǎng)結(jié)合良性循環(huán)的關(guān)鍵[4-6]。歐美發(fā)達(dá)國(guó)家畜禽養(yǎng)殖業(yè)規(guī)模化發(fā)展較早,他們大多采取集約化、農(nóng)場(chǎng)化的種養(yǎng)結(jié)合模式,養(yǎng)殖場(chǎng)通過(guò)配套足夠的土地或與周邊種植業(yè)農(nóng)場(chǎng)合作以消納畜禽糞便,并且根據(jù)土地面積嚴(yán)格控制養(yǎng)殖規(guī)模[7],這種農(nóng)場(chǎng)間的種養(yǎng)合作被認(rèn)為是經(jīng)濟(jì)和效率優(yōu)化的結(jié)果,可以實(shí)現(xiàn)畜禽養(yǎng)殖和農(nóng)作物生產(chǎn)的畜禽糞便交易成本最小化以及肥料利用效率最大化[8]。Logistic回歸分析模型被廣泛運(yùn)用于畜禽養(yǎng)殖污染治理愿意的影響因素研究。張暉等[9]基于長(zhǎng)三角地區(qū)207戶(hù)生豬養(yǎng)殖戶(hù)的實(shí)地調(diào)查數(shù)據(jù),運(yùn)用多元Logisitic回歸模型分析了影響農(nóng)戶(hù)參與畜禽糞便無(wú)害化處理意愿的影響因素;王鵬等[10]基于福建南平地區(qū)286戶(hù)農(nóng)戶(hù)的問(wèn)卷調(diào)查,運(yùn)用多元Logisitic回歸模型對(duì)農(nóng)戶(hù)參與畜禽養(yǎng)殖污染整治的意愿及其影響因素進(jìn)行實(shí)證研究;張郁等[11]通過(guò)建立多元有序Logistic模型,研究了環(huán)境風(fēng)險(xiǎn)感知對(duì)湖北省280戶(hù)規(guī)模生豬養(yǎng)殖戶(hù)環(huán)境行為的影響。文獻(xiàn)回顧發(fā)現(xiàn),在政策層面推動(dòng)區(qū)域種養(yǎng)結(jié)合,實(shí)現(xiàn)畜禽養(yǎng)殖生產(chǎn)與生態(tài)協(xié)調(diào)發(fā)展,既要考察區(qū)域種養(yǎng)平衡狀況,為種植與養(yǎng)殖產(chǎn)業(yè)空間規(guī)劃布局做好支撐,又要量化分析養(yǎng)殖業(yè)主實(shí)施種養(yǎng)結(jié)合意愿的影響因素,為制定激勵(lì)政策提供科學(xué)依據(jù)。當(dāng)前,在種養(yǎng)平衡和畜禽養(yǎng)殖場(chǎng)污染治理意愿方面的研究已十分豐富,但綜合運(yùn)用這2種方法更有助于系統(tǒng)研究特定區(qū)域種養(yǎng)結(jié)合問(wèn)題。為此,筆者試圖估算江蘇全省及各地市種養(yǎng)平衡空間狀況,并考察畜禽養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合意愿的影響因素,以期為江蘇省種養(yǎng)產(chǎn)業(yè)規(guī)劃布局和畜禽養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合政策出臺(tái)提供參考。
(1)數(shù)據(jù)來(lái)源。畜禽養(yǎng)殖存出欄數(shù)據(jù)主要來(lái)源于2001—2013年《江蘇統(tǒng)計(jì)年鑒》和《江蘇農(nóng)村統(tǒng)計(jì)年鑒》,部分來(lái)源于《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》和《中國(guó)畜牧業(yè)年鑒》。參照文獻(xiàn)[9,12]和畜禽養(yǎng)殖生產(chǎn)實(shí)踐,生豬的飼養(yǎng)周期一般為199 d,以年出欄量作為當(dāng)年飼養(yǎng)總量;牛按用途分為肉牛、役用牛和奶牛,肉牛一般1 a內(nèi)出欄,役用牛和奶牛飼養(yǎng)周期一般大于1 a,故肉牛采用年出欄量作為當(dāng)年飼養(yǎng)總量,役用牛和奶牛采用年末存欄量作為當(dāng)年飼養(yǎng)總量;羊、馬、驢和騾的飼養(yǎng)周期一般大于1 a,故采用年末存欄量作為當(dāng)年的飼養(yǎng)總量;家禽的飼養(yǎng)周期一般小于1 a,故采用年出欄量作為當(dāng)年的飼養(yǎng)總量。
(2)估算方法。農(nóng)田生產(chǎn)系統(tǒng)的氮磷平衡狀況是決定土壤肥力、農(nóng)作物產(chǎn)量和農(nóng)業(yè)環(huán)境影響的重要因素,系統(tǒng)氮磷缺損會(huì)降低土壤肥力,造成作物減產(chǎn);系統(tǒng)氮磷盈余會(huì)增加土壤氮磷的流失風(fēng)險(xiǎn),威脅地表水和地下水環(huán)境[13]。鑒于農(nóng)田氮磷流失監(jiān)測(cè)困難,隨機(jī)性強(qiáng),運(yùn)用土壤表觀(guān)養(yǎng)分平衡理論,核算農(nóng)田生產(chǎn)系統(tǒng)氮磷等養(yǎng)分盈余或缺損狀況,是分析農(nóng)業(yè)環(huán)境影響的重要方法。土壤表觀(guān)養(yǎng)分平衡理論的核心是農(nóng)田系統(tǒng)養(yǎng)分守恒,在農(nóng)田系統(tǒng)中化肥和畜禽糞便是養(yǎng)分輸入的主要來(lái)源,作物移走和養(yǎng)分損失則是養(yǎng)分輸出的主要途徑,良性的農(nóng)田系統(tǒng)養(yǎng)分輸入量與輸出量應(yīng)相等,即:化肥養(yǎng)分輸入量+畜禽糞便養(yǎng)分輸入量=農(nóng)作物移走養(yǎng)分量(假定農(nóng)作物秸稈全部還田)[14]。為評(píng)估特定區(qū)域畜禽養(yǎng)殖污染水平,假定畜禽糞便不跨區(qū)域流動(dòng)。計(jì)算公式為
M=F/S,
(1)
f=YN/P/(Ym+YN/P),
(2)
F=Ym+YN/P-Yc,
(3)
(4)
(5)
式(1)~(5)中,M為區(qū)域農(nóng)田土壤氮(磷)盈余強(qiáng)度,kg·hm-2;F為農(nóng)田土壤氮(磷)盈余量,萬(wàn)t;S為農(nóng)作物播種面積,hm2;f為化肥氮(磷)投入占農(nóng)田氮(磷)投入比例,%;Ym為畜禽糞便含氮(磷)量,萬(wàn)t;YN/P為化肥中氮(磷)折純量,萬(wàn)t;Qi為第i種畜禽的存欄或出欄量,萬(wàn)頭或萬(wàn)只;ri為第i種畜禽的糞便排泄系數(shù),t·a-1;pi為第i種畜禽的糞便的氮(磷)含量,%;Cj為第j類(lèi)農(nóng)作物年產(chǎn)量,萬(wàn)t;θj為第j類(lèi)農(nóng)作物100 kg產(chǎn)量所需的氮(磷)比例,%;Yc為農(nóng)作物移走的氮(磷)量,萬(wàn)t。采用M值估算區(qū)域農(nóng)田土壤氮(磷)環(huán)境風(fēng)險(xiǎn),f表示化肥投入對(duì)農(nóng)田氮(磷)盈余的貢獻(xiàn)。
畜禽糞便和農(nóng)作物氮磷含量參照文獻(xiàn)[12,15-16]獲得(表1~2)。
表1 一個(gè)飼養(yǎng)周期內(nèi)單位畜禽的糞便排泄系數(shù)及糞便氮磷含量[12,15]
Table 1 Manure excretion coefficient of an animal per breeding cycle and nitrogen and phosphorus contents in the manure
畜禽種類(lèi)糞便排泄系數(shù)/(t·a-1)w(總氮)/%w(總磷)/%豬1.05470.2380.074役用牛10.10.3510.082肉牛7.70.3510.082奶牛19.40.3510.082馬5.90.3780.077驢/騾5.00.3780.077羊0.871.0140.220家禽0.0321.2500.940
該研究中養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合指畜禽養(yǎng)殖場(chǎng)將養(yǎng)殖與種植相結(jié)合,利用畜禽糞便就近還田用于農(nóng)作物種植。根據(jù)已有研究,畜禽養(yǎng)殖場(chǎng)實(shí)施種養(yǎng)結(jié)合受多方面因素的影響,主要有:
(1)養(yǎng)殖場(chǎng)業(yè)主對(duì)養(yǎng)殖污染的認(rèn)知。與傳統(tǒng)農(nóng)戶(hù)散養(yǎng)不同,規(guī)?;B(yǎng)殖場(chǎng)畜禽糞便的隨意排放對(duì)環(huán)境造成嚴(yán)重污染,養(yǎng)殖業(yè)主對(duì)養(yǎng)殖場(chǎng)畜禽糞便排放對(duì)周邊環(huán)境危害的認(rèn)知程度越高,通過(guò)種養(yǎng)結(jié)合治理糞污的意愿越強(qiáng)[17]。
表2 農(nóng)作物氮磷含量[16]
Table 2 N and P contents in crop
農(nóng)作物 w(氮)/%w(磷)/%谷物2.741.18水果0.510.18蔬菜0.400.19豆類(lèi)5.151.33薯類(lèi)0.430.19棉花5.001.80花生6.801.30油菜籽5.802.50芝麻8.232.07甘蔗0.380.04甜菜0.400.15煙葉4.100.70茶葉6.402.00
(2)政策引導(dǎo)作用的影響?,F(xiàn)行的畜禽養(yǎng)殖污染治理政策主要分為命令控制型和經(jīng)濟(jì)激勵(lì)型2類(lèi)。環(huán)保部門(mén)對(duì)養(yǎng)殖場(chǎng)環(huán)境監(jiān)管越嚴(yán)格,養(yǎng)殖場(chǎng)通過(guò)種養(yǎng)結(jié)合治理糞污的可能性越高;業(yè)主對(duì)畜禽糞便資源化利用經(jīng)濟(jì)激勵(lì)政策了解越多,開(kāi)展種養(yǎng)結(jié)合治理糞污的意愿可能越強(qiáng);已獲得過(guò)相關(guān)畜禽糞便資源化利用政策補(bǔ)貼的養(yǎng)殖場(chǎng),實(shí)施種養(yǎng)結(jié)合治理糞污的意愿可能越強(qiáng)[18]。
(3)社會(huì)網(wǎng)絡(luò)行為影響。養(yǎng)殖場(chǎng)業(yè)主處在特定社會(huì)網(wǎng)絡(luò)之中,業(yè)主的種養(yǎng)結(jié)合意愿受到周邊社會(huì)網(wǎng)絡(luò)行為的影響[19]。采用種養(yǎng)結(jié)合的同行越多,可能意味著實(shí)施種養(yǎng)結(jié)合的經(jīng)營(yíng)風(fēng)險(xiǎn)越小,業(yè)主實(shí)施種養(yǎng)結(jié)合的意愿可能越高。養(yǎng)殖場(chǎng)因環(huán)保問(wèn)題影響到與周邊村民、村委會(huì)或政府的關(guān)系,其來(lái)自社會(huì)的環(huán)保壓力越大,業(yè)主實(shí)施種養(yǎng)結(jié)合的意愿可能越高。
(4)養(yǎng)殖業(yè)主個(gè)體特征。養(yǎng)殖業(yè)主行為決策與其個(gè)體特征密切相關(guān)[20],受教育程度越高的業(yè)主環(huán)保意識(shí)越高,實(shí)施種養(yǎng)結(jié)合的意愿可能越高;養(yǎng)殖場(chǎng)養(yǎng)殖規(guī)模越大,畜禽糞便排放量越大,環(huán)保壓力也就越大,實(shí)施種養(yǎng)結(jié)合開(kāi)展糞污治理的意愿可能越高;養(yǎng)殖場(chǎng)通過(guò)流轉(zhuǎn)或與周邊農(nóng)民合作的方式,擁有的配套土地越多,通過(guò)種養(yǎng)結(jié)合開(kāi)展糞污治理的意愿可能越高;養(yǎng)殖場(chǎng)經(jīng)營(yíng)效益越好,越有能力開(kāi)展環(huán)保治理,實(shí)施種養(yǎng)結(jié)合開(kāi)展糞污治理的意愿可能越高;同時(shí),種養(yǎng)結(jié)合不僅是消納養(yǎng)殖糞污的一種手段,也是一種經(jīng)營(yíng)模式,需要養(yǎng)殖業(yè)主具備生態(tài)農(nóng)業(yè)經(jīng)營(yíng)能力,養(yǎng)殖業(yè)主對(duì)種養(yǎng)結(jié)合經(jīng)營(yíng)難度的認(rèn)知,也會(huì)影響其實(shí)施種養(yǎng)結(jié)合的意愿。
模型數(shù)據(jù)來(lái)源于課題組對(duì)畜禽養(yǎng)殖密集的環(huán)洪澤湖流域淮安(淮陰區(qū)、淮安區(qū)、清浦區(qū))和宿遷(泗陽(yáng)縣)184家規(guī)?;B(yǎng)豬場(chǎng)的調(diào)研。環(huán)洪澤湖地區(qū)是江蘇畜禽養(yǎng)殖密集區(qū),2013年淮安和宿遷2市生豬、肉牛和家禽出欄量分別占全省總量的23.68%、26.92%和24.34%,且緊鄰湖區(qū)的洪澤縣等6縣(區(qū))生豬、肉牛和家禽出欄量分別占全省的12.92%、22.86%和11.98%。中國(guó)科學(xué)院南京地理與湖泊研究所于2011年12月發(fā)布的《中國(guó)五大淡水湖泊保護(hù)與發(fā)展報(bào)告》指出,洪澤湖每年入湖污染物總量已經(jīng)大大超過(guò)了湖泊的自?xún)裟芰?且對(duì)洪澤湖周邊產(chǎn)業(yè)調(diào)查發(fā)現(xiàn),湖區(qū)萬(wàn)元產(chǎn)值污染水平以畜禽養(yǎng)殖業(yè)的排污水平最高。從一定程度上而言,環(huán)洪澤湖地區(qū)畜禽養(yǎng)殖業(yè)發(fā)展所面臨的污染問(wèn)題正是現(xiàn)階段江蘇畜禽養(yǎng)殖業(yè)發(fā)展面臨環(huán)境約束的縮影。模型相關(guān)變量含義及描述性統(tǒng)計(jì)見(jiàn)表3。
表3 畜禽養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合意愿影響因素變量及描述性統(tǒng)計(jì)
Table 3 Variables affecting willingness of livestock farms to combine breeding with planting and descriptive statistics
類(lèi)別變量變量含義變量賦值均值標(biāo)準(zhǔn)差預(yù)期方向因變量 Y實(shí)施種養(yǎng)結(jié)合的意愿0=不愿意采納;1=愿意采納0.4460.498養(yǎng)殖業(yè)主環(huán)境認(rèn)知X1認(rèn)為養(yǎng)殖場(chǎng)對(duì)周邊環(huán)境的污染程度0=沒(méi)有污染;1=污染很小;2=一般污染;3=污染很大;4=污染嚴(yán)重1.3100.995正向政府引導(dǎo)行為X2對(duì)畜禽養(yǎng)殖糞便資源化利用政策的了解程度1=不了解;2=了解一點(diǎn);3=一般了解;4=比較了解;5=熟悉3.3861.007正向X3是否獲得過(guò)養(yǎng)殖環(huán)保政策補(bǔ)貼0=否;1=是0.3860.488正向X4環(huán)保部門(mén)的監(jiān)管壓力0=沒(méi)有壓力;1=壓力很小;2=有一定壓力;3=壓力較大;4=壓力很大1.7501.352正向社會(huì)行為X5是否得到同行的推薦采納種養(yǎng)結(jié)合模式0=否;1=是0.0870.283正向X6是否因環(huán)保問(wèn)題影響到與周邊農(nóng)民、村委會(huì)或政府的關(guān)系0=否;1=是0.6470.479正向養(yǎng)殖場(chǎng)業(yè)主個(gè)體特征X7業(yè)主文化程度1=文盲;2=小學(xué);3=初中;4=高中或中專(zhuān);5=大專(zhuān)及以上3.6361.299正向X8年存欄規(guī)模100~299、300~499、500~999、1000~2999、3000~4999、5000~9999、1萬(wàn)頭及以上分別賦值1、2、3、4、5、6、73.5712.050不確定X9配套土地規(guī)?!?.67、>6.67~33.35、>33.35~66.7、>66.7~200.1、>200.1hm2分別賦值1、2、3、4、51.8591.229正向X10經(jīng)濟(jì)效益情況1=很不好;2=不太好;3=一般;4=比較好;5=很好3.2450.775正向X11種養(yǎng)結(jié)合模式經(jīng)營(yíng)難度是否較大0=否;1=是0.4350.497負(fù)向
從農(nóng)田氮磷盈余強(qiáng)度的變化趨勢(shì)(圖1)來(lái)看,2001—2013年江蘇全省及各地市氮磷盈余強(qiáng)度總體上呈現(xiàn)出上升—下降的2階段特征,且相對(duì)于2001年,2013年全省及多數(shù)地市氮磷盈余強(qiáng)度下降幅度較大,全省農(nóng)田氮磷平衡狀況正在改觀(guān)。在農(nóng)田氮盈余強(qiáng)度方面,全省和南京、無(wú)錫、常州、連云港、淮安、鹽城、宿遷7地市在2003年達(dá)到峰值后呈現(xiàn)下降趨勢(shì),蘇州和揚(yáng)州2地市在2004年達(dá)到峰值后呈現(xiàn)下降趨勢(shì),鎮(zhèn)江、南通和徐州在2007年達(dá)到峰值后呈現(xiàn)下降趨勢(shì),泰州農(nóng)田氮盈余強(qiáng)度呈現(xiàn)持續(xù)下降趨勢(shì),徐州和宿遷氮盈余強(qiáng)度雖然呈現(xiàn)下降趨勢(shì),但2013年的氮盈余強(qiáng)度明顯高于2001年。農(nóng)田磷盈余強(qiáng)度與氮盈余強(qiáng)度變化趨勢(shì)基本一致,全省以及南京、無(wú)錫、連云港、淮安、鹽城及宿遷6地市在2003年達(dá)到峰值后呈現(xiàn)下降趨勢(shì),蘇州和揚(yáng)州在2004年達(dá)到峰值后呈現(xiàn)下降趨勢(shì),鎮(zhèn)江和南通在2007年達(dá)到峰值后呈現(xiàn)下降趨勢(shì),泰州磷盈余強(qiáng)度呈現(xiàn)持續(xù)下降趨勢(shì),而常州和徐州磷盈余強(qiáng)度總體上呈現(xiàn)上升趨勢(shì),但徐州磷盈余強(qiáng)度在2011年后有所回落。
從農(nóng)田氮磷盈余強(qiáng)度水平(圖1)來(lái)看,各地市普遍存在農(nóng)田氮磷盈余現(xiàn)象,蘇北各地市氮磷盈余強(qiáng)度普遍高于蘇南、蘇中地區(qū)。2013年,全省農(nóng)田氮盈余強(qiáng)度為113 kg·hm-2,其中4地市高于全省總體水平,從高到低依次為徐州(253.91 kg·hm-2)、宿遷(180.78 kg·hm-2)、連云港(120.29 kg·hm-2)和淮安(117.77 kg·hm-2);9地市低于全省總體水平,從高到低依次為鹽城(83.12 kg·hm-2)、常州(75.11 kg·hm-2)、南通(70.00 kg·hm-2)、揚(yáng)州(52.90 kg·hm-2)、蘇州(49.96 kg·hm-2)、無(wú)錫(48.98 kg·hm-2)、南京(24.42 kg·hm-2)、泰州(4.88 kg·hm-2)和鎮(zhèn)江(-8.26 kg·hm-2),鎮(zhèn)江自2012年已出現(xiàn)農(nóng)田氮虧損,泰州已接近農(nóng)田氮平衡。2013年,全省農(nóng)田磷盈余強(qiáng)度為23.13 kg·hm-2,其中3地市高于全省總體水平,從高到低依次為徐州(80.78 kg·hm-2)、宿遷(43.63 kg·hm-2)和常州(35.21 kg·hm-2);10地市低于全省總體水平,從高到低依次為淮安(20.01 kg·hm-2)、南通(16.87 kg·hm-2)、鹽城(16.27 kg·hm-2)、連云港(12.69 kg·hm-2)、無(wú)錫(4.88 kg·hm-2)、蘇州(2.97 kg·hm-2)、揚(yáng)州(-0.81 kg·hm-2)、南京(-1.85 kg·hm-2)、鎮(zhèn)江(-15.38 kg·hm-2)和泰州(-18.55 kg·hm-2),南京、鎮(zhèn)江、揚(yáng)州和泰州分別自2012、2009、2011和2007年起出現(xiàn)農(nóng)田磷虧損。
圖1 2001—2013年江蘇省及各地市農(nóng)田氮、磷盈余強(qiáng)度以及化肥氮、磷投入比例 Fig.1 Farmland N,P surplus intensity and proportion of N,P fertilizer input of Jiangsu Province and local cities in 2001-2013
總體而言,全省各地市農(nóng)田氮盈余強(qiáng)度顯著高于磷盈余強(qiáng)度,考慮到未將人糞尿、生物固氮、復(fù)合肥投入、作物秸稈還田和農(nóng)田氮磷流失等因素納入估算框架,可推測(cè)全省僅鎮(zhèn)江和泰州2地市農(nóng)田氮磷已基本實(shí)現(xiàn)平衡。
從農(nóng)田氮磷的投入結(jié)構(gòu)來(lái)看,化肥是江蘇全省及各地市農(nóng)田氮磷投入的主要來(lái)源。基于農(nóng)田化肥氮投入比例的觀(guān)察(圖1),2013年全省化肥氮投入占農(nóng)田投入的79%,其中連云港、揚(yáng)州、淮安、泰州和宿遷的化肥氮投入比例高于全省水平,連云港化肥氮投入占農(nóng)田投入的89%,為全省最高;無(wú)錫、蘇州、鎮(zhèn)江、鹽城、徐州、南京、南通和常州低于全省水平,常州化肥氮投入占農(nóng)田投入的58%,為全省最低;13 a來(lái),除揚(yáng)州、連云港、淮安和鹽城外,全省及各地市化肥氮投入比例均有所下降。基于農(nóng)田化肥磷投入比例的觀(guān)察(圖1),2013年全省化肥磷投入占農(nóng)田投入的63%,其中連云港、淮安、泰州、揚(yáng)州和宿遷5地市的化肥磷投入比例高于全省水平,連云港化肥磷投入占農(nóng)田投入的80%,為全省最高;蘇州、無(wú)錫、鎮(zhèn)江、鹽城、徐州、南京、南通和常州低于全省水平,常州化肥磷投入占農(nóng)田投入的35%,為全省最低;13 a來(lái),除揚(yáng)州和淮安外,全省及其他各地市化肥磷投入比例均有所下降??傮w而言,化肥的過(guò)量投入是江蘇全省及各地市農(nóng)田氮磷盈余的主要驅(qū)動(dòng)因素,迫切需要發(fā)展種養(yǎng)結(jié)合,促進(jìn)畜禽糞肥還田,減少化肥投入。
運(yùn)用SPSS 19.0統(tǒng)計(jì)軟件對(duì)184份有效樣本數(shù)據(jù)進(jìn)行多元Logistic回歸分析。采用Wald值向后逐步回歸方法,逐步剔除Wald值最小的變量,直到所有的解釋變量都達(dá)到顯著水平為止。在該逐步回歸模型中共進(jìn)行了6次回歸分析,選擇第1步和第6步的回歸分析結(jié)果進(jìn)行對(duì)比分析,根據(jù)模型最終估計(jì)結(jié)果(表4中第6步),模型整體擬合優(yōu)度較高,具有較強(qiáng)的可信度。
表4 畜禽養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合意愿影響因素估計(jì)結(jié)果
Table 4 Reckoning of the factors affecting willingness of livestock farms to combine breeding with planting
模型解釋變量第1步第6步系數(shù)標(biāo)準(zhǔn)差Wald值顯著度系數(shù)標(biāo)準(zhǔn)差Wald值顯著度X1-0.8570.4972.9750.085*-0.9530.3686.6980.010***X2-0.5460.6000.8260.363X32.0800.9005.3360.021**2.3030.58015.7930.000***X40.8230.5911.9400.1640.5360.2713.9020.048**X50.4301.4100.0930.761X60.5190.7890.4330.511X7-0.5850.4741.5240.217-0.5800.2684.6900.030**X80.5390.3971.8490.1740.7660.2449.9030.002***X91.6370.49311.0270.001***1.3950.30021.6110.000***X100.6230.6001.0810.299X110.4221.0630.1580.691常量-4.1863.0211.9200.166-2.0921.3462.4150.120
*、**、***分別表示統(tǒng)計(jì)值在10%、5%和1%置信水平上顯著。第1步和第6步的-2對(duì)數(shù)似然值、Cox&SnellR2、NagelkerkeR2和預(yù)測(cè)準(zhǔn)確率分別為53.654、0.560、0.748、87.5%和58.574、0.537、0.717、88.5%。
回歸結(jié)果顯示,養(yǎng)殖業(yè)主對(duì)養(yǎng)殖場(chǎng)對(duì)周邊環(huán)境污染程度的認(rèn)知、是否獲得過(guò)養(yǎng)殖環(huán)保政策補(bǔ)貼、環(huán)保部門(mén)的監(jiān)管壓力、業(yè)主文化程度、養(yǎng)殖規(guī)模(年存欄規(guī)模)和配套土地面積6個(gè)因素對(duì)養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合意愿具有顯著影響,而養(yǎng)殖場(chǎng)業(yè)主對(duì)畜禽養(yǎng)殖糞便資源化利用政策的了解程度、是否得到同行的推薦、是否因環(huán)保問(wèn)題影響到與周邊農(nóng)民、村委會(huì)或政府的關(guān)系以及養(yǎng)殖場(chǎng)經(jīng)濟(jì)效益情況這4個(gè)因素不具有顯著影響。具體分析如下:
(1)養(yǎng)殖業(yè)主對(duì)周邊環(huán)境污染程度的認(rèn)知對(duì)其種養(yǎng)結(jié)合的意愿具有顯著負(fù)影響,與前文預(yù)期假設(shè)不一致。模型中變量X1的系數(shù)為負(fù),即養(yǎng)殖業(yè)主認(rèn)為養(yǎng)殖場(chǎng)對(duì)周邊環(huán)境影響越小,開(kāi)展種養(yǎng)結(jié)合的意愿越高,與預(yù)期假設(shè)不一致,原因可能來(lái)自2個(gè)方面:一是因養(yǎng)殖污染問(wèn)題十分敏感,養(yǎng)殖業(yè)主在填寫(xiě)調(diào)查問(wèn)卷時(shí)有意回避該問(wèn)題,刻意淡化養(yǎng)殖場(chǎng)對(duì)周邊環(huán)境污染的程度;二是種養(yǎng)結(jié)合能實(shí)現(xiàn)畜禽養(yǎng)殖糞便的資源化利用,經(jīng)營(yíng)得當(dāng)能取得良好的經(jīng)濟(jì)效益,養(yǎng)殖場(chǎng)即便能夠較好地處理好養(yǎng)殖糞污,未對(duì)周邊環(huán)境未造成污染,仍愿意開(kāi)展種養(yǎng)結(jié)合。
(2)是否獲得過(guò)養(yǎng)殖環(huán)保財(cái)政補(bǔ)貼對(duì)養(yǎng)殖業(yè)主的種養(yǎng)結(jié)合意愿具有顯著正影響,與前文預(yù)期假設(shè)一致。模型中變量X3的系數(shù)為正,養(yǎng)殖業(yè)主獲得過(guò)養(yǎng)殖環(huán)保財(cái)政補(bǔ)貼,開(kāi)展種養(yǎng)結(jié)合的意愿越高,與預(yù)期假設(shè)一致。種養(yǎng)結(jié)合是治理養(yǎng)殖污染的一種手段,在環(huán)保領(lǐng)域享受過(guò)財(cái)政補(bǔ)貼的養(yǎng)殖業(yè)主,相對(duì)于未享受過(guò)補(bǔ)貼的業(yè)主,其對(duì)政策導(dǎo)向的認(rèn)知更清晰,會(huì)更愿意相信實(shí)施種養(yǎng)結(jié)合會(huì)得到政策的支持,由此開(kāi)展種養(yǎng)結(jié)合的意愿就越高。
(3)是否有來(lái)自環(huán)保部門(mén)的監(jiān)管壓力對(duì)養(yǎng)殖業(yè)主的種養(yǎng)結(jié)合意愿具有顯著正影響,與前文預(yù)期假設(shè)一致。模型中變量X4的系數(shù)為正,養(yǎng)殖業(yè)主面臨的環(huán)境監(jiān)管壓力越大,采納種養(yǎng)結(jié)合模式的意愿越高,與預(yù)期假設(shè)一致。新《環(huán)境保護(hù)法》、《畜禽規(guī)模養(yǎng)殖污染防治條例》等法律法規(guī)的實(shí)施,尤其是禁養(yǎng)區(qū)的劃分,使養(yǎng)殖場(chǎng)面臨前所未有的環(huán)保壓力,養(yǎng)殖場(chǎng)環(huán)保不達(dá)標(biāo)將面臨高昂的處罰,環(huán)保部門(mén)的監(jiān)管壓力倒逼養(yǎng)殖場(chǎng)開(kāi)展種養(yǎng)結(jié)合與污染整治。
(4)業(yè)主文化程度對(duì)其種養(yǎng)結(jié)合意愿具有顯著負(fù)影響,與前文預(yù)期假設(shè)不一致。模型中變量X7的系數(shù)為負(fù),養(yǎng)殖業(yè)主文化程度越低,開(kāi)展種養(yǎng)結(jié)合的意愿越高,與預(yù)期假設(shè)不一致。已有研究中,養(yǎng)殖業(yè)主環(huán)保意識(shí)會(huì)隨著文化程度的提高而提高,與預(yù)期不一致的原因可能是文化程度較高的養(yǎng)殖業(yè)主多是年輕一代,對(duì)種植業(yè)缺乏經(jīng)驗(yàn),從專(zhuān)業(yè)化的角度理解,認(rèn)為自身經(jīng)營(yíng)好種養(yǎng)結(jié)合的難度較大,同時(shí)種養(yǎng)結(jié)合更是一種傳統(tǒng)理念,文化程度低的業(yè)主一般年齡較大,對(duì)種養(yǎng)結(jié)合的認(rèn)同感更強(qiáng)。
(5)養(yǎng)殖場(chǎng)年存欄生豬規(guī)模對(duì)其種養(yǎng)結(jié)合意愿具有顯著正影響,與前文預(yù)期假設(shè)一致。模型中變量X8的系數(shù)為正,隨著養(yǎng)殖規(guī)模的擴(kuò)大,養(yǎng)殖業(yè)主的種養(yǎng)結(jié)合意愿越高。養(yǎng)殖規(guī)模越大,畜禽糞便處理的環(huán)保壓力越大,養(yǎng)殖場(chǎng)更愿意開(kāi)展種養(yǎng)結(jié)合治理畜禽糞污。
(6)養(yǎng)殖場(chǎng)的農(nóng)地配套規(guī)模對(duì)其種養(yǎng)結(jié)合意愿具有顯著正影響,與前文預(yù)期假設(shè)一致。模型中變量X9的系數(shù)為正,養(yǎng)殖場(chǎng)配套的土地面積越多,養(yǎng)殖業(yè)主開(kāi)展種養(yǎng)結(jié)合的意愿越高。農(nóng)地是種養(yǎng)結(jié)合的必備條件,養(yǎng)殖場(chǎng)能夠流轉(zhuǎn)到足夠的農(nóng)地或周邊農(nóng)民愿意與其達(dá)成糞肥消納協(xié)議,養(yǎng)殖場(chǎng)才具備開(kāi)展種養(yǎng)結(jié)合的條件,實(shí)施種養(yǎng)結(jié)合的意愿就會(huì)越高。
基于對(duì)省域農(nóng)田氮磷平衡狀況的考察,2001—2013年間江蘇全省及各地市普遍存在農(nóng)田氮磷盈余現(xiàn)象,化肥的過(guò)量投入是造成農(nóng)田氮磷盈余主要驅(qū)動(dòng)因素,全省及各地市氮磷盈余強(qiáng)度總體上呈現(xiàn)出上升—下降的2階段特征,蘇北各地市氮磷盈余強(qiáng)度普遍高于蘇南、蘇中,但相對(duì)于2001年,2013年全省及多數(shù)地市氮磷盈余強(qiáng)度下降幅度較大,全省農(nóng)田氮磷平衡狀況正在改觀(guān)。運(yùn)用Logistic模型研究養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合意愿的影響因素,得出養(yǎng)殖場(chǎng)是否獲得過(guò)養(yǎng)殖環(huán)保政策補(bǔ)貼、環(huán)保部門(mén)的監(jiān)管壓力、養(yǎng)殖規(guī)模和配套土地面積等因素,對(duì)養(yǎng)殖場(chǎng)種養(yǎng)結(jié)合意愿具有顯著正向影響。
種養(yǎng)結(jié)合是促進(jìn)畜禽糞肥還田、防治畜禽養(yǎng)殖污染和實(shí)現(xiàn)化肥減量化的重要手段,在省域?qū)用嫱苿?dòng)畜禽養(yǎng)殖業(yè)種養(yǎng)結(jié)合,應(yīng)在執(zhí)行中央政府出臺(tái)的相關(guān)法規(guī)政策的基礎(chǔ)上,做好農(nóng)田養(yǎng)分管理和種養(yǎng)平衡空間規(guī)劃,科學(xué)合理地確定各地畜禽養(yǎng)殖規(guī)模,實(shí)現(xiàn)畜禽糞便排放與農(nóng)業(yè)種植需肥在時(shí)間與空間上的匹配。加強(qiáng)畜禽養(yǎng)殖環(huán)境監(jiān)管,并對(duì)實(shí)施種養(yǎng)結(jié)合的畜禽養(yǎng)殖場(chǎng)給予扶持,鼓勵(lì)大中型養(yǎng)殖場(chǎng)流轉(zhuǎn)土地,開(kāi)展農(nóng)牧一體化經(jīng)營(yíng),引導(dǎo)小規(guī)模畜禽養(yǎng)殖場(chǎng)、養(yǎng)殖小區(qū)和養(yǎng)殖散戶(hù)采取多元化的集中治理模式。
[1] 楊飛,楊世琦,諸云強(qiáng),等.中國(guó)近30年畜禽養(yǎng)殖量及其耕地氮污染負(fù)荷分析[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(5):1-11.[YANG Fei,YANG Shi-qi,ZHU Yun-qiang,etal.Analysis on Livestock and Poultry Production and Nitrogen Pollution Load of Cultivated Land During Last 30 Years in China[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29(5):1-11.]
[2] 耿維,胡林,崔建宇,等.中國(guó)區(qū)域畜禽糞便能源潛力及總量控制研究[J].農(nóng)業(yè)工程學(xué)報(bào),2013,29(1):171-179;295.[GENG Wei,HU Lin,CUI Jian-yu,etal.Biogas Energy Potential for Livestock Manure and Gross Control of Animal Feeding in Region Level of China[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29(1):171-179.]
[3] 閻波杰,潘瑜春,閆靜杰.安徽省縣域耕地畜禽養(yǎng)殖廢棄物養(yǎng)分負(fù)荷時(shí)空演變特征[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào),2016,32(3):466-472.[YAN Bo-jie,PAN Yu-chun,YAN Jing-jie.Spatio-Temporal Evolution of Livestock Manure Nutrient Load in Farmland on a County Scale:A Case Study of Anhui Province[J].Journal of Ecology and Rural Environment,2016,32(3):466-472.]
[4] 馬永喜,王穎.規(guī)?;竽琉B(yǎng)殖廢棄物處理的環(huán)境經(jīng)濟(jì)優(yōu)化研究:基于生態(tài)經(jīng)濟(jì)模型的分析[J].農(nóng)業(yè)現(xiàn)代化研究,2014,35(3):340-344.[MA Yong-xi,WANG Ying.Optimization of Environmental and Economic Effects in Waste Treatment of Intensive Livestock Production:An Analysis Based on "Ecological-E-c-o-n-o-m-i-c" Model[J].Research of Agricultural Modernization,2014,35(3):340-344.]
[5] 黃治平,郝利,高尚賓,等.西柏店村級(jí)養(yǎng)殖種植園區(qū)氮素流和能量流分析[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2011,30(2):395-403.[HUANG Zhi-ping,HAO Li,GAO Shang-bin,etal.Material Flow Analysis on Nitrogen Cycle and Energy Flow Analysis of Planting and Breeding Park of Village Scale at Xibaidian Village[J].Journal of Agro-Environment Science,2011,30(2):395-403.]
[6] 侯勇,高志嶺,馬文奇,等.京郊典型集約化“農(nóng)田-畜牧”生產(chǎn)系統(tǒng)氮素流動(dòng)特征[J].生態(tài)學(xué)報(bào),2012,32(4):1027-1036.[HOU Yong,GAO Zhi-ling,MA Wen-qi,etal.Nitrogen Flows in Intensive "Crop-Livestock" Production Systems Typically for the Peri-Urban Area of Beijing[J].Acta Ecologica Sinica,2012,32(4):1027-1036.]
[7] ASAI M,LANGER V.Collaborative Partnerships Between Organic Farmers in Livestock-Intensive Areas of Denmark[J].Organic Agriculture,2014,4(1):63-77.
[8] ANDERSSON H,LARSéN K,LAGERKVIST C J,etal.Farm Cooperation to Improve Sustainability[J].Ambio,2005,34(4/5):383-387.
[9] 張暉,虞祎,胡浩.基于農(nóng)戶(hù)視角的畜牧業(yè)污染處理意愿研究:基于長(zhǎng)三角生豬養(yǎng)殖戶(hù)的調(diào)查[J].農(nóng)村經(jīng)濟(jì),2011(10):92-94.[ZHANG Hui,YU Yi,HU Hao.Study on the Willingness to Treat the Pollution of Animal Husbandry Based on the Farmer′s Perspective:Based on the Survey of Pigs in the Yangtze River Delta[J].Rural Economy,2011(10):92-94.]
[10]王鵬,于瀟.農(nóng)戶(hù)參與畜禽養(yǎng)殖污染整治意愿及其影響因素:基于福建南平地區(qū)286份調(diào)查問(wèn)卷[J].湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版),2013,14(6):44-49.[WANG Peng,YU Xiao.Desire and Factors Influencing Farmers Participating in Livestock and Poultry Breeding Pollution Treatment:Based on 286 Questionnaires of Nanping,Fujian[J].Journal of Hunan Agricultural University(Social Sciences),2013,14(6):44-49.]
[11]張郁,江易華.環(huán)境規(guī)制政策情境下環(huán)境風(fēng)險(xiǎn)感知對(duì)養(yǎng)豬戶(hù)環(huán)境行為影響:基于湖北省280戶(hù)規(guī)模養(yǎng)殖戶(hù)的調(diào)查[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì),2016(11):76-86.[ZHANG Yu,JIANG Yi-hua.Impact of Enviromnental Risk Awareness on Environmental Behavior of Pig Household in the Environmental reguation policy context:Based on the Survey of 280 Households in Hubei Province[J].Journal of Agrotechnical Economics,2016(11):76-86.]
[12]王方浩,馬文奇,竇爭(zhēng)霞,等.中國(guó)畜禽糞便產(chǎn)生量估算及環(huán)境效應(yīng)[J].中國(guó)環(huán)境科學(xué),2006,26(5):614-617.[WANG Fang-hao,MA Wen-qi,DOU Zheng-xia,etal.The Estimation of the Production Amount of Animal Manure and Its Environmental Effect in China[J].China Environmental Science,2006,26(5):614-617.]
[13]陳敏鵬,陳吉寧.中國(guó)區(qū)域土壤表觀(guān)氮磷平衡清單及政策建議[J].環(huán)境科學(xué),2007,28(6):1305-1310.[CHEN Min-peng,CHEN Ji-ning.Inventory of Regional Surface Nutrient Balance and Policy Recommendations in China[J].Environmental Science,2007,28(6):1305-1310.]
[14]胡浩,郭利京.農(nóng)區(qū)畜牧業(yè)發(fā)展的環(huán)境制約及評(píng)價(jià):基于江蘇省的實(shí)證分析[J].農(nóng)業(yè)技術(shù)經(jīng)濟(jì),2011(6):36-42.[HU Hao,GUO Li-jing.Environmental Restriction and Evaluation of Animal Husbandry Development in Agricultural Region:An Empirical Analysis Based on Jiangsu Province[J].Journal of Agrotechnical Economics,2011(6):36-42.]
[15]李飛,董鎖成.西部地區(qū)畜禽養(yǎng)殖污染負(fù)荷與資源化路徑研究[J].資源科學(xué),2011,33(11):2204-2211.[LI Fei,DONG Suo-cheng.Pollution From Livestock and Poultry and Its Resource Strategy in West China[J].Resources Science,2011,33(11):2204-2211.]
[16]孟祥海,周海川,張俊飚.中國(guó)畜禽污染時(shí)空特征分析與環(huán)境庫(kù)茲涅茨曲線(xiàn)驗(yàn)證[J].干旱區(qū)資源與環(huán)境,2015,29(11):104-108.[MENG Xiang-hai,ZHOU Hai-chuan,ZHANG Jun-biao.Spatial and Temporal Characteristics and EKC Verification for Livestock Pollution[J].Journal of Arid Land Resources and Environment,2015,29(11):104-108.]
[17]孔凡斌,張維平,潘丹.養(yǎng)殖戶(hù)畜禽糞便無(wú)害化處理意愿及影響因素研究:基于5省754戶(hù)生豬養(yǎng)殖戶(hù)的調(diào)查數(shù)據(jù)[J].農(nóng)林經(jīng)濟(jì)管理學(xué)報(bào),2016,15(4):454-463.[KONG Fan-bin,ZHANG Wei-ping,PAN Dan.Factors Influencing Farmers′ Willingness to Adopt Harmless Disposal of Livestock Manure:Based on 754 Pig Farmers in Five Provinces[J].Journal of Agro-Forestry Economics and Management,2016,15(4):454-463.]
[18]馮淑怡,羅小娟,張麗軍,等.養(yǎng)殖企業(yè)畜禽糞尿處理方式選擇、影響因素與適用政策工具分析:以太湖流域上游為例[J].華中農(nóng)業(yè)大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版),2013(1):12-18.[FENG Shu-yi,LUO Xiao-juan,ZHANG Li-jun,etal.Livestock Enterprises′ Choice of Waste Treatment Modes,Its Influencing Factors and Possible Policy Instruments:A Case Study in Upper Reaches of Taihu Lake Basin[J].Journal of Huazhong Agricultural University,2013(1):12-18.]
[19]ZHENG C H,BLUEMLING B,LIU Y,etal.Managing Manure From China′s Pigs and Poultry:The Influence of Ecological Rationality[J].Ambio,2014,43(5):661-672.
[20]劉雪芬,楊志海,王雅鵬.畜禽養(yǎng)殖戶(hù)生態(tài)認(rèn)知及行為決策研究:基于山東、安徽等6省養(yǎng)殖戶(hù)的實(shí)地調(diào)研[J].中國(guó)人口·資源與環(huán)境,2013,23(10):169-176.[LIU Xue-fen,YANG Zhi-hai,WANG Ya-peng.Ecological Cognition and Behavior Decision of Poultry Farmers:Based on Field Research of Farmers in Six Provinces Such as Shandong and Anhui[J].China Population,Resources and Environment,2013,23(10):169-176.]