虛擬水貿(mào)易的可計算非線性動態(tài)投入產(chǎn)出分析模型

馬超+許長新+田貴良+吳丹+石常峰+吳兆丹

摘要 虛擬水是水資源需求管理的創(chuàng)新領(lǐng)域。投入產(chǎn)出法是研究虛擬水的重要方法，但既有研究多以線性靜態(tài)投入產(chǎn)出模型為主，存在較強(qiáng)的比例性假設(shè)，在解釋現(xiàn)實(shí)經(jīng)濟(jì)活動中的虛擬水貿(mào)易時存在一定缺陷和不足。為更加貼近現(xiàn)實(shí)狀況，本文利用可計算一般均衡（CGE）思想，設(shè)計了一種新的虛擬水測算思路，對傳統(tǒng)的線性靜態(tài)投入產(chǎn)出模型進(jìn)行了非線性和動態(tài)化的拓展，旨在對區(qū)域經(jīng)濟(jì)系統(tǒng)中的虛擬水貿(mào)易進(jìn)行更加科學(xué)合理的計算和考察。首先，參照國家統(tǒng)計機(jī)構(gòu)常用的42產(chǎn)業(yè)部門劃分方式，根據(jù)一般均衡理論，圍繞區(qū)域經(jīng)濟(jì)系統(tǒng)的生產(chǎn)模塊、價格模塊和供需平衡模塊三個部分，定義有關(guān)變量和參數(shù)，對區(qū)域經(jīng)濟(jì)系統(tǒng)非線性動態(tài)投入產(chǎn)出模型進(jìn)行了詳細(xì)的方程列寫。在此基礎(chǔ)上，將虛擬水流動的因素與一般形式的非線性動態(tài)投入產(chǎn)出模型進(jìn)行嵌套，采用“母表”（價值型流量表）和“子表”（水資源流量表）相結(jié)合的形式，給出了一種全新的區(qū)域水資源投入產(chǎn)出表的設(shè)計思路與編制方法，將可計算非線性動態(tài)產(chǎn)出模型從一般形式擴(kuò)展至水資源領(lǐng)域，構(gòu)建了區(qū)域經(jīng)濟(jì)系統(tǒng)中虛擬水貿(mào)易的可計算非線性動態(tài)投入產(chǎn)出分析框架。分析同時指出，通過對該模型進(jìn)行求解，可以推導(dǎo)和計算平衡增長路徑和最優(yōu)增長路徑下的區(qū)域經(jīng)濟(jì)系統(tǒng)產(chǎn)出結(jié)構(gòu)和用水結(jié)構(gòu)，并可結(jié)合虛擬水貿(mào)易對區(qū)域經(jīng)濟(jì)的作用機(jī)理，分析最優(yōu)增長路徑下實(shí)施虛擬水貿(mào)易對區(qū)域經(jīng)濟(jì)增長的貢獻(xiàn)。

關(guān)鍵詞 虛擬水貿(mào)易；可計算；非線性；動態(tài)；投入產(chǎn)出

中圖分類號 TV213.4 文獻(xiàn)標(biāo)識碼 A 文章編號 1002-2104（2016）11-0160-10

doi：10.3969/j.issn.1002-2104.2016.11.020

虛擬水貿(mào)易是無形的“引調(diào)水工程”[1-2]?！疤摂M水”的概念由Tony Allan教授于1993年首先提出[3-5]，與實(shí)體水相對應(yīng)，是指生產(chǎn)商品或服務(wù)中所需要的水資源量[6-7]，以無形的方式蘊(yùn)含在產(chǎn)品和服務(wù)中。虛擬水貿(mào)易，則是指缺水國家或地區(qū)通過貿(mào)易的方式從富水國家和地區(qū)購買水密集產(chǎn)品來緩解本國或本地區(qū)的水資源壓力，實(shí)現(xiàn)當(dāng)?shù)氐乃Y源安全及食品安全。其理論的提出與發(fā)展拓展了過去人們對于水與糧食問題上從問題發(fā)生范圍內(nèi)尋求解決方案的傳統(tǒng)思路，而是強(qiáng)調(diào)從系統(tǒng)的角度出發(fā)，基于系統(tǒng)思考的方法從問題發(fā)生的范圍之外尋找解決問題的對策，為水資源管理特別是需求側(cè)管理開辟了新的視角。虛擬水的量化分析和政策研究，也由此成為水問題政策專家與學(xué)者的關(guān)注焦點(diǎn)與創(chuàng)新領(lǐng)域。我們試圖對傳統(tǒng)的線性靜態(tài)水資源投入產(chǎn)出模型加以擴(kuò)展和改進(jìn)，引用當(dāng)前較為前沿的可計算非線性動態(tài)投入產(chǎn)出模型，基于可計算非線性動態(tài)投入產(chǎn)出的思想，建立區(qū)域經(jīng)濟(jì)系統(tǒng)非線性動態(tài)水資源投入產(chǎn)出模型，對區(qū)域經(jīng)濟(jì)系統(tǒng)中的虛擬水貿(mào)易進(jìn)行更加科學(xué)合理的計算和考察。本文將重點(diǎn)給出可計算非線性動態(tài)水資源I-O模型的構(gòu)建思路。

1 文獻(xiàn)綜述

