Imaging of reflection seismic energy for mapping shallow fracture zones in crystalline rocks

Geophysics ◽  
1994 ◽  
Vol 59 (5) ◽  
pp. 753-765 ◽  
Author(s):  
J. S. Kim ◽  
Wooil M. Moon ◽  
Ganpat Lodha ◽  
Mulu Serzu ◽  
Nash Soonawala

The high‐resolution reflection seismic technique is being used increasingly to address geologic exploration and engineering problems. There are, however, a number of problems in applying reflection seismic techniques in a crystalline rock environment. The reflection seismic data collected over a fractured crystalline rock environment are often characterized by low signal‐to‐noise ratios (S/N) and inconsistent reflection events. Thus it is important to develop data processing strategies and correlation schemes for the imaging of fracture zones in crystalline rocks. Two sets of very low S/N, high‐resolution seismic data, previously collected by two different contractors in Pinawa, Canada, and the island of Äspö, Sweden, were reprocessed and analyzed, with special emphasis on the shallow reflection events occurring at depths as shallow as 60–100 m. The processing strategy included enhancing the signals hidden behind large‐amplitude noise, including clipped ground roll. The pre‐ and poststack processing includes shot f-k filtering, residual statics, careful muting after NMO correction, energy balance, and coherency filtering. The final processed seismic sections indicate that reflected energy in these data sets is closely related to rock quality in Äspö data and fracturing in Atomic Energy of Canada, Ltd. (AECL) data. The lithologic boundaries are not clearly mappable in these data. When thickness of the reflection zone is of the order of a wavelength, the top and bottom of the zone may be resolved. The major fracture zones in crystalline rocks correlate closely with the well‐log data and are usually characterized by very low velocity and produce low‐acoustic‐impedance contrasts compared to those of surrounding rocks. Because the incidence angles vary rapidly for shallow‐reflection geometries, segments of major fracture zones can effectively be analyzed in terms of reflectivity. Reflection images of each fracture zone were investigated in the common‐offset section, where each focused event was associated with a consistent incidence angle on the reflectivity map. The complex attributes of the data indicate that strong reflectors at shallow depth coincide with intensely fractured zones. These correlate well with instantaneous amplitude plots and instantaneous frequency plots. The instantaneous phase plot also identifies the major and minor fractures.

Geophysics ◽  
2012 ◽  
Vol 77 (1) ◽  
pp. A5-A8 ◽  
Author(s):  
David Bonar ◽  
Mauricio Sacchi

The nonlocal means algorithm is a noise attenuation filter that was originally developed for the purposes of image denoising. This algorithm denoises each sample or pixel within an image by utilizing other similar samples or pixels regardless of their spatial proximity, making the process nonlocal. Such a technique places no assumptions on the data except that structures within the data contain a degree of redundancy. Because this is generally true for reflection seismic data, we propose to adopt the nonlocal means algorithm to attenuate random noise in seismic data. Tests with synthetic and real data sets demonstrate that the nonlocal means algorithm does not smear seismic energy across sharp discontinuities or curved events when compared to seismic denoising methods such as f-x deconvolution.


2013 ◽  
Vol 167 ◽  
pp. 72-83 ◽  
Author(s):  
J.S. L'Heureux ◽  
M. Long ◽  
M. Vanneste ◽  
G. Sauvin ◽  
L. Hansen ◽  
...  

2020 ◽  
Author(s):  
Hao Zhang ◽  
Jianguang Han ◽  
Heng Zhang ◽  
Yi Zhang

<p>The seismic waves exhibit various types of attenuation while propagating through the subsurface, which is strongly related to the complexity of the earth. Anelasticity of the subsurface medium, which is quantified by the quality factor Q, causes dissipation of seismic energy. Attenuation distorts the phase of the seismic data and decays the higher frequencies in the data more than lower frequencies. Strong attenuation effect resulting from geology such as gas pocket is a notoriously challenging problem for high resolution imaging because it strongly reduces the amplitude and downgrade the imaging quality of deeper events. To compensate this attenuation effect, first we need to accurately estimate the attenuation model (Q). However, it is challenging to directly derive a laterally and vertically varying attenuation model in depth domain from the surface reflection seismic data. This research paper proposes a method to derive the anomalous Q model corresponding to strong attenuative media from marine reflection seismic data using a deep-learning approach, the convolutional neural network (CNN). We treat Q anomaly detection problem as a semantic segmentation task and train an encoder-decoder CNN (U-Net) to perform a pixel-by-pixel prediction on the seismic section to invert a pixel group belongs to different level of attenuation probability which can help to build up the attenuation model. The proposed method in this paper uses a volume of marine 3D reflection seismic data for network training and validation, which needs only a very small amount of data as the training set due to the feature of U-Net, a specific encoder-decoder CNN architecture in semantic segmentation task. Finally, in order to evaluate the attenuation model result predicted by the proposed method, we validate the predicted heterogeneous Q model using de-absorption pre-stack depth migration (Q-PSDM), a high-resolution depth imaging result with reasonable compensation is obtained.</p>


Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1278-1289 ◽  
Author(s):  
William J. Lutter ◽  
Rufus D. Catchings ◽  
Craig M. Jarchow

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.


Geophysics ◽  
2012 ◽  
Vol 77 (6) ◽  
pp. B269-B285 ◽  
Author(s):  
Mai-Britt Jensen ◽  
Artem Kashubin ◽  
Christopher Juhlin ◽  
Sten-Åke Elming

Potential weakness zones due to mining-related fracture development under the town of Kiruna, Sweden, have been investigated by integration of seismic, gravity, and petrophysical data. Reflection seismic data were acquired along two subparallel 2D profiles within the residential area of the town. The profiles of [Formula: see text], each oriented approximately east–west, nearly perpendicular to the general strike of the local geology, crossed several contact zones between quartz-bearing porphyries, a sequence of interchanging sedimentary rocks (siltstone, sandstone, conglomerate, and agglomerate), and metabasalt. The resulting reflection seismic sections revealed a strong east-dipping reflectivity that is imaged down to approximately 1.5 km. The location and orientation of major features agree well between the profiles and with the surface geology and known contact zones between the different rock types. Our imaging results, supported by traveltime modelling, indicate that the contact zones dip 40°–50° to the east. The deepest and the weakest reflections are associated with a [Formula: see text] dipping structure that is presumably related to the Kiirunavaara iron mineralization. Tomographic inversion of refracted arrivals revealed a more detailed image of the velocity distribution in the upper 100–200 m along the profiles, enabling us to identify near-surface low velocity zones. These could be possible weakness zones developed along the lithological contacts and within the geologic units. The structural image obtained from the seismic data was used to constrain data inversion along a 28 km long east–northeast to west–southwest-oriented gravity profile. The resulting density model indicates that the quartz-bearing porphyry in the hanging wall of the Kiirunavaara mineralization can be separated into two blocks oriented parallel to the ore body. One block has an unexpected low density, which could be an indication of extensive fracturing and deformation.


Sign in / Sign up

Export Citation Format

Share Document