Soil Cd recognition and remediation method based on mulberry root exudates and microorganisms

Abstract

The invention relates to the technical field of soil cadmium remediation, in particular to a soil Cd recognition and remediation method based on mulberry root exudates and microorganisms. According to the method, spatial microscopic image frames of rhizosphere microflora are divided through a community structure, community node increment updating calculation of adjacent time sequences is executed, the transition recombination frequency of the rhizosphere microflora is obtained, and the transition color of the actual concentration content of root exudates is rendered based on a metal pollution level system and the transition recombination frequency. Generating a stress level heat model of the soil cadmium gradient in the target restoration area; and according to the stress level popularity model, carrying out popularity color analysis and identification on the soil cadmium pollution of the target restoration area, and generating a global stress identification result of the soil cadmium. According to the method, rapid identification of cadmium pollution enrichment can be realized by utilizing the synergistic response between mulberry root exudates and rhizosphere microorganisms, and a repair regulation and control plan is made, so that the cadmium pollution degree is reduced, the soil structure and fertility are optimized, and a safer and healthier environment is provided for plant growth.

Description

TECHNICAL FIELD

[0001] The application relates to the technical field of soil cadmium remediation, in particular to a soil Cd identification and remediation method based on mulberry root exudates and microorganisms. BACKGROUND

[0002] With the rapid development of industrialization and urbanization, cadmium (Cd) as a highly toxic heavy metal has become one of the main risk factors of soil pollution. Cadmium in soil has strong migration and bioavailability, and is easily accumulated through the food chain, causing serious threats to crop yield and food safety. The existing soil Cd identification method mainly judges the cadmium pollution degree through the enrichment or passivation effect of the whole plant, lacks a precise identification mechanism based on the interaction of mulberry root exudates and Cd metals, and is difficult to assist the above-mentioned identification mechanism in judging the precise pollution of soil Cd enrichment by means of the change of rhizosphere microbial community structure, and usually relies on manual operation tools for field identification, which is complicated and requires a large amount of labor cost, so that the final identification result of soil Cd enrichment has large deviation and inaccuracy. Secondly, the existing method usually relies on plant absorption or rhizosphere microbial transformation, and the efficiency is limited, which is difficult to cope with different concentration gradients of cadmium contaminated soil, and lacks a means of synergistically regulating mulberry root exudates, rhizosphere microbial community and soil cadmium migration and transformation process. In addition, the traditional method needs to rely on artificial experience to regulate the rhizosphere environment to reduce the biological availability of soil Cd, and the remediation efficiency is low, the cycle is long and the adaptability to different pollution environments is limited. Therefore, it is urgent to develop an innovative method and system combining the identification and remediation functions of mulberry root exudates and rhizosphere microorganisms to improve the accuracy of soil Cd pollution identification and the efficiency of remediation. SUMMARY

[0003] The application overcomes the shortcomings of the prior art and provides a soil Cd identification and remediation method based on mulberry root exudates and microorganisms.

[0004] To achieve the above purpose, the technical scheme adopted by the application is as follows: The application provides a soil Cd identification and remediation method based on mulberry root exudates and microorganisms, which comprises the following steps: S102: Determining the concentration content data of root exudates and the metabolic product data of rhizosphere microorganisms in the target remediation area, tracking and analyzing the migration enrichment channel of the metabolic product data based on the complexation of the concentration content data and Cd ions, and obtaining a soil Cd enrichment layout map of the target remediation area; S104: Obtain the transition reorganization frequency of the rhizosphere microbial community by dividing the spatial micrograph frames of the rhizosphere microbial community through community structure and performing adjacent time community node incremental update calculation, render the transition color of the actual concentration content of root exudates based on the metal pollution level system and the transition reorganization frequency, and generate a stress level heat model of the soil Cd gradient in the target remediation area; S106: According to the stress level heat model, analyze and identify the soil Cd pollution in the target remediation area by heat color, and generate a global stress identification result of soil Cd; S108: Obtain the resistance coefficient of rhizosphere microbial species to soil Cd according to the global stress identification result, and perform Gaussian modeling and analysis planning and deployment of the actual decline performance of rhizosphere microbial species based on the resistance coefficient to realize the optimal community recruitment deployment of rhizosphere microorganisms; S110: Establish a multiscale coupling model of rhizosphere environment-dissolved oxygen interface, introduce the oxidation and reduction characteristics of the optimal community recruitment deployment into the multiscale coupling model to calculate the patch pattern, and regulate the rhizosphere environment of the target remediation area according to the oxidation-reduction patch pattern to remediate soil Cd pollution.

[0005] More specifically, the S102 specifically includes the following steps: Obtain the target remediation area of soil Cd, and use biological test method to determine the root exudates and rhizosphere microorganisms of mulberry planted in the target remediation area; Through determination, obtain the species information of different mulberry root exudates and the actual concentration content of various information, and obtain the metabolic product distribution of each metabolic pathway formed by different rhizosphere microorganisms; Collect real-time environmental parameters of the target remediation area and obtain prior ionization characteristics of Cd ions under different preset environmental conditions, and retrieve the established migration efficiency of Cd ions in the target remediation area based on the prior ionization characteristics in the heavy metal knowledge graph according to the real-time environmental parameters; Based on the actual concentration content, construct a complex model of soil Cd-mulberry root exudates, and drive the upper and lower migration interaction hierarchical structure established by the candidate migration vector in the complex model through the established migration efficiency, to track and calculate the migration mode of the metabolic product distribution when the Cd ion is complexed with the root exudates, and determine the real-time migration channel of soil Cd; Detect the false touch signal flux of metal transport proteins misabsorbed by Cd ions in each real-time migration channel, if the false touch signal flux is greater than the preset threshold, mark the real-time migration channel as a dense enrichment channel; if the false touch signal flux is less than the preset threshold, mark the real-time migration channel as a sparse enrichment channel; According to the dense enrichment channel and the sparse enrichment channel, identify and plan the target remediation area, and obtain the soil Cd enrichment layout of the target remediation area.

[0006] More specifically, the complexation model of soil Cd-root exudate is constructed based on the actual concentration content, and the migration interaction layered structure of the upper and lower layers is established by driving the candidate migration vectors in the complexation model with the given migration efficiency, to track the migration mode of the metabolic product distribution when Cd ions are complexed with root exudates, and determine the real-time migration channel of soil Cd, which specifically includes the following steps: The complexation data of soil Cd and various types of root exudates corresponding to the actual concentration content in the historical time period is extracted through the soil Cd treatment log, and the complexation model of soil Cd-root exudate is constructed according to the complexation data, and a number of candidate migration vectors are preset based on the inhibition remodeling behavior of the complexation compound for Cd ions mapped by the complexation model; The candidate migration directions are divided into candidate upper and lower layers of migration interaction layered structure, and the given migration efficiency is injected into the candidate upper layer as a layered driving item to perform trend calculation; During the injection process, the elite neighborhood trend wind direction factor of the candidate upper layer is transmitted to the local framework of the metabolic product distribution deployment to guide the difference variation calculation of the candidate lower layer under the premise of the metabolic product distribution, to generate a potential migration test mode based on the metabolic regulation of Cd ions and rhizosphere microorganisms before complexation with root exudates, and determine the real-time migration channel of soil Cd according to the potential migration test mode.

[0007] More specifically, the S104 specifically includes the following steps: The rhizosphere microorganisms in the target repair area are captured and observed by a confocal laser microscope to obtain spatial microscopic image frames of the rhizosphere microbial community in time sequence; The static LPA algorithm is introduced to divide the community nodes of the spatial microscopic image frames at the current time sequence and the next time sequence, to obtain the current community network and the alternation community network, detect the incremental update of the current community network relative to the alternation community node, and output the disturbance index of the original community node; If the disturbance index is greater than the preset disturbance index, the original community node is replaced according to the community label with the highest appearance frequency in the neighbor community node at this time, and the above community node replacement step is repeated to perform local label propagation until the minimum disturbance coefficient is approached; The label results of all community nodes are merged after the propagation operation, the community division of the global spatial community is continuously tracked, the dynamic evolution result of the rhizosphere microbial community at different time steps is generated, and the transition reorganization frequency of the rhizosphere microbial community is determined according to the dynamic evolution result; The frequency of rhizosphere microbial community changes in different enrichment sub-regions in the soil Cd enrichment distribution map was obtained, and the empirical concentration of different root exudates to prevent negative interference from soil Cd pollution was obtained based on planting cases. The deviation between the actual concentration and the empirical concentration is calculated to obtain the concentration gradient difference. By identifying the concentration gradient difference in the heavy metal knowledge graph, the soil Cd metal pollution level system and the corresponding color scale registration rules that meet the concentration gradient difference are determined. Based on the frequency of change and recombination, the temporal structural changes of each rhizosphere microbial community are divided into rendering level intervals that follow time step jumps. Based on the color level registration rule, the actual concentration content is mapped to the transition color of the interpolated gradient of the rendering level interval, and a stress level heat model of soil Cd gradient in the target remediation area is generated.

[0008] More specifically, S106 includes the following steps: The heat color value range of soil Cd stress level in different enrichment sub-regions is obtained through the stress level heat model, and the toxicity heat color level is preset according to the soil Cd identification requirements. If the toxicity heat value scale does not cover any color value in the heat value range, then the enriched sub-region is marked as a low stress region of soil Cd. If the number of any color values ​​in the toxicity heat value range is less than a preset threshold, then the enriched sub-region is marked as a medium to low stress region for soil Cd. If the number of any color values ​​in the toxicity heat value range exceeds a preset threshold, then the enriched sub-region is marked as a medium to high stress region for Cd in the soil. If the toxicity heat value range completely covers any color value in the heat value range, then the enriched sub-region is marked as a high-stress region for soil Cd; and the global stress identification result for soil Cd is generated.

[0009] More specifically, S108 includes the following steps: The remediation requirements for soil Cd in the target remediation area are obtained. Based on the global stress identification results, multiple rhizosphere microbial species that meet the remediation requirements and the resistance coefficients of various species to soil Cd stress levels are obtained through big data retrieval. Based on the pre-defined expected stress identification results and expected decline rate of soil Cd activity according to the remediation requirements, the hash drift function between the global stress identification results and the expected stress identification results is calculated, and the actual decline rate of current soil Cd pollution is determined based on the hash drift function. The Gaussian process algorithm is introduced, and the actual decline rate is modeled in Gaussian based on the resistance coefficient. This generates a Gaussian regression model of the actual decline in soil Cd activity for rhizosphere microbial species. The resistance estimate and prediction distribution of different rhizosphere microbial species to the actual decline rate are calculated through the Gaussian regression model. By balancing the standard deviation and mean of the predicted distribution of resistance, the limiting reduction rate of soil Cd activity stress for each rhizosphere microbial species at different enrichment locations was constructed. If the limiting rate of decline is higher than the expected rate of decline, the rhizosphere microorganisms at the enrichment sites corresponding to the limiting rate of decline are considered suitable species; if the limiting rate of decline is lower than the expected rate of decline, the rhizosphere microorganisms at the enrichment sites corresponding to the limiting rate of decline are considered unsuitable species. By combining suitable and unsuitable species, the rhizosphere microbial community in the target remediation area was planned and analyzed to obtain the optimal community recruitment and deployment of rhizosphere microorganisms.

[0010] More specifically, S110 includes the following steps: Based on big data retrieval, the relative reduction characteristics of rhizosphere microbial species recruited and deployed in the target remediation area rhizosphere environment were obtained, as well as the relative oxidation characteristics of root exudates in the local environment were obtained. We introduce a dissolved oxygen flow model and a soil-specific redox model, defining the dissolved oxygen flow model as the macroscale and the soil-specific redox model as the microscale. At the same time, we construct a soil Cd environment map of the target remediation area based on the soil Cd enrichment distribution map. Based on the sampling criterion of transforming macroscopic conditions into microscopic constraints, the macroscopic and microscopic scales are co-mapped onto the soil Cd environment map, resulting in a multi-scale coupled model of the redox dynamics of the rhizosphere environment-dissolved oxygen interface. The interaction mechanism between root exudates and rhizosphere microorganisms was obtained. Based on the interaction mechanism, a coupling interaction strategy for dissolved oxygen in the soil environment during redox was constructed. Based on the coupling interaction strategy, a multi-scale coupling model was introduced to perform multi-scale state calculations on the relative oxidation and relative reduction characteristics from limited oxygen supply to rhizosphere microorganisms to severe oxygen consumption. Several particle behavior traces were output. During the multi-scale calculation, the effective response step size when each particle behavior trace representation is formed is continuously recorded, and the residence index of each particle behavior trace representation on the enriched and intersecting micro-region is planned and formulated according to the span of the effective response step size. If the residence index is less than the preset residence index, the enriched interlaced micro-region is ignored; if the residence index is greater than the preset residence index, patches are created on the enriched interlaced micro-region by stacking based on the characterization boundary of particle behavior trace characteristics, and finally the redox patch pattern of the rhizosphere environment is obtained. Based on the interaction mechanism, the rhizosphere environment of the target remediation area is regulated by planting according to the redox patch pattern, thereby regulating the Cd remediation ecology between root exudates and rhizosphere microorganisms and realizing soil Cd remediation.

