金玨雯，穆月英，楊 鑫

（1.中國農(nóng)業(yè)大學(xué)經(jīng)濟管理學(xué)院，北京 100083； 2.中國社會科學(xué)院農(nóng)村發(fā)展研究所，北京 100732）

中國糧食安全在取得巨大成就的同時，結(jié)構(gòu)性問題突出［1-2］。2020年，中國糧食生產(chǎn)實現(xiàn)歷史性的“十七連豐”，全國糧食總產(chǎn)量13 390億斤（1斤等于500 g），較2019年增長0.9%，連續(xù)6年保持在1.3萬億斤以上，完全有能力保障糧食和重要農(nóng)產(chǎn)品供給，中國人的飯碗牢牢端在自己手中。然而，中國糧食安全面臨供求結(jié)構(gòu)失衡的重大挑戰(zhàn)，主要表現(xiàn)為區(qū)域結(jié)構(gòu)失衡和品種結(jié)構(gòu)失衡。一方面，糧食生產(chǎn)持續(xù)向糧食主產(chǎn)區(qū)集中，糧食主銷區(qū)糧食供求缺口大，糧食區(qū)域貿(mào)易量不斷加大，進一步增加了市場調(diào)控壓力［3］。糧食主產(chǎn)區(qū)糧食產(chǎn)量占全國產(chǎn)量比重由2004年的72.7%上升至2020年的78.6%，而同期糧食主銷區(qū)、產(chǎn)銷平衡區(qū)的產(chǎn)量占比卻分別從20.0%、7.4%下降至17.1%、4.3%。另一方面，伴隨居民食物消費結(jié)構(gòu)的轉(zhuǎn)型升級，糧食品種結(jié)構(gòu)供給與需求匹配錯位問題突出。作為口糧的水稻和小麥的國內(nèi)產(chǎn)出始終維持在較高水平，但國內(nèi)消費增長乏力；玉米政策性去庫存結(jié)束，玉米市場出現(xiàn)供應(yīng)缺口；大豆進口量持續(xù)快速增長，嚴重依賴國際市場［4-6］。因此，聚焦于與糧食安全直接相關(guān)的糧食種植結(jié)構(gòu)調(diào)整問題具有重要的現(xiàn)實意義。

與此同時，以總量少、空間分布不均為特征的水資源短缺問題日益加重，成為糧食安全的關(guān)鍵制約因素。一方面，中國人均水資源占有量少，約為世界人均水平的1/4，被聯(lián)合國列為貧水國。另一方面，經(jīng)濟發(fā)展和居民收入提高使得工業(yè)用水和生活用水不斷增加，農(nóng)業(yè)用水空間不斷被擠占，農(nóng)業(yè)用水比例從2004年的65%下降至2020年的62%。此外，水資源空間分布嚴重不均［7］。2020年，作為中國主要糧食生產(chǎn)基地的北方地區(qū)，僅依靠36%的水資源生產(chǎn)了全國59%的糧食。

為緩解水資源短缺和空間分布不均問題，中國已修建了20余項大型跨流域調(diào)水工程，進行實體水調(diào)運，總長度超過7 200 km［8］。以南水北調(diào)工程為例，截至2019年底，南水北調(diào)東、中線一期工程累計調(diào)水量303.88億m3。對水資源空間布局的調(diào)控，勢必會緩解缺水地區(qū)的生活和產(chǎn)業(yè)用水壓力，與此同時，一個地區(qū)水的總供給的改變，又會對農(nóng)業(yè)產(chǎn)生何種影響，也是各界關(guān)注的問題［9-11］。

除了跨流域調(diào)水為代表的實體水流動，以區(qū)域間商品和服務(wù)貿(mào)易為載體的虛擬水流動也被認為會改變空間上的水資源配置［12-13］。Allan［14］首次提出虛擬水概念，指的是在生產(chǎn)產(chǎn)品和服務(wù)中所需要的水資源數(shù)量。他指出，因為糧食商品的生產(chǎn)需要消耗水資源，糧食貿(mào)易的背后隱含著看不見的水資源交易。中國糧食生產(chǎn)越來越向北方糧食主產(chǎn)區(qū)集中，北方地區(qū)每年通過糧食貿(mào)易向南方地區(qū)輸出了大量虛擬水，且越來越多的省份加入糧食輸入行列或擴大輸入行列，加劇了糧食生產(chǎn)水資源空間分布不均問題［15-16］。

跨流域調(diào)水和糧食區(qū)域貿(mào)易均會影響區(qū)域水資源供給與利用結(jié)構(gòu)，進而影響糧食種植結(jié)構(gòu)。那么，基于“實體水-虛擬水”流動視角，即以跨流域調(diào)水“南水北調(diào)”和糧食區(qū)域貿(mào)易“北水南運”為表現(xiàn)的實體水流動與虛擬水流動互為對照、共同作用，對區(qū)域水資源約束的影響如何？這種影響又會給糧食種植結(jié)構(gòu)帶來怎樣的作用？糧食輸出區(qū)和糧食輸入?yún)^(qū)水資源約束變化，對糧食種植結(jié)構(gòu)的影響是否存在異質(zhì)性以及非對稱的空間溢出效應(yīng)？為了給出水資源約束與糧食種植結(jié)構(gòu)調(diào)整的新解釋，為保障水資源和糧食安全可持續(xù)發(fā)展的政策制定提供決策參考，文章將對上述問題作進一步的研究。

1 文獻綜述

圍繞水資源約束與糧食種植結(jié)構(gòu)的關(guān)鍵議題，已有研究主要從實體水和虛擬水兩個視角展開，前者主要分析實體水短缺和實體水資源管理對糧食種植結(jié)構(gòu)影響問題；后者主要從糧食貿(mào)易與消費出發(fā)，重點分析糧食虛擬水流動驅(qū)動因素及效應(yīng)問題。已有研究表明，實體水短缺對不同糧食作物種植結(jié)構(gòu)的影響是不同的。李玉敏等［17］通過對全國10個省調(diào)研數(shù)據(jù)的分析發(fā)現(xiàn)，水資源越短缺，農(nóng)戶越傾向于種植對灌溉依賴程度低、需水量小的作物；水資源短缺會降低水稻和小麥的種植比例，提高玉米和豆類的種植比例。楊鑫等［18］研究表明，灌溉水壓力增加，促進了玉米產(chǎn)量比例大幅增加、豆類和薯類產(chǎn)量比例顯著下降。也有研究從緩解實體水短缺角度出發(fā)，分析了灌溉基礎(chǔ)設(shè)施對于糧食種植結(jié)構(gòu)的影響，發(fā)現(xiàn)農(nóng)田水利設(shè)施建設(shè)改善了灌溉條件，保障了土地集約經(jīng)營，提高了中國糧食種植比例［19］。在實體水資源管理方面，已有研究主要分析了農(nóng)業(yè)水價對糧食種植結(jié)構(gòu)的影響。目前農(nóng)業(yè)水價尚處于試點階段，大部分地區(qū)水價低廉，難以起到價格杠桿作用，對農(nóng)戶種植結(jié)構(gòu)的影響并不顯著［20-21］。易福金等［22］認為只有當農(nóng)業(yè)灌溉用水單價持續(xù)上漲到一定水平時，農(nóng)戶才可能調(diào)整農(nóng)作物種植結(jié)構(gòu)。許朗等［23］認為在一些較為干旱且農(nóng)業(yè)水價定價較高的地區(qū)，農(nóng)業(yè)水價才會發(fā)揮節(jié)水作用，促使農(nóng)戶調(diào)整種植結(jié)構(gòu)。在虛擬水視角下，以往研究多關(guān)注區(qū)域間糧食虛擬水流動格局、驅(qū)動因素以及所產(chǎn)生的資源、社會、經(jīng)濟等影響［24-26］。或者以虛擬水為工具，分析糧食種植結(jié)構(gòu)調(diào)整的節(jié)水效應(yīng)，認為糧食種植結(jié)構(gòu)調(diào)整節(jié)約了水資源消耗，但不同區(qū)域間存在異質(zhì)性［27-28］。

