郭聰 徐濤 陳江平

摘要：旭龍水電站前期工程勘察中面臨壩庫(kù)區(qū)高陡邊坡、滑坡體及超厚覆蓋層等復(fù)雜地質(zhì)條件，準(zhǔn)確查明該區(qū)域的隱蔽構(gòu)造、裂縫、斷層及覆蓋層厚度十分關(guān)鍵。引入新型大地電磁測(cè)深探測(cè)系統(tǒng)，并基于優(yōu)化改進(jìn)二維快速松弛反演（RRI）算法，得到了收斂快、曲線平滑的大地電磁測(cè)深剖面，提取了豐富的深部地層信息。結(jié)果表明：采用低頻磁棒、改進(jìn)的RRI反演算法，有效數(shù)據(jù)頻帶延伸到1～0.1 Hz，有效探測(cè)深度近1 000 m。該技術(shù)在超厚覆蓋層及深部構(gòu)造探測(cè)中獲得了較好的應(yīng)用效果，克服了旭龍水電站前期勘察技術(shù)難點(diǎn)，可為新型大地電磁系統(tǒng)在類似工程的應(yīng)用提供參考。

關(guān)鍵詞：地球物理勘察； 大地電磁法； 二維RRI反演； 隱蔽構(gòu)造； 旭龍水電站

中圖法分類號(hào)：P631

文獻(xiàn)標(biāo)志碼：A

DOI：10.15974/j.cnki.slsdkb.2024.05.007

文章編號(hào)：1006-0081（2024）05-0036-07

0 引言

旭龍水電站是金沙江上游河段“一庫(kù)十三級(jí)”開(kāi)發(fā)方案的第12級(jí)，也是國(guó)家西電東送骨干電源點(diǎn)之一，2022年6月獲得國(guó)家發(fā)展和改革委員會(huì)核準(zhǔn)批復(fù)，成為“十四五”期間核準(zhǔn)開(kāi)工建設(shè)的國(guó)家重大工程之一。旭龍水電站建成投產(chǎn)后，預(yù)計(jì)年發(fā)電量約105億kW·h，節(jié)約標(biāo)準(zhǔn)煤315萬(wàn)t，將有效緩解區(qū)域電力供應(yīng)緊張，對(duì)地區(qū)經(jīng)濟(jì)社會(huì)發(fā)展意義重大。

早在2004年，旭龍水電站預(yù)可研階段勘察工作已經(jīng)啟動(dòng)。前期資料表明，水電站所在的金沙江上游河段地質(zhì)條件復(fù)雜，隱蔽構(gòu)造多，區(qū)域裂縫及斷層發(fā)育，覆蓋層較厚，庫(kù)區(qū)邊坡穩(wěn)定性差，地質(zhì)勘察難度較大。多年來(lái)，國(guó)內(nèi)專家學(xué)者對(duì)此開(kāi)展大量研究工作，楊鵬等［1］對(duì)旭龍水電站庫(kù)區(qū)格亞頂堆積體變形及穩(wěn)定性進(jìn)行分析，劉沖平等［2］對(duì)旭龍水電站庫(kù)壩區(qū)邊坡危巖體進(jìn)行評(píng)價(jià)。

為進(jìn)一步調(diào)查分析旭龍水電站庫(kù)壩區(qū)隱蔽構(gòu)造特征及深部不良地質(zhì)體，國(guó)內(nèi)學(xué)者充分借鑒類似工程的物探勘察方法開(kāi)展研究工作，郭一兵等［3］研究高密度電法在特大型滑坡勘察中的應(yīng)用，趙驚濤等［4］采用地震繞射波探測(cè)方法對(duì)煤礦隱蔽致災(zāi)地質(zhì)體進(jìn)行分析，張建清等［5］也在同類水電站建設(shè)中開(kāi)展過(guò)綜合物探方法研究工作，都取得了較好的效果，給旭龍水電站隱蔽構(gòu)造及深部地質(zhì)構(gòu)造探測(cè)提供了思路。采用高密度電法能較好地查明測(cè)區(qū)淺層覆蓋層厚度及斷層分布情況，但無(wú)法探測(cè)超厚覆蓋層及深部斷裂信息；采用瞬變電磁、地震法等能滿足探測(cè)深度要求，但容易受到環(huán)境條件影響，不易達(dá)到勘察精度要求。因此，復(fù)雜地質(zhì)條件及不利環(huán)境因素成為了旭龍水電站隱蔽構(gòu)造探測(cè)的難點(diǎn)。

