File Estimating Dewatering in An Underground Mine by Using 3d.en - Es
File Estimating Dewatering in An Underground Mine by Using 3d.en - Es
File Estimating Dewatering in An Underground Mine by Using 3d.en - Es
com
MÁSUNO
ARTÍCULO DE INVESTIGACIÓN
Litang HuIDENTIFICACIÓN
1,2 *, Menglin Zhang1,2, Zhengqiu Yang1,2, Yongfan3*, Jixiu Li4,
Hongliang Wang5,6, Celestin Lubale7
1Facultad de Ciencias del Agua, Universidad Normal de Beijing, Beijing, República Popular China,2Centro de Investigación de Ingeniería para el
Control y Remediación de la Contaminación de las Aguas Subterráneas del Ministerio de Educación, Universidad Normal de Beijing, Beijing,
República Popular China,3China ENFI Engineering Corporation, Beijing, RP China,4Jinchuan Group Co., Ltd, Beijing Road, distrito de Jinchuan,
ciudad de Jinchang, provincia de Gansu, República Popular China,5Instituto de Investigación de Ingeniería del Norte de China Co., Ltd., provincia
de Hebei, Shijiazhuang, República Popular China,6Centro de Innovación Tecnológica para la Seguridad de las Aguas Subterráneas Mineras de la
a1111111111 Provincia de Hebei, Shijiazhuang, República Popular China,7Ruashi Mining SAS, Sitio minero de Ruashi, Lubumbashi, Lualaba, República
a1111111111 Democrática del Congo
a1111111111
a1111111111 * [email protected] (LH);[email protected] (YF)
a1111111111
Abstracto
La entrada de agua subterránea a una mina subterránea afectará gravemente su plan minero y la seguridad de la
ACCESO ABIERTO ingeniería geológica. Los modelos de aguas subterráneas son herramientas poderosas que se utilizan
Citación:Hu L, Zhang M, Yang Z, Fan Y, Li J, Wang H, et al. comúnmente en las minas para desarrollar estrategias de deshidratación. Muchas minas en el área de Kolwezi han
(2020) Estimación de la deshidratación en una mina
estado presentes desde la década de 1950, y los patrones de flujo de agua subterránea se han visto
subterránea mediante el uso de un modelo 3D de
significativamente influenciados por las actividades mineras. Se desarrolla un plan minero para una mina
elementos finitos. MÁS UNO 15(10): e0239682.https://doi.
org/10.1371/journal.pone.0239682 subterránea con estratos sinclinales volcados en Kolwezi, Congo. Los modelos de aguas subterráneas anteriores
que utilizaban medios homogéneos en capas redujeron la precisión del modelo. Se formula un nuevo modelo
Editor:Jianguo Wang, Universidad de Minería y
Tecnología de China, CHINA tridimensional de agua subterránea que utiliza FEFLOW, que consiste en modelos geológicos regionales y locales
combinados que integran 16 secciones transversales hidrogeológicas y datos de registro de pozos, para predecir
Recibió:19 de abril de 2020
la deshidratación subterránea en el área de estudio. Se realizan pruebas de bombeo de 31 días con 3 pozos de
Aceptado:11 de septiembre de 2020
bombeo y 28 pozos de observación para estimar las propiedades hidrogeológicas. Los datos simulados del nivel
Publicado:29 de octubre de 2020
del agua coinciden bastante bien con los datos observados. Bajo 8 escenarios de posibles diseños de pozos, el
Derechos de autor:©2020 Hu et al. Este es un artículo de modelo predice una posible capacidad de drenaje superior a 23.900 m3/d en la etapa inicial de minería. El
acceso abierto distribuido bajo los términos delLicencia
concepto del modelo y su aplicación puede ser una referencia para otras minas con geología compleja para la
de atribución Creative Commons , que permite el uso,
seguridad minera en la región de interés.
distribución y reproducción sin restricciones en cualquier
medio, siempre que se cite al autor y la fuente originales.
Introducción
Declaración de disponibilidad de datos:El conjunto de datos de prueba de
La explotación de minerales y las actividades mineras desempeñan un papel vital en el desarrollo económico de la
financian investigaciones de hidrogeología y llevaron a Congo representó el 51% de la producción minera de cobre de África, y el sector de minería y procesamiento de minerales
cabo análisis de datos, metodología y administración de en la República Democrática del Congo representó aproximadamente el 22,1% de su producto interno bruto (PIB). Muchas
proyectos. Norte de China
de las operaciones mineras más importantes en la República Democrática del Congo consisten en la producción de cobre y
Engineering Investigation Institute Co., Ltd es
cobalto en la región del cinturón cuprífero de la provincia sureña de Katanga. El cinturón se extiende a lo largo de 250 km
responsable de las investigaciones de campo y llevó a
entre Kolwezi y Lubumbashi. La producción industrial de cobre comenzó en 1911 con la Union Minière du Haut Katanga
cabo la recopilación y el análisis de datos. China ENFI
Engineering Corporation es responsable de la (que se convirtió en Gécamines), y la mina Musonoi es conocida desde hace muchos años en la parte occidental del
diseño de deshidratador ng plan y análisis de datos. Este cinturón de cobre de Katanga.
La investigación también cuenta con el apoyo financiero parcial Durante la explotación de depósitos minerales, la deshidratación de los estratos del subsuelo es una de las tareas más
del Proyecto de la Fundación Nacional de Ciencias Naturales de
importantes para garantizar una operación minera geológicamente segura y eficiente.1 ]. Bajo condiciones geológicas
China (Número de subvención: 41877173) y el Proyecto principal
complejas, solución analítica [2 ,3 ] y los modelos de elementos analíticos no están disponibles y, por lo tanto, los modelos
de ciencia y tecnología de tratamiento y control de la
numéricos son herramientas eficientes para proporcionar un plan de deshidratación en áreas de medio poroso [4–7] y
contaminación del agua (Subvención No.
2018NX07109-003). Los financiadores no tuvieron ningún papel regiones fracturadas y kársticas [8–11]. La zonificación de los parámetros hidrogeológicos generalmente se realiza en base
en el diseño del estudio, la recopilación y análisis de datos, la a datos de pozo e interpolación de datos [12 ]. A diferencia de los suelos en capas [13 ], la zonificación es difícil de lograr
decisión de publicar o la preparación del manuscrito. mediante la interpolación de datos para entornos geológicos complejos, como sinclinales volcados. Los modelos
Conflicto de intereses:YF está afiliada a China ENFI geológicos 3D, que integran muchos tipos de datos de investigación de campo, son herramientas necesarias y
Engineering Corporation, JL está afiliada a Jinchuan generalmente utilizadas para demostrar la distribución de formaciones en las actividades mineras.14 ,15 ]. Es importante
Group Ltd, HW está afiliada a North China Engineering destacar que la heterogeneidad geológica está asociada con huecos mineros abiertos, fallas y fracturas tectónicas, y las
Investigation Institute Co., Ltd y CL está afiliada a Ruashi
redes estructuradas tradicionalmente tienen grandes dificultades para representar una alta heterogeneidad. Los modelos
Mining SAS. No existen patentes, productos en desarrollo
geológicos, a partir de software de modelado geológico existente como Petrel, GeoModeller y Leapfrog, proporcionan una
ni productos comercializados que declarar. Esto no altera
nuestro cumplimiento de las políticas de PLOS ONE sobre demostración visual y conocimiento de la distribución de formaciones, especialmente controlando acuíferos o fallas. Sin
el intercambio de datos y materiales. embargo, pocas investigaciones se preocupan por combinar modelos de flujo de agua subterránea con modelos
geológicos 3D en acuíferos complicados. MOD-FLUJO-USG [dieciséis ] y FEFLOW 7 [17 ] han sido desarrollados para
describir condiciones geológicas complejas con sistemas de cuadrícula no estructurados. Mientras tanto, todavía existen
incertidumbres en la geometría de los acuíferos y sus características debido a las investigaciones limitadas y la complejidad
de las condiciones geológicas. Hassen y Gibsonet al. [18 ] construir un modelo geológico 3D utilizando el software
GeoModeller y evaluar la influencia de la conectividad del acuífero y el papel de las fallas y la geología en el flujo de agua
subterránea dentro del sistema acuífero; sin embargo, hasta donde sabemos, se reportan pocos modelos de flujo de agua
subterránea con un modelo geológico acoplado. El acoplamiento del modelo de flujo y el modelo geológico es eficaz para
mejorar la precisión del modelo numérico porque proporciona una visualización visual de las condiciones geológicas y del
patrón de flujo del agua subterránea y, por lo tanto, una estimación precisa de la deshidratación minimizará los riesgos
La megabrecha de Kolwezi contiene depósitos de Cu-Co alojados en estructuras plegadas y frágiles fracturadas
del Subgrupo de Minas. Como lo indican el mapa de superficie y las secciones transversales [19 –21 ], estas
estructuras están alineadas más o menos paralelas a la orientación EW/60—SW-NE de la napa de empuje [22 ,23 ].
La mina Musonoi es una mina subterránea recientemente planificada en el depósito de cobre de Kolwezi,
República Democrática del Congo. Los estratos sinclinales de la mina Musonoi están volcados, con condiciones
geológicas e hidrogeológicas complejas. Varias minas están presentes en el área de Kolwezi nappe, incluido el tajo
abierto KOV (Kamoto, Oliveira y Virgule) [24 ], la mina subterránea Kamoto, Kamoto East, el antiguo tajo abierto
Musonoi, el tajo abierto Dikuluwe, los tajos abiertos Mashamaba West y East, la mina Sicomines y la mina de cobre
Kolwezi. Las intensas actividades mineras alteran el camino natural del flujo de agua subterránea y fortalecen las
conexiones de agua subterránea entre diferentes minas en la mina de cobre Kolwezi. Las intensas actividades
mineras alteran el camino natural del flujo de agua subterránea y fortalecen las conexiones de agua subterránea
entre diferentes minas en el área de Kolwezi Nappe. Se han desarrollado modelos de flujo de agua subterránea
para explorar las estrategias de deshidratación en el área de Kolwezi nappe. El primer modelo informático público
del flujo de aguas subterráneas [25 ] se realiza para los tajos Dikuluwe y Mashamaba West e indica que el progreso
de la deshidratación en 1985 no sería satisfactorio para los planes mineros existentes. Africa Geo-Environmental
Services (Pty) Ltd. desarrolla modelos de flujo de agua subterránea que consideran la deshidratación de KOV,
Kamoto East, Kamoto Underground, Dikuluwe, Mashamaba East y Mashamaba West en 2006. SLR Consulting
(Sudáfrica) (Pty) Ltd. toma la iniciativa El área de Kolwezi Nappe como área modelo para estimar la deshidratación
en la mina Musonoi recientemente planificada en 2013 y 2017 utilizando el software FEFLOW. Dichos modelos
simplifican los estratos en propiedades en capas, y las simplificaciones afectan la precisión del modelo hasta cierto
punto.
Aunque se llevan a cabo muchas investigaciones de campo en el área minera local, pocas están presentes a
escala regional. Sin embargo, se requiere que el límite del modelo se extienda más allá de la mina.
área para evitar las influencias de las actividades mineras adyacentes y los límites creados por el hombre en la
precisión del modelo de flujo de agua subterránea. Entonces, los principales objetivos de esta investigación son
combinar el modelo de flujo de agua subterránea con el modelo geológico regional y local, y luego explorar la
deshidratación en minas subterráneas utilizando el modelo establecido. Primero, se diseña un modelo conceptual
hidrogeológico para las actividades mineras subterráneas. En segundo lugar, el modelo geológico 3D se construye
a escalas regional y local y así se llevará a cabo la parametrización. Luego, se prepararán los datos y luego se
calibrará el modelo de flujo de agua subterránea. Por último, el modelo establecido se utiliza para predecir la
deshidratación en la mina Musonoi. Los resultados de este artículo se pueden utilizar eficazmente como referencia
para desarrollar un programa de deshidratación eficiente.
material y métodos
Área de estudio
El área de estudio está ubicada al sur del ecuador en la meseta de Katanga en la República Democrática del Congo (Figura
1 ). El área de estudio tiene un clima de sabana. La temperatura media anual es de aproximadamente 21,2 ˚C, la
temperatura más alta es de aproximadamente 31,9 ˚C y la temperatura más baja es de aproximadamente 31,9 ˚C.
Fig 1. Ubicación de la zona de estudio.Imagen creada por los autores en QGIS 3.10 y Visio 2016; no se utilizó ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g001
Serie (De joven a mayor) Formación Nombre local Breve descripción Espesor aproximado (m)
Kundelungu kú Sedimentos 3000–5000
Nguba ng Arenisca, esquisto 200–500
Alto Roan R4 Mwashya esquisto, limolita, arenisca, dolomitas 50-100
R3-2 dipeta esquistos arenosos alrededor de 1000
Modelo conceptual
Debido a las actividades mineras a largo plazo, los caminos naturales del flujo de agua subterránea en cada área de la
mina se alteran significativamente. Las conexiones hidráulicas del sistema de agua subterránea en cada área minera están
mucho más cerca del proceso minero debido a la intensa actividad minera. Por lo tanto, no es razonable establecer un
modelo numérico del flujo de agua subterránea en una determinada área de la mina y, basándose en el estudio del
modelo de elevación digital del área de estudio, se utilizan las cajas de herramientas de hidrología en el software QGIS
para dividir el límite de la cuenca. El límite del modelo, mostrado enFigura 1 , considera la estructura del cuerpo de la napa
límites artificiales. Los límites oeste y sur se establecen como divisiones naturales de cuencas
hidrográficas y se representan como límites sin flujo. Los límites norte y este del área del modelo se
establecen como límites de flujo y se representan en el modelo como límites de Neumann. Como
indican las perforaciones, la profundidad de la fractura y el desarrollo tectónico se extiende hasta
1.000 m por debajo de la superficie. El fondo del modelo se establece a una profundidad de 2000 m
debajo de la superficie y, por lo tanto, el fondo puede representarse como un límite sin flujo. Para
tener en cuenta la complejidad del cambio de estrato, el área de estudio se divide en 18 capas de
estructura hidrogeológica. El espesor total de la capa estudiada es de 2000 m, y los espesores de las
capas divididas son, en orden desde la superficie del suelo, 50 m, 50 m, 50 m, 50 m, 50 m, 50 m, 50
m, 50 m. , 50m, 50m, 100m, 100m, 100m, 100m, 100m, 200m, 300m y 500m. El límite superior recibe
principalmente la infiltración de la precipitación y los flujos de agua subterránea hacia las minas.
Modelo de desarrollo
Diagrama de flujo del estudio.
El diagrama de flujo se muestra enFigura 2 . En primer lugar, se define el área del modelo y se establece el modelo
conceptual. En segundo lugar, el área del modelo se discretiza usando FEFLOW [26 ] software en sistemas de red y se
utiliza un modelo de flujo de agua subterránea en capas. En tercer lugar, se preparan todos los datos, incluidos los
registros de pozos, los datos de las pruebas de bombeo, el modelo geológico local y los informes hidrogeológicos, y luego
se desarrollará un modelo geológico 3D a nivel regional. A continuación, se recopilan y utilizan otros datos del modelo
para llevar a cabo la calibración del modelo mediante pruebas de bombeo en el área de Musonoi. Finalmente, el modelo
Fig 2. Diagrama de flujo del desarrollo del modelo.Imagen creada por los autores en Visio 2016; no se utilizó ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g002
Fig. 3.Cambios litológicos en la llanura (A) y las secciones transversales (B) en el área de la napa de Kolwezi. Imagen creada por los autores en QGIS 3.10, Datamine Studio RM 1.2.45.0 y Visio 2016;
no se utilizó ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g003
Fuentes y sumideros
En el límite superior del sistema, el acuífero recibe recarga principalmente de las precipitaciones, los ríos y las
aguas residuales de los depósitos de relaves. Los elementos de descarga de aguas subterráneas incluyen
principalmente drenaje minero, drenaje a pozos abandonados, drenaje al río, fugas de pozos y evaporación de
aguas subterráneas. La evaporación del agua subterránea ocurre en las áreas de poca profundidad del agua
subterránea. La recarga del acuífero proviene principalmente de las precipitaciones durante la temporada de
lluvias. El coeficiente de infiltración de la precipitación está estrechamente relacionado con la precipitación, la
estructura litológica de la zona de aireación y la profundidad del nivel freático. Según estudios previos, el
coeficiente de recarga de infiltración de lluvia generalmente está entre el 10% y el 25%, pero es solo del 3% en la
zona de roca RAT, del 75% en algunos tajos abiertos y hasta del 25% en otras areniscas, cráteres de desechos y
estanques de relaves. . En el modelo, la tasa de recarga de infiltración se calcula multiplicando la infiltración de
precipitación promedio anual por el coeficiente de infiltración. Los datos de observación de precipitaciones son
monitoreados por DRC Huakan Geological Services Co., Ltd. Las extracciones de agua subterránea de diferentes
minas son elementos importantes del drenaje de agua subterránea y tienen influencias significativas en el nivel
del agua subterránea. El flujo de deshidratación promedio fue de aproximadamente 5,52×104metro3/ d en el tajo
Dikuluwe en septiembre de 2018, y 4,22×104metro3/d en el tajo oeste de Mashambo en septiembre de 2018. La
deshidratación en estos tajos abiertos se representa como el límite de altura variable en el tiempo en el software
FEFLOW, y la altura se establece según la profundidad actual hasta la superficie en cada área de la mina.
Los pozos de observación se distribuyen en el área modelo (Figura 1 ). Para analizar las condiciones
hidrogeológicas del área de la mina Musonoi y estimar con precisión la deshidratación de la mina Musonoi, DRC
Huakan Geological Services Co., Ltd lleva a cabo una prueba de bombeo a gran escala. Las pruebas de bombeo
comienzan a las 8:00 am del mes de noviembre. 22 de diciembre de 2018, y parará a las 8:00 am del 18 de
diciembre de 2018. El nivel freático restablecido se observa desde las 8:00 pm del 23 de diciembre de 2018. Se
cuenta con 3 pozos de bombeo, PW01, PW02 y PW03, y las producciones de cada pozo de bombeo son 1232.4,
3532.32 y 2790.64 m3/d, respectivamente (figura 4 ). Hay 28 pozos de observación (incluidos 3 pozos de bombeo)
en el área de interés. El 18 de diciembre de 2018 se observa la máxima reducción entre los pozos,
aproximadamente 61 m en el pozo PW01, 58 m en el pozo PW03 y 45 m en el pozo PW02. Las profundidades de
perforación de cada pozo están marcadas enfigura 4 , y todos los pozos tienen varias capas.
El dominio modelo tiene un área de aproximadamente 623,63 km2. El área de estudio se divide en 12.740
nodos y 24.988 elementos, con un total de 18 capas de modelo (el número total de nodos es
242.060 y el número de elementos es 449.784). En el área de estudio, se refina la región de Kolwezi nappe y
se refina nuevamente el área de la mina Musonoi. Los pozos de observación y los pozos de bombeo se
consideran puntos fijos y la malla se divide según la forma de la red fluvial. La estación seca en el área de
estudio es de mayo a septiembre, y el período de calibración del modelo se establece de septiembre a
diciembre. La distribución inicial del nivel del agua se calcula basándose en la simulación del flujo de agua
subterránea en estado estacionario y luego se ajusta después de la entrada de datos del modelo [27 ]. La
calibración del modelo se realiza con base en los parámetros hidrogeológicos estimados, el método de
prueba y error y la herramienta FePEST. Los pozos multicapa se representaron como el módulo de pozo
multicapa en FEFLOW 7. El modelo adopta un control automático de paso de tiempo a través del esquema
de integración de tiempo hacia adelante Adams-Bashforth/trapezoide hacia atrás (AB/TR).
Fig 4. Distribuciones de pozos de observación en el área de interés, el número alrededor del pozo significa la profundidad de los pozos.Imagen creada por los autores en QGIS 3.10 y Visio 2016; no se utilizó ningún
material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g004
modelo de geología 3D
Los datos de elevación del terreno utilizados en este modelo se obtienen mediante datos de elevación de alta resolución
del software BIGEMAP (http://www.bigemap.com/ ) con una resolución horizontal de 7 metros, que están disponibles
públicamente. Para el modelo de geología regional, el rango de profundidad se toma desde la superficie del suelo hasta los
1. Establecer un modelo geológico regional en el área de estudio. En primer lugar, se digitalizan el mapa geológico plano
y 16 secciones transversales. Luego, el cuerpo geológico se genera mediante la integración de datos 3D y el software
Datamine RM de todos los mapas geológicos. Se derivan los clips de horizonte de los perfiles en las 18 profundidades
Figura 5.Demostración del modelo geológico en el área del modelo y el área de la mina Musonoi, (A), (B), (C) y (D) significa malla 3D en el área de estudio, secciones transversales en el área de estudio, malla
3D en el área de interés y cortes transversales en el área de interés, respectivamente. Imagen creada por los autores en Tecplot 360 Ex 2014R1 y Visio 2016; no se utilizó ningún material protegido por
derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g005
Tabla 2. Estimación de la conductividad hidráulica a lo largo del eje principal de anisotropía larga en la llanura para cada parámetro de zonificación de cada capa del modelo en el área de
estudio (unidad: m/d).
litología Nizilo_Block kú ng Mwashya dipeta CMN Dakota del Sur kamoto RATA Brecha
Capa 1 0.0050 0.1307 0.0200 0.0200 0.0579 0.4803 0.3194 0.7500 0.0500 2.4000
Capa 2 0.0050 0.1307 0.0200 0.0200 0.0579 0.4803 0.3194 0.7500 0.0500 2.4000
Capa 3–4 0.0100 0.1307 0.0200 0.0200 0.1000 1.2000 0.4500 1.2000 0,2000 3.0000
Capa 5–8 0.0100 0.1307 0.0200 0.0200 0.0800 0.6500 0.3500 1.0500 0.1200 1.2000
Capa 9-11 0.0100 0.0984 0.0100 0.0100 0.0434 0.3600 0.2396 0.5625 0.0150 1.2000
Capa 12-13 0.0050 0.0660 0.0050 0.0050 0.0289 0.2396 0.1597 0.3750 0.0150 1.2000
Capa 14-15 0.0050 0.0660 0.0050 0.0050 0.0289 0.2396 0.1597 0.3750 0.0250 1.2000
Capa 16-18 0.0050 0.0050 0.0050 0.0050 0.0289 0.2396 0.1597 0.3750 0.0250 1.2000
https://doi.org/10.1371/journal.pone.0239682.t002
El pozo de bombeo del PW01 y el pozo de observación están a 130 m del pozo PW01 al pozo GM11 y a 44 m
del pozo PW01 al pozo TB01, lo que sugiere que la anisotropía aún está presente en el área de estudio y
debe ser considerada. De acuerdo con la distribución de la formación, el ángulo entre la dirección del eje
principal de anisotropía de la conductividad hidráulica y el eje del sistema de coordenadas cartesiano se
generaliza como el mismo valor en el área de estudio, y los ángulos de los ejes principales son 60, 20 y 20
grados desde elX,yyzejes, respectivamente. La conductividad hidráulica del eje principal de anisotropía
larga en la llanura se estima mediante los datos de las pruebas de bombeo de los pozos de ingeniería en
las áreas mineras vecinas, y el valor para cada capa del modelo se muestra enTabla 2 . La conductividad
hidráulica del eje principal de anisotropía corta en la llanura se establece igual que la del eje largo. La
conductividad hidráulica del eje principal de anisotropía vertical se establece en 1/10 de la del eje mayor. El
rendimiento específico varía de 0,01 a 0,05 y se establece en 0,02, 0,02, 0,04 y 0,05 para las zonas CMN, SD,
Kamoto y Breccia, respectivamente. La capacidad de almacenamiento varió de 0,00001 a 0,0001 m-1y se
establece en 0,0002, 0,0002, 0,0002 y 0,0001 m-1para las zonas CMN, SD, Kamoto y Breccia (zona
fracturada), respectivamente.
Resultados
Conductividad hidráulica (k) y los límites de recarga, como la precipitación y la infiltración de agua
superficial, son factores críticos para la calibración del modelo, y es necesario evaluar sus incertidumbres.
En este estudio, la sensibilidad local se calcula con base en el método de variación de un solo factor. Los
principales factores incluyen las conductividades hidráulicas en la zona Ku, la zona RAT, la zona Kamoto,
la zona SD, la zona CMN y la zona Dipeta, y la infiltración de ríos y embalses. El escenario base (S0) se
establece como parámetros estimados iniciales. Los escenarios se establecen mediante un factor con el
25%, 75%, 125% y 150% de los parámetros iniciales (que se muestran enTabla 3 ). la simulación
https://doi.org/10.1371/journal.pone.0239682.t003
Fig 6. Sensibilidad local de diferentes zonas geológicas y condiciones de recarga ante cambios en el nivel del agua.Imagen creada por los autores en Python y
Visio 2016; no se utilizó ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g006
El período se fija desde las 8:00 am del 10 de septiembre de 2018, hasta las 8:00 am del 18 de diciembre de 2018. El índice
de sensibilidad (Eq (1 )) se define como la raíz del error cuadrático medio de la diferencia simulada del nivel del agua.
1PAGnorte
SAsi¼ j¼1 dhsidjÞ hs0djÞÞ2 d1Þ
norte
El período de calibración se establece desde las 8:00 am del 10 de septiembre de 2018 hasta las 8:00 am del 18 de
diciembre de 2018, para un total de 104 días. Las pruebas de bombeo se realizan desde las 8:00 am del 22 de noviembre de
2018 hasta las 8:00 am del 18 de diciembre de 2018. Los principales objetivos de ajuste son los datos del nivel freático.
Xng X
GRAMOdpag1;pag2;���;pagnorteÞ ¼ ohhdi; jÞ � ½Sdi; jÞ Sgramodi; jÞ�2 d2Þ
i¼1 j¼ntc
Tabla 4. Estadísticas de las diferencias absolutas entre la reducción observada y simulada entre 28 pozos.
Diferencia absoluta entre los niveles de agua subterránea ΔS(m) en la mina Musonoi ΔS�0,5 0,5<ΔS�1.0 1.0<ΔS�5.0 ΔS>5.0 total
Fig 7. Comparación de la reducción simulada y observada en 12 pozos de observación.Imagen creada por los autores en Excel 2016 y Visio 2016; no se utilizó ningún material
protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g007
Discusión
Cambio del nivel freático y distribución del abatimiento máximo.
Por lo general, la altura inicial de cada nodo se interpola a partir de las alturas observadas en los pozos de
observación. Hay 18 capas de modelo en el área de estudio y, por lo tanto, es difícil obtener la altura inicial
para todas las capas de modelo porque los pozos no están bien distribuidos en cada capa y sobre el área
de estudio. La altura inicial se simula mediante el método de parámetros y cabeza acoplados [27 ]. El cono
de depresión del nivel freático está claramente simulado y visto en el área del modelo. En el área de
interés, el agua subterránea fluye desde los límites este y sur. Debido a que la distancia entre el área de
interés y las minas adyacentes es corta, el patrón de flujo del agua subterránea probablemente se verá
alterado después de las intensas actividades mineras en el área de interés, y no está claro si la división del
agua subterránea está presente entre las dos minas.
Fig 8. Comparación de la reducción máxima durante la prueba de bombeo para 6 pozos.Imagen creada por los autores en Excel 2016; no se utilizó
ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g008
La distribución de la reducción máxima durante las pruebas de bombeo puede proporcionar una comprensión
profunda de las vías dominantes del flujo de agua subterránea. El desarrollo del cono de depresión en respuesta al
bombeo está sujeto principalmente al área de las zonas de fractura controladas por fallas (aproximadamente a
una profundidad de aproximadamente 150 a 350 m). Aunque los pozos de observación de múltiples capas no
tienen las mismas profundidades que los pozos de sondeo, las reducciones máximas de todos los pozos se
demuestran en la llanura (figura 9 ). La distribución de la reducción máxima es muy desigual. El eje largo de
máxima caída es de aproximadamente 45̊ al norte del este, y la longitud de influencia es de aproximadamente
1,50 km. El eje corto de máxima caída está aproximadamente a 45° al norte del oeste y la longitud de influencia es
de aproximadamente 1,0 km. Los pozos de observación a una profundidad vertical de aproximadamente 150 a 350
m tienen la mayor reducción y luego los pozos a poca profundidad tienen la menor reducción.
En el área de interés, la reducción inducida por la deshidratación es tan grande que generalmente ocurre el
secado del acuífero. Durante la calibración se observa el fenómeno de secado y rehumedecimiento de la capa
modelo (Capa 1) (Higo 10 ). Bajo la influencia del pozo de bombeo PW02, el cambio del nivel del agua subterránea
con el tiempo es normal en el sitio P1 pero anormal en el sitio P2 (aquí las ubicaciones de P1 y P2 están muy cerca
de PW02). La interfaz de capa entre la capa 1 del modelo y
Tabla 5. Conductividad hidráulica calibrada a lo largo del eje principal de anisotropía larga en la llanura para cada parámetro de zonificación de cada capa del modelo en el área de estudio
(unidad: m/d).
litología Nizilo_Block kú ng Mwashya dipeta CMN Dakota del Sur kamoto RATA Brecha
Capa 1 0.0015 0.0172 0.0029 0.0021 0.0132 0.1112 0.0784 0.1078 0.0463 0.1012
Capa 2 0.0009 0.0122 0.0022 0.0016 0.0083 0.0554 0.0420 0.0557 0.0240 0,1284
Capa 3–4 0.0018 0.0127 0.0041 0.0030 0.0136 0.0788 0.0068 0.0169 0.0199 0.1836
Capa 5 0.0012 0.0167 0.0030 0.0023 0.0202 0.0866 0.0217 0.0182 0.0252 0.2823
Capa 6 0.0012 0.0182 0.0030 0.0023 0.0639 0,1232 0.0358 0.0272 0.0365 0.4288
Capa 7–8 0.0012 0.0162 0.0032 0.0023 0.0361 0.0503 0.0339 0.0226 0.0272 0,1849
Capa 9-11 0.0012 0.0120 0.0021 0.0013 0.0145 0.0371 0.0283 0.0298 0.0126 0.0630
Capa 12-13 0.0007 0.0104 0.0014 0.0007 0.0158 0.0214 0.0180 0.0198 0.0126 0.0654
Capa 14-15 0.0006 0.0081 0.0007 0.0006 0.0084 0.0113 0.0130 0.0129 0.0048 0.0463
https://doi.org/10.1371/journal.pone.0239682.t005
Fig 9. Mapa de contorno de máxima caída en la zona de interés.Imagen creada por los autores en QGIS 3.10 y Visio 2016; no se utilizó ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g009
2 está marcado enHigo 10 ; Se presentan dos saltos bruscos del nivel del agua al inicio y al final de las
pruebas de bombeo, lo que suele indicar problemas de convergencia y baja eficiencia de cálculo.
Fig 10. Secado y rehumedecimiento observados de la capa modelo.Imagen creada por los autores en Excel 2016; no se utilizó ningún material
protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g010
Según el borrador del plan minero, la primera fase durante el período de servicio de la mina del área de interés es
a una profundidad de 140 ml y completamente deshidratada mediante pozos de bombeo. Durante el período de
prueba de bombeo se perforaron tres pozos de bombeo, a saber, PW01, PW02 y PW03. El área del cono de
depresión bajo los pozos de bombeo actuales es relativamente pequeña y se deben agregar pozos adicionales
para ampliar el área. El plan de deshidratación de 140 ml se logrará en 1 año y, por lo tanto, se consideran 8
escenarios (enumerados enTabla 6 ). La profundidad de los pozos PW01, PW02 y PW03 es de 310,20 m, 251,51 m y
325,00 m, respectivamente. El diámetro de los tres pozos varía con la profundidad desde 0,780 m en la parte
superior hasta 0,378 m en la parte inferior. El diámetro de la tubería de la pared del pozo en el fondo es de 0,273
m para tres pozos de bombeo. En el modelo, el radio del pozo se establece uniformemente en 0,189 m para tres
pozos de bombeo. En diferentes escenarios para la deshidratación del pozo de 140 ml, el radio de todos los pozos
recién agregados se establece en 0,189 m y la profundidad de estos pozos se establece en 300 m. Los escenarios
S1 y S2 suman los pozos PB01, NW101, NW102 y PB02, y las tasas de bombeo son como se muestran en la tabla.
Los escenarios S3 y S4 suman los pozos PB01, NW201, NW202, NW203 y PB02 con caudales de bombeo de 2000 m.
3/día, 2500 m3/día, 2500 m3/d y 2500 m3/d y 1500 m3/día, 1500 m3/día, 1500 m3/d y 2000 m3/d, respectivamente.
Los pozos PB01, NA01, NA02, NA03, NA04, NA05 y PB02 con caudales de bombeo de 2000 m3/día, 2500 m3/día,
2500 m3/día, 2500 m3/día, 2500 m3/día, 1500 m3/día, 1500 m3/día, 2500 m3/día, 1500 m3/d y 2000 m3/d, se añaden
en los escenarios S5 y S6. Los escenarios S7 y S8 suman los pozos PB01, NA01, NW201, NA02, NW202, NW203 y
PB02. La distribución de los pozos de achique en los seis escenarios se muestra enFigura 11 .
Los resultados simulados bajo los ocho escenarios se muestran enFigura 12 . Se requiere un cono de
depresión con una caída de 140 m para cubrir el área del área de la mina (como lo indica el color rojo en
Figura 12 ). Las áreas del cono de depresión en los ocho escenarios se encuentran básicamente dentro del
área del yacimiento, pero las áreas del cono son diferentes para los escenarios. El mapa de contorno del
patrón de descenso es una forma elíptica con un eje principal de noreste a suroeste. El escenario S5 tiene
el área más grande del cono, mientras que el escenario S1 tiene el área más pequeña. Como se puede ver
en la figura, los escenarios S1, S2, S3, S4, S5, S6, S7 y S8 pueden formar un área de descenso de cono
continuo con 80 m, 80 m, 90 m, 80 m, 150 m, 140 m, 120 m y 110 m, respectivamente. El flujo total de
deshidratación para el escenario S1, S2, S3, S4, S5, S6, S7 y S8, es de 16.400 m3/día, 17.100 m3/día, 18.900
m3/día, 15.900 m3/día, 28.400 m3/día, 22.800 m3/día, 23.900 m3/d y 18.900 m3/d, respectivamente. El
escenario S5 tiene el flujo de deshidratación más alto, mientras que el escenario S4 tiene el flujo más bajo.
Aunque el área del cono de depresión formado por el escenario S5 es grande, la inversión para perforación
será mayor. Por lo tanto, se recomienda el escenario S7 al considerar el área de deshidratación y el costo
de perforación y operación. Se puede ver en la figura que el escenario S7 (Figura 12 ) puede formar un
cono de depresión con una caída de 120 m en el medio de un pozo de 140 ml después de un año con una
caída máxima de 320 m, y una caída de 200 m en las regiones locales.
Con base en los patrones simulados de cambios de reducción para el escenario S7, se puede encontrar a partir
de los resultados que el cono de depresión en la sexta capa del modelo es el más grande y la reducción máxima es
de casi 320 m. La reducción máxima disminuye desde la sexta capa de modelo a la primera capa de modelo, y
también desde la sexta capa de modelo a la octava. De las pruebas de bombeo existentes, la tasa máxima de
producción de los pozos PW03, PW01 y PW02 es de 2700 m3/ día, 1200 m3/d y 3500 m3/d, y la cantidad de tasa de
producción para pozos adicionales enTabla 6 puede ser razonable. El volumen total de agua de todos los pozos en
la sección media del pozo de 140 ml es de aproximadamente 23.900 m3/d para el escenario S7. Así, la capacidad de
deshidratación con 23.900 m3Se requerirá /d al menos para el plan de deshidratación.
Las pruebas de bombeo se realizan durante la estación seca, donde la infiltración es de aproximadamente 2,19 mm/
día. Las diferentes estaciones, como la estación normal y la estación húmeda, también tendrán una influencia importante
en la deshidratación. El modelo también se utiliza para evaluar la influencia de la estación húmeda y normal en el cambio
del nivel del agua subterránea. La precipitación es de 1952 mm/año y 1200 mm/año para la temporada húmeda y la
temporada normal, respectivamente, y luego se utiliza como entrada del modelo. Una mayor infiltración aumentará el
nivel del agua subterránea. Para el escenario S7, en relación con la estación seca, la diferencia de descenso es de
actividades de deshidratación alterarán las condiciones de presión del agua y, entonces, puede ocurrir un colapso desde el
punto de vista de la seguridad geológica. Se debe prestar mucha atención a estos problemas de seguridad geológica
Fig 11. Diseño de pozos de drenaje diseñados para pozo de 140 mL.Imagen creada por los autores en QGIS 3.10 y Visio 2016; no se utilizó ningún
material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g011
Figura 12.Mapa de contorno de reducción después de 1 año bajo un eje de 140 ml para 8 escenarios: (A), (B), (C), (D), (E), (F), (G) y (H) representan el escenario 1 , 2,
3, 4, 5, 6, 7 y 8. Imagen creada por los autores en QGIS 3.10 y Visio 2016; no se utilizó ningún material protegido por derechos de autor.
https://doi.org/10.1371/journal.pone.0239682.g012
deshidratación, el cambio de los niveles de agua subterránea debe ser monitoreado durante la implementación del proyecto de
deshidratación.
Debido a la alta heterogeneidad y anisotropía de las formaciones geológicas, los resultados simulados
tienen incertidumbres. Los principales efectos de las incertidumbres son los siguientes:
1. Incertidumbre de las formaciones geológicas. Aunque este artículo establece el modelo geológico 3D
acoplando el modelo geológico 3D regional basado en secciones transversales geológicas con un modelo
geológico refinado local basado en datos de pozo, pueden existir fracturas locales y las distribuciones de las
fracturas aún se desconocen debido a las limitadas investigaciones de campo.
2. Incertidumbre de los parámetros hidrogeológicos. Los datos de las pruebas de bombeo en el área de interés se utilizan
para calibrar el modelo. Sin embargo, los tres pozos de bombeo están ubicados a profundidades de 250 a 320 m
debajo de la superficie del suelo, y los pozos de observación no están bien distribuidos en el área y en dirección
vertical. La duración limitada de una prueba de bombeo (menos de 3 meses en este estudio) generalmente no es
suficiente para rastrear la estacionalidad y los cambios en las cabezas de agua subterránea. A partir de la calibración de
los datos de las pruebas de bombeo en el área de interés, las propiedades del acuífero obtenidas pueden no ser
representativas de toda la extensión de las capas y el área del modelo.
Conclusiones
Este estudio recopila muchos datos, como precipitación, escorrentía de ríos, nivel de agua subterránea y los
resultados de pruebas de bombeo, y construye un modelo tridimensional de flujo de agua subterránea en un
área de interés con un modelo geológico 3D utilizado para ayudar en la predicción de la deshidratación. de minas
subterráneas. La calibración del modelo se realiza utilizando datos de pruebas de bombeo y se analizan las
características del flujo de agua subterránea. Se obtienen las siguientes conclusiones.
El interés es pequeño, el área del modelo se extiende hasta los límites de la cuenca para minimizar el
impacto de las condiciones de contorno del modelo de agua subterránea en la precisión debido a la
influencia de intensas actividades mineras. Verticalmente, la profundidad de 2000 m bajo la superficie
del suelo se considera la capa impermeable inferior del modelo.
3. La calibración del modelo se realiza utilizando datos de prueba de bombeo. A partir de las curvas de ajuste de los niveles de agua
subterránea, los resultados del modelo son relativamente buenos y podrían proporcionar cambios en el nivel de agua
Este estudio proporciona un método científico para la construcción de modelos y la predicción inicial de la
deshidratación. Sin embargo, los resultados del modelo pueden tener algunas incertidumbres. Debido a la
privacidad de las minas, es difícil obtener datos detallados, como datos precisos de drenaje fuera del área de
interés, y esto puede reducir la precisión del modelo hasta cierto punto. Durante el desarrollo actual de la rampa
Musonoi, se han recopilado datos de fracturas y vetas locales desde el frente o la pared y estos parámetros serán
útiles para contribuir al modelo geológico local. Especifique la conductividad hidráulica de la capa RSC, que
generalmente tiene, junto con el CMN, la transmisividad y la abundancia de agua más fuertes debajo de las
zonas de meteorización. El modelo geológico puede mejorarse a escalas mucho más finas si se finalizan las
investigaciones en curso a pequeña escala.
Información de soporte
Archivo S1.
(XL)
Archivo S2.
(DOCX)
Expresiones de gratitud
Contribuciones de autor
Curación de datos:Menglin Zhang, Hongliang Wang, Celestin Lubale.
Investigación:Hongliang Wang.
Metodología:Litang Hu, Zhengqiu Yang, Yong Fan, Jixiu Li, Celestin Lubale.
Referencias
1.Yihdego Y, Salem HS, Ayongaba B, Veljkovic Z. Desafíos del sector minero en los países en desarrollo,
Tigray, Etiopía e inspiradoras historias de éxito de Australia. Revista Internacional de Minería e
Ingeniería Mineral 2018; 9(4): 321–367.
2.Yihdego Y, Paffard A. Predicción del flujo de entrada y la profundidad de recuperación de una mina a cielo abierto en el soum
Durvuljin, provincia de Zavkhan, Mongolia. Agua de Mina y Medio Ambiente 2017; 36: 114-123.
5.Murdoch LC, Germanovich LN, Wang H, et al. Hidrogeología de las proximidades de la mina Homestake, Dakota del
Sur, EE.UU. Revista de Hidrogeología 2012; 20: 27–43.
6.Jiang SM, Kong XL, Ye H, et al. Optimización de la deshidratación de aguas subterráneas en la mina de carbón a cielo abierto
Shengli número 1, Mongolia Interior, China. Ciencias Ambientales de la Tierra 2013; 69: 187–196.
7.Bahrami S, Ardejani FD, Aslani S, Baafi E. Modelado numérico del flujo de agua subterránea a una mina a cielo
abierto en avance: pozo Kolahdarvazeh, Irán central. Monitoreo y Evaluación Ambiental 2014; 186: 8573–
8585.https://doi.org/10.1007/s10661-014-4025-x PMI:25186026
8.Yihdego Y, Danis C, Paffard A. Simulación numérica del flujo de agua subterránea en 3D para discontinuidades
geológicas en la cuenca de Unkheltseg, Mongolia. Ciencias Ambientales de la Tierra 2015; 73 (8): 4119–4133.
9.Hu LT, Jiao JJ. Modelado de las influencias de la recuperación de tierras en los sistemas de aguas subterráneas: un estudio de
caso en la península de Shekou, Shenzhen, China. Ingeniería Geológica 2010; 114: 144-153.
10.Schrader A, Winde F, Erasmus E. Uso de los impactos de la minería a profundidad para investigar la hidrología kárstica: un
enfoque basado en Darcy para predecir el futuro de los manantiales dolomíticos secos en el campo aurífero de Far West
Rand (Sudáfrica). Parte 1: un modelo conceptual de recarga y flujo intercompartimental. Ambiental
Ciencias de la Tierra 2014; 72: 3549–3565.
11.Bodeux S, Pujades E, Orban P, et al. Interacciones entre el agua subterránea y la cavidad de una antigua mina de pizarra utilizada
como depósito inferior de una UPSH (Hidroelectricidad de almacenamiento por bombeo subterráneo): un enfoque de
modelización. Ingeniería Geológica 2017; 217: 71–80.
12.Mu WP, Wu Q, Xing Y, et al. Uso de simulación numérica para la predicción de la deshidratación de minas de un sistema de agua
kárstica subyacente a la veta de carbón en la cuenca de Yuxian, en el norte de China. Ciencias Ambientales de la Tierra 2018;
77: 215.
13.Pujades E, SunordesteEV, Carrera J, et al. Deshidratación de una excavación profunda realizada en un suelo estratificado.
Ingeniería Geológica 2014; 178: 15-27.
15.Garzonio CA, Piccinini L, Gargini A. Modelado de aguas subterráneas de acuíferos fracturados en minas: el estudio
de caso de Gavorrano (Toscana, Italia). Mecánica de Rocas e Ingeniería de Rocas 2014; 47: 905–921.
18.Hassen I, Gibson H, Hamzaouiazaza F, et al. Modelado geológico 3D del sistema acuífero de Kasserine, Túnez central: nuevos
conocimientos sobre la geometría del acuífero y las interconexiones para una mejor evaluación de los recursos de agua
subterránea. Revista de Hidrología 2016; 539: 223–236.
19.Mendelsohn. La geología del cinturón de cobre de Rodesia del Norte. Londres: Macdonald. 1961, págs. 351–405.
20.François A. La extremidad occidental del arco cuprifmire shabien Etude geologique. Bureau d'études
géologiques. Gécamines-Exploitation, Likasi, Zaire. Aulhenlie consultoría de inversiones (China) Lo. Ltd.
Traducción en 2006 (en chino), 1973.
21.François A. Carte géoogique de la region de Kolwezi-Kalukundi (Shaba, República de Zaire). Bulletin de
la société belge de Geologie 1980; 89: 141-143.
22.Cailteux JLH, Muchez P, Cuyper JD, et al. Origen de las megabrechas en el cinturón de cobre de Katanga (RD Congo).
Revista de Ciencias de la Tierra Africanas 2018; 140: 76–93.
23.Nkulu CBL, Casas L, Haufroid V, et al. Sostenibilidad de la minería artesanal de cobalto en la República Democrática del Congo.
Sostenibilidad de la Naturaleza 2018; 1: 495–504.https://doi.org/10.1038/s41893-018-0139-4 PMI:30288453
24.Mukonki PM. Planificación estratégica de minas: un análisis FODA aplicado a la mina a cielo abierto KOV en
la República Democrática del Congo. En Academia Mundial de Ciencia, Ingeniería y Tecnología Revista
Internacional de Ingeniería Geológica y Ambiental. 2017; Vol: 11, Núm. 11.
25.Straskraba V, Placet J, Holubec M. Aplicación de modelado por computadora para el diseño de deshidratación de minas a cielo
abierto. Procedimiento IMWA de 1985; 579–597.
26.Diersch HJG. FEFLOW: modelado de elementos finitos de flujo, masa y transporte de calor en medios porosos
y fracturados. vol. XXXV. Springer, Berlín Heidelberg, 2014, págs.996.https://doi.org/10.1111/gwat. 12285
PMI:25393965
27.Chen CX, Pei SP, Jiao JJ. Hundimiento del terreno causado por la explotación de aguas subterráneas en la ciudad de Suzhou, China.
Revista de Hidrogeología 2003; 11: 275–287.