摘要：山東省是經(jīng)濟和人口大省，同時也是水資源緊缺省份。水足跡是認(rèn)識水資源的新視角，能夠揭示對水資源的真實占用?；谒阚E視角構(gòu)建生產(chǎn)水足跡測度模型和評價體系能對2010—2021年山東省水資源利用情況進行評價。結(jié)果表明：2010—2021年山東省生產(chǎn)水足跡波動上升，表現(xiàn)為西部和中南部高、東部和中北部低，藍(lán)、綠、灰水足跡均呈增長趨勢。在所有核算賬戶中，農(nóng)業(yè)占比最大。從綜合評價角度看，生產(chǎn)水足跡的經(jīng)濟效益逐漸增長，但單位生產(chǎn)水足跡的人口承載力較差，水資源壓力指數(shù)較大。山東省總體表現(xiàn)為水資源不可持續(xù)，水資源壓力較大，應(yīng)及時采取調(diào)整產(chǎn)業(yè)結(jié)構(gòu)、提高節(jié)水技術(shù)、增強污染防治等措施以緩解用水緊張的現(xiàn)狀。

關(guān)鍵詞：生產(chǎn)水足跡；時空格局；結(jié)構(gòu)差異；綜合評價；山東省

中圖分類號：TV213.4 文獻(xiàn)標(biāo)識碼：A 文章編號：1001-9235（2024）04-0107-09

Comprehensive Evaluation of Production Water Footprint in Shandong Province

ZHANG Min¹，ZHANG Yu^{1，2，3*}，TIAN Qing¹，WANG Mengfan¹

（1.School of Geography and Tourism，Qufu Normal University，Rizhao 276826，China;

2.Institute of Yellow River Ecology，Qufu Normal University，Qufu 273165，China;

3.Institute of Geographic Sciences and Natural Resources Research，Beijing 100101，China）

Abstract： As a significant economic and populous region，Shandong Province also grapples with water scarcity.The water footprint concept provides a fresh understanding of water resources，revealing the true utilization extent of the resources.A model and evaluation system for measuring the production water footprint is established to assess water resource utilization in the province from 2010 to 2021.The results indicate that the production water footprint in the studied region experiences fluctuations and presents an overall increasing trend during this period，with higher levels in the western and southern regions and lower ones in the eastern and northern areas.The blue，green，and grey water footprints all display a rising trend.Agriculture accounts for the largest proportion among all the sectors.From a comprehensive perspective，the economic benefits of the production water footprint steadily increase，but the per capita carrying capacity of the production water footprint is relatively poor，with a relatively large water resource pressure index.Generally，Shandong faces unsustainable water resources and considerable water pressure，which highlights the need for timely measures to adjust the industrial structure，enhance water-saving technologies，and improve pollution prevention and control to alleviate current water scarcity.

Keywords：production water footprint;spatio-temporal pattern;structural differences;comprehensive evaluation;Shandong Province

隨著近年來水資源供需矛盾日益加重，如何對水資源利用情況進行測度和評價成為了學(xué)者們關(guān)注的話題。水足跡包括實體水和虛擬水的特性，使其成為國家或區(qū)域水資源可持續(xù)利用評價、水資源壓力評價的新視角。山東省水資源嚴(yán)重短缺，嚴(yán)重制約經(jīng)濟社會可持續(xù)發(fā)展，故基于生產(chǎn)水足跡對山東省進行水資源利用狀況進行測度和評價，有利于為山東省水資源管理提供決策參考。

