本研究選取強地動觀測計畫(Taiwan Strong Motion Instrumental Program, TSMIP)的集集地震事件序列之垂直分量、南北分量、東西分量等三分量合成之總量記錄，對於收錄到集集主震的四個分區：台中分區(TCU)、花蓮分區(HWA)、嘉南分區(CHY)與台東分區(TTN)進行分析。使用複數解載(CMD)與三維Qp值模型解算出4~10 Hz內各測站之震源加速度函數；先做同一分區的場址效應平均化，再做各個分區間比較。以同樣的方式再做集集地震序列之一餘震各分區之平均震源加速度函數，同樣進行各分區間比較，並進行主震與餘震事件的比較。
研究結果有以下發現：
集集主震與餘震的震源加速度函數，去除異常紀錄之台中測站，各分區的震源加速度時間函數，都有共同的狀態。首先從0至4秒之間僅有微量的水平振幅，約在0.5~1 gal之間，4秒以後開始有所增幅，至第7秒多達到峰值，然後振幅逐漸減低，之後因不同分區而有不同的減低幅度，直至原來前3~4秒時的振幅。50個測站中有16個測站其S-P波走時差在10.24秒以內的花蓮分區、與35個測站中有5個測站其S-P波走時差在10.24秒以內的嘉南分區，震源加速度時間函數最大振幅皆達到3 gal。全部測站S-P走時差皆在10.24秒以上之台東分區，震源加速度時間函數最大約1 gal。
集集餘震的震源加速度函數，因所有分區走時差在6秒以上，且起始與結束部分波動幅度差不多大，這代表了所有分區的震源加速度函數只有P波到達，或是有部分S波滲入，並未像集集主震台中分區之震源加速度時間函數，受到表面波滲入影響，在時間末端的加速度時間函數有極大的振幅。將集集餘震所有分區的震源加速度函數進行比較，會發現此次餘震事件中，台中分區與花蓮分區的測站走時差皆在10.24秒內，花蓮分區走時差中位數約在7.711秒、台中分區走時差中位數約7.716秒，也表示震央距中位數兩者相差不大，但台中分區的震源加速度函數最大振幅卻仍達到±15 gal，花蓮分區的震源加速度函數最大振幅卻僅有±5 gal，由於分區裏所有測站之S-P波走時差皆在6秒以上，可以排除是表面波到達所影響，造成兩個分區震源加速度函數差異甚大原因，推測台中分區與同一事件其他分區相比，場址效應明顯大於同一事件其他分區。而餘震事件中摒除較大場址效應之台中分區不計，其震源加速度時間函數之變化，在0~5秒以內，花蓮分區有3~4 gal、台東分區有2~3 gal、嘉南分區有0.5~1.5 gal的起伏，5秒過後開始出現反向振幅，在7秒時達到峰值，直到10.24秒結束為止；花蓮分區最大振幅可達到5gal、台東分區最大振幅可達到3gal、嘉南分區最大振幅僅有2gal；由以上結果發現，震央距越遠，震源加速度函數之最大振幅越小。

