焦士興， 王安周， 李青云， 李中軒，趙榮欽， 張 曼， 林璐霜

（1.安陽(yáng)師范學(xué)院資源環(huán)境與旅游學(xué)院，河南安陽(yáng) 455002； 2.洛陽(yáng)市第十九中學(xué)，河南洛陽(yáng) 471000；3.許昌學(xué)院城市與環(huán)境學(xué)院，河南許昌 461000； 4.華北水利水電大學(xué)測(cè)繪與地理信息學(xué)院，鄭州 450046）

水土資源是人類生產(chǎn)生活的基礎(chǔ)資源，其協(xié)調(diào)發(fā)展有利于實(shí)現(xiàn)區(qū)域糧食安全和可持續(xù)發(fā)展［1］. 伴隨著社會(huì)經(jīng)濟(jì)的發(fā)展，工業(yè)化、城鎮(zhèn)化進(jìn)程加快，導(dǎo)致農(nóng)業(yè)用地持續(xù)減少、水環(huán)境惡化等現(xiàn)象，并嚴(yán)重威脅農(nóng)村耕地的可持續(xù)利用［2-3］. 國(guó)外研究在水土資源生態(tài)占用、土地利用對(duì)水資源影響、糧食生產(chǎn)對(duì)水資源影響等方面較多，如AL-Mulali等采用生態(tài)足跡方法探討了可再生能源生產(chǎn)對(duì)水土占用帶來(lái)的影響［4］；Yira等依據(jù)水文模擬模型研究了土地利用變化對(duì)水資源產(chǎn)生的影響［5］；孫世坤等采用虛擬水流模型分析了糧食進(jìn)出口對(duì)區(qū)域水資源壓力［6］. 國(guó)內(nèi)學(xué)者研究集中在構(gòu)建模型、時(shí)空分異、耦合協(xié)調(diào)等方面，如杜捷構(gòu)建了寧夏水土資源評(píng)價(jià)指標(biāo)體系，研究其耦合關(guān)系［7］；姜秋香構(gòu)建了水土資源優(yōu)化配置體系，分析了三江平原水土資源空間分布及其特征［8］；馬歷等通過構(gòu)建水資源與耕地利用效益指標(biāo)體系，并以四川省為例探討了二者的耦合關(guān)系［9］.

綜上，有關(guān)水資源與耕地利用的相關(guān)研究主要集中于土地對(duì)水資源的影響［5-6］、時(shí)空差異分析［8］以及區(qū)域研究［4，7，9］等，但對(duì)二者耦合協(xié)調(diào)、相互影響研究較少. 水資源系統(tǒng)改善與耕地利用效益增長(zhǎng)存在雙向作用，開展二者相互影響、相互作用關(guān)系的研究更具有現(xiàn)實(shí)意義. 河南省作為農(nóng)業(yè)大省和糧食生產(chǎn)核心區(qū)，農(nóng)業(yè)耗水量多且受到“非農(nóng)業(yè)”用水?dāng)D壓，農(nóng)業(yè)水資源相對(duì)短缺. 本研究構(gòu)建了河南省耕地利用效益與水資源壓力的綜合評(píng)價(jià)體系，采用耦合協(xié)調(diào)度等模型，研究其耦合協(xié)調(diào)狀況及空間分異，提出了針對(duì)性的政策建議.

1 研究方法及數(shù)據(jù)來(lái)源

1.1 耕地利用效益指標(biāo)體系和評(píng)價(jià)模型

結(jié)合相關(guān)研究［10-12］，構(gòu)建了河南省經(jīng)濟(jì)等3方面14項(xiàng)指標(biāo)體系（表1）. 具體公式如下：

表1 河南省耕地利用效益指標(biāo)Tab.1 Indexes of cultivated land use benefit in Henan Province

式中：f(t)為耕地利用效益指數(shù)；f(a)、f(b)、f(c)分別表示耕地利用社會(huì)效益、經(jīng)濟(jì)效益和生態(tài)效益；Wj表示第j個(gè)指標(biāo)權(quán)重；i表示第i年；j表示第j個(gè)指標(biāo)；Yij表示i年第j個(gè)指標(biāo)的標(biāo)準(zhǔn)化值.

1.2 水資源綜合壓力評(píng)價(jià)模型

1.2.1 標(biāo)準(zhǔn)差法賦權(quán)重

為消除指標(biāo)單位差異，采用標(biāo)準(zhǔn)化處理指標(biāo)層［13］. 根據(jù)各指標(biāo)標(biāo)準(zhǔn)化值，計(jì)算指標(biāo)權(quán)重［14］. 具體公式如下：

