徐 陽

(山東省煤田地質(zhì)局第一勘探隊(duì)，山東青島 266555)

0 引言

在多金屬礦勘探過程中，往往采用磁法勘探和激電中梯的方法進(jìn)行工作區(qū)域的勘探(崔東鄆和王偉建，2009；管志寧，2015；王猛和梁中朋，2015)，或運(yùn)用綜合物探方法分區(qū)域探測(cè)多金屬礦在測(cè)區(qū)的分布狀況(李榮亮等，2017)。多金屬礦往往分布在地理環(huán)境較為復(fù)雜的地域，因此尋找創(chuàng)新型的找礦方法成為了重要的研究對(duì)象(任為望等，2012；郭世強(qiáng)，2013；王曉軍等，2013；林斌，2018；邰文星等，2019；李俊杰等，2019)。而地質(zhì)勘探往往處在人煙稀少的地段，環(huán)境比較惡劣，嚴(yán)重威脅到野外工作人員的生命安全，因此需要尋找合適的方法以大大降低野外工作環(huán)境的危險(xiǎn)程度(苑立清和趙金凱，2009)。因而，施工前的工作部署圖的合理規(guī)劃以及坐標(biāo)點(diǎn)的布置成為了重中之重(曾群望，1985；劉勵(lì)慎和汪津華，1987；臧忠淑和蔣樹一，2003；艾斯卡爾和李華，2010；王寶江，2010；張建恩等，2010)，而將測(cè)繪遙感等技術(shù)合理應(yīng)用到地質(zhì)勘探工作部署圖的規(guī)劃中可以大大提高工作效率(包世泰等，2004；康世英和張宏偉，2006；胡志文等，2012；張維寬等，2012；何輝等，2014；徐力峰，2014；狄廣禮和劉劍英，2015；周朝憲等，2015；李飛等，2018；王浩峰等，2019)。在眾多測(cè)繪軟件中，CASS具有非常強(qiáng)大的地質(zhì)圖編輯功能，及坐標(biāo)點(diǎn)導(dǎo)入導(dǎo)出的功能，因此運(yùn)用CASS軟件對(duì)多金屬礦的工作部署圖進(jìn)行轉(zhuǎn)換具有十分強(qiáng)大的優(yōu)點(diǎn)(李得基等，2006；韋晶珠，2010；李丙法，2017)。

結(jié)合以往的研究發(fā)現(xiàn)，本文首次提出運(yùn)用測(cè)繪程序CASS 8.0來直接繪制地質(zhì)勘查類的工作部署圖，并利用測(cè)繪學(xué)科軟件與地質(zhì)學(xué)科知識(shí)的交叉應(yīng)用優(yōu)勢(shì)，使數(shù)據(jù)導(dǎo)入導(dǎo)出操作簡(jiǎn)便化，易于掌握，并且繪圖精確度提高，使點(diǎn)坐標(biāo)直接呈現(xiàn)在圖形上，便于后期修改調(diào)整。

針對(duì)大面積比例尺坐標(biāo)點(diǎn)導(dǎo)出導(dǎo)入的困難問題，本文嘗試直接運(yùn)用測(cè)繪軟件CASS 8.0中對(duì)高程點(diǎn)的繪制和輸出，直接轉(zhuǎn)換為帶有測(cè)點(diǎn)線號(hào)的點(diǎn)坐標(biāo)。該點(diǎn)坐標(biāo)無須經(jīng)過再次變換，經(jīng)Excel調(diào)整變換即可按線號(hào)批量導(dǎo)入手持GPS定位儀中。

CASS 8.0在地貌、地物繪制，點(diǎn)坐標(biāo)的繪制、提取等方面具有明顯的優(yōu)勢(shì)。為了有效解決坐標(biāo)點(diǎn)批量生成及導(dǎo)入問題，本文首次提出運(yùn)用測(cè)繪軟件CASS 8.0的點(diǎn)陣排列工具及有編碼高程點(diǎn)的批量輸出的功能，輸出帶有線號(hào)以及序號(hào)的點(diǎn)坐標(biāo)。此方法操作簡(jiǎn)單快捷，無需通過編程或應(yīng)用Excel表格方程式，即可得出直觀并且方便測(cè)繪儀器使用的數(shù)據(jù)及工作部署圖，可極大提高坐標(biāo)點(diǎn)提取的可操作性及簡(jiǎn)便性。

1 技術(shù)方法

1.1 測(cè)線坐標(biāo)點(diǎn)的生成

在勘探野外工作中，采用手持GPS定位儀所需要的測(cè)點(diǎn)坐標(biāo)包括四個(gè)部分：橫坐標(biāo)、縱坐標(biāo)、線號(hào)、點(diǎn)號(hào)。文件格式為(李榮亮等，2017)：

