Based on the waveguide mode (WGM) method, coupling interaction of electromagnetic wave in a groove doublet configuration is studied. The formulation obtained by WGM method for a single groove [Prog. Electromagn. Res. 18, 1–17 (1998)] is extended to two grooves. By exploring the total scattered field of the configuration, coupling interaction ratios are defined to describe the interaction between grooves quantitatively. Since each groove in this groove doublet configuration is regarded as the basic unit, the effects of coupling interaction on the scattered fields of each groove can be investigated respectively. Numerical results show that an oscillatory behavior of coupling interaction is damped with increasing groove spacing. The incident and scattering angle dependence of coupling interaction is symmetrical when the two grooves are the same. For the case of two subwavelength grooves, the coupling interaction is not sensitive to the incident angle and scattering angle. Although the case of two grooves is discussed for simplicity, the formulation developed in this article can be generalized to arbitrary number of grooves. Moreover, our study offers a simple alternative to investigate and design metallic gratings, compact directional antennas, couplers, and other devices especially in low frequency regime such as THz and microwave domain.