$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Fizik tabanlı sinirsel operatörlerin geçerliliğini test etmek amacıyla, doğrusal olmayan çok fazlı akış sisteminin dinamiğini modellemede FNO ile sinirsel operatör modelinin performansını karşılaştırmak için ek performans testleri gerçekleştirilmiştir. Sonuçlar, FNO'nun küresel uzaysal tutarlılıkta daha iyi performans gösterdiğini, sinirsel operatör modelinin ise heterojen girdilere uyumlulukta daha başarılı olduğunu göstermiştir.
Önerilen dijital ikiz çerçevesinin etkili bir şekilde eğitilmesini ve değerlendirilmesini sağlamak amacıyla tüm deneyler yüksek performanslı bir hesaplama ortamında gerçekleştirilmiştir. Uygulama için temel programlama dili olarak Python kullanılmıştır. Fizik tabanlı sinir operatörleri, uzaysal-zamansal grafik sinir ağları ve pekiştirmeli öğrenme modülleri gibi derin öğrenme modellerini oluşturmak için CUDA ile GPU hızlandırması etkinleştirilmiş popüler bir derin öğrenme çerçevesi ve optimize edilmiş bir derin sinir ağı kütüphanesi kullanılmıştır. Büyük ölçekli çok değişkenli zaman serisi veri işlemini ve model eğitimini kolaylaştırmak için testler; çok çekirdekli bir CPU, özel belleğe sahip yüksek performanslı bir grafik işlem birimi (GPU) ve en az 64 GB sistem belleği içeren bir iş istasyonunda yürütülmüştür. Derin öğrenme ekosistemiyle uyumluluğu sağlamak için işletim sistemi ortamı bir Linux dağıtımı üzerine kurulmuştur.
Modelleri eğitmek için mini-batch gradyan tabanlı optimizasyon kullanılmış ve istatistiksel güvenilirliği garanti altına almak adına deneylerin çok sayıda bağımsız çalışması gerçekleştirilmiştir. Sağlanan hesaplama kurulumu sayesinde, önerilen metodoloji benzer yüksek performanslı hesaplama ortamlarında tekrarlanabilir.
Sonuçların tekrarlanan deneyler yoluyla doğrulanmasını sağlamak için, tüm deneyler farklı rastgele başlatma değerleri ile birden fazla bağımsız çalışma ile tekrarlanmıştır. Simüle edilmiş verilerin kullanıldığı deneylerde, model (N) bağımsız deneme için çalıştırılmış ve değerler ortalama performans ve standart sapma ile raporlanmıştır. SWAT ve WADI karşılaştırma verilerinin kullanıldığı deneylerde ise, sonuçları istatistiksel olarak güçlü kılmak için MAE, RMSE, F1 Score ve AUROC değerlendirme metrikleri test verileri üzerinde hesaplanmış ve tekrarlanan deneyler üzerinden ortalaması alınmıştır. Bu işlem, modelin performansına dair güvenilir bir tahmin elde etmek ve performans artışlarının tek bir denemeden kaynaklanmadığından emin olmak için yapılmıştır.
Buna rağmen, fizik bilgisine dayalı öğrenme, uzaysal-zamansal modelleme ve veri asimilasyon tekniklerinin entegre bir dijital ikiz sisteminde birleştirilmesinin, geleneksel yöntemlerle karşılaştırıldığında, değerlendirilen koşullar altında tutarlı bir şekilde daha iyi izleme performansı, daha hızlı yakınsama ve gerçek zamanlı optimal kararlar sağladığı doğrulanmıştır. Tablo 4, önerilen çalışmanın simülasyon ortamını göstermektedir.
Simülasyon ortamı, endüstrideki gerçek zamanlı işletimi doğru şekilde yansıtacak şekilde modellenmiştir. Akış veri ortamını temsil etmek için sensörlerden gelen veriler kısa aralıklarla alınmış ve dayanıklılığı test etmek amacıyla kontrollü miktarlarda gürültü eklenmiştir. Anomali tespiti ve kontrol yeteneklerini test etmek için çeşitli hata ve bozulma durumları simüle edilmiştir. Değerlendirme için kullanılan kayan pencere yaklaşımı; tahminler, durum senkronizasyonları ve kontrol eylemlerinin üretilmesi için önerilen dijital ikizin dinamik güncellemelerine uygun şekilde, gerçek zamanlı işletimin sürekli bir biçimde temsil edilmesini sağlar.
Kontrol farkındalıklı dijital ikiz modelinin tasarımında, eş zamanlı değerlendirme için dikkate alınması gereken çeşitli değerlendirme kriterleri kullanılır. İlk olarak, dijital ikiz mimarisinin izleme ve durum tahmin gücü için Ortalama Mutlak Hata (MAE), Kök Ortalama Kare Hata (RMSE) ve Ortalama Mutlak Yüzde Hata (MAPE) gibi regresyon kriterleri kullanılır. MAE ve MAPE; basınç, akış hızı ve sıcaklık gibi temel sistem parametreleri üzerinde nöral operatörlerin ve uzamsal-zamansal grafik modellerinin doğruluğunu hesaplar ve petrol ile gaz üretim sürecinde kullanılan sensörlerin sürekli çok değişkenli verileri için geçerlidir.
İkinci olarak, durum kestirim hatası ve tahmin kayması giderimiyle ilgili olarak dijital ikizin kendi senkronizasyon süreci de değerlendirilir. Durum kestirim hatası, gerçek şema durumu ile senkronize edilmiş dijital ikiz durumu arasındaki hata derecesini ölçer; bu, geliştirilen nöral 4D-Var modülünün, sistem dinamiği tahminlerini iyileştirmek için gerçek sistem durumu gözlemlerini ne kadar etkili kullandığının bir göstergesidir. Son olarak, tahmin kayması azalması, önerilen yaklaşımın gürültülü veya seyrek sistem gözlemlerine dayanma yeteneği sayesinde, senkronize edilmemiş dijital ikiz sistemine karşı elde ettiği performans kazancını gösterir.
Üçüncü olarak, anomali tespiti ve hata teşhisine gelince; Kesinlik (Precision), Duyarlılık (Recall), F1 ve ROC Eğrisi Altındaki Alan (AUC) gibi sınıflandırma odaklı değerlendirme parametreleri kullanılır. Duyarlılık, herhangi birini gözden kaçırmadan gerçek hataları doğru şekilde teşhis edebilme yeteneğini ifade ederken; kesinlik, tespit edilen bir anomaliyle ilgili geçerliliğin bir ölçüsü olarak görülebilir. Performans, uygun bir şekilde F1 metriği kullanılarak ölçülebilir. Değerlendirme parametreleri genel olarak bir SWAT ve WADI sisteminde benimsenmektedir.
Tablo 5, önerilen nöral 4D-Var senkronizasyon modülünün fiziksel ve siber modeller arasındaki tahmin kaymasını azaltmadaki etkinliğini üç veri seti için değerlendirmek amacıyla kullanılmıştır. Başlangıç durum hatası, senkronizasyon uygulanmamış modeller arasındaki hataya karşılık gelirken; nihai durum hatası, senkronizasyon uygulandıktan sonraki model hatasına karşılık gelmektedir. Her üç veri seti için de durum hatasında belirgin bir iyileşme olduğu görülmektedir; en yüksek iyileşme (%72,2), SWAT (%63,6) ve WADI (%57,7) ile karşılaştırıldığında göreceli basitliği nedeniyle simüle edilmiş petrol ve gaz sistemi modeli için gerçekleşmiştir. Bu okumalar, ilgili karmaşıklık ve "gürültü" seviyeleriyle tutarlıdır. Simüle edilmiş petrol ve gaz veri setindeki kayma azalma oranının daha yüksek olması, simüle edilmiş ortamın SWAT ve WADI gibi diğer gerçekçi siber-fiziksel sistemlerle karşılaştırıldığında genellikle daha düşük gürültü seviyesi, daha az rastgele pertürbasyon ve iyi tanımlanmış sistem davranışı içermesiyle açıklanabilir. Ancak, bu argüman bir hipotez olarak değerlendirilmeli ve bu çalışmanın sonuçlarına dayanarak nicel olarak kanıtlanmamalıdır. Gelecek çalışmalarda, bu argümanı doğrulamak için entropi, gürültü varyansı ve sistem boyutu gibi karmaşıklık ölçütleri kullanılabilir.

