AU2020103130A4 - Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology - Google Patents

Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology Download PDF

Info

Publication number
AU2020103130A4
AU2020103130A4 AU2020103130A AU2020103130A AU2020103130A4 AU 2020103130 A4 AU2020103130 A4 AU 2020103130A4 AU 2020103130 A AU2020103130 A AU 2020103130A AU 2020103130 A AU2020103130 A AU 2020103130A AU 2020103130 A4 AU2020103130 A4 AU 2020103130A4
Authority
AU
Australia
Prior art keywords
fish
model
habitat
stocks
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
AU2020103130A
Inventor
Quan QUAN
Bing Shen
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Original Assignee
Xian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Technology filed Critical Xian University of Technology
Priority to AU2020103130A priority Critical patent/AU2020103130A4/en
Application granted granted Critical
Publication of AU2020103130A4 publication Critical patent/AU2020103130A4/en
Ceased legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G09EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
    • G09BEDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
    • G09B25/00Models for purposes not provided for in G09B23/00, e.g. full-sized devices for demonstration purposes
    • G09B25/08Models for purposes not provided for in G09B23/00, e.g. full-sized devices for demonstration purposes of scenic effects, e.g. trees, rocks, water surfaces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • AHUMAN NECESSITIES
    • A01AGRICULTURE; FORESTRY; ANIMAL HUSBANDRY; HUNTING; TRAPPING; FISHING
    • A01KANIMAL HUSBANDRY; AVICULTURE; APICULTURE; PISCICULTURE; FISHING; REARING OR BREEDING ANIMALS, NOT OTHERWISE PROVIDED FOR; NEW BREEDS OF ANIMALS
    • A01K61/00Culture of aquatic animals
    • A01K61/10Culture of aquatic animals of fish
    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02BHYDRAULIC ENGINEERING
    • E02B3/00Engineering works in connection with control or use of streams, rivers, coasts, or other marine sites; Sealings or joints for engineering works in general
    • EFIXED CONSTRUCTIONS
    • E02HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
    • E02BHYDRAULIC ENGINEERING
    • E02B8/00Details of barrages or weirs ; Energy dissipating devices carried by lock or dry-dock gates
    • E02B8/08Fish passes or other means providing for migration of fish; Passages for rafts or boats
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/25Design optimisation, verification or simulation using particle-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/60Ecological corridors or buffer zones
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/80Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in fisheries management
    • Y02A40/81Aquaculture, e.g. of fish

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Educational Technology (AREA)
  • Business, Economics & Management (AREA)
  • Geology (AREA)
  • Educational Administration (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Farming Of Fish And Shellfish (AREA)

Abstract

The invention discloses a fish habitat diagnosis method based on fish individual dynamic simulation technology, which includes the following steps: (1) observation and experiment of fish individual behavior; (2) simulation of fish individual behavior; (3) habitat diagnosis analysis method of target fish stocks. Based on the field parallel observation and experimental data of target fish stocks, the correlation between fish behavior and environmental factors was obtained, and it was transformed into a mathematical function to realize the simulation of fish behavior based on individual movement trajectory of target fish stocks, and the habitat ecology and environment of target fish stocks were diagnosed according to the comparison of individual trajectory tracking and simulation results of target fish stocks. With the help of a small-volume and high-precision underwater acoustic tracking system the invention uses dynamic mathematical model and physical entity model for auxiliary and verification, and uses fish behavior as the judgment basis for fish habitat diagnosis, thus making the habitat diagnosis method more rigorous, reasonable and has more ecological significance.

Description

Habitat Identification Method Based on Fish Individual Dynamic Simulation
Technology
TECHNICAL FIELD
[01] The invention relates to a distinguishing method of fish habitat, in particular to a fish habitat identification method based on dynamic simulation of fish individuals.
BACKGROUND
[02] The construction of water conservancy project may change the river's hydrological situation, hydrodynamic force and natural state of water environment, and make the spatial distribution of river aquatic organisms habitat change to different extent. River aquatic organisms habitat can provide the complete habitat space and the environmental conditions needed for the complete ecological and life-history process for fish in the basin, and lay a solid foundation for maintaining the integrity of river ecological processes in the region, maintaining a certain amount of aquatic biological population resources, and maintaining the species and genetic diversity of aquatic organisms. In aquatic biological system, the relationship between fish and human is the most close, and fish is also the top living creature of river ecosystem. The change of fish stocks is closely related to river ecosystem. Many studies show that the change degree of fish diversity can reflect the health status of river system.
[03] In China, the protection of aquatic organisms and their habitats is mainly in the form of nature reserve construction or fish germplasm resource reserve construction, as well as corresponding protection and management schemes. In recent years, with the need of social and economic development, a large number of water conservancy and hydropower projects have been built, which affect fish habitat and make many wild fish in a state of stress. In view of this situation, people also protect fish habitat through various methods, and the establishment of fish nature reserve is the best measure to protect the fish habitat. Therefore, it is necessary to predict the change of fish habitat before and after the construction of water conservancy project reasonably and scientifically.
[04] Because of the high complexity and non-linearity of the ecosystem, the research on its evolution mechanism is often semi-empirical, while the ecological monitoring data are generally sparse. Therefore, the traditional ecological models are mainly based on semi-empirical aggregates, which usually ignore the behavior of individuals and their adaptation to the environment. Therefore, we should strengthen the investment in the key monitoring and research of fish behavior, systematically carry out long-term monitoring work, and use advanced acoustic marker positioning technology to track the fish migration trajectory, and real-time monitor the changes of river ecological factors and fish physiological and ecological behavior. However, with the help of small volume and high precision underwater acoustic tracking system, the water dynamic mathematical model and physical entity model are used for auxiliary and verification, and the results of fish habitat research based on fish behavior have not been reported. In view of the above problems, on the basis of the traditional ecological model to diagnose the habitat quality of target fish stocks, the present invention introduces the method of simulating and verifying the habitat distribution range of fish by the combination of fish individual behavior simulation and fish individual marker, and provides a more effective and scientific habitat diagnosis method.
SUMMARY
[05] The purpose of the invention is to make up for the lack offish physiology and behavior in the existing habitat identification methods, thus providing a habitat identification method based on the dynamic simulation technology of fish individuals, which can scientifically evaluate the habitat changes quantitatively in combination with the elements of fish life history, and provide a basis for the protection and restoration of fish habitat, and improve the rationality of habitat identification.
[06] In order to achieve the above purpose, the invention adopts the following technical scheme:
[07] One kind of fish habitat diagnosis method based on fish individual dynamic simulation technology, including the following steps:
[08] (1) Observation and experiment of fish individual behavior;
[09] (2) Simulation of fish individual behavior;
[010] (3) Habitat diagnosis and analysis of target fish stocks.
[011] As a further scheme of the invention, the step (1) observation and experiment of fish individual behavior includes the following steps:
[012] (1-1) Determine the target fish stocks: First, investigate the status quo of fish resources in the study area, then determine the target fish stocks in the study area based on the historical data and visiting investigation, and obtain the variation law of the population and quantity of the target fish stocks;
[013] (1-2) Trajectory tracking: Based on the fish acoustic telemetry tracer system, a monitoring network is set up to study the three-dimensional motion trajectory of the target fish in different flow fields, pressure fields (water depth) and temperature fields through the monomer calibration (sensor implantation) and tracking during releasing, and establish the relationship between the occurrence frequency, velocity, residence time of the target fish and water temperature, water velocity, water pressure, longitude and latitude, circadian rhythm environment factors.
[014] (1-3) Relation establishment: Based on the outstanding ability of generalized additive model in solving the highly nonlinear and non-monotonic relationship between response variables and prediction factors, the response model of target fish stocks and water environment factors is established, and the relationship between data response variables and prediction factors is analyzed.
[015] As a further scheme of the invention, the step (2) simulation of fish individual behavior includes the following steps:
[016] (2-1) Model construction: Using shallow water equation to construct two dimensional hydrodynamic model of the river reach studied, simulating water depth, vertical average velocity and flow field velocity. Based on the two-dimensional water hydrodynamic model, the convection-diffusion equation is added in two-dimensional water environment model to consider the heat exchange process and simulate the change of key environmental factors, and the result of the model are used as input. According to the response relationship of fish to water environment factors established in step (1-3), the Eulerian-Lagrangian method is used to establish the simulation model of the target fish, and the real-time tracking trajectory of the target fish is used to provide the velocity at any time for the model to determine the number of the model parameters, calculate the position of the target fish in the next moment, simulate the individual behavior and growth state of fish, and propose an improved simulation model of fish individual behavior. On this basis, the particle swarm is established and the different properties of each particle are given, so that all the individuals can move through the above rules, the spatial distribution of the fish in the whole river section can be obtained with the dynamic change of the water environment condition, and the precise simulation of the growth, survival and reproduction behavior of the fish stocks in the study area can be realized, and the relationship between the real fish individual behavior and the stocks distribution can be reproduced accurately. The governing equations of the model mainly include:
[017] Water flow continuity equation:
S+ +Ow S
[018] ax ay Oz (1);
[019] Flow momentum equation (xdirection):
au au2 avu awu aq I OP, _ ga17OPd _+ + + =fv-g 1 dz at ax ay az ax p0 ax P0 ax 1 (as "?+ asJ, ''-F±+a (au vI+uS
[020] Poh ax ay " z ' z (2);
[021] Flow momentum equation (ydirection):
av _+ a2 + auv + awv =-fu-g a 7 A IOa 1 gaapa J dz at ay ax az ay p 0 ay P0 a
I asV as av
[022] Poh ax ay uz az (3);
[023] Temperature convection diffusion equation:
[024] or 2 Our ± Ov+ ± Ow r + a nDw+H =F,± O +s 0 (4); +
at Ox DyV Oz az ( Oz )
Dh Fr z= a a Dh +T_
[025] aTa aj (5);
[026] Water quality control equation:
ac ac ac a2C a 2C - +u +v- D +D
[027] at ax ay Y (6);
[028] Water quality degradation equation:
-9C KC
[029] 9t (7);
[030] In the above formula, t is time; x, y, and z are Cartesian coordinates; iis the
height of the water surface; d is the still water depth; h = 7 + d is the total water depth;
u, v and w are velocity components in the directions of x, y and z, respectively; f = 2Q sin# is the Coriolis parameter (Q is the angular velocity of rotation, # is the
latitude); g is the acceleration of gravity; p is the density of water; Sxx, Sxy, Syx and Syy are components of radiation stress tensor; vt is the vertical eddy viscosity coefficient; p, is the atmospheric pressure; po is the reference density of water; S is the flow rate
of the point source, and (us, vs) is the flow rate of the source and sink items; (F,, F,) is a horizontal stress term, which is described by pressure gradient correlation and T is water temperature.(Dh, D,) are the temperature diffusion coefficients in the horizontal
direction and the vertical direction respectively; H is the source term from atmospheric heat exchange; So is the other temperature source term; C: concentration in mg/l; D, Dy: is diffusion coefficients in x and y directions, K: attenuation coefficients, in units of s-1;
[031] (2-2) Model applicability: The model constructed by step (2-1) is applicable to most of the fish stocks in the river currently studied, and it is not limited to the target fish stocks. The reason is that the building module of the fish individual dynamic model is based on the combination of the basic model and the fish individual model. The basic model is the simulation construction of terrain, hydrology, water environment and hydraulic field, and the fish individual model is based on the Eulerian-Lagrangian framework. In order to obtain the growth, foraging and breeding behavior of the target fish stocks with different fish attributes, the applicability of the dynamic model of the individual fish stocks is mainly embodied in the transformation of the living habits of the target fish stocks into a mathematical function, and the rule is summarized. The growth function of the fish is derived by mathematical statistics, and the process of the fish movement is induced by fluid hydrodynamics and river dynamics.
[032] The growth of fish includes body weight growth and body length growth. Since the growth of fish exists in non-uniform linear in the life cycle, in order to describe the growth process of fish, mathematical equations can be used to describe the growth characteristics of fish, and von Bertalanffy growth equation based on metabolism theory is adopted.
It = L.[1 - e-k(t- t o )]
[033] (8);
Wt = W.[1 - e-k(t- to)]3
[034] (9);
[035] Where t is the age, the unit is day (d), It and Wt are the average body length
(cm) and body weight (kg) at time t, I is the average progressive body length (cm) and
W 1 is the average progressive body weight (kg), k is the growth coefficient (1/d), and to is the hypothetical theoretical growth starting age;
[036] The motion process in the fish individual model is actually the embodiment of Lagrangian algorithm in the basic model. It has accurate position (x, y, z), and is independent of the whole model network structure. The direction of motion, velocity and state variables (volume, mass) can be obtained by assignment. At the same time, the rules such as birth and death can be defined, and the feedback information between Euler and Langrange can be obtained.
[037] As a further scheme of the invention, the step (3) diagnosis and analysis of the habitat of the target fish stocks includes: reproducing the habitat distribution of the target fish stocks according to the simulation model of the fish individual behavior, and studying the data distribution characteristics from the data sample itself with the acoustic telemetry trajectory within the range of the predetermined "habitat", combining the kernel density estimation function to divide the high density and low density area, and then dividing the high density area into small area. The chi-square test theory is used to test whether the data change with the regional characteristics, so as to judge the specific location of the real habitat; at the same time, based on the psychological Q-matrix theory, an improvement is made and a diagnosis scheme is put forward to evaluate the small-scale behavior of the target fish stocks. Based on the on the-spot fishing of "habitat" and the behavior trajectory of target fish, the relationship between the probability of individual occurrence and the sex, body size, date, season and environment of the target fish is quantitatively analyzed, and a habitat diagnostic classification model is constructed to determine the physical attributes of the "habitat" of the native fish, including the spawning ground, wintering ground and feeding ground.
[038] The invention has the following advantages: the invention focuses on the behavior characteristics of various high-level biological populations including fish in the river ecosystem, adopts a variety of methods, such as physiological ecology, fluid hydrodynamics, river dynamics, fuzzy mathematics, computational fluid dynamics, etc. to reasonably and scientifically determine the suitable habitat of the target fish stocks, and fills the blank of the technology of fish habitat protection and recognition, improves the research system and has great significance for river ecosystem protection.
[039] The fish individual dynamic model of the invention is universal to river fish, and different fish movement trajectories will be obtained by changing the growth and movement functions of the particle fish in the model, which is of great significance to perfect the research of river ecological flow and habitat.
[040] With the help of a small-volume and high-precision underwater acoustic tracking system the invention uses dynamic mathematical model and physical entity model for auxiliary and verification, and uses fish behaviour as the judgment basis for fish habitat diagnosis, thus making the habitat diagnosis method more rigorous, reasonable and has more ecological significance.
BRIEF DESCRIPTION OF THE FIGURES
[041] Figure 1 is a flowchart of the invention;
[042] Figure 2 shows a grid division diagram for the implementation embodiment study area model.
[043] Figure 3 is a simulation graph of individual fish trajectories (in which a is a graph of the fish aggregation effect; b is a graph of shoreline aggregation and sidewall migration trajectory of fish; c is a consistency graph of fish swimming tendency);
[044] Figure 4 is a distribution map of suitable spawning grounds by combining fish tracer map and simulated habitat map.
DESCRIPTION OF THE INVENTION
[045] The invention is further described in detail in combination with the attached drawing and the specific embodiment.
[046] As shown in the figure, a habitat diagnosis method based on dynamic simulation of individual fish in this embodiment includes the following steps:
[047] 1. Fish individual behavior observation and experimental study methods;
[048] (1-1) Determine the target fish stocks: First, investigate the status quo of the fish resources in the study area, then determine the target fish stocks in the study area, and obtain the variation law of the population and quantity of the target fish stocks;
[049] On the basis of on-the-spot investigation, visit to local residents and literature, the amount of fish resources, composition structure, distribution location and living habits in the study area were statistically summarized. Finally, the target stocks in the study area were Gymnocypris eckloni eckloni Herzenstein, which began to lay eggs in May each year. The spawning ground was located at the gravel bottom of the main stream of the Yellow River. Generally, it was located in the section where the water is clear and the water flow is relatively fast. They're demersal eggs and mature in the pit; both flowing water and still water can live, but most of the time they live in flowing water, usually scattered or concentrated in small groups in habitat, and during the breeding season they gather in large groups to the larger tributaries leading to the main stream, reservoir or lake. After hatching, the young fish are fed in the shallow water along with the flowing water into the main stream bay fork or lake and the banks of the reservoir; the main food is aquatic multi-vertebrates and they also eat small-scale loaches.
[050] (1-2) Trajectory tracking: Based on the fish acoustic telemetry tracer system, a monitoring network is set up to study the three-dimensional trajectory of the target fish in different flow field, pressure field (water depth) and temperature field; establish the relationship between the frequency, velocity and residence time of the target fish and water temperature, water velocity, water pressure, longitude, latitude and circadian rhythm.
[051] (1-3) Relationship establishment: Based on the outstanding ability of generalized additive model in solving the highly nonlinear and non-monotonic relationship between response variables and prediction factors, the response model of target fish stocks and water environment factors is established, and the relationship between data response variables and prediction factors is analyzed.
[052] By using the generalized additive model (GAM) coupling function, the functional relationship between the mathematical expectation value of the response variable and the prediction variable is established. Taking the presence or absence of the labeled sample as the observed variable, the response variable follows a binomial distribution with a value of 1 and no value of 0. The self-velocity and trajectory distribution density of the labeled samples are the dependent variables of the model, and the independent variables of the model include three sets of data: three-dimensional spatial variables, time variables and environment variables.
[053] 2. Research methods of fish individual behavior simulation;
[054] (2-1) Model construction: Using shallow water equation to build two dimensional hydrodynamic model of the river reach studied, simulating water depth, vertical average velocity and flow field velocity. Based on the two-dimensional water hydrodynamic model, the convection-diffusion equation is added in two-dimensional water environment model to consider the heat exchange process and simulate the change of key environmental factors, and the result of the model are used as input, using the established fish-water response relationship, calculate the position of the target fish in the next moment, i. e. the velocity u of the target fish at the time of t is synchronized into the Eulerian-Lagrangian formula by using the time, longitude and latitude, and the parameters are calculated and calibrated to obtain the position at the next moment:
[055] According to the response relationship of fish to water environment factors established in step (1-3), the Eulerian-Lagrangian method is used to establish the simulation model of the target fish, and the real-time tracking trajectory of the target fish is used to provide the velocity at any time for the model to determine the number of the model parameters, calculate the position of the target fish in the next moment, simulate the individual behavior and growth state of fish, and propose an improved simulation model of fish individual behavior. On this basis, the particle swarm is established and the different properties of each particle are given, so that all the individuals can move through the above rules, the spatial distribution of the fish in the whole river section can be obtained with the dynamic change of the water environment condition, and the precise simulation of the growth, survival and reproduction behavior of the fish stocks in the study area can be realized, and the relationship between the real fish individual behavior and the stocks distribution can be reproduced accurately. The governing equations of the model mainly include:
[056] Water flow continuity equation:
S+ +Ow S
[057] Ox 0y Oz (1);
[058] Flow momentum equation (x direction):
au + au-+2 - avu -+ awu -= fv-g aq 1 ap g~aJ--dz - fZ7O
at ax ay az ax p0 ax p0 ax
H-+F±+ a (au 1 (as "?+asJ, vI- uS
[059] Poh ax ay " z ' z (2):
[060] Flow momentum equation (ydirection):
[061] or + ouT + + = Fr + D + H + T, (4); at ax ay Oz az 'az
aaK aT FT = Dh + Dh
[062] ax ax ay ay (5);
[063] Water quality control equation:
ac ac ac a 2C a 2C +u-+v =D - +D
[064] at ax ay ay (6);
[065] Water quality degradation equation:
-9C KC
[066] 9t (7);
[067] In the above formula, t is time; x, y, and z are Cartesian coordinates; q is the
height of the water surface; d is the still water depth; h = q + d is the total water depth;
u, v and w are velocity components in the directions of x, y and z, respectively; f = 2Q sin# is the Coriolis parameter (Q is the angular velocity of rotation, # is the
latitude); g is the acceleration of gravity; p is the density of water; Sxx, Sxy, Syx and Syy are components of radiation stress tensor; vt is the vertical eddy viscosity coefficient; p, is the atmospheric pressure; po is the reference density of water; S is the flow rate
of the point source, and (us, vs) is the flow rate of the source and sink items; (F,, F,) is a horizontal stress term, which is described by pressure gradient correlation and T is water temperature.(Dh, D,) are the temperature diffusion coefficients in the horizontal
direction and the vertical direction respectively; H is the source term from atmospheric heat exchange; So is the other temperature source term; C: concentration in mg/l; D, Dy: is diffusion coefficients in x and y directions, K: attenuation coefficients, in units of s-1;
[068] (2-2) Model applicability: The model constructed by step (2-1) is applicable to most of the fish stocks in the river currently studied, and is not limited to the target fish stocks. The reason is that the building module of the fish individual dynamic model is based on the combination of the basic model and the fish individual model. The basic model is the simulation construction of terrain, hydrology, water environment and hydraulic field. The fish individual model is based on Eulerian-Lagrangian framework, which gives different fish attributes to the particle swarm to obtain the behavior of growth, foraging and reproduction of the target stocks. The applicability of the individual dynamic model of fish is mainly embodied in the transformation of the life habits of the target stocks into mathematical functions, summarizing the rules, the growth function of fish derived from mathematical statistics, and the process of fish movement induced by fluid hydrodynamics and river dynamics.
[069] The growth of fish includes body weight growth and body length growth. Since the growth of fish exists in non-uniform linear in the life cycle, in order to describe the growth process of fish, mathematical equations can be used to describe the growth characteristics of fish, and von Bertalanffy growth equation based on metabolism theory is adopted.
It = L [1 - e-k(t- to)]
[070] (8);
Wt = W.[1 - e k(t- to) 3
[071] (9);
[072] Where t is the age, the unit is day (d),it and Wt are the average body length
(cm) and body weight (kg) at time t, I is the average progressive body length (cm) and
W 1 is the average progressive body weight (kg), k is the growth coefficient (1/d), and to is the hypothetical theoretical growth starting age;
[073] The motion process in the fish individual model is the embodiment of Lagrangian algorithm in the basic model. It has accurate position (x, y, z), and is independent of the whole model network structure. The direction, velocity and state variables (volume, mass, etc.) can be obtained by assignment. At the same time, the rules such as birth, death, and the feedback information between Euler and Langrange can be defined.
[074] 3. Analysis methods for habitat diagnosis of target stocks:
[075] Reproducing the habitat distribution of the target fish stocks according to the simulation model of the fish individual behavior, and studying the data distribution characteristics from the data sample itself with the acoustic telemetry trajectory within the range of the predetermined "habitat", combining the kernel density estimation function to divide the high density and low density area, and then dividing the high density area into small area. The chi-square test theory is used to test whether the data change with the regional characteristics, so as to judge the specific location of the real habitat. On the basis of Q-matrix theory, an improved "habitat attribute" diagnosis scheme is proposed. The method of evaluating the small-scale behavior of the target stocks: Q matrix is Kxm matrix, K means the attributes of the target stocks, which is defined here as three behaviors of bait, breeding and overwintering; m means the items of the test, and the distribution of 16 tracks of the target stocks is defined here. Based on the on-the-spot fishing and small-scale behavior assessment of target fish, the relationship between the probability of individual occurrence and the sex, body size, date, season and environment of target fish is quantitatively analyzed. A classification model of habitat diagnosis is constructed to determine the physical attributes of the "habitat" of native fish (spawning ground, wintering ground and feeding ground).
[076] Although the invention has been described with reference to specific examples, it will be appreciated by those skilled in the art that the invention may be embodied in many other forms, in keeping with the broad principles and the spirit of the invention described herein.
[077] The present invention and the described embodiments specifically include the best method known to the applicant of performing the invention. The present invention and the described preferred embodiments specifically include at least one feature that is industrially applicable

Claims (4)

THE CLAIMS DEFINING THE INVENTION ARE AS FOLLOWS:
1. A fish habitat diagnosis method based on fish individual dynamic simulation technology is characterized in that following steps:
(1) Observation and experiment of fish individual behavior;
(2) Simulation of fish individual behavior;
(3) Habitat diagnosis and analysis of target fish stocks.
2. The method of fish habitat diagnosis based on fish individual dynamic simulation technology, as described in claim 1, is characterized by the following steps: (1) The observation and experiment of fish individual behavior includes the following steps:
(1-1) Determine the target fish stocks: First, investigate the status quo of fish resources in the study area, then determine the target fish stocks in the study area based on the historical data and visiting investigation, and obtain the variation law of the population and quantity of the target fish stocks;
(1-2) Trajectory tracking: Based on the fish acoustic telemetry tracer system, a monitoring network is set up to study the three-dimensional motion trajectory of the target fish in different flow fields, pressure fields (water depth) and temperature fields through the monomer calibration (sensor implantation) and tracking during releasing, and establish the relationship between the occurrence frequency, velocity, residence time of the target fish and water temperature, water velocity, water pressure, longitude and latitude, circadian rhythm environment factors.
(1-3) Relation establishment: Based on the outstanding ability of generalized additive model in solving the highly nonlinear and non-monotonic relationship between response variables and prediction factors, the response model of target fish stocks and water environment factors is established, and the relationship between data response variables and prediction factors is analyzed.
3. The fish habitat diagnosis method based on fish individual dynamic simulation technology, as described in claim 2, is characterized by the following steps: (2) The fish individual behavior simulation includes the following steps:
(2-1) Model construction: Using shallow water equation to construct two dimensional hydrodynamic model of the river reach studied, simulating water depth, vertical average velocity and flow field velocity. Based on the two-dimensional water hydrodynamic model, the convection-diffusion equation is added in two-dimensional water environment model to consider the heat exchange process and simulate the change of key environmental factors, and the result of the model are used as input. According to the response relationship of fish to water environment factors established in step (1-3), the Eulerian-Lagrangian method is used to establish the simulation model of the target fish, and the real-time tracking trajectory of the target fish is used to provide the velocity at any time for the model to determine the number of the model parameters, calculate the position of the target fish in the next moment, simulate the individual behavior and growth state of fish, and propose an improved simulation model of fish individual behavior. On this basis, the particle swarm is established and the different properties of each particle are given, so that all the individuals can move through the above rules, the spatial distribution of the fish in the whole river section can be obtained with the dynamic change of the water environment condition, and the precise simulation of the growth, survival and reproduction behavior of the fish stocks in the study area can be realized, and the relationship between the real fish individual behavior and the stocks distribution can be reproduced accurately. The governing equations of the model mainly include:
Water flow continuity equation:
au 0V __ + Ow =s ax ay az (1);
Flow momentum equation (x direction):
Bu -+-Bu2 +-aVu +---=fv-g- Owu 8q 1 1 a& g fZ7O@ -dz at ax ay az ax po ax P0 ax 1 as.+ asX+, a , au +uS __+ +±F±+ v Poh ax y z(' az (2);
Flow momentum equation (ydirection):
Dv Dy2 Duv Dwv aDq 1 DPa g fI17 a + + + =-fu-g D 1 Pdz at ay ax az ay po ay po ay 1 sY asY a av p+h+Dx+ " +F + - v +vS poh Ox ay z ' z(3);
Temperature convection diffusion equation:
DT7 DufT DvT DwT aD a DT~ + + =F,+ ,+ D + H+TS (4); at ax Dy az az V az)
FT= -QDh D-+- DhD T ax ax ay
( ay ()
Water quality control equation:
aC - +u aC -+v aC -= D 8 2C +D 0 2C at ax ay X ax2 Y aY 2 (6);
Water quality degradation equation:
C= -KC at (7);
In the above formula, t is time; x, y, and z are Cartesian coordinates;iis the
height of the water surface; d is the still water depth; h = q + d is the total water depth; u, v and w are velocity components in the directions of x, y and z, respectively; f = 2Q sin# is the Coriolis parameter (Q is the angular velocity of rotation, # is the latitude); g is the acceleration of gravity; p is the density of water; Sxx, Sxy, Syx and Syy are components of radiation stress tensor; vt is the vertical eddy viscosity coefficient; pa is the atmospheric pressure; po is the reference density of water; S is the flow rate of the point source, and (us, vs) is the flow rate of the source and sink items; (Fu, Fv) is a horizontal stress term, which is described by pressure gradient correlation and T is water temperature.(Dh, D) are the temperature diffusion coefficients in the horizontal direction and the vertical direction respectively; H is the source term from atmospheric heat exchange; So is the other temperature source term; C: concentration in mg/l; D,
Dy: is diffusion coefficients in x and y directions, K: attenuation coefficients, in units of s-1;
(2-2) Model applicability: The model constructed by step (2-1) is applicable to most of the fish stocks in the river currently studied, and it is not limited to the target fish stocks. The reason is that the building module of the fish individual dynamic model is based on the combination of the basic model and the fish individual model. The basic model is the simulation construction of terrain, hydrology, water environment and hydraulic field, and the fish individual model is based on the Eulerian-Lagrangian framework. In order to obtain the growth, foraging and breeding behavior of the target fish stocks with different fish attributes, the applicability of the dynamic model of the individual fish stocks is mainly embodied in the transformation of the living habits of the target fish stocks into a mathematical function, and the rule is summarized. The growth function of the fish is derived by mathematical statistics, and the process of the fish movement is induced by fluid hydrodynamics and river dynamics.
The growth of fish includes body weight growth and body length growth. Since the growth of fish exists in non-uniform linear in the life cycle, in order to describe the growth process of fish, mathematical equations can be used to describe the growth characteristics of fish, and von Bertalanffy growth equation based on metabolism theory is adopted.
It = LO.[1 - e--ktt- to)] (8):
Wt = W.[1 - e-k(t- to)1 3 (9);
Where t is the age, the unit is day (d), l and Wt are the average body length
(cm) and body weight (kg) at time t, 11 is the average progressive body length (cm) and
W 1 is the average progressive body weight (kg), k is the growth coefficient (1/d), and to is the hypothetical theoretical growth starting age;
The motion process in the fish individual model is actually the embodiment of
Lagrangian algorithm in the basic model. It has accurate position (x, y, z), and is independent of the whole model network structure. The direction of motion, velocity and state variables (volume, mass) can be obtained by assignment. At the same time, the rules such as birth and death can be defined, and the feedback information between Eulerian and Langrange can be obtained.
4. The method of fish habitat diagnosis based on the technique of fish individual dynamic simulation technology, as described in claim 2, is characterized by the following steps: (3) The diagnostic analysis of the habitat of the target fish stocks includes: reproducing the habitat distribution of the target fish stocks according to the simulation model of fish individual behavior, and studying the data distribution characteristics from the data sample itself with the acoustic telemetry trajectory within the range of the predetermined "habitat", combining the kernel density estimation function to divide the high density and low density area, and then dividing the high density area into small area. The chi-square test theory is used to test whether the data change with the regional characteristics, so as to judge the specific location of the real habitat; at the same time, based on the psychological Q-matrix theory, an improvement is made and a diagnosis scheme is put forward to evaluate the small-scale behavior of the target fish stocks. Based on the on-the-spot fishing of "habitat" and the behavior trajectory of target fish, the relationship between the probability of individual occurrence and the sex, body size, date, season and environment of the target fish is quantitatively analyzed, and a habitat diagnostic classification model is constructed to determine the physical attributes of the "habitat" of the native fish, including the spawning ground, wintering ground and feeding ground.
AU2020103130A 2020-10-30 2020-10-30 Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology Ceased AU2020103130A4 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AU2020103130A AU2020103130A4 (en) 2020-10-30 2020-10-30 Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
AU2020103130A AU2020103130A4 (en) 2020-10-30 2020-10-30 Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology

Publications (1)

Publication Number Publication Date
AU2020103130A4 true AU2020103130A4 (en) 2021-01-07

Family

ID=74041717

Family Applications (1)

Application Number Title Priority Date Filing Date
AU2020103130A Ceased AU2020103130A4 (en) 2020-10-30 2020-10-30 Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology

Country Status (1)

Country Link
AU (1) AU2020103130A4 (en)

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112753623A (en) * 2021-01-22 2021-05-07 交通运输部天津水运工程科学研究所 Experimental device and testing method for brave and exploration behaviors of fishes
CN113176328A (en) * 2021-04-25 2021-07-27 中国水产科学研究院珠江水产研究所 Method for evaluating influence of hydro-junction on fish habitat based on otolith micro-area in-situ trace elements
CN113326743A (en) * 2021-05-10 2021-08-31 大连海洋大学 Fish shoal movement behavior parameter extraction and analysis method under breeding background condition
CN113475431A (en) * 2021-07-16 2021-10-08 四川大学 Method for constructing branch habitat substitution suitability model matched with fish life history
CN113598099A (en) * 2021-08-10 2021-11-05 中国电建集团华东勘测设计研究院有限公司 Design method of artificial fish spawning site with water depth-flow velocity coupling
CN113807017A (en) * 2021-09-23 2021-12-17 西安理工大学 Method for determining fish preference habitat and terminal equipment
CN114304035A (en) * 2021-12-10 2022-04-12 大连海洋大学 Research method for influence of fish swarm movement of recirculating aquaculture system on hydrodynamics and flow field in aquaculture pond
CN114463675A (en) * 2022-01-11 2022-05-10 北京市农林科学院信息技术研究中心 Underwater fish group activity intensity identification method and device
CN114529170A (en) * 2022-01-26 2022-05-24 中国长江三峡集团有限公司 Fish salinity tolerance assessment method and device, storage medium and electronic equipment
CN114532260A (en) * 2022-01-18 2022-05-27 中国长江三峡集团有限公司 Accurate regulation and control system for breeding fish producing drifting eggs
CN114698576A (en) * 2022-05-23 2022-07-05 中国长江三峡集团有限公司 Natural propagation incubation protection device for special rare fishes at upper reaches of Yangtze river and use method
CN115081201A (en) * 2022-06-15 2022-09-20 三峡大学 Prediction method applied to assessment of culvert fish passing capacity
CN115152683A (en) * 2022-07-29 2022-10-11 广东海洋大学 Simulation test device for research on water body pressure for fish coping
CN116050304A (en) * 2023-03-15 2023-05-02 重庆交通大学 Intelligent fish flow field simulation control method, system, equipment and storage medium
CN116432902A (en) * 2023-03-31 2023-07-14 中国水利水电科学研究院 Species long-term viability assessment method considering water environment factor variation
WO2023165299A1 (en) * 2022-03-03 2023-09-07 石河子大学 Intelligent analysis system used for ethology of various high-density tiny insects
CN117541078A (en) * 2023-11-21 2024-02-09 交通运输部规划研究院 Ecological protection strategy customizing method based on artificial canal development
CN117710176A (en) * 2024-02-05 2024-03-15 中国水产科学研究院南海水产研究所 Ecological construction and restoration method and system for ecological sea wall oyster reef
CN117892980A (en) * 2024-03-14 2024-04-16 长江水资源保护科学研究所 Ecological scheduling method and device for round-mouth copper fish

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112753623B (en) * 2021-01-22 2022-03-29 交通运输部天津水运工程科学研究所 Experimental device and testing method for brave and exploration behaviors of fishes
CN112753623A (en) * 2021-01-22 2021-05-07 交通运输部天津水运工程科学研究所 Experimental device and testing method for brave and exploration behaviors of fishes
CN113176328A (en) * 2021-04-25 2021-07-27 中国水产科学研究院珠江水产研究所 Method for evaluating influence of hydro-junction on fish habitat based on otolith micro-area in-situ trace elements
CN113326743A (en) * 2021-05-10 2021-08-31 大连海洋大学 Fish shoal movement behavior parameter extraction and analysis method under breeding background condition
CN113326743B (en) * 2021-05-10 2023-10-13 大连海洋大学 Method for extracting and analyzing fish school movement behavior parameters under cultivation background condition
CN113475431B (en) * 2021-07-16 2022-05-10 四川大学 Method for constructing branch habitat substitution suitability model matched with fish life history
CN113475431A (en) * 2021-07-16 2021-10-08 四川大学 Method for constructing branch habitat substitution suitability model matched with fish life history
CN113598099A (en) * 2021-08-10 2021-11-05 中国电建集团华东勘测设计研究院有限公司 Design method of artificial fish spawning site with water depth-flow velocity coupling
CN113807017A (en) * 2021-09-23 2021-12-17 西安理工大学 Method for determining fish preference habitat and terminal equipment
US11645354B2 (en) 2021-09-23 2023-05-09 Xi'an University Of Technology Determination method for preferred habitat of fish and terminal device
CN113807017B (en) * 2021-09-23 2022-05-24 西安理工大学 Method for determining fish preference habitat and terminal equipment
CN114304035A (en) * 2021-12-10 2022-04-12 大连海洋大学 Research method for influence of fish swarm movement of recirculating aquaculture system on hydrodynamics and flow field in aquaculture pond
CN114304035B (en) * 2021-12-10 2023-04-25 大连海洋大学 Method for researching influence of shoal motion of circulating water culture system on hydrodynamic force and flow field in culture pond
CN114463675A (en) * 2022-01-11 2022-05-10 北京市农林科学院信息技术研究中心 Underwater fish group activity intensity identification method and device
CN114532260A (en) * 2022-01-18 2022-05-27 中国长江三峡集团有限公司 Accurate regulation and control system for breeding fish producing drifting eggs
CN114532260B (en) * 2022-01-18 2022-12-09 中国长江三峡集团有限公司 Accurate regulation and control system for breeding fish producing drifting eggs
CN114529170A (en) * 2022-01-26 2022-05-24 中国长江三峡集团有限公司 Fish salinity tolerance assessment method and device, storage medium and electronic equipment
CN114529170B (en) * 2022-01-26 2023-07-11 中国长江三峡集团有限公司 Fish salinity tolerance assessment method, device, storage medium and electronic equipment
US11967182B2 (en) 2022-03-03 2024-04-23 Shihezi University Intelligent analysis system applied to ethology of various kinds of high-density minimal polypides
WO2023165299A1 (en) * 2022-03-03 2023-09-07 石河子大学 Intelligent analysis system used for ethology of various high-density tiny insects
CN114698576A (en) * 2022-05-23 2022-07-05 中国长江三峡集团有限公司 Natural propagation incubation protection device for special rare fishes at upper reaches of Yangtze river and use method
CN115081201B (en) * 2022-06-15 2024-07-05 三峡大学 Prediction method applied to evaluation of culvert fish passing capability
CN115081201A (en) * 2022-06-15 2022-09-20 三峡大学 Prediction method applied to assessment of culvert fish passing capacity
CN115152683A (en) * 2022-07-29 2022-10-11 广东海洋大学 Simulation test device for research on water body pressure for fish coping
CN116050304B (en) * 2023-03-15 2024-03-26 重庆交通大学 Intelligent fish flow field simulation control method, system, equipment and storage medium
CN116050304A (en) * 2023-03-15 2023-05-02 重庆交通大学 Intelligent fish flow field simulation control method, system, equipment and storage medium
CN116432902A (en) * 2023-03-31 2023-07-14 中国水利水电科学研究院 Species long-term viability assessment method considering water environment factor variation
CN116432902B (en) * 2023-03-31 2023-09-12 中国水利水电科学研究院 Species long-term viability assessment method considering water environment factor variation
CN117541078A (en) * 2023-11-21 2024-02-09 交通运输部规划研究院 Ecological protection strategy customizing method based on artificial canal development
CN117541078B (en) * 2023-11-21 2024-05-28 交通运输部规划研究院 Ecological protection strategy customizing method based on artificial canal development
CN117710176A (en) * 2024-02-05 2024-03-15 中国水产科学研究院南海水产研究所 Ecological construction and restoration method and system for ecological sea wall oyster reef
CN117710176B (en) * 2024-02-05 2024-05-03 中国水产科学研究院南海水产研究所 Ecological construction and restoration method and system for ecological sea wall oyster reef
CN117892980A (en) * 2024-03-14 2024-04-16 长江水资源保护科学研究所 Ecological scheduling method and device for round-mouth copper fish
CN117892980B (en) * 2024-03-14 2024-05-24 长江水资源保护科学研究所 Ecological scheduling method and device for round-mouth copper fish

Similar Documents

Publication Publication Date Title
AU2020103130A4 (en) Habitat Identification Method Based on Fish Individual Dynamic Simulation Technology
CN109271694B (en) Habitat identification method based on fish individual dynamic simulation technology
Berryman Principles of population dynamics and their application
Quan et al. Assessment of the sustainability of Gymnocypris eckloni habitat under river damming in the source region of the Yellow River
Jongejans et al. Dispersal, demography and spatial population models for conservation and control management
Williams et al. Estimating habitat carrying capacity for migrating and wintering waterfowl: considerations, pitfalls and improvements
Metaxas et al. Quantifying the “bio-” components in biophysical models of larval transport in marine benthic invertebrates: advances and pitfalls
Doerr et al. Dispersal range analysis: quantifying individual variation in dispersal behaviour
Spear et al. Disease transmission models for public health decision making: toward an approach for designing intervention strategies for Schistosomiasis japonica.
Neuenfeldt et al. Analysing migrations of Atlantic cod Gadus morhua in the north‐east Atlantic Ocean: then, now and the future
Wang et al. Integrated modelling of the life cycle and aeroecology of wind-borne pests in temporally-variable spatially-heterogeneous environment
Botsford Individual state structure in population models
Michel et al. Climate change, hydrology, and fish morphology: predictions using phenotype-environment associations
Perry et al. Application of the stream salmonid simulator (S3) to Klamath River fall Chinook Salmon (Oncorhynchus tshawytscha), California—parameterization and calibration
CN110414088B (en) Space fuzzy evaluation method for suitability of birds-involved habitat by combining hydrodynamic model
Cowled et al. A review of geospatial and ecological factors affecting disease spread in wild pigs: considerations for models of foot-and-mouth disease spread
Perry et al. Model structure of the stream salmonid simulator (S3)—A dynamic model for simulating growth, movement, and survival of juvenile salmonids
Hamda et al. Applying a simplified energy-budget model to explore the effects of temperature and food availability on the life history of green sturgeon (Acipenser medirostris)
Özyalın Interpretation of volcanic magnetic anomalies using differential search algorithm: case study from the Kula volcanic park, western Türkiye
Goodwin et al. Simulating mobile populations in aquatic ecosystems
Sobocinski et al. Trophic transfer in seagrass systems: estimating seasonal production of an abundant seagrass fish, Bairdiella chrysoura, in lower Chesapeake Bay
Delaney et al. Predicting regional spread of non-native species using oceanographic models: validation and identification of gaps
Perry et al. Application of the Stream Salmonid Simulator (S3) to the restoration reach of the Trinity River, California—Parameterization and calibration
Gross et al. Hatching distribution, abundance, and losses to freshwater diversions of longfin smelt inferred using hydrodynamic and particle-tracking models
Gentleman Factors controlling the seasonal abundance and distribution of Calanus finmarchicus in the Gulf of Maine/Georges Bank region

Legal Events

Date Code Title Description
FGI Letters patent sealed or granted (innovation patent)
MK22 Patent ceased section 143a(d), or expired - non payment of renewal fee or expiry