Vol. 18, No.

1: 108-116
Enero - junio 2022
ISSN 1680-8894
Portal de revistas: http://revistas.utp.ac.pa ISSN-L 2219-6714

Combinando el método de refracción sísmica y análisis multicanal

de ondas superficiales para la detección de niveles freáticos someros
en el noroccidente del Istmo de Panamá
Combining the seismic refraction method and multichannel analysis of surface waves
for detection of shallow water-table in the northwestern Isthmus of Panama
Alexis Mojica 1,2* , Osiris Cubilla 3 , Eduardo Silva 3 , José Ladrón-De Guevara 4 , Alberto Ruíz 5 , Carlos A. Ho 1 ,
Jaime Toral 6
1
Universidad Tecnológica de Panamá, Laboratorio de Investigación en Ingeniería y Ciencias Aplicadas, Centro Experimental de
Ingeniería, Panamá
2
Sistema Nacional de Investigación (SNI), SENACYT, Panamá.
3
3P Sociedad Ltda., Panamá.
4
Universidad Tecnológica de Panamá, Facultad de Ciencias y Tecnología, Panamá.
5
Fundación INDICRI, Panamá.
6
Universidad Nacional Autónoma de Chiriquí, Panamá.

alexis.mojica@utp.ac.pa; osiris.cubilla.1986@gmail.com; silvasan@cwpanama.net; cheoldg.13@gmail.com; einstein40@gmail.com;

carlos.ho@utp.ac.pa; jaime.toral@gmail.com
Fecha de recepción: 27 de abril de 2022; Fecha de aprobación: 6 de mayo de 2022.
*Autor de correspondencia: Alexis Mojica (alexis.mojica@utp.ac.pa)

RESUMEN. En el año 2017, se realizaron una serie de estudios geofísicos en la comunidad de Finca 30, provincia de Bocas del
Toro, noroccidente de Panamá y cuyo objetivo se focalizó en evaluar el uso combinado de los resultados obtenidos a través de la
refracción sísmica (velocidad de onda P) y el análisis multicanal de ondas superficiales (velocidad de onda S) vía el cociente de
Poisson para la detección de capas freáticas someras en el sitio de estudio. En esta comunidad se estableció un perfil de 30 m de
longitud en donde se utilizó un sismógrafo para estudios de onda P y S, con 16 geófonos de alta frecuencia para el primero y 12 de
baja frecuencia para el segundo y adicionalmente, un prototipo de inyección de corriente alterna de 400 V con un arreglo tipo
Schlumberger a fin de verificar los resultados de la prospección sísmica; el procedimiento para estudios de onda P se basaron en la
norma ASTM D5777-18, para el de onda S en el documento sobre el proyecto InterPACIFIC, y para el sondeo eléctrico, la norma
ASTM D6431-18. El rango de valores del cociente de Poisson obtenido (0.47-0.49), y 10.7 ohm.m a una profundidad de 1.8 m
sugieren la existencia de una zona saturada la cual fue corroborada a través de una perforación desarrollada en el sitio. Gracias a
estos resultados, los acuíferos superficiales podrían ser aprovechados para abastecimiento de la población y para evaluar el nivel de
afectación de la capacidad portante para futuros proyectos de ingeniería civil.

Palabras clave. Análisis multicanal de ondas superficiales, Bocas del Toro, coeficiente de Poisson, nivel freático, refracción
sísmica, resistividad eléctrica.

ABSTRACT. In 2017, a series of geophysical surveys were carried out in the community of Finca 30, province of Bocas del
Toro, north-western Panama and whose aim was focused on evaluating the combined use of the results obtained through seismic
refraction (P-wave velocity) and multichannel analysis of surface waves (S-wave velocity) via Poisson's ratio for the detection of
shallow water-table at the site. In this community, a 30-m-long profile was established where a seismograph was used for P and S
wave studies, with 16 high-frequency geophones for the first and 12 low-frequency geophones for the second, and additionally, a
AC electrical current injection prototype of 400 V with a Schlumberger type array in order to verify the results of the seismic survey;
the procedure for P-wave studies was based on the ASTM D5777-18 standard, for the S-wave study in document of InterPACIFIC
project, and for the electrical sounding, the ASTM D6431-18 standard. The range of Poisson's ratio values obtained (0.47-0.49), and
10.7 ohm.m at 1.8 m of depth suggest the existence of a saturated zone which was corroborated through a borehole developed at the
 Alexis Mojica | Osiris Cubilla | Eduardo Silva |José Ladrón-De Guevara | Alberto Ruíz | Carlos A. Ho | Jaime Toral

site. Thanks to these results, the shallow aquifers could be used to supply the population and to evaluate the level of affectation of
the bearing capacity for future civil engineering projects.

