Three-dimension geological faults modeling rely on multi-sources faults data heavily. A series of approaches are adopted to triangulate geological faults by analyzing different data types. For multi faults, the priority of each fault is set according to its intersecting or restricting relationships. Based on the highest priority of faults within research region, a new Recursion algorithm dividing the geological bodies with multi faults is presented. By using the dividing algorithm, the process how to construct constraining surfaces and the method how to close constraining surfaces into entities are also been elaborated.