馬 圣 敏,夏 峰2，張 建 清

(1.長江勘測規(guī)劃設(shè)計研究有限責(zé)任公司，湖北 武漢 430010； 2.黔南州水利水電勘測設(shè)計研究院，貴州 黔南 588000)

在城市軌道交通巖溶探測時，通常會在預(yù)開挖的隧道兩側(cè)布置多條地質(zhì)鉆孔測線，每條測線按一定間隔布置鉆孔，每對鉆孔之間采用電磁波CT技術(shù)對隧道區(qū)間巖溶發(fā)育情況進行探測(見圖1)。由于鉆孔相對較密，形成的電磁波CT剖面較多，所以在制作相關(guān)圖件時耗費時間較長，以剖面圖中的視吸收系數(shù)異常投影到平面圖最為典型。

圖1 軌道交通電磁波CT剖面布置示意Fig.1 Schematic diagram of electromagnetic wave CT profile layout for subway

目前主要采用人工的方式將縱剖面圖中視吸收系數(shù)異常投影到平面圖的測線上。具體過程是:① 將電磁波CT縱剖面圖視吸收系數(shù)異常復(fù)制到平面圖中(以鉆孔孔口位置為復(fù)制基點)；② 沿平面圖中視吸收系數(shù)異常區(qū)域左、右兩側(cè)分別繪制鉛垂線，并任意繪制一條水平線與兩條鉛垂線相交；③ 從水平線中點引出一條垂線與平面圖上的測線相交，交點為A；④ 計算中點與兩條鉛錘線的距離D1(左)，D2(右);⑤ 以交點A為中心，沿平面圖中測線左側(cè)方向移動距離D1并繪制一個點P1，沿平面圖中測線右側(cè)方向移動距離D2并繪制一個點P2，連接點P1、A、P2，形成一條多段線，如圖2所示。

圖2 電磁波CT剖面圖視吸收系數(shù)異常平面投影常規(guī)方法Fig.2 The usual method of projecting profile abnormal apparent absorption coefficient curve to plane diagram of electromagnetic wave CT

由此可見，人工方式較為繁瑣，效率低，而且人工投影容易犯錯，一旦發(fā)生錯漏，很難找到問題所在，影響到數(shù)據(jù)處理成果的質(zhì)量控制。因此，急需一種新的電磁波CT剖面圖中視吸收系數(shù)異?？焖倨矫嫱队暗姆椒?。

查閱國內(nèi)外相關(guān)文獻，尚無與軌道交通CT探測剖面布置特點相適應(yīng)，且能將剖面圖中視吸收系數(shù)異常投影到平面圖的簡單易行的計算機實現(xiàn)方法。如龐林軍[1]提到了將電磁波剖面CT圖像基巖界面及以下主要視衰減系數(shù)異常進行水平面投影，但未提及具體方法。狄衛(wèi)民等[2-7]只提及平剖面圖計算機輔助繪圖方法，可有效提高生產(chǎn)效率。劉柏根等[8-10]提出了將剖面圖中二維曲線轉(zhuǎn)換為三維曲線的方法，為剖面圖平面投影提供了一種途徑。江龍劍[11]提出了礦山地質(zhì)圖的平剖面圖自動轉(zhuǎn)換系統(tǒng)，先切割礦條段(三維)，然后通過自動連圖的方式形成平剖面圖，這與電磁波CT的剖面圖“切割”地下不良地質(zhì)體類似，不同點在于通過CT剖面圖中視吸收異常形成平面圖異常時事先并不清楚不良地質(zhì)體的三維形態(tài)。楊麗[12]采用空間坐標(biāo)變換實現(xiàn)地質(zhì)平/剖面圖自動成圖，但在軌道交通CT探測剖面梅花形布置條件下，將剖面中視吸收系數(shù)異常區(qū)域坐標(biāo)轉(zhuǎn)換成三維坐標(biāo)再進行投影，算法相對比較復(fù)雜。文獻[13]為實現(xiàn)在CAD中自動投影提供了思路。

本文結(jié)合軌道交通電磁波CT剖面布置特征，以鉆孔位置為參照，形成剖面視吸收系數(shù)異常的準(zhǔn)確、唯一編碼，并制定相應(yīng)的解碼規(guī)則，從而實現(xiàn)大數(shù)據(jù)量電磁波CT剖面視吸收系數(shù)異常的快速平面投影。

1 問題分析

1.1 電磁波CT剖面布置與剖面圖形特征

如圖1所示，武漢軌道交通7號線北延線工程共布置電磁波CT測線5條，即L1-L1′、M1-M1′、R1-R1′、L1-M1、R1-M1測線，每條測線形成的剖面圖如圖3所示。在剖面圖中，上下兩條白色的平行多段線圍成的區(qū)域為CT探測范圍，緊鄰探測范圍頂部的多段線為覆蓋層與基巖的分界線?；鶐r界線以下的虛線為較完整基巖分界線。閉合圓圈為視吸收系數(shù)異常，分別位于MK165孔與MK167孔下。

