Abstract： Using the Helmholtz resolution method, the solidmatrix, matrix pore and fracture pore coupling equations were decoupled. A series of elastic waves equations were obtained by using the wave function expansion method. Then, with the help of the transformation theory for spherical coordinates and the superposition principle, the scattering problem of double spherical inclusions were expressed in the same coordinate system. Finally, examples of calculation and analysis of different parameters were concerned to determine their impacts on first order wave amplitude with variation to frequency. The results show that different inclusions distance and inclusion radius have significant impact on P1 to P3 wave and Swave amplitude. The value of P1 to P3 wave and Swave amplitude decrease as the distance between the inclusions increases. This impact is more significant at higher frequency. Besides, the value of P1 to P3 wave and Swave amplitude increase as the radius of the other inclusion increases. The radius has more impact on P1 and P3 wave but has a rather small influence on P2 wave and S wave amplitude.