Abstract: A method and system for estimating a full elastic tensor. The method may comprise taking a measurement for compressional wave sonic data and cross-dipole shear data with a sonic logging tool at a first location as cross-dipole data, processing the compressional wave sonic data to produce a compressional wave slowness (P), and processing the cross-dipole shear data to produce a fast horizontal polarized shear wave slowness (SH) and a slow quazi-vertical shear wave slowness (qSV) as a function of depth. The method may further comprise setting an initial guess for at least five constants of the full elastic tensor for Vertical Transversely Isotropy (VTI) symmetry, determining a modeled slowness surfaces from the full elastic tensor, and comparing the modeled slowness surfaces with measured values of the P, the SH, and the qSV. The method may be performed by a system comprising a sonic logging tool and an information handling system.

Abstract: Systems, methods, and apparatuses for generating a subsurface image using diffraction energy information are disclosed. The systems, methods, and apparatuses may include converting a shot gather into one or more plane-wave gather using a Radon transform. The plane-wave gathers may be extrapolated into source-side wavefields and receiver-side wavefields and further generate a pseudo dip-angle gather. The diffraction energy information may be extracted from the pseudo dip-angle gather, and an image containing subsurface features may be generated from the extracted diffraction energy information. The receiver-side wavefields may be decomposed using a recursive Radon transform.

Abstract: A hydrophone array tool as described herein is configured to locate leakages throughout a borehole with improved accuracy using acoustic beamforming techniques with acoustic velocity estimation. An acoustic beamforming processor generates an initial beamforming map and corresponding initial estimated leakage location using acoustic measurements throughout the borehole. The acoustic beamforming processor generates additional beamforming maps at the initial estimated leakage location, each additional beamforming map corresponding to an acoustic velocity within a range of anticipated acoustic velocities. An acoustic velocity estimator determines an acoustic velocity corresponding to a beamforming map with spatial statistics that indicate a most prominent leakage location. The acoustic beamforming processor updates the leakage location according to this beamforming map for improved accuracy.

Abstract: A seismic time-frequency analysis method based on generalized Chirplet transform with time-synchronized extraction, which has higher level of energy aggregation in the time direction and can better describe and characterize the local characteristics of seismic signals, and is applicable to the time-frequency characteristic representation of both harmonic signals and pulse signals, comprising the steps of processing generalized Chirplet transform with time-synchronized extraction for each seismic signal to obtain a time spectrum by: carrying out generalized Chirplet transform, calculating group delay operator and carrying out time-synchronized extraction on seismic signals, thereby the boundary and heterogeneity structure of the rock slice are more accurately and clearly shown and subsequence seismic analysis and interpretation are facilitated.

Abstract: A system and method can be used for to calibrating time-lapse seismic volumes by cross-migration rescaling and reorientation for use in determining optimal wellbore placement or production in a subsurface environment. Certain aspects include methods for cross-migration of data sets processed using different migration techniques. Pre-processing of the data sets, optimization of rescaling and reorientation, and identification of adjustment parameters associated with minimum global error can be used to achieve a time-dependent formation data set that addresses error in all input data sets.

Abstract: A method can include acquiring imagery of an exposed surface of the Earth; generating a multi-dimensional model based at least in part on the imagery; generating synthetic seismic data utilizing the multi-dimensional model; acquiring seismic data of a subsurface region of the Earth; performing a search that matches a portion of the acquired seismic data and a portion of the synthetic seismic data; and characterizing the subsurface region of the Earth based at least in part on the portion of the synthetic seismic data.

Abstract: Methods and devices for seismic exploration of an underground structure apply W2-based full-wave inversion to transformed synthetic and seismic data. Data transformation ensures that the synthetic and seismic data are positive definite and have the same mass using an adaptive normalization. This approach yields superior results particularly when the underground structure includes salt bodies.

