# Comparison of Methods for Assessing the Assimilation Capacity of the Kazakhstani Sector of the Ili River

Authors: A. T. Nurseitova; J. K. Jamalov; A. A. Azimov; D. B. Nurseitov; E. A. Tursunov
Journal: Modelling and Simulation in Engineering
 (2021)

---

## Abstract

A mixed inverse problem for determining the biochemical oxygen demand of water (L0) and the rate of biochemical oxygen consumption (k0), which are important indicators of water quality, has been formulated and numerically solved based on real experimental data. The inverse problem is reduced to the optimization problem consisting in minimization of the deviation of the calculated values from the experimental data, which is solved numerically using the Nelder–Mead method (zero order) and the gradient method (first order). A number of examples of processing both model experimental data and field experimental data provided by hydrological stations monitoring pollutants in the Kazakhstani part of the Ili River basin are presented. A mathematical model that adequately describes the processes in the river system has been constructed.


---

## Body

## 1. Introduction

Many models describing the processes of pollution and purification of water resources are based on the work of Streeter and Phelps published in 1925 [1], where a simple model (considered classical in our time) was proposed, which for many years satisfied the needs of engineers and other practitioners in the field of monitoring of water pollution in many countries. In the course of the development of this direction, it became clear that there are many rather important factors that have a great influence on the studied processes but are not taken into account by the classical Streeter-Phelps model. A reasonably good analysis of the models proposed over 40 years, taking into account new input data, was made in [2], but, as the author himself noted, all these models are semiempirical and require a lot of input data that are not always available for measurements. At first glance, the classical Streeter-Phelps model is simple to implement, and even its more complicated versions allow an analytical solution, which is undoubtedly an advantage for engineers, i.e., for practical use. Indeed, even with the advent of computers until the end of the 20th century, the issue of economical use of computer time in calculations had a priority position in the choice of models and methods of their implementation. Therefore, there are many works that consider different factors of oxygen consumption and recovery in water (dispersion, photosynthesis, reaeration, etc.), but all of them are based either on a large amount of input data, or on assumptions that are inconvenient in practice, but this is not the main weakness of these studies. For example, in [3], to calculate the missing (inaccessible for measurement) input data, where reaeration is taken into account, the original Streeter-Phelps model is synthesized with the shallow water equation, and here, the scientists already face the problem of cumbersome calculations. As shown in the work of Gotovtsev [4], the classical Streeter-Phelps model does not guarantee physically correct solutions for any input data, i.e., a rather difficult process of data calibration leads to a separately posed problem that requires separate consideration [2]. In [4], a modification of the classical Streeter-Phelps model (a closed Streeter-Phelps system) is proposed, which does not require additional hydrochemical analyzes; moreover, it is shown that it is physically correct. To obtain an analytical solution of the Streeter-Phelps closed-loop system, certain conditions are again set on the input data, for example, the frequency of incubation times [5], which is not applicable in practice. Thus, a positive aspect of these models is the explicit form of the desired value, which greatly facilitates the work of practitioners, but the models are complicated by introduction of additional factors imposing new conditions and restrictions on the input data, which creates a number of new problems; the solution of which already requires other resources. With development of more powerful computers, the capabilities of calculators have increased significantly and works have appeared where application of the modified Streeter-Phelps models is considered solutions to some inverse problems, for example [6]. That is, the general tendency to consider some problems as inverse to direct ones, which appeared in the 20th century [7, 8], was also reflected in the implementations of models describing the processes of pollution and purification of water resources.In this paper, the sought BOD (L0) and k0 are considered solutions to an inverse problem for which additional information is given at a fixed moment of time T>0. More exactly, a closed Streeter-Phelps model is considered, for which there are no data restrictions (taken from real sources), which is achieved by using iterative methods instead of searching for an analytical solution. It is shown that the inverse problem is ill-posed [7]; therefore, two optimization methods for its solution are considered: the semiheuristic Nelder–Mead [9] and the gradient (also known as regularizing [10]). The Nelder–Mead method is an unconstrained zero-order optimization method, i.e., the method that does not require computation of either the first or the second derivatives of the objective function, which is an undoubted advantage in the implementation of this iterative algorithm, especially in the case of many variables. However, the convergence of this method was investigated only for strictly convex functions in the one-dimensional case; for multidimensional objective functions, the convergence has not been proven (although for the two-dimensional case weak convergence was established under certain restrictions on the algorithm) [11]. The use of gradient methods (of the first order) in the case of a multidimensional objective function is complicated by the calculation of the gradient, but the convergence of these iterative algorithms has been sufficiently well studied [12]. This paper presents a comparative analysis of the numerical results of realization of the considered modified Streeter-Phelps model by the Nelder–Mead method and the steepest descent method for synthesized data. More exactly, two aspects of the methods are compared: the complexity of realization and investigation of convergence. The fact that for the Nelder–Mead method the relative “efficiency at the best case” outweighs the absence of a theory of convergence, noticed in [11], was also confirmed by the results of our comparative analysis. As this method is still a favorite among other zero-order optimization methods, it was used for testing the model under consideration based on real experimental data (provided by hydrological stations monitoring pollutants in the Kazakhstani part of the Ili River basin). The analysis of the obtained numerical results was carried out on the basis of the report on the research work of the JSC Institute of Geography and Water Security [13].


