Next Article in Journal
Compatibility Analysis Between RedCap Non-Public Networks and 5G NR in TDD FR1 and FR2 Bands
Previous Article in Journal
Improving the Thermal Efficiency of Gasket Plate Heat Exchangers Used in Vegetable Oil Processing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Control of Irrigation Systems Using Predictive Methods and Disturbance Rejection

by
Jose Carreño-Zagarra
1,*,
Diana Poveda-Rodriguez
2 and
Marco Flórez
1
1
Programa de Ingeniería Electrónica, Universidad de Investigación y Desarrollo, Bucaramanga 680001, Colombia
2
Escuela de Ingenierías Eléctrica, Electrónica y de Telecomunicaciones, Universidad Industrial de Santander, Bucaramanga 680002, Colombia
*
Author to whom correspondence should be addressed.
Inventions 2025, 10(1), 11; https://doi.org/10.3390/inventions10010011
Submission received: 30 October 2024 / Revised: 1 December 2024 / Accepted: 6 December 2024 / Published: 31 January 2025

Abstract

:
Ensuring that the world’s population meets its food needs, despite water restrictions, can be significantly improved by increasing irrigation efficiency and productivity. Achieving this goal necessitates technological advancements in control systems. Therefore, the implementation of effective control systems across the entire irrigation water distribution chain is crucial and requires technological modernization. This paper presents a control scheme that combines the benefits of model predictive control (MPC) and active disturbance rejection by using generalized proportional integral (GPI) observers. The proposed control scheme was applied to a three-canal irrigation system. The simulation results confirm that the proposed controller is robust to disturbances and ensures accurate tracking for all reference levels. The controller’s performance is highlighted by the improvement in response time and considerable reduction in overshoot compared with the optimized proportional integral (PI) controllers. Additionally, the use of GPI observers allows for the precise estimation of nonlinear disturbances and phase variables, enhancing the robustness of the system. The efficiency of the observer is due to its ability to adequately estimate global additive disturbances, including unknown parameters and external disturbances in the input–output dynamics.

1. Introduction

Irrigation consumes large amounts of water. It is estimated that irrigation represents 80% of the world’s total freshwater consumption, and that only 30% of this water is delivered to plants [1]. This low efficiency in water distribution can be attributed to several factors. Significant losses occur owing to evaporation, infiltration, and surface runoff during water transport through the canals. Additionally, improper maintenance of irrigation infrastructure, such as cracked or unlined canals, exacerbates seepage losses [2]. Inefficiencies are also linked to outdated or poorly managed water allocation practices that fail to match crop water needs with the delivery schedules. Furthermore, over irrigation, often caused by a lack of precise control systems, leads to water wastage and uneven distribution across fields. These factors collectively contribute to the substantial gap between the volume of water extracted and the amount effectively utilized by plants.
Current irrigation systems face several challenges, including water scarcity, inefficient water distribution, and the need for sustainable agricultural practices. These challenges are exacerbated by climate change and increasing global food demands. The necessity for improved control of irrigation systems is critical for enhancing water use efficiency and ensuring reliable water delivery to crops. Despite advancements in control technologies, there remain gaps in the literature regarding the integration of robust control strategies that can handle the nonlinearities and uncertainties inherent in irrigation systems [3,4,5].
Automatic control techniques are widely used in irrigation canals, most of which are based on local gate control using classical approaches such as PI controllers [6]. These decentralized approaches provide reasonable performances in many cases. However, because the coupling effect between different local controllers is not considered, there is often significant loss in the performance of the control system [7]. Studies such as [8,9,10,11,12,13] have used new techniques for the optimal tuning of the parameters of these controllers.
A global controller, or centralized control, is an excellent alternative for addressing the problem of dynamic coupling in multivariate systems. Model-based predictive control (MPC) approaches are a family of control algorithms whose common properties are state and output predictions using an internal model and performing optimization with current and forecasted future data [14,15]. These controllers have shown excellent results in open channels [16,17,18,19,20,21]. However, the inclusion of disturbance models is a typical prerequisite in any standard industrial MPC implementation because the sources of disturbances are known [16]. The difficulty in tuning these disturbance models lies in the unknown nature of the disturbances.
Some studies, such as [18], add a feedforward control component to a predictive control scheme for irrigation systems. This approach aims to counteract the effects of known disturbances by adjusting the control actions in advance, thereby enhancing system performance. However, the use of such a stage requires accurate measurement of the external disturbances, which can be quite challenging in real-world applications owing to the dynamic and unpredictable nature of environmental factors, such as weather conditions, soil moisture variability, and water flow inconsistencies.
Disturbance observers are often employed to adress this issue because they provide an effective estimate of disturbances without the need for direct measurement [22,23]. Disturbance observers monitor the system output and estimate the impact of disturbances based on deviations from the expected behavior. This estimation is then used to compensate for disturbances in real-time, thereby improving the control accuracy. Moreover, the use of disturbance observers has an added benefit: not only can they handle external disturbances, but they can also account for uncertainties in the plant model parameters or unmodeled dynamics. By treating these uncertainties as part of the disturbance, the control scheme can significantly improve the robustness and adaptability of the system to varying operating conditions, ensuring a more reliable and efficient irrigation process.
This paper proposes an active disturbance rejection approach combined with model-based predictive control to improve the performance of the irrigation channels. The proposed control scheme combines the benefits of MPC and active disturbance rejection using generalized proportional integral (GPI) observers and is applied to a three-canal irrigation system. The simulation results confirm that the proposed controller is robust to disturbances and ensures the accurate tracking of all reference levels. The controller’s performance is highlighted by the improvement in response time and considerable reduction in overshoot compared with traditional proportional integral (PI) controllers. Additionally, the use of GPI observers allows for the precise estimation of nonlinear disturbances and phase variables, enhancing the robustness of the system. The contributions of this study include the development of a robust control scheme, application to a practical system, and demonstration of improved performance metrics.
The document is organized as follows: Section 2 presents the Saint-Venant equations and dynamic model of the experimental channel. Section 3 discusses the design of the predictive controller and the disturbance observer for the experimental channel. Section 4 presents a comparison of the results of the proposed robust controller with those of a classical controller. Finally, Section 5 presents conclusions based on the experimental results.

