輸水隧洞瞬變電磁法勘察的三維可視化分析

渠繼明

（山西省水利水電勘測設計研究院　山西　太原　030024）

?

輸水隧洞瞬變電磁法勘察的三維可視化分析

渠繼明

（山西省水利水電勘測設計研究院山西太原030024）

摘要瞬變電磁法勘察所得到的二維斷面圖，對于復雜的異常體較難圈定其連通性及空間分布狀態，針對這一特定問題，利用Voxler軟件實現瞬變電磁法勘察剖面的三維可視化，使得瞬變電磁法勘察對地下復雜異常體的分析解釋得到明顯提高。

關鍵詞輸水隧洞；采空區；瞬變電磁法；三維可視化

1　序言

山西是一個產煤大省，隨著煤炭資源的大量開采，留下了大量的煤礦采空區。山西又是一個多山的省份，有很多長距離的輸水隧洞經過煤礦采空區，顯然煤礦采空區對輸水隧洞構成巨大威脅，尤其是無圖紙、無規律的私挖亂采形成的采空區對輸水隧洞具有更大的安全隱患。因此煤礦采空區的勘察就顯得格外重要。瞬變電磁法作為煤礦采空區勘察的較有效方法，在煤礦采空區勘察中發揮著重要作用，但瞬變電磁法的二維勘察形成的是獨立的斷面，斷面之間異常的連通性在二維斷面上得不到反映，通過Voxler軟件實現瞬變電磁法勘察剖面的三維可視化，使得二維斷面上的異常在空間展布情況得到直觀的反映。

2　方法與技術

2.1基本原理

瞬變電磁法（簡稱TEM）是一種時間域的電磁感應方法，其探測的基本原理是以不接地的發射線框中加以脈沖電流，在一次場的激勵下，地質體產生渦流，其大小與地質體的導電程度有關。在一次場消失后，該渦流不會立即消失，它是一個逐漸衰減的過程，該衰減過程又產生衰減的二次磁場向地表傳播，由在地表布設的接收線框接收二次磁場，此二次磁場的變化反映了地下地質體的電性特征。瞬變電磁法就是基于上述原理來測量地質體的電導率來達到解決實際問題的目的。

渦旋場在大地中擴散傳播，在傳播過程中，二次電磁場的能量直接在導電介質中傳播而逐漸消耗。由于趨膚效應，低頻部分傳播到深部，而高頻部分主要集中在淺部，且其分布范圍是源下面的局部，其傳播速度和深度為：

式中t為傳播時間；σ為介質電導率；μ0為介質中的磁導率。

依據上述原理，根據觀測的感應電壓值隨時間的特定的變化關系，便可反演得到不同深度的地層的視電阻率值。

2.2地球物理特征

當地層中的煤層未被采動時，地層呈現縱向的成層性和橫向的連續性，反映在視電阻率斷面圖上為視電阻率等值線在橫向上呈近似水平或呈一定傾角的特征，在縱向上呈近似層狀的特征。煤層被采空后，破壞了原有地層的連續性、完整性，反映在視電阻率斷面圖上表現為視電阻率等值線出現扭曲和變形，出現圈閉等現象。如果采空區被水充填，則在視電阻率等值線圖中表現為局部低阻異常圈，如果采空區未被充填，則表現為局部高阻異常圈。

2.3數據預處理與反演

資料處理采用EMRS瞬變電磁處理軟件。處理流程為：數據傳輸→資料預處理→數據轉換→繪制圖件→定性解釋、定量解釋。在資料處理之前須進行資料的預處理，包括：

（1）關斷時間的影響及其校正。理論上發射波形應為階躍波，但在實際工作中由于發射線框與大地間存在電容和線框存在電感，導致發射電流關斷時不為階躍波，而為一個有一定后延時間的關斷波形，因此須將其校正到理想波形狀態。

（2）曲線的圓滑。由于各類干擾的影響，往往使實測V/I曲線尾支出現波動，為此必須對實測曲線進行圓滑處理。

（3）將經過預處理的實測數據反演成與深度對應的視電阻率值，通過視電阻率二維等值線斷面圖或三維可視化分析進行定性和定量的解釋。

2.4三維數據處理與可視化

傳統的瞬變電磁法成果為二維視電阻率等值線圖，它反映的是同一測線上不同深度的視電阻率變化特征，在進行資料解釋時，僅靠二維視電阻率等值線圖很難對異常的空間展布及不同剖面上的異常變化特征作出直觀而準確的判定。Voxler軟件的出現使得物探數據的三維可視化分析成為可能，Voxler主要針對三維數據的顯示和體積的渲染，利用離散點數據以及二維DEM網格化數據和二維圖像進行三維成圖，可以形成矢量圖、等值線圖、切片圖、三維散點圖和體積渲染圖等。使用Voxler輸出的三維模型，可以通過鼠標操作改變視角，也可以被任意的角度切片或顯示。

