A method for evaluating physical and biological synergistic prevention and control effects based on apple snails

Abstract

The application discloses a kind of based on apple snail physical and biological synergistic prevention and treatment effect evaluation method, including collecting the experimental monitoring data and associated data of preset apple snail area, the experimental monitoring data and the associated data are pretreated;According to the biological control data, the resistance gene frequency and the natural enemy satiation threshold are obtained by the analysis of the measures sustainability of controlling snail, the natural enemy satiation threshold is revised using predator and prey dynamic model, the resistance gene frequency is revised using selective pressure;According to the physical control data, the interception efficiency is obtained by the analysis of intercepting trapping, the interception efficiency is adjusted using snail behavior response, temperature and humidity coupling, according to the synergistic prevention and treatment data, obtain governance efficiency data;According to the natural enemy satiation threshold, the resistance gene frequency, the interception efficiency and the governance efficiency data, construct three-dimensional apple snail physical and biological synergistic prevention and treatment effect evaluation model, and output evaluation result.

Description

Technical Field

[0001] This invention relates to the field of evaluating the synergistic effect of physical and biological control of golden apple snails, and particularly to a method for evaluating the synergistic effect of physical and biological control of golden apple snails. Background Technology

[0002] Golden apple snails, as a typical invasive alien species, have a strong reproductive capacity and wide adaptability, posing a serious threat to my country's farmland ecosystems, water conservancy facilities, and biodiversity. They feed on rice and other crop seedlings, leading to reduced crop yields, while also disrupting water balance, causing ecological chain disorder, and resulting in enormous losses to agricultural production and the ecological environment.

[0003] Current methods for controlling golden apple snails often employ a synergistic approach that combines physical interception and trapping with biological control by natural enemies. However, existing methods for evaluating control effectiveness have significant limitations: they only consider the physical interception rate or the amount of snails controlled by biological control, lacking a systematic analysis of the synergistic effect of the two; and they do not fully incorporate key variables such as environmental factors, snail resistance evolution, and socio-economic costs, leading to a large discrepancy between the evaluation results and the actual control effectiveness.

[0004] Furthermore, golden apple snails are prone to developing resistance gene mutations, and the effectiveness of natural enemies in controlling snails is significantly affected by environmental conditions such as temperature and prey substitution. Physical interception efficiency also fluctuates due to snail behavior responses and changes in temperature and humidity. Traditional assessment methods struggle to quantify the impact of these dynamic factors, failing to provide precise evidence for optimizing control programs and hindering the promotion and application of collaborative control technologies. Therefore, there is an urgent need to construct a multi-dimensional and refined method for evaluating the effectiveness of collaborative control, comprehensively considering core dimensions such as ecology, technology, and economics, scientifically quantifying control effectiveness, and providing technical support for precise control of golden apple snails. Summary of the Invention

[0005] The purpose of this invention is to provide a method for evaluating the synergistic effect of physical and biological control of golden apple snails.

[0006] To achieve the above objectives, the present invention is implemented according to the following technical solution:

[0007] This invention includes the following steps:

[0008] Experimental monitoring data and related data were collected from a pre-defined area of ​​*Pomacea canaliculata*, and the experimental monitoring data and related data were preprocessed. The experimental monitoring data included physical control data, biological control data, and synergistic control data. The related data included environmental data, ecological data, and economic data. The synergistic control data included synergistic gain coefficient, combination of control measures, and implementation cycle data.

[0009] Based on the biological control data, a sustainability analysis of snail control measures was conducted to obtain the predator saturation threshold and resistance gene frequency. The predator saturation threshold was corrected using a predator-prey dynamic model, and the resistance gene frequency was corrected using selective pressure.

[0010] Interception efficiency is obtained by interception and trapping analysis based on the physical control data. The interception efficiency is adjusted by snail behavior response and temperature and humidity coupling. Governance effectiveness data is obtained based on the collaborative control data.

