Elektryczna tomografia pojemnościowa

Sekwencja obrazów ECT – obiekt wykonany ze szkła akrylowego, obracający się wewnątrz sondy pomiarowej
System pomiarowy – Tomograf pojemnościowy[1] z podłączoną 16-elektrodową sondą

Elektryczna tomografia pojemnościowa (ang. electrical capacitance tomography – ECT) – technika obrazowania pozwalająca na wizualizację przestrzennego rozkładu przenikalności elektrycznej w badanej objętości. Obraz przestrzennego rozkładu przenikalności elektrycznej jest rekonstruowany na podstawie pomiarów pojemności wzajemnych elektrod powierzchniowych otaczających badaną objętość[2].

Technika ta na tle innych tomografii wyróżnia się bardzo dużą częstością pomiarów, dzięki czemu umożliwia obrazowanie zjawisk dynamicznych[3]. Duża liczba obrazów na sekundę, skalowalność sondy pojemnościowej, jak również możliwość wbudowania sondy w instalację przemysłową sprawia, że technika ta znajduje zastosowanie w monitorowaniu procesów dynamicznych i jest zaliczana do technik tomografii procesowej. Ze względu na bardzo małą rozdzielczość przestrzenną, ECT nie znalazła jak dotąd zastosowania w diagnostyce medycznej. Potencjalnie ECT może mieć podobne medyczne zastosowania co elektryczna tomografia impedancyjna, takie jak monitorowanie pracy płuc[4] lub wykrywanie niedokrwienia lub krwotoku w mózgu[5].

Chociaż w XX wieku metody pomiaru pojemności były w powszechnym użyciu, pomysł wykorzystania ich jako metody obrazowania przypisuje się zespołowi badawczemu pod przewodnictwem Maurica Becka z Instytutu Naukowo-Technicznego Uniwersytetu w Manchesterze (dziś Uniwersytet Manchesterski). Grupa ta, którą współtworzył również polski naukowiec Andrzej Pląskowski, w latach 80 opublikowała pierwszy artykuł naukowy na temat elektrycznej tomografii pojemnościowej[2].

Pomiar małych pojemności

Pojemność wzajemna elektrod mierzona w elektrycznej tomografii pojemnościowej jest na poziomie ułamków femtofaradów. Pomiar tak małych wartości jest dużym wyzwaniem dla układów elektronicznych. Powinny one być nieczułe na pojemności pasożytnicze i mieć szum własny zredukowany do minimum. Wskazane jest również ekranowanie układów, aby ograniczyć wpływ zewnętrznego pola elektrycznego[6].

Metoda z pobudzeniem sinusoidalnym

W metodzie z pobudzeniem sinusoidalnym, nazywanej również metodą AC, na mierzoną pojemność, podłączoną szeregowo, podawane jest napięcie przemienne. Pomiar odbywa się przy pomocy wzmacniacza prądu zmiennego. Wzmocniony sygnał jest następnie mieszany z sygnałem o częstotliwości napięcia pobudzającego (detekcja synchroniczna). Sygnały o tej samej częstotliwości i będące w zgodnej fazie dają po zmieszaniu napięcie stałe. Uzyskany sygnał wygładzany jest filtrem dolnoprzepustowym, w efekcie uzyskuje się stałe napięcie o amplitudzie wprost proporcjonalnej do mierzonej pojemności. Wolne zmiany napięcia są traktowane jako zmiany pojemności. Detekcja synchroniczna charakteryzuje się wysokim stosunkiem sygnału do szumu. Wadą tej metody jest względnie mała prędkość pomiaru, zależna od częstotliwości sygnału pobudzającego oraz częstotliwości granicznej filtru dolnoprzepustowego[7].

Metoda z pobudzeniem impulsowym

Metoda z pobudzeniem impulsowym polega na cyklicznym ładowaniu i rozładowaniu mierzonej pojemności. Impulsy prądu ładowania i rozładowania są całkowane w dwóch układach całkujących przełączanych za pomocą kluczy analogowych synchronicznie z prostokątnym sygnałem wymuszającym. Napięcie na integratorach jest proporcjonalne do mierzonej pojemności. Obie wartości są sumowane za pomocą wzmacniacza pomiarowego. Metoda ta jest bardzo szybka i minimalizuje wpływ zakłóceń zewnętrznych. Dla uzyskania mierzalnego sygnału stosuje się uśrednianie[8].