在虛擬水的量化計算中，國內(nèi)外學(xué)者常用的方法有兩種：一種是單位水量法，主要來自Chapagain和Hoekstra提出的研究不同產(chǎn)品生產(chǎn)樹的方法[7]，以及Zimmer和Renault基于對不同產(chǎn)品類型的區(qū)分的計算方法[8]。比如，在計算農(nóng)產(chǎn)品的虛擬水含量時，普遍采用聯(lián)合國糧農(nóng)組織（FAO）推薦并修正的標(biāo)準(zhǔn)彭曼-孟蒂斯公式，計算作物在生長發(fā)育期間蒸發(fā)蒸騰所消耗的全部水資源量ETc（Plant Evapotranspiration），并對農(nóng)作物產(chǎn)量YC進(jìn)行統(tǒng)計，即可得到單位質(zhì)量農(nóng)產(chǎn)品的虛擬水含量[9-11]。另一種是投入產(chǎn)出法。與生產(chǎn)樹法相比，后者因為其所依據(jù)的投入產(chǎn)出表覆蓋了國民經(jīng)濟(jì)系統(tǒng)中各種產(chǎn)品的生產(chǎn)過程和使用去向，可以在產(chǎn)業(yè)普遍聯(lián)系之中，更為直觀、準(zhǔn)確的追蹤產(chǎn)品貿(mào)易背后的水資源流動效應(yīng)，因而成為研究社會生產(chǎn)系統(tǒng)中水資源問題最重要的分析工具之一[12-14]。近年來，大量學(xué)者使用投入產(chǎn)出模型對虛擬水貿(mào)易進(jìn)行計算分析。Yu等用投入產(chǎn)出模型計算了英國虛擬水貿(mào)易格局[15]。Antonelli等利用投入產(chǎn)出法計算了地中海地區(qū)“藍(lán)水”和“綠水”的流動情況[16]。田貴良和許長新基于投入產(chǎn)出分析方法，構(gòu)建起虛擬水貿(mào)易理論框架下的水價敏感性模型，分析了水價上漲對各產(chǎn)業(yè)產(chǎn)品價格帶來的影響[17]。Guan等用投入產(chǎn)出法，計算了我國北方6省和廣東省的虛擬水流動情況[18]。譚圣林提出運(yùn)用單區(qū)域投入產(chǎn)出表和多區(qū)域投入產(chǎn)出表研究虛擬水，并以廣東省為例進(jìn)行了量化分析[19]。李方一對山西省與國內(nèi)其他地區(qū)的虛擬水流動關(guān)系進(jìn)行了實(shí)證研究[20]。雷玉桃依據(jù)投入產(chǎn)出分析原理，對中國17個行業(yè)的各種用水系數(shù)和虛擬水進(jìn)出口量進(jìn)行了計算分析[21]。周姣等建立起價值-實(shí)物混合型區(qū)域間投入產(chǎn)出表對華北地區(qū)的虛擬水貿(mào)易進(jìn)行了量化分析[22]。黃曉榮等運(yùn)用投入產(chǎn)出法計算了1998-2002年寧夏地區(qū)的虛擬水狀況[23]。馬忠等利用投入產(chǎn)出法對甘肅省張掖市各產(chǎn)業(yè)間的虛擬水流動進(jìn)行了分析[24]?？偟目磥恚S著虛擬水概念的逐漸普及和研究領(lǐng)域的不斷擴(kuò)大，投入產(chǎn)出法以其行業(yè)覆蓋面廣、計算精確度高的特點(diǎn)，成為當(dāng)前虛擬水貿(mào)易研究中的主流方法。

但需要看到，既有研究中的投入產(chǎn)出模型多以線性靜態(tài)模型為主，存在著很強(qiáng)的比例性假設(shè)，即假設(shè)在一定的值域內(nèi)，某種產(chǎn)品的產(chǎn)出量與包括水資源在內(nèi)的投入要素之間是成線性比例的，各投入要素間必須是固定比例而不能變化的[25]。而事實(shí)上，在實(shí)際生產(chǎn)過程中，各產(chǎn)業(yè)內(nèi)部投入要素的比例是在不斷變化的。比如，發(fā)展中國家處于經(jīng)濟(jì)起步階段時，往往投入較多的勞動力資源而使用較少的固定資本，但隨著經(jīng)濟(jì)發(fā)展，固定資本占投入要素的比例不斷增加，而投入的勞動力總量基本變化不大，導(dǎo)致其投入比例的減小，由此引起了兩種投入要素其比例隨著生產(chǎn)的進(jìn)行在不斷變化。由此看來，基于比例性假設(shè)的線性靜態(tài)投入產(chǎn)出模型在描述現(xiàn)實(shí)經(jīng)濟(jì)活動中的投入產(chǎn)出關(guān)系時不夠準(zhǔn)確。為了使投入產(chǎn)出模型更加切近現(xiàn)實(shí)中的區(qū)域經(jīng)濟(jì)系統(tǒng)，近年來，許多經(jīng)濟(jì)學(xué)家都致力于將線性投入產(chǎn)出模型推廣到非線性的情況，并取得一定進(jìn)展，但尚未有研究將這一思路應(yīng)用到水資源及虛擬水貿(mào)易領(lǐng)域。

2 模型構(gòu)建的思路