Keywords. Multichannel analysis of surface waves, Bocas del Toro, Poisson ratio, water-table, seismic refraction, electrical
resistivity.

1. Introducción
El desarrollo de proyectos de ingeniería civil en áreas empleando el time-term inversion method, y la inversión
de suelos compuestos por sedimentos no consolidados y de la curva de dispersión de velocidad de fase de las
saturados, podrían representar un riesgo significativo ya ondas superficiales a través del algoritmo neighborhood;
que pueden verse afectados de forma natural debido no de forma paralela, un sondeo eléctrico fue realizado en el
solamente a la amplificación sísmica y licuación que sitio de interés.
experimentan al encontrarse en una región sísmicamente
activa sino también por la presencia de acuíferos 2. Materiales y métodos
superficiales. Ante este tipo de escenarios, los métodos 2.1. Descripción del sitio de interés
geofísicos pueden jugar un papel importante debido a la El sector noroccidental del Istmo de Panamá
rápida adquisición de los datos en campo, la cobertura en (específicamente en la provincia de Bocas del Toro) se
extensión y su alto nivel de confiabilidad, así como el caracteriza por presentar desde finales del siglo
bajo costo que representa; hecho que se traduce en una antepasado, una alta inversión y desarrollo de la
optimización de recursos económico y humano. En productividad agrícola bananera y gran parte de las
estudios hidrogeológicos, los métodos eléctricos y comunidades circundantes a estas zonas se construyeron
electromagnéticos han sido los más empleados [1], [2] y sobre rellenos constituidos por sedimentos,
en Panamá, se han desarrollado algunos proyectos que específicamente por limo.
contemplan la aplicación de estas técnicas a estudios Además, la existencia de acuíferos superficiales en
hidrogeológicos e hidrológicos [3-10]; sin embargo, los estas zonas constituye un factor natural que podría afectar
métodos sísmicos no han sido ampliamente utilizados en la capacidad portante de este tipo de suelos para futuros
este tipo de estudios. proyectos de ingeniería civil. Las características físicas
En muchos casos, el uso independiente de la del suelo y la fácil accesibilidad al sitio fueron factores
refracción sísmica para la obtención de modelos de que influyeron en el desarrollo de este trabajo. Un mapa
terreno en base a las velocidades de onda de compresión de ubicación del sitio de estudio se presenta en la Figura
u ondas P (Vp), y del análisis multicanal de ondas 1, con coordenadas 328,908mE y 1,042,979mN.
superficiales MASW para obtener modelos de terreno en
función de las velocidades de las ondas de cizalla u ondas 2.2. Aspectos geológicos
S (Vs), no conllevan a un amplio conocimiento de los A pesar de que la zona de estudio se ubica dentro de
parámetros que definen un entorno hidrogeológico dado; la Formación Almirante Lajas, la cual domina los
no obstante, el análisis de las razones Vp/Vs y Poisson, sectores de media montaña y el litoral costero, sobre yace
podrían ser clave para la detección del nivel freático de a formaciones sedimentarias Gatún Uscari, en
acuíferos superficiales [11], [12], [13]. complemento con algunos cuerpos intrusivos de la
El objetivo de este trabajo se centra en aplicar por formación Guayabito sobre la cordillera esta zona
primera vez esta nueva técnica de métodos combinados a costera, la zona se caracteriza por presentar aluviones,
fin de detectar el nivel freático en un área ubicada al sedimentos no-consolidados, areniscas, corales,
sureste de la ciudad de Changuinola, provincia de Bocas manglares, conglomerados, lutitas carbonosas,
del Toro. Para esto, se utilizarán los resultados de los deposiciones tipo delta, propios de la Formación Las
procesos de inversión correspondiente al análisis de los Lajas, del período Cuaternario (ver la Figura 1) además
tiempos de llegada de las ondas refractadas y directas

RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022. 109
 Combinando el método de refracción sísmica y análisis multicanal de ondas superficiales para la detección de niveles freáticos someros en el noroccidente del Istmo de Panamá
Combining the seismic refraction method and multichannel analysis of surface waves for detection of shallow water-table in the northwestern Isthmus of Panama

de ser zona de influencia y depósito de material eruptivo y durante la jornada geofísica el nivel freático se detectó
reciente del cercano volcán Barú. a una profundidad de 1.80m.
Por otro lado, los ensayos de humedad natural del
suelo presentan un incremento porcentual para las
muestras obtenidas entre 0 y 0.6m (46.12%), para las
muestras obtenidas entre 1.5 y 2.1m de profundidad, la
humedad natural del suelo aumentó a 48.46%. Las
muestras obtenidas entre las profundidades de 3.0 y 3.6m
presentaron por su parte una disminución en el valor de
este parámetro (45.41%) y finalmente, para el conjunto
de muestras obtenidas entre 4.5 y 5.1m, la humedad
natural reportada fue de 46.12%. La Figura 2(c) presenta
un gráfico sobre este resultado.