1.2 實現(xiàn)思路

由于剖面圖是垂直向下的，視吸收系數(shù)異常投影到平面圖上的表現(xiàn)形式為覆蓋在探測線上的一段折線?？梢砸糟@孔為位置參照，為每個剖面視吸收系數(shù)異常制定唯一的位置編碼，如MK165孔下的視吸收系數(shù)異常可表達為“MK165孔左xm至MK165孔左ym”；MK167孔下的視吸收系數(shù)異常可表達為“MK167孔左xm至MK167孔左ym”，表1為示例數(shù)據(jù)。

在進行平面投影時，以平面圖中鉆孔的平面坐標(biāo)為基點，沿視吸收系數(shù)異常所在剖面對應(yīng)的測線方向分別移動一定的距離，形成視吸收系數(shù)異常平面投影折線的兩個端點，也就是位置解碼。如表1中2001號視吸收系數(shù)異常，沿LK048～LK049剖面測線大樁號方向，分別移動0.5 m和4.5 m，形成視吸收系數(shù)異常平面投影折線的兩個端點，將兩個端點連接起來，便完成了剖面視吸收系數(shù)異常的平面投影。

圖3 某軌道交通電磁波CT剖面示意Fig.3 Schematic diagram of electromagnetic wave CT profile of a subway

異常號存在剖面高程范圍/m水平范圍/m2001LK048～LK049-1.5～2.5LK048孔右0.5～4.5m2002LK048～LK049-12.5～-9.0LK048孔右0.3～4.6m2003LK048～LK049～LK050-14.9～-9.2LK049孔左5.8m至LK049孔右2.6m2004LK051～LK052～LK0534.0～7.4LK052孔左4.7m至LK052孔右1.1m2005LK052～LK053～LK0540.1～4.9LK053孔左5.9m至LK053孔右4.3m

該示例是單個視吸收系數(shù)異常位于2個鉆孔之間的情況。對于視吸收系數(shù)異常與1個或多個鉆孔相交、與鉆孔不相交的情況，也可以制定相應(yīng)的編碼、解碼規(guī)則?；谠撍悸?，便可實現(xiàn)剖面視吸收系數(shù)異常投影到平面圖的計算機方法。

2 具體方案

具體而言，基于計算機方法實現(xiàn)大數(shù)據(jù)量電磁波CT剖面圖中視吸收系數(shù)異常在平面圖中的投影分為以下幾個步驟。

2.1 位置編碼

選擇剖面圖中任一視吸收系數(shù)異常，以與該視吸收系數(shù)異常最為鄰近的鉆孔為參照，根據(jù)鉆孔與視吸收系數(shù)異常的位置關(guān)系，遵照一定的編碼規(guī)則，形成每個視吸收系數(shù)異常投影為折線的唯一位置編碼。

(1) 如圖4(a)所示，對于視吸收系數(shù)異常區(qū)域位于兩個鉆孔m、n之間的情況，通過視吸收系數(shù)異常與鉆孔n的空間位置關(guān)系，獲取視吸收系數(shù)異常最左端、最右端與鉆孔n的水平距離參數(shù)a、b，形成位置編碼為“n孔左am至n孔左bm”。

(2) 如圖4(b)所示，對于視吸收系數(shù)異常與一個鉆孔m相交的情況，通過視吸收系數(shù)異常和穿過該視吸收系數(shù)異常的鉆孔m的空間位置關(guān)系，獲取視吸收系數(shù)異常最左端、最右端與鉆孔m的水平距離參數(shù)a、b，形成位置編碼為“m孔左am至m孔右bm”。

(3) 如圖4(c)所示，對于視吸收系數(shù)異常與兩個鉆孔m，n相交的情況，通過視吸收系數(shù)異常和穿過該視吸收系數(shù)異常的兩個鉆孔m、n的空間位置關(guān)系，獲取視吸收系數(shù)異常最左端與鉆孔m、最右端與鉆孔n的水平距離參數(shù)a、b，形成位置編碼為“m孔左am至n孔右bm”。

圖4 電磁波CT剖面圖鉆孔與視吸收系數(shù)異常位置關(guān)系示意Fig.4 Schematic diagram of relationship between borehole and apparent absorption coefficient abnormal in electromagnetic wave CT profile

2.2 解譯與投影

在進行電磁波CT剖面視吸收系數(shù)異常平面圖投影時，逐一取出視吸收系數(shù)異常與鉆孔位置關(guān)系的編碼記錄，根據(jù)記錄中視吸收系數(shù)異常與鉆孔的關(guān)系類型，選擇相應(yīng)的計算機解碼規(guī)則進行解譯，并完成平面圖投影。

(1) 對于視吸收系數(shù)異常位于2個鉆孔之間的情況，解譯編碼記錄，獲取鉆孔孔號n，空間關(guān)系“孔左”，以及距離參數(shù)a、b。從鉆孔平面坐標(biāo)記錄中檢索鉆孔孔號n所對應(yīng)的平面坐標(biāo)。計算平面圖中沿測線方向距鉆孔孔號n左邊a、bm的2個點的平面坐標(biāo)。將其與鉆孔n的平面坐標(biāo)連成折線，完成投影過程。

