本研究主要應用連續表面波試驗(Continuous Surface Wave Test，簡稱CSWT)之現地震測試驗，求得現地之頻散曲線(dispersion curve)。從實驗中，我們可以輸入一固定頻率(f )之震波，將得到之震測資料經過處理，可得到波長(L )。利用Vr=f*L ，即可得到雷利波速。雷利波速和波長的關係圖，即是頻散曲線。為了方便進行反算，則將左右兩側震測資料，共同進行多項式二次或三次迴歸。
在反算的方法上，本研究採用有限差分法來進行數值運算。利用有限差分法的原理，將2-D波傳理論，在平面應變的條件下，改寫成有限差分方程式。在本研究中，所採用的格點為600x600，邊界為鉸接(hinge)，以震源位置垂直自由表面之軸線為對稱軸，並且考慮初始條件為所有應力狀態皆為零。震源以正弦波輸入。格點距離範圍為0.25∼2m，時步大小範圍為0.0001∼0.0005sec。反算的流程，則是先利用簡化法，猜測初步評估之地盤剪力波速與地盤厚度，代入有限差分法中，建立正算模式，以求得有限差分頻散曲線；再將所得之有限差分頻散曲線跟現地迴歸之頻散曲線比較，進行試誤法，直到有限差分所得之頻散曲線與現地頻散曲線吻合，則猜測之地盤剪力波速與地盤厚度，即為反算結果。
利用有限差分法對雷利波波傳所進行的反算，針對台大校園、士林百齡橋、台北民族公園SCPT比較結果，大致都相當符合。對花蓮和平電廠之跨孔法，亦有良好的結果。將各場址反算結果與簡化法加以比較，則簡化法較易有高估的現象。

The objective of this research is to get dispersion curve by using in situ test of the Continuous Surface Wave Test. In the test, we can get some data by inputting seismic waves of a specific frequency(f ) and then get wavelength(L ) by analyzing the data. By using Vr=f*L , we get the Rayleigh Wave Velocity. The graphic relationship of the Rayleigh Wave Velocity and wavelength is dispersion curve. In order to process inversion, we use 2 or 3 orders polynomial to deal with the new data which are combined with the data from left and right side seismic test, and then make a fitting curve.
As to the inversion, this research applies Finite Difference Method(FDM) to process numerical calculation. Based on FDM theory and the plane strain condition, we rewrite the 2-D wave propagation theory into finite difference equations. In this research, the mesh is 600x600 grids, the boundary condition is hinge, and symmetry axis is a vertical line which goes through the seismic source. Moreover, the initial condition of the stress state is zero. This research assumes that the seismic source is inputted by sine waves. The range of the grid size is 0.25m to 2m and the time step is from 0.0001 to 0.0005 second. The first step of the process of inversion is applying simplify method to assume the data of the shear wave velocity and thickness of layers, and input the assumed data to the finite difference equations for accomplishing forward model to get finite difference dispersion curve. Applying the Try and Error Method to compare the results of FDM dispersion curve and in situ dispersion curve. If the FDM dispersion curve matches the in situ dispersion curve, the results of this study indicate that the ground shear wave velocity and thickness of the layers are the final output of the inversion.
This research compares the inversion results and SCPT data of the National Taiwan University campus, Pai-Ling Bridge in Shi-Ling, and Ming-Tsu Park in Taipei, and the results indicate that the inversion results meet the SCPT. Furthermore, the comparison of the inversion result and the crosshole method test in Ho-Ping Power Plant in Hua-Lien is satisfying. Comparing the inversion results of different locations to the simplify method, the latter is likely to be over-estimation.

第一章 導論1
1.1 前言1
1.2 研究動機與目的2
1.3 研究內容3
第二章 文獻回顧5
2.1 震波之相關背景5
2.2 雷利波與波傳反算相關研究5
2.3 有限差分法於波傳反算之相關研究7
第三章 波傳理論與有限差分法10
3.1 無限均質等向性線彈性介質中之波傳理論10
3.2 半無限空間均質線彈性之波傳理論14
3.3 波傳理論之有限差分分析法23
3.3.1 有限差分方程式23
3.3.2內部格點25
3.3.3 外力作用格點27
3.3.4 自由表面格點28
3.3.5 數值穩定條件29
3.3.6 數值震盪之分析29
第四章 連續表面波試驗(CSWT)設備與試驗方法35
4.1 試驗設備35
4.2 試驗方法36
第五章 試驗場址介紹39
5.1 台大校園39
5.2 五股工業區39
5.3 士林百齡橋39
5.4 台北市民族公園40
5.5 花蓮和平電廠40
5.6 蘆洲抽水站40
第六章 連續表面波分析方法47
6.1 試驗資料分析方法47
6.2 以簡化法初步反算土層動態參數48
第七章 速度-應力有限差分法之正反算模式61
7.1 正算模式61
7.2 反算模式62
第八章 分析結果與討論68
8.1 台大校園68
8.1.1 連續表面波試驗結果68
8.1.2 反算分析69
8.2 五股工業區70
8.2.1 連續表面波試驗結果70
8.2.1 反算分析70
8.3 士林百齡橋71
8.3.1 連續表面波試驗結果71
8.3.2 反算分析72
8.4 台北市民族公園72
8.4.1 連續表面波試驗結果72
8.4.2 反算分析73
8.5 花蓮和平電廠73
8.5.1 連續表面波試驗結果73
8.5.2 反算分析74
8.6 蘆洲抽水站75
8.6.1 連續表面波試驗結果75
8.6.2 反算分析75
8.7 綜合討論75
第九章 結論與建議106
9.1 結論106
9.2 建議107

