馬天華，廖祿云，陳國輝，劉泊雷，石元帥

基于三維傾斜攝影及GIS技術(shù)的華鎣山康養(yǎng)城地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)

馬天華1，2，廖祿云1，2，陳國輝2，劉泊雷2，石元帥2

（1.四川省天晟源環(huán)保股份有限公司，成都 610037；2.四川省地質(zhì)工程勘察院集團(tuán)有限公司，成都 610072）

本文在分析三維傾斜攝影成果及研究區(qū)地質(zhì)環(huán)境條件的基礎(chǔ)上，利用GIS技術(shù)，選取了高程、坡向、坡度、曲率、地層巖性、與采空區(qū)的距離、與公路的距離、地質(zhì)災(zāi)害面密度等8個(gè)評(píng)價(jià)指標(biāo)，采用信息量法對(duì)研究區(qū)進(jìn)行地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)，將研究區(qū)劃分為低危險(xiǎn)區(qū)、中危險(xiǎn)區(qū)、高危險(xiǎn)區(qū)，并對(duì)各危險(xiǎn)區(qū)進(jìn)行評(píng)價(jià)，在此基礎(chǔ)上，提出了防治措施建議。結(jié)論可作為華鎣山康養(yǎng)城規(guī)劃建設(shè)的地質(zhì)基礎(chǔ)支撐。

華鎣山康養(yǎng)城；傾斜攝影；GIS；信息量法；危險(xiǎn)性評(píng)價(jià)

為貫徹落實(shí)習(xí)近平總書記關(guān)于“推動(dòng)成渝地區(qū)雙城經(jīng)濟(jì)圈建設(shè)”的部署，四川省廣安市擬加快推進(jìn)華鎣山康養(yǎng)城的建設(shè)。華鎣山康養(yǎng)城定位是打造資源枯竭型城市轉(zhuǎn)型的發(fā)展典范、高品位深度休閑體驗(yàn)的旅游典范、傳承歷史文脈的人文典范、構(gòu)建投資樂土的宜業(yè)典范，集文化朝圣體驗(yàn)、避暑養(yǎng)生度假、生態(tài)觀光探奇、娛樂休閑運(yùn)動(dòng)等功能為一體的國際立體山地休閑度假旅游目的地①基于三維傾斜攝影技術(shù)的華鎣山康養(yǎng)城地質(zhì)環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)報(bào)告. 四川省地質(zhì)工程勘察院集團(tuán)有限公司，2019.。

然而，華鎣山康養(yǎng)城規(guī)劃區(qū)由于歷史采礦等一系列人類工程活動(dòng)，引發(fā)了較多的地質(zhì)環(huán)境問題，對(duì)華鎣山康養(yǎng)城規(guī)劃建設(shè)帶來一定的影響。三維傾斜攝影技術(shù)能為研究區(qū)提供較為全面直觀的地形數(shù)據(jù)、地質(zhì)災(zāi)害發(fā)育情況等地質(zhì)環(huán)境要素，結(jié)合GIS技術(shù)對(duì)研究區(qū)進(jìn)行地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)，為華鎣山康養(yǎng)城規(guī)劃建設(shè)提供地質(zhì)基礎(chǔ)支撐。

圖1 控制點(diǎn)布設(shè)及研究范圍圖

1 研究區(qū)概況

華鎣山康養(yǎng)城位于廣安市華鎣市及鄰水縣中部，地處亞熱帶濕潤季風(fēng)氣候區(qū)，四季分明，多年平均氣溫17.2℃，多年平均降雨量1082.2mm。地貌類型為低山—低中山地貌，華鎣山中段的西翼大部分是低山，只在海拔超過1000m的華鎣山背脊一帶，分布有低中山地貌。地層巖性主要有二疊系、三疊系及第四系之黏性土、碎塊石土、泥巖、砂巖、頁巖、灰?guī)r等，巖性巖相變化較大。研究區(qū)處于華鎣山復(fù)式背斜中段，地質(zhì)構(gòu)造較為復(fù)雜?；?、崩塌、地面塌陷等地質(zhì)災(zāi)害較為發(fā)育，地質(zhì)環(huán)境問題較為突出②華鎣山景區(qū)李子埡村—704轉(zhuǎn)播站片區(qū)建設(shè)用地地質(zhì)環(huán)境危險(xiǎn)性評(píng)估報(bào)告. 四川煤田一四一建設(shè)投資有限公司，2018.③華鎣山景區(qū)李子—椿木坪片區(qū)建設(shè)用地地質(zhì)環(huán)境危險(xiǎn)性評(píng)估報(bào)告. 四川煤田一四一建設(shè)投資有限公司，2013.。

