Unstructured grids are automatically constructed in a domain containing complex internal boundaries. Simulation grids are constructed for reservoirs or fields which contain complex fault planes. Reconciling among generated fault grid-points and other reservoir/field grid-points is performed, enabling the use of unconstrained Delaunay triangulation. High-quality orthogonal unstructured grids are provided with good convergence properties for reservoir simulation.