本研究利用重力異常資料與地層層位資料來建構苗栗地區地下地質構造。在重力修正的部份，過高或過低的修正密度將導入地形所引起之修正誤差，本研究採用N-密度測定法、g-H關係圖法、布蓋密度與地形密度法等數值方法決定修正密度，並比較其優劣與適用性，以提高布蓋修正與地形修正之精度。
　　本研究也利用地質統計Kriging推估法來提高重力資料準確性。由於Kriging法具有最小變異數與不偏推估的特性，希望藉由該方法來降低網格插值過程所造成的誤差，並解決部份地區重力測點稀少、繁複或重疊的情形。此外，由於重力測勘所求得之地下構造並不具唯一性，本研究亦利用Kriging法來進行苗栗地區地層推估，選擇幾條垂直主要構造的剖面，利用地層資料作為控制，來模擬地下岩層之側向密度變化，並與重力測勘結果相互比對，以期獲得既能反應重力異常分佈亦有實際地層層位控制之地下構造。
　　研究結果顯示，將重力測點資料庫建立後，可隨時取用任意範圍的資料，並配合數值地形資料計算出該範圍的岩層平均密度，不但省去直接測量密度的麻煩，所求得之修正密度較具代表性。此外，透過剖面的模擬，顯示在研究區西北側的地層有些微向上抬升，這可能與前人認為台灣西北部中新世地層表層滑動及基底逆衝有關。至於在研究區的東側，整個沈積地層受到麓山帶邊緣造山運動的影響，構造呈現上揚並產生不少高角度的斷層。利用地層資料模擬的結果與重力觀測的結果大致相同，只有在山區無法吻合。整體上理論重力值與觀測重力值的差值皆在 ±5 mgal以內。

　　It is well known that the use of an incorrect density in constructing Bouguer gravity map will cause abnormal results. To reconstruct a Bouguer gravity map in the Miao-Li area, average rock density was investigated in this study. For optimization of gravity result, the density for the Bouguer correction and terrain correction of gravity was derived from N-method, g-H relations, and the method of Bouguer density and terrain density in this study.
　　In addition, gravity data are often acquired over long periods of time using different instruments and various techniques, resulting in data sets of non-uniform accuracy. In order to obtain more precise gravity data, the Kriging technique of geological statistics was employed to determine Bouguer anomaly value exactly. With the feature of least variance and the best linear unbiased estimator, Kriging technique could take into account station accuracy and decrease interpolation error. The Kriging technique was also applied to the simulation of geological structure.
　　The general procedure to determine geological structure is as fallows. Density structures of three cross-sections perpendicular to strike of the surface geology were obtained from modeling of well strata. We compared calculation result from modeling with observed values, which were judged by traditional correction of gravity. Moreover, the density structures of the model were fixed until the anomaly curves trend to coincide each other.
　　After processing, however, the new method for determination of density has two advantages. First, since the density was computed with the function of the latitude, longitude, elevation, and observed gravity value of the station, the new method can save more time and has smaller error than the previously direct measured method. Second, the density of any region can be obtained form database immediately without reconstructing station data. Besides, the surface density structures obtained from modeling of three cross-sections indicated that the formations in west-northwestern side are uplift. This feature can be recognized by a step in decollement of a deep thrust in the base of Miocene formations. The fault-bend folding mountain belt in eastern side shows that deformation of the Miao-Li area was caused by the northwestward compressive force which was possible the results of the orogeny of early Pleistocene. In general, the gravity results between observed and modeling are consonance except for mountain area.