2 數(shù)據(jù)采集及處理

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

數(shù)據(jù)采集采用電動(dòng)垂降固定翼無人機(jī)CW-15，搭載五鏡頭傾斜云臺(tái)，布設(shè)12個(gè)像控點(diǎn)、8個(gè)檢查點(diǎn)及2個(gè)水泥樁控制點(diǎn)（圖1）。正射影像時(shí)，旁向重疊度40%，航向重疊度70%；傾斜影像時(shí)，旁向重疊度75%，航向重疊度85%。完成了38km2（其中，研究區(qū)27.78km2）的航飛影像分辨率0.1m的傾斜攝影測(cè)量數(shù)據(jù)采集，形成了正射影像（DOM）、高程模型（DEM）以及全景三維模型等成果（表1）。

2.2 數(shù)據(jù)提取及處理

（1）根據(jù)傾斜攝影的成果，解譯出地形地貌、地層巖性、水系、人類工程活動(dòng)、地質(zhì)災(zāi)害（滑坡、不穩(wěn)定斜坡、崩塌、地面塌陷）等地質(zhì)環(huán)境條件及地質(zhì)環(huán)境問題。

表1 研究區(qū)傾斜攝影工作量

（2）基于傾斜攝影技術(shù)的基礎(chǔ)上，利用GIS技術(shù)，提取高程、坡向、坡度、曲率、地層巖性、與采空區(qū)的距離、與公路的距離、地質(zhì)災(zāi)害面密度等數(shù)據(jù)。

3 地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)

3.1 評(píng)價(jià)方法

本次評(píng)價(jià)采用信息量法，信息量法的理論基礎(chǔ)是信息論，用于地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)的主要思路是根據(jù)對(duì)已知災(zāi)害的現(xiàn)實(shí)情況提供的信息，把影響因素的實(shí)測(cè)值轉(zhuǎn)化為反映地質(zhì)環(huán)境的信息量值。地質(zhì)現(xiàn)象()受多種因素(x,=1,2,的影響，各種因素所起作用的大小、性質(zhì)是不相同的。因此，對(duì)于區(qū)域地質(zhì)問題要綜合研究其影響因素類別和具體狀態(tài)的組合，而不是停留在單個(gè)因素上。地質(zhì)產(chǎn)生與否是與預(yù)測(cè)過程中所獲取的信息的數(shù)量和質(zhì)量有關(guān)的，用信息量來表示即為：

圖3 危險(xiǎn)性評(píng)價(jià)體系圖

根據(jù)條件概率運(yùn)算，上式可進(jìn)一步寫成：

I(y,xx…x)=I(y,x)+Ix(y,x)+…+Iy,xx…x(y,x)（2）

式中：I(y,xx…x)—因素組合xx…x對(duì)地災(zāi)所提供的信息量；P(y,xx…x)—因素xx…x組合條件下滑坡發(fā)生的概率；I(y,x)—因素x1存在時(shí)，因素x對(duì)地災(zāi)提供的信息量；—地災(zāi)發(fā)生的概率。

公式（2）說明：因素組合xx…x對(duì)地質(zhì)環(huán)境所提供的信息量等于因素x1提供的信息量，加上因素x1確定后因素x對(duì)地質(zhì)環(huán)境提供的信息量，直至因素xx…x確定后，x對(duì)地質(zhì)環(huán)境提供的信息量。從而說明區(qū)域地質(zhì)信息預(yù)測(cè)是充分考慮因素組合的共同影響與作用。

3.2 評(píng)價(jià)指標(biāo)體系的建立

根據(jù)華鎣山康養(yǎng)城地質(zhì)環(huán)境條件及人類工程活動(dòng)的影響，遵循評(píng)價(jià)因子簡明性、代表性和可操作性的原則，并考慮現(xiàn)有數(shù)據(jù)與GIS技術(shù)功能的局限性等條件，最終選取了地形地貌、地質(zhì)環(huán)境、人類工程活動(dòng)、災(zāi)害特征作為評(píng)價(jià)指標(biāo)（白玉明，2019；袁茂珂等，2021；陳國輝等，2019；賀鵬等，2016；詹慶明等，2016；朱良峰等，2002）（圖2）。

圖2 重分類對(duì)話框界面

3.3 危險(xiǎn)性分析與計(jì)算

（1）各評(píng)價(jià)指標(biāo)的重分類

在實(shí)際情況中，某些因素對(duì)地質(zhì)環(huán)境的影響并不完全是按要素量的增加而成正比增加的，而往往在一個(gè)數(shù)量范圍內(nèi)對(duì)地質(zhì)環(huán)境的影響是穩(wěn)定的，如距離道路2m和距離道路5m的地方其對(duì)地質(zhì)環(huán)境的影響可能幾乎是一樣的，而150m則可能與5m的地方有較大的差異。因此，需要對(duì)部分連續(xù)數(shù)據(jù)重新分類。首先Spatial Analyst Tools中Reclass模塊中的Reclassify工具，將各因子?xùn)鸥駡D層進(jìn)行重分類，隨后計(jì)算出i等級(jí)所占的面積與研究區(qū)總面積比值，及分級(jí)面積比。以高程分級(jí)面積比計(jì)算為例，根據(jù)研究區(qū)地形條件高程分級(jí)為：0～200、200～400、400～600、600～800、800～1000、1000～1200、>1200；分類主要在ArcGIS空間分析模塊中利用reclassify命令完成（圖3）。

