基于Ann-CA-Markov模型的生態(tài)空間預測模擬: 以重慶市萬州區(qū)為例

幸瑞燊,周啟剛,4①

(1.生態(tài)環(huán)境空間信息數(shù)據(jù)挖掘與大數(shù)據(jù)集成重慶市重點實驗室，重慶 401320；2.重慶財經(jīng)學院訊飛人工智能學院，重慶 401320；3.重慶工商大學環(huán)境與資源學院，重慶 400067；4.重慶工商大學公共管理學院，重慶 400067)

生態(tài)空間指能夠提供重要生態(tài)服務功能、生態(tài)產(chǎn)品并保證國家及區(qū)域生態(tài)安全的空間[1]，其不僅是國土空間重要的一部分，而且是生態(tài)系統(tǒng)的核心載體[2]。隨著城市化進程的快速發(fā)展，“三生”空間之間的利用矛盾日益突出，我國對生態(tài)系統(tǒng)服務功能保護、恢復以及生態(tài)空間的可持續(xù)發(fā)展與保護進入關鍵時期，生態(tài)空間成為熱點問題。目前國內(nèi)外學者對于生態(tài)空間的研究主要集中于生態(tài)空間的演化特征及其梯度效應[3-4]、生態(tài)空間的識別[5]、生態(tài)空間的重構(gòu)優(yōu)化[6]以及生態(tài)空間的保護和利用[7]，但是鮮有對生態(tài)空間模擬預測的研究。未來的生態(tài)空間規(guī)模和格局直接影響著國土空間的生態(tài)安全，故構(gòu)建一個科學、合理的模型對生態(tài)空間進行模擬預測，對生態(tài)環(huán)境的保護以及國土空間的優(yōu)化管控具有重要的理論和現(xiàn)實意義。

在對于空間規(guī)模和格局的模擬預測研究中，國內(nèi)外學者主要使用各種模型對土地利用格局進行模擬預測，使用最為廣泛的模型包括元胞自動機(CA)[8-9]、人工神經(jīng)網(wǎng)絡(Ann)[10]、土地利用變化及效應(CLUE-S)[11-12]等模型，許多學者還構(gòu)建復合模型進行模擬預測[13-15]。這些模型在模擬預測過程中都有著各自的優(yōu)勢，但是同樣存在著不足之處：CA模型在模擬過程中對各種宏觀因素的考慮明顯不足；Ann模型對網(wǎng)絡結(jié)構(gòu)選擇不一，且存在局部極小化等問題；CLUE-S 模型對于自然和政策等要素的影響考慮不足，過多強調(diào)經(jīng)濟效益的影響。因此，尋找一種能夠快速、科學、合理的模型對生態(tài)空間進行模擬預測成為亟待解決的問題。CA-Markov模型是國內(nèi)外學者構(gòu)建最多的模擬預測模型，其不僅能夠?qū)Ω鞣N生態(tài)空間類型的轉(zhuǎn)移總量以及轉(zhuǎn)移概率矩陣進行預測，而且能夠通過鄰域之間的分析對生態(tài)空間類型的空間分布情況進行模擬預測[16]，但是生態(tài)空間與其影響因子之間是非常復雜的非線性關系，CA-Markov模型很難確定它們之間的關系確定轉(zhuǎn)換規(guī)則，Ann模型則極為擅長通過學習的方式處理復雜的非線性關系，因此，將Ann與CA-Markov模型結(jié)合，彌補了各個模型之間的缺點，對于生態(tài)空間的模擬預測具有較大的研究價值。

筆者以萬州區(qū)作為研究區(qū)，對萬州區(qū)生態(tài)空間的變化特征進行分析，并通過構(gòu)建Ann-CA-Markov耦合模型對萬州區(qū)未來生態(tài)空間的規(guī)模及格局進行模擬預測，同時對萬州區(qū)生態(tài)空間的生境質(zhì)量進行評估，以期為相關的研究提供方法參考，同時也為生態(tài)環(huán)境保護政策制定以及國土空間的優(yōu)化管控提供參考依據(jù)。