Szczególnym przykładem pomiaru impulsowego jest metoda polegająca na pobudzaniu elektrod napięciem rzędu kilkuset woltów. Sygnał pobudzenia poprzez pojemność mierzoną trafia na wzmacniacz pomiarowy. Mierzone napięcie jest proporcjonalne do sygnału pobudzającego i badanej pojemności. Zastosowanie stosunkowo wysokiego dla elektroniki pomiarowej napięcia wymuszenia zwiększa mierzony sygnał i poprawia SNR[6].

Sekwencja pomiarowa

Schemat przykładowej sondy pomiarowej
Pierścień elektrod zamocowanych do rury PCW

W elektrycznej tomografii pojemnościowej typowo stosuje się wielokanałowe systemy pomiarowe, w których do każdej elektrody pomiarowej przyłączony jest osobny układ odbiorczy. Jedna z elektrod pełni funkcję elektrody nadawczej, a na wszystkich pozostałych elektrodach jednocześnie odbywa się pomiar. Tomograficzna sonda pod względem elektrycznym jest układem kondensatorów, w którym rolę okładki dodatniej pełni elektroda nadawcza, a pozostałe elektrody wspólnie są okładkami ujemnymi. Trudnością tak dokonywanych pomiarów jest duży zakres mierzonych wartości. Pojemności wzajemne elektrod przyległych do elektrody nadawczej sięgają pikofaradów, podczas gdy dla elektrod przeciwległych do elektrody nadawczej nie przekraczają ułamków femtofaradów. Z tego powodu przed każdym pomiarem we wszystkich układach odbiorczych w zależności od ułożenia względem elektrody nadawczej ustawiane są odpowiednie poziomy wzmocnień[9]. Każda z elektrod po kolei pełni funkcję elektrody nadawczej, aż do uzyskania wszystkich dostępnych kombinacji. Liczba pomiarów jest równa dwuelementowej kombinacji z L-elementowego zbioru gdzie to liczba elektrod pomiarowych. Jeżeli wykorzystywane są również powtórzenia pomiarów dla par elektrod ze zamianą elektrody wymuszającej, liczbę uzyskiwanych pomiarów mnoży przez liczbę permutacji, czyli zwiększa się dwa razy[10].

Przy założeniu, że w obrazowanej przestrzeni zachodzi przepływ laminarny o znanej prędkości, możliwe jest składanie wielu zrekonstruowanych obrazów dwuwymiarowych z następujących po sobie pomiarów w jeden obraz objętościowy. Takie działanie określa się mianem tomografii dwuipółwymiarowej[11]. W ECT możliwe jest również uzyskiwanie obrazów trójwymiarowych. W tym celu stosuje się sondy wyposażone w wiele równoległych pierścieni elektrod, między którymi odbywa się pomiar[12].

Elektryczna tomografia pojemnościowa 2,5D
Obiekt przepływający przez sondę pomiarową i uzyskane obrazy 2D
Złożony obraz objętościowy

Rekonstrukcja obrazu

Wartości mierzonych pojemności wzajemnych są silnie związane ze stosunkiem powierzchni elektrod do średnicy sondy pomiarowej. Zmniejszanie szerokości elektrod nie może być równoważone zwiększaniem ich wysokości, gdyż pogarszałoby to zdolność rozdzielczą tomografu w osi Z cylindra pomiarowego. Ze względu na ograniczone możliwości układów pomiarowych (najmniejsze mierzalne pojemności są na poziomie ułamków femtofaradów), niemożliwe jest dowolne zwiększanie liczby elektrod. W praktyce nie stosuje się więc więcej niż 16 elektrod pomiarowych w jednym pierścieniu[13]. Tak mocno ograniczona liczba elektrod skutkuje małą liczbą uzyskiwanych pomiarów.

