نوع مقاله : مقاله پژوهشی
نویسندگان
1 استاد دانشکده محیط زیست دانشگاه تهران .
2 دانشیار دانشکده محیط زیست دانشگاه تهران
3 دانشیار گروه مهندسی آب دانشکدة کشاورزی، دانشگاه بوعلی سینا همدان
4 دانشجوی کارشناسی ارشد پردیس ارس دانشگاه تهران*(مسئول مکاتبات).
چکیده
کلیدواژهها
فصلنامه انسان و محیط زیست، شماره 21، تابستان 91
پیش بینی تغییرات پارامتر رسانایی الکتریکی در آب زیر زمینی شهر تهران
با کمک شبکه عصبیمصنوعی
ناصر مهردادی [1]
غلام رضا نبی بید هندی 1
اکبر باغوند[2]
حمیدزارعابیانه[3]
پویان عباسی مائده[4]*
جهت بررسی کیفیت آب زیر زمینی تهران با توجه به برداشت 10 سال گذشته از پارامتر های هیدروشیمیایی مربوط به 71 نقطه مختلف از شهر تهران سه مدل مختلف شبکه عصبی مصنوعی با تعداد پارامتر های مختلف ورودی و خروجی پارامتر رسانایی الکتریکی تعریف گردید. مشاهده می شود که به جهت تخمین پارامتر رسانایی الکتریکی با سعی و خطای فراوان تابع محرک تانژانت با تابع آموزش مومنتم دارای خطای کمی خواهند بود. با کمتر کردن مقدار پارامتر های ورودی و ثابت نگه داشتن تابع آموزش به مومنتوم و تابع محرک به تانژانت و تغییرات در مقدار پردازنده های داخلی شبکه عصبی میزان خطا کاسته خواهد شد و در نهایت سه مدل با یک لایه پنهان و تابع آموزش مومنتوم و تابع محرک تانژانت ساخته شد.حداکثر خطای نتایج به دست آمده نشان دهنده میزان بیشینه ضریب تعیین برابر 0.986می باشد که مربوط به مدل های 1 و 3 با تعداد پارامتر بیشتر می باشد. همچنین با توجه به شبکه عصبی ساخته شده در یک لایه کمترین میزان کمینه ی نرمال ریشه ی میانگین مربعات خطا برابر 0.110 در مدل های 1 و 3 خواهد بود، با توجه به تعداد ورودی های کمتر مدل شماره 2 و خطاهای با تقریب بسیار نزدیک به این دو مدل با میزان بیشینه ضریب تعیین برابر 96/0 و میزان کمینه ی نرمال ریشه ی میانگین مربعات خطا برابر 176/0 تقریب بسیار مناسب و قابل قبول خواهد بود و از میزان ورودی و آزمایشات جهت اندازه گیری پارامتر ورودی کاسته خواهد شد. در بحث تاثیر گذاری پارامتر ها در پیش بینی با شبکه عصبی مشخص می شود که میزان اهمیت دو یون سولفات و کلراید با اهمیت تر از پارامتر سدیم می باشند.
کلمات کلیدی: شبکه عصبی مصنوعی، تهران، آب زیر زمینی، رسانایی الکتریکی، پیش بینی.
پیش بینی کیفیت آب زیرزمینی برای مدیریت تأمین و بهره برداری از منابع آب، امری لازم است. نوسان های کیفیت آب زیرزمینی، اثریافته از عوامل زیادی است که قطعی نبودن آن ها، به دلیل تعدد متغیرهای موثر بر کیفیت آب حایز اهمیت می باشد. نامتجانس بودن محیط، سبب پیچیدگی پیش بینی کیفیت و تلاش پژوهشگران در بررسی تغییرات کیفیت مخازن شده است (1). پارامتر رسانایی الکتریکی (EC) جزو پارامترهای اصلی در پایش کیفیت آب به لحاظ شرب و کشاورزی می باشد. این پارامتر ارتباط مستقیم با میزان شوری آب، میزان جذب سدیم و میزان کیفیت شرب آب را دارد، به همین جهت در این تحقیق جهت پیش بینی مورد بررسی قرار گرفته است (2، 3، 4). یکی از رویکردهای مناسب در مطالعه رفتار آب های زیرزمینی، استفاده از مدل های رایانه ای است بنابراین، شناخت سازوکار تغییرات سطح کیفیت در طول زمان و پیش بینی آن با استفاده از روند حاکم برای اطلاع از وضعیت سفره و میزان آب قابل دسترسی، ضروری است (5). در این راستا، مدل های شبکه عصبی مصنوعی، به دلیل نبود درک فیزیکی از ماهیت مسئله، تنها براســـاس آموزش قــادر به مدل سازی رفتار دینامیکی فرایند غیرخطّی هستند. ویژگی آموزش سبب انعطاف پذیری شبکه ی عصبی مصنوعی در برابر خطاهای ناخواسته و بی تأثیر شدن آن ها در روند کار است(6،7). شبکه های عصبی مصنوعی با قابلیّت اعتماد بیشتر زمانی و سازگاری با تغییرات پیش بینی نشده، جایگزین مناسب مدل های تجربی و رگرسیونی برای پیش بینی رفتار منابع آب به شمار می روند.کاربرد شبکه هـــای عصبـــی نه تنها در پیش بینی وضعیت کیفی بلکه حتی در پیش بینی وضعیت سطوح آب های زیر زمینی و میزان احجام آن ها به کار می رود و مدل سازی های فراوانی در این زمینه ها صورت گرفته است (8-11). در رابطه با پیش بینی با کمک شبکه عصبی نیز مهردادی و همکاران در سال 2012 به پیش بینی کردن پارامتر مواد معلق محلول (TDS) توسط شبکه عصبی در تصفیه خانه فاضلاب فجر در جنوب کشور پرداخته اند (4). در بررسی های مشابه به جهت پیش بینی پارامتر نیترات با مدل سازی مشابه این تحقیق زارع ابیانه و همکاران دارای تجارب موفقیت آمیز بوده اند (16). همچنین در تحقیق مشابهی مدل سازی به جهت پیش بینی سطح آب زیر زمینی دشت ملایر صورت پذیرفته است که نتایج قابل قبول با مدل سازی مشابه حاصل شده است (17).
معرفی محل مورد مطالعه :
استان تهران به مرکزیت شهر تهران، با وسعتی حدود ۱۲۹۸۱ کیلومتر مربع بین ۳۴ تا 36.5 درجه عرض شمالی و ۵۰ تا ۵۳ درجه طول شرقی واقع شدهاست. تهران بزرگترین شهر استان و پایتخت کشور ایران است که موقعیت آن در شکل (1) نمایش داده شده است. جمعیت آن ۸،۴۲۹،۸۰۷ نفر است و هجدهمین شهر پرجمعیت جهان به شمار میآید. این استان همراه توابع خود (استان تهران)، جمعیتی برابر ۱۳،۴۲۲،۳۶۶ نفر و مساحتی برابر ۱۸،۸۱۴ کیلومتر مربع دارد. ارتفاع شهر در بلندترین نقاط شمال به ۲۰۰۰ متر و در جنوبیترین نقاط به ۱۰۵۰ متر از سطح دریا میرسد. تهران از شمال به نواحی کوهستانی و از جنوب به نواحی کویری منتهی شده در نتیجه در جنوب و شمال دارای آب و هوایی متفاوت است. نواحی شمالی از آب و هوای سرد و خشک و نواحی جنوبی از آب و هوای گرم و خشک برخوردارند.ساختار اداری ایران در تهران متمرکز شدهاست. طی یک دوره 30 ساله متوسط میزان بارندگی محدوده مورد بررسی 230 میلی متر گزارش گردیده است.وجود رودخانههای همیشگی مانند رودخانه کرج، رودخانه جاجرود، رود لار، حبله رود، رود شور یا ابهررود و طالقانرود موجب شده تا استان تهران از لحاظ منابع آب کمبودی نداشته باشد. بیشترین رودخانههای استان از کوه های البرز سرچشمه میگیرند. در استان تهران قناتهای زیادی وجود داشته که در گذشته نه چندان دور در تأمین آب مورد نیاز مناطق شهری و روستایی سهم بسزایی داشتهاند. ولی امروزه با استفاده از امکانات آب لولهکشی که از سدهایی چون سد امیرکبیر، لتیان و لار تأمین می شود، آب قناتها و چشمهها فقط برای کشاورزی و آبیاری مصرف می شود و فقط بعضی چشمهها، به ویژه چشمههای آب معدنی که بیشتر در شمال شرقی استان متمرکزند، اهمیت سابق خود را حفظ کردهاند. مهمترین این چشمهها عبارتاند از: چشمه اعلا دماوند، چشمه قلعه دختر، چشمه آب علی هراز، چشمه وله در گچسر، چشمه شاه دشت کرج، چشمه علی در شهرری، چشمه تیزآب، چشمه گله گیله(12، 13) .
شکل 1 -موقعیت شهر تهران در کشور ایران(12)
مواد و روش ها
پارامترهایی نظیر درجه حرارت، pH، هدایت الکتریکی(EC) و اکسیژن محلول(DO) در محل نمونه برداری توسط دستگاه های قابل حمل مورد ارزیـابی قرار گرفته و سایر پارامترها در آزمایشـگاه آنالیز گردید. اندازه گیری نیترات، نتیریت، سولفات و فلوئور با استفاده از دستگاه هک انجام شده که به ترتیب از روش های شماره 8039، 8507، 8051 و 8029 استفاده شده است. کلیه کاتیونها توسط روش 3005 EPA- با دستگاه پلاسمای جفت شونده القایی اندازه گیری شده است جهت اندازه گیری آنیون های کربنـات و بی کربنـات از روش اسـتاندارد مـتد با شمـاره 4500 استفاده شـــده و اندازه گیری آنیون کلرور به روش آرگنتومتریک (اسـتاندارد مـتد با شمـاره2330 ) صورت گرفت(14،15،2). مقادیر مختلف بیشینه، کمینه، میانگین و انحراف از معیار با کمک نرم افزار SPSS 19 در جدول 1 نمایش داده شده است.
جدول 1 - میانگین و مقدیر کمترین و بیشترین پارامتر ها در محدوده مطالعات
پارامتر |
تعداد |
کمینه |
بیشینه |
میانگین |
انحراف از معیار |
|
Statistic |
Statistic |
Statistic |
Statistic |
Std. Error |
Statistic |
|
EC |
1652 |
285 |
9630 |
96/1465 |
42/33 |
58/1358 |
TDS |
1652 |
164 |
6652 |
95/917 |
94/21 |
03/892 |
SO4 |
1652 |
51/0 |
54/56 |
76/5 |
18/0 |
61/7 |
Cl |
1652 |
17/0 |
60/45 |
7/4 |
14/0 |
93/5 |
HCO3 |
1652 |
05/1 |
75/12 |
24/4 |
52/0 |
13/2 |
TH |
1652 |
35 |
2472 |
27/394 |
53/8 |
78/346 |
SAR |
1652 |
234/0 |
39/28 |
35/3 |
07/0 |
13/3 |
K |
1652 |
01/0 |
31/0 |
036/0 |
00/0 |
03/0 |
Na |
1652 |
33/0 |
5/58 |
89/6 |
19/0 |
11/8 |
Mg |
1652 |
.16/0 |
84/29 |
95/2 |
08/0 |
41/3 |
Ca |
1652 |
4/0 |
9/31 |
93/4 |
09/0 |
88/3 |
جهت تحلیل و بررسی ارتباط میان آلاینده ها بر ایجاد برقـــراری رابطه براساس ضریب همبستگی پیرسون با نرم افزارSPSS 19 پرداخته می شود. با مطالعه ضرایب همبستگی به وجــود آمده مشخص می شــود که میــزان تاثیر پذیری هر پارامتر نسبت به پارامتر های دیگر چقدر است. در انتها نیز به بررسی مدل های پیشنهادی با شبکه عصبی طراحی شده جهت پیش بینی پارامتر رسانایی الکتریکی پرداخته شد. در مرحله اول با توجه به جدول 3 همبستگی های به وجود آمده از داده های آماری تهیه شده از سال 1381 تا 1390 در فصول خشک و تر مربوط به چاه های نمونه گیری شده به تعداد 71 حلقه در سطــح شهر تهــران بیشترین ارتباط های معنی دار با پارامتر رسانایی الکتریکی در پارامتر های جدول 3 که همبستگی به روش پیرسون در نرم افزار SPSS 19 به دست آمده است دیده می شوند. به جهت پیش بینی با کمترین خطا جهت هر کدام از مدل ها ساختار یکسان برای شبکه عصبی پیشنهاد شده است که دارای کمترین خطاها باشند و تغییرات در تعداد پردازنده های لایه های شبکه ایجاد شده است.
ساختار شبکه عصبی
برای بررسی نوسان های سطح کیفیت آب زیر زمینی دشت تهران، از روش شبکه عصبی مصنوعی نوع پرسپترون چندلایه استفاده شده است. این روش امکان بررسی تغییرات کیفی سفره ی آب زیرزمینی را در گام های مختلف زمانی با لحاظ اطلاعات گزینش شده در ورودی شبکه، فراهم می کند طراحی ساختارهای مختلف شبکه عصبی با توجه به مطالعات انجام شده در این زمینه، از نوع پرسپتترون چندلای (MLP ) با الگوریتم پس انتشارخطا (BP) استفاده شده است (17). در شبکه ی عصبی چندلایه، براساس شکل ارتباط بین عناصر ورودی در لایه ی اول( xi) و خروجی در لایه آخر (y )، به کمک وزن های نرونی (W)، بایاس ( b ) و تابع فعّالیّت ((f(x) در لایه (های) میانی است(شکل 2). طراحی شبکه براساس ترکیبی از اطلاعات پارامتر های تاثیر گذاردر کیفیت سفره در زمان های قبل، در قالب ساختارهای متفاوت از اطلاعات در لایه ی ورودی انجام شده است(16،17،18). در هر یک از ساختارها، اطلاعات ورودی پس از پردازش، از راه خروجی نرون های لایه ی اول، به نرون های لایه (های) بعدی و درنهایت، در صورت قابل قبول بودن به خروجی شبکه منتقل می شود. در غیر این صورت با انتشار خطای محاسبــاتی به لایه های پیش، محاسبه ها دوباره تکرار می شوند. این روند تا کسب نتیجه ی قابل قبول ادامه می یابد. برای افزایش سرعت پردازش اطلاعات و عدم توقف شبکه در حداقل های محلی، از اطلاعات نرمال شده به عنوان ورودی شبکه استفاده شده. از مزایای دیگر این نرم افزار ، وجودتوابع مختلف با الگوریتم های گوناگون در بانک نرم افزار است(16،17،18).
شکل2- مدل شبکه عصبی به کار رفته در تحقیق
در این مطالعه، داده ها به صورت تصادفی و براساس تجربه های پژوهشگران دیگر، به دو گروه داده آموزش70 درصد و داده آزمون30% تقسیم بندی شده اند. و کلیه مراحل به صورت آزمون و خطا بوده است.در اجرای شبکه عصبی مصنوعی، تابع محرک تانژانت استفاده شد. همچنین برای هر تابع محرک، قوانین مختلف آموزشی مومنتوم به کار رفته است برای دستیابی به ساختار بهینه، تلاش شد تا تمامی قوانین آموزشی و توابع محرک تعریف شده در محیط نرم افزار فوق، به شکل سعی و خطا مورد آزمون قرار گیرد. برای تعیین مقدار بهینه تعداد تکرار محاسباتی شبکه از روش سعی و خطا استفاده و خطای پیش بینی آن با مقادیر مختلفی از تکرارهای محاسباتی برای شبکه بهینه، محاسبه شد. گفتنی است که توابع محرک نرون های ورودی، میانی و خروجی یکسان در نظر گرفته شدند. در این راستا، مطالعات نیز نشان داد که یکسان در نظر گرفتن توابع محرک نسبت به متفاوت بودن توابع محــرک لایه هـــای مختلف، منجــر به نتــایج بهتـــری می شـــود (5،16،17،18). تعدادکل داده های مورد بررسی سری زمانی 1652 مورد می باشد. که 1156 مورد جهت آموزش شبکه و 496 پارامتر جهت بررسی و آزمون نهایی استفاده شده اند. جهت ساخت مدل، 3 مدل پیشنهادی به صورت حذفی با توجه به درصد اهمیت وابستگی که در جدول3 حاصل آزمون پیرسون مشخص شده است. پارامتر های ورودی هر مدل نیز در جدول(2) نمایش داده شده است.
جدول 2 - ساختار مدل های پیشنهادی در شبکه عصبی
مدل 1 |
مدل 2 |
مدل 3 |
|||
پارامتر ورودی |
پارامتر خروجی |
پارامتر ورودی |
پارامتر خروجی |
پارامتر ورودی |
پارامتر خروجی |
SO4 |
Ec |
SO4 |
Ec |
SO4 |
Ec |
Cl |
Cl |
Cl |
|||
HCo3 |
Na |
HCO3 |
|||
TH |
|
Na |
|||
SAR |
|
|
|||
K |
|
|
|||
Na |
|
|
|||
Mg |
|
|
|||
Ca |
|
|
جهت ورودی های مدل اول از کلیه پارامتر های موجود اندازه گیری شده جهت این تحقیق به عنوان ورودی استفاده شده است، در تلاش بعدی تعداد ورودی ها را به حداقل کاهش داده و پارامتر هایی که بیشترین مقدار همبستگی با پارامتر هدایت الکتریکی را دارا بوده انتخاب نموده و در تلاش آخر (مدل 3) با اضافه کردن پارامتر بی کربنات عملاً توان سیستم را در ارزیابی این پارامتر جهت بهبود پیش بینی چک می کنیم، در حقیقت هدف از اضافه کردن این پارامتر این است که آیا با اضافه کردن یک پارامتر اضافی تر سیستم قدرت پیش بینی بالاتری پیدا خواهد کرد و یا تاثیر محسوسی در روند تخمین نخواهد داشت.
با توجه به توضیحات بالا و ورود اطلاعات کامل هر مـدل در نرم افزار NeuroSolutions 5 و استفاده از سیستم تحت نرم افزار EXCEL و تعاریف مربــوط بــه شبکه عصبی در نرم افزار خروجــی هـای جــدول(4) ظاهر می شود. معیارهای پذیرش نتایج پیش بینی شده در هر یک از ســاختـارهــای شبکه ی عصبی، برمبنای کمتر بودن خطای کمی حاصل از مقادیر محاسباتی و مقادیر مشاهداتی شامل، بیشینه ضریب تعیین()،کمینه ریشه میانگین مربعـات خطا (RMSE) ، میانگین قدر مطلق خطا (MAE)، و خطای نسبی با روابط ذیل خواهد بود (6،8،10،19).
1)
2)
3)
در روابط فوق
در مقادیر فوق بیشینه ضریب تعیین( ) ،کمینه ریشه میانگین مربعات خطا(RMSE) ، نرمال کمینه ریشه میـانگین مربعــــات خطا (NRMSE). در مقــادیر فــوقEC actual میزان EC به دست آمده توسط آزمایش و نمونه گیری،EC forecast میزان EC محاسبه شده توسط شبکه عصبی ،EC average میزان میانگین EC به دست آمده از آزمایشات ،n تعداد پارامتر های مورد بررسی خواهد بود بود. ساختن پارامتر هــــای هر مدل با توجــه به همبستگی های جدول 3 خواهد بود و به ترتیب کمتر بودن اهمیت در مدل ها حذف می شوند.
باتوجه به نقشه شکل 3 خروجی نرم افزار ARC GIS می باشــد که مشخص کننده پراکنــدگی نقــاط نمونه گیری می باشد، معلوم می گردد می شود که بیشترین میزان رسانایی الکتریکی در مناطق جنوبی و شرق خواهد بود، واضح است با بررسی پارامتر هــای مختلف در ایـن مناطــق، مشخص می شود که بیشترین افزایش و کاهش همراه با پارامتر رسانایی الکتریکی، مربوط به پارامتر های کلراید، سولفات و سدیم می باشدکه تاکیدی بر درجه اهمیت مدل های ساخته شده می باشد(20،12،13).
شکل 3- نقشه پراکندگی پارامتر رسانایی الکتریکی در شهر تهران
جدول 3- میزان همبستگی پارامترهای ورودی شبکه عصبی با پارامتر رسانایی الکتریکی
K |
Ca |
Mg |
Na |
Cl |
HCo3 |
So4 |
SAR |
TH |
TDS |
Ec |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
Pearson Correlation |
Ec |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
993/0** |
Pearson Correlation |
TDS |
|
|
|
|
|
|
|
|
1 |
923/0** |
921/0** |
Pearson Correlation |
TH |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
345/0** |
617/0** |
633/0** |
Pearson Correlation |
AR |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
573/0** |
904/0** |
959/0** |
958/0** |
Pearson Correlation |
So4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
460/0** |
539/0** |
515/0** |
548/0** |
560/0** |
Pearson Correlation |
HCo3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
935/0** |
845/0** |
545/0** |
866/0** |
930/0** |
939/0** |
Pearson Correlation |
Cl |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
893/0** |
542/0** |
917/0** |
785/0** |
778/0** |
947/0** |
957/0** |
Pearson Correlation |
Na |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
790/0** |
814/0** |
583/0** |
889/0** |
401/0** |
943/0** |
902/0** |
900/0** |
Pearson Correlation |
Mg |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
804/0** |
693/0** |
830/0** |
407/0** |
833/0** |
264/0** |
956/0** |
855/0** |
853/0** |
Pearson Correlation |
Ca |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
853/0** |
853/0** |
853/0** |
853/0** |
853/0** |
853/0** |
853/0** |
853/0** |
853/0** |
853/0** |
Pearson Correlation |
K |
|
|
|
|
|
|
|
|
|
|
|
|
یافته های تحقیق
با توجه به ساخت سه مدل مختلف شبکه عصبی با ساختار های ذکر شده در جدول 4 و در نظر گرفتن، ثابت بودن توابع آموزش و محرک نتایج ذیل با سعی و خطا حاصل گردیده است.
جدول 4- ساختار شبکه های عصبی و مقادیر خطاهای به دست آمده
مدل |
تعداد پارامتر ورودی |
تعداد پارامتر خروجی |
تعداد لایه پنهان |
درصد پارامتر های آموزشی |
درصد پارامترهای پاسخی |
تعداد پردازنده لایه پنهان |
قانون آموزش |
تابع لایه پنهان |
قانون آموزش لایه خروجی |
تابع محرک لایه خروجی |
1 |
9 |
1 |
1 |
70 |
30 |
11 |
مومنتوم |
تانژانت |
مومنتوم |
تانژانت |
2 |
3 |
1 |
1 |
70 |
30 |
38 |
مومنتوم |
تانژانت |
مومنتوم |
تانژانت |
3 |
4 |
1 |
1 |
70 |
30 |
28 |
مومنتوم |
تانژانت |
مومنتوم |
تانژانت |
خطای حاصل از اجرای هر یک از سه مدل شبکه عصبی مصنوعی در پیش بینی پارامتر رسانایی الکتریکی در جدول 5 آورده شد.
جدول 5 -خطاهای به دست آمده از مدل های شبکه عصبی
مدل |
فرایند |
RMSE |
NRMSE |
تعداد ورودی |
|
1 |
آموزش |
986/0 |
21/161 |
110/0 |
1156 |
تست |
985/0 |
170 |
113/0 |
496 |
|
2 |
آموزش |
964/0 |
76/261 |
179/0 |
1156 |
تست |
963/0 |
62/262 |
176/0 |
496 |
|
3 |
آموزش |
986/0 |
12/160 |
110/0 |
1156 |
تست |
986/0 |
95/159 |
110/0 |
496 |
با بررسی مــدل هــا مشخص می شود که جهت پیش بینی پارامتر رسانایی الکتریکی در مدل اول و سوم با توجه به تعـــداد پارامتر های ورودی بیشتر دارای خطای بیشتری نسبت به مدل دوم می باشند، با بررسی دقیق تر این موضوع مشاهده می شود که با کم کردن پارامتر های ورودی و ثابت نگه داشتن توابع محرک و آموزشی میزان خطا های مدل بیشتر خواهد شد در صورتی که در صورت کم کردن تعداد پارامتر های ورودی و بیشتر کردن عوامل پردازشی لایه پنهان تا حد مشخص از مقدار خطا کاسته می شود. در مدل های اول تا سوم با ثابت نگه داشتن توابع آموزش و محرک میزان خطا کاسته می شود. همچنین معلوم می شود که با کاهش تعداد نرون های ورودی و افزایش پردازنده های لایه پنهان در کمترین حالت قرار خواهد گرفت. در مقایسه مدل های اول و سوم با مدل دوم ،این مدل با توجه به تعداد کمتر نرون ورودی دارای خطای قابل قبول خواهد بود(جدول 5). در نمودار شکل (4) پراکندگی داده های پیش بینی شده (محـور عمودی) و داده های مشاهده ای(محور افقی) در مرحله آزمون مدل اول تا سوم خواهد بود .در شکل(4) ارتباط بین داده های ورودی و خروجی شبکه عصبی به شکل یک معادله درجه یک نشان داده شده است و میزان انحراف از معیار از نیمساز ربع اول مشخص شده است.لازم به ذکر است هر چقدر داده ها به نمودار یک به یک نزدیک تر باشند ،نشان دهنده توانایی بیشتر مدل در سنجش میزان رسانایی الکتریکی خواهد بود. در شکل محور افقی نمایان گر میزان هدایت الکتریکی پیش بینی شده توسط مدل مربوطه خواهد بود و محور عمودی نشان دهنده مقدار واقعی پارامتر هدایت الکتریکی که توسط آزمایش به دست آمده است می باشد(20،21).
شکل 4 - نمودارهای پراکندگی میزان رسانایی الکتریکی پیش بینی شده و مشاهده شده
همچنین در بررسی شکل5 به طور کامل میزان همپوشانی مقادیر محاسبه شده توسط مدل های مختلف شبکه عصبی و میزان مشاهده شده توسط نتایج آزمایشات نمایش داده شده است. واضح است که در هر سه مدل بالای 90% مقــادیر پیش بینی و واقعی با هم مساوی خواهد بود و دو نمودار بر هم منطبق شده اند.در شکل محور افقی نمایانگر شماره نمونه بوده و محور عمودی مقدار آن را جهت پیش بینی توسط مدل و نتایج آزمایشات با توجه به همپوشانی نمودار های روی هم مشخص می نماید. بدیهی است هر چقدر این همپوشانی دارای مقدار بیشتری باشد میزان پیش بینی توسط مدل دقیق تر خواهد بود.
شکل 5- نمودارهای همپوشانی میزان رسانایی الکتریکی پیش بینی شده و مشاهده شده
با بررسی شکل6 معلوم می شود که میزان حساسیت پارامتر های مختلف در پیش بینی شبکه عصبی مدل سوم به ترتیب برای پارامترهای سولفات و کلراید و سدیم می باشد. علت این تغییر اهمیت نسبت به جدول همبستگی ممکن است تغییرات کمتر نسبت به دو پارامتر دیگر باشد.در حقیقت به دلیل این که سدیم تا حد زیادی جزو بافت خاک منطقه مورد مطالعه می باشد (20)، تقریبا دارای میزان ثابت تری نسبت سولفات و کلراید می باشد که در پیش بینی تاثیرات کمتری را دارا خواهد بود(20).
شکل 6- میزان حساسیت پارامتر های مختلف درتخمین شبکه عصبی مدل 2
نتیجه گیری
با توجه به نتایج به دست آمده از ساخت مدل های مختلف شبکه عصبی مشاهده می شود که مدل دوم با کمترین تعداد داده ها و در نتیجه کم ترین مقدار آزمایشات مربوط به بدست آوردن پارامتر های مختلف دارای کمترین هزینه و میزان خطای قابل قبول به جهت پیش بینی پارامتر رسانایی الکتریکی آب زیر زمینی شهر تهران خواهد بود. این مدل با تابع آموزش مومنتوم و تابع محرک تانژانت دارای نتایج مطلوب به شرح ذیل می باشد.از دلایل بهتر شدن نتایج می توان بازه باز تر تابع تانژانت (1و1-) بوده را نام برد(16،17،18). نتایج به دست آمده نشان دهنده حداکثر خطای نتایج به دست آمده نشان دهنده میزان بیشینه ضریب تعیین برابر 0.96 می باشد که با توجه به مقدار پارامتر ورودی به نسبت مدل های دیگر دارای مقدار خطای کمتری خواهد بود .همچنین با توجه به شبکه عصبی ساخته شده در یک لایه میزان کمینه ی نرمال ریشه میانگین مربعات خطا برابر 0.176 خواهد بود که با توجه به مقدار ورودی ها و شبکه های عصبی ساخته شده در مدل های دیگر دارای تقریب بسیار مناسب و قابل قبول می باشد. در حقیقت این مدل 98% از بیشترین مقدار پیش بینی که در مدل اول می باشد را دارا می باشد و با تقریب 98% پیش بینی را انجام خواهد داد. با توجه به مدل به دست آمده می توان با کنترل پارامتر سولفات و کلراید و سدیم در پیش بینی پارامتر رسانایی الکتریکی در آینده تقریب درستی را ارائه نمود و در جهت بهبود کیفیت و مدیریت منابع آب زیر زمینی توان کافی پیش بینی را دارا بود. همچنین در بررسی میزان حساسیت های پارامتر های ورودی و تاثیرات آن بر نتایج شبکه عصبی مدل پنجم مشخص می شود که به ترتیب پارامتر های سولفات و کلراید و سدیم دارای تاثیرات خواهند بود .جدول همبستگی پیرسون نیز تاکیدی بر صلاحیت شکل به دست آمده از شبکه عصبی خواهد بود. تفاوت در میزان حساسیت پارامتر های سدیم و کلـر نیز می تواند جزو نکات مهم بررسی کیفــی آب زبر زمینی تهران باشد. یکی از دلایل مهم در این مطلب زمین ساخت بودن پارامتر سدیم که در خصوصیات زمین شناسی منطقه به دلیل وجود دریاچه های نمک در جنوب شهر دارای اهمیت می باشد و دارای تغییرات کمتری در ورود سفـره آب زیر زمینی خواهد بود(2،3)، و انسان ساخت بودن پارامتر کلرکه حاصل فعالیت های شهری و صنعتی و ورود مواد شوینده و ضد عفونی کننده می باشد که میزان ورود بیشتری به سفـره را در فصــول و مکان های مختلف دارا می باشد(2).
منابع
1- استاد دانشکده محیط زیست دانشگاه تهران.
2- دانشیار دانشکده محیط زیست دانشگاه تهران.
3- دانشیار گروه مهندسی آب دانشکدة کشاورزی، دانشگاه بوعلی سینا همدان
4- دانشجوی کارشناسی ارشد پردیس ارس دانشگاه تهران*(مسئول مکاتبات).