遙感與GIS技術(shù)支持下的機(jī)場選址與工程地質(zhì)分析

李天華,楊武年

(1.成都空軍勘察設(shè)計(jì)院,成都 610041;2.成都理工大學(xué) 遙感與GIS研究所,成都 610059)

遙感與GIS技術(shù)支持下的機(jī)場選址與工程地質(zhì)分析

李天華1,楊武年2

(1.成都空軍勘察設(shè)計(jì)院,成都 610041;2.成都理工大學(xué) 遙感與GIS研究所,成都 610059)

研究將“3S”技術(shù)(RS、GIS、GPS)應(yīng)用到機(jī)場地址勘察中,主要是采用現(xiàn)代化遙感技術(shù),以最新遙感圖像LANDSAT ETM(最高空間分辨率15m),SPOT圖像(2.5m)和QU ICKB IRD圖像(0.61m)等作為信息源,通過遙感圖像數(shù)據(jù)融合、數(shù)字鑲嵌和影像增強(qiáng)等數(shù)字處理,制作不同類型、不同波段和不同分辨率的正射遙感影像及正射影像地圖,為機(jī)場工程地質(zhì)勘察提供最佳圖像資料。通過遙感圖像增強(qiáng)處理、信息提取及GIS空間分析等,對研究區(qū)的地質(zhì)構(gòu)造,地質(zhì)災(zāi)害等進(jìn)行解譯。并借助GPS進(jìn)行野外現(xiàn)場調(diào)研,查明研究區(qū)地質(zhì)構(gòu)造的特征,尤其是隱伏斷裂、活動性斷裂的分布,以及巖溶漏斗、洼地等水文地質(zhì)現(xiàn)象,為小哨機(jī)場工程地質(zhì)勘察工作提供了基礎(chǔ)資料。

“3S”技術(shù);地質(zhì)勘察;機(jī)場選址;正射影像

1 機(jī)場概況

昆明小哨地區(qū)位于云貴高原中部,處于云南省境內(nèi)自西北向東南遞降地勢的第二大階梯上。大部份地區(qū)海拔在1 500m～2 800m之間,最高點(diǎn)是北部祿勸縣拱王山脈的馬鬃嶺,海拔4 247.7m;最低點(diǎn)是祿勸縣則黑鄉(xiāng)小河坪子?xùn)|北1 km處普渡河與金沙江交匯點(diǎn),海拔746m。該區(qū)內(nèi)有二條主要山脈:①川西魯南山脈越過金沙江南下的拱王山脈,蜿蜒盤旋于祿勸、富民、西山、安寧、晉寧等地,最高峰馬鬃嶺,海拔4 247.7m,為區(qū)內(nèi)最高峰;②與滇東北烏蒙山脈連接的梁王山脈,分布于嵩明、官渡、呈貢、宜良、石林等縣區(qū),最高峰海拔2 833m。受主體構(gòu)造線方向的控制,二條山脈走向近南北向。研究區(qū)大部份地形呈波狀起伏,平均海拔2 000m左右,平均自然坡度5°～8°,表現(xiàn)為丘陵略有起伏溶蝕平原和沖洪積平原。

作者在本次研究的方案與技術(shù)思路,如圖1(見下頁)所示。

2 遙感圖像數(shù)字處理

考慮要進(jìn)行機(jī)場選址設(shè)計(jì)和工程勘察,作者選用了QuickB ird圖像、SPOT圖像、TM圖像,作為主要原始圖像數(shù)據(jù)。

通過對原始遙感圖像進(jìn)行多波段數(shù)據(jù)融合,圖像增強(qiáng),正射精校正,高精度DEM生成等處理,并調(diào)用DEM中的高程信息和遙感圖像展現(xiàn)的形狀、色彩、質(zhì)地等信息,將二種信息融合到一起,形成對研究區(qū)地貌的三維可視化。從而提供可供反復(fù)使用的真實(shí)、客觀、信息連續(xù)的宏觀分析地面景觀影像,使機(jī)場建設(shè)工作者在一種接近現(xiàn)實(shí)的虛擬環(huán)境中,進(jìn)行機(jī)場選址設(shè)計(jì),場區(qū)的構(gòu)造地質(zhì)、工程地質(zhì)、環(huán)境地質(zhì)等特征信息提取和綜合評價(jià)。

2.1 遙感圖像正射精校正

由于正射校正可以消除投影差的影響,因此是衛(wèi)星區(qū)域鑲嵌圖和標(biāo)準(zhǔn)地理分幅遙感圖像制作的基礎(chǔ)。正射精校正的原理和幾何精校正差不多,但正射精校正有其自身的特點(diǎn):①用于正射精校正的地面控制點(diǎn),必須具有高程信息;②進(jìn)行正射精校正時(shí),必須利用數(shù)字高程模型(DEM)。遙感影像正射校正是遙感正射影像地圖的重要信息源,同時(shí)它也是進(jìn)行其它遙感圖像處理(遙感圖像融合、鑲嵌等)的基礎(chǔ)。