[0011] A three-dimensional physical and biological synergistic control effect evaluation model for golden apple snails is constructed based on the predator saturation threshold, the resistance gene frequency, the interception efficiency, and the control effectiveness data. The data to be evaluated is input into the golden apple snail physical and biological synergistic control effect evaluation model, and the evaluation results are output. The three dimensions include ecological dimension, technological synergy dimension, and socio-economic adaptability dimension.

[0012] Furthermore, the method for obtaining the predator satiety threshold and the resistance gene frequency includes:

[0013] Data on daily snail predation counts, predator numbers, snail population density, average daily temperature, and prey substitute biomass were obtained. The baseline satiety threshold was calculated using a Michaelis-Menten dynamics model.

[0014] ;

[0015] in The basic predator satiety threshold, For the attack rate of natural enemies, Time for prey processing The effective predation time per day, The number of snails preyed upon per day. The number of natural enemies;

[0016] Snail samples are collected quarterly, and resistance-related genes are amplified using polymerase chain reaction. The number of resistance alleles and susceptible alleles is counted, and the initial resistance gene frequency is calculated.

[0017]

[0018] in The initial resistance gene frequency, It is a resistance allele. It is a sensitive allele. The number of homozygous resistant individuals. The number of heterozygous individuals. This represents the number of samples.

[0019] Furthermore, the method for correcting the predator satiety threshold using a predator-prey dynamic model includes:

[0020] Introducing a two-factor correction for the baseline satiety threshold, the expression is:

[0021] ;

[0022] in This is a temperature correction factor. For prey modification factor, To replace food biomass with prey This is the corrected predator satiety threshold. The basic predator satiety threshold;

[0023] When the measured snail density is less than the corrected saturation threshold for three consecutive days, it is determined that the saturation threshold has been triggered.

[0024] Furthermore, the method of selectively adjusting the frequency of the resistance gene using selective pressure includes:

[0025] Determine the initial resistance gene frequency, record the population size, and calculate the selection pressure index:

[0026] ;

[0027] in To select the stress index, Let be the snail population size in generation t. The population size of snails in generation t+1. The frequency of the resistance gene in generation t. This represents the contribution coefficient of resistance genes to population growth.

[0028] The selection coefficients for different genotypes are adjusted based on the selection pressure index, expressed as follows:

[0029] ；

[0030] ;

[0031] ;

[0032] in Sensitive homozygote The selection coefficient, Heterozygote The selection coefficient, resistant homozygote Selection coefficient;

[0033] Calculate the average fitness and the corrected gene frequency:

[0034]

[0035]

[0036] in The frequency of the resistance gene in generation t. To correct the resistance gene frequency in generation t+1, Let be the average fitness of generation t.

[0037] Furthermore, the method for obtaining interception efficiency by performing interception and trapping analysis based on the physical prevention and control data includes:

[0038] Obtain information on the mesh size of the interception net, the deployment density, the number of trapping devices, the net material, the amount of snails coming from upstream, the amount of snails intercepted, the snail movement rate, the flow velocity, the water temperature, and the pH value.

[0039] Calculate the mechanical interception rate based on the physical barrier capability of the mesh:

[0040] ;

[0041] Among them The maximum body width of a mature golden apple snail, For standard deployment density, For mechanical interception rate, To intercept the mesh aperture, For deployment density;

[0042] Based on the active attraction effect of the trap on snails, the contribution of the trapping device is calculated:

[0043] ;

[0044] in The daily snail capture rate of the trapping device. Standard trap density, The trapping density gain coefficient, For the density of the traps, This refers to the amount of snails coming from upstream. Contribute to the trapping device;

[0045] The base interception rate is calculated based on the mechanical interception rate and the contribution of the trapping device:

[0046] ;

[0047] in This represents the maximum mechanical interception rate. This is the weighting coefficient for the mechanical interception rate. The weighting coefficients contributed to the trapping device. Based on the interception rate.

[0048] Furthermore, the method of adjusting the interception efficiency by employing snail behavior response and temperature and humidity coupling includes:

[0049] By monitoring the behavioral characteristics of golden apple snails after encountering obstacles using high-speed video, a three-dimensional correction term was introduced to calculate the snail behavior response correction coefficient:

[0050] ;

