萬生新,王悅泰

?

基于DPSIR模型的沂河流域水生態(tài)安全評價方法

萬生新1,王悅泰2

1. 咸陽師范學院經(jīng)濟與管理學院, 陜西 咸陽 712000 2. 山東省新泰市泉溝鎮(zhèn)農(nóng)技站, 山東 新泰 271207

準確對流域水生態(tài)安全進行定量評價，明確流域水生態(tài)安全的限制因素。對DPSIR(驅動力-壓力-狀態(tài)-影響-響應)模型中“狀態(tài)”和“影響”指標的選取進行辨析后建立了流域水生態(tài)安全評價多層次指標體系，采用改進的群組AHP法確定各層指標權重并最終計算水生態(tài)安全指數(shù)。以沂河流域臨沂段為例，對該流域2015年水生態(tài)安全狀況進行評價。結果表明：流域內水生態(tài)安全狀態(tài)良好區(qū)域占45.2%，一般區(qū)域占54.8%，整體處在一般區(qū)間頂部，其水生態(tài)安全指數(shù)為0.573。水生態(tài)系統(tǒng)結構和功能仍較穩(wěn)定，但面臨著流域內經(jīng)濟水平亟須提高的威脅以及施肥強度過大的壓力，在生態(tài)投入與措施采取中，恢復水生態(tài)系統(tǒng)沿岸濕地是重點工作。

水生態(tài)安全; 評價; DPSIR模型; 沂河流域

20世紀80年代以來，各領域學者從全球安全、國家安全、以及人類自身安全的角度認識到了生態(tài)安全的重要性[1]，并不斷發(fā)展生態(tài)安全評價的理論與方法，其評價指標體系也從單因子發(fā)展到了多指標綜合體系[2-6]。90年代以來，我國學者在區(qū)域、流域、資源等生態(tài)安全評價方面取得了重大成果[7-10]，形成了以區(qū)域為核心，流域、生態(tài)脆弱區(qū)為輔的生態(tài)安全評價研究格局[11]。但目前為止，由于對生態(tài)安全概念及評價內容等方面沒有達成統(tǒng)一，所以生態(tài)安全評價的研究依然處在指標體系的構建和評價方法的探究階段[12]。水循環(huán)是流域內物質、能量循環(huán)的重要承載者，對流域的植被、氣候產(chǎn)生重要影響。水體中溶解性物質的侵蝕是造成流域內水土流失的重要因素。水力、物種多樣性等為流域內居民提供了巨大服務功能。但伴隨社會、經(jīng)濟的發(fā)展，流域水生態(tài)正在遭受人為的強烈干擾，另一方面，關于流域水生態(tài)的安全評價鮮有報道。

水源不僅對于是動植物的生存必需，同時也是生產(chǎn)生活的必需品，在人類社會生活方面起著無可替代的作用，對于江河等水流域而言，其生態(tài)安全直接制約著人類社會生活的方方面面[13-15]，人類活動或多或少都會對水域產(chǎn)生一定的影響，無論是有利影響還是不利影響，都難以避免地對水域生態(tài)結構產(chǎn)生一定作用，只有水域能夠穩(wěn)定地為人類社會生產(chǎn)生活提供足夠的支撐，才能說明其具有穩(wěn)定的可持續(xù)發(fā)展能力，其對外界威脅的應對能力正是其生態(tài)安全的充分體現(xiàn)，也是高于生態(tài)健康的體現(xiàn)，這也是生態(tài)學研究中常用的流域水生態(tài)安全。這種分析方法側重于水域的生態(tài)安全，利用水的因果聯(lián)系DPSIR來進行多層次的評價；之后通過AHP法對指標權重進行確定，然后對指標數(shù)據(jù)進行標準化賦值，并在權重的基礎上對影響水生態(tài)安全的因素予以明晰；此外，通過實例分析的方式對沂河流域臨沂段開展水域安全方面的探究，從而為提升水域安全水平提供有益參考和借鑒。

1 評價方法

1.1 DPSIR模型

