VBA在計算重力固體潮理論值中的應用

張永旺 趙麗娟 趙婧 蘇鵬飛

摘要：重力固體潮改正是重力數(shù)據(jù)處理的基礎，在重力工作中重力儀的儀器試驗、重力基點網(wǎng)聯(lián)測、重力測點絕對重力值計算等主要重力數(shù)據(jù)處理中，重力固體潮理論值的計算都必不可少。但重力固體潮的計算公式繁雜，為提高計算效率，簡化計算過程，筆者應用EXCEL中的Visual Basic編輯器，通過VBA代碼實現(xiàn)了重力固體潮理論值的一鍵計算。

關鍵詞：重力固體潮；VBA；重力計算

Abstract： the theory value of Gravity solid earth tides is the basis of gravity data processing， in the gravity test instrument， gravity basepointsnet conjunction， absolute gravity value of observation point and so on， the theory value of Gravity solid earth tides is necessary of gravity data-processing. But Gravity solid earth tides correction is complicated. In order to improve the computational efficiency， simplify the calculation process， the author take full advantage of VBA editor in the EXCEL， through VBA code calculating the theory value of Gravity solid earth tides by one key.

Key words： Gravity solid earth tides；VBA；Gravity Calculation

前言

隨著礦產(chǎn)資源的勘查與開發(fā)日趨向深部發(fā)展，深入了解大地構造，劃分斷裂構造及其控礦特征，圈定隱伏巖體，對于開展深部找礦具有決定性的作用，而重力工作對于解決以上地質(zhì)問題具有明顯優(yōu)于其他物化探方法的優(yōu)勢，所以重力工作對于進一步推動地質(zhì)礦產(chǎn)勘查的轉型和發(fā)展都具有重要意義。

重力工作中數(shù)據(jù)處理的真實準確是開展解釋推斷的基礎，而重力固體潮理論值的計算又是重力數(shù)據(jù)處理的前提，所以準確無誤的獲得重力固體潮的理論值顯得尤為重要。

目前，重力固體潮理論值的計算有幾款商業(yè)軟件可以實現(xiàn)，但不同軟件輸入格式多樣，其計算結果應用于下級計算步驟繁雜，不利于資料的二次利用。通過VBA實現(xiàn)重力固體潮的計算可以簡化重力資料處理，提高工作效率。希望文章對重力工作者有所幫助，同時也歡迎廣大重力工作者批評指正。

1. 概念

1.1 重力固體潮

由于地球和太陽及月亮的相對位置不斷變化而引起的重力變化稱為重力固體潮。

1.2 重力固體潮理論值

月球和太陽對地球上某一點作用值之和即為重力固體潮的理論值。

2. 格式要求

數(shù)據(jù)格式依次為點號、日期（年-月-日）、時間（60進制，h：mm：ss）、經(jīng)度（60進制，DDDMMSS.S）、緯度（60進制，DDMMSS.S）

運行VBA代碼后會在F列顯示點對應的重力固體潮理論值。

3. VBA代碼及解釋說明

4. VBA代碼的應用

實際使用時，新建一個Excel表格，在Sheet1中輸入數(shù)據(jù)，然后在工具選項卡中打開宏-Visual Basic編輯器，在代碼窗口中，粘貼上述代碼，點擊運行（如圖1紅圈所示）按鈕即可在數(shù)據(jù)尾列獲得重力固體潮理論值（單位毫伽）。

5. 算例與分析

為確保VBA代碼的準確性和可靠性，筆者針對性的選擇以下特殊情況的數(shù)據(jù)（見表2）進行驗算，并與“區(qū)域重力信息系統(tǒng)RGIS2006”、“跨平臺金維地學信息處理研究應用系統(tǒng)GeoIPASV3.2”計算成果進行對比。

從表2可知，由于計算過程中截取位數(shù)的差別，筆者編寫的VBA代碼與RGIS2006的計算成果最大偏差0.002毫伽，與GeoIPAS V3.2最大偏差0毫伽，證明在不同特殊情況下，筆者編寫的VBA代碼計算成果準確可靠。經(jīng)過大量檢查計算（此處不一一列舉），VBA代碼的計算成果與RGIS 2006的最大偏差僅0.003毫伽，與GeoIPASV3.2基本沒有偏差。

6. 結束語

利用VBA代碼進行重力固體潮的改正不僅方便快捷，便于進行數(shù)據(jù)的下級計算，而且精度可靠，能夠滿足大比例尺重力工作，具有較高的應用推廣價值，同時歡迎廣大重力工作者批評指正。

參考文獻：

[1] DZ/T0082-2006，區(qū)域重力調(diào)查規(guī)范[S].

[2] 馬玄龍，劉永生.重力固體潮在Excel工作表中的計算方法[A].資源環(huán)境與工程.2010，24（2）：203-206.

[3] 宋偉，吳建國等.中文Visual Basic6.0編程基礎[M].北京：清華大學出版社，1999.