Theoretical analysis method and apparatus for simulating the transition distribution of multi-size microcystis colonies

Abstract

To solve the technical problem of the difficulty in accurately describing transition distribution process of multi-sized Microcystis colonies.SOLUTION: The present disclosure provides a theoretical analysis method and device for simulating transition distribution of multi-sized Microcystis colonies. The method comprises: classifying all Microcystis colonies into multiple groups based on size; constructing an initial concentration control equation for each group and simplifying it into a target concentration control equation; setting initial and boundary conditions for a concentration; converting these into a dimensionless form; and further analytically calculating a dimensionless concentration control equation based on target quantity moments and central moments, additional conditions, initial conditions, and boundary conditions for each group, to obtain an analytical expression for a spatial concentration distribution of all Microcystis colonies. Consequently, it is possible to precisely describe a transition distribution process of multi-sized Microcystis colonies.SELECTED DRAWING: Figure 1

Description

[Technical Field]

[0001] This disclosure relates to the field of aquatic ecological environment technology, and more particularly to a theoretical analysis method and apparatus for simulating the transition distribution of multi-sized Microcystis colonies. [Background technology]

[0002] Microcystis can be divided into single cells and colonies. In natural environments, microcystis often exists in colonial form and varies in size. The spatial distribution of microcystis determines the dynamics of the occurrence and disappearance of water blooms. In eutrophic water bodies, microcystis can usually reach high biomass, so the succession of microcystis in water bodies is a crucial process for water bloom generation. Currently, research methods on the succession distribution of multi-sized microcystis colonies include field observations, laboratory experiments, and numerical simulations, but theoretical analysis methods are lacking. Therefore, it is difficult to accurately describe the succession distribution process of multi-sized microcystis colonies and to provide reliable support for water bloom prediction and forecasting. [Overview of the Initiative] [Problems that the invention aims to solve]

[0003] To solve the above technical problems, or at least partially solve them, this disclosure provides a theoretical analysis method and apparatus for simulating the transition distribution of multi-size microcystis colonies, which differs from theoretical analysis methods for the transition distribution of single-size microcystis colonies. This disclosure can accurately describe the transition distribution process of multi-size microcystis colonies, reflecting the influence of different environmental gradients on the transition characteristics and transport features of microcystis colony patches, and provides reliable support for the prediction and forecasting of water blossoms. [Means for solving the problem]

[0004] The embodiments of this disclosure provide a theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies, the method being: The steps include dividing all microcystis colonies into multiple groups of microcystis colonies according to their size, and obtaining the quantity percentage of microcystis colonies in each group relative to the total quantity of all microcystis colonies, The steps include constructing an initial concentration control equation for each group of microcystis colonies, performing a simplification process on the initial concentration control equation based on a two-dimensional wind flow velocity pattern to obtain a target concentration control equation for each group of microcystis colonies, Based on the quantity percentage of microcystis colonies in each group and the total quantity, the initial concentration conditions for the microcystis colonies in each group are determined, and the concentration boundary conditions for the microcystis colonies in each group are obtained. The steps include: transforming the target concentration control equation, the initial concentration conditions, and the concentration boundary conditions based on the introduction of dimensionless parameters to obtain the dimensionless concentration control equation, the initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis colonies; The steps include: processing the dimensionless concentration control equation for each group of microcystis colonies based on the target number of-order concentration moment analysis equation, the target quantity center concentration moment analysis equation, the preset concentration distribution addition conditions, the dimensionless initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis colonies to obtain the target number of-order concentration moments and the target quantity center concentration moments for each group of microcystis colonies; The method includes the step of calculating, based on the target number-order concentration moment and target number-center concentration moment of each group of microcystis colonies, to obtain a spatial concentration distribution analysis formula for all microcystis colonies.

[0005] Selectively construct the following initial concentration control equations for each group of microcystis colonies: JPEG0007923361000001.jpg24170 Here, C iis the concentration of the i-th group Microcystis colonies, i is a positive integer, t is time, μ is the vector flow velocity, and μ b i is the vector transition velocity of the i-th group Microcystis colonies, λ is the diffusion coefficient of water flow, ∇ is the gradient operator, simplifying the initial concentration control equation based on a two-dimensional wind-driven flow velocity pattern to obtain the target concentration control equation of Microcystis colonies of each group as follows, JPEG0007923361000002.jpg28170Here, μ is the longitudinal flow velocity, and w b i is the transition velocity of the i-th group Microcystis colonies themselves, and λ b ixx is the diffusion coefficient along the x-direction of the i-th group Microcystis colonies, and λ b izz is the diffusion coefficient along the z-direction of the i-th group Microcystis colonies.

