張靜 胡尊平 陳實(shí) 陶鵬飛 牛輝

摘? ?要：利用Surfer軟件實(shí)現(xiàn)物探數(shù)據(jù)的三維可視化。通過(guò)數(shù)據(jù)網(wǎng)格化利用成果數(shù)據(jù)制作等值線圖、據(jù)高程數(shù)據(jù)制作表面圖、采用點(diǎn)位坐標(biāo)制作張貼圖，最后通過(guò)覆蓋地圖功能將物探異常數(shù)據(jù)投影在三維地形圖上，實(shí)現(xiàn)三維可視化。同時(shí)介紹了一些成圖細(xì)節(jié)，為物探異常推斷解釋提供參考。

關(guān)鍵詞：Surfer;三維可視化;物探異常;物探數(shù)據(jù);立體地形投影

物化探常用軟件成圖效果單一，往往是用鮮明顏色區(qū)分和突出異常成果，多為二維圖形，無(wú)法清晰直觀地表達(dá)地形和異常數(shù)據(jù)關(guān)系。常用地質(zhì)軟件流程復(fù)雜，不具備將異常成果圖和實(shí)際點(diǎn)位圖投影在三維地形圖上的功能。Surfer軟件較好的解決了此問(wèn)題，其三維成圖功能強(qiáng)大，操作簡(jiǎn)便，已廣泛應(yīng)用于教學(xué)、科研、測(cè)繪、地質(zhì)等領(lǐng)域[1]。由于測(cè)區(qū)地形起伏使物探成果數(shù)據(jù)在激發(fā)和采集測(cè)量過(guò)程中發(fā)生了畸變，原來(lái)適用于水平地形條件下的異常解釋方法，如不考慮起伏地形影響，直接用來(lái)解釋物探異常，將導(dǎo)致較大誤差。本文將物探異常直接套合在地形圖上，為異常推斷解釋提供直觀參考。

1? Surfer軟件介紹

Surfer軟件是美國(guó)Golden Software公司出品的一款著名地理數(shù)據(jù)網(wǎng)格化繪圖軟件，擁有強(qiáng)大插值功能和繪制圖件能力，是地質(zhì)工作者必備的繪圖軟件之一（圖1）。該軟件簡(jiǎn)單易學(xué)，可輕松制作基面圖、數(shù)據(jù)點(diǎn)位圖、分類數(shù)據(jù)圖、等值線圖、線框圖、地形地貌圖、趨勢(shì)圖、矢量圖及三維表面圖等。Surfer被廣泛用于地形建模、測(cè)深建模、景觀可視化、表面分析、等高線測(cè)繪、流域和3D曲面的繪制等[2-8]。

2? 數(shù)據(jù)準(zhǔn)備

在開(kāi)始制作三維可視化立體地形投影圖之前，需將投影成果數(shù)據(jù)準(zhǔn)備好。本次繪圖之前準(zhǔn)備了調(diào)查區(qū)點(diǎn)位數(shù)據(jù)、磁法測(cè)量數(shù)據(jù)（ΔT）、視極化率數(shù)據(jù)（M1）、視電阻率數(shù)據(jù)和高程（地形）數(shù)據(jù)，即在excel表格中將橫坐標(biāo)、縱坐標(biāo)、成果數(shù)據(jù)分別成列，保存為成果數(shù)據(jù)“.csv”格式。之所以不保存成“.xls”格式，是因?yàn)镾urfer 8軟件無(wú)法較好的識(shí)別其格式。而“.csv”格式文件會(huì)將其中公式及格式全部清除，方便Surfer軟件識(shí)別。

3? 繪圖流程

3.1? 將成果數(shù)據(jù)進(jìn)行網(wǎng)格化處理

