CN115017727A - Sewage collection simulation method based on Masjing root method - Google Patents

Sewage collection simulation method based on Masjing root method Download PDF

Info

Publication number
CN115017727A
CN115017727A CN202210746387.XA CN202210746387A CN115017727A CN 115017727 A CN115017727 A CN 115017727A CN 202210746387 A CN202210746387 A CN 202210746387A CN 115017727 A CN115017727 A CN 115017727A
Authority
CN
China
Prior art keywords
section
equation
flow
concentration
masjing
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.)
Granted
Application number
CN202210746387.XA
Other languages
Chinese (zh)
Other versions
CN115017727B (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN202210746387.XA priority Critical patent/CN115017727B/en
Publication of CN115017727A publication Critical patent/CN115017727A/en
Application granted granted Critical
Publication of CN115017727B publication Critical patent/CN115017727B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • 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/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • 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
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152—Water filtration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a sewage collection simulation method based on a Masjing root method, which comprises the steps of (1) collecting river channel section flow data and initial moment head section concentration, and setting various parameters in the Masjing root method; (2) dispersing the basic equation of the one-dimensional water quality model by a Maskikyo method to obtain a dispersion equation; (3) solving the discrete equation to obtain the flow of each time interval and each section of the river channel; (4) calculating the area of each time interval and each section according to the flow, and each parameter in the discrete equation; (5) and substituting various parameters in the Masjing root method back into the discrete equation, and solving the migration and transformation process of the pollutant concentration in the river channel along with the time. The method reduces the data requirement of pollutant migration and transformation simulation, and expands the application scope of the Mas Jing root method.

Description

Sewage collection simulation method based on Masjing root method
Technical Field
The invention relates to a water quality simulation method, in particular to a sewage collection simulation method based on a Masjing root method.
Background
The migration and transformation process of the pollutant concentration in the river channel along with time and space is generally solved by using a one-dimensional water quality model basic equation, but the method needs a large amount of hydrodynamic data. The data is obtained by solving a one-dimensional Saint-Venn equation set, and river section data, boundary water level, flow and other data are generally needed, and the data is difficult to obtain and is often accompanied with data loss, so that the result precision is influenced. Meanwhile, the process of solving the Saint-Vietnam partial differential equation set is complicated, the data size is large, and the possibility of artificial errors in the actual working process is increased.
The Masjing root method is a hydrology method for river flood calculation, is different from the method for solving the Saint-Vietnam equation set, is simple, convenient and easy to operate, and is suitable for confluence simulation in hilly areas, but the method can only forecast section flow and cannot solve the problem of migration and transformation simulation of pollutants.
Disclosure of Invention
The purpose of the invention is as follows: aiming at the problems, the invention provides a sink pollution simulation method based on the Masjing root method, which reduces the data requirement of pollutant migration and conversion simulation and expands the application range of the Masjing root method in water environment simulation.
The technical scheme is as follows: the technical scheme adopted by the invention is a dirt collection simulation method based on the Masjing root method, which comprises the following steps:
(1) collecting river channel section flow data and initial moment first section pollutant concentration, and setting various parameters in the Mas Jing root method;
(2) dispersing the basic equation of the one-dimensional water quality model by a Masjing root method to obtain a dispersion equation; the discrete equation is as follows:
Figure BDA0003717466350000011
wherein:
Figure BDA0003717466350000012
in the formula, theta is a spatial term weighting coefficient, a variable superscript n represents a current time value, and j is a section number; Δ t is the time step, K is the empirical coefficient of Mas Jing, β, d 1 ，d 2 ，d 3 ，d 4 For calculating parameters, A is the area of the section; q is the flow; Δ x is the section distance; c is the concentration of the contaminant; s is a summary item; alpha is a degradation coefficient; q is a side single wide inflow; c. C 0 Is the side inflow concentration, K c W is the tank storage amount, and W is KQ.
(3) Solving the discrete equation to obtain the flow of each time interval and each section of the river channel; solving the discrete equation to obtain the flow of each time interval and each section of the river channel, namely solving the parameters beta and d in the discrete equation 1 ，d 2 ，d 3 ，d 4 The method is divided into the sum of two parts of an invariant parameter and a variable parameter. When the variable parameter is expressed as γ, the split discrete equation is:
Figure BDA0003717466350000021
wherein:
Figure BDA0003717466350000022
Figure BDA0003717466350000023
q′＝qΔx
Figure BDA0003717466350000024
S′ q ＝q′C 0
wherein gamma is a calculation parameter, S' q The convergence rate of the side inflow pollutants is represented, theta is a space term weighting coefficient, a variable superscript n represents a current time value, and j is a section number; Δ t is the time step, K is the empirical coefficient of Mas Jing, β, d 1 ，d 2 ，d 3 ，d 4 For calculating parameters, A is the area of the section; q is the flow; Δ x is the section distance; c is the concentration of the contaminant; s is a summary item; alpha is a degradation coefficient; q is a side single wide inflow; c. C 0 Is the side inflow concentration, K c W is the tank storage amount, and W is KQ.
(4) Calculating the area of each time interval and each section according to the flow and each parameter in the discrete equation;
(5) and substituting various parameters in the Masjing root method back into the discrete equation, and solving the pollutant concentration and the migration and transformation process of the pollutant concentration in the river channel along with time. Wherein the pollutant concentration is the pollutant concentration of the main trunk, and the calculation formula is as follows:
Figure BDA0003717466350000031
in the formula, D in Is the crossover point pollutant discharge rate.
Has the advantages that: compared with the prior art, the invention has the following advantages: the data requirement of pollutant migration and conversion simulation is reduced, the pollutant concentration simulation of different sections in the river channel is realized, and the dynamic process of pollutant migration in the river channel is reflected.
Drawings
FIG. 1 is a process flow diagram of a dirt collection simulation method based on the Masjing root method according to the present invention;
fig. 2 is a graph showing the simulation result of the present invention.
Detailed Description
The technical solution of the present invention is further described below with reference to the accompanying drawings and examples.
The invention relates to a dirt collection simulation method based on a Masjing root method, a flow chart of which is shown in figure 1, and the method comprises the following steps:
first, the convergent flow is processed by Masjing's method
Figure BDA0003717466350000032
Assuming that the river reach channel storage amount W has the following linear relationship with the flow rate Q and the inflow rate I:
W＝K[XI++(1-X)Q] (1)
in the formula, K and X are empirical coefficients, and X is more than or equal to 0 and less than or equal to 0.5.
Substituting the formula into the water balance equation can obtain:
Q 2 ＝C 1 I′ 1 +C 2 I′ 2 +C 3 Q 1 (2)
in the formula:
Figure BDA0003717466350000033
Figure BDA0003717466350000034
Figure BDA0003717466350000035
Figure BDA0003717466350000036
I i ′＝I i +qΔxi＝1，2 (7)
calculating the river channel dirt accumulation:
the basic equation of the one-dimensional water quality model is as follows:
Figure BDA0003717466350000041
wherein A is the cross-sectional area, m 2 (ii) a Q is the flow, m 3 S; x is the length of the river channel, m; c is the concentration of the contaminant, g/m 3 (ii) a S is a summary item, g/m/S; alpha is degradation coefficient 1/s; q is a side single wide inflow, m 2 /s；c 0 Is the side inflow concentration g/m 3 。
Figure BDA0003717466350000042
The two ends of the equation are divided by the area A at the same time, and the simplification is as follows:
Figure BDA0003717466350000043
the method comprises the following steps of dispersing a one-dimensional water quality model basic equation according to a differential format of a MasJinggen method, wherein a space term weighting coefficient is theta, a time term adopts a central differential format, a variable superscript n represents a current time value, a variable superscript n +1 represents a next time value, and j is a section number, so that the following form is obtained:
Figure BDA0003717466350000044
thus, there are:
Figure BDA0003717466350000045
order:
Figure BDA0003717466350000046
multiplying both sides of the equation by K c Δ t, and finishing to obtain:
Figure BDA0003717466350000047
in the formula:
Figure BDA0003717466350000051
wherein, β, d 1 ，d 2 ，d 3 ，d 4 For calculating the parameter, W is the tank storage amount, and W ═ KQ.
Considering that the concentration of pollutants contained in the side inflow influences the concentration of pollutants in the main trunk, the calculation is carried out by establishing a cross point:
Figure BDA0003717466350000052
in the formula, D in Is the crossover point pollutant discharge rate, g/s. i denotes different river numbers, c i Concentration of contaminant in the ith river, Q i The flow of the ith river channel. The physical meaning is as follows: and comparing the sum of the cross point pollutant discharge rate and all the side river pollutant discharge rates with the sum of all the side river flows, and finally calculating to obtain the pollutant concentration in the main road.
Beta, d is 1 ，d 2 ，d 3 ，d 4 The split was divided into the sum of two parts, the first half consisting of invariant parameters and the second half consisting of variable parameters and denoted γ.
Figure BDA0003717466350000053
In the formula:
Figure BDA0003717466350000054
Figure BDA0003717466350000061
q′＝qΔx (19)
Figure BDA0003717466350000062
S′ q ＝q′C 0 (21)
S′ q expressed as sidestream contaminant sink flow rate, g/s. Here only as a calculation parameter.
As discussed above, A' of different time intervals and sections is subject to flow
Figure BDA0003717466350000063
The influence of (c). Therefore, the flow calculated in the step (1) is substituted into the formula (17), gamma is calculated, and the pollutant concentration of the next section is solved through the formula (16)
Figure BDA0003717466350000064
Case description
The method is used for carrying out water quality simulation calculation on a certain river reach. The time interval Δ t was 3600s, and the cross-section interval Δ x was 9000m, for a total of 9 cross-sections. Wherein, the empirical parameter K is 3600s, X is 0.3, the spatial term weighting coefficient is θ 0.3, and the degradation coefficient α is 7.0 × 10 -7 1/s, side single wide inflow q is 0.001m 2 /s，c 0 Side inflow concentration c 0 ＝5mg/L。
And (3) obtaining a pollutant migration and transformation process line through a differential solution of a MaskAccu method, wherein the result is shown in figure 2.

Claims (5)

1. A sink pollution simulation method based on the Masjing root method is characterized by comprising the following steps:
(1) collecting river channel section flow data and initial moment first section pollutant concentration, and setting various parameters in the Mas Jing root method;
(2) dispersing the basic equation of the one-dimensional water quality model by a Masjing root method to obtain a dispersion equation;
(3) solving the discrete equation to obtain the flow of each time interval and each section of the river channel;
(4) calculating the area of each time interval and each section according to the flow and each parameter in the discrete equation;
(5) and substituting various parameters in the Masjing root method back into the discrete equation, and solving the pollutant concentration and the migration and transformation process of the pollutant concentration in the river channel along with time.
2. The sink contamination simulation method based on the masjing root method according to claim 1, characterized in that: the discrete equation in the step (2) is as follows:
Figure FDA0003717466340000011
wherein:
Figure FDA0003717466340000012
in the formula, theta is a spatial term weighting coefficient, a variable superscript n represents a current time value, and j is a section number; Δ t is the time step, K is the empirical coefficient of Mas Jing, β, d 1 ，d 2 ，d 3 ，d 4 For calculating parameters, A is the area of the cross section; q is the flow; Δ x is the section distance; c is the concentration of the contaminant; s is a sink item; alpha is a degradation coefficient; q is a side single wide inflow; c. C 0 Is the side inflow concentration, K c W is the tank storage amount, and W is KQ.
3. The sink contamination simulation method based on the masjing root method according to claim 1, characterized in that: in the step (3), the discrete equation is solved to obtain the flow of each time interval and each section of the river channel, namely, the parameters beta and d in the discrete equation are used 1 ，d 2 ，d 3 ，d 4 The method is divided into the sum of two parts of an invariant parameter and a variable parameter.
4. The method of claim 3 for simulating the sinking of pollutants based on the Masjing's methodIs characterized in that: the parameters beta, d in the discrete equation 1 ，d 2 ，d 3 ，d 4 Splitting the variable parameter into the sum of two parts of an invariant parameter and a variable parameter, and expressing the variable parameter as gamma, wherein the split discrete equation is as follows:
Figure FDA0003717466340000013
wherein:
Figure FDA0003717466340000021
Figure FDA0003717466340000022
q′＝qΔx
Figure FDA0003717466340000023
S′ q ＝q′C 0
wherein gamma is a calculation parameter, S' q Expressing the convergence rate of the side inflow pollutants, theta is a space term weighting coefficient, a variable superscript n expresses a current time value, and j is a section number; Δ t is the time step, K is the empirical coefficient of Mas Jing, β, d 1 ，d 2 ，d 3 ，d 4 For calculating parameters, A is the area of the section; q is the flow; Δ x is the section distance; c is the concentration of the contaminant; s is a summary item; alpha is a degradation coefficient; q is a side single wide inflow; c. C 0 Is the side inflow concentration, K c W is the tank storage amount, and W is KQ.
5. The sink contamination simulation method based on the masjing root method according to claim 1, characterized in that: the pollutant concentration in the step (5) is the pollutant concentration of the main trunk, and the calculation formula is as follows:
Figure FDA0003717466340000024
in the formula, D in Is the crossover point pollutant discharge rate.
CN202210746387.XA 2022-06-28 2022-06-28 Sewage collection simulation method based on Ma Sijing root method Active CN115017727B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210746387.XA CN115017727B (en) 2022-06-28 2022-06-28 Sewage collection simulation method based on Ma Sijing root method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210746387.XA CN115017727B (en) 2022-06-28 2022-06-28 Sewage collection simulation method based on Ma Sijing root method

Publications (2)

Publication Number Publication Date
CN115017727A true CN115017727A (en) 2022-09-06
CN115017727B CN115017727B (en) 2023-04-21

Family

ID=83076252

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210746387.XA Active CN115017727B (en) 2022-06-28 2022-06-28 Sewage collection simulation method based on Ma Sijing root method

Country Status (1)

Country Link
CN (1) CN115017727B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106815448A (en) * 2017-02-07 2017-06-09 长江水资源保护科学研究所 A kind of river attenuation type pollutant analogy method
CN111241758A (en) * 2020-01-10 2020-06-05 燕山大学里仁学院 Evaluation method based on transport and diffusion model of soluble pollutants in water environment
CN112685688A (en) * 2021-01-12 2021-04-20 长江水资源保护科学研究所 Long-distance water delivery open channel water quality simulation and prediction method and system
CN112861360A (en) * 2021-02-19 2021-05-28 河海大学 Maskyo flow calculation error correction method based on system response theory
CN113836758A (en) * 2021-07-23 2021-12-24 中国市政工程华北设计研究总院有限公司 Water quantity and water quality simulation calculation method for low-influence development facility
CN114139259A (en) * 2021-11-29 2022-03-04 深圳市水务科技发展有限公司 Method for constructing river channel water quality model

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106815448A (en) * 2017-02-07 2017-06-09 长江水资源保护科学研究所 A kind of river attenuation type pollutant analogy method
CN111241758A (en) * 2020-01-10 2020-06-05 燕山大学里仁学院 Evaluation method based on transport and diffusion model of soluble pollutants in water environment
CN112685688A (en) * 2021-01-12 2021-04-20 长江水资源保护科学研究所 Long-distance water delivery open channel water quality simulation and prediction method and system
CN112861360A (en) * 2021-02-19 2021-05-28 河海大学 Maskyo flow calculation error correction method based on system response theory
CN113836758A (en) * 2021-07-23 2021-12-24 中国市政工程华北设计研究总院有限公司 Water quantity and water quality simulation calculation method for low-influence development facility
CN114139259A (en) * 2021-11-29 2022-03-04 深圳市水务科技发展有限公司 Method for constructing river channel water quality model

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
SAID M. EASA, MENG, PHD, FCAE, FEIC, FCSCE: "New and improved four-parameter non-linear Muskingum model" *
张清寰等: "湘江捞刀河流域氨氮和总磷水环境容量计算" *
杨程炜;李订芳;刘德地;: "基于模拟优化模式的梯级水库群突发水污染事件应急调度研究" *
潘沛: "太湖丘陵地区农业小流域非点源污染研究" *
边凯旋;黄强;杨元园;刘登峰;刘东;韩知明;万新;: "河流突发水污染事件的水库应急调度流量研究――以珠江上游红水河为例" *

Also Published As

Publication number Publication date
CN115017727B (en) 2023-04-21

Similar Documents

Publication Publication Date Title
CN108664647B (en) Basin fine management system of integrated water environment model
CN109558973B (en) Water quality and water ecology integrated early warning system, control equipment and storage medium
CN111241758B (en) Evaluation method based on transport and diffusion model of soluble pollutants in water environment
CN113569438B (en) Urban flood model construction method based on multisource rainfall fusion and real-time correction
CN109711095B (en) A method of section local inflow is obtained based on hydrological model
CN114117848B (en) Basin water environment simulation prediction method and device based on multi-model coupling
CN105825230A (en) Soil moisture status prediction method and system based on DBN (Deep Belief Network) model
CN110119590B (en) Water quality model particle filtering assimilation method based on multi-source observation data
CN111985106A (en) Distributed hydrological model parameter calibration method based on multipoint parallel correction
CN110765676B (en) Watershed water quality simulation method based on stable flow field
CN112149314A (en) Multi-sand reservoir capacity silt flushing simulation method based on virtual reservoir capacity correction
CN114580316B (en) Small reservoir flood level forecasting method based on two-dimensional-zero-dimensional coupling model
CN113435630B (en) Basin hydrological forecasting method and system with self-adaptive runoff yield mode
CN112113146B (en) Synchronous self-adaptive check method for roughness coefficient and node water demand of water supply pipe network pipeline
CN111597732B (en) River network water flow numerical simulation method using water surface gradient of branch of a river point influence area
CN115796075A (en) River network hydrodynamic process simulation method based on branch point water quantity conservation
CN111104746B (en) River flood beach elevation determination method based on wavelet analysis
CN114881544A (en) Reservoir flow regulation and control method and device, electronic equipment and storage medium
CN110334851B (en) Hybrid cascade reservoir combined flood control optimal scheduling method considering branch flood storage civil application
CN113343601A (en) Dynamic simulation method for water level and pollutant migration of complex water system lake
CN113553782B (en) Downscaling method for forecasting wind speed
CN112541611B (en) Rain-raising agricultural ground source pollution emission prediction method and system
CN115017727A (en) Sewage collection simulation method based on Masjing root method
CN110847112B (en) River flood discharge early warning method based on hydraulics simulation
CN112651185A (en) Method for dealing with sudden water pollution event based on water quantity and water quality coupling model

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
GR01 Patent grant
GR01 Patent grant