[0006] Optionally, an initial time is set, the Microcystis colony patches are concentrated at x=0 and z=0, and based on the quantity percentage of Microcystis colonies of each group and the total quantity, the initial concentration conditions of Microcystis colonies of each group are determined as follows, JPEG0007923361000003.jpg25170Here, N i is the quantity percentage of the i-th group Microcystis colonies, Q is the total quantity of Microcystis colonies, δ is the Dirac function, At x=±∞, z=0, and z=H (where H is water depth), the concentration boundary conditions of Microcystis colonies of each group are JPEG0007923361000004.jpg50170It is used to indicate that the flux of Microcystis colonies at the free surface and the bottom is zero.

[0007] Optionally, the dimensionless parameters introduced are Based on the introduction of dimensionless parameters, the target concentration control equation, the initial concentration conditions, and the concentration boundary conditions are transformed to obtain the dimensionless concentration control equation, the initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis colonies, which are as follows: JPEG0007923361000006.jpg94170i reflects the relative strength of longitudinal transition and total longitudinal diffusion in the microcystis colony of group i.

[0008] Selectively, the target-order-level concentration moment analysis formula for each group of microcystis colonies is: The filename is JPEG0007923361000007.jpg24170. The target quantity central concentration moment analysis formula for each of the aforementioned groups of microcystis colonies is: The second-order central concentration moment, third-order central concentration moment, and fourth-order central concentration moment obtained from JPEG0007923361000008.jpg83170 are as follows: The filename is JPEG0007923361000009.jpg36170.

[0009] The pre-set concentration distribution conditions can be selected. The filename is JPEG0007923361000010.jpg21170. Based on the target order-order concentration moment analysis formula for each group of microcystis colonies, the pre-set concentration distribution addition conditions, the dimensionless initial concentration conditions, and the dimensionless concentration boundary conditions, the dimensionless concentration control equation for each group of microcystis colonies is processed. JPEG0007923361000011.jpg109170 p>0 for the i-th group of Microcystis colonies as follows: JPEG0007923361000012.jpg77170

[0010] Optionally, the step of calculating based on the target number-order concentration moment and target number central concentration moment of the Microcystis colonies of each group to obtain an analytical expression for the spatial concentration distribution of all Microcystis colonies is: Based on the target number-order concentration moment of the Microcystis colonies of each group, applying the Chatwin asymptotic expansion formula to obtain an analytical expression for the spatial concentration distribution of all Microcystis colonies as follows: JPEG0007923361000013.jpg70170where the polynomial coefficient a is determined according to the target number central concentration moment n .

[0011] Embodiments of the present disclosure further provide a theoretical analysis apparatus for simulating the transition distribution of multi-size Microcystis colonies, the apparatus comprising: a grouping and obtaining module, configured to divide all Microcystis colonies into a plurality of groups of Microcystis colonies according to size, and obtain the quantity percentage of each group of Microcystis colonies relative to the total quantity of all Microcystis colonies; a constructing and simplifying module, configured to construct an initial concentration control equation for each group of Microcystis colonies, perform simplification processing on the initial concentration control equation based on a two-dimensional wind-induced flow velocity pattern, and obtain a target concentration control equation for each group of Microcystis colonies; a determining and obtaining module, configured to determine initial concentration conditions for each group of Microcystis colonies based on the quantity percentage of each group of Microcystis colonies and the total quantity, and obtain concentration boundary conditions for each group of Microcystis colonies; a conversion module, configured to convert the target concentration control equation, the initial concentration conditions and the concentration boundary conditions based on introduction of dimensionless parameters, and obtain dimensionless concentration control equation, dimensionless initial concentration conditions and dimensionless concentration boundary conditions for each group of Microcystis colonies; a processing module configured to process the dimensionless concentration control equations of Microcystis colonies in each group based on the target number-order concentration moment analytical expression, the target number-center concentration moment analytical expression of the Microcystis colonies in each group, the preset additional concentration distribution conditions, the initial dimensionless concentration condition and the dimensionless concentration boundary condition, so as to obtain the target number-order concentration moment and the target number-center concentration moment of the Microcystis colonies in each group; a calculation module configured to perform calculation based on the target number-order concentration moment and the target number-center concentration moment of the Microcystis colonies in each group, so as to obtain an analytical expression of the spatial concentration distribution of all Microcystis colonies.

[0012] Embodiments of the present disclosure provide an electronic device, the electronic device comprising a processor and a memory for storing executable instructions executable by the processor, wherein the processor is configured to read the executable instructions from the memory and execute the instructions to implement the theoretical analysis method for simulating the transition distribution of multi-size Microcystis colonies provided in the embodiments of the present disclosure.

