Our study proposes a three-dimensional analytical solution for the self-potential (SP) method to detect seafloor massive sulfide deposits. Using the mirror image method, we derived formulas for the SP generated by polarized bodies in layered media. We conducted experiments with a simulated  environment and compared the results with our theoretical predictions. The findings confirmed a high degree of accuracy, demonstrating the reliability of our approach for faster and precise SP modeling.