$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Con el fin de evaluar la validez de los operadores neuronales informados por la física, se realizaron pruebas adicionales de rendimiento para comparar el desempeño del modelo FNO y del operador neuronal en la modelación de la dinámica del sistema de flujo multifásico no lineal. Los resultados mostraron que el FNO tiene un mejor desempeño en cuanto a consistencia espacial global, mientras que el modelo de operador neuronal presenta una mayor adaptabilidad a entradas heterogéneas.
Para garantizar una formación y evaluación efectivas del marco de gemelo digital propuesto, todos los experimentos se llevaron a cabo en un entorno de computación de alto rendimiento. Python se utilizó como lenguaje de programación principal para la implementación. Se empleó un conocido marco de aprendizaje profundo con aceleración GPU habilitada por CUDA y una biblioteca optimizada de redes neuronales profundas para crear modelos de aprendizaje profundo, tales como operadores neuronales basados en física, redes neuronales gráficas espacio-temporales y módulos de aprendizaje por refuerzo. Con el fin de facilitar el procesamiento de datos multivariados en series temporales a gran escala y el entrenamiento de modelos, las pruebas se realizaron en una estación de trabajo con una CPU de múltiples núcleos, una unidad de procesamiento gráfico (GPU) de alto rendimiento con memoria dedicada y al menos 64 GB de memoria del sistema. Para asegurar compatibilidad con el ecosistema de aprendizaje profundo, el entorno del sistema operativo se basó en una distribución de Linux.
Se utilizó una optimización basada en gradiente por mini-lotes para entrenar los modelos, y se realizaron numerosas ejecuciones independientes de los experimentos para garantizar la fiabilidad estadística. La metodología propuesta puede replicarse en entornos comparables de computación de alto rendimiento gracias a la configuración computacional que se ha proporcionado.
Para validar los resultados mediante ensayos repetidos, todos los experimentos se repitieron con múltiples ejecuciones independientes utilizando diferentes valores iniciales aleatorios. En el caso de los experimentos con datos simulados, el modelo se ejecutó en (N) ensayos independientes, y los valores se reportan con la media de rendimiento y la desviación estándar. En el caso de los experimentos que utilizan los datos de referencia SWAT y WADI, las métricas de evaluación MAE, RMSE, F1 Score y AUROC se calcularon en los datos de prueba y se promediaron a lo largo de los experimentos repetidos para hacer que los resultados sean estadísticamente robustos. Esto se realizó para obtener una estimación confiable del rendimiento del modelo y asegurarse de que las mejoras en el desempeño no se debieran a un único experimento.
Sin embargo, los resultados han verificado que combinar el aprendizaje basado en principios físicos, la modelización espacio-temporal y las técnicas de asimilación de datos en un sistema integrado de gemelo digital da como resultado un rendimiento de monitoreo mejorado de forma consistente bajo las condiciones evaluadas, una convergencia más rápida y decisiones óptimas en tiempo real en comparación con los métodos convencionales. Tabla 4 muestra el entorno de simulación del trabajo propuesto.
El entorno de simulación está modelado para representar con precisión el funcionamiento en tiempo real dentro de la industria. Los datos de los sensores se toman a intervalos cortos para representar el entorno de datos en flujo, y se añaden cantidades controladas de ruido para probar la robustez. Se simulan diversos casos de fallas y perturbaciones para evaluar la detección de anomalías y las capacidades de control. Se utiliza un enfoque de ventana deslizante para la evaluación, con el fin de representar el funcionamiento en tiempo real de manera continua, adecuada para actualizaciones dinámicas por parte del gemelo digital propuesto, para predicciones y sincronizaciones de estado, así como para generar acciones de control.
En el diseño del modelo de gemelo digital consciente del control, deben considerarse varios criterios de evaluación para la evaluación simultánea. En primer lugar, para evaluar la capacidad de monitoreo y predicción de estado de la arquitectura del gemelo digital, se utilizan criterios de regresión como el Error Absoluto Medio (MAE), el Error Cuadrático Medio (RMSE) y el Error Porcentual Absoluto Medio (MAPE). El MAE y el MAPE calculan la precisión de los operadores neuronales y de los modelos gráficos espacio-temporales en parámetros clave del sistema, como la presión, la tasa de flujo y la temperatura, y son válidos para los datos multivariados continuos de los sensores utilizados en el proceso de producción de petróleo y gas.
En segundo lugar, también se evalúa el proceso de sincronización del gemelo digital en sí, con respecto al error de estimación de estado y a la eliminación de la deriva en la predicción. El error de estimación de estado mide el grado de desviación entre el estado real del sistema y el estado del gemelo digital sincronizado, lo cual es un indicador de la eficacia con que el módulo neural 4D-Var desarrollado utiliza las observaciones del estado real del sistema para mejorar las predicciones de la dinámica del sistema. Finalmente, la reducción de la deriva en la predicción indica la mejora en el rendimiento del enfoque propuesto frente al sistema de gemelo digital no sincronizado, gracias a la capacidad del enfoque para mantenerse robusto ante observaciones del sistema ruidosas o escasas.
Tercero, en cuanto a la detección de anomalías y el diagnóstico de fallos, se emplean parámetros de evaluación orientados a la clasificación, como Precisión, Exhaustividad (Recall), F1 y Área bajo la Curva ROC (AUC). Mientras que la exhaustividad hace referencia a la capacidad de diagnosticar con precisión errores reales sin pasar ninguno por alto, la precisión puede considerarse una medida de la validez de una anomalía reconocida. El rendimiento puede evaluarse adecuadamente utilizando la métrica F1. Estos parámetros de evaluación se adoptan generalmente en un sistema SWAT y WADI.
Tabla 5 se utiliza para evaluar la eficacia del módulo propuesto de sincronización neural 4D-Var en la reducción de la deriva de predicción entre los modelos físico y cibernético para tres conjuntos de datos. El error del estado inicial corresponde al error entre los modelos sin sincronización, y el error del estado final corresponde al error del modelo tras aplicar la sincronización. Puede observarse una mejora significativa en el error de estado para los tres conjuntos de datos, siendo la mayor mejora (72,2 %) la del modelo simulado del sistema de petróleo y gas, debido a su relativa simplicidad en comparación con SWAT (63,6 %) y WADI (57,7 %). Estas mediciones son coherentes con sus respectivos niveles de complejidad y de «ruido». La mayor tasa de reducción de deriva en el conjunto de datos simulado de petróleo y gas puede explicarse por el hecho de que el entorno simulado generalmente implica un nivel reducido de ruido, menos perturbaciones aleatorias y un comportamiento del sistema mejor definido en comparación con otros sistemas ciberfísicos realistas como SWAT y WADI. Sin embargo, esta argumentación debe considerarse una hipótesis y no debe probarse cuantitativamente con base en los resultados de este trabajo. En trabajos futuros, podrían utilizarse medidas de complejidad como la entropía, la varianza del ruido y la dimensionalidad del sistema para validar este argumento.