2002年荷蘭學(xué)者Hoekstra將實體水與虛擬水相聯(lián)系，提出“水足跡”的概念，即在一定物質(zhì)標(biāo)準(zhǔn)下，生產(chǎn)一定區(qū)域的人群所消耗的產(chǎn)品和服務(wù)需要的水資源數(shù)量^［1］。近幾年國內(nèi)外對于水足跡的研究較豐富，研究尺度涉及全球^［2-3］、國家^［4］、流域^［5］、省域或特定地區(qū)^［6-7］等。國內(nèi)主要集中在水足跡時空演變及影響因素分析^［8］、區(qū)域農(nóng)業(yè)用水評價^［9-10］、農(nóng)業(yè)生產(chǎn)結(jié)構(gòu)調(diào)整^［11］、糧食作物水足跡核算^［12］和區(qū)域間差異^［13］等方面。近幾年省域尺度的研究成果不斷涌現(xiàn)，如劉巍等^［14］、王月英等^［15］分別對黑龍江省、陜西省進行了研究。在基于水足跡理論的水資源評價方面，主要從水足跡的經(jīng)濟效益、生態(tài)指標(biāo)和水資源可持續(xù)性等角度構(gòu)建指標(biāo)體系評價水資源利用狀況，如商慶凱等^［16］衡量了青海省的進口、出口、農(nóng)業(yè)、工業(yè)、居民生活和生態(tài)水足跡，選取人均水足跡、水資源匱乏度、水資源壓力指數(shù)、水資源進口依賴度等指標(biāo)對水資源利用狀況進行評價。

本文從生產(chǎn)水足跡理論出發(fā)，在已有研究基礎(chǔ)上架構(gòu)生產(chǎn)水足跡測度模型和區(qū)域生產(chǎn)水足跡評價體系，以山東省為研究對象，從省域、市域、分賬戶水足跡、不同水資源類型等角度分析山東省2010—2021年生產(chǎn)水足跡的時空格局和結(jié)構(gòu)差異，并從效益、生態(tài)和可持續(xù)性3個指標(biāo)層面綜合評價山東省水資源利用情況，以期為山東省水資源科學(xué)管理和可持續(xù)發(fā)展提供有益參考。

1 研究區(qū)與數(shù)據(jù)

1.1 研究區(qū)域

山東省是工農(nóng)業(yè)和人口大省，經(jīng)濟發(fā)展居全國前列，其中農(nóng)業(yè)用水常年占總用水量的60%以上。但山東省水資源年際、區(qū)域分布不均，水資源緊缺是基本省情，2021年全省水資源總量為525.33億m³，僅占全國水資源總量的1.77%。同時，旱澇災(zāi)害頻發(fā)、嚴(yán)重依賴黃河水、地下水超采等問題^［17］也是制約工農(nóng)業(yè)發(fā)展和生活用水的重要瓶頸。所以本文從生產(chǎn)水足跡視角出發(fā)綜合評價山東省2010—2021年水資源利用情況，以期為緩解山東省水資源緊張及促進水資源合理分配提供理論和數(shù)據(jù)支持。

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

研究時段為2010—2021年，農(nóng)產(chǎn)品產(chǎn)量、工業(yè)用水、城鎮(zhèn)公共用水、居民生活用水、生態(tài)用水、污染物排放量、農(nóng)業(yè)總產(chǎn)值、人口等數(shù)據(jù)來自《山東省統(tǒng)計年鑒》《山東省水資源公報》、各市統(tǒng)計年鑒和水資源公報。農(nóng)業(yè)生產(chǎn)水足跡涉及種植業(yè)、畜牧業(yè)、林業(yè)和漁業(yè)四大農(nóng)業(yè)部門共34種產(chǎn)品，單位產(chǎn)品水足跡數(shù)據(jù)參考Mekonnen等^［2-3］有關(guān)中國的研究成果。部分畜牧業(yè)產(chǎn)品數(shù)據(jù)在此基礎(chǔ)上通過山東省飼料用谷物水足跡與全國平均比值以及全國畜牧產(chǎn)品平均水足跡得出，牛奶和淡水產(chǎn)品參考程雨菲等^［18］的研究。此外，因2019年萊蕪劃入濟南，所以本文統(tǒng)一將萊蕪各項數(shù)據(jù)納入濟南計算。

2 研究方法

2.1 生產(chǎn)水足跡測度

本文構(gòu)建了生產(chǎn)水足跡測度模型，包括農(nóng)業(yè)、工業(yè)、居民生活以及城鎮(zhèn)公共用水和生態(tài)環(huán)境用水，即式（1）：

WF_prod=WF_agriculture+WF_industry+WF_domestic+WF_public+WF_environment（1）

式中 WF_prod——總生產(chǎn)水足跡，m³；WF_agriculture——農(nóng)業(yè)生產(chǎn)水足跡，m³；WF_industry——工業(yè)生產(chǎn)水足跡，m³；WF_domestic——居民生活水足跡，m³；WF_public——城鎮(zhèn)公共用水量，m³；WF_environment——生態(tài)環(huán)境用水量，m³。