Figura 1. Ubicación geográfica del sector de Finca 30 en el mapa geológico y

alrededores, provincia de Bocas del Toro.
Fuente: Dirección Nacional de Recursos Minerales [14].

Otro elemento importante en la geología es la

presencia de grandes fallas regionales de corrimiento que
afectan la superficie con una extensión Noroeste y
buzamientos predominantes hacia el Noreste producto de
las tensiones que ejerce la placa del Caribe en el
desarrollo del grabben Centroamericano del cual forma
parte.

2.3. Información estratigráfica del sitio

En la Figura 2(a) se presenta una imagen aérea de la
zona bajo estudio con el perfil de 30m en donde se Figura 2. (a) Imagen aérea de la zona de interés y el perfil establecido, (b)
desarrollaron todas las pruebas geofísicas (dirección detalle estratigráfico de la perforación desarrollada cercano al centro del perfil
y (c) contenido natural de humedad de las muestras obtenidas en dicha
Noroeste-Sureste). Por su parte, la Figura 2(b) muestra la perforación.
descripción estratigráfica en el centro del perfil Fuente: Imagen aérea tomada de Google Earth.
establecido. Un conjunto de muestras fue tomado de esta
perforación para su posterior análisis de laboratorio. Los 2.4. Estrategia de adquisición
resultados mostraron dos horizontes: el primero 2.4.1. El sondeo eléctrico – inversión 1D
correspondiente a un limo elástico de plasticidad media, La configuración electródica tipo Schlumberger ha
pH neutro y consistencia medianamente firme hasta una sido empleada en diversas investigaciones relacionadas
profundidad de 2.10m, seguido de un segundo horizonte con la caracterización estratigráfica y la búsqueda de
compuesto por un limo gris oscuro, caracterizado por fuentes de aguas subterráneas [14], [15], [16]. Para el
presentar un elevado contenido de humedad, baja sondeo eléctrico desarrollado en el punto central del
plasticidad, pH neutro y consistencia medianamente perfil en el área de estudio, se empleó un prototipo de
firme. La perforación concluyó a la profundidad de 6m, inyección de corriente alterna; los valores de resistividad
eléctrica aparente se obtuvieron a través de la medición
de los voltajes registrados, de la intensidad de la corriente

110 RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022.
 Alexis Mojica | Osiris Cubilla | Eduardo Silva |José Ladrón-De Guevara | Alberto Ruíz | Carlos A. Ho | Jaime Toral

inyectada y la geometría del arreglo electródico una placa metálica ubicada en un determinado punto a lo
empleado. Un total de 10 pares de datos (resistividad largo del perfil. El incremento en el número de golpes en
eléctrica aparente-distancia AB/2) fueron obtenidos para un punto dado aumenta la razón señal/ruido; un total de
dos diferentes separaciones MN/2. La Figura 3(a) 5 puntos de impacto fueron establecidos a lo largo del
presenta la distribución de los electrodos de corriente y perfil tal como lo muestra la Figura 3(b).
voltaje para la separación electródica MN/2 = 0.3m El sistema fue configurado para: (a) una duración
(arriba) y MN/2 = 1.0m (abajo), ubicándose el punto temporal de las trazas sísmicas de 250ms, (b) 50ms de
central del sondeo cercano a la perforación. Del total de pre-evento (delay) y (c) un intervalo de muestreo de
datos medidos en campo, 9 fueron sometidos al proceso 0.250ms. La Figura 4 ilustra las 16 trazas sísmicas
de inversión ya que para la distancia AB/2 = 4m, se obtenidas para cada uno de los 5 shot point establecidos
registró el voltaje para las posiciones MN/2 = 0.3m y en el perfil bajo estudio.
1.0m.

Figura 3. Bosquejo de la configuración geométrica para la adquisición de

datos de (a) resistividad eléctrica y sísmicos para estudios de onda (b) P y (c)
S. En la posición 15m se encuentra la perforación.

Los datos fueron invertidos empleado el software

IX1D de Interpex y cuyo código de inversión está basado
en el método Ridge Regression [17], [18]; su estimador
interpreta las curvas de resistividad eléctrica aparente
para modelos de terreno horizontalmente estratificados,
considerando una pequeña cantidad de eigenvalores
contenidos en el esquema de inversión y compuestos por
las matrices Jacobiana e identidad.

2.4.2. Prospección de onda P – refracción símica

