Akademik araştırmalar, anlaşılır dil

Verianla | Akademik Araştırmalardan Türkçe Ekonomi ve Bilim İçerikleri

04 Kasım 2026, Çarşamba
VERİANLA Bağımsız bilim yayıncılığı
Verianla on Google Play
Menüyü aç veya kapat
...
Home / Uygulamalı Bilimler / Mühendislik / Rüzgâr Türbini Arızalarında Yorumlanabilir Erken Uyarı için Transfer Entropisi Tabanlı Dinamik Ağ Belirteçleri
Mühendislik

Rüzgâr Türbini Arızalarında Yorumlanabilir Erken Uyarı için Transfer Entropisi Tabanlı Dinamik Ağ Belirteçleri

Transfer Entropisi Tabanlı Yönlü Ağ (Transfer-Entropy-Based Directed Network, TEDN), rüzgâr türbini SCADA değişkenleri arasındaki yönlü öngörücü bilgi aktarımını her kayar zaman penceresinde transfer entropisiyle temsil eden bir ağ modelidir.

23/10/2026  Veri Anla  191 görüntüleme
Rüzgâr Türbini Arızalarında Yorumlanabilir Erken Uyarı için Transfer Entropisi Tabanlı Dinamik Ağ Belirteçleri

Transfer Entropisi Tabanlı Yönlü Ağ (Transfer-Entropy-Based Directed Network, TEDN), rüzgâr türbini SCADA değişkenleri arasındaki yönlü öngörücü bilgi aktarımını her kayar zaman penceresinde transfer entropisiyle temsil eden bir ağ modelidir. Çalışma, yalnız tek tek sıcaklık, güç veya akım değerlerinde belirgin bir sınır aşımı beklemek yerine, arızadan önce çok değişkenli sistemin etkileşim yapısında meydana gelen yeniden örgütlenmeyi izler. Ardışık TEDN ağları arasındaki yönlü bağ değişimleri; düğüm dalgalanması, birinci derece komşuluk içi değişim ve birinci/ikinci derece komşuluklar arası değişim ile birleştirilerek düğüm düzeyinde Transfer Entropy Dynamic Network Marker (TE-DNM) skoru hesaplanır. Düğüm skorlarının toplamı ise türbin için Comprehensive Warning Indicator (CWI), yani Kapsamlı Uyarı Göstergesi'ni oluşturur.

Yöntem EDP Renewables tarafından sağlanan, 2017 yılına ait 10 dakikalık SCADA kayıtları üzerinde değerlendirilmiştir. Aynı model üç türbin kullanılmıştır: T01, T06 ve T07. Normal durum kalibrasyonunda T01'in 8–12 Ocak 2017 tarihleri arasındaki 500 ardışık örneği kullanılmış; 100 örneklik, yani 1000 dakikalık pencereler ve bir örneklik kayma adımı uygulanmıştır. Normal durumda CWI ağırlıklı olarak 0–3 aralığında kaldığı için 3 değeri vaka çalışmalarında alarm eşiği olarak kullanılmıştır. Tek bir eşik aşımı yeterli kabul edilmemiş; herhangi beş ardışık pencerenin en az ikisinde CWI'nin eşiği aşması durumunda erken uyarı oluşturulmuştur.

T01 transformatör arızasında gürültüsüz temel durumda ilk geçerli öncül pencere 199, kayıtlı arıza penceresi 281 olarak belirlenmiş; bu fark 82 pencere veya yaklaşık 13,67 saate karşılık gelmiştir. Gaussian gürültüsü eklenen 15 tekrar deneyinde ilk uyarı pencereleri 183 ile 206 arasında kalmış ve kaynak, test edilen gürültü aralığında uyarı zamanının genel olarak kararlı kaldığını bildirmiştir. T06 dişli kutusu rulman arızasında TE-DNM, kayıtlı arızadan 17 pencere yani yaklaşık 2,83 saat önce; T07 hidrolik bileşen yağ kaçağında ise 141 pencere yani yaklaşık 23,5 saat önce geçerli uyarı üretmiştir.

Karşılaştırmada 3DV-DNM yöntemi T06 için 193. pencerede, TE-DNM ise 189. pencerede; T07 için 3DV-DNM 260. pencerede, TE-DNM ise 241. pencerede uyarı vermiştir. Bununla birlikte bu sonuçlar yalnız incelenen iki gerçek arıza vakası ve üç türbin için geçerlidir. Transfer entropisi bağlantıları fiziksel arıza nedenselliğini kanıtlamaz; yönlü öngörücü bilgi aktarımındaki istatistiksel değişimleri gösterir. Yüksek skorlu düğümler bakım personeline hangi ölçümlerin ve bileşenlerin önce incelenebileceğine ilişkin aday kanıt sunabilir, ancak gerçek arıza mekanizmasının belirlenmesi için bakım kayıtları ve doğrudan fiziksel inceleme gerekir.

Rüzgâr türbinlerinde erken uyarı problemi neden çok değişkenlidir?

