黃佳銘，劉羊，樓康明，譚建輝，陳凌偉

(1.廣州市城市規(guī)劃勘測(cè)設(shè)計(jì)研究院，廣東 廣州 510060； 2.廣東省城市感知與監(jiān)測(cè)預(yù)警企業(yè)重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510060)

1 引 言

花崗巖因球狀風(fēng)化作用形成的微風(fēng)化-中風(fēng)化似球狀堅(jiān)硬塊體通常被稱為孤石[1，2]，根據(jù)其埋藏情況可分為兩種，一種是裸露于地表或局部埋入殘坡積土層中的坡面孤石，另一種是埋藏于全-強(qiáng)風(fēng)化巖或殘積土中的地下孤石[3]。孤石在我國(guó)沿海一帶分布較為廣泛，特別是廣東、福建地區(qū)[4]，地下孤石具有分布不規(guī)律性和不可視性，它的存在對(duì)基坑開挖、隧道掘進(jìn)等重大地下工程建設(shè)產(chǎn)生較大的影響，眾多學(xué)者對(duì)其開展了持續(xù)性研究工作[5～7]，也取得了較為豐富的研究成果。針對(duì)坡面孤石的形成和分布特征，國(guó)內(nèi)外學(xué)者開展了一系列研究工作，國(guó)內(nèi)很多學(xué)者認(rèn)為花崗巖孤石的形成是物理風(fēng)化和化學(xué)風(fēng)化共同作用的結(jié)果[8，9]，國(guó)外很多學(xué)者則認(rèn)為形成花崗巖孤石的主要因素有巖石體積變化、巖石裂縫、卸荷及外部應(yīng)力因素[10～14]；劉治軍等[4]指出花崗巖地區(qū)坡面孤石通常呈集群或線狀分布，亦有少量孤石成孤立零星分布；王平等[1]將坡面孤石分為原生孤石和滾石，原生孤石常以單體形式成群或線狀分布，滾石則多以堆疊形式存在于溝谷或山腳；王浩等[7]闡述了花崗巖孤石的發(fā)育特征，指出了其在工程建設(shè)中可能引發(fā)的地質(zhì)問(wèn)題。

廣州市花崗巖廣泛分布，其中孤石非常發(fā)育，近年來(lái)，因城市建設(shè)和極端強(qiáng)降雨的作用，坡面孤石崩塌災(zāi)害頻繁發(fā)生，嚴(yán)重威脅了廣州市的經(jīng)濟(jì)發(fā)展和城市安全。本文以廣州市黃埔區(qū)為例，基于充足的野外調(diào)查成果，采用信息量模型[15～19]結(jié)合ArcGIS的空間分析功能，對(duì)黃埔區(qū)的坡面孤石發(fā)育程度進(jìn)行分區(qū)評(píng)價(jià)，并對(duì)評(píng)價(jià)結(jié)果進(jìn)行檢驗(yàn)，為研究區(qū)坡面孤石崩塌災(zāi)害的防治工作提供理論依據(jù)和技術(shù)支撐。

2 研究區(qū)概況

黃埔區(qū)位于廣州市中部核心位置，珠江三角洲北緣，行政區(qū)域總面積約 484 km2，占廣州市總面積的6.51%。區(qū)內(nèi)中北大部丘陵山地區(qū)域主要出露三疊系、侏羅系、白堊系地層，巖性從石英閃長(zhǎng)巖到花崗巖，以二長(zhǎng)花崗巖為主，微風(fēng)化花崗巖天然抗壓強(qiáng)度強(qiáng)度大，屬較硬～堅(jiān)硬巖類，在丘陵山地區(qū)域巖石風(fēng)化強(qiáng)烈，殘坡積層較厚，地表和殘坡積層中因風(fēng)化不均勻易發(fā)育花崗巖球中風(fēng)化體。北部鎮(zhèn)龍一帶出露前震旦系云開巖組地層，巖性以片麻巖為主，微風(fēng)化片麻巖天然扛壓強(qiáng)度較大，屬較硬巖類，個(gè)別達(dá)到堅(jiān)硬。南部區(qū)域局部出露古近系地層，巖性以砂礫巖、砂巖、粉細(xì)砂巖為主，微風(fēng)化砂巖天然抗壓強(qiáng)度一般，屬較軟巖，層理較發(fā)育，其中粉細(xì)砂巖、泥巖易風(fēng)化崩解。