橫坐標(biāo) 縱坐標(biāo) 線號(hào) 點(diǎn)號(hào)

4654445.415 48860.885 100 001

4654485.957 48860.609 100 002

4654524.640 48861.241 100 003

……

首先使用CASS 8.0繪制出需要的測(cè)點(diǎn)布陣圖，打開CASS軟件，點(diǎn)擊“繪圖處理”中的“展野外測(cè)點(diǎn)點(diǎn)號(hào)”，利用GPS定位儀測(cè)出的工區(qū)邊界點(diǎn)坐標(biāo)繪制出工區(qū)邊界。運(yùn)用CASS中的對(duì)象捕捉功能，將邊界點(diǎn)連成線，繪制首條測(cè)線，并運(yùn)用陣列復(fù)制功能繪制工區(qū)的所有測(cè)線，在每條測(cè)線端點(diǎn)處輸入測(cè)線號(hào)，輸入命令g，捕捉每條測(cè)線與邊界點(diǎn)的交點(diǎn)，并輸入相應(yīng)的測(cè)線號(hào)來作為高程點(diǎn)高程，選中所繪制高程點(diǎn)，點(diǎn)擊“陣列復(fù)制”輸入角度、間距、數(shù)目等信息，刪除工區(qū)邊界外高程點(diǎn)。繪制一條閉合的邊界線將所有高程點(diǎn)包含在內(nèi)，點(diǎn)擊“工程應(yīng)用”——“高程點(diǎn)生成數(shù)據(jù)文件”——“有編碼高程點(diǎn)”，此時(shí)操作平臺(tái)上會(huì)自動(dòng)彈出“輸入數(shù)據(jù)文件名”的對(duì)話框，選擇相應(yīng)文件并點(diǎn)擊“打開”按鈕，輸入數(shù)據(jù)文件名后出現(xiàn)請(qǐng)選擇(1)選取高程點(diǎn)范圍(2)直接選取高程點(diǎn)或控制點(diǎn)，輸入1，點(diǎn)擊回車按鈕，此時(shí)命令欄出現(xiàn)“請(qǐng)選取建模區(qū)域邊界”的命令，點(diǎn)擊繪制的閉合邊界，即可生成所需坐標(biāo)文件(*.dat格式)(林斌，2018)。文件格式為：

測(cè)點(diǎn)序號(hào) 縱坐標(biāo) 橫坐標(biāo) 高程點(diǎn)

1， ，48860.885 4654445.415 1

2， ，48860.609 4654485.957 1

3， ，48861.241 4654524.640 1

……

對(duì)于CASS 8.0導(dǎo)出的坐標(biāo)文件，利用Microsoft Excel對(duì)數(shù)據(jù)進(jìn)行處理，即可得到野外勘探工作的坐標(biāo)文件。具體步驟如下：①打開生成的測(cè)點(diǎn)坐標(biāo)；②對(duì)數(shù)據(jù)進(jìn)行分裂；③按高程點(diǎn)的數(shù)值對(duì)整個(gè)表格進(jìn)行重新排序，并改為測(cè)線序號(hào)；④將測(cè)點(diǎn)序號(hào)調(diào)整到最后一列，并重新按照測(cè)線的線號(hào)更改對(duì)應(yīng)的序號(hào)；將縱坐標(biāo)列與橫坐標(biāo)列互換。具體操作流程圖如圖1所示：

圖1 CASS 8.0生成坐標(biāo)流程Fig.1 Flowchart of CASS 8.0 coordinate generation

1.2 工作部署圖的繪制

將生成的坐標(biāo)文件按照以下格式重新修改并保存，文件格式為：

測(cè)點(diǎn)序號(hào) 縱坐標(biāo) 橫坐標(biāo) 高程點(diǎn)

1， ，48860.885 4654445.415 1

2， ，48860.609 4654485.957 1

3， ，48861.241 4654524.640 1

……

具體步驟如下：①利用Microsoft Excel打開生成的測(cè)點(diǎn)坐標(biāo)；②將測(cè)點(diǎn)序號(hào)調(diào)到第一列，并在后面插入新的一列；③將橫坐標(biāo)與縱坐標(biāo)對(duì)調(diào)；④將文件保存為*.csv格式；⑤將*.csv格式的文件修改為*.dat格式文件。