表2 高程評(píng)價(jià)因子重分類分級(jí)

表3 坡度評(píng)價(jià)因子重分類分級(jí)

表4 坡向評(píng)價(jià)因子重分類分級(jí)

參照區(qū)域地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)分級(jí)賦值形式，并依據(jù)誘發(fā)滑坡、崩塌、地面塌陷、不穩(wěn)定斜坡等地質(zhì)環(huán)境問題的可能性，根據(jù)《礦山地質(zhì)環(huán)境保護(hù)與恢復(fù)治理方案編制規(guī)范》附錄 E、2003年全國山洪災(zāi)害防治規(guī)劃組定制的《山洪災(zāi)害防治規(guī)劃地形地質(zhì)區(qū)劃技術(shù)細(xì)則》等技術(shù)標(biāo)準(zhǔn)對(duì)各個(gè)評(píng)價(jià)因子分級(jí)，各評(píng)價(jià)因子重分類情況見表2～9。

表5 曲率評(píng)價(jià)因子重分類分級(jí)

表6 地層巖性評(píng)價(jià)因子重分類分級(jí)

表7 與公路的距離評(píng)價(jià)因子重分類分級(jí)

表8 與采空區(qū)的距離評(píng)價(jià)因子重分類分級(jí)

表9 地質(zhì)災(zāi)害面密度評(píng)價(jià)因子重分類分級(jí)

（2）信息量值計(jì)算

將各圖層因子的屬性表和對(duì)應(yīng)的與地質(zhì)災(zāi)害分布作“乘”運(yùn)算后得到的各圖層因子的屬性表存為Excel表格，在Excel中按如下公式計(jì)算出每個(gè)圖層各類的信息量值：

式中：為研究區(qū)像素總數(shù)；為研究區(qū)內(nèi)含有滑坡分布的像素總數(shù)；S為研究區(qū)內(nèi)含有評(píng)價(jià)因素x的像素?cái)?shù)；N為分布在因素x內(nèi)特定類別內(nèi)的災(zāi)害點(diǎn)像素?cái)?shù)。

最終得到8個(gè)圖層因子內(nèi)各類別的信息量值（表10～17）。然后按照各類別與對(duì)應(yīng)的信息量之間的關(guān)系，對(duì)各類別賦其對(duì)應(yīng)的信息量值，得到8個(gè)信息量圖層。

表10 高程和災(zāi)害點(diǎn)的關(guān)系及其對(duì)應(yīng)的信息量值

表11 坡度和災(zāi)害點(diǎn)的關(guān)系及其對(duì)應(yīng)的信息量值

表12 坡向和災(zāi)害點(diǎn)的關(guān)系及其對(duì)應(yīng)的信息量值

表13 曲率和災(zāi)害點(diǎn)的關(guān)系及其對(duì)應(yīng)的信息量值

表14 地層巖性和災(zāi)害點(diǎn)的關(guān)系及其對(duì)應(yīng)的信息量值

表15 與采空區(qū)的距離、災(zāi)點(diǎn)關(guān)系及對(duì)應(yīng)的信息量值

