بسم الله الرحمن الرحیم –
فرمت فایل :word (با قابلیت ویرایش) – تعداد صفحه : 23 صفحه -
مقدمه
فرسایش بادی یک معضل جدی در بیشتر مناطق خشک و نیمه خشک دنیا و ایران است . توانائی پیش بینی دقیق فرسایش بادی خاک برای بسیاری منظورها ، از جمله برنامه های حفاظتی ، منابع طبیعی ، و کاهش آلودگی هوا ناشی از طوفان ضروری است ( 3 ) . از آنجایی که نیروی باد در طول سال ، ماه وحتی روز تا حد زیادی تغییر می نماید ، و همچنین قدرت فرسایندگی باد بستگی به توان سوم سرعت باد دارد . به منشور پیش بینی و کنترل فرسایش بادی در هر منطقه توزیع سرعت باد حائز اهمیت میباشد . همچنین علاوه برسرعت باد ، دانستن چگونگی تغییرات جهت باد در منطقه نیز امری ضروری است . زیرا نسبت جهت باد به جهت اضلاع زمین ، بادشکن ها ، ردیف کاشت گیاهان ، و شخم زمین ، نقش مهمی را در پیش بینی مقدار و جهت فرسایش بادی ایفا می کند ( 4) . مدل های مختلفی برای نشان دادن توزیع سرعت باد استفاده شده است . بی شک توزیع ویبل یکی از گسترده ترین توزیع هایی است که تا بحال برای نشان دادن پراکندگی سرعت باد مورد استفاده قرار گرفته است ( 5 ) . اهداف این پژوهش عبارت بودند از : 1 ) شبیه سازی ساعتی سرعت و جهت باد به روش استوکاستیک با استفاده از توزیع ویبل ، به منظور استفاده در مدل WEPS ، برای پیش بینی فرسایش بادی در شهرهای مختلف ایران ، 2) آزمون اعتبار سنجی توزیع ویبل و مدل کامپیوتری windpred ، در پیش بینی ساعتی سرعت و جهت باد ، . 3) ترسیم نقشه های سرعت و جهت بادهای فرساینده در ایران .
مواد و روشها
ابتدا 38 شهر که دارای حداقل 10 سال آمار ساعتی سرعت و جهت باد بودند ، انتخاب گردیدند . در مرحله بعد تعداد سال آماری هر شهری به دو دوره برابر تقسیم گردید ، بطوریکه از دوره اول برای شبیه سازی و از دوره دوم برای آزمون اعتبار سنجی مدل ( با استفاده از معنی دار بودن و نبودن ضرائب همبستگی ) استفاده شد . سپس با استفاده از توزیع ویبل ، شبیه سازی سرعت و جه باد بصورت ساعتی توسط برنامه Windpred ( 2 و 1 ) انجام گرفت . تابع توزیع تجمعی ویبل F (U) به صورت زیر میباشد :
(1) [-(u/c)k]F(u)=1-exp
که در این معادله u سرعت باد ( متر بر ثانیه ) ، c، پارامتر مقیاس ( با واحد سرعت ) ، و k پارامتر شکل ( بدون واحد ) ، میباشند ( 6 ) . در هرمرحله بعد ، دوره های باد آرام حذف و فراوانی باد در هر گروه سرعتی نرمالیزه گردیدند . بنابراین :
(2) [-(u/c)k]= 1-exp [ (F(u)-F0 ) / (1-F0)] F1(u) =
که در آن F1(u) توزیع تجمعی در حالتی است که دوره های باد آرام حذف شده است ، و F0 فراوانی دوره های باد آرام میباشد . پارامترهای k, c به روش حداقل مربعات و بکارگیری تابع توزیع تجمعی محاسبه شدند ( معادله 2 ) .
با استفاده از پارامترهای توزیع ویبل (c,k) فراوانی سرعت باد در هر ماه و در سال بصورت تجمعی و نرمال شده بدست آمد . به منظور شبیه سازی جهت باد ، اعداد بین صفر و یک بصورت تصادفی انتخاب ، و با جدول توزیع تجمعی جهت باد مقایسه گردیدند . برای شبیه سازی سرعت باد براساس جهت باد تعیین شده ، پارامترهای c،k توزیع ویبل برای آن جهت خاص از جدول های تعیین شده قبلی ، بدست امد و از معادله زیر استفاده گردید :
تحقیق در موردپیش بینی سرعت و جهت بادهای فرساینده در ایران