1 研究區(qū)概況

萬州區(qū)屬重慶市轄區(qū)，位于長江上游地區(qū)、重慶東北部，處三峽庫區(qū)腹心，幅員面積3 457 km2，城區(qū)面積100 km2，城區(qū)人口102萬人，戶籍總?cè)丝?73.56萬人，常住人口為165.01萬人。其境內(nèi)山丘起伏，少平壩和臺地，經(jīng)濟發(fā)展十分迅速，其地區(qū)生產(chǎn)總值在2019年達920.91億元，不僅是成渝城市群沿江城市帶區(qū)域中心城市、成渝經(jīng)濟區(qū)的東向開放門戶，還是“一帶一路”和長江經(jīng)濟帶重要節(jié)點城市。

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

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

采用的主要數(shù)據(jù)是自三峽工程修建蓄水以來萬州區(qū)2000、2006、2012、2018年4期Landsat TM/ETM遙感影像、DEM影像、交通道路、河流、行政中心、生態(tài)保護紅線和土地利用數(shù)據(jù)，以GB/T 21010—2007《土地利用現(xiàn)狀分類》為標準解譯，通過外業(yè)調(diào)查的方式保證解譯的準確性，如表1所示。

2.2 生態(tài)空間界定

綜合國內(nèi)外對于生態(tài)空間概念的研究發(fā)現(xiàn)，國內(nèi)外學者對于生態(tài)空間的概念界定的最大爭議在于是否將農(nóng)用地納入生態(tài)空間，部分學者認為農(nóng)用地以經(jīng)濟產(chǎn)出為主，不應該納入生態(tài)空間之內(nèi)，而另外一部分學者則認為農(nóng)用地雖然以經(jīng)濟產(chǎn)出為主，但是其同樣具有生態(tài)服務功能，應該納入生態(tài)空間之內(nèi)[7]。2017年中共中央辦公廳、國務院辦公廳印發(fā)了《關于劃定并嚴守生態(tài)保護紅線的若干意見》，指出生態(tài)空間是具有自然屬性、以提供生態(tài)服務或生態(tài)產(chǎn)品為主體功能的國土空間，包括森林、草原、濕地、河流、湖泊、灘涂、岸線、海洋、荒地、荒漠、戈壁、冰川、高山凍原、無居民海島等。筆者以《關于劃定并嚴守生態(tài)保護紅線的若干意見》為依據(jù)，界定該研究的生態(tài)空間包括林地、草地、水域和未利用地。

表1 數(shù)據(jù)來源

2.3 研究方法

2.3.1景觀格局指數(shù)

景觀格局指數(shù)主要是對景觀自身的結(jié)構(gòu)和分布特征進行定量的描述，通過計算景觀格局指數(shù)對生態(tài)空間的空間變化特征進行研究[17-18]，主要選取斑塊數(shù)量(對生態(tài)空間格局進行描述)、斑塊密度(描述某一生態(tài)空間類型的單位面積上的斑塊個數(shù))、最大斑塊指數(shù)(描述各個生態(tài)空間類型斑塊間的連通度)、分離度指數(shù)(描述某一生態(tài)空間類型的分化程度)、香農(nóng)多樣性指數(shù)(描述各個生態(tài)空間類型的均衡化程度)對生態(tài)空間格局的空間變化特征進行描述。

2.3.2InVEST模型

InVEST模型指生態(tài)系統(tǒng)服務評估與權衡模型，主要應用于陸地、淡水和海洋3類生態(tài)系統(tǒng)服務評估，從而支持生態(tài)系統(tǒng)管理和決策的一套模型系統(tǒng)[19]。InVEST生境質(zhì)量模型主要是通過計算生境退化程度進而計算生境質(zhì)量指數(shù)對生境質(zhì)量進行評估。

生境質(zhì)量指數(shù)指生境為某一個體或物種提供生存條件的潛力指數(shù)，通過各個生境的敏感性和威脅因素對生境的威脅強度進行計算[20]，計算公式為