Rekonstrukcja obrazów z projekcji, gdy mamy stosunkowo mały zbiór pomiarów jest bardzo trudna, problem odwrotny jest niedookreślony (nieoznaczony układ równań). W efekcie rozdzielczość przestrzenna tej techniki jest słaba w porównaniu do innych metod tomograficznych[10].

Problem prosty

Problem prosty polega na wyznaczeniu wyniku obserwacji w danym modelu dla założonych danych wejściowych. W ECT, znając model dla przyjętego rozkładu przestrzennego przenikalności elektrycznej, można wyznaczyć wartości pomiarów pojemności.

Pomiar można opisać równaniem: gdzie to wektor pomiarów pojemności, a oznacza nielinową, wektorową funkcję rozkładu przenikalności opisywanego wektorem Korzystając z rozwinięcia w szereg Taylora, nielinową funkcję można przybliżyć modelem liniowym gdzie

jest macierzą pochodnych cząstkowych. jest nazywana macierzą wrażliwości, ponieważ jej współczynniki określają wpływ zmiany wartości przenikalności elektrycznej w małym elemencie objętości sondy pomiarowej na pojemności wzajemne jej elektrod. Macierzy wrażliwości można wyznaczyć metodą zaburzeń eksperymentalnie lub w symulacji numerycznej.


Dla założonego rozkładu przenikalności elektrycznej wyznacza się rozkład potencjału elektrycznego w sondzie pomiarowej. Korzysta się w tym celu z jednej z metod numerycznego rozwiązywania problemów brzegowych, takich jak metoda różnic skończonych, metoda elementów skończonych lub metoda objętości skończonych[14]. Korzystając z prawa Gaussa, wyznacza się ładunek na jednej z elektrod pomiarowych. Znając różnicę potencjału pomiędzy elektrodami i można wyznaczyć pojemność wzajemną -tej pary elektrod Dla rozkładu przenikalności, w którym wprowadza się małe zaburzenie w małym -tym elemencie objętości, powtarza się obliczenia numeryczne dla wyznaczenia nowej wartości pojemności i stąd wyznacza się zmianę pojemności Powtarzając symulacje numeryczne dla kolejnych par elektrod i kolejnych elementów objętości uzyskuje się pochodne cząstkowe

W praktyce ze względu na konieczność powtórzenia symulacji numerycznych, macierzy wrażliwości nie wyznacza się tą metodą, lecz korzysta się ze wzoru Geselowitza[15]:

gdzie:

– pojemność mierzona między elektrodami i
– wartość przenikalności elektrycznej w punkcie o indeksie
– różnica potencjału między elektrodami i
– wektor opisujący rozkład potencjału w punktach siatki dyskretnej, gdy elektroda jest źródłem dodatniego wymuszenia,
– wektor opisujący rozkład potencjału, gdy elektroda jest źródłem dodatniego wymuszenia.

Korzystając z powyższego wzoru, wyprowadzonego z teorii impedancji wzajemnej, rozkład potencjału jest wyznaczany -razy, gdzie jest liczbą elektrod sondy.

Problem odwrotny

Problem odwrotny w ECT polega na wyznaczeniu rozkładu przenikalności elektrycznej na podstawie zmierzonych pojemności W rekonstrukcji obrazów można posłużyć się metodami rozwiązującymi problem nieliniowy albo przybliżony model liniowy.

W przypadku przybliżenia liniowego nieznaną wartość wektora wyznacza się na podstawie wektora pomiarów przy założeniu wartości modelu W rzeczywistości macierz wrażliwości nie jest znana, gdyż wyznaczenie jej współczynników wymaga znajomości rozkładu przenikalności Zakładając, że szukany rozkład przenikalności elektrycznej nieznacznie różni się od rozkładu równomiernego, możliwe jest zastosowanie w modelu liniowym przybliżenia macierzy wrażliwości wyznaczonej dla rozkładu równomiernego.