3 坡面孤石發(fā)育特征

本次調(diào)查發(fā)現(xiàn)研究區(qū)坡面孤石共計(jì)1 765處，巖性均為花崗巖，中等風(fēng)化-微風(fēng)化，基座大多由殘積土和不同風(fēng)化程度的風(fēng)化巖組成，形狀多樣，埋深和尺寸不一。

3.1 坡面孤石及其基座巖土特征

研究區(qū)坡面孤石主要表現(xiàn)為中微風(fēng)化，花崗巖球狀風(fēng)化由外至內(nèi)進(jìn)行，無(wú)結(jié)構(gòu)面時(shí)風(fēng)化難以在孤石內(nèi)部進(jìn)行，因而，坡面孤石風(fēng)化程度以微風(fēng)化為主，中等風(fēng)化巖次之，強(qiáng)風(fēng)化花崗巖多呈半巖半土狀，難以形成孤石，中等風(fēng)化孤石共399處，占孤石總數(shù)的22.61%，微風(fēng)化孤石共 1 366處，占孤石總數(shù)的77.39%。孤石基座組成可分為殘積土和風(fēng)化巖兩類，在原位風(fēng)化或經(jīng)搬運(yùn)后孤石大概率會(huì)坐落于殘積土之上，少部分坐落于基巖之上，基座為殘積土的坡面孤石共 1 324處，占孤石總數(shù)的75.14%，基座為風(fēng)化巖的坡面孤石共441處，占孤石總數(shù)的24.86%。

3.2 坡面孤石規(guī)模特征

研究區(qū)坡面孤石體積大多不超過(guò) 50 m3，根據(jù)《廣東省地質(zhì)災(zāi)害特征認(rèn)定和分級(jí)標(biāo)準(zhǔn)》中崩塌地質(zhì)災(zāi)害規(guī)模等級(jí)劃分標(biāo)準(zhǔn)，結(jié)合不同尺寸孤石可能造成的危害程度，將其規(guī)模劃分為大、中、小和微型四個(gè)等級(jí)，建立規(guī)模等級(jí)劃分標(biāo)準(zhǔn)表，并依據(jù)標(biāo)準(zhǔn)表對(duì)研究區(qū)坡面孤石規(guī)模等級(jí)進(jìn)行統(tǒng)計(jì)分析(表1)。

表1 坡面孤石規(guī)模等級(jí)劃分標(biāo)準(zhǔn)

從坡面孤石的體積來(lái)看，中型孤石數(shù)量最多，共678處，占坡面孤石總數(shù)的38.41%，微型和小型孤石次之，分別為543處和385處，占坡面孤石總數(shù)分別為30.76%和21.81%，大型孤石最少，共159處，占坡面孤石總數(shù)的9.01%。從坡面孤石的形狀來(lái)看，條塊狀孤石數(shù)量最多，共840處，占孤石總數(shù)的47.59%，橢球體孤石次之，共514處，占孤石總數(shù)的29.12%，不規(guī)則和球體孤石較少，分別為255處和156處，占比分別為14.45%和8.84%。從坡面孤石的埋深來(lái)看，淺埋孤石(出露介于50%～85%)的數(shù)量最多，共 1 565處，占孤石總數(shù)的88.67%，深埋孤石(出露小于50%)和無(wú)埋深孤石(出露大于85%)數(shù)量均較小，分別為47處和153處，占比分別為2.66%和8.67%。

4 坡面孤石分布規(guī)律