[0051] in This is a correction factor for the behavior response of snails. For turning frequency, Based on the interception rate, To explore the proportion of time, Due to the cluster effect, The weighting coefficient for turning frequency. To explore the weighting coefficients for time percentages, This represents the weighting coefficient for the clustering effect;

[0052] Introducing a temperature and humidity coupling compensation coefficient, the expression is:

[0053] ;

[0054] in This is the temperature and humidity coupling compensation coefficient. For temperature sensitivity coefficient, Humidity sensitivity coefficient This is the actual measured water temperature;

[0055] The base interception rate is corrected based on the snail behavior response correction coefficient and the temperature and humidity coupling compensation coefficient, expressed as follows:

[0056] ;

[0057] in This is the corrected interception rate.

[0058] Furthermore, the method for constructing the evaluation model for the synergistic physical and biological control effect of the golden apple snail includes:

[0059] Ecological dimension assessment indicators are calculated based on water health, crop damage rate, and non-target organism impact:

[0060]

[0061] in The safe threshold for water health The appropriate threshold for water health For water body health and safety values, For the optimal value for water health, For crop damage rate, Corrected resistance gene frequencies As an ecological dimension assessment indicator, The impact rate on non-target organisms;

[0062] The evaluation index for technical synergy is calculated based on interception efficiency, resistance gene frequency, predator satiety threshold, and synergy gain coefficient.

[0063]

[0064] in As an evaluation indicator for the dimension of technological collaboration, For the cooperative gain, This is the corrected interception rate. This is the corrected predator satiety threshold. This is the synergistic effect coefficient;

[0065] The socioeconomic suitability assessment indicators are calculated based on prevention and control costs, reduction in rice loss, and ecosystem service value.

[0066]

[0067] in For the loss of prevention, To prevent future losses, For the cost of physical prevention and control measures, For the cost of biological control measures, For ecological service value, For total investment, As an indicator for assessing socioeconomic adaptability, Weight for ecological services For prevention and control weight;

[0068] The ecological dimension assessment indicators, the technological synergy dimension assessment indicators, and the socio-economic adaptability dimension assessment indicators are weighted and combined to obtain a comprehensive assessment index.

[0069] The beneficial effects of this invention are:

[0070] This invention is a method for evaluating the synergistic effect of physical and biological control of golden apple snails. Compared with existing technologies, this invention has the following technical advantages:

[0071] This invention constructs a three-dimensional assessment model encompassing ecology, technological synergy, and socioeconomic adaptability, overcoming the limitations of single-dimensional assessments and achieving comprehensive quantification of control effects. By modifying key parameters through predator-prey dynamic models and selective pressure, and by adjusting interception efficiency in conjunction with snail behavior responses and temperature / humidity coupling, the accuracy of the assessment is enhanced. It integrates experimental monitoring data and correlated data, taking into account the sustainability of biological snail control, the effectiveness of physical interception, and synergistic gains, adapting to complex environmental changes. The assessment results can directly guide the optimization of control programs, reducing costs and ecological losses, combining scientific rigor with practicality, and providing reliable technical support for the precise control of golden apple snails. Attached Figure Description

[0072] Figure 1 This is a flowchart illustrating the steps of an evaluation method for the synergistic physical and biological control of golden apple snails according to the present invention. Detailed Implementation

[0073] The present invention will be further described below through specific embodiments. The illustrative embodiments and descriptions herein are used to explain the present invention, but are not intended to limit the present invention.

[0074] The present invention provides a method for evaluating the synergistic physical and biological control effects of golden apple snails, comprising the following steps:

[0075] like Figure 1 As shown, this embodiment includes the following steps:

[0076] Experimental monitoring data and related data were collected from a pre-defined area of ​​*Pomacea canaliculata*, and the experimental monitoring data and related data were preprocessed. The experimental monitoring data included physical control data, biological control data, and synergistic control data. The related data included environmental data, ecological data, and economic data. The synergistic control data included synergistic gain coefficient, combination of control measures, and implementation cycle data.

