辛 超, 黃 超, 付榮翔, 龐 詠, 黃 輝, 葉 雷, 劉 達(dá)

(1.江西省地質(zhì)調(diào)查勘查院礦產(chǎn)勘查所,南昌 330000;2.中國(guó)地質(zhì)調(diào)查局 昆明自然資源綜合調(diào)查中心,昆明 650100)

0 引言

富源縣位于云南省地質(zhì)災(zāi)害次重點(diǎn)防治區(qū)中的富源-師宗滑坡次重點(diǎn)防治亞區(qū),區(qū)內(nèi)地質(zhì)災(zāi)害以滑坡為主[1]。半坡村滑坡位于富源縣恩洪-老廠煤層勘查區(qū)內(nèi),地處大則勒煤礦礦界附近,該滑坡地裂縫每年持續(xù)擴(kuò)張,多次破壞道路,對(duì)半坡村產(chǎn)生多處房屋裂縫、地面裂縫、地面塌陷等問(wèn)題。雖然半坡村于2018年6月整村搬遷,但半坡村滑坡仍然處于活動(dòng)狀態(tài),威脅其東部的S204省道,探明其地下空間形態(tài)、滑面位置、成因機(jī)制對(duì)以后的災(zāi)害治理具有重要的指導(dǎo)意義。

滑坡的基本特征主要包括平面和剖面形態(tài)特征、滑體組成與結(jié)構(gòu)及滑面特征,是研究其變形特征和成因機(jī)制的重要依據(jù)[2]。長(zhǎng)期以來(lái),傳統(tǒng)的滑坡調(diào)查手段多以觀測(cè)水位變化、進(jìn)行位移監(jiān)測(cè)、形變分析、調(diào)查誘發(fā)原因以及降雨量等因素為主,而滑移面形態(tài)特征主要以通過(guò)鉆孔資料繪制的地質(zhì)剖面為依據(jù)[3]。鉆孔資料能夠清楚地反映地層結(jié)構(gòu),但由于施工成本大,無(wú)法大面積布設(shè),只能局部布設(shè)。通過(guò)地球物理勘察,能夠真實(shí)反映滑坡體的地下空間結(jié)構(gòu)及滑面位置。 加以鉆孔進(jìn)行驗(yàn)證,既能保證工作精度又能控制施工成本。

隨著現(xiàn)在地球物理理論的研究不斷深入及儀器設(shè)備的勘察精度不斷提高,不少學(xué)者將地球物理勘察方法引入到滑坡調(diào)查當(dāng)中,對(duì)滑坡的基本特征進(jìn)行勘察。能運(yùn)用于滑坡勘察的地球物理方法有地震法(地震反射法、地震折射法、地震 面波法、微動(dòng)勘查和陸地聲納法)、電法(高密度電阻率法、自然電位法和激發(fā)極化法)及電磁法(地質(zhì)雷達(dá))[4]。楊成林[5]等運(yùn)用探地雷達(dá)對(duì)趙子秀山滑坡進(jìn)行勘探;程慶[6]等運(yùn)用高密度電阻率發(fā)對(duì)四川高川茶園溝滑坡進(jìn)行勘探;許新剛等運(yùn)用瑞利面波法對(duì)秦嶺一小型滑坡進(jìn)行勘察;金維民[8]等運(yùn)用等淺層地震勘探方法對(duì)下垢坪斜坡進(jìn)行勘察。近10年來(lái),高密度電法是滑坡勘察主要工作手段[9-15],利用高密度電法查明滑移面形態(tài)、分析滑坡工程地質(zhì)特點(diǎn)和滑坡的穩(wěn)定性評(píng)價(jià)[3]。

首先通過(guò)res2dmod軟件構(gòu)建滑坡模型進(jìn)行數(shù)值模擬,使用res2dinv軟件對(duì)正演數(shù)據(jù)進(jìn)行反演計(jì)算。通過(guò)反演結(jié)果與模型進(jìn)行對(duì)比,從理論上驗(yàn)證高密度電法在滑坡勘察中的有效性。以半坡村滑坡為例,在滑坡范圍內(nèi)布設(shè)高密度電法剖面,查明滑坡體地層結(jié)構(gòu)、滑床形態(tài)、滑面埋深等信息,通過(guò)工程鉆孔驗(yàn)證反演結(jié)果的可靠性,構(gòu)建工程地質(zhì)剖面,研究誘發(fā)形成機(jī)理,分析致災(zāi)成因。

1 高密度電法工作原理

高密度電法屬于直流電阻率法的一種,通過(guò)A、B供電電極向地下供電,用M、N測(cè)量電極測(cè)量電位差,以計(jì)算出視電阻率ρs。

(1)