該模型評價建立在以往的PSR模型之上，這種壓力-狀態(tài)-響應關系鏈的評價方法在環(huán)境問題分析中效果明顯，該模型在此基礎上融入了“驅動力”和“影響”指標，從而更加全面地分析生態(tài)環(huán)境社會問題；在這種因果關系鏈分析的方式之下，人類社會活動對生態(tài)環(huán)境的相互作用及影響能夠多層次的加以分析評價[16]，不僅僅社會生產(chǎn)生活能夠顯著影響水域結構及安全等方面，而水域生態(tài)的變化也會直接對社會生活產(chǎn)生反作用，二者之間具有明顯的關系鏈，在該模型分析之下通過“狀態(tài)”分析能夠對水域當前狀態(tài)有效識別，同時結合“驅動力”、“壓力”指標分析來對其面臨的威脅及壓力加以識別，在“影響”分析機制下判別其對相應生態(tài)及社會活動等方面的制約效果，“響應”識別能夠分析維持生態(tài)安全所需保障；綜合來看，該模型能夠多方面、深層次地開展水域生態(tài)安全分析。

在DPSIR模型的應用中，“狀態(tài)”和“影響”指標的選取一直較難區(qū)分。例如，在流域生態(tài)安全評價方面：徐藝揚等[17]將森林覆蓋率定義為流域生態(tài)的“狀態(tài)”指標，而張松等[11]認為森林覆蓋率是河流健康狀態(tài)造成的影響，所以將其定義為“影響”指標；姚遠[18]將生物多樣性定義為流域生態(tài)的“狀態(tài)”指標，而趙振亞等[19]認為生物多樣性為植被蓋度等狀態(tài)造成的“影響”。在水資源安全評價方面：肖新成等[20]將水質指標定義為水資源的“狀態(tài)”，而陳洋波等[21]認為水質的下降是水資源超載影響的結果，所以將水質指標列入“影響”范疇。究其原因，這是由于評價對象內部因果關系的多重性和復雜性導致的，即評價對象的水質指標既是水資源超載影響的結果，同時水質又作為狀態(tài)對水生生物、植被覆蓋等要素產(chǎn)生影響。為明確流域水生態(tài)安全評價中“狀態(tài)”與“影響”指標的選擇，該評價方法對于水生態(tài)系統(tǒng)內部指標僅關注評價時刻的靜態(tài)狀態(tài)，即忽略水生態(tài)系統(tǒng)內部指標因果關系的作用過程，皆定義為“狀態(tài)”指標。“影響”則主要關注水生態(tài)系統(tǒng)此刻狀態(tài)對與其緊密聯(lián)系的生態(tài)系統(tǒng)以及社會經(jīng)濟、人類健康造成的影響。

圖 1 DPSIR模型示意

1.2 評價單元

在進行流域水生態(tài)整體評價的同時也要考慮其空間差異性，需將評價對象劃分為次級單元進行評價。以行政區(qū)劃作為次級評價單元具有可操作性強、易于將評價結果反饋于管理的優(yōu)勢。依據(jù)水系和流域邊界確定次級評價單元在體現(xiàn)空間差異性的同時不會破壞流域水生態(tài)的連續(xù)性[13]。在評價時可從可操作性等方面選取合適的次級評價單元。

1.3 評價指標體系的構建

對水域生態(tài)分析建立在全面、可靠的指標分析基礎上，前提是準確把控生態(tài)安全內涵，之后在此基礎上開展相應的研究分析；在水域生態(tài)分析過程中不僅僅具有“目標”、“準則”層次，還要有“要素”、“指標”層次，在合理、全面的基礎上選取的指標共22個，具體如表1所示。

1.4 權重的確定

各指標在評價過程中的重要程度不同，該評價方法采用改進的群組AHP法確定各層指標權重。主要步驟為：①各位專家構建兩兩判斷矩陣獲得個體排序向量；②根據(jù)專家所構造判斷矩陣的一致性程度，并結合聚類分析思想得到每位專家的權重；③加權平均得到綜合排序向量，即指標權重。詳細計算過程見文獻[22]，所得權重見表1。

表 1 水生態(tài)安全評價指標體系及權重

1.5 評價指標標準化處理

考慮到指標的取值性難度以及研究的可操作性，同時考慮到不同指標之間的量綱差異較大[23]，因此在本研究過程中對流域的綜合評析采取量化分級的形式，根據(jù)實際測量值的優(yōu)劣分為優(yōu)秀、良好、一般、較差、差五級，并進行相應的標準化賦值，具體定義等如表2所示。

（1）建設用地比例：在人類長期影響下形成的人為景觀格局，可表征流域內物質、能量循環(huán)的阻礙。以陸域面積中建設用地的比例表示。

