夏 浩，苑韶峰*，楊麗霞，王 瑤

（1. 浙江工商大學公共管理學院，浙江?杭州 310018；2. 浙江財經大學城鄉(xiāng)規(guī)劃與管理學院，浙江?杭州 310018）

基于空間自相關分析的浙南城市群建設用地利用效益評價與分區(qū)管制

夏 浩1，苑韶峰1*，楊麗霞2，王 瑤1

（1. 浙江工商大學公共管理學院，浙江?杭州 310018；2. 浙江財經大學城鄉(xiāng)規(guī)劃與管理學院，浙江?杭州 310018）

結合建設用地綜合效益評價和空間自相關分析，為浙西南城市群建設用地分區(qū)管制和精細化管理提供依據(jù)。分析結果表明：浙南城市群各縣市之間的建設用地利用綜合效益差異較大，且存在明顯的空間溢出作用；浙南城市群建設用地效益呈現(xiàn)區(qū)域集聚性；建設用地管制區(qū)可劃分為調整優(yōu)化區(qū)、協(xié)調發(fā)展區(qū)、重點發(fā)展區(qū)、適度發(fā)展區(qū)和生態(tài)保護區(qū)。研究認為：建設用地分區(qū)管制要依據(jù)建設用地效益以及區(qū)域發(fā)展特征，并制定差異化政策。

建設用地；效益評價；分區(qū)管制；空間自相關分析；浙南城市群

城市化推動城市建設用地面積不斷擴張，但單一偏重建設用地擴張的城鎮(zhèn)化模式使得土地資源粗放利用嚴重[1]，土地利用沖突明顯[2]，建設用地水土污染問題突出[3]。因此有必要對城市建設用地效益進行綜合評價，反映城市建設用地利用效率的高低，為城市建設用地的發(fā)展進行管制和調控提供有益借鑒。

當前城市建設用地效益評價涉及村鎮(zhèn)[4]、城市[5,6]、省域[7,8]、區(qū)域[9]和國家[10,11]等多個尺度，而鮮有針對城市群建設用地效益評價的研究。目前建設用地效益評價研究已涉及經濟效益、生態(tài)效益、集約水平、綜合效益等多方面的評價內容，如王筱明和鄭新奇以二、三產業(yè)產值與建成區(qū)面積之比作為衡量城市用地效益的指標，對濟南市城市用地效益和擴張合理性進行分析[12]，吳振華等人對江浙滬地區(qū)城市建設用地生態(tài)效率、空間格局演變進行了實證分析[13]，任飛陽和韋仕川對三亞市建設用地集約利用水平進行定量評價[14]，趙勝男等人對西安市建設用地擴張綜合效益進行評價[15]。建設用地效益評價結果多與建設用地效益區(qū)域差異研究[11,16]、城市合理用地規(guī)模研究[12]以及建設用變化特征[9]等研究相關聯(lián)，而將評價結果用于分區(qū)管制的研究較少。分區(qū)管制對土地利用產生積極效果，能有效控制建設用地規(guī)模，也有利于實現(xiàn)人口、經濟、資源和環(huán)境*通訊作者: 苑韶峰(博士后/教授):shaofengyuan1975@163.com的空間均衡[17]。

鑒于此，本文以浙南城市群為研究區(qū)域，從縣域視角出發(fā)，通過建設用地利用綜合效益的評價和空間自相關分析，對區(qū)域建設用地進行更精細的分區(qū)管制。

1 區(qū)域概況及數(shù)據(jù)來源

1.1 區(qū)域概況

浙南城市群包括溫州市、麗水市和臺州市三市，下轄25個縣（市）。2014年地區(qū)總面積為38503km2，其中建設用地面積約為1272.6km2。浙南地區(qū)民營經濟活躍，“溫州模式”和“臺州模式”指導下區(qū)域社會經濟發(fā)展良好，2005~2014年區(qū)域GDP從 3154.11億元增至8742.18億元，總人口從1561.52萬增長到1676.44萬。

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

研究以浙南25個縣（市、區(qū)）為評價單元，各縣市數(shù)據(jù)主要來源于《2015年浙江省統(tǒng)計年鑒》、《2015年溫州市統(tǒng)計年鑒》、《2015年臺州市統(tǒng)計年鑒》以及《2015年麗水市統(tǒng)計年鑒》，部分建設用地面積缺失的縣市數(shù)據(jù)，通過查詢當?shù)貒敛块T公布的建設用地面積進行填補。

