綜合物探方法在城市地質調查中的應用研究：以湖南省常德市鼎城區隱伏基巖探測為例

賀轉利，何 禹

（1. 湖南省水文地質環境地質調查檢測所，湖南 長沙 410014；2. 湖南省地質調查所，湖南 長沙 410114）

0 引言

自2017年自然資源部發布《關于加強城市地質調查工作的通知》以來[1]，從國家層面的自然資源部和中國地質調查局到省級層面的自然資源廳以及地質行業管理部門部署和開展了一系列城市地質調查相關工作。2017年底自然資源部發布了城市地質調查規范[2]，涉及到的地球物理探測工作主要有以下幾個方面：城市基礎地質條件探測、城市自然資源調查與評價、城市地質環境調查與評價、城市地質災害調查等內容。

城市地質調查是城市開發的先行工作，是城市空間有效綜合利用規劃的基礎。自1999年開展國土資源大調查以來[3]，城市地質工作從粗到細，呈現出快速發展趨勢，由國土資源部、中國地質調查局牽頭先后試點完成了全國主要城市（群）環境地質調查評價，以及與上海、北京、天津、廣州、南京、杭州等市政府合作，開展三維城市地質調查，系統建立了城市地下三維結構，建立了三維可視化城市地質信息管理決策平臺和面向公眾的城市地質信息服務系統。2004—2012年，國土資源部和中國地質調查局統一規劃和部署，完成了全國主要城市環境地質調查，初步查明了滑坡崩塌泥石流、地面沉降、水土污染、活動斷裂、礦山地質環境問題等各類城市環境地質問題，摸清了地下水、地熱、礦泉水、地質景觀等地質資源狀況。2009年開始，在國土資源部的統一部署下，采用部、省、市多方合作模式，完成了福州、廈門、泉州、蘇州、鎮江、嘉興、合肥、石家莊、唐山、秦皇島、濟南等28個城市地質調查工作。2010年以來，為服務國家區域戰略和主體功能區劃的需求，組織開展了京津冀、長三角、珠三角、海峽西岸、北部灣、長江中游、關中、中原、成渝等重點城市群綜合地質調查工作。2008年7月以來，湖南省地質調查院承擔了長株潭城市群地質環境綜合調查與區劃、長株潭城市群水工環地質調查等項目，目前，正在實施的有常德市城市地質環境調查、益陽市城市地質環境調查和岳陽市城市地質調查項目[3]。

隨著城市地質調查工作逐漸深入，地球物理探測技術也取得了一定科技成果。吳儀芳等[4]以銅陵市長江路大面積出現的地裂縫及塌陷為例，利用電法、磁法、人工地震等方法進行多方法探測試驗，得出采用淺層地震勘探為主的方法是有效探測手段；陳相府等[5]以廣州市城市地質調查為例，指出利用地震橫波勘探結合波速測井進行第四紀松散層層序劃分及厚度探測是可行的；李曉斌等[6]對常見的重、磁、電、震、放射性勘探實例進行論證，證實了物探方法在城市地質調查中具有良好的應用前景；楊歧焱等[7]采用夯源為人工震源的淺層地震勘探方法對日喀則城市活斷層進行了評價；韓術合等[8]利用三維地震、可控源音頻大地電磁法和瞬變電磁法對赤峰市煤礦采空區進行了綜合探測，效果良好；陳實等[9-10]利用高密度電阻率法和天然源面波微動勘探法在烏魯木齊進行了調查應用研究，對后期的城市地下空間資源探測工作具有指導意義；蔣波等[11]利用可控源音頻大地電磁、高密度電阻率、瞬變電磁和地質雷達等方法技術對徐州城市主要斷裂帶分布位置和古河道等進行了探測；王亞輝等[12]對大西安西咸新區利用綜合物探方法對城市地下空間進行了探測和建模；周磊等[13]在湖南郴州市某城區進行了利用等值反磁通瞬變電磁法對城鎮有限場地條件下物探方法找水試驗研究。以上成果都不同程度推動了城市區物探方法應用的進步，本文針對湖南特有的湖湘盆地城市電性的低阻覆蓋層下伏基巖結構進行探測試驗，以期為城市地質調查工作中的物探應用工作提供經驗。

1 地質概況