其中，農(nóng)業(yè)生產(chǎn)水足跡測度見式（2）：

式中 P_i——第i種農(nóng)產(chǎn)品產(chǎn)量；VW_i——單位產(chǎn)品水足跡。

工業(yè)生產(chǎn)水足跡包括工業(yè)用水和工業(yè)灰水足跡。工業(yè)灰水足跡指生產(chǎn)過程中排放的污染物達(dá)到自然濃度或水質(zhì)標(biāo)準(zhǔn)所需吸收淡水的理論體積，本文采用廢水中主要污染物（COD和氨氮排放量）為衡量指標(biāo)，取兩者最大值作為灰水足跡，見式（3）—（6）：

WF_industry=GWF_industry+Water_industry（3）

其中：

GWF_industry=max（GWF_{industry（COD）}，GWF_{industyr（NH}⁺_4－N））（4）

式中 WF_industry——工業(yè)生產(chǎn)水足跡，m³；Water_industry——工業(yè)用水量，m³；GWF_industry——工業(yè)灰水足跡，m³；GWF_{industry（COD）}——COD的灰水足跡，m³；GWF⁺_{industry（NH4-N）}——NH⁺₄-N的灰水足跡，m³；L_{industry（COD）}——COD的排放負(fù)荷，kg；L⁺_{industry（NH4-N）}——NH⁺₄-N的排放負(fù)荷，kg；C_max——環(huán)境水質(zhì)標(biāo)準(zhǔn)允許的污染物濃度，kg/m³，根據(jù)GB 8979—1996《污染物綜合排放標(biāo)準(zhǔn)》，C_max（COD）和C⁺_{max（NH4-N）}分別取0.060、0.015 kg/m³；根據(jù)《水足跡評價手冊》，C_nat取0 kg/m³。

居民生活水足跡包括居民生活用水和生活灰水足跡，計算方法同工業(yè)。

2.2 綜合評價體系

參考前人研究成果^［19］，從效益、生態(tài)和可持續(xù)性層面構(gòu)建了基于生產(chǎn)水足跡的水資源利用評價體系（表1）。水資源可持續(xù)指數(shù)等級見表2。

3 結(jié)果與討論

3.1 生產(chǎn)水足跡時空差異

3.1.1 年際變化

2010—2021年山東省生產(chǎn)水足跡為27 722.39億m³，整體呈波動上升趨勢（圖1）。農(nóng)業(yè)生產(chǎn)水足跡變化趨勢與總生產(chǎn)水足跡保持一致，反映了農(nóng)業(yè)作為用水大戶的主導(dǎo)作用。2011—2016年波動下降，2017年回升。2016年生產(chǎn)水足跡下降為2 219.76億m³，原因可能是2016年山東省的低溫受凍、風(fēng)雹、旱澇等自然災(zāi)害造成了農(nóng)產(chǎn)品大量減產(chǎn)^［23］，生產(chǎn)水足跡大幅減少。具體來看，煙臺、臨沂、濱州、菏澤上升幅度較大，而濟南、濟寧、泰安和濟寧顯著減少。

3.1.2 空間格局

本文利用ArcGIS繪制了山東省生產(chǎn)水足跡空間分布（圖2）。生產(chǎn)水足跡分布不均，整體呈西部、中南部高，東部、中北部低的特征。年均生產(chǎn)水足跡為2 310.20億m³，超200億m³的有菏澤、濰坊、臨沂、德州。這些地區(qū)均為農(nóng)業(yè)生產(chǎn)大市，種植業(yè)發(fā)達(dá)，農(nóng)業(yè)生產(chǎn)水足跡多。而日照、威海、東營和淄博較少，其農(nóng)業(yè)種植面積和產(chǎn)量較低，未消耗較多水資源。與之相對，其工業(yè)和居民生活生產(chǎn)水足跡占比較高。

