Method for purifying and cultivating SPF duck closed colony

Abstract

The present application belongs to the technical field of data processing, and particularly relates to a purification and cultivation method for SPF duck closed group, comprising: moving the to-be-purified basic duck group into an isolated environment with independent ventilation, positive pressure protection and full-in full-out management, sterilizing and purifying the incoming goods and personnel; collecting multi-dimensional physical and chemical signals and biochemical detection signals, generating physical and discrete sequences containing time stamps, time matching according to relative time difference and air pressure deviation, calculating comprehensive deviation degree in combination with environmental offset characteristics and duck group genetic quality parameters, and adjusting the ventilation equipment according to the comprehensive deviation degree and air flow rate; obtaining the microbiological detection results of the ex vivo samples of the to-be-purified duck group, excluding individuals not meeting the preset microbiological standards from the breeding group, making the qualified individuals reproduce internally under the condition of not introducing external ducks, and continuously rechecking the offspring for several generations; and updating the early warning threshold of the next cultivation cycle according to the hatching performance of the breeding batch. The present application improves the stability and controllability of purification and cultivation.

Description

Technical Field

[0001] This invention relates to the field of data processing technology. More specifically, this invention relates to a method for the purification and breeding of SPF duck closed flocks. Background Technology

[0002] my country is the world's largest duck producer. However, farming efficiency has long been hampered by bottlenecks in disease control. Farmed duck flocks are prone to carrying various infectious pathogens, which not only cause significant economic losses to the farming industry but also seriously threaten public health and safety.

[0003] The establishment of a disease prevention and control system requires standardized laboratory animals as raw materials for research. Conventionally farmed duck flocks, carrying pathogens, can severely interfere with scientific experimental results, leading to a series of problems such as delays in vaccine development, false positives / false negatives in diagnostic reagents, and contamination of biological products. Therefore, establishing SPF (Specific Pathogen Free) duck populations with clearly defined microbial control is a crucial foundation for overcoming technical challenges in disease prevention and control and enhancing industrial capabilities.

[0004] SPF ducks are artificially bred ducks that do not carry zoonotic pathogens or highly contagious animal diseases that pose serious threats to animal and / or human health, nor pathogens that significantly interfere with scientific research. They have a clear genetic background or origin and are used for scientific research, teaching, production, testing, and other scientific experiments. SPF ducks and SPF duck embryos are key raw materials for the research and production of veterinary biological products such as duck plague vaccines and duck hepatitis vaccines. Establishing SPF duck flocks effectively avoids exogenous viral contamination, ensuring the safety and effectiveness of biological products. SPF ducks are highly sensitive to major duck disease pathogens, making them ideal models for studying the pathogenesis, immune response patterns, and vaccine efficacy evaluation of diseases such as duck viral hepatitis, duck influenza, and duck reovirus infection. Compared to ordinary ducks, SPF ducks do not have maternal antibody interference or latent infection, accurately reflecting the true pathogenicity of pathogens and the true protective efficacy of vaccines. SPF ducks are essential experimental materials in key stages such as pathogen pathogenesis research, diagnostic reagent development and quality control, and vaccine development and efficacy evaluation.

[0005] As standardized laboratory animals, ensuring the reliability and reproducibility of experiments is crucial. The breeding and preservation methods for SPF ducks in closed flocks generally include several steps: artificial domestication, purification and selection, family selection, and closed flock selection. In the purification and selection process of SPF duck flocks, especially in the daily operation of closed breeding and breeding isolation, implementing isolation environment early warning and control based on multi-source data can effectively protect the health of the duck flock and improve the success rate of selection. Its main task is to continuously analyze multi-dimensional physicochemical and biochemical detection signals collected from the environment to promptly detect environmental deviations or pathogen invasions and adjust equipment, thereby avoiding duck flock damage and breeding failure caused by environmental deterioration or disease spread.

[0006] Existing isolated environment monitoring typically uses independent sensors to collect environmental data and relies on fixed upper and lower threshold methods to trigger the operation of variable frequency motors or ventilation regulating valves.

[0007] However, in actual continuous monitoring and control, early warning and control face challenges such as misalignment of data time points and the cross-cycle decline of the controlled object's physical condition. Firstly, the sampling frequencies of multidimensional physicochemical signals and discrete sequences differ, causing gaps and misalignments in the underlying data over time. Traditional monitoring methods rely on the extreme values ​​of data from a single sensor as independent judgment criteria, making it difficult to effectively correlate data from different sources. The lack of a time-axis matching mechanism between timestamped physical sequences and discrete sequences easily leads to missed anomalies or misjudging routine environmental fluctuations as risks.

[0008] Secondly, existing technologies do not take into account the decline in physical condition of SPF duck flocks during cross-cycle breeding. Directly applying fixed parameters will lead to a disconnect between safety standards and actual breeding conditions. At the same time, when faced with the need for regulation and intervention when environmental parameters exceed limits, a single step control command is difficult to ensure the stability of equipment operation. Directly outputting large adjustment commands can easily cause airflow impact, thereby causing secondary stress damage to the duck flock. Summary of the Invention

[0009] To address the technical problems of inaccurate early warnings due to time misalignment of multi-source data, fixed safety limits, and stress caused by excessive equipment movements in duck flocks, this invention provides solutions in the following aspects.

[0010] In a first aspect, the present invention provides a method for purifying and breeding an SPF duck closed flock, comprising: transferring the basic duck flock to be purified into an isolation environment for rearing, wherein the isolation environment adopts a high-efficiency filtration device to achieve independent ventilation, positive pressure protection, and all-in / all-out; sterilizing items entering the isolation environment, and managing the changing clothes, bathing, and disinfection of personnel entering the isolation environment; monitoring the isolation environment using multi-source data, and adjusting the ventilation equipment in the isolation environment according to the monitoring results; obtaining the microbiological test results of in vitro samples of the basic duck flock to be purified, isolating individuals whose microbiological test results do not meet the preset microbiological standards from the breeding flock, and not including such individuals in the next generation of the breeding flock; not introducing external ducks during the breeding period, allowing individuals whose microbiological test results meet the preset microbiological standards to reproduce within the flock, and retesting the offspring for microbiology, and identifying duck flocks that meet the preset microbiological standards for consecutive preset generations as an SPF duck closed flock.

[0011] Preferably, the step of monitoring the isolation environment using multi-source data and adjusting the ventilation equipment in the isolation environment based on the monitoring results includes: acquiring historical concentration sequences to determine concentration limits; collecting biochemical detection signals and multidimensional physicochemical signals including ambient temperature, relative humidity, air pressure, ammonia concentration, and air velocity; calculating the arithmetic mean of each signal in the multidimensional physicochemical signals within a preset sliding time window to obtain the corresponding smoothed mean, and generating a timestamped physical sequence based on the smoothed mean; generating a discrete sequence based on the biochemical detection signals and concentration limits; obtaining the time series deviation based on the relative time difference between the physical sequence and the discrete sequence, and the air pressure deviation between the current transient sample value of the air pressure signal in the physical sequence and its smoothed mean; and performing operations on the physical sequence and the discrete sequence based on the time series deviation. Time axis matching is used to construct a multi-source data array. If no step abrupt changes are detected in the multi-source data array, the maximum environmental offset rate is obtained based on the smoothed average values ​​of ambient temperature, relative humidity, air pressure, and ammonia concentration within the multi-source data array. The comprehensive deviation is then calculated by combining the target heterozygosity with the actual heterozygosity. When the comprehensive deviation exceeds the currently effective safety threshold, a warning signal is generated. The wind speed parameter is obtained by calculating the smoothed average value of airflow velocity in the physical sequence, and the smoothing suppression amount is obtained by combining it with the comprehensive deviation to adjust the operating frequency of the variable frequency motor and the step opening of the ventilation regulating valve. The actual hatching rate and the preset design hatching rate are obtained. A correction value is obtained based on the hatching deviation rate extracted from the actual hatching rate and the design hatching rate, and the safety threshold for the next breeding cycle is updated based on the correction value.

[0012] Preferably, the step of obtaining historical concentration sequences to determine concentration limits includes: during the empty house purification period before the basic duck flock to be purified is moved into the isolation environment, driving the bioaerosol detection device deployed in the isolation environment to perform air sampling, continuously obtaining the pathogen concentration quantification value in the empty house state, and defining the pathogen concentration quantification value obtained within a preset time window as the historical concentration sequence; calculating the arithmetic mean of the historical concentration sequence to obtain the noise floor mean, and calculating the standard deviation of the historical concentration sequence to obtain the noise dispersion; multiplying the noise dispersion by a preset threshold coefficient to obtain the fluctuation tolerance term, and adding the fluctuation tolerance term to the noise floor mean to obtain the concentration limit.

