Um die Gültigkeit der physikbasierten neuronalen Operatoren zu überprüfen, wurden zusätzliche Leistungstests durchgeführt, um die Leistung von FNO und dem neuronalen Operator-Modell bei der Modellierung der Dynamik des nichtlinearen Mehrphasen-Flusssystems zu vergleichen. Die Ergebnisse zeigten, dass FNO in der globalen räumlichen Konsistenz besser abschneidet, während das neuronale Operator-Modell eine bessere Anpassungsfähigkeit an heterogene Eingaben aufweist.
Um eine effektive Schulung und Bewertung des vorgeschlagenen digitalen Zwillings-Frameworks sicherzustellen, wurden alle Experimente in einer Hochleistungsrechenumgebung durchgeführt. Als primäre Programmiersprache für die Implementierung wurde Python verwendet. Zum Erstellen von Deep-Learning-Modellen wie physikbasierten neuronalen Operatoren, räumlich-zeitlichen Graph-Neuronalen-Netzen und Verstärkungslern-Modulen kam ein gängiges Deep-Learning-Framework mit durch CUDA ermöglichter GPU-Beschleunigung sowie eine optimierte Bibliothek für tiefe neuronale Netze zum Einsatz. Um die Verarbeitung umfangreicher multivariater Zeitreihendaten und das Modelltraining zu ermöglichen, wurden die Tests auf einer Workstation mit einer Mehrkern-CPU, einer leistungsstarken Grafikverarbeitungseinheit (GPU) mit dediziertem Speicher und mindestens 64 GB Systemspeicher durchgeführt. Zur Gewährleistung der Kompatibilität mit dem Deep-Learning-Ökosystem basierte die Betriebssystemumgebung auf einer Linux-Distribution.
Die Modelle wurden mithilfe einer Mini-Batch-Gradienten-basierten Optimierung trainiert, und zahlreiche unabhängige Durchläufe der Experimente wurden durchgeführt, um die statistische Zuverlässigkeit sicherzustellen. Die vorgeschlagene Methodik kann aufgrund des bereitgestellten rechnerischen Aufbaus in vergleichbaren Hochleistungsrechnungs-Umgebungen repliziert werden.
Um die Ergebnisse durch wiederholte Versuche zu validieren, wurden alle Experimente mit mehreren unabhängigen Durchläufen und unterschiedlichen zufälligen Initialisierungswerten wiederholt. Bei Experimenten mit simulierten Daten wurde das Modell in (N) unabhängigen Versuchen ausgeführt, und die Werte wurden mit der mittleren Leistung und der Standardabweichung angegeben. Bei den Experimenten mit den Benchmark-Daten SWE und WADI wurden die Evaluierungsmetriken MAE, RMSE, F1-Score und AUROC an den Testdaten berechnet und über die wiederholten Experimente gemittelt, um die Ergebnisse statistisch robust zu machen. Dies geschah, um eine verlässliche Schätzung der Modellleistung zu erhalten und sicherzustellen, dass die Leistungsverbesserungen nicht auf ein einzelnes Experiment zurückzuführen sind.
Trotzdem haben die Ergebnisse bestätigt, dass die Kombination physikbasierter Lernverfahren, räumlich-zeitlicher Modellierung und Datenassimilationstechniken in einem integrierten digitalen Zwillingsystem unter den evaluierten Bedingungen durchgängig zu einer verbesserten Überwachungsleistung, schnellerer Konvergenz und optimalen Echtzeitentscheidungen im Vergleich zu herkömmlichen Methoden führt. Tabelle 4 zeigt die Simulationsumgebung der vorgeschlagenen Arbeit.
Die Simulationsumgebung wird so modelliert, dass sie den Echtzeitbetrieb in der Industrie genau widerspiegelt. Daten von Sensoren werden in kurzen Intervallen erfasst, um die Stream-Daten-Umgebung darzustellen, und gezielt kontrollierte Mengen an Rauschen werden hinzugefügt, um die Robustheit zu testen. Verschiedene Fehler- und Störungsfälle werden simuliert, um die Fähigkeiten zur Anomalieerkennung und Regelung zu überprüfen. Zur Bewertung wird ein Sliding-Window-Ansatz verwendet, der den Echtzeitbetrieb kontinuierlich darstellt und für dynamische Aktualisierungen durch den vorgeschlagenen digitalen Zwilling hinsichtlich Vorhersagen und Zustandssynchronisierungen sowie zur Erzeugung von Steuermaßnahmen geeignet ist.
Bei der Auslegung des steuerungsbewussten Digitalen Zwillingsmodells müssen mehrere Bewertungskriterien für eine gleichzeitige Beurteilung berücksichtigt werden. Erstens werden zur Bewertung der Überwachungs- und Zustandsvorhersagefähigkeit der Digitalen-Zwilling-Architektur Regressionskriterien wie der mittlere absolute Fehler (MAE), der Wurzelmittelquadratische Fehler (RMSE) und der mittlere absolute prozentuale Fehler (MAPE) verwendet. MAE und MAPE bewerten die Genauigkeit der neuronalen Operatoren und der raumzeitlichen Graphmodelle hinsichtlich zentraler Systemparameter wie Druck, Durchflussrate und Temperatur und sind für die kontinuierlichen multivariaten Sensordaten geeignet, die im Öl- und Gasförderprozess eingesetzt werden.
Zweitens wird auch der Synchronisierungsprozess des digitalen Zwillings selbst hinsichtlich des Zustandsschätzfehlers und der Beseitigung von Vorhersagedrift bewertet. Der Zustandsschätzfehler misst das Ausmaß der Abweichung zwischen dem tatsächlichen Zustand des Systems und dem synchronisierten Zustand des digitalen Zwillings und ist somit ein Indikator dafür, wie effektiv das entwickelte neuronale 4D-Var-Modul die Beobachtungen des tatsächlichen Systemzustands nutzt, um die Vorhersagen der Systemdynamik zu verbessern. Schließlich zeigt die Verringerung der Vorhersagedrift den Leistungsgewinn des vorgeschlagenen Ansatzes im Vergleich zu einem nicht synchronisierten digitalen Zwillingssystem auf, bedingt durch die Fähigkeit des Ansatzes, gegenüber verrauschten oder spärlichen Systembeobachtungen robust zu sein.
Drittens werden bei der Erkennung von Anomalien und der Fehlerdiagnose klassifikationsbasierte Bewertungsparameter wie Präzision, Rückruf (Recall), F1 und die Fläche unter der ROC-Kurve (AUC) verwendet. Während der Rückruf die Fähigkeit beschreibt, reale Fehler korrekt zu diagnostizieren, ohne dabei welche zu übersehen, kann die Präzision als Maß für die Gültigkeit einer erkannten Abweichung angesehen werden. Die Leistung kann angemessen anhand des F1-Kennwerts beurteilt werden. Die Bewertungsparameter werden im Allgemeinen in einem SWAT- und WADI-System angewendet.
Tabelle 5 dient der Bewertung der Wirksamkeit des vorgeschlagenen neuronalen 4D-Var-Synchronisationsmoduls hinsichtlich der Reduzierung des Vorhersagefehlers zwischen physikalischem und Cyber-Modell für drei Datensätze. Der anfängliche Zustandsfehler entspricht dem Fehler zwischen den Modellen ohne Synchronisation, während der finale Zustandsfehler dem Modellfehler nach Anwendung der Synchronisation entspricht. Es ist ersichtlich, dass sich der Zustandsfehler bei allen drei Datensätzen deutlich verbessert, wobei die größte Verbesserung (72,2 %) für das simulierte Öl- und Gassystemmodell auftritt, da es im Vergleich zu SWAT (63,6 %) und WADI (57,7 %) relativ einfach ist. Diese Ergebnisse stimmen mit den jeweiligen Komplexitätsgraden und Rauschniveaus überein. Die höhere Driftreduktionsrate beim simulierten Öl- und Gas-Datensatz lässt sich dadurch erklären, dass simulierte Umgebungen in der Regel ein geringeres Rauschniveau, weniger zufällige Störungen und ein besser definiertes Systemverhalten aufweisen als realistische cyber-physikalische Systeme wie SWAT und WADI. Diese Argumentation sollte jedoch als Hypothese betrachtet und anhand der Ergebnisse dieser Arbeit nicht quantitativ belegt werden. In zukünftigen Arbeiten könnten Komplexitätsmaße wie Entropie, Rauschvarianz und Systemdimensionalität verwendet werden, um diese Aussage zu validieren.