[0013] Embodiments of the present disclosure further provide a computer-readable storage medium storing a computer program for executing the theoretical analysis method for simulating the transition distribution of multi-size Microcystis colonies provided in the embodiments of the present disclosure. [Effects of the Invention]

[0014] The technical solutions provided in the embodiments of this disclosure have the following advantages compared to the conventional technology. The theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies provided in the embodiments of this disclosure divides all microcystis colonies into multiple groups of microcystis colonies according to their size, obtains the quantity percentage of microcystis colonies in each group relative to the total quantity of all microcystis colonies, constructs an initial concentration control equation for each group of microcystis colonies, performs a simplification process on the initial concentration control equation based on a two-dimensional wind flow velocity pattern to obtain a target concentration control equation for each group of microcystis colonies, determines the initial concentration conditions for each group of microcystis colonies based on the quantity percentage and total quantity of microcystis colonies in each group, obtains the concentration boundary conditions for each group of microcystis colonies, and sets dimensionless parameters Based on the introduction, the target concentration control equation, initial concentration conditions, and concentration boundary conditions are transformed to obtain the dimensionless concentration control equation, initial concentration conditions, and boundary conditions for each group of microcystis colonies. Based on the target number of-order concentration moment analysis equation, target quantity-centered concentration moment analysis equation, pre-set concentration distribution addition conditions, dimensionless initial concentration conditions, and dimensionless concentration boundary conditions for each group of microcystis colonies, the dimensionless concentration control equation for each group of microcystis colonies is processed to obtain the target number of-order concentration moments and target quantity-centered concentration moments for each group of microcystis colonies. Based on the target number of-order concentration moments and target quantity-centered concentration moments for each group of microcystis colonies, the spatial concentration distribution analysis equation for all microcystis colonies is calculated.Regarding the transition distribution of multi-sized Microcystis colonies, the lack of theoretical analysis methods makes it difficult to accurately describe the transition distribution process of multi-sized Microcystis colonies, thus making it difficult to provide reliable support for predicting and forecasting water blooms. This technical problem is solved by the embodiment of this disclosure, which, unlike the theoretical analysis method for the transition distribution of single-sized Microcystis colonies, reflects the influence of different environmental gradients on the transition characteristics and transport features of Microcystis colony patches, enabling an accurate description of the transition distribution process of multi-sized Microcystis colonies and providing reliable support for predicting and forecasting water blooms. [Brief explanation of the drawing]

[0015] The above and other features, advantages and aspects of each embodiment of this disclosure will become clearer when viewed in conjunction with the drawings and the following specific embodiments. Throughout the drawings, the same or similar reference numerals represent the same or similar elements. Note that the drawings are schematic and the dimensions of actual objects and elements are not necessarily depicted proportionally.

[0016] [Figure 1] This is an illustrative flowchart of a theoretical analysis method for simulating the transition distribution of multi-sized microcystis colonies provided in the embodiments of this disclosure. [Figure 2] This is a schematic diagram illustrating the transition of microcystis colonies of different sizes in a wind-flow system, as provided in the embodiments of this disclosure. [Figure 3] This is a schematic diagram of the spatial concentration distribution of the microcystis colony provided in the embodiments of this disclosure. [Figure 4] This is a schematic diagram of the spatial concentration distribution of another microcystis colony provided in the embodiments of this disclosure. [Figure 5] This is a schematic diagram of the spatial concentration distribution of the microcystis colony provided in the embodiments of this disclosure. [Figure 6] This is a schematic diagram of the spatial concentration distribution of further microcystis colonies provided in the embodiments of this disclosure. [Figure 7]This is a schematic diagram of the structure of a theoretical analysis device for simulating the transition distribution of multi-sized microcystis colonies provided in the embodiments of this disclosure. [Modes for carrying out the invention]

[0017] The embodiments of this disclosure will be described in more detail below with reference to the drawings. Although some embodiments of this disclosure are shown in the drawings, this disclosure can be realized in various forms and should not be construed as being limited to the embodiments described herein. Rather, these embodiments are provided for a more thorough and complete understanding of this disclosure. The drawings and embodiments of this disclosure are for illustrative purposes only and do not limit the scope of protection of this disclosure.

[0018] It should be understood that the steps described in the embodiments of the method of this disclosure may be performed in a different order and / or in parallel. Furthermore, embodiments of the method may include additional steps and / or omission of the performed steps shown. The scope of this disclosure is not limited in this respect.

[0019] As used herein, the terms “including” and their variations mean “including, but not limited to.” The term “based on” means “at least partially based on….” The term “one embodiment” means “at least one embodiment,” the term “another embodiment” means “at least one other embodiment,” and the term “several embodiments” means “at least several embodiments.” Relevant definitions of other terms are given below.

