نوع مقاله : مقاله پژوهشی
نویسندگان
1 کارشناس ارشد اقلیمشناسی، دانشگاه زنجان، ایران. *(مسوول مکاتبات)
2 دانشیار، گروه جغرافیا، دانشگاه پیام نور، تهران، ایران.
چکیده
کلیدواژهها
فصلنامه انسان و محیط زیست، شماره 47، زمستان 97
واسنجی تابش دریافتی سطح زمین در پهنه شمالشرق ایران
مهناز رستمیان [1] *
امیرحسین حلبیان[2]
تاریخ دریافت: 29/07/1396 |
تاریخ پذیرش: 16/12/1396 |
چکیده
تابش خورشیدی رسیده به سطح زمین (Rs)، در بسیاری از زمینههای مهندسی کشاورزی، آبیاری و هیدرولوژی کاربرد گستردهای دارد و به این دلیل روشهای مختلفی برای تخمین آن ارایه شده است که معادله آنگستروم – پرسکات یکی از مهمترین معادلات موجود در این زمینه است. این معادله دارای ضرایب تجربی است که برای هر موقعیت متفاوت میباشند. به دلیل اهمیت Rs، روشهای مختلفی برای تخمین آن توسعه داده شد. پژوهش حاضر با استفاده از دادههای میانگین روزانه ساعت آفتابی ایستگاههای سینوپتیک وابسته به اداره کل هواشناسی کشور در گستره خراسانشمالی، خراسانرضوی و خراستان جنوبی، طی بازهی 2016-1988 به برآورد تابش رسیده به سطح زمین پرداخته است. جهت برآورد تابش خورشیدی در سطح زمین از روش انگستروم_ پرسکات بهره گرفته شد. در این پژوهش با استفاده از روش فرر ضرایب a و b محاسبه و سپس در معادله آنگسترم جایگذاری گردید، در نهایت با استفاده از مدل درونیابی کریجینگ نقشههای همارزش تابش گستره مورد مطالعه ترسیم و تحلیل گردید. نتایج این پژوهش نشان داد که مناطق جنوب، جنوب شرق منطقه مورد مطالعه، بهخصوص ایستگاه نهبندان پتانسیل دریافت تابش خورشیدی، بیشتری دارند. کمینه مقادیر تابش دریافتی ماهانهی در ماههای آذر و دی در شمال و شمالغرب منطقه مورد مطالعه قرار دارد، که به ترتیب مقادیر 36/7، 98/7 مگاژول بر متر مربع میباشد. بیشینه مقدار تابش دریافتی، در ماه تیر با مقدار 24/30 مگاژول بر مترمربع اتفاق میافتد و مناطق جنوب و جنوبشرق تا مرکز پهنه بهخصوص ایستگاه نهبندان و به ترتیب فصول تابستان و بهار از نظر دریافت تابش خورشیدی، پتانسیل بیشتری دارند.
کلمات کلیدی: تابش خورشیدی، ساعت آفتابی، مدل فرر، آنگستروم- پرسکات، شمالشرق ایران.
|
Calibration of earth surface received radiation in north east of Iran
Mahnaz Rostamian [3]*(Corresponding Author)
Amir Hossein Halabian [4]
Abstract
The solar radiation reached the earth surface (RS) has wide applications in many areas of agriculture engineering, irrigation and hydrology and that’s why various methods have been represented for its estimation, which Angstrom – Prescott equation is one of the most important equation in this field. This equation has experimental coefficients that are different for any situation. Because of the importance of Rs, various methods have developed for its estimation. The radiation reached to earth surface is evaluated by using average daily sunny hours data of synoptic affiliated to the general state metrological office in the north of Khorasan, Khorasan Razavi and southern Khorasan 1988- 2016. The Angstrom – prescient method has been used for estimating solar radiation of earth surface in the studied area. In this research, the a & b coefficients were calculated by using Ferrer's method fastly and then have been placed in Angstrom equation. Finally, the radiation equivalent map of the studied area has been drawn and analyzed by using kriging interpolation model. The research outcome indicated that southern, south-east regions and at the next level, the center of the studied area especially Nehbndan station have more potential for receiving solar radiation. The minimum monthly received radiation of area during January and December month in the north and northwest of this region, which are 7.98 and 7.36 megajoules per square meter respectively. The maximum received radiation of area occurs in July about 30.24 megajoules. South and south-east regions up to the center of the of the area Nehbandan station have more potential for receiving in terms of radiation during summer & spring seasons respectively.
Key Words: Solar Radiation, Sundial, Ferrer Model, Angstrom- Prescott, North East of Iran.
مقدمه
انرژی تابشی خورشید سرمنشا انرژیهای روی کره زمین است که یکی از منابع مهم انرژی پاک به ویژه در کشور ما است. انرژی خورشید از معروفترین انرژیهای پاک و انرژیهای غیر فسیلی است و استفاده از آن کمترین آلودگی را در محیط زیست به وجود میآورد. دریافت اطلاعات انرژی خورشید به دو روش امکانپذیر است. روش اول اندازهگیری مستقیم انرژی خورشید توسط دستگاههای اندازهگیری مانند پیرانومتر و پیرهیلومتر میباشند. در روش دوم که با استفاده از روشهای تجربی به برآورد تابش میپردازد. مهمترین پارامتر که در بیشتر معادلات تجربی برآورد تابش خورشیدی از آن استفاده میشود، ساعات آفتابی است. تحقیقات نشان داده است که ساعات آفتابی مهمترین عامل تعیین کنندهی میزان تابش خورشیدی است. هدف از پژوهش حاضر واسنجی و برآورد مقدار انرژی دریافتی در سطح زمین در سه استان خراسانشمالی، رضوی و جنوبی، با استفاده از روش تجربی آنگستروم[5] و نیز با کمک داده های ساعت آفتابی و عرض جغرافیایی می باشد، از آنجا که هزینه دستگاههای اندازه گیری بالا بوده؛ با استفاده از روشهای فیزیکی و محاسبات سعی می شود تا مقدار تابش با کمترین مقدار خطا برآورد گردد.
پیشینه پژوهش
در ایران و جهان تاکنون پژوهشهای زیادی توسط محققان مختلف در زمینه برآورد انرژی خورشیدی با روشهای مختلفی صورت گرفته است. از جمله این مطالعات می توان به موارد زیر اشاره نمود: چو و لی[6] با استفاده از روش انگستروم به برآورد تابش خورشیدی در مالزی پرداختند. نتایج مطالعه نشان داد که در طول یک ماه همبستگی خطی وابستگی رخداد فرینهای بیشینه و کمینه روزانه طول ساعت آفتابی بسیار ضعیف است. برآوردهای بدست آمده برای شهرهای نزدیک تایلند و سنگاپور موافقت بیشتری را با دادههای این کشورها نشان دادند (1). لوپز و همکاران[7] به برآورد تابش خورشیدی از دادههای مودیس[8] در جنوب اسپانیا پرداختند. آنها یک پدیده دیگر به مدل پارامتری برای تخمین تابش خورشیدی تحت شرایط ابرناکی اضافه نمودند. محتوای بخار آب قابل بارش جو تنها ورودی تجربی مدل میباشد (2). ماکوچی و همکاران[9] مدلی جدید برای برآورد و ارزیابی تابش خورشیدی در مراکش شمالی ارایه دادند نتایجشان نشان میدهد که مدل پیشنهادی برای برآورد تابش خورشیدی در تمام فصلهای سال برای موقعیت مورد مطالعه و برای روز مورد نظر با استفاده از ورودیهای ارتفاع (متر)، طول جغرافیایی (درجه) و عرض جغرافیایی (درجه) مورد استفاده قرار می گیرد (3). ژانک و همکاران[10] به مرور انتقادی بر مدلهای مورد استفاده برآورد تابش خورشیدی طی یک دهه قبل پرداختند. در این مطالعه روشهای مورد استفاده در مقالات مرتبط با تابش خورشید مورد بررسی قرار گرفتهاند. نتایج نشان میدهد که مدلهای کسری زمانی تابش آفتاب و شبکههای عصبی مصنوعی عملکرد مشابهی دارند که برای برآورد ماهانه متوسط رزوانه تابش جهانی استفاده میشود. در حالی که بیشتر کارها برای مطالعه روش برآورد به فواصل زمانی کوتاهتر و مکانیزمهای اتمسفر برای تابش خورشید نیاز دارند (4). ربیع و همکاران[11] با استفاده از دادههای چند سال اخیر و مدل شبکه های عصبی مصنوعی برای برآورد متوسط روزانه تابش خورشیدی در کویت پرداختند. مدل شبکههای عصبی مصنوعی قادر بوده با استفاده از روش غیرخطی به پیشبینی بپردازد. توسعه پیشبینی وابسته به این است که دیدگاههای مقامات کشوری، سرمایهگذاران و نیروی مهندسی سیستم، و تعیین مکانهای مناسب برای تاسیسات خورشید و شایستگی فنی-اقتصادی با یکپارچه سازی انرژی خورشیدی در مقیاس طولانی مدت باشد (5). یالدریم و همکاران[12] به برآورد روزانه تابش جهانی خورشید با استفاده از رابطه کاربری گرافیکی در منطقه مدیترانه شرقی ترکیه در یک دوره 11 ساله پرداختند. آنها روشهای مختلف مدل های مختلف تابش خورشیدی را تحلیل نمودند. همچنین یک مدل بهینه شده بر اساس شبکه عصبی مصنوعی و مدل آنگستروم - پرسکات برای برآورد تابش خورشیدی جهانی روزانه ارایه دادند. و یک مدل بهینه شده بر اساس شبکه عصبی مصنوعی و مدل آنگستروم - پرسکات [13]برای برآورد تابش خورشیدی جهانی روزانه ارایه دادند (6) . آچور[14] مدلی ترکیبی برای برآورد تابش خورشیدی ماهانه در جنوب الجزایر برای دوره 5 ساله ارایه دادند. نتایج به دست آمده آنها نشان داد که مدل جدید هیبرید میتواند قادر به پیش بینی تابش خورشیدی جهانی با دقت عالی در این مکان باشد (7). رودریگو و همکاران[15] با استفاده از تصاویر ماهوارهای و الگوریتم ژنتیکی به برآورد تابش جهانی خورشید پرداختند. آنها از تصاویر با یاختههای 3 کیلومتر و در بازهی ژانویه تا دسامبر 2008 بهره بردند. برآیند مطالعه ضریب همبستگی 99 درصد را نشان داد. برآوردهای تولید شده نسبتا دقیق بوده و خطای فضایی نسبتا یکنواختی را نشان دادند. مدل فوق حتی برای روزهای ابری نیز برآوردهای قابل اعتمادی را نشان داده است (8). خلیلی و رضایی صدر به برآورد تابش کلی خورشید در گستره ایران بر مبنای دادههای اقلیمی در دوره (1991-1966) باروش انگستروم پرداختند و نتایج آنها نشان داد که به جز ایستگاههای جاسک و زاهدان که آنالیزها منجر به رابطه معنیداری برای آنها نشده است، در سایر ایستگاهها روابط در سطح 01/0 معنی دار میباشد. ضریب A معادله آنگستروم در شبکه تابشسنجی ایران بین 21/0 در یزد تا 39/0 در کرمانشاه متفاوت میباشد (9). صفائی و همکاران در تحقیقی به برآورد پتانسیل تابش خورشیدی در ایران و تهیه اطلس تابشی آن در یک دوره 5 تا 11 ساله با استفاده از سه مدل ریاضی: آنگستروم، روش حداکثر احتمال و مدل هیبرید اقدام کردند و مدل هیبرید بهترین برازش را با این دادهها داشت و حاصل کار آنها 12 نقشه مربوط به میانگین روزانه تابش کلی خورشیدی در هر ماه است (10).علیزاده و خلیلی در پژوهشی به تعیین ضرایب معادله آنگستروم و توسعه یک معادله رگرسیونی برآورد تابش خورشیدی در مشهد در سالهای 1380 الی 1378 ضرایب رابطه آنگسترم – پرسکات واسنجی پرداختند و همچنین یک معادله رگرسیونی منطقهای ارائه دادند نتایج آنالیز آماری روابطی که بررسی کردند نشان داد که دقت معادله منطقه ای چند متغیره با رابطه آنگستروم – پرسکات اختلاف معنی داری نداشته است (11). یزدان پناه و همکاران در تحقیقی برآورد تابش کلی خورشید در سطح افقی زمین در اصفهان با استفاده از آمار هواشناسی ایستگاه ازن سنجی اصفهان در دوره آماری 1384 تا 1387 نسبت به یافتن مدل مناسب برآورد شدت تابش خورشیدی رسیده به سطح از روی آمار ساعات آفتابی با استفاده از مدل آنگستروم اقدام کردند ونتایج این تحقیق نشان داد مدل تهیه شده بهویژه برای ماههای می و آگوست دقت بسیار بالایی دارد. و مقدار ضریب رگرسیون مدل ساخته شده ما معادل 91/0 است، ضرایب رگرسیون برای مدلهای، فرر و ام سی کلوچ 6/0 و 7/0 است (12). معینی و همکاران در تحقیقی به برآورد تابش خورشیدی در ایران با استفاده از یک مدل بهینه پرداختند و از مدل آنگستروم برای برآورد میزان تابش خورشید روی سطح افقی بهره گرفتند و تابش کل را برای نقاط مختلف ایران وهمچنین محاسبه طول روز و تابش فراجو را برآورد کردند (13). سبزیپرور و بیاتورکشی به ارزیابی دقت روشهای شبکه عصبی مصنوعی و عصبی – فازی در شبیهسازی تابش کل خورشیدی طی سالهای 1992 تا 2006 در ایستگاههای کرمان، ارومیه، شیراز و اصفهان پرداختند و و نتیجه گرفتند روابط تجربی مورد استفادشان در تخمین تابش خورشیدی موید ارتباط مناسب بین پارامترهای هواشناسی و مقدار تابش خورشیدی میباشد و بهترین مدل شبکه عصبی مصنوعی با تابع فعالیت تانژانت هایپربولیک در ایستگاه کرمان میباشد (14). موسویبایگی و اشرف در تحقیقی به شناسایی مناطق کمترین میزان ابرناکی به منظور پهنهبندی نواحی پرتابش کشور در دوره 1387-1368 با روش درون یابی اسپیلاین spline)) اقدام نمودند و نتایج آنها نشان داد که بیشترین روزهای آفتابی در فروردین ماه متعلق به شهرستان بافت در اردیبهشت ماه متعلق به شهرستان خاش بوده است و همچنین اغلب مناطق مناطق ایران از استعداد و قابلیت بالایی برای استفاده از انرزی تجدیدپذیر خورشیدی برخوردارند (15). پیری و همکاران در پژوهشی به مدلسازی تابش خورشید رسیده به زمین با استفاده از ANFIS و مدلهای تجربی ایستگاههای زاهدان و بجنورد در دوره آماری به ترتیب 5 تا 7 ساله پرداختند و دو روش تجربی را جهت برآورد تابش روزانه خورشید واسنجی و ارزیابی کردند و نتیجه گرفتند که مدل نرو فازی برآورد بهتری نسبت به دو مدل تجربی در برآورد تابش دارد (16). زمانی محی آبادی در تحقیقی به پیش بینی لحظهای تابش کل خورشید در شهر رفسنجان توسط شبکه عصبی برای یک سال پرداخت و نتیجه گرفت که شاخصهای آمار میانگین خطای بایاس [16]و ریشه میانگین مربعات خطا [17]و t برای شبکه عصبی بترتیب 0037/0، 5366/36 و 0232/0 میباشند که نمایانگر عملکرد خوب مدل طراحی شده در شبکه عصبی است (17). عرفانیان و بابایی حصار در پژوهشی به ارزیابی مدل هیبرید در تخمین تابش خورشیدی روزانه در تعدادی از ایستگاههای تابشسنجی ایران پرداختند. و مقدار کم معیارهای خطا برای مدل هیبرید بیانگر کارایی نسبتا بالای این مدل در تخمین تابش خورشیدی ایستگاههاس انتخابی در مقیاس روزانه میباشد (18). هوشنگی و همکاران در تحقیقی به بررسی منطقهای پتانسیل تابش خورشیدی با ارزیابی و بهینه سازی روشهای درون یابی در سطح کشور ایران و نتایج حاصل با اعتبار سنجی متقابل مقایسه شدند. و بررسی نتایجشان نشان داد که بهینه سازی عوامل موثر در روشهای درون یابی موجب افزایش همگرایی 3/8 درصدی مقادیر نتایج میگردد (19).کرمپور و همکاران در پژوهشی توازن تابش در سطح شهر یزد در ماه ژانویه 2006 تا 2008 و میانگین ماهانه ساعتی سپیدایی، شار تابش خالص در شرایط هوای صاف و ابری واکاوی کردند و به این نتیجه رسیدند که میزان سپیدایی غرب شهر به دلیل ساخت شهری بیشتر از شرق شهر است. و مقدار آن در روزهای که سطح شهر مرطوب بوده تا 2 درصد کاهش پیدا میکند(20). مجرد و همکاران با به کارگیری مدل بهینه شدهی برد و هول استورم در دورهی زمانی 2009-1990 اقدام به برآورد تابش خورشیدی دریافتی سطح زمین در استان کرمانشاه نمودند. نتایج تحقیق آنها نشان داد بین مقادیر تابش اندازه گیری شده با پیرانومتر و مقادیر برآورد شده توسط مدل اختلاف کمی وجود دارد که نشان دهندهی توانایی مدل در برآورد تابش است. و در فصول پاییز و زمستان به دلیل نقش تابش پراکنده، اختلاف تابش در مناطق مختلف استان کاهش مییابد (21). فلاح قالهری و شاکری با به کارگیری مدلهای آنگستروم – پرسکات به واسنجی ضرایب رابطه آنگستروم – پرسکات در ایستگاههای منتخب استان خراسان رضوی نمودند. نتایج آنها نشان داد که بین مقادیر تابش اندازهگیری شده با پیرانومتر و مقادیر تابش برآورد شده توسط رابطه بلکه اختلاف کمی وجود دارد. و همچنین نتایج حاصل از پهنهبندی تابش نشان داد مناطق جنوبی استان از نظر دریافت تابش خورشیدی، پتانسیل بیشتری دارند (22).
مواد و روش
منطقه مورد مطالعه
منطقه مورد مطالعه در این تحقیق، پهنهی شمالشرق متشکل از استانهای خراسانجنوبی، خراسان رضوی وخراسان شمالی واقع در طول جغرافیایی از 55 درجه و 55 دقیقه و53 ثانیه غربی تا 61 درجه و 17 دقیقه و 2 ثانیه شرقی و عرض جغرافیایی از 30 درجه و 29 دقیقه و 33 ثانیه جنوبی تا 38 درجه و 18 دقیقه و 3 ثانیه جنوبی است مساحت پهنه 5/ 239730 کیلومتر مربع میباشد(شکل 1). این استان ها تحت تاثیر توده هوای سرد قطبی، سیبری ، توده هوای مدیترانهایی و توده هوای گرم و خشک (پرفشار جنب حاره ) قرار دارند که شرایط اقلیمی متفاوتی را برای این استانها رقم می زنند.
شکل1-موقعیت جغرافیایی منطقه مورد مطالعه به همراه ایستگاهای منتخب
همانطور که در جدول 1 مشاهده میشود کمترین میانگین مجموع ساعات آفتابی به ترتیب ایستگاه بجنورد در استان خراسانشمالی با 6/7 و در مرتبه بعدی ایستگاههای قوچان، گلمکان، سرخس، مشهد در استان خراسانرضوی با مقادیر 7/7، 0/8،0/8، 1/8 دارا میباشند و بیشترین میانگین مجموع ساعات آفتابی در جنوبشرق و مرکز پهنه به ترتیب در ایستگاههای نهبندان و قائن و بیرجند و فرودس با مقادیر 2/9، 1/9 به خود اختصاص دادهاند. انحراف معیار میزان پراکندگی دادهها را نسبت به میانگین میسنجد.کمترین نوسانات مجموع ساعات آفتابی در طول پهنه در ایستگاههای نهبندان و بیرجند با مقادیر (7/1، 9/1) و بیشترین نوسانات در ایستگاههای سرخس و مشهد و گلمکان به ترتیب با مقادیر (8/2، 7/2، 5/2) مشاهده شد.
جدول 1- آمار توصیفی ساعات آفتابی ایستگاههای سینوپتیک منطقه مورد مطالعه طی سالهای (2016-1988)
بیشینه |
کمینه |
واریانس یا پراش |
انحراف معیار |
میانگین |
ایستگاه |
استان |
8/11 |
3/5 |
8/4 |
2/2 |
3/8 |
بشرویه |
خراسان جنوبی
|
9/11 |
6/6 |
6/3 |
9/1 |
1/9 |
بیرچند |
|
1/12 |
0/6 |
2/4 |
1/2 |
1/9 |
فردوس |
|
3/12 |
2/6 |
6/4 |
1/2 |
1/9 |
قائن |
|
8/11 |
8/6 |
8/2 |
7/1 |
2/9 |
نهبندان |
|
3/12 |
8/5 |
5/5 |
3/2 |
9/8 |
تربت حیدریه |
خراسان رضوی
|
5/11 |
6/5 |
7/4 |
2/2 |
5/8 |
سبزوار |
|
0/12 |
3/4 |
9/7 |
8/2 |
0/8 |
سرخس |
|
3/11 |
6/4 |
0/6 |
4/2 |
7/7 |
قوچان |
|
9/11 |
8/5 |
7/4 |
2/2 |
8/8 |
کاشمر |
|
6/11 |
0/5 |
2/6 |
5/2 |
0/8 |
گلمکان |
|
0/12 |
0/6 |
3/4 |
1/2 |
8/8 |
گناباد |
|
1/12 |
6/4 |
3/7 |
7/2 |
1/8 |
مشهد |
|
9/10 |
6/4 |
5/5 |
3/2 |
6/7 |
بجنورد |
خراسان شمالی |
روش تحقیق
در این تحقیق برای برآورد میزان تابش دریافتی سطح زمین در استانهای خراسانرضوی، خراسانجنوبی و خراسان شمالی از میانگین دادههای روزانه ساعات آفتابی 13 ایستگاه سینوپتیک منطقه که دارای آمار طولانی، طی یک بازهی زمانی 28 ساله (2016-1988) بودند از سازمان هواشناسی کشور دریافت شد و سپس استفاده گردید، البته اندازهگیری مستقیم انرژی خورشید توسط دستگاههای اندازهگیری مانند پیرانومتر [18] پیرهیلومتر[19] میباشد. در منطقه مورد بررسی، اکثر ایستگاهها مجهز به دستگاههای اندازهگیری تابش خورشید نبودند، یا اندازهگیری تابش در آن نقاط به صورت پراکنده میباشد لذا به برآورد مقدار تابش دریافتی سطح زمین در شمالشرق ایران پرداخته شد، در این تحقیق ابتدا با استفاده از مدل فرر ضرایب a و b محاسبه گردید سپس با جایگذاری این ضرایب در معادله آنگستروم – پرسکات انرژی خورشید برآورد گردید. (روابط 6-1). مقدار تابش دریافتی سطح زمین در خراسانهای رضوی،شمالی، جنوبی برای12 ماه براساس روش کریجینگ[20] نقشه آنها در Arc GIS ترسیم شد شکل (2).
مدل انگستروم :
پریسکات از طریق رابطه زیر برآورد می گردد (21):
(1) |
|
در رابطه فوق بیانگر میانگین تابش کل روزانه درهر ماه میباشد.
(2)
|
بیانگر تابش اندازهگیری شده در خارج از جو میباشد مقدار از رابطه (2) بدست میآید: که در آن شماره ژوولیوسی ، ثابت خورشید که معادل 1367، عرض جغرافیایی محل و زاویه میل که از رابطه (3) زیر به دست میآید: ( 12).
(3) |
که در آن ( زاویه ساعتی خورشید در هنگام غروب) از رابطه زیر به دست میآید:
(4) |
بیانگر میانگین ماهیانه ساعات آفتابی روزانه و میانگین ماهیانه ماکزیمم ساعات آفتابی (طول روز) میباشد و مقدار آن از رابطه (5) محاسبه میشود.
(5) |
از مدل فرر برای بدست آوردن مقادیر (a,b) در برآورد تابش کلی ( درسطح افقی) استفاده شده است( رابطه6). مدل فرر به شرح زیر است: (12).
(6) |
a= -0.27+1.75 b= 1.32-2.90 |
یافته ها
همان طور که در جدول 2 مشاهده میشود، بیشترین میانگین تابش خورشیدی ماهانه در جنوب شرق پهنه مورد نظر در ایستگاه نهبندان بامقدار 4/20 و در مرتبه های بعدی ایستگاه های بیرجند و قائن ، فردوس ، گناباد و تربت حیدریه به ترتیب با مقادیر 9/19، 9/19، 7/19، 3/19، 3/19 مگاژول بر متر مربع را به خود اختصاص دادهاند. اکثر ایستگاهها نوسانات تابش خورشیدی درون سالانه نسبتا زیادی را نشان دادند. ایستگاه سرخس با انحراف معیار 5/7 بیشینه نوسانات را نشان داده است. و کمینه نوسانات در ایستگاه نهبندان و با انحراف معیار 6 مشاهده شده است. در این بین ایستگاه های سرخس و بجنورد حداقل متوسط تابش خورشیدی دریافتی را با حدود 5/7 و 6/7 مگاژول بر متر مربع را نشان دادند. بیشینه متوسط روزانه تابش را ایستگاه قائن و تربت حیدریه به ترتیب با حدود 3/30 و 30 مگاژول بر متر مربع را تجربه کرده است.
جدول2- آمارتوصیفی تابش خورشیدی ماهانه (2016-1988) برحسب مگاژول بر متر مربع در روز
ایستگاه |
میانگین |
انحراف معیار |
واریانس |
کمینه |
بیشینه |
بیرجند |
9/19 |
3/6 |
3/40 |
0/11 |
3/29 |
بشرویه |
9/18 |
7/6 |
0/45 |
2/9 |
9/28 |
نهبندان |
4/20 |
0/6 |
0/36 |
9/11 |
1/29 |
قائن |
9/19 |
9/6 |
9/47 |
3/10 |
3/30 |
فردوس |
7/19 |
7/6 |
6/44 |
3/10 |
6/29 |
قوچان |
3/17 |
9/6 |
2/47 |
8/7 |
2/27 |
گلمکان |
9/17 |
0/7 |
0/49 |
2/8 |
1/28 |
گناباد |
3/19 |
6/6 |
7/43 |
0/10 |
3/29 |
کاشمر |
1/19 |
8/6 |
3/46 |
5/9 |
9/28 |
مشهد |
0/18 |
4/7 |
8/54 |
1/8 |
2/29 |
سبزوار |
5/18 |
7/6 |
5/45 |
9/8 |
8/27 |
سرخس |
9/17 |
5/7 |
9/56 |
5/7 |
0/29 |
تربت حیدریه |
3/19 |
2/7 |
6/51 |
4/9 |
0/30 |
بجنورد |
1/17 |
7/6 |
7/44 |
6/7 |
3/26 |
باتوجه به نقشههای شکل 2 مشاهده میشود کمترین مقادیر تابش دریافتی ماهانهی پهنه در ماههای دی و آذر در شمال و شمال غرب پهنه رخ میدهد که به ترتیب مقادیر 98/7، 36/7 مگاژول بر متر مربع میباشد. بیشترین مقدار تابش دریافتی پهنه در ماه تیر با مقدار 24/30 مگاژول بر مترمربع اتفاق میافتد ودر مرتبه های بعدی به ترتیب ماههای خرداد، مرداد، اردیبهشت و فروردین 81/28، 78/27، 72/25 و 24/22 مگازول بر متر مربع در جنوبشرق و مرکز پهنه مورد نظررا به خود اختصاص دادهاند. افزایش مقادیر تابش در ماههای گرم سال، در مرحله نخست به خاطر افزایش زاویه خورشید از فصل سرد به گرم میشود. از طرفی دیگر در ماههای سرد سال ابرناکی آسمان بیشتر است و باعث کاهش دریافتی خواهد بود در فصل زمستان به دلیل نقش تابش پراکنده، اختلاف تابش در پهنه مورد نظر کاهش مییابد. فصل تابستان و در مرتبه بعدی فصل بهار بیشترین مقدار تابش خورشیدی را به خود اختصاص دادهاند. . با مشاهده نقشههای فصل پاییز مشاهده میشود که همچنان بیشینه تابش دریافتی در جنوبشرق پهنه، ولی دامنه آن به طرف مرکز در حال گسترش است و کمینه تابش دریافتی درهمچنان در شمال پهنه، متمرکز شده است. کمترین مقدار تابش در ماه آبان با مقدار57/9 مگاژول و بیشترین مقدار تابش در ماه مهر با مقدار 49/19 مگاژول را به خود اختصاص دادهاند. در فصل زمستان بیشینه تابش دریافتی در جنوبشرق پهنه، در ایستگاه نهبندان و کمینه تابش دریافتی در شمال پهنه، در ایستگاههای قوچان و بجنورد متمرکز شده است. در بین ماههای فصل زمستان کمینه تابش در ماه دی با 98/7
مگاژول و در ماه اسفند بیشینه تابش با 08/18 مگاژول را به خود اختصاص داده است. در فصل بهار بیشینه تابش در جنوبشرق پهنه، و درماههای اردیبهشت و خرداد تا مرکز پهنه و ایستگاههای بیرجند و قائن گسترش یافته اند. و کمترین مقدار تابش همچنان در شمال پهنه، قرار دارد. کمترین مقدار تابش در ماه فروردین با مقدار 60/18 مگاژول و بیشینه تابش در ماه خرداد با مقدار 81/28 مگاژول متمرکز شده است. در فصل تابستان کمینه تابش در ایستگاه بجنورد در شمال پهنه محدود شده است. و بیشینه تابش از جنوبشرق پهنه جابجا شده است به طوریکه در ماههای مرداد و تیر بیشترین تابش در شرق و مرکز پهنه متمرکز شده است. کمترین تابش دریافتی در ماه شهریور با مقدار 54/17 مگاژول و بیشترین تابش دریافتی در ماه تیر با مقدار 24/30 مگاژول را به خود اختصاص دادهاند. جدول3 نتایج محاسبات تابش خورشیدی را نشان میدهد. و ملاحظه میشود ماه دی کمترین میانگین تابش خورشیدی دریافتی با مقدار 8/9 مگاژول برمتر مربع در پهنه مطالعاتی دارد و بیشینه میانگین تابش دریافتی به ترتیب در ماههای تیر، خرداد ، مرداد و اردیبهشت با مقادیر 8/28، 8/27، 2/27، 4/24 مگاژول بر متر مربع متمرکز شده است. همچنین با توجه به نتایج جدول 3 مشاهده میشود کمترین نوسانات در طول پهنه به ترتیب در ماههای مهر، خرداد و اردیبهشت با مقادیر 5/0، 8/0، 0/1 مگاژول بر متر مربع و بیشترین نوسات در طول پهنه در ماه آبان به مقدار 4/1 مگاژول و در ماههای آذر، اسفند و بهمن با مقدار 3/1 مگا ژول را به خود اختصاص دادهاند.
جدول3-آمار توصیفی متوسط ماهانه تابش خورشیدی طی سالهای (2016-1988 (برحسب مگاژول بر متر مربع در روز
آذر |
آبان |
مهر |
شهریور |
مرداد |
تیر |
خرداد |
اردیبهشت |
فروردین |
اسفند |
بهمن |
دی |
ایستگاه |
3/9 |
6/11 |
4/18 |
6/19 |
2/27 |
8/28 |
8/27 |
4/24 |
4/20 |
8/15 |
7/12 |
8/9 |
میانگین |
9/11 |
3/14 |
5/19 |
6/21 |
9/28 |
3/30 |
8/28 |
7/25 |
3/22 |
1/18 |
0/15 |
1/12 |
بیشینه |
5/7 |
7/9 |
7/17 |
7/17 |
2/25 |
3/26 |
0/26 |
7/22 |
8/18 |
9/13 |
8/10 |
1/8 |
کمینه |
3/1 |
4/1 |
5/0 |
2/1 |
1/1 |
1/1 |
8/0 |
0/1 |
1/1 |
3/1 |
3/1 |
2/1 |
انحرافمعیار |
7/1 |
9/1 |
3/0 |
4/1 |
2/1 |
1/1 |
6/0 |
1/1 |
1/1 |
8/1 |
7/1 |
5/1 |
واریانس |
|
|
شکل 2- مقدار تابش دریافتی سطح زمین در خراسان رضوی، شمالی، جنوبی (2016-1988) با مدل آنگستروم بر حسب مگاژول بر متر مربع |
|
|
ادامه شکل 2- مقدار تابش دریافتی سطح زمین در خراسان رضوی، شمالی، جنوبی (2016-1988) با مدل آنگستروم بر حسب مگاژول بر متر مربع
بحث و نتیجه گیری
تابش خورشیدی مهمترین عنصر هواشناسی است که به صورت مستقیم و غیر مستقیم بر تمام فرآیندهای آب و هوایی اثر میگذارد . در بسیاری از ایستگاههای همدید پهنه مورد نظر، دستگاه پیرانومتر که مستقیما مقدار تابش خورشیدی را اندازهگیری میکند وجود نداشت. در این پژوهش برای برآورد تابش خورشیدی از روش تجربی آنگستروم – پرسکات استفاده شد و ضرایب aو b با استفاده از مدل فرر محاسبه گردید. با توجه به مدل مورد استفاده و نقشههای حاصل از مدل درون یابی کریجینگ، شکل 2 مشخص شد که به ترتیب از شمال به جنوب بر میزان تابش دریافتی سطح گستره مورد مطالعه افزوده میشود، کمینه تابش خورشیدی در گستره شمالی و بیشینه آن در بخش جنوبی گستره مطالعاتی که ایستگاه نهبندان نماینده آن است مشاهده گردید. به میزان نزدیک شدن به دوره گرم سال بر میزان تابش خورشیدی گستره مطالعاتی افزوده میگردد به همان نسبت منطقه بیشینه تابش خورشیدی در سمت شرق و مرکز گستره مطالعاتی و در محدوده ایستگاههای قائن، بیرجند و تربت حیدریه متمرکز شده ولی کماکان کمینه ساعت آفتابی در شمال گستره و محدوده ایستگاه بجنورد و بیشینه آن در ایستگاه نهبندان متمرکز است. ماه های تیر، خرداد و مرداد بیشینه میزان تابش خورشیدی را نشان دادند.
از آنجا که در جنوب منطقه مورد مطالعه ارتفاع خورشید بیشتر است، لذا تابش بیشتری نسبت به نواحی شمالی دریافت می دارد. ماه آذر و دی نیز کمینه میزان تابش خورشیدی را در پهنه مورد مطالعه نشان داده است. اکثر ایستگاهها نوسانات تابش خورشیدی زیادی مشاهده میشود. نتایج بدست آمده حاکی از آن است که در مناطق مورد بررسی ساعات آفتابی و عرض جغرافیایی نقش زیادی در تابش بیشتر مناطق مورد نظر داشته است و مناطقی که ساعات آفتابی کمتری و هوای سردی داشتهاند تابش کمتری را به خود اختصاص داده است، لذا بهتر است مدل های مربوطه را توسعه داد تا از متغیر های اقلیمی بیشتری جهت برآورد تابش استفاده گردد.
منابع
1- کارشناس ارشد اقلیمشناسی، دانشگاه زنجان، ایران. *(مسوول مکاتبات)
2- دانشیار، گروه جغرافیا، دانشگاه پیام نور، تهران، ایران.
1- MSc in Climatology, University of Zanjan, Iran. *(Corresponding Author)
2-Associate Professor, Department of Geography, Payame Noor University, Tehran, Iran.
1-Angstrom
2-S.L. Lee
3-Lópeza et al
4-MODIS
5-Mghouchi et al
6-Zhang et al
7-Rabee et al
1-Yıldırım et al
2-Angstrom- Prescott
3-Achour et al
4-Rodriguez et al
[16]-MBE
[17]-RMSE
1-Pyranometer
2-Pyridomometer
3-Kriging
فصلنامه انسان و محیط زیست، شماره 47، زمستان 97
واسنجی تابش دریافتی سطح زمین در پهنه شمالشرق ایران
مهناز رستمیان [1] *
امیرحسین حلبیان[2]
تاریخ دریافت: 29/07/1396 |
تاریخ پذیرش: 16/12/1396 |
چکیده
تابش خورشیدی رسیده به سطح زمین (Rs)، در بسیاری از زمینههای مهندسی کشاورزی، آبیاری و هیدرولوژی کاربرد گستردهای دارد و به این دلیل روشهای مختلفی برای تخمین آن ارایه شده است که معادله آنگستروم – پرسکات یکی از مهمترین معادلات موجود در این زمینه است. این معادله دارای ضرایب تجربی است که برای هر موقعیت متفاوت میباشند. به دلیل اهمیت Rs، روشهای مختلفی برای تخمین آن توسعه داده شد. پژوهش حاضر با استفاده از دادههای میانگین روزانه ساعت آفتابی ایستگاههای سینوپتیک وابسته به اداره کل هواشناسی کشور در گستره خراسانشمالی، خراسانرضوی و خراستان جنوبی، طی بازهی 2016-1988 به برآورد تابش رسیده به سطح زمین پرداخته است. جهت برآورد تابش خورشیدی در سطح زمین از روش انگستروم_ پرسکات بهره گرفته شد. در این پژوهش با استفاده از روش فرر ضرایب a و b محاسبه و سپس در معادله آنگسترم جایگذاری گردید، در نهایت با استفاده از مدل درونیابی کریجینگ نقشههای همارزش تابش گستره مورد مطالعه ترسیم و تحلیل گردید. نتایج این پژوهش نشان داد که مناطق جنوب، جنوب شرق منطقه مورد مطالعه، بهخصوص ایستگاه نهبندان پتانسیل دریافت تابش خورشیدی، بیشتری دارند. کمینه مقادیر تابش دریافتی ماهانهی در ماههای آذر و دی در شمال و شمالغرب منطقه مورد مطالعه قرار دارد، که به ترتیب مقادیر 36/7، 98/7 مگاژول بر متر مربع میباشد. بیشینه مقدار تابش دریافتی، در ماه تیر با مقدار 24/30 مگاژول بر مترمربع اتفاق میافتد و مناطق جنوب و جنوبشرق تا مرکز پهنه بهخصوص ایستگاه نهبندان و به ترتیب فصول تابستان و بهار از نظر دریافت تابش خورشیدی، پتانسیل بیشتری دارند.
کلمات کلیدی: تابش خورشیدی، ساعت آفتابی، مدل فرر، آنگستروم- پرسکات، شمالشرق ایران.
|
Calibration of earth surface received radiation in north east of Iran
Mahnaz Rostamian [3]*(Corresponding Author)
Amir Hossein Halabian [4]
Abstract
The solar radiation reached the earth surface (RS) has wide applications in many areas of agriculture engineering, irrigation and hydrology and that’s why various methods have been represented for its estimation, which Angstrom – Prescott equation is one of the most important equation in this field. This equation has experimental coefficients that are different for any situation. Because of the importance of Rs, various methods have developed for its estimation. The radiation reached to earth surface is evaluated by using average daily sunny hours data of synoptic affiliated to the general state metrological office in the north of Khorasan, Khorasan Razavi and southern Khorasan 1988- 2016. The Angstrom – prescient method has been used for estimating solar radiation of earth surface in the studied area. In this research, the a & b coefficients were calculated by using Ferrer's method fastly and then have been placed in Angstrom equation. Finally, the radiation equivalent map of the studied area has been drawn and analyzed by using kriging interpolation model. The research outcome indicated that southern, south-east regions and at the next level, the center of the studied area especially Nehbndan station have more potential for receiving solar radiation. The minimum monthly received radiation of area during January and December month in the north and northwest of this region, which are 7.98 and 7.36 megajoules per square meter respectively. The maximum received radiation of area occurs in July about 30.24 megajoules. South and south-east regions up to the center of the of the area Nehbandan station have more potential for receiving in terms of radiation during summer & spring seasons respectively.
Key Words: Solar Radiation, Sundial, Ferrer Model, Angstrom- Prescott, North East of Iran.
مقدمه
انرژی تابشی خورشید سرمنشا انرژیهای روی کره زمین است که یکی از منابع مهم انرژی پاک به ویژه در کشور ما است. انرژی خورشید از معروفترین انرژیهای پاک و انرژیهای غیر فسیلی است و استفاده از آن کمترین آلودگی را در محیط زیست به وجود میآورد. دریافت اطلاعات انرژی خورشید به دو روش امکانپذیر است. روش اول اندازهگیری مستقیم انرژی خورشید توسط دستگاههای اندازهگیری مانند پیرانومتر و پیرهیلومتر میباشند. در روش دوم که با استفاده از روشهای تجربی به برآورد تابش میپردازد. مهمترین پارامتر که در بیشتر معادلات تجربی برآورد تابش خورشیدی از آن استفاده میشود، ساعات آفتابی است. تحقیقات نشان داده است که ساعات آفتابی مهمترین عامل تعیین کنندهی میزان تابش خورشیدی است. هدف از پژوهش حاضر واسنجی و برآورد مقدار انرژی دریافتی در سطح زمین در سه استان خراسانشمالی، رضوی و جنوبی، با استفاده از روش تجربی آنگستروم[5] و نیز با کمک داده های ساعت آفتابی و عرض جغرافیایی می باشد، از آنجا که هزینه دستگاههای اندازه گیری بالا بوده؛ با استفاده از روشهای فیزیکی و محاسبات سعی می شود تا مقدار تابش با کمترین مقدار خطا برآورد گردد.
پیشینه پژوهش
در ایران و جهان تاکنون پژوهشهای زیادی توسط محققان مختلف در زمینه برآورد انرژی خورشیدی با روشهای مختلفی صورت گرفته است. از جمله این مطالعات می توان به موارد زیر اشاره نمود: چو و لی[6] با استفاده از روش انگستروم به برآورد تابش خورشیدی در مالزی پرداختند. نتایج مطالعه نشان داد که در طول یک ماه همبستگی خطی وابستگی رخداد فرینهای بیشینه و کمینه روزانه طول ساعت آفتابی بسیار ضعیف است. برآوردهای بدست آمده برای شهرهای نزدیک تایلند و سنگاپور موافقت بیشتری را با دادههای این کشورها نشان دادند (1). لوپز و همکاران[7] به برآورد تابش خورشیدی از دادههای مودیس[8] در جنوب اسپانیا پرداختند. آنها یک پدیده دیگر به مدل پارامتری برای تخمین تابش خورشیدی تحت شرایط ابرناکی اضافه نمودند. محتوای بخار آب قابل بارش جو تنها ورودی تجربی مدل میباشد (2). ماکوچی و همکاران[9] مدلی جدید برای برآورد و ارزیابی تابش خورشیدی در مراکش شمالی ارایه دادند نتایجشان نشان میدهد که مدل پیشنهادی برای برآورد تابش خورشیدی در تمام فصلهای سال برای موقعیت مورد مطالعه و برای روز مورد نظر با استفاده از ورودیهای ارتفاع (متر)، طول جغرافیایی (درجه) و عرض جغرافیایی (درجه) مورد استفاده قرار می گیرد (3). ژانک و همکاران[10] به مرور انتقادی بر مدلهای مورد استفاده برآورد تابش خورشیدی طی یک دهه قبل پرداختند. در این مطالعه روشهای مورد استفاده در مقالات مرتبط با تابش خورشید مورد بررسی قرار گرفتهاند. نتایج نشان میدهد که مدلهای کسری زمانی تابش آفتاب و شبکههای عصبی مصنوعی عملکرد مشابهی دارند که برای برآورد ماهانه متوسط رزوانه تابش جهانی استفاده میشود. در حالی که بیشتر کارها برای مطالعه روش برآورد به فواصل زمانی کوتاهتر و مکانیزمهای اتمسفر برای تابش خورشید نیاز دارند (4). ربیع و همکاران[11] با استفاده از دادههای چند سال اخیر و مدل شبکه های عصبی مصنوعی برای برآورد متوسط روزانه تابش خورشیدی در کویت پرداختند. مدل شبکههای عصبی مصنوعی قادر بوده با استفاده از روش غیرخطی به پیشبینی بپردازد. توسعه پیشبینی وابسته به این است که دیدگاههای مقامات کشوری، سرمایهگذاران و نیروی مهندسی سیستم، و تعیین مکانهای مناسب برای تاسیسات خورشید و شایستگی فنی-اقتصادی با یکپارچه سازی انرژی خورشیدی در مقیاس طولانی مدت باشد (5). یالدریم و همکاران[12] به برآورد روزانه تابش جهانی خورشید با استفاده از رابطه کاربری گرافیکی در منطقه مدیترانه شرقی ترکیه در یک دوره 11 ساله پرداختند. آنها روشهای مختلف مدل های مختلف تابش خورشیدی را تحلیل نمودند. همچنین یک مدل بهینه شده بر اساس شبکه عصبی مصنوعی و مدل آنگستروم - پرسکات برای برآورد تابش خورشیدی جهانی روزانه ارایه دادند. و یک مدل بهینه شده بر اساس شبکه عصبی مصنوعی و مدل آنگستروم - پرسکات [13]برای برآورد تابش خورشیدی جهانی روزانه ارایه دادند (6) . آچور[14] مدلی ترکیبی برای برآورد تابش خورشیدی ماهانه در جنوب الجزایر برای دوره 5 ساله ارایه دادند. نتایج به دست آمده آنها نشان داد که مدل جدید هیبرید میتواند قادر به پیش بینی تابش خورشیدی جهانی با دقت عالی در این مکان باشد (7). رودریگو و همکاران[15] با استفاده از تصاویر ماهوارهای و الگوریتم ژنتیکی به برآورد تابش جهانی خورشید پرداختند. آنها از تصاویر با یاختههای 3 کیلومتر و در بازهی ژانویه تا دسامبر 2008 بهره بردند. برآیند مطالعه ضریب همبستگی 99 درصد را نشان داد. برآوردهای تولید شده نسبتا دقیق بوده و خطای فضایی نسبتا یکنواختی را نشان دادند. مدل فوق حتی برای روزهای ابری نیز برآوردهای قابل اعتمادی را نشان داده است (8). خلیلی و رضایی صدر به برآورد تابش کلی خورشید در گستره ایران بر مبنای دادههای اقلیمی در دوره (1991-1966) باروش انگستروم پرداختند و نتایج آنها نشان داد که به جز ایستگاههای جاسک و زاهدان که آنالیزها منجر به رابطه معنیداری برای آنها نشده است، در سایر ایستگاهها روابط در سطح 01/0 معنی دار میباشد. ضریب A معادله آنگستروم در شبکه تابشسنجی ایران بین 21/0 در یزد تا 39/0 در کرمانشاه متفاوت میباشد (9). صفائی و همکاران در تحقیقی به برآورد پتانسیل تابش خورشیدی در ایران و تهیه اطلس تابشی آن در یک دوره 5 تا 11 ساله با استفاده از سه مدل ریاضی: آنگستروم، روش حداکثر احتمال و مدل هیبرید اقدام کردند و مدل هیبرید بهترین برازش را با این دادهها داشت و حاصل کار آنها 12 نقشه مربوط به میانگین روزانه تابش کلی خورشیدی در هر ماه است (10).علیزاده و خلیلی در پژوهشی به تعیین ضرایب معادله آنگستروم و توسعه یک معادله رگرسیونی برآورد تابش خورشیدی در مشهد در سالهای 1380 الی 1378 ضرایب رابطه آنگسترم – پرسکات واسنجی پرداختند و همچنین یک معادله رگرسیونی منطقهای ارائه دادند نتایج آنالیز آماری روابطی که بررسی کردند نشان داد که دقت معادله منطقه ای چند متغیره با رابطه آنگستروم – پرسکات اختلاف معنی داری نداشته است (11). یزدان پناه و همکاران در تحقیقی برآورد تابش کلی خورشید در سطح افقی زمین در اصفهان با استفاده از آمار هواشناسی ایستگاه ازن سنجی اصفهان در دوره آماری 1384 تا 1387 نسبت به یافتن مدل مناسب برآورد شدت تابش خورشیدی رسیده به سطح از روی آمار ساعات آفتابی با استفاده از مدل آنگستروم اقدام کردند ونتایج این تحقیق نشان داد مدل تهیه شده بهویژه برای ماههای می و آگوست دقت بسیار بالایی دارد. و مقدار ضریب رگرسیون مدل ساخته شده ما معادل 91/0 است، ضرایب رگرسیون برای مدلهای، فرر و ام سی کلوچ 6/0 و 7/0 است (12). معینی و همکاران در تحقیقی به برآورد تابش خورشیدی در ایران با استفاده از یک مدل بهینه پرداختند و از مدل آنگستروم برای برآورد میزان تابش خورشید روی سطح افقی بهره گرفتند و تابش کل را برای نقاط مختلف ایران وهمچنین محاسبه طول روز و تابش فراجو را برآورد کردند (13). سبزیپرور و بیاتورکشی به ارزیابی دقت روشهای شبکه عصبی مصنوعی و عصبی – فازی در شبیهسازی تابش کل خورشیدی طی سالهای 1992 تا 2006 در ایستگاههای کرمان، ارومیه، شیراز و اصفهان پرداختند و و نتیجه گرفتند روابط تجربی مورد استفادشان در تخمین تابش خورشیدی موید ارتباط مناسب بین پارامترهای هواشناسی و مقدار تابش خورشیدی میباشد و بهترین مدل شبکه عصبی مصنوعی با تابع فعالیت تانژانت هایپربولیک در ایستگاه کرمان میباشد (14). موسویبایگی و اشرف در تحقیقی به شناسایی مناطق کمترین میزان ابرناکی به منظور پهنهبندی نواحی پرتابش کشور در دوره 1387-1368 با روش درون یابی اسپیلاین spline)) اقدام نمودند و نتایج آنها نشان داد که بیشترین روزهای آفتابی در فروردین ماه متعلق به شهرستان بافت در اردیبهشت ماه متعلق به شهرستان خاش بوده است و همچنین اغلب مناطق مناطق ایران از استعداد و قابلیت بالایی برای استفاده از انرزی تجدیدپذیر خورشیدی برخوردارند (15). پیری و همکاران در پژوهشی به مدلسازی تابش خورشید رسیده به زمین با استفاده از ANFIS و مدلهای تجربی ایستگاههای زاهدان و بجنورد در دوره آماری به ترتیب 5 تا 7 ساله پرداختند و دو روش تجربی را جهت برآورد تابش روزانه خورشید واسنجی و ارزیابی کردند و نتیجه گرفتند که مدل نرو فازی برآورد بهتری نسبت به دو مدل تجربی در برآورد تابش دارد (16). زمانی محی آبادی در تحقیقی به پیش بینی لحظهای تابش کل خورشید در شهر رفسنجان توسط شبکه عصبی برای یک سال پرداخت و نتیجه گرفت که شاخصهای آمار میانگین خطای بایاس [16]و ریشه میانگین مربعات خطا [17]و t برای شبکه عصبی بترتیب 0037/0، 5366/36 و 0232/0 میباشند که نمایانگر عملکرد خوب مدل طراحی شده در شبکه عصبی است (17). عرفانیان و بابایی حصار در پژوهشی به ارزیابی مدل هیبرید در تخمین تابش خورشیدی روزانه در تعدادی از ایستگاههای تابشسنجی ایران پرداختند. و مقدار کم معیارهای خطا برای مدل هیبرید بیانگر کارایی نسبتا بالای این مدل در تخمین تابش خورشیدی ایستگاههاس انتخابی در مقیاس روزانه میباشد (18). هوشنگی و همکاران در تحقیقی به بررسی منطقهای پتانسیل تابش خورشیدی با ارزیابی و بهینه سازی روشهای درون یابی در سطح کشور ایران و نتایج حاصل با اعتبار سنجی متقابل مقایسه شدند. و بررسی نتایجشان نشان داد که بهینه سازی عوامل موثر در روشهای درون یابی موجب افزایش همگرایی 3/8 درصدی مقادیر نتایج میگردد (19).کرمپور و همکاران در پژوهشی توازن تابش در سطح شهر یزد در ماه ژانویه 2006 تا 2008 و میانگین ماهانه ساعتی سپیدایی، شار تابش خالص در شرایط هوای صاف و ابری واکاوی کردند و به این نتیجه رسیدند که میزان سپیدایی غرب شهر به دلیل ساخت شهری بیشتر از شرق شهر است. و مقدار آن در روزهای که سطح شهر مرطوب بوده تا 2 درصد کاهش پیدا میکند(20). مجرد و همکاران با به کارگیری مدل بهینه شدهی برد و هول استورم در دورهی زمانی 2009-1990 اقدام به برآورد تابش خورشیدی دریافتی سطح زمین در استان کرمانشاه نمودند. نتایج تحقیق آنها نشان داد بین مقادیر تابش اندازه گیری شده با پیرانومتر و مقادیر برآورد شده توسط مدل اختلاف کمی وجود دارد که نشان دهندهی توانایی مدل در برآورد تابش است. و در فصول پاییز و زمستان به دلیل نقش تابش پراکنده، اختلاف تابش در مناطق مختلف استان کاهش مییابد (21). فلاح قالهری و شاکری با به کارگیری مدلهای آنگستروم – پرسکات به واسنجی ضرایب رابطه آنگستروم – پرسکات در ایستگاههای منتخب استان خراسان رضوی نمودند. نتایج آنها نشان داد که بین مقادیر تابش اندازهگیری شده با پیرانومتر و مقادیر تابش برآورد شده توسط رابطه بلکه اختلاف کمی وجود دارد. و همچنین نتایج حاصل از پهنهبندی تابش نشان داد مناطق جنوبی استان از نظر دریافت تابش خورشیدی، پتانسیل بیشتری دارند (22).
مواد و روش
منطقه مورد مطالعه
منطقه مورد مطالعه در این تحقیق، پهنهی شمالشرق متشکل از استانهای خراسانجنوبی، خراسان رضوی وخراسان شمالی واقع در طول جغرافیایی از 55 درجه و 55 دقیقه و53 ثانیه غربی تا 61 درجه و 17 دقیقه و 2 ثانیه شرقی و عرض جغرافیایی از 30 درجه و 29 دقیقه و 33 ثانیه جنوبی تا 38 درجه و 18 دقیقه و 3 ثانیه جنوبی است مساحت پهنه 5/ 239730 کیلومتر مربع میباشد(شکل 1). این استان ها تحت تاثیر توده هوای سرد قطبی، سیبری ، توده هوای مدیترانهایی و توده هوای گرم و خشک (پرفشار جنب حاره ) قرار دارند که شرایط اقلیمی متفاوتی را برای این استانها رقم می زنند.
شکل1-موقعیت جغرافیایی منطقه مورد مطالعه به همراه ایستگاهای منتخب
همانطور که در جدول 1 مشاهده میشود کمترین میانگین مجموع ساعات آفتابی به ترتیب ایستگاه بجنورد در استان خراسانشمالی با 6/7 و در مرتبه بعدی ایستگاههای قوچان، گلمکان، سرخس، مشهد در استان خراسانرضوی با مقادیر 7/7، 0/8،0/8، 1/8 دارا میباشند و بیشترین میانگین مجموع ساعات آفتابی در جنوبشرق و مرکز پهنه به ترتیب در ایستگاههای نهبندان و قائن و بیرجند و فرودس با مقادیر 2/9، 1/9 به خود اختصاص دادهاند. انحراف معیار میزان پراکندگی دادهها را نسبت به میانگین میسنجد.کمترین نوسانات مجموع ساعات آفتابی در طول پهنه در ایستگاههای نهبندان و بیرجند با مقادیر (7/1، 9/1) و بیشترین نوسانات در ایستگاههای سرخس و مشهد و گلمکان به ترتیب با مقادیر (8/2، 7/2، 5/2) مشاهده شد.
جدول 1- آمار توصیفی ساعات آفتابی ایستگاههای سینوپتیک منطقه مورد مطالعه طی سالهای (2016-1988)
بیشینه |
کمینه |
واریانس یا پراش |
انحراف معیار |
میانگین |
ایستگاه |
استان |
8/11 |
3/5 |
8/4 |
2/2 |
3/8 |
بشرویه |
خراسان جنوبی
|
9/11 |
6/6 |
6/3 |
9/1 |
1/9 |
بیرچند |
|
1/12 |
0/6 |
2/4 |
1/2 |
1/9 |
فردوس |
|
3/12 |
2/6 |
6/4 |
1/2 |
1/9 |
قائن |
|
8/11 |
8/6 |
8/2 |
7/1 |
2/9 |
نهبندان |
|
3/12 |
8/5 |
5/5 |
3/2 |
9/8 |
تربت حیدریه |
خراسان رضوی
|
5/11 |
6/5 |
7/4 |
2/2 |
5/8 |
سبزوار |
|
0/12 |
3/4 |
9/7 |
8/2 |
0/8 |
سرخس |
|
3/11 |
6/4 |
0/6 |
4/2 |
7/7 |
قوچان |
|
9/11 |
8/5 |
7/4 |
2/2 |
8/8 |
کاشمر |
|
6/11 |
0/5 |
2/6 |
5/2 |
0/8 |
گلمکان |
|
0/12 |
0/6 |
3/4 |
1/2 |
8/8 |
گناباد |
|
1/12 |
6/4 |
3/7 |
7/2 |
1/8 |
مشهد |
|
9/10 |
6/4 |
5/5 |
3/2 |
6/7 |
بجنورد |
خراسان شمالی |
روش تحقیق
در این تحقیق برای برآورد میزان تابش دریافتی سطح زمین در استانهای خراسانرضوی، خراسانجنوبی و خراسان شمالی از میانگین دادههای روزانه ساعات آفتابی 13 ایستگاه سینوپتیک منطقه که دارای آمار طولانی، طی یک بازهی زمانی 28 ساله (2016-1988) بودند از سازمان هواشناسی کشور دریافت شد و سپس استفاده گردید، البته اندازهگیری مستقیم انرژی خورشید توسط دستگاههای اندازهگیری مانند پیرانومتر [18] پیرهیلومتر[19] میباشد. در منطقه مورد بررسی، اکثر ایستگاهها مجهز به دستگاههای اندازهگیری تابش خورشید نبودند، یا اندازهگیری تابش در آن نقاط به صورت پراکنده میباشد لذا به برآورد مقدار تابش دریافتی سطح زمین در شمالشرق ایران پرداخته شد، در این تحقیق ابتدا با استفاده از مدل فرر ضرایب a و b محاسبه گردید سپس با جایگذاری این ضرایب در معادله آنگستروم – پرسکات انرژی خورشید برآورد گردید. (روابط 6-1). مقدار تابش دریافتی سطح زمین در خراسانهای رضوی،شمالی، جنوبی برای12 ماه براساس روش کریجینگ[20] نقشه آنها در Arc GIS ترسیم شد شکل (2).
مدل انگستروم :
پریسکات از طریق رابطه زیر برآورد می گردد (21):
(1) |
|
در رابطه فوق بیانگر میانگین تابش کل روزانه درهر ماه میباشد.
(2)
|
بیانگر تابش اندازهگیری شده در خارج از جو میباشد مقدار از رابطه (2) بدست میآید: که در آن شماره ژوولیوسی ، ثابت خورشید که معادل 1367، عرض جغرافیایی محل و زاویه میل که از رابطه (3) زیر به دست میآید: ( 12).
(3) |
که در آن ( زاویه ساعتی خورشید در هنگام غروب) از رابطه زیر به دست میآید:
(4) |
بیانگر میانگین ماهیانه ساعات آفتابی روزانه و میانگین ماهیانه ماکزیمم ساعات آفتابی (طول روز) میباشد و مقدار آن از رابطه (5) محاسبه میشود.
(5) |
از مدل فرر برای بدست آوردن مقادیر (a,b) در برآورد تابش کلی ( درسطح افقی) استفاده شده است( رابطه6). مدل فرر به شرح زیر است: (12).
(6) |
a= -0.27+1.75 b= 1.32-2.90 |
یافته ها
همان طور که در جدول 2 مشاهده میشود، بیشترین میانگین تابش خورشیدی ماهانه در جنوب شرق پهنه مورد نظر در ایستگاه نهبندان بامقدار 4/20 و در مرتبه های بعدی ایستگاه های بیرجند و قائن ، فردوس ، گناباد و تربت حیدریه به ترتیب با مقادیر 9/19، 9/19، 7/19، 3/19، 3/19 مگاژول بر متر مربع را به خود اختصاص دادهاند. اکثر ایستگاهها نوسانات تابش خورشیدی درون سالانه نسبتا زیادی را نشان دادند. ایستگاه سرخس با انحراف معیار 5/7 بیشینه نوسانات را نشان داده است. و کمینه نوسانات در ایستگاه نهبندان و با انحراف معیار 6 مشاهده شده است. در این بین ایستگاه های سرخس و بجنورد حداقل متوسط تابش خورشیدی دریافتی را با حدود 5/7 و 6/7 مگاژول بر متر مربع را نشان دادند. بیشینه متوسط روزانه تابش را ایستگاه قائن و تربت حیدریه به ترتیب با حدود 3/30 و 30 مگاژول بر متر مربع را تجربه کرده است.
جدول2- آمارتوصیفی تابش خورشیدی ماهانه (2016-1988) برحسب مگاژول بر متر مربع در روز
ایستگاه |
میانگین |
انحراف معیار |
واریانس |
کمینه |
بیشینه |
بیرجند |
9/19 |
3/6 |
3/40 |
0/11 |
3/29 |
بشرویه |
9/18 |
7/6 |
0/45 |
2/9 |
9/28 |
نهبندان |
4/20 |
0/6 |
0/36 |
9/11 |
1/29 |
قائن |
9/19 |
9/6 |
9/47 |
3/10 |
3/30 |
فردوس |
7/19 |
7/6 |
6/44 |
3/10 |
6/29 |
قوچان |
3/17 |
9/6 |
2/47 |
8/7 |
2/27 |
گلمکان |
9/17 |
0/7 |
0/49 |
2/8 |
1/28 |
گناباد |
3/19 |
6/6 |
7/43 |
0/10 |
3/29 |
کاشمر |
1/19 |
8/6 |
3/46 |
5/9 |
9/28 |
مشهد |
0/18 |
4/7 |
8/54 |
1/8 |
2/29 |
سبزوار |
5/18 |
7/6 |
5/45 |
9/8 |
8/27 |
سرخس |
9/17 |
5/7 |
9/56 |
5/7 |
0/29 |
تربت حیدریه |
3/19 |
2/7 |
6/51 |
4/9 |
0/30 |
بجنورد |
1/17 |
7/6 |
7/44 |
6/7 |
3/26 |
باتوجه به نقشههای شکل 2 مشاهده میشود کمترین مقادیر تابش دریافتی ماهانهی پهنه در ماههای دی و آذر در شمال و شمال غرب پهنه رخ میدهد که به ترتیب مقادیر 98/7، 36/7 مگاژول بر متر مربع میباشد. بیشترین مقدار تابش دریافتی پهنه در ماه تیر با مقدار 24/30 مگاژول بر مترمربع اتفاق میافتد ودر مرتبه های بعدی به ترتیب ماههای خرداد، مرداد، اردیبهشت و فروردین 81/28، 78/27، 72/25 و 24/22 مگازول بر متر مربع در جنوبشرق و مرکز پهنه مورد نظررا به خود اختصاص دادهاند. افزایش مقادیر تابش در ماههای گرم سال، در مرحله نخست به خاطر افزایش زاویه خورشید از فصل سرد به گرم میشود. از طرفی دیگر در ماههای سرد سال ابرناکی آسمان بیشتر است و باعث کاهش دریافتی خواهد بود در فصل زمستان به دلیل نقش تابش پراکنده، اختلاف تابش در پهنه مورد نظر کاهش مییابد. فصل تابستان و در مرتبه بعدی فصل بهار بیشترین مقدار تابش خورشیدی را به خود اختصاص دادهاند. . با مشاهده نقشههای فصل پاییز مشاهده میشود که همچنان بیشینه تابش دریافتی در جنوبشرق پهنه، ولی دامنه آن به طرف مرکز در حال گسترش است و کمینه تابش دریافتی درهمچنان در شمال پهنه، متمرکز شده است. کمترین مقدار تابش در ماه آبان با مقدار57/9 مگاژول و بیشترین مقدار تابش در ماه مهر با مقدار 49/19 مگاژول را به خود اختصاص دادهاند. در فصل زمستان بیشینه تابش دریافتی در جنوبشرق پهنه، در ایستگاه نهبندان و کمینه تابش دریافتی در شمال پهنه، در ایستگاههای قوچان و بجنورد متمرکز شده است. در بین ماههای فصل زمستان کمینه تابش در ماه دی با 98/7
مگاژول و در ماه اسفند بیشینه تابش با 08/18 مگاژول را به خود اختصاص داده است. در فصل بهار بیشینه تابش در جنوبشرق پهنه، و درماههای اردیبهشت و خرداد تا مرکز پهنه و ایستگاههای بیرجند و قائن گسترش یافته اند. و کمترین مقدار تابش همچنان در شمال پهنه، قرار دارد. کمترین مقدار تابش در ماه فروردین با مقدار 60/18 مگاژول و بیشینه تابش در ماه خرداد با مقدار 81/28 مگاژول متمرکز شده است. در فصل تابستان کمینه تابش در ایستگاه بجنورد در شمال پهنه محدود شده است. و بیشینه تابش از جنوبشرق پهنه جابجا شده است به طوریکه در ماههای مرداد و تیر بیشترین تابش در شرق و مرکز پهنه متمرکز شده است. کمترین تابش دریافتی در ماه شهریور با مقدار 54/17 مگاژول و بیشترین تابش دریافتی در ماه تیر با مقدار 24/30 مگاژول را به خود اختصاص دادهاند. جدول3 نتایج محاسبات تابش خورشیدی را نشان میدهد. و ملاحظه میشود ماه دی کمترین میانگین تابش خورشیدی دریافتی با مقدار 8/9 مگاژول برمتر مربع در پهنه مطالعاتی دارد و بیشینه میانگین تابش دریافتی به ترتیب در ماههای تیر، خرداد ، مرداد و اردیبهشت با مقادیر 8/28، 8/27، 2/27، 4/24 مگاژول بر متر مربع متمرکز شده است. همچنین با توجه به نتایج جدول 3 مشاهده میشود کمترین نوسانات در طول پهنه به ترتیب در ماههای مهر، خرداد و اردیبهشت با مقادیر 5/0، 8/0، 0/1 مگاژول بر متر مربع و بیشترین نوسات در طول پهنه در ماه آبان به مقدار 4/1 مگاژول و در ماههای آذر، اسفند و بهمن با مقدار 3/1 مگا ژول را به خود اختصاص دادهاند.
جدول3-آمار توصیفی متوسط ماهانه تابش خورشیدی طی سالهای (2016-1988 (برحسب مگاژول بر متر مربع در روز
آذر |
آبان |
مهر |
شهریور |
مرداد |
تیر |
خرداد |
اردیبهشت |
فروردین |
اسفند |
بهمن |
دی |
ایستگاه |
3/9 |
6/11 |
4/18 |
6/19 |
2/27 |
8/28 |
8/27 |
4/24 |
4/20 |
8/15 |
7/12 |
8/9 |
میانگین |
9/11 |
3/14 |
5/19 |
6/21 |
9/28 |
3/30 |
8/28 |
7/25 |
3/22 |
1/18 |
0/15 |
1/12 |
بیشینه |
5/7 |
7/9 |
7/17 |
7/17 |
2/25 |
3/26 |
0/26 |
7/22 |
8/18 |
9/13 |
8/10 |
1/8 |
کمینه |
3/1 |
4/1 |
5/0 |
2/1 |
1/1 |
1/1 |
8/0 |
0/1 |
1/1 |
3/1 |
3/1 |
2/1 |
انحرافمعیار |
7/1 |
9/1 |
3/0 |
4/1 |
2/1 |
1/1 |
6/0 |
1/1 |
1/1 |
8/1 |
7/1 |
5/1 |
واریانس |
|
|
شکل 2- مقدار تابش دریافتی سطح زمین در خراسان رضوی، شمالی، جنوبی (2016-1988) با مدل آنگستروم بر حسب مگاژول بر متر مربع |
|
|
ادامه شکل 2- مقدار تابش دریافتی سطح زمین در خراسان رضوی، شمالی، جنوبی (2016-1988) با مدل آنگستروم بر حسب مگاژول بر متر مربع
بحث و نتیجه گیری
تابش خورشیدی مهمترین عنصر هواشناسی است که به صورت مستقیم و غیر مستقیم بر تمام فرآیندهای آب و هوایی اثر میگذارد . در بسیاری از ایستگاههای همدید پهنه مورد نظر، دستگاه پیرانومتر که مستقیما مقدار تابش خورشیدی را اندازهگیری میکند وجود نداشت. در این پژوهش برای برآورد تابش خورشیدی از روش تجربی آنگستروم – پرسکات استفاده شد و ضرایب aو b با استفاده از مدل فرر محاسبه گردید. با توجه به مدل مورد استفاده و نقشههای حاصل از مدل درون یابی کریجینگ، شکل 2 مشخص شد که به ترتیب از شمال به جنوب بر میزان تابش دریافتی سطح گستره مورد مطالعه افزوده میشود، کمینه تابش خورشیدی در گستره شمالی و بیشینه آن در بخش جنوبی گستره مطالعاتی که ایستگاه نهبندان نماینده آن است مشاهده گردید. به میزان نزدیک شدن به دوره گرم سال بر میزان تابش خورشیدی گستره مطالعاتی افزوده میگردد به همان نسبت منطقه بیشینه تابش خورشیدی در سمت شرق و مرکز گستره مطالعاتی و در محدوده ایستگاههای قائن، بیرجند و تربت حیدریه متمرکز شده ولی کماکان کمینه ساعت آفتابی در شمال گستره و محدوده ایستگاه بجنورد و بیشینه آن در ایستگاه نهبندان متمرکز است. ماه های تیر، خرداد و مرداد بیشینه میزان تابش خورشیدی را نشان دادند.
از آنجا که در جنوب منطقه مورد مطالعه ارتفاع خورشید بیشتر است، لذا تابش بیشتری نسبت به نواحی شمالی دریافت می دارد. ماه آذر و دی نیز کمینه میزان تابش خورشیدی را در پهنه مورد مطالعه نشان داده است. اکثر ایستگاهها نوسانات تابش خورشیدی زیادی مشاهده میشود. نتایج بدست آمده حاکی از آن است که در مناطق مورد بررسی ساعات آفتابی و عرض جغرافیایی نقش زیادی در تابش بیشتر مناطق مورد نظر داشته است و مناطقی که ساعات آفتابی کمتری و هوای سردی داشتهاند تابش کمتری را به خود اختصاص داده است، لذا بهتر است مدل های مربوطه را توسعه داد تا از متغیر های اقلیمی بیشتری جهت برآورد تابش استفاده گردد.
منابع
1- کارشناس ارشد اقلیمشناسی، دانشگاه زنجان، ایران. *(مسوول مکاتبات)
2- دانشیار، گروه جغرافیا، دانشگاه پیام نور، تهران، ایران.
1- MSc in Climatology, University of Zanjan, Iran. *(Corresponding Author)
2-Associate Professor, Department of Geography, Payame Noor University, Tehran, Iran.
1-Angstrom
2-S.L. Lee
3-Lópeza et al
4-MODIS
5-Mghouchi et al
6-Zhang et al
7-Rabee et al
1-Yıldırım et al
2-Angstrom- Prescott
3-Achour et al
4-Rodriguez et al
[16]-MBE
[17]-RMSE
1-Pyranometer
2-Pyridomometer
3-Kriging