表16 與公路的距離、災(zāi)點(diǎn)關(guān)系及對(duì)應(yīng)的信息量值

表17 災(zāi)害點(diǎn)和災(zāi)害點(diǎn)的關(guān)系及其對(duì)應(yīng)的信息量值

根據(jù)高程、坡度、坡向等8個(gè)指標(biāo)的信息量值在-3.16～21.5之間，信息量值與地災(zāi)害危險(xiǎn)性呈正比，信息量值越大，表明該區(qū)域發(fā)生地災(zāi)的概率越大，危險(xiǎn)性等級(jí)也越高。

3.4 危險(xiǎn)性評(píng)價(jià)結(jié)果及區(qū)劃

完成所有評(píng)論因子信息量值柵格圖層以后，運(yùn)用Spatial Analyst Tools中的Map Algebra模塊的RASTER Calculatr 工具，將評(píng)價(jià)因子信息值疊加，得到研究區(qū)總信息量值。并運(yùn)用Reclassify工具對(duì)其進(jìn)行重分類為低危險(xiǎn)區(qū)、中危險(xiǎn)區(qū)、高危險(xiǎn)區(qū)三個(gè)區(qū)。

表18 研究區(qū)地質(zhì)環(huán)境危險(xiǎn)性程度分區(qū)特征表

高危險(xiǎn)區(qū)以滑坡、地面塌陷較為發(fā)育，不穩(wěn)定斜坡、崩塌（危巖）次之，發(fā)育密度為1.52個(gè)/km2；地形坡度一般在35°～45°之間；巖性多為軟硬相間的紫紅色砂巖、礫巖與泥巖、頁巖不等厚互層，且?guī)r層傾角多小于地形坡度；褶皺、斷裂發(fā)育；在暴雨或連續(xù)降水條件下易引發(fā)地質(zhì)災(zāi)害發(fā)生。按重點(diǎn)防治區(qū)進(jìn)行地質(zhì)災(zāi)害風(fēng)險(xiǎn)管控，不應(yīng)進(jìn)行人工切坡，防止人類工程活動(dòng)引發(fā)災(zāi)害發(fā)生，盡量避開地質(zhì)災(zāi)害危險(xiǎn)區(qū)域或地段，無法避免時(shí)應(yīng)進(jìn)行工程治理。

中危險(xiǎn)區(qū)主要發(fā)育不穩(wěn)定斜坡為主，次為崩塌（危巖），發(fā)育密度為0.59個(gè)/km2；地形坡度一般在20°～30°之間；巖性以灰?guī)r、白云質(zhì)灰?guī)r、鈣質(zhì)泥巖為主，且?guī)r層傾角多與地形坡度接近；褶皺、斷裂較為發(fā)育；在暴雨或連續(xù)降水條件下較易引發(fā)地質(zhì)災(zāi)害發(fā)生。按次重點(diǎn)防治區(qū)進(jìn)行地質(zhì)災(zāi)害風(fēng)險(xiǎn)管控，應(yīng)盡量避免人工切坡，對(duì)現(xiàn)有隱患點(diǎn)采取合適的防治措施，加強(qiáng)邊坡穩(wěn)定性監(jiān)測(cè)。

低危險(xiǎn)區(qū)主要發(fā)育滑坡、崩塌（危巖）為主，發(fā)育密度為0.28個(gè)/km2；地形坡度一般15°～30°；巖性以砂巖、泥巖為主且?guī)r層傾角多大于地形坡度；受褶皺、斷裂的影響較大；在暴雨或連續(xù)降水條件下可引發(fā)地質(zhì)災(zāi)害發(fā)生。按一般區(qū)進(jìn)行地質(zhì)災(zāi)害風(fēng)險(xiǎn)管控，加強(qiáng)群策群防。

圖４ 研究區(qū)地質(zhì)環(huán)境危險(xiǎn)性程度分區(qū)圖

3.5 防治措施建議

本評(píng)價(jià)結(jié)果可作為研究區(qū)規(guī)劃建設(shè)的地質(zhì)依據(jù)，根據(jù)前述評(píng)價(jià)結(jié)果，提出如下建議：①應(yīng)結(jié)合地質(zhì)環(huán)境條件進(jìn)一步優(yōu)化相關(guān)規(guī)劃內(nèi)容；②建設(shè)時(shí)應(yīng)加強(qiáng)地質(zhì)災(zāi)害危險(xiǎn)性評(píng)估；③必要時(shí)，應(yīng)對(duì)地質(zhì)災(zāi)害進(jìn)行治理；④建立健全群測(cè)群防體系，落實(shí)地質(zhì)災(zāi)害防治責(zé)任制。