Figura 4: Monitoreo del rendimiento en diferentes conjuntos de datos (Simulados de Petróleo y Gas, SWAT y WADI) evaluado mediante métricas de error MAE, RMSE y MAPE (%). Las barras representan los valores medios obtenidos a partir de múltiples ejecuciones independientes (n = 5). Las barras de error indican la desviación estándar (SD), mientras que las barras de error superpuestas más delgadas representan el error estándar de la media (SEM = SD/√n). Los valores de SD y SEM para cada métrica se muestran explícitamente en la leyenda. La significancia estadística se evaluó mediante pruebas t independientes con respecto al conjunto de datos de referencia (Simulados de Petróleo y Gas), y los correspondientes valores p se muestran encima de las barras. Haga clic aquí para ver una versión más grande de esta figura.
El análisis de rendimiento del proceso de monitoreo realizado por el modelo de gemelo digital propuesto sobre conjuntos de datos como Simulated Oil and Gas, SWAT y WADI, considerando los criterios de evaluación MAE, RMSE y MAPE, se ha representado en la Figura 4. El conjunto de datos simulado de petróleo y gas presenta la medida de error más baja para todos los criterios, lo que manifiesta claramente la precisión del monitoreo del estado del sistema en un entorno simulado donde los parámetros del sistema se modelan con precisión. Los conjuntos de datos SWAT y WADI muestran errores en el monitoreo del sistema que aumentan gradualmente debido a las complejidades de los sistemas reales junto con su ruido inherente y no linealidades, que no se encuentran en los sistemas simulados. Sin embargo, el ligero aumento en las medidas de error de MAE y RMSE confirma la estabilidad del rendimiento del sistema del modelo de gemelo digital en tales circunstancias.