本模型的構(gòu)建思路源于可計算一般均衡（CGE）思想。這一思想源于經(jīng)濟(jì)學(xué)家瓦爾拉斯的一般均衡理論，該理論有別于數(shù)據(jù)導(dǎo)向的計量經(jīng)濟(jì)學(xué)模型，而是依賴堅實(shí)的經(jīng)濟(jì)理論基礎(chǔ)，將整個經(jīng)濟(jì)系統(tǒng)作為研究對象，全面考察系統(tǒng)中各種商品和要素供給、需求和供求的變化關(guān)系，成為研究經(jīng)濟(jì)復(fù)雜系統(tǒng)和政策分析的有效工具[26-28]。在CGE框架中，各經(jīng)濟(jì)主體和市場都將價格視為參數(shù)，并通過價格相互作用，既體現(xiàn)經(jīng)濟(jì)主體彼此之間的聯(lián)系，又包含對市場機(jī)制的描述[29]，一般框架如圖1所示。而水資源作為資源環(huán)境可持續(xù)發(fā)展的基礎(chǔ)性資源，與經(jīng)濟(jì)社會關(guān)系密切、相互依存，幾乎所有行業(yè)的生產(chǎn)和消費(fèi)都直接或間接與水有關(guān)。在接下來可計算非線性動態(tài)水資源投入產(chǎn)出模型的構(gòu)建中，本文將結(jié)合CGE的思想進(jìn)行闡釋。

國內(nèi)學(xué)者曾給出6部門區(qū)域經(jīng)濟(jì)系統(tǒng)非線性動態(tài)投入產(chǎn)出模型的計算方法[25]。從本質(zhì)上講，如果采用上述原理與方法進(jìn)一步構(gòu)造更多部門的可計算非線性動態(tài)投入產(chǎn)出模型，除了在數(shù)據(jù)采集及參數(shù)估計方面要求較大的工作量，在技術(shù)上是完全可行的?？紤]建模、求解的復(fù)雜程度以及統(tǒng)計數(shù)據(jù)的可得性，本研究擬擴(kuò)展到投入產(chǎn)出分析中更為具體的42產(chǎn)業(yè)部門，建立42部門可計算非線性動態(tài)水資源I-O模型。需要說明的是，本文對非線性投入產(chǎn)出模型中相關(guān)參數(shù)的設(shè)置，與投入產(chǎn)出表的元素結(jié)構(gòu)基本一致。

3 區(qū)域經(jīng)濟(jì)系統(tǒng)非線性動態(tài)投入產(chǎn)出模型方程列寫

根據(jù)一般均衡理論，區(qū)域經(jīng)濟(jì)系統(tǒng)由生產(chǎn)模塊、價格模塊和供需平衡模塊三部分組成。本小節(jié)分別構(gòu)建前兩模塊的非線性動態(tài)方程組，其中生產(chǎn)模塊包括總產(chǎn)出方程、增加值生產(chǎn)函數(shù)、固定資產(chǎn)生產(chǎn)函數(shù)、存貨生產(chǎn)函數(shù)和資本存量平衡方程；價格模塊包括產(chǎn)品價格方程、要素價格方程、增加值價格模型。

參照國家統(tǒng)計局慣用的部門劃分方式，將國民經(jīng)濟(jì)劃分為42個部門，其中第一產(chǎn)業(yè)包括1個部門，即農(nóng)林牧漁業(yè)，簡稱農(nóng)業(yè)，代碼為01；第二產(chǎn)業(yè)包括25個部門，代碼為02到26；第三產(chǎn)業(yè)包括16個部門，代碼為27到42。為體現(xiàn)區(qū)域經(jīng)濟(jì)系統(tǒng)固定資本的實(shí)際變化，我們在建立非線性動態(tài)投入產(chǎn)出模型時，將建筑業(yè)部門（代碼26）從工業(yè)部門中單獨(dú)分離進(jìn)行分析。

3.1 生產(chǎn)模塊方程

關(guān)于區(qū)域經(jīng)濟(jì)系統(tǒng)生產(chǎn)模塊的數(shù)學(xué)描述，主要包括各產(chǎn)業(yè)部門的生產(chǎn)動態(tài)方程、固定資本累計與投資方程兩部分。下面，首先給出區(qū)域經(jīng)濟(jì)系統(tǒng)生產(chǎn)動態(tài)方程，定義變量見表1：

需要強(qiáng)調(diào)的是，表中所定義的變量與投入產(chǎn)出表中的定義完全一致。另外，在以上42個部門或者42種產(chǎn)品中，只有工業(yè)產(chǎn)品及建筑業(yè)產(chǎn)品可以作為固定資本，由于每個部門其使用的固定資本總額Kj中包括工業(yè)品固定資本Kj′和建筑業(yè)產(chǎn)品固定資本Kj″，一般的，有如下平衡關(guān)系：

對于第j個部門，其投入產(chǎn)出關(guān)系可用圖2來表示：

對于圖2所示的第一層關(guān)系，根據(jù)宏觀經(jīng)濟(jì)學(xué)的基本知識，以及區(qū)域經(jīng)濟(jì)系統(tǒng)投入產(chǎn)出經(jīng)濟(jì)流量表的平衡關(guān)系，有：

產(chǎn)業(yè)部門j中間投入的i產(chǎn)品數(shù)量xij由j部門的經(jīng)濟(jì)產(chǎn)出Yj和中間投入系數(shù)aij決定，后者可以直接從投入產(chǎn)出表直接獲取。對于j部門的增加值Vj，根據(jù)投入產(chǎn)出表的平衡方程，則有：

