基于InSAR技術(shù)的內(nèi)蒙古巴彥淖爾市地面沉降演化特征及成因分析

葛偉麗 ，李元杰 ，張春明 ，張紅霞 ，王志超 ，楊紅磊

（1.內(nèi)蒙古自治區(qū)地質(zhì)調(diào)查研究院, 內(nèi)蒙古 呼和浩特 010020；2.內(nèi)蒙古煤炭建設(shè)生態(tài)環(huán)境研究院有限責(zé)任公司, 內(nèi)蒙古 呼和浩特 010020；3.內(nèi)蒙古自治區(qū)測(cè)繪地理信息中心, 內(nèi)蒙古 呼和浩特 010020；4.中國(guó)地質(zhì)大學(xué)（北京）土地科學(xué)技術(shù)學(xué)院, 北京 100083）

地面沉降是指由于自然因素或人類工程活動(dòng)引發(fā)的松散層固結(jié)壓縮并導(dǎo)致一定區(qū)域范圍內(nèi)地面高程降低的地質(zhì)現(xiàn)象，是一種不可補(bǔ)償?shù)挠谰眯原h(huán)境和資源損失[1-2]。地震、采礦、石油、天然氣或地下水的過(guò)度開(kāi)采以及構(gòu)造運(yùn)動(dòng)等各種因素都可能導(dǎo)致地面沉降[3-4]。長(zhǎng)期地下水超采導(dǎo)致松散沉積物中顆粒間孔隙流體壓力降低從而引發(fā)的地面沉降，已是世界范圍內(nèi)各大城市面臨的主要地質(zhì)環(huán)境問(wèn)題之一[5-7]。

在地面沉降測(cè)量方面，傳統(tǒng)的水準(zhǔn)測(cè)量手段和GPS監(jiān)測(cè)更適合離散點(diǎn)的監(jiān)測(cè)，對(duì)于地域范圍跨度較大、沉降變化較為緩慢且時(shí)間較長(zhǎng)的城市地面沉降監(jiān)測(cè)，這兩種方法都具有其局限性[8]。合成孔徑雷達(dá)干涉測(cè)量技術(shù)（Interferometric Synthetic Aperture Radartechnique, InSAR）以其全天候、全天時(shí)、精度高、覆蓋范圍大的優(yōu)點(diǎn)被廣泛應(yīng)用[9-10]，其中永久散射體干涉法（Permanent Scatterer Intervention Method，PS-InSAR）和小基線集方法（Small Baseline Set Method，SBAS-InSAR）作為先進(jìn)的InSAR技術(shù)方法，可以克服大氣延遲和時(shí)空非相干性等缺陷，獲得高精度的沉降反演結(jié)果，已被許多國(guó)家用于城市地面沉降監(jiān)測(cè)[11-12]。此外，高分辨率SAR圖像也得到廣泛應(yīng)用，如搭載相控陣型L波段合成孔徑雷達(dá)（PALSAR）傳感器的ALOS衛(wèi)星（ALOSPALSAR）和搭載C波段傳感器的Radarsat-2衛(wèi)星的SAR圖像[13-14]。

巴彥淖爾市位于內(nèi)蒙古自治區(qū)西部，近年來(lái)隨著人口激增和城市建設(shè)加快，地下水降落漏斗范圍不斷擴(kuò)大，且該地區(qū)固有的第四系粉砂層具有發(fā)育廣、厚度大的特點(diǎn)，易產(chǎn)生不均勻地面沉降，且沉降持續(xù)時(shí)間較長(zhǎng)[15]，因此巴彥淖爾市具備地面沉降產(chǎn)生的基礎(chǔ)條件，然而目前該市地面沉降的區(qū)域性研究尚未開(kāi)展，已有的水平面測(cè)點(diǎn)近年來(lái)也遭到嚴(yán)重破壞。因此，本文利用PS-InSAR和SBAS-InSAR技術(shù)首次對(duì)巴彥淖爾市2007—2011年和2015—2016年的地面沉降演化特征進(jìn)行定量研究分析，通過(guò)現(xiàn)場(chǎng)實(shí)地調(diào)查驗(yàn)證InSAR技術(shù)監(jiān)測(cè)結(jié)果的可靠性，并進(jìn)一步分析研究區(qū)地面沉降成因及發(fā)展趨勢(shì)，以期為地面沉降防災(zāi)減災(zāi)提供技術(shù)支撐。