## 2. Materials and Methods

### 2.1. Statement of the Inverse Problem

In [4] and in the previous studies, the correctness of the physical formulation of the Streeter-Phelps closed system, describing the process of biochemical oxygen consumption, was shown:(1)dLdt=−k0Cs∙Ct∙Lt,t>0,dCdt=−k0Cs∙Ct∙Lt+k2Cs−Ct,t>0.Here,t is the time, Lt is the concentration of dissolved organic substance, Ct is the concentration of dissolved oxygen, k0 is the rate of biochemical oxygen consumption, Cs is the concentration of oxygen saturation, and k2 is the rate of reaeration.In this work, it is assumed thatk2=0, i.e., the process of decomposition of organic matter, which occurs in a water sample, is placed in a sealed flask or in the ice-covered river channels and reservoirs. For a given T>0 in the space L20,T, we have the following direct problem:(2)dLdt=−k0Cs∙Ct∙Lt,t∈0,T,dCdt=−k0Cs∙Ct∙Lt,t∈0,T,(3)L0=L0,C0=C0,where using the given values of constants k0,L0,C0,Cs∈ℝ+, it is necessary to determine functions L,C∈L20,T satisfying system (2) and (3).Definition 1.
For someT>0, a pair of functions L,C∈L20,T will be called a solution of direct problems (2) and (3) if for any functions ω∈H10,T such that ωT=0, the following equalities hold:(4)L0∙ω0+∫0TLt∙ω′tdt=k0Cs∫0TCt∙Lt∙ωtdt,C0∙ω0+∫0TCt∙ω′tdt=k0Cs∫0TCt∙Lt∙ωtdt.Note 1.
Direct problems (2) and (3) are correct in the sense of Hadamard in the space L20,T for T>0.It is obvious that problems (2) and (3) are reduced to an equivalent system of Volterra integral equations of the second kind, which are known to be correct.To determine the most important indicators of water qualityL0 and k0, the problem inverse to ((2)) and ((3)) is considered: for some T≥5 determine the values L0,k0∈ℝ+, where the functions L,C∈L20,T are a solution to direct problems (2) and (3), using the given values of C0,Cs∈ℝ+ and the additional information:(5)LT=LT,C5=C5.Note 2.
Inverse problems (2) and (5) are incorrect.Indeed, inverse problems (2) and (5) are incorrect, since direct problems (2) and (3) are correct [10].Note 3.
Let for someT≥5 functions L,C∈L20,T be a solution to direct problems (2) and (3), then, there are traces of LT,C5∈ℝ and the following estimates:(6)LT≤L0+k0CsCL20,T∙LL20,T,C5≤C0+k0CsCL20,T∙LL20,T,are correct.


### 2.2. Statement of the Optimization Problem and Methods for Its Solution

