1. Resumen — todo el problema en una página
Un pozo en Cumaná puede fracasar por exactamente tres causas, y no hay una cuarta: no da agua, da agua salada, o da agua contaminada. Las tres están gobernadas por dos escalares del subsuelo: la conductividad hidráulica K (¿sale agua?) y la conductividad eléctrica del agua de poro σw (¿sirve?).
El núcleo
Los dos escalares se leen desde la superficie, sin perforar, en una sola magnitud: la resistividad eléctrica del terreno, ρb.
La ley de Archie los une: la resistividad del terreno es el producto de la geometría del poro (que es K) por la resistividad del fluido que lo llena (que es la salinidad).
Pero hay una trampa, y es el verdadero núcleo del problema: los dos modos de fracaso producen la misma señal. Una arcilla saturada y una arena llena de agua de mar tienen ambas resistividades bajísimas, de 1 a 10 Ω·m. Un mínimo de resistividad en el sondeo puede significar «aquí no hay acuífero» o «aquí hay un acuífero excelente lleno de agua de mar»: conclusiones opuestas, lectura idéntica.
Y en Cumaná esto es letal, porque el subsuelo tiene las dos cosas a la vez: la llanura es limo-arcillosa y hay cuña salina en el delta.
La solución, en una frase
La resistividad dice «aquí conduce».
La cargabilidad (polarización inducida) dice por qué conduce.
La arcilla polariza —tiene doble capa eléctrica—; el agua salada no —conduce por iones libres, es óhmica pura—. Ese par de observables resuelve el problema de Cumaná.
A eso se suma un segundo número, que decide dónde no perforar. Por Ghyben-Herzberg, cada metro de carga hidráulica dulce sostiene 40 metros de agua dulce bajo el nivel del mar. Es decir: si bombeas y bajas la carga un metro, la interfaz salina sube cuarenta. Esa amplificación de 40× es lo que convierte al puente Gómez Rubio —un sitio con estación de bombeo ociosa, tentador— en una trampa irreversible.
La consecuencia práctica, que cabe en una línea
No pedir «sondeos eléctricos verticales». Pedir «sondeos de resistividad
con polarización inducida, calibrados contra el pozo existente del INTI».
Mismo equipo, mismo despliegue, mismo día de trabajo. Y es la diferencia entre un estudio
que distingue arcilla de sal y uno que no.
2. El problema, formulado como física
| Modo de fracaso | Magnitud que lo gobierna | |
|---|---|---|
| A | No da agua (o da un hilo) | Conductividad hidráulica K → transmisividad T = K·b |
| B | Da agua salada | Contraste de densidad ρs vs ρf → posición de la interfaz |
| C | Da agua contaminada | Velocidad advectiva v = K·i / ne desde las cloacas |
Obsérvese qué tienen en común: en A, B y C aparece K, o aparece la salinidad del fluido de poro. No hay una cuarta variable escondida. Todo el problema se reduce a esos dos escalares, y de ahí que la prospección eléctrica no sea una elección arbitraria entre métodos: es el método que mide precisamente esas dos cosas.
Nota de clase — ¿por qué K y no la porosidad?
Es un error frecuente. Una arcilla puede tener porosidad del 50 % —más que una grava— y sin embargo no entregar agua: sus poros son diminutos y el agua queda retenida por tensión superficial y adsorción. Lo que importa no es cuánta agua hay, sino cuán fácilmente circula. Eso es K, y depende del cuadrado del tamaño de grano (Kozeny-Carman). Por eso K barre diez órdenes de magnitud: ~10−9 m/s en arcilla, ~10−2 m/s en grava limpia.1
3. La resistividad lo ve todo — ley de Archie
Para un medio granular saturado y limpio de arcilla, la resistividad del terreno obedece la relación empírica de Archie (1942):2
con a ≈ 1 y m ≈ 1,3 – 2,0 (el exponente de cementación crece con la
consolidación; en arena suelta, m ≈ 1,3 – 1,8).
Léase la estructura de la fórmula, porque ahí está todo: la resistividad del terreno es el producto de dos factores independientes:
- φ−m — la geometría del poro. Eso es K.
- ρw — la resistividad del agua que lo llena. Eso es la salinidad.
Un solo número observable contiene ambas incógnitas.
Los órdenes de magnitud que hay que tener en la cabeza
| Agua | σw | ρw = 10⁴/σw |
|---|---|---|
| Lluvia / muy dulce | ~50 µS/cm | ~200 Ω·m |
| Potable típica | 300 – 800 µS/cm | 12 – 35 Ω·m |
| Salobre | ~5.000 µS/cm | ~2 Ω·m |
| Mar | ~50.000 µS/cm | ~0,2 Ω·m |
El agua barre tres órdenes de magnitud. Por eso la salinidad domina la señal eléctrica: es, con diferencia, el factor más ruidoso y el más informativo.
| Material del subsuelo | ρb | |
|---|---|---|
| Grava / arena limpia, agua dulce | 300 – 2.000 Ω·m | ✅ el objetivo |
| Arena, agua potable | 100 – 400 Ω·m | ✅ |
| Limo | 10 – 50 Ω·m | — |
| Arcilla saturada | 1 – 10 Ω·m | ❌ no da agua |
| Arena con agua salada | 1 – 10 Ω·m | ❌ agua inservible |
Comprobémoslo con números
Arena, φ = 0,30, m = 1,8, a = 1 →
F = 0,30−1,8 ≈ 8,7
con agua de mar (ρw = 0,2 Ω·m) → ρb = 8,7 × 0,2 ≈ 1,7 Ω·m la misma arena · idéntica geometría de poro · factor 150× de diferencia
La misma arena, con la misma porosidad y la misma permeabilidad, cambia su resistividad en un factor de ~150 solo por el fluido que la llena. La medida es exquisitamente sensible a lo que queremos saber.
4. 🔴 La degeneración — aquí está el verdadero núcleo
Vuelva a mirar las dos filas rojas de la tabla anterior:
El problema
Arcilla saturada: 1 – 10 Ω·m
Arena con agua salada: 1 – 10 Ω·m
Los dos modos de fracaso producen la misma señal. La resistividad, por sí sola, no los distingue.
Esto no es un detalle de laboratorio: es una degeneración del observable. Un mínimo de ρb en el sondeo admite dos interpretaciones opuestas:
- «Aquí no hay acuífero» — porque es arcilla. Descarta el sitio.
- «Aquí hay un acuífero magnífico… lleno de agua de mar» — también descarta el sitio, pero por una razón completamente distinta, y con implicaciones distintas para el resto de la campaña.
Por qué en Cumaná esto es letal, y no académico
Porque el subsuelo de Cumaná tiene las dos cosas a la vez. El MOP (1969) documenta que la llanura donde está la ciudad es limo-arcillosa. Y el mismo proyecto documenta la cuña salina penetrando por el delta.
Una campaña de resistividad sin más, en este subsuelo, devolverá mínimos que nadie sabrá interpretar. Y el riesgo no es simétrico: confundir «arena salada» con «arcilla» hace que descartes una zona que quizá sea recuperable; confundir «arcilla» con «arena salada» te hace perforar donde no hay nada.
Nota de clase — degeneración, equivalencia y supresión
En la inversión de sondeos eléctricos existen tres ambigüedades clásicas, y conviene no confundirlas:
- Equivalencia: una capa delgada y muy conductora es indistinguible de una gruesa y menos conductora — solo se resuelve el producto conductancia = σ·h (o la resistencia transversal ρ·h). El espesor y la resistividad por separado no se determinan.
- Supresión: una capa intermedia delgada puede no aparecer en la curva de campo.
- Degeneración física (la nuestra): dos materiales físicamente distintos dan la misma resistividad. No es un problema de inversión: es un problema del observable. No se arregla con más datos de resistividad, ni con mejores algoritmos. Solo se arregla midiendo OTRA cosa.3
5. Cómo se rompe la degeneración — polarización inducida
La corriente eléctrica atraviesa el subsuelo por dos mecanismos físicamente distintos, y esa distinción es la llave:
| Mecanismo | Dónde ocurre | ¿Almacena carga? |
|---|---|---|
| Conducción electrolítica iones libres moviéndose en el fluido de poro |
agua salada en los poros | NO — es óhmica pura, sin memoria |
| Conducción de superficie doble capa eléctrica sobre las láminas de arcilla |
arcillas (y sulfuros metálicos) | SÍ — la doble capa se carga y se descarga |
La polarización inducida mide exactamente eso. Se inyecta corriente, se corta, y se observa la tensión residual que decae: la carga que el medio había almacenado.
| Material | ρb | Cargabilidad M | Diagnóstico |
|---|---|---|---|
| Arcilla saturada | 1 – 10 Ω·m | 10 – 50 mV/V | ALTA → doble capa → es arcilla |
| Arena con agua salada | 1 – 10 Ω·m | < 5 mV/V | BAJA → conducción iónica → es sal |
| Arena/grava, agua dulce | 300 – 2.000 Ω·m | < 5 mV/V | alta ρ + baja M → el objetivo |
La frase que resume el método
La resistividad dice «aquí conduce».
La cargabilidad dice por qué conduce.
Un par de observables independientes — (ρb, M) — separa inequívocamente las dos causas que la resistividad sola confunde.
Nota de clase — el origen físico de la polarización
Las láminas de arcilla tienen carga superficial negativa neta (por sustitución isomórfica en la red cristalina). Esa carga atrae una nube de cationes del fluido, formando la doble capa eléctrica (Stern + difusa). Cuando se aplica un campo, esos iones adsorbidos se desplazan pero no pueden escapar: se acumulan en un extremo. Al cortar la corriente, relajan por difusión, y esa relajación es la tensión residual que medimos.
El agua salada en un poro no tiene superficie que polarizar: sus iones circulan libres y, al cortar la corriente, se detienen sin memoria. Conduce mucho, pero no almacena nada. Ahí está toda la diferencia, y es una diferencia de mecanismo, no de grado.
Aclaratoria — esto NO es un sobrecosto
Un resistivímetro moderno mide ρ y M en la misma inyección de corriente: el mismo equipo, los mismos electrodos, el mismo despliegue, la misma jornada de campo. Pedir IP no encarece la campaña de forma significativa: cambia la especificación, no el presupuesto.
Pero si el pliego no lo pide, no lo miden. Y entonces el estudio nace ciego a la mitad del problema.
6. El anclaje — un pozo lo vuelve todo cuantitativo
La geofísica de superficie es, por naturaleza, relativa: entrega contrastes. Un pozo existente la vuelve absoluta.
En el pozo del INTI se puede medir ρw directamente, con un simple conductímetro sumergido. Y con ρw conocido, la ley de Archie deja de tener dos incógnitas y pasa a tener una:
Lo que realmente aporta el pozo del INTI
No es solo «agua que ya está perforada».
Es el punto que convierte toda la campaña geofísica de cualitativa en cuantitativa.
Sin él, un sondeo dice «aquí es más resistivo que allá».
Con él, dice «aquí hay arena con agua de tantos µS/cm».
⚠️ Corrección a una versión anterior de este análisis — y es importante
En una versión previa se afirmó que «con 10 o 15 aljibes se levanta el primer mapa piezométrico del acuífero de Cumaná». Eso era optimista y hay que retirarlo.
Revisada la memoria oral con frialdad, los puntos concretos, con dirección y probablemente accesibles son dos: la casa del Dr. Urosa y la casa de Jairo Díaz, ambas en San Francisco. Todo lo demás son sectores (el casco histórico, La Casimba, Bebedero) o sitios históricos que muy probablemente estén cegados. Sumando el pozo del INTI, tenemos a lo sumo dos o tres puntos de control — no catorce.
Y esto cambia la física del planteamiento, no solo la logística:
- No se puede construir un campo de carga h(x,y) con dos puntos. Sin campo de carga no hay gradiente ∇h, y sin gradiente no se conoce la dirección del flujo — es decir, no se sabe de dónde viene el agua ni hacia dónde viaja la contaminación.
- Por tanto, la carga de la prueba recae íntegramente sobre la geofísica. Y con un único punto de calibración, la polarización inducida deja de ser una mejora y pasa a ser IMPRESCINDIBLE: es el único modo de romper la degeneración sin una red de control densa.
- El campo piezométrico hay que construirlo, no encontrarlo: con piezómetros de pequeño diámetro perforados a propósito. Son baratos comparados con un pozo de producción, y sin ellos el modelo hidrogeológico no cierra.
7. El número 40 — por qué bombear cerca de la cuña la acelera
Equilibrio hidrostático entre una lente de agua dulce (ρf = 1.000 kg/m³) flotando sobre agua salada (ρs = 1.025 kg/m³):4
Cada metro de carga hidráulica dulce por encima del nivel del mar sostiene cuarenta metros de agua dulce por debajo. La lente flota, y flota con un lastre de 40 a 1.
La consecuencia
Si bombeas y bajas la carga 1 metro, la interfaz salina sube 40 metros.
Un abatimiento modesto produce un ascenso salino catastrófico. La amplificación es de 40×. Ese es el número que convierte al puente Gómez Rubio —con su estación de bombeo ociosa, sus tableros ya instalados y su conexión al alimentador de 36″— de sitio conveniente en trampa.
Y el daño no es reversible en escalas humanas. Una vez que el agua salada entra en la matriz porosa, difunde hacia los poros muertos y hacia las láminas de arcilla, desde donde vuelve a salir lentamente durante años. El lavado por advección tarda décadas, incluso si se detiene el bombeo.
Nota de clase — «upconing» y por qué el pozo profundo es peor
Localmente, bajo un pozo que bombea, la interfaz no sube de forma plana: forma un cono invertido que asciende hacia la rejilla (upconing). Existe un abatimiento crítico: superado, el cono se vuelve inestable y la sal irrumpe bruscamente en el pozo. No hay aviso gradual: el pozo pasa de dulce a salobre casi de golpe.
Corolario contraintuitivo: en la costa, un pozo MÁS PROFUNDO no es más seguro — es más peligroso, porque su rejilla está más cerca de la interfaz. La estrategia correcta en zona costera es la contraria: pozos SOMEROS, de bajo caudal y muchos, en vez de uno profundo y potente.5
Aclaratoria — los límites del modelo
Ghyben-Herzberg supone interfaz abrupta (en la realidad hay una zona de mezcla por dispersión), régimen estático y flujo esencialmente horizontal (hipótesis de Dupuit). Es una idealización. Pero su valor no está en la precisión del número: está en el orden de magnitud del acoplamiento. La relación es 40 a 1, no 2 a 1. Aunque el modelo se equivocara por un factor 2, la conclusión no cambia.
8. El paleocauce corta para los dos lados — Darcy
Ley de Darcy y velocidad real del agua en el poro:
Un paleocauce —un río enterrado, un cuerpo de arena y grava limpias sepultado bajo la
ciudad— con K ≈ 10⁻³ m/s, gradiente i ≈ 10⁻³ y ne ≈ 0,25:
y en grava limpia (K ≈ 10⁻² m/s) → ≈ 1,3 km/año
El activo y el pasivo son la misma cosa
El mismo cuerpo permeable que te da caudal te trae la contaminación de las cloacas desde un kilómetro de distancia, en un año.
Alta K es simultáneamente lo que buscas y lo que te mata. El gradiente hidráulico decide el signo.
De ahí que el paleocauce sea un objetivo condicionado: excelente si está aguas arriba de la fuente de contaminación; pésimo si está aguas abajo. Y de ahí, también, que el piedemonte gane: no solo tiene K alta (gravas coluviales) — está aguas arriba de todo: de las cloacas, de la cuña salina, de la ciudad entera.
Nota de clase — el «perímetro de protección» sale de esta fórmula
Con v y un tiempo de tránsito objetivo (típicamente 50 días, el tiempo que tardan en morir los patógenos bacterianos en el subsuelo), se define el radio mínimo que debe quedar libre de fuentes de contaminación alrededor del pozo:
r ≈ v · t → con v = 0,35 m/día y t = 50 d
→ r ≈ 17,5 m (arena) · pero en grava limpia (v = 3,5 m/día) → r ≈ 175 m
El perímetro de protección no es un número reglamentario arbitrario: sale de la permeabilidad del terreno. Cuanto mejor es el acuífero, más grande debe ser el perímetro que hay que proteger.
9. La campaña, ordenada por información obtenida por bolívar
| Fase | Medida | Observable que entrega | Costo |
|---|---|---|---|
| 0a | Pozo del INTI: σw (conductímetro), Fe, Mn, Cl⁻, caudal, nivel dinámico | ρw → el ancla de Archie. Convierte toda la geofísica posterior en cuantitativa. | frascos y una tarde |
| 0b | Los 2 aljibes vivos (casa del Dr. Urosa, casa de Jairo Díaz): nivel de agua, σw, profundidad | Segundo y tercer punto de control. Insuficientes para un campo de carga — ver la corrección del §6. | cinta, sonda y frascos |
| 1 | Perfil VERTICAL de salinidad del Manzanares (por estratos, con CTD) | La condición de contorno del problema. La cuña es corriente de fondo: un muestreo de superficie es ciego por construcción. | una lancha |
| 2 | Resistividad + POLARIZACIÓN INDUCIDA (SEV y/o tomografía) en Jagüey de Luna y sobre los paleocauces | El par (ρb, M) → rompe la degeneración arcilla/sal → mapa de K y de salinidad | campaña geofísica |
| 3 | Piezómetros de pequeño diámetro (3–5 puntos) | El campo de carga h(x,y) → ∇h → dirección de flujo. Hay que construirlo: los aljibes no bastan. | perforación menor |
| 4 | Pozo exploratorio + ensayo de bombeo | T y S — los únicos parámetros que predicen el caudal sostenible | perforación |
El ensayo de bombeo — Cooper-Jacob
Aproximación logarítmica de la solución de Theis, válida para tiempos largos:
Representando s frente a log t se obtiene una recta. De su pendiente sale T; de su intercepto sale S.
La consecuencia que hay que decir en voz alta
Antes del ensayo de bombeo, NADIE puede decir cuántos litros por segundo dará el pozo. Cualquier cifra anterior es una invención. Si en una reunión alguien promete un caudal sin haber hecho el ensayo, está adivinando — y conviene decirlo.
10. Consecuencia práctica — qué pedir en el pliego
La línea que hay que escribir
❌ NO pedir: «sondeos eléctricos verticales».
✅ SÍ pedir: «sondeos de resistividad con polarización inducida (IP), calibrados contra el pozo existente del INTI, con medición in situ de la conductividad del agua de poro».
Es la misma campaña, el mismo equipo, el mismo día de trabajo. Y es la diferencia entre un estudio que distingue arcilla de sal y uno que no.
En un subsuelo que es limo-arcilloso y salino a la vez —y el de Cumaná lo es, está documentado desde 1969— un sondeo de resistividad sin IP puede devolver un mínimo y nadie sabrá cuál de las dos cosas es. El informe llegará, se pagará, y no habrá resuelto la pregunta.
«Yo escucho el susurro que viene de las entrañas del suelo, del corazón de la tierra, así encontrar el agua: es un tiro fijo.»
El tiro fijo tiene ecuación.
Es Archie, es Darcy, es Ghyben-Herzberg.
Y el susurro se mide en ohm·metro.
Notas y aclaratorias
- Kozeny-Carman.
K ∝ (φ³ / (1−φ)²) · d². La dependencia con el cuadrado del diámetro de grano es lo que produce el barrido de diez órdenes de magnitud. Una arcilla con más porosidad que una grava puede tener una K un millón de veces menor, simplemente porque sus granos son mil veces más pequeños. Porosidad y permeabilidad no son la misma cosa, y confundirlas es el error clásico. - Los límites de Archie. La ley es empírica y se dedujo para areniscas limpias. En presencia de arcilla falla, porque aparece la conducción de superficie: la arcilla conduce además por su doble capa, no solo por el fluido. Los modelos de Waxman-Smits o de doble agua corrigen esto añadiendo un término de conductividad superficial. Y aquí hay una ironía útil: el hecho mismo de que Archie falle en presencia de arcilla es lo que hace posible detectarla — porque esa conducción de superficie es, precisamente, la que produce la polarización inducida. El defecto del modelo es la señal que buscamos.
- Degeneración ≠ mal condicionamiento. Vale la pena insistir en la distinción, porque tiene consecuencias prácticas. La equivalencia y la supresión son problemas del inverso: se mitigan con mejor cobertura, más aperturas, regularización o información a priori. La degeneración física es un problema del observable: dos estados del mundo distintos producen exactamente el mismo dato. Ningún algoritmo, ninguna cantidad de datos y ninguna regularización pueden distinguirlos. Solo se rompe midiendo una magnitud física diferente. Es la misma lógica por la que en física de partículas se añade un detector de otro tipo cuando dos canales dan la misma firma.
- Ghyben-Herzberg, 1888-1901. Descubierto independientemente por el
neerlandés Badon Ghyben y el alemán Herzberg, ambos observando lentes de agua dulce en dunas
costeras. El resultado se obtiene igualando presiones a la profundidad de la interfaz:
ρf·g·(h+z) = ρs·g·z→z = ρf·h/(ρs−ρf). Es hidrostática elemental — y sin embargo es el resultado que más pozos costeros ha salvado y arruinado en la historia. - Sobre el upconing. El criterio de estabilidad (Dagan-Bear) sitúa el ascenso crítico en torno a un tercio de la distancia inicial entre la rejilla y la interfaz. Superado ese umbral, el cono se vuelve inestable y la irrupción salina es abrupta. La consecuencia de diseño es contraintuitiva y conviene retenerla: en zona costera, muchos pozos someros de bajo caudal son muy superiores a un pozo profundo y potente, aunque el caudal total sea el mismo. La geometría del bombeo importa tanto como su magnitud.
- Fuentes de los datos geológicos citados. MOP, Proyecto del Aliviadero y Control de Inundaciones del río Manzanares (Ing. Víctor Silva Bermúdez, 1969) — litología limo-arcillosa, cuña salina, paleocauces. BID/IDOM/IH-Cantabria (2017), Informe 7.1 — nivel freático, licuefacción, piedemonte. Martínez, Senior & Márquez (2005), Bol. Inst. Oceanogr. Venezuela 44(2) — penetración de la cuña salina. Detalle completo en el expediente.