Chen＇s technique of computing synthetic seismograms, which decomposes every vector with a set of basis of orthogonality and completeness before applying the Luco-Apsel-Chen （LAC） generalized reflection and transmission coefficients method, is confirmed to be efficient in dealing with elastic waves in multi-layered media and accurate in any frequency range. In this article, we extend Chen＇s technique to the computation of coupled seismic and electromagnetic （EM） waves in layered porous media. Expanding the involved mechanical and electromagnetic fields by a set of scalar and vector wave-function basis, we obtain the fundamental equations which are subsequently solved by using a recently developed version of the LAC generalized reflection and transmission coefficients method. Our approach and corresponding program is validated by reciprocity tests. We also show a numerical example of a two-layer model with an explosion source. The P-to-EM conversion waves radiated from the interface may have potential application.