Un total de 16 geófonos de 14Hz fueron separados 2m
entre ellos obteniéndose una longitud total de perfil de
30m tal como lo muestra la Figura 3(b), en donde las
Figura 4. Sismogramas registrados por los geófonos de 14Hz y obtenidos para
posiciones del primer y último geófono fueron de 0 y los 5 shot points o fuentes establecidas a lo largo del perfil.
30m, respectivamente.
Estos geófonos se interconectaron a un sismógrafo tipo La secuencia del procesamiento de estas señales de
ES-3000 de Geometrics Inc. a través de un sistema multi- refracción sísmica consistió en la determinación del
cable, los cuales fueron empleados también para los tiempo de llegada de las ondas P a cada geófono; para
estudios de onda S; por su parte las ondas sísmicas fueron esto, se empleó el módulo Pickwin de SeisImager2D
generadas a través de 5 golpes de un mazo de 20lb sobre (Geometrics). Seguidamente se empleó el módulo

RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022. 111
 Combinando el método de refracción sísmica y análisis multicanal de ondas superficiales para la detección de niveles freáticos someros en el noroccidente del Istmo de Panamá
Combining the seismic refraction method and multichannel analysis of surface waves for detection of shallow water-table in the northwestern Isthmus of Panama

Plotrefa para definir las dromocronas y luego obtener un fase de las ondas Rayleigh a diferentes frecuencias, lo
modelo 2D de terreno a través de un proceso de inversión que conlleva a la generación de un espectro de velocidad
basado en el método time-term inversión [19]. Este es un de fase del cual se extrae una curva de dispersión.
método simple de inversión de tiempos de trayecto de Tomando en cuenta que el centro del arreglo de los
una onda sísmica, el cual fue desarrollado por [20]. geófonos ubicados a lo largo del perfil corresponderá a la
posición del sondeo sísmico, se procedió con la
2.4.3. Prospección de onda S – (Análisis Multicanal de extracción de la curva de dispersión normalizando la
Ondas Superficiales 1D o MASW-1D) amplitud del espectro para cada frecuencia dada y para
Sobre el mismo perfil, se utilizaron 12 geófonos de cada uno de los dos puntos de impacto establecidos;
4.5Hz los cuales se separaron 2m entre ellos, posteriormente, se apilaron los dos espectros de
obteniéndose una longitud del arreglo de 22m; las dispersión normalizados y se procedió en identificar los
posiciones del primer y último geófono fueron 4 y 26m, puntos de la curva correspondiente al modo fundamental.
respectivamente. Un total de dos puntos de impacto se Este primer procesamiento se realizó a través del
establecieron en los extremos de dicho perfil (-1 m y 31 programa de código abierto Geopsy para pruebas de
m). La configuración del sistema consistió en establecer: fuente activa (MASW-1D).
(a) una extensión temporal de los registros de 1000ms, En este mismo punto, la selección de la frecuencia
(b) un intervalo de muestreo de 0.5ms y (c) un total de 5 mínima en el espectro apilado depende tanto de la
golpes en cada punto de impacto de forma análoga a la longitud máxima del arreglo de geófonos (L) como de la
prueba precedente. longitud de onda máxima. [21] y [22] recomiendan que
La Figura 3(c) presenta la distribución de los esta longitud de onda máxima debe ser inferior a L debido
geófonos y puntos de impacto para el desarrollo de esta al efecto de campo cercano y la sobreestimación del valor
prueba MASW. La Figura 5 muestra las 12 trazas de la velocidad a baja frecuencia. En este estudio, la
sísmicas obtenidas para cada uno de los puntos de frecuencia mínima obtenida fue de 11.6Hz la cual se
impacto sobre el perfil dado. El análisis multi-canal de encuentra asociada con una longitud máxima de 10.6m.
ondas superficiales se basa en el estudio de la propiedad De acuerdo con [21] la profundidad de investigación
de dispersión de las ondas Rayleigh en las trazas sísmicas corresponde a la mitad del valor de la longitud de onda
obtenidas y de ahí, la posibilidad de determinar un perfil máxima (λmax/2 = 5.3m). Por otro lado, la frecuencia
de velocidad de onda S en función de la profundidad. máxima depende del efecto aliasing o distorsiones que se
presentan en el espectro de dispersión normalizado. Para
este estudio, dicha frecuencia fue de 24.4Hz.
La inversión del conjunto de datos de velocidad de
fase de las ondas Rayleigh vs la frecuencia se desarrolló
empleando el algoritmo neighborhood el cual fue
desarrollado por [23] y emplea los constructores
geométricos o celdas de Voronoi para derivar la
búsqueda en el espacio de los parámetros; dicho
algoritmo fue adaptado por [24] y [25] a estudios
subsuperficiales a través de la herramienta Dinver
incorporado en Geopsy.

