Abstract: A method for estimating a fluid pressure required to stimulate a subsurface formation includes using seismic signals detected by a plurality of seismic sensors disposed proximate the subsurface formation. Spatial positions and times of origin (“hypocenters”) of each of a plurality of microseismic events induced by pumping fluid into the subsurface formation are estimated. Magnitudes and directions of principal stresses are estimated from the hypocenters and from amplitude and phase of the detected seismic signals for each of the microseismic events. Shear and normal stresses of induced fractures are from the estimated principal stresses. A fluid pressure required to cause formation failure on each fracture is estimated using the estimated shear and normal stresses.

Abstract: A method for determining hypocenters of microseismic events includes entering as input to a computer seismic signals recorded by a plurality of seismic sensors disposed proximate a volume of subsurface to be evaluated. For each point in space in the volume, and for a plurality of preselected origin times, a seismic energy arrival time at each seismic sensor is determined. Event amplitudes for each arrival time are determined. A synthetic event amplitude is calculated for each arrival time. A semblance between the determined event amplitudes and the synthetic event amplitudes is determined. Existence of an actual microseismic is determined event when the semblance exceeds a selected threshold.

Abstract: A method for determining a stimulated rock volume includes determining a position of a plurality of seismic events from seismic signals recorded in response to pumping fracturing fluid into a formation penetrated by a wellbore. The signals generated by recording output of a plurality of seismic receivers disposed proximate a volume of the Earth's subsurface to be evaluated. A source mechanism of each seismic event is determined and is used to determine a fracture volume and orientation of a fracture associated with each seismic event. A volume of each fracture, beginning with fractures closest to a wellbore in which the fracturing fluid was pumped is subtracted from a total volume of proppant pumped with the fracture fluid until all proppant volume is associated with fractures. A stimulated rock volume is determined from the total volume of fractures associated with the volume of proppant pumped.

Abstract: A method for imaging microseismic events includes determining a hypocenter of microseismic events generated by at least one stage of a hydraulic fracturing procedure from recorded signals detected by seismic sensors disposed above a wellbore in the subsurface. Spatial position of the microseismic events occurring sequentially in the fracturing procedure is determined with reference to a center of fracturing procedure. Each microseismic event is assigned to one of a plurality of selected size bins defined positionally with reference to the center of the fracturing procedure. A property of each microseismic event assigned to each bin is aggregated and an image of the aggregated property is generated with respect to position referenced to the center of the fracturing procedure.

Abstract: A method for determining spatial distribution of proppant incudes using signals detected by seismic sensors disposed proximate a formation treated by pumping fracturing fluid containing the proppant. Origin time and spatial position of seismic events induced by pumping the fracturing fluid are determined. Volume and orientation of at least one fracture in the subsurface formation associated with each induced seismic event are determined. Spatial distribution of a volume of the pumped fracturing fluid is determined using the volume and orientation of each fracture. A length of ellipsoidal axes is selected using a surface defined by a selected fractional amount of the total volume of frac fluid pumped into the formation. Spatial distribution of the proppant is determined using proppant mass, specific gravity and expected proppant porosity in the fractures, and spatially distributing a volume of the fractures within an ellipsoid defined by the ellipsoidal axes.

Abstract: 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.

Abstract: A method for determining spatial distribution of proppant incudes using signals detected by seismic sensors disposed proximate a formation treated by pumping fracturing fluid containing the proppant. Origin time and spatial position of seismic events induced by pumping the fracturing fluid are determined. Volume and orientation of at least one fracture in the subsurface formation associated with each induced seismic event are determined. Spatial distribution of a volume of the pumped fracturing fluid is determined using the volume and orientation of each fracture. A length of ellipsoidal axes is selected using a surface defined by a selected fractional amount of the total volume of frac fluid pumped into the formation. Spatial distribution of the proppant is determined using proppant mass, specific gravity and expected proppant porosity in the fractures, and spatially distributing a volume of the fractures within an ellipsoid defined by the ellipsoidal axes.

Abstract: A method for determining hypocenters of microseismic events includes entering as input to a computer seismic signals recorded by a plurality of seismic sensors disposed proximate a volume of subsurface to be evaluated. For each point in space in the volume, and for a plurality of preselected origin times, a seismic energy arrival time at each seismic sensor is determined. Event amplitudes for each arrival time are determined. A synthetic event amplitude is calculated for each arrival time. A semblance between the determined event amplitudes and the synthetic event amplitudes is determined. Existence of an actual microseismic is determined event when the semblance exceeds a selected threshold.

Abstract: A method for determining positions and origin times of seismic events occurring in the Earth's subsurface includes accepting as input to the method signals recorded from a plurality of seismic sensors deployed above a subsurface volume of interest. The recorded signals are a representation of seismic amplitude with respect to time. Origin time and location of each of a plurality of subsurface seismic events are determined from the recorded signals. The origin times and locations of each event are inverted to obtain Thomsen's parameters in formations in the volume of interest. Depths of each of the events are determined by individually searching the depth of each event, the inversion with each incorporated new depth including updating the Thomsen parameters and setting as a limit a minimum value of RMS error.

Abstract: A method for determining a stimulated rock volume includes determining a position of a plurality of seismic events from seismic signals recorded in response to pumping fracturing fluid into a formation penetrated by a wellbore. The signals generated by recording output of a plurality of seismic receivers disposed proximate a volume of the Earth's subsurface to be evaluated. A source mechanism of each seismic event is determined and is used to determine a fracture volume and orientation of a fracture associated with each seismic event. A volume of each fracture, beginning with fractures closest to a wellbore in which the fracturing fluid was pumped is subtracted from a total volume of proppant pumped with the fracture fluid until all proppant volume is associated with fractures. A stimulated rock volume is determined from the total volume of fractures associated with the volume of proppant pumped.