[0011] A second aspect of the present invention also provides a soil Cd identification and remediation system based on mulberry root exudates and microorganisms. The system includes: a memory, a processor, and a communication interface. The memory includes a soil Cd identification and remediation method program based on mulberry root exudates and microorganisms. The communication interface is used for data connection and communication between the memory and the processor. When the soil Cd identification and remediation method program based on mulberry root exudates and microorganisms is executed by the processor, it implements the steps of the soil Cd identification and remediation method described in any one of the present invention.

[0012] This invention addresses the technical deficiencies in the prior art, and its beneficial technical effects are as follows: Concentration data of root exudates and rhizosphere microbial metabolites in the target remediation area were obtained by measuring and acquiring these data. Based on the complexation of Cd ions with the concentration data as a migration prior, the migration and enrichment channels forming the metabolite data were tracked and analyzed to obtain a soil Cd enrichment distribution map of the target remediation area. Spatial microscopic image frames of the rhizosphere microbial community were divided by community structure, and incremental update calculations of adjacent time-series community nodes were performed to obtain the frequency of change and recombination of the rhizosphere microbial community. Transitional colors of the actual concentration of root exudates were rendered based on the metal pollution level system and the frequency of change and recombination to generate a stress level heat model of the soil Cd gradient in the target remediation area. Based on the stress level heat model... The model analyzes and identifies soil Cd pollution in the target remediation area using thermal color analysis, generating a global Cd stress identification result. Based on this result, it obtains the resistance coefficients of rhizosphere microorganisms to Cd in the soil. Using these resistance coefficients, it performs Gaussian modeling and analysis to plan and deploy rhizosphere microorganisms based on their actual Cd activity reduction performance, obtaining the optimal community recruitment and deployment. A multi-scale coupling model of the rhizosphere environment-dissolved oxygen interface is established. The oxidation and reduction characteristics of the optimal community recruitment and deployment are imported into the multi-scale coupling model to calculate patch patterns. Based on these redox patch patterns, the rhizosphere environment in the target remediation area is regulated to remediate soil Cd pollution. This invention utilizes the connection between mulberry root exudates and rhizosphere microorganisms to achieve high-accuracy collaborative identification of soil Cd pollution in the target area, significantly improving identification efficiency. Based on the identification results, it formulates corresponding rhizosphere environment remediation directions and plans, thereby reducing soil Cd pollution levels and improving soil health and the plant growth environment. Attached Figure Description

[0013] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other embodiments can be obtained from these drawings without creative effort.

[0014] Figure 1 A flowchart of the first method for soil Cd identification and remediation based on mulberry root exudates and microorganisms is shown; Figure 2 A flowchart of a second method for soil Cd identification and remediation based on mulberry root exudates and microorganisms is shown. Figure 3 A system framework diagram of a soil Cd identification and remediation system based on mulberry root exudates and microorganisms is shown. Detailed Implementation

[0015] To better understand the above-mentioned objectives, features, and advantages of the present invention, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. It should be noted that, unless otherwise specified, the embodiments and features described in these embodiments can be combined with each other.

[0016] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and therefore the scope of protection of the invention is not limited to the specific embodiments disclosed below.

[0017] The first aspect of this invention provides a method for soil Cd identification and remediation based on mulberry root exudates and microorganisms, such as... Figure 1 As shown, it includes the following steps: S102: Measure and obtain the concentration data of root exudates and rhizosphere microbial metabolite data in the target remediation area. Based on the concentration data and the complexation of Cd ions as the migration prior, track and analyze the migration and enrichment channels that form the metabolite data to obtain the soil Cd enrichment layout map of the target remediation area. S104: The spatial microscopic image frames of the rhizosphere microbial community are divided by community structure and the incremental update calculation of community nodes in adjacent time series is performed to obtain the frequency of change and recombination of the rhizosphere microbial community. The transition color of the actual concentration of root exudates is rendered based on the metal pollution level system and the frequency of change and recombination to generate a stress level heat model of soil Cd gradient in the target remediation area. S106: Based on the stress level heat model, analyze and identify the heat color of Cd pollution in the soil of the target remediation area to generate global stress identification results for soil Cd. S108: Based on the global stress identification results, obtain the resistance coefficient of rhizosphere microorganisms to soil Cd. Based on the resistance coefficient, perform Gaussian modeling and analysis of the actual decline performance of soil Cd activity of rhizosphere microorganisms to plan and deploy them, and obtain the optimal community recruitment and deployment of rhizosphere microorganisms. S110: Establish a multi-scale coupling model of the rhizosphere environment-dissolved oxygen interface, import the oxidation and reduction characteristics of optimal community recruitment and deployment into the multi-scale coupling model to calculate the patch pattern, and regulate the rhizosphere environment of the target remediation area according to the redox patch pattern to remediate soil Cd pollution.

[0018] More specifically, step S102 includes the following steps: S102 specifically includes the following steps: The target remediation area for soil Cd was obtained, and the root exudates and rhizosphere microorganisms of mulberry trees planted in the target remediation area were measured using a bioassay. By measuring, we obtained information on the types of secretions from different mulberry roots and the actual concentration of each type, as well as the distribution of metabolites from different rhizosphere microorganisms in various metabolic pathways. Real-time environmental parameters of the target remediation area are collected and the prior free characteristics of Cd ions under different preset environmental conditions are obtained. Based on the prior free characteristics, the predetermined migration efficiency of Cd ions in the target remediation area is retrieved in the heavy metal knowledge graph based on the real-time environmental parameters. Based on the actual concentration, a complexation model of soil Cd-mulberry root exudates was constructed. The migration interaction layer structure of the upper and lower layers was established by driving the candidate migration vectors in the complexation model through a predetermined migration efficiency. This was used to track and calculate the migration pattern of the distribution of metabolites generated when Cd ions complex with root exudates, and to determine the real-time migration channels of soil Cd. The false trigger signal flux of metal transport proteins being mistakenly absorbed by Cd ions in each real-time migration channel is detected. If the false trigger signal flux is greater than a preset threshold, the real-time migration channel is marked as a dense enrichment channel; if the false trigger signal flux is less than the preset threshold, the real-time migration channel is marked as a sparse enrichment channel. Based on the dense enrichment channels and sparse enrichment channels, the target remediation area is identified and the soil Cd enrichment distribution map of the target remediation area is obtained.