3. Resultados y Discusión
3.1. Resultado de la inversión eléctrica 1D
Figura 5. Sismogramas registrados por los geófonos de 4.5Hz y obtenidos para Para este primer ensayo geofísico, las observaciones
los 2 shot points o fuentes establecidas a lo largo del perfil.
o datos de campo (valores de resistividad eléctrica
Las ondas sísmicas de baja frecuencia pueden llegar a aparente medidos en campo) y la respuesta del modelo
agitar la columna de suelo de mayor espesor cercana a la (valores de resistividad eléctrica aparente calculados) se
superficie, hecho que permite determinar la velocidad de presentan en el gráfico log-log de la Figura 6(a).

112 RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022.
 Alexis Mojica | Osiris Cubilla | Eduardo Silva |José Ladrón-De Guevara | Alberto Ruíz | Carlos A. Ho | Jaime Toral

Como puede observarse en el gráfico bilogarítmico, la

respuesta del modelo se ajusta a los datos de campo al
establecer un total de 2 capas; la resolución del problema
inverso 1D dio por resultado el modelo unidimensional
que se presenta en la Figura 6(b) en conjunto con los
modelos equivalentes, obteniéndose un error RMS de
2.73%. La primera capa se encuentra caracterizada por
presentar un valor de resistividad eléctrica calculada
igual a 37.9ohm.m y un espesor de 1.7m, mientras que el
segundo horizonte presenta una resistividad más baja
(10.7ohm.m); como puede observarse, el contraste
existente entre ambos horizontes se aproxima al nivel Figura 7. (a) Curvas de tiempo de arribo de las ondas refractas y (b) modelo
freático detectado durante la jornada de trabajo (1.8m). 2D de terreno obtenido con el método de inversión del tiempo de recorrido,
n.f. corresponde al nivel freático.

Un análisis de los resultados obtenidos en la Figura

7(a) sugieren la presencia de dos capas: un horizonte
superficial de un espesor que varía entre 1.9 y 1.5m, y
con una velocidad de onda P calculada de 163m/s,
seguido de un horizonte caracterizado por presentar un
valor de velocidad de onda P calculado de 715m/s. En la
posición 15m a lo largo del perfil, muy cercano a la
perforación, el horizonte superficial presentó un espesor
de 1.8m (ver Figura 7(b)).
Figura 6. (a) Gráfico bilogarítmico de las observaciones o datos de campo,
respuesta del modelo o datos de campo sintéticos versus el espaciamiento 3.3. Resultado de prospección de onda S
medio (AB/2) entre los electrodos de corriente eléctrica, (b) modelo de terreno
obtenido posterior a la resolución del problema inverso (línea negra continua)
El proceso de apilamiento de los espectros de dispersión
con 6 modelos equivalentes (líneas continuas en tono gris), e (c) información normalizados obtenidos a partir de los puntos de
estratigráfica del sitio de interés. impactos desarrollados en las posiciones -1m y 31m
(Figura 8(a)), se presenta en la Figura 8(b). La inversión
3.2. Resultado de prospección de onda P 1D de la curva de dispersión extraída de la Figura 8(b) se
La Figura 7(a) muestra las curvas de tiempos de ilustra a través de la línea negra continua de la Figura 8(c)
llegada en función de la distancia para todos los 5 puntos la cual representa el menor valor de discrepancia
de impacto establecidos a lo largo del perfil, mientras que (0.0118) de los 2550 modelos calculados. Dicho modelo
la Figura 7(b) muestra el modelo de capas del terreno presenta 2 horizontes: el superficial de 3.15m de espesor
obtenido mediante la resolución del problema inverso. El y de 85m/s de velocidad de onda S calculada, seguido por
error RMS en la inversión de la matriz obtenido para esta un horizonte de mayor velocidad (177m/s). La
prueba fue de 1.82ms. profundidad de investigación (DOI) obtenida en esta
prueba fue de 5.3m, tal como se mencionó con
anterioridad.

RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022. 113
 Combinando el método de refracción sísmica y análisis multicanal de ondas superficiales para la detección de niveles freáticos someros en el noroccidente del Istmo de Panamá
Combining the seismic refraction method and multichannel analysis of surface waves for detection of shallow water-table in the northwestern Isthmus of Panama

Figura 9. (a) Modelos 1D de velocidad de onda P extraído del perfil de