Rüzgâr türbinleri değişken rüzgâr yükleri, çevre koşulları, mekanik yaşlanma, elektriksel yük değişimleri ve kontrol sistemi müdahaleleri altında çalışan karmaşık sistemlerdir. Erken aşamadaki bozulma, tek bir SCADA değişkenini hemen klasik alarm sınırının dışına çıkarmayabilir. Sıcaklıklar, akımlar veya dönme hızları tek tek normal kabul edilen aralıkta kalırken bu değişkenler arasındaki dinamik ilişkiler değişmeye başlayabilir.

Çalışmanın temel yaklaşımı bu nedenle “hangi sensör değeri anormal?” sorusundan önce “sistemin değişkenleri arasındaki etkileşim yapısı değişiyor mu?” sorusunu sorar. Bu yaklaşım kritik yavaşlama olarak bilinen dinamik sistem kavramıyla ilişkilidir. Sistem bir kritik geçişe yaklaştığında bozucu etkilerden toparlanma hızının azalması, dalgalanmaların kalıcı hale gelmesi ve değişkenler arasındaki etkileşim yapısının güçlenmesi veya yeniden örgütlenmesi beklenebilir.

Geleneksel dinamik ağ belirteçlerinin sınırı

Dinamik ağ belirteçleri, kritik geçiş öncesinde dalgalanması ve birlikte hareketi güçlenen değişkenlerden bir geçiş-duyarlı alt ağ oluşturabilir. Ancak birçok geleneksel DNM yaklaşımında ağ kenarları Pearson korelasyonu gibi simetrik ilişki ölçülerine dayanır. Böyle bir ağ A ile B'nin birlikte değiştiğini gösterebilir fakat A'nın geçmişinin B'nin geleceği hakkında ek öngörücü bilgi sağlayıp sağlamadığını ve ilişkinin yönünü ayırt edemez.

Bu çalışmanın temel yeniliği, simetrik korelasyon kenarının yerine transfer entropisine dayalı yönlü bir ağ kullanmasıdır. Böylece yalnız ilişkinin kuvveti değil, tahmine katkı sağlayan bilgi akışının yönü de modele dahil edilir.

Transfer Entropisi Tabanlı Yönlü Ağ (TEDN) Nedir?

TEDN, her kayar SCADA zaman penceresinde izlenen değişkenleri düğüm, değişkenler arasındaki anlamlı transfer entropisi değerlerini ise yönlü ve ağırlıklı kenar olarak kullanan bir ağdır; böylece sistemin yalnız birlikte değişimini değil, hangi değişkenin geçmiş bilgisinin başka bir değişkenin geleceğine ek öngörücü bilgi taşıdığını temsil eder.

Transfer entropisinin matematiksel anlamı

Kaynak, A değişkeninden B değişkenine yönlü transfer entropisini koşullu karşılıklı bilgi biçiminde tanımlar:

\[ TE_{A\rightarrow B}^{(w)} = I\left(X_t^{(A)};X_{t+\mathrm{lag}}^{(B)}\mid X_t^{(B)}\right) \]

Burada \(TE_{A\rightarrow B}^{(w)}\), \(w\) zaman penceresinde A'dan B'ye yönlü transfer entropisidir. \(X_t^{(A)}\), A'nın mevcut/geçmiş durumunu; \(X_t^{(B)}\), B'nin mevcut durumunu; \(X_{t+\mathrm{lag}}^{(B)}\) ise B'nin gecikme süresi kadar sonraki değerini temsil eder. Koşullu karşılıklı bilgi \(I(X;Y\mid Z)\), B'nin kendi geçmişi bilinirken A'nın geçmişinin B'nin geleceği hakkında ne kadar ek bilgi sağladığını ölçer.

Çalışma transfer entropisini parametrik doğrusal bir model yerine KSG-kNN tabanlı koşullu karşılıklı bilgi tahmincisiyle kestirir. Her örnek için kullanılan üçlü yapı şöyledir:

\[ (x_n,y_n,z_n) = \left( X_n^{(A)}, X_{n+\mathrm{lag}}^{(B)}, X_n^{(B)} \right) \]

Komşuluk uzaklığı ortak \((X,Y,Z)\) uzayında Chebyshev, yani maksimum norm ile hesaplanır:

\[ L_\infty(u_a,u_b) = \max_r |u_{a,r}-u_{b,r}| \]

Kaynakta verilen KSG tabanlı koşullu karşılıklı bilgi tahmincisi şu yapıya sahiptir:

\[ I_{KSG}(X;Y\mid Z) = \psi(k)+\psi(N) - \frac{1}{N} \sum_{a=1}^{N} \left[ \psi(n_{xz}(a)+1) + \psi(n_{yz}(a)+1) - \psi(n_z(a)+1) \right] \]

\(\psi\) digamma fonksiyonunu; \(k\) en yakın komşu derecesini; \(N\) pencere içindeki kullanılabilir örnek sayısını; \(n_{xz}\), \(n_{yz}\) ve \(n_z\) ise belirlenen komşuluk yarıçapı içinde ilgili alt uzaylarda bulunan örnek sayılarını gösterir.

