The implicit Zbus Gauss method is widely used in distribution network power flow calculation for radial and weakly meshed networks because of its efficiency and reliability. This method is appropriate for handling PQ or PQ(V) nodes in power flow calculation but it encounters divergence problems when distributed generations (DGs) modeled as PV nodes are incorporated in the system. In order to solve the problem, this paper proposes an enhanced implicit Zbus Gauss method on the basis of compensation method. The main idea of which is the reactive power of PV nodes is revised through node-impedance-matrix and deviation of voltage amplitude after each iteration, so that the real value of reactive power of PV nodes can be obtained and the convergent ones of power flow solutions can further be achieved. This method is applied to IEEE-33 bus distribution system. Compared with the homotopy-enhanced implicit Zbus Gauss method, it demonstrates that the proposed method has better convergence and faster convergent speed. In addition, the number of PV nodes exert little influence on convergence and convergent speed of this method, which manifests the proposed method is more suitable for power flow calculation with DG modeled as PV nodes.