在前人研究基礎(chǔ)上，本文針對(duì)旭龍水電站隱蔽構(gòu)造及超厚覆蓋層開(kāi)展大地電磁法研究工作。大地電磁法探測(cè)深度大，不易受到環(huán)境干擾。音頻大地電磁法（AMT）或大地電磁法（MT）起步較早，數(shù)據(jù)處理及解釋較為成熟，反演方法種類繁多。近年來(lái)，隨著技術(shù)發(fā)展，高精度大地電磁法反演方法也有所突破。目前運(yùn)用較成熟的有奧可姆法（OCCAM）反演及簡(jiǎn)化基奧可姆法（REBOCC）反演 ［6］，這些反演方法能反映出地層電性結(jié)構(gòu)特征。Smith等［7］于1999年提出的SBI反演方法，改變了邊界位置的收斂度，同年開(kāi)發(fā)的快速松弛法算法（RRI）成功引入到國(guó)內(nèi)，成為了國(guó)內(nèi)MT工作者經(jīng)常使用的反演方法，并內(nèi)嵌到各種軟件中。該反演方法通過(guò)建立數(shù)據(jù)擾動(dòng)和模型參數(shù)擾動(dòng)之間的線性積分方程，使求解過(guò)程變得迅速，且更容易收斂。

本文從AMT的方法原理入手，介紹一種新型大地電磁測(cè)深系統(tǒng)（GEM-3D）。對(duì)儀器裝備技術(shù)總結(jié)、反演方法改進(jìn)、數(shù)據(jù)處理與解釋、實(shí)例分析等開(kāi)展了研究，并在旭龍水電站前期物探勘察中取得了良好的應(yīng)用效果。

1 音頻大地電磁（AMT）方法原理

音頻大地電磁（AMT）探測(cè)基本原理：當(dāng)天然交變電磁場(chǎng)入射大地，在地下以波的形式傳播時(shí)，由于電磁感應(yīng)，接收器會(huì)接收到地面電磁場(chǎng)信號(hào)，這些信號(hào)會(huì)包含地下介質(zhì)的電阻率分布信息。由于不同頻率的電磁場(chǎng)信號(hào)具有不同穿透深度，因此可以通過(guò)研究地表采集的大地電磁數(shù)據(jù)來(lái)反演出地下不同深度介質(zhì)電阻率分布的信息［8］。

郭聰 等大地電磁探測(cè)技術(shù)在旭龍水電站隱蔽構(gòu)造勘察中的應(yīng)用

通過(guò)傅立葉變換將時(shí)間域的電磁信號(hào)變成頻譜信號(hào)，進(jìn)而計(jì)算卡尼亞電阻率Z［9］：

2 二維快速松弛反演（RRI）

2.1 RRI反演算法改進(jìn)

快速松弛反演（Rapid Relax Inversion，RRI）是通過(guò)求解與一維相近的反演問(wèn)題來(lái)計(jì)算測(cè)點(diǎn)的視電阻率矩陣，從原理上可以把二維問(wèn)題轉(zhuǎn)成一維問(wèn)題［11］?；舅悸啡缦?。

圖1頻率與電阻率交會(huì)量版圖

式中：Q（yi）為各測(cè)點(diǎn)的拉普拉斯（Laplace）范數(shù)方程；f（z）為控制標(biāo)度尺長(zhǎng)度，來(lái)度量不同深度的構(gòu)造；z0表示模型最高頻率情況下的趨膚深度；g（z）為控制水平方向構(gòu)造的懲罰因子［12-13］。