[0013] Preferably, the calculation of the arithmetic mean of each signal in the multidimensional physicochemical signal within a preset sliding time window to obtain the corresponding smoothed mean, and generating a physical sequence with timestamps based on the smoothed mean; generating a discrete sequence based on the biochemical detection signal and concentration limit includes: calculating the arithmetic mean of ambient temperature, relative humidity, air pressure, ammonia concentration and air velocity within a continuous sliding time window to obtain the smoothed mean of each multidimensional physicochemical signal; concatenating the smoothed mean of each multidimensional physicochemical signal within the same sampling period into a multidimensional vector and arranging them according to the chronological order of sampling time to generate a physical sequence; generating a first state value representing the microbiological safety status when the biochemical detection signal is less than the concentration limit; generating a second state value representing the microbiological risk status when the biochemical detection signal is greater than or equal to the concentration limit; and arranging the continuously generated first state value or second state value according to the chronological order of recording time to generate a discrete sequence.

[0014] Preferably, the step of obtaining the time series deviation based on the relative time difference between the physical sequence and the discrete sequence, and the pressure deviation between the current transient sample value and its smoothed mean of the pressure signal in the physical sequence, includes: subtracting the current timestamp corresponding to the current acquisition point in the physical sequence from the entry timestamp corresponding to the latest data entry in the discrete sequence and taking the absolute value to obtain the relative time difference; subtracting the current transient sample value and its smoothed mean of the pressure signal in the physical sequence and taking the absolute value to obtain the pressure deviation; performing dimensionless processing on the relative time difference using a preset time scale parameter to obtain a dimensionless relative time difference, and performing dimensionless processing on the pressure deviation using a preset pressure scale parameter to obtain a dimensionless pressure deviation; performing exponential mapping calculation with the natural constant as the base and the dimensionless pressure deviation as the exponent to obtain a pressure deviation mapping term; and multiplying the pressure deviation mapping term by the dimensionless relative time difference to obtain the time series deviation.

[0015] Preferably, the step of performing timeline matching between the physical sequence and the discrete sequence based on the time series deviation to construct a multi-source data matrix includes: determining a matching time window between the physical sequence and the discrete sequence based on the time series deviation; within the matching time window, aligning the state values ​​in the discrete sequence with the physical sequence data whose time positions are closest; filling the data gaps in the physical sequence using an effective value decay preservation method, and filling the data gaps in the discrete sequence using a state preservation method; and splicing the aligned and filled physical sequence with the discrete sequence to construct a multi-source data matrix.

[0016] Preferably, the step of obtaining the maximum environmental offset rate feature based on the smoothed average values ​​of ambient temperature, relative humidity, air pressure, and ammonia concentration within the multi-source data array, and calculating the comprehensive deviation by combining the target heterozygosity and the actual heterozygosity, includes: obtaining preset safe intervals corresponding to ambient temperature, relative humidity, air pressure, and ammonia concentration respectively; calculating the ambient temperature offset rate, relative humidity offset rate, air pressure offset rate, and ammonia concentration offset rate respectively based on the degree of deviation of the smoothed average values ​​of ambient temperature, relative humidity, air pressure, and ammonia concentration from the corresponding preset safe intervals; extracting the maximum value among the ambient temperature offset rate, relative humidity offset rate, air pressure offset rate, and ammonia concentration offset rate as the maximum environmental offset rate feature; dividing the target heterozygosity by the actual heterozygosity to obtain the heterozygosity gain ratio; and multiplying the heterozygosity gain ratio by the maximum environmental offset rate feature to obtain the comprehensive deviation.

[0017] Preferably, the wind speed parameter is obtained by calculating the smoothed mean value of air velocity in the physical sequence, and the smoothing suppression amount is obtained by combining the comprehensive deviation. This includes: performing dimensionless processing on the wind speed parameter using a preset air velocity scale parameter to obtain a dimensionless wind speed parameter; performing exponential mapping calculation with the natural constant as the base and the negative of the dimensionless wind speed parameter as the exponent to obtain a wind speed attenuation mapping term; calculating the difference between the constant 1 and the wind speed attenuation mapping term to obtain an airflow cutoff gain term; and multiplying the comprehensive deviation and the airflow cutoff gain term to obtain the smoothing suppression amount.

[0018] Preferably, the adjustment of the operating frequency of the variable frequency motor and the step opening of the ventilation regulating valve includes: determining the adjustment direction of the variable frequency motor and the ventilation regulating valve based on the abnormal physicochemical signal that causes an increase in the overall deviation; using the smoothing suppression amount as the amplitude correction coefficient to limit the single-time operating frequency change of the variable frequency motor and the single-time step opening of the ventilation regulating valve; when the ammonia concentration is higher than the corresponding preset safety range, increasing the exhaust volume after high-efficiency filtration or disinfection according to the amplitude correction coefficient; when the air velocity is higher than the corresponding preset safety range, limiting the frequency increase of the variable frequency motor or reducing the single-time opening change of the ventilation regulating valve according to the amplitude correction coefficient; when the positive pressure gradient of the isolation environment is lower than the preset positive pressure gradient, increasing the supply air volume or decreasing the exhaust air volume according to the amplitude correction coefficient.

[0019] Preferably, the step of obtaining a correction value based on the hatching deviation rate extracted from the actual hatching rate and the designed hatching rate, and updating the safety limit threshold for the next breeding cycle based on the correction value, includes: subtracting the designed hatching rate from the actual hatching rate to obtain the hatching deviation rate; performing an exponential mapping calculation with the natural constant as the base and the hatching deviation rate as the exponent to obtain a deviation mapping term; calculating the difference between the deviation mapping term and constant 1 to obtain an error residual term; multiplying a preset learning step size term and the error residual term to obtain a correction value; calculating the difference between constant 1 and the correction value to obtain an attenuation gain term; multiplying the safety limit threshold effective in the current breeding cycle with the attenuation gain term to obtain a candidate updated safety limit threshold; limiting the candidate updated safety limit threshold to between a preset lower threshold and a preset upper threshold to obtain an updated safety limit threshold; and overwriting the judgment threshold parameter in the control terminal with the updated safety limit threshold as the safety limit threshold for the next breeding cycle.

[0020] The beneficial effects of this invention are as follows: This invention integrates isolation environment purification, item sterilization, personnel purification, and animal microbial quality control into the same breeding process. By excluding individuals that do not meet preset microbiological standards, internal reproduction within the group, and continuous retesting of offspring, the purification results can be extended to subsequent breeding generations. This invention calculates the time-series deviation using relative time difference and air pressure deviation, and performs time matching on continuous multidimensional physicochemical signals and discontinuous biochemical detection signals based on the time-series deviation, reducing data misalignment, false alarms, and missed alarms caused by different sampling frequencies. This invention utilizes environmental offset characteristics and current duck flock genetic quality parameters to calculate the comprehensive deviation, enabling isolation environment early warning to adapt to the genetic quality status of duck flocks in different breeding generations, rather than relying solely on... Fixed single-item environmental thresholds; This invention calculates smoothing inhibition based on comprehensive deviation and airflow velocity, and uses the smoothing inhibition to limit the single-action amplitude of variable frequency motors and ventilation regulating valves, which can reduce the secondary stress caused by abrupt changes in wind pressure and airflow velocity to duck flocks; This invention distinguishes between microbial risk control mode and environmental deviation control mode, and switches independent filtration or disinfection ventilation branches when microbial risks occur, while maintaining the minimum number of air exchanges and positive pressure gradient, avoiding ammonia accumulation or pressure gradient failure caused by simply closing ventilation; This invention updates the safety limit threshold for the next breeding cycle based on hatching performance that meets the effectiveness conditions, and limits the amplitude through upper and lower thresholds, so that the environmental warning standard can be stably adjusted according to changes in breeding status. Attached Figure Description

[0021] Figure 1 This is a flowchart of an SPF isolation environment early warning method based on multi-source data, used in the purification and breeding method of an SPF duck closed flock according to the present invention.

[0022] Figure 2This is a functional block diagram of an SPF isolation environment early warning system based on multi-source data used in the purification and breeding method of an SPF duck closed flock of the present invention;

[0023] Figure 3 This is a comparison chart of airflow velocity parameters and actual operating load rate control strategies of variable frequency motors;

[0024] Figure 4 It is a state diagram of the process of timing alignment of multi-source heterogeneous signals and construction of multi-source data arrays;

[0025] Figure 5 It is a mapping diagram of the smoothing suppression amount driven by the combined deviation and wind speed parameters;

[0026] Figure 6 A pedigree diagram of the breeding of SPF ducks in a closed colony;

[0027] Figure 7 This is a diagram showing the phenotypic breeding results of the SPF duck closed group. Detailed Implementation

[0028] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0029] The specific embodiments of the present invention will now be described in detail with reference to the accompanying drawings.

[0030] Example 1

[0031] This invention discloses a method for purifying and cultivating an SPF duck closed colony, comprising the following steps:

[0032] The ducks to be purified were moved into an isolation environment for rearing. The isolation environment uses a high-efficiency filtration device to achieve independent ventilation, positive pressure protection and all-in-all-out.