[0019] It should be noted that the Chinese name for soil Cd is soil cadmium. Information on the types of exudates from different mulberry root systems includes citric acid, malic acid, oxalic acid, amino acids, and polyphenols. Rhizosphere microorganisms include phosphate-solubilizing bacteria, nitrogen-fixing bacteria, and Bacillus. Because these mulberry root exudates can promote the continuous migration of divalent cadmium ions in the soil through complexation, forming different enrichment patterns, cadmium ions are usually in a free state. However, existing identification methods often struggle to accurately track the migration trajectory of soil Cd, leading to a significant decrease in the accuracy of soil Cd enrichment identification results. To address this, this method retrieves the corresponding migration efficiency from a heavy metal knowledge graph based on the prior free characteristics of cadmium ions. This predetermined migration efficiency measures the efficiency and performance of spontaneous migration of cadmium ions under target remediation conditions and is a key indicator for tracking cadmium ion migration trajectories. Since the migration of cadmium ions in soil is driven by the complexation effect of mulberry root exudates, a model of soil Cd combined with root exudates to form a complex is constructed as the underlying navigation logic. This makes the tracking of cadmium ion migration more in line with the natural laws of the soil root environment, significantly improving the reliability and credibility of cadmium ion migration exploration. Subsequently, the candidate migration vectors mapped by this complexation model are divided into an interactive hierarchical structure of upper and lower layers. Because certain rhizosphere microorganisms (such as phosphate-solubilizing bacteria, nitrogen-fixing bacteria, and Bacillus) induce the synthesis of metal-binding proteins, glutathione, or extracellular polysaccharides when they encounter cadmium ions migrating in the soil driven by mulberry root exudates, this method uses a predetermined migration efficiency setting as the driving force for the input interaction hierarchy to guide the migration wind direction from the upper layer to the lower layer to transfer the distribution of metabolites in a fixed differential variation. The migration test pattern generated by the differential variation is the wandering channel created when soil cadmium ions migrate in the rhizosphere environment of the target remediation area, which constitutes the precise migration vector field of soil cadmium, providing an accurate and reliable location basis for the enrichment location and direction of soil cadmium.

[0020] It should be noted that mulberry roots can sense the presence of Cd in the soil environment through ion channels and transport proteins. This sensing is mainly triggered by the misabsorption of Cd ions by metal transport proteins on plant cell membranes. Therefore, if the misabsorption signal flux is greater than a preset threshold, it indicates that the metal transport proteins can sense cadmium ions at a high frequency within this real-time migration channel, suggesting a dense enrichment trend of Cd ions in the soil, which is a dense enrichment phenomenon of Cd ions. Thus, this real-time migration channel is labeled as a dense enrichment channel. Conversely, if the flux is lower, it indicates a sparse enrichment trend of Cd ions, and this real-time migration channel is labeled as a sparse enrichment channel. This method can utilize the complexation metabolic characteristics between mulberry root exudates, rhizosphere microorganisms, and soil Cd ions to track their migration, thereby accurately tracing the enrichment and distribution of Cd in the soil and providing high-precision location guidance for subsequent soil Cd pollution identification.

[0021] More specifically, the method of constructing a soil Cd-root exudate complexation model based on actual concentration content, and using a predetermined migration efficiency to drive the upper and lower migration interaction hierarchical structure established by candidate migration vectors in the complexation model, is used to track and calculate the migration pattern of metabolites generated when Cd ions complex with root exudates, thereby determining the real-time migration channels of soil Cd. This method specifically includes the following steps: By extracting the complexation data of soil Cd and various types of root exudates corresponding to the actual concentration of soil Cd in historical time periods from soil Cd remediation logs, a soil Cd-root exudate complexation model was constructed based on the complexation data. Based on the inhibitory remodeling behavior of complexes on Cd ions mapped by the complexation model, several candidate migration vectors were preset. Several candidate migration vectors are divided into a hierarchical structure of migration interaction between candidate upper layer and candidate lower layer. The predetermined migration efficiency is injected into the candidate upper layer to perform trend inference. During the injection process, trend direction factors are transmitted from the elite neighborhood of the candidate upper layer to the local framework of metabolite distribution deployment to guide the differential variation calculation of the candidate lower layer under the premise of metabolite distribution. This generates a potential migration test model based on the metabolic regulation of Cd ions and rhizosphere microorganisms under the precursor of root exudate complexation. The real-time migration channel of soil Cd is determined based on the potential migration test model.

[0022] It should be noted that, firstly, this method establishes a complexation model. The overall architecture of this model maps the migration activity of root exudates after they inhibit the morphological shaping behavior of Cd ions by complexing with soil. Based on this, candidate guide individuals, i.e., candidate migration vectors, are randomly generated in the migration solution space, providing a diverse exploration starting network for global migration channel tracking. Next, the candidate migration vectors are divided into a candidate upper-layer exploration framework that stores global elite individuals and a candidate lower-layer development framework that accommodates local ordinary individuals. The candidate upper layer is mainly responsible for performing a global calculation of the migration trend of soil Cd ions based on a predetermined migration efficiency. This allows the exploration of migration trajectories to move as far as possible towards the objective function that satisfies the migration efficiency characteristics of soil Cd ions, thereby exploring potential elite guides that conform to the soil Cd migration trend within the target remediation area, improving the accuracy and coverage of Cd ion migration tracking, and avoiding enrichment identification omissions and orientation errors. During the global extrapolation of predetermined migration efficiency, the elite domain continuously output by the candidate upper layer generates a trend indicator factor. This trend indicator factor is the implicit vector information of the global migration trend, containing the trend direction, magnitude, and scalar indicator. It is used as an elite guide representative to be transmitted to the candidate lower layer. Based on the phenomenon that rhizosphere microorganisms induce the synthesis of metal-binding proteins, glutathione, or extracellular polysaccharides when they encounter cadmium ions driven by root exudates in the soil, this method deploys a local architecture in the candidate lower layer that follows the premise of metabolite distribution. This allows the distribution of metabolites to provide accurate and known migration enrichment guidance for the differential variation of the candidate lower layer led by the trend indicator factor. This can effectively accelerate the local convergence of the candidate lower layer in the reasonable enrichment area, further ensuring the precision of soil Cd enrichment identification and the accuracy of traceability. This method innovatively establishes a hierarchical driven tracking structure that transmits from the upper layer, which has the ability to explore the migration trend across the entire region, to the lower layer, which has the ability to guide development locally. This addresses the technical shortcoming of existing methods in lacking real-time tracking and identification of soil Cd, and achieves dynamic enrichment and localization of soil Cd. Compared with the cumbersome process of traditional methods that calculate step by step, this method can effectively enhance the efficiency of computation and identification, and ensure the accuracy and reliability of subsequent identification results.

