A method for dynamically adjusting the concentration of hydrogen-rich water

Abstract

The application provides a hydrogen-rich water concentration dynamic adjustment method, relates to the field of concentration self-adaptive adjustment, and comprises the following steps: S1, data acquisition; S2, data preprocessing; S3, key influence factor extraction; S4, hydrogen-rich water preliminary concentration acquisition; and S5, hydrogen-rich water optimized concentration acquisition. The dynamic concentration adjustment amount, the extreme environment adjustment amount and the implicit stress adjustment amount are respectively obtained through reinforcement learning dynamic optimization, extreme environment adaptation and implicit stress early warning, and the final concentration of the hydrogen-rich water is obtained. According to the growth state of the plant and environmental factors, the hydrogen-rich water concentration can be self-adaptively adjusted, the accuracy of the hydrogen-rich water concentration is ensured, and the problem that the hydrogen-rich water concentration is too high or too low to affect the normal growth of the plant is avoided.

Description

Technical Field

[0001] This invention relates to the field of concentration adaptive adjustment technology, and specifically to a method for dynamically adjusting the concentration of hydrogen-rich water. Background Technology

[0002] Adding hydrogen-rich water to irrigation water can regulate plant physiological metabolism by utilizing the antioxidant properties of hydrogen in the water. Specifically, hydrogen-rich water selectively neutralizes excess reactive oxygen species (ROS, such as hydroxyl radicals and superoxide anions) produced during plant growth, reducing oxidative damage to plants and protecting the function of organelles such as chloroplasts and mitochondria. At the same time, hydrogen promotes the synthesis of growth hormones such as auxin (IAA) and cytokinin (CTK) in plants and inhibits the accumulation of stress hormones such as abscisic acid (ABA), thereby effectively regulating plant cell division, growth, and nutrient synthesis. In addition, hydrogen-rich water can also regulate the soil microbial community, improve soil permeability, and promote root development, thereby enhancing the plant's ability to absorb water and nutrients and improving plant survival rate and growth.

[0003] While hydrogen-rich water can effectively enhance plant growth and increase agricultural yields, excessively high concentrations—meaning the concentration is mismatched with the plant's growth status—can disrupt the plant's redox balance (overly high concentrations of hydrogen-rich water excessively remove ROS, affecting cell division and gene expression, leading to root stagnation and hindering pollen germination). It can also inhibit plant growth, reduce photosynthetic efficiency, and cause nutrient imbalances. Currently, the concentration of hydrogen-rich water is primarily adjusted based on operator experience and extensive experimental data, making it highly susceptible to human error and prone to problems such as excessively high or low concentrations. (Excessive concentration leads to widespread plant wilting or yield reduction; insufficient concentration fails to achieve the desired effect, requiring repeated application, which is time-consuming and labor-intensive). This not only wastes hydrogen-rich water but also easily causes plant wilting or death, ultimately impacting crop yields and economic value. Summary of the Invention

[0004] To address the problems existing in the prior art, the present invention aims to provide a method for dynamically adjusting the concentration of hydrogen-rich water. This method belongs to the process control of crop cultivation. It adaptively adjusts the concentration of hydrogen-rich water according to the growth status of the plant and environmental factors, effectively ensuring the accuracy of the applied concentration of hydrogen-rich water and avoiding the problems of excessively high or low concentrations.

[0005] The objective of this invention is achieved through the following technical solution:

[0006] A method for dynamically adjusting the concentration of hydrogen-rich water, comprising:

[0007] Step S1: Data Acquisition: Based on the crop type and field conditions, deploy data acquisition sensors to collect plant growth trend data and field environmental factor data respectively.

[0008] Step S2, Data Preprocessing: The collected data undergoes noise removal, abnormal data detection and correction, and data normalization preprocessing.

[0009] Step S3: Extraction of key influencing factors: Explore plant growth trend data and field environmental factor data to identify key factors affecting crop growth;

[0010] Step S4: Obtaining the initial concentration of hydrogen-rich water: Input the data of key influencing factors into the fuzzy logic system, define the fuzzy subsets of each key influencing factor, and obtain the initial concentration of hydrogen-rich water.

[0011] Step S5, Obtaining the optimized concentration of hydrogen-rich water: The dynamic concentration adjustment amount, extreme environment adjustment amount, and latent stress adjustment amount are obtained through reinforcement learning dynamic optimization, extreme environment adaptation, and latent stress early warning, respectively, to obtain the final concentration of hydrogen-rich water.

[0012] Based on further optimization of the above scheme, the growth trend data includes chlorophyll fluorescence, canopy multispectral data, stem microdeformation, leaf infrared temperature, root growth, etc.; field environmental factor data includes light intensity, air temperature and humidity, soil temperature and humidity, CO2 concentration, soil pH, etc.

[0013] Based on further optimization of the above scheme, the noise removal in step S2 is to remove random noise using the Kalman filter algorithm.

[0014] Based on further optimization of the above scheme, the abnormal data detection and correction in step S2 specifically includes:

[0015] For input feature data X =[ x 1, x 2,…, x n The Isolation Forest algorithm was used to detect anomalies in the data after removing random noise, and anomaly scores were obtained.

[0016] ;

[0017] In the formula: h(X) Indicates sample X Path length in an isolated tree; E[h(X)] Indicates the expected path length; Represents Euler's constant;

[0018] Preset abnormal score threshold SX d ,like If the value is 0, then the corresponding sample is determined to be an outlier.

[0019] For outliers, Huber robust regression is used for correction. The Huber loss function is:

[0020] ;

[0021] ;

[0022] In the formula: The parameters are constants; by optimizing the model parameters Improve tolerance to extreme environments.

[0023] Based on further optimization of the above scheme, in step S3, key factors are extracted using a stress-metabolism coupling network, a dynamic time-delay network, a multi-scale feature fusion network, and a stress state adaptive network, respectively. Specifically:

[0024] The preprocessed data variables from step S2 are sequentially input into four networks: a stress-metabolism coupling network, a dynamic time-delay network, a multi-scale feature fusion network, and a stress state adaptive network. If a variable is identified as a key element by one of the four network structures, the corresponding variable is incremented by 1 in the key factor determination. If the total value of a variable is not less than 3, the variable is identified as a key influencing factor. If the total value of a variable is 2, it is considered a secondary influencing factor (when no key influencing factor is obtained, the secondary influencing factor is used as the key influencing factor).

[0025] Based on further optimization of the above scheme, the specific method for identifying key elements in the stress-metabolism coupling network is as follows:

[0026] First, the preprocessed data variables from step S2 X’ Based on the effects of ROS (reactive oxygen species) on plants, the plant is divided into an environmental layer (CeX), a physiological layer (CeY), and an effect layer (CeZ). The environmental layer (CeX) corresponds to environmental indicators that directly induce ROS generation (i.e., ROS generation drivers, including data such as light intensity, air temperature and humidity, soil temperature and humidity, CO2 concentration, and soil pH). The physiological layer (CeY) reflects indicators of ROS damage to plant function (i.e., ROS damage targets, including data such as chlorophyll fluorescence and leaf infrared temperature). The effect layer (CeZ) reflects indicators of the final traits of plants affected by ROS (i.e., ROS regulation results, including data such as root growth, canopy multispectral data, and stem microdeformation).

[0027] Then, two interlayer unidirectional influence pathways were established, including the first pathway CeX→CeY (the influence of environmental factors on physiological state) and the second pathway CeY→CeZ (the influence of physiological state on growth effect).

[0028] Then, the path coefficients of the two paths are obtained using the partial least squares path model (PLS-PM): , ,

[0029] Finally, preset path coefficient threshold ,like Then the stress-metabolism coupling network will identify the corresponding variables as key elements.

[0030] Based on further optimization of the above scheme, the specific method for identifying key elements using the dynamic time-delay network is as follows:

[0031] First, obtain the preprocessed data variables from step S2. X’(t) (This indicates that the characteristic variables were collected at time t and preprocessed in step S2) and concentration requirements. In time lag Mutual information under TDMI for:,

[0032] ;

[0033] In the formula: P represents the probability distribution; Representing variables X’ In time t Pick a Value and concentration requirement B in time Pick b The probability of the value; Representing variables X’ In time t Pick a The probability of the value This indicates that the concentration requirement B is in time. Pick b The probability of the value;

[0034] Plant growth trend data and field environmental factor data were obtained at time lags. The TDMI is set below, with a preset TDMI peak threshold. max If the mutual information TDMI of a certain variable is greater than TDMI max Then the corresponding variable is determined to be a key element by the dynamic time-delay network.

[0035] Based on further optimization of the above scheme, the specific method for identifying key elements using the multi-scale feature fusion network is as follows:

[0036] For time series variables X’(t)A three-level wavelet packet decomposition was performed (sampling frequency was once per hour, with a cumulative sampling period of 4 weeks to ensure data coverage of the "hour-day-week" scale); the db4 wavelet basis function was used for decomposition, dividing the time series variables into three levels: instantaneous (hourly), short-term (daily), and long-term (weekly), with a frequency band of 8; the 8 frequencies were sorted from high to low. j =1~2 represents high frequency, corresponding to the hour pole. j =3 to 6 represent medium frequency, corresponding to short-term day-level frequencies. j =7~8 are low frequencies, corresponding to long-term cycles), obtaining the energy values ​​of eight frequency bands. E i,j :

[0037] ;

[0038] In the formula: x i,j (t) Indicates the first i The first variable j Sub-signals in each frequency band; t 0、 t n These represent the start and end times of the acquisition cycle corresponding to the sub-signal, respectively.

[0039] The sub-signals of each frequency band are input as independent features into a random forest to obtain the Gini importance of each sub-signal; weights are assigned to different frequency bands based on the metabolic cycle caused by hydrogen-rich water. To obtain the overall importance of the corresponding variables:

[0040] ;

[0041] In the formula: Gini i,j Indicates the first i The first variable j The importance of the Gini coefficient in frequency band characteristics;

[0042] Preset comprehensive importance threshold I d If the overall importance of a variable is greater than the overall importance threshold, then the corresponding variable is determined to be a key element by the multi-scale feature fusion network.

[0043] Based on further optimization of the above scheme, the specific method for identifying key elements by the stress state adaptive network is as follows:

[0044] Introducing the Stress Index SI:

[0045] ;

[0046] In the formula: Tleaf Indicates the infrared temperature of the blade; T opt This indicates the optimal temperature for the leaves; T max Indicates the critical high temperature; Indicates soil volumetric moisture content. Indicates field holding capacity; K 4. K 5 represents the corresponding weight coefficients;

[0047] Preset stress index threshold SI d :

[0048] when When the time is right, it indicates that the plant's metabolism is basically normal. Variables related to photosynthesis (such as canopy multispectral, light intensity, CO2 concentration, chlorophyll fluorescence, root growth, etc.) are used as input variables of the Bayesian network. The Bayesian network is used to calculate the posterior probability of the variables and the concentration of hydrogen-rich water. The variables are arranged from largest to smallest and the top three are taken as key elements.

[0049] when When the plant's photosynthesis / respiration is inhibited, stress-related variables (such as the degree of stem microdeformation, leaf infrared temperature, soil temperature and humidity, air temperature and humidity, soil pH, etc.) are used as input variables of the Bayesian network. The posterior probability of the variables and the concentration of hydrogen-rich water is calculated through the Bayesian network, and the variables are arranged from largest to smallest. The top three variables are taken as key elements.

[0050] A Bayesian network is a probabilistic graphical model that outputs posterior probabilities based on the causal relationship between "field environmental factors / plant growth trends" and "hydrogen-rich water concentration requirements".

[0051] Based on further optimization of the above scheme, step S4 specifically includes:

[0052] Define the membership function of the fuzzy subset of key influencing factors, including a small number of memberships. U dr (S) Suitable membership degree U su (S) Excessive membership U we (S) :

[0053] ;

[0054] In the formula: S Indicates the current measured values ​​of key influencing factors; S min This indicates the critical value when the key influencing factor is too low.S a Indicates the appropriate central value of key influencing factors; S max This indicates the critical value when key influencing factors are too high.

[0055] A fuzzy rule base is constructed by integrating multi-factor correlation and expert practical experience, and the centroid method is used to defuzzify and obtain the initial concentration.

[0056] ;

[0057] In the formula: U i Indicates the first i The membership degree of a rule; C i Indicates the first i The reference concentration for this rule; m This indicates the number of rules that are active.

[0058] Based on the above scheme, the reinforcement learning dynamic optimization specifically includes:

[0059] Using the vectors of key influencing factors (initial concentration and real-time monitoring) as the state space Y, and the concentration adjustment of hydrogen-rich water as the action space A, a reward function is constructed as follows:

[0060] ;

[0061] In the formula: This indicates the change in photosynthetic efficiency; This indicates the change in root growth. This indicates the adjustment range of the hydrogen-rich water concentration; These represent the corresponding weight coefficients;

[0062] Update the Q value using the Q-Learning algorithm:

[0063] ;

[0064] In the formula: This represents the learning rate, typically between 0.1 and 0.3. This represents the discount factor, typically ranging from 0.7 to 0.9. Q(Y t , A t ) Represents the state at the current time t. Y t -action A t The corresponding Q value; Indicates the next state Y t+1Next, iterate through all possible actions. a Choose the one with the largest Q value;

[0065] When the change in Q value is less than a preset threshold, Q value converges; after Q value convergence, for any current state... Y t Choose the action that maximizes the Q value. A * :

[0066] ;

[0067] action A * This refers to the adjustment amount for dynamic optimization of hydrogen-rich water concentration in the current state using reinforcement learning. .