Dodatkowo problem odwrotny jest problemem źle postawionym. Ze względu na małą liczbę elektrod liczba pomiarów jest mała, mniejsza od liczby punktów przestrzeni, w których wyznaczana jest wartość przenikalności. Zwyczajowo przyjmuje się, że matryca obrazu reprezentującego dwuwymiarowy rozkład przenikalności nie powinna być mniejsza niż 32×32 piksele. W takiej sytuacji układ równań jest niedookreślony, a macierz wrażliwości nie jest kwadratowa. Ponieważ w modelu liniowym występującym w ECT ten układ równań jest nieoznaczony i ma wiele rozwiązań, konieczne jest zastosowanie dodatkowego kryterium pozwalającego uzyskać unikalne rozwiązanie[10].

Podstawową metodą rozwiązanie tego problemu jest wybór rozwiązania o najmniejszej normie:

spełniającego równanie Rozwiązanie tego problemu optymalizacyjnego dane jest wzorem (prawa pseudoodwrotność Moore’a-Penrose’a):

Bardzo dużą różnicę wrażliwości pomiarów dla par elektrod sąsiednich i par elektrod przeciwległych (bardzo duży stosunek wartości szczególnych macierzy ) powoduje, że układ jest źle uwarunkowany numerycznie. W takim przypadku małe zmiany wartości pomiarowych mogą powodować bardzo duże zmiany w estymowanym obrazie. W przypadku, w którym rozwiązanie jest nieciągłą funkcją danych wejściowych konieczne jest dodanie dodatkowego warunku, który stabilizuje rozwiązanie, co matematycznie nazywa się regularyzacją. Zastosowanie metod regularyzacji, takich jak metoda Tichonowa albo metoda odcięcia wartości szczególnych (osobliwych) (ang. truncated singular value decomposition – TSVD), pozwala uzyskać dobre rezultaty pod warunkiem wybrania prawidłowej wartości parametru regularyzacji. Wartość tego parametru decyduje, jak bardzo rekonstruowany obraz będzie wygładzony. Wyznaczenie wartości parametru regularyzacji może być oddzielnym problemem optymalizacyjnym. Znane metody wyznaczania optymalnej wartości parametru regularyzacji, które były proponowane dla ECT to metoda krzywej L (ang. L-curve) i metoda walidacji krzyżowej (ang. cross validation method – CGV)[16].

Praktyczne zastosowanie pseudoodwrotności w ECT jest ograniczone wrażliwością tej metody na odstępstwo modelu numerycznego (macierzy wrażliwości) od modelu rzeczywistego. Dlatego w praktyce najczęściej stosowaną metodą rekonstrukcji obrazów jest metoda linowej projekcji wstecznej (ang. linear back-projection – LBP)[10]. W metodzie tej zakłada się przybliżenie macierzy odwrotnej macierzą transponowaną:

gdzie jest macierzą o wartościach znormalizowanych w kolumnach (suma współczynników każdej kolumny musi mieć jednakową wartość dla wszystkich kolumn). Normalizacja koryguje nierównomierność wrażliwości w przestrzeni wewnątrz sondy, tak aby suma wrażliwości wszystkich pomiarów dla pikseli w środku obrazu była taka sama, jak dla pikseli znajdujących się na brzegu obrazu. Algorytm LBP nie pozwala uzyskać obrazów o dobrej jakości, ze względu na mocne duże rozmycie obrazu[10].

W ECT znalazła również zastosowanie metoda optymalizacji średniokwadratowej, gdzie rozwiązaniem problemu liniowego jest wektor dany równaniem:

Minimalizacja tej normy może zostać przeprowadzona za pomocą algorytmu iteracyjnego. Najczęściej stosowanym algorytmem jest algorytm Landwebera, który jest prostą wersją algorytmu najszybszego spadku[10].

W przypadku rekonstrukcji nieliniowej stosuje się iteracyjny algorytm Levenberga-Marquardta[10].

Zastosowanie

Tomografię pojemnościową stosuje się w badaniach mieszanin materiałów elektrycznie nieprzewodzących o różnej przenikalności elektrycznej. ECT umożliwia wizualizację i ocenę hydrodynamiki przepływów wielofazowych[17], które odgrywają istotną rolę w procesach technologicznych przemysłu chemicznego[18], petrochemicznego[19], spożywczego[20].

Tomografia procesowa znajduje zastosowanie w przemyśle energetycznym między innymi przy monitorowaniu wodno-parowych układów chłodzenia w reaktorach jądrowych[21], jak również w przemyśle farmaceutycznym przy kontrolowaniu pracy reaktorów fluidalnych[22].

