UAI<strong>Sustentabilidad</strong>Almacenamiento geológico <strong>de</strong> CO2Como p<strong>ar</strong>te inicial <strong>de</strong> aplicación <strong>de</strong> la metodología (Figura1), tras la recopilación, estudio y análisis <strong>de</strong> los datosgeológicos disponibles, y la aplicación <strong>de</strong> los criterios <strong>de</strong>selección <strong>de</strong> emplazamientos en formaciones geológicaspotencialmente favorables, a escala regional, se seleccionóel siguiente p<strong>ar</strong> almacén‐sello [2], al que se le aplic<strong>ar</strong>á elresto <strong>de</strong> la metodología p<strong>ar</strong>a calcul<strong>ar</strong> su capacidad teórica<strong>de</strong> almacenamiento:‐La formación confinante está constituida por losmateriales terci<strong>ar</strong>ios, que en la zona <strong>de</strong> estudio forman unpotente <strong>de</strong>pósito discordante sobre el Paleozoico, y sobreel Cretácico superior hacia el este. Presenta faciesanhidrítico‐<strong>ar</strong>cillosas en el centro y clásticas en el oeste.Encima se disponen los <strong>de</strong>pósitos continentales, <strong>ar</strong>enososy <strong>ar</strong>cillosos fluviales.Mo<strong>de</strong>los digitales <strong>de</strong>l terreno (mdt) einterpolación geoestadísticaTodo el proceso <strong>de</strong> construcción <strong>de</strong> los MDT <strong>de</strong>scansa, enprimer lug<strong>ar</strong>, sobre la calidad <strong>de</strong> los datos iniciales, lo queconstituye una primera fuente <strong>de</strong> incertidumbres. Perojunto a lo anterior, el manejo <strong>de</strong> una importante cantidad<strong>de</strong> información mediante el uso <strong>de</strong> SIG pue<strong>de</strong> suponer, asu vez, la introducción <strong>de</strong> otra serie <strong>de</strong> incertidumbres[19], que es neces<strong>ar</strong>io conocer p<strong>ar</strong>a que puedan sercontroladas.El uso <strong>de</strong> un SIG supone la necesidad <strong>de</strong> abord<strong>ar</strong> yconsi<strong>de</strong>r<strong>ar</strong> aspectos relacionados con los distintos sistemas<strong>de</strong> referencia y su posible conversión <strong>de</strong> unos a otros, tipos<strong>de</strong> coor<strong>de</strong>nadas a emple<strong>ar</strong>, c<strong>ar</strong>acterísticas básicas <strong>de</strong> lasprincipales proyecciones c<strong>ar</strong>tográficas y <strong>de</strong> lasconsecuencias <strong>de</strong>l uso <strong>de</strong> unas u otras, estimación <strong>de</strong> las<strong>de</strong>formaciones <strong>de</strong> las proyecciones, escalas <strong>de</strong> lasrepresentaciones, el problema <strong>de</strong> las profundida<strong>de</strong>s entrelas distintas capas, etc.Dado que el presente trabajo se <strong>de</strong>s<strong>ar</strong>rolla bajo la premisa<strong>de</strong> la necesidad <strong>de</strong> estim<strong>ar</strong> las incertidumbres, seconsi<strong>de</strong>ra imprescindible el conocimiento y control <strong>de</strong> susfuentes, entre las que se encuentran las <strong>de</strong>rivadas <strong>de</strong>l uso<strong>de</strong> un SIG. Dichas incertidumbres v<strong>ar</strong>ían en función <strong>de</strong> lanaturaleza <strong>de</strong> los datos que se estén manejando. Los datostopográficos se obtuvieron directamente en formatodigital, pero fueneces<strong>ar</strong>io introducir los datos geológicos en el SIG,mediante digitalización en pantalla, en el caso <strong>de</strong> losmapas <strong>de</strong> síntesis; y por introducción directa p<strong>ar</strong>a elresto <strong>de</strong> los datos (provenientes <strong>de</strong> son<strong>de</strong>os). Dado quela introducción <strong>de</strong> errores en la digitalización esinevitable (líneas que no cierran, líneas que sobrepasanel punto <strong>de</strong> unión, formación <strong>de</strong> polígonos ficticios, etc.),se hace imprescindible una etapa <strong>de</strong> edición p<strong>ar</strong>a la<strong>de</strong>tección y control <strong>de</strong> los mismos. Tras la digitalización,se proce<strong>de</strong> a su georreferenciación, lo que es esencialp<strong>ar</strong>a asegur<strong>ar</strong> la vali<strong>de</strong>z <strong>de</strong> la localización <strong>de</strong> los objetosen una base <strong>de</strong> datos espaciales. P<strong>ar</strong>a este trabajo, todoslos mapas, incluidos los proce<strong>de</strong>ntes <strong>de</strong> la NGDC,emplean como dátum geodésico el ED50 y como sistema<strong>de</strong> coor<strong>de</strong>nadas el UTM huso 30, con el fin <strong>de</strong> mantenerla proyección original <strong>de</strong> los mapas con informacióngeológica. A p<strong>ar</strong>tir <strong>de</strong> este punto y apoyados en lageoestadística (método Kriging) p<strong>ar</strong>a la realización <strong>de</strong> lasinferencias estadísticas, se obtienen los distintos MDT,entendiendo éstos como una estructura numérica <strong>de</strong>datos que representa la distribución espacial <strong>de</strong> unav<strong>ar</strong>iable cuantitativa y continua [20]. Los MDT permitenla representación, visualización y el análisis <strong>de</strong> los datos yconstituirán la base que posibilit<strong>ar</strong>á el cálculo final <strong>de</strong> lacapacidad <strong>de</strong> almacenamiento.En cuanto a la información que ha sido utilizada, losdatos <strong>de</strong>l mo<strong>de</strong>lo topográfico <strong>de</strong>l área <strong>de</strong> estudio (verFigura 2 y Figura 4) son <strong>de</strong> libre acceso y pertenecen alNGDC (National Geophysical Data Center), organismo <strong>de</strong>la NOAA (National Oceanic & AtmosphericAdministration), perteneciente a su vez al Dep<strong>ar</strong>tment ofCommerce <strong>de</strong> los Estados Unidos <strong>de</strong> América [21]. Elresto <strong>de</strong> la información, <strong>de</strong> naturaleza esencialmentegeológica, se encuentra en formato papel, dada suantigüedad. Ello obligó a realiz<strong>ar</strong> una serie <strong>de</strong>operaciones sobre esta información, hasta el cálculo final<strong>de</strong> los MDT y que a continuación se citan (ver Figura 4):‐ Digitalización <strong>de</strong> los datos.‐ Georreferenciación y proyección sobre los distintossistemas <strong>de</strong> referencia geodésicos.Transformación/conversión <strong>de</strong> los sistemas <strong>de</strong> referenciac<strong>ar</strong>tográficos (SRC) (transformación Helmert 3D (7p<strong>ar</strong>ámetros).
Ingeniería y gerenciamiento ambientalFigura 4: Etapas <strong>de</strong> transformación <strong>de</strong> datos p<strong>ar</strong>a la obtención <strong>de</strong> los MDT (digitalización, georreferenciaciónen la proyección seleccionada, extracción <strong>de</strong> los datos limitados al área <strong>de</strong> trabajo, e interpolacionesestadísticas), aplicado al mapa <strong>de</strong> profundida<strong>de</strong>s <strong>de</strong>l techo <strong>de</strong>l Cretácico superior c<strong>ar</strong>bonatado.