[0033] Items entering the isolation environment are sterilized, and staff entering the isolation environment undergo changing clothes, showering, and disinfection management.

[0034] Multi-source data monitoring of the isolation environment, and adjustment of ventilation equipment in the isolation environment based on the monitoring results;

[0035] Obtain the microbiological test results of in vitro samples from the basic duck flock to be purified, isolate individuals whose microbiological test results do not meet the preset microbiological standards from the breeding flock, and do not include such individuals in the next generation of the breeding flock;

[0036] During the breeding period, no outside ducks are introduced. Individuals whose microbiological test results meet the preset microbiological standards are allowed to reproduce within the group. The offspring are then retested for microbiology. Duck flocks that meet the preset microbiological standards for consecutive generations are identified as SPF duck closed flocks.

[0037] Among them, a new SPF isolation environment early warning method based on multi-source data is proposed in the step of monitoring the isolation environment with multi-source data and adjusting the ventilation equipment in the isolation environment according to the monitoring results, so as to improve the stability and controllability of purification cultivation. To achieve this goal, specific references are made. Figure 1 The isolation and early warning method includes steps S1-S4:

[0038] S1. Acquire multidimensional physical and chemical signals and biochemical and electrical signals, generate physical sequences and discrete sequences with timestamps, and obtain the time deviation to perform time axis matching.

[0039] It should be noted that under the closed breeding conditions of duck flocks, environmental physicochemical indicators will exhibit high-frequency continuous fluctuations, while the occurrence of pathogens is a low-frequency, sporadic event. This difference in sampling frequency makes it difficult to obtain aligned data at the same moment. Furthermore, the inherent background noise in multi-source sensor networks when capturing weak signals can cause deviations in fixed threshold judgments. Therefore, this invention primarily measures the temporal misalignment of different data collections through temporal deviation, laying a solid foundation for accurate alignment of subsequent underlying data.

[0040] Specifically, during the clean period before the duck flock is moved into the isolation chamber, the bioaerosol sensor network deployed in the isolation chamber is driven to perform air sampling, continuously capture pathogen concentration values ​​in the empty chamber state, and define the pathogen concentration values ​​within the preset time window as the historical concentration sequence.

[0041] It should be further explained that, in this embodiment, the steps for obtaining the quantitative value of the pathogen concentration are as follows: The bioaerosol detection device is used to perform continuous enrichment sampling of the target pathogenic microorganisms in a unit volume of air within the isolated environment; the enriched aerosol samples are then subjected to liquid-phase elution and nucleic acid lysis extraction sequentially; subsequently, reverse transcription and real-time quantitative PCR amplification reactions are performed based on the extracted nucleic acid products, and fluorescence signals are collected simultaneously; finally, the collected fluorescence signals are converted into the corresponding pathogenic nucleic acid equivalent value using a preset standard curve. In specific industrial implementations, the bioaerosol detection device specifically adopts a fully automated online bioaerosol nucleic acid analyzer that integrates automated aerosol sampling and microfluidic PCR amplification functions.

[0042] The arithmetic mean of the historical concentration series is calculated to obtain the noise floor mean; the standard deviation of the historical concentration series is calculated to obtain the noise dispersion; and the concentration limit is obtained based on the noise dispersion, the noise floor mean, and the sensitivity coefficient.

[0043] Specifically, the concentration limit satisfies the expression:

[0044] ;

[0045] In the formula, For concentration limits, This represents the average noise level. The sensitivity coefficient, This represents the noise dispersion.

[0046] in, The larger the value, the more intense the underlying voltage level oscillations of the multi-source sensor network in its pure state, causing... The larger the value, the wider the fault tolerance buffer reserved for background noise becomes, thus making it more likely to achieve the filtering effect of suppressing false blocking actions caused by hardware interference. The smaller the value, the more stable the baseline output of the multi-source sensor network, making it more stable. The smaller the value, the closer the dynamically generated safety boundary becomes to the background signal, thus improving the sensitivity of pathogen mutation response monitoring.

[0047] It should be further clarified that the proprietary technical terms involved in the preliminary data collection and analysis of this invention are defined as follows: the pure period is the completely enclosed and absolutely empty stage in which the isolation chamber is completely closed and free from interference from living biological sources after the previous breeding batch of ducks is removed, after cleaning, disinfection and empty house rest, until the current batch of ducks is moved in; the pathogen concentration value is the equivalent concentration electrical signal output by the bioaerosol sensing network configured inside the isolation chamber after continuous sampling and quantitative analysis of target pathogenic microorganisms in a unit volume of air; the historical concentration sequence is the time series of pathogen concentration values ​​continuously recorded by the sensing network during the baseline period without interference from biological metabolism within the above-mentioned pure period.

[0048] Furthermore, the ambient temperature, relative humidity, air pressure, ammonia concentration, and air velocity signals inside the isolation chamber are collected in real time. Within a continuously set sliding time window, multiple instantaneous sample values ​​of each signal are extracted, and the arithmetic mean of each signal within the sliding time window is calculated to generate a smoothed mean of each signal. The smoothed mean values ​​of the ambient temperature, relative humidity, air pressure, ammonia concentration, and air velocity signals within the same sampling period are then spliced ​​into a multi-dimensional vector to generate a physical sequence arranged in chronological order.

[0049] Simultaneously, real-time biochemical electrical signals for target pathogen detection are acquired, and these signals are compared with concentration limits.

[0050] When the biochemical electrical signal is less than the concentration limit, a value of 0 representing a safe state is generated; when the biochemical electrical signal is greater than or equal to the concentration limit, a value of 1 representing a warning state is generated; the continuously generated values ​​of 0 or 1 are combined into a discrete sequence according to the order of their entry time.

[0051] Extract the current timestamp corresponding to the current acquisition point in the physical sequence and the entry timestamp corresponding to the latest data entry in the discrete sequence. Subtract the current timestamp from the entry timestamp and take the absolute value to obtain the relative time difference.

[0052] Extract the current transient sample value of the air pressure signal in the physical sequence, subtract it from the smoothed mean of the air pressure signal, and take the absolute value to obtain the air pressure deviation.

[0053] The time sequence deviation is obtained by measuring the air pressure deviation and the relative time difference.

[0054] Specifically, the timing deviation satisfies the expression:

[0055] ;

[0056] In the formula, This refers to the time series deviation. The relative time difference This represents the air pressure deviation.

[0057] in, The larger the value, the wider the natural time interval between the physical sequence and the discrete sequence, resulting in a greater time series deviation. The larger the value, the stronger the alignment penalty cardinality becomes for asynchronous features, thus becoming more likely to trigger interpolation calculations to smooth out order mismatches and improve alignment accuracy. The smaller the value, the greater the overlap between the physical sequence and the discrete sequence on the time axis, resulting in a higher degree of temporal deviation. The smaller the value, the weaker the alignment penalty for asynchronous features becomes, making it more likely to preserve the original temporal mapping relationship and maintain the accuracy of the judgment.

[0058] at the same time, The larger the value, the more drastic the changes in the aerodynamic field or transient pressure fluctuations inside the isolation chamber, causing... The more it tends to increase, the more it will eventually lead to... Enlarging causes the timing alignment containment window for asynchronous states to tighten accordingly, thereby making it more likely to achieve the goal of amplifying the penalty weight of the current section and rejecting mismatched data input filtering; The smaller the value, the more stable the aerodynamic environment inside the isolation chamber, making it easier to operate. The more it tends to decrease and converge to 1, the more it makes... The output falls back to the baseline, making the timing alignment tolerance for asynchronous states more relaxed, thus achieving the control effect of widening the timing alignment tolerance window to maintain the continuity of data alignment operations.

[0059] For example, the preset time window determines the completeness of the coverage of day-night temperature drift and background environmental fluctuations by the noise baseline extraction. The empirical range is [24, 72] hours, and it is set to 48 hours in this embodiment. The implementers can adjust it according to the length of the isolation shed's rest period.

[0060] For example, the preset sliding time window has an empirical range of [5, 30] minutes, and is set to 15 minutes in this embodiment. Implementers can adjust this parameter according to the inherent sampling frequency of various sensors and the severity of dynamic fluctuations in the isolation environment.

[0061] For example, the empirical range of the sensitivity coefficient is [1.5, 3.5], and it is set to 2.5 in this embodiment. Implementers can set the sensitivity coefficient according to the level of external pathogen exposure risk faced by the duck flock during the breeding season. For example, when an epidemic occurs in the surrounding area or during the peak winter and spring season, the sensitivity coefficient can be appropriately increased to enhance the initial high-sensitivity prevention capability; when the isolation chamber has an extremely high level of closed protection and the sensor network is susceptible to non-biological interference such as dust, the sensitivity coefficient can be appropriately reduced to suppress false blocking actions.