在CASS 8.0中打開相應(yīng)的地形地質(zhì)圖，點(diǎn)擊“繪圖處理”——“展野外測(cè)點(diǎn)點(diǎn)號(hào)”，此時(shí)操作平臺(tái)上會(huì)自動(dòng)彈出“輸入坐標(biāo)數(shù)據(jù)文件名”的對(duì)話框，選擇以上生成的*.dat格式文件，并點(diǎn)擊“打開”，即可將野外測(cè)點(diǎn)點(diǎn)號(hào)展到圖上。繼續(xù)點(diǎn)擊“繪圖處理”——“展高程點(diǎn)”，此時(shí)操作平臺(tái)上會(huì)自動(dòng)彈出“輸入坐標(biāo)數(shù)據(jù)文件名”的對(duì)話框，再次選擇以上生成的*.dat格式文件并點(diǎn)擊“打開”，即可將測(cè)線的線號(hào)展到地形地質(zhì)圖上。具體操作流程圖如圖2所示。

圖2 CASS 8.0繪制工作部署圖流程Fig.2 Flowchart of drawing work deployment map by CASS 8.0

2 應(yīng)用實(shí)例

2.1 1∶1萬高精度磁測(cè)工作中坐標(biāo)點(diǎn)生成及工作部署

在內(nèi)蒙古自治區(qū)額濟(jì)納旗黑鷹山南區(qū)鉛多金屬礦普查項(xiàng)目中，勘查區(qū)位于內(nèi)蒙古自治區(qū)西部額濟(jì)納旗的達(dá)來呼布鎮(zhèn)西偏北約200 km處，地理極值坐標(biāo)為：東經(jīng)98°51′56″—98°59′36″，北緯42°01′31″—42°04′01″，勘查區(qū)面積20.14 km2。

勘查過程中，首先在該區(qū)進(jìn)行了1∶1萬地質(zhì)草測(cè)，1∶1萬高精度磁法測(cè)量，然后選擇位于F1、F2斷層附近的成礦有利部位，進(jìn)行了1∶1萬激電測(cè)量、槽探、激電測(cè)深工作。

(1)GPS衛(wèi)星定位儀參數(shù)的設(shè)置

磁測(cè)掃面測(cè)網(wǎng)以本地區(qū)國家三角點(diǎn)為起算點(diǎn)，施測(cè)GPS(E級(jí))控制網(wǎng)作為該工區(qū)的首級(jí)控制。

測(cè)網(wǎng)基線在首級(jí)控制網(wǎng)的基礎(chǔ)上發(fā)展，采用動(dòng)態(tài)GPS(RTK)進(jìn)行，其相應(yīng)的參數(shù)設(shè)置如下：

坐標(biāo)系統(tǒng)：選擇西安80坐標(biāo)系，高程采用1985高程系。

(2)測(cè)網(wǎng)布置

測(cè)區(qū)內(nèi)測(cè)線按照0°方位角布置。

點(diǎn)、線號(hào)的編排方法原則是：測(cè)線號(hào)與測(cè)點(diǎn)號(hào)從南到、北從西向東由小至大。

(3)測(cè)量網(wǎng)度(100×40)

每條測(cè)線的測(cè)點(diǎn)插衛(wèi)生筷子標(biāo)記，并用記號(hào)筆清楚寫明點(diǎn)線號(hào)，以便后續(xù)工作時(shí)能在野外現(xiàn)場(chǎng)識(shí)別。測(cè)地工作依據(jù)執(zhí)行GB/T18314-2009《全球定位系統(tǒng)(GPS)測(cè)量規(guī)范》、GB/T18314-2001《地質(zhì)礦產(chǎn)勘查測(cè)量規(guī)范》、DZ/T0153-95《物化探工程測(cè)量規(guī)范》、《物化探測(cè)量規(guī)范》，嚴(yán)格遵守有關(guān)規(guī)范和設(shè)計(jì)要求。對(duì)外業(yè)測(cè)量采用不定時(shí)的自檢、互檢、確保測(cè)量質(zhì)量。為保證點(diǎn)位相對(duì)誤差小于5 m，工作中衛(wèi)星數(shù)量不少于5顆。

普查區(qū)高精度磁測(cè)測(cè)線方位角0°點(diǎn)、線號(hào)的編排方法原則是：測(cè)線號(hào)與測(cè)點(diǎn)號(hào)從南到北、從西向東由小至大。測(cè)網(wǎng)網(wǎng)度100 m×40 m，每條測(cè)線的測(cè)點(diǎn)插衛(wèi)生筷子標(biāo)記，并用記號(hào)筆清楚寫明點(diǎn)線號(hào)，以便后續(xù)工作時(shí)能在野外現(xiàn)場(chǎng)識(shí)別。該工區(qū)邊界4個(gè)點(diǎn)的坐標(biāo)分別為：

4657687.170，488859.530

4657685.430，499360.740

4654445.420，488860.890

4654445.280，499361.040

按照1.1的方法得到的前10個(gè)點(diǎn)的坐標(biāo)如下：

橫坐標(biāo) 縱坐標(biāo) 線號(hào) 點(diǎn)號(hào)

4654445.42, 488860.89, 1, 1

