Dissertations / Theses on the topic '3D seismic interpretation'

Extracting fault, unconformity, and horizon surfaces from a seismic image is useful for interpretation of geologic structures and stratigraphic features. Although interpretation of these surfaces has been automated to some extent by others, significant manual effort is still required for extracting each type of these geologic surfaces. I propose methods to automatically extract all the fault, unconformity, and horizon surfaces from a 3D seismic image. To a large degree, these methods just involve image processing or array processing which is achieved by efficiently solving partial differential equations.

For fault interpretation, I propose a linked data structure, which is simpler than triangle or quad meshes, to represent a fault surface. In this simple data structure, each sample of a fault corresponds to exactly one image sample. Using this linked data structure, I extract complete and intersecting fault surfaces without holes from 3D seismic images. I use the same structure in subsequent processing to estimate fault slip vectors. I further propose two methods, using precomputed fault surfaces and slips, to undo faulting in seismic images by simultaneously moving fault blocks and faults themselves.

For unconformity interpretation, I first propose a new method to compute a unconformity likelihood image that highlights both the termination areas and the corresponding parallel unconformities and correlative conformities. I then extract unconformity surfaces from the likelihood image and use these surfaces as constraints to more accurately estimate seismic normal vectors that are discontinuous near the unconformities. Finally, I use the estimated normal vectors and use the unconformities as constraints to compute a flattened image, in which seismic reflectors are all flat and vertical gaps correspond to the unconformities. Horizon extraction is straightforward after computing a map of image flattening; we can first extract horizontal slices in the flattened space and then map these slices back to the original space to obtain the curved seismic horizon surfaces.

The fault and unconformity processing methods above facilitate automatic flattening and horizon extraction by providing an unfaulted image with continuous reflectors across faults and unconformities as constraints for an automatic flattening method. However, human interaction is still desirable for flattening and horizon extraction because of limitations in seismic imaging and computing systems, but the interaction can be enhanced. Instead of picking or tracking horizons one at a time, I propose a method to compute a volume of horizons that honor interpreted constraints, specified as sets of control points in a seismic image. I incorporate the control points with simple constraint preconditioners in the conjugate gradient method used to compute horizons.

Integrated analysis of well and geophysical data can provide detailed geologic interpretation of the subsurface in Osage County, Oklahoma. Systems tracts and depositional system successions can be interpreted at marginal seismic resolution using well log motif with seismic reflector character within a depositional context. Shelf-prism and subaqueous, delta-scale clinoforms of Missourian age observed in 3D seismic were interpreted with greater sequence stratigraphic detail when coupled with wireline well logs. The Late Pennsylvanian Midcontinent Sea was thought to be approximately 150 feet average depth across the southern Midcontinent during the Missourian Stage, and deepen towards the Arkoma and Anadarko Basins to the south. Here we show that the Late Pennsylvanian Midcontinent Sea floor was in water depths greater than 600 feet and sloped to the southeast, toward major, southern basins, during the Missourian Stage in Osage County. Shelf-prism and delta scale clinoforms up to 600 and 300 feet of relief, respectively, were observed in paired seismic and well log cross sections, thickness maps, and structure maps dipping northwest at 052° strike, upon a basin floor dipping southeast at 253° strike. Lithologic and sequence stratigraphic interpretation revealed a mixed carbonate-siliciclastic system comprising of delta, offshore shelf, and carbonate buildup depositional systems of mesothem, 3rd order sequence magnitude. The observed succession included: 1) falling stage to lowstand, sand-prone, subaqueous delta, 2) transgressive to highstand offshore shelf and carbonate bank, and 3) falling stage delta. The depositional sucession demonstrates how carbonate banks related spatially to terrigenous sediment input in northeastern Oklahoma during the Late Pennsylvanian because of glacio-eustasy and possible tectonism.

Detection of Morrow A sandstones is a major problem in the exploration of new fields and the characterization of existing fields because they are very thin and laterally discontinuous. The present research shows the advantages of S-wave data in detecting and characterizing the Morrow A sandstone. Full-waveform modeling is done to understand the sandstone signature in P-, PS- and S-wave gathers. The sandstone shows a distinct high-amplitude event in pure S-wave reflections as compared to the weaker P- and PS-wave events. Modeling also helps in understanding the effect of changing sandstone thickness, interbed multiples (generated by shallow high-velocity anhydrite layers) and sidelobe interference effect (due to Morrow shale) at the Morrow A level.