Abstract: Optimizing seismic to depth conversion to enhance subsurface operations including measuring seismic data in a subsurface formation, dividing the subsurface formation into a training area and a study area, dividing the seismic data into training seismic data and study seismic data, wherein the training seismic data corresponds to the training area, and wherein the study seismic data corresponds to the study area, calculating target depth data corresponding to the training area, training a machine learning model using training inputs and training targets, wherein the training inputs comprise the training seismic data, and wherein the training targets comprise the target depth data, computing, by the machine learning model, output depth data corresponding to the study area based at least in part on the study seismic data; and modifying one or more subsurface operations corresponding to the study area based at least in part on the output depth data.

Abstract: Systems and computer readable media are described that actuate at least one marine seismic source according to a constrained sequence. The sequence exhibits an actuation time or distance interval between each actuation. The actuation time or distance interval corresponds to the sum of a nominal time or distance and a respective dither time or dither distance for each actuation. The sequence is constrained such that differences between consecutive dither times or dither distances remain within a threshold dither difference. Constraining the sequence according to the threshold dither difference enables increased bottom speeds for the source (i.e., increased speeds of the source relative to the seafloor), while still maintaining at least a minimum actuation time or distance interval for the source, taking into account both the nominal time or distance and the respective dither times or dither differences.

Abstract: A method for subsurface fault extraction using undirected graphs is provided. Extracting faults in the subsurface may assist in various stages of geophysical prospecting. To that end, an undirected graph may be used in order to identify distinctive fault branches in the subsurface. Fault probability data, from seismic data, may be used to establish connections in the undirected graph. Thereafter, some of the connections in the undirected graph may be removed based on analyzing one or more attributes, such as dip, azimuth, or context, associated with the connections or nodes associated with the connections. After which, the undirected graph may be analyzed in order to extract the faults in the subsurface.

Abstract: System and methods for predicting time-dependent rock properties are provided. Seismic data for a subsurface formation is acquired over a plurality of time intervals. A value of at least one rock property of the subsurface formation is calculated for each of the plurality of time intervals, based on the corresponding seismic data acquired for that time interval. At least one of a trend or a spatio-temporal relationship in the seismic data is determined based on the value of the at least one rock property calculated for each time interval. A value of the at least one rock property is estimated for a future time interval, based on the determination. The estimated value of the at least one rock property is used to select a location for a wellbore to be drilled within the subsurface formation. The wellbore is then drilled at the selected location.

Abstract: A method for predicting a pressure window for drilling a borehole in a formation includes: obtaining a pore pressure related data value of the formation using a data acquisition tool; predicting pore pressure uncertainty from the pore pressure related data value of the formation using a processor; estimating uncertainty of a pressure window for drilling fluid using the predicted pore pressure uncertainty using a processor; and applying the estimated uncertainty to the pressure window to provide a modified pressure window using a processor.

Abstract: Methods, apparatuses, and computer-readable media utilize data stacking to facilitate identification and/or correction of data interpretation conducted for a subsurface formation. Related data sets, such as well logs, may be displayed along with markers representing a common entity in the related data sets, such as formation features in a surface formation, and a visualization of stacked data may be generated and centered on the markers to highlight mis-alignment of any of the markers.

Abstract: Systems and methods are provided for reducing error in data compression and decompression when data is transmitted over low bandwidth communication links, such as satellite links. Embodiments of the present disclosure provide systems and methods for variable block size compression for gridded data, efficiently storing null values in gridded data, and eliminating growth of error in compressed time series data.

Abstract: Disclosed are systems and methods for processing broadband acoustic signals acquired by a plurality of acoustic sensors, using an array-signal-processing technique to compute fused-signal maps in the frequency domain for a plurality of frequency bins. In accordance with various embodiments, the fused-signal maps are combined across the frequency bins, with respective weightings that are based on eigenvalues of covariance matrices computed for the plurality of frequency bins. The combined maps can be used to locate an acoustic source in a wellbore.