1 研究區(qū)概況

巴彥淖爾市位于黃河沖積平原上，總面積65 139 km2，地面開(kāi)闊平坦，海拔1 045～1 209 m，地勢(shì)由西南向東北傾斜。地表出露多為第四系全新統(tǒng)沖積砂及砂質(zhì)粉土，厚6～51 m。第四系厚約1 000 m，推測(cè)最大厚度為1 200～1 500 m。第四系沉積層厚度大、含水率高、固結(jié)程度差，在重力等作用下，易發(fā)生地層自然壓實(shí)、固結(jié)。構(gòu)造上屬于河套盆地，西界為郎山前緣斷裂，北界為色爾騰山斷裂，南界為鄂爾多斯北緣斷裂?？傮w而言，區(qū)域地殼處于相對(duì)穩(wěn)定狀態(tài)。研究區(qū)內(nèi)有2條斷裂帶：1條分布在東北方向，從晨光村到高峰村；另1條分布在民豐村—友誼村—扶豐村，沿東西方向分布。區(qū)內(nèi)淺層地下水分布普遍，為主要開(kāi)采層[16]。含水層埋藏厚度受構(gòu)造控制明顯，盆地邊緣較薄，深拗陷帶較厚，臨河區(qū)一般為55～60 m，粉砂層厚度自東向西、自南向北逐漸增加，埋深相應(yīng)加深。粉砂層厚度在北部杭錦后旗地區(qū)最大，達(dá)280 m[17]（圖1—2）。

圖1 研究區(qū)地理位置及第四系粉砂層厚度圖Fig.1 Geographical location of the study area and depth of the Quaternary silt strata

圖2 I—I′水文地質(zhì)剖面圖Fig.2 Hydrogeological profile of section I—I′

2 數(shù)據(jù)來(lái)源及研究方法

2.1 數(shù)據(jù)來(lái)源

本研究采用ALOS PALSAR（2007-01-13—2011-03-11，98景）和Radarsat-2（2015-06-28—2016-12-31，10景）圖像。ALOS PALSAR數(shù)據(jù)較豐富，數(shù)據(jù)采樣較均勻，L波段數(shù)據(jù)波長(zhǎng)較長(zhǎng)、穿透性較好，在非建筑區(qū)仍有相干性。目前在軌衛(wèi)星有Radarsat-2、TerraSAR-X、COSMOSkyMed和Sentinel-1，其中Radarsat-2從2015年開(kāi)始，在研究區(qū)有拍攝數(shù)據(jù)，TerraSAR-X和COSMO-SkyMed在研究區(qū)無(wú)存檔數(shù)據(jù)。綜合對(duì)比已有數(shù)據(jù)情況、數(shù)據(jù)幅寬和波長(zhǎng)等指標(biāo)，最終選用ALOS PALSAR和Radarsat-2數(shù)據(jù)分別對(duì)巴彥淖爾市2007—2011年和2015—2016年的地面歷史沉降情況進(jìn)行分析研究。采集參數(shù)如表1所示。

表1 巴彥淖爾市ALOS PALSAR 和Radarsat-2 數(shù)據(jù)參數(shù)Table 1 Parameters of ALOS PALSAR data and Radarsat-2 SAR data in the city of Bayannur

2.2 研究方法

采用PS-InSAR方法對(duì)2007-01-13—2011-03-11期間巴彥淖爾市的地面沉降進(jìn)行識(shí)別。PS-InSAR技術(shù)本質(zhì)上屬于雷達(dá)差分干涉法，將測(cè)量區(qū)域的序列SAR圖像對(duì)應(yīng)的地面目標(biāo)分為永久散射（Permanent Scatterer，PS）點(diǎn)和非永久散射點(diǎn)兩類，然后提取PS點(diǎn)，采用差分干涉法和建模方法實(shí)現(xiàn)信號(hào)相位的分離，最后提取PS變形信息并進(jìn)行插值得到研究區(qū)域的地面形變情況[18]，如圖3（a）所示。選取2009年9月22日的圖像為主影像，與其他影像組合24個(gè)干涉像對(duì)，得到地表形變信息。采用10景Radarsat-2數(shù)據(jù)對(duì)2015-06-02—2016-12-31期間巴彥淖爾市的地面沉降情況進(jìn)行簡(jiǎn)單研究，而PS-InSAR方法要求SAR數(shù)據(jù)量大于25景，因此采用SBAS-InSAR方法研究2015-06-02—2016-12-31期間的地面沉降情況。以時(shí)間基線不超過(guò)350 d、空間基線不長(zhǎng)于560 m為原則，組成干涉對(duì)，得到該地區(qū)的地表沉降速率，具體處理過(guò)程如圖3（b）所示。