研究區位于洞庭盆地西南緣[14]，地理位置位于常德市鼎城區工業開發區，試驗剖面東部為太陽山隆起，西部為武陵隆起，南部為沅江，研究區內水系較為發育，出露地層較為簡單。

1.1 地層

根據區域地質資料，本次工作研究區地表出露的地層主要為前第四紀基巖（包括青白口系Qb～古近系E）、第四系下更新統常德組、第四系中更新統馬王堆組、第四系全新統（圖1）。主要地層情況如下所述。

圖1 研究區地質及物探試驗剖面位置圖Fig. 1 Location map of geological and geophysical test sections in the study area

1）前第四紀基巖（Qb-E）。因常德市區及外圍大面積被第四系覆蓋，本次研究把元古界、古生界、中生界和新生界的古近系統一為前第四系基巖。其中，青白口系為研究區內出露的最老地層，包括冷家溪群、板溪群，集中分布于研究區北側常德鼎城區太陽山一帶。南華系和震旦系僅分布于太陽山呈帶狀展布。寒武系為研究區內分布最廣泛的地層，分布于研究區北西側，集中于石板灘鎮以東的連臺-槐樹崗-鐵家咀-觀獅堰、灌溪鎮以西的馬頭寺-馮家橋-樟樹灣、太平-龔家灣-李家-南堰，陬市鎮劉家坡一帶亦有少量分布。奧陶系分布于研究區內北西角，主要呈北西向帶狀分布于石板灘鎮以西的金雞山-趙家灣、灌溪鎮周家-寨子山-長嶺崗水庫。白堊系集中于研究區內西側灌溪鎮北西的桫木塘-福家山、吳家灣-寶鉦灣一帶，局部第四系分布區也有零星露底。古近系集中于研究區北側常德柳葉湖北北東側黃家沖-十字灣-肖伍鋪一帶。

2）第四系下更新統常德組（Qp1cd）。主要分布于丘崗-崗狀平原地貌區，呈面狀或帶狀展布，具體分布于河伏-灌溪鎮和丁家港-斗姆湖-草坪地區。地層出露標高一般60～120 m，沉積物一般厚10～30 m，最厚可達70 m。巖性主要為礫石層夾少量透鏡狀砂層、黏土層等。由下往上沉積物粗細變化頻繁，內部沖刷充填、侵蝕切割等現象十分普遍，為以辮狀河沉積為主的沖洪積扇狀體，內部偶夾少量泥石流沉積。

3）第四系中更新統馬王堆組（Qp2mw）。在丘陵地貌，崗狀平原、波狀平原等地區都有分布，出露標高40～150 m，地層厚度橫向變化較大，厚度3～20 m，一般在地貌發生截變處厚度變大。巖性主要為棕紅色、棕黃色網紋狀黏土、含（弱）網紋狀含鐵錳質黏土，底部常見網紋狀礫質黏土。其與下部地層（常德組、白砂井組等）為不整合接觸關系，常呈小角度披覆或大角度斜切下部砂礫石層。

4）第四系全新統（Qh）。發育于地勢相對較低的現代河流沖溝中，地層分布較為局限，厚度較小，厚2～9 m；多具河流二元結構，下部為礫石層，上部為粉砂質黏土層。

研究區內深部有奧陶系白水溪組和寒武系探溪組、污泥塘組、牛蹄塘組等地層，主要地層情況如下所述。

1）奧陶系白水溪組（O1bs）。本組巖性簡單，以灰黃-淺灰綠色薄-中層狀粉砂質泥巖為主，局部夾含粉砂質泥巖、泥巖、條帶狀泥巖，底部夾泥質灰巖透鏡體，中部夾粉砂質頁巖。體，中部夾粉砂質頁巖。

2）寒武系探溪組（∈3-4t）。以砂屑灰巖、粉屑灰巖、泥晶含泥質灰巖、微晶灰質白云巖、板狀頁巖等為主要巖性組合，整合于下伏地層污泥塘組之上。

3）寒武系污泥塘組（∈2-3w）。該組以中層狀、薄層狀條帶泥晶泥質灰巖、條帶狀粉晶灰巖、粉屑灰巖夾鈣質泥巖、泥巖和鈣質頁巖為特征，與下伏牛蹄塘組呈整合接觸，在灌溪鎮廖家灣偶見星點露頭。

