The theory of complex variables methods and multi-polar coordinate systems are used to solve the problem of scatter and dynamic stress concentration of elastic waves around double elliptic cavities in saturated soil. The 2D steady state Biot wave equations are decoupled by introducing the potentials and reduced to three Helmholtz equations that the potentials satisfy. The solving of the Helmholtz equations by separating variable methods obtains general solutions of the potentials. The expressions of displacements

stresses and pore pressure can be obtained by determining the unknown coefficients in the potentials and utilizing the boundary conditions of soil and pore pressure. Numerical examples are provided to show the effect of wave number

distance between centers of cavities upon the dynamic stress concentration and pore pressure concentration around the cavities under incident steady elastic wave.