$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
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 Operatormodell 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 Operatormodell 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 kam Python zum Einsatz. Mit einem verbreiteten Deep-Learning-Framework, das durch CUDA GPU-Beschleunigung ermöglicht und eine optimierte Bibliothek für tiefe neuronale Netze nutzt, wurden Deep-Learning-Modelle wie physikbasierte neuronale Operatoren, räumlich-zeitliche graphenbasierte neuronale Netze und Verstärkungslern-Module erstellt. Um die Verarbeitung umfangreicher multivariater Zeitreihendaten und das Modelltraining zu erleichtern, 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.
Zur Schulung der Modelle wurde eine auf Mini-Batches basierende Gradientenoptimierung verwendet, 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 Setups in vergleichbaren Hochleistungsrechnenumgebungen 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 für (N) unabhängige Versuche ausgeführt, und die Werte werden mit mittlerer Leistung und Standardabweichung angegeben. Bei den Experimenten mit den Benchmark-Datensätzen SWAT 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 zuverlässige Schätzung der Modellleistung zu erhalten und sicherzustellen, dass 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, um den Echtzeitbetrieb kontinuierlich darzustellen, was dynamische Aktualisierungen durch den vorgeschlagenen digitalen Zwilling für Vorhersagen und Zustandssynchronisierungen sowie für die Erzeugung von Steuermaßnahmen ermöglicht.
Bei der Konzeption des steuerungsbewussten digitalen Zwillingsmodells müssen mehrere Bewertungskriterien für eine gleichzeitige Beurteilung berücksichtigt werden. Zunächst werden zur Bewertung der Überwachungs- und Zustandsvorhersagefähigkeit der digitalen Zwillingsarchitektur Regressionskriterien wie der mittlere absolute Fehler (MAE), der quadratische Mittelwert des mittleren quadratischen Fehlers (RMSE) und der mittlere absolute prozentuale Fehler (MAPE) verwendet. MAE und MAPE berechnen die Genauigkeit neuronaler Operatoren und räumlich-zeitlicher Graphenmodelle hinsichtlich zentraler Systemparameter wie Druck, Durchflussrate und Temperatur und gelten für die kontinuierlichen multivariaten Daten der Sensoren, 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 Systemzustand 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 zum 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, echte Fehler korrekt zu diagnostizieren, ohne 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-Werts beurteilt werden. Die Bewertungsparameter werden im Allgemeinen in einem SWAT- und WADI-System angewandt.
Tabelle 5 dient zur Beurteilung der Wirksamkeit des vorgeschlagenen neuronalen 4D-Var-Synchronisationsmoduls hinsichtlich der Verringerung des Vorhersagedrifts zwischen dem physikalischen und dem Cyber-Modell für drei Datensätze. Der Anfangszustandsfehler entspricht dem Fehler zwischen den Modellen ohne Synchronisation, und der Endzustandsfehler entspricht dem Modellfehler nach Anwendung der Synchronisation. 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 Werte stimmen mit den jeweiligen Komplexitätsgraden und Rauschniveaus überein. Die höhere Driftreduktionsrate beim simulierten Öl- und Gasdatensatz 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 realistischere 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 Argumentation zu validieren.

Abbildung 4: Überwachung der Leistung über verschiedene Datensätze hinweg (simuliertes Öl und Gas, SWAT und WADI), bewertet anhand der Fehlermetriken MAE, RMSE und MAPE (%). Die Balken stellen die aus mehreren unabhängigen Durchläufen (n = 5) ermittelten Mittelwerte dar. Die Fehlerbalken zeigen die Standardabweichung (SD), während dünnere überlagerte Fehlerbalken den Standardfehler des Mittelwerts (SEM = SD/√n) darstellen. Die SD- und SEM-Werte für jede Metrik sind explizit in der Legende angegeben. Die statistische Signifikanz wurde mittels unabhängiger t-Tests im Vergleich zum Baseline-Datensatz (simuliertes Öl und Gas) ermittelt, und die entsprechenden p-Werte sind über den Balken dargestellt. Bitte klicken Sie hier, um eine größere 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 Gasnetz, SWAT und WADI), bewertet anhand des F1-Scores und der AUC. Die Balken zeigen die Mittelwerte aus mehreren unabhängigen Durchläufen (n = 5). Die Fehlerbalken geben die SD an, während dünnere ü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, weist die beste Leistung auf, 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 offenbart und das vorgeschlagene System zu einem zuverlässigen System für die frühzeitige Erkennung von Fehlern und Eindringversuchen macht. Das vorgeschlagene System zeichnet sich durch eine hervorragende Leistung bei der Erkennung von Anomalien in den drei Datensätzen aus.
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, während hohe Rückrufwerte darauf hinweisen, 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: ROC-Kurven (Receiver Operating Characteristic) zur Anomalieerkennung anhand der Datensätze Simulated Oil and 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 Linie stellt die Leistung einer zufälligen Klassifizierung dar. Bitte klicken Sie hier, um eine größere 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 zufälliger Klassifikation entfernen, garantiert, dass das entwickelte Digital-Twin-Modell für die präzise Anomalieerkennung deutlich besser abschneidet als der Zufallslevel. 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 Verstärkungslernmodul zur geschlossenen Regelung, das in den digitalen Zwilling integriert ist. Die erreichte kumulative Belohnung repräsentiert die Gesamtleistung der Regelung basierend auf den betrieblichen Zielvorgaben, während die Nebenbedingungsverletzungen die Gesamtanzahl der Verstöße während des Betriebs darstellen. Die kumulativen Belohnungen bei null Nebenbedingungsverletzungen sind für unsere Modellumgebung im Öl- und Gassektor 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 das vorgeschlagene Digital-Twin-Framework auch in dynamischen Umgebungen wirksam für sichere und optimale Steuerung bleibt.
Die Echtzeit-Eignung des vorgeschlagenen, durch Verstärkungslernen (RL) gesteuerten Regelungssystems 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 Verstärkungslernen-Entscheidungsfindung einen sehr geringen Rechenaufwand. Die kumulative Ende-zu-Ende-Inferenzzeit aller in dem vorgeschlagenen System verwendeten Module liegt im Bereich typischer industrieller Regelzyklen. Insbesondere wurde eine durchschnittliche Ende-zu-Ende-Inferenzlatenz von 50–150 Millisekunden pro Regelzyklus beobachtet, was mit den Anforderungen standardmäßiger industrieller Regeltechnik 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 physikgeleiteter neuronaler Operatoren weiter zu untermauern. Die Bewertung aus Tabelle 9 konzentriert sich darauf, wie gut jedes Modell die dynamischen und nichtlinearen Dynamiken des Systems erfasst. Zur Leistungsbeurteilung werden Standardmetriken wie mittlerer quadratischer Fehler (MSE), mittlerer absoluter Fehler (MAE) und relativer Fehler zwischen erwarteten und realen Systemzuständen verwendet. Laut experimentellen Ergebnissen übertreffen sowohl FNO als auch das Neural-Operator-Modell die Baseline-Ansätze bei der Modellierung komplexer Mehrphasendynamiken, da sie geringere 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 neuronalen Operatormodelle wird durch die quantitativen Vergleichsergebnisse in Tabelle 9 belegt. Im Vergleich zu Baseline-Numeriklösern und aktuellen Methoden des Operatoren-Lernens erreichen sowohl FNO als auch das neuronale Operatormodell 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. Zudem zeigt das neuronale Operatormodell signifikante Verbesserungen, insbesondere beim Erlernen nichtlinearer Operatorenabbildungen unter verschiedenen Bedingungen. Unabhängige t-Tests werden verwendet, um die statistische Signifikanz zu bewerten; p-Werte unter 0,05 deuten auf eine beträchtliche Verbesserung gegenüber dem Baseline-Numeriklöser hin. Die von beiden vorgeschlagenen Modellen erzielten p-Werte liegen unter 0,001, was darauf hindeutet, dass die beobachteten Leistungssteigerungen statistisch signifikant sind und nicht zufällig zustande gekommen sind. 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 Vorhersagedrift und kumulierter Belohnung ermöglicht eine integrierte Bewertung der Überwachungsgenauigkeit, Anomalieerkennungsfähigkeit, Synchronisationseffizienz und Steuerungsleistung 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 verbessern. Darüber hinaus werden Skalierbarkeit und rechnerische Effizienz hinsichtlich Inferenzlatenz, GPU-Auslastung und Echtzeitfähigkeit des vorgeschlagenen Ansatzes in Tabelle 12 bewertet.
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 sind aufgrund statischer Ansätze weniger anpassungsfähig und weisen 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 hat die besten Werte mit einem RMSE von 0,20, einer Driftreduzierung von 72,2 %, einem AUROC von 0,94, einem F1-Score von 0,92 und einer kumulativen Belohnung von 95,5 erreicht. Die oben genannten Ergebnisse zeigen eindeutig, dass die effektive Kombination physikbasierter neuronaler Operatoren, graphbasierter Modellierung von Raum und Zeit, neuronaler Datenassimilation und verstärkenden Lernens ein genaues, robustes und regelungsbewusstes System für digitale Zwillinge erheblich verbessern kann.