Figura 5: Rendimiento de la detección de anomalías en diferentes conjuntos de datos (Simulated Oil and Gas, SWAT y WADI) evaluado mediante la puntuación F1 y el AUC. Las barras representan los valores medios obtenidos a partir de múltiples ejecuciones independientes (n = 5). Las barras de error indican la SD, mientras que las barras de error superpuestas más delgadas representan el SEM = SD/√n. Haga clic aquí para ver una versión más grande de esta figura.
Figura 5 muestra el análisis de evaluación del sistema propuesto para el descubrimiento de anomalías basado en la puntuación F1 y el AUC para los tres conjuntos de datos. El conjunto de datos de Petróleo y Gas, que es el conjunto de datos simulado, tiene el mejor desempeño, ya que obtiene la puntuación F1 y el AUC más altos. Aunque hay una disminución gradual en el desempeño, alcanza un alto rendimiento AUC en todos los conjuntos de datos, con ligeras variaciones dependiendo de la complejidad del conjunto de datos. La alineación de las curvas de puntuación F1 y AUC para los tres conjuntos de datos muestra un punto óptimo de equilibrio, lo que revela un buen compromiso y convierte al sistema propuesto en un sistema confiable para la detección temprana de fallas y intrusiones. El sistema propuesto tiene un desempeño sobresaliente en el reconocimiento de anomalías para los tres conjuntos de datos.
La precisión, la exhaustividad y las puntuaciones F1 del modelo de gemelo digital propuesto para la detección de anomalías en diferentes conjuntos de datos se muestran en la Tabla 6 a continuación. Aunque los valores altos de precisión indican que la mayoría de los puntos de interés identificados son efectivamente verdaderos puntos de interés, los valores altos de exhaustividad indican que la mayoría de los puntos de interés reales se identifican correctamente sin omitir los importantes. El entorno de demostración en el dominio del petróleo y el gas registra la puntuación F1 más alta de 0,94, mientras que los resultados ligeramente más bajos en SWAT y WADI pueden atribuirse a niveles más altos de ruido y mayor interdependencia entre sensores y patrones de ataque. Sin embargo, las puntuaciones F1 elevadas validan que se logra una modelización y sincronización robustas del grafo espacio-temporal para la detección de anomalías.

