Abstract: In this paper we discuss the boundary element method for solving the free boundary problems of non-uniform media. A boundary integral equation with integral kernals of weak singularity is presented with its application to the seepage tree boundary problem in which the eartbdams are piece-wise homogeneous.