Numerical simulation method and system of rapid gonadal maturation of Lateolabrax japonicus in complex marine system model

Abstract

The present invention discloses a method and system for numerical simulation of striped sea bass with rapid gonadal maturation in a complex sea system model, and relates to the technical field of numerical simulation of fish. The method comprises: randomly grouping striped sea bass of the same age, raising each group of test samples under different salinity water samples, and keeping other environmental parameters consistent; randomly detecting the key hormone levels and key gene expression levels of the striped sea bass in each group of test samples, and establishing a prediction model with salinity as input and gonadal development-related indicators as output; inputting the salinity of the sea area to be raised into the model, and obtaining the predicted values ​​of the key hormone levels and key gene expression levels of the striped sea bass gonadal development, and generating a preliminary gonadal development index in combination with environmental parameters; correcting the gonadal development index based on the sea area disturbance parameters, and obtaining an accurate gonadal development index to complete the simulation of the gonadal development degree of the striped sea bass. The method can accurately predict the gonadal development degree of the striped sea bass in a complex sea environment, and provide guidance for scientific breeding.

Description

Technical Field

[0001] The invention relates to the technical field of fish numerical simulation, and in particular to a method and system for numerical simulation of Lateolabrax japonicus with rapid gonadal maturation in a complex sea system model. Background Art

[0002] Lateolabrax japonicus is one of the important economic marine fish species in China. It is widely distributed in coastal waters of my country and is valued by farmers for its delicious meat, fast growth and strong adaptability. However, with the expansion of the scale of Lateolabrax aquaculture and the complexity of the ecological environment in natural sea areas, how to promote the gonadal development of Lateolabrax japonicus and improve its reproductive efficiency has become a core issue in the current aquaculture and fishery resource management. During the gonadal development of Lateolabrax japonicus, salinity, temperature, light, nutrition and other environmental disturbances (such as sea surface wind speed, ocean current circulation, etc.) have a significant impact on its gonadal maturation and reproductive results. For example, salinity fluctuations can significantly affect the secretion of key hormones such as gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone by regulating the endocrine system and gene expression of fish. Changes in these hormone levels further lead to the upregulation or downregulation of the expression levels of genes related to gonadal development (such as aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor), ultimately affecting the maturation process and quality of gonadal development.

[0003] At present, research on gonadal development in marine ecology and aquaculture, especially for fish with complex growth cycles such as striped seabass, lacks a systematic modeling method that can accurately reflect the impact of environmental changes on individual organisms. Most existing studies focus on the role of a single environmental factor, while ignoring the comprehensive effects of multiple factor interactions on the growth and gonadal development of striped seabass. Therefore, a new method is urgently needed to monitor and simulate the gonadal development process of striped seabass in real time, especially the changes under different salinities and sea conditions. This requires combining environmental changes with the growth characteristics of individual organisms to propose a method that can quickly and efficiently predict the level of gonadal development of striped seabass.

[0004] As a powerful data processing and prediction tool, neural network models also have wide application potential in the field of aquaculture, but they are currently less used, especially in the process of gonadal development, which is still in its infancy. Therefore, building a gonadal development prediction model for Lateolabrax japonicus based on neural network technology and combining it with environmental factors for real-time monitoring and simulation can effectively improve the accuracy and application scope of gonadal development prediction, thereby providing strong support for marine ecological protection and aquaculture management.

[0005] In the prior art, the announcement number CN117350068B discloses a method and device for numerical simulation of herbivorous fish in a complex water ecosystem model, which collects the basic data required for constructing a complex water ecosystem; based on the PCLake+ model, according to the collected basic data, a PCLake+Gras numerical simulation model of a complex water ecosystem is constructed; according to the collected physiological and biochemical process data of herbivorous fish, various parameters of herbivorous fish in the PCLake+Gras model are calibrated to obtain a calibrated PCLake+Gras model; the collected measured lake water quality and water ecology data are input into the calibrated PCLake+Gras model, and the water ecosystem simulation result of the simulation object is output. This scheme requires a large amount of accurate basic data, including water quality, water ecology, physiological and biochemical processes of herbivorous fish in complex water ecosystems, etc. However, in practical applications, data collection may not cover the spatial heterogeneity of the entire lake ecosystem (such as differences in water quality and fish distribution in different regions), or reflect the trend of long-term dynamic changes. At the same time, in complex aquatic ecosystems, many key ecological processes are highly nonlinear or random, but this scheme may not fully consider these factors, which limits the model's predictive ability. Therefore, the accuracy and effectiveness of the simulation results are reduced.

[0006] The above information disclosed in this Background section is only for enhancement of understanding of the background of the present disclosure and therefore it may contain information that does not constitute the prior art that is already known to one of ordinary skill in the art. Summary of the invention

[0007] The object of the present invention is to provide a method and system for numerical simulation of rapid gonadal maturation of Lateolabrax japonicus in a complex marine system model, so as to solve the problems raised in the above-mentioned background technology.

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

[0009] A numerical simulation method for the rapid maturation of gonads of Lateolabrax japonicus in a complex marine system model, the specific steps comprising:

[0010] Select several striped seabass of the same age of the species to be raised, randomly group the selected striped seabass groups of the same age to form several groups of test samples, set different salinity water samples for different groups of test samples, and keep the other environmental parameters of each group of test samples consistent except salinity, and raise the striped seabass of each group of test samples, where the raising time length is K;

[0011] After the feeding experiment is completed, the key hormone level of gonad development and the key gene expression ability level of each seabass in each group of test samples are randomly tested, wherein the key hormone level of gonad development is calculated by calculating the content of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone, and the key gene expression ability level is calculated by calculating the content of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor;