Şekil 4: MAE, RMSE ve MAPE (%) hata metrikleri kullanılarak değerlendirilen farklı veri setleri (Simüle edilmiş Petrol ve Gaz, SWAT ve WADI) genelindeki performans izleme. Çubuklar, birden fazla bağımsız çalışmadan (n = 5) elde edilen ortalama değerleri temsil etmektedir. Hata çubukları standart sapmayı (SD) gösterirken, üzerine bindirilmiş daha ince hata çubukları ortalamanın standart hatasını (SEM = SD/√n) temsil etmektedir. Her metrik için SD ve SEM değerleri lejantta açıkça gösterilmiştir. İstatistiksel anlamlılık, baz veri setine (Simüle edilmiş Petrol ve Gaz) göre bağımsız t-testleri kullanılarak değerlendirilmiş ve ilgili p-değerleri çubukların üzerinde gösterilmiştir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Önerilen dijital ikiz modelinin Simulated Oil and Gas, SWAT ve WADI gibi veri kümeleri üzerinde gerçekleştirdiği izleme sürecinin performans analizi; MAE, RMSE ve MAPE değerlendirme kriterleri dikkate alınarak Şekil 4'te sunulmuştur. Simüle edilmiş Petrol ve Gaz (Simulated Oil and Gas) veri kümesi, tüm kriterler için en düşük hata ölçümüne sahiptir; bu durum, sistem parametrelerinin hassas şekilde modellendiği simüle bir ortamda sistem durum izlemenin doğruluğunu açıkça ortaya koymaktadır. SWAT ve WADI veri kümeleri, simüle sistemlerde bulunmayan, gerçek sistemlerin karmaşıklıkları ile doğal gürültü ve doğrusal olmama durumları nedeniyle sistem izlemede kademeli olarak artan hatalar sergilemektedir. Buna rağmen, MAE ve RMSE hata ölçümlerindeki küçük artış, dijital ikiz modelinin bu tür koşullar altındaki sistem performans kararlılığını doğrulamaktadır.

