A method for assessing fish resources in shallow water fisheries

Abstract

The invention discloses a method for assessing the fish resource quantity in a shallow-water fishery water body, and relates to the field of fishery acoustic technology. The method specifically comprises the steps of designing a sampling route, collecting acoustic data, determining density conditions, analyzing approach-avoidance effects, determining an acoustic target intensity-body length conversion equation, and estimating specifications and biomass. The method of the invention is applicable to various types of shallow-water fishery waters with a water depth of 2 to 5 m, including lakes, reservoirs, ponds, ditches, etc., and can assess all fish distributed above the bottom of the water body. The method is particularly applicable to rapid assessment of fish resource quantity in shallow-water aquaculture and aquaculture fishery waters, and has the outstanding advantages of being non-destructive, precise, and efficient.

Description

Technical Field

[0001] The invention relates to the technical field of fishery acoustics, and in particular to a method for evaluating the amount of fish resources in a shallow water fishery body. Background Art

[0002] Among common media such as sound, light, and electromagnetism, ultrasound is the only effective medium that can propagate over long distances in water. Therefore, the use of hydroacoustic methods to observe and study fish and other living organisms or water environment characteristics in water bodies is a widely used technical means. In existing applications, the cruise detection method using a hydroacoustic transducer pointing vertically to the bottom of the water is widely used in the assessment of fish resources in the ocean and deep lakes and reservoirs, and can achieve good assessment results. However, in shallow (water depth between 2 and 5 meters) lakes, reservoirs, ponds, and ditches, the vertical cruise detection method has a very small sampled water volume due to the influence of the near-field effect, and its applicability is relatively poor.

[0003] However, in shallow lakes, reservoirs, ponds, ditches and other aquaculture waters, the evaluation of fish behavior, density, size and resource quantity is of great significance for the scientific management of the aquaculture process, precise feeding, disease prevention and control, and benefit prediction, which can significantly improve production efficiency and benefits. At present, the traditional evaluation methods for fish resources in such shallow waters mainly include biological sampling and optical observation, which have problems such as damage, high cost, limited observation range and low efficiency. The hydroacoustic method has outstanding advantages such as non-destructiveness, accuracy and efficiency. Among them, the horizontal hydroacoustic method has a good application effect in shallow water systems, but its technical requirements are quite different from those of vertical hydroacoustics. There is an urgent need for an accurate evaluation method for fish resources in shallow waters. Summary of the invention

[0004] The purpose of the present invention is to provide a method for assessing the amount of fish resources in shallow water fishery water bodies, so as to solve the problems existing in the above-mentioned prior art and enable the amount of fish resources in shallow water aquaculture fishery water bodies to be accurately assessed.

[0005] To achieve the above object, the present invention provides the following solutions:

[0006] The present invention provides a method for assessing the amount of fish resources in a shallow water fishery body, which specifically comprises the following steps:

[0007] S1, design the sampling route, determine the trajectory, number of sections and spacing of the assessment route according to the coefficient of variation of the resource assessment results that can be accepted within the estimated assessment requirements;

[0008] S2, collect acoustic data, use the horizontal acoustic detection method, and conduct navigation acoustic survey sampling according to the designed sampling route;

[0009] S3, determining the density status, determining whether the fish density distribution status is in a discrete distribution state, when it is determined that the target is in a discrete distribution state, the fish target intensity observed in situ is used for fish body length conversion;

[0010] S4, analyze the approach-avoidance effect, use the theoretical model to correct the sampling bias caused by the fish approach-avoidance effect, and estimate the actual average density of fish;

[0011] S5, determine the acoustic target intensity-body length conversion equation, obtain the equation using field measurements, or directly use the existing empirical formula;

[0012] S6, estimate the size and biomass, use the empirical formula of horizontal acoustic target intensity-body length and the body length-weight relationship to calculate the average size of fish, and use the cross-section method to estimate the total biomass in the assessed waters.

[0013] Preferably, in S1, the sampling route adopts a square waveform or a triangular waveform route, and the spacing between adjacent parallel routes depends on the requirements of the assessment time and assessment accuracy; to ensure the assessment accuracy, the coefficient of variation of the resource assessment results acceptable within the assessment requirements is estimated according to the following formula, so as to determine the number and spacing of the sections of the assessment route:

