棉花貿(mào)易中虛擬土地資源進口能否緩解國內(nèi)棉花用地緊張？

曹　沖　張雅茜　馬惠蘭　宋玉蘭

（新疆農(nóng)業(yè)大學(xué)經(jīng)濟與貿(mào)易學(xué)院，新疆烏魯木齊830052）

棉花貿(mào)易中虛擬土地資源進口能否緩解國內(nèi)棉花用地緊張？

曹沖張雅茜馬惠蘭宋玉蘭

（新疆農(nóng)業(yè)大學(xué)經(jīng)濟與貿(mào)易學(xué)院，新疆烏魯木齊830052）

［摘要］文章針對虛擬土地資源進口能否緩解國內(nèi)棉花用地緊張這一現(xiàn)實問題，以虛擬土地資源進口和國內(nèi)棉花用地變動趨勢為研究基礎(chǔ)，運用土地警度判定模型分析國內(nèi)棉花用地緊張程度，運用協(xié)整檢驗和誤差修正模型來說明二者之間的長短期均衡關(guān)系。研究結(jié)果顯示：國內(nèi)棉花用地處于一種緊張狀態(tài)；虛擬土地資源進口是緩解國內(nèi)棉花用地緊張的原因，而國內(nèi)棉花進口不是引起虛擬土地資源進口的原因；在長期內(nèi)，虛擬土地資源進口和國內(nèi)棉花用地之間存在長期均衡關(guān)系，并且虛擬土地資源進口能夠緩解國內(nèi)棉花用地緊張格局；在短期內(nèi)，國內(nèi)棉花用地對虛擬土地資源進口有負的顯著影響，但是虛擬土地資源進口對緩解國內(nèi)棉花用地緊張的影響不顯著。

［關(guān)鍵詞］棉花貿(mào)易；虛擬土地資源進口；緩解；棉花用地

一、引言

面對當前經(jīng)濟快速發(fā)展和糧食對土地資源的剛性需求等因素，國內(nèi)土地資源供求矛盾日益突顯。然而，在農(nóng)產(chǎn)品中占有重要地位的棉花，其用地面積已由棉花流通體制改革（1999年）時的372.56萬公頃增加到2013年的435萬公頃，但是面對日益增長的需求，國內(nèi)棉花供給還是不能滿足國內(nèi)及世界對中國棉紡織品的迫切需要，使得國內(nèi)棉花安全性逐步降低。那么，如何在有限的耕地資源下能夠保證棉花產(chǎn)業(yè)的安全，虛擬土地資源進口或許是一個不錯的新選擇，但是虛擬土地資源進口能否緩解國內(nèi)棉花用地資源緊張到目前為止還是一個值得商榷的問題。

關(guān)于虛擬土地資源進口對土地資源的貢獻，羅貞禮等指出，在耕地資源缺乏的國家，如果國家或者地區(qū)間通過合理的分工和交換，那么進口的糧食和食品是能夠補充國內(nèi)資源的耕地短缺，同時也能夠提高土地資源的利用效率和產(chǎn)出效率［1］。成麗等提出糧食的凈進口對耕地資源的持續(xù)發(fā)展有正的顯著影響，虛擬土地資源進口能夠?qū)χ袊找婢o張的耕地資源起到一定的緩解作用［2］。劉紅梅等和強文麗等研究表明虛擬土地貿(mào)易對中國耕地資源壓力的緩解起到了積極的貢獻［3］［4］。嚴志強等認為虛擬土地貿(mào)易對糧食安全、土地安全、土地資源利用效益和土地資源配置等方面具有重要意義，同時也認為虛擬土地貿(mào)易是糧食安全和土地安全得以實現(xiàn)的強有力工具，是土地資源利用效益得以提高的調(diào)節(jié)工具［5］。閆麗珍等認為虛擬土地資源戰(zhàn)略的實施能夠使耕地資源和糧食安全得到有效的保障［6］。馬博虎等通過虛擬耕地資源貿(mào)易量的計算方法，估計了中國的小麥、大豆、玉米和大米等4種主要糧食作物，認為虛擬耕地戰(zhàn)略有助于使耕地資源管理得到加強，能夠使耕地資源配置得以優(yōu)化，如果合理實施虛擬耕地戰(zhàn)略，能夠使得國內(nèi)的耕地資源壓力得以緩解［7］。關(guān)于土地資源對虛擬土地資源進口的影響，劉紅梅等表明了糧食播種面積對虛擬土地資源進口量有負的顯著影響［3］。針對虛擬土地資源和耕地資源之間的關(guān)系，鄒君等認為農(nóng)產(chǎn)品虛擬耕地與耕地資源之間有密切關(guān)系，有一定波動，但波動幅度不大，且還在一個較為合理的范圍內(nèi)（基尼系數(shù)亦保持在0.3－0.4）［8］。

綜上所述，諸多學(xué)者對耕地資源短缺國家的研究表明虛擬土地資源進口能夠緩解土地資源短缺。然而，虛擬土地資源進口和國內(nèi)棉花用地之間到底是一個怎樣的關(guān)系，目前研究較少，因此本文主要對目前虛擬土地資源進口同國內(nèi)棉花用地之間的長短期均衡關(guān)系進行研究。

二、國內(nèi)棉花用地和虛擬土地資源進口狀況

（一）“先增后減”趨勢的國內(nèi)棉花用地

棉花流通體制改革以來的15年間，國內(nèi)棉花用地面積總體呈緩慢波動增加趨勢，由372.56萬公頃增加到435萬公頃，累計增加62.44萬公頃，增加了16.76%，年均增加1.11%。期間棉花用地面積呈現(xiàn)“波動增加—逐步減少”兩個階段特征：

1.第一階段（1999－2007年）。國內(nèi)棉花用地面積呈現(xiàn)波動增加態(tài)勢。由1999年的372.56萬公頃波動增加到2007年的592.61萬公頃，9年間累計增加220.05萬公頃，增加了59.06%，其中2007年是這一時期棉花用地面積最高的一年，為592.61萬公頃，此階段棉花用地面積隨著時間變動的線性回歸方程為：Y＝26.897T＋358.4906

①1999－2007年的時間序列把1999年作為1，2000年作為2，以此類推。②2008－2013年的時間序列把2008年作為1，2009年作為2，以此類推。（R2＝0.831573）①。

2.第二階段（2008－2013年）。國內(nèi)棉花用地面積逐步減少。由2008年的575.41萬公頃減少到2013年的435萬公頃，6年間累計減少140.41萬公頃，減少了24.40%，年均減少5.44%，其中2013年是這一階段棉花用地面積最少的年份，為435萬公頃，此階段棉花用地面積隨著時間變動的線性回歸方程為：Y＝－21.752T＋569.922 （R2＝0.758070）②。

（二）“增加—減少—回升”趨勢的虛擬土地資源進口

1999年以來，中國棉花貿(mào)易中的虛擬土地資源進口總體呈波動增加趨勢，由1999年的42.14萬公頃增加到2013年的273.13萬公頃，15年間累計增加230.99萬公頃，增加了5.48倍，年均增加14.28%，番了2翻有余，其虛擬土地資源進口隨時間的變動線性回歸方程為：Y＝18.36879T＋32.90971（R2＝0.664117）。期間中國棉花貿(mào)易中的虛擬土地資源進口呈現(xiàn)“增加—減少—回升”三時期特征。

1.第一時期（1999－2006年）。棉花貿(mào)易中虛擬土地資源進口快速增加。由1999年的42.14萬公頃增加到2006年的309.13萬公頃，8年間虛擬土地資源進口累計增加266.99萬公頃，平均每年增加33.37萬公頃，翻了三翻，其中2006年是虛擬土地資源進口的歷史最高點，為309.13萬公頃。

2.第二時期（2007－2009年）。棉花貿(mào)易中虛擬土地資源進口逐步減少。2007年虛擬土地資源進口為212.96萬公頃，2008年減少為177.55萬公頃，到2009年進一步減少為140.01萬公頃，3年時間虛擬土地資源進口累計減少72.95萬公頃。

3.第三時期（2010－2013年）。棉花貿(mào)易中虛擬土地資源進口逐步回升，由2010年的254.05萬公頃增加到2013年的273.13萬公頃，累計增加19.08萬公頃，增加了7.51%，其中2012年是這一時期的新高，為304.13萬公頃，2013年低于2012年大約10個百分點，為273.13萬公頃。

（三）“異步—同步—異步”趨勢的虛擬土地資源進口與國內(nèi)棉花用地

