Dynamic simulation method for the structural stability of aquatic ecosystem networks

Abstract

This invention provides a dynamic simulation method for the structural stability of aquatic ecosystem networks, which improves the accuracy of simulation analysis results by dynamically simulating changes in the aquatic environment. [Solution] The dynamic simulation method constructs an aquatic ecosystem network based on ecological information in the aquatic environment, removes nodes in the aquatic ecosystem network according to an importance index of the aquatic ecosystem network, dynamically adjusts the structure of the aquatic ecosystem network to simulate the temporal changes in the aquatic environment, determines the destruction resistance threshold of the aquatic ecosystem network based on this, and uses it as the result of a simulated analysis of the structural stability of the aquatic ecosystem network.

Description

[Technical Field]

[0001] The present invention relates to the field of ecological environments, and more specifically, to the structural stability of aquatic ecosystem networks. This concerns dynamic simulation methods. [Background technology]

[0002] Currently, in the field of environmental science, assessing the stability of aquatic ecosystem networks in aquatic environments is important. This is the direction of research. Ecological networks describe the relationships between biological populations in aquatic environments. It is a mathematical model established based on ecological principles, and it describes the functional structure and properties of aquatic ecosystems. It can effectively reflect quality cycles. Conventional methods for evaluating ecological network structures are mainly We measure the stability of a system through the number of biological populations, species composition, and interaction relationships. This method depends on the physical distribution of species and population dynamics, and encompasses the ecosystem network structure. A thorough analysis is lacking, and the spatial distribution of aquatic organisms in the water environment, pollutants, and environmental factors are not being analyzed for aquatic ecosystems. The comprehensive impact on [the subject] is often overlooked, resulting in insufficient accuracy in the analysis results. [Overview of the project]

