Volume-phase holographic gratings are suitable for use in greenhouse gas detection imaging spectrometers, enabling the detection instruments to achieve high spectral resolution, high signal-to-noise ratios, and high operational efficiency. However, when utilized in the infrared wavelength band with high dispersion requirements, gratings struggle to meet the demands for low polarization sensitivity due to changes in diffraction performance caused by phase delays in the incidence of light waves with distinct polarization states, and current methods for designing bulk-phase holographic gratings require a large number of calculations that complicate the balance of diffraction properties. To overcome this problem, a design method for transmissive bulk-phase holographic gratings is proposed in this study. The proposed method combines two diffraction theories (namely, Kogelnik coupled-wave theory and rigorous coupled-wave theory) and establishes a parameter optimization sequence based on the influence of design parameters on diffraction characteristics. Kogelnik coupled-wave theory is employed to establish the initial Bragg angle range, ensuring that the diffraction efficiency and phase delay of the grating thickness curve meet the requirements for incident light waves in various polarization states. Utilizing rigorous coupled-wave theory, we optimize grating settings based on criteria such as a center wavelength diffraction efficiency greater than 95%, polarization sensitivity less than 10%, maximum bandwidth, and spectral diffraction efficiency exceeding 80%. The ideal grating parameters are ultimately determined, and the manufacturing tolerances for various grating parameters are analyzed. The design results show that the grating stripe frequency is 1067 lines per millimeter, and the diffraction efficiencies of TE and TM waves are 96% and 99.89%, respectively. The diffraction efficiency of unpolarized light is more than 88% over the whole spectral range with an average efficiency of 94.49%, an effective bandwidth of 32 nm, and a polarization sensitivity of less than 7%. These characteristics meet the performance requirements for dispersive elements based on greenhouse gas detection, the spectral resolution of the detection instrument is up to 0.1 nm, and the signal-to-noise ratio and working efficiency are improved by increasing the transmittance of the instrument.