Indian Journal of Marine Sciences Vol. 43(7), July 2014, pp. Thermal pollution modelling of cooling water discharge into a closed creek system S. Buvaneshwari1*, Vijaya Ravichandran2 & B. V. Mudgal3 1 Indian Institute of Science, Department of Civil Engineering, Bangalore, India. 2 Coastal and Environmental Engineering Division, National Institute of Ocean Technology, Chennai, India. 3 Centre for water resources, Anna University, Chennai, India. *[Email: [email protected]] Received 26 august 2013; revised 16 October 2013 This paper is focused on the behaviour of the heated effluent discharged at an elevated temperature into the Ennore creek by North Chennai Thermal Power Station (NCTPS).Functioning of thermal power plants near tidal creeks and estuaries, due to the feasibility of intake and discharge of water for cooling poses serious environmental concern from heated water discharge due to the flow characteristics of the creek. Discharge of the heated effluent and its dispersion is a function of ebb and flood flow into the Ennore creek as the creek mouth remains closed for most part of the year. This study aims to investigate the dispersion characteristics of the heated effluent in the near-field and far-field under various discharge and dynamic conditions of the creek using calibrated hydrodynamic models and comparisonwith the field data measured during varying tidal conditions. [Keywords:Thermal Pollution; Hydrodynamics; CORMIX; MIKE21] Introduction Urban sprawl, industrialization and development, fuel the insatiable demand for energy. Thermal power contributes to the major portion of the power demand. It uses the fossil fuel to generate power. The heat that is generated by the power plants needs to be taken out of the system and is generally done through cooling waters. If the cooling water is recycled, then cooling towers are adopted, else in a once through cooling systems the heated water is safely discharged back into the water bodies adhering to the norms. Here heated water becomes the pollutant. Tamil Nadu Electricity board operates a thermal power station of 600 MW capacity at Ennore, North Chennai Thermal Power Station (NCTPS).The power plant uses sea water to meet its cooling water requirement and the condenser reject heated water which is about 8°C higher than the intake waters is let into the Ennore creek. The heated water discharge is about 27.5 m3/s. Ennore creek is a dynamic and ecologically sensitive water body receiving waste water from numerous sources including industrial and municipal effluents. In open coastal waters effluent discharged from outfalls undergoes rapid initial dilution, typically within a few hundred meters of the outfall, before reaching either a level of neutral buoyancy or the ocean surface. It undergoes a number of physical, chemical and biological processes in the water column before being transported out of the region or settling down into the sediments or being accumulated by biota. Mixing and dilution of the effluent occur due to turbulence that causes entrainment and diffusion.The source of this turbulence changes with distance from the discharge point. Near the diffuser, the primary source of turbulence is shear, induced by the discharge momentum and buoyancy. Local currents and local boundaries may modify this “self induced” turbulence. This region is often called near field or initial mixing region. Farther from the diffuser, this self-induced turbulence decays, and mixing is primarily due to turbulence naturally present in the ocean. This region is called far field1,2. A prerequisite for estimating the changes to the water environment is the understanding of the movement of the water mass as defined by the waves and tides, otherwise called wave and tidal hydrodynamics and current circulation patterns. Behaviour of plumes from ocean outfalls can be predicted using numerical models and assessment needs to be made of the three different stages, viz., initial dilution (near field modeling), transport (advection and dispersion modeling) and decay (fate) of pollutants (water quality modeling). INDIAN N J. MAR. SCI.,, VOL. 43, NO. 7, JULY 2014 The mixinng and transport process inn the creek iss highly h depen ndent on the tidal t ebbing and floodingg 3,4 conditions c .Shoreline chhanges along the Ennoree coast c and cloosure of the creek mouth due to sandd bar b formation n results in loow circulatioon of flow inn the t creek. Thhis study aim ms at modelingg the fate andd transport t of thermal disccharges into the Ennoree Creek C in the near field and a far field. The modell results r are vaalidated against real time measurements m s observed o in the field durring various seasons of a year y and different timess of the day given thee influence i of ambient a temp perature on th he advection-dispersion d chharacteristics of o the creek. Materials M an nd Methods Study S area Ennore E Creekk The southheast coast of o India is an a importantt stretch of coaastline, wheree many majorr rivers drainn into i the Bay of Bengal annd they are also a richer inn marine m fauna than the western coaast of India. Ennore E creek k (13°13'54.48" N, 80°19' 26.60" E) iss located l in the t northeastt coast of metropolitann Chennai C city y, Thiruvalluur district, Tamil T Nadu, India I as show wn in Fig. 1. N Fig. 1.Ennoore Creek-Enviroons Most of the t area conssists of alluviial tracts andd beach b dunes,, tidal flats and creek in n the easternn part. p Ennoree comprises of lagoonss, with saltt marshes m andd backwaterss, which aree submergedd under u water during d high tide and form an a arm of thee sea opening in i to the Bay of Bengal. The T total areaa of o the creek is 4 km2 annd is nearly 400 4 m wide. Creek C channeels connect to the Pulicatt Lake in thee north n and to the t Kortalaiyaarriver in the south. Ennore crreek was oncce a place wiith rich floraa and a fauna. It was not onlyy the nature’ss gift but alsoo a predominannt source for supporting th he traditionall fishermen f co ommunity seettled in this Creek. Thee creek c once fiilled with richh biodiversityy contributedd an a excellent green belt which w is now totally takenn out o by the untreated seewage from Royapuram m sewage outfall, unttreated and treated induustrial effluentts from Manali Industrial Belt, B which hoouses huge nuumber of cheemical indusstries. The naatural wealth of the creeek is given place to seewage channel and the bioloogical producctivity of the coast me down. Thee northern seection of the creek has com is the Kosasthalaiyar K r backwater which w is connnected to the Pulicat P Lake, with North Chennai Theermal Power Plant and the t satellite port of Ennnore. Wastew water enters the creeek through the Buckinggham canal inn the north an nd south of creek, c Kortalaiiyarriver andd Amullavoyyal canal inn the south. Ennore Theermal Poweer Station ETPS E withdraw ws cooling water from m the creek and disposess warm watter into seaa through marine m outfall. CTPS) Nortth Chennai Thhermal Powerr Station (NC unit disccharges the cooling c waterr into Buckinggham canal which w joins wiith Ennore Creek. C It has aopen a pre cooling channel having widthh up to 130 m for about 2.5 2 km, afterr flowing thrrough the exiisting warm water w channel for about 2km m, the warm water w travels a distance off 4.5 km to reach the Ennnore creek annd thereby miixing with thee creek waterr. The objectivve of such aan arrangemeent was to bring down thhe thermal polllution. The frequent closure of the Ennore E mouthh has resultedd in insufficcient tidal inflow and thus reducedd cooling waters for the thermal power p plants. Hence NCTP PS has opted d to dischargge the warm water w back into the crreek throughh the Buckinggham Canal, in an attemppt to maintaiin the water qu uantity in thee creek. However, it was found f that the warm water w would flow directly d back to t the intake structure wiith minimum m retention time, resultingg in the pow wer plant witthdrawing waarmer water for their operations, it has stoopped wing water from creek and now intake withdraw water foor cooling is ttaken inside Ennore E port. Ennore Creek Mouthh The currents in thhe coastal watters are influeenced by the seasonal circculation in thhe Bay of Beengal. Currents move alongg the coast towards t the North N M to Octoober and moove towards south from March from Noovember to Feebruary. The flood and ebbb tide is a significant s ccomponent of o the dischharge characteeristics in the Ennore creek k The dynamic Ennnore creek moouth influencees the current and tide ppattern insidee the creek, the currentss are generallly tidally inflluenced depennding on the mouth m condittion. As a parrt of the studdy the width of o the creek mouth was measured inn the month of o April and it was foundd that it was 72 m wide. BUVANESHWARI et al.: THERMAL POLLUTION MODELLING OF COOLING WATER DISCHARGE Field Investigation The temperature observation of NCTPS plant effluent is divided into two parts as shown in the Fig. 2. In Fig. 2 I indicates the observations inside the NCTPS plant of about 2km stretch channel and the measurements with respect to discharge, velocity profiles, depth, temperature are taken using Q-Liner instrument which works under acoustic Doppler principle. The region II (Fig. 2) indicates the observations of the temperature inside the creek for spring tide and neap tide and it is covered for both flood and ebb tide conditions using the YSI probe. Fig. 2. Layout of NCTPS Thermal Discharge NCTPS Effluent Dispersion- Choice of model Dispersion and mixing of pollutant concentrations in receiving waters is dependent on the ambient conditions of the receiving environment and discharge characteristics of the effluent. The ambient conditions are defined by the water body's geometry, currents, as well as its dynamic characteristics. The discharge conditions are a function of outfall geometry (diameter, height above bed, orientation) and flux characteristics (discharge rate, density, momentum and buoyancy)5, 6. Immediately after the effluent exits the diffuser into the ambient waters, the initial jet characteristics of momentum flux, buoyancy flux, and outfall geometry influences the jet trajectory and mixing. This region will be referred to as the "near-field", where outfall design can control the initial mixing characteristics through adjustment of design variables. In this region, the effluent jets are turbulent fluid zones in which vigorous mixing takes place. The USEPA length scale rule based empirical CORMIX model is particularly useful in the near field as turbulence modeling using deterministic hydrodynamic approaches is very complex7. Beyond the near field, ambient currents and density stratification gradually deflect and mix the plume into the predominant flow. As the turbulent plume travels further away from the source, the source characteristics become less important. The ambient environmental conditions will control trajectory and dilution of the turbulent plume through buoyant spreading motions and passive diffusion due to ambient turbulence. This region will be referred to here as the "far-field"8.Thus, the designer must keep in mind that minor geometric changes in the diffuser block may alter plume characteristics in the near field, while changes in the far field may need alterations in the significant geometric changes such as wider spacing between risers/ports9. CORMIX predictions in the near field are often benchmarked by some modellers, yet the model is limited in its ability to predict far field plume characteristics. Another model used for diffuser design is USEPA’s PLUMES (DOS version) or Visual Plumes (Windows version). For far field predictions, deterministic 2-D and 3-D hydrodynamic models like USEPA’s HSCTM2D and EFDC, Danish Hydraulic Institute’s MIKE21 and MIKE3 are considered to be more suitable. In nearshore coastal waters, where ambient conditions can be turbulent due to waves, currents and winds, CORMIX’s far field formulations are likely to be conservative. Thus, for large discharges, models such as MIKE21 should also be considered. For cooling water effluents from power plants and brine effluents from desalination plants, buoyancy flux is important10. As the effect of the relative density difference between the effluent discharge and ambient conditions in combination with the gravitational acceleration determines the tendency for the effluent flow to rise (i.e. positive buoyancy) or to fall (i.e. negative buoyancy). This study uses CORMIX3 as well as MIKE21. The modeling of NCTPS effluent dispersion into Ennore creek is done for the present discharge rate27.5 m3/sec by varying the creek mouth conditions. Indicates the present plant capacity (600MW) discharge that is let into the creek after travelling a distance of about 4.5 km inside the NCTPS plant and the hydrodynamic conditions comprise flooding and ebbing of spring and neap tides under creek mouth open and closed condition.The hydrodynamic dispersion INDIAN J. MAR. SCI., VOL. 43, NO. 7, JULY 2014 characteristics of the thermal effluent discharge by NCTPS at Ennore creek is modelled using CORMIX3 for near-field dispersion characteristics and MIKE21for the far-field dispersion characteristics. Model results are compared with the field observations to ensure the applicability of the hydrodynamic model to this problem. The temperature observations in the Ennore creek for the spring and neap tide are generated as temperature plots for the flood and ebb tide using the software surfer. Further the MIKE21 hydrodynamic model is used to predict the dispersion characteristics for the excessive discharge. The scenarios simulated for the heated effluent dispersion are given in the Table 1 Table 1.Simulation Scenarios of MIKE21 Scenarios Creek Mouth Creek Mouth Fig. 3. Measurements in the NCTPS Channel Bridges Table 2. Temperature Measurement Section Location Width of the Mouth (m) Discharge rate (m3/sec) Location Latitude Longitude Remarks 210 72 27.5 27.5 Location 1 80° 19' 40" 13° 15' 38" starting point Location 2 80° 19' 43" 13° 15' 47" Bridge 1 Location 3 80° 19' 48 " 13° 15' 50" Bridge 2 Location 4 80° 20' 03 " 13° 15' 39" Bridge 3 Location 5 80° 20' 02 " 13° 15' 25 " Bridge 4 Location 6 80° 20' 07" 13° 15' 24" Bridge 5 Location 7 80° 20' 06" 13° 15' 19" Diffuser point Results Field Observation inside NCTPS North Chennai thermal power station a 600MW capacity plant discharges 27.5 m3/sec of condenser reject water at about 5°C to 8°C above the ambient, travelling a distance of about 2km in a lined channel and then into a very shallow water body of about 2.5km distance inside the NCTPS plant and reaches the Ennore creek. The first set of measurements started from the NCTPS condenser reject into the channel. Since the channel sides were covered by bushes taking temperature measurements with the Q-Liner was difficult hence seven bridge sections as given in table 2 were selected to deploy the instrument to measure the discharge, depth, surface temperature and velocity profiles. The temperature was also measured with respect to depth using YSI probe and it was observed there is no thermal stratification along the depth of 2m. The results also indicated that the temperature drop was negligible over a distance of 4.5 km. Q-liner output Fig. 4. Q-Liner Output i. ii. Flow rate Velocity = 27.13 m3/sec = 1.2 to 0.9 m/sec BUV VANESHWARII et al.: THERM MAL POLLUTIO ON MODELLIN NG OF COOLIN NG WATER DIS SCHARGE iii. iv. Depth = 2.33 m Temperatu ure = 36°C The condenser rej eject temperatture dependss on o the intaake sea surfface temperaature. Thesee measurement m ts were taken in February 2011. 2 Field F Observvation inside the t Creek-Floood Tide In order to t measure th he extent of dispersion d inn the t creek, temperatures were w observeed using YSII probe p and thhe locations are a noted usinng hand heldd GPS. G The obbservations were w made onn spring andd neap n tide day ys during the month of Appril and Mayy in i each day temperature t w measuredd in the creekk was starting from m creek mouth h up to Buckiingham canall during d flood tide t and ebb tide. t The tempeerature disperrsion plots foor the presentt discharge d conndition is gen nerated using surfer for thee tidal t flow conndition in the creek. Towarrds Buckinngham canal Fig. 6. Therm mal Plume Trace during Flood Tiide The temperature oobservation inn the Ennore creek is coverred for duratioon of 3hours before the start of flood tiide. The therrmal effluentt dispersed inn the creek iss travelling a distance of about a 500 m from the disccharge point to the creekk mouth, as it is pushed by the creek water to the open o sea as shhown in the Fiig. 7 Fig. 5.Floood Tide Temperature Plots- Aprril 2011Spring annd Neap The tempeerature observ vation in the Ennore E creekk is i covered foor 3hours beffore the start of ebb tide, when w almostt the whole creek is fillled with thee seawater as the maximum m flood tide condition iss reached. r It was observed that t the tempeerature of thee NCTPS N effluuent entering the creek varied 8°C too 9°C 9 above thhe ambient inn April whichh is the mostt critical c periodd as the ambiient temperatuure is high inn summer wheen compared to other seassons and thee thermal t plum me travels thee entire creekk distance off about a 1km before disp persing to the t ambientt conditions. c This T is observved for the sppring tide off April A and Maay flood conddition as shoown in Fig. 6 and a Fig. 7. The temperrature plot of o neap tidee indicateslikel i ly, rise in tem mperature of about 0.5°C to t 1°C which h may be due to the solar radiation r andd the t time of the day duuring summeer when thee measurement m t was made. It needs to be b noted thatt daytime d of summer is thee most critical period forr temperature t r rise in the ambient a as well w as in thee effluent e water. Figg. 7. Ebb Tide Teemperature Plotss- April 2011Sprring and Neap The thermal plum me movement towards the sea s at o 5°C abovee the a high temperature of about to ambientt and this phenomenon n is predom minant during April. A The sollar radiation is i observed too play a vital role in raising the tempperature of about 0.5°C too 1°C during noon which is spread oveer the entire creek c causing negligible variation inn the temperaature profile. INDIAN N J. MAR. SCI.,, VOL. 43, NO. 7, JULY 2014 Thermal Dispersion (27.5 m3/s) for Creek Mouth M Closed Condition C usiing Mike21 Creek mouth opening to sea Fig g. 8. Thermal Plu ume Trace during g Ebb Tide Itindicatess the thermall plume beinng pushed byy the t creek water w from Buckingham m canal andd Kortalaiyar K R River, therebyy the plume gets g deflectedd towards t the sea. s The obseervations werre also madee along a the traajectory to trrace out the temperaturee variation v alonng the thermaal plume trajecctory. Figg. 11. Thermal D Dispersion from Low Tide to Higgh Tide Thermal Dispersion (27.5 m3/s) for Creek Mouth M Open Condition C using Mike21 Thermal Disppersion (27.55 m3/s) of NC CPTS Effluentt during d Flood d Tide using Cormix3 C Fig. 9..Plots of Presentt Discharge durinng Flood Tide It is founnd that the teemperature at a the end off near-field n is 4°C above the t ambient and it is nott completely c d dispersed even n after a distaance travel off about a 625 m towards the Buckingham B canal. c Fiig. 12.Thermal D Dispersion from Low Tide to Higgh tide Fig. 10.Plots of Presen nt Discharge durring Ebb Tide It is founnd that the teemperature at a the end off near-field n is 5°C above the t ambient and it is nott completely c dispersed d eveen after a traavel of aboutt 750 7 m towaards sea as the thermal plume getss deflected d by the t creek watter. Indiccate the thhermal plum me being puushed towardss sea and traavelling towaards port andd the effect iss felt to a disttance of 500m m at 5°C abovve the ambientt and enters the Buckinghham canal att 4°C above thhe ambient annd thermal plu ume width off both tide is leess when com mpare to moutth close condiition Comparrison of Modeel Results withh Field Data The Thermal plum me simulated by the CORM MIX3 and MIK KE21for the present dischharge under creek BUVANESHWARI et al.: THERMAL POLLUTION MODELLING OF COOLING WATER DISCHARGE mouth closed condition is compared with the surfer plots generated by the field data for various tides.The model results and surfer plots are compared at the entrance of Buckingham canal and at the creek mouth point opening to sea. As MIKE21 gives a time series data the average temperature of flood and ebb tide is taken for comparison. In CORMIX3 the temperature at the downstream coinciding with the points of comparison is taken into account.The temperature at the Buckingham canal during flood tide is 4°C above the ambient and at the creek mouth 5°C above the ambient during ebb tide, are matching with the model results and the temperature plots. Thus the results showed a good argument with the field data emphasizing the proper simulation of the model. plume is felt to a distance of 500 m towards Ennore port during ebb and the thermal plume travels a distance of 1.5 km and reaches the Buckingham canal at 4°C above the ambient. Acknowledgements I would like to express my earnest gratitude to Dr. Purnima Jalihal, Dr. Prince Prakash Jebakumar and Mr. Rajkumar Mr. Raju Abraham, Mr. Prasad Dudhgaonkar, Scientist National Institute of Ocean Technology, Mr. E. Arun Babu, Centre for Water Resources, for their valuable suggestions. References Grace R A, Marine outfall systems: planning, design, 1. 2. 3. Conclusions The near-field dispersion characteristics of the NCTPS heated effluent discharge during flood tide simulated using CORMIX3 surface buoyant discharge of about 27.5 m3/sec indicates a plume travel distance of about 625 m towards Buckingham canal and the temperature at thenear field region end is 4°C above the ambient. This indicates the temperature is not dispersed in the Ennore creek of 2 km stretch. Further the plume enters the Buckingham canal at a temperature rise of about 3°C above the ambient and during ebb tide condition plume travels a distance of about 750 m towards sea and the temperature at the near field region end is 5.5°C above the ambient. As the NCTPS outfall is close to the creek mouth, the flow getting reversed during ebb tide condition deflects the plume at a shorter distance and reaches open sea at an elevated temperature. The far-field dispersion behaviour of the heated effluent is simulated using MIKE21 Advection Dispersion Model for the creek mouth closed and open condition for a discharge rate of 27.5m3/sec. The results showed a good variation of temperature with respect to time, a clear picture of the thermal plume movement with respect to flood and ebb condition and it is found that the effect of thermal 4. 5. 6. 7. 8. 9. 10. and construction. 1978: Prentice-Hall. Fischer H B, Mixing in Inland and Coastal Waters. 1979: Academic Press. Webb T and Tomlinson R, Design Procedures for Effluent Discharge to Estuaries during Ebb Tide.Journal of Environmental Engineering, 1992. 118(3):338-362. Lewis R, Dispersion in Estuaries and Coastal Waters. 1997: Wiley. Jirka G and Akar P, Hydrodynamic Classification of Submerged Multiport‐Diffuser Discharges.Journal of Hydraulic Engineering, 1991. 117(9):1113-1128. Jirka G and Domeker R, Hydrodynamic Classification of Submerged Single‐Port Discharges.Journal of Hydraulic Engineering, 1991a. 117(9):1095-1112. USEPA., Technical Bulletin No. 624, User's Guide For The Cornell Mixing Zone Expert System (CORMIX), USEPA, Center for Exposure Assessment Modeling, Athen, Georgia, USA. 1992. Doneker R and Jirka G, Expert Systems for Mixing‐Zone Analysis and Design of Pollutant Discharges.Journal of Water Resources Planning and Management, 1991b. 117(6):679-697. Bleninger T, User’s Manual for CORHYD: An Internal Diffuser Hydraulics Model, Institute of Hydromechanik, Universtat Karlsruhe. 2005. Schreiner S P, Krebs T A, Strebel D E, and Brindley A, Testing the CORMIX model using thermal plume data from fourMaryland power plants.Environmental Modelling & Software, 2002. 17:321-331.
© Copyright 2017 ExploreDoc