同期，虛擬土地資源進口同國內(nèi)棉花用地之間的變動趨勢在時間序列上存在“異步—同步—異步”三個階段特征（見圖1）。

圖1　虛擬土地資源進口量與國內(nèi)棉花用地變動趨勢

1.“異步”階段（1999－2001年）。虛擬土地資源進口逐年減少，由42.14萬公頃減少到38.29萬公頃；國內(nèi)棉花用地面積逐年增加，由372.56萬公頃增加到480.98萬公頃。此階段出現(xiàn)“異步”的原因很有可能是由于棉花流通體制改革處于初期，在市場機制下的宏觀調(diào)控對棉花資源合理配置還沒有發(fā)揮出其應(yīng)有的作用。

2.“同步”階段（2002－2006年）。虛擬土地資源進口和國內(nèi)棉花用地資源二者同步增加。虛擬土地資源進口由51.96萬公頃波動增加到309.13萬公頃，累計增加257.17萬公頃，其中2006年是棉花貿(mào)易中虛擬土地資源進口最多的一年；國內(nèi)棉花用地面積由418.42萬公頃增加到581.57萬公頃，累計增加163.15萬公頃。此階段虛擬土地資源進口和國內(nèi)棉花用地面積不斷增加，這一方面主要歸因于棉花流通體制改革在市場機制作用下逐步發(fā)揮作用，對資源配置產(chǎn)生影響；另一方面因為中國加入世界貿(mào)易組織，充分利用“兩種資源”和“兩個市場”為國內(nèi)棉花生產(chǎn)服務(wù)；同時世界棉紡織品配額的取消也是刺激棉花需求不斷高漲的主要因素，這必然促使虛擬土地資源進口和國內(nèi)棉花用地面積不斷增加。

3.“異步”階段（2007－2013年）。虛擬土地資源進口波動增加，國內(nèi)棉花用地面積逐步減少。虛擬土地資源進口由212.96萬公頃增加到273.13萬公頃，累計增加60.17萬公頃，增加了28.26%；國內(nèi)棉花用地面積由592.61萬公頃波動減少到435萬公頃，累計減少157.61萬公頃，減少了26.60%。此階段的“異步”趨勢主要是因為世界經(jīng)濟危機的波及，降低了世界市場對原棉產(chǎn)品的需求，再加上國內(nèi)棉花成本不斷上升，棉農(nóng)會根據(jù)產(chǎn)品的比較收益生產(chǎn)收益較高的產(chǎn)品，致使國內(nèi)棉花用地面積不斷減少。同時中國作為一個紡織品制造大國，為了滿足國內(nèi)和國外兩個市場的需求，在國內(nèi)棉花價格高居的情況下，國外低價棉花就成為企業(yè)選擇對象，那么虛擬土地資源進口也相對增加。

三、國內(nèi)棉花用地緊張程度判定

根據(jù)葛向東等提出的土地資源警度判定模型［9］，嘗試用來測算國內(nèi)棉花用地的緊張程度，模型如下：

式（1）中，A為土地警度判定系數(shù)；為土地資源下所產(chǎn)的實際增量，其中Ti和Tj分別為預(yù)警年份和基準年；△T＝△X·N為計劃增量，△X＝Xi·Xj為基準年人口增加量，N為基準年的人均產(chǎn)量。土地資源安全警度判定如表1所示。

表1　土地資源安全警度判定

根據(jù)（1）式和國內(nèi)棉花改革體制改革的起始點，國內(nèi)棉花用地安全度的計算基準年份定為1999年。由表2國內(nèi)棉花用地計算結(jié)果可以看出：國內(nèi)棉花用地均處于中警以上水平（14年中1年為重警、6年為巨警和7年為中警），其中，2003年國內(nèi)棉花用地為重警，2002年、2005年、2008年、2009年、2010年和2013年這6年為巨警級別，其他年份為中警級別。可見，隨著國內(nèi)人均棉花用地面積和后備土地資源數(shù)量的不斷減少，再加上全球環(huán)境惡化和自然災(zāi)害的嚴重加劇，國內(nèi)棉花用地安全問題會頻繁突出，將引發(fā)更加緊張的局勢，這充分說明了國內(nèi)棉花用地處于一種不安全狀態(tài)即緊張狀態(tài)。