現(xiàn)有研究或單獨探討區(qū)域?qū)嶓w水資源約束對糧食種植結(jié)構(gòu)的影響，忽略了區(qū)域間水資源的流動性；或僅僅分析基于糧食區(qū)域貿(mào)易的虛擬水流動所帶來的影響，使得水資源約束對糧食種植結(jié)構(gòu)影響估計產(chǎn)生偏誤，原因如下：第一，未同時考慮實體水和虛擬水流動過程。糧食種植結(jié)構(gòu)不僅受到實體水資源約束的影響，也受到產(chǎn)后貿(mào)易與消費隱含的虛擬水流動的影響。第二，沒有考慮實體水和虛擬水區(qū)域流動產(chǎn)生的空間異質(zhì)性。實體水和虛擬水流動對于糧食輸出區(qū)和糧食輸入?yún)^(qū)糧食種植結(jié)構(gòu)的影響可能是不同的。

探究水資源約束對糧食種植結(jié)構(gòu)的影響，需要統(tǒng)籌考慮“實體水-虛擬水”流動過程。吳普特等［29-30］從水資源的自然、經(jīng)濟和社會屬性出發(fā)，提出了統(tǒng)籌考慮“實體水-虛擬水”流動過程的概念，包含三個層次：第一，水作為自然物質(zhì)，受自然規(guī)律、人工調(diào)配等作用，形成以“天然降水-廣義水資源-自然蒸散-經(jīng)濟社會利用-退水”路徑為特征的實體水流動過程；第二，水作為生產(chǎn)資料，形成“實體水攝取-實體水耗散-虛擬水嵌入”的基本路徑，即實體水向虛擬水轉(zhuǎn)化并嵌入產(chǎn)品的過程；第三，嵌入產(chǎn)品中的虛擬水隨區(qū)域間商品流通形成“虛擬水流動（商品流通）-虛擬水回用（商品回收再生產(chǎn)）”的虛擬水流動路徑。作為實體水調(diào)運工程之一的跨流域調(diào)水工程，已成為中國解決水資源時空分配不均勻問題和提高水資源配置能力的重要手段，而相較于其他實體水調(diào)運工程，跨流域調(diào)水工程具有投資成本高、建設(shè)周期長、影響范圍廣等特點［31-32］。與此同時，基于糧食區(qū)域貿(mào)易，每年有大量的糧食虛擬水從北方地區(qū)輸送到南方地區(qū)，同實體水的南水北調(diào)形成互補［16］。因此，文章借鑒吳普特等［29-30］提出的統(tǒng)籌考慮“實體水-虛擬水”流動過程概念，重點關(guān)注由于人工調(diào)配和產(chǎn)后市場貿(mào)易對水資源約束的改變，進一步將“實體水”流動聚焦于跨流域調(diào)水，“虛擬水”流動聚焦于糧食區(qū)域貿(mào)易，從“實體水-虛擬水”流動視角探究跨流域調(diào)水、糧食區(qū)域貿(mào)易對糧食種植結(jié)構(gòu)的影響。

文章基于“實體水-虛擬水”流動視角，聚焦于跨流域調(diào)水、糧食區(qū)域貿(mào)易及其共同作用對區(qū)域水資源約束的改變，進而探究其對糧食種植結(jié)構(gòu)的影響。相較于以往研究，文章有兩個方面的擴展：第一，統(tǒng)籌考慮實體水與虛擬水，從區(qū)域間水資源流動視角考察跨流域調(diào)水和糧食區(qū)域貿(mào)易對糧食種植結(jié)構(gòu)的影響；第二，采用糧食區(qū)域貿(mào)易空間權(quán)重矩陣的空間杜賓模型考察“實體水-虛擬水”流動對糧食輸出區(qū)與糧食輸入?yún)^(qū)糧食種植結(jié)構(gòu)的差異化影響，以及相應(yīng)的空間溢出效應(yīng)。文章的第二部分將闡述“實體水-虛擬水”流動影響糧食種植結(jié)構(gòu)的理論分析；第三部分介紹模型設(shè)定、變量選擇和數(shù)據(jù)來源；第四部分為實證結(jié)果與分析；第五部分為穩(wěn)健性檢驗；最后，第六部分為結(jié)論與啟示。

2 理論與事實

2.1 “實體水-虛擬水”流動過程的理論邏輯：實體水“南水北調(diào)”與虛擬水“北水南運”

根據(jù)虛擬水貿(mào)易理論，缺水地區(qū)可以通過貿(mào)易的方式從富水區(qū)進口水資源密集型產(chǎn)品緩解本地區(qū)的水資源短缺，實現(xiàn)社會經(jīng)濟發(fā)展的水資源和糧食供應(yīng)安全［14］。而中國糧食區(qū)域貿(mào)易“北糧南運”的虛擬水流動格局是水資源分布的逆向配置，即由缺水的北方地區(qū)向富水的南方地區(qū)進行轉(zhuǎn)移。這是由于虛擬水貿(mào)易的出發(fā)點和立足點僅注重了水資源要素稟賦方面，忽略了其他資源（如耕地）、生產(chǎn)技術(shù)水平等因素。從虛擬水流動和區(qū)域水資源分布角度來看，“北糧南運”和“南水北調(diào)”共存仍是今后發(fā)展的方向［24］。有必要從邏輯和事實認清“北糧南運”和“南水北調(diào)”格局的成因、長期性，以及其對糧食安全和水資源安全的影響。

根據(jù)比較優(yōu)勢理論，糧食區(qū)域貿(mào)易“北糧南運”格局是資源稟賦、社會經(jīng)濟、技術(shù)進步、宏觀政策等因素共同作用的結(jié)果。耕地資源“南少北多”是糧食增長中心“北上”的首要因素［33-34］。其次，南北方地區(qū)社會經(jīng)濟發(fā)展水平和發(fā)展模式差異，導(dǎo)致產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化調(diào)整方向不同［35］。相對于北方地區(qū)，南方水熱資源充足，生產(chǎn)條件較好，種植作物類型選擇較多，加上價格政策因素影響，種植經(jīng)濟作物收益高于種植糧食作物，導(dǎo)致糧食播種面積下降較多［36］。另一方面，南方地區(qū)相較于北方非農(nóng)產(chǎn)業(yè)發(fā)展迅速，導(dǎo)致耕地、資金、勞動等農(nóng)業(yè)資源向非農(nóng)產(chǎn)業(yè)轉(zhuǎn)移較多［37］。第三，農(nóng)業(yè)技術(shù)進步減輕了糧食生產(chǎn)對硬約束資源的依賴。地膜技術(shù)的廣泛應(yīng)用，改善了北方很多地區(qū)的耕作制度，復(fù)種指數(shù)明顯提高。良種研發(fā)、化肥和農(nóng)藥等投入要素的研發(fā)和推廣、節(jié)水灌溉技術(shù)采用、有效灌溉面積增加等減輕了北方各地區(qū)糧食生產(chǎn)資源稟賦不足的約束［38］。此外，國家出臺了一系列政策穩(wěn)定糧食主產(chǎn)區(qū)糧食生產(chǎn)與供給。劃定13個省份為糧食主產(chǎn)區(qū)，堅守18億畝（1畝=666.67 m2）耕地紅線，建立健全糧食安全省長責任制，劃定糧食生產(chǎn)功能區(qū)等，促使區(qū)域分工格局進一步固化，糧食生產(chǎn)持續(xù)向北方集中［39］?；谏鲜鲈颍谙喈旈L時間內(nèi)糧食區(qū)域貿(mào)易格局“北糧南運”特征不會改變，且越來越多省份成為糧食輸入?yún)^(qū)或擴大輸入量［15］。