[0023] More specifically, S104 includes the following steps: The rhizosphere microorganisms in the target repair area were captured and observed using a confocal laser microscope to obtain spatial microscopic image frames of the rhizosphere microbial community in a continuous time sequence. The static LPA algorithm is introduced to perform structural partitioning of community nodes in spatial microscopic image frames located in the current time series and the next time series, to obtain the current community network and the iterative community network. The incremental update of the current community network relative to the iterative community network is detected, and the perturbation index of the original community nodes is output. If the disturbance index is greater than the preset disturbance index, then the original community node is replaced with the community label corresponding to the highest frequency of appearance among the neighboring community nodes. The above community node replacement steps are repeated to perform local label propagation until the minimum disturbance coefficient is approached. By merging the label results of all community nodes after the propagation operation, the community division of the global spatial community is continuously tracked, and the dynamic evolution results of the rhizosphere microbial community at different time steps are generated. The frequency of change and recombination of the rhizosphere microbial community is determined based on the dynamic evolution results. The frequency of rhizosphere microbial community changes in different enrichment sub-regions in the soil Cd enrichment distribution map was obtained, and the empirical concentration of different root exudates to prevent negative interference from soil Cd pollution was obtained based on planting cases. The deviation between the actual concentration and the empirical concentration is calculated to obtain the concentration gradient difference. By identifying the concentration gradient difference in the heavy metal knowledge graph, the soil Cd metal pollution level system and the corresponding color scale registration rules that meet the concentration gradient difference are determined. Based on the frequency of change and recombination, the temporal structural changes of each rhizosphere microbial community are divided into rendering level intervals that follow time step jumps. Based on the color level registration rule, the actual concentration content is mapped to the transition color of the interpolated gradient of the rendering level interval, and a stress level heat model of soil Cd gradient in the target remediation area is generated.

[0024] It should be noted that the community structure of rhizosphere microorganisms undergoes recombination changes with varying frequencies under specific Cd-stress pollution in soil. Therefore, this method acquires continuous spatial microscopic images of the rhizosphere microorganism community using confocal laser microscopy. These adjacent spatial microscopic images are then divided into community structures, forming a community network of the rhizosphere microorganism community structure distribution, namely the current community network and the iterative community network. This community networking provides an effective and stable focal reference and starting point for the spatial community evolution process, enhancing the interpretability and recombination continuity of community changes. Subsequently, incremental changes in the community network structure between adjacent time series are detected during the temporal evolution process, including the relative addition or deletion of community nodes or edges, suitable for both minute rhizosphere microbial activities and large-scale microbial communities in the target remediation area. This allows for the localization or locking of only the areas (nodes and their neighbors) affected by the incremental changes, uncovering the global perturbation effect of the temporal changes in the rhizosphere microorganism community. If the disturbance index is greater than the preset disturbance index, it indicates that the fluctuation range of the rhizosphere microbial community structure on the original community node is large, indicating that the local community structure is frequently reorganized and changed. Therefore, the community labels on the original nodes where the local community structure is located are updated according to the label distribution of neighboring community nodes, so that newly added community nodes in the next time series can be quickly assigned to suitable community communities, and the community instability caused by the deletion of nodes or edges can be quickly repaired, thereby ensuring the dynamic adaptability of time series community division, improving the accuracy of the dynamic evolution of rhizosphere microbial community, and thus effectively and accurately capturing the frequency of rhizosphere microbial community change and reorganization. This frequency of change and reorganization reflects the pollution response of rhizosphere microbial community to the soil Cd environment in the target remediation area. Therefore, this method uses it as a reliable biological indicator for soil Cd identification.

[0025] It should be noted that, since root exudates exhibit concentration variations and species shifts under specific soil Cd concentrations, this method calculates the concentration gradient difference between the actual concentration of root exudates in the current target remediation area and the empirical concentration of different root exudates unaffected by soil Cd pollution. This concentration gradient difference represents the quality change of root exudates and is a key identification signal for inferring the pollution level. Therefore, this method further identifies a soil Cd metal pollution level system that satisfies the concentration gradient difference and the corresponding color scale registration rules to determine the enrichment of soil Cd pollution in the target remediation area. Specifically, the transition and recombination frequency reveals the transient rate of soil Cd enrichment at the rhizosphere microbial level. Therefore, based on the transition and recombination frequency, a rendering level interval following the time step jump (transient rate) is constructed, and the color scale registration rules are further used to assign values ​​to the transition colors of the actual concentration within the rendering level interval. Finally, the enrichment pollution level of soil Cd is visualized in the form of a gradient color gamut, realizing rapid identification and visualization of soil Cd enrichment pollution. This method enables precise identification of soil Cd enrichment levels based on the synergistic response of mulberry root exudate quality and rhizosphere microbial community. It replaces the cumbersome steps and errors of traditional manual identification tools, saving time and effort, significantly improving identification efficiency, and ensuring the reliability and accuracy of the identification results.

[0026] More specifically, S106 includes the following steps: The heat color value range of soil Cd stress level in different enrichment sub-regions is obtained through the stress level heat model, and the toxicity heat color level is preset according to the soil Cd identification requirements. If the toxicity heat value scale does not cover any color value in the heat value range, then the enriched sub-region is marked as a low stress region of soil Cd. If the number of any color values ​​in the toxicity heat value range is less than a preset threshold, then the enriched sub-region is marked as a medium to low stress region for soil Cd. If the number of any color values ​​in the toxicity heat value range exceeds a preset threshold, then the enriched sub-region is marked as a medium to high stress region for Cd in the soil. If the toxicity heat value range completely covers any color value in the heat value range, then the enriched sub-region is marked as a high-stress region for soil Cd; and the global stress identification result for soil Cd is generated.

[0027] It should be noted that if the toxicity heat value scale does not cover any color value within the heat value range, it indicates that the heat value within the enriched sub-region is not covered by the toxicity heat value scale range. This means that the soil Cd stress level within the enriched sub-region does not meet or reach the high toxicity standard for plant growth, and therefore the enriched sub-region is displayed as a low-stress area on the soil Cd identification terminal. If the number of any color value within the heat value range in the toxicity heat value scale is less than a preset threshold, it indicates that a small portion of the heat value within the enriched sub-region is covered by the toxicity heat value scale range. This means that the soil Cd pollution in the enriched sub-region is trending towards toxicity but the process or amount is small, resulting in a low level of pollution. Therefore, the enriched sub-region is marked as a moderate to low-stress area. The remaining analysis cases are the opposite but similar, and will not be elaborated here. This method can identify and analyze different pollution levels of soil Cd in the target remediation area based on the stress level heat model, thereby achieving accurate identification of soil Cd enrichment, achieving full coverage of pollution location, and providing planters with a highly reliable basis for Cd remediation planning.

[0028] More specifically, S108, as Figure 2 As shown, the specific steps include: The remediation requirements for soil Cd in the target remediation area are obtained. Based on the global stress identification results, multiple rhizosphere microbial species that meet the remediation requirements and the resistance coefficients of various species to soil Cd stress levels are obtained through big data retrieval. Based on the pre-defined expected stress identification results and expected decline rate of soil Cd activity according to the remediation requirements, the hash drift function between the global stress identification results and the expected stress identification results is calculated, and the actual decline rate of current soil Cd pollution is determined based on the hash drift function. The Gaussian process algorithm is introduced, and the actual decline rate is modeled in Gaussian based on the resistance coefficient. This generates a Gaussian regression model of the actual decline in soil Cd activity for rhizosphere microbial species. The resistance estimate and prediction distribution of different rhizosphere microbial species to the actual decline rate are calculated through the Gaussian regression model. By balancing the standard deviation and mean of the predicted distribution of resistance, the limiting reduction rate of soil Cd activity stress for each rhizosphere microbial species at different enrichment locations was constructed. If the limiting rate of decline is higher than the expected rate of decline, the rhizosphere microorganisms at the enrichment sites corresponding to the limiting rate of decline are considered suitable species; if the limiting rate of decline is lower than the expected rate of decline, the rhizosphere microorganisms at the enrichment sites corresponding to the limiting rate of decline are considered unsuitable species. By combining suitable and unsuitable species, the rhizosphere microbial community in the target remediation area was planned and analyzed to obtain the optimal community recruitment and deployment of rhizosphere microorganisms.

