深圳市基本生態(tài)控制線劃定對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響

吳健生,黃　喬,曹祺文

1 北京大學(xué)城市規(guī)劃與設(shè)計(jì)學(xué)院,城市人居環(huán)境科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室, 深圳　518055 2 北京大學(xué)城市與環(huán)境學(xué)院,地表過(guò)程與模擬教育部重點(diǎn)實(shí)驗(yàn)室, 北京　100871

生態(tài)系統(tǒng)服務(wù)(Ecosystem services)是人類從生態(tài)系統(tǒng)中所獲得的各種惠益[1],其價(jià)值評(píng)估是生態(tài)環(huán)境保護(hù)、生態(tài)功能區(qū)劃、環(huán)境經(jīng)濟(jì)核算和生態(tài)補(bǔ)償決策的重要依據(jù)和基礎(chǔ)[2]。1997年Costanza等人發(fā)表在Nature上關(guān)于全球生態(tài)系統(tǒng)服務(wù)價(jià)值核算的文章[3]引發(fā)了研究者對(duì)生態(tài)系統(tǒng)服務(wù)的關(guān)注,生態(tài)系統(tǒng)服務(wù)迅速成為生態(tài)學(xué)、地理學(xué)和環(huán)境科學(xué)的研究熱點(diǎn)和前沿[4]。國(guó)內(nèi)外關(guān)于生態(tài)系統(tǒng)服務(wù)價(jià)值(Ecosystem services value, ESV)評(píng)估的方法的研究很多,但尚未形成一套統(tǒng)一的評(píng)估體系[5- 8]。其中,謝高地等[9]在Costanza的研究成果上建立的中國(guó)陸地生態(tài)系統(tǒng)單位面積服務(wù)價(jià)值表的研究成果憑借直觀易用、數(shù)據(jù)量需求少等特點(diǎn)為廣大學(xué)者引用。土地利用變化蘊(yùn)含大量人類活動(dòng)的信息,可以通過(guò)改變生態(tài)系統(tǒng)類型、空間格局以及生態(tài)過(guò)程直接影響生態(tài)系統(tǒng)服務(wù)[10],其格局的變化對(duì)生態(tài)系統(tǒng)服務(wù)有著顯著的影響[11]。隨著土地利用土地覆被變化(LUCC)成為全球變化的熱點(diǎn),其對(duì)生態(tài)系統(tǒng)服務(wù)的影響也受到越來(lái)越多的關(guān)注[10- 16]。

《基本生態(tài)控制線》(后簡(jiǎn)稱“控制線”)是深圳市于2005年頒布的一項(xiàng)生態(tài)政策。政策頒布后,深圳市的城市建設(shè)與發(fā)展多了一項(xiàng)強(qiáng)制性的約束條件,土地利用格局也因此發(fā)生了很大的變化,與之密切相關(guān)的生態(tài)系統(tǒng)服務(wù)價(jià)值可能也發(fā)生了變化。

倍差法最早由Meyer[17]于1995年在《Natural and Quasi-Experiments in Economics》中提出,由于操作簡(jiǎn)單且邏輯清晰,該方法已被廣泛應(yīng)用在公共政策分析和工程評(píng)估中,其優(yōu)點(diǎn)在于直接剔除公共政策實(shí)施前后的相同影響因素的影響,留下公共政策本身的影響,在經(jīng)濟(jì)學(xué)領(lǐng)域應(yīng)用廣泛[18- 20],近年來(lái)也逐漸被引用到別的領(lǐng)域,如環(huán)境科學(xué),研究高速路收費(fèi)與空氣污染之間的關(guān)系[21]。

為探究控制線對(duì)深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值的影響,本研究基于深圳市2002、2005和2010年3期土地利用數(shù)據(jù)和控制線數(shù)據(jù),分析了政策頒布前后深圳市土地利用變化情況,進(jìn)而估算深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值,探討生態(tài)政策與生態(tài)系統(tǒng)服務(wù)價(jià)值之間的關(guān)系,以期為深圳市土地資源的合理優(yōu)化配置,生態(tài)環(huán)境保護(hù),生態(tài)政策的制定、調(diào)整及優(yōu)化提供決策依據(jù)。

1　研究區(qū)與研究方法

1.1　研究區(qū)概況