從整體來(lái)看，研究區(qū)坡面孤石的分布具有明顯的區(qū)域性，集中分布在聯(lián)和街、蘿崗街和云埔街，分別為561處、458處和370處，分別占坡面孤石總數(shù)的31.78%、25.95%和20.96%，長(zhǎng)嶺街、永和街和穗東街次之，分別為160處、104處和104處，占比分別為9.07%、5.89%和5.89%，其他鎮(zhèn)街零星分布或無(wú)分布。坡面孤石主要發(fā)育于以花崗巖山地區(qū)域，研究區(qū)南部以三角洲平原為主，坡面孤石數(shù)量很少；中部地區(qū)基巖主要為花崗巖，地貌以丘陵、臺(tái)地為主，區(qū)域內(nèi)聯(lián)和街、蘿崗街、云埔街、長(zhǎng)嶺街、永和街等街道坡面孤石很發(fā)育；北部地區(qū)基巖主要為花崗巖和變質(zhì)巖，地貌類型以丘陵為主，但區(qū)域內(nèi)花崗巖成分及年代不易形成孤石，使得研究區(qū)北部較少出現(xiàn)坡面孤石。

從局部來(lái)看，坡面孤石的分布表現(xiàn)出三個(gè)特征：一是坡面孤石密集區(qū)，孤石往往成群分布，表現(xiàn)出明顯的群體分布特征；二是坡面孤石密度較小區(qū)域，由雨水沖刷、人類活動(dòng)等因素影響，導(dǎo)致其分布呈現(xiàn)出線狀或帶狀分布特征；三是坡面孤石被外力搬運(yùn)至此，僅零星分布，坡面孤石以點(diǎn)狀形式存在。研究區(qū)坡面孤石的分布以群體分布為主，點(diǎn)狀分布次之，線狀或帶狀分布較少出現(xiàn)。坡面孤石分布特征如圖1所示。

圖1 坡面孤石分布特征

5 發(fā)育程度評(píng)價(jià)

坡面孤石發(fā)育程度是指在一定自然條件和人類活動(dòng)影響下，坡面孤石形成的難易程度，是通過(guò)孤石形成以及出露于坡面的諸多影響因素進(jìn)行綜合分析，確定坡面孤石在研究區(qū)內(nèi)的形成概率、分布密度等，評(píng)價(jià)坡面孤石的發(fā)育程度，坡面孤石形成及出露于坡面受諸多因素影響，因而，坡面孤石發(fā)育程度是多因子綜合評(píng)價(jià)的結(jié)果。

5.1 評(píng)價(jià)模型

根據(jù)坡面孤石發(fā)育特征及發(fā)育程度評(píng)價(jià)內(nèi)容，綜合考慮選取信息量模型開展研究區(qū)坡面孤石發(fā)育程度評(píng)價(jià)。孤石的形成及出露于坡面受多個(gè)因素的共同影響，信息量模型反映了一定地質(zhì)環(huán)境條件下最易形成坡面孤石及其細(xì)分區(qū)間的組合，通過(guò)特定評(píng)價(jià)單元內(nèi)某種因素作用下坡面孤石形成的概率與區(qū)域坡面孤石調(diào)查成果相比較實(shí)現(xiàn)的，對(duì)應(yīng)某種因素特定狀態(tài)下的坡面孤石發(fā)育的信息量公式可表示為：

式中：IAj-B—對(duì)應(yīng)因素A、j狀態(tài)(或區(qū)間)下坡面孤石發(fā)育的信息量；

Nj—對(duì)應(yīng)因素A、j狀態(tài)(或區(qū)間)下坡面孤石分布的單元格；

N—研究區(qū)已知有坡面孤石分布的單元格；

Sj—對(duì)應(yīng)因素A、j狀態(tài)(或區(qū)間)分布的單元格；

S—工作區(qū)單元總數(shù)。

當(dāng)IAj-B>0時(shí)，反映了對(duì)應(yīng)因素A在j狀態(tài)(或區(qū)間)下坡面孤石發(fā)育的信息量較大，坡面孤石形成的可能性較大，或者說(shuō)利于坡面孤石發(fā)育；當(dāng)IAj-B<0時(shí)，說(shuō)明對(duì)應(yīng)因素A在j狀態(tài)(或區(qū)間)下坡面孤石發(fā)育的信息量較小，坡面孤石形成的可能性較小，或者說(shuō)不利于坡面孤石發(fā)育；當(dāng)IAj-B=0時(shí)，表明對(duì)應(yīng)因素A在j狀態(tài)(或區(qū)間)下不提供坡面孤石形成與否的任何信息，即因素A在j狀態(tài)(或區(qū)間)可以剔除，排除其作為坡面孤石發(fā)育的預(yù)測(cè)因子。