[0029] It should be noted that some rhizosphere microorganisms can alter the physicochemical properties of soil by releasing their metabolites into soil locations with high Cd concentrations, thereby affecting the speciation and activity of high-concentration cadmium. They can also directly adsorb Cd ions through polysaccharides, proteins, and functional groups on their cell walls. Therefore, recruiting suitable rhizosphere microorganisms at different enrichment locations plays a crucial role in the comprehensive remediation of Cd in soil. To this end, this method retrieves multiple different rhizosphere microbial species and their corresponding resilience coefficients based on global stress identification results. These resilience coefficients reflect the remediation performance of rhizosphere microbial species against different levels of Cd stress in the soil. Then, based on remediation requirements, the expected stress identification results and the expected reduction rate of soil Cd activity are preset. This expected reduction rate is an objective function, serving as the supporting anchor point for rhizosphere microbial suitability screening. The expected stress identification results represent the reasonable stress distribution identified after remediating the soil Cd pollution in the target remediation area to the expected state. The actual reduction magnitude represents the order of magnitude and mileage by which rhizosphere microorganisms need to reduce soil Cd activity, clarifying the target span of the remediation process. Subsequently, a Gaussian process was used to model the actual reduction magnitude. Based on this Gaussian inference, the predicted mean and standard deviation of rhizosphere microorganisms for remediating Cd pollution in soil were derived under different resistance coefficient constraints. The predicted mean represents the Gaussian model's estimate of the actual reduction magnitude, while the predicted standard deviation represents the uncertainty of whether the rhizosphere microorganism can completely remediate the pollution. This method not only obtains the remediation achievement predictions for different individual rhizosphere microorganisms in advance but also quantifies the uncertainty of the predictions, effectively achieving a balance between remediation exploration and recruitment and utilization. This clarifies the maximum reduction benefit of each rhizosphere microorganism species for Cd activity stress in soil at different enrichment locations, i.e., the limit reduction rate. If the limit reduction rate is higher than the expected reduction rate, it indicates that the rhizosphere microorganism species can efficiently reduce soil Cd activity at the current enrichment location, demonstrating a high degree of suitability for Cd remediation at that location, and thus it is recruited as a suitable species. Conversely, if the limit reduction rate is lower, it indicates low reduction efficiency and clearly cannot effectively inhibit Cd activity at the current enrichment location, thus it is considered an unsuitable species and is subject to replacement or elimination recruitment planning. This method enables the selective recruitment of suitable root microorganisms for the remediation of Cd pollution in the target remediation area, thereby reducing the bioavailability of Cd in the soil and further improving the remediation benefits of rhizosphere microorganisms on the soil.

[0030] More specifically, S110 includes the following steps: Based on big data retrieval, the relative reduction characteristics of rhizosphere microbial species recruited and deployed in the target remediation area rhizosphere environment were obtained, as well as the relative oxidation characteristics of root exudates in the local environment were obtained. We introduce a dissolved oxygen flow model and a soil-specific redox model, defining the dissolved oxygen flow model as the macroscale and the soil-specific redox model as the microscale. At the same time, we construct a soil Cd environment map of the target remediation area based on the soil Cd enrichment distribution map. Based on the sampling criterion of transforming macroscopic conditions into microscopic constraints, the macroscopic and microscopic scales are co-mapped onto the soil Cd environment map, resulting in a multi-scale coupled model of the redox dynamics of the rhizosphere environment-dissolved oxygen interface. The interaction mechanism between root exudates and rhizosphere microorganisms was obtained. Based on the interaction mechanism, a coupling interaction strategy for dissolved oxygen in the soil environment during redox was constructed. Based on the coupling interaction strategy, a multi-scale coupling model was introduced to perform multi-scale state calculations on the relative oxidation and relative reduction characteristics from limited oxygen supply to rhizosphere microorganisms to severe oxygen consumption. Several particle behavior traces were output. During the multi-scale calculation, the effective response step size when each particle behavior trace representation is formed is continuously recorded, and the residence index of each particle behavior trace representation on the enriched and intersecting micro-region is planned and formulated according to the span of the effective response step size. If the residence index is less than the preset residence index, the enriched interlaced micro-region is ignored; if the residence index is greater than the preset residence index, patches are created on the enriched interlaced micro-region by stacking based on the characterization boundary of particle behavior trace characteristics, and finally the redox patch pattern of the rhizosphere environment is obtained. Based on the interaction mechanism, the rhizosphere environment of the target remediation area is regulated by planting according to the redox patch pattern, thereby regulating the Cd remediation ecology between root exudates and rhizosphere microorganisms and realizing soil Cd remediation.

[0031] It should be noted that in soil areas near root hairs or aerenchyma, plants can transport oxygen outward through their roots, creating a relatively oxidizing local environment. Simultaneously, certain rhizosphere microorganisms rapidly consume oxygen, creating a locally enhanced reducing environment in areas of high oxygen consumption but insufficient oxygen supply. This causes Cd to be reduced to a more stable form or combine with sulfur to precipitate. Therefore, this method establishes a multi-scale coupled model of the redox dynamics at the rhizosphere environment-dissolved oxygen interface. Specifically, using a soil Cd environment map of the target remediation area as the underlying framework, a mapping relationship is established between the macroscopic scale of the dissolved oxygen flow model and the microscopic scale of the soil environment-specific redox model. This achieves the interaction between the macroscopic flow of dissolved oxygen and the microscopic redox dynamics, ensuring the coupling and communication between macroscopic airflow and microscopic redox reactions (micro-reactions – diffusion). The specific mapping relationship refers to sampling criteria that transform macroscopic conditions into microscopic constraints. For example, the theoretical parameters of the macroscopic flow medium of dissolved oxygen are transformed into particle distribution constrained by the microscopic redox exchange between the root system and the rhizosphere environment. The coupling interaction strategy constructed by the interaction mechanism determines the interaction mode between models at different scales, ensuring that the multi-scale calculation between macroscopic and microscopic scale models fully follows the interaction effect of root exudates and rhizosphere microorganisms, and ensuring the stability and convergence of iterative coupling calculations. Ultimately, the multi-scale coupling model can infer the behavioral traces of dissolved oxygen particles from the redox level. For example, it shows a short-term retention of dissolved oxygen on the soil surface during the process from root oxygen supply to rapid oxygen consumption by rhizosphere microorganisms, providing a reliable planning basis for redox patch patterns.

