Finite difference time domain (FDTD) method is used for the simulation and analysis of electromagnetic field in the top coupling layer of GaAs/AlGaAs quantum well infrared photodetector (QWIP). Simulation results demonstrated the coupling efficiencies and distributions of electromagnetic (EM) field in a variety of 2D photonic crystal coupling layer structures. A photonic crystal structure for bi-color-QWIP is demonstrated with high coupling efficiency for two wavelengths.