張亮霞　劉利霞

摘 要：地理學的研究內(nèi)容復雜，其要素相互聯(lián)系，過程多為不可控制，不同于物理學和化學。定量的解釋地理問題是地理學的發(fā)展趨勢，運用統(tǒng)計相關(guān)模型來進行地理建模已經(jīng)成為非常普遍的研究手段。該文通過概述植物響應氣候變化、空氣污染物與氣象因子之間的關(guān)系、人口數(shù)據(jù)統(tǒng)計分析3個方面的應用，大致介紹了統(tǒng)計相關(guān)模型在地理建模方面的發(fā)展，同時由于統(tǒng)計模型優(yōu)于數(shù)據(jù)分析而弱于邏輯推理，仍需不斷改進完善。

關(guān)鍵詞：統(tǒng)計模型；地理建模；動態(tài)變化；回歸分析；多格網(wǎng)尺度

中圖分類號 P208 文獻標識碼 A 文章編號 1007-7731（2018）05-0012-03

1 引言

地理學主要研究人地關(guān)系、空間格局、地理現(xiàn)象的時空演變等問題。這些問題覆蓋多種因素，且多為不可控制，明顯區(qū)別于物理學和化學[1]。針對這些問題，對應于具體實際問題，建立可解釋的模型，對地理現(xiàn)象、地理機理與過程進行總結(jié)與表達，以實現(xiàn)地理模擬，為社會、經(jīng)濟的決策提供服務(wù)。地理建模是人類固有的、自然的、不可避免的活動[2-3]。地理建模通過運用邏輯分析法、物理學方法、統(tǒng)計分析法、模糊數(shù)學法、動態(tài)數(shù)據(jù)分析、空間數(shù)據(jù)分析等方法對不同問題進行分析，能夠得到有效、合理的解釋。統(tǒng)計相關(guān)模型是通過用實驗數(shù)據(jù)和調(diào)查來替代具有隨機不確定性的地理現(xiàn)象特征，運用統(tǒng)計學的數(shù)學方法，呈現(xiàn)出可視化、易分析的模型，結(jié)合實際問題，進行現(xiàn)象解釋、變量間研究、預測及數(shù)據(jù)、結(jié)構(gòu)簡化等工作。統(tǒng)計相關(guān)模型通過定量的方法能夠?qū)σ恍┑乩韱栴}進行建模和解釋。隨著研究的不斷深入，定量模擬已成為地理建模的發(fā)展趨勢。

2 統(tǒng)計相關(guān)建模應用

2.1 植物響應氣候變化模擬 根據(jù)器測記錄和過去一段時間的溫度重建可知，全球變暖是目前人們公認的全球溫度變化。關(guān)于溫室氣體的增溫效應及幅度大小一直是研究的重點，其中植物對氣候變化的響應成為了學者關(guān)注的熱點[4-7]。目前的研究多數(shù)集中在植物生長狀態(tài)及生物量、種群的構(gòu)成成分及動態(tài)、群落結(jié)構(gòu)及功能、植被分布范圍和界限變化模擬等方面，主要的定量模擬方法包括統(tǒng)計模型、生物地理相關(guān)模型、生態(tài)響應面模型、立地模型、植物生理模型等。其中，統(tǒng)計模型中發(fā)展較為成熟的模型有元胞自動機及馬爾科夫轉(zhuǎn)移矩陣。第一種是元胞自動機，元胞自動機以動態(tài)形式聯(lián)接離散時間、多維空間的空間位置，并同步其變化狀態(tài)。不同元胞間的狀態(tài)和轉(zhuǎn)移函數(shù)將引起單個元胞發(fā)生改變。植物種及其周圍種空間分布格局及其時間變化受到廣泛重視，利用元胞自動機模型的原理和特點，可以模擬這種時空變化[8]。例如Dytham于1995年模擬了2種植物之間的共存；Loh和Hsieh于1995年模擬了薩王納景觀的演替。從過去的例子中，人們發(fā)現(xiàn)元胞自動機模型的轉(zhuǎn)移函數(shù)過于簡單，缺乏生物學和地學意義，因此還需進行反復改進。第二種是馬爾科夫鏈模型，馬爾科夫鏈模型可以通過計算生態(tài)系統(tǒng)成分，包括植被類型、年齡組、空間位置等的狀態(tài)矩陣，來對草地、森林、珊瑚礁、景觀及土地利用動態(tài)等動態(tài)變化進行模擬[9]。但由于單一的馬爾科夫鏈模型僅適用于轉(zhuǎn)移概率不隨時間變化和狀態(tài)不發(fā)生變化的系統(tǒng)，而群落演替過程中物種的種類和數(shù)量往往是不斷發(fā)生變化的，因此，一些學者開始將空間維引入馬爾科夫鏈形成時空的馬爾科夫鏈模型（STMC），考慮了鄰域狀態(tài)對元胞動態(tài)的影響，這樣就可以使數(shù)據(jù)同時具有時間和空間特征，解決了單一馬爾科夫鏈模型不能模擬一些群落演替和預測變化的弊端[10]。