SCADA verisi zaman pencerelerine nasıl ayrılıyor?

Ön işlenmiş veri, uzunluğu \(l\), kayma adımı \(h\) olan pencerelere ayrılır. Toplam \(T\) zaman noktası için kaynakta pencere sayısı:

\[ K=\left(\frac{T-l}{h}\right)+1 \]

biçiminde verilir. Her pencere, \(l\) zaman noktası ve \(m\) izleme değişkeninden oluşan \(l\times m\) boyutlu bir veri matrisi taşır. Her değişken çifti için transfer entropisi hesaplandıktan sonra pencereye özgü yönlü komşuluk matrisi elde edilir:

\[ TE^{(w)} = \left[ TE_{A\rightarrow B}^{(w)} \right]_{m\times m} \]

Her hesaplanan transfer entropisi kenar olarak tutuluyor mu?

Hayır. Çalışma önce her pencere içinde transfer entropisi bağlantılarının anlamlılığını permütasyon testiyle değerlendirir. Her çift için 50 permütasyon kullanılmış, çoklu test problemi için False Discovery Rate (FDR) kontrolü uygulanmış ve FDR eşiği 0,2 olarak ayarlanmıştır.

ParametreDeğer
Permütasyon sayısı50
FDR eşiği0,2
Ortalama anlamlı yönlü kenar oranı%30,914
Ortalama anlamlı kenar sayısı170,65
Komşuluk hesaplarına alınan bağlantılarEn yüksek TE değerine sahip %5'lik kesim

En güçlü %5'lik kenarın seçilmesi, düğüm komşuluklarını daha kuvvetli tahmini bağımlılıklar üzerine yoğunlaştırmak için kullanılan çalışma-parametresidir. Kaynak bunun başka türbinler veya veri setleri için optimal olduğunun gösterilmediğini açıkça belirtir.

TE-DNM ve CWI Erken Uyarıyı Nasıl Üretiyor?

TE-DNM, ardışık TEDN pencerelerindeki yönlü bağ değişimini düğüm ve komşuluk dalgalanmasıyla birleştirir; her düğümün TE-DNM skorları toplandığında CWI elde edilir ve CWI normal durumdan öğrenilen eşiği tekrarlı biçimde aştığında arıza öncülü uyarısı oluşturulur.

1. Ardışık pencerelerde yönlü bağ değişimi: STEC

A ile B arasındaki transfer entropisinin ardışık iki pencere arasındaki değişimi kaynakta Specific Transfer Entropy Coefficient (STEC) yapısıyla temsil edilir:

\[ STE_{A\rightarrow B}^{(w)} = TE_{A\rightarrow B}^{(w)} - TE_{A\rightarrow B}^{(w-1)} \]

İşaret, yönlü bağın kuvvetlenip zayıfladığını; mutlak değer ise değişimin büyüklüğünü gösterir:

\[ \left|STE_{A\rightarrow B}^{(w)}\right| = \left| TE_{A\rightarrow B}^{(w)} - TE_{A\rightarrow B}^{(w-1)} \right| \]

2. Düğüm ve birinci derece komşuluk dalgalanması: DCVI

Dynamic Directed Volatility Intensity (DCVI), merkez düğüm ile doğrudan yönlü bağlantılı birinci derece komşularının göreli değişkenliğini bir araya getirir. Kaynak, standart sapmanın ortalamaya göre normalize edildiği katsayıları ve sıfıra bölünmeyi önlemek için küçük pozitif \(\varepsilon\) sabitini kullanır:

\[ DCVI_A(w) = \frac{ \frac{\sigma_A(w)}{\mu_A(w)+\varepsilon} + \sum_{B\in N_1(A,w)} \frac{\sigma_B(w)}{\mu_B(w)+\varepsilon} }{ n_1(A,w)+1 } \]

\(N_1(A,w)\), A düğümünün birinci derece yönlü komşuluk kümesi; \(n_1(A,w)\) bu kümedeki düğüm sayısıdır. DCVI'nin görevi, yalnız ağ kenarındaki değişime değil, düğüm ve yakın komşularındaki ölçüm dalgalanmasına da duyarlılık kazandırmaktır.

3. Birinci derece komşuluk içindeki değişim: TEin

Birinci derece komşularla olan yönlü transfer entropisi değişimlerinin ortalaması:

\[ TEin_A(w) = \frac{1}{n_1(A,w)} \sum_{B\in N_1(A,w)} \left| TE_{A\rightarrow B}^{(w)} - TE_{A\rightarrow B}^{(w-1)} \right| \]

şeklinde hesaplanır. TEin, A düğümünün doğrudan yönlü çevresindeki bağ yeniden yapılanmasını ölçer.

4. Birinci ve ikinci derece komşuluklar arasındaki değişim: TEout

İkinci derece komşuluk \(N_2(A,w)\), A'nın birinci derece komşularına bağlı fakat A ve birinci derece komşuluk kümesi içinde olmayan düğümleri içerir. Kaynak, birinci ve ikinci derece komşuluklar arasındaki ortalama yönlü bağ değişimini:

\[ TEout_A(w) = \frac{1}{n_1(A,w)n_2(A,w)} \sum_{B\in N_1(A,w)} \sum_{C\in N_2(A,w)} \left| TE_{B\rightarrow C}^{(w)} - TE_{B\rightarrow C}^{(w-1)} \right| \]

ile tanımlar.

Bir düğüm için birinci veya ikinci derece komşuluk boşsa ilgili TE-DNM düğüm skoru sıfır kabul edilir. Kaynak bu sıfırın “normal çalışma” anlamına gelmediğini; yeterli ağ-komşuluk bilgisi bulunmadığını ifade ettiğini özellikle belirtir. Bir pencerede hiçbir geçerli düğüm yoksa geçerli CWI hesaplanmaz ve pencere uyarı değerlendirmesinden çıkarılır.

5. Düğüm düzeyindeki TE-DNM skoru

Üç yapı şu biçimde birleştirilir:

\[ TE\text{-}DNM_A(w) = \frac{ DCVI_A(w)\times TEin_A(w) }{ TEout_A(w)+\varepsilon } \]

Bu oran, merkez düğüm ve yakın komşulardaki dalgalanma ile yakın çevredeki yönlü bağ yeniden örgütlenmesini birlikte yükseltirken, daha geniş komşuluk çevresiyle bir karşılaştırma/normalizasyon sağlar.

6. Sistem düzeyindeki Kapsamlı Uyarı Göstergesi

Türbinin tüm geçerli düğüm skorları toplanarak:

\[ CWI(w) = \sum_{A=1}^{m} TE\text{-}DNM_A(w) \]

elde edilir. CWI, tek bir değişkenin genliğini değil, ağın farklı bölgelerinde ortaya çıkan yönlü dinamik yeniden yapılanmanın toplam işaretini temsil eder.

Yöntemin Verianla süreç görünümü

AşamaİşlemÜretilen bilgi
1SCADA verisini seç, hatalı/aykırı örnekleri işle ve değişkenleri standartlaştır.Karşılaştırılabilir çok değişkenli zaman serisi.
2100 örneklik kayar zaman pencereleri oluştur.Her biri 1000 dakikayı kapsayan yerel veri blokları.
3Her değişken çifti için yönlü transfer entropisini hesapla.Pencereye özgü TE matrisi.
4Permütasyon testi ve FDR kontrolü uygula; güçlü bağlantılardan TEDN komşuluklarını oluştur.Anlamlı yönlü bilgi aktarımı ağı.
5Ardışık pencereler arasında STEC, DCVI, TEin ve TEout hesapla.Yerel dalgalanma ve yönlü bağ yeniden örgütlenmesi.
6Her geçerli düğüm için TE-DNM skoru hesapla.Değişken düzeyinde yorumlanabilir katkı.
7Düğüm skorlarını toplayarak CWI oluştur.Türbin düzeyinde uyarı göstergesi.
8CWI'nin normal durumdan kalibre edilen eşiği tekrarlı biçimde geçip geçmediğini kontrol et.Normal çalışma veya arıza öncülü uyarısı.

Transfer Entropisi Fiziksel Nedenselliği Kanıtlar mı?

Hayır. Bu çalışmada transfer entropisi “nedensel” sözcüğünü bilgi kuramsal, yönlü öngörücü bağımlılık anlamında kullanır; A'dan B'ye yüksek veya değişen bir TE bağlantısı, A'nın fiziksel olarak B'deki arızaya neden olduğunu kanıtlamaz ve gerçek arıza mekanizmasının doğrulanması için çalışma kayıtları, bakım bulguları ve doğrudan bileşen incelemesi gerekir.

Bu ayrım yöntemin yorumlanabilirlik iddiası açısından kritiktir. Yüksek düğüm skoru, bakım personeline hangi değişkenlerin veya ilişkilerin dikkat çekici biçimde yeniden örgütlendiğini gösterir. Ancak “yüksek skor alan değişken arızanın sebebidir” sonucu çıkarılamaz.

Çalışmanın Yöntemi ve Bulguları

Kullanılan SCADA veri seti

Çalışma, EDP Renewables tarafından yayımlanan Wind farm 1 - Wind turbine SCADA signals 2017 veri setini kullanır. Ölçümler 10 dakikalık aralıklarla kaydedilmiştir. Aynı model T01, T06 ve T07 türbinleri seçilmiştir. Türbinlerin nominal gücü 2 MW, nominal rüzgâr hızı 12 m/s, devreye girme hızı 4 m/s ve devreden çıkma hızı 25 m/s olarak raporlanmıştır.

İzlenen 24 sürekli SCADA değişkeni