[0012] Based on the key hormone levels and key gene expression levels of gonadal development of each striped sea bass in the test sample, a neural network prediction model was established. The salinity of water samples from different groups of test samples was used as the model input, and the model was trained with the corresponding key hormone levels and key gene expression levels of gonadal development as labels to obtain a gonadal development prediction model.

[0013] The salinity of the sea area where the japonica sea bass is to be reared is collected, and the salinity is input into the trained gonad development prediction model to obtain the predicted values ​​of the key hormone levels and key gene expression ability levels of the japonica sea bass gonad development in this sea area, and the japonica sea bass gonad development index is generated by combining the environmental parameters of the sea area where the japonica sea bass is to be reared;

[0014] The disturbance parameters in the sea area where the striped seabass to be raised are obtained, and the striped seabass gonad development index is corrected based on the disturbance parameters to obtain the accurate striped seabass gonad development index. The simulation of the gonad development degree of the striped seabass is completed according to the accurate striped seabass gonad development index, wherein the disturbance parameters include the maximum wind speed on the sea surface and the sea circulation speed.

[0015] Furthermore, different salinity water samples are set for different groups of test samples, and other environmental parameters except salinity are kept consistent among the test samples of each group, wherein the other environmental parameters include water temperature, dissolved oxygen concentration, average daily light intensity and average daily nutrient content;

[0016] Each group of test samples of Leptotris japonicus was reared for a period of time K, wherein the period of time K was calibrated as [T 1 , T 1 +K], where T 1 Indicates the start of the feeding experiment, K indicates the length of the time period, the unit is day, and 75≤K≤180, K is a positive integer.

[0017] Furthermore, based on the key hormone levels and key gene expression levels of gonadal development of each sea bass in the test sample, a neural network prediction model was established. A neural network prediction model was established based on the LSTM model, and the activation function and optimization algorithm were selected. The tanh function was selected as the activation function, and Adam was selected as the optimization algorithm of the LSTM model. The formula of the Tanh function is:

[0018]

[0019] In the formula, f(r) represents the tanh function, and the independent variable r represents the weighted sum of the neuron's input, that is, the result of the weighted summation of the input received by the neuron from the previous layer;

[0020] At the same time, the hyperparameters of the LSTM model are set, and the hyperparameters of the LSTM model include: the number of network layers, the number of iterations, the learning rate, the batch size, the number of training times, the batch processing number, and the number of hidden layer neurons;

[0021] The number of network layers is set to a 3-layer network structure, the number of iterations is set to 200, the learning rate is set to 0.001, the batch size is set to 32, the number of training times is set to 100, the batch size is set to 256, and the number of hidden layer neurons is 32.

[0022] Furthermore, according to the predicted values ​​of the key hormone levels and key gene expression ability levels of the striped sea bass gonad development in this sea area, combined with the environmental parameters of the sea area where the striped sea bass is to be raised, the striped sea bass gonad development index is generated, wherein the specific formula for calculating the striped sea bass gonad development index is:

[0023]

[0024] Where, ZS represents the gonadal development index of Lateolabrax japonicus, G h and G g are the predicted values ​​of key hormone levels and key gene expression levels in the gonadal development of Leptospermum japonicum, T S To detect the sea water temperature, T 0 is the culture water temperature for the feeding experiment, L 0 is the average daily light intensity of the feeding experiment, L S is the average daily light intensity of the detected sea area, DO is the dissolved oxygen concentration of the detected sea area, and ya is the average daily nutrient content of the detected sea area;

[0025] Among them, the key hormone levels and key gene expression levels of L. japonicus gonadal development h and G g Characterized by the produced hormone content and related proteins, the specific calculation formula is:

[0026]

[0027] Wherein, GnRH is the predicted value of gonadotropin-releasing hormone content after the same time period as the feeding experiment, FSH is the predicted value of follicle-stimulating hormone content within the same time period as the feeding experiment, LH is the predicted value of luteinizing hormone content within the same time period as the feeding experiment, CYP, FSHR and LHR are the predicted values ​​of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor contents within the same time period as the feeding experiment, respectively.

[0028] Furthermore, the specific calculation formula for the average daily nutrient content ya of the detected sea area is:

[0029]

[0030] In the formula, and is the mean nitrate content and phosphate content of the detected sea area, SL is the number of biological species in the detected sea area, and Z pl is the abundance of zooplankton, α and β are the weight coefficients of the number of biological species in the detected sea area and the sum of the mean nitrate content and the mean phosphate content, respectively, where β>α, and α and β are both greater than 0, and the abundance of zooplankton Z pl The calculation is based on the formula:

[0031]

[0032] Where W total is the total mass of zooplankton in the sample water, V fil is the volume of water sample filtered with zooplankton, E net is the capture efficiency correction factor.

[0033] Furthermore, the gonad development index of the striped seabass was corrected based on the disturbance parameters to obtain the accurate gonad development index of the striped seabass, wherein the formula for calculating the accurate gonad development index of the striped seabass is:

[0034]

[0035] Where ZS′ is the precise gonadal development index of Lateolabrax japonicus, S wind is the maximum wind speed on the sea surface, HL is the sea circulation speed, ω 1 and ω 2 are the weight coefficients of the maximum wind speed on the sea surface and the sea circulation speed, respectively, where ω 2 ≥ω 1 And ω 1 and ω 2 All are greater than 0;

[0036] Set the judgment threshold of gonadal development of L. japonicus to yz, where the logic for judging the gonadal development of L. japonicus is:

[0037] When ZS′≥1.0*yz, it is judged that the gonad development of L. japonicus is excellent, which means that in this sea area, the gonad development speed of the L. japonicus species to be reared is fast and suitable for rearing;

[0038] When 0.4*yz≤ZS′<1.0*yz, it is judged that the gonad development of L. japonicus is good, which means that the gonad development speed of the L. japonicus species to be reared in this sea area is good;

[0039] When 0≤ZS′<0.4*yz, it is judged that the gonad development of the Lateolabrax japonicus is poor, indicating that in this sea area, the gonad development speed of the Lateolabrax japonicus species to be reared is slow and is not suitable for rearing.

[0040] The present invention also provides a numerical simulation system for a striped sea bass with rapid gonadal maturation in a complex sea system model, and the numerical simulation system for a striped sea bass with rapid gonadal maturation in a complex sea system model is used to execute the numerical simulation method for a striped sea bass with rapid gonadal maturation in a complex sea system model, comprising:

[0041] The test sample preparation module is used to select a number of striped sea bass of the same age of the species to be raised, randomly group the selected striped sea bass groups of the same age to form a number of test samples, set different salinity water samples for different groups of test samples, and keep the environmental parameters of each group of test samples consistent except salinity, and raise the striped sea bass of each group of test samples, where the raising time length is K;

[0042] The gene and hormone analysis module is used to randomly detect the key hormone level of gonadal development and the key gene expression ability level of each striped seabass in each group of test samples after the feeding experiment is completed, wherein the key hormone level of gonadal development is obtained by calculating the content of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone, and the key gene expression ability level is obtained by calculating the content of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor;

[0043] The prediction model training module is used to establish a neural network prediction model based on the key hormone levels and key gene expression levels of gonadal development of each striped sea bass in the test sample, take the salinity of water samples of different groups of test samples as model input, and use the corresponding key hormone levels and key gene expression levels of gonadal development as labels to train the model, and obtain a gonadal development prediction model;

[0044] The gonad development index calculation module is used to collect the salinity of the sea area where the japonica sea bass is to be raised, input the salinity into the trained gonad development prediction model, obtain the predicted values ​​of the key hormone levels and key gene expression ability levels of the japonica sea bass gonad development in the sea area, and generate the japonica sea bass gonad development index in combination with the environmental parameters of the sea area where the japonica sea bass is to be raised;

[0045] The precise development index simulation module is used to obtain the disturbance parameters in the sea area where the striped seabass to be raised grows, correct the striped seabass gonad development index based on the disturbance parameters to obtain the precise striped seabass gonad development index, and complete the simulation of the striped seabass gonad development degree according to the precise striped seabass gonad development index. The disturbance parameters include the maximum wind speed on the sea surface and the sea circulation speed.

[0046] Compared with the prior art, the present invention has the following beneficial effects:

[0047] First, through random grouping and multi-salinity experimental design, the effects of salinity on the key hormone levels and gene expression during the gonadal development of L. japonicus, including aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor, were systematically studied, and the multi-dimensional regulatory mechanism of salinity on gonadal development was fully revealed. Through deep learning and fitting of experimental data by neural network model, the complex nonlinear relationship between salinity and gonadal development indicators can be effectively captured, overcoming the limitation of traditional analysis methods that can only evaluate simple linear relationships. Secondly, the model can not only make predictions under different salinity conditions, but also accurately evaluate the levels of key hormones and gene expression capabilities of gonadal development when other environmental factors (such as temperature, nutrients, etc.) remain unchanged. Finally, the dynamic effects of complex marine environment on the gonadal development of L. japonicus were considered. In traditional studies, the effects of disturbance factors such as wind speed and circulation on reproductive physiology are often ignored, resulting in large deviations from the actual environment. This scheme uses disturbance parameters in natural sea areas such as maximum sea surface wind speed and ocean circulation speed as correction factors to make secondary corrections to the gonadal development index, which significantly improves the applicability and accuracy of model prediction. It provides technical reference for aquaculture and fishery management, and promotes the efficient and sustainable development of the sea bass industry. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] Figure 1 It is a schematic diagram of the overall method flow of the present invention;

[0049] Figure 2 It is a schematic diagram of the overall system structure of the present invention. DETAILED DESCRIPTION

[0050] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with specific embodiments.

[0051] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the present invention should be understood by people with ordinary skills in the field to which the present invention belongs. The words "first", "second" and similar words used in the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0052] Example:

[0053] See also Figure 1 , the present invention provides a technical solution:

[0054] A numerical simulation method for the rapid maturation of gonads of Lateolabrax japonicus in a complex marine system model, the specific steps comprising:

[0055] Step 1: Select several striped seabass of the same age of the species to be raised, randomly group the selected striped seabass groups of the same age to form several groups of test samples, set different salinity water samples for different groups of test samples, and keep other environmental parameters of each group of test samples consistent except salinity, and raise the striped seabass of each group of test samples, where the raising time is K.

[0056] Different salinity water samples were set for different groups of test samples. Except for salinity, other environmental parameters of each group of test samples were kept consistent, wherein the other environmental parameters included water temperature, dissolved oxygen concentration, average daily light intensity and average daily nutrient content;

[0057] Each group of test samples of Leptotris japonicus was reared for a period of time K, wherein the period of time K was calibrated as [T 1 , T 1 +K], where T 1 represents the start time of the feeding experiment, K represents the length of the time period, the unit is day, and 75≤K≤180, K is a positive integer.