Abbildung 4: Überwachung der Leistung über verschiedene Datensätze hinweg (simuliertes Öl- und Gas, SWAT und WADI), bewertet anhand der Fehlerkennzahlen MAE, RMSE und MAPE (%). Die Balken stellen die aus mehreren unabhängigen Durchläufen (n = 5) ermittelten Mittelwerte dar. Fehlerbalken zeigen die Standardabweichung (SD), während schmalere überlagerte Fehlerbalken den Standardfehler des Mittelwerts (SEM = SD/√n) darstellen. Die SD- und SEM-Werte für jede Kennzahl sind explizit in der Legende angegeben. Die statistische Signifikanz wurde mittels unabhängiger t-Tests im Vergleich zum Referenzdatensatz (simuliertes Öl- und Gas) ermittelt, und die entsprechenden p-Werte sind über den Balken angezeigt. Klicken Sie bitte hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Die Leistungsanalyse des Überwachungsprozesses, der durch das vorgeschlagene Digital-Twin-Modell anhand der Datensätze Simulated Oil and Gas, SWAT und WADI unter Berücksichtigung der Bewertungskriterien MAE, RMSE und MAPE durchgeführt wurde, ist in Abbildung 4 dargestellt. Der simulierte Datensatz „Oil and Gas“ weist für alle Kriterien die geringsten Fehlerwerte auf, was deutlich die Genauigkeit der Systemzustandsüberwachung in einer simulierten Umgebung zeigt, in der die Systemparameter präzise modelliert sind. Die Datensätze SWAT und WADI zeigen aufgrund der Komplexität realer Systeme sowie deren inhärentem Rauschen und Nichtlinearitäten, die in simulierten Systemen nicht vorhanden sind, schrittweise zunehmende Fehler in der Systemüberwachung. Die geringfügige Erhöhung der Fehlermaße MAE und RMSE bestätigt jedoch die Stabilität der Systemleistung des Digital-Twin-Modells unter solchen Bedingungen.

