SUB-SURFACE WELL LOCATION DETERMINATION

A method of estimating the current position of a drill bit or other noise source within a subsurface formation of the Earth. The method comprises collecting seismic data from a plurality of seismic sensors spread across a region of the surface of the Earth above the drill bit or other noise source and / or located in one or more nearby wells, pre-processing the seismic data to enhance a contribution of drill bit or noise source generated noise, and defining a set of points on a grid in 3-dimensional space that includes an expected position of the drill bit or noise source. Travel times for seismic waves from each said point to each seismic sensor location are located and, for each said point, the pre-processed seismic data, sensor location data, and computed travel times are computed to compute a semblance stack of travel-time corrected seismic data in a time window. The grid location of the maximum semblance is determined and a 3-dimensional polynomial fitted around this grid location. The location of a maximum of the 3-dimensional polynomial is located and used as an estimate of the current position of the drill bit or other noise source. La présente invention concerne un procédé d'estimation de la position actuelle d'un trépan ou d'une autre source de bruit dans une formation souterraine de la Terre. Le procédé consiste à collecter des données sismiques provenant d'une pluralité de capteurs sismiques déployés sur une région de la surface de la Terre au-dessus du trépan ou d'une autre source de bruit et/ou situés dans un ou plusieurs puits proches, à prétraiter les données sismiques pour améliorer une contribution du bruit généré par un trépan ou une source de bruit, et à définir un ensemble de points sur une grille dans un espace tridimensionnel qui comprend une position attendue du trépan ou de la source de bruit. Des temps de déplacement pour des ondes sismiques depuis chaque dit point jusqu'à chaque emplacement de capteur sismique sont localisés et, pour chaque dit point, les données sismiques prétraitées, les données d'emplacement de capteur et les temps de déplacement calculés sont calculés pour calculer une pile de semblance de données sismiques corrigées par le temps de déplacement dans une fenêtre temporelle. L'emplacement de grille de la semblance maximale est déterminé et un polynôme tridimensionnel est ajusté autour de cet emplacement de grille. L'emplacement d'un maximum du polynôme tridimensionnel est localisé et utilisé comme estimation de