Este estudio no involucra a ningún participante humano ni sujetos animales. Los conjuntos de datos utilizados (SWaT y WADI) son conjuntos de datos de referencia disponibles públicamente, y no se incluyen datos sensibles ni personales.
El enfoque propuesto establece un gemelo digital consciente del control, que permite incorporar algoritmos de aprendizaje profundo de vanguardia, asimilación de datos y evaluación continua y optimización de los sistemas y el control de extracción de petróleo y gas. En primer lugar, se recopilan y modelan datos multimodales en tiempo real procedentes de infraestructuras de producción distribuidas geográficamente, utilizando una estructura de grafos que define la conectividad y las interacciones entre los componentes de producción. Luego, se aprovechan operadores neuronales basados en principios físicos para construir aproximaciones a nivel de operador de las ecuaciones que rigen el flujo multifásico no lineal, facilitando así simulaciones escalables dentro del gemelo digital. En esta investigación, el término «modelo de operador neuronal» se utiliza para referirse a la red profunda de operadores (Deep Operator Network). Con el fin de incorporar dependencias espaciales y temporales entre activos de producción distribuidos geográficamente, se entrenan redes neuronales gráficas espacio-temporales dentro de representaciones basadas en principios físicos. Se aplica un enfoque diferenciable de asimilación de datos que utiliza la red neuronal 4DVAR para mantener las predicciones de simulación del gemelo digital coherentes y sincronizadas con los datos de monitoreo en tiempo real, incluso en presencia de mediciones ruidosas e irregulares en el espacio. Sobre la base del gemelo digital sincronizado en tiempo real, se entrena una herramienta de optimización basada en aprendizaje por refuerzo para el control en tiempo real. Finalmente, el gemelo digital permanece preparado para la actualización continua de modelos y una escalabilidad adaptada durante períodos prolongados, basándose en actualizaciones continuas a partir de nuevas mediciones disponibles.
La arquitectura de gemelo digital propuesta mejora la arquitectura convencional de gemelo digital al incorporar inteligencia, adaptabilidad y conciencia del control. La canalización completa consta de cuatro módulos: (i) adquisición de datos/construcción del grafo, (ii) estimación del estado mediante operadores neuronales y ST-GNN, (iii) detección de anomalías mediante análisis basado en residuos, y (iv) control basado en aprendizaje por refuerzo.