Abbildung 5: Leistung der Anomalieerkennung über verschiedene Datensätze hinweg (simuliertes Öl und Gas, SWAT und WADI), bewertet anhand des F1-Scores und der AUC. Die Balken stellen Mittelwerte dar, die aus mehreren unabhängigen Durchläufen gewonnen wurden (n = 5). Die Fehlerbalken zeigen die SD, während dünner überlagerte Fehlerbalken das SEM = SD/√n darstellen. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
Abbildung 5 zeigt die Evaluierungsanalyse des vorgeschlagenen Systems zur Anomalieerkennung basierend auf dem F1-Score und der AUC für die drei Datensätze. Der Öl- und Gas-Datensatz, welcher der simulierte Datensatz ist, schneidet am besten ab, da er die höchsten Werte für F1-Score und AUC erreicht. Obwohl eine allmähliche Leistungsabnahme zu verzeichnen ist, erzielt das System über alle Datensätze hinweg eine hohe AUC-Leistung, mit geringfügigen Schwankungen abhängig von der Komplexität des jeweiligen Datensatzes. Die Übereinstimmung der F1-Score- und AUC-Kurven für die drei Datensätze zeigt einen optimalen Schnittpunkt, der einen guten Kompromiss zwischen den beiden Metriken aufzeigt und das vorgeschlagene System zu einem zuverlässigen System für die frühzeitige Erkennung von Fehlern und Eindringversuchen macht. Das vorgeschlagene System leistet hervorragende Ergebnisse bei der Erkennung von Anomalien in den drei Datensätzen.
Präzision, Rückruf und F1-Scores des vorgeschlagenen Digital-Twin-Modells bei der Erkennung von Anomalien in verschiedenen Datensätzen sind in Tabelle 6 unten dargestellt. Hohe Präzisionswerte zeigen zwar, dass die meisten identifizierten interessanten Punkte tatsächlich echte interessante Punkte sind, doch hohe Rückrufwerte zeigen auch, dass die meisten tatsächlichen interessanten Punkte korrekt erkannt werden, ohne wichtige zu übersehen. Die Demonstrationsumgebung im Bereich Öl und Gas weist den höchsten F1-Score von 0,94 auf, während die leicht niedrigeren Ergebnisse bei SWAT und WADI auf höhere Rauschanteile sowie stärkere Abhängigkeiten zwischen Sensoren und Angriffsmustern zurückzuführen sind. Dennoch bestätigen hohe F1-Scores, dass eine robuste räumlich-zeitliche Graphmodellierung und -synchronisation zur Anomalieerkennung erreicht wurde.