山東省2010—2021年人均生產(chǎn)水足跡為2 335.49 m³，與前人研究差距較大^［24］，原因是研究年份不同、不同學(xué)者測算農(nóng)業(yè)生產(chǎn)水足跡方法和范圍不同，且不同學(xué)者也可能會采取不同污染物濃度標(biāo)準(zhǔn)測算工業(yè)灰水。各市人均生產(chǎn)水足跡由高到低為德州、濱州和聊城，均在3 000 m³以上，原因在于其地形平坦，耕地資源豐富，氣候及土地條件適宜農(nóng)業(yè)生產(chǎn)^［6］，第一產(chǎn)業(yè)比重大。較低的為淄博、青島和濟南。淄博和濟南地處魯中山區(qū)，年際、年內(nèi)降水分配不均，不利于作物生產(chǎn)，且相對耕地較少。青島多山地丘陵，耕地少，農(nóng)作物種植面積小。另外，青島和濟南是山東經(jīng)濟中心，經(jīng)濟重心為第二、三產(chǎn)業(yè)，且人口較多，所以人均生產(chǎn)水足跡較低。

3.2 生產(chǎn)水足跡結(jié)構(gòu)差異

3.2.1 水足跡結(jié)構(gòu)

山東省藍(lán)、綠、灰水足跡與生產(chǎn)水足跡變化趨勢相同（圖3），呈波動上升趨勢，占比分別為13.41%、61.91%和24.68%。各市生產(chǎn)水足跡中藍(lán)、綠、灰水足跡的構(gòu)成比例與山東省一致，均為綠水足跡＞藍(lán)水足跡＞灰水足跡（圖4）。綠水和灰水足跡較高的依次為濰坊、菏澤、臨沂，較低的依次為日照、東營和威海。藍(lán)水足跡較高的是菏澤、德州、濟寧，威海最低。

綠水足跡占主體且增速最快。綠水全部由農(nóng)業(yè)貢獻(xiàn)，主要來自種植業(yè)和林業(yè)消耗的有效降水及畜牧業(yè)和漁業(yè)飼料中的綠水^［25］。山東省綠水足跡波動增加，表明山東省對綠水的利用更加充分，有效緩解了藍(lán)水壓力。

藍(lán)水足跡雖增速較慢，但目前山東省地下水超采問題仍然存在，地下漏斗得不到有效恢復(fù)^［17］。水資源問題仍較嚴(yán)重。山東省應(yīng)大力發(fā)展多水源利用技術(shù)，提升節(jié)水潛力，做好“開源節(jié)流”^［26］，一方面要充分利用土壤水和雨洪水資源；另一方面要大范圍推廣節(jié)水灌溉和新型工農(nóng)業(yè)節(jié)水技術(shù)，建設(shè)高標(biāo)準(zhǔn)節(jié)水工程，保證藍(lán)水可持續(xù)利用。

灰水足跡波動增加但增速較緩。農(nóng)業(yè)和生活灰水波動增長，工業(yè)灰水逐年減少。農(nóng)業(yè)灰水增加表明農(nóng)藥化肥及其他廢棄物導(dǎo)致水污染加重。因此需加快建設(shè)污染處理系統(tǒng)，加大對農(nóng)業(yè)污染物處理的監(jiān)管力度。生活灰水主要來自廚房、洗滌等。面對經(jīng)濟發(fā)展和人口增加的用水需求和污水排放壓力，山東省需進一步通過水價改革、提高生活污水處理技術(shù)等方式提高生活用水使用率和回收率，降低生活灰水。工業(yè)灰水逐年減少，原因在于山東省增加了工業(yè)污染治理投資和治污設(shè)施運行費用投入^［27］，工業(yè)水污染逐漸減少。未來山東省應(yīng)繼續(xù)提高工業(yè)廢物處理技術(shù)，進一步減少工業(yè)污染。

3.2.2 核算賬戶結(jié)構(gòu)

山東省各核算賬戶生產(chǎn)水足跡占比由高到低為農(nóng)業(yè)、生活、工業(yè)、生態(tài)和城鎮(zhèn)公共用水（表3），占比分別為92.86%、4.22%、2.17%、0.43%和0.32%。

