Numerical modeling of a North American hydraulic fracture experiment is done to demonstrate the accuracy with which the volume containing proppant could be estimated when electrically conductive proppant is used. An electromagnetic (EM) acquisition system with surface electric and magnetic field receivers and a grounded electric dipole source is simulated. The source has one electrode on the surface and one down a steel-cased lateral well that is adjacent to the lateral well that is being hydraulically fractured. The simulations are performed using measured EM noise at the site during hydraulic fracturing. A 3D OcTree finite-volume code is used that allows very fine meshing around the wells and fractures that expands rapidly toward the boundaries keeping memory requirements within available resources. The effect of the steel casings is modeled in the forward and inverse solutions. Possible scenarios for source-receiver configurations, proppant conductivity, number of perforations per frac stage, variations in the steel casing properties, as well as geometric errors in the locations of receivers and in the placement of lateral wells are considered. Hydraulic fracture stages are modeled as 3D geobodies with variability in the direction perpendicular to the well. Frac stages are embedded in a layered background model built from logged resistivities. The inversion of the EM data starts with the pre-frac data to recover the anisotropic layered background conductivity, steel casing conductivity, and susceptibility. Data differencing between the frac stage and the background or between successive frac stages is used for inversion of frac stage properties. A parametric box model is fit to each stage to estimate the volume, length, height, and mean stage conductivity. Hundreds of inversions starting from random parameter values are run to calculate parameter means and standard deviations. The mean values of the recovered volume, length, and height are all within 20% of the true values.