Surface to borehole electromagnetic (SBEM) is a potential tool for reservoir monitoring. However, it is difficult to give a good resolution in reservoir imaging in the horizontal direction due to the fact that receivers are located only inside the borehole. Multiple sources are a natural choice to overcome this problem, but the resultant required memory and CPU time are large due to multiple sources and complex inhomogeneous media. A two-dimensional biconjugate gradient fast Fourier transform (BCGS-FFT) method is developed for multiple sources in a layered medium for vertical electromagnetic profiling of the surface to borehole method. With the help of the presented solver, vertical electromagnetic profiling for reservoir monitoring would be possible.