Abbildung 6: Empfänger-Charakteristik-Kurven (ROC-Kurven) zur Anomalieerkennung anhand der simulierten Datensätze für Öl und Gas, SWAT und WADI. Die Leistung wird anhand der Fläche unter der Kurve (AUC) bewertet. Die Standardabweichung (SD) und der Standardfehler des Mittelwerts (SEM = SD/√n, n = 5) der AUC-Werte sind in der Legende für jeden Datensatz angegeben. Die gestrichelte Diagonale repräsentiert die Leistung bei zufälliger Klassifizierung. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
Abbildung 6 zeigt die Receiver-Operating-Characteristic-(ROC)-Kurven zusammen mit den zugehörigen AUROC-Maßen für die drei verschiedenen Datensätze, wodurch die Genauigkeit der Anomalieerkennung unabhängig von der Wahl des Schwellenwerts bewertet wird. Bei der ROC-Kurve für den simulierten Öl- und Gas-Datensatz nähert sich der Bereich bemerkenswert dem Ursprung an, was einen AUROC-Wert von 0,99 ergibt und somit die nahezu fehlerfreie Klassifikationsfähigkeit bestätigt. Bei den SWAT- und WADI-Datensätzen liegen die AUROC-Werte bei 0,97 bzw. 0,95 und belegen damit eine gute Klassifikationsleistung, trotz höherer Unsicherheiten. Die Tatsache, dass sich die ROC-Kurven deutlich von den Linien unterscheiden, die eine zufällige Klassifikation darstellen, garantiert, dass das entwickelte Digital-Twin-Modell für die präzise Anomalieerkennung deutlich besser abschneidet als der Zufall. Die Receiver-Operating-Characteristic-(ROC)-Kurven in Abbildung 6 zeigen die Leistung des vorgeschlagenen Digital-Twin-Modells bei der Anomalieerkennung anhand der drei Datensätze. Die AUROC-Werte betragen 0,99 für den simulierten Öl- und Gas-Datensatz, 0,97 für den SWAT-Datensatz und 0,95 für den WADI-Datensatz und weisen somit auf die Diskriminierungsfähigkeit des Modells sowohl in simulierten als auch in etablierten industriellen Steuersystemen hin.
Tabelle 7 enthält eine kurze Übersicht über die kumulativen Belohnungen und Nebenbedingungsverletzungen für das Modul des Verstärkungslernens zur geschlossenen Regelung, das in den digitalen Zwilling integriert ist. Die erzielte kumulative Belohnung repräsentiert die Gesamtleistung der Regelung basierend auf den betrieblichen Zielvorgaben, während die Nebenbedingungsverletzungen die Gesamtanzahl der während des Betriebs aufgetretenen Verstöße darstellen. Die kumulativen Belohnungen bei null Nebenbedingungsverletzungen sind für unsere Modellumgebung im Bereich Öl und Gas am höchsten (96,5), was eine optimale Regelung und betriebliche Leistung belegt.
Obwohl die SWaT- und WADI-Datensätze leicht niedrigere kumulative Belohnungen und einige Verletzungen von Nebenbedingungen aufweisen, zeigt die insgesamt verbesserte Einsparung betrieblicher Kosten, dass der vorgeschlagene Digital-Twin- Rahmen effektiv für sichere und optimale Regelung in dynamischen Umgebungen bleibt.
Die Echtzeit-Eignung des vorgeschlagenen, durch Verstärkendes Lernen (RL) gesteuerten Kontrollsystems wurde durch die Analyse der Ende-zu-Ende-Inferenzlatenz aller in dem vorgeschlagenen System verwendeten Rechenmodule bewertet. Die experimentellen Beobachtungen in Tabelle 8 zeigen, dass das Neuronale-Operatoren-Modul, das FNO-/Neural-Operator-Modell und das räumlich-zeitliche Graph-Neuronale-Netzwerk-Modul (ST-GNN) aufgrund der während des Prozesses verwendeten zeitlichen Merkmale einen moderaten Rechenaufwand aufweisen. Demgegenüber hat das Modul für die Entscheidungsfindung mittels Verstärkendem Lernen einen sehr geringen Rechenaufwand. Die kumulative Ende-zu-Ende-Inferenzzeit aller in dem vorgeschlagenen System verwendeten Module liegt im Bereich typischer industrieller Steuerungszyklen. Insbesondere wurde eine durchschnittliche Ende-zu-Ende-Inferenzlatenz von 50–150 Millisekunden pro Steuerungszyklus beobachtet, was mit den Anforderungen standardmäßiger industrieller Regelungssysteme kompatibel ist. Dies zeigt, dass das vorgeschlagene System für Echtzeit- und nahezu Echtzeitanwendungen geeignet ist. Beachten Sie jedoch, dass die Ende-zu-Ende-Inferenzzeit je nach verwendeter Hardware variieren kann. In dieser Hinsicht kann das vorgeschlagene System für Echtzeitanwendungen eingesetzt werden.
Ein quantitativer Vergleich der FNO-, Neural-Operator-Modell- und Baseline-Ansätze, wie traditionelle numerische Löser und aktuelle Operator-Lernverfahren, wird durchgeführt, um die Wirksamkeit physikbasierter neuronaler Operatoren weiter zu untermauern. Die Bewertung anhand von Tabelle 9 konzentriert sich darauf, wie gut jedes Modell die dynamischen und nichtlinearen Systemdynamiken erfasst. Zur Leistungsbeurteilung werden Standardmetriken wie mittlerer quadratischer Fehler (MSE), mittlerer absoluter Fehler (MAE) und relativer Fehler zwischen erwarteten und realen Systemzuständen herangezogen. Laut experimenteller Ergebnisse übertrumpfen sowohl das FNO- als auch das Neural-Operator-Modell die Baseline-Ansätze bei der Modellierung komplexer Mehrphasendynamiken, da sie niedrigere Fehlerwerte erreichen. Insbesondere erfasst das Neural-Operator-Modell nichtlineare Operatormapping unter verschiedenen Eingangsbedingungen erfolgreich, während FNO gut darin ist, globale räumliche Abhängigkeiten zu lernen.
Die Wirksamkeit der vorgeschlagenen Neural-Operator-Modelle wird durch die quantitativen Vergleichsergebnisse in Tabelle 9 belegt. Im Vergleich zu basislinearen numerischen Lösern und aktuellen Operator-Lernverfahren erreichen sowohl FNO als auch das Neural-Operator-Modell deutlich niedrigere Fehlermetriken. Insbesondere weist FNO den geringsten relativen Fehler und den niedrigsten mittleren quadratischen Fehler auf, was seine überlegene Fähigkeit zeigt, die Dynamik komplexer nichtlinearer Systeme abzubilden. Außerdem zeigt das Neural-Operator-Modell signifikante Fortschritte, insbesondere beim Erlernen nichtlinearer Operator-Abbildungen unter verschiedenen Bedingungen. Zur Bewertung der statistischen Signifikanz werden unabhängige t-Tests verwendet; p-Werte unter 0,05 deuten auf eine wesentliche Verbesserung gegenüber dem basislinearen numerischen Löser hin. Die von beiden vorgeschlagenen Modellen erzielten p-Werte liegen unter 0,001 und zeigen damit, dass die beobachteten Leistungssteigerungen statistisch signifikant sind und nicht zufällig entstanden. Darüber hinaus weisen die vorgeschlagenen Verfahren deutlich kürzere Inferenzzeiten auf, was ihre Eignung für Echtzeit-Digital-Twin-Anwendungen unterstreicht.
Um die Klarheit der experimentellen Bewertung des vorgeschlagenen steuerungsbewussten Digitalen Zwillingsmodells weiter zu verbessern, wird ein umfassender quantitativer Vergleich des vorgeschlagenen Ansatzes mit anderen Methoden in einem integrierten tabellarischen Format bereitgestellt. Dieser Vergleich verschiedener Leistungsparameter wie mittlerer absoluter Fehler (MAE), Wurzel des mittleren quadratischen Fehlers (RMSE), F1-Score, Fläche unter der Kurve (AUC), Reduktion der Vorhersageabweichung und kumulierte Belohnung ermöglicht eine integrierte Bewertung der Überwachungsgenauigkeit, der Anomalieerkennungsfähigkeit, der Synchronisationseffizienz und der Regelungsleistung der Methoden in Tabelle 10. Zusätzlich werden verschiedene Parameter des Trainingsverfahrens wie Batch-Größe, Lernrate, Anzahl der Epochen und modulweise Trainingskonfiguration explizit genannt, um die Transparenz der experimentellen Bewertung in Tabelle 11 zu erhöhen. Darüber hinaus werden Skalierbarkeit und Recheneffizienz hinsichtlich Inferenz-Latenz, GPU-Auslastung und Echtzeitfähigkeit des vorgeschlagenen Ansatzes in Tabelle 12 evaluiert.
Der vorgeschlagene digitale Zwilling kann umfassend mit anderen Methoden wie der physikbasierten Synchronisierungsmethode, LSTM- und GNN-Modellen hinsichtlich Überwachungsgenauigkeit, Synchronisationsfähigkeit, Anomalieerkennung und Regelungsleistung verglichen werden, wie in Tabelle 13 dargestellt. Physikbasierte digitale Zwillinge weisen aufgrund statischer Ansätze eine geringere Anpassungsfähigkeit und eine höhere Fehlerquote auf. Neuronale Ansätze wie LSTM und GNN sind besser, kombinieren Physik, Synchronisation und Regelung jedoch nicht effektiv. Der vorgeschlagene digitale Zwilling erreichte die besten Werte mit einem RMSE von 0,20, einer Driftreduktion von 72,2 %, einem AUROC von 0,94, einem F1-Score von 0,92 und einer kumulativen Belohnung von 95,5. Die oben genannten Ergebnisse zeigen deutlich, dass die effektive Kombination physikbasierter neuronaler Operatoren, graphenbasierter Modellierung von Raum und Zeit, neuronaler Datenassimilation und Verstärkendes Lernen ein genaueres, robustes und regelungsbewusstes System für digitale Zwillinge erheblich verbessern kann.

