Identification of unknown electric charges or sources distributed in space is made from the data observed over the field boundary using dual reciprocity boundary element models. The inhomogeneous term of the Poisson field can equivalently be expressed as the linear combination of the functions associated with the particular solutions to transform into Laplace equation. For the solution procedure, the variational formulation is employed, in which the regular boundary integral approach is incorporated to avoid the singularity. Numerical examples are presented to demonstrate the availability and the capability.