مدل سازي فرونشست دشت ايذه با استفاده از كد رياضي MODFLOW

Izeh Plain Subsidence Modeling Using MODFLOW Mathematical Code

فرونشست ساختار دشت‌ها، مي‌تواند مستقيما بر اثر افت سطح آب زيرزميني و تخريب بافت آبرفتي آبخوان صورت پذيرد. در اين مطالعه در دشت ايذه مدل رياضي MODFLOW، با استخراج نقشه رستري فرونشست آبخوان اشباع در دوره ده ساله از 2008 تا 2019 و دوره پيش بيني سه سناريو اقليمي بين سال‌هاي 2019 تا 2029 انجام شد. مدل فرونشست با توسعه پكيج SUB در ساختار كد عددي تفاضل محدود MODFLOW صورت پذيرفت. مدل كمي جريان آب زيرزميني با استفاده از كد خودكار PEST واسنجي و تحليل حساسيت شد. نتايج بدست آمده حاكي از آن است كه مدل رياضي استفاده شده جهت شبيه‌سازي آبخوان ايذه داراي خطاي نسبي NRMSE مطلوب مي‌باشد كه مويد مدل سازي ايده آل پس از بررسي فرآيند صحت سنجي است. بررسي تغييرات عمودي ساختار زمين نشان داد كه در بلند مدت 20 ساله با فرضيات اقليمي مطرح شده سطح آبخوان مجموعا حداكثر تا 5/1 متر و حداقل 9/0 متر فرونشست رخ مي‌دهد. در صورت تخريب ساختار آبرفتي آبخوان اين مقدار ممكن است كمتر از اين ارقام بوده اما در مقابل برابر با شرايط بدون بازگشت تعذيه سفره آب زيرميني باشد.

The subsidence of the plains structure can be directly due to the drop in groundwater level and the destruction of aquifer alluvial tissue. In this study, in the Izeh plain, the MODFLOW mathematical model was performed by extracting a raster map of saturated aquifer subsidence in the ten-year period from 2008 to 2019 and the forecast period for three climatic scenarios between 2019 and 2029. The subsidence model was developed by developing the SUB package in the limited numerical code structure structure of MODFLOW. The quantitative groundwater flow model was calibrated using PEST automation code and sensitivity analysis. The results show that the mathematical model used to simulate Izeh aquifer has a relative relative error NRMSE, which confirms the ideal modeling after reviewing the validation process. An examination of the vertical changes in the structure of the land showed that in the long term of 20 years, with the climatic assumptions raised, the aquifer level will occur a total of up to 1.5 meters and a subsidence of at least 0.9 meters. If the aquifer alluvial structure is destroyed, this amount may be less than these figures, but in return, it is equal to the conditions without returning the groundwater aquifer feed.