周申蓓，黃媛媛，呂玲玲

(河海大學 商學院，江蘇 南京 211100）

小浪底水利樞紐工程位于黃河中游峽谷出口,控制流域面積69.4萬km2,集防洪、調(diào)沙、發(fā)電和灌溉等功能于一體,是黃河綜合治理的控制性工程。工程的興建和使用,一方面極大改善了黃河下游防洪形勢,另一方面顯著影響了黃河下游河岸帶的開發(fā)和保護狀態(tài)。已有研究評估了小浪底工程運行對黃河下游河道泥沙輸送、生態(tài)環(huán)境、浮游生物及魚類、農(nóng)業(yè)景觀、沿線濕地和河口保護等方面的影響。有研究顯示,水庫運行減少了黃河下游泥沙淤積,對生態(tài)環(huán)境產(chǎn)生了有利影響[1］。水庫運行總體對黃河下游水生生物資源的數(shù)量和種類產(chǎn)生了不利影響,對河流生態(tài)功能服務的正面影響體現(xiàn)在資源開發(fā)和水質(zhì)凈化等方面[2-4］。目前的研究從不同方面評估了小浪底工程及人類活動對黃河下游生態(tài)環(huán)境的影響,但評估的整體性尚有不足。為此,以黃河河岸帶為研究對象,將小浪底工程影響區(qū)域劃分為庫區(qū)上游100 km段、庫區(qū)下游焦作—菏澤段、聊城—濱州段以及入海河口段4段,探討小浪底工程建設和運營對上下游河岸帶生態(tài)系統(tǒng)服務價值的影響,以期為黃河流域生態(tài)系統(tǒng)空間格局優(yōu)化及可持續(xù)管理提供參考。

1 研究區(qū)概況

1.1 研究區(qū)地理位置

研究區(qū)為小浪底工程所在位置上游100 km處至黃河河口,共包括17個市。工程上游100 km以內(nèi)河岸帶位于山西省運城市,河南省三門峽市和洛陽市；工程下游河岸帶位于河南省濟源市、焦作市、鄭州市、開封市、新鄉(xiāng)市與濮陽市,山東省聊城市、泰安市、濟南市、德州市、濱州市、淄博市與菏澤市；入海口位于山東省東營市。

1.2 小浪底工程基本情況

小浪底水利樞紐工程坐落于河南省孟津縣與濟源市交界,位于三門峽水利樞紐下游130 km、洛陽市以北40 km的黃河干流,是黃河干流三門峽以下唯一能取得較大庫容的控制性工程。小浪底工程于1991年4月開工,2001年12月主體工程完工,2002年正式投入運營。小浪底水庫正常蓄水位275 m,死水位230 m,設計洪水位274 m,校核洪水位275 m,總庫容126.5億m3。

2 數(shù)據(jù)來源與研究方法

2.1 數(shù)據(jù)來源與處理

小浪底工程建設并使用至今共30余a,本文土地利用數(shù)據(jù)取自中國科學院資源環(huán)境科學與數(shù)據(jù)中心中國陸地生態(tài)系統(tǒng)宏觀結(jié)構(gòu)數(shù)據(jù)中1990年、1995年、2000年、2005年、2010年、2015年及2020年7期分辨率為1 km的遙感監(jiān)測數(shù)據(jù)。上述數(shù)據(jù)覆蓋了小浪底工程建設前、建設中和建成后運營期不同階段。根據(jù)《土地利用現(xiàn)狀分類》(GB/T21010—2017)及小浪底工程影響區(qū)域?qū)嶋H情況,將土地利用類型劃分為農(nóng)田、林地、草地、水域、濕地和建設用地6類。研究區(qū)的社會發(fā)展、經(jīng)濟等數(shù)據(jù)來源于《中國統(tǒng)計年鑒》《中國環(huán)境年鑒》《全國農(nóng)產(chǎn)品成本收益資料匯編》《中國環(huán)境統(tǒng)計年鑒》等。

2.2 研究方法

2.2.1 河岸帶范圍確定

河岸帶包括非永久被淹沒的河床與周圍新生的和殘余的洪泛平原,是受水陸相互作用的地區(qū),其范圍的確定受河流徑流量、土壤以及其他自然環(huán)境影響。蘇晨等[5］研究表明,小浪底工程運行20 a后,黃河下游影響帶寬度為5～25 km,以新鄉(xiāng)—開封一帶最寬,黃河入海口處最窄,分別約25、5 km?？紤]研究區(qū)空間跨度以及各省市土地利用情況等,河岸帶影響區(qū)域為河流兩岸向外擴展10 km,即研究區(qū)包括河流及其河岸帶,總寬度接近25 km。研究時期為1990—2020年,其中1990年為工程建設前的初始狀態(tài),1995年、2000年為工程建設期,2005年、2010年為工程運營前期,2015年、2020年為工程運營后期。期間研究范圍保持一致,其中工程上游研究范圍按庫區(qū)蓄水后的水面及其河岸帶面積界定,總面積為2561 km2。

2.2.2 河岸帶生態(tài)系統(tǒng)服務價值評估

基于研究區(qū)土地利用實際情況,將河岸帶土地利用類型分為建設用地和非建設用地,其中非建設用地包括農(nóng)田、林地、草地、水域、濕地。建設用地生態(tài)系統(tǒng)服務價值采用替代成本法和防治成本法進行評估,非建設用地生態(tài)系統(tǒng)服務價值采用價值當量進行評估。

(1)非建設用地生態(tài)系統(tǒng)服務價值評估。參照謝高地等[6-8］修訂的基準年2010年全國單位面積生態(tài)系統(tǒng)服務價值當量因子表,根據(jù)河岸帶生態(tài)系統(tǒng)特征及當?shù)亟?jīng)濟社會系統(tǒng)實際情況,考慮產(chǎn)量和經(jīng)濟因素,對非建設用地生態(tài)系統(tǒng)服務價值當量進行修正。

產(chǎn)量因素修正系數(shù):

式中:α為農(nóng)田生態(tài)系統(tǒng)服務當量因子修正系數(shù)；σt為t時期研究區(qū)單位面積糧食產(chǎn)量；σ為t時期全國單位面積糧食產(chǎn)量。

根據(jù)支付意愿,結(jié)合不同時期研究區(qū)城鎮(zhèn)化水平構(gòu)建社會發(fā)展因子調(diào)整模型進行經(jīng)濟因素修正。

式中:lt為研究區(qū)社會發(fā)展系數(shù)；l為全國社會發(fā)展系數(shù)；l0為極富社會發(fā)展階段的支付意愿,取值為1；h為t時期研究區(qū)城鎮(zhèn)化水平；H為t時期全國城鎮(zhèn)化水平；β為t時期研究區(qū)支付意愿系數(shù)；En為t時期研究區(qū)恩格爾系數(shù)。

考慮價格變動因素,以2010年為基準年,采用可比價格計算農(nóng)田生態(tài)系統(tǒng)單位面積食物生產(chǎn)服務價值,具體計算公式:

式中:Di為第i年研究區(qū)1個標準當量因子的生態(tài)系統(tǒng)服務價值,元/hm2；Sri、Swi、Sci、Sbi分別為第i年研究區(qū)稻谷、小麥、玉米和大豆的播種面積占4種作物播種總面積的百分比,%；Fi為第i年研究區(qū)糧食的單位面積可比現(xiàn)金純收益,可分為Fri、Fwi、Fci和Fbi,分別為第i年研究區(qū)稻谷、小麥、玉米和大豆的單位面積可比現(xiàn)金純收益,元/hm2；Pi為第i年研究區(qū)稻谷、小麥、玉米和大豆的當年現(xiàn)金純收益,元/kg；APIi為第i年稻谷、小麥、玉米和大豆的生產(chǎn)價格指數(shù)；r為某類農(nóng)產(chǎn)品可比價格指數(shù)；i10為2010年。

采用修正后的各省(區(qū))當年單位面積生態(tài)系統(tǒng)服務價值當量因子表,計算非建設用地生態(tài)系統(tǒng)服務價值。

式中:ESVki為i時期k類生態(tài)系統(tǒng)服務價值；ESVie為i時期研究區(qū)非建設用地生態(tài)系統(tǒng)服務價值；rki為i時期修正后的研究區(qū)k類生態(tài)系統(tǒng)服務標準當量因子；A為k類生態(tài)系統(tǒng)土地利用面積。

(2)建設用地生態(tài)系統(tǒng)服務價值評估。河岸帶建設用地指因人類活動而實際被占用的土地,其內(nèi)的少量水域或綠地因沒有形成完整的生態(tài)功能而忽略不計。建設用地生活或生產(chǎn)活動廢棄物及有害物質(zhì)的排放往往對生態(tài)系統(tǒng)服務價值產(chǎn)生較強的負面影響,影響主要體現(xiàn)在氣體調(diào)節(jié)、水文調(diào)節(jié)和土壤保持三方面。水文調(diào)節(jié)價值采用替代成本法評估,氣體調(diào)節(jié)和土壤保持價值采用防治成本法評估,計算公式為

式中:Pg為建設用地氣體調(diào)節(jié)的單位面積生態(tài)系統(tǒng)服務價值；Qg為廢氣排放量；Cg為廢氣排污費單價；Sc為建設用地總面積；Pw為建設用地水文調(diào)節(jié)的單位面積生態(tài)系統(tǒng)服務價值；Qw為廢水排放量；Cw為廢水排污費單價；Pl為建設用地土壤保持的單位面積生態(tài)系統(tǒng)服務價值；Ql為固體廢棄物的排放量；Cl為固體廢棄物的排污費單價；ESVif為i時期建設用地生態(tài)系統(tǒng)服務價值。

研究區(qū)生態(tài)系統(tǒng)服務總價值即非建設用地生態(tài)系統(tǒng)服務價值和建設用地生態(tài)系統(tǒng)服務價值之和。

3 結(jié)果分析

3.1 研究區(qū)土地利用面積變化

典型河岸帶土地利用變化情況見圖1。在小浪底工程建設期及運營期,各類土地利用面積稍有變化。

圖1 典型河岸帶土地利用類型

(1)工程建設期農(nóng)田、草地面積增大,但工程建設期及運營期總體呈下降趨勢。農(nóng)田面積下降反映了研究區(qū)工業(yè)化及城鎮(zhèn)化發(fā)展與農(nóng)田生態(tài)系統(tǒng)之間的土地流轉(zhuǎn)關(guān)系。河口段農(nóng)田面積顯著減小,可能是小浪底工程建設導致水沙條件變化所致[9-10］。草地面積總體減小,2010年以后隨著環(huán)保意識的提高及土地管控政策趨嚴,草地面積變化不大[11-12］。

(2)水域面積變化符合大型水利工程建設規(guī)律[13-14］。工程建設期水域面積呈減小趨勢,其中1990—1995年減小速度最快；工程竣工后水域面積呈增大趨勢,2000年以后水域面積變化不大。

(3)庫區(qū)下游段建設用地面積逐漸增大,小浪底工程建設推動了下游沿岸行政區(qū)的社會發(fā)展。小浪底工程建設后,洪澇風險進一步降低,城鎮(zhèn)居民點及工礦用地、交通運輸用地及水利設施用地面積增大,研究區(qū)建設用地面積呈穩(wěn)步增大趨勢,2000年至2020年建設用地面積總體增加了857 km2,面積占比增加了4.14%。

(4)林地面積保持穩(wěn)定。黃河流域下游林地面積占比較小,1990—2020年共增加了12 km2,處于較穩(wěn)定狀態(tài),說明林地保護效果較好。

(5)濕地面積呈下降趨勢。河口區(qū)人類活動強度不斷增大,使得天然濕地面積減小,濕地系統(tǒng)結(jié)構(gòu)失調(diào)[15］。

3.2 研究區(qū)單位面積生態(tài)系統(tǒng)服務價值

在工程運營前期,研究區(qū)單位面積生態(tài)系統(tǒng)服務價值呈增大趨勢,其中水域生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務價值增幅最明顯。工程運營后期,研究區(qū)單位面積生態(tài)系統(tǒng)服務價值減小,原因主要是非建設用地生態(tài)系統(tǒng)產(chǎn)生的正向影響增幅不大,建設用地面積增加對生態(tài)系統(tǒng)產(chǎn)生的負面影響陡增。典型年研究區(qū)單位面積生態(tài)系統(tǒng)服務價值見表1。

3.3 各河段生態(tài)系統(tǒng)服務總價值

庫區(qū)上游河岸帶建設期生態(tài)系統(tǒng)服務總價值呈上升趨勢,運營前期生態(tài)系統(tǒng)服務價值增長幅度較大,見表2。相對于2000年,2010年上游河岸帶生態(tài)系統(tǒng)服務價值增長率為303%。運營后期庫區(qū)河岸帶生態(tài)系統(tǒng)服務總價值下降,主要是小浪底工程上游土地開發(fā)強度增大。

表1 典型年研究區(qū)單位面積生態(tài)系統(tǒng)服務價值 千元/hm2

表2 庫區(qū)上游河岸帶生態(tài)系統(tǒng)服務價值 百萬元

焦作—菏澤段河岸帶生態(tài)系統(tǒng)服務價值增長幅度與河岸帶距離小浪底工程所在位置正相關(guān)。不同于一般工程影響區(qū)域主要在工程緊鄰下游范圍,小浪底調(diào)水調(diào)沙將黃河下游主河槽過流能力提高到4000 m3/s,對黃河下游生態(tài)系統(tǒng)服務價值提升具有顯著促進作用,見表3、表4。這種促進作用在工程運營前期表現(xiàn)尤為突出,與2000年相比,2010年庫區(qū)下游焦作—菏澤段河岸帶生態(tài)系統(tǒng)服務總價值增長率為208%,庫區(qū)下游聊城—濱州段增長率為350%。

表3 庫區(qū)下游（焦作—菏澤段）生態(tài)系統(tǒng)服務價值 百萬元

表4 庫區(qū)下游（聊城—濱州段）生態(tài)系統(tǒng)服務價值 百萬元

小浪底工程運營對河口段河岸帶產(chǎn)生了顯著影響,見表5。2000—2010年河口段生態(tài)系統(tǒng)服務價值快速增大,與2000年相比,2005年河口段生態(tài)系統(tǒng)服務價值增長率為138%；2015年河口段生態(tài)系統(tǒng)服務價值下降,小浪底工程后續(xù)調(diào)水調(diào)沙動力不足[16-17］導致此時期草地、水域等面積下降。

表5 河口段生態(tài)系統(tǒng)服務價值 百萬元

研究區(qū)生態(tài)系統(tǒng)服務總價值30余a呈波浪變化趨勢,共增加了64199萬元,見表6。小浪底工程運營前期研究區(qū)生態(tài)系統(tǒng)服務價值增大,運營后期小浪底工程調(diào)水調(diào)沙動力下降和人類活動導致研究區(qū)生態(tài)系統(tǒng)服務價值下降,說明小浪底工程運營前期,工程建設對研究區(qū)生態(tài)系統(tǒng)服務價值產(chǎn)生了較大正向影響,工程運營后期沿岸土地開發(fā)強度增大導致工程產(chǎn)生的正向影響被削弱。

表6 研究區(qū)生態(tài)系統(tǒng)服務總價值 百萬元

4 結(jié)論

采用當量因子法利用遙感監(jiān)測數(shù)據(jù)對黃河小浪底工程建設前后30余a庫區(qū)及黃河下游兩岸各10 km岸線的生態(tài)系統(tǒng)服務價值進行了核算,得出以下結(jié)論。

(1)小浪底工程建設期河岸帶農(nóng)田和草地面積減小,但總體呈下降趨勢,研究區(qū)生態(tài)系統(tǒng)服務總價值呈增大趨勢,說明工程建設對河岸帶生態(tài)系統(tǒng)產(chǎn)生的負面影響不大。

(2)工程運營前期,小浪底工程對生態(tài)系統(tǒng)產(chǎn)生了較大的正向影響,研究區(qū)生態(tài)系統(tǒng)服務價值增大。其中,小浪底工程對庫區(qū)上游河岸帶生態(tài)系統(tǒng)服務價值產(chǎn)生的正向影響最大,對下游的影響程度與河岸帶和小浪底工程的距離正相關(guān),主要原因是調(diào)水調(diào)沙提高了黃河下游主河漕過流能力,對下游生態(tài)系統(tǒng)恢復具有正向影響。

(3)工程運營后期,研究區(qū)兩岸人類活動的負面影響削弱了小浪底工程對生態(tài)系統(tǒng)服務價值的正向影響。小浪底工程建設保障了下游防洪安全,也刺激了研究區(qū)開發(fā),進而使得黃河生態(tài)壓力增大。聊城—濱州段、河口段人類活動對生態(tài)系統(tǒng)服務價值的負面影響尤為顯著。