No.SCADA değişkeniNo.SCADA değişkeni
1Ortalama dönme hızı13Ortalama jeneratör slip-ring sıcaklığı
2Ortalama rulman sıcaklığı14Ortalama rotor tarafı konvertör sıcaklığı 1
3Ortalama jeneratör faz-1 sıcaklığı15Ortalama rotor tarafı konvertör sıcaklığı 2
4Ortalama jeneratör faz-2 sıcaklığı16Ortalama rotor tarafı konvertör sıcaklığı 3
5Ortalama jeneratör faz-3 sıcaklığı17Şebeke bağlantı noktasındaki ortalama aktif güç
6Ortalama hidrolik yağ sıcaklığı18Şebeke bağlantı noktasındaki ortalama faz-1 gerilimi
7Ortalama dişli kutusu yağ sıcaklığı19Şebeke bağlantı noktasındaki ortalama faz-2 gerilimi
8Ortalama dişli kutusu rulman sıcaklığı20Şebeke bağlantı noktasındaki ortalama faz-3 gerilimi
9Ortalama yüksek gerilim transformatörü faz-A sıcaklığı21Şebeke bağlantı noktasındaki ortalama faz-1 akımı
10Ortalama yüksek gerilim transformatörü faz-B sıcaklığı22Şebeke bağlantı noktasındaki ortalama faz-2 akımı
11Ortalama yüksek gerilim transformatörü faz-C sıcaklığı23Şebeke bağlantı noktasındaki ortalama faz-3 akımı
12Gerçek konvertör sıcaklığı24Ortalama jeneratör rulman sıcaklığı

Veri ön işleme ve normal durum kalibrasyonu

Ön işleme aşamasında Local Outlier Factor (LOF) yöntemi aykırı örnekleri belirlemek için kullanılmış, çıkarılan değerler KNN imputasyonu ile tamamlanmış ve ardından değişkenler standardize edilmiştir.

Normal durum referansı için T01 türbininin 8–12 Ocak 2017 arasındaki 500 ardışık örneği seçilmiştir. Her pencere 100 örnek içerdiğinden bir pencere 1000 dakikalık çalışma süresini kapsar. Kayma adımı \(h=1\)'dir. 500 örnekten 401 ağ penceresi oluşturulmuş; ağ değişim hesabı ardışık ağ çiftlerini gerektirdiğinden 400 geçerli CWI değeri elde edilmiştir.

ParametreDeğer
SCADA örnekleme aralığı10 dakika
Normal durum örnek sayısı500
Pencere uzunluğu \(l\)100 örnek
Bir pencerenin fiziksel süresi1000 dakika
Kayma adımı \(h\)1 örnek
Oluşturulan ağ penceresi401
Geçerli ardışık CWI karşılaştırması400
Normal referans aralığıAğırlıklı olarak [0, 3]
Uyarı eşiğiCWI = 3
Onay kuralıHerhangi 5 ardışık pencerenin en az 2'sinde eşik aşımı

Kaynak, 3 eşiğinin farklı türbin modellerine veya farklı çalışma ortamlarına doğrudan taşınmaması gerektiğini özellikle belirtir. Yeni bir türbin veya işletme koşulunda eşik, o birimin kendi normal durum dağılımından yüzdelik veya ekstrem-değer ölçütü gibi yöntemlerle yeniden kalibre edilmelidir. Beş pencere içinde iki aşım kuralının da evrensel optimum olduğu test edilmemiştir.

T01 transformatör arızası: temel erken uyarı sonucu

T01 transformatör arızası 11 Ağustos 2017 saat 13:10'da kaydedilmiştir. Analiz edilen dizide kayıtlı arıza 281. pencereye karşılık gelmektedir. Gürültü eklenmeyen durumda ilk geçerli öncül eşik aşımı 199. pencerede belirlenmiştir. Böylece arızayla ilk öncül arasında 82 pencere, 10 dakikalık örnekleme aralığında yaklaşık 13,67 saat bulunmuştur.

Kaynak burada önemli bir zamanlama ayrımı yapar: retrospektif “öncül başlangıç” zamanı, iki-of-beş doğrulama kriterinin daha sonra tamamlanmasına yol açan en erken eşik aşımıdır. Dolayısıyla 13,67 saatlik değer gerçek zamanlı bir sistemin alarmı operatöre kesin olarak verdiği andan hesaplanan çevrim içi alarm süresi olarak yorumlanmamalıdır.

Gaussian gürültü dayanıklılığı

T01 vakasında orijinal veriye \(\sigma=0{,}1\), \(0{,}2\) ve \(0{,}5\) olmak üzere üç Gaussian gürültü seviyesi eklenmiş ve her seviye için beş bağımsız gerçekleştirim oluşturulmuştur. Böylece gürültüsüz temel durum dışında toplam 15 gürültülü tekrar değerlendirilmiştir.

Gürültü seviyesi \(\sigma\)Tekrarİlk uyarı penceresi ortalamasıStandart sapmaBeş gerçekleşimdeki pencereler
0,15197,45,083191, 199, 197, 194, 206
0,25197,60,800197, 199, 197, 197, 198
0,55193,05,060196, 196, 183, 196, 194