四、虛擬土地資源進口與國內(nèi)棉花用地之間的協(xié)整關(guān)系及誤差修正

（一）虛擬土地資源進口與國內(nèi)棉花用地數(shù)據(jù)的平穩(wěn)性檢驗

為了減緩對非平穩(wěn)序列直接進行回歸可能產(chǎn)生的“偽回歸”問題，一般采用ADF檢驗方法。在ADF檢驗中，為了保證方程中的是白噪聲，在方程的右邊加入一些滯后項。于是單位根檢驗方程如下所示：

模型1：

模型2（模型1＋常數(shù)項）：

模型3（模型2＋時間趨勢項）：

在實際檢驗時，首先從模型3開始，然后是模型2，最后是模型1，什么時候拒絕零假設(shè)，什么時候檢驗停止。這里需要注意的是，零假設(shè)為H0：p＝1。否則就要繼續(xù)進行檢驗。只有通過對模型3、模型2和模型1的檢驗，才能得到各自相應(yīng)的臨界值。

模型中用于分析虛擬土地資源進口與國內(nèi)棉花用地之間關(guān)系的變量分別標記為V和D，為了克服時間序列中存在的自相關(guān)現(xiàn)象，對原始變量采用自然對數(shù)轉(zhuǎn)換，分別記為lnV和lnD。同時，借助Eviews 7.2軟件，對變量lnV和lnD進行平穩(wěn)性檢驗，檢驗結(jié)果由表3可知：變量lnV和lnD的水平檢驗結(jié)果均存在單位根，具有不穩(wěn)定性，而一階差分檢驗結(jié)果均不存在單位根，具備穩(wěn)定性，并且變量lnV在1%顯著水平上拒絕零假設(shè)，lnD在5%顯著水平上拒絕零假設(shè)。

表2　國內(nèi)棉花用地警度變動

表3　lnV和lnD單位根檢驗結(jié)果

表4　lnV和lnD的因果關(guān)系檢驗

（二）虛擬土地資源進口與國內(nèi)棉花用地之間的因果關(guān)系檢驗

由于對變量的數(shù)值變換并不會改變變量之間的因果關(guān)系，因此利用變量對數(shù)值對虛擬土地資源進口和國內(nèi)棉花用地之間進行Granger因果關(guān)系檢驗，結(jié)果如表4所示?？梢钥闯?，虛擬土地資源進口不是國內(nèi)棉花用地的Grange原因的F值為1.36165，P值為0.3550，在10%的顯著水平下不顯著，故國內(nèi)棉花用地不是引起虛擬土地資源進口的原因；而國內(nèi)棉花用地緊張不是虛擬土地資源進口的Grange原因的F值為8.84863，P值為0.0201，在5%的顯著水平下顯著，故虛擬土地資源進口是緩解國內(nèi)棉花用地緊張的原因?？梢?，虛擬土地資源進口是有效緩解國內(nèi)棉花用地緊張的一個原因，而國內(nèi)棉花用地緊張并不是引起虛擬土地資源進口的原因，二者并不存在明顯的互饋關(guān)系。

（三）虛擬土地資源進口與國內(nèi)棉花用地之間的長期均衡模型

協(xié)整檢驗主要包括兩種：一種是基于回歸系數(shù)的協(xié)整檢驗，如Johansen協(xié)整檢驗；另一種是基于回歸殘差的協(xié)整檢驗，如CRDW檢驗、DF檢驗和ADF檢驗。這里主要運用Engle和Granger （1987）提出的協(xié)整檢驗方法（回歸殘差的協(xié)整檢驗）。檢驗的主要步驟如下：

（1）若k個序列y1t和y2t，y3t，…，ykt都是1階單整序列，建立回歸方程：

模型估計的殘差為：

Engle和Granger（1987）提出，如果兩個或者兩個以上變量之間呈現(xiàn)平穩(wěn)性，或者它們的一階差分具有相同的階數(shù)，并且也都具有平穩(wěn)性，那么變量之間可能具有協(xié)整關(guān)系。如果變量之間不存在平穩(wěn)性，或者在顯著水平下不存在相同的單整階數(shù)，變量之間的協(xié)整關(guān)系就很難確定。因此，根據(jù)上文的單位根檢驗結(jié)果可以看出，變量lnV 和lnD在5%的顯著水平下，單整階數(shù)相同，具有平穩(wěn)性?？梢?，二者之間很有可能存在協(xié)整關(guān)系。

