Brevet FR2648562A1 — publié le 19901221

L'invention est relative à un procédé et à un système de détermination de l'anisotropie azimutale dans un domaine d'un milieu matériel. L'étude de l'isotropie des milieux matériels présente actuellement un intérêt croissant. Dans le domaine technique de la prospection géophysique ou minière, il est particulièrement intéressant de pouvoir déterminer les paramètres d'anisotropie de milieux tels que les couches géologiques sousjacentes à la surface de l'écorce terrestre. Des travaux récents ont été effectués, en France notamment, afin d'établir une méthodologie permettant la détermination des caractéristiques d'anisotropie des milieux matériels à partir d'ondes de cisaillement transversales, ondes S ou "shearwaves" en langage anglo-saxon. Les méthodologies proposées permettent pour un milieu matériel donné de déterminer les matrices de transfert du milieu aux ondes de cisaillement précitées. Pour une description plus détaillée de la méthodologie proposée, on pourra se reporter à la demande de brevet français n" 2 599 508 publiée le 04.12.1987 au nom de la demanderesse. La méthodologie précédemment citée donne toute satisfaction. Cependant, elle nécessite la mise en oeuvre de sources génératrices d'ondes de cisaillement et de récepteurs de ces ondes en fonction de leur direction de polarisation, par exemple des géophones omnidirectionnels tri axes. Enfin, la méthodologie précitée est avant tout avantageuse en sismique de puits pour la détermination de la fonction de transfert d'une zone déterminée de terrain. La présente invention a pour but la mise en oeuvre d'un procédé et d'un système de détermination de l'anisotropie azimutale d'un milieu matériel à. partir d'ondes de pression, donc longitudinales. Un autre objet de la présente invention est la mise en oeuvre d'un procédé et d'un système de détermination de l'anisotropie azimutale d'un milieu matériel pouvant être mis en oeuvre en sismique terrestre de surface à partir de matériel de relevé de profils d'acquisition en sismique tridimentionnelle par exemple. Un autre objet de la présente invention est enfin la mise en oeuvre d'un procédé et d'un système de détermination de l'anisotropie azimutale d'un milieu matériel susceptible de permettre l'établissement d'une carte tridimensionnelle d'anisotropie par un processus de type tomographie. Le procédé de détermination de l'anisotropie azimutale dans un domaine d'un milieu matériel présentant au moins une surface de séparation sensiblement plane, objet de la présente invention, est remarquable en ce qu'il consiste, pour une acquisition pour un domaine élémentaire, à disposer sur la surface de séparation une pluralité de couples d'émetteur-récepteur d'une onde de pression dirigée vers le milieu matériel à étudier. Chaque couple d'émetteur-récepteur est disposé selon une direction déterminée par rapport à une direction de référence autour d'un centre commun et un émetteur est décalé d'un récepteur correspondant d'une distance de décalage. Une onde de pression est engendrée à partir d'un émetteur et les échos ou ondes de pression réfléchies sont enregistrés au niveau du récepteur qui lui est associé dans une fenêtre de temps de propagation définissant le domaine élémentaire. L'étape précédente d'émission et d'enregistrement est répétée pour chacun des couples d'émetteur-récepteur pour le domaine élémentaire considéré. La loi de variation de propagation des ondes de pression réfléchies ou échos en fonction de l'azimut correspondant est établie à partir des fluctuations d'amplitude et de temps de propagation pour chaque écho reçu. Le système de détermination de l'anisotropie azimutale dans un domaine constitué de domaines élémentaires d'un milieu matériel présentant au moins une surface de séparation sensiblement plane, conforme à l'objet de la présente invention, est remarquable en ce qu'il comprend un réseau de couples d'émetteur-récepteur disposés selon une pluralité d'azimuts et des moyens de commande des couples d'émetteur récepteur interconnectés aux couples d'émetteur-récepteur pour assurer la commande en émission des émetteurs et en réception des récepteurs définissant le domaine élémentaire. Des moyens de calcul sont prévus pour effectuer un calcul des fluctuations d'amplitude et de temps de propagation des échos reçus en fonction de l'azimut. Le procédé et le système, objets de l'invention, trouvent application à la détermination de l'anisotropie azimutale des milieux matériels de toute nature, en particulier les milieux matériels stratifiés. Une application particulièrement avantageuse du procédé et du système objets de l'invention concerne la sismique tridimensionnelle terrestre en tectonique .calme. Une description plus détaillée du procédé et du système, objets de l'invention, sera donnée de façon plus détaillée à la lecture de la description et à l'observation des dessins dans lesquels la figure I représente un diagramme illustratif de la mise en oeuvre du procédé objet de l'invention, la figure 2 représente la loi de variation du temps de propagation entre un émetteur Ei et un récepteur Ri en fonction de l'orientation de ces derniers par rapport à une direction de référence, la figure 3 représente la loi de variation de l'amplitude des échos reçus par un récepteur Ri en fonction de l'angle d'orientation par rapport à une direction de référence, la figure 4 représente un diagramme tridimensionnel des signaux d'échos reçus par un récepteur Ri déterminé en fonction de l'axe vertical gradué en temps de propagation, I'axe horizontal gradué en distance de décalage entre émetteur Ei et récepteur Ri, une graduation supplémentaire en angle d'orientation de û à 3600 par rapport à une direction de référence étant en outre portée sur l'axe horizontal pour différentes valeurs discrètes de distances de décalage. la figure 5 représente un diagramme illustratif de la mise en oeuvre du procédé selon l'invention pour différentes profondeurs d'analyse, c'est-à-dire pour différents domaines élémentaires successifs de susceptibles de constituer un domaine total D, .la figure 6 représente un mode de réalisation avantageux non limitatif d'un système permettant la mise en oeuvre du procédé objet de l'invention, . Ia figure 7 représente l'organigramme d'un programme intégré à un calculateur de commande et de traitement permettant la mise en oeuvre du procédé objet de l'invention pour un ensemble de domaines élémentaires de constituant un domaine D. Le procédé, objet de l'invention, sera tout d'abord décrit en liaison avec le diagramme illustratif de la figure I. Selon la figure précitée, et ainsi que représenté sur celle-ci, le procédé, objet de l'invention, consiste, afin de déterminer l'anisotropie azimutale d'un milieu matériel présentant au moins une surface de séparation sensiblement plane, notée S sur cette figure, à disposer sur cette surface de séparation une pluralité de couples d'émetteur-récepteur El, E2, Ei, respectivement RI, R2, Ri d'une onde de pression dirigée vers le milieu matériel à étudier. Afin d'effectuer une acquisition sur un domaine élémentaire de sur la figure I, chaque couple d'émetteur-récepteur Ei, Ri est disposé selon une direction a i déterminée par rapport à une direction de référence N. Les couples d'émetteur-récepteur Ei, Ri sont disposés autour d'un centre commun CK. A titre d'exemple non limitatif, la direction de référence N peut être quelconque mais elle peut, dans le cas où le procédé, objet de l'invention, est utilisé pour l'étude des couches géologiques sousjacentes à la surface terrestre, être constituée par la direction du nord. Dans tous les cas, et quelle que soit la direction de référence utilisée, I'anisotropie mise en évidence pour le milieu matériel étudié sera définie par l'anisotropie azimutale, I'angle a i étant alors considéré comme l'azimut de chaque couple d'émetteur-récepteur Ei, Ri par rapport à la direction de référence N, que celle-ci soit constituée par la direction du nord géographique ou non. o' Ainsi qu'on l'a représenté sur la figure 1, un émetteur Ei est séparé d'un récepteur correspondant Ri d'une distance de décalage Oi. On notera que les distances de décalage des différents couples d'émetteurrécepteur peuvent a priori être quelconques. Cependant, pour la commodité de mise en oeuvre du procédé, objet de l'invention, une valeur prédéterminée de décalage sera utilisée, ainsi qu'il sera décrit ultérieurement dans la description. Conformément à la figure I précitée, le procédé, objet de l'invention, consiste alors à engendrer à partir d'un émetteur Ei une onde de pression Pi et à enregistrer au niveau du récepteur Ri qui lui est associé les échos ou ondes de pression réfléchies dans une fenêtre de temps de propagation définissant un domaine élémentaire de dans le domaine D à étudier. On comprendra,. bien sûr, que l'enregistrement par le récepteur Ri assdcié à l'émetteur Ei correspondant des échos ou ondes de pression réfléchies sur une couche sousjacente à la surface de séparation S, cette réflexion étant, bien entendu, partielle, et l'enregistrement précité, étant effectué pendant une fenêtre de temps de largeur 6 t décalée par rapport à l'instant d'émission de l'onde de pression Pi par l'émetteur Ei d'un temps t, permet d'obtenir en fait l'information d'amplitude et de phase ou de retard de propagation de l'onde réfléchie sur une tranche d'épaisseur o Z correspondant à la largeur de la fenêtre temporelle d'enregistrement 6 t. La tranche d'épaisseur 6 Z sur une surface de réflexion élé.mentaire a déterminée permet de constituer le volume élementaire ou domaine élémentaire , un ensemble de ces domaines élémentaires de constituant le domaine d'étude D précité. On comprendra, bien sûr, qu'en fonction de la largeur de la fenêtre temporelle 6 t, il est possible de rendre l'épaisseur de la tranche de domaine élémentaire de fl Z aussi petite que possible de façon à obtenir un domaine élémentaire de de taille appropriée. En ce qui concerne la surface élémentaire de réflexion a celle-ci est déterminée par des conditions d'incidence de l'onde de pression Pi, ainsi qu'il sera décrit ultérieurement dans la description. Selon une autre caractéristique du procédé, objet de l'invention, I'étape précédente, consistant à engendrer à partir d'un émetteur Ei une onde de pression Pi et à enregistrer au niveau du récepteur Ri associé les échos ou ondes de pression réfléchies, est répétée pour chacun des couples d'émetteur-récepteur Ei,Ri pour le domaine élémentaire de considéré. Les étapes précédentes ayant été répétées pour chaque couple d'émetteur-récepteur précité, le procédé, objet de l'invention, consiste alors à déterminer pour chaque écho reçu les fluctuations d'amplitude et de temps de propagation pour établir une loi de variation de propagation des ondes de pression réfléchies ou échos en fonction de l'azimut ou angle d'orientation correspondant par rapport à la direction de référence. Dans le cas où la surface élémentaire a est située au-dessus de terrains ou de couches matérielles présentant l'anisotropie autour de l'axe vertical passant par le centre commun CK, I'anisotropie précitée se traduit par une déformation en temps des réflexions correspondant à des trajets traversant les terrains anisotropes suivant une sinusolde de période 1800 et par une variation de l'amplitude de ces réflexions suivant une sinusotde également de période 1800. Compte tenu de l'observation précitée, le procédé, objet de l'invention, dans son étape de détermination de la loi de variation de. propagation des ondes de pression réfléchies ou échos en fonction de l'azimut ou angle d'orientation ai par rapport à la direction de référence N consiste alors à mémoriser les valeurs.de temps de propagation Ti pour chaque valeur d'azimut o i, puis à comparer pour une valeur d'azimut a i donnée la valeur du temps de propagation Ti correspondant à la valeur moyenne des temps de propagation T au point considéré. Selon une autre caractéristique du procédé, objet de l'invention, I'étape de détermination de la loi de variation de propagation des ondes de pression réfléchies ou échos en fonction de l'azimut correspondant peut égalemént consister à mémoriser les valeurs d'ampli tude Ai des ondes de pression réfléchies ou échos pour chaque valeur d'azimut a i ou angle d'orientation par rapport à la direction N, puis à comparer pour une valeur d'azimut donnée la valeur de l'amplitude Ai correspondante à la valeur moyenne des amplitudes A au point considéré. Selon une caractéristique avantageuse non limitative du procédé, objet de l'invention, et afin de faciliter l'analyse mathématique des phénomènes précédemment décrits, I'émission et la réception des ondes de pression réfléchies ou échos sont effectuées pour une plage d'azimuts ou d'angles d'orientation a i comprise entre 0 et 1800 à incidence constante. Dans ce cas, ainsi que représenté notamment en figure 1, les couples d'émetteur-récepteur Ei,Ri présentent une même distance de décalage Oi et sont disposés en un motif formant une même circonférence par rapport au centre commun CK. A titre d'exemple non limitatif, les paramètres d'émission, tels que le décalage Oi entre chaque émetteur Ei et chaque récepteur Ri correspondant et, en définitive, I'angle d'incidence des ondes de pression émise Pi sur les surfaces réfléchissantes élémentaires a peuvent être choisies égales de façon que Oi = 2 x Z où Z désigne la profondeur de la surface réfléchissante élémentaire 'a par rapport à la surface de séparation S. Sur la figure 2, on a représenté la loi de variation du temps de propagation des ondes de pression réfléchies Pi pour différents angles d'azimut ou d'orientation a i pour une plage de valeurs de cet angle d'orientation comprise entre 00, correspondant sensiblement à la direction de référence N, par rapport au centre commun CK, et 1800 correspondant à la direction opposée par rapport à ce même centre commun CK. A titre d'exemple non limitatif, les incréments d'angles d'orientation a i ont été pris égaux à 150. Sur la figure 2, I'axe vertical est gradué à une valeur de temps de propagation par rapport à un temps de propagation moyen, et on constate que les valeurs de temps de propagation, en fonction de l'angle d'orientation ou d'azimut a i, suivent sensiblement une loi sinusoldale. Ainsi, la détermination du temps de propagation le plus faible per.met de déterminer, d'une part, la direction de l'axe d'anisotropie dit "axe rapide", cette direction étant notée a r sur la figure 2, et, d'autre part, pour le temps de propagation maximum la direction de l'axe d'anisotropie dit "axe lent" correspondant à la direction a I sur ia figure 2. Les directions rapide et lente précitées constituent les directions principales d'anisotropie du milieu matériel considéré au centre commun CK précité. De la même manière, ainsi qu'on l'a représenté en figure 3, la loi de variation de l'amplitude Ai des ondes de pression Pi en fonction de l'angle d'orientation ou azimut a i des couples d'émetteur-récepteur Ei,Ri a été représentée en figure 3. Sur la figure 3, I'axe vertical représente des valeurs d'amplitude par rapport à l'amplitude moyenne et l'axe horizontal des valeurs d'angle d'orientation par rapport à la direction de référence N. La loi de variation d'amplitude des échos en fonction de l'angle d'orientation a i ou d'azimuts précités suit également une variation sinusoldale. La. loi,de variation d'amplitude, telle que représentée en figure 3, a été obtenue pour des incréments d'orientation a i égaux à 15"'Comme dans le cas de la détermination de la loi de variation de temps de propagation. La loi de variation d'amplitude, telle que représentée en figure 3, permet de déterminer l'atténuation différentielle d'amplitude selon les axes lent et rapide précités pour les valeurs,a r, al précédemment déterminées, c'est-à-dire pour les directions rapide et lente constituant les directions principales d'anisotropie ou coordonnées naturelles du milieu matériel considéré. Bien entendu, en ce qui concerne les figures 2 et 3, les axes de valeurs de temps et d'amplitude peuvent être gradués en valeur relative par rapport à la valeur de temps de propagation. moyenne, respectivement la valeur d'amplitude moyenne. En ce qui concerne la conduite du procédé, objet de l'invention, en particulier, afin d'obtenir des valeurs d'amplitude stables des échos ou ondes de pression réfléchies au cours des essais, il est souhaitable de travailler à incidence constante, I'amplitude des échos réfléchis dépendant, bien entendu, de incidence de l'onde de pression Pi engendrée par chaque émetteur Ei. Ainsi qu'on a pu le constater au cours d'essais, de faibles incidences, c'est-à-dire de faibles valeurs de distance de décalage entre émetteurs Ei et récepteurs associés Ri, ont pour effet de provoquer de très faibles écarts de temps de propagation. En conséquence, la valeur d'incidence précédemment indiquée, correspondant à un angle d'incidence des ondes de pression Pi par rapport à la verticale sensiblement égale à 450, a été retenue. Sur la figure 4, on a représenté un diagramme tridimentionnel des signaux d'échos reçus par un récepteur Ri déterminé en fonction de l'axe vertical gradué en temps de propagation T, cet axe étant gradué en milliseconde. L'axe horizontal est gradué en distance de décalage entre émetteurs Ei et récepteurs Ri pour des distances de décalage exprimées en pied (un pied = 0,3048 mètre), une graduation supplémentaire en angle d'orientation ou d'azimut de 0 à 3600 par rapport à une direction, de référence étant en outre portée sur ce même axe horizontal pour différentes valeurs discrètes de distance de décalage. Sur le diagramme de la figure 4 précité, on peut constater, d'une part, la périodicité selon une loi sinusoîdale des fluctuations du temps de propagation T en fonction de l'azimut pour chaque valeur de décalage correspondante, ainsi qu'une variation périodique de période 1800 des valeurs d'amplitude des échos ou ondes de pression réfléchies, ces valeurs périodiques d'amplitude apparaissant du fait de l'intensité de coloration en noir pour un temps de propagation moyen donné. Sur la figure 4, on a également porté les courbes d'isoincidence des ondes de pression Pi, ces courbes ayant été calculées pour les valeurs d'incidence 200, 300 et 45 , pour l'ensemble des valeurs de décalage entre chaque émetteur Ei et chaque récepteur associé Ri pour des valeurs de profondeur et finalement de temps de propagation T correspondantes. On constate finalement que l'angle d'incidence le plus favorable pour la conduite des essais est compris entre 300 et 45". Le diagramme de la figure 4 justifie pour l'étude du phénomène d'anisotropie, conformément au procédé objet de l'invention, l'étude des fluctuations des temps -de propagation et des amplitudes des échos en fonction de l'angle d'orienta- tion ou d'azimut par rapport à la direction de référence sur une valeur d'orientation ou d'azimut comprise entre 0 et 1800 pour déterminer les directions principales et les atténuations différentielles selon ces directions. Une telle constatation est justifiée par le fait que l'anisotropie de propagation est identique par rotation de 1800 d'un couple d'émetteur Ei récepteur Ri considéré autour d'un centre commun CK, la distance de décalage Oi correspondante étant conservée. Une variante de réalisation ou de mise en oeuvre du procédé, objet de l'invention, sera décrite en liaison avec la figure 5. La variante de réalisation ou de mise en oeuvre du procédé, objet de l'invention, décrite en liaison avec la figure 5 permet d'effectuer une mesure ou une acquisition dans une succession de domaines élémentaires de de profondeur successive incrémentée. Dans ce but, ainsi que représenté de façon schématique en figure 5, le procédé, selon l'invention, consiste à effectuer une pluralité d'acquisitions, le motif formant une même circonférence par rapport au centre commun CK étant transformé par une suite d'homothéties de centre CK et de rapport Kr. Ainsi, le motif initial peut être transformé de façon que la distance de décalage Oir prenne des valeurs successives Oir Kr x Oir Les valeurs du coefficient ou rapport d'homothétie Kr peuvent être choisies de façon que Oir et Oi r-I représentent par exemple les distances de décalage des couples émetteur-récepteu,r Ei,Ri pour un .incrément de profondeur Pf à incidence constante, ainsi que représenté sur la figure 5 et correspondant à un pas de décalage en profondeur du domaine élémentaire de considéré. Sur la figure 5, on a représenté à titre d'exemple non limitatif par rapport à un centre commun CK un motif constitué par une pluralité de couples émetteur-récepteur Ei,R 1, Rital issue d'un motif initial non représenté et suivi successivement d'un motif constitué par des couples émetteur-récepteur Eir, Rir puis Ei +1 > Ri r+l Les domaines élementaires successifs sont notés alors de 1. der, der+ 1 Ils correspondent, bien entendu, à différents domaines élémentaires de surface élémentaire a décalés d'un incrément de profondeur Pf, ainsi que représenté sur la figure 5. Bien entendu, en fonction de l'incidence des ondes de pression Pi sur chaque surface élémentaire a et de la largeur de la feriêtre tempor.elle 6 t, il est ainsi possible d'assurer un recouvrement successif des domaines élémentaires successifs de , der, de 1 pour constituer tout ou partie d'un domaine d'analyse D, ainsi qu'il sera décrit ci-après dans la description. Ainsi qu'on l'a en outre schématisé en figure 5, en vue d'effectuer une acquisition en une succession de domaines élémentaires de décalés en position selon une première x ou deuxième y direction parallèle à la surface de séparation s; le motif, formant une même circonférence par rapport au centre commun CK, ainsi que ce dernier, sont alors décalés en translation sur la surface S d'un ou plusieurs incréments px, py respectivement. Ainsi, et à titre purement illustratif, on a représenté sur la figure 5 les directions x et y précitées, celles-ci étant a priori purement arbitraires sur la surface S. Le centre commun CK peut ainsi être amené en une position ultérieure CK+l à la figure 5, le centre commun ayant ainsi été déplacé d'un pas py selon la direction y puis d'un pas px selon la direction x par exemple. Il est ainsi possible, conformément au procédé, objet de l'invention, d'effectuer en nouvelle position du centre commun CK+I une nouvelle acquisition, telle que précédemment décrite pour la position précédente CK. Un recouvrement de l'ensemble des domaines élémentaires de I à de 1 et des domaines élémentaires connexes obtenus à partir du centre commun CK+1 I pour obtenir une acquisition dans un domaine D constitué par l'ensemble des domaines élémentaires précités. Le procédé, objet de l'invention, est ainsi analogue et particulièrement remarquable en ce qu'il s'analyse en un procédé comparable à un procédé de tomographie permettant une analyse fine de la matière en domaines élémentaires. Bien entendu, le recouvrement des domaines élémentaires successifs, tant dans la direction de profondeur Z que dans les directions de translation latérale, selon les directions x ou y, peut être effectué et déterminé par produits d'intercorrélation d'images de deux domaines successifs, afin de déterminer les limites moyennes en direction Z, X ou Y correspondantes selon les techniques de corrélation d'images utilisées habituellement. Des indications complémentaires seront données relativement à la mise en oeuvre du procédé, objet de l'invention, en ce qui concerne la mise en oeuvre de l'émission de chaque onde de pression Pi considérée. Selon un premier aspect particulièrement avantageux du procédé, objet de l'invention, la répétition pour chaque couple d'émetteurrécepteur Ei,Ri peut être effectuée par une émission-enregistrement séquentielle des ondes de pression réfléchies pour les directions d'azimut ou d'orientation a i considérées. Ainsi, Si l'on se reporte à nouveau à la figure I, le procédé, objet de l'invention, peut ainsi consister à effectuer tout d'abord l'émission par l'émetteur El de l'onde de pression P1 et. conjointement la réception de l'onde de pression réfléchie par la surface élémentaire a du domaine de et la mesure du temps de propagation TI et des valeurs d'amplitude Al de.l'onde de pression réfléchie ou écho, puis successivement pour l'émetteur-récepteur E2, R2, puis l'émetteur-récepte,ur Ei,Ri et ainsi de suite jusqu'à épuisement des couples d'émetteur-récepteur. Le traitement des valeurs enregistrées peut alors être effectué de façon à déterminer les directions principales ou azimuts d'axes lent et rapide, ainsi que mentionné précédemment. Au contraire, dans le cas où le procédé, objet de l'invention, doit être mis en oeuvre par exemple dans le cadre d'une acquisition en sismique tridimentionnelle terrestre, le nombre des couples d'émetteurrécepteur, ainsi qu'il sera décrit ultérieurement dans la description, peut être très important bien que le procédé, objet de l'invention, puisse être mis en oeuvre à partir d'un nombre de couples d'émetteur-récepteur supérieur ou égal à 3. Dans le cas précité toutefois, la répétition pour chaque couple d'émetteur-récepteur Ei,Ri peut, de manière plus avantageuse, être effectuée par une émission-enregistrement simultanée des ondes de pression réfléchies. Dans ce cas, l'émission de l'onde de pression Pi par chacun des émetteur Ei consiste en une onde de pression périodique codée, chaque émetteur Ei émettant une onde de pression Pi codée selon un code qui lui est propre. A titre d'exemple non limitatif, le codage de l'émission de l'onde de pression Pi par chaque émetteur Ei peut consister en un codage de phase. Le codage de phase peut alors être avantageusement représentatif de l'azimut a i ou angle d'orientation correspondant du couple Ei,Ri d'émetteur-récepteur correspondant. A titre d'exemple non limitatif, le codage de phase peut être réalisé selon un code binaire du numéro d'ordre i précité. Le principe de codage de phase à l'émission de sources, de vibrations ne sera pas décrit plus en détail car il correspond à des principes de codage parfaitement connus de l'homme de métier. Une description plus détaillée d'un système de détermination de l'anisotropie azimutale dans un domaine D constitué de domaines élémentaires de d'un milieu matériel, conformément à l'objet de l'invention, ce milieu matériel présentant au moins une surface de séparation sensiblement plane sera maintenant donné en liaison avec les figures 6 et 7. Conformément à la figure 6 notamment, le système, objet de l'invention, comporte un réseau de couples d'émetteur-récepteur El, RI, Ei, Ri à E6, R6 sur cette figure. Bien entendu, ainsi que déjà décrit en liaison avec les figures 1 et 3, illustratives de la mise en oeuvre du procédé, objet de l'invention, les couples d'émetteur-récepteur précités sont disposés selon une pluralité de directions d'orientation a i par rapport à la direction de référence N autour du centre commun CK. Sur cette même figure 6, le réseau de couples d'émetteur-récepteur est noté R. Ce réseau R est matérialisé par le cercle en trait mixte sur la figure 6. A titre d'exemple non limitatif, les émetteurs Ei peuvent être constitués par des sources de vibrations sismique et les récepteurs Ri par des géophones omnidirectionnels. En outre, ainsi qu'on l'a représenté sur cette figure des moyens de commande sont prévus pour effectuer la commande des couples d'émetteur-récepteur Ei,ri interconnectés aux couples d'émetteur-récepteur. Les moyens précités permettent d'assurer la commande en émission des émetteurs Ei et en réception des récepteurs. Ri dans une fenêtre de temps Gt de propagation définissant le domaine élémentaire de, ainsi que décrit précédemment. Sur la figure tridimensionnelle. Ces sources peuvent également être éventuellement constituées par des sources explosives. Dans ce dernier cas, I'émission des ondes de pression Pi peut être réalisée séquentiellement, ainsi que décrit précédemment en liaison avec le procédé, objet de l'invention. Les lignes de commande et de transmission. précitées sont alors reliées à partir de chaque couple Ei,Ri aux moyens de commande 1 par un système dinterfaçage connu en lui-même, par exemple aux moyens calculateurs 100 de façon classique. Le système, objet de l'invention, peut en outre, ainsi que représenté en figure 6, être mis en oeuvre de manière particulièrement avantageuse dans le cas de l'existence d'une grille G d'acquisition de données sismiques en sismique tridimensionnelle terrestre, ainsi qu'il sera décrit ci-après. Un tel cas apparaît particulièrement avantageux dans la mesure où la grille G, ainsi que représenté en figure 6, étant déjà implantée, le système, objet de l'invention, peut alors être mis en oeuvre par configuration d'un réseau R sur la grille G, ainsi qu'il sera décrit ci-après. Dans le cas précité, le réseau R de couples d'émetteurrécepteur Ei,Ri peut alors être avantageusement formé sur une grille G de lignes d'émetteurs El-lk et de récepteurs Rlk, espacés d'un pas Px par exemple. Deux lignes successives sont alors formées par un alignement d'émetteurs respectivement de récepteurs, ces deux lignes étant espacées d'un même pas Py, égal par exemple au pas Px de façon à constituer dans une direction perpendiculaire aux lignes, c'est-à-dire la direction x, un arrangement en colonnes similaire d'émetteurs El Ik-récepteurs Rlk alternés. Chaque couple d'émetteur-récepteur Ei,Ri est alors choisi par adressage de couples d'émetteur-récepteur Ei,lk, Rlk pour constituer un réseau R selon un motif formant sensiblement une circonférence. Bien entendu, dans le'cas précité, le réseau R est également représenté sur la figure 6 par le cercle en trait mixte,.les couples d'émetteur-récepteur El, RI à E6, R6 étant constitués par des émetteurs-récepteurs El.lk, Rlk de la grille G. Dans ce dernier cas et ainsi qu'on l'a également représenté en figure 6, les moyens de commande I des couples d'émetteur-récepteur Ei,Ri interconnectés aux couples d'émetteur-récepteur précités et les moyens de calcul peuvent être constitués respectivement par un circuit d'interface 200 assurant l'interconnexion des lignes d'émetteurs et de-récepteurs. Le circuit d'interface permet d'assurer le déclenchement de l'émission des ondes.de pression Pi par adressage des émetteurs Ei considérés, c'est-à-dire par adressage dans la grille C des émetteurs El I k et des récepteurs Rlk d'adresses choisies en conséquence. Le circuit d'interface 200 permet d'assurer également la réception et l'enregistrement par adressage des récepteurs Ri considérés des ondes de pression réfléchies ou échos dans la fenêtre temporelle 8 t définissant le domaine élémentaire de considéré. Bien entendu, l'adréssage des récepteurs Ri est effectué à partir des adresses Rlk des lignes de. récepteurs considérés. Enfin, le circuit d'interface 200 permet d'effectuer la conversion analogique numérique des signaux d'échos des ondes de pression Pi et la mémorisation des signaux d'échos numérisés obtenus. En outre, un calculateur 100 de conduite d'acquisition et de traitement est interconnecté au circuit d'interface 200 précité afin de permettre la commande par ce dernier des émetteurs Ei et des récepteurs Ri. Le calculateur 100 de conduite d'acquisition et de traitement permet également d'effectuer le traitement des échos enregistrés et numérisés pour déterminer, pour un centre CK considéré, la direction des axes lent et rapide par rapport à la direction de référence N ainsi que l'atténuation différentielle d'amplitude selon les axes lent et rapide précités. Le calculateur de conduite d'acquisition et de traitement peut être constitué de manière classique lorsque le système, objet de l'invention, correspond à un système d'acquisition de données sismiques tridimentionnelles par exemple par un calculateur VAX SYSTEM par exemple muni de ses ressources périphériques. Afin d'effectuer le traitement des échos enregistrés, le calculateur 100 de conduite d'acquisition et de traitement peut comprendre un sous programme de tri des valeurs des temps de propagation Ti permettant de déterminer les azimuts ou angles d'orientation a 1 et a r correspondant, d'une part, au temps de propagation maximum et, d'autre part, au temps de propagation minimum pour le centre commun CK considéré. Ces directions correspondent et définissent respectivement les directions des axes lent et rapide, ainsi que décrit précédemment. En outre, le calculateur de conduite d'acquisition et de traitement comprend également un sous programme de tri des valeurs des amplitudes et de calcul de l'amplitude crête à crête de la variation d'amplitude des échos reçus en fonction de l'azimut ou angle d'orientation a i. Le sous programme précité permet de calculer également l'atténuation différentielle selon les axes lent et rapide des amplitudes des échos réfléchis. Les deux sous programmes de tri précités peuvent être constitués par des modules de tri de type classique, tels que des modules de tri à bulle permettant par. exemple de déterminer la plus grande valeur d'une série de valeurs de temps de propagation mémorisées, ainsi que la plus petite valeur de temps de propagation mémorisée et les valeurs correspondantes d'angle d'orientation ou d'azimut a i. Il en va de même pour les valeurs d'amplitude mémorisées, les valeurs d'amplitude pour les valeurs d'orientation des axes lent et rapide oL 1 a r étant déterminées par simple lecture des amplitudes correspondantes. Enfin, le calcul de l'atténuation différentielle, selon les axes lent et rapide, peut être effectué par simple calcul du rapport des valeurs des amplitudes aux valeurs a 1 et a r précitées par rapport à la valeur d'amplitude moyenne. Ces programmes correspondent à des programmes dont les algorithmes sont parfaitement connus de l'homme de métier. Une description plus détaillée de l'afgorithme d'un programme permettant une conduite automatique du procédé, objet de l'invention, au moyen du système décrit en figure 6 sera donnée en liaison avec la figure 7. Conformément à la figure précitée, en vue d'effectuer une acquisition sur un domaine D constitué par un ensemble de domaines élémentaires de, le calculateur 100 de conduite d'acquisition et de traitement peut comporter un programme de conduite d'acquisition. Ce programme de conduite d'acquisition peut comprendre, ainsi que représenté sur la figure 7 précitée, un module d'initialisation 1000 permettant la définition sur la grille G d'émetteur-récepteur El lk Rlk du centre commun CK du réseau R de couples d'émetteur-récepteur Ei,Ri. Le module d'initialisation 1000 précité permet par exemple l'introduction des coordonnées Ik du centre CK commun précité ou la lecture de celles-ci dans un mode de réalisation qui sera décrit ultérieurement. il permet également l'introduction de la distance de décalage Oi de chaque couple émetteur-récepteur Ei,Ri, ce décalage constituant le diamètre du motif circulaire ou réseau R de centre CK. Le module d'initialisation 1000 permet également l'introduction de la profondeur initiale Z et de l'incidence de l'onde de pression Pi engendrée par chaque émetteur Ei. Le programme de conduite d'acquisition, représenté en figure 7, comporte également un module de calcul 1001 des adresses k, I sur la grille G des émetteurs El Ik-récepteurs Rlk de cette grille G permettant de constituer les couples d'émetteur-récepteur Ei,Ri du réseau R. On comprendra, bien entendu, que le module d'initialisation 1000 et le module 1001 de calcul des adresses permettent de définir sur la grille G un motif initial correspondant au réseau R. A titre d'exemple non limitatif, le module de calcul des adresses ki sur la grille G, module 1001, permet, à partir de la lecture des valeurs initiales entrées ou, respectivement lues en ce qui concerne le centre commun CK par le module 1000, de calculer les adresses correspondantes kl sur la grille G constituant grille d'acquisition en sismique tridimensionnelle. A titre d'exemple non limitatif, les coordonnées du centre commun CK ayant été entrées, ainsi que la distance de décalage Oi, le module de calcul d'adresses peut comporter un sous module de calcul du rayon du motif correspondant au réseau R, ce rayon étant égal à Q , suivi d'un module de calcul de 2 la distance.au centre commun CK de tous les' émetteurs El Ik et les récepteurs Rlk compris dans un domaine délimité par un carré Cl circonscrit au réseau R à un pas de distance Px ou Py correspondant, ce carré étant noté Cl sur la figure 6 et entre un domaine en croix C2 inscrit dans le réseau R, ainsi que représenté sur la figure 6. On comprendra, bien entendu, que le domaine, ainsi délimité par les courbes Cl et CZ, peut très facilement être défini à partir des adresses en lignes I et en colonnes K des émetteurs et des récepteurs El I k Rlk de la grille G. Un calcul de la distance au centre commun CK de chaque émetteur correspondant situé entre les courbes Cl et C2 permet alors de choisir les émetteurs-récepteurs El Ikw Rlk dont la différence de la distance au centre commun CK au rayon du motif R est minimum. L'ensemble de ces émeneurs-récepteurs constitue ainsi le réseau R constitué en motif circulaire. Afin de faciliter les calculs, de façon avantageuse mais non limitative, il sera possible de 'choisir le centre CK comme centre de symétrie d'une maille constituée par deux récepteurs d'une première ligne et deux émetteurs d'une ligne successive. Dans ce cas, il est alors avantageux de choisir un décalage ou distance de décalage entre couples d'émetteur-récepteur Ei,Ri égal à un multiple impair de pas Px ou Py de façon à obtenir un maximum d'émetteur-récepteur formant couples au voisinage du motif constitué par la circonférence. Sur la figure 6, ces émetteurs-récepteurs ont été notés respectivement El, R2, E4, ES, A6 et R1, E2, E3, R4, R5, E6. On comprendra, en effet, que dans ces conditions, du point de vue de la seule symétrie, les émetteurs-récepteurs El,E2, R1,R2, E4,ES, et RS,R4 jouent sensiblement le même rôle par rapport au centre commun CK. Le calcul des adresses correspondantes des émetteurs et des récepteurs ayant été effectué, ainsi que décrit précédemment, il est alors possible d'attribuer aux émetteurs d'une demi-circonférence du réseau R, ainsi déterminé, un numéro d'ordre qui correspond au numéro d'ordre i, les récepteurs d'ordre correspondants associés se trouvant automatiquement symétriquement par rapport au centre CK. Ainsi qu'on le remarquera en outre sur la figure 7, le programme de conduite d'acquisition comporte egalement ensuite un module 1002 de transformation des adresses kl des émetteurs El Ik- récepteurs Rlk, précédemment déterminés pour constituer le réseau R constituant le motif circulaire d'origine, pour choisir une succession de nouveaux motifs obtenus par homothétie du motif initial. Les nouveaux motifs successifs permettent une acquisition variable en profondeur d'un. pas Pf détermine. Ainsi le module 1002 précité permet-il l'entrée du paramètre de pas de profondeur Pf, lequel permet, en fonction de l'incidence choisie précédemment, de déterminer la nouvelle valeur de décalage 01r+I par rapport à la valeur de décalage ir précédente à partir d'un simple calcul trigonométrique et par appel à nouveau du module de calcul des adresses pour déterminer les nouvelles adresses des émetteurs-récepteurs correspondant aux nouveaux motifs, ainsi que représenté sur la figure 5, les adresses kl nouvelles des nouveaux émetteurs-récepteurs constituant le nouveau motif pouvant être obtenues à partir de courbes C'I et C'2 obtenues par homothétie par rapport aux courbes Cl et C2 précitées. Bien entendu, le module de transformation des adresses 1002 est itératif afin d'obtenir le nombre de domaines élémentaires de , de , de espacés selon l'axe des profondeur Z, ainsi que précédemment décrit sur la figure 5. En outre, ainsi qu'on l'a représenté sur la figure 7, les modules d'initialisation 1000 de calcul des adresses 1001 et de transformation des adresses 1002 sont intégrés dans une boucle itérative permettant d'assurer une translation selon les directions x et y de la surface S du milieu matériel. On a représenté les positions successives CK-I, CK, CK+I de centres communs successifs, décalées en translation successivement selon la direction y, puis sur la direction x de deux pas Py, respectivement Px successivement. Il va de soi que tout chemin de déplacement du centre commun CK peut être envisagé pour constituer un balayage approprié de la surface S par l'intermédiaire des émetteurs Ellk, Rlk de la grille G pour obtenir une acquisition dans un domaine D de dimension sensiblement quelconque. Ainsi qu'on le remarquera sur la figure 7, la boucle itérative comprend un module d'initialisation, lequel est noté 1005, permettant d'entrer notamment les coordonnées sur la grille G d'une'suite d'un nombre NK de positions décalées en translation d'un pas Px ou Py ou d'un multiple de ces pas, des centres communs CK de coordonnées Ik. Les positions des centres communs CK peuvent être entrées par leurs coordonnées Ik exprimées par rapport à la grille G. Ce module 1005 permet notamment l'entrée du nombre NK 'de 'positions de centres communs nécessaires à assurer le balayage de la surface S précitée, ainsi que l'initialisation d'un paramètre de comptage K à la valeur 1. La bducle itérative précitée comporte également une boucle de comptage à partir d'une variable de comptage K précitée initialisée précédemment à 1, cette boucle de comptage comportant un test 1003 de la variable de comptage K comparé à la valeur NK entrée précédemment. Une réponse négative à ce test a pour effet d'incrémenter la variable de comptage d'une valeur I et de retourner au module d'initialisa tion d'entrée, lequel permet, par lecture des coprdonnées du centre commun C de valeurs K correspondantes, d'appeler ensuite le module 1000 d'initialisation, lequel permet à nouveau la conduite du processus, ainsi que précédemment décrit, grâce aux modules 1000, 1001 et 1002. Un balayage complet des couchés sousjacentes de la surface S est ainsi effectué. Sur réponse positive au test 1003 précité, un module 1006 permet d'appeler les sous programmes de tri et de valeurs de temps de propagation et d'amplitude, tels que précédemment décrits dans la description, afin d'assurer le traitement de ces valeurs et la détermination des fluctuations de temps de propagation et d'amplitude correspondants pour les différents centres communs CK répartis sur la surface S. On a ainsi décrit un procédé et un système de détermination de l'anisotropie azimutale d'un domaine D d'un milieu matériel particulièrement performants dans la mesure où un tel système et un tel procédé peuvent être utilisés dans des domaines aussi variés que l'étude des matériaux et, de manière très significative, dans le domaine de la sismique tridimensionnelle. On comprendra, en particulier, que dans ce dernier cas le procédé et le système, objets de l'invention, apparaissent particulièrement avantageux dans la mesure où une grille d'acquisition en sismique tridimensionnelle peut facilement être utilisée, celle-ci étant déjà implantée grâce à un logiciel approprié, lequel a fait l'objet de la présente description. Dans le cas précité, selon la fréquence des signaux reçus, le domaine élémentaire de peut présenter une dimension de quelques mètres à quelques dizaines de mètres, la surface de réflexion élémentaire a peut être constituée par un carré ou un rectangle de quelques dizaines de mètres de côté correspondant à une cellule élémentaire désignée en langage anglo-saxon, par "bin" dans le domaine technique de la sismique tridimensionnelle. REVENDICATIONS I. Procédé de détermination de l'anisotropie azimutale dans un domaine (D) d'un milieu matériel présentant au moins une surface de séparation sensiblement plane, caractérisé en ce que ledit procédé consiste pour une acquisition pour un domaine élémentaire de - à disposer sur ladite surface de séparation une pluralité de couples d'émetteur-récepteur (Ei,Ri) d'une onde de pression dirigée vers le milieu matériel à étudier, chaque couple d'émetteur-récepteur étant disposé selon une direction (a i) déterminée par rapport à une direction de référence (N) autour d'un centre commun (CK) et un émetteur (Ei) étant séparé d'un récepteur correspondant (Ri) d'une distance (Oi) de décalage, - à engendrer à partir d'un émetteur (Ei) une onde de pression (Pi) et à enregistrer au niveau du récepteur (Ri) qui lui est associé les échos ou ondes de pression réfléchies dans une fenêtre de temps de propagation définissant ledit domaine élémentaire (de), - à répéter l'étape précédente pour chacun desdits couples d'émetteur récepteur (Ei,Ri) pour le domaine élémentaire (de) considéré, - à déterminer pour chaque écho reçu les fluctuations d'amplitude et de temps de propagation pour établir une loi de variation de propagation desdites ondes de pression réfléchies ou échos en fonction de l'azimut correspondant. 2. Procédé selon la revendication 1, caractérisé en ce que l'étape de détermination de la loi de variation de propagation desdites ondes de pression réfléchies ou échos en fonction de l'azimut correspondant consiste - à mémoriser les valeurs de temps de propagation (Ti) pour chaque valeur d'azimut (a i), - à comparer pour une valeur d'azimut (ai) donnée la valeur du temps de propagation (Ti) correspondant à la valeur moyenne des temps de propagation (T) au point considéré. 3. Procédé selon l'une des revendications I ou 2, caractérisé en ce que l'étape de détermination de la loi de variation de propagation desdites ondes de pression réfléchies ou échos en fonction de l'azimut correspondant consiste - à mémoriser les valeurs d'amplitude (Ai) des ondes de pression réfléchies ou écHos pour chaque valeur d'azimut ( a i), - à comparer pour une valeur d'azimut (ai) donnée la valeur de l'amplitude (Ai) correspondante à la valeur moyenne des amplitudes (A) au point considéré. 4. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'émission et la réception des ondes de pression réfléchies ou échos sont effectuées pour une plage d'azimuts ( ai) comprise entre 0 et 1800 à incidence constante, lesdits couples d'émetteurrécepteur (Ei,Ri) présentant une même distance de décalage (Oi) et étant disposés en un motif formant une même circonférence par rapport audit centre commun (CK). 3. Procédé selon la revendication 4, caractérisé en ce que, en vue d'effectuer une mesure en une succession de domaines élémentaires (de) de profondeur successive incrémentée, ledit procédé consiste à effectuer une pluralité d'acquisitions, ledit motif, formant une même circonférence par rapport audit centre commun (CK), étant transformé par une suite d'homothéties de centre (CK) et de rapport (Kr) telles que Oir = K r x Oi,,l où Oir et ir I représentent les distances de décalage des couples émetteur-récepteur (Ei,Ri) pour un incrément de profondeur (Pf) correspondant à un pas de décalage en profondeur du domaine élémentaire (de) considéré. 6. Procédé selon la revendication 4 ou 5, caractérisé en ce que, en vue d'effectuer une acquisiton en une succession de domaines élémentaires (de) décalés en position selon une première (x) ou deuxième (y) direction parallèle à la surface de séparation, ledit motif formant une meme circonférence par rapport audit centre commun (CK), ainsi que ce dernier, sont décalés en translation, sur ladite surface d'un incrément (Px, Py) respectivement. 7. Procédé selon l'une des revendications précédentes, caractérisé en ce que la répétition pour chaque couple d'émetteur-récepteur (Ei,Ri) est effectuée par une émission-enregistrement séquentielle des ondes de pression réfléchies pour les directions d'azimut (a i) considérées. 8. Procédé selon l'une des revendications 1 à 6, caractérisé en ce que la répétition pour chaque couple d'émetteur-récepteur (Ei,Ri) est effectuée par une émission-enregistrement simultanée des ondes de pression réfléchies, l'émission de l'onde de pression (Pi) par chacun des émetteurs (Ei) consistant en une onde de pression périodique codée, chaque émetteur (Ei) émettant une onde de pression (Pi) codée selon un code qui lui est propre. 9. Procédé selon la revendication 8, caractérisé en ce que le codage consiste en un codage de phase, ledit codage étant représentatif du numéro d'ordre (i) de l'azimut ( a i) ou angle d'orientation correspondant. 10. Système de détermination de l'anisotropie azimutale dans un domaine (D) constitué de domaines élémentaires (de) d'un milieu matériel présentant au moins une surface de séparation sensiblement plante, caractérisé en ce qu'il comporte - un réseau de couples d'émetteur-récepteur (Ei,Ri) disposés selon une pluralité de directions d'orientation ( a i) par rapport à une direction de référence (N), autour d'un centre commun (CK), - des moyens (I) de commande desdits couples d'émetteur-récepteur (Ei,Ri) interconnectés auxdits couples d'émetteur-récepteur pour assurer la commande en émission des émetteurs (Ei) et en réception des récepteurs (Ri) dans une fenêtre de temps de propagation définissant ledit domaine élémentaire (de), - des moyens (100) de calcul des fluctuations d'amplitude et de temps de propagation des échos reçus en fonction de la direction d'orientation ( a i). Il. Système selon la revendication 10, caractérisé en ce que ledit réseau (R) de couples d'émetteur-récepteur (Ei,Ri) est formé par une grille (G) de lignes d'émetteurs (El Ik) et de récepteurs (Rlk), espacés d'un pas (Px), deux lignes successives étant formées par un alignement d'émetteurs respectivement de récepteurs et espacées d'uh même pas (Py), égal à (Px), de façon à constituer dans une direction perpendiculaire auxdites lignes un arrangement en colonnes similaire d'émetteurs (El Ik) récepteurs (Rlk) alternés, chaque couple d'émetteur-récepteur (Ei,Ri) étant choisi par adressage de couples (El.lk), (Rlk) d'émetteur-récepteur pour constituer un réseau (R) selon un motif formant sensiblement une circonférence. 12. Système selon la revendication 10 ou Il, caractérisé en ce que lesdits moyens de commande desdits couples d'émetteur-récepteur (Ei,Ri) interconnectés auxdits couples d'émetteur-récepteur et lesdits moyens de calcul sont . constitués par - un circuit d'interface (200) assurant l'interconnexion des. lignes d'émetteurs et de récepteurs et permettant d'assurer le déclenchement de l'émission des ondes de pression (Pi) par adressage des émetteurs (Ei) considérés, . la réception et l'enregistrement, par adressage des récepteurs (Ri) considérés, des ondes de pression réfléchies ou échos dans une fenêtre temporelle définissant le domaine élémentaire (de) considéré, . la conversion analogique numérique des signaux d'échos et la mémorisa tion des signaux d'échos numérisés obtenus, - un calculateur (100) de conduite d'acquisition et de traitement permettant la commande par le circuit d'interface (i00) des émetteurs (Ei) et récepteurs (Ri), et le traitement des échos enregistrés pour déterminer, pour un centre (CK) considéré, la direction des axes lent et rapide par rapport à la direction de référence et l'atténuation différentielle d'amplitude selon les axes lent et rapide. 13. Système selon la revendication 12, caractérisé en ce que, afin d'effectuer le traitement des échos enregistrés, ledit calculateur de conduite d'acquisiton comprend - un sous programme de tri des valeurs des temps de propagation permettant de déterminer les azimuts (a I) et (a r) correspondant, d'une part, au temps de propagation maximum et, d'autre part, au temps de propagation minimum, lesdites directions ou azimuts ( a 1), ( a r) définissant respectivement les directions des axes lent et rapide, - un sous programme de tri des valeurs des amplitudes et de calcul de l'amplitude crête à crête de la variation d'amplitude des échos reçus en fonction de l'azimut ( a i) et de l'atténuation différentielle selon lesdits axes lent et rapide. 14. Système selon l'une des revendications précédentes 10 à 13, caractérisé en ce que, en vue d'effectuer une acquisition sur un domaine (D) constitué par un ensemble de domaines élémentaires (de), ledit calculateur de conduite d'acquisition comporte un programme de conduite d'acquisition, ledit programme de conduite d'acquisition comprenant - un module d'initialisation (1000) permettant la définition, sur la grille (G) d'émetteur-récepteur (El Ik) (Rlk) du centre commun (CK) du réseau de couples d'émetteur-récepteur (Ei,Ri) de la distance de décalage (Oi) diamètre du motif circulaire de centre (CK), de la profondeur initiale (Z) et de l'incidence de l'onde de pression (Pi) engendrée par chaque émetteur (Ei), - un module de calcul (1001) des adresses (k, I) sur la grille (G) des émetteurs (E, Ik) récepteurs (Rlk) permettant de constituer les couples d'émetteur-récepteur (Ei,Ri), lesdits modules d'initialisation et de calcul des adresses permettant de définir sur la grille (G) un motif initial, - un module de transformation desdites adresses (k, I) des émetteurs (Ei,1 k) récepteurs (Rlk) pour choisir une succession.de nouveaux motifs obtenus par homothétie du motif initial, lesdits nouveaux motifs successifs permettant une acquisition variable en profondeur d'un pas (Pf) déterminé. 15. Système selon la revendication 14, caractérisé en ce que lesdits modules d'initialisation, de calcul des adresses, et de transformation des adresses sont intégrés dans une boucle itérative permettant d'assurer une translation selon les directions (x) et (y) de la surface (S) du milieu matériel d'un pas ou d'un nombre de pas déterminé de la position successive du centre commun (CK-I, CK, CK+I) du réseau (R) constituant le motif initial, de façon à assurer une acquisition pour l'ensemble des domaines élémentaires (de) constituant le domaine (D). 16. Système selon la revendication 15, caractérisé en ce que ladite boucle itérative comprend - un module d'initialisation (1005) permettant d'entrer les coordonnées sur la grille (G) d'une suite d'un nombre (NK) de positions décalées en translation d'un pas (Px), (Py) ou d'un multiple de ces pas, - une boucle de comptage dans laquelle lesdits modules sont intégrés.