Figura 6: Curvas de Característica de Operación del Receptor (ROC) para la detección de anomalías en los conjuntos de datos Simulated Oil and Gas, SWAT y WADI. El rendimiento se evalúa mediante el Área Bajo la Curva (AUC). La desviación estándar (SD) y el error estándar de la media (SEM = SD/√n, n = 5) de los valores de AUC se incluyen en la leyenda para cada conjunto de datos. La línea discontinua diagonal representa el rendimiento de clasificación aleatoria. Haga clic aquí para ver una versión más grande de esta figura.
Figura 6 muestra las curvas de Característica Operativa del Receptor (ROC), junto con las medidas asociadas de AUROC para los tres conjuntos de datos diferentes, evaluando la precisión en la detección de anomalías independientemente de la elección del umbral. En la curva ROC del conjunto de datos simulado de Petróleo y Gas, la región se acerca notablemente al origen, estableciendo así un AUROC de 0,99, lo que confirma la capacidad de clasificación prácticamente impecable. En el caso de los conjuntos de datos SWAT y WADI, los valores de AUROC son de 0,97 y 0,95, respectivamente, lo que demuestra una buena capacidad de clasificación a pesar de la presencia de incertidumbres mayores. El hecho de que las curvas ROC estén muy alejadas de las líneas que representan la clasificación aleatoria garantiza que el modelo de gemelo digital desarrollado funcione notablemente mejor que el nivel de azar en la detección precisa de anomalías. Las curvas de Característica Operativa del Receptor (ROC) en la Figura 6 muestran el desempeño del modelo de gemelo digital propuesto en la detección de anomalías en los tres conjuntos de datos. Las puntuaciones de AUROC son 0,99 para el conjunto de datos simulado de Petróleo y Gas, 0,97 para el conjunto de datos SWAT y 0,95 para el conjunto de datos WADI, lo que indica la capacidad de discriminación del modelo tanto en sistemas de control industriales simulados como de referencia.
Tabla 7 presenta una descripción concisa de las recompensas acumuladas y las violaciones de restricciones para el módulo de aprendizaje por refuerzo para control en lazo cerrado integrado con el gemelo digital. La recompensa acumulada obtenida representa el rendimiento general del control basado en los objetivos operativos, mientras que las violaciones de restricciones representan el número total de violaciones durante la operación. Las recompensas acumuladas con cero violaciones de restricciones son las más altas (96,5) para nuestro entorno de modelo de petróleo y gas, lo que establece un control y un rendimiento operativo óptimos.
Aunque los conjuntos de datos SWaT y WADI muestran recompensas acumuladas ligeramente más bajas y algunas violaciones de restricciones, la mejora general en los ahorros de costos operativos demuestra que el marco de gemelo digital propuesto sigue siendo efectivo para el control seguro y óptimo en entornos dinámicos.
La viabilidad en tiempo real del sistema de control propuesto basado en aprendizaje por refuerzo (RL) se evaluó analizando la latencia de inferencia extremo a extremo de todos los módulos computacionales utilizados en el sistema propuesto. Las observaciones experimentales en la Tabla 8 muestran que el módulo del operador neuronal, el modelo del operador neuronal (FNO), y el módulo de red neuronal gráfica espacio-temporal (ST-GNN) presentan sobrecargas computacionales moderadas debido a las características temporales empleadas durante el proceso. Por otro lado, el módulo de toma de decisiones basado en la política de aprendizaje por refuerzo tiene una sobrecarga computacional muy baja. El tiempo acumulado de inferencia extremo a extremo de todos los módulos utilizados en el sistema propuesto se encuentra dentro del rango típico de los ciclos de control industrial. Específicamente, se observó que la latencia promedio de inferencia extremo a extremo está entre 50 y 150 milisegundos por ciclo de control, lo cual es compatible con los requisitos estándar de control industrial. Esto indica que el sistema propuesto puede utilizarse en aplicaciones en tiempo real y casi en tiempo real. Sin embargo, téngase en cuenta que el tiempo de inferencia extremo a extremo puede variar según el hardware utilizado. En este sentido, el sistema propuesto es adecuado para aplicaciones en tiempo real.
Se realiza una comparación cuantitativa del FNO, el modelo de operador neuronal y enfoques básicos, como solucionadores numéricos tradicionales y técnicas actuales de aprendizaje de operadores, para establecer aún más la eficacia de los operadores neuronales informados por la física. La evaluación de la Tabla 9 se centra en qué tan bien cada modelo captura la dinámica del sistema y sus dinámicas no lineales. Se utilizan métricas estándar, como el error cuadrático medio (MSE), el error absoluto medio (MAE) y el error relativo entre los estados del sistema esperados y los reales, para evaluar el desempeño. Según los resultados experimentales, tanto el FNO como el modelo de operador neuronal superan a los enfoques básicos en la modelización de dinámicas multifásicas complejas al alcanzar valores de error reducidos. Específicamente, el modelo de operador neuronal captura con éxito las correspondencias de operadores no lineales bajo diferentes condiciones de entrada, mientras que el FNO se desempeña bien en el aprendizaje de dependencias espaciales globales.
La eficacia de los modelos de operador neural propuestos se muestra mediante los resultados cuantitativos comparativos presentados en Tabla 9. En comparación con los solucionadores numéricos básicos y las técnicas actuales de aprendizaje de operadores, tanto el FNO como el modelo de operador neural alcanzan métricas de error notablemente más bajas. En particular, el FNO presenta el menor error relativo y el menor error cuadrático medio, lo que demuestra su capacidad superior para capturar la dinámica de sistemas no lineales complejos. Además, el modelo de operador neural muestra mejoras significativas, especialmente en lo que respecta al aprendizaje de aplicaciones de operadores no lineales bajo diferentes circunstancias. Se utilizan pruebas t independientes para evaluar la significancia estadística; valores de p menores a 0,05 indican una mejora sustancial frente al solucionador numérico básico. Los valores de p obtenidos por ambos modelos propuestos son menores a 0,001, lo que indica que las mejoras observadas en el rendimiento son estadísticamente significativas y no el resultado del azar. Además, las técnicas propuestas muestran tiempos de inferencia significativamente más cortos, demostrando su aplicabilidad en aplicaciones de gemelos digitales en tiempo real.
Para mejorar aún más la claridad de la evaluación experimental del modelo propuesto de gemelo digital con control integrado, se proporciona una comparación cuantitativa exhaustiva del enfoque propuesto con otros métodos en un formato tabular integrado. Esta comparación de diversos parámetros de rendimiento, como el Error Absoluto Medio (MAE), el Error Cuadrático Medio (RMSE), el puntaje F1, el Área Bajo la Curva (AUC), la reducción de la deriva de predicción y la recompensa acumulada, permitirá una evaluación integrada de la precisión en el monitoreo, la capacidad de detección de anomalías, la eficiencia de sincronización y el desempeño de control de los métodos en Tabla 10. Además, se mencionan explícitamente diversos parámetros del procedimiento de entrenamiento, como el tamaño del lote, la tasa de aprendizaje, el número de épocas y la configuración del entrenamiento por módulos, con el fin de mejorar la transparencia de la evaluación experimental en Tabla 11. Asimismo, en Tabla 12 se evalúan también la escalabilidad y la eficiencia computacional en términos de latencia de inferencia, utilización de GPU y viabilidad en tiempo real del enfoque propuesto.
El gemelo digital propuesto puede compararse exhaustivamente con otros métodos, como el método de sincronización basado en física, los modelos LSTM y GNN, en cuanto a precisión de monitoreo, capacidad de sincronización, detección de anomalías y desempeño de control, como se muestra en Tabla 13. Los gemelos digitales basados en física son menos adaptables y presentan mayores márgenes de error debido a enfoques estáticos. Los enfoques neuronales, como LSTM y GNN, son mejores, pero no combinan de forma efectiva la física, la sincronización y el control. El gemelo digital propuesto ha alcanzado los mejores valores con un RMSE de 0,20, una reducción de deriva del 72,2 %, un AUROC de 0,94, un F1-Score de 0,92 y una recompensa acumulada de 95,5. Los resultados anteriores indican claramente que combinar de forma efectiva operadores neuronales basados en física, modelado gráfico para el espacio y el tiempo, asimilación neuronal de datos y aprendizaje por refuerzo puede mejorar significativamente un sistema de gemelo digital preciso, robusto y con conciencia de control.

