To solve seismic wave scattering by a large-scale three-dimensional (3-D) hill topography, a fast parallel indirect boundary element method (IBEM) is developed by proposing a new construction method for the wave field, modifying the generalized minimum residual (GMRES) algorithm and constructing an OpenMP plus MPI parallel model. The validations of accuracy and efficiency show that this method can solve 3-D seismic response of a large-scale hill topography for broadband waves, and overcome the weakness of large storage and low efficiency of the traditional IBEM. Based on this new algorithm architecture, taking the broadband scattering of plane SV waves by a large-scale Gaussian-shaped hill of thousands-meters height as an example, the influence of several important parameters is investigated, including the incident frequency, the incident angle and the height-width and length-width ratio of the hill. The numerical results illustrate that the amplification effect on the ground motion by a near-hemispherical hill is more significant than the narrow hill. For low-frequency waves, the scattering effect of the higher hill is more pronounced, and there is only a single peak near the top of the hill. However, for high-frequency waves, rapid spatial variation of displacement amplitude appears on the hill surface.

为了解决大规模三维（3-D）丘陵地形的地震波散射问题，提出了一种新的波场构造方法，并修改了广义最小残差（ GMRES）算法并构建OpenMP加MPI并行模型。精度和效率的验证表明，该方法可以解决宽带山丘大尺度地形对宽带波的3D地震响应，克服了传统IBEM的存储量大，效率低的缺点。基于这种新的算法架构，以数千米高的大规模高斯形山丘对平面SV波的宽带散射为例，研究了几个重要参数的影响，包括入射频率，入射角和山的高宽比。数值结果表明，近半球形的山丘对地震动的放大作用比狭窄的山丘更为显着。对于低频波，较高丘陵的散射效应更为明显，并且在丘陵顶部附近只有一个峰。但是，对于高频波，在山丘表面会出现位移幅度的快速空间变化。