1. |
- Abreu, Rafael, et al.
(författare)
-
Micropolar modelling of rotational waves in seismology
- 2017
-
Ingår i: Geophysical Journal International. - : OXFORD UNIV PRESS. - 0956-540X .- 1365-246X. ; 210:2, s. 1021-1046
-
Tidskriftsartikel (refereegranskat)abstract
- In this contribution we study elastic wave propagation via the introduction of the micropolar theory. As a generalization of a classical linear elastic medium, a micropolar medium allows each particle to have intrinsic rotational degrees of freedom (spin). We perform numerical experiments using the Pseudospectral method. We find analytical harmonic micropolar solutions for different problem configurations, which result in waveform differences between the classical linear elastic and micropolar media. In contrast to linear elastic media, wave propagation in micropolar media is dispersive. We study how the spin waveform depends on the micropolar elastic parameters and frequency content of the simulation. The micropolar effect on numerical seismograms has a direct implication on the phase, amplitude and arrival time. For frequencies lower than the cut-off frequency, the spin waveform has the same amplitude as the macrorotation field. For frequencies higher than the cut-off frequency, the amplitude of the spin waveform decreases with increasing frequency, so that then it is no longer comparable to the amplitude of macroscopic rotations. When both frequencies are equal there is no wave propagation. This work attempts to clarify the theory of micropolar media for its applications in seismology. We argue that micropolar theory should be further investigated for its potential uses in seismology to, for example, describe energy dissipation, seismograms recorded with rotational seismometers and rupture processes.
|
|
2. |
- Abtahi, Sayyed Mohammad, et al.
(författare)
-
A new reference model for 3D inversion of airborne magnetic data in hilly terrain : A case study from northern Sweden
- 2018
-
Ingår i: Geophysics. - : SOC EXPLORATION GEOPHYSICISTS. - 0016-8033 .- 1942-2156. ; 83:1, s. B1-B12
-
Tidskriftsartikel (refereegranskat)abstract
- The inherent nonuniqueness in modeling magnetic data can be partly reduced by adding prior information, either as mathematical constructs or simply as bounds on magnetization obtained from laboratory measurements. If a good prior model can be used as a reference model, then the quality of estimated models through an inverse approach can be greatly improved. But even though data on magnetic properties of rocks might exist, their distribution may often be quite irregular on local and regional scales, so that it is difficult to define representative classes of rock types suitable for constraining geophysical models of magnetization. We have developed a new way of constructing a reference model that varies only laterally and is confined to the part of the terrain that lies above the lowest topography in the area. To obtain this model, several estimated 2D magnetization distributions were constructed by data inversion as a function of the iteration number. Then, a suitable 2D model of the magnetization in the topography was chosen as a starting point for constructing a 3D reference model by modifying it with a vertical decay such that its average source depth was the same for all horizontal positions. The average source depth of the reference model was chosen to satisfy the average source depth obtained from analyzing the radial power spectrum of the area studied. Finally, the measured magnetic data were inverted in three dimensions using the given reference model. For a selected reference model, shallow structures indicated a better overall correlation with large remanent magnetizations measured on rock samples from the area. Throughout the entire model, the direction of magnetization was allowed to vary freely. We found that the Euclidean norm of the estimated model was reduced compared with the case where the magnetization direction was fixed.
|
|
3. |
|
|
4. |
- Abtahi, Sayyed Mohammad
(författare)
-
Airborne Gravity Gradient, Magnetic and VLF datasets : Case studies of modelling, inversion and interpretation
- 2016
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Northern Sweden is one of the largest hosts for mineral resources in Europe and always has been an interesting area for researchers from various disciplines of Earth sciences. This dissertation is a comprehensive summary of three case study papers on airborne VLF, gravity gradient and magnetic data in the area.In the first paper, tensor VLF data is extracted from an old data set which contains only the total and the vertical magnetic components. The anomalous part of the horizontal magnetic field components is computed by a Hilbert transform of the vertical magnetic field. The normal part of the horizontal magnetic field component is computed as a function of total, vertical and anomalous part of horizontal magnetic fields. The electric field is also calculated for TE mode and impedance tensor and apparent resistivity are computed. In addition tippers are calculated for two transmitters and inverted by a 3D inversion algorithm. Comparison of the estimated model and geology map of bedrock shows that lower resistivity zones are correlated with mineralizations.The second paper deals with the internal consistency of airborne gravity gradient data. The six components of the data are estimated from a common potential function. It is shown that the data is adequately consistent but at shorter land clearances the difference between the estimated data and the original data is larger. The technique is also used for computing the Bouguer anomaly from terrain corrected FTG data. Finally the data is inverted in 3D, which shows that the estimated density model in shallow depth is dominated by short wave length features.Inversion of TMI data is the topic of the third paper where a new type of reference model for 3D inversion of magnetic data is proposed by vertically extending the estimated magnetization of a 2D terrain magnetization model. The final estimated 3D result is compared with the magnetization model where no reference model is used. The comparison shows that using the reference model helps the high magnetization zones in the estimated model at shallow depths to be better correlated with measured high remanent magnetization from rock samples. The high magnetization zones are also correlated with gabbros and volcanic metasediments.
|
|
5. |
- Abtahi, Sayyed Mohammad, et al.
(författare)
-
Consistency investigation, vertical gravity estimation and inversion of airborne gravity gradient data – A case study from northern Sweden
- 2016
-
Ingår i: Geophysics. - : Society of Exploration Geophysicists. - 0016-8033 .- 1942-2156. ; 81:3, s. B65-B76
-
Tidskriftsartikel (refereegranskat)abstract
- For airborne gravity gradient data, it is a challenge to distinguish between high-frequency intrinsic and dynamically produced noise caused by the aircraft and small-scale effects from shallow density variations. To facilitate consistent interpretation, techniques that include all of the measured gravity gradient components are particularly promising. We represented the measurements by a common potential function accounting for lateral and height variations. Thus, it was possible to evaluate the internal consistency between the measured components and to identify components with bias or particularly strong noise. As an extra benefit for data sets that contain terrain-corrected and nonterrain-corrected gravity gradient measurements at flight altitude, we estimated terrain-corrected anomalies on the topographic relief using downward continuation and retrieved nonterrain-corrected gravity gradient data suitable for inversion using upward continuation. For a field data set from northern Sweden, the largest differences (up to 50 eotvos) between the measured and estimated components of the gravity gradient data were found in areas of high topographical relief. But the average residual standard deviations of the individual components were between 3.6 and 7.4 eotvos, indicating that the components were consistent in an average sense. We have determined the successful conversion of terrain-corrected airborne gravity gradient data to Bouguer gravity data on the topographic relief using ground-based vertical gravity data as a reference. A 3D inverse model computed from the nonterrain-corrected data clearly showed the depth extent of the geologic structures observed at the surface, but it only produced a weak representation of the shallow structure. In contrast, a 2D surface density model in which only lateral variations of density in the topographic relief was allowed exhibited more realistic density distributions in fair correlation with geology.
|
|
6. |
- Abtahi, Sayyed Mohammad, et al.
(författare)
-
Extracting geoelectrical maps from vintage very-low-frequency airborne data, tipper inversion, and interpretation : A case study from northern Sweden
- 2016
-
Ingår i: Geophysics. - 0016-8033 .- 1942-2156. ; 81:5, s. B135-B147
-
Tidskriftsartikel (refereegranskat)abstract
- In 1985, the mining company Luossavaara-Kiirunavaara Aktiebolag collected airborne very-low-frequency (VLF) data in northern Sweden. The operators stored only the vertical component and the total magnetic field, which at that time were believed to be sufficient for qualitative interpretation. Therefore, the data could not be directly used for quantitative tensor VLF processing and inversion. To avoid the costs of resurveying, we have developed a novel technique to estimate the tippers from the measured VLF data by computing anomalous and normal parts of the horizontal components of the magnetic field from two transmitters separately. Retrieval of the normal horizontal components was possible because one component of the horizontal magnetic field was used as the phase reference during the measurements. Additionally, we have determined how the approximate apparent resistivity suitable for data visualization can be computed from the components of the magnetic field assuming an average normal resistivity of the subsurface. Maps of apparent resistivity combined with topography show a clear correlation between high topography and high resistivity, whereas conductive zones are found in valleys in between. More importantly, the 3D model inverted from the calculated tippers shows excellent agreement with a map of the surface geology. Based on this comparison, some less resistive zones can be related to fluids in fractures and others can be related to mineralized contact zones. We suggest to focus further exploration on conductive zones surrounding areas with basaltic composition.
|
|
7. |
- Braun, Martina, et al.
(författare)
-
Simultaneous inversion of magnetic resonance sounding in terms of water content, resistivity and decay times
- 2009
-
Ingår i: Near Surface Geophysics. - 1569-4445 .- 1873-0604. ; 7:5-6, s. 589-598
-
Tidskriftsartikel (refereegranskat)abstract
- Magnetic resonance sounding (MRS) or surface nuclear magnetic resonance (SNMR) is used for direct groundwater exploration and for an improved aquifer characterization. Currently, it is the only geophysical method that is capable of directly determining the free water content and estimating the pore sizes of the aquifer in the subsurface. However, MRS is basically an electromagnetic method. Therefore, it is sensitive to the resistivity of the subsurface. The water content is the main target of investigation, therefore first inversion routines focused on the water content. Later on, inversion routines determining water content and decay times became available. Very recently, MRS inversion for water content and resistivity has been realized. We present here a simultaneous inversion of MRS in terms of determining the three inversion parameters - water content, resistivity and decay time - within one single inversion routine. Within the iterative inversion scheme, the extrapolated initial values are determined on the basis of the physical effective decay times in the subsurface, which are estimated within the inversion scheme. Due to an instrumental dead time, the initial values for amplitude and phase, which are related to water content and resistivity, cannot be measured directly. Therefore, the initial amplitude must be extrapolated using the decay time of the signal. The standard approach is a mono-exponential decay curve; implicitly, the phase is assumed to be time-invariant. However, multi-exponential signals are natural when considering relaxation behaviour in the underground. It originates from multi-modal pore size distributions or simply a number of differently relaxing signal contributions from the various lithological units.
|
|
8. |
- Cherevatova, M., et al.
(författare)
-
Crustal structure beneath southern Norway imaged by magnetotellurics
- 2014
-
Ingår i: Tectonophysics. - : Elsevier BV. - 0040-1951 .- 1879-3266. ; 628, s. 55-70
-
Tidskriftsartikel (refereegranskat)abstract
- We use data from two magnetotelluric profiles, ToSca10 and ToSca'09, over the Scandinavian Mountains to study the crustal structure in southern Norway. The profiles cross the major tectonic structures of the Caledonian orogen as well as the western margin of the Precambrian Baltica. Dimensionality and strike analyses indicate generally 3-D behavior of the data. However, the majority of the used data distinguishes a preferable strike direction, which is supported by the geology of the region. Hence, we employ 2-D inversion and choose to invert the determinant of the impedance tensor to mitigate 3-D effects in the data on our 2-D models. Magnetotelluric data from both profiles are inverted using a damped least squares solution based on a singular value decomposition. We improved the solution by defining the inverse model covariance matrix through gradient or Laplacian smoothing operators. The two-dimensional inversion models of the ToSca'09 and ToSca'10 field data from southern Norway derived from the damped least squares scheme with the Laplacian inverse model covariance matrix are presented. Resistive rocks, extending to the surface, image the autochthonous Southwest Scandinavian Domain and the allochthonous Western Gneiss Region. Near-surface conductors, which are located between the resistive Caledonian nappes and Precambrian basement, delineate highly conductive shallow-sea sediments, so called alum shales. They exhibit a decollement along which the Caledonian nappes were overthrust. A deeper, upper to midcrustal conducting layer in the Southwest Scandinavian Domain may depict the remnants of closed ocean basins formed during the accretions and collisions of various Sveconorwegian terranes. In ToSca'10, the Caledonian nappes, the conducting alum shales and the deeper conductor are terminated in the west by the Faltungsgraben shear complex which represents a crustal scale boundary between the Western Gneiss Region in the west and the Southwest Scandinavian Domain in the east.
|
|
9. |
- Cherevatova, M., et al.
(författare)
-
Magnetotelluric array data analysis from north-west Fennoscandia
- 2015
-
Ingår i: Tectonophysics. - : Elsevier BV. - 0040-1951 .- 1879-3266. ; 653, s. 1-19
-
Tidskriftsartikel (refereegranskat)abstract
- New magnetotelluric (MT) data in north-west Fennoscandia were acquired within the framework of the project "Magnetotellurics in the Scandes" (MaSca). The project focuses on the investigation of the crustal and upper mantle lithospheric structure in the transition zone from stable Precambrian cratonic interior to passive continental margin beneath the Caledonian orogen and the Scandinavian Mountains in western Fennoscandia. An array of 59 synchronous long period and 220 broad-band MT sites was occupied in the summers of 2011 to 2013. We estimated MT transfer functions in the period range from 0.003 to 10(5) s. The Q-function multi-site multi-frequency analysis and the phase tensor were used to estimate strike and dimensionality of MT data. Dimensionality and strike analyses indicate generally 2-D behaviour of the data with 3-D effects at some sites and period bands. In this paper we present 2-D inversion of the data, 3-D inversion models are shown in the parallel paper. We choose to invert the determinant of the impedance tensor to mitigate 3-D effects in the data on our 2-D models. Seven crustal-scale and four lithospheric-scale 2-D models are presented. The resistive regions are images of the Archaean and Proterozoic basement in the east and thin Caledonian nappes in the west. The middle and lower crust of the Svecofennian province is conductive. The southern end of the Kittila Greenstone Belt is seen in the models as a strong upper to middle crustal conductor. In the Caledonides, the highly conductive alum shales are observed along the Caledonian Thrust Front. The thickest lithosphere is in the Palaeoproterozioc Svecofennian Domain, not in the Archaean. The thickness of the lithosphere is around 200 km in the north and 300 km in the south-west.
|
|
10. |
- Juhojuntti, Niklas, et al.
(författare)
-
Joint inversion of seismic refraction and resistivity data using layered models : Applications to groundwater investigation
- 2015
-
Ingår i: Geophysics. - 0016-8033 .- 1942-2156. ; 80:1, s. EN43-EN55
-
Tidskriftsartikel (refereegranskat)abstract
- We developed a method for joint inversion of seismic refraction and resistivity data, using sharp-boundary models with few layers (typically three). We demonstrated the usefulness of the approach via examples from near-surface case studies involving shallow groundwater exploration and geotechnical investigations, although it should also be applicable to other types of layered environments, e.g., sedimentary basins. In our model parameterization, the layer boundaries were common for the resistivity and velocity distributions. Within the layers, only lateral variations in the material parameters (resistivity and velocity) were allowed, and we assumed no correlation between these. The inversion was performed using a nonlinear least-squares algorithm, using lateral smoothing to the layer boundaries and to the materialparameters. Depending on the subsurface conditions, the smoothing can be applied either to the depth of the layer boundaries or to the layer thicknesses. The forward responses and Jacobian for refraction seismics were calculated through ray tracing. The resistivity computations were performed with finite differences and a cell-to-layer transform for the Fréchet derivatives. Our method performed well in synthetic tests, and in the case studies, the layer boundaries were in good agreement with in situ tests and seismic reflection data, although minimum-structure inversion generally has a better data fit due to more freedom to introduce model heterogeneity. We further found that our joint inversion approach can provide more accurate thickness estimates for seismic hidden layers.
|
|