英文參考文獻
1、 Abbiss C. P. " Shear wave measurements of elasticity of the ground" Geotechnique, 31(1), p91-104, 1981.
2、 Addo, K.O. and Robertson, P. K. " Shear-wave velocity measurment of soil using Rayleigh waves ", Can. Geotech. J. 29, p558-568 , 1992.
3、 Addo, K.O. " Development of Computerized SASW Technique ", Ph.D. thesis, The University of Alberta., 1992.
4、 Bolt, B.A. " Nuclear explosions and earthquakes ", W.H. Freeman and Company, San Francisco, California, 309pp., 1976.
5、 Dunkin, J.W. " Computation of modal solution in layered, elastic media at high frequencies ", Bulletin of the Seismological Society of America, Vol. 55, p335-358, 1965.
6、 Ganji, V., Gucunski, N. and Maher, A. " Detection of underground obstacles by SASW method numerical aspects ", Journal of Geotechnical and Geoenviromental Engineering, p212-219, March 1997.
7、 Gazetas G. " Vibrational characteristics of soil deposits with variable velocity ", J. Num. Anal. Meth. Geomech., 6, p1-20, 1982.
8、 Gucunski, N. and Woods, R. D. " Instrumentation for SASW testing ", Recent advances in instrumentation data acquisition and testing in soil dynamics,Proceedings of sessions sponsored by the Geotechnical Engineering Division of the American Society of Civil Engineers in conjunction with the ASCE convention Orlando, Florida, p1-16, 1991.
9、 Gucunski, N., Ganji, V. and Maher, M. H. " Effects of obstacles on Rayleigh wave dispersion obtained from the SASW test ", Soil Dynamics and Earthquake Engineering 15, p223-231, 1996.
10、 Haskell, N. A. " The dispersion of surface waves on multilayered media ", B.S.S.A., Vol. 43, No. 1,p17-34, 1953.
11、 Heukolom W. and Foster C. R. " Dynamic testing of pavements ", Trans. ASCE, 127(1), p425-456, 1962.
12、 Knopoff, L. " A matrix method for elastic wave problems ", Bulletin of the Seismological Society of America, Vol. 54, No. 1, p431-438, 1964.
13、 Kramer, S. L. " Geotechnical earthquake engineering ", University of Washington, 1996.
14、 Lamb, H. " On the propagation of tremors over the surface of an elastic solid ", Philosophical Transaction of the Royal Society, London, Ser. A, Vol. 203, p1-42, 1904.
15、 Matthews, M. C., Hope, V. S. and Clayton R.I. " The use of surface waves in the determination of ground stiffness profiles ", Proc. Instn Civ. Engrs Geotech. Engng.119, p84-95, Apr. 1996.
16、 Menzies, B.K. and Matthews, M.C. " The continuous surface-wave system: A modern technique for site investigation ", Special lecture: Indian Geotechnical Conference Madras, Dec 11-14,1996.
17、 Nazarian, S. " In situ determination of elastic moduli of soil deposits and pavement systems by spectral-analysis-of-surface-waves method ", PhD thesis, The University of Texas at Austin, Austin, Tex, 1984.
18、 Nazarian, S., Stokoe, II, K. H. " Evaluation of moduli and thickness of pavement system by SASW methods ", Report No. 256-4, Center for Transportation Research, The Univ. of Texas, Austin, Tx., p.123, 1982.
19、 Nazarian, S., Stokoe, II, K. H. " Use of surface waves in pavement evaluation ", T. R. R. 1070, T. R. B., p132-144, 1986.
20、 Rayleigh, Lord, " On waves propagated along the plane surface of an elastic solid ", London Mathematical Society Proc., 17, p4-11, 1885.
21、 Richart, F. E., Hall, J. R. and Woods, R. D. " Vibrations of soil and foundations ", Englewood Cliffs, New Jersey, 1970.
22、 Thomson, W.T. " Transmisson of elastic waves through a stratified solid medium ", Journal of Applied Physics, Vol. 21, No. 2, p89-93, 1950.
23、 Timoshenko, S. P., and Goodier, J. N. " Theory of elasticity ", McGraw-Hill Book Co., New York, 506pp., 1951.
24、 Virieux, J. " SH-wave propagation in hetergeneous media: Velocity-stress finite-difference method ", Geophysics, 49, p1933-1957, 1984.
25、 Virieux, J. " P-SV wave propagation in hetergeneous media: Velocity-stress finite-difference method ", Geophysics, 51,p889-901, 1986.
中文參考文獻
26、 朱武男，"以電子震測錐試驗（SCPT）評估土層動態參數之研究"，國立台灣大學土木工程研究所，碩士論文 ，民國87年。
27、 李咸亨，"震波之量測"，地工技術雜誌，第17期，第57-69頁，民國84年。
28、 吳偉特，"土壤動力學與大地工程"，地工技術雜誌，第9期，第5-19頁，民國74年。
29、 常正之，"應用雷利波散射曲線反算土層動態參數之研究"，國立成功大學土木工程研究所，博士論文，民國 82年。
30、 鮑亦興，"波傳反算震測法之研究"，國道新建工程局研究報告，民國82年。
31、 賴政興，"介質正反算之初步研究"，國立台灣大學應用力學研究所，碩士論文，民國80年。
32、 龔俊宏，"以彈性波檢測介面與內含物之初步研究"，國立台灣大學應用力學研究所，碩士論文，民國81年。