在實(shí)際反演中，為了得到較好的反演效果和有效深度，需要對(duì)RRI反演方法進(jìn)行改進(jìn)。具體做法如下：① 先對(duì)數(shù)據(jù)進(jìn)行預(yù)處理，剔除異常畸點(diǎn)；② 對(duì)數(shù)據(jù)初次反演，并進(jìn)行平滑處理，得到初始模型；③ 將初次反演結(jié)果作為初始模型進(jìn)行RRI反演，同時(shí)控制好f（z）參數(shù)。一般情況下，設(shè)置合適的首層厚度，選取適當(dāng)?shù)脑鲩L(zhǎng)因子和水平構(gòu)造懲罰因子g（z），使反演邊界達(dá)到趨膚深度。經(jīng)過(guò)改進(jìn)后的RRI算法具有收斂速度快、反演剖面平滑、深部信息豐富等特點(diǎn)，適用于隱蔽構(gòu)造的反演成像［14-15］。

2.2 RRI反演應(yīng)用

為驗(yàn)證RRI反演算法可靠性，選取旭龍水電站壩址區(qū)的一段測(cè)線進(jìn)行反演分析。試驗(yàn)測(cè)線長(zhǎng)1 200 m，點(diǎn)距50 m，按照2.1節(jié)中反演步驟，選取首層厚度為13 m，增長(zhǎng)因子1.1，為了達(dá)到一定的趨膚深度，層數(shù)選取23層。最終得到反演結(jié)果如圖2顯示。

從圖2中看出，二維反演后的地層構(gòu)造電性特征較為明顯，表層覆蓋層低阻特征明顯，盲區(qū)也相對(duì)小。從深部看，在樁號(hào)800～1 000 m之間，深度范圍2 000～2 400 m見(jiàn)一處低阻異常，兩側(cè)電阻差異明顯，結(jié)合地質(zhì)資料和鉆孔資料推斷為巖體破碎?？傮w來(lái)看，覆蓋層與基巖分界面清晰，邊界輪廓無(wú)毛刺現(xiàn)象，深部地質(zhì)構(gòu)造明顯，低阻異常區(qū)反應(yīng)突出，表明改進(jìn)后的反演算法對(duì)隱蔽性強(qiáng)的地質(zhì)構(gòu)造反演效果好。

3 旭龍水電站實(shí)例應(yīng)用

3.1 工程地質(zhì)概況

旭龍水電站位于云南省德欽縣與四川省得榮縣交界的金沙江干流上游河段，地跨川、滇、藏兩省一區(qū)，屬于青藏高原強(qiáng)烈隆起區(qū)。水電站選址區(qū)地處橫斷山山地，呈明顯臺(tái)階狀下降趨勢(shì)，屬于峽谷地貌特征，總體地勢(shì)為西北高東南低。選址區(qū)河谷為深切型的下切峽谷地形。

選址區(qū)建基巖體以花崗巖和斜長(zhǎng)角閃片巖為主，兩岸斷裂構(gòu)造較發(fā)育，但構(gòu)造組合相對(duì)簡(jiǎn)單，斷層規(guī)模小，查明這些斷裂構(gòu)造和斷層是地質(zhì)勘察的關(guān)鍵。

庫(kù)區(qū)沿線山坡陡峻、為高山峽谷。兩岸山體呈臺(tái)階下降，覆蓋層以第四系沖洪積、崩坡積為主，覆蓋層厚度大，基巖基本以大理巖、花崗巖為主。覆蓋層與基巖電性特征差異較大，為大地電磁法勘探提供了必要的地球物理?xiàng)l件。

3.2 新型大地電磁系統(tǒng)介紹

