In this paper, we propose a NMPC scheme based on NARX-Laguerre model in its discrete-time structure in order to control the pH neutralization process. This model is the result of the development of discrete-time NARX model parameters on five independent Laguerre bases introducing the benefit of an essentially diminished number of parameters compared to the classical NARX model. The decrease in parametric complexity actually depends on choosing the ideal poles that characterize the Laguerre bases. In this paper, we develop an analytical technique to optimize the NARX-Laguerre poles. The parameters of the NARX-Laguerre model are reached using a recursive technique. The proposed model is utilized to integrate a nonlinear model predictive control, where we create a j-step ahead predictor on the prediction horizon [k+1,k+Np] and formulate the optimization problem using a performance criterion taking into consideration the restrictions imposed on the input and output of the process. The proposed nonlinear model predictive control strategy is validated on pH neutralization process.