Abstract: A method of underwater tsunami detection includes detecting a trigger event using disruption of at least one of a plurality of hard disk drives (HDDs), each different one of the plurality of HDDs in a different one of a plurality of autonomous underwater vehicles (AUVs). A time and location of each of the at least one HDD for the trigger event is logged. Based on at least one of the HDD disruptions, times, and locations of the at least one HDD of the plurality of HDDs, a size, strength, and direction of a tsunami caused by the trigger event is determined. Information regarding the tsunami is transmitted to a monitoring station.

Abstract: An inspection tool and associated methods for testing physical and electrical properties of a perforating gun and sending the perforating gun to a wellbore site with at least one of an electrical property and a dimension that has been previously verified. The perforating gun may be received in a perforating gun holder positioned between a first connecting portion and second connecting portion. The first connecting portion may be moved from a first position to a second position at which the at least one of the electrical property and the dimension is measured. Upon verification that the at least one of the electrical property and the dimension is within a predetermined specification, the perforating gun may be sent to the wellbore site.

Abstract: A direction-finding and positioning system of electromagnetic emission of coal or rock fracture includes a three-axis electromagnetic sensor array, a signal acquisition module and a direction-finding and positioning terminal; the three-axis electromagnetic sensor array is composed of at least four three-axis electromagnetic sensors configured to synchronously sense magnetic field strength in three-axis direction based on a tunnel magneto resistance technology, and obtain a real magnetic field vector in space by measuring; the signal acquisition module is configured to acquire magnetic field vector variable information of multiple measuring points in real-time, and after extracting magnetic field vector variable parameters, transmitting the magnetic field vector variable parameters to the direction-finding and positioning terminal; the direction-finding and positioning terminal is configured to perform direction-finding and positioning calculations according to the magnetic field vector variable parameters rec

Abstract: A method for providing a calibrated rock-physics model of a subsoil. First, a geological model of the subsoil comprising a grid made of cells, associated with a rock-physics parameter is obtained. A group of cells forming a calibration body is selected in the grid. The calibration body corresponds to a region of the subsoil having substantially homogenous rock-physics parameter values. Finally, an adjustable constant parameter in a physical equation expressing a relationship between the petro-physical parameter and a petro-elastic parameter in the calibration body is calibrated so as to reduce a mismatch between the petro-elastic parameter estimated using the physical equation and a petro-elastic parameter value determined from inverted seismic data, the calibrated physical equation providing a calibrated rock-physics model of the subsoil in the calibration body.

Abstract: Delay calibration for digital signal chains of SFCW systems is disclosed. An example calibration method includes receiving a burst with a test pulse, the burst having a duration of L clock cycles; receiving a trigger indicative of a time when the burst was transmitted; generating a digital signal indicative of the received burst; for each of L clock cycles, computing a moving average of a subset of digital samples and an amplitude for each average; identifying one moving average for which the computed amplitude is closest to an expected amplitude; identifying the clock cycle of the identified moving average; and updating at least one delay to be applied in digital signal processing of received bursts based on a difference between the trigger and the identified clock cycle. The delay may be used for selecting digital samples of the received signal that contain valid data for performing further data processing.

Abstract: A seabed object detection system is provided. The system can include a receiver array including streamers. The system can include a plurality of receivers coupled with the streamers. The system can include a receiver array cross-cable to couple with the first streamer and to couple with the second streamer. The receiver array cross-cable can be disposed at a first depth of a body of water. The system can include a first diverter and a second diverter coupled with the receiver array cross-cable. The system can include a source array including a first source and a second source. The source array can be coplanar to the receiver array. The system can include a source array cross-cable to couple with the first source and to couple with the second source, the source array cross-cable disposed at a second depth of the body of water.

Abstract: Methods are taught for creating training data for a learning algorithm, training the learning algorithm with the training data and using the trained learning algorithm to suggest domain names to users. A domain name registrar may store activities of a user on a registrar website. Preferably, domain name searches, selected suggested domain names and domain names registered to the user are stored as the training data in a training database. The training data may be stored so that earlier activities act as inputs to the learning algorithm while later activities are the expected outputs of the learning algorithm. Once trained, the learning algorithm may receive activities of other users and suggest domain names to the other users based on their activities.