打開(kāi)Surfer軟件，點(diǎn)擊菜單欄中的網(wǎng)格→數(shù)據(jù)，選擇準(zhǔn)備好的“成果數(shù)據(jù).csv”，將XYZ軸分別對(duì)應(yīng)設(shè)置。在高級(jí)選項(xiàng)中可根據(jù)實(shí)際測(cè)量點(diǎn)線距設(shè)置搜索橢圓半徑和角度，其他參數(shù)可根據(jù)需要自行設(shè)置。本次實(shí)測(cè)數(shù)據(jù)采用線距100 m，點(diǎn)距20 m，方位角20°，將搜索半徑設(shè)置為點(diǎn)線距的2～3倍既能反應(yīng)出細(xì)小異常，又能使網(wǎng)格化數(shù)據(jù)較為平滑。選擇“確認(rèn)”則很快創(chuàng)建了“磁測(cè)數(shù)據(jù)（△T）”網(wǎng)格化數(shù)據(jù)。以此類推，再分別創(chuàng)建“視極化率M1（%）”、“視電阻率（Ω·m）”和“高程（m）”的網(wǎng)格化數(shù)據(jù)。

3.2? 利用創(chuàng)建的成果數(shù)據(jù)網(wǎng)格? ? ? ? ?化文件，制作等值線圖

在菜單欄中選擇地圖→等值線圖→新建等值線圖，選擇打開(kāi)“磁測(cè)數(shù)據(jù)（△T）.grd”，則在圖形顯示區(qū)生成了磁測(cè)等值線圖[10]。右鍵點(diǎn)擊等值線圖，則可以進(jìn)行標(biāo)注和屬性的設(shè)置。在屬性設(shè)置中，選擇“常規(guī)”菜單欄，將填充等值線和平滑等值線勾選，在“等級(jí)”菜單欄選擇填充→前景色，進(jìn)行顏色譜的設(shè)置。點(diǎn)擊確定，則完成了等值線圖的“上色”（圖2）。制作其他等值線圖以此類推，分別制作“視極化率M1（%）”、“視電阻率（Ω·m）”等值線圖。

3.3? 利用創(chuàng)建的高程網(wǎng)格化? ? ? ? ?文件（高程.grd），制作? ? ? ? ? ? ? “表面圖”

在菜單欄中選擇地圖→表面圖，選擇打開(kāi)“高程（m）.grd”文件，則在圖形顯示區(qū)生成三維地形圖（圖3）。

在對(duì)象管理器中，右鍵點(diǎn)擊“3D Surface”，選擇屬性，可查看或調(diào)整圖形顯示顏色、網(wǎng)格、光照角度等。在繪圖工具欄中，選擇“3D軌跡球”或菜單欄“地圖→軌跡球”可對(duì)圖形進(jìn)行旋轉(zhuǎn)、翻動(dòng)等操作，使圖形更直觀。

3.4? 利用“成果數(shù)據(jù).csv”，制作“張貼圖”

在菜單欄中選擇地圖→張貼圖→新建張貼圖，選擇“成果數(shù)據(jù).csv”，則在圖形顯示區(qū)生成張貼圖（圖4）。在對(duì)象管理器中，右鍵點(diǎn)擊“Post”，選擇屬性，可查看“X坐標(biāo)”和“Y坐標(biāo)”，同時(shí)可重新設(shè)置標(biāo)注的符號(hào)類型、顯示頻率及符號(hào)大小等。

3.5? 覆蓋地圖，完成三維可視化立體地形投影圖的? ? ? ? ? 制作

參考以上步驟，我們已經(jīng)完成“等值線圖”、“表面圖”和“張貼圖”的制作。接下來(lái)我們需將3幅圖進(jìn)行覆蓋處理，讓各元素呈現(xiàn)在一張圖上。在對(duì)象管理器中，同時(shí)選中以上3幅圖（利用“ctrl+”功能），然后在菜單欄中選擇地圖→覆蓋地圖，則在圖形顯示區(qū)生成“三維可視化立體地形投影圖”（圖5）。在對(duì)象管理器中，我們?nèi)钥梢杂益I點(diǎn)擊修改各類圖件的屬性，這為我們后期調(diào)整一些細(xì)節(jié)提供了方便。最后，如對(duì)某一幅圖形效果不滿意，可在對(duì)象管理器中選定該圖形，然后在菜單欄中選擇地圖→拆解覆蓋，即可將該圖形取消覆蓋。