由于每個(gè)評(píng)價(jià)單元受多個(gè)因素的綜合影響，各因素又存在若干狀態(tài)，各因素組合情況下坡面孤石發(fā)育的總信息量可用下式確定：

式中：I—對(duì)應(yīng)特定單元坡面孤石發(fā)育的總信息量，代表坡面孤石發(fā)育的可能性，可作為坡面孤石發(fā)育程度指數(shù)；

Ii—對(duì)應(yīng)特定單元、特定因素第i狀態(tài)(或區(qū)間)條件下坡面孤石發(fā)育的信息量；

Ni—對(duì)應(yīng)特定因素第i狀態(tài)(或區(qū)間)條件下分布有坡面孤石的單元數(shù)；

Si—對(duì)應(yīng)特定因素第i狀態(tài)(或區(qū)間)的單元總數(shù)；

N—調(diào)查區(qū)分布有坡面孤石的單元總數(shù)；

S—調(diào)查區(qū)單元總數(shù)。

5.2 單元?jiǎng)澐?/h3>

調(diào)查發(fā)現(xiàn)，出露型坡面孤石僅發(fā)育于未經(jīng)開發(fā)的山地范圍內(nèi)，平原和建成區(qū)均未發(fā)現(xiàn)有坡面孤石，為評(píng)價(jià)的簡(jiǎn)便性和準(zhǔn)確性，研究區(qū)坡面孤石發(fā)育程度評(píng)價(jià)僅考慮山地范圍，依據(jù)《地質(zhì)災(zāi)害風(fēng)險(xiǎn)調(diào)查評(píng)價(jià)技術(shù)要求(1∶5萬(wàn))》，綜合考慮孤石的分布規(guī)律和發(fā)育特征，將黃埔區(qū)建成區(qū)和平原以外的區(qū)域劃分成 10 m×10 m的評(píng)價(jià)網(wǎng)格單元，共計(jì) 2 737 623個(gè)網(wǎng)格單元。

5.3 評(píng)價(jià)因子

坡面孤石是指獨(dú)立于母巖存在于斜坡之上的單個(gè)花崗巖孤石，發(fā)育程度定義為地質(zhì)環(huán)境條件對(duì)坡面孤石形成及出露于坡面的影響程度，基巖成因類型及時(shí)代是能否形成孤石的必要條件，經(jīng)風(fēng)化剝蝕或重力搬運(yùn)形成的孤石以埋藏型為主，能否出露受諸多因素的影響，包括所處地形坡度、坡向、坡面曲率、高程、斷裂距離、植被覆蓋率等，綜合考慮將巖性、坡度、坡向、坡面曲率、高程、斷裂距離、植被覆蓋率等共計(jì)7個(gè)因素作為評(píng)價(jià)因子。

(1)巖性因子

孤石指花崗巖因差異風(fēng)化而形成的球狀風(fēng)化物，花崗巖或經(jīng)變質(zhì)作用的花崗巖是孤石形成的必要條件，但不同成分、不同年代的花崗巖形成孤石的難易程度不盡相同，根據(jù)工作區(qū)基巖地質(zhì)圖將區(qū)內(nèi)基巖按照成因、年代、組成等進(jìn)行分類，基巖類型可劃分為志留紀(jì)片麻狀花崗巖(S)、元古界云開巖群(Pt1)、元古代片麻雜巖(Pt2)、第三系沉積巖(R)、晚三疊細(xì)中?；◢弾r(T1)、晚三疊中細(xì)粒花崗巖(T2)、中侏羅世中粗?；◢弾r(J1)、晚侏羅世細(xì)粒花崗巖(J2)和白堊世花崗巖(K)。

(2)坡度因子