Abstract: A method for determining a volume of a fracture network includes detecting seismic signals deployed over an area of the subsurface during pumping of fracturing fluid into at least one wellbore drilled through the area. A hypocenter of each fracture induced by the pumping is determined using the seismic signals. A facture network and associated fracture volume is determined using the determined hypocenters and seismic moments determined from the detected seismic signals. A maximum value of a scaling factor is determined based on a subset of the hypocenters having a highest cumulative seismic moments. The scaling factor is determined by relating a pumped volume of the fracturing fluid with respect to the determined fracture volume. Dimensions of each fracture are scaled using the maximum value of the scaling factor. The fracture volumes are recalculated using the scaled dimensions.

Abstract: A method for estimating moment magnitude of a seismic event occurring in subsurface formations includes measuring seismic signals at each of a plurality of seismic sensors disposed in a selected pattern proximate a subsurface area in which the seismic event occurs. Amplitude events corresponding to the seismic event from the signals detected by each receiver are time aligned. Corrections are applied to the aligned events for density, for the formation velocity, for the radiation pattern, for propagation effects and instrument response. The corrected events are summed. Seismic moment is determined from the summed, corrected events. A moment magnitude is estimated from the seismic moment.

Abstract: A method for determining fracture plane orientation from seismic signals detected above a subsurface formation of interest includes detecting seismic signals using an array of seismic sensors deployed above the subsurface formation during pumping of a hydraulic fracture treatment of the subsurface formation. A time of origin and a spatial position of origin (hypocenter) of microseismic events resulting from the hydraulic fracture treatment are determined. Time consecutively occurring ones of the hypocenters falling within a selected temporal sampling window are selected. A best fit line through the selected hypocenters using a preselected linear regression coefficient is determined. The selecting hypocenters and determining best fit lines is repeated for a selected number of windows.

Abstract: A method for mapping a fracture network that includes determining a source of at least one seismic event from features in recorded seismic signals exceeding a selected amplitude (“visible seismic event”). The signals are generated by a plurality of seismic receivers disposed proximate a volume of subsurface to be evaluated. The signals are electrical or optical and represent seismic amplitude. A source mechanism of the at least one visible seismic event is determined. A fracture size and orientation are determined from the source mechanism. Seismic events are determined from the signals from features less than the selected amplitude (“invisible seismic events”) using a stacking procedure. A source mechanism for the invisible seismic events is determined by matched filtering. At least one fracture is defined from the invisible seismic events. A fracture network model is generated by combining the fracture determined from the visible seismic event with the fracture determined from the invisible seismic events.

Abstract: 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.

Abstract: A method for determining a volume of a fracture network includes detecting seismic signals deployed over an area of the subsurface during pumping of fracturing fluid into at least one wellbore drilled through the area. A hypocenter of each fracture induced by the pumping is determined using the seismic signals. A facture network and associated fracture volume is determined using the determined hypocenters and seismic moments determined from the detected seismic signals. A maximum value of a scaling factor is determined based on a subset of the hypocenters having a highest cumulative seismic moments. The scaling factor is determined by relating a pumped volume of the fracturing fluid with respect to the determined fracture volume. Dimensions of each fracture are scaled using the maximum value of the scaling factor. The fracture volumes are recalculated using the scaled dimensions.

Abstract: A method for determining fracture plane orientation from seismic signals detected above a subsurface formation of interest includes detecting seismic signals using an array of seismic sensors deployed above the subsurface formation during pumping of a hydraulic fracture treatment of the subsurface formation. A time of origin and a spatial position of origin (hypocenter) of microseismic events resulting from the hydraulic fracture treatment are determined. Time consecutively occurring ones of the hypocenters falling within a selected temporal sampling window are selected. A best fit line through the selected hypocenters using a preselected linear regression coefficient is determined. The selecting hypocenters and determining best fit lines is repeated for a selected number of windows.

Abstract: A method for passive seismic surveying includes deploying seismic sensors in a plurality of spatially distributed wellbores disposed above a volume of subsurface formations to be evaluated. The sensors in each wellbore form a line of sensors. Each sensor generate optical or electrical signals in response to seismic amplitude. The seismic signals from each sensor are recorded for a selected period of time. The response of the seismic sensor recordings is beam steered to at least one of a selected point and a selected volume in the subsurface. At least one microseismic event is identified in the beam steered response.

Abstract: A method for evaluating fractures induced in subsurface formations includes determining hypocenters and origin times of induced seismic events from seismic signals detected by seismic sensors deployed above a volume of the subsurface formations to be evaluated and recording the detected signals for a selected time. A modified Gutenberg-Richter exponent is determined for each induced seismic event having at least a selected number of other induced seismic events within a selected radial distance thereof. A spatial distribution of the modified Gutenberg-Richter exponents is used to determine at least one property of the induced fractures.

Abstract: A method of imaging the Earth's subsurface using passive seismic interferometry tomography includes detecting seismic signals from within the Earth's subsurface over a time period using an array of seismic sensors, the seismic signals being generated by seismic events within the Earth's subsurface. The method further includes adaptively velocity filtering the detected signals. The method further includes cross-correlating the velocity filtered seismic signals to obtain a reflectivity series at a position of each of the seismic sensors.