FR 2614997 A1 19881110 FR 8706544 A 19870507 - 1 - L'invention a pour objet une méthode de prospection sismique permettant une connaissance améliorée de l'emplacement et de la nature des discontinuités géologiques du sous-sol. Les méthodes de prospection sismique comportent généralement une étape d'émission d'ondes sismiques, une étape de réception en un ou plusieurs emplacements des ondes qui se sont propagées dans le terrain à étudier, une étape d'enregistrement des signaux reçus et une étape d'application aux signaux enregistrés, de traitements de natures très diverses avant de dresser des coupes sismiques o les discontinuités du sous-sol sont restituées à leur vraie place. Une méthode connue de prospection sismique comporte l'émission d'ondes sismiques par une source disposée soit à la surface du sol soit dans un autre puits et la réception des ondes après propagation en un ou plusieurs emplacements d'un puits foré dans le sol. La réception est effectuée au moyen d'un dispositif de réception comportant au moins une sonde de faible section, suspendue par un câble électro-porteur à un ensemble de support disposé en surface, et munie de bras d'ancrage escamotables pour plaquer le corps de la sonde contre la paroi du puits. Dans un compartiment du corps de sonde sont disposés un ou plusieurs capteurs sismiques directionnels. Ces capteurs sont généralement associés par trois et leurs axes disposés suivant trois axes orthogonaux. -2- Le dispositif de réception peut comporter aussi plusieurs sondes de puits suspendues les unes au-dessous des autres et contenant chacune au moins un capteur triaxial, de manière-à multiplier les emplacements de réception des ondes sismiques. La sonde ou l'ensemble de sondes est descendu(e) dans le puits et remonté(e) par paliers successifs. A chaque pause, on effectue un cycle d'émission d'ondes sismiques. Les ondes qui se sont propagées jusqu'à la sonde ou l'ensemble des sondes, sont reçues par les capteurs triaxiaux et les signaux sismiques sont collectés par un ensemble d'acquisition qui les transmet à un laboratoire de surface. De telles sondes de puits sont décrites par exemple dans les brevets français 2.501.380 et 2.548.727 et les demandes de brevets français EN. 87/04.365 et 87/04.677. Les signaux, reçus sont en général complexes. Les géophysiciens ont cherché à déterminer la courbe décrite par le centre de gravité de chaque récepteur tri-axial et désignée ci-après par trajectoire, à partir des signaux reçus par les trois capteurs élémentaires d'un même récepteur et ont trouvé que la surface qui pouvait le mieux la représenter, au moins pendant un certain temps à partir du début de la réception, était un ellipsolde généralement très aplati. On a jusqu'ici considéré que cet ellipsoide aplati modélisant la trajectoire réelle, était coritenu dans un plan sensiblement vertical, les écarts pouvant être négligés, et que ce plan vertical était celui qui passait par la source sismique et l'emplacement de chaque récepteur. La sonde ou chacune d'elles étant susceptible de tourner plus ou moins sur son axe pendant les manoeuvres de descente ou de remontée, le - 3 trièdre constitué des axes respectifs des trois capteurs élémentaires de chaque récepteur, peut avoir une orientation quelconque par rapport au plan vertical passant par la source et l'emplacement de réception. Dans l'hypothèse admise, on avait à comparer et à combiner les signaux délivrés par les capteurs d'un même récepteur, de manière à déterminer l'angle dièdre du plan vertical de la trajectoire avec un plan vertical de référence formé par deux axes du trièdre et l'angle entre les axes de symétrie de la trajectoire et le plan horizontal. Ces deux angles que, par analogie on peut désigner par longitude et latitude, suffisaient pour caractériser la position de la trajectoire. La connaissance de ces deux angles et d'autres paramètres tels que l'ellipticité de la trajectoire permettait dans l'hypothèse retenue, de déterminer certaines caractéristiques géophysiques et géologiques des terrains traversés par les ondes. La méthode de prospection sismique selon l'invention permet d'améliorer de façon très sensible la connaissance des caractéristiques géologiques des terrains traversés. Elle comporte l'émission d'ondes acoustiques dans le sol, la réception des ondes acoustiques qui se sont propagées dans le sol par au moins un récepteur sismique triaxial et l'enregistrement des ondes reçues. Elle est caractérisée en ce qu'elle comporte la détermination à partir des différents signaux délivrés par chaque récepteur, de paramètres définissant la position réelle dans l'espace et la forme de la trajectoire suivie par le centre de gravité de chaque récepteur, en réponse aux ondes reçues, la détermination en fonction de la distance entre chaque récepteur et le lieu d'émission, des variations desdits paramètres et la localisation des discontinuités géologiques aux profondeurs o l'on a mesuré des variations significatives desdits paramètres. La trajectoire étant modélisée par son "ellipsolde équivalent" on combine par exemple les ondes reçues par les trois capteurs de chaque - 4 récepteur, pour déterminer sa position dans l'espace et des paramètres indicatifs de sa forme, et on établit les courbes représentatives des variations des différents paramètres en fonction de la distance entre chaque récepteur et le lieu d'émission des ondes sismiques. Chaque récepteur est par exemple appliqué contre la paroi d'un forage et le lieu d'émission est situé à la surface du sol. Dans ce cas, sa distance au lieu d'émission est sensiblement sa profondeur, et c'est en fonction de celle-ci que lesdits paramètres sont représentés. D'autres caractéristiques et avantages de la méthode apparaîtront à la lecture de la description d'un mode de réalisation donné à titre d'exemple non limitatif, en se référant aux dessins annexés o: - la figure 1 montre un ensemble de réception comportant plusieurs sondes, descendu dans un forage; - la figure 2 montre les projections sur trois plans orthogonaux, d'une trajectoire ayant une polarisation privilégiée dans un plan vertical, pendant un intervalle de temps limité. - la figure 3 montre des projections sur trois plans orthogonaux d'une trajectoire relevée pendant un intervalle de temps limité, - la figure 4 montre une représentation de la trajectoire du centre de gravité des trois capteurs d'un même récepteur triaxial sous l'action des ondes sismiques captées, polarisée dans un plan vertical, o représente la latitude et, la longitude. - la figure 5 représente également des projections sur trois plans orthogonaux d'une trajectoire relevée pendant un intervalle de temps limité, et illustre le cas général o le plan de polarisation de la trajectoire est orienté dans l'espace de manière quelconque; et -5 - - la figure 6 montre une trajectoire quelconque dans l'espace et ses paramètres représentatifs. Des opérations de prospection sismique du sous-sol peuvent être réalisées en descendant dans un forage 1 (fig. 1) à l'extrémité d'un câble électroporteur multi-fonctions 2, un dispositif de réception comportant une sonde unique 3 sous laquelle éventuellement sont suspendues des sondes satellites 4. Dans chacune des sondes, est disposé au moins un récepteur sismique 5 comportant au moins un capteur directif. On utilise généralement des récepteurs triaxiaux comportant trois capteurs directifs dont les axes sont dirigés suivant trois directions orthogonales les unes par rapport aux autres et forment un trièdre trirectangle. Les sondes sont munies de bras d'ancrage et peuvent être immobilisées dans le forage à toute profondeur o des relevés sismiques doivent être effectués. Le câble éLectro-porteur 2 relie la sonde 3 ou l'ensemble de sondes (3, 4) à un ensemble de surface comportant une structure de support 6 et un treuil 7 disposé sur un véhicule 8. Les signaux captés par chaque récepteur sont transmis par des conducteurs électriques du câble, à un système central de commande et d'enregistrement 9 disposé sur le véhicule. La connexion est directe si le nombre de récepteurs est faible. Elle est assurée dans un autre cas par l'intermédiaire d'un appareil d'acquisition de données disposé dans une des sondes et adapté à collecter, numériser et coder les signaux sismiques à transmettre. De tels ensembles de réception pour forage sont décrits dans les demandes de brevets ou brevet français précités. La sonde ou l'ensemble de sondes étant plaqué contre la paroi du forage à une profondeur choisie, on actionne une source sismique 10 disposée au contact du sol et on procède à l'enregistrement des ondes qui se sont propagées jusqu'aux différents récepteurs. Le déplacement des particules de terrain dans la paroi du forage - 6 contre laquelle sont appliqués les récepteurs sismiques, est en général fort complexe. Le mouvement du centre de gravité de chaque récepteur sismique triaxial peut être représenté comme le montre la figure 2, par sa projection sur trois plans orthogonaux, le plan vertical X 02 contenant l'emplacement d'un récepteur sismique et celui de la source sismique 10, le plan vertical y'O? perpendiculaire au précédent et le plan horizontal Xoy. On désignera ci-après par trajectoire ce mouvement du centre de gravité. Si on limite la durée d'enregistrement des signaux, ta trajectoire décrite peut être représentée de manière plus simple par une courbe dans l'espace qui peut, avec une bonne approximation, être considérée comme un ellipsoide très aplati. On voit sur la figure 3 que la projection de la trajectoire sur le plan X OZ est quasiment une ellipse alors que les projections sur les deux plans orthogonaux au premier (yOZtet XOy) sont presque réduites à des segments de droite. L'exploitation des enregistrements obtenus a été faite jusqu'ici en partant de l'hypothèse que l'ellipsolde quasi-plan parcouru par la trajectoire, était contenu dans un plan vertical et que ce plan passait par l'emplacement de réception et celui de la, source sismique S (10). Du fait de la torsion du câble électro-porteur durant les manoeuvres de montée, la position du trièdre (H1, H2, H3) formé par les axes des capteurs d'un récepteur triaxial (fig. 4) peut être quelconque. L'axe H3 pouvant être considéré comme toujours vertical, l'indétermination porte en fait sur la position du repère (H1, H2). Dans l'hypothèse en vigueur o la trajectoire était considérée comme quasiment contenue dans le plan vertical XOZ passant par le point d'émission S (fig. 4), un seul angle A qui correspond par analogie à une longitude, suffisait à caractériser le plan de polarisation de la trajectoire. Le traitement des données sismiques lues sur les enregistrements, - 7 - consistait dans ce cas à combiner les trois signaux de chaque récepteur de manière à calculer L'angle k et aussi l'angle Y, analogue à une latitude, représentant l'inclinaison du grand axe de l'ellipse par rapport à l'horizontale. La méthode selon l'invention permet, par une prise en compte plus réaliste de la complexité des trajectoires détectées, de mesurer tout au long du forage un nombre plus important de paramètres traduisant généralement des variations lithologiques et géologiques des différentes couches traversées et donc d'obtenir une représentation plus détaillée du sous-sol. La figure 5 montre en effet que la trajectoire ne se projette pas de façon privilégiée dans le plan XOE comme dans le cas représenté à la figure 3 mais dans les deux plans complémentaires yC>i et XOY et que c'est la projection dans le plan XOy qui est ici réduite presque à un segment. On le vérifie aisément par la constitution et l'examen d'un couple stéréoscopique de la trajectoire qui montre que le plan moyen n'est plus vertical. Ceci s'explique par le fait que dans la plupart des cas, le plan de polarisation privilégié de la trajectoire n'est pas situé dans le plan vertical passant par chaque récepteur et l'emplacement d'émission 10 mais qu'il présente une certaine obliquité par rapport à la verticale. L'ellipsoîde T représentant le mieux la trajectoire observée doit donc être défini par un plus grand nombre de paramètres, comme le montre la figure 6. 1 désigne l'axe principal de polarisation de l'ellipsoide; V 2, son second axe orthogonal au premier et t3 ' son troisième axe perpendiculaire au plan formé par les deux précédents. Pour définir la position de ces trois axes, on peut utiliser par exemple les angles d'Euler. L'angle O dit de nutation est celui que fait l'axe principal V1 avec la vertical OQ. Le plan normal à l'axe principal 1 coupe le plan horizontal suivant une ligne hh' qui fait un angle YJ dit angle de précession, avec 2 6 1 4 9 9 7 - 8 - l'axe H1 d'un des deux capteurs horizontaux de chaque récepteur sismique triaxial. Cette même ligne hh' 'fait avec le troisième axe 3 un angle C dit de rotation propre. On définit aussi un quatrième angle 3 ditd'inclinaison latérale. C'est l'angle qui fait le troisième vecteur V3 avec le plan horizontal. On définit une portion de trajectoire par un certain nombre n d'échantillons de signaux à trois composantes. A partir des valeurs des signaux enregistrés par les trois capteurs d'un même récepteur triaxial et par une méthode de transformation connue en soi, on calcule la matrice de covariance des échantillons sur un intervalle de temps déterminé et l'on obtient trois valeurs propres 1, 2, l 3 et les trois vecteurs propres V1l V2' I /3 définis chacun par trois cosinus-directeurs, par rapport aux trois axes H1, H2, i du trièdre lié à chaque récepteur sismique. Au vecteur propre V. (i = 1,2 ou 3) correspondent les cosinus- directeurs Vi i2/et V. On peut alors calculer les racines carrées des trois valeurs propres précédemment déterminées, racines qui représentent les longueurs a, b et c des trois demi-axes de l'ellipsoîde équivalent, à un coefficient constant près. Les valeurs des quatres angles précédemment définis, se calculent alors par les formules suivantes: = cos 1( V1) Y = tan-1( V11/vH2) (D = tan-1( V3 / V2) 3 = sin1 ( V 3 On détermine aussi d'autres paramètres de polarisation tels que notamment les coefficients d'ellipticité e21, e31, e32 en faisant les - 9 - les rapports entre les valeurs des trois demi-axes a, b, c précédemment calculés, ou bien encore un coefficient de planéité p en combinant les valeurs des trois demi-axes a, b, c. Les valeurs des différents paramètres ayant été calculées pour chaque récepteur sismique triaxial, on peut donc déterminer leurs courbes représentatives en fonction de la profondeur de réception et en faire des enregistrements. La lecture de ces enregistrements de paramètres facilite beaucoup l'interprétation des enregistrements sismiques dans la mesure o il existe une corrélation manifeste entre les profondeurs o se situent des zones de discontinuités géologiques et celles o l'on observe des variations importantes dans les courbes représentatives des paramètres. On peut observer en particulier: - qu'un accroissement de l'ellipticité principale e21 de l'ellipsoîde va de pair avec une atténuation accrue; et - qu'une décroissance de l'angle de nutation O est associée à une diminution locale de la vitesse de propagation des ondes sismiques. En tenant compte des courbes de variation des paramètres, on peut donc préciser la position des horizons sur les enregistrements sismiques. Dans l'exemple considéré, les récepteurs triaxiaux sont disposés dans un forage. Il va de soi que la méthode selon l'invention pourrait être appliquée en disposant les récepteurs en surface. - 10 - R E V E N D I C A T I 0 N S 1. - Méthode de prospection sismique permettant de connaître avec une meilleure précision la nature des discontinuités géologiques du sous-sol, comportant l'émission d'ondes acoustiques dans le sol, la réception des ondes acoustiques qui se sont propagées dans le sol par au moins un récepteur sismique triaxial et l'enregistrement des ondes reçues, caractérisée en ce qu'elle comporte la détermination à partir des différents signaux délivrés par chaque récepteur, de paramètres définissant la position réelle dans l'espace et la forme de la trajectoire suivie par le centre de gravité de chaque récepteur, en réponse aux ondes reçues, la détermination en fonction de la distance entre chaque récepteur et le lieu d'émission, des variations desdits paramètres et la localisation des discontinuités géologiques aux profondeurs o l'on a mesuré des variations significatives desdits paramètres. 2. - Méthode selon la revendication 1, caractérisée en ce que ladite trajectoire étant modélisée par un ellipsoide, on combine les amplitudes des ondes reçues par les trois capteurs de chaque récepteur pour déterminer sa position dans L'espace et des paramètres indicatifs de sa forme, et on établit les courbes représentatives des variations des différents paramètres en fonction de la distance entre chaque récepteur et le lieu d'émission des ondes sismiques. 3. - Méthode selon la revendication 1 ou 2, caractérisée en ce que chaque récepteur est appliqué contre la paroi d'un forage et le lieu d'émission est situé à la surface du sol, la distance entre chaque récepteur et le lieu d'émission étant sensiblement sa profondeur. 4. - Méthode selon la revendication 1 ou 2, caractérisée en ce que les récepteurs sont disposés en surface, ladite distance étant leur déport longitudinal par rapport au lieu d'émission. - 11 - 5. - Méthode selon la revendication 1, caractérisée en ce que l'on définit La trajectoire par n points en prélevant n échantillons sur les signaux délivrés respectivement par les trois capteurs de chaque récepteur triaxial, on calcule la matrice de covariance de ces n échantillons à trois composantes, de manière à déterminer les valeurs propres et vecteurs propres définis par leurs cosinus-directeurs respectifs et on combine les valeurs propres et les cosinus directeurs obtenus de manière à déterminer les angles repérant la position réelle dans l'espace de l'ellipsoîde modélisant la trajectoire et des paramètres indicatifs de sa forme. 6. - Méthode selon la revendication 5, caractérisé en ce que l'on détermine notamment des coefficients définissant l'ellipticité et la planéité de la trajectoire, on détermine leurs variations en fonction de la distance entre le lieu d'émission et l'emplacement de chaque récepteur et on localise des discontinuités géologiques pour des valeurs de la distance o l'on observe des variations importantes des valeurs desdits coefficients.