深圳市位于中國(guó)南部海濱,珠江三角洲東岸,113°46′—114°37′E,22°27′—22°52′N。改革開放以來(lái),深圳市從一個(gè)以農(nóng)業(yè)為主的邊陲小鎮(zhèn)轉(zhuǎn)變成高度城市化的大都市,景觀格局經(jīng)歷了巨大的變化。城市建設(shè)用地不斷擴(kuò)張,生態(tài)用地不斷被割裂和蠶食,生態(tài)環(huán)境出現(xiàn)了惡化的趨勢(shì)[22],生態(tài)系統(tǒng)服務(wù)不可避免受到重大影響。

為此,深圳市于2005年頒布了《深圳市基本生態(tài)控制線管理規(guī)定》,成為全國(guó)第一個(gè)將生態(tài)保護(hù)納入法規(guī)的城市。該政策要求將一級(jí)水源保護(hù)區(qū)、風(fēng)景名勝區(qū)、自然保護(hù)區(qū)、集中成片的基本農(nóng)田保護(hù)區(qū)、森林及郊野公園；坡度大于25%的山地、林地以及特區(qū)內(nèi)海拔超過(guò)50 m、特區(qū)外海拔超過(guò)80 m的高地等納入了控制線范圍,最終約占全市面積一半的土地被劃入了控制線內(nèi)[23]?？刂凭€內(nèi)嚴(yán)格限制各種建設(shè)活動(dòng),且已建合法建筑不得擅自改建和擴(kuò)建。

1.2　數(shù)據(jù)

本文基于深圳市2002、2005和2010年3期土地變更調(diào)查數(shù)據(jù),參照國(guó)土資源部2002年實(shí)行的“土地利用分類系統(tǒng)”進(jìn)行土地利用重分類,并依據(jù)研究需要將其合并為耕地、園地、林地、草地、建設(shè)用地、水域和其他用地七類。控制線數(shù)據(jù)由深圳市城市規(guī)劃委員會(huì)(http：//www.szplan.gov.cn/szupb/)于2005年公布,并根據(jù)公布面積對(duì)其進(jìn)行修正的數(shù)據(jù)。此外,還用到了深圳市2002、2005和2010年統(tǒng)計(jì)年鑒上糧食產(chǎn)量相關(guān)社會(huì)經(jīng)濟(jì)數(shù)據(jù)等。

1.3　研究方法

本文擬通過(guò)土地利用動(dòng)態(tài)動(dòng)度指標(biāo)定量描述深圳市控制線政策頒布前后土地利用的變化特征,在此基礎(chǔ)上利用價(jià)值系數(shù)法對(duì)深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行估算,最后利用倍差法定量分析控制線政策對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響。

1.3.1土地利用動(dòng)態(tài)度

土地利用動(dòng)態(tài)度是用來(lái)描述某一地區(qū)一定時(shí)間段內(nèi)各類土地利用類型的數(shù)量變化及變化速率情況[24]。計(jì)算公式如下：

(1)

式中,K表示研究時(shí)段內(nèi)某一土地利用類型的動(dòng)態(tài)度；Ua,Ub分別表示研究初期及研究末期某一土地利用類型的面積；T表示研究時(shí)段長(zhǎng)度,當(dāng)T的單位設(shè)定為年時(shí),K的值就是研究區(qū)某地類的年變化率[25]。后文所用的即為土地利用類型的年變化率。

1.3.2生態(tài)系統(tǒng)服務(wù)價(jià)值估算

本文通過(guò)謝高地等[9]提出的我國(guó)陸地生態(tài)系統(tǒng)服務(wù)的價(jià)值系數(shù)來(lái)計(jì)算深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值,并根據(jù)深圳市的實(shí)際情況,對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子表進(jìn)行部分修正：耕地與農(nóng)田對(duì)應(yīng),園地取森林和草地的平均值,林地與森林對(duì)應(yīng),水域與水體對(duì)應(yīng),其他用地與裸地對(duì)應(yīng),建設(shè)用地賦值為0,得到深圳市生態(tài)系統(tǒng)單位面積生態(tài)服務(wù)價(jià)值當(dāng)量因子。將深圳市2002、2005和2010年3年的糧食平均產(chǎn)量(5915 kg/hm2)以及稻谷平均收購(gòu)價(jià)(2.28 元/kg)作為基準(zhǔn)單產(chǎn)和價(jià)格確定深圳市平均糧食單產(chǎn)市場(chǎng)價(jià)值,進(jìn)而以其價(jià)值的1/7[9]確定1個(gè)生態(tài)系統(tǒng)服務(wù)價(jià)值當(dāng)量因子的經(jīng)濟(jì)價(jià)值量為13486.2 元/hm2,利用公式(2)計(jì)算出深圳市不同土地利用類型單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(表1)。

ESV=∑(Ak×VCk)

(2)

式中,ESV為生態(tài)系統(tǒng)服務(wù)價(jià)值；Ak為土地利用類型k的面積；VCk為生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)。

