Ontwikkeling van een controle-bewust digital twin-framework voor real-time monitoring en optimalisatie van olie- en gasproductiesystemen

Samenvatting

Deze studie presenteert een controle-bewust digital twin-framework voor het monitoren en optimaliseren van olie- en gassystemen met behulp van machine learning en natuurkundig geïnformeerde modellen, geëvalueerd op gesimuleerde en benchmark-datasets.

Samenvatting

In dit artikel wordt een controle-bewust digital twin-framework voorgesteld voor het volgen en verbeteren van olie- en gasproductiesystemen. De methode combineert spatiotemporele graph neural networks, neurale 4D-Var data-assimilatie, reinforcement learning-gebaseerde controle en physics-informed neural operators (FNO en Neural operator model). Benchmark-datasets van cyber-fysische systemen (SWaT en WADI) en een gesimuleerde olie- en gasdataset worden gebruikt om het framework te beoordelen. De resultaten, die zijn bevestigd door meerdere onafhankelijke runs met behulp van gemiddelde maten, tonen een verbeterde prestatie in toestandschatting, anomaliedetectie en controle-optimalisatie. De beoordeling is echter beperkt tot benchmark- en gesimuleerde datasets; aanvullende validatie met werkelijke industriële gegevens is nodig om de praktische toepasbaarheid te verifiëren.

Inleiding

De behoefte aan veilige en effectieve transportnetwerken voor olie en gas is toegenomen als gevolg van de groeiende wereldwijde energievraag. Echter wordt het handhaven van de betrouwbaarheid van pijpleidingen ernstig belemmerd door verslechterende infrastructuur en toenemende systeemcomplexiteit. In dit verband heeft digital twin (DT)-technologie interesse gewekt als een datagestuurde strategie voor levenscyclusbeheer die risicobeoordeling, voorspellend onderhoud en systeemoptimalisatie in industriële omgevingen mogelijk maakt^1,2,3,4,5. De digitale transformatie in industriële systemen, van reactief onderhoud naar intelligente, datagestuurde besluitvorming, heeft de anomaliedetectie, faalvoorspelling en systeembetrouwbaarheid verbeterd^6,7,8. DT-toepassingen in procesoptimalisatie, voorspellend onderhoud en smart manufacturing zijn het onderwerp geweest van talrijke onderzoeken. Zo zijn er DT-gebaseerde modellen voorgesteld voor voorspellende faalanalyse¹⁰ en evaluatie van groene prestaties¹¹, alsook productieoptimalisatie op basis van deep learning. DT's zijn in de olie- en gasindustrie gebruikt om monitoring en anomalie-identificatie te verbeteren door middel van systeemmodellering, prestatiebeoordeling en integratie van multisensorgegevens^12,13,14,15.

Ondanks deze ontwikkelingen concentreren de meeste huidige DT-technieken zich op discrete kenmerken zoals monitoring, modellering of anomaliedetectie, met weinig integratie van real-time synchronisatie en controle. Veel systemen missen bidirectionele interactie en adaptieve leervermogens, waardoor ze functioneren als digitale schaduwen met een unidirectionele datastroom. In olie- en gassystemen, waar dynamische omstandigheden en complexe foutpatronen voortdurend leren en controle vereisen, is deze beperking bijzonder relevant. Bovendien wordt real-time closed-loop controle in complexe industriële processen niet ondersteund door huidige cyber-fysieke systeemtechnieken, die grotendeels beperkt zijn tot offline analyse¹⁶. Daarnaast vertrouwen huidige benaderingen voor het beschrijven van interacties tussen gedistribueerde sensoren hoofdzakelijk op traditionele machine learning- en filteringtechnieken, met weinig gebruik van geavanceerde modellen zoals spatio-temporele graph neural networks. Verder is er een gebrek aan onderzoek naar het gebruik van operator learning-technieken om niet-lineaire fysische dynamiek vast te leggen, vooral in olie- en gasproductiesystemen¹⁷. Hierdoor zijn bestaande DT-systemen niet in staat om real-time data-assimilatie en adaptieve controle te integreren terwijl ze complexe spatio-temporele verbindingen modelleren.

Omdat ze industriële besturingsinstellingen nabootsen, worden benchmarkdatasets zoals SWaT en WADI frequent gebruikt om DT-gebaseerde technieken in cyber-fysische systemen te valideren¹⁸. Daarnaast bieden DT's meer algemene voordelen, zoals een verhoogde productiviteit, een betere kennisuitwisseling en veilige gegevensoverdracht gedurende industriële processen^19,20. Vooral in risicovolle industrieën zoals olie en gas^21,22,23 wordt hun potentieel in veiligheidsmanagement, real-time monitoring en industriële automatisering steeds breder erkend. Desondanks zijn voor praktische implementatie nog completere en geïntegreerde oplossingen nodig die modellerings-, monitoring- en besturingsmogelijkheden combineren. Door de toepassingen en voordelen van DT's aan te tonen bij het verbeteren van de operationele efficiëntie en de wijze van besluitvorming, hebben evaluatiestudies zoals die geciteerd in^24,25 de basis gelegd. Tabel 1 toont de vergelijking van bestaande methoden met betrekking tot olie- en gasproductie en besturingsmanagement^26,27,28,29. Vanwege de hoge operationele risico's en de complexe dynamiek van olie- en gasproductiesystemen schieten traditionele onderhoudsmethoden vaak tekort voor een tijdige anomalie-identificatie en reactie^30,31,32,33. Hoewel veiligheidsmanagement en monitoring op afstand recentelijk zijn verbeterd, zijn intelligente systemen die continu in real-time kunnen waarnemen, leren en besturen nog steeds vereist.

Dit artikel stelt een controle-bewuste digital twin-architectuur voor die ruimtetijd-graafneurale netwerken, neurale 4D-Var data-assimilatie, natuurkundig geïnformeerde neurale operatoren en reinforcement learning-gebaseerde controle combineert om deze moeilijkheden te overwinnen. Het voorgestelde systeem vestigt een enkele pijplijn voor gegevensverzameling, toestandschatting, anomaliedetectie en controle-optimalisatie. Terwijl ruimtetijd-graafneurale netwerken interacties tussen verspreide sensoren registreren, worden natuurkundig geïnformeerde neurale operatoren, zoals de Fourier Neural Operator en het Deep Operator Network (Neuraal operator-model), gebruikt om complexe niet-lineaire systeemdynamiek weer te geven. Reinforcement learning maakt adaptieve gesloten-luscontrole mogelijk, terwijl neurale data-assimilatie de nauwkeurigheid van de toestandschatting verbetert. De generaliseerbaarheid van de voorgestelde methode over industriële omgevingen wordt aangetoond door deze te evalueren met benchmark-datasets van cyber-fysische systemen (SWaT en WADI) en gesimuleerde olie- en gasgegevens. Dit werk overwint belangrijke tekortkomingen van huidige DT-systemen en biedt een schaalbare en aanpasbare oplossing voor real-time monitoring en optimalisatie in olie- en gasproductiesystemen door modellering, leren en controle te combineren in een uniform raamwerk.

Protocol

Deze studie omvat geen menselijke deelnemers of dierproefpersonen. De gebruikte datasets (SWaT en WADI) zijn publiekelijk beschikbare benchmark-datasets en er is geen sprake van gevoelige of persoonlijke gegevens.

De voorgestelde aanpak vestigt een controlebewuste digitale tweeling die de integratie van geavanceerde deep learning-algoritmen, data-assimilatie en continue evaluatie en optimalisatie van olie- en gasextractiesystemen en hun aansturing mogelijk maakt. Eerst worden real-time multimodale sensorgegevens van geografisch verspreide productie-infrastructuur verzameld en gemodelleerd met behulp van een graafstructuur die de connectiviteit en interacties tussen productiecomponenten definieert. Vervolgens worden physics-informed neural operators ingezet voor de constructie van benaderingen op operator-niveau van de bepalende vergelijkingen die niet-lineaire meerfasenstroming beheersen, wat schaalbare simulatie binnen de digitale tweeling faciliteert. In dit onderzoek is de term 'Neural operator model' gebruikt om te verwijzen naar het Deep Operator Network. Om ruimtelijke en temporele afhankelijkheden tussen geografisch verspreide productieactiva te integreren, worden spatio-temporele graaf-neurale netwerken getraind binnen de physics-informed representaties. Een differentieerbare data-assimilatiebenadering die gebruikmaakt van neural 4DVAR wordt toegepast om de simulatievoorspellingen van de digitale tweeling consistent en gesynchroniseerd te houden met real-time monitoringsgegevens, zelfs bij ruisachtige en ruimtelijk onregelmatige metingen. Bovenop de real-time gesynchroniseerde digitale tweeling wordt een optimalisatietool op basis van reinforcement learning getraind voor real-time aansturing. Tot slot blijft de digitale tweeling uitgerust voor de continue bijwerking van modellen en aangepaste schaalbaarheid over een langere periode op basis van continue updates uit nieuwe beschikbare metingen.

De voorgestelde digitale tweelinarchitectuur verbetert de conventionele digitale tweelinarchitectuur door hier intelligentie, adaptiviteit en controlebewustzijn aan toe te voegen. De volledige pijplijn bestaat uit vier modules: (i) data-acquisitie/graafconstructie, (ii) toestandschatting middels neurale operatoren en ST-GNN, (iii) anomaliedetectie middels residu-gebaseerde analyse, en (iv) controle op basis van reinforcement learning.

Figuur 1: Algemene control-aware digital-twin architectuur. Anomaliedetectie, spatiotemporele graaf-neurale netwerken, neurale 4D-Var synchronisatie, physics-informed neurale operatoren en reinforcement learning-gebaseerde closed-loop regeling voor olie- en gasproductiesystemen zijn allemaal geïntegreerd in het voorgestelde control-aware digital twin-model. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 1 stelt het functionele proces voor van een digitale tweeling die is bewust van de besturing van olie- en gasproductie als een cyberfysisch systeem in een gesloten lus. Het functionele proces begint op het niveau van het olie- en gasproductiesysteem, waarbij verschillende sensoren die zijn verspreid over putten, pijpleidingen, kleppen en verwerkingseenheden gegevens vastleggen over druk, debiet, temperatuur en besturingen op een specifiek tijdstip. De diverse sensoren leggen gegevens vast in verschillende formaten, en deze gegevens worden naar de digitale tweeling gestuurd, wat de basis vormt voor het organisatorische bewustzijn op dat moment.

De voorgestelde methodologie omvat de volgende stappen:

Gegevensverzameling

Real-time informatie van sensoren met betrekking tot druk, debiet, temperatuur en regelingsinputs wordt verzameld uit olie- en gasproductieomgevingen. Andere beschikbare benchmark-datasets, zoals SWAT en WADI, worden gebruikt voor de robuustheid en generalisatiecapaciteiten. Om reproduceerbaarheid te waarborgen, wordt elke stap van de voorgestelde methodologie gespecificeerd met technische details.

[Raadpleeg Supplementary File (S1) voor de generatie van de dataset en de simulatie.]