2. Modeling of Hydraulic Channels

2.1. The Saint-Venant Equations

In addition to issues related to inefficient water management and unsatisfactory service to farmers, the hydraulic limitations of irrigation canals also significantly contribute to water loss and low efficiency. Factors such as canal design, wall roughness, inadequate slopes, and flow conditions directly impact infiltration and evaporation losses, as well as conduit discharge capacity [24]. These hydraulic characteristics, which are often interdependent with control and management strategies, must be considered when seeking solutions to improve the performance of irrigation systems. Therefore, the integration of models that include hydraulic aspects is essential to achieve more robust and effective control.
The one-dimensional equations governing unsteady open channel flow are widely known as the Saint-Venant equations. These equations were derived based on several key assumptions [24,25]:
  • The flow is considered one-dimensional, implying that the velocity is uniform across any cross-section, and that the free-surface profile in the transverse direction is horizontal.
  • The curvature of the streamlines was assumed to be minimal, and the vertical acceleration of the fluid was negligible. Consequently, the pressure distribution throughout the flow is hydrostatic.
  • The flow resistance and turbulent losses are consistent with those of a steady uniform flow, regardless of depth variations, as long as the flow depth and velocity remain the same.
  • The slope of the channel bed is small enough to allow the following approximations:
    cos θ 1 , sin θ tan θ .
  • The density of the water is constant throughout the flow.
The Saint-Venant equations have traditionally been considered a complete tool for modeling irrigation system sections [24]. These processes are modeled by employing hyperbolic equations for the conservation of mass and momentum that represent the process of one-dimensional flow in a free sheet:
A w ( x , t ) t + Q x = q l
Q t + ( Q 2 / A w ) x + g A w Z x = g A w S f + k q l V ( x , t )
where A w ( x , t ) is the wetted area [ m 2 ] ; S f is the friction slope; q l ( x ) is the lateral discharge; Q ( x , t ) is the discharge [ m 3 / s ] across section A w ( x , t ) ; Z ( x , t ) is the water depth [ m ] ; V ( x , t ) is the average velocity [ m / s ] ; I is the bed slope; and g is the gravitational acceleration [ m / s 2 ] [26]. The initial and boundary conditions of the system are as follows.
Z ( x , 0 ) = Z 0 ( x ) , Q ( x , 0 ) = Q 0 ( x )
Z ( 0 , t ) = Z i ( t ) , Q ( 0 , t ) = Q i ( t ) , Z ( X f , t ) = Z f ( t ) , Q ( X f , t ) = Q f ( t )

2.2. Simplified Linear Model for a Stretch of Channel