[0062] For example, the empirical range for the average noise floor is [5, 15] micrograms per milliliter, and in this embodiment it is set to 10 micrograms per milliliter. Implementers can set the average noise floor based on the hardware's own temperature drift characteristics or dark current reference. For instance, when sensor element aging causes a significant increase in zero-point dark current, the average noise floor can be appropriately increased to prevent inherent noise from causing system misjudgments. When using a high-precision, low-noise sensor chip and the space is extremely clean, the average noise floor can be appropriately reduced to improve the extraction accuracy of initial signals at extremely low concentrations.

[0063] S2. Based on the maximum amplitude characteristics of the ambient temperature signal, relative humidity signal, air pressure signal, and ammonia concentration signal after extreme value normalization processing, and combined with the alignment results of the aforementioned output, calculate the comprehensive deviation at the current moment.

[0064] It should be noted that in closed breeding conditions, even minor environmental degradation can cause chronic physiological damage to duck flocks. Conventional monitoring methods often rely solely on environmental data, lacking consideration of the tolerance characteristics of the controlled species. Therefore, this invention, based on the temporal alignment of the underlying data, further measures the superposition of environmental drift and genetic characteristics through comprehensive deviation, weakening the one-sided interference of random fluctuations from a single sensor, and providing a unique basis for subsequent triggering of global early warning and wind speed cutoff adjustment.

[0065] Specifically, based on the time-series deviation obtained from the assessment, time-axis matching is performed on the physical sequence and the discrete sequence, and the two are aligned and stitched together to construct a multi-source data matrix;

[0066] When the latest data bit of the multi-source data array undergoes a step change from a value of 0 to a value of 1, an early warning signal is generated and the action of closing the ventilation regulating valve of the isolation chamber is triggered.

[0067] In a normal production scenario where no step change is detected, the smoothed mean values ​​of the ambient temperature signal, relative humidity signal, air pressure signal, and ammonia concentration signal extracted from the multi-source data array are subjected to extreme value normalization algorithms. The term with the largest value in the normalization result is defined as the maximum environmental offset feature.

[0068] Obtain the batch identifier of the current duck flock, retrieve the preset standard target heterozygosity based on the batch identifier, and simultaneously retrieve the actual heterozygosity of the duck flock that has been pre-determined and recorded by gene sequencing.

[0069] The overall deviation is obtained based on the target hybridity, the actual hybridity, and the maximum environmental offset rate.

[0070] Specifically, the overall deviation satisfies the expression:

[0071] ;

[0072] In the formula, For the overall deviation, For the target heterozygosity, For actual heterozygosity, This represents the maximum environmental offset rate.

[0073] in, The larger the value, the greater the deviation of a certain core physicochemical indicator inside the isolation chamber from the safety baseline, making... The greater the increase, the more likely it is to trigger an early entry into emergency control levels to protect the safety of the duck flock. The smaller the value, the more stable and safe the aerodynamic parameters inside the isolation chamber are, ensuring that... The more it decreases, the more it tends to achieve the goal of extending the regulation sleep cycle to reduce the overall energy consumption of the underlying execution hardware.

[0074] also, The smaller the value, the higher the degree of inbreeding and degeneration within the duck flock. The more it tends to increase, the more it makes The larger the increase, the more likely it is to raise the protection threshold for vulnerable duck flocks in order to implement precise intervention; The larger the size, the richer the genetic diversity of the duck population. The more it tends to decrease, the more it makes The more it tends to decrease, the more likely it is to achieve the control effect of relaxing the tolerance for sudden physical and chemical degradation response to improve the environmental drift tolerance.

[0075] It should be noted that the proprietary technical terms involved in the multi-source data deep fusion analysis of this invention are defined as follows: the maximum environmental offset rate characteristic is the deviation value of the multidimensional physicochemical signal with the largest value after the smoothed mean values ​​of the environmental temperature signal, relative humidity signal, air pressure signal, and ammonia concentration signal within the multi-source data array are respectively subjected to the extreme value normalization algorithm; the target heterozygosity is the health standard baseline value of the population allele heterozygosity status set according to the poultry breeding and preservation technical specifications to maintain the basic stress resistance and population vitality of the current closed breed duck flock; the actual heterozygosity is the true heterozygosity ratio of the current batch of duck flock samples after specific genotyping determination and transfer recording by gene sequencing analysis equipment.

[0076] In this embodiment, the actual heterozygosity obtained is 0.6. For example, the empirical range of the target heterozygosity is [0.5, 0.9], and in this embodiment it is set to 0.7. Implementers can set the target heterozygosity based on the breeding generation and basic stress resistance baseline of the duck strain. For example, when facing sensitive inbred breeding lines with high conservation value, the target heterozygosity can be appropriately increased to enhance the protective sensitivity of vulnerable lines, ensuring that even minor environmental drifts trigger high-level defenses. When the breeding population is a commercial population with strong stress resistance and extremely low degradation risk, the target heterozygosity can be appropriately reduced to relax the degradation response tolerance and improve the energy-saving effect of the environmental control dormancy cycle.

[0077] S3. In response to the calculated comprehensive deviation exceeding the limit, trigger the on-site early warning action and calculate the smoothing suppression amount of the device based on the moving average characteristics of the airflow signal.

[0078] It should be noted that in the air circulation control of the isolation chamber, the fans and regulating valves have inherent mechanical inertia. If a large adjustment command is issued instantaneously when an anomaly is detected, it can easily lead to a sudden increase in air pressure and turbulent airflow within the chamber. This sudden change in airflow can severely stimulate the sensitive duck flock, causing secondary stress injuries. Therefore, this invention, while triggering the early warning defense line based on the comprehensive deviation, simultaneously introduces a smoothing suppression amount to measure the instantaneous airflow impact of the adjustment command on the wind field, providing a safety limit for the stable control of the underlying equipment.

[0079] It should be noted that the proprietary technical terms involved in the smooth frequency reduction flexible control of this invention are defined as follows: the wind speed parameter is the smoothed flow velocity characteristic value obtained by tracing back to the historical time axis and extracting the air flow velocity signal from the multi-source data array and performing an arithmetic average operation within a set continuous sliding time window; the smoothing suppression amount is the amplitude adjustment correction coefficient used to constrain the mechanical output amplitude of the underlying environmental control hardware, which is generated by multiplying the global comprehensive deviation value under the current cross section with the dynamic attenuation characteristic of the air flow velocity signal.

[0080] Specifically, the current safety limit threshold that is effective within the current control cycle of the system is obtained. The safety limit threshold is taken as the preset initial judgment benchmark when the first breeding batch of the system is run, and the safety limit threshold after cross-cycle feedback update is taken in subsequent breeding batches.

[0081] When the overall deviation exceeds the safety limit threshold, an early warning signal containing the current overall deviation and abnormal physical and chemical indicators is generated, driving the audible and visual alarm device configured on the field industrial control node to perform high-frequency flashing and whistling early warning actions, and simultaneously pushing the early warning signal to the host computer monitoring terminal for interface alarm display;

[0082] While performing the above-mentioned warning actions, the smoothed mean of the air velocity signal of the physical sequence within the set time window is traced back to the historical time axis and used as the smoothed wind speed parameter.

[0083] The smoothing suppression amount is obtained based on the wind speed parameter and the overall deviation.

[0084] Specifically, the smoothing suppression amount satisfies the expression:

[0085] ;

[0086] In the formula, To smooth out the suppression amount, For the overall deviation, This refers to the wind speed parameter.

[0087] in, The larger the value, the stronger the mechanically driven gas movement or pipeline air pressure inside the isolation chamber, making... The more it tends to increase, the more it leads to The greater the pressure, the stronger the force that compels the fan to reduce its frequency and the valve to open in stages to calm the high-pressure cyclone inside the shed, thus preventing secondary stress injuries to the duck flock caused by sudden changes in wind pressure. The smaller the value, the more it indicates that the airflow replacement rate inside the isolation chamber is maintained in a slow-reducing state, thus... As it decreases and approaches zero, the final amount of smoothing inhibition decreases. The synchronous convergence decreases, thus becoming more and more likely to allow the conventional environmental control regulation loop to perform fine-tuning actions without hindrance, thereby avoiding the maintenance purpose of large-scale control dead zones caused by regulation overload.

[0088] at the same time, The larger the value, the closer the overall physical and chemical deterioration of the space is to the environmental safety limit. The larger the value, the more likely it is to reach the constraint target of forcibly tightening the instantaneous step opening change tolerance of the underlying aerodynamic field control hardware, thereby defending against sudden mechanical impacts. The smaller the value, the more stable the current environmental health indicators are at normal levels. The more it decreases, the closer it gets to the normal control state where the limiting clamp on the wind farm actuators is released, enabling the underlying equipment to respond efficiently to regular climate regulation commands.

[0089] Furthermore, the smoothing suppression amount is used as the amplitude correction coefficient to adjust the operating frequency of the variable frequency motor and the step opening of the ventilation regulating valve, so as to achieve smooth and flexible control of the airflow in the warehouse without sudden changes.