“北糧南運”格局下，糧食輸出區(qū)水資源短缺問題日趨嚴峻，是糧食安全可持續(xù)發(fā)展面臨的重大挑戰(zhàn)。糧食虛擬水從糧食生產(chǎn)灌溉用水效率較高的地區(qū)流向較低地區(qū)，從缺水地區(qū)流向豐水地區(qū)，使得糧食輸出區(qū)水資源壓力指數(shù)由1.61上升至2.04［26］。虛擬水資源的流動緩解了南方地區(qū)的耕地短缺，提高了全國的水土資源生產(chǎn)效率，卻加大了水土不平衡的差距［40］。糧食虛擬水輸出區(qū)在保障國家整體利益及對輸入?yún)^(qū)的糧食安全作出重要貢獻的同時，因糧食生產(chǎn)擠占綠水資源而對陸地自然生態(tài)系統(tǒng)造成嚴重威脅［41］。

基于糧食區(qū)域貿(mào)易的虛擬水流動格局短期內(nèi)很難改變，而糧食虛擬水進口面臨的國際環(huán)境復(fù)雜多變，南水北調(diào)等跨流域調(diào)水工程的作用就顯得尤為突出［24］。盡管南水北調(diào)等跨流域調(diào)水工程并不是直接服務(wù)于農(nóng)業(yè)，但一定程度上彌補了北方糧食輸出區(qū)的虛擬水消耗，緩解了區(qū)域水資源壓力，對改進區(qū)域水土資源平衡、優(yōu)化糧食生產(chǎn)資源配置等起到積極的作用［40-41］?？缌饔蛘{(diào)水能夠促進區(qū)域生態(tài)環(huán)境改善，通過水量置換遏制受水區(qū)地下水嚴重超采的趨勢，并逐步退還被長期大量擠占的農(nóng)業(yè)及生態(tài)用水，保障國家糧食安全與水安全［42］。

根據(jù)上述理論邏輯，文章認為基于跨流域調(diào)水“南水北調(diào)”和糧食區(qū)域貿(mào)易“北水南運”的“實體水-虛擬水”流動格局將長期存在，其對區(qū)域水資源約束的改變不容忽視，進而顯著作用于糧食種植結(jié)構(gòu)。

2.2 “實體水-虛擬水”流動視角下跨流域調(diào)水、糧食區(qū)域貿(mào)易對糧食種植結(jié)構(gòu)的影響

跨流域調(diào)水和糧食區(qū)域貿(mào)易是分別通過實體水調(diào)運和虛擬水進口的方式依賴其他地區(qū)水資源，達到緩解本地區(qū)水資源短缺的目的。因此，文章用區(qū)域總用水量中的跨流域調(diào)水量以及糧食虛擬水凈進口量作為衡量區(qū)域“實體水-虛擬水”流動影響的基本元素。

“實體水-虛擬水”流動對區(qū)域水資源約束的改變，可以通過糧食生產(chǎn)誘致性技術(shù)進步程度、資源稟賦約束以及糧食區(qū)域生產(chǎn)與分工強度三個途徑影響糧食種植結(jié)構(gòu)，且對糧食輸出區(qū)和糧食輸入?yún)^(qū)的影響具有異質(zhì)性（圖1）。對于糧食輸出區(qū)來說，跨流域調(diào)水量減少和糧食虛擬水凈出口量增加，會加強水資源稟賦約束，對糧食生產(chǎn)產(chǎn)生一定的限制作用，從而糧食種植比例降低。其次，根據(jù)誘致性技術(shù)變遷理論，水資源稀缺性增加會誘致節(jié)約水資源的技術(shù)進步，如抗旱種子、節(jié)水灌溉技術(shù)等。由于技術(shù)選擇的非均衡性，主要糧食作物的技術(shù)創(chuàng)新率明顯更快［43-44］，使得水資源稀缺性增加情況下，糧食輸出區(qū)糧食作物相對于經(jīng)濟作物技術(shù)進步更快，從而糧食種植比例提升。第三，根據(jù)產(chǎn)業(yè)布局理論，區(qū)域之間貿(mào)易自由化可以整合國內(nèi)市場，強化區(qū)域?qū)I(yè)化程度，逐漸形成產(chǎn)業(yè)集聚，充分發(fā)揮比較優(yōu)勢，促使各地區(qū)產(chǎn)業(yè)結(jié)構(gòu)差異化發(fā)展［45-46］。糧食輸出區(qū)糧食虛擬水凈出口量增加，會進一步強化糧食區(qū)域生產(chǎn)與分工格局，提高糧食輸出區(qū)糧食種植比例。

圖1 “實體水-虛擬水”流動對不同區(qū)域糧食種植結(jié)構(gòu)的影響機制

對于糧食輸入?yún)^(qū)來說，跨流域調(diào)水量減少同樣會加強水資源稟賦約束，一方面對糧食生產(chǎn)具有限制作用，降低糧食種植比例；另一方面促進糧食作物誘致性技術(shù)變遷，提高糧食種植比例。而糧食輸入?yún)^(qū)的糧食虛擬水凈進口量增加，會進一步強化糧食區(qū)域生產(chǎn)與分工格局，使得糧食輸入?yún)^(qū)糧食種植比例進一步下降。

由此可知，“實體水-虛擬水”流動對不同區(qū)域糧食種植比例的影響取決于誘致性技術(shù)進步程度、資源稟賦約束以及糧食區(qū)域生產(chǎn)與分工強度的影響大小。

從糧食內(nèi)部結(jié)構(gòu)來看，“實體水-虛擬水”流動通過影響不同糧食作物的邊際用水成本，進而影響糧食內(nèi)部結(jié)構(gòu)?！皩嶓w水-虛擬水”流動影響降低，加強本地糧食生產(chǎn)水資源約束，提高邊際用水成本。邊際用水成本增加一方面會降低高耗水作物種植比例，提高低耗水作物種植比例；另一方面，會促進節(jié)水技術(shù)的創(chuàng)新和推廣，反而會提高高耗水作物種植比例，即“灌溉效率悖論”［47-48］。對不同糧食作物種植結(jié)構(gòu)的差異性影響取決于這兩種作用的大小。

由于糧食區(qū)域貿(mào)易關(guān)系，“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)的影響存在空間溢出效應(yīng)，且對于糧食輸出區(qū)與糧食輸入?yún)^(qū)的影響具有異質(zhì)性。糧食輸入?yún)^(qū)“實體水-虛擬水”流動影響提高，一方面，加劇了糧食輸出區(qū)的水資源壓力，制約糧食生產(chǎn)；另一方面，可能會引致糧食輸出區(qū)糧食生產(chǎn)技術(shù)進步，提高糧食種植比例；此外，糧食區(qū)域生產(chǎn)與分工的強化使得糧食輸出區(qū)糧食種植比例提高。糧食輸入?yún)^(qū)“實體水-虛擬水”流動對糧食輸出區(qū)糧食種植比例的空間溢出效應(yīng)取決于這三者的大小。糧食輸出區(qū)“實體水-虛擬水”流動影響降低，一方面，由于糧食區(qū)域生產(chǎn)與分工的強化會降低糧食輸入?yún)^(qū)糧食種植比例；另一方面，當糧食輸出區(qū)水資源壓力較大，對糧食生產(chǎn)的限制作用大于誘致性技術(shù)進步時，可能會促使糧食輸入?yún)^(qū)提高糧食種植比例，從而滿足本地糧食需求。糧食輸出區(qū)“實體水-虛擬水”流動對糧食輸入?yún)^(qū)糧食種植比例的空間溢出效應(yīng)取決于這兩者的大小。

3 納入“實體水-虛擬水”流動影響的研究方法與數(shù)據(jù)說明