In some studies a method was proposed that linearizing the Saint-Venant equations around a state of equilibrium transforms them into the Laplace domain [1,27,28]. For a channel section in which there are no lateral retreats, the equations that relate the y levels at the origin and the end of the section with the q flows are
y ( 0 , s ) y ( 1 , s ) = p 11 ( s ) p 12 ( s ) p 21 ( s ) p 22 ( s ) q ( 0 , s ) q ( 1 , s )
where the transfer functions in the Laplace domain p i j ( s ) are the integrator, delay, and zero (IDZ) models. The parameters that define the IDZ models are obtained analytically from the physical parameters of the section, such as the Manning coefficient, geometry, and slope, and they capture most of the flow dynamics [29].
The IDZ models corresponding to the end of tranche level y ( L , s ) are
p 21 ( s ) = 1 A d s + b d e τ d s
p 22 ( s ) = 1 A d s + b ¯ d
where A d is the equivalent area, b d and b ¯ d are static gains, and τ d is the delay. These values are functions of the physical parameters of the section; their exact expressions can be found in [29]. In distribution channels, it is recommended to use gate openings or closures as control variables rather than for manipulating flow rates. Therefore, Equation (5) must be transformed into another equivalent system of equations in which the inputs are the variations in the gates at the origin c 1 and the end of section c 2 . The flow that circulates through a semi-submerged gate is defined as
Q = C d L c C 2 Y 1 Y 2
L c is the gate width, C is the opening, C d is the discharge coefficient, Y 1 and Y 2 are the upstream and downstream levels of the gate, respectively, and g is the gravitational acceleration. Linearizing this expression to consider small deviations in the flow around an equilibrium point, the following expression is obtained,
q = K 1 Y 1 + K 2 Y 2 + K c C
where K c = Q / L c , K 1 = Q / 2 ( Y 1 Y 2 ) and K 2 = K 1 . By substituting into Equation (5) the flow in Equation (9), we introduce the general expressions of the IDZ models and obtain the equation for the end-of-span level. In this context, the “end-of-span level” refers to the water level at the downstream boundary of a canal section or span, which is crucial for maintaining the desired hydraulic conditions and ensuring water delivery accuracy.

2.3. PAC-UPC Experimental Channel

The control algorithm testing laboratory of the Polytechnic University of Catalonia (PAC-UPC for its acronym in Spanish) is a laboratory-scale model of an irrigation canal, situated at the Polytechnic University of Catalonia. Specifically designed for conducting both fundamental and applied research in the domain of irrigation canal control, this channel integrates various subfields, such as instrumentation, modeling, and water measurement of irrigation canals.
Originally conceived to replicate the significant transport delays encountered in real irrigation control scenarios, the design features a long, zero-slope, rectangular cross-sectional canal. However, owing to spatial limitations within the laboratory, a serpentine configuration was adopted to maximize the length of the canal within the confined area.
This study utilizes the PAC-UPC channel, which measures 220 m in length and 0.44 m in width, maintains a zero slope, and exhibits a transport delay of 100 s, as illustrated in Figure 1. This delay represents the time taken for water to traverse from the start to the end of the canal. The PAC-UPC includes three motorized gates and several extraction points, and is arranged into three pools, each demarcated by one of the gates, with gravity intakes positioned at the downstream end of each pool [16].
The equations that relate the final levels of the section with their corresponding origin and end gates for the three study pools of the PAC-UPC are
h 1 ( s ) = G 11 ( s ) u 1 ( s ) + G 12 ( s ) u 2 ( s ) h 2 ( s ) = G 22 ( s ) u 2 ( s ) + G 23 ( s ) u 3 ( s ) h 3 ( s ) = G 33 ( s ) u 3 ( s )
By employing parametric identification with step-type inputs, like in ref. [28], the following dynamic model for the three pools of the UPC-PAC channel were obtained.
G 11 ( s ) = 0.5 ( s + 0.007 ) s e 24 s G 12 ( s ) = 0.5 ( s + 0.011 ) s e 4 s
G 22 ( s ) = 0.6 ( s + 0.02 ) s e 45 s G 23 ( s ) = 0.4 ( s + 0.004 ) s e 4 s
G 33 ( s ) = 0.45 375 s + 1 e 25 s
As observed, all models conform to the integrator-delay-zero (IDZ) type, with the exception of the G 33 ( s ) model, which is a first-order model featuring a time delay. This variation in modeling approaches stems from the presence of a spillway at the conclusion of the third segment of the canal, rather than a semi-submerged gate.

3. Robust Controller Design

The development of a model that integrates robust strategies for irrigation systems is a complex yet feasible task that uses advanced control and modeling approaches. Tools such as model-based predictive control (MPC) allow for the incorporation of uncertainties into a model through disturbance observers and active rejection techniques. Furthermore, the reliability of these models can be assessed both quantitatively and qualitatively. Quantitative methods include performance metrics, such as the mean squared error (MSE), settling time, and overshoot reduction. On the other hand, qualitative assessments can be conducted through comparative studies with traditional controllers, sensitivity analysis to critical parameter changes, and experimental validation under dynamic conditions.
This section presents a control approach that combines model-based predictive control and a disturbance observer for active rejection of disturbances in an experimental irrigation channel.

3.1. Predictive Control Overview