Abstract: Accordingly, there are disclosed herein geologic modeling methods and systems employing reference-based inversion of seismic image volumes. An illustrative method embodiment includes: (a) obtaining a measured seismic image volume; (b) determining a reference seismic image volume based on a reference model; (c) deriving a synthesized seismic image volume from a geologic model; (d) detecting at least one geologic model region where the synthesized seismic image volume and the measured seismic image volume are mismatched; (e) finding a reference model region where the reference seismic image volume best matches the measured seismic image volume; (f) replacing content of the at least one geologic model region with content of the reference model region to obtain an improved geologic model; and (g) outputting the improved geologic model.

Abstract: A computer-based method for real-time communication authorization includes receiving, from a first communication device, a communication request, verifying, with a verification engine, a pre-approval status of the communication request, storing the communication request in an approval queue if the communication authorization the pre-approval status is set to false, issuing an alert to the authorization device, and receiving one or more authorization parameters from an authorization device.

Abstract: Reducing casing wave effects on sonic logging data by positioning two or more receivers in a borehole in a subsurface formation; receiving, at two or more receivers in a borehole in a subsurface formation, a first signal associated with a first acoustic signal originating from a first transmitter position; receiving, at the two or more receivers, a second signal associated with a second acoustic signal originating from a second transmitter position; creating a dataset based on the first signal and the second signal; identifying casing wave signals in the dataset based at least in part on the second signal; calculating inverse-phase casing wave signals based at least in part on the casing wave signals and the second signal; and reducing effects of the casing wave signals on the dataset using the inverse-phase casing wave signals.

Abstract: In a first step, a defined scope value is selected for each of a plurality of hydrodynamic input parameters. A simulated topographical result is generated using the selected scope values and a forward model. A detailed seismic interpretation is generated to represent specific seismic features or observed topography. A calculated a misfit value representing a distance between the simulated topographical result and a detailed seismic interpretation is minimized. An estimated optimized sand ratio and optimized hydrodynamic input parameters are generated. In a second step, a genetic algorithm is used to determine a proportion of each grain size in the estimated optimized sand ratio. A misfit value is used that is calculated from thickness and porosity data extracted from well data and a simulation result generated by the forward model to generate optimized components of different grain sizes. Optimized hydrodynamic input parameters and optimized components of different grain sizes are generated.

Abstract: Methods and related systems for charging a refrigerant into a climate control system. In an embodiment, the method includes (a) coupling a storage tank to a refrigerant loop of the climate control system through a charging valve; (b) opening and closing the charging valve in a plurality of cycles; and (c) flowing refrigerant from the storage tank to the refrigerant loop through the charging valve when the charging valve is open, during (b). In addition, the method includes (d) determining a detected saturated temperature of the refrigerant within the refrigerant loop after each cycle of the plurality of cycles; and (e) adjusting an amount of time that the charging valve is open during each cycle of the plurality of cycles during (b) as a function of the detected saturated temperature from a previous cycle of the plurality of cycles.

Abstract: A method for evaluating accuracy in positioning a receiver point, which is associated with at least one shot point, and for which a derived position data is obtained, wherein a pair of the receiver point and a respective shot point is associated with a characteristic parameter which includes an offset, a velocity of first arrival wave and a first arrival time, the method comprising: step S10, constructing a residual vector associated with the receiver point and the respective shot point based on the characteristic parameter; step S12, determining a characterization parameter of the derived position data based on the residual vector; and step S14, evaluating accuracy of the derived position data based on the characterization parameter. An apparatus for evaluating accuracy in positioning a receiver point is also provided.