Abbildung 7: Trainingskosten und Ressourcennutzung verschiedener Komponenten des vorgeschlagenen Systems, einschließlich FNO, neuronales Operatormodell, GNN, RL-Modul und des gesamten Frameworks. Zu den Metriken gehören die Trainingszeit (Stunden) und die GPU-Speichernutzung (GB). Die Werte stellen Mittelwerte aus mehreren unabhängigen Durchläufen dar (n = 5). Fehlerbalken zeigen die SD an, während dünnere überlagerte Fehlerbalken das SEM = SD/√n darstellen und die Variabilität über die Durchläufe widerspiegeln. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
In Abbildung 7 sind die Trainingskosten und die Ressourcennutzung des vorgeschlagenen Ansatzes in ihren Hauptkomponenten 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 Graph-Neural-Network- und die Verstärkendes-Lernen-Module geringere Trainingskosten haben, 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 das vorgeschlagene System erreicht eine verbesserte Genauigkeit und Robustheit, 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 adressiert. Im Gegensatz zu herkömmlichen numerischen Lösungsverfahren, die häufig rechenintensiv 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 zwar verwendet werden, lernen jedoch in der Regel keine zugrundeliegenden physikalischen Operatoren und können nicht auf verschiedene Systemkonfigurationen verallgemeinern. Abhängigkeiten zwischen verteilten Sensoren, die in Öl- und Gasförderanlagen inhärent miteinander verknüpft sind, werden mithilfe des räumlich-zeitlichen Graph-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 Steuerung, die auf dem Verstärkungslernmodul basiert, ermöglicht es dem System, betriebliche Entscheidungen in dynamischen und unsicheren Umgebungen zu optimieren. Regelbasierte oder statische Steuerungsstrategien hingegen sind starr und nicht in der Lage, sich an wechselnde Systembedingungen oder unvorhergesehene Störungen anzupassen. Obwohl es Alternativen mit geringerer Komplexität gibt, behandeln diese meist nur bestimmte Aspekte des Problems und reichen nicht aus, um eine integrierte Lösung zu erreichen, die präzise Modellierung, Echtzeitüberwachung und adaptive Steuerung vereint. Durch die sicher gestellte Beitragsleistung jedes einzelnen Elements zur Gesamtsystemleistung 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 Zwillingssystem mit virtueller Sensorik basierend auf DeepONet27 | Virtuelle Sensoren nach 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 & nichtlinearem MPC | Verbesserte Regelungsleistung; geringere Regeldifferenz und adaptive Anpassung des Zwillings | Anwendung auf Schiffe beschränkt; weitere Anwendung auf Produktionssysteme notwendig |
| KI-gestützter Digitaler Zwilling in der industriellen Automatisierung29 | Überwachtes/Unüberwachtes maschinelles Lernen (LSTM, CNN, RF, etc.) | Hohe Genauigkeit bei Fehlerprognose & -überwachung | Überwiegend klassisches maschinelles Lernen; fortschrittliche räumlich-zeitliche Deep-Learning-Modelle fehlen |
| Übersichtsarbeit zum Digitalen Zwilling in Öl- & Gasindustrie16 | Verschiedene Ansätze aus KI/ML & mechanistischer Modellierung | Umfassende Taxonomie der Modellierungsstrategien | Überwiegend konzeptionell; minimale praktische Implementierungen mit fortschrittlichem Deep Learning |
Tabelle 1: Vergleichende Analyse der Öl- und Gasförderung sowie der Steuerungsmanagementverfahren mit bestehenden Methoden. Ein Überblick über relevante Digital-Twin- und Regelstrategien unter Betonung der eingesetzten Methoden, erzielter Ergebnisse und wesentlicher Nachteile.
| Datensatzname | Bereich / System | Anzahl der Merkmale | Dauer & Größe | Dateneigenschaften |
| Simuliertes Öl & Gasproduktionsdaten | Öl & Gaserzeugungssystem (Simulation) | Multivariat (Druck, Durchfluss, Temperatur, etc.) | Für diese Forschung konzipiert | Simulierter Echtzeitsensor & Kontrollsignale, normal & Fehlerzustände |
| SWaT | Sichere Wasseraufbereitung CPS | ~51 Sensoren + Aktuatoren | ~11 Tage (~450.000+ Proben) | Normal & Anomalie (36 Angriffsszenarien) |
| WADI | Verteilung von Wasserversorgungs-CPS | ~123 Sensoren + Aktuatoren | ~16 Tage (~1 Mio. + Proben) | Normal & Anomalie (15 Anfälle) |
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 | Schichten des Policy-Netzwerks | [128, 64] |
Tabelle 3: Modellarchitektur und Hyperparameter. Wichtige architektonische Konfigurationen und Hyperparameter-Einstellungen für die neuronalen Operatoren, die Datenassimilation, das Verstärkungslernen und 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 der Vorhersageabweichung über Datensätze hinweg sowie Fehler in den Schätzungen des Anfangs- und Endzustands zeigt.
| Datensatz | Präzision | Recall | F1-Wert |
| Simuliertes Öl & 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 |
Tabelle 6: Leistung der Anomalieerkennung. Ergebnisse der Anomalieerkennung an simulierten Öl- und Gas-, SWAT- und WADI-Datensätzen, einschließlich Präzision, Sensitivität und F1-Score.
| Datensatz | Kumulierte Belohnung ↑ | Nebenbedingungsverletzungen ↓ | Kostenreduzierung (%) ↑ |
| Simuliertes Öl & Gas | 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 Reduzierung 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 | Eingabenormalisierung und Merkmalsvorbereitung |
| Gesamte Latenz der Steuerschleife | 55–82 ms | Ausführungszeit für die Steuerungsaktion von Ende zu Ende |
Tabelle 8: Bewertung der Echtzeit-Regelungsleistung. Die Latenzwerte stellen die durchschnittlichen Inferenzzeiten dar, die pro Regelungsschritt unter einer Standard-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 steuerungsbewussten Digital-Twin-Modells mit anderen Baseline-Ansätzen hinsichtlich Überwachung, Anomalieerkennung, Synchronisation und Steuerung. Bei der Vorhersagegenauigkeit werden niedrigere Werte für MAE und RMSE bevorzugt. Bei der Leistung in Anomalieerkennung und Synchronisation werden höhere Werte für den F1-Score, die AUC, die Driftreduktion und die 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 |
| Zeitfenster | 10–20 Zeitschritte |
| Neuronales 4D-Var | Assimilationsfenster | Fest (gleitendes Fenster) |
| Optimierung | Gradientenabstieg |
| 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 digitalen Zwillings-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 (Std) | 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-Controller | 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 wird verwendet, um die Echtzeit-Inferenzleistung darzustellen, GPU repräsentiert die Rechenleistung, und Skalierbarkeit beschreibt die Fähigkeit des vorgeschlagenen Modells, für industrielle Anwendungen im großen Maßstab eingesetzt zu werden.
| Methode | Überwachung des RMSE | Reduzierung der Synchronisationsdrift (%) | AUROC | F1-Score | Kumulierte Belohnung |
| 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 Digitalen Zwillingsmodells 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 | Umgang mit Beobachtungslücken |
| Normal | 100% | Niedrig | Hoch | Stabil |
| Teilweiser Verlust | 50–70% | Mäßig | Hoch | Robuste Interpolation |
| Schwerer Verlust | 20–30% | Erhöht | Mäßig | Geordneter Rückgang |
| Sensorausfall | Regionenweise 0% | Höher | Mäßig | Räumliche Generalisierung |
Tabelle 15: Robustheitsanalyse unter Datenknappheit. Die folgende Tabelle zeigt die Leistung des vorgeschlagenen neuronalen 4D-Var-Frameworks unter verschiedenen Bedingungen der Datenverfügbarkeit. Die Tabelle veranschaulicht die Robustheit des Modells, indem sie darstellt, wie sich Fehler bei teilweisem Verlust, hoher Sparsamkeit und vollständigem Verlust von Sensordaten verändern.