式中：E(Xj)、σ(Xi)、Wj分別表示xij標(biāo)準(zhǔn)化后的均值、標(biāo)準(zhǔn)差和指標(biāo)權(quán)重；xij為第i年第j個(gè)指標(biāo)值；N表示指標(biāo)數(shù)量；Zij為經(jīng)過標(biāo)準(zhǔn)化處理的數(shù)值.

1.2.2 水資源綜合壓力評(píng)價(jià)模型

結(jié)合相關(guān)研究成果，構(gòu)建了包括水資源人口等4個(gè)壓力指標(biāo)的綜合壓力評(píng)價(jià)模型［9］.具體公式如下：

式中：P為水資源綜合壓力指數(shù)，該值越大，壓力越大；Wα、Wβ、Wγ、Wδ分別為水資源人口、經(jīng)濟(jì)、生態(tài)、技術(shù)壓力指數(shù)權(quán)重；Zi、Ze、Zc、Zt則分別為水資源的人口、經(jīng)濟(jì)、生態(tài)、技術(shù)壓力指數(shù)標(biāo)準(zhǔn)化值［13］. 水資源人口等壓力指數(shù)計(jì)算見參考文獻(xiàn)［9］.

1.3 耦合協(xié)調(diào)度計(jì)算

根據(jù)耕地利用效益和水資源綜合壓力指數(shù)值，計(jì)算耦合度及耦合協(xié)調(diào)度. 具體公式如下：

式中：O為耦合度；I為協(xié)調(diào)指數(shù)；X為兩個(gè)系統(tǒng)的耦合協(xié)調(diào)度；Z1和Z2分別是耕地利用效益和水資源綜合壓力評(píng)價(jià)值；a、b為待定系數(shù)，假設(shè)耕地利用效益與水資源壓力系統(tǒng)同等重要，則a=b=0.5［15］. 耦合度反映了耕地利用效益與水資源綜合壓力相互作用程度，耦合協(xié)調(diào)度則探討了其動(dòng)態(tài)演化過程. 耦合度、耦合協(xié)調(diào)度類別標(biāo)準(zhǔn)見表2［9，16］.

表2 耦合度與耦合協(xié)調(diào)類別標(biāo)準(zhǔn)Tab.2 Standards for coupling degree and coupling coordination category

1.4 數(shù)據(jù)來(lái)源

數(shù)據(jù)包括有效灌溉面積、農(nóng)民人均可支配收入等18個(gè)指標(biāo). 其中人均水資源量、有效灌溉面積、水資源開發(fā)量、萬(wàn)元GDP用水量、水資源利用量5個(gè)數(shù)據(jù)來(lái)源于《河南省水資源公報(bào)》（2005—2018年）；農(nóng)民人均可支配收入、第一產(chǎn)業(yè)產(chǎn)值、人均GDP、農(nóng)業(yè)增加值、糧食總產(chǎn)量、耕地總面積、全年農(nóng)作物播種總面積、總?cè)丝跀?shù)、非農(nóng)業(yè)人口數(shù)、第一產(chǎn)業(yè)從業(yè)人口、農(nóng)業(yè)機(jī)械總動(dòng)力、土地總面積化肥使用量、農(nóng)藥使用量13個(gè)數(shù)據(jù)來(lái)源于《河南省統(tǒng)計(jì)年鑒》（2006—2019年）. 個(gè)別缺失數(shù)據(jù)依據(jù)相關(guān)資料整理獲得.

2 耕地利用效益與水資源綜合壓力指數(shù)時(shí)間演變特征

利用公式（1）～（5）計(jì)算河南省18市耕地利用效益指數(shù)和水資源綜合壓力指數(shù)，具體如圖1、圖2所示.

由圖可知，河南省耕地利用效益發(fā)展態(tài)勢(shì)較好，均呈現(xiàn)增長(zhǎng)趨勢(shì)，但增長(zhǎng)速度不同. 依據(jù)耕地利用效益指數(shù)，運(yùn)用SPSS軟件進(jìn)行聚類分析，將其分為快速增長(zhǎng)和緩慢增長(zhǎng)兩類（圖1）. 快速增長(zhǎng)的城市為鄭州、鶴壁、濮陽(yáng)、南陽(yáng)、信陽(yáng). 其中鄭州耕地利用效益指數(shù)由2005年的0.21增長(zhǎng)到2018年的0.769，增速最快，原因在于在經(jīng)濟(jì)與科技的推動(dòng)下，社會(huì)、經(jīng)濟(jì)和生態(tài)效益均呈現(xiàn)高速增長(zhǎng)，增速分別為185.6%、758.6%、184.7%，從而使耕地利用效益指數(shù)不斷提高. 緩慢增長(zhǎng)的城市有商丘、開封、洛陽(yáng)、平頂山、安陽(yáng)、新鄉(xiāng)、焦作、許昌、漯河、三門峽、周口、駐馬店、濟(jì)源. 其中商丘市耕地利用效益指數(shù)由2005年的0.478增長(zhǎng)到2018年的0.629，增速最慢，原因在于近年來(lái)商丘市耕地大量被占用轉(zhuǎn)變?yōu)樯虡I(yè)用地［17］，生態(tài)效益降低（降低354.9%），但社會(huì)效益和經(jīng)濟(jì)效益有所增加，耕地利用效益指數(shù)小幅度增長(zhǎng)，但增速較慢.

