Ocean is an important part of the global life support system and also a natural resource for human survival. Its huge water and good hydrodynamic conditions has become an ideal place for acceptance and purification of tail water. Coastal waters, as the most closely related to human activities in the ocean, are deeply affected by human activities. In recent years, as the world energy crisis intensifies, vigorously develop Liquefied Natural Gas (LNG) project in the coastal areas is becoming an important part of China's energy strategy. However, in the process of the production, the sea water is changed into the ocean after heat transfer as the heating medium, and has a certain degree of influence on the water environment in the coastal area . Therefore, it is very important to study the hydrodynamic characteristics of the coastal waters and the 2-D temperature distribution model of the relevant sea areas. It will make a reasonable planning of marine resources protection for the environmental protection departments and provide strong scientific support for marine engineering decision.
Yaohua Dong established a 2-D water flow, temperature and concentration transport numerical model of river, processed the determination of parameters and coefficient, selection of boundary definite condition, the treatment of shore boundary and moving boundary, and proved the feasibility and applicability of the model . The finite nodal method was based on the triangular grid and the sliding boundary condition was used to deal with the change of water level by Bin Luo . Zhenhai Li from China Institute of Water Resources and Hydropower Research used the windward format to deal with the convective term of the momentum equation, used the central difference scheme to deal with the diffusion term, used the control volume method to deal with the continuity equation and heat transfer equation. Then, the flow field and temperature field of discharge from power plant are solved . Haijie Wu predicted for the heated discharge from coastal power plant. He used dry wet unit hydraulic model to deal with the calculation of complicated boundaries, established a second-order Osher format water-temperature model, which provided a reliable technical data for marine environmental impact assessment . But the study above has not forecast the maximum possible influence range of temperature.
In this paper, a low temperature drainage problem from the LNG project was used as a research example, and hydrological conditions of summer is generally larger than that of winter. The tide current field in shallow sea area of Jiangsu coastal area was simulated through 2-D hydrodynamic model, and the time variation and spatial distribution characteristics of pollutant migration and transformation in the receiving sea area were simulated based on these hydrodynamic conditions. The heat of the 2-D convection diffusion model was established, and which was used to forecast the temperature field of temporal and spatial distribution characteristics considering the premise of water heat transfer effect. In the forecast, maximum influence range of the low temperature drainage was the focus, which meant that the value of the temperature envelope line selected the maximum value of the whole simulated period at one point in the position, so as to provide the basis for the maximum ecological environment security of the adjacent water area of the LNG project.
2. Hydrodynamic Simulation
A 2-D hydrodynamic model was used to simulate the tide current field in the shallow sea area of the radial sand ridges in the Jiangsu coast (see Figure 1). The range from south to north is 31˚58'N - 33˚14'N, the spacing is 152 km; the range from east to west is 120˚46'E - 121˚51'E, width across 102 km, and the total area of the calculation area is about 9850 km2.
Different combination of sequence boundary points on the tide level process was taken as the input condition of hydrology. Due to the location of the tide stations in the area was not conducive to the model boundary conditions in the study area, the tide data couldn’t be used as a tide of input data to simulate the moving boundary. Therefore, the ocean tide assimilation model (TPXO) was used to obtain the boundary of tide level process, which is one of the 7 ocean tide models in the world. In order to verify the accuracy of the model, we simulated the hydrodynamic conditions during April 24, 2008 to May 6, 2008, in which period contained more measured data. This method is mature, many literatures have discussed, so it would not repeat in this article  .
2.1. Hydrodynamic Model Validation
According to the measured data report of the relevant ocean hydrological observation, took the tide observation data to validate the hydrodynamic model, and compared the consistency between the calculation of tide level and tide duration and the measured value. The longitude and latitude of the tide observation stations are shown in Table 1. CW1-CW4 tide stations results shown in Figures 2-4.
Visible from above, the simulation results of tide level and tide duration are in good agreement with observation data. Based on the tide level and tide duration
Figure 1. Sea area terrain around LNG project.
Figure 2. Verification process of the mid tide on April 2008 24-25.
Figure 3. Verification process of the neap tide on April 2008 28-29
Table 1. The longitude and latitude of the tide observation stations.
verification above, it shows that the model could reflect the current situation of the hydrodynamic conditions well. Therefore, the grid layout, boundary conditions and calculation parameters could be used to simulate the hydrodynamic conditions in the predicted area well.
2.2. Hydrodynamic simulation results
Considering the seasonal difference of tidal current in the predicted area, the different flow field directly affects the migration and diffusion of pollutants, so as to influence the prediction results. Select summer and winter two seasons as a typical period of time to simulate the hydrodynamic conditions in the region. Flow field distribution of summer and winter are shown in Figure 5 and Figure 6.
2.3. Hydrodynamic Characteristics Analysis
The results show that, the main trend of the flow was similar to NE-SW during the flood tide, which is perpendicular to the coastline, and when the ebb tide is opposite. The regional flow rate is generally larger, the flow rate of the summer
Figure 4. Verification process of the spring tide on May 2008 5-6.
Figure 5. Summer flow field map.
Figure 6. Winter flow field map.
was significantly greater than that in winter, and the water velocity of the flood tide is higher than the ebb tide.
In the area near the shore, affected by the shoreline and topography, tidal current flows paralleled to the shore, and the local water flow is different from the main tidal current due to the change of coastline. Overall, during the flood tide, the trend in the northern coast area is southward, the trend in southern coast area is in the direction of the northeast, and when the ebb tide it is opposite. The velocity of water flow in the near shore water is less than the velocity of the water.
The drainage outlet is located on the artificial island which has a certain distance from the shore boundary, so the velocity and direction is less affected by the shore, and the water flow condition is relatively simple. Results shows that, the flow direction is approximately west at the sewage outfall during the flood tide, and the flow direction is east during the ebb tide. The average velocity at the outfall was 0.253 m/s in spring tide and 0.182 m/s in neap tide.
3. Chlorine and Water Temperature Distribution Simulation
3.1. 2-D Heat Diffusion Equation
Form in: T――receiving water temperature, ˚C; t――time, s;,――ver- tical and horizontal vertical average velocity, m/s; Kx, Ky――thermal diffusion coefficient of water body, m2/s; - source terms of water gas heat exchange; -source term temperature, ˚C; -source term flow, m3/s.
3.2. Definite Solution Condition
Form in: is boundary normal.
All initial conditions are 0.
3.3. Determination of Parameters and Source Strength
3.3.1. Diffusion Coefficient
The velocity field in natural sea area is generally not uniform. Diffusion coefficient is an important parameter which is not easy to be determined, and it often has the property of empirical coefficient. The diffusion coefficient formula used in this paper was as follows:
Form in: h――water depth, m; ――frictional velocity ( , j is hydralic gradient), m/s; ， ―― empirical coefficient, respectively taken as 4.0 and 0.5.
3.3.2. Water Surface Integrated Heat Dissipation Coefficient
The water surface comprehensive heat dissipation coefficient Ks was calculated by Gunnerberg formula of foreign experience formula, and the influence of wind and sea water temperature on thermal diffusion was considered:
Form in: t――water surface temperature.
3.3.3. Pollutant Source Concentration
According to the data, the original concentration of the drainage of a LNG receiving station project was 10.35 m3/s during operation, the water temperature was 3 degrees lower than the normal water.
3.4. Forecast Result
In this paper, we used the above model and its solution, simulated the full tide cycle under two typical hydrological conditions in winter and summer by considering the most adverse effects, and predicted the temperature increment and range. The equivalent envelope line of the water temperature drop value near the drains in winter and summer are shown in Figure 7 and Figure 8, and the spatial distribution statistical parameters of contour lines are shown in Table 2 and Table 3.
Figure 7. Equivalent envelope line of maximum temperature drop near drains in summer.
Figure 8. Equivalent envelope line of maximum temperature drop near drains in winter.
Table 2. Spatial distribution statistical parameters of water temperature drop near drains in summer.
Table 3. Spatial distribution statistical parameters of water temperature drop near drains in winter.
4. The Effect of Low Temperature Drainage on Environment
Based on the predicted results, the low temperature water formed an extension of the strip area around the drain under the action of tide current. And due to the local influence of artificial island, the pollutant zone extended along the direction of the shore boundary attaching to the bank near the drain. In the summer hydrological conditions, the velocity of tide current field is larger, and the diffusion condition is better, so the influence range of the low temperature water is larger than winter, but large temperature drop covered by a smaller area. The typical pollution load values are represented by the size of the pollution belt as follows:
1) In summer, the pollution zone near the drains that the maximum temperature drop greater than 1˚C was approximated to a band range of the length of 230 m, width of 100 m, and the area of the 0.023 km2; and the zone that greater than 0.1˚C was approximated to the length of 6900 m, width of 1400 m, and the area of 6.68 km2.
2) In winter, the pollution zone near the drains that the maximum temperature drop greater than 1˚C was approximated to a band range of the length of 220 m, width of 130 m, and the area of the 0.027 km2; and the zone that greater than 0.1˚C was approximated to the length of 7400 m, width of 1850 m, and the area of 9.73 km2.
In summary, the influence range of low temperature drainage was small, and pollution near the drain was more serious. In a full tide current, the average value of the pollution load in the water was much lower than the maximum value, because the equivalent line was selected as the maximum value at different time, which appeared only in an instant.
In addition, according to Pollution control standards for marine sewage disposal project (in China) and Seawater quality standard, the range allowed by the sea area is 3.0 km2, and the water temperature change cannot exceed 4˚C. So it can be concluded that the low temperature drainage of LNG project has a less impact on the water environment in the research sea area.
5. Conclusion and Prospect
1) Simulation results show that, the low temperature water formed an extension of the strip area around the drain under the action of tide current. And due to the local influence of artificial island, the pollutant zone extended along the direction of the shore boundary attaching to the bank near the drain. In the summer hydrological conditions, the velocity of tide current field is larger, and the diffusion condition is better, so the influence range of the low temperature water is larger than winter, but large temperature drop covered by a smaller area.
2) In summer, the pollution zone near the drains that the maximum temperature drop greater than 1˚C was approximated to a band range of the length of 230 m, width of 100 m, and the area of the 0.023 km2. In winter, the zone was approximated to a band range of the length of 220 m, width of 130 m, and the area of the 0.027 km2.
3) The predicted result was the maximum increment of pollutant concentration in adverse hydrological condition, that is, the concentration of pollutant was selected as the maximum concentration at different time in the whole simulation period. This maximum occurs only in an instant. Therefore, in a full tide current, the average value of the pollution load in the water was much lower than the maximum value. The simulation results also show that the maximum increment of the pollutant concentration in the project area was also relatively small, therefore the low temperature drainage of the project was relatively small impact on the ocean environment nearby.
1) This paper was based on the numerical simulation of 2-D model. Actual water flow has three dimensional characteristics, so 3-D model is more appropriate to reflect the actual flow of water. The future research may consider using the 3-D model to simulate the sea water transport characteristics and pollutants.
2) The influence of residual current and wave was not considered in this paper, so the simulation of the water flow couldn’t show the nonlinear characteristics of the actual coastal water. These factors should be considered in further research in the future.
3) The heat exchange mechanism between the low temperature water discharge and the natural water body is different from that of warm water discharge, and this paper has not carried out in-depth research. These factors can be considered in further research in the future in order to better simulate imitate the transport characteristics of temperature field after the cold water enters the natural water body.
Work partially supported by grant of the special funds from Ocean and Fisheries Bureau of Jiangsu Province (key technical research on the discharge of sewage in deep sea area).