[0020] Furthermore, the concepts of “First,” “Second,” etc., as used in this disclosure are used solely to distinguish between different devices, modules, or units, and are not intended to limit the order or interdependence of the functions performed by these devices, modules, or units.

[0021] Furthermore, the modifiers "one" and "multiple" as used in this disclosure are merely illustrative and not limiting; those concerned will understand that unless otherwise explicitly stated in the context, they should be understood as "one or more."

[0022] The names of messages or information exchanged between multiple devices in embodiments of this disclosure are for illustrative purposes only and are not used to limit the scope of such messages or information.

[0023] Conventional methods have a technical problem in that they cannot fully reflect the state of existence of a single-size microcystis colony in its natural environment. However, the theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies provided in the embodiments of this disclosure can overcome the aforementioned shortcomings. By constructing and analytically solving a multi-size microcystis colony concentration control equation that takes into account the different sizes of the microcystis colonies, the transition distribution status of multi-size microcystis colonies can be accurately reflected, thereby improving the prediction and forecasting accuracy of water blossoms.

[0024] Figure 1 is a schematic diagram of the flow of a theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies provided in the embodiments of this disclosure, which can be implemented in software and / or hardware and is generally integrable into electronic devices. As shown in Figure 1, the method comprises steps 101 to 106.

[0025] In step 101, all microcystis colonies are divided into multiple groups of microcystis colonies according to their size, and the quantity percentage of microcystis colonies in each group relative to the total quantity of all microcystis colonies is obtained.

[0026] Specifically, the transport process of microcystis colonies in water is very complex and involves convection-induced entrainment and turbulent diffusion of the water, longitudinal transition and self-diffusion of microcystis colonies, and interactions between microcystis colonies and water. In the embodiments of this disclosure, the convection-diffusion effect of the water and the transition characteristics of the microcystis colonies themselves are mainly considered, while the effects of changes in the density of the microcystis colonies, interactions between microcystis colonies, and the hydrodynamic transition characteristics of the microcystis colonies themselves are ignored.

[0027] In the embodiments of this disclosure, each microcystis colony has various different shapes in the actual application scene, and pretreatment can simplify the various different shapes of microcystis colonies into circles. Furthermore, it can be understood that all microcystis colonies contain microcystis colonies of different sizes, where size can be understood as the size of the equivalent diameter, and microcystis colonies of different sizes can be understood as circles with different equivalent diameters, and such circles contain multiple microcystis.

[0028] In the embodiments of this disclosure, all microcystis colonies can be divided into multiple groups of microcystis colonies according to their size. For example, different size ranges such as 15-20 micrometers and 20-25 micrometers can be set, and the microcystis colonies belonging to different size ranges can be divided into multiple groups, with each group containing microcystis colonies of the same or similar size.

[0029] Furthermore, by statistically analyzing the quantity of each microcystis within each group of microcystis colonies, the total quantity of all microcystis colonies can be obtained, thereby calculating the percentage of microcystis colonies in each group relative to the total quantity of all microcystis colonies. Subsequently, multiple groups of microcystis colonies can be processed simultaneously to quickly and effectively obtain theoretical analysis of the transition distribution of microcystis colonies of various sizes.

[0030] Specifically, in natural environments, there are generally microcystis colonies of different sizes, and all microcystis colonies are divided into n groups according to their size, where the equivalent diameter is D. i The proportion of the i-th group of microcystis colonies to the total colonies is N i The concentration C of multi-size microcystis colonies can be calculated from the following formula (1): JPEG0007923361000014.jpg29151 Here, C i This represents the concentration of microcystis colonies in group i.

[0031] As an example, consider the transition process of microcystis colony patches in a wind-flowing water system in Cartesian coordinates, where x is parallel to the water flow direction, the z axis points vertically upward, and the origin is set at the bottom. Figure 2 shows a schematic diagram illustrating the transition of microcystis colonies of different sizes in a wind-flowing water system.

[0032] For example, as shown in Figure 2, every microcystis colony contains three different sizes of microcystis colonies, so they can be divided into three groups, and subsequent treatments can be performed simultaneously on the microcystis colonies of all three groups.

[0033] In step 102, the initial concentration control equations for each group of microcystis colonies are constructed, and the initial concentration control equations are simplified based on the two-dimensional wind flow velocity pattern to obtain the target concentration control equations for each group of microcystis colonies.

[0034] In several examples, the initial concentration control equations for the microcystis colonies of each group were constructed and are shown in equation (2). JPEG0007923361000015.jpg29157