圖1 研究方案與技術(shù)思路Fig.1 Research p rogram s and technical ideas

正射精校正的精度與DEM、GCP或DRG的精度有關(guān),還與衛(wèi)星影像自身的分辨率有關(guān)。作者對研究區(qū)進(jìn)行了正射精校正,針對不同的圖像采用不同的正射校正方法。對于Spot5圖像,利用的Erdas中Spo t5模型對圖像進(jìn)行了正射校正。對于QuickB ird圖像,作者采用了RPC(有理多項(xiàng)式)校正方法實(shí)現(xiàn)對圖像的正射校正。正射校正的精度一般控制在半個(gè)像元內(nèi),其地面控制點(diǎn)(GCP)的高程值通過DEM自動提取。校正好的正射影像圖與等高線套合,效果較好,可滿足1∶5 000區(qū)調(diào)填圖和遙感解譯使用。

2.2 遙感圖像三維可視化實(shí)現(xiàn)

遙感圖像三維可視化實(shí)現(xiàn)的技術(shù)和思想,就是依據(jù)DEM建立表面模型來顯示真實(shí)地形,然后再疊加遙感影像或航片等數(shù)字影像進(jìn)行紋理疊加,來顯示地表細(xì)節(jié)。充分發(fā)揮計(jì)算機(jī)圖示技術(shù)及虛擬現(xiàn)實(shí)技術(shù)的優(yōu)勢,利用遙感圖像、地理要素和文字符號標(biāo)注等多種數(shù)據(jù),生成三維地形影像[5]。將這些可視化數(shù)據(jù)集成套合成三維地形影像,必須做到不同數(shù)據(jù)間的坐標(biāo)配準(zhǔn),目的是將不同來源的同一地區(qū)的圖像、地理要素,以及文字符號等數(shù)據(jù),轉(zhuǎn)換到同一坐標(biāo)系中[6]。作者在本文中是以地形圖的地理坐標(biāo)作為配準(zhǔn)參考,其它數(shù)據(jù)均需轉(zhuǎn)換到地形圖所在的坐標(biāo)系中。

遙感圖像地質(zhì)解譯三維可視化,與影像動態(tài)分析系列動畫制作工序主要包括以下內(nèi)容:

(1)遙感圖像數(shù)字處理。遙感圖像數(shù)字處理是數(shù)字區(qū)調(diào)工作中信息提取的重要環(huán)節(jié)。圖像處理的目的,是對原始遙感圖像進(jìn)行輻射校正、幾何校正和投影差改正、地理編碼、圖像鑲嵌、圖像增強(qiáng),以及不同時(shí)相不同傳感器不同分辨率的多波段遙感圖像數(shù)據(jù)融合等,最終制作出統(tǒng)一規(guī)格標(biāo)準(zhǔn)的高質(zhì)量遙感圖像,提高數(shù)字區(qū)調(diào)地質(zhì)解譯應(yīng)用效果[2～4]。

多類型遙感影像數(shù)據(jù)融合,是一種通過高級影像處理技術(shù)來復(fù)合多源遙感影像數(shù)據(jù)的技術(shù),其目的是將單一傳感器的多波段信息或不同類傳感器所提供的信息加以綜合,消除多傳感器信息之間可能存在的冗余和矛盾,加以互補(bǔ),降低其不確定性,減少模糊度,以增強(qiáng)影像中信息透明度,提高圖像解像力和解譯精度[5、6]。

(2)高精度DEM生成。高精度三維立體圖像需要高精度的DEM來支撐,在遙感圖像的正射處理中,各像點(diǎn)的投影差改正也需要對應(yīng)點(diǎn)的高精度DEM[4]。高精度DEM生成常用有三種方法:①利用地形圖數(shù)字化求得;②利用立體像對(具有立體效果的航片、衛(wèi)星圖像等)使用專門的軟件求得;③利用干涉雷達(dá)圖像處理技術(shù)求得。

(3)三維飛行路線選取。根據(jù)工作區(qū)的地質(zhì)構(gòu)造復(fù)雜程度和工作需要,按照一定的規(guī)則進(jìn)行飛行勘察路線部署。

(4)三維可視化系列動畫產(chǎn)品制作。首先進(jìn)行三維飛行的參數(shù)設(shè)置,如航高、時(shí)速、夸大系數(shù)、屏幕大小、視角設(shè)置及背景效果等;然后根據(jù)布置的飛行路線完成三維動畫制作。