3.6? 細(xì)節(jié)調(diào)整

顯示高程差? 在對(duì)象管理器中右鍵點(diǎn)擊3D Surface→屬性→勾選顯示基面，可用所需顏色顯示高程差。

增加圖例單位? 在圖形顯示區(qū)右鍵點(diǎn)擊圖例，選擇屬性→格式，在“后綴”欄輸入對(duì)應(yīng)的單位。

增加X(jué)YZ軸坐標(biāo)? 在對(duì)象管理器相應(yīng)的軸屬性中可顯示、隱藏、調(diào)整、加密刻度值。

成果圖件的輸出? 點(diǎn)擊菜單欄→文件→輸出，可將成果圖件保存成20多種格式的文件。

白化文件的建立與應(yīng)用? 成圖后，圖形有時(shí)會(huì)出現(xiàn)邊界有毛刺現(xiàn)象，這是因?yàn)樵陂_(kāi)始數(shù)據(jù)網(wǎng)格化時(shí)，選擇的網(wǎng)格搜索半徑過(guò)大或過(guò)小，這時(shí)就需進(jìn)行白化處理，即將某范圍內(nèi)或以外線條進(jìn)行剪裁。方法如下：①利用Excel（如用WPS軟件，可能出現(xiàn)無(wú)法復(fù)制數(shù)據(jù)現(xiàn)象）建立白化節(jié)點(diǎn)表格，節(jié)點(diǎn)要有順序（表1）;②在菜單欄點(diǎn)擊文件→新建→工作表，將Excel中的數(shù)據(jù)復(fù)制到上表中，文件→另存為，在保存類型中選擇“*.bln”格式;③在菜單欄點(diǎn)擊網(wǎng)格→白化，選擇需要白化的“*.grd”文件，再打開(kāi)相應(yīng)的“*.bln”文件，修改文件名，保存;④在菜單欄點(diǎn)擊地圖→等值線圖→新建等值線圖，選擇并打開(kāi)上述完成的文件，即可得到白化后的圖形。

4? 成圖效果分析

利用以上流程，可將各項(xiàng)異常成果進(jìn)行三維可視化處理，并按需排列放置在一個(gè)Surfer文件中，利于成果對(duì)比分析研究。圖6-a反應(yīng)了調(diào)查區(qū)地形變化情況，同時(shí)也反應(yīng)了物探測(cè)量點(diǎn)位在不同海拔高度上的位置情況;圖6-b為磁測(cè)異常在調(diào)查區(qū)空間位置上的反應(yīng)，由此可較清晰判斷磁異常形態(tài)分布特征是否與地形變化有關(guān);圖6-c，e分別為視極化率和視電阻率特征，它們主要分布在地形較低的南部，異常走向?yàn)镋W向。同時(shí)，調(diào)查區(qū)東部云海調(diào)查區(qū)均能看到弱異常反應(yīng)。圖6-d，f分別為布格重力和剩余重力異常特征。圖6-e顯示布格重力異常值由南向北，由低到高逐漸增大;圖6-f顯示剩余重力異常在云海礦區(qū)及海拔不同地區(qū)均有所反應(yīng)。

綜上，三維立體地形投影圖清晰、直觀的反應(yīng)了各類物探異常的空間位置及分布形態(tài)，為物探解釋和地質(zhì)推斷工作提供了助力。

5? 結(jié)語(yǔ)