(2) 對于視吸收系數(shù)異常與一個鉆孔相交的情況，解譯編碼記錄，獲取鉆孔孔號m，空間關(guān)系“孔左、孔右”，以及參數(shù)a、b。從鉆孔平面坐標(biāo)記錄中檢索鉆孔孔號m所對應(yīng)的平面坐標(biāo)。計算平面圖中沿測線方向距鉆孔孔號m左邊a、右邊bm的兩個點的平面坐標(biāo)。將其與鉆孔m、n孔的平面坐標(biāo)一周繪制成點并連成折線，完成投影過程；

(3) 對于視吸收系數(shù)異常與2個鉆孔相交的情況，解譯編碼記錄，獲取鉆孔孔號m、n，空間關(guān)系“孔左、孔右”，以及參數(shù)a、b；從鉆孔平面坐標(biāo)記錄中檢索鉆孔孔號m、n所對應(yīng)的平面坐標(biāo)。計算圖中沿測線方向距鉆孔孔號m孔左邊am、鉆孔孔號n孔右邊bm的兩個點的平面坐標(biāo)。將其與鉆孔m、n孔的平面坐標(biāo)連成折線，完成投影過程。

完成投影過程后，可依據(jù)測線的空間位置關(guān)系，將平面圖中多條具有相關(guān)關(guān)系的視吸收系數(shù)異常投影線的兩端端點逆時針或順時針連接起來，形成視吸收系數(shù)異常平面圖，如圖5所示。

圖5 平面圖中視吸收系數(shù)異常投影線合成示意Fig.5 Sketch diagram of apparent absorption coefficient abnormal synthesis on plan map

3 工程應(yīng)用

該方法在武漢軌道交通7號線北延線工程中應(yīng)用效果良好，以下以該工程為例進行說明。

武漢軌道交通7號線北延線工程起于黃陂前川，先后經(jīng)過黃陂前川新城區(qū)、武漢臨空經(jīng)濟區(qū)、黃陂盤龍城區(qū)、東西湖區(qū)，止于7號線一期工程園博園北站(不含)。工程全長36.47 km，其中高架段長約22.60 km，U型槽段1.08 km，地下段約12.79 km。

電磁波CT探測范圍包括橫臨高架區(qū)段及臨空北路站、巨龍大道站-馬池站和馬池站-園博園北站區(qū)段，共完成電磁波CT剖面1 256對。

通過數(shù)據(jù)分析與整理，形成如表2所示電磁波CT剖面視吸收系數(shù)異常表，表中水平范圍列為視吸收系數(shù)異常的水平位置編碼。采用文中所述方法，在平面圖中形成投影線段如圖6所示。為了區(qū)分不同的視吸收系數(shù)異常，圖中視吸收系數(shù)異常投影的線段采用了不同的顏色與線寬。投影線段正確反映了L1-L1′、R1-R1′、L1-M1、R1-M1測線上如圖4(a)～(b)所示類型的異常。

表2 前川線巨湯區(qū)段電磁波CT剖面視吸收系數(shù)異常結(jié)果(部分?jǐn)?shù)據(jù))Tab.2 Apparent absorption coefficient abnormal of electromagnetic wave CT profile of Qianchuan Subway Project

圖6 前川線平面圖中視吸收系數(shù)異常投影線示意Fig.6 Sketch diagram of apparent absorption coefficient abnormal synthesis on plan map of Qianchuan subway project

4 結(jié) 論

(1) 本文提出了一種新的電磁波CT剖面圖中視吸收系數(shù)異常區(qū)域平面投影方法，可為大數(shù)據(jù)量電磁波CT剖面圖視吸收系數(shù)異?？焖偻队暗狡矫鎴D提供基礎(chǔ)。

該方法以視吸收系數(shù)異常鄰近鉆孔為參照，按照鉆孔與視吸收系數(shù)異常的各種空間位置關(guān)系，遵照一定的編碼規(guī)則，形成視吸收系數(shù)異常投影為折線的唯一位置編碼。在平面投影時，按剖面視吸收系數(shù)異常與鉆孔的位置關(guān)系解碼位置編碼記錄，沿著剖面所對應(yīng)的平面測線方向，投影形成多個平面坐標(biāo)點，連接各坐標(biāo)點形成折線，實現(xiàn)視吸收系數(shù)異常的平面投影。

(2) 剖面圖視吸收系數(shù)異常的快速平面投影只是第一步，如何在平面圖中將具有關(guān)聯(lián)性的平面投影視吸收系數(shù)異常(折線)組合起來，形成閉合的區(qū)域，是下一步需要解決的問題。不良地質(zhì)體形態(tài)錯綜復(fù)雜，視吸收系數(shù)異常平面投影線的合成需要制定相應(yīng)的準(zhǔn)則，才能為快速電磁波CT視吸收系數(shù)異常平面圖制作打下基礎(chǔ)。