[0058] The striped seabass of the same age are generally selected as one-year-old striped seabass. One-year-old striped seabass has important scientific significance and practical value as the research object. One-year-old striped seabass refers to individual striped seabass that grows for about one year after birth. The fish population at this stage is relatively stable in physiological and ecological characteristics, and it is a key stage for studying the dynamic changes of gonadal development of striped seabass. One-year-old striped seabass is usually in a stage of rapid growth, and its gonadal development begins to enter the initial stage, which is an important time window for studying the laws of reproductive development. From a physiological point of view, the reproductive system of fish at this stage is not yet fully mature, but the gonads have already shown certain developmental characteristics.

[0059] From an ecological perspective, one-year-old japonicus usually have good adaptability and respond significantly to changes in environmental factors (such as salinity, temperature, and light). Therefore, choosing one-year-old japonicus helps to clearly observe the impact of the external environment on its gonadal development.

[0060] Step 2: After the feeding experiment is completed, the key hormone levels for gonadal development and the key gene expression ability levels of each striped seabass in each group of test samples are randomly tested, wherein the key hormone levels for gonadal development are calculated by calculating the contents of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone, and the key gene expression ability levels are calculated by calculating the contents of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor.

[0061] The key hormone levels of gonadal development and the key gene expression ability levels of each striped seabass in each group of test samples were randomly tested. The key hormone levels of gonadal development were characterized by the produced hormone content, and the key gene expression ability levels were characterized by the content of related proteins of the parameters, including aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor.

[0062] In order to accurately determine the key hormone levels and key gene expression capacity of the gonadal development of sea bass, it is necessary to randomly select individuals from the experimental samples to collect tissue and blood samples. The detection method is enzyme-linked immunosorbent assay (ELISA), and the specific steps are: select a commercial ELISA kit for the content of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone to ensure specificity and sensitivity, dilute and process the serum sample according to the kit instructions, add it to a 96-well plate, and perform antibody binding, washing, color development and other steps in sequence, use an enzyme reader (450nm) to read the absorbance value, and calculate the hormone content according to the standard curve.

[0063] The expression level of key genes is characterized by the content of related proteins of the parameters, including aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor. The content of key proteins (aromatase, follicle-stimulating hormone receptor, luteinizing hormone receptor) reflects the final expression product level of the gene. Commonly used techniques are Western Blot or enzyme-linked immunosorbent assay (ELISA), and the specific steps include: extracting total protein from gonadal tissue, measuring protein concentration (such as BCA method), separating proteins by SDS-PAGE gel electrophoresis, transferring to PVDF membrane, detecting with specific antibodies against target proteins, developing with chemiluminescence (ECL), recording signal band intensity, and calculating relative protein levels by normalization with internal reference (such as β-actin).

[0064] Step 3: Based on the key hormone levels and key gene expression levels of gonadal development of each striped seabass in the experimental samples, a neural network prediction model was established. The salinity of water samples in different groups of experimental samples was used as the model input, and the model was trained with the corresponding key hormone levels and key gene expression levels of gonadal development as labels to obtain the gonadal development prediction model.

[0065] Based on the key hormone levels and key gene expression levels of gonadal development of each sea bass in the test sample, a neural network prediction model was established. A neural network prediction model was established based on the LSTM model, and the activation function and optimization algorithm were selected. The tanh function was selected as the activation function, and Adam was selected as the optimization algorithm of the LSTM model. The formula of the Tanh function is:

[0066]

[0067] In the formula, f(r) represents the tanh function, and the independent variable r represents the weighted sum of the neuron's input, that is, the result of the weighted summation of the input received by the neuron from the previous layer;

[0068] At the same time, the hyperparameters of the LSTM model are set, and the hyperparameters of the LSTM model include: the number of network layers, the number of iterations, the learning rate, the batch size, the number of training times, the batch processing number, and the number of hidden layer neurons;

[0069] The number of network layers is set to a 3-layer network structure, the number of iterations is set to 200, the learning rate is set to 0.001, the batch size is set to 32, the number of training times is set to 100, the batch size is set to 256, and the number of hidden layer neurons is 32.

[0070] LSTM can capture the dynamic regulation of hormone and gene expression during gonadal development. This time-based dynamic modeling capability can make up for the shortcomings of traditional statistical models (such as linear regression or multivariate regression) in dealing with time correlation. LSTM neural networks can fit highly complex nonlinear relationships through multiple layers of neurons and weight adjustment. The effect of salinity on key indicators of gonadal development may be regulated by multiple feedback mechanisms, and LSTM can capture these complex relationships more accurately.

[0071] Finally, the input is the salinity of the water in the growth environment of the striped seabass, and the output is the key hormone levels and key gene expression ability levels of the striped seabass' gonadal development within a specific time range.

[0072] Step 4: Collect the salinity of the sea area where the striped seabass to be raised grows, and input the salinity into the trained gonad development prediction model to obtain the predicted values ​​of the key hormone levels and key gene expression ability levels of the striped seabass gonad development in this sea area. Combined with the environmental parameters of the sea area where the striped seabass to be raised grows, generate the striped seabass gonad development index.

[0073] According to the predicted values ​​of the key hormone levels and key gene expression levels of the striped perch gonad development in this sea area, combined with the environmental parameters of the sea area where the striped perch will grow, the striped perch gonad development index is generated. The specific formula for calculating the striped perch gonad development index is:

[0074]

[0075] Where ZS represents the gonadal development index of Lateolabrax japonicus, G h and G g are the predicted values ​​of key hormone levels and key gene expression levels in the gonadal development of Leptospermum japonicum, T S To detect the sea water temperature, T 0 is the culture water temperature for the feeding experiment, L 0 is the average daily light intensity of the feeding experiment, L S is the average daily light intensity of the detected sea area, DO is the dissolved oxygen concentration of the detected sea area, and ya is the average daily nutrient content of the detected sea area;

