Formation damage in petroleum reservoirs is a combination of many complicated phenomena. The formation damage reduces the productivity index significantly. The damage mechanisms may occur from the stage of drilling process to tertiary oil recovery stage. This study investigates the formation damage of petroleum wells due to shooting in well completion. The study is conducted with numerical and experimental investigations. The experimental set is a prototype. Two techniques are considered for numerical simulations: one is ANSYS-CFX and other one is a computational methodology which is based on weighted residual collocation method. Single phase and two-phase flows in porous media have been taken into consideration to explore the productivity index, which is affected by formation damage. Therefore, a better understanding of damage mechanisms for various reservoir conditions can lead to optimize oil recovery rate.