根據(jù)協(xié)整理論，得出反映二者關(guān)系的兩個不同的方程：lnD＝1.22185lnV和lnD＝5.650834＋0.109173lnV。然而，根據(jù)本文的研究目標，認為后者較為正確，主要是因為前者在構(gòu)造變量之間的長期均衡關(guān)系時沒有把具有顯著影響的常數(shù)項包含進去，同時后者也是對前者誤差的一種修訂，并且對殘差進行單位根檢驗時，檢驗形式中也包含了常數(shù)項和時間趨勢項。因此，本文估計變量之間的長期均衡穩(wěn)定關(guān)系建立的回歸方程如下所示：

根據(jù)公式（7），可以看出變量之間存在協(xié)整關(guān)系，但是二者是否真正存在協(xié)整，需要進一步對殘差進行檢驗。

通過對回歸方程的殘差進行單位根檢驗，結(jié)果顯示（見表5）：上述方程的殘差在5%的顯著水平拒絕原假設(shè)，說明其是平穩(wěn)序列。因此，可以斷定虛擬土地資源進口和國內(nèi)棉花用地面積之間存在協(xié)整關(guān)系，同時也進一步說明了上述選擇的模型是正確的。

綜上所述，虛擬土地資源進口同國內(nèi)棉花用地面積之間存在一個長期的均衡穩(wěn)定關(guān)系。在其他條件不變的情況下，當虛擬土地資源進口量每增加1%時，將會緩解國內(nèi)棉花用地需求0.11%。同時也應(yīng)看到，該方程回歸系數(shù)的符號和大小與經(jīng)濟理論的期望值相符合。

（四）虛擬土地資源進口與國內(nèi)棉花用地之間的誤差修正

根據(jù)上述協(xié)整方程得到虛擬土地資源進口與國內(nèi)棉花用地之間存在長期的均衡穩(wěn)定關(guān)系，但是短期內(nèi)的均衡關(guān)系還不知曉，為此需要建立反映二者關(guān)系的短期均衡模型，所以嘗試用誤差修正模型來測算二者是否存在短期均衡，同時也衡量虛擬土地資源進口對國內(nèi)棉花用地緊張的緩解程度。

表5　殘差平穩(wěn)性的ADF檢驗

表6　誤差修正模型估計結(jié)果

根據(jù)Grange（1969）提出的涉及變量短期關(guān)系的因果檢驗方法，文中采用Engle和Grange （1987）提出的誤差修正模型。該模型的基本方程如下所示：

式中，αD和αV分別是調(diào)整系數(shù)，△lnDt和△lnVt分別表示虛擬土地資源進口和國內(nèi)棉花用地的滯后項。當αD和β1i等于0時，那么△lnVt不是引起△lnDt的Grange原因；同理，當αV和β2i等于0時，那么△lnDt不是引起△lnVt的Grange原因。

由表6的誤差修正模型可知，根據(jù)公式（8）所得出的結(jié)果是△lnVt－1和△lnVt－2在10%的顯著水平下對△lnDt不顯著，但是其誤差修正項系數(shù)為負值，卻符合反向修正機制。根據(jù)公式（9）所得出的結(jié)果是誤差修正項的系數(shù)在1%的顯著水平下為－0.685，符合誤差修正模型的反向修正機制，并且根據(jù)△lnDt－1和△lnDt－2對△lnVt的參數(shù)分別在1%和5%的顯著水平下顯著，并且系數(shù)為負值，說明了在短期內(nèi)國內(nèi)棉花用地對虛擬土地資源進口產(chǎn)生負的影響，并且影響較大，這與經(jīng)濟理論的假設(shè)和現(xiàn)實經(jīng)濟情況較為符合。

五、結(jié)論及啟示