Om de capaciteit van de conceptuele digitale tweeling aan te tonen met betrekking tot complexe en hoogdimensionale tijdreeksgegevens alsook ongebruikelijke gebeurtenissen, worden twee openbare benchmark-datasets voor industriële controlesystemen uit het onderzoeksdomein voor industriële controlesystemen gebruikt: de Secure Water Treatment (SWaT) en de Water Distribution (WADI) benchmark-datasets³⁴. Beide datasets zijn afgeleid van realistische cyber-fysieke testopstellingen en bevatten multivariate sensor- en actuatorbedieningsgegevens voor zowel normale scenario's als anomalie-/aanvalscenario's die op dergelijke systemen zijn uitgelokt. Door gebruik te maken van de SWaT- en WADI-datasets is het mogelijk om cross-domeintests uit te voeren en de generaliseerbaarheid aan te tonen voorbij cyber-fysieke testopstellingen naar bredere gedistribueerde industriële systemen met locatie- en tijdsafhankelijke relaties. De beschrijving van de datasets wordt weergegeven in Tabel 2.

Het voorgestelde onderzoek maakt gebruik van drie verschillende datasets voor de uitgebreide evaluatie van de ontwikkelde digital twin-methodologie, weergegeven in Tabel 2. De basisdataset bestaat uit gesimuleerde gegevens met betrekking tot de productieorganisatie in de olie- en gasindustrie en is bedoeld om werkelijke operationele metingen van het systeem, alsook controle- en/of foutcondities te reproduceren. Daarnaast wordt de Secure Water Treatment (SWaT) dataset geïntroduceerd als voorbeeld van een fysiek systeem op gereduceerde schaal, bestaande uit metingen van 51 sensoren en actuatoren gedurende ongeveer 11 dagen, inclusief zowel normale als kunstmatig gecreëerde anomalieën en aanvallen. Deze dataset wordt algemeen beschouwd als de standaard benchmark voor onderzoek. Een andere set is de Water Distribution (WADI) dataset, waarvan de metingen het gedrag vastleggen van het waterdistributiesysteem op gereduceerde schaal, bestaande uit ongeveer 123 variabelen over een periode van 16 dagen, inclusief verschillende aanvallen en normale bedrijfsomstandigheden. De dataset bevat in totaal 15 verschillende aanvallen. De twee bovenstaande datasets draaien hoofdzakelijk om hoogdimensionale tijdreeksen en ruimtelijke relaties, zoals gebruikelijk in industriële controlesystemen.

Preprocessing van gegevens

Voorafgaand aan de modeltraining ondergaan alle datasets een grondige voorbehandelingsprocedure om de datakwaliteit, consistentie en reproduceerbaarheid te garanderen. De preprocessingsprocedures omvatten onder meer het beheer van ontbrekende waarden, ruisfiltering, normalisatie en feature-schaling. De gedetailleerde preprocessingspipeline is beschreven in Aanvullend Bestand (S2).

Constructie van de systeemgrafiek

In de voorgestelde methodologie wordt het fysieke productiesysteem voor olie en gas weergegeven in de vorm van een graaf, die dient als wiskundige structuur voor de beschrijving van ruimtelijke correlaties tussen activa, samen met de dynamiek van de toestand van het systeem over tijd (raadpleeg het Aanvullend Bestand (S3) voor de formules voor de constructie van de systeemgraaf).

Natuurkunde leren met behulp van neurale operatoren

Voor een correcte representatie van de niet-lineaire en multidimensionale kenmerken van meerfasenstroomsystemen maakt de voorgestelde aanpak gebruik van fysica-geïnformeerde operatoren in de vorm van Fourier Neural Operators (FNO) en een neuraal operatormodel. In tegenstelling tot de traditionele methode waarbij neurale netwerken als puntgewijze mappingen worden weergegeven, maakt het concept van operatoren generalisatie mogelijk met betrekking tot verschillende randvoorwaarden. Deze eigenschap maakt deze methode geschikter voor meerfasenstroomsystemen die worden gebruikt bij de productie van olie en gas, aangezien de fysische fenomenen continu optreden. Raadpleeg aanvullend bestand (S4) voor de wiskundige vergelijkingen van de neurale operatoren.

Om de validiteit van de physics-informed neurale operatoren verder te evalueren, werden aanvullende experimenten uitgevoerd om hun vermogen te beoordelen om niet-lineaire operatorafbeeldingen te benaderen onder verschillende stromingstoestanden en controle-inputs. De resultaten tonen aan dat de operatoren de afbeelding van de niet-lineaire operator nauwkeurig kunnen benaderen, waardoor het vermogen van de operatoren wordt gevalideerd om de dynamiek van het systeem nauwkeurig te modelleren, vooral bij een hoge dimensionaliteit, tegen lagere computationele kosten dan andere solvers. De kwantitatieve evaluatie van de nauwkeurigheid van het operatorleren wordt uitgevoerd met behulp van de gemiddelde kwadratische fout tussen voorspelde en gesimuleerde systeemtoestanden, wat een consistente benaderingsprestatie over verschillende stromingscondities aantoont.

Figuur 2: Workflow van het physics-informed neural operator-raamwerk. Een voorbeeld van de physics-informed neural operator-workflow demonstreert hoe regelinputs, sensorgegevens en randvoorwaarden worden gebruikt om niet-lineaire meerfasenstroomdynamica te leren en snelle systeempredicties op operatorniveau te genereren. Klik hier om een grotere versie van deze figuur te bekijken.

In het PINO-werkproces dat wordt gepresenteerd in Figuur 2, vindt het leerproces van fysiek complexe niet-lineaire meerfasestromingen in olie- en gasprocessen op een continue manier plaats. Dit begint bij de invoerlaag, waar het model heterogene gegevens ontvangt over het fysiek evoluerende systeem. Deze invoer bestaat uit randvoorwaarden, variërend van druk, temperatuur en snelheid, evenals realtimegegevens van sensoren die in het olie- en gassysteem zijn geplaatst, en regelvariabelen zoals de opening van kleppen, de snelheid van pompen of de opening van chokes. Al deze inputs definiëren de werking van het fysiek evoluerende systeem.

Het sleutelcomponent binnen de workflow is de PINO, gerealiseerd via FNO en een Neural operator-model. De voorgestelde methode is verschillend en onderscheidend in die zin dat er neurale operatoren worden gebruikt om operatoren tussen functieruimtes te leren, in tegenstelling tot traditionele neurale netwerken, die zich doorgaans richten op puntgewijze mapping tussen inputs en outputs. Op dit niveau is FNO in staat om globale ruimtelijke relaties vast te leggen met behulp van het Fourier-domein, en is het Neural operator-model in staat om operatorrepresentaties te leren met behulp van zowel branch- als trunk-netwerken. Het geleerde model wordt geleid door fysische beperkingen om te garanderen dat de voorspelde dynamiek zich gedraagt zoals verwacht in de meerfasenstroming. De uitkomst van de neurale operator neemt de vorm aan van de continue veldoplossing van belangrijke fysische variabelen zoals druk, debieten en fasedistributie voor het gehele productie-netwerk. Bovendien zijn de outputwaarden accuraat voor het gehele veld, inclusief de putten, pijpleidingen en andere aangesloten eenheden, in plaats van dat ze op specifieke punten zoals sensoren worden gemeten. Op deze manier biedt de output van de neurale operator de mogelijkheid om snel zeer nauwkeurige simulaties te genereren die de oplossing van partiële differentiaalvergelijkingen kunnen benaderen zonder deze noodzakelijkerwijs op te lossen. Deze oplossing fungeert als de snelle oplossing in de digitale tweeling.

Het leerproces van de neurale operator maakt gebruik van kenmerken op functieniveau en systeemniveau, die zijn geëxtraheerd uit olie- en gasproductieprocessen. De kenmerken bestaan grotendeels uit de fysische toestandsvariabelen voor olie- en gasproductie, welke continue velden zijn: temperatuur, debiet, druk en fasenfracties voor olie, gas en water. De kenmerken omvatten daarnaast randvoorwaarden, waaronder inlaatdrukken, debieten en klepposities. Controle-inputs, zoals posities van verstikkleppen en pompsnelheden, worden eveneens meegenomen om responsen onder verschillende operationele scenario's te voorspellen. Het neurale operator-model ontvangt de fysische kenmerken en controleparameters, die worden verwerkt door de trunk- en branch-netwerken. De kenmerken bieden voldoende informatie voor neurale operators om nauwkeurige operatorrelaties te voorspellen die afhankelijk zijn van zowel fysische processen als operationele factoren voor snelle digital twin-simulaties.

Fourier neurale operator (FNO)

In dit onderzoek wordt een Fourier Neural Operator (FNO)³⁶ gebruikt om effectief de operator-mappings te schatten die de niet-lineaire meerfasige stromingsdynamica van olie- en gasproductiesystemen beschrijven. In tegenstelling tot traditionele neurale netwerken, die input-output-mappings benaderen voor een specifieke set punten, benadert de voorgestelde FNO mappings tussen functieruimtes, waardoor generalisatie over diverse schalen en configuraties van een systeem mogelijk wordt. In dit geval werkt de FNO door de representaties van de systeemtoestand van de oorspronkelijke ruimte naar het Fourier-domein te verplaatsen, waar spectrale convoluties effectief kunnen worden toegepast om de globale ruimtelijke afhankelijkheden binnen de systemen vast te leggen. Hierbij wordt de representatie van de inputtoestand geprojecteerd naar een hoger-dimensionale latente ruimte, waarbij transformaties en activaties iteratief worden toegepast met behulp van Fourier-transformaties, spectrale kernen en activaties. Raadpleeg het aanvullende bestand (S4.1) voor de wiskundige evaluaties van de Fourier neural operator.

Neuraal operator-model

Het neurale operator-model³⁷ dient als een alternatieve benadering binnen het neurale operator-proces voor het benaderen van een niet-lineaire fysieke operator die de dynamiek van olie- en gasproductie definieert. De gedetailleerde beschrijving van het neurale operator-model staat beschreven in Aanvullend Bestand (S4.2) (Raadpleeg Aanvullend Bestand S10 voor Algoritme 1: Physics Learning with Neural Operators for Digital Twin Algoritme).

De voorgestelde aanpak met gebruik van neurale operatoren begint daarom met het verzamelen van gepaarde voorbeelden van inputs en outputs, waarbij de inputs voorbeelden zijn van systeemcondities in de vorm van systeemgrenswaarden, regelvariabelen of fysieke kenmerken. De neurale operator wordt getraind in een iteratieve fase waarbij responsen op een inputfunctie worden voorspeld, met daaropvolgende verbetering van de afwijkingen tussen voorspellingen en werkelijke responsen met behulp van een gemiddelde kwadratische foutfunctie. Deze aanpak leidt tot optimalisatie op basis van gradiënten, waarbij een neurale operator systemische natuurkundige operatoren ontdekt die het systeemgedrag definiëren. Het getrainde systeem kan zo efficiënt systeemresponsen op verschillende inputcondities voorspellen.