圖3 數(shù)據(jù)處理流程Fig.3 Workflow of InSAR data processing

3 研究結(jié)果

對(duì)2007—2011年98幅ALOS PALSAR數(shù)據(jù)的SAR圖像進(jìn)行干涉處理，最終獲取PS點(diǎn)。巴彥淖爾市的地面沉降主要表現(xiàn)為緩慢性沉降，沉降速率集中在0～10 mm/a。所有PS點(diǎn)中最大沉降速率達(dá)23.7 mm/a，位于臨河區(qū)喬和圪旦的空地上。以InSAR影像解譯結(jié)果為基礎(chǔ)，選擇InSAR解譯成果沉降速率較大的PS點(diǎn)集中分布區(qū)，結(jié)合城市地下水開(kāi)采情況、區(qū)域地下水水位動(dòng)態(tài)變化特征、地下水降落漏斗形態(tài)特征、第四系粉砂層發(fā)育特征等條件，并進(jìn)行野外驗(yàn)證，綜合確定地面沉降區(qū)的位置及范圍。結(jié)果表明2007—2011年巴彥淖爾市存在3處地面沉降集中區(qū)，即杭錦后旗沉降區(qū)（A）、臨河區(qū)沉降區(qū)（B）、開(kāi)發(fā)區(qū)水源地及北部沉降區(qū)（C），如表2和圖4（a）所示。2015—2016年巴彥淖爾市的沉降也為緩慢性沉降，沉降速率集中在0～2.6 mm/a；沉降速率最大值為5.634 mm/a，位于臨河區(qū)光明村中央路口以北0.3 km處；且存在1處地面沉降區(qū)，即臨河區(qū)及北部沉降區(qū)（D），如表2和圖4（b）所示。巴彥淖爾市地面沉降整體表現(xiàn)為緩慢性沉降，按《地面沉降干涉雷達(dá)數(shù)據(jù)處理技術(shù)規(guī)程》（DD 2014—11）中地面沉降嚴(yán)重程度分級(jí)標(biāo)準(zhǔn)，巴彥淖爾市地面沉降嚴(yán)重程度屬于低級(jí)。

圖4 巴彥淖爾市地面形變速度Fig.4 Ground deformation velocity in the city of Bayannur

表2 巴彥淖爾市地面沉降集中區(qū)特征Table 2 Characteristics of land subsidence concentration area in the city of Bayannur

實(shí)地調(diào)查發(fā)現(xiàn)，在杭錦后旗陜壩鎮(zhèn)西二街村房屋東墻發(fā)現(xiàn)多處裂縫，該點(diǎn)距離杭錦后旗地下水超采區(qū)約1 km；同時(shí)在第一水源地4號(hào)水源井、第二水源地6號(hào)水源井發(fā)現(xiàn)井管相對(duì)抬升、井周存在輻射狀裂縫，其他地區(qū)未見(jiàn)明顯地面沉降明顯跡象（圖5）。這是由于巴彥淖爾市地面沉降尚處于初步發(fā)生階段，地面沉降嚴(yán)重程度屬于低級(jí)，地面沉降發(fā)展較緩，因此未出現(xiàn)較大范圍的區(qū)域性地面沉降跡象。

圖5 巴彥淖爾市地面沉降特征野外調(diào)查照片F(xiàn)ig.5 Field survey photos of land subsidence characteristics in the city of Bayannur