2 建設用地利用效益評價

2.1 評價指標體系構建

建設用地利用效益的評價是對土地的一個綜合性評價，包含土地的經濟效益、社會效益以及環(huán)境效益。評價指標應當要體現(xiàn)出建設用地的利用特征，并且考慮其功能和產生的效益。本文遵循科學性、層次性、可比性、目的性等原則，并結合實際情況以及指標數(shù)據(jù)的可獲得性，借鑒相關研究成果[7,9,11,16]，最終構建包含經濟效益、社會效益與環(huán)境效益三個層次的共13個指標。

2.2 數(shù)據(jù)計算

表1 建設用地效益評價指標體系Table 1 The evaluation index system of construction land benefit

（1）指標標準化

由于原始數(shù)據(jù)量綱不同，因此對數(shù)據(jù)進行無量綱化處理：

式中：Xij是指標標準化值；Zij為某一指標屬性值，maxZij、minZij分別為某一指標的最大值和最小值。

（2）權重確定與建設用地利用效益計算

本研究采用改進熵值法[18]確定各指標權重。用改進熵值法確定指標權重，能夠減少主觀賦權法無法避免的臆斷性問題，并可有效解決多指標變量間信息重疊問題，以減少極端值對評價結果的影響，使得到的熵值非常接近原始熵值。采用改進熵值法得到各項指標權重，進而計算各部分效益水平，公式如下：

式中，wij為各指標權重，Ui為各縣市的單項效益。

式中，Wj為準則層的權重，E為各縣市的綜合效益。

2.3 空間分析法

（1）變異系數(shù)

根據(jù)建設用地利用效益分布情況，采取變異系數(shù)進行空間分析。變異系數(shù)用來衡量離散程度的大小，其公式為：

式中，V為變異系數(shù)；E為建設用地利用效益的平均值；m為區(qū)域內城市個數(shù)；Ei區(qū)域內第i個城市的建設用地利用效益。

（2）空間自相關分析

空間自相關反映區(qū)域建設用地效益與鄰近區(qū)域建設用地效益的關聯(lián)程度。本研究采用Moran’s I來衡量全局空間自相關，公式為：

式中，I為全局空間自相關指標；Sij2為各縣（市）建設用地效益的方差；n是空間單元數(shù)目；Di、Dj表示空間單元和 的建設用地效益；Wij是空間權重系數(shù)矩陣。

采用LISA值局部空間自相關相關程度，公式為：

式中，Ii為第i個縣（市）的局部空間自相關指標。

3 建設用地利用效益的空間差異及分區(qū)管制

3.1 建設用地利用效益空間差異

通過計算得到各縣市建設用地利用效益，結合表2和圖1可知，浙南建設用地綜合效益和各分項效益均有顯著的區(qū)域差異。

表2 各市建設用地利用效益Table 2 Construction land use benefit of cities

就綜合效益而言，溫州市綜合效益最高，麗水市次之，臺州市綜合效益最低。從空間上來看，綜合效益高值區(qū)位于東南部溫州市的文成縣、瑞安市、平陽縣、蒼南縣、永嘉縣和樂清市。低值區(qū)位于北部臺州市的天臺縣、三門縣、臨海市、溫嶺市等縣市。溫州市各縣市間綜合效益變異系數(shù)達0.448，表明各縣市之間綜合效益差異顯著，而臺州和麗水兩市變異系數(shù)較小，縣市間綜合效益差異較小。

經濟效益衡量單位建設用地面積上經濟投入和產出狀況[19]，經濟效益指數(shù)以東南部永嘉縣、平陽縣、樂清市、瑞安市等縣市較高，該地區(qū)二三產業(yè)相對發(fā)達，地均二三產業(yè)產值大，土地集約利用水平高。而北部仙居縣、天臺縣和三門縣等縣市建設用地利用相對粗放，地均經濟投入產出小，經濟效益較低。溫州市各縣市間經濟效益變異系數(shù)最大，臺州市次之，麗水市最小，反映出溫州市各縣市間經濟效益區(qū)域差異較顯著，而臺州市和麗水市相對差異較小。