由式2和式3可得到區(qū)域經(jīng)濟(jì)系統(tǒng)投入產(chǎn)出關(guān)系第一層關(guān)系的數(shù)學(xué)方程，列寫如下：

第二層關(guān)系表示了生產(chǎn)活動中部門增加值與固定資本投資投入產(chǎn)出的數(shù)量關(guān)系，其數(shù)學(xué)表達(dá)式可用學(xué)者慣用的CES類型的生產(chǎn)函數(shù)來描述。方程如下：

其中，αj、dj、σj是j產(chǎn)業(yè)部門CES函數(shù)的參數(shù)。

下面，考慮圖2的第三層投入產(chǎn)出關(guān)系，該部分用以反映固定資本投資的構(gòu)成比例?？疾焱度氘a(chǎn)出表的結(jié)構(gòu)不難發(fā)現(xiàn)，在實(shí)際的固定資本積累中只有工業(yè)品與建筑業(yè)品可作為固定資本投資，對于固定投資中的其他部分由于無法在投入產(chǎn)出結(jié)構(gòu)中明顯體現(xiàn)，本文不予考慮。不妨假設(shè)部門j工業(yè)產(chǎn)品投資占固定資本的比例為βj，則建筑業(yè)品所占比例為1-βj。其關(guān)系方程很容易由下式表達(dá)：

用列昂惕夫型的函數(shù)表示，上式轉(zhuǎn)化為：

同樣的，βj值也可以依據(jù)區(qū)域經(jīng)濟(jì)系統(tǒng)歷年的統(tǒng)計數(shù)據(jù)進(jìn)行參數(shù)估算。以上方程構(gòu)成了區(qū)域經(jīng)濟(jì)系統(tǒng)各部門的生產(chǎn)方程。

下面，仍以產(chǎn)業(yè)部門j為例，對固定資本累計與投資方程進(jìn)行描述。為體現(xiàn)固定資本的累計效應(yīng)，需設(shè)定基準(zhǔn)年t及基準(zhǔn)年的后續(xù)年t+1，此外還需設(shè)定固定資產(chǎn)折舊率δ′和δ″，分別代表工業(yè)品固定資本及建筑業(yè)固定資本的折舊率，部門j在第t年使用的工業(yè)品及建筑業(yè)品固定資本數(shù)量分別為Ij′（t）和Ij″（t）。由此，我們對部門j的固定資本積累與投資方程表述如下：

至此，區(qū)域經(jīng)濟(jì)系統(tǒng)第j個產(chǎn)業(yè)部門的生產(chǎn)函數(shù)、固定資本投資及資本積累的非線性動態(tài)方程推導(dǎo)完畢，聯(lián)立方程組如下：

3.2 價格模塊方程

區(qū)域經(jīng)濟(jì)系統(tǒng)價格方程可由生產(chǎn)函數(shù)方程進(jìn)行演推，根據(jù)公式（9）的第一式，即列昂惕夫型的生產(chǎn)函數(shù)，不妨假設(shè)j部門產(chǎn)品的產(chǎn)出Yj量為1，反映在方程上，此時中間投入xij等于中間投入系數(shù)aij，即：

對于增加值Vj，有：

類似的，我們首先對價格方程推導(dǎo)中涉及的諸多變量進(jìn)行定義，詳見表2。

其中，pj和qj分別為j部門的產(chǎn)品價格及增加值價格。應(yīng)當(dāng)注意，該價格并非現(xiàn)實(shí)中產(chǎn)品價格的絕對值，而是價格指數(shù)的概念。這是因為，對于價值型投入產(chǎn)出表，每單位的經(jīng)濟(jì)產(chǎn)出并非數(shù)量概念而是代表價值概念，所以這里的pj和qj均指調(diào)價幅度。如pj=1.2，則對于j部門1億元的經(jīng)濟(jì)產(chǎn)出而言，是指當(dāng)前的產(chǎn)品價格被低估，應(yīng)調(diào)整為1.2億元。在市場平衡時，一種產(chǎn)品的市場價格應(yīng)和實(shí)際投入的成本相等。

接下來，考慮生產(chǎn)函數(shù)非線性方程的第三式，即：

用價格方程來表示生產(chǎn)函數(shù)中固定資本的構(gòu)成比例，首先應(yīng)明確工業(yè)固定資本與建筑業(yè)品固定資本的價格指數(shù)。對于后者，由于在投入產(chǎn)出表中由唯一的部門與之對應(yīng)，可直接由該部門的價格指數(shù)予以表示，而工業(yè)固定資本則包含了投入產(chǎn)出表中的多個部門，由于模型中的參數(shù)設(shè)置與投入產(chǎn)出表結(jié)構(gòu)完全一致，因此不能用某一參數(shù)直接給出。這里，我們考慮用區(qū)域經(jīng)濟(jì)系統(tǒng)各工業(yè)部門的平均價格指數(shù)來表示工業(yè)固定資本的價格指數(shù)。如，對于42部門投入產(chǎn)出表，其工業(yè)部門包含了代碼由02至25總共24個經(jīng)濟(jì)部門，因此其固定資本價格指數(shù)可表示為：p′=p2+p3+…+p2524，而建筑業(yè)固定資本的價格指數(shù)為p″=p26。簡明表示起見，使用p′和p″分別表示兩種固定資本的價格指數(shù)。根據(jù)生產(chǎn)方程的第三式，以及前面給出的參數(shù)，我們得到部門j使用兩種固定資本分別付出的成本（租金），包括資本增長率及折舊，如下式所示：