Şekil 5: Farklı veri kümeleri (Simüle Petrol ve Gaz, SWAT ve WADI) genelindeki anomali tespit performansının F1-skoru ve AUC kullanılarak değerlendirilmesi. Çubuklar, birden fazla bağımsız çalıştırmadan elde edilen ortalama değerleri temsil eder (n = 5). Hata çubukları SD'yi, üzerindeki daha ince hata çubukları ise SEM = SD/√n değerini göstermektedir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Şekil 5, önerilen anomali keşif sisteminin üç veri seti için F1-skoru ve AUC tabanlı değerlendirme analizini göstermektedir. Simülasyon veri seti olan Petrol ve Gaz veri seti, en yüksek F1-skoru ve AUC değerlerini elde ederek en iyi performansı sergilemiştir. Performansta kademeli bir düşüş olsa da, veri seti karmaşıklığına bağlı olarak küçük farklılıklar göstermekle birlikte tüm veri setleri genelinde yüksek AUC performansı elde edilmiştir. Üç veri seti için hem F1-skoru hem de AUC eğrilerinin uyumu, iyi bir dengenin mevcut olduğunu ortaya koyan optimal bir kesişim noktasını göstermekte ve önerilen sistemi erken hata tespiti ve sızma için güvenilir bir sistem haline getirmektedir. Önerilen sistem, üç veri setinde anomalileri tanıma konusunda üstün performans göstermektedir.
Önerilen dijital ikiz modelinin farklı veri setlerinde anomali tespitindeki kesinlik (precision), geri çağırma (recall) ve F1 skorları aşağıdaki Tablo 6'da gösterilmiştir. Yüksek kesinlik değerleri, belirlenen ilgi noktalarının çoğunun gerçekten doğru ilgi noktaları olduğunu gösterirken, yüksek geri çağırma değerleri ise gerçek ilgi noktalarının çoğunun önemli olanları atlamadan doğru şekilde tanımlandığını göstermektedir. Petrol ve gaz alanı üzerindeki demonstrasyon ortamı 0,94 ile en yüksek F1 skorunu kaydederken, SWAT ve WADI üzerindeki biraz daha düşük sonuçlar, sensörlerdeki ve saldırı modellerindeki daha yüksek gürültü seviyelerine ve karşılıklı bağımlılıklara bağlanabilir. Bununla birlikte, yüksek F1 skorları, anomalileri tespit etmek için sağlam bir uzay-zamansal grafik modellemesi ve senkronizasyonun sağlandığını doğrulamaktadır.

Şekil 6: Simüle Edilmiş Petrol ve Gaz, SWAT ve WADI veri kümelerinde anomali tespiti için Alıcı İşletim Karakteristiği (ROC) eğrileri. Performans, Eğri Altındaki Alan (AUC) kullanılarak değerlendirilmiştir. AUC değerlerinin standart sapması (SD) ve ortalamanın standart hatası (SEM = SD/√n, n = 5), her veri kümesi için lejandda belirtilmiştir. Diyagonal kesikli çizgi, rastgele sınıflandırma performansını temsil etmektedir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Şekil 6, eşik değer seçiminden bağımsız olarak anomali tespit doğruluğunu değerlendiren, üç farklı veri seti için ilişkili AUROC ölçümleriyle birlikte Alıcı İşletim Karakteristiği (ROC) eğrilerini göstermektedir. Simüle edilmiş Petrol ve Gaz veri setinin ROC eğrisinde, bölge başlangıç noktasına oldukça yakındır ve böylece neredeyse kusursuz sınıflandırma yeteneğini doğrulayan 0,99'luk bir AUROC değeri saptanmıştır. SWAT ve WADI veri setleri durumunda AUROC değerleri sırasıyla 0,97 ve 0,95 olup, yüksek belirsizliklerin varlığına rağmen iyi bir sınıflandırma yeteneği olduğunu ortaya koymaktadır. ROC eğrilerinin rastgele sınıflandırmayı temsil eden çizgilerden uzak olması, geliştirilen dijital ikiz modelinin doğru anomali tespiti için şans seviyesinden belirgin şekilde daha iyi performans gösterdiğini garanti etmektedir. Şekil 6'daki Alıcı İşletim Karakteristiği (ROC) eğrileri, önerilen dijital ikiz modelinin üç veri seti üzerindeki anomali tespit performansını göstermektedir. AUROC skorları Petrol ve Gaz simüle veri seti için 0,99, SWAT veri seti için 0,97 ve WADI veri seti için 0,95 olup, modelin hem simüle edilmiş hem de kıyaslama amaçlı endüstriyel kontrol sistemlerindeki ayırt etme yeteneğini göstermektedir.
Tablo 7, dijital ikiz ile entegre edilmiş kapalı döngü kontrol için takviyeli öğrenme modülüne ait kümülatif ödüllerin ve kısıt ihlallerinin kısa bir açıklamasını sunmaktadır. Elde edilen kümülatif ödül, operasyonel hedeflere dayalı genel kontrol performansını temsil ederken, kısıt ihlalleri operasyon sırasındaki toplam ihlal sayısını temsil eder. Sıfır kısıt ihlali ile elde edilen kümülatif ödüller, petrol ve gaz model ortamımız için en yüksek (96,5) seviyededir ve bu durum optimal kontrol ve operasyonel performansı kanıtlamaktadır.
SWaT ve WADI veri setleri biraz daha düşük kümülatif ödüller ve bazı kısıtlama ihlalleri gösterse de, operasyonel maliyet tasarruflarındaki genel iyileşme, önerilen dijital ikiz çerçevesinin dinamik ortamlarda güvenli ve optimal kontrol için etkili olmaya devam ettiğini kanıtlamaktadır.
Önerilen pekiştirmeli öğrenme (RL) tabanlı kontrol sisteminin gerçek zamanlı uygulanabilirliği, sistem içinde kullanılan tüm hesaplama modüllerinin uçtan uca çıkarım gecikmesi analiz edilerek değerlendirilmiştir. Tablo 8'deki deneysel gözlemler; nöral operatör modülünün, FNO/Nöral operatör modelinin ve uzay-zamansal graf nöral ağ (ST-GNN) modülünün, işlem sırasında kullanılan zamansal özellikler nedeniyle orta düzeyde hesaplama yüküne sahip olduğunu göstermektedir. Öte yandan, pekiştirmeli öğrenme politika karar modülünün hesaplama yükü oldukça düşüktür. Önerilen sistem içinde kullanılan tüm modüllerin kümülatif uçtan uca çıkarım süresi, tipik endüstriyel kontrol döngüleri aralığındadır. Spesifik olarak, ortalama uçtan uca çıkarım gecikmesinin kontrol döngüsü başına 50–150 milisaniye arasında olduğu gözlemlenmiş olup bu değer standart endüstriyel kontrol gereksinimleriyle uyumludur. Bu durum, önerilen sistemin gerçek zamanlı ve gerçek zamanlıya yakın uygulamalar için kullanılabileceğini göstermektedir. Bununla birlikte, uçtan uca çıkarım süresinin kullanılan donanıma bağlı olarak değişebileceği unutulmamalıdır. Bu bağlamda, önerilen sistem gerçek zamanlı uygulamalar için kullanılabilir.
Fizik tabanlı sinirsel operatörlerin etkinliğini daha detaylı olarak ortaya koymak amacıyla; FNO, Sinirsel operatör modeli ile geleneksel sayısal çözücüler ve mevcut operatör öğrenme teknikleri gibi temel yaklaşımların nicel bir karşılaştırması gerçekleştirilmiştir. Tablo 9'daki değerlendirme, her bir modelin sistem dinamiğini ve doğrusal olmayan dinamikleri ne kadar iyi yakaladığına odaklanmaktadır. Performansı değerlendirmek için ortalama kare hata (MSE), ortalama mutlak hata (MAE) ve beklenen sistem durumları ile gerçek durumlar arasındaki bağıl hata gibi standart metrikler kullanılmıştır. Deneysel sonuçlara göre, hem FNO hem de Sinirsel operatör modeli, daha düşük hata değerleri elde ederek karmaşık çok fazlı dinamiklerin modellenmesinde temel yaklaşımlardan daha iyi performans göstermiştir. Özellikle, Sinirsel operatör modeli farklı giriş koşulları altında doğrusal olmayan operatör eşlemelerini başarıyla yakalarken, FNO küresel uzaysal bağımlılıkların öğrenilmesinde başarılı olmuştur.
Önerilen nöral operatör modellerinin etkinliği, Tablo 9'da sunulan nicel karşılaştırmalı sonuçlarla gösterilmiştir. Temel sayısal çözücüler ve mevcut operatör öğrenme teknikleriyle karşılaştırıldığında, hem FNO hem de nöral operatör modeli fark edilir derecede daha düşük hata metrikleri elde etmiştir. Özellikle FNO, karmaşık doğrusal olmayan sistemlerin dinamiklerini yakalamadaki üstün kapasitesini kanıtlayarak en düşük bağıl hata ve ortalama kare hata değerlerini sergilemiştir. Ek olarak, nöral operatör modeli, özellikle farklı koşullar altında doğrusal olmayan operatör eşlemelerini öğrenme konusunda önemli kazanımlar göstermektedir. İstatistiksel anlamlılığı değerlendirmek için bağımsız t-testleri kullanılmıştır; 0,05'ten küçük p-değerleri, temel sayısal çözücüye kıyasla önemli bir iyileşmeyi ifade eder. Önerilen her iki model tarafından elde edilen p-değerleri 0,001'den küçüktür; bu durum, performansta gözlemlenen iyileşmelerin istatistiksel olarak anlamlı olduğunu ve tesadüf sonucu olmadığını göstermektedir. Ayrıca, önerilen teknikler önemli ölçüde daha kısa çıkarım süreleri sergileyerek gerçek zamanlı dijital ikiz uygulamalarına uygunluklarını kanıtlamıştır.
Önerilen kontrol farkındalıklı dijital ikiz modelinin deneysel değerlendirmesinin netliğini daha da artırmak amacıyla, önerilen yaklaşımın diğer yöntemlerle kapsamlı nicel karşılaştırması entegre bir tablo formatında sunulmuştur. Ortalama Mutlak Hata (MAE), Kök Ortalama Kare Hata (RMSE), F1-skoru, Eğri Altındaki Alan (AUC), tahmin sapmasının azaltılması ve kümülatif ödül gibi çeşitli performans parametrelerinin bu karşılaştırması; yöntemlerin izleme doğruluğu, anomali tespit yeteneği, senkronizasyon verimliliği ve kontrol performansının entegre bir şekilde değerlendirilmesine Tablo 10'da olanak sağlayacaktır. Ayrıca, deneysel değerlendirmenin şeffaflığını artırmak için grup büyüklüğü, öğrenme oranı, epok sayısı ve modül bazlı eğitim yapılandırması gibi eğitim prosedürünün çeşitli parametreleri Tablo 11'de açıkça belirtilmiştir. Ayrıca, önerilen yaklaşımın çıkarım gecikmesi, GPU kullanımı ve gerçek zamanlı uygulanabilirliği açısından ölçeklenebilirliği ve hesaplama verimliliği de Tablo 12'de değerlendirilmiştir.
Önerilen dijital ikiz; izleme doğruluğu, senkronizasyon yeteneği, anomali tespiti ve kontrol performansı açısından Tablo 13'te gösterildiği gibi fizik tabanlı senkronizasyon yöntemi, LSTM ve GNN modeli gibi diğer yöntemlerle kapsamlı bir şekilde karşılaştırılabilir. Fiziksel dijital ikizler, statik yaklaşımlara dayandıkları için daha az uyarlanabilirdir ve daha yüksek hata payına sahiptir. LSTM ve GNN gibi nöral yaklaşımlar daha iyidir ancak fizik, senkronizasyon ve kontrolü etkili bir şekilde birleştirmemektedir. Önerilen dijital ikiz; 0.20 RMSE, drift azaltmada %72.2, 0.94 AUROC, 0.92 F1-Skoru ve 95.5 kümülatif ödül ile en iyi değerlere ulaşmıştır. Yukarıdaki sonuçlar; fizik tabanlı nöral operatörlerin, zaman ve mekan için graf modellemenin, nöral veri asimilasyonunun ve pekiştirmeli öğrenmenin etkili bir şekilde birleştirilmesinin, doğru, sağlam ve kontrol farkındalığına sahip bir dijital ikiz sistemini önemli ölçüde geliştirebileceğini açıkça göstermektedir.

Şekil 7: FNO, Sinirsel operatör modeli, GNN, RL modülü ve toplam çerçeve dahil olmak üzere önerilen sistemin farklı bileşenlerinin eğitim maliyeti ve kaynak kullanımı. Metrikler eğitim süresini (saat) ve GPU bellek kullanımını (GB) içermektedir. Değerler, birden fazla bağımsız çalışmanın (n = 5) ortalama sonuçlarını temsil etmektedir. Hata çubukları SD değerini, üzerine binen daha ince hata çubukları ise çalışmalar arasındaki değişkenliği yansıtan SEM = SD/√n değerini göstermektedir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.
Şekil 7'de, önerilen yaklaşımın eğitim maliyeti ve kaynak kullanımı ana bileşenleri şeklinde sunulmuştur. FNO ve Nöral operatör modeli gibi nöral operatör modellerinin, genel sisteme kıyasla önemli eğitim maliyetlerine ve GPU bellek kullanımına sahip olduğu gözlemlenebilir. Bu yüksek eğitim maliyeti, operatörlerin yüksek boyutlu öğrenme yeteneğinden kaynaklanmaktadır. Graf nöral ağ ve pekiştirmeli öğrenme modüllerinin daha düşük eğitim maliyetine sahip olduğu, ancak mekânsal bağımlılıkların öğrenilmesi ve adaptif kontrol için önemli oldukları tespit edilmiştir. Önerilen sistem daha yüksek eğitim maliyetlerine sahip olsa da, çıkarım süresinin verimli olduğu bulunmuştur. Bağımsız LSTM veya GNN modelleri gibi temel yöntemlerle karşılaştırıldığında, önerilen sistem çoklu modül entegrasyonu nedeniyle daha yüksek eğitim maliyeti oluşturmakta; ancak Tablo 13'te gösterildiği üzere, geliştirilmiş doğruluk ve sağlamlık sağlamaktadır.
Genel olarak, önerilen dijital ikiz yapısı, her biri karmaşık endüstriyel sistemlerin modelleme ve kontrolündeki farklı bir sorunu ele alan çeşitli elementleri birleştirir. Sıklıkla hesaplama maliyeti yüksek olan ve gerçek zamanlı uygulamalar için daha az uygun olan geleneksel sayısal çözücülerin aksine, doğrusal olmayan ve yüksek boyutlu sistem dinamiklerini daha etkili bir şekilde yakalamak için Fourier Neural Operator ve Deep Operator Network gibi fizik destekli nöral operatörler kullanılır. İleri beslemeli veya tekrarlayan nöral ağlar gibi daha basit makine öğrenimi modelleri kullanılabilir, ancak bunlar genellikle temel fiziksel operatörleri öğrenemezler ve farklı sistem konfigürasyonları genelinde genelleme yapamazlar. Petrol ve gaz üretim sistemlerinde doğası gereği birbirine bağlı olan dağıtılmış sensörler arasındaki bağımlılıklar, uzamsal-zamansal grafik nöral ağ kullanılarak modellenir. Bağımsız zaman serisi modelleri veya basit tekrarlayan mimariler gibi geleneksel yöntemler, uzamsal bağlantıları açıkça yansıtmadıkları için durum kestirimi ve anomali tespitinde daha düşük doğruluğa yol açabilir. Takviyeli öğrenme modülünden yararlanan adaptif kontrol, sistemin dinamik ve belirsiz ortamlarda operasyonel seçimleri maksimize etmesini sağlar. Öte yandan, kural tabanlı veya statik kontrol stratejileri katıdır ve değişen sistem koşullarına veya beklenmedik aksaklıklara uyum sağlayamazlar. Daha az karmaşıklığa sahip alternatifler mevcut olsa da, bunlar genellikle sorunun yalnızca belirli kısımlarıyla ilgilenirler ve hassas modelleme, gerçek zamanlı izleme ve adaptif kontrolü içeren entegre bir çözüm sağlamak için yetersizdirler. Önerilen modüler tasarım, her bir parçanın genel sistem performansına katkıda bulunmasını sağlayarak, zorlu endüstriyel ortamlarda doğruluğu, ölçeklenebilirliği ve dayanıklılığı artırır.
| Çalışma / İş | Kullanılan Teknikler | Sonuçlar | Kısıtlamalar |
| Dijital İkiz için Derin Sinirsel Operatör26 | Vekil model olarak DeepONet sinir operatörü | Gerçek zamanlı tahmin; fizik simülatörlerinden katlarca daha hızlı | Sensör yerleşiminin dikkatle yapılması gerekmektedir; değerlendirme yöntemlerinin geliştirilmesi gerekir |
| DeepONet ile Sanal Algılama Destekli Dijital İkiz27 | DT'de DeepONet sanal sensörleri | Hızlı tahminler; gerçek zamanlı çıkarım; kapsamlı durum kestirimi | Spektral yanlılık yüksek frekanslı fenomenleri etkiler; hibrit modelleme gereklidir |
| Kontrolde DRL + NMPC ile Dijital İkiz28 | Derin Takviyeli Öğrenme & Doğrusal Olmayan MPC entegrasyonu | İyileştirilmiş kontrol performansı; daha düşük izleme hatası ve adaptif ikiz güncellemeleri | Uygulama kaplarla sınırlıdır; üretim sistemlerine daha fazla uygulama yapılması gerekmektedir |
| Endüstriyel Otomasyonda Yapay Zeka Güdümlü Dijital İkiz29 | Gözetimli/Gözetimsiz Makine Öğrenmesi (LSTM, CNN, RF, vb.) | Hata tahmininde yüksek doğruluk & izleme | Öncelikle klasik ML; gelişmiş uzaysal-zamansal derin modellerden yoksundur |
| Yağ & Gaz Dijital İkiz İncelemesi16 | Çeşitli YZ/MO & mekanistik modelleme | Modelleme stratejilerinin kapsamlı taksonomisi | Çoğunlukla kavramsal; ileri düzey derin öğrenme ile minimal pratik uygulamalar |
Tablo 1: Petrol ve gaz üretimi ile kontrol yönetiminin mevcut yöntemlerle karşılaştırmalı analizi. Kullanılan yöntemlere, elde edilen sonuçlara ve önemli dezavantajlara vurgu yapan ilgili dijital ikiz ve kontrol stratejilerinin genel görünümü.
| Veri Seti Adı | Alan / Sistem | Öznitelik Sayısı | Süre & Boyut | Veri Özellikleri |
| Simüle Edilmiş Petrol & Gaz Üretim Verileri | Petrol & Gaz Üretim Sistemi (Simülasyon) | Çok Değişkenli (Basınç, Akış, Sıcaklık, vb.) | Bu araştırma için tasarlanmıştır | Simüle edilmiş gerçek zamanlı sensör & kontrol sinyalleri, normal & hata durumları |
| SWaT | Güvenli Su Arıtma Siber-Fiziksel Sistemi (CPS) | ~51 sensör + aktüatör | ~11 gün (~450k+ örnek) | Normal & anomali (36 saldırı senaryosu) |
| WADI | Su Dağıtım Siber-Fiziksel Sistemi (CPS) | ~123 sensör + aktüatör | ~16 gün (~1M+ örnek) | Normal & anomali (15 saldırı) |
Tablo 2: Veri seti açıklaması.
Veri Seti Açıklaması. Alan, özellik sayısı, süre ve veri öznitelikleri dahil olmak üzere çalışmadaki veri setlerinin özeti.
| Bileşen | Hiperparametre | Değer |
| ST-GNN | Katman sayısı | 3 |
| ST-GNN | Gizli boyut | 64 |
| ST-GNN | Zamansal pencere uzunluğu | 12 |
| Sinyiral Operatör (FNO) | Fourier modlarının sayısı | 16 |
| Sinyiral Operatör (FNO) | Katman sayısı | 4 |
| DeepONet | Dal ağ boyutu | [128, 128] |
| DeepONet | Gövde ağ boyutu | [128, 128] |
| Veri Asimilasyonu | Asimilasyon penceresi | 10 zaman adımı |
| Optimize Edici | Adam | |
| Öğrenme oranı | 0.001 | |
| Yığın boyutu | 64 | |
| RL Denetleyicisi | İndirim faktörü (γ) | 0.99 |
| RL Denetleyicisi | Politika ağı katmanları | [128, 64] |
Tablo 3: Model Mimarisi ve Hiperparametreler. Sinirsel operatörler, veri asimilasyonu, pekiştirmeli öğrenme ve uzay-zamansal grafik sinir ağı bileşenleri için önemli mimari konfigürasyonlar ve hiperparametre ayarları.
| Parametre | Açıklama |
| Simülasyon platformu | PyTorch (v2.1) içeren Python (v3.10) |
| Destekleyici kütüphaneler | NumPy (v1.24), SciPy (v1.10), Pandas (v1.5) |
| Donanım | NVIDIA RTX serisi GPU (örneğin RTX 3090, 24 GB VRAM), Intel Core i7/i9 CPU, 32–128 GB RAM |
| GPU hızlandırma | CUDA Toolkit (v11.8), cuDNN (v8.x) |
| İşletim sistemi | Linux (Ubuntu 20.04) / Windows 10 |
| Geliştirme ortamı | Jupyter Notebook / VS Code |
| Örnekleme aralığı | 1–10 saniye (gerçek zamanlı simülasyon için yapılandırılabilir) |
| Sensör gürültü modellemesi | Gauss gürültüsü (σ = 0.01–0.05) |
| Hata senaryoları | Sensör hataları, akış bozulmaları ve kontrol anomalileri |
| Değerlendirme modu | Kayan pencere tabanlı gerçek zamanlı simülasyon |
| Kod erişilebilirliği | Depo bağlantısı, tekrarlanabilirlik için talep üzerine / yayınlandığında sağlanacaktır |
Tablo 4: Simülasyon Ortamı. Deneysel düzendeki donanım konfigürasyonu, yazılım platformu, örnekleme aralığı, gürültü modellemesi ve hata senaryoları hakkında bilgiler.
| Veri Seti | Başlangıç Durumu Hatası | Final Durum Hatası | Sürüklenme Azaltımı (%) |
| Simüle Edilmiş Yağ & Gaz | 0.9 | 0.25 | 72.2 |
| SWAT | 1.1 | 0.4 | 63.6 |
| WADI | 1.3 | 0.55 | 57.7 |
Tablo 5: Dijital İkiz Senkronizasyonunun deneysel sonucu. Veri setleri genelindeki tahmin sapması azalmasını ve başlangıç ile son durum kestirim hatalarını gösteren senkronizasyon performansının nicel değerlendirmesi.
| Veri Seti | Hassasiyet | Geri çağırma | F1-Skoru |
| Simüle Petrol & Gaz | 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 |
Tablo 6: Anomali tespit performansı. Hassasiyet, geri çağırma ve F1-skoru dahil olmak üzere simüle edilmiş petrol ve gaz, SWAT ve WADI veri setleri üzerindeki anomali tespiti sonuçları.
| Veri Kümesi | Kümülatif Ödül ↑ | Kısıtlama İhlalleri ↓ | Maliyet Azalması (%) ↑ |
| Simüle Petrol & Gaz | 96.5 ± 0.8 | 0 | 18.7 |
| SWAT | 92.1 ± 1.2 | 1 | 14.3 |
| WADI | 89.7 ± 1.5 | 2 | 11.6 |
Tablo 7: Kapalı Döngü kontrol ve Optimizasyonun deneysel sonuçları. Kapalı döngü kontrol performansını değerlendirmek için kümülatif ödül, kısıtlama ihlalleri ve operasyonel maliyet azaltımı kullanılmıştır.
| Modül | Ortalama Gecikme (ms) | Açıklama |
| Nöral Operatör (FNO/DeepONet) | 18-25 ms | Operatör öğrenimi ve sistem durum tahmini |
| ST-GNN | 22–35 ms | Uzay-zamansal bağımlılık modellemesi |
| RL Politika Ağı | 5–10 ms | Kontrol kararının oluşturulması |
| Veri Ön İşleme | 8–12 ms | Girdi normalizasyonu ve özellik hazırlama |
| Toplam Kontrol Döngüsü Gecikmesi | 55–82 ms | Uçtan uca kontrol eylemi yürütme süresi |
Tablo 8: Gerçek zamanlı kontrol performans değerlendirmesi. Gecikme değerleri, standart bir donanım konfigürasyonu altında kontrol adımı başına ölçülen ortalama çıkarım sürelerini temsil etmektedir.
| Model | MSE ↓ | MAE ↓ | Bağıl Hata (%) ↓ | Çıkarım Süresi (ms) ↓ | p-değeri (Bazal değere karşı) |
| Sayısal Çözücü (Temel Hat) | 0.0125 ± 0.0012 | 0.089 ± 0.006 | 8.75 ± 0.54 | 120.5 ± 5.2 | — |
| Geleneksel ML Modeli | 0.0098 ± 0.0009 | 0.072 ± 0.005 | 6.42 ± 0.48 | 85.3 ± 4.7 | 0.021 |
| Mevcut Operatör Öğrenimi | 0.0076 ± 0.0007 | 0.061 ± 0.004 | 5.18 ± 0.41 | 42.8 ± 3.1 | 0.008 |
| DeepONet (Önerilen) | 0.0052 ± 0.0005 | 0.044 ± 0.003 | 3.67 ± 0.32 | 18.6 ± 2.4 | < 0.001 |
| FNO (Önerilen) | 0.0047 ± 0.0004 | 0.039 ± 0.002 | 3.21 ± 0.28 | 15.2 ± 2.1 | < 0.001 |
Tablo 9: Temel ve önerilen modellerin tahmin doğruluğu (MSE, MAE, Bağıl Hata) ve hesaplama verimliliği (Çıkarım Süresi) açısından karşılaştırması. Değerler, birden fazla çalışma üzerinden ortalama ± standart sapma olarak raporlanmıştır. Daha düşük değerler daha iyi performansı gösterir. p-değerleri, temel sayısal çözücüye kıyasla istatistiksel anlamlılığı belirtir.
| Yöntem | MAE ↓ | RMSE ↓ | F1-Skoru ↑ | AUC ↑ | Sürüklenme Azalması (%) ↑ | Kümülatif Ödül ↑ |
| Fizik tabanlı DT | 0.45 | 0.62 | 0.78 | 0.82 | 35.4 | 65.2 |
| LSTM | 0.32 | 0.48 | 0.85 | 0.88 | 48.7 | 74.5 |
| GNN | 0.28 | 0.41 | 0.88 | 0.91 | 55.3 | 81.3 |
| Önerilen Çatı | 0.2 | 0.3 | 0.94 | 0.97 | 72.2 | 95.5 |
Tablo 10: Kapsamlı performans karşılaştırması. Önerilen kontrol farkındalıklı dijital ikiz modelinin; izleme, anomali tespiti, senkronizasyon ve kontrol açısından diğer temel yaklaşımlarla kapsamlı performans karşılaştırması. Tahmin doğruluğu için MAE ve RMSE'nin düşük değerleri tercih edilir. Anomali tespiti ve senkronizasyon performansı için F1-skoru, AUC, kayma azalması ve ödülün yüksek değerleri tercih edilir.
| Bileşen | Parametre | Değer/Ayar |
| Sinirsel Operatörler (FNO/DeepONet) | Epoch Sayısı | 100–200 |
| Batch Boyutu | 32–64 |
| Öğrenme Oranı | 0.001 |
| ST-GNN | Katmanlar | 3 |
| Gizli Birimler | 64 |
| Zamansal Pencere | 10–20 zaman adımı |
| Sinirsel 4D-Var | Asimilasyon Penceresi | Sabit (kayan pencere) |
| Optimizasyon | Gradyan İnişi |
| Pekiştirmeli Öğrenme | İndirim Faktörü (γ) | 0.95 |
| Ödül Ağırlıkları (α,β,γ) | 0.6, 0.25, 0.15 |
| Genel Eğitim | Optimize Edici | Adam |
| Donanım | NVIDIA RTX GPU |
Tablo 11: Önerilen dijital ikiz çerçevesinin farklı bileşenleri için kullanılan eğitim konfigürasyonu ve uygulama parametreleri. Tablo, farklı bileşenler için önerilen sistemin tekrarlanabilirliğini sağlamak amacıyla kullanılan farklı hiperparametrelerin ve eğitim parametrelerinin genel bir görünümünü sunmaktadır.
| Bileşen | Eğitim Süresi (saat) | Çıkarım Gecikmesi (ms) | GPU Kullanımı (%) | Ölçeklenebilirlik |
| Sinirsel Operatörler | 6–8 | 25–40 | 80–90 | Yüksek |
| ST-GNN | 3–4 | 15–25 | 60–70 | Yüksek |
| Nöral 4D-Var | 2–3 | 20–30 | 65–75 | Orta düzeyde |
| RL Denetleyici | 1–2 | 10–20 | 50–60 | Yüksek |
| Genel Sistem | 12–16 | <100 ms | 70–85 | Ölçeklenebilir |
Tablo 12: Önerilen dijital ikiz modelinin farklı bileşenler için hesaplama performansı ve ölçeklenebilirlik analizi. Gecikme süresi gerçek zamanlı çıkarım performansını, GPU hesaplama performansını ve ölçeklenebilirlik ise önerilen modelin endüstriyel ölçekli uygulamalarda kullanılabilme yeteneğini temsil etmektedir.
| Yöntem | İzleme RMSE | Senkronizasyon Kayma Azalması (%) | AUROC | F1-Skoru | Kümülatif Ödül |
| Fizik tabanlı Dijital İkiz | 0.45 | 28.6 | 0.78 | 0.74 | 42.3 |
| LSTM tabanlı Model | 0.31 | 41.2 | 0.86 | 0.83 | 61.7 |
| GNN tabanlı Model | 0.26 | 53.4 | 0.91 | 0.88 | 74.5 |
| Önerilen Dijital İkiz | 0.2 | 72.2 | 0.94 | 0.92 | 95.5 |
Tablo 13: Önerilen çalışmanın mevcut yöntemlerle karşılaştırmalı değerlendirmesi. Önerilen dijital ikiz modelinin izleme, senkronizasyon, anomali tespiti ve kontrol metrikleri açısından performansı, temel ve öğrenme tabanlı yaklaşımlarla karşılaştırılmıştır.
| Model Konfigürasyonu | İzleme Hatası (RMSE) | Anomali Tespiti (F1) | Kontrol Performansı (Ödül) |
| Tam Model (Tüm Modüller) | 0.2 | 0.92 | 95.5 |
| – Veri Asimilasyonu Olmadan | 0.29 | 0.88 | 89.2 |
| – ST-GNN Olmadan | 0.27 | 0.85 | 87.1 |
| – Pekiştirmeli Öğrenme Olmadan | 0.21 | 0.9 | 78.4 |
| – Sürekli Öğrenme Olmadan | 0.23 | 0.91 | 92.3 |
Tablo 14: Deneysel sonuçlar ve Ablasyon Çalışması. Önerilen dijital ikiz modelinin her bir parçasının sistemin genel performansına nasıl katkıda bulunduğunu inceleyen ablasyon çalışması.
| Senaryo | Veri Kullanılabilirliği | Kök Ortalama Kare Hata (RMSE) | Doğruluk | Gözlem Boşluğu Yönetimi |
| Normal | 100% | Düşük | Yüksek | Stabil |
| Kısmi Kayıp | 50–70% | Orta derece | Yüksek | Robust interpolasyon |
| Ağır Kayıp | 20–30% | Artırılmış | Orta düzeyde | Kademeli bozulma |
| Sensör Arızası | Bölge bazlı %0 | Daha yüksek | Orta düzeyde | Uzamsal genelleme |
Tablo 15: Veri kıtlığı altında dayanıklılık analizi. Aşağıdaki tablo, önerilen nöral 4D-Var çerçevesinin farklı veri kullanılabilirliği durumlarındaki performansını göstermektedir. Tablo; sensör verilerinin kısmi kaybı, yüksek seyreklik ve toplam kaybı durumlarında hataların nasıl değiştiğini göstererek modelin dayanıklılığını ortaya koymaktadır.