表1　深圳市不同土地利用類型單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)/(元/hm2)

ESV: Ecosystem serices value

1.3.3倍差法

分析控制線劃定是否會(huì)對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生影響,重點(diǎn)在于比較控制線劃定前后兩個(gè)時(shí)間段的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化。但是直接進(jìn)行前后比較不可能得到準(zhǔn)確結(jié)論,因?yàn)榭刂凭€劃定前后生態(tài)系統(tǒng)服務(wù)價(jià)值的差異可能不只是政策作用的結(jié)果,還包含了大量其他因素,如自然環(huán)境的變化、社會(huì)經(jīng)濟(jì)因素的影響等。本文擬采用倍差法(Difference-in-Differences approach)來(lái)解決這個(gè)問(wèn)題。

具體而言,構(gòu)造劃定控制線的處理組和沒(méi)有劃定控制線的對(duì)照組,通過(guò)控制其他因素,比較控制線劃定前后的處理組和對(duì)照組之間的差異,從而檢驗(yàn)控制線劃定對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響。

將研究區(qū)劃分為500 m×500 m的網(wǎng)格,根據(jù)每個(gè)網(wǎng)格中各土地利用類型的面積及單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)計(jì)算得到每個(gè)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值。將所有的網(wǎng)格劃分為4組：控制線劃定前的處理組和對(duì)照組、控制線劃定后的處理組對(duì)照組。本研究將控制線以內(nèi)的網(wǎng)格作為處理組,控制線以外的網(wǎng)格作為對(duì)照組,用2005年表示控制線劃定前,2010年表示控制線劃定后。設(shè)置兩個(gè)虛擬變量來(lái)度量樣本的劃分,一個(gè)是du,處理組取值為1,對(duì)照組為0；另一個(gè)是dt,控制線劃定前取值為0,劃定后為1。具體回歸方程設(shè)定如下：

fit=β0+β1duit+β2dtit+β3duit×dtit+εit

(3)

式中,i和t分別表示網(wǎng)格編號(hào)和時(shí)間,f和ε分別為生態(tài)系統(tǒng)服務(wù)價(jià)值和擾動(dòng)項(xiàng)。顯然,系數(shù)β3度量了控制線劃定對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響,即控制線的政策效應(yīng),解釋如下：

在對(duì)照組,即duit=0,由式(3)可知,控制線劃定前后生態(tài)系統(tǒng)服務(wù)價(jià)值分別記為：

可見,控制線劃定前后,對(duì)照組的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化為(β0+β2+εit)―(β0+δit),即Δf1=β2。

相應(yīng)的,在處理組,即duit=1,由式(3)可知,控制線劃定前后生態(tài)系統(tǒng)服務(wù)價(jià)值分別記為：

可見,控制線劃定前后,處理組的生態(tài)系統(tǒng)服務(wù)價(jià)值的變化為(β0+β1+β2+εit)―(β0+β1+δit),即Δf2=β2+β3。

因此,控制線劃定對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的凈影響為Δf2-Δf1=β3,即交叉項(xiàng)duit×dtit的系數(shù)β3。

由于控制線劃定需要政府批準(zhǔn),可視為外生變量,可以假定與回歸方程的誤差項(xiàng)不相關(guān),即E(εit|duit×dtit)=0。

在實(shí)證分析過(guò)程中,在上述回歸方程中加入一個(gè)新的變量,建設(shè)用地面積占比xit(每個(gè)網(wǎng)格中建設(shè)用地面積與該網(wǎng)格面積的比值),得到以下回歸方程：

fit=α+β1duit+β2dtit+β3duit×dtit+β4xit+εit

(4)

本文關(guān)注參數(shù)β3,當(dāng)β3為正時(shí),表示控制線劃定提高了研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值,反之,控制線對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值無(wú)作用。

2　結(jié)果與分析

2.1　土地利用變化分析

整體來(lái)看土地利用類型變化表現(xiàn)為：建設(shè)用地面積占比最大,且不斷增大；耕地、園地、草地、水域面積不斷減小,林地和其他用地面積先減小后增大(表2)。具體從年變化率來(lái)看,控制線劃定前后,面積不斷減小的耕地、園地、草地、水域,后期的減小速率均小于前期減小的速率；前期建設(shè)用地是唯一年變化率為正值的地類,其面積前期增大的速率遠(yuǎn)大于后期,建設(shè)用地?cái)U(kuò)張遇到明顯的阻力；其他用地面積減小的速率最大,林地面積變化的趨勢(shì)由減小轉(zhuǎn)變成增大,林地作為主要的生態(tài)用地,受到了生態(tài)政策明顯的保護(hù)和恢復(fù)作用。

表3給出了控制線內(nèi)/外各土地利用類型面積比例以及動(dòng)態(tài)度變化的情況?？梢钥吹?園地和林地的絕大部分面積都劃分到控制線以內(nèi),耕地、草地、水域面積次之,建設(shè)用地面積占比最小。隨著時(shí)間的推移,控制線內(nèi)耕地、園地、水域面積比例不斷提高,林地面積比例變化不大,建設(shè)用地面積比例增加速率變小,草地和其他用地面積比例先增大后減小?？刂凭€外耕地、園地水域面積比例不斷降低,林地面積比例變化不大,建設(shè)用地面積比例略有減小,草地和其他用地面積比例先減小后增大。從線內(nèi)/外各土地利用類型的動(dòng)態(tài)度來(lái)看,控制線劃定前,除建設(shè)用地外各地類面積變化情況呈現(xiàn)線內(nèi)小于線外的趨勢(shì)；控制線劃定后普遍存在控制線內(nèi)變化小于線外變化的現(xiàn)象,且各地類面積變化幅度明顯小于控制線劃定前。

表2　深圳市土地利用類型面積比例及動(dòng)態(tài)度

表3　控制線內(nèi)/外土地利用類型面積比例及動(dòng)態(tài)度

2.2　生態(tài)系統(tǒng)服務(wù)價(jià)值估算

整體來(lái)看,2002—2010年間,深圳市生態(tài)系統(tǒng)服務(wù)總價(jià)值不斷降低,從2002年的376.41億元下降到2005年的333.56億元,降幅11.38%,再降到2010年的314.81億元,降幅5.62%(表4),生態(tài)系統(tǒng)服務(wù)總價(jià)值降低的速率減小。從各土地利用類型來(lái)看,耕地、園地、草地和水域的價(jià)值均不斷降低,而林地和其他用地的價(jià)值出現(xiàn)了先降低后升高的變化,控制線劃定后生態(tài)系統(tǒng)服務(wù)價(jià)值回升,但林地和其他用地增加的價(jià)值并不能彌補(bǔ)其他地類價(jià)值的損失,因此,研究時(shí)段內(nèi)深圳市生態(tài)系統(tǒng)服務(wù)總價(jià)值是不斷減小的。從生態(tài)系統(tǒng)服務(wù)總價(jià)值的構(gòu)成來(lái)看,林地和水域的價(jià)值是總價(jià)值的主要組成部分,兩種地類價(jià)值之和占比超過(guò)80%。園地的價(jià)值占比較小,約為14%—16%,草地和其他用地的服務(wù)價(jià)值占比最小。盡管各土地利用類型對(duì)總生態(tài)系統(tǒng)服務(wù)價(jià)值的重要程度分異明顯,2002—2010年各土地利用類型服務(wù)價(jià)值的比例變化并不顯著。從控制線內(nèi)/外生態(tài)系統(tǒng)服務(wù)價(jià)值來(lái)看,雖然線內(nèi)/外價(jià)值均在不斷降低,但線內(nèi)的價(jià)值均遠(yuǎn)高于線外,并且兩者之間的差異在不斷增大,各土地利用類型變化表現(xiàn)為耕地、水域和其他用地價(jià)值線內(nèi)低于線外,園地和林地價(jià)值線內(nèi)高于線外,尤其線內(nèi)林地價(jià)值遠(yuǎn)遠(yuǎn)高于線外林地價(jià)值。

