Abstract：The work of the first part of this paper is continued by describing the numerical simulation technique of steady state flow in three-dimensional artificial fracture network. Every fracture disk is taken as a two-dimensional flow domain with a uniform transmissivity，in which the flow is determined by the Laplace equation. The boundary element method is used to calculate the steady flow of single fracture. For flow between two intersecting fractures，additional internal boundary elements are introduced to represent the line of intersection between the fractures. For a network of intersecting fractures，with consideration of mass balance and equivalence of heads along the interboundaries，a large matrix can be formed. All the unknown head and fluxes can be determined by solving this matrix. In addition，the fundamental solutions of the boundary element method for the governing function of flow problem is discussed and the equivalence of three typical fundamental solutions is proved.