[0003] The present invention is as follows: In a first embodiment, the present invention provides a method for dynamically simulating the structural stability of aquatic ecosystem networks. This method includes the following: Based on ecological information in aquatic environments, we will construct aquatic ecosystem networks. Here, ecological information includes species information, environmental factors, and pollutant information, and species information is present in the aquatic environment. This includes the spatial distribution of multiple species of aquatic organisms, including algae, fish, and aquatic plants. Environmental factors include temperature, pH value, and dissolved oxygen of the aquatic environment, and pollutant information is water This includes the concentrations of multiple types of pollutants present in the local environment, which include heavy metals and organic pollutants. And containing nutrients, the aquatic ecosystem network connects multiple nodes and multiple nodes Including multiple consecutive edges, Based on pre-set removal ratios and importance indicators for each node in the aquatic ecosystem network and at least one node and at least one node in the aquatic ecosystem network Remove at least one edge connected to the dot, and after each removal, the aquatic ecosystem network Calculate the stability index value, After removal, until the stability index value of the aquatic ecosystem network falls below the stability threshold, the aquatic ecosystem Number of nodes removed from the system network and the number of nodes in the aquatic ecosystem network before repeated removals Based on the total number of units, the threshold for resistance to destruction of aquatic ecosystem networks is determined, and the aquatic ecosystem network The fracture resistance threshold of the workpiece is used as a result of a simulated analysis of the structural stability of the aquatic ecosystem network. ru. In one embodiment of the first aspect, multiple nodes of the aquatic ecosystem network are multiple species of aquatic It has a one-to-one correspondence with organisms, environmental factors, and multiple types of pollutants, and is one of the aquatic ecosystem networks. The edges represent the ecological interaction relationships between two nodes connected by one edge. The relationships include predation, symbiosis, competition, survival threat, and mutual influence, where ecological relationships are defined as water Spatial distribution of multiple aquatic organisms present in the local environment, environmental factors, and multiple organisms included in the aquatic environment It is determined based on the concentration of contaminants in the species. In one embodiment of the first aspect, a predetermined removal ratio and in the aquatic ecosystem network Based on the importance indicator of each node, at least one node in the aquatic ecosystem network Remove at least one edge connected to at least one node, after each removal. The step of calculating the stability index value of the aquatic ecosystem network includes: calculating an importance index of each node in the aquatic ecosystem network, wherein the importance index satisfies the following formula: JPEG2026143293000027.jpg21122wherein, I i represents the importance index of the i-th node in the aquatic ecosystem network, k i represents the node degree of the i-th node in the aquatic ecosystem network, B i represents the node centrality of the i-th node in the aquatic ecosystem network, C i represents the node correlation degree of the i-th node in the aquatic ecosystem network, α1 represents the weight of the node degree, α 2 represents the weight of the node centrality, and α3 represents the weight of the node correlation degree; sorting the removal order of nodes in the aquatic ecosystem network according to the magnitude of the importance index, and iteratively removing at least one node matching a preset removal ratio and at least one edge connected to the at least one node from the aquatic ecosystem network, after each removal operation, calculating a stability index value of the aquatic ecosystem network after each removal, wherein the stability index value of the aquatic ecosystem network after each removal is the network efficiency of the aquatic ecosystem network after each removal and the maximum connectivity of the aquatic ecosystem network after each removal, wherein the preset removal ratio is preset as the ratio of the number of nodes removed in each removal to the total number of nodes of the aquatic ecosystem network before iterative removal, and the preset removal ratio is 1%, 2%, 3% or 5% ; the network efficiency of the aquatic ecosystem network after each removal satisfies the following formula: ​JPEG2026143293000028.jpg39122where E represents the network efficiency of the water ecosystem network after each removal, and N is the number of iterations resents the total number of nodes in the aquatic ecosystem network before repeated removal, d ji represents the aquatic ecosystem after each removal represents the shortest path length between node i and node j in the system network, the maximum connectivity of the aquatic ecosystem network after each removal satisfies the following formula: JPEG2026143293000029.jpg3873where S represents the maximum connectivity of the aquatic ecosystem network after each removal, |C max | represents the total number of nodes in the maximum connectivity subgraph of the aquatic ecosystem network after each removal. In an embodiment of the first aspect, the stability threshold comprises a network efficiency threshold and a maximum connectivity threshold the value range of the network efficiency threshold is 0.10 to 0.15, and the maximum connectivity threshold has a value range of 0.4 to 0.6. In an embodiment of the first aspect, the stability index value of the aquatic ecosystem network after removal is the stability the number of removed nodes of the aquatic ecosystem network and the number before repeated removal until the value is below the threshold based on the total number of nodes of the aquatic ecosystem network, the destruction resistance threshold of the aquatic ecosystem network is determined the step of determining the value comprises: the network efficiency of the aquatic ecosystem network after removal is less than or equal to the network efficiency threshold , or when the maximum connectivity of the aquatic ecosystem network after removal is less than or equal to the maximum connectivity threshold , the number of removed nodes of the aquatic ecosystem network and the aquatic ecosystem network before repeated removal the destruction resistance threshold of the aquatic ecosystem network is calculated from the total number of nodes of the work, and the aquatic ecosystem the destruction resistance threshold of the network satisfies the following formula: JPEG2026143293000030.jpg3053 Here, fa represents the fracture resistance threshold, and Nr represents the removed noise of the aquatic ecosystem network. This represents the number of nodes, where N represents the total number of nodes in the aquatic ecosystem network before repeated removal. In one embodiment of the first aspect, the method further includes the following steps: When the stability index value of the aquatic ecosystem network after removal falls below the stability threshold, the aquatic ecosystem The resilience of the aquatic ecosystem network is determined based on the number of nodes removed from the network. By combining the resilience of aquatic ecosystem networks and the threshold for their resistance to disruption, water As a result of a simulated analysis of the structural stability of regional ecosystem networks, the resilience of aquatic ecosystem networks is It is used to describe the resilience and recovery capacity of aquatic ecosystem networks to changes in the ecological environment. The resilience of the aquatic ecosystem network is satisfied by the following equation: JPEG2026143293000031.jpg33110 Here, R represents the resilience of the aquatic ecosystem network, and Nr represents the aquatic ecosystem network This represents the number of removed nodes, d ji This refers to node i in the aquatic ecosystem network after removal. This represents the shortest path length between node j and node j. [Effects of the Invention]

