Cross-equalization and prestack seismic inversion of time-lapse land seismic data, Central Thrace Basin, Turkey

2017 ◽  
Vol 36 (7) ◽  
pp. 588-597
Author(s):  
Sait Baytok ◽  
Burcu Selek ◽  
Hüseyin Özdemir
2014 ◽  
Vol 2 (4) ◽  
pp. T255-T271 ◽  
Author(s):  
Roderick Perez Altamar ◽  
Kurt Marfurt

Differentiating brittle and ductile rocks from surface seismic data is the key to efficient well location and completion. Brittleness average estimates based only on elastic parameters are easy to use but require empirical calibration. In contrast, brittleness index (BI) estimates are based on mineralogy laboratory measurements and, indeed, cannot be directly measured from surface seismic data. These two measures correlate reasonably well in the quartz-rich Barnett Shale, but they provide conflicting estimates of brittleness in the calcite-rich Viola, Forestburg, Upper Barnett, and Marble Falls limestone formations. Specifically, the BI accurately predicts limestone formations that form fracture barriers to be ductile, whereas the brittleness average does not. We used elemental capture spectroscopy and elastic logs measured in the same cored well to design a 2D [Formula: see text] to brittleness template. We computed [Formula: see text] and [Formula: see text] volumes through prestack seismic inversion and calibrate the results with the [Formula: see text] template from well logs. We then used microseismic event locations from six wells to calibrate our prediction, showing that most of the microseismic events occur in the brittle regions of the shale, avoiding more ductile shale layers and the ductile limestone fracture barriers. Our [Formula: see text] to brittleness template is empirical and incorporates basin- and perhaps even survey-specific correlations of mineralogy and elastic parameters through sedimentation, oxygenation, and diagenesis. We do not expect this specific template to be universally applicable in other mudstone rock basins; rather, we recommend interpreters generate similar site-specific templates from logs representative of their area, following the proposed workflow.


2012 ◽  
Vol 466-467 ◽  
pp. 400-404
Author(s):  
Jin Zhang ◽  
Huai Shan Liu ◽  
Si You Tong ◽  
Lin Fei Wang ◽  
Bing Xu

Elastic impedance (EI) inversion is one of the prestack seismic inversion methods, which can obtain P-wave and S-wave velocity, density, Poisson ratio, Lame coefficients and other elastic parameters. But there have been many EI formulas nowadays, so which formula should be used in inversion is an urgent problem. This paper divides these formulas into two categories, and use several forward modeling to test the accuracy of these EI formulas. It shows that using the first kind of EI formulas in near offset seismic data can get high precision results.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. M73-M87 ◽  
Author(s):  
Alvaro Rey ◽  
Eric Bhark ◽  
Kai Gao ◽  
Akhil Datta-Gupta ◽  
Richard Gibson

We have developed an efficient approach of petroleum reservoir model calibration that integrates 4D seismic surveys together with well-production data. The approach is particularly well-suited for the calibration of high-resolution reservoir properties (permeability) because the field-scale seismic data are areally dense, whereas the production data are effectively averaged over interwell spacing. The joint calibration procedure is performed using streamline-based sensitivities derived from finite-difference flow simulation. The inverted seismic data (i.e., changes in elastic impedance or fluid saturations) are distributed as a 3D high-resolution grid cell property. The sensitivities of the seismic and production surveillance data to perturbations in absolute permeability at individual grid cells are efficiently computed via semianalytical streamline techniques. We generalize previous formulations of streamline-based seismic inversion to incorporate realistic field situations such as changing boundary conditions due to infill drilling, pattern conversion, etc. A commercial finite-difference flow simulator is used for reservoir simulation and to generate the time-dependent velocity fields through which streamlines are traced and the sensitivity coefficients are computed. The commercial simulator allows us to incorporate detailed physical processes including compressibility and nonconvective forces, e.g., capillary pressure effects, while the streamline trajectories provide a rapid evaluation of the sensitivities. The efficacy of our proposed approach was tested with synthetic and field applications. The synthetic example was the Society of Petroleum Engineers benchmark Brugge field case. The field example involves waterflooding of a North Sea reservoir with multiple seismic surveys. In both cases, the advantages of incorporating the time-lapse variations were clearly demonstrated through improved estimation of the permeability heterogeneity, fluid saturation evolution, and swept and drained volumes. The value of the seismic data integration was in particular proven through the identification of the continuity in reservoir sands and barriers, and by the preservation of geologic realism in the calibrated model.