高密度電法是一種陣列勘探方法,能夠一次布設(shè)幾十至上百根電極,利用程控電極轉(zhuǎn)換開關(guān),能夠?qū)崿F(xiàn)數(shù)據(jù)的自動(dòng)采集。高密度電法裝置已多達(dá)十多種,其中常用的有溫納裝置和斯倫貝謝裝置(圖1)。這兩種視電阻率計(jì)算公式分別為:

圖1 溫納裝置、施倫貝謝裝置電極排列示意圖Fig.1 Diagram of electrode arrangement of Wenner device and Schlumberger device

(2)

(3)

2 數(shù)值模擬

為驗(yàn)證高密度電法在滑坡勘察中的有效性,利用res2dmod軟件進(jìn)行數(shù)值模擬,構(gòu)建滑坡模型,剖面長(zhǎng)度為100 m,剖面布設(shè)101根電極,電極距為1 m。采用有限差分法進(jìn)行網(wǎng)格剖分:網(wǎng)格列間距采用每相鄰兩電極為4個(gè)網(wǎng)格單元等間距剖分;網(wǎng)格行間距采用不等間距剖分,共分17層,分別為0.25 m、0.5 m、0.75 m、1 m、1.5 m、2 m、2.5 m、3 m、3.5 m、4 m、4.5 m、5.5 m、10.4 m、15.4 m、30.4 m、60.4 m、120.4 m。

構(gòu)建如圖2所示模型:上部為地表風(fēng)化層,電阻率為120 Ω·m,位置在5.5 m以淺;下部為基巖層,電阻率為300 Ω·m,位置在5.5 m以深;上部存在三條滑面,電阻率為10 Ω·m,底部埋深分別為1 m、2.5 m和5.5 m;在滑面上為松散層,電阻率為100 Ω·m。

圖2 模型剖面及電性參數(shù)Fig.2 Model profile and electrical parameters

利用res2dinv軟件,采用最小二乘法對(duì)正演數(shù)據(jù)(圖3)進(jìn)行反演計(jì)算,從反演擬斷面圖(圖4)可知,高密度電法能夠有效識(shí)別多期次滑坡的多條滑面位置,并且隨著滑面埋深越深,反演滑面埋深位置越接近實(shí)際模型(表1)。

表1 滑面埋深誤差對(duì)比表Tab.1 Comparison table of buried depth error of sliding surface

圖3 模型正演擬斷面圖Fig.3 Simulation section of the forward model

圖4 模型反演擬斷面圖Fig.4 Simulation section of the inversion model

3 應(yīng)用實(shí)例

如圖5所示,半坡村滑坡地理位置為,位于富源縣營(yíng)上鎮(zhèn)塊抉河右岸一級(jí)支流左岸半坡村,東部為S204省道,受構(gòu)造活動(dòng)及地層巖性的控制,周圍地形地貌極為復(fù)雜地形地貌復(fù)雜,總體以溶蝕-侵蝕地貌(中切山地)和侵蝕構(gòu)造地貌(中山剝蝕-侵蝕地貌)為主,山脈和河流多呈北北東向展布。滑坡后部山脊高程為1 896.1 m,滑坡后緣高程為1 883 m,滑坡前緣高程為1 673 m,滑坡高差為210 m,主滑方向?yàn)?40°,滑坡體長(zhǎng)約500 m,寬約670 m。通過(guò)測(cè)量估算,滑體表面積為33×104m2,體積約396×104m3。

圖5 半坡村滑坡平面圖Fig.5 Landslide plan of Banpo Village

3.1 測(cè)區(qū)地質(zhì)概況

根據(jù)1996年出版的《云南省巖石地層》和2000年《中國(guó)地層指南及中國(guó)地層指南說(shuō)明書》,結(jié)合測(cè)區(qū)不同構(gòu)造單元內(nèi)部的構(gòu)造特征、地層層序及生物組合特征,測(cè)區(qū)地層隸屬于華南地層大區(qū)之揚(yáng)子地層區(qū)的曲靖地層小區(qū)。根據(jù)地質(zhì)剖面資料揭示,測(cè)區(qū)內(nèi)出露的主要地層為:峨眉山組玄武巖、茅口組灰?guī)r、宣威組砂巖、泥巖。

3.2 測(cè)區(qū)地球物理概況

物性測(cè)量工作采用露頭測(cè)量、標(biāo)本測(cè)量及水文鉆孔三側(cè)向電阻率測(cè)量等3種物性測(cè)定方式,共測(cè)量4種巖性電阻率,露頭測(cè)量123個(gè)點(diǎn)(圖6),測(cè)量標(biāo)本120塊(圖7)。從測(cè)量結(jié)果來(lái)看(圖8),各測(cè)定方法測(cè)量的各巖性電阻率高低分布特征基本一致,泥巖、砂巖電阻率基本一致,玄武巖電阻率略大于砂巖、泥巖電阻率,灰?guī)r電阻率遠(yuǎn)大于其他巖石電阻率。綜合考慮,各電阻率從小到大以此為泥巖、砂巖、玄武巖、灰?guī)r。

