شبيه سازي عددي و بررسي آزمايشگاهي انحلال كارست نمكي در مخازن سدها

از جمله مهمترين منابع آلوده كننده در محيط ­هاي آبي بويژه در مخازن سدها، سامانه ­هاي كارست نمكي مي­ باشد كه انحلال آن در آب مخزن باعث كاهش شديد كيفيت آب مخزن مي­ شود. در تحقيق حاضر براي بررسي رفتار كيفي مخزن سد در حضور سازند كارست نمكي، ابتدا بررسي ­هاي ميداني از مخزن انجام گرفت. سپس براي كاليبراسيون مدل عددي و تعيين ضريب انحلال سازند كارست نمكي، مدل آزمايشگاهي بصورت كانال ذوزنقه ­اي طراحي و ساخته شد. سپس مدل آزمايشگاهي با استفاده از مدل عددي MIKE شبيه ­سازي شد. با مقايسه نتايج مربوط به مدل آزمايشگاهي و مدل عددي، مقدار تقريبي ضريب انحلال سازند نمكي موجود در مخزن تعيين و مدل عددي كاليبره شد. در ادامه با تعميم ضريب انحلال بدست­ آمده به مخزن واقعي و شبيه ­سازي عددي مخزن، رفتار كيفي آن مورد تحليل و بررسي قرارگرفت. همچنين تحليل حساسيت پارامترهاي موثر بر فرآيند انحلال و توزيع شوري انجام شد و نحوه اثرگذاري و ميزان تأثير مقادير ضريب انحلال نمك، دبي، دما، تبخير و رطوبت نسبي بر توزيع شوري و شرايط كيفي مخزن در اثر انحلال سازند نمكي مورد مطالعه قرار گرفت. نتايج اين مطالعه مي­ تواند در طراحي و بهره ­برداري از مخازن داراي سازند نمكي مفيد باشد.

The presence of halite karst systems in dam reservoirs can be one of the crucial sources of pollutants for the water resources. Karstic regions are those with dissolving caves. If the dam is constructed in the regions with reactive rocks it is likely to be endangered by karst phenomenon, dissolution and corrosion of current waters. construction of dams and reservoirs in these formations are associated with high risk. A comprehensive and accurate understanding of the salt dissolution processes facilitates the description of salt karst formations dissolution process in dam reservoirs. This research was performed based on three bases: field studies, experimental investigation and numerical modeling. For determination of dissolution coefficient, the physical model of the reservoir was built in the form of the trapezoidal channel and samples of salt rock in the reservoir were used in the physical model of karst formation. This model was implemented in different conditions and the results of dissolution amount and salt concentration were measured for different time scales. At the next stage, the physical model was simulated in the numerical model for determining of dissolution coefficient. The numerical model ran with different dissolution coefficients and results were compared with experimental results. With this comparison and adaptation of results, the dissolution rate of halite karst determined and then used in the numerical simulation of the main reservoir. With simulation of main reservoir using dissolution coefficient obtained from laboratory study, the effect of salt karst formation on the reservoir’s water quality was analyzed. The influences of effective parameters such as dissolution rate, air temperature, discharge and relative humidity on the water quality of dam reservoir were studied along with carrying out a wide range of sensitivity numerical analysis. To this purpose, a versatile finite volume tool ‘MIKE’ was used. Using the experimental study and numerical simulation and implementation of results, the dissolution coefficient of salt karst formation was calculated about 0.2 cm/h. Numerical model results showed that early months of operation were not critical to the project and salt concentration in reservoir was in the range of 0.8 to 0.9 of allowable limit. But considering the changes occurring in the reservoir’s conditions, the process of salinity distribution increases and salt concentration increases with time. By comparing the salt concentration in reservoir and numerical model, it was found that the average error percentage of concentration values was about 13.65 percent and the numerical model showed the values lower than the reservoir in most times. The sensitivity analysis of parameters affecting the salinity distribution showed that these parameters have a significant impact on salinity distribution and salt concentration in the reservoir. It was found that with 10 percent changes in studied parameters, dissolution rate with 21%, air temperature with18%, discharge with 21.1% and relative humidity with 6.1% caused changes in salt concentration of reservoir. According to the results, it is evident that dissolution rate and discharge were the most important factors that influence the behavior of salinity distribution in the dam reservoir.