[0068] Based on further optimization of the above scheme, the extreme environment adaptation specifically includes:

[0069] Collect historical extreme environment data to construct a meta-task set:

[0070] ;

[0071] In the formula: This represents a small amount of training data for the k-th extreme scenario; This represents the test data for the k-th extreme scenario;

[0072] Use MAML (Model-Agnostic Meta-Learning) to pre-train general initial parameters :

[0073] For each metatask k ,use Data is updated within the task:

[0074] ;

[0075] In the formula: Indicates the learning rate of the inner loop; This represents the training loss for task k; This represents the gradient operator, i.e., with respect to parameters. Find the partial derivative;

[0076] use Data is updated meta-level and initial parameters are optimized:

[0077] ;

[0078] In the formula: This represents the outer loop learning rate;

[0079] The novelty of environmental states is tested using isolated forests, and a novelty score is output.

[0080] A preset novelty score threshold is established: if the novelty score exceeds the threshold, it is considered a new extreme environment; in the new extreme environment, real-time data is collected for meta-learning updates.

[0081] ;

[0082] In the formula: The loss function represents the loss function for the new environmental data;

[0083] ;

[0084] In the formula: Indicates the model uses the current parameters And based on the characteristics of the new environment X new Predicted demand for hydrogen-rich water concentration; y new This represents the optimal concentration requirement actually obtained in the corresponding extreme environment (obtained through historical experience data). N new Indicates the number of samples in the new environmental data;