Abstract: A method and apparatus for generating a fluid saturation model for a subsurface region. One example method generally includes obtaining a model of the subsurface region; for each of a plurality of fluid types: flooding the subsurface region model with the fluid type to generate a flood model; and running a trial petrophysical inversion with the flood model to generate a trial petrophysical model; identifying potential fluid contact regions in the trial petrophysical models; partitioning the subsurface region model at the identified potential fluid contact regions; and constructing the fluid saturation model from the partitioned subsurface region model.

Abstract: Various aspects described herein relate to a machine learning based signal recovery. In one example, a computer-implemented method of noise contaminated signal recovery includes receiving, at a server, a first signal including a first portion and a second portion, the first portion indicative of data collected by a plurality of sensors, the second portion representing noise; performing a first denoising process on the first signal to filter out the noise to yield a first denoised signal; applying a machine learning model to determine a residual signal indicative of a difference between the first signal and the first denoised signal; and determining a second signal by adding the residual signal to the first denoised signal, the second signal comprising (i) signals of the first portion with higher magnitudes than the noise in the second portion, and (ii) signals of the first portion having lower magnitudes than the noise in the second portion.

Abstract: A method for accurately measuring one or more network performance statistics during one or more sampling periods comprises processing a received data packet to identify flow information and time of arrival; transmitting the data packet, along with metadata comprising the flow information and the time of arrival, to a sample processing module; replaying the reception of the data packet based on the relative speed of the transmission link and the databus and the time of arrival; incrementing a number of bytes received for the data packet for a first sampling period; starting a second sampling period at a predetermined time and incrementing a number of bytes received for the data packet for the second sampling period, if the replaying is not completed. The number of bytes received for the data packet for the first and second sampling periods are processed into first and second sampling period statistics, respectively.

Abstract: In one embodiment, a method of processing petroleum extraction data includes receiving petroleum extraction data from one or more data sources, wherein the petroleum extraction data includes structured data and unstructured data, filtering the petroleum extraction data by applying one or more statistical methods to generate filtered petroleum extraction data including transactional data and non-transactional data, and verifying the filtered petroleum extraction data according to one or more attributes. The method further includes calculating one or more parameters from the transactional data at a first frequency and one or more parameters from the non-transactional data at a second frequency, and correlating the one or more parameters from the transactional data and the non-transactional data based on one or more attributes to generate correlated data. The method also includes aggregating and storing the correlated data in a data structure, and performing one or more analytic processes on the aggregated data.

Abstract: A method of determining an arrival-time of a first seismic event in a seismic data set including, obtaining the seismic data set and an initial seismic velocity model, and determining an updated seismic velocity model based on the seismic data set. Furthermore, the method includes determining a simulated arrival-time of the first seismic event based on the updated seismic velocity model and defining a predicted time-window based on the simulated arrival-time of the first seismic event, and picking the arrival-time of the first seismic event in the seismic data set based on the predicted time-window.

Abstract: Processes and systems for generating images of a subterranean formation from recorded seismic data obtained in a marine survey are described. Processes and systems compute reverse-time receiver-motion-corrected upgoing and downgoing pressure wavefields at different locations of corresponding upgoing and downgoing observation levels based on the recorded seismic data. The reverse-time receiver-motion-corrected upgoing and downgoing pressure wavefields are time forwarded and extrapolated to obtain a corresponding receiver-motion-corrected upgoing and downgoing pressure wavefields at locations of a static observation level. An image of the subterranean formation is generated based at least in part on the receiver-motion-corrected upgoing pressure wavefield and the receiver-motion-corrected downgoing pressure wavefield.

Abstract: Systems and methods for interpreting and visualizing multi-Z horizons from seismic data are disclosed. A two-dimensional (2D) representation of seismic data is displayed via a graphical user interface (GUI). User input is received via the GUI for interpreting a multi-Z horizon within a portion of the displayed 2D representation. The user's input is tracked relative to displayed 2D representation within the GUI. Based on the tracking, each of a plurality of surfaces for the multi-Z horizon is determined. At least one intersection point between the multi-Z horizon surfaces is identified. A depth position for each surface relative to other surfaces is determined. The 2D representation of the seismic data is dynamically updated to include visual indications for the plurality of surfaces and the intersection point(s), based on the depth position of each surface, where the visual indications use different visualization styles to represent the surfaces and intersection point(s).

