Summary: The truncated pluri-Gaussian model is a powerful tool for representing realistic spatial distributions of facies in reservoir characterization. It is suitable for generating stochastic three-dimensional facies realizations with complex vertical and lateral relationships such as are observed in algal mound behavior. Truncated pluri-Gaussian realizations account for anisotropies and relative proportions of the facies. Despite their advantages, truncated pluri-Gaussian models have not been extensively used in data assimilation techniques such as ensemble-Kalman-based algorithms. One of the major limitations encountered in the existing implementations is the difficulty of preserving facies observations at well locations through the data assimilation procedure arising even for weakly correlated data. In this work, the problem of maintaining consistency of realizations with facies is solved by merging the data assimilation algorithm (Levenberg-Marquardt ensemble randomized maximum likelihood) with an interior-point method suitable of inequality constraints. The iterative ensemble smoother is effective at assimilating highly nonlinear production data, while the interior-point method takes into account the inequality constraints on the Gaussian model variables during the data assimilation. The formulation uses an objective function that includes a model-mismatch term, a data-mismatch term and a boundary penalization function. The method allows for its approximate version neglecting model-mismatch. The algorithm was tested on a three-dimensional synthetic reservoir mimicking the algal mounds shapes of an outcrop in the Paradox Basin of Utah and includes a large number of strongly correlated facies data and significant change of petrophysical properties for different facies. It resulted in good decrease in the data-mismatch while preserving the mound structure realism and variability in the ensemble.