refracción sísmica de la Figura 7(b) y de velocidad de onda S obtenido en la
Figura 8(b); (b) razón de Poisson computada a partir de los modelos 1D de
Figura 8. (a) Espectros de dispersión normalizados obtenidos en los puntos de onda P y S; (c) información estratigráfica del sitio y (d) contenido natural de
impacto establecidos en las posiciones -1.0m y 31.0m. (b) Espectro y curva de agua; n.f. corresponde al nivel freático.
dispersión obtenidos a partir del apilamiento de los dos puntos de impacto, y
(c) resultado de la inversión 1D de los datos de velocidad de fase vs frecuencia
extraídos de la Figura 8(a) empleando el algoritmo neighborhood Por su parte, el resultado obtenido en el análisis multi-
implementado por [24] en el programa Geopsy. DOI corresponde a la canal de ondas superficiales (1D) no presenta una
profundidad de investigación y n.f., el nivel freático.
correlación muy significativa entre el límite de los dos
horizontes identificados en esta prueba (3.15m de
Con los modelos 1D de velocidad de las ondas P y S
profundidad) y el nivel freático (1.80m); sin embargo, el
en la posición 15m, se procedió en calcular los valores de
primer horizonte de velocidad de onda S de 85m/s se
la razón de Poisson en función de la profundidad para un
asocia con el limo elástico de plasticidad media, y el
DOI de 5.3m. La Figura 9(a) presenta ambos modelos 1D
segundo de 177m/s, al limo de plasticidad baja. Para el
desplegados en un mismo gráfico.
coeficiente de Poisson (Figura 9(b)) se presentan 3
Para la prospección de onda P, el primer horizonte de
horizontes: con un valor de 0.31 y un espesor de 1.8m
velocidad calculada de 163m/s se encuentra asociado al
para el horizonte superficial; este es un típico valor que
limo elástico de plasticidad media identificado en la
de acuerdo con [26] y [27] corresponde a sedimentos no
perforación, mientras que el segundo horizonte con
consolidados y no saturados, el cual es el caso del limo
velocidad de 715m/s se asocia con un limo de plasticidad
elástico de plasticidad media identificada en la
baja. Como se mencionó antes, el nivel freático se
perforación.
correlaciona bien con el resultado de esta prueba; a partir
Entre 1.8m y 3.15m de profundidad se presenta un
de 1.80m de profundidad, la velocidad de onda P se
segundo horizonte con un valor del coeficiente de
mantiene constante hasta los 5.3m.
Poisson igual a 0.49, el cual es típico de los sedimentos
saturados ([28], [29], [30], [31]) y corresponde al limo de
plasticidad baja. El límite entre estos dos horizontes se
asocia con el nivel freático identificado en la perforación.
El tercer horizonte posee un valor de coeficiente de
Poisson cercano al identificado en el horizonte
precedente (0.47), por lo que se asocia al mismo tipo de
limo de baja plasticidad, pero con un menor nivel de
saturación.

4. Conclusiones
En la actualidad, los métodos sísmicos no han sido
extensamente utilizados en estudios hidrogeológicos en

114 RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022.
 Alexis Mojica | Osiris Cubilla | Eduardo Silva |José Ladrón-De Guevara | Alberto Ruíz | Carlos A. Ho | Jaime Toral

