陳 曦，葉可陌，李 坤，金 陽

（1.浙江農(nóng)林大學(xué) 風(fēng)景園林與建筑學(xué)院，浙江 杭州 311300；2.浙江工業(yè)大學(xué) 設(shè)計(jì)與建筑學(xué)院，浙江 杭州 310012）

土地資源作為人類生存發(fā)展的基本要素和載體，根據(jù)功能可劃分為生產(chǎn)空間、生活空間和生態(tài)空間(“三生空間”）[1]。資源型城市與其他類型城市相比，受到城鎮(zhèn)化發(fā)展與資源開采的雙重壓力，生態(tài)空間被生活和生產(chǎn)空間擠壓的時(shí)段更長，“三生空間”沖突更為強(qiáng)烈[2]。并且，由于土地功能不斷變更，土地利用敏感度增加，景觀格局劇烈變化，面臨的生態(tài)安全問題也更為復(fù)雜與緊迫[3]。因此，科學(xué)提升資源型城市生態(tài)安全是促進(jìn)土地可持續(xù)發(fā)展的關(guān)鍵。資源型城市邊緣區(qū)作為城市建設(shè)中最富變化的區(qū)域，不僅涉及城鄉(xiāng)之間的過渡，礦區(qū)與非礦區(qū)之間的矛盾也尤為突出[4]，需給予重點(diǎn)關(guān)注。生態(tài)系統(tǒng)服務(wù)價(jià)值(ESV）和景觀生態(tài)風(fēng)險(xiǎn)(ERI）作為評(píng)估區(qū)域生態(tài)環(huán)境的重要內(nèi)容，與維護(hù)生態(tài)安全密切相關(guān)[5]?，F(xiàn)階段，多數(shù)學(xué)者將生態(tài)系統(tǒng)服務(wù)價(jià)值與景觀生態(tài)風(fēng)險(xiǎn)作為2 個(gè)獨(dú)立的科學(xué)問題，分別從格局、過程、功能等方面開展大量研究[6-7]。但是，近期研究發(fā)現(xiàn)，生態(tài)系統(tǒng)服務(wù)價(jià)值的時(shí)空異質(zhì)性可以體現(xiàn)景觀格局與生態(tài)過程的最終結(jié)果，是修正景觀生態(tài)風(fēng)險(xiǎn)的最佳指示指標(biāo)[8-9]。因此，將生態(tài)系統(tǒng)服務(wù)價(jià)值與景觀生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)合研究，有助于精確分析區(qū)域生態(tài)環(huán)境變化過程[10-11]，同時(shí)也可開展生態(tài)分區(qū)等相關(guān)研究[12-13]，為提升資源型城市邊緣區(qū)生態(tài)安全提供新的切入點(diǎn)。

本研究以典型的資源型城市大慶城市邊緣區(qū)為例，運(yùn)用地學(xué)信息圖譜、時(shí)空立方體模型、空間自相關(guān)分析等方法研究1980—2020 年“三生空間”土地利用、景觀風(fēng)險(xiǎn)生態(tài)系統(tǒng)服務(wù)價(jià)值，深入探討資源型城市人類活動(dòng)對(duì)土地利用變化的影響、“三生空間”耦合協(xié)調(diào)發(fā)展趨勢(shì)及未來發(fā)展路徑。

1 研究方法

1.1 研究區(qū)域

大慶市(45°23′～47°29′N，123°45′～125°47′E）地處黑龍江省西部，松嫩平原中部。東部與綏化市、哈爾濱市相接，南部與吉林省隔江(松花江）相望，西部、北部與齊齊哈爾市接壤。大慶市蘊(yùn)含豐富的石油資源，大慶油田曾是中國最大的油田。本研究以大慶城市總體規(guī)劃(2011—2020 年）為依據(jù)，劃定城市邊緣區(qū)與油田生產(chǎn)區(qū)(簡稱礦區(qū)，圖1）。

圖1 研究區(qū)示意圖Figure 1 Location of the study area

1.2 數(shù)據(jù)來源及處理

1980、1990、2000、2010、2020 年大慶城市邊緣區(qū)土地利用現(xiàn)狀數(shù)據(jù)來源于地理監(jiān)測云平臺(tái)，空間分辨率為30 m。行政區(qū)劃數(shù)據(jù)源于國家基礎(chǔ)地理信息中心(http：//www.ngcc.cn/ngcc/）。糧食數(shù)據(jù)來自《中國農(nóng)村統(tǒng)計(jì)年鑒 1981—2021 年》和黑龍江省2020 年國民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。本研究根據(jù)研究區(qū)土地利用特點(diǎn)，參考關(guān)于“三生空間”的相關(guān)文獻(xiàn)[14-15]，將土地利用分類系統(tǒng)分為3 個(gè)一級(jí)土地利用空間與8 個(gè)二級(jí)土地利用空間(表1）。

表1 土地利用主導(dǎo)功能分類Table 1 Land use function classification

1.3 研究方法

1.3.1 評(píng)價(jià)單元?jiǎng)澐?為了將研究數(shù)據(jù)進(jìn)行可視化表達(dá)，從空間上更直觀地分析研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值的變化情況，借鑒已有研究和結(jié)合研究區(qū)實(shí)際情況，利用ArcGIS 漁網(wǎng)工具，將研究區(qū)劃分成1 km×1 km 的單元網(wǎng)格(評(píng)價(jià)單元），共得到2 342 個(gè)評(píng)價(jià)單元。

1.3.2 土地利用轉(zhuǎn)移矩陣 以 1980、1990、2000、2010 和2020 年“三生空間”數(shù)據(jù)為基本地理單元，利用ArcGIS 構(gòu)建 1980—1990、1990—2000、2000—2010、2010—2020 年2 位數(shù)編碼地學(xué)信息圖譜。計(jì)算公式如下：Nc=Na×10+Nb。其中，Na和Nb分別表示研究初期a年和末期b年土地“三生空間”二級(jí)類型圖譜單元編碼；Nc為土地轉(zhuǎn)移新生圖譜編碼。

1.3.3 景觀生態(tài)風(fēng)險(xiǎn)估算 景觀生態(tài)風(fēng)險(xiǎn)可以衡量生態(tài)系統(tǒng)受到的外界干擾強(qiáng)度和內(nèi)部抵抗能力[16]。計(jì)算公式如下：

其中：ERIk為第k個(gè)評(píng)價(jià)單元土地利用景觀生態(tài)風(fēng)險(xiǎn)；Aki為第k個(gè)評(píng)價(jià)單元土地利用類型i的面積；Ak為第k個(gè)評(píng)價(jià)單元的總面積；Ri為土地利用類型i的景觀損失度指數(shù)；Ei為土地利用類型i的景觀干擾度指數(shù)；Ci、Si和Di分別為景觀破碎度指數(shù)、景觀分離度指數(shù)和景觀優(yōu)勢(shì)度指數(shù)；a1、b1、c1為各景觀指數(shù)權(quán)重，且三者相加為1，根據(jù)參考文獻(xiàn)[17-18]，將3 個(gè)指數(shù)分別賦值為0.502、0.301、0.197；Fi為景觀i類型的景觀脆弱度指數(shù)，根據(jù)參考文獻(xiàn)[19-21]以及結(jié)合研究區(qū)土地利用類型的變化情況，對(duì)8 個(gè)土地利用空間空間賦分(表1），并進(jìn)行歸一化處理。

1.3.4 生態(tài)系統(tǒng)服務(wù)價(jià)值估算 生態(tài)系統(tǒng)服務(wù)功能是人們從生態(tài)系統(tǒng)獲取的效益，包含直接和間接的、有形和無形的效益[22]。計(jì)算公式為：。其中，ESV為研究區(qū)生態(tài)系統(tǒng)服務(wù)功能總價(jià)值(元）；ESVk為第k個(gè)評(píng)價(jià)單元生態(tài)系統(tǒng)服務(wù)功能價(jià)值(元）；Ai為評(píng)價(jià)單元內(nèi)第i類土地利用類型面積(hm2）；n為土地利用類型數(shù)；VCi為單位面積土地利用類型i的生態(tài)系統(tǒng)服務(wù)價(jià)值(元·hm-2）。其中，生態(tài)系統(tǒng)服務(wù)價(jià)值以謝高地等[23]的中國陸地生態(tài)系統(tǒng)服務(wù)單位面積價(jià)值為基礎(chǔ)，輔以文獻(xiàn)[24-25]所提供的比例關(guān)系確定。在確定農(nóng)業(yè)生活空間生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)，從《大慶統(tǒng)計(jì)年鑒》中獲取 1980、1990、2000、2010、2020 年糧食作物種植面積(hm2）及糧食作物總產(chǎn)值(元）數(shù)據(jù)，得出大慶市農(nóng)田提供食物生產(chǎn)服務(wù)的單價(jià)平均值為1 079.49 元·hm-2。按照自然斷點(diǎn)法分類并采用相對(duì)指標(biāo)[26]，從景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值上劃分為低、較低、中、較高、高共5 類區(qū)域。

1.3.5 時(shí)空熱點(diǎn)分析 以時(shí)空立方體模型為基礎(chǔ)[27]，計(jì)算每個(gè)立方體條柱的Getis-Ord Gi*統(tǒng)計(jì)量，探測某特征在時(shí)空尺度的熱點(diǎn)或冷點(diǎn)，結(jié)合M-K 檢驗(yàn)法對(duì)時(shí)空尺度的熱點(diǎn)分析結(jié)果進(jìn)行趨勢(shì)檢驗(yàn)，從而識(shí)別數(shù)據(jù)聚集程度的時(shí)空趨勢(shì)[28]。

1.3.6 空間自相關(guān)分析 空間自相關(guān)分析可以探索某屬性值是否具有空間集聚性和異質(zhì)性[29]。運(yùn)用GEODA 空間分析工具，通過局部空間自相關(guān)分析[30]確定未來“三生空間”的管控區(qū)域與發(fā)展方向。

2 結(jié)果與分析

2.1 “三生空間”數(shù)量變化及類型轉(zhuǎn)移

2.1.1 “三生空間”分布特征 研究區(qū)總面積為2 342 km2，其中礦區(qū)面積為639.52 km2，占研究面積的23.86%。根據(jù)1980—2010 年土地面積統(tǒng)計(jì)可知(表2），研究區(qū)生態(tài)空間面積較廣，各類型面積從大到小依次為草地、其他、水域、林地。生產(chǎn)空間次之，各類型面積從大到小依次為農(nóng)業(yè)、工業(yè)。生活空間最少，各類型面積從大到小依次為鄉(xiāng)村、城市。根據(jù)礦區(qū)內(nèi)外“三生空間”比例分布可知，農(nóng)業(yè)生產(chǎn)空間、草地生態(tài)空間礦區(qū)內(nèi)外分布較為均等；城市生活空間、水域生態(tài)空間、林地生態(tài)空間分布礦區(qū)外大于礦區(qū)內(nèi)；鄉(xiāng)村生活空間、工業(yè)生產(chǎn)空間、其他生態(tài)空間分布礦區(qū)外小于礦區(qū)內(nèi)。

表2 1980—2020 年研究區(qū)“三生空間”面積與分布比例統(tǒng)計(jì)Table 2 Statistics on the area and distribution ratio of the “production-living-ecological space” in the study area in 1980-2020

2.1.2 “三生空間”土地利用轉(zhuǎn)移 由表3 可知：1980—2020 年，生產(chǎn)空間與生態(tài)空間之間土地轉(zhuǎn)移最為活躍，面積為418.44 km2，生活與生態(tài)空間之間轉(zhuǎn)移次之，面積為88.11 km2，生產(chǎn)與生活之間轉(zhuǎn)移最少，面積為25.56 km2。其中，1980—1990 年，主要為礦外生態(tài)與生產(chǎn)空間之間相互轉(zhuǎn)移；1990—2000 年，主要為礦外生態(tài)空間轉(zhuǎn)移為生產(chǎn)空間；2000—2010 年，主要為礦區(qū)外生產(chǎn)與生態(tài)空間之間相互轉(zhuǎn)移，及礦區(qū)內(nèi)生活空間轉(zhuǎn)移為生態(tài)空間；2010—2020 年，主要為礦外生產(chǎn)與生態(tài)空間相互轉(zhuǎn)移，及礦內(nèi)生態(tài)空間轉(zhuǎn)移為生活空間。

表3 1980—2020 年研究區(qū)“三生空間”面積轉(zhuǎn)移與分布比例統(tǒng)計(jì)Table 3 Statistics on the area transfer and distribution ratio of the “production-living-ecological space” in the study area in 1980-2020

2.2 風(fēng)險(xiǎn)-價(jià)值時(shí)空分析

2.2.1 風(fēng)險(xiǎn)-價(jià)值總體變化趨勢(shì) 研究區(qū)總體景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值指數(shù)呈波動(dòng)上升態(tài)勢(shì)(圖2）。研究期初風(fēng)險(xiǎn)與價(jià)值指數(shù)較低，1980—1990 年變化幅度最大，呈現(xiàn)大幅增長趨勢(shì)。從“三生空間”角度分析可知，景觀生態(tài)風(fēng)險(xiǎn)指數(shù)：生活空間與生產(chǎn)空間呈上升趨勢(shì)，生態(tài)空間呈下降趨勢(shì)?！叭臻g”生態(tài)系統(tǒng)服務(wù)價(jià)值指數(shù)均呈上升趨勢(shì)。

圖2 1980—2020 年研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)和生態(tài)系統(tǒng)服務(wù)價(jià)值變化趨勢(shì)Figure 2 Changing trend of landscape ecological risk and ecosystem service value in 1980-2020

2.2.2 風(fēng)險(xiǎn)和價(jià)值空間分布特征 研究區(qū)1980—2020 年景觀生態(tài)風(fēng)險(xiǎn)的分布格局圖(圖3）與生態(tài)系統(tǒng)服務(wù)價(jià)值分布格局圖(圖4）顯示：低風(fēng)險(xiǎn)區(qū)主要由生活與生態(tài)空間組成，較低風(fēng)險(xiǎn)區(qū)主要由生產(chǎn)與生態(tài)空間組成；高與較高風(fēng)險(xiǎn)區(qū)主要位于礦區(qū)內(nèi)部，其中較高風(fēng)險(xiǎn)區(qū)主要由生產(chǎn)與生態(tài)空間組成，高風(fēng)險(xiǎn)區(qū)主要由生態(tài)空間組成。中價(jià)值區(qū)主要由礦區(qū)外生產(chǎn)用地組成，低與較低價(jià)值區(qū)主要由礦區(qū)內(nèi)部生態(tài)空間與生產(chǎn)空間組成，高與較高價(jià)值區(qū)主要由礦區(qū)外部生態(tài)空間組成。

圖3 1980—2020 年研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)分布示意圖Figure 3 Landscape ecological risk distribution map in 1980-2020

圖4 1980—2020 年研究區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值分布示意圖Figure 4 Ecosystem service value distribution map in 1980-2020

2.2.3 風(fēng)險(xiǎn)和價(jià)值時(shí)空熱點(diǎn)分析 研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空熱點(diǎn)分析圖(圖5）顯示：冷、熱點(diǎn)分別由連續(xù)、持續(xù)、加強(qiáng)、減弱、新增、分散共6 個(gè)部分組成。根據(jù)研究區(qū)景觀生態(tài)風(fēng)險(xiǎn)時(shí)空熱點(diǎn)分析可知：持續(xù)與連續(xù)冷點(diǎn)主要為水域生產(chǎn)空間、城市生活空間與礦區(qū)內(nèi)部工業(yè)生產(chǎn)空間組成。新增冷點(diǎn)主要為礦區(qū)外部其他生態(tài)空間向農(nóng)村生活空間轉(zhuǎn)移。分散冷點(diǎn)主要為生態(tài)用地與生產(chǎn)用地的過渡區(qū)，減弱冷點(diǎn)為礦區(qū)外生態(tài)空間向生產(chǎn)空間轉(zhuǎn)移，加強(qiáng)冷點(diǎn)為礦區(qū)外生產(chǎn)空間向生活空間轉(zhuǎn)移。持續(xù)與連續(xù)熱點(diǎn)主要為其他生態(tài)空間中的鹽堿地組成，新增熱點(diǎn)為礦區(qū)外生態(tài)用地向生產(chǎn)用地轉(zhuǎn)移，分散熱點(diǎn)為生態(tài)用地與生產(chǎn)用地的過渡區(qū)，減弱熱點(diǎn)為礦區(qū)外草地生態(tài)空間內(nèi)部轉(zhuǎn)移。

圖5 景觀生態(tài)風(fēng)險(xiǎn)和生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空熱點(diǎn)聚類示意圖Figure 5 Spatial and temporal hotspot analysis of landscape ecological risk and ecosystem service

根據(jù)生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空熱點(diǎn)分析可知：持續(xù)與連續(xù)熱點(diǎn)主要為生態(tài)空間，由高覆蓋度草地與水域組成。新增與加強(qiáng)的熱點(diǎn)主要為礦區(qū)外部生態(tài)空間內(nèi)部轉(zhuǎn)移，體現(xiàn)為草地生態(tài)空間向水域生態(tài)空間轉(zhuǎn)移。分散熱點(diǎn)主要為生態(tài)用地與生產(chǎn)用地的過渡區(qū)，減弱熱點(diǎn)為礦區(qū)外生態(tài)空間向生產(chǎn)空間轉(zhuǎn)移。持續(xù)與連續(xù)冷點(diǎn)主要由其他生態(tài)用地中的鹽堿地組成，新增冷點(diǎn)為礦區(qū)內(nèi)生態(tài)用地向生產(chǎn)用地的轉(zhuǎn)移，分散冷點(diǎn)為生態(tài)用地與生產(chǎn)用地的過渡區(qū)，加強(qiáng)冷點(diǎn)為礦區(qū)內(nèi)生態(tài)空間內(nèi)部轉(zhuǎn)移體現(xiàn)為草地生態(tài)空間向其他生態(tài)空間轉(zhuǎn)移。

2.2.4 風(fēng)險(xiǎn)和價(jià)值關(guān)聯(lián)分析 利用局部空間自相關(guān)分析，將研究區(qū)空間關(guān)系分為高風(fēng)險(xiǎn)-高價(jià)值區(qū)、低風(fēng)險(xiǎn)-低價(jià)值區(qū)、低風(fēng)險(xiǎn)-高價(jià)值區(qū)、高風(fēng)險(xiǎn)-低價(jià)值區(qū)和不顯著區(qū)共5 種(圖6），結(jié)果顯示：低風(fēng)險(xiǎn)-高價(jià)值區(qū)域主要由礦區(qū)外水域生態(tài)空間組成，是未來生態(tài)保護(hù)的重點(diǎn)區(qū)域。高風(fēng)險(xiǎn)-低價(jià)值區(qū)域主要由礦區(qū)內(nèi)其他生態(tài)空間(鹽堿地）組成，是未來生態(tài)修復(fù)的重點(diǎn)區(qū)域。這2 類空間是維護(hù)區(qū)域生態(tài)安全穩(wěn)定的核心區(qū)域，需要重點(diǎn)管控。研究區(qū)空間關(guān)系主要變化年限為1980—1990 年，其原因主要為礦區(qū)外草地生態(tài)空間向水域生態(tài)空間的轉(zhuǎn)移，以及礦區(qū)內(nèi)部草地生態(tài)空間向其他生態(tài)空間的轉(zhuǎn)移。高風(fēng)險(xiǎn)-高價(jià)值區(qū)域與低風(fēng)險(xiǎn)-低價(jià)值區(qū)域主要為“三生空間”過渡區(qū)域，其中高風(fēng)險(xiǎn)-高價(jià)值區(qū)域主要位于礦區(qū)外，是水域生態(tài)空間與農(nóng)業(yè)生產(chǎn)空間的過渡區(qū)域，低風(fēng)險(xiǎn)-低價(jià)值區(qū)域主要由礦區(qū)內(nèi)部生產(chǎn)與生態(tài)過渡區(qū)域組成，這2 類空間雖然占比較少，變化較為穩(wěn)定，但是是影響區(qū)域未來生態(tài)安全的主要風(fēng)險(xiǎn)源。因此，未來應(yīng)將生產(chǎn)空間作為主要風(fēng)險(xiǎn)源進(jìn)行科學(xué)防范。

圖6 1980—2020 年景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)價(jià)值關(guān)系聚類示意圖Figure 6 Cluster map of the relationship between landscape ecological risk and ecosystem value in 1980-2020

3 討論

3.1 人類擾動(dòng)因子對(duì)“三生空間”土地利用的影響

資源型城市“三生空間”土地利用變化的影響機(jī)制與人類活動(dòng)關(guān)系密切，通過探討人類擾動(dòng)因子對(duì)“三生空間”土地利用的影響，有助于未來土地利用的決策制定。

研究區(qū)生產(chǎn)空間變化主要受到產(chǎn)業(yè)產(chǎn)值的影響。生活空間變化主要與人口因素相關(guān)，由于城市邊緣區(qū)經(jīng)濟(jì)發(fā)展較為落后，與中心城區(qū)距離較近，使得農(nóng)村人口多去城市務(wù)工，邊緣區(qū)農(nóng)村人口數(shù)量流失明顯。生態(tài)空間變化主要與油田開采相關(guān)，油田開采破壞地表植被使得土地裸露，并且開采過程中消耗大量地下水資源，引起地表沉降形成多個(gè)湖泊的同時(shí)也加重土地鹽堿程度。這與其他資源型城市研究結(jié)果相似[31]。相比之下，非資源型城市的生產(chǎn)空間變化主要受到經(jīng)濟(jì)發(fā)展水平的影響，生活空間變化主要與城市規(guī)劃和管理相關(guān)，生態(tài)空間變化則主要與環(huán)境保護(hù)政策相關(guān)[32]。

3.2 “三生空間”耦合特征

景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值的耦合性能夠間接反映出區(qū)域生態(tài)系統(tǒng)中“三生空間”的協(xié)調(diào)性，探討生態(tài)、生產(chǎn)、生活空間之間的耦合協(xié)調(diào)關(guān)系對(duì)于維護(hù)資源型城市邊緣區(qū)生態(tài)安全具有重要意義。

資源型城市受到資源不可再生性影響，一般要經(jīng)歷興起期、成長期、繁榮期、衰退期(再生期） 4 個(gè)階段。與非資源型城市相比，其“三生空間”耦合協(xié)調(diào)性不僅受到城市規(guī)劃建設(shè)的影響[33]，也受資源發(fā)展周期影響[34]。1980—2020 年間，研究區(qū)“三生空間”呈由失調(diào)向協(xié)調(diào)轉(zhuǎn)變的趨勢(shì)，礦區(qū)內(nèi)“三生空間”耦合協(xié)調(diào)程度低于礦區(qū)外，這與其他資源型城市研究結(jié)果[35]相似。2000 年前，大慶市處于興起期與成長期，這一時(shí)段大慶市進(jìn)行大規(guī)模的城市建設(shè)與油田開采活動(dòng)，使得“三生空間”呈現(xiàn)“高耦合、低協(xié)調(diào)”的極度失調(diào)特征。2000 年后，大慶市進(jìn)入繁榮期與衰退期，這一時(shí)期城市建設(shè)逐漸穩(wěn)步，同時(shí)石油儲(chǔ)備降低、開采活動(dòng)減弱，“三生空間”呈現(xiàn)協(xié)調(diào)發(fā)展特征。

3.3 區(qū)域生態(tài)安全格局構(gòu)建

1980—2020 年，研究區(qū)景觀生態(tài)系統(tǒng)服務(wù)價(jià)值雖然有所改善，但景觀生態(tài)風(fēng)險(xiǎn)處于上升趨勢(shì)，說明生態(tài)安全仍面臨一定壓力。根據(jù)風(fēng)險(xiǎn)-價(jià)值的關(guān)系分析結(jié)果，結(jié)合大慶市國土空間總體規(guī)劃(2021—2035 年），從“三生空間”視域協(xié)調(diào)土地關(guān)系推動(dòng)構(gòu)建生態(tài)安全格局：①劃定城市生態(tài)修復(fù)區(qū)，將礦區(qū)內(nèi)其他生態(tài)空間(鹽堿地）作為重點(diǎn)修復(fù)區(qū)域。要特別注重礦區(qū)內(nèi)部生態(tài)環(huán)境保護(hù)，做到石油開采與環(huán)境建設(shè)兼顧，提高礦區(qū)生態(tài)環(huán)境質(zhì)量。②保護(hù)城市水源地與濕地，將礦區(qū)外部水域生態(tài)空間作為核心保護(hù)區(qū)域，對(duì)于該區(qū)生態(tài)保護(hù)策略應(yīng)側(cè)重原有景觀的維護(hù)，將其作為城市重要的生態(tài)斑塊進(jìn)行管理，助力構(gòu)建“一核一區(qū)一帶、多廊多點(diǎn)”的城市生態(tài)安全格局。③生產(chǎn)空間作為大慶城市邊緣區(qū)主要的風(fēng)險(xiǎn)源，禁止農(nóng)業(yè)生產(chǎn)空間侵占生態(tài)空間，防止工業(yè)生產(chǎn)空間無序擴(kuò)張，建立風(fēng)險(xiǎn)預(yù)警機(jī)制，強(qiáng)化生態(tài)安全格局屏障。

4 結(jié)論

本研究基于 1980、1990、2000、2010 和 2020 年的土地利用現(xiàn)狀分類圖，從“三生空間”的視域，以大慶市邊緣區(qū)為例，定量研究資源型城市景觀生態(tài)風(fēng)險(xiǎn)指數(shù)及生態(tài)系統(tǒng)服務(wù)價(jià)值時(shí)空變化特征，并利用空間回歸模型進(jìn)一步分析兩者之間的關(guān)系，主要結(jié)論如下：①根據(jù)“三生空間”土地分布與轉(zhuǎn)移分析可知：研究區(qū)生態(tài)空間面積較廣，生產(chǎn)空間次之，生活空間最少。礦區(qū)外土地轉(zhuǎn)移活躍度大于礦區(qū)內(nèi)，其中，生產(chǎn)空間與生態(tài)空間之間轉(zhuǎn)移最為活躍；②根據(jù)風(fēng)險(xiǎn)與價(jià)值的時(shí)空特征分析可知，1980—2020 年，研究區(qū)“三生空間”景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值呈上升趨勢(shì)。根據(jù)時(shí)空熱點(diǎn)分析可知，景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值的高低分布有明顯的區(qū)域性特點(diǎn)且較為穩(wěn)定。其中，礦區(qū)內(nèi)風(fēng)險(xiǎn)高于礦區(qū)外，礦區(qū)內(nèi)價(jià)值低于礦區(qū)外。③根據(jù)景觀生態(tài)風(fēng)險(xiǎn)與生態(tài)系統(tǒng)服務(wù)價(jià)值關(guān)聯(lián)分析可知：風(fēng)險(xiǎn)與價(jià)值之間存在顯著相關(guān)性。其中，低風(fēng)險(xiǎn)高價(jià)值區(qū)域主要由礦區(qū)外水域生態(tài)空間組成，高風(fēng)險(xiǎn)低價(jià)值區(qū)域主要由礦區(qū)內(nèi)其他生態(tài)空間組成。

綜上所述，未來大慶城市邊緣區(qū)在生態(tài)環(huán)境建設(shè)過程中應(yīng)以縮小礦區(qū)內(nèi)外差異為首要任務(wù)，降低人為因素(特別是油田開采）對(duì)“三生空間”的擾動(dòng)，加強(qiáng)各區(qū)政府與油田管理局之間的協(xié)作機(jī)制，共同統(tǒng)籌區(qū)域“三生空間”建設(shè)，重點(diǎn)保護(hù)礦區(qū)外水域生態(tài)空間，修復(fù)礦區(qū)內(nèi)鹽堿地，并防止生產(chǎn)空間進(jìn)一步擴(kuò)張。