根據(jù)理論分析，“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)的影響具有空間溢出效應(yīng)，因此采用空間計量模型進行分析。納入“實體水-虛擬水”流動影響的研究思路與方法如圖2所示。首先對糧食種植結(jié)構(gòu)和“實體水-虛擬水”流動影響進行空間自相關(guān)檢驗，初步判斷空間相關(guān)性。進一步地，同時考慮空間滯后因變量和自變量，以及糧食種植結(jié)構(gòu)可能存在的“時間慣性”，構(gòu)建更一般化的動態(tài)空間杜賓模型分析“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)的影響。然后，從糧食區(qū)域貿(mào)易關(guān)系說明空間權(quán)重矩陣的設(shè)定。最后，對核心自變量跨流域調(diào)水、糧食區(qū)域貿(mào)易以及“實體水-虛擬水”流動影響指標構(gòu)建進行說明。

圖2 納入“實體水-虛擬水”流動影響的研究思路與方法

3.1 空間自相關(guān)檢驗

在決定使用普通面板模型或者空間面板模型之前，需要對糧食種植結(jié)構(gòu)的空間相關(guān)性進行檢驗，采用Moran’sI指數(shù)進行檢驗，計算公式如下：

其中：Yi為樣本觀測值（即為不同糧食作物播種面積占比及“實體水-虛擬水”流動影響指標）為樣本平均值；S2為樣本方差；Wij為空間權(quán)重矩陣W的(i，j)元素。其中，I∈[-1，1]。當I＞0時，說明存在正向的空間相關(guān)性，即在空間權(quán)重矩陣設(shè)定下，糧食種植比例、三大主糧種植比例以及“實體水-虛擬水”流動影響呈現(xiàn)高值與高值相鄰、低值與低值相鄰的空間分布；當I＜0時，說明存在負向的空間相關(guān)性，即糧食種植比例、三大主糧種植比例以及“實體水-虛擬水”流動影響呈現(xiàn)高值與低值相鄰的空間分布；當I=0時，說明空間分布是隨機的，不存在空間相關(guān)性。

3.2 空間模型設(shè)定

確定了空間相關(guān)性存在之后，考慮空間因素建立空間計量模型進行分析。根據(jù)理論分析，選擇更一般化的空間杜賓模型，同時考慮空間滯后因變量和自變量；此外，糧食種植結(jié)構(gòu)可能存在“時間慣性”，因此需要考慮時間維度上的動態(tài)性。最終，選擇利用動態(tài)空間杜賓模型進行實證分析。與靜態(tài)空間面板模型相比，動態(tài)空間面板模型既考慮了糧食種植結(jié)構(gòu)的動態(tài)變化和空間溢出效應(yīng)［49-50］，又可以避免“雞蛋相生”的內(nèi)生性問題［51］，從而使得模型估計結(jié)果更加準確和可靠。首先，考察“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)的影響，模型設(shè)定為：

進一步地，考察跨流域調(diào)水和糧食區(qū)域貿(mào)易對糧食種植結(jié)構(gòu)的影響，模型設(shè)定為：

其中：Yit為i省份第t年某種糧食作物播種面積占比，δit為i省份第t年“實體水-虛擬水”流動影響，δtran_it為i省份第t年跨流域調(diào)水影響，δvw_it為i省份第t年糧食區(qū)域貿(mào)易影響，Zit為其他影響糧食種植結(jié)構(gòu)的控制變量。W為空間權(quán)重矩陣；ρ和ρ'為空間自回歸系數(shù)，表示相關(guān)省份糧食種植結(jié)構(gòu)對本省糧食種植結(jié)構(gòu)的影響；τ和τ'表示種植結(jié)構(gòu)調(diào)整慣性，表示上一期糧食種植結(jié)構(gòu)對當期糧食種植結(jié)構(gòu)的影響。動態(tài)空間杜賓模型同時考慮了糧食種植結(jié)構(gòu)的時間和空間滯后性，不僅可以得到“實體水-虛擬水”流動影響的直接效應(yīng)和間接效應(yīng)，還可以得到兩種效應(yīng)的短期效應(yīng)和長期效應(yīng)。

直接效應(yīng)指一個省份中的解釋變量發(fā)生一個單位的變化對本省份糧食種植結(jié)構(gòu)的影響，等于模型系數(shù)與反饋效應(yīng)之和，即其他省份的糧食種植結(jié)構(gòu)反過來影響該省份的糧食種植結(jié)構(gòu)。間接效應(yīng)指相關(guān)省份的解釋變量發(fā)生一個單位的變化對本省份糧食種植結(jié)構(gòu)的影響。直接效應(yīng)與間接效應(yīng)之和為總效應(yīng)，可解釋為所有省份解釋變量發(fā)生一個單位的變化對本省份糧食種植結(jié)構(gòu)的影響。以式（2）為例，直接效應(yīng)和間接效應(yīng)表示為：

3.3 空間權(quán)重矩陣設(shè)定

空間權(quán)重矩陣反映了空間單元之間的相互依賴性與關(guān)聯(lián)程度，合理恰當?shù)剡x擇空間權(quán)重矩陣對于空間計量模型的估計結(jié)果和解釋力至關(guān)重要［52］。根據(jù)研究問題與視角不同，已有研究發(fā)現(xiàn)相對于地理鄰接、地理距離矩陣，根據(jù)經(jīng)濟距離、貿(mào)易距離等構(gòu)建的空間權(quán)重矩陣有時能更好地反映空間單元之間的關(guān)系。如林光平等［53］使用地區(qū)間人均GDP差額構(gòu)建經(jīng)濟權(quán)重矩陣，發(fā)現(xiàn)能更好地擬合中國地區(qū)經(jīng)濟的發(fā)展狀況。徐春華等［54］采用貿(mào)易份額指標構(gòu)建貿(mào)易權(quán)重矩陣，更有力地解釋了FDI和政府消費對CO2排放的影響。此外，兩個空間單元之間的相互影響在強度上可能是不對稱的，如經(jīng)濟發(fā)展水平較高的空間單元會對經(jīng)濟發(fā)展水平較低地區(qū)產(chǎn)生更強的空間影響與輻射作用，合理地反映這種非對稱性十分必要［55］。楊鑫［39］根據(jù)糧食省際貿(mào)易量分別構(gòu)建了糧食省際調(diào)出矩陣和糧食省際調(diào)入矩陣，檢驗了糧食生產(chǎn)用水效率通過糧食省際貿(mào)易對糧食調(diào)入?yún)^(qū)和糧食調(diào)出區(qū)經(jīng)濟水平影響的非對稱空間溢出效應(yīng)。

根據(jù)理論分析，“實體水-虛擬水”流動對于糧食輸出區(qū)（包括山西、內(nèi)蒙古、吉林、黑龍江、江蘇、安徽、江西、河南、陜西、甘肅、寧夏和新疆12個省份）與糧食輸入?yún)^(qū)（包括北京、天津、河北、遼寧、上海、浙江、福建、山東、湖北、湖南、廣東、廣西、海南、重慶、四川、貴州、云南、西藏和青海19個省份）糧食種植結(jié)構(gòu)的影響路徑存在差異，且兩者之間的空間溢出效應(yīng)具有非對稱性。因此，分別從生產(chǎn)者與消費者兩個角度核算省際糧食貿(mào)易的虛擬水貿(mào)易量，并構(gòu)建兩個空間權(quán)重矩陣，以捕獲上述區(qū)域異質(zhì)性和空間溢出效應(yīng)的非對稱關(guān)系。以生產(chǎn)者角度計算的糧食虛擬水貿(mào)易量，指的是糧食輸出區(qū)由于糧食貿(mào)易真實消耗的虛擬水量；以消費者角度計算的糧食虛擬水貿(mào)易量，指的是假如糧食輸入?yún)^(qū)本地生產(chǎn)輸入的糧食所要消耗的虛擬水量。可以看出，糧食輸入?yún)^(qū)對糧食輸出區(qū)糧食種植結(jié)構(gòu)影響的空間溢出效應(yīng)主要通過生產(chǎn)者角度糧食虛擬水貿(mào)易作用（糧食輸出區(qū)為滿足糧食輸入?yún)^(qū)的糧食需求而實際消耗的糧食虛擬水），而糧食輸出區(qū)對糧食輸入?yún)^(qū)糧食種植結(jié)構(gòu)影響主要通過消費者角度糧食虛擬水貿(mào)易作用（糧食輸入?yún)^(qū)因為從糧食輸出區(qū)進口糧食而節(jié)約的本地糧食虛擬水）。