Model-based predictive control (MPC) is a sophisticated control technique that is widely recognized and adopted in various industrial and academic settings. The appeal of this method lies in its ability to handle multivariable systems with constraints, thereby providing a future-oriented control strategy that adjusts actions based on predicted outcomes.
MPC operates by forecasting future outputs over a prediction horizon N, with these outputs anticipated at each time step t using a process model. The predicted outputs y ( t + k t ) for k = 1 N are dependent on the known values at time t and projected future control signals u ( t + k t ) for k = 1 N 1 . These future control signals are computed and implemented in the system.
The calculation of the control sequence u ( t + k t ) for k = 1 N 1 involves an optimization process aimed at minimizing the deviation of the process from a predefined reference path w ( t + k ) . Optimization typically minimizes a quadratic cost function based on the discrepancy between the predicted output and the reference trajectory.
The immediate control signal executed during the process is u ( t t ) . Signals from other times such as y ( t ) are disregarded because by the next time step, y ( t + 1 ) , the outcome is already known, and thus, the control strategy updates using this new information. Accordingly, the control signal for the subsequent time step, u ( t + 1 t + 1 ) , is calculated using freshly updated data.
To operationalize this strategy, MPC uses a foundational structure as depicted in Figure 2. This structure utilizes a predictive model to forecast the future behavior of the plant based on past and current values, along with optimal predictions derived from control actions. The optimizer, considering the defined cost function and any imposed constraints, computes these control actions.
One type of predictive control widely accepted in the industry and academia is generalized predictive control (GPC) [30]. The formulation of GPC using the transfer function structure is as follows:
A ( z 1 ) y ( t ) = B ( z 1 ) z d u ( t 1 ) + C ( z 1 ) e ( t ) Δ
Here, Δ = 1 z 1 , e ( t ) is zero-mean white noise, A, B, and C are polynomials in the delay operator z 1 , and d is the dead time of the system and u ( t ) and y ( t ) are the input and output variables, respectively. Polynomial C can be used for the optimal rejection of disturbances, although its usefulness is greater in improving the robustness.
The GPC algorithm involves applying an optimization problem that minimizes the following cost function:
J = j = N 1 N 2 δ ( j ) [ y ^ ( t + j t ) w ( t + j ) ] 2 + j = 1 N u λ ( j ) [ Δ u ( t + j 1 ) ] 2
where δ ( j ) and λ ( i ) are weighting sequences that are generally constant. y ^ ( t + j t ) is the optimal prediction of the output j steps forward, obtained with known data at time t. N 1 , N 2 , and N u are the minimum prediction horizon, maximum prediction horizon, and control horizon, respectively, and finally w ( t + j ) is the future reference trajectory. The optimal prediction y ^ ( t + j t ) is obtained by means of a recursive algorithm derived from the solution of a Diophantine equation [31]. The solution of this algorithm indicates that the future values of the output can be expressed as
y ( t + j ) = F j y ( t ) + E j B ( z 1 ) Δ u ( t + j d 1 ) + E j e ( t + j )
Polynomials E j and F j are derived from the Diophantine equation. The noise terms of Equation (13) will be investigated in the future. Therefore, the best prediction is given by
y ^ ( t + j t ) = G j ( z 1 ) Δ u ( t + j d 1 ) + F j ( z 1 ) y ( t )
where G j ( z 1 ) = E j ( z 1 ) B ( z 1 ) .
The solution to the GPC problem involves obtaining the sequence of control signals u ( t ) , u ( t + 1 ) , ⋯, u ( t + N ) that optimizes the cost function. Since the process has a dead time of d sampling periods, the system output will only be influenced by the control signal u ( t ) after d + 1 sampling periods. Therefore, the horizons of the GPC problem can be redefined as N 1 = d + 1 , N 2 = d + N , and N u = N . Considering the following set of predictions of j steps, it is noted that
y ^ ( t + d + 1 ) = G d + 1 Δ u ( t ) + F d + 1 y ( t ) y ^ ( t + d + 2 ) = G d + 2 Δ u ( t + 1 ) + F d + 2 y ( t ) y ^ ( t + d + N ) = G d + N Δ u ( t + N 1 ) + F d + N y ( t )
which can be written as
y = G u + F ( z 1 ) y ( t ) + G ( z 1 ) Δ u ( t 1 )
with
y = [ y ^ ( t + d + 1 t ) y ^ ( t + d + N t ) ] T
u = [ Δ u ( t ) Δ u ( t + 1 ) Δ u ( t + N 1 ) ] T
G = g 0 0 0 g 1 g 0 0 g N 1 g N 2 g 0
G ( z 1 ) = ( G d + 1 ( z 1 ) g 0 ) z ( G d + 2 ( z 1 ) g 0 g 1 ( z 1 ) ) z 2 ( G d + N ( z 1 g N 1 ( z ( N 1 ) ) z N )
F ( z 1 ) = F d + 1 ( z 1 ) F d + 2 ( z 1 ) F d + N ( z 1 )
The last two terms in the equation depend only on the past; therefore, these terms can be grouped in vector f of the free response. By rearranging the previous equation, the following expression is obtained,
y = Gu + f
Under null initial conditions, the free response also becomes zero. By applying a unit step to the input of the system at time t, we have Δ u ( t ) = 1 ,   Δ u ( t + 1 ) = 0 , , Δ u ( t + N 1 ) = 0 . The output sequence [ y ^ ( t + 1 ) ,   y ^ ( t + 2 ) , , y ^ ( t + N ) ] T is the same as the first column of the matrix G .
Therefore, the first column of the matrix G can be calculated as the response of the system to the unit step in the manipulated variable. The free response can be calculated recursively as follows:
f j + 1 = z ( 1 A ˜ ( z 1 ) ) f j + B ( z 1 ) Δ u ( t d + j )
where f 0 = y ( t ) and Δ u ( t + j ) = 0 .
Finally, the cost function is expressed in a matrix scheme, as follows:
J = ( Gu + f w ) T ( Gu + f w ) + λ u T u
where w = w ( t + d + 1 ) w ( t + d + N ) T .

3.2. Disturbance Observer

Consider the following n-dimensional dynamical system:
y ( n ) = K ( t , y ) u + ψ ( t )
where K ( t , y ) is known, uniformly bounded, and far from zero, and the function ψ ( t ) can be unknown and is uniformly bounded in an absolute manner, along with all of its time derivatives up to a finite order m. To follow a smooth known trajectory y * ( t ) , t [ 0 , ) by means of a feedback control law u, the following observer-based linearizing controller is proposed,
u = 1 K ( t , y ) y * ( t ) ( n ) i = 0 n 1 k i y i y * ( t ) ( i ) ψ ^ ( t )
where y i , i = 0 , 1 , , n 1 , and ψ ^ = z 1 are the variables generated by the following extended Luenberger observer, also known as the GPI observer:
y ˙ 0 = y 1 + λ m + n 1 ( y y 0 ) y ˙ j = y j + 1 + λ m + n j 1 ( y y 0 ) , j = 1 , , n 2 y ˙ n 1 = K ( t , y ) u + z 1 + λ m ( y y 0 ) z ˙ i = z i + 1 + λ m 1 ( y y 0 ) , i = 2 , , m 1 z ˙ m = λ 0 ( y y 0 )
The controller coefficients k j are selected in such a way that the following polynomial in the complex variable s is Hurwitz:
p c l ( s ) = s n + k n 1 s n 1 + + k 1 s + k 0 = 0
Similarly, the coefficients λ 0 , , λ m + n 1 of the observer are selected such that the following polynomial in the complex variable s is also Hurwitz.
p o b s ( s ) = s n + λ n 1 s n 1 + + λ 1 s + λ 0 = 0

3.3. Integration of GPC and GPI Observer

The implementation of the generalized predictive control (GPC) integrated with the generalized proportional integral (GPI) observer follows the steps outlined below.
1.
Formulate the Transfer Function Model of the System:
  • Develop the mathematical model of the irrigation canal system.
  • Represent the system dynamics using transfer functions.
2.
Derive the Polynomials A, B, and C:
  • Identify the system parameters and use them to derive the polynomials.
  • Ensure that these polynomials accurately capture the system’s behavior.
3.
Set Up the Cost Function for the GPC Algorith:
  • Define the prediction horizon N and control horizon.
  • Construct the cost function, which is typically a quadratic function that balances tracking performance and control effort.
4.
Solve the Diophantine Equation to Obtain E j and F j :
  • Use the system model and cost function to derive the Diophantine equation.
  • Solve for polynomials E j and F j , which are used to predict future outputs.
5.
Implement the GPC Algorithm to Optimize the Control Signals:
  • The derived polynomials and cost functions were used to compute the optimal control inputs at each time step.
  • Ensure the control signals minimize the cost function while maintaining system stability and performance.
6.
Design the GPI Observer:
  • Formulate the observer model to estimate the system states and disturbances.
  • Tune the observer gains to ensure accurate and fast disturbance estimation.
7.
Use the Disturbance Estimates from the GPI Observer to Adjust the Control Inputs in the GPC Algorithm:
  • Integrate the disturbance estimates into the GPC framework.
  • The control inputs are adjusted dynamically based on the observer’s output to improve disturbance rejection.

3.4. Robust GPC with GPI Observer Design

This study aimed to regulate the water level at the downstream end of each canal pool by adjusting the discharge (both upstream and downstream) in response to both known and unknown disturbances. These disturbances include deviations such as the discharge extracted from the canal at the downstream end of the pools.
The parameters for the predictive controller were selected based on an analysis described in [16]. The weight parameter δ in Equation (13) was calculated as the reciprocals of the squares of the maximum permissible values:
δ i = 1 M V P e 2
where MVPe is the expected maximum allowed value of the water level error signal and δ i is the i-th element of the δ weight matrix. Similarly, the parameters for weight λ of the cost function were established as the reciprocals of the squares of the maximum allowable deviations of the control variable. The formula used for the calculations is as follows:
λ i = 1 M V P u 2
The final tuned values of these parameters were δ i = 1111 and λ i = 4444 for all the sections. When designing the disturbance observer, the planar outputs of the multivariate system described in Equation (10) were used to formulate the dynamic system shown in Equation (25). Utilizing these balanced outputs, a GPI observer was developed for each controlled variable, with its characteristic equation provided by
p o b s ( s ) = s 3 + λ 2 s 2 + λ 1 s + λ 0 = 0
For Section 1 and Section 2 of the canal the following characteristic polynomial was used for the observation error:
p o b s ( s ) = s + 0.8 3
where the gains of both observers were given by λ 0 = 0.512 , λ 1 = 1.92 y λ 2 = 2.4 .
Owing to its slow dynamics, in Section 3, the following characteristic polynomial for the observation error was determined.
p o b s ( s ) = s + 1 3
whereby the observer’s gains were given by λ 0 = 1 , λ 1 = 3 and λ 2 = 3 .

4. Numerical Validation Discussion

This section presents the simulation results of the GPC scheme with the proposed disturbance observer, applied to the linear model of the PAC-UPC experimental canal described in Section 2.3. The prediction horizons were set to 15 samples for segment 1 and 2 and 40 samples for segment 3. A sampling period of 10 s was used for the measurements and control, with gate adjustments performed at the same interval.
Figure 3, Figure 4 and Figure 5 present the simulation outcomes of the control tests conducted on the UPC-PAC channel. In these charts, the red line illustrates the performance of the tuned PI controller, whereas the blue line depicts the results of the proposed control scheme. The PI controllers were designed by optimizing the integral time squared error (ITSE) for each section, that is, by minimizing the following objective function:
J ( t ) = 0 t e t 2 d t
In Figure 3, a reference change is implemented for the level described in Section 1. It is observed that the GPC controller with disturbance rejection rapidly attains the set reference compared with the optimized PI controller. While the PI controller experiences an overshoot greater than 25% and takes approximately 12 min to adjust to the 10 cm reference change, the robust controller achieves the target in approximately 3 min. However, the cost of ensuring accurate tracking of the reference involves a slight fluctuation in the other two controlled variables: a maximum deviation of 2 cm for y 2 and approximately 4 mm for y 3 .
In Figure 4, a reference change of 10 cm is displayed for the level in the second segment. Once again, the proposed controller demonstrates rapid reference tracking in comparison to the tuned PI controller. The MPC control with observer stabilizes within 5 min, while the PI controller requires approximately 13.5 min and exhibits significant overshoot. Analyzing the behavior of the level in span 1, the observer-based control swiftly rejects external disturbances, whereas the PI controller shows prolonged oscillations. For the variable y 3 , a maximum fluctuation of 5 mm is observed, which is considered negligible.
In Figure 5, a reference change of 5 cm is shown for the level in the third segment. Furthemore, predictive control utilizing active disturbance rejection significantly outperforms the classic tuned controller. The PI controller requires 20.5 min to stabilize, whereas the robust control accomplishes this in only 7 min. Additionally, it should be noted that the levels in the other two sections are immediately regulated in response to disturbances arising from the reference change in Section 3.
The experimental test flume used in this study lacks a gradient and has a serpentine shape, which may influence the hydraulic conditions, including the flow patterns and hydraulic losses. These characteristics were considered in the design and analysis of the experiments. The absence of gradient simplifies the system dynamics, allowing a focus on the controller’s ability to handle disturbances and uncertainties, whereas the serpentine shape introduces realistic flow complexities that are common in irrigation systems with varying geometries. Although these features may not fully replicate all hydraulic conditions in real-world irrigation canals, they provide a controlled environment for evaluating the proposed control strategy under consistent and reproducible conditions.

5. Conclusions

This study introduces a novel control scheme that seamlessly merges the advantages of model predictive control with the precision of active disturbance rejection, utilizing GPI observers. This innovative control framework finds practical application in a three-canal irrigation system, serving as a testament to its real-world viability.
The simulation results, derived from a comprehensive analysis, unequivocally validate the robustness of the proposed controller when faced with disturbances, while also affirming its exceptional ability to accurately track the desired setpoints across various levels. This level of performance is underpinned by the efficacy of the observers, which plays a pivotal role in the control system. The key findings from this research indicate a significant improvement in response time and considerable reduction in overshoot compared with traditional PI controllers, highlighting the effectiveness of the proposed scheme.
Future work should focus on extending this control strategy to larger and more complex irrigation networks, incorporating real-time data analytics and machine learning techniques to enhance the adaptability and efficiency of the control system. Additionally, exploring the integration of IoT devices and sensors can provide more precise and timely data, further improving control accuracy and resource management.
The broader implications of this study in the field of irrigation and water resource management are substantial. Implementing such advanced control strategies can lead to more sustainable agricultural practices by optimizing water usage and reducing waste. This contributes not only to better crop yields but also to the conservation of vital water resources, which is crucial in the face of global climate change and increasing water scarcity.
By addressing both the technical and practical aspects of irrigation control, this research paves the way for future innovations in smart irrigation systems, ultimately contributing to the advancement of precision agriculture and efficient water management practices.

Author Contributions

Conceptualization, J.C.-Z. and D.P.-R.; methodology, J.C.-Z. and M.F.; software, J.C.-Z.; validation, J.C.-Z., D.P.-R. and M.F.; formal analysis, J.C.-Z.; investigation, J.C.-Z.; resources, J.C.-Z.; data curation, M.F.; writing—original draft preparation, J.C.-Z.; writing—review and editing, J.C.-Z.; visualization, J.C.-Z.; supervision, J.C.-Z.; project administration, J.C.-Z.; funding acquisition, J.C.-Z. and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Universidad de Investigación y Desarrollo (UDI).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
MPCModel Predictive Control
ADRCActive Disturbance Rejection Control
GPIGeneralized Proportional Integral
GPCGeneralized Predictive Control

References

  1. Litrico, X.; Fromion, V. Modeling and Control of Hydrosystems; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  2. Plusquellec, H.; Burt, C.; Wolter, H.W. Modern Water Control in Irrigation: Concepts, Issues, and Applications; World Bank: Washington, DC, USA, 1994. [Google Scholar]
  3. Prawarisa, P.; Kuijk, A.; Kaijser, H. Digital Twin for Irrigation Canal System Control: A Case Study. Comput. Electron. Agric. 2023, 192, 106556. [Google Scholar]
  4. Garcia, L.; Ocampo-Martinez, C.; Puig, V. A Review of Modern Control Techniques for Irrigation Systems. Control Eng. Pract. 2018, 74, 107–124. [Google Scholar]
  5. Torres, M.; Perez, A.; Hernandez, C. Modeling and Control of Irrigation Canals: Advances and Applications. J. Irrig. Drain. Eng. 2023, 149, 04022043. [Google Scholar]
  6. Malaterre, P.O.; Rogers, D.C.; Schuurmans, J. Classification of canal control algorithms. J. Irrig. Drain. Eng. 1998, 124, 3–10. [Google Scholar] [CrossRef]
  7. Álvarez, A.; Ridao, M.; Ramirez, D.; Sánchez, L. Constrained predictive control of an irrigation canal. J. Irrig. Drain. Eng. 2013, 139, 841–854. [Google Scholar] [CrossRef]
  8. Clemmens, A.; Schuurmans, J. Simple optimal downstream feedback canal controllers: Theory. J. Irrig. Drain. Eng. 2004, 130, 26–34. [Google Scholar] [CrossRef]
  9. van Overloop, P.J.; Schuurmans, J.; Brouwer, R.; Burt, C.M. Multiple-model optimization of proportional integral controllers on canals. J. Irrig. Drain. Eng. 2005, 131, 190–196. [Google Scholar] [CrossRef]
  10. Litrico, X.; Fromion, V.; Baume, J.P. Tuning of robust distant downstream PI controllers for an irrigation canal pool. II: Implementation issues. J. Irrig. Drain. Eng. 2006, 132, 369–379. [Google Scholar] [CrossRef]
  11. Litrico, X.; Malaterre, P.O.; Baume, J.P.; Vion, P.Y.; Ribot-Bruno, J. Automatic tuning of PI controllers for an irrigation canal pool. J. Irrig. Drain. Eng. 2007, 133, 27–37. [Google Scholar] [CrossRef]
  12. Piao, X.; Burt, C. Tuning Algorithms for Automated Canal Control; ITRC ARI 04-3-005; Cal Poly: San Luis Obispo, CA, USA, 2005. [Google Scholar]
  13. Lozano, D.; Arranja, C.; Rijo, M.; Mateos, L. Simulation of automatic control of an irrigation canal. Agric. Water Manag. 2010, 97, 91–100. [Google Scholar] [CrossRef]
  14. Maciejowski, J.M. Predictive Control: With Constraints; Pearson Education: London, UK, 2002. [Google Scholar]
  15. Rossiter, J.A. Model-Based Predictive Control: A Practical Approach; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  16. Horváth, K.; Galvis, E.; Valentín, M.G.; Rodellar, J. New offset-free method for model predictive control of open channels. Control Eng. Pract. 2015, 41, 13–25. [Google Scholar] [CrossRef]
  17. Horváth, K.; van Overloop, P.J.; Galvis, E.; Gómez, M.; Rodellar, J. Multivariable model predictive control of water levels on a laboratory canal. In Advances in Hydroinformatics; Springer: Berlin/Heidelberg, Germany, 2014; pp. 77–92. [Google Scholar]
  18. Aguilar, J.V.; Langarita, P.; Linares, L.; Rodellar, J. Automatic control of flows and levels in an irrigation canal. IEEE Trans. Ind. Appl. 2009, 45, 2198–2208. [Google Scholar] [CrossRef]
  19. Aguilar, J.V.; Langarita, P.; Linares, L.; Rodellar, J.; Soler, J. Adaptive predictive expert control of levels in large canals for irrigation water distribution. Int. J. Adapt. Control Signal Process. 2012, 26, 945–960. [Google Scholar] [CrossRef]
  20. Van Overloop, P.J. Model Predictive Control on Open Water Systems; IOS Press: Amsterdam, The Netherlands, 2006. [Google Scholar]
  21. Zafra-Cabeza, A.; Maestre, J.; Ridao, M.A.; Camacho, E.F.; Sánchez, L. Hierarchical distributed model predictive control for risk mitigation: An irrigation canal case study. In Proceedings of the 2011 American Control Conference, San Francisco, CA, USA, 29 June–1 July 2011; IEEE: Piscataway, NJ, USA, 2011; pp. 3172–3177. [Google Scholar]
  22. Han, J. From PID to active disturbance rejection control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  23. Gao, Z. Active disturbance rejection control: A paradigm shift in feedback control system design. In Proceedings of the 2006 American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; IEEE: Piscataway, NJ, USA, 2006; pp. 2399–2405. [Google Scholar]
  24. Akan, A.O. Open Channel Hydraulics; Elsevier: Amsterdam, The Netherlands, 2011. [Google Scholar]
  25. Chanson, H. Hydraulics of Open Channel Flow; Elsevier: Amsterdam, The Netherlands, 2004. [Google Scholar]
  26. Litrico, X.; Fromion, V. Infinite dimensional modelling of open-channel hydraulic systems for control purposes. In Proceedings of the 41st IEEE Conference on Decision and Control, Las Vegas, NV, USA, 10–13 December 2002; IEEE: Piscataway, NJ, USA, 2002; Volume 2, pp. 1681–1686. [Google Scholar]
  27. Litrico, X.; Fromion, V. Frequency modeling of open-channel flow. J. Hydraul. Eng. 2004, 130, 806–815. [Google Scholar] [CrossRef]
  28. Horváth, K.; Galvis, E.; Rodellar, J.; Gómez, M. Experimental comparison of canal models for control purposes using simulation and laboratory experiments. J. Hydroinform. 2014, 16, 1390–1408. [Google Scholar] [CrossRef]
  29. Litrico, X.; Fromion, V. Analytical approximation of open-channel flow for controller design. Appl. Math. Model. 2004, 28, 677–695. [Google Scholar] [CrossRef]
  30. Clarke, D.W. Application of generalized predictive control to industrial processes. IEEE Control Syst. Mag. 1988, 8, 49–55. [Google Scholar] [CrossRef]
  31. Clarke, D.W.; Mohtadi, C.; Tuffs, P. Generalized predictive control—Part I. The basic algorithm. Automatica 1987, 23, 137–148. [Google Scholar] [CrossRef]
Figure 1. PAC-UPC three-pool configuration.
Figure 1. PAC-UPC three-pool configuration.
Inventions 10 00011 g001
Figure 2. Model-based predictive control (MPC).
Figure 2. Model-based predictive control (MPC).
Inventions 10 00011 g002
Figure 3. Controller validation for reference tracking for segment 1.
Figure 3. Controller validation for reference tracking for segment 1.
Inventions 10 00011 g003
Figure 4. Controller validation for reference tracking for segment 2.
Figure 4. Controller validation for reference tracking for segment 2.
Inventions 10 00011 g004
Figure 5. Controller validation for reference tracking for segment 3.
Figure 5. Controller validation for reference tracking for segment 3.
Inventions 10 00011 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Carreño-Zagarra, J.; Poveda-Rodriguez, D.; Flórez, M. Robust Control of Irrigation Systems Using Predictive Methods and Disturbance Rejection. Inventions 2025, 10, 11. https://doi.org/10.3390/inventions10010011

AMA Style

Carreño-Zagarra J, Poveda-Rodriguez D, Flórez M. Robust Control of Irrigation Systems Using Predictive Methods and Disturbance Rejection. Inventions. 2025; 10(1):11. https://doi.org/10.3390/inventions10010011

Chicago/Turabian Style

Carreño-Zagarra, Jose, Diana Poveda-Rodriguez, and Marco Flórez. 2025. "Robust Control of Irrigation Systems Using Predictive Methods and Disturbance Rejection" Inventions 10, no. 1: 11. https://doi.org/10.3390/inventions10010011

APA Style

Carreño-Zagarra, J., Poveda-Rodriguez, D., & Flórez, M. (2025). Robust Control of Irrigation Systems Using Predictive Methods and Disturbance Rejection. Inventions, 10(1), 11. https://doi.org/10.3390/inventions10010011

Article Metrics

Back to TopTop