[0077] In the actual assessment, a rice-growing area in southern China (latitude 28°N, longitude 114°E) was selected as the pre-set monitoring area for golden apple snails. This area covered 50 mu (approximately 3.3 hectares) and was divided into five 10-mu (approximately 0.67 hectares) experimental plots. The rice variety was conventional indica rice, with a planting density of 300,000 plants per mu (approximately 20,000 plants per hectare). The average initial density of golden apple snails in the area was 12 snails per m². 2 The main control method adopted is a synergistic control model of physical interception nets and the regulation of natural enemies of grass carp, with a control period from April to October 2024 (the entire growth period of rice).

[0078] Physical control data: The mesh size of the interception net is 1.5cm, the maximum body width of a mature golden apple snail is 2.5cm, the deployment density is 8 nets / mu, and the standard deployment density is 10 nets / mu; the number of trapping devices is 5 / mu, the standard trapping device density is 4 / mu, and the trapping density gain coefficient is 0.3; the average daily snail arrival from upstream is 35, and the daily snail capture by the trapping devices is 12; the net material is nylon net, the snail movement speed is 0.8m / h, the water flow velocity is 0.3m / s, the measured water temperature is 28℃, and the pH value is 7.2;

[0079] Biological control data: The number of natural predators, grass carp, was 20 per mu (approximately 0.067 hectares), and the daily predation rate of snails was 15 per fish. The predator attack rate was 0.6. Through controlled indoor experiments, five grass carp density gradients (5-25 fish / m²) were set up, and the predation rate of grass carp on golden apple snails at different densities was recorded. The predator attack rate was found to be 0.6 using the Michaelis-Menten equation. The prey treatment time was 0.8 h per snail, and the effective predation time per day was 8 h. The prey substitute food biomass was 1.2 kg / mu (approximately 0.067 hectares). 300 snail samples were collected each quarter, including 18 homozygous resistant individuals and 45 heterozygous individuals.

[0080] Synergistic control data: Synergistic gain was 15%; loss without control was 80 kg / mu, loss after control was 12 kg / mu; cost of physical control measures was 120 yuan / mu, cost of biological control measures was 80 yuan / mu;

[0081] Environmental data: The average daily temperature is 26℃, the relative humidity is 82%, the temperature sensitivity coefficient is 0.2, and the humidity sensitivity coefficient is 0.15; the total nitrogen in the water is 1.2 mg / L, and the safe threshold for water health is 1.5 mg / L; the dissolved oxygen is 6.8 mg / L, and the suitable threshold for water health is 7.0 mg / L.

[0082] Ecological data: Crop damage rate was 3.5%; the corrected frequency of resistance genes was 0.12; the ecological service value was 200 yuan / mu, and the total input was 450 yuan / mu;

[0083] Economic data: Prevention and control weight is 0.6, ecosystem service weight is 0.4; the contribution coefficient of resistance genes to population growth is 0.15; the population size of snails in generation t is 12 individuals / m². 2 The snail population size in generation t+1 was 8 individuals / m². 2 ;

[0084] Outlier removal and standardization were performed on the collected data: the ratio of upstream incoming screw volume to intercepted screw volume was 35 / (12×5) = 0.58 < 1.2, indicating no outlier data; physical parameters such as the aperture and deployment density of the interception net have been standardized according to standard values.

[0085] Based on the biological control data, a sustainability analysis of snail control measures was conducted to obtain the predator saturation threshold and resistance gene frequency. The predator saturation threshold was corrected using a predator-prey dynamic model, and the resistance gene frequency was corrected using selective pressure.

[0086] Interception efficiency is obtained by interception and trapping analysis based on the physical control data. The interception efficiency is adjusted by snail behavior response and temperature and humidity coupling. Governance effectiveness data is obtained based on the collaborative control data.

[0087] In actual assessments, the data on the effectiveness of the treatment include the rate of change in snail density before and after treatment, the reduction in crop damage rate, and the degree of improvement in the aquatic environment.

