With the development of electronic technology, time-domain electromagnetic (TDEM) transmitters are now capable of emitting various waveforms such as trapezoidal, triangular, and semi-sinusoidal waveforms, as well as their combinations. However, the application of multi-waveform combination detection methods is still immature, and modeling methods are required to analyze the characteristics of different waveforms and complete data interpretation. Based on the Crank-Nicolson finite difference time-domain (CN-FDTD) method for discretizing the electromagnetic wave equations, the current waveform can be converted into current surface density for source-included calculations, enabling high-precision modeling of the responses to different current waveforms. The responses to multiple transmission current waveforms were compared and analyzed using various models, and the results showed that selecting targeted excitation parameters based on the characteristics of the target can improve detection sensitivity. We applied the multi-waveform combination detection method to the actual exploration in the field, obtained the pure secondary field response by subtracting the primary field of the actual current waveform from the measured data, converted the responses of different transmission current waveforms into step wave responses by deconvolution method, and performed resistivity imaging. The results of modeling and data interpretation demonstrate the effectiveness of the multi-waveform combination detection method in improving geological exploration accuracy, and promoted the development of fine-scale geophysical exploration.