(5)三維可視化產(chǎn)品輸出。根據(jù)選擇的飛行路線,逐條生成遙感圖像地質(zhì)解譯三維可視化及影像動態(tài)分析系列動畫,并把這些產(chǎn)品轉(zhuǎn)化成通用動畫所支持的格式打包,如QuickTim eM ovie格式等,刻盤提交給有關(guān)的工作人員使用。

(6)信息提取。地質(zhì)構(gòu)造、巖溶地貌解譯和影像判讀等。

作者在本次工作中共制作了四條飛行路線,分別位于測區(qū)的西部、中部和東部,其目的是為了了解測區(qū)的總體區(qū)域地質(zhì)構(gòu)造格局以及巖溶發(fā)育、地貌特征等情況(見圖2)。

3 基于“3S”技術(shù)機(jī)場工程地質(zhì)分區(qū)

在上述所有工作基礎(chǔ)上,綜合運(yùn)用遙感與GIS對遙感圖像、機(jī)場數(shù)字DEM、三維可視化(見圖2),以及前期工作中提取的地貌、地層、構(gòu)造等信息進(jìn)行分析,將小哨機(jī)場地區(qū)試驗(yàn)方案劃分為四個(gè),并以方案三區(qū)為例,對該區(qū)工程地質(zhì)條件進(jìn)行了評價(jià)。

3.1 場地工程地質(zhì)分析

方案三出露泥盆系宰格組白云巖、?？诮M上段灰?guī)r及下段砂巖、雙龍?zhí)督M白云巖、陡坡寺組砂巖及泥巖,未見斷層。但裂隙發(fā)育,碳酸鹽巖和碎屑巖相間分布,成為走向北東傾南東的單斜狀構(gòu)造,極有利于巖溶的發(fā)育。單個(gè)漏斗的長軸方向受結(jié)構(gòu)面控制,以近北東向、東西向及北西向?yàn)橹?但漏斗的空間分布很大程度上具有層控性(巖性段控制)。地質(zhì)構(gòu)造對工程的影響主要是:①傾斜互層狀碳酸鹽巖與砂泥巖單斜構(gòu)造導(dǎo)致區(qū)內(nèi)巖溶發(fā)育;②產(chǎn)狀平緩,次級褶皺發(fā)育的宰格組內(nèi)巖溶受橫張節(jié)理控制,漏斗星羅棋布,地下巖溶管道發(fā)育。

3.2 工程地質(zhì)分區(qū)

綜合“3S”技術(shù)、綜合物探、區(qū)域地質(zhì)、工程地質(zhì)勘察等綜合分析,擬建場地的工程地質(zhì)巖組可分以下幾個(gè)亞區(qū)(見下頁圖3):

(1)構(gòu)造剝蝕丘陵。地貌為深丘,位于試驗(yàn)區(qū)北緣,面積0.096 km2,海拔2 120.4 m～2 137m,坡度3°～10°,為非可溶性砂泥巖分布區(qū),僅出現(xiàn)二處洼地,單個(gè)面積均達(dá)到4 900 m2以上。地形波狀起伏,殘積層厚,為工程挖方區(qū),巖石單軸抗壓強(qiáng)度可以達(dá)到12M Pa～15M Pa。

圖2 昆明小哨機(jī)場遙感圖像三維可視化系列圖片之一幀(右上角箭頭指向南)Fig.2 One of the remo te sensing im agesof three-dim ensional visualization of X iao-shao airport in Kunm ing(the arrow points to South)

(2)剝蝕溶蝕丘陵。地形地貌為深丘,面積約0.027 km2,海拔2 095m～2 135m,坡度3°～10°,巖溶地貌漏斗和洼地發(fā)育,局部漏斗和洼地有落水洞。地層巖性為泥盆系中統(tǒng)?？诮M(D 2h2),由灰黃色中、厚層狀亮晶砂屑灰?guī)r、礫屑灰?guī)r夾瀝青質(zhì)灰?guī)r頂部為砂巖、雙龍?zhí)督M白云巖組成。上覆紅色粘土。巖體節(jié)理裂隙發(fā)育,巖石破碎,巖溶漏斗發(fā)育。巖質(zhì)堅(jiān)硬,屬堅(jiān)硬巖,巖石單軸抗壓強(qiáng)度可達(dá)到50M Pa～90M Pa,是良好的基礎(chǔ)持力層。

③構(gòu)造剝蝕丘陵。地形地貌為深丘,面積約0.14 km2,海拔2 095m～2 135m,坡度3°～10°,為?？诮M下段砂巖分布區(qū),上覆厚度較薄的殘積。該區(qū)域出現(xiàn)八處洼地,單個(gè)面積在320 m2～1 400m2之間,少部份達(dá)到5 100m2以上。巖石單軸抗壓強(qiáng)度高,開挖后為良好的持力層。

