The new meshless method based on the cover least square approximation, which utilize cover instead of

points that are used in the moving least square approximation, is proposed. For the cover distribution, multi

scale voxel data is used for the ease of analysis in 3D solid. Several 3D examples are shown for

demonstration. The method is applied to adaptive analysis of seepage flow problem of rock including

complicated cracks. The rock model with 3 cracks and 1000 cracks are analyzed