1. Chai Y. and Hinze W. J., 1988, Gravity inversion of an interface above which the density contrast varies exponentially with depth, Geophysics, Vol. 53, pp. 837-845.
2. Cheeney R. F., 1983, Statistical methods in geology for field and lab decisions, George Allen & Unwin Ltd, London, 169 p.
3. Cordell L., 1973, Gravity analysis using an exponential density depth function, San Jacinto graben, California, Geophysics, Vol. 38, pp. 684-690.
4. Davis J. C., 1973, Statistics and data analysis in geology, John Wiley & Sons, New York, pp. 239-248, 383-405.
5. Dehlincer P., 1978, Marine gravity, Elsevier Sci. Pub., Netherlands, 137 p.
6. Edward H. I. and Srivastava R. M., 1989, Applied geostatistics, Oxford University Press, New York, 561 p.
7. Edzer J. P. and Cees G. W., 1998, Gstat: a program for geostatistical modelling, prediction and simulation, Computers & Geosciences, Vol. 24, No. 1, pp. 17-31.
8. Edzer J. P., 1999, Gstat user''s manual, Internet, Available: http://www.frw.ruu.nl/gstat/manual/gstat.html.
9. Englund E. and Sparks A., 1991, GEO-EAS 1.2.1 user''s guide, Environmental Monitoring Systems Laboratory Press, U.S. Environmental Protection Agency, Las Vegas, 186 p.
10. Environmental Modeling Research Laboratory, 1999, Groundwater modeling system tutorials, Brigham Young University Press, pp. 2-1~3-18.
11. Fukao Y., Yamamoto A. and Nozaki K., 1981, A method of density determination for gravity correction, Journal of Physical Earth, Vol. 29, pp. 163-166.
12. Garcia-Abdeslem J., 1992, Gravitational attraction of a rectangular prism with depth-dependent density, Geophysics, Vol. 57, pp. 470-473.
13. Grant F. S. and West G. F., 1965, Interpretation theory in applied geophysics, McGraw-Hill, New York.
14. Hagiwara Y., 1967, Analysis of gravity values in Japan, Bull. Earthq. Res. Inst., Tokyo University, Vol. 45, pp. 1091-1228.
14. Hagiwara Y., 1967, Analysis of gravity values in Japan, Bull. Earthq. Res. Inst., Tokyo University, Vol. 45, pp. 1091-1228.
16. Hammer S., 1939, Terrain corrections for gravity stations, Geophysics, Vol. 4, pp. 184-194.
17. Howell B. F., 1971, Introduction to geophysics, 214 p.
17. Howell B. F., 1971, Introduction to geophysics, 214 p.
19. Jeffreys H., 1952, The earth, 3rd ed., Cambridge University Press.
19. Jeffreys H., 1952, The earth, 3rd ed., Cambridge University Press.
21. Keating P., 1995, Error estimation and optimization of gravity surveys, Geophysical Prospecting, Vol. 43, pp. 569-580.
22. Keckler D., 1995, Surfer 6 user''s guide, Golden Software, Inc., pp. 5-32-5-58.
23. Marsal D. and Daniel F. M., 1987, Statistics for geoscientists, Pergamon Books Ltd., New York, 176 p.
23. Marsal D. and Daniel F. M., 1987, Statistics for geoscientists, Pergamon Books Ltd., New York, 176 p.
25. Maus S., 1999, Variogram analysis of magnetic and gravity data, Geophysics, Vol. 64, No. 3, pp. 776-784.
25. Maus S., 1999, Variogram analysis of magnetic and gravity data, Geophysics, Vol. 64, No. 3, pp. 776-784.
27. Nettleton L. L., 1942, Gravity and magnetic calculations, Geophysics, Vol. 7, pp. 293-310.
28. Nettleton L. L., 1976, Gravity and magnetics in oil prospecting, McGraw-Hill, New York, pp. 144-147.
28. Nettleton L. L., 1976, Gravity and magnetics in oil prospecting, McGraw-Hill, New York, pp. 144-147.
30. Rao D. B., 1990, Analysis of gravity anomalies of sedimentary basins by an asymmetrical trapezoidal model with quadratic density function, Geophysics, Vol. 55, pp. 226-231.
30. Rao D. B., 1990, Analysis of gravity anomalies of sedimentary basins by an asymmetrical trapezoidal model with quadratic density function, Geophysics, Vol. 55, pp. 226-231.
30. Rao D. B., 1990, Analysis of gravity anomalies of sedimentary basins by an asymmetrical trapezoidal model with quadratic density function, Geophysics, Vol. 55, pp. 226-231.
33. Telford W. M., Geldart L. P., and Sheriff R. E., 1990, Applied geophysics, Cambridge University Press, New York, pp 6-61.
33. Telford W. M., Geldart L. P., and Sheriff R. E., 1990, Applied geophysics, Cambridge University Press, New York, pp 6-61.
35. Turcotte D. L. and Schubert G., 1983, Geodynamics, 206 p.
36. Wen W. B. and Fan W. C., 1989, Applied geophysics, Chinese Geophysical Society, Beijing, 159 p.
37. Yamamoto A., Nozaki K., Fukao Y., Furumoto M., Shichi R. and Ezaki T., 1982, Gravity survey in the central ranges, Honshu, Japan, Journal of Physical Earth, Vol. 30, pp. 201-243.
37. Yamamoto A., Nozaki K., Fukao Y., Furumoto M., Shichi R. and Ezaki T., 1982, Gravity survey in the central ranges, Honshu, Japan, Journal of Physical Earth, Vol. 30, pp. 201-243.
39. 中國石油學會勘採組，1982，地球物理探勘概論，中國石油學會，共386頁。
40. 王謙身，張赤軍，周文虎，蔣福珍，1995，微重力測量-理論、方法與應用，科學出版社，北京，共154頁。
41. 史天元，邱元宏，1998，台灣地區數值密度模型之建立，測量工程，第40卷，第3期，第3-22頁。
42. 台灣地質圖，1986，經濟部中央地質調查所編印。
43. 江新春，1970，苗栗永和山背斜之震測解釋，台灣石油地質，第7號，第65-72頁。
44. 何春蓀，1986，台灣地質概論－台灣地質圖說明書，經濟部中央地質調查所，共164頁。
45. 吳榮章，梅文威，童英銓，蔡錦樁，林國安，1992，新竹外海沈積中心遷移分析及油氣潛能研究，台灣石油地質，第27號，第103-127頁。
46. 呂崇嘉，1988，三義、銅鑼地區之地球物理探勘，國立中央大學地球物理研究所碩士論文，共93頁。
47. 李維峰，1997，濁水溪沖積扇之地電阻測勘研究，國立中央大學應用地質研究所碩士論文，共74頁。
48. 李錫堤，1996，數值地形及遙測影像在個人電腦上的應用，國立中央大學應用地質研究所講義，共55頁。
49. 邱華燈，1967，台灣西北部地下地層對比，台灣石油地質，第5號，第33-44頁。
50. 邱華燈，1968，台灣西北部苗栗細道邦附近之地質構造，台灣石油地質，第6號，第33-43頁。
51. 胡錦城，邱華燈，1984，台灣西北部卓蘭附近之深部地質構造，台灣石油地質，第20號，第21-33頁。
52. 張忠富，1998，利用數值地形模型從事重力之地形修正研究，國立中央大學地球物理研究所碩士論文，共75頁。
53. 張錫齡，1974，台灣苗栗區域地下地質研究，台灣石油地質，第11號，第1-25頁
54. 張徽正，林啟文，陳勉銘，盧詩丁，1998，臺灣活動斷層概論－臺灣活動斷層分析圖說明書，經濟部中央地質調查所，共103頁。
55. 陳培源，1992，野外及礦業地質學，正中書局，第381-392頁。
55. 陳培源，1992，野外及礦業地質學，正中書局，第381-392頁。
57. 楊潔豪，1987，苗栗縣三義及銅鑼地區地球物理探勘，經濟部中央地質調查所專題研究成果報告，共157頁。
58. 劉發祥，蔡清吉，1999，苗栗縣細道邦構造重力補測(細道邦-烏石坑)報告，中國石油股份有限公司台灣油礦探勘總處，未發表。
59. 謝世雄，1964，台灣苗栗出磺坑油田出磺坑層內儲油層之研究，台灣石油地質，第3號，第75-92頁。
60. 謝世雄，1970，台灣苗栗出磺坑構造南段地下地質之研究，台灣石油地質，第7號，第103-119頁。
61. 謝世雄，胡錦城，1972，台灣重力與磁力之研究，台灣石油地質，第10號，第283-322頁。
62. 藍生傑，1981，台灣西部麓山帶之構造：苗栗及新竹區之南部，台灣石油地質，第18號，第31-51頁。
63. 藍生傑，1983，台灣西部麓山帶之構造：苗栗及新竹區之中央部份，台灣石油地質，第19號，第51-76頁。
64. 顏宏元，1992，台灣地區重力異常分佈及其在地體構造上之含意，國立中央大學地球物理研究所博士論文，共98頁。