InSAR解譯影像顯示，部分地區(qū)存在地面抬升現(xiàn)象，本研究在野外核查工作中有針對(duì)性地對(duì)地面抬升特征區(qū)域進(jìn)行了實(shí)地核查，核查點(diǎn)分布情況見(jiàn)圖4（a）。結(jié)果表明，除BN1調(diào)查點(diǎn)位于臨河區(qū)團(tuán)結(jié)路萬(wàn)豐小區(qū)外，其余調(diào)查點(diǎn)均位于農(nóng)田、林地、草地、蘆葦?shù)?，受作物種類、長(zhǎng)勢(shì)影響，對(duì)InSAR解譯成果造成一定干擾。而B(niǎo)N1調(diào)查點(diǎn)處的萬(wàn)豐小區(qū)始建于2010年，當(dāng)時(shí)地表人工開(kāi)挖擾動(dòng)較大，因此2007—2011年InSAR解譯影像顯示的地表回升現(xiàn)象與此有關(guān)。此外，2015—2016年黃河西南地區(qū)遙感解譯影像顯示，研究區(qū)出現(xiàn)大面積的地面沉降集中點(diǎn)；隨后野外驗(yàn)證表明，該區(qū)域位于庫(kù)布齊沙漠，地表變形易受水土風(fēng)蝕等自然因素的影響，因此，推測(cè)這些沉降點(diǎn)為地面沉降干擾點(diǎn)；而且由于該區(qū)域隸屬于鄂爾多斯市，不在本研究區(qū)范圍內(nèi)，因此未進(jìn)行詳細(xì)分析研究。

4 地面沉降成因及趨勢(shì)分析

4.1 地面沉降成因分析

4.1.1 自然因素

研究表明地面沉降量主要來(lái)自含水層頂板自身的固結(jié)壓縮[19]。2007—2011年巴彥淖爾市A、B、C等3處地面沉降集中區(qū)沉降面積較大，其中B沉降面積達(dá)139.98 km2；2015—2016年該沉降區(qū)繼續(xù)向周邊及北部擴(kuò)展，沉降面積增至221.43 km2。由于上述地面沉降區(qū)均位于研究區(qū)粉砂層厚度較大的區(qū)域，因此認(rèn)為巴彥淖爾市沉積物的固結(jié)壓實(shí)是導(dǎo)致該地區(qū)地面沉降的重要因素。為了進(jìn)一步探討巴彥淖爾市沉積物固結(jié)壓實(shí)對(duì)地面沉降的影響，本研究比較分析了2007—2011年3個(gè)地面沉降區(qū)平均沉降速率與相應(yīng)粉砂層厚度的相關(guān)性。結(jié)果表明，粉砂層厚度最薄的C區(qū)沉降速率最低，為4.615 mm/a；中等粉砂層厚度的B區(qū)沉降速率最高，為8.336 mm/a；而A區(qū)粉砂層厚度最大，達(dá)260～280 m，沉降速率卻較小。這說(shuō)明沉積物固結(jié)是引發(fā)該區(qū)地面沉降的原因之一，地面沉降是多種因素共同作用的結(jié)果。在本研究中，巴彥淖爾市兩條斷裂兩側(cè)未見(jiàn)明顯差異變形，斷裂兩側(cè)也未見(jiàn)地面沉降的帶狀分布。此外，由構(gòu)造活動(dòng)引起的地表高程差異發(fā)生在漫長(zhǎng)的歷史時(shí)期，速度極慢，影響很小[20]。因此，不考慮構(gòu)造活動(dòng)對(duì)研究區(qū)的影響。

4.1.2 人為因素