4654485.96, 488860.61, 1, 2

4654524.64, 488861.24, 1, 3

4654565.46, 488860.98, 1, 4

4654605.28, 488860.93, 1, 5

4654645.43, 488861.03, 1, 6

4654685.39, 488861.1, 1, 7

4654725.23, 488861.01, 1, 8

4654765.97, 488861.09, 1, 9

4654805.61, 488861.25, 1, 10

將該方法產(chǎn)生的坐標(biāo)文件生成*.csv格式后導(dǎo)入GPS高精度手持儀器(合眾思?jí)袽G838)即可進(jìn)行野外定點(diǎn)。該方法的優(yōu)點(diǎn)是能利用CASS 8.0的專業(yè)性繪圖工具直觀、精準(zhǔn)、簡(jiǎn)潔地繪制工區(qū)坐標(biāo)點(diǎn)的位置。在勘探工區(qū)中若存在不可逾越的障礙物時(shí)，如河流、湖泊、陡崖等，亦可根據(jù)地形地質(zhì)圖及時(shí)調(diào)整測(cè)點(diǎn)位置，合理規(guī)劃路線，提高工作效率，降低野外工作風(fēng)險(xiǎn)，迅速準(zhǔn)確輸入輸出測(cè)點(diǎn)坐標(biāo)。

將以上文件轉(zhuǎn)變?yōu)镃ASS 8.0使用的*.dat文件如圖3所示。

圖3 1∶1萬磁測(cè)測(cè)點(diǎn)轉(zhuǎn)換后的dat文件格式Fig.3 Dat file of 1∶10000 magnetic measurement points after transformation

將以上文件展到CASS打開的地形地質(zhì)圖中，調(diào)整坐標(biāo)點(diǎn)的格式、大小等特性，便可得到相應(yīng)的工作部署圖，如圖4所示。本次應(yīng)用共生成測(cè)線106條，測(cè)點(diǎn)坐標(biāo)8693個(gè)，使用天寶R4 GPS差分定位系統(tǒng)復(fù)測(cè)，精確度達(dá)到mm級(jí)別。

2.2 1∶1萬高精度激電中梯工作中坐標(biāo)點(diǎn)生成及工作部署

普查區(qū)激電區(qū)測(cè)線方位角45°，測(cè)網(wǎng)網(wǎng)度100 m×20 m，線號(hào)、點(diǎn)號(hào)均以1 m為單位編號(hào)，線號(hào)由南向北、點(diǎn)號(hào)由西向東增大。每條測(cè)線的測(cè)點(diǎn)插衛(wèi)生筷子標(biāo)記，并用記號(hào)筆清楚寫明點(diǎn)線號(hào)，以便后續(xù)工作時(shí)能在野外現(xiàn)場(chǎng)識(shí)別。該工區(qū)邊界6個(gè)點(diǎn)的坐標(biāo)分別為：

圖4 黑鷹山南區(qū)鉛多金屬礦普查1∶1萬磁測(cè)測(cè)點(diǎn)坐標(biāo)投影到地圖上效果Fig.4 Map projection of coordinates of 1∶10000 magnetic survey points for lead polymetallic ore survey in the south of Heiyingshan district

4656292.976，490727.507

4656294.945，494406.776

4657708.181，494406.058

4657708.181，494406.058

4654455.492，495071.401

4654455.191，490729.890

按照1.1的方法得到的前10個(gè)點(diǎn)的坐標(biāo)如下：

橫坐標(biāo)， 縱坐標(biāo)， 線號(hào)， 點(diǎn)號(hào)

4655799.580， 490730.022, 1300.000， 1

4655813.722， 490744.164, 1300.000， 2

4655827.864， 490758.306， 1300.000， 3

4655842.006， 490772.448， 1300.000， 4

4655856.148， 490786.590， 1300.000， 5

4655870.290， 490800.732， 1300.000， 6

4655884.433， 490814.875， 1300.000， 7

4655898.575， 490829.017， 1300.000， 8

4655912.717， 490843.159， 1300.000， 9

4655926.859， 490857.301， 1300.000， 10

將該方法產(chǎn)生的坐標(biāo)文件生成*.csv格式后導(dǎo)入GPS高精度手持儀器(合眾思?jí)袽G838)即可進(jìn)行野外定點(diǎn)，將以上文件轉(zhuǎn)變?yōu)镃ASS 8.0使用的*.dat文件,如圖5所示。

圖5 1∶1萬電測(cè)測(cè)點(diǎn)轉(zhuǎn)換后的dat文件格式Fig.5 dat file format of 1∶10000 electric measuring points after transformation

將以上文件展到CASS打開的地形地質(zhì)圖中，

