Several industrial systems are characterized by high nonlinearities with wide operating ranges and large set point changes. Identification and representation of these systems represent a challenge, especially for control engineers. Multimodel technique is one effective approach that can be used to describe nonlinear systems through the combination of several submodels, where each is contributing to the output with a certain degree of validity. One major concern in this technique, especially for systems with unknown operating conditions, is the partitioning of the system's operating space and thus the identification of different submodels. This paper proposes a three-stage approach to obtain a multimodel representation of a nonlinear system. A reinforced combinatorial particle swarm optimization and hybrid K-means are used to determine the number of submodels and their respective parameters. The proposed method automatically optimizes the number of submodels with respect to the submodel complexity. This allows operating space partition and generation of a parsimonious number of submodels without prior knowledge. The application of this approach on several examples, including a continuous stirred tank reactor, demonstrates its effectiveness.