Abstract: Methods and systems for optimizing drilling operations in a wellbore using a drill string are provided. The methods and systems include measuring a first formation characteristic with at least one sensor, measuring a second formation characteristic by means of a hydraulic test, the at least one second formation characteristic being different from the at least one first formation characteristic, generating a model to represent a formation around the wellbore, the model incorporating the first formation characteristic and the second formation characteristic, and performing a drilling operation based on the generated model.

Abstract: A method for partial differential equation inversion, the method including receiving measured data do; selecting an objective function having first and second measures N1 and N2, wherein the objective function depends on three independent variables V, u, and f, V being a perturbation of a wave equation operator L from a background operator L0, u being a wavefield that satisfies the wave equation operator L, and f being a source function that describes the source of the waves; optimizing with a processor the objective function by finding a minimum or a maximum using the inversion; calculating with the processor solutions V*, u*, and f* of the three independent variables V, u, and f; and generating with the processor an image of an object based on the solutions V*, u*, and f*.

Abstract: An interface autonomously collects input data to control the operation or travel path of an object, includes a plurality of independently operating sensors that detect the occurrence of an event. A master fractal controller initiates a space-variant fractalization protocol in order to collect the input data related to the detected event. In response to the input, the master fractal controller sends control signals to a servient controller to control the operation or travel path of the object. The master fractal controller initiates the space-variant fractalization protocol by identifying or defining a fractal space, and by further initiating sub-servient fractals in sub-servient spaces, within which an intersection of a graphical function layer that includes an executable grid with a graphical data layer causes a command or function be executed within the sub-servient fractal space.

Abstract: In an example implementation, first seismic energy is generated using first seismic sources positioned on an earth's surface. First data including measurements of the first seismic energy is obtained from first geophones positioned at a first depth below the earth's surface. Second data including measurements of the first seismic energy is obtained from second geophones positioned on the earth's surface. Second seismic energy is generated using second seismic sources positioned on an earth's surface and proximal to the second geophones. Third data including measurements of the second seismic energy is obtained from third geophones positioned at the first depth below the earth's surface. A propagation of the first seismic energy along a first path is estimated based on the first, second and third data. One or more characteristics of the target are determined based on the estimate.

Abstract: Systems and methods are provided for reducing error in data compression and decompression when data is transmitted over low bandwidth communication links, such as satellite links. Embodiments of the present disclosure provide systems and methods for variable block size compression for gridded data, efficiently storing null values in gridded data, and eliminating growth of error in compressed time series data.

Abstract: A method and apparatus for forecasting output power of wind turbine in a wind farm. The present invention provides a method for forecasting output power of a wind turbine in a wind farm, including: generating a corrected data set based on environmental data collected from at least one sensor in the wind farm; correcting a weather forecasting model by using the corrected data set; obtaining a forecast value of wind information at the wind turbine based on the corrected weather forecasting model; and forecasting the output power of the wind turbine based on the forecast value and a power forecasting model.

Abstract: Disclosed are method and apparatus for obtaining residual gravity anomaly. The method comprises obtaining Bouguer gravity anomaly of a target region, determining a first pre-set range corresponding to each sampling point in the target region, obtaining a first regional field value of each sampling point in the first pre-set range through a surface fitting method based on the coordinate and field value of the sampling point, traversing the target region to obtain first regional gravity anomaly of the target region according to the first regional field values of the sampling points in the target region, and obtaining first residual gravity anomaly of the target region according to the Bouguer gravity anomaly and the first regional gravity anomaly. The residual gravity anomaly thus obtained is more accurate, thus enabling an accurate prediction of an underground geological body.

