本研究主要是參考美、日兩國之經驗，以地震儀陣列量測微地動資料來推估淺層沈積層之剪力波速度構造之方法。包括從野外施測步驟，反覆驗證後取得測區之微地動資料，到室內進行頻率-波速分析求取表面波頻散曲線，然後執行多次基因演算法逆推得淺層剪力波速度構造，最後以傳統表面波頻散曲線逆推技術，進一步修正此速度模型。
本研究主要施測的地區為台北關渡(關渡大橋北岸右側數百公尺空地上)及921受創嚴重的埔里盆地，在研究中發現近地表之剪力波速度應呈梯度增加之形態；以關渡地區為例，在深度約400-700公尺處應為基盤之所在位置，而埔里盆地地區，則在不同深度(10-25公尺)會碰到約1.0-1.4km/s之高速層(礫石層)。但詳細深度必須有更低頻之訊號才能解析得出。目前結果只能確定該處應有-『高速基盤』，如此計算的頻散曲線才能與陣列觀測的相速度吻合。

Follow Japan and United States seismologists` experience, We developed a method to estimate the near surface sediment`s shear wave velocity structure by using microtremor data which measured from a portable seismic array. This method includes the data acquisition procedures for the field works and the analysis in the laboratory. The analysis procedures can be categorized as (1) Frequency-Wavenumber analysis to get the dispersion curve, (2) to search the good velocity models using genetic algorithm, and (3) to refine the GA searched model using traditional surface wave inversion technique.
In this study, we applied our method to Guendu, Taipei and Puli basin. Puli basin is seriously damaged the 921 chi-chi earthquake. From our analysis, we found that the velocity is gradually increasing for the near surface sediments. The results show that the high speed basement appear in the depth of 400-700m for Guendu, and roughly about 10-25m deep for Puli basin.

Aki, K. and P. G. Richard (1980). Quantitative Seismology: Theory and Methods .
Brune, J. N., J. E. Nafe, and J. Oliver (1960). A simplified method too the analyis and synthesis of dispersed wave trains, J. Geophys. Res. 65, 287-304.
Capon, J., R. J. Greenfield, R. J. Koller, and R. T. LaCoss (1968). Short-period signal processing results for the large aperture seismic array, Geophysics. 33, 452-472.
Capon, J., R. J. Greenfield, and R. T. LaCoss (1969). Long-period signal processing results for the large aperture seismic array, Geophysics.
34, 305-329.
Capon, J. (1974). Characterization of crust and upper mantle structure under LASA as a random medium, Bull. Seism. Soc. Am. 64, No.1, 235-266.
Capon, J. (1969). Investigation of long-period noies at the large aperture seismic array, J. Geophys. Res. 74, 3182-3194.
Chang, T. M. (1997). Evaluation of suface-wave waveform modeling for lith sphere velocity structure, Ph. D Dissertation, Saint Louis University.
Constable, S. C., R. L. Parker and C. G. Constable (1987). Occam’s inversion: A practical algorithm for generating smooth model from
electromagnetic sounding data, Geophysics. 52, 289-300.
Dziewonski, A. M., S. Bloch, and M. Landisman (1969). A technique for the analysis of transiend seismic signal, Bull. Seism. Soc. Am. 54, 427-444.
Goldberg D. E. (1989). Genetic algorithms in search, Optimization and machine learning, Reading, MA:Addison-Wesley.
Herrin, E., and T. Goforth (1977). Phase-matched filters: application to the study of rayleith waves, Bull. Seism. Soc. Am. 67, 1259-1275.
Herrmann, R. B. (1973). Someaspects of band-pass filtering of surface wave, Bull. Seism. Soc. Am. 63, 663-671.
Herrmann, R. B. and Russel(1985). Computer program in seismology, Saint Louis University.
Holland, J. H. (1992). Genetic algorithms, Sci. Am., 66-72, July.
Horike, M.(1985). Inversion of phase velocity of long-period microtremors to the S-wave-velocity strusture down to the basement in urbanized areas, J. Phys. Earth 33, 59-96
Kanai, K. (1962). On the spectrum of strong earthquake motions, Bull. Earthq. Res. Inst.40, 71- 90.
Kanai, K. (1962). On the predominant period of earthquake motions, Bull. Earthq. Res. Inst. 40, 855- 860.
Katz, L. J. (1976). Microtremors analysis of local geological conditions, Bull. Seism. Soc. Am. 66, No.1, 45- 60.
Marquardt, D. W. (1970). Generalized inversion, ridge regression, biased linear estimation and nonlinear estimation, Technometrics. 12, 591-612
Matsushima, T., and H. Okada (1990). Determination of deep geological structures under urban areas using long-periods microtremors, Butsuri Tansa 43, 21-33.
McMechan, G. A., and M. J. Yedlin (1981). Analysis of dispersive wave by wave field transformation, Geophysics, 46,NO. 6, 869- 874, 7 FIGS.
Menke, W. (1989). Geophysical data analysis: Discrete inversion theory, Academic, san diego, 266 pp.
Russell, D. R. (1987). Multi-channel processing of dispersed surface wave, Ph. D Dissertation, Saint Louis University.
Sambridge, M., and G. Drijkoningen(1992). Genetic algorithms in seismic waveform inversion, Geophys. J. Int. 109, 323-324.
Sambridge, M., and K. Gallagher(1993). Earthquake hypocenter location using genetic algorithms, Bull. Seis. Soc. Am. 83 , 1467-1491.
Sasatani, T., K. Yoshida, H. Okada, O. Nakano, T. Kobayashi, and S. Ling(2001). Estimation of deep subsurface structures and observation of strong ground in Sapporo urban districts, J. Jsnds 20-3, 325-342.
Satoh, T., H. Kawase, T. Iwata, S. Higashi, T. Sato, K. Irikura, and H. C. Huang (2001). S-wave veolocity structure of the Taichung Basin, Taiwan, Estimated form Array and single-station records of microtremors, Bull .Seism. Soc. Am. 91, No.5, 1267- 1282.
Satoh, T., H. Kawase, and S. Matsushima (2001). Structures in and around the Sendai Basin, Japan, Using array records of micretremors, Bull. Seism. Soc. Am. 91, No.2, 206- 218.
Satoh, T., H. Kawase, and S. Matsushima (2001). Different between site characteristics obtained from microtremors, S-waves, P-waves, Bull. Seism. Soc. Am. 91, No.2, 313- 334.
Satoh, T., H. Kawase, M. Matsui, and S. Kataoki (1991). ARRAY measurement of high frequency microtremors for underground structure estimation, Proceeding Fonrth International Conference on Seismic Zonation.2, 409- 416.
Twomey, S. (1977). Introduction to the mathematics of inversion in remote sensing and indirect measurement, Elsevier Scientific Pub. Co.
Wang, C. Y. (1981). Wave theory for seismogram synthesis, Ph. D Dissertation, Saint Louis University.
周鵬程 (2001)。遺傳演算法原理與應用。全華出版。
唐昭榮 (2001)。探討埔里地區之場址特性。國立中正大學應用地球物理所碩士論文。
黃文紀 (1986)。羅東強震儀陣列區微地動之來源與特性。國立中央大學地球物理所碩士論文。
黃紹揚 (2000)。頻散曲線之有限差分法反算土層剪力波速。國立台灣大學土木工程學研究所碩士論文。
黃瑞德 (1999)。由表面波頻散推研中國大陸之上部地函構造。立中央大學地球物理所博士論文。
黃鑑水、謝凱旋、陳勉銘 (2000)。五萬分之一埔里地質圖及說明書。經濟部中央地質調查所。