基于上式，再根據(jù)工業(yè)品及建筑業(yè)品分別占固定資本的構(gòu)成比例，部門j每單位的固定資本投資需要付出的成本（租金）可表示為：

容易理解的是，當(dāng)市場達(dá)到平衡時，產(chǎn)業(yè)部門使用的各種固定資本其統(tǒng)一回報率相等，即r=r′=r″。則上式轉(zhuǎn)化為：

不難發(fā)現(xiàn)，對于42部門的區(qū)域經(jīng)濟(jì)系統(tǒng)投入產(chǎn)出模型，上述價格動態(tài)方程共有126個方程，而變量卻共有128個變量。按照方程自身的經(jīng)濟(jì)學(xué)特點(diǎn)，變量個數(shù)大于方程個數(shù)，且方程相互之間并非全部線性獨(dú)立，如果按照數(shù)學(xué)上非線性的方程組算法，求解過程將變得十分復(fù)雜，具體方法我們將后續(xù)研究中詳細(xì)解釋。

3.3 供需平衡模塊方程

在構(gòu)建了區(qū)域經(jīng)濟(jì)系統(tǒng)生產(chǎn)與價格動態(tài)方程后，下面主要根據(jù)列昂惕夫動態(tài)I-O系統(tǒng)總供給流量與總需求流量的關(guān)系，通過對區(qū)域經(jīng)濟(jì)系統(tǒng)的供需模塊的動態(tài)平衡進(jìn)行數(shù)學(xué)描述和方程的列寫。

3.3.1 數(shù)學(xué)描述

列昂惕夫動態(tài)投入產(chǎn)出模型屬于動態(tài)一般均衡分析的范疇，在區(qū)域經(jīng)濟(jì)系統(tǒng)的平衡增長中，首先應(yīng)滿足如下關(guān)系：

總產(chǎn)出=中間產(chǎn)品+固定資產(chǎn)投資品+

存貨投資品+最終凈產(chǎn)品（24）

對其進(jìn)行數(shù)學(xué)方程的列寫，上式可表述為：

這里，將區(qū)域經(jīng)濟(jì)系統(tǒng)中的實(shí)物平衡方程表示成了與傳統(tǒng)的列昂惕夫動態(tài)投入產(chǎn)出模型完全一致的形式。其中，δ′=δ″=0，即不考慮固定資本折舊，A為中間投入系數(shù)矩陣，根據(jù)歷年投入產(chǎn)出表可以直接查找得到，B（p，w，r）中的元素是價格p、工資率ω和利潤率r的非線性函數(shù)，反映了區(qū)域經(jīng)濟(jì)生產(chǎn)過程中資本與勞動相互替代的情況，這也是與線性投入產(chǎn)出模型的根本區(qū)別。此外，這些非線性函數(shù)的數(shù)學(xué)表達(dá)式由上一部分與本節(jié)前面部分構(gòu)件的區(qū)域經(jīng)濟(jì)系統(tǒng)三大模塊動態(tài)方程組確定，其中的參數(shù)也可根據(jù)實(shí)際數(shù)據(jù)進(jìn)行參數(shù)估計，因而是可計算的。

為便于對實(shí)物平衡方程進(jìn)行求解，還應(yīng)當(dāng)相應(yīng)對消費(fèi)需求向量c（t）的形式進(jìn)行等價變換。用于最終消費(fèi)需求的是消費(fèi)者收入，即勞動力的工資總額，而工資總額又等于生產(chǎn)過程中由于勞動力投入帶來的增加值。在實(shí)際的消費(fèi)中，假設(shè)消費(fèi)者用工資總額wL購買i中產(chǎn)品的支出比例是固定的，令i產(chǎn)品的價格為pi，購買量為ci，支付比例為gi，g1+g2+…+gi=1，則有如下關(guān)系：

4 基于非線性動態(tài)投入產(chǎn)出的虛擬水貿(mào)易分析框架

建立了區(qū)域經(jīng)濟(jì)系統(tǒng)的非線性動態(tài)投入產(chǎn)出模型之后，接下來，將虛擬水的因素考慮進(jìn)去，建立相應(yīng)的非線性動態(tài)水資源I-O模型。為了更加直觀的表示水資源在各部門之間的投入產(chǎn)出關(guān)系，這里引入一種全新的區(qū)域水資源投入產(chǎn)出表的設(shè)計思路與編制方法。

如表3、表4所示，用于可計算非線性動態(tài)水資源I-O模型求解的投入產(chǎn)出表，采用“母表”與“子表”的形式。其中，“母表”為區(qū)域經(jīng)濟(jì)價值型流量表，其結(jié)構(gòu)與區(qū)域經(jīng)濟(jì)投入產(chǎn)出表的一般形式完全一致，“子表”為區(qū)域水資源流量表，其結(jié)構(gòu)布局與“母表”是一致的，表中的元素與“母表”中的元素一一對應(yīng)，即“子表”中所有數(shù)據(jù)均為“母表”同一位置價值量所對應(yīng)的直接用水量。