表4　深圳市2002—2010年生態(tài)系統(tǒng)服務(wù)價(jià)值/億元

從表5可以看到四大類服務(wù)價(jià)值、各項(xiàng)服務(wù)產(chǎn)生的價(jià)值及其比例?？傮w來(lái)看,從2002年到2010年,深圳市四大類服務(wù)價(jià)值均減小,但控制線劃定后價(jià)值減小的速率有所減緩。各項(xiàng)服務(wù)價(jià)值結(jié)構(gòu)保持不變,不過(guò)各項(xiàng)服務(wù)價(jià)值之間的差距有所減小,生態(tài)系統(tǒng)服務(wù)價(jià)值差異極大值由2002年的90.25億元減小至2010年的70.60億元。各項(xiàng)生態(tài)系統(tǒng)服務(wù)按照服務(wù)價(jià)值大小排序依次為：水源涵養(yǎng)>廢物處理>土壤形成與保護(hù)>生物多樣性保護(hù)>氣體調(diào)節(jié)>氣候調(diào)節(jié)>原材料>娛樂(lè)文化>食物生產(chǎn)。雖然各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值排序沒(méi)有發(fā)生變化,但各項(xiàng)服務(wù)價(jià)值的變化趨勢(shì)有所不同,其中氣體調(diào)節(jié)、氣候調(diào)節(jié)、土壤形成與保護(hù)、生物多樣性保護(hù)和原材料幾項(xiàng)服務(wù)價(jià)值增大,而水源涵養(yǎng)、廢物處理、食物生產(chǎn)和娛樂(lè)文化價(jià)值減小。

表5　深圳市2002—2010年各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值

“總計(jì)”是各大類服務(wù)價(jià)值的總和；“排序”是各項(xiàng)服務(wù)產(chǎn)生的價(jià)值對(duì)總價(jià)值的貢獻(xiàn)排序；“趨勢(shì)”表示的是研究時(shí)段內(nèi)各項(xiàng)服務(wù)價(jià)值的變化趨勢(shì),“↑”表明該項(xiàng)服務(wù)價(jià)值增大,“↓”表明該項(xiàng)服務(wù)價(jià)值減小[26]

將研究區(qū)劃分為500 m×500 m的網(wǎng)格,根據(jù)每個(gè)網(wǎng)格中各土地利用類型的面積計(jì)算每個(gè)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值,采用自然斷點(diǎn)法將其劃分為低、較低、中、較高和高5個(gè)等級(jí),最終得到2002—2010年深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值等級(jí)空間分布(圖1)。從空間上來(lái)看,深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值整體呈現(xiàn)下降的趨勢(shì),價(jià)值低的范圍不斷擴(kuò)大,蠶食價(jià)值高的區(qū)域,尤其以西部沿海、寶安北部、光明新區(qū)、龍華北部以及龍崗南部地區(qū)價(jià)值降低最為顯著,大鵬新區(qū)存在部分零星價(jià)值升高的區(qū)域?？刂凭€劃定之前,生態(tài)系統(tǒng)服務(wù)價(jià)值的高值區(qū)快速流失,而控制線劃定之后,情況相對(duì)好轉(zhuǎn),除了西部沿海地區(qū)水域面積縮減造成生態(tài)系統(tǒng)服務(wù)價(jià)值明顯降低之外其他區(qū)域變化相對(duì)不顯著。

圖1　2002—2010年深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值Fig.1　ESV in Shenzhen from 2002 to 2010ESV: 生態(tài)系統(tǒng)服務(wù)價(jià)值Ecosystem serices value

控制線內(nèi)/外生態(tài)系統(tǒng)服務(wù)價(jià)值的空間變化分別如圖2、3所示。對(duì)比兩張圖可以發(fā)現(xiàn),控制線以內(nèi)生態(tài)系統(tǒng)服務(wù)價(jià)值以中、較高為主,高值零星分布,整體價(jià)值較高；控制線以外生態(tài)系統(tǒng)服務(wù)價(jià)值以低、較低為主,高值主要分布在西部沿海地區(qū)。主要原因是線內(nèi)以林地、園地等生態(tài)用地類型為主,單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值較高；線外主要以建設(shè)用地為主,單位生態(tài)系統(tǒng)服務(wù)價(jià)值很低。由于水域單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)是各土地類型的極大值,因此水域分布的地區(qū)價(jià)值最高。

