نوع مقاله : مقاله پژوهشی
نویسندگان
1 کارشناسی ارشد آبیاری و زهکشی، گروه علوم و مهندسی آب، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران
2 استادیار گروه علوم و مهندسی آب، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران*(مسئول مکاتبات)
3 دانشیار گروه علوم و مهندسی آب، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران
چکیده
کلیدواژهها
موضوعات
علوم و تکنولوژی محیط زیست، دورهبیست و یکم، شماره هفت، مهر ماه 98
مدلسازی کیفیت آب رودخانه با استفاده از مدل QUAL2Kw
(مطالعه موردی: رودخانه شاهرود)
رزیتا آریایی نژاد[1]
مهدی سرائی تبریزی[2]*
M.sarai@srbiau.ac.ir
حسین بابازاده[3]
تاریخ دریافت: 29/06/96 |
تاریخ پذیرش:08/09/96 |
چکیده
زمینه و هدف: مدلسازی کیفیت آب رودخانهها میتواند بهعنوان یکی از مؤثرترین ابزارها در جهت مدیریت کیفیت آب رودخانهها و کاهش تبعات مسایل محیط زیستی ناشی از ورود آلایندهها به آن مورد استفاده قرار گیرد. هدف از مقاله حاضر استفاده از مدل معتبر QUAL2Kw.5.1 در راستای مدلسازی کیفیت آب در رودخانه شاهرود میباشد.
روش بررسی: در مقاله حاضر ازمدل QUAL2Kw.5.1 برای مدلسازی هفت پارامتر کیفیت آب شامل اکسیژن محلول، اکسیژنخواهی بیوشیمیایی، قلیاییت، کل مواد جامد محلول، هدایت الکتریکی، فسفر کل و نیتروژن کل در چهار دوره زمانی در رودخانه شاهرود استفاده شده است. به منظور واسنجی مدل از دادههای آبان 1386 و تیر 1387 و برای صحتسنجی مدل از دادههای شهریور 1387 و مهر 1387 استفاده شده است. واسنجی خودکار ضرایب مدل با استفاده از الگوریتم ژنتیک موجود در مدل انجام شده است. به منظور مقایسه نتایج شبیهسازی شده مدل با دادههای مشاهداتی از دو پارامتر آماری شامل ضریب تبیین و میانگین قدرمطلق خطا استفاده شده است.
یافتهها: مهمترین ضرایب واسنجی شده در مدل مربوط به مواد جامد معلق، اکسیژنخواهی بیوشیمیایی، نیتروژن و فسفر میباشد. مدل در شبیهسازی پارامترهای pH و هدایت هیدرولیکی با مقادیر میانگین قدرمطلق خطای 19/0 و 89/163 در مرحله صحتسنجی بهترتیب بیشترین و کمترین دقت را از خود نشان داد. به طور میانگین حداقل و حداکثر اکسیژن محلول اندازهگیری شده در رودخانه شاهرود در شهریور 1387 و مهر 1387 و بهترتیب برابر 93/6 و 99/9 میلیگرم بر لیتر اندازهگیری شده است. همچنین بیشترین و کمترین دقت مدل در شبیهسازی این پارامتر مربوط به تیر و مهر 1387 با میانگین قدرمطلق خطای برابر 86/0 و 29/1 میباشد. همچنین نتایج نشان داده که مدلسازی هیدرولیکی دقیق از تغییرات پارامترهای هیدرولیکی رودخانه و جریان در طول رودخانه تأثیر زیادی در مدلسازی کیفیت آب رودخانه دارد.
بحث و نتیجهگیری: نتایج مقاله حاضر دقت مناسب مدل QUAL2Kw.5.1را در شبیهسازی پارامترهای کیفیت آب رودخانه شاهرود نشان میدهد. از طرف دیگر دقت شبیهسازی هر پارامتر با توجه به مقدار تغییرات آن در طول رودخانه متغیر میباشد. به نحویکه هر چه تغییرات یک پارامتر در طول رودخانه و در بازههای مختلف کمتر باشد، دقت مدل در شبیهسازی این پارامتر بیشتر خواهد بود.
واژههای کلیدی:رودخانه شاهرود، مدل QUAL2Kw.5.1مدلسازی کیفیت آب، صحتسنجی، واسنجی.
|
Modeling Water Quality of Rivers Using QUAL2Kw Model (Case Study: Shahroud River)
Rozita Aryaee Nezhad[4]
Mahdi Sarai Tabrizi[5]*
Hossein Babazadeh[6]
Accepted:2017.09.13 |
Received:2017.11.29
|
Abstract
Background and Objective: Modeling water quality of rivers can be used as one of the most effective tools for water quality management in rivers and reducing the environmental impacts of entering pollutants. The purpose of this paper is to use the valid QUAL2Kw 5.1 model to model water quality in Shahrood River.
Method: In this paper, seven parameters of water quality have been used including dissolved oxygen (DO), biochemical demand oxygen (BOD), pH, total dissolved solids (TDS), total phosphorus, and total nitrogen four times in Shahrood River. Data from October of 2007 and July of 2008 were used to calibrate and data from September and October of 2008 were used to verify the model. Auto-calibration of model coefficients was done using genetic algorithm of the model. In order to compare simulated results with the observed data, determination coefficient and mean absolute error were used.
Findings: The most important calibration coefficients of the model were related to TSS, BOD, total nitrogen and phosphorus. This model in simulation of pH and EC with mean absolute error of 0.19 and 163.89 during verification stage showed the most and the least accuracy, respectively. On average the minimum and maximum DO were measured 6.93 and 9.99 mg/L in September and October of 2008 respectively in Shahroud River. Also the highest and lowest accuracy of the model in simulating these parameters were related to July and October of 2008 with mean absolute error of 0.86 and 1.29, respectively. In addition the results showed accurate hydraulic modeling of hydraulic parameters changes of the river along the river had a great influence on modeling of the water river quality.
Discussion and Conclusion: The results of this paper show the accuracy of the QUAL2Kw model in simulating water quality parameters of Shahrood River. On the other hand, the accuracy of the simulation of each parameter varies with the amount of its variation along the river so that the less the changes in a parameter along the river and at different intervals, the higher the accuracy of the model in simulating this parameter will be.
Keywords: Shahrood River, QUAL2Kw Model, Water Quality Modeling, Verification, Calibration
مقدمه
ورود انواع مختلف آلایندههای کشاورزی، صنعتی و پسابهای شهری به رودخانهها که یکی از منابع اصلی تأمین آب آشامیدنی، صنعتی و کشاورزی میباشند، باعث ایجاد بحران محیطزیستی در عصر حاضر شده است. واحدهای صنعتی، کشاورزی و مسکونی واقع در حوضه رودخانه از نظر مصرف آب به رودخانه وابسته میباشند. زهکشهای کشاورزی، پسابهای صنایع و خانگی و غیره وارد رودخانه شده و آلودگی رودخانه ها را همراه با مصارف آب تشدیدکرده است (1).
شکری و همکاران (1393) به شبیهسازی کیفی آمونیوم و نیترات در طول رودخانه گرگر با استفاده از مدل QUAL2KW پرداختند. نتایج حاصل از صحتسنجی ایشان نشان داد که مدل تطابق خوبی با واقعیت دارد و با توجه به مطالعات به عمل آمده بیشتر منابع آلوده کننده رودخانه گرگر را میتوان فاضلابهای روستاهای مسیر، پسابها و فاضلابهای ماهیسراها و حوضچههای پرورش ماهی و زهکشهای کشاورزی دانست(2). خداممحمدی و همکاران (1395) به ارزیابی توان خودپالایی و نقش اکسیژن محلول در کیفیت آب رودخانه کر (مطالعه موردی: پایین دست سد درودزن تا دریاچه طشک- بختگان) به کمک مدل کیفی QUAL2K پرداختند. نتایج تحقیق ایشان نشان داد که وجود پلها و تأسیسات آببند مانند بند امیر، پل رحمت آباد و بند حسنآباد میتواند باعث افزایش تلاطم و افزایش هوادهی و قدرت خودپالایی رودخانه کر و در نتیجه افزایش میزان اکسیژن محلول در رودخانه شود (3).
حسینی و همکاران (1395) با بررسی کارایی مدل QUAL2Kw در خودپالایی رودخانه (مطالعه موردی رودخانه کارون در بازه زرگان-کوت امیر) دریافتند که این مدل کارایی خوبی برای بررسی خودپالایی رودخانه دارد. در این تحقیق پارامترهایEC، NO3، TSS، BOD، DO،pH و دما توسط مدل شبیهسازی و واسنجی شده است. خروجی مدل به طور کلی حکایت از عدم خودپالایی رودخانه در پارامترهای pH وBOD داشته و خودپالایی کم در پارامترهای TSS و NO3با توجه به بار آلودگی ورودی در این بازه از رودخانه داشته است (4). حسینی و همکاران (1396) به بررسی تغییرات توان خودپالایی رودخانه کارون در سالهای 87 و 92 با استفاده از مدل QUAL2KW در محدوده شهر اهوازرودخانه کارون پرداختند. نتایج تحقیق ایشان نشاندهنده عدم خودپالایی پارامتر EC در رودخانه کارون برای سال 87 و عدم خودپالایی پارامتر BOD در سال 92 بوده است. همچنین بهترین شبیهسازی مدل برای پارامتر pH بدست آمد و در رتبههای بعدی بهترتیب پارامتر BOD و EC قرار گرفتند (5).
عاشقمعلا و همکاران (2015) در تحقیقی دو رودخانهی سبزکوه در چهارمحالوبختیاری و رودخانه قشلاق در کردستان را با مدل QUAL2KW شبیهسازی کردند. نتیجه این شبیهسازی نشان داد که رعایت استانداردهای تخلیه پساب، شرایط کیفی مناسب رودخانه قشلاق را به خوبی حفظ میکند. اما در رودخانه سبزکوه به دلیل میزان توان خودپالایی و تعداد زیاد منابع آلاینده موجود، رعایت استانداردهای تخلیه پساب، نه تنها کمکی به حفظ کیفیت آب نمیکند، بلکه تاثیر بسیار منفی بر کیفیت آب رودخانه خواهد گذاشت. سپس، با استفاده از شبیه سازی کیفی رودخانه، حد مجاز مناسب هر رودخانه بر اساس توان خودپالایی رودخانه و تعداد و نوع منابع آلاینده موجود برآورد شد و نتایج نشان داد که در تدوین استانداردهای دقیق و کارآمد باید به توان خودپالایی رودخانه، تعداد منابع آلاینده حاشیه رودخانه و میزان بارگذاری مواد زاید و دیگر شرایط متفاوت رودخانهها نیز توجه نمود (6). مهراسبی و همکاران (2015) کیفیت آب رودخانه کینهورز را با استفاده از مدل QUAL2KW شبیهسازی کردند. نتایج تحقیق ایشان نشان داد که نیتروژن کل و فسفر کل بهترتیب باید 76٪ و 93٪ کاهش یابد تا به هدف کیفیت آب مطلوب در رودخانه دست یافت (7).
ذوالفقاریپور و همکاران (2016) به منظور مدیریت کیفی رودخانه و حداقلسازی هزینههای حفظ محیط زیست و تامین معیارهای کیفیت آب به مدلسازی با استفاده از نرمافزار QUAL2Kw پرداختند. ایشان به منظور لحاظ کردن نقش تخلیهکنندگان بار آلودگی در تصمیمگیری و افزایش قابلیت اجرا برای برنامههای پیشنهادی تخصیص بار آلودگی در رودخانه از رویکرد گزینش اجتماعی استفاده کردند. برای هر تخلیهکننده، سناریوهای مختلف تصفیه فاضلاب مشخص و گزینههای ترکیبی تصفیه تعریف شدند. سپس با استفاده از مدل شبیه سازی کیفی QUAL2Kw، برای هر یک از گزینههای تصفیه، مقادیر جریمه تخطی از استاندارد کیفی محاسبه شد. نتایج حاصل از کاربرد مدل پیشنهادی در یک مطالعه موردی که قسمتی از رودخانه زرجوب در شمال ایران بود، نشان داد که تدوین سیاستهای تخصیص در کاهش بار آلودگی رودخانهها کارآیی مناسبی دارد (8). محققانی در سال 2017 با مدل کردن کیفیت آب رودخانه یامیونا در هند به کمک مدلQUAL2Kw نشان دادند که پارامترهایTN ، CBODf،DO و TC نسبت به جریان بالادست و کیفیت و موقعیت ورودی منبع آلاینده بسیار حساس است (9).
هدف اصلی پژوهش حاضر شبیهسازی پارامترهای کیفیت آب شامل هفت پارامتر اکسیژن محلول[7]، اکسیژنخواهی بیوشیمیایی[8]، قلیاییت[9]، کل مواد جامد محلول[10]، هدایت الکتریکی[11]، فسفر کل[12] و نیتروژن کل[13] در چهار دوره زمانی (آبان 1386، تیر 1387، شهریور 1387 و مهر 1387) در رودخانه شاهرود با استفاده از مدل QUAL2Kw 5.1 میباشد.
روش بررسی
معرفی رودخانه مورد مطالعه
حوضه آبریز رودخانه شاهرود بخشی از حوضه سفیدرود بزرگ بوده و با مشخصات طول جغرافیایی 30-49 الی 10-59 و عرض جغرافیایی 7-36 الی 45-36 قرار دارد که فاصله بین حد شرقی و غربی آن حدود 160 کیلومتر و عرض آن قریب به 25 کیلومتر میباشد. رودخانه شاهرود مهمترین و پر آبترین رودخانه استان قزوین است. رودخانههای الموترود و طالقانرود در محل روستای شیرکوه به هم پیوسته و رودخانه شاهرود را تشکیل میدهند. شکل 1 موقعیت رودخانه شاهرود را به همراه ایستگاههای نمونهبرداری نشان میدهد.
شکل 1- موقعیت ایستگاههای نمونه برداری آب
Figure 1- Location of water sampling stations
بازهبندی رودخانه شاهرود بر اساس 10 ایستگاههای نمونهبرداری شامل 4 ایستگاه روی شاخه اصلی رودخانه شاهرود و 6 ایستگاه رویدو شاخه ورودی به این رودخانه (2 ایستگاه روی الموترود و 4 ایستگاه روی طالقانرود) انجام شده است. موقعیت ایستگاههای نمونهبرداری در شکل 1 و مشخصات بازههای رودخانه شاهرود در جدول 1 نشان داده شده است. مطابق با جدول 1، طول کل مسیر مورد بررسی در رودخانه شاهرود برابر 13/189 کیلومتر مییاشد.
جدول 1- مشخصات بازههای رودخانه شاهرود مورد استفاده در مدل کیفی
Table 1- Characteristics of Shahrood River range used in qualitative model
ردیف |
رودخانه |
کد ایستگاه |
عرض جغرافیایی |
طول جغرافیایی |
فاصله از پاییندست(کیلومتر) |
طول بازه (کیلومتر) |
ارتفاع (متر) |
شیب بازه |
۱ |
طالقان رود |
SS1 |
18/36 |
88/50 |
13/189 |
- |
۱۹۴۰ |
- |
۲ |
طالقان رود |
SS2 |
17/36 |
75/50 |
59/176 |
54/12 |
۱۷۹۶ |
09/0 |
۳ |
طالقان رود |
SS3 |
20/36 |
63/50 |
10/164 |
49/12 |
۱۶۹۳ |
12/0 |
۴ |
الموت رود |
SS4 |
40/36 |
68/50 |
48/129 |
62/34 |
۱۶۳۶ |
61/0 |
۵ |
الموت رود |
SS5 |
39/36 |
39/50 |
49/102 |
99/26 |
۱۰۷۸ |
05/0 |
۶ |
طالقان رود |
SS6 |
39/36 |
39/50 |
98/101 |
50/0 |
۱۰۵۸ |
03/0 |
۷ |
شاهرود |
SS7 |
40/36 |
17/50 |
46/85 |
52/16 |
۹۲۰ |
12/0 |
۸ |
شاهرود |
SS8 |
61/36 |
75/49 |
50/26 |
96/58 |
۴۸۶ |
14/0 |
۹ |
شاهرود |
SS9 |
62/36 |
52/49 |
89/0 |
61/25 |
۳۲۶ |
16/0 |
۱۰ |
شاهرود |
SS10 |
62/36 |
51/49 |
00/0 |
89/0 |
۲۷۹ |
02/0 |
نمونهبرداری از رودخانه شاهرود طی چهار دوره در آبانماه 1386، تیرماه 1387، شهریورماه 1387 و مهر ماه 1387 توسط سازمان حفاظت محیط زیست انجام شده است. زمانهای نمونهبرداری با توجه به تغییرات دبی رودخانه و نوع منبع آلاینده و زمانهای تاثیرپذیری آنها انتخاب شده است. اکسیژن محلول، اکسیژن
خواهی بیوشیمیایی[14]، قلیاییت، کل مواد جامد محلول، هدایت الکتریکی، کدورت، فسفر کل و نیتروژن کل از مهمترین پارامترهای اندازهگیری شده در نمونهها میباشد. جدول 2 لیست دادههای مورد استفاده و مشخصات آماری آنها را نشان میدهد.
جدول 2- خلاصه آماری دادههای رودخانه شاهرود (سازمان حفاظت محیط زیست، 1389)
Table 2- Summary of Statistical data of Shahrood River (Environmental Protection Agency, 2010)
تاریخ نمونهبرداری |
آماره |
Q (m3/s) |
T (0C) |
DO (mg/l) |
BOD (mg/l) |
pH |
EC (dS/m) |
TS (mg/l) |
Turbidity (NTU) |
TP (mg/l) |
TN (mg/l) |
آبان 1386 |
مینیمم |
88/0 |
3/6 |
1/6 |
0/9 |
7/7 |
۰.۳۹ |
۳۴۴ |
۱۰ |
02/0 |
65/0 |
ماکزیمم |
20/9 |
8/9 |
7/9 |
0/38 |
2/8 |
۰.۸۷ |
۱۴۵۰ |
۸۸۰ |
05/0 |
04/2 |
|
میانگین |
92/3 |
5/8 |
6/7 |
3/18 |
9/7 |
۰.۶۸ |
۷۷۵ |
۳۱۳ |
03/0 |
25/1 |
|
انحراف معیار |
67/2 |
2/1 |
2/1 |
7/8 |
2/0 |
۰.۲۰ |
۴۱۸ |
۳۲۲ |
01/0 |
45/0 |
|
تیر 1387 |
مینیمم |
54/0 |
3/13 |
4/6 |
5/18 |
0/8 |
۰.۶۰ |
۳۸۰ |
۲ |
04/0 |
69/0 |
ماکزیمم |
10/5 |
3/27 |
9/7 |
0/84 |
5/8 |
۱.۳۴ |
840 |
۱۸۰ |
06/0 |
34/1 |
|
میانگین |
70/2 |
5/19 |
4/7 |
5/59 |
2/8 |
92/0 |
۶۱۴ |
۴۶ |
05/0 |
04/1 |
|
انحراف معیار |
61/1 |
1/5 |
5/0 |
9/28 |
2/0 |
27/0 |
۲۱۲ |
۵۹ |
01/0 |
26/0 |
|
شهریور 1387 |
مینیمم |
21/0 |
0/12 |
0/6 |
5/11 |
8/7 |
61/0 |
۳۸۵ |
۲ |
05/0 |
36/0 |
ماکزیمم |
76/4 |
0/27 |
8/7 |
5/22 |
2/8 |
30/1 |
824 |
90 |
13/0 |
05/1 |
|
میانگین |
84/1 |
1/20 |
9/6 |
8/18 |
9/7 |
00/1 |
۶۴۶ |
۳۹ |
۰.۰۷ |
62/0 |
|
انحراف معیار |
39/1 |
8/4 |
6/0 |
5/3 |
2/0 |
26/0 |
۱۷۷ |
۴۱ |
03/0 |
21/0 |
|
مهر 1387 |
مینیمم |
83/2 |
5/13 |
8/7 |
4/2 |
9/7 |
46/0 |
۲۲۷ |
۳ |
89/0 |
13/1 |
ماکزیمم |
19/9 |
5/23 |
3/11 |
8/42 |
4/8 |
74/0 |
۴۲۸ |
92 |
32/1 |
38/8 |
|
میانگین |
29/5 |
0/19 |
0/10 |
8/9 |
1/8 |
62/0 |
۳۳۲ |
۴۱ |
07/1 |
37/3 |
|
انحراف معیار |
34/2 |
1/4 |
0/1 |
1/13 |
1/0 |
10/0 |
۷۴ |
۴۰ |
17/0 |
14/2 |
* در جدول 2، واحد درجه حرارت آب درجه سانتیگراد، دبی مترمکعب بر ثانیه، اکسیژن محلول و اکسیژنخواهی بیوشیمیایی میلیگرم بر لیتر و هدایت الکتریکی دسی زیمنس بر متر میباشد.
معادله هندرسون– هاسلبالخ [15] برای محاسبه pH یک محلول با استفاده از ثابت تفکیک اسیدی بر اساس رابطه (1) استفاده میشود.
(1)
در رابطه (1)، [HA] غلظت مولی اسید نامحلول ضعیف،[A⁻]
غلظت باز مزدوج این اسید و pKaثابت تفکیک اسیدی میباشد.
مدلسازی کیفیت آب رودخانه شاهرود با استفاده از نرمافزار QUAL2Kw 5.1
در مقاله حاضر مدلسازی کیفیت آب رودخانه شاهرود با استفاده از نرمافزار QUAL2Kw 5.1، در پنج مرحله کلی انجام شده است.
این مراحل در شکل 2 به صورت فلوچارت نشان داده شده است.
شکل 2- فلوچارت مراحل مدلسازی کیفیت آب رودخانه با استفاده از مدل QUAL2Kw
Figure 2- Flowchart of river water quality modeling process using QUAL2Kw model
در مقاله حاضر برای واسنجی مدل از دادههای آبان 1386 و تیر 1387 و برای صحتسنجی مدل از دادههای شهریور 1387 و مهر 1387 استفاده شده است. واسنجی خودکار ضرایب مدل با استفاده از الگوریتم ژنتیک موجود در مدل انجام شده است. به منظور مقایسه نتایج شبیهسازی شده مدل با دادههای مشاهداتی از دو پارامتر آماری شامل ضریب تبیین (R2) و میانگین قدرمطلق خطا (MAE) مطابق روابط 1 و 2 استفاده شده است.
(2) |
|
(3) |
که در روابط (2) و (3)،، ، ، بهترتیب مقادیر مشاهداتی، شبیهسازی شده، میانگین مشاهداتیها و
میانگین شبیهسازی شدهها در ایستگاه iام و N تعداد دادهها میباشد.
یافتهها
در نرمافزار QUAL2Kw.5.1 تعداد زیادی ضریب و پارامتر
قابلیت انتخاب جهت واسنجی شدن را دارند. مقادیر تعدادی از مهمترین پارامترهای واسنجی شده در مدل مقاله حاضر با استفاده از روش واسنجی خودکار توسط الگوریتم ژنتیک در جدول3 ارایه شده است.
جدول 3- مقادیر پارامترهای واسنجی شده در مدل QUAL2Kw.5.1 با استفاده از الگوریتم ژنتیک
Table 3 - The values of calibrated parameters in the QUAL2Kw model using the genetic algorithm
ردیف |
دستهبندی پارمتر |
ضریب |
واحد |
حداقل |
حداکثر |
مقدار واسنجیشده |
1 |
جامدات معلق غیرآلی |
سرعت تهنشینی |
متر بر روز |
0 |
2 |
06/0 |
2 |
اکسیژنخواهی بیوشیمیایی کربنی |
نرخ هیدرولیز |
یک بر روز |
0 |
5 |
93/1 |
3 |
نرخ اکسیداسیون |
یک بر روز |
0 |
5 |
54/0 |
|
4 |
نیتروژن آلی |
نرخ هیدرولیز |
یک بر روز |
0 |
5 |
84/0 |
5 |
نرخ تهنشینی |
متر بر روز |
0 |
2 |
25/0 |
|
6 |
آمونیوم |
نرخ نیتریفیکاسیون |
یک بر روز |
0 |
10 |
15/2 |
7 |
نیترات |
نرخ دنیتریفیکاسیون |
یک بر روز |
0 |
2 |
03/1 |
8 |
فسفر آلی |
نرخ هیدرولیز |
یک بر روز |
0 |
5 |
44/3 |
9 |
فسفر غیرآلی |
نرخ تهنشینی |
متر بر روز |
0 |
2 |
014/0 |
جدول 3 نشان میدهد که مهمترین ضرایب واسنجی شده در مدل مربوط به مواد جامد معلق، اکسیژنخواهی بیوشیمیایی، نیتروژن و فسفر میباشد. از جمله مهمترین پارامترهایی که اثر زیادی بر محیط زیست سالم در یک رودخانه دارد، میتوان به دمای آب، اکسیژن محلول، اکسیژنخواهی بیوشیمیایی و نیتروژن کل اشاره کرد که در ادامه نمودارهای مقادیر مشاهداتی و شبیهسازی شده این پارامترها ارایه شده است.شکل 3 مقادیر مشاهداتی و شبیهسازی شده پارامتر دمای آب را در چهار دوره اندازهگیری شده نشان میدهد.
شکل 3- مقادیر مشاهداتی و شبیهسازی شده دمای آب در چهار دوره واسنجی و صحتسنجی
Figure 3 - Observed and Simulated Water Temperature in Four Calibration and Validation Periods
دمای آب یکی از مهمترین پارامترهایی است که روی سایر پارامترهای کیفیت آب رودخانه تأثیرگذار است. به طور معمول هر چه دمای آب بالاتر باشد، سرعت واکنشهای شیمیایی در آب بیشتر میشود. شکل 3 نشان میدهد که به طور میانگین کمترین و بیشترین دمای آب اندازهگیری شده در رودخانه شاهرود مربوط به آبان 1386 و شهریور 1387 و بهترتیب برابر 71/8 و 56/19 درجه سانتیگراد میباشد. همچنین با توجه به میانگین قدرمطلق خطا (MAE)، بیشترین و کمترین دقت شبیهسازی مدل مربوط به آبان 1386 و مهر1387 و بهترتیب برابر 15/1 و 16/4 میباشد. شکل 4 مقادیر مشاهداتی و شبیهسازی شده پارامتر اکسیژن محلول را در چهار دوره مدلسازی شده نشان میدهد.
شکل 4- مقادیر مشاهداتی و شبیهسازی شده اکسیژن محلول در چهار دوره واسنجی و صحتسنجی
Figure 4 - Observed and Simulated Dissolved Oxygen (DO) Values in Four Calibration and Validation Periods
پارامتر اکسیژن محلول آب با توجه به اینکه برای هر نوع واکنش شیمیایی و همچنین حیات موجودات زنده در آب ضروری است، یکی از پارامترهای بسیار مهم در بین سایر پارامترهای کیفیت آب میباشد. غلظت اکسیژن محلول در آب به طور مستقیم به برخی از سایر پارامترهای کیفی از جمله دمای آب وابسته است. با توجه به شکل 4، به طور میانگین حداقل و حداکثر اکسیژن محلول اندازهگیری شده در رودخانه شاهرود در شهریور 1387 و مهر 1387 و بهترتیب برابر 93/6 و 99/9 میلیگرم بر لیتر اندازهگیری شده است. همچنین بیشترین و کمترین دقت مدل در شبیهسازی پارامتر اکسیژن محلول مربوط به تیر 1387 و مهر 1387 با میانگین قدرمطلق خطای برابر 86/0 و 29/1 میباشد. شکل 5 مقادیر مشاهداتی و شبیهسازی شده پارامتر اکسیژنخواهی بیوشیمیایی را در چهار دوره مدلسازی شده نشان میدهد.
شکل 5- مقادیر مشاهداتی و شبیهسازی شده اکسیژنخواهی بیوشیمیایی در چهار دوره مدلسازی شده
Figure 5 - Observed and simulated biochemical oxygen demand (BOD) Values in Four Calibration and Validation Periods
غلظت پارامتر اکسیژنخواهی بیوشیمیایی (BOD) در آب به طور مستقیم به میزان فعالیت میکروارگانیسمهای موجود در آب وابسته است به طوریکه هر چه فعالیتهای موجودات زنده و واکنشهای شیمیایی در آب بیشتر باشد، مقدار مصرف اکسیژن محلول بیشتر و اکسیژنخواهی بیوشیمیایی بیشتر میشود. شکل 5 نشان میدهد که به طور میانگین بیشترین و کمترین مقادیر اکسیژنخواهی بیوشیمیایی در تیر 1387 و مهر 1387 و بهترتیب برابر 78/52 و 80/9 میلیگرم بر لیتر اندازهگیری شده است. مدل بیشترین و کمترین دقت شبیهسازی اکسیژنخواهی بیوشیمیایی را بهترتیب در شهریور 1387 و تیر 1387 و با میانگین قدرمطلق خطای برابر 20/5 و77/36 داشته است. شکل 6 مقادیر مشاهداتی و شبیهسازی شده پارامتر نیتروژن کل را در چهار دوره مدلسازی شده نشان میدهد.
شکل 6- مقادیر مشاهداتی و شبیهسازی شده نیتروژن کل در چهار دوره مدلسازی شده
Figure 6 - Observed and simulated total nitrogen Values in Four Calibration and Validation Periods
پارامتر نیتروژن کل از مجموع چهار فرم نیتروژن شامل نیتروژن آلی، نیتروزن آمونیاکی، نیتریت و نیترات بدست میآید. با توجه به شکل 6 بیشترین و کمترین غلظت نیتروژن کل بهترتیب در مهر 1387 و شهریور 1387 برابر 37/3 و 65/0 میلیگرم بر لیتر اندازهگیری شده است. همچنین بیشترین و کمترین دقت مدل در شبیهسازی نیتروژن کل در تیر 1387 و مهر 1387 با میانگین قدرمطلق خطای برابر 59/0 و 23/5 بدست آمده است. نتایج آماری مراحل واسنجی و صحتسنجی مدل QUAL2Kw.5.1 برای پارامترهای مختلف در جدول 4 ارایه شده است.
جدول 4- نتایج آماری مراحل واسنجی و صحتسنجی مدل QUAL2Kw
Table 4- The statistical results of the calibration and verification procedures
of the QUAL2Kw model
ردیف |
زمان نمونهبرداری |
آماره |
|
پارامترهای مدلسازی شده |
||||||
T |
pH |
EC |
TN |
TP |
CBODu |
DO |
||||
1 |
آبان 1386(واسنجی) |
MAE |
15/1 |
43/0 |
140 |
92/0 |
0 |
19/6 |
04/1 |
|
ضریب همبستگی |
12/0 |
3/0 |
23/0 |
05/0 |
45/0 |
03/0 |
13/0 |
|||
2 |
تیر 1387(واسنجی) |
MAE |
99/3 |
2/0 |
33/233 |
59/0 |
0 |
77/36 |
86/0 |
|
ضریب همبستگی |
33/0 |
12/0 |
12/0 |
32/0 |
31/0 |
33/0 |
14/0 |
|||
3 |
شهریور 1387(صحتسنجی) |
MAE |
86/3 |
23/0 |
67/246 |
59/0 |
06/0 |
20/5 |
14/1 |
|
ضریب همبستگی |
01/0 |
03/0 |
18/0 |
45/0 |
04/0 |
01/0 |
04/0 |
|||
4 |
مهر 1387(صحتسنجی) |
MAE |
16/4 |
14/0 |
11/81 |
23/5 |
07/1 |
21/7 |
29/1 |
|
ضریب همبستگی |
81/0 |
26/0 |
27/0 |
61/0 |
11/0 |
23/0 |
23/0 |
|||
دادههای واسنجی |
MAE |
57/2 |
32/0 |
67/186 |
75/0 |
0 |
48/21 |
96/0 |
||
ضریب همبستگی |
23/0 |
21/0 |
18/0 |
19/0 |
38/0 |
18/0 |
13/0 |
|||
دادههای صحتسنجی |
MAE |
02/4 |
19/0 |
89/163 |
91/2 |
565/0 |
18/4 |
12/3 |
||
ضریب همبستگی |
41/0 |
15/0 |
23/0 |
53/0 |
075/0 |
12/0 |
14/0 |
مقدار میانگین قدرمطلق خطا هر چه به صفر نزدیکتر و مقدار ضریب همبستگی هر چه به یک نزدیکتر باشند حاکی از دقت بیشتر نتایج مدل به دادههای اندازهگیری شده میباشد. جدول 4 نشان میدهد که بیشترین همبستگی برای پارامتر دمای آب در مهر 1387 و برابر 81/0 بدست آمده است. همچنین نتایج جدول 4 نشان میدهد که دقت مدل در شبیهسازی مقادیر هدایت الکتریکی دقت کمتری نسبت به سایر پارامترها داشته است.
بحث و نتیجهگیری
در مقاله حاضر مدل QUAL2Kw.5.1دقت بالایی در
شبیهسازی پارامتر pH از خود نشان داد که با نتایج حسینی و همکاران (1396) مطابقت دارد. همچنین نتایج تحقیق حاضد نشان داد که مهمترین ضرایب واسنجی شده در مدل مربوط به مواد جامد معلق، اکسیژنخواهی بیوشیمیایی، نیتروژن و فسفر میباشد که با نتایج شارما و همکاران (2017) همخوانی دارد. نتایج مقاله حاضر دقت مناسب مدل QUAL2Kw.5.1را در شبیهسازی پارامترهای کیفیت آب رودخانه شاهرود نشان میدهد که شکری و همکاران (1393) چنین نتیجهای را در مورد شبیهسازی کیفی آمونیوم و نیترات در طول رودخانه گرگر با استفاده ازمدل QUAL2Kw.5.1عنوان کرده بودند. دقت شبیهسازی هر پارامتر با توجه به مقدار تغییرات آن در طول رودخانه متغیر میباشد. به نحویکه هر چه تغییرات یک پارامتر در طول رودخانه و در بازههای مختلف کمتر باشد، دقت مدل در شبیهسازی این پارامتر بیشتر خواهد بود. داشتن دادههای دقیق از تغییرات پارامترهای هیدرولیکی رودخانه و جریان در طول رودخانه تأثیر زیادی در مدلسازی هیدرولیکی و بالتبع در مدلسازی کیفیت آب رودخانه دارد. همچنین دقت مدل در مرحله صحتسنجی به نوع دادههای مورد استفاده در مرحله واسنجی، مشخصات الگوریتم ژنتیک مورد استفاده در مرحله واسنجی تعداد ضرایب واسنجی شده، و تغییرات داههای مرحله صحتسنجی بستگی دارد.
تشکر و قدردانی
دادهها و اطلاعات مورد استفاده در این مقاله از سازمان حفاظت محیطزیست تهیه شده است که به این وسیله تشکر و قدردانی میشود.
منابع
[1]- کارشناسی ارشد آبیاری و زهکشی، گروه علوم و مهندسی آب، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران
[2]- استادیار گروه علوم و مهندسی آب، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران*(مسئول مکاتبات)
[3]- دانشیار گروه علوم و مهندسی آب، دانشگاه آزاد اسلامی واحد علوم و تحقیقات تهران
[4]- MSc. Student Department of Water Engineering, Science and Reserch Branch, Islamic Azad University, Tehran, Iran.
[5]- Assistant Professor, Department of Water Engineering, Science and Reserch Branch, Islamic Azad University, Tehran, Iran.* (Corresponding Author)
[6]- Associate Professor, Department of Water Engineering, Science and Reserch Branch, Islamic Azad University, Tehran, Iran.
[7]- Dissolved Oxygen (DO)
[8]- Biochemical Oxygen Demaon (BOD)
[9]- Potential of Hydrogen (pH)
[10]- Total Dissolved Solids (TDS)
[11]- Electrical Conductivity(EC)
[12]- Total Phosphorus (TP)
[13]- Total Nitrogen (TN)
[14]-Biochemical Oxygen Demand
[15]- Henderson–Hasselbalch equation
Sharma, D., Kansal, A., Pelletier, G., 2017. Water quality modeling for urban reach of Yamuna River, India (1999–2009), using QUAL2Kw. Applied Water Science; 7(3): 1535–155