社會效益反映建設用地社會承載力和提供公共服務的能力，永嘉縣、文成縣和平陽縣三市社會效益最高，而仙居縣、玉環(huán)縣、三門縣和天臺縣等縣市社會效益較低。社會效益市域內差異以溫州市最大，麗水市次之，臺州市最小。

環(huán)境效益反映建設用地利用對城市生態(tài)環(huán)境影響，環(huán)境效益最高值位于西部麗水市龍泉縣，而次高值區(qū)域由溫州市各縣市和臺州市區(qū)在東部地區(qū)形成帶狀集聚，低值區(qū)位于臺州市仙居縣、三門縣、天臺縣、玉環(huán)縣和溫嶺縣。環(huán)境效益市域內部差異以麗水和臺州較大，而溫州具有顯著較小的差異性。

圖1 浙南城市群建設用地效益空間分異Fig.1 The spatial differentiation of construction land benefit in the urban agglomeration of south Zhejiang

3.2 建設用地利用效益空間關聯(lián)性分析

利用GeoDA軟件進行空間分析，從而進一步分析建設用地利用效益的空間關聯(lián)性并繪制出Moran散點圖。Moran散點圖由四個象限構成，每一個象限代表不同的空間自相關類型[19]：第Ⅰ象限和第Ⅲ象限代表正相關，即該區(qū)域建設用地效益值較高（低），相鄰區(qū)域建設用地效益值也較高（低），第II象限和第Ⅳ象限代表負相關，即該區(qū)域建設用地效益值較低（高），相鄰區(qū)域建設用地效益值較高（低）。

根據(jù)公式計算得到浙南建設用地利用綜合效益的全局空間自相關指數(shù)Moran’s I=0.25155 (顯著度為0.05)，表明浙南建設用地利用綜合效益存在較顯著的正向空間自相關關系，各縣市建設用地利用綜合效益表現(xiàn)出空間集聚。通過如Moran圖所示，15個縣（市）位于第Ⅰ象限和第Ⅲ象限，占城市總數(shù)的60%，且位于第Ⅲ象限的縣（市）多于第Ⅰ象限的縣（市）。

通過LISA圖進一步進行區(qū)域集聚的分析，處于顯著HH類型的縣市為平陽縣，處于顯著LL類型的縣市為縉云縣、臨海市和洞頭縣，處于LH類型的縣市為溫州市區(qū)，其他城市均為不顯著類型。從以上Moran圖與LISA圖的分析中反映出浙南各縣市之間的建設用地利用綜合效益存在明顯的空間溢出作用，且土地利用經濟效益“兩極化”空間顯著。

3.3 基于利用效益空間分異的建設用地分區(qū)管制

從縣域角度對各縣市進行劃分，從而進行精細化的分區(qū)管制。采用幾何間隔分類法將25個縣的建設用地利用綜合效益水平分成高效、中效和低效三類。再結合Moran圖中的空間關聯(lián)性類型，對浙西南城市群建設用地進行分區(qū)管制。

圖2 Moran散點圖Fig.2 Moran scatter plots

圖3 建設用地利用效益LISA集聚圖Fig.3 LISA gather diagram of construction land use benefit

（1）調整優(yōu)化區(qū)：建設用地總面積為77.4km2，占浙南建設用地總面積的6.08%。這一類型地區(qū)的建設用地利用綜合效益水平高，同時對周邊地區(qū)建設用地利用綜合效益具有溢出效應。因這一類型地區(qū)以及周邊地區(qū)建設用地利用綜合效益都已很高，因此該類地區(qū)對周邊地區(qū)的溢出效應邊際不斷減少。管制措施：①優(yōu)化建設用地結構和布局；②加強與低效益地區(qū)的產業(yè)聯(lián)系，可以將部分產業(yè)遷移至協(xié)同發(fā)展區(qū)與重點發(fā)展區(qū)進而置換部分產業(yè)。

（2）協(xié)同發(fā)展區(qū)：建設用地總面積為239.57km2，占浙南建設用地總面積的18.83%。這一類型地區(qū)的建設用地利用綜合效益水平從浙南整體上來看較高，但與周邊地區(qū)相比較低，形成被包圍的一個低值區(qū)，受周邊地區(qū)溢出效益影響不顯著。管制措施：①推動區(qū)域經濟、社會和環(huán)境效益的協(xié)調發(fā)展，同時優(yōu)化區(qū)域土地利用方式，促進建設用地的集約和節(jié)約利用；②加強與周邊區(qū)域在產業(yè)結構上的連接，以強化周邊地區(qū)對本區(qū)域的溢出作用的影響，促進區(qū)域社會經濟的協(xié)調發(fā)展。