圖1 2005—2018年河南省耕地利用效益態(tài)勢(shì)Fig.1 Cultivated land use benefit index of Henan Province from 2005 to 2018

河南省水資源壓力區(qū)域差異較大，但整體呈波動(dòng)增長(zhǎng)態(tài)勢(shì). 依據(jù)水資源綜合壓力指數(shù)增長(zhǎng)速度，運(yùn)用SPSS軟件進(jìn)行聚類分析，將其分為快速增長(zhǎng)和緩慢增長(zhǎng)兩類（圖2）. 快速增長(zhǎng)的城市有鄭州、開封、洛陽(yáng)、安陽(yáng)、新鄉(xiāng)、平頂山、鶴壁等7市，其中鄭州水資源壓力指數(shù)由2005年的0.04增長(zhǎng)到2018年的1.10，增長(zhǎng)幅度較大，原因是近年來(lái)經(jīng)濟(jì)高速發(fā)展，水資源經(jīng)濟(jì)、人口和生態(tài)壓力指數(shù)均不斷增大，導(dǎo)致綜合壓力指數(shù)波動(dòng)增長(zhǎng)，水資源利用效益有待提高. 緩慢增長(zhǎng)的城市有焦作、濮陽(yáng)、許昌、漯河、三門峽、南陽(yáng)、商丘、信陽(yáng)、周口、駐馬店、濟(jì)源等11市，其中濟(jì)源市水資源壓力指數(shù)從2005年的0.083增長(zhǎng)到2018年的0.120，增長(zhǎng)幅度相對(duì)較小，原因是濟(jì)源市水資源技術(shù)、人口、生態(tài)壓力均有小幅度的上升，而經(jīng)濟(jì)壓力指數(shù)則大幅下降，水資源壓力指數(shù)變化幅度較小.

圖2 2005—2018年河南省水資源綜合壓力態(tài)勢(shì)Fig.2 Comprehensive pressure index of water resources in Henan Province from 2005 to 2018

3 耕地利用效益與水資源壓力的耦合分析

3.1 耕地利用效益與水資源壓力綜合評(píng)價(jià)指數(shù)動(dòng)態(tài)分析

利用公式（1）～（5），計(jì)算河南省耕地利用效益指數(shù)和水資源綜合壓力指數(shù)，具體如圖3、4所示.

由圖3可知，河南省耕地利用效益與水資源壓力綜合評(píng)價(jià)指數(shù)變化趨勢(shì)相同，相關(guān)系數(shù)為0.973，表明二者具有較強(qiáng)的正向相關(guān)性，說(shuō)明耕地利用效益增長(zhǎng)導(dǎo)致水資源壓力不斷增長(zhǎng)，且水資源壓力指數(shù)波動(dòng)幅度相對(duì)較大.

圖3 2005—2018年河南省耕地利用效益與水資源綜合壓力評(píng)價(jià)指數(shù)Fig.3 Evaluation index of cultivated land use efficiency and comprehensive pressure on water resources in Henan Province from 2005 to 2018

河南省耕地利用效益指數(shù)呈波動(dòng)增加趨勢(shì)，但波動(dòng)幅度較小. 耕地利用的社會(huì)、經(jīng)濟(jì)和生態(tài)效益權(quán)重依次為0.338、0.275、0.387，表明生態(tài)效益對(duì)耕地利用效益的影響最大. 耕地利用效益指數(shù)由2005 年的0.236波動(dòng)增長(zhǎng)至2018年的0.692，年均增長(zhǎng)0.035. 2006年后，隨著國(guó)家“十一五”計(jì)劃的貫徹執(zhí)行，河南省加強(qiáng)了耕地資源的管理，耕地利用水平不斷提高，利用方式不斷優(yōu)化，復(fù)種指數(shù)、灌溉指數(shù)均有大幅增加［18］.