Spatio-temporele systeemmodellering

In de voorgestelde digitale-tweelingarchitectuur³⁸ is het olie- en gasproductieproces weergegeven door middel van een spatio-temporele graaf. Dit is gedaan om zowel de connectiviteit als het tijdsevolutieproces van de gedistribueerde productiemiddelen te conceptualiseren.

Het productieproces van olie en gas heeft inherent een gedistribueerde structuur wat betreft putten, pijpleidinginfrastructuur en sensoren. Om ruimtelijke/temporele relaties voor een dergelijk proces efficiënt vast te leggen, maakt de digital twin-benadering die in dit artikel is ontwikkeld gebruik van spatio-temporele grafische neurale netwerken (ST-GNN). In deze benadering worden systeementiteiten weergegeven als graafknopen en worden ruimtelijke/functionele systeemrelaties aangeduid door graafranden. Temporele relaties worden uitgedrukt in termen van window-functies om systeemtoestandsovergangen en anomalieën efficiënt te voorspellen. Raadpleeg Aanvullend Bestand (S5) voor een gedetailleerde beschrijving van het Spatio-Temporal Graph Neural Network (ST-GNN).

Digital twin-synchronisatie met behulp van Neural 4D VAR

Synchronisatie van digitale tweelingen kan worden gedefinieerd als de continue synchronisatie van de virtuele tweeling met een voortdurend veranderende fysieke olie- en gasproductieomgeving op basis van real-time sensormetingen. In de praktijk doen digitale tweelingen voorspellingen over toestanden in een systeem, zoals druk, debiet en temperatuur, gebaseerd op een begrip van de dynamiek van het fysieke systeem. Desondanks ontstaat er door onzekerheden en sensorruis een onvermijdelijk verschil tussen voorspellingen en metingen. Om deze reden wordt in dit werk, om voorspellingen met werkelijke metingen te synchroniseren, een andere methode geïntroduceerd, bekend als 'data-assimilatie', om voorspellingen consistent te maken met fysieke metingen. Om de synchronisatie van het fysieke proces en de digitale tweeling te waarborgen tijdens ruizige en schaarse sensorobservaties, is een neurale 4D-Var data-assimilatiemethode geïntegreerd. De toegevoegde component is erop gericht de kloof tussen sensormetingen en modelvoorspellingen binnen een assimilatievenster te verkleinen. Door gebruik te maken van differentieerbaarheid biedt de neurale 4D-Var-methode end-to-end learning en maakt deze het mogelijk om de toestand van de digitale tweeling online te corrigeren (Raadpleeg Aanvullend Bestand (S6) voor Neurale 4D-Var Data-assimilatie, Aanvullend Bestand (S10) voor Algoritme 2: Synchronisatie van Digitale Tweelingen met behulp van Neurale 4D-Var en raadpleeg Aanvullend Bestand (S11) voor hyperparameters en trainingsdetails).

Detectie en diagnose van anomalieën

In de gesimuleerde olie- en gasdataset werden anomalieën geïnjecteerd door gecontroleerde perturbaties te introduceren, zoals drukvalstappen die 15% overschrijden, abnormale stroompieken en klefstoringen. Datapunten die een ingestelde drempelwaarde overschreden of afweken van een statistische basisverdeling worden geïdentificeerd als anomalieën. Voor de wiskundige formulering van het Anomaly Detection Framework (Aanvullend bestand, S7).

Closed-loop optimalisatie en controle

In het voorgestelde systeem kan gesloten-lusregeling worden bereikt door de harmonieuze DT en de reinforcement learning-module voor optimalisatie te integreren. Wat betreft de digitale tweelingen, voorspelt de digitale tweeling het toekomstige gedrag van het olie- en gasproductieproces onder verschillende omstandigheden, zoals de positionering van verschillende kleppen, pompsnelheden of debietregeling. Dit maakt de ontwikkeling van een veiligheidsgebied mogelijk waarin de regelagent het resultaat van gewenste acties kan evalueren zonder dat deze in het proces in de echte wereld hoeven te worden gerealiseerd. Door zo gebruik te maken van de voorspellingen van de digitale tweeling, kan de reinforcement learning-module operationele doelen optimaliseren, waaronder de efficiëntie van het proces, het energieverbruik en de processtabiliteit, zonder beperkt te worden door procesbeperkingen. Raadpleeg aanvullend bestand (S8) voor de wiskundige beschrijving van de gesloten-lusoptimalisatie en -regeling.

Continu leren en aanpassen

Industriële systemen worden geconfronteerd met niet-stationaire omgevingen met verouderde apparatuur, variërende bedrijfsomstandigheden en onvoorziene verstoringen. Om hiervoor een aanpak te ontwerpen, bevat het voorgestelde digitale tweelingmodel een continu lerend onderdeel dat neurale operatoren, grafmodellen en beleid upgrade op basis van incrementeel leren. Wanneer er grote variaties en niet geïdentificeerde patronen worden waargenomen, wordt het leerproces gefaciliteerd door nieuwe gegevens in te voegen en parameters geassocieerd met het digitale tweelingmodel aan te passen.

Continu leren verbetert ook de robuustheid van zowel de besturingsstrategieën als de detectie van anomalieën, aangezien het het systeem helpt om nieuwe of onbekende gedragspatronen te identificeren. Feedback van anomalieën in de besturingsresultaten wordt toegepast om modellen aan te passen en een reinforcement learning-strategie te optimaliseren, wat effectieve besluitvormingsprocessen waarborgt terwijl het systeem zich blijft veranderen en ontwikkelen. Dit adaptieve leerproces transformeert een digitale tweeling van een passieve waarnemer naar een intelligent systeem dat in staat is om real-time optimalisatieprocessen van complexe olie- en gasproductiesystemen over een langere periode in stand te houden, gebruikmakend van een technologische infrastructuromgeving die rigoureus en effectief is. In Tabel 3 worden de modelarchitectuur en hyperparameters weergegeven.

Raadpleeg Aanvullend Bestand (S9) voor de wiskundige formulering van Continuous Learning Adaptation.

De module voor gegevensvoorbewerking is verantwoordelijk voor het transformeren van ruwe sensorgegevens naar gestructureerde en genormaliseerde inputs die geschikt zijn voor de daaropvolgende leermodules. Specifiek aggregeert deze module historische sensorwithmetingen, verwijdert ruis en standaardiseert kenmerken om een consistente datakwaliteit te waarborgen. De bewerkte gegevens worden, samen met randvoorwaarden en controle-inputs, vervolgens aangeboden als input aan de module van de physics-informed neural operator (FNO/Neural operator-model), die de onderliggende niet-lineaire fysische relaties leert en een continue schatting van de systeemtoestand genereert. Deze geschatte systeemtoestand wordt vervolgens door de module voor controle-optimalisatie gebruikt om controlestrategieën op een dynamische manier bij te werken. De geoptimaliseerde controle-outputs worden vervolgens doorgestuurd naar de besluitvormingsmodule, die de resultaten interpreteert en uitvoerbare beslissingen genereert voor de systeemoperatie. Er is een feedbackmechanisme geïntegreerd om het systeem continu bij te werken op basis van nieuw waargenomene gegevens, waardoor de aanpasbaarheid en prestaties in de loop van de tijd worden verbeterd. Om de duidelijkheid en reproduceerbaarheid te vergroten, wordt een architectuur op systeemniveau gepresenteerd om de interactie tussen de modules te illustreren. Het voorgestelde systeem bestaat uit vier hoofdcomponenten, namelijk gegevensvoorbewerking, physics-informed neural operator-modellering, controle-optimalisatie en besluitvorming. Deze modules zijn onderling verbonden via goed gedefinieerde input-output-interfaces, wat een naadloze datastroom, iteratief leren en gecoördineerde werking over het gehele digital twin-systeem garandeert.

Figuur 3: Stroomdiagram van de architectuur voor de optimalisatie van de regeling op systeembreid niveau. Het toont de algemene systeemsarchitectuur van het voorgestelde werk. Klik hier om een grotere versie van deze figuur te bekijken.

De systeemarchitectuur geeft een compleet beeld van het model dat kan worden gebruikt voor het modelleren en optimaliseren van complexe niet-lineaire fysische dynamiek in olie- en gasproductiesystemen. Dit begint met de invoergegevens, die de parameters bevatten die voor de simulatie worden gebruikt en de randvoorwaarden die worden gehanteerd om de fysische omgeving weer te geven. Deze gegevens passeren eerst de module voor gegevensvoorbehandeling, waar normalisatie en feature-extractie plaatsvinden. Dit bereidt de gegevens voor op de volgende module, waarin de neurale operator, FNO, en het neurale operator-model worden gebruikt. Deze modellen worden ingezet om de systeemtoestanden te voorspellen, en deze informatie wordt gebruikt om de systeemparameters te optimaliseren met behulp van de optimalisatiemodule. De geoptimaliseerde resultaten worden vervolgens gebruikt om beslissingen te nemen via de besluitvormingsmodule. Deze module helpt bij het bepalen van de efficiëntie van het systeem. Het raamwerk omvat ook de uitvoerresultaten, waarbij de resultaten van het optimalisatieproces worden getoond in Figuur 3. Daarnaast is er een feedbacklus opgenomen, waarbij de gegevens na de uitvoerresultaten opnieuw door de module voor voorbehandeling passeren.

Resultaten

Om de geldigheid van de natuurkundig geïnformeerde neurale operatoren te testen, werden aanvullende prestatietests uitgevoerd om de prestaties van FNO en het neurale operator-model te vergelijken bij het modelleren van de dynamiek van het niet-lineaire meerfasenstromingssysteem. De resultaten toonden aan dat FNO beter presteert in termen van globale ruimtelijke consistentie, terwijl het neurale operator-model beter presteert in de aanpasbaarheid aan heterogene inputs.

Om een effectieve training en evaluatie van het voorgestelde digitale tweeling-framework te garanderen, werden alle experimenten uitgevoerd in een high-performance computing-omgeving. Python werd gebruikt als de primaire programmeertaal voor de implementatie. Een populair deep learning-framework met GPU-acceleratie via CUDA en een geoptimaliseerde deep neural network-bibliotheek werden gebruikt voor het creëren van deep learning-modellen, zoals physics-informed neural operators, spatiotemporele grafische neurale netwerken en reinforcement learning-modules. Om grootschalige verwerking van multivariate tijdreeksgegevens en modeltraining te faciliteren, werden de tests uitgevoerd op een workstation met een multi-core CPU, een high-performance graphics processing unit (GPU) met dedicated geheugen en minimaal 64 GB systeemgeheugen. Om compatibiliteit met het deep learning-ecosysteem te waarborgen, was de besturingsomgeving gebaseerd op een Linux-distributie.

Mini-batch gradient-gebaseerde optimalisatie werd gebruikt om de modellen te trainen, en er werden talrijke onafhankelijke runs van de experimenten uitgevoerd om de statistische betrouwbaarheid te garanderen. De voorgestelde methodologie kan in vergelijkbare high-performance computing-omgevingen worden gereproduceerd dankzij de verstrekte computationele opzet.

