مدلسازی خطر جابهجایی زمین با رویکرد چند معیاره، یک خط تحقیق برای ارائه چارچوب روششناختی برای نقشهبرداری ریسک است. هدف از این کار ایجاد مدلی برای ترسیم «خطر جابجایی زمین» با رویکرد GIS بر اساس تحلیل چند معیاره است. روش ایجاد نقشه های موضوعی با ترکیب در GIS، عوامل تعیین کننده (شیب، سنگ شناسی، آب، …) در محرک پدیده زمین لغزش و انقباض – تورم خاک با استفاده از روش مدل جمع وزنی (WSM) است. سپس این نقشهها با هم ترکیب میشوند تا نقشهای از حرکات زمین را ارائه دهند. کاربرد این روش امکان توزیع مکانی معیارها و پدیده های مختلف در منطقه برهوم، منطقه حوضه هدنه، شرق الجزایر را فراهم می کند. در نتیجه، منطقه مورد مطالعه به چهار منطقه مختلف تقسیم شده است: 1) مناطق بدون خطر جابجایی زمین،
کلید واژه ها
حرکات زمینی ، خطر ، MCDA ، GIS ، برهوم
1. مقدمه
خطرات طبیعی با شدت کم و بیش فعالیت های انسان را تهدید می کند. الجزایر مکرراً تحت تأثیر بلایای طبیعی قرار گرفته است. این آخری ها جان افراد را به خطر می اندازند و خسارات غیرقابل برآوردی به بار می آورند.
مناطق شمالی الجزایر که با تراکم جمعیت و سرعت رشد شهری مشخص می شوند، شاهد مشکلاتی در ساخت و ساز ناشی از مشکلات مربوط به جابجایی زمین هستند. علت این امر ضعف شناخت پدیده هایی مانند زمین لغزش و چروکیدگی و متورم شدن خاک است.
حرکات زمین شامل تمام موارد میدانی جابجایی توده های زمین تحت اثر گرانش است. آنها می توانند بسیار کند یا بسیار سریع باشند.
حرکات زمین را می توان به دو نوع طبقه بندی کرد: حرکات شیب: حرکات اریب هستند که هم شیب های طبیعی و هم مصنوعی را تحت تاثیر قرار می دهند، فرو ریختن ها، لغزش ها و غیره، حرکات عمودی که شامل فرونشست، فروریختن، نشست و غیره است که ناشی از پیش رونده است. تحکیم رسوبات تحت تأثیر وزن خود … آنها می توانند ناشی از اضافه بار موضعی، سطح خاک در نتیجه ساخت یک کار یا استخراج مواد زیر سطحی، مانند بهره برداری از نفت، آب باشند [ 1 ] .
بر اساس این تعریف، حرکات زمین با دو مؤلفه اساسی مرتبط است: زمین لغزش و خاک در حال چروکیدگی. زمین لغزش و خاک در حال انقباض و متورم شدن بزرگترین گروه خطرات را تشکیل می دهند.
ارزیابی خطر حرکت زمین، از جمله زمین لغزش و خاک در حال انقباض و متورم شدن، موضوع چندین مقاله تحقیقاتی [ 2 ] – [ 11 ] بوده است.
چندین رابطه و مدل های ریاضی خطر ایجاد شده است، بر اساس یک استاندارد واحد یا چندین استاندارد، و متخصص در پدیده های خاص، مانند زمین لغزش یا کوچک شدن و متورم شدن خاک است.
آیا راهی برای جمع آوری و ارزیابی این دو پدیده با هم وجود دارد؟
این را می توان با ادغام تجزیه و تحلیل تصمیم چند معیاره (MCDA) و قابلیت های تحلیلی GIS [ 12 ] به دست آورد.
چندین روش برای حل و ارزیابی مسائل انتخاب در MCDM وجود دارد [ 13 ] [ 14 ]. تکنیک اصلی که نزدیک به این مقاله است، مدل جمع وزنی (WSM) است. WSM را می توان به عنوان یک ترکیب خطی وزن دار یا روش امتیازدهی توصیف کرد [ 15 ] [ 16 ]. WSM نمایشی از ترجیحات تصمیم گیرنده در یک تابع افزودنی خطی است [ 17 ]. Triantaphyllou و Mann [ 18 ] اثربخشی روش های تصمیم گیری را بررسی کردند و WSM را به عنوان یکی از ساده ترین روش ها برای حل مسائل MCDM معرفی کردند [ 19 ].
به طور کلی، فرض کنید که یک مسئله MCDA معین بر روی m جایگزین و n معیار تصمیم تعریف شده است، علاوه بر این، فرض کنیم که همه معیارها معیارهای سود هستند، یعنی هر چه مقادیر بالاتر باشد، بهتر است. در مرحله بعد فرض کنید که w j نشان دهنده وزن نسبی اهمیت معیار C j است و a ij مقدار عملکرد جایگزین A i است زمانی که بر اساس معیار C j ارزیابی می شود. سپس، کل (یعنی وقتی همه معیارها به طور همزمان در نظر گرفته شوند) اهمیت جایگزین A i ، به صورت آدبلیواسم– s c o r eمنAiWSM−score، به صورت زیر تعریف می شود:
آدبلیواسم– s c o r eمن=∑nj = 1آمن جwj, برای n = 1 , 2 , 3 , ⋯ , mAiWSM−score=∑j=1naijwj,forn=1,2,3,⋯,m(1)
برای حالت بیشینه سازی، بهترین جایگزین گزینه ای است که حداکثر مقدار عملکرد کل را به دست آورد [ 20 ].
از (1)، ماهیت شهودی و سهولت پرداختن به مسائل تصمیم گیری چند معیاره قابل مشاهده است و این سادگی، WSM را به یکی از مهم ترین روش ها برای بررسی راه حل های ممکن برای یک مسئله تصمیم گیری چند معیاره و ارائه راه حل در مقایسه تبدیل می کند. به روش های دیگر [ 13 ] .
هدف این کار استفاده از مدلی از خطر جابجایی زمین بر اساس رویکرد تحلیل چند معیاره و ارزیابی توسط GIS است. این مدل با استفاده از روش (WSM) در چارچوب GIS مورد آزمایش قرار گرفته و نتایج به دست آمده مورد بحث قرار گرفته است. نتایج اطلاعات مفیدی را در مورد کاهش خطرات زمین لغزش و کاهش خطرات ناشی از تورم-کوچک شدن ارائه می دهد و ممکن است به عنوان دستورالعملی برای برنامه ریزی کاربری زمین در منطقه برهوم باشد.
این مطالعه بخشی از کمک به مدیریت منطقه برهوم است که در برابر خطرات طبیعی آسیب پذیر است و همچنین در زمینه تولید اطلاعات و اسناد آگاهی. بنابراین، روش های پیشگیرانه بهترین راه برای غلبه بر این مشکلات یا کاهش آنها هستند. یکی از ابزارهایی که امکان اتخاذ روشهای پیشگیرانه مؤثر را فراهم میکند، درک پدیده و پیشبینی مکانهایی است که در آن رخ میدهد.
انتخاب محل تحقیق بر روی منطقه برهوم، واقع در جنوب شرقی الجزایر و بخشی از حوضه الهدنه در شمال شرقی الجزایر متمرکز است ( شکل 1 ). این انتخاب به دلیل تنوع خواص طبیعی و چگالی تجهیزات آن است.
2. روش شناسی
پس از شناسایی پارامترهای دخیل در حرکت زمین، یک مدل ریاضی به خوبی ایجاد شده است.
برای به دست آوردن مدل خطر جابجایی زمین در منطقه مورد مطالعه از روش مدل جمع وزنی (WSM) استفاده شد ( شکل 2 ).
در تجزیه و تحلیل، مقادیر وزن به لایه ها و مقادیر رتبه به طبقات هر لایه اختصاص داده شد. برای هر طبقه از لایه، رتبه بندی بر اساس اهمیت آنها داده شد. پس از تخصیص رتبه بندی به طبقات هر لایه، وزن ها با توجه به اهمیت آنها به لایه ها اختصاص یافت. تعامل بین لایه ها در نظر گرفته نشد زیرا لایه ها مستقل از یکدیگر فرض می شدند [ 15 ] [ 16 ].
شکل 1 . موقعیت منطقه برهوم (منبع: Google Earth).
شکل 2 . نمودار سازمانی یک سیستم اطلاعات جغرافیایی برای نقشه برداری خطر حرکت زمین.
مقادیر وزن و رتبه لایه ها و طبقات هر لایه به منظور به دست آوردن یک واحد بدون بعد مشترک استاندارد شد. پس از آن، خروجی با ضرب مقدار وزن اختصاص داده شده به هر لایه در مقدار رتبه داده شده به طبقات آن لایه و در نهایت با جمع کردن محصولات ایجاد شد [ 22 ].
این مراحل را برای به دست آوردن نقشه زمین لغزش و همچنین نقشه خاک در حال انقباض و متورم شدن اعمال می کنیم. در نتیجه نقشه خطر حرکت زمین منطقه برهوم در زون های مختلف گروه بندی شد.
کاربرد مدل پیشنهادی، از رویکرد SIG استفاده میکند. آن را به صورت زیر بررسی می کنیم:
2.1. فرضیه ها
برای دستیابی به این مدل، باید فرضیاتی مطرح شود. در این زمینه دو فرضیه مطرح می شود:
· حرکات زمین پدیده ای است که با پدیده های زمین لغزش و چروکیدگی-تورم شدن خاک قابل ارزیابی است.
・ خطر مجموع وزنی این عوامل است.
2.2. مدل سازی خطر حرکت زمین
خطر حرکت زمین نیز مجموع وزنی (WSM) عوامل تعیین کننده پدیده است. خطر جابجایی زمین با رابطه زیر ارائه می شود:
اچL M=∑n1wمن⋅سیمنHLM=∑1nwi⋅Ci(2)
جایی که:
w i : وزن تخصیص یافته به هر تعیین کننده، C i : معیار.
طبق تعریف خطر جابجایی زمین، عبارت است از مجموع زمین لغزش و چروکیدگی – متورم شدن خاک. اثر ترکیبی این دو پارامتر با معادله زیر بیان می شود:
اچL M= 0.5اچال اس+ 0.5اچاساسHLM=0.5HLS+0.5HSS(3)
جایی که:
H LM : خطر حرکت زمین، H LS : خطر زمین لغزش، H SS : خطر انقباض و متورم شدن خاک.
2.2.1. مدل سازی خطر زمین لغزش
بر اساس اصل جمع وزنی و در رابطه با شیب، آب و خاک، خطر زمین لغزش با فرمول زیر به دست میآید:
اچال اس= 0.4 P+ 0.4 Tاس+ 0.1 Pنپ+ 0.1 EپHLS=0.4P+0.4TS+0.1PNP+0.1EP(4)
جایی که:
P: زاویه شیب، TS: انواع خاک، PNP: عمق سطح آب زیرزمینی، EP: باران.
2.2.2. مدلسازی خطر انقباض – تورم خاک
بر همین اساس و در رابطه با انواع خاک و آب، خطر انقباض و متورم شدن خاک با فرمول زیر ارائه می شود:
اچاساس= 0.6 Tاس+ 0.2 Pنپ+ 0.2 EپHSS=0.6TS+0.2PNP+0.2EP(5)
جایی که:
TS: انواع خاک، PNP: عمق سطح آب زیرزمینی، EP: باران.
2.2.3. معیارهای مورد استفاده برای انقباض و تورم خاک و خطرات زمین لغزش
زاویه شیب (P): بیشتر زمین لغزش های بزرگ مربوط به زوایای شیب است که مقدار آنها به طور کلی بیش از 45٪ است. مقادیر زاویه شیب به چهار دسته تقسیم شده است ( جدول 1 ).
انواع خاک (TS): بازگشت به طبقهبندی LPC، خاکهای ریز [خاک پلاستیکی کمی (Ap)، خاک رس بسیار پلاستیکی (At)، سیلت پلاستیکی کمی (Lp) و سیلت بسیار پلاستیکی (Lt)]، حساسترین خاکها به تورم هستند. – پدیده انقباض خاک خاک های دانه ای [Gravel (G) و Sand (S)] و سنگ ها به دور از تأثیر این پدیده هستند. انواع خاک به چهار دسته تقسیم شدند ( جدول 2 ).
وجود آب های زیرزمینی و محتوای آب در تغییرات خاک به دلیل تغییرات اقلیمی فصلی، پدیده های تورم- انقباض و رانش زمین را پیچیده می کند.
عمق سطح آب زیرزمینی (PNP): عمق سطح آب زیرزمینی به چهار دسته تقسیم شد ( جدول 3 ).
باران (EP): باران به چهار دسته تقسیم شد ( جدول 4 ).
2.2.4. تقاطع بین پارامترها و طبقه بندی
برای ایجاد تلاقی بین پارامترها، بهتر است آنها را با توجه به درجه تأثیر آنها به کلاس ها تقسیم کنید. نتایج در جدول 5 و جدول 6 فهرست شده است.
ترکیب نتایج وزنی معیارهای مختلف، یک به یک و با در نظر گرفتن عوامل مشترک در هر زمان، سطح هر پدیده را به دست می دهد که در قالب جداول ارائه کرده ایم.
ماتریس های تقاطع زمین لغزش به صورت زیر آورده شده است ( جدول 7 ).
ماتریس های تقاطع برای خاک در حال انقباض و متورم شدن به صورت زیر آورده شده است ( جدول 8 ).
2.2.5. طبقه بندی برای خطرات زمین لغزش و انقباض و تورم خاک
طبقه بندی خطر زمین لغزش در جدول 9 نشان داده شده است.
طبقه بندی خاک در حال انقباض و متورم شدن در جدول 10 نشان داده شده است.
2.2.6. مدل سازی خطر حرکت زمین
برای ترکیبی از لغزش و انقباض – تورم، از فرمول (3) برای مدلسازی نقشه خطر حرکت زمین استفاده میشود. رتبه وزن ها و مقادیر در جدول 11 به شرح زیر است .
ماتریس های تقاطع به صورت جداول زیر آورده شده اند ( جدول 12 و جدول 13 را ببینید).
3. فضاسازی با رویکرد GIS: کاربرد در منطقه برهوم
رویکردی که ما برای تهیه نقشه توزیع فضایی برای ارزیابی خطر حرکت زمین اتخاذ کردهایم بر اساس سیستمهای اطلاعات جغرافیایی است.
مجموعه دادههای مورد استفاده در این مطالعه را میتوان در چهار منبع داده اصلی بهعنوان نقشه سنگشناسی (نقشه سنگها)، بررسی میدانی و دادههای حفاری ژئوتکنیکی، نقشه توپوگرافی و نقشه بارش گروهبندی کرد. این اطلاعات به دست آمده از نقشه توپوگرافی، نقشه بارش، نقشه سنگ شناسی و یا از بررسی میدانی و حفاری ژئوتکنیکی، در برنامه GIS (Arc Gis 9.3) برای ترسیم نقشه های اولیه پیاده سازی شده است.
پس از ساخت بانک اطلاعاتی مرجع جغرافیایی، سازماندهی و ساختاربندی آن با استفاده از برنامه ها، نقشه حرکات زمین خطر را نشان می دهد که بهترین مناطق برای ساخت و ساز را نشان می دهد.
3.1. منابع داده
3.1.1. نقشه سنگ شناسی
نقشه سنگ شناسی برگرفته از نقشه زمین شناسی منطقه مورد مطالعه [ 23 ] . پس از شمارهگذاری نقشه زمینشناسی، پایگاه دادههای صخرهای دارای ساختمان است. برای منطقه مورد مطالعه، نه واحد سنگ شناسی شناسایی شد ( شکل 3 ). این واحدها شامل دولومیت، سنگ آهک، ماسه سنگ، کنگلومرا، مارن، رس، آبرفت درشت شنی-رسی و آبرفت ریز می باشد. در منطقه مورد مطالعه، سازندهای دولومیت، سنگ آهک، ماسه سنگ و کنگلومرا به عنوان سنگ طبقه بندی شدند.
3.1.2. داده های بررسی میدانی و گمانه ژئوتکنیکی
در واقع، از سه تکنیک متفاوت اما مکمل برای تکمیل پایگاه داده استفاده شده است:
・ داده های مربوطه 49 بررسی ژئوتکنیکی مطالعه بررسی و ارزیابی شد.
・ بررسی میدانی خوب بر اساس مشاهدات مستقیم،
・ و همچنین ضبط و نمونه برداری میدانی (GPS، تصاویر و غیره).
انواع خاک (TS) و لایه های سطح آب زیرزمینی با استفاده از نقشه سنگ شناسی و گمانه ژئوتکنیکی تهیه شده است ( شکل 4 و شکل 5 ).
3.1.3. نقشه توپوگرافی
برای بدست آوردن مدل رقومی ارتفاع (DEM) منطقه مورد مطالعه از نقشه توپوگرافی در مقیاس 1/50000 استفاده شده است. دومی با دیجیتالی کردن خطوط کانتور نقشه توپوگرافی به دست آمد. این نرم افزار Arc Gis 9.3 است که برای درونیابی استفاده شده است. نقشه شیب ها بر اساس درجه، با استفاده از DEM منطقه مورد مطالعه تهیه شده است ( شکل 6 ).
3.1.4. نقشه بارش
برای بدست آوردن لایه آب باران منطقه برهوم از نقشه بارندگی در مقیاس 1/50000 استفاده شده است ( شکل 7 ).
3.2. ارزیابی داده ها
این داده ها برای به دست آوردن نقشه های خطرات زمین لغزش و انقباض- تورم خاک تجمیع و ارزیابی شدند ( شکل 8 و شکل 9 ). برای اجرای این مراحل، ما از شبکهسازی ماشینابزار ماشینابزار (Calculator Raster) نوار ابزار تحلیل فضایی (Suitability Special) در Arc Gis 9.3 [ 24 ] استفاده میکنیم.
3.3. شناسایی و طبقه بندی خطر حرکت زمین
ترکیب نتایج زمین لغزش و انقباض – متورم شدن خاک (فرمول 3)، درجه خطر حرکت (H LM ) را به دست می دهد که در شکل 10 ارائه شده است.
4. اعتبار سنجی نقشه خطر حرکت زمین
اعتبار سنجی نتایج در دو سطح انجام شد. از یک طرف، نقشه های موضوعی به دست آمده با داده های میدانی قبل از تلاقی آنها مقایسه شد. تطابق کامل بین این مشاهدات میدانی ( شکل 11 ) و مکان های خطر متوسط و بالا که توسط نقشه ساخته شده از GIS نشان داده شده است، امکان اعتبارسنجی نقشه نهایی را تا حدودی فراهم می کند.
شکل 3 . نقشه سنگ شناسی منطقه برهوم.
شکل 4 . لایه انواع خاک
شکل 5 . لایه سطح آب زیرزمینی
شکل 6 . لایه شیب.
شکل 7 . لایه آب باران.
شکل 8 . نقشه خطر زمین لغزش
شکل 9 . نقشه خطر خاک در حال انقباض و متورم شدن
شکل 10 . نقشه خطر حرکات زمینی منطقه برهوم.
شکل 11 . عکس های صحرایی که مشکلات انقباض و متورم شدن خاک را در منطقه مورد مطالعه نشان می دهد.
از سوی دیگر، در طول ماموریت های میدانی، نظرسنجی انجام شده (چه از طریق دفاتر فنی شهرداری ها، چه توسط آژانس بازرسی فنی ساختمان ها (CTC-EST، Agence de M’SILA) و چه به طور مستقیم در میدان) هیچ گونه اسلاید نمایشی را ثبت نکرد. در منطقه مورد مطالعه، و بنابراین هیچ گونه خسارت به ساختمان های ناشی از رانش زمین ثبت نشده است. همچنین وجود برخی آسیبهای محدود به ساختمانها را ثبت کردیم (از طریق گزارشهای آژانس بازرسی فنی ساختوسازها)، و همان مشاهداتی است که [ 25 ] در تحقیقات خود در هر یک از شهرداریهای برهوم و اولد آدی گوبالا به دلیل کاهش مشکلات خاک متورم [ 26 ].
5. بحث
گسترش شهری و افزایش جمعیت شهرنشینی را به سمت استفاده از فضاهای «در معرض خطر» سوق داد. ادغام و برخورد با حرکات اراضی مشخصه مهمی از تعادل ایجاد شده توسط طبیعت است [ 11 ].
برای حل این مشکلات، نیاز به ارائه نقشه خطر جابجایی زمین با استفاده از سیستم های اطلاعات جغرافیایی (GIS) با اتخاذ مدل (WSM) بود. این ابزارها ما را قادر می سازند تا درک فضا و پردازش داده ها را بیان کنیم و در نتیجه کارتوگرافی به روشی بهینه انجام می شود. این نقشه خطر جابجایی زمین در تمامی سطوح مطالعه به عنوان اطلاعات، کاغذ کار، هشدار و به ویژه ابزاری برای کمک به تصمیم گیری با بیان گرایش ها و جهت گیری ها از اهمیت بالایی برخوردار است.
نقشه خطر جابجایی زمین منطقه مورد مطالعه ( شکل 10 ) و یک پایگاه داده با تمام اطلاعات ژئوتکنیکی گردآوری شده (GEODATA) کمک بسیار مهمی به برنامه ریزی و مدیریت زمین منطقه مورد مطالعه و تعریف مناطق خطر طبیعی است. این یک منبع داده های ژئوتکنیکی عالی را فراهم می کند تا زمانی که پروژه های جدید برای منطقه مورد مطالعه برنامه ریزی و طراحی می شوند، مورد استفاده قرار گیرد.
یک پایگاه داده باز، با اطلاعات مربوط به گمانهها، مکان جدول آب، دادههای آزمایش آزمایشگاهی، مکانهای نمونه و عکسهای نمونهها، باید به نقشههای سنگشناسی، انواع خاک و جدولهای آب زیرزمینی اضافه شود. چنین پایگاه داده ای امکان ارزیابی سریع شرایط خطر جابجایی زمین را فراهم می کند.
دادهها از گزارشهای ژئوتکنیکی ارائهشده توسط مؤسسات دولتی و همچنین از کار میدانی و تحقیقات خود ما به جمعآوری دادهها اجازه میدهد تا نقشههای زمین لغزش و کوچک شدن و متورم شدن خاک را به دست آورد ( شکل 8 و شکل 9 ). این داده ها برای محاسبه خطر حرکت زمین ( شکل 10 ) با رابطه (3) استفاده شده است.
مدل پیشنهادی ساده است و به راحتی توسط مهندس ژئوتکنیک یا موارد دیگر استفاده می شود. مدل رزوه بسیار جالب است زیرا بر اساس تشخیص کارشناسان بر روی تعداد کمی از ویژگی های ژئوتکنیکی ساخته می شود. به طور مشابه، این مدل را می توان به راحتی در سایر مناطق مشابه از نظر چشم انداز محیطی و آب و هوا اعمال کرد. با این حال، این مدل ثابت نیست. گاهی اوقات می توان با توجه زیاد به یک عامل تغییر داد و آن را به عامل دیگر کاهش داد.
6. نتیجه گیری
این کار در نقشه کشی خطر جابجایی زمین، با استفاده از تکنیک های سیستم های اطلاعات جغرافیایی (GIS)، با مدل وزنی سام (WSM)، نشان دهنده علاقه ناشی از این روش ها است که به نظر می رسد به اشکال مختلف، راه را برای پیشرفت بسیار حساس باز می کند. پیشگیری از خطرات طبیعی این ابزارها ما را قادر می سازند تا درک فضا و پردازش داده ها را بیان کنیم و در نتیجه کارتوگرافی به روشی بهینه انجام می شود.
این نقشههای خطر جابجایی زمین به عنوان اطلاعات، کاغذ کار، هشدار و به ویژه ابزاری برای کمک به تصمیمگیری با بیان گرایشها و جهتگیریها، در تمامی سطوح مطالعه از اهمیت بالایی برخوردار است.
در این کار نقشه های زمین لغزش و انقباض – تورم خاک برای منطقه برهوم تهیه شده است. مقادیر وزن ها، تجزیه و تحلیل و طبقه بندی توسط مدل Weighted Somme اختصاص داده شده است. در نتیجه و با توجه به کارایی پی، محدوده مورد مطالعه به چهار منطقه مختلف تقسیم شده است: 1) مناطق بدون خطر جابجایی زمین، 2) مناطق با خطر جابجایی زمین متوسط، 3) مناطق با خطر جابجایی زمین بالا، و 4) مناطق با خطر جابجایی زمین بسیار بالا. بر اساس این نقشه، تراکم کم هر نوع حرکت منطقه مورد مطالعه شناسایی شد.
نقشه خطر حرکات زمین وسیله مبارزه نیست، بلکه به ایجاد طرح شهرنشینی کمک می کند و کنترل کافی، علاوه بر مدیریت مشکل شهرنشینی در مناطق پرخطر، امکان پذیر و قابل کنترل می شود.
روش مورد استفاده در این مطالعه برای منطقه برهوم می تواند در مکان های دیگر، روش های دیگر برای انتخاب مکان و بازسازی استانداردهای لازم به طور مناسب اعمال شود.
منابع
بدون دیدگاه