農(nóng)業(yè)各部門生產(chǎn)水足跡由高到低依次是種植業(yè)、畜牧業(yè)、漁業(yè)和林業(yè)。種植業(yè)、林業(yè)、畜牧業(yè)波動上升，漁業(yè)波動下降。種植業(yè)呈波動增長，但農(nóng)業(yè)用水量占比持續(xù)降低，原因在于近幾年山東省大力發(fā)展節(jié)水農(nóng)業(yè)，采用節(jié)水工藝和機械化作業(yè)，節(jié)水灌溉面積增加^［28］，用水效率提升。畜牧業(yè)年均占比為30.99%，畜產(chǎn)品單位產(chǎn)品水足跡較高使得畜牧業(yè)以約7.62%的產(chǎn)量貢獻(xiàn)了近1/3的水足跡。畜牧業(yè)生產(chǎn)水足跡增加是山東省近年來畜牧業(yè)機械化的成果，而且隨著經(jīng)濟發(fā)展水平的提高，人民生活水平提升，對生活品質(zhì)有了更高要求，增加了畜產(chǎn)品等高水足跡產(chǎn)品需求。漁業(yè)和林業(yè)水足跡年均占比分別為3.23%和0.57%，且年際變化較小。所以減少農(nóng)業(yè)水足跡應(yīng)主要從種植業(yè)和畜牧業(yè)入手。如優(yōu)化種植結(jié)構(gòu)，減少高耗水作物種植，種植蔬菜、水果等經(jīng)濟產(chǎn)值高且耗水少的產(chǎn)品；充分利用虛擬水貿(mào)易進口高耗水的畜禽肉、蛋奶產(chǎn)品來緩解水資源短缺壓力。

工業(yè)生產(chǎn)水足跡呈下降趨勢。工業(yè)藍(lán)水足跡波動上升，工業(yè)灰水足跡大幅下降，反映了近年來山東省工業(yè)轉(zhuǎn)型的巨大成效，工業(yè)灰水減少量大于工業(yè)藍(lán)水增加量，使得山東省工業(yè)生產(chǎn)水足跡緩慢降低。工業(yè)藍(lán)水足跡增加表明隨著工業(yè)發(fā)展消耗了較多水資源，所以山東省需要在大力發(fā)展工業(yè)生產(chǎn)的同時大力應(yīng)用工業(yè)節(jié)水技術(shù)。工業(yè)灰水足跡由2010年的49%下降至2021年的33%，表明隨著工業(yè)產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化，山東省加強了工業(yè)水污染管理，將灰水足跡控制在較低水平。

居民生活水足跡波動上升。藍(lán)水足跡平穩(wěn)增加，2010—2015年較穩(wěn)定，2016年后增長較快，主要原因在于經(jīng)濟發(fā)展和城鎮(zhèn)化進程推動了用水人口和需求增加，居民生活用水量增加。生活灰水足跡波動增加，灰水足跡約為藍(lán)水足跡的2.5倍，所以相對于節(jié)水，控制廢水中污染物的排放和治理是重中之重。

城鎮(zhèn)公共用水和生態(tài)用水均波動增加。城鎮(zhèn)公共用水年均增長率為1.04%，這是山東省城鎮(zhèn)化加快的結(jié)果^［29］，城鎮(zhèn)化帶來了城市人口膨脹和城鎮(zhèn)用地規(guī)模擴張，對城鎮(zhèn)公共建設(shè)用水需求增大。雖然生態(tài)用量總量較少，但增長最快，年均增長13.81%，因生態(tài)文明建設(shè)提出之后，山東省提高了生態(tài)意識，保證了生態(tài)環(huán)境蓄水量，注重生態(tài)建設(shè)與經(jīng)濟發(fā)展的協(xié)調(diào)。在城鎮(zhèn)化進程中，山東省應(yīng)重視發(fā)展節(jié)水技術(shù)，減少空間延展導(dǎo)致的水量損失，提高水資源利用率，同時注意避免城鎮(zhèn)化建設(shè)中生態(tài)環(huán)境用水被擠占，保障生態(tài)用水供應(yīng)量。

3.3 綜合評價

3.3.1 效益指標(biāo)

首先，2010—2021年山東省生產(chǎn)水足跡經(jīng)濟效益整體上升（表4），說明生產(chǎn)水足跡在經(jīng)濟發(fā)展中的作用不斷提高，即生產(chǎn)水足跡經(jīng)濟產(chǎn)出逐年增加。青島和濟南經(jīng)濟效益較高，菏澤、德州和聊城較低，這與各市產(chǎn)業(yè)結(jié)構(gòu)有關(guān)。但山東省年均僅為27.4元/m³，明顯低于其他地區(qū)水足跡的經(jīng)濟效益^［30］，所以不斷提高水資源利用率和單位經(jīng)濟價值是未來努力的方向。其次，山東省生產(chǎn)水足跡人口效益略有波動，保持在4.11～4.48人/萬m³，年均4.28人/萬m³，說明山東省生產(chǎn)水足跡人口承載力較穩(wěn)定，2010—2016年波動上升，2017—2021年波動下降。但劉燦輝等^［5］計算的太湖流域生產(chǎn)水足跡人口密度值為16.46～19.23人/萬m³，約為山東省的4倍，即山東省生產(chǎn)水足跡的人口效益較低，單位水足跡人口承載力較差。山東省的低人口承載力表明生產(chǎn)消耗的水資源在區(qū)域內(nèi)起到的有效作用力較差^［19］，有關(guān)部門需引起重視并采取有效措施提高用水效率，從而提升人口效益。最后，水資源負(fù)載指數(shù)均大于10，呈先波動上升后波動下降的趨勢，即山東省水資源開發(fā)利用水平較高，日后開發(fā)困難，需要依靠跨流域調(diào)水補充水資源。

3.3.2 生態(tài)指標(biāo)

山東省生產(chǎn)水足跡壓力指數(shù)年均值為9.11（表4），表明山東省生產(chǎn)耗水量遠(yuǎn)超地區(qū)可利用水資源量，即水資源長期處于負(fù)均衡的超載狀態(tài)，水資源消耗量遠(yuǎn)超可利用水資源量，在一定程度上反映了山東省地下水超采問題依然存在。從變化趨勢角度看，山東省總體呈現(xiàn)先增加后波動降低的趨勢，隨著農(nóng)業(yè)節(jié)水灌溉面積不斷擴大、工業(yè)節(jié)水技術(shù)進步和南水北調(diào)工程的推動^［31-32］，地下水超采問題得以緩解，水資源緊張狀態(tài)有所緩和，但仍面臨巨大壓力。2021年的水資源壓力指數(shù)最低，2021年受大氣環(huán)流異常等的影響，山東省降水量大幅增加，比多年平均偏多24.5%，屬于豐水年份^［33］。2014年水資源壓力最高，原因在于2014年屬于大旱之年^［34］。即山東省水資源壓力不僅與本年度生產(chǎn)水足跡有關(guān)，也與降水補給程度有極大的關(guān)系，所以山東省在采取措施減少生產(chǎn)水足跡的同時，還應(yīng)通過植樹造林、恢復(fù)植被、修建水庫和水利工程、退田還湖等人類活動增加空氣濕度，進而增加降水量。

3.3.3 可持續(xù)性指標(biāo)

山東省生產(chǎn)水足跡增長指數(shù)年均為0.54%，其中2012—2013、2015—2016、2018年為負(fù)值?？衫盟Y源增長指數(shù)年均為12.41%，年際波動較大，原因在于山東省降水量年際變化較大，導(dǎo)致地表水和地下水資源量變動幅度較大。山東省水資源可持續(xù)指數(shù)年均為12.31%，年際變化為先下降后上升，除2021年為中不可持續(xù)外，其余年份均為強不可持續(xù)。這表明山東省水資源可持續(xù)性處于危機狀態(tài)，水資源壓力和水資源緊張狀況必須引起重視。所以山東省必須因地制宜，合理調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，降低高耗水工業(yè)和農(nóng)業(yè)的生產(chǎn)水足跡比重，提高節(jié)水能力和污水處理技術(shù)，同時完善水利設(shè)施，保障調(diào)水和蓄水，緩解地下水開發(fā)壓力，發(fā)展節(jié)水技術(shù)。

4 結(jié)論

本文基于水足跡理論，對山東省及其16地市生產(chǎn)水足跡年際變化、空間格局和結(jié)構(gòu)差異進行研究，從效益、生態(tài)和可持續(xù)性3個層面構(gòu)建了區(qū)域生產(chǎn)水足跡評價體系并對山東省水資源利用狀況進行實證研究。