Om de resultaten te valideren via herhaalde proeven, werden alle experimenten herhaald in meerdere onafhankelijke runs met verschillende willekeurige initialisatiewaarden. In het geval van experimenten met gesimuleerde gegevens is het model uitgevoerd voor (N) onafhankelijke proeven, waarbij de waarden worden gerapporteerd als gemiddelde prestatie en standaarddeviatie. In het geval van de experimenten met de benchmarkgegevens, SWAT en WADI, zijn de evaluatiemetrieken MAE, RMSE, F1 Score en AUROC berekend op de testgegevens en gemiddeld over de herhaalde experimenten om de resultaten statistisch robuust te maken. Dit is gedaan om een betrouwbare schatting van de prestaties van het model te verkrijgen en om ervoor te zorgen dat de prestatieverbeteringen niet te wijten zijn aan een enkel experiment.

Desalniettemin hebben de resultaten bevestigd dat het integreren van physics-informed learning, ruimtelijk-temporele modellering en data-assimilatietechnieken in een geïntegreerd digital twin-systeem consistent leidt tot verbeterde monitoringprestaties onder de geëvalueerde omstandigheden, een snellere convergentie en optimale realtime beslissingen in vergelijking met conventionele methoden. Tabel 4 toont de simulatieomgeving van het voorgestelde werk.

De simulatieomgeving is gemodelleerd om een nauwkeurige weergave te vormen van real-time operaties binnen de industrie. Gegevens van sensoren worden met korte tussenpozen verzameld om een streaming-dataomgeving weer te geven, waarbij gecontroleerde hoeveelheden ruis worden toegevoegd om de robuustheid te testen. Verschillende fout- en storingsscenario's worden gesimuleerd om de mogelijkheden voor anomaliedetectie en aansturing te testen. Voor de evaluatie wordt een sliding window-benadering gebruikt om real-time operaties op een continue wijze weer te geven, geschikt voor dynamische updates door de voorgestelde digital twin voor voorspellingen en toestandsynchronisaties en voor het genereren van regelacties.

Bij het ontwerp van het besturingsbewuste digitale tweemaalmodel moeten verschillende evaluatiecriteria in overweging worden genomen voor een gelijktijdige beoordeling. Ten eerste worden voor de bewakings- en toestandsvoorspellingscapaciteit van de digitale tweelingsarchitectuur regressiecriteria zoals de Mean Absolute Error (MAE), Root Mean Squared Error (RMSE) en Mean Absolute Percentage Error (MAPE) gebruikt. MAE en MAPE berekenen de nauwkeurigheid van neurale operatoren en spatio-temporele grafiekmodellen op belangrijke systeemparameters zoals druk, debiet en temperatuur, en zijn geldig voor de continue multivariabele gegevens van sensoren die worden gebruikt in het olie- en gasproductieproces.

Ten tweede wordt het synchronisatieproces van de digitale tweeling zelf beoordeeld, met betrekking tot de fout in de toestandschatting en het verwijderen van de predictiedrift. De fout in de toestandschatting meet de mate van afwijking tussen de werkelijke toestand van het schema en de gesynchroniseerde toestand van de digitale tweeling; dit is een indicator voor hoe effectief de ontwikkelde neurale 4D-Var-module de werkelijke systeemtoestandswaarnemingen gebruikt om de voorspellingen van de systeemdynamica te verbeteren. Ten slotte geeft de reductie van de predictiedrift de prestatiewinst van de voorgestelde aanpak ten opzichte van het niet-gesynchroniseerde digitale tweeling-systeem aan, dankzij het vermogen van de methode om bestand te zijn tegen ruisachtige of schaarse systeemwaarnemingen.

Ten derde worden bij het detecteren van anomalieën en foutdiagnose classificatiegerichte evaluatieparameters zoals precisie, recall, F1 en de oppervlakte onder de ROC-curve (AUC) gebruikt. Terwijl recall verwijst naar het vermogen om werkelijke fouten nauwkeurig te diagnosticeren zonder er een over het hoofd te zien, kan precisie worden beschouwd als een maat voor de validiteit van een herkende abnormaliteit. De prestaties kunnen op passende wijze worden gemeten met behulp van een F1-metriek. De evaluatieparameters worden over het algemeen toegepast in een SWAT- en WADI-systeem.

Tabel 5 wordt gebruikt om de effectiviteit van de voorgestelde neurale 4D-Var synchronisatiemodule te beoordelen voor de reductie van voorspellingsdrift tussen de fysieke en cybermodellen voor drie datasets. De initiële toestandsfout komt overeen met de fout tussen modellen zonder synchronisatie, en de uiteindelijke toestandsfout komt overeen met de modelfout na toepassing van synchronisatie. Er is een significante verbetering in de toestandsfout te zien voor alle drie de datasets, waarbij de grootste verbetering (72,2%) optreedt voor het gesimuleerde olie- en gassysteemmodel vanwege de relatieve eenvoud in vergelijking met SWAT (63,6%) en WADI (57,7%). Deze resultaten zijn consistent met hun respectievelijke complexiteitsniveaus en "ruis"-niveaus. De verhoogde mate van driftreductie voor de gesimuleerde olie- en gasdataset kan worden verklaard door het feit dat de gesimuleerde omgeving over het algemeen een lager ruisniveau, minder willekeurige perturbaties en een beter gedefinieerd systeemgedrag vertoont in vergelijking met andere realistische cyber-fysische systemen zoals SWAT en WADI. Deze argumentatie moet echter als een hypothese worden beschouwd en is op basis van de resultaten van dit werk niet kwantitatief bewezen. In toekomstig werk kunnen complexiteitsmaten zoals entropie, ruisvariantie en systeemdimensionaliteit worden gebruikt om dit argument te valideren.

Figuur 4: Monitoring van de prestaties over verschillende datasets (Simulated Oil and Gas, SWAT en WADI), geëvalueerd aan de hand van de foutmetrieken MAE, RMSE en MAPE (%). De staven vertegenwoordigen de gemiddelde waarden verkregen uit meerdere onafhankelijke runs (n = 5). De foutenbalken geven standaarddeviatie (SD), terwijl dunnere, overliggende foutenbalken de standaardfout van het gemiddelde (SEM = SD/√n). De SD- en SEM-waarden voor elke metriek worden expliciet weergegeven in de legenda. Statistische significantie wordt beoordeeld met behulp van onafhankelijke t-toetsen ten opzichte van de basisdataset (Simulated Oil and Gas), en overeenkomstige p-waarden worden boven de balken weergegeven. Klik hier om een grotere versie van deze figuur te bekijken.

De prestatieanalyse van het monitoringproces dat is uitgevoerd door het voorgestelde digitale tweelingmodel op datasets zoals Simulated Oil and Gas, SWAT en WADI, rekening houdend met de evaluatiecriteria MAE, RMSE en MAPE, is weergegeven in Figuur 4. De gesimuleerde Oil and Gas dataset vertoont de laagste foutmeting voor alle criteria, wat duidelijk de nauwkeurigheid van de systeemtoestandsmonitoring in een gesimuleerde omgeving aantoont waarin de systeemparameters nauwkeurig zijn gemodelleerd. De SWAT- en WADI-datasets vertonen fouten in de systeemmonitoring die geleidelijk toenemen vanwege de complexiteit van echte systemen, samen met hun inherente ruis en nonlineariteiten, die niet voorkomen in gesimuleerde systemen. De geringe stijging van de foutmetingen van MAE en RMSE bevestigt echter de stabiliteit van de systeemprestaties van het digitale tweelingmodel onder dergelijke omstandigheden.

Figuur 5: Prestaties van anomaliedetectie over datasets (Simulated Oil and Gas, SWAT en WADI), geëvalueerd met behulp van de F1-score en AUC. De staven vertegenwoordigen de gemiddelde waarden verkregen uit meerdere onafhankelijke runs (n = 5). Foutenbalken geven de SD aan, terwijl dunnere overliggende foutenbalken de SEM = SD/√n vertegenwoordigen. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5 toont de evaluatieanalyse van het voorgestelde systeem voor anomaliedetectie op basis van de F1-score en AUC voor de drie datasets. De Oil and Gas-dataset, de gesimuleerde dataset, presteert het beste met de hoogste F1-score en AUC. Hoewel er een geleidelijke afname in prestaties is, wordt er een hoge AUC-prestatie behaald over de datasets heen, met lichte variaties afhankelijk van de complexiteit van de dataset. De uitlijning van zowel de F1-score- als de AUC-curves voor de drie datasets toont een optimaal omslagpunt, wat wijst op een goede trade-off en het voorgestelde systeem een betrouwbaar systeem maakt voor vroege foutdetectie en intrusie. Het voorgestelde systeem presteert uitmuntend bij het herkennen van anomalieën voor de drie datasets.

Precisie, recall en F1-scores van het voorgestelde digital twin-model bij het detecteren van anomalieën in verschillende datasets worden getoond in Tabel 6 hieronder. Hoewel hoge precisiewaarden aantonen dat de meeste geïdentificeerde interessante punten inderdaad ware interessante punten zijn, tonen hoge recall-waarden aan dat de meeste werkelijke interessante punten correct worden geïdentificeerd zonder belangrijke punten te missen. De demonstratieomgeving in het olie- en gasdomein registreert de hoogste F1-score van 0,94, terwijl iets lagere resultaten voor SWAT en WADI kunnen worden toegeschreven aan hogere niveaus van ruis en interdependencies tussen sensoren en aanvalspatronen. Desondain valideren de hoge F1-scores dat een robuuste spatio-temporele graafmodellering en synchronisatie voor het detecteren van anomalieën is bereikt.

Figuur 6: Receiver Operating Characteristic (ROC)-curves voor anomaliedetectie in de Simulated Oil and Gas-, SWAT- en WADI-datasets. De prestaties worden geëvalueerd aan de hand van de Area Under the Curve (AUC). De standaarddeviatie (SD) en de standaardfout van het gemiddelde (SEM = SD/√n, n = 5) van de AUC-waarden zijn opgenomen in de legenda voor elke dataset. De diagonale stippellijn vertegenwoordigt de prestaties bij willekeurige classificatie. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6 toont de Receiver Operating Characteristic (ROC)-curves, samen met de bijbehorende AUROC-maten voor de drie verschillende datasets, waarbij de nauwkeurigheid van de anomaliedetectie onafhankelijk van de keuze van de drempelwaarde wordt beoordeeld. In de ROC voor de gesimuleerde Oil and Gas-dataset komt de regio opvallend dicht bij de oorsprong, waardoor een AUROC van 0,99 wordt vastgesteld, wat de vrijwel foutloze classificatiecapaciteit bevestigt. In het geval van de SWAT- en WADI-datasets liggen de AUROC-waarden respectievelijk op 0,97 en 0,95, waarmee een goede classificatiecapaciteit wordt aangetoond ondanks de aanwezigheid van hogere onzekerheden. Omdat de ROC-curves ver verwijderd zijn van de lijnen die een willekeurige classificatie vertegenwoordigen, is gegarandeerd dat het ontwikkelde digital twin-model aanzienlijk beter presteert dan het toevalsniveau voor nauwkeurige anomaliedetectie. De Receiver Operating Characteristic (ROC)-curves in Figuur 6 tonen de prestaties van het voorgestelde digital twin-model bij anomaliedetectie op de drie datasets. De AUROC-scores zijn 0,99 voor de gesimuleerde Oil and Gas-dataset, 0,97 voor de SWAT-dataset en 0,95 voor de WADI-dataset, wat wijst op het onderscheidend vermogen van het model in zowel gesimuleerde als benchmark industriële controlesystemen.

