scholarly journals Imaging crustal structures through a passive seismic imaging approach in a mining area in Saxony, Germany

Solid Earth ◽  
2021 ◽  
Vol 12 (12) ◽  
pp. 2703-2715
Author(s):  
Hossein Hassani ◽  
Felix Hloušek ◽  
Stefan Buske ◽  
Olaf Wallner

Abstract. We have used several flooding-induced microseismic events that occurred in an abandoned mining area to image geological structures close to the hypocentres in the vicinity of the mine. The events have been located using a migration-based localization approach. We used the recorded full waveforms of these localized microseismic events and have processed these passive source data as if they resulted from active sources at the known hypocentre location and origin time defined by the applied location approach. The imaging was then performed using a focusing 3D prestack depth migration approach for the secondary P-wave arrivals. The needed 3D migration velocity model was taken from a recent 3D active (controlled-source) seismic survey in that area. We observed several clear and pronounced reflectors in our obtained 3D seismic image cube, some of them related to a major fault zone in that area and some correlating well with information from the nearby mining activities. We compared our results to the 3D seismic image cube obtained directly from the 3D active seismic survey and have found new structures with our approach that were not known yet, probably because of their steep dips which the 3D active seismic survey had not illuminated. The location of the hypocentres at depth with respect to the illumination angles of those structures proved to be favourable in that case, and our 3D passive image complements the 3D active seismic image in an elegant way, thereby revealing new structures that cannot be imaged otherwise with surface seismic configurations alone.

2021 ◽  
Author(s):  
Hossein Hassani ◽  
Felix Hloušek ◽  
Stefan Buske ◽  
Olaf Wallner

Abstract. We have used several flooding induced microseismic events that occurred in an abandoned mining area to image geological structures close to the hypocentres in the vicinity of the mine. The events have been located using a migration-based localization approach. We used the recorded full waveforms of these localized microseismic events and have processed these passive source data as if they resulted from active sources at the known hypocentre location and origin time defined by the applied location approach. The imaging was then performed by using a focusing 3D prestack depth migration approach for the secondary P-wave arrivals. The needed 3D migration velocity model was taken from a recent 3D active (controlled-source) seismic survey in that area. We observed several clear and pronounced reflectors in our obtained 3D seismic image cube, some of them related to a major fault zone in that area and some correlating well with information from the nearby mining activities. We compared our results to the 3D seismic image cube obtained directly from the 3D active seismic survey and have found new structures with our approach that were not know yet, probably because of their steep dips which the 3D active seismic survey had not illuminated. The location of the hypocentres at depth with respect to the illumination angles of those structures proved to be favourable in that case, and our 3D passive image complements the 3D active seismic image in an elegant way thereby revealing new structures that cannot be imaged otherwise with surface seismic configurations alone.


2021 ◽  
Author(s):  
Myriam Lajaunie ◽  
Céleste Broucke ◽  
Jean-Philippe Malet ◽  
Clément Hibert ◽  
Guy Sénéchal ◽  
...  