孤石形成后能否出露于坡面或能否存在于坡面極大程度受所處地形坡度的影響，坡度越大，表層土體的剝蝕越容易，孤石越容易出露于坡面，但當(dāng)坡度過(guò)大時(shí)，孤石出露于坡面后容易失穩(wěn)滾落或滑落，難以存在于坡面，因而，坡面孤石的形成很大程度受坡度控制。根據(jù)工作區(qū)地形起伏程度，將工作區(qū)的坡度劃分為 ≤5°、5°～15°、15°～25°、25°～35°和>35°共五個(gè)等級(jí)。

(3)坡向因子

坡向定義為坡面法線在水平面上的投影方向，坡向?qū)ι降厣鷳B(tài)有較大的作用，它能夠規(guī)律地反映基本的地形特征，坡向?qū)θ照?、太?yáng)輻射和降雨沖刷等有較大影響，一般而言太陽(yáng)輻射南向坡多于北向坡、東向坡多于西向坡，太陽(yáng)輻射間接影響植被的生長(zhǎng)發(fā)育，對(duì)山體表層土體的剝蝕速度有很大影響。研究區(qū)坡向可劃分為平地(-1)、北向(0°～22.5°和337.5°～360°)、北東向(22.5°～67.5°)、東向(67.5°～112.5°)、南東向(112.5°～157.5°)、南向(157.5°～202.5°)、南西向(202.5°～247.5°)、西向(247.5°～292.5°)、北西向(292.5°～337.5°)共九個(gè)級(jí)別。

(4)坡面曲率因子

斜坡表面一定扭曲變化程度定量化的度量因子即為地形曲率，垂直和水平方向上的分量分別稱為平面曲率和剖面曲率，剖面曲率是對(duì)斜坡坡度沿最大坡降方向地面高程變化率的度量，調(diào)查分析表明剖面曲率對(duì)坡面孤石的發(fā)育程度有著重要的影響，剖面曲率一方面影響著坡面土體剝蝕沖刷的速率，另一方面極大地影響著坡面孤石的穩(wěn)定性，根據(jù)研究區(qū)斜坡形態(tài)特征，將邊坡曲率劃分為≤-5、-5～-1.5、-1.5～1.5、1.5～5和>5共五個(gè)等級(jí)。

(5)高程因子

高程的變化是影響坡面孤石發(fā)育的重要因素之一，一方面，高程影響著區(qū)內(nèi)的降雨量、氣溫、植被以及巖層風(fēng)化程度，另一方面，高程的變化影響著人類活動(dòng)強(qiáng)度，地質(zhì)環(huán)境條件破壞程度會(huì)對(duì)坡面孤石的發(fā)育產(chǎn)生影響，因而高程對(duì)坡面孤石發(fā)育程度起到一定的控制作用。由于研究區(qū)高程普遍不大，可將工作區(qū)高程劃分為 <20 m、20 m～50 m、50 m～100 m、100 m～150 m、150 m～200 m和>200 m六個(gè)等級(jí)。

(6)斷裂距離因子

地質(zhì)構(gòu)造對(duì)區(qū)域地形地貌的塑造起著一定的控制作用，對(duì)巖體節(jié)理裂隙的發(fā)育程度有著重要影響，不同構(gòu)造單元對(duì)孤石發(fā)育的難易程度、規(guī)模等具有一定影響。研究區(qū)地質(zhì)構(gòu)造主要體現(xiàn)為斷裂的影響，斷裂對(duì)坡面孤石發(fā)育的影響主要表現(xiàn)在斷裂帶及其附近一定范圍內(nèi)的巖土體結(jié)構(gòu)遭到破壞，降低巖體完整性，加速巖體風(fēng)化剝蝕。根據(jù)斷裂的延伸性、規(guī)模和類型等地質(zhì)特征分析其影響范圍，采用距斷裂距離來(lái)分析其對(duì)坡面孤石發(fā)育程度的影響。研究區(qū)坡面孤石距斷裂的距離可劃分為 <100 m、100 m～300 m、300 m～500 m和>500 m四個(gè)級(jí)別。