表3 浙南城市群建設用地利用管制分區(qū)Table 3 Regulated zoning of construction land use in the urban agglomeration of south Zhejiang

圖4 浙南城市群建設用地分區(qū)管制圖Fig.4 Map of regulated zoning of construction land use in the urban agglomeration of south Zhejiang

（3）重點發(fā)展區(qū)：建設用地總面積為202.05km2，占浙南建設用地總面積的15.88%。這一類型的地區(qū)建設用地利用綜合效益水平較高，受周邊地區(qū)帶動作用小，同時自身對周邊地區(qū)的拉動作用也不顯著。管制措施：①推動本地區(qū)產業(yè)結構的優(yōu)化和升級，加快轉變經濟發(fā)展方式，加快“退二進三”步伐；②加強與協(xié)同發(fā)展區(qū)與調整優(yōu)化區(qū)之間的聯(lián)系，以加強自身極化作用。

（4）適度發(fā)展區(qū)：建設用地總面積為297.44km2，占浙南建設用地總面積的23.37%。這類地區(qū)建設用地利用綜合效益低，但周邊地區(qū)綜合效益較高。管制措施：①承接部分調整優(yōu)化去與協(xié)同發(fā)展區(qū)的產業(yè)轉移，利用自身優(yōu)勢發(fā)展特色經濟；②在發(fā)展過程中要協(xié)調好經濟、社會和環(huán)境三者之間的關系，促進三者的耦合協(xié)調發(fā)展，以促進地區(qū)綜合效益水平的提高。

（5）生態(tài)保護區(qū)：建設用地總面積為456.14km2，占浙南建設用地總面積的35.84%。這類地區(qū)建設用地利用綜合效益低，周邊地區(qū)綜合效益也低，周邊地區(qū)基本無溢出效應。其中，西南區(qū)塊的松陽縣、云和縣和麗水市區(qū)屬于浙西南山地生態(tài)環(huán)境保護區(qū)；東北區(qū)塊的天臺縣屬于丘陵盆地生態(tài)環(huán)境治理區(qū)，三門縣和臨海市屬于浙東南沿海港灣平原生態(tài)環(huán)境治理區(qū)，玉環(huán)縣屬于浙東海洋島嶼生態(tài)環(huán)境保護區(qū)。這一類型地區(qū)建設用地因生態(tài)敏感性高，建設用地利用的綜合效益低。管制措施：①以保護地區(qū)生態(tài)環(huán)境為重點，在維持生態(tài)系統(tǒng)平衡的基礎上發(fā)展地區(qū)經濟，限制或遷移高污染、高能耗企業(yè)。②加強地區(qū)生態(tài)服務價值研究，加強補償資金管理。

4 結論與討論

本文利用改進熵值法計算浙南各縣市的建設用地利用經濟效益、社會效益、環(huán)境效益和綜合效益，并結合空間自相關分析對建設用地進行利用分區(qū)，最終得到以下結論：

（1）浙南城市群各縣市建設用地效益區(qū)域差距顯著，效益高值區(qū)主要分布于東南部溫州各縣市，低值區(qū)主要分布于北部臺州市各縣市。

（2）通過對建設用地利用綜合效益進行空間自相關分析，利用Moran指數(shù)反映浙南城市群建設用地利用效益空間相關性特征，發(fā)現(xiàn)各縣市之間的建設用地利用綜合效益存在明顯的空間溢出作用。同時，浙南城市群建設用地效益呈現(xiàn)區(qū)域集聚性，南部平陽縣為High-High集聚性輻射熱點，北部縉云縣和臨?？h以及海島洞頭縣為Low-Low集聚型輻射冷點。

（3）結合綜合效益與空間自相關類型，對浙南城市群建設用地進行分區(qū)管制，分為調整優(yōu)化區(qū)、協(xié)調發(fā)展區(qū)、重點發(fā)展區(qū)、適度發(fā)展區(qū)和生態(tài)保護區(qū)，以進一步對地區(qū)建設用地利用進行精細化和差異化管理。

References)