In view of Note3, formulated inverse problems (2) and (5) naturally reduce to an equivalent optimization problem associated with minimization of the following objective functional:(7)IL0,k0=LT−LT2+C5−C52.Here,LTandC5 are given values and Lt and Ct are solutions of direct problems (2) and (3), which correspond to parameters L0 and k0.Thus, parametersL0,k0∈ℝ+ must be determined from the condition of minimum of functional (7).A convenient way to solve the optimization problems is the Nelder–Mead method [9], which is an unconditional optimization of the functional in several variables without using its gradient. Since the convergence of the Nelder–Mead method has not been established for the two-dimensional case, for comparison, we also consider the regularizing gradient method, which gives us conditional convergence [10]. The gradient of functional (7) is equal to(8)I′=IL0′,Ik0′,where(9)IL0′=p10,Ik0′=−1Cs∫0TCtLtp1t+p2tdt,functions p1t,p2t are solutions for adjoint system(10)dp1dt=k0CsCtp1t+p2t−2Lt−LTδt−T,t∈0,T,dp2dt=k0CsLtp1t+p2t−2Ct−C5δt−5,t∈0,T,p1T=0,p2T=0.


## 2.1. Statement of the Inverse Problem

In [4] and in the previous studies, the correctness of the physical formulation of the Streeter-Phelps closed system, describing the process of biochemical oxygen consumption, was shown:(1)dLdt=−k0Cs∙Ct∙Lt,t>0,dCdt=−k0Cs∙Ct∙Lt+k2Cs−Ct,t>0.Here,t is the time, Lt is the concentration of dissolved organic substance, Ct is the concentration of dissolved oxygen, k0 is the rate of biochemical oxygen consumption, Cs is the concentration of oxygen saturation, and k2 is the rate of reaeration.In this work, it is assumed thatk2=0, i.e., the process of decomposition of organic matter, which occurs in a water sample, is placed in a sealed flask or in the ice-covered river channels and reservoirs. For a given T>0 in the space L20,T, we have the following direct problem:(2)dLdt=−k0Cs∙Ct∙Lt,t∈0,T,dCdt=−k0Cs∙Ct∙Lt,t∈0,T,(3)L0=L0,C0=C0,where using the given values of constants k0,L0,C0,Cs∈ℝ+, it is necessary to determine functions L,C∈L20,T satisfying system (2) and (3).Definition 1.
For someT>0, a pair of functions L,C∈L20,T will be called a solution of direct problems (2) and (3) if for any functions ω∈H10,T such that ωT=0, the following equalities hold:(4)L0∙ω0+∫0TLt∙ω′tdt=k0Cs∫0TCt∙Lt∙ωtdt,C0∙ω0+∫0TCt∙ω′tdt=k0Cs∫0TCt∙Lt∙ωtdt.Note 1.
Direct problems (2) and (3) are correct in the sense of Hadamard in the space L20,T for T>0.It is obvious that problems (2) and (3) are reduced to an equivalent system of Volterra integral equations of the second kind, which are known to be correct.To determine the most important indicators of water qualityL0 and k0, the problem inverse to ((2)) and ((3)) is considered: for some T≥5 determine the values L0,k0∈ℝ+, where the functions L,C∈L20,T are a solution to direct problems (2) and (3), using the given values of C0,Cs∈ℝ+ and the additional information:(5)LT=LT,C5=C5.Note 2.
Inverse problems (2) and (5) are incorrect.Indeed, inverse problems (2) and (5) are incorrect, since direct problems (2) and (3) are correct [10].Note 3.
Let for someT≥5 functions L,C∈L20,T be a solution to direct problems (2) and (3), then, there are traces of LT,C5∈ℝ and the following estimates:(6)LT≤L0+k0CsCL20,T∙LL20,T,C5≤C0+k0CsCL20,T∙LL20,T,are correct.


## 2.2. Statement of the Optimization Problem and Methods for Its Solution

In view of Note3, formulated inverse problems (2) and (5) naturally reduce to an equivalent optimization problem associated with minimization of the following objective functional:(7)IL0,k0=LT−LT2+C5−C52.Here,LTandC5 are given values and Lt and Ct are solutions of direct problems (2) and (3), which correspond to parameters L0 and k0.Thus, parametersL0,k0∈ℝ+ must be determined from the condition of minimum of functional (7).A convenient way to solve the optimization problems is the Nelder–Mead method [9], which is an unconditional optimization of the functional in several variables without using its gradient. Since the convergence of the Nelder–Mead method has not been established for the two-dimensional case, for comparison, we also consider the regularizing gradient method, which gives us conditional convergence [10]. The gradient of functional (7) is equal to(8)I′=IL0′,Ik0′,where(9)IL0′=p10,Ik0′=−1Cs∫0TCtLtp1t+p2tdt,functions p1t,p2t are solutions for adjoint system(10)dp1dt=k0CsCtp1t+p2t−2Lt−LTδt−T,t∈0,T,dp2dt=k0CsLtp1t+p2t−2Ct−C5δt−5,t∈0,T,p1T=0,p2T=0.