3.3.1 糧食省際貿(mào)易量估計方法

參考丁雪麗等［16］、楊鑫［39］，分為口糧消費、飼料消費、工業(yè)消費、種子用糧和糧食損耗五大類對糧食消費量進行估計。口糧消費包括家庭糧食消費和家庭外糧食消費，由各省份城鎮(zhèn)居民人均糧食消費量和農(nóng)村居民人均糧食消費量分別乘以城鎮(zhèn)人口和農(nóng)村人口相加得到家庭糧食消費，城鎮(zhèn)居民在外糧食消費比例設(shè)為15%，農(nóng)村設(shè)為10%；飼料消費根據(jù)各省畜禽產(chǎn)品產(chǎn)量按照飼料轉(zhuǎn)換率進行轉(zhuǎn)換，其中豬肉1∶2.5，牛肉1∶0.7，羊肉1∶0.5，禽肉1∶1.8，禽蛋1∶1.7，水產(chǎn)品1∶1；工業(yè)消費按照一定的折算比例對工業(yè)用糧主要用途酒精、白酒、啤酒、味精和其他用途的用糧情況進行估算，其中酒精1∶3，白酒1∶2.3；啤酒1∶0.172，味精1∶2.4，其他工業(yè)用糧按占比25%計算；種子消費依據(jù)各省分品種糧食播種面積，并按照水稻2.78 kg/畝、小麥15 kg/畝、玉米15 kg/畝和大豆5.26 kg/畝進行折算，薯類種子用量按產(chǎn)量的10%計算；糧食損耗按糧食消費的5%計算。根據(jù)2004—2020年各省份糧食產(chǎn)量與消費量均值，構(gòu)建省際間糧食貿(mào)易線性規(guī)劃模型，以運輸成本最小化為約束目標，估計省際間糧食貿(mào)易的平均情況。具體線性規(guī)劃模型為：

其中：Z為糧食省際貿(mào)易的運輸總費用；i為糧食輸入省份（總數(shù)為N），j為糧食輸出省份（總數(shù)為31-N），tradeij為省份j向省份i輸出的糧食貿(mào)易量；cij為省份j到省份i的單位糧食運輸成本，rij為省份j到省份i的單位糧食鐵路運輸成本（鐵路運輸成本計算方式為省會間鐵路距離乘以鐵路貨運價格，水路運輸成本計算方式為碼頭間航線距離乘以平均貨運成本。鐵路運價來自2015年國家發(fā)展改革委發(fā)布的《各類貨物鐵路運輸基準運價率表》，糧食運價號為4），wij為省份j到省份i的單位糧食水路運輸成本（水路運價價格變化頻繁，為增加計算可操作性，將水運價格固定為0.041 5元/（km·t））；Si為省份i的糧食短缺量；Rj為省份j的糧食過剩量。通過線性規(guī)劃求解，即得到最優(yōu)條件下的糧食省際貿(mào)易量tradeij。

將糧食省際貿(mào)易量tradeij作為糧食虛擬水貿(mào)易矩陣的基本元素?？臻g權(quán)重矩陣具體設(shè)置如下。

3.3.2 生產(chǎn)者角度糧食虛擬水貿(mào)易矩陣

根據(jù)生產(chǎn)者角度糧食虛擬水貿(mào)易矩陣計算的直接效應(yīng)是糧食輸出區(qū)的本地效應(yīng)，間接效應(yīng)是糧食輸入?yún)^(qū)對糧食輸出區(qū)的空間溢出效應(yīng)，矩陣設(shè)置為：

其中：i為糧食輸出區(qū)，j為糧食輸入?yún)^(qū)，tradeij為2004—2020年i省份向j省份輸出的糧食貿(mào)易量均值（萬t），vwcsi為糧食輸出區(qū)2004—2020年單位糧食生產(chǎn)水足跡均值（m3/t）（利用2000年國際水足跡網(wǎng)絡(luò)數(shù)據(jù)庫各省份品種糧食生產(chǎn)單位水足跡數(shù)據(jù)（WaterStat）（http：//waterfootprint.org/en/resources/waterstat），采用“速達法”對2004—2020年份省單位糧食生產(chǎn)水足跡進行推算，由小麥、玉米、水稻、高粱、大麥、谷子、大豆、馬鈴薯及紅薯9種作物加權(quán)計算得出?！八龠_法”是基于作物單位生產(chǎn)水足跡與單產(chǎn)水平呈負線性相關(guān)關(guān)系的原理，以特定年份作物生產(chǎn)水足跡數(shù)據(jù)庫推算多年份作物生產(chǎn)水足跡量級的快速方法，計算公式為：VWC(t)=[VWC(T)×y(T)]/y(t)。其中，VWC(t)與VWC(T)（m3/t）分別為t年與T年單位糧食生產(chǎn)水足跡，y(t)與y(T)（t/ha）分別為t年與T年的單位面積糧食產(chǎn)量。Tuninetti等［58］通過全球尺度不確定分析驗證該方法誤差在10%以內(nèi)。因具有計算成本低、結(jié)果可靠性高的特點，該方法在國內(nèi)外已得到廣泛應(yīng)用［59-61］）。

3.3.3 消費者角度糧食虛擬水貿(mào)易矩陣

根據(jù)消費者角度糧食虛擬水貿(mào)易矩陣計算的直接效應(yīng)是糧食輸入?yún)^(qū)的本地效應(yīng)，間接效應(yīng)是糧食輸出區(qū)對糧食輸入?yún)^(qū)的空間溢出效應(yīng)，矩陣設(shè)置為：

其中：i為糧食輸入?yún)^(qū)，j為糧食輸出區(qū)，tradeij為2004—2020年i省份從j省份輸入的糧食貿(mào)易量均值（萬t），vwcdi為糧食輸入?yún)^(qū)2004—2020年單位糧食生產(chǎn)水足跡均值（m3/t）。

3.4 “實體水-虛擬水”流動、跨流域調(diào)水以及糧食區(qū)域貿(mào)易影響的指標構(gòu)建

參考Zhao等［62］、Ye等［63］，從水資源壓力指數(shù)出發(fā)，構(gòu)建“實體水-虛擬水”流動、跨流域調(diào)水以及糧食區(qū)域貿(mào)易影響指標。水資源壓力指數(shù)（WSI）是指由于人類活動，從本地可用水源取水引起的水資源壓力。不考慮跨流域調(diào)水的各省份水資源壓力指數(shù)（WSI）為：

假設(shè)沒有跨流域調(diào)水以及糧食區(qū)域貿(mào)易，即各省份需要利用本地水資源滿足其糧食消費需求，則各省份的水資源壓力指數(shù)（WSI*）為：

則兩者之間的差距，為“實體水-虛擬水”流動產(chǎn)生的影響：

進一步地，可分解為跨流域調(diào)水影響以及糧食區(qū)域貿(mào)易影響：

其中：δ為“實體水-虛擬水”流動影響；δtran為跨流域調(diào)水影響；δvwt為糧食區(qū)域貿(mào)易影響；WU為用水總量（億m3）；Wtransfer為跨流域調(diào)水量（億m3）；VWim為糧食虛擬水凈進口量（億m3），C為糧食消費量（t），P為糧食產(chǎn)量（t），VWC為單位糧食生產(chǎn)水足跡（億m3/t）；Q為水資源稟賦量（億m3），由當?shù)厮Y源總量和入境水資源量（不包括省際跨流域調(diào)水量）組成，考慮生態(tài)用水，定義水資源稟賦量的20%為可利用水資源量［26］；B為調(diào)整系數(shù)（消除豐水和枯水年影響），其數(shù)值為各省份歷年水資源稟賦量除以對應(yīng)所有年份的平均水資源稟賦量［18］。

