用于不规则界面多层介质的地震波方法比较研究：不规则地形问题

通过研究SH波入射时各种地形的散射，从计算结果的精确性、稳定性和计算效率三个方面研究Aki Larner方法、Bouchon Campillo方法、以及全局广义反-透射矩阵方法（Chen's方法），比较这三种计算具有不规则界面多层介质中理论地震图的离散波数方法. 由于利用了FFT算法，Aki Larner方法速度最快，但由于方法本身受到Rayleigh假设的限制，结果精度最差. Bouchon Campillo方法非常稳定，结果较精确；但其计算效率最低，并且不能正确地处理非常陡峭的地形问题. 全局广义反-透射矩阵方法非常稳定，结果最精确，能够很好地处理非常陡峭的地形问题，并具有适中的计算效率；因此是三种方法中的最佳方法，为计算复杂介质中的地震波动问题提供了十分有效的工具.

A COMPARATIVE STUDY ON SEISMIC WAVE METHODS FOR MULTI-LAYERED MEDIA WITH IRREGULAR INTERFACES:IRREGULAR TOPOGRAPHY PROBLEM

Through studies on the scattering from various topographies due to an incident SH wave, three methods for seismogram synthesis (Aki Larner method, Bouchon Campillo method, and global generalized R/T matrices method) are evaluated in terms of the accuracy of the results, numerical stability and computational efficiency. Using the FFT algorithm, the Aki Larner method is the most efficient of the three methods, but is least accurate because of the limitation of the Rayleigh hypothesis. The Bouchon Campillo method is quite accurate and stable, but has the lowest efficiency;moreover it cannot treat very steep topography precisely. The global generalized R/T matrices method not only is very accurate, stable and with a modest computational efficiency, but also can treat the very steep topography well. Therefore, the global generalized R/T matrices method is the best method among the three, and provides a very effective tool to synthesize seismograms for complex media.