Научная статья на тему 'Results of computing experiments for water ecological system mathematical modelin'

Results of computing experiments for water ecological system mathematical modelin Текст научной статьи по специальности «Математика»

CC BY
75
19
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
MATHEMATICAL MODEL / MATHEMATICAL MODELING OF AQUATIC ECOSYSTEMS / COMPUTING EXPERIMENT

Аннотация научной статьи по математике, автор научной работы — Petrov Y. S., Raspopov V. E.

The point-wise imitation and one-dimensional mathematical models of aquatic ecosystems have been overlooked. The developed models are intended for studying ecosystems in the Krasnoyarsk aquatic basin and in separate locations on the Yenisei River. The results of the computing experiments are presented.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Results of computing experiments for water ecological system mathematical modelin»

Disadvantages:

- algorithms are complicated;

- there are high system requirements.

Except for the algorithm of segmentation itself in this work there were described some other new ideas:

- the method of multi scale analysis with extraction of frequency information;

- the algorithm of two-dimensional function minimization which uses morphological filling;

- combination of areas growing and areas uniting;

- the criterion of the definition of the optimal moment to stop uniting.

References

1. Gonzalez R. C., Woods R. E. Digital Image Processing. N. Y. : Prentice Hall, 2002.

2. Steinbrecher R. Bildverarbeitung in der Praxis

[Electronic resource]. URL: http://www.rst-

software.de/dbv/download.html.

© Palamar I. N., Sizov P. V., 2010

Y. S. Petrov, V. E. Raspopov Siberian Federal University, Russia, Krasnoyarsk

RESULTS OF COMPUTING EXPERIMENTS FOR WATER ECOLOGICAL SYSTEM

MATHEMATICAL MODELING

The point-wise imitation and one-dimensional mathematical models of aquatic ecosystems have been overlooked. The developed models are intended for studying ecosystems in the Krasnoyarsk aquatic basin and in separate locations on the Yenisei River. The results of the computing experiments are presented.

Keywords: mathematical model, mathematical modeling of aquatic ecosystems, computing experiment.

Environmental issues have a designated place in the general list of issued for which mathematical modeling is used. The increase of the anthropogenic environmental impact, caused by intense exploitation of natural resources and the growth of industry leads to an ecological balance disruption. This is happening both on local (in separate areas of globe) and on planetary scale. The importance of struggling against anthropogenic eutrophication of reservoirs and their pollution is understood everywhere in the world. There had been a great amount of researches in limnology, mathematical modeling, and economy, connected with problems of preservation, restoration, and the effective exploitation of natural resources, such as lakes and manmade reservoirs. The ecological condition of the water bodies depends on a number of various factors and processes: hydrophysical, hydrobiological, hydrochemical, meteorological, and anthropogenic. Hydrophysical processes appreciably form a habitat of hydrobionts, define the transferred and sedimentation of substances, the intensity of pollution, and the self-cleaning of reservoirs.

The problem of water quality is complicated. Water bodies are complex physical, biochemical and ecological systems. To be able to predict the consequence of one decision or another, the corresponding tool by dint of which it is possible to analyze the sufficiency of information is required. Such a tool is the computing experiment based on mathematical modeling and numerical methods. An effective means of the arising problem objective analysis in the field of hydrobiology problems are the methods based on constructing and studying mathematical models of water ecosystems. The using of mathematical modeling and carrying out computing experiments allows us to predict the dynamics

of water ecosystem development, and also to estimate the consequences of realizing various projects, connected with influence on the ecosystem.

A number of general claims to each mathematical model are known: the corresponding system of the equations should be closed and consistent; the model should describe a variety of physical phenomena and suppose the designing of realized numerical algorithm.

In the given work, some results of the calculations, carried out with a mathematical model of the water ecosystem (being an improvement of the model considered in [1]) are presented. The model is modified by the separation of green algae as independent components of a mathematical model and the introduction of an additional equation, describing the change in algae concentration.