(1)

式(1)中，Qxj為j類型生境柵格x的生境質(zhì)量；Hj為j類型的生境適宜程度；Dxj為j類型生境的柵格x的總威脅水平指數(shù)；z和k為比例因子。

2.3.3Cramer′s V指數(shù)值

在進行生態(tài)空間模擬預測的過程中需要使用影響因子，故通過計算各個影響因子與生態(tài)空間之間的Cramer′s V指數(shù)值(Cv)來測定各個影響因子對生態(tài)空間的影響力，以為生態(tài)空間的模擬預測提供依據(jù)，計算公式為

(2)

式(2)中，X2為卡方檢驗數(shù)值；r為行數(shù)；s為列數(shù)；n為樣本數(shù)。

2.3.4Ann-CA-Markov耦合模型

Ann模型是一種通過大量簡單的神經(jīng)元相互連接來模擬人腦功能的基本特征所形成的一個復雜的非線性動力學習網(wǎng)絡系統(tǒng)，擁有較強的并行、分布式存儲、處理、自組織、自適應和自學能力[21]。BP人工神經(jīng)網(wǎng)絡主要包括輸入層、隱藏層以及輸出層在內(nèi)的3層結(jié)構(gòu)，其本質(zhì)就是通過對訓練樣本的不斷訓練從而獲取最小誤差的網(wǎng)絡參數(shù)值，而后通過輸入與訓練樣本相似的數(shù)據(jù)輸出結(jié)果誤差最小的數(shù)據(jù)[22]。CA模型是一種基于元胞空間以及某種特定的規(guī)則在空間和時間上進行演化模擬的動力學系統(tǒng)[23]。CA模型主要由6部分組成，包括元胞、元胞空間、鄰域和規(guī)則4大最基礎的部分以及狀態(tài)和離散時間2個部分。Markov模型以一類隨機過程為基礎的一種用來預測事件的先驗概率和條件概率的數(shù)學方法[24]。Markov模型由4個部分所組成，分別為狀態(tài)、狀態(tài)轉(zhuǎn)移過程、狀態(tài)轉(zhuǎn)移概率和狀態(tài)轉(zhuǎn)移矩陣。

Ann-CA-Markov耦合模型是以CA作為基礎模型并結(jié)合Markov和Ann模型進行構(gòu)建的，模型的整體結(jié)構(gòu)為CA模型，總共包括元胞、元胞空間、狀態(tài)、鄰域、規(guī)則和循環(huán)次數(shù)6個部分，其中最重要的規(guī)則部分通過Markov和Ann模型計算得到，具體構(gòu)建情況如下。

(1)元胞的定義：生態(tài)空間結(jié)構(gòu)類型圖的大小為30 m×30 m的柵格。

(2)元胞空間的定義：生態(tài)空間結(jié)構(gòu)類型圖的所有柵格。

(3)狀態(tài)的定義：元胞的狀態(tài)就是柵格的屬性，劃分為耕地、林地、草地、水域、未利用地和建設用地。

(4)鄰域的定義：定義5×5的擴展性摩爾鄰域進行研究，即周邊有24個元胞對元胞的屬性產(chǎn)生影響。

(5)規(guī)則的定義：元胞的轉(zhuǎn)換規(guī)則主要包括空間和數(shù)量上的轉(zhuǎn)換規(guī)則，數(shù)量上的轉(zhuǎn)換規(guī)則主要是通過Markov模型進行計算(具體詳見3.2.1節(jié))，空間上的轉(zhuǎn)換規(guī)則通過Ann模型進行計算(具體詳見3.2.2和3.2.3節(jié))。

(6)確定循環(huán)次數(shù)：因使用的每期基礎數(shù)據(jù)之間的間隔均為6 a，故設置確定循環(huán)次數(shù)的年限間隔為6 a。

3 結(jié)果與討論

3.1 研究區(qū)生態(tài)空間變化特征

