Abstract: In 3D rock-mass structure modeling,when reconstructing non-continuous stratum surface was cut by complex fault system,complexity and difficulty were unavoidable in application of the partition method.In order to solve this problem,integrative reconstruction of non-continuous stratum surface was proposed.In this method,firstly,a primary stratum surface was formed without considering faults through continuous interpolation,and surfaces of faults were also built.Secondly,the intersecting lines between all faults and the primary stratum surface were found,and regarding these intersecting lines as control,the initial stratum surface was rebuilt through re-mesh and re-interpolation.Then,by using the iterative arithmetic of intersection,the true intersecting lines between faults and stratum surface were found.Finally,after determining the affected region of faults,the coordinates of nodes on affected region of the initial stratum surface were adjusted using discrete smooth interpolation(DSI),and the final stratum surface was reconstructed.