[0085] Use the adapted parameters Combined with reinforcement learning dynamic optimization to achieve Q ( Y t ,A t ; The biggest action A 2 * :

[0086] ;

[0087] action A 2 * That is, the adjustment amount of hydrogen-rich water concentration to adapt to extreme environments. .

[0088] Based on further optimization of the above scheme, the implicit stress early warning specifically refers to:

[0089] Microscopic features of early latent stress were extracted from chlorophyll fluorescence and Raman spectroscopy data to construct a stress feature vector. Z x ;

[0090] A variational autoencoder (VAE) is used to learn the feature distribution under normal conditions to obtain early latent stress: the VAE encoder will learn the feature distribution under normal conditions. Z xMapping to latent space Output mean and variance; the VAE decoder will output the latent spatial variables. Reconstructed as features (in, and Z x (For vehicles with the same dimensions), the loss function for VAE is:

[0091] ;

[0092] In the formula: This represents the "conditional probability distribution" defined by the decoder, which describes the reconstructed features. The probability distribution of the original input feature Z; E[·] represents the expected value of the reconstruction loss; This represents the approximate posterior distribution of the encoder output; represents the prior distribution of the latent variable; KL[·] represents the KL divergence; This represents the balance coefficient (usually 1).

[0093] Calculate the reconstruction error REW:

[0094] ;

[0095] In the formula: Representing the eigenvector Z x The dimension;

[0096] Preset error threshold If the reconstruction error REW is greater than the error threshold, then the feature is determined. Z x This is early-stage, latent stress;

[0097] For characteristics identified as early latent stress Z x Calculate its relationship with the historical coercion database { Z 1, Z 2,…, Z n The Euclidean distance of}

[0098] ;

[0099] Choose the one with the smallest distance. h 1 sample, to obtain the initial adjustment amount:

[0100] ;

[0101] In the formula: C h,i Representation and Features Z x The smallest distancei The historical adjustment amount corresponding to each historical sample;

[0102] Combining reinforcement learning dynamic optimization, that is, the action space A in reinforcement learning dynamic optimization is... If the fine-tuning interval is centered, the reward function is updated as follows:

[0103] ;

[0104] In the formula: This indicates the change in photosynthetic efficiency; This indicates the change in root growth. This indicates the adjustment range of the hydrogen-rich water concentration; These represent the corresponding weight coefficients;

[0105] Latent stress adjustment for:

[0106] ;

[0107] In the formula: It represents the gradient of the Q-value with respect to the reward function, reflecting the impact of changes in the adjustment amount on the reward; Indicates the fine-tuning coefficient;

[0108] The specific method for obtaining it is as follows:

[0109] First, the Q-value is defined as "the long-term cumulative reward for performing adjustment action A in state Y," and its expression is the classic form of reinforcement learning, namely...

[0110] ;

[0111] but:

[0112] ;

[0113] In the formula: This represents the discount factor.

[0114] The following are the technical effects of the present invention:

[0115] This invention identifies key elements through four networks: a stress-metabolism coupling network, a dynamic time-delay network, a multi-scale feature fusion network, and a stress state adaptive network. This allows for the output of key influencing factors for hydrogen-rich water. The identification process not only anchors the core physiological mechanisms, supplements the dynamics of the time dimension, covers the completeness of the time scale, and adapts to the variability of plant states, but also avoids problems such as detachment from physiological mechanisms, time mismatch, and inability to adapt to complex and changing field environments caused by outputting key influencing factors through a single network. This overcomes the limitations of traditional static analysis, ultimately identifying key influencing factors strongly correlated with hydrogen-rich water concentration requirements, achieving precise monitoring and regulation, eliminating redundant factor detection, and improving monitoring and control efficiency. Meanwhile, this invention dynamically adjusts the initial concentration of hydrogen-rich water output by the fuzzy logic system through reinforcement learning dynamic optimization, extreme environment adaptation, and implicit stress early warning. This effectively solves problems such as static rules being unable to adapt to real-time changes, conventional rules being unable to cover extreme scenarios, and explicit rules being unable to capture early implicit stress damage. It ensures that the hydrogen-rich water concentration not only matches the real-time state but also covers extreme situations and can predict implicit risks, achieving real-time and precise adjustment of the hydrogen-rich water concentration and avoiding problems of excessively low or high hydrogen-rich water concentration.

[0116] This invention can adjust the concentration of hydrogen-rich water in real time, continuously and accurately according to the growth status of plants and environmental factors, thereby improving the accuracy, adaptability and stability of hydrogen-rich water irrigation, thus meeting the needs of complex and ever-changing field scenarios, and avoiding problems such as reduced yield and death of plants due to excessively high hydrogen-rich water concentration or poor irrigation effect and increased cost due to excessively low hydrogen-rich water concentration. Attached Figure Description

[0117] Figure 1 This is a structural block diagram illustrating the dynamic adjustment of hydrogen-rich water concentration in an embodiment of the present invention. Detailed Implementation

[0118] The technical solutions in the embodiments of the present invention will be clearly and completely described below. In the following description, specific details such as specific system structures and technologies are presented for illustration rather than limitation, so as to provide a thorough understanding of the embodiments of the present invention.

[0119] Example 1:

[0120] A method for dynamically adjusting the concentration of hydrogen-rich water, comprising:

[0121] Step S1, Data Acquisition: Based on the crop type and field conditions, deploy acquisition sensors to collect plant growth trend data and field environmental factor data respectively; growth trend data includes chlorophyll fluorescence, canopy multispectral data, stem microdeformation degree, leaf infrared temperature, root growth, etc.; field environmental factor data includes light intensity, air temperature and humidity, soil temperature and humidity, CO2 concentration, soil pH value, etc.

[0122] Step S2, Data Preprocessing: The collected data undergoes preprocessing including Kalman filtering to remove random noise (using existing conventional Kalman filtering methods), outlier detection and correction, and data normalization (using existing conventional normalization methods). Specifically, outlier detection and correction are performed as follows:

[0123] For input feature data X =[ x 1, x 2,…, x n The Isolation Forest algorithm was used to detect anomalies in the data after removing random noise, and anomaly scores were obtained.

[0124] ;

[0125] In the formula: h(X) Indicates sample X Path length in an isolated tree; E[h(X)] Indicates the expected path length; Represents Euler's constant;

[0126] Preset abnormal score threshold SX d (Usually 0.7), if If the value is 0, then the corresponding sample is determined to be an outlier.

[0127] For outliers, Huber robust regression is used for correction. The Huber loss function is:

[0128] ;

[0129] ;

[0130] In the formula: This is a constant, typically taken as 1.35; the model parameters are optimized. Improve tolerance to extreme environments.

[0131] Step S3: Extraction of key influencing factors: Explore plant growth trend data and field environmental factor data to identify key factors affecting crop growth;

[0132] Key factors were extracted using stress-metabolism coupling networks, dynamic time-delay networks, multi-scale feature fusion networks, and stress state adaptive networks, respectively. Specifically:

[0133] The preprocessed data variables from step S2 are sequentially input into four networks: a stress-metabolism coupling network, a dynamic time-delay network, a multi-scale feature fusion network, and a stress state adaptive network. If a variable (i.e., any variable in the growth trend data and field environmental factor data) is identified as a key element by one of the four network structures, the corresponding variable's value is incremented by 1 in the key factor determination. If the total value of a variable is not less than 3, the variable is considered a key influencing factor. If the total value of a variable is 2, it is considered a secondary influencing factor (in the absence of a key influencing factor, the secondary influencing factor is used as the key influencing factor). For example, if leaf infrared temperature is identified as a key element by the stress-metabolism coupling network, the dynamic time-delay network, and the stress state adaptive network, but not by the multi-scale feature fusion network, then the leaf infrared temperature value is 3, making it a key influencing factor.

[0134] The specific method for identifying key elements in stress-metabolism coupling networks is as follows:

[0135] First, the preprocessed data variables from step S2 X’ Based on the effects of ROS (reactive oxygen species) on plants, the plant is divided into an environmental layer (CeX), a physiological layer (CeY), and an effect layer (CeZ). The environmental layer (CeX) corresponds to environmental indicators that directly induce ROS generation (i.e., ROS generation drivers, including data such as light intensity, air temperature and humidity, soil temperature and humidity, CO2 concentration, and soil pH). The physiological layer (CeY) reflects indicators of ROS damage to plant function (i.e., ROS damage targets, including data such as chlorophyll fluorescence and leaf infrared temperature). The effect layer (CeZ) reflects indicators of the final traits of plants affected by ROS (i.e., ROS regulation results, including data such as root growth, canopy multispectral data, and stem microdeformation).

[0136] Then, two interlayer unidirectional influence pathways were established, including the first pathway CeX→CeY (the influence of environmental factors on physiological state) and the second pathway CeY→CeZ (the influence of physiological state on growth effect).

[0137] Then, the path coefficients of the two paths are obtained using the partial least squares path model (PLS-PM): , :

[0138] If the stress-metabolism coupling network (i.e., the partial least squares path model) is a mapping of "latent variables to manifest variables", then the sign of the latent variables in the environmental layer is... The corresponding explicit variable is the light intensity after preprocessing in step S2. air temperature and humidity Soil temperature and humidity CO2 concentration Soil pH Etc., that is, the explicit variables of the environment layer are ( j =1, 2, 3, 4, 5); similarly, the symbols for latent variables in the physiological layer are... The corresponding dominant variable is chlorophyll fluorescence. Blade infrared temperature etc., that is ( j =1, 2); the latent variables in the effect layer are , and the corresponding manifest variables are root growth potential. Canopy Multispectral Degree of slight deformation of stem etc., that is ( j =1, 2, 3);

[0139] For example: the latent variable score for the environmental layer is:

[0140] ;

[0141] In the formula: In the CeX environment layer, the first j The first environment variable, the first k The weight coefficients for the next iteration;

[0142] Simultaneous generation of latent variable scores for the physiological layer:

[0143] ;

[0144] PLS-PM iterates through a loop of "external weights → internal weights → path coefficients" until the weight change is less than a threshold (e.g., 10). -4 That is, the weight difference between two adjacent iterations. Stop the iteration and output the final path coefficients; the path coefficients represent the strength and direction of the influence between latent variables, and their essence is the "standardized regression coefficient of the target latent variable on the source latent variable".

[0145] in:

[0146] External weights of the environment layer:

[0147] ;

[0148] Internal weights:

[0149] ;

[0150] In the formula: Cov() Represents covariance;Var() Indicates variance;

[0151] because Standardization has been performed, and the variance of the standardized variables can be simplified to 1, i.e.:

[0152] ;

[0153] ;

[0154] In the formula: n represents the sample size (i.e., the total number of plants collected in the field, such as n=100), and X and Y represent the corresponding variables (such as the latent variable scores of the environmental layer). Physiological layer latent variable score wait); X i Representing variables X The i Observed values ​​of plant samples; Representing variables X The sample mean;

[0155] Similarly, the second path coefficients are obtained. ;

[0156] Finally, preset path coefficient threshold (Usually 0.5), if In this case, the stress-metabolism coupling network identifies the corresponding variables as key elements; for example: (That is, in this path, CeX is the input layer variable and CeY is the output layer variable), and the corresponding variable is the data in the environment layer CeX, then the variable in the environment layer is determined to be a key element; if (That is, CeY is the input layer variable and CeZ is the output layer variable in this path), and the corresponding variable is the data in the physiological layer CeY. Then the variable in the physiological layer is determined to be the key element.

[0157] The specific method for identifying key elements using dynamic time-delay networks is as follows:

[0158] First, obtain the preprocessed data variables from step S2. X’(t) (This indicates that the characteristic variables were collected at time t and preprocessed in step S2) and concentration requirements. In time lag Mutual information under TDMI for:,

[0159] ;

[0160] In the formula: P represents the probability distribution; Representing variables X’ In time t Pick aValue and concentration requirement B in time Pick b The probability of the value; Representing variables X’ In time t Pick a The probability of the value This indicates that the concentration requirement B is in time. Pick b The probability of the value; where concentration requirement B is the target concentration of hydrogen-rich water required by the plant under the corresponding growth state; (the probability can be obtained using the frequency probability approximation method, for example: based on the time series data collected in step S1: the growth trend / environmental factor variable X is collected once per hour for a cumulative period of 4 weeks; the time series data of concentration requirement B comes from the matching records of "key factors-optimal concentration" in the historical database, and the continuous variables after preprocessing in step S2 are...) X’ The concentration requirement B was divided into 5–10 equally spaced intervals, with the number of intervals adjusted according to the data distribution density; for the variable… X’ The result is obtained by counting the frequency of the value 'a' within a certain interval in hour t and dividing by the total amount of data. For concentration requirement B, its statistical value is... Take the value within the interval within the hour b The frequency percentage, that is, Statistical variables X’ Take value a within hour t, and the concentration requirement B is... Take values ​​within an hour b The frequency of synchronous occurrences divided by the total amount of data yields the result. (Data interval partitioning and frequency statistics can be performed using Python's pandas library, and probability calculations can be performed using the scipy.stats library.)

[0161] Plant growth trend data and field environmental factor data were obtained at time lags. TDMI under the range of 0–48 h, with a preset TDMI peak threshold. max (Typically 0.8 bits), if the mutual information TDMI of a certain variable is greater than TDMI max Then the corresponding variable is determined to be a key element by the dynamic time-delay network.

[0162] The specific method for identifying key elements using multi-scale feature fusion networks is as follows:

[0163] For time series variables X’(t)A three-level wavelet packet decomposition was performed (sampling frequency was once per hour, with a cumulative sampling period of 4 weeks to ensure data coverage across hourly, daily, and weekly scales; the db4 wavelet basis function was used for decomposition, dividing the time-series variables into three levels: instantaneous (hourly), short-term (daily), and long-term (weekly). The first level, instantaneous response: sudden increase in leaf temperature (direct sunlight) → hydrogen-rich water needs to quickly alleviate transpiration stress, with 2 frequency bands; the second level, short-term response: slight stem deformation (water deficit) → hydrogen-rich water needs to continuously regulate turgor pressure, with 4 frequency bands; the third level, long-term response: root growth (nutrient absorption) → hydrogen-rich water needs to promote metabolic repair, with 8 frequency bands; the 8 frequencies were sorted from high to low. j =1~2 represents high frequency, corresponding to the hour pole. j =3 to 6 represent medium frequency, corresponding to short-term day-level frequencies. j =7~8 are low frequencies, corresponding to long-term cycles), obtaining the energy values ​​of eight frequency bands. E i,j :

[0164] ;

[0165] In the formula: x i,j (t) Indicates the first i The first variable j Sub-signals in each frequency band; t 0、 t n These represent the start and end times of the acquisition cycle corresponding to the sub-signal, respectively.

[0166] The sub-signals of each frequency band are used as independent features and input into a random forest (the random forest model chosen is RandomForestRegressor from the Python scikit-learn library, using 8 frequency band sub-signals as 8 independent feature inputs and "hydrogen-rich water concentration requirement" as the label). A label vector is formed by matching the probability of "8 sub-signal states → optimal irrigation hydrogen-rich water concentration" in field experiments. The preprocessed feature matrix is ​​then used as the input. X’ The label vectors are input into a random forest model for training. After the random forest is trained, the Gini importance of each of the eight frequency band sub-signals is obtained directly through the model's built-in feature_importance_ attribute. Based on the metabolic cycle caused by hydrogen-rich water, weights are assigned to different frequency bands. (Generally high frequency) Mid-frequency low frequency ), to obtain the overall importance of the corresponding variables:

[0167] ;

[0168] In the formula: Gini i,j Indicates the first i The first variable j The importance of the Gini coefficient in frequency band characteristics;

[0169] Preset comprehensive importance threshold I d If the overall importance of a variable is greater than the overall importance threshold, then the corresponding variable is determined to be a key element by the multi-scale feature fusion network.

[0170] The specific method for identifying key elements using a stress-state adaptive network is as follows:

[0171] Introducing the Stress Index SI:

[0172] ;

[0173] In the formula: T leaf Indicates the infrared temperature of the blade; T opt This indicates the optimal temperature for the leaves; T max Indicates the critical high temperature; Indicates soil volumetric moisture content. Indicates field holding capacity; K 4. K 5 represents the corresponding weight coefficients;

[0174] Preset stress index threshold SI d (Usually 0.3):

[0175] when When the plant's metabolism is basically normal, variables related to photosynthesis (such as canopy multispectral density, light intensity, CO2 concentration, chlorophyll fluorescence, root growth, etc.) are used as input variables of the Bayesian network. The Bayesian network calculates the posterior probability of the variables and the concentration of hydrogen-rich water, and the variables are arranged from largest to smallest. The top three variables are selected as key elements (e.g., canopy multispectral density, light intensity, CO2 concentration).

[0176] when When the plant's photosynthesis / respiration is inhibited, stress-related variables (such as the degree of stem microdeformation, leaf infrared temperature, soil temperature and humidity, air temperature and humidity, soil pH, etc.) are used as input variables of a Bayesian network. The Bayesian network calculates the posterior probability of the variables and the concentration of hydrogen-rich water, and the variables are arranged from largest to smallest. The top three variables are selected as key elements (e.g., the degree of stem microdeformation, leaf infrared temperature, and soil temperature and humidity).

[0177] A Bayesian network is a probabilistic graphical model that outputs a posterior probability based on the causal relationship between "field environmental factors / plant growth trends" and "hydrogen-rich water concentration requirements," specifically:

[0178] First, growth trend data and corresponding variables of field environmental factors are used as input variables for the parent node, and hydrogen-rich water concentration requirements are used as output variables for the child nodes. Based on the mechanism of action of hydrogen-rich water and past empirical experimental data, a directed acyclic graph (DAG) is drawn (the structure can be defined using the BayesianNetwork class of the Python pgmpy library, for example: light intensity ↑ → concentration requirement ↓, soil moisture content ↓ → concentration requirement ↑, chlorophyll fluorescence ↓ → concentration requirement ↑, stem microdeformation degree ↑ → concentration requirement ↑). Then, the conditional probability of the child node (concentration requirement) under different parent node states is calculated to form a conditional probability table (CPT, which can be implemented using the MaximumLikelihoodEstimator class of the pgmpy library). Finally, the measured values ​​of variables at the current moment in the field are collected and converted into network-recognizable states (consistent with the discretization / continuation rules in the data preparation stage). Based on Bayes' theorem, the posterior probability of each variable and the current concentration requirement is inferred (using the VariableElimination inferencer of the pgmpy library).

[0179] ;

[0180] In the formula: These represent the prior probability of the variable (obtained by the frequency percentage of each state of the variable occurring naturally in the field), the conditional probability of the concentration requirement (obtained by statistically analyzing the corresponding frequency percentages based on historical matching data of "variable state - optimal concentration requirement"), and the marginal probability of the concentration requirement (obtained by directly statistically analyzing the frequency percentages of each state of the concentration requirement in historical data).

[0181] Step S4: Obtaining the initial concentration of hydrogen-rich water: Input the data of key influencing factors into the fuzzy logic system, define the fuzzy subsets of each key influencing factor, and obtain the initial concentration of hydrogen-rich water.

[0182] Specifically:

[0183] Define the membership function (range [0,1]) of the fuzzy subset of key influencing factors, including a small number of memberships. U dr (S) Suitable membership degree U su (S) Excessive membership U we (S) :

[0184] ;

[0185] In the formula: S Indicates the current measured values ​​of key influencing factors; S min This indicates the critical value when the key influencing factor is too low. S a Indicates the appropriate central value of key influencing factors; S max This indicates the critical value when key influencing factors are too high.

[0186] A fuzzy rule base is constructed by integrating multi-factor correlations and expert practical experience. (When constructing the fuzzy rule base, each rule is a combination of the membership degrees of multiple factors; for example, a rule could be: if the light intensity is "excessive", the soil moisture is "slight", and the chlorophyll fluorescence is "suitable", then the reference concentration of hydrogen-rich water is...) C i The initial concentration was 1.5 mg / L, and the centroid method was used to defuzzify the sample to obtain the preliminary concentration.

[0187] ;

[0188] In the formula: U i Indicates the first i The membership degree of a rule; C i Indicates the first i The reference concentration for this rule; m This indicates the number of rules that are active.

[0189] Step S5, Obtaining the optimized concentration of hydrogen-rich water: The dynamic concentration adjustment amount, extreme environment adjustment amount, and latent stress adjustment amount are obtained through reinforcement learning dynamic optimization, extreme environment adaptation, and latent stress early warning, respectively, to obtain the final concentration of hydrogen-rich water.

[0190] Specifically, reinforcement learning dynamic optimization includes:

[0191] The initial concentration and the vector of key influencing factors in real-time monitoring are used as the state space Y (state Y = {initial concentration of hydrogen-rich water}). C st +Key influencing factor vector}), using the concentration adjustment of hydrogen-rich water as the action space A, construct the reward function:

[0192] ;

[0193] In the formula: This indicates the change in photosynthetic efficiency; This indicates the change in root growth. This indicates the adjustment range of the hydrogen-rich water concentration; These represent the corresponding weight coefficients;

[0194] The Q-value is updated using the Q-Learning algorithm (Q-Learning is a classic value-based, model-free algorithm in reinforcement learning. Its core objective is to learn a "state-action" value function, or Q-function, through trial and error between the agent and the environment, thereby guiding the agent to choose the optimal action in any state to maximize long-term cumulative reward).

[0195] ;

[0196] In the formula: This represents the learning rate, typically ranging from 0.1 to 0.3. This represents the discount factor, typically ranging from 0.7 to 0.9. Q(Y t , A t ) Represents the state at the current time t. Y t -action A t The corresponding Q value (to the left of the arrow) Q(Y t ,A t ) The updated, more optimal value estimate, i.e., the "corrected new understanding"; the right side of the arrow Q(Y t ,A t ) This represents the value estimate of action A performed in state Y before the update, i.e., the old knowledge; both the left and right sides of the formula show... Q(Y t ,A t ) This is because the Q-Learning algorithm uses an "incremental update" logic, which means it does not directly replace the old Q value, but rather makes corrections based on the original estimate and new experience. Indicates the next state Y t+1 Next, iterate through all possible actions. a Select the one with the largest Q-value (a core feature of the Q-Learning algorithm, demonstrated through offline learning, representing the new state). Y t+1 (the optimal expected reward);

[0197] When the change in Q value is less than a preset threshold (typically 0.01), Q value converges; after Q value convergence, for any current state... Y tChoose the action that maximizes the Q value. A * :

[0198] ;

[0199] action A * This refers to the adjustment amount for dynamic optimization of hydrogen-rich water concentration in the current state using reinforcement learning. .

[0200] Extreme environment adaptation specifically refers to:

[0201] Collect historical extreme environment data to construct a meta-task set:

[0202] ;

[0203] In the formula: This represents a small amount of training data for the k-th extreme scenario; This represents the test data for the k-th extreme scenario;

[0204] Use MAML (Model-Agnostic Meta-Learning) to pre-train general initial parameters :

[0205] For each metatask k ,use Data is updated within the task:

[0206] ;

[0207] In the formula: Indicates the learning rate of the inner loop; This represents the training loss for task k; This represents the gradient operator, i.e., with respect to parameters. Find the partial derivative;

[0208] use Data is updated meta-level and initial parameters are optimized:

[0209] ;

[0210] In the formula: This represents the outer loop learning rate;

[0211] The novelty of environmental states is assessed using isolated forests, and a novelty score is output, which includes:

[0212] Collect at least three types of historical extreme scenario data (such as high temperature and drought, heavy rain and high humidity, extreme pH stress, etc.), with no fewer than 50 data samples for each scenario, and remove completely duplicate samples to serve as the training set; collect real-time data of the current field environment as the test set, with its feature dimensions completely consistent with the training set (e.g., if the training set contains 5 features, the test set cannot be increased or decreased), and the collection frequency of the real-time data is consistent with the training set (e.g., 1 data sample per hour); then, use the data preprocessing method of step S2 (Kalman filtering algorithm and Huber robust regression to preprocess the collected data), and map the features of the training set and test set to the [0,1] interval; use the IsolationForest class of the Python scikit-learn library as the isolated forest model, and train the model through the training set; the novelty score of the isolated forest is based on the "average path length of the sample in the isolated tree" transformation, and the score range is [0,1]. For example, for real-time test samples... x test The isolated forest model was used to calculate its performance. N tree Average path length in an isolated tree Thus, a novelty score is obtained. S xyx :

[0213] ;

[0214] In the formula: E[h(x)] This represents the average expected path length of the training set samples (automatically calculated by the model and related to the sample size and number of features).

[0215] A preset novelty score threshold (typically 0.9) is set: if the novelty score exceeds the threshold, it is considered a new extreme environment; in the new extreme environment, real-time data is collected for meta-learning updates.

[0216] ;

[0217] In the formula: The loss function represents the data from the new environment;

[0218] ;

[0219] In the formula: Indicates the model uses the current parameters And based on the characteristics of the new environment X new Predicted demand for hydrogen-rich water concentration; y new This represents the optimal concentration requirement actually obtained in the corresponding extreme environment (obtained through historical experience data). The model parameters are: At that time, in the new extreme environment, the first i The predicted output value for each sample; N new Indicates the number of samples in the new environmental data;

[0220] Use the adapted parameters Combined with reinforcement learning dynamic optimization to achieve Q ( Y t ,A t ; The biggest action A 2 * (That is, when the change in Q value is less than a preset threshold of 0.01, Q value converges; after Q value convergence, for any current state) Y t Choose the action that maximizes the Q value. A 2 * (Similar to the reinforcement learning dynamic optimization above)

[0221] ;

[0222] action A 2 * That is, the adjustment amount of hydrogen-rich water concentration to adapt to extreme environments. .

[0223] The latent stress early warning mechanism is a mechanism for identifying potential stresses (such as mild drought, latent nutrient imbalance, and micro-abiotic stress) in the early stages of dynamic regulation of hydrogen-rich water concentration, before any visible damage is observed in the plant. Its core is to accurately determine the presence and severity of stress by capturing abnormal changes in the plant's microscopic physiological characteristics, before the stress manifests as visible symptoms (such as wilting, yellowing, or growth stagnation). This provides a "predictive basis" for adjusting the hydrogen-rich water concentration, preventing the stress from developing into visible damage. The latent stress early warning mechanism is the perception mechanism for "early detection of potential stresses," while the hydrogen-rich water adjustment amount is the execution parameter for "precisely correcting the concentration for potential stresses." The two are causally linked.

[0224] Microscopic features of early latent stress were extracted from chlorophyll fluorescence and Raman spectroscopy data to construct a stress feature vector. Z x ;

[0225] A variational autoencoder (VAE) is used to learn the feature distribution under normal conditions to obtain early latent stress: the VAE encoder will learn the feature distribution under normal conditions. Z x Mapping to latent space Output mean and variance; the VAE decoder will output the latent spatial variables. Reconstructed as features (in, and Z x (For vehicles with the same dimensions), the loss function for VAE is:

[0226] ;

[0227] In the formula: This represents the "conditional probability distribution" defined by the decoder, which describes the reconstructed features. The probability distribution of the original input feature Z; E[·] represents the expected value of the reconstruction loss; This represents the approximate posterior distribution of the encoder output; represents the prior distribution of the latent variable; KL[·] represents the KL divergence; This represents the balance coefficient (usually 1).

[0228] Calculate the reconstruction error REW:

[0229] ;

[0230] In the formula: Representing the eigenvector Z x The dimension;

[0231] Preset error threshold (Generally 0.8), if the reconstruction error REW is greater than the error threshold, then the feature is determined. Z x This is early-stage, latent stress;

[0232] For characteristics identified as early latent stress Z x Calculate its relationship with the historical coercion database { Z 1, Z 2,…, Z n The Euclidean distance of}

[0233] ;

[0234] Choose the one with the smallest distance. h 1 sample, to obtain the initial adjustment amount:

[0235] ;

[0236] In the formula: C h,i Representation and Features Z x The smallest distance i The historical adjustment amount corresponding to each historical sample;

[0237] Combining reinforcement learning dynamic optimization, that is, the action space A in reinforcement learning dynamic optimization is... If the fine-tuning interval is centered, the reward function is updated as follows:

[0238] ;

[0239] In the formula: This indicates the change in photosynthetic efficiency; This indicates the change in root growth. This indicates the adjustment range of the hydrogen-rich water concentration; These represent the corresponding weight coefficients;

[0240] Latent stress adjustment for:

[0241] ;

[0242] In the formula: It represents the gradient of the Q-value with respect to the reward function, reflecting the impact of changes in the adjustment amount on the reward; This represents the fine-tuning coefficient (typically 0.01).

[0243] The specific method for obtaining it is as follows:

[0244] First, the Q-value is defined as "the long-term cumulative reward for performing adjustment action A in state Y," and its expression is the classic form of reinforcement learning, namely...

[0245] ;

[0246] but:

[0247] ;

[0248] In the formula: This represents the discount factor (typically 0.7 to 0.9).

[0249] The final concentration of hydrogen-rich water obtained is:

[0250] ;

[0251] In the formula: This represents the correction factor (i.e., the correction factor based on the physiological characteristics of the crop at its current growth stage); when the plant is in the flowering stage, When the plant is in the seedling stage, .

[0252] Example 2:

[0253] As a further optimization of the scheme in this application, based on the scheme in Example 1, the chlorophyll fluorescence parameters include the maximum optical efficiency ( F v / F m The potential activity of the photosystem II reaction center is reflected in the photosystem II protein complex located on the thylakoid membrane during photosynthesis. Its main function is to absorb light energy and catalyze the photolysis of water to produce oxygen, protons, and electrons. The photochemical quenching coefficient (qP, reflecting the open proportion of the photosystem II reaction center) is also considered.

[0254]

[0255] In the formula: F 0 indicates the initial fluorescence intensity; F m Indicates the maximum fluorescence intensity; F v Indicates variable fluorescence intensity; F t Indicates the steady-state fluorescence intensity under light adaptation; Indicates the maximum fluorescence intensity under light adaptation; Indicates the initial fluorescence intensity under light adaptation;

[0256] Canopy multispectral parameters include the Normalized Difference Vegetation Index (NDVI) and Red Edge Position (REP, a sensitive indicator reflecting changes in chlorophyll content):

[0257]

[0258] In the formula: Indicates near-infrared reflectance; Indicates the reflectivity in the red light band; , , These represent the reflectivity at wavelengths of 680nm, 705nm, and 750nm, respectively.

[0259] Degree of micro-deformation of the stem (measured by strain sensor):

[0260]

[0261] In the formula: This indicates the change in stem length; Indicates the initial length of the stem;

[0262] Leaf infrared temperature (canopy temperature, measured by an infrared sensor):

[0263]

[0264] In the formula: C1 and C2 represent the parameters related to Planck's constant (C1 is...). C2 is ); This represents the spectral radiance received by the sensor; Indicates the wavelength of the infrared band;

[0265] Root growth parameters (obtained through 3D image scanning) include root length density (RLD, reflecting the root absorption range) and root surface area index (RSI, reflecting absorption capacity):

[0266]

[0267] In the formula: L 根 Indicates the root system length; V 土壤 This indicates the corresponding soil volume; r 根 Indicates the root radius.

[0268] Field environmental data, such as light intensity, air temperature and humidity, soil temperature and humidity, CO2 concentration, and soil pH, can be obtained using corresponding light intensity sensors, air temperature and humidity sensors, soil temperature and humidity sensors, CO2 concentration sensors, and soil pH sensors.

[0269] Example 3:

[0270] As a further optimization of the present application, a dynamic adjustment system for hydrogen-rich water concentration is also included, which uses any one of the methods in Example 1 or Example 2 to adjust the concentration of hydrogen-rich water during the crop planting process.

Claims

1. A method for dynamically adjusting the concentration of hydrogen-rich water, characterized in that: include: Step S1: Data Acquisition: Based on the crop type and field conditions, deploy data acquisition sensors to collect plant growth trend data and field environmental factor data respectively. Step S2, Data Preprocessing: The collected data undergoes noise removal, abnormal data detection and correction, and data normalization preprocessing. Step S3: Extraction of key influencing factors: Explore plant growth trend data and field environmental factor data to identify key factors affecting crop growth; Key factors were extracted using a stress-metabolism coupling network, a dynamic time-delay network, a multi-scale feature fusion network, and a stress state adaptive network, respectively. Specifically, the preprocessed data variables from step S2 were input into the four networks in sequence. If a variable was identified as a key element by one of the four networks, the corresponding variable was incremented by 1 in the key factor determination. If the total value of a variable is not less than 3, then the variable is considered a key influencing factor; if the total value of a variable is 2, then it is considered a minor influencing factor. Step S4: Obtaining the initial concentration of hydrogen-rich water: Input the data of key influencing factors into the fuzzy logic system, define the fuzzy subsets of each key influencing factor, and obtain the initial concentration of hydrogen-rich water. Step S5, Obtaining the optimized concentration of hydrogen-rich water: The dynamic concentration adjustment amount, extreme environment adjustment amount, and latent stress adjustment amount are obtained through reinforcement learning dynamic optimization, extreme environment adaptation, and latent stress early warning, respectively, to obtain the final concentration of hydrogen-rich water.

2. The method for dynamically adjusting the concentration of hydrogen-rich water according to claim 1, characterized in that: In step S2, noise removal is performed by using the Kalman filter algorithm to remove random noise.

3. A method for dynamically adjusting the concentration of hydrogen-rich water according to claim 1 or 2, characterized in that: The abnormal data detection and correction in step S2 specifically involves: For input feature samples X =[ x 1, x 2,…, x n The Isolation Forest algorithm was used to detect anomalies in the data after removing random noise, and anomaly scores were obtained. ； In the formula: h(X) Representing feature samples X Path length in an isolated tree; E[h(X)] Indicates the expected path length; Represents Euler's constant; Preset abnormal score threshold SX d ,like If the corresponding feature sample is an outlier, then it is determined that the sample is an outlier. For outliers, Huber robust regression is used for correction. The Huber loss function is: ； ； In the formula: The parameters are constants; by optimizing the model parameters Improve tolerance to extreme environments.

4. The method for dynamically adjusting the concentration of hydrogen-rich water according to claim 3, characterized in that: The specific method for identifying key elements in the stress-metabolism coupling network is as follows: First, the preprocessed data variables from step S2 X’ Based on the effects of ROS on plants, the plants are divided into an environmental layer (CeX), a physiological layer (CeY), and an effect layer (CeZ). Among them, the environmental layer (CeX) corresponds to environmental indicators that directly induce ROS generation, the physiological layer (CeY) reflects indicators of ROS damage to plant function, and the effect layer (CeZ) reflects indicators of the final traits of plants affected by ROS. Then, two inter-layer unidirectional influence paths are established, including the first path CeX→CeY and the second path CeY→CeZ; Then, the path coefficients of the two paths are obtained using a partial least squares path model: , , Finally, preset path coefficient threshold ,like Then the stress-metabolism coupling network will identify the corresponding variables as key elements.

5. The method for dynamically adjusting the concentration of hydrogen-rich water according to claim 3, characterized in that: The specific method for identifying key elements using the dynamic time-delay network is as follows: First, obtain the preprocessed data variables from step S2. X'(t) Concentration requirements In time lag Mutual information under TDMI for:, ； In the formula: P represents the probability distribution; Representing variables X’ In time t Pick a Value and concentration requirement B in time Pick b The probability of the value; Representing variables X’ In time t Pick a The probability of the value This indicates that the concentration requirement B is in time. Pick b The probability of the value; Plant growth trend data and field environmental factor data were obtained at time lags. The TDMI is set below, with a preset TDMI peak threshold. max If the mutual information TDMI of a certain variable is greater than TDMI max Then the corresponding variable is determined to be a key element by the dynamic time-delay network.

6. The method for dynamically adjusting the concentration of hydrogen-rich water according to claim 3, characterized in that: The specific method for identifying key elements using the multi-scale feature fusion network is as follows: For time series variables X'(t) Perform three-level wavelet packet decomposition to obtain energy values ​​for eight frequency bands. E i,j : ； In the formula: x i,j (t) Indicates the first i The first variable j Sub-signals in each frequency band; t 0、 t n These represent the start and end times of the acquisition cycle corresponding to the sub-signal, respectively. The sub-signals of each frequency band are input as independent features into a random forest to obtain the Gini importance of each sub-signal; weights are assigned to different frequency bands based on the metabolic cycle caused by hydrogen-rich water. To obtain the overall importance of the corresponding variables: ； In the formula: Gini i,j Indicates the first i The first variable j The importance of the Gini coefficient in frequency band characteristics; Preset comprehensive importance threshold I d If the overall importance of a variable is greater than the overall importance threshold, then the corresponding variable is determined to be a key element by the multi-scale feature fusion network.

7. The method for dynamically adjusting the concentration of hydrogen-rich water according to claim 3, characterized in that: The specific method for identifying key elements by the stress state adaptive network is as follows: Introducing the Stress Index SI: ； In the formula: T leaf Indicates the infrared temperature of the blade; T opt This indicates the optimal temperature for the leaves; T max Indicates the critical high temperature; Indicates soil volumetric moisture content. Indicates field holding capacity; K 4. K 5 represents the corresponding weight coefficients; Preset stress index threshold SI d : when When the time is right, it indicates that the plant's metabolism is basically normal. Variables related to photosynthesis are used as input variables of the Bayesian network. The posterior probability of the variables and the concentration of hydrogen-rich water is calculated through the Bayesian network. The variables are arranged from largest to smallest and the top three variables are taken as key elements. when When the plant's photosynthesis / respiration is inhibited, stress-related variables are used as input variables of a Bayesian network. The posterior probability of the variables and the concentration of hydrogen-rich water is calculated through the Bayesian network. The variables are arranged from largest to smallest and the top three variables are selected as key elements. A Bayesian network is a probabilistic graphical model that outputs posterior probabilities based on the causal relationship between "field environmental factors / plant growth trends" and "hydrogen-rich water concentration requirements".

8. The method for dynamically adjusting the concentration of hydrogen-rich water according to claim 1, characterized in that: Step S4 specifically involves: Define the membership function of the fuzzy subset of key influencing factors, including a small number of memberships. U dr (S) Suitable membership degree U su (S) Excessive membership U we (S) : ； In the formula: S Indicates the current measured values ​​of key influencing factors; S min This indicates the critical value when the key influencing factor is too low. S a Indicates the appropriate central value of key influencing factors; S max This indicates the critical value when key influencing factors are too high. A fuzzy rule base is constructed by integrating multi-factor correlation and expert practical experience, and the centroid method is used to defuzzify and obtain the initial concentration. ； In the formula: U i Indicates the first i The membership degree of a rule; C i Indicates the first i The reference concentration for this rule; m This indicates the number of rules that are active.

Images