4）寒武系牛蹄塘組（∈1-2n）。該組總體為一套深灰色、灰黑色（含）炭質頁巖、條帶狀含炭質頁巖、炭質泥巖，薄層-中厚層狀含炭質泥巖、含炭質硅質泥巖、含硅質炭質泥巖，夾粉砂質泥巖，局部夾鈣質頁巖和透鏡狀、似層狀硅質巖。

1.2 地質構造

研究區前第四紀斷裂構造主要分布于北西部太陽山-馬頭山一帶，主要生成于印支期、燕山期和古近紀。鑒于前第四紀地層分布在研究區北部，所有前第四紀斷裂都分布在研究區北部，其中，北北東向斷裂主要分布在太陽山凸起的兩側，具有規模大，構造行跡醒目，多期活動的特征，為不同時期區域應力場或局限應力場作用的產物；北東東-北東向斷裂規模較小，被北北東斷裂切割、限制較明顯。試驗剖面附近的主要以研究區內F3斷裂為代表（圖1）。河伏-石板灘斷裂（F3斷裂）為太陽山西側規模較大的斷裂，大體沿漸河延伸于河伏至石板灘一帶，研究區內延伸長度18.5 km，走向NE10°～15°，斷面向南東傾斜，傾角50°～70°，呈正斷裂兼具走滑特征。斷裂在地貌上表現為平直的深溝或谷地，兩側地形及地層存在明顯差異，斷裂東側為太陽山西麓殘坡積物，水系不發育，斷裂西側為丘崗地貌，為早更新世砂礫石層組成的沖積扇體，發育垂直斷裂走向的羽狀水系。

2 地球物理方法技術及地球物理特征

2.1 地球物理方法技術

本次應用試驗研究工作在前期調研的基礎上，遵循“由表及里、由淺入深、循序漸進”的工作思路，首先，進行高密度電阻率法淺部探測；其次，在其基礎上進行可控源音頻大地電磁法中深部探測；然后，為了對比試驗增加了廣域電磁法進行大深度探測；最后，利用天然場面波微動新方法在擬布設鉆探位置從地震波速物性層面進行了多參數探測分析。

2.1.1 高密度電阻率法

眾所周知，高密度電阻率法值陣列布設電極的高效率電法探測方法，是傳統電法測深的改良，具有較高的野外數據采集效率，目前，在城市地質調查工作中應用最為廣泛。本次研究中，為了探測近地表淺層100 m以內的地層結構，點距選取10 m，一個排列共布設60個電極，共計590 m，實行連續滾動方式，采集的參數為視電阻率，因視電阻率是地層信息的整體反應，采用Res2Dinv反演軟件對其反演，其中，正演利用有限差分法進行計算，反演利用基于圓滑約束的最小二乘法進行計算，利用重慶奔騰數控技術研究所研制的WGMD-4電阻率數據采集系統。本次研究是以探測地層分層為主要目的，因此，選擇縱向分辨率較高的溫納裝置[15]。

2.1.2 可控源音頻大地電磁法

因探測深度未穿透低阻覆蓋層，在高密度電法完成的基礎上，在探測剖面基礎上完成了可控源音頻大地電磁法。可控源音頻大地電磁法（以下簡稱“CSAMT法”）是在大地電磁探測基礎上發展而來的一種人工源頻率域測深方法，其信號強度較大，探測深度最深可達1.5 km，但其目前處理野外實測數據時一般采用“遠區”數據，因此，收發距的布設是方法有效使用的前置條件之一。剖面與高密度電阻率法一致，點距20 m，利用曹創華等[16]以長株潭城市群地質調查項目成果的新認識，把供電源AB布設在測線剖面的西南部陬市鎮楊家巷村附近，供電線AB長1.04 km，AB線與測線剖面平行，收發距為9.8 km，觀測裝置為TM模式，采用儀器是美國Zonge公司研發的GDP32系列多功能電法工作站，發電機發射功率為30 kW，發射電流呈現出高頻低電流，低頻高電流的特點，電流最大為16 A，最小為3.5 A，采集數據頻率為0.125～8192 Hz。數據處理軟件Scs2D，具體處理過程有飛點剔除、靜態校正、地形改正和反演等步驟，采用以一維Bostick反演結果為初始模型的二維Occam反演方法進行反演。

2.1.3 廣域電磁法