Akaike, H. (1971).Autoregressive model fitting for control, Annals of the Institute of Statistical. Mathematics, 23, 163-180.
Akaike, H. (1979). A bayesian extension of the minimum AIC procedure of autoregressive model fitting, Biometrika, 66, 237-242.
Allen, R.M., Kanamori, H.(2003) The Potential for Earthquake Early Warning in Southern California, Science, 300, 786-789.
Azakawa, T. (2004a). Automatic Detection Procedure of Onset Time of P- and S-Phase in Seismic Records. Architectural Institude of Japan, 21246, 491-492. [Text in Japanese]
Azakawa, T. (2004b). A Technique for Automatic Detection of Onset Time of P- and S-Phase in Strong Motion Records. 13th World Conference on Earthquake Engineering, Symposium held at Vancouver, B. C., Canada.
Banks, R. J. (1975). Complex demodulation of geomagnetic data and the estimation of transfer function. Geophys. J. R. astr. Soc., 43, 87-101.
Brigham, E. O. (1974). The fast fourier transform (2nd ed.). Englewood Cliffs, NJ: Prentice-Hall International.
Brigham, E. O. （1985）。快速傅立葉變換（The Fast Fourier Transform）（黎文明譯）。台南市：復漢出版社。（原作1974年出版）
Chen, K. J. (1998). S-Wave Attenuation Structure in the Taiwan Area and its Correlation to Seismicity. Terrestrial, Atmospheric and Oceanic Sciences, 9(1), 97-118.
Chen, K. J., Yeh, Y. H., & Shyu, C. T. (1996). Qp Structure in the Taiwan Area and Its Correlation to Seismicity. Terrestrial, Atmospheric and Oceanic Sciences, 7(4), 409-429.
Chen, K.J., Yeh, Y. H., & Tsai, Y. B. (1985). A study on complex demodulation of the geomagnetic total intensity data of Taiwan. Bull. Inst. Earth Sci., Academia Sinica, 5, 125-136.
Chen, K. J., Ho, Y. R., Chiu B., Wang, J.S., Sun, R., Lin, C. H. (2007). Correlation Between Time Change in Modulus of Short-Period Geomagnetic Variation and Seismicity in Taiwan. Terrestrial, Atmospheric and Oceanic Sciences, 18(3), 577-591.
Chopra, S., Kumar, D., & Rastogi, B. K. (2011). Attenuation of high frequency P and S waves in the Gujarat region, India. Pure and applied geophysics, 168(5), 797-813.
Geo SIG Presentation Seismic Signals and Sensors, Retrieved August 07, 2014, from
http://www.geosig.com/files/GS_Presentation_Seismic_Signals_and_Sensors.pps
Hough, S. E., Anderson, J. G., Brune, J., Vernon III. F., Berger, J., Fletcher, J., Haar, L., Hanks, T., & Baker L. (1988). Attenuation near Anza, California. Bulletin of the Seismological Society of America, 78(2), 672-691.
Kanamori, H. (1967). Spectrum of Short-period core phases in relation to the attenuation in the mantle. Journal of Geophysical Research, 72(8), 2181-2186.
Maeda, N., (1985). A method for reading and checking phase times in auto-processing system of seismic wave data. Zisin Ser. 2, 38(3), 365-379. [Text in Japanese]
Ma, K. F., Wang, J. H., Zhao, D. (1996). Three-Dimensional Seismic Velocity Structure of the Crust and Uppermost Mantle beneath Taiwan. Journal of Physics of the Earth, 44(2), 85-105.
Rau, R. J.,Wu. Y. M. (1995). Tomographic imaging of lithospheric structures under Taiwan. Earth and Planetary Science Letters, 133(3-4), 517-532.
Roecker, S. W., Yeh, Y. H., & Tsai, Y. B. (1987). Three-dimensional P and S wave velocity structures beneath Taiwan: deep structure beneath an arc-continent collision. Journal of Geophysical Research, 92(B10), 10547-10570.
Sleeman, R., Torild van Eck (1999). Robust automatic P-phase picking: an on-line implementation in the analysis of broadband seismogram recordings. Physics of the Earth and Planetary Interiors, 113, 265-275.
WPSMap - Purchase-Build a Seismometer, Retrived Novenber 23, from http://www.wpsmap.com/menu/seismometer.php
Waves, Seismometers, and Seismograms, Retrived August 07, 2017, from http://eqseis.geosc.psu.edu/~cammon/HTML/Classes/IntroQuakes/Notes/seismometers.html
Yoshida, M., Chiba, H., Hagiwara, H., Ogata, Y., and Mizoue, M. (1994). Summary of Real Time Processing Methods for Micro-earthquake Automatic Processing System (1981-1993). Bull. Earthq. Res. Inst., 69, 121-138. [Text in Japanese]
工程地震學研究室(2017)。工程地震學研究室。 國立中正大學地球與環境科學系。上網日期：2017年9月，檢自： http://www.eq.ccu.edu.tw/lab/lab520/introduction.html
中央氣象局地球物理資料管理系統。自由場強地動觀測網。上網日期：2014年9月，檢自：
http://gdms.cwb.gov.tw/gdms-freefield.php
中央氣象局地震測報中心(2009)。中央氣象局地震測報中心20週年專刊. 中央氣象局地震測報中心. Retrieved from http://scweb.cwb.gov.tw/Twenty.aspx?ItemId=8&loc=tw
王仁聖(2011)。以時頻分析探討地動加速度譜之異常。國立臺灣師範大學 地球科學系碩士論文。
地震活動彙整，中央氣象局。上網日期：2014年9月，檢自
http：//www.cwb.gov.tw/V7/earthquake/rtd_eq.htm
林孝維(2004)。利用傅氏振幅譜比法分析全台灣強震站的場址。國立中央大學 地球物理研究所碩士論文。
何英銳(2003)。地磁複變解調與地震活動之相關性：以台灣崙坪全磁力資料為例。國立臺灣師範大學 地球科學研究所碩士論文。
李錫堤、邱宏智、林柏伸、鄭錦桐(2006)。強震資料庫建置與維護及使用者平台開發(一)。國立中央大學財團法人中興工程顧問社 技術報告。
梁文宗(2016)。氣象局新一代地球物理管理資料系統規劃。交通部中央氣象局地震類委託研究計畫期末報告。
陳光榮(1993)。台灣地區Qp值之空間分佈及其特性。國立台灣大學 海洋研究所博士論文。
陳光榮(2010)。台灣北部地區中大型地震震波時頻即時分析系統。交通部中央氣象局101年度地震類委託研究計劃。
張議仁(2004)。 大屯火山區三維Vp、Vp/Vs及Qp構造反演。國立中正大學 地震研究所碩士論文。
劉貞伶(2009)。以複變解調法分析竹子湖強震站最大地動加速度的異常放大現象。國立臺灣師範大學 地球科學系碩士論文。