[0035] Based on the two-dimensional wind flow velocity pattern, the initial concentration control equation was simplified to obtain the target concentration control equation for each group of microcystis colonies, which is shown in equation (3). JPEG0007923361000017.jpg23170

[0036] JPEG0007923361000018.jpg45170

[0037] In step 103, the initial concentration conditions for microcystis colonies in each group are determined based on the quantity percentage and total quantity of microcystis colonies in each group, and the concentration boundary conditions for microcystis colonies in each group are obtained.

[0038] In some examples, an initial time is set, the microcystis colony patches are concentrated at x=0 and z=0, and the initial concentration conditions for the microcystis colonies in each group are determined based on the quantity percentage and total quantity of microcystis colonies in each group, as shown in equation (4). JPEG0007923361000019.jpg19170

[0039] JPEG0007923361000020.jpg54170

[0040] JPEG0007923361000021.jpg33164 Here, equation (6) is used to show that the flux of the microcystis colony at the free surface and bottom is 0.

[0041] In step 104, based on the introduction of dimensionless parameters, the target concentration control equation, initial concentration conditions, and concentration boundary conditions are transformed to obtain the dimensionless concentration control equation, initial concentration conditions, and boundary conditions for each group of microcystis colonies.

[0042] In some embodiments, the dimensionless parameter to be introduced is: JPEG0007923361000022.jpg15170

[0043] JPEG0007923361000023.jpg35170

[0044] Based on the introduction of dimensionless parameters, the target concentration control equation, initial concentration conditions, and concentration boundary conditions are transformed to obtain the dimensionless concentration control equation, initial concentration conditions, and boundary conditions for each group of microcystis colonies, which are as follows. JPEG0007923361000024.jpg87170

[0045] In other words, by substituting the dimensionless parameters into equations (3), (4), (5), and (6), we obtain equations (7), (8), (9), and (10).

[0046] JPEG0007923361000025.jpg54170

[0047] In step 105, the dimensionless concentration control equation for each group of microcystis colonies is processed based on the target number of-order concentration moment analysis equation, the target quantity center concentration moment analysis equation, pre-set concentration distribution addition conditions, dimensionless initial concentration conditions, and dimensionless concentration boundary conditions for each group of microcystis colonies to obtain the target number of-order concentration moments and the target quantity center concentration moments for each group of microcystis colonies.

[0048] In the embodiments of this disclosure, the target quantity can be set as needed, for example, four floors.

[0049] In several examples, the Aris concentration moment method was developed to analyze the overall transition characteristics of microcystis colonies in windy flow, and the p-th order concentration moment and p-th order central concentration moment of the i-th group of microcystis colonies were defined. The target number order concentration moment analysis formulas for each group of microcystis colonies are as follows. JPEG0007923361000026.jpg29170

[0050] The target quantity central concentration moment analysis formula for each group of microcystis colonies is as follows: JPEG0007923361000027.jpg48170

[0051] JPEG0007923361000028.jpg44170 shows the position of the centroid of a Microcystis colony patch.

[0052] For example, the obtained second-order central concentration moment, third-order central concentration moment, and fourth-order central concentration moment are, respectively, The filename is JPEG0007923361000029.jpg36170.

[0053] In some examples, the pre-set conditions for adding concentration distribution are as follows. JPEG0007923361000030.jpg28170

[0054] Based on the target order-order concentration moment analysis equation for each group of microcystis colonies, the pre-set concentration distribution addition conditions, the dimensionless initial concentration conditions, and the dimensionless concentration boundary conditions, the dimensionless concentration control equation for each group of microcystis colonies is processed. JPEG0007923361000031.jpg181170

[0055] In step 106, the spatial concentration distribution analysis formula for all microcystis colonies is obtained by calculating based on the target number-order concentration moment and the target number-center concentration moment of each group of microcystis colonies.

[0056] In some embodiments, the step of obtaining a spatial concentration distribution analysis formula for all microcystis colonies by calculating based on the target number of order concentration moments and the target number of central concentration moments for each group of microcystis colonies is: The process includes the step of applying the Chatwin asymptotic expansion equation based on the target order of concentration moments of each group of microcystis colonies to obtain the following spatial concentration distribution analysis equation for all microcystis colonies: JPEG0007923361000032.jpg63170 Here, the polynomial coefficient a depends on the target quantity and central concentration moment. n To decide.

[0057] It should be understood that the target quantity can be set as needed. In this disclosure, using the 4th order as an example, and based on the concentration moments from the 0th to the 4th order calculated previously, the Chatwin asymptotic expansion, i.e., the Edgeworth expansion, is applied to obtain the spatial concentration distribution analysis formula for a microcystis colony in a wind-flow system. Therefore, the polynomial coefficients are: The filename is JPEG0007923361000033.jpg19170.