(7)植被覆蓋率

植被對(duì)表層土體風(fēng)化剝蝕具有很大的影響，其根系可以提高土體抗剪強(qiáng)度，起到固土護(hù)坡的作用，葉片蒸騰作用可以促進(jìn)地下水的排泄，因而，一般而言植被越茂盛的斜坡表層土體的剝蝕速度越慢。植被指數(shù)是一種利用衛(wèi)星探測(cè)數(shù)據(jù)的線性或非線性組合來(lái)反映植被的存在、數(shù)量、質(zhì)量、狀態(tài)及時(shí)空分布特點(diǎn)的指數(shù)，用歸一化植被指數(shù)(NDVI)表征工作區(qū)植被覆蓋程度，根據(jù)研究區(qū)植被發(fā)育特征，將NDVI劃分為≤0、0～0.1、0.1～0.2、0.2～0.3、0.3～0.4和>0.4共六個(gè)等級(jí)。

5.4 發(fā)育程度評(píng)價(jià)

從研究區(qū)內(nèi)各行政區(qū)實(shí)際出發(fā)選擇合適的評(píng)價(jià)因子，并對(duì)其進(jìn)行相應(yīng)的因子分級(jí)，為保證所選取的評(píng)價(jià)因子之間相互獨(dú)立，采用ArcGIS的空間分析功能，對(duì)選取的因子進(jìn)行獨(dú)立性檢驗(yàn)，通過(guò)計(jì)算相關(guān)系數(shù)R來(lái)衡量各指標(biāo)之間的相關(guān)程度，R的取值范圍為[-1，1]，|R|越接近1，則表示兩兩指標(biāo)相關(guān)性越高，相關(guān)程度分級(jí)如表2所示。根據(jù)相關(guān)性計(jì)算結(jié)果，確定評(píng)價(jià)因子，采用信息量法計(jì)算信息量值，運(yùn)用自然斷點(diǎn)法結(jié)合研究區(qū)坡面孤石調(diào)查成果，劃分坡面孤石發(fā)育程度分區(qū)，研究區(qū)可劃分為坡面孤石強(qiáng)發(fā)育區(qū)、中等發(fā)育區(qū)、弱發(fā)育區(qū)和不發(fā)育區(qū)。

表2 |R|表征相關(guān)程度表

通過(guò)ArcGIS軟件提取各評(píng)價(jià)因子圖層的屬性數(shù)據(jù)，運(yùn)用空間分析工具計(jì)算得到各評(píng)價(jià)因子之間的相關(guān)性矩陣如表3所示。

表3 評(píng)價(jià)因子相關(guān)性矩陣

根據(jù)表2所示的|R|表征相關(guān)程度表，結(jié)合表3計(jì)算得到的各評(píng)價(jià)因子相關(guān)性矩陣表，本次選取的7個(gè)評(píng)價(jià)因子兩兩相關(guān)性系數(shù)絕對(duì)值均小于0.4，即任意兩評(píng)價(jià)因子之間均屬低度相關(guān)性，因而，可將初步選取的7個(gè)因子全部作為本次坡面孤石發(fā)育程度的評(píng)價(jià)因子。

根據(jù)研究區(qū)坡面孤石調(diào)查統(tǒng)計(jì)結(jié)果，按照信息量法計(jì)算坡面孤石發(fā)育程度各評(píng)價(jià)因子的信息量，各評(píng)價(jià)因子對(duì)應(yīng)的信息量計(jì)算結(jié)果如表4所示。

表4 各級(jí)評(píng)價(jià)因子信息量統(tǒng)計(jì)表

