Brevet FR2753802A1 — publié le 19980327

INTERPRETATION DE COUCHES SUPERIEURES DE FORMATION La présente invention concerne l'interprétation de couches supérieures de formations à partir de courbes de diagraphie de puits afin d'identifier les caractéristiques identiques parmi les puits. La collecte de données de diagraphie de puits à partir de trous de sondage forés dans la terre fournit des informations qui peuvent être analysées concernant la structure en profondeur de formations de couches sous la surface ou de couches supérieures (par exemple, une limite entre deux formations) dans les terrains pétrolifères et les gisements de gaz. Les informations déduites indiquent le type de roche de la couche sous la surface. L'instrumentation nucléaire, à rayon gamma, électromagnétique, acoustique, magnétique ou d'autres sources est descendue dans les trous de sondage pour générer des signaux de source qui sondent les formations souterraines. Les formations ou les couches supérieures de formations modifient les signaux de source ou réagissent à ceux-ci et des détecteurs sont placés avec les instruments de source dans les trous de sondage pour surveiller les signaux de réponse résultants ou modifiés. Les caractéristiques du signal de réponse, par exemple, son amplitude, varient avec les différents types de signaux de source et dépendront également du type de formation ou de couche supérieure observée. Cés données sont collectées dans le temps et sont appelées, collectivement, "courbes de diagraphie de puits " ou "diagraphie de puits". Les diagraphies de puits sont enregistrées, de manière typique, comme une fonction de la profondeur dans les trous de forage, une courbe enregistrée (ou trace) pour chaque type de source. Plusieurs trous de forage peuvent être utilisés pour collecter des données provenant de puits multiples. Le sondage de puits multiples permet de suivre de manière spatiale les diverses formations ou couches supérieures testées. Des puits, par exemple, des douzaines de puits, peuvent être forés à des miles les uns des autres pour cette analyse. Des milliers de puits peuvent être forés pour des études plus détaillées, espacés d'un quart à un demi-mile les uns des autres. Le forage à une profondeur typique de 10000 pieds permet d'observer de multiples formations dans les traces de réponse. Le défi des spécialistes des sciences de la terre est de cartographier l'emplacement des formations de couches sous la surface à partir de nombreuses diagraphies. Cela peut impliquer de tenter de faire correspondre les caractéristiques dans des traces enregistrées à partir d'une diagraphie de puits aux caractéristiques similaires dans les traces d'autres (par exemple, beaucoup d'autres) diagraphies, un procédé appelé corrélation. L'interprétation des couches supérieures dans les terrains pétrolifères et les gisements de gaz peut, cependant, être une activité difficile et de longue durée. EXPOSE DE L'INVENTION L'invention peut présenter l'avantage d'améliorer la précision de la sélection d'événements fins dans l'interprétation des sommets des couches ou des couches supérieures. L'invention peut également présenter l'avantage de faire gagner beaucoup de temps aux spécialistes des sciences de la terre lors de l'interprétation de sommets de couches ou de couches supérieures. L'invention peut, de plus, présenter l'avantage d'apporter une fiabilité accrue à l'interprétation des sommets des couches ou des couches supérieures. En général, selon un premier aspect, l'invention concerne un procédé d'interprétation des courbes de diagraphie de référence et cibles enregistrées à partir de différents puits pour aider à déterminer la profondeur approximative de caractéristiques correspondantes aux emplacements de puits. Le procédé peut comprendre : le prétraitement des courbes de diagraphie de référence et cibles selon des paramètres entrés par l'utilisateur, l'application d'une fonction de guidage aux courbes qui limite l'interprétation à des parties des courbes, le calcul d'une indication de correspondance entre les parties et l'évaluation de l'indication de correspondance pour déterminer si des événements sur les parties sont liés. En général, selon un autre aspect, l'invention concerne un appareil pour interpréter les courbes de diagraphie de référence et cibles enregistrées à partir de différents puits pour aider à déterminer la profondeur approximative de caractéristiques correspondantes aux emplacements de puits. L'appareil peut comprendre : une mémoire qui stocke les instructions pour interpréter les données de diagraphie de puits et un processeur couplé à la mémoire qui est chargé par les instructions : de prétraiter les données de diagraphie de référence et cibles selon des paramètres entrés par l'utilisateur, d'appliquer une fonction de guidage aux données de diagraphie de puits préretraitées qui limite l'interprétation à des parties des données, de calculer une indication de correspondance entre les parties et d'évaluer l'indication de correspondance pour déterminer si des événements sur les parties sont liés. En général, selon un autre aspect, l'invention concerne un système pour analyser les données de courbes de diagraphie enregistrées à partir des puits répartis géographiquement pour déterminer l'emplacement des sommets correspondants de couches dans les puits. Le système peut comprendre : des dispositifs d'entrée/sortie (E/S) qui reçoivent les données de diagraphie enregistrées à partir de puits et un processeur couplé aux dispositifs d'entrée/sortie (E/S) qui est chargé : de prétraiter les données de diagraphie de puits selon des paramètres entrés par l'utilisateur, de déterminer si une fonction de guidage doit être appliquée aux données, de déterminer si les événements dans les données sont corrélés, en se basant en partie sur les paramètres entrés par l'utilisateur et, si les événements sont liés, de convertir un emplacement associé à la corrélation en un emplacement approximatif des sommets des couches dans les puits. En général, selon un autre aspect, l'invention concerne un procédé d'évaluation de la correspondance entre des événements sur des courbes de diagraphie de référence et cible enregistrées dans les puits de référence et cibles pour identifier un emplacement approximatif de couches connexes dans les puits. Le procédé peut comprendre : la détermination du fait qu'un indicateur de correspondance entre les événements sur les diagraphies de référence et cible a ou non une valeur supérieure à un paramètre de coupure, l'application d'un guide géologique prescrit pour déterminer l'emplacement d'un extremum de l'indicateur de correspondance sur la courbe cible de diagraphie de puits et la conversion de l'emplacement de l'extremum en un emplacement absolu approximatif à l'intérieur du puits cible. En général, selon un autre aspect, l'invention concerne un programme informatique, résidant sur un support susceptible d'être lu par un ordinateur. Le programme informatique peut comprendre des instructions pour entraîner un ordinateur à : prétraiter des données de diagraphie de référence et cible selon des paramètres entrés par l'utilisateur, appliquer une fonction de guidage aux données de diagraphie qui dirige le traitement des données de diagraphie, calculer une indication de correspondance entre des parties des données de diagraphie de puits et évaluer l'indication de correspondance pour déterminer si les parties sont liées. D'autres caractéristiques et avantages de l'invention seront évidents à partir de la description qui suit, y compris les dessins, et à partir des revendications. BREVE DESCRIPTION DES FIGURES La figure 1 est un schéma fonctionnel d'un système informatique selon un mode de réalisation de l'invention. La figure 2 montre un résultat de calcul fonctionnel selon un mode de réalisation de l'invention. La figure 3 montre une courbe de diagraphie de référence et une courbe de diagraphie cible selon un mode de réalisation de I'invenion. La figure 4 montre une mémoire comportant des blocs d'instructions et de données selon un mode de réalisation de l'invention. La figure 5 montre où sont attendues les positions des maxima du calcul fonctionnel sur une courbe de diagraphie cible par rapport à une courbe de diagraphie de référence pour différentes règles selon un mode de réalisation de l'invention. La figure 6 montre une partie de règle de la figure 4 plus en détail selon un mode de réalisation de l'invention. La figure 7 montre l'application d'une des règles selon le mode de réalisation de l'invention. La figure 8 montre comment les sélections faites par un utilisateur ont été affinées à un emplacement plus précis selon un mode de réalisation de l'invention. La figure 9 montre un exemple de la manière d'appliquer un paramètre d'utilisateur selon un mode de réalisation de l'invention. La figure 10 montre l'application d'une procédure à un ensemble de courbes de coupes de sonde selon un mode de réalisation de l'invention. La figure 11 est un organigramme de la procédure de la figure 10 selon un mode de réalisation de 1' invention. La figure 12 montre l'application d'une autre procédure à un ensemble de courbes de diagraphie selon un mode de réalisation de l'invention. La figure 13A est un organigramme d'un procédé selon un mode de réalisation de l'invention. Les figures 13B-E sont des organigrammes plus détaillés du procédé de la figure 13A selon un mode de réalisation de l'invention. La figure 14 est un organigramme d'un autre procédé selon un mode de réalisation de l'invention. La figure 15 montre les effets d'un paramètre de taille trop petit sur une diagraphie de référence selon un mode de réalisation de l'invention. La figure 16 montre les effets d'un paramètre de taille (de la figure 15) trop large sur une diagraphie de référence selon un mode de réalisation de l'invention. La figure 17 montre les effets de la dimension des paramètres de courbe de diagraphie de puits de référence de la figure 15 sélectionnée avec soin selon un mode de réalisation de l'invention. La figure 18 montre comment un procédé est appliqué selon un mode de réalisation de l'invention. La figure 19 montre le résultat de l'application du procédé de la figure 18 selon un mode de réalisation de l'invention. La figure 20 montre le résultat de l'application d'un autre procédé selon un mode de réalisation de 1 'invention. DESCRIPTION DETAILLEE DE MODES DE REALISATION L'invention peut être utilisée pour interpréter des corrélations de sommets de couches géologiques (qui échantillonnent des horizons géologiques) entre des courbes de diagraphie échantillonnées. Avec référence à la figure 1, un système informatique 10 est illustré selon un mode de réalisation de l'invention. Le système informatique 10 peut être un ordinateur personnel (PC), une station de travail, une unité centrale, etc. Des exemples de stations de travail possibles comprennent une station de travail Silicon Graphics Indigo 2 et une station de travail Sun SPARC. Le système informatique 10 stocke et exécute des instructions utilisées pour interpréter des corrélations de sommets de couches entre deux courbes de diagraphie échantillonnées ou plus. Le système informatique 10 comprend un processeur programmable 12 et une mémoire 14 (par exemple, un support susceptible d'être lu par ordinateur) couplés à un bus de système 16. Le processeur 12 peut être, par exemple, un microprocesseur, un micro-contrôleur, ou un processeur d'unité centrale ou de station de travail. La mémoire 14 peut être, par exemple, un disque dur, une mémoire morte (ROM), un disque optique compact (CD ROM), une mémoire vive dynamique (DRAM) ou une autre mémoire vive (RAM), une mémoire flash, un élément de stockage magnétique, un élément de stockage optique, des registres, ou d'autres mémoires volatiles et/ou non volatiles. Un algorithme de corrélation à marqueur automatique (AMCA) comporte des instructions qui peuvent être stockées dans une zone ou un bloc de mémoire AMCA 14a dans la mémoire 14. Ces instructions sont utilisées pour charger le processeur 12 de déterminer une corrélation (un type d'indication de correspondance ou indicateur) entre les courbes de diagraphie échantillonnées. Les données de diagraphie peuvent être reçues à travers un bloc de données reçues 18 qui peut comprendre, par exemple, des dispositifs d'entrée/sortie, tels qu'un clavier d'utilisateur (pour les paramètres entrés par l'utilisateur), un modem, un téléphone, un adaptateur de réseau intégré de service numérique (ISDN), des dispositifs ou systèmes sans fil, ou d'autres lignes ou services de communication numérique ou analogique quelconques. Les données de diagraphie qui sont reçues (qui peuvent être déjà numérisées (échantillonnées) ou converties d'analogique en numérique [non montré] lors de leur réception) par le système 10 peuvent être acheminées depuis le bloc de données reçues 18, qui est couplé au bus de système 16, vers une zone 14b de données de diagraphie dans la mémoire 14 à travers le bus de système 16 ou, alternativement, peuvent déjà être stockées dans la zone 14b de données de diagraphie. Avec référence à la figure 2, les instructions stockées dans la zone de mémoire 14a chargent le processeur 12 de générer la fonction AMCA, Cmax(X) 30, entre une courbe de diagraphie de référence et une ou plusieurs courbes cibles qui peuvent être stockées dans la zone 14b de données de diagraphie. La fonction AMCA 30 peut être stockée dans la mémoire 14 dans le bloc de mémoire 14c. La fonction AMCA 30 est utilisée pour faire correspondre (corréler) les événements connexes enregistrés sur les courbes de diagraphie de référence et cibles pour les mêmes formations géologiques ou couches supérieures qui échantillonnent des horizons géologiques. Elle quantifie la qualité atteinte d'une correspondance et produit une valeur numérique qui peut être interprétée comme une mesure de la fiabilité de la correspondance. Un nombre de "fiabilité" pourrait être attribué sur la base de la corrélation produite par la fonction AMCA 30. La fonction AMCA 30 peut être déterminée pour n'importe quel nombre de puits et n'importe quel nombre d'événements dans chaque puits. La fonction AMCA 30 peut être affichée sur un support d'enregistrement de sortie 20 couplé au bus de système 16. Le support d'enregistrement de sortie 20 peut être, par exemple, un écran d'ordinateur, un enregistreur de diagrammes, un traceur, un élément de stockage à mémoire tel que décrit ci-dessus, ou tout autre format d'enregistrement de sortie tangible. Des exemples d'une courbe de diagraphie de référence échantillonnée x 21 et d'une courbe de diagraphie cible échantillonnée y 23 pour calculer la fonction AMCA 30, toutes deux enregistrées comme une fonction de la profondeur du sondage 22 sous la surface de la terre, sont montrées sur la figure 3. Une fenêtre 26 de longueur w et une fenêtre 28 de longueur Wy peuvent être définies, respectivement, sur les données de courbe de diagraphie de puits x 21 et de diagraphie de puits y 23. La fonction AMCA 30 entre les deux courbes de diagraphie échantillonnées 21 et 23 est calculée en faisant glisser la fenêtre 26 de manière incrémentale (illustrée de manière schématique sur la figure 3) sur la fenêtre 28 et en additionnant les termes selon pour la position de corrélation T = Wy - w/2 à Wy + w/2, où CT est la fonction de covariance Dans l'équation de CT, les termes Xw et Yw sont les valeurs moyennes et les termes vxw et oy, sont les écarts types, des données dans les fenêtres 26 et 28, des courbes de diagraphie échantillonnées x 21 et y 23, respectivement, pour chaque T. Ca/Gb est le rapport de ces écarts types calculé à chaque position T et Sa/Eb est le rapport des sommes des amplitudes des courbes de diagraphie de référence et cibles (par exemple, la diagraphie x 21 et la diagraphie y 23 de la figure 3) dans les fenêtres 26 et 28 également calculé à chaque position de corrélation 1. Bien que cela ne soit pas nécessaire, il est commode de maintenir 0 Pour cela, pour chaque calcul de la fonction AMCA entre une courbe de coupe de diagraphie de référence et une courbe de coupe de diagraphie cible : (1) le terme d'écart type "b" est déduit de n'importe laquelle des deux courbes qui a la plus grande valeur d'écart type et le terme d'écart type "a" est déduit de l'autre courbe. (Cela est fait de sorte que le rapport d'écart type soit inférieur ou égal à un et supérieur ou égal à zéro (une approche similaire est utilisée pour déterminer le rapport de somme d'amplitude de "a" et "b") et (2) les paramètres a, ss et y sont ajustés de sorte que a + + t r = 1. La fonction AMCA 30 est exprimée comme une fonction de la position de corrélation T et a une longueur Wy+w comme montré sur les figures 2 et 3. Les longueurs w et Wy (c' est-à- dire, la quantité de données échantillonnées) peuvent toutes deux être sélectionnées au choix de l'utilisateur pour opérer le calcul AMCA entre les courbes de référence et cibles. Les termes xi et Y sont les données d'amplitude échantillonnées provenant des courbes de coupes de diagraphie de référence et cibles, respectivement. L'index i va de -w/2 à w/2 et est incrémenté, pour exécuter l'addition, de la fréquence d'échantillonnage, des courbes de diagraphie échantillonnées, par exemple, 6 pouces ou 1 pied. La position de corrélation T a une profondeur intégrée et représente une différence de profondeur entre les courbes de diagraphie de référence et cibles. Parce que la position de corrélation T a une profondeur intégrée, la profondeur absolue approximative par rapport à la surface de la terre est connue et donnée par T + une profondeur constante connue (par exemple, connue lorsque les courbes de diagraphie sont enregistrées). La position de corrélation T est également incrémentée de la fréquence d'échantillonnage. Avec référence à la figure 4, selon un mode de réalisation de l'invention, des règles géologiques 50 peuvent être utilisées pour établir des paramètres pour le calcul AMCA et pour analyser la fonction AMCA 30 résultante pour une courbe de référence et une ou plusieurs courbes cibles. Les règles géologiques 50 peuvent être stockées comme des instructions pour le processeur 12 dans la zone de mémoire 14d dans la mémoire 14. Les règles géologiques 50 peuvent comprendre des horizons de guidage 50A, des règles de séquence 50B et une fonction de guidage par l'utilisateur 50C dont les codes correspondants peuvent être stockés dans la zone de mémoire 14d correspondante de la mémoire 14, comme montré sur la figure 4. Avec référence à la figure 5, des horizons de guidage 50A peuvent être définis ou appliqués comme des limites (par exemple, les horizons 52A et 52B) sur les courbes de diagraphie échantillonnées 42 et 44 (identiques aux courbes de diagraphie 21 et 23, respectivement) qui peuvent être utilisées pour générer la fonction AMCA 30. Les horizons 50A (52A, 52B) sont utilisés pour limiter l'étendue des courbes de diagraphie échantillonnées 42, 44 sur laquelle la fonction AMCA 30 doit être calculée. Ces limites sont connues comme des horizons de guidage parce qu'elles identifient ce que l'utilisateur sait ou croit être la même couche supérieure (limite ou horizon échantillonné) sur chacune des courbes. La connaissance ou la conviction de l'utilisateur est basée sur d'autres informations connues de l'utilisateur concernant la formation ou les données de courbe de diagraphie. Par exemple, on peut savoir que les signaux enregistrés de sommet de couche de formation 53A et 53B identifiés sur la courbe de diagraphie de référence 42 sont liés aux signaux enregistrés de sommets de couche 55A et 55B (c'est-à-dire, qu'ils sont déduits des mêmes sommets de la formation), respectivement , sur la courbe de diagraphie cible 44 sur la base d'autres données régionales ou géologiques (par exemple, sismiques préexistantes) A la suite d'une procédure automatique selon un mode de réalisation de l'invention, l'utilisateur sélectionne (sélectionne ou marque) les courbes échantillonnées (c'est-à-dire, sur les courbes de diagraphie de référence et cibles) pour définir les horizons de guidage (par exemple, avec un dispositif d'entrée, tel qu'une souris, couplé au support d'enregistrement de sortie 20). On doit comprendre que toute sélection ou tout marquage, tel que décrit ici, peut impliquer l'identification d'emplacements ou d'adresses de mémoire spécifiques dans la mémoire 14 associés à des données particulières des courbes échantillonnées. Le support de sortie 20 peut être, par exemple, un écran qui montre les courbes de diagraphie et les horizons de guidage sélectionnés par l'utilisateur sur celles-ci. Les horizons de limite 50A (52A, 52B) peuvent également être sélectionnés à partir d'interprétations géologiques dont l'utilisateur a connaissance et, également, à partir de réseau cartographique de ces interprétations. Alternativement, avant de commencer le calcul de la fonction AMCA 30, les horizons de guidage peuvent être marqués automatiquement à partir de données stockées au préalable dans les horizons de guidage 50A du bloc de mémoire 14b de la mémoire 14. Lorsque aucun horizon de guidage n'est prévu par l'utilisateur ou marqué automatiquement pour le calcul AMCA, les horizons sont supposés être plats. Les horizons de guidage plats sont alignés à la même profondeur à travers les courbes de diagraphie de référence et cibles similaires à l'horizon de guidage 52A montré sur la figure 5. Avec référence à la figure 6, les règles de séquence 50B des règles géologiques 50 qui peuvent être utilisées pour aider au calcul de la fonction AMCA 30 comprennent, par exemple, des définitions de séquence superposée 54A, de séquence tronquée 543, de séquence conforme 54C et de séquence non structurée 54D. Le code correspondant pour ces règles de séquence peut être stocké dans la zone de règles de séquence 50B du bloc de mémoire 14d de la mémoire 14, comme montré sur la figure 4. Les règles de séquence 50B aident à limiter l'étendue générale sur la courbe de diagraphie cible sur laquelle la fonction AMCA 30 est calculée lorsqu'on cherche une correspondance à un événement choisi par l'utilisateur sur la courbe de diagraphie de référence (par exemple, sur la courbe 42). Les règles de séquence 50B peuvent être choisies par l'utilisateur. Lorsqu'elles sont appliquées, elles fournissent une indication de la position supposée, sur la courbe diagraphie cible (par exemple, sur la courbe 44), de l'emplacement de l'extremum ou du maximum de la corrélation (par exemple, le maximum de la fonction AMCA 30 identique au maximum 30A de la figure 3). Avec les règles de séquence 50B, la fenêtre 28 définie sur la courbe de diagraphie cible est débutée (c'est-à-dire que la position (bord) initiale de la fenêtre 28 est déterminée) en fonction des horizons de guidage, avec la longueur Wy de la fenêtre définie par l'utilisateur et ne s'étend pas au-delà des limites des horizons de guidage 50A (52A, 52B). L'utilisateur peut choisir la règle de séquence 50B pour chaque intervalle de diagraphie limité par deux horizons de guidage 50A (52A, 52B) à travers les courbes de diagraphie 42 et 44 montrées sur la figure 5. Le choix peut être réalisé parmi ce qui suit 1. Séquence superposée 54A - la position maximale de corrélation de la fonction AMCA 30 sur la courbe de diagraphie cible est supposée être parallèle et en dessous de l'horizon de limite supérieure (par exemple, l'horizon 52A) à la même position de profondeur (par exemple, la position 62 sur la figure 5) que le marqueur sélectionné par l'utilisateur (par exemple, le marqueur A 60 sur la figure 5) à travers les courbes 42 et 44 (les sélections ou les marqueurs peuvent etre montrés sur les dessins sur les courbes de diagraphie ou déplacés par rapport aux courbes, comme dans le cas du marqueur 60). Dans l'un ou l'autre cas, on doit comprendre que chaque marqueur marque un emplacement spécifique sur la courbe correspondant à un emplacement d'adresse de données de courbe dans la mémoire 14b, où, par exemple, la ligne en traits pointillés croise la courbe de référence, comme sur la figure 5, ou à la même position de profondeur que le marqueur même si une ligne en traits pointillés peut ne pas être illustrée); 2. Séquence tronquée 54B - la position maximale de corrélation de la fonction AMCA 30 sur la courbe de diagraphie cible est supposée être parallèle et audessus de l'horizon de limite inférieure (par exemple, l'horizon 52B) à la meme position de profondeur relative (par exemple, la position 64) par rapport à l'horizon de limite inférieure que le marqueur sélectionné par l'utilisateur (par exemple, 60) à travers les courbes 42 et 44 3. Séquence conforme 54C (examinée plus en détail ci-dessous) - la position maximale de corrélation de la fonction AMCA 30 sur la courbe de diagraphie cible est supposée être décalée vers une position (par exemple, la position 66) qui est proportionnelle au changement d'épaisseur de la couche (changement d'épaisseur de la formation tel qu'enregistré par des événements sur les courbes de diagraphie de puits) en allant de la courbe de diagraphie de référence vers la courbe de diagraphie cible pour chaque point d'évaluation sur la courbe de diagraphie cible (c'est-à-dire que chaque couche évaluée pour la corrélation peut présenter une épaisseur différente sur la courbe de diagraphie cible de son épaisseur sur la courbe de diagraphie de référence) ; et 4. Séquence non structurée 54D - la position maximale de corrélation de la fonction AMCA 30 sur la courbe de diagraphie cible est supposée être n'importe où dans l'intervalle défini Wy sur la courbe de diagraphie cible (voir, par exemple, la figure 3). La figure 7 fournit plus de détails à propos de la séquence conforme 54C. Pour chaque marqueur sélectionné par l'utilisateur A 60 dans la fenêtre 68 (similaire à la fenêtre 26), la fonction AMCA 30 peut être calculée avec la fenêtre 68 de longueur w de la courbe de diagraphie de référence x 42 d'abord étirée (ou compressée) proportionnellement (par exemple, en utilisant une fonction d'interpolation en sinx/x) pour faire correspondre l'épaisseur (minceur) d'un intervalle correspondant 70 sur la courbe de diagraphie de cible d'évaluation y 44. La figure 7 montre schématiquement une section (c'est-à-dire, un intervalle) 43 dans la fenêtre 68 de la courbe de référence 42. La section 43 doit être évaluée en calculant la fonction AMCA 30 et est montrée en cours d'étirement proportionnellement à sa dimension (w) dans la fenêtre 68 sur la courbe 42 pour correspondre à l'intervalle de dimension proportionnel correspondant de la fenêtre 70 sur la courbe 44. L'étirement (ou la compression) est réalisé avant le calcul de la fonction AMCA 30. De nouveau avec référence à la figure 4, une procédure ou une règle 50 supplémentaire, selon un mode de réalisation de l'invention, peut être appliquée comme une fonction de guidage de l'utilisateur ,(également appelée procédé semi-automatique) 50C dont le code peut être stocké dans le bloc de mémoire 14d de la mémoire 14. La fonction de guidage par l'utilisateur 50C peut être efficace, par exemple, dans l'environnement de station de travail, lorsque : (1) les courbes de diagraphie (par exemple, 42, 44) sont compressées à des fins d'affichage (par exemple, pour afficher une courbe de diagraphie de référence et un certain nombre de courbes de diagraphie cibles) et (2) un interprétateur (utilisateur) peut facilement réaliser des corrélations visuelles entre des puits, mais ne peut pas marquer la profondeur de manière précise (par exemple, sur l'échelle de profondeur 22 de la figure 3) sur la courbe de diagraphie cible à l'échelle de visualisation de l'écran. Dans de telles situations, des instructions sont prévues dans le bloc de mémoire 14d qui donnent l'ordre au processeur 12 de permettre que la fonction AMCA 30 soit plutôt guidée par la sélection "sommaire" de l'utilisateur sur la courbe de diagraphie cible pour trouver une profondeur exacte pour un maximum (corrélation) (par exemple, 30A de la figure 2) de la fonction AMCA 30. La fonction de guidage par l'utilisateur 50C peut permettre une sélection rapide de couches supérieures de formations tout en permettant à l'interprétateur un affichage régional (à grande échelle) des données. Par exemple, l'utilisateur doit d'abord sélectionner et marquer un événement sur une courbe de diagraphie de référence affichée (par exemple, affichée sur le support d'enregistrement de sortie 20) avec le degré de précision souhaité, mais peut ne pas être capable de le faire si facilement sur la courbe de diagraphie cible lorsque la courbe de diagraphie cible est affichée en une petite dimension (par exemple, lorsqu'elle est affichée dans une vue compressée ou sans zoom) . En utilisant la fonction de guidage de l'utilisateur 50C, l'utilisateur doit seulement sélectionner (par exemple, avec une souris ou un clavier connecté au système 10) un emplacement approximatif sur la courbe de diagraphie cible pour évaluation par le calcul de la fonction AMCA 30. Avec la fonction de guidage par l'utilisateur 50C, le maximum de corrélation est supposé être proche du point approximatif sélectionné par l'utilisateur sur les courbes cibles. Ensuite, le système 10 trouvera le maximum de corrélation (fonction AMCA 30) qui dépasse le paramètre de coupure AMCA défini par l'utilisateur (examiné ci-dessous) qui est le plus proche de la sélection de l'utilisateur pour la profondeur exacte. Le procédé semi-automatique permet à l'utilisateur de sélectionner de nombreux trous de sondage et de nombreux événements à une échelle réduite (c' est-à- dire, sans zoom), avec une faible définition, avec tous les trous de sondage et tous les événements présentés sur l'écran en même temps. Cela peut épargner beaucoup de temps de repositionnement, de réalisation de zoom, de sortie du mode zoom et de défilement des données. Un affichage compact d'une pseudo-coupe (par exemple, sur le support d'enregistrement de sortie 20 qui peut être un écran d'ordinateur) peut être utilisé pour effectuer la sélection sur un groupe de courbes de trou de sondage cibles, typiquement jusqu'à 40 sur un écran. La figure 18 (examinée ci-dessous) montre un exemple d'un événement (indicatif d'un sommet de couche) sélectionné sur la diagraphie de référence en utilisant le procédé de corrélation semi-automatique. Une fois que l'événement 94a est sélectionné, les événements correspondants 94b-g dans les courbes cibles sont trouvés et marqués (une fois que l'utilisateur a effectué ses sélections sommaires sur les courbes de diagraphie cibles) . La figure 18 montre que même lorsqu'on utilise une très petite échelle, il est possible de sélectionner des profondeurs précises en se basant sur l'événement de référence sélectionné. Sur la figure 8, la ligne 98 montre comment les sélections de commande (par exemple, la sélection de référence 98a sur la courbe de diagraphie de référence 21 et la sélection de l'utilisateur 98b sur la courbe de diagraphie cible 23) effectuées par l'utilisateur ont été affinées pour l'emplacement précis des marqueurs 98c sur les courbes de diagraphie cibles. La courbe cible 23 sur la figure 8 apparaît en zoom ou agrandie pour les besoins de l'illustration du procédé semiautomatique. Un utilisateur peut, cependant, choisir d'afficher la courbe 23 sur l'écran d'ordinateur soit en zoom (comme sur la figure 8) avec peu d'autres courbes affichées, soit sans zoom avec beaucoup d'autres courbes affichées, selon ses besoins. De nouveau avec référence à la figure 4, les règles géologiques 50 sont montrées comprenant également un paramètre de coupure AMCA 50D, dont le code est stocké dans le bloc de mémoire 14d de la mémoire 14. Le paramètre de coupure 50D est un critère, susceptible d'etre sélectionné par l'utilisateur, d'acceptation du fait qu'une correspondance valide a été trouvée entre les courbes de diagraphie de référence et cible (par exemple, 42 et 44). Autrement dit, c'est une valeur de la fonction AMCA 30 à laquelle, ou au-dessus de laquelle, l'utilisateur est prêt à accepter qu'une corrélation valide a été trouvée. Le paramètre de coupure 50D (valeur) est exprimé comme une valeur en pourcentage, avec une correspondance identique représentée par 100 % (c'est-à-dire, une correspondance identique a une valeur de 1 pour le maximum de corrélation). Lors de l'exécution du calcul AMCA, il peut être souhaitable de fixer le paramètre de coupure 50D à une valeur élevée (par exemple, 0,7 - 0,8) afin d'assurer une fiabilité améliorée que des correspondances ont été trouvées (entre des événements sur les courbes de diagraphie) pour des sommets de couches connexes. Un exemple de la manière d'appliquer le paramètre de coupure 50D est montré sur la figure 9. Il peut être nécessaire d'ajuster la valeur du paramètre de coupure 50D fixé par l'utilisateur, en fonction de la qualité des données de courbe diagraphie. Par exemple, la valeur initiale 31 (ligne A de la figure 9) peut être fixée à une valeur élevée (par exemple, à 0,8) pour des données de puits cibles plus près du puits de référence et abaissée, ensuite, à une valeur 32 (ligne B de la figure 9) pour des courbes de diagraphie cibles de puits cibles plus éloignés du puits de référence. Il peut etre nécessaire ou avantageux de s'adapter, par exemple, à tous les changements qui peuvent survenir dans la même formation à des distances éloignées du puits de référence. Quelques exemples de changements pouvant survenir seront examinés ci-dessous. Pour les puits cibles qui sont plus éloignés, la fonction AMCA 30 peut apparaltre comme sur la figure 9 par rapport aux valeurs de paramètre de coupure 31 et 32. Avec référence de nouveau à la figure 4, des géologies complexes 50E, telles qu'un croisement, des sections répétées et des sections inversées, peuvent etre adaptées ou traitées par l'ensemble de règles 50 pour le calcul de la fonction AMCA 30. Si l'utilisateur suppose que des géologies complexes doivent être examinées, il ou elle peut positionner des indicateurs appropriés. Les géologies complexes 50E peuvent être stockées dans le bloc de mémoire 14d de la mémoire 14. Un croisement se produit lorsque des couches de formation successives inversent leur ordre de succession dans une région. Sans croisement, pour calculer la fonction AMCA 30, on suppose que les strates (formations) entre les horizons de guidage (profondeurs) sont déposées dans le même ordre aux emplacements de puits de référence et de puits cible. De là, un croisement n'est pas autorisé. Si deux maximums de corrélation successifs dans la fonction AMCA réalisent un croisement (signifiant que deux couches successives se sont croisées), cependant, celle avec la valeur maximale de corrélation la plus élevée de Cmax(z) est sélectionnée comme l'événement correct. L'autre maximum de corrélation est indiqué par un point d'interrogation ("?") pour indiquer à l'utilisateur un conflit possible ou une "erreur de sélection" (expliquée ci-dessous) dans l'interpolation. Des sections répétées sont simplement là où des couches de formation successives se répètent en profondeur. Une section répétée peut se produire dans des zones de faille inverse. Si des sections répétées sont supposées, elles peuvent être traitées sous la forme d'une multiplication de surfaces définies par marqueurs (ou horizons) définies en profondeur. Dans ce cas, le calcul de la fonction AMCA 30 sélectionnera plusieurs maximums de corrélation pour le nombre de sections répétées supposé. Des sections inversées ou renversées (c'est-à-dire que la structure en profondeur des sections est physiquement inversée) peuvent survenir dans des zones de plis couchés. Si des sections inversées sont supposées, elles peuvent être traitées en autorisant des corrélations inversées entre les régions limitées des courbes de diagraphie de référence et cibles. Avec les corrélations inversées, la section approximative d'un des ensembles des données de diagraphie de référence ou cible est renversée (inversée) et l'autre n'est pas inversée avant que le calcul de la fonction AMCA 30 soit exécuté. Pour le test de cohérence, une corrélation de retour (calculs de la fonction AMCA 30 "de retour") 56 (pour des coupes de strates de formation) et une fermeture de boucle (calculs de fermeture de boucle de la fonction AMCA 30) 58 (pour l'interprétation de carte de strates de formation régionales), dont les codes sont stockés dans le bloc de mémoire 14d de la mémoire 14 (figure 4), peuvent être exécutés pour validation. Par exemple, un calcul de corrélation de retour peut être exécuté et peut informer l'utilisateur de sélections de sommets de couches de formation incompatibles. On se réfère à la figure 10 qui illustre schématiquement comment la corrélation de retour est exécutée et à la figure 11 qui est un organigramme d'un procédé selon un mode de réalisation de l'invention pour exécuter la corrélation de retour. Une courbe de diagraphie de référence 420 déterminée auparavant (par exemple, une partie d'une courbe de diagraphie comme la courbe 21 de la figure 3 que l'utilisateur considère comme une bonne courbe de diagraphie de puits) est sélectionnée (252) comme la courbe de diagraphie de référence initiale. Les autres courbes de diagraphie 422, 424 et 426 sont prises à partir de puits forés le long d'une direction linéaire (ou n'importe quelle direction dans le cas de la fermeture de boucle, tel que décrit ci-dessous) du puits à partir duquel la courbe 420 a été déduite. La courbe de diagraphie cible est choisie (254) pour être la courbe de diagraphie 422 et un calcul de la fonction AMCA 30 (corrélation) est exécuté (256) entre les courbes 420 et 422. Ensuite, la courbe 422 est choisie (258) pour être la nouvelle courbe de diagraphie de référence et un nouveau calcul de corrélation est exécuté (260) entre la courbe 422 et la courbe 424. Le procédé est répété jusqu'à ce que la dernière courbe de diagraphie 426 soit atteinte (262). Ensuite, le programme inverse la direction (264), la courbe 426 est choisie (266) comme la nouvelle courbe de diagraphie de référence initiale et la courbe 424 est choisie comme la nouvelle courbe de diagraphie cible initiale. La corrélation (calcul de la fonction AMCA 30) entre les courbes 426 et 424 est déterminée (268), et ainsi de suite, jusqu'à ce que la dernière corrélation (calcul de la fonction AMCA 30) soit déterminée entre les courbes 422 et 420. Si le résultat de la dernière corrélation correspond (270) au résultat du calcul de corrélation initial (tout premier), alors, tous les résultats de corrélation sont supposés être valides, signifiant que les événements connexes sur les courbes de diagraphie ont été identifiés correctement (272). Autrement, l'identification n'est pas valide (274). Un procédé de validation de corrélation similaire pourrait être décrit pour la fermeture de boucle qui est illustré schématiquement sur la figure 12 en commençant avec la partie de courbe de diagraphie 320 et en poursuivant par les parties de courbes de diagraphie 322, 324 et 326 et, ensuite, de retour vers la partie 320. Avec la fermeture de boucle, les courbes de diagraphie sont déduites de puits qui sont répartis de manière spatiale (c'est-à-dire que les puits ne sont pas limités à être forés décalés linéairement les uns des autres le long de la même direction) en deux ou trois dimensions. Une description fonctionnelle plus détaillée du fonctionnement du système 10 commandé par le processeur 12 et les instructions stockées dans la mémoire 14 est présentée ci-dessous. Avec référence aux figures 13A-E, un procédé automatique et le procédé semi-automatique selon des modes de réalisation de l'invention sont montrés en tant qu'organigrammes de modules de logiciel. En commençant au bloc de début 102 de la figure 13A, les données d'une courbe de diagraphie de référence (par exemple, 21 sur la figure 3) et d'une courbe de diagraphie cible (par exemple, 23 sur la figure 3) sont prétraitées (104). Les instructions 35 pour le prétraitement (104) peuvent etre stockées dans le bloc de mémoire 14c de la mémoire 14 et exécutées afin d'améliorer la recherche de détermination d'une bonne corrélation entre les événements sur les courbes de diagraphie de référence et cibles. Le prétraitement (104) peut comprendre le rétrécissement de la fenêtre, la correspondance spectrale et un procédé de dérivation dont les instructions sont stockées dans les blocs 36, 37 et 38, respectivement, des instructions 35 dans le bloc de mémoire 14c de la mémoire 14. Les instructions 36 pour le rétrécissement de fenêtre peuvent impliquer une "pondération centrale", par exemple, l'application d'une fonction de pondération triangulaire bien connue aux données sur la diagraphie de référence et sur la diagraphie cible, comme cela sera examiné davantage cidessous. La pondération centrale est utilisée pour centrer les données de diagraphie entre les diagraphie de référence et cibles de préférence vers le centre des fenêtres respectives (par exemple, les fenêtres 26 et 28) avant de calculer la fonction AMCA 30. Par exemple, avec la pondération triangulaire, pour chaque fenêtre 26 ou 28 de la courbe de diagraphie de référence ou de la courbe de diagraphie cible, respectivement, en commençant à la valeur de données d'amplitude de centre (non montrée), chaque valeur de données d'amplitude est multipliée par un nombre qui décale toutes les valeurs vers le centre de la fenêtre. En particulier, la valeur de centre de la fenêtre 26 (ou 28) peut etre multipliée par 1 et les valeurs aux bords de la fenêtre 26 (ou 28) peuvent être multipliées par 0, avec les autres valeurs entre le centre et les bords multipliées par des valeurs proportionnelles à la pente des parties de rampe de la fonction de pondération triangulaire (non montré). Plus la fenêtre 26 ou 28 est grande (c' est-à- dire, plus la valeur de w ou de Wy est grande), plus la pente (parties de rampe de) de la fonction de pondération triangulaire est grande, et vice versa. Les instructions 37 pour la correspondance spectrale impliquent la transformation des données de diagraphie de référence et cible dans le domaine des fréquences pour filtrage (par exemple, filtrage de limitation de fréquences). Le filtrage est une procédure pour enlever le bruit des données afin d'améliorer la probabilité de détermination d'une bonne corrélation. Les instructions 38 pour le procédé de dérivation impliquent l'exécution du calcul de la fonction AMCA 30 avec les dérivées (c'est-à-dire, après différentiation) des données de diagraphie de référence et cible plutôt qu'avec les données de diagraphie elles-mêmes. Cette procédure peut contribuer à la détermination de bonnes corrélations lorsque des pics de forte amplitude (par rapport au contexte) ou des diagraphies moyennées par blocs sont impliquées. Le calcul de moyenne par blocs est un procédé de lissage, où les données de courbe de diagraphie sont prémoyennées en blocs avant d'être échantillonnées. La figure 13B montre davantage de détails du prétraitement (104). Dans le prétraitement (104), les paramètres de l'utilisateur sont lus (106). Les exemples de paramètres d'utilisateur comprennent les valeurs de a, D et y fixées par l'utilisateur, des valeurs du paramètre de coupure 50D allant de 0,0 à 1,0 et le fait que l'utilisateur souhaite ou non fixer des indicateurs pour une géologie complexe et pour l'ordre de traitement. Une fois que les paramètres de l'utilisateur sont lus, les courbes de diagraphie et les marqueurs (s'il y en a) provenant de la base de données (par exemple, les marqueurs d'horizon de guidage qui sont stockés dans le bloc de mémoire 14b) sont lus (108) . Il est ensuite déterminé si un filtrage des diagraphie aura lieu ou non (110), par exemple, en questionnant l'utilisateur. Si oui, un filtre approprié est appliqué (112), tel que pour la correspondance spectrale (une procédure pour filtrer une courbe de diagraphie de sorte que ses composantes de fréquence correspondent à celles d'une autre courbe de diagraphie), la dérivée d'une courbe de diagraphie, ou la normalisation (une procédure destinée à normaliser les courbes de diagraphie de référence et cibles à une valeur centrale sur les courbes respectives). Sinon, ou si le filtre a été appliqué, il est déterminé si les paramètres de l'utilisateur indiquent que des indicateurs devraient être positionnés pour une géologie complexe 114. Si oui, des indicateurs appropriés de géologie complexe sont positionnés (116) pour des croisements, des sections répétées et/ou des sections inversées. Si l'utilisateur ne voulait pas fixer des indicateurs pour une géologie complexe ou s'ils ont été fixés, il est déterminé si les paramètres de l'utilisateur indiquent que des indicateurs devraient être positionnés pour l'ordre de traitement (118). Si oui, des indicateurs appropriés sont positionnés pour l'ordre de traitement (120). L'ordre de traitement peut être fixé parmi, mais non limité à, la corrélation de retour, la fermeture de boucle et la référence fixe ou "roulante" (décrite ci-dessous). Si des indicateurs d'ordre de traitement ne sont pas fixés ou une fois qu'ils sont fixés, le programme continue par le traitement de fonction de guidage (122) (voir les figures 13A et 13C). Dans le traitement de fonction de guidage (122), la règle géologique (par exemple, provenant des règles géologiques 50) est obtenue pour le marqueur (s'il y en a) (124). La fenêtre cible est ensuite calculée (12fez (par exemple, la dimension Wy et la position de la fenêtre cible 28 sur la courbe de diagraphie cible sont déterminées), laquelle peut être basée sur la règle géologique, s'il y en a une. Cela peut impliquer la détermination du fait que l'utilisateur souhaite lancer ou non l'option automatique (128). Si l'utilisateur veut l'option automatique, la position de début de la fenêtre (Wy) par rapport à l'horizon de guidage (fonction) et la règle de séquence sont déterminées (130). Si, d'autre part, l'utilisateur souhaite actionner l'option semi-automatique (132), la position de début de la fenêtre (Wy) par rapport à une sélection de l'utilisateur sur la courbe cible est déterminée (134). Si la règle de séquence conforme doit être appliquée (136), l'étirement ou la compression proportionnellement à la courbe de diagraphie de référence est exécuté (138). Si la règle de séquence conforme n'est pas appliquée, ou une fois qu'elle est appliquée, il est déterminé si une section inversée est supposée (140). Si c'est le cas, la courbe de diagraphie de référence (ou la section appropriée de la courbe de diagraphie de référence) est inversée (142). Si une courbe ou une section inversée n'est pas supposée, ou une fois que la courbe ou la section de diagraphie de référence est inversée, le programme exécute le traitement AMCA (144) (voir les figures 13A et 13D). Dans le traitement AMCA (144), pour calculer la fonction AMCA 30, un index de boucle est fixé (146) pour la position de corrélation T (c'est-à-dire que l'index de boucle est fixé de manière à correspondre aux positions sur la courbe de diagraphie cible allant de Wy(début)-w/2 à Wy(fin)Tw/2). Si l'option automatique a été choisie (148), la position de début de la fenêtre (Wy) (voir par exemple, les positions 92B sur la figure 18) est fixée (150) par rapport à la fonction de guidage (par exemple, A sur la figure 18) et à la règle de séquence. Si, d'autre part, l'option semi-automatique a été choisie (152), le début de la fenêtre (Wy) (par exemple, le début de la fenêtre (Wy) 98d sur la courbe cible 23 sur la figure 8) par rapport à la sélection de l'utilisateur (par exemple, la sélection de l'utilisateur 98b sur la figure 8) sur la courbe de diagraphie cible (par exemple, la courbe cible 23 sur la figure 8) est fixée (154). (Le début de la fenêtre (w) 98e est le début de la fenêtre sur la courbe de référence 21 sur la figure 8). La pondération centrale, par exemple, la pondération triangulaire (décrite ci-dessus), est ensuite exécutée (156) dans les fenêtres w et Wy au retard courant (T) . Cmax(T) (fonction AMCA 30) est ensuite calculé (158) pour l'alignement courant (T) en : (1) calculant (160) la covariance C#; (2) calculant (162) les écarts types 6w, (c' est-à-dire, 6xw) et 6wy (c' est-à-dire, oyw) ; (3) calculant (164) les sommes des amplitudes dans les fenetres (Ew et Swy) ; et (4) déterminant (166) Cmax(T) (voir l'équation ci-dessus), où ab (axa) est le plus grand (ou le plus petit) de Sw et Swy comme examiné cidessus. Si le dernier T n a pas encore été atteint (168) dans le traitement, le programme retourne à l'étape 156 jusqu'à ce que Cmax(X) ait été déterminé pour tous les T. Ensuite, la fonction AMCA 30 (Cmax(T)) est stockée (170) dans la zone 51 Cmax(n) 30 de la zone de mémoire 14e de la mémoire 14 et le programme exécute le traitement d'évaluation (172) (voir les figures 13A et 13E). Dans le traitement d'évaluation (172), le maximum de la fonction AMCA 30 est trouvé (174) qui satisfait aux règles géologiques et au paramètre de coupure. Le maximum qui est trouvé peut être, par exemple, le maximum le plus proche (non montré) de la fonction AMCA 30 qui correspond à une position (par exemple, la sélection AMCA 98c sur la courbe de diagraphie cible 23 sur la figure 8 dans le procédé semi-automatique) par rapport à la position supposée (voir les figures 5 et 7) ou la sélection de l'utilisateur (par exemple, la sélection de l'utilisateur 98b sur la figure 8 dans le procédé semi-automatique), sur la base de la règle géologique particulière appliquée. Le maximum est également un maximum qui dépasse le paramètre de coupure (valeur). La position du maximum (c'est-à-dire, la position T), par exemple, la valeur 3de T de 30A sur la figure 2, est ensuite convertie (176) en une profondeur absolue approximative sur la courbe de diagraphie cible et, donc, en la profondeur absolue approximative dans le puits cible foré dans lequel la courbe de données cible est enregistrée. La conversion peut être effectuée en prenant la profondeur connue approximative du bord supérieur 33 (voir figure 3), par exemple, de la fenêtre 28 (c'est-à-dire, la position de début Wy) de la diagraphie de puits cible et en lui ajoutant la valeur de profondeur de la position du maximum de la fonction AMCA 30 (par exemple, 30A à 6, comme montré sur la figure 2) qui est connue par rapport au bord supérieur 34 (figure 3). Un nouveau marqueur est ensuite créé (ou un marqueur existant est mis à jour) (178) sur la diagraphie cible (par exemple, stocké comme un emplacement de marqueur 25 dans le bloc 14b de la mémoire 14 pour marquer la courbe de diagraphie cible) à cette profondeur absolue approximative, qui correspond à la profondeur absolue approximative dans le puits cible pour les données cibles marquées. S'il y d'autres diagraphies cibles à traiter, (180), le programme retourne ensuite à l'étape 146 dans le traitement AMCA (144) (figures 13A et 13D) pour la diagraphie cible suivante et la sélection de l'utilisateur suivante (si le procédé semi-automatique est lancé). Si, cependant, il n'y a plus de diagraphie cible, il est déterminé s'il faut davantage de marqueurs ou si l'utilisateur souhaite davantage de marqueurs (182). Si oui, le programme retourne à l'étape 124 dans le traitement de fonction de guidage (122) (figures 13A et 13C). S'il ne doit pas y avoir davantage de marqueurs, cependant, le programme avance à fin (184) (figure 13A). A partir des procédures cidessus (option automatique ou option semi-automatique), la fonction AMCA 30 peut être obtenue ainsi que son maximum (par exemple, 30A sur la figure 2) pour une courbe de diagraphie de référence et une ou plusieurs courbes de diagraphie cibles et, de ce fait, peut identifier des caractéristiques connexes (sommet de couche) entre les courbes. Selon un mode de réalisation de l'invention, les courbes de diagraphie de référence et cibles avec la caractéristique connexe identifiée (c'est-à-dire que les signaux enregistrés sont identifiés à partir de la même couche supérieure de formation) peuvent être affichées sur le support d'enregistrement de sortie 20 (par exemple, un écran d'ordinateur) pour un utilisateur. Un tel affichage peut être d'apparence similaire à la figure 18 avec les emplacements de marqueurs 94a-g de sommet de formation (événement) (similaires aux marqueurs 80A-F) identifiant la caractéristique connexe parmi les courbes. Les caractéristiques d'un tel affichage peuvent s'avérer utiles pour l'utilisateur pour l'identification de l'emplacement (c'est-à-dire que les emplacements des marqueurs peuvent être convertis en profondeurs absolues approximatives dans les puits) de sources potentielles d'hydrocarbures (par exemple, du pétrole) dans un horizon (par exemple, dans des couches supérieures de formation particulières) et/ou à des emplacements de failles. L'horizon et/ou les failles sont ceux qui ont été sondés dans les puits forés pour collecter les données de courbes de diagraphie de référence et cibles par lesquelles on doit commencer. L'utilisateur est capable d'utiliser les informations affichées pour identifier des emplacements pour forer des puits supplémentaires vers des sommets de couches de formation à des profondeurs absolues approximatives particulières par rapport à la surface de la terre où des hydrocarbures sont susceptibles d'être trouvés pour être récupérés. Selon d'autres modes de réalisation de l'invention, le procédé semi-automatique peut être utilisé en même temps que le procédé automatique pour procurer une flexibilité à l'utilisateur pour des réglages géologiques divers. Par exemple, si les marqueurs les plus superficiels sont sélectionnés en utilisant le procédé semi-automatique, ensuite les marqueurs suivants (en profondeur) peuvent être contraints par une règle de non-croisement et la fenêtre de diagraphie de référence ne s'étendra pas (ne croisera pas) au-delà des marqueurs ci-dessus. Avec référence à la figure 14, un procédé selon un mode de réalisation de l'invention est montré. Ce procédé peut être utilisé pour évaluer la fonction AMCA 30 une fois qu'elle a été calculée (une opération de traitement ultérieure). La fonction AMCA 30 peut être évaluée avec ce procédé entre deux horizons de guidage 50A (par exemple, 52A et 52B de la figure 5) comme suit : (1) si toutes les valeurs de Cmax(X) 30 sont (352) inférieures au paramètre de coupure 50D, l'événement présentant un intérêt n'est pas trouvé (354) ; (2) une section de la courbe de diagraphie cible (par exemple, de la courbe 23 sur la figure 3) correspondant à une section d'une courbe de diagraphie de référence (par exemple, de la courbe 21 sur la figure 3) est supposée (356) manquante et le support d'enregistrement de sortie 20 (par exemple, un écran) est marqué (358) (par exemple, par un point d'interrogation ["?"]) ; (3) d'autre part, si toutes les valeurs de Cmax(7) 30 sont (360) supérieures au paramètre de coupure 50D, plusieurs positions de corrélation T peuvent dépasser (362) le paramètre de coupure 50D ; et (4) la position de corrélation choisie T peut être déterminée (364) en appliquant (déterminées par l'utilisateur) les règles géologiques 50 comme suit (autrement, l'utilisateur reçoit l'instruction de remettre à zéro le paramètre de coupure (361)) ; (5) si une séquence superposée 54A doit être appliquée (366), la valeur de Cmax(X) 30 dépassant le paramètre de coupure 50D et la plus proche de la position parallèle au-dessous de l'horizon de surface de limite supérieure 52A est sélectionnée (368) comme le maximum de corrélation (par exemple, 30A de la figure 2) ; (6) si une séquence tronquée 54B doit être appliquée (370), la valeur de Cmax(l) 30 dépassant le paramètre de coupure 50D et la plus proche de la position parallèle audessus de l'horizon de surface de limite inférieure 52B est sélectionnée (372) comme le maximum de corrélation ; (7) si une séquence conforme 54C doit être appliquée (374), la valeur de Cmax(r) 30 dépassant le paramètre de coupure 50D et la plus proche de la distance proportionnelle par rapport à l'un ou l'autre horizon de surface de limite 52A ou 52B est sélectionnée (376) comme le maximum de corrélation (8) si une fonction de séquence non structurée 54D doit être appliquée (378), la valeur de Cmax(T) 30 avec la valeur la plus haute dépassant le paramètre de coupure 50D est sélectionnée (380) comme le maximum de corrélation ; et (9) si une fonction de guidage de l'utilisateur 50C doit être appliquée (382), la valeur de Cmax(T) 30 avec la valeur la plus haute dépassant le paramètre de coupure 50D la plus proche d'un emplacement sélectionné par l'utilisateur est sélectionnée (384) comme le maximum de la fonction AMCA 30. Sans se soucier de laquelle parmi les règles de séquence 50B ci-dessus ou la fonction de guidage par l'uti l'utilisateur 50C n'est pas appliquée, alors, une erreur s'est produite (386) qui peut être indiquée à l'utilisateur. Il existe des variables qui affectent la valeur maximale de la fonction AMCA 30. Ces variables peuvent comprendre la longueur w de la fenêtre de référence, la longueur Wy de la fenêtre cible et des critères pour la correspondance de forme et de dimension en plus de la covariance. L'utilisateur peut contrôler la sensibilité de la fonction AMCA 30 à la covariance, à la forme et à la dimension en contrôlant les paramètres a, ss, 7. La longueur w contrôle quelle proportion de l'événement sur la courbe de référence est utilisée pour la comparaison avec la cible. Si la longueur w pour la fenêtre 26 est trop petite comparée à la largeur de l'événement sur la courbe de diagraphie de référence, l'événement présentant un intérêt ne sera pas trouvé (c'est-à-dire qu'une bonne corrélation ne peut pas être trouvée et qu'une erreur de sélection se produira, par exemple, aucun marqueur ou un marqueur incorrect sera identifié sur la courbe cible), comme montré schématiquement sur la figure 15. Un marqueur incorrect (par exemple, le marqueur 81 sur la figure 15) peut être reconnu par un spécialiste des sciences de la terre expérimenté dans de nombreux cas. De même, si la longueur w de la fenêtre 26 est trop grande, également, l'événement présentant un intéret ne sera pas trouvé (par exemple, il y aura une erreur de sélection), comme montré schématiquement sur la figure 16, parce que la valeur du maximum de corrélation devient moins distincte par rapport aux valeurs du fond de la fonction AMCA 30. De plus, si l'événement de référence comporte essentiellement des composantes basse fréquence, le calcul de la fonction de covariance 30 peut s'exécuter mieux avec des valeurs plus grandes de w ; pour des composantes essentiellement haute fréquence, de plus petites valeurs de w peuvent être plus appropriées. La figure 17 montre un exemple d'une fenêtre 26 dont la longueur w a été sélectionnée avec soin pour l'événement présentant un intérêt contenu à l'intérieur de celle-ci. La figure 17 montre la courbe de diagraphie de référence 21, les courbes de diagraphie cibles 23A-F (chacune similaire à la courbe cible 23 sur la figure 3) provenant de puits différents et les marqueurs 82 (sélectionnés par l'utilisateur) et 80A-F (trouvés selon un mode de réalisation de l'invention, par exemple, le procédé semi-automatique). Sur la figure 17, aucun des marqueurs 80A-F n'est mal sélectionné, à la différence des figures 15 et 16. Selon un mode de réalisation de l'invention, tel que montré sur la figure 18, un événement A a été sélectionné auparavant et les marqueurs ont été sauvegardés (soit à partir du procédé automatique, soit à partir du procédé semi-automatique décrit ci-dessus) Il est maintenant disponible pour l'utilisateur comme horizon de guidage 94 (similaire à l'horizon de guidage 50A). Le début 92B des fenêtres (similaires à la fenêtre 28) pour sélectionner l'événement B est ensuite commandé en utilisant la profondeur de l'horizon de guidage A 94 et la règle géologique sélectionnée. Pour cet exemple, la règle de superposition est appliquée de sorte que la fenêtre (Wy) pour chaque courbe cible soit parallèle à l'horizon de guidage A 94. Cela peut être accompli en réglant le centre 95 (non montré sur toutes les courbes de diagraphie cibles) de la fenêtre (Wy) à la même profondeur au-dessous de l'horizon de guidage A 94 sur les diagraphies cibles que le centre 93 de la fenêtre w sur la diagraphie de référence. Une pondération plus élevée peut être appliquée au terme de dimension (le terme multiplié par y) de la fonction AMCA 30 en augmentant la valeur y. Cela produira un plus grand maximum de corrélation (fonction AMCA 30) (non montré) lorsque les amplitudes des échantillons entre la diagraphie de référence et la diagraphie cible sont les plus similaires (comme avant, cc + P + 7 = 1,0). Les nouvelles sélections pour l'événement B sont montrées sur la figure 19. D'autre part, sur la figure 20, les résultats provenant d'une position de début constante (non montrée) des fenetres de recherche (Wy) sont montrés. Ici, l'horizon A 94' (similaire à l'horizon de guidage 50A) est plat. A partir des figures 19 et 20, il est évident que la position de début Wy a un effet sur les emplacements des maximums de corrélation (fonction AMCA 30) qui correspondent aux marqueurs 96 (similaires aux marqueurs 8OA-F sur la figure 17 et 94b-g sur la figure 18) sur les courbes cibles. Certains marqueurs n'apparaissent même pas (c'est-à-dire qu'il peut y avoir des erreurs de sélection) dans les données de diagraphie cibles dans les courbes cibles 23G et 23H (toutes deux similaires à la courbe cible 23) sur la figure 20. Les maximums de la fonction AMCA 30 peuvent également être affectés par la correspondance de forme et de dimension qui peut être utilisée pour pondérer la similarité de forme et de dimension d'événements sur les courbes de diagraphie de référence et cibles. Si des critères de correspondance de forme et de dimension doivent être utilisés, ils peuvent être appliqués en modifiant les valeurs des paramètres a, ss, 7 du calcul de Cmax(X) 30. La pondération qui affecte la forme est le terme ss et une pondération élevée produira ici un plus grand maximum de corrélation (fonction AMCA 30) lorsque les diagraphies de référence et cibles ont une dispersion de valeurs (écarts types) similaire. Outre une courbe de diagraphie de référence fixe (par exemple, la courbe de diagraphie X 21 de la figure 3), selon le mode de réalisation semi-automatique de l'invention, cependant, un procédé a été développé dans lequel la courbe de diagraphie de référence peut varier ("rouler") du puits de sondage sélectionné par l'utilisateur à la courbe de diagraphie de trou de sondage qui précède immédiatement la courbe de diagraphie de trou de sondage cible courante en cours d'évaluation. Cela peut être avantageux dans des environnements géologiques qui varient systématiquement dans une certaine direction (par exemple, une formation qui s'épaissit continuellement dans la direction du nord). . L'utilisateur peut sélectionner plusieurs puits cibles dans un ordre tel que les sommets des formations (enregistrés comme des événements sur les courbes de diagraphie) présentant un intérêt seront vues comme changeant graduellement de la diagraphie de référence initiale à la diagraphie cible finale en passant par les diagraphies cibles intermédiaires, où une comparaison directe des diagraphies de référence d'origine et cible finale serait assez différente. Avec référence aux figures 10 et 12, la corrélation de retour et la fermeture de boucle qui ont été décrites sont des exemples du procédé de roulement. Le tableau 1 montre des exemples de résultats obtenus pour la covariance maximale (fonction AMCA 30) pour les références fixes et roulantes, ainsi que des profondeurs sélectionnées. Tableau 1 Référence Référence Référence Référence fixe fixe "roulante" "roulante" Trou de Covariance Profondeur Covariance Profondeur sondage maximale sélectionnée maximale sélectionnée 2 0,994385 5855 0,994385 5855 3 0,998366 5853 0,993058 5855 4 0,995838 5835 0,996672 5802 5 0,996072 5854 0,999435 5747 6 0,993388 5820 0,999432 5725 7 0,995725 5819 0,998902 5760 A partir du Tableau 1, il apparaît que le procédé de référence "roulante" a produit de meilleures corrélations que l'approche à référence fixe pour cet exemple. Le procédé de guidage par l'utilisateur (semiautomatique) décrit ci-dessus peut utiliser soit l'option de diagraphie de référence fixe, soit l'option de diagraphie à référence "roulante". D'autres modes de réalisation sont compris dans l'étendue des revendications qui suivent. REVENDICATIONS 1. Procédé d'interprétation de courbes de diagraphie de référence et cibles enregistrées à partir de différents puits pour aider à la détermination de la profondeur approximative de caractéristiques correspondantes aux emplacements de puits, comportant le prétraitement des courbes de diagraphie de référence et cibles selon des paramètres entrés par l'utilisateur ; l'application d'une fonction de guidage aux courbes qui limite l'interprétation à des parties des courbes ; le calcul d'une indication de correspondance entre les parties ; et l'évaluation de 11 indication de correspondance pour déterminer si des événements sur les parties sont liés. 2. Procédé selon la revendication 1 comportant de plus le calcul d'une fenêtre cible pour les données cibles - 3. Procédé selon la revendication 1, comportant de plus l'application d'une règle de séquence conforme aux données de référence, dans laquelle un étirement proportionnel ou une compression proportionnelle est appliqué. 4. Procédé selon la revendication 1, comportant de plus l'inversion des données de référence. 5. Procédé selon la revendication 1 comportant de plus l'application de pondérations aux parties des données de référence et cibles. 6. Procédé selon la revendication 1 comportant de plus l'affichage des données de référence et cibles sur un écran. 7. Appareil pour l'interprétation de courbes de diagraphie de référence et cibles enregistrées à partir de différents puits pour aider à déterminer la profondeur approximative de caractéristiques correspondantes, aux emplacements de puits, comportant une mémoire (14) qui stocke des instructions pour interpréter les données de diagraphie; et un processeur (12) couplé à la mémoire (14) qui est chargé par les instructions de prétraiter les données de diagraphie de référence et cibles selon des paramètres entrés par l'utilisateur, appliquer une fonction de guidage aux données de diagraphie prétraitées qui limite l'interprétation à des parties de données, calculer une indication de correspondance entre les parties, et évaluer l'indication de correspondance pour déterminer si des événements sur les parties sont liés. 8. Système pour l'analyse de données de courbes de diagraphie enregistrées à partir de puits répartis géographiquement pour déterminer l'emplacement de sommets de formation correspondant dans les puits, comportant des dispositifs d'entrée/sortie (E/S) qui reçoivent des données de diagraphie enregistrées à partir de puits ; et un processeur (12) couplé aux dispositifs d'entrée/sortie (E/S) qui est chargé de prétraiter les données de diagraphie selon des paramètres entrés par l'utilisateur, déterminer si une fonction de guidage doit être appliquée aux données, effectuer une détermination quant à savoir si les événements dans les données sont corrélés, en se basant en partie sur les paramètres entrés par l'utilisateur, et si les événements sont liés, convertir un emplacement associé à la corrélation en un emplacement approximatif des sommets de formation dans les puits. 9. Procédé d'évaluation de la correspondance entre des événements sur une courbe de diagraphie de référence et une courbe de diagraphie cible enregistrées dans des puits de référence et cibles pour identifier un emplacement approximatif de couches correspondantes dans les puits, comportant déterminer si un indicateur de correspondance entre les événements sur les courbes de diagraphie de référence et cibles à une valeur supérieure à un paramètre de coupure appliquer un guide géologique prescrit pour déterminer l'emplacement d'un extremum de l'indicateur de correspondance sur la courbe de diagraphie cible ; et convertir l'emplacement de l'extremum en un emplacement absolu approximatif dans le puits cible. 10. Programme informatique, résidant sur un support susceptible d'être lu par ordinateur, comportant des instructions pour amener un ordinateur à : prétraiter des données de diagraphie de référence et cibles selon des paramètres entrés par l'utilisateur; appliquer une fonction de guidage aux données de diagraphie qui dirige le traitement des données de diagraphie; calculer une indication de correspondance entre les parties des données de diagraphie; et évaluer l'indication de correspondance pour déterminer si les parties sont liées.