3.5 變量選取

3.5.1 被解釋變量

參考楊進等［64］、錢龍等［65］，使用各類作物占作物總種植面積的比例來表示糧食種植結(jié)構(gòu)的變化。具體包括兩個層面的種植結(jié)構(gòu)，首先是糧食作物種植面積占農(nóng)作物總種植面積的比例，反映糧食生產(chǎn)在整個農(nóng)作物生產(chǎn)當中的結(jié)構(gòu)問題。第二是糧食內(nèi)部生產(chǎn)結(jié)構(gòu)的調(diào)整。對中國而言，糧食主要指三大主糧，即水稻、小麥和玉米，因而使用三類主糧種植面積占總種植面積比例變化來表示糧食內(nèi)部生產(chǎn)結(jié)構(gòu)調(diào)整。

3.5.2 核心解釋變量

核心解釋變量即為根據(jù)上文公式測算出的“實體水-虛擬水”流動影響、跨流域調(diào)水影響以及糧食區(qū)域貿(mào)易影響。

對于糧食輸出區(qū)來說，“實體水-虛擬水”流動對水資源壓力的影響取決于跨流域調(diào)水和糧食區(qū)域貿(mào)易的影響大小，當影響為正時，糧食區(qū)域貿(mào)易發(fā)揮作用越大，反之，則跨流域調(diào)水影響更大。對于糧食輸入?yún)^(qū)來說，“實體水-虛擬水”流動影響為正，跨流域調(diào)水量和糧食虛擬水輸入量越大，對緩解水資源壓力的作用越大。“實體水-虛擬水”流動對全國層面水資源壓力的影響取決于其對糧食輸出區(qū)和糧食輸入?yún)^(qū)的影響大?。ㄋY源壓力分級參考Hoekstra等［66］劃分：低水資源壓力（水資源壓力指數(shù)WSI＜1.0）、中等水資源壓力（1.0≤WSI＜1.5）、高水資源壓力（1.5≤WSI＜2.0）和嚴重水資源壓力（WSI＞2.0）。任一區(qū)域的水資源壓力指數(shù)按該區(qū)域各省級行政區(qū)的國土面積加權(quán)后進行計算［26］）。

由圖3可知，從糧食輸出區(qū)來看，2004—2020年“實體水-虛擬水”流動影響均為負值，且持續(xù)降低，說明糧食輸出區(qū)糧食虛擬水輸出大于跨流域調(diào)水的補充，且程度不斷加大，糧食輸出區(qū)水資源壓力不斷加劇。從糧食輸入?yún)^(qū)來看，2004—2020年“實體水-虛擬水”流動影響均為正值，說明“實體水-虛擬水”流動格局緩解了糧食輸入?yún)^(qū)的水資源壓力，糧食輸入?yún)^(qū)水資源壓力呈下降趨勢。從全國來看，2004—2006年“實體水-虛擬水”流動影響為正，2007—2020年“實體水-虛擬水”流動影響變?yōu)樨撝?，說明隨著糧食生產(chǎn)不斷向中西部和北方轉(zhuǎn)移和集中，糧食輸出區(qū)水資源壓力遠大于糧食輸入?yún)^(qū)，跨流域調(diào)水作用有限，全國水資源壓力總體上呈上升趨勢。

圖3 2004—2020年“實體水-虛擬水”流動影響的變化趨勢

3.5.3 控制變量

參考楊鑫等［18］、楊宗輝等［67］以及楊進等［64］，從自然基礎(chǔ)條件、非農(nóng)就業(yè)機會、糧食及替代作物生產(chǎn)效益以及糧食生產(chǎn)政策四個方面引入控制變量。包括全省平均年降水量、全省平均年氣溫、有效灌溉面積占耕地面積比重，農(nóng)村居民消費指數(shù)平減后的農(nóng)村居民人均工資性收入，分品種生產(chǎn)價格指數(shù)平減后的水稻、小麥、玉米、大豆、油菜的畝均凈利潤（根據(jù)《農(nóng)產(chǎn)品成本收益資料匯編》，水稻畝均凈利潤主要考察23個水稻生產(chǎn)省份，小麥畝均凈利潤考察15個小麥生產(chǎn)省份，玉米畝均凈利潤考察20個玉米生產(chǎn)省份。糧食作物之間的替代關(guān)系選擇，參考林大燕等［68］，主要根據(jù)主產(chǎn)區(qū)糧食作物的耕作制度和種植規(guī)模，分別選取水稻、玉米、小麥的不同替代作物。水稻在黑龍江的替代作物主要是大豆；小麥在江蘇、安徽、湖北的替代作物主要是油菜；玉米在黑龍江、吉林、遼寧、內(nèi)蒙古的替代作物主要是大豆），小麥最低收購價格、水稻最低收購價格以及玉米臨時收儲價格等變量。小麥和水稻最低收購價政策的實施時間分別為2006—2020年和2004—2020年，實施省份參照歷年《小麥和稻谷最低收購價執(zhí)行預(yù)案》第一條；玉米臨時收儲政策的實施時間為2008—2015年，具體實施省份為內(nèi)蒙古、遼寧、吉林和黑龍江。

3.6 數(shù)據(jù)來源與說明

研究對象為31個省份（未涉及香港、澳門、臺灣）的糧食整體及水稻、小麥、玉米三大主糧作物，研究時間段為2004—2020年。分品種糧食播種面積、產(chǎn)量、有效灌溉面積、農(nóng)村居民人均工資性收入、農(nóng)產(chǎn)品生產(chǎn)價格指數(shù)以及農(nóng)村居民消費價格指數(shù)來自國家統(tǒng)計局分省年度數(shù)據(jù)庫。水資源總量、用水總量、跨流域調(diào)水量、入境水資源量數(shù)據(jù)來自歷年各省《水資源公報》。降水量、溫度數(shù)據(jù)來自歷年《中國統(tǒng)計年鑒》。畝均凈利潤數(shù)據(jù)來自歷年《農(nóng)產(chǎn)品成本收益資料匯編》。估計糧食消費量的基礎(chǔ)數(shù)據(jù)主要來自歷年《中國統(tǒng)計年鑒》、各省統(tǒng)計年鑒、《中國輕工業(yè)年鑒》以及《中國食品工業(yè)年鑒》。變量的詳細說明和統(tǒng)計性指標見表1。

表1 變量說明與描述性統(tǒng)計

4 實證結(jié)果與分析

4.1 糧食種植結(jié)構(gòu)與“實體水-虛擬水”流動影響的空間自相關(guān)檢驗

通過式（1）計算得到2004—2020年中國糧食種植結(jié)構(gòu)及“實體水-虛擬水”流動影響的Moran’sI指數(shù)。由表2可得，2004—2020年糧食種植比例具有顯著的負向空間相關(guān)性，說明糧食輸出區(qū)與糧食輸入?yún)^(qū)的糧食種植比例呈現(xiàn)高值與低值“相鄰”的空間分布；玉米種植比例具有顯著的正向空間相關(guān)性，說明糧食輸出區(qū)與糧食輸入?yún)^(qū)的玉米種植比例呈現(xiàn)高值與高值、低值與低值“相鄰”的空間分布；“實體水-虛擬水”影響、糧食區(qū)域貿(mào)易影響具有顯著的負向空間相關(guān)性，說明糧食輸出區(qū)與糧食輸入?yún)^(qū)的“實體水-虛擬水”流動影響、糧食區(qū)域貿(mào)易影響呈現(xiàn)高值與低值“相鄰”的空間分布，跨流域調(diào)水影響的空間相關(guān)性不明顯。