[0090] For example, the empirical range of the safety threshold is [2, 5]. In this embodiment, the initial safety threshold for the first breeding batch is set to 3. Implementers can set the initial safety threshold based on the age and stress vulnerability of the live ducks in the isolation chamber. For instance, when the breeding chamber is in the stage of young ducklings prone to death due to environmental disturbances, the initial safety threshold can be appropriately lowered to activate early agile protection; when the breeding batch belongs to the stage of adult breeding ducks with strong resistance and tolerance, the initial safety threshold can be appropriately increased to reduce the frequency of alarms caused by occasional physical and chemical disturbances and avoid mechanical fatigue of the environmental control actuator.

[0091] S4. Integrate the batch hatching rate characteristics of the number of hatched chicks and the number of qualified hatching eggs, and reverse generate correction values ​​for adjusting the trigger threshold of early warning actions.

[0092] It should be noted that in long-term operation, a fixed warning threshold often fails to adapt to the microscopic decline in the tolerance of the controlled species. Without a feedback and update mechanism, when the duck flock experiences a decline in its resilience due to intergenerational accumulation, the environmental control terminal may stubbornly maintain the original high threshold without triggering an alarm, leading to delayed warning actions and a decrease in the breeding success rate. Therefore, this step, as the advanced macro-control loop of the entire method, measures the feedback deviation of batch-level hatching data on the control cycle, and autonomously lowers the safety threshold for triggering warning actions, forming an increasingly sensitive and realistic dynamic warning closed-loop mechanism.

[0093] Specifically, after the breeding batch is completed, the number of chicks and the number of qualified hatching eggs carried over from that batch of duck flock are obtained. The number of chicks is divided by the number of qualified hatching eggs to obtain the actual hatching rate of that batch. The preset design hatching rate is obtained, and the design hatching rate is subtracted from the actual hatching rate to obtain the hatching deviation rate.

[0094] The correction value is obtained based on the incubation deviation rate and the preset learning step size.

[0095] Based on the correction value and the current safety limit threshold effective in the current control cycle, the updated safety limit threshold for parameter overwriting in the next control cycle is obtained, and the judgment threshold parameter at the bottom layer of the system control terminal is overwritten with the updated safety limit threshold. Thus, when the actual hatching rate is lower than the design target, the proportion of the red line for triggering the early warning action is reduced, so that even a small environmental degradation in the next cycle can quickly exceed the limit and trigger the audible and visual early warning action, realizing cross-cycle adaptive iteration of highly sensitive early warning.

[0096] Specifically, the correction value and the updated safety limit threshold satisfy the following relationship:

[0097] ;

[0098] In the formula, This is the correction value for the deviation. To learn the step length item; Calculate the carry-over hatching deviation rate for the current breeding batch;

[0099] ;

[0100] In the formula, The current safety limit threshold that is in effect during the current control cycle; This is the updated safety limit threshold used for parameter overwriting in the next control cycle.

[0101] Among them, when The larger the value, the lower the actual hatching rate is compared to the designed hatching rate, resulting in a higher error residual term. The more it tends to increase, in Under constraints, make the correction value The corresponding increase makes The more it tends to decrease, the more it passes the current safety threshold. Multiplication updates the safety limit threshold. The more it tends to decrease, the more it tends to achieve an adaptive control effect that forces the judgment defense line to actively tighten towards the high-sensitivity defense zone in order to improve the initial interception performance;

[0102] when The smaller the value, the closer the actual hatching rate is to or higher than the designed hatching rate. The more it tends to decrease and become negative, the more the correction value... Synchronization tends to zero or becomes negative, making The closer it gets to or greater than 1, the more it will cause... The more stable or moderately rising the parameters are, the more likely they are to reach the stability goal of locking the judgment boundary to prevent random environmental fluctuations from causing high-frequency blind adjustments to the parameters.

[0103] also, The larger the value, the greater the residual of biological symptom decline carried over from the previous production stage, making it more likely to cause further decline. The more it tends to decrease, the more it makes The more it tends to decrease, the more it tends to lower the warning threshold to ensure a sensitive protection effect against minor degradation that can quickly exceed the limit; The smaller the value, the better the biological function of the production system is maintained, making it more efficient. As it tends to increase and approach 1, it makes... The closer they are to each other and equivalent to This causes the underlying parameter update algorithm to enter a dormant, locked state, thereby increasing the likelihood of achieving the goal of enhancing the stability of the timing continuity of the entire industrial control network.

[0104] It should be further explained that the proprietary technical terms involved in the cross-cycle adaptive parameter evolution of this invention are defined as follows: the actual hatching rate is the percentage of the total number of healthy chicks hatched after the completion of the entire production process of the current single breeding batch, calculated as the final carry-over of the controlled duck flock in the isolation shed, relative to the total number of qualified hatching eggs; the designed hatching rate is the expected hatching rate control benchmark index preset based on the production average level of the current closed strain's long-term historical breeding big data; the hatching deviation rate is the algebraic deviation residual term generated after subtracting the actual hatching rate from the designed hatching rate benchmark value; and the correction value is a dimensionless feedback parameter generated by performing nonlinear exponential mapping and step size scaling on the hatching deviation rate, used to perform proportional gain reduction on the safety limit threshold at the beginning of the next control cycle.

[0105] In this embodiment, the actual hatching rate obtained is 88%.

[0106] For example, the empirical value range of the learning step size term is [0.01, 0.1], and it is set to 0.05 in this embodiment. Implementers can adjust this parameter according to the intensity of random fluctuations in the actual hatching rate of multiple batches. For example, when the hatching results show high-frequency irregular fluctuations due to seasonal changes, the learning step size term can be appropriately reduced to increase the damping of the adaptive loop and avoid the oscillation and divergence of the control red line. When the trend of the decline in the genetic tolerance of the strain across generations is clear and the current batch experiences a severe decline in production, the learning step size term can be appropriately increased to accelerate the downward adjustment of the response and force the judgment defense line to tighten rapidly across cycles.

[0107] For example, the empirical range of the designed hatching rate is [80%, 95%], and it is set to 90% in this embodiment. Implementers can adjust it according to the long-term historical average production level of the base or refined benefit indicators. For example, when implementing a strict zero-defect quality and benefit assessment and maximizing the breeding advantages of the strain, the designed hatching rate can be appropriately increased to enhance the residual mapping base and force the safety limit threshold to always remain in a high-sensitivity protection state. When the basic hatching indicators of the base experience a global force majeure decline due to seasonal extreme heat stress, the designed hatching rate can be appropriately reduced to accelerate the residual convergence of the feedback loop and maintain baseline consistency.

[0108] For example, Figure 3 This figure compares the control strategies for airflow velocity parameters and the actual operating load rate of the variable frequency motor. As shown, the existing technology uses a step control strategy, whose control commands are greatly affected by threshold judgments. When the operating load rate crosses the judgment boundary, it jumps and directly reaches the physical upper limit of the inverter's full load, which can easily cause mechanical shock to the actuator and fluctuations in environmental parameters. In contrast, this invention constructs a buffer mechanism by extracting smoothing suppression, ensuring that the actual operating load rate of the variable frequency motor maintains a continuous and smooth transition with changes in airflow velocity parameters, avoiding truncation and jumps in command output. Operational results show that this invention effectively balances exhaust regulation needs with mechanical operational stability, suppresses command jumps across operating conditions, and enhances the flexible control of the underlying fan equipment while ensuring a smooth transition of environmental parameters.

[0109] For example, Figure 4 This diagram illustrates the state of the process for time-series alignment of heterogeneous signals and construction of a multi-source data array. As shown, the multi-dimensional physicochemical signals such as environmental temperature, humidity, and air pressure in the timestamped physical sequence differ from the discrete sequence in terms of sampling frequency and time point misalignment, resulting in gaps in the underlying data. This invention, after performing time-axis matching, employs an effective value attenuation preservation strategy for the data gaps in the multi-dimensional physicochemical signals and performs state-preserving operations on the discrete sequence to fill the gaps. The processing results demonstrate that this invention effectively overcomes the data asynchrony problem of heterogeneous sensors, eliminates abnormal gaps in the characteristic mapping amplitude of the multi-source data array, and improves the data reliability of subsequent comprehensive deviation assessment while ensuring the coherence and consistency of the underlying logic.

[0110] For example, Figure 5This is a mapping diagram of the smoothing suppression amount driven by the combined deviation and wind speed parameters. As shown in the figure, the output load of the system actuator is constrained by the combined environmental conditions. The smoothing suppression amount increases non-linearly with the changes in the combined deviation and wind speed parameters, resulting in a smooth control surface transition without local jumps or drops. Operational results demonstrate that this invention, based on multi-parameter fusion to construct an amplitude correction mechanism, limits the action step amplitude of the underlying adjustment commands, suppresses command saturation under fluctuating operating conditions, and achieves smooth operation of the variable frequency motor and ventilation control valves while ensuring environmental correction effects.