河南省水資源壓力綜合評(píng)價(jià)指數(shù)呈增加趨勢(shì)，但波動(dòng)幅度較大. 水資源利用壓力不斷增加，從2005年的0.338 增至2018 年的0.738，年均增長(zhǎng)0.031，且水資源人口、經(jīng)濟(jì)、生態(tài)、技術(shù)壓力的權(quán)重依次為0.244、0.280、0.257、0.219，表明水資源受經(jīng)濟(jì)影響較大. 隨著用水結(jié)構(gòu)優(yōu)化和管理水平的提升，水資源綜合利用的水平也不斷提高，但持續(xù)增長(zhǎng)的經(jīng)濟(jì)、生態(tài)、人口等用水導(dǎo)致水資源壓力不斷增大. 2013年水資源綜合壓力指數(shù)達(dá)到低值（0.410），原因是2013年水資源生態(tài)、技術(shù)壓力均達(dá)到2005—2018年中最高值，分別為0.751和0.591，且綜合用水效率在2005—2017年最低［19］，因此水資源壓力較高，其綜合壓力指數(shù)則相應(yīng)較低.

3.2 耕地利用效益與水資源綜合壓力耦合度及耦合協(xié)調(diào)度動(dòng)態(tài)分析

利用公式（6）～（8）計(jì)算河南省耕地利用效益與水資源綜合壓力耦合度及耦合協(xié)調(diào)度，具體如圖4所示.

由圖4可知，河南省耕地利用效益與水資源綜合壓力耦合度及耦合協(xié)調(diào)度均呈增長(zhǎng)態(tài)勢(shì). 作為農(nóng)業(yè)大省和糧食生產(chǎn)核心區(qū)，伴隨著耕地利用效益的提高，河南農(nóng)業(yè)用水需求增大，二者耦合度從2005年低水平耦合（0.262）發(fā)展到2018年磨合時(shí)期（0.661），耦合協(xié)調(diào)度則從中度失調(diào)（0.274）發(fā)展到低度耦合協(xié)調(diào)（0.687）.2012—2013年，河南省的耕地利用效益與水資源綜合壓力耦合度及耦合協(xié)調(diào)度呈現(xiàn)下降趨勢(shì)，原因是2012年是“十二五”規(guī)劃的開局之年，伴隨著經(jīng)濟(jì)新常態(tài)調(diào)整以及農(nóng)業(yè)方面的供給側(cè)改革，經(jīng)濟(jì)發(fā)展變動(dòng)幅度較大，耕地利用效益略有下降. 但隨著水資源保護(hù)政策不斷發(fā)展完善，水污染狀況趨于減緩，且2012—2013年河南省水資源總量最大（256×108m3），因此2012—2013年水資源壓力有所下降.

圖4 2005—2017年河南省耕地利用效益與水資源綜合壓力耦合度及耦合協(xié)調(diào)度Fig.4 Coupling degree and coupling coordination degree between cultivated land use efficiency and comprehensive pressure of water resources in Henan Province from 2005 to 2017

3.3 耕地利用效益與水資源綜合壓力耦合協(xié)調(diào)度的時(shí)空分異

利用公式（6）～（8）計(jì)算河南18市耕地效益與水資源綜合壓力耦合度及耦合協(xié)調(diào)度，具體如表3所示.

由表3可以看出，鄭州、洛陽(yáng)、鶴壁和南陽(yáng)耕地利用效益與水資源綜合壓力耦合協(xié)調(diào)度由2005年的嚴(yán)重失調(diào)轉(zhuǎn)為2018年的低度失調(diào)；開封、安陽(yáng)、周口、駐馬店、許昌、漯河、商丘、三門峽和濟(jì)源由中度失調(diào)過渡到低度失調(diào)；新鄉(xiāng)和濮陽(yáng)由嚴(yán)重失調(diào)轉(zhuǎn)為中度失調(diào)；平頂山、焦作則均處于中度失調(diào)；信陽(yáng)由嚴(yán)重失調(diào)轉(zhuǎn)為弱度失調(diào). 河南省18市耕地利用效益與水資源綜合壓力耦合協(xié)調(diào)度整體呈增長(zhǎng)趨勢(shì)，原因是伴隨著經(jīng)濟(jì)新常態(tài)穩(wěn)步發(fā)展，河南省注重保護(hù)和利用耕地，耕地利用效益不斷增加，同時(shí)加大對(duì)水資源的經(jīng)濟(jì)、生態(tài)、技術(shù)等方面投資，水資源壓力不斷減小，促使二者耦合協(xié)調(diào)度不斷提升.