表2 糧食種植結(jié)構(gòu)與“實體水-虛擬水”流動的空間自相關(guān)檢驗

在空間自相關(guān)檢驗的基礎(chǔ)上，進一步根據(jù)空間面板數(shù)據(jù)Moran檢驗、LM檢驗以及LR檢驗結(jié)果，并結(jié)合擬合優(yōu)度等指標，最終選擇個體固定效應(yīng)的空間杜賓模型進行實證分析，糧食輸出區(qū)和糧食輸入?yún)^(qū)模型估計結(jié)果見表3和表4。進一步通過效應(yīng)分解，分析“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)的本地效應(yīng)和空間溢出效應(yīng)，結(jié)果整理見表5和表6。

表3 糧食輸出區(qū)“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的回歸結(jié)果

表4 糧食輸入?yún)^(qū)“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的回歸結(jié)果

4.2 “實體水-虛擬水”流動對糧食輸出區(qū)糧食種植結(jié)構(gòu)的本地效應(yīng)

表5空間權(quán)重矩陣為Wpij部分報告了“實體水-虛擬水”流動對糧食輸出區(qū)糧食種植結(jié)構(gòu)的短期和長期直接效應(yīng)，即本地效應(yīng)。

對于糧食輸出區(qū)來說，“實體水-虛擬水”流動影響是降低的，“實體水-虛擬水”流動影響降低1%，糧食種植比例在短期和長期分別提高0.013%和0.182%。說明“實體水-虛擬水”流動對于糧食輸出區(qū)水資源壓力的增加，對糧食生產(chǎn)誘致性技術(shù)變遷的促進作用與區(qū)域生產(chǎn)分工的強化作用大于限制作用。進一步來看，跨流域調(diào)水量降低和糧食虛擬水出口增加對糧食輸出區(qū)的糧食生產(chǎn)促進作用均大于限制作用。

從糧食內(nèi)部結(jié)構(gòu)來看，“實體水-虛擬水”流動影響對糧食輸出區(qū)水稻和小麥種植比例不顯著，主要表現(xiàn)為對玉米種植比例的影響?！皩嶓w水-虛擬水”流動影響降低1%，玉米種植比例在短期和長期分別下降0.080%和0.157%，且糧食區(qū)域貿(mào)易的作用大于跨流域調(diào)水的影響。即當用水邊際成本提高時，抗旱能力較強、相較于水稻和小麥屬于低耗水作物的玉米種植比例反而會下降，說明節(jié)水灌溉技術(shù)進步促使糧食內(nèi)部結(jié)構(gòu)向更加耗水的方向發(fā)展，證明存在“灌溉效率悖論”。

4.3 “實體水-虛擬水”流動對糧食輸入?yún)^(qū)糧食種植結(jié)構(gòu)的本地效應(yīng)

表5空間權(quán)重矩陣為Wcij部分報告了“實體水-虛擬水”流動對糧食輸入?yún)^(qū)糧食種植結(jié)構(gòu)的短期和長期直接效應(yīng)，即本地效應(yīng)。

表5 “實體水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的本地效應(yīng)

對于糧食輸入?yún)^(qū)來說，“實體水-虛擬水”流動影響是提高的，“實體水-虛擬水”流動影響提高1%，糧食種植比例在短期和長期分別下降0.016%和0.129%。一方面，糧食虛擬水進口增加強化區(qū)域糧食生產(chǎn)與分工，使得糧食輸入?yún)^(qū)糧食種植比例下降。另一方面，跨流域調(diào)水量和糧食虛擬水進口增加放松了糧食輸入?yún)^(qū)的水資源約束，而由于種植經(jīng)濟作物的收益高于糧食作物，經(jīng)濟作物播種面積提高得更多，從而導(dǎo)致糧食種植比例下降。

從糧食內(nèi)部結(jié)構(gòu)來看，“實體水-虛擬水”流動影響主要表現(xiàn)為對水稻和玉米種植比例的影響?！皩嶓w水-虛擬水”流動影響提高1%，水稻種植比例短期提高0.014%，長期則下降0.025%；玉米種植比例短期提高0.012%，長期則下降0.029%。即當“實體水-虛擬水”流動影響提高時，短期糧食輸入?yún)^(qū)高耗水作物水稻種植比例的提高幅度大于低耗水作物玉米種植比例，長期來看低耗水作物玉米種植比例的下降幅度更大，說明無論從短期還是長期來看，糧食輸入?yún)^(qū)糧食內(nèi)部結(jié)構(gòu)均向耗水方向發(fā)展。

4.4 “實體水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的空間溢出效應(yīng)

在分析了“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)本地效應(yīng)的基礎(chǔ)上，進一步分析其對糧食輸出區(qū)和糧食輸入?yún)^(qū)糧食種植結(jié)構(gòu)的空間溢出效應(yīng)。

4.4.1 糧食輸入?yún)^(qū)對糧食輸出區(qū)的空間溢出效應(yīng)

表6空間權(quán)重矩陣為Wpij部分報告了糧食輸入?yún)^(qū)對糧食輸出區(qū)的空間溢出效應(yīng)。從糧食整體結(jié)構(gòu)來看，糧食輸入?yún)^(qū)“實體水-虛擬水”流動影響提高1%，糧食輸出區(qū)的糧食種植比例在短期和長期分別下降0.012%和0.213%。說明隨著糧食輸入?yún)^(qū)“實體水-虛擬水”流動影響的提高，糧食輸出區(qū)水資源壓力增加對糧食生產(chǎn)的制約作用要大于誘致性技術(shù)進步程度以及糧食區(qū)域生產(chǎn)與分工強度帶來的促進作用，使得糧食輸出區(qū)糧食種植比例下降。進一步來看，糧食輸入?yún)^(qū)對糧食輸出區(qū)的空間溢出效應(yīng)主要是由于糧食區(qū)域貿(mào)易，跨流域調(diào)水的間接效應(yīng)不顯著。

從糧食內(nèi)部結(jié)構(gòu)來看，糧食輸入?yún)^(qū)“實體水-虛擬水”流動主要表現(xiàn)為對糧食輸出區(qū)玉米種植比例的空間溢出效應(yīng)。“實體水-虛擬水”流動影響提高1%，糧食輸出區(qū)的玉米種植比例在短期和長期分別下降0.142%和0.207%。說明糧食輸入?yún)^(qū)加劇了糧食輸出區(qū)的“灌溉效率悖論”情況。

4.4.2 糧食輸出區(qū)對糧食輸入?yún)^(qū)的空間溢出效應(yīng)

表6空間權(quán)重矩陣為Wcij部分報告了糧食輸出區(qū)對糧食輸入?yún)^(qū)的空間溢出效應(yīng)。從糧食整體結(jié)構(gòu)來看，糧食輸出區(qū)“實體水-虛擬水”流動對糧食輸入?yún)^(qū)糧食種植比例沒有影響。從糧食內(nèi)部結(jié)構(gòu)來看，糧食輸出區(qū)“實體水-虛擬水”流動對糧食輸入?yún)^(qū)的影響，主要表現(xiàn)為對水稻和玉米種植比例的空間溢出效應(yīng)。糧食輸出區(qū)“實體水-虛擬水”流動影響降低1%，短期糧食輸入?yún)^(qū)水稻和玉米種植比例分別降低0.051%和0.074%，長期分別提高0.017%和0.140%，即短期高耗水作物水稻種植比例下降更多，長期低耗水作物玉米種植比例提高更多，說明無論長期還是短期，糧食輸出區(qū)“實體水-虛擬水”流動影響的溢出效應(yīng)使得糧食輸入?yún)^(qū)糧食內(nèi)部結(jié)構(gòu)向節(jié)水方向發(fā)展。

表6 “實體水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的空間溢出效應(yīng)

5 穩(wěn)健性檢驗