為了建立虛擬水貿(mào)易的可計算非線性動態(tài)投入產(chǎn)出分析框架，這里需要厘清區(qū)域經(jīng)濟(jì)系統(tǒng)一般形式的投入產(chǎn)出模型和水資源系統(tǒng)的投入產(chǎn)出模型之間的關(guān)系，重點(diǎn)作出兩點(diǎn)說明。首先，前面建立的區(qū)域經(jīng)濟(jì)系統(tǒng)實(shí)物平衡方程，具有與傳統(tǒng)的列昂惕夫投入產(chǎn)出模型相同的形式，并且與區(qū)域經(jīng)濟(jì)系統(tǒng)中的價值型投入產(chǎn)出表是相互對應(yīng)的，即方程中的元素Y（t）是指社會總產(chǎn)品的價值，c（t）代表區(qū)域經(jīng)濟(jì)系統(tǒng)的最終消費(fèi)需求對應(yīng)的價值量。其次，區(qū)域經(jīng)濟(jì)系統(tǒng)的平衡關(guān)系適用于接下來的水資源投入產(chǎn)出分析，也就是說，利用標(biāo)準(zhǔn)型的可計算非線性動態(tài)投入產(chǎn)出模型，根據(jù)水資源投入產(chǎn)出母表與子表的對應(yīng)關(guān)系，以水資源的流量對模型中的經(jīng)濟(jì)變量進(jìn)行一一替代，就可以建立起可計算非線性動態(tài)水資源I-O模型，據(jù)此可對區(qū)域經(jīng)濟(jì)系統(tǒng)中的虛擬水貿(mào)易進(jìn)行分析測算。對于第二點(diǎn)說明，我們基于前面建立的區(qū)域經(jīng)濟(jì)系統(tǒng)平衡方程給出簡要證明，為便于閱讀，將該方程再度列寫如下：

根據(jù)前面的基本假定，由于產(chǎn)品類型不同及其價值的差別懸殊決定了使用單位質(zhì)量產(chǎn)品生產(chǎn)所需的水資源數(shù)量來表示產(chǎn)品對水資源依賴程度大小，因此使用單位價值產(chǎn)品生產(chǎn)中所需要的水資源數(shù)量來對區(qū)域經(jīng)濟(jì)系統(tǒng)的虛擬水含量進(jìn)行表征。這里，我們將n個產(chǎn)業(yè)部門產(chǎn)品的虛擬水含量分別以vw1，vw2，…vwn來表示，不妨假設(shè)各產(chǎn)品虛擬水含量在一段時間內(nèi)保持不變。

上面的方程中，Y（t）是一個n行列向量，代表了n個產(chǎn)業(yè)部門的產(chǎn)出，因此對于產(chǎn)業(yè)部門i來說，有：

上式即為可計算非線性動態(tài)水資源I-O模型的標(biāo)準(zhǔn)形式?？梢姡捎嬎惴蔷€性動態(tài)水資源I-O模型依然保持了與傳統(tǒng)列昂惕夫投入產(chǎn)出模型基本一致的形式，但兩者的一個重要區(qū)別是，靜態(tài)投入產(chǎn)出模型中的最終需求包括了最終消費(fèi)和投資需求，盡管將固定資產(chǎn)的投資與需求也包含在內(nèi)，但分不清投資結(jié)構(gòu)與消費(fèi)結(jié)構(gòu)，而可計算非線性動態(tài)水資源I-O模型則定義了固定資產(chǎn)使用系數(shù)，將最終需求中的消費(fèi)需求與投資需求分開，從理論上對后者進(jìn)行了完善，更為清晰的表述區(qū)域經(jīng)濟(jì)系統(tǒng)中的投入產(chǎn)出關(guān)系。

與一般性的非線性動態(tài)投入產(chǎn)出模型一樣，非線性動態(tài)水資源投入產(chǎn)出模型考慮了區(qū)域經(jīng)濟(jì)發(fā)展中資本與勞動的相互替代規(guī)律，更為貼近區(qū)域經(jīng)濟(jì)的發(fā)展實(shí)際。式中，A矩陣中的元素與一般性的非線性動態(tài)水資源投入產(chǎn)出模型完全相同，系數(shù)矩陣B和T中的元素仍是價格的函數(shù)，可以依照實(shí)際數(shù)據(jù)進(jìn)行參數(shù)估計，因而是可計算的投入產(chǎn)出模型。需要說明的是，盡管我們在數(shù)學(xué)方程的表達(dá)上對兩個投入產(chǎn)出模型采用了相同的參數(shù)表示形式，但模型中參數(shù)的現(xiàn)實(shí)意義卻并不相同?；蛘哒f，區(qū)域經(jīng)濟(jì)系統(tǒng)和水資源系統(tǒng)的非線性動態(tài)I-O模型擁有相同的數(shù)學(xué)表達(dá)形式，但參數(shù)的量綱是明顯存在區(qū)別的。特別的，對于可計算非線性動態(tài)水資源I-O模型而言，盡管這里對其參數(shù)的現(xiàn)實(shí)意義難以立刻給出合理的解釋，但這并不影響模型的計算分析。至此，我們已經(jīng)將可計算非線性動態(tài)產(chǎn)出模型從一般形式擴(kuò)展至水資源領(lǐng)域，根據(jù)該模型，可以對區(qū)域經(jīng)濟(jì)系統(tǒng)中各產(chǎn)業(yè)部門間的虛擬水流動，進(jìn)行更加客觀、精確的測算。

5 小 結(jié)

