CN106844857B - Method and device for simulating time-space distribution of water pollution accident pollution zone - Google Patents

Method and device for simulating time-space distribution of water pollution accident pollution zone Download PDF

Info

Publication number
CN106844857B
CN106844857B CN201611190316.7A CN201611190316A CN106844857B CN 106844857 B CN106844857 B CN 106844857B CN 201611190316 A CN201611190316 A CN 201611190316A CN 106844857 B CN106844857 B CN 106844857B
Authority
CN
China
Prior art keywords
section
state variable
auxiliary
river
variable
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.)
Expired - Fee Related
Application number
CN201611190316.7A
Other languages
Chinese (zh)
Other versions
CN106844857A (en
Inventor
张波
郝千婷
闫其年
秦宇
郝莹
陈瀛
徐昊
李诗瑶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Information Center Of Ministry Of Ecological Environment
Original Assignee
Information Center Of Ministry Of Ecological Environment
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Information Center Of Ministry Of Ecological Environment filed Critical Information Center Of Ministry Of Ecological Environment
Priority to CN201611190316.7A priority Critical patent/CN106844857B/en
Publication of CN106844857A publication Critical patent/CN106844857A/en
Application granted granted Critical
Publication of CN106844857B publication Critical patent/CN106844857B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The application provides a method and a device for simulating the time-space distribution of a pollution zone of a water pollution accident, a system for executing the method is provided with a graphical modeling interface of a one-dimensional hydrodynamic model, and the method comprises the following steps: receiving state variables and auxiliary variables input by a user through the graphical modeling interface; establishing an incidence relation between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model; calculating the value of each state variable at the corresponding time of each section by using the one-dimensional hydrodynamic model by taking the obtained auxiliary variable parameter value and the estimated parameter value of other auxiliary variables in the credibility range as input parameters; and generating the simulation time for the water pollution accident pollution zone to reach each section according to the value of the corresponding time of each state variable on each section. Through the means, the problem that the time-space distribution of the water pollution accident pollution zone cannot be predicted timely and accurately in the prior art can be effectively solved.

Description

Method and device for simulating time-space distribution of water pollution accident pollution zone
Technical Field
The application relates to the technical field of river environment monitoring and forecasting, in particular to a method and a device for simulating time-space distribution of a water pollution accident pollution zone.
Background
In recent years, serious water pollution accidents in China frequently occur, and serious influences are caused on ecological environment, human health and social safety. For example, the water pollution accident of Songhua river, which occurs in month 11 of 2005, and the cadmium pollution accident of Guangdong and North river, which occurs in month 12, are typical cases of major water pollution accidents. In the emergency process of water pollution accidents, the time-space distribution of pollution zones and the concentration of pollutants are urgently needed to be mastered, so that an emergency response scheme can be quickly and effectively formulated. In order to meet the requirement of rapid disposal of water pollution accidents, a water quality mathematical model (called a water quality model for short) is used for timely and accurately predicting the time, the influence range and the like of a pollution group reaching each important section in combination with actual monitoring data of an accident site, the effects of various emergency strategies and schemes are simulated, a basis for quantitative reference is provided for emergency decision of accidents, and the method is an important way for improving the emergency capacity of environmental accidents.
However, the disadvantages of the existing water quality model system are: the regulation and control function of the model is weak, and the basis of quantitative reference is provided for prediction of water pollution accidents and optimization of emergency measures and schemes by adjusting various parameters difficultly; the model parameters are numerous, the data quantity required by parameter calibration is large, and the problem of difficult calibration of the model parameters caused by lack of monitoring data in a water pollution accident site is difficult to solve; the model operation is complicated, a user needs to spend a certain time to know the model, and a numerical solution meeting certain precision requirements cannot be quickly obtained through simple operation after an accident occurs.
Disclosure of Invention
The application provides a method and a device for simulating time-space distribution of a water pollution accident pollution zone, which are used for solving the problems that the existing water quality model system has numerous parameters, large data quantity to be acquired and complex operation, and cannot timely and accurately predict the time-space distribution of the water pollution accident pollution zone.
The application discloses a method for simulating time-space distribution of a water pollution accident pollution zone, a system for executing the method is provided with a graphical modeling interface of a one-dimensional hydrodynamic model, and the method comprises the following steps: receiving state variables and auxiliary variables input by a user through the graphical modeling interface; establishing an incidence relation between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model; calculating the value of each state variable at the corresponding time of each section by using the one-dimensional hydrodynamic model by taking the obtained auxiliary variable parameter value and the estimated parameter value of other auxiliary variables in the credibility range as input parameters; and generating the simulation time for the water pollution accident pollution zone to reach each section according to the value of the corresponding time of each state variable on each section.
The application discloses a device of water pollution accident pollution area space-time distribution is simulated, includes: the graphical modeling interface is used for receiving state variables and auxiliary variables required by building the one-dimensional hydrodynamic model; the model establishing module is used for establishing the incidence relation between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model; the operation execution module is used for calculating the value of each state variable at the corresponding time of each section by using the one-dimensional hydrodynamic model by taking the acquired auxiliary variable parameter value and the estimated parameter value of other auxiliary variables in the credibility range as input parameters; and the result generation module is used for generating simulation time for the water pollution accident pollution zone to reach each section according to the value of the corresponding time of each state variable on each section.
Compared with the prior art, the method has the following advantages: the method has the advantages that the one-dimensional hydrodynamic model is established in the preferred embodiment of the application, the parameter insensitivity is realized, when the influence factor (namely auxiliary variable) data of the state variable is insufficient, the same system trend, behavior mode and fluctuation period can still be shown as long as the parameter estimation value can be within the credibility range, the method can be well suitable for the condition that real-time monitoring data are difficult to obtain due to the sudden water pollution accident, and the problem that the time-space distribution of the water pollution accident pollution zone cannot be timely and accurately predicted in the prior art can be effectively solved.
Drawings
The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the application. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a flow chart of a method for simulating a time-space distribution of a pollution zone of a water pollution accident according to an embodiment of the present application;
FIG. 2 is a schematic block diagram of a hydrodynamic model of the system used in the embodiment shown in FIG. 1;
FIG. 3 is a comparison of an actual cross-section of a river and a simulated cross-section;
FIG. 4 is a schematic diagram of the gradient of a river bed;
FIG. 5 is a schematic structural diagram of an apparatus for simulating spatial and temporal distribution of a pollution zone in a water pollution accident according to an embodiment of the present disclosure.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, the present application is described in further detail with reference to the accompanying drawings and the detailed description.
In the description of the present application, it is to be understood that the terms "first", "second" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implying any number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. The meaning of "plurality" is two or more unless specifically limited otherwise. The terms "comprising," including, "and the like are to be construed as open-ended terms, i.e.," including/including but not limited to. The term "based on" is "based, at least in part, on". The term "an embodiment" means "at least one embodiment"; the term "another embodiment" means "at least one additional embodiment". Relevant definitions for other terms will be given in the following description.
Referring to fig. 1, a flow of a method for simulating a time-space distribution of a pollution zone of a water pollution accident according to an embodiment of the present application is shown, in order to facilitate a modeling operation of a user, a system for executing the method is provided with a graphical modeling interface of a one-dimensional hydrodynamic model, through which the user can add icons of state variables and auxiliary variables in a dragging manner and establish a connection relationship between the variables. The method comprises the following steps:
step S101: and receiving the state variable and the auxiliary variable input by the user through the graphical modeling interface.
Step S102: and establishing an incidence relation between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model.
Step S103: and calculating the value of each state variable at the corresponding time of each section by using the one-dimensional hydrodynamic model by taking the acquired auxiliary variable parameter value and the estimated parameter value of other auxiliary variables in the credibility range as input parameters.
Step S104: and generating the simulation time for the water pollution accident pollution zone to reach each section according to the value of the corresponding time of each state variable on each section.
In another embodiment, when the user is presentThe selected state variables are river depth H and river flow Q (the corresponding initial values are the water depth and the flow of the section where the water pollution accident occurs), and the auxiliary variables are the unit branch inflow QinUnit substream outflow qoutWidth of river surface BwFlow velocity unevenness coefficient β, river cross section area A, gravity acceleration g, and bed gradient SxRiverbed height BhRiverbed frictional resistance SfMean velocity of flow and dispersion coefficient of fracture ExHydraulic radius R, Manning roughness coefficient n and average branch flow velocity of branch flow velocity in river direction
Figure GDA0001257337970000032
The relationship between the state variable H, Q and the auxiliary variable may be established using the following formula:
Figure GDA0001257337970000031
Figure GDA0001257337970000041
in the above formula, x represents the distance between the cross section and the place where the water pollution accident occurs, Δ x represents the distance between adjacent cross sections, and t represents the time length after the water pollution accident occurs; the process of establishing the association relationship between the auxiliary variables specifically includes:
in general, the length L of the river is much greater than the width Bw (i.e., Bw)<<L). For example, Songhua river is about 500km long (L) from Harbin city to Betula section, while river width does not exceed 1km (B)w) The narrowest part (canada river reach) is only 405m, and most of the narrowest part is 600-900 m. In addition, most of Chinese rivers are shallow rivers, the water depth is generally 2-3 m, and the water depth H is far less than the river width Bw(i.e. H)<<Bw) Therefore, the cross section of the river can be assumed to be rectangular (as shown in FIG. 3), and the cross section A of the river cross section can be expressed as the river depth H and the river surface width BwThe form of the product, namely:
A=H·Bw
the slope angle of the river bed is generally very small, as shown in fig. 4, the river bed is taken as an x axis and points to the flowing direction of the river, a y axis is the width direction of the river, a z axis is vertical to the river bed and belongs to the water depth direction, the slope angle α of the river bottom is less than 1, sin α is equal to tg αxSlope S of riverbedxCan be expressed in the following form:
Figure GDA0001257337970000042
the water flow velocity of the river section is unevenly distributed due to the turbulence effect, and the dispersion coefficient E can be adoptedxThe nonuniformity of the water velocity caused by the turbulent flow is corrected. The application adopts an approximate estimation formula for irregular natural rivers proposed by Feichel:
Figure GDA0001257337970000043
in the formula u*Represents the flow resistance (in m/s) by friction,
Figure GDA0001257337970000044
in specific implementation, the gravity acceleration g can be 9.8m/s2Coefficient of dispersion ExThe values can be: e is not less than 0x≤100km2D (i.e., square kilometers per day).
The flow velocity u is different at each point on the cross section, and the non-uniformity coefficient β of the flow velocity can reflect the non-uniformity of the flow velocity u on the cross section, and can be specifically calculated by the following formula:
Figure GDA0001257337970000051
wherein rho represents the density of water and can be taken as value rho in specific implementationWater (W)=1000kg/m3
In the case of non-constant flow, the resistance to which the flow is subjected is not much different from that of constant flow, whereas in the case of constant flow, the frictional resistance SfThe river water push is triggered by the gradient of the riverbedThe power f is balanced. The frictional resistance S is preferably calculated by the following Manning formulaf
Figure GDA0001257337970000052
In the formula (I), the compound is shown in the specification,
mean velocity of cross section
Figure GDA0001257337970000053
The larger the Manning coefficient n is, the larger the riverbed friction is, and the relative flow velocity value is reduced correspondingly. The value of the Mannich coefficient n is generally measured by experimental data, and can be selected by looking up a table when in use.
The hydraulic radius R is the section area A (A is H.B)w) Divided by the wet perimeter (2H + 2B)w) Namely:
Figure GDA0001257337970000054
the form of the one-dimensional hydrodynamic model created by the above process is shown in fig. 2. In the model shown in fig. 2, two branch flows are provided for the inflow and two branch flows for the outflow, and if the number of branch flows is increased, it is possible to easily expand the model. Wherein:
inite 1 and inite 2 represent the positions of two inflowing branches, qin1 and qin2 represents the unit flow rate of the two branches (the water discharge of the reservoir causes the flow rate to increase in a certain period); q. q.sin1=Qin1/Bw;qin2=Qin2/Bw
outnite 1 and outnite 2 represent the positions of two branches flowing out of the main river reach respectively, and q isout1 and qout2 represents the unit flow rate of the outflow of the two branches; q. q.sout1=Qout1/Bw;qout2=Qout2/Bw
In a further preferred embodiment, the simulation time for the water pollution accident pollution zone to reach each section is calculated by adopting the following formula:
Figure GDA0001257337970000061
wherein, when x is 0, tx=0。
In another further preferred embodiment, the one-dimensional hydrodynamic model may further include a rate control parameter of a state variable, wherein:
Figure GDA0001257337970000062
Figure GDA0001257337970000063
Figure GDA0001257337970000064
in the above equation, hr represents the rate control parameter for the state variable H, slope represents the first rate control parameter for the state variable Q, and qr represents the second rate control parameter for the state variable Q.
For simplicity of description, the foregoing method embodiments are described as a series of acts or combination of acts, but those skilled in the art will appreciate that the present application is not limited by the order of acts described, as some steps may, in accordance with the present application, occur in other orders and concurrently; further, those skilled in the art should also appreciate that the above-described method embodiments are preferred embodiments and that the acts and modules involved are not necessarily required for the application.
Referring to fig. 5, a schematic structural diagram of an apparatus for simulating a spatial-temporal distribution of a pollution zone in a water pollution accident according to an embodiment of the present application is shown, the apparatus including:
and the graphical modeling interface 51 is used for receiving the state variables and the auxiliary variables required by establishing the one-dimensional hydrodynamic model.
And the model creating module 52 is configured to create an association relationship between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model.
When the state variables input by the user are river depth H and river flow Q, the auxiliary variable is unit branch inflow QinUnit substream outflow qoutWidth of river surface BwFlow velocity unevenness coefficient β, river cross section area A, gravity acceleration g, and bed gradient SxRiverbed height BhRiverbed frictional resistance SfMean velocity of flow and dispersion coefficient of fracture ExHydraulic radius R, Manning roughness coefficient n and average branch flow velocity of branch flow velocity in river direction
Figure GDA0001257337970000071
Model creation module 52 may then establish the relationship between state variable H, Q and the auxiliary variables using the following formula:
Figure GDA0001257337970000072
Figure GDA0001257337970000073
wherein x represents the distance between the cross section and the water pollution accident site, Δ x represents the spacing distance between adjacent cross sections, and t represents time; the relationship between the auxiliary variables is:
A=H·Bw
Figure GDA0001257337970000074
Figure GDA0001257337970000075
the diffusion coefficient ExThe value range of (a) can be set to 0-100 square kilometers per day.
Figure GDA0001257337970000076
Figure GDA0001257337970000077
Where ρ represents the density of water and u represents the point flow velocity on the cross section.
And the operation execution module 53 is configured to use the acquired auxiliary variable parameter values and estimated parameter values of other auxiliary variables within the reliability range as input parameters, and calculate values of each state variable at corresponding time of each section by using the one-dimensional hydrodynamic model.
And the result generation module 54 is used for generating the simulation time for the water pollution accident pollution zone to reach each section according to the value of the corresponding time of each state variable in each section.
In specific implementation, the result generating module 54 may specifically calculate the simulation time of the water pollution accident pollution zone reaching each section by using the following formula:
Figure GDA0001257337970000081
wherein, when x is 0, tx=0。
In another embodiment, the one-dimensional hydrodynamic model may further include a rate control parameter of the state variable, wherein:
Figure GDA0001257337970000082
Figure GDA0001257337970000083
Figure GDA0001257337970000084
in the above equation, hr represents the rate control parameter for the state variable H, slope represents the first rate control parameter for the state variable Q, and qr represents the second rate control parameter for the state variable Q.
It should be noted that the above device embodiments belong to preferred embodiments, and the units and modules involved are not necessarily essential to the present application.
The embodiments in the present specification are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the device embodiments of the present application, since they are substantially similar to the method embodiments, the description is relatively simple, and for the relevant points, reference may be made to the description of the method embodiments.
The method and the device for simulating the space-time distribution of the pollution zone of the water pollution accident are introduced in detail, specific examples are applied in the method for explaining the principle and the implementation mode of the method, and the description of the examples is only used for helping to understand the method and the core idea of the method; meanwhile, for a person skilled in the art, according to the idea of the present application, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present application.

