A Simultaneous Dynamic Optimization Approach For Parameter Estimation In Water Quality Models

VANINA G. ESTRADA; ELISA PARODI; M. SOLEDAD DIAZ

Salt Lake City, EEUU

Congreso; AIChE Annual Meeting 2007; 2007

AIChE

In this work, we formulate a parameter estimation problem with a least-squares objective function subject to a partial differential algebraic equations (PDE) model resulting from temporal and spatial (one dimensional, along water column) dynamic mass balances in phytoplankton (in the form of diatoms, green algae and cyanobacteria), dissolved oxygen and nutrients, such as nitrate, ammonium, organic nitrogen, silica, phosphate and organic phosphorus. Algebraic equations represent profiles for temperature, solar radiation and river inflows, in addition to the calculation of most factors that affect rate equations, such as effect of solar radiation, nutrients, etc. The PDE is transformed into an ordinary differential equation system by spatially discretizing the PDE into sets of ordinary differential-algebraic equations (DAE) (Rodriguez and Diaz, 2006). The DAE optimization problem is solved through a simultaneous approach by transforming it into a large-scale nonlinear programming (NLP) problem by representing state and control variables profiles by polynomial functions over finite elements in time. The discretized NLP model has more than 10,000 variables. Data sets from an entire year, with a frequency of twice a week have been included. The present eutrophication model has been formulated for Lake Paso de las Piedras, an artificial lake which supplies drinking water for more than 400,000 inhabitants in Bahia Blanca (Argentina) and for industrial purposes at a petrochemical complex nearby. The high content of phosphorus and nitrogen in Paso de las Piedras Lake is consequence of agricultural activities. The discretized NLP problem has been solved with a reduced successive quadratic programming algorithm (Biegler et al., 2002). Numerical results show good agreement with values from the literature. The model is currently being validated with recently obtained additional data from the lake.