[0032] It should be noted that the effective response step length during the formation of each particle behavior trace reflects the duration of the soil redox response of dissolved oxygen in the target remediation area, while its span reveals the residence characteristics of dissolved oxygen redox in the soil. If the residence index is less than the preset residence index, it indicates that dissolved oxygen is unlikely to form redox patches in this enriched and intertwined micro-region, belonging to an unreasonable patch planning area; conversely, it indicates that patches are easy to form, ensuring the bioavailability of Cd in the soil at this location. Finally, the rhizosphere environment of the target remediation area is regulated according to the redox patch pattern. This method can determine the appropriate redox patch area for the optimal community recruitment and deployment of selectively recruited rhizosphere microorganisms, thus clarifying the remediation targets of different rhizosphere microbial species in soil Cd enrichment, thereby reducing the bioavailability of soil Cd in the rhizosphere environment, enhancing the speciation reduction performance, fixation and migration control capabilities of cadmium, improving soil health and plant growth environment, and achieving green and sustainable soil Cd remediation benefits.

[0033] The second aspect of this invention also provides a soil Cd identification and remediation system based on mulberry root exudates and microorganisms, such as... Figure 3As shown, the system includes: a memory 301, a processor 302, and a communication interface 303. The memory 301 includes a soil Cd identification and remediation method program based on mulberry root exudates and microorganisms. The communication interface 303 is used for data connection and communication between the memory 301 and the processor 302. When the soil Cd identification and remediation method program based on mulberry root exudates and microorganisms is executed by the processor 302, it implements any of the steps of the soil Cd identification and remediation method described above.

[0034] The above are merely specific embodiments of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A soil Cd identification and remediation method based on mulberry root exudates and microorganisms, characterized in that, The method comprises the following steps: S102: Determine the concentration content data of root exudates in the target remediation area and the rhizosphere microbial metabolite data, track and analyze the migration enrichment channels of the metabolite data based on the complexation of the concentration content data and Cd ions, and obtain the soil Cd enrichment layout of the target remediation area; S104: Divide the spatial micrograph frames of the rhizosphere microbial community by community structure and perform adjacent time sequence community node incremental update calculation to obtain the transition reorganization frequency of the rhizosphere microbial community, render the transition color of the actual concentration content of the root exudates based on the metal pollution level system and the transition reorganization frequency, and generate the stress level heat model of the soil Cd gradient in the target remediation area; S106: According to the stress level heat model, analyze and identify the heat color of the soil Cd pollution in the target remediation area, and generate the global stress identification result of the soil Cd; S108: According to the global stress identification result, obtain the resistance coefficient of the rhizosphere microbial species to soil Cd, and perform Gaussian modeling and analysis planning and deployment of the soil Cd activity actual decline performance of the rhizosphere microbial species based on the resistance coefficient, to obtain the optimal community recruitment deployment of the rhizosphere microorganisms; S110: Establish a multi-scale coupling model of the rhizosphere environment-dissolved oxygen interface, introduce the oxidation and reduction characteristics of the optimal community recruitment deployment into the multi-scale coupling model to calculate the patch pattern, and regulate the rhizosphere environment of the target remediation area according to the oxidation-reduction patch pattern to repair the soil Cd pollution.

2. The method for identifying and remediating soil Cd based on mulberry root exudates and microorganisms according to claim 1, characterized in that, The S102 specifically comprises the following steps: Obtain the target remediation area of soil Cd, and use biological testing method to determine the root exudates and rhizosphere microorganisms of mulberry planted in the target remediation area; Through determination, obtain the type information of different mulberry root exudates and the actual concentration content of various type information, and obtain the metabolic product distribution of each metabolic pathway formed by different rhizosphere microorganisms; Collect real-time environmental parameters of the target remediation area and obtain the prior free characteristics of Cd ions under different preset environmental conditions, and search for the established migration efficiency of Cd ions in the target remediation area in the heavy metal knowledge graph based on the prior free characteristics and real-time environmental parameters; Based on the actual concentration content, a complexation model of soil Cd-mulberry root exudates is constructed, and the upper and lower migration interaction hierarchical structures established by the candidate migration vectors in the complexation model are driven by the established migration efficiency to track and calculate the migration mode of the metabolic product distribution when Cd ions are complexed with root exudates, so as to determine the real-time migration channel of soil Cd; Detect the false touch signal flux of the metal transport protein misabsorbed by Cd ions in each real-time migration channel, if the false touch signal flux is greater than the preset threshold, mark the real-time migration channel as a dense enrichment channel; if the false touch signal flux is less than the preset threshold, mark the real-time migration channel as a sparse enrichment channel; According to the dense enrichment channel and the sparse enrichment channel, the target remediation area is identified and planned to obtain the soil Cd enrichment layout of the target remediation area.

3. The method according to claim 2, wherein the soil Cd identification and remediation method based on mulberry root exudates and microorganisms is characterized by, The complexation model of soil Cd-root exudate is constructed based on the actual concentration content, and the migration interaction layered structure of the upper layer and the lower layer is established by driving the candidate migration vector in the complexation model with the given migration efficiency, so as to track the migration mode of the metabolic product distribution when the Cd ion is complexed with the root exudate, and determine the real-time migration channel of the soil Cd, and the specific steps include the following steps: The complexation data of the actual concentration content of the soil Cd and various types of root exudate in the historical time period is extracted through the soil Cd treatment log, and the complexation model of the soil Cd-root exudate is constructed according to the complexation data, and a plurality of candidate migration vectors are preset based on the inhibition remodeling behavior of the complexation product to the Cd ion mapped by the complexation model; A plurality of candidate migration directions are divided into a migration interaction layered structure of a candidate upper layer and a candidate lower layer, and a given migration efficiency is injected into the candidate upper layer as a layered driving item to perform trend calculation; During the injection process, the elite neighborhood of the candidate upper layer transmits the trend wind direction factor to the local framework of the metabolic product distribution to guide the difference variation calculation of the candidate lower layer under the premise of the metabolic product distribution, generate a potential migration test mode based on the metabolic regulation of the Cd ion and the rhizosphere microorganism before the complexation of the root exudate, and determine the real-time migration channel of the soil Cd according to the potential migration test mode.

