^{1}

^{1}

^{1}

Taking Liyang City of Jiangsu Province as the research object, a water pollution accident diffusion simulation and analysis system, taking WebGIS as the core, was established by taking finite difference method to do the system dynamics water quality model solution. This paper introduces the construction process of the model, the integration of the model with GIS, the interactive application of the system and users and the visual dynamic expression of water pollution accidents. By combining the statistical analysis of the spatial and environmental database in the system, the scientific simulation of the temporal and spatial relationship of the development of water pollution accidents on the client side is realized and providing precise decision support for early warning and control of water pollution diffusion.

Water resources are not only the basic resources of social and economic development, but also the source of human life. However, in the past decade, due to the acceleration of industrialization and urbanization, sudden water pollution incidents have occurred frequently in China, seriously threatening the surface water ecosystem and urban water supply security. However, the water environment system is a complex dynamic evolution system, which continuously exchanges material and energy with outside world. By using the geographic information system (GIS) based on the geospatial database, this system not only completes the real-time query and analysis of water quality basics, but also realizes the visual expression of sudden water pollution accidents through its spatial display ability to illustrate the extent and scope of the accidents, making it easier to make a scientific quantitative assessment and prediction so as to achieve the role of water quality early warning and timely control of pollution.

The research on water pollution diffusion started early in developed countries. The STORET-COGENT (Storage and Retrieval Database System & Construction Grants Evaluation and Network Tracking System) [

In the mid-1980s, China began to attach importance to water pollution accidents. In order to study and explore the impact of water environmental safety accidents on social economy, the Songliao River Basin Management Committee developed a water quality model and water environment information management system [

In summary, under the new normal development of the Internet, the new subject Geographic Information System (GIS), which integrates geographic data processing, can integrate environmental management and ecological protection into the GIS system. By utilizing the efficient spatial analysis ability and powerful data display ability of GIS, it has injected new ideas into the traditional low-efficiency and high-cost water environmental protection work. However, the above water environment information management system still has many areas to be improved: it cannot simulate one-dimensional and two-dimensional water pollution accidents at the same time; the model has weak pertinence and strong limitations in application; the interface is complex or the model is obscure and difficult to popularize. Taking Liyang City of Jiangsu Province as an example, this paper chooses the finite difference water quality model which is most suitable for the hydrological characteristics of the basin, develops a surface water quality monitoring and water pollution diffusion analysis and simulation system suitable for county and municipal environmental protection applications, realizes the simulation visualization of the temporal and spatial distribution of sudden water pollution accidents in rivers and reservoirs and constitutes basic information , decision-making and human-computer interactions supporting the emergency treatment of sudden water pollution accidents scientifically, efficiently and accurately.

The main watershed of Liyang is the Shahe Reservoir of Tianmu Lake and Daxi Reservoir and its tributaries. The Daxi Reservoir is located in the upper reaches of the Daxi River in the hilly area of the southwest of Liyang City. The Shahe Reservoir is composed of the Dasha River originating from Tianlu Mountain and the Fushi River originating from Xiafang Village. Both belong to the Tianmu Lake watershed of the Nanhe River system in the western part of the Taihu Lake Basin and are closely related to the ecological stability of the Taihu Basin.

The water quality in this area mainly was influenced by agricultural pollution caused by farming and livestock farming, with ammonia nitrogen and P as the main elements. Because it is located in the subtropical monsoon climate with abundant precipitation, these pollutants are discharged into reservoirs from paddy fields with irrigation, resulting in water eutrophication.

Therefore, combined with the hydrological characteristics and pollution source types of the study area, the system chooses a suitable water environment quality prediction model for the watershed to quantify the sudden water events in the area.

Based on system dynamics, this system establishes a dynamic simulation model of sudden water pollution diffusion simulation, which mainly includes system dynamics hydrodynamic model and water quality model and integrates the final calculation results with GIS system to realize the visualized dynamic expression of water pollution diffusion simulation.

The one-dimensional water quality model is a relatively simple environmental water model. When rivers or lakes are polluted, this water pollution prediction model is often used when the difference between the concentration distribution of cross sections and the average concentration of cross sections is small. It mainly studies variations of pollutant concentration distribution along the river and the variations over time in each section. The one-dimensional model only considers the transportation and distribution characteristics of pollutant concentration and its diffusion parameters in the longitudinal direction, i.e., the flow direction (generally set as the x-axis). When pollutants are transported along the river, not only the flow but also the self-purification of the river should be considered. The corresponding water quality model [

∂ C ∂ t = E x ∂ 2 C ∂ x 2 − U x ∂ C ∂ x − K C (1)

Here, C denotes the concentration of pollutants in rivers, t denotes time, E_{x} denotes the diffusion coefficient in the X direction, U_{x} denotes the average discharge, and K denotes the decay rate coefficient of the pollutants.

Because the finite difference method allows for easy construction of the system dynamics model, the finite difference method is chosen to calculate the numerical difference of the model. Obtained by deduction and rearrangement:

C i j + 1 = C i j + ( E i C i + 1 j + C i − 1 j Δ x 2 + u i C i j Δ x ) Δ t − ( E i 2 C i j Δ x 2 + u i C i + 1 j Δ x + k i C i j ) Δ t (2)

In the formula, E_{i} is the longitudinal dispersion coefficient of section I; u_{i} is the average velocity of section I; k_{i} is the pollutant decay rate coefficient of section I; Δx are steps in the X direction; Δt are the time steps; i is the section I; and j is the time, j.

From the formula, the concentration value at time (j + 1), section I can be calculated iteratively from the concentration value at time j, section (i + 1), (i − 1) and other parameters (diffusion coefficient E_{i}, average velocity u_{i}, time steps Δt, position steps Δx and pollutant decay rate coefficient k_{i}).

When pollutants enter the water body, it is impossible to achieve uniform concentration mixing of the whole section in a short time. When there are significant differences in both the vertical and horizontal directions, the one-dimensional model cannot meet the needs of the system and the two-dimensional simulation model should be used [

∂ C ∂ t = E x ∂ 2 C ∂ x 2 + E y ∂ 2 C ∂ y 2 − U x ∂ C ∂ x − K C (3)

In the formula, the dispersion coefficients E_{x} and E_{y} are in the x and y directions, respectively. U_{x} is the velocity component in the X direction, and K is the decay rate coefficient of the pollutants. Since the essence of the system dynamics model is a set of differential equations, the finite difference solution of the two-dimensional river water quality model is as follows:

C i , j n + 1 = C i , j n + ( E i , j x C i + 1 , j n + E i , j x C i − 1 , j n Δ x 2 + E i , j y C i , j + 1 n + E i , j y C i , j − 1 n Δ y 2 + u i , j C i , j n Δ x ) Δ t − ( 2 E i , j x C i , j n Δ x 2 + 2 E i , j y C i , j n Δ y 2 + u i , j C i + 1 , j n Δ x + k i , j C i , j n ) Δ t (4)

Here, E i , j x is the longitudinal dispersion coefficient of the I and J sections, E i , j y is the transverse dispersion coefficient of the I and J sections, u_{i}_{,j} is the average flow velocity of the I and J sections, and k_{i}_{,j} is the decay rate coefficient of pollutants in the I and J sections.

In order to realize the space-time simulation and visualization of water pollution accidents, it is necessary to display the model calculation results in the browser with the help of webGIS so as to facilitate interaction with users. The research area is gridded in two-dimensional plane map by two programming languages: Java and JavaScript. Then the pollutant concentration in the grid is calculated by the above-mentioned water pollution diffusion analysis model. And then the grid is rendered dynamically according to the pollutant concentration and the display effect is enhanced to achieve accurate visualization. Specific implementation steps (

Detailed steps are as follows:

1) Obtaining the coordinates of river basin nodes by webGIS: divide the acquisition area into 10 m * 10 m grid on average and get the coordinates of each grid in the corresponding coordinate system.

The geometric attributes of x_{min} and y_{min} of the corresponding watershed are accessed by webGIG. The cell grid with 10 m step length is assigned by for cycle, stored in the array in turn and the maximum outer rectangle of the watershed is obtained at the same time.

2) By setting the initial parameters for the model calculation under the simulated state, the pollutant concentration values at each time in each grid are obtained.

According to the arithmetic provided by the water model, the concentration of pollutants in each grid can be obtained by iteration of the concentration of pollutants in the preceding moment near the top, left, bottom and right of the grid.

3) In order to enhance the visualization of pollution diffusion, graded color settings were used to convert the water pollution concentration values into RGB true color images.

4) According to the set unit time, the dynamic effect of water pollution diffusion is generated.

Through the statistical analysis of the system spatial basic information database and environmental basic information database, the water quality status of each section can be grasped macroscopically and the main pollutant sources, flow rate and direction of the corresponding watershed can be analyzed so as to obtain more accurate initial parameters and simulate the diffusion process of water pollution accidents.

In addition, by using GIS, we can get the geographic coordinates of the incident site according to the place names and locations provided by the eyewitnesses and make full use of remote sensing technology to obtain water environment and water pollution information as an integrated information resource for the establishment and inspection of water pollution model [

When the user visits the browser, by specifying the pollution point source in the corresponding River basin and setting the initial state parameters such as the average flow velocity, the pollutant concentration of the pollution source, the diffusivity coefficient in all directions and the attenuation coefficient of the pollutant, the user can start the analysis and then calculate the concentration of each point iteratively through the above model. The client parameter setting interface is shown in

Visualization of water pollution accident diffusion will simulate the results of calculation after inputting the corresponding parameters in different time periods and realize the dynamic expression of the spatial distribution of pollutants at different time in each area of the client through the map layer released by ArcGIS Sever so as to realize the simulation of the whole process of water pollution accident and determine its diffusion range and impact degree (

The analysis of pollutant concentration change includes the analysis of pollution source concentration and pollution reference point concentration. It can analyze the trend of concentration change with time at the designated point and display in broken line graph (as shown in

This system takes Liyang water area of Jiangsu Province as the research object, adopts B/S architecture, chooses the finite difference method which is more suitable for expressing its diffusion state according to the characteristics of the research area, couples the water quality model with the simulation and prediction model of water environment quality of the river basin and GIS, realizes the visual expression of one-dimensional and two-dimensional water quality mathematical models on the server side and can simulate space-time relationship of water pollution development of rivers and lakes respectively. Combining the attribute database of spatial geographic and environmental basic information and eyewitness feedback mechanism, we can analyze and summarize the initial parameters which are most in line with the real-time state of the study area, accurately simulate and control the water pollution diffusion and timely establish

the spatial relationship between the pollution source and the object of water environmental protection. The model provides decision support for the early warning and control of water pollution diffusion, effectively improving the efficiency of pollution control and making the integrated management of water environment more precise and quantitative.

At the end of this paper, I would first like to thank my tutor, Professor Du Jinglong, for his careful guidance and strict requirements, which enabled me to overcome the difficulties and finally complete the research task. Mr. Du’s noble personality, meticulous academic attitude, inspired me in the spirit of my work, gave me the motivation to continue to work hard on the road of scientific research.

In addition, I would also like to thank my team members, Li Wei and Wan Genyi, for their great help and encouragement during the study. We have jointly optimized the implementation of complex water quality model algorithms through the exchange discussion process and improved system operating efficiency.

Secondly, I would also like to thank the project fund: Changzhou Social Development Project “Study of County (City) Level Intelligent Environmental Information Platform―Taking Fuyang City as an Example” to provide data and financial support.

Finally, I am especially grateful to my parents, who taught me to be a human being, often spurred and encouraged me, and silently supported me to complete the research.

The authors declare no conflicts of interest regarding the publication of this paper.

Kuang, T.Q., Du, J.L. and Zhang, W.W. (2019) Simulation of Water Pollution Diffusion Based on WebGIS. Journal of Computer and Communications, 7, 28-38. https://doi.org/10.4236/jcc.2019.77004