[0014]

[0015] Where, CV is the coefficient of variation of the resource assessment results that the assessment requires to be acceptable, that is, the ratio of the standard deviation to the mean value; Λ is the monitoring coverage;

[0016]

[0017] Where D is the length of the evaluation route, in meters; A is the area of ​​the evaluation area, in square meters. According to the requirements, CV is less than 0.25 and Λ is greater than 4, and then the length of the evaluation route is determined. Then, the number of sections and spacing of the evaluation route are calculated from the area of ​​the evaluation area.

[0018] Preferably, in S2, two transducers are symmetrically arranged on both sides of the ship's side, the transducers are immersed in water to a depth of at least 0.5m, and the horizontal directions of the transducers are opposite and perpendicular to the direction of the ship's travel; before the evaluation begins, the scientific fish finder system is acoustically calibrated using the standard sphere method with reference to the instrument operation manual, and the calibration of the split-beam scientific fish finder system includes the system transceiver signal gain and beam parameter calibration.

[0019] Preferably, in S2, during the acoustic data underway sampling process, the set sampling parameters need to remain consistent, the speed is lower than 5km / h, the pulse length is set according to the distribution characteristics and water depth of the assessment object, and the data collection and recording range is set according to the water depth of the assessment water area; the pulse length selects a short pulse less than 0.2ms, and the data collection range is greater than 10m.

[0020] Preferably, in S3, when the fish cluster distribution in the assessment water area leads to too high density, an acoustic shading effect will be generated, resulting in a nonlinear relationship between the volume scattering intensity and the fish density, which in turn leads to analysis errors. The Nv value is calculated using the following formula (3) to determine whether the target is in a discrete distribution state:

[0021]

[0022] Where c is the speed of sound in meters per second; τ is the pulse length in seconds; Ψ is the decibel representation of the effective beam solid angle in decibels; R is the distance between the target and the transducer surface in meters; ρ is v —Volume fish density, ind. / m 3 , calculated by the following formula (4);

[0023]

[0024] In the formula, s v —Volume scattering coefficient, in m 2 / m 3 ; σ bs —Backscattering area, in square meters, S v —Volume scattering intensity, in decibels, derived from the acoustic data post-processing software in S2; TS—Fish target intensity, in decibels, derived from the acoustic data post-processing software in S2;

[0025] When Nv<0.1, the target is in a discrete distribution state, and the in situ observed fish target intensity is used to convert the fish body length; when Nv>0.1, the target is in a non-discrete distribution state, and the relevant analysis units or areas need to be excluded from the analysis.

[0026] Preferably, during the acoustic sampling survey, due to the avoidance behavior of fish towards detection activities, the spatial distribution density of fish will deviate from the normal state. It is necessary to use a theoretical model to correct the deviation according to the fish density distribution to obtain the correct fish density value; first, the water layer is divided symmetrically on both sides with the acoustic sampling track as the center, starting from 2m away from the track, and several water layers are divided with a spacing of 1m. It is analyzed whether the fish density in each water layer is uniform or there is a density gradient. When the fish distribution in each water layer is uniform, it is confirmed that the fish avoidance effect is not obvious, and the fish biomass in the water area is directly estimated using the density in S3.

[0027] Preferably, when there is a density gradient in the distribution of fish in each water layer, the following theoretical model is used to correct the deviation caused by the sampling process and estimate the actual average density of fish:

[0028]

[0029] Where, d is the distance from the water layer to the transducer, in meters; δa d —The sound density obtained in the water layer corresponding to this distance, ind. / km 2 ; l—the total number of water layers divided, which is a positive integer; δ max —The maximum density recorded in each water layer, ind. / km 2 ;

[0030] Therefore, the actual average density of fish in the water body δe is:

[0031] δe=δa / P (6)

[0032] Preferably, the field measurement method for determining the acoustic target intensity-fish body length conversion equation in S5 is as follows:

[0033] First, a net cage with an opening on the top is made. The net is made of polyethylene monofilament and the mesh size is 1cm-2cm. A float is installed on the top of the net cage and a sinker is installed on the bottom, so that the net cage can float on the water surface and has a regular square shape after being stretched. According to the biological characteristics of the fish to be studied, the appropriate sample specifications and quantity are determined, and samples are taken at equal intervals between 10cm and the maximum recorded body length of the fish. The body length interval is between 2cm-20cm, and the number of samples is between 10 and 30. The acoustic target intensity of fish samples of different specifications is measured horizontally from the side of the net cage, and the number of clear and identifiable single target echo signals obtained from the measurement of each sample should be more than 300. A regression equation between the horizontal acoustic target intensity of fish and the body length of the fish is established, and the standard form of the equation is: TS=20.0×log(BL)-b20, where b20 is a constant.

[0034] The existing empirical formulas for determining the acoustic target intensity-body length conversion equation include: the conversion formula for the horizontal acoustic target intensity-fish body length of the common freshwater aquaculture fish carp is:

[0035] TS = 24.63 × log (SL) - 93.97, the operating frequency of the fish finder is 200kHz;

[0036] TS = 23.61 × log (SL) - 93.78, the working frequency of the fish finder is 430kHz,

[0037] Among them, SL is the standard body length of fish.

[0038] Preferably, according to the horizontal acoustic target intensity-fish body length conversion equation determined in step S5, the fish body length is converted from the target intensity, and then the fish body weight can be converted according to the fish body length and weight relationship:

[0039] W=a×Lb (7)

[0040] Where, W is the weight of the fish, in grams; L is the length of the fish, in millimeters; a, b are constants, obtained through measurement and regression analysis of fish samples, or using existing empirical values ​​for the research object;

[0041] Then, the total fish biomass in the water area is calculated based on the fish density obtained in S4 and the water volume of the assessment water area.

[0042] Preferably, if the water body is small and the fish distribution density is relatively uniform, the total resource amount is estimated directly from the fish density and water volume obtained in step 4; if the water area is large and the fish distribution density varies greatly, the cross-section method or square area method is required to estimate the fish biomass of each sub-area, and then the total fish biomass of the entire water area is accumulated:

[0043] In the cross-section method, the cross-section observation value represents the average value of the water area between the two half cross-sections on both sides of the cross-section. The sum of the resource volume of the water area represented by each cross-section is the total resource volume within the monitoring range. Suppose the resource tail number and biomass of the assessed species i in the water area represented by a given cross-section are:

[0044]

[0045]

[0046] In the formula, —The average integral value of assessment category i within the section, in square meters / square kilometer; —The average acoustic scattering cross section of the assessment category i within the section, in square meters; D—section length, in kilometers; S—section spacing, in kilometers; - mean body weight of assessment category i, in grams;

[0047] The square area method divides the entire monitoring area into several small square areas, and the calculation is performed with the square area as the unit. The sum of the resource quantity in each area is the total resource quantity in the monitoring area. Suppose the resource tail number and biomass of the assessed species i in the waters represented by a given square area are:

[0048]

[0049]

[0050] In the formula, —The average integral value of assessment category i in the area, in square meters / square kilometer; —The average acoustic scattering cross section of the assessment type i in the square area, in square meters; A—the water area represented by the square area, in square kilometers; - the mean body weight of assessment species i, in grams.

[0051] Compared with the prior art, the present invention has achieved the following technical effects:

[0052] The method for assessing the amount of fish resources in shallow water fishery water bodies of the present invention can be applied to various types of shallow water fishery waters with a water depth of 2 to 5 m, including lakes, reservoirs, ponds, ditches, etc., and can assess all fish distributed above the bottom of the water body. It is particularly suitable for rapid assessment of the amount of fish resources in shallow water aquaculture and breeding fishery waters, and has the outstanding advantages of being non-destructive, accurate, and efficient. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] 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 work.

[0054] Figure 1 A schematic diagram of a route of a method for assessing fish resources in a shallow-water fishery water body according to an embodiment of the present invention;

[0055] Figure 2 It is another route schematic diagram of the method for assessing fish resources in shallow water fisheries in an embodiment of the present invention;

[0056] Figure 3 It is a schematic diagram of the principle of the water layer division step of the approach and avoidance effect analysis in an embodiment of the present invention;

[0057] Figure 4 Schematic diagram of the principle of the step of measuring the intensity of a horizontal acoustic target in an embodiment of the present invention;

[0058] Among them: 1-transducer, 2-net cage, 3-float, 4-sinker. DETAILED DESCRIPTION

[0059] 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.