On beş gürültülü çalışmada ilk uyarı penceresi 183–206 aralığında kalmıştır. Kaynak bu sonucu test edilen gürültü aralığı ve T01 transformatör vakası için gösterge kararlılığına destek olarak yorumlar; bütün sensör gürültüsü türleri, bütün çalışma koşulları veya bütün türbinler için genel dayanıklılık kanıtı olarak sunmaz.

Ablasyon: DCVI çıkarıldığında ne oluyor?

DCVI'nin katkısını test etmek için TEDN korunmuş fakat \(DCVI=1\) sabitlenerek düğüm dalgalanması modülasyonu kaldırılmıştır. Bu konfigürasyonun kendi normal durum ölçeğinde CWI değerleri [0,1] aralığında kalmış ve 1 değeri alarm eşiği olarak kullanılmıştır. T01 arıza verisinde CWI hiçbir geçerli pencerede bu eşiğe ulaşmamış ve doğrulanmış alarm üretilememiştir.

Bu sonuç, incelenen T01 vakasında düğüm ve birinci derece komşuluk dalgalanmasını taşıyan DCVI bileşeninin erken uyarı duyarlılığına katkıda bulunduğunu göstermektedir. Ancak farklı ablasyonların CWI ölçekleri ve alarm eşikleri aynı olmadığından yalnız mutlak CWI büyüklüklerine bakılarak yöntem üstünlüğü çıkarılamaz.

Ablasyon: Transfer entropisi yerine Pearson korelasyonu kullanıldığında ne oluyor?

İkinci ablasyonda yönlü transfer entropisi kenarları Pearson korelasyon katsayısı (PCC) ile değiştirilmiştir. PCC tabanlı normal durum için [0, 0,25] referans aralığı ve 0,25 eşik değeri belirlenmiştir.

T01 arıza verisinde PCC tabanlı CWI yalnız 49. ve 56. pencerelerde eşiğe ulaşmıştır. Bu iki pencere aynı beş pencerelik aralıkta bulunmadığı için iki-of-beş uyarı kuralı karşılanmamış ve doğrulanmış alarm oluşmamıştır. Kaynak bu vakada TE tabanlı yönlü ağın başlangıç arıza işaretlerine PCC tabanlı yönsüz ilişki ağından daha duyarlı olduğunu değerlendirmektedir.

T06 dişli kutusu rulman arızası

T06 için kayıtlı olay 17 Ekim 2017 saat 08:38'de gerçekleşen dişli kutusu rulman arızasıdır. CWI eşik aşımı 189, 190, 198 ve 201. pencerelerde meydana gelmiştir. Kayıtlı arıza 206. pencerededir.

ÖlçütTE-DNM3DV-DNMOperasyon kaydı
İlk geçerli uyarı / arıza189. pencere193. pencere206. pencere
TE-DNM öncül süresi17 pencere ≈ 2,83 saatKaynak TE-DNM için hesaplamıştır.

Karşılaştırmada 3DV-DNM de arızadan önce uyarı vermiş, ancak TE-DNM'nin ilk geçerli uyarısı dört pencere daha erken gerçekleşmiştir. Bununla birlikte iki yöntemin alarm doğrulama kuralları aynı değildir: TE-DNM iki-of-beş, 3DV-DNM ise kendi yöntem tanımına göre üç-of-beş eşik aşımı kullanır. Bu nedenle sonuç “her koşulda mutlak üstünlük” değil, kaynakta uygulanan yöntem-özel kurallar altındaki karşılaştırmadır.

T07 hidrolik bileşen yağ kaçağı

T07'de 17 Haziran 2017 saat 11:35'te hidrolik bileşende yağ kaçağı kaydedilmiştir. İlk CWI eşik aşımı 241. pencerede ortaya çıkmış; 245, 259, 260, 307, 308, 329 ve 330. pencerelerde ek aşım olayları gözlenmiştir. Kayıtlı arıza 382. pencerededir.

ÖlçütTE-DNM3DV-DNMOperasyon kaydı
İlk geçerli uyarı / arıza241. pencere260. pencere382. pencere
TE-DNM öncül süresi141 pencere ≈ 23,5 saatKaynak TE-DNM için hesaplamıştır.

İki gerçek arıza vakasının ikisinde de TE-DNM, kaynakta kullanılan 3DV-DNM karşılaştırmasından daha erken ilk uyarı üretmiştir. Ancak karşılaştırma yalnız tek benchmark yöntemi ve iki gerçek arıza türüyle sınırlıdır.

Üç ana vaka tek tabloda

TürbinOlayİlk TE-DNM uyarı penceresiKayıtlı arıza penceresiÖncül farkı
T01Transformatör arızası19928182 pencere ≈ 13,67 saat
T06Dişli kutusu rulman arızası18920617 pencere ≈ 2,83 saat
T07Hidrolik bileşen yağ kaçağı241382141 pencere ≈ 23,5 saat

Önemli yorum sınırı: T01, T06 ve T07 için verilen bu süreler retrospektif vaka analizlerinde elde edilmiştir. Kaynakta özellikle T01 için erken “öncül başlangıç” zamanının, çevrim içi sistemin alarm kriterini fiilen tamamladığı andan daha erken olabileceği belirtilir. Bu nedenle değerler gerçek saha alarm sisteminde garanti edilmiş operasyonel reaksiyon süreleri olarak kullanılmamalıdır.

