A method for determining a favorable time window of an infill well of an unconventional oil and gas reservoir, which comprises the following steps: S1, establishing a three-dimensional geological model with physical properties and geomechanical parameters; S2, establishing a natural fracture network model in combination with indoor core-logging-seismic monitoring; S3, calculating complex fractures in hydraulic fracturing of parent wells; S4, establishing an unconventional oil and gas reservoir model and calculating a current pore pressure field; S5, establishing a dynamic geomechanical model and calculating a dynamic geostress field; S6, calculating complex fractures in horizontal fractures of the infill well in different production times of the parent wells based on pre-stage complex fractures and the current geostress field; S7, analyzing a microseismic event barrier region and its dynamic changes in infill well fracturing; and S8, analyzing the productivity in different infill times, and determining an infill time window.

A method and device of identifying a fracture are provided in the embodiments of the present application. The method comprises: determining three components of structure quantification for each data point in a seismic data volume; constructing a structure quantification matrix of the data point according to the three components of structure quantification for each of the data points; determining feature value and feature vector of the structure quantification matrix of each of the data points; determining fracture attribute value of the data point according to the feature value and feature vector of the structure quantification matrix of each of the data points; constructing a data volume of the fracture attribute according to the fracture attribute values of respective data points; and performing a fracture extraction for the data volume of the fracture attribute according to the feature vectors of the structure quantification matrix of the respective data points. The embodiments of the present application can improve the accuracy of identifying a minor fracture, so as to realize an effective identification of the minor fracture.

Systems and methods for identifying potential hydrocarbon traps in a subterranean region can include: receiving seismic data of the subterranean region, the seismic data acquired by at least one seismic sensor, the seismic data indicating positions of physical barriers to hydrocarbon flow in the subterranean region and using anisotropic lateral and upward erosion to identify possible locations of hydrocarbons.

Systems and methods for generating a fractured reservoir model include: receiving a seismic dataset of a surveyed subsurface; identifying a dynamic response of each parameter; selecting a subset of parameters from the set of parameters based on the dynamic response of each parameter; sampling an outer boundary of a parameter uncertainty domain; adjusting the range of values associated with each parameter of the subset based on the sampling; generating a geo-model based on the adjusted range of values associated with each parameter of the subset; generating a discrete fracture network model based on the geo-model; generating a scenario of a simulation model based on the discrete fracture network model and the geo-model; performing a forward simulation based on the scenario of the simulation model; determining that a misfit of the forward simulation is below a threshold by evaluating an objective function; and producing a model based on the forward simulation.

A seismic dataset and a task to be performed with the seismic dataset may be received. A representative seismic line representative of the seismic dataset may be generated. The representative seismic line may include pixel data representative of the seismic dataset. Based on the representative seismic line, the task may be performed. The task may include at least finding an analogous geological region by searching for an analogous seismic dataset existing in a seismic database by comparing the representative seismic line with the analogous seismic dataset's representative seismic line.

The invention pertains to a method for estimating faults in a three-dimensional seismic image block. Directrices are generated within respective first cross-sections of the seismic image block based on points selected by a user. Similarly, generatrices are generated within respective second cross-sections of the block based on points selected by the user. The user inputs relationships between directrices and generatrices. A fault is estimated within the seismic image block as a surface including at least one directric and at least one generatrix having a relationship therebetween.

Methods, systems, and apparatus, including computer programs encoded on a computer storage medium, to generate a custom seismic surface and volume attribute. In one aspect, a method includes receiving a seismic cube and a seismic surface, and the seismic cube includes traces recorded at receivers deployed to collect seismic data. The seismic surface is picked on the seismic cube. Seismic wavelets are extracted with a selected length from the seismic cube along an intersection with the seismic surface for each spatial coordinate associated with the seismic surface. A reference wavelet is determined. A surface attribute map is generated based on comparing each of the seismic wavelets to the reference wavelet. A productivity of the seismic surface is evaluated using the surface attribute map.

A system and method is provided for restoring a 3D tomographic model of the Earth's subsurface geology from the present-day to a past restoration time. Whereas at the present time all faults represent active discontinuities, at a past restoration time some faults have not yet formed. Accordingly, the restored model divides the fault network into -active faults (discontinuous surfaces for faults that intersect the layer deposited at the past restoration time) and -inactive faults (continuous surfaces for faults that do not intersect the layer deposited at the past restoration time). A new 3D restoration transformation is also provided that uses linear geological constraints to process the restoration model in less time and generate more accurate geological images.

A system and method is provided for restoring a 3D tomographic model of the Earth's subsurface geology from the present-day to a past restoration time. Whereas at the present time all faults represent active discontinuities, at a past restoration time some faults have not yet formed. Accordingly, the restored model divides the fault network into -active faults (discontinuous surfaces for faults that intersect the layer deposited at the past restoration time) and -inactive faults (continuous surfaces for faults that do not intersect the layer deposited at the past restoration time). A new 3D restoration transformation is also provided that uses linear geological constraints to process the restoration model in less time and generate more accurate geological images.

A system and method is provided for restoring a 3D tomographic model of the Earth's subsurface geology from the present-day to a past restoration time. Whereas at the present time all faults represent active discontinuities, at a past restoration time some faults have not yet formed. Accordingly, the restored model divides the fault network into -active faults (discontinuous surfaces for faults that intersect the layer deposited at the past restoration time) and -inactive faults (continuous surfaces for faults that do not intersect the layer deposited at the past restoration time). A new 3D restoration transformation is also provided that uses linear geological constraints to process the restoration model in less time and generate more accurate geological images.