在考察跨流域調(diào)水和糧食區(qū)域貿(mào)易對不同糧食作物種植結(jié)構(gòu)的影響時，空間關(guān)系不同設(shè)定可能會導(dǎo)致估計結(jié)果的差異。糧食種植結(jié)構(gòu)受地理位置、自然資源稟賦等影響較大，因此在糧食虛擬水貿(mào)易空間權(quán)重矩陣的基礎(chǔ)上，考慮地理距離因素，采用地理距離加權(quán)的糧食虛擬水貿(mào)易矩陣①穩(wěn)健性檢驗的生產(chǎn)者角度糧食虛擬水貿(mào)易矩陣設(shè)置為，其中i為糧食輸出區(qū)，j為糧食輸入?yún)^(qū)，tradeij為2004—2020年i省份向j省份輸出的糧食貿(mào)易量均值（萬t），vwcsi為糧食輸出區(qū)2004—2020年單位糧食生產(chǎn)水足跡均值（m3/t）。穩(wěn)健性檢驗的消費者角度糧食虛擬水貿(mào)易矩陣設(shè)置為，其中i為糧食輸入?yún)^(qū)，j為糧食輸出區(qū)，tradeij為2004—2020年i省份從j省份輸入的糧食貿(mào)易量均值（萬t），vwcdi為糧食輸入?yún)^(qū)2004—2020年單位糧食生產(chǎn)水足跡均值（m3/t）。d為根據(jù)經(jīng)緯度計算的省會城市之間的地表距離（km）。，對“實體水-虛擬水”流動對糧食種植結(jié)構(gòu)的影響及空間溢出效應(yīng)再次進行估計和檢驗?！皩嶓w水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的本地效應(yīng)和空間溢出效應(yīng)分別見表7和表8。本地效應(yīng)的符號和顯著性與采用不考慮地理距離的糧食虛擬水貿(mào)易矩陣的估計結(jié)果基本一致。

表7 “實體水-虛擬水”流動對糧食種植結(jié)構(gòu)影響的本地效應(yīng)穩(wěn)健性檢驗

表8 “實體水-虛擬水”流動對糧食種植結(jié)構(gòu)空間溢出效應(yīng)的穩(wěn)健性檢驗

從糧食輸入?yún)^(qū)對糧食輸出區(qū)的空間溢出效應(yīng)來看，糧食輸入?yún)^(qū)“實體水-虛擬水”流動影響提高，對糧食輸出區(qū)糧食種植比例具有負向作用，但短期間接效應(yīng)不再顯著，說明長時間的糧食省際貿(mào)易，即糧食區(qū)域生產(chǎn)與分工的穩(wěn)定，使得糧食輸入?yún)^(qū)“實體水-虛擬水”流動影響對糧食輸出區(qū)糧食種植比例產(chǎn)生負的空間溢出效應(yīng)。糧食輸出區(qū)對糧食輸入?yún)^(qū)糧食種植比例的影響不顯著，與原始估計結(jié)果相似。從糧食內(nèi)部結(jié)構(gòu)來看，跨流域調(diào)水影響對糧食輸出區(qū)水稻種植比例和糧食輸入?yún)^(qū)玉米種植比例的間接效應(yīng)不再顯著，也進一步說明“實體水-虛擬水”流動影響對糧食內(nèi)部結(jié)構(gòu)的空間溢出效應(yīng)主要是由于糧食區(qū)域貿(mào)易的作用。其他影響方向和顯著性與原始回歸相似，說明主要研究結(jié)論具有一定的穩(wěn)健性。

6 結(jié)論與啟示

文章基于2004—2020年中國31個省份的面板數(shù)據(jù)，從“實體水-虛擬水”流動視角測算出全國和糧食輸出區(qū)、糧食輸入?yún)^(qū)跨流域調(diào)水和糧食區(qū)域貿(mào)易對于地區(qū)水資源壓力的影響，并利用動態(tài)空間杜賓模型分別估計了“實體水-虛擬水”流動、跨流域調(diào)水以及糧食區(qū)域貿(mào)易對糧食整體種植結(jié)構(gòu)和水稻、小麥、玉米三大主糧作物種植結(jié)構(gòu)的影響及空間溢出效應(yīng)。主要結(jié)論如下。

（1）跨流域調(diào)水和糧食區(qū)域貿(mào)易對水資源壓力的影響研究表明，2004—2006年，跨流域調(diào)水和糧食區(qū)域貿(mào)易一定程度上減輕了全國層面的水資源壓力，2007—2020年跨流域調(diào)水和糧食區(qū)域貿(mào)易加劇了全國層面的水資源壓力。隨著糧食生產(chǎn)不斷向中西部和北方轉(zhuǎn)移和集中，糧食輸出區(qū)水資源壓力遠大于糧食輸入?yún)^(qū)，跨流域調(diào)水作用有限，全國水資源壓力總體呈上升趨勢。

（2）“實體水-虛擬水”流動短期和長期均提高了糧食輸出區(qū)糧食種植比例，降低了玉米種植比例，說明“實體水-虛擬水”流動對于糧食輸出區(qū)水資源壓力的增加，對糧食生產(chǎn)誘致性技術(shù)變遷的促進作用與區(qū)域生產(chǎn)分工的強化作用大于限制作用，且存在“灌溉效率悖論”現(xiàn)象，即節(jié)水技術(shù)進步導(dǎo)致高耗水作物種植增加?！皩嶓w水-虛擬水”流動短期和長期均降低了糧食輸入?yún)^(qū)糧食種植比例；短期水稻和玉米種植比例提高，且高耗水作物水稻種植比例提高幅度更大，長期水稻和玉米種植比例下降，且低耗水作物玉米種植比例下降程度更大，說明糧食輸入?yún)^(qū)糧食內(nèi)部結(jié)構(gòu)向耗水方向發(fā)展。

（3）“實體水-虛擬水”流動作用下，糧食輸入?yún)^(qū)對糧食輸出區(qū)糧食種植比例和玉米種植比例具有負向的空間溢出效應(yīng)，說明對糧食輸出區(qū)的糧食生產(chǎn)具有負向影響，且加劇了糧食內(nèi)部結(jié)構(gòu)的“灌溉效率悖論”現(xiàn)象，不利于長期的糧食安全與水資源安全。糧食輸出區(qū)對糧食輸入?yún)^(qū)的糧食種植比例不具有顯著影響，短期降低了水稻和玉米的種植比例，且高耗水作物水稻種植比例下降更多，長期提高了水稻和玉米種植比例，且低耗水作物玉米種植比例提高更多，說明糧食輸出區(qū)促進了糧食輸入?yún)^(qū)糧食內(nèi)部結(jié)構(gòu)向節(jié)水方向發(fā)展。

基于以上結(jié)論，提出如下政策啟示：①實施區(qū)域“實體水-虛擬水”統(tǒng)籌管理，綜合考慮“實體水-虛擬水”流動的地區(qū)效應(yīng)與空間溢出效應(yīng)，促進糧食種植結(jié)構(gòu)適水調(diào)整。②在當前“實體水-虛擬水”流動作用下，糧食輸出區(qū)糧食種植比例在提高的同時，糧食內(nèi)部結(jié)構(gòu)趨向耗水，不利于長期糧食安全和水資源安全，亟須大力促進糧食輸出區(qū)的水利工程建設(shè)和節(jié)水農(nóng)業(yè)發(fā)展。③當前“實體水-虛擬水”流動的空間溢出效應(yīng)，對糧食輸出區(qū)的糧食生產(chǎn)具有負向作用，且使得糧食輸出區(qū)糧食內(nèi)部結(jié)構(gòu)趨向耗水，糧食輸入?yún)^(qū)糧食內(nèi)部結(jié)構(gòu)趨向節(jié)水。應(yīng)完善糧食輸出區(qū)糧食生產(chǎn)利益補償機制，建立區(qū)域間糧食虛擬水補償制度，提高糧食輸出區(qū)水資源利用效率，促進糧食安全和水資源安全的可持續(xù)發(fā)展。