4. The method for identifying and remediating Cd in soil based on mulberry root exudates and microorganisms according to claim 1, characterized in that, The S104 specifically includes the following steps: The rhizosphere microorganism in the target repair area is captured and observed by a confocal laser microscope to obtain spatial microscopic image frames of the rhizosphere microbial community in time sequence; A static LPA algorithm is introduced to divide the community nodes of the spatial microscopic image frames at the current time sequence and the next time sequence to obtain the current community network and the alternation community network, detect the incremental update of the current community network relative to the alternation community node, and output the disturbance index of the original community node; If the disturbance index is greater than the preset disturbance index, the community node is replaced according to the community label corresponding to the neighbor community node with the highest appearance frequency at this time, and the above community node replacement step is repeated to perform local label propagation until the minimum disturbance coefficient is approached; The label results of all community nodes are merged after the propagation operation, the community division of the global spatial community is continuously tracked, the dynamic evolution result of the rhizosphere microbial community at different time steps is generated, and the transition reorganization frequency of the rhizosphere microbial community is determined according to the dynamic evolution result; The transition reorganization frequency of the rhizosphere microbial community in different enrichment sub-regions in the soil Cd enrichment layout is obtained, the experienced concentration content of different root exudates to prevent negative interference of soil Cd pollution is obtained based on the planting case; The deviation between the actual concentration content and the experienced concentration content is calculated to obtain the concentration content gradient difference, the concentration content gradient difference is identified in the heavy metal knowledge graph, and the soil Cd metal pollution grade system and the corresponding color scale registration rule that meet the concentration content gradient difference are determined; The stress level heat model of the soil Cd gradient in the target repair area is generated by dividing the time sequence structure change of each rhizosphere microbial community into a rendering level interval following a time step jump based on transition recombination frequency, and mapping the actual concentration content to a transition color of a rendering level interval interpolation gradient based on a color scale registration rule.

5. The method for identifying and remediating soil Cd based on mulberry root exudates and microorganisms according to claim 1, characterized in that, The S106 specifically includes the following steps: A heat color value range of the soil Cd stress level in different enrichment sub-regions is obtained through the stress level heat model, and a toxicity heat color scale is preset according to the identification requirement of the soil Cd in synchronization; If the toxicity heat color scale does not cover any color value in the heat color value range, the enrichment sub-region is marked as a low stress area of the soil Cd; If the number of any color value in the heat color value range existing in the toxicity heat color scale is less than a preset number threshold, the enrichment sub-region is marked as a medium-low stress area of the soil Cd; If the number of any color value in the heat color value range existing in the toxicity heat color scale is greater than the preset number threshold, the enrichment sub-region is marked as a medium-high stress area of the soil Cd; If the toxicity heat color scale completely covers any color value in the heat color value range, the enrichment sub-region is marked as a high stress area of the soil Cd; and a global stress identification result of the soil Cd is generated.

6. The method for identifying and remediating Cd in soil based on mulberry root exudates and microorganisms according to claim 1, characterized in that, The S108 specifically includes the following steps: A repair requirement for the soil Cd in the target repair area is obtained, and a plurality of rhizosphere microbial species meeting the repair requirement and a resistance coefficient of each species to the stress level of the soil Cd are obtained through big data retrieval according to the global stress identification result; An expected stress identification result and an expected decline rate of the activity of the soil Cd are preset according to the repair requirement, a hash drift function between the global stress identification result and the expected stress identification result is calculated, and an actual decline amplitude of the current soil Cd pollution is determined according to the hash drift function; A Gaussian process algorithm is introduced, the actual decline amplitude is Gaussian modeled in the Gaussian process algorithm based on the resistance coefficient, a Gaussian regression model of the rhizosphere microbial species for the actual activity decline of the soil Cd is generated, and a resistance estimation prediction distribution of different rhizosphere microbial species reaching the actual decline amplitude is calculated through the Gaussian regression model; A balance degree between the standard deviation and the mean of the resistance estimation prediction distribution is analyzed by trade-off, and a limit decline rate of each rhizosphere microbial species to the activity stress of the soil Cd at different enrichment positions is constructed; If the limit decline rate is higher than the expected decline rate, the rhizosphere microbial species at the enrichment position corresponding to the limit decline rate is marked as a suitable species; if the limit decline rate is lower than the expected decline rate, the rhizosphere microbial species at the enrichment position corresponding to the limit decline rate is marked as a non-suitable species; The rhizosphere microbial community of the target repair area is planned and analyzed in combination with the suitable species and the non-suitable species, and an optimal community recruitment deployment of the rhizosphere microorganism is obtained.

7. The method according to claim 1, wherein the method is characterized by, The S110 specifically includes the following steps: Based on big data retrieval, the relative reduction characteristics of the rhizosphere microbial species in the optimal community recruitment deployment in the rhizosphere environment of the target repair area are obtained, and the relative oxidation characteristics of the root exudates in the local environment are retrieved. The dissolved oxygen flow model and the specific oxidation-reduction model of the soil environment are introduced, the dissolved oxygen flow model is defined as a macro scale, the specific oxidation-reduction model is defined as a micro scale, and a soil Cd environment map of a target remediation area is constructed based on a Cd enrichment layout map of the soil; The macro scale and the micro scale are cooperatively mapped to the soil Cd environment map based on a sampling criterion for converting macro conditions into micro constraints, to obtain a multi-scale coupling model of rhizosphere environment-dissolved oxygen interface oxidation-reduction dynamics; The interaction mechanism of root exudates-rhizosphere microorganisms is obtained, the coupling interaction strategy of the dissolved oxygen in the soil environment during oxidation-reduction is constructed according to the interaction mechanism, and the multi-scale coupling model is introduced to calculate the multi-scale state of the relative oxidation characteristics and the relative reduction characteristics under the condition of the limited supply of oxygen to rhizosphere microorganisms by root exudates, and a plurality of particle behavior trace representations are output; During the multi-scale calculation process, the effective response step length at which each particle behavior trace representation is formed is continuously recorded, and the residence index of each particle behavior trace representation on the enrichment staggered micro area is planned according to the span of the effective response step length; If the residence index is less than a preset residence index, the enrichment staggered micro area is ignored; if the residence index is greater than the preset residence index, a patch is created in the enrichment staggered micro area based on the representation boundary of the particle behavior trace characteristics in a stacking manner, and finally an oxidation-reduction patch pattern of the rhizosphere environment is obtained. Based on the interaction mechanism, the rhizosphere environment of the target remediation area is regulated according to the oxidation-reduction patch pattern, so as to adjust the Cd treatment ecology between the root exudates and the rhizosphere microorganisms, and realize the remediation of the soil Cd.

8. A soil Cd identification and remediation system based on mulberry root exudates and microorganisms, characterized in that, The system comprises a memory, a processor and a communication interface, the memory comprises a soil Cd identification and remediation method program based on mulberry root exudates and microorganisms, the communication interface is used for data connection and communication between the memory and the processor, and the soil Cd identification and remediation method program based on mulberry root exudates and microorganisms is executed by the processor to realize the soil Cd identification and remediation method steps of any one of claims 1-7.

Images