Claims (8)

1. A method of simulating the spatiotemporal distribution of the pollution zone of a water pollution accident, characterized in that the system implementing said method is provided with a graphical modeling interface of a one-dimensional hydrodynamic model, said method comprising:
receiving state variables and auxiliary variables input by a user through the graphical modeling interface;
establishing an incidence relation between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model;
calculating the value of each state variable at the corresponding time of each section by using the one-dimensional hydrodynamic model by taking the obtained auxiliary variable parameter value and the estimated parameter value of other auxiliary variables in the credibility range as input parameters;
generating simulation time for the water pollution accident pollution zone to reach each section according to the value of each state variable at the corresponding time of each section;
the state variables comprise river water depth H and river flow Q, and the initial values are water depth and flow of a section where the water pollution accident occurs respectively;
the auxiliary variable comprising the unit branch inflow qinUnit substream outflow qoutWidth of river surfaceBwFlow velocity unevenness coefficient β, river cross section area A, gravity acceleration g, and bed gradient SxRiverbed height BhRiverbed frictional resistance SfMean velocity of flow and dispersion coefficient of fracture ExHydraulic radius R, Manning roughness coefficient n and average branch flow velocity of branch flow velocity in river direction
Figure 573086DEST_PATH_IMAGE002
The relationship between the state variable H, Q and the auxiliary variables is established using the following formula:
Figure 490226DEST_PATH_IMAGE004
Figure 332280DEST_PATH_IMAGE006
wherein x represents the distance between the section and the water pollution accident site, Δ x represents the distance between adjacent sections, and t represents time; the relationship between the auxiliary variables is:
A=H·Bw
Figure 727490DEST_PATH_IMAGE008
Figure 541862DEST_PATH_IMAGE010
Figure 364324DEST_PATH_IMAGE012
Figure 942942DEST_PATH_IMAGE014
Figure 204159DEST_PATH_IMAGE016
Figure 873038DEST_PATH_IMAGE018
Figure 866402DEST_PATH_IMAGE020
(ii) a Where ρ represents the density of water and u represents the point flow velocity on the cross section.
2. The method of claim 1, wherein the simulated time for the water pollution accident zone to reach each section is calculated by the following formula:
t(x+∆x)=tx+ ∆x/ūx
wherein, when x =0, tx=0。
3. The method of claim 1, wherein the one-dimensional hydrodynamic model further comprises a rate control parameter for a state variable, wherein:
Figure 433780DEST_PATH_IMAGE022
Figure 170792DEST_PATH_IMAGE024
Figure 694178DEST_PATH_IMAGE026
in the above equation, hr represents the rate control parameter for the state variable H, slope represents the first rate control parameter for the state variable Q, and qr represents the second rate control parameter for the state variable Q.
4. The method of claim 1, wherein the diffusion coefficient ExThe value range of (A) is 0-100 square kilometers per day.
5. A device for simulating the space-time distribution of a water pollution accident pollution zone is characterized by comprising:
the graphical modeling interface is used for receiving state variables and auxiliary variables required by building the one-dimensional hydrodynamic model;
the model establishing module is used for establishing the incidence relation between the state variable and the auxiliary variable to form a one-dimensional hydrodynamic model;
the operation execution module is used for calculating the value of each state variable at the corresponding time of each section by using the one-dimensional hydrodynamic model by taking the acquired auxiliary variable parameter value and the estimated parameter value of other auxiliary variables in the credibility range as input parameters;
the result generation module is used for generating simulation time for the water pollution accident pollution zone to reach each section according to the value of each state variable at the corresponding time of each section;
the state variables comprise river water depth H and river flow Q, and the initial values are water depth and flow of a section where the water pollution accident occurs respectively;
the auxiliary variable comprising the unit branch inflow qinUnit substream outflow qoutWidth of river surface BwFlow velocity unevenness coefficient β, river cross section area A, gravity acceleration g, and bed gradient SxRiverbed height BhRiverbed frictional resistance SfMean velocity of flow and dispersion coefficient of fracture ExHydraulic radius R, Manning roughness coefficient n and average branch flow velocity of branch flow velocity in river direction
Figure 124022DEST_PATH_IMAGE002
The model creation module establishes the association between the state variable H, Q and the auxiliary variable using the following formula:
Figure 427964DEST_PATH_IMAGE004
Figure 703088DEST_PATH_IMAGE006
wherein x represents the distance between the section and the water pollution accident site, Δ x represents the distance between adjacent sections, and t represents time; the relationship between the auxiliary variables is:
A=H·Bw
Figure 346559DEST_PATH_IMAGE008
Figure 681725DEST_PATH_IMAGE010
Figure 410647DEST_PATH_IMAGE012
Figure 818624DEST_PATH_IMAGE014
Figure 316601DEST_PATH_IMAGE016
Figure 88248DEST_PATH_IMAGE018
Figure 101204DEST_PATH_IMAGE020
(ii) a Where ρ represents the density of water and u represents the point flow velocity on the cross section.
6. The apparatus of claim 5, wherein the result generation module calculates the simulated time for the water pollution accident zone to reach each section by using the following formula:
t(x+∆x)=tx+ ∆x/ūx
wherein, when x =0, tx=0。
7. The apparatus of claim 5, wherein the one-dimensional hydrodynamic model further comprises a rate control parameter for a state variable, wherein:
Figure 983709DEST_PATH_IMAGE022
Figure 149242DEST_PATH_IMAGE024
Figure 91790DEST_PATH_IMAGE026
in the above equation, hr represents the rate control parameter for the state variable H, slope represents the first rate control parameter for the state variable Q, and qr represents the second rate control parameter for the state variable Q.
8. The device of claim 5, wherein the diffusion coefficient ExThe value range of (A) is 0-100 square kilometers per day.
CN201611190316.7A 2016-12-21 2016-12-21 Method and device for simulating time-space distribution of water pollution accident pollution zone Expired - Fee Related CN106844857B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611190316.7A CN106844857B (en) 2016-12-21 2016-12-21 Method and device for simulating time-space distribution of water pollution accident pollution zone

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611190316.7A CN106844857B (en) 2016-12-21 2016-12-21 Method and device for simulating time-space distribution of water pollution accident pollution zone