調(diào)整坐標(biāo)點(diǎn)的格式、大小等特性，便可得到相應(yīng)的工作部署圖，如圖6所示。本次應(yīng)用共生成測(cè)線44條，測(cè)點(diǎn)坐標(biāo)4511個(gè)，使用天寶R4 GPS差分定位系統(tǒng)復(fù)測(cè)，精確度達(dá)到mm級(jí)別。

2.3 工作部署圖生成后測(cè)網(wǎng)敷設(shè)及質(zhì)量評(píng)述

2.3.1 測(cè)網(wǎng)敷設(shè)及質(zhì)量評(píng)述

物探網(wǎng)測(cè)線放樣全部采用科力達(dá)風(fēng)云K-9雙頻接收機(jī)進(jìn)行實(shí)時(shí)定位放樣，投入一個(gè)基準(zhǔn)站、四個(gè)移動(dòng)站。這樣既保證了工作效率的提高，又保證了測(cè)量精度。

動(dòng)態(tài)RTK測(cè)量精度指標(biāo)：

平面精度：±2 cm+2 ppm

高程精度：±3 cm+1 ppm

儀器通過山東省計(jì)量科學(xué)研究院質(zhì)量檢測(cè)，其精度指標(biāo)符合作業(yè)要求，可以投入生產(chǎn)。

(1)作業(yè)方法

本次物探網(wǎng)的布設(shè)采用網(wǎng)絡(luò)RTK模式，即在工區(qū)內(nèi)只架設(shè)一個(gè)基準(zhǔn)站，通過GPRS網(wǎng)絡(luò)與通訊控制中心連接，流動(dòng)站也通過GPRS將自己的初始化位置傳輸給控制中心，并接受中心發(fā)送的差分信號(hào)，生成厘米級(jí)的位置信息。網(wǎng)絡(luò)RTK的作業(yè)半徑可達(dá)20 km，能覆蓋整個(gè)工區(qū)。

圖6 黑鷹山南區(qū)鉛多金屬礦普查1∶1萬電測(cè)測(cè)點(diǎn)坐標(biāo)投影到地圖上效果Fig.6 Map projection of coordinates of 1∶10000 electric survey points for lead polymetallic ore survey in the south of Heiyingshan district

選取通視開闊，設(shè)置相對(duì)制高點(diǎn)作為基準(zhǔn)站(基站架設(shè)在未知點(diǎn)上)，在其上安置一臺(tái)接收機(jī)，連續(xù)跟蹤所有的可見衛(wèi)星，連接通訊控制中心，發(fā)送實(shí)時(shí)改正信息；另四臺(tái)有接收GPRS網(wǎng)絡(luò)信號(hào)的GPS接收機(jī)作為流動(dòng)站，流動(dòng)站的移動(dòng)手簿通過在已知控制點(diǎn)上進(jìn)行點(diǎn)校正，校正后便可在事先輸入好的測(cè)線上進(jìn)行測(cè)點(diǎn)的布設(shè)。在測(cè)點(diǎn)位置，壓有寫有點(diǎn)線號(hào)的紅紙條以作標(biāo)記，并且每到達(dá)一個(gè)點(diǎn)都在電子手薄中記錄下當(dāng)前的坐標(biāo)位置。

(2)RTK作業(yè)要求

①至少同步觀測(cè)4顆以上分布良好的衛(wèi)星，并在運(yùn)動(dòng)過程中保持連續(xù)的跟蹤。

②在實(shí)時(shí)動(dòng)態(tài)狀態(tài)下顯示移動(dòng)站的三維坐標(biāo)，以及其相應(yīng)的坐標(biāo)分量的中誤差，通過設(shè)置最大點(diǎn)位中誤差實(shí)現(xiàn)對(duì)測(cè)量精度的限制，符合要求的記錄，否則拒絕記錄。

2.3.2 高精度磁測(cè)工作及其質(zhì)量評(píng)述

普查區(qū)高精度磁測(cè)測(cè)線方位角0°點(diǎn)、線號(hào)的編排方法原則是：測(cè)線號(hào)與測(cè)點(diǎn)號(hào)從南到北、從西向東由小至大。測(cè)網(wǎng)網(wǎng)度100 m×40 m，每條測(cè)線的測(cè)點(diǎn)插衛(wèi)生筷子標(biāo)記，并用記號(hào)筆清楚寫明點(diǎn)線號(hào)，以便后續(xù)工作時(shí)能在野外現(xiàn)場(chǎng)識(shí)別。

磁法測(cè)量使用5臺(tái)GSM19T型質(zhì)子磁力儀。在開工前和工作結(jié)束后對(duì)儀器性儀器一致性進(jìn)行校驗(yàn)。