Figura 1: Arquitectura general del gemelo digital consciente del control. El modelo propuesto de gemelo digital consciente del control integra la detección de anomalías, redes neuronales gráficas espacio-temporales, sincronización neuronal mediante 4D-Var, operadores neuronales informados por principios físicos y control en bucle cerrado basado en aprendizaje por refuerzo para sistemas de producción de petróleo y gas. Haga clic aquí para ver una versión más grande de esta figura.
Figura 1 representa el proceso funcional de un gemelo digital con conocimiento del control de producción de petróleo y gas como un sistema ciberfísico en bucle cerrado. El proceso funcional comienza a nivel del sistema de producción de petróleo y gas, en el que diferentes sensores distribuidos en pozos, tuberías, válvulas y unidades de procesamiento capturan datos sobre presión, caudal, temperatura y controles en un instante determinado. Los diversos sensores capturan los datos en distintos formatos, y esta información se envía al gemelo digital, lo que constituye la base del conocimiento organizacional en ese instante.
La metodología propuesta incluye los siguientes pasos:
Recolección de datos
Se obtiene información en tiempo real de sensores sobre presión, caudal, temperatura y entradas de control procedentes de entornos de producción de petróleo y gas. Se utilizan otros conjuntos de datos de referencia disponibles, como SWAT y WADI, para evaluar las capacidades de robustez y generalización. Para garantizar la reproducibilidad, cada paso de la metodología propuesta se especifica utilizando detalles técnicos.
[Consulte el Archivo Suplementario (S1) para la Generación del Conjunto de Datos y la Simulación.]
Se utilizan dos conjuntos de datos industriales de referencia públicos del ámbito de investigación para sistemas de control industrial con el fin de demostrar la capacidad del gemelo digital conceptual frente a datos temporales complejos y de alta dimensionalidad, así como frente a eventos inusuales: los conjuntos de datos de referencia Secure Water Treatment (SWaT) y Water Distribution (WADI)34. Ambos conjuntos de datos se derivan de bancos de pruebas ciberfísicos realistas e incluyen datos multivariantes de sensores y actuadores correspondientes a operaciones normales y anomalías/ataques provocados en dichos sistemas. Mediante el uso de los conjuntos de datos SWaT y WADI, es posible realizar pruebas interdominio y demostrar la generalización más allá de los bancos de pruebas ciberfísicos hacia sistemas industriales distribuidos más amplios, con dependencias relacionadas con la ubicación y el tiempo. La descripción del conjunto de datos se muestra en Tabla 2.
La investigación propuesta utiliza tres conjuntos de datos diferentes para la evaluación exhaustiva de la metodología de gemelo digital desarrollada que se muestra en Tabla 2. El conjunto de datos base está compuesto por datos simulados relacionados con la organización de producción en la fabricación de petróleo y gas, y tiene como objetivo reproducir mediciones operativas reales detectadas por el sistema, así como condiciones de control y/o fallos. Además, se introduce el conjunto de datos de Tratamiento Seguro de Agua (SWaT) como ejemplo de un sistema físico a escala reducida, que comprende mediciones de 51 sensores y actuadores durante aproximadamente 11 días, incluyendo tanto condiciones normales como anomalías y ataques creados artificialmente. Este conjunto de datos es ampliamente considerado como el punto de referencia estándar para la investigación. Otro conjunto es el conjunto de datos de Distribución de Agua (WADI), cuyas mediciones capturan el comportamiento del sistema de distribución de agua a escala reducida, que comprende aproximadamente 123 variables durante un período de 16 días, incluyendo diferentes ataques y condiciones normales de operación. El conjunto de datos contiene un total de 15 ataques diferentes. Los dos conjuntos de datos anteriores giran principalmente en torno a series de alta dimensionalidad a lo largo del tiempo y relaciones espaciales, como en los sistemas de control industrial.
Preprocesamiento de datos
Antes del entrenamiento del modelo, todos los conjuntos de datos pasan por un procedimiento exhaustivo de pretratamiento para garantizar la calidad, consistencia y reproducibilidad de los datos. Entre los procedimientos de preprocesamiento se incluyen la gestión de valores ausentes, el filtrado de ruido, la normalización y la escalación de características. La descripción detallada de la tubería de preprocesamiento se encuentra en el Archivo Suplementario (S2).
Construcción del grafo del sistema
En la metodología propuesta, el sistema de producción física de petróleo y gas se representa en forma de grafo, que sirve como una estructura matemática para describir las correlaciones espaciales entre los activos junto con la dinámica del estado del sistema a lo largo del tiempo (véase el Archivo Suplementario (S3) para las Formulaciones de Construcción del Grafo del Sistema).
Aprendizaje de la física mediante operadores neuronales
Para la representación correcta de las características no lineales y multidimensionales de los sistemas de flujo multifásico, el enfoque propuesto utiliza operadores basados en principios físicos en forma de Operadores Neuronales de Fourier (FNO) y modelo de operador neuronal. A diferencia del método tradicional de representar redes neuronales como aplicaciones puntuales, el concepto de operadores permite la generalización en términos de diferentes condiciones de contorno. Esta característica hace que este método sea más adecuado para sistemas de flujo multifásico utilizados en la producción de petróleo y gas, ya que los fenómenos físicos ocurren de manera continua. Consulte el Archivo Suplementario (S4) para ver las ecuaciones matemáticas de los operadores neuronales.
Para evaluar más a fondo la validez de los operadores neuronales informados por la física, se realizaron experimentos adicionales para evaluar su capacidad de aproximar aplicaciones operadoras no lineales bajo diferentes estados de flujo y entradas de control. Los resultados muestran que los operadores pueden aproximar con precisión la aplicación del operador no lineal, validando así su capacidad para modelar con exactitud la dinámica del sistema, especialmente en presencia de alta dimensionalidad y con un costo computacional menor que otros solucionadores. La evaluación cuantitativa de la precisión del aprendizaje del operador se realiza mediante el error cuadrático medio entre los estados del sistema predichos y simulados, demostrando un rendimiento de aproximación consistente en diferentes condiciones de flujo.