Tabel 7 geeft een beknopte beschrijving van de cumulatieve beloningen en schendingen van restricties voor de reinforcement learning-module voor closed-loop controle, geïntegreerd met de digitale tweeling. De behaalde cumulatieve beloning vertegenwoordigt de algehele controleprestatie op basis van operationele doelen, terwijl schendingen van restricties het totaal aantal schendingen tijdens de werking vertegenwoordigen. De cumulatieve beloningen met nul schendingen van restricties zijn het hoogst (96,5) voor onze olie- en gasmodelomgeving, wat duidt op een optimale controle en operationele prestatie.

Hoewel de SWaT- en WADI-datasets iets lagere cumulatieve beloningen en enkele schendingen van beperkingen vertonen, bewijst de algemene verbetering in operationele kostenbesparingen dat het voorgestelde digital twin-framework effectief blijft voor veilige en optimale besturing in dynamische omgevingen.

De real-time haalbaarheid van het voorgestelde besturingssysteem op basis van reinforcement learning (RL) werd beoordeeld door de end-to-end inferentielatentie van alle gebruikte computationele modules binnen het voorgestelde systeem te analyseren. De experimentele waarnemingen in Tabel 8 laten zien dat de neurale operator-module, het FNO/neurale operator-model en de spatio-temporele graph neural network (ST-GNN)-module een matige computationele overhead hebben vanwege de temporele kenmerken die tijdens het proces worden gebruikt. Daarentegen heeft de beslissingsmodule voor de reinforcement learning-policy een zeer lichte computationele overhead. De cumulatieve end-to-end inferentietijd van alle modules die in het voorgestelde systeem zijn gebruikt, valt binnen het bereik van typische industriële besturingscycli. Specifiek is waargenomen dat de gemiddelde end-to-end inferentielatentie binnen 50–150 milliseconden per besturingscyclus ligt, wat compatibel is met standaard industriële besturingsvereisten. Dit toont aan dat het voorgestelde systeem kan worden gebruikt voor real-time en near real-time toepassingen. Houd er echter rekening mee dat de end-to-end inferentietijd kan variëren afhankelijk van de gebruikte hardware. In dat opzicht kan het voorgestelde systeem worden ingezet voor real-time toepassingen.

Er wordt een kwantitatieve vergelijking uitgevoerd tussen de FNO, het Neural operator-model en baseline-benaderingen, zoals traditionele numerieke solvers en huidige operator-leertechnieken, om de effectiviteit van physics-informed neural operators verder vast te stellen. De evaluatie uit Tabel 9 richt zich op de mate waarin elk model de dynamiek en niet-lineaire dynamiek van het systeem vastlegt. Standaardmetrieken, waaronder de gemiddelde kwadratische fout (MSE), de gemiddelde absolute fout (MAE) en de relatieve fout tussen de verwachte en de ground-truth systeattoestanden, worden gebruikt om de prestaties te evalueren. Op basis van de experimentele resultaten presteren zowel de FNO als het Neural operator-model beter dan de baseline-benaderingen bij het modelleren van complexe meerfasendynamiek door lagere foutwaarden te behalen. Specifiek legt het Neural operator-model succesvol niet-lineaire operator-mappings vast onder verschillende inputcondities, terwijl de FNO goed presteert in het leren van globale ruimtelijke afhankelijkheden.

De effectiviteit van de voorgestelde neurale operator-modellen wordt aangetoond door de kwantitatieve vergelijkende resultaten in Tabel 9. In vergelijking met basisnumerieke solvers en huidige operator-leermethoden behalen zowel FNO als het neurale operator-model merkbaar lagere foutmetrieken. In het bijzonder vertoont FNO de laagste relatieve fout en gemiddelde kwadratische fout, wat de superieure capaciteit aantoont om de dynamiek van complexe niet-lineaire systemen vast te leggen. Daarnaast vertoont het neurale operator-model significante winst, vooral bij het leren van niet-lineaire operator-mappings onder verschillende omstandigheden. Onafhankelijke t-toetsen worden gebruikt om de statistische significantie te beoordelen; p-waarden kleiner dan 0,05 duiden op een substantiële verbetering ten opzichte van de basisnumerieke solver. De p-waarden verkregen voor beide voorgestelde modellen zijn kleiner dan 0,001, wat aangeeft dat de waargenomen prestatieverbeteringen statistisch significant zijn en niet het resultaat van toeval. Bovendien vertonen de voorgestelde technieken aanzienlijk kortere inferentietijden, wat hun toepasbaarheid voor real-time digital twin-toepassingen demonstreert.

Om de duidelijkheid van de experimentele evaluatie van het voorgestelde control-aware digital twin-model verder te verbeteren, wordt in een geïntegreerd tabelformaat een uitgebreide kwantitatieve vergelijking van de voorgestelde benadering met andere methoden gegeven. Deze vergelijking van diverse prestatieparameters, zoals de Mean Absolute Error (MAE), Root Mean Squared Error (RMSE), F1-score, Area Under the Curve (AUC), vermindering van predictiedrift en cumulatieve beloning, maakt een geïntegreerde evaluatie mogelijk van de monitoringsnauwkeurigheid, de anomaliedetectiecapaciteit, de synchronisatie-efficiëntie en de controleprestaties van de methoden in Tabel 10. Daarnaast worden diverse parameters van de trainingsprocedure, zoals batchgrootte, leersnelheid, aantal epochs en de module-specifieke trainingsconfiguratie, expliciet vermeld in Tabel 11 om de transparantie van de experimentele evaluatie te vergroten. Verder worden de schaalbaarheid en computationele efficiëntie van de voorgestelde benadering, uitgedrukt in inferentielatentie, GPU-benutting en real-time haalbaarheid, geëvalueerd in Tabel 12.

De voorgestelde digitale tweeling kan uitgebreid worden vergeleken met andere methoden, zoals de op fysica gebaseerde synchronisatiemethode, LSTM- en GNN-modellen, op basis van bewakingsnauwkeurigheid, synchronisatievermogen, anomaliedetectie en regelprestaties, zoals weergegeven in Tabel 13. Fysische digitale tweelingen zijn minder adaptief en hebben een hogere foutmarge vanwege hun statische benaderingen. Neurale benaderingen, zoals LSTM en GNN, presteren beter, maar combineren fysica, synchronisatie en regeling niet effectief. De voorgestelde digitale tweeling heeft de beste waarden behaald met een RMSE van 0,20, 72,2% in driftreductie, 0,94 in AUROC, 0,92 in F1-Score en 95,5 in cumulatieve beloning. De bovenstaande resultaten geven duidelijk aan dat het effectief combineren van physics-informed neurale operatoren, graafmodellering voor ruimte en tijd, neurale data-assimilatie en reinforcement learning een nauwkeurig, robuust en regelbewust digitaal tweelingsysteem aanzienlijk kan verbeteren.

Figuur 7: Trainingskosten en resourcegebruik van verschillende componenten van het voorgestelde systeem, inclusief FNO, het neurale operatormodel, GNN, de RL-module en het totale framework. De metrieken omvatten de trainingstijd (uren) en het GPU-geheugengebruik (GB). De waarden vertegenwoordigen de gemiddelde resultaten over meerdere onafhankelijke runs (n = 5). De foutenbalken geven de SD aan, terwijl de dunnere overlappende foutenbalken de SEM = SD/√n vertegenwoordigen, wat de variabiliteit tussen de runs weerspiegelt. Klik hier om een grotere versie van deze figuur te bekijken.

In Figuur 7 worden de trainingskosten en de resourcebenutting van de voorgestelde aanpak weergegeven op basis van hun hoofdcomponenten. Er kan worden waargenomen dat de neurale operator-modellen, zoals FNO en het neurale operator-model, aanzienlijke trainingskosten en GPU-geheugenbenutting hebben in vergelijking met het totale systeem. De aanzienlijke trainingskosten zijn het gevolg van het hoog-dimensionale leervermogen van de operators. De grafen-neurale netwerk- en reinforcement-learningmodules blijken lagere trainingskosten te hebben, maar ze zijn essentieel voor het aanleren van ruimtelijke afhankelijkheden en adaptieve controle. Hoewel het voorgestelde systeem hogere trainingskosten heeft, blijkt de inferentietijd efficiënt te zijn. Vergeleken met referentiemethoden, zoals stand-alone LSTM- of GNN-modellen, brengt het voorgestelde systeem hogere trainingskosten met zich mee door de integratie van meerdere modules, maar het bereikt een verbeterde nauwkeurigheid en robuustheid, zoals aangetoond in Tabel 13.

Over het geheel genomen combineert de voorgestelde digitale tweelingstructuur verschillende elementen, waarvan elk een specifiek probleem aanpakt bij het modelleren en besturen van complexe industriële systemen. In tegenstelling tot traditionele numerieke solvers, die vaak rekenintensief zijn en minder geschikt voor realtime-toepassingen, worden physics-informed neural operators, zoals de Fourier Neural Operator en het Deep Operator Network, gebruikt om niet-lineaire en hoogdimensionale systeemdynamiek effectiever vast te leggen. Eenvoudigere machine learning-modellen, zoals feedforward- of recurrente neurale netwerken, kunnen worden ingezet, maar zij leren doorgaans geen fundamentele fysische operatoren en kunnen niet generaliseren over verschillende systeemconfiguraties. Afhankelijkheden tussen gedistribueerde sensoren, die inherent verbonden zijn in olie- en gasproductiesystemen, worden gemodelleerd met behulp van het spatio-temporal graph neural network. Conventionele methoden, zoals onafhankelijke tijdreeksmodellen of eenvoudige recurrente architecturen, kunnen leiden tot een lagere nauwkeurigheid bij toestandschatting en anomaliedetectie, aangezien zij ruimtelijke koppelingen niet expliciet weerspiegelen. Adaptieve regeling, die gebruikmaakt van de reinforcement learning-module, stelt het systeem in staat om operationele keuzes te maximaliseren in dynamische en onzekere omgevingen. Regelgebaseerde of statische regelstrategieën zijn daarentegen rigide en niet in staat zich aan te passen aan veranderende systeemcondities of onvoorziene verstoringen. Hoewel er minder complexe alternatieven bestaan, behandelen deze doorgaans slechts bepaalde delen van het probleem en zijn ze onvoldoende om een geïntegreerde oplossing te bereiken die nauwkeurige modellering, realtime monitoring en adaptieve regeling combineert. Door ervoor te zorgen dat elk onderdeel bijdraagt aan de algehele systeemprestatie, verbetert het voorgestelde modulaire ontwerp de nauwkeurigheid, schaalbaarheid en robuustheid in uitdagende industriële omgevingen.