Düğüm skorları yöntemi nasıl yorumlanabilir hale getiriyor?

Sistem düzeyindeki CWI yalnız “uyarı var/yok” bilgisi üretmez. Düğüm TE-DNM skorları hangi SCADA değişkenlerinin ağ yeniden örgütlenmesine daha fazla katkı verdiğini de gösterir. Kaynak bu sıralamayı bakım incelemesini önceliklendirmek için aday kanıt olarak değerlendirir.

VakaÖne çıkan değişkenlerKaynakta verilen yorum
T01 transformatör9 en yüksek; ayrıca 5, 6, 12, 14 ve 15Transformatör sıcaklığı doğrudan bileşen termal durumuyla; jeneratör, akım/konvertör ve yardımcı sıcaklıklar ise yük ve işletme durumu değişimleriyle birlikte incelenmelidir.
T06 dişli kutusu21, 22 ve 17 belirgin; ayrıca 2, 14 ve 6Akım ve aktif güç düğümleri tahrik sistemi yük değişimleriyle ilişkili olabilir; rulman ve sıcaklık ölçümleri destekleyici kanıttır.
T07 hidrolik yağ kaçağı22 en yüksek; 7 ve 14 de yüksekElektriksel yük, dişli kutusu işletimi ve sıcaklık değişimleri öncül durum değişikliklerini gösterebilir; ancak yağ kaçağının yerini doğrudan tanımlamaz.

T07 için kaynak özellikle basınç, yağ seviyesi veya doğrudan kaçak akışı gibi ölçümlerin daha doğrudan hidrolik arıza yorumu sağlayacağını belirtir. Kullanılan açık SCADA veri setinde ayrıntılı bakım raporları ve bileşen inceleme bulguları mevcut olmadığı için düğüm skorlarından fiziksel arıza yayılım yolu çıkarılmamıştır.

Uyarı sonrasında önerilen mühendislik yaklaşımı

Çalışma CWI artışının otomatik olarak “parçayı değiştir” komutu olması gerektiğini savunmaz. Kaynakta önerilen yaklaşım katmanlıdır:

  1. Önce sensör geçerliliği ve veri kalitesi kontrol edilir.
  2. En yüksek skorlu düğümler ve yönlü bağlantılar incelenir.
  3. Mevcut örüntü önceki normal çalışma dönemiyle karşılaştırılır.
  4. Güç, akım, dönme hızı ve ilgili sıcaklık kayıtları birlikte değerlendirilir.
  5. İşaret kalıcıysa ilgili bileşenlerin hedefli fiziksel muayenesi planlanır.
  6. Bakım, yük azaltma veya derating kararı mevcut operasyon kanıtlarıyla birlikte verilir.
  7. Yanlış alarm sayıları ve sonraki fiziksel inceleme sonuçları kaydedilerek yöntemin gerçek mühendislik değeri izlenir.

Çalışmanın desteklediği sonuçlar

  • Yönlü transfer entropisi, test edilen SCADA verilerinde arıza öncesi çok değişkenli ilişki yeniden örgütlenmesini temsil etmek için kullanılabilmiştir.
  • Normal T01 referans döneminde CWI ağırlıklı olarak [0,3] aralığında kalmıştır.
  • T01 transformatör vakasında test edilen Gaussian gürültü seviyeleri altında ilk uyarı penceresi genel olarak kararlı kalmıştır.
  • DCVI'nin kaldırıldığı ablasyonda T01 için doğrulanmış uyarı üretilememiştir.
  • PCC tabanlı yönsüz ağ ablasyonunda T01 için iki-of-beş alarm kriteri karşılanmamıştır.
  • T06 ve T07 gerçek arıza vakalarında TE-DNM, kayıtlı olaylardan sırasıyla yaklaşık 2,83 ve 23,5 saat önce öncül üretmiştir.
  • T06 ve T07'de TE-DNM, kaynakta kullanılan 3DV-DNM benchmark'ından daha erken ilk uyarı vermiştir.
  • Düğüm skorları hangi ölçümlerin uyarıya daha fazla katkı verdiğini görünür hale getirerek bakım incelemesini yönlendirebilecek aday kanıt sağlamaktadır.

Çalışmanın desteklemediği sonuçlar

  • TE bağlantıları fiziksel nedensellik veya arıza yayılım yolunu kanıtlamaz.
  • CWI=3 eşiğinin bütün rüzgâr türbinlerine doğrudan uygulanabileceği gösterilmemiştir.
  • 100 örneklik pencere, %5 bağlantı seçimi, FDR=0,2 veya iki-of-beş doğrulama kuralının bütün veri setleri için optimal olduğu gösterilmemiştir.
  • Üç türbin ve iki gerçek arıza tipi üzerinden bütün arıza türlerinin kapsandığı sonucu çıkarılamaz.
  • Test edilen Gaussian gürültü seviyeleri bütün sensör hataları veya gürültü dağılımları için genel dayanıklılık kanıtı değildir.
  • Yüksek düğüm skoru tek başına ilgili fiziksel bileşenin arızanın nedeni olduğunu göstermez.
  • Retrospektif öncül zamanı çevrim içi sahada operatöre verilecek onaylanmış alarm süresiyle otomatik olarak aynı değildir.
  • Yanlış alarm oranı uzun normal çalışma dönemlerinde henüz kapsamlı biçimde nicelleştirilmemiştir.
  • 3DV-DNM dışındaki erken uyarı yöntemlerine göre genel üstünlük gösterilmemiştir.