（1） 利用Surfer軟件分別生成等值線圖、張貼圖、表面圖。通過(guò)調(diào)整圖形屬性及覆蓋地圖功能，很好地實(shí)現(xiàn)了將物探異常成果投影在三維地形圖上的功能。且可在軟件中任意角度旋轉(zhuǎn)圖件，成圖過(guò)程簡(jiǎn)單易學(xué)，成圖效果美觀易懂，是物化探工作者展示成果的又一得力助手。存在的不足之處，主要有：生成的成果文件無(wú)法直接使用“MapGIS”軟件打開(kāi)，不能進(jìn)一步加工，需進(jìn)行相應(yīng)的格式轉(zhuǎn)換。

（2） 同一巖性相同地質(zhì)條件下，地形凹陷與地形隆起處對(duì)應(yīng)的物探異常，認(rèn)為是由地形參與引起的變化，這為排除干擾異常，消減多解性，研究異常原因，做出推斷解釋具重要作用。

（3） 三維地形圖的制作對(duì)高程數(shù)據(jù)要求較高，普通的GPS精度可能無(wú)法滿足數(shù)據(jù)質(zhì)量的要求，建議采用RTK（高精度測(cè)量?jī)x）進(jìn)行高程數(shù)據(jù)的采集。

參考文獻(xiàn)

[1]? ? 馬肖波.綜合物探數(shù)據(jù)的聯(lián)合成圖技術(shù)及應(yīng)用[D].成都理工大學(xué).

[2]? ? 朱慶俊， 李風(fēng)哲. Surfer軟件在地球物理資料處理中的應(yīng)用[J]. 物? ?探與化探，2007（3）：250-251.

[3]? ? 陳少?gòu)?qiáng)， 宋利好， 姚敬金，等. 可視化技術(shù)在物化探找礦中的應(yīng)用及 前景[J]. 物探與化探， 2002.

[4]? ? 任印國(guó)， 魏永強(qiáng). 使用Surfer軟件繪制地質(zhì)圖件和處理地質(zhì)數(shù)據(jù)的方法[J]. 測(cè)繪技術(shù)裝備， 2006， 8（1）：34-36.

[5]? ? 劉宏偉， 徐明， 萇沙沙，等. ArcGIS、MapGIS與Surfer軟件聯(lián)合構(gòu)建三維地形及其可視化[J]. 地理空間信息， 2012， 10（3）：67-69.

[6]? ? 蔣偉. Surfer在繪制物探圖件中的新方法技術(shù)[J]. 礦產(chǎn)勘查， 2019.

[7]? ? 密蓓蓓， 徐濤， 陸凱，等. 基于MapGIS和Surfer軟件的三維地質(zhì)圖制作[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì)， 2014， 34（6）：181-185.

[8]? ? 陳歡歡， 李星， 丁文秀. Suffer 8.0等值線繪制中的十二種插值方法[J]. 工程地球物理學(xué)報(bào)， 2007， 4（1）：52-57.

[9]? ? 劉兆平， 楊進(jìn)， 武煒. 地球物理數(shù)據(jù)網(wǎng)格化方法的選取[J].物探與化探， 2010（1）：97-101.

[10]? 宋明藝， 張春灌. 借助Suffer軟件實(shí)現(xiàn)快速繪制平面等值線圖[J].工程地球物理學(xué)報(bào)， 2009.

Application of Surfer Software in Making 3D Visualization Stereo Topographic Projection Map

——Take Geophysical Data of Yunhai as an Example

Zhang Jing，Hu Zunping，Chen Shi，Tao Pengfei，Niu Hui

（Geological Survey Institute of Xinjiang，Urumqi，Xinjiang，830000，China）

Abstract：Surfer software is used to realize the three-dimensional visualization of geophysical data. Through data gridding， contour map is made by using results data， surface map is made by using elevation data， and post map is made by using point coordinates. Finally， the geophysical anomaly data is mapped on the three-dimensional topographic map by covering map function to realize three-dimensional visualization. At the same time， some drawing details are introduced to make the results figures more beautiful and high-definition. It provides a reference for geophysical anomaly inference and interpretation.

Key words：Surfer;Three-dimensional visualization;Geophysical anomaly;Geophysical data;Stereo terrain projection