[0058] Therefore, by obtaining the first four coefficients and using them as an approximate solution, computational efficiency is guaranteed while ensuring computational efficiency.

[0059] This disclosure establishes a two-dimensional model of multi-size microcystis colony transport in a shallow lake environment, taking into account the transition characteristics of the microcystis colonies themselves, based on a multi-size microcystis colony concentration distribution model. By obtaining analytical solutions of concentration moments from the 0th to the 4th level using the Aris concentration moment method, the overall transition characteristics of the microcystis colonies are obtained. Different concentration moments can reflect different characteristics of the transition of the microcystis colonies. Based on the analytical solutions of concentration moments from the 0th to the 4th level, a two-dimensional concentration distribution of the microcystis colonies is obtained by combining it with the Chatwin asymptotic expansion. By simultaneously calculating the theoretical analysis of the transition distribution of multiple groups of microcystis colonies of different sizes, the theoretical analysis of the transition distribution of all microcystis colonies of different sizes is obtained, enabling rapid and accurate simulation of the transition distribution process of multi-size microcystis colonies and providing reliable support for water flower prediction and forecasting.

[0060] To summarize the above, the theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies provided in the embodiments of this disclosure divides all microcystis colonies into multiple groups of microcystis colonies according to their size, obtains the quantity percentage of microcystis colonies in each group relative to the total quantity of all microcystis colonies, constructs an initial concentration control equation for each group of microcystis colonies, performs a simplification process on the initial concentration control equation based on a two-dimensional wind flow velocity pattern to obtain a target concentration control equation for each group of microcystis colonies, determines the initial concentration conditions for each group of microcystis colonies based on the quantity percentage and total quantity of microcystis colonies in each group, obtains the concentration boundary conditions for each group of microcystis colonies, and performs a dimensionless process Based on the introduction of the lameter, the target concentration control equation, initial concentration conditions, and concentration boundary conditions are transformed to obtain the dimensionless concentration control equation, initial concentration conditions, and boundary conditions for each group of microcystis colonies. Based on the target number of-order concentration moment analysis equation, target quantity-center concentration moment analysis equation, pre-set concentration distribution addition conditions, dimensionless initial concentration conditions, and dimensionless concentration boundary conditions for each group of microcystis colonies, the dimensionless concentration control equation for each group of microcystis colonies is processed to obtain the target number of-order concentration moments and target quantity-center concentration moments for each group of microcystis colonies. Based on the target number of-order concentration moments and target quantity-center concentration moments for each group of microcystis colonies, the spatial concentration distribution analysis equations for all microcystis colonies are calculated. Regarding the transition distribution of multi-sized Microcystis colonies, the lack of theoretical analysis methods makes it difficult to accurately describe the transition distribution process of multi-sized Microcystis colonies, thus making it difficult to provide reliable support for predicting and forecasting water blooms. This technical problem is solved by the embodiment of this disclosure, which, unlike the theoretical analysis method for the transition distribution of single-sized Microcystis colonies, reflects the influence of different environmental gradients on the transition characteristics and transport features of Microcystis colony patches, enabling an accurate description of the transition distribution process of multi-sized Microcystis colonies and providing reliable support for predicting and forecasting water blooms.

[0061] As an example, we study the spatial concentration distribution of microcystis colony patches under different environmental gradient conditions based on a transport model of multi-size microcystis colonies, and further interpret the differences in microcystis colony concentration distribution under different environmental gradient conditions. Tables 1-3 show the proportion and ascent velocity of microcystis colonies of different sizes under three environmental gradient conditions.

[0062] Table 1 The proportion of microcystis colonies of different sizes and their ascent rates in an environmental gradient of 1. JPEG0007923361000034.jpg75170

[0063] Table 2 The proportion of microcystis colonies of different sizes and their ascent rates in environmental gradient 2. JPEG0007923361000035.jpg77170

[0064] Table 3 The proportion of microcystis colonies of different sizes and their ascent rates in an environmental gradient of 3. JPEG0007923361000036.jpg65170