2007—2011年巴彥淖爾市存在的3個(gè)地面沉降集中區(qū)中，A區(qū)位于陜壩鎮(zhèn)水源地的東北方向，該區(qū)地面沉降呈小范圍小幅度均勻沉降，沉降中心與杭錦后旗地下水超采區(qū)及陜壩鎮(zhèn)水源地在區(qū)域上較為接近，呈現(xiàn)出與地下開(kāi)采較為相關(guān)的特性；B區(qū)內(nèi)有第一水廠、第二水廠、鐵路水廠，企事業(yè)單位自備井較多，目前區(qū)內(nèi)已形成臨河區(qū)小型孔隙淺層地下水超采區(qū)，該沉降區(qū)位于城區(qū)地下水降落漏斗中部，區(qū)內(nèi)地面沉降呈大范圍小幅度均勻沉降，地面沉降中心與地下水降落漏斗中心在空間分布上具有較高的一致性，但沉降中心與地下水漏斗中心并非完全吻合，沉降漏斗有向東北偏移的趨勢(shì)，如圖6（a）所示。其主要原因可能是軟土層固結(jié)速度滯后于地下水水頭變化，并且不同地區(qū)的軟土層厚度對(duì)該區(qū)域地面沉降的發(fā)生、發(fā)展也有顯著影響[21-22]。2004—2015年地下水降落漏斗范圍不斷擴(kuò)大，可能導(dǎo)致地面沉降范圍逐步擴(kuò)大，本次對(duì)巴彥淖爾市2015—2016年的地面沉降調(diào)查也恰好驗(yàn)證了這一發(fā)展趨勢(shì)：B區(qū)在原有沉降基礎(chǔ)上繼續(xù)向周邊擴(kuò)展，特別是東北部，北至白虎高圪旦，東到長(zhǎng)豐村，沉降區(qū)面積達(dá)221.43 km2，如圖6（b）所示；C區(qū)東南緊鄰開(kāi)發(fā)區(qū)水源地，由于2010年后開(kāi)發(fā)區(qū)工業(yè)產(chǎn)值增加較快，導(dǎo)致區(qū)內(nèi)企業(yè)用水增多，地下水位下降相對(duì)較快，因此C區(qū)地面沉降速率相對(duì)較高，選取巴彥淖爾市2007—2011年兩個(gè)具有完整水位監(jiān)測(cè)數(shù)據(jù)的形變特征點(diǎn)，觀測(cè)其地面形變速度和地下水位，并對(duì)二者的相關(guān)性進(jìn)行回歸分析。兩個(gè)特征點(diǎn)分別位于臨河區(qū)城關(guān)鎮(zhèn)五四十社村和臨河區(qū)聯(lián)邦制藥（內(nèi)蒙古）有限公司，對(duì)應(yīng)的水位監(jiān)測(cè)點(diǎn)分別為B120和B128-1，具體位置如圖6（a）所示。結(jié)果表明，地下水位與地面形變速度具有一定的相關(guān)性，水位下降時(shí)地面沉降，水位上升時(shí)地面抬升，地面形變與地下水位變化密切相關(guān)，各點(diǎn)兩者的相關(guān)系數(shù)均大于0.8（圖7）。因此，地下水開(kāi)采是巴彥淖爾市地面沉降的主要影響因素。

圖6 地面沉降環(huán)境影響因素圖Fig.6 Diagram showing the environmental influencing factors in land subsidence areas

圖7 地面形變速度與承壓水水位相關(guān)性Fig.7 Correlation of the ground deformation velocity and artesian groundwater level elevation

2007—2011年，巴彥淖爾市城市建設(shè)進(jìn)入快速增長(zhǎng)階段，新增建筑物較多。地基建設(shè)往往需要在開(kāi)挖過(guò)程中抽取地下水，造成地基下沉。本研究查明的幾個(gè)地面沉降集中區(qū)，不僅是地下水開(kāi)發(fā)利用程度較高的區(qū)域，同樣為居民聚居區(qū)（如A區(qū)和B區(qū)）或新增建筑較多區(qū)域（如C區(qū)）。分別在3個(gè)地面沉降集中區(qū)選取沉降點(diǎn)M1、M2、M3，基于歷史衛(wèi)星圖像，分析2007—2011年沉降點(diǎn)及沉降點(diǎn)周圍建筑物變化情況（圖8）。結(jié)果表明，3個(gè)典型的地面沉降點(diǎn)及周圍均為新增建筑物較多的地區(qū)，因此地表荷載增加也是巴彥淖爾市地面沉降的影響因素之一。

圖8 各沉降區(qū)研究期始末歷史衛(wèi)星影像Fig.8 Historical satellite images at the beginning and the end of the study period in each land subsidence area

4.2 地面沉降發(fā)展趨勢(shì)分析

