بررسی پدیده فرونشست دشت ابهر با استفاده از مدل ریاضی modflow و بر مبنای توسعه بسته عددی sub

در این مقاله مدل‌سازی فرونشست دشت ابهر با مدل شبیه‌ساز modflow و استخراج نقشه رستری فرونشست آبخوان اشباع در یک دوره ده ‌ساله انجام شده است. مدل فرونشست با توسعه زیرمدل sub در ساختار کد عددی تفاضل محدود modflow انجام شد. دوره شبیه‌سازی به دو مرحله جهت واسنجی مدل کمی جریان در شرایط غیرماندگار از روش تحلیل حساسیت تلفیقی و با بکارگیری مدل pest و جهت صحت‌سنجی تقسیم شد. نتایج نشان داد شبیه‌سازی آبخوان دارای 8% خطای نسبی می‌باشد که موید مدل‌سازی ایده‌آل است. بررسی تغییرات عمودی ساختار لایه های زمین نشان داد در دوره ده ساله آبخوان 34 سانتی‌متر فرونشست داشته است. در بررسی مدل فازی تکمیلی با استفاده از همپوشانی گاما بین لایه‌های موثر با رخداد بیشینه همبستگی رگرسیون خطی لایه نقاط فرونشست مشخص شد که تحلیل اثر کاربری اراضی با تفسیر تغییرات فرونشست مدل فازی مطابقت دارد. با وجود آنکه کاربری شهری تنها 4 درصد از سطح آبخوان را شامل می‌شود ولی 26 درصد از واقعه فرونشست و در مقابل کاربری کشاورزی که 42 درصد از سطح می باشد 56 درصد از فرونشست و زمین‌های بایر با 54 درصد از سطح آبخوان اشباع تنها 19 درصد از فرونشست دشت را به خود اختصاص داده‌اند. با به‌کارگیری همپوشانی وزنی لایه‌های ایجاد شده با استفاده از عملگر گاما در مدل فازی جهت پهنه بندی بعنوان روشی نوین در محاسبات نرخ فرونشست زمین به تفکیک نوع کاربری اراضی، می‌توان مناطق مستعد فرونشست زمین را شناسایی کرد تا با مدیریت صحیح، از وقوع فرونشست و تاثیرات مخرب آن جلوگیری نمود.

Studying land subsidence phenomenon of Abhar plain using MODFLOW mathematical model based on developing SUB numeric package

Introduction The reduced groundwater levels of plains have led to increased water extraction cost, increased energy consumption, reduced water quality, and the appearance of subsidence. Structural subsidence in plains could be directly resulted from reduced level of groundwater as well as destruction of alluvial texture in aquifers. Although it results from the compacted underlying layers of soil, subsidence represents unpleasant outcomes in the future. The compacted underlying layers and reduced water table in groundwater basins reduce the water storage space. In other words, reduced groundwater storage is a rational consequence of subsidence (Sharifikia et al., 2014). Scientific investigation and various instruments, including mathematical models, are required to overcome the problems of groundwater resources (Gaura et al., 2011; Hu et al., 2010). Mathematical models allow for investigating changes in the current and future situations of groundwater tables by incorporating different influential factors (Yaoutia et al., 2008; Zhang, 2010). Therefore, using mathematical models, the situation of an aquifer can be simulated by collecting information on the inputs and outputs of the groundwater system at a small cost in a short time. Such simulation models can provide the mutual effects of surface water and groundwater in short and longterm periods. MODFLOW is among the most important mathematical models (Lachaal et al., 2012). In this current paper the approaches to minimize subsidence in Abhar plain is evaluated that based on an empirically derived relationship between cumulative subsidence rates and groundwater levels is used with common groundwater model software MODFLOW. Also, the vertical changes in the ground structure were nonlinearly modeled based on the cellular changes of the finite difference method. To determine the essential factor of the changes, the fuzzy model approach and spatialstatistical analysis were adopted. The largest effect on the appearance of subsidence was studied by the regression relationships of the corresponding points in the spreadsheet setting. Material and methods Abhar Plain is located in the northwest of the Namak Lake Basin. It occupies an area of 1926.5 km2, 1040.92 km2 of which is plain, while the remaining area is composed of mountains. The maximum and minimum elevations of Abhar Plain are 2166 and 749 m, respectively. The saturated area of the plain is approximately 657.9 km2 – the corrected area is 657.4 km2. The groundwater extraction area of the plain includes 1359 water wells with an annual discharge of 233.36 million m3, 169 fountains with an annual discharge of 4.88 million m3, and 101 aqueducts with an annual discharge of 1.3 million m3. Considering that parameters involved in the calculation of the final subsidence layer had reduction or enhancement effects on each other and given the descriptions on the overlapping functions, the current study employed the gamma operator. Computation was performed by the MODFLOW2005 engine in GMS v.10. The study period was selected to be 119 months, which could be made more accurate and rebuilt based on the maximum available data. The computation engine of PCG2 with 100 outer and inner iterations, a critical convergence variation limit of 0.01 m, and a critical convergence error limit of 0.01 m3 per day was selected. 75% of the interval length was used for calibration in nonstable conditions. After seven executions of the calibration model with a certain number of internal iterations, the optimal final number of surface feeding, horizontal hydraulic conductivity and horizontal hydraulic conductivity anisotropy, specific yield in the form of pilot extraction points, and transferability parameters in the boundaries, and the waterway network in the form of the bunches of lines. The layers were also used to develop the fuzzy subsidence model. According to the described fuzzy theory, each of the basic parameters influencing or representing subsidence was transformed into a standard raster map in the range of 01 by a linear equation. The fuzzified layers included water level reduction, the difference between the initial and final water levels, horizontal hydraulic conductivity anisotropy, horizontal hydraulic conductivity, saturated aquifer depth, surface feeding, extraction flow rate, and specific yield. Direct fuzzification was applied to all the parameters, except for surface feeding, to which inverse fuzzification was applied. Table 2 shows the zoning errors of the selected layers. The parameters of Table 2 were used to develop the fuzzy model. The coefficients of the aquifer were merely extracted from the last calibration round of the groundwater flow model. In the next step, to investigate the spatial variations of subsidence occurrence, the landuse layer was utilized as the statistical analysis basis of the subsidence raster output. Conclusion The results of the current study indicated that subsidence modeling can be carried out by linear regression analysis with a determination correlation coefficient of 70%. Among the influential parameters, merely the spotted aquifer depth layer had a correlation of 30%. The fuzzified variations of the layers with the enhancementreduction gamma combination had the highest correlations with the spotted subsidence layer of the MODFLOW output. However, the correlation was noticeably lower than expected. In the subsidence model, the final raster layer was transferred to the GIS setting using dispersed points. After zoning using the regional analysis command in two manners, the separate indicators of landuses and a basic landuse set were extracted. These results revealed that the urban landuse accounted for 26% of subsidence, even though it occupied merely 4% of the aquifer’s area. The agricultural landuse (including gardens), which accounted for 42% of the aquifer’s area, involved 56% of subsidence events. Finally, occupying 54% of the aquifer’s area, the idle landuse accounted for merely 19% of the vertical changes of the aquifer’s structure. The high subsidence in urban areas was noticeable. Considering the alluvial groundwater yield, the effect of water level reduction can be seen with small spatial distances.