目前瞬變電磁法數據的采集本質上是二維采集，然后疊加形成三維數據，因此前期的數據處理與反演仍沿用傳統的處理解釋方法。在二維測線的數據處理時選用統一的處理方法和參數，這樣在后期的三維可視化處理時測線與測線之間才會保持一致性，當每條測線反演完畢之后即可進行三維化處理。首先，把每條測線上各個測點對應的樁號或距離轉換成測網對應的平面坐標(x,y)，把測線上各點的深度經過高程換算轉換為高程數據，再把測網內的所有測線的數據合并在一起，便形成一個瞬變電磁法的三維數據體。

三維數據體形成后便可進行網格化，數據經網格化后會形成一個三維模型，然后便可利用Voxler軟件對三維模型進行可視化分析。

3　工程實例

山西省某引水工程是山西省“十二五規劃”大水網建設中一項重要的工程，其中有很大一部分輸水管線穿越煤層或煤礦采空區域，為了確保輸水管線的安全，須對煤礦采空區進行詳細勘察并加以治理。受某引水工程管理局的委托，對該引水工程的某段輸水管線區域進行了瞬變電磁法物探勘察。

根據地質資料，勘察區內地層自上而下為第四系、第三系、石炭系和奧陶系地層。正常情況下，縱向上第四系表層視電阻率較高；第三系、石炭系砂巖等煤系地層一般為表現為低阻，奧陶系地層以灰巖為主，視電阻率值較高。整套地層的電阻率在縱向上呈現為高～低～高的特征反映。該區域采空區埋深約80m～90m，未被水充填，屬于私挖亂采形成的古空區，因此采空區無任何規律性。為了對采空區進行詳細的勘察，在勘察區域內，瞬變電磁點剖面采用小點距、小線距布設，為了基于Voxler輔助條件下對采空區三維可視化分析的方便，在外業測試過程中保持采集參數不變。

3.1勘察成果二維分析：

圖1、圖2、圖3所示為在勘察區域內a2～a2'剖面、a5～a5'剖面、a8～a8'剖面3個TEM剖面的視電阻率等值線圖，從圖中可看出：在a2～a2'剖面的深度80m～90m、水平距離25m～35m出現高阻異常區；在a5～a5'剖面的深度80m～90m、水平距離15m～35m出現高阻異常區，在a8～a8'剖面的深度80m～90m、水平距離16m～23m呈現高阻異常。從選取的3個視電阻率等值線斷面圖可看出：在3個剖面上均存在高阻異常區，結合地質資料推斷為煤礦采空區引起的，但3個剖面之間采空區的空間展布如何，異常之間是否連通，在二維剖面上顯然無法反映出來。

圖1　a2- a2'TEM視電阻率等值線斷面圖

圖2　a5- a5'TEM視電阻率等值線斷面圖

圖3　a8- a8'TEM視電阻率等值線斷面圖

3.2三維可視化分析

在二維圖像上，三個剖面均存在高阻異常區域，但異常僅獨立展現在三個二維剖面上，異常的空間形態及其連通性難以確定，經過三維可視化處理后，三個剖面的異常的連通性在三維等值面圖下的空間展布形態有了清晰的展示，且可進行水平方向不同深度或垂直方向不同位置的切片，便于對采空區綜合分析。因此，經過Voxler軟件三維可視化處理后的資料，突出了異常在三維空間的展布情況，具有較好的實用性。

圖4　三維等值面圖與a2～a2'斷面的切片

圖5　三維等值面圖與a5～a5'斷面的切片

圖6　三維等值面圖與a8～a8'斷面的切片

3.3應用效果

為了對采空區域進行治理并驗證Voxler輔助條件下的三維可視化數據處理的有效性，在物探勘察區域布置鉆孔進行驗證并對采空區范圍布設注漿孔進行注漿處理。

圖7為利用Voxler輸出的三維模型在異常體的中心的水平切片圖（圖中深色部分為推斷的煤礦采空區），圖8為鉆孔驗證與注漿孔確認的采空區平面分布圖，從二者的對比來看，利用Voxler輸出的三維模型劃定的采空區與鉆孔確認的范圍基本吻合，說明了在Voxler輔助條件下的三維可視化數據處理的有效性和實用性。

圖7　三維模型異常體中心水平切片圖

圖8　經鉆孔驗證后的采空區分布圖

4　結語

（1）本文通過具體的工程實例對瞬變電磁數據的二維數據處理與基于Voxler輔助條件下的三維數據分析進行了對比分析；傳統的二維等值線斷面圖僅能反映沿剖面方向的異常體的分布；而經過Voxler軟件處理后的瞬變電磁數據，不僅能夠反映異常體沿剖面方向展布情況，而且能夠直觀地反映剖面之間異常體的空間展布，并可人機交互實現對異常體的全方位觀測。

（2）雖然利用Voxler軟件三維可視化分析可直觀地反映瞬變電磁法勘察的異常體三維空間的展布情況，但處理效果仍依賴于二維資料的精度，并不能從根本上提高資料的解釋精度。陜西水利

參考文獻：

[1]李永銘，楊新明.電磁波井間層析三維可視化分析[J]《工程物探》.2014年第2期

[2]蘆安貴，李法濱.EH4三維數據采集處理與可視化[J]《工程物探》.2015年第2期

（責任編輯：暢妮）

中圖分類號：P631

文獻標識碼：A