Abstract: Strata surfaces can be identified in well logs and correlated across the well logs taking into account manual corrections. For example, a computing device can receive well logs. The computing device can determine multiple stratum-surface correlations based on the well logs. Then, the computing device can receive user input indicating a correction to a particular stratum-surface correlation. Based on the correction to the particular stratum-surface correlation, the computing device can update some or all of the other stratum-surface correlations.

Abstract: Systems and methods are provided to correct seismic data for the undesired effects caused by near surface velocity variations. In one embodiment, a method includes receiving travel time data for a near surface region and estimating an initial velocity model for the near surface region using the travel time data. The method can include updating the velocity model by performing an inversion iteration of including inversion of travel times to estimate slowness. The process can also include calculating at least one long-wavelength static for the near surface region. The long-wavelength statics may be used to correct for undesired effects caused by near surface velocity variations.

Abstract: Methods for full waveform inversion (FWI) in the midpoint-offset domain include using a computer system to sort seismic traces into common midpoint-offset bins (XYO bins). For each XYO bin, a linear moveout correction is applied to a collection of seismic traces within the XYO bin. The collection of seismic traces is stacked to form a pilot trace. The computer system determines a surface-consistent residual static correction for each seismic trace. The computer system determines that the surface-consistent residual static correction for each seismic trace is less than a threshold. Responsive to the determining that the surface-consistent residual static correction is less than the threshold, the computer system stacks the collection of seismic traces to provide the pilot trace. The computer system groups the pilot traces for the XYO bins into a set of virtual shot gathers. The computer system performs one-dimensional FWI based on each virtual shot gather.

Abstract: A method and apparatus for machine learning for use with automated seismic interpretation include: obtaining input data; extracting patches from a pre-extraction dataset based on the input data; transforming data of a pre-transformation dataset based on the input data and geologic domain knowledge and/or geophysical domain knowledge; and generating augmented data from the extracted patches and the transformed data. A method and apparatus for machine learning for use with automated seismic interpretation include: a data input module configured to obtain input data; a patch extraction module configured to extract patches from a pre-extraction dataset that is based on the input data; a data transformation module configured to transform data from a pre-transformation dataset that is based on the input data and geologic domain knowledge and/or geophysical domain knowledge; and a data augmentation module configured to augment data from the extracted patches and the transformed data.

Abstract: Systems and methods develop a three-dimensional model of a subterranean formation based on vertical seismic profiles at a plurality of well locations.

Abstract: A new method for iteratively picking the seismic first breaks and conducting imaging of the near-surface velocity structures in an iterative fashion is provided that the first-break picks of the input seismic data are applied to image the near-surface velocity structures and the calculated travel times associated with the updated velocity structures are applied to help refine the first-break picks in the first break picking process until first-break picks satisfy a number of quality control criteria, statics solutions are optimized, and the near surface imaging reaches an acceptable data misfit. This invention produces a velocity model that can be used for near surface statics corrections or for the prestack depth migration.

Abstract: Provided in this disclosure are systems and methods for selection of petrophysical techniques to model reservoirs. Reservoir properties were calculated using two techniques—a deterministic technique and an optimizing petrophysics technique, and simulation models were developed. The technique that yielded a simulation model that aligned more accurately and consistently with the reservoir data was selected to further develop the reservoir model.

Abstract: Systems and methods for editing multi-Z horizons interpreted from seismic data are provided. A multi-Z horizon having a plurality of surfaces is visualized within a two-dimensional (2D) representation of seismic data displayed via a graphical user interface (GUI) of an application executable at a computing device of a user. Input is received via the GUI from the user for editing one or more of the plurality of surfaces of the multi-Z horizon within a current view of the displayed 2D representation of the seismic data. A location of the received input relative to each of the plurality of surfaces within the current view is determined. The one or more surfaces of the multi-Z horizon are modified based on the location of the received input within the current view. The visualization of the multi-Z horizon within the GUI is updated, based on the modified one or more surfaces.