[0111] Example 2

[0112] In one embodiment, using the Muscovy duck as the base breeding material, the complete application process of the purification and breeding method for an SPF duck closed colony described in this invention is illustrated. The following descriptions of artificial domestication, family construction, and cyclic mating are further implementation methods and do not constitute all the necessary conditions for implementing the isolation purification, multi-source data monitoring, and microbiological evaluation steps described in claim 1.

[0113] I. Artificial domestication of basic duck flocks

[0114] The Mountain Duck is characterized by its small size, early onset of egg production, high egg yield, high feed conversion rate, and strong adaptability, making it suitable as the basic breeding material in this embodiment.

[0115] The initial acclimatization process for the duck flock involved a gradual transition from hydroponics to dryland rearing. In the early stages, shallow water areas were maintained, but the daily surface activity time was limited to gradually reduce the ducks' dependence on open water. During the middle stages, open water was eliminated, and intermittent spraying was used to meet the ducks' needs for feather cleaning and normal behavior. In the later stages, dryland rearing was implemented throughout, and the ducks were gradually moved to an isolated environment with high-efficiency filtration, independent ventilation, and positive pressure protection.

[0116] In terms of breeding methods, the basic duck flock is successively subjected to a free-range adaptation period, a cage-raising transition period, and a quarantine period. During the free-range adaptation period, the ducks adapt to artificial lighting, timed feeding, and human management. During the cage-raising transition period, cages that allow the ducks to stand, turn, stretch their necks, eat, and drink normally are used, and each duck is numbered to record egg production, feeding, and feces. During the quarantine period, the ducks are moved into a high-efficiency filtered isolation environment with independent ventilation, positive pressure protection, and all-in, all-out.

[0117] Regarding drinking methods, in the early stages of domestication, both open and closed drinking devices were used to guide the ducks to adapt to closed drinking; in the middle stages of domestication, open drinking was gradually reduced; and in the later stages of domestication, closed drinking devices were used for water supply.

[0118] Regarding egg-laying behavior, both open nests and closed egg boxes were set up in the early stage of domestication; in the middle stage of domestication, the open nests were gradually removed, and only the closed egg boxes were retained to guide the duck flock to form fixed-point egg-laying behavior.

[0119] After domestication, the ducks can adapt to dry-rearing, cage-rearing, and isolation environments without open water surfaces. They can drink water through closed drinking devices and develop the behavior of laying eggs at fixed points in closed egg boxes.

[0120] II. Isolation Environment Configuration and Purification Control

[0121] The acclimatized and preliminarily screened duck flocks were transferred to isolation pens for rearing. The isolation pens employed high-efficiency filtration systems to achieve independent ventilation, positive pressure protection, and all-in / all-out ventilation. The ambient temperature within the isolation pens was controlled between 16 and 35°C, relative humidity between 30% and 70%, with a maximum daily temperature difference not exceeding 4°C, a pressure gradient not less than 50 Pa, ammonia concentration not exceeding 14 mg / m³, airflow velocity at the cages not exceeding 0.20 m / s, air exchanges not less than 20 times, working illumination not less than 150 Lx, and animal illumination 5–10 Lx. Relevant environmental monitoring indicators are shown in Table 1.

[0122] Table 1. Monitoring Indicators for the Isolation Environment

[0123]

[0124] Items entering the isolation chamber are sterilized using either high-pressure steam sterilization, ultraviolet disinfection, or ozone sterilization, depending on their material. The high-pressure steam sterilization parameters are 121℃, 0.1MPa, and 30min; the ultraviolet disinfection irradiation time is no less than 15min; and the ozone concentration for ozone sterilization is no less than 60mg / m³, with a sterilization time of no less than 30min.

[0125] Staff entering the isolation ward must undergo changing clothes, shoes, showering, and disinfection procedures. Unauthorized personnel are not allowed to enter the clean area.

[0126] III. Multi-source data environment early warning and equipment adjustment

[0127] The isolation chamber is equipped with a bioaerosol detection device, an ambient temperature sensor, a relative humidity sensor, an air pressure sensor, an ammonia concentration sensor, and an air flow rate sensor, as well as a variable frequency motor and ventilation regulating valves.

[0128] According to the method in the foregoing embodiments of the present invention, the historical concentration sequence within the cleanliness cycle of the air ventilator is obtained to determine the concentration limit; biochemical detection signals and multi-dimensional physicochemical signals including ambient temperature, relative humidity, air pressure, ammonia concentration and air velocity are collected; the smoothed mean is calculated to generate a physical sequence with timestamps; a discrete sequence is generated based on the biochemical detection signals and concentration limits; the time series deviation is obtained and a multi-source data array is constructed accordingly; the comprehensive deviation is calculated in the case where no step change is detected; when the comprehensive deviation is greater than the currently effective safety limit threshold, an early warning signal is generated, and the smoothing suppression amount is obtained to adjust the variable frequency motor and ventilation regulating valve.

[0129] Specifically, when the ammonia concentration exceeds the preset safe range, the exhaust volume after high-efficiency filtration or disinfection is increased; when the air velocity exceeds the preset safe range, the frequency increase of the variable frequency motor is limited or the single opening change of the ventilation regulating valve is reduced; when the positive pressure gradient is below the preset lower limit, the supply air volume is increased or the exhaust air volume is reduced. By using the smoothing suppression amount as the amplitude correction coefficient, drastic step changes in wind speed and air pressure within the isolation chamber are avoided, thereby mitigating the stress impact on the duck flock.

[0130] IV. Microbiological evaluation, breeding access and maintenance of closed colonies

[0131] Regularly collect at least one of the following in vitro tests from the ducks in the isolation pens: blood, throat swabs, cloacal swabs, feces, or body surface samples. Tests may include Salmonella, Pasteurella multocida, Riemerella anatipestifer, avian influenza virus, Newcastle disease virus, duck enteritis virus, duck hepatitis A virus, duck reovirus, Tembusu virus, duck circovirus, avian reticuloendotheliosis virus, egg drop syndrome virus, duck parvovirus, Chlamydia, ectoparasites, and coccidia. Refer to Table 2 for specific test items and methods.

[0132] Table 2 Microbiological Detection Items and Methods for SPF Ducks

[0133]

[0134] For bacterial tests, pathogen isolation, culture, or nucleic acid detection can be used; for viral tests, nucleic acid detection, antigen detection, or serological detection can be used; for chlamydia, nucleic acid detection or other suitable in vitro detection methods can be used; for ectoparasites and coccidia, surface examination, fecal microscopy, or other parasitological detection methods can be used.

[0135] The determination of whether an individual meets the preset microbiological standards is based on the test results corresponding to different test items, rather than uniformly stating that all test items are negative for both antigens and antibodies.

[0136] Individuals whose test results do not meet the preset microbiological standards are excluded from the breeding group. Individuals that meet the preset microbiological standards are bred within the group, and their offspring are retested. In actual breeding management, individuals excluded from the breeding group may be subject to isolation observation, drug treatment, immunization, or culling, depending on the pathogen type and breeding management standards. Such treatments do not constitute the necessary technical steps for obtaining the SPF duck closed flock of this invention.

[0137] No new ducks should be introduced during the breeding period. After the duck flock meets the preset microbiological standards for consecutive generations, it can be identified as a qualified microbiologically purified basic flock and transferred to the routine monitoring stage; thereafter, the frequency of routine testing can be once per quarter, or adjusted according to the surrounding epidemiological risk and the operation of the isolation environment.

[0138] V. Further family structure and cyclical mating

[0139] In a further embodiment, to maintain the genetic quality and reproductive performance of the closed colony, multiple families can be established by selecting breeding ducks from a microbiologically purified and qualified base colony.

[0140] Artificial insemination can be used during the family breeding stage to maintain a clear paternal lineage. After collecting and incubating the eggs, the ducklings are individually tagged, and the paternal and maternal parent numbers are recorded.

[0141] Egg production, single egg weight, eggshell strength, survival rate, and hatchability were measured for each family. The measurement items, trait types, and methods for each trait are detailed in Table 3. Continuous recording and statistical analysis of these traits provide fundamental data support for inter- and intra-family evaluations.

[0142] Table 3. Selection Trait Determination Table

[0143]

[0144] Interfamilial and intrafamilial evaluations were conducted based on phenotypic and pedigree data from each family. Statistical analysis software or breeding value assessment software could be used to process the phenotypic and pedigree data, retaining families with high overall evaluations and satisfactory microbiological testing results from the initial families.

[0145] After the family line was established, closed breeding was carried out without introducing any outside ducks. From the offspring of the preserved families, selection was carried out at the hatching, rearing, and breeding entry stages, and individuals for propagation were determined based on individual phenotype, health status, microbiological test results, and reproductive capacity.