Figura 2: Flujo de trabajo del marco de trabajo del operador neuronal informado por la física. Un ejemplo del flujo de trabajo del operador neuronal informado por la física muestra cómo se utilizan las entradas de control, los datos de sensores y las condiciones de contorno para aprender la dinámica de flujo multifásico no lineal y generar predicciones rápidas del sistema a nivel de operador. Haga clic aquí para ver una versión más grande de esta figura.
En el proceso de trabajo de PINO presentado en la Figura 2, el proceso de aprendizaje de flujos multifásicos no lineales físicamente complejos en procesos de petróleo y gas tiene lugar de manera continua. Esto comienza en la capa de entrada, donde el modelo recibe datos heterogéneos sobre el sistema físicamente cambiante. Esta entrada adopta la forma de condiciones de contorno, que incluyen presión, temperatura y velocidad, así como datos en tiempo real procedentes de sensores ubicados dentro del sistema de petróleo y gas, además de variables de control como la apertura de válvulas, la velocidad de las bombas o la apertura de estranguladores. Todas estas entradas definen el funcionamiento del sistema físicamente cambiante.
El componente clave dentro del flujo de trabajo es el PINO, implementado mediante FNO y el modelo de operador neuronal. El método propuesto es diferente y distinto en el sentido de que utiliza operadores neuronales para aprender operadores entre espacios de funciones, a diferencia de las redes neuronales tradicionales, que suelen centrarse en mapeos puntuales entre entradas y salidas. En este nivel, el FNO es capaz de capturar las relaciones espaciales globales utilizando el dominio de Fourier, y el modelo de operador neuronal puede aprender representaciones de operadores mediante redes de tipo rama y troncal. El modelo aprendido se guía mediante restricciones físicas para asegurar que la dinámica predicha se comporte como se espera en el flujo multifásico. El resultado del operador neuronal adopta la forma de una solución de campo continuo de variables físicas importantes, como la presión, los caudales y la distribución de fases para toda la red de producción. Además, los valores de salida son precisos en todo el campo, incluyendo pozos, tuberías y otras unidades conectadas, en lugar de medirse en puntos específicos como los sensores. De este modo, la salida del operador neuronal proporciona la capacidad de generar rápidamente simulaciones altamente precisas que pueden aproximar la solución de ecuaciones diferenciales parciales sin necesidad de resolverlas explícitamente. Esta solución actúa como la solución rápida en el gemelo digital.
El proceso de aprendizaje del operador neural toma características a nivel de función y a nivel de sistema, que son extracciones de los procesos de producción de petróleo y gas. Las características consisten principalmente en variables de estado físico para la producción de petróleo y gas, que son campos continuos, temperatura, caudal, presión y fracciones de fase para petróleo, gas y agua. Las características también incluyen condiciones de contorno, que involucran presiones de entrada, caudales y posiciones de válvulas. Las entradas de control, que incluyen posiciones de válvulas estranguladoras y velocidad de bombas, también se consideran para predecir respuestas bajo diferentes escenarios operativos. El modelo de operador neural recibe las características físicas y los parámetros de control, los cuales son procesados por las redes troncal y de ramificación. Las características proporcionan información suficiente para que los operadores neuronales predigan relaciones operadoras precisas que dependen tanto de los procesos físicos como de los factores operativos, para simulaciones rápidas de gemelos digitales.
Operador neural de Fourier (FNO)
En esta investigación, se utiliza un operador neural de Fourier (FNO)36 para estimar eficazmente las aplicaciones operacionales que describen la dinámica no lineal de flujo multifásico en sistemas de producción de petróleo y gas. A diferencia de las redes neuronales tradicionales, que aproximan aplicaciones entrada-salida para un conjunto específico de puntos, el FNO propuesto aproxima aplicaciones entre espacios funcionales, lo que permite la generalización sobre diversas escalas y configuraciones de un sistema. En este caso, el FNO actúa trasladando las representaciones del estado del sistema desde el espacio original al dominio de Fourier, donde pueden aplicarse eficazmente convoluciones espectrales con el fin de capturar las dependencias espaciales globales existentes dentro de los sistemas. Aquí, la representación del estado de entrada se proyecta en un espacio latente de mayor dimensión, y las transformaciones y activaciones se aplican de forma iterativa utilizando transformadas de Fourier, núcleos espectrales y funciones de activación. Consulte el Archivo Suplementario (S4.1) para las evaluaciones matemáticas del operador neural de Fourier.
Modelo de operador neural
El modelo de operador neuronal37 actúa como un enfoque alternativo dentro del proceso de operador neuronal para aproximar un operador físico no lineal que define la dinámica de producción de petróleo y gas. La descripción detallada del modelo de operador neuronal se encuentra en el Archivo Suplementario (S4.2) (consulte el Archivo Suplementario S10 para el Algoritmo 1: Aprendizaje Físico con Operadores Neuronales para el Algoritmo del Gemelo Digital).
Por lo tanto, el enfoque propuesto que utiliza operadores neuronales comienza con la recopilación de ejemplos emparejados de entradas y salidas, en los que las entradas son ejemplos de condiciones del sistema en forma de valores de frontera del sistema, variables de control o características físicas. El operador neuronal se entrena en una fase iterativa que consiste en predecir respuestas a una función de entrada y luego mejorar las desviaciones entre las predicciones y las respuestas reales mediante una función de error cuadrático medio. Este enfoque conduce a una optimización basada en gradientes, mediante la cual un operador neuronal descubre operadores de física sistémica que definen el comportamiento del sistema. Así, el sistema entrenado puede predecir eficientemente las respuestas del sistema a diferentes condiciones de entrada.
Modelado del sistema espacio-temporal
En la arquitectura de gemelo digital propuesta38, el proceso de producción de petróleo y gas ha sido representado mediante un grafo espacio-temporal. Esto se ha realizado para conceptualizar la conectividad, así como el proceso de evolución temporal de los recursos productivos distribuidos.
El proceso de producción de petróleo y gas tiene inherentemente una estructura distribuida en cuanto a pozos, infraestructura de tuberías y sensores. Para capturar de manera eficiente las relaciones espaciales y temporales de dicho proceso, el enfoque de gemelo digital desarrollado en este artículo utiliza redes neuronales gráficas espacio-temporales (ST-GNN). En este enfoque, las entidades del sistema se representan mediante nodos del grafo, y las relaciones espaciales o funcionales del sistema se indican mediante aristas del grafo. Las relaciones temporales se expresan mediante funciones de ventana para predecir de forma eficiente las transiciones de estado del sistema y las anomalías. Consulte el Archivo Suplementario (S5) para obtener una descripción detallada de la red neuronal gráfica espacio-temporal (ST-GNN).
Sincronización de gemelo digital mediante Neural 4D VAR
La sincronización del gemelo digital puede definirse como la sincronización continua del gemelo virtual con un entorno físico de producción de petróleo y gas en constante cambio, basada en lecturas en tiempo real de sensores. En la práctica, los gemelos digitales realizan predicciones sobre estados de un sistema, como presión, caudal y temperatura, basándose en el conocimiento de la dinámica del sistema físico. Sin embargo, debido a incertidumbres y ruido en los sensores, surge inevitablemente una discrepancia entre las predicciones y las mediciones. Por esta razón, con el fin de sincronizar las predicciones con las mediciones reales, en este trabajo se introduce otro método conocido como «asimilación de datos». Para garantizar la sincronización entre el proceso físico y el gemelo digital durante observaciones de sensores ruidosas y dispersas, se incorpora un método de asimilación de datos neural 4D-Var. El componente añadido tiene como objetivo reducir la brecha entre las mediciones de los sensores y las predicciones del modelo dentro de una ventana de asimilación. Aprovechando la diferenciabilidad, el método neural 4D-Var permite un aprendizaje de extremo a extremo y posibilita la corrección en línea del estado del gemelo digital (véase el Archivo Suplementario (S6) para la Asimilación de Datos Neural 4D-Var, el Archivo Suplementario (S10) para el Algoritmo 2: Sincronización del Gemelo Digital Usando Neural 4D-Var y el Archivo Suplementario (S11) para los hiperparámetros y detalles del entrenamiento).
Detección y diagnóstico de anomalías
En el conjunto de datos simulados de petróleo y gas, se introdujeron anomalías mediante perturbaciones controladas, como caídas de presión superiores al 15 %, picos anormales de flujo y fallos en válvulas. Los puntos de datos que superaron un valor umbral establecido o que se desviaron de una distribución estadística de referencia se identifican como anomalías. Para la formulación matemática del marco de detección de anomalías (Archivo Suplementario, S7).
Optimización y control en bucle cerrado
En el sistema propuesto, se puede lograr un control en lazo cerrado mediante la incorporación del DT armónico y el módulo de aprendizaje por refuerzo para la optimización. En la parte de los gemelos digitales, el gemelo digital predice el comportamiento futuro del proceso de producción de petróleo y gas bajo diversas condiciones, como la posición de diferentes válvulas, tasas de bombeo o control de flujo. Esto permite el desarrollo de un área de seguridad en la que el agente de control puede evaluar el resultado de las acciones deseadas sin que estas se materialicen en el proceso real. Así, aprovechando las predicciones ofrecidas por el gemelo digital, el módulo de aprendizaje por refuerzo puede optimizar objetivos operativos, incluyendo la eficiencia del proceso, el consumo de energía y la estabilidad del proceso, sin estar limitado por restricciones del proceso. Consulte el Archivo Suplementario (S8) para la descripción matemática de la optimización y control en lazo cerrado.
Aprendizaje y adaptación continuos
Los sistemas industriales enfrentan entornos no estacionarios con equipos envejecidos, condiciones operativas variables en el tiempo y perturbaciones imprevistas. Para diseñar un enfoque para esto, el modelo de gemelo digital propuesto incluye un componente de aprendizaje continuo que actualiza operadores neuronales, modelos gráficos y políticas basadas en aprendizaje incremental. Al observar variaciones importantes y patrones no identificados, el aprendizaje se facilita mediante la incorporación de nuevos datos y el ajuste de los parámetros asociados al modelo de gemelo digital.
El aprendizaje continuo también está mejorando la robustez de las políticas de control y la detección de anomalías, ya que ayuda al sistema a identificar nuevos patrones o comportamientos novedosos. La retroalimentación proveniente de las anomalías en los resultados de control se aplica para ajustar los modelos y optimizar una política de aprendizaje por refuerzo que garantice procesos de toma de decisiones eficaces a medida que el sistema sigue cambiando y evolucionando. Este proceso de aprendizaje adaptativo transforma al gemelo digital de un mero observador pasivo a un sistema inteligente capaz de mantener procesos de optimización en tiempo real de sistemas complejos de producción de petróleo y gas durante períodos prolongados, utilizando un entorno de infraestructura tecnológica riguroso y eficaz. En Tabla 3, se muestra la arquitectura del modelo y los hiperparámetros.
Consulte el Archivo Suplementario (S9) para la formulación matemática de la adaptación del aprendizaje continuo.
El módulo de preprocesamiento de datos se encarga de transformar los datos brutos de los sensores en entradas estructuradas y normalizadas adecuadas para los módulos de aprendizaje posteriores. En concreto, agrega mediciones históricas de sensores, elimina el ruido y estandariza las características para garantizar una calidad de datos consistente. Los datos procesados, junto con las condiciones de contorno y las entradas de control, se proporcionan entonces como entrada al módulo del operador neuronal con base física (modelo FNO/operador neuronal), que aprende las relaciones físicas no lineales subyacentes y genera una estimación continua del estado del sistema. Este estado del sistema estimado es utilizado posteriormente por el módulo de optimización del control para actualizar las estrategias de control de forma dinámica. Las salidas de control optimizadas se envían entonces al módulo de toma de decisiones, el cual interpreta los resultados y genera decisiones ejecutables para la operación del sistema. Se incorpora un mecanismo de retroalimentación para actualizar continuamente el sistema según nuevos datos observados, mejorando así la adaptabilidad y el rendimiento con el tiempo. Para mejorar la claridad y la reproducibilidad, se presenta una arquitectura a nivel de sistema que ilustra la interacción entre los módulos. El sistema propuesto consta de cuatro componentes principales, a saber, preprocesamiento de datos, modelado del operador neuronal con base física, optimización del control y toma de decisiones. Estos módulos están interconectados mediante interfaces de entrada-salida bien definidas, asegurando un flujo de datos continuo, un aprendizaje iterativo y una operación coordinada en todo el sistema de gemelo digital.