測(cè)點(diǎn)使用GSM19T型質(zhì)子磁力儀作總場(chǎng)觀測(cè)，每次觀測(cè)時(shí)探頭的高度保持一致，日變和野外磁測(cè)儀器時(shí)間調(diào)成一致。觀測(cè)人員嚴(yán)格“去磁”，未攜帶任何鐵磁性物品。在生產(chǎn)期間，避開了干擾。通過每天對(duì)日變數(shù)據(jù)的分析，未遇到磁暴。變觀測(cè)點(diǎn)位于磁測(cè)總基點(diǎn)上。

日變觀測(cè)使用了性能較穩(wěn)定的GSM19T質(zhì)子磁力儀，進(jìn)行日變觀測(cè)，每20秒觀測(cè)一次。每天的日變觀測(cè)早于各儀器的早校正點(diǎn)觀測(cè)之前開始，遲于晚校正點(diǎn)觀測(cè)之后結(jié)束。

當(dāng)天原始記錄傳入計(jì)算機(jī)之后進(jìn)行100%自檢驗(yàn)收，并填寫磁法資料日驗(yàn)收記錄表,確保記錄數(shù)據(jù)完整齊全，合格后保存原始觀測(cè)文件。內(nèi)業(yè)計(jì)算主要包括儀器一致性精度計(jì)算、測(cè)點(diǎn)ΔT值計(jì)算、質(zhì)量檢查精度計(jì)算等。

質(zhì)量檢查采用“一同三不同”(同點(diǎn)位、不同儀器、不同日期、不同操作員)的質(zhì)量檢查原則進(jìn)行。本區(qū)1∶10000高精度磁測(cè)工作共完成34.55 km2，全區(qū)掃面共完成物理點(diǎn)8693個(gè)，質(zhì)檢點(diǎn)490個(gè)，檢查率5.6%，野外觀測(cè)精度：±1.3 nT。野外觀測(cè)質(zhì)量滿足《地面高精度磁測(cè)勘查技術(shù)規(guī)程》和設(shè)計(jì)的要求，質(zhì)量合格。

ΔDi—第i點(diǎn)觀測(cè)值與放樣坐標(biāo)的點(diǎn)位差；n—觀測(cè)總點(diǎn)數(shù)，i=1,2,3,……n。

平面點(diǎn)位中誤差：m位(剖)=±0.57 m;

m允(剖)=±6.25 m

高程點(diǎn)位中誤差：m高(剖)=±0.41 m;

m允(剖)=±5 m

從檢查結(jié)果來看，本區(qū)的測(cè)量工作精度可靠，精度指標(biāo)滿足設(shè)計(jì)及規(guī)范要求，為野外物探數(shù)據(jù)的采集提供了可靠的保障。

2.3.3 激電中梯工作及其質(zhì)量評(píng)述

普查區(qū)激電區(qū)測(cè)線方位角45°，測(cè)網(wǎng)網(wǎng)度100 m×20 m，線號(hào)、點(diǎn)號(hào)均以1 m為單位編號(hào)，線號(hào)由南向北、點(diǎn)號(hào)由西向東增大。

使用儀器為重慶地質(zhì)儀器廠生產(chǎn)的DFZ10-2型大功率智能發(fā)射系統(tǒng)和與之配套的DJS-8型數(shù)字直流激電接收機(jī)，供電電源選擇為本田發(fā)電機(jī)ET12000交流發(fā)電機(jī)，配以直流穩(wěn)壓裝置。

在野外施工過程中嚴(yán)格按照《時(shí)間域激發(fā)極化法技術(shù)規(guī)程(DZ/T0070-93)》、(DZ/T0153-2014)《物化探工程測(cè)量規(guī)范》執(zhí)行。

激電異常是明顯高于背景值的ηa變化。為減少背景值本身的波動(dòng)和觀測(cè)誤差給異常圈定帶來的影響，一般通過異常下限來劃定異常，異常下限取決于背景值的大小、穩(wěn)定性和ηa的觀測(cè)精度。在激電中梯所測(cè)量的全部極化率數(shù)據(jù)中，剔除畸變、突變數(shù)據(jù)后去掉大于3.0%的數(shù)據(jù)，將剩余數(shù)據(jù)取算術(shù)平均值作為該測(cè)區(qū)極化率背景值。

經(jīng)計(jì)算本區(qū)極化率背景值為1.45%。根據(jù)公式ηa下限=(1+3Mηa)×計(jì)算異常下限，其中Mηa=3.34%，ηa背景=1.45%，計(jì)算得ηa下限=1.60%。