Multicomponent data need proper care while processing, especially the S-wave data which are aected by the near-surface complexity. Cross-spread geometry and 3D FK filtering are effective in removing the low-velocity noise trends. The S-wave data obtained after stripping the S-wave splitting in the overburden show improvement for imaging and reservoir property determination. Individual P- and S-wave attributes as well as their combinations have been analyzed to predict the A sandstone thickness. A multi-attribute map and collocated cokriging procedure is used to derive the seismic-guided isopach of the A sandstone.

Postle Field is undergoing CO₂ flooding and it is important to understand the characteristics of the reservoir for successful flood management. Density can play an important role in finding and monitoring high-quality reservoirs, and to predict reservoir porosity. prestack P- and S-wave AVO inversion and joint P- and S-wave inversion provide density estimates along with the P- and S-impedance for better characterization of the Morrow A sandstone. The research provides a detailed multicomponent processing, inversion and interpretation work flow for reservoir characterization, which can be used for exploration in other parts of the world as well.

Hydrocarbon recovery rates within the Niobrara Shale are estimated as low as 2–8%. These recovery rates are controlled by the ability to effectively hydraulic fracture stimulate the reservoir using multistage horizontal wells. Subsequent to any mechanical issues that affect production from lateral wells, the variability in production performance and reserve recovery along multistage lateral shale wells is controlled by the reservoir heterogeneity and its consequent effect on hydraulic fracture stimulation efficiency. Using identical stimulation designs on a number of wells that are as close as 600ft apart can yield variable production and recovery rates due to inefficiencies in hydraulic fracture stimulation that result from the variability in elastic rock properties and in-situ stress conditions.

As a means for examining the effect of the geological heterogeneity on hydraulic fracturing and production within the Niobrara Formation, a 3D geomechanical model is derived using geostatistical methods and volumetric calculations as an input to hydraulic fracture stimulation. The 3D geomechanical model incorporates the faults, lithological facies changes and lateral variation in reservoir properties and elastic rock properties that best represent the static reservoir conditions pre-hydraulic fracturing. Using a 3D numerical reservoir simulator, a hydraulic fracture predictive model is generated and calibrated to field diagnostic measurements (DFIT) and observations (microseismic and 4D/9C multicomponent time-lapse seismic). By incorporating the geological heterogeneity into the 3D hydraulic fracture simulation, a more representative response is generated that demonstrate the variability in hydraulic fracturing efficiency along the lateral wells that will inevitability influence production performance.

Based on the 3D hydraulic fracture simulation results, integrated with microseismic observations and 4D/9C time-lapse seismic analysis (post-hydraulic fracturing & post production), the variability in production performance within the Niobrara Shale wells is shown to significantly be affected by the lateral variability in reservoir quality, well and stage positioning relative to the target interval, and the relative completion efficiency. The variation in reservoir properties, faults, rock strength parameters, and in-situ stress conditions are shown to influence and control the hydraulic fracturing geometry and stimulation efficiency resulting in complex and isolated induced fracture geometries to form within the reservoir. This consequently impacts the effective drainage areas, production performance and recovery rates from inefficiently stimulated horizontal wells.

The 3D simulation results coupled with the 4D seismic interpretations illustrate that there is still room for improvement to be made in optimizing well spacing and hydraulic fracturing efficiency within the Niobrara Formation. Integrated analysis show that the Niobrara reservoir is not uniformly stimulated. The vertical and lateral variability in rock properties control the hydraulic fracturing efficiency and geometry. Better production is also correlated to higher fracture conductivity. 4D seismic interpretation is also shown to be essential for the validation and calibration hydraulic fracture simulation models. The hydraulic fracture modeling also demonstrations that there is bypassed pay in the Niobrara B chalk resulting from initial Niobrara C chalk stimulation treatments. Forward modeling also shows that low pressure intervals within the Niobrara reservoir influence hydraulic fracturing and infill drilling during field development.