2.2 北京地區(qū)PM2.5模型 城市污染問題是城市生態(tài)發(fā)展中的一個重要部分。城市污染物中的汽車尾氣、化工等排放的氮氧化物（NOX）、一氧化碳（CO）、二氧化硫（SO2）等都會危害人們的身體健康。污染物的維持、擴散與氣象條件有密切的物理化學關(guān)系：根據(jù)大氣擴散理論和研究表明，不同的氣象條件下，同一污染源排放所造成的近地層污染物體積分數(shù)可相差幾十倍乃至幾百倍[11]；北京地區(qū)NOX、SO2、CO和O3污染氣體濃度變化氣溶膠垂直分布廓線及其濃度峰值層次與邊界層大氣熱力狀況結(jié)構(gòu)（逆溫層特征）存在顯著聯(lián)系[12]。周麗等則利用北京白石橋的污染物觀測資料和北京地區(qū)的氣象常規(guī)資料，運用逐步回歸的統(tǒng)計方法探討了氣象因子與氣溶膠細粒子相互關(guān)系。在這項研究中指出，氣象因子對氣溶膠的濃度分布有很大的影響[13]；北京地區(qū)PM10污染狀況的變化與氣象條件密切相關(guān)，其濃度的變化在穩(wěn)定天氣時與能見度、風速和氣象條件呈反相關(guān)[14]；風速和PM2.5的反相關(guān)特征十分明顯，PM2.5高濃度的數(shù)值發(fā)生在小風或靜風的時候，隨著風速的增大，PM2.5濃度含量減??；氣象因子與PM2.5有較好的相關(guān)關(guān)系。在建立回歸模型過程中，自變量為風速（x1）、氣壓（x2）、氣溫（x3）、濕度（x4），樣本數(shù)量為100個，進行8點平滑，擬合對象為PM2.5，擬合99個時次的PM2.5，并對4個自變量進行檢驗。結(jié)果表明，對污染物影響最大的是氣壓，其次是風速，兩者與PM2.5的變化反相關(guān)，溫度和濕度對PM2.5的影響不顯著，氣壓和風速可以作為預測氣溶膠細粒子濃度的重要參考因子；獨立樣本預報試驗結(jié)果顯示雖然預測的氣溶膠粒子濃度與觀測值不能完全吻合，但擬合曲線與觀測曲線變化趨勢相同，說明PM2.5粒子濃度的變化與大氣動力結(jié)構(gòu)及其近地層風場擴散能力密切相關(guān)[15-18]。