Publications (2)

Publication Number Publication Date
CN106844857A CN106844857A (en) 2017-06-13
CN106844857B true CN106844857B (en) 2020-06-19

Family

ID=59135260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611190316.7A Expired - Fee Related CN106844857B (en) 2016-12-21 2016-12-21 Method and device for simulating time-space distribution of water pollution accident pollution zone

Country Status (1)

Country Link
CN (1) CN106844857B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109558652B (en) * 2018-11-15 2020-11-27 同济大学 Auxiliary correction method for river network water quality model calculation result
CN110008632B (en) * 2019-04-18 2023-04-18 江西师范大学 Method for rapidly calculating emergency pollution event of river channel type drinking water source
CN110866367B (en) * 2019-08-13 2021-07-27 河海大学 Real-time tracking method for conventional sudden water pollution clusters under complex river network water flow condition
CN112132264B (en) * 2020-09-11 2023-04-07 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) Regional exhaust emission prediction method and system based on space-time residual perception network
CN116258403B (en) * 2023-01-19 2023-09-15 河南省新乡水文水资源勘测局 Big data-based water quality monitoring and early warning method and system

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7254564B2 (en) * 2001-10-22 2007-08-07 Coppola Jr Emery J Neural network based predication and optimization for groundwater/surface water system
CN102608282B (en) * 2011-01-21 2015-03-25 环境保护部信息中心 Water quality information computing device and method
CN106168991B (en) * 2016-06-24 2019-03-15 珠江水利委员会珠江水利科学研究院 A kind of tidal river network tide prediction method based on hydrodynamic simulation