a）時空格局。年際變化角度，2010—2021年山東省生產(chǎn)水足跡整體呈波動上升趨勢，各市變化特征與山東省保持一致?？臻g格局角度，生產(chǎn)水足跡整體呈現(xiàn)西部、中南部高，東部、中北部低的分布特征；西北部和西南部人均生產(chǎn)水足跡較高；山東省年均生產(chǎn)水足跡為2 310.20億m³。

b）結(jié)構(gòu)差異。水資源結(jié)構(gòu)角度，山東省藍(lán)、綠、灰水足跡均呈增長趨勢，占比分別為13.41%、61.91%和24.68%，綠水足跡增長速率最快；各市藍(lán)綠灰水足跡比例與全省基本一致。核算賬戶結(jié)構(gòu)上，由高到低為農(nóng)業(yè)、生活、工業(yè)、生態(tài)和城鎮(zhèn)公共用水，農(nóng)業(yè)占比高達(dá)92.86%。

c）綜合評價。2010—2021年山東省生產(chǎn)水足跡經(jīng)濟效益逐漸增長，單位水足跡人口承載力較差，水資源壓力指數(shù)均較大，山東省水資源利用總體表現(xiàn)為強不可持續(xù)，水資源壓力較大，應(yīng)及時采取調(diào)整產(chǎn)業(yè)結(jié)構(gòu)、提高節(jié)水技術(shù)、增強污染防治等措施以緩解農(nóng)業(yè)用水緊張的現(xiàn)狀。

參考文獻(xiàn)：

［1］HOEKSTRA A，HUNG P Q.Virtual water trade：A quantification of virtual water flows between nations in relation to international crop trade［J］.Water Science and Technology，2002，49：203-209.

［2］MEKONNEN M M，HOEKSTRA A Y.The green，blue and grey water footprint of crops and derived crop products，Value of Water Research Report Series No.47［R］.Netherlands：Delft，2010.

［3］MEKONNEN M M，HOEKSTRA A Y.The green，blue and grey water footprint of farm animals and animal products，Value of Water Research Report Series No.48［R］.Netherlands：Delft，2010.

［4］孫世坤，王玉寶，劉靜，等.中國主要糧食作物的生產(chǎn)水足跡量化及評價［J］.水利學(xué)報，2016，47（9）：1115-1124.

［5］劉燦輝，金文龍，孫英杰，等.水足跡視角下的太湖流域水資源利用評價［J］.中國農(nóng)村水利水電，2022（8）：70-77.

［6］余灝哲，韓美.基于水足跡的山東省水資源可持續(xù)利用時空分析［J］.自然資源學(xué)報，2017，32（3）：474-483.

［7］鄭亨，杜海龍，黃靖.杭州市2010—2016年水足跡分析研究［J］.人民珠江，2019，40（5）：69-75.

［8］郭相平，高爽，吳夢洋，等.中國農(nóng)作物水足跡時空分布與影響因素分析［J］.農(nóng)業(yè)機械學(xué)報，2018，49（5）：295-302.

［9］王旭，高翔，曹蕾，等.水足跡視角下的寧夏中衛(wèi)市農(nóng)業(yè)水資源利用評價［J］.蘭州大學(xué)學(xué)報（自然科學(xué)版），2015，51（5）：619-624.

［10］操信春，任杰，吳夢洋，等.基于水足跡的中國農(nóng)業(yè)用水效果評價［J］.農(nóng)業(yè)工程學(xué)報，2018，34（5）：1-8.

［11］侯慶豐.基于水足跡的甘肅省農(nóng)作物種植結(jié)構(gòu)優(yōu)化分析［J］.中國沙漠，2013，33（6）：1921-1927.

［12］田園宏，諸大建，王歡明，等.中國主要糧食作物的水足跡值：1978—2010［J］.中國人口·資源與環(huán)境，2013，23（6）：122-128.

［13］孫世坤，王玉寶，吳普特，等.小麥生產(chǎn)水足跡區(qū)域差異及歸因分析［J］.農(nóng)業(yè)工程學(xué)報，2015，31（13）：142-148.