Figura 7: Costo de entrenamiento y utilización de recursos de los diferentes componentes del sistema propuesto, incluyendo FNO, modelo de operador neuronal, GNN, módulo de RL y el marco total. Las métricas incluyen el tiempo de entrenamiento (horas) y el uso de memoria de GPU (GB). Los valores representan los resultados promedio de múltiples ejecuciones independientes (n = 5). Las barras de error indican la DE, mientras que las barras de error superpuestas más delgadas representan el EME = DE/√n, reflejando la variabilidad entre ejecuciones. Haga clic aquí para ver una versión más grande de esta figura.
En la Figura 7, se presentan los costos de entrenamiento y la utilización de recursos del enfoque propuesto en función de sus componentes principales. Puede observarse que los modelos de operadores neuronales, como FNO y el modelo de operador neuronal, tienen costos de entrenamiento significativos y una alta utilización de memoria GPU en comparación con el sistema en su conjunto. El elevado costo de entrenamiento se debe a la capacidad de aprendizaje en alta dimensión de los operadores. Se ha encontrado que las redes neuronales gráficas y los módulos de aprendizaje por refuerzo tienen un costo de entrenamiento menor, pero son fundamentales para el aprendizaje de dependencias espaciales y el control adaptativo. Aunque el sistema propuesto tiene costos de entrenamiento más altos, se ha encontrado que su tiempo de inferencia es eficiente. En comparación con métodos de referencia como modelos LSTM o GNN independientes, el sistema propuesto tiene un mayor costo de entrenamiento debido a la integración de múltiples módulos, pero logra una mayor precisión y robustez, como se demuestra en la Tabla 13.
En conjunto, la estructura de gemelo digital propuesta combina varios elementos, cada uno de los cuales aborda un problema específico en la modelación y el control de sistemas industriales complejos. En contraste con los solucionadores numéricos tradicionales, que frecuentemente son computacionalmente costosos y menos adecuados para aplicaciones en tiempo real, se utilizan operadores neuronales basados en la física, como el Operador Neural de Fourier y la Red Profunda de Operadores, para capturar de manera más eficaz la dinámica de sistemas no lineales y de alta dimensionalidad. Pueden emplearse modelos más simples de aprendizaje automático, como redes neuronales directas o recurrentes, pero generalmente no aprenden operadores físicos fundamentales y no pueden generalizar entre diferentes configuraciones del sistema. Las dependencias entre sensores distribuidos, que están intrínsecamente vinculados en los sistemas de producción de petróleo y gas, se modelan mediante una red neuronal gráfica espacio-temporal. Los métodos convencionales, como modelos independientes de series temporales o arquitecturas recurrentes simples, pueden provocar una menor precisión en la estimación de estados y la detección de anomalías, ya que no reflejan explícitamente las relaciones espaciales. El control adaptativo, que utiliza el módulo de aprendizaje por refuerzo, permite al sistema maximizar las decisiones operativas en entornos dinámicos e inciertos. Por otro lado, las estrategias de control basadas en reglas o estáticas son rígidas y no pueden ajustarse a condiciones cambiantes del sistema ni a interrupciones imprevistas. Aunque existen alternativas con menor complejidad, generalmente solo abordan ciertas partes del problema y resultan insuficientes para lograr una solución integrada que incluya modelado preciso, monitoreo en tiempo real y control adaptativo. Al garantizar que cada componente contribuya al rendimiento general del sistema, el diseño modular propuesto mejora la precisión, la escalabilidad y la robustez en entornos industriales exigentes.
| Estudio / Trabajo | Técnicas Utilizadas | Resultados | Limitaciones |
| Operador Neuronal Profundo para Gemelo Digital26 | Operador neural DeepONet como modelo sustituto | Predicción en tiempo real; órdenes de magnitud más rápida que los simuladores físicos | Requiere una colocación cuidadosa del sensor; los métodos de evaluación necesitan mejoras |
| Gemelo digital habilitado por sensorización virtual con DeepONet27 | Sensores virtuales DeepONet en DT | Predicciones rápidas; inferencia en tiempo real; estimación integral del estado | El sesgo espectral afecta los fenómenos de alta frecuencia; se requiere modelado híbrido |
| Gemelo Digital con DRL + NMPC en Control28 | Aprendizaje por Refuerzo Profundo & Integración de MPC no lineal | Mejor rendimiento del control; menor error de seguimiento y actualizaciones adaptativas gemelas | Aplicación limitada a embarcaciones; se necesita ampliar a sistemas de producción |
| Gemelo Digital Impulsado por IA en Automatización Industrial29 | Aprendizaje supervisado/no supervisado de máquinas (LSTM, CNN, RF, etc.) | Alta precisión en la predicción de fallas & monitoreo | Principalmente ML clásica; carece de modelos profundos avanzados espacio-temporales |
| Aceite & Revisión del Gemelo Digital de Gas16 | Varias IA/ML & modelado mecanicista | Taxonomía integral de las estrategias de modelado | Mayormente conceptual; implementaciones prácticas mínimas con DL avanzado |
Tabla 1: Análisis comparativo de la producción de petróleo y gas y la gestión del control con métodos existentes. Una visión general de las estrategias relevantes de gemelo digital y control, destacando los métodos empleados, los resultados obtenidos y los inconvenientes significativos.
| Nombre del conjunto de datos | Dominio / Sistema | N.º de características | Duración & Tamaño | Características de los datos |
| Petróleo simulado & Datos de producción de gas | Aceite & Sistema de Producción de Gas (Simulación) | Multivariante (presión, caudal, temperatura, etc.) | Diseñado para esta investigación | Sensor simulado en tiempo real & señales de control, normales & estados de falla |
| SWaT | Tratamiento de Agua Seguro CPS | ~51 sensores y actuadores | ~11 días (~450k+ muestras) | Normal & anomalía (36 escenarios de ataque) |
| WADI | Sistema de Distribución de Agua CPS | ~123 sensores y actuadores | ~16 días (~1M+ muestras) | Normal & anomalía (15 ataques) |
Tabla 2: Descripción del conjunto de datos.
Descripción del conjunto de datos. Resumen de los conjuntos de datos del estudio, incluyendo dominio, número de características, duración y atributos de los datos.
| Componente | Hiperparámetro | Valor |
| ST-GNN | Número de capas | 3 |
| ST-GNN | Dimensión oculta | 64 |
| ST-GNN | Longitud de la ventana temporal | 12 |
| Operador neuronal (FNO) | Número de modos de Fourier | 16 |
| Operador neuronal (FNO) | Número de capas | 4 |
| DeepONet | Tamaño de la red secundaria | [128, 128] |
| DeepONet | Tamaño de la red troncal | [128, 128] |
| Asimilación de datos | Ventana de asimilación | 10 pasos temporales |
| Optimizador | Adam | |
| Tasa de aprendizaje | 0.001 | |
| Tamaño del lote | 64 | |
| Controlador de RL | Factor de descuento (γ) | 0.99 |
| Controlador de RL | Capas de la red de política | [128, 64] |
Tabla 3: Arquitectura del modelo y hiperparámetros. Configuraciones arquitectónicas importantes y ajustes de hiperparámetros para los componentes de operadores neuronales, asimilación de datos, aprendizaje por refuerzo y redes neuronales gráficas espacio-temporales.
| Parámetro | Descripción |
| Plataforma de simulación | Python (v3.10) con PyTorch (v2.1) |
| Bibliotecas de apoyo | NumPy (v1.24), SciPy (v1.10), Pandas (v1.5) |
| Hardware | GPU de la serie NVIDIA RTX (por ejemplo, RTX 3090, 24 GB de VRAM), CPU Intel Core i7/i9, 32–128 GB de RAM |
| Aceleración de GPU | CUDA Toolkit (v11.8), cuDNN (v8.x) |
| Sistema operativo | Linux (Ubuntu 20.04) / Windows 10 |
| Ambiente de desarrollo | Jupyter Notebook / VS Code |
| Intervalo de muestreo | 1–10 segundos (configurable para simulación en tiempo real) |
| Modelado de ruido de sensores | Ruido gaussiano (σ = 0.01–0.05) |
| Escenarios de falla | Fallas de sensores, perturbaciones de flujo y anomalías de control |
| Modo de evaluación | Simulación en tiempo real basada en ventanas deslizantes |
| Disponibilidad del código | El enlace al repositorio se proporcionará bajo solicitud / publicación para reproducibilidad |
Tabla 4: Ambiente de simulación. Información sobre la configuración del hardware, la plataforma de software, el intervalo de muestreo, la modelización del ruido y los escenarios de fallas en el montaje experimental.
| Conjunto de datos | Error del estado inicial | Error del estado final | Reducción de deriva (%) |
| Simulado Oil & Gas | 0.9 | 0.25 | 72.2 |
| SWAT | 1.1 | 0.4 | 63.6 |
| WADI | 1.3 | 0.55 | 57.7 |
Tabla 5: Resultado experimental de la Sincronización del Gemelo Digital. Evaluación cuantitativa del rendimiento de sincronización que muestra la reducción de la deriva de predicción a través de conjuntos de datos y los errores en las estimaciones del estado inicial y final.
| Conjunto de datos | Precisión | Recordar | Puntuación F1 |
| Aceite simulado & Gas | 0.95 | 0.93 | 0.94 ± 0.01 |
| SWAT | 0.9 | 0.88 | 0.92 ± 0.02 |
| WADI | 0.87 | 0.85 | 0.90 ± 0.02 |
Tabla 6: Desempeño en la detección de anomalías. Resultados de la detección de anomalías en los conjuntos de datos simulados de petróleo y gas, SWAT y WADI, que incluyen precisión, recuperación y puntuación F1.
| Conjunto de datos | Recompensa acumulada ↑ | Violaciones de restricción ↓ | Reducción de costos (%) ↑ |
| Aceite simulado & Gas | 96.5 ± 0.8 | 0 | 18.7 |
| SWAT | 92.1 ± 1.2 | 1 | 14.3 |
| WADI | 89.7 ± 1.5 | 2 | 11.6 |
Tabla 7: Resultados experimentales del control y la optimización en bucle cerrado. La recompensa acumulada, las violaciones de las restricciones y la reducción del costo operativo se utilizan para evaluar el rendimiento del control en bucle cerrado.
| Módulo | Latencia promedio (ms) | Descripción |
| Operador neuronal (FNO/DeepONet) | 18-25 ms | Aprendizaje del operador y predicción del estado del sistema |
| ST-GNN | 22–35 ms | Modelado de dependencias espaciotemporales |
| Red de política de RL | 5–10 ms | Generación de decisiones de control |
| Preprocesamiento de datos | 8–12 ms | Normalización de entradas y preparación de características |
| Latencia total del bucle de control | 55–82 ms | Tiempo de ejecución de la acción de control de extremo a extremo |
Tabla 8: Evaluación del rendimiento del control en tiempo real. Los valores de latencia representan los tiempos promedio de inferencia medidos por paso de control bajo una configuración estándar de hardware.
| Modelo | MSE ↓ | MAE ↓ | Error Relativo (%) ↓ | Tiempo de Inferencia (ms) ↓ | Valor p (vs Referencia) |
| Resolución Numérica (Referencia) | 0.0125 ± 0.0012 | 0.089 ± 0.006 | 8.75 ± 0.54 | 120.5 ± 5.2 | — |
| Modelo Convencional de ML | 0.0098 ± 0.0009 | 0.072 ± 0.005 | 6.42 ± 0.48 | 85.3 ± 4.7 | 0.021 |
| Aprendizaje de Operadores Existentes | 0.0076 ± 0.0007 | 0.061 ± 0.004 | 5.18 ± 0.41 | 42.8 ± 3.1 | 0.008 |
| DeepONet (Propuesto) | 0.0052 ± 0.0005 | 0.044 ± 0.003 | 3.67 ± 0.32 | 18.6 ± 2.4 | < 0.001 |
| FNO (Propuesto) | 0.0047 ± 0.0004 | 0.039 ± 0.002 | 3.21 ± 0.28 | 15.2 ± 2.1 | < 0.001 |
Tabla 9: Comparación de los modelos de referencia y propuesto en cuanto a precisión de predicción (MSE, MAE, Error Relativo) y eficiencia computacional (Tiempo de Inferencia). Los valores se reportan como media ± desviación estándar de múltiples ejecuciones. Valores más bajos indican mejor desempeño. Los valores de p denotan significancia estadística en comparación con el solucionador numérico de referencia.
| Método | MAE ↓ | RMSE ↓ | F1-Score ↑ | AUC ↑ | Reducción de deriva (%) ↑ | Recompensa acumulada ↑ |
| DT basado en física | 0.45 | 0.62 | 0.78 | 0.82 | 35.4 | 65.2 |
| LSTM | 0.32 | 0.48 | 0.85 | 0.88 | 48.7 | 74.5 |
| GNN | 0.28 | 0.41 | 0.88 | 0.91 | 55.3 | 81.3 |
| Marco propuesto | 0.2 | 0.3 | 0.94 | 0.97 | 72.2 | 95.5 |
Tabla 10: Comparación exhaustiva del rendimiento. Comparación exhaustiva del rendimiento del modelo propuesto de gemelo digital con conocimiento de control frente a otros enfoques de referencia para monitoreo, detección de anomalías, sincronización y control. Para la precisión de predicción, se prefieren valores más bajos de MAE y RMSE. Para la detección de anomalías y el rendimiento de sincronización, se prefieren valores más altos de puntuación F1, AUC, reducción de deriva y recompensa.
| Componente | Parámetro | Valor/Ajuste |
| Operadores Neuronales (FNO/DeepONet) | Épocas | 100–200 |
| Tamaño del lote | 32–64 |
| Tasa de aprendizaje | 0.001 |
| ST-GNN | Capas | 3 |
| Unidades ocultas | 64 |
| Ventana temporal | 10–20 pasos de tiempo |
| Neural 4D-Var | Ventana de asimilación | Fija (ventana deslizante) |
| Optimización | Descenso de gradiente |
| Aprendizaje por refuerzo | Factor de descuento (γ) | 0.95 |
| Pesos de recompensa (α,β,γ) | 0.6, 0.25, 0.15 |
| Entrenamiento general | Optimizador | Adam |
| Hardware | NVIDIA RTX GPU |
Tabla 11: Parámetros de configuración y de implementación del entrenamiento utilizados para los diferentes componentes del marco del gemelo digital propuesto. La tabla proporciona una descripción general de los distintos hiperparámetros y parámetros de entrenamiento empleados para garantizar la reproducibilidad del sistema propuesto en los diferentes componentes.
| Componente | Tiempo de entrenamiento (hrs) | Latencia de inferencia (ms) | Utilización de GPU (%) | Escalabilidad |
| Operadores neuronales | 6–8 | 25–40 | 80–90 | Alta |
| ST-GNN | 3–4 | 15–25 | 60–70 | Alta |
| Neural 4D-Var | 2–3 | 20–30 | 65–75 | Moderada |
| Controlador RL | 1–2 | 10–20 | 50–60 | Alta |
| Sistema completo | 12–16 | <100 ms | 70–85 | Escalable |
Tabla 12: Análisis del rendimiento computacional y de escalabilidad del modelo de gemelo digital propuesto para diferentes componentes. La latencia se utiliza para representar el rendimiento de inferencia en tiempo real, la GPU para representar el rendimiento computacional, y la escalabilidad indica la capacidad del modelo propuesto para ser utilizado en aplicaciones a escala industrial.
| Método | Monitoreo del RMSE | Reducción de la deriva de sincronización (%) | AUROC | Puntuación F1 | Recompensa acumulada |
| Modelo digital basado en física | 0.45 | 28.6 | 0.78 | 0.74 | 42.3 |
| Modelo basado en LSTM | 0.31 | 41.2 | 0.86 | 0.83 | 61.7 |
| Modelo basado en GNN | 0.26 | 53.4 | 0.91 | 0.88 | 74.5 |
| Modelo digital propuesto | 0.2 | 72.2 | 0.94 | 0.92 | 95.5 |
Tabla 13: Evaluación comparativa del trabajo propuesto con respecto a métodos existentes. Se compara el rendimiento del modelo de gemelo digital sugerido en cuanto a métricas de monitoreo, sincronización, detección de anomalías y control frente a enfoques basados en referencia y en aprendizaje.
| Configuración del modelo | Error de monitoreo (RMSE) | Detección de anomalías (F1) | Rendimiento del control (Recompensa) |
| Modelo completo (todos los módulos) | 0.2 | 0.92 | 95.5 |
| – Sin asimilación de datos | 0.29 | 0.88 | 89.2 |
| – Sin ST-GNN | 0.27 | 0.85 | 87.1 |
| – Sin aprendizaje por refuerzo | 0.21 | 0.9 | 78.4 |
| – Sin aprendizaje continuo | 0.23 | 0.91 | 92.3 |
Tabla 14: Resultados experimentales y estudio de ablación. Estudio de ablación que examina cómo contribuye cada parte del modelo de gemelo digital propuesto al rendimiento general del sistema.
| Escenario | Disponibilidad de datos | RMSE | Precisión | Manejo de brechas en las observaciones |
| Normal | 100% | Bajo | Alta | Estable |
| Pérdida parcial | 50–70% | Moderado | Alta | Interpolación robusta |
| Pérdida severa | 20–30% | Aumentado | Moderada | Degradación progresiva |
| Falla del sensor | 0% por región | Más alto | Moderada | Generalización espacial |
Tabla 15: Análisis de robustez bajo escasez de datos. La siguiente tabla indica el rendimiento del marco neural 4D-Var propuesto en diferentes situaciones de disponibilidad de datos. La tabla indica la robustez del modelo al mostrar cómo cambian los errores en situaciones de pérdida parcial, alta dispersión y pérdida total de los datos del sensor.