Underground hydrogen storage has garnered interest in recent years owing to the considerable demand for clean energy. Hydrogen is more diffusive than air and has relatively low density and viscosity. These unique properties induce distinct hydrodynamic phenomena during hydrogen storage. Cushion gas has been proven to be a potential remedy for attenuating the adverse impacts of underground hydrogen storage. To investigate the influence of the cushion gas, a microscopic numerical simulation was performed with Fluent software using the Eulerian multi-fluid model. Carbon dioxide, nitrogen, and methane are usually used as the preferred candidates for cushion gases in underground hydrogen storage systems. In this study, nitrogen was used as the cushion gas and was injectedalong with hydrogen into heterogeneous porous media with volume fractions ranging from 0% to 70%. A parameterization study was then performed to elucidate the influences of the injection rate and viscosity of the fluid on the fingering pattern. Two representative types offingering, viscous fingering and capillary fingering, were observed under different gas mixtures and boundary conditions. After the simulation, an image analysis was performed to capture the evolution of the fingering pattern. The specific fingering area, number of branches, and fractal dimensions are proposed as geometric indices to describe the shape of the fingering pattern. The results showed that there was a remarkable enhancement in saturation due to the injection of the cushion gas, depending on the concentration of the gas mixture. This study offers insight on the design of gas mixture injection in underground hydrogen storage and can be furtherextended to the hydrochemo-mechanical coupled numerical simulation of multiphase gas injection in porous media.