3.1.1研究區(qū)生態(tài)空間數(shù)量變化特征

根據(jù)表2和圖1可知，萬州區(qū)的生態(tài)空間的面積超過了萬州區(qū)總面積的53%，萬州區(qū)的生態(tài)空間以林地為主，其所占面積超過了萬州區(qū)總面積的49%，其次依次為水域、草地和未利用地。

表2 2000—2018年萬州區(qū)生態(tài)空間結(jié)構(gòu)類型面積統(tǒng)計表

2000—2018年間，萬州區(qū)的生態(tài)空間面積持續(xù)增加，非生態(tài)空間用地中的耕地面積持續(xù)減少、建設用地面積持續(xù)增加。生態(tài)空間面積增加了846 45 hm2，增加幅度占萬州區(qū)總面積的0.24%，其中2006—2012年生態(tài)空間面積增加幅度最為明顯，增加幅度占萬州區(qū)總面積的0.22%。從生態(tài)空間的各個用地類型來看，2000—2018年間林地、草地和未利用地面積持續(xù)減少，其減少幅度分別占萬州區(qū)的總面積的0.36%、0.02%和0.01%；水域面積持續(xù)增加，增加面積2 157.0 5 hm2，增加幅度占萬州區(qū)總面積的0.63%。究其原因在于城鎮(zhèn)以及經(jīng)濟化發(fā)展對林地、草地和未利用地的占用以及三峽庫區(qū)蓄水增加了三峽庫區(qū)的水域面積。

3.1.2研究區(qū)生態(tài)空間在空間上的變化特征分析

為了對生態(tài)空間在空間上的變化特征進行研究，計算萬州區(qū)的斑塊數(shù)量(NP)、斑塊密度(PD)、最大斑塊指數(shù)(LPI)、分離度指數(shù)(SPLIT)和香農(nóng)多樣性指數(shù)(SHDI)(圖2)，在2000—2018年間萬州區(qū)的NP和PD整體上呈增加趨勢，生態(tài)空間的破碎化程度越來越大；LPI持續(xù)減少且幅度較大，最大生態(tài)空間斑塊不斷被侵蝕、破壞；SPLIT持續(xù)增加，各個生態(tài)空間斑塊之間的分化程度越來越大；SHDI持續(xù)增加，各個生態(tài)空間斑塊之間的面積趨于均勻化。由此可見，萬州區(qū)的生態(tài)空間的不穩(wěn)定性正在逐步擴大。

3.2 Ann-CA-Markov耦合模型轉(zhuǎn)換規(guī)則

3.2.1生態(tài)空間類型轉(zhuǎn)移概率矩陣

生態(tài)空間類型的轉(zhuǎn)移概率矩陣是Ann-CA-Markov耦合模型的重要轉(zhuǎn)換規(guī)則之一，為2個不同時期生態(tài)空間類型流轉(zhuǎn)變化的轉(zhuǎn)移概率矩陣。使用萬州區(qū)2000、2006、2012、2018年4期土地利用現(xiàn)狀數(shù)據(jù)，通過馬爾科夫模型分別計算萬州區(qū)在2000—2006、2006—2012和2012—2018年3個時間段內(nèi)的生態(tài)空間類型的轉(zhuǎn)移概率矩陣，如表3～5所示。

表3 2000—2006年萬州區(qū)生態(tài)空間類型轉(zhuǎn)移概率矩陣

3.2.2影響因子的影響力分析

在制作適宜性圖集之前需要對生態(tài)空間變化影響因子的影響力進行分析。通過國內(nèi)外相關研究結(jié)合萬州區(qū)的實際情況，選取高程、道路交通、行政中心、河流、生態(tài)紅線5個影響因子作為Ann-CA-Markov耦合模型的主要影響因子(圖3)，其中高程、道路交通、行政中心、河流為條件因子，生態(tài)紅線為限制因子。