廣域電磁法由中南大學何繼善院士發明，是21世紀初發展起來的電磁探測技術。該方法采用人工場源，從電磁場精確的表達式出發，嚴格定義了廣域電磁法視電阻率參數，改善了非遠區的畸變效應，使得測深能在廣大的、不局限于“遠區”的區域進行，在同等收發距條件下勘探深度更大。該方法繼承了可控源音頻大地電磁法使用人工場源的優點，以及磁偶源頻率測深法（以下簡稱“MELOS”法）非遠區測量的優勢；改良了CSAMT法遠區信號微弱的劣勢，拓展了觀測適用的范圍，同時摒棄了MELOS方法的校正辦法；用適合于全域的公式計算視電阻率，保留了計算公式的高次項，大大拓展了人工源電磁法的觀測范圍，提高了觀測速度、精度和野外效率[17]。本次研究利用可控源音頻大地電磁法采集數據，而不使用磁道數據，僅利用電場標量和地電場觀測參數進行全區電阻率反演，利用頻率為0.125～8192 Hz全頻。

2.1.4 天然源面波微動

因地球近地表人文活動區域都有一種微弱的連續天然面波信號場源，自從20世紀50年代就有學者利用其探測計算地層的速度結構。目前，主要利用低頻信號如由潮汐、海浪、氣壓變動產生的信號和人文活動產生的機械振動高頻信號，根據這些瑞雷面波特性處理反演得到的頻散曲線進行推斷地下橫波速度模型[10]。本次研究利用北京水電物探研究所資助研發的WD-1型智能天然源面波數據采集系統，數據采集利用2 Hz檢波器10通道等邊三角形內嵌布極方式進行數據采集，最大邊長為90 m、采樣間距為5 ms、迭代次數為600次。

2.2 研究區地球物理特征

2.2.1 常德市區域地球物理特征

區域1∶50萬重力測量資料顯示：研究區位于麻陽-常德-湘潭-衡陽弧形重力高的東部和南西邊緣，布格重力異常值為（-50～-10）×10-5m/s2。該弧形重力高與相應的中生代紅層盆地輪廓基本吻合，相對基底有較大的質量虧損，主要反映了變質結晶基底隆起及莫霍面的抬升。

研究區內航磁背景場為△T正異常，△T異常值在15～40 nT之間變化；航磁△T異常總體上在北西部和北東部小，從北部往南部，△T異常梯度變化由大變小；航磁異常等值線由中部向北西方向和北東方向凸出。研究區內航磁異常等值線方向與地層和構造走向線不相協調，反映測區出露的第四紀地層具有低緩磁場或無磁性的特征。

研究區位于常德幔隆區，該幔隆區包括整個洞庭湖地區，常德-安仁轉換斷裂帶將其分為沅江、桃源兩個走向完全不同的局部幔隆區。研究區位于兩個次級幔隆區之間，常德-安仁轉換斷裂帶從研究區中部呈北西向通過，莫霍面埋深31～32 km，兩個幔隆區對應于中生代洞庭盆地。

2.2.2 研究區巖石物理特征

本次試驗的四種探測方法涉及到兩種參數，即電性參數電阻率和波速。

研究區內的不同地層巖石電阻率統計見表1，其電性特征為：第四系全新統、中更新統馬王堆組、下更新統常德組呈相對低阻特征，其均值一般都小于300 Ω·m，若含礫石較多則電阻率相對高一些，而粉砂質黏土成分多一些則呈現較低阻特征；此地區奧陶系白水溪組粉砂質頁巖電阻率均值為1000.3 Ω·m，在此區域埋藏較深，含水量較少呈現出中等電阻率特征；寒武系探溪組白云質灰巖是本區探測深度范圍內電阻率最高的地層，均值為2788.5 Ω·m；寒武系污泥塘組帶泥晶泥質灰巖相對電阻率變低，均值為1757.3 Ω·m，因其含有泥質成分相對探溪組有所降低；而寒武系牛蹄塘組因含炭質，相比區內寒武系其他地層電阻率急劇降低，均值為798.2 Ω·m；當隱伏的地層內有斷裂構造時，因上部覆蓋的第四系具有天然的湖湘盆地地表水滲透通道而呈現局部低阻特征。不同地層的物性具有一定差異，可以通過電阻率參數探測進行研究區內隱伏前第四系基巖的探測，用于圈定構造等地質異常體。

表1 研究區不同地層巖石電阻率特征表Table 1 Resistivity characteristics of different strata in the study area