[0076] In the formula, the gonad development index ZS of striped seabass characterizes the gonad development of striped seabass by integrating the predicted values ​​of key hormone levels and key gene expression ability levels of striped seabass gonad development and environmental factors. The larger the gonad development index ZS of striped seabass is, the higher the maturity of the gonad development of the striped seabass is and the faster the gonad development is.

[0077] G h This reflects the direct regulatory effect of the endocrine system of L. japonicus on gonadal development. By measuring the levels of gonadotropins (gonadotropin-releasing hormone, follicle-stimulating hormone, and luteinizing hormone), the hormone drive of gonadal development is revealed. g Characterizing the transcription and expression capacity of genes related to gonadal development (such as aromatase, gonadotropin receptor, etc.) is the molecular basis of gonadal development. These are key biological indicators reflecting the physiological process of gonadal development in striped seabass. The progress of gonadal development is closely related to the levels of sex hormones (such as gonadal hormones, sex hormone receptors, etc.), and the level of gene expression capacity can also affect gonadal development. The combination of the two is used to quantify the level of gonadal development in striped seabass. Hormone levels and gene expression capacity are the core driving forces of gonadal development, so G h and G g It is proportional to the gonad development index ZS of striped seabass. The two are introduced in the form of product, which reflects the importance of their synergistic effect. That is, only when the hormone level and gene expression reach a high level at the same time can they support good gonad development.

[0078] Water temperature is a key environmental factor in the development of the gonads of Lateolabrax japonicus, directly affecting its metabolic rate, enzyme activity, and hormone secretion. S -T 0 The larger the value of |, the less likely the Lateolabrax will adapt, and the gonadal development will be inhibited. S -T 0 |It is inversely proportional to the zS index of gonadal development of striped seabass. The water temperature difference is added to the denominator, which plays a penalty role, indicating that the greater the water temperature deviation, the lower the zS index of gonadal development of striped seabass.

[0079] Light affects the circadian rhythm of organisms, and then regulates the levels of hormones related to gonadal development (such as melatonin secreted by the pineal gland). If the actual light intensity is significantly different from the experimental conditions, it may cause biological rhythm disorders, thereby affecting gonadal development. 0 -L S |It is inversely proportional to the zS index of gonadal development of striped bass. The light difference is added to the denominator, which plays a penalty role similar to the water temperature difference.

[0080] Among them, the culture water temperature T 0 The average daily light intensity L 0 Generally, these are the optimal temperature and optimal daily light intensity for the growth and development of striped sea bass. The specific values ​​can be set by referring to relevant materials and combining expert experience.

[0081] DO represents the dissolved oxygen concentration in the tested sea area, which is an important environmental factor for the metabolism, energy consumption and gonad development of the striped seabass. The higher the dissolved oxygen level, the more sufficient the oxygen supply in the water body, and the more active the physiological functions of the striped seabass (such as oxidative metabolism and germ cell development), thereby promoting gonad development. Therefore, the dissolved oxygen concentration DO in the tested sea area is proportional to the gonad development index ZS of the striped seabass. The introduction of the logarithmic form ln(1+DO) weakens the impact of extremely high dissolved oxygen values, while ensuring that the negative impact of low dissolved oxygen concentrations is significantly amplified.

[0082] Nutrients are an important material basis for the development of the gonads of Lateolabrax japonicus. In particular, the formation of yolk protein and the growth of gonadal tissues have high energy and nutrient requirements. Nutritional status is crucial to the gonadal development of Lateolabrax japonicus. The average daily nutrient content directly affects the growth rate and gonadal development of fish. The reproductive health of fish is closely related to its nutritional status. Therefore, the average daily nutrient content ya in the tested sea area is proportional to the gonadal development index ZS of Lateolabrax japonicus. The index is expressed as The introduction of can highlight the positive effect of nutrient content while avoiding excessive amplification of extreme high values.

[0083] Among them, the key hormone levels and key gene expression levels of L. japonicus gonadal development h and G g Characterized by the produced hormone content and related proteins, the specific calculation formula is:

[0084]

[0085] Wherein, GnRH is the predicted value of gonadotropin-releasing hormone content after the same time period as the feeding experiment, FSH is the predicted value of follicle-stimulating hormone content within the same time period as the feeding experiment, LH is the predicted value of luteinizing hormone content within the same time period as the feeding experiment, CYP, FSHR and LHR are the predicted values ​​of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor contents within the same time period as the feeding experiment, respectively.

[0086] The specific calculation formula for the daily average nutrient content ya in the tested sea area is:

[0087]

[0088] In the formula, and is the mean nitrate content and phosphate content of the detected sea area, SL is the number of biological species in the detected sea area, and Z pl is the abundance of zooplankton, α and β are the weight coefficients of the number of biological species in the detected sea area and the sum of the mean nitrate content and the mean phosphate content, respectively, where β>α, and α and β are both greater than 0.

[0089] Among them, the daily average nutrient content ya is characterized by inorganic nutrients, zooplankton and species richness in the sea area. The larger the value, the richer the resources in the sea area, the more suitable it is for the growth of striped seabass, and the higher the degree of gonadal development of striped seabass.

[0090] Inorganic nutrients include and They are the main nutrient source for the growth of primary producer phytoplankton, and phytoplankton is the direct nutrient source for the bait organisms of the striped sea bass. Therefore, the nitrate and phosphate content in the tested sea area is proportional to the daily average nutrient content ya of the tested sea area.