Also Published As

Publication number Publication date
CN106844857A (en) 2017-06-13

Similar Documents

Publication Publication Date Title
CN106844857B (en) Method and device for simulating time-space distribution of water pollution accident pollution zone
Bertoldi et al. Planform dynamics of braided streams
Gaudio et al. Non-universality of von Kármán's κ in fluvial streams
Dhar et al. Optimal operation of reservoirs for downstream water quality control using linked simulation optimization
Constantinescu LE of shallow mixing interfaces: A review
An et al. Three-dimensional modeling of turbid density currents in Imha Reservoir, South Korea
Bazin et al. Computing flooding of crossroads with obstacles using a 2D numerical model
Becker et al. Modelling the effects of thermal stratification on the oxygen budget of an impounded river
Momplot et al. Typology of the flow structures in dividing open channel flows
Gulley et al. Large values of hydraulic roughness in subglacial conduits during conduit enlargement: implications for modeling conduit evolution
CN106446438B (en) A kind of nearly far field approach of coupled numerical simulation of power plant&#39;s warm water discharge deep water discharge
Sauida Calibration of submerged multi-sluice gates
Hohermuth et al. Numerical simulation of air entrainment in uniform chute flow
KR20090093425A (en) Method for monitoring and modeling reservoir turbidity in real time
Ferrara et al. Flow field in shallow reservoir with varying inlet and outlet position
Xie et al. Large-eddy simulation of turbulent free surface flow over a gravel bed
Chen et al. Experimental and numerical investigation on the aerodynamics of isolated high-rise building and phenomenology of twisted wind field
Padhi et al. Conditional turbulence characteristics in water-worked and screeded gravel-bed flows
CN104657783B (en) For the section sliding-model control method of Mixing Coefficient in Rectangular Channels Release and dispersion prediction
Lopes et al. Numerical and experimental study of the fundamental flow characteristics of a 3D gully box under drainage
Wang et al. Model development in OpenFOAM to predict spillway jet regimes
Rao et al. Coarse-grained modeling of sheared granular beds
Al-Khatib et al. Development of empirical regression-based models for predicting mean velocities in asymmetric compound channels
Mark et al. An analytical model for solute mixing in surcharged manholes
Witt et al. Predicting oxygen transfer efficiency at low-head gated sill structures

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 100029 No. 1 Yu Hui South Road, Beijing, Chaoyang District

Applicant after: Information center of the Ministry of ecological environment

Address before: 100029 No. 1 Yu Hui South Road, Beijing, Chaoyang District

Applicant before: MINISTRY OF ENVIRONMENTAL PROTECTION INFORMATION CENTER

GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200619

Termination date: 20201221