[0088] A three-dimensional physical and biological synergistic control effect evaluation model for golden apple snails is constructed based on the predator saturation threshold, the resistance gene frequency, the interception efficiency, and the control effectiveness data. The data to be evaluated is input into the golden apple snail physical and biological synergistic control effect evaluation model, and the evaluation results are output. The three dimensions include ecological dimension, technological synergy dimension, and socio-economic adaptability dimension.

[0089] In this embodiment, the method for obtaining the predator satiety threshold and the resistance gene frequency includes:

[0090] Data on daily snail predation counts, predator numbers, snail population density, average daily temperature, and prey substitute biomass were obtained. The baseline satiety threshold was calculated using a Michaelis-Menten dynamics model.

[0091] ;

[0092] in The basic predator satiety threshold, For the attack rate of natural enemies, Time for prey processing The effective predation time per day, The number of snails preyed upon per day. The number of natural enemies;

[0093] Snail samples are collected quarterly, and resistance-related genes are amplified using polymerase chain reaction. The number of resistance alleles and susceptible alleles is counted, and the initial resistance gene frequency is calculated.

[0094]

[0095] in The initial resistance gene frequency, It is a resistance allele. It is a sensitive allele. The number of homozygous resistant individuals. The number of heterozygous individuals. The number of samples;

[0096] In the actual assessment, 300 snail samples were collected each quarter, with a baseline satiety threshold of 3.15 snails / m². 2 The initial resistance gene frequency was 0.135.

[0097] In this embodiment, the method for correcting the predator satiety threshold using a predator-prey dynamic model includes:

[0098] Introducing a two-factor correction for the baseline satiety threshold, the expression is:

[0099] ;

[0100] in This is a temperature correction factor. For prey modification factor, To replace food biomass with prey This is the corrected predator satiety threshold. The basic predator satiety threshold;

[0101] When the measured snail density is less than the corrected saturation threshold for 3 consecutive days, it is determined that the saturation threshold is triggered.

[0102] In the actual assessment, the temperature correction factor was 0.861, the prey correction factor was 1.024, and the corrected predator satiety threshold was 3.09 predators / m². 2 The actual snail density monitored for three consecutive days was 2.8 snails / m³. 2 2.6 birds / m 2 2.9 birds / m 2 All values ​​are less than the modified satiety threshold, indicating that the satiety threshold has been triggered.

[0103] In this embodiment, the method of selectively adjusting the frequency of the resistance gene using selective pressure includes:

[0104] Determine the initial resistance gene frequency, record the population size, and calculate the selection pressure index:

[0105] ;

[0106] in To select the stress index, Let be the snail population size in generation t. The population size of snails in generation t+1. The frequency of the resistance gene in generation t. This represents the contribution coefficient of resistance genes to population growth.

[0107] The selection coefficients for different genotypes are adjusted based on the selection pressure index, expressed as follows:

[0108] ；

[0109] ;

[0110] ;

[0111] in Sensitive homozygote The selection coefficient, Heterozygote The selection coefficient, resistant homozygote Selection coefficient;

[0112] Calculate the average fitness and the corrected gene frequency:

[0113]

[0114]

[0115] in The frequency of the resistance gene in generation t. To correct the resistance gene frequency in generation t+1, The average fitness of generation t; the initial value of the resistance gene frequency in generation t is equal to the initial resistance gene frequency. ;

[0116] In the actual assessment, a stress index of 0.68 was selected. It is 0.136. It is 0.068. The value was 0.932, the average fitness was 0.3372, and the corrected gene frequency was 0.43.

[0117] In this embodiment, the method for obtaining interception efficiency by interception and trapping analysis based on the physical prevention and control data includes:

[0118] Obtain information on the mesh size of the interception net, the deployment density, the number of trapping devices, the net material, the amount of snails coming from upstream, the amount of snails intercepted, the snail movement rate, the flow velocity, the water temperature, and the pH value.

[0119] Calculate the mechanical interception rate based on the physical barrier capability of the mesh:

[0120] ;

[0121] Among them The maximum body width of a mature golden apple snail, For standard deployment density, For mechanical interception rate, To intercept the mesh aperture, For deployment density;

[0122] Based on the active attraction effect of the trap on snails, the contribution of the trapping device is calculated:

[0123] ;

[0124] in The daily snail capture rate of the trapping device. Standard trap density, The trapping density gain coefficient, For the density of the traps, This refers to the amount of snails coming from upstream. Contribute to the trapping device;

[0125] The base interception rate is calculated based on the mechanical interception rate and the contribution of the trapping device:

[0126] ;

[0127] in This represents the maximum mechanical interception rate. This is the weighting coefficient for the mechanical interception rate. The weighting coefficients contributed to the trapping device. Base interception rate;

[0128] In actual assessment, when the ratio of the number of incoming snails to the number of intercepted snails in a single group is greater than 1.2, it is judged as data anomaly and corrected by the 3-day moving average method; standardization processing: the aperture of different specifications of nets is uniformly converted into relative aperture ratio; the mechanical interception rate is 0.32, the contribution of the trapping device is 0.471, and the basic interception rate is 0.4284.

[0129] In this embodiment, the method of adjusting the interception efficiency using snail behavior response and temperature and humidity coupling includes:

[0130] By monitoring the behavioral characteristics of golden apple snails after encountering obstacles using high-speed video, a three-dimensional correction term was introduced to calculate the snail behavior response correction coefficient:

[0131] ;

[0132] in This is a correction factor for the behavior response of snails. For turning frequency, Based on the interception rate, To explore the proportion of time, Due to the cluster effect, The weighting coefficient for turning frequency. To explore the weighting coefficients for time percentages, This represents the weighting coefficient for the clustering effect;

[0133] Introducing a temperature and humidity coupling compensation coefficient, the expression is:

[0134] ;

[0135] in This is the temperature and humidity coupling compensation coefficient. For temperature sensitivity coefficient, Humidity sensitivity coefficient This is the actual measured water temperature;

[0136] The base interception rate is corrected based on the snail behavior response correction coefficient and the temperature and humidity coupling compensation coefficient, expressed as follows:

[0137] ;

[0138] in This is the corrected interception rate;

[0139] In the actual assessment, the turning frequency is the proportion of time that snails spend exploring meshes, derived from infrared counters and high-speed cameras; the clustering effect is the snail aggregation density in the same area, derived from quadrat method surveys; and the exploration time proportion is the proportion of time that golden apple snails spend exploring meshes, derived from behavioral observation experiments.

[0140] , , The value ranges from 0.01 to 0.05 and can be determined through orthogonal experiments based on the species of snails.

[0141] When the water temperature is below 15℃ or above 35℃, the activity intensity of the golden apple snail decreases significantly, and the influence of humidity decreases. Therefore, the humidity coupling compensation coefficient is automatically corrected to 0.7. When the cluster density is greater than 2 birds / m 2 hour, Increased to 0.05; turning frequency was 0.3 times / min, exploration time percentage was 0.25, and swarm effect was 1.8 individuals / m². 2 ; It is 0.03. It is 0.04. The coefficient for snail behavior response is 0.02, the coefficient for temperature and humidity coupling compensation is 1.082, and the corrected interception rate is 0.489.

[0142] In this embodiment, the method for constructing the evaluation model for the synergistic physical and biological control effect of the golden apple snail includes:

[0143] Ecological dimension assessment indicators are calculated based on water health, crop damage rate, and non-target organism impact:

[0144]

[0145] in The safe threshold for water health The appropriate threshold for water health For water body health and safety values, For the optimal value for water health, For crop damage rate, Corrected resistance gene frequencies As an ecological dimension assessment indicator, The impact rate on non-target organisms;

[0146] The evaluation index for technical synergy is calculated based on interception efficiency, resistance gene frequency, predator satiety threshold, and synergy gain coefficient.

[0147]

[0148] in As an evaluation indicator for the dimension of technological collaboration, For the cooperative gain, This is the corrected interception rate. This is the corrected predator satiety threshold. This is the synergistic effect coefficient;

[0149] The socioeconomic suitability assessment indicators are calculated based on prevention and control costs, reduction in rice loss, and ecosystem service value.

[0150]