<div> <p>Bedrock geometry, geological discontinuities, geotechnical units and shear surfaces/bands control the deformation patterns and the mechanisms of slope instabilities. Seismic P-wave refraction tomography is useful to detect these features because P-wave velocity significantly decreases in fractured and weathered rocks relative to consolidated ones, and because lateral changes of velocity can highlight alternation of dipping fracture zones and consolidated rocks. Acquiring this information at high spatial resolution is of paramount importance to model landslide behaviour. </p> </div><div> <p>The Viella slope instability (Hautes-Pyrénées, France) is a complex and deep-seated (> 80 m) landslide which has reactivated in Spring 2018 as a consequence of both a 100-yr return period flash flood (Bastan torrent) which affects the lower part of the slope, and a major rockslide (> 100.000 m<sup>3</sup>) modifying the stress conditions in the upper part. The landslide, which covers an area of ca. 650 000 m², is primarily composed of schists with different degrees of weathering, forming several kinematic units with surface velocities in the range [0.5 – 5] mm.month<sup>-1</sup>. Many buildings and infrastructures (roads, bridge) are progressively damaged (cracks, progressive tilting) and scarps and lobes develop at the surface delineating the kinematic units.  </p> </div><div> <p>In order to model the evolution of the landslide and design possible mitigation measures (drainage, slope reprofiling), a 3D seismic survey has been carried out in summer 2020. The survey was designed to provide a highly detailed velocity model untill 100 m depth, highlighting possible lithological and mechanical contrasts as well as water preferential flow paths. The acquisition was carried out using 71 3C miniaturized seismic sensors buried at ca. 30 cm in the ground and spaced with an average intertrace of 70 m in accordance with slope topography. IGU16HR-3C 5Hz SmartSolo geophones of the DENSAR service (EOST) were used. The seismic array was recording continuously from June, 22nd to July, 21st 2020 at a sampling rate of 500 Hz. 370 controlled seismic sources were triggered at 122 locations using blank 12-gauge shotgun cartridges, Seismic Impulse Source Systemshots, 90-kg Propelled Energy Generator shots and a Mechatronics Lightning source generating P and S-waves with mono-frequency and sweep signals between 5 and 60 Hz of maximum 80 s length.  </p> </div><div> <p>We present the results of this active P-wave traveltime tomography. We first discuss the quality of the recorded signals related to each different type of source, given the noise and attenuation conditions at Viella. Because the signals were challenging to detect a methodology based on manual picking was used, supported by automatic detection tools and considerations regarding the network geometry in an a priori velocity model.  </p> </div><div> <p>The P-wave model was obtained using the inversion library pyGIMLI, which permits an accurate description of the topography, and provides a spatial discretization adapted to the problem. To supplement and constrain the interpretation of the P-wave velocity model, borehole information as well as a 3D resistivity model of the zone are available. With regards to these data, the geometric features and physical parameters of the main geological structures of the landslide are discussed. </p> </div>


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B41-B57 ◽  
Author(s):  
Himanshu Barthwal ◽  
Mirko van der Baan

Microseismicity is recorded during an underground mine development by a network of seven boreholes. After an initial preprocessing, 488 events are identified with a minimum of 12 P-wave arrival-time picks per event. We have developed a three-step approach for P-wave passive seismic tomography: (1) a probabilistic grid search algorithm for locating the events, (2) joint inversion for a 1D velocity model and event locations using absolute arrival times, and (3) double-difference tomography using reliable differential arrival times obtained from waveform crosscorrelation. The originally diffusive microseismic-event cloud tightens after tomography between depths of 0.45 and 0.5 km toward the center of the tunnel network. The geometry of the event clusters suggests occurrence on a planar geologic fault. The best-fitting plane has a strike of 164.7° north and dip angle of 55.0° toward the west. The study region has known faults striking in the north-northwest–south-southeast direction with a dip angle of 60°, but the relocated event clusters do not fall along any mapped fault. Based on the cluster geometry and the waveform similarity, we hypothesize that the microseismic events occur due to slips along an unmapped fault facilitated by the mining activity. The 3D velocity model we obtained from double-difference tomography indicates lateral velocity contrasts between depths of 0.4 and 0.5 km. We interpret the lateral velocity contrasts in terms of the altered rock types due to ore deposition. The known geotechnical zones in the mine indicate a good correlation with the inverted velocities. Thus, we conclude that passive seismic tomography using microseismic data could provide information beyond the excavation damaged zones and can act as an effective tool to complement geotechnical evaluations.


Geophysics ◽  
2019 ◽  
Vol 84 (3) ◽  
pp. KS59-KS69 ◽  
Author(s):  
Chao Song ◽  
Zedong Wu ◽  
Tariq Alkhalifah

Passive seismic monitoring has become an effective method to understand underground processes. Time-reversal-based methods are often used to locate passive seismic events directly. However, these kinds of methods are strongly dependent on the accuracy of the velocity model. Full-waveform inversion (FWI) has been used on passive seismic data to invert the velocity model and source image, simultaneously. However, waveform inversion of passive seismic data uses mainly the transmission energy, which results in poor illumination and low resolution. We developed a waveform inversion using multiscattered energy for passive seismic to extract more information from the data than conventional FWI. Using transmission wavepath information from single- and double-scattering, computed from a predicted scatterer field acting as secondary sources, our method provides better illumination of the velocity model than conventional FWI. Using a new objective function, we optimized the source image and velocity model, including multiscattered energy, simultaneously. Because we conducted our method in the frequency domain with a complex source function including spatial and wavelet information, we mitigate the uncertainties of the source wavelet and source origin time. Inversion results from the Marmousi model indicate that by taking advantage of multiscattered energy and starting from a reasonably acceptable frequency (a single source at 3 Hz and multiple sources at 5 Hz), our method yields better inverted velocity models and source images compared with conventional FWI.