在電腦上對(duì)測(cè)區(qū)所有測(cè)線起始點(diǎn)坐標(biāo)和重復(fù)觀測(cè)點(diǎn)坐標(biāo)等資料數(shù)據(jù)進(jìn)行統(tǒng)一整理記錄。由技術(shù)人員進(jìn)行100%自檢及互檢驗(yàn)收，確保記錄數(shù)據(jù)完整齊全。對(duì)已檢查野外數(shù)據(jù)進(jìn)行室內(nèi)綜合處理，保證每個(gè)數(shù)據(jù)采集對(duì)畸變點(diǎn)進(jìn)行分析并剔除，保證處理結(jié)果的高保真度和電性分布特征的客觀真實(shí)性。結(jié)合數(shù)據(jù)整理激電中梯視電阻率和視極化率等值線平面圖和激電測(cè)深成果圖件，根據(jù)所得資料對(duì)工作區(qū)電性特征進(jìn)行分析解釋，并編寫物探成果報(bào)告點(diǎn)與實(shí)際坐標(biāo)相對(duì)應(yīng)，計(jì)算視電阻率數(shù)值，根據(jù)實(shí)際數(shù)據(jù)采集點(diǎn)位整理物探工作實(shí)際材料圖。

本區(qū)1∶10000激電中梯工作共完成9.00 km2，共計(jì)4644個(gè)有效物理點(diǎn)。質(zhì)量檢查按照“同點(diǎn)位、不同接收機(jī)、不同操作員、不同電流”的“一同三不同”原則進(jìn)行，共計(jì)檢查150個(gè)點(diǎn)，質(zhì)檢比例為3.1%(設(shè)計(jì)及規(guī)范要求不低于3.0%)。經(jīng)統(tǒng)計(jì)計(jì)算，質(zhì)量檢查視極化率實(shí)達(dá)均方相對(duì)誤差為±3.34%(設(shè)計(jì)要求低于4.0%)；視電阻率實(shí)達(dá)均方相對(duì)誤差為±4.40%(設(shè)計(jì)要求低于7.0%)，表明質(zhì)量滿足設(shè)計(jì)及規(guī)范要求，整體工作完成情況良好。

3 結(jié)論

基于對(duì)工作中測(cè)線點(diǎn)坐標(biāo)批量提取及輸入的需求，本文提出了一種運(yùn)用CASS 8.0更加簡(jiǎn)便得出工作部署圖以及測(cè)線點(diǎn)坐標(biāo)的方法，并批量輸入到GPS手持定位儀中。對(duì)比結(jié)合VB6.0語言或Visual C++編程等得出坐標(biāo)點(diǎn)數(shù)據(jù)的方法，此操作方式更加簡(jiǎn)單、直觀。通過直接繪制點(diǎn)和提取點(diǎn)坐標(biāo)，無需坐標(biāo)轉(zhuǎn)換即可獲得能夠批量導(dǎo)入GPS手持儀中的數(shù)據(jù)。同時(shí)，在繪制點(diǎn)陣的過程中，運(yùn)用CASS 8.0帶來的優(yōu)勢(shì)，可以在施工前結(jié)合地形地質(zhì)圖合理刪除一些無效點(diǎn)，甚至可以繪制曲線型的測(cè)線來規(guī)避危險(xiǎn)區(qū)域或者障礙物(懸崖、峭壁等)，以提高工作效率，保障工作人員的生命安全。將這些技術(shù)運(yùn)用到內(nèi)蒙古自治區(qū)額濟(jì)納旗黑鷹山南區(qū)鉛多金屬礦普查工作中，得到了較好的效果。

[附中文參考文獻(xiàn)]

艾斯卡爾，李華.2010.關(guān)于磁法勘探測(cè)點(diǎn)坐標(biāo)的生成與工作部署圖繪制[J].物探與化探，34(2):250-252.

包世泰，夏斌，蔣鵬，付文生.2004.基于GIS的地質(zhì)勘察信息系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)[J].地理與地理信息科學(xué)，20(4):31-35.

崔東鄆，王偉建.2009.激電中梯測(cè)量工作在多金屬找礦中的應(yīng)用[J].貴州地質(zhì)，26(3):218-220.

狄廣禮，劉劍英.2015.基于GPS-RTK技術(shù)在地質(zhì)勘探工程測(cè)量中的應(yīng)用研究[J].科技資訊，13(5):51-52.

管志寧.2005.地磁場(chǎng)與磁力勘探[M].北京:地質(zhì)出版社:1-309.

郭世強(qiáng).2013.淺談在地質(zhì)找礦中地質(zhì)勘探的應(yīng)用[J].科技創(chuàng)業(yè)家，(6):198,200.

何輝，薛典軍，王林飛.2014.地球物理數(shù)據(jù)處理解釋系統(tǒng)中坐標(biāo)轉(zhuǎn)換的設(shè)計(jì)及實(shí)現(xiàn)[J].物探與化探，38(5):1060-1063.

胡志文，歐陽燕，羅湘.2012.水工環(huán)地質(zhì)勘察及遙感技術(shù)在地質(zhì)工作中的應(yīng)用[J].江西建材，(5):187-188.

