FR 3064370 A1 20180928 FR 1752521 A 20170327 1 L'invention concerne le domaine des procédés géophysiques d'exploration et de caractérisation du sous-sol. En particulier, la présente invention concerne un procédé pour construire une image du sous-sol, représentative de la résistivité électrique du sous-sol et/ou de la polarisation provoquée du sous-sol. En effet, les caractéristiques du sous-sol peuvent-être appréhendées par l'étude des propriétés électriques, à savoir la résistivité électrique et la chargeabilité des formations en présence. La chargeabilité est mesurée simultanément avec la résistivité, sur un même signal d'injection. De telles images de propriétés électriques sont utilisées dans de nombreux procédés 10 industriels nécessitant la mise en évidence d'anomalies de résistivité ou de chargeabilité présentes dans le sous-sol, tels que la recherche d'eau souterraine, la recherche de cavités, la géothermie, ou encore l'exploration minière. On appelle « résistivité électrique» d'un matériau, sa capacité à s'opposer à la circulation du courant électrique. 15 On appelle « polarisation provoquée » du sous-sol le phénomène qui consiste en la polarisation de ses matériaux constitutifs sous l'effet d'un courant induit. La « chargeabilité » est la mesure de la quantité d'énergie potentielle électrique stockée par le matériau durant sa polarisation. Cette énergie est restituée au milieu lors de la coupure de l'injection et se traduit par un retour progressif à zéro du signal de réception selon une 20 courbe dite de décharge. Ces propriétés sont fonction, entre autres, de la nature de la roche, de sa structure, de la perméabilité, du diamètre des pores, de la conductivité des fluides d'imbibition et de son contenu en sulfures métalliques. Pour réaliser de telles images, un procédé connu consiste à utiliser un appareil qui 25 introduit un courant continu ou alternatif à très basse fréquence dans le sol par l'intermédiaire de deux électrodes « d'injection » (A, B) et mesure la différence de potentiel entre deux autres électrodes « de potentiel » (M, N). Ce type de procédé utilise des équipements de mesures du type « centralisé », c'est-à-dire qu'il existe une même centrale d'émission et de mesures (centrale d'acquisition) 30 raccordée à des électrodes qui peuvent être à la fois des électrodes d'injection (A, B) et de réception (P1, P2, Pn). Pour des applications à faible profondeur (dédiées à la recherche en eau ou à la surveillance de l'environnement), les fonctions d'émetteur et de récepteur sont assurées par une même centrale. 35 En revanche, pour des applications nécessitant des profondeurs d'investigation plus grandes, et donc des puissances d'injection plus élevées, les fonctions d'émetteur et de récepteur sont assurées par deux instruments distincts. 3064370 2 Les équipements de mesures de l'état de l'art posent des contraintes techniques de mise en oeuvre de la mesure, principalement lors de la mise en place des câbles. L'investigation de grands volumes demande habituellement des espacements maximaux entre les électrodes pouvant aller jusqu'au kilomètre. La couverture de telles zones exige 5 de tirer des câbles multibrins, lourds et longs, dont la mise en place peut être complexifiée par la présence d'une forte topographie ou de végétation. Les temps d'installation sont donc très longs et l'opération est couteuse, à la fois en temps et en main d'ceuvre. Pour s'affranchir de ces contraintes liées au système de mesures du type « centralisé », on connaît des procédés utilisant des enregistreurs de résistivité électrique 10 et de chargeabilité du type « distribué » comportant : des récepteurs communiquant avec la centrale d'acquisition pour lui envoyer les données en temps réel ; des centrales nécessitant que toutes les mesures de différence de potentiel soient faites par rapport à une électrode de référence ; ceci impose un cordon entre chaque 15 récepteur monovoie, et les récepteurs sont communicants et transmettent les informations à la centrale d'acquisition. Ainsi les dispositifs et méthodes d'analyse associées comportent au moins l'un des inconvénients suivants : - une potentielle difficulté de communication en cas de forte topographie ou de 20 zone investiguée trop grande ; - une profondeur d'investigation limitée ; une perte de temps lors de l'installation des câbles nécessaires pour relier les différents éléments du système de mesure ; la restitution d'un signal traité sans possibilité de pouvoir revenir à la donnée 25 brute. L'invention a pour but de remédier à ces inconvénients en fournissant un procédé pour construire une image d'une zone du sous-sol, représentative de la résistivité et/ou de la chargeabilité de ladite zone, à partir de systèmes d'acquisition, dans lequel les systèmes d'acquisition sont indépendants les uns des autres et positionnés de façon distribuée sur 30 la zone. Le procédé comporte les étapes suivantes : on enregistre simultanément des mesures d'intensité du courant et des mesures de différences de potentiel entre deux dipôles orthogonaux ; on stocke lesdites mesures dans lesdits systèmes d'acquisition en y associant un temps d'enregistrement et une position géographique ; 35 on organise lesdites mesures selon un temps d'enregistrement, puis on filtre et on accumule lesdites mesures ; et 3064370 3 - on construit l'image de la zone du sous-sol en déterminant des données apparentes de résistivité et de chargeabilité, à partir desdites mesures accumulées et filtrées. Le procédé comporte ainsi les avantages suivants : 5 - Grande profondeur d'investigation Les systèmes d'acquisition indépendants n'ayant pas besoin de communiquer (de façon filaire ou sans fil) pendant l'acquisition, la distance entre eux n'est pas limitée. Le déploiement du dispositif sur une surface de grande dimension permet d'accéder à une profondeur d'investigation élevée. L'état de l'art est limité quant à lui par la dimension des câbles. Cet atout est particulièrement prisé dans le secteur de l'exploration minière dont les cibles sont de plus en plus profondes. - Gain de temps lors du levé Un autre avantage de l'invention réside dans le temps gagné dans les opérations de levés que ce soit en 2D ou en 3D. Les câbles multibrins utilisés dans l'état de l'art pour les acquisitions 3D à grande profondeur sont nombreux, très lourds et souvent difficiles à mettre en place du fait de la topographie et de la végétation. L'absence de câble permet d'installer très rapidement les systèmes d'acquisition en se rendant d'un point d'installation à l'autre en véhicule motorisé quand cela est possible. - Une meilleure qualité de mesure Les systèmes n'étant pas connectés les uns aux autres, les longueurs de câbles utilisées sont petites par rapport à celles des appareils de l'état de l'art. Cela contribue à réduire fortement les effets d'antenne et les potentiels couplages électromagnétiques entre câbles et avec le sol, améliorant significativement la qualité des formes d'ondes mesurées. Alors que les appareils de l'état de l'art interprètent directement les signaux mesurés en calculant et enregistrant la résistivité et la chargeabilité, les systèmes d'acquisition de courant et de potentiel utilisés par l'invention enregistrent toutes les données en fonction du temps (séries temporelles). Cela permet à l'utilisateur d'améliorer le traitement a posteriori en focalisant le traitement sur une période moins bruitée, de sélectionner une méthode de traitement plus adaptée que celle qui aurait été utilisée avec l'appareil de l'état de l'art ou bien de rejeter la mesure du fait d'un mauvais rapport signal/bruit dans le cas d'un signal indétectable à l'oeil nu par exemple. D'autre part, compte tenu de la diversité des directions possibles du courant de retour du dispositif et donc du champ électrique qu'il provoque, le fait de mesurer en tout point les signaux dans deux directions plus ou moins orthogonales permet de disposer d'un signal mesurable sur au moins une composante, également de caractériser le champ électrique par l'amplitude de la résultante de ses deux composantes. 3064370 4 Le procédé peut en outre comporter l'une ou plusieurs des caractéristiques suivantes, prises seules ou en combinaison : on utilise des dispositifs de mesure de données de différences de potentiel et de 5 données de courant, dans lequel on utilise au moins un émetteur disposé en série avec au moins deux électrodes pour injecter un courant dans ladite zone ; on enregistre simultanément les mesures d'intensité du courant induits par ledit émetteur à partir du premier système d'acquisition, et l'on enregistre les mesures de différences de potentiel entre deux dipôles orthogonaux induits par ledit émetteur à 10 partir du second système d'acquisition ; on associe le temps d'enregistrement et la position géographique au moyen d'un système GPS et/ou d'une horloge interne du système d'acquisition ; on réalise également un contrôle qualité au moyen d'un logiciel ; les premier et second systèmes d'acquisition sont autonomes en termes 15 d'alimentation électrique et en capacité de stockage des mesures ; les premier et second systèmes d'acquisition enregistrent des séries temporelles de données ; on réalise, préalablement à l'étape de positionnement distribués des premier second systèmes d'acquisition, une étape de détermination de la position des électrodes 20 d'injection et des dipôles orthogonaux de réception, permettant à l'opérateur de disposer les systèmes d'acquisition en tenant compte des contraintes de surface de ladite zone et d'une zone illuminée en profondeur ; on organise les mesures d'intensité et les mesures de différence de potentiel, au moyen d'un logiciel ; 25 on construit l'image de la zone du sous-sol en calculant une résultante de couples de résistivité apparentes et une résultante de couples de chargeabilité apparentes. on construit l'image de la zone du sous-sol en utilisant toutes les valeurs de résistivité apparente et de chargeabilité apparente. on construit l'image de la zone du sous-sol en réalisant une inversion des données 30 apparentes de résistivité et de chargeabilité, permettant de construire des images de résistivité vraies et de chargeabilité vraies. L'invention sera mieux comprise à la lecture des figures annexées, qui sont fournies à titre d'exemples et ne présentent aucun caractère limitatif, dans lesquelles : 35 la figure 1 illustre les principales étapes du procédé selon l'invention ; la figure 2 illustre le principe d'injection et la diffusion de courant dans la zone d'étude ; 3064370 5 la figure 3A illustre le positionnement de façon distribuée des systèmes d'acquisitions indépendants utilisés par le procédé ; la figure 3B illustre, en comparaison, un système de l'art antérieur avec un positionnement centralisé ; 5 la figure 4 illustre l'implantation des électrodes de courant et de potentiel raccordées aux enregistreurs de courant et de potentiels ; la figure 5 illustre (en haut) un enregistrement de l'intensité de courant par le premier système d'acquisition lors de l'injection d'un signal carré ; la figure 4, illustre (en bas) un enregistrement du potentiel sur la voie 1 d'un second système d'acquisition durant 10 la même période de temps ; la figure 6 représente les données issues d'une campagne de mesure avec le premier système d'acquisition (en haut) et le deuxième système d'acquisition (ici trois, trois images en bas), pendant une période de temps (t) ; la figure 7 illustre le résultat de l'étape de sommation sur une période, des différences 15 de potentiel mesurées pendant toute une injection; la figure 8 illustre le résultat de l'étape de sommation de la courbe de décharge (mV/V), pendant toute une injection; Comme illustré sur la figure 1, le procédé selon l'invention permet de construire une 20 image d'une zone du sous-sol, représentative de la résistivité et/ou de la polarisation provoquée de ladite zone, à partir de dispositifs de mesure de données temporelles de potentiel et de courant. Le procédé utilise au moins un émetteur placé en série d'au moins deux électrodes pour injecter un courant dans ladite zone. La figure 2, illustre le principe d'injection et de répartition du courant dans la zone d'étude : un appareil (I) introduit un 25 courant continu ou alternatif à très basse fréquence dans le sol (S) par l'intermédiaire de deux électrodes « d'injection » (A, B) et mesure la différence de potentiel (V) entre deux autres électrodes « de potentiel » (M, N). Les lignes courbes en pointillé représentent les lignes de courant, et les lignes courbes continues représentent les surfaces équipotentielles. 30 Le procédé comporte les étapes suivantes, décrites plus en détail ultérieurement : i. on réalise une étape d'acquisition (ACQ) de mesures (levé) dans laquelle : on positionne (POS) sur la zone de travail, de façon distribuéee, un premier et un second système d'acquisition, ces systèmes d'acquisition étant indépendants les uns des autres sur la zone ; 35 on enregistre (ENR) des mesures d'intensité du courant à partir du premier système d'acquisition, et l'on enregistre des mesures de potentiels à partir du second système d'acquisition ; 3064370 6 on stocke (STO) chaque enregistrement dans les systèmes d'acquisition en y associant un temps d'enregistrement et une position géographique, au moyen d'un système GPS et/ou d'une horloge interne du système d'acquisition ; puis ii. on réalise une étape de traitement (TRA) des mesures enregistrées dans laquelle : 5 on récupère (REC) les mesures stockées dans chaque système d'acquisition ; on organise (ORG) les mesures d'intensité et les mesures de potentiel, selon le temps d'enregistrement ; et on réalise, au moyen d'un logiciel, un contrôle qualité (QC), un filtrage et une accumulation de la mesure ; 10 on détermine (DET) les paramètres de résistivité et de chargeabilité apparentes ; enfin, on restitue (IMA) une image d'une zone du sous-sol après inversion des données apparentes. Etape de préparation du levé de terrain (PREPA) 15 Selon un mode de réalisation, on réalise, préalablement à l'étape de positionnement distribués des premier et second systèmes d'acquisition, une étape de détermination de la position de l'émetteur et des systèmes d'acquisition, permettant à l'opérateur de disposer les systèmes d'acquisition en tenant compte des contraintes de surface de ladite zone et de simuler la zone illuminée en profondeur (points de report). 20 L'utilisateur peut préparer le levé de terrain à l'aide d'un algorithme qui fait partie intégrante de l'invention. Dans un premier temps, on définit des positions (DEFP) sur la zone à investiguer pour les systèmes d'acquisition et leurs électrodes. Les coordonnées correspondantes à ces positions sont enregistrées et leurs positions 25 automatiquement reportées (AFF) dans un environnement numérique de cartographie, de préférence satellitaire, de type Google Earth. On vérifie alors l'adéquation de ces positions avec les conditions de surface (habitat, infrastructures, obstacles ...). Les positions sont modifiées (MOD) jusqu'à trouver des positions compatibles avec les conditions de surface. 30 Puis, l'opérateur simule la zone illuminée (atteint par le courant injecté) en profondeur selon la disposition des systèmes d'acquisition en surface. Les positions sont à nouveau modifiées jusqu'à trouver des positions compatibles avec l'illumination requise. 35 Etape d'acquisition des mesures L'opérateur injecte, au moyen de l'émetteur, un courant dans le sol à plusieurs endroits à des périodes de temps différentes. Ce courant est injecté sous la forme de créneaux. 3064370 7 L'opérateur enregistre simultanément l'intensité de ces courants grâce au premier système d'acquisition disposé en série et les potentiels aux électrodes grâce aux seconds systèmes d'acquisition. De préférence, l'enregistrement dure plusieurs périodes afin de procéder lors du 5 traitement à l'amélioration du rapport signal sur bruit grâce à la sommation des signaux enregistrés pendant une période (accumulation ou « stacking » en anglais). Comme l'illustre la figure 2, les électrodes de l'émetteur dédiées à l'injection de courant (A et B) sont disposées selon la configuration souhaitée (pôle ou dipôle) par rapport aux électrodes de potentiels (M, N). Ces électrodes de courant (A, B) sont 10 déplacées pendant le levé dans la zone d'étude. Comme l'illustre la figure 3A, au cours de cette étape, on positionne de façon distribuée le premier système d'acquisition (IFW) et le second système d'acquisition (VFW). Une fois le premier système distribué implanté, on procède à l'injection de courant par l'intermédiaire d'une de ses électrodes E-IFW qui sera déplacée au cours du levé. 15 Ces systèmes d'acquisition VFW et IFW sont indépendants les uns des autres sur la zone, c'est-à-dire qu'il n'existe aucune connexion ni communication entre eux. La figure 3B, illustre, en comparaison, un système de l'art antérieur avec un positionnement centralisé du système d'acquisition (SA). On enregistre des mesures d'intensité du courant créé par l'émetteur (Tx) à partir du 20 premier système d'acquisition (IFW) disposé en série (Figure 4). On peut par exemple utiliser un système d'acquisition du type 1-Fullwavers® (Iris Instruments, France). On enregistre des mesures de potentiels aux électrodes à partir du second système d'acquisition (VFW). On peut par exemple utiliser un système d'acquisition du type VFullwavers® (Iris Instruments, France). Les systèmes V-Fullwaver ont été spécifiquement 25 développés pour des mesures de chargeabilité et de résistivité. Chaque système est entièrement indépendant, intégrant sa propre source d'alimentation, un module GPS et une mémoire numérique pour un enregistrement en continu jusqu'à 3 mois. Le V-Fullwaver est un enregistreur de différences de potentiel mesurée sur ses 2 voies de mesures (voie 1 mesurée entre Pl et P2 et voie 2 entre P2 et P3). L'enregistrement 30 continu à une fréquence d'échantillonnage de 10 ms fournit un enregistrement complet des séries temporelles de la différence de potentiel entre l'électrode P2 et les électrodes Pl et P3 respectivement (Figure 4). Ces trois électrodes seront de préférence disposées orthogonalement. L'enregistrement peut être synchronisé sur le signal GPS PPS, permettant un traitement facile du signal. 35 Les systèmes d'acquisition 1-Fullwaver (IFW) sont connectés en série avec l'émetteur de courant (Tx). N'importe quel émetteur peut être utilisé avec ce système. On note qu'un système d'acquisition 1-Fullwaver (IFW) est nécessaire pour l'ensemble des systèmes 3064370 8 d'acquisition V-Fullwaver (VFW) installés sur le terrain. Le I-Fullwaver est disposé en série avec l'émetteur (Tx) entre deux électrodes A et B. Les systèmes d'acquisition de courant (I-Fullwaver) et de potentiel (V-Fullwaver) utilisés par l'invention enregistrent les séries temporelles complète (systèmes 5 d'acquisition « full waveform ») en fonction du temps. Cela permet à l'utilisateur d'améliorer le traitement a posteriori en focalisant le traitement sur une période moins bruitée, de sélectionner une méthode de traitement plus adaptée que celle qui aurait été utilisée avec l'appareil de l'état de l'art ou bien de rejeter la mesure du fait d'un mauvais rapport signal/bruit dans le cas d'un signal indétectable à l'oeil nu par exemple. 10 Les mesures acquises par les systèmes d'acquisition distribués sont stockées dans les systèmes d'acquisition eux-mêmes. Ainsi, de par sa conception distribuée, c'est-à-dire que les mesures ne sont pas concentrées en un point fixe (la centrale d'acquisition), les mesures sont enregistrées par des systèmes d'acquisition distribués sur la zone à investiguer. 15 Lors de ce stockage, on associe à chaque mesure un temps d'enregistrement et une position géographique, au moyen d'un système GPS et/ou d'une horloge interne du système d'acquisition. Les systèmes d'acquisition distribués sont autonomes en termes d'alimentation électrique. Par ailleurs, ils comportent une capacité de stockage des mesures permettant 20 un stockage de mesures sur plusieurs semaines ou plusieurs mois selon la fréquence d'échantillonnage choisie. De tels systèmes d'acquisition distribués, du type 1-Fullwavers® et V-Fullwavers® comportent un écran pour visualiser en temps réel les données en cours d'enregistrement. Ils peuvent réaliser des enregistrements jusqu'à une fréquence de 100 25 Hz ce qui permet d'assurer une bonne représentativité des évolutions du milieu en particulier de leur constante de temps. Selon le procédé, les seconds systèmes d'acquisition (V-Fullwaver) sont disposés en ligne ou en nappe. En utilisant un nombre suffisant de ces systèmes pour la zone d'étude, il n'est pas nécessaire de modifier leurs places du début à la fin du levé de terrain. 30 La figure 5, illustre (en haut) un enregistrement, en fonction du temps t, de l'intensité de courant (i) par le premier système d'acquisition (IFW) lors de l'injection d'un signal carré. La figure 5, illustre (en bas) un enregistrement, en fonction du temps t, de la différence de potentiel (V) entre les électrodes P1 et P2 (voie 1) d'un second système d'acquisition (VFW) durant la même période de temps. L'intensité du courant injecté est 35 exprimée en mA et le potentiel en mV. Sur ces figures, la flèche noire du haut indique la période, et la flèche grise indique la faculté de synchroniser manuellement les signaux de courant et les signaux de tension si nécessaire. 3064370 9 La figure 6 représente les données issues d'une campagne de mesure pendant une période de temps (t) : l'intensité i avec le premier système d'acquisition (en haut), et la différence de potentiel avec les deuxièmes systèmes d'acquisition (ici trois deuxièmes systèmes d'acquisition ont été utilisés, correspondant aux trois images du bas). Le 5 courant est injecté dans le sol entre deux électrodes pendant un certain temps et est enregistré (figure 6 en haut). L'enregistreur de courant est assimilé à un ampèremètre. Sur cette mesure, on distingue 11 injections. L'intensité du courant injecté est enregistrée par le I-Fullwaver alors que la différence de potentiel (ici la voie 1 notée V01) est enregistrée simultanément par trois V-Fullwaver situés à des positions différentes (3 10 figures les plus basses). La proximité des électrodes de potentiels par rapport aux électrodes de courant et la variation de résistivité électrique du terrain se traduit par une différence de potentiel plus ou moins élevée. Etape de traitements des mesures 15 Au cours de cette étape, on calcule des résistivités et des chargeabilités à partir des données temporelles de potentiel et de courant. Pour ce faire, on récupère les mesures stockées dans chaque système d'acquisition, on organise les mesures d'intensité et les mesures de potentiel, selon le temps d'enregistrement, et on réalise, un contrôle qualité et un filtrage des mesures, puis l'on 20 construit les images de résistivités et de polarisations provoquées à partir des mesures ainsi ordonnées. Lors de cette dernière étape, on restitue une image de résistivité apparente et une image de polarisation provoquée apparente. L'utilisation d'un logiciel d'inversion 3D disponible sur le marché permet d'obtenir des images de résistivité et de polarisation 25 provoquée « vraies » ou « interprétées ». A l'issu du levé, les données enregistrées sont extraites de chaque récepteur (V et IFullwaver). Les fichiers du I-Fullwaver comportent au moins les informations suivantes : coordonnées géographiques x, y, z ; 30 temps (date d'acquisition) : t ; et intensité du courant injecté : i. Les fichiers des V-Fullwaver comportent au moins les informations suivantes : coordonnées géographiques x, y, z ; temps d'acquisition : t ; 35 différence de potentiel de la voie 1 : V1 ; et différence de potentiel de la voie 2 : V2. relatif au courant injecté : i. Au cours d'une première étape, on ordonne tous les enregistrements d'intensité (I- 3064370 10 Fullwaver) et de potentiels (V-Fullwaver) en fonction de leur position géographique (x, y, z) et du temps de leur acquisition (t). Le fichier résultant de cette fusion comportent donc au moins les informations suivantes, pour chaque position géographique (x, y, z) : coordonnées géographiques x, y, z ; 5 temps d'acquisition : t ; intensité du courant injecté : i ; différence de potentiel de la voie 1 : V1 ; et différence de potentiel de la voie 2 : V2. Cette étape peut être réalisée au moyen d'un logiciel dédié, tel que FullwaveViewer® 10 (Iris Instruments, France). On réalise ensuite une étape facultative de contrôle qualité des mesures ainsi enregistrées et ordonnées. Cette étape peut être réalisée au moyen d'un logiciel dédié, de préférence à partir du même logiciel précité. Lors de cette étape, l'opérateur peut soit procéder à des corrections automatiques (synchronisation sur le temps GPS ou sur 15 l'horloge interne de l'instrument, rejet de pics de forte amplitude ...), soit sélectionner lui- même les fenêtres temporelles sur lesquelles sont fait les traitements. Le traitement consiste en la réalisation d'une moyenne (AV) des différences de potentiel et de courant sur une période. Les résultats de l'étape de sommation de plusieurs périodes successives sont illustrés sur la figure 7, qui montre la moyenne AV 20 en fonction du temps t. Ceci permet de calculer la résistivité apparente p pour chaque rapport potentiel/intensité. La courbe de décharge (D en mV) illustrée sur la figure 8, correspond à la décroissance de potentiel (V1) enregistré sur la voie 1 (en mV) après la fin de l'injection. La chargeabilité m correspond à l'intégrale de cette courbe de décharge (mV.ms) divisée 25 par l'amplitude du voltage mesuré pendant l'injection (en V) et par la période d'intégration (ms). On réalise ensuite une étape de traitement au sens strict, au cours de laquelle on somme la différence de potentiel sur le nombre de périodes sélectionnées pour chacune 30 des 2 voies. Deux solutions s'offrent à l'utilisateur : - Soit il calcule alors pour chaque V-Fullwaver et pour chaque injection deux résistivités apparentes pl et p2 et deux chargeabilités apparentes mi et m2 : pl= k Vi/i et p2 = k V2/i où k est le facteur géométrique (m). 35 mi et m2 correspond à l'intégrale de chacune des courbes de décharge (enregistrement de la voie 1 et enregistrement de la voie 2) divisée 3064370 11 par l'amplitude du voltage mesuré pendant l'injection (en mV/V) et par la période d'intégration. - Soit il calcule une résistivité et une chargeabilité résultante pour chaque VFullwaver et pour chaque injection, qui dans le cas où les deux voies de mesure sont orthogonales entre elles et de même longueur, peut s'exprimer de la façon suivante : N/172 + V 2 P = k 2 et m = \k/n1171)2 + (m2172)2 .\/1712 + V2 2 15 L'opérateur dispose donc d'un ou de deux sets de données de résistivité et de chargeabilité dites apparentes. Afin d'obtenir des valeurs de résistivité et de chargeabilité « interprétées », on peut compléter le procédé en utilisant un logiciel d'inversion, connu des spécialistes. Le ou les sets de données (p, et p2, rn, et m2 ou p' et m') peuvent être introduits directement dans le logiciel d'inversion. En choisissant d'introduire la résistivité 20 et la chargeabilité résultantes, l'utilisateur réduira le nombre de données par deux mais aura pour effet de régulariser les données, en supprimant les minimisant le poids de mesure faites sur de très faibles différences de potentiel. 5 10 REVENDICATIONS1. Procédé pour construire une image d'une zone du sous-sol, représentative de la résistivité et de la chargeabilité de ladite zone, à partir de systèmes d'acquisition (IFW, VFW), caractérisé en ce que lesdits systèmes d'acquisition (IFW, VFW) sont indépendants les uns des autres et positionnés de façon distribuée sur ladite zone ; et en ce que l'on réalise les étapes suivantes : on enregistre (ENR) simultanément des mesures d'intensité (i) du courant et des mesures de différences de potentiel (V) entre deux dipôles orthogonaux ; on stocke (STO) lesdites mesures (i, V) dans lesdits systèmes d'acquisition en y associant un temps d'enregistrement et une position géographique ; on organise (ORG) lesdites mesures (i, V) selon un temps d'enregistrement, puis on filtre et on accumule lesdites mesures (i, V) ; et on construit (IMA) l'image de la zone du sous-sol en déterminant (DET) des données apparentes de résistivité et de chargeabilité, à partir desdites mesures accumulées et filtrées. 2. Procédé selon la revendication 1, dans lequel on utilise des dispositifs de mesure de données de différences de potentiel et de données de courant, dans lequel on utilise au moins un émetteur disposé en série avec au moins deux électrodes pour injecter un courant dans ladite zone. 3. Procédé selon la revendication précédente, dans lequel on enregistre (ENR) simultanément les mesures d'intensité (i) du courant induits par ledit émetteur à partir du premier système d'acquisition (IFW), et l'on enregistre les mesures de différences de potentiel (V) entre deux dipôles orthogonaux induits par ledit émetteur à partir du second système d'acquisition (VFW). 4. Procédé selon l'une des revendications précédentes, dans lequel on associe le temps d'enregistrement et la position géographique au moyen d'un système GPS et/ou d'une horloge interne du système d'acquisition. 5. Procédé selon l'une des revendications précédentes, dans lequel on réalise également un contrôle qualité (QC) au moyen d'un logiciel. 6. Procédé selon l'une des revendications précédentes, dans lequel les premier (IFW) et second (VFW) systèmes d'acquisition sont autonomes en termes d'alimentation électrique et en capacité de stockage des mesures. 7. Procédé selon l'une des revendications précédentes, dans lequel les premier (IFW) et second (VFW) systèmes d'acquisition enregistrent des séries temporelles de données. 8. Procédé selon l'une des revendications précédentes, dans lequel on réalise, 3064370 13 préalablement à l'étape de positionnement (POS) distribués des premier (IFW) et second (VFW) systèmes d'acquisition, une étape de détermination de la position (DEFD) des électrodes d'injection et des dipôles orthogonaux de réception, permettant à l'opérateur de disposer les systèmes d'acquisition (IFW, VFW) en 5 tenant compte des contraintes de surface de ladite zone et d'une zone illuminée en profondeur. 9 Procédé selon l'une des revendications précédentes, dans lequel on organise (ORG) les mesures d'intensité (i) et les mesures de différence de potentiel (V) , au moyen d'un logiciel. 10 10 Procédé selon l'une des revendications précédentes, dans lequel on construit l'image de la zone du sous-sol (IMA) en calculant une résultante de couples de résistivité apparentes et une résultante de couples de chargeabilité apparentes. 11 Procédé selon l'une des revendications 1 à 9, dans lequel on construit l'image de la zone du sous-sol (IMA) en utilisant toutes les valeurs de résistivité apparente et de 15 chargeabilité apparente. 12. Procédé selon l'une des revendications précédentes, dans lequel on construit l'image de la zone du sous-sol (IMA) en réalisant une inversion des données apparentes de résistivité et de chargeabilité, permettant de construire des images de résistivité vraies et de chargeabilité vraies. 20