本文在傳統(tǒng)靜態(tài)水資源投入產(chǎn)出模型的基礎(chǔ)上，引入可計算非線性動態(tài)投入產(chǎn)出理論，通過區(qū)域經(jīng)濟(jì)系統(tǒng)生產(chǎn)、價格以及實(shí)物平衡方程的列寫，闡述了將區(qū)域經(jīng)濟(jì)系統(tǒng)非線性動態(tài)投入產(chǎn)出模型從一般性分析向水資源領(lǐng)域的轉(zhuǎn)換思路，建立了測算區(qū)域經(jīng)濟(jì)系統(tǒng)中虛擬水貿(mào)易的投入產(chǎn)出分析框架。需要說明的是，由于聯(lián)立方程組涉及的變量及參數(shù)較多，大大超過了方程的數(shù)量，使得模型的求解過程變得十分復(fù)雜，限于篇幅限制這里不再進(jìn)行專門闡述。

在后續(xù)工作中，我們將重點(diǎn)討論該模型的參數(shù)估計思路及求解方法，利用庇隆-弗羅賓紐斯定理（P-F定理）、Jordan分解定理及GaussSeidel迭代運(yùn)算思想，對經(jīng)濟(jì)系統(tǒng)平衡增長路徑的存在性和唯一性進(jìn)行分析，對基于可計算非線性動態(tài)I-O模型的虛擬水測算方法展開探討，并通過著名的“快車道”定理，推導(dǎo)最優(yōu)增長路徑下區(qū)域經(jīng)濟(jì)系統(tǒng)產(chǎn)出結(jié)構(gòu)與用水結(jié)構(gòu)的求解思路。同時，結(jié)合虛擬水貿(mào)易對區(qū)域經(jīng)濟(jì)的作用機(jī)理[30]，回答如何基于以上模型計算對最優(yōu)增長路徑下實(shí)施虛擬水貿(mào)易對區(qū)域經(jīng)濟(jì)增長的貢獻(xiàn)。

（編輯：李 琪）

參考文獻(xiàn)（References）

[1]馬靜， 汪黨獻(xiàn)， HOEKSTRA A Y. 虛擬水貿(mào)易與跨流域調(diào)水[J]. 中國水利， 2004（13）： 37-39.[MA Jing，WANG Dangxian，HOEKSTRA A Y. Virtual water trade and interbasin water transfer[J]. China water resources，2004（13）：37-39.]

[2]馬超. 糧食生產(chǎn)重心北移進(jìn)程中的“虛擬水”流動效應(yīng)分析[J]. 水利發(fā)展研究，2014（4）：21-27.[MA Chao. Analysis of the flow effect of virtual water in the shift of the center of gravity of grain production[J]. Water resources development research，2014（4）：21-27.]

[3]ALLAN J A. Fortunately there are substitutes for water otherwise our hydropolitical futures would be impossible[A]// Priorities for water resources allocation and management. London： ODA， 1993：13-26.

[4]ALLAN J A. Virtual water： a strategic resource. Global solutions to regional deficits[J].Ground water， 1998， 36（4）： 545-546.

[5]ALLAN J A. Virtual water-the water， food and trade nexus： useful concept or misleading metaphor？[J]. Water int， 2003， 28 （1）： 106-113.

[6]HERMAN B. Integrated water management： emerging issues and challenges[J]. Agricultural water management， 2000， 45： 217-228.

[7]HOEKSTRA A Y， HUNG P Q. Virtual water trade： a quantification of virtual water flows between nations in relation to international crop trade[A]//Value of water research report series No. 11. Delft， Netherlands：IHE， 2002.

[8]ZIMMER D， RENAULT D. Virtual water in food production and global trade. Review of methodological issues and preliminary results[A]//Hoekstra A Y. Virtual water trade： proceedings of the international expert meeting on virtual water trade， research report series No.12. Delft， Netherlands： IHE，2003.

[9]CHANPAGAIN A K， HOEKSTRA A Y. The water needed to have the Dutch drink tea[A]//Value of Water Research Report Series No.15.Delft，Netherlands： UNESCOIHE Institute for Water Education， 2003.

[10]羅貞禮， 黃璜. 紅三角地區(qū)農(nóng)產(chǎn)品虛擬水的計算分析[J].冰川凍土， 2005， 27（3）： 46-430.[ LUO Zhenli， HUANG Huang. Calculation of virtual water of agricultural products in the Red Triangular Region analysis[J]. Journal of glaciology and geocryology， 2005， 27 （3）： 46-430.]

[11]王紅瑞， 韓兆興， 韓魯杰， 等. 虛擬水理論與方法的研究進(jìn)展[J]. 中國水利水電科學(xué)研究院學(xué)報， 2008， 6（1）： 66-73.[WANG Hongrui， HAN Zhaoxing， HAN Lujie， et al. Research progress of virtual water theory and method[J]. Journal of China institute of water resources and hydropower research， 2008， 6 （1）： 66-73.]

[12]許新宜，楊中文，王紅瑞，等.水資源與環(huán)境投入產(chǎn)出研究進(jìn)展及關(guān)鍵問題[J].干旱區(qū)地理，2013，36（5）：818-830.[ XU Xinyi， YANG Zhongwen， WANG Hongrui， et al. Research progress and key problems of water resources and environmental input output[J]. Arid area geography， 2013，36 （5）： 818-830.]

[13]劉秀麗，陳錫康.投入產(chǎn)出分析在我國九大流域水資源影子價格計算中的應(yīng)用[J].管理評論，2003，15（1）：49-53.[ LIU Xiuli， CHEN Xikang. Application of input output analysis in the calculation of the shadow price of water resources in the nine basins of China[J]. Management review， 2003， 15 （1）： 49-53.]