As dynamic variables of model, the concentrations of green algae (CA0), of blue-green algae (CA1), of diatoms (CA2), of zooplankton (CZ), of bacteria (CB), of detritus (CD), of the inorganic phosphorus dissolved in water (PS), of the inorganic nitrogen dissolved in water (NS), of the organic matter dissolved in water (POB), and of the oxygen dissolved in water (02) are taken.

In model the following processes are considered:

- growth of microorganisms;

- outflow of products of a metabolism;

- death rate of microorganisms;

- processes of settling;

- transitions on a trophic chain;

- decomposition processes;

- atmospheric reaeration (isolation of oxygen from water);

- denitrification (process of restoration of nitrates to the molecular nitrogen, caused by bacteria);

- limiting factors (illumination, temperature);

- water aeration (saturation of water by oxygen of air).

The main feature of the given model is the division of

blue-green algae into two species: greens and blue-greens; this isn’t presented in many models, but is of great importance for the research of the reservoirs’ ecology, for their development is various.

The model allows predicting the dynamics of water ecosystem development; including the transformations of nitrogen and phosphorus, as basic biogenous elements, defining the efficiency and water quality in reservoirs.

The structure of model describing the functioning of an ecosystem is given in the flow chart (fig. 1).

Fig. 1. Flow chart of the model. Arrows designate directions of substance streams between ecosystem components

On the basis of the flow chart the imitating model, describing the considered processes in the ecosystem is constructed. The mathematical model represents the following system of the ordinary differential equations with corresponding initial data:

dCA0

dt

dCA1

dt

dCA2

= (mA0 - RA0 - MA0) ■ CA0 + a0 • CA1- CA0, = (mA1 - RA1 - MAY) ■ CA1 - a1 ■ CA1 ■ CA0,

m7 ■ C7

-----= (mA2 - RA2 - SA2 -MA2) ■ CA2--------------,

dt Y1

dC7

dt

= (m7 - R7 - M7) ■ C7,

dCB = (mB - RB - MB) ■ CB - m7 ^C7

dt

Y 2

dCD

dt

= MA0 ■ CA0 + MA1 ■ CA1 + MA2 ■ CA2 + M7 ■ C7 +

+MB■ CB -SA3*CD-, (1)

Y 3

Y 4

dPS

dt

= -(mA0 - RA0) ■ PP0 ■ CA0 - (mA1 - RA1);

< PP1 ■ CA1 -(mA2 - RA2) ■ PP2 ■ CA2 +

+ R7 ■ C7 ■ PP3 + RB ■ CB ■ PP4,

dNS

dt

= RA0 ■ PN0 ■ CA0 + RA1 ■ PN1 ■ CA1 - (mA2 - RA2) >

x PN2 ■ CA2 + R7 ■ C7 ■ PN3 + RB ■ CB ■ PN4, dPOB mB ■CB + h0 ■ RA0 ■ CA0 + hb RAb CA1 +

dt Y 5

+h2 ■ RA2 ■ CA2 + h3 ■ R7 ■ C7 + h4 ■ RB ■ CB, dO2

dt

• = K1 ■ (O2O - O2) + Kacc ■ (mA0 ■ CA0 + mA1x

x CA1 + mA2 ■ CA2) -alf ■ (RA1 ■ CA1 + RA2 ■ CA2 +

+ R7 ■ C7 + RB ■ CB) - Bb m7 ■ C7.

where mAi are functions describing growth; coefficients RAi are breath; MAi are the death rate; SAi is the settling; Yi are proportionality coefficients; T in temperature in C0; t is time.

In the live description the incoming and proceeding streams’ components are taken into account. Also is included the share of received resources (food) spent for growth and reproduction; proceeding is the consumption of species from given components; predators and death rate depending on every other possible reason. Meanwhile, the influence on the stream speed of the environment (temperature, etc.) is considered.

In microbiological systems as a rule, the growth rate is limited by a concentration of substrates. We have applied the hyperbolic dependence offered by 7 mono for the description of the limitation process.

It is supposed, that the growth of green and blue-green algae is limited by phosphorus, while the growth of diatoms - by nitrogen and phosphorus. The growth functions, death rate, illumination, and temperature dependence, as well as all entrance data are included according to researches [1-3].