## 3. Results and Discussions

### 3.1. Comparison of Numerical Results for Synthetic Data

To carry out computational experiments, as sought values, we takeL0=20,k0=0.23. Let alsoC0=Cs=9,T=5.Solving direct problems (2) and (3), we find functions Lt,Ct, from which we obtain the values C5=C5,LT=LT, which we take as additional information (5). And, using these data, we will reconstruct the parameters L0,k0 from the condition of functional minimum (7). Using the Nelder–Mead (NM) and gradient (GM) methods, the results presented in Table 1 were obtained. Here, the exact values of C5 and LT are taken as additional information. Figure 1 shows the graphs of the functions Lt,Ct, corresponding to the values of parameters L0,k0 obtained by the NM and GM methods.Table 1
Results of calculations with exact values ofC5 andLT.
Nelder–Mead methodGradient methodL0k0L0k0Reconstructed value20.0000360.23000119.99850.229947Absolute error0.0000360.0000010.0014140.000053Relative error0.00018%0.00043%0.00707%0.02304%Functional value6.3·10−101.0·10−6Figure 1
FunctionsLt and Ct calculated with exact values of LT and C5.As can be seen from Table1, in the case when the exact values of LT and C5 are known, the values of L0 and k0 are reconstructed quite accurately. In particular, the relative NM error was less than a hundredth of a percent.Let us consider the case when the additional information (5) is given with some error α; that is, the values of LT and C5 can randomly deviate from their true value by less than α. After performing a number of computational experiments, it is possible to estimate the range of variation of the sought parameters and the corresponding functions Lt,Ct. Tables 2–5 show the corresponding calculation results for different α. Note that these tables present the maximum error values.Table 2
Results of calculations with an errorα=5%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error0.550.020.550.02Maximum relative error2.75%7.35217%2.76%7.33%Table 3
Results of calculations with an errorα=10%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error1.100.031.100.03Maximum relative error5.50%15.58%5.51%15.55%Table 4
Results of calculations with an errorα=15%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error1.650.051.650.05Maximum relative error8.25%24.85%8.25%24.82%Table 5
Results of calculations with an errorα=20%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error2.200.082.200.08Maximum relative error11.00%35.36%11.00%35.32%Figure2 shows the graphs of the functions Lt,Ct, corresponding to the values of parameters L0,k0 obtained by the NM and GM methods. The solid lines show the functions constructed by the results of solving the inverse problem with exact data for LT and C5. Dashed lines indicate maximum deviations when additional information (5) is specified with 5% error. As can be seen from Table 2 and Figure 2, the numerical results obtained by the NM and GM methods almost coincide.Figure 2
FunctionsLtandCt, α=5%.Similar results forα=10%, α=15%, and α=20% are presented in Tables 3–5 and in Figures 3–5, respectively.Figure 3
FunctionsLt and Ct, α=10%.Figure 4
FunctionsLt and Ct, α=15%.Figure 5
FunctionsLt andCt, α=20%.The results of Tables2–5 show that an increase in the input error causes an increase in the error of the results. Moreover, it can be noted that the relative error for k0 is about 3 times greater than that for L0, while the value of the objective functional always remains sufficiently small. It is also easy to see that the function Ct as a whole is reconstructed much more accurately: in the case when the additional information (5) is given with some error, the deviation of the graph of function Ct from the true one is insignificant. This is likely to be explained by the fact that the initial condition C0 for Ct in the studied problem is known, whereas the initial condition for Ltis not known.Thus, the optimization problem, equivalent to inverse problems (2) and (5), associated with minimization of functional (7), was solved by the Nelder–Mead method and by the gradient method, and the sought parameters L0,k0 were reconstructed quite accurately with the same maximum deviations for data errors. A series of calculations were carried out using field experimental data provided by the Institute of Geography. Oil products were considered pollutants, as they are included in the list of parameters, which must be obligatory determined according to the mandatory program of monitoring the quality of surface waters by hydrochemical and hydrological indicators.