[14]王雪妮. 基于區(qū)域間投入產(chǎn)出模型的中國虛擬水貿(mào)易格局及趨勢研究[J].管理評論，2014，26（7）：46-54.[WANG Xueni. Research on the pattern and trend of Chinas virtual water trade based on regional inputoutput model[J]. Management review， 2014，26 （7）： 46-54.]

[15]YU Y，HUBACEK K，F(xiàn)ENG K，et al. Assessing regional and global water footprints for the UK[J]. Ecological economics，2010，69（5）：1140-1147.

[16]ANTONELLI M， ROSON R，SARTORI M. Systemic inputoutput computation of green and blue virtual water ‘flows with an illustration for the Mediterranean region[J].Water resource manage，2012，26（14）：4133-4146.

[17]田貴良，許長新.社會產(chǎn)品價格對水價的敏感性——基于虛擬水貿(mào)易框架的分析[J].中國人口·資源與環(huán)境，2011，21（2）：41-48.[ TIAN Guiliang， XU Changxin. Sensitivity of social product price to water price： an analysis based on the framework of virtual water trade[J]. China population， resources and environment， 2011，21 （2）： 41-48.]

[18]GUAN D， HUBACEK K. Assessment of regional trade and virtual water flows in China[J]. Ecological economics，2007，61（1）：159-170.

[19]譚圣林，邱國玉，雄育久.投入產(chǎn)出法在虛擬水消費(fèi)與貿(mào)易應(yīng)用中的新應(yīng)用[J].自然資源學(xué)報，2014，29（2）：355-364.[ TAN Shenglin， QIU Guoyu， XIONG Yujiu. New application of input output method in the application of virtual water consumption and trade[J]. Journal of natural resources， 2014， 29 （2）： 355-364.]

[20]李方一，劉衛(wèi)東，劉紅光.區(qū)域間虛擬水貿(mào)易模型及其在山西省的應(yīng)用[J].資源科學(xué)，2012，34（5）：802-810.[ LI Fangyi， LIU Weidong， LIU Hongguang. The model of inter regional virtual water trade and its application in Shanxi Province[J]. Resources science， 2012，34 （5）：802-810.]

[21]雷玉桃，蔣璐.中國虛擬水貿(mào)易的投入產(chǎn)出分析[J].經(jīng)濟(jì)問題探索，2012（3）：116-120.[ LEI Yutao， JIANG Lu. The input and output of China virtual water trade[J]. Inquiry into economic issues， 2012 （3）： 116-120.]

[22]周姣，史安娜.區(qū)域虛擬水貿(mào)易計算方法及實(shí)證[J].中國人口·資源與環(huán)境，2008，18（4）：184-188.[ ZHOU Jiao， SHI Anna. Regional virtual water trade calculation method and empirical[J]. China population，resources and environment， 2008，18 （4）： 184-188.]

[23]黃曉榮，裴源生，梁川.寧夏虛擬水貿(mào)易計算的投入產(chǎn)出方法[J].水科學(xué)進(jìn)展，2005，16（4）：564-568.[ HUANG Xiaorong， PEI Yuansheng， LIANG Chuan. The input and output method of virtual water trade in Ningxia[J]. Advances in water science， 2005，16 （4）： 564-568.]

[24]馬忠，張繼良.張掖市虛擬水投入產(chǎn)出分析[J].統(tǒng)計研究，2008（5）：65-70.[MA Zhong， ZHANG Jiliang. Analysis on the input output of virtual water in Zhangye city[J]. Statistical research， 2008 （5）： 65-70.]

[25]張金水. 可計算非線性動態(tài)投入產(chǎn)出模型[M]. 北京： 清華大學(xué)出版社， 2000.[ ZHANG Jinshui. Computable nonlinear dynamic input output model[M]. Beijing： Tsinghua University Press， 2000.]

[26]趙永， 王勁峰， 蔡煥杰. 水資源問題的可計算一般均衡模型研究綜述[J]. 水科學(xué)進(jìn)展， 2008， 19（5）： 756-762.[ ZHAO Yong， WANG Jinfeng， CAI Huanjie. Review of CGE models on water resources[J].Advances in water science， 2008， 19 （5）： 756-762.]

[27]PARTRIDGE M D， RICKMAN D S. Regional computable general equilibrium modeling： a survey and critical appraisal[J].International regional science review，1998， 21（3）： 205-248.

[28]趙永， 王勁峰. 經(jīng)濟(jì)分析CGE模型與應(yīng)用[M]. 北京： 中國經(jīng)濟(jì)出版社， 2008.[ ZHAO Yong， WANG Jinfeng. Economic analysis CGE model and application[M]. Beijing： China Economic Publishing House， 2008.]

[29]PETERSEN T W. An introduction to CGEmodelling and an illustrative application to eastern European intergration wth the EU[EB/OL]. http：//www.dreammodel.dk， 1997.

[30]許長新， 馬超， 田貴良，等. 虛擬水貿(mào)易對區(qū)域經(jīng)濟(jì)的作用機(jī)理及貢獻(xiàn)份額研究[J]. 中國軟科學(xué)， 2011（12）： 110-119.[ XU Changxin， MA Chao， TIAN Guiliang， et al. Study on the mechanism of virtual water trade and its contribution to regional economy[J]. China soft science， 2011 （12）： 110-119.]