A method for estimating uncertainties in determining hypocenters of seismic events occurring in subsurface formations according to one aspect includes determining estimates of event locations by choosing local peaks in summed amplitude of seismic energy detected by an array of sensors disposed above an area of the subsurface to be evaluated. For each peak, the following is performed: recomputing the summed amplitude response for a selected set of points of comprising small perturbations in time and space from the estimated event locations; computing second derivatives of log likelihood function from the stacked responses at the estimated location and the perturbed locations; assembling the second derivatives into a Fisher information matrix; computing an inverse of the Fisher information matrix; determining variances of estimated parameters from the elements from the diagonal of the inverted matrix; and computing standard deviations of the estimated parameters by calculating a square root of the variances.

Correlation matrices may be used to simultaneously correlate multiple wells. A correlation matrix may be generated for individual pairs of multiple wells. The values of elements of the correlation matrices may be determined based on matching between segments of the multiple wells and segments of one or more computational stratigraphic models. An N-dimensional space including an axis for individual wells may be generated. Directed walk may be performed within the N-dimensional space to generate paths representing scenarios of correlations for segments of the multiple wells.

A scintillation crystal can include a cesium halide that is co-doped with thallium and another element. In an embodiment, the scintillation crystal can include CsX:Tl, Me, where X represents a halogen, and Me represents a dopant selected from the group consisting of chromium (Cr), zirconium (Zr), cobalt (Co), manganese (Mn), cadmium (Cd), dysprosium (Dy), thulium (Tm), tantalum (Ta), and erbium (Er), the dopant concentration of the element selected from the group consisting of chromium (Cr), zirconium (Zr), cobalt (Co), manganese (Mn), cadmium (Cd), dysprosium (Dy), thulium (Tm), tantalum (Ta), and erbium (Er) in the scintillation crystal is in a range of 1×10.sup.−7 mol % to 0.5 mol %. In a particular embodiment, the scintillation crystal may have a cesium iodide host material, a first dopant including a thallium cation, and a second dopant including a cation.

A distributed acoustic sensing cable package having a polymer composite extruded over an optical waveguide to encase the waveguide and to form a crystalline matrix layer acoustically coupled to the waveguide. The crystalline matrix includes reinforcement fibers to further enhance transmission of a cable strain to the optical waveguide. During manufacture of the cable, the polymer composite may be extruded over the optical waveguide and subsequently subjected to heat treatment to increase the crystallinity of the polymer composite and increase the elastic modulus. Both axial and radial strain fields are effectively interact with cased fiber waveguide for producing measurable phase shift signal for distributed acoustic noise detection.

A distributed acoustic sensing cable package having a polymer composite extruded over an optical waveguide to encase the waveguide and to form a crystalline matrix layer acoustically coupled to the waveguide. The crystalline matrix includes reinforcement fibers to further enhance transmission of a cable strain to the optical waveguide. During manufacture of the cable, the polymer composite may be extruded over the optical waveguide and subsequently subjected to heat treatment to increase the crystallinity of the polymer composite and increase the elastic modulus. Both axial and radial strain fields are effectively interact with cased fiber waveguide for producing measurable phase shift signal for distributed acoustic noise detection.

A method of monitoring hydraulic fracturing using DAS sensors in a treatment well and/or observation well is described. The raw data is transformed using a low pass filter (≦0.05 Hz) and down-sampled to show the signal as the stimulation progresses. The resulting data can be used to optimize the hydraulic fracturing or improve reservoir models for other reservoirs.

A self-adapting digital twin of a wellbore environment can be created. The self-adapting digital twin incorporates a standalone Ensemble Kalman Filter (EnKF) module with a constant parameter digital twin developed for a fixed environment. The standalone EnKF module receives streaming measurement data from multiple sensors and prediction data from the digital twin and executes the standalone EnKF module using the streaming measurement data and the prediction data from the digital twin. The results of executing the standalone EnKF module are input parameter corrections for the digital twin that are communicated to the digital twin. Output predictions of the digital twin are used to modify operational parameters of an oil or gas recovery process.

An exploration method starts from cuttings associated with sampling intervals and well data for a well in a subsurface formation. The cuttings are prepared and analyzed to extract textural and chemical/mineralogical data for plural fragments in each sample that is made of the cuttings in one sampling interval. The method then includes matching lithotypes of rock defined according to the textural and chemical/mineralogical data for each fragment with segments of the well data in the corresponding sampling interval to obtain correspondences between the lithotypes and depth ranges. The correspondences between the lithotypes and the depth ranges may be used as constraints for seismic data inversion.

An exploration method starts from cuttings associated with sampling intervals and well data for a well in a subsurface formation. The cuttings are prepared and analyzed to extract textural and chemical/mineralogical data for plural fragments in each sample that is made of the cuttings in one sampling interval. The method then includes matching lithotypes of rock defined according to the textural and chemical/mineralogical data for each fragment with segments of the well data in the corresponding sampling interval to obtain correspondences between the lithotypes and depth ranges. The correspondences between the lithotypes and the depth ranges may be used as constraints for seismic data inversion.

Disclosed are tools, systems, and methods for detecting one or more underground acoustic sources and localizing them in depth and radial distance from a wellbore, for example, for the purpose of finding underground fluid flows, such as may result from leaks in the well barriers. In various embodiments, acoustic-source detection and localization are accomplished with an array of at least three acoustic sensors disposed in the wellbore, in conjunction with array signal processing.