[0060] The purpose of the present invention is to provide a method for assessing the amount of fish resources in shallow water fishery water bodies, so as to solve the problems existing in the prior art and enable the amount of fish resources in shallow water aquaculture fishery water bodies to be accurately assessed.

[0061] 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.

[0062] like Figures 1 to 4 As shown: This embodiment provides a method for assessing the amount of fish resources in a shallow water fishery body, which specifically includes the following steps:

[0063] S1, design the sampling route, determine the trajectory, number of sections and spacing of the assessment route according to the coefficient of variation of the resource assessment results that can be accepted within the estimated assessment requirements;

[0064] S2, collect acoustic data, use the horizontal acoustic detection method, and conduct navigation acoustic survey sampling according to the designed sampling route;

[0065] S3, determining the density status, determining whether the fish density distribution status is in a discrete distribution state, when it is determined that the target is in a discrete distribution state, the fish target intensity observed in situ is used for fish body length conversion;

[0066] S4, analyze the approach-avoidance effect, use the theoretical model to correct the sampling bias caused by the fish approach-avoidance effect, and estimate the actual average density of fish;

[0067] S5, determine the acoustic target intensity-body length conversion equation, obtain the equation using field measurements, or directly use the existing empirical formula;

[0068] S6, estimate the size and biomass, use the empirical formula of horizontal acoustic target intensity-body length and the body length-weight relationship to calculate the average size of fish, and use the cross-section method to estimate the total biomass in the assessed waters.

[0069] As an optional solution, in S1 of this embodiment, the sampling route adopts a square wave or triangular wave route, and the spacing between adjacent parallel routes depends on the requirements of the assessment time and assessment accuracy; the cross-section spacing of the sampling route adopting a triangular wave route is half of the spacing between the starting and ending points of a broken line. To ensure the assessment accuracy, the coefficient of variation of the resource assessment results that can be accepted within the assessment requirements is estimated according to the following formula, so as to determine the number and spacing of the sections of the assessment route:

[0070]

[0071] Where, CV is the coefficient of variation of the resource assessment results that the assessment requires to be acceptable, that is, the ratio of the standard deviation to the mean value; Λ is the monitoring coverage;

[0072]

[0073] In the formula, D is the length of the assessment route, in meters; A is the area of ​​the assessment area, in square meters; according to the requirements, CV is less than 0.25 and Λ is greater than 4, and then the length of the assessment route is determined, and then the number of sections and spacing of the assessment route are calculated from the area of ​​the assessment area. Generally, according to the requirement that Λ is greater than 4, the length of the assessment sample line must be greater than 400m in 1 hectare of water.

[0074] As an optional solution, in S2 of this embodiment, two transducers 1 are symmetrically arranged on both sides of the ship's side, the depth of the transducer 1 immersed in water is at least 0.5m, and the horizontal directions of the transducer 1 are opposite and perpendicular to the direction of the ship's travel; before the evaluation begins, the scientific fish finder system is acoustically calibrated using the standard sphere method with reference to the instrument operation manual, and the calibration of the split-beam scientific fish finder system includes the system transceiver signal gain and beam parameter calibration.

[0075] As an optional solution, in S2 of this embodiment, during the acoustic data cruise sampling process, the set sampling parameters must remain consistent, the speed is less than 5km / h, the pulse length is set according to the distribution characteristics and water depth of the evaluation object, and the data collection and recording range is set according to the water depth of the evaluation waters; the pulse length selects a short pulse less than 0.2ms, and the data collection range is greater than 10m. The data collection range is the range of data collection and recording in the horizontal direction. After the sound wave is emitted in the horizontal direction, it will form a reflected wave when it encounters the fish body. Since the beam shape is conical like a flashlight, the data collection and recording range should not be too small, and should be appropriately larger to increase the volume of the sampled water body; but it cannot be too large, because the expansion of the sound wave beam will encounter the water surface and the bottom of the water to form a reflected wave, thereby interfering with the fish body signal. In the case of a water depth of 2-5m, it is generally recommended that the data collection range is greater than 10m. Of course, if the water body is deeper, the data collection and recording range can be appropriately larger.

[0076] As an optional solution, in S3 of this embodiment, when the fish cluster distribution in the evaluation water area leads to too high density, an acoustic shading effect will be generated, resulting in a nonlinear relationship between the volume scattering intensity and the fish density, which in turn leads to analysis errors. The following formula (3) is used to calculate the Nv value to determine whether the target is in a discrete (i.e., non-clustered) distribution state:

[0077]

[0078] Where c is the speed of sound in meters per second; τ is the pulse length in seconds; Ψ is the decibel representation of the effective beam solid angle in decibels; R is the distance between the target and the surface of the transducer 1 in meters; ρ v —Volume fish density, ind. / m 3 , calculated by the following formula (4);

[0079]

[0080] In the formula, s v —Volume scattering coefficient, in m 2 / m 3 ; σ bs —Backscattering area, in square meters, S v —Volume scattering intensity, in decibels, derived from the acoustic data post-processing software in S2; TS—Fish target intensity, in decibels, derived from the acoustic data post-processing software in S2;

[0081] When Nv<0.1, that is, the distribution density of fish clusters is not high and the targets are in a discrete distribution state, the in-situ observed fish target intensity is used to convert the fish body length; when Nv>0.1, that is, the distribution density of fish clusters is high and the targets are in a non-discrete distribution state, the relevant analysis units or areas need to be excluded from the analysis.

[0082] As an optional solution, in the acoustic sampling survey process of this embodiment, due to the avoidance behavior of fish to the detection activity, the spatial distribution density of fish will deviate from the normal state. It is necessary to use a theoretical model to correct the deviation according to the distribution of fish density to obtain the correct fish density value; first, the water layer is divided symmetrically on both sides with the acoustic sampling track as the center, starting from 2m away from the track, and several water layers are divided with a spacing of 1m. Analyze whether the fish density of each water layer is uniform or there is a density gradient. When the fish distribution in each water layer is uniform, it is confirmed that the fish avoidance effect is not obvious, and the fish biomass in the water area is directly estimated by the density in S3. The fish density is analyzed from the fishery sonar acoustic data post-processing software, that is, after the fish finder collects the data, there is a dedicated data analysis software, and the fish density can be obtained through software analysis. By variance analysis, it is judged whether there is a difference in the fish density of each water layer. If there is no difference, it is considered that the fish density is uniform, and the total biomass can be directly calculated based on this mean. On the contrary, if the variance analysis shows that the density of each water layer is different, it is necessary to use formula 5 to calibrate the model and estimate the total biomass by the fish density after the model correction. The density gradient is the density difference. Generally, the water layer close to the transducer has a lower density due to the influence of the approach-avoidance effect, and the farther the water layer is, the less interference it receives, so the density is higher.

[0083] As an optional solution, in this embodiment, when there is a density gradient in the distribution of fish in each water layer, the following theoretical model is required to correct the deviation caused by the sampling process and estimate the actual average density of fish:

[0084]

[0085] Where, d is the distance from the water layer to the transducer 1, in meters; δa d —The sound density obtained in the water layer corresponding to this distance, ind. / km 2 , δα can be obtained from the d value; l—the total number of water layers divided, which is a positive integer; δ max —The maximum density recorded in each water layer, ind. / km 2 ;

[0086] Therefore, the actual average density of fish in the water body δe is:

[0087] δe=δa / P (6)

[0088] As an optional solution, the field measurement method for determining the acoustic target intensity-fish body length conversion equation in S5 in this embodiment is as follows:

[0089] First, make a net box 2 with an open top. The net is made of polyethylene monofilament with a mesh size of 1cm-2cm. A float 3 is installed on the top of the net box 2 and a sinker 4 is installed on the bottom, so that the net box 2 can float on the water surface and has a regular cube shape after stretching. The net box 2 preferably has a length×width×depth of 1.5m×1.5m×1.5m. The float and sinker can be cylinders with holes in the middle for ropes, or they can be spheres. According to the biological characteristics of the studied fish, the appropriate sample size and quantity are determined, and samples are taken at equal intervals between 10 cm and the maximum recorded body length of the fish, and the body length interval is between 2 cm and 20 cm, and the sample quantity is between 10 and 30 fish; the acoustic target intensity of fish samples of different sizes is measured horizontally from the side of the cage 2, and the number of clear and identifiable single target echo signals obtained from the measurement of each sample should be more than 300; a regression equation between the horizontal acoustic target intensity of fish and the body length of the fish is established, and the standard form of the equation is: TS = 20.0 × log (BL) - b20, where b20 is a constant;