為采集數(shù)據(jù)質(zhì)量好的大地電磁原始資料。本文介紹一種新型大地電磁Geode-EM3D系統(tǒng)（簡(jiǎn)稱GEM-3D系統(tǒng)），該系統(tǒng)采用分布式網(wǎng)格化設(shè)計(jì)，配置多個(gè)接收機(jī)（Station），可以實(shí)現(xiàn)3D采集，單臺(tái)接收機(jī)有6個(gè)通道，可以配置3個(gè)磁通道或3個(gè)電通道。經(jīng)測(cè)試，單站可以同時(shí)采集5個(gè)測(cè)點(diǎn)（標(biāo)量），效率大大增加。在進(jìn)行AMT采集過(guò)程時(shí)，通過(guò)單站與微機(jī)進(jìn)行網(wǎng)絡(luò)信息鏈接（Network Information Build）NIB連接，可實(shí)現(xiàn)數(shù)據(jù)采集動(dòng)態(tài)監(jiān)視，提高數(shù)據(jù)采集質(zhì)量。在野外工作布置中，通常采用單站式采集模式。圖3是單站式GEM-3D系統(tǒng)的工作布置示意。

為了更好地采集深部地層信息，在實(shí)際探測(cè)中采用GEM-3D系統(tǒng)的G20K低頻磁棒，它的頻率范圍為0.1～20 kHz，不僅探測(cè)深，且噪聲干擾小。在實(shí)際測(cè)試中，對(duì)比分析了兩種不同頻率范圍的磁棒，圖4為采用不同頻率磁棒采集的兩組標(biāo)量數(shù)據(jù)。從圖4維Bostick反演曲線看出，探測(cè)深度達(dá)到數(shù)千米。此外，低頻采集對(duì)“死頻帶”有一定的壓制作用。

3.3 覆蓋層勘察應(yīng)用

在旭龍水電站庫(kù)區(qū)移民安置點(diǎn)和金沙江兩岸堆積體中選擇2處測(cè)點(diǎn)開(kāi)展覆蓋層厚度和滑坡體地質(zhì)狀況物探勘察工作。旭龍水電站移民安置點(diǎn)A點(diǎn)位于樞紐區(qū)下游4～5 km，海拔高程2 300～3 000 m。覆蓋層為冰水堆積碎石土層，主要由碎塊石、礫石及大量粉土組成。在該安置點(diǎn)上布置AMT測(cè)線，點(diǎn)距為50 m。經(jīng)過(guò)RRI反演后，最終得到視電阻率剖面如圖5所示。

從圖5看出，剖面整體成層性較好，淺層視電阻率低，中深部視電阻率有升高趨勢(shì)，橫向上看電性特征變化均勻，近地表電性差異明顯。覆蓋層與基巖界面清晰，覆蓋層厚度從小樁號(hào)到大樁號(hào)逐漸增大，符合該地區(qū)深切峽谷地質(zhì)特征。厚度變化范圍130～300 m；在樁號(hào)1 260 m左右有一鉆孔ZK2孔，鉆孔揭露在239 m尚未鉆穿覆蓋層，推斷該處的覆蓋層深度在250～300 m之間，反演結(jié)果與鉆孔結(jié)果相吻合。

在旭龍水電站庫(kù)區(qū)金沙江兩岸分布很多滑坡體和堆積體，如格亞頂1號(hào)堆積體、格亞頂2號(hào)滑坡體、木擁橋堆積體、莫丁1號(hào)堆積體、拿榮1號(hào)堆積體、吉布卡堆積體等。這些滑坡體、堆積體地形陡峭，巖體疏松，經(jīng)常發(fā)生小規(guī)模的崩塌、掉塊現(xiàn)象。查明其覆蓋層厚度和基巖面情況，對(duì)庫(kù)區(qū)高邊坡防護(hù)有重要意義。

在一堆積體上布置沿江方向測(cè)線M1，點(diǎn)距為25 m。通過(guò)數(shù)據(jù)采集、處理、反演得到M1測(cè)線AMT視電阻率反演剖面如圖6所示。

從圖6看出，剖面視電阻率縱向呈由低到高的漸變趨勢(shì)，淺部局部電性分布不均勻，地層電性差異較明顯，覆蓋層與基巖界面清晰，厚度變化范圍30～80 m。覆蓋層以砂卵石、碎塊石為主。從橫向上看，基巖埋深總體呈現(xiàn)均勻變化，推測(cè)該堆積體由于山體巖體風(fēng)化、巖屑滑坡逐漸形成厚度不均的覆蓋層，具有不穩(wěn)定性，會(huì)隨邊坡發(fā)生走滑和落石。