康世英，張宏偉.2006.GPS測(cè)量在石油物探中的應(yīng)用探討[J].礦產(chǎn)與地質(zhì)，20(4-5):552-555.

李榮亮，田建榮，劉洋，任杰，白順寶，劉曉峰.2017.綜合物探方法在甘肅梧桐井鐵銅多金屬礦勘查中的應(yīng)用[J].地質(zhì)與勘探，53(4):755-764.

李丙法.2017.CASS 8.0應(yīng)用于大比例尺地形圖繪制及使用技巧[J].工程技術(shù)，(3):307.

李得基，梅貴福，田淑梅.2006.“南方測(cè)繪軟件CASS6.1”在工程測(cè)量中的應(yīng)用[J].青海國土經(jīng)略，(2):41-43.

李飛，丁志強(qiáng)，崔志強(qiáng)，胥值禮，李軍峰.2018.CH-3無人機(jī)航磁測(cè)量系統(tǒng)在我國新疆不同地形區(qū)的應(yīng)用示范[J].地質(zhì)與勘探,54(4):735-746.

李俊杰，張紅綱，王偉，何建設(shè)，李劍強(qiáng)，郭佳豪.2019.綜合物探技術(shù)在灰?guī)r地區(qū)隧洞超前預(yù)報(bào)中的應(yīng)用[J].地質(zhì)與勘探,55(6):1452-1462.

林斌.2018.關(guān)于找礦地質(zhì)勘探技術(shù)創(chuàng)新的思考[J].世界有色金屬，(3):84-85.

劉勵(lì)慎,汪津華.1987.PC-1500微型計(jì)算機(jī)處理鉆孔在空間坐標(biāo)位置的實(shí)用方法[J].地質(zhì)與勘探23(12):41-45.

任為望，宋子萱，耿博.2012.關(guān)于找礦地質(zhì)勘探技術(shù)創(chuàng)新的探討[J].科學(xué)技術(shù)創(chuàng)新，(12):17.

邰文星，金尚剛，何彥南，楊梅珍，劉萬龍.2019.綜合物探方法在湖北大冶雞冠咀銅金礦深部及外圍找礦中的應(yīng)用[J].地質(zhì)與勘探,55(4):1010-1025.

王寶江.2010.VB6.0物探線平面坐標(biāo)推算程序設(shè)計(jì)[J].山西建筑，36(9):358-359.

王浩鋒，劉波，陳霜，陳鳳興，薛文浩.2019.綜合物探測(cè)量在二連盆地沉積巖型鈾礦勘查中應(yīng)用[J].地質(zhì)與勘探,55(1):127-134.

王猛，梁中朋.2015.大功率激電在內(nèi)蒙古多金屬礦勘查的應(yīng)用[J].工程地球物理學(xué)報(bào)，12(2):176-182.

王曉軍，陳孝強(qiáng)，傅揚(yáng).2013.草原覆蓋區(qū)綜合物探找礦工作方法—以內(nèi)蒙古阿吉勒銀鉛鋅多金屬礦區(qū)為例[J].高校地質(zhì)學(xué)報(bào)，19(Z):376.

韋晶珠.2010.淺談如何利用CASS軟件繪制地形圖和斷面圖[J].現(xiàn)代企業(yè)文化，(8):119-120.

徐力峰.2014.物化探測(cè)點(diǎn)批量導(dǎo)入手持GPS的方法[J].地球，(1):141-142.

苑立清，趙金凱.2009.地質(zhì)勘探安全生產(chǎn)工作的思考與對(duì)策[J].安全與環(huán)境工程，26(3):87-91.

臧忠淑,蔣樹一.2003.ArcView GIS中圖件比例系數(shù)及坐標(biāo)位置的設(shè)置技巧[J].地質(zhì)與勘探,39(1):58-61.

曾群望.1985.繪制幾種地質(zhì)坐標(biāo)圖的BASIC程序―適用于SHARP PC-1500袖珍機(jī)[J].地質(zhì)與勘探，21(9):32-35.

張建恩，劉亞寧，李康虎，賈傳杰.2010.RTK物探測(cè)量放樣坐標(biāo)轉(zhuǎn)換參數(shù)求取方法的探討[J].物探裝備，20(4):264-268.

張維寬，張永勝，孟紅.2012.全站儀極坐標(biāo)法測(cè)設(shè)物探網(wǎng)的一種新方法[J].價(jià)值工程，(11):36-37.

周朝憲，董少波，和志軍，姜高珍，燕丹晨.2015.自定義坐標(biāo)系的建立及其坐標(biāo)變換實(shí)現(xiàn)[J].地質(zhì)與勘探,51(4):699-704.