（2）水資源開發(fā)利用強度：水資源的過度開發(fā)會影響水文循環(huán)進而導致流域內濕地面積減少、過渡帶退化、湖泊萎縮等負面反應[29]。

區(qū)域水資源開發(fā)強度=年平均農(nóng)業(yè)、工業(yè)以及生態(tài)用水量/區(qū)域水資源的總量×100% (1)

（3）工業(yè)廢水排放強度：表征流域內點源污染壓力，以評價年單位面積的工業(yè)廢水排放量表示。

（4）施肥強度：流域內的水資源是農(nóng)業(yè)發(fā)展的有利條件，但較低的化肥利用率使農(nóng)業(yè)成為非點源污染的罪魁禍首。以單位農(nóng)作物播種面積(公頃)的化肥施用量(折純、千克)表示。

（5）河道連通性：水利工程以其隔斷效應干擾了水生生物的洄游及不同水域水生生物群體間的遺傳交流，致使生物多樣性降低[30-31]。以每公里的閘壩、水電站等水利工程個數(shù)表示。

（6）平水期水量狀況：以平水期水流淹沒河道范圍的百分比來表示評價目標的生態(tài)需水滿足度。

（7）消落帶濕地退化指數(shù)：消落帶濕地可有效攔截面源污染。

消落帶濕地退化指數(shù)=(歷史濕地面積數(shù)據(jù)-當前濕地面積)/目前濕地面積×100% (2)

（8）底泥重金屬風險指數(shù)：重金屬是富集性污染物，底泥是水中重金屬遷移的歸宿，也是造成水中重金屬二次污染的源頭。

式中：為潛在生態(tài)風險指數(shù)；為第種重金屬的潛在生態(tài)風險系數(shù)；為第種重金屬的毒性響應系數(shù)；為第種重金屬的污染系數(shù)；為第種重金屬的實測值；為第種重金屬的參比值，可選用當?shù)刂亟饘俦尘爸怠?/p>

（9）水生生物可以綜合反映水生態(tài)系統(tǒng)的健康狀況。該評價方法選擇魚類物種多樣性綜合指數(shù)，詳細計算見《流域健康評估指南》[23]。

（10）森林覆蓋率：森林是豐富的物質、能量、資源庫，是陸域生態(tài)系統(tǒng)的重要組成，以其表征水生態(tài)對陸域生態(tài)系統(tǒng)的影響。

森林覆蓋率=森林面積/陸域總面積×100% (4)

（11）水生態(tài)公眾滿意度：流域內生活居民對水生態(tài)現(xiàn)狀、措施、宣傳教育的滿意程度，通過調查問卷的方式實施。

（12）退耕還濕地率：濕地是水生態(tài)系統(tǒng)的屏障，能有效控制面源污染，退耕還濕具有重要意義。以評價對象當年水生態(tài)系統(tǒng)沿岸退耕面積比例表示。

表 2 流域水生態(tài)安全評價指標量化標準

1.6 水生態(tài)安全指數(shù)計算

式中，X表示第個準則層的第個指標的標準化值，B為第個準則層的計算標準化值，W為第個指標的權重。

式中，W為第個準則層的權重。

1.7 水生態(tài)安全等級劃分

按照優(yōu)秀、良好、一般、較差、差5個等級劃分水生態(tài)安全等級（表3）

表 3 水生態(tài)安全等級

本研究過程中涉及的統(tǒng)計數(shù)據(jù)取自于臨沂市統(tǒng)計年鑒及職能部門數(shù)據(jù)；通過遙感影像來獲取遙感數(shù)據(jù)，并采取高分一號；課題組2015年采樣數(shù)據(jù)作為環(huán)境監(jiān)測數(shù)據(jù)。

2 沂河流域(臨沂段)水生態(tài)安全評價

2.1 研究區(qū)域

臨沂因地臨沂河得名，位于北緯34°22′~36°13′，東經(jīng)117°24′~119°11′之間。臨沂境內沂河流經(jīng)沂水縣、沂南縣、蒙陰縣、平邑縣、費縣、蘭山區(qū)、河東區(qū)、羅莊區(qū)、郯城縣，全長約333 km，流域面積11820 km2。近年來，臨沂市社會經(jīng)濟發(fā)展迅速，水利興建、耕地向建設用地轉變、水資源大量開發(fā)以及點源、面源污染導致水體出現(xiàn)富營養(yǎng)化，水生態(tài)安全面臨威脅。評價對象及其次級評價單元劃分見圖2。