[0151] in For the loss of prevention, To prevent future losses, For the cost of physical prevention and control measures, For the cost of biological control measures, For ecological service value, For total investment, As an indicator for assessing socioeconomic adaptability, Weight for ecological services For prevention and control weight;

[0152] The ecological dimension assessment indicators, the technological synergy dimension assessment indicators, and the socio-economic adaptability dimension assessment indicators are weighted and combined to obtain a comprehensive assessment index;

[0153] In actual assessments, a comprehensive assessment index greater than or equal to 0.8 indicates highly efficient and synergistic control effects; a comprehensive assessment index greater than or equal to 0.6 and less than 0.8 indicates basically synergistic control effects; a comprehensive assessment index greater than or equal to 0.4 and less than 0.6 indicates control effects that need optimization; and a comprehensive assessment index less than 0.4 indicates synergistic failure in control effects.

[0154] The ecological dimension assessment index is 0.5748, the technological synergy dimension assessment index is 1.738, the socio-economic adaptability dimension assessment index is 0.4, and the comprehensive assessment index is 0.6858. The current regional physical and biological synergistic control effect of golden apple snail is basically synergistic.

[0155] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A method for evaluating the physical and biological synergistic control effect of apple snails, characterized in that, Includes the following steps: Experimental monitoring data and related data were collected from a pre-defined area of ​​*Pomacea canaliculata*, and the experimental monitoring data and related data were preprocessed. The experimental monitoring data included physical control data, biological control data, and synergistic control data. The related data included environmental data, ecological data, and economic data. The synergistic control data included synergistic gain coefficient, combination of control measures, and implementation cycle data. Based on the biological control data, a sustainability analysis of snail control measures was conducted to obtain the predator saturation threshold and resistance gene frequency. The predator saturation threshold was corrected using a predator-prey dynamic model, and the resistance gene frequency was corrected using selective pressure. Interception efficiency is obtained by interception and trapping analysis based on the physical control data. The interception efficiency is adjusted by snail behavior response and temperature and humidity coupling. Governance effectiveness data is obtained based on the collaborative control data. A three-dimensional physical and biological synergistic control effect evaluation model for golden apple snails is constructed based on the predator saturation threshold, the resistance gene frequency, the interception efficiency, and the control effectiveness data. The data to be evaluated is input into the golden apple snail physical and biological synergistic control effect evaluation model, and the evaluation results are output. The three dimensions include ecological dimension, technological synergy dimension, and socio-economic adaptability dimension.

2. The method for evaluating the synergistic physical and biological control effect of golden apple snails according to claim 1, characterized in that, A method for obtaining the predator satiety threshold and the resistance gene frequency includes: Data on daily snail predation counts, predator numbers, snail population density, average daily temperature, and prey substitute biomass were obtained. The baseline satiety threshold was calculated using a Michaelis-Menten dynamics model. ； in The basic predator satiety threshold, For the attack rate of natural enemies, Time for prey processing The effective predation time per day, This represents the number of snails consumed daily. The number of natural enemies; Snail samples are collected quarterly, and resistance-related genes are amplified using polymerase chain reaction. The number of resistance alleles and susceptible alleles is counted, and the initial resistance gene frequency is calculated. ； in The initial resistance gene frequency, It is a resistance allele. It is a sensitive allele. The number of homozygous resistant individuals. The number of heterozygous individuals. This represents the number of samples.

3. The method according to claim 1, wherein the physical and biological synergistic control effect of the apple snail is evaluated. A method for correcting the predator satiety threshold using a predator-prey dynamic model includes: Introducing a two-factor correction for the baseline satiety threshold, the expression is: ; in This is a temperature correction factor. For prey modification factor, To replace food biomass with prey This is the corrected predator satiety threshold. The basic predator satiety threshold; When the measured snail density is less than the corrected saturation threshold for three consecutive days, it is determined that the saturation threshold has been triggered.