④剝蝕溶蝕殘丘。地形地貌為淺丘,面積約0.377 km2,海拔2 075m～2 121m,坡度3°～8°,巖溶地貌漏斗和洼地特別發(fā)育,共出現(xiàn)洼地四十五處,面積大致在200m2～3 000 m2之間,其中以Kh323最大,面積達(dá)到21 418m2。

圖3 方案三工程地質(zhì)分區(qū)示意圖Fig.3 Schem atic diagram of 3 engineering geo logy division diagram

4 結(jié)論

作者在本文中將“3S”技術(shù)應(yīng)用到西部高原區(qū)機(jī)場選址和工程勘察中,虛擬了擬選機(jī)場場址的高精度三維遙感動態(tài)模型,從圖像的三維可視化和動態(tài)分析角度,進(jìn)行了擬選機(jī)場場址的勘察和工程地質(zhì)概況等遙感信息的提取。這改變了傳統(tǒng)的機(jī)場勘探工程中的思維方式,為機(jī)場勘察提供了新的技術(shù)平臺,為機(jī)場建設(shè)工程開拓了新思路,為決策者(各級領(lǐng)導(dǎo)、工程負(fù)責(zé)人等)從宏觀把握擬選機(jī)場場址提供了很好的工程評價(jià)平臺。同時(shí),對加快機(jī)場建設(shè)工作現(xiàn)代化進(jìn)程,具有重要意義和實(shí)用價(jià)值。

[1]楊武年,廖崇高,濮國梁,等.數(shù)字區(qū)調(diào)新技術(shù)新方法—遙感圖像地質(zhì)解譯三維可視化及影像動態(tài)分析[J].地質(zhì)通報(bào).2003,22(1):60.

[2]賈永紅,李德仁,孫家柄.多源遙感影像數(shù)據(jù)融合[J].遙感技術(shù)與應(yīng)用,2001,15(1):41.

[3]張瑞軍,楊武年,劉漢湖,等.數(shù)字高程模型(DEM)的構(gòu)建及其應(yīng)用[J].工程勘察,2005,196:(5):61.

[4]YANGWU-N IAN,CAUNEAU F,PAR IS J-P,et al.Fusion of SAR and SPOT im ages for the detection of geo logical featuresover the Three GorgesDam site,China.In Proceedingsof the third conference"Fusion of Earth data:m erging po in tm easurem en ts,rasterm ap s and remotely sensed im ages"[C].Sophia Antipo lis,France,January 26-28,2000,T.Ranchin and L.W ald Editors,pub lished by SEE/UR ISCA,N ice,France,2000:137.

[5]曹代勇,李青元,朱小弟,等.地質(zhì)構(gòu)造三維可視化模型探討[J].地質(zhì)與勘探,2001,37(4):60.

[6]龔建華,林琿,肖樂斌,等.地學(xué)可視化探討[J].遙感學(xué)報(bào),1999,3(3):236.

[7]劉志麗,陳曦.基于ERDAS IMAGING軟件的TM影像幾何精校正方法初探[J].干旱區(qū)地理,2001,12:353.

[8]濮國梁,楊武年,劉漢超.正射遙感影像地圖制作技術(shù)在岷江上游滑坡研究中的應(yīng)用[J].成都理工學(xué)院學(xué)報(bào),2000,27(3):312.

[9]劉曉龍.基于影像匹配接邊糾正的數(shù)字正射影像的鑲嵌技術(shù)[J].遙感學(xué)報(bào),2001,3:104.

[10]曹建君,苗天寶,蔡喜琴,等.SPOT5衛(wèi)星影像在退耕精確監(jiān)測的應(yīng)用研究[J].遙感技術(shù)與應(yīng)用,2004:131.

[11]楊武年,FRANO ISCAUNEAU,TH IERRYRANCH IN,etal.應(yīng)用多種數(shù)據(jù)融合方法對多時(shí)相SPOT圖像處理提取三峽工程大壩區(qū)地面動態(tài)變化信息[J].物探化探計(jì)算技術(shù),2005,27(3):267.

[12]蔣勇軍,袁道先,謝世友.典型巖溶農(nóng)業(yè)區(qū)地下水質(zhì)與土地利用變化分析～以云南小江流域?yàn)槔齕J].地理學(xué)報(bào),2006,61(5):471.

[13]李天華,廖崇高,楊武年.遙感數(shù)字圖像處理關(guān)鍵技術(shù)及其在青藏高原機(jī)場建設(shè)工程中的應(yīng)用[J].工程勘察,2006,(12):51.

TP 75

A

1001—1749(2011)01—0079—04

國家自然基金項(xiàng)目(60841006);云南政府重點(diǎn)建設(shè)項(xiàng)目(20070616009)

2010-06-30 改回日期:2010-11-03

李天華(1961-),男,博士,高級工程師,研究方向:“3S”技術(shù)及應(yīng)用。