Dzięki bardzo wysokiej liczbie uzyskiwanych obrazów na sekundę tomografia pojemnościowa używana jest do obrazowania procesów spalania zarówno w kotłach energetycznych, jak i w komorach spalania silników lotniczych i rakietowych[23].

Zobacz też

  • Elektryczna tomografia impedancyjna
  • Tomografia procesowa
  • Metoda elektrooporowa

Przypisy

  1. J. Kryszyn i inni, Architecture of EVT4 data acquisition system for electrical capacitance tomography, „Measurement”, 101, Elsevier, 2017, s. 28–39, DOI10.1016/j.measurement.2017.01.020 (ang.).
  2. a b S.M. Huang i inni, Capacitance-based tomographic flow imaging system, „Electron. Lett.”, 7, 24, 1988, s. 418–419, DOI10.1049/el:19880283 (ang.).
  3. W.Q. Yang, A.L. Stott, M.S. Beck, High frequency and high resolution capacitance measuring circuit for process tomography, „IEE Proceedings – Circuits, Devices and Systems”, 3, 141, 1994, s. 215–219, DOI10.1049/ip-cds:19941019 (ang.).
  4. D. Wanta i inni, Numerical Evaluation of Complex Capacitance Measurement Using Pulse Excitation in Electrical Capacitance Tomography, „Electronics”, 11, MDPI, 2022, DOI10.3390/electronics11121864 (ang.).
  5. Y.D. Jiang, M. Soleimani, Capacitively Coupled Electrical Impedance Tomography for Brain Imaging, „Transactions on medical imaging”, 38, IEEE, 2019, s. 2104–2113, DOI10.1109/TMI.2019.2895035 (ang.).
  6. a b W.T. Smolik i inni, Single shot high voltage circuit for electrical capacitance tomography, „Measurement Science & Technology”, 28 (2), 2016, s. 1–12, DOI10.1088/1361-6501/aa50e1 (ang.).
  7. W.Q. Yang, T.A. York, New AC-based capacitance tomography system, „IEE Proceedings – Science, Measurement and Technology”, 1, 146, 1999, s. 47–53, DOI10.1049/ip-smt:19990008 (ang.).
  8. A. Plaskowski i inni, Imaging Industrial Flows: Applications of Electrical Process Tomography, CRC Press, 1995, ISBN 978-0-7503-0296-8 (ang.).
  9. J. Kryszyn, D. Wanta, W.T. Smolik, Gain Adjustment for Signal-to-Noise Ratio Improvement in Electrical Capacitance Tomography System EVT4, „IEEE Sensors Journal”, 2017, DOI10.1109/JSEN.2017.2744985 (ang.).
  10. a b c d e f g W. Smolik, Rekonstrukcja obrazów w elektrycznej tomografii pojemnościowej, Warszawa: Oficyna Wydawnicza Politechniki Warszawskiej, 2013, ISBN 978-83-7814-097-9.
  11. R. Banasiak i inni, Wizualizacja 4D ECT w czasie rzeczywistym przemysłowych procesów przepływu grawitacyjnego materiałów sypkich, „Automatyka”, 12 (3), 2008, s. 863–869 (ang.).
  12. W. Warsito, Q. Marashdeh, L.-S. Fan, Electrical Capacitance Volume Tomography, „IEEE Sensors Journal”, 1, 2007, s. 525–535, DOI10.1109/jsen.2007.891952 (ang.).
  13. L. Peng i inni, Evaluation of Effect of Number of Electrodes in ECT Sensors on Image Quality, „IEEE Sensors Journal”, 5, 12, 2012, s. 1554–1565, DOI10.1109/JSEN.2011.2174438 (ang.).
  14. D. Wanta i inni, A Finite Volume Method using a Quadtree Non-Uniform Structured Mesh for Modeling in Electrical Capacitance Tomography, „Proceedings of the National Academy of Sciences, India Section A: Physical Sciences”, Springer, 2021, DOI10.1007/s40010-021-00748-7 (ang.).
  15. D.B. Geselowitz, An application of Electrocardiographic Lead Theory to Impedance Plethysmography, „IEEE Trans. Biomed. Eng.”, BME-18 (1), 1971, s. 38–41, DOI10.1109/TBME.1971.4502787 (ang.).
  16. M. Soleimani, W.R.B. Lionheart, Nonlinear image reconstruction in electrical capacitance tomography using experimental data, „Meas. Sci. Technol.”, 16, 2005, s. 1987–1996, DOI10.1088/0957-0233/16/10/014 (ang.).
  17. A.J. Jaworski, T. Dyakowski, Application of electrical capacitance tomography for measurement of gas–solid flow characteristics in a pneumatic conveying system, „Measurement Science and Technology”, 12, 2001, s. 1109–1119, DOI10.1088/0957-0233/12/8/317 (ang.).
  18. T. Rymarczyk i inni, Electrical Capacitance Tomography and Optical Detection in Quality Control System, „Przegląd Elektrotechniczny”, 93 (12/2017), DOI10.15199/48.2017.12.53, ISSN 0033-2097 (ang.).
  19. T.A. York, Status of electrical tomography in industrial applications, „Journal of Electronic Imaging”, 10 (3), 2001, DOI10.1117/1.1377308 (ang.).
  20. M.S. Beck i inni, Principles and Industrial Applications of Electrical Capacitance Tomography, „Measurement and Control”, 30 (7), 1997, DOI10.1177/002029409703000702 (ang.).
  21. A. Dupré, G. Ricciardi, S. Bourennane, Development of a prototype device for imaging two-phase flows via Electrical Impedance Tomography, „Specialist Workshop on Advanced Instrumentation and Measurement Techniques for Nuclear Reactor Thermal Hydraulics (SWINTH-2016)”, Livorno, Italy, 17 czerwca 2016 (ang.).
  22. H. Wang, W. Yang, Scale-up of an electrical capacitance tomography sensor for imaging pharmaceutical fluidized beds and validation by computational fluid dynamics, „Measurement Science and Technology”, 22 (10), 2011, DOI10.1088/0957-0233/22/10/104015 (ang.).
  23. A. Hunt, Industrial Applications of High-speed Electrical Capacitance Tomography, „Proceedings of 9th WORLD CONGRESS ON INDUSTRIAL PROCESS TOMOGRAPHY”, Bath, UK, 6 października 2018, ISBN 978-0-853-16356-5 (ang.).

