A physical model has been established previously to simulate the cathodic protection (CP) system of pipeline passing through the thermocline environment in deep water and the cathodic potential distribution was numerically calculated by finite element method (FEM), which has been verified by laboratory results. In order to extend the application of the numerical model to CP design and optimization, new various situations of anode placement, such as different number and positions, were simulated by the same model. The potential distributions on the pipeline surface and also in the whole domain of the physical model were predicted to evaluate the sacrificial anode throwing power in the heterogeneous environment. The optimum anode positions were determined for different anode number. The results showed that the anode number and also the positions had a significant influence on the potential distributions.