CN114547992A - Method and device for constructing tidal estuary brine invasion model and electronic equipment - Google Patents
Method and device for constructing tidal estuary brine invasion model and electronic equipment Download PDFInfo
- Publication number
- CN114547992A CN114547992A CN202210087093.0A CN202210087093A CN114547992A CN 114547992 A CN114547992 A CN 114547992A CN 202210087093 A CN202210087093 A CN 202210087093A CN 114547992 A CN114547992 A CN 114547992A
- Authority
- CN
- China
- Prior art keywords
- inflow
- tidal
- sea
- invasion
- brine
- 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.)
- Pending
Links
- 230000009545 invasion Effects 0.000 title claims abstract description 147
- 239000012267 brine Substances 0.000 title claims abstract description 82
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 50
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims abstract description 73
- 239000011780 sodium chloride Substances 0.000 claims abstract description 73
- 230000035945 sensitivity Effects 0.000 claims abstract description 60
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 51
- 150000003839 salts Chemical class 0.000 claims abstract description 20
- 230000015654 memory Effects 0.000 claims description 21
- 238000003860 storage Methods 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 4
- 238000004519 manufacturing process Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000007726 management method Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 2
- 239000013505 freshwater Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02B—HYDRAULIC ENGINEERING
- E02B1/00—Equipment or apparatus for, or methods of, general hydraulic engineering, e.g. protection of constructions against ice-strains
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Human Resources & Organizations (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Tourism & Hospitality (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Public Health (AREA)
- Life Sciences & Earth Sciences (AREA)
- Structural Engineering (AREA)
- Evolutionary Biology (AREA)
- Civil Engineering (AREA)
- Mechanical Engineering (AREA)
- Primary Health Care (AREA)
- Probability & Statistics with Applications (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Quality & Reliability (AREA)
- Computing Systems (AREA)
Abstract
The invention discloses a method for constructing a tidal estuary saline water invasion model, which comprises the following steps: acquiring data information of a target river mouth, wherein the data information comprises tidal ranges and maximum salt water invasion areas corresponding to turning points of different sea inflow rates; determining a first sensitivity coefficient and a first relation between the sea inflow and the tidal range according to the sea inflow, the tidal range and the maximum saline invasion area; determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area; and constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point. The invention can predict the entry state of the saline water at the river mouth and estimate the critical flow of the saline water entry by determining the tidal river mouth saline water entry model, can make quick judgment on the saline water entry in the practical production application, and provides data support for water supply safety and water resource management.
Description
Technical Field
The invention relates to the technical field of water power, in particular to a method and a device for constructing a tidal estuary brine invasion model and electronic equipment.
Background
The estuary region is the intersection of the ocean and the river and is one of the most economically active and densely populated regions in the world. Along with the rapid development of social economy, the resource development and utilization of the estuary region are deepened more and more, and the problem of the invasion of estuary brine comes along with the development of the estuary region, so that the fresh water supply of the residential life, the farmland irrigation, the urban industrial production and the like of the estuary region is directly influenced. Estuary brine invasion is an important and common problem faced by human activities, and as a global problem, the brine invasion is aggravated by the increase of fresh water demand in coastal areas, the rise of sea level, climate change and the like. Therefore, the mechanism and the rule of the brine invasion are determined to have important practical significance.
The saline invasion related research and technology is mainly based on field observation, numerical simulation, mathematical analysis and the like, the field observation usually needs a large amount of manpower and material resources, the numerical simulation needs a large amount of basic data, the parameter calibration and verification also needs a large amount of time, and a common mathematical model method is lack of a physical mechanism; resulting in long time consumption and low efficiency in the whole research process.
Disclosure of Invention
In view of this, the embodiment of the invention provides a method for constructing a tidal estuary brine invasion model, so as to solve the problems of long time consumption and low efficiency in a brine invasion research process in the prior art.
In order to achieve the purpose, the invention provides the following technical scheme:
the embodiment of the invention provides a method for constructing a tidal estuary saline water invasion model, which comprises the following steps:
acquiring data information of a target river mouth, wherein the data information comprises tidal ranges and maximum salt water invasion areas corresponding to turning points of different sea inflow rates;
determining a first sensitivity coefficient and a first relation between the sea inflow and the tidal difference according to the sea inflow, the tidal difference and the maximum saline invasion area;
determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area;
and constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point.
Optionally, the determining a first sensitivity coefficient and a first relationship between the sea inflow and the tidal range according to the sea inflow, the tidal range and the maximum saline invasion area includes:
obtaining a first sensitivity coefficient between the tidal difference and the saline invasion area by carrying out calibration fitting on the maximum saline invasion area and the tidal difference;
and fitting the sea inflow rate and the tidal range to obtain a first relation between different sea inflow rates and tidal ranges.
Optionally, the obtaining a first relationship between different sea inflow rates and tidal ranges by fitting the sea inflow rates and the tidal ranges includes:
calculating the average value of the sea inflow;
comparing the different sea inflow rates with the average value to obtain a plurality of dimensionless sea inflow rate data;
fitting the dimensionless inflow volume data and the corresponding tidal range to obtain a first relation between different inflow volumes and tidal ranges.
Optionally, the determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area includes:
obtaining a middle turning point of the inflow flow by performing calibration fitting on the inflow flow and the maximum saline water invasion area;
and calculating to obtain a second sensitivity coefficient between the sea inflow and the saline invasion area according to the sea inflow, the maximum saline invasion area and the middle turning point.
Optionally, the obtaining of the middle turning point of the sea inflow by performing calibration fitting on the sea inflow and the maximum saline invasion area includes:
calculating the average value of the sea inflow;
comparing the different sea inflow rates with the average value to obtain a plurality of dimensionless sea inflow rate data;
carrying out calibration fitting on the maximum saline water invasion area and the dimensionless inflow volume data to obtain a fitting curve;
and determining the abscissa value corresponding to the turning point in the fitting curve as the middle turning point of the sea inflow.
Optionally, the constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the intermediate turning point includes:
constructing a relationship between saline invasion area and tidal range based on the maximum saline invasion area, the first sensitivity coefficient and the first relationship;
constructing a relation between the saline invasion area and the inflow flow based on the maximum saline invasion area, the second sensitivity coefficient and the middle turning point;
and obtaining a tidal estuary brine invasion model according to the relationship among the maximum brine invasion area, the brine invasion area and the tidal range and the relationship among the brine invasion area and the inflow flow.
Optionally, the method further includes:
acquiring sea inflow data of a target estuary in real time;
and substituting the inflow volume data into the tidal estuary brine invasion model to obtain a brine invasion area.
The embodiment of the invention also provides a tidal estuary saline invasion model device, which comprises:
the acquisition module is used for acquiring data information of a target estuary, wherein the data information comprises tidal ranges and maximum salt water invasion areas corresponding to turning points of different sea inflow rates;
a first analysis module for determining a first sensitivity coefficient and a first relationship between the sea inflow and the tidal range according to the sea inflow, the tidal range and the maximum saline invasion area;
the second analysis module is used for determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area;
and the construction module is used for constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point.
An embodiment of the present invention further provides an electronic device, including:
the device comprises a memory and a processor, wherein the memory and the processor are mutually connected in a communication mode, the memory stores computer instructions, and the processor executes the computer instructions so as to execute the method for constructing the tidal estuary brine invasion model provided by the embodiment of the invention.
Embodiments of the present invention also provide a computer-readable storage medium storing computer instructions for causing a computer to execute the method for constructing the tidal estuary brine intrusion model provided by the embodiments of the present invention.
The technical scheme of the invention has the following advantages:
the invention provides a method for constructing a tidal estuary brine invasion model, which comprises the steps of acquiring data information of a target estuary, wherein the data information comprises tidal ranges and maximum brine invasion areas corresponding to turning points with different sea inflow rates; determining a first sensitivity coefficient and a first relation between the sea inflow and the tidal range according to the sea inflow, the tidal range and the maximum saline invasion area; determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area; and constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point. The invention can predict the entry state of the saline water at the river mouth and estimate the critical flow of the saline water entry by determining the tidal river mouth saline water entry model, can make quick judgment on the saline water entry in the practical production application, and provides data support for water supply safety and water resource management.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a flow chart of a method of constructing a tidal estuary brine intrusion model in an embodiment of the invention;
FIG. 2 is a flow chart of determining a first sensitivity factor and a first relationship between sea inflow and tidal range according to an embodiment of the present invention;
FIG. 3 is a flow chart of determining a first relationship in an embodiment in accordance with the invention;
FIG. 4 is a flow chart of determining a second sensitivity factor and an intermediate inflection point in accordance with an embodiment of the present invention;
FIG. 5 is a flow chart of determining an intermediate turning point according to an embodiment of the present invention;
FIG. 6 is a flow chart of a method for constructing a tidal estuary saltwater intrusion model in accordance with an embodiment of the present invention;
FIG. 7 is a flow chart of an embodiment of the invention using a tidal estuary saltwater intrusion model;
FIG. 8 is a schematic structural diagram of a device for constructing a tidal estuary brine invasion model according to an embodiment of the invention;
FIG. 9 is a schematic structural diagram of an electronic device in an embodiment of the invention;
FIG. 10 is a schematic view of an S-shaped curve in an embodiment of the present invention;
FIG. 11 is a schematic view of a curve fitting tidal range to saline invasion area at different flow rates in an embodiment of the present invention;
FIG. 12 shows the dimensionless sea inflow and tDA schematic of the relationship;
FIG. 13 is a graphical illustration of a plot of dimensionless sea inflow rate as a function of maximum salt water intrusion area overproof rate in an embodiment of the present invention;
FIG. 14 is a diagram illustrating fitting verification according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In accordance with an embodiment of the present invention, there is provided an embodiment of a method for constructing a tidal estuary brine intrusion model, wherein the steps illustrated in the flow chart of the drawings may be performed in a computer system, such as a set of computer executable instructions, and wherein the logical sequence is illustrated in the flow chart, in some cases the steps illustrated or described may be performed in a different order than presented herein.
In this embodiment, a method for constructing a tidal estuary brine intrusion model is provided, which can be used in an occasion where a estuary region needs to predict brine intrusion, as shown in fig. 1, the method for constructing the tidal estuary brine intrusion model includes the following steps:
step S1: and acquiring data information of the target river mouth, wherein the data information comprises tidal ranges and maximum salt water invasion areas corresponding to turning points of different sea inflow rates. Specifically, the tidal range and the maximum salt water invasion area corresponding to the turning points with different sea inflow rates are obtained by performing statistical calculation based on the estuary salinity data, the sea inflow rate and the tidal change.
Step S2: a first sensitivity coefficient and a first relationship between the sea inflow and the tidal range are determined based on the sea inflow, the tidal range, and the maximum saltwater intrusion area.
Step S3: and determining a second sensitivity coefficient and an intermediate turning point according to the inflow of the sea and the maximum saline invasion area.
Step S4: and constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point.
Through the steps S1 to S4, the method for constructing the tidal estuary brine invasion model according to the embodiment of the present invention can predict the estuary brine invasion state and estimate the critical flow of brine invasion by determining the tidal estuary brine invasion model, and can make a quick judgment for the brine invasion in actual production application, so as to provide data support for water supply safety and water resource management.
Specifically, in an embodiment, as shown in fig. 2, the step S2 includes the following steps:
step S21: and obtaining a first sensitivity coefficient between the tidal difference and the saline invasion area by carrying out calibration fitting on the maximum saline invasion area and the tidal difference.
Step S22: and fitting the sea inflow rate and the tidal range to obtain a first relation between different sea inflow rates and tidal ranges. Specifically, the first relationship refers to a linear relationship between the tidal range corresponding to the turning point and the dimensionless inflow sea current.
Specifically, the general form of the tidal estuary salt water invasion diffusion curve is as follows: salt water intrusion is simultaneously influenced by the effects of tides as driving force and ocean currents as resistance. The saline invasion intensity changes in a monotonically increasing S-shaped change of "slow-fast-slow". A, D, B are the three turning points of the S-shaped curve, respectively, as shown in FIG. 10.
According to the basic form of the saline invasion diffusion curve, selecting a curve tangent function as a basic function, and constructing a conceptual model through translation and generalization treatment, wherein the conceptual model comprises the following steps:
wherein: x represents the driving or impedance factor during the salt water intrusion, alpha represents the sensitivity coefficient, S represents the maximum limit value during the seawater intrusion (selected seawater intrusion area), xDAnd represents the corresponding numerical value of the turning point in the curve.
According to the central symmetry of the S-curve, by performing a third order derivation on equation 1, the critical salinity and the sensitivity coefficient corresponding to the A, D, B three turning points can be calculated as follows:
based on equation (1), the relationship between the saline invasion area and tidal range is constructed:
during the construction process, a first sensitivity factor alpha is required1First relation t between sum tidal range and inflowDObtaining tidal range and maximum saline water invasion area under different flow rates by obtaining data information of a target estuary, substituting the data information into a formula (4), carrying out rating fitting on parameters, and determining alpha through multiple rating fitting1The value of (a) is shown in figure 11, and when there is a difference in the inflow, the function curve fitting is detailed, and the linear fitting is performed on the tidal range and the inflow, so that t can be obtainedD。
For example: TABLE 1 shows the tidal ranges (t) corresponding to the turning points at different sea inflow ratesD) And the maximum saline invasion area (A)imax)。
TABLE 1
According to Table 1, the tidal difference (t) is associated with the turning pointD) And dimensionless inflow to the sea
A linear fit is performed, as in fig. 12, determining the linear relationship between:
coefficient of sensitivity alpha 12 and linear relation
Substituting equation (4) yields the area of saline invasion as a function of tidal range as follows:
specifically, in an embodiment, as shown in fig. 3, the step S22 includes the following steps:
step S221: the average value of the sea inflow was calculated.
Step S222: and comparing the different sea inflow rates with the average value to obtain a plurality of dimensionless sea inflow rate data.
Step S223: fitting the dimensionless inflow data and the corresponding tidal range to obtain a first relation between different inflow and tidal range.
Specifically, in an embodiment, as shown in fig. 4, the step S3 includes the following steps:
step S31: and obtaining the middle turning point of the sea inflow by performing calibration fitting on the sea inflow and the maximum saline water invasion area.
Step S32: and calculating to obtain a second sensitivity coefficient between the sea inflow and the saline invasion area according to the sea inflow, the maximum saline invasion area and the middle turning point.
Specifically, based on formula 1, a relationship between the maximum area of the brine invasion and the inflow flow is constructed:
Amax＝max(A(xi)) (7)
during the construction process, a second coefficient of sensitivity is requiredα2And intermediate turning point
Tidal ranges and maximum salt water invasion areas under different flow rates can be obtained by obtaining data information of a target river mouth, the data information is substituted into the formula (5), and the parameters are subjected to calibration fitting, such as the data information shown in fig. 13, so that the tidal ranges and the maximum salt water invasion areas under different flow rates can be obtained
Then alpha is determined according to the sea inflow, the maximum saline water invasion area and the middle turning point2The value of (a) is,
for example, based on the formula (6) and the formula (7), the data in the above table 1 are substituted into the calculation to obtain
α2＝5.0,Amax1. The relationship between the maximum saline invasion area and the inflow is therefore as follows:
wherein x isi: the sea inflow rate;
the middle turning point (the middle turning point of the dimensionless sea inflow corresponding to the relation between the maximum saline water invasion area and the sea inflow).
Specifically, in an embodiment, as shown in fig. 5, the step S31 includes the following steps:
step S311: the average value of the inflow to the sea is calculated.
Step S312: and comparing the different sea inflow rates with the average value to obtain a plurality of dimensionless sea inflow rate data.
Step S313: and carrying out calibration fitting on the maximum saline water invasion area and the dimensionless inflow volume data to obtain a fitting curve.
Step S314: will be fitted withAnd determining the abscissa value corresponding to the turning point in the curve as the middle turning point of the sea inflow. Specifically, the abscissa value corresponding to the turning point D in the fitting curve is XD。
Specifically, in an embodiment, as shown in fig. 6, the step S4 includes the following steps:
step S41: and constructing a relation between the saline invasion area and the tidal range based on the maximum saline invasion area, the first sensitivity coefficient and the first relation.
Step S42: and obtaining a tidal estuary brine invasion model according to the maximum brine invasion area, the relation between the brine invasion area and the tidal range and the relation between the brine invasion area and the inflow flow.
Specifically, the formula (4), the formula (6) and the formula (7) are integrated, processed and calculated, so that a tidal estuary brine invasion model can be obtained:
wherein x isi: the sea inflow rate; a. themax: area of maximum saline invasion (in x)iMaximum area reached by saline invasion under conditions); t: tidal range; alpha is alpha1: first coefficient of sensitivity, α2: a second coefficient of sensitivity; t is tD: a first relationship;
dimensionless sea inflow (measured/average);
the function of the maximum salt water invasion area and the sea inflow rate corresponds to the middle turning point of the dimensionless sea inflow rate.
According to the above formula and related parameters, in practical application, further according to the simulation data (or measured value), the calculation determination A is performedmax、tD、
α1、α2Substituting the parameter values into the tidal estuary brine invasion model to obtain the final tidal estuary brine invasion model.
For example, according to the above process, the parameters and the relations that can be determined are: alpha is alpha1＝2.0,
α2＝5.0、A max1 and
therefore, a final tidal estuary salt water intrusion model can be obtained, as shown in equation 10.
By the above calculation, under different sea inflow rates, the salinity exceeding rate is always increased in a 'slow-fast-slow' S shape along with the increase of tidal range, and then is gradually reduced when the maximum area of salt water invasion is reached. According to the formula (10), the flow rate of the incoming sea is less than 10000m3In the time of/s, the salt water invasion exceeding standard gradually approaches to the peak value, and the change process tends to be gentle; when the flow is more than 10000m3At the time of/s, the invasion intensity of the saline water is continuously weakened; the flow rate reaches 25000m3And when the water content is more than the second, the salt water invasion intensity is very small, and the region basically has no salt water invasion trouble and accords with the model simulation result.
Specifically, in an embodiment, as shown in fig. 7, the method further includes the following steps:
step S51: and acquiring the sea inflow data of the target estuary in real time.
Step S52: and substituting the inflow volume data into a tidal estuary brine invasion model to obtain a brine invasion area.
Specifically, by determining a tidal estuary brine invasion model, the estuary brine invasion state can be predicted according to real-time inflow volume data, the critical flow of brine invasion is estimated, and in practical production application, the method can make a quick judgment for the brine invasion and provide data support for water supply safety and water resource management.
Specifically, in an embodiment, as shown in fig. 14, the method further includes:
and processing the historical data information of the target estuary to obtain a data scatter diagram.
And carrying out fitting verification on the tide estuary saline water invasion model based on the data scatter diagram.
And judging whether the verification result meets the expectation, if not, optimizing the model until the model meets the expectation standard.
Specifically, tidal range plays a key role in promoting the over-standard salinity, and sea inflow plays a key role in resisting the over-standard salinity. Through verifying and optimizing the tidal estuary salt water invasion model, the model has higher accuracy and more accurate result, and provides more favorable data support for users of the model.
In this embodiment, a device for constructing a tidal estuary brine invasion model is further provided, and the device is used for implementing the above embodiments and preferred embodiments, and the description of the device is omitted. As used below, the term "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware, or a combination of software and hardware is also possible and contemplated.
The embodiment provides a device for constructing a tidal estuary brine invasion model, as shown in fig. 8, comprising:
the obtaining module 101 is configured to obtain data information of the target river mouth, where the data information includes a tidal range and a maximum saline invasion area corresponding to turning points of different sea inflow rates, and details of the data information refer to relevant description of step S1 in the foregoing method embodiment, and are not described herein again.
The first analysis module 102 is configured to determine a first sensitivity coefficient and a first relationship between the sea inflow rate and the tidal range according to the sea inflow rate, the tidal range, and the maximum saline invasion area, for details, refer to the related description of step S2 in the foregoing method embodiment, and details are not repeated here.
The second analysis module 103 is configured to determine a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area, for details, refer to the related description of step S3 in the foregoing method embodiment, and details are not repeated herein.
A building module 104, configured to build a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relationship, the second sensitivity coefficient, and the middle turning point, for details, refer to the related description of step S4 in the foregoing method embodiment, and details are not repeated here.
The means for constructing the tidal estuary brine intrusion model in this embodiment is in the form of a functional unit, where the unit is an ASIC circuit, a processor and memory executing one or more software or fixed programs, and/or other devices capable of providing the above-described functionality.
Further functional descriptions of the modules are the same as those of the corresponding embodiments, and are not repeated herein.
There is also provided an electronic device according to an embodiment of the present invention, as shown in fig. 9, the electronic device may include a processor 901 and a memory 902, where the processor 901 and the memory 902 may be connected by a bus or in another manner, and fig. 9 takes the example of being connected by a bus as an example.
The memory 902, which is a non-transitory computer readable storage medium, may be used for storing non-transitory software programs, non-transitory computer executable programs, and modules, such as program instructions/modules corresponding to the methods in the method embodiments of the present invention. The processor 901 executes various functional applications and data processing of the processor by executing non-transitory software programs, instructions and modules stored in the memory 902, that is, implements the methods in the above-described method embodiments.
The memory 902 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application program required for at least one function; the storage data area may store data created by the processor 901, and the like. Further, the memory 902 may include high speed random access memory, and may also include non-transitory memory, such as at least one magnetic disk storage device, flash memory device, or other non-transitory solid state storage device. In some embodiments, the memory 902 may optionally include memory located remotely from the processor 901, which may be connected to the processor 901 via a network. Examples of such networks include, but are not limited to, the internet, intranets, local area networks, mobile communication networks, and combinations thereof.
One or more modules are stored in the memory 902, which when executed by the processor 901 performs the methods in the above-described method embodiments.
The specific details of the electronic device may be understood by referring to the corresponding related descriptions and effects in the above method embodiments, and are not described herein again.
It will be understood by those skilled in the art that all or part of the processes of the methods of the embodiments described above can be implemented by hardware related to instructions of a computer program, and the program can be stored in a computer readable storage medium, and when executed, can include the processes of the embodiments of the methods described above. The storage medium may be a magnetic Disk, an optical Disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a Flash Memory (Flash Memory), a Hard Disk (Hard Disk Drive, abbreviated as HDD) or a Solid State Drive (SSD), etc.; the storage medium may also comprise a combination of memories of the kind described above.
Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the invention, and such modifications and variations fall within the scope defined by the appended claims.
Claims (10)
1. A method for constructing a tidal estuary brine invasion model is characterized by comprising the following steps:
acquiring data information of a target river mouth, wherein the data information comprises tidal ranges and maximum salt water invasion areas corresponding to turning points of different sea inflow rates;
determining a first sensitivity coefficient and a first relation between the sea inflow and the tidal difference according to the sea inflow, the tidal difference and the maximum saline invasion area;
determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area;
and constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point.
2. The method of constructing a tidal estuary saltwater intrusion model of claim 1, wherein said determining a first sensitivity factor and a first relationship between sea inflow and tidal range from said sea inflow, said tidal range and said maximum saltwater intrusion area comprises:
obtaining a first sensitivity coefficient between the tidal difference and the saline invasion area by carrying out calibration fitting on the maximum saline invasion area and the tidal difference;
and fitting the sea inflow rate and the tidal range to obtain a first relation between different sea inflow rates and tidal ranges.
3. The method of constructing a tidal estuary brine intrusion model of claim 2, wherein the obtaining a first relationship between different sea inflow and tidal differences by fitting the sea inflow and the tidal differences comprises:
calculating the average value of the sea inflow;
comparing the different sea inflow rates with the average value to obtain a plurality of dimensionless sea inflow rate data;
fitting the dimensionless inflow volume data and the corresponding tidal range to obtain a first relation between different inflow volumes and tidal ranges.
4. The method of constructing a tidal estuary saltwater intrusion model of claim 1, wherein said determining a second sensitivity factor and an intermediate turning point according to said sea inflow and said maximum saltwater intrusion area comprises:
obtaining a middle turning point of the sea inflow rate by rating fitting of the sea inflow rate and the maximum saline water invasion area;
and calculating to obtain a second sensitivity coefficient between the sea inflow and the saline invasion area according to the sea inflow, the maximum saline invasion area and the middle turning point.
5. The method of constructing a tidal estuary brine intrusion model of claim 4, wherein said obtaining the intermediate turning point of the sea inflow by ratiometric fitting of the sea inflow and the maximum brine intrusion area comprises:
calculating the average value of the sea inflow;
comparing the different sea inflow rates with the average value to obtain a plurality of dimensionless sea inflow rate data;
carrying out calibration fitting on the maximum saline water invasion area and the dimensionless inflow volume data to obtain a fitting curve;
and determining the abscissa value corresponding to the turning point in the fitting curve as the middle turning point of the sea inflow.
6. The method of constructing a tidal estuary brine intrusion model according to claim 1, wherein the constructing a tidal estuary brine intrusion model based on the maximum brine intrusion area, the first sensitivity coefficient, the first relationship, the second sensitivity coefficient and the intermediate turning point comprises:
constructing a relationship between saline invasion area and tidal range based on the maximum saline invasion area, the first sensitivity coefficient and the first relationship;
constructing a relation between the saline invasion area and the inflow flow based on the maximum saline invasion area, the second sensitivity coefficient and the middle turning point;
and obtaining a tidal estuary brine invasion model according to the relationship among the maximum brine invasion area, the brine invasion area and the tidal range and the relationship among the brine invasion area and the inflow flow.
7. The method of constructing a tidal estuary brine intrusion model of claim 1, further comprising:
acquiring sea inflow data of a target estuary in real time;
and substituting the inflow volume data into the tidal estuary brine invasion model to obtain a brine invasion area.
8. A device for constructing a tidal estuary brine invasion model is characterized by comprising:
the acquisition module is used for acquiring data information of a target estuary, wherein the data information comprises tidal ranges and maximum salt water invasion areas corresponding to turning points of different sea inflow rates;
a first analysis module for determining a first sensitivity coefficient and a first relationship between the sea inflow and the tidal range according to the sea inflow, the tidal range and the maximum saline invasion area;
the second analysis module is used for determining a second sensitivity coefficient and a middle turning point according to the sea inflow and the maximum saline invasion area;
and the construction module is used for constructing a tidal estuary brine invasion model based on the maximum brine invasion area, the first sensitivity coefficient, the first relation, the second sensitivity coefficient and the middle turning point.
9. An electronic device, comprising:
a memory and a processor, the memory and the processor being communicatively connected to each other, the memory having stored therein computer instructions, the processor executing the computer instructions to perform the method of constructing a tidal estuary brine intrusion model of any one of claims 1-7.
10. A computer-readable storage medium storing computer instructions for causing a computer to perform the method of constructing a tidal estuary brine intrusion model of any one of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210087093.0A CN114547992A (en) | 2022-01-25 | 2022-01-25 | Method and device for constructing tidal estuary brine invasion model and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210087093.0A CN114547992A (en) | 2022-01-25 | 2022-01-25 | Method and device for constructing tidal estuary brine invasion model and electronic equipment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114547992A true CN114547992A (en) | 2022-05-27 |
Family
ID=81672101
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210087093.0A Pending CN114547992A (en) | 2022-01-25 | 2022-01-25 | Method and device for constructing tidal estuary brine invasion model and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114547992A (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102590469A (en) * | 2012-01-16 | 2012-07-18 | 河海大学 | Flat plate tracing simulating device system and simulating method for salt water intrusion under tidal action |
| CN110110942A (en) * | 2019-05-20 | 2019-08-09 | 水利部交通运输部国家能源局南京水利科学研究院 | A kind of salinity effect method of riverine getting water from water head site mouth |
-
2022
- 2022-01-25 CN CN202210087093.0A patent/CN114547992A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102590469A (en) * | 2012-01-16 | 2012-07-18 | 河海大学 | Flat plate tracing simulating device system and simulating method for salt water intrusion under tidal action |
| CN110110942A (en) * | 2019-05-20 | 2019-08-09 | 水利部交通运输部国家能源局南京水利科学研究院 | A kind of salinity effect method of riverine getting water from water head site mouth |
Non-Patent Citations (1)
| Title |
|---|
| ZHI XU等: "Saltwater Intrusion Function and Preliminary Application in the Yangtze River Estuary, China", 《INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH》, pages 1 - 19 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2020103356A4 (en) | Method and device for building river diatom bloom warning model | |
| CN104318325B (en) | Many basin real-time intelligent water quality prediction methods and system | |
| Adnan et al. | Pan evaporation estimation by relevance vector machine tuned with new metaheuristic algorithms using limited climatic data | |
| US20180181876A1 (en) | Unsupervised machine learning to manage aquatic resources | |
| CN112288193A (en) | Ocean station surface salinity prediction method based on GRU deep learning of attention mechanism | |
| CN107292383A (en) | The variation water quality interval prediction method being combined based on deep learning algorithm with MILP | |
| CN109002888A (en) | A kind of red tide prewarning method | |
| Shang | Lake surface area method to define minimum ecological lake level from level-area-storage curves | |
| Liu et al. | Automatic and accurate prediction of key water quality parameters based on SRU deep learning in mariculture | |
| Bi et al. | Daily runoff forecasting based on data-augmented neural network model | |
| Zhang et al. | Fuzzy logic method for evaluating habitat suitability in an estuary affected by land reclamation | |
| CN113668479B (en) | Method and system for adjusting filter screen for sea surface cleaning ship | |
| Liang | Water pollution prediction based on deep belief network in big data of water environment monitoring | |
| CN114547992A (en) | Method and device for constructing tidal estuary brine invasion model and electronic equipment | |
| CN117253344A (en) | Seawater acidification early warning and forecasting method, system and electronic equipment | |
| Jiang et al. | COD optimization prediction model based on CAWOA-ELM in water ecological environment | |
| CN116559059A (en) | Tidal energy reservoir water dike chloride ion erosion detection method and device and electronic equipment | |
| Calderón-Vega et al. | Analysis of extreme waves with seasonal variation in the Gulf of Mexico using a time-dependent GEV model | |
| CN114971354A (en) | Water quality grade evaluation method, device and medium | |
| CN111027202B (en) | Digital city prediction method, device, equipment and storage medium | |
| CN114897319A (en) | Method for establishing bearing capacity evaluation and early warning index system of water resource environment | |
| CN114862062A (en) | Early warning method and device for environmental bearing capacity, electronic equipment and readable storage medium | |
| Eze et al. | Developing a Novel Water Quality Prediction Model for a South African Aquaculture Farm. Water 2021, 13, 1782 | |
| CN111859707B (en) | Pipeline diameter determining method, urban waterlogging early warning method and device | |
| Li et al. | Assessment and Analysis of Water Environmental Risk Using Pressure-state-response Model in Ningbo, China. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |