Deze studie omvat geen menselijke deelnemers of dierproeven. De gebruikte datasets (SWaT en WADI) zijn publiek toegankelijke benchmark-datasets, en er is geen sprake van gevoelige of persoonlijke gegevens.
De voorgestelde aanpak realiseert een controlebewuste digitale tweeling waarmee de integratie van geavanceerde deep learning-algoritmen, data-assimilatie en continue evaluatie en optimalisatie van de olie- en gasextractiesystemen en -besturing mogelijk wordt gemaakt. Ten eerste worden real-time multi-modale 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 operatorniveau van de bepalende vergelijkingen voor niet-lineaire meerfasestroming, wat schaalbare simulaties 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 graafneurale netwerken getraind binnen de physics-informed representaties. Een differentieerbare data-assimilatiebenadering op basis van neurale 4DVAR wordt toegepast om de simulatievoorspellingen van de digitale tweeling consistent en gesynchroniseerd te houden met real-time monitoringsgegevens, zelfs bij ruisgevoelige en ruimtelijk onregelmatige metingen. Bovenop de real-time gesynchroniseerde digitale tweeling wordt een door reinforcement learning aangestuurde optimalisatietool getraind voor real-time besturing. Ten slotte blijft de digitale tweeling uitgerust voor het continu bijwerken van modellen en aangepaste schaalbaarheid over een langere periode, gebaseerd op continue updates vanuit nieuwe beschikbare metingen.
De voorgestelde digitale tweeling-architectuur verbetert de conventionele digitale tweeling-architectuur door deze aan te vullen met intelligentie, adaptiviteit en controlebewustzijn. 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 graph neurale netwerken, neurale 4D-Var synchronisatie, physics-informed neurale operatoren en op reinforcement learning gebaseerde closed-loop controle 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 controle-bewuste digitale tweeling voor olie- en gasproductie als een gesloten-lus cyber-fysiek systeem. Het functionele proces begint op het niveau van het olie- en gasproductiesysteem, waarbij verschillende sensoren die verspreid zijn over putten, pijpleidingen, kleppen en verwerkingseenheden op een specifiek tijdstip gegevens vastleggen over druk, debiet, temperatuur en besturingen. De diverse sensoren leggen gegevens vast in verschillende formaten, en deze gegevens worden naar de digitale tweeling gestuurd, wat de basis vormt voor de organisatiebewustheid op dat moment.
De voorgestelde methodologie omvat de volgende stappen:
Gegevensverzameling
Real-time informatie van sensoren met betrekking tot druk, debiet, temperatuur en regelingingangen wordt verzameld uit olie- en gasproductieomgevingen. Andere beschikbare benchmark-datasets, zoals SWAT en WADI, worden gebruikt voor de robuustheid en generalisatievermogens. Om reproduceerbaarheid te waarborgen, wordt elke stap van de voorgestelde methodologie gespecificeerd met technische details.
[Raadpleeg Aanvullend Bestand (S1) voor de generatie van de dataset en de simulatie.]
Twee openbare benchmarkdatasets voor industriële controlesystemen uit het onderzoeksdomein voor industriële controlesystemen worden gebruikt om de capaciteit van de conceptuele digitale tweeling aan te tonen met betrekking tot complexe en hoogdimensionale tijdreeksgegevens, evenals ongebruikelijke gebeurtenissen: de Secure Water Treatment (SWaT) en de Water Distribution (WADI) benchmarkdatasets34. Beide datasets zijn afgeleid van realistische cyber-fysieke testopstellingen en bevatten multivariate sensor- en actuatorbedieningsgegevens voor zowel normale als anomale/aanvalsscenario's die op dergelijke systemen zijn geactiveerd. Door gebruik te maken van de SWaT- en WADI-datasets is het mogelijk om cross-domeintesting uit te voeren en de generaliseerbaarheid aan te tonen voorbij cyber-fysieke testopstellingen naar bredere gedistribueerde industriële systemen met locatie- en tijdsgebonden afhankelijkheden. De beschrijving van de dataset wordt weergegeven in Tabel 2.
Het voorgestelde onderzoek maakt gebruik van drie verschillende datasets voor de uitgebreide evaluatie van de ontwikkelde digital twin-methodologie, zoals weergegeven in Tabel 2. De basisdataset bestaat uit gesimuleerde gegevens met betrekking tot de productieorganisatie in de olie- en gasindustrie en is bedoeld om zowel echte operationele metingen die door het systeem zijn waargenomen als regel- 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 over een periode van ongeveer 11 dagen, inclusief zowel normale als kunstmatig gecreëerde anomalieën en aanvallen. Deze dataset wordt algemeen beschouwd als de standaardbenchmark 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 bovengenoemde datasets concentreren zich voornamelijk op hoogdimensionale tijdreeksen en ruimtelijke relaties, zoals gebruikelijk bij industriële controlesystemen.
Preprocessing van gegevens
Vóór de training van het model ondergaan alle datasets een grondige voorbehandelingsprocedure om de datakwaliteit, consistentie en reproduceerbaarheid te garanderen. De voorbewerkingsstappen omvatten onder meer het beheer van ontbrekende waarden, ruisfiltering, normalisatie en feature scaling. De gedetailleerde voorbewerkingspipeline wordt 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 systeemtoestand in de loop van de tijd (zie het Supplementary File (S3) voor de formuleringen voor de constructie van de systeemgraaf).
Fysica leren met behulp van neurale operatoren
Voor een correcte weergave van de niet-lineaire en multidimensionale kenmerken van meerfasenstroomsystemen maakt de voorgestelde aanpak gebruik van natuurkundig geïnformeerde operatoren in de vorm van Fourier Neural Operators (FNO) en een Neural operator-model. In tegenstelling tot de traditionele methode waarbij neurale netwerken worden weergegeven als puntgewijze afbeeldingen, 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 neural operators verder te evalueren, werden aanvullende experimenten uitgevoerd om hun vermogen te beoordelen om niet-lineaire operatorafbeeldingen te benaderen onder verschillende stromingstoestanden en regelingsinputs. De resultaten tonen aan dat de operators de afbeelding van de niet-lineaire operator nauwkeurig kunnen benaderen, waarmee het vermogen van de operators om de dynamica van het systeem nauwkeurig te modelleren wordt gevalideerd, vooral bij een hoge dimensionaliteit en 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 de voorspelde en gesimuleerde systeemtoestanden, wat een consistente benaderingsprestatie over verschillende stromingscondities aantoont.

Figuur 2: Workflow van het physics-informed neural operator-framework. Een voorbeeld van de physics-informed neural operator-workflow demonstreert hoe regelinputs, sensorgegevens en randvoorwaarden worden gebruikt om niet-lineaire meerfasestromingsdynamica aan te leren en snelle systeemvoorspellingen 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 multifasestromingen 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 is in de vorm van randvoorwaarden, variërend van druk, temperatuur en snelheid, evenals realtime gegevens 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 belangrijkste onderdeel 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 deze neurale operatoren 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 gestuurd door fysische beperkingen om te garanderen dat de voorspelde dynamiek zich gedraagt zoals verwacht in de multifasestroom. Het resultaat van de neurale operator is de continue veldoplossing van de belangrijke fysische variabelen, zoals druk, debieten en fasedistributie voor het gehele productienetwerk. Bovendien zijn de outputwaarden nauwkeurig voor het gehele veld, inclusief de putten, pijpleidingen en andere verbonden 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, welke extracties zijn van olie- en gasproductieprocessen. De kenmerken bestaan grotendeels uit de fysische toestandsvariabelen voor olie- en gasproductie, wat continue velden zijn: temperatuur, debiet, druk en fasefracties voor olie, gas en water. De kenmerken omvatten ook randvoorwaarden, waaronder inlaatdrukken, debieten en klepposities. Besturingsinputs, zoals posities van verstikkleppen en pompsnelheid, worden eveneens meegenomen om responsen onder verschillende operationele scenario's te voorspellen. Het neurale operator-model ontvangt de fysische kenmerken en besturingsparameters, die worden verwerkt door de trunk- en branch-netwerken. De kenmerken leveren voldoende informatie om de neurale operators in staat te stellen 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)36 gebruikt om effectief de operator-mappings te schatten die de niet-lineaire multifase-stroomdynamica 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 verschillende schalen en configuraties van een systeem mogelijk wordt. In dit geval werkt de FNO door de systeemtoestandsrepresentaties 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 input-toestandsrepresentatie geprojecteerd in een hoger-dimensionale latente ruimte, waarna transformaties en activaties iteratief worden toegepast met behulp van Fourier-transformaties, spectrale kernels en activaties. Raadpleeg het Aanvullend Bestand (S4.1) voor de wiskundige evaluaties van de Fourier neural operator.
Neuraal operator-model
Het neurale operator-model37 dient als een alternatieve benadering binnen het neurale operator-proces voor het benaderen van een niet-lineaire fysische operator die de dynamiek van olie- en gasproductie definieert. De gedetailleerde beschrijving van het neurale operator-model staat in Aanvullend bestand (S4.2) (Raadpleeg Aanvullend bestand S10 voor Algoritme 1: Physics Learning with Neural Operators for Digital Twin Algorithm).
De voorgestelde aanpak waarbij neurale operatoren worden gebruikt, begint daarom met het verzamelen van gepaarde voorbeelden van inputs en outputs, waarbij de inputs voorbeelden zijn van systeemcondities in de vorm van systeemgrenswaarden, controlevariabelen of fysieke kenmerken. De neurale operator wordt getraind in een iteratieve fase waarin responsen op een inputfunctie worden voorspeld, met een daaropvolgende verbetering van de afwijkingen tussen voorspellingen en werkelijke responsen met behulp van een mean squared error-functie. Deze aanpak leidt tot optimalisatie op basis van gradiënten, waarbij een neurale operator systemische fysica-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 architectuur van de digitale tweeling38 is het olie- en gasproductieproces weergegeven als een spatio-temporele graaf. Dit is gedaan om zowel de connectiviteit als het tijdsverloop van de gedistribueerde productiemiddelen te conceptualiseren.
Het olie- en gasproductieproces 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 in dit artikel ontwikkelde digital twin-benadering gebruik van spatio-temporele graph neural networks (ST-GNN). Bij deze benadering worden systeementiteiten weergegeven als graafknopen en worden ruimtelijke/functionele systeemrelaties aangeduid door graafranden. Temporele relaties worden uitgedrukt in termen van windowing-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 sensorwaarden. In de praktijk maken 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 onvermijdelijk een verschil tussen voorspellingen en metingen. Om die reden wordt in dit werk, om voorspellingen met feitelijke metingen te synchroniseren, een andere methode geïntroduceerd die bekend staat als 'data-assimilatie' om voorspellingen consistent te maken met fysieke metingen. Om synchronisatie van het fysieke proces en de digitale tweeling te waarborgen tijdens ruisgevoelige en schaarse sensorwaarnemingen, 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 het mogelijk om de toestand van de digitale tweeling online te corrigeren (Raadpleeg Supplementary File (S6) voor Neurale 4D-Var Data-assimilatie, Supplementary File (S10) voor Algoritme 2: Synchronisatie van Digitale Tweelingen middels Neurale 4D-Var en raadpleeg Supplementary File (S11) voor hyperparameters en trainingsdetails).
Detectie en diagnose van anomalieën
In de gesimuleerde olie- en gasdataset werden anomalieën geïnduceerd door gecontroleerde perturbaties te introduceren, zoals drukval over 15%, abnormale flowpieken en klefdefecten. 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 (zie Aanvullend bestand, S7).
Closed-loop optimalisatie en controle
In het voorgestelde systeem kan gesloten-lusregeling worden gerealiseerd door de harmonieuze DT en de module voor reinforcement learning voor optimalisatie te integreren. In het gedeelte van 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 werkelijke proces worden gerealiseerd. Door gebruik te maken van de voorspellingen van de digitale tweeling kan de module voor reinforcement learning dus 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 hebben te maken met niet-stationaire omgevingen met verouderende apparatuur, tijdsvariërende bedrijfsomstandigheden en onvoorziene verstoringen. Om hiervoor een aanpak te ontwerpen, bevat het voorgestelde digital twin-model een continu lerend component dat neurale operatoren, grafmodellen en beleidsregels upgrade op basis van incrementeel leren. Na het waarnemen van grote variaties en niet-geïdentificeerde patronen wordt het leerproces gefaciliteerd door nieuwe gegevens te incorporeren en parameters geassocieerd met het digital twin-model aan te passen.
Continu leren verbetert ook de robuustheid van zowel de regelingsstrategieë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 regelingsresultaten wordt toegepast om modellen aan te passen en een reinforcement learning-strategie te optimaliseren, wat effectieve besluitvormingsprocessen waarborgt terwijl het systeem blijft veranderen en evolueren. Dit adaptieve leerproces transformeert een digitale tweeling van een passieve waarnemer naar een intelligent systeem dat real-time optimalisatieprocessen van complexe olie- en gasproductiesystemen over een langere periode kan ondersteunen, gebruikmakend van een rigoureuze en effectieve technologische infrastructuromgeving. In Tabel 3 worden de modelarchitectuur en de hyperparameters weergegeven.
Raadpleeg Aanvullend Bestand (S9) voor de wiskundige formulering van de Continuous Learning Adaptation.
De module voor datapreprocessing 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 sensormetingen, verwijdert ruis en standaardiseert kenmerken om een consistente datakwaliteit te waarborgen. De verwerkte gegevens worden, samen met randvoorwaarden en regelinputs, vervolgens aangeboden als input aan de physics-informed neural operator module (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 regeloptimalisatie gebruikt om regelstrategieën op een dynamische manier bij te werken. De geoptimaliseerde regeloutputs worden vervolgens doorgegeven aan de besluitvormingsmodule, die de resultaten interpreteert en actiegerichte beslissingen genereert voor de systeemwerking. Er is een feedbackmechanisme geïntegreerd om het systeem continu bij te werken op basis van nieuw geobserveerde gegevens, waardoor de aanpassingscapaciteit 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 datapreprocessing, physics-informed neural operator modellering, regeloptimalisatie en besluitvorming. Deze modules zijn met elkaar verbonden via goed gedefinieerde input-outputinterfaces, wat zorgt voor een naadloze datastroom, iteratief leren en gecoördineerde werking over het gehele digital twin-systeem.

Figuur 3: Stroomschema van de architectuur voor optimalisatie van de besturing op systeemniveau. Het toont de algehele systeemarchitectuur van het voorgestelde werk. Klik hier om een grotere versie van deze figuur te bekijken.
De systeemarchitectuur geeft een volledig beeld van het model dat kan worden gebruikt voor het modelleren en optimaliseren van complexe niet-lineaire fysieke dynamiek in olie- en gasproductiesystemen. Dit begint bij de inputgegevens, die de parameters voor de simulatie en de randvoorwaarden die de fysieke omgeving representeren, omvatten. Deze gegevens passeren eerst de module voor data-preprocessing, waar normalisatie en feature-extractie worden uitgevoerd. Dit bereidt de gegevens voor op de volgende module, waarin de neurale operator, FNO, en het neurale operator-model worden gebruikt. Deze modellen worden aangewend om de systeemtoestanden te voorspellen, en deze informatie wordt gebruikt om de systeemparameters te optimaliseren middels 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 outputresultaten, waarbij de resultaten van het optimalisatieproces worden getoond in Figuur 3. Daarnaast is er een feedbackloop opgenomen, waarbij de gegevens na de outputresultaten opnieuw via de preprocessingmodule lopen.