Brevet FR2716006A1 — publié le 19950811

i La présente invention se rapporte à la détection et à la cartographie de réservoirs d'hydrocarbures en utilisant des techniques radar. Jusqu'à présent, en général, deux types prédomi- nants de techniques de détection et de cartographie de réservoirs d'hydrocarbures ont existé qui utilisaient l'énergie des ondes. Une de ces techniques, et actuellement la plus largement utilisée, était basée sur les ondes d'énergie acoustique ou sismique, soit pour des relevés par réflexion, soit pour des relevés par réfraction. L'exploration sismique est généralement onéreuse et prend beaucoup de temps, nécessitant de nombreux arrangements de géophones de détection ou de transducteurs devant être placés sur des surfaces considérables. L'autre technique de détection par énergie d'onde destinée aux réservoirs d'hydrocarbures est basée sur l'utilisation d'énergie électromagnétique ou ondes à radar. Une masse considérable de travaux de l'art antérieur existe et est basée sur l'hypothèse que la présence de réservoirs d'hydrocarbures souterrains provoque la présence identifiable discernable de gaz s'échappant dans l'atmosphère à la surface de la terre. Une réflexion soutenue et de nombreux efforts se sont appliqués à détecter et à identifier ces gaz s'échappant. Il était supposé que la présence et le type de gaz s'échappant de la terre pouvaient être détectés à partir des caractéristiques moléculaires uniques, identifiables par radar, des différentes espèces de gaz. Des exemples dans ce domaine comprennent les brevets des Etats-Unis d'Amérique n 3 651 395; 4 100 481; 4 132 943; 4 385 516; et 4 504 833. Brièvement, la présente invention réalise un procédé nouveau et amélioré de détection et de cartographie de réservoirs d'hydrocarbures à l'aide d'ondes radar. Grâce au procédé de la présente invention, une modification géochimique proche de la surface de la roche et du sol est explorée à l'aide d'énergie électromagnétique afin de détecter la présence de réservoirs d'hydrocarbures souterrains, sous la surface de la terre. L'onde électromagnétique ou l'énergie radar est transmise conformément à la présente invention dans la surface de la terre afin de pénétrer les matériaux proches de la surface de la terre. L'énergie radar est adaptée pour pénétrer toute accumulation de carbone, silice ou équivalents dans les matériaux proches de la terre. A l'aide de la présente invention, la présence de ces carbonates, silices ou équivalents a été considérée comme résultant de modifications géochimiques à, ou proches de, la surface de la terre. Ces modifications géochimiques surviennent en raison des microfiltrations des hydrocarbures dissous dans l'eau. La source d'hydrocarbures généralement concernée peut être un réservoir de pétrole ou de gaz enterré à plusieurs centaines de mètres sous terre. Une zone anomale de résistivité électromagnétique est formée dans les matériaux proches de la surface de la terre en raison de ces modifications géochimiques. L'énergie d'onde électromagnétique qui pénètre la terre, à proximité de la surface, est déviée vers le bas par la zone anomale de résistivité puis est réfléchie en vue de sa détection par un récepteur d'énergie électromagnétique ou un récepteur radar. Les ondes d'énergie radar réfléchies sont ensuite collectées en tant que données d'analyse afin d'indiquer la présence de réservoirs d'hydrocarbures souterrains. D'autres objets, caractéristiques et avantages de l'invention ressortent de la description qui suit, donnée à titre d'exemple non limitatif, en regard des dessins annexés sur lesquels: - les figures 1 et 2 sont des vues en coupe de caractéristiques géologiques et de résistivité de la surface de la terre au voisinage d'un réservoir d'hydrocarbures souterrain; - les figures 3A et 3B sont des diagrammes de forme d'onde de réponses radar possibles correspondant à des caractéristiques géologiques souterraines telles que celles montrées sur les figures 1 et 2; - la figure 4 est un diagramme schématique d'une onde électromagnétique se déplaçant à travers les caractéristiques géologiques telles que celles des figures 1 et 2; - la figure 5 est un diagramme schématique d'une géométrie de transmetteur/récepteur radar pour les opérations de relevé selon la présente invention; - la figure 6 est un bloc-diagramme schématique simplifié d'un système transmetteur/récepteur radar selon la présente invention; - la figure 7 est un tracé du paramètre y transversal radar en fonction de l'angle 0 de dépression ou de brillance destiné aux ondes radar transmises en direction des surfaces de la terre; - les figures 8 et 9 sont des tracés de coefficients de réflexion radar de la terre qui sont fonction de l'angle de brillance ou de pénétration d'une antenne polarisée verticalement et horizontalement, respectivement; - la figure 10 est un diagramme schématique d'une géométrie de transmetteur/récepteur radar alternative selon la présente invention par rapport à celle de la figure 5. Selon la présente invention, il s'est avéré que l'énergie électromagnétique ou les ondes radar pouvaient être utilisées pour détecter et cartographier une modification de rocher et de sol géochimiques à proximité de la surface. Un système de transmetteur et récepteur radar est utilisé pour transmettre et recevoir les ondes radar qui pénètrent l'accumulation de carbonate et de silice à la surface de la terre au-dessus de la microfiltration d'hydrocarbures dissous dans l'eau. La source d'hydrocarbures concernée est typiquement un réservoir de pétrole ou de gaz enterré à plusieurs centaines de mètres, bien que dans certains cas, il puisse s'agir d'un réservoir proche de la surface, ayant des caractéristiques artificielles. L'épaisseur et la résistivité de l'accumulation sont liées à la taille du réservoir. Un écho radar provenant du fond de l'accumulation offre un outil de détection et de cartographie de réservoir selon la présente invention. Dans les dessins, un réservoir d'hydrocarbures R enterré à une certaine distance par exemple, à plusieurs centaines de mètres, dans un grand bassin B sédimentaire sous la surface de la terre 11 est monté (figure 1). Les rocs du bassin B contiennent des microfractures. L'eau provenant des parties les plus profondes du bassin B continue à progresser vers le haut comme indiqué par les lignes de mouvement 14 s'étendant comme indiqué par les flèches 16. Le mouvement de l'eau vers le haut à travers les microfractures est dû à l'action de pompage naturelle des changements diurnes de l'atmosphère de la terre. Généralement, cette eau possède des valeurs Eh et Ph qui la place du côté réducteur d'une échelle de réaction d'oxydo-réduction. L'eau qui passe par le réservoir d'hydrocarbures R sur son trajet la conduisant à la surface 11 devient plus réductrice en raison des gaz d'hydrocarbures dissous en elle. L'eau se déplaçant vers le haut voit sa nature réductrice se modifier au point qu'elle provoque une modification supplémentaire de la roche par rapport aux lieux terrestres proches. Il en résulte un panache ou cheminée 18 réducteur isolé qui finalement progresse vers l'eau souterraine oxydante dans la zone proche de la surface indiquée par 20 et provoque une réaction chimique réductrice sur des matériaux de roche normalement oxydés. A une profondeur comprise entre plusieurs dizaines de mètres et plusieurs centaines de mètres, le résultat de cette réaction chimique réductrice est la formation d'argile. Depuis la surface 11 jusqu'à quelques milliers de mètres de profondeur, le carbonate de calcium et la silice sont précipités comme montré en 22 dans les pores de la roche du matériau proche de surface 20. En passant au-dessus de cette zone et en regardant cette caractéristique géologique à partir d'un avion, on verrait une forme circulaire au sein de laquelle existeraient des états différents de réactions d'oxydo-réduction, selon la quantité d'hydrocarbures dissous dans l'eau. Dans un certain état d'oxydation, habituellement proche de la surface, le carbonate et la silice précipitent et remplissent les espaces de pores afin de former diverses quantités de matériau plus dur indiqué en 22, possédant l'apparence visuelle du calcaire. Sur la figure 2, les caractéristiques géologiques identiques sont montrées, mais la résistivité de la roche est indiquée. La strate dure 22 possède une résistivité grande ou élevée qui parfois dépasse 100 ohmsmètre. Directement sous la strate 22, se trouve une zone de plus basse résistivité en raison de la formation d'argile indiquée par 20. Habituellement, la résistivité moyenne intermédiaire du bassin B est de l'ordre de deux fois supérieure à celle de la zone d'argile 18, mais est bien inférieure à celle de la zone de la strate dure 22. A l'aide de la présente invention, un système radar, (figure 5), habituellement dans la bande X des fréquences (c'est-à-dire 8000 à 12 500 MHz) est utilisé de façon qu'une réponse radar anomale puisse être détectée en provenance de la zone de strate dure 22. La géométrie du système radar est sélectionnée de façon que l'épaisseur de la strate 22 et l'intensité du contraste de résistivité dans la zone d'argile 20 donnent des informations pouvant être utilisées pour déterminer la position d'un trou de forage en vue de tester le réservoir R en cas d'hydrocarbures éventuels. Les figures 3A et 3B sont des diagrammes de formes d'ondes de réponses de réflexion radar possibles, selon la présente invention. La figure 3A représente une réponse dans la zone de la strate 22 à forte résistivité, tandis que la figure 3B représente une réflexion radar dans une zone o il n'y a pas de strate 22 à haute résistivité présente. Comme cela est visible sur la figure 3A sur la strate 22, on trouve en plus de l'amplitude A, de réflexion radar ordinaire une amplitude supplémentaire A2 plus tardive pour le retour vers le récepteur radar. Cette amplitude supplémentaire A2 est une anomalie de réflexion. La hauteur de cette anomalie de réflexion et la zone placée sous sa courbe A2 sont liées directement au contraste de résistivité autour et sous la strate de résistivité 22. L'anomalie A2 supplémentaire lorsqu'elle est ajoutée au retour de la réponse de fond normal Al, apparaît dans le récepteur radar comme si une phosphorescence existait dans l'air situé au-dessus de la strate à forme résistivité 11. Plus la réponse indiquée par l'amplitude A2 est forte, plus la microfiltration d'hydrocarbures provenant du réservoir R est grande. Comme cela est visible sur la figure 4, une impulsion radar incidente se déplaçant le long d'un trajet comme indiqué par une flèche 24 et entrant dans la strate à forte résistivité 22 sous la surface 11 de la terre selon un angle 0 de brillance ou de dépression, pénètre en réalité dans la zone à forte résistivité 22, comme indiqué par une flèche 26. L'énergie radar se déplace vers le bas à travers la strate à forte résistivité 22, une partie d'une telle énergie étant réfléchie en tant qu'impulsion de rétrodiffusion 28 provenant de l'argile 20 et retournant au transmetteur et récepteur X de radar (figure 5). A l'aide de la présente invention, la détection de l'impulsion de rétrodiffusion est possible lorsque l'angle de brillance 0 à la surface 11 est petit comparé à l'angle de dépression existant dans la zone à forte résistivité 22. Ceci rend la section transversale radar, a, en retrait de l'interface 30 entre la strate 22 et l'argile 20, bien supérieure à la section radar, ac, de défaut des échos parasites à la surface (a >> ac). Même lorsque l'angle de brillance > sur la strate 22 est petit, l'impulsion réfractée est presque normale à l'interface argile/strate 30. Ceci se traduit par une augmentation significative ou poussée de la section transversale radar de l'interface a comparée à la surface de la section transversale radar ac. L'importance d'une telle augmentation détermine à quelle profondeur le radar peut effectivement pénétrer ou "voir" dans la couche souterraine de la terre, lorsque l'écho radar est diminué d'environ 18 dB à chaque interface ou transition avec un matériau de résistivité différente. Egalement, il existe un angle de brillance optimal et une polarisation selon laquelle ces échos radar peuvent être observés. Cet angle de brillance s'est avéré être d'environ 19 degrés pour une constante diélectrique E = 8E0, la polarisation étant verticale. De plus, une comparaison entre les échos polarisés verticalement et horizontalement peut être utilisée pour aider à fournir une indication sur ces positions o la pénétration survient. Comme exemple, on considère un système radar X en bande X conventionnel muni d'une antenne 42 (figure 5) montée sur un véhicule V et possédant les caractéristiques suivantes: Puissance de pointe du transmetteur = 5 kW Longueur de l'impulsion rayonnée = 0,1016 gsec = T Largeur de bande passante du radar = 1/T = Bn Diamètre de l'antenne 42 = 1,22 m = D Fréquence de fonctionnement = 010 Hz Hauteur de l'antenne 42 au-dessus d'une surface plate terrestre 11 = 3,04 m Distance à la cible = 304 m Il est préféré que le système radar fonctionne en bande X (8000-12 500 MHz). Il doit être compris, cependant, que des niveaux de fréquence légèrement inférieurs peuvent également être utilisés, le facteur limitant étant le diamètre ou la taille de l'antenne 42 pouvant être éventuellement transportable. NIVEAU DE BRUIT D'ENTREE Un diagramme schématique de système récepteur de radar simplifié, est montré sur la figure 6, afin de permettre le calcul du niveau de bruit d'entrée équivalent destiné au système radar X. L'antenne 42 transmet l'énergie d'onde radar et détecte ou collecte en tant que données l'énergie d'onde électromagnétique réfléchie provenant des zones souterraines telles que l'interface 30. Chaque élément récepteur est supposé posséder la même bande passante de bruit, Bn, égale à la bande passante du récepteur de 3 dB. Seul un transducteur actif unique 40, ou amplificateur, est montré bien que les récepteurs possèdent souvent plusieurs étages amplificateurs. Cependant, si le gain du premier amplificateur est grand, comme c'est typiquement le cas, les étages successifs ayant habituellement peu d'effet lors de la détermination de l'énergie du bruit d'entrée équivalent. La température de l'antenne, Ta, représente l'énergie du bruit de sortie disponible, Pa, aux bornes de l'antenne réceptrice selon la relation: Pa = K Ta Bn Ta dépend des sources de rayonnement diverses (selon la loi de rayonnement de Planck) dans la configuration de réception de l'antenne 42 qui comprend des lobes latéraux et des lobes arrière. La ligne de transmission 44 entre l'antenne 42 et l'amplificateur 40 possède une perte utile Lr et une température de bruit efficace Tri tandis que l'amplificateur 40 possède un gain utile GA et une température de bruit efficace TA. Le bruit d'entrée équivalent, N, peut être montré, qui doit satisfaire la relation: N - bruit d'entrée équivalent = K Ts B. K étant la constante de Boltzmann et Ts, la température du système, est approximativement donnée par l'expression: T, = Ta + Tr + Lr + TA Pour le présent exemple, les valeurs du cas le meilleur ou du pire cas sont supposées être la température d'antenne Ta. Lorsqu'on considère une antenne 42 au, ou proche du, sol 11, Ta peut être considéré être typiquement d'environ 300 K. La ligne de transmission 44 possède une température de bruit d'entrée efficace qui est donnée par l'équation: Tr = Ttr (Lr- 1) dans laquelle Ttr est la température thermique de la ligne de transmission 44, typiquement de 290 K. Il serait raisonnable de s'attendre à ce que Lr - 2 dB. TA étant donné par la relation suivante: TA = To (Fn - 1) dans laquelle To = 290 K et Fn étant le chiffre du bruit, typiquement environ 10 dB pour une conception simple, comme dans cet exemple. La température du système est d'environ Ts = 4607 K Cela conduit à la valeur suivante pour le niveau de bruit d'entrée équivalent N: N = 62575 x 1013 watts. MESURE DE NIVEAU DES ECHOS PARASITES L'étape suivante consiste à déterminer le niveau des échos parasites. Le niveau d'écho des échos parasites, C, est donné par l'équation suivante: C PtG2V2oC (4n) 3R4L dans laquelle Pt = puissance délivrée à l'antenne 42 = puissance du transmetteur/2 G = gain de l'antenne A = c/f = 3 x 10-2m Gamme R = 304,8 m L = 3 dB (blocage d'ouverture) ac = section transversale de radar des échos parasites Il est visible qu'il existe un besoin de connaître les valeurs G et ac avant que la puissance des échos parasites puisse être calculée. Pour y parvenir, on doit d'abord estimer la largeur du faisceau de 3 dB, 03dB- 63dB's 1.01 = 0,0249 rad D A partir de cela, il est possible d'estimer la directivité GD le long de l'axe de l'antenne 42 GD - 4 3dB f3dB On s'attendrait à ce que l'antenne ait une perte très petite en raison du chauffage ohmique, c'est-à- dire G = GD = 43 dB. La portée de la gamme du défaut des échos parasites est c,/2 1,52 m et la portée latérale du défaut des échos parasites est d'environ: 1000 x 03dB = 7,62 m. La surface du défaut des échos parasites peut ainsi être considérée dans cet exemple comme étant environ de 115,66 m2. La figure 7 est un diagramme provenant d'un livre de référence connu, "Introduction to Radar Systems", Skolnick, p. 523. Elle est utilisée pour estimer la section transversale radar y en bande X en supposant une voie accidentée. Une telle surface serait approximativement semblable à la courbe descendante d'un terrain ayant très peu de végétation. La section transversale y est également définie par la réflectivité de surface sans dimensions a' et l'angle d'incidence y, comme suit: GO Y= sin4 A partir de la figure 7, une valeur de y = -36 dB = a /sin (angle de brillance y), avec a' étant la section transversale des échos parasites par unité de surface, est utilisée. Dans le cas présent, ao est: a = Y sin 0,573 = 2,512 x 10-6 et ac est donné par: ac = aO [RO3B] [C,/2] = 2,9 x 10-4 m2 - -35,37 dB En utilisant les formules ci-dessus, la puissance des échos parasites C peut être trouvée comme étant égale à : C = 7,6 x 10-9 watts A partir de cela, il est clairement visible que la puissance des échos parasites C domine le bruit du récepteur N. De ce fait, c'est le rapport signal-écho parasite qui détermine la performance du radar dans des intervalles de l'ordre de 300 mètres. SIGNAL PROVENANT DE L'INTERFACE Le niveau de puissance du signal S provenant de l'interface 30 est donné par l'équation suivante: Pt G2A2 (4I) 3R'L dans laquelle a est la section transversale à radar d'une interface 30 sous la surface 11. Dans cet exemple, une hypothèse est faite qui consiste à obtenir cette valeur de a grâce à une formule de section transversale radar de plaque plate modifiée: 4nA2 cos2O(Rc) 2 dans laquelle A est la surface identique à celle du défaut des échos parasites, 6 étant l'angle de brillance de pénétration indiqué par la flèche 26 (figure 4) par rapport à la normale, et RC étant le coefficient de réflexion depuis l'interface d'argile 30. Lorsque l'argile 20 est fortement conductrice, Rc - 1. La perte L doit inclure la perte d'épaisseur de peau ainsi que les pertes de transmission à travers l'interface air-terre puis en retour, à travers l'interface terre-air. Dans tout ce qui suit, on a fait l'hypothèse que E = 8 E0, g = g0 et que la conductivité se modifie selon quatre valeurs 1/10, 1/100, 1/250 et 1/5000. L'épaisseur de peau de ces valeurs de conductivité est donnée, comme suit: Conductivité (Äm)-' Epaisseur de peau (mètres) 0,1 0,15 m 0,01 1,50 m 0,004 3,75 m 0,002 7,50 m Il existe une perte à deux voies de 17,37 dB associée à une pénétration de l'épaisseur de peau. L'angle de rayon 28 sous la surface, est de 20,7 par rapport à la verticale. Il n'est pas exactement vertical car les conductivités sont trop faibles. A partir de cette valeur d'angle de rayon, on peut alors calculer: a = 47,54 dB, Rc = 0,2 étant le coefficient de réflexion supposé provenant de l'interface d'argile 30. Cette valeur de a est une valeur très grande par rapport à ac. Les échos radar provenant de la surface 10 sont observés selon des angles de faible brillance 0, mais la surface souterraine 30 est observée selon des angles presque verticaux. Ceci s'est avéré conduire à une différence frappante entre les valeurs. Pour être complet, on doit encore connaître les pertes de transmission allant vers le bas et allant vers le haut à travers l'interface air-terre 11. La figure 8 et la figure 9 montrent le coefficient d'une réflexion provenant de la terre pour une antenne polarisée verticalement et horizontalement, respectivement. Pour être simple, l'analyse suivante suppose une polarisation verticale. A partir de la figure 8, on peut voir que la perte à travers l'interface air-terre 11 est d'environ de 26 dB. La perte à travers l'interface terre-air s'est avérée être de 20,9 dB, en utilisant les mêmes techniques de calcul mentionnées ci-dessus. Le résultat est une perte totale de transmission de 46,9 dB. Si l'on suppose que le coefficient de transmission est de un diminué du coefficient de réflexion, il est possible de calculer le signal par: 3,037 x 10-5 Lskin dans laquelle Lskin, la perte de peau, est conservée en tant que paramètre. Le rapport signal-écho parasite S/C est donné par S/C = 36 dB Lskin(dB) Deux profondeurs de peau correspondraient environ au point o S/C = 1, et il s'agirait grossièrement d'une profondeur pour laquelle on pourrait raisonnablement attendre que l'énergie d'onde radar pénètre la géométrie supposée décrite. Une géométrie alternative est montrée sur la figure 10, la gamme oblique R étant maintenue à 304 m pour les besoins de l'analyse. Il est à noter à partir de la figure 8, que toute l'énergie radar pénètre dans le sol lorsque l'angle de brillance 0 est d'environ 19 . Dans ces circonstances, il n'y a pas d'énergie de perdue lors de la pénétration dans la terre, l'angle de réfraction étant d'environ 6,61 , augmentant ainsi a de 0,52 dB. Seule une perte de 5,3 dB est rencontrée lors du retour à travers l'interface terre-air 11. Cela donne une capacité de pénétration supplémentaire de 41,6 dB. Naturellement, le niveau des échos parasites augmente lorsque l'angle de brillance 0 est accru jusqu'à 19 . En revenant à nouveau à la figure 7 concernant le paramètre de section transversale radar y, il est visible que y s'accroît d'environ 10 dB. Egalement, le défaut des échos parasites augmente d'environ 0,25 dB. De ce fait, S/C total augmente d'environ 32 dB, ce qui est presque assez pour voir quatre épaisseurs de peau dans la terre. Il est à noter que lorsque la couche d'argile 20 est fortement conductrice, on peut obtenir 14 dB supplémentaires (Rc 1), le radar pénètrant facilement ou "voyant" jusqu'à quatre profondeurs de peau. Même si l'angle de brillance n'est pas exactement de 19 mais possède une valeur de l'ordre de 5 , on anticiperait seulement environ 3 dB de perte lors de la pénétration de l'interface air-terre. A partir de ce qui précède, il est visible que le système X transmetteur/récepteur radar transmet et reçoit les ondes radar qui pénètrent l'accumulation de carbonate et de silice 22 à la surface de la terre 11 vers l'interface 30 au-dessus des microfiltrations d'hydrocarbures dissous dans l'eau, et se déplaçent vers le haut à partir du réservoir R. A l'aide des techniques précédentes, il est possible d'utiliser l'énergie électromagnétique ou les ondes radar pour détecter et cartographier la modification géochimique proche de la surface de la roche et du sol. L'épaisseur et la résistivité de l'accumulation de la couche 22 sont liées à la taille du réservoir souterrain R. En outre, l'écho radar provenant de l'interface inférieure 30 fournit un outil de détection de réservoir et de cartographie. La divulgation et la description précédentes de l'invention ont pour but d'illustrer et d'expliquer cette invention, divers changements quant à la taille, la forme, les matériaux, les éléments, les composants de circuit, les connexions de câbles et les contacts, et aussi relatifs aux détails des circuits illustrés et à la construction et au procédé du fonctionnement pouvant être faits sans sortir du cadre de l'invention. REVENDICATIONS 1. Procédé d'exploration à l'aide d'énergie d'onde électromagnétiques afin de détecter la présence de réservoirs d'hydrocarbures souterrains (R) sous la surface de la terre (11), caractérisé par les étapes consistant a: transmettre l'énergie d'onde électromagnétique dans la surface de la terre afin de pénétrer des zones de résistivité anomale dans les matériaux terrestres proches de la surface (20) ; et collecter les données indiquant la présence de réflexions d'énergie d'onde électromagnétique provenant des zones de résistivité anomale en vue de l'analyse pour indiquer la présence de réservoirs d'hydrocarbures souterrains (R). 2. Procédé selon la revendication 1, caractérisé en ce que les zones de résistivité anomale sont provoquées par des modifications géochimiques des matériaux terrestres proches de la surface (20) en raison de la microfiltration de l'eau entraînée par les gaz d'hydrocarbures. 3. Procédé selon la revendication 2, comprenant en outre l'étape consistant à : dévier l'énergie d'onde électromagnétique vers le bas à travers les matériaux terrestres proches de la surface (20) lorsque l'énergie entre dans les matériaux terrestres proches de la surface (20). 4. Procédé selon la revendication 2, caractérisé en ce que l'eau contenant les gaz d'hydrocarbures entraînés provoque une action réductrice sur les matériaux terrestres proches de la surface (20), formant une argile. 5. Procédé selon la revendication 4, caractérisé en ce que l'argile formée est un matériau tel que du carbonate de calcium ou de la silice, qui est précipité dans les pores des roches et du sol, dans les matériaux terrestres proches de la surface (20). 6. Procédé selon la revendication 5, caractérisé en ce que le matériau précipité à la surface proche de la surface terrestre forme une strate (22) de matériau à forte résistivité proche de la surface terrestre, audessus d'une couche (20) de matériau à faible résistivité composé d'argile. 7. Procédé selon la revendication 6, comprenant en outre l'étape consistant à positionner une antenne de transmission d'onde d'énergie électromagnétique à une hauteur et une distance à partir d'une zone examinée, de sorte que l'énergie arrive à la surface de la terre selon un angle de brillance (0) qui est faible comparé à l'angle d'incidence au sein de la strate (22). 8. Procédé selon la revendication 1, caractérisé en ce que l'énergie d'onde électromagnétique transmise durant ladite étape de transmission se trouve dans la bande en X ou légèrement dessous.