comparación con los métodos geoeléctricos y misma. J.L.G: redacción parcial y revisión del
electromagnéticos; como método de verificación de los manuscrito. A.R.: análisis del contexto geológico del
límites de saturación del suelo, el sondeo eléctrico mostró área de estudio. C.H.: configuración y análisis de los
una buena correlación entre el nivel freático reportado y datos de onda S. J.T.: estudio del contexto geológico y
el límite estratigráfico de los limos identificados en la tectónico de la zona de interés y revisión del artículo.
perforación. Por su parte, la combinación de perfiles de Todos los autores afirmamos que se leyó y aprobó la
velocidad de onda P y S a través del cálculo de la razón versión final de este artículo.
de Poisson, podría jugar un rol importante en la detección
de los niveles freáticos en los suelos sedimentarios, y no
REFERENCIAS
solamente a nivel de detección de los límites de los [1] S.S. Hubbard y Y. Rubin. “Introduction to
diferentes estratos sino también a la extensión del tamaño Hydrogeophysics”, en: Hydrogeophysics (Eds. Y. Rubin y
ya que los estratos saturados se caracterizan por presentar S.S. Hubbard), Holanda: Springer, pp. 3-22. 2005
valores cercanos al límite de 0.50. En el área de interés, [2] S.S. Hubbard y N. Linde. “Hydrogeophysics”. Treatise on
el análisis de refracción sísmica presentó una buena Water Science, vol. 2, pp. 401-434. 2011
[3] A. Mojica A. “Geoelectrical Sounding and Imaging over the
correlación con: el resultado del sondeo eléctrico, el nivel Central Zone of Panama”, en: Geophysics (Ed. A. Okiwelu),
freático y una diferencia aceptable de 30cm con el límite Londres: IntechOpen, pp. 136-147. 2018
estratigráfico entre ambos limos; sin embargo, el modelo [4] A. Mojica, I. Díaz, C.A. Ho, F. Ogden, R. Pinzón, J.
1D obtenido a través del análisis multicanal de ondas Fábrega, D. Vega y J. Hendrickx. “Study of seasonal rainfall
infiltration via time – lapse surface electrical resistivity
superficiales presentó discrepancias con el resultado del tomography. A case study of Gamboa area, Panama Canal
sondeo eléctrico, la lectura del nivel freático y con el Watershed”, Air, Soil and Water Research Journal, vol. 6,
límite estratigráfico. El cálculo del coeficiente de Poisson pp. 131-139. 2013
obtenido en el área de estudio permitió detectar el nivel [5] F. Ogden, A. Mojica, E. Kempema, J. Briceño y J. Regina.
freático a la profundidad reportada y por la magnitud “Diagnosing Hydrologic Flow Paths in Forest and Pasture
Land Uses within the Panama Canal Watershed Using
(0.49), y manteniéndose casi constante (con una leve Simulated Rainfall and Electrical Resistivity Tomography”,
diferencia de 0.02) hasta 3.15m y luego hasta el máximo en AGU Fall Meeting Suppl., San Francisco, California,
de profundidad alcanzado en este trabajo. 2014
[6] E. Kempema, A. Mojica, G. Litt, F. Ogden y A. Carey, A.
“Effect of forest age on rainwater infiltration in the lowlands
AGRADECIMIENTOS humid tropics”, en AGU Fall Meeting Suppl., San
Los autores desean agradecer a la Universidad Francisco, California, 2015
Tecnológica de Panamá y a la Secretaría Nacional de [7] M. Castrellón, J. Fábrega, L. Foglia, A. Mojica, A. Ruíz y
S. Saavedra. “Groundwater model of a fractured rock
Ciencia, Tecnología e Innovación (SENACYT) a través system as a tool for groundwater management: The Estibaná
del SNI, por todo el apoyo brindado para el desarrollo de Sub-Catchment, Azuero Peninsula, Panama”, en E-
este proyecto. proceedings of the 38th IAHR World Congress, Panama,
2019
[8] G. Litt, F. Ogden, A. Mojica, E. Kempema, C. Gardner, J.
CONFLICTO DE INTERESES Regina, J. Bruce, J. Harrison, J. Hendrickx, M. Bretfeld, Y.
Los autores declaran no tener algún conflicto de Cheng, R. Stallard, J. Hall y W. Lyons. “Land Cover Effects
interés. on Soil Infiltration Capacity Measured Using Plot Scale
Rainfall Simulation in Steep Tropical Lowlands of Central
Panama”. Hydrological Processes, vol. 34, no. 4, pp. 878-
CONTRIBUCIÓN Y APROBACIÓN DE LOS 897, 2020.
AUTORES [9] F. Rejiba, A. Mojica, M. Schmutz, A. Ruíz, M. Castrellón,
J. Ladrón De Guevara, S. Saavedra, J. Fábrega, C. Finco, C.
A.M.: medición de parámetros geofísicos en campo, Schamper, L.H. Cavalcante y M. Llubes. “Geophysical
análisis e interpretación de los datos geofísicos, characterization of near-surface formations in the La Villa
redacción del artículo. O.C.: medición de parámetros River catchment (Los Santos, Panama)”, en Conference
geofísicos en campo, análisis de las muestras de suelo Proceedings, NSG2021 1st Conference on
Hydrogeophysics, Bordeaux, Francia, 2021, pp. 1-5
obtenidas en la perforación . E.S.: coordinación de la [10] A. Mojica, B. Duarte, F. Vergara, M. Pinto, R. Pinzón, J.
perforación y el análisis de las muestras obtenidas en la Pérez, M. Castrellón, C.A. Ho y M. Gómez. “Time-lapse
electrical resistivity tomography for assessment of seasonal

RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022. 115
 Combinando el método de refracción sísmica y análisis multicanal de ondas superficiales para la detección de niveles freáticos someros en el noroccidente del Istmo de Panamá
Combining the seismic refraction method and multichannel analysis of surface waves for detection of shallow water-table in the northwestern Isthmus of Panama