從各評(píng)價(jià)因子的信息量統(tǒng)計(jì)結(jié)果可以看出各評(píng)價(jià)因子在孤石易發(fā)性程度評(píng)價(jià)中的權(quán)重大小。①坡度因子信息量統(tǒng)計(jì)結(jié)果顯示，隨著坡度的增加，坡面孤石發(fā)育程度信息量先增大后減小，僅在15°～35°區(qū)間內(nèi)表現(xiàn)為正值，表明坡度在此區(qū)間易于坡面孤石形成和存在于坡面，其他坡度均不利于形成坡面孤石；②從坡向因子信息量統(tǒng)計(jì)結(jié)果可以看出，平地調(diào)查未發(fā)現(xiàn)孤石，表明平地不利于坡面孤石的形成，其余8個(gè)坡向?qū)?yīng)的信息量均在0左右浮動(dòng)，僅南東向、南向和南西向?yàn)檎?，說(shuō)明南向坡有利于坡面孤石發(fā)育；③從高程因子的信息量統(tǒng)計(jì)結(jié)果可以看出，高程小于 50 m和大于 150 m時(shí)信息量均為負(fù)，介于 50 m～150 m時(shí)信息量為正，說(shuō)明高程在 50 m～150 m之間有利于坡面孤石的形成；④曲率因子的信息量值浮動(dòng)較小，且曲率絕對(duì)值過(guò)大或過(guò)小時(shí)，信息量值偏小，說(shuō)明曲率對(duì)坡面孤石的形成影響不大，且直線形坡和大曲率坡均不利于坡面孤石的形成；⑤NDVI表示植被發(fā)育程度，其值越大表示植被越發(fā)育，建成區(qū)以外植被裸露區(qū)域極小，調(diào)查未發(fā)現(xiàn)坡面孤石，故信息量很小，隨著NDVI的增大，植被因子信息量先增大后減小，且僅在0.2～0.4之間時(shí)信息量為正，說(shuō)明植被較茂盛時(shí)最有利于坡面孤石形成，植被裸露和植被茂盛區(qū)均不利于坡面孤石的形成；⑥斷裂距離因子的信息量整體較小，說(shuō)明其對(duì)坡面孤石形成的影響較小，但從其變化趨勢(shì)來(lái)看，隨著斷裂距離的增加，信息量逐漸增大，說(shuō)明斷裂對(duì)基巖的作用不利于坡面孤石的形成；⑦巖性因子的信息量浮動(dòng)較大，說(shuō)明基巖巖性是坡面孤石形成的決定性因素，調(diào)查僅在晚三疊和中侏羅發(fā)現(xiàn)坡面孤石，說(shuō)明該時(shí)代花崗巖最有利于坡面孤石的形成，其他成因類型和形成時(shí)代的基巖均不利于形成坡面孤石。

6 發(fā)育程度分區(qū)

基于各評(píng)價(jià)因子信息量統(tǒng)計(jì)結(jié)果，利用ArcGIS的空間分析功能將7個(gè)評(píng)價(jià)因子圖層的數(shù)據(jù)進(jìn)行空間疊加分析，得到每個(gè)評(píng)價(jià)單元的總信息量，用該信息量值表征各單元格的坡面孤石發(fā)育程度指數(shù)，信息量頻率分布如圖2所示，運(yùn)用自然斷點(diǎn)法并結(jié)合研究區(qū)坡面孤石發(fā)育特征及分布規(guī)律，將研究區(qū)劃分為坡面孤石強(qiáng)發(fā)育區(qū)、中等發(fā)育區(qū)、低發(fā)育區(qū)和不發(fā)育區(qū)(表5)，并生成研究區(qū)坡面孤石發(fā)育程度分區(qū)圖。

圖2 坡面孤石信息量頻率分布圖

表5 坡面孤石發(fā)育程度劃分表

