Brevet FR2784756A1 — publié le 20000421

Méthode de détection de ruptures dans des signaux diagraphiques concernant une zone d'un milieu La présente invention concerne une méthode de détection de ruptures dans des signaux diagraphiques concernant une zone d'un milieu, les signaux diagraphiques étant constitués de logs de natures différentes enregistrés pour ladite zone en fonction de la profondeur, et l'application de cette méthode à un recalage en profondeur, desdits logs. Dans de nombreux domaines, il est nécessaire de coITéler rapi- dement deux ou plusieurs courbes représentant les variations d'une première grandeur en fonction d'une seconde grandeur, pour des besoins de comparaison, d'ajustement, etc. Les courbes à comparer peuvent être de même nature, c'est-à-dire représenter les variations d'une même première grandeur en fonction d'une même deuxième grandeur, ou de nature différente. Il peut par exemple s'agir des enregistrements d'un même phénomène physique mais décalés dans le temps ou l'espace, ou des enregistrements relatifs à des phénomènes physiques différents ou encore des enregistrements relatifs à un même phénomène physique enregistré par exemple par des méthodes différentes de sorte que leur contenu firéquentiel est différent. Les conrrélations peuvent s'effectuer numériquement. Le résultat obtenu est en généraI global et peu fiable si on ne fait pas d'hypothèses contraignantes sur les signaux, la méthode consistant alors à choisir entre plusieurs pics d'auto-corrélation. La corrélation peut être effectuée visuel- lement, par décalage manuel d'une des courbes par rapport à l'autre le long de l'axe de la deuxième grandeur. On recherche ainsi par décalages successifs la similitude optimale sur une ou plusieurs portions de la courbe. Cette méthode permet de tenir compte des connaissances a priori. C'est celle qui est couramment employée en géophysique pour le calage en profondeur ou en temps des horizons sismiques ou pour la corrélation d'enregistrements effectués dans un puits et d'enregistrements sismiques. 2 2784756 Le principal inconvénient d'une telle méthode réside dans la difficulté qu'il y a à comparer des signaux de formes qui peuvent être ntrès différentes, par exemple si leur contenu firéquentiel est différent. On connaît un procédé d'analyse d'un signal dit procédé d'analyse par ondelettes, permettant de décomposer ledit signal comme une somme de fonctions ondelettes élémentaires aab. qui vibrent chacune comme des sinusoïdes sur une plage dont la position sur un axe est liée au paramètre b et dont la largeur est liée au paramètre a (fréquence centrale), et qui s'amolrtissent très follrtement à l'extérieur de cette plage. La décomposition d'un signal à l'aide d'une famille de ces ondelettes constitue ce qu'on appelle une analyse "temps-fiéquence" car les premières et plus courantes décom- positions ont été effectuées sur des enregistrements des variations d'une première grandeur en fonction du temps (la deuxième grandeur). Dans ce cas, la dimension du paramètre b est celle d'un temps et la dimension du paramètre a est la dimension de l'inverse d'un temps, donc d'une firéquence temporelle. Pour plus d'informations sur la décomposition en ondelettes ou les analyses "temps-fiéquence", on se référera à l'article "L'analyse par ondelette" de Yves MNIEYER et al., paru dans "Pour la Science" de Septembre 1987, à l'ouvrage "Wavelets" de J.N. COMBES et al. aux éditions Springer-Verlag, ou encore à la demande de brevet internationale publiée sous le n WO 92/18941, documents qui sont intégrés à la présente demande. Plusieurs types de fonctions peuvent être utilisés, permettant de définir de nombreuses familles d'ondelettes ayant des propriétés différentes. Celles-ci peuvent être par exemple des gaussiennes, des fonctions créneaux ou triangles, des fonctions réelles ou complexes, orthogonales ou non entre elles. On se référera à l'article cité précédemment pour connaître les contraintes applicables à ces diverses fonctions et à d'autres pour générer des familles d'ondelettes. Pour une famille d'ondelettes 'Fab déterminée, on définit la "transformée en ondelettes" à deux dimensions z et x associée à un enregistrement s(z) selon l'axe z comme la suite des coefficients Cab correspondant chacun à l'intégrale du produit de l'enregistrement s(z) à analyser par l'ondelette élémentaire d'analyse 'Pa b selon les valeurs de b 3 2784756 suivant l'axe z et les valeurs de a suivant un axe x. Dans le cas o l'on a choisi des ondelettes complexes pour effectuer l'analyse tempsfréquence d'un enregistrement ou d'un signal, il devient possible de définir la partie réelle, la partie imaginaire, le module ou encore la phase de la transformée en ondelettes. Les coefficients Ca,b sont calculés par la formnule bien connue: Ca,b = _ s(:)Y+,b ()c& Des méthodes et des dispositifs pour identifier des structures géologiques en utilisant des transformées en ondelettes sont notamment décrits dans les brevets US 5 673 191, US 5 740 036 et US 5 757 309 et dans l'article intitulé "Detection of non stationarites in geological time series: wavelet transform of chlaotic and cyclic sequences", de Andreas PROKOPH et al, paru dans Computers and Geosciences, Vol. 22, No. 10, pages 1097-1108, 1996. Toutefois, ces derniers documents concernent soit des mesures magnétiques et gravitationnelles pour faire la distinction entre les structures géologiques relativement profondes et les structures peu profondes, soit des moyens pour simuler la succession des structures. La présente invention concerne une méthode de détection de ruptures dans des signaux diagraphiques, qui utilise une analyse par ondelettes desdits signaux. On sait que l'analyse des signaux diagraphiques obtenus à l'aide de dispositifs bien connus, permet de détermnniner la nature minéralogique, la texture, le type de réseau poreux et le contenu fluide des formnations traversées par des forages. L'évolution des signaux suivant la profondeur reflète l'évolution des propriétés des formations et permet de retracer leur histoire sédimentaire diagénétique et structurale. Sur les signaux diagraphiques, il est possible de distinguer des ruptures qui correspondent à des modifications importantes de la nature des formations qui s'effectuent sur un faible intervalle de profondeur. Des bancs électrofaciologiques peuvent être caractérisés à partir des ruptures repérées sur au moins un des canaux du signal diagraphique. A l'intérieur d'un banc, chaque canal du signal diagraphique montre une variation continue, à une échelle de résolution de profondeur donnée. Les ruptures remarquables sont utilisées par le géologue à des fins de corréla- tions lithostratigraphiques. Dans certains cas, des corrélations chronostrati- graphiques sont possibles en effectuant une interprétation à partir d'un modèle conceptuel de l'évolution des dépôts sédimenitaires. Les spécialistes effectuant l'analyse des signaux diagraphiques, utilisent les ruptures remarquables, dans la première étape de l'inter- prétation, pour recaler en profondeur les différents signaux enregistrés par les capteurs du dispositif diagraphique qui ne se trouvent pas tous au même moment devant une même formation. De part et d'autre des ruptures, le signal diagraphique est affecté par un effet d'épontes ("shoulder effect" en anglais) sur un intervalle qui dépend de la résolution des dispositifs diagra- phiques et du contraste des réponses diagraphiques caractéristiques des fonnrmations. Cet effet d'épontes est une source d'erreurs et d'incertitude dans les interprétations. Les processus actuels d'interprétation des signaux diagraphiques sont basés sur un traitement de chaque échantillon du signal diagraphique indépendamment des échantillons situés au-dessus et au-dessous de l'échan- tillon traité, la notion de profondeur n'intervenant pas. En conséquence, l'information portée par l'évolution du signal avec la profondeur n'est pas prise en compte. Pour que cette information soit prise en compte, il est nécessaire de définir des ruptures sur le signal diagraphique et des évolutions à l'intérieur des ruptures. La détermination des ruptures est actuellement effectuée manuel- lement et elle nécessite un opérateur expérimenté. Le résultat est à la fois subjectif et difficilement reproductible à l'identique. Or, ces ruptures qui correspondent aux limites de bancs ou de formations sont nécessaires à un recalage en profondeur. Le recalage en profondeur est une étape fondamentale de toute interprétation des logs, car il consiste à remettre à la méême profondeur des mesures effectuées par les différents capteurs des dispositifs diagraphiques, qui ne passent pas simultanément devant le même point du puits. On distingue deux types de recalages - les recalages "intra run" qui concernent des mesures enregistrées au cours d'une même remontée d'un ensemble de capteurs liés mécaniquement entre eux; - les recalages "entre runs", qui sont facilités en enregistrant toujours une diagraphie ou log commun dans les différents runs, ce log commun étant généralement le log "Gamma ray" qui selt de référence de profondeur. Un premier recalage est effectué au moment de l'acquisition et ne concerne que les mesures effectuées au cours du même enregistrement. Il n'est satisfaisant que dans les meilleurs cas et nécessite d'être toujours vérifié. Il existe des stations d'analyse de logs à l'aide desquelles il est possible de faire des recalages. Toutefois, les opérations de recalage restent manuelles ou lorsqu'elles sont automatiques, elles ne concernent que des logs de même nature, les stations d'analyse ne pouvant pas analyser automatiquement des logs de nature différente. De ce fait, seuls les recalages entre runs sont possibles. De plus, les procédés actuels, basés sur des corrélations, ne permettent pas d'identifier, de hiérarchiser et d'affecter un index de qualité aux corrélations. La présente invention a pour but de proposer une méthode qui permet de détecter automatiquement les ruptures dans les signaux diagra- phiques ou logs et qui est susceptible d'êtr-e appliquée pour un recalage en profondeur des logs enregistrés. La présente invention a pour objet une méthode de détection de ruptures dans des signaux diagraphiques concernant une zone d'un milieu et constitués de logs de natures différentes enregistrés pour ladite zone en fonction de la profondeur, du type consistant à: - sélectionner une portion dans chacun desdits logs de telle sorte que toutes les portions sélectionnées aient en commun un même intervalle profondeur, une des portions sélectionnées étant considérée comme portion de référence, déterminer une suite de fréquences spatiales d'analyse, - sélectionner une fonction ondelette mère et construire, à partir de ladite fonction mère, une famille de fonctions ondelettes d'analyse dépen- dant de la firéquence spatiale (ou nombre d'onde) et de la profondeur, 6 2784756 - calculer une transformée en ondelettes de la portion sélectionnée de chaque log et pour chaque fréquence d'analyse, - choisir une grandeur caractéristique de la transformée en ondelettes et à utiliser cette grandeur comme représentation de la transformée en ondelettes, ladite méthode étant caractérisée en ce qu'elle consiste en outre à - calculer pour chaque portion de log sélectionné et pour chaque cote profondeur, la valeur absolue du gradient moyen de la grandeur caractéris- tique de la transformée en ondelettes pour les différentes fréquences d'analyse, - sélectionner, pour chaque portion de log traitée, les pics de la valeur absolue du gradient moyen de la grandeur caractéristique, chaque pic correspondant à une rupture, - déterminer sur la portion de log de référence les ruptures correspondantes, - définir une fenêtre d'analyse centrée sur chaque rupture de la portion de log de référence, et sélectionner les ruptures des aunes portions de logs qui sont situées dans la fenêtre d'analyse. Un avantage de la présente invention réside dans le fait que toutes les courbes représentatives des différents logs enregistrés sont prises en compte et traitées rapidement (simultanément ou séquentiellement les unes après les autres). De plus, la méthode selon l'invention permet de s'affranchir des effets d'épontes liés à la résolution (efficacité) des dispositifs de diagraphie pour découper en bancs les intervalles fournissant les logs et pour pouvoir analyser l'évolution verticale de leurs différentes caractéristiques géologiques. Selon une autre caractéristique, le résultat de la transformée en ondelettes est un nombre complexe et la grandeur caractéristique de la transformée en ondelettes est la partie réelle dudit nombre complexe. La fonction mère peut être par exemple une fonction du type f(z) = (1 - z2)exp (-z2/2). Selon une autre caractéristique, la valeur absolue du gradient moyen est normalisée, les pics de la valeur absolue du gradient moyen qui sont sélectionnés sont supérieurs ou égaux à un seuil prédéterminé. En 7 2784756 particulier, la valeur absolue du gradient moyen est normalisée et le seuil de sélection des pics est égal ou supérieur à 0,2. Selon une autre caractéristique, le log fournissant la portion de référence est obtenu par gammagraphie, ledit log (gamma ray) étant une excellente référence profondeur puisqu'il peut être enregistré dans tous les types de boue de forage et même à travers un tubage. Selon une autre caractéristique, chaque portion de log traitée est incluse dans un intervalle d'étude contenant un nombre donné d'échan- tillons N. En particulier, lorsque le nombre d'échantillons à traiter dans une portion de log est soit inférieur ou bien supérieur au nombre d'échantillons N de l'intervalle d'étude, le log est soit centré dans ledit intervalle et les parties vides de ce dernier sont remplies par des échantillons ayant une valeur égale à la valeur moyenne du log, ou bien il est divisé en au moins deux parties comprenant chacune un nombre d'échantillons inférieur à N, de manière à traiter chaque partie comme indiqué ci-dessus. Selon une autre caractéristique, la succession des firéquences spatiales d'analyse utilisées pour le calcul de la transformée en ondelettes, a comme limites une fréquence correspondant à une longueur d'onde de 4 m et une fréquence correspondant à une longueur de 200 m. La succession desdites fréquences est par exemple une progression géométrique. De préférence, on sélectionne dix firéquences spatiales dont les limites correspondent à des longueurs d'onde de 10 m et 100 ni. Chaque firéquence spatiale est analysée indépendamment des autres, sans filtrages successifs. De plus, le choix des fréquences permet d'avoir un nombre de coefficients d'ondelettes qui est suffisant pour réaliser une étude de leur organisation spatiale, dans le plan profondeurfréquence. De cette manière, on obtient une information à trois dimensions reliant la profondeur, les fréquences présentes dans le signal diagraphique de départ et l'amplitude de la transformée. Ainsi, la méthode selon l'invention permet d'étudier, pour chaque cote profondeur, le signal diagraphique à différentes échelles c'est-à-dire sur des intervalles de profondeur dont la taille est différente. La méthode selon l'invention permet de hiérarchiser les ruptures suivant différents critères diagraphiques et de donner des priorités (ou des critères de qualité) dans les références qui sont utilisées notamment dans le recalage profondeur. Ainsi, il est possible de traiter automatiquement aussi bien le recalage "intra run" que le recalage "entre runs". D'autres avantages et caractéristiques ressortiront mieux à la lecture d'un mode de réalisation préféré de la méthode selon l'invention, ainsi que des dessins annexés sur lesquels: - la figure 1 est une représentation des ruptures sur un ensemble de logs non recalés, comportant un log gamma ray conmmoe référence - la figure 2 est une représentation de la patile réelle des coefficients de la transfolrmée en ondelettes du log de référence - la figure 3 est une représentation du gradient de la partie réelle des coefficients de la transformnnée en ondelettes du log de référence - la figure 4 est une représentation de la valeur absolue du gradient moyen normalisé pour le log de référence. Pour un milieu donné à explorer, divers logs sont enregistrés en fonction de la profondeur et correspondent à une zone dudit milieu. Dans l'exemple de la figure 1, six logs référencés I à 6 ont été enregistrés. Le log 1 est obtenu par gammagraphie et constitue un log gamma ray; le log 2 correspond à un enregistrement représentant l'index d'hydrogène et il est usuellement référencé par NPHI. Le log 3 correspond à la densité de la roche en place dans la zone considérée, et il est référencé par RHOB; le log 4 correspond à la lenteur Dt dans ladite zone considérée; les logs 5 et 6 correspondent à des résistivités désignées par LLS pour représenter la résistivité peu profonde et par LLD pour représenter la résistivité profonde. Dans chaque log 1 à 6, on sélectionne une portion qui concerne un même intervalle profondeur, puis on choisit parmi ces portions de log une portion qui est utilisée comme portion de référence. Dans l'exemple représenté sur la figure 1, la portion de référence est celle du gamma ray. Dans une autre étape, on sélectionne une fonction ondelette mère, par exemple la fonction ondelette de NIORLET ou mieux encore la fonction appelée fonction ondelette du "Chapeau Mexicain", du type f(z) = (l-z2) exp (-z2/2) qui est équivalente à un lissage "gaussien" dont on prend la dérivée seconde. A partir de la fonction ondelette mère, on construit, ainsi que cela est connu et rappelé ci-dessus, une famille de fonctions d'analyse qui dépend de la firéquence spatiale et de la profondeur. Dans le cas présent, on choisit dix fréquences d'analyse dont les limites sont comprises entre une fréquence F1 correspondant à une longueur d'onde X de mètres et une fréquence Flo correspondant à une longueur d'onde d'environ 100 m, la succession des fréquences d'analyse F1 à Flo décrois- sante suivant une progression géométrique, par exemple de raison 1,24 Chaque log 1 à 6 ou plutôt chaque portion de log sélectionnée est ainsi traitée à l'aide de la transformée en ondelettes et ce, pour chaque firéquence d'analyse. Le résultat de la transformée en ondelettes est un coefficient représenté par un nombre complexe dont on ne conserve, comme grandeur caractéristique, que la patile réelle. La figure 2 représente la patllie réelle des coefficients de la transformée en ondelettes pour la portion de référence du log gamma ray, les valeurs étant croissantes de la gauche vers la droite. La colonne située immédiatement après le log gamma ray 1 sur la figure 2 correspond à la partie réelle des coefficients de la transforniée en ondelettes pour la première fi-équence F1, la colonne suivante correspondant à la partie réelle des coefficients de la transformée en ondelettes pour la fréquence F2, et ainsi de suite en allant vers la droite jusqu'à la fiéquence FI0. Pour chaque cote de chaque portion de log analysée, on calcule la valeur absolue du gradient moyen de la partie réelle du nombre complexe du résultat de la transformée en ondelettes. Sur la figure 3, on a représenté le gradient de la partie réelle des coefficients de la transfonrmée en onde- lettes pour la portion de référence sélectionnée dans le gamma ray, comme à propos de la figure 2. C'est ainsi que la colonne consécutive au gamma ray situé le plus à gauche, à côté de l'échelle des profondeurs en mètres, correspond au gradient de la patile réelle représentée sur la figure 2 et correspondant à la firéquence initiale F1. A chaque firéquence F2, F3... Fo10 correspond un gradient de la patllie réelle. On procède de même pour chacun des logs 2 à 6, c'est-à-dire que chaque log est traité à l'aide de la famille de fonctions ondelettes d'analyse issue de la même fonction ondelette mère, de manière qu'à chaque log corresponde une représentation du gradient de la patlie réelle de la tr-ansformée en ondelettes pour chacune des fréquences d'analyse F1 à FiO. Dans une autre étape, on calcule pour chaque portion de log sélectionnée et pour chaque cote profondeur, la valeur absolue du gradient moyen de la partie réelle choisie comme grandeur caractéristique de la transformée en ondelettes, à partir des valeurs absolues des gradients déterminées pour l'ensemble des fréquences d'analyse F1 à F0lo. En se référant à la figure 3 et pour chaque cote profondeur, par exemple la cote profondeur 1613, on calcule le gradient moyen en prenant, par exemple, la moyenne arithmétique des valeurs du gradient pour chaque cote et pour chacune des fréquences F1 à F0lo. Puis on sélectionne, pour chaque portion de log sélectionnée, les pics de la valeur absolue du gradient moyen. Sur la figure 4 représentant le gradient moyen normalisé qui correspond au seul log gamma ray, on peut voir qu'entre 1610 et 1623 m, le gradient moyen pour chaque cote profondeur comprise entre 1610 et 1623 m est relati- vement faible, le premier pic apparaissant pour la cote 1623 mni. De proche en proche, on constate qu'on a un pic significatif pour la cote 1625 nim, d'autres pics à 1638 mni, 1650 mni, 1653 mn et ainsi de suite. A chaque valeur du gradient moyen correspond une rupture recherchée. De préférence, et pour chaque cote profondeur, seuls les pics de la valeur absolue du gradient moyen supérieurs à un seuil prédéterminé sont sélectionnés. Par exemple, on normalise la valeur absolue du gradient moyen et on ne retient que les pics supérieurs au seuil 0,2. La figure I fournit, à côté de chaque portion de log sélectionnée 1 à 6, une représentation des pics sélectionnés pour chaque cote profondeur et donc des ruptures qui ont été déterminées grâce à la présente invention. La suite des ruptures dans la colonne référencée 1' correspond à la série des pics supérieurs à 0,2, pour la portion de log gamma ray. Comme on peut le constater, pour l'intervalle de la profondeur replrésentée, une rupture a été déterminée pour la cote 1623, une autre pour la cote 1625 mn, une autre pour la cote 1626 m et ainsi de suite. La colonne 2' représente la série de ruptures déterminées pour la portion de log NPHI, la colonne 3' représentant la série de ruptures déterminées pour la portion de log RHOB et ainsi de suite. Pour des raisons de clarté et d'encombrement, les potilions de log sélectionnées sont comprises entre la cote profondeur 1610 ni et la cote profondeur 1679 mni. En réalité et de préférence, les portions de log sélec- tionnées sont incluses dans un intervalle d'étude qui comprend un nombre donné N d'échantillons diagraphiques, par exemple 4096 échantillons. Si on considère que le pas d'échantillonnage est de 15,24 cm (1/2 pied), cela l l 2784756 représente une longueur de log de 624 mètres au lieu des 69 mètres indiqués sur les figures 1 à 3. Comme indiqué précédemment, on traite N d'échantillons dans chaque portion de log sélectionnée. Lorsque la portion de log à traiter comprend moins du nombre d'échantillons N fixé, 4096 dans le cas présent, la portion de log est recentrée dans l'ensemble des 4096 échantillons. Par exemple, si la portion de log sélectionnée comprend 2 000 échantillons, on crée un log de 4096 échantillons dont les 1048 premiers et 1048 derniers sont affectés de la valeur moyenne du log, les 2000 autres échantillons correspondant au log d'origine. Lorsque la portion de log sélectionnée comprend plus d'échan- tillons que le nombre fixé N, par exemple supérieur à 4096, on coupe ladite portion en deux et on traite chaque moitié comme si elle comportait moins d'échantillons que 4096, de la manière indiquée ci-dessus. Selon l'invention, il est préférable de rechercher les ruptures pertinentes dans chaque log. A cet effet, on définit une fenêtre de pertinence qui est centrée successivement sur chaque rupture du log considéré et pour chaque position de ladite fenêtre de pertinence, on calcule un coefficient de pertinence défini par l'expression: (lOgmax 10gmiiii)ielatre / (10ginax - loigmin)portion dans laquelle: logmax et logrinjl dans la fenêtre représentent les amplitudes maximale et minimale du signal diagraphique situé dans la fenêtre, logmnax et log, 11in, dans la portion de log traitée représentent les amplitudes maximale et minimale de la portion traitée du signal diagraphique. Une telle opération contribue à l'élimination des problèmes liés au bruit sur les signaux diagraphiques enregistrés. Pour chaque log, on conserve seulement les ruptures pour lesquelles le coefficient de pertinence est supérieur à une valeur donnée. Les colonnes 1' à 6' représentent les ruptures conservées à l'aide du calcul précédent pour les logs I à 6 respectivement. 12 2784756 De manière à ne conserver ou à ne sélectionner que les ruptures intéressantes qui se correspondent d'un log à l'autre, la présente invention préconise de hiérarchiser les ruptures conservées sur chacune des portions de log sélectionnées. Pour cela, on définit une fenêtre d'analyse qui recouvre un nombre donné M d'échantillons, par exemple égal à 5. Pour chaque rupture du log de référence, on centre la fenêtre d'analyse sur cette rupture et on recherche, sur les autres logs, les ruptures qui sont situées dans ladite fenêtre. Chaque rupture déterminée sur un ou plusieurs logs, dans la fenêtre considérée, est associée à la rupture du log de référence sur laquelle la fenêtre a été centrée. Après quoi, on peut hiérarchiser les ruptures entre elles, en opérant par exemple comme suit. Dans les colonnes des ruptures, on peut voir qu'une rupture déterminée pour la cote 1623 m dans la colonne 1' n'est présente également que dans la colonne 4' et absente dans les autres colonnes 2', 3', 5' et 6'. Comme elle n'est présente que dans une colonne 4' autre que la colonne de référence 1', le coefficient un lui est affecté. La rupture déterminée à 1626 m environ dans la colonne 1', se retrouve présente dans les colonnes 2', 3', 5' et 6'; le coefficient quatre lui est affecté. La rupture déterminée à 1631 m dans la colonne 1' se retrouve également présente dans chacune des autres colonnes 2' à 6': le coefficient cinq lui est affecté on répète cette opération pour chacune des ruptures de la colonne 1'. La hiérarchisation est représentée dans la colonne 7' de la figure 1. Dans cette colonne 7', on peut constater qu'il y a quatre ruptures ayant le coefficient 5, cinq ruptures ayant le coefficient quatre, et ainsi de suite. Parce qu'elles correspondent à des événements géologiques remarquables, les ruptures sont utilisées pour un recalage en profondeur des diagraphies. Un recalage profondeur des diagraphies consiste à remettre à la même profondeur des mesures effectuées par différents capteurs qui nepassent pas simultanément devant le même point du puits. Pour effectuer un recalage en profondeur automatique, on détecte sur chaque diagraphie des ruptures, de la manière indiquée ci-dessus. De même, on recherche les ruptures existant sur les autres logs visibles dans une fenêtre de profondeur centrée sur chaque rupture repérée sur la diagraphie de référence. Lorsqu'une rupture est détectée sur un ou plusieurs logs, elle est associée à la rupture détectée sur le log de réfé-rence. Les valeurs des logs sont 13 2784756 -1 recalculées de telle sorte que les ruptures associées à chacune des ruptures du log de référence apparaissent à la même cote profondeur que la rupture considérée du log de référence. 14 1 2784756 REVENDICATIONS 1. Méthode de détection de ruptures dans des signaux diagra- phiques concernant une zone d'un milieu et constitués de logs de natures différentes enregistrés pour ladite zone en fonction de la profondeur, du type consistant à - sélectionner une portion dans chacun desdits logs de telle sorte que toutes les portions sélectionnées aient en commun un même intervalle profondeur, une des portions sélectionnées étant considérée comme portion de référence, - déterminer une suite de fréquences spatiales d'analyse, - sélectionner une fonction ondelette mère et construire, à partir de ladite fonction mère, une famille de fonctions ondelettes d'analyse dépendant de la fréquence spatiale et de la profondeur, calculer une transformée en ondelettes de la portion sélectionnée de chaque log et pour chaque fréquence spatiale d'analyse, - choisir une grandeur caractéristique de la transformée en ondelettes et à utiliser cette grandeur comme représentation de la transformée en ondelettes, ladite méthode se caractérisant en ce qu'elle consiste en outre à calculer pour chaque portion de log sélectionnée et pour chaque cote profondeur, la valeur absolue du gradient moyen de la grandeur caractéristique de la transformée en ondelettes pour les différentes fréquences spatiales d'analyse, - sélectionner, pour chaque portion de log traité, les pics de la valeur absolue du gradient moyen de la grandeur caractéristique, chaque pic correspondant à une rupture, - déterminer sur la portion de log de référence les ruptures correspondantes, - définir une fenêtre d'analyse centrée sur chaque rupture de la portion de log de référence, et - sélectionner les ruptures des autres portions de logs qui sont situées dans la fenêtre d'analyse. j5 2784756 2. Méthode selon la revendication 1, caractérisée en ce que le résultat de la transformée en ondelettes est un nombre complexe et en ce que la grandeur caractéristique de la transfonnée en ondelettes est la partie réelle dudit nombre complexe. 3. Méthode selon la revendication 1, caractérisée en ce que les pics de la valeur absolue du gradient moyen sélectionnés sont supérieurs ou égaux à un seuil prédéterminé. 4. Méthode selon la revendication 3, caractérisée en ce que la valeur absolue du gradient moyen est normalisée et en ce que ledit seuil est supérieur ou égal à 0,2. 5. Méthode selon la revendication 1, caractérisée en ce que le log fournissant la portion de référence est un log obtenu par gammagraphie. 6. Méthode selon la revendication 1, caractérisée en ce que chaque portion de log est incluse dans un intervalle d'étude contenant un nombre prédéterminé d'échantillons N. 7. Méthode selon la revendication 6, caractérisée en ce que lorsqu'un log destiné à fournir une portion sélectionnée est constitué par un nombre d'échantillons inférieur au nombre N, ledit log est centré dans l'intervalle d'étude et les parties vides dudit intervalle sont remplies par des échantillons ayant une valeur égale à la valeur moyenne du log. 8. Méthode selon les revendications 6 et 7, caractérisée en ce que lorsqu'un log destiné à fournir une portion sélectionnée est constitué par un nombre d'échantillons supérieur au nombre N, ledit log est divisé en au moins deux parties ayant chacune un nombre d'échantillons inférieur au nombre N, et en ce que chaque partie est traitée selon la revendication 7. 9. Méthode selon la revendication 1, caractérisée en ce que les fréquences spatiales d'analyse sélectionnées pour le calcul de la transformée en ondelettes, ont des fréquences limites correspondant à une longueur d'onde de 4 m et à une longueur d'onde de 200 nim. 10. Méthode selon la revendication 9, caractérisée en ce qu'on sélec- tionne dix firéquences spatiales dont les firéquences limites correspondent à des longueurs d'onde de 10 m et 100 m. 11. Méthode selon la revendication 9 ou 10, caractérisée en ce que la succession des fréquences spatiales d'analyse est une progression géométrique. 12. Méthode selon l'une des revendications I à 11, caractérisée en ce que la fonction ondelette mère est une fonction du type f(z) = (1 - z) exp(-z2/2) dans laquelle z est la profondeur. 13. Méthode de recalage, en profondeur, de logs de natures différentes enregistrés pour une zone d'un milieu en fonction de la profondeur caractérisée en ce qu'elle consiste à déterminer des ruptures sur une portion d'un log choisie comme référence et sur des portions correspondantes des autres logs en appliquant la méthode selon l'une des revendications 1 à 12 et à recalculer les valeurs desdites autres portions de log de telle sorte que les ruptures desdites portions, qui sont associées à chacune des ruptures de la portion de log choisie comme référence, apparaissent à la même cote profondeur que la rupture considérée du log de référence.