Estimating The Rate Of Post-Mining Filling Of Pit Lakes

Introduction Deep open-pit mines invariably affect the local and regional hydrologic systems. Pit dewatering, occurring during mining operations, puts an obvious hydrologic stress on these hydrologic systems. However, post-mining hydrologic impacts resulting from the pit refilling with groundwater following the cessation of mining activity can also be significant. The prediction of the rate at which the post-mining pit will fill with groundwater is a critical aspect of assessing the long-term hydrologic impacts. Numerical groundwater flow modeling provides a method for predicting the groundwater refilling rate of the pit. The rate at which pit "lakes" fill depends on several factors: •the rate and duration of pit dewatering; •the depth and size of the ultimate pit and •the pre-mining hydrologic regime. These factors can be incorporated into a detailed numerical groundwater flow model that can then be used to assess the effects of dewatering and post-mining recovery on the local and regional hydrologic systems. A sufficiently detailed, numerical groundwater model provides the oportunity to: •account for complex geology near the pit; •assess the impact of active pit dewatering and •predict the long-term impacts of post-mining groundwater flow into the pit. A detailed groundwater model incorporating these items has been developed and applied at an operating open-pit mine. Developed by Durbin and O'Brien (1987), the three-dimensional, finite-element, groundwater flow model was used to represent the hydrologic system of an approximately 253-km2 (98-sq mile) area surrounding the pit. Historical groundwater elevation data, stream flows and meteorologic, geologic and geophysical data were used to establish the dimensions and initial conditions for the model. Steady-state conditions, representing the pre-mining local and regional hydrologic systems, were simulated using the initial conditions incorporated into the groundwater model. The groundwater model was then utilized to simulate various dewatering programs, to predict the filling rate and the groundwater depth in the ultimate pit once mining activities are complete and to assess the long-term impacts on the regional groundwater flow system. Development of pit lake model Groundwater modeling efforts were completed in two phases. The first focused on pit dewatering activities, while the second phase concentrated on the post-mining effects on the hydrologic system. The final estimates of groundwater elevations calculated during the pit dewatering simulations were used in predicting the post-mining recovery of the hydrologic system. The groundwater model was also modified prior to the second phase to account for the volume of rock removed during mining activities. To account for the actual volume of rock mined, the geometry of the post-dewatering model grid was modified to approximate the final pit geometry. The depth and width of the ultimate pit were divided into eight idealized stages that represented significant changes in the bench geometries. These eight stages were then introduced sequentially into the model according to the predicted water elevations within the pit. In this way, changes in the volume and depth of water within the pit were accounted for through time. Once the ultimate pit geometry was accounted for in the model, it was necessary to assign new hydraulic characteristics to those parts of the model grid (elements) that represented excavated rock. The solution of the numerical model requires that finite hydraulic conductivity values be assigned to the portion of the groundwater model that represents excavated rock. Therefore, the calculated groundwater elevations differ, somewhat, between the edges and the center of the open pit. These model-calculated water elevations at the edge and in the middle of the open pit represent the elevation of water that would occur in the pit lake. To minimize the error in the estimated level of water within the pit lake, the hydraulic conductivity was increased to a value that would: •minimize the predicted difference between the groundwater elevations across the open pit and •produce a numerically stable solution. Specific storage is the hydrologic parameter that accounts for the water produced by compaction of the aquifer matrix. To predict the groundwater volume that would flow into the ultimate pit, this parameter was assigned a value equivalent to the compressibility of water. This value of specific storage reflects the post-mining groundwater storage occurring as an open body of water. Additionally, a specific yield of 1.0 was assigned to the pit elements to represent the 100% porosity of the open pit. In