Brevet FR2520117A1 — publié le 19830722

La présente invention concerne la spectroscopie par résonance magnétique nucléaire utilisant la transforma- tion de Fourier, et elle porte plus particulièrement sur la suppression ou la compensation de phénomènes d'oscillations acoustiques dans des données de résonance excitées par impul- sions. La pratique moderne en spectroscopie par résonance magnétique nucléaire (RMN) utilise une excitation par impul- sions des systèmes résonnants, suivie par l'acquisition de la réponse dans le domaine des temps, c'est-à-dire d'un signal. L'irnformation contenue dans le signal donne un spectre de fréquence par l'application d'une transformation de Fourier appropriée Un certain nombre d'effets parasites accompagnent la mise en oeuvre de la résonance magnétique nucléaire avec transformation de Fourier, et on peut minimiser les artefacts résultants dans le spectre par des mesures appropriées por- tant sur la cause de l'artefact, ou bien on peut souvent an- nuler l'artefact en exploitant la différence de propriétés de cohérence entre le signal de résonance nucléaire et l'ar- tefact. On sait tque des oolurants de haute fréquence circu- Lant dans des conducteurs situés dans des champs magnétiques statiques peuvent donner naissance a des ondes ultrasonores stationnaires qui, par un mécanisme inverse, produisent un signal ayant l'allure d'une réssonance Cet effet parasite a souvent été observé et il est étudié longuement par Fukushima et col dans la revue Journal of ïlagnetic Resonance, Vol 33, pages 199-205 ( 1979) Cet article décrit également des mesu- res visant à minimiser l'effet par le choix de matières et par la configuration du matériel. L'oscillation acoustique constitue un phénomène semblable à une résonance large qui traduit le couplage élevé entre le courant de haute fréquence et un réseau cristallin dans le champ magnétique L'onde stationnaire résultante est donc caractérisée par des temps de décroissance assez courts. Néanmoins, l'amplitude est élevée et des signaux de RM Nde courte durée peuvent être complètement perdus du fait de la présence de ce transitoire. On a déterminé dans l'art antérieur qu'il était possible de faire disparaître des spectres les résonances à courte durée de vie par le moyen simple consistant à interca- ler un retard entre l'excitation et la détection de la ré- senance On observe sans difficulté des signaux de RIA rela- tivement longs en présence d'une oscillation acoustique cour- te en retardant l'acquisition pendant une durée appropriée pour permettre l'atténuation de l'oscillation acoustique Ma- lheureusement, cette stratégie introduit un déphasage fonc- tion de la fréquence qu'il faudra finalement compenser et, d'autre part, le rapport signal à bruit est réduit par la dé- croissance du niveau du signal pendant la durée de retard. Dans une version plus élaborée de cette technique de temporisation, cn a utilisé une mesure d'écho de spin pour préserver pratiquement le signal de résonance pendant la dé- croissance de l'oscillation due à l'instrumentation Le dé- phasage introduit par le procédé simple de retard et acqui- sition est éliminé par la mesure d'écho de spin, mais on doit supporter un facteur d'atténuation du signal, égal à e-2 /T 2. En outre, des mesures d'écho de spin introduisent une oscil- lation acoustique avec l'impulsion d'inversion aussi bien qu'à avec l'impulsion d'observation, et le choix du temps de refocalisation doit nécessairement être suffisamment long pour permettre l'atténiation de l'oscillation qui correspond à l'impulsion d'inversion. L'invention a pour but de compenser la présence d'une oscillation acoustique superposée sur des signaux de résonance magnétique nucléaire excités par impulsions. Selon une caractéristique d'un premier mode de réa- lisation, on obtient un ensemble de mesures d'écho de spin dans lesquelles la phase relative d'impulsions d'inversion correspondantes, par rapport à chaque impulsion parmi quatre impulsions de 90 (ces dernières ayant une phase située dans un plan z-x donné),est soumise à une variation cyclique entre des demi-plans x, -x, y et -y contenant l'axe de polarisation commun z, et on effectue l'accumulation du signal en procé- dant respectivement de façon additive, additive, soustractive et à nouveau soustractive, ce qui a pour effet d'annuler les composantes d'oscillation. Selon une caractéristique d'un autre mode de réa- lisation de l'invention, on fait tourner les vecteurs de ma- gnétisation d'un angle sélectionné 0, par rapport au champ polarisant, et l'intervalle entre la rotation initiale des vecteurs de magnétisation et l'acquisition des données com- prend, dans l'ordre, une unité de retard de déphasage suivie par une impulsion de 1800, suivie à son tour par deux unités de retard pour effectuer une opération progressive de mise en phase/déphasage, suivie à nouveau par une autre impulsion de 1800 et une autre unité de retard pour le rétablissement progressif de la phase, grâce à quoi le signal est préservé sous la forme d'un double écho de spin sur l'intervalle to- tal entre la rotation initiale et l'observation, tandis que la composante d'oscillation acoustique décro t sur le mê- me intervalle, et le système est rétabli dans son orienta- tion d'origine. Selon une autre caractéristique des mesures par double écho de spin pour la compensation de l'oscillation acoustique, on utilise des mesures consécutives par double écho de spin avec des impulsionsd'inversion correspondantes dont la phase, par rapport aux impulsions d'excitation de phase constante, prend de façon cyclique les valeurs de 0 et 1800, avec une mémorisation additive des signaux respec- tifs, et prend de façon cyclique des valeurs de 900 et 2700 avec accumulation des compléments des signaux correspondants qui en résultent. Selon une caractéristique d'un autre mode de réali- sation de l'invention, on fait tourner le vecteur de magné- tisation d'un angle sélectionné G par rapport au champ pola- risant, et on déclenche l'acquisition du signal de résonance à une première phase, suivie par une période de retard cons- tante; on procède ensuite à l'inversion de la résonance par une impulsion de 1800 et on termine un retard-sélectionné en orientant à nouveau la magnétisation selon le premier axe et en procédant à l'acquisition du signal avec une phase opposée à la première phase, grâce à quoi l'oscillation résultant de chacune des deux impulsions de rotation s'annule, tandis que le signal de résonance résultant des impulsions de rotation est effectivement additifs Selon une autre caractéristique d'un autre mode de réalisation encore de l'invention, le plan dans lequel on fait tourner les résonateurs par rapport à l'axe z est sou- mis à une rotation séquentielle autour de cet axe z, de fa- çon à projeter le vecteur de magnétisation dans des plans contenant les axes x et -x (ou y et -y), gr 4 ce à quoi l'os- cillation acoustique due aux impulsions de 1800 dans les plans correspondants s'annule également sur une séquence de quatre excitations. L'invention sera mieux comprise à la lecture de la description qui va suivre de modes de réalisation et en se référant aux dessins annexés sur lesquels: La figure 1 représente une séquence d'écho de spin simple de l'art antérieur. Les figures 2 A-2 D montrent une séquence d'écho de spin à variation cyclique de la phase conforme à-l'invention. la figure 3 représente une séquence d'impulsions d'écho de spin double conforme à l'invention. Les figures 4 A et 4 B comparent des spectres obtenus par l'utilisation de la séquence de la figure 3 et un spectre de l'art antérieur. Les figures 5 A-5 D représentent une-séquence d'écho de spin double à variation cyclique de phase conforme à l'in - vention. Ia figure 6 représente un autre mode de réalisation de l'invention. Les figures 7 A et 7 B représentent une variante du mode de réalisation de la figure 6, pour la suppression de l'oscillation acoustique des impulsions de 1800. Les figures 8 A-80 comparent des spectres obtenus en utilisant le procédé des figures 7 A et 7 B, avec un procédé de l'art antérieur et un procédé à écho spin double. On va maintenant considérer la figure 1 qui montre une séquence correspondant à la mesure par écho de spin classique de l'art antérieur Une impulsion de 900, 10, fait tourner la composante de magnétisation z dans le plan x-y. Conformément à la convention habituelle, l'axe z est pris parallèle au champ magnétique polarisant La fréquence de précession du vecteur de magnétisation est répartie sur le spectre de fréquence dans un intervalle de fréquence qui est déterminé par des décalages chimiques de l'échantillon Au bout d'ame durée 12 de longueur Z, au cours de laquelle la distribution de la fréquence de précession se disperse sur un intervalle de fréquence qui est déterminé par les décala- ges chimiques, une impulsion de 180 , 14, est appliquée pour inverser lamagnéUtisaion pour l'axe z Ceci est fonctionnelle- ment équivalent à une inversion de la dispersion de phase de tous les vecteurs accomplissant un mouvement de précession. Après une durée de refocalisation 16, à nouveau de longueur V, on observe le signal 18 Pendant ce processus, le signal manifeste une atténuation exponentielle avec une constante de temps 2 t/T 2, en désignant par T 2 le temps de relaxation transverse Lorsque l'oscillation acoustique diminue avec une constante de temps Tr T 2, le signal résultant est pratiquemrrent exempt d'artefact Si au contraire Tr D T 2, le spectre sera affecté par 1 'oscillation acoustique ette dernitère est introduite par les deux impulsions 10 et 14, et elle régit l'intervalle de recoca lsation Outre la contrain- te concernant le temps de rela Tation T 2, araetéristlque des résonances observées, ce procédé exige des impulsions de 90 pour faire tourner l' ensemble de l a magnétisation dans le plan x-y La magnétisation rémanente le long de l'axe z est laissée inversée à l'achèvement de la mesure élémentaire. Les mesures par écho de spin simple correspondant à la description faite ci-dessus peuvent également etre clas- sées d'après la phase relative de l'impl-sion de retournement de 90 et de l'impulsion de refocalisation de 1800 En utili- sant la notation classique, une impulsion 90 x fait tourner la magnétisation de 900 par rapport à l'axe z, la rotation s'effectuant autour de l'axe x Une rotation autour de l'axe -x a le sens opposé, c'està-dire qu'elle se produit dans le demi-plan opposé La séquence de MeiboomGill est représentée par les impulsions 900 x -T -1800 acquisition, tandis que la séquence de Carr-Purcell est représentée par le pro- gramme d'impulsions 90 x Z 180 c acquisition. x x L'impulsion de 90 comme l'impulsion de 180 d'une mesure par écho de spin simple sont des sources d'oscillation acoustique. Les figures 2 A-2 D montrent un exemple d'un procédé préféré pour compenser l'oscillation dans une mesure du type à écho de spin simple La sous-séquence de la figure 2 A représente une mesure de Carr-Purcell caractéristique Celle-ci est sui- vie par une seconde sous-séquence (figure 2 B), à nouveau de la forme Carr-Purcell, qui se distingue par la relation de phase opposée pour l'impulsion de 180 O 14 x par rapport à l'impulsion de 1800 correspondante 14 de la figure 2 A Des x C signaux cumulatifs dans le domaine des temps pour ces deux sous-séquences additionnent l'oscillation due aux impulsions de 900, 1 Ox, tandis qu'ils annulent l'oscillation due aux impulsions de 180 , 14-x et 14 x Les sous-séquences des figu- res 2 C et 2 D, correspondant à la séquence de Meiboom-Gill, accomplissent la même annulation entre les impulsions de 180 respectives 14 y et 14 _y, tandis que l'oscillation provenant des impulsions respectives de 90 , 10 x, s'ajoute à nouveau. Du fait que les acquisitions sont maintfnant effectivement soustraites du signal cumulatif pour les sous-séquences des figures 2 C et 2 D, les composantes d'oscillation à 90 sommées des figures 2 A et 2 B arnnulenlt conjointement les composantes d'oscillation à 90 sommées des sous-séquences des figures 20 et 2 Do La soustraction effective peut être réalisée en modi- fiant la phase du récepteur, comme représenté, ou par sous- traction du signal converti sous forme numérique qui est acquis au moment considéré, par rapport au signal cumulatif mémorisé, en employant des moyens numériques simples Le si- gnal, dont la phase a été inversée pour les deux dernières sous-séquences,est maintenant additif et toute l'oscillation a été annulée Du coté gauche des figures 2 A-2 D, le diagramme vectoriel montre la rotation 900 x qui fait tourner la magnéti- sation de l'axe z vers le plan x-y, par l'application d'une impulsion 10 dirigée le long de l'axe +x Dans chacune des quatre sous-séquences, la magnétisation transverse Mt demeure animée d'un mouvement de précession, puis elle est soumise à l'action d'une impulsion de 180 dirigée le long desaxes +x, -x, +y, -y, ce qui produit les rotations autour des axes correspondants, comme il est représenté. Les techniques d'écho de spin, décrites ci-dessus, utilisent des impulsions de 900, alors qu'il est bien connu qu'on peut optimiser la sensibilité expérimentale en utili- sant des angles de basculement plus faibles, à une cadence de répétition plus élevée On peut réaliser ceci d'après le procédé de la figure 3, correspondant à un autre mode de réalisation de l'invention, en utilisant des impulsions 60, pour produire une projection dans le plan x-y du vecteur de magnétisation, de valeur suffisante pour satisfaire un critère de rapport signal à bruit désiré Les composantes animées d'un mouvement de précession dans le plan x-y se dispersent pendant l'intervalle 62 et l'inversion pour l'axe z est accomplie par l'impulsion de 1800 f, 64 Les composan- tes de vecteurs qui sont animées d'un mouvement de précession dans le plan x-y se refocalisent pendant l'intervalle 66 A. Si l'acquisition devait commencer à ce point (écho de spin simple), l'expérience serait confrontée à une composante de magnétisation -z élevée (inversée) Par conséquent, on permet aux vecteurs animés d'un mouvement de précession de se dis- perser à nouveau sur l'intervalle 66 B, puis une impulsion de réinversion 68 rétablit la composante de magnétisation z à son sens d'origine, ce qui inverse une fois de plus la phase des composantes animées d'un mouvement de précession, dans le plan x-y La refocalisation se produit à nouveau sur l'in- tervalle 70, après quoi on déclenche l'acquisition du signal 72. On va maintenant considérer les figures 4 A-et 4 B qui comparent un spectre de l'art antérieur, obtenu par re- tard et acquisition, et un spectre obtenu par écho de spin double Le pic est dû à 15 N, obtenu à partir d'une solution à 90 % de CHONH 2, en utilisant un appareil Varian XL-200 fonctionnant à 20,3 M Hz, avec une largeur spectrale de 10 k Hz Les spectres représentés sur la figure 4 A correspondent à l'interposition d'un retard de 20 ps entre l'impulsion et l'acquisition des données, tandis que le spectre à écho de spin double de la figure 4 B est obtenu avec un temps caracté- ristique t = 500 >s, pour un intervalle total de 2 ms en- tre l'excitation et l'acquisition Les spectres sont tracés sur la même échelle, après un calcul de moyenne sur 512 ex- citations Les impulsions d'observation sont à 80 et un élar- gissement de la raie à 20 Hz est appliqué pour faire ressor- tir la variation de la ligne de base Dans cette comparaison, le temps T 2 pour l'échantillon étudié est très supérieur à la durée du signal d'oscillation acoustique. Le procédé à écho de spin double qui est décrit de- meure soumis à des contraintes qui consistent en ce que l'os- cillation qui résulte des impulsions de 180 s'atténue en un temps Tr extension évidente de celle qu'on a envisagée pour les figu- res 2 A-2 D Comme dans ce mode de réalisation, les diverses impulsions de 1800 correspondantes qui interviennent dans les sous-séquences sont également soumises ici à une variation cyclique de la phase sous la forme +x, -x +y, -y (ou en or- dre inverse) Les impulsions de réinversion 68 portent les indices + p + q pour faire ressortir la variation cyclique de phase qui est appliquée séparément à ces impulsions la seule exigence pour la relation entre ces phases est: Ip-q I = lr /2 On peut maintenant donner aux intervalles d'évolution/ refocalisation des durées extrêmement courtes, sans contrain- tes résultant des longueurs relatives du temps de-décroissan- ce de l'oscillation acoustique et des temps de relaxation transverse qui sont caractéristiques des résonances observées. On va maintenant considérer la figure 6 qui montre que la séquence d'impulsions du mode de réalisation préféré utilise une impulsion 8 x, 80, et une acquisition de signal immédiate, 82 Le signal est accumulé en mémoire de façon additive, ce qui est indiqué par le symbole "+" Une durée de retard fixe 84 permet au signal de récupérer un niveau acceptable après l'application d'une impulsion de 180 , 86, après quoi la magnétisation est inversée pour l'axe z et un autre intervalle 88 de longueur W permet la décroissance de l'oscillation qui est due à l'impulsion à 180 , 86 Une im- pulsion 8 x, 90, identique à l'impulsion 80, est appliquée et une autre acquisition avec signe, 92, est déclenchée et est soustraite du signal cumulatif en mémoire, ce qui est désigné symboliquement par 1 "-" et correspond par exemple au signe op- posé à celui de la première paire d'acquisitions Le résultat conjoint de l'inversion de 180 et des phases d'accumulation opposées consiste en ce que les signaux acquis s'additionnent, tandis que la partie d'oscillation qui résulte des signaux Gx, bien que cohérente, n'est pas modifiée par l'inversion spa- tiale, ce qui conduit à l'annulation de la contribution d'os- cillation acoustique entre des paires d'excitations G x Conformément à une autre forme du mode de réalisa- tion décrit ci-dessus, les figures 7 A et 7 B illustrent l'ap- plication d'lu cycle de phase à la séquence de la figure 6, dans le but de compenser la contribution d'oscillation des impulsions d'inversion 86 la séquence de la figure 6 est sui- vie mqr une séquence similaire dans laquelle les impulsions de 180 sont soumises à une variation cyelique alterrnative enltre +x et -x (ou +y et y) Il en résulte que la contrainte concernant R et Tr (de l'impl Jsion de 1800 prècédente) est r assouplie Par exemple, Z peut etre assez court et/ou '1 peut Ctre assez long, sans effet défavorable. Dans la forme décrite ci-dessus, il faut évidemment prendre garde au fait que la magnétisation z est inversée à l'achèvement d'une séquence ne comportant que deux excitations. Si le retard 84 est suffisamment long par rapport à T 1, les deux excitations d'une séquence élémentaire sont pratiquement indépendantes, dans la mesure o la magnétisation z inversée restante retourne à sa valeur d'équilibre avant l'excitation suivante On peut assouplir cette contrainte s'il existe une impulsion de 180 supplémentaire dans une séquence quelconque, pour rétablir la magnétisation inversée Ceci peut prendre la forme de l'insertion d'une telle impulsion d'inversion supplémentaire 94 à l'achèvement de l'acquisition 92 Il est préférable que l'impulsion de réinversion 94 soit soumise à une variation cyclique de phase L'insertion d'impulsions d'inversion 94 avec variation cyclique de phase garantit donc l'annulation des effets de magnétisation rémanente et de l'os- cillation acoustique qui résulte de telles impulsions L'os- cillation qui est introduite par les impulsions de réiriversion 94 peut être compensée sans affecter le signal, si les phases relatives d'impulsions de réinversion alternées sont égales à q T /2 La phase relative de ces impulsions est arbitraire lorsque T 1 ou T 2 est relativement court et, bien entendu, lors- que les deux temps T 1 et T 2 sont relativement courts, les im- pulsions de réinversion sont inutileso Les figures 8 A-80 montrent des exemples de spectres obtenus par le procédé de retard de l'art antérieur et par les séquences des figures 3 et 6 On a examiné une solution contenant 25 % de D 20 et 75 % d'acétone sur un spectromètre Varian XL-200, à une fréquence de 27,1 M Hz La largeur spec- trale était de 30,030 k Hz et les spectres résultants repré- sentent la moyenne sur 1024 excitations, à une cadence de ré- pétition de 20 ms Les deux pics sont dus aux éléments 170 de l'eau lourde et de l'acétone la figure 8 C montre lmun exem- ple de spectres obtenus avec le mode de réalisation de la fi- gure 6 La figure 8 A montre l'exem-le concernant le procédé d'acquisition retardée de l'art nvtérieur, pour un certain nombre de durées de retard intercalées entre l'impulsion d'observation et l'acquisition Des artefacts spectraux nota- bles sont présents, sauf pour la plus longue durée de retard. Les six spectres tracés avec la même échelle sont paramétrés par valeurs successives de temps de retard ou de temps de re- focalisation La différence entre les spectres pour les va- leurs respectives de V est minimale, tandis que la différen- ce avec les spectres des figures 8 A et 8 B se manifeste clai- rement par l'absence d'artefacts spectraux pour les valeurs les plus courtes de D ou du temps de retard Contrairement à l'exemple de la figure 4 B dans lequel T 2 >> Tr, la mesure par écho de spin double de la figure 8 B est un exemple dans lequel T 2 n'est plus aussi grand par rapport à Tro Sur les figures 8 B et 8 C, l'impulsion d'observation et la phase du 1 1 récepteur ont été soumises à une variation cyclique de 900, en plus de la variation de phase cyclique qui fait l'objet de l'invention Conformément à la technique connue, les im- pulsions de 1800 sont des impulsions composites, pour minimi- ser les effets de défauts d'étalonnage et de couplages hors résonance. L'homme de l'art notera que bien qu'on ait envisagé un effet parasite particulier, les procédés décrits ici peu- vent compenser tout transitoire similaire d'origine instru- mentale. Il va de soi que de nombreuses autres modifications peuvent être apportées auxprocédés décrits et représentés, sans sortir du cadre de l'invention. REVENDICATIONS 1 Procédé de spectroscopie par résonance gyromagné- tique agissant sur un ensemble de résonateurs gyromagnétiques immergés dans un champ magnétique polarisant, ce champ étant dirigé selon un premier axe (z), dans lequel les vecteurs de magnétisation des résonateurs gyromagnétiques accomplissent un mouvement de précession autour du champ magnétique polari- sant, caractérisé en ce qu'il comprend les opérations suivan- tes: (a) on fait tourner les vecteurs magnétiques de 900 à partir du premier axe, par une rotation du vecteur magnétique autour d'un second axe (x), ce second axe étant orthogonal au premier axe et la rotation s'effectuant dans un sens sélection- né, grace à quoi les vecteurs magnétiques tournent dans un plan orthogonal au premier axe; (b) on permet la perte pro- gressive de la cohérence de phase entre ces vecteurs magnéti- ques tournants, pendant un intervalle de temps T; (c) on inverse le sens de rotation des vecteurs magnétiques en in- versant le sens du premier axe, par une rotation de 1800 par rapport au premier axe, grace à quoi la perte progressive de cohérence de phase est inversée; (d) on permet aux vecteurs magnétiques tournants de s'approcher progressivement de la cohérence de phase; (e) on détecte un signal de résonance d'écho de spin des résonateurs, ce signal pouvant contenir des transitoires d'origine instrumentale qui-résultent des opérations (a) et (c); (f) on mémorise le signal de résonan- ce; (g) on répète les opérations (a) et (b) et on inverse le sens du premier axe par rotation de 1800 du système de coordoÈn- nées, 'par rapport au premier axe, cette rotation s'effec- tuant autour du second axe, avec un sens de rotation opposé au sens sélectionné, et on répète les opérations (d) et (e); (h) on additionne le signal mémorisé à l'opération (f) avec le signal détecté à l'opération (g) pour former un signal cumulé exempt de transitoires résultant des rotations de 1800 des opérations (c) et de la réinversion de l'opération (g); (i) on répète les opérations (a) et (b) et on inverse le sens du premier axe par rotation du vecteur magnétique autour d'un troisième axe (y), cet axe étant orthogonal au premier axe (z) et au second axe (x), et cette rotation de 1800 ayant un second sens sélectionné, gr Cce à quoi la perte de cohé- rence de phase est inversée, et on répète à nouveau les opé- rations (d) et (e); (j) on soustrait le signal détecté à l'opération (i) du signal cumulé obtenu à l'opération (h) et on mémorise le signal cumulé; (k) on répète les opérations (a) et (b) et on inverse le sens du premier axe par rotation du vecteur magnétique sur 1800 autour du troisième axe, cet- te rotation de 180 s'effectuant en sens opposé au second sens sélectionné, et on répète les opérations (d) et (e); et ( 1) on soustrait le signal détecté à l'opération (k) du signal cumulé obtenu et mémorisé à l'opération (j) et on mé- morise le nouveau signal cumulé, grace à quoi l'oscillation transitoire qui résulte de l'opération (a), répétée à l'opé- ration (g), est annulée par l'oscillation transitoire qui se manifeste dans les opérations (i) et (k). 2 Procédé de spectroscopie par résonance gyroma- gnétique agissant sur un ensemble de résonateurs gyromagné- tiques immergés dans un champ magnétique polarisant orienté selon un premier axe (z), dans lequel les vecteurs de magné- tisation des résonateurs gyromagnétiques accomplissent un meuvement de précession autour du champ ma nétique polarisant, caractérisé en ce que: (a) on fait tolurner les vecteurs de ziagnétisation d'uni angle S 1 lectionné par rapport a-u premier axe (z), cette rotation s'effectuant dans un sens sélectionné autour d'un second axe (x) orthogonal au premier axe; (b) on permet la perte progressive de la coherence de phase entre les résonateurs gyromagnétiques pendant un intervalle de temps T; (c) on inverse le sens du premier axe, par rotation du système sur 180 par rapport au premier axe, cette rota- tion s'accompagnant d'une oscillation transitoire; (d) on permet l'augmentation de la cohérence de phase parmi les ré- sonateurs gyromagnétiques pendant un autre intervalle de temps; (e) on répète l'opération (c); et (f) on détecte la résonance d'écho de spin double des résonateurs, grâce à quoi la composante de magnétisation qui demeure orientée selon le premier axe après l'opération (a) et qui est inversée par l'opération (b) est rétablie dans son orientation initiale. 3 Procédé selon la revendication 2, caractérisé en ce que la rotation de l'opération (c) est soumise à des variations de phase cycliques, entre des phases séparées de , au cours de répétitions suivantes des opérations (a) à (f) 4 Procédé de spectroscopie par résonance magnéti- que de résonateurs gyromagnétiques, caractérisé en ce que: (a) on polarise les vecteurs magnétiques des résonateurs se- lon l'axe de polarisation (z), ce qui fait que les vecteurs représentatifs de ces résonateurs sont soumis à un mouvement de précession autour de l'axe de polarisation; (b) on fait tourner les vecteurs magnétiques de façon cohérented'un angle désiré par rapport à l'axe de polarisation; (c) on détecte et on mémorise un premier signal de résonance résul- tant ( 82); (d) on observe un retard correspondant à un in- tervalle de temps fixe; (e) on inverse les vecteurs de ma- gnétisation de 180 par rapport à la direction de polarisa- tion; (f) on fait à nouveau tourner les vecteurs magnétiques de façon cohérente de l'angle désiré par rapport à l'axe de polarisation, et on obtient un second signal de résonance ( 92); et (g) on soustrait le second signal du premier pour former un signal cumulé et on mémorise ce signal cumulé, grace à quoi les parties non résonsntes du signal qui accompagnent les opérations consistant à faire tourner les vecteurs de fa- çon cohérente et à faire tourner a nouveau les vecteurs de faon cohérente s'annulent, tandis que les pares résonnrmantes des signaux s'ajoutent. 5 Procédé selon la revendication 4, caractérisé en ce qu'on répète les opérations (a) à (g) incluses et on ef- fectue des opérations d'inversion alternées avec un écart de phase de 180 0. 6 Procédé selon la revendication 5, caractérisé en ce que l'opération de soustraction et de mémorisation comprend en outre une réinversion des vecteurs magnétiques soumis à un mouvement de précession. 7 Procédé selon la revendication 6, caractérisé en ce que des opérations de réinversion alternées présentent un écart de phase de 90 .