[0146] Each family lineage constructs the next generation through cyclical mating to reduce the probability of inbreeding. The cyclical mating relationships can be found in Table 4.

[0147] Table 4. Closed-group cyclic mating relationship table

[0148]

[0149] For example, Figure 6 This is a pedigree diagram of the breeding of SPF ducks in a closed colony. Figure 6 This illustrates the pedigree evolution of a base population formed using Muscovy ducks as the primary breeding material, through artificial domestication, microbiological purification, family construction, family selection, and continuous generation in a closed colony. Specifically, after the base population passed microbiological purification, multiple families were established. Egg production, single egg weight, eggshell strength, survival rate, and hatchability were measured in each family. Families with superior overall performance and passing microbiological tests were retained through inter- and intra-familial evaluations. Based on these retained families, the next generation of families was constructed through a cyclical mating process without introducing foreign ducks. Microbiological testing, individual markers, breeding entry screening, and offspring retesting were continuously implemented across each breeding generation to maintain the genetic quality and microbiological stability of the population. Thus, Figure 6 It can reflect the overall breeding path from the basic group to the closed SPF duck group and the closed breeding organization relationship in this embodiment.

[0150] Microbiological testing, individual markers, breeding records, genetic quality testing, and multi-source data environmental early warning are continuously implemented in each breeding generation. Passage hatching, offspring retesting, breeding admission screening, in-group mating, and next generation construction are repeated until a stable and maintainable SPF duck closed flock is obtained.

[0151] For example, Figure 7 This image shows the phenotypic results of a closed colony of SPF ducks. The left side of the image represents male individuals, and the right side represents female individuals. Male individuals have broad chests and backs, relatively long bodies, yellowish-green beaks, brown irises, orange-red shanks and webs, grayish-brown back feathers, and grayish-white belly feathers. Female individuals have slender, well-proportioned, and elegant bodies, yellow beaks, brown irises, orange-red shanks and webs, and light mottled feathers. These SPF ducks exhibit rapid growth and development, high survival rates, significant weight gain during the brooding period, high egg production rates and long laying periods, robust physiques, symmetrical development, and stable reproductive performance.

[0152] VI. Evaluation of Cultivation Results and Cross-Period Feedback

[0153] Through the above-mentioned artificial domestication, isolation and purification, breeding access control, internal reproduction within the group, continuous re-examination of offspring, and further cyclic mating, SPF duck closed flocks with clear microbial quality, relatively stable genetic quality, and the ability to be maintained in an isolated environment for a long time can be obtained.

[0154] After completing a breeding batch, the following statistics are recorded: the number of qualified hatching eggs, the number of fertilized eggs, the number of healthy chicks hatched, the egg storage time, and the incubation equipment operation records. First, it is determined whether the fertilization rate is within the preset range, whether the egg storage time exceeds the preset period, and whether there are any abnormalities in temperature, humidity, or egg turning in the incubation equipment. If there is a significant decrease in fertilization rate, excessive egg storage time, or incubation equipment malfunction, the hatching result of this batch is considered to be mainly affected by non-isolated environmental factors, and the hatching rate of this batch will not be used to update the safety threshold. If all the above parameters are within the preset effective range, the actual hatching rate is calculated based on the number of healthy chicks hatched and the number of qualified hatching eggs, and the deviation between the actual hatching rate and the designed hatching rate is used to update the safety threshold for the next breeding cycle.

[0155] The updated safety thresholds are stored in the isolation environment early warning control device and take effect at the start of the next breeding batch. This effectiveness assessment reduces erroneous feedback caused by egg fertilization status, egg storage conditions, or incubation equipment malfunctions affecting environmental threshold updates.

[0156] Example 3

[0157] In the purification and breeding method of SPF duck closed flock described in this invention, the multi-source data monitoring, early warning and equipment adjustment process can be realized by an SPF isolation environment early warning system based on multi-source data.

[0158] For example, Figure 2 This diagram illustrates the functional modules of a multi-source data-based SPF isolation environment early warning system used in the purification and breeding method for SPF duck colonies according to the present invention, encompassing multi-source data monitoring, early warning, and equipment adjustment processes. This multi-source data-based SPF isolation environment early warning system can be installed in an electronic device. Depending on the functions implemented, the multi-source data-based SPF isolation environment early warning system 100 may include: a bottom-level data acquisition module 101, a sequence generation and alignment module 102, a time-series alignment and array construction module 103, a comprehensive deviation evaluation module 104, an early warning and flexible suppression module 105, and a cross-period adaptive module 106.

[0159] The module described in this invention, also referred to as a unit, refers to a series of computer program segments that can be executed by the processor of an electronic device and perform a fixed function, and are stored in the memory of the electronic device. In this embodiment, the functions of each module / unit are as follows:

[0160] The underlying data acquisition module 101 is used to acquire historical concentration sequences to determine concentration limits, collect biochemical electrical signals, and multidimensional physicochemical signals including ambient temperature, relative humidity, air pressure, ammonia concentration, and air velocity.

[0161] The sequence generation and comparison module 102 is used to calculate the arithmetic mean of each item in the multidimensional physicochemical signal within a preset sliding time window to obtain a smoothed mean, and generate a physical sequence with timestamps accordingly; and generate a discrete sequence based on the biochemical electrical signal and concentration limit.

[0162] The time-series alignment and array building module 103 is used to obtain the time-series deviation degree based on the relative time difference between the physical sequence and the discrete sequence, and the pressure deviation between the current transient sample value and its smoothed mean value of the pressure signal in the physical sequence, and to perform time axis matching on the physical sequence and the discrete sequence accordingly to build a multi-source data array.

[0163] The comprehensive deviation assessment module 104 is used to obtain the maximum environmental deviation rate characteristics based on the smoothed average values ​​of ambient temperature, relative humidity, air pressure and ammonia concentration in the array when no step change is detected in the multi-source data array, and to calculate the comprehensive deviation by combining the obtained target heterozygosity and the actual heterozygosity.

[0164] The warning and flexible suppression module 105 is used to generate a warning signal when the comprehensive deviation is greater than the currently effective safety limit threshold, and to calculate the smoothed average value of the air velocity in the physical sequence to obtain the wind speed parameter, and to obtain the smoothing suppression amount in combination with the comprehensive deviation, so as to adjust the variable frequency motor and the ventilation regulating valve.

[0165] The cross-cycle adaptive module 106 is used to obtain the actual hatching rate and the preset design hatching rate, obtain the correction value based on the hatching deviation rate extracted from the actual hatching rate and the design hatching rate, and update the safety limit threshold of the next control cycle using the correction value.

[0166] In the several embodiments provided by this invention, it should be understood that the disclosed methods and systems can be implemented in other ways. For example, the system embodiments described above are merely illustrative; for instance, the division of modules is only a logical functional division, and other division methods may be used in actual implementation.

[0167] The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0168] Furthermore, the functional modules in the various embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or in a combination of hardware and software functional modules. It will be apparent to those skilled in the art that the present invention is not limited to the details of the above exemplary embodiments, and that the present invention can be implemented in other specific forms without departing from its spirit or essential characteristics.

[0169] The system also includes other components well known to those skilled in the art, such as communication buses and communication interfaces, the settings and functions of which are known in the art and will not be described in detail here.

[0170] Although preferred embodiments of the present application have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of this application.

[0171] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.

[0172] The above provides a detailed description of the purification and cultivation method for SPF duck closed colonies provided in this application. Specific examples have been used to illustrate the principles and implementation methods of this application. The description of the above embodiments is only for the purpose of helping to understand the method and core ideas of this application. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application areas based on the ideas of this application. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for purifying and cultivating an SPF duck closed flock, characterized in that, include: The ducks to be purified were moved into an isolation environment for rearing. The isolation environment uses a high-efficiency filtration device to achieve independent ventilation, positive pressure protection and all-in-all-out. Items entering the isolation environment are sterilized, and staff entering the isolation environment undergo changing clothes, showering, and disinfection management. Multi-source data monitoring of the isolation environment, and adjustment of ventilation equipment in the isolation environment based on the monitoring results; Obtain the microbiological test results of in vitro samples from the basic duck flock to be purified, isolate individuals whose microbiological test results do not meet the preset microbiological standards from the breeding flock, and do not include such individuals in the next generation of the breeding flock; During the breeding period, no outside ducks are introduced. Individuals whose microbiological test results meet the preset microbiological standards are allowed to reproduce within the group. The offspring are then retested for microbiology. Duck flocks that meet the preset microbiological standards for consecutive generations are identified as SPF duck closed flocks.