2022 ◽  
Vol 12 (1) ◽  
Author(s):  
Manzar Fawad ◽  
Nazmul Haque Mondol

AbstractTo mitigate the global warming crisis, one of the effective ways is to capture CO2 at an emitting source and inject it underground in saline aquifers, depleted oil and gas reservoirs, or in coal beds. This process is known as carbon capture and storage (CCS). With CCS, CO2 is considered a waste product that has to be disposed of properly, like sewage and other pollutants. While and after CO2 injection, monitoring of the CO2 storage site is necessary to observe CO2 plume movement and detect potential leakage. For CO2 monitoring, various physical property changes are employed to delineate the plume area and migration pathways with their pros and cons. We introduce a new rock physics model to facilitate the time-lapse estimation of CO2 saturation and possible pressure changes within a CO2 storage reservoir based on physical properties obtained from the prestack seismic inversion. We demonstrate that the CO2 plume delineation, saturation, and pressure changes estimations are possible using a combination of Acoustic Impedance (AI) and P- to S-wave velocity ratio (Vp/Vs) inverted from time-lapse or four-dimensional (4D) seismic. We assumed a scenario over a period of 40 years comprising an initial 25 year injection period. Our results show that monitoring the CO2 plume in terms of extent and saturation can be carried out using our rock physics-derived method. The suggested method, without going into the elastic moduli level, handles the elastic property cubes, which are commonly obtained from the prestack seismic inversion. Pressure changes quantification is also possible within un-cemented sands; however, the stress/cementation coefficient in our proposed model needs further study to relate that with effective stress in various types of sandstones. The three-dimensional (3D) seismic usually covers the area from the reservoir's base to the surface making it possible to detect the CO2 plume's lateral and vertical migration. However, the comparatively low resolution of seismic, the inversion uncertainties, lateral mineral, and shale property variations are some limitations, which warrant consideration. This method can also be applied for the exploration and monitoring of hydrocarbon production.


2016 ◽  
Vol 4 (4) ◽  
pp. T543-T556
Author(s):  
Sait Baytok ◽  
Şeref Arzu Aktepe ◽  
Muhlis Ünaldi

The Thrace Basin that is located in northwestern Turkey contains sandstone and carbonate reservoirs of Eocene and Oligocene age. Production and exploration activities are still underway. Mapping undrained sweet spots from seismic data is currently a challenge, so time lapse (4D) seismic is used to reduce the risk for new production and development drilling. We have evaluated the normalization and amplitude variation with offset (AVO) analysis of 3D-4D land seismic data in a gas producing field from which baseline and monitor surveys were acquired in 2002 and 2011, respectively. Through AVO analysis, intercept (A) and gradient (B) analysis was conducted, and fluid factor (FF) attribute maps were generated for the assessment of the remaining potential areas. Synthetic gathers were created for simulation of the AVO response, drained and undrained stages and compared with the corresponding 4D seismic data. The drainage of gas from the reservoir interval is evident from the difference maps between 2002 and 2011 seismic data. Both data sets were processed using an amplitude friendly processing sequence. This parallel processing was followed a mild data conditioning and crossequalization for reliable 4D interpretation. The 4D seismic data, especially land data, has low repeatability and requires conditioning to reduce the 4D noise. The 4D noise can be described as nonrepeatable noise, and any difference outside the reservoir zone is not related to production. A so-called crossequalization was applied to the base and the monitor data to bring out similarities so that they cancel out when differences of seismic data and its attributes indicated only the production results over the reservoir zones. As the available 4D data crossequalization software was implemented for stack data only, we created angle band stacks and crossequalized each angle band stack from the base and the monitor data cubes. Five angle band stacks from the base and the monitor prestack data cubes 0°–55° (0°–15°, 15°–25°, 25°–35°, 35°–45°, and 45°–55°) were crossequalized individually. The crossequalized angle band stacks were used in AVO analysis and AVO inversion to generate pore fill identifiers such as FF to map possible undrained zones after 10 years of production.


Sign in / Sign up

Export Citation Format

Share Document