Methods for dolomite mapping using multiscale fracture characterization include using a computer system to receive seismic data for a geographical area. The computer system identifies one or more macroscale fractures located within the geographical area based on a three-dimensional (3D) visualization of the seismic data. The computer system identifies one or more mesoscale fractures located within the geographical area based on a curvature map generated from the seismic data. The computer system identifies one or more microscale fractures located within the geographical area based on an amount of chaotic seismic reflections indicated by the seismic data. The computer system identifies a dolomite distribution of the geographical area based on the one or more macroscale fractures, the one or more mesoscale fractures, and the one or more microscale fractures. A display device of the computer system generates a graphical representation of the dolomite distribution.

The present disclosure describes methods and systems, including computer-implemented methods, computer program products, and computer systems, for generating a velocity model for a subsurface structure. One computer-implemented method for determining velocity model for a subsurface structure includes generating, by at least one hardware processor, a first velocity model for the subsurface structure by performing a refraction traveltime tomography procedure based on an initial velocity model; and generating, by the at least one hardware processor, a first refined velocity model based on the first velocity model and a structure skeleton model, wherein the structure skeleton model is determined based on reflection seismic data of the subsurface structure.

Seismic data processing may include computing the travel time shift between two seismic signals or the depth shift between two seismic images. In Full Waveform Inversion (FWI), the travel time difference between an observed trace and a simulated trace may be computed such that the two traces match after the travel time shift is applied to the observed trace. The travel time shift may be computed based on a constrained optimization that maximizes the windowed cross-correlation between the two seismic traces by constraining the time derivative of the travel time shift to be less than a constant while maximizing the windowed cross-correlation. Further, the travel time shift may be computed during the model line search in FWI by computing a plurality of travel time shifts where a first travel time shift is dependent on the observed trace and a second travel time shift is independent of the observed trace.

A DAS VSP technique is used to determine the induced fracture height and fracture density of an induced fracture region. The DAS VSP technique obtains pre-hydraulic fracturing DAS VSP survey time-lapse data to establish a baseline reference for the direct acoustic wave travel time. The DAS VSP technique obtains one or more time-lapse data corresponding to the subsequent monitor surveys conducted after each hydraulic fracturing stage along the well. Forward modeling is used to determine a theoretical acoustic wave travel time difference. The forward modeling uses seismic anisotropy to describe the behavior of seismic waves traveling through the induced fracture regions. An inversion scheme is then used to invert for the induced fracture height and the fracture density using the forward modeling. The two extracted induced fracture characteristics may then be used to determine optimal hydraulic fracturing parameters.

Methods and systems including computer programs encoded on a computer storage medium, for utilizing equivalent linear velocity for first arrival picking of seismic refraction. In one aspect, a method includes receiving data for the shot gather record, generating a diving wave equation curve for a particular parameter pair of multiple parameter pairs, and integrating the shot gather record data corresponding to the diving wave equation curve over a selected range of offsets of the shot gather to generate an equivalent linear velocity value for the particular parameter pair and the shot gather record data, selecting, from the equivalent linear velocity values for the plurality of parameter pairs, a greatest equivalent linear velocity value of the equivalent linear velocity values, the greatest equivalent linear velocity value corresponding to a first-arrival parameter pair, and determining, using the first-arrival parameter pair, a set of first-arrival onsets for the selected sub-range of offsets.

Embodiments of the disclosure relate to distributed storage methods and apparatus for travel time tables. One method includes determining a quantity of server sides, determining a server side from all server sides as a primary server side, storing travel time tables corresponding to any seismic work area grids at each moment or travel time tables corresponding to multiple seismic work area grids at each moment, establishing a travel time table storage index in real time according to identifiers of the travel time tables corresponding to the any seismic work area grids at each moment or identifiers of the travel time tables corresponding to the multiple seismic work area grids at each moment, and storing the travel time table storage index in the primary server side.

A retrofit valve shutoff device is provided that comprises a coupling key for coupling with an actuator of a shutoff valve on a fluid supply line, an inertial measurement unit for generating one or more signals in response to arrival of seismic waves, a motor for rotating the coupling key and the actuator of the shutoff valve, and a processing unit for receiving the one or more signals from the inertial measurement unit, analyzing the received signals to determine whether to close the shutoff valve, and sending a signal to the motor to rotate the coupling key and the actuator of the shutoff valve to close the shutoff valve based on the analysis of the received signals.

A method is described for seismic imaging that will produce a seismic image with correctly focused and positioned reflectors. This is accomplished by adding physical geological information to a beam tomography process to generate an updated earth model for the seismic imaging. The method may be executed by a computer system.

Disclosed is a method of characterizing a subsurface volume. The method comprises: extracting a geobody from seismic data arranged within a discretized volume comprising a plurality of cells, the geobody comprising a subset of the plurality of cells, each cell of the subset having one or more properties indicative of a particular fluid phase. The extraction of the geobody comprises: determining a propagation probability value for each cell indicative of the probability that a front will propagate through the cell; beginning from a source within the discretized volume, using the propagation probability value to calculate a traveltime for each cell, the travel time describing the time the front takes to travel from the source point to the cell; and using the traveltimes to extract the geobody from the seismic data.

An exemplary embodiment of the invention relates a method for determining hydraulic permeability of rocks in a subsurface region, the method comprising: in-situ measuring and/or calculating at least one of independent seismic velocities of rocks at different locations in said subsurface region; determining at least one lithological unit in said subsurface region based on the measured seismic velocities; for the at least one lithological unit, acquiring at least one rock sample at an at least one location of said subsurface region; measuring the rock-sample porosity and permeability as functions of stress; measuring and/or calculating of at least one of independent seismic velocities of said at least one rock sample as functions of stress; computing the stiff and compliant porosity of said at least one rock sample; computing numerical coefficients of a given analytical permeability model based on the stiff and compliant porosities; computing coefficients of a given analytical model of an elastic-compliance characteristic of the rock based on the stiff and compliant porosity; computing the stiff and compliant porosity for a plurality of other locations in said subsurface region; and computing the permeability of rocks belonging to said at least one lithological unit, for said at least one location in said subsurface region and for said plurality of other locations in said subsurface region.