Esta es la versión para lectura y consulta. El portal completo, con el recorrido narrado, las visualizaciones y el glosario interactivo, está en ecprym.org.
Este capítulo es una guía didáctica de ingeniería de sistemas. Expone las siete decisiones que dan forma al clúster ECPRYM y las deriva de primeros principios. Cada decisión se presenta con la misma estructura: el reto de partida, la razón por la que la solución convencional no lo resuelve, el fundamento físico de la solución adoptada y la ecuación o el parámetro que la cuantifica, con su respaldo normativo y bibliográfico.
El océano es la fuente del agua, del viento más denso y constante y de los nutrientes del afloramiento de Benguela. El clima costero del desierto del Namib, sin embargo, degrada casi todo lo que no necesita estar junto al mar: los equipos se corroen, la radiación de la mañana llega atenuada y las sales higroscópicas no secan.
La práctica habitual concentra todo el complejo —plantas, campo solar, invernaderos— junto al mar para acortar las conducciones de agua. Esa elección minimiza una única partida, el transporte de agua, y traslada su coste a todas las demás: estructuras y electrónica en la categoría de corrosividad más severa, irradiación matinal imprevisible y superficies de evaporación más extensas.
El aire litoral circula sobre el agua de afloramiento, a unos 13 °C de media, con una humedad relativa superior al 85 % y un punto de rocío próximo a 14,5 °C. En contacto con el mar se enfría por debajo de ese punto y condensa: es la niebla de advección de Benguela, que la brisa introduce en la franja costera cada mañana. Sus consecuencias son tres:
Más allá del kilómetro 25, el aire diurno tiene una humedad inferior al 35 %, la niebla no alcanza la meseta y la categoría de corrosividad baja a C2–C3. La irradiación global horizontal supera los 2.200 kWh/m²·año sin atenuación matinal, y la evaporación potencial supera los 3.200 mm/año. El parque fotovoltaico ocupa por eso la franja de los kilómetros 25 a 40: el suelo desértico seco devuelve un albedo de 0,32 que la cara posterior de los módulos bifaciales convierte en producción, y el conjunto rinde un factor de capacidad del 30,6 % con el perfil físico del emplazamiento. Esos dos números son los que sostienen los 200,92 MW del parque en la matriz canónica y en el gemelo digital.
La contrapartida es el transporte: el agua desalada debe subir al interior por el acueducto DN1000. La toma y la impulsión consumen 0,84 kWh/m³ de los 3,00 kWh/m³ del consumo global. Ese coste se paga en la ventana solar, cuando sobra energía, y se compensa con la reducción de la superficie de evaporación, la vida útil de los equipos del interior y la predictibilidad del recurso solar. El aire más cálido del interior eleva la temperatura de célula y resta algunos puntos al rendimiento fotovoltaico frente a la costa; el balance es favorable al interior por corrosividad, predictibilidad y albedo del suelo seco para los módulos bifaciales, no por el rendimiento de célula.
T_d = \frac{243{,}04\,\gamma}{17{,}625-\gamma},\quad \gamma = \ln\frac{HR}{100} + \frac{17{,}625\,T}{243{,}04+T}\qquad A_{camas} = \frac{V_{evap}}{E_{neta}}
Con una tasa neta de 3.200 mm/año en lugar de 2.500 mm/año, la superficie de camas para el mismo caudal se reduce en torno al 22 %. Respaldo: Olivier (1995), Lancaster et al. (1984), Koloutsou-Vakakis y Rood (1994), Greenspan (1977), ISO 9223:2012 e ISO 12944-2:2017.
Una red con más de 290 MW de generación por convertidores —fotovoltaica y eólica— debe mantener la frecuencia y la tensión cuando un generador se desconecta de golpe. Sin masa rotante, la frecuencia cae en milisegundos.
La respuesta habitual es instalar acumuladores electroquímicos (BESS) o un bombeo hidroeléctrico reversible (PHES). Los acumuladores de litio pierden capacidad con los ciclos y con la temperatura, que en el interior del Namib supera los 40 °C en verano, exigen reposición de celdas a lo largo de la vida del proyecto y dependen de cadenas de suministro de materiales críticos. Su respuesta a una perturbación, además, es sintética: llega después de un retardo de medida y control. El bombeo reversible exige elevar agua para devolverla después, y con ella un ciclo de pérdidas que el clúster no necesita pagar: el agua que mueve tiene ya un destino propio —regar y abastecer— y desciende una sola vez, de modo que lo que procede es recuperar energía de ese descenso y no fabricar un ciclo de ida y vuelta.
La energía cinética de una masa rotante responde a una perturbación en el instante t = 0, por la propia física del acoplamiento electromagnético, sin medida ni decisión. El clúster dispone de 494 MW·s de energía cinética real: 480 MW·s de dos condensadores síncronos de 40 MVA con volante físico y 14 MW·s de los turbogeneradores de vapor de recuperación (HRSG). Con una potencia base de 107,9 MVA, la constante de inercia es H = 4,58 s.
La inercia sólo es suficiente si la perturbación es acotada. De ahí la regla topológica R1: el parque fotovoltaico de 200,92 MW se divide en seis circuitos de evacuación independientes, de modo que la mayor pérdida creíble al mediodía es de 18,56 MW. Con cuatro circuitos sería de 27,83 MW y el RoCoF superaría el límite de 1,0 Hz/s. Esa pérdida corresponde al despacho de 111,4 MW que se aplica en isla: con la irradiancia física del mediodía cada circuito llegaría a 33,49 MW, y por eso, cuando el enlace de 132 kV está abierto, el despacho fotovoltaico se limita a ese valor (Apéndice 19).
La regla R2 completa el diseño por el lado de la energía. La inercia resuelve los segundos, pero no las horas: en una noche en calma no hay viento ni sol que sostengan la carga firme de procesos continuos del clúster, 28,14 MW. Sin acumuladores ni bombeo reversible, esas horas exigen una máquina despachable, y el clúster la tiene en el grupo de biometano de 11,0 MW alimentado por su propio digestor, que aporta 4,81 GWh al año —unas 437 horas equivalentes a plena carga— y suma inercia por ser síncrono. Con él, y con la parada nocturna autorizada de la electrólisis, el modelo horario no importa energía en las noches ordinarias. A ese grupo se suma la central de pie de presa del sistema nodal de agua, que durante las catorce horas de descarga nocturna entrega 1,137 MW firmes y 15,92 MWh, con inercia propia por ser su alternador síncrono. Una calma prolongada —catorce horas con la eólica al 3 %— deja un hueco de 22,48 MW que ninguna máquina de 11 MW cubre: con el enlace cerrado se importa y en isla el clúster baja al nivel 3 de supervivencia; para esa noche el grupo tiene el combustible asegurado desde la Fase 1 con 28 t de biometano presurizado a 250 bar (unos 39.000 Nm³), el consumo exacto de una calma de catorce horas a plena carga, y es además de doble combustible: desde el Año 20 se le suma una reserva de 69,62 t de metanol de la planta de química verde. El combustible deja de ser la restricción; la potencia de la máquina sigue siéndolo. El enlace de 132 kV queda para lo que es: evacuación comercial del excedente y anclaje de tensión con la red regional, con más de 10 MW de margen rodante permanente.
El parque de 291,54 MW genera 898,50 GWh al año con los factores de capacidad físicos del emplazamiento: 30,6 % en la fotovoltaica del interior y 45,4 % en la eólica del litoral. De esa generación, 631,88 GWh —el 70,3 %— alimentan los procesos del clúster, y el resto es excedente. Lo que decide el destino del excedente es la capacidad de evacuación: el enlace de doble circuito admite hasta 39,52 MW de exportación, el mayor valor cuya pérdida simultánea sostiene la inercia del clúster (Apéndice 19), y con él se exportan 113,26 GWh al año y se vierten 153,37 por falta de salida; con 120 MW de evacuación se exportarían 265,22 GWh, el 29,5 % de la generación, y el vertido se anularía, pero esa ampliación exige circuitos adicionales independientes y su propio estudio de estabilidad. Ése es el sentido físico de la regla 70/30: un reparto anual entre consumo propio y red nacional, no una cuota horaria plana. La energía que sale por el doble circuito de 132 kV entra en la red troncal nacional. Desde ella, la distribuidora regional abastece a la población de Lüderitz. Ésa es la frontera entre la microrred industrial, que es privada y alimenta procesos, y el consumo municipal, que es un servicio público: el clúster entrega energía a la red y agua potable en alta, y no distribuye electricidad a los hogares. De día, con la fotovoltaica a pleno rendimiento, la exportación instantánea supera con holgura esa media; de noche, cuando el viento sólo cubre la carga firme, no hay excedente que exportar.
\frac{2H}{f_0}\frac{df}{dt} = \frac{P_{gen}-P_{dem}}{S_{base}}\qquad \mathrm{RoCoF} = \frac{f_0\,\Delta P}{2\,E_{cin}} = \frac{50 \times 18{,}56}{2 \times 494} = 0{,}94\ \mathrm{Hz/s}
El nadir de frecuencia ante cualquier contingencia N-1 del día de operación de diseño se mantiene por encima de 49,52 Hz, con el umbral de deslastre en 49,0 Hz. Respaldo: Kundur (1994); IEEE Std 2800-2022, que establece los requisitos de desempeño de los recursos basados en convertidores conectados a redes de transporte.
La generación solar es diurna y la eólica, variable; la demanda de agua es continua y tiene una punta estacional en el verano austral, cuando se riega el Tecnosol. Hay que casar tres ritmos distintos sin verter energía renovable y sin racionar agua.
La solución estándar almacena kilovatios hora eléctricos para mantener la desalación a caudal constante las 24 horas. Convierte dos veces la energía —a química y de vuelta a eléctrica—, con pérdidas en cada conversión, para alimentar una carga que no necesita ser constante.
El agua ya desalada e impulsada contiene la energía que costó producirla y elevarla. Si la planta produce cuando sobra energía y el agua se almacena, la demanda eléctrica se traslada en el tiempo sin reconversión alguna. La planta de ósmosis inversa y su impulsión operan como una batería virtual: hasta 12 MW en la ventana solar de 10:00 a 20:00 h y, fuera de ella, un único módulo de los diez a plena carga (1,2 MW), para no someter las membranas a ciclos de presión. El 84,7 % de la energía de bombeo se consume dentro de la ventana, y la producción diaria de agua se mantiene íntegra.
El embalse nodal S04 no es un depósito de regulación diaria: se dimensiona por la curva de masas del balance horario de un año completo —el acumulado de producción menos consumo—, que es la única forma de absorber el desfase estacional entre el recurso y la punta de riego. Su volumen elimina las 479 horas anuales de déficit que existirían con una regulación sólo diaria.
El coste relativo ilustra la decisión: el embalse cuesta 171,88 MUSD para 9.428 MWh de demanda desplazable, del orden de 18 USD por kWh, una fracción del coste por kWh de un acumulador electroquímico de gran escala. Además, el embalse sería necesario de todos modos para la regulación del agua: su función eléctrica es una economía de alcance.
E_{bv} = V_{S04}\cdot SEC = 3.140.521\ \mathrm{m^3}\times 3{,}0022\ \mathrm{kWh/m^3} = 9.428\ \mathrm{MWh}\qquad \frac{E_{bv}}{P_{max}} = \frac{9.428}{11{,}99} \approx 786\ \mathrm{h}
La relación energía/potencia de 786 horas define su función: desplazamiento intradiario y estacional, no regulación en segundos, que corresponde a la inercia. Respaldo: Atia y Fthenakis (2019), sobre la ósmosis inversa como carga flexible; Elimelech y Phillip (2011), sobre el consumo energético de la desalación.
La industria del clúster necesita 16.498 t/año de carbonato sódico para el vidrio y la extracción de alginato, y 17,71 t/h de vapor a 12 bar para la papelera, los evaporadores de salmuera y la digestión, sin sobredimensionar reactores y sin quemar combustibles.
La carbonatación cáustica —absorber CO₂ en sosa— exigiría 15.866 t/año de NaOH. Producir esa sosa por electrodiálisis bipolar obligaría a multiplicar el tren de electrodiálisis por once y generaría una cantidad equivalente de ácido clorhídrico sin destino. Para el vapor, la solución estándar es una caldera de gas o de fuelóleo, que el clúster no puede alimentar sin importar combustible fósil.
La ruta Solvay aprovecha tres corrientes que el clúster ya produce: la salmuera saturada de las aguas madres, el CO₂ geogénico de la calcinación de la caliza en el horno de cal y el amoníaco, que circula en bucle cerrado y sólo requiere una reposición de 41,25 t/año. La ruta no consume sosa cáustica. La sosa que necesitan otros procesos, 3.409 t/año, se cubre con la electrodiálisis bipolar (1.432 t/año) y con un suministro de proximidad declarado (1.978 t/año).
El horno de vidrio carga el carbonato cálcico crudo —caliza de cantera y conchas de molusco—, sin calcinación previa: el carbonato se descompone dentro del propio fundido, y la cal viva del horno S02 queda para el hempcrete, la regeneración del amoníaco y el Tecnosol.
El anillo de vapor, de 2,99 km a 12,0 bar, recoge el calor de los humos de los hornos de cal y de vidrio (7,01 t/h) y lo completa con una caldera de electrodos de 7,11 MW, operada como carga firme no modulable: modularla con el excedente renovable despresurizaría el colector cada noche. El enclavamiento homeostático prohíbe la parada simultánea de los dos hornos, porque la inercia de la red caería a 4,45 s, por debajo del mínimo de 4,50 s.
\mathrm{2\,NaCl + CaCO_3 \rightarrow Na_2CO_3 + CaCl_2}\qquad \dot Q = \dot m\,(h_g - h_{fw}),\quad \Delta h \approx 2{,}345\ \mathrm{MJ/kg}
Estequiometría: 16.498 t de Na₂CO₃ equivalen a 155.656 kmol/año, que exigen 6.850 t/año de CO₂ frente a las 8.804 t/año que libera la calcinación. Balance de vapor: 7,01 t/h recuperadas y 10,70 t/h aportadas por la caldera, con el mismo salto entálpico en los dos lados. Respaldo: Steinhauser (2008); Tong y Elimelech (2016).
El cáñamo industrial produce cada año decenas de miles de toneladas de cañamiza leñosa. Un material que se acumula se degrada, emite y deja de ser un sumidero.
Vender la biomasa como materia prima agrícola sin transformar traslada fuera del territorio su valor y su carbono. Fabricar material de construcción sin un destino asegurado produce acopios en campa, donde la cañamiza se descompone y el producto pierde calidad.
La cañamiza destinada a construcción se transforma en la planta local de prefabricados en bloques de hempcrete, un biocompuesto de cañamiza y cal aérea hidratada del horno propio. La producción se fija por contrato a la colocación: 96.277 m³/año, a razón de 70,02 m³ por vivienda, en 1.375 viviendas bioclimáticas al año del programa social de bioconstrucción. La cola de material sin colocar es nula por diseño. Una vez en el muro, la cal recarbonata absorbiendo CO₂ del aire durante la vida útil del edificio.
El cierre biológico tiene su frontera toxicológica: los lodos de la depuradora no entran en el suelo ni en ninguna cadena alimentaria. El protocolo TEC-02 los piroliza a 550–650 °C y confina el biochar al 100 % en matrices de geopolímero incorporadas a la obra civil.
\mathrm{Ca(OH)_2 + CO_2 \rightarrow CaCO_3 + H_2O}\qquad V_{hempcrete} = 1.375 \times 70{,}02 = 96.277\ \mathrm{m^3/año}
La recarbonatación del hempcrete colocado aporta −7.822,49 tCO₂e/año, con un descuento de permanencia del 30 %, y el biochar confinado −2.822,00 tCO₂e/año; junto con el carbono orgánico del Tecnosol y el balance de proceso, el balance neto del clúster es de −25.123,56 tCO₂e/año. Respaldo: Amziane y Arnaud (2013); EN 15804:2012+A2:2019, que exige informar el carbono biogénico en su indicador propio; Davidovits (2008).
Construir una infraestructura agroindustrial de 2.824,88 MUSD sin exponer el proyecto a una deuda que su operación temprana no puede servir.
La estructura convencional financia la obra completa desde el primer día con un crédito único de gran tamaño. Sobre la misma base física del clúster, un esquema de ese tipo alcanza un pico de deuda de 1.941,94 MUSD en el Año 10 y acumula caja ociosa mientras paga intereses sobre ella. Además, exige a un único prestamista asumir el riesgo de construcción de todo el sistema antes de que exista un solo ingreso.
El clúster se despliega como una sucesión ecológica: los estratos basales, de alto retorno y sin dependencias de otros procesos productivos, entran primero y financian a los siguientes. La halita lavada AG-10 (Q01, 47,86 MUSD/año de EBITDA) y la microalga Dunaliella salina (Q05, 16,74 MUSD/año) operan desde el trimestre 12 y transforman el primer residuo del sistema, la salmuera, en su primera caja.
El ticket de entrada de la Fase 1 es de 639,49 MUSD: 191,85 MUSD de capital promotor y 447,65 MUSD del tramo concesional A1. Hasta el primer ingreso, en el trimestre 8, sólo se necesitan 281,76 MUSD. El segundo tramo concesional, de 533,29 MUSD, se contrata en el Año 6 con las primeras factorías ya en facturación. El 58,0 % del CAPEX —1.549,51 MUSD— se paga con EBITDA reinvertido; el pico de deuda viva es de 718,05 MUSD, el DSCR anual no baja de 1,33 en operación y la deuda queda amortizada en el Año 21.
Ese 1,33 se lee junto con la regla que lo hace posible: la inversión sólo se ejecuta cuando la caja la sostiene, y lo que no cabe se aplaza (Apéndice 17, §17.3). Las pruebas de estrés del mismo apéndice muestran por qué importa. Con el calendario de obra rígido, un EBITDA un 15 % menor, un CAPEX un 15 % mayor, ambos a la vez o un retraso de cuatro trimestres en los ingresos hunden el DSCR anual mínimo hasta valores entre 0,06 y 0,23 (§17.2). La viabilidad del caso base depende, por tanto, de esa flexibilidad del calendario: el plan de obra es un techo, no una obligación.
El modelo distingue con precisión dos proposiciones. La primera es física —un sistema sin fugas de materia se financia a sí mismo una vez arrancado—, y las cuentas del modelo muestran que es coherente: según ellas, la fase final se autofinancia sin un dólar externo desde el Año 14. La rentabilidad del arranque para el capital es otra proposición: la rentabilidad del activo sin apalancamiento es del 4,91 %, y la del accionista alcanza el 7,05 % porque la deuda concesional cuesta un 3 %.
Autofin(t) = \min\bigl(CAPEX(t),\ \max(0,\ Pre(t) - C_{min})\bigr),\qquad C_{min} = 40\ \mathrm{MUSD}
La caja operativa por encima del colchón es la primera fuente de cualquier inversión; la financiación externa sólo cubre el remanente. Respaldo: Yescombe (2014), sobre la estructura de la financiación de proyectos; OCDE (2018), principios de la financiación combinada para el desarrollo.
Un clúster de 37 nodos acoplados propaga cualquier parada: si la desaladora se detiene, se detienen el riego, el suelo, la papelera y el vidrio; si un horno dispara, el anillo de vapor pierde presión en dos minutos. Y sin acumuladores químicos ni bombeo reversible, que la regla de diseño excluye, falta la herramienta habitual para absorber esos golpes.
La respuesta convencional es una reserva grande y única que cubra cualquier déficit. Concentra el riesgo en un solo equipo, cuesta lo mismo tanto si el golpe llega como si no, y no aprovecha los excesos: la energía que sobra una noche de viento se vierte o se almacena con pérdidas.
La ley de la variedad requerida (Ashby, 1956) establece que un regulador sólo absorbe tantas perturbaciones distintas como estados distintos puede adoptar. ECPRYM no tiene una reserva grande: tiene muchas palancas pequeñas e independientes —electrólisis, desaladora, caldera de electrodos, modulación del horno de cal, grupo despachable de biometano, enlace de 132 kV— y diez existencias físicas. Cada perturbación se asigna a la palanca cuya constante de tiempo le corresponde: la masa rotante a los primeros 250 milisegundos, los variadores al segundo, los relés de subfrecuencia al décimo de segundo siguiente, la caldera a los dos minutos y los búferes a los días.
Un sistema es antifrágil cuando gana más con los picos de lo que pierde con los valles (Taleb, 2012). Sin almacenamiento electroquímico, esa convexidad se fabrica con existencias: cada megavatio hora que sobra se convierte en caliza machacada, cal viva, sal lavada, hipoclorito o agua embalsada, y cada existencia devuelve después potencia en las horas escasas, porque permite parar su proceso sin parar a los procesos que dependen de él. El silo de cal viva incorpora 11,4 GWh: es el mayor almacén energético implícito del clúster. La microalga Dunaliella salina ofrece la forma biológica del mismo principio: el β-caroteno que se vende se sintetiza precisamente como respuesta al exceso de luz y de salinidad (Ben-Amotz y Avron, 1983).
El tercer principio es de red. La carga firme de procesos continuos que hay que sostener en todo momento es de 28,14 MW, y con el bombeo en su mínimo escalonado el mínimo técnico nocturno consolidado del clúster es de 29,34 MW. Mientras el enlace de 132 kV está cerrado, la frecuencia la impone el sistema interconectado y cualquier pérdida de generación dentro del margen del enlace se absorbe sin excursión. La inercia local —494 MW·s— sólo gobierna cuando el enlace se abre. Por eso el clúster opera en dos modos, conectado e isla, con reglas distintas para cada uno. El propio enlace se construye con dos circuitos: la pérdida de uno no interrumpe la exportación, y la de ambos, que es un evento extremo, fija la exportación máxima en 39,52 MW.
Autonomía de un búfer ante la parada de su proveedor:
t_{aut} = \frac{V_{útil}}{\sum_j \dot m_j}
Con la prioridad de consumidores, la autonomía de un hijo pequeño se multiplica: la torre Solvay, que usa el 1,06 % de las aguas madres, dispone de 2.826 días de un estanque que la halita agota en 30.
Mayor escalón de potencia admisible en isla y mayor exportación por el enlace de doble circuito, con el límite extremo de 2,0 Hz/s:
\Delta P_{max} = \frac{2\,E_{cin}\,\mathrm{RoCoF}_{max}}{f_0} = \frac{2 \times 494 \times 0{,}94}{50} = 18{,}57\ \mathrm{MW}\qquad P_{exp,max} = \frac{2 \times 494 \times 2{,}0}{50} = 39{,}52\ \mathrm{MW}
Capacitancia del anillo de vapor, que fija el tiempo disponible para reponer el vapor de un horno:
C = V\,\frac{d\rho}{dP} = 211\ \mathrm{m^3} \times 0{,}49\ \mathrm{kg/(m^3\,bar)} \approx 104\ \mathrm{kg/bar}
Respaldo: Ashby (1956), sobre la variedad requerida; Taleb (2012), sobre la convexidad ante la perturbación; Kundur (1994), sobre la ecuación de oscilación de la máquina síncrona.
Las siete decisiones comparten una lógica: asignar cada función al recurso cuyo principio físico la resuelve con la menor conversión de energía y de materia. El mar hace lo que sólo el mar puede hacer, y el interior, lo que exige aire seco. La masa rotante sostiene los milisegundos; el agua almacenada, las horas y las estaciones. Las corrientes residuales se convierten en reactivos, vapor y muros. Y la caja sigue el mismo principio que la materia: lo que un estrato produce alimenta al siguiente. La continuidad sigue la misma lógica que el diseño: ninguna reserva única, sino muchas existencias pequeñas que convierten los excesos de hoy en la autonomía de mañana.
Las siete decisiones no son siete elecciones independientes: son aplicaciones de un mismo criterio de diseño, el criterio de cierre metabólico, que conviene dejar escrito porque es lo que permite discutirlas con física y no con preferencias.
El marco es el de los sistemas abiertos lejos del equilibrio. La entropía de un sistema abierto se descompone en el intercambio con el entorno y la producción interna irreversible:
\frac{dS}{dt} = \frac{d_eS}{dt} + \frac{d_iS}{dt}\qquad \frac{d_iS}{dt} \ge 0
Sólo el primer término puede ser negativo. Un sistema que recibe energía de alta calidad y devuelve calor degradado importa orden y puede sostener una organización estable: es una estructura disipativa (Nicolis y Prigogine, 1977), lo mismo que una célula (Schrödinger, 1944) o un bosque. La economía no escapa a esa contabilidad: el proceso económico no crea ni destruye materia ni energía, sólo degrada su disponibilidad (Georgescu-Roegen, 1971), y la calidad de cada corriente puede medirse (Odum, 1996).
Criterio central de cierre metabólico, en lenguaje llano. Un bosque en el desierto no vive porque alguien lo riegue cada día: vive porque cada estrato se alimenta de lo que el anterior desecha y porque el sol paga la factura. Si los procesos industriales se ordenan con esa lógica, el sistema necesita del exterior, en lo esencial, energía solar y agua de mar, además del capital y el trabajo que lo construyen y lo operan, y su estabilidad no la da un almacén de electricidad, sino la masa que ya gira en sus máquinas y el agua que ya está elevada.
Enunciado formal. Un sistema de n procesos con matriz de requerimientos internos A \ge 0 admite producción no negativa y finita para toda demanda final no negativa si y sólo si su radio espectral cumple \rho(A) < 1: es la condición de Hawkins y Simon (1949), que hace al sistema productivo siempre que reciba sus insumos primarios —energía, agua de mar, capital y trabajo—, que quedan fuera de la matriz. La productividad es necesaria para el cierre, pero no suficiente: para hablar de cierre hay que contabilizar también los residuos, como hace el análisis insumo-producto ampliado a residuos (Nakamura y Kondo, 2002), y comprobar que lo único que entra del exterior es sol, agua de mar y capital. Si además la variedad de control iguala o supera la variedad de las perturbaciones creíbles, y la energía cinética de las máquinas síncronas mantiene el RoCoF y el nadir de la mayor pérdida creíble dentro de los límites de operación, el cierre es estable sin acumulación electroquímica.
Qué aporta el criterio. La diferencia entre productividad y cierre se entiende con una fábrica muy eficiente: puede atender todo lo que se le pide y, a la vez, seguir llenando un vertedero. Por eso la matriz no basta: el cierre se comprueba, además, con el balance de residuos de cada nodo. Las tres condiciones del criterio proceden de campos distintos —el análisis insumo-producto, la cibernética (Ashby, 1956) y la estabilidad de los sistemas eléctricos— y ninguna es nueva. Lo propio de ECPRYM es reunirlas en un criterio de diseño y exigir que el clúster las cumpla a la vez.
Leídas desde aquí, las siete decisiones se ordenan solas: las decisiones 1, 4 y 5 construyen la matriz A con \rho(A) < 1, que hace productiva la red, y el balance de residuos que exige el cierre, en el que cada residuo es insumo de otro proceso; las decisiones 2 y 3 aportan la tercera condición, la inercia de los segundos y el agua de las horas; y las decisiones 6 y 7 aportan la segunda, la variedad de control: caja que se reinvierte antes de endeudarse y existencias que absorben lo que no se puede prever.
Balance de materia en la forma de Leontief (1986), con x el vector de producción de los n procesos, A = (a_{ij}) la matriz de requerimientos —a_{ij} es lo que el proceso j consume del producto i por unidad de su propio producto— y d la demanda final:
x = A\,x + d \quad \Longleftrightarrow \quad (I - A)\,x = d \quad \Longrightarrow \quad x = (I - A)^{-1} d
Que el sistema sea productivo exige que (I-A)^{-1} exista y sea no negativa. Por el teorema de Perron-Frobenius, toda matriz no negativa tiene un autovalor real igual al máximo módulo de sus autovalores, su radio espectral \rho(A), y la serie de Neumann
(I - A)^{-1} = I + A + A^2 + A^3 + \cdots
converge si y sólo si \rho(A) < 1; al converger, todos sus términos son no negativos, luego la inversa lo es (Horn y Johnson, 2013). El significado físico es directo: A^k es lo que el sistema consume de sí mismo en la vuelta k del bucle de reutilización, y \rho(A) < 1 dice que cada vuelta pide menos que la anterior, de modo que la suma es finita. Con \rho(A) \ge 1 la serie diverge y el bucle se realimenta: si todos los procesos están conectados entre sí —una matriz irreducible—, ninguna demanda final no negativa y no nula admite un plan de producción no negativo; las soluciones algebraicas que puedan existir exigen producciones negativas, sin sentido físico. Ese sistema no es productivo —y, por tanto, tampoco puede cerrar—, por mucho que se le añada capacidad. Es la razón por la que la clausura no se consigue agregando plantas a un polígono convencional.
Condición de gobierno. La ley de la variedad requerida (Ashby, 1956) establece que sólo la variedad absorbe variedad: si el entorno presenta V_e estados distinguibles y el regulador dispone de V_r, la variedad del resultado no baja de V_e / V_r. De ahí que cada perturbación creíble necesite su propio grado de libertad, y no un margen genérico: cargas modulables para la variabilidad del recurso, existencias para la interrupción de suministro, dos modos de red para la pérdida del enlace. En este estudio la condición funciona como guía cualitativa: todavía no se define cómo medir la variedad del clúster, de modo que se aplica como esa regla de diseño —un grado de libertad por perturbación creíble— y no como una desigualdad calculada.
Condición dinámica. Ante una pérdida de generación \Delta P, la pendiente inicial de la frecuencia y la constante de inercia son:
\mathrm{RoCoF} = \frac{f_0\,\Delta P}{2\,E_{cin}} \qquad H = \frac{E_{cin}}{S_{base}}
Con E_{cin} = 494 MW·s y S_{base} = 107{,}9 MVA resulta H = 4{,}58 s. Según el modelo, la mayor pérdida creíble al mediodía, 18,56 MW, produce un RoCoF de 0,94 Hz/s y un nadir de 49,52 Hz: ambos dentro de los límites de operación. El transitorio se sostiene, por tanto, con masa rotante.
Las tres condiciones son independientes, y el criterio las exige a la vez: la primera hace que el metabolismo tenga solución, la segunda que pueda gobernarse y la tercera que sobreviva al primer segundo de la perturbación. Sólo la primera es un resultado matemático; la segunda y la tercera se comprueban con el modelo del clúster, y el cierre mismo, con el balance de masa de cada nodo. Ninguna de las tres invoca acumulación electroquímica, y por eso el diseño puede prescindir de ella.
Toda magnitud de esta obra procede de dos fuentes y de ninguna más: la matriz canónica de parámetros del clúster y el gemelo digital que la resuelve. La matriz es un único fichero de datos, congelado en memoria, que ningún componente de la plataforma puede alterar mientras la página se usa. El gemelo digital resuelve el sistema físico completo —oscilación de la máquina síncrona, balance de agua, metabolismo de materia y cascada de caja— y comprueba que cada cifra de la matriz es solución de ese sistema y no una elección editorial. El orden es siempre el mismo: una decisión de diseño cambia la matriz, el gemelo la vuelve a resolver, las herramientas de gobierno contrastan cifra a cifra lo que la plataforma publica y sólo entonces cambian los textos. No existe el camino inverso.
El libro de cálculo de prefactibilidad que fijó el dimensionado inicial es un documento fundacional. Se conserva íntegro y se abre en modo de sólo lectura; cada herramienta verifica su huella digital al empezar y al terminar, de modo que ninguna ejecución puede modificar una sola celda. Su función es dejar constancia del punto de partida, no alimentar lo que se publica: cuando el gemelo digital y el libro fundacional dan valores distintos para la misma magnitud, gobierna el gemelo, porque es el que resuelve el sistema con la física y el perímetro vigentes.
Esta jerarquía tiene una consecuencia práctica que conviene enunciar: la trazabilidad de una cifra publicada no termina en un documento, sino en una ecuación. Cualquier magnitud de la obra puede seguirse hasta la celda del libro fundacional que la origina o hasta el documento externo que la declara, y desde ahí hasta la ecuación del gemelo que la reproduce. Lo que no puede hacerse es al revés: escribir una cifra en un texto y esperar que el sistema se ajuste a ella.
Esa cadena es una verificación de consistencia interna, la propia de un estudio de prefactibilidad conceptual: comprueba que las cuentas cuadran, que los balances cierran y que cada cifra remite a su celda y a su ecuación. Verificar no es validar. Validar exige contrastar el modelo con observaciones independientes de él —datos de plantas de ósmosis inversa en operación, series meteorológicas de Lüderitz, rendimientos agrícolas publicados y, en último término, un piloto—, y ese contraste está pendiente. Es el paso que convertiría la coherencia de las cuentas en evidencia sobre el terreno.
Los decimales piden la misma lectura. Sirven para la trazabilidad —permiten encontrar cada cifra en su celda—, no expresan la precisión de una medida. En fase conceptual, una estimación de coste de clase 5 según AACE International (2020) puede desviarse entre un −20 % y un −50 % a la baja y entre un +30 % y un +100 % al alza; una de clase 4, entre un −15 % y un −30 % y entre un +20 % y un +50 %. Los 2.824,88 MUSD de CAPEX dicen de dónde sale la estimación, no cuánto costará la obra al dólar.