This paper proposes a two-stage analytical methodology for predicting the response of pile groups in layered fluid-saturated soils induced by tunnel excavation. The Loganathan analytical approach is adopted for estimating the free-field response induced by tunneling. At the second phase, the solution to the surrounding soil of pile groups is employed as the kernel function based on the analytical layer-element method (ALEM). Subsequently, the analysis of soil-structure interaction is presented using the boundary element method (BEM). The Timoshenko-beam model is utilized to simulate each monopile combined with finite element method (FEM). Finally, the responses of adjacent pile groups are derived by the coordination conditions and coupling of BEM and FEM. Comparisons are made between the presented analytical solutions and existing experimental and numerical results. Parametric analyses are performed to evaluate the influence of soil permeability, material stratification, tunnel-pile distance, tunnel buried depth, and pile slenderness ratio on the mechanical behavior of grouped piles. The results show that the pile group in hard top and soft bottom saturated soils demonstrates better load-bearing capacity than that in soft upper and hard lower soils. Then, with the increment of tunnel-pile distance, the displacement and bending moment are reduced by 22.2 %-50 % and 41.2 %-66.7 %, respectively, and the increased rate of mechanical response over time slows down. Moreover, increasing the burial depth reduces the displacement by 8.6 %-33.3 % and the bending moment by 33.4 %-50 %, respectively, and the locations of peak responses increase from 0.6 times the pile length to 0.8 times the pile length.