隧道含水構造三維瞬變電磁場響應特征及突水災害源預報研究.pdf

1、針對隧道含水構造的瞬變電磁場響應特征和隧道突水災害瞬變電磁探測與預報問題進行了兩方面的研究:第一，采用時域有限差分法實現(xiàn)了針對任意復雜模型瞬變電磁探測的三維正演，完成了隧道含水構造瞬變電磁場響應的三維建模并針對典型含水構造進行了響應特征分析;第二，提出了隧道突水災害瞬變電磁三維多分量陣列式探測方法，并就其解釋技術進行了探索。論文針對上述兩方面的幾個關鍵問題進行了研究。
　　 地面和航空瞬變電磁的三維正反演技術在歐美、澳洲等發(fā)達國

2、家已經(jīng)開始應用于工業(yè)生產(chǎn)，然而，隧道內(nèi)的瞬變電磁探測存在一個復雜的偽均勻全空間環(huán)境，現(xiàn)有的正演技術均不適用于隧道內(nèi)的瞬變電磁探測，其主要原因是采用現(xiàn)有方法對隧道探測施加激勵源存在難度。論文針對這一問題，同時考慮有源和無源區(qū)域的Maxwell方程組，引入激勵源電流密度項到安培定理中，從而實現(xiàn)回線激勵源的注入。理論上，這種施加激勵源的方法適用于任意復雜的三維模型，但其缺點是增加了計算成本。為了提高計算效率，采用了CPU+GPU異構并行計算方

3、案。為了驗證本文給出的新激勵源施加方法的正確性，首先采用均勻半空間模型解析解與三維FDTD解進行對比，結果顯示衰減曲線誤差小于5％，視電阻率誤差小于3.5％;然后采用四類三層模型(A、H、K、Q)的線性數(shù)字濾波解與三維FDTD計算結果進行對比;之后采用均勻半空間包含低阻塊體的模型進行計算并與Wang&Hohmann(1993)文獻中的FDTD解和積分方程解進行對比。以地面驗證的理論為基礎，針對隧道探測進行了專門改進，提出了隧道超前探測瞬

4、變電磁三維正演建模方法，針對直立斷層、傾斜斷層、充水溶洞等隧道開挖過程中可能遇到的典型含水構造進行了正演模擬，考察了含水構造與掌子面的距離、含水構造規(guī)模、含水構造與背景圍巖電阻率差異等參數(shù)對瞬變電磁響應的影響，得到了部分典型含水構造的瞬變電磁響應特征。針對TBM開挖的隧道，計算和分析了TBM機對瞬變電磁超前探測的影響規(guī)律，為隧道瞬變電磁超前地質(zhì)預報用于TBM施工的隧道進行了數(shù)值的探索。分析數(shù)值算例后，新發(fā)現(xiàn)隧道超前地質(zhì)預報中尚未采用的水

5、平感應電動勢響應對含水構造具有明顯的反應，對于判識大規(guī)模含水構造具有重要作用。
　　 現(xiàn)有的隧道瞬變電磁超前探測一般采用中心回線的方式并在掌子面布置1條測線，得到電阻率斷面的二維分布，對含水構造的空間形態(tài)與分布規(guī)律反映較差，本文借鑒地面大回線中心裝置的工作方式，提出了隧道突水災害瞬變電磁三維多分量陣列式探測方法，給出了詳細的裝置布置形式與觀測方法。以一維瞬變電磁基本理論為基礎，研究了小回線半徑下瞬變電磁回線內(nèi)帶偏移距觀測點的水平

6、和垂直磁場頻率域響應中雙重貝塞爾函數(shù)強振蕩型積分的計算方法，以中心點頻率域響應多項式為基函數(shù)，定義多項式形式的帶偏移距觀測點頻率域響應表達式，通過最小二乘擬合獲得多項式系數(shù)。將得到的頻率域響應方程按照反傅里葉變換轉換到時間域，得到瞬變電磁回線內(nèi)帶偏移距觀測點水平和垂直感應電動勢響應的時間域方程，在瞬變電磁場晚期近似條件下，給出了與隧道多點陣列式探測方法配套的水平和垂直感應電動勢響應晚期視電阻率公式，從而使隧道多分量陣列式瞬變電磁探測采集

7、的數(shù)據(jù)能夠轉換為視電阻率圖像進行解釋。以層狀地層的橫向和縱向電阻率合成方法給出了綜合利用水平和垂直響應的電阻率矢量合成方法。采用均勻半空間模型和四類三層模型(A、H、K、Q)進行對比驗證，并采用現(xiàn)場實測數(shù)據(jù)進行分析，結果得到了鉆孔驗證。
　　 上述本文的兩部分研究工作為瞬變電磁正演理論提供了研究基礎，尤其是隧道突水災害瞬變電磁探測的響應特征研究提供了基礎。建立的隧道三維探測方法能夠獲得三維數(shù)據(jù)體，連同正演理論為進一步的三維反演提