شبیه سازی عددی زمین لغزش در اثر بارش های با شدت های متفاوت

نوع مقاله : پژوهشی اصیل (کامل)

نویسنده
گروه مهندسی آب، دانشگاه تبریز
چکیده
هر ساله وقوع باران­‌های شدید باعث ایجاد شکست در تعداد زیادی از شیب‌های خاکی می­شود .در طول این باران‌ها، سطح آب زیرزمینی بالا آمده و باعث افزایش فشار آب منفذی و کاهش پایداری شیب می­گردد. هدف از این تحقیق، بررسی تغییرات فشار آب منفذی و ضریب اطمینان پایداری شیب است. جریان آب درون خاک، به صورت دو بعدی، اشباع- غیراشباع و غیردائمی مدلسازی عددی می­شود. از روش اجزای محدود برای محاسبه فشار آب منفذی و از روش تعادل حدی برای تعیین ضریب اطمینان پایداری شیب استفاده می­شود. برای شبیه­سازی جریان غیردائمی، زمان برابر 4 روز، با گام 1/0 روز استفاده می­شود. مدت بارش 24 ساعت فرض شد. شدت بارش­های مختلف به صورت شار ورودی روی سطح خاک، معرف یکی از شرایط مرزی استفاده شد. منحنی رطوبتی خاک که نشان دهنده رابطه مکش محتوی رطوبتی و مکش- هدایت هیدرولیکی است، استفاده شد. نتایج نشان داد که شدت بارش­های مختلف، اثرات متفاوتی روی پروفیل رطوبتی خاک دارند. با افزایش شدت بارش، ضریب اطمینان پایداری شیب کاهش یافت. این کاهش برای 24 ساعت اول بارندگی شیب بیشتری داشت و در سه روز بعدی بعد از قطع بارش، شیب آن کمتر شد. به عنوان مثال، برای شدت بارش 04/2 میلی­متر بر ساعت و تدوام یک روز، نتایج نشان داد که در شروع بارش روی سطح خاک، ضریب اطمینان پایداری شیب برابر 853/1 بود که در پایان روز چهارم، به مقدار 743/1 کاهش یافت. به بیان دیگر، ضریب اطمینان پایداری شیب حدود 9/5 درصد کاهش یافت.

کلیدواژه‌ها

موضوعات


عنوان مقاله English

Numerical simulation of landslides caused by rainfall of different intensities

نویسنده English

Farzin Salmasi
Department of Water Engineering, Faculty of Agriculture, University of Tabriz
چکیده English

Foothills, plains, alluvial areas, and sloping areas whose soil is geologically susceptible to landslides, can become unstable and dangerous. In Iran, because of their soil type, Mazandaran, Gilan, Lorestan, Golestan and Khuzestan provinces are more vulnerable to landslides than other provinces. But taking into account factors such as earthquakes, extreme weather, and human interference, other parts of the country can also be prone to landslides. In case of rainfall and absorption of water by clay layers, the possibility of landslides increases. If the slope of the land is suitable, the slope will move on the clay layer and the thrust will occur downwards. In many mountains and steep areas, the conditions for slope failure may be available in terms of the slope angle, the type of soil and the presence of clay layers. But in the absence of sufficient soil moisture, this phenomenon is not observed. Landslides occur whenever rainfall happens and water penetrates the clay layer. Cutting down forest trees, destroying vegetation and pastures, along with possible rainfall and soil moisture, cam cause landslides on steep slopes. In fact, land use changes contribute to landslides. By examining the statistics over the last three decades in the northern provinces of Iran, it can be seen that landslides were far less likely when there was proper vegetation. The occurrence of heavy rains can cause failure in a large number of soil slopes. During these heavy rainfalls, the underground water level rises and increases the pore water pressure and reduces the stability of the soil. The purpose of this study is to investigate changes in pore water pressure and the factor of safety for slope stability. The flow of water inside the soil is numerically modeled as a two-dimensional, saturated-unsaturated, unsteady flow. The finite element method (FEM) is used to calculate the pore water pressure and the limit equilibrium method is used to determine the factor of safety for slope stability. To simulate the unsteady flow, time duration of 4 days is used with a time step of 0.1 days. The rainfall duration was assumed to be 24 hours. The intensities of different rainfalls are used as the input flux on the soil surface. The soil moisture curve, which shows the relationship between suction-moisture content and suction-hydraulic conductivity, is used. The results show that different rainfall intensities have different effects on soil moisture profile. As the intensity of precipitation increase, the factor of safety of slope stability decreases. This decrease is steeper for the first 24 hours of rainfall and over the three days after the rain stopped, its slope decreased. For example, for a rainfall intensity of 2.04 mm/hour and the duration of one day, resulted in a factor of safety of slope stability equal to 1.853, and at the end of the fourth day, it was 1.743. In other words, the factor of safety decreased by 5.9%.

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

Pore water pressure
soil slope stability
factor of safety
rainfall intensity
unsteady flow
unsaturated flow
[1] Rahardjo H, Ong TH, Rezaur RB, Leong EC, 2007. Factors controlling instability of homogeneous soil slopes under rainfall. J Geotech Geoenviron Eng, 133:1532–1543 20.
[2] Dhanai P, Singh VP, Soni P, 2022. Rainfall Triggered Slope Instability Analysis with Changing Climate, Indian Geotech J, 52(2):477–492
[3] Salmasi F, Hosseinzadeh Dalir A, Norouzi Sarkarabad R, 2019. Investigation of the performance of horizontal drains in increasing slope stability in intense rainfall conditions by numerical simulation, Amirkabir Journal of Civil Engineering, Amirkabir J. Civil. Eng., 51(3), 153- 156, (In Persian with English abstract)
[4] Salmasi F, Pradhan B, Nourani B, 2019. Prediction of the sliding type and critical factor of safety in homogeneous finite slopes. Applied Water Science, 9:158, 1-11
[5] Moharrami A, Hassanzadeh Y, Salmasi F, Moradi G, Moharrami G (2014). Performance of the horizontal drains in upstream shell of earth dams on the upstream slope stability during rapid drawdown conditions. Arabian Journal of Geosciences. 7(5): 1957–1964
[6] Lau KC, Kenney TC, 1984. Horizontal drains to stabilize clay slopes, Canadian Geotechnical 21(2), 241-249.
[7] Ahmadi H, Verdinejad VR, Hemati M, (2013). Determination of stable side slops in irrigation canals, Journal of Water and Soil, Publisher: Ferdowsi Mashad University, 27 (3), 603-612, (In Persian with English abstract)
[8] Subramaniya K., Madhav R.M., and Govinda C.M. 1973. Studies on Seepage from canals with partial lining. Journal of Hydraulic Division, 99(12), 2333-2351
[9] Taheri Aghdam, A., Salmasi, F., Abraham, J., Arvanaghi, H. 2021. Effect of drain pipes on uplift force and exit hydraulic gradient and the design of gravity dams using the finite element method. Geotechnical and Geological Engineering 5 (3), 12-29
[10] Sartipi N, Salmasi F, Abraham J, Hosseinzadeh Dalir A, 2021. Investigation of the effect of depth and distance between cutoff walls on uplift force for gravity dams, International Journal of Environmental Science and Technology, 18, 1361–1378
[11] Van Genuchten, M.T. 1980. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal 44(5): 892-898.
[12] Nourani B., Salmasi F., Abraham, J. 2025. Investigating Seepage and Stability of Unsaturated Embankment Dams Under Transient Conditions. Transportation Infrastructure Geotechnology. 12, 17 (2025).