The immiscible displacement behavior in porous media is crucial for oil recovery and subsurface remediation, yet how wettability influences this process across different pore structures remains unclear. Using a color-gradient lattice Boltzmann model, this study investigates immiscible displacement dynamics in porous media. Heterogeneous porous structures with various degrees of permeability were reconstructed using the quarter structure generation set algorithm, and wettability effects were analyzed with the contact angle set in the range of ${30}^{\circ }$ to ${150}^{\circ }$. The numerical results showed that porous heterogeneity greatly affects the displacement efficiency via permeability-dependent flow pathway optimization. Enhanced efficiency was observed in high-permeability media through low-tortuosity channels, whereas low-permeability systems exhibited reduced efficiency due to capillary trapping in pores with lower flow capacity. Wettability alters displacement patterns via capillary forces – under hydrophilic condition, the displacing fluid preferentially enters smaller pores. Fractal dimension and Euler number were used to quantify flow heterogeneity, revealing that increased permeability reduces flow complexity and improves connectivity. Moreover, permeability heterogeneity and wettability interact to disrupt classical linear flow responses, leading to non-monotonic efficiency trends in low-permeability systems. These findings highlight the importance of pore-scale multiphase flow in heterogeneous media and offer new insights for predicting wettability effects in subsurface flows.