FR 3087642 A1 20200501 FR 1859826 A 20181024 Procédé et système d'analyse spectrale et de détermination d'un marqueur permettant d'assurer la sécurité d'interventions d'ultrasons thérapeutiques La présente invention concerne un procédé d'analyse spectrale de la réponse acoustique d'un tissu biologique et de détermination d'un marqueur de sécurité permettant d'assurer la sécurité d'interventions d'ultrasons thérapeutiques. L'invention concerne également un système correspondant à l'analyse spectrale de la réponse acoustique d'un tissu biologique et de détermination d'un marqueur de sécurité permettant d'assurer la sécurité d'interventions d'ultrasons thérapeutiques. Le tissu biologique est un tissu biologique mou vascularisé quelconque, par exemple un tissu compris dans l'ensemble des tissus du cerveau, du foie, du coeur, des muscles, des seins, des reins, des yeux, de la thyroïde, de la prostate, de l'utérus, des tendons, du pancréas, de la peau, et de préférence un tissu cérébral. Malgré l'accroissement du nombre de médicaments actifs et l'avènement des thérapies ciblées en cancérologie, les progrès thérapeutiques pour les maladies cérébrales (cancer inclus) restent encore modestes. L'un des obstacles majeurs réside dans l'absence de délivrance spécifique et contrôlée des molécules thérapeutiques au sein des tissus. En effet, les parois des vaisseaux sanguins cérébraux forment une barrière endothéliale très efficace appelée barrière hémato-encéphalique. Elle limite le passage de molécules du sang vers les cellules à traiter. Les méthodes actuelles pour administrer les agents thérapeutiques sont invasives, non localisées ou présentent un risque élevé pour le patient. En outre, la libre circulation des substances thérapeutiques dans l'organisme provoque des effets indésirables sur les tissus sains. L'acheminement efficace, spécifique et localisé de molécules thérapeutiques constitue donc un défi majeur. Depuis 2000, de nombreuses études ont démontré que les ultrasons focalisés peuvent être utilisés pour accomplir cette tâche. Combinés à l'injection par voie intraveineuse de microbulles de gaz, les ultrasons peuvent être utilisés afin d'induire une ouverture localisée et réversible des barrières biologiques. En effet, les forces mécaniques (i.e., micro-flux et oscillations) 2 issues des interactions bulles-ultrasons (cavitation) fragilisent la barrière et favorisent le passage de molécules vers le tissu cérébral en général, et en particulier la région pathologique que l'on cherche à traiter si cette dernière est correctement visée par le faisceau d'ultrasons. La durée pouvant aller 5 jusqu'à plusieurs heures, l'efficacité de l'ouverture de la barrière en termes de passage de larges molécules ainsi que la sécurité de la technique peuvent être contrôlées en modifiant les paramètres ultrasonores. La « dose de cavitation » délivrée joue un rôle prépondérant dans l'efficacité et la sécurité de cette technologie. La pression acoustique au sein 10 du tissu traité doit être suffisante pour entraîner une oscillation maîtrisée des microbulles (régime de cavitation stable) et engendrer une perméabilisation réversible et non lésionnelle des parois vasculaires. En revanche, soumises à de trop fortes pressions acoustiques, les microbulles entrent alors dans un régime de cavitation inertielle impliquant 15 des phénomènes physiques localement violents (i.e., onde de chocs, micro- jet, implosion locale de la bulle) pouvant conduire à une détérioration des tissus et en l'apparition d'effets secondaires graves (e.g., inflammation, hémorragie). L'écart entre dose efficace et dose lésionnelle est faible et nécessite 20 de développer de nouvelles méthodes précises de dosimétrie in situ. Le contrôle de la dose de cavitation en temps réel est donc un enjeu majeur à l'heure où cette technologie entame son passage en essais cliniques. Dans le cadre de l'ouverture de la barrière hémato-encéphalique par ultrasons, l'objectif consiste à conserver un haut degré de cavitation stable 25 (efficacité) pendant le traitement tout en maintenant la cavitation inertielle à un faible niveau (sécurité). Lors d'une thérapie par ultrasons transcrâniens, l'hétérogénéité du crâne peut avoir une influence indésirable sur l'efficacité et la sécurité de la technique. En effet, l'épaisseur du crâne variant selon les zones, l'atténuation du 30 faisceau ultrasonore s'en trouve modifiée. L'amplitude de l'onde ultrasonore peut aisément varier d'un facteur 2 d'un point à l'autre chez l'Homme. Pour les gros animaux et plus particulièrement les primates non humains, la présence de tissus (e.g.. des muscles) entre la peau et le crâne peut également modifier le faisceau ultrasonore. 3 Ainsi, lorsque les ultrasons sont transmis au travers une région plus épaisse qu'attendu, l'amplitude de l'onde ultrasonore dans la zone d'intérêt sera plus faible et le traitement potentiellement inefficace (i.e., la pression acoustique ne sera pas suffisante pour permettre une oscillation notoire des 5 microbulles). Au contraire, si la zone est moins épaisse, l'amplitude de l'onde ultrasonore sera sous-estimée et des problèmes de sécurité résultant de la cavitation inertielle des bulles pourront se produire. La détection passive de cavitation peut être réalisée à l'aide de 10 transducteurs ultrasonores disposés autour de la zone traitée. Ces capteurs permettent de mesurer la réponse fréquentielle issue des microbulles et ainsi de déterminer en temps réel le régime de cavitation induit, comme décrit dans l'article de Coussios et al., intitulé « Role of acoustic cavitation in the delivery and monitoring of cancer treatment by high-intensity focused 15 ultrasound (HIFU) » et publié dans Int. J.Hyperthermia, March 2007, 23(2), pp.105-120. Un défi actuel consiste à discriminer la réponse des microbulles de celle des tissus environnants puis d'assurer que cette réponse correspond à une cavitation stable. Par convention, la dose de cavitation stable est 20 déterminée par le signal harmonique des microbulles : (n+1) fo (avec n E N* et fo la fréquence d'émission des ultrasons) alors que la dose de cavitation inertielle correspond à la mesure du signal émis à toutes les fréquences (signal large bande), comme le décrit l'article de Konofagou et al. (2012), intitulé « Ultrasound-Induced Blood-Brain Barrier Opening" et publié dans 25 Curr. Parma.Biotechnol. ;13(7), pp. 1332-1345. Toutefois, ces marqueurs de cavitation manquent de robustesse et de sensibilité. Des évènements indésirables (e.g., oedèmes) peuvent survenir malgré les précautions. D'autres indicateurs plus fiables et reproductibles sont donc nécessaires. 30 De manière connue, les microbulles sont généralement composées d'un gaz lourd (perfluorocarbone, hexafluorure de soufre) afin de réduire leur cinétique de dissolution dans le sang et d'augmenter ainsi la durée de l'examen échographique qui atteint plusieurs minutes. Ce gaz est entouré d'une coque dont la fonction est de protéger la bulle. L'épaisseur de la paroi 4 varie de quelques nm à plusieurs centaines de nm. Elle est généralement composée de protéines, de phospholipides, de surfactants ou de polymères. De récents travaux in vitro, décrits dans l'article de Shekar et al. (2014), intitulé « The delayed onset of subharmonic and ultraharmonic 5 emissions from a phospholipid contrast agent » et publié dans Ultrasound in Med.& Biol., 40(4), pp. 727-738, ont montré l'apparition d'oscillations spécifiques des microbulles lorsque ces dernières sont soumises à des conditions particulières. En effet, au cours du temps, le gaz contenu dans la bulle se diffuse dans le milieu environnant résultant sur un état de flambage 10 de son enveloppe (e.g., excès de lipide sur l'enveloppe de la bulle). Cet état est associé à l'apparition de fréquences spécifiques (sub-harmonique et ultra-harmonique : (2 n2+ 1) fo avec n E N et fo la fréquence d'émission des ultrasons) et peut conduire à la disparition de la microbulle. Ce phénomène est accentué par l'application de longues séquences 15 ultrasonores qui favorisent la diffusion du gaz, comme décrit dans l'article de O'Brien et al. (2013), intitulé « Surfactant shedding and gas diffusion during pulsed ultrasound through a microbubble contrast agent suspension » et publié dans J. Acoust. Soc. Am., 134, 1416-27, et le flambage de l'enveloppe de la microbulle comme décrit dans l'article de Kooiman et al. 20 (2017), intitulé "Focal areas of increased lipid concentration on the coating of microbubbles during short tone-burst ultrasound insonification" et publié dans PLoS ONE 12(7):e0180747. Cet état de déstabilisation des microbulles apparaît au cours de l'excitation ultrasonore et n'est pas à ce jour utilisé comme marqueur du 25 risque d'effets indésirables lors d'une thérapie ultrasonore assistée par microbulles de gaz. Il semblerait donc que l'on puisse en réalité définir trois régimes d'activité de microbulles : cavitation stable (rayonnement harmonique uniquement), cavitation inertielle (rayonnement harmonique, sub- 30 harmonique, ultra-harmonique et large-bande) et un régime intermédiaire correspondant à une déstabilisation de la coque des microbulles (rayonnement harmonique, sub-harmonique et ultra-harmonique). Aujourd'hui, les travaux mesurant l'activité de cavitation s'intéressent principalement aux composantes harmoniques et au spectre large-bande. 5 Actuellement, aucun consensus au niveau de la communauté scientifique n'est établi sur l'utilisation des composantes sub/ultraharmoniques comme marqueur d'activité stable ou de cavitation inertielle, ainsi qu'en atteste l'article de Haqshenas et al. (2015), intitulé « Multi- 5 resolution analysis of passive cavitation detector signais, et publié dans Journal of Physics, doi :10.1088/1742-6596/581/1/012004. Afin d'éviter les effets délétères liés à l'utilisation d'une trop forte pression acoustique, les travaux actuels déterminent un seuil de cavitation « tolérable » basé sur l'exploitation du signal large bande ou sur l'analyse 10 d'une ou plusieurs composantes fréquentielles spécifiques (e.g., harmoniques, sub-harmonique, ultra-harmoniques). Parmi ces travaux actuels, peuvent être cités les documents suivants : .- l'article de O'Reilly et al. (2012), intitulé « Blood-Brain Barrier : Realtime Feedback-controlled Focused Ultrasound Disruption by Using an 15 Acoustic Emissions-based Controller », et publié dans Radiology, 263(1), pp. 96-106 ; .- l'article de Tsai et al. (2016), intitulé "Real-time monitoring of focused ultrasound bloodbrain barrier opening via subharmonic acoustic emission detection: implementation of confocal dual-frequency piezoelectric 20 transducers", et publié dans Phys.Med. Biol., 61, pp. 2926-2946; .- l'article de Kamimura et al. (2018), intitulé "Feedback control of microbubble cavitation for ultrasound-mediated blood-brain barrier disruption in non human primates under magnetic resonance guidance", et publié dans J. Cereb. Blood Flow Metab. 1:271678X17753514; 25 .- la demande de brevet W02012042423 A1; et .- la demande de brevet W02008062342 A3. Qu'ils traitent l'information en temps-réel ou pas, avec une rétroaction ou pas sur le tir ultrasonore en cours, et qu'ils s'intéressent au spectre de rayonnement absolu ou relatif, ces travaux actuels ont comme point commun 30 de calculer des doses de cavitation tir après tir, une valeur pour chaque tir. Ces marqueurs de sécurité souffrent d'un manque de sensibilité et n'empêchent pas dans certains cas l'apparition d'effets délétères. Le problème technique est de fournir un marqueur de sécurité plus sensible, plus robuste et plus fiable qui permette d'éviter l'apparition d'effets 35 délétères sur un tissu soumis à une exposition d'ultrasons thérapeutiques, 6 et de fournir un procédé permettant de déterminer ce marqueur thérapeutique. A cet effet, l'invention a pour objet un procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, représentatif d'un état de 5 déstabilisation de microbulles contenues dans une zone d'un tissu biologique mou vascularisé, lesdites microbulles étant soumises à un signal d'excitation ultrasonore à une fréquence d'émission fo prédéterminée pour induire une ouverture localisée et réversible des barrières biologiques dans ladite zone, et ledit état de déstabilisation des microbulles étant dommageable pour le 10 tissu biologique, et ledit signal d'excitation ultrasonore étant formé par une séquence ultrasonore composée d'un nombre entier prédéterminé Nb, supérieur ou égal à 1, de train(s) d'onde, appelés « tir(s) ». Le procédé de détection et de détermination d'un marqueur de sécurité est caractérisé en ce qu'après le déclenchement de chaque tir Bb, b compris entre 1 et Nb, un 15 système d'analyse spectrale et de détermination d'un marqueur de sécurité : .- mesure régulièrement au cours du tir Bb en une série d'instants ta, l'évolution temporelle des raies spectrales correspondant aux fréquences subharmoniques et ultra-harmoniques du signal de réponse acoustique des 20 microbulles reçu, le signal de réponse reçu étant détecté par un détecteur passif de cavitation ayant une bande passante de détection prédéterminée, et .- détermine en la quantifiant l'évolution temporelle sur les instants ta d'un marqueur de sécurité, défini à chaque instant ta par un nombre MDDa 25 égal au rapport de la somme des aires des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal reçu de réponse acoustique des microbulles, sur la somme des aires des raies spectrales, mesurées au premier instant t1 et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du 30 signal de réponse acoustique des microbulles dans leur état initial. Suivant des modes particuliers de réalisation, le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité comprend l'une ou plusieurs des caractéristiques suivantes : .- la mesure à chaque instant ta des raies spectrales correspondant 35 aux fréquences subharmoniques et/ou ultra-harmoniques du signal de 7 réponse acoustique des microbulles reçu exploite le signal de réponse reçu dans une fenêtre d'observation wa contenant l'instant ta et incluse dans l'intervalle temporel de réception correspondant au tir considéré ; .- les fenêtres d'observation ou d'analyse wa sont adjacentes ou 5 séparées ou se recouvrent deux à deux partiellement ; .- le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité décrit ci-dessus comprend pour un tir donné Bb, b compris entre 1 et Nb, une première étape de mesure et de segmentation dans laquelle le signal de réponse acoustique des microbulles au train d'onde Bb, reçu et 10 mesuré par le détecteur passif de cavitation, est décomposé suivant un nombre entier prédéterminé k, supérieur ou égal à 2, de fenêtres temporelles wa, a variant de 1 à k, d'égales durées qui permettent de déterminer l'évolution des composantes fréquentielles au cours du train d'onde ultrasonore Bb ; 15 .- le nombre k de fenêtres et leurs tailles tw dépendent directement de la durée et de la fréquence d'excitation fo du tir ultrasonore Bb, la durée du tir ultrasonore étant comprise entre quelques microsecondes et plusieurs centaines de millisecondes ; .- la durée des fenêtres wa est comprise entre la durée de 8 cycles 20 du signal d'excitation et la moitié de la durée d'un tir, et/ou le nombre k de fenêtres wa est supérieur ou égal à 2 et inférieur ou égal au huitième du produit de la durée d'un tir TB par la fréquence d'excitation ultrasonore fo ; .- le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité décrit ci-dessus comprend pour un tir Bb donné, b compris entre 25 1 et Nb, une deuxième étape de calcul de spectres, exécutée après la première étape, dans laquelle pour chaque fenêtre wa du tir Bb, a variant de 1 à k, le système d' analyse spectrale et de détermination d'un marqueur de sécurité calcule le spectre fréquentiel de la portion du signal de réponse acoustique des microbulles au train d'onde Bb, contenue dans ladite fenêtre 30 wa ; .- la méthode de calcul des spectres fréquentiels utilise une transformée de Fourier ; .- le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, défini ci-dessus, comprend pour un tir Bb donné, b compris entre 35 1 et Nb, une troisième étape de calcul de l'évolution au cours du tir d'un 8 signal de cavitation s(a), exécutée après la deuxième étape, dans laquelle pour chaque fenêtre temporelle wa, a variant de 1 à k, le système d'analyse spectrale et de détermination d'un marqueur de sécurité calcule le signal de cavitation s(a) comme la somme des aires des raies spectrales, mesurées à 5 l'instant ta et correspondant aux fréquences subharmoniques et/ou ultra- harmoniques du signal de réponse acoustique des microbulles ; .- le nombre des composantes ultra-harmoniques et/ou subharmoniques considéré pour le calcul du signal de cavitation s(a) dépend de la bande passante du ou des transducteur(s) utilisé(s) pour la détection 10 de cavitation et formant le détecteur passif de cavitation ; .- l'amplitude des pics mesurés pour les composantes ultraharmonique(s) et/ou subharmonique(s) comprises dans la bande passante du détecteur passif de cavitation sont utilisés en complément ou en remplacement dans le calcul du signal de cavitation s(a) ; 15 .- le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, décrit ci-dessus, comprend pour un tir Bb donné, b compris entre 1 et Nb, une quatrième étape de calcul de l'évolution au cours du tir d'un marqueur de cavitation s(a), exécutée après la troisième étape (108), dans laquelle pour chaque fenêtre temporelle wa du tir Bb, a variant de 1 à k, le 20 système d'analyse spectrale et de détermination d'un marqueur de sécurité calcule un marqueur de sécurité , appelé « dose de cavitation MDDa », défini par un nombre MDDa, égal au rapport du signal de cavitation s(a) de la aième fenêtre wa sur le signal de cavitation s(1) de la première fenêtre temporelle w1, la dose de cavitation MDDa étant exprimé dans une échelle 25 linéaire ou logarithmique ; .- le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, décrit ci-dessus, comprend pour un tir Bb donné, b compris entre 1 et Nb, une cinquième étape de calcul de l'évolution au cours du tir d'un premier paramètre d'alerte A11(a) et/ou d'un deuxième paramètre d'alerte 30 Al2(a), exécutée après la quatrième étape, dans laquelle le premier paramètre d'alerte Al1 est mis dans un état actif lorsque le marqueur de sécurité MDDa dépasse une première valeur seuil de sécurité Th1 prédéterminée, et le deuxième paramètre d'alerte Al2 est mis dans un état actif lorsque le nombre de fois nf où le marqueur de sécurité MDDa a 9 dépassé la première valeur de seuil Th1 a dépassé une deuxième valeur de seuil Th2 prédéterminée ; .- le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, décrit ci-dessus, comprend pour un tir Bb donné, b compris entre 5 1 et Nb, une sixième étape, exécutée après la quatrième ou la cinquième étape, dans laquelle le système d'analyse spectral et de détermination d'un marqueur de sécurité transmet à un dispositif de commande et de contrôle qui intervient dans une boucle de rétroaction des paramètres de tir : les doses de cavitation MDDa évoluant au cours du tir fournies dans la quatrième lo étape, et/ou les états Al1 (a), Al2(a) du premier paramètre d'alerte et/ou du deuxième paramètre d'alerte déterminés dans la cinquième étape. L'invention a également pour objet un système d'analyse spectrale et de détermination d'un marqueur de sécurité, représentatif d'un état de déstabilisation de microbulles, contenues dans une zone d'un tissu 15 biologique, lesdites bulles étant soumises à un signal d'excitation ultrasonore émis à une fréquence d'émission fo prédéterminée pour induire une ouverture localisée et réversible des barrières biologiques dans ladite zone du tissu, et ledit état de déstabilisation des microbulles étant dommageable pour le tissu biologique, ledit signal d'excitation ultrasonore 20 étant formé par une séquence ultrasonore composée d'un nombre entier prédéterminé Nb, supérieur ou égal à 1 de train(s) d'onde, appelés « tir(s) Le système d'analyse spectrale et de détermination d'un marqueur de sécurité est caractérisé en ce qu'il est configuré pour, après le déclenchement de chaque tir Bb, b compris entre 1 et Nb : 25 .- mesurer régulièrement au cours du tir Bb en une série d'instants ta, l'évolution temporelle des raies spectrales correspondant aux fréquences subharmoniques et/ou ultra-harmoniques d'un signal de réponse acoustique des microbulles reçu, le signal de réponse reçu étant détecté par un détecteur passif de cavitation ayant une bande passante de détection 30 prédéterminée, et .- déterminer en la quantifiant l'évolution temporelle sur les instants ta de la série temporelle un marqueur de sécurité, défini à chaque instant ta par un nombre MDDa égal au rapport de la somme des aires des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences 35 subharmoniques et/ou ultra-harmoniques du signal reçu de réponse 10 acoustique des microbulles, sur la somme des aires des raies spectrales, mesurées au premier instant tl et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles. 5 Suivant des modes particuliers de réalisation, le système d'analyse spectrale et de détermination d'un marqueur de sécurité comprend l'une ou plusieurs des caractéristiques suivantes : .- à chaque instant ta, la mesure des raies spectrales correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de 10 réponse acoustique des microbulles reçu, exploite le signal de réponse reçu dans une fenêtre d'observation wa contenant l'instant ta et incluse dans l'intervalle temporel de réception correspondant au tir considéré ; .- les fenêtres d'observation ou d'analyse wa sont adjacentes ou séparées ou se recouvrent deux à deux partiellement et faiblement par leurs 15 bords. L'invention a également pour objet un produit ou programme d'ordinateur comprenant un ensemble d'instructions configurées pour mettre en oeuvre le procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, défini ci-dessus, lorsqu'elles sont chargées dans et exécutées 20 par un calculateur ou plusieurs calculateurs implémentés dans le système d'analyse spectrale et de détermination d'un marqueur de sécurité défini ci-dessus. L'invention a également pour objet un système d'assistance ultrasonore d'un traitement thérapeutique visant une zone d'un tissu 25 biologique mou vascularisé contenant des microbulles, comprenant : .- un dispositif d'excitation et d'émission d'une séquence thérapeutique de tir(s) d'excitation à une fréquence d'émission prédéterminée fo, focalisés sur la zone à traiter du tissu biologique, .- un capteur passif de cavitation pour détecter et mesurer la réponse 30 des microbulles contenues dans la zone en réponse aux tirs de la séquence, .- un système d'analyse spectrale et de détermination d'un marqueur de sécurité, représentatif d'un état de déstabilisation des microbulles, décrit ci-dessus, .- un dispositif de commande et de contrôle des paramètres de tir(s) 35 du dispositif d'excitation ultrasonore, 11 le capteur passif de cavitation, le système d'analyse spectrale et de détermination d'un marqueur de sécurité, le dispositif de commande et de contrôle et le dispositif d'excitation ultrasonore étant mis en série selon un chaîne de sorte à former une boucle rétroactive de sécurité. 5 L'invention sera mieux comprise à la lecture de la description de plusieurs formes de réalisation qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins sur lesquels : .- la Figure 1 est une vue schématique d'un système d'assistance ultrasonore d'un traitement thérapeutique visant une zone d'un tissu 10 biologique contenant des microbulles, un système d'analyse spectrale et de détermination d'un marqueur de sécurité étant intégré dans ledit système d'assistance ultrasonore ; .- la Figure 2 est une vue d'une séquence ultrasonore pour la thérapie avec microbulles de gaz (e.g., ouverture de la barrière hémato-encéphalique, 15 sono-perméabilisation) ; .- la Figure 3 est une vue du contenu fréquentiel de la réponse des microbulles en réponse à l'application de séquences de tirs ultrasonores ; .- la Figure 4 est une vue du procédé d'analyse spectrale et de détermination d'un marqueur de sécurité, ledit procédé étant exécuté lors de 20 chaque tir d'une séquence ; .- la Figure 5 est une vue de la décomposition en plusieurs fenêtres temporelles wa, a variant de 1 à k, du signal de cavitation reçu, issu des microbulles pendant un tir ultrasonore ; .- la Figure 6 est une vue d'un exemple de réponses spectrales avant 25 et après déstabilisation des microbulles au cours d'un tir d'excitation ultrasonore transmis à 0,5 MHz dans le cerveau de primate ; .- la Figure 7 est une vue de l'évolution de la dose de composantes ultra-harmoniques (1,5f0; 2,5f0; 3,5f0) issues du signal des microbulles (SonoVue) au cours d'un tir de 10 ms à 500 kHz appliqué dans le cerveau de 30 primate, un premier exemple aboutissant à l'apparition d'une hémorragie et un deuxième exemple sans survenue d'hémorragie ; .- les Figures 8A et 8B sont des vues de mesure du signal ultra-harmonique au cours des séquences ultrasonores chez le primate dans un cas non hémorragique pour la Figure 8A et un cas hémorragique pour la 35 Figure 8B. 12 Suivant la Figure 1, un système d'assistance ultrasonore 2 d'un traitement thérapeutique visant une zone 4 d'un tissu biologique 6 contenant des microbulles 8 comprend : .- un dispositif 12 d'excitation et d'émission d'une séquence 5 thérapeutique de trains d'onde d'excitation à une fréquence d'émission prédéterminée fo, appelés « tirs », focalisés sur la zone 4 à traiter du tissu biologique 6 ; .- un capteur passif de cavitation 14, réalisé à l'aide d'un ou plusieurs transducteurs de réception, pour détecter et mesurer la réponse des 10 microbulles contenues dans la zone 4 en réponse aux tirs de la séquence ; .- un système selon l'invention 16 d'analyse spectrale et de détermination d'un marqueur de sécurité, représentatif d'un état de déstabilisation des microbulles 8 ; et .- un dispositif 18 de commande et de contrôle des paramètres de 15 tir(s) du dispositif d'excitation 12. La fréquence d'excitation fo est choisie de sorte à permettre d'induire une ouverture localisée et réversible des barrières biologiques dans la zone traitée 4 du tissu 8. Le capteur passif de cavitation 14 a une bande passante de réception 20 qui dépend de la ou des bandes passantes des transducteurs électroacoustiques de réception. Le système selon l'invention 16 d'analyse spectrale et de détermination d'un marqueur de sécurité est réalisé, par exemple, par un ou plusieurs calculateurs électroniques. 25 Le système selon l'invention 16 d'analyse spectrale et de détermination d'un marqueur de sécurité est configuré pour, après le déclenchement de chaque tir Bb d'une séquence d'un nombre Nb prédéterminé de tirs, Nb étant un nombre entier supérieur ou égal à 1 et b étant un indice d'ordre dans la séquence compris entre 1 et Nb : 30 .- mesurer régulièrement au cours du tir Bb en une série d'instants ta, l'évolution temporelle des raies spectrales correspondant aux fréquences subharmoniques et ultra-harmoniques d'un signal de réponse acoustique des microbulles reçu, le signal de réponse reçu étant détecté par un détecteur passif de cavitation ayant une bande passante de détection prédéterminée, 35 et 13 .- déterminer en la quantifiant l'évolution temporelle sur les instants ta de la série temporelle un marqueur de sécurité, défini à chaque instant ta par un nombre MDD, égal au rapport de la somme des aires des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences 5 subharmoniques et/ou ultra-harmoniques du signal reçu de réponse acoustique des microbulles, sur la somme des aires des raies spectrales, mesurées au premier instant t1 et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles. 10 Le dispositif de 18 de commande et de contrôle des paramètres de tir(s) du dispositif d'excitation 12 est configuré pour recevoir en temps réel l'information du marqueur de sécurité et/ou des informations d'alerte et à partir de ces informations inhiber le dispositif d'excitation (arrêt du ou des tirs, modulation des paramètres de la séquence de tir(s)) lorsque la mise en 15 sécurité du tissu l'exige. Une boucle rétroactive de sécurité 20 automatique peut être mise en oeuvre par la mise en série selon une chaine du capteur de cavitation 14, du système d'analyse spectrale et de détermination d'un marqueur de sécurité 16, du dispositif de commande et de contrôle 18 des paramètres de tir(s), et 20 du dispositif d'excitation et d'émission 12 d'une séquence thérapeutique de tir(s) ultrasonore(s), comme illustré sur la Figure 1. Un signal d'excitation ultrasonore est formé par une séquence ultrasonore composée d'un nombre entier prédéterminé Nb, supérieur ou égal à 1, de train(s) d'onde, appelés « tir(s) »; et 25 Suivant la Figure 2, la séquence ultrasonore thérapeutique 52 est composée de plusieurs trains d'onde 54 (en anglais « burst »), appelés par la suite « tirs », répétés et séparés par des temps de repos. Typiquement, des tirs d'une durée comprise entre 3 et 10 ms sont répétés à une fréquence de répétition, de valeur comprise typiquement entre 5 et 10Hz, pendant une 30 durée de séquence ou d'exposition de 30 secondes à 10 minutes. Comme il a été déjà décrit en préambule, actuellement, la détection de cavitation (stable ou inertielle) est réalisée sur l'analyse spectrale de la réponse acoustique globale des microbulles à un train d'onde, c'est à dire la moyenne sur la durée totale de chaque tir. Comme mentionné 14 précédemment, les microbulles peuvent émettre des fréquences spécifiques (sub-harmoniques et ultra-harmoniques) liées à leur déstabilisation. Suivant la Figure 3, l'apparition dans le contenu spectral 62 d'un tir pris dans sa globalité de raies à des fréquences subharmoniques et ultra- s harmoniques est illustrée. Ici, une raie subharmonique 64 à 0,5 fo est observé et trois raies ultra-harmoniques 66, 68, 70, respectivement à 1,5 fo ; 2,5 fo ; 3,5 fo sont observées. L'apparition des raies subharmoniques et ultra-harmoniques est augmentée lorsque des tirs de longues durée, c'est-à-dire de plusieurs cycles d'onde à fo sont transmis. 10 Dans l'invention, les signaux de cavitation sont utilisés et exploités afin de déterminer un marqueur de sécurité en temps réel de la thérapie ultrasonore. L'invention repose sur la mesure de l'évolution des fréquences subharmoniques et ultra-harmoniques au cours d'un tir ultrasonore (un train 15 d'onde d'excitation émis) afin d'assurer la sécurité de la séquence thérapeutique. Les données ayant conduit à l'invention indiquent que l'apparition de raies à des fréquences subharmoniques et ultra-harmoniques dans le contenu spectral est un phénomène soudain qui se produit au cours du tir ultrasonore et qui perdure ensuite le temps de celui-ci. La détection de 20 ce phénomène peut être utilisée immédiatement pour arrêter le tir ultrasonore avant de reprendre le tir suivant une fois des nouvelles bulles intactes entrées dans le volume de tir. La séquence peut aussi être ajustée dynamiquement en terme d'amplitude ultrasonore et/ou de durée des tirs afin d'éviter la répétition du phénomène. 25 L'idée d'observer l'apparition de ce phénomène au cours du tir permet donc de définir précisément le temps à partir duquel ledit phénomène se produit et ainsi d'ajuster en conséquence la durée des tirs suivants de la séquence ultrasonore soit de façon fixe, soit préférentiellement de façon dynamique. 30 La mesure de ce phénomène au cours du tir suivant le procédé de l'invention est plus sensible que les outils de détection conventionnels basés sur l'analyse fréquentielle d'un tir complet. Dans le cas de l'analyse conventionnelle des données mesurées par le capteur de cavitation, ces phénomènes d'apparition soudaine peuvent se 35 produire sans être détectés. En effet, l'analyse conventionnelle revient à 15 réaliser une moyenne du contenu fréquentiel sur le tir complet. Il devient évident que dans le cas où le phénomène apparaît à la fin du tir, le moyennage masquera les composantes subharmoniques et ultra-harmoniques qui n'apparaîtront pas (ou trop peu) lors de l'analyse 5 fréquentielle. En plus de ce phénomène de dilution du paramètre utile par moyennage, il faut noter que le suivi tir à tir observe des réalisations indépendantes sur des nuages de bulles différents à chaque tir. Les stratégies conventionnelles, actuellement proposées, qui analysent le signal 10 de réponse sur le tir complet, utilisent comme référence un signal acquis avant injection des microbulles. Cette acquisition permet de prendre entre autre en considération le signal réfléchi par l'os et la propagation non linéaire des ultrasons dans le milieu. Toutefois, cette acquisition est réalisée plusieurs secondes à plusieurs minutes avant le traitement. Cette acquisition 15 est contraignante car elle doit être répétée pour toutes les amplitudes ultrasonores potentiellement évaluées lors de l'injection des microbulles et les positions du transducteur pendant la thérapie. Elle nécessite une parfaite immobilisation du système ultrason et du patient pendant tout le traitement. A contrario, dans l'invention, ce qui est observé est l'évolution du 20 signal des mêmes bulles au sein du tir, c'est-à-dire la dynamique de leur déstabilisation sous ultrasons, en comparant les bulles à elles-mêmes en début de tir. La référence étant prise sur les même bulles quelques cycles auparavant, avec comme seule variable leur exposition récente aux ultrasons, il est logique d'exhiber des paramètres plus physiquement 25 pertinents d'une telle analyse. Avec le procédé de l'invention, les difficultés évoquées précédemment, rencontrées avec les stratégies conventionnelles disparaissent, puisque le signal de référence est établi sur la première fenêtre w1 du tir pendant le traitement et l'application des ultrasons. Il en résulte un gain de temps et l'absence de perturbation. 30 De manière générale, un procédé d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'invention est mis en oeuvre par le système 16 d'analyse spectrale et de détermination d'un marqueur de sécurité. Le marqueur de sécurité selon l'invention est représentatif d'un état de 35 déstabilisation de microbulles contenues dans la zone 4 du tissu biologique 16 6, lesdites microbulles 8 étant soumises au signal d'excitation ultrasonore à la fréquence d'émission fo prédéterminée pour induire une ouverture localisée et réversible des barrières biologiques dans la zone 4, et ledit état de déstabilisation des microbulles 8 étant dommageable pour le tissu 5 biologique 6. Le signal d'excitation ultrasonore est formé par une séquence ultrasonore composée d'un nombre entier prédéterminé Nb, supérieur ou égal à 1 de tirs. Le procédé d'analyse spectrale et de détermination d'un marqueur de 10 sécurité est caractérisé en ce qu'après le déclenchement de chaque tir Bb, b compris entre 1 et Nb, le système d'analyse spectrale et de détermination d'un marqueur de sécurité : .- mesure régulièrement au cours du tir Bb en une série d'un nombre entier prédéterminée k d'instants ta, a variant de 1 à k, l'évolution temporelle 15 des raies spectrales correspondant aux fréquences subharmoniques et ultra- harmoniques du signal de réponse acoustique des microbulles reçu, le signal de réponse reçu étant détecté par un détecteur passif de cavitation ayant une bande passante de détection prédéterminée, et .- détermine en la quantifiant l'évolution temporelle sur les instants ta 20 d'un marqueur de sécurité, défini à chaque instant ta par un nombre MDDa égal au rapport de la somme des aires des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal reçu de réponse acoustique des microbulles, sur la somme des aires des raies spectrales, mesurées au premier instant t1 et 25 correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles. La mesure à chaque instant ta des raies spectrales correspondant aux fréquences subharmoniques et ultra-harmoniques du signal de réponse acoustique des microbulles reçu exploite le signal de réponse reçu dans une 3o fenêtre d'observation wa contenant l'instant ta et incluse dans l'intervalle temporel de réception correspondant au tir considéré. Les fenêtres d'observation ou d'analyse wa sont adjacentes ou séparées ou se recouvrent deux à deux partiellement et faiblement par leurs bords. 17 Il est à remarquer que la première fenêtre w1 qui sert de référence peut d'une part être légèrement décalée par rapport au début du tir pour attendre la stabilisation du signal et d'autre part la durée de la fenêtre w1 sur laquelle on calcule le spectre de référence peut être plus grande que la 5 durée d'une fenêtre d'observation consécutive ou la durée de la somme de plusieurs fenêtres consécutives. Cela permet d'obtenir un spectre de référence moins bruité. Suivant la Figure 4, un procédé 102 d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'invention, exécuté lors de o chaque tir Bb d'une séquence, b étant compris entre 1 et Nb, comprend un ensemble d'étapes 106, 108, 110, 112, 114 et 116. Dans une première étape 106 de mesure et de segmentation, le signal de réponse temporelle des microbulles 8 au tir ultrasonore Bb, reçu et détecté au préalable par le détecteur passif de cavitation 14, est décomposé 15 suivant un nombre entier prédéterminé k, supérieur ou égal à 2, de fenêtres temporelles wa, a variant de 1 à k, d'égales durées qui permettent de déterminer l'évolution des composantes fréquentielles au cours du train d'onde ultrasonore Bb. Le nombre k de fenêtres wa et leurs tailles tw vont directement 20 dépendre de la longueur TB et de la fréquence fo du tir ultrasonore utilisés pour la thérapie. Typiquement, les travaux actuels montrent une ouverture réversible de la barrière pour des tirs ultrasonores dont la longueur varie de quelques microsecondes à plusieurs dizaines de millisecondes. D'un point de vue 25 théorique, deux fenêtres temporelles de taille équivalentes suffisent pour réaliser la méthode. La première fenêtre temporelle w1 est alors utilisée comme référence, la seconde w2 permet d'observer l'évolution du signal reçu des microbulles. Toutefois, si le nombre de fenêtres k n'est pas suffisamment grand, il 30 sera impossible de déterminer avec précision, la précision étant déterminée par le nombre et la taille de la fenêtre, le temps d'apparition des fréquences subharmoniques et ultra-harmoniques, et donc d'ajuster efficacement les temps et amplitude de tir ultrasonore dans la séquence. D'autre part, l'utilisation de fenêtres temporelles trop longues induira un moyennage du 35 contenu fréquentiel du signal reçu rendant difficile la détection des 18 fréquences sub-harmoniques et ultra-harmoniques. Il est donc recommandé de subdiviser le signal récupéré en une multitude de fenêtres temporelles k consécutives, relativement courtes (si possible quelques dizaines de cycles ultrasonores) afin d'augmenter la sensibilité de la méthode. 5 Dans une deuxième étape 108 de calcul de spectres, exécutée après la première étape 106, pour chaque fenêtre temporelle wa (a, 1 8 cycles ultrasonores) afin d'éviter le recouvrement des composantes fréquentielles lors de l'analyse spectrale. 15 En résumé, et prenant en compte les exigences concernant la première étape 106 et la deuxième étape 108, la durée des fenêtres wa est comprise entre 8 cycles et la moitié de la durée TB d'un tir. Le nombre k de fenêtres wa est supérieur ou égal à 2 et inférieur ou égal au huitième du produit de la durée d'un tir TB par la fréquence d'excitation ultrasonore fo. 20 Dans une troisième étape 110 de l'évolution au cours du tir Bb d'un signal de cavitation s(a), exécutée après la deuxième étape 108, pour chaque fenêtre temporelle wa, a variant de 1 à k, le signal de cavitation s(a) est calculé à partir de la représentation du signal dans le domaine fréquentiel, fournie dans la deuxième étape 108, comme la somme des aires 25 des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles. Le nombre de composantes ultra-harmoniques à considérer peut varier selon les dispositifs et la bande passante du ou des transducteur(s) 30 utilisé(s) pour la détection de cavitation. C'est une question de rapport signal sur bruit de la détection des différentes composantes. L'analyse peut être réalisée sur une composante, de préférence fo/2 ou 3 f0/2 qui sont les plus intensément rayonnées par les bulles, ou plusieurs composantes, par exemple f0/2, 3 f0/2, 5 f0/2, 7 f0/2 et plus si le capteur passif de cavitation 14 35 le permet. 19 Il est à remarquer que l'amplitude des pics à ces composantes ultra-harmoniques et/ou subharmoniques peut également être considérée dans ce calcul du signal de cavitation, en complément ou remplacement de l'aire des raies spectrales sous la courbe de bandes spectrales. Les bandes spectrales 5 Af utilisées pour calculer ces composantes fréquentielles dépendent de la fréquence d'émission fo et de la fréquence d'échantillonnage du dispositif de mesure et de la durée des fenêtres temporelles wa. Les travaux concernant cette invention montrent que la méthode peut être appliquée efficacement pour des bandes spectrales Af allant de 0,04 * fo à 0,5 * fo avec une 10 préférence autour de 0,1* fo. La largeur de la bande spectrale doit être adaptée au nombre de cycle de la fenêtre temporelle afin de minimiser la mesure du bruit. Lorsque le nombre de cycles ne est grand, la bande spectrale Af doit être petite, et inversement lorsque le nombre de cycles ne est petit, la bande spectrale Af doit être grande. 15 Il est à remarquer qu'en variante on peut calculer séparément le signal de cavitation pour chacune des composantes ultra-harmoniques et/ou subharmoniques et considérer de manière associée des stratégies de rétroaction plus complexes dans un arbre de décision. Par exemple, selon la bande passante du capteur de cavitation et le bruit électronique ambiant 20 (plus présent à basse fréquence), si l'amplitude de la composante à f0/2 cc explose », c'est à dire est excessive, mais que l'amplitude de la composante à 3 fo/2 reste admissible, alors dans une tentative de rétroaction l'amplitude du faisceau est diminuée sans le couper puis si l'amplitude de la composante à 3 f0/2 finit par « exploser » alors le tir est vraiment coupé. 25 Dans une quatrième étape 112 de calcul de l'évolution au cours du tir Bb d'un marqueur de cavitation s(a), exécutée après la troisième étape 110, pour chaque fenêtre wa du tir Bb, a variant de 1 à k, le système d'analyse spectrale et de détermination d'un marqueur de sécurité calcule un marqueur de sécurité MDDa, appelé cc dose de cavitation MDDa » (en anglais 30 « Microbubble Destabilization Dose »), défini par un nombre MDDa, égal au rapport du signal de cavitation s(a) de la a-ième fenêtre wa sur le signal de cavitation s(1) de la première fenêtre temporelle w1, soit MDDa , la dose de cavitation MDDa étant exprimée dans une échelle linéaire ou logarithmique. 20 De manière générale, les fenêtres d'observation ou d'analyse wa sont adjacentes ou séparées ou se recouvrent deux à deux partiellement. L'exécution des première, deuxième, troisième et quatrième peut être mise en oeuvre dans enchainement en mode série ou en mode pipeline. 5 Par exemple, dans un premier mode de réalisation, l'analyse du signal de réponse temporelle est conduite sur des fenêtres temporelles adjacentes successives. Dans un deuxième mode de réalisation, l'analyse du signal de réponse temporelle est conduite sur des fenêtres temporelles glissantes avec 10 toujours une normalisation par rapport au signal de cavitation s(1) de la première fenêtre. Dans une cinquième étape 114 de calcul de l'évolution au cours du tir Bb d'un premier paramètre d'alerte Al1 et/ou d'un deuxième paramètre d'alerte Al2, exécutée par le système d'analyse spectrale et de détermination 15 d'un marqueur de sécurité après la quatrième étape 112, le premier paramètre d'alerte Al1 est mis dans un état actif lorsque le marqueur de sécurité MDDa dépasse une première valeur seuil de sécurité Th1 prédéterminée, et le deuxième paramètre d'alerte Al2 est mis dans un état actif lorsque le nombre de fois nf où le marqueur de sécurité MDDa a 20 dépassé la première valeur de seuil Th1 a dépassé une deuxième valeur de seuil Th2 prédéterminée. Ce nombre de fois nf peut être comptabilisé soit pour des tirs consécutifs ou en variante pour des tirs non nécessairement consécutifs. Dans une sixième étape 116, exécutée après la quatrième 112 ou la 25 cinquième étape 115, le système d'analyse spectrale et de détermination d'un marqueur de sécurité transmet les données suivantes au dispositif de commande et de contrôle qui intervient dans une boucle de rétroaction sur les paramètres de tir : .- les doses de cavitation MDDa évoluant au cours du tir fournies dans 3o la quatrième étape ; et/ou .- les états AI1(a), Al2(a) du premier paramètre d'alerte et/ou du deuxième paramètre d'alerte déterminés dans la cinquième étape. Par exemple, si la dose de cavitation MDDa dépasse un premier seuil de sécurité Th1 qui peut varier entre 6 et 22 dB selon les données et leurs 35 post traitements avec une valeur préférentielle autour de 8 dB, la séquence 21 ultrasonore de tirs peut être soit ajustée en temps réel par une boucle de rétrocontrôle, les paramètres de contrôle étant compris parmi la diminution de la pression acoustique, la réduction de la durée de tirs, et éventuellement l'espacement des tirs, ou arrêtée pour éviter tout effet indésirable lié aux 5 ultrasons. Il est aussi possible de définir un deuxième seuil de sécurité Th2 comme un nombre maximum de tirs tolérables pouvant dépasser la dose MDD seuil, i.e. le premier seuil de sécurité Thl , tirs consécutifs ou non consécutifs. A ce jour, aucun procédé et système d'assistance ultrasonore d'un traitement thérapeutique utilisant le marqueur MDDa de l'invention pour contrôler la sécurité d'une séquence ultrasonore thérapeutique n'est connu. Actuellement, l'analyse fréquentielle du signal de cavitation est effectuée sur la globalité du signal de réponse émis par les microbulles en réponse à chaque tir. 15 La décomposition du signal de réception de réponse à un tir sur au moins deux fenêtres temporelles permet d'observer l'apparition des raies correspondantes à composantes subharmoniques et/ou ultra-harmoniques au cours du tir. Ce phénomène d'apparition ne peut être ni observé ni quantifié correctement sans décomposition du signal de cavitation reçu sur 20 au moins deux fenêtres temporelles. De manière avantageuse, le bio-marqueur de sécurité MDDa selon l'invention, ses méthodes de calcul ainsi que son utilisation dans le cadre d'une boucle de rétrocontrôle sur les propriétés acoustiques du faisceau de tir permettent d'éviter la déstabilisation de la bulle engendrée par la 25 modification de son enveloppe et par la diffusion du gaz contenu initialement à l'intérieur de celle-ci, laquelle déstabilisation peut entraîner des effets physiques localement violents comme cause possible d'effets biologiques délétères. De manière avantageuse encore et bien que la motivation majeure de 30 l'invention soit l'amélioration de la sécurité du protocole de traitement, le fait d'éviter la déstabilisation des microbulles permet de les maintenir actives dans le flux sanguin plus longtemps et ainsi d'augmenter la durée totale pendant laquelle on peut solliciter lesdites microbulles pour ouvrir la barrière hémato-encéphalique BHE et ainsi d'augmenter l'efficacité du traitement. 22 L'utilisation du marqueur de sécurité MDD selon l'invention peut permettre d'allonger les tirs, augmenter l'amplitude du faisceau ou augmenter le rapport cyclique des tirs (en anglais « duty cycle ») tant que le seuil délétère établi, par exemple un seuil de 8 dB, n'est pas atteint en cours 5 de tir. En garantissant une moindre déstabilisation des microbulles, il allonge leur temps de circulation dans le sang et donc leur temps total d'interaction possible avec les parois de vaisseaux. Le marqueur de sécurité MDD selon l'invention peut donc servir d'outil d'optimisation d'efficacité, cet outil pouvant être combiné à d'autres outils plus classiques comme celui de la dose de 10 cavitation harmonique. Suivant la Figure 5, un exemple d'un signal de réponse acoustique 202 d'un milieu formé de microbulles, de tissu osseux et de tissus mous à un train d'onde, i.e un tir, de 10 ms transmis chez le primate à une fréquence fo de 0,5 MHz. La réponse 202 des microbulles au tir ultrasonore a été 15 décomposée en 75 fenêtres temporelles 208 équivalentes wa, a variant de 1 à 75, de 128 ps de durée. Il est à remarquer que la première fenêtre w1 ne débute pas nécessairement au temps de début du signal de réponse acquis. Ici sur la Figure 5, un délai de 160 ps a été pris. Cela permet au signal électronique de 20 mesure de se stabiliser et rend possible la comparaison avec les autres fenêtres. La Figure 6 montre un exemple de l'évolution 212 de la dose des composantes ultra-harmoniques 214, 216, 218 aux fréquences respectives 1,5f0; 2,5f0; 3,5f0 entre une première réponse spectrale 222 en trait noir 25 avant déstabilisation des microbulles et une deuxième réponse spectrale 224 en trait gris après déstabilisation des microbulles. Les première et deuxième réponses spectrales 222, 224 ont été mesurées dans deux fenêtres d'observations temporelles consécutives au cours d'un tir de 10 ms à 500 kHz appliqué dans le cerveau de primate pour l'ouverture de la barrière 30 hémato-encéphalique et à partir du signal de réponse des microbulles reçu par le capteur passif de cavitation. Le signal de réponse reçu est mesuré avec un transducteur acoustique, formant le capteur de cavitation, dont la bande passante de réception est centrée à 1,5 MHz. 23 L'intervalle de temps séparant la mesure de la première réponse spectrale et la mesure de la deuxième réponse spectral est égale à la durée d'une fenêtre d'observation temporelle, soit 128 ps. Le signal ultra-harmonique aux fréquences ultra-harmoniques 1,5f0, 5 2,5f0, 3,5f0 correspondant à une déstabilisation des microbulles apparait nettement dans le deuxième spectre par rapport au premier spectre avec une augmentation nette du signal de cavitation s(a) associé. Ainsi, à partir de l'analyse spectrale réalisée sur chaque fenêtre temporelle du signal de réponse au tir, la déstabilisation des microbulles peut lo être détectée. Suivant la Figure 7, l'évolution 252 de la dose de composantes ultra-harmoniques (1,5f0; 2,5f0 ; 3,5f0) issues du signal des microbulles (SonoVue) au cours d'un tir de 10ms à 500 kHz appliqué dans le cerveau de primate est illustrée, un premier exemple aboutissant à l'apparition d'une hémorragie 15 étant décrit par une première courbe 254 en trait noir, et un deuxième exemple sans survenue d'hémorragie étant décrit par une deuxième courbe 256 en trait gris. Les première et deuxième courbes 254, 256 d'évolution démontrent que l'apparition des composantes ultra-harmoniques est un phénomène 20 soudain se produisant à partir d'un certain temps en fonction des paramètres ultrasonores appliqués, du patient et du milieu soumis aux ultrasons. En effet, une excitation longue entraîne une oscillation répétée des bulles conduisant à leur potentielle déstabilisation à partir de plusieurs cycles d'excitation. Dans les conditions ultrasonores présentes, c'est-à-dire des 25 trains d'ondes de 10ms à 500 kHz chez le primate et à 650 kHz chez le rat, les résultats d'étude concernant l'invention indiquent que le temps moyen d'apparition du phénomène est de 5,4 ms dans le cerveau de primate et de 3,2 ms dans le cerveau de rat. En pratique, un seuil de MDD de 8 dB a été considéré comme marqueur d'un évènement de déstabilisation des 30 microbulles. Les Figures 8A et 8B montrent des exemples d'un cas (Figure 8A) où aucun effet délétère n'a été visualisé et d'un cas (Figure 8B) où une hémorragie a été observée dans le cerveau du primate. La déstabilisation des microbulles induite par la séquence ultrasonore 35 est caractérisée par l'apparition de fréquences ultra-harmoniques, en clair 24 sur la Figure 8B au cours du tir ultrasonore. La visualisation de ce phénomène spécifique est associée à l'observation d'une hémorragie chez le primate, confirmée par IRM (Imagerie par Résonance Magnétique). Pour ces exemples, sur une période d'excitation ultrasonore de 30s, la 5 fréquence d'apparition du phénomène est de 73,8% dans le cas hémorragique et de 0,7% dans le cas non hémorragique. L'apparition répétée de ce phénomène, c'est-à-dire plus de deux évènements pour des tirs consécutifs, est associée à la présence d'effets indésirables (oedèmes et hémorragies) chez l'animal. 10 Il est à noter que cet exemple se limite à l'étude des fréquences ultra- harmoniques (trois premières ultra-harmoniques) car la bande passante limitée en réception du capteur ultrasonore de cavitation ne permet pas une mesure fiable de la composante subharmonique dans cet exemple. Toutefois, l'utilisation d'un transducteur de détection de cavitation adapté à la 15 mesure du signal subharmonique permet de réaliser une étude similaire sur la composante subharmonique. Sur la Figure 8B, on peut mettre en évidence deux temps distincts : verticalement, le nombre de tirs avant que les bulles ne commencent à émettre des ultra-harmoniques pendant le tir. Puis horizontalement à partir 20 de ce temps, la durée au début du tir avant qu'elles ne commencent à en émettre. Ces deux temps sont importants et renseignent sur la qualité des bulles dans le milieu. 25 REVENDICATIONS.1- Système d'analyse spectrale et de détermination d'un marqueur de sécurité, représentatif d'un état de déstabilisation de microbulles (8) contenues dans une zone (4) d'un tissu biologique (6) mou vascularisé, lesdites microbulles(8) étant soumises à un signal d'excitation ultrasonore à une fréquence d'émission fo prédéterminée pour induire une ouverture localisée et réversible des barrières biologiques dans ladite zone (4), et ledit état de déstabilisation des microbulles (8) étant dommageable pour le tissu biologique (6), ledit signal d'excitation ultrasonore étant formé par une séquence ultrasonore composée d'un nombre entier prédéterminé Nb, supérieur ou égal à 1 de train(s) d'onde, appelés « tir(s) »; et le système d'analyse spectrale et de détermination d'un marqueur de sécurité étant caractérisé en ce qu'il est configuré pour, après le déclenchement de chaque tir Bb, b compris entre 1 et Nb : .- mesurer (106, 108) régulièrement au cours du tir Bb en une série d'instants ta, l'évolution temporelle des raies spectrales correspondant aux fréquences subharmoniques et ultra-harmoniques du signal de réponse acoustique des microbulles reçu, le signal de réponse reçu étant détecté par un détecteur passif de cavitation ayant une bande passante de détection prédéterminée, et .- déterminer (110, 112) en la quantifiant l'évolution temporelle sur les instants ta d'un marqueur de sécurité, défini à chaque instant ta par un nombre Mn), égal au rapport de la somme des aires des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal reçu de réponse acoustique des microbulles, sur la somme des aires des raies spectrales, mesurées au 3o premier instant t1 et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles dans leur état initial. 2- Système d'analyse spectrale et de détermination d'un marqueur 35 de sécurité selon la revendication 1, dans lequel la mesure à chaque instant 26 ta des raies spectrales correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles reçu exploite le signal de réponse reçu dans une fenêtre d'observation wa contenant l'instant ta et incluse dans l'intervalle temporel de réception 5 correspondant au tir considéré. 3- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon la revendication 2, dans lequel les fenêtres d'observation ou d'analyse wa sont adjacentes ou séparées ou se recouvrent deux à deux 10 partiellement. .4- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'une quelconque des 1 à 3, configuré pour un tir donné Bb, b compris entre 1 et Nb, 15 dans une première étape (104) de mesure et de segmentation, décomposer le signal de réponse acoustique des microbulles au train d'onde Bb, reçu et mesuré par le détecteur passif de cavitation, suivant un nombre entier prédéterminé k, supérieur ou égal à 2, de fenêtres temporelles wa, a variant de 1 à k, d'égales durées qui permettent de déterminer l'évolution des 20 composantes fréquentielles au cours du train d'onde ultrasonore Bb. .5- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon la revendication 4, dans lequel le nombre k de fenêtres et leurs tailles tw dépendent directement de 25 la durée et de la fréquence d'excitation fo du tir ultrasonore Bb; la durée du tir ultrasonore étant comprise entre quelques microsecondes et plusieurs centaines de millisecondes. .6- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon la revendication 4 ou la revendication 5, dans lequel la durée des fenêtres wa est comprise entre la durée de 8 cycles du signal d'excitation et la moitié de la durée d'un tir ; et/ou le nombre k de fenêtres wa est supérieur ou égal à 2 et inférieur ou égal au huitième du produit de la durée d'un tir TB par la fréquence d'excitation ultrasonore fo. 27 .7- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'une quelconque des revendications 4 à 6, configuré pour, un tir Bb étant donné avec b compris entre 1 et Nb, 5 dans une deuxième étape (106) de calcul de spectres, exécutée après la première étape (104), calculer pour chaque fenêtre wa du tir Bb, a variant de 1 à k, le spectre fréquentiel de la portion du signal de réponse acoustique des microbulles au train d'onde Bb, contenue dans ladite fenêtre wa. 10 .8- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon la revendication 7, dans lequel le calcul des spectres fréquentiels utilise une transformée de Fourier. 15 .9- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'une quelconque des revendications 7 à 8, configuré pour, un tir Bb étant donné avec b compris entre 1 et Nb, dans une troisième étape (108) de calcul de l'évolution au cours du tir d'un signal de cavitation s(a), exécutée après la deuxième étape (106), 20 calculer pour chaque fenêtre temporelle wa, a variant de 1 à k, le signal de cavitation s(a) comme la somme des aires des raies spectrales, mesurées à l'instant ta et correspondant aux fréquences subharmoniques et/ou ultra-harmoniques du signal de réponse acoustique des microbulles. 25 .10- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon la revendication 9, dans lequel le nombre des composantes ultra-harmoniques et/ou subharmoniques considéré pour le calcul du signal de cavitation s(a) dépend de la bande passante d'un ou plusieurs transducteur(s) utilisé(s) pour la 30 détection de cavitation et formant le détecteur passif de cavitation. .11- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'une des revendications 9 à 10, dans lequel l'amplitude des pics mesurés pour les composantes ultra-35 harmonique(s) et/ou subharmonique(s) comprises dans la bande passante 28 du détecteur passif de cavitation sont utilisés en complément ou en remplacement dans le calcul du signal de cavitation s(a). .12- Système d'analyse spectrale et de détermination d'un marqueur 5 de sécurité selon l'une quelconque des revendications 9 à 11, configuré pour, un tir Bb étant donné avec b compris entre 1 et Nb, dans une quatrième étape (110) de calcul de l'évolution au cours du tir d'un signal de cavitation s(a), exécutée après la troisième étape (108), calculer pour chaque fenêtre temporelle wa du tir Bb, a variant de 1 à k, un 10 marqueur de sécurité, appelé « dose de cavitation MDDa », défini par un nombre MDDa, égal au rapport du signal de cavitation s(a) de la a-ième fenêtre wa sur le signal de cavitation s(1) de la première fenêtre temporelle wl, la dose de cavitation MDDa étant exprimé dans une échelle linéaire ou logarithmique. 15 .13- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon la revendication 12, configure pour, un tir Bb étant donné; avec b compris entre 1 et Nb, dans une cinquième étape (112) de calcul de l'évolution au cours du tir 20 d'un premier paramètre d'alerte AI1(a) et/ou d'un deuxième paramètre d'alerte Al2(a), exécutée après la quatrième étape (110), mettre le premier paramètre d'alerte Al1 dans un état actif lorsque le marqueur de sécurité MDDa dépasse une première valeur seuil de sécurité Th1 prédéterminée, et 25 mettre le deuxième paramètre d'alerte Al2 dans un état actif lorsque le nombre de fois nf où le marqueur de sécurité MDDa a dépassé la première valeur de seuil Thl a dépassé une deuxième valeur de seuil Th2 prédéterminée. 30 .14- Système d'analyse spectrale et de détermination d'un marqueur de sécurité selon l'une quelconque des revendications 12 à 13, configuré pour, un tir Bb étant donné avec b compris entre 1 et Nb, dans une sixième étape (114), exécutée après la quatrième (110) ou la cinquième étape (112), transmettre à un dispositif de commande et de 29 contrôle qui intervient dans une boucle de rétroaction des paramètres de tir : .- les doses de cavitation MDDa(k) évoluant au cours du tir fournies dans la quatrième étape ; et/ou 5 .- les états AI1(a). Al2(a) du premier paramètre d'alerte et/ou du deuxième paramètre d'alerte déterminés dans la cinquième étape. .15- Système d'assistance ultrasonore d'un traitement 10 thérapeutique visant une zone (4) d'un tissu biologique (6) mou vascularisé contenant des microbulles (8), comprenant .- un dispositif (12) d'excitation et d'émission d'une séquence thérapeutique de tir(s) d'excitation à une fréquence d'émission prédéterminée fo, focalisés sur la zone (4) à traiter du tissu biologique (6), 15 .- un capteur passif de cavitation (14) pour détecter et mesurer la réponse des microbulles (8) contenues dans la zone (4) en réponse aux tirs de la séquence, .- un système d'analyse spectrale et de détermination (16) d'un marqueur de sécurité, représentatif d'un état de déstabilisation des 20 microbulles (8), défini selon l'une des revendications 1 à 14, .- un dispositif de commande et de contrôle (18) des paramètres de tir(s) du dispositif d'excitation ultrasonore (12), le capteur passif de cavitation (14), le système d'analyse spectrale et de détermination d'un marqueur de sécurité (16), le dispositif de commande 25 et de contrôle (18) et le dispositif d'excitation ultrasonore (12) étant mis en série selon un chaîne de sorte à former une boucle rétroactive de sécurité (20). 30 35