4. The method according to claim 1, wherein the physical and biological synergistic control effect of the apple snail is evaluated. A method for selectively adjusting the frequency of the resistance gene using selective pressure includes: Determine the initial resistance gene frequency, record the population size, and calculate the selection pressure index: ； in To select the stress index, Let be the snail population size in generation t. The population size of snails in generation t+1. The frequency of the resistance gene in generation t. This represents the contribution coefficient of resistance genes to population growth. The selection coefficients for different genotypes are adjusted based on the selection pressure index, expressed as follows: ； ； ； in Sensitive homozygote The selection coefficient, Heterozygote The selection coefficient, resistant homozygote Selection coefficient; Calculate the average fitness and the corrected gene frequency: ； ； in The frequency of the resistance gene in generation t. To correct the resistance gene frequency in generation t+1, Let be the average fitness of generation t.

5. The method according to claim 1, wherein the method is characterized by, A method for obtaining interception efficiency by interception and trapping analysis based on the physical prevention and control data includes: Obtain information on the mesh size of the interception net, the deployment density, the number of trapping devices, the net material, the amount of snails coming from upstream, the amount of snails intercepted, the snail movement rate, the flow velocity, the water temperature, and the pH value. Calculate the mechanical interception rate based on the physical barrier capability of the mesh: ； wherein the maximum body width of the adult golden apple snail, is the standard laying density, is the mechanical interception rate, is the interception mesh size, is the laying density; Based on the active attraction effect of the trap on snails, the contribution of the trapping device is calculated: ； wherein is the daily capture of moths by the trap, is the standard trap density, is the trap density gain factor, is the trap density, is the upstream moth flow, is the contribution of the trap. The base interception rate is calculated based on the mechanical interception rate and the contribution of the trapping device: ； wherein is the maximum value of the mechanical interception rate, is the weight coefficient of the mechanical interception rate, is the weight coefficient contributed by the trapping device, is the base interception rate.

6. The method for evaluating the synergistic physical and biological control effect of golden apple snails according to claim 1, characterized in that, The method for adjusting the interception efficiency by employing snail-like behavior response and temperature and humidity coupling includes: By monitoring the behavioral characteristics of golden apple snails after encountering obstacles using high-speed video, a three-dimensional correction term was introduced to calculate the snail behavior response correction coefficient: ； in This is a correction factor for the behavior response of snails. For turning frequency, Based on the interception rate, To explore the proportion of time, Due to the cluster effect, The weighting coefficient for turning frequency. To explore the weighting coefficients for time percentages, This represents the weighting coefficient for the clustering effect; Introducing a temperature and humidity coupling compensation coefficient, the expression is: ； wherein is a temperature and humidity coupling compensation coefficient, is a temperature sensitive coefficient, is a humidity sensitive coefficient, is a measured water temperature; The base interception rate is corrected based on the snail behavior response correction coefficient and the temperature and humidity coupling compensation coefficient, expressed as follows: ； wherein is the corrected interception rate.

7. The method for evaluating the synergistic physical and biological control effect of golden apple snails according to claim 1, characterized in that, The method for constructing the evaluation model for the synergistic physical and biological control effect of the golden apple snail includes: Ecological dimension assessment indicators are calculated based on water health, crop damage rate, and non-target organism impact: ； in The safe threshold for water health The appropriate threshold for water health For water body health and safety values, For the optimal level of water health, For crop damage rate, Corrected resistance gene frequencies As an ecological dimension assessment indicator, The impact rate on non-target organisms; The evaluation index for technical synergy is calculated based on interception efficiency, resistance gene frequency, predator satiety threshold, and synergy gain coefficient. ； in As an evaluation indicator for the dimension of technological collaboration, For the cooperative gain, This is the corrected interception rate. This is the corrected predator satiety threshold. This is the synergistic effect coefficient; The socioeconomic suitability assessment indicators are calculated based on prevention and control costs, reduction in rice loss, and ecosystem service value. ； in For the loss of prevention, To prevent future losses, For the cost of physical prevention and control measures, For the cost of biological control measures, For ecological service value, For total investment, As an indicator for assessing socioeconomic adaptability, Weight for ecological services For prevention and control weight; The ecological dimension assessment indicators, the technological synergy dimension assessment indicators, and the socio-economic adaptability dimension assessment indicators are weighted and combined to obtain a comprehensive assessment index.

Images