We investigate the leading power corrections to the decay rates and distributions in the decay B-->X(s)e(+)e(-) in the standard model (SM) using heavy quark expansion (HQE) in (1/m(b)) and a phenomenological model implementing the Fermi motion effects of the b quark bound in the B hadron. In the HQE method, we find that including the leading power corrections the decay width Gamma(B-->X(s)e(+)e(-)) decreases by about 4% and the branching ratio B(B-->X(s)e(+)e) by about 1.5% from their (respective) parton model values. The dilepton invariant mass spectrum is found to be stable against power corrections over a good part of this spectrum. However, near the high-mass end point this distribution becomes negative with the current value of the nonperturbative parameter lambda(2) (the lambda(1)-dependent corrections are found to be innocuous), implying the breakdown of the HQE method in this region. Our results are at variance with the existing ones in the literature in both the decay rate and the invariant dilepton mass distribution calculated in the HQE approach. As an alternative, we implement the nonperturbative effects in the decay B-->X(s)e(+)e(-) using a phenomenologically motivated Gaussian Fermi motion model. We find small corrections to the branching ratio, but the nonperturbative effects are perceptible in both the dilepton mass distribution and the forward-backward asymmetry in the high dilepton mass region. Using this model for estimating the nonperturbative effects, modeling the dominant long distance contributions from the decays B-->X-s+(J/psi,psi',...)-->X(s)e(+)e(-), and taking into account the next-to-leading order perturbative QCD corrections in b-->se(+)e(-), we present the decay rates and distributions for the inclusive process B-->X(s)e(+)e(-) in the SM.