A fracture model for a hydraulic fracture in a wellbore can be updated and calibrated. Information about a microseismic event can be received from a sensor that is monitoring a subterranean formation. The information can be received subsequent to a fracking fluid being introduced into the formation. An observed geometry of a hydraulic fracture can be determined based on the information and a predicted geometry of the fracture can be determined based on properties of the fracking fluid and a fracture model. The fracture model can be updated using the information about the microseismic event where it is determined that an uncertainty value of the observed geometry does not exceed a pre-set maximum. The uncertainty value can be based on the predicted geometry of the hydraulic fracture.

A symmetry attribute is described that may be used for determining seismic stratigraphic features in a formation. In one example, seismic input data from a formation is processed to determine an attribute by selecting a center trace, assigning a first cluster of the traces to a left image and a second cluster of the traces to a right image, and determining symmetry about the center trace between the left and the right images.

The present invention is a method to determine an azimuthal orientation of a borehole seismometer sensor performed by a computing device using a control server having a database and an arithmetic function, the computing device performing a method to determine the azimuthal orientation of a borehole seismometer sensor using long-period surface waves in microseisms, including step S100 in which a data collection unit 100 collects continuous waveform data recorded by a borehole seismometer and a reference seismometer; step S200 in which a frequency band setting unit 200 sets a frequency band to be analyzed in the collected continuous waveform data; step S300 in which a filtering unit 300 performs bandpass filtering on the frequency band to be analyzed; step S400 in which a waveform dividing unit 400 divides seismic waveform into waveform segments with preset time units; step S500 in which a phase shift unit 500 shifts the phase of the divided vertical component waveforms by 90°; step S600 in which a waveform calculation unit 600 combines the divided N′ and E′ component seismic waveforms to calculate horizontal components for rotation angles waveform between 0 and 360° from the N′ orientation; step S700 in which a correlation calculation unit 700 calculates a correlation coefficient between the horizontal and vertical component waveforms; step S800 in which a Rayleigh wave orientation determination unit 800 repeats steps S500 to S700 for each divided time domain; step S900 in which an orientation comparison unit 900 performs steps S400 to S800, respectively, with respect to the borehole seismometer data for which the sensor orientation is to be determined and the reference seismometer data for which the sensor orientation is already known; and step S1000 in which a result calculation unit 1000 averages 0 determined for each time period to calculate a final result.

Uncertainty in microseismic monitoring sensor data can be reduced. A computing device can receive information about at least one sensor that is monitoring a subterranean formation, including a location, after a fracturing fluid is introduced into the formation. The computing device can also receive information about a microseismic event and determine a seismic ray bath between a location of the event and the at least one sensor, and an uncertainty value of the location based on information about the formation and the information about the event. The computing device can determine a total uncertainty value associated with the locations of a plurality of microseismic events, including the microseismic event. The computing device can determine a solution to an objective function based on the total uncertainty value and a number of sensors. The computing device can determine a new location of the at least one sensor based on the solution.

A rock burst hazard prediction method based on vibration wave energy attenuation characteristics of a mine earthquake cluster is provided. The rock burst hazard prediction method comprehensively considers the static load and dynamic load effects of the vibration waves of the mine earthquake cluster based on vibration wave energy attenuation characteristics of the mine earthquake cluster. The static load strength index and the dynamic load strength index involved in the method have clear physical meanings. A comprehensive prediction index calculation model proposed based on the dynamic and static load superposition principle of rock burst occurrence is clear, and the method has a firm theoretical support as well as strong universality and operability. Meanwhile, the updating and adjustment of weights are rapid and the objective judgment and prediction of the final comprehensive prediction results are efficient, and the high-energy mine earthquake and impact behavior area can be effectively predicted.

A system (100) for monitoring a subterranean structure comprises an array (10) with n acoustic sensors capable of detecting P-waves and/or S-waves from the subterranean structure and a central controller (120) for receiving a signal (X) from the sensors. The system further comprises a lookup table (20) comprising a pre-computed travel time curve (24) expressed as relative arrival times of a signal from a location (L.sub.m) to each of the sensors (1−n); a comparison unit for comparing the received signal (X) with the pre-computed travel time curve (24), and means for raising an alarm if the received signal (X) matches the precomputed travel time curve (24). Preferably, the alarm is raised if a computed semblance value (26, 27) exceeds a predefined threshold. The system may monitor several locations (L.sub.m) in parallel using a fraction of the computer resources and time required by prior art techniques.

A method for improved data-driven estimation of a stimulated reservoir volume may generate an optimized surface that encloses a set of data points including microseismic event data corresponding to a treatment of a subterranean formation. A Delaunay triangulation may be performed on the set of data points to generate a set of polytopes. A Voronoi polygon may be generated about each data point and used to obtain a local density measure that is locally and adaptively determined for each data point. Based on the local density measure, polytopes in the set of polytopes may be discriminated for inclusion in the optimized surface.

A data acquisition program, which includes core, image log, microseismic, DAS, DTS, and pressure data, is described. This program can be used in conjunction with a variety of techniques to accurately monitor and conduct well stimulation.

An embodiment of a method of stimulating an earth formation includes: disposing a stimulation device at a borehole in an earth formation, the earth formation having been stimulated by an initial stimulation operation; subsequent to the stimulation operation, performing a probe operation configured to cause movement of existing fractures in the formation; and measuring microseismic events occurring in the formation by one or more seismic receivers. The method further includes: identifying one or more target zones in the formation based on the measuring, the one or more target zones exhibiting a reduced micro seismicity relative to another zone in the formation; and designing a re-stimulation operation configured to stimulate the one or more target zones to increase hydrocarbon production from the formation.

A buried hydrophone may be configured to identify real-time changes in pipeline acoustics and vibrations during a flood event. Monitoring vibrations may allow predicting the likelihood of pipeline exposure and suspension.