Author(s):  
D. Hollis ◽  
C. Cox ◽  
R. Clayton ◽  
F. Lin ◽  
D. Li ◽  
...  

2020 ◽  
Author(s):  
Vera Lay ◽  
Stefan Buske ◽  
Sascha Barbara Bodenburg ◽  
Franz Kleine ◽  
John Townend ◽  
...  

<p>The Alpine Fault along the West Coast of the South Island (New Zealand) is a major plate boundary that is expected to rupture in the next 50 years, likely as a magnitude 8 earthquake. The Deep Fault Drilling Project (DFDP) aims to deliver insight into the geological structure of this fault zone and its evolution by drilling and sampling the Alpine Fault at depth.  </p><p>Here we present results from a 3D seismic survey around the DFDP-2 drill site in the Whataroa Valley where the drillhole penetrated almost down to the fault surface. Within the glacial valley, we collected 3D seismic data to constrain valley structures that were obscured in previous 2D seismic data. The new data consist of a 3D extended vertical seismic profiling (VSP) survey using three-component receivers and a fibre optic cable in the DFDP-2B borehole as well as a variety of receivers at the surface.</p><p>The data set enables us to derive a reliable 3D P-wave velocity model by first-arrival travel time tomography. We identify a 100-460 m thick sediment layer (average velocity 2200±400 m/s) above the basement (average velocity 4200±500 m/s). Particularly on the western valley side, a region of high velocities steeply rises to the surface and mimics the topography. We interpret this to be the infilled flank of the glacial valley that has been eroded into the basement. In general, the 3D structures implied by the velocity model on the upthrown (Pacific Plate) side of the Alpine Fault correlate well with the surface topography and borehole findings.</p><p>A reliable velocity model is not only valuable by itself but it is also required as input for prestack depth migration (PSDM). We performed PSDM with a part of the 3D data set to derive a structural image of the subsurface within the Whataroa Valley. The top of the basement identified in the P-wave velocity model coincides well with reflectors in the migrated images so that we can analyse the geometry of the basement in detail.</p>


2014 ◽  
Vol 2 (2) ◽  
pp. T79-T88 ◽  
Author(s):  
Adam D. Halpert ◽  
Robert G. Clapp ◽  
Biondo Biondi

Although it is a crucial component of seismic velocity model building, salt delineation is often a major bottleneck in the interpretation workflow. Automatic methods like image segmentation can help to alleviate this bottleneck, but issues with accuracy and efficiency can hinder their effectiveness. However, a new graph-based segmentation algorithm can, after modifications to account for the unique nature of seismic data, quickly and accurately delineate salt bodies on 3D seismic images. In areas where salt boundaries are poorly imaged, limited manual interpretations can be used to guide the automatic segmentation, allowing for interpreter insight to be combined with modern computational capabilities. A successful 3D field data example demonstrates that this method could become an important tool for interactive interpretation tasks.


Geophysics ◽  
1995 ◽  
Vol 60 (4) ◽  
pp. 1178-1186 ◽  
Author(s):  
M. Reza Daneshvar ◽  
Clarence S. Clay ◽  
Martha K. Savage

We have developed a method of processing seismic signals generated by microearthquakes to image local subsurface structure beneath a recording station. This technique uses the autocorrelation of the vertically traveling earthquake signals to generate pseudoreflection seismograms that can be interpreted for subsurface structure. Processed pseudoreflection data, from microearthquakes recorded in the island of Hawaii, show consistent reflectivity patterns that are interpreted as near‐surface horizontal features. Forward modeling of the pseudoreflection data results in a P‐wave velocity model that shows reasonable agreement with the velocity model derived from a refraction study in the region. Usable signal‐to‐noise ratio is obtained down to 2 s. A shear‐wave velocity model was also generated by applying this technique to horizontal component data.


2014 ◽  
Vol 6 (1) ◽  
pp. 985-1021
Author(s):  
I. Janutyte ◽  
E. Kozlovskaya ◽  
M. Majdanski ◽  
P. H. Voss ◽  
M. Budraitis ◽  
...  