本文在對棉花流通體制改革以來的國內(nèi)棉花用地和虛擬土地資源進口狀況分析的基礎(chǔ)上，通過警度判定模型說明國內(nèi)棉花安全警度級別，運用協(xié)整關(guān)系和誤差修正模型來驗證二者之間的長短期均衡關(guān)系，進而得出以下幾點結(jié)論：（1）國內(nèi)棉花用地同虛擬土地資源進口存在“異步—同步—異步”關(guān)系；（2）國內(nèi)棉花用地處于一種不安全狀態(tài)，即緊張狀態(tài)；（3）虛擬土地資源進口是緩解國內(nèi)棉花用地緊張的原因，而國內(nèi)棉花用地緊張不是引起虛擬土地資源進口的原因；（4）虛擬土地資源進口同國內(nèi)棉花用地存在長期穩(wěn)定均衡關(guān)系；（5）在短期內(nèi)，虛擬土地資源進口對緩解國內(nèi)棉花用地資源緊張不顯著，但是國內(nèi)棉花用地卻對虛擬土地資源進口產(chǎn)生負的顯著影響。

通過以上結(jié)論發(fā)現(xiàn)，虛擬土地資源進口在長期內(nèi)可以緩解國內(nèi)棉花用地緊張，可以保障國內(nèi)棉花產(chǎn)業(yè)安全，但是在短期內(nèi)如何把握棉花安全，可能會出現(xiàn)一系列問題。如果能夠兼顧國內(nèi)和國外政策等因素的一致性，適度運用棉花貿(mào)易中的虛擬土地資源戰(zhàn)略可能會使國內(nèi)棉花安全得到保障、資源配置得到優(yōu)化、產(chǎn)業(yè)結(jié)構(gòu)得到調(diào)整。

參考文獻

［1］羅貞禮，龍愛華，等.虛擬土戰(zhàn)略與土地資源可持續(xù)利用的社會化管理［J］.冰川凍土，2004，（05）：624－631.

［2］成麗，吳迪，等.糧食貿(mào)易對耕地資源可持續(xù)利用影響分析［J］.農(nóng)業(yè)經(jīng)濟，2011，（01）：41－43.

［3］嚴志強，顏章雄，等.虛擬土地、虛擬土地戰(zhàn)略與區(qū)域土地資源優(yōu)化配置管理的理論探討［J］.廣西社會科學(xué)，2007，（10）：70－74.

［4］劉紅梅，王克強，等.中國糧食虛擬土地資源進口的實證分析［J］.中國農(nóng)村經(jīng)濟，2007，（11）：26－33＋51.

［5］強文麗，劉愛民，等.中國農(nóng)產(chǎn)品貿(mào)易的虛擬土地資源量化研究［J］.自然資源學(xué)報，2013，（08）：1289－1297.

［6］閆麗珍，成升魁，等.玉米南運的虛擬耕地資源流動及其影響分析［J］.中國科學(xué)院研究生院學(xué)報，2006，（03）：342－348.

［7］馬博虎，張寶文.中國糧食對外貿(mào)易中虛擬耕地貿(mào)易量的估算與貢獻分析——基于1978－2008年中國糧食對外貿(mào)易數(shù)據(jù)的實證分析［J］.西北農(nóng)林科技大學(xué)學(xué)報（自然科學(xué)版），2010，（06）：115－119，126.

［8］鄒君，劉文潔，等.中國農(nóng)產(chǎn)品虛擬耕地與資源環(huán)境經(jīng)濟要素的時空匹配分析［J］.長江流域資源與環(huán)境，2012，（04）：477－481.

［9］葛向東，彭補拙，等.耕地總量動態(tài)平衡的監(jiān)測和預(yù)警研究［J］.自然資源學(xué)報，2002，（01）：35－41.

（責任編輯：胡宜挺）

［作者簡介］曹沖（1988－），男，安徽臨泉人，碩士研究生，研究方向：貿(mào)易經(jīng)濟和區(qū)域經(jīng)濟；張雅茜（1990－），女，山西運城人，碩士研究生，研究方向：農(nóng)業(yè)經(jīng)濟管理；宋玉蘭（1979－），女，山西運城人，博士，碩士生導(dǎo)師，副教授，研究方向：農(nóng)業(yè)經(jīng)濟管理；通訊作者：馬惠蘭（1962－），女，博士，博士生導(dǎo)師，教授，研究方向：區(qū)域經(jīng)濟和農(nóng)業(yè)經(jīng)濟。

［基金項目］新疆維吾爾自治區(qū)產(chǎn)學(xué)研聯(lián)合培養(yǎng)研究生示范基地資助項目（xjaucxy－yjs－20131038）；干旱區(qū)農(nóng)村發(fā)展研究中心資助項目（XJEDU030113Y10）。