2.3 人口數(shù)據(jù)多尺度空間化模型 人類活動和人口增長是造成生態(tài)環(huán)境發(fā)生變化的重要驅(qū)動力之一。人口空間分布可以直觀、便捷的反映人類活動。隨著研究的進步和空間分析技術(shù)的發(fā)展，基于行政單元的統(tǒng)計型人口數(shù)據(jù)對行政單元內(nèi)部的人口分布空間差異描述有限，時間可比性差，難以滿足實際研宄和應用的需求。地理對象能夠呈現(xiàn)出特有的時間、空間乃至時空尺度特征，應采用一定尺寸的地理格網(wǎng)描述統(tǒng)計型人口數(shù)據(jù)[19]。目前發(fā)展的人口數(shù)據(jù)空間化方法包括基于城市地理學人口分布理論規(guī)律的人口密度模型、利用面插值技術(shù)的人口空間化方法、以遙感和GIS空間分析技術(shù)為手段的人口數(shù)據(jù)空間化方法研究、基于空間數(shù)據(jù)挖掘理論的智能化人口數(shù)據(jù)空間化方法等方法[20-23]。有學者以安徽省為研究區(qū)域，利用遙感和GIS技術(shù)，通過多源數(shù)據(jù)融合處理，對該省2010年縣級統(tǒng)計型人口普查數(shù)據(jù)進行了多格網(wǎng)尺度的空間化模型研究，依據(jù)人口特征二次分區(qū)結(jié)果分別進行基于土地利用類型的統(tǒng)計分析建模，引入夜間燈光數(shù)據(jù)進行分級，實現(xiàn)城鎮(zhèn)用地再分類，使用多元統(tǒng)計回歸分析方法和地理加權(quán)回歸方法，進行了人口空間化模型研究[24]。獲取了不同格網(wǎng)尺度的人口空間化數(shù)據(jù)，并結(jié)合縣級、部分鄉(xiāng)鎮(zhèn)人口統(tǒng)計數(shù)據(jù)，運用統(tǒng)計分析和空間分析，對人口數(shù)據(jù)的空間化結(jié)果進行了誤差比較和模型精度分析。研究結(jié)果表明：以1km格網(wǎng)尺度進行的分區(qū)建模的整體平均相對誤差比整體建模誤差相比下降了16.28%，大大提高了建模精度；引入夜間燈光數(shù)據(jù)，可以很好的表現(xiàn)經(jīng)濟發(fā)展水平和人口聚集程度的關(guān)系，使人口空間分布差異性難以區(qū)分的情況得到改善；研究分別獲取了1km、5km、10km 3種格網(wǎng)尺度的人口空間化數(shù)據(jù)，通過直接精度驗證表明，對于多元統(tǒng)計回歸方法，模型精度隨著尺度的增加而增加[25-29]。

3 結(jié)語

在地理學定量發(fā)展的驅(qū)使下，對實際問題的地理建模中，統(tǒng)計相關(guān)模型能夠?qū)⒌乩憩F(xiàn)象、地理過程等定量的表述以便進行分析和解釋。但單純的統(tǒng)計模型則會使某些地理要素的作用弱化，更多情況下這種模型只適合一時一地一事的情況，在有些實例中，不能解釋機理，參數(shù)之間缺乏邏輯性。對于不同地區(qū)、不同條件，往往可以得出多種統(tǒng)計規(guī)律，不能進行時間和空間的外推。可見，統(tǒng)計相關(guān)模型還需要不斷完善，才能更好地利用在地理建模實例中，更合理的用定量方法解釋地理現(xiàn)象和地理問題。

參考文獻

[1]韋玉春，陳鎖忠.地理建模原理與方法[M].北京：科學出版社，2005.

[2]陳旻，盛業(yè)華，溫永寧，等.語義引導的圖標式地理概念建模環(huán)境初探[J].地理研究，2009，03：705-715.

[3]Westervelt J.流域管理的模擬建模[M].程國棟，李新，王書功，譯.鄭州：黃河水利出版社，2004.

[4]Box E O，Crumpacker D W，Hardin E D.Predicted effect s of climatic change on distribution of ecologically important native tree and shrub species in Florida[J].Climatic Change，1999，41：213-248.

[5]Sykes M T，Prentice I C.Quantifying the impact of global climate change on potential natural vegetation[J].Climate Change，1999，41：37-52.

[6]Channell R，Lomolino M V.Dynamic biogeography and conservation of end angered species[J].Nature，2000，403：84-86.

