In this article, we have numerically investigated an intense terahertz (THz) pulses generation in gaseous plasma based on the third-order nonlinear effect, four-wave mixing rectification (FWMR). We have proposed that the fundamental fields and second-harmonic field of ultra-short pulse lasers are combined and focused into a very small gas chamber to induce a gaseous plasma, which intense THz pulse is produced. To understand the THz generation process, the first-order multiple-scale perturbation method (MSPM) has been utilized to derive a set of nonlinear coupled-mode equations for interacting fields such as two fundamental fields, a second-harmonic field, and a THz field. Then, we have simulate the intense THz-pulse generation by using split step-beam propagation method (SS-BPM) and calculated output THz intensities. Finally, the output THz intensities generated from induced air, nitrogen, and argon plasma have been compared.