2007—2011年和2015—2016年巴彥淖爾市地面沉降速率分別集中在0～10 mm/a和0～2.6 mm/a。由地面沉降階段性預(yù)測(cè)標(biāo)準(zhǔn)[23]可知巴彥淖爾市地面沉降尚屬于發(fā)生階段。盡管本次調(diào)查顯示2015—2016年A區(qū)沒(méi)有繼續(xù)擴(kuò)展或加重的趨勢(shì)，但據(jù)《巴彥淖爾市城市總體規(guī)劃（2011—2030年）》，陜壩鎮(zhèn)規(guī)劃保留現(xiàn)狀自來(lái)水廠，并新建二水廠，因此若不合理開(kāi)采地下水，未來(lái)該沉降區(qū)很有可能向東南部擴(kuò)展至陜壩鎮(zhèn)水源地保護(hù)區(qū)。B區(qū)內(nèi)設(shè)有第一水廠、第二水廠、鐵路水廠，企事業(yè)單位自備井較多，據(jù)《巴彥淖爾市城市總體規(guī)劃（2011—2030）》，中心城區(qū)保留現(xiàn)狀一水廠、二水廠、三水廠及鐵路水廠，并在規(guī)劃期末達(dá)到設(shè)計(jì)供水規(guī)模，新建四水廠；未來(lái)隨著地下水開(kāi)采強(qiáng)度增加，中心城區(qū)降落漏斗可能繼續(xù)擴(kuò)大，致使該沉降區(qū)向東部擴(kuò)展至第三水廠，沉降速率有所加快。C區(qū)東南緊鄰第三水廠，近年由于臨河?xùn)|部巴彥淖爾經(jīng)濟(jì)開(kāi)發(fā)區(qū)建設(shè)和開(kāi)發(fā)區(qū)最大企業(yè)“聯(lián)邦（內(nèi)蒙古）制藥公司”的入駐，開(kāi)發(fā)區(qū)已新建水源地（35眼井）和聯(lián)邦制藥自備水源地（58眼井），未來(lái)區(qū)內(nèi)企事業(yè)單位新增用水量仍會(huì)增加，且已有研究表明，隨著建筑荷載的增大沉降范圍也會(huì)增大，當(dāng)建筑荷載足夠大且建筑物間距過(guò)小時(shí)，地面沉降漏斗可能連成一片。因此，未來(lái)隨著地下水用水量增加及開(kāi)發(fā)區(qū)大規(guī)模城鎮(zhèn)建設(shè)，該沉降區(qū)可能呈輻射狀向周邊擴(kuò)展，西南部與臨河區(qū)沉降區(qū)相連通，沉降速率也有所增加。

5 結(jié)論及建議

（1）巴彥淖爾市地面沉降呈緩慢性沉降。2007—2011年沉降速率集中在0～10 mm/a，存在3處地面沉降集中區(qū)，即杭錦后旗沉降區(qū)、臨河區(qū)沉降區(qū)、開(kāi)發(fā)區(qū)水源地及北部沉降區(qū)；2015—2016年，沉降速率主要是集中在0 ～2.6 mm/a，且存在1處地面沉降區(qū)，即臨河區(qū)及北部沉降區(qū)，該沉降區(qū)是在2007—2011年形成的臨河區(qū)沉降區(qū)基礎(chǔ)上繼續(xù)向周邊及北部擴(kuò)展形成的。

（2）巴彥淖爾市緩慢性地面沉降的主要影響因素為地下水的開(kāi)發(fā)利用引發(fā)的地下水水位下降，沉積固結(jié)壓實(shí)也是該地區(qū)地面沉降產(chǎn)生的重要因素，地表荷載增加也產(chǎn)生了一定推動(dòng)作用。

（3）巴彥淖爾市地面沉降屬于發(fā)生階段。未來(lái)杭錦后旗沉降區(qū)在原有沉降基礎(chǔ)上，有可能向東南方向擴(kuò)展至陜壩鎮(zhèn)水源地保護(hù)區(qū)。此外，臨河區(qū)沉降區(qū)可能呈現(xiàn)出向東部擴(kuò)展至第三水廠，沉降速率有加快的趨勢(shì)。開(kāi)發(fā)區(qū)水源地及北部沉降區(qū)可能呈輻射狀向周邊擴(kuò)展，西南部與臨河區(qū)沉降區(qū)相連通，沉降速率也會(huì)有所增加。

建議將巴彥淖爾市地面沉降的下一步研究重點(diǎn)放在長(zhǎng)時(shí)間序列自然因素與短時(shí)間序列人為因素各自影響的區(qū)分分析上，并建立地面沉降監(jiān)測(cè)系統(tǒng)，開(kāi)發(fā)區(qū)域地面沉降預(yù)測(cè)模型，為地質(zhì)環(huán)境保護(hù)、城市基礎(chǔ)設(shè)施建設(shè)、城市發(fā)展規(guī)劃、重要工程選址提供地質(zhì)環(huán)境依據(jù)和參考。