3.4 隱蔽構(gòu)造勘察應(yīng)用

壩址區(qū)兩岸雖然巖性相對(duì)單一，但區(qū)域斷裂、斷層發(fā)育，全區(qū)分布著70余條大小斷層，絕大多數(shù)斷層走向與金沙江岸坡交角較大。區(qū)域處于地質(zhì)活躍帶，受構(gòu)造運(yùn)動(dòng)影響造成褶皺、巖體局部破裂，如不及時(shí)查明并處置這些隱蔽構(gòu)造及不良地質(zhì)體，將對(duì)大壩施工造成安全隱患。以過(guò)壩址的交通隧洞線路G2探測(cè)成果作為實(shí)例，分析隧洞高程上下巖體的完整性，查明隧洞線路上地層地質(zhì)構(gòu)造發(fā)育特征，為進(jìn)一步勘察設(shè)計(jì)提供依據(jù)。G2測(cè)線點(diǎn)距50 m，經(jīng)數(shù)據(jù)采集、處理、反演得到G2測(cè)線AMT視電阻率反演剖面如圖7所示。

從圖7看出，剖面總體視電阻率較高，橫向電阻率變化較大。在樁號(hào)1 500 m左右，視電阻率等值線不連續(xù)，存在電阻率差異，推測(cè)為一斷層F，傾向小樁號(hào)方向，傾角約85°，為一逆斷層。在樁號(hào)400～1 200 m、高程2 100 m以下，有一處低阻異常帶，引起上部地層電阻率在橫向上變化較大，推斷為巖體破裂帶。經(jīng)鉆孔驗(yàn)證，該區(qū)域基巖為中、弱風(fēng)化帶，破碎充填物多為鐵錳質(zhì)礦物，導(dǎo)致電阻率偏低。此外，在樁號(hào)1 600～2 600 m地表面視電阻率較高，電阻率值與基巖相似，推斷該區(qū)段除有局部較薄的崩坡積層揭露外，其余為基巖出露。

從G2測(cè)線整體來(lái)看，反演剖面不僅反映出該測(cè)線宏觀地質(zhì)構(gòu)造，而且對(duì)小尺度的巖性變化、裂隙等有很好的探測(cè)效果。

4 結(jié)論

（1） 大地電磁探測(cè)技術(shù)對(duì)旭龍水電站測(cè)區(qū)超厚覆蓋層和斷裂帶、堆積體、滑坡體等不良地質(zhì)體勘察有良好的應(yīng)用效果，對(duì)進(jìn)一步研究三維觀測(cè)系統(tǒng)布置、提升深部勘探成果的準(zhǔn)確性具有重要意義。

（2） 隱蔽構(gòu)造大地電磁探測(cè)方法在旭龍水電站勘察的運(yùn)用，彌補(bǔ)了常規(guī)物探方法的不足，解決了旭龍水電站隱蔽構(gòu)造、超厚覆蓋層、邊坡堆石體等勘察難點(diǎn)，為水電站后續(xù)建設(shè)提供了技術(shù)支撐。

（3） 本文對(duì)新型大地電磁技術(shù)裝備應(yīng)用進(jìn)行總結(jié)，優(yōu)化改進(jìn)大地電磁探測(cè)技術(shù)，拓寬了物探思路，對(duì)類似工程物探技術(shù)應(yīng)用具有參考價(jià)值。

參考文獻(xiàn)：

［1］楊鵬，鄭光，黃金成，等.旭龍水電站庫(kù)區(qū)格亞頂堆積體變形及穩(wěn)定性分析［J］.人民長(zhǎng)江，2023，54（9）：128-135.

［2］劉沖平，王吉亮，黃孝泉，等.環(huán)境邊坡危巖體穩(wěn)定性定性評(píng)價(jià)研究［J］.水利水電快報(bào)，2023，44（1）：43-47.

［3］郭一兵，姜鑫，郭富赟，等.高密度電法在特大型滑坡應(yīng)急勘查中的應(yīng)用［J］.甘肅地質(zhì)，2023，32（3）：82-87，93.