4 結(jié)論

（1）滿足精度要求的傾斜攝影技術(shù)可以較高效率的情況下提供評(píng)價(jià)所需各項(xiàng)因子。

（2）在傾斜攝影成果的基礎(chǔ)上，采用GIS技術(shù)提取了高程、坡度、坡向、曲率、與采空區(qū)的距離、與公路的距離、地質(zhì)災(zāi)害面密度等8個(gè)評(píng)價(jià)指標(biāo)，并進(jìn)行了分析與計(jì)算，結(jié)果與實(shí)際較為符合。

（3）采用信息量法對(duì)研究區(qū)進(jìn)行地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)，評(píng)價(jià)結(jié)果與實(shí)際情況較為接近。

（4）本地質(zhì)環(huán)境危險(xiǎn)性評(píng)價(jià)結(jié)果可作為華鎣山康養(yǎng)城規(guī)劃建設(shè)的地質(zhì)基礎(chǔ)支撐。

白玉明. 2019. 桂林恭城巖口小學(xué)地質(zhì)災(zāi)害危險(xiǎn)性分析及防治[J]. 四川地質(zhì)學(xué)報(bào)，39(2)：308-314.

袁茂珂，宋偉，杜杰，鄭小敏，卓曼他. 2021. 基于3S技術(shù)的九寨溝8·8地震后生態(tài)地質(zhì)環(huán)境敏感性研究[J]. 四川地質(zhì)學(xué)報(bào)，41(2)：289-294.

陳國輝，吳森，潘國耀. 2019. 基于貢獻(xiàn)率模型的管道河(溝)道水毀災(zāi)害危險(xiǎn)性評(píng)價(jià)指標(biāo)體系[J]. 中國地質(zhì)災(zāi)害與防治，31(4)：128-123.

賀鵬，童立強(qiáng)，郭兆成，劉春玲，涂杰楠，王珊珊，許杰敏. 2016. GIS支持下基于層次分析法的西藏札達(dá)地區(qū)滑坡災(zāi)害易發(fā)性評(píng)價(jià)研究[J]. 科學(xué)技術(shù)與工程，16(25)：193-200

詹慶明，張維思，鄒芳，劉慧民，陳暢. 2016. 基于GIS的神農(nóng)架林區(qū)崩滑流災(zāi)害風(fēng)險(xiǎn)評(píng)價(jià)研究[J]. 地理信息世界，23(5)：1-9

朱良峰，殷坤龍，張梁，李閩. 2002. 基于GIS技術(shù)的地質(zhì)災(zāi)害風(fēng)險(xiǎn)分析系統(tǒng)研究[J]. 工程地質(zhì)學(xué)報(bào)，10(4)：428-433,348.

Geological Environment Risk Assessment of Kangyang City of Huaying Mountain Based on 3D Oblique Photography and GIS Technology

MA Tian-hua1,2LIAO Lu-yunn1,2CHEN Guo-hui2LIU Bo-lei2SHI Yuan-shuai2

(1-Sichuan Tianshengyuan Environmental Protection Co., LTD., Chengdu 610037; 2-Sichuan Geological Engineering Survey Institute Group Co., LTD., Chengdu 610072)

Based on the analysis of 3D oblique photography results and geological environment conditions of the study area, this paper uses GIS technology to select eight evaluation indexes, including elevation, slope direction, slope, curvature, stratum lithology, distance from goaf, distance from highway, and geological hazard surface density, and uses information volume method to evaluate geological environmental risk of the study area. The study area is divided into low risk area, medium risk area and high risk area, and each risk area is evaluated. On this basis, the prevention and control measures are proposed. It can be used as the geological base support for the planning and construction of Huaying Mountain health care city.

Huaying Mountain health city; oblique photography; GIS; information volume method; risk assessment

P237

A

1006-0995（2022）04-0695-06

10.3969/j.issn.1006-0995.2022.04.030

2022-11-04

馬天華（1976— ），男，重慶石柱人，本科，高級(jí)工程師，主要從事環(huán)境保護(hù)方面的工作.

陳國輝（1977— ），男，福建漳浦人，正高級(jí)工程師，主要從事水工環(huán)地質(zhì)方面的研究工作