Temel sınırlılıklar ve gelecek doğrulama gereksinimleri

Kaynak, mevcut doğrulamanın bir normal referans dönemi, üç türbin ve iki gerçek arıza tipiyle sınırlı olduğunu belirtir. Yük değişimleri, start-stop geçişleri ve kontrol sistemi eylemleri hem ölçüm değişkenliğini hem de yönlü öngörücü bağımlılıkları etkileyebilir. Bu nedenle CWI artışı yalnız gösterge değerine bakılarak arıza öncülü kabul edilmemeli; güç, akım, dönme hızı ve işletme kayıtlarıyla birlikte değerlendirilmelidir.

Gelecek çalışmalar için uzun normal işletme dönemlerinde yanlış alarm oranlarının ölçülmesi, daha fazla türbin ve arıza tipinde doğrulama, TE ve komşuluk parametrelerine duyarlılık analizi, güç düzeyi/yük değişimi/işletme durumuna göre katmanlı değerlendirme ve yüksek katkılı kenarların bakım kayıtlarıyla karşılaştırılması önerilmektedir.

Kaynak ve Yöntem Notu

Özgün çalışma: Transfer-Entropy-Based Dynamic Network Markers for Interpretable Early Warning of Wind Turbine Failures.

Yazarlar: Tongding Yu, Ruiming Fang, Nan Ying ve Yihan Huang.

Sorumlu yazar: Ruiming Fang.

Kurum: College of Information Science and Engineering, Huaqiao University, Xiamen 361021, China.

Kaynak türü: Araştırma preprinti / working-paper niteliğinde teknik çalışma.

Hakemlik durumu: PDF her sayfada çalışmanın henüz hakem değerlendirmesinden geçmediğini belirtmektedir.

Platform: SSRN.

SSRN kimliği: abstract=7625808.

DOI: Yüklenen preprint içinde çalışma için doğrulanabilir bir DOI verilmemiştir; bibliyografik boşluk kaynak dışı tahminle doldurulmamıştır.

Veri kaynağı: EDP Open Data, Wind farm 1 - Wind turbine SCADA signals 2017. Kaynak çalışma veri setine 7 Temmuz 2025 erişim tarihi vermektedir.

Temel yöntem: Kayar pencere içinde KSG-kNN ile koşullu karşılıklı bilgi tabanlı transfer entropisi kestirimi; permütasyon testi ve FDR; yönlü TEDN oluşturulması; STEC, DCVI, TEin ve TEout ile düğüm düzeyinde TE-DNM skoru; düğüm skorlarının CWI içinde toplanması; normal durum verisinden eşik kalibrasyonu.

Yöntem sınıfı: Gözetimsiz normal durum kalibrasyonlu, çok değişkenli SCADA tabanlı arıza öncülü tespit yaklaşımı. Çalışma büyük etiketli arıza kütüphanesi üzerinden denetimli sınıflandırıcı eğitmemektedir.

Finansman: National Natural Science Foundation of China, hibe 52477048; Xiamen Key Innovation Platforms Major Science and Technology Projects, hibe 3502Z20261010; Zhangzhou Municipal Science and Technology Program, Open Competition Project, hibe ZZ2025JB02.

Çıkar çatışması: Yazarlar çalışmayı etkileyebilecek bilinen rekabet eden mali çıkar veya kişisel ilişki olmadığını beyan etmektedir.

Veri erişilebilirliği: Kaynak SCADA veri seti EDP Open Data üzerinde açık olarak belirtilmiştir. İşlenmiş veri alt kümeleri ve analiz kodu ise corresponding author'dan makul talep üzerine sağlanabileceği şeklinde raporlanmıştır.

Yazar katkıları: Tongding Yu: kavramsallaştırma, veri kürasyonu, biçimsel analiz, araştırma, metodoloji, yazılım, doğrulama, görselleştirme ve ilk taslak. Ruiming Fang: kavramsallaştırma, finansman edinimi, metodoloji, proje yönetimi, kaynaklar, danışmanlık ve gözden geçirme/düzenleme. Nan Ying: doğrulama. Yihan Huang: araştırma.


Paylaş:

Yorumlar incelendikten sonra yayımlanır. Gönderdiğiniz yorum onay sürecine alınır ve uygun bulunduğunda görünür hâle gelir.

Bir yorum bırakın

E-posta adresiniz yayınlanmayacaktır. Gerekli alanlar * ile işaretlenmiştir