就限制因子而言，生態(tài)紅線是生態(tài)空間內(nèi)具有特殊重要生態(tài)功能、必須強制性嚴格保護的區(qū)域，是保障和維護國家生態(tài)安全的底線和生命線，原則上按照禁止開發(fā)區(qū)域的要求進行管理，嚴禁不符合主體功能定位的各類開發(fā)活動，嚴禁任意改變用途。故在該研究中生態(tài)紅線是生態(tài)空間的關鍵區(qū)域。

就條件因子而言，通過計算Cramer′s V指數(shù)值對各個影響因子的影響力進行解釋，計算結(jié)果如表6所示。根據(jù)表6可知，2000—2018年間道路交通因子和河流因子對萬州區(qū)生態(tài)空間變化的影響力持續(xù)增加，高程因子和其影響力先減小后增加，究其原因在于近年來大量農(nóng)村人口涌入城市，致使城鎮(zhèn)化發(fā)展進程加快且出現(xiàn)大量農(nóng)村山區(qū)地帶耕地撂荒等問題，從而使大量生態(tài)空間被侵占。與此同時，土地復墾、土地整治、退耕還林、退耕還草以及各類扶貧政策、生態(tài)環(huán)境保護政策的實施致使道路交通建設加快、生態(tài)空間內(nèi)的各種自然資源被開發(fā)利用，進而使高程、道路交通、行政中心、河流等條件因子對生態(tài)空間變化的影響力持續(xù)受影響。

表6 各個影響因子影響下的Cramer′s V指數(shù)值

3.2.3生態(tài)空間類型適宜性圖集

生態(tài)空間類型的適宜性圖集是Ann-CA-Markov耦合模型的另外一個重要的轉(zhuǎn)換規(guī)則，主要通過MLP_ANN模型計算生成。選取高程、距道路交通的距離、距行政中心的距離、距河流的距離、現(xiàn)有土地利用類型、鄰近各個土地利用類型的柵格數(shù)量、生態(tài)保護紅線圖等共17個影響因子作為MLP_ANN模型中的輸入層的17個神經(jīng)元，并設置隱藏層的神經(jīng)元為12個，因未利用地本身基數(shù)太小，故設置輸出層的神經(jīng)元一共有5個，分別對應耕地、林地、草地、水域、建設用地的轉(zhuǎn)換概率。通過MLP_ANN模型的學習計算精確度達到了96.44%，故生成的適宜性圖集完全可以用于預測模擬。由于Ann-CA-Markov模型是基于CA模型進行構(gòu)建，最后將修正過后的適宜性圖集進行拉伸，拉伸值為0～255，拉伸后生成的適宜性圖集如圖4所示。

3.3 研究區(qū)生態(tài)空間模擬預測

3.3.1Ann-CA-Markov耦合模型模擬精度檢驗

在完成Ann-CA-Markov耦合模型的構(gòu)建后，對萬州區(qū)2012和2018年的生態(tài)空間進行模擬預測以檢驗Ann-CA-Markov耦合模型的模擬精度，模擬結(jié)果如圖5所示。

主要通過Kappa系數(shù)對Ann-CA-Markov耦合模型的模擬精度進行檢驗，Kappa系數(shù)主要進行圖件之間的一致性驗證，如果計算的Kappa系數(shù)值≥0.75，則認為模擬精度高；如果Kappa系數(shù)值<0.75且≥0.4，則認為模擬精度一般；如果Kappa系數(shù)值<0.4，則認為模擬精度低。

分別計算了萬州區(qū)2012和2018年的標準、數(shù)量、位置和隨機Kappa系數(shù)值，2012年分別為0.983 6、0.987 6、0.992 7、0.992 7，2018年分別為0.985 0、0.988 7、0.996 3、0.996 3，由此可見，筆者所構(gòu)建的Ann-CA-Markov耦合模型具有較高的模擬精度，可以用于萬州區(qū)2024年生態(tài)空間的模擬預測。

3.3.2研究區(qū)2024年生態(tài)空間模擬預測

