Integrated D-InSAR and Ground-based Radar for Open Pit Slope Stability Monitoring and Implications for Rock Mass YoungAos Modulus Reduction \\\\\\\\ Faris Ridwan Maulana.
Ridho Kresna Wattimena & Budi Sulistianto Mining Engineering Department.
Faculty of Mining and Petroleum Engineering.
Institut Teknologi Bandung.
Jalan Ganesha 10.
Bandung 40132.
Indonesia Corresponding author: farisridma01@gmail.
Abstract Excavation and material stockpiling activities in the mining process cause a change in the distribution of forces and stresses in the material.
As a result, the material will seek a new equilibrium by releasing the load through a landslide.
As part of mitigation, it is necessary to monitor displacements occurring on slopes using high-accuracy devices.
Ground-based radar is a technique considered to have good ability to detect displacements in real-time.
However, ground-based radar has a weakness in that it cannot detect vertical displacement.
One of the emerging technologies that is used for monitoring vertical displacements as LoS displacements is D-InSAR analysis.
With the integration of both methods, displacements that occur horizontally and vertically on a slope can be detected properly.
In addition, the decrease in rock mass strength due to displacement can be predicted based on numerical analysis using the finite element method.
Monitoring was carried out from December 10th, 2021 to April 9th, 2022.
The monitoring results from the beginning to the end of the period showed that the horizontal and vertical displacements that occurred in the low wall area were 1247.
34 mm and 292.
5 mm.
Correlated with these conditions, the YoungAos modulus value of the rock mass decreased between 3% and 35%.
Keywords: D-InSAR.
ground-based radar.
rock mass.
YoungAos modulus.
Introduction Mining activities, whether carried out using an open pit or underground mining, cannot be separated from the activities of excavating and dumping materials.
The consequences of these activities can cause changes in the distribution of forces and stresses in the material.
As a result, the material will seek a new equilibrium by releasing the load through rock mass avalanches.
Instability in openings and slopes can cause various problems, from negative impact on the environment, damage to mining equipment, threatening the safety of workers, to mine closure due to economic aspects that are not on target .
As a part of a mitigation system, it is necessary to have good hazard risk management.
One form of management is by integrating various methods to monitor areas that are considered critical.
One of the parameters of concern is monitoring rock mass displacement.
Various methods are used in monitoring displacement, ranging from conventional methods by direct survey of critical areas, to methods that use high-tech devices to monitor rock mass displacement.
The conventional method is carried out by surveying the critical area.
A limitation of this method is that it is very dependent on weather conditions and can only be done during the day.
Therefore, many companies have now switched to using more modern devices, one of which is Ground-based radar.
This is one of the most effective and accurate tools for detecting rock mass displacement.
Displacement information is provided in a few minute intervals with accuracy up to the millimeter scale with wide area coverage .
Although it has many advantages, it has one major weakness, which is that it cannot detect vertical displacement.
Sensors ground-based radar can only detect the movement of rock masses that tend to approach or move away from the sensor.
Based on these conditions, it is necessary to integrate other devices or technologies that can detect vertical displacement with high accuracy.
One of the developing technologies is utilizing the Sentinel-1 satellite to determine changes in the earthAos surface condition through D-InSAR analysis by LoS displacement.
Another Copyright A2023 Published by IRCS - ITB
ISSN: 2337-5779
Eng.
Technol.
Sci.
Vol.
No.
3, 2023, 247-260
DOI: 10.
5614/j.
Research Paper Journal of Engineering and Technological Sciences Faris Ridwan Maulana et al.
advantage of D-InSAR is that it can detect rock mass displacements less than 5 mm .
In addition, the waves emitted by the satellite to detect displacements do not require a reflector device, so the coverage area in the DInSAR analysis is very wide.
With the integration of both techniques, the displacement values that occur horizontally and vertically can be detected properly.
Displacement in a rock mass is a form of decreasing the strength of the rock mass.
By using the rock mass displacement value, a decrease in rock mass strength can be obtained based on numerical analysis with the finite element method.
Information on displacements and the magnitude of the decrease in rock mass strength can be used as a reference in determining policies related to the safety and the operation of mining activities.
Literature Review Types of Slope Movement As a part of monitoring, it is important to know the mechanism of the rock mass displacement.
This information will be useful to determine which equipment to use and the interpretation of the results of monitoring devices.
Landslides that occur on slopes do not just happen at once.
The occurrence of rock mass landslides is preceded by a slow rock mass movement, which is then followed by a landslide event.
Therefore, it is important to know the characteristics of the rock mass movement occurring slowly.
The relationship between rock mass displacement against time can be divided into three types .
The first type of slope movement is called regressive movement.
The slope will tend to experience regressive movement, where the velocity of the rock mass displacement slows down until it reaches stability (Figure .
The slope will become more stable as the ratio between the driving force and the resisting force decreases.
The slope movement will continue to slow down until there are factors that can increase its velocity of movement.
These factors can be in the form of blasting or changes in groundwater pressure due to high rainfall intensity.
The progressive deformation trend is characterized by an increase in the velocity of the rock mass displacement.
However, not all increases in velocity will lead to progressive deformation trend.
The slope will only experience progressive if the discontinuity plane angle is steeper than the friction angle, which will cause the rock shear strength to decrease (Figure .
As shown in Figure 1, with increasing time, the C curve, which initially experienced a regressive deformation trend, will turn into a progressive one.
These changes in movement will cause a landslide.
These changes can be triggered by an increase in groundwater pressure, rock excavation on the toe of the slope, or other factors that cause instability of the slope.
Figure 1 Rock mass displacement curve against time .
Differential Interferometry Synthetic Aperture Radar (D-InSAR) Although radar technology has developed over fifty years, the components that support interferometric radar processing have been developed only about ten years ago.
One of the uses of remote sensing technology, in the world of mining, especially D-InSAR, is as a means of monitoring slope stability.
The principle used is to look for differences in the condition of the earthAos surface at the same location in different acquisition data (Figure .
this case, the radar sensor will read the line of sight (LoS) component that appears when a displacement occurs.
Integrated D-InSAR and ground-based radar for Open Pit Slope Stability Monitoring DOI: 10.
5614/j.
Figure 2 Illustration of D-InSAR analysis for surface displacements using line of sight (LoS) .
Information on displacements that occur on the earthAos surface can be known based on imagery produced by In the images, there is information on the phase difference of the electromagnetic waves emitted and then received back by the satellite.
Phase differences in electromagnetic waves occur due to changes in the conditions of the earthAos surface, as shown in Figure 2.
With targets that are only surface conditions, the coverage area of the monitoring method with D-InSAR is very wide.
In addition, electromagnetic waves emitted by satellites can penetrate clouds and are not affected by weather conditions in the local area .
Although it has several advantages, the D-InSAR monitoring method, has one major disadvantage.
Sentinel-1 data as raw data for D-InSAR analysis is available only at intervals of several days.
Therefore, takes several days to obtain displacement information from the satellite-based data acquisition process.
In addition, as a principle of monitoring with SAR interferometry, sensors on satellites can only find displacement information based on the trend of surface changes that are moving away from or approaching the sensor.
Thus, the displacements obtained by D-InSAR analysis are only displacements that occur vertically.
Methodology PT Adaro Indonesia is a coal mining company located in Tabalong Regency.
South Kalimantan (Figure .
The present study was conducted in the low wall area of pit South Tutupan, which is one of the areas included in the concession area of PT Adaro Indonesia.
Two other locations that were also included in the PT Adaro Indonesia concession area were Paringin and Wara.
The main product of PT Adaro Indonesia is Envirocoal, which is included in the sub-bituminous coal type with a calorific value in the range of 4000 to 5000 kcal/kg.
In 2021.
PT Adaro IndonesiaAos coal reserves amounted to 731 million tons with resources reaching 3.
3 billion ton.
In the same year.
PT Adaro IndonesiaAos coal production reached 43.
2 million tons with more than 70% coming from Tutupan.
Figure 3 Area of South Tutupan min: .
plane view.
camera directed to the north.
Because of high activities in the South Tutupan pit area, the company has determined various policies with high standards so that operations can run well and effectively.
One of the policies made is related to the aspect of slope stability.
Based on a survey the slope height on the low wall side of the South Tutupan pit exceeds 200 m.
This makes the slope stability of the South Tutupan pit low wall area one of the main focuses of the company.
Various methods and devices have been installed to monitor the stability of the low wall side slopes, one of Faris Ridwan Maulana et al.
which is a rock mass displacement monitoring device.
The displacement aspect is important because it is one of the indicators in determining policies related to mining operations.
One of the displacement monitoring devices that has been used and has good accuracy is ground-based radar because it can provide information on the horizontal displacement of rock masses in real-time with measurement accuracy up to millimeters .
Data from ground-based radar monitoring is presented in Table 1 for the period from December 10th, 2021 to April 9th, 2022 with a monitoring interval of two hours.
The graph in Figure 4 is a visualization of the horizontal displacement that occurred in the low wall pit area of South Tutupan.
Meanwhile, to identify LoS displacements in the low wall area.
D-InSAR analysis was used based on Sentinel-1 data processing.
The Sentinel-1 data used was open-source data with an interval of twelve days.
The monitoring period was based on the availability of data between the ground-based radar data and the Sentinel-1 data.
Therefore, monitoring was carried out from December 10th, 2021 to April 9th, 2022.
In general, the concept of monitoring with D-InSAR is carried out by comparing two Sentinel-1 data from successive acquisition periods.
there is a change in the topographic conditions, there will be a phase difference in the electromagnetic wave which is interpreted as LoS displacement.
To get the LoS displacement value, some processing must be done in the D-InSAR analysis using the Sentinel Application Platform (SNAP) software as illustrated in Figure 5 .
Figure 4 Correlation value and velocity of horizontal displacement in the low wall area.
Figure 5 Sentinel-1 processing for D-InSAR analysis .
With the Sentinel-1 data used in the D-InSAR analysis, the LoS displacement value is obtained by processing the data as illustrated in Figure 5.
One of the purposes of doing this processing is to remove noise arising from the Sentinel-1 data acquisition process.
Noise that appears during data acquisition comes from interactions between electromagnetic waves and atmospheric conditions, which can produce ambiguity in the SAR data.
This noise consists of two types, namely regional noise, and single-data anomalies.
Regional noise is noise that appears on a wide scale and usually occurs in several data acquisition periods.
One of the factors that causethis type of noise Integrated D-InSAR and ground-based radar for Open Pit Slope Stability Monitoring DOI: 10.
5614/j.
is the difference in thermal conditions on the surface due to changes in weather and season.
To detect the presence of this noise, a coherence value parameter is used.
Pixels that show high deformation values and low coherence values (< 0.
indicate that the detected deformation is affected by the presence of noise .
The type of noise that appears in one of the data .
ingle-data anomal.
can be identified based on the approach taken to regional noise.
If there are three or more data sets affected by high atmospheric conditions, the data can be replaced with data that has a low impact.
One of the rock mechanics parameters that affect rock mass displacement is YoungAos modulus.
A decreasing YoungAos modulus value correlates with decreasing rock mass strength.
To determine the decrease percentage, this study used finite element numerical modeling.
The modelling was done by reducing the value of YoungAos modulus gradually until the displacement value obtained from the numerical analysis was close to the actual displacement value from the monitoring results.
To carry out this modelling, data are needed such as slope geometry, subsurface lithology sections, mechanical parameters of rock mass, and load borne by the slope.
The slope geometry was created based on updated topographical conditions in the South Tutupan pit area, which were then made into a two-dimensional .
D) cross-section.
This cross-section was made transversely from the low wall to the high wall are with a length of more than 2500 m.
Based on these conditions, it was found that the slope height reached 250 m with an overall slope of 20A.
Then the subsurface lithology arrangement was input into the previously made cross-section.
The load experienced by the slope is divided into two types, namely the seismic load and the pore pressure coefficient (Ru coefficien.
Meanwhile, the rock types found in the South Tutupan pit area, especially in the low wall area, are three types of rock sandstone, mudstone, and coal with mechanical parameters as presented in Table 1.
When all the modeling parameters were known, the two-dimensional slope model presented in Figure 6 was obtained.
Table 1 Mechanical parameters of lithology in the South Tutupan pit area.
Mechanical Parameters Unit weight .
N/mA) UCS (MP.
YoungAos modulus (MP.
Poisson ratio Tensile strength (MP.
Cohesion (MP.
Friction angle (A) Sandstone 25 - 22.
80 - 11.
81 - 21209.
217 - 0.
0005 Ae 0.
18 Ae 0.
62 - 32.
Mudstone 00 - 24.
90 - 19.
51 - 15422.
261 - 0.
0014 Ae 0.
17 - 0.
68 - 27.
Coal 45 Ae 21.
25 Ae 12.
81 Ae 5284.
217 Ae 0.
0006 Ae 0.
17 Ae 0.
73 Ae 30.
AAo Figure 6 2D Cross-section A-AAo in low wall area.
Results and Discussion Horizontal Displacements Monitoring In this study, horizontal displacement monitoring was carried out using ground-based radar.
Monitoring results in the form of displacement values were obtained in real-time.
To get the same time interval as D-InSAR monitoring, a recap of the monitoring results was made with twelve-day intervals.
The results of monitoring with ground-based radar are shown in Table 2.
Faris Ridwan Maulana et al.
Table 2 Horizontal displacement recap based on ground-based radar monitoring.
Date
12/10/2021
12/22/2021
01/03/2022
01/15/2022
01/27/2022
02/08/2022
02/20/2022
03/04/2022
03/16/2022
03/28/2022
04/09/2022
Horizontal Displacement .
Based on the monitoring results, there was an increase in the horizontal displacement from December 10th, 2021 to January 15th, 2022.
On December 10th, 2021, the horizontal displacement value was 202.
28 mm, while the horizontal displacement value on January 15th, 2022 was 938.
53 mm.
This shows an increase in the displacement value of 736.
25 mm in an interval of 36 days.
In the following period, the increase in the displacement value tended to be linear and did not experience a significant increase.
The increase in the horizontal displacement value is shown in Figure 4.
Furthermore, to determine the effect of the horizontal displacement on the rock mass characteristics, a numerical analysis was carried out using the finite element method.
The analysis was carried out by varying the YoungAos modulus parameter for each type of material and applying the strength factor reduction (SRF) method.
Varying the value of YoungAos modulus is considered good if the horizontal displacement from the numerical analysis is close to the horizontal displacement from ground-based radar monitoring.
Based on these assumptions, the horizontal displacement value from the numerical analysis on December 10 th, 2021 was 281 mm, 205 mm in the low wall area and 76 mm in the high wall area with a reduction of the YoungAos modulus value of 3% compared to the initial value (Figure .
E-3%
E-35%
Figure 7 Horizontal displacement based on numerical analysis by reduction of YoungAos modulus and applying the SRF method.
Integrated D-InSAR and ground-based radar for Open Pit Slope Stability Monitoring DOI: 10.
5614/j.
In the same period, the horizontal displacement value from the monitoring of the ground-based radar was known to be 202.
28 mm.
This method was used to determine the decrease of the YoungAos modulus value based on the displacement that occurred in the next monitoring period.
The results obtained are presented in Table Table 3 Comparison of horizontal displacement based on ground-based radar monitoring and numerical analysis by reduction of the YoungAos modulus value and applying the SRF method.
Horizontal Displacement Monitoring .
Date
GBR
12/10/2021
12/22/2021
01/03/2022
01/15/2022
01/27/2022
02/08/2022
02/20/2022
03/04/2022
03/16/2022
03/28/2022
04/09/2022
Horizontal Displacement Numerical Critical Analysis .
SRF Total YoungAos modulus Reduction E-3% E-15% E-18% E-25% E-27% E-28% E-29% E-29.
E-30%
E-33%
E-35%
The regression results of the two horizontals displacement value produced an R-square value of 0.
9992 (Figure This shows that the horizontal displacement value of the results of numerical analysis and the results of monitoring were mutually correlated at 99.
In addition to the R-square value, from the regression results, other results were found by applying the following equation:
yccEaycuycyco = 1.
0034yccEayayaAycI where, yccEaycuycyco is the horizontal displacement value as a result of numerical analysis, and yccEayayaAycI is the horizontal displacement value based on ground-based radar monitoring.
Figure 8 Regression graph of horizontal displacement based on ground-based radar monitoring and numerical analysis by YoungAos modulus reduction and applying the SRF method.
Faris Ridwan Maulana et al.
LoS Displacements Monitoring To find the value of LoS displacement, the monitoring method used was D-InSAR analysis.
Like the monitoring period to determine the value of horizontal displacement, the monitoring period to find the amount of LoS displacement was December 10th, 2021 to April 9th, 2022 with an interval of twelve days.
Validation of the displacement values from the D-InSAR analysis was carried out on mining sequences.
This was done because prism monitoring data was not available.
Based on the mining sequence data during the monitoring period no mining was carried out in the area that was the monitoring target.
Thus, the displacement detected by D-InSAR was not affected by the mining sequences.
The LoS displacement at the beginning of the period was known to be 3 mm.
This value was obtained by comparing the topographic conditions in the Sentinel-1 data from December 10th, 2021 with the results from the previous data acquisition on November 28th, 2021.
The same was done for each twelve-day time interval until the end of the monitoring period.
On April 9th, 2022, the LoS displacement value was 50 mm.
The results of the D-InSAR analysis for the monitoring period of December 10th, 2021 to April 9th, 2022 are presented in Figure The LoS displacement value of the D-InSAR monitoring results is not an accumulation of displacement values.
To compare them to the horizontal horizontal displacement values, the LoS displacement values were All results of LoS displacement monitoring with D-InSAR are given in Table 4.
Figure 9 LoS displacement monitoring based on D-InSAR analysis from December 10th, 2021 to April 9th, 2022.
Integrated D-InSAR and ground-based radar for Open Pit Slope Stability Monitoring DOI: 10.
5614/j.
Figure 9 Continued.
LoS displacement monitoring based on D-InSAR analysis from December 10th, 2021 to April 9th, 2022.
Faris Ridwan Maulana et al.
Table 4 LoS displacement values based on D-InSAR analysis from December 10th, 2021 to April 9th, 2022.
Date
12/10/2021
12/22/2021
01/03/2022
01/15/2022
01/27/2022
02/08/2022
02/20/2022
03/04/2022
03/16/2022
03/28/2022
04/09/2022
LoS Displacements .
D-InSAR Cumulative To determine the reduction of the YoungAos modulus values based on the LoS displacement values, numerical analysis was carried out.
This numerical analysis was carried out with the same assumptions, namely variation of YoungAos modulus and the SRF method.
Based on the results of the variation of YoungAos modulus, the LoS displacement value was found by reducing the YoungAos modulus value by 3% from the actual value of 18 mm.
Variation of YoungAos modulus was continuously carried out until the LoS displacement values obtained from the numerical analysis were close to the results of D-InSAR monitoring.
Based on these assumptions, the largest LoS displacement value obtained from the results of numerical analysis with a 35% reduction in the value of YoungAos modulus was 284 mm (Figure .
The LoS displacement values resulting from the numerical analysis with the assumption of a complete reduction of YoungAos modulus are given in Table 5, while the results of the overall analysis in the form of a model are presented in Figure 10.
Furthermore, a comparison graph of the LoS displacement values obtained from the numerical analysis and D-InSAR monitoring is presented in Figure 11.
E-3%
E-35%
Figure 10 LoS displacement based on numerical analysis by reduction of YoungAos modulus and applying the SRF Integrated D-InSAR and ground-based radar for Open Pit Slope Stability Monitoring DOI: 10.
5614/j.
Table 5 Comparison of LoS displacement based on D-InSAR analysis and numerical analysis by reducing YoungAos modulus and applying the SRF method.
LoS Displacement Monitoring .
LoS Displacement SRF Date D-InSAR Cumulative Numerical Analysis .
Displacement 12/10/2021 12/22/2021 01/03/2022 01/15/2022 01/27/2022 02/08/2022 02/20/2022 03/04/2022 03/16/2022 03/28/2022 04/09/2022 Critical YoungAos Modulus
SRF
Reduction (%) E-3% E-15% E-18% E-25% E-27% E-28% E-29% E-29.
E-30%
E-33%
E-35%
Based on the LoS displacement values obtained from the results of the D-InSAR analysis and the numerical analysis, a regression analysis was performed to determine the correlation between both analyses.
The regression results are presented in Figure 11.
The regression results show an R-square value of 0.
This value indicates that the LoS displacement value from the numerical analysis and from D-InSAR monitoring were correlated at 98.
In addition, other results came from a regression in the form of the following equation:
yccycycuycyco = 1.
0051yccycyayaycuycIyaycI Where, yccycycuycyco are the LoS displacement values resulting from the numerical analysis and yccycyayaycuycIyaycI are the LoS displacement values based on the D-InSAR analysis.
Figure 11 Regression graph of LoS displacement based on D-InSAR analysis and numerical analysis by reducing YoungAos modulus and applying the SRF method.
To determine the level of accuracy of the D-InSAR analysis, the trend formed from the displacement values of the D-InSAR monitoring results was compared to the trend of the ground-based radar monitoring results as shown in Figure 12.
Overall, the two methods showed a constant or linear increase.
The displacement trend based on the ground-based radar monitoring showed a significant increase at the beginning of the period.
This may have been affected by noise in the ground-based radar data.
Furthermore, correlation of the two methods showed an R-square value of 85% based on the following equation:
yaycuycIyccycnycycy = 0.
Where the ground-based radar registered a horizontal displacement of around 1 mm, the D-InSAR analysis gave an LoS displacement of around 0.
1419 mm with an R-square of 85%.
Faris Ridwan Maulana et al.
Figure 12 .
Horizontal and LoS displacement charts.
correlation of horizontal and LoS displacements.
Discussion Based on the results of horizontal and LoS displacement monitoring, it was found that the displacement values In line with these conditions, the results of the numerical analysis also showed an increase in the displacement values, one of which was caused by a decrease in the YoungAos modulus value of the material.
the numerical analysis that was carried out, to obtain a displacement value that is close to the monitoring results, the YoungAos modulus value for each material was reduced by 3% to 35%.
As one of the parameters of rock mechanics, changes in YoungAos modulus are influenced by the heating-cooling cycles of the rock.
This can cause rock to experience deterioration so that it experiences a decrease in strength .
The intensity of light and high rainfall in areas located around the equator make this cycle faster.
In this study, the reduction in the YoungAos modulus value was calculated based on a monitoring period for horizontal and LoS displacement from December 10th, 2021 to April 9th, 2022, as fully described in Table 6.
the beginning of the monitoring period, the YoungAos modulus value had decreased by 3%, as is known from the results of the numerical analysis.
This value continued to increase until the end of the monitoring period, where the percentage of YoungAos modulus reduction was 35%.
Based on these assumptions, this study also conducted a regression to determine the decrease in the YoungAos modulus value as a function of time.
The results of the regression are presented in Figure 13, which shows an R-square value of 0.
In addition, an equation was obtained to determine the percentage of reduction in YoungAos modulus of the material in the South Tutupan pit at the following monitoring time intervals:
yaycyceycc (%) = Oe6.
Oe 0.
where, yaycyceycc (%) is the percentage of YoungAos modulus reduction, and yccycayc is the day of monitoring as the time.
Table 6 Percentage of YoungAos modulus reduction based on horizontal displacement and LoS displacement monitoring from December 10th, 2021 to April 9th, 2022.
Date
Day
10/12/2021
22/12/2021
03/01/2022
15/01/2022
27/01/2022
08/02/2022
20/02/2022
04/03/2022
16/03/2022
28/03/2022
09/04/2022
Horizontal Displacement .
LoS Displacement .
Numerical Analysis D-InSAR Analysis Numerical GBR Low wall High wall Sum D-InSAR Cumulative Analysis YungAos Modulus Reduction E-3% E-15% E-18% E-25% E-27% E-28% E-29% E-29.
E-30%
E-33%
E-35%
Integrated D-InSAR and ground-based radar for Open Pit Slope Stability Monitoring DOI: 10.
5614/j.
Figure 13 Regression graph of YoungAos modulus reduction against time.
Based on the YoungAos modulus reduction values in Table 6 and the regression results in Figure 13, it is known that the YoungAos modulus values decreased significantly at the beginning of the monitoring period.
In line with increasing time, the decreasing trend in YoungAos modulus was not significant until the end of the monitoring period with a decrease percentage value of 5%.
There was a decrease in the YoungAos modulus value in the South Tutupan Pit low wall area, caused among others by heating-cooling cycles.
This phenomenon causes several components in rocks to change.
One component that changes due to the heating-cooling cycle is the mineral composition of the rocks.
A type of mineral that can be found in large quantities in almost all types of rocks is Rocks with a predominantly quartz composition can experience a rapid change in YoungAos modulus at a temperature interval of 25AC to 520AC .
This shows when temperatures are not too high, changes in YoungAos modulus can occur.
When correlated with the conditions in the study area, the increase in temperature due to heat from sunlight could be one of the factors causing a decrease in the value of YoungAos modulus in the rocks.
In addition to changes the mineral contents, the heating-cooling cycle can also affect the presence of microcracks in rocks.
In sedimentary rocks, microcracks appear in the gaps between the grains that make up the If the rock undergoes a rapid heating-cooling cycle, it will cause the constituent grains to shrink.
This shrinkage will be the initiator of the formation of microcracks.
In line with the formation of more and more microcracks, the YoungAos modulus value of the rock will also decrease.
Conclusion To determine the displacement value, two methods were used, namely ground-based radar monitoring for horizontal displacement and D-InSAR analysis for LoS displacement.
Meanwhile, to determine the YoungAos modulus value of the rock mass, finite element numerical analysis was carried out by reducing the value of YoungAos modulus to obtain displacement value that are close to those from the monitoring results.
In line with the increase in the horizontal and LoS displacement values that were known based on monitoring, the results of the numerical analysis also showed a reduction of the YoungAos modulus value for each material.
Based on the results of the analysis, it was found that the YoungAos modulus value was reduced by 3% on December 10th, 2021 and continued decrease to 35% on April 9th, 2022.
Based on the results of the analysis, it is known that the decrease in the value of YoungAos modulus is related to heating-cooling cycles in the research area.
References