Linki zewnętrzne

Media użyte na tej stronie

ModelSondyECT.svg
Autor: Damian Wanta, Licencja: CC BY-SA 4.0
Schemat przykładowej sondy pomiarowej wykonanej z użyciem rury PCW
EctImageMovie.ogg
Autor: Damian Wanta, Licencja: CC BY-SA 4.0
Sekwencja obrazów ECT - obiekt wykonany ze szkła akrylowego, obracający się wewnątrz sondy pomiarowej
EVT4 i sonda.jpg
Autor: Damian Wanta, Licencja: CC BY-SA 4.0
System pomiarowy - Wielokanałowy tomograf pojemnościowy z podłączoną 16-elektrodową sondą
Obraz2i5dECT.jpg
Autor: Damian Wanta, Licencja: CC BY-SA 4.0
Złożony obraz objętościowy
ElekrodyECT.jpg
Autor: WolfgangNihil, Licencja: CC BY-SA 4.0
Pierścień elektrod pojemnościowych zamocowanych do rury pcw
SondaFantom2D.gif
Autor: Damian Wanta, Licencja: CC BY-SA 4.0
Obrazowanie dynamiczne w Elektrycznej Tomografii Pojemnościowej.
Po lewej stronie animacja przedstawiająca siedem cylindrycznych obiektów przemieszczających się wzdłuż sondy pomiarowej. Po prawej stronie seria obrazów przekroju sondy zrekonstruowanych metodą linowej projekcji wstecznej z następujących po sobie pomiarów pojemności. Kolor reprezentuje wartość względnej przenikalności elektrycznej ε w danym punkcie przestrzeni. Animacja stworzona z użyciem oprogramowania Blender oraz Matlab.