通過構(gòu)建的Ann-CA-Markov耦合模型(選取的主要影響因子包括高程、道路交通、河流、行政中心以及生態(tài)紅線因子，詳見3.2.2和3.2.3節(jié))對萬州區(qū)2024年的生態(tài)空間格局進行模擬預測，其模擬結(jié)果如圖5和表7所示。

表7 2024年萬州區(qū)生態(tài)空間模擬結(jié)果統(tǒng)計表

根據(jù)圖5和表7的萬州區(qū)2024年生態(tài)空間模擬結(jié)果來看，到2024年萬州區(qū)生態(tài)空間的面加增加，增加幅度占整個萬州區(qū)總面積的0.15%，其中林地、水域和未利用地面積增加，林地和水域面積增加幅度分別占萬州區(qū)總面積的0.18%和0.01%，未利用地增加面積可以忽略不計；草地面積減少，減少幅度占萬州區(qū)總面積的0.04%。

在對2024年萬州區(qū)生態(tài)空間的數(shù)量變化特征進行分析后，進一步對萬州區(qū)2024年生態(tài)空間的空間變化特征進行分析。根據(jù)2024年的生態(tài)空間模擬結(jié)果計算的景觀格局指數(shù)(表8)可知，萬州區(qū)生態(tài)空間的斑塊數(shù)量、斑塊密度、最大斑塊指數(shù)、分離度指數(shù)和香農(nóng)多樣性指數(shù)到2024年都進一步增加，說明萬州區(qū)生態(tài)空間的破碎化程度以及各個生態(tài)空間斑塊之間的分化程度和面積均勻化程度越來越大，萬州區(qū)生態(tài)空間的穩(wěn)定性進一步下降。

表8 2018—2024年萬州區(qū)生態(tài)空間景觀格局指數(shù)表

3.4 研究區(qū)生態(tài)空間生境質(zhì)量評估

為了進一步對萬州區(qū)的生態(tài)空間進行研究，通過InVEST模型對萬州區(qū)2000—2018年以及2024年模擬的生態(tài)空間的生境質(zhì)量進行計算。根據(jù)圖6可知，在整個萬州區(qū)范圍內(nèi)，高生境質(zhì)量區(qū)域(敏感區(qū)域)分布較廣，在萬州區(qū)的各個地方均有分布，集中分布于長江流域以及萬州區(qū)的東南部區(qū)域；而萬州區(qū)的低生境質(zhì)量主要集中分布于萬州區(qū)城區(qū)。

通過ArcGIS 10.1軟件平臺分別計算萬州區(qū)2000、2006、2012、2018和2024年模擬的生態(tài)空間的平均生境質(zhì)量(圖7)。

根據(jù)圖7可知，無論是2000—2018年還是通過模擬的2024年，萬州區(qū)生態(tài)空間的生境質(zhì)量呈現(xiàn)持續(xù)下降的趨勢，其中2000—2006年萬州區(qū)生態(tài)空間的生境質(zhì)量下降最為顯著，2018—2024年下降幅度較小。綜上所述，萬州區(qū)對于生態(tài)空間以及生態(tài)環(huán)境的保護并未使萬州區(qū)的生態(tài)環(huán)境得到較大程度的改善，故需要進一步加強對萬州區(qū)生態(tài)空間的優(yōu)化管控及生態(tài)環(huán)境的保護。

3.5 研究區(qū)生態(tài)空間重構(gòu)策略