Abstract. The presented study is a part of the passive seismic experiment PASSEQ 2006–2008 which took place around the Trans-European Suture Zone (TESZ) from May 2006 to June 2008. The dataset of 4195 manually picked arrivals of teleseismic P waves of 101 earthquakes (EQs) recorded in the PASSEQ seismic stations deployed to the east of the TESZ was inverted using the non-linear teleseismic tomography algorithm TELINV. Two 3-D crustal models were used to estimate the crustal travel time (TT) corrections. As a result, we obtained a model of P wave velocity variations in the upper mantle beneath the TESZ and the EEC. In the study area beneath the craton we observed 5 to 6.5% higher and beneath the TESZ about 4% lower seismic velocities compared to the IASP91 velocity model. We found the seismic lithosphere-asthenosphere boundary (LAB) beneath the TESZ at a depth of about 180 km, while we observed no seismic LAB beneath the EEC. The inversion results obtained with the real and the synthetic datasets indicated a ramp shape of the LAB in the northern TESZ where we observed values of seismic velocities close to those of the craton down to about 150 km. The lithosphere thickness in the EEC increases going from the TESZ to the NE from about 180 km beneath Poland to 300 km or more beneath Lithuania. Moreover, in western Lithuania we possibly found an upper mantle dome. In our results the crustal units are not well resolved. There are no clear indications of the features in the upper mantle which could be related with the crustal units in the study area. On the other hand, at a depth of 120–150 km we possibly found a trace of a boundary of proposed palaeosubduction zone between the East Lithuanian Domain (EL) and the West Lithuanian Granulite Domain (WLG). Also, in our results we may have identified two anorogenic granitoid plutons.


2021 ◽  
Author(s):  
Vera Lay ◽  
Stefan Buske ◽  
Franz Kleine ◽  
John Townend ◽  
Richard Kellett ◽  
...  

<p>The Alpine Fault at the West Coast of the South Island (New Zealand) is a major plate boundary that is expected to rupture in the next 50 years, likely as a magnitude 8 earthquake. The Deep Fault Drilling Project (DFDP) aimed to deliver insight into the geological structure of this fault zone and its evolution by drilling and sampling the Alpine Fault at depth. Here we present results from a seismic survey around the DFDP-2 drill site in the Whataroa Valley where the drillhole almost reached the fault plane. This unique 3D seismic survey includes several 2D lines and a 3D array at the surface as well as borehole recordings. Within the borehole, the unique option to compare two measurement systems is used: conventional three-component borehole geophones and a fibre optic cable (heterodyne Distributed Vibration Sensing system (hDVS)). Both systems show coherent signals but only the hDVS system allowed a recording along the complete length of the borehole.</p><p>Despite the challenging conditions for seismic imaging within a glacial valley filled with sediments and steeply dipping valley flanks, several structures related to the valley itself as well as the tectonic fault system are imaged. The pre-processing of the seismic data also includes wavefield separation for the zero-offset borehole data. Seismic images are obtained by prestack depth migration approaches.</p><p>Within the glacial valley, particularly steep valley flanks are imaged directly and correlate well with results from the P-wave velocity model obtained by first arrival travel-time tomography. Additionally, a glacially over-deepened trough with nearly horizontally layered sediments is identified about 0.5 km south of the DFDP-2B borehole.</p><p>With regard to the expected Alpine fault zone, a set of several reflectors dipping 40-56° to the southeast are identified in a ~600 m wide zone between depths of 0.2 and 1.2 km that is interpreted to be the minimum extent of the damage zone. Different approaches image one distinct reflector dipping at 40°, which is interpreted to be the main Alpine Fault reflector. This reflector is only ~100 m ahead from the lower end of the borehole. At shallower depths (z<0.5 km), additional reflectors are identified as fault segments and generally have steeper dips up to 56°. About 1 km south of the drill site, a major fault is identified at a depth of 0.1-0.5 km that might be caused by the regional tectonics interacting with local valley structures. A good correlation is observed among the separate seismic data sets and with geological results such as the borehole stratigraphy and the expected surface trace of the fault.</p><p>In conclusion, several structural details of the fault zone and its environment are seismically imaged and show the complexity of the Alpine Fault at the Whataroa Valley. Thus, a detailed seismic characterization clarifies the subsurface structures, which is crucial to understand the transpressive fault’s tectonic processes.</p>


Sign in / Sign up

Export Citation Format

Share Document