moisture variations in a tropical regolith”. Hydrological [22] A. O’Neill. “Full-waveform reflectivity for modelling,
Research Letters, vol. 16, no. 1, pp. 18-24, 2022 inversion and appraisal of seismic surface wave dispersion
[11] S. Pasquet, L. Bodet, A. Dhemaied, A. Mouhri, Q. Vitale, in shallow site investigation”, PhD thesis, The University of
F. Rejiba, N. Flipo y R. Guérin. “Detecting different water Western Australia, 2003
table levels in a shallow aquifer with combined P-, surface [23] D. Zywicki y G. Rix. “Mitigation of near-field effects for
and SH-wave surveys: Insights from Vp/Vs or Poisson’s seismic Surface wave velocity estimation with cylindrical
ratios”. Journal of Applied Geophysics, vol. 113, pp. 38-50, beamformers”. Journal of Geotechnical and
2015 Geoenvironmental Engineering, vol. 131, no. 8, pp. 970-
[12] S. Pasquet, L. Bodet, L. Longuevergne, A. Dhemaied, C. 977, 2005
Camerlynck, F. Rejiba y R. Guérin. “2D characterization of [24] M. Sambridge. “Geophysical inversion with a
near-surface Vp/Vs: surface-wave dispersion inversion neighborhood algorithm-I. Searching a parameter space”.
versus refraction tomography”. Near Surface Geophysics, Geophysical Journal International, vol. 138, no. 2, pp. 479-
vol. 13, pp. 315-331, 2015 494, 1999
[13] S. Pasquet y L. Bodet. “SWIP: An integrated workflow for [25] M. Wathelet, D. Jongmans y M. Ohrnberger. “Surface-wave
surface-wave dispersion inversion and profiling”. inversion using a direct search algorithm and its application
Geophysics, vol. 82, no. 6, pp. WB47-WB61, 2017 to ambient vibration measurements”. Near Surface
[14] Atlas Nacional de la República de Panamá. Instituto Geophysics, vol. 2, no. 4, pp. 211-221, 2004
Geográfico Nacional Tommy Guardia. Colombia: [26] M. Wathelet M. “An improved neighborhood algorithm:
Impresiones Carpal, pp. 36-37, 2016. parameter conditions and dynamic scaling”. Geophysical
[15] M. Karous y T.K. Pernu. “Combined sounding-profiling Research Letters, vol. 35, no. 9, pp. L09301, 2008
resistivity measurements with the three-electrode arrays”. [27] K.H. Stokoe y R.D. Woods. “In situ shear wave velocity by
Geophysical Prospecting, vol. 33, no. 3, pp. 447-459, 1985 cross-hold method”. Journal of Soil Mechanics &
[16] M.C. Stadelhofen. Applications de la géophysique aux Foundations Div., vol. 98, pp. 443-460, 1972
recherches d'eau. Francia: Technique et Documentation- [28] D. Tiab y E.C. Donaldson. Petrophysics. Theory and
Lavoisier, 1991 practice of measuring reservoir rock and fluid transport
[17] Y.L. Ekinci y A. Demirci. “A Damped Least-Squares properties. Houston: Gulf Publishing, 1996
Inversion Program for the Interpretation of Schlumberger [29] H. Stuempel S. Kaehler, R. Meissner y B. Milkereit. “The
Sounding Curves”. Journal of Applied Sciences, vol. 8, no. use of seismic shear waves and compressional waves for
22, pp. 4070-4078, 2008 lithological problems of shallow sediments”. Geophysical
[18] D.W. Marquardt. “Generalized Inverses, Ridge Regression, Prospecting, vol. 32, no. 4, pp. 662-675., 1984
Biased Linear Estimation, and Nonlinear Estimation”. [30] R. Meissner, H. Stuempel y F. Theilen. “Shear wave studies
Technometrics, vol. 12, no. 3, pp. 519-612, 1970 in shallow sediments. Seismic shear waves”, en Handbook
[19] J.R. Inman. “Resistivity inversion with Ridge Regression”. of geophysical exploration. Applications, (Eds. K. Helbig y
Geophysics, vol. 40, no. 5, pp. 798-817, 1975 S. Treiter), London: Geophysical Press, pp. 225-253, 1985
[20] T. Iwasaki. “Extended time-term method for identifying [31] H.S. Salem. “Poisson's ratio and the porosity of surface soils
lateral structural variations from seismic refraction data”. and shallow sediments, determined from seismic
Earth Planets Space, vol. 54, pp. 663-677, 2002 compressional and shear wave velocities”. Géotechnique,
[21] A. Scheidegger y P.L. Willmore. “The use of a least square vol. 50, no. 4, pp. 461-463, 2000
method for the interpretation of data from seismic surveys”. [32] O. Uyanik. “The porosity of saturated shallow sediments
Geophysics, vol. 22, no. 1, pp. 9-21, 1957 from seismic compressional and shear wave velocities”.
Journal of Applied Geophysics, vol. 73, no. 1, pp. 16-24,
2011.

116 RIDTEC (ISSN 1680-8894 - ISSN-L 2219-6714) | Vol. 18, n.° 1, enero - junio 2022.