研究區內的不同地層橫波波速范圍見表2，場地內出露地層主要為第四系，對于百米以內的地表黏土、礫巖和基底的前第四系基巖具有明顯的波速差異，特別是第四系內的含礫石層可以有效區分，說明此區內具有開展天然源面波的前提條件。

表2 研究區不同地層巖石橫波速度范圍Table 2 S-wave velocity range of different strata in the study area

3 成果解釋及鉆探驗證

3.1 成果解釋

對設計完成的G4剖面高密度電法的原始數據按照Res2Dinv軟件的輸入格式要求，拼接整理成一個獨立的數據文件（圖2）。利用基于圓滑約束的最小二乘法進行反演，經過地形校正得到圖2（a）中的結果。整個剖面在海拔-50 m以上電阻率值都較低，一般小于2000 Ω·m，結合地表地質填圖調查結果，推斷大部分為第四系，且低阻特征比較連續，綜合表1結果，推斷在測點1920前和測點2450后的近地表均為中更新統馬王堆組網狀黏土，其下伏地層為第四系全新統礫石層。

圖2 研究區G4試驗剖面探測結果集地質解釋Fig. 2 Geological interpretation of G4 test profile in the study area

圖2（b）為可控源音頻大地電磁法數據處理結果，整體上地層分層較為明顯，可分為3層，第一層為海拔-100 m以上的低阻層，電阻率一般小于2000 Ω·m；第二層為海拔-550 m以上的中阻層，但呈現類似斷陷盆地的兩個凹陷，第一個位于測點1650～測點2120之間，另一個位于測點2320～測點2920之間，第一個凹陷電阻率無明顯的帶狀分布特征且電阻率梯級帶變化較緩，推斷其由地層引起，第二個凹陷有一個朝向東東南的低阻帶延伸且梯級帶變化明顯，推斷是深部基巖內存在構造，因淺部發育的水系滲漏而形成低阻，深部-650～-950 m海拔之間呈現出高阻特征，推斷為寒武系灰巖、白云質灰巖等地層。但值得注意的是可控源音頻大地電磁法反演剖面存在若干個“公牛眼”形態的閉合小異常，且分布毫無規律，例如，測點2520～測點2620海拔-200 m處、測點2820～測點2920海拔-150 m處等位置，側面證實了可控源音頻大地電磁法在城區嚴重受磁場干擾而引起電磁場信號跳變，給解釋帶來了不便。

為了加強對試驗剖面的地質認識，利用GDP32采集的電場數據，由中南大學劉春明老師處理完成得到廣域電磁法數據處理結果，如圖2（c）所示。由圖2（c）可知，成圖數據較為圓滑穩定，分層較為明顯；結合上述兩種方法對剖面進行了綜合解釋，如圖2（d）所示。由圖2（d）可知，淺部-100 m以上為第四系，僅在大號測點3220～測點3720之間存在一個高阻，經實地調查可知，因在工業開發區路邊的人文活動密集區采集，可能淺部受到了工業游散電流的干擾；覆蓋于第四系下的白堊系地層呈現出相對中等電阻率特征，電阻率在600 Ω·m左右，中間部位相對缺失，推斷在測點2120前和測點2720后存在白堊系地層；中間測點2120～測點2520之間被剖面的高阻體拱起，推斷為寒武系地層，其頂部與第四系呈不整合接觸；在深部的寒武系地層中結合區域地質調查研究結果，推斷在測點2520的深部存在一個寬約80 m向東南傾的斷層破碎帶，其中，寒武系探溪組白云質灰巖和污泥塘組帶泥晶泥質灰巖被構造運動改造；在斷裂深部的東南部存在寒武系牛蹄塘組含炭質地層，產狀傾向西北，傾角約為30°。

電阻率類探測方法對識別大的地質構造、地層分布具有良好的分辨率，但對于第四系內進行更加細化的分層則效果不佳，這一點在以湖湘盆地的典型水系發育的南方洞庭盆地常德地區顯得更為突出。為了更好地為基巖探測提供多參數信息，本次試驗在擬選的sk09驗證鉆孔（測點2120）附近、sk10驗證鉆孔（測點3700）附近分別完成了天然面波微動探測，驗證鉆孔探測結果分別如圖3和圖4所示。

