مدلسازی عددی فرونشست زمین ناشی از استخراج آب های زیرزمینی با لحاظ کردن اثرات غیر اشباع و به روش بدون المان گالرکین (efg)

در پژوهش حاضر، به کمک یک مدل عددی کاملاً همبسته‌ی هیدرومکانیکال، پدیده‌ی فرونشست ناشی از استخراج آب‌های زیرزمینی از آبخوان‌ها با درنظرگرفتن آثار غیراشباع به روش بدون‌المان‌گالرکین (efg) شبیه‌سازی شده است. به کمک مدل اخیر می‌توان یک پیش‌بینی مناسب برای تغییرات هد هیدرولیکی و میزان فرونشست زمین ارائه کرد. پس از حل مسائل محک، مدل نهائی برای حل مسئله‌ی‌ فرونشست استفاده شد. نتایج پژوهش حاضر با نتایج روش المان‌محدود (fem) برای یک آبخوان همسان و با درنظرگرفتن شرایط غیراشباع مقایسه و خطای فرونشست زمین 5.51٪ و خطای افت هدهیدرولیکی 13.35٪ محاسبه شد، که نشان از تطابق مناسب بین نتایج دو روش اخیر است. نتایج نشان می‌دهند که روش efg قابلیت مدل‌سازی پدیده‌ی فرونشست را دارد. همچنین آثار غیرخطی استخراج آب‌های زیرزمینی در کاهش هد هیدرولیکی و افزایش فرونشست قابل تعیین است. درنهایت، یک مطالعه‌ی پارامتری به‌منظور بررسی میزان تاثیر ویژگی‌های ژئوتکنیکی آبخوان در میزان فرونشست انجام شد. پارامترهای مدول کشسان و ضریب پواسون بیشترین تاثیر را در میزان فرونشست زمین دارند و لازم است پیش از اجرای پروژه‌های برداشت آب‌های زیرزمینی، پارامترهای مذکور با دقت قابل قبول درنظر گرفته ‌شوند.

numerical modeling of land subsidence induced by groundwater extraction considering unsaturated effects and using element free galerkin (efg) method

this study aims to provide a coupled flow deformation model for simulating land subsidence associated with groundwater extraction in aquifers. for this simulation, we have adopted the element free galerkin (efg) method and considered the unsaturated effects in the aquifers based on the aquifer’s hydrologic and geotechnical characteristics. this model gives us a better understanding of the aquifer’s hydrogeological characteristics, enabling us to forecast changes in the hydraulic head and land subsidence. to ensure the credibility of our model and to verify the code, we modeled unsaturated hydromechanical benchmark problems. then, using the efg method as a numerical tool, we modeled an isotropic aquifer to investigate the effects of groundwater pumping on land subsidence and hydraulic changes in the aquifer. to ascertain the reliability of the modeling, we compared the results obtained from the efg method with those from the finite element method (fem). the comparative analysis of efg and fem models demonstrates discrepancies of 5.51% in land subsidence and 13.35% in hydraulic head reduction, which are satisfying. the land subsidence and hydraulic head profiles demonstrate that the efg method is capable of land subsidence simulation caused by water pumping. furthermore, our findings highlight the nonlinear correlation between groundwater extraction and the subsequent decrease in hydraulic head and land subsidence augmentation. finally, we conducted a parametric study to better understand the effect of various characteristics of aquifers and observe the effect of the aquifer’s parameters, such as hydraulic conductivity, elastic modulus, and poisson’s ratio. we investigated the effect of each parameter on land subsidence increase and hydraulic head decline. the results show that elastic modulus and poisson’s ratio have the most significant effect on land subsidence. although hydraulic conductivity controls the hydraulic decrement and land subsidence increase time, it slightly affects the ultimate hydraulic head and land subsidence at the steady state stage. these results highlight the importance of in situ measurement of elastic modulus and poisson’s ratio paraeters with acceptable accuracy for groundwater extraction projects, as these parameters play a significant role in the feasibility studies.