Breathing estimation method and system in water-free environment in fish keep-alive process
By monitoring and calculating the oxygen consumption rate and carbon dioxide release rate of fish in anhydrous environment, the problem of lack of research on fish respiration rates in anhydrous environments is solved, and the survival rate of fish during anhydrous transportation is improved.
Patent Information
- Application Number
- CN202411407248.X
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-10
- Publication Date
- 2025-01-10
AI Technical Summary
In the prior art, there are few researches on the respiration rate of fish in anhydrous environments, and there is a lack of effective monitoring and estimation methods, making it difficult to judge the vitality characteristics of fish during invigorating transportation.
Provide a respiration estimation method in anhydrous environment of fish in the process of maintaining fish. By obtaining the weight, volume, volume of sealed packaging bags and initial oxygen volume of fish individuals, collecting timing parameters in anhydrous environment, and calculating correlation coefficients and rate of change Coefficients, calculate the oxygen consumption rate and carbon dioxide release rate, and then estimate the respiration rate per unit time of fish.
This method helps to optimize the respiration correlation parameters of anhydrous environment, improve the resilience rate of fish during anhydrous transportation, and provides an effective method for monitoring and estimating respiration rate in anhydrous environment of fish.
Smart Images

Figure CN119293446A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the technical field of logistics and transportation monitoring, and in particular to a method and system for estimating respiration in a waterless environment during a fish keeping alive process. Background Art
[0002] Respiration is an important indicator of fish's physiological metabolism and vitality level. Fish obtain oxygen through breathing to maintain life activities. The respiration rate can reflect the oxygen consumption of fish's life activities. At the same time, the stress level and vitality state of live fish can be judged through respiration. Therefore, monitoring and estimating the respiration rate of fish can help better judge the vitality characteristics of fish during the process of keeping them alive and transporting them.
[0003] There are many factors that affect the respiratory state of aquatic organisms. It is affected by internal factors such as aquatic product category and genes, as well as external environmental stress factors. Environmental temperature, gas, stress level, etc. all have important effects on the respiratory state of fish. Current research focuses more on the estimation of the respiration of aquatic organisms in water environments, and there are fewer studies on the estimation of respiration in waterless environments. Therefore, it is necessary to propose a new method for estimating the respiratory rate of individual fish in waterless environments. Summary of the invention
[0004] In order to overcome the deficiencies of the prior art, the purpose of the present invention is to provide a method and system for estimating respiration in a waterless environment during the fish keeping alive process, so as to solve the problem of lack of research on the respiration rate of fish in a waterless environment in the prior art.
[0005] To achieve the above object, the present invention provides the following solutions:
[0006] A method for estimating respiration in a waterless environment during fish preservation, comprising:
[0007] Obtain the weight of the individual fish to be tested, the volume of the fish body, the volume of the sealed packaging bag, and the initial amount of oxygen in the sealed packaging bag, and collect the timing parameters associated with breathing in a waterless environment with a preset sampling period and a preset sampling step length during the survival process of the individual fish in the preset temperature zone:
[0008] Divide the preset sampling period into multiple time periods, and calculate the correlation coefficient and change rate coefficient of the timing parameters of each time period under the preset temperature zone;
[0009] Calculate the oxygen consumption rate and carbon dioxide release rate in the anhydrous environment during the survival process of individual fish in each period of the preset temperature zone;
[0010] The respiration rate of fish per unit time is calculated based on the oxygen consumption rate and carbon dioxide release rate in each period under the preset temperature zone and preset cycle.
[0011] Preferably, the timing parameters include: temperature, oxygen, carbon dioxide and gill amplitude.
[0012] Preferably, the preset sampling period is divided into multiple time periods, including:
[0013] Obtain standard reference values of parameters of temperature, oxygen, carbon dioxide and gill amplitude within a preset sampling period; the calculation formula of the standard reference value is:
[0014]
[0015] Where X represents the standard reference value of the parameter within the preset sampling period, X i (max) indicates the maximum value of the parameter within the preset sampling period, X i (min) represents the minimum value of the parameter within the preset sampling period, n represents the number of timing parameters, and i represents the timing parameter sequence;
[0016] Obtain the fluctuation deviation of the parameters of temperature, oxygen, carbon dioxide and gill amplitude within a preset sampling period; the calculation formula of the fluctuation deviation is:
[0017]
[0018] Among them, γ is the fluctuation deviation value of each parameter within the preset sampling period, X i is the sampling value within the preset sampling period;
[0019] The boundary of the fluctuation deviation is used as a time period division point to perform multiple time period divisions.
[0020] Preferably, the step of calculating the correlation coefficient comprises:
[0021] Calculate the correlation coefficients between the temperature, oxygen, carbon dioxide, and gill amplitude parameters of each time period under the preset temperature zone, and compare the orientation of the correlation coefficients between the parameters of each time period. If the correlation coefficients between the parameters of each time period are in opposite orientations, check the data consistency and identify abnormalities.
[0022] Preferably, the calculation formula of the change rate coefficient is:
[0023]
[0024] in, A represents the average value of the parameter in the Δt period, i (Δt)max represents the maximum value of the parameter within the Δt period. is the average value within the period Δt+1.
[0025] Preferably, the calculation formulas for the oxygen consumption rate and the carbon dioxide release rate are as follows:
[0026] OCR(O2)=[k·(V P -V W )·W N ] / (G(B a -B i ))×100%
[0027] OCR(CO2)=[k·(V P -V W )·W N ] / (G(C a -C i ))×100%
[0028] Wherein, OCR(O2) is the oxygen consumption rate, OCR(CO2) is the carbon dioxide release rate, k is a parameter related to the individual difference characteristics of fish, k∈(0,1], V P V is the volume of the sealed packaging bag during the process of keeping fish alive without water. W is the fish body volume, W N is the weight of the fish at each time period, G(B a ), G(B i ) are the starting concentration and ending concentration of oxygen consumption of fish in each period, G(C a ), G(C i ) are the starting concentration and ending concentration of carbon dioxide released by fish in each time period.
[0029] Preferably, the calculation formula of the respiratory rate is:
[0030] BR N =(θ·(OCR(O2))+(1-θ)·(OCR(CO2)))·ε ci
[0031] Wherein, θ represents the first weight coefficient, which is comprehensively determined according to the change rate coefficient of the timing parameter. When OCR(O2)>OCR(CO2), θ>(1-θ); when OCR(O2)≤OCR(CO2), θ<(1-θ), ε ci is the correction coefficient of the change rate coefficient of the gill amplitude, if δ ci is not 0, then ε ci The value is 1, if δ ci is 0, then ε ci The value is 0, BR N is the respiration rate.
[0032] Preferably, after calculating the respiration rate of the fish per unit time according to the oxygen consumption rate and carbon dioxide release rate in each time period in the preset temperature zone and preset cycle, the method further comprises:
[0033] The estimated value of the comprehensive respiratory rate of the individual fish in the preset temperature zone and preset period is determined according to the respiratory rate; the calculation formula of the estimated value is:
[0034]
[0035] Among them, w N is the second weight coefficient, and the second weight coefficient is based on BR N To determine, if BR N ≥BR N-1 , then w N ≥w N-1 , w1+w2+……w N =1, and is the estimated value, BR min and BR max are the minimum and maximum respiration rates, respectively.
[0036] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:
[0037] The present invention provides a method and system for estimating respiration in an anhydrous environment during the fish preservation process. The method includes: obtaining the weight of the individual fish to be tested, the volume of the fish body, the volume of the sealed packaging bag, and the initial oxygen content in the sealed packaging bag, and collecting the timing parameters associated with the respiration in an anhydrous environment with a preset sampling period and a preset sampling step length during the fish individual preservation process in a preset temperature zone: dividing the preset sampling period into multiple time periods, and calculating the correlation coefficient and the change rate coefficient of the timing parameters of each time period in the preset temperature zone; calculating the oxygen consumption rate and carbon dioxide release rate in an anhydrous environment during the fish individual preservation process in each time period in the preset temperature zone; calculating the fish respiration rate per unit time according to the oxygen consumption rate and carbon dioxide release rate in each time period in the preset temperature zone and the preset period. The present invention is helpful to optimize the parameters associated with respiration in an anhydrous environment and improve the survival rate of fish during anhydrous transportation. BRIEF DESCRIPTION OF THE DRAWINGS
[0038] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.
[0039] Figure 1 A flow chart of a method provided by an embodiment of the present invention;
[0040] Figure 2 A schematic diagram of steps provided for an embodiment of the present invention;
[0041] Figure 3 A schematic diagram of the system structure provided by an embodiment of the present invention;
[0042] Figure 4 A schematic diagram of the structure of a device provided in an embodiment of the present invention. DETAILED DESCRIPTION
[0043] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.
[0044] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.
[0045] Figure 1 A flow chart of a method provided by an embodiment of the present invention, such as Figure 1 As shown, the present invention provides a method for estimating respiration in a waterless environment during a fish keeping alive process, comprising:
[0046] Step 100: Obtain the weight of the individual fish to be tested, the volume of the fish body, the volume of the sealed packaging bag, and the initial amount of oxygen in the sealed packaging bag, and collect the timing parameters associated with breathing in a waterless environment with a preset sampling period and a preset sampling step length during the survival process of the individual fish in the preset temperature zone:
[0047] Step 200: Divide the preset sampling period into multiple time periods, and calculate the correlation coefficient and change rate coefficient of the timing parameters of each time period in the preset temperature zone;
[0048] Step 300: Calculate the oxygen consumption rate and carbon dioxide release rate in the anhydrous environment during the fish individual survival process in each period in the preset temperature zone;
[0049] Step 400: Calculate the respiration rate of the fish per unit time according to the oxygen consumption rate and carbon dioxide release rate in each time period in the preset temperature zone and preset cycle.
[0050] Figure 2 is a flow chart of a multi-stage evaluation method for the storage and transportation of raw aquatic products provided by an embodiment of the present invention, such as Figure 2 As shown, including:
[0051] S1, obtain the weight W and volume V of the fish to be tested W , Volume of sealed packaging bag V PThe initial oxygen volume B0 in the sealed packaging bag and the breathing-related timing parameters of the fish individuals collected in the preset temperature zone Tem in the waterless environment with the same preset sampling period ΔT and the same preset sampling step length t: temperature a i 、Oxygen b i 、CO2 i , gill amplitude d i Timing parameters of
[0052] Specifically, the fish individuals to be tested are sea bass and salmon, and the preset temperature zones are set to 5 temperature zones: the first temperature zone [0-1.5℃], the second temperature zone [1-2.5℃], the third temperature zone [1.5-3℃], the fourth temperature zone [2-3.5℃], and the fifth temperature zone [2.5-4℃]. The respiratory-related timing parameters with the same preset sampling period ΔT and the same preset sampling step t (t<ΔT) in the 5 temperature zones are obtained respectively: the timing parameters of temperature, oxygen, carbon dioxide, and fish gill amplitude.
[0053] S2, dividing the preset sampling period ΔT into N time periods;
[0054] Specifically, the standard reference values of temperature, oxygen, carbon dioxide, and gill amplitude parameters within the preset sampling period ΔT are obtained and processed by the following formula:
[0055]
[0056] Where X represents the standard reference value of the parameter within the preset sampling period ΔT, X i (max) represents the maximum value of the parameter within the preset sampling period ΔT, X i (min) represents the minimum value of the parameter within the preset sampling period ΔT, and n represents the number of data;
[0057] Obtain the fluctuation deviation of temperature, oxygen, carbon dioxide, and gill amplitude parameters within the preset sampling period ΔT and process them using the following formula:
[0058]
[0059] Among them, γ∈[0,1] is the fluctuation deviation value of each parameter within the preset sampling period ΔT, X i is the sampling value within the preset sampling period ΔT, and X is the standard reference value within the preset sampling period ΔT.
[0060] The deviation limit is used as the time period division point to divide the time period into N time periods.
[0061] In the actual process of fish transportation without water, γ>0.5 is specifically set as the deviation limit for time period division.
[0062] S3, calculate the correlation coefficient r of temperature, oxygen, carbon dioxide, and gill amplitude parameters in each period under the preset temperature zone i , specifically including:
[0063] The correlation analysis between the time series parameters of each divided period is carried out. i ∈[-1,1], and compare the correlation coefficients of each time period based on historical experience knowledge. For example, oxygen and carbon dioxide are negatively correlated, and oxygen and respiratory amplitude are positively correlated. If the correlation coefficients between the parameters of a certain time period are opposite, the data needs to be verified and abnormal judgment processed, and the correlation judgment is performed between the processing until they are consistent.
[0064] The consistency of the correlation coefficients in each period can indicate the consistency of the data and provide reliability for the estimation of respiratory rate.
[0065] S4, calculating the change rate coefficient of two adjacent divided time periods in the preset temperature zone is specifically:
[0066]
[0067] in, A represents the average value of the parameter in the Δt period, i (Δt)max represents the maximum value of the parameter within the Δt period. is the average value within the period Δt+1.
[0068] It should be pointed out that δ i ∈[0,1), during the waterless live transport process, the rate of change of the parameters at each time period reflects the respiratory state of the fish, δ i The greater the rate of change, the greater the impact on the respiratory state of fish.
[0069] It is necessary to further explain that the oxygen consumption of fish is greatly affected by temperature, oxygen and carbon dioxide parameters. The delta values of each parameter are set according to the preset knowledge base. i (First temperature zone)≤δ i (Second temperature zone)≤δ i (Third temperature zone)≤δ i (Fourth temperature zone)≤δ i (Fifth temperature zone).
[0070] S4: Calculate the oxygen consumption rate OCR (O2) and carbon dioxide release rate OCR (CO2) in the anhydrous environment during the fish individual survival process in each period of the preset temperature zone and preset period ΔT. The specific calculation is:
[0071] OCR(O2)=[k·(V P -V W )·W N ] / (G(Ba -B i ))×100%
[0072] OCR(CO2)=[k·(V P -V W )·W N ] / (G(C a -C i ))×100%
[0073] Among them, k is related to the individual difference characteristics of fish, K∈(0,1],V P V is the volume of the sealed packaging bag during the process of keeping fish alive without water. W is the fish body volume, W N is the fish body weight at each time period, G(B a ), G(B i ) are the initial and final oxygen concentrations of fish in each period, assuming that in the first period G(B a )=B0,G(C a ), G(C i ) are the starting and ending concentrations of carbon dioxide for fish in each time period.
[0074] It should be noted that it is assumed that the changes in fish weight and volume during the waterless transportation of fish are very small, that is, there is no change treatment.
[0075] In addition, the respiration rates of different fish species are not only affected by the environment, but also related to their own individual differences. For example, the suitable temperature zone for sea bass is different from that for salmon. The oxygen consumption of salmon obtained based on the knowledge base has a greater impact, so the k value of salmon is larger than that of sea bass.
[0076] S5: Calculate the oxygen consumption rate OCR (O2) and carbon dioxide release rate OCR (CO2) in each period of the preset temperature zone and preset period ΔT, and calculate the respiratory rate of the fish individual in each period by the following formula:
[0077] BR=(θ·(OCR(O2))+(1-θ)·(OCR(CO2)))·ε ci
[0078] Among them, θ represents the weight coefficient, which is comprehensively determined based on the change rate coefficient of temperature, oxygen, carbon dioxide and gill amplitude. When OCR(O2)>OCR(CO2), θ>(1-θ), otherwise θ<(1-θ).
[0079] ε ci is the coefficient of change of gill amplitude δ ci The relevant correction coefficient of ci is not 0, then εci The value is 1, if δ ci is 0, then ε ci The value is 0.
[0080] The embodiment of the present invention estimates the respiratory rate of individual fish in a waterless environment, which helps to optimize the respiratory-related parameters of the waterless environment and improve the survival rate of fish during waterless transportation.
[0081] The following is a description of a system for estimating respiration in an anhydrous environment for a fish keeping alive process provided by an embodiment of the present invention. The system for estimating respiration in an anhydrous environment for a fish keeping alive process described below and the method for estimating respiration in an anhydrous environment for a fish keeping alive process described above can be referred to each other. Figure 3 As shown, the system includes:
[0082] The acquisition unit acquires the weight of the individual fish to be tested, the volume of the fish body, the volume of the sealed packaging bag, and the initial amount of oxygen in the sealed packaging bag, and collects the timing parameters associated with the breathing in the waterless environment with the preset sampling period and the preset sampling step length during the survival process of the individual fish in the preset temperature zone:
[0083] A first calculation unit, used to divide the preset sampling period into a plurality of time periods, and calculate the correlation coefficient and the change rate coefficient of the timing parameters of each time period in the preset temperature zone;
[0084] The second calculation unit is used to calculate the oxygen consumption rate and carbon dioxide release rate of the anhydrous environment during the fish individual survival process in each period of the preset temperature zone;
[0085] The third calculation unit is used to calculate the respiration rate of the fish per unit time according to the oxygen consumption rate and carbon dioxide release rate in each time period in the preset temperature zone and preset cycle.
[0086] In addition, this embodiment also provides a schematic diagram of the structure of a multi-stage evaluation device for the storage and transportation process of raw aquatic products, such as Figure 4 As shown, it includes: an acquisition module 41, a monitoring module 42, a first processing module 43 and a second processing module 44; wherein:
[0087] The acquisition module 41 is used to obtain basic data of the individual fish to be tested, including fish weight, fish volume, sealed packaging bag volume and initial oxygen content in the sealed packaging bag;
[0088] Monitoring module 42, used for real-time collection of temperature, oxygen, carbon dioxide, gill amplitude and other data in a water-free environment under a preset period, and for time division;
[0089] The first processing module 43 is used to calculate the correlation coefficient and the change rate coefficient of the temperature, oxygen, carbon dioxide, and gill amplitude parameters in each time period in the preset temperature zone;
[0090] The second processing module 44 is used to calculate the oxygen consumption rate and carbon dioxide release rate in the anhydrous environment during the fish individual survival process in each period of the preset temperature zone, and establish a fish individual breathing rate estimation model in the anhydrous environment.
[0091] The embodiment of the present invention estimates the respiratory rate in different time periods by considering multiple key influencing parameters of respiration, which helps to optimize the respiration-related parameters in an anhydrous environment and improve the survival rate of fish during anhydrous transportation.
[0092] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.
[0093] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of the present invention. At the same time, for those skilled in the art, according to the idea of the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.
Claims
1. A method for estimating respiration in a waterless environment during fish preservation, characterized in that: include: Obtain the weight of the individual fish to be tested, the volume of the fish body, the volume of the sealed packaging bag, and the initial amount of oxygen in the sealed packaging bag, and collect the timing parameters associated with breathing in a waterless environment with a preset sampling period and a preset sampling step length during the survival process of the individual fish in the preset temperature zone: Divide the preset sampling period into multiple time periods, and calculate the correlation coefficient and change rate coefficient of the timing parameters of each time period under the preset temperature zone; Calculate the oxygen consumption rate and carbon dioxide release rate in the anhydrous environment during the survival process of individual fish in each period of the preset temperature zone; The respiration rate of fish per unit time is calculated based on the oxygen consumption rate and carbon dioxide release rate in each period under the preset temperature zone and preset cycle.
2. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 1, characterized in that: The timing parameters include: temperature, oxygen, carbon dioxide and fish gill amplitude.
3. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 2, characterized in that: The preset sampling period is divided into multiple time periods, including: Obtain standard reference values of parameters of temperature, oxygen, carbon dioxide and gill amplitude within a preset sampling period; the calculation formula of the standard reference value is: Where X represents the standard reference value of the parameter within the preset sampling period, X i (max) indicates the maximum value of the parameter within the preset sampling period, X i (min) represents the minimum value of the parameter within the preset sampling period, n represents the number of timing parameters, and i represents the timing parameter sequence; Obtain the fluctuation deviation of the parameters of temperature, oxygen, carbon dioxide and gill amplitude within a preset sampling period; the calculation formula of the fluctuation deviation is: Among them, γ is the fluctuation deviation value of each parameter within the preset sampling period, X i is the sampling value within the preset sampling period; The boundary of the fluctuation deviation is used as a time period division point to perform multiple time period divisions.
4. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 2, characterized in that: The calculation steps of the correlation coefficient include: Calculate the correlation coefficients between the temperature, oxygen, carbon dioxide, and gill amplitude parameters of each time period under the preset temperature zone, and compare the orientation of the correlation coefficients between the parameters of each time period. If the correlation coefficients between the parameters of each time period are in opposite orientations, check the data consistency and identify abnormalities.
5. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 2, characterized in that: The calculation formula of the change rate coefficient is: Among them, A i (Δt) represents the average value of the parameter within the Δt period, A i (Δt)max represents the maximum value of the parameter within the Δt period. is the average value within the period Δt+1.
6. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 2, characterized in that: The calculation formulas for the oxygen consumption rate and the carbon dioxide release rate are as follows: OCR(O2)=[k·(V P -V W )·W N ] / (G(B a -B i ))×100% OCR(CO2)=[k·(V P -V W )·W N ] / (G(C a -C i ))×100% Wherein, OCR(O2) is the oxygen consumption rate, OCR(CO2) is the carbon dioxide release rate, k is a parameter related to the individual difference characteristics of fish, k∈(0,1], V P V is the volume of the sealed packaging bag during the process of keeping fish alive without water. W is the fish body volume, W N is the weight of the fish at each time period, G(B a ), G(B i ) are the starting concentration and ending concentration of oxygen consumption of fish in each period, G(C a ), G(C i ) are the starting concentration and ending concentration of carbon dioxide released by fish in each time period.
7. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 6, characterized in that: The calculation formula of the respiratory rate is: BR N =(θ·(OCR(O2))+(1-θ)·(OCR(CO2)))·ε ci Wherein, θ represents the first weight coefficient, which is comprehensively determined according to the change rate coefficient of the timing parameter. When OCR(O2)>OCR(CO2), θ>(1-θ); when OCR(O2)≤OCR(CO2), θ<(1-θ), ε ci is the correction coefficient of the change rate coefficient of the gill amplitude, if δ ci is not 0, then ε ci The value is 1, if δ ci is 0, then ε ci The value is 0, BR N is the respiration rate.
8. The method for estimating the respiratory rate of fish in an anhydrous environment during the process of keeping fish alive according to claim 7, characterized in that: After calculating the respiration rate of fish per unit time according to the oxygen consumption rate and carbon dioxide release rate in each period of the preset temperature zone and preset cycle, it also includes: The estimated value of the comprehensive respiratory rate of the individual fish in the preset temperature zone and preset period is determined according to the respiratory rate; the calculation formula of the estimated value is: Among them, w N is the second weight coefficient, and the second weight coefficient is based on BR N To determine, if BR N ≥BR N-1 , then w N ≥w N-1 , w1+w2+……w N =1, and is the estimated value, BR min and BR max are the minimum and maximum respiration rates, respectively.
9. A breathing estimation system in an anhydrous environment during fish preservation, characterized in that: include: The acquisition unit acquires the weight of the fish to be tested, the volume of the fish body, the volume of the sealed packaging bag, and the initial oxygen content in the sealed packaging bag, and collects the timing parameters associated with the breathing in the waterless environment with the preset sampling period and the preset sampling step length during the survival process of the fish in the preset temperature zone: A first calculation unit, used to divide the preset sampling period into a plurality of time periods, and calculate the correlation coefficient and the change rate coefficient of the timing parameters of each time period in the preset temperature zone; The second calculation unit is used to calculate the oxygen consumption rate and carbon dioxide release rate of the anhydrous environment during the fish individual survival process in each period of the preset temperature zone; The third calculation unit is used to calculate the respiration rate of the fish per unit time according to the oxygen consumption rate and carbon dioxide release rate in each time period in the preset temperature zone and preset cycle.