圖 2 評價區(qū)域

2.2 數(shù)據(jù)來源

研究區(qū)數(shù)據(jù)包括統(tǒng)計數(shù)據(jù)、遙感影像數(shù)據(jù)和環(huán)境監(jiān)測數(shù)據(jù)。其中統(tǒng)計數(shù)據(jù)來源于山東省臨沂市各相關職能部門及統(tǒng)計年鑒，遙感影像為覆蓋全流域的高分一號數(shù)據(jù)，環(huán)境監(jiān)測數(shù)據(jù)來源于課題組2015年的采樣監(jiān)測。

3 結果分析

3.1 流域整體分析

沂河流域(臨沂段)水生態(tài)安全指數(shù)為0.573，處在“一般”區(qū)間頂部，水生態(tài)系統(tǒng)結構尚完整、功能可保證，但已出現(xiàn)一定威脅。繪制準則層“驅動力、壓力、狀態(tài)、影響、響應”的得分雷達圖，如圖3。由圖3可知：①“驅動力”得分為0.523，處在“一般”區(qū)間，在準則層中得分較低。究其原因，較低的人均GDP(標準化值為0.290)將會對該流域水生態(tài)產(chǎn)生較大威脅。②“壓力”、“狀態(tài)”得分分別為0.646和0.600，處在“良好”區(qū)間，水生態(tài)面臨的直接壓力較小、處在較健康狀態(tài)。施肥強度(標準化值為0.520)是目前該流域水生態(tài)面臨的較大直接壓力。周期性淹水以及人類活動的影響使得消落帶內濕地數(shù)量減少，消落帶濕地退化指數(shù)(標準化值為0.327，處在較差區(qū)間)為“狀態(tài)”的限制因素。③“影響”的計算得分為0.459，處在“一般”區(qū)間，是準則層中的得分最低項。除水災害損失占GDP比例的標準化值處在良好區(qū)間，其余各指標標準化值皆處在一般區(qū)間。

圖 3 準則層評價結果

3.2 空間差異性分析

次級評價單元間具有差異性，水生態(tài)安全狀況“良好”區(qū)域面積占45.2%，“一般”區(qū)域面積占54.8%，詳見圖4。從空間差異性來看：①處在上游的沂水縣、蒙陰縣，政策上限制高強度工業(yè)、推行生態(tài)化農(nóng)業(yè)，在較小壓力下，水生態(tài)系統(tǒng)結構、功能穩(wěn)定，水生態(tài)安全狀況處在“良好”區(qū)間。②流域中上游的費縣、平邑縣、沂南縣，水環(huán)境狀態(tài)與生物種群狀態(tài)成為水生態(tài)安全的限制因素，水生態(tài)安全狀況處在“一般”區(qū)間。③蘭山區(qū)與河東區(qū)處在流域中游，近年來政府部門加大生態(tài)投入，嚴格執(zhí)行環(huán)保措施，水生態(tài)安全狀態(tài)明顯改善，處在“良好”區(qū)間。④位于中下游的羅莊區(qū)，重工業(yè)為其支柱產(chǎn)業(yè)，較大的水資源利用及工業(yè)廢水排放強度對水生態(tài)造成了較大壓力，并在此壓力下造成了水環(huán)境狀態(tài)和生物種群狀態(tài)的下降，水生態(tài)安全狀況處在“一般”區(qū)間。⑤郯城縣位于羅莊區(qū)下游，水環(huán)境狀態(tài)受到上游影響，水生態(tài)安全狀況處在“一般”區(qū)間。各次級評價單元的準則層得分及水生態(tài)安全指數(shù)見表4。

圖 4 次級評價單元水生態(tài)安全評價結果

表 4 各次級評價單元準則層得分及水生態(tài)安全指數(shù)

4 結論

1) DPSIR模型不僅可以識別水生態(tài)目前的健康狀態(tài)，又可預見水生態(tài)面臨的潛在威脅，以及在人類保護下水生態(tài)可持續(xù)發(fā)展的能力。DPSIR模型是進行水生態(tài)安全評價的適合模型。