Abbildung 7: Trainingskosten und Ressourcennutzung der verschiedenen Komponenten des vorgeschlagenen Systems, einschließlich FNO, neuronales Operatormodell, GNN, RL-Modul und des gesamten Frameworks. Die Metriken umfassen die Trainingsdauer (Stunden) und die GPU-Speicherbelegung (GB). Die Werte stellen Mittelwerte aus mehreren unabhängigen Durchläufen dar (n = 5). Die Fehlerbalken zeigen die SD, während überlagerte, dünnere Fehlerbalken das SEM = SD/√n darstellen und somit die Variabilität zwischen den Durchläufen widerspiegeln. Klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
In Abbildung 7 sind die Trainingskosten und die Ressourcennutzung des vorgeschlagenen Ansatzes in ihren Hauptbestandteilen dargestellt. Es ist ersichtlich, dass neuronale Operatormodelle wie FNO und das neuronale Operatormodell im Vergleich zum Gesamtsystem erhebliche Trainingskosten und eine hohe GPU-Speichernutzung aufweisen. Die hohen Trainingskosten ergeben sich aus der Fähigkeit der Operatoren, hochdimensionale Zusammenhänge zu lernen. Es zeigt sich, dass die graphbasierten neuronalen Netzwerke und die Module des Verstärkungslernens geringere Trainingskosten verursachen, jedoch wesentlich für das Erlernen räumlicher Abhängigkeiten und einer adaptiven Steuerung sind. Obwohl das vorgeschlagene System höhere Trainingskosten verursacht, erweist sich die Inferenzzeit als effizient. Im Vergleich zu Baseline-Methoden wie eigenständigen LSTM- oder GNN-Modellen entstehen durch die Integration mehrerer Module höhere Trainingskosten, doch wird eine verbesserte Genauigkeit und Robustheit erreicht, wie in Tabelle 13 gezeigt.
Insgesamt kombiniert die vorgeschlagene Digital-Twin-Struktur mehrere Elemente, von denen jedes ein spezifisches Problem bei der Modellierung und Steuerung komplexer industrieller Systeme angeht. Im Gegensatz zu herkömmlichen numerischen Lösern, die häufig rechenaufwendig sind und weniger für Echtzeitanwendungen geeignet, werden physikbasierte neuronale Operatoren wie der Fourier-Neuronale-Operator und das Deep-Operator-Netzwerk eingesetzt, um nichtlineare und hochdimensionale Systemdynamiken effizienter abzubilden. Einfachere maschinelle Lernmodelle, wie Vorwärts- oder rekurrente neuronale Netze, können verwendet werden, lernen jedoch meist keine grundlegenden physikalischen Operatoren und können nicht auf verschiedene Systemkonfigurationen verallgemeinern. Abhängigkeiten zwischen verteilten Sensoren, die in Öl- und Gasförderanlagen intrinsisch verknüpft sind, werden mittels eines räumlich-zeitlichen graphbasierten neuronalen Netzwerks modelliert. Konventionelle Ansätze wie unabhängige Zeitreihenmodelle oder einfache rekurrente Architekturen können zu geringerer Genauigkeit bei der Zustandsschätzung und Anomalieerkennung führen, da sie räumliche Verbindungen nicht explizit berücksichtigen. Die adaptive Regelung, die das Modul des Verstärkungslernens nutzt, ermöglicht es dem System, betriebliche Entscheidungen in dynamischen und unsicheren Umgebungen zu optimieren. Regelbasierte oder statische Regelstrategien hingegen sind starr und nicht in der Lage, sich an veränderte Systembedingungen oder unvorhergesehene Störungen anzupassen. Obwohl es Alternativen mit geringerer Komplexität gibt, befassen diese sich typischerweise nur mit bestimmten Teilaspekten des Problems und reichen nicht aus, um eine integrierte Lösung zu erreichen, die präzise Modellierung, Echtzeitüberwachung und adaptive Regelung vereint. Durch die Sicherstellung, dass jeder Bestandteil zur Gesamtsystemleistung beiträgt, verbessert das vorgeschlagene modulare Design die Genauigkeit, Skalierbarkeit und Robustheit in anspruchsvollen industriellen Umgebungen.
| Studie / Arbeit | Verwendete Techniken | Ergebnisse | Einschränkungen |
| Tiefer Neuronaler Operator für Digital Twin26 | Tiefer Neuronaler Operator DeepONet als Ersatzmodell | Echtzeitvorhersage; um Größenordnungen schneller als physikbasierte Simulationen | Benötigt sorgfältige Sensorplatzierung; Evaluationsmethoden müssen verbessert werden |
| Digitales Zwillingsmodell mit virtueller Sensorik basierend auf DeepONet27 | Virtuelle Sensoren mittels DeepONet im Digitalen Zwilling | Schnelle Vorhersagen; Echtzeit-Inferenz; umfassende Zustandsschätzung | Spektrale Verzerrung beeinträchtigt hochfrequente Phänomene; hybride Modellierung erforderlich |
| Digitaler Zwilling mit DRL + NMPC in der Regelung28 | Integration von Deep Reinforcement Learning und nichtlinearer modellprädiktiver Regelung (NMPC) | Verbesserte Regelungsleistung; geringere Regelabweichung und adaptive Aktualisierung des Zwillings | Anwendung bisher auf Schiffe beschränkt; Erweiterung auf Fördersysteme notwendig |
| KI-gestützter Digitaler Zwilling in der industriellen Automatisierung29 | Überwachtes/Unüberwachtes maschinelles Lernen (LSTM, CNN, RF, etc.) | Hohe Genauigkeit bei Fehlerprognose und -überwachung | Überwiegend klassische ML-Methoden; fortschrittliche räumlich-zeitliche Deep-Learning-Modelle fehlen |
| Übersichtsarbeit zum Digitalen Zwilling in der Öl- und Gasindustrie16 | Verschiedene Ansätze aus KI/ML und mechanistischer Modellierung | Umfassende Taxonomie der Modellierungsstrategien | Überwiegend konzeptionell; kaum praktische Implementierungen mit fortschrittlichem Deep Learning |
Tabelle 1: Vergleichende Analyse der Öl- und Gasförderung und -steuerung im Vergleich zu bestehenden Methoden. Ein Überblick über relevante Digital-Twin- und Steuerungsstrategien unter Betonung der eingesetzten Methoden, erzielten Ergebnisse und wesentlichen Nachteile.
| Datensatzname | Domäne / System | Anzahl der Merkmale | Dauer und Größe | Datencharakteristiken |
| Simulierte Öl- und Gasförderungsdaten | Öl- und Gasförderungssystem (Simulation) | Multivariat (Druck, Durchfluss, Temperatur usw.) | Für diese Studie konzipiert | Simulierte Echtzeitsensor- und Steuersignale, normale und Fehlerzustände |
| SWaT | Sichere Wasserbehandlungs-CPS | ~51 Sensoren + Aktoren | ~11 Tage (~450.000+ Messwerte) | Normal- und Anomaliezustände (36 Angriffsszenarien) |
| WADI | Wasser-Versorgungs-CPS | ~123 Sensoren + Aktoren | ~16 Tage (~1 Mio.+ Messwerte) | Normal- und Anomaliezustände (15 Angriffe) |
Tabelle 2: Datensatzbeschreibung.
Beschreibung des Datensatzes. Zusammenfassung der Datensätze der Studie, einschließlich Domäne, Anzahl der Merkmale, Dauer und Datenattribute.
| Komponente | Hyperparameter | Wert |
| ST-GNN | Anzahl der Schichten | 3 |
| ST-GNN | Versteckte Dimension | 64 |
| ST-GNN | Länge des zeitlichen Fensters | 12 |
| Neuronaler Operator (FNO) | Anzahl der Fourier-Modi | 16 |
| Neuronaler Operator (FNO) | Anzahl der Schichten | 4 |
| DeepONet | Größe des Branch-Netzwerks | [128, 128] |
| DeepONet | Größe des Trunk-Netzwerks | [128, 128] |
| Datenassimilation | Assimilationsfenster | 10 Zeitschritte |
| Optimierer | Adam | |
| Lernrate | 0,001 | |
| Batch-Größe | 64 | |
| RL-Regler | Abzinsungsfaktor (γ) | 0,99 |
| RL-Regler | Anzahl der Schichten im Policy-Netzwerk | [128, 64] |
Tabelle 3: Modellarchitektur und Hyperparameter. Wichtige architektonische Konfigurationen und Hyperparameter-Einstellungen für die neuronalen Operatoren, die Datenassimilation, das Verstärkungslernen sowie die Komponenten des räumlich-zeitlichen Graph-Neuronalen-Netzwerks.
| Parameter | Beschreibung |
| Simulationsplattform | Python (v3.10) mit PyTorch (v2.1) |
| Unterstützende Bibliotheken | NumPy (v1.24), SciPy (v1.10), Pandas (v1.5) |
| Hardware | NVIDIA RTX-Serie GPU (z. B. RTX 3090, 24 GB VRAM), Intel Core i7/i9 CPU, 32–128 GB RAM |
| GPU-Beschleunigung | CUDA Toolkit (v11.8), cuDNN (v8.x) |
| Betriebssystem | Linux (Ubuntu 20.04) / Windows 10 |
| Entwicklungsumgebung | Jupyter Notebook / VS Code |
| Abtastintervall | 1–10 Sekunden (konfigurierbar für Echtzeitsimulation) |
| Modellierung von Sensorrauschen | Gaußsches Rauschen (σ = 0,01–0,05) |
| Fehlerszenarien | Sensorfehler, Strömungsstörungen und Steuerungsanomalien |
| Evaluierungsmodus | Echtzeitsimulation basierend auf gleitendem Fenster |
| Verfügbarkeit des Codes | Repository-Link wird auf Anfrage / nach Veröffentlichung zur Reproduzierbarkeit bereitgestellt |
Tabelle 4: Simulationsumgebung. Informationen zur Hardwarekonfiguration, Softwareplattform, Abtastintervall, Rauschmodellierung und Fehler-Szenarien im experimentellen Aufbau.
| Datensatz | Anfangsfehler | Endfehler | Driftreduzierung (%) |
| Simuliertes Öl & Gas | 0.9 | 0.25 | 72.2 |
| SWAT | 1.1 | 0.4 | 63.6 |
| WADI | 1.3 | 0.55 | 57.7 |
Tabelle 5: Experimentelles Ergebnis der Digital Twin-Synchronisation. Quantitative Bewertung der Synchronisationsleistung, die die Verringerung des Vorhersagefehlers über verschiedene Datensätze hinweg sowie die Schätzfehler des Anfangs- und Endzustands zeigt.
| Datensatz | Präzision | Recall | F1-Score |
| Simuliertes Öl- und Gasfeld | 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 |
Tabelle 6: Leistung der Anomalieerkennung. Ergebnisse der Anomalieerkennung an simulierten Öl- und Gasdatensätzen sowie den SWAT- und WADI-Datensätzen, einschließlich Präzision, Sensitivität und F1-Score.
| Datensatz | Kumulierte Belohnung ↑ | Nebenbedingungsverletzungen ↓ | Kostenreduktion (%) ↑ |
| Simuliertes Öl- und Gasfeld | 96.5 ± 0.8 | 0 | 18.7 |
| SWAT | 92.1 ± 1.2 | 1 | 14.3 |
| WADI | 89.7 ± 1.5 | 2 | 11.6 |
Tabelle 7: Experimentelle Ergebnisse der geschlossenen Regelung und Optimierung. Kumulierte Belohnung, Verletzungen von Nebenbedingungen und Reduktion der Betriebskosten werden zur Bewertung der Leistung der Regelung im geschlossenen Kreis verwendet.
| Modul | Durchschnittliche Latenz (ms) | Beschreibung |
| Neuraler Operator (FNO/DeepONet) | 18-25 ms | Lernen von Operatoren und Vorhersage des Systemzustands |
| ST-GNN | 22–35 ms | Modellierung räumlich-zeitlicher Abhängigkeiten |
| RL-Richtliniennetzwerk | 5–10 ms | Erzeugung von Steuerungsentscheidungen |
| Datenvorverarbeitung | 8–12 ms | Eingangsnormalisierung und Merkmalsvorbereitung |
| Gesamtlatenz der Regelungsschleife | 55–82 ms | Ausführungszeit für die End-to-End-Steuerungsaktion |
Tabelle 8: Bewertung der Echtzeit-Regelungsleistung. Die Latenzwerte stellen die durchschnittlichen Inferenzzeiten dar, die pro Regelungsschritt unter einer standardmäßigen Hardwarekonfiguration gemessen wurden.
| Modell | MSE ↓ | MAE ↓ | Relativer Fehler (%) ↓ | Inferenzzeit (ms) ↓ | p-Wert (gegenüber Baseline) |
| Numerischer Löser (Baseline) | 0.0125 ± 0.0012 | 0.089 ± 0.006 | 8.75 ± 0.54 | 120.5 ± 5.2 | — |
| Konventionelles ML-Modell | 0.0098 ± 0.0009 | 0.072 ± 0.005 | 6.42 ± 0.48 | 85.3 ± 4.7 | 0.021 |
| Bekanntes Operator-Lernen | 0.0076 ± 0.0007 | 0.061 ± 0.004 | 5.18 ± 0.41 | 42.8 ± 3.1 | 0.008 |
| DeepONet (vorgeschlagen) | 0.0052 ± 0.0005 | 0.044 ± 0.003 | 3.67 ± 0.32 | 18.6 ± 2.4 | < 0.001 |
| FNO (vorgeschlagen) | 0.0047 ± 0.0004 | 0.039 ± 0.002 | 3.21 ± 0.28 | 15.2 ± 2.1 | < 0.001 |
Tabelle 9: Vergleich der Baseline- und der vorgeschlagenen Modelle hinsichtlich Vorhersagegenauigkeit (MSE, MAE, relativer Fehler) und rechnerischer Effizienz (Inferenzzeit). Die Werte sind als Mittelwert ± Standardabweichung über mehrere Durchläufe angegeben. Niedrigere Werte weisen auf eine bessere Leistung hin. Die p-Werte bezeichnen die statistische Signifikanz im Vergleich zum Baseline-Numerik-Löser.
| Methode | MAE ↓ | RMSE ↓ | F1-Score ↑ | AUC ↑ | Drift-Reduktion (%) ↑ | Kumulierte Belohnung ↑ |
| Physikbasiertes DT | 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 |
| Vorgeschlagener Ansatz | 0.2 | 0.3 | 0.94 | 0.97 | 72.2 | 95.5 |
Tabelle 10: Vergleich der umfassenden Leistung. Vergleich der umfassenden Leistung des vorgeschlagenen steuerungsbezogenen Digital-Twin-Modells mit anderen Basisverfahren hinsichtlich Überwachung, Anomalieerkennung, Synchronisation und Steuerung. Bei der Vorhersagegenauigkeit werden niedrigere Werte von MAE und RMSE bevorzugt. Bei der Leistung von Anomalieerkennung und Synchronisation werden höhere Werte von F1-Score, AUC, Drift-Reduktion und Belohnung bevorzugt.
| Komponente | Parameter | Wert/Einstellung |
| Neuronale Operatoren (FNO/DeepONet) | Epochen | 100–200 |
| Batch-Größe | 32–64 |
| Lernrate | 0,001 |
| ST-GNN | Schichten | 3 |
| Versteckte Einheiten | 64 |
| Zeitliches Fenster | 10–20 Zeitschritte |
| Neurales 4D-Var | Assimilationsfenster | Fest (gleitendes Fenster) |
| Optimierung | Gradientabstieg |
| Verstärkendes Lernen | Abzinsungsfaktor (γ) | 0,95 |
| Bewertungsgewichte (α,β,γ) | 0,6, 0,25, 0,15 |
| Allgemeines Training | Optimierer | Adam |
| Hardware | NVIDIA RTX GPU |
Tabelle 11: Trainingskonfiguration und Implementierungsparameter, die für verschiedene Komponenten des vorgeschlagenen Digital-Twin-Frameworks verwendet wurden. Die Tabelle bietet einen Überblick über verschiedene Hyperparameter und Trainingsparameter, die verwendet wurden, um die Reproduzierbarkeit des vorgeschlagenen Systems für verschiedene Komponenten sicherzustellen.
| Komponente | Trainingszeit (Stunden) | Inferenz-Latenz (ms) | GPU-Auslastung (%) | Skalierbarkeit |
| Neural Operators | 6–8 | 25–40 | 80–90 | Hoch |
| ST-GNN | 3–4 | 15–25 | 60–70 | Hoch |
| Neural 4D-Var | 2–3 | 20–30 | 65–75 | Mäßig |
| RL-Regler | 1–2 | 10–20 | 50–60 | Hoch |
| Gesamtsystem | 12–16 | <100 ms | 70–85 | Skalierbar |
Tabelle 12: Rechenleistung und Skalierbarkeitsanalyse des vorgeschlagenen Digital-Twin-Modells für verschiedene Komponenten. Die Latenz steht für die Echtzeit-Inferenzleistung, GPU für die Rechenleistung, und Skalierbarkeit beschreibt die Fähigkeit des vorgeschlagenen Modells, in industriellen Anwendungen im großen Maßstab eingesetzt zu werden.
| Methode | RMSE-Überwachung | Synchronisationsdrift-Reduktion (%) | AUROC | F1-Score | Kumulativer Reward |
| Physikbasiertes digitales Zwillingsmodell | 0.45 | 28.6 | 0.78 | 0.74 | 42.3 |
| LSTM-basiertes Modell | 0.31 | 41.2 | 0.86 | 0.83 | 61.7 |
| GNN-basiertes Modell | 0.26 | 53.4 | 0.91 | 0.88 | 74.5 |
| Vorgeschlagenes digitales Zwillingsmodell | 0.2 | 72.2 | 0.94 | 0.92 | 95.5 |
Tabelle 13: Vergleichende Bewertung der vorgeschlagenen Arbeit mit bestehenden Methoden. Die Leistung des vorgeschlagenen Digital Twin-Modells hinsichtlich Überwachung, Synchronisation, Anomalieerkennung und Steuerungsmetriken wird mit Basisverfahren und lernbasierten Ansätzen verglichen.
| Modellkonfiguration | Überwachungsfehler (RMSE) | Anomalieerkennung (F1) | Regelleistung (Belohnung) |
| Vollständiges Modell (alle Module) | 0.2 | 0.92 | 95.5 |
| – Ohne Datenassimilation | 0.29 | 0.88 | 89.2 |
| – Ohne ST-GNN | 0.27 | 0.85 | 87.1 |
| – Ohne Verstärkendes Lernen | 0.21 | 0.9 | 78.4 |
| – Ohne kontinuierliches Lernen | 0.23 | 0.91 | 92.3 |
Tabelle 14: Experimentelle Ergebnisse und Ablationsstudie. Ablationsstudie zur Untersuchung, wie jeder Teil des vorgeschlagenen Digital-Twin-Modells zur Gesamtleistung des Systems beiträgt.
| Szenario | Datenverfügbarkeit | RMSE | Genauigkeit | Behandlung von Beobachtungslücken |
| Normal | 100% | Niedrig | Hoch | Stabil |
| Teilweiser Verlust | 50–70% | Mäßig | Hoch | Robuste Interpolation |
| Starker Verlust | 20–30% | Erhöht | Mäßig | Sanfte Degradation |
| Sensorausfall | Regional 0% | Höher | Mäßig | Räumliche Generalisierung |
Tabelle 15: Robustheitsanalyse unter Datenknappheit. Die folgende Tabelle zeigt die Leistung des vorgeschlagenen neuronalen 4D-Var-Verfahrens unter verschiedenen Bedingungen der Datenverfügbarkeit. Die Tabelle veranschaulicht die Robustheit des Modells, indem sie darstellt, wie sich die Fehler bei teilweisem Verlust, hoher Sparsamkeit und vollständigem Verlust von Sensordaten verändern.