圖3 sk09天然源面波微動探測頻散波速反演巖性分層結果Fig. 3 Lithographic layering results of sk09 natural source surface wave fretting detection of dispersive wave velocity

圖4 sk10天然源面波微動探測頻散波速反演巖性分層結果Fig. 4 Lithographic layering results of sk10 natural source surface wave fretting detection of dispersive wave velocity

在反演后的面波頻散曲線中，利用曲線拐點、曲線的切線突變點、波速變化等信息進行分層，結合表2波速差異信息，綜合地質認識sk09可分為5層，sk10可分為8層。其中，sk09推斷第一層為在鉆孔地下18 m以上，波速在159～446 m/s之間，為粉質黏土和粉砂；第二層在地表下40～18 m之間，波速在547～581 m/s之間，為土夾卵礫石層；第三層在地表下75～40 m之間，波速在575～740 m/s之間，為一卵石層；第四層在地表下145～75 m之間，波速在825～930 m/s之間，為頁巖層；第五層在地表下145 m以下，波速在992～1400 m/s之間，為灰巖層。sk10推斷第一層在地表下10 m以上，波速為306 m/s，為粉質黏土層；第二層在地表下70～10 m，波速為504 m/s，為土夾卵礫石層；第三層在地表下110～70 m，波速為1093 m/s，為粉砂質泥巖層；第四層在地表下180～110 m，波速為1400 m/s，為粉砂巖；第五層在地表下210～180 m，波速為1303 m/s，為泥質粉砂巖；第六層在地表下270～210 m，波速為1473 m/s，為粉砂巖；第七層在地表下310～270 m，波速為1510 m/s，為泥巖；第八層在地表下320～310 m，波速為1653 m/s，為灰巖。

3.2 鉆探驗證

依據綜合探測解釋結果，對圖1和圖2中的sk09和sk10分別進行鉆探驗證，其中，sk09終孔深度為202.18 m（圖5），sk10終孔深度為302.80 m（圖6）。sk09地表下89.50 m以上為第四系，其中，3.70 m以上為耕植土，3.70～18.30 m為粉質黏土，18.30～89.50 m為砂卵石層；89.50～147.80 m為奧陶系白水溪組粉砂質頁巖；147.80～202.18 m為寒武系探溪組白云質灰巖。sk10地表下9.70 m為中更新統馬王堆組網狀黏土，地表下9.70～71.52 m為第四系全新統礫石層；地表下71.52～275.74 m為白堊系紅層；275.74 m以下為寒武系污泥塘組晶泥質灰巖。上述鉆探結果與微動探測結果對應較好，特別是呈現低阻特征的第四系和白堊系紅層的內部分層較為準確。

圖5 sk9鉆探巖芯編錄結果Fig. 5 Cataloging results of sk9 drilling core

圖6 sk10鉆探巖芯編錄結果Fig. 6 Cataloging results of sk10 drilling core

4 結 論

1）高密度電阻率法在湖湘盆地地區利用常規10 m電極距的電纜布極，在常德市區地層結構條件下難以穿透地阻層；可控源音頻大地電磁法探測深度較大，但其視電阻率計算時必須考慮磁場信息，在強干擾環境下會對反演結果造成誤差放大，會對深部地層結構信息的判斷帶來一定的干擾；廣域電磁法具有信號強，利用單個分量進行處理，從根本上降低了信號的干擾，最大限度地提高了信噪比，且可利用處理反演的頻帶較寬，探測深度較大；天然面波微動有效地提高了低阻覆蓋區近地表探測的分辨率，可對第四系和白堊系紅層進行分層反演，探測效果明顯。

2）在城市地質調查工作中需結合多種方法優勢，利用多參數方法進行綜合評價，遵循由宏觀到微觀，由表及里，仔細分析判別各種異常產生的原因，綜合地質調查研究結果進行物探資料解釋，是提高城市地球物理勘探效果應遵循的基本思路。

致謝：本次試驗研究工作歷時近兩年，湖南省地質調查院物化遙地質分院和水工環分院（現湖南省地質調查所物化遙地質技術中心和水工環地質研究中心）眾多同事參與了野外工作；天然源面波技術工作得到了北京水電物探研究所范宏強工程師的耐心指導和幫助，在此一致向他們表示衷心的感謝！范宏強工程師的耐心指導和幫助，在此一致向他們表示衷心的感謝！