The constructed mathematical model represents the Cauchy problem for a system of ten ordinary differential equations. For the numerical solution of the Cauchy problem, the Runge-Kutta method of the fourth approximation order is applied:

yB+1 =y„+ 1 T( K1 + 2 K 2 + 2 K3 +K4),

6

K1 = F(tn,yn),

K 2 = F (tn + 2, yn + T yX K3 = F (tn + 2, y„ + T K2l),

K4 = F (tn + t, yn + тKз), n = ^ 1, ^,

where y is a vector function of unknown; F is the right part of system (1); t is a step in time; y0 is specified.

Let’s note that the set of components in the model considerably complicates the problem, both the modeling, and in studying the model; as it is required to specify its value for each coefficient (fig. 2).

A complex of the programs is written, allowing the inputting of entrance data in an interactive mode. The calculation results can be received numerically, presented

graphically, and transferred outside for subsequent processing. For the management of graphic representation of calculation results, a corresponding menu is provided.

August decreases to 1.2-2 mg/l. Model calculation has also shown a falling in values of biomass during this period (fig. 3).

Fig. 2. Coefficients of modeling and initial system values

The program’s complex is realized in Visual C ++ 6,0 with the use of MFC (Microsoft Foundation Classes) which is one of the most convenient and powerful tools among Windows’ applications. The software product has a friendly interface, it is convenient in work, and is intended not only for mathematicians, but also for researchers who are not experts in programming.

The first numerical experiments have been devoted to the comparative analysis of results, obtained by means of a working mathematical model [1] and by its updating means (1) with same input data [4]. The calculation results have shown that the received concentrations of diatoms, bacteria, and detritus for a working model [1] poorly correspond to experimental data, in comparison with the results obtained from the model aforementioned.

Thus, the computing experiments that have been carried out have shown the effectiveness of separating green algae as an independent variable of the mathematical model for the reservoir ecosystem.

The following calculations with an improved mathematical model are meant for researching general tendencies of seasonal dynamic variable change for a model, using field data from the Novoselovsky reach of the Krasnoyarsk impoundment for 1998-2000. Notice that for the comparison of calculations results, we have used only the field data, the time moments of which are precisely known. During other time periods, due to the incompleteness of existing information, the comparative analysis of average data also shows qualitative calculation coincidence.

Particularly seasonal dynamics of diatoms demonstrate a qualitative picture of two “flowering” peaks: summer - with the maximum biomass of 5.9 mg/l, and autumn with the maximum biomass of 2.27 mg/l; this corresponds with the supervision data [5]. The total biomass of diatoms according to supervisions in July and

Fig. 3. Comparison of experimental data with numerical results for diatoms. Natural data is marked by triangles

The comparison of a seasonal course model for green and blue-green algae and experimental data [6], has shown that in a general understating of model concentration (approximately by 1.5 times) the relative time course had been precisely reconstructed (fig. 4). The annual course of biomass for zooplankton has a single peak and corresponds to the maximum values for green and blue-green algae biomasses, which also correspond with theoretical representations.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Fig. 4. Comparison of experimental data with numerical results for blue-green algae. Natural data are marked by triangles

For a seasonal course of chemical combinations of nitrogen and phosphorus concentrations, the calculation has shown the falling of values of nitrogen compound concentration during seasons, when intensive plankton growth occurs, and the maximum values at minima of plankton biomass.

The annual course of organic matter concentration has two expressed maxima with values of 0.29 mg/l in the beginning of the summer; and 0.54 mg/l in the autumn. These periods correspond to the maximum values of

phytoplankton and zooplankton biomasses. For seasonal dynamics of bacteria and detritus, two peaks of development are typical: the first falls in the middle of the summer, the second - in the beginning of autumn. Thus, detritus influence the growth of bacteria and stimulate their productivity, which also corresponds to natural data.

The obtained calculation model coincides with experimental data, which testifies the adequacy of the examined model.

The one-dimensional model of the aquatic ecosystem. Along with point-wise model the mathematical model allowing mass transfer along the length of a reservoir (one-dimensional model in a horizontal plane) is examined. The mathematical model represents a differential equation system in partial derivatives of the first degree:

dU dU ~ ~

+ V----------= F (t, x,U) (2)

dt dx

with corresponding initial and edge conditions:

U (t ,0) = U0(t),

U (0, x) = Uj(x).

where V is the current speed of the reservoir; x is the spatial variable corresponding to the length of a reservoir; t is time; the right part of the equations (2) corresponds to the right part (1). In such modeling it is supposed, that the substance is evenly distributed along the width of the stream and moves with the average speed of the stream. So, the data is averaged for the depth and width of a

reservoir. We will notice that the given model is

expedient for using in the case when the length of a

reservoir is greater than its width.

The introduced mathematical model is also realized numerically by means of an implicit difference scheme:

an insignificant role. Bacterial mass cumulates which here, undergoes intense decomposition only at the inlet stream of the Angara.

> n+1

yj - y

j + Vyj -yj-1

T

h

= Fj.

For initial data, the data from the point-wise model is taken. For edge conditions - the solutions obtained from the point-wise model are used.

Initial distribution of all system components is considered uniform. The calculation was carried out for time from t = 0 to t = 365 (one year), for the Novoselovsky reach of the Krasnoyarsk impoundment, and also for sites on the Yenisei River downstream from the Krasnoyarsk Hydroelectric Power Plant with average current speed of the Yenisei of 1.2 km/h on a distance from ten to one hundred kilometers. In fig. 5 particularly, the results for bacteria calculations on a river site are given. It is visible, that under the specified conditions a concentration of bacteria in the chosen part of a reservoir changes considerably (we suppose that this difference of values is caused by the current).

Notice that the ecosystem of the Yenisei on the site adjoining the power plant is strongly impoverished because of the destroying action of the plant’s turbines and the low water temperature. The self-cleaning process of water in this heavily polluted site is weakened. The zooplankton in process of substance decomposition plays

Fig. 5. Results of calculations for the bacteria carried out in a one-dimensional mathematical model

Calculations have shown that the further the distance is from the Krasnoyarsk Hydroelectric Power Plant, the more change occur in concentration, selected as a component; in particular, there is a shift of peaks for diatoms, bacteria, and detritus. At a distance from the power plant of over one hundred kilometers, the model depicts the dynamics of water ecosystem development in space less adequately. This is caused by the studying of quality water characteristics; it is necessary to take in account the more difficult and full processes: hydrodynamic (convective stirring, pressure, wind, and deep currents), heat transfer and illumination. The use of the models based on various variants of mechanic equations for liquid and heat transfer, and corresponding boundary conditions is necessary [7].

Let’s note the basic results obtained in the work:

Point-wise and one-dimensional mathematical models of water ecosystems consisting of ten differential equations had been constructed. The use of models allows revealing of processes progressing dynamics in difficult ecological water systems, to predict a system status in time and in space (the one-dimensional model in a horizontal plane) for a distance up to one hundred kilometers, and to analyze problematic situations. The models in particular, make it possible to describe the change of hydrobionts and the basic biogenic elements, and also to reproduce occurrence situations, both for one and for two peaks of phytoplankton flowering during the vegetative season, depending on external conditions. It is necessary to notice, however, that the offered mathematical models are very sensitive to changes in parameters and demand a meticulous selection of coefficients for each specific aquatic ecosystem.

The results obtained by means of the models described above, can as well be used as well for estimating ecological risks.

A complex of computer programs has been produced, and the numerical simulation of some processes progressing in the ecosystem of the Novoselovsky reach of the Krasnoyarsk impoundment and in the Yenisei River had been carried out.

The authors express their gratitude to Professor Z. G. Gold for his useful consultations.

References

1. Mathematical simulation of reaches of the Krasnoyarsk impoundment / V. A. Sapozhnikov [et al.] // Association of subjects of the Russian Federation and a wildlife management problem in Prienisejsky Siberia : theses and materials of reports of inter-regional scientifically-practical conference / KSU. Krasnoyarsk, 2005. P. 296-298.

2. Borodin, A. L., Raspopov V. E. Numerical identification of coefficients of mathematical model of an ecosystem of a reservoir // Joint issue. Computing technologies. T. 13. Herald of the KazSU of AL-FARABI. A series of the mathematics, the mechanics, the informatics. Vol. 3 (58). Almaty-Novosibirsk, 2008. P. 302-306.

3. Gubanov V. G. Biotic circulation and interaction of trophic links in artificial and natural biosystems : diss. dr. phys.-math. sciences. Krasnoyarsk, 2004.

4. Petrov J. S. Special-purpose software for carrying out of computing experiments at mathematical modeling of water ecosystems // YOUTH AND THE SCIENCE: the XXI-st century BEGINNING : Materials of the All-Russia scientific and technical conf. of students, postgraduate students and young scientists. In 4 p. P. 1 / SFU. Krasnoyarsk, 2009. P. 78-80.

5. Kozhevnikova N. A., Phytoplankton of a deep-water part of a Krasnoyarsk impoundment // Alkologia. № 2. 2002. P. 39-40.

6. ShChur L. A. Structure and functional characteristics of bacterial plankton and phytoplankton in ecosystems of reservoirs of different type : diss. dr. biol. sciences. Krasnoyarsk, 2006.

7. Belolipetsky V. M., Genova S. N., Gurevich K. J. Platform for research of dynamics of hydrophysical and radio ecological characteristics of river system // Computing technologies / the Siberian Branch of the Russian Academy of Science. Vol. 6, № 2. 2001. P. 14-24.

© Petrov Y. S., Raspopov V. E., 2010

V. V. Podkolzin, V. O. Osipyan Kuban State University, Russia, Krasnodar

ON PROPERTIES OF KNAPSACK SYSTEMS OF INFORMATION PROTECTION WITH THE OPEN KEY IN Zp

Properties of sequences of numbers expressed through components of a knapsack vector are investigated. Properties of isomorphic and similar knapsack systems of information protection are analyzed. Methods of increasing cryptographic security of knapsack systems of information protection with an open key are presented.

Keywords: a knapsack vector, isomorphism, cryptoanalysis, density, injectivity.

Let’s express a set of natural numbers {0, 1, ..., p-1} through Zp and a set of all numerical sets of length n with

components from Zp. through ZUp .

A knapsack problem for set w e N and vector A = (a1,

a2, ..., an), where aieN, I = 1...n, if there is an equation solution

AxT = w

has the solution in Zp

(1)

j-1

• >Y, (p -1) ai

i-1

j

2... n

(2)

A knapsack vector A = (a1, a2, ..., an) is a nondecreasing one if its components are ordered according to the rule ai-1 < a„ I = 2. n. Accordingly, the vector is increasing if its components are ordered according to the rule ai-1 <a„ I = 2. n.

Definition. Let’s call vector AA = (8 i, 82, ..., 8p) a variation of vector A = (a1, a2, ..., an) (ai e N, I = 1...n) in Zp, For its components the following correlations are carried out:

we will call vector A of equations (1) a knapsack vector.

A knapsack vector A = (a1, a2, ..., an) is an injective one if for any natural w the equation (1) has not more than one solution. A knapsack vector which has two elements ai = aj, I £ j, is not injective. Injectivity of a knapsack vector allows to speak about uniqueness of restoration of the original text according to the cryptogram. Supergrowing knapsack vectors are the simplest injective knapsack vectors from the point of view of understanding and algorithmization. For their components in Zp the following relationships are carried out:

j-1

81 - a1, 8j = aj (p - 1)a<, j = 2-n.

(3)

On the basis of vector AA it is possible to define a knapsack vector A in Zp corresponding to it:

a1 - 81, i-1 i-1

at - 8,- + (p - 1)X aj - 8,- + (p -1)£ p‘-j-18j,

j-1 j-1

I = 2-n. (4)

Let’s express a set of various values w for which equation (1) has the solution through |a (p, A). Capacity ^ (p, A) does not exceed pn since the quantity of various

i-1

i Надоели баннеры? Вы всегда можете отключить рекламу.