Page 24 - Rev_SEA_50_1_2_web
P. 24

Efecto del método de definición de las propiedades acústicas de cráneo humano…






             Comenzando por el modelo heterogéneo de densi-    donde el foco no llega a las tonalidades rojizas corres-
          dad aparente, se emplea la información obtenida a partir   pondientes a los 0.5 MPa (que es el umbral de apertura
          del CT-scan y se aplica la expresión siguiente para obte-  de la barrera hematoencefálica) alcanzando en torno a
          ner los mapas de velocidad de propagación, densidad y   los 0.32 MPa de presión máxima focal, como es obser-
          absorción:                                           vable también en las curvas que cortan cada eje indivi-
                                           (
                 c x, y, z) = c 0 + c max − c 0)⋅ ρ CT x, y, z)  dualmente (Figura 2 g,h,i).
                  (
                               (
                                                        (5)       También se han estudiado las características del foco
                               (
                 ρ x, y, z) = ρ 0 + ρ max − ρ 0)⋅ ρ CT x, y, z)
                  (
                                            (
                                                               acústico generado en el interior de la cavidad craneal:
          donde ρ CT  es la densidad aparente obtenida a partir del   atenuación del máximo del foco, desviación axial y des-
          CT-scan normalizada, ρ 0 es la densidad del agua, ρ max es   viaciones transversales. Bien, si aplicamos un requisito
          la densidad máxima en el cráneo, c  es la velocidad de   fundamental en este tipo de estudios de propagación
                                         0
          propagación en agua, y c max  es la velocidad de propaga-  transcraneal que asegura una localización precisa del
          ción máxima en el cráneo.                            foco, es decir, permitiéndonos cometer un error de como
             Otro modelo heterogéneo es el proporcionado por el   máximo un milímetro en la colocación del foco, el método
          toolbox k-Wave de Matlab, basado en los datos experi-  se considera preciso y factible. Por tanto, se establece
          mentales de Schneider junto con la relación empírica de   este criterio, y tal y como se observa en la Tabla 2, los
          Mast. Esta herramienta proporciona la función “houns-  modelos heterogéneos ofrecen un foco de calidad acús-
          field2density”, que convierte las unidades Hounsfield del   tica muy similar, mientras que el modelo homogéneo, y
          CT-scan en unidades de impedancia acústica.          muy importante para la zona del cráneo de incidencia del
                                                               haz elegida, ofrece una mayor pérdida energética y el
             Por último, el modelo homogéneo se ha definido de   foco presenta una mayor atenuación (en torno a 4 dB
          manera muy simple, ya que en lugar de asignar valores a   mayor), la desviación axial es muy notable (el foco se
          cada elemento del CT-scan, se ha asignado el mismo   acerca al transductor unos 4 mm), mientras que la des-
          valor de velocidad de propagación o densidad a todos   viación transversal entra dentro del umbral de error por
          los elementos del mapa craneal.                      debajo del milímetro, por lo que se considera similar a los
             La Tabla 1 muestra los valores de densidad y veloci-  casos heterogéneos.
          dad del sonido utilizados en el modelo heterogéneo de
          densidad aparente y en el modelo homogéneo, además   4.  Discusión
          de los parámetros de la ley de potencia para la absorción
          acústica.                                               Los dos modelos heterogéneos se comportan de la
                                                               misma manera, ofreciendo resultados prácticamente
          3.  Resultados                                       idénticos, mientras que el modelo homogéneo presenta
                                                               mayores pérdidas y una desviación axial muy notable.
             Los tres métodos estudiados se han comparado me-  Esto ocurre concretamente para la zona de incidencia
          diante la distribucíon del campo acústico obteniendo di-  sobre el cráneo estudiada (zona parietal superior), pues-
          versos cortes en axial y transversal de la presión máxima   to que en el caso heterogéneo la densidad media de esa
          simulada numéricamente. Los dos métodos heterogé-    zona es menor y por tanto ofrece menores pérdidas en
          neos son idénticos (Figura 2 a,b,d,e,g,h,i); los mapas   la propagación. La gran desviación axial, acercando el
          del campo acústico en 2D (Figura 2 a,b,d,e) son prác-  foco unos 4 mm hacia el transductor se debe a que la
          ticamente iguales, observando también y con más deta-  curvatura craneal en el caso homogéneo, al presentar
          lle el mismo resultado en las curvas 1D (Figura 2 g,h,i).   una impedancia acústica mayor, concentra el haz en una
          El método que difiere de los demás es el homogéneo,   posión del espacio más cercana a dicha curvatura, el
          presentando menor energía en la zona focal, aunque sí   efecto convergente de esa pared craneal aumenta cuan-
          que preserva la forma del campo acústico, apreciándose   to mayor sea su impedancia acústica, cuanto más duro
          tanto en las distribuciones 2D del campo (Figura 2 c,f),   sea ese material.

                                           Tabla 1.  Parámetros acústicos del estudio numérico
                                      Densidad (kg/m )     Velocidad (m/s)        Ley de potencia
                                                 3
                        Medio       Máxima     Media    Máxima     Media     α  (db/cm/MHz )    y
                                                                                        Y
                                                                              0
                    Agua (37º)        –        1000        –       1524  (1)  139.36.10 –5 (1)  2
                    Cráneo           2482      1617      3171       2203        11.17  (2)     1.1  (2)
                   (1) Valores obtenidos de funciones de la herramienta k-Wave; (2) [17].

          [ 22]  revista de acústica  |  Vol. 50  |  N.  1 y 2
                                            os
   19   20   21   22   23   24   25   26   27   28   29