從時(shí)間變化上看,生態(tài)系統(tǒng)服務(wù)價(jià)值在控制線以內(nèi)的變化不明顯,而在控制線以外,表現(xiàn)出降低的趨勢(shì),尤其以2002—2005年降低最為顯著,價(jià)值降低的區(qū)域主要分布在寶安北部、光明新區(qū)、龍華北部以及龍崗南部地區(qū)；2005—2010年西南部地區(qū)水域面積收縮,造成生態(tài)系統(tǒng)服務(wù)價(jià)值的降低?？梢?控制線以內(nèi)限制開發(fā)的區(qū)域生態(tài)系統(tǒng)服務(wù)的價(jià)值得以維持和保護(hù),而控制線以外非限制性區(qū)域,由于城市的發(fā)展,價(jià)值繼續(xù)降低。

圖2　控制線內(nèi)生態(tài)系統(tǒng)服務(wù)價(jià)值Fig.2　ESV within the control line

圖3　控制線外生態(tài)系統(tǒng)服務(wù)價(jià)值Fig.3　ESV outside the control line

2.3　控制線政策對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的效應(yīng)分析

通過(guò)倍差法,構(gòu)建了生態(tài)系統(tǒng)服務(wù)價(jià)值與控制線政策之間的回歸模型。模型樣本共計(jì)22292個(gè)網(wǎng)格,將每個(gè)網(wǎng)格的生態(tài)系統(tǒng)服務(wù)價(jià)值作為因變量,du、dt、du×dt、建設(shè)用地占比作為自變量,利用SPSS軟件進(jìn)行多元線性回歸,得到模型回歸結(jié)果如表6。

從模型回歸結(jié)果來(lái)看,最終多元回歸模型回歸方程為ESV=0.26du-0.33dt+0.016du×dt-0.769x。模型通過(guò)F檢驗(yàn),模型整體有效；擬合優(yōu)度為0.617,自變量對(duì)因變量的解釋程度較高,擬合較好；自變量的回歸系數(shù)均通t-檢驗(yàn)。建設(shè)用地占比系數(shù)值為負(fù),說(shuō)明建設(shè)用地面積的比例對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值存在負(fù)效應(yīng),即建設(shè)用地面積占比越大,生態(tài)系統(tǒng)服務(wù)價(jià)值越小,與事實(shí)吻合,進(jìn)一步說(shuō)明模型的有效性。本文重點(diǎn)關(guān)注的du×dt的系數(shù)β3的值為0.016,值為正,說(shuō)明控制線的劃定能夠提高且大約提高1.6%的生態(tài)系統(tǒng)服務(wù)價(jià)值。由此可見,控制線的劃定確實(shí)對(duì)區(qū)域生態(tài)系統(tǒng)服務(wù)具有促進(jìn)作用,其政策效應(yīng)得以定量表現(xiàn),而前文關(guān)于兩者關(guān)系的描述性分析具有可靠性。

表6　倍差法結(jié)果

***表示P<0.01

3　結(jié)論與討論

3.1　結(jié)論