[0065] From this, we can see that the concentration distribution of microcystis colonies differs under different environmental gradient conditions. Specifically, as shown in Figures 3 to 6, these are the spatial concentration distributions of microcystis colonies at different times (τ = 0.1, 0.5, 1.0, 5.0) under different environmental gradient conditions. Figure 3 shows the changes in Ω with ξ and ζ, and Pe x =1.0, τ=0.1, and in Figure 3, (a) is environmental gradient 1 from Table 1, (b) is environmental gradient 2 from Table 2, and (c) is environmental gradient 3 from Table 3. Figure 4 shows the changes of Ω with respect to ξ and ζ, and Pe x =1.0, τ=0.5, and in Figure 4, (a) is environmental gradient 1 from Table 1, (b) is environmental gradient 2 from Table 2, and (c) is environmental gradient 3 from Table 3. Figure 5 shows the changes of Ω with respect to ξ and ζ, and Pe x=1.0, τ=1.0, and in Figure 5, (a) is environmental gradient 1 from Table 1, (b) is environmental gradient 2 from Table 2, and (c) is environmental gradient 3 from Table 3. Figure 6 shows the changes of Ω with respect to ξ and ζ, and Pe x =1.0 and τ=5.0. In Figure 6, (a) is environmental gradient 1 from Table 1, (b) is environmental gradient 2 from Table 2, and (c) is environmental gradient 3 from Table 3.

[0066] Figure 7 is a schematic diagram of the structure of a theoretical analysis device for simulating the transition distribution of multi-size microcystis colonies provided in the embodiments of this disclosure, the device can be implemented in software and / or hardware and is generally integrable into electronic equipment. As shown in Figure 7, the device is A grouping and acquisition module 701 for dividing all microcystis colonies into multiple groups of microcystis colonies according to their size, and obtaining the quantity percentage of microcystis colonies in each group relative to the total quantity of all microcystis colonies, A construction / simplification module 702 is used to construct initial concentration control equations for each group of microcystis colonies, and to perform a simplification process on these initial concentration control equations based on a two-dimensional wind flow velocity pattern to obtain target concentration control equations for each group of microcystis colonies. Based on the quantity percentage of microcystis colonies in each group and the total quantity, a determination / acquisition module 703 is used to determine the initial concentration conditions of the microcystis colonies in each group and to obtain the concentration boundary conditions of the microcystis colonies in each group. A transformation module 704 for transforming the target concentration control equation, the initial concentration conditions, and the concentration boundary conditions based on the introduction of dimensionless parameters to obtain the dimensionless concentration control equation, the initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis colonies, A processing module 705 processes the dimensionless concentration control equation for each group of microcystis groups based on the target number of-order concentration moment analysis equation, the target quantity center concentration moment analysis equation, the preset concentration distribution addition conditions, the dimensionless initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis groups to obtain the target number of-order concentration moments and the target quantity center concentration moments for each group of microcystis groups. The system includes a calculation module 706 for obtaining spatial concentration distribution analysis formulas for all microcystis colonies by calculating based on the target number-order concentration moment and target number-center concentration moment of each group of microcystis colonies.

[0067] The theoretical analysis device for simulating the transition distribution of multi-sized microcystis colonies provided in the embodiments of this disclosure can perform a theoretical analysis method for simulating the transition distribution of multi-sized microcystis colonies provided in any embodiment of this disclosure, and includes a corresponding functional module for performing the method, and has beneficial effects associated with the method.

[0068] The embodiments of this disclosure provide a computer program product including a computer program / instruction, which, when executed by a processor, performs a theoretical analysis method for simulating the transition distribution of multi-sized microcystis colonies provided in any embodiment of this disclosure.

[0069] According to one or more embodiments of this disclosure, this disclosure provides electronic equipment, Processor and The processor includes a memory for storing executable instructions, The processor is used to read the executable instructions from the memory and execute the instructions in order to carry out a theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies as described in any one of the items provided in this disclosure.

[0070] According to one or more embodiments of the present disclosure, the present disclosure provides a computer-readable storage medium that stores a computer program for performing a theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies as described in any one of the items provided in the present disclosure.

[0071] The above description is merely an explanation of preferred embodiments and applied technical principles of the present disclosure. Those skilled in the art will understand that the scope of the present disclosure is not limited to technical solutions formed by specific combinations of the above technical features, but should include other technical solutions formed by arbitrary combinations of the above technical features or equivalent features, without departing from the above-disclosed concepts. For example, technical solutions formed by replacing the above features with technical features having similar functions disclosed in this disclosure (but not limited to these).

[0072] Furthermore, although each operation is described in a specific order, it should not be understood that these operations must be performed in a specific order or sequence as indicated. In certain environments, multitasking and parallel processing may be advantageous. Similarly, although the above description includes some specific implementation details, these should not be construed as limiting the scope of this disclosure. Some features described in the context of individual embodiments may be implemented in a single embodiment in a combined manner. Conversely, various features described in the context of a single embodiment may be implemented in multiple embodiments individually or in any appropriate sub-combination.

[0073] While this subject matter has been described using language specific to structural features and / or methodological operations, it should be understood that the subject matter as defined in the attached claims is not necessarily limited to the specific features or operations described above. Rather, the specific features and operations described above are merely exemplary forms of implementing the claims.