根據(jù)表5所示坡面孤石發(fā)育程度劃分標(biāo)準(zhǔn)，可將研究區(qū)劃分為坡面孤石強(qiáng)發(fā)育區(qū)、中等發(fā)育區(qū)、弱發(fā)育區(qū)和不發(fā)育區(qū)，其中，強(qiáng)發(fā)育區(qū)包括兩個(gè)亞區(qū)(A1、A2)，地貌類型以丘陵為主，面積共 58.40 km2，主要分布在聯(lián)和街西部天鹿湖森林公園一帶，長(zhǎng)嶺、聯(lián)和、蘿崗、永和交界處暹崗山-耙田山-香雪公園一帶，大沙、文沖、云埔交界的大田山一帶，紅山、穗東、云埔交界的龍頭山森林公園-丹水坑風(fēng)景區(qū)一帶；中等發(fā)育區(qū)包括兩個(gè)亞區(qū)(B1、B2)，地貌類型以丘陵和山間谷地為主，面積 100.07 km2，主要分布在廣汕公路-永順大道以北的聯(lián)和東部、長(zhǎng)嶺北部、永和北部、新龍南部丘陵區(qū)，聯(lián)和南部、大沙北部猛虎山-將軍地一帶、大沙南部丘陵區(qū)；弱發(fā)育區(qū)包括四個(gè)亞區(qū)(C1、C2、C3和C4)，地貌類型包括丘陵、臺(tái)地、山間谷地和三角洲平原，面積共 113.58 km2，主要分布在新龍北部、九佛和龍湖街丘陵、臺(tái)地區(qū)，聯(lián)和街西北部天鹿湖森林公園一帶，南部云埔街、南崗街、穗東街、紅山街、黃埔街、魚珠街和長(zhǎng)洲街的丘陵和剝蝕殘丘也有零星分布；不發(fā)育區(qū)包括兩個(gè)亞區(qū)(D1和D2)，地貌類型以丘陵、山間谷地和三角洲平原為主，面積共 212.12 km2，主要分布在黃埔區(qū)北部九佛街、龍湖街和新龍鎮(zhèn)的建成區(qū)，南部建成區(qū)和三角洲平原區(qū)。

將研究區(qū)坡面孤石分布情況與坡面孤石發(fā)育程度分區(qū)圖進(jìn)行疊加(圖3)可知，坡面孤石發(fā)育程度與實(shí)際調(diào)查中坡面孤石的分布具有高度一致性，充分證明了坡面孤石發(fā)育程度評(píng)價(jià)的準(zhǔn)確性。

圖3 坡面孤石發(fā)育程度分區(qū)圖

7 結(jié) 論

(1)坡面孤石是花崗巖球狀風(fēng)化產(chǎn)物，主要分布于以花崗巖為基巖的山地丘陵區(qū)，研究區(qū)坡面孤石以中小微型、淺埋、微風(fēng)化為主要發(fā)育特征，分布特征主要表現(xiàn)為集群分布、線狀分布和點(diǎn)狀分布。

(2)基于ArcGIS工具和信息量模型，選取7個(gè)影響坡面孤石發(fā)育的主要因子，對(duì)黃埔區(qū)進(jìn)行坡面孤石發(fā)育程度評(píng)價(jià)。計(jì)算各因子對(duì)坡面孤石發(fā)育的信息量，通過(guò)各因子信息量的空間疊加，結(jié)合坡孤石在研究區(qū)的分布特征，將研究區(qū)劃分為坡面孤石強(qiáng)發(fā)育區(qū)、中等發(fā)育區(qū)、弱發(fā)育區(qū)和不發(fā)育區(qū)4類。將坡面孤石發(fā)育程度分區(qū)與調(diào)查樣本分布進(jìn)行對(duì)比分析，結(jié)果表明，坡面孤石發(fā)育程度分區(qū)具有較高的準(zhǔn)確性。

(3)通過(guò)各因子信息量計(jì)算結(jié)果可知，對(duì)坡面孤石發(fā)育程度影響較大的因素為基巖巖性、地形坡度、坡向和高程，晚三疊和中侏羅的花崗巖有利于坡面孤石的形成，50 m～150 m高程范圍內(nèi)坡度介于15°～35°的南向坡有利于坡面孤石的發(fā)育。

(4)隨著坡面孤石發(fā)育程度等級(jí)的提升，面積占比逐漸減小，研究區(qū)坡面孤石強(qiáng)發(fā)育區(qū)面積共 58.40 km2，占研究區(qū)總面積的12.00%；中等發(fā)育區(qū)面積共 100.07 km2，占比為20.56%；弱發(fā)育區(qū)面積共 113.58 km2，占比為23.34%；不發(fā)育區(qū)面積共 214.58 km2，占比為44.10%。