2) 提出了DPSIR模型在水生態(tài)安全評價過程中“狀態(tài)”和“影響”指標的選擇原則：對于水生態(tài)系統(tǒng)內部指標，忽略其因果關系的作用過程，皆定義為“狀態(tài)”指標，“影響”指標則主要關注水生態(tài)系統(tǒng)此刻狀態(tài)對與其緊密聯(lián)系的生態(tài)系統(tǒng)以及社會經(jīng)濟、人類健康造成的影響。

3) 該方法以DPSIR模型為基礎構建包含22個指標的流域水生態(tài)安全評價指標體系，采用改進的群組AHP法確定各層指標權重，并建立指標標準化標準，形成系統(tǒng)的流域水生態(tài)安全評價方法。

4) 沂河流域臨沂段水生態(tài)安全指數(shù)為0.579。其中“狀態(tài)”得分為0.621，水生態(tài)處于較健康狀態(tài)，但人均GDP標準化值低至0.290，存在流域內經(jīng)濟水平亟待提高的威脅。在水生態(tài)當前健康狀態(tài)下，陸域生態(tài)系統(tǒng)中存在水土流失比例較大(標準化值為0.42)現(xiàn)象，水生態(tài)對與其緊密聯(lián)系的生態(tài)系統(tǒng)造成了一定消極影響。

[1] 肖篤寧,陳文波,郭福良.論生態(tài)安全的基本概念和研究內容[J].應用生態(tài)學報,2002,13(3):354-358

[2] Hacho? J, El?bieta Bondar-Nowakowska. Aquatic plants–based risk model for assessment of ecological safety of rivers[J]. Human and Ecological Risk Assessment: An International Journal, 2016,22(4):1065-1077

[3] Dos Santos FB, Esteves KE. A Fish-Based Index of Biotic Integrity for the Assessment of Streams Located in a Sugarcane-Dominated Landscape in Southeastern Brazil[J]. Environ Manage, 2015,56(2):532-548

[4] 高遠,蘇宇祥,亓樹財.沂河流域浮游植物與水質評價[J].湖泊科學,2008(4):544-548

[5] 凌紅波,徐海量,史薇,等.新疆瑪納斯河流域綠洲生態(tài)安全評價[J].應用生態(tài)學報,2009(9):2219-2224

[6] Barbara Sowińska-?wierkosz. Application of surrogate measures of ecological quality assessment: The introduction of the Indicator of Ecological Landscape Quality (IELQ)[J]. Ecological Indicators, 2017,7(3):224-234

[7] 萬利,陳佑啟,譚靖,等.北京郊區(qū)生態(tài)安全動態(tài)評價與分析[J].地理科學進展,2009(2):238-244

[8] 張銳,劉友兆.我國耕地生態(tài)安全評價及障礙因子診斷[J].長江流域資源與環(huán)境,2013(7):945-951

[9] 許妍,高俊峰,郭建科.太湖流域生態(tài)風險評價[J].生態(tài)學報,2013,33(9):2896-2906

[10] 金相燦,王圣瑞,席海燕.湖泊生態(tài)安全及其評估方法框架[J].環(huán)境科學研究,2012(4):357-362

[11] 張松,郭懷成,盛虎,等.河流流域生態(tài)安全綜合評估方法[J].環(huán)境科學研究,2012(7):826-832

[12] 高吉喜,張向暉,姜昀,等.流域生態(tài)安全評價關鍵問題研究[J].科學通報,2007(S2):216-224

[13] 張遠,高欣,林佳寧,等.流域水生態(tài)安全評估方法[J].環(huán)境科學研究,2016(10):1393-1399

[14] 王耕,王利,吳偉.區(qū)域生態(tài)安全概念及評價體系的再認識[J].生態(tài)學報,2007(4):1627-1637

[15] Song Y, Chang RY, Yu ZD,. The Relationship between Ecosystem Health and Ecological Safety[C]. Information Engineering Research Institute, USA. Proceedings of 2012 Asian Pacific Conference on Energy, Environment and Sustainable Development(APEESD 2012 V260-261). Information Engineering Research Institute, USA: 2012:781-785

[16] 曹紅軍.淺評DPSIR模型[J].環(huán)境科學與技術,2005(S1):110-111,126

[17] 徐藝揚,錢敏蕾,李響,等.基于DPSIR的太平湖流域(黃山區(qū))生態(tài)安全綜合評估[J].復旦學報:自然科學版,2015(4):407-415