通過對萬州區(qū)的生態(tài)空間格局以及生境質(zhì)量的模擬結(jié)果可知，萬州區(qū)的生態(tài)空間減少最主要的原因是城鎮(zhèn)化以及工業(yè)化進程的快速發(fā)展，因此對于控制人口密度以及合理規(guī)劃城市和工業(yè)的發(fā)展是對生態(tài)空間進行合理重構(gòu)的重要策略之一。萬州區(qū)的生境質(zhì)量呈現(xiàn)持續(xù)下降的趨勢，最主要的原因在于高生境質(zhì)量的生態(tài)空間類型向低生境質(zhì)量的生態(tài)空間類型轉(zhuǎn)換，而對于整個萬州區(qū)而言，高生境質(zhì)量生態(tài)空間主要包括生態(tài)紅線以及自然保護區(qū)，低生境質(zhì)量主要包括耕地和未利用地，故對生態(tài)空間進行合理重構(gòu)的另一重要策略是加強高生境質(zhì)量生態(tài)空間的保護(生態(tài)紅線、自然保護區(qū)、林地、水域等)，并在保證區(qū)域糧食安全的情況下將低生境質(zhì)量生態(tài)空間類型向高生境質(zhì)量生態(tài)空間類型轉(zhuǎn)換。因此，生態(tài)空間的管控及優(yōu)化最主要在于對生態(tài)空間總量進行控制、對生態(tài)空間進行分類管制、對生態(tài)空間進行合理布局以及對生態(tài)空間保護實施保障機制。

3.6 討論

筆者所構(gòu)建Ann-CA-Markov耦合模型(Kappa系數(shù)0.983 6)比其他相關研究所使用模型的模擬精度(表9)更高，模擬預測結(jié)果的可信度更高，研究結(jié)果有利于掌握未來生態(tài)空間的變化趨勢，并為相關研究提供方法參考，同時為生態(tài)環(huán)境保護政策的制定以及國土空間的優(yōu)化管控提供參考依據(jù)。但是研究仍然存在著以下2個方面的問題，需要進一步改進。

表9 相關模擬預測研究模型精度

第一方面，該研究使用的TM影像精度較低，且解譯只到一級地類，未精確到二級地類。故在未來的研究中，將盡可能使用高分辨率的遙感影像，獲取更加精細化的生態(tài)空間類型，加強對生態(tài)空間變化特征的研究。

第二方面，Ann-CA-Markov耦合模型創(chuàng)建轉(zhuǎn)換規(guī)則時未設置動態(tài)轉(zhuǎn)換參數(shù)，通過人工神經(jīng)網(wǎng)絡模型創(chuàng)建適宜性圖集時僅設置了靜態(tài)參數(shù)，而未設置動態(tài)參數(shù)，在生態(tài)空間發(fā)生變化的過程中各個影響因子也是發(fā)生變化的。在未來的研究中盡可能地獲取動態(tài)轉(zhuǎn)換的參數(shù)需進一步優(yōu)化Ann-CA-Markov耦合模型，提高模擬預測結(jié)果的可信度。

4 結(jié)論

通過對2000—2018年萬州區(qū)的生態(tài)空間演變特征進行分析，并構(gòu)建Ann-CA-Markov耦合模型對萬州區(qū)2024年的生態(tài)空間進行預測模擬，對萬州區(qū)生態(tài)空間的生境質(zhì)量進行評估，得出以下結(jié)論：

(1)構(gòu)建的Ann-CA-Markov耦合模型具有較高的模擬精度，通過Markov模型和Ann模型對CA模型的主要規(guī)則進行構(gòu)建，解決了CA模型在模擬過程中對各種宏觀因素的考慮不足的問題，使模擬的結(jié)果更加科學、合理、可信。

(2)根據(jù)預測模擬結(jié)果，2000—2024年萬州區(qū)生態(tài)空間面積持續(xù)增加，但是生態(tài)空間的破碎化程度和穩(wěn)定性持續(xù)下降，且生態(tài)空間各個斑塊的分化程度持續(xù)增大，需要從區(qū)域整體上對生態(tài)空間進行優(yōu)化管控。

(3)2000—2024年萬州區(qū)生態(tài)空間的生境質(zhì)量持續(xù)下降，雖然萬州區(qū)生態(tài)空間的數(shù)量不斷增加，但是其生境質(zhì)量仍然持續(xù)下降，生態(tài)空間的優(yōu)化管控不僅需要體現(xiàn)在生態(tài)空間的數(shù)量以及穩(wěn)定性上，更應該關注生態(tài)空間的生境質(zhì)量。