以锥体断面法直接做重力地形改正.PDF

以錐體斷面法直接做重力地形改正 王建立 國立成功大學地球科學系暨研究所 摘 要 重力異常（gravity anomaly）一般定義為大地水準面的實測重力值與正常地球 橢球面正常重力值的差值。但大地水準面的重力值並無法直接測量，所以需將地 表測量的重力值歸算至大地水準面上。重力歸算主要包含二步驟：一為自由空間 改正；一為去除測站與大地水準面層間質量引力的重力地形改正。傳統上重力地 形改正是先作布格改正（將測站附近視為無限延伸的平坦地形），再移去及補上 高出測站或低於水平面地形質量的引力（圖1）。本研究為捨棄傳統方式，採扇 形錐體斷面法直接計算大地水準面上層間質量對測站的引力。 扇形錐體斷面法為將測站視為原點，將地形劃分為等圓心角的扇形，並以各 扇形中心線的斷面（剖面）代表各扇形區塊的地形變化，圓心角愈小表現的地形 也愈真實。扇形區塊依斷面數值（相對測站的距離及高程）再細分為更小扇形圓 錐區塊，再計算及累計各扇形錐體區塊對測站的引力。扇形錐體的內半徑及外半 徑並不等高，並且以大地水準面為基準面直接計算。有別於漢默法將各扇形視為 平坦且以測站水平面為基準面計算方式。 本研究先以圓錐模型測試第三章推導的扇形錐體引力公式（圖2、式1），再 選擇4 點地形變化較大，並能代表不同類型地形特徵，分別為山頂2 點、山谷1 點及山腰1 點進行實際地形測試，以圓心角π/6（相差π/12 建制之斷面資料） 之扇形錐體作測試，誤差最大為青山（山腰）1.576mGal（表1、圖3）。在真實 地形測試方面與傳統漢默法比較，除玉山（山頂）約13.4 mGal，其它3 測試點 誤差皆在4 mGal以內（表2、圖4）。由測試結果得到，在圓錐模型或真實地形 測試方面，錐體斷面算法的引力大於傳統漢默法，這是因漢默法將修正的扇形區 塊視為平坦地形（取平均高程）的關係。 參考書目 李坤哲，2001，重力地形修正之探討，國立中央大學地球物理研究所碩士論文， 桃園，72頁。 陳國華，2004，整合 TWVD2001 水準及 GPS 資料改近台灣區域大地水準面模式已 應用於 GPS 高程測量，國立成功大學測量及空間資訊學系博士論文，台南， 154頁。 蕭宇伸，2002，嚴密正高改正，國立交通大學土木工程研究所碩士論文，新竹， 67頁。 胡明城，2003，現代大地測量學的理論及其應用，測繪出版社，北京，476 頁。 孔祥元、郭际明、劉宗泉，2002，大地測量學基礎，武漢大學出版社，武漢，265 頁。 圖 1 傳統地形改正及直接計算地形改正的斷面示意圖 圖 2 環狀錐體模板對 P點引力示意圖 錐體模板引力公式： ⎡ ( )2 2 ⎤ + + + + 2π δ a p q a p G ⎛ 2 2 2 2 ⎞ i+1 i+1 ⎢ ( + + −) ( + +) − ⋅ln ⎥ FZ ⎜ ai+1 p q ai p q ⎟ p (m 2 +1)⎢⎝ ⎠ ( )2 2