[18] 姚遠.基于DPSIR模型的流域環(huán)境變遷與生態(tài)安全指標體系研究[J].安徽農(nóng)業(yè)科學,2012(12):7345-7348

[19] 趙振亞,姬寶霖,宋小園,等.基于DPSIR-TOPSIS模型的黃河內蒙古段十大孔兌綜合治理生態(tài)環(huán)境效應評估[J].水土 保持通報,2016(1):143-150

[20] 肖新成,何丙輝,倪九派,等.農(nóng)業(yè)面源污染視角下的三峽庫區(qū)重慶段水資源的安全性評價——基于DPSIR框架的 分析[J].環(huán)境科學學報,2013(8):2324-2331

[21] 陳洋波,陳俊合,李長興,等.基于DPSIR模型的深圳市水資源承載能力評價指標體系[J].水利學報,2004(7):98-103

[22] 夏萍,汪凱,李寧秀,等.層次分析法中求權重的一種改進[J].中國衛(wèi)生統(tǒng)計,2011(2):151-154,157

[23] 中國環(huán)境科學研究院.流域生態(tài)健康評估技術指南[R].北京:中國環(huán)境科學研究院,2013

[24] 解雪峰,吳濤,肖翠,等.基于PSR模型的東陽江流域生態(tài)安全評價[J].資源科學,2014(8):1702-1711

[25] 毛志剛,谷孝鴻,陸小明,等.太湖東部不同類型湖區(qū)疏浚后沉積物重金屬污染及潛在生態(tài)風險評價[J].環(huán)境科 學,2014(1):186-193

[26] 劉欽普,林振山,周亮.山東省化肥使用時空分異及潛在環(huán)境風險評價[J].農(nóng)業(yè)工程學報,2015(7):208-214

[27] 金中彥,鄭彥強,趙海生.基于DPSIR模型的嵐漪河流域生態(tài)系統(tǒng)安全評估[J].人民黃河,2012(3):54-56

[28] 陳華偉,黃繼文,張欣,等.基于DPSIR概念框架的水生態(tài)安全動態(tài)評價[J].人民黃河,2013(9):34-37,45

[29] 張軍民.瑪納斯河流域水文循環(huán)二元分化生態(tài)效應研究[J].干旱區(qū)資源與環(huán)境,2008(1):133-136

[30] 周小愿.水利水電工程對水生生物多樣性的影響與保護措施[J].中國農(nóng)村水利水電,2009(11):144-146

[31] 蘇國歡,沙永翠,熊鷹,等.大壩截流前后金沙江觀音巖水電站魚類群落功能多樣性的變化[J].長江流域資源與環(huán) 境,2015(6):965-970

The Method Assessing Water Eco-security in Yi River Based on DPSIR Model

WAN Sheng-xin1, WANG Yue-tai2

1.712000,2.271207,

In order to quantitatively assess the water eco-security and determine the limiting factors of it, this method established a multi-level index system by discriminating the selection of the "state" and "impact" in the DPSIR (Drives-Pressures-State-Impact-Responses) model. After establishing the index system, the improved group AHP method was used to determine the weight of each index and water eco-security index was calculated finally. Taking the Linyi section of Yi River watershed as an example, the water eco-security in 2015 was assessed. The results show that the good regions of assessment object accounted for 45.2%, the general areas accounted for 54.8%, and the water eco-security index of the assessment object is 0.573, at the top of the general interval. The structure and function of the water ecosystem are stable, but it is facing the threat of increasing the economic level and the pressure of excessive fertilization. The restoration of wetland in water ecosystem is priority for the process of ecological investment and measures taking.

water ecological security; assessment; DPSIR model; Yi River Watershed

X826

A

1000-2324(2019)03-0502-07

10.3969/j.issn.1000-2324.2019.03.032

2018-03-24

2018-04-11

教育部人文社會科學研究青年基金:社會資本、組織結構與農(nóng)民用水合作組織發(fā)展研究(13YJC790135);陜西省社會科學基金:基于成員退出行為的關中灌區(qū)農(nóng)民用水合作組織發(fā)展研究(2016D026)

萬生新(1975-),男,博士,副教授.主要研究方向:農(nóng)業(yè)與農(nóng)村發(fā)展. E-mail:287405457@qq.com