### 3.2. Numerical Results for Experimental Data

Kazakhstani part of the Ili River, which is the main tributary of lake Balkhash basin, was chosen as the object of modeling. In the system of the economy of Kazakhstan, the basin is a diversified economic complex, which has environmentally hazardous enterprises of the extractive industry and nonferrous metallurgy. About 30% of the water resources of the Ili River are formed on the territory of Kazakhstan [14–16]. The Ili River is a transboundary river, and in the recent decades, Kazakhstan has faced the growing water deficit, and one of the reasons for this is the policy of China to increase unilaterally the water intake from transboundary rivers Irtysh and Ili ignoring herewith the interests of the Kazakh side. The chemical composition of the Ili River on the territory of the Kazakhstan to the Kapshagai water reservoir is formed under the influence of pollutants coming from the territory of China as well as contaminated surface runoff and washout from the farmland adjacent to the basin. [17] The data of chemical analysis of water from 4 hydrological stations for the period from 2000 to 2014 were obtained. Table 6 presents the data for these hydrological stations. Figure 6 shows their location on the map.Table 6
The main points of quality control of surface waters in the Ili River basin.
Waterbody codeLocationPost codeDistance from the outfall (km)Catchment area (square km)113200001Pier Dobyn1400372364388113200001164 km above the Kapshagai hydroelectric station140046078540011320000137 km below the Kapshagai hydroelectric station14011434111000113200001Ushzharma village14014264129000Figure 6
Location of hydrological posts on the map (yellow triangles denote hydrological posts),To study the transfer of pollution, as well as to assess the assimilation capacity of the Kazakhstani part of the Ili River, еру data from two hydrological stations were considered: Dobyn pier (the first hydrological station on the Kazakhstani part of the Ili River) and the village of Ushzharma (the last hydrological station before the Ili delta). Figures7 and 8 show the results of usage of two methods for numerical calculations of the sought parameters L0 and k0 at hydrological stations for the period from 2000 to 2014.Figure 7
Initial concentrationL0 and biochemical oxidation rate k0 of oil products at the hydrological station the Dobyn pierFigure 8
Initial concentrationL0 and biochemical oxidation rate k0 of oil products at the hydrological station in the village of UshzharmaFigures7 and 8 show that the solutions of the inverse problem by the Nelder–Mead method and the gradient method give very close results. It can be argued that the constructed mathematical model adequately describes the processes in the Ili River system, as evidenced by the results of a series of numerical experiments. Intra-annual hydrochemical regime of pollution of the Ili River in the section of the Dobyn pier is, in general, consistent with the theoretical description of its natural flow. There is a decrease in pollution in the March, an increase in the July low-water period with an improvement in quality during the period of precipitation in August. Further, during the period of rains, an improvement in quality should be observed, but its quality is deteriorated, which can be explained by the washout of pollutants from the adjacent territories. The increased water content in 2001 and 2008 led to an increase in the oxygen content in the water due to a more rapid flow of water in the river, and as a result of increased aeration. High concentrations of both mineral substances and pollutants, some of which are characterized by high chemical activity, led to more active oxidative reactions, as a result of which, in Figures 7 and 8, extremely high values of L0 and k0 are observed oil products in these years. These data are confirmed by the high values of CIWP (a complex index of water pollution, obtained by averaging all substances that exceed the maximum permissible concentration) in these years. Figure 7 shows that in 2006 the minimum value of the concentration of the pollutant after that the value increases. According to the space images [18], the area of irrigated land in China’s Xinjiang Uygur Autonomous Region (XUAR) had grown to 465,500 hectares. The increase in water consumption affects the river inflow into the territory of Kazakhstan which is one of the factors that affect the increase of concentration of pollutants [16, 18].


## 3.1. Comparison of Numerical Results for Synthetic Data

