FR 2547062 A1 19841207 FR 8308914 A 19830530 La présente invention se rapporte à la géophysique et a notamment pour objet un procédé de détection des mouvements géodynamiques actuels dans un massif montagneux L'invention peut être utilisée avec le plus de succès dans la reconnais5 sance du sol de fondation Elle peut trouver son application là ou l'on se heurte à des problèmes scientifiques et pratiques que l'on peut résoudre à l'aide des méthodes de levé et d'établissement de cartes structurales et géodynamiques: dans le génie minier, la seismologie, la recherche des gisements 10 de minéraux utiles. L'un des facteurs essentiels réduisant le risque géologique lors de l'élaboration des projets et la construction des ouvrages de génie civil modernes (surtout des ouvrages uniques en leur genre), est la présence de cartes structurales 15 et géodynamiques détaillées qui fournissent des informations sur la structure tectonique de la région à étudier, ainsi que sur le caractère et l'intensité des processus géodynamiques actuels Ces informations sont surtout nécessaires dans le cas de régions de séismicité élevée, de zones caractérisées par 20 un développement karstique, de zones de glissements, de structures tectoniques actives ainsi que dans les régions de forte influence de travaux de génie civil sur le milieu montagneux. Ce sont les données obtenues lors des travaux de forage 25 qui sont à la base des cartes structurales Mais, il faut le reconnaître, ces travaux exigent des dépenses considérables et, de par leur nature, ne peuvent pas être directement utilisés pour obtenir des informations sur le caractère des mouvements géologiques actuels. Les méthodes géophysiques classiques utilisées pour l'établissement des cartes structurales sont basées sur l'analyse des caractéristiques des champs d'origine naturelle ou artificielle Selon la présence et le caractère des perturbations apportées dans ces champs, on peut juger des modifi- 2 2547062 cations tectoniques et de leur disposition, ce qui permet d'obtenir des données pour l'établissement des cartes structurales En principe, tous les champs géophysiques renferment des informations sur le caractère des mouvements géodynamiques actuels Toutefois, il est pratiquement impossible d'obtenir ces informations sous une forme nette, du fait des perturbations apportées dans les champs à étudier par l'hétérogénéité du milieu montagneux (cette hétérogénéité rend plus compliqués les autres travaux utilisant les métho10 des géophysiques: établissement de cartes structurales, recherche de gisements de minéraux utiles) Par exemple, les procédés d'établissement de cartes structurales utilisent des champs électromagnétiques induits et qui sont largement employés (voir, par exemple, Ao B Velikin, G S Frantov 15 "Elektromagnitnye polya, primenyaemye v induktsionnykh metodakh elektromagnitnoi razvedki", Leningrad, 1962), nécessitent souvent l'utilisation d'un sondage profond du milieu montagneux pour que le champ électromagnétique puisse atteindre l'objet de la recherche(la modification tectonique)masqué 20 par les dépôts de couverture, l'intensité de ce champ étant suffisante pour enregistrer la perturbation due à cet objet. Dans le cas d'un tel sondage 4 toute l'assise des roches de couverture se trouvant au-dessus de l'objet de la recherche 25 joue le rôle d'un écran relativement à cet objet, et on voit se manifester dans le champ électromagnétique enregistré l'hétérogénéité de ladite assise, les perturbations de ce champ dues à ladite hétérogénéité étant parfois identiques (de par leurs intensité et leur morphologie) à celles dues 30 à l'objet de la recherche lui-même. On assiste à un phénomène essentiellement identique en cas de sondage circulaire pour l'établissement de cartes et la recherche de minéraux utiles C'est un type de sondage profond o le champ électromagnétique induit dans le massif 35 montagneux est étudié dans des directions azimutales différentes à partir d'une source de rayonnement électromagné- 3 -2547062 tique à une distance de dizaines de mètres (voir V Ya. Khmelevsky, "Osnovnoi kurs elektrorazvedki", vol II, éditions de l'Université de Moscou, 1971, p 256). Dans le cas d'un tel rayon de la zone à étudier, on se trouve en présence de roches non homogènes et non isotropes caractériséespar une période de relaxation prolongée (plusieurs décennies), ce qui ne permet pas d'utiliser ledit procédé pour la détection des mouvements géodynamiques actuels Les études en profondeur se caracté10 risent par la présence probable d'inclusions magnétiques dans le milieu à étudier, ce qui nécessite une technique de mesure plus compliquée. Le fait que l'assise des roches de couverture joue le rôle d'un écran relativement à l'objet de la recherche, 15 analogue à celui décrit cidessus, se fait également sentir lors de l'étude des champs de force naturels (gravitationnel, magnétique, gazeux, etc) à l'aide de méthodes de prospection géophysiques. Les informations concernant le caractère des mouvements 20 géophysiques actuels peuvent être obtenues grâce aux méthodes instrumentales classiques (nivellements de précision de retour, inclinométrie, déformographie, extensométrie, etc) Toutefois, toutes ces méthodes présentent un caractère local et nécessitent des mesures compliquées et onéreuses visant à neutrali25 ser l'influence de facteurs de déséquilibre tels que la température, l'humidité, la pression,etc;sur les résultats des mesures, ce qui affecte l'efficacité de ces méthodes. Il y a lieu de noter que l'ensemble des méthodes classiques énumérées peuvent être utilisées, en principe, pour 30 l'établissement de cartes structurales et géodynamiques, mais lorsqu'on se trouve devant des problèmes pratiques pour la solution desquels on doit disposer d'informations détaillées concernant un territoire relativement étendu, les dépenses nécessaires à l'établissement de cartes de ce genre sont excessives. La première tentative de trouver une solution valable au: point de vue économique a abouti à la mise au point d'un procédé de détection des mouvements géodyna5 miques actuels dans les structures tectoniques à l'aide d'un levé d'émanation (voir le certificat d'auteur URSS N 396659). Ledit procédé est basé sur l'étude de la capacité d'émanation des dépôts de couverture et consiste à prélever de l'air tellurique, à l'aide d'un échantillonneur de gaz et 10 d'une pompe à vide, à partir de puits dont la profondeur va jusqu'à 0,8 m et qui sont préparés sur les piquets des profils, et à déterminer, à l'aide d'analyseurs-, la concentration et la composition des gaz radio-actifs (radon et thoron) Selon la variation de la concentration desdits gaz par rapport à celle 15 de fond, on peut juger des mouvements géodynamiques actuels dans les structures tectoniques. La possibilité d'un tel emploi du procédé connu de levé d'émanation de la zone adjacente à la surface est basée sur le fait, établi expérimentalement, qu'il existe une dépendance, 20 dans l'espace et dans le temps, entre les anomalies d'émanation et les zones des mouvements géodynamiques actuels. Le procédé décrit ne demande pas un appareillage onéreux et consommant beaucoup d'énergie et diffère avantageusement des procédés connus par sa rapidité, ce qui permet, notamment, 25 de l'utiliser en vue de l'étude des processus dynamiques qui se déroulent relativement vite et/ou sont finement différenciés dans l'espace, processus inaccessibles jusqu'à présent à l'étude. Toutefois, la nécessité de pénétrer dans le milieu à étudier (enfoncement de l'échantillonneur de gaz) conduit à un accroissement de la quantité de travail (d'habitude, des mesures sont effectuées par une équipe de deux ou trois personnes), à une augmentation du coût des travaux, à un endommagement des revêtements artificiels (routes, pavés, trottoirs, etc), et, lorsque les travaux sont effectués dans une ville, on doit prévoir des mesures spéciales visant à éviter tout risque de détérioration des canalisations (câbles, conduits d'eau et de gaz, etc) En plus de cela, les travaux de terrassement rendent ledit procédé moins rapide. Etant donné que ce procédé est basé sur l'étude du champ naturel, on ne peut pas modifier arbitrairement les propriétés qui sont inhérentes à ce champ En particulier, en ce qui concerne le radon, qui est le plus informatif des deux gaz indiqués, le domaine d'étude du milieu est constitué par un demi-espace inférieur d'un rayon de 10 m environ, et il est impossible de régler la profondeur d'étude à une autre valeur correspondant à la profondeur de gisement du milieu à étudier Ceci conduit à une altération des résultats des 15 mesures, due à la superposition des variations du champ d'émanation causées par l'hétérogénéité du milieu à étudier. La certitude des informations obtenues est sensiblement diminuée du fait de l'influence des régimes barohydrothermiques des sols, de leur type et de leur état initial (par 20 exemple, sol rendu meuble) sur les résultats du levé d'émanation. L'inconvénient essentiel de ce procédé réside dans le caractère peu informatif de la caractéristique physique du milieu qui est choisie en tant qu'objet des mesures et dont 25 la variation permet d'apprécier la variabilité des champs de tensions tectoniques actuels et les mouvements géologiques dans le massif montagneux, mouvements correspondant à cette variabilité Les essais de laboratoire effectués à l'Institut de Prospection Géologique (Moscou) ont montré que la capa30 cité d'émanation des sols est déterminée non seulement par l'intensité des mouvements actuels, mais aussi par leur dynamique Il est établi que des mouvements même insignifiants; mais susceptibles de varier rapidement, augmentent fortement la capacité d'émanation des sols, tandis que les mouvements 35 intenses, mais susceptibles de varier lentement, peuvent conduire, dans des conditions déterminées,à un affaiblis- sement de cette capacité De ce fait, la vitesse de développement des variations de la capacité d'émanation des sols, déterminée lors des mesures de retour, permet de déterminer l'une des caractéristiques quantitatives des mouvements actuels: le degré de leur variabilité (dynamique) dans le temps pour un terrain déterminé (car ce degré de variabilité se fait sentir le plus nettement au-dessus des zones de perturbations techniques et des zones soumises à leur influence), mais ne permet pas de 10 déterminer quantitativement l'intensité des mouvements détectés. Les informations obtenues lors des mesures décrites ci-dessus, portant sur la capacité d'émanation des sols, présentent un caractère scalaire, le caractère des varia15 tions de cette capacité, détectées lors des observations de retour, étant lui aussi scalaire De ce fait, ledit procédé ne permet pas de déterminer le caractère de la déformation déterminée par les processus géodynamiques détectés, ce qui, compte tenu de ce qui précède, affecte 20 le caractère informatif du procédé. En outre, le domaine d'application d'un tel procédé est limité du fait que les mesures de l'émanation des gaz du sol, décrites ci-dessus, ne peuvent être effectuées qu'en présence de roches de couverture telles que: argiles, 25 sols argileux, terres noires. Le but de l'invention est d'éliminer les inconvénients mentionnés. L'invention vise par conséquent un procédé de détection des mouvements géodynamiques actuels dans un massif 30 montagneux, assurant une amélioration de la capacité informative, tout en permettant de réduire les frais pour sa mise en oeuvre grâce au choix, en tant qu'objet de mesures, d'une caractéristique physique du milieu à étudier telle qu'elle soit douée d'une capacité informative élevée quant I au champ de tensions tectoniques et qu'elle soit susceptible d'un contrôle rapide sans contact, ainsi que grâce à une modification de la technique des mesures telle qu'elle corresponde de manière optimale à la nature de l'objet de mesures choisi. Le problème ainsi posé est résolu à l'aide d'un procédé de détection des mouvements géodynamiques dans un massif montagneux, consistant à enregistrer et à apprécier les variations des paramètres d'un milieu adjacent à la surface, 10 ledit procédé étant caractérisé, selon l'invention, en ce que l'on engendre préalablement dans ledit milieu un champ électromagnétique à l'aide d'une source de rayonnement électromagnétique, on enregistre les variations des paramètres du milieu exposé au rayonnement en comparant les résultats des mesures périodiques des séries de valeurs de la fonction d'au moins l'un des paramètres électromagnétiques de ce milieu en des points disposés suivant une circonférence de rayon allant jusqu'à 5 m par rapport à un centre situé en un point d'un profil, et on apprécie les variations des paramètres du milieu exposé au rayonnement d'après les variations de la caractéristique qualitative de l'anisotropie azimutale du paramètre électromagnétique à étudier du milieu exposé au rayonnement, caractéristique représentée par une série de valeurs de la fonction dudit paramètre enregistrées lors 25 desdites mesures périodiques, la présence de mouvements géodynamiques actuels provoquant une déformation du type compression-traction des roches étant constatée en présence de ladite caractéristique qualitative stable, tandis que la présence de mouvements géodynamiques actuels conduisant à une déformation du type cisaillement des roches est constatée en présence de ladite caractéristique qualitative variable. L'objet des mesures, dans le procédé revendiqué, est constitué par une caractéristique physique du milieu qui diffère fondamentalement de celle qui est choisie dans le procédé connu Cette caractéristique physique est l'anisotropie azimutale des propriétés électromagnétiques, qui est observée dans un milieu soumis à l'action d'un champ magnétique et faisant partie du massif montagneux ou entrant en contact, même temporaire, avec celui-ci Les auteurs de la présente invention ont établi expérimentalement (en traitant les résultats de plus de 300 000 mesures) que cette anisotropie se fait sentir à tout niveau contrôlé, y compris 10 la zone qui se trouve à proximité de la surface du sol (tant en profondeur que du côté extérieur), et qu'elle correspond au champ de tensions tectonique actuel et aux mouvements géodynamiques actuels qui sont liés à ce champ. Il convient de noter que l'anisotropie azimutale des 15 propriétés électromagnétiques des milieux montagneux et des revêtements artificiels, dont la présence a été mise en évidence par les auteurs de l'invention, présente la même nature géodynamique et n'est pas identique à l'anisotropie due à l'hétérogénéité naturelle des milieux montagneux. 20 Cette dernière anisotropie est déterminée à l'aide du procédé connu, décrit ci-dessus, de sondage circulaire du milieu montagneux, dans lequel la méthode de mesures et l'appareillage mis en oeuvre ne permettent en aucun cas de détecter l'anisotropie de nature géodynamique et, par consé25 quent, les mouvements géodynamiques actuels, car aux distances pour lesquelles ledit procédé est applicable, ces informations sont neutralisées par les informations sur l'hétérogénéité naturelle du milieu. Lorsqu'on procède aux mesures, suivant le procédé conforme à l'invention, en des points situés sur une circonférence de rayon-allant jusqu'à 5 m, ce sont les informations relatives à l'anisotropie due aux champs de tensions tectoniques actuels qui prédominent Ceci s'explique par le fait que, dans les limites de surfaces et de profondeurs peu importantes, le milieu à étudier est généralement amagnétique, il présente une homogénéité et une isotropie naturelles des propriétés électromagnétiques, ainsi qu'une vitesse de relaxation élevée (plusieurs heures), ce qui- offre la possibilité de contrôler la dynamique des variations de ses propriétés physiques et, notamment, des propriétés électromagnétiques, variations dues à celles du champ de tensions tectonique actuel. Du fait que le champ engendré dans le milieu à étudier o 10 est artificiel, on a la possibilité de régler la profondeur de recherche en modifiant les paramètres des moyens techniques mis en jeu pour la réalisation du procédé Ceci permet de limiter la profondeur d'étude du milieu à des valeurs pour lesquelle le milieu est de nature homogène et isotropeet d'augmenter 15 par conséquent la fiabilité des résultats des mesures effectuées conformément au procédé revendiqué. Le fait d'éliminer la nécessité d'une pénétration matérielle dans le milieu à étudier permet de diminuer sensiblement les frais pour la réalisation des travaux nécessaires à 20 la détection des mouvements géodynamiques actuels grâce à la diminution de la quantité de travail (au lieu des 2 ou 3 personnes nécessaires à la réalisation des mesures suivant le procédé connu, on a besoin d'une seule personne dans le cas du procédé revendiqué), ainsi que grâce à l'élimination 25 du risque d'endommagement lors des travaux de terrassement, mentionné ci-dessus. La nature de la caractéristique physique choisie en tant quaobjet des mesures suivant le procédé revendiqué admet l-emploi de méthodes d'analyse non paramétriques, ce qui élimine pratiquement les facteurs de déséquilibre intervenant lors de l'appréciation des résultats des mesures et augmente sensiblement la fiabilité des conclusions. Le choix de l'anisotropie azimutale des propriétés électromagnétiques du milieu à étudier en tant qu'objet des mesures, ainsi que la technique de mesure de la fonction des paramètres électromagnétiques dans les différentes directions, permettent d'apprécier les mouvements géodynamiques actuels en grandeurs vectorielles et de déter5 miner le caractère des déformations dues à ces mouvements. Une telle appréciation qualitative des déformations constitue l'un des avantages du procédé revendiqué et, comme on le verra ci-dessous, offre la possibilité d'améliorer encore plus la capacité informative du procédé. Pour l'appréciation quantitative des déformations dues aux mouvements géodynamiques actuels, le milieu à exposer au rayonnement peut avantageusement être porté à un état quasi-isotrope avant chacune des mesures successives, si les valeurs de la fonction du paramètre électromagnétique à étudier, 15 obtenues lors des mesures précédentes, sont sensiblement différentes. Le fait de porter le milieu à l'état quasi-isotrope avant les mesures assure, dans tous les cas, la possibilité d'apprécier la vitesse de développement de l'anisotropie a20 zimutale du paramètre à contrôler dans le milieu à exposer au rayonnnement et, par çóla-même, d'apprécier l'intensité des mouvements géodynamiques actuels, qui est fonction de ladite vitesse. Pour l'appréciation quantitative des déformations du 25 type cisaillement, il est avantageux, dans chaque quadrant du cercle formé par les points en lesquels on effectue les mesures périodiques des séries de valeurs de la fonction du paramètre électromagnétique du milieu exposé au rayonnement, de choisir une valeur maximale de ladite fonction pour la 30 série considérée, de comparer les positions azimutales des points auxquels correspondent, dans le même quadrant, les valeurs maximales choisies, obtenues lors des mesures espacées dans le temps, et, d'après le rapport entre la différence des azimuts des points à comparer et l'intervalle de temps entre les mesures à comparer, d'apprécier l'intensité des mouvements qui provoquent ladite déformation des roches. Pour l'appréciation quantitative des déformations du 5 type compression- traction, il est avantageux, dans chaque quadrant du cercle formé par les points en lesquels on effectue les mesures périodiques des séries de valeurs de la fonction du paramètre électromagnétique du milieu exposé au rayonnement, de choisir les valeurs maximale et minimale 10 de ladite fonction pour la série considérée, de comparer le rapport de ces valeurs obtenues lors des mesures espacées dans le temps, et, d'après la vitesse de variation desdits rapports, d'apprécier l'intensité des mouvements géodynamiques actuels qui provoquent ladite déformation des roches. Une telle technique de mesures permet d'éliminer efficacement les facteurs aléatoires qui influencent l'appréciation des résultats des mesures. En tant que milieu exposé au rayonnement on peut utiliser la couche superficielle du massif montagneux L'avantage 20 d'un tel milieu réside en ce qu'il présente un caractère informatif et qu'il est facilement accessible A le différence du procédé connu, le procédé revendiqué impose des impératifs moins sévères à l'indice d'homogénéité et d'isotropie qui doit caractériser cette couche, par suite du caractère non 25 paramétrique du traitement des résultats des observations prévues par le procédé revendiqué, ce qui découle de ce qui précède Pour la même raison, dans le procédé envisage, on peut négliger le facteur barohydrothermique Néanmoins, l'utilisation de la couche superficielle du massif montagneux 30 en tant que milieu exposé au rayonnement n'est pas toujours désirable et possible, elle ne l'est que dans les cas o cette couche, de par sa nature, est homogène et isotrope et possède une faible mémoire, c'est-àdire si elle présente une période de relaxation peu considérable (dans les limites de 20 à 30 heures, mais il est désirable que cette période soit bien moindre) d'après la caractéristique physique étudiée Ces propriétés sont typiques notamment pour les argiles, les sols argileux, les terres noires, la glace dont sont recouverts les bassins d'eau, ainsi que pour certains revêtements artificiels tels que l'asphalte. La couche superficielle du massif montagneux, qui est constituée par les milieux rocheux, incohérents ou aqueux, est peu informative pour le procédé revendiqué Dans ce 10 cas, il est avantageux d'utiliser, en tant que milieu exposé au rayonnement, un milieu de contrôle mis en contact avec la couche superficielle du massif montagneux et constitué par une matière amagnétique présentant l'homogénéité et l'isotropie du paramètre électromagnétique à étudier et dont le point de transformation de la phase solide en phase liquide se trouve dans un intervalle de changement naturel des températures. Mis en phase solide, de tels milieux sont de bons porteurs d'information, tandis que la mise en phase liquide 20 de ceux-ci assure un effacement fiable de l'information accumulée En particulier, en qualité de tels milieux on peut utiliser la glace, les sols argileux humides ou l'asphalte L'avantage de la glace réside en ce qu'elle est presque toujours disponible ou qu'elle peut être obtenue en 25 présence de températures correspondantes; l'asphalte, quant à lui, peut être également obtenu d'une manière relativement simple à partir de matériaux largement répandus; enfin, les sols argileux constituent un milieu qui est classé parmi les milieux le plus répandus. L'utilisation des milieux de contrôle élargit amplement les possibilités technologiques du procédé en rendant possible son emploi dans tous les cas, indépendamment du type des sols qui constituent la couche superficielle du massif montagneux. 1-3 L'anisotropie des propriétés électromagnétiques du milieu exposé au rayonnement peut être représentée de la manière la plus complète par les valeurs de tous les paramètres du milieu (conductibilité éléctrique, permittivité, perméabilité magnétique et résistance de pertes magnétiques) , mais dans ce but il faut avoir une représentation complète des composantes magnétique ou électrique à l'aide de six caractéristiques du champ (soit elliptiques, soit orthogonales, soit leur combinaison). Du fait que l'appareillage assurant l'enregistrement de six caractéristiques du champ est compliqué au point de vue conception et exploitation et se caractérise par une faible rapidité de fonctionnement, on est obligé d'avoir recours, pour les besoins pratiques, non pas à une repré15 sentation complète du champ, mais à une représentation partielle. Un compromis optimal entre les impératifs imposés au procédé revendiqué quant à ses qualités économiques et à la faibilité des résultats obtenus est assuré par 20 un mode de réalisation selon lequel on enregistre, lors des mesures périodiques, des séries d'amplitudes de la composante magnétique du champ électromagnétique, qui sont les valeurs de la fonction de l'admittance complexe du milieu exposé au rayonnement. Dans le cas d'un tel mode de réalisation du procédé revendiqué, l'appareillage utilisé pour les mesures est très simple et ne demande pas une formation spéciale. En plus de cela, comme le montrent les essais, la fonction déterminée est très informative sur le plan de l'anisotropie des propriétés électromagnétiques du milieu exposé au rayonnement, et elle assure d'une manière assez efficace l'obtention des données fiables nécessaires à l'établissement des cartes structuraless et géodynamiques. L'invention sera mieux comprise et d'autres buts, dé35 tails et avantages de celle-ci apparaîtront mieux à la lumière de la description explicative, qui va suivre, de plusieurs modes de réalisation donnés uniquement à titre d'exemples non limitatifs, avec références aux dessins annexés dans lesquels: la figure 1 est une représentation schématique du dispositif à induction conçu pour la réalisation du procédé conforme à l'invention; les figures 2 a à 2 c représentent les diagrammes polaires des résultats- de trois mesures des séries d'ampli10 tudes de la composante magnétique du champ électromagnétique induit, effectuées sur le même piquet selon le procédé conforme à l'invention et espacées dans le temps. les figures 3 a et 3 b représentent les diagrammes polaires des résultats de deux mesures des séries d'ampli15 tudes de la composante magnétique du champ électromagnétique, effectuées sur un autre piquet selon le procédé conforme à l'invention et espacées dans le temps; la figure 4 représente une carte structurale et géodynamique du terrain avec indication des limites des zones des mouvements géodynamiques actuels de différents types, cette carte étant établie sur la base de données de mesures effectuées selon le procédé conforme à l'invention; la figure 5 est une vue fragmentaire de la même carte; la figure 6 est une vue fragmentaire d'une carte d'intensité des mouvements géodynamiques actuels, utilisant des isoclines et établie sur la base de données de mesures effectuées selon le procédé conforme à l'invention; la figure 7 est une vue fragmentaire d'une carte d'intensité des mouvements géodynamiques actuels, utilisant des isolignes et établie sur la base de données de mesures effectuées selon le procédé conforme à l'invention; les figures 8 a à 8 c représentent les diagrammes de variation du champ d'activité géodynamique actuelle pour trois secteurs de glissement actif, sur la base de données de mesures effectuées à l'aide du procédé instrumental connu, du procédé connu basé sur l'appréciation de la capacité d'émanation des sols et du procédé conforme à l'inven5 tion. Pour détecter les mouvements géodynamiques actuels dans un massif montagneux à l'aide du procédé conforme à l'invention, on induit un champélectromagnétique dans le milieu attenant à la surface du massif montagneux, à l'aide 10 d'un dispositif à induction montré schématiquement sur la figure 1. Ce dispositif comporte une source de rayonnement électromagnétique 1, qui peut être constitué par n'importe quel générateur harmonique, et un récepteur 2, dis15 posés aux extrémités d'une barre diélectrique 3 La barre 3 est articulée par l'une de ses extrémités sur une plateforme tournante 4 montée de manière à pouvoir tourner autour de l'axe vertical d'une plate-forme d'appui 5 L'extrémité libre de la barre 3 porte un manche 6 prévu pour faciliter 20 la manoeuvre, et un appui amortisseur 7 La barre 3 est réalisée en plusieurs pièces (non représentées et peut être réglée en longueur En vue de l'orientation initiale du dispositif, il est prévu un compas magnétique 8 fixé sur la barre 3. Lorsque le massif montagneux à étudier comporte une couche superficielle d'un milieu à faible temps de relaxation (plusieurs heures), on crée un champ électromagnétique directement dans cette couche, à une profondeur qui ne dépasse pas la profondeur de sa disposition naturelle Parmi 30 les milieux de ce type on trouve plusieurs roches de couverture (argile, sol argileux, terre noire), la glace des nappes d'eau, certains types de revêtements artificiels (asphalte, béton, bitume) La profondeur d'étude est, généralement, de quelques décimètres à 1 mètre, et peut atteindre I 6 2,5 mètres au maximum Dans chaque cas particulier, on choisit cette profondeur selon les conditions réelles et on la règle d'une manière connue en soi en faisant varier la fréquence du générateur ou, lorsque la fréquence est constante, en faisant varier la distance entre la source de rayonnement électromagnétique 1 et le récepteur 2. Si la couche superficielle du massif montagneux ne se prête pas, ou est tout à fait inapte, de par ses propriétés, aux mesures des caractéristiques du champ élec10 tromagnétique (c'est le cas d'une couche superficielle constituée par des rochers, dont le temps de relaxation se chiffre par plusieurs années, ou par des roches incompétentes, par exemple des roches friables), on ménage, sur le piquet du profil, un palier constitué par le milieu 15 à contrôler mis en contact avec la surface du massif montagneux, et on créé un champ magnétique dans ce milieu de contrôle à une profondeur correspondant à la couche déposée. On choisit le milieu de contrôle parmi les substances amagnétiques, homogènes et isotropes dont le point de transfor20 mation de la phase solide en phase liquide se trouve dans l'intervalle de changement natureldes températures On accorde la préférence à ceux des matériaux qui sont disponibles ou sont faciles à transporter Ce sont, en particulier, les sols argileux humides, la glace et l'asphalte On peut également utiliser le béton, le bitume et d'autres matériaux analogues. On place le dispositif sur le profil du secteur à étudier de manière que l'axe géométrique de rotation de la plate-forme tournante soit confondu avec le piquet La source de rayonnement électromagnétique I se trouve également sur ce même axe, ou bien, lorsqu'on utilise un enregistreur automatique, on place sur cet axe le récepteur 2 tandis qu'on dispose la source 1 à l'extrémité opposée de la barre 3 On oriente, à l'aide du compas magnétique 8, la barre 3 en position de départ de manière que son extrémité libre marque la direction Nord On règle la longueur de la barre 3 en fonction des conditions données, en particulier en fonction de la profondeur d'étude requise La 5 longueur maximale de la barre étant atteinte, la base du dispositif (c'est-à-dire la distance entre la source de rayonnement électromagnétique 1 et le récepteur 2) constitue m L'utilisation de dispositifs dont la base dépasse 5 m n'est pas avantageuse, car la porfondeur d'étude)qui augmente 10 avec l'accroissement de la base (si cette profondeur n'est pas limitée par celle de la couche de couverture), atteint, lorsque la limite indiquée est dépassée, une valeur pour laquelle les résultats des études peuvent être affectés par l'hétérogénéité naturelle du milieu à étudier En plus de cela, la rapidité du procédé diminue fortement et, si l'on a recours à un milieu de contrôleles frais et la quantité de travail pour transporter et placer ce milieu sur le piquet augmentent. En principe, dans tous les cas, on tâche de réduire 20 la base à un minimum Mais il faut tenir compte de ce que, pour certains types de milieux à étudier, la valeur minimale de la profondeur d'étude doit dépasser la profondeur de la couche adjacente à la surface, qui s'avère incapable de jouer le rle d'un porteur d'informations géodynamiques par suite 25 de l'érosion et/ou de l'ameublissement. Apres avoir réglé les parties de la barre 3 à la longueur désirée, on les fixe et on fait tourner la barre 3 à l'aide du manche 6 en faisant tourner le récepteur 2 (ou la source de rayonnement électromagnétique) suivant une 30 circonférence, avec arrêts périodiques (dans le cas de l'enregistreur automatique, la rotation se fait d'une manière continue) Le pas de rotation angulaire, déterminé à l'aide d'un limbe fixé sur la plate-forme d'appui 5 dudispositif, est généralement de 5 à 30 O Aux points de la circonférence 35 qui correspondent auxdits arrêts, on mesure à l'aide de l'enregistreur du récepteur 2 l'amplitude de la composante magnétique du champ électromagnétique du milieu exposé au rayonnement. Les résultats des mesures simultanées de séries d'am5 plitutes pour chaque piquet du profil sont enregistrée sous la forme d'un tableau ou d'un diagramme polaire Sur ce diagramme (figures 2), chacune des amplitudes obtenues est représentée sous forme d'un vecteur dirigé du centre des coordonnées vers la périphérie, sous un angle par rapport 10 à l'axe des ordonnées, cet angle étant égal à l'azimut du point de mesure correspondant La ligne brisée reliant, les extrémités de tous les vecteurs d'un tel diagramme correspond approximativement à la caractéristique qualitative de l'anisotropie azimutale des propriétés électromagnétiques du milieu En réalité, le diagramme représente l'anisotropie de l'admittance complexe de cemilieu, admittance dont la fonction est l'amplitude de la composante magnétique du champ électromagnétique Comme le montrent les essais, l'admittance complexe est le paramètre le plus informatif en ce qui concerne 20 les propriétés électromagnétiques du milieu, tandis que l'amplitude de la composante magnétique du champ électromagnétique est la caractéristique qui se prête le mieux aux mesures sans nécessiter un appareillage compliqué et onéreux. Cependant, il y a lieu de noter qu'en principe il est impossible 25 de déterminer l'anisotropie des propriétés électromagnétiques par l'intermédiaire d'autres caractéristiques du champ En particulier, la mesure des amplitudes de la composante électrique du champ électromagnétique n'est pas moins informative. C'est l'absence d'un appareillage de mesure suffisamment efficace qui empêche l'utilisation d'un tel mode de réalisation du procédé. On répète les mesures décrites ci-dessus dans un intervalle de temps donné Le nombre de mesures de retour et l'appréciation des résultats obtenus sont dictés par l'appli- cation concrète des recherches. Si les recherches ne visent qu'à différencier les. mouvements géodynamiques actuels suivant le type de déformations des roches dues à ces mouvements, il suffit d'effectuer les mesures deux fois avec un intervalle dépassant le temps de relaxation du milieu à étudier et de comparer les caractéristiques qualitatives détectées représentées par les séries d'amplitudes de la composante magnétique du champ électromagnétique (mises en tableau 10 ou représentées sous forme de diagrammes polairee) En cas de constance de la caractéristique Qualitative, qui S 1 exprime graphiquement par la similitude géométrique des figures obtenues par des mesures espacées dans le temps et effectuées sur les mêmes piquets (voir figures 2 a et 15 2 c), on constate la présence de mouvements géodynamiques correspondant à une déformation du type compression-traction des roches, et lorsque la caractéristique qualitative est variable, ce qui s'exprime par une altération de la similitude géométrique (figures 3 a et 3 b), on constate la pré20 sence de mouvements géodynamiques actuels correspondant à une déformation du type cisaillement des roches. Les conclusions concernant la constance ou les variations de la caractéristique qualitative peuvent être tirées sur la base d'un traitement des valeurs de mesures par ordi25 nateur à l'aide des méthodes de corrélation de phase paire. Comme le montrent les nombreux essais effectués, le maximum d'hétérogénéité des mouvements du type déformation des roches est observé là o ont lieu des altérations tectoniques, ainsi que dans les zones subissant leur influence, 30 dans les zones d'évolution karstique, de glissements, et dans les secteurs de mouvements technogènes De ce fait, la différenciation des mouvements suivant le type des déformations des roches assure un levé efficace de cartes avec indication des altérations tectoniques, karstiques, des 35 glissements, des zones de mouvements technogènes. Sur la base des mesures effectuées deux fois de la manière décrite, on peut apprécier l'intensité relative des mouvements géodynamiques actuels dus au cisaillement des roches (d'après l'angle de rotation du vecteur de l'amplitude maximale dans un intervalle de temps compris entre deux mesures). La figure 4 représente une carte d'un secteur réel, établie sur la base des résultats de deux mesures effectuées de la manière décrite ci-dessus Les parties non hachurées 10 sur cette carte correspondent aux mouvements géodynamiques actuels liés à une déformation du type compression-traction des couches, la hachuration inclinées indiquant les zones de cisaillement le plus intense. L'établissement de telles cartes peut constituer le 15 but final des recherches, mais dans plusieurs cas cette carte n'est qu'une étape préliminaire des travaux visant la détermination de l'intensité des mouvements géodynamiques actuels Cette carte permet au chercheur de concentrer les études ultérieures sur les zones les plus informatives 20 de la surface à étudier, ce qui rend les travaux plus économiques et efficaces. Lorsqu'il est nécessaire d'apprécier l'intensité des mouvements géodynamiques, on procède à des mesures multiples, le milieu exposé au rayonnement étant porté à un état 25 quasi-isotrope avant chacune des mesures successives lors des mesures périodiques si les résultats de la mesure précédente sont sensiblement différents, c'est-à-dire s'ils différent l'un de l'autre d'une valeur dépassant celle imposée, cette dernière étant choisie en fonction des propriétés physiques du milieu considéré (en moyenne 10 % environ) Ceci assure l'effacement des informations accumulées et l'obtention d'un étalon servant à apprécier les variations de l'anisotropie de l'admittance complexe, d'après lesquelles on peut avoir une idée de l'intensité des processus géo35 dynamiques actuels. La mise a l'état quasi-isotrope de la couche superficielle ou du milieu de contrôle est effectuée, en fonction des propriétés du milieu, à l'aide d'un ameublissement ou d'une action thermique de courte durée. Pour l'appréciation quantitative des mouvements géodynamiques actuels, liés à une déformation du type cisaillement des roches, on choisit l'amplitude maximale dans chaque quadrant du diagramme polaire correspondant au quadrant de la circonférence formée par les points de mesure, et on compare les angles formés par les vecteurs choisis (mesures espacées dans le temps) avec l'axedes ordonnées. Ces angles d'inclinaison correspondent à des positions azimutales des points de mesure, Le rapport entre la différence des angles à comparer et l'intervalle de temps séparant les mesures à comparer constitue la vitesse de rotation des vecteurs choisis, qui caractérise l'intensité des mouvements géodynamiques actuels correspondant au cisaillement des roches. Pour l'appréciation quantitative des mouvements géo20 dynamiques actuels liés à une déformation du type compression-traction des roches, -on choisit les amplitudes maximale et minimale dans chacun des quadrants du diagramme polaire et on compare les rapports de ces valeurs obtenues par les mesures espacées dans le temps L'intensité d'un tel genre de mouvements géodynamiques actuels est appréciée d'après la vitesse de variation desdits rapports, vitesse obtenue sous forme d'un quotient résultant de la division desdits rapports (pour les mesures à comparer espacées dans le temps) par l'intervalle de temps compris entre ces mesures. L'appréciation qualitative des mouvements géodynamiques actuels revêt une importance particulière lors de l'élaboration des projets d'ouvrages de génie civil et est étroitement liée à la prévision des phénomènes séismiques, des coupes de roche et des éruptions de puits. Pour une meilleure compréhension du procédé qui vient d'être décrit, on donne ci-dessous plusieurs exemples concrets mais non limitatifs de réalisation de l'invention. Exemple 1 Les études ont été effectuées dans la zone d'une structure tectonique profonde détectée lors de travaux de forage sous la couverture de roches sédimentaires de 500 m d'épaisseur, y compris 70 m de dépôts incohérents Le secteur à étudier présentait une nappe d'eau d'une épaisseur de 6 à 10 8 mètres Les travaux ont été effectués en hiver, sur de la glace d'une épaisseur de 200 à 300 mm. Aux piquets d'un réseau de profils géophysiques d'un pas de 10 m, on a étudié le champ électromagnétique engendré dans la glace et dans la couche superficielle de l'eau à l'aide d'un générateur harmonique faisant partie d'un'dispositif à induction présentant une faible base, décrit ci- dessus et représenté sur la figure 1 La distance entre le générateur et le récepteur était de 5 m, la fréquence de travail du générateur était de 2 m Hz, ce qui assurait une profondeur 20 du sondage électromagnétique de 2 m, la sensibilité du récepteur étant de 0,5 m V/m. Pour enregistrer les amplitudes de la composante magnétique du champ électromagnétique induit, on à utilisé un récepteur que l'on a déplacé suivant une circonférence par rapport au générateur disposé au point du profil à pas angulaire de 30 Cette opération a été effectuée par un seul opérateur pour 100 piquets durant la journée. Le levé de retour a été effectué au bout de vingt-quatre heures, ce qui dépassait le temps de relaxation du milieu à 30 étudier. Les résultats des mesures pour trois piquets d'un profil, référencés en 73, 76, 79, sont repris dans le tableau 1,o les chiffres romains I et II indiquent le premier levé et le deuxième levé, respectivement. Tableau 1 15 20 25 30 Levé I II Piquets 73 76 79 73 76 79 Azimut O 17 38 21 14 20 18 40 29 27 33 13 24 20 17 50 17 11 47 23 14 18 18 10 16 31 10 12 30 12 14 27 12 30 22 38 27 35 40 42 33 42 39 210 25 45 23 31 34 12 240 32 30 11 27 14 18 2 _ 70 16 15 17 14 12 16 300 42 14 30 37 11 20 330 30 30 20 29 28 18 360 16 36 21 14 21 17 Tableau 2 Piquets 73 76 79 Rang I II I II I II 1 300 300 210 180 60 60 2 30 30 180 150 180 180 3 180 180 O 210 150 150 4 240 210 240 330 300 30 120 120 120 330 30 300 6 330 330 30 240 210 O 7 150 240 60 30 O 240 8 210 150 270 120 330 330 9 90 90 150 270 90 90 60 60 90 60 270 270 11 O O 300 300 120 120 12 270 270 120 90 240 210 2.4 Dans le tableau 2, les mêmes résultats sont représentés sous une forme commode pour la comparaison avec introduction d'un élément de corrélation de rang. Lors de l'analyse et la comparaison des caractéristiques qualitatives de l'anisotropie azimutale de l'admittance complexe du milieu à étudier, obtenues à l'aide de relevés espacés et représentées par les séries d'amplitudes de la composante magnétique du champ électromagnétique pour chaque piquet, ce sont les zones dans lesquelles ont été détectés des 10 indices nettement marqués de rotation des axes d'anisotropie, détectée par le premier levé, qui ont été rapportées aux zones de mouvements géodynamiques actuels liés à une déformation du type cisaillement des roches, tandis que les autres zones, dans lesquelles une telle rotation n'avait pas eu lieu, 15 ont été rapportées aux zones de mouvements géodynamiques liés à une déformation du type compression-traction des roches. Il résulte du tableau 2, en particulier, que le déséquilibre qualitatif de l'anisotropie azimutale lors du levé de retour est nettement marqué dans le cas du piquet 76. La figure 5 représente une vue fragmentaire de la carte du secteur à étudier, établie sur la base de mesures effectuées de la manière décrite ci-dessus Sur ce fragment on voit les piquets 73, 76 et 79, ainsi que les limites entre la zone o les mouvements géodynamiques actuels correspondent 25 aux déformations de cisaillement des roches (hachure inclinée) et la zone o les déformations ont un caractère compressiontraction L'axe de la structure tectonique détectée lors des travaux de forage est indiqué par un trait fort Comme le montre ce fragment, l'hétérogénéité la plus marquée des mou30 vements quant au caractère des déformations résultantes des roches est observée dans la zone de la structure tectonique. Ainsi, la différenciation des mouvements géodynamiques actuels suivant le type des déformations des roches permet de trouver les structures tectoniques. Exemple 2 Les études ont été effectuées dans un secteur caractérisé par une faille formant un chevauchement étendu dans la direction Nord-Est et dont la disposition a été précisée grâce aux travaux miniers L'épaisseur des dépôts incohérents se chiffrait par des dizaines de mètres La partie superficielle des dépôts était constituée par des limons sableux avec de faibles inclusions de cailloux arrondis et de gros galets Les études ont été effectuées sur les piquets d'un 10 réseau de profils géophysiques avec un pas de 5 m à l'aide du même dispositif a induction que dans l'exemple 1 La base du dispositif, lors des essais de tous les profils (à l'exception du profil 3-3 '), était de 2,0 m, et la fréquence de travail était de 3,0 m Hz Le profil 3-3 ' passait partiel15 lement le long de l'axe d'une route asphaltée Le reste du même profil a été asphalté sur une largeur de 2,5 m avant de procéder aux essais La base du dispositif, lors de l'étude de ce profil, constituait 0,8 m, la fréquence de travail était de 3 m Hz La profondeur du sondage dans ce profil était 20 réduit à quelques décimètres, ce qui était suffisant, car ce n'était que la couche d'asphalte qui était porteuse d'informations. Sur la base des résultats de mesures opérées pour chaque piquet avec un pas angulaire de 30 de la même manière que 25 dans l'exemple 1, et des résultats de l'analyse des données obtenues, analyse identique à celle décritedans l'exemple 1, on a établi une carte (schéma) du secteur étudié, qui est représentée dans la figure 4 Les parties hachurées de ce schéma correspondent aux mouvements géodynamiques actuels liés aux déformations du type cisaillement La hachure inclinée correspond aux zones o ces mouvements sont particulièrement intenses, ce qui peut être déterminé d'après la valeur absolue de la différence des nombres donnés dans les colonnes I et II du rang correspondant dans un tableau établi comme le tableau 2, au stade d'étude auquel se réduit l'exemple décrit ci-dessus. En se servant de cette carte, on a précisé les limites du secteur considéré comme valable pour les recherches ulté5 rieures (sur la carte, c'est la zone avoisinant le profil référencé en 3-3 '), après quoi on a procédé aux mesures de travail dans les limites du secteur indiqué, en utilisant la technique décrite dans l'exemple 1 Sur la base des résultats de ces mesures, se rapportant au piquet 17 du profil 3-3 ', on 10 a établi le diagramme polaire montré sur la figure 3 a La ligne brisée fermée reliant les extrémités de tous les vecteurs est identique à la caractéristique qualitative de l'anisotropie azimutale de l'admittance complexe du revêtement asphaltique Dans chacun des quatre quadrants du diagramme polaire 15 ont été révélés les axes d'anisotropie correspondant aux directions des vecteurs des amplitudes présentant la valeur absolue la plus grande (sur la figure 3 a, ils sont indiqués en traits forts) Lesdits axes sont, en l'occurrence, les directions suivant les azimuts de 90 , 240 , 330 (l'axe 20 suivant l'azimut de 90 est commun aux deux quadrants). Après avoir effectué les mesures, on a procédé à un chauffage de courte durée de l'asphalte sur le piquet 17 dans un rayon de 1 m Tout de suite après cette opération on a effectué les mesures de contrôle Ces mesures ont révélé un écart des valeurs absolues des amplitudes mesurées de l'ordre de 10 %, ce qui a donfirmé l'état quasi-isotrope du revêtement asphaltique dans la zone des mesures. Quelque_ temps après (plus exactement, 16 heures 45 minutes après les mesures de travail), on a répété les mesures de travail, dont les résultats (pour le piquet 17) sont reflétés sur le diagramme polaire (figure 3 b). L'analyse de ce diagramme a permis de déceler quatre axes d'anisotropie correspondant aux azimuts suivants: , 120 , 210 et 300 (indiqués en traits forts). La comparaison des diagrammes polaires représentés sur les figures 3 a et 3 b montre que trois axes d'anisotropie sur quatre sont des axes décelés lors des mesures précédentes, tournés d'un angle de 30 dans le sens antihoraire. L'intensité de cisaillement des roches a été déterminée d'après la vitesse de rotation des axes d'anisotropie, obtenue par division de ladite valeur de l'angle de rotation des axes d'anisotropie par l'intervalle de temps compris entre lesdites mesures de travail Cette vitesse était de 1,85 degrés/h, ce qui témoignait de la faible intensité des mouvements géodynamiques actuels correspondant au cisaillement des roches (comme le montrent les recherches faisant appel au procédé revendiqué, la vitesse de rotation des axes peut atteindre 10 degrés/h). Sur la base des mesures effectuées sur les piquets de tout le secteur étudié, on a établi une carte d'intensité des mouvements géodynamiques actuels correspondant au cisaillement des roches avec utilisation des isoclines Un fragment de cette carte, o sont montrés le profil 3-3 ' et le piquet 20 17, est représenté sur la figure 6 Les chiffres indiqués sur les insoclines correspondent à la vitesse de rotation des axes d'anisotropie en degrés/h Le trait fort indique l'axe d'une structure tectonique détectée lors des travaux miniers Comme le montre ce fragment, un maximum d'intensité 25 des mouvements géodynamiques actuels est observé dans la zone de la structure tectonique, dont la disposition peut donc être déterminée, à l'aide du procédé décrit, avec un degré de précision élevé. Exemple 3 Le secteur à étudier se trouvait dans la zone d'une faille L'épaisseur des dépôts incohérents dépassait 10 m. La couche superficielle était constituée par des sols argileux humides du type loesseo Les essais ont été effectués sur un réseau de profils géophysiques pour des piquets d'un pas de 10 m à l'aide d'un dispositif à induction identique à celui utilisé dans l'exemple 1 o La base du dispositif constituait 1,6 m, la fréquence de travail était de 100 k Hzo Un sondage efficace de la couche superficielle dudit milieu, en présence des paramètres indiqués du dispositif, a été assuré sur une profondeur de 0,4 m. Sur la base des résultats des mesures effectuées deux fois à l'aide du procédé décrit dans l'exemple 1, on a établi une carte (schéma) identique à celle qui est représentée sur 10 la figure 4 En utilisant cette carte, on a choisi dans le secteur une zone o prédominaient les mouvements géodynamiques actuels liés aux déformations du type compression-traction des roches Sur chaque piquet, dans les limites de la zone indiquée, a été creusé un évidement de 2,0 m de rayon et d'une profondeur 15 de 0,6 m Apres ameublissementle sol a été placé dans l'évidement (dans le même ordre de succession des couches) et soigneusement damé Les mesures effectuées immédiatement après cette opération ont montré que l'écart maximal des valeurs des amplitudes était de 6 %, ce qui était inférieur à la limite 20 admissible ( 10 %) pour l'état quasi-isotrope des sols argileux humides, compte tenu des erreurs de mesure et des propriétés du milieu Le diagramme polaire établi sur la base des résultats des mesures pour le piquet 6 du profil 14 14 ' est représenté sur la figure 2 b. Quelques temps après (plus exactement, 17,5 heures après), on a répété les mesures dont les résultats sont reflétés (pour le même piquet) sur le diagramme polaire de la figure 2 c. On a choisi, dans chaque quadrant du diagramme polaire, l'amplitude maximale Hmax et l'amplitude minimal Hmin) et on a comparé les rapports max obtenus par des mesures espacées H min dans le temps L'intensité des mouvements géodynamiques actuels a été appréciée suivant la vitesse de variation dudit rapport. On donne ci-dessous les résultats des deux mesures décrites et 10 les résultats de l'analyse pour le quatrième quadrant, qui est le plus illustratif, du diagramme polaire du piquet 6 du profil 14-14 ' (tableau 3). Tableau 3 Intervalle Intensité Numéro H max de temps du mouvedela HH de la H Max Hmin H entre les ment géomesure mn mesures,h dynamit que, h-1 H max min 1 18 17 1,06 17,5 1,81 = O 1 17,5 2 26 14,33 1,81 La valeur obtenuelégale à 0,1,témoigne d'une faible 15 intensité des mouvements. Sur la base des données et de l'analyse des résultats des mesures effectuées sur tous les piquets de la zone étudiée, on a établi une carte d'intensité des mouvements géodynamiques actuels correspondant aux déformations du type compression20 traction des roches (avec utilisation des isolignes) Une vue fragmentaire de cette carte dans laquelle est indiqué le piquet du profil 14-14 ', est représentée sur la figure 7. Les chiffres indiqués sur les isolignes expriment l'intensité par la méthode décrite ci-dessus Le trait fort correspond à 25 l'axe d'une structure tectonique détectée lors de travaux miniers Comme dans l'exemple 2, la zone d'influence de la structure tectonique se caractérisait par l'intensité de mouvements géodynamiques actuels la plus élevée. La fiabilité des conclusions obtenues à l'aide du procédé conforme à l'invention est vérifiée par la comparaison des résultats des observations dans le temps avec utilisation de trois procédés différents, ce qui est illustré par les figures 8 a à 8 c montrant les courbes de variations des champs géodynamiques pour trois secteurs de glissement actif sur la base de résultats de mesures directes instrumentales (trait interrompu mixte), de données de mesures des émanations de l'air tellurique (trait discontinu) et de l'anisotropie de l'admittance complexe du milieu (traitcontinu) Les courbes sont établies sur la base des résultats du traitement des données par ordinateur (à l'aide des méthodes de la corrélation paire) La comparaison révèle un parfait synchronisme de la variation des champs géodynamiques et du signe d'accroissement avec le déroulement des déformations du massif montagneux. 10 Le procédé de détection des mouvements géodynamiques actuels qui vient d'être décrit est économique, simple, très informatif, il peut trouver une large application L'image de l'activité des mouvements géodynamiques actuels dans le champ de l'admittance azimutale des milieux est plus nette, 15 la gamme de processus géodynamiques étudiés est plus étendue et, en plus de cela, l'interprétation est plus concrète en comparaison du procédé connu basé sur la détermination de la capacité d'émanation des gaz telluriques Il faut noter que les études sont effectuées dans des milieux qui sont d'un 20 accès facile et largement répandus, y compris les milieux peu informatifs dans le cas du procédé connu: les sols gelés, rocheux, la glace dont sont recouvertes les nappes d'eau, les revêtements artificiels tels que l'asphalte et le béton. Les techniques de traitement des données des mesures, admises 25 par le procédé conforme à l'invention, suppriment les facteurs de déséquilibre. Le procédé proposé permet l'étude de la déformation en laboratoire aussi bien qu'en campagne En outre, à l'aide de ce procédé, on peut prévoir les phénomènes géodynamiques dans 30 les carrières et les mines (éruption, formation d'une pente raide, etc), ainsi que les phénomènes séismiques. * Sur la base des données obtenues à l'aide du procédé revendiqué on peut établir des cartes structurales et géodynamiques jusqu'à l'échelle 1:500, ce qui facilite le rattachement géodésique et l'élaboration de projets d'ouvrages de génie civil R E V E N D I C A T I O N S 1 Procédé de détection des mouvements géodynamiques dans un massif montagneux, du type consistant à enregistrer et à apprécier les variations des paramètres du milieu attenant à la surface du massif montagneux, caractérisé en ce que l'on crée préalablement dans ledit milieu un champ électromagnétique à l'aide d'une source de rayonnement électromagnétique ( 1), on enregistre les variations des paramètres du milieu exposé audit rayonnement en comparant les résultats de mesures périodiques des séries de valeurs de la fonction 10 d'au moins l'un des paramètres électromagnétiques de ce milieu en des points disposés suivant une circonférence de rayon allant jusqu'à 5 m par rapport à son centre situé en un point d'un profil, et on apprécie les variations des paramètres du milieu exposé audit rayonnement d'après les variations de la 15 caractéristique qualitative de l'anisotropie azimutale du paramètre électromagnétique à étudier du milieu exposé au rayonnement, caractéristique représentée par une série de valeurs de la fonction dudit paramètre enregistrées lors desdites mesures périodiques, la présence des mouvements 20 géodynamiques actuels provoquant une déformation du type compression-traction des roches étant constatée en présence d'une caractéristique qualitative stable, tandis que la présence de mouvements géodynamiques actuels conduisant à une déformation du type cisaillement des roches est cons25 tatée en présence d'une caractéristique qualitative variable. 2 Procédé suivant la revendication 1, caractérisé en ce que l'on porte le milieu exposé au rayonnement à un état quasi-isotrope avant chacune des mesures successives faisant 30 partie des mesures périodiques, si les valeurs de la fonction du paramètre électromagnétique à étudier du milieu exposé au rayonnement, obtenues lors des mesures précédentes, sont sensiblement différentes. 3 Procédé suivant l'une des revendications 1 et 2, caractérisé en ce que dans chaque quadrant du cercle formé par les points en lesquels on effectue les mesures périodiques des séries de valeurs de la fonction du paramètre électromagnétique du milieu exposé au rayonnement, on choisit la valeur maximale de ladite fonction pour la série considérée, on compare les positions azimutales des points auxquels cor10 respondent, dans un même quadrant, les valeurs maximales choisies, obtenues lors de mesures espacées dans le temps, et, d'après le rapport entre la différence des azimuts des -points à comparer et l'intervalle de temps entre les mesures à comparer, on apprécie l'intensité des mouvements qui provo15 quent ladite déformation des roches. 4 Procédé suivant l'une des revendications 1 et 2, caractérisé en ce que dans chaque quadrant du cercle formé par les points en lesquels on effectue les mesures périodiques des séries de valeurs de la fonction du paramètre électro20 magnétique du milieu exposé au rayonnement, ou choisit les valeurs maximale et minimale de ladite fonction pour la série considérée, on compare les rapports de ces valeurs obtenues lors des mesures espacées dans le temps, et, d'après la vitesse de variation desdits rapports, on apprécie l'inten25 sité des mouvements géodynamiques actuels qui provoquent la déformation du type compression-traction des roches. Procédé suivant l'une des revendications 1 à 4, caractérisé en ce que l'on utilise, en tant que milieu à exposer au rayonnement, la couche superficielle du massif 30 montagneux. 6 Procédé suivant l'une des revendications 1 à 4, caractérisé en ce que l'on utilise, en tant que milieu à exposer au rayonnement, un milieu de contrôle mis en contact avec la couche superficielle du massif montagneux et constitué par une substance amagnétique présentant l'homogénéité et l'isotropie du paramètre électromagnétique à étudier et dont le point de transformation de la phase solide en phase liquide se trouve dans l'intervalle de variation naturelle des tempé5 ratures. 7 Procédé suivant la revendication 6, caractérisé en ce que l'on utilise, en tant que milieu de contrôle, de la glace. 8 Procécé suivant la revendication 6, caractérisé en ce 10 que l'on utilise, en tant que milieu de contrôle, l'asphalte. 9 Procédé suivant l'une des revendications 1 à 8, caractérisé en ce que l'on enregistre, lors des mesures périodiques, les séries d'amplitudes de la composante magnétique du champ électromagnétique, qui sont les valeurs de la fonction de 15 l'admittance complexe du milieu exposé au rayonnement.