[0091] The abundance of zooplankton Z pl The larger the value, the more zooplankton there are. Zooplankton is the direct food source of sea bass. Therefore, the abundance of zooplankton Z pl It is directly proportional to the average daily nutrient content ya in the tested sea area.

[0092] The number SL of biological species in the detected sea area can be obtained through field surveys, and its value represents the number of secondary producers. A high abundance of SL will consume the abundance of primary productivity phytoplankton and zooplankton, thereby reducing the nutrients available to the sea bass and producing a strong competitive effect. Therefore, the number SL of biological species in the detected sea area is inversely proportional to the average daily nutrient content ya of the detected sea area.

[0093] Among them, the higher the content of inorganic nutrients, the more primary producer phytoplankton it can produce and the faster the production rate, which has a greater impact on the nutrient supply of the detected sea area. Therefore, β>α is set, and both α and β are greater than 0.

[0094] The abundance of zooplankton is Z pl The calculation is based on the formula:

[0095]

[0096] Where W total is the total mass of zooplankton in the sample water, V fil is the volume of water sample filtered with zooplankton, E net is the capture efficiency correction factor. In the natural environment, the capture efficiency of zooplankton may vary depending on the mesh size and operation method. To improve the calculation accuracy, the capture efficiency correction factor is added, where the capture efficiency correction factor E net Generally it is 0.6-0.9, depending on the type of net and experimental conditions.

[0097] Step 5: Obtain disturbance parameters in the sea area where the striped seabass to be raised grows, correct the striped seabass gonad development index based on the disturbance parameters to obtain an accurate striped seabass gonad development index, and complete the simulation of the striped seabass gonad development degree according to the accurate striped seabass gonad development index, wherein the disturbance parameters include the maximum sea surface wind speed and the sea circulation speed.

[0098] The gonad development index of striped sea bass was corrected based on the disturbance parameters to obtain the accurate gonad development index of striped sea bass. The formula for calculating the accurate gonad development index of striped sea bass is:

[0099]

[0100] Where ZS′ is the precise gonadal development index of Lateolabrax japonicus, S wind is the maximum wind speed on the sea surface, HL is the sea circulation speed, ω 1 and ω 2 are the weight coefficients of the maximum wind speed on the sea surface and the sea circulation speed, respectively, where ω 2 ≥ω 1 And ω 1 and ω 2 Both are greater than 0.

[0101] Among them, the precise gonad development index ZS′ of striped seabass is a disturbance parameter that integrates the complex marine environment and describes the precise value of the gonad development index of striped seabass under the disturbance parameters. The larger the value, the higher the maturity of the gonad development of the striped seabass and the faster the gonad development.

[0102] The maximum wind speed on the sea surface is an important parameter affecting ocean disturbances. High wind speeds can cause water mixing, affect the vertical transport of nutrients, and affect the mixing degree of the water area and the distribution of nutrients, which in turn affects the growth and gonadal development of L. japonicus. Therefore, the maximum wind speed on the sea surface is inversely proportional to the precise gonadal development index of L. japonicus. The exponential function shows that the effect on gonadal development increases significantly when the wind speed is too high.

[0103] The sea circulation velocity HL can significantly affect the stability of local habitats. Excessive circulation velocity may lead to the dispersion of bait organisms and habitat destruction, thereby inhibiting gonadal development. Therefore, the sea circulation velocity HL is inversely proportional to the precise gonadal development index of Lateolabrax japonicus. The logarithmic function shows that when the circulation velocity is too high, the impact on gonadal development increases significantly.

[0104] The maximum sea surface wind speed refers to the maximum value of the wind speed on the sea surface within a specified time. It is usually used to assess the environmental dynamic conditions and potential meteorological impacts of the sea area. The specific method of obtaining it is: deploy anemometers (such as ultrasonic anemometers or mechanical anemometers) in specific sea areas to directly measure the sea surface wind speed; the time series of wind speed needs to be recorded during measurement, and the maximum wind speed value is obtained through statistical analysis.

[0105] The ocean circulation velocity can be measured by deploying ADCP equipment in the ocean and using the acoustic Doppler effect to measure the flow rate and direction of different water layers.

[0106] Since the ocean circulation speed can directly change the distribution of bait organisms and nutrients, setting ω 2 ≥ω 1 And ω 1 and ω 2 Both are greater than 0.

[0107] Set the judgment threshold of gonadal development of L. japonicus to yz, where the logic for judging the gonadal development of L. japonicus is:

[0108] When ZS′≥1.0*yz, it is judged that the gonad development of L. japonicus is excellent, which means that in this sea area, the gonad development speed of the L. japonicus species to be reared is fast and suitable for rearing;

[0109] When 0.4*yz≤ZS′<1.0*yz, it is judged that the gonad development of L. japonicus is good, which means that the gonad development speed of the L. japonicus species to be reared in this sea area is good;

[0110] When 0≤ZS′<0.4*yz, it is judged that the gonadal development of the japonica seabass is poor, indicating that the gonadal development of the japonica seabass species to be reared in this sea area is slow and unsuitable for rearing. The judgment threshold of the gonadal development of the japonica seabass is set by consulting relevant materials and combining expert experience.

[0111] See also Figure 2The present invention also provides a numerical simulation system for the rapid maturation of gonads of striped sea bass in a complex sea system model, and the numerical simulation system for the rapid maturation of gonads of striped sea bass in a complex sea system model is used to execute the numerical simulation method for the rapid maturation of gonads of striped sea bass in a complex sea system model, comprising:

[0112] The test sample preparation module is used to select a number of striped sea bass of the same age of the species to be raised, randomly group the selected striped sea bass groups of the same age to form a number of test samples, set different salinity water samples for different groups of test samples, and keep the environmental parameters of each group of test samples consistent except salinity, and raise the striped sea bass of each group of test samples, where the raising time length is K;