[1] 劉新衛(wèi),張定祥,陳百明. 快速城鎮(zhèn)化過程中的中國城鎮(zhèn)土地利用特征[J]. 地理學報,2008,63(3):301-310.

Liu X W, Zhang D X, Chen B M. Characteristics of China's townlevel land use in rapid urbanization stage[J]. Acta Geographica Sinica, 2008,63(3):301-310.

[2] 周德,徐建春,王莉. 環(huán)杭州灣城市群土地利用的空間沖突與復雜性[J]. 地理研究,2015,34(9):1630-1642.

Zhou D, Xu J C, Wang L. Land use spatial conflicts and complexity: A case study of the urban agglomeration around Hangzhou Bay, China[J]. Geographical Research, 2015,34(9):1630-1642.

[3] 羅啟仕. 我國城市建設用地水土污染治理現(xiàn)狀與問題分析[J]. 上海國土資源,2015,36(4):59-63.

Luo Q S. Status of remediation of contaminated urban development land in China and analysis of existing problems[J]. Shanghai Land & Resources, 2015,36(4):59-63.

[4] 鄭沃林,鄭榮寶,唐曉蓮,等. 村鎮(zhèn)建設用地再開發(fā)后評估指標體系探討[J]. 中國土地科學,2016,30(4):70-78.

Zheng W L, Zheng R B, Tang X L, et al. Evaluation index system of village and town construction land redevelopment[J]. China Land Sciences, 2016,30(4):70-78.

[5] 李志強,趙凱,夏顯力. 西安市建設用地綜合效益評價[J]. 國土資源科技管理,2010,27(4):102-107.

Li Z Q, Zhao K, Xia X L. Comprehensive evaluation of construction land in Xi’an[J]. Scientific and Technological Management of Land and Resources, 2010,27(4):102-107.

[6] 錢競,羅罡輝,彭云飛. 基于因子分析的深圳市建設用地經濟效益評價[J]. 特區(qū)經濟,2013,(3):32-34.

Qian J, Luo G H, Peng Y F. Evaluation on the economic benefit of construction land in Shenzhen city based on factor analysis[J]. Special Zone Economy, 2013,(3):32-34.

[7] 吳斌,郭杰,殷爽,等. 江蘇省建設用地利用效益區(qū)域差異及分區(qū)管制[J]. 中國土地科學,2013,27(12):25-31.

Wu B, Guo J, Yin S, et al. Regional differences of construction land use efficiency and zoning regulation in Jiangsu province[J]. China Land Sciences, 2013,27(12):25-31.

[8] 劉桂菊,吳明發(fā). 經濟發(fā)達地區(qū)城市建設用地利用效益評價——以廣東省為例[J]. 江西農業(yè)學報,2015,27(11):99-102.

Liu G J, Wu M F. Benefit evaluation of urban construction land use in developed area: A case study of Guangzhou province[J]. Acta Agriculturae Jiangxi, 2015,27(11):99-102.

[9] 曹銀貴,周偉,喬陸印,等. 中國東部地區(qū)城市建設用地變化與利用效益分析[J]. 地理科學進展,2012,31(7):869-877.

Cao Y G, Zhou W, Qian L Y, et al. Research on urban construction land changes and using benefit in eastern China[J]. Progress in Geography, 2012,31(7):869-877.

[10] 強真,杜艦,吳尚昆. 我國城市建設用地利用效益評價[J]. 中國人口?資源與環(huán)境,2007,17(1):92-95.

Qiang Z, Du J, Wu S K. Appraisal on integrated benefit of urban construction land use in China[J]. China Population, Resources and Environment, 2007,17(1):92-95.

[11] 朱珠,張琳,葉曉雯,等. 基于TOPSIS方法的土地利用綜合效益評價[J]. 經濟地理,2012,32(10):139-144.

Zhu Z, Zhang L, Ye X W, et al. Evaluation of comprehensive land use efficiency based on TOPSIS[J]. Economic Geography, 2012,32(10):139-144.

[12] 王筱明,鄭新奇. 基于效益分析的濟南市城市合理用地規(guī)模研究[J]. 中國人口?資源與環(huán)境, 2010,20(6):160-165.

Wang X M, Zheng X Q. Study on rational urban land scale in Jinan based on benefits analysis[J]. China Population, Resource and Environment, 2010,20(6):160-165.