圖6 露頭測(cè)量各巖性電阻率分布圖Fig.6 Resistivity distribution of each lithology measured on outcrop

圖7 標(biāo)本測(cè)量各巖性電阻率分布圖Fig.7 Resistivity distribution of each lithology measured by specimen

圖8 物性測(cè)量各巖性電阻率平均值分布圖Fig.8 Average resistivity distribution of each lithology in physical property measurement

3.3 野外數(shù)據(jù)采集

為查明半坡村滑坡的基本特征及誘發(fā)形成機(jī)理,沿半坡村滑坡走向布設(shè)1條長(zhǎng)1 195 m的高密度電法剖面L1,使用深圳賽盈地脈公司生產(chǎn)的GD-10高密度電阻率法測(cè)量系統(tǒng),采用溫納裝置方式進(jìn)行測(cè)量,最小隔離系為1,最大隔離系數(shù)為39。為滿足測(cè)量分辨率要求,采用5 m電極距進(jìn)行測(cè)量。由于測(cè)線較長(zhǎng),采用滾動(dòng)排列模式進(jìn)行測(cè)量,第一次布設(shè)108個(gè)電極,每次滾動(dòng)12個(gè)電極,共布設(shè)240個(gè)電極(圖9)。

圖9 高密度電法視電阻率數(shù)據(jù)剖面圖Fig.9 Apparent resistivity data profile of high density electrical method

3.4 成果資料解釋

從WT205線反演圖(圖11)可知,剖面視電阻率值主要表現(xiàn)為西高東低的特征,變化范圍在3.8 Ω·m～3 507 Ω·m之間,電阻率幅值較為分散,反映巖性較為復(fù)雜。通過(guò)與地質(zhì)剖面對(duì)比,高密度電法能夠準(zhǔn)確的反映出各巖性界線位置,各巖性視電阻率與物性測(cè)量結(jié)果基本一致,玄武巖電阻率在43 Ω·m～77 Ω·m之間,灰?guī)r電阻率在430 Ω·m～2 400 Ω·m之間,砂巖、泥巖電阻率在7 Ω·m～77 Ω·m之間。

根據(jù)WT205剖面反演圖( 圖11),結(jié)合地質(zhì)剖面(圖10)及物性資料,共推測(cè)4個(gè)巖溶塌陷區(qū)、1個(gè)采空區(qū)和3條滑面。巖溶塌陷區(qū)電性特征表現(xiàn)為中等電阻率異常,視電阻率在129 Ω·m～201 Ω·m之間,由于塌陷導(dǎo)致巖石較為松散,含水率較高,相對(duì)灰?guī)r視電阻率降低。該區(qū)域的煤礦為地下開采模式,多數(shù)煤礦未進(jìn)行回填或少數(shù)采取煤矸石進(jìn)行回填,采空區(qū)電性特征表現(xiàn)為高阻異常,視電阻率在312 Ω·m～3 507 Ω·m之間。h1、h2、h3滑面的電性特征都表現(xiàn)為串珠狀低阻異常,視電阻率在3.8 Ω·m～14.3 Ω·m之間,由于滑面含水量較高,相對(duì)滑坡體及滑床視電阻率較低,滑面最大埋深分別為25.2 m、19.9 m、19.1 m。

圖10 地質(zhì)剖面圖Fig.10 Geological profile

圖11 高密度電法反演擬斷面圖Fig.11 Pseudo-section diagram of high-density electrical inversion

根據(jù)高密度電法測(cè)量結(jié)果,在剖面上布設(shè)zk2010、zk2011、zk2012三個(gè)工程鉆孔(圖12)對(duì)高密度電法的推測(cè)結(jié)果進(jìn)行驗(yàn)證。工程鉆孔zk2010在21.09 m以深,可見巖溶填充物,主要巖性為礫巖,成分以玄武巖、灰?guī)r為主,基質(zhì)為灰色粘土、灰黃色砂土等,基底式膠結(jié),與推測(cè)的巖溶塌陷區(qū)基本一致。工程鉆孔zk2011在16.24 m和28.88 m處見滑帶粘土,驗(yàn)證了該位置推測(cè)的兩條滑帶的存在,但滑帶埋深推測(cè)結(jié)果與實(shí)測(cè)結(jié)果有所差異。工程鉆孔zk2012在7.22 m處見滑帶粘土,驗(yàn)證了該位置推測(cè)的一條滑帶的存在,但滑帶埋深推測(cè)結(jié)果與實(shí)測(cè)結(jié)果有所差異。通過(guò)對(duì)比可以發(fā)現(xiàn),隨著滑面埋深越深,推測(cè)結(jié)果與實(shí)測(cè)結(jié)果的誤差比越小,這與數(shù)值模擬結(jié)果一致(表2)。