To carry out computational experiments, as sought values, we takeL0=20,k0=0.23. Let alsoC0=Cs=9,T=5.Solving direct problems (2) and (3), we find functions Lt,Ct, from which we obtain the values C5=C5,LT=LT, which we take as additional information (5). And, using these data, we will reconstruct the parameters L0,k0 from the condition of functional minimum (7). Using the Nelder–Mead (NM) and gradient (GM) methods, the results presented in Table 1 were obtained. Here, the exact values of C5 and LT are taken as additional information. Figure 1 shows the graphs of the functions Lt,Ct, corresponding to the values of parameters L0,k0 obtained by the NM and GM methods.Table 1
Results of calculations with exact values ofC5 andLT.
Nelder–Mead methodGradient methodL0k0L0k0Reconstructed value20.0000360.23000119.99850.229947Absolute error0.0000360.0000010.0014140.000053Relative error0.00018%0.00043%0.00707%0.02304%Functional value6.3·10−101.0·10−6Figure 1
FunctionsLt and Ct calculated with exact values of LT and C5.As can be seen from Table1, in the case when the exact values of LT and C5 are known, the values of L0 and k0 are reconstructed quite accurately. In particular, the relative NM error was less than a hundredth of a percent.Let us consider the case when the additional information (5) is given with some error α; that is, the values of LT and C5 can randomly deviate from their true value by less than α. After performing a number of computational experiments, it is possible to estimate the range of variation of the sought parameters and the corresponding functions Lt,Ct. Tables 2–5 show the corresponding calculation results for different α. Note that these tables present the maximum error values.Table 2
Results of calculations with an errorα=5%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error0.550.020.550.02Maximum relative error2.75%7.35217%2.76%7.33%Table 3
Results of calculations with an errorα=10%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error1.100.031.100.03Maximum relative error5.50%15.58%5.51%15.55%Table 4
Results of calculations with an errorα=15%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error1.650.051.650.05Maximum relative error8.25%24.85%8.25%24.82%Table 5
Results of calculations with an errorα=20%.
Nelder–Mead methodGradient methodL0k0L0k0Maximum absolute error2.200.082.200.08Maximum relative error11.00%35.36%11.00%35.32%Figure2 shows the graphs of the functions Lt,Ct, corresponding to the values of parameters L0,k0 obtained by the NM and GM methods. The solid lines show the functions constructed by the results of solving the inverse problem with exact data for LT and C5. Dashed lines indicate maximum deviations when additional information (5) is specified with 5% error. As can be seen from Table 2 and Figure 2, the numerical results obtained by the NM and GM methods almost coincide.Figure 2
FunctionsLtandCt, α=5%.Similar results forα=10%, α=15%, and α=20% are presented in Tables 3–5 and in Figures 3–5, respectively.Figure 3
FunctionsLt and Ct, α=10%.Figure 4
FunctionsLt and Ct, α=15%.Figure 5
FunctionsLt andCt, α=20%.The results of Tables2–5 show that an increase in the input error causes an increase in the error of the results. Moreover, it can be noted that the relative error for k0 is about 3 times greater than that for L0, while the value of the objective functional always remains sufficiently small. It is also easy to see that the function Ct as a whole is reconstructed much more accurately: in the case when the additional information (5) is given with some error, the deviation of the graph of function Ct from the true one is insignificant. This is likely to be explained by the fact that the initial condition C0 for Ct in the studied problem is known, whereas the initial condition for Ltis not known.Thus, the optimization problem, equivalent to inverse problems (2) and (5), associated with minimization of functional (7), was solved by the Nelder–Mead method and by the gradient method, and the sought parameters L0,k0 were reconstructed quite accurately with the same maximum deviations for data errors. A series of calculations were carried out using field experimental data provided by the Institute of Geography. Oil products were considered pollutants, as they are included in the list of parameters, which must be obligatory determined according to the mandatory program of monitoring the quality of surface waters by hydrochemical and hydrological indicators.


## 3.2. Numerical Results for Experimental Data

