An algorithm and implementation of a dual-percolation model based on bond-percolation theory was introduced which combines pore-network with fracture-network. Monte-Carlo numerical simulation was conducted based on statistical analysis of key parameters,i. e. fracture population densities,Fracture angles,fracture lengths,and pore connectivity possibilities. Fractal dimension was presented to quantify fractal connectivity. Numerical results were in agreement with those of Balberg and Berkowitz et al. Finally,appropriate sample scale was discussed considering both calculation precision and computing time.