Integral equation method for 3D modeling of electromagnetic fields in complex structures with inhomogeneous background conductivity