[0113] The gene and hormone analysis module is used to randomly detect the key hormone level of gonadal development and the key gene expression ability level of each striped seabass in each group of test samples after the feeding experiment is completed, wherein the key hormone level of gonadal development is obtained by calculating the content of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone, and the key gene expression ability level is obtained by calculating the content of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor;

[0114] The prediction model training module is used to establish a neural network prediction model based on the key hormone levels and key gene expression levels of gonadal development of each striped sea bass in the test sample, take the salinity of water samples of different groups of test samples as model input, and use the corresponding key hormone levels and key gene expression levels of gonadal development as labels to train the model, and obtain a gonadal development prediction model;

[0115] The gonad development index calculation module is used to collect the salinity of the sea area where the japonica sea bass is to be raised, input the salinity into the trained gonad development prediction model, obtain the predicted values ​​of the key hormone levels and key gene expression ability levels of the japonica sea bass gonad development in the sea area, and generate the japonica sea bass gonad development index in combination with the environmental parameters of the sea area where the japonica sea bass is to be raised;

[0116] The precise development index simulation module is used to obtain the disturbance parameters in the sea area where the striped seabass to be raised grows, correct the striped seabass gonad development index based on the disturbance parameters to obtain the precise striped seabass gonad development index, and complete the simulation of the striped seabass gonad development degree according to the precise striped seabass gonad development index. The disturbance parameters include the maximum wind speed on the sea surface and the sea circulation speed.

[0117] The above embodiments may be implemented in whole or in part by software, hardware, firmware or any other combination thereof. When implemented by software, the above embodiments may be implemented in whole or in part in the form of a computer program product. Those skilled in the art may appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein may be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed by hardware or software methods depends on the specific application and design constraints of the technical solution.

[0118] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, and may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0119] The above description is only a specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application.

Claims

1. A numerical simulation method for the rapid maturation of gonads of Lateolabrax japonicus in a complex marine system model, characterized in that: The specific steps include: Select several striped seabass of the same age of the species to be raised, randomly group the selected striped seabass groups of the same age to form several groups of test samples, set different salinity water samples for different groups of test samples, and keep the other environmental parameters of each group of test samples consistent except salinity, and raise the striped seabass of each group of test samples, where the raising time length is K; After the feeding experiment is completed, the key hormone level of gonad development and the key gene expression ability level of each seabass in each group of test samples are randomly tested, wherein the key hormone level of gonad development is calculated by calculating the content of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone, and the key gene expression ability level is calculated by calculating the content of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor; Based on the key hormone levels and key gene expression levels of gonadal development of each striped sea bass in the test sample, a neural network prediction model was established. The salinity of water samples from different groups of test samples was used as the model input, and the model was trained with the corresponding key hormone levels and key gene expression levels of gonadal development as labels to obtain a gonadal development prediction model. The salinity of the sea area where the japonica sea bass is to be reared is collected, and the salinity is input into the trained gonad development prediction model to obtain the predicted values ​​of the key hormone levels and key gene expression ability levels of the japonica sea bass gonad development in this sea area, and the japonica sea bass gonad development index is generated by combining the environmental parameters of the sea area where the japonica sea bass is to be reared; The disturbance parameters in the sea area where the striped seabass to be raised are obtained, and the striped seabass gonad development index is corrected based on the disturbance parameters to obtain the accurate striped seabass gonad development index. The simulation of the gonad development degree of the striped seabass is completed according to the accurate striped seabass gonad development index, wherein the disturbance parameters include the maximum wind speed on the sea surface and the sea circulation speed.

2. The method for numerically simulating the rapid gonadal maturation of Lateolabrax japonicus in a complex marine system model according to claim 1, characterized in that: Different salinity water samples were set for different groups of test samples. Except for salinity, other environmental parameters of each group of test samples were kept consistent, wherein the other environmental parameters included water temperature, dissolved oxygen concentration, average daily light intensity and average daily nutrient content; Each group of test sample sea bass was raised, wherein the raising time length was K, wherein the raising time length K was calibrated as [T1, T1+K], wherein T1 represents the start time of the raising experiment, K represents the length of the time period, the unit is day, and 75≤K≤180, K is a positive integer.

3. The method for numerically simulating the rapid gonadal maturation of Lateolabrax japonicus in a complex marine system model according to claim 2, characterized in that: Based on the key hormone levels and key gene expression levels of gonadal development of each sea bass in the test sample, a neural network prediction model was established. A neural network prediction model was established based on the LSTM model, and the activation function and optimization algorithm were selected. The tanh function was selected as the activation function, and Adam was selected as the optimization algorithm of the LSTM model. The formula of the Tanh function is: In the formula, f(r) represents the tanh function, and the independent variable r represents the weighted sum of the neuron's input, that is, the result of the weighted summation of the input received by the neuron from the previous layer; At the same time, the hyperparameters of the LSTM model are set, and the hyperparameters of the LSTM model include: the number of network layers, the number of iterations, the learning rate, the batch size, the number of training times, the batch processing number, and the number of hidden layer neurons; The number of network layers is set to a 3-layer network structure, the number of iterations is set to 200, the learning rate is set to 0.001, the batch size is set to 32, the number of training times is set to 100, the batch size is set to 256, and the number of hidden layer neurons is 32.