［4］趙驚濤，彭蘇萍，陳宗南，等.煤礦隱蔽致災(zāi)地質(zhì)體地震繞射波探測(cè)方法［J］.礦業(yè)科學(xué)學(xué)報(bào)，2022，7（1）：1-8.

［5］張建清，陳敏，蔡加興，等.綜合物探檢測(cè)技術(shù)在烏東德水電站建設(shè)中的應(yīng)用［J］.人民長(zhǎng)江，2014，45（20）：59-63.

［6］CONSTABLE S C，PARKER R L，CONSTABLE C G .Occam′s inversion，a practical algorithm for generating smooth models from electromagnetic sounding data［J］.Geophysics，1987，52（3）：289-300.

［7］SMITH J T，BOOKER J R.Rapid inversion of two- and three-dimensional magnetotelluric data［J］.Journal of Geophysical Research Solid Earth，1991，96（B3）：3905-3922.

［8］尹劍，張建清.基于音頻大地電磁測(cè)深法數(shù)據(jù)的不良地質(zhì)體探測(cè)研究——以西藏某水利工程為例［J］.水資源與水工程學(xué)報(bào) 2020，31（3）：180-185.

［9］湯井田，何繼善.可控源音頻大地電磁法及其應(yīng)用［M］.長(zhǎng)沙：中南大學(xué)出版社，2005.

［10］陳樂(lè)壽，王光鍔.大地電磁測(cè)深法［M］.北京：地震出版社，1990.

［11］譚捍東，余欽范，BOOKER J，等.大地電磁法三維快速松弛反演［J］.地球物理學(xué)報(bào)，2003，46（6）：850-854.

［12］葉益信，鄧居智，楊海燕，等.SBI與RRI結(jié)合的MT二維反演方法研究及應(yīng)用［J］.東華理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，36（2）：204-209.

［13］高才坤，湯井田，王燁，等.基于RRI反演的高頻大地電磁測(cè)深在深邊部礦產(chǎn)勘探中的試驗(yàn)研究［J］.地球物理學(xué)進(jìn)展，2009，24（1）：309-314.

［14］劉戰(zhàn).音頻大地電磁法在探測(cè)斷層發(fā)育區(qū)中的應(yīng)用［J］.工程地球物理學(xué)報(bào)，2019，16（5）：730-736.

［15］陳玉玲，韓凱，陳貽祥，等.可控源音頻大地電磁法在巖溶塌陷勘察中的應(yīng)用［J］.地球物理學(xué)進(jìn)展，2015，30（6）：2616-2622.

編輯：李慧

Study on application of AMT in geophysical exploration of hidden geological structure in Xulong Hydropower Station

GUO Cong1，2，XU Tao1，2，CHEN Jiangping1，2

（1.Changjiang Geophysical Exploration & Testing Co.，Ltd.，Wuhan 430010，China；2.National Dam Safety Engineering Technology Research Center，Wuhan 430010，China）

Abstract：

Xulong Hydropower Station is faced with complex geological conditions such as high and steep slopes，landslides and ultra-thick overburden in the dam and reservoir area，and it is very critical to accurately identify the hidden structures，cracks，faults and thickness of the overburden in the area.A new Audio-frequency Magnetotelluric method was introduced.By improving the algorithm of two-dimensional Rapid Relax Inversion（RRI），the AMT apparent resistivity profile with fast convergence and smooth curve was obtained，which contained information of deep formation.The results showed that the effective data band extended from 1.0 to 0.1 Hz and the effective detection depth was nearly 1 000 m with low frequency coil data acquisition and improved RRI inversion.The method used in the super-thick cover layer and deep hidden structure geophysical exploration was effective，which solved the difficulties in the preliminary exploration and provided a reference for geophysical exploration of similar hydropower stations.

Key words：

geophysical exploration； AMT method； 2D RRI Inversion； hidden geological structure； Xulong Hydropower Station

收稿日期：2023-11-08

作者簡(jiǎn)介：郭聰，男，工程師，碩士，主要從事水利工程物探技術(shù)研究與應(yīng)用工作。E-mail：xjqxgc@qq.com