Figura 3: Diagrama de flujo de la arquitectura de optimización del control a nivel del sistema. Muestra la arquitectura general del sistema del trabajo propuesto. Haga clic aquí para ver una versión más grande de esta figura.
La arquitectura del sistema muestra una imagen completa del modelo que puede utilizarse para modelar y optimizar dinámicas físicas no lineales complejas en sistemas de producción de petróleo y gas. Esto comienza con los datos de entrada, que incluyen los parámetros utilizados para la simulación y las condiciones de contorno empleadas para representar el entorno físico. Estos datos primero pasan por el módulo de preprocesamiento de datos, donde se realiza la normalización y la extracción de características. Esto prepara los datos para el siguiente módulo, donde se utilizan el operador neuronal, FNO y el modelo de operador neuronal. Estos modelos se emplean para predecir los estados del sistema, y esta información se utiliza para optimizar los parámetros del sistema mediante el módulo de optimización. Los resultados optimizados se utilizan luego para tomar decisiones mediante el módulo de toma de decisiones. Este módulo ayuda a determinar la eficiencia del sistema. El marco también incluye los resultados de salida, donde se muestran los resultados del proceso de optimización en Figura 3. Además, se incluye un bucle de retroalimentación, en el cual los datos pasan nuevamente por el módulo de preprocesamiento tras los resultados de salida.