Las cotas PDVSA por nivelación geodésica convencional 2004, 2007 y 2012 [PDVSA, 2013], [Acuña, 2013] de 2089 BMs de control de subsidencia en la COLM-Edo. Zulia, y aquellas determinadas recientemente en esos sitios con altimetría satelital SWOT, época 2026.1 [Acuña, 2026], fueron combinadas óptimamente por GGenLUZ usando -fast least-squares prediction FLSP- para generar un modelo regional de alta-resolución 450x450m que representa la velocidad media vertical de la deformación del terreno por subsidencia en el periodo 2004-2026 y sus errores estimados, ver Figuras 1 y 2.
Figura 1. Velocidades medias verticales de deformación del terreno por subsidencia (cm/año) en BMs PDVSA según MDTSubCOLM0426v1.0 [Acuña, 2026]: Modelo GGenLUZ de Deformación de Terreno por Subsidencia en la COLM-Edo. Zulia, periodo 2004-2026, versión 1.0, de resolución 450x450m, con filtraje gaussiano de 900 m 'full-width'.
En términos de magnitudes de velocidades medias verticales de deformación, MDTSubCOLM0426v1.0 muestra tasas de subsidencia extremas entre -12 cm/año y +1 cm/año, con valor medio en la región de -1.2 cm/año. Para el periodo 2004-2026, el modelo señala la zona de mayor deformación en Lagunillas y otras significantes de deformación moderada en Tía Juana, Bachaquero y San Lorenzo (ver Figura 1).
Respecto a los errores de las velocidades asociados a MDTSubCOLM0426v1.0, estos se acumulan con valores elevados en la zona de mayor deformación de Lagunillas, donde se presentan en forma de incertidumbres de hasta ±3.7 cm/año para el referido periodo. En general, como valor medio para la región, las velocidades muestran errores estándar de ±0.5 cm/año (ver Figura 2).
MDTSubCOLM0426v1.0 y su software de aplicación MDTSubCOLM0426v1.0_intp.yab/.exe [Acuña, 2026] permitirán estimar, o más precisamente, propagar en el tiempo, con calidad de pocos centímetros y en la zona de validez del modelo (ver Figura 1), la cota de cualquier BM o estación de control geodésico a una época especifica de interés en el periodo 2004-2026 (o próximo a éste), conocida su posición ITRF/GRS80 y cota en otra época, previa o posterior.
Sobre el uso de MDTSubCOLM0426v1.0:
El modelo de deformación es aplicable a través de su software de interpolación MDTSubCOLM0426v1.0_intp.yab/.exe [Acuña, 2026], ver Figura 3. Éste posibilita, de nuevo, para BMs o estaciones de control geodésico en la zona del modelo, con cotas niveladas o alturas elipsoidales GNSS, respectivamente, propagar en el tiempo 'corrigiendo por subsidencia' esos valores de elevación, conocida la posición ITRF/GRS80 de tales sitios. A continuación se presenta un ejemplo típico de la aplicación de MDTSubCOLM0426v1.0 con MDTSubCOLM0426v1.0_intp.yab/.exe.
Figura 3. Software GGenLUZ MDTSubCOLM0426v1.0_intp.yab/.exe [Acuña, 2026].
El ejemplo refiere al BM PDVSA 212A de control de subsidencia (también estación GPS/GNSS pasiva de medición episódica bienal), ubicado en Lagunillas, muy cerca de la zona de máxima deformación en la COLM por tal fenómeno, ver Figura 4.
El BM 212A fue medido con GPS de punto-preciso (PPP) durante varias sesiones en la campaña de subsidencia PDVSA de 2007; sus coordenadas resultantes refieren al ITRF2020/GRS80 en la época 2007.9 presentando calidad centimétrica (mejor a ±1-cm) en las 3 componentes de la posición [Acuña, 2025].
A partir de su altura elipsoidal observada (h=-18.7563 m), la cota del BM para la época de medición 2007.9 (H=-4.2508 msnl) fue determinada con nivelación-GPS y el geoide nacional VGM25v1.0 (N=-14.5055 m) desde los BMs profundos PDVSA 9201, 9202 y 9203. Ahora bien, aquí se emplea MDTSubCOLM0426v1.0_intp.yab/.exe para estimar su valor (cota) en la época actual, i.e., 2026.1 (aprox.).
El resultado de la determinación se muestra en detalle en las Figuras 4 y 5. Considerando un periodo de transformación de 18.2 años, y una velocidad media vertical anual de deformación por subsidencia de -1.76±0.5 cm/año (dada por el modelo MDTSubCOLM0426v1.0), la cota del BM fue propagada corrigiéndola por subsidencia a la época de interés 2026.1. Su valor estimado final, H=-4.5715 msnl, obtenido con un error de ±0.0919 m, refleja el impacto de -32.07 cm de hundimiento por subsidencia del sitio, entre las épocas 2007.9 y 2026.1. El valor propagado de la cota del BM 212A con MDTSubCOLM0426v1.0_intp.yab/.exe difiere sólo -5.99 cm del obtenido con altimetría satelital SWOT y -1.68 cm respecto al valor sintético predicho por regresión lineal según cuadrados mínimos utilizando el registro histórico de la cota del BM. Tales diferencias permanecen dentro del rango de error de la estimación con el software empleado. Esto es un indicador veraz del muy satisfactorio nivel de confiabilidad de la herramienta computacional objeto de esta nota técnica.
Figura 4. Aplicación de MDTSubCOLM0426v1.0_intp.yab/.exe en la estación GPS/BM PDVSA 212A para propagar su cota y corregirla por subsidencia desde la época de medición 2007.9 a la actualidad (2026.1 aprox.).
Figura 5. Resultados completos (archivo ASCII .out) de MDTSubCOLM0426v1.0_intp.yab/.exe para la estación GPS/BM PDVSA 212A [Acuña, 2026].
Una versión DEMO de MDTSubCOLM0426v1.0_intp.yab/.exe está disponible en el siguiente enlace:
Más información en:
Acuña, G. (2026): MDTSubCOLM0426v1.0 - modelo GGenLUZ de deformación de terreno por subsidencia 2004-2026 para la COLM-Edo. Zulia. Tópicos de Geodesia Geométrica. Mayo 10-16, 2026. Cátedra Geodesia Geométrica (GGenLUZ). Dpto. de Geodesia Superior. Esc. de Ingeniería Geodésica. Fac. de Ingeniería. Universidad del Zulia (LUZ). Maracaibo, Venezuela.
No hay comentarios:
Publicar un comentario