[13] 吳振華,唐芹,王亞蓓. 江浙滬地區(qū)城市建設用地生態(tài)效率評價——基于三階段DEA與Bootstrap-DEA方法[J]. 生態(tài)經濟,2016,32(4):105-110.

Wu Z H, Tang Q, Wang Y B. Evaluation of ecological efficiency of urban construction land in Jiangsu-Zhejiang-Shanghai Area: Based on three-stage DEA and Bootstrap-DEA methods[J]. Ecological Economy, 2016,32(4):105-110.

[14] 任飛陽,韋仕川. 區(qū)域建設用地節(jié)約集約利用評價研究——以三亞市為例[J]. 上海國土資源,2016,37(3):35-39.

Ren F Y, Wei S C. Evaluating regional economic and intensive land use: The case of Sanya city[J]. Shanghai Land & Resources, 2016,37(3):35-39.

[15] 趙勝男,任志遠,張晗,等. 2000-2013年西安市建設用地擴張?zhí)卣骷靶б娣治鯷J]. 陜西師范大學學報(自科版),2015,43(4):89-95.

Zhao S N, Ren Z Y, Zhang H, et al. Analysis of characteristics and benefits of urban sprawl in Xi’an city from 2000 to 2013[J]. Journal of Shanxi Normal University (Natural Science Edition), 2015,43(4):89-95.

[16] 王志寶,曹廣忠. 中國各省市區(qū)建設用地空間差異分析:1996-2004年[J]. 中國土地科學,2010,24(6):20-25.

Wang Z B, Cao G Z. Analysis on the spatial differences of construction land at provincial levels in China: From 1996 to 2004[J]. China Land Sciences, 2010,24(6):20-25.

[17] 柴志陽,范琰,董淼. 基于土地利用的縣域主體功能區(qū)劃分研究——以湖南省北湖區(qū)為例[J]. 上海國土資源,2013,34(4):62-66.

Chai Z Y, Fan Y, Dong M. A study of the division of major functional zones in counties based on land use: Beihu District, Hunan Province[J]. Shanghai Land & Resources, 2013,34(4):62-66.

[18] Deng X, Xu Y, Han L, et al. Assessment of river health based on an improved entropy-based fuzzy matter-element model in the Taihu Plain, China[J]. Ecological Indicators, 2015,57:85-95.

[19] 楊麗霞,夏浩,苑韶峰. 基于耦合協(xié)調度的土地利用經濟效益空間差異分析--以浙江省為例[J]. 中國土地科學,2015,29(11):83-88.

Yang L X, Xia H, Yuan S F. Spatial analysis on economic benefit of land utilization based on coupling coordination degree: A case study of Zhejiang Province[J]. China Land Sciences, 2015,29(11):83-88.

The regulated zoning of construction land in the urban agglomeration of south Zhejiang based on spatial autocorrelation analysis

XIA Hao1, YUAN Shao-Feng1, YANG Li-Xia2, WANG Jing-Yao1
(1. School of Public Administration, Zhejiang Gongshang University, Zhejiang Hangzhou 310018, China; 2. School of Urban and Rural Planning Management, Zhejiang University of Finance & Economics, Zhejiang Hangzhou 310018,, China)

Combined with a comprehensive evaluation of the benefits of construction land and spatial autocorrelation analysis, the purpose of this study is to provide a basis for regional construction land regulated zoning and precise management. The results show the following. (1) There is a significant difference between the counties that benefit and there are spatial spillover effects. (2) The benefits of construction land use show regional agglomeration. (3) The construction land of the urban agglomeration of southern Zhejiang can be divided into adjusting and optimizing zones: a coordinated development zone, key development zone, moderate development zone, and ecological reserve zone. The study concludes that regulated zoning of construction land should be based on the characteristics of land and regional conditions, and strategies need to be differentiated.

construction land; benefit evaluation; regulated zoning; spatial autocorrelation analysis; the urban agglomeration of south Zhejiang

F301.24

A

2095-1329(2017)01-0009-05

10.3969/j.issn.2095-1329.2017.01.003

2016-11-28

修回日期: 2017-01-30

夏浩(1994-),男,主要研究方向為土地利用效益評價.

電子郵箱: ryan-1h@163.com

聯(lián)系電話: 0571-28008341基金項目: 國家自然科學基金項目(41171151, 41371188);浙江工商大學教育部創(chuàng)新訓練計劃(201510353030)