[0090] The existing empirical formulas for determining the acoustic target intensity-body length conversion equation include: the conversion formula for the horizontal acoustic target intensity-fish body length of the common freshwater aquaculture fish carp is:

[0091] TS = 24.63 × log (SL) - 93.97, the operating frequency of the fish finder is 200kHz;

[0092] TS = 23.61 × log (SL) - 93.78, the working frequency of the fish finder is 430kHz,

[0093] Among them, SL is the standard body length of fish.

[0094] As an optional solution, in this embodiment, according to the horizontal acoustic target intensity-fish body length conversion equation determined in step S5, the fish body length is converted from the target intensity, and then the fish body weight can be converted according to the fish body length and weight relationship:

[0095] W=a×Lb (7)

[0096] Where, W is the weight of the fish, in grams; L is the length of the fish, in millimeters; a, b are constants, obtained through measurement and regression analysis of fish samples, or using existing empirical values ​​for the research object;

[0097] Then, the total fish biomass in the water area is calculated based on the fish density obtained in S4 and the water volume of the assessment water area.

[0098] As an optional solution, in this embodiment, if the water body is small and the fish distribution density is relatively uniform, the total resource amount can be directly estimated from the fish density and water volume obtained in step 4; if the water area is large and the fish distribution density varies greatly, the cross-section method or square area method is required to estimate the fish biomass of each sub-area, and then the total fish biomass of the entire water area is accumulated:

[0099] In the cross-section method, the cross-section observation value represents the average value of the water area between the two half cross-sections on both sides of the cross-section. The sum of the resource volume of the water area represented by each cross-section is the total resource volume within the monitoring range. Suppose the resource tail number and biomass of the assessed species i in the water area represented by a given cross-section are:

[0100]

[0101]

[0102] In the formula, —The average integral value of assessment category i within the section, in square meters / square kilometer; —The average acoustic scattering cross section of the assessment category i within the section, in square meters; D—section length, in kilometers; S—section spacing, in kilometers; - mean body weight of assessment category i, in grams;

[0103] The square area method divides the entire monitoring area into several small square areas, and the calculation is performed with the square area as the unit. The sum of the resource quantity in each area is the total resource quantity in the monitoring area. Suppose the resource tail number and biomass of the assessed species i in the waters represented by a given square area are:

[0104]

[0105]

[0106] In the formula, —The average integral value of assessment category i in the area, in square meters / square kilometer; —The average acoustic scattering cross section of the assessment type i in the square area, in square meters; A—the water area represented by the square area, in square kilometers; - the mean body weight of assessment species i, in grams.

[0107] Take the evaluation of fish stocks in a rectangular pond (100m×50m×2m) as an example, its surface area is 5000m 2 , volume is 10000m 3 , the main breeding target is largemouth bass:

[0108] 1) Design the cross-sectional observation sample line length to be 300m so that the monitoring coverage (Λ) is greater than 4;

[0109] 2) Use horizontal acoustic detection method to conduct underway acoustic survey sampling;

[0110] 3) The fish density distribution is determined to be Nv<0.1, and the target is determined to be in a discrete distribution state. The in-situ observed fish target intensity can be used to convert the fish body length;

[0111] 4) The theoretical model was used to correct the sampling bias caused by the fish avoidance effect, and the actual average density of fish was estimated to be 4.9×10 6 ind. / km 2 ;

[0112] 5) The average size of fish was calculated to be 987.1 g using the empirical formula of horizontal acoustic target intensity-body length and the relationship between body length and weight;

[0113] 6) The total biomass in the pond was estimated to be 24,700 kg using the cross-section method.

[0114] The method for assessing the amount of fish resources in shallow fishery water bodies of this embodiment can be applied to various types of shallow fishery waters with a water depth of 2 to 5 meters, including lakes, reservoirs, ponds, ditches, etc., and can assess all fish distributed above the bottom of the water body. It is particularly suitable for rapid assessment of the amount of fish resources in shallow water aquaculture and fishery waters, and has outstanding advantages such as being non-destructive, accurate, and efficient.