Claims

1. A theoretical analysis method for simulating the transition distribution of multi-size microcystis colonies, The above method is performed using a theoretical analysis system that simulates the transition distribution of multi-sized microcystis colonies. The processor divides all microcystis colonies into multiple groups of microcystis colonies according to their size, and obtains the quantity percentage of each group of microcystis colonies relative to the total quantity of all microcystis colonies. The processor constructs an initial concentration control equation for each group of microcystis colonies, and performs a simplification process on the initial concentration control equation based on a two-dimensional wind flow velocity pattern to obtain a target concentration control equation for each group of microcystis colonies. Here, the initial concentration control equation is: And, C ｉ is the concentration of the i-th group of microcystis colonies, where i is a positive integer, t is time, and μ is the vector velocity. b i λ is the vector transition velocity of the i-th group of microcystis colonies, λ is the diffusion coefficient of the water flow, and ∇ is the gradient operator. The aforementioned target concentration control equation is: And, μ is the longitudinal flow velocity, w b i is the transition velocity of the i-th group of Microcystis colonies themselves, λ b ixx is the diffusion coefficient along the x-direction of the i-th group of Microcystis colonies, λ b izz is the diffusion coefficient along the z-direction of the i-th group of Microcystis colonies, After setting the initial time, the microcystis colony patches are concentrated at x=0 and z=0. The processor determines the initial concentration conditions for the microcystis colonies in each group based on the quantity percentage and total quantity of the microcystis colonies in each group, and obtains the concentration boundary conditions for the microcystis colonies in each group. The initial conditions for the concentration of the aforementioned microcystis colony are: And, N ｉ is the quantity percentage of the i-th group of microcystis colonies, Q is the total quantity of microcystis colonies, and δ is the Dirac function. At x=±∞, z=0, z=H (where H is the water depth), the concentration boundary conditions for the microcystis colonies of each group are: This is used to indicate that the flux of the microcystis colony at the free surface and bottom is 0. The processor transforms the target concentration control equation, the initial concentration conditions, and the concentration boundary conditions based on the introduction of dimensionless parameters to obtain the dimensionless concentration control equation, the initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis colonies. The dimensionless parameters in the first half are: And, Based on the introduction of the dimensionless parameter, the target concentration control equation, the initial concentration conditions, and the concentration boundary conditions are transformed to obtain the dimensionless concentration control equation, the initial concentration conditions, and the dimensionless concentration boundary conditions for each group of microcystis colonies. The dimensionless concentration control equation, the dimensionless concentration initial condition, and the dimensionless concentration boundary condition are, They are as follows: Reflecting the relative strength of longitudinal transition and total longitudinal diffusion in the i-th group of microcystis colonies, The processor processes the dimensionless concentration control equation for each group of microcystis groups based on the target number of order concentration moment analysis equations, the target quantity center concentration moment analysis equations, preset concentration distribution addition conditions, the dimensionless initial concentration conditions, and the dimensionless concentration boundary conditions, in order to obtain the target number of order concentration moments and the target quantity center concentration moments for each group of microcystis groups. Here, the target number-order concentration moment analysis formula is: And, The aforementioned target quantity central concentration moment analysis formula is: And, Here, And, Here, the pre-set concentration distribution conditions are: And, Based on the target order-order concentration moment analysis formula for each group of microcystis colonies, the pre-set concentration distribution addition conditions, the dimensionless initial concentration conditions, and the dimensionless concentration boundary conditions, the dimensionless concentration control equation for each group of microcystis colonies is processed. If p = 0, If p > 0, And, The process includes the step of the processor calculating, based on the target number of order concentration moments and the target number of central concentration moments of each group of microcystis colonies, to obtain a spatial concentration distribution analysis formula for all microcystis colonies, Based on the target order-order concentration moments of the microcystis colonies in each of the aforementioned groups, applying the Chatwin asymptotic expansion equation yields the equation for analyzing the spatial concentration distribution of all microcystis colonies: Here, the polynomial coefficient a depends on the target quantity central concentration moment. ｎ to decide A theoretical analysis method for simulating the transition distribution of multi-sized microcystis colonies, characterized by the following features.

2. Processor and The processor includes a memory for storing executable instructions, The processor is used to read the executable instructions from the memory and execute the instructions in order to carry out the theoretical analysis method for simulating the transition distribution of multi-size microcystis groups described in claim 1 above. An electronic device characterized by the following features.

3. A computer-readable storage medium characterized by storing a computer program for executing a theoretical analysis method for simulating the transition distribution of multi-sized microcystis colonies as described in claim 1 above.

Images