4. The method for numerically simulating the rapid gonadal maturation of Lateolabrax japonicus in a complex marine system model according to claim 3, characterized in that: According to the predicted values ​​of the key hormone levels and key gene expression levels of the striped perch gonad development in this sea area, combined with the environmental parameters of the sea area where the striped perch will grow, the striped perch gonad development index is generated. The specific formula for calculating the striped perch gonad development index is: Where ZS represents the gonadal development index of Lateolabrax japonicus, G h and G g are the predicted values ​​of key hormone levels and key gene expression levels in the gonadal development of Leptospermum japonicum, T S To detect the sea water temperature, T0 is the culture water temperature of the feeding experiment, L0 is the average daily light intensity of the feeding experiment, L S is the average daily light intensity of the detected sea area, DO is the dissolved oxygen concentration of the detected sea area, and ya is the average daily nutrient content of the detected sea area; Among them, the key hormone levels and key gene expression levels of L. japonicus gonadal development h and G g Characterized by the produced hormone content and related proteins, the specific calculation formula is: Wherein, GnRH is the predicted value of gonadotropin-releasing hormone content after the same time period as the feeding experiment, FSH is the predicted value of follicle-stimulating hormone content within the same time period as the feeding experiment, LH is the predicted value of luteinizing hormone content within the same time period as the feeding experiment, CYP, FSHR and LHR are the predicted values ​​of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor contents within the same time period as the feeding experiment, respectively.

5. The method for numerically simulating the rapid gonadal maturation of Lateolabrax japonicus in a complex marine system model according to claim 4, characterized in that: The specific calculation formula for the daily average nutrient content ya in the tested sea area is: In the formula, and is the mean nitrate content and phosphate content of the detected sea area, SL is the number of biological species in the detected sea area, and Z pl is the abundance of zooplankton, α and β are the weight coefficients of the number of biological species in the detected sea area and the sum of the mean nitrate content and the mean phosphate content, respectively, where β>α, and α and β are both greater than 0, and the abundance of zooplankton Z pl The calculation is based on the formula: Where W total is the total mass of zooplankton in the sample water, V fil is the volume of water sample filtered with zooplankton, E net is the capture efficiency correction factor.

6. The method for numerically simulating the rapid maturation of gonads of Lateolabrax japonicus in a complex marine system model according to claim 4, characterized in that: The gonad development index of striped sea bass was corrected based on the disturbance parameters to obtain the accurate gonad development index of striped sea bass. The formula for calculating the accurate gonad development index of striped sea bass is: Where ZS′ is the precise gonadal development index of Lateolabrax japonicus, S wind is the maximum wind speed on the sea surface, HL is the sea circulation speed, ω1 and ω2 are the weight coefficients of the maximum wind speed on the sea surface and the sea circulation speed, respectively, where ω2≥ω1 and ω1 and ω2 are both greater than 0; Set the judgment threshold of gonadal development of L. japonicus to yz, where the logic for judging the gonadal development of L. japonicus is: When ZS′≥1.0*yz, it is judged that the gonad development of L. japonicus is excellent, which means that in this sea area, the gonad development speed of the L. japonicus species to be reared is fast and suitable for rearing; When 0.4*yz≤ZS′<1.0*yz, it is judged that the gonad development of L. japonicus is good, which means that the gonad development speed of the L. japonicus species to be reared in this sea area is good; When 0≤ZS′<0.4*yz, it is judged that the gonad development of the Lateolabrax japonicus is poor, indicating that in this sea area, the gonad development speed of the Lateolabrax japonicus species to be reared is slow and is not suitable for rearing.

7. A numerical simulation system for the rapid maturation of gonads of Lateolabrax japonicus in a complex marine system model, characterized by: The numerical simulation system for the rapid maturation of gonads of Lateolabrax japonicus in a complex sea system model is used to execute the numerical simulation method for the rapid maturation of gonads of Lateolabrax japonicus in a complex sea system model according to any one of claims 1 to 6, comprising: The test sample preparation module is used to select a number of striped sea bass of the same age of the species to be raised, randomly group the selected striped sea bass groups of the same age to form a number of test samples, set different salinity water samples for different groups of test samples, and keep the environmental parameters of each group of test samples consistent except salinity, and raise the striped sea bass of each group of test samples, where the raising time length is K; The gene and hormone analysis module is used to randomly detect the key hormone level of gonadal development and the key gene expression ability level of each striped seabass in each group of test samples after the feeding experiment is completed, wherein the key hormone level of gonadal development is obtained by calculating the content of gonadotropin-releasing hormone, follicle-stimulating hormone and luteinizing hormone, and the key gene expression ability level is obtained by calculating the content of aromatase, follicle-stimulating hormone receptor and luteinizing hormone receptor; The prediction model training module is used to establish a neural network prediction model based on the key hormone levels and key gene expression levels of gonadal development of each striped sea bass in the test sample, take the salinity of water samples of different groups of test samples as model input, and use the corresponding key hormone levels and key gene expression levels of gonadal development as labels to train the model, and obtain a gonadal development prediction model; The gonad development index calculation module is used to collect the salinity of the sea area where the japonica sea bass is to be raised, input the salinity into the trained gonad development prediction model, obtain the predicted values ​​of the key hormone levels and key gene expression ability levels of the japonica sea bass gonad development in the sea area, and generate the japonica sea bass gonad development index in combination with the environmental parameters of the sea area where the japonica sea bass is to be raised; The precise development index simulation module is used to obtain the disturbance parameters in the sea area where the striped seabass to be raised grows, correct the striped seabass gonad development index based on the disturbance parameters to obtain the precise striped seabass gonad development index, and complete the simulation of the striped seabass gonad development degree according to the precise striped seabass gonad development index. The disturbance parameters include the maximum wind speed on the sea surface and the sea circulation speed.

Images