[0004] The present invention has the following beneficial effects. This invention constructs an aquatic ecosystem network based on ecological information in the aquatic environment, and pre-configures Based on the determined removal ratio and the importance indicator of each node in the aquatic ecosystem network, repeated removal Through the above method, at least one node in the aquatic ecosystem network and at Remove at least one edge connected to a node several times, and after each removal, the water body after removal Calculate the stability index value of the ecosystem network. Stability of the aquatic ecosystem network after removal. Until the indicator value falls below the stability threshold, the ecology of the aquatic environment corresponding to the aquatic ecosystem network The overall function of the system is impaired and it cannot maintain normal operation, i.e., the aquatic ecosystem network It is determined that the network becomes unstable. The number of nodes removed and the water volume before repeated removal are also determined. The threshold for the destruction resistance of the aquatic ecosystem network is calculated from the total number of nodes in the ecosystem network. The fracture resistance threshold is used as the result of a simulated analysis of the structural stability of the aquatic ecosystem network. The above method is: During the process of constructing aquatic ecosystem networks, species information, environmental factors, and pollutants in the aquatic environment Considering the information, and based on this, at least one node in the aquatic ecosystem network By removing at least one edge connected to at least one node several times , the changes in the aquatic environment caused by natural or human activities in actual situations ( It can dynamically simulate changes in species, environmental factors, and pollutants, and therefore This allows for more accurate simulations of the stability of aquatic ecosystem network structures. [Brief explanation of the drawing]

[0005] [Figure 1] This is a schematic diagram (part 1) of the dynamic simulation method for the structural stability of aquatic ecosystem networks provided by the examples. [Figure 2] This is a schematic diagram (part 2) of the dynamic simulation method for the structural stability of aquatic ecosystem networks provided by the examples. [Figure 3] This is a schematic diagram of the node removal process in the aquatic ecosystem network in the example. [Figure 4] This is a schematic diagram of indicators related to the aquatic ecosystem network after node removal in the example, where (a) shows connectivity, (b) shows the relative size of the maximum connectivity subgraph, (c) shows network efficiency, and (d) shows the ratio of maximum connectivity to node attacks. [Modes for carrying out the invention]

[0006] The methods provided by embodiments of this application relate to aquatic ecosystem networks, and aquatic ecosystem By removing nodes from the network, the water area corresponding to the aquatic ecosystem network It can dynamically simulate changes in the environment over time and the disruption of fluctuating aquatic ecosystem networks. By using the collapse resistance threshold as the result of a simulated analysis of the stability of the aquatic ecosystem network structure, It can reflect the stability of the environment. Naturally, an aquatic ecosystem network is defined as a network of a specific aquatic environment at a spatial scale. Identify ecological factors in aquatic bodies (including environmental factors, pollutants, and aquatic organisms), and identify different ecological factors It is a complex system that connects offspring and reflects the interactions and correlations between each ecological element. The stability of aquatic ecosystem networks is affected by natural or human activities in waters. When ecological changes occur in ecosystems corresponding to ecosystem networks, the normal functioning of these networks is maintained. It reflects the sustainable capacity it possesses. Studying the resilience of ecological network structures is important for ecosystems. In ensuring stability and sustainability, external disturbances such as environmental changes, species invasion, or pollution are particularly important. Even when faced with turmoil, it is extremely important for maintaining the normal functioning of ecosystems. The spatial distribution of aquatic organisms refers to the distribution patterns and locations of aquatic organisms within aquatic environments. Reflecting the aggregation, dispersion, or movement patterns of aquatic organisms at a given spatial scale, The spatial distribution of is typically influenced by environmental factors, biointeractions, and pollutants. To solve the aforementioned problems in the background technology, the embodiment describes the structure of aquatic ecosystem networks. It provides a dynamic simulation method for stability, and during the process of constructing aquatic ecosystem networks, in the aquatic environment Considering species information, environmental factors, and pollutant information, and based on these, the aquatic ecosystem network At least one node in the network and at least one node connected to it By removing one side several times, the result is drawn by natural or human action in a real-world situation. The changes in the aquatic environment caused by these changes (including species changes, environmental factor changes, and pollutant changes) It can simulate the structure of aquatic ecosystem networks more accurately, and the results of the simulated analysis of the stability of aquatic ecosystem networks can be more precise. It can be done. As shown in Figure 1, the dynamics of the structural stability of the aquatic ecosystem network provided by the embodiment The target simulation method includes S101 to S103. S101: Construct an aquatic ecosystem network based on ecological information in the aquatic environment. Here, ecological information includes species information, environmental factors, and pollutant information. Species information is present in the aquatic environment. This includes the spatial distribution of multiple species of aquatic organisms, including algae, fish, and aquatic plants. Includes substances. Environmental factors include temperature, pH value, and dissolved oxygen of the aquatic environment. Pollutant information is: This includes concentrations of multiple types of pollutants present in the aquatic environment, and multiple types of aquatic organisms, including algae, fish, and The contaminants, including aquatic plants, contain multiple types of heavy metals, organic pollutants, and nutrients. Selectively, the above multiple types of aquatic organisms may further include microorganisms, and the above environmental factors may be hardness The above-mentioned multiple types of contaminants may further include other contaminants, and this application In the embodiment, the above-mentioned multiple types of aquatic organisms, the above-mentioned environmental factors and the above-mentioned multiple types of pollutants are further This is not limited to. Specifically, in order to acquire the above ecological information, remote sensing Through technology, field sampling, and other means, we can gather species information and data on the airspace of multiple aquatic organisms in aquatic environments. Inter-species distribution data and species density may be acquired, and environmental factors can be continuously monitored using water quality monitoring equipment. Alternatively, sampling may be performed periodically, and multiple types can be identified using water quality monitoring equipment and laboratory analysis. You may obtain the concentration of contaminants. Regarding the above aquatic ecosystem network, the above aquatic ecosystem network consists of multiple nodes It also includes multiple edges that connect multiple nodes. Multiple nodes in an aquatic ecosystem network It has a one-to-one correspondence with multiple species of aquatic organisms, environmental factors, and multiple species of pollutants. Aquatic ecosystem One edge of a network represents the ecological relationship between two nodes connected by that edge. This represents ecological relationships, including predation, symbiosis, competition, survival threats, and mutual influence. Here, Ecological relationships involve the spatial distribution of multiple aquatic organism species present in aquatic environments, environmental factors, and the aquatic environment itself. It is determined based on the concentrations of multiple types of pollutants present in the environment. For example, for a single edge, there are six possible ecological relationships that it represents. ru. Case 1: If two nodes connected by one edge are both aquatic organisms, The ecological relationship represented by this side is one of the relationships between predator, symbiosis, and competition. The interaction is determined based on the spatial distribution of multiple aquatic organism species present in the aquatic environment. Case 2: One of two nodes connected by one edge is an aquatic organism. If the other node is an environmental factor, the ecological relationship represented by that node is a survival crisis, and above The environmental factors described above mean that they can lead to the extinction of the aquatic organisms. The assignment is determined based on environmental factors. Case 3: One of two nodes connected by one edge is an aquatic organism. If the other node is polluted, the ecological relationship represented by that node is a survival crisis, as described above. This means that the pollutants will cause the extinction of the above-mentioned aquatic organisms. The ecological relationship that this represents is This is determined based on the concentration of the above-mentioned pollutants in multiple types of pollutants contained in the aquatic environment. Case 4: If two nodes connected by one edge are both environmental factors, The ecological relationship represented by that side is mutual influence. The ecological relationship represented by that side is based on environmental factors. It will be decided. Case 5: If two nodes connected by one edge are both contaminants, The ecological relationships represented by the edges are mutual influences. These ecological relationships are included in the aquatic environment. It is determined based on the concentrations of the two aforementioned contaminants in multiple types of contaminants. Case 6: One of two nodes connected by one edge is an environmental factor. If the other node is a pollutant, the ecological relationship represented by that node is mutual influence. Furthermore, in the above aquatic ecosystem network, each node further includes node characteristics, A key feature is that it can be used as supplementary information during the simulation process of aquatic ecosystem networks. Please note the following. Specifically, if a node is an aquatic organism, the nodal features of this node are The node may be the species name, spatial distribution, and species density of an aquatic organism, and the node may be an environmental factor. The node features of the node may be the values ​​of the environmental factor from continuous or periodic sampling. Often, if a node is contaminated, the node characteristics of that node are continuous or periodic sampling. The concentration of the contaminant may be the same as the concentration of the contaminant in the rind. In one application scenario, the implementation uses graph theory to construct an aquatic ecosystem network. The resulting aquatic ecosystem network structure is shown in Table 1.

[0007] Table 1: Ecological Network Configuration Information Table JPEG2026143293000032.jpg238168 In the above application scenario, after obtaining the aquatic ecosystem network, The adjacency matrix is ​​constructed using the connection weights w between the nodes, and the Floyd-Washer algorithm is used. The shortest path matrix W between each node is determined using the Zoom. It is shown as follows. JPEG2026143293000033.jpg2659 Here, w ij This is the path length between nodes i and j. S102, Pre-set removal ratios and importance indicators for each node in the aquatic ecosystem network Based on this, at least one node in the aquatic ecosystem network and at least one Iteratively remove at least one edge connected to the node, and after each removal, the aquatic ecosystem network Calculate the stability index value of the workpiece. Repetition refers to the action of repeating a feedback process, where each repetition of the process is counted as one repetition. It is called a repeat, and it is understood that the result obtained in one iteration is used as the initial value for the next iteration. In embodiments of this application, repeated removal of nodes and edges in an aquatic ecosystem network is performed in the aquatic ecosystem. This is an operation that repeatedly removes nodes and edges in an ecosystem network, and after each removal operation, The aquatic ecosystem network that was previously identified will be used as the target for repeated removal in the next cycle. Selectable, as shown in Figure 2 in conjunction with Figure 1, S102 includes S1021 to S1022. nothing. S1021 calculates the importance index for each node in the aquatic ecosystem network. Selectable, the above importance indicator is the node order of the node in the aquatic ecosystem network, node It can be calculated from centrality and node correlation. The above importance index satisfies the following formula: JPEG2026143293000034.jpg21122 Here, Ii represents the importance index of the i-th node in the aquatic ecosystem network, and k i represents the node degree of the i-th node in the aquatic ecosystem network, and B i represents the aquatic ecosyste m node centrality of the i-th node in the network, and C i represents the node correlation degree of the i-th no de in the aquatic ecosystem network, α1 represents the weight of node degree, and α2 represents the weight of node centrality, and α3 represents the weight of node correlation degree. In the above formula, node degree (Degree) is the number of edges through which one node is directly connected to other nodes. Node centrality (Betweenness Centrality ) measures the frequency of a single node appearing in all shortest paths, and is used in a network to reflect intermediate effects. Node correlation degree (Closeness Cen trality) measures the average shortest path length from one node to all other nodes , and is used to reflect the proximity between a node and other nodes. In graph theory, node degree , node centrality and node correlation degree are used to measure node importance and network structure as commonly used conventional indicators, and the definitions and applications of node degree, node centrality and node correlation degree all belong to the prior art. The calculation formula of node degree k i is as follows. JPEG2026143293000035.jpg4768Here, N represents the total number of nodes in the aquatic ecosystem network before removal, and a ij represents the adjacency matr ix of the aquatic ecosystem network. When node i is directly connected to node j, aij =1; when node i is not directly connected to node j, aij=0. Node centrality B i The formula for calculating this is as follows: JPEG2026143293000036.jpg3991 Here, n jk This represents the number of shortest paths between node j and node k, and n jk (i) is node i This represents the number of shortest paths passing between node j and node k via [a specific point]. Node correlation C i The formula for calculating this is as follows: JPEG2026143293000037.jpg6177 Here, C ij represents the degree of correlation between node i and node j, and d ji The aquatic ecosystem after removal This represents the shortest path length between node i and node j in the system network. In some cases, d max =max(d ij ) and d max This is the maximum value of all valid shortest path lengths in the network. Therefore, if there is no path connecting node i and node j within the network, The shortest path length between node i and node j is d max That is the case. S1022, Order of removal of nodes in the aquatic ecosystem network according to the magnitude of the importance index. Sort the order and match the predetermined removal ratio from the aquatic ecosystem network to at least Iterate through one node and at least one edge connected to at least one node. Remove the material, and after each removal operation, calculate the stability index value of the aquatic ecosystem network after each removal. Here, the pre-set removal ratio is determined by the number of nodes removed and the aquatic ecosystem before repeated removal. This is the ratio of the total number of nodes in the network. The above pre-set removal ratio is based on the actual aquatic environment. It can be set according to the situation. Selectable, pre-set removal ratios are 1%, It may be 2%, 3%, or 5%, or any other reasonable value, as in the embodiments of this application. This is not limited to this. For example, if the predetermined removal ratio is set to 1%, when the first removal is performed, The magnitude of the importance index value at the nodes of the aquatic ecosystem network before removal is the first 1% Remove the node, calculate the stability index value of the aquatic ecosystem network after the first removal, and When removal is carried out twice, the water ecosystem network after the first removal is repeated, and the water before the removal is repeated. The magnitude of the importance index value in the nodes of the ecosystem network is the first 1-2% of nodes ( The total number of removed nodes is 2% of the total number of aquatic ecosystem network nodes before iterative removal. Remove the material and calculate the stability index value of the aquatic ecosystem network after the second removal. The number of nodes removed was 3%, 4%, and 5% of the total number of aquatic ecosystem network nodes before iterative removal. ...If the result is 100%, the stability index value of the aquatic ecosystem network after removal is calculated. Please refer to Figure 3 for the removal process. As another example, if the predetermined removal ratio is set to 5%, the first removal will be performed. The magnitude of the importance index value at the nodes of the aquatic ecosystem network before repeated removal was initially Remove 5% of the nodes and calculate the stability index value of the aquatic ecosystem network after the first removal. And when the second removal is carried out, the aquatic ecosystem network after the first removal will be replaced with the pre-removal version. The magnitude of the importance index value at the nodes of the aquatic ecosystem network is the first 5-10% Nodes (At this time, the total number of nodes removed is the total number of nodes in the aquatic ecosystem network before iterative removal) (10% of the total) were removed, and the stability index value of the aquatic ecosystem network after the second removal was determined. Calculate. By analogy, the number of nodes removed is the total number of nodes in the aquatic ecosystem network before iterative removal. If the number is 15%, 20%, 25%...100%, then the aquatic ecosystem network after removal Calculate the stability index value. In the two examples above, removing a node also removes the edges connected to that node. I want you to understand this. Furthermore, the stability indicator values ​​of the aquatic ecosystem network after each removal are as follows: Network efficiency of the system network and maximum connectivity of each post-removal aquatic ecosystem network. Includes sex. The importance index value reflects the importance of each node in the aquatic ecosystem network. The higher the importance index value, the more important the node is in the aquatic ecosystem network. Please understand that removing the ne will have a significant impact on the aquatic ecosystem network. Network efficiency refers to the number of nodes in a network. Used to measure the efficiency of information transmission between devices. Maximum Connectivity (Maximum Conn Connectivity) is the size of the maximum connectivity subgraph in the network (i.e., the number of connections). It is the number of (codes). The network efficiency of the aquatic ecosystem network after each of the above removals satisfies the following equation. JPEG2026143293000038.jpg39122 Here, E represents the network efficiency of the aquatic ecosystem network after each removal, and N is , represents the total number of nodes in the aquatic ecosystem network before repeated removal, d ji The water areas after each removal This represents the shortest path length between node i and node j in the ecosystem network. The maximum connectivity of the aquatic ecosystem network after the above removal satisfies the following equation. JPEG2026143293000039.jpg3873 Here, S represents the maximum connectivity of each aquatic ecosystem network after removal, |C max | is, N is the total number of nodes in the maximum connectivity subgraph of each aquatic ecosystem network after removal, where N is the total number of nodes. , the total number of nodes in the aquatic ecosystem network before repeated removal. ru. example For example, early aquatic ecosystem networks consist of numerous nodes (fish, dissolved oxygen, nutrients, etc.) Assuming it is composed of these, these nodes are interconnected in a network Form a network. A key node in the aquatic ecosystem network (e.g., dissolved oxygen node). Removing the node will result in the direct correlations (e.g., fish survival and pollutant decomposition) being connected to that node. Relationships are destroyed, the connectivity of the network structure weakens, and the overall function of the ecosystem is impaired. It can be done. S103, until the stability index value of the aquatic ecosystem network after removal falls below the stability threshold. , the number of nodes removed from the aquatic ecosystem network and the number of nodes removed from the aquatic ecosystem network before repeated removals Based on the total number of nodes in the network, the threshold for the destruction resistance of the aquatic ecosystem network is determined, and the aquatic ecosystem Determining the threshold for the disruption resistance of the ecosystem network based on the results of a simulated analysis of the structural stability of aquatic ecosystem networks. do. Here, the stability threshold includes the network efficiency threshold and the maximum connectivity threshold. The range of the efficiency threshold value is 0.10 to 0.15, and the range of the maximum connectivity threshold value is 0.4 It is approximately 0.6. Selectively, S103 above includes the following: The network efficiency of the aquatic ecosystem network after removal is below the network efficiency threshold. or, if the maximum connectivity of the aquatic ecosystem network after removal is below the maximum connectivity threshold. In total, the number of nodes removed from the aquatic ecosystem network and the aquatic ecosystem network before repeated removals. From the total number of nodes in the work, the destruction resistance threshold of the aquatic ecosystem network is calculated, and the aquatic ecosystem The network's failure resistance threshold satisfies the following equation: JPEG2026143293000040.jpg3053 Here, fa represents the fracture resistance threshold, and Nr represents the removed noise of the aquatic ecosystem network. This represents the number of nodes, where N represents the total number of nodes in the aquatic ecosystem network before repeated removal. The above network efficiency threshold may be 0.14, and the above maximum connectivity threshold may be optional. The network efficiency threshold may be 0.4 or 0.6, and the embodiments of this application are as described above. This does not limit the maximum connectivity threshold mentioned above. In one application scenario, the aquatic environment is a waters area in Bishan, Chongqing, China. To analyze the stability of the ecological network. Regarding aquatic organisms contained in the aquatic environment, algae are... It includes Lorera, cyanobacteria, diatoms, and green algae, and the fish include grass carp, carp, crucian carp, and mackerel. The plants include water hyacinth, orchard grass, aquatic plants, goldfish weed, and reeds. Regarding environmental factors, water temperature is The temperature is 25℃, the pH is 7.5, the dissolved oxygen is 8 mg / L, and the hardness is 50 mg / L. It is L. Of the multiple types of pollutants contained in the aquatic environment, heavy metals include cadmium, lead and It contains copper, organic pollutants include organochlorine pesticides, and nutrients include nitrogen and phosphorus, in addition to Furthermore, it also contains zinc, petroleum-based pollutants, and polycyclic aromatic hydrocarbons. Specifically, the nitrogen concentration is 10m³. The concentration is g / L, the phosphorus concentration is 1.2 mg / L, and the concentration of organochlorine pesticides is 0.05 mg / L. The concentration is g / L, the concentration of petroleum-based pollutants is 2.3 mg / L, and the concentration of polycyclic aromatic hydrocarbons is 2.3 mg / L. The concentration was 3 mg / L, the cadmium concentration was 0.2 μg / L, and the lead concentration was 0.05 μg The concentration of copper is 0.15 μg / L, and the concentration of zinc is 0.3 μg / L. The concentration of heavy metal chromium is 0.1 μg / L. Applying method S101 provided by the embodiment to the above application scenario results in the aquatic environment Obtain the aquatic ecosystem network of the aquatic environment. The culprit is repeatedly removed, and S103 determines the threshold resistance to disruption of the aquatic ecosystem network in the aquatic environment. To decide. In the above process, the maximum connectivity of the aquatic ecosystem network after removal is the maximum connectivity threshold (0 5) If the amount is less than or equal to 1%, and the predetermined removal rate is 1%, then the aquatic ecosystem network From the number of removed nodes and the total number of nodes in the aquatic ecosystem network before repeated removal, It is configured to calculate the disruption resistance threshold of the regional ecosystem network. The removal rate, maximum connectivity of the aquatic ecosystem network, and network efficiency are shown in Table 1. It will be done.

[0008] Table 1: Removal Ratio, Maximum Connectivity, and Network Efficiency Diagram As can be seen from Table 1, when the removal rate reaches 24%, the aquatic ecosystem network after removal When the maximum connectivity of the system reaches 0.5001 and the removal rate reaches 26%, the water body after removal The maximum connectivity of the system network is 0.33435. In other words, the rejection rate reaches 24%. When this happens, the maximum connectivity of the aquatic ecosystem network after removal reaches the maximum connectivity threshold (0.5). At this point, the network function of the aquatic ecosystem network after removal is on the verge of collapse. This indicates that when the removal rate reaches 26%, the aquatic ecosystem network after removal When the maximum connectivity becomes less than the maximum connectivity threshold, the aquatic ecosystem network after removal This indicates that the network functionality of the system has collapsed, and as can be seen, the removal rate is 24%. The network function of aquatic ecosystem networks collapses when the level is around 26%. To more accurately indicate the stability of the aquatic ecosystem network, the removal rate reached 24%. When the aquatic ecosystem network reaches a destruction resistance threshold of 0.24 and a removal ratio of 26% Aquatic ecosystem network destruction resistance threshold range of 0.24~0 0.26 will be used as the result of a simulated analysis of the structural stability of aquatic ecosystem networks. In one embodiment, the above method further includes S104. S104, When the stability index value of the aquatic ecosystem network after removal falls below the stability threshold, Based on the number of nodes removed from the aquatic ecosystem network, the resilience of the aquatic ecosystem network Determine the resilience of aquatic ecosystem networks and the threshold for their resistance to disruption. This, combined, represents the results of a simulated analysis of the structural stability of aquatic ecosystem networks. Here, the resilience of aquatic ecosystem networks is the resistance of aquatic ecosystem networks to changes in the ecological environment. It is used to describe the resilience and recovery of the aquatic ecosystem network. The following conditions must be met: JPEG2026143293000042.jpg33110 Here, R represents the resilience of the aquatic ecosystem network, and Nr represents the resilience of the aquatic ecosystem network This represents the number of nodes removed from the twerk, d ji In the aquatic ecosystem network after removal This represents the shortest path length between node i and node j. In summary, the method of this application relates to a method that is brought about by natural or human action in actual circumstances. This system dynamically simulates changes in the aquatic environment (changes in species, environmental factors, and pollutants). This allows for the simulation analysis results of the stability of aquatic ecosystem network structures to be more accurate. It is possible.

Claims

1. A method for dynamically simulating the structural stability of aquatic ecosystem networks, comprising the following steps: Based on ecological information in aquatic environments, we will construct aquatic ecosystem networks. Here, the ecological information includes species information, environmental factors, and pollutant information, and the species information is the water This includes the spatial distribution of multiple species of aquatic organisms present in the local environment, and these multiple species of aquatic organisms include algae and fish. This includes species and aquatic plants, and the environmental factors are the temperature, pH value and dissolved oxygen of the aquatic environment. The pollutant information includes the concentrations of multiple types of pollutants contained in the aquatic environment, and the Multiple types of pollutants include heavy metals, organic pollutants, and nutrients, and the aquatic ecosystem network The element includes multiple nodes and multiple edges connecting the said multiple nodes. Based on the pre-set removal ratio and the importance index of each node in the aquatic ecosystem network Based on this, at least one node in the aquatic ecosystem network and at least The aquatic ecosystem after each removal is repeatedly removed from at least one edge connected to another node. Calculate the stability index value of the system network, After removal, until the stability index value of the aquatic ecosystem network falls below the stability threshold, Number of nodes removed from the aquatic ecosystem network and the aquatic ecosystem network before repeated removals Based on the total number of nodes in the work, the threshold for the destruction resistance of the aquatic ecosystem network is determined. The threshold for resistance to disruption of the aquatic ecosystem network is determined by the structural stability of the aquatic ecosystem network. The structural stability of aquatic ecosystem networks, characterized by being used as a result of simulated analysis. Dynamic simulation method.

2. Multiple nodes of the aquatic ecosystem network include multiple species of aquatic organisms and environmental factors. and correspond one-to-one with the aforementioned multiple types of pollutants, and one side of the aquatic ecosystem network is, This represents the ecological interaction relationship between two nodes connected by the aforementioned edge, and the ecological interaction relationship The relationships include predation, symbiosis, competition, survival threat, and mutual influence, where the ecological relationships are , the spatial distribution of multiple species of aquatic organisms present in the aquatic environment, the environmental factors and the aquatic environment The determination is based on the concentrations of multiple types of pollutants contained in the boundary, as described in 1. Method of description.

3. The aforementioned pre-set removal ratio and the importance indicator of each node in the aquatic ecosystem network Based on the target, at least one node in the aquatic ecosystem network and the few At least one edge connected to at least one node is repeatedly removed, and the water body after each removal The steps for calculating the stability index value of the ecosystem network include: The importance index for each node in the aquatic ecosystem network is calculated, and the importance index is as follows: The equation satisfies: Here, I ｉ This represents the importance index of the i-th node in the aquatic ecosystem network, k ｉ This represents the node order of the i-th node in the aquatic ecosystem network, and B ｉ teeth, This represents the node centrality of the i-th node in the aquatic ecosystem network, C ｉ is the water This represents the node correlation of the i-th node in the regional ecosystem network, and α １ , node degree Represents the weight of α ２ α represents the weight of node centrality, ３ This represents the weight of the node correlation. 、 The order in which nodes in the aquatic ecosystem network are removed is determined according to the magnitude of the aforementioned importance index. And, from the aquatic ecosystem network, at least the amount that matches the predetermined removal ratio. Iterate through one more node and at least one edge connected to at least one node. Remove the material, and after each removal operation, calculate the stability index value of the aquatic ecosystem network after each removal. Here, the stability index value of the aquatic ecosystem network after each removal is, Network efficiency of the regional ecosystem network and the aquatic ecosystem network after each removal This includes the maximum connectivity, where the preset removal ratio is proportional to the number of nodes in each removal. This is the ratio of the total number of nodes in the aquatic ecosystem network before removal, and the predetermined The removal ratio is 1%, 2%, 3%, or 5%. The network efficiency of the aquatic ecosystem network after each removal satisfies the following equation: Here, E represents the network efficiency of the aquatic ecosystem network after each removal, and N is , represents the total number of nodes in the aquatic ecosystem network before repeated removal, d ｊｉ This is after each removal. This represents the shortest path length between node i and node j in the aforementioned aquatic ecosystem network. The maximum connectivity of the aquatic ecosystem network after each removal satisfies the following equation: Here, S represents the maximum connectivity of the water ecosystem network after each removal, and |C ｍａｘ | This represents the total number of nodes in the maximum connectivity subgraph of each post-removal aquatic ecosystem network. The method according to feature 1.

4. The stability threshold includes a network efficiency threshold and a maximum connectivity threshold, and the network The range of the efficiency threshold value is 0.10 to 0.15, and the range of the maximum connectivity threshold value is 0 The method according to claim 3, characterized in that the value is 4 to 0.

6.

5. Until the stability index value of the aquatic ecosystem network after the removal falls below the stability threshold, The number of nodes removed from the aquatic ecosystem network and the number of nodes removed from the aquatic ecosystem network before repeated removal Based on the total number of nodes in the network, the threshold for the destruction resistance of the aquatic ecosystem network is determined. The steps include: The network efficiency of the aquatic ecosystem network after the removal is the network efficiency threshold. The value is less than or equal to the maximum connectivity of the aquatic ecosystem network after removal. If the number of nodes removed from the aquatic ecosystem network is below the large connectivity threshold, From the total number of nodes in the aquatic ecosystem network before repeated removal, The destruction resistance threshold of the aquatic ecosystem network is calculated, and the destruction resistance threshold of the aquatic ecosystem network satisfies the following equation. Tashi: Here, fa represents the fracture resistance threshold, and Nr represents the removal of the aquatic ecosystem network. This represents the number of nodes, where N represents the total number of nodes in the aquatic ecosystem network before repeated removal. The method according to feature 4.

6. The method described above further includes the following steps: When the stability index value of the aquatic ecosystem network after removal falls below the stability threshold, the water Based on the number of nodes removed from the regional ecosystem network, the aquatic ecosystem network Determine the toughness of the aquatic ecosystem network and the destruction of the aquatic ecosystem network. By adjusting the resistance threshold, the results of the simulated analysis of the structural stability of the aquatic ecosystem network are obtained, and the water The resilience of a regional ecosystem network is the tolerance of aquatic ecosystem networks to changes in the ecological environment. Used to describe resilience, the toughness of the aquatic ecosystem network is given by the following equation: Satisfied: Here, R represents the resilience of the aquatic ecosystem network, and Nr represents the resilience of the aquatic ecosystem network This represents the number of nodes removed from the twerk, d ｊｉ The aquatic ecosystem network after removal The method according to claim 1, characterized in that it represents the shortest path length between node i and node j in the middle. 。

Images