electrical impedance tomography inverse prob lems
play

Electrical Impedance Tomography, inverse prob- lems, material - PowerPoint PPT Presentation

Electrical Impedance Tomography, inverse prob- lems, material characterization & structural health monitoring Aku Sepp anen Department of Applied Physics University of Eastern Finland Kuopio, Finland Finnish Inverse Problems Summer


  1. Electrical Impedance Tomography, inverse prob- lems, material characterization & structural health monitoring Aku Sepp¨ anen Department of Applied Physics University of Eastern Finland Kuopio, Finland Finnish Inverse Problems Summer School 2019 UEF // University of Eastern Finland

  2. Inverse problems

  3. Linnanm¨ aki

  4. Contents Electrical Impedance Tomography(EIT) Ill-posedness of the EIT inverse problem EIT-imaging of concrete EIT-based sensing skin

  5. Electrical Impedance Tomography (EIT) In EIT electric currents I are applied to electrodes on the surface of the object and the resulting potentials V are measured using the same electrodes. The conductivity distribution σ = σ ( x ) is reconstructed based on the potential measurements.

  6. The forward model in EIT ∇ · ( σ ∇ u ) = 0 , x ∈ Ω u + z ℓ σ∂ u ∂ν = U ( ℓ ) , x ∈ e ℓ , ℓ = 1 , 2 , . . . , L e ℓ σ∂ u ∂ν d S = I ( ℓ ) , � ℓ = 1 , 2 , . . . , L σ∂ u x ∈ ∂ Ω \ ∪ L ∂ν = 0 , ℓ =1 e ℓ

  7. Finite element approximation of the EIT forward model FE-approximation of the complete electrode model ⇒ V = U ( σ ) where σ ∈ R N is a finite dimensional approximation of the conductivity.

  8. Solution of the forward problem

  9. Solution of the forward problem

  10. Solution of the forward problem

  11. Solution of the forward problem

  12. Solution of the forward problem

  13. Solution of the forward problem

  14. Solution of the forward problem

  15. Two different targets & electrode potentials

  16. Inverse problem of EIT

  17. Two different targets & electrode potentials

  18. Two different targets & electrode potentials

  19. Two different targets & electrode potentials

  20. Two different targets & electrode potentials

  21. Two different targets & electrode potentials

  22. Two different targets & electrode potentials

  23. Two different targets & electrode potentials

  24. Two different targets & electrode potentials

  25. Inverse problem of EIT

  26. Modeling errors in EIT Example 1: Modeling error due to unknown contact impedances z (true z = 1, assumed z = 0.01). Left: True conductivity distribution. Middle: EIT reconstruction based on correct model (z=1). Right: EIT reconstruction based on incorrect model (z=0.01).

  27. Modeling errors in EIT Example 2: Modeling error due to inaccuracy in the injected current I (error level in I: 0.5%). Left: True conductivity distribution. Middle: EIT reconstruction based on correct model. Right: EIT reconstruction based on incorrect model.

  28. Modeling errors in EIT Example 3: Modeling error due to unknown boundary shape. Left: Photograph of a target. Middle: EIT reconstruction based on correct geometry. Right: EIT reconstruction based on circular model geometry. ◮ Nissinen et al 2010

  29. Inverse problem of EIT

  30. Solution of the inverse problem of EIT

  31. Solution of the inverse problem of EIT Reconstructing the conductivity σ based on noisy observations V obs is an ill-posed inverse problem . The solution of the inverse problem is typically written in the form σ> 0 {� L n ( V obs − U ( σ )) � 2 + A ( σ ) } σ MAP = arg min The functional A ( σ ) models the prior information on the conductivity distribution σ . Non-linear, constrained optimization problem

  32. Iteration step 1 Left: estimated conductivity distribution. Right: Measured vs. computed potentials.

  33. Iteration step 2 Left: estimated conductivity distribution. Right: Measured vs. computed potentials.

  34. Iteration step 3 Left: estimated conductivity distribution. Right: Measured vs. computed potentials.

  35. Iteration step 4 Left: estimated conductivity distribution. Right: Measured vs. computed potentials.

  36. Iteration step 5 Left: estimated conductivity distribution. Right: Measured vs. computed potentials.

  37. Final estimate Figure: Left: Photo of the true target; Right: estimated conductivity distribution.

  38. The resolution of EIT is usually not very high...

  39. ”Blobology”

  40. ”Blobology” However, if feasible prior information on the resistivity is available, the resolution can be improved...

  41. Concrete The most extensively used construction material in the world About 7.5 cubic kilometers of concrete made each year In the United States, more than 55,000 miles of highways paved with concrete $ 35-billion industry. In the United States, 2 million workers Evaluation, repair and restoration: 35 % of the total volume work in building industry

  42. Concrete, need for evaluation/testing/monitoring On-site testing & evaluation Crack detection Prediction of rebar corrosion risk, etc... Material characterization in lab scale Evaluation of transport properties – esp. the ability of concrete to impede the ingress of water

  43. EIT imaging of concrete

  44. EIT imaging of 3D moisture flow in concrete

  45. ECT imaging of 3D moisture flow in concrete

  46. ECT imaging of 3D moisture flow in concrete

  47. ECT imaging of 3D moisture flow in concrete

  48. Results with X-ray CT...

  49. Imaging of cracks

  50. EIT-based sensing skin for damage detection Electrically conductive material (e.g. copper tape, CNT film, copper/silver paint) is applied on the surface of concrete The cracking of concrete breaks also the sensing skin Detecting of cracks in the surface material with EIT

  51. EIT-based sensing skin for damage detection We choose the painted sensing skin (Easy to apply & applicable to a large scale). 2D EIT imaging problem

  52. Case 1: Sensing skin on plexi-glass Sensing skin painted on plexi-glass 16 electrodes for EIT Synthetic cracks made by scratching the paint

  53. Case 1

  54. Case 1

  55. Case 1

  56. Case 1

  57. Case 1

  58. Case 1

  59. Case 1

  60. Case 1

  61. Case 1

  62. No blobology!

  63. Case 1: difference vs absolute reconstructions

  64. How? We fit homogeneous conductivity distribution σ ref to reference EIT data V ref Denote the discrepancy between V ref and the modeled data by ǫ ǫ = V ref − U ( σ ref ) This error is mostly due to inhomogeneity of the sensing skin. An approximative modeling error correction; observation model V = U ( σ ) + ǫ + n

  65. How? MAP estimate 0 <σ<σ ref { 1 2 � L n ( V − U ( σ ) − ǫ ) � 2 + A ( σ ) } σ MAP = arg min where A ( σ ) is a potential function related to a total variation prior � A ( σ ) = α �∇ σ � d r Ω A ( σ ) promotes sparsity of ∇ σ .

  66. Case 2: Notched concrete beam in 4-point bending

  67. Case 2: Notched concrete beam in 4-point bending

  68. Case 2: Notched concrete beam in 4-point bending

  69. Case 2: Notched concrete beam in 4-point bending

  70. Case 2: Notched concrete beam in 4-point bending

  71. Case 2: Photo vs. EIT reconstruction

  72. Case 2: Photo vs. EIT reconstruction

  73. Case 2: Photo vs. EIT reconstruction

  74. Case 2: Photo vs. EIT reconstruction

  75. Case 2: Photo vs. EIT reconstruction

  76. Case 2: reconstructions, denser FE mesh

  77. Temperature sensing experiment Sensing skin was exposed to temperature changes by contact with a heat source. Temperature of the heat source could be controlled within 2 ◦ C, when in contact with the temperature sensor. Reconstructed conductivities were converted to temperature maps based on an experimentally determined T vs. σ curve.

  78. Local temperature change 77 ◦ C

  79. Local temperature changes 37 ◦ C and 77 ◦ C

  80. H2020 project, Science for clean energy

  81. H2020 project, Science for clean energy

  82. Thank you!

Recommend


More recommend