2. The method for purifying and cultivating an SPF duck closed flock according to claim 1, characterized in that, The method of performing multi-source data monitoring of the isolation environment and adjusting the ventilation equipment in the isolation environment based on the monitoring results includes: Historical concentration sequences were obtained to determine concentration limits, and biochemical detection signals and multidimensional physicochemical signals including ambient temperature, relative humidity, air pressure, ammonia concentration and air velocity were collected. Calculate the arithmetic mean of each signal in the multidimensional physicochemical signal within a preset sliding time window, obtain the corresponding smoothed mean, and generate a physical sequence with timestamps based on the smoothed mean; generate a discrete sequence based on the biochemical detection signal and concentration limit. The time series deviation is obtained based on the relative time difference between the physical sequence and the discrete sequence, and the pressure deviation between the current transient sample value and its smoothed mean in the physical sequence; the time series deviation is used to perform time axis matching between the physical sequence and the discrete sequence to construct a multi-source data matrix. In the absence of detected step abrupt changes in the multi-source data array, the maximum environmental offset rate characteristics are obtained based on the smoothed average values ​​of ambient temperature, relative humidity, air pressure, and ammonia concentration within the multi-source data array, and the comprehensive deviation is calculated by combining the target heterozygosity with the actual heterozygosity. When the overall deviation exceeds the currently effective safety limit threshold, an early warning signal is generated. The wind speed parameter is obtained by calculating the smoothed mean value of the airflow velocity in the physical sequence and combining it with the overall deviation to obtain the smoothing suppression amount, so as to adjust the operating frequency of the variable frequency motor and the step opening of the ventilation regulating valve. Obtain the actual hatching rate and the preset design hatching rate, obtain the correction value based on the hatching deviation rate extracted from the actual hatching rate and the design hatching rate, and update the safety limit threshold for the next breeding cycle based on the correction value.

3. The method for purifying and cultivating an SPF duck closed flock according to claim 2, characterized in that, The process of obtaining historical concentration sequences to determine concentration limits includes: During the empty house purification period before the basic duck flock to be purified is moved into the isolation environment, the bioaerosol detection device deployed in the isolation environment is driven to perform air sampling, continuously obtain the pathogen concentration quantification value in the empty house state, and define the pathogen concentration quantification value obtained within the preset time window as the historical concentration sequence. Calculate the arithmetic mean of the historical concentration series to obtain the noise floor mean, and calculate the standard deviation of the historical concentration series to obtain the noise dispersion; Multiply the noise dispersion by the preset threshold coefficient to obtain the fluctuation tolerance term, and add the fluctuation tolerance term to the mean noise floor value to obtain the concentration limit.

4. The method for purifying and cultivating an SPF duck closed flock according to claim 2, characterized in that, The arithmetic mean of each signal in the multidimensional physicochemical signal within the preset sliding time window is calculated to obtain the corresponding smoothed mean, and a physical sequence with timestamps is generated based on the smoothed mean. A discrete sequence is generated based on the biochemical detection signal and the concentration limit, including: Within a continuous sliding time window, the arithmetic mean of ambient temperature, relative humidity, air pressure, ammonia concentration, and air velocity are calculated respectively to obtain the smoothed mean of each multidimensional physicochemical signal. The smoothed mean values ​​of various multidimensional physicochemical signals within the same sampling period are concatenated into multidimensional vectors and arranged in chronological order of sampling time to generate a physical sequence. When the biochemical detection signal is less than the concentration limit, a first state value representing the microbiological safety status is generated; when the biochemical detection signal is greater than or equal to the concentration limit, a second state value representing the microbiological risk status is generated. Arrange the continuously generated first or second state values ​​in chronological order of their entry time to generate a discrete sequence.

5. The method for purifying and cultivating an SPF duck closed flock according to claim 2, characterized in that, The step of obtaining the time series deviation based on the relative time difference between the physical sequence and the discrete sequence, and the air pressure deviation between the current transient sample value and its smoothed mean value of the air pressure signal in the physical sequence, includes: The relative time difference is obtained by subtracting the current timestamp corresponding to the current acquisition point in the physical sequence from the entry timestamp corresponding to the latest data entry in the discrete sequence and taking the absolute value. The pressure deviation is obtained by subtracting the current transient sample value of the air pressure signal in the physical sequence from its smoothed mean and taking the absolute value. The relative time difference is dimensionless by using preset time scale parameters to obtain a dimensionless relative time difference, and the pressure deviation is dimensionless by using preset pressure scale parameters to obtain a dimensionless pressure deviation. Using the natural constant as the base and the dimensionless pressure deviation as the exponent, an exponential mapping calculation is performed to obtain the pressure deviation mapping term; Multiplying the pressure deviation mapping term by the dimensionless relative time difference yields the time series deviation.

6. The method for purifying and cultivating an SPF duck closed flock according to claim 5, characterized in that, The step of performing time-axis matching between the physical sequence and the discrete sequence based on the time-series deviation to construct a multi-source data matrix includes: Determine the matching time window between the physical sequence and the discrete sequence based on the time series deviation; Within the matching time window, the state values ​​in the discrete sequence are aligned with the physical sequence data whose time position is closest. The data gaps in the physical sequence are filled using an effective value decay preservation method, and the data gaps in the discrete sequence are filled using a state preservation method. The aligned and filled physical sequences are spliced ​​together with the discrete sequences to construct a multi-source data matrix.

7. The method for purifying and cultivating an SPF duck closed flock according to claim 2, characterized in that, The maximum environmental offset rate is obtained based on the smoothed average values ​​of ambient temperature, relative humidity, air pressure, and ammonia concentration within the multi-source data array. The comprehensive deviation is then calculated by combining the target heterozygosity with the actual heterozygosity, including: The preset safe ranges corresponding to ambient temperature, relative humidity, air pressure, and ammonia concentration are obtained respectively. Based on the degree of deviation of the smoothed mean values ​​of ambient temperature, relative humidity, air pressure, and ammonia concentration from the corresponding preset safe range, the ambient temperature deviation rate, relative humidity deviation rate, air pressure deviation rate, and ammonia concentration deviation rate are calculated respectively. The maximum value among the ambient temperature offset rate, relative humidity offset rate, air pressure offset rate, and ammonia concentration offset rate is extracted as the maximum environmental offset rate feature. The ratio of the target heterozygosity to the actual heterozygosity is used to obtain the heterozygosity gain ratio. The overall deviation is obtained by multiplying the hybridity gain ratio by the maximum environmental offset characteristic.

8. The method for purifying and cultivating an SPF duck closed flock according to claim 2, characterized in that, The wind speed parameter is obtained by calculating the smoothed mean value of airflow velocity in the physical sequence, and the smoothing suppression amount is obtained by combining the comprehensive deviation, including: The wind speed parameter is dimensionless by using a preset air velocity scale parameter to obtain a dimensionless wind speed parameter. Using the natural constant as the base and the negative of the dimensionless wind speed parameter as the exponent, an exponential mapping calculation is performed to obtain the wind speed attenuation mapping term; The difference between constant 1 and the wind speed attenuation mapping term is calculated to obtain the airflow cutoff gain term; Multiplying the overall deviation by the airflow cutoff gain term yields the smoothing suppression amount.

9. The method for purifying and cultivating an SPF duck closed flock according to claim 8, characterized in that, The method of adjusting the operating frequency of the variable frequency motor and the step opening of the ventilation regulating valve includes: Based on the abnormal physicochemical signals that cause an increase in overall deviation, determine the adjustment direction of the variable frequency motor and the ventilation regulating valve; The smoothing suppression amount is used as the amplitude correction coefficient to limit the single-run frequency change of the variable frequency motor and the single-step opening of the ventilation regulating valve. When the ammonia concentration is higher than the corresponding preset safety range, the exhaust volume after high-efficiency filtration or disinfection is increased according to the amplitude correction factor. When the air velocity is higher than the corresponding preset safety range, the frequency increase of the variable frequency motor is limited or the single opening change of the ventilation regulating valve is reduced according to the amplitude correction coefficient. When the positive pressure gradient of the isolated environment is lower than the preset positive pressure gradient, the air supply volume is increased or the exhaust volume is decreased according to the amplitude correction coefficient.

10. The method for purifying and cultivating an SPF duck closed flock according to claim 2, characterized in that, The step of obtaining a correction value based on the hatching deviation rate extracted from the actual hatching rate and the designed hatching rate, and updating the safety limit threshold for the next breeding cycle based on the correction value, includes: The hatching deviation rate is obtained by subtracting the designed hatching rate from the actual hatching rate. Using the natural constant as the base and the hatching deviation rate as the exponent, an exponential mapping calculation is performed to obtain the deviation mapping term; Calculate the difference between the deviation mapping term and the constant 1 to obtain the error residual term; Multiply the preset learning step size term by the error residual term to obtain the correction value; The difference between constant 1 and the correction value is calculated to obtain the attenuation gain term; Multiply the safety threshold that is effective in the current breeding cycle by the decay gain term to obtain the candidate updated safety threshold. The candidate update safety threshold is restricted to between the preset lower threshold and the preset upper threshold to obtain the update safety threshold. The decision threshold parameter in the control terminal is overwritten with the updated safety limit threshold, which will then serve as the safety limit threshold for the next breeding cycle.

Images