Kazakhstani part of the Ili River, which is the main tributary of lake Balkhash basin, was chosen as the object of modeling. In the system of the economy of Kazakhstan, the basin is a diversified economic complex, which has environmentally hazardous enterprises of the extractive industry and nonferrous metallurgy. About 30% of the water resources of the Ili River are formed on the territory of Kazakhstan [14–16]. The Ili River is a transboundary river, and in the recent decades, Kazakhstan has faced the growing water deficit, and one of the reasons for this is the policy of China to increase unilaterally the water intake from transboundary rivers Irtysh and Ili ignoring herewith the interests of the Kazakh side. The chemical composition of the Ili River on the territory of the Kazakhstan to the Kapshagai water reservoir is formed under the influence of pollutants coming from the territory of China as well as contaminated surface runoff and washout from the farmland adjacent to the basin. [17] The data of chemical analysis of water from 4 hydrological stations for the period from 2000 to 2014 were obtained. Table 6 presents the data for these hydrological stations. Figure 6 shows their location on the map.Table 6
The main points of quality control of surface waters in the Ili River basin.
Waterbody codeLocationPost codeDistance from the outfall (km)Catchment area (square km)113200001Pier Dobyn1400372364388113200001164 km above the Kapshagai hydroelectric station140046078540011320000137 km below the Kapshagai hydroelectric station14011434111000113200001Ushzharma village14014264129000Figure 6
Location of hydrological posts on the map (yellow triangles denote hydrological posts),To study the transfer of pollution, as well as to assess the assimilation capacity of the Kazakhstani part of the Ili River, еру data from two hydrological stations were considered: Dobyn pier (the first hydrological station on the Kazakhstani part of the Ili River) and the village of Ushzharma (the last hydrological station before the Ili delta). Figures7 and 8 show the results of usage of two methods for numerical calculations of the sought parameters L0 and k0 at hydrological stations for the period from 2000 to 2014.Figure 7
Initial concentrationL0 and biochemical oxidation rate k0 of oil products at the hydrological station the Dobyn pierFigure 8
Initial concentrationL0 and biochemical oxidation rate k0 of oil products at the hydrological station in the village of UshzharmaFigures7 and 8 show that the solutions of the inverse problem by the Nelder–Mead method and the gradient method give very close results. It can be argued that the constructed mathematical model adequately describes the processes in the Ili River system, as evidenced by the results of a series of numerical experiments. Intra-annual hydrochemical regime of pollution of the Ili River in the section of the Dobyn pier is, in general, consistent with the theoretical description of its natural flow. There is a decrease in pollution in the March, an increase in the July low-water period with an improvement in quality during the period of precipitation in August. Further, during the period of rains, an improvement in quality should be observed, but its quality is deteriorated, which can be explained by the washout of pollutants from the adjacent territories. The increased water content in 2001 and 2008 led to an increase in the oxygen content in the water due to a more rapid flow of water in the river, and as a result of increased aeration. High concentrations of both mineral substances and pollutants, some of which are characterized by high chemical activity, led to more active oxidative reactions, as a result of which, in Figures 7 and 8, extremely high values of L0 and k0 are observed oil products in these years. These data are confirmed by the high values of CIWP (a complex index of water pollution, obtained by averaging all substances that exceed the maximum permissible concentration) in these years. Figure 7 shows that in 2006 the minimum value of the concentration of the pollutant after that the value increases. According to the space images [18], the area of irrigated land in China’s Xinjiang Uygur Autonomous Region (XUAR) had grown to 465,500 hectares. The increase in water consumption affects the river inflow into the territory of Kazakhstan which is one of the factors that affect the increase of concentration of pollutants [16, 18].


## 4. Conclusions

On the basis of the closed Streeter-Phelps system, an inverse problem that determines the power of point sources of pollutants, at which the concentration of pollutants at hydrological posts is equal to the maximum permissible, has been formulated. A comparative analysis of optimization methods (NM and GM) for solving the formulated inverse problem was carried out, which showed the effectiveness of application of the Nelder–Mead method and the gradient method. The results of numerical experiments made it possible to assess the assimilation capacity of the Kazakhstani part of the Ili River basin. The resulting estimate determines the upper limit of the assimilation capacity of the basin (the highest seasonal value of the maximum permissible load), since the calculations for the BOD were carried out atk=0.23 day-1. This value corresponds to a water temperature of 20°С, which is usually observed during the summer low-water period. It is obvious that in winter, when the rate of decomposition of pollutants is significantly lower, the calculated value of the maximum permissible load will decrease. Also, the obtained results can be used to predict changes in the concentration of pollutants in case a decrease in the river inflow into the territory of Kazakhstan.



---

*Source: 9387827-2021-08-06.xml*