[7]王紹強，羅承文.中國陸地自然植被碳量空間分布特征探討[J].地理科學進展，1999，18（3）：238-244.

[8]Woodward F I，Thomps on G F，McKee I F.The effects of elevated concentrated concentrations of carbon dioxide on individual plants，populations，communities and ecosystems[J].Ann.Bot.，1991，67：23-28.

[9]Childress W M，Crisafulli C M et al.Comparison of Markovian matrix models of a primary successional plant community[J].Ecological Modelling，1998，107：93-102.

[10]李雙成.植物響應氣候變化模型模擬研究進展[J].地理科學進展，2001，03：216-225.

[11]李宗愷，潘云仙，孫潤橋.空氣污染氣象學原理及應用[M].北京：氣象出版社，1985.

[12]徐祥德，丁國安，周麗，等.北京地區(qū)冬季大氣污染動力——化學過程區(qū)域性三維結(jié)構(gòu)特征[J].科學通報，2003，48（1）：1-6.

[13]楊東貞，于海青.北京北郊冬季低空大氣氣溶膠分析[J].應用氣象學報，2002，13（特刊）：113-126.

[14]王淑英，張小玲.北京地區(qū)PM10污染的氣象特征[J].應用氣象學報，2002，13（特刊）：177-184.

[15]徐祥德.城市化環(huán)境大氣污染模型動力學問題[J].應用氣象學報，2002，13（特刊）：1-12.

[16]黃嘉佑.氣象統(tǒng)計分析與預報方法[M].北京：氣象出版社，1990，84-86.

[17]魏風英.現(xiàn)代氣候統(tǒng)計診斷與預測技術(shù)[M].北京：氣象出版社，2007.

[18]周麗，徐祥德，丁國安，等.北京地區(qū)PM2.5統(tǒng)計—相關(guān)模型[A]//中國氣象學會.新世紀氣象科技創(chuàng)新與大氣科學發(fā)展——中國氣象學會2003年年會“大氣氣溶膠及其對氣候環(huán)境的影響”分會論文集[C].中國氣象學會，2003：4.

[19]李雙成，蔡運龍.地理尺度的轉(zhuǎn)換若干問題的初步探討[J].地理研究，2005，24（1）：11-18.

[20]Gallopin G-C.Human dimensions of global change：linking the global and local processes[J].International Social Science Journal.1991，43（4）：707-718.

[21]Clark C.Urban population densities[J].Journal of the Royal Statistical Society.Series A（General），1951，114（4）：490-496.

[22]潘志強，劉高煥.面插值的研宄進展[J].地理科學進展，2002，21（2）：146-152.

[23]王雪梅，李新，馬明國.基于遙感和GIS的人口數(shù)據(jù)空間化研宄進展及案例分析[J].遙感技術(shù)與應用，2004，19（5）：320-327.

[24]黃河清，王有亮，胡寶清，等.基于神經(jīng)網(wǎng)絡(luò)和GIS的廣西都安縣人口數(shù)據(jù)空間化研究[J].測繪與空間地理信息，2009，32（6）：46-49.

[25]曹麗琴，李平湘，張良培.基于DMSP/OLS夜間燈光數(shù)據(jù)的城市人口估算——以湖北省各縣市為例[J].遙感信息，2009，1：83-87.

[26]董春，劉紀平，趙榮，等.地理因子與空間人口分布的相關(guān)性研究[J].遙感信息，2002，4：61-64.

[27]Harvey J T.Estimating census district populations from satellite imagery：some approaches andlimitations[J].International Journal of Remote Sensing，2002，23（10）：2071-2095.

[28]Martin D.Mapping population date from zone centroid locations[J].Transactions of the Institute of British Geographers，1989，14（1）：90-97.

[29]王珂靖.統(tǒng)計型人口數(shù)據(jù)多尺度空間化模型研究[D].上海：華東師范大學，2015.

[30]仝慧杰.森林生物量遙感反演建?；A(chǔ)與方法研究[D].北京：北京林業(yè)大學，2007. （責編：張宏民）