表2 實(shí)測(cè)數(shù)據(jù)滑面埋深誤差對(duì)比表Tab.2 Comparison of buried depth error of sliding surface in measured data

圖12 工程鉆孔柱狀圖Fig.12 Borehole column diagram

3.5 滑坡形態(tài)特征與形成機(jī)制初步探討

3.5.1 滑坡形態(tài)特征

根據(jù)高密度電法的實(shí)測(cè)結(jié)果,結(jié)合地質(zhì)剖面、工程鉆孔資料,獲取半坡村滑坡工程地質(zhì)剖面圖(圖13)。分析圖圖13可知:①半坡村滑坡在高程、地貌上自上而下分為后緣的圍巖、中部的滑坡體和前緣的松散堆積物;②從物質(zhì)組成上來(lái)看,滑床由茅口組灰?guī)r和宣威組砂巖、泥巖組成,宣威組砂巖、泥巖352°∠18°;滑坡體由茅口組灰?guī)r、峨眉山組玄武巖和第四系碎石土組成;滑面由高至低逐漸平緩,埋深在8 m～26 m。

圖13 工程地質(zhì)圖Fig.13 Engineering geological map

3.5.2 滑坡形成機(jī)制探討

通過(guò)高密度電法剖面獲得的半坡村滑坡的電法數(shù)據(jù),有效揭示了該滑坡的剖面形態(tài)特征、滑體組成與結(jié)構(gòu)及滑面特征,對(duì)滑坡形成機(jī)制分析提供了全面準(zhǔn)確的參考依據(jù)。通過(guò)高密度電法測(cè)量工作,結(jié)合地質(zhì)剖面及鉆孔資料,筆者認(rèn)為半坡村滑坡形成機(jī)制分為兩部分:①滑坡后緣的地裂縫最為密集,并在地裂縫密集區(qū)觀察到有巖溶塌陷區(qū)的存在,結(jié)合通過(guò)高密度電法推測(cè)的地下巖溶塌陷區(qū)的位置,埋深較淺的巖溶塌陷區(qū)位于地裂縫密集區(qū)內(nèi),推測(cè)該區(qū)域的地裂縫由巖溶塌陷與滑坡綜合因素所形成;②滑坡中部和后緣的砂巖、泥巖的傾向與滑坡體的滑動(dòng)方向相反,地層較為穩(wěn)定,不易產(chǎn)生滑坡,滑坡后緣的滑坡剪出口位于高密度電法推測(cè)的采空區(qū)位置,采空區(qū)規(guī)模110×75 m2,推測(cè)該區(qū)域的滑坡主要由于采空區(qū)所形成。

半坡村滑坡推測(cè)存在3個(gè)滑面,形成機(jī)制較為復(fù)雜,并非單一因素影響所致。推測(cè)該滑坡的地質(zhì)歷史演化過(guò)程為:半坡村滑坡后緣的灰?guī)r在水蝕作用下形成多個(gè)不連續(xù)的巖溶塌陷區(qū),導(dǎo)致地表存在多條地裂縫及地表塌陷;又由于灰?guī)r地層較為穩(wěn)定,只形成地層局部變形。后期由于大則勒煤礦開采至半坡村滑坡山腳,在坡腳附近形成采空區(qū),局部發(fā)生塌陷失穩(wěn),致使坡腳失去支撐,引起坡體結(jié)構(gòu)失去平衡,加劇地層形變,形成三階滑坡。

4 結(jié)論

1)通過(guò)數(shù)值模擬及應(yīng)用實(shí)例,驗(yàn)證了高密度電法在多期次滑坡勘察中的有效性,能夠有效識(shí)別多條滑面。

2)半坡村多期次滑坡其滑坡體主要表現(xiàn)為局部不連續(xù)高阻異常電性特征,3條滑面主要表現(xiàn)為串珠狀低阻異常電性特征。

3)通過(guò)高密度電法剖面、地質(zhì)剖面及工程鉆孔資料,構(gòu)建半坡村滑坡工程地質(zhì)剖面,揭示滑坡剖面形態(tài)特征、滑體組成與結(jié)構(gòu)及滑面特征,分析滑坡形成機(jī)制,得出半坡村滑坡是由于巖溶塌陷及后期煤礦開采形成采空區(qū)所致。