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

نوع مقاله : مقاله پژوهشی

نویسندگان

1 دانشجوی دکتری مهندسی منابع آب، پردیس ارس، دانشگاه تهران.

2 استاد گروه مهندسی آبیاری و آبادانی، دانشکده مهندسی و فناوری کشاورزی، دانشگاه تهران کرج ایران

10.22034/gmpj.2020.122222

چکیده

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

کلیدواژه‌ها


عنوان مقاله [English]

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

نویسندگان [English]

  • Mohammadreza Asgari 1
  • |KUMARS EBRAHIMI 2
1 Water Engineering Department, ' Aras International Campus, University of Tehran
2 Department of Irrigation & Reclamation Engineering Faculty of Agricultural Engineering and Technology College of Agriculture & Natural Resources, University of Tehran, Karaj, Iran
چکیده [English]

-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 long-term 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 spatial-statistical 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 MODFLOW-2005 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 non-stable 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 0-1 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 land-use 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 enhancement-reduction 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 land-uses and a basic land-use set were extracted. These results revealed that the urban land-use accounted for 26% of subsidence, even though it occupied merely 4% of the aquifer’s area. The agricultural land-use (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 land-use 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.

کلیدواژه‌ها [English]

  • Modeling
  • fuzzy model
  • Abhar Aquifer
  • subsidence
شریفی کیا، م.، افضلی، ع و شایان، س.، 1394. استخراج و ارزیابی اثرات پدیده‌های ژئومورفولوژیک ناشی از فرونشست در دشت دامغان، فصلنامه پژوهش‌های ژئومورفولوژی کمّی، سال چهارم، شماره 2، پاییز 1394، صص.74-60.
صفاری، ا.، جعفری، ف و توکلی صبور، س، م.، 1395. پایش فرونشست زمین و ارتباط آن با برداشت آب‌های زیرزمینی مطالعه موردی: دشت کرج-شهریار، فصلنامه پژوهش‌های ژئومورفولوژی کمّی، سال پنجم، شماره 2، پاییز 1395، صص.93-82.
قدرتی، م و ثعبانی، ع.، 1391. مدل‌های ریاضی آبهای زیرزمینی ،آموزش کاربردی مدل GMS، چاپ اول، انتشارات سیمای دانش، 210ص.
Ahmed, W, A., Kalkan, E. Guzy, A. Alacali, M. Malinowska, A., 2020, Modeling of land subsidence caused by groundwater withdrawal in Konya Closed Basin, Turkey, Tenth International Symposium on Land Subsidence (TISOLS), 382, 397–401.
Al-Sittawy, M. Gad, S. Fouad, R. Nofal, E.; 2019, Assessment of soil subsidence due to long-term dewatering, Esna city, Egypt. Journal Water Science Pages 40-53 | Received 05 Mar 2019, Accepted 30 May 2019, Published online: 25 Jul 2019.
Bear, J., 1979. Hydraulics of Groundwater, McGraw-Hill, New York, 569p.
Chenini, I. Mammou, A, B., 2010, Groundwater recharge study in arid region: Anapproach using GIS techniques and numerical modeling, Computers & Geosciences, Pages 801–817.
Cho, J. Barone, V, A. Mostaghimi, S., 2009, Simulation of land use impacts on groundwater levels and streamflow in a Virginia watershed, agricultural water management, Volume 96, pages 1–11.
Garfias, J. Richard, M. Angus, C., 2020, Assessment of groundwater exploitation and land subsidence development in the Toluca aquifer system, Mexico, Repositorio Instrucional, 290-302.
Gaura, S. Chahar, B, R. Graillota, D., 2011, Combined use of groundwater modeling and potential zone analysis for management of groundwater, International Journal of Applied Earth Observation and Geoinformation, Volume 13, Pages 127–139.
Hu, Y. Moiwo, J, P. Yang, Y. Han, S. Yang, Y., 2010, Agricultural water-saving and sustainable groundwater management in Shijiazhuang Irrigation District, North China Plain, Journal of Hydrology, Volume 393, Pages 219–232.
Kooi, H. Erkens, G., 2020, Creep consolidation in land subsidence modelling; integrating geotechnical and hydrological approaches in a new MODFLOW package (SUB-CR), Proceedings of the International Association of Hydrological Sciences, Proc. IAHS, 382, 499–503.
Lachaal, F. Mlayah, A. Bédir, M. Tarhouni, J. Leduc, C., 2012, Implementation of a 3-D groundwater flow model in a semi-arid region using MODFLOW and GIS tools: The Zéramdine–Béni Hassen Miocene aquifer system (east-central Tunisia), Computers & Geosciences, Volume 48, Pages 187–198.
Le Page, M, L. Berjamy, B. Fakir, Y. Bourgin, F. Jarlan, L. Abourida, A. Benrhanem, M. Jacob, G. Huber, M. Sghrer, F. Simonneaux, V. Chehbouni, G., 2012, An Integrated DSS for Groundwater Management Based on Remote Sensing. The Case of a Semi-Arid Aquifer in Morocco, Water Resour Manage, Volume 26, Pages 3209–3230.
Malczewski, J., 1999. GIS and Multicriteria Decision Analysis. United States, 408 p.
Ranjbar, A. Ehteshami; M., 2019, Development of an Uncertainty Based Model to Predict Land Subsidence Caused by Groundwater Extraction (Case Study: Tehran Basin). Geotechnical and Geological Engineering, Volume 37, Issue 4, pp 3205–3219
Rashvand M, Li, J. Liu, Y., 2019, Coupled Stress-Dependent Groundwater Flow-Deformation Model to Predict Land Subsidence in Basins with Highly Compressible Deposits. Journal of Hydrology 2019, 6(3), 78.
Yang, F, R. Lee, C, H. Kung, W, J. Yeh, H, F., 2009, The impact of tunneling construction on the hydrogeological environment of “Tseng-Wen Reservoir Transbasin Diversion Project” in Taiwan, Engineering Geology, Volume 103, Issues 1–2, Pages 39–58.
Yaoutia, F, E. Mandourb, A, E. Khattacha, D. Kaufmannc, O., 2008, Modelling groundwater flow and advective contaminant transport in the Bou-Areg unconfined aquifer (NE Morocco), Journal of Hydro-environment Research, Volume 2, Issue 3, 12 December 2008, Pages 192–209.
Zhang, H. Hiscock, K, M., 2010, Modelling the impact of forest cover on groundwater resources: A case study of the Sherwood Sandstone aquifer in the East Midlands, UK, Journal of Hydrology, Journal of Hydrology, Volume 392, Issues 3–4, Pages 136–149.