Studie / Werk	Gebruikte technieken	Resultaten	Beperkingen
Deep Neural Operator voor Digital Twin26	DeepONet neurale operator als surrogaatmodel	Realtime voorspelling; grootteordes sneller dan physics-simulatoren	Vereist zorgvuldige sensorplaatsing; evaluatiemethoden moeten worden verbeterd
Virtual Sensing-enabled Digital Twin met DeepONet27	DeepONet virtuele sensoren in DT	Snelle voorspellingen; realtime inferentie; uitgebreide toestandschatting	Spectrale bias beïnvloedt hoogfrequente fenomenen; hybride modellering is noodzakelijk
Digital Twin met DRL + NMPC in Control28	Integratie van Deep Reinforcement Learning & Nonlinear MPC	Verbeterde regelprestaties; lagere trackingfout en adaptieve twin-updates	Toepassing beperkt tot vaatuigen; verdere toepassing op productiesystemen is nodig
AI-gestuurde Digital Twin in industriële automatisering29	Supervised/Unsupervised ML (LSTM, CNN, RF, etc.)	Hoge nauwkeurigheid bij foutvoorspelling & monitoring	Voornamelijk klassieke ML; gebrek aan geavanceerde spatio-temporele deep-modellen
Review van Digital Twin in Olie & Gas16	Diverse AI/ML & mechanistische modellering	Uitgebreide taxonomie van modelleringsstrategieën	Grotendeels conceptueel; minimale praktische implementaties met geavanceerde DL

Tabel 1: Vergelijkende analyse van olie- en gasproductie en beheersingsmanagement met bestaande methoden. Een overzicht van relevante digitale tweelingen en beheersingsstrategieën, met nadruk op de gehanteerde methoden, behaalde resultaten en significante tekortkomingen.

Naam dataset	Domein / Systeem	Aantal kenmerken	Duur & Grootte	Datakenmerken
Gesimuleerde olie- & gasproductiegegevens	Olie- & gasproductiesysteem (simulatie)	Multivariaat (druk, debiet, temperatuur, etc.)	Ontworpen voor dit onderzoek	Gesimuleerde real-time sensorsignalen & controlesignalen, normale toestanden & fouttoestanden
SWaT	Secure Water Treatment CPS	~51 sensoren + actuatoren	~11 dagen (~450k+ monsters)	Normaal & anomalie (36 aanvalscenario's)
WADI	Water Distribution CPS	~123 sensoren + actuatoren	~16 dagen (~1M+ monsters)	Normaal & anomalie (15 aanvallen)

Tabel 2: Beschrijving van de dataset.
Beschrijving van de dataset. Samenvatting van de datasets van de studie, inclusief domein, aantal kenmerken, duur en data-attributen.

Component	Hyperparameter	Waarde
ST-GNN	Aantal lagen	3
ST-GNN	Verborgen dimensie	64
ST-GNN	Lengte van het temporele venster	12
Neurale Operator (FNO)	Aantal Fourier-modi	16
Neurale Operator (FNO)	Aantal lagen	4
DeepONet	Grootte van het vertakkingsnetwerk	[128, 128]
DeepONet	Grootte van het stamnetwerk	[128, 128]
Data-assimilatie	Assimilatievenster	10 tijdstappen
Optimizer	Adam
Leersnelheid	0.001
Batchgrootte	64
RL-controller	Disconteringsfactor (γ)	0.99
RL-controller	Lagen van het beleidsnetwerk	[128, 64]

Tabel 3: Modelarchitectuur en hyperparameters. Belangrijke architecturale configuraties en hyperparameterinstellingen voor de neurale operatoren, data-assimilatie, reinforcement learning en de componenten van het spatiotemporele grafische neurale netwerk.

Parameter	Beschrijving
Simulatieplatform	Python (v3.10) met PyTorch (v2.1)
Ondersteunende bibliotheken	NumPy (v1.24), SciPy (v1.10), Pandas (v1.5)
Hardware	NVIDIA RTX serie GPU (bijv. RTX 3090, 24 GB VRAM), Intel Core i7/i9 CPU, 32–128 GB RAM
GPU-acceleratie	CUDA Toolkit (v11.8), cuDNN (v8.x)
Besturingssysteem	Linux (Ubuntu 20.04) / Windows 10
Ontwikkelomgeving	Jupyter Notebook / VS Code
Bemonsteringsinterval	1–10 seconden (configureerbaar voor real-time simulatie)
Modellering van sensorruis	Gaussische ruis (σ = 0,01–0,05)
Foutscenario's	Sensorfouten, stromingsverstoringen en controle-anomalieën
Evaluatiemodus	Real-time simulatie op basis van sliding window
Beschikbaarheid van code	Repository-link wordt op aanvraag / publicatie verstrekt voor reproduceerbaarheid

Tabel 4: Simulatieomgeving. Informatie over de hardwareconfiguratie, het softwareplatform, het bemonsteringsinterval, de ruismodellering en de foutscenario's in de experimentele opstelling.

Dataset	Initiële Toestandsfout	Eindtoestandsfout	Driftreductie (%)
Gesimuleerde Olie & Gas	0.9	0.25	72.2
SWAT	1.1	0.4	63.6
WADI	1.3	0.55	57.7

Tabel 5: Experimenteel resultaat van Digital Twin-synchronisatie. Kwantitatieve beoordeling van de synchronisatieprestaties, die de reductie van de voorspellingsdrift over datasets en de fouten in de initiële en uiteindelijke toestatsschatting weergeeft.

Dataset	Precisie	Herinnering	F1-score
Gesimuleerde olie & 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

Tabel 6: Prestaties van anomaliedetectie. Resultaten voor anomaliedetectie op gesimuleerde olie- en gas-, SWAT- en WADI-datasets, inclusief precisie, recall en F1-score.

Dataset	Cumulatieve Beloning ↑	Schendingen van Beperkingen ↓	Kostenreductie (%) ↑
Gesimuleerde Olie & Gas	96.5 ± 0.8	0	18.7
SWAT	92.1 ± 1.2	1	14.3
WADI	89.7 ± 1.5	2	11.6

Tabel 7: Experimentele resultaten van gesloten-lusregeling en optimalisatie. Cumulatieve beloning, schendingen van randvoorwaarden en vermindering van de operationele kosten worden gebruikt om de prestaties van de gesloten-lusregeling te evalueren.

Module	Gemiddelde latentie (ms)	Beschrijving
Neurale operator (FNO/DeepONet)	18-25 ms	Operator-leren en voorspelling van de systeemtoestand
ST-GNN	22–35 ms	Spatio-temporele afhankelijkheidsmodellering
RL-beleidnetwerk	5–10 ms	Genereren van regelbeslissingen
Gegevensvoorbehandeling	8–12 ms	Normalisatie van invoer en voorbereiding van kenmerken
Totale latentie van de regelkring	55–82 ms	End-to-end uitvoeringstijd van de regelactie

Tabel 8: Evaluatie van de prestaties van de real-time aansturing. Latentiewaarden vertegenwoordigen de gemiddelde inferentietijden gemeten per aansturingsstap onder een standaard hardwareconfiguratie.

Model	MSE ↓	MAE ↓	Relatieve fout (%) ↓	Inferentietijd (ms) ↓	p-waarde (vs. baseline)
Numerieke oplosser (referentiewaarde)	0.0125 ± 0.0012	0.089 ± 0.006	8.75 ± 0.54	120.5 ± 5.2	—
Conventioneel ML-model	0.0098 ± 0.0009	0.072 ± 0.005	6.42 ± 0.48	85.3 ± 4.7	0.021
Bestaand Operator Learning	0.0076 ± 0.0007	0.061 ± 0.004	5.18 ± 0.41	42.8 ± 3.1	0.008
DeepONet (Voorgesteld)	0.0052 ± 0.0005	0.044 ± 0.003	3.67 ± 0.32	18.6 ± 2.4	< 0.001
FNO (Voorgesteld)	0.0047 ± 0.0004	0.039 ± 0.002	3.21 ± 0.28	15.2 ± 2.1	< 0.001

Tabel 9: Vergelijking van de baseline- en voorgestelde modellen op basis van voorspellingsnauwkeurigheid (MSE, MAE, relatieve fout) en computationele efficiëntie (inferentietijd). Waarden worden gerapporteerd als gemiddelde ± standaarddeviatie over meerdere runs. Lagere waarden duiden op betere prestaties. p-waarden geven de statistische significantie aan ten opzichte van de numerieke baseline-solver.

Methode	MAE ↓	RMSE ↓	F1-score ↑	AUC ↑	Reductie van drift (%) ↑	Cumulatieve beloning ↑
Natuurkundige DT	0.45	0.62	0.78	0.82	35.4	65.2
LSTM (Long Short-Term Memory)	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
Voorgesteld kader	0.2	0.3	0.94	0.97	72.2	95.5

Tabel 10: Uitgebreide prestatievergelijking. Uitgebreide prestatievergelijking van het voorgestelde controle-bewuste digital twin-model met andere basisbenaderingen voor monitoring, anomaliedetectie, synchronisatie en controle. Voor de voorspellingsnauwkeurigheid hebben lagere waarden van MAE en RMSE de voorkeur. Voor de prestaties van anomaliedetectie en synchronisatie hebben hogere waarden van F1-score, AUC, driftreductie en beloning de voorkeur.

Component	Parameter	Waarde/Instelling
Neurale Operatoren (FNO/DeepONet)	Epochs	100–200
	Batchgrootte	32–64
	Leersnelheid	0.001
ST-GNN	Lagen	3
	Verborgen eenheden	64
	Tijdsvenster	10–20 tijdstappen
Neuraal 4D-Var	Assimilatievenster	Vast (schuivend venster)
	Optimalisatie	Gradiëntafdaling
Versterkend leren	Disconteringsfactor (γ)	0.95
	Beloningsgewichten (α,β,γ)	0.6, 0.25, 0.15
Algemene training	Optimalisator	Adam
	Hardware	NVIDIA RTX GPU

Tabel 11: Trainingsconfiguratie en implementatieparameters gebruikt voor verschillende componenten van het voorgestelde digital twin-raamwerk. De tabel geeft een overzicht van de verschillende hyperparameters en trainingsparameters die zijn gebruikt om de reproduceerbaarheid van het voorgestelde systeem voor verschillende componenten te waarborgen.

Component	Trainingstijd (u)	Inferentielatentie (ms)	GPU-benutting (%)	Schaalbaarheid
Neurale Operatoren	6–8	25–40	80–90	Hoog
ST-GNN	3–4	15–25	60–70	Hoog
Neuraal 4D-Var	2–3	20–30	65–75	Matig
RL-controller	1–2	10–20	50–60	Hoog
Algeheel systeem	12–16	<100 ms	70–85	Schaalbaar

Tabel 12: Analyse van de computationele prestaties en schaalbaarheid van het voorgestelde digitale tweelingmodel voor verschillende componenten. Latentie wordt gebruikt om de real-time inferentieprestaties weer te geven, GPU wordt gebruikt om de computationele prestaties weer te geven, en schaalbaarheid staat voor het vermogen van het voorgestelde model om te worden gebruikt voor toepassingen op industriële schaal.

Methode	Monitoring RMSE	Reductie synchronisatiedrift (%)	AUROC	F1-score	Cumulatieve beloning
Op natuurkunde gebaseerde Digital Twin	0.45	28.6	0.78	0.74	42.3
Op LSTM gebaseerd model	0.31	41.2	0.86	0.83	61.7
Op GNN gebaseerd model	0.26	53.4	0.91	0.88	74.5
Voorgestelde Digital Twin	0.2	72.2	0.94	0.92	95.5

Tabel 13: Vergelijkende evaluatie van het voorgestelde werk met bestaande methoden. De prestaties van het voorgestelde digital twin-model op het gebied van monitoring, synchronisatie, anomaliedetectie en controlestatistieken worden vergeleken met baseline- en leergebaseerde benaderingen.

Modelconfiguratie	Monitoringfout (RMSE)	Anomaliedetectie (F1)	Controleprestatie (beloning)
Volledig model (alle modules)	0.2	0.92	95.5
– Zonder data-assimilatie	0.29	0.88	89.2
– Zonder ST-GNN	0.27	0.85	87.1
– Zonder reinforcement learning	0.21	0.9	78.4
– Zonder continu leren	0.23	0.91	92.3

Tabel 14: Experimentele resultaten en ablatiestudie. Ablatiestudie die onderzoekt hoe elk onderdeel van het voorgestelde digitale tweelingmodel bijdraagt aan de algehele prestaties van het systeem.

Scenario	Beschikbaarheid van gegevens	RMSE (root-mean-square error)	Nauwkeurigheid	Omgang met hiaten in observaties
Normaal	100%	Laag	Hoog	Stabiel
Gedeeltelijk verlies	50–70%	Matig	Hoog	Robuuste interpolatie
Ernstig verlies	20–30%	Verhoogd	Matig	Gecontroleerde degradatie
Sensorfout	Regiogewijs 0%	Hoger	Matig	Ruimtelijke generalisatie

Tabel 15: Robuustheidsanalyse bij gegevensschaarste. De volgende tabel geeft de prestaties van het voorgestelde neurale 4D-Var-raamwerk weer in verschillende situaties van gegevensbeschikbaarheid. De tabel geeft de robuustheid van het model aan door te laten zien hoe fouten veranderen in situaties van gedeeltelijk verlies, hoge schaarsheid en totaal verlies van sensorgegevens.

Discussie

De resultaten van dit onderzoek tonen de effectiviteit aan van de ontwikkelde controle-bewuste digitale tweeling wat betreft het vermogen om realtime monitoring, detectie en controle-optimalisatie te verbeteren voor een breed scala aan CPS. Vergeleken met traditionele monitoringmodellen is de combinatie van de PINO en de graafneurale netwerken effectiever wat betreft de nauwkeurigheid van de toestandschatting, zelfs voor een ruisgevoelig systeem. Deze bevinding is correct omdat een geavanceerd machine learning-model in de digitale tweeling effectiever is wat betreft voorspelling en detectie in het systeem³⁹.

In vergelijking met huidige oplossingen voor anomaliedetectie, zoals digital twin-anomaliedetectoren die gebruikmaken van curriculum learning of causale representatiemodellen, is op te merken dat de voorgestelde methode een concurrerende prestatie kan behalen bij vergelijking van de in de literatuur behaalde resultaten. Zo maakten digital twin-gebaseerde anomaliedetectoren die curriculum learning toepasten een bescheiden verbetering van de F1-score mogelijk voor CPS-testbeds, wat aangeeft dat leerparadigma's de efficiëntie van detectieprocessen kunnen beïnvloeden⁴⁰. Tegelijkertijd leggen de huidige populaire digital twin-modellen een sterke nadruk op efficiënte redenering en het onderdrukken van fout-positieven via causale inferentie, maar dit gaat ten koste van een grotere modelcomplexiteit in het causale model. In dit opzicht biedt onze voorgestelde gesynchroniseerde reinforcement learning-benadering zowel een efficiënte operationele optimalisatie als een betrouwbare anomaliedetectie binnen een enkel framework. Er moet echter worden opgemerkt dat de vergelijking van de voorgestelde methode met andere benaderingen niet is gebaseerd op de resultaten van directe vergelijkende experimenten, maar enkel op de door andere onderzoekers gerapporteerde resultaten; daarom dient de vergelijking met deze overweging te worden gemaakt.

Bovendien illustreert de toevoeging van anomaliedetectie en de optimalisatie van regelacties binnen de voorgestelde digitale tweeling het groeiende belang van beslissingsondersteuning en actuatievermogens geboden door digitale tweelingen, een huidige richting van onderzoeks- en ontwikkelingsactiviteiten in verband met CPS. Onderzoeksartikelen⁴¹ met betrekking tot de kalibratie van digitale tweelingenen met behulp van reinforcement learning, zoals adaptieve benaderingen voor het sturen van sensoren, laten zien dat het gebruik van leeragenten extra voordelen kan bieden, met name wat betreft nauwkeurigheid en aanpassingsvermogen, vooral voor dynamische systemen. Dit bevestigt de bevindingen over de verbetering van de effectiviteit van closed-loop optimalisatie, waardoor de hypothese wordt versterkt dat intelligente integratie de robuustheid en efficiëntie van CPS kan verhogen. De benadering is, zoals eerder vermeld, ontworpen voor de olie- en gasindustrie, maar de evaluatie van de benadering is gebaseerd op simulaties van de olie- en gasindustrie, evenals op benchmarkdatasets van CPS, namelijk SWaT en WADI. Er moet echter worden opgemerkt dat de huidige validatie geen operationele gegevens uit de echte olie- en gaswereld bevat, waardoor de toepasbaarheid op scenario's van implementatie in de praktijk voorzichtig moet worden geïnterpreteerd. In de toekomst zal een vergelijkende evaluatie van de voorgestelde curriculum learning-methode en de digitale tweelingmethode op basis van causale inferentie worden uitgevoerd binnen een uniform experimenteel kader om de twee benaderingen op een gecontroleerde manier te vergelijken. Hoewel de integratie van deze modules, waaronder neurale operatoren, graafnetwerken en reinforcement learning, de computationele kosten verhoogt, heeft elke module een specifieke functie. De ablatiestudie valideert daarnaast de prestaties van de voorgestelde benadering, aangezien de prestaties afnemen wanneer een van de modules wordt verwijderd.

Theoretisch gezien levert dit onderzoek aanzienlijke bijdragen aan de toevoeging van modelgebaseerde en datagestuurde methoden binnen de digital twin-structuur. De integratie van physics-informed neural operators en neurale data-assimilatietechnieken in het proces van toestandschatting biedt een hybride modelleringsbenadering voor CPS, waardoor een aanzienlijk gat wordt gedicht in de literatuur over digital twins, die voor CPS-simulaties tot nu toe ofwel vertrouwde op puur natuurkundige simulaties, ofwel op black-box-modellen, maar niet op beide. Theoretisch gezien sluit dit onderzoek aan bij de vooruitgang van digital twin-onderzoek voor autonome CPS-modellen die continu kunnen leren en zich kunnen aanpassen, zoals in recente literatuur is benadrukt.

Wat betreft de toepasbaarheid kan de voorgestelde aanpak de werking van CPS in een industriële setting aanzienlijk verbeteren, aangezien deze real-time betrouwbaarheid, snelle en nauwkeurige anomalie-identificatie en optimale regelstrategieën mogelijk maakt, resulterend in een directe verhoging van de veiligheid en efficiëntie bij de bedrijfsvoering van CPS. Bijvoorbeeld, in een industriële setting waar de omgeving waarschijnlijk gepaard gaat met ruis in sensoren en een niet-lineair systeemmodel, waarborgt de gesynchroniseerde detectieaanpak een snelle identificatie van abnormaliteiten terwijl het voorkomen van valse alarmen wordt onderdrukt. In tegenstelling tot bestaande oplossingen voor intrusiedetectie in CPS-systemen op basis van digital twins, die hoofdzakelijk gericht zijn op classificatiemetrics, maakt de voorgestelde aanpak zowel identificatie- als correctiebenaderingen mogelijk, wat bijdraagt aan het vermogen van systeemoperators om de operationele omgeving van de CPS proactief aan te passen voordat er downtime optreedt. Daarnaast wijst de capaciteit van de voorgestelde aanpak om universeel toegepast te worden in verschillende settings, zoals vastgesteld door de prestaties op de benchmarkdataset, op de algemene toepasbaarheid van de methode in diverse CPS-omgevingen, zoals slimme productie, het energiedomein en beheersystemen voor kritieke infrastructuurnetwerken.

De schaalbaarheid van het voorgestelde framework voor de digitale tweeling wordt beoordeeld op basis van het vermogen om groeiende systeemomvang en temporele complexiteit te beheren. Door sensorinteracties weer te geven met behulp van ijle grafstructuren, schaalt het spatio-temporele graf neurale netwerk efficiënt naarmate het aantal sensoren en gemonitorde componenten toeneemt, waardoor een exponentiële groei in rekencomplexiteit wordt voorkomen. Door operator-mappings te leren onafhankelijk van de discretiseringsresolutie, verbetert de toepassing van physics-informed neurale operators de schaalbaarheid verder en maakt het effectieve modellering van grootschalige systemen mogelijk zonder merkbare toename van de computationele kosten. Het framework maakt gebruik van sliding window-methoden en sequentiegebaseerd leren om temporele afhankelijkheden over grotere tijdshorizons te beheren zonder de computationele levensvatbaarheid in gevaar te brengen. Daarnaast verminderen incrementele leerstechnieken de overhead bij grootschalige implementaties door continue updates mogelijk te maken zonder dat volledige hertraining nodig is. Bovendien wordt effectieve real-time verwerking van hoogdimensionale datastromen mogelijk gemaakt door het gebruik van parallelle GPU-gebaseerde berekeningen. De voorgestelde architectuur is geschikt voor grootschalige industriële toepassingen, zoals olie- en gasproductiesystemen met gedistribueerde sensorennetwerken en continue monitoringvereisten, aangezien deze is ontworpen om op te schalen met zowel het aantal sensoren als de operationele tijd.

De ablatiestudie helpt bij het begrijpen van de individuele bijdrage van elke module, namelijk de neurale operatoren, ST-GNN, data-assimilatie en reinforcement learning. Er is een ablatiestudie uitgevoerd om de bijdrageniveaus van diverse belangrijke individuele modules in deze voorgestelde digitale tweeling-architectuur te meten. Deze individuele modules omvatten data-assimilatie, graafmodellering met ruimte en tijd, reinforcement learning en continu leren. De effecten van het verwijderen van deze individuele modules zijn onderzocht in termen van de effectniveaus die zij creëren voor de bewakingsnauwkeurigheid, anomaliedetectie en controle. De uiteindelijke resultaten valideren duidelijk dat het altijd beter is om gebruik te maken van de volledige architectuur; het verwijderen van individuele modules resulteert dan ook in variërende niveaus van prestatievermindering, wat aangeeft dat verschillende componenten op verschillende wijze bijdragen aan de algehele systeemprestatie. Het effectniveau van het verwijderen van de data-assimilatiemodule laat duidelijk een toename in voorspellingsfouten zien en beïnvloedt de betrouwbaarheid op een gelijkwaardig niveau wat betreft anomaliedetectie. Het effectniveau van het verwijderen van reinforcement learning laat duidelijk zien dat de controle significant is beïnvloed, vooral wanneer er enkel rekening wordt gehouden met bewaking.

De ablatieresultaten in Tabel 14 laten zien dat verschillende modules een verschillende bijdrage leveren aan de prestaties van het systeem. De module voor data-assimilatie heeft de grootste bijdrage aan de nauwkeurigheid van de monitoringprestaties, gevolgd door de ST-GNN-module. De module voor continu leren levert een matige bijdrage aan de nauwkeurigheid van de monitoringprestaties, en de module voor reinforcement learning heeft de kleinste bijdrage aan de nauwkeurigheid van de monitoringprestaties, hoewel deze de grootste bijdrage levert aan de controleprestaties in termen van de daling in de beloningswaarde van 95.5 naar 78.4.

Om de robuustheid van het voorgestelde neurale 4D-Var-raamwerk te valideren, werden experimenten uitgevoerd onder steeds verslechterde observatieomstandigheden. Specifiek omvatten de experimenten de volgende scenario's: (i) gedeeltelijke sensoruitval, waarbij het percentage ontbrekende gegevens varieerde van 30% tot 70%; (ii) volledige sensoruitval in specifieke regio's; en (iii) communicatievertraging die resulteerde in temporeel scharse observaties. De experimentele resultaten tonen aan dat het voorgestelde model een robuuste prestatie vertoont bij matig dataverlies en geleidelijk degradeert bij extreem dataverlies, wat wijst op een hoge mate van aanpassingsvermogen. In Tabel 15 worden de resultaten van de robuustheidsanalyse weergegeven.

Een beperking van deze studie is verbonden aan de reikwijdte van de validatie van het voorgestelde model. Er dient te worden opgemerkt dat hoewel het model bedoeld is voor olie- en gasproductiesystemen, de evaluatie van het voorgestelde model gebaseerd was op gesimuleerde olie- en gaswaarden en cyber-fysieke benchmarkproblemen zoals SWaT en WADI. Het is daarnaast belangrijk op te merken dat hoewel de benchmarkproblemen worden geassocieerd met olie- en gasprocessen, de problemen in dit geval betrekking hebben op waterbehandeling en -distributie.

Ondanks deze optimale prestaties zijn er echter enkele beperkingen. Zo is de computationele complexiteit hoger in vergelijking met sommige implementaties van digitale tweelingen, vooral bij de integratie van neurale operatoren, graafnetwerken en reinforcement learning. Dit zou een grote belemmering kunnen vormen voor de adoptie op platforms met beperkte rekencapaciteit. Bovendien is er uitsluitend gefocust op simulaties en benchmarksets, die niet de mate van diversiteit kunnen bevatten die in een operationeel CPS zou worden aangetroffen. Dit is te zien in simulatiesets zoals SWAT en WADI, die over het algemeen worden gebruikt bij CPS-beoordelingen. Deze kunnen de diversiteit missen die in echte CPS-omgevingen wordt ervaren. Er wordt daarnaast uitgegaan van een adequate aanvoer van sensorgegevens, wat bij oudere CPS niet altijd haalbaar is. De belangrijkste beperking van het systeem is het gebrek aan real-world data uit de olie- en gasindustrie om het systeem te evalueren, aangezien de evaluatie is gebaseerd op simulaties van de olie- en gasindustrie en op benchmarkdatasets van CPS, namelijk SWaT en WADI, die de real-world complexiteiten van de olie- en gasindustrie mogelijk niet nauwkeurig weerspiegelen. Het voorgestelde systeem heeft veelbelovende prestaties getoond in voorspellende vermogens. Het voorgestelde systeem vereist echter hogere computationele kosten tijdens het trainingsproces. In de toekomst kunnen lichtgewicht alternatieven voor het voorgestelde systeem worden overwogen, waaronder reduced-order neurale operatoren, model pruning en hybride physics-ML-modellen.

Schaalbaarheid en implementatie in een praktijomgeving, en de integratie in bestaande legacy-CPS-infrastructuren waarbij de sensordichtheid en de kwaliteit van de waargenomen informatie sterk kunnen variëren, moeten in toekomstig onderzoek naar modelcompressietechnieken worden bestudeerd. Bovendien zou de uitbreiding van dit systeem om uitlegbaarheid en feedbackloops voor operators gunstig kunnen zijn voor het verbeteren van de transparantie en het opbouwen van vertrouwen in autonome besluitvorming. Daarnaast zou de integratie van causale inferentiemethoden binnen de huidige processen voor anomaliedetectie, vergelijkbaar met die in causale digitale tweelingen, kunnen bijdragen aan een verbeterde root-cause-analyse en het minimaliseren van valse alarmen. Ten slotte zal dit raamwerk, om robuustheid te creëren, gevalideerd moeten worden met werkelijke veldgegevens.

Dit onderzoek heeft het concept van een intelligent, controlebewust digital twin-model gedemonstreerd dat erop gericht is om real-time systeemmonitoring, anomaliedetectie en besluitvormingsmogelijkheden in complexe cyber-fysische systemen te verbeteren. In tegenstelling tot traditionele digital twins, die voornamelijk passieve waarnemers zijn en worden gebruikt voor real-time monitoring, beoogt dit onderzoek de integratie van physics-informed neural operators, spatio-temporele graph neural networks, neurale data-assimilatie, reinforcement learning-gebaseerde controle en leeralgoritmen binnen één enkel model van een digital twin-systeem. Het systeem wordt eerst in real-time gemonitord door een digital twin-systeem op basis van een physics-informed leermodel van een niet-lineair stelsel vergelijkingen. Ten tweede wordt een neuraal 4DVar-systeem toegepast voor een data-assimilatieprocedure die zowel de fysische als de virtuele systemen synchroniseert om de drift van systeempredicties als gevolg van ruis en systeemveranderingen te minimaliseren. Anomaliedetectie wordt uitgevoerd op een spatio-temporeel model van ofwel de ruimtelijke relaties of de systeemevolutie over tijd van een systeem. Ten slotte creëert een reinforcement learning-agent systeemcontroleacties vanuit een gesynchroniseerde systeemtoestand die geldig is vanuit het perspectief van een systeemveiligheidsbeperking, terwijl een controlebeperking wordt meegewogen. Uitgebreide experimenten uitgevoerd op een gesimuleerd olie- en gasproces en de benchmark CPS-scenario's SWAT en WADI demonstreerden de potentiële effectiviteit van de voorgestelde aanpak onder gesimuleerde en benchmarkcondities. Het voorgestelde model resulteerde in verminderde monitoringsfouten (MAE en RMSE), uitstekende resultaten bij anomaliedetectie en een consistent gedrag voor closed-loop controle. Deze bevindingen hebben bevestigd dat de synchronisatie van anomaliedetectie en controle binnen de digital twin een belangrijke factor is in relatie tot het behoud van de betrouwbaarheid en efficiëntie van de operaties. Een ander gebied voor potentiële ontwikkeling van digital twin-technologie is de uitbreiding hiervan met kunstmatige intelligentie en causaliteit. De evaluatie van het framework zal worden uitgevoerd door het framework te valideren met real-world data uit de olie- en gasindustrie.

Het is belangrijk om te benadrukken dat de conclusies in deze studie zijn gebaseerd op gesimuleerde en benchmark-datasets, en dat de toepasbaarheid op echte olie- en gassystemen daarom met voorzichtigheid moet worden geïnterpreteerd.

DATAbeschikbaarheid:

De SWaT- en WADI-datasets die in deze studie zijn gebruikt, zijn openbaar beschikbare datasets van industriële controlesystemen en kunnen worden geraadpleegd via de datasetrepository van het iTrust Lab: https://itrust.sutd.edu.sg/itrust-labs_datasets/dataset_info/.

Materialen

Lijst van materialen gebruikt in dit artikel

Naam	Bedrijf	Catalogusnummer	Opmerkingen
CPU (Intel Xeon / AMD Ryzen)	Intel / AMD	Xeon / Ryzen-serie	Gebruikt voor systeemsimulatie, data-preprocessering en algemene berekeningen
Deep Learning Framework (PyTorch 2.0)	Meta AI (PyTorch)	https://pytorch.org RRID:SCR_018536	Gebruikt voor het implementeren van neurale netwerken, optimalisatie en trainingspipelines
Deep Neural Network Library (cuDNN 8.x)	NVIDIA cuDNN	cuDNN 8.x	Biedt geoptimaliseerde primitieven voor deep learning-operaties op de GPU
DeepONet	Niet van toepassing	https://arxiv.org/abs/1910.03193	In het script wordt hier naar verwezen als "neural operator model"
GPU (NVIDIA RTX 4090)	NVIDIA	RTX 4090	Gebruikt voor het trainen van deep learning-modellen, inclusief neurale operatoren, ST-GNN en componenten voor reinforcement learning
GPU-acceleratie (CUDA 11.8)	NVIDIA CUDA	CUDA 11.8 RRID:SCR_018131	Maakt GPU-gebaseerde parallelle berekeningen mogelijk voor high-performance modeltraining
Linux-distributie (versie 22.04 LTS)	Canonical (Ubuntu)	https://ubuntu.com/download; RRID:SCR_018127	Besturingssysteem gebruikt voor het uitvoeren van experimenten en modelimplementatie
Besturingssysteem (Ubuntu 22.04 LTS)	Canonical (Ubuntu)	https://ubuntu.com/download RRID:SCR_018127	Biedt de uitvoeringsomgeving voor alle experimenten en modelimplementatie
Optimalisatiealgoritme (Adam Optimizer)	Niet van toepassing	https://arxiv.org/abs/1412.6980	Gebruikt voor het trainen van modellen met adaptieve leer-snelheid optimalisatie
geoptimaliseerde deep neural network library (versie 8.x)	NVIDIA cuDNN	https://developer.nvidia.com/cudnn RRID:SCR_018131	Biedt geoptimaliseerde, GPU-versnelde primitieven voor deep neural network-operaties
Programmeertaal (Python 3.10)	Python Software Foundation	https://www.python.org/downloads RRID:SCR_008394	Kernprogrammeertaal gebruikt voor implementatie en experimentatie
Systeemgeheugen (≥ 64 GB RAM)	Corsair / Kingston / Crucial	https://www.corsair.com / https://www.kingston.com / https://www.crucial.com	Ondersteunt het verwerken van hoogdimensionale multivariate tijdreeksdata en grootschalige modeltraining