［14］劉巍，李雨晴，付強，等.黑龍江省糧食水足跡時空分布規(guī)律與影響因素分析［J］.農(nóng)業(yè)機械學(xué)報，2020，51（8）：214-222，335.

［15］王月英，劉科偉.陜西省農(nóng)業(yè)水足跡量化及對水資源的影響分析［J］.河北農(nóng)業(yè)大學(xué)學(xué)報，2019，42（4）：109-116.

［16］商慶凱，陰柯欣，米文寶.基于水足跡理論的青海省水資源利用評價［J］.干旱區(qū)資源與環(huán)境，2020，34（5）：70-77.

［17］董亞楠，邢立亭，張欣慧，等.地下水超采區(qū)評價法及其應(yīng)用［J］.地質(zhì)與勘探，2019，55（2）：641-648.

［18］程雨菲，曹升樂，楊裕恒，等.山東省農(nóng)業(yè)水足跡周期及趨勢變化分析［J］.南水北調(diào)與水利科技，2018，16（4）：169-175.

［19］戚瑞，耿涌，朱慶華.基于水足跡理論的區(qū)域水資源利用評價［J］.自然資源學(xué)報，2011，26（3）：486-495.

［20］史利潔，吳普特，王玉寶，等.基于作物生產(chǎn)水足跡的陜西省水資源壓力評價［J］.中國生態(tài)農(nóng)業(yè)學(xué)報，2015，23（5）：650-658.

［21］劉獻(xiàn)，袁旦，韓軼菲，等.基于水足跡理論的甘肅省水資源利用與經(jīng)濟增長脫鉤分析研究［J］.人民珠江，2022，43（7）：52-60.

［22］甘容，李旖旎，楊峰，等.基于三維水足跡模型的河南省水資源可持續(xù)利用水平評價［J］.水資源與水工程學(xué)報，2023，34（1）：30-39.

［23］于小兵，陳虹，吉中會，等.山東省農(nóng)業(yè)氣象災(zāi)害與主要農(nóng)作物產(chǎn)量的灰色關(guān)聯(lián)分析［J］.災(zāi)害學(xué)，2019，34（3）：103-108.

［24］何治韋.基于水足跡的生態(tài)補償機制研究［J］.水利科技與經(jīng)濟，2019，25（11）：26-32.

［25］賈浩，李碩，宋英杰，等.滄州市農(nóng)業(yè)水足跡變化特征與驅(qū)動因素分析［J］.中國農(nóng)業(yè)大學(xué)學(xué)報，2018，23（6）：1-14.

［26］杜貞棟，王昕，周榮星.山東省農(nóng)業(yè)適水發(fā)展對策研究［J］.山東水利，2021（1）：1-4，11.

［27］宋雪青.高技術(shù)產(chǎn)業(yè)技術(shù)創(chuàng)新能力及影響因素研究［D］.廣州：廣東財經(jīng)大學(xué)，2022.

［28］李明輝.山東糧食生產(chǎn)水資源配置及優(yōu)化策略研究［D］.濟南：山東農(nóng)業(yè)大學(xué)，2019.

［29］王迎婧.城鎮(zhèn)居民采用節(jié)水器具的影響因素研究［D］.西安：西安理工大學(xué)，2022.

［30］焦士興，郭力菲，王安周，等.基于水足跡理論的河南省水資源利用評價［J］.人民黃河，2021，43（11）：87-91，96.

［31］馬承新，鄧海燕，王維平.山東省與其他區(qū)域間虛擬水貿(mào)易定量研究［J］.灌溉排水學(xué)報，2015，34（7）：54-56.

［32］楊亞鋒，高媛媛，許新宜，等.水資源生態(tài)安全的時空分異特征：以南水北調(diào)東線山東段為例［J］.南水北調(diào)與水利科技（中英文），2022，20（5）：833-841.

［33］山東省水利廳.山東省水資源公報［M］.濟南：山東出版社，2010—2021.

［34］盧書梅，楊傳國，谷黃河.1959—2018年山東省降水時空分布特征與洪旱演變規(guī)律［J］.水資源與水工程學(xué)報，2020，31（6）：60-66.

（責(zé)任編輯：程 茜）