[0115] The present specification uses specific examples to illustrate the principles and implementation methods of the present invention. The above examples are 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 assessing fish resources in shallow water fisheries, characterized in that: The specific steps include: S1, design the sampling route, determine the trajectory, number of sections and spacing of the assessment route according to the coefficient of variation of the resource assessment results that can be accepted within the estimated assessment requirements; S2, collect acoustic data, use the horizontal acoustic detection method to conduct navigation acoustic survey sampling according to the designed sampling route; S3, determining the density status, determining whether the fish density distribution status is in a discrete distribution state, when it is determined that the target is in a discrete distribution state, the fish target intensity observed in situ is used for fish body length conversion; S4, analyze the approach-avoidance effect, use the theoretical model to correct the sampling deviation caused by the fish approach-avoidance effect, and estimate the actual average density of fish; when there is a density gradient in the distribution of fish in each water layer, the following theoretical model is needed to correct the deviation caused by the sampling process and estimate the actual average density of fish δe. The fish density correction coefficient P needs to be calculated first: Where, d is the distance from the water layer to the transducer, in meters; n is the distance from the farthest water layer from the transducer, in meters; δa d —The fish density in the water layer corresponding to the distance, ind. / km 2 ; l—the total number of water layers divided, which is a positive integer; δ max —The maximum density recorded in each water layer, ind. / km 2 ; Therefore, the actual average density of fish in the water body δe is: δe=δa / P (6) Where P is the fish density correction coefficient calculated based on the fish density values ​​of each divided water layer; δa is the fish density of each water layer detected, ind. / km 2 ; δe is the fish density after theoretical calculation and correction, ind. / km 2 ; S5, determine the acoustic target intensity-body length conversion equation, obtain the equation using field measurements, or directly use the existing empirical formula; S6, estimate the size and biomass, use the empirical formula of horizontal acoustic target intensity-body length and the body length-weight relationship to calculate the average size of fish, and use the cross-section method or square area method to estimate the total biomass in the assessed waters.

2. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: In S1, the sampling route adopts a square wave or a triangular wave route, and the spacing between adjacent parallel routes depends on the requirements of the assessment time and assessment accuracy. To ensure the assessment accuracy, the coefficient of variation of the resource assessment results acceptable within the assessment requirements is estimated according to the following formula, so as to determine the number and spacing of the sections of the assessment route: Where, CV is the coefficient of variation of the resource assessment results that the assessment requires to be acceptable, that is, the ratio of the standard deviation to the mean value; Λ is the monitoring coverage; Where D is the length of the evaluation route, in meters; A is the area of ​​the evaluation area, in square meters. According to the requirements, CV is less than 0.25 and Λ is greater than 4, and then the length of the evaluation route is determined. Then, the number of sections and spacing of the evaluation route are calculated from the area of ​​the evaluation area.

3. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: In S2, two transducers are symmetrically arranged on both sides of the ship's side, the transducers are immersed in water to a depth of at least 0.5m, and the horizontal directions of the transducers are opposite and perpendicular to the direction of the ship's travel; before the evaluation begins, the scientific fish finder system is acoustically calibrated using the standard sphere method with reference to the instrument operation manual. The calibration of the split-beam scientific fish finder system includes the system transceiver signal gain and beam parameter calibration.

4. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: In S2, during the acoustic data underway sampling process, the set sampling parameters must remain consistent, the speed is lower than 5km / h, the pulse length is set according to the distribution characteristics and water depth of the assessment object, and the data collection and recording range is set according to the water depth of the assessment water area; the pulse length is selected as a short pulse less than 0.2ms, and the data collection range is greater than 10m.

5. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: In S3, when the fish cluster distribution in the assessment water area leads to too high density, an acoustic shading effect will occur, resulting in a nonlinear relationship between the volume scattering intensity and the fish density, which in turn leads to analysis errors. The Nv value is calculated using the following formula (3) to determine whether the target is in a discrete distribution state: Where c is the speed of sound in meters per second; τ is the pulse length in seconds; Ψ is the decibel representation of the effective beam solid angle in decibels; R is the distance between the target and the transducer surface in meters; ρ is v —Volume fish density, ind. / m 3 , calculated by the following formula (4); In the formula, s v —Volume scattering coefficient, in m 2 / m 3 ; σ bs —Backscattering area, in square meters, S v — volume scattering intensity in decibels, derived from the acoustic data post-processing software in S2; TS—fish target strength, in decibels, derived from the acoustic data post-processing software in S2; When Nv<0.1, the target is in a discrete distribution state, and the in situ observed fish target intensity is used to convert the fish body length; when Nv>0.1, the target is in a non-discrete distribution state, and the relevant analysis units or areas need to be excluded from the analysis.

6. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: During the acoustic sampling survey, the fish avoidance behavior towards detection activities will cause the spatial distribution density of fish to deviate from the normal state. It is necessary to use a theoretical model to correct the deviation according to the fish density distribution to obtain the correct fish density value. First, the water layer is divided symmetrically on both sides with the acoustic sampling track as the center. Starting from 2m away from the track, several water layers are divided with a spacing of 1m. It is analyzed whether the fish density in each water layer is uniform or there is a density gradient. When the fish distribution in each water layer is uniform, it is confirmed that the fish avoidance effect is not obvious, and the fish biomass in the water area is directly estimated using the density in S3.

7. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: The field measurement method for determining the acoustic target intensity-fish body length conversion equation in S5 is as follows: First, a net cage with an opening on the top is made. The net is made of polyethylene monofilament and the mesh size is 1cm-2cm. A float is installed on the top of the net cage and a sinker is installed on the bottom, so that the net cage can float on the water surface and has a regular square shape after being stretched. According to the biological characteristics of the fish to be studied, the appropriate sample specifications and quantity are determined, and samples are taken at equal intervals between 10cm and the maximum recorded body length of the fish. The body length interval is between 2cm-20cm, and the number of samples is between 10 and 30. The acoustic target intensity of fish samples of different specifications is measured horizontally from the side of the net cage, and the number of clear and identifiable single target echo signals obtained from the measurement of each sample should be more than 300. A regression equation between the horizontal acoustic target intensity of fish and the body length of the fish is established, and the standard form of the equation is: TS=20.0×log(BL)-b20, where b20 is a constant. The existing empirical formulas for determining the acoustic target intensity-body length conversion equation include: the conversion formula for the horizontal acoustic target intensity-fish body length of the common freshwater farmed fish carp is: TS = 24.63 × log (SL) - 93.97, the operating frequency of the fish finder is 200kHz; TS = 23.61 × log (SL) - 93.78, the working frequency of the fish finder is 430kHz, Among them, SL is the standard body length of fish.

8. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: According to the horizontal acoustic target intensity-fish body length conversion equation determined in step S5, the fish body length is converted from the target intensity, and then the fish body weight can be converted according to the fish body length and weight relationship: W=a×L b (7) Where, W is the weight of the fish, in grams; L is the length of the fish, in millimeters; a, b are constants, obtained through measurement and regression analysis of fish samples, or using existing empirical values ​​for the research object; Then, the total fish biomass in the water area is calculated based on the fish density obtained in S4 and the water volume of the assessment water area.

9. The method for assessing fish resources in shallow water fisheries according to claim 1, characterized in that: If the water body is small and the fish distribution density is relatively uniform, the total resource volume can be directly estimated from the fish density and water volume obtained in step 4; if the water area is large and the fish distribution density varies greatly, the cross-section method or square area method is needed to estimate the fish biomass of each sub-area, and then the total fish biomass of the entire water area is accumulated: In the cross-section method, the cross-section observation value represents the average value of the water area between the two half cross-sections on both sides of the cross-section. The sum of the resource volume of the water area represented by each cross-section is the total resource volume within the monitoring range. Suppose the resource tail number and biomass of the assessed species i in the water area represented by a given cross-section are: In the formula, —The average integral value of assessment category i within the section, in square meters / square kilometer; —The average acoustic scattering cross section of the assessment category i within the section, in square meters; D—section length, in kilometers; S—section spacing, in kilometers; - mean body weight of assessment category i, in grams; The square area method divides the entire monitoring area into several small square areas, and the calculation is performed with the square area as the unit. The sum of the resource quantity in each area is the total resource quantity in the monitoring area. Suppose the resource tail number and biomass of the assessed species i in the waters represented by a given square area are: In the formula, —The average integral value of assessment category i in the area, in square meters / square kilometer; —The average acoustic scattering cross section of the assessment type i in the square area, in square meters; A—the water area represented by the square area, in square kilometers; - the mean body weight of assessment species i, in grams.

Images