本研究基于2002年、2005年和2010年深圳市控制線劃定前后的土地變更調(diào)查數(shù)據(jù),分析了控制線政策頒布前后深圳市土地利用的變化情況,運(yùn)用基于單位面積價(jià)值當(dāng)量因子法對(duì)深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行評(píng)估并分析,最后利用倍差法定量研究了控制線政策對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響。結(jié)果表明：(1)深圳市林地和水域的價(jià)值最大,且價(jià)值之和超過(guò)總價(jià)值的80%?？刂凭€劃定之前,隨著城市化進(jìn)程的推進(jìn),生態(tài)用地被建設(shè)用地蠶食,造成生態(tài)系統(tǒng)服務(wù)總價(jià)值減少42.9億元；而控制線劃定之后,生態(tài)用地受到政策的保護(hù),面積有所增加,但并不足以彌補(bǔ)建設(shè)用地?cái)U(kuò)張帶來(lái)的價(jià)值損失,所以總價(jià)值繼續(xù)減少18.74億元,減少的速率明顯下降。生態(tài)系統(tǒng)服務(wù)價(jià)值線內(nèi)遠(yuǎn)高于線外,價(jià)值降低速率線內(nèi)遠(yuǎn)低于線外。各項(xiàng)服務(wù)按照價(jià)值大小排序依次為：水源涵養(yǎng)>廢物處理>土壤形成與保護(hù)>生物多樣性保護(hù)>氣體調(diào)節(jié)>氣候調(diào)節(jié)>原材料>娛樂(lè)文化>食物生產(chǎn)。(2)空間上,低價(jià)值區(qū)域面積不斷擴(kuò)大,高價(jià)值區(qū)域被蠶食,尤其以西部沿海、寶安北部、光明新區(qū)、龍華北部以及龍崗南部地區(qū)價(jià)值降低最明顯?？刂凭€劃定之前生態(tài)系統(tǒng)服務(wù)價(jià)值快速流失,而控制線劃定之后,情況稍有好轉(zhuǎn)。控制線內(nèi)生態(tài)系統(tǒng)服務(wù)價(jià)值以較高、高為主,整體價(jià)值較高,并且隨時(shí)間的變化不明顯；控制線外生態(tài)系統(tǒng)服務(wù)價(jià)值以低、較低為主,價(jià)值較低,且呈現(xiàn)不斷降低的趨勢(shì),尤其控制線劃定之前價(jià)值降低顯著。(3)倍差法得到控制線政策對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響為正,控制線的劃定大約提高1.6%的生態(tài)系統(tǒng)服務(wù)價(jià)值。

3.2　倍差法適用性

本文所采用的倍差法是經(jīng)濟(jì)學(xué)、公共管理學(xué)中廣泛使用的政策效應(yīng)評(píng)估工具,目前在生態(tài)政策效應(yīng)研究中應(yīng)用較少。由于政策具有難量化的特點(diǎn),對(duì)政策效應(yīng)的一般定性描述往往不能排除政策之外其他因素的影響。倍差法操作簡(jiǎn)單且能較好解決政策難定量化、難剝離的問(wèn)題,得到越來(lái)越廣泛的運(yùn)用。

控制線政策是深圳市頒布的一項(xiàng)強(qiáng)制性空間生態(tài)政策,任何城市開發(fā)活動(dòng)都需要遵守控制線的規(guī)定。研究控制線政策對(duì)深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值的影響同樣屬于政策的效應(yīng)研究范疇,倍差法具有一定適用性；而文章使用倍差法的結(jié)果也表明控制線政策對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值確實(shí)是存在影響的,證實(shí)了倍差法在生態(tài)政策效應(yīng)評(píng)估中的適用性。

3.3　不足與展望

由于土地利用的分類與自然生態(tài)系統(tǒng)類型并非完全對(duì)應(yīng),所以各種土地利用類型單位面積服務(wù)價(jià)值的取值只是一種近似[26- 27]。使用不同的方法計(jì)算得到的深圳市生態(tài)系統(tǒng)服務(wù)總價(jià)值可能有較大的偏差。本研究采用基于單位面積價(jià)值當(dāng)量因子的方法對(duì)深圳市生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行評(píng)估,雖然已經(jīng)根據(jù)深圳市的實(shí)際情況對(duì)價(jià)值當(dāng)量因子表進(jìn)行部分修正,但發(fā)現(xiàn)單位面積水域的水源涵養(yǎng)服務(wù)價(jià)值當(dāng)量因子偏高,導(dǎo)致計(jì)算得到的水源涵養(yǎng)服務(wù)價(jià)值系數(shù)偏高,水域價(jià)值計(jì)算存在一定的誤差。不過(guò),本研究側(cè)重分析生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空變化情況,可以消除由生態(tài)系統(tǒng)服務(wù)價(jià)值單價(jià)取值不同帶來(lái)的影響,也能夠客觀反映研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的變化,因此,研究的結(jié)論依舊具有可靠性。

本研究在進(jìn)行倍差法分析時(shí),所考量的生態(tài)系統(tǒng)服務(wù)價(jià)值變化僅基于單一時(shí)間段,今后應(yīng)進(jìn)一步考慮做長(zhǎng)時(shí)間序列上的動(dòng)態(tài)變化,以更為準(zhǔn)確地衡量城市生態(tài)空間政策的時(shí)間效應(yīng),進(jìn)而輔助決策者及時(shí)調(diào)整生態(tài)管控政策,也為建設(shè)可持續(xù)發(fā)展的生態(tài)城市提供理論和實(shí)踐基礎(chǔ)。