由表3 還可看出，河南省18 市耕地利用效益與水資源綜合壓力耦合協(xié)調(diào)度差異較小. 鄭州、許昌、洛陽(yáng)、三門峽、漯河、開封、鶴壁、商丘、周口、南陽(yáng)、駐馬店、信陽(yáng)、平頂山、安陽(yáng)、新鄉(xiāng)、焦作、濮陽(yáng)17市耦合協(xié)調(diào)度年均值為0.207～0.284，處于中度失調(diào)；濟(jì)源市則為0.309，處于低度失調(diào). 中度失調(diào)，說(shuō)明耕地利用效率提高對(duì)水資源壓力的影響較小，但人口壓力、生態(tài)壓力較大，如鄭州2018 年人口壓力、生態(tài)指數(shù)值分別為0.980 和0.955. 濟(jì)源的原因是隨著經(jīng)濟(jì)的發(fā)展，耕地利用社會(huì)和生態(tài)效益迅速增加［20］，耕地有效灌溉用水面積增加了69.4%，且水資源總量在河南18市中最少，水資源壓力指數(shù)較大.

表3 河南省18市耕地利用效益與水資源綜合壓力的耦合協(xié)調(diào)度Tab.3 Coupling coordination degree of cultivated land use benefit and comprehensive pressure of water resources in 18 cities of Henan Province

4 結(jié)論與政策建議

4.1 結(jié)論

1）河南省耕地利用效益發(fā)展態(tài)勢(shì)較好但增速度不同，水資源壓力區(qū)域差異較大且整體呈現(xiàn)波動(dòng)增長(zhǎng)趨勢(shì).

2）河南省耕地利用效益增長(zhǎng)導(dǎo)致水資源壓力不斷增長(zhǎng)，二者具有正向相關(guān)性，且水資源壓力波動(dòng)幅度相對(duì)較大.

3）河南省耕地利用效益與水資源壓力均呈波動(dòng)增長(zhǎng)態(tài)勢(shì)，耦合度從低水平耦合發(fā)展到磨合階段，耦合協(xié)調(diào)度則從中度失調(diào)發(fā)展到低度耦合協(xié)調(diào).

4）河南省18 市耕地利用效益與水資源壓力的耦合協(xié)調(diào)度均呈增長(zhǎng)趨勢(shì)，但增長(zhǎng)幅度不同且區(qū)域差異較小.

4.2 政策建議

1）加大技術(shù)和資金投入，完善灌溉設(shè)施，發(fā)展節(jié)水、節(jié)地農(nóng)業(yè)；適當(dāng)壓縮農(nóng)業(yè)用水規(guī)模，優(yōu)化水資源配置結(jié)構(gòu)，保障區(qū)域用水安全；在提高耕地利用效率的同時(shí)，強(qiáng)化水資源的再生循環(huán)與重復(fù)利用，提高公民節(jié)約用水意識(shí)，減少水資源浪費(fèi)，緩解水資源壓力.

2）推進(jìn)農(nóng)業(yè)綠色轉(zhuǎn)型和產(chǎn)業(yè)化發(fā)展，優(yōu)化農(nóng)業(yè)結(jié)構(gòu)，建設(shè)高標(biāo)準(zhǔn)農(nóng)田，提高農(nóng)業(yè)發(fā)展水平，確保耕地利用效益持續(xù)向好. 因地制宜，制定相應(yīng)的水資源保護(hù)政策，促進(jìn)水資源利用；多采用滴灌噴灌技術(shù)，提高農(nóng)業(yè)用水效率，促進(jìn)水土資源優(yōu)化配置.

3）統(tǒng)籌耕地資源和水資源開發(fā)利用機(jī)制，加強(qiáng)耕地?cái)?shù)量保護(hù)，提高科學(xué)耕作水平，依托科技支撐，高質(zhì)量和集約化發(fā)展農(nóng)業(yè)；不斷優(yōu)化水資源管理和利用，加強(qiáng)用水管理，提高水資源的綜合利用，促進(jìn)水資源系統(tǒng)進(jìn)一步優(yōu)化.

4）依據(jù)耦合協(xié)調(diào)類型，因地制宜制定區(qū)域政策，積極發(fā)展節(jié)水耐旱作物，培養(yǎng)蓄水聚肥技術(shù)，提高土壤自身的持水、蓄水能力，發(fā)揮“土壤水庫(kù)”的功能，促進(jìn)水資源的高效利用，多措并舉，以實(shí)現(xiàn)耕地利用效益和水資源的高效利用二者的耦合協(xié)調(diào)發(fā)展.