首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 27 毫秒
1.
According to recent estimates, the continental mid-crust contains 35–40 per cent amphibolites. Heating of the crust by an underlying mantle plume, for example beneath continental rifts, high plateaus, and areas of intraplate volcanic activity, releases water. Dehydration of amphibole-bearing rocks at depths of 20–40  km occurs mainly in the temperature range 650–700 °C, and this releases about 0.4  wt per cent of water.
  Seismic tomography studies of the crust in the Kirgyz Tien Shan Range, where the age of the tectonic activity is less than 30  Ma, revealed a low-velocity zone in the mid-crust. The velocity of P waves was 0.4  km  s−1 lower than in normal crust. MT sounding data in the region show the existence of a low-resistivity layer with an average resistivity of about 25  Ω  m at the depth of the low-velocity layer. The spatial correlation of the observed anomalous layers and calculated effect of fluid phase on seismic and electric parameters of rocks suggests the presence of aqueous fluids released by the heating of the mid-crust.  相似文献   

2.
About 50 000 P and S arrival times and 25 000 values of t * recorded at seismic arrays operated in the Central Andes between 20°S and 25°S in the time period from 1994 to 1997 have been used for locating more than 1500 deep and crustal earthquakes and creating 3-D P , S velocity and Qp models. The study volume in the reference model is subdivided into three domains: slab, continental crust and mantle wedge. A starting velocity distribution in each domain is set from a priori information: in the crust it is based on the controlled sources seismic studies; in slab and mantle wedge it is defined using relations between P and S velocities, temperature and composition given by mineral physics. Each iteration of tomographic inversion consists of the following steps: (1) absolute location of sources in 3-D velocity model using P and S arrival times; (2) double-difference relocation of the sources and (3) simultaneous determination of P and S velocity anomalies, P and S station corrections and source parameters by inverting one matrix. Velocity parameters are computed in a mesh with the density of nodes proportional to the ray density with double-sided nodes at the domain boundaries. The next iteration is repeated with the updated velocity model and source parameters obtained at the previous step. Different tests aimed at checking the reliability of the obtained velocity models are presented. In addition, we present the results of inversion for Vp and Vp/Vs parameters, which appear to be practically equivalent to Vp and Vs inversion. A separate inversion for Qp has been performed using the ray paths and source locations in the final velocity model. The resulting Vp , Vs and Qp distributions show complicated, essentially 3-D structure in the lithosphere and asthenosphere. P and S velocities appear to be well correlated, suggesting the important role of variations of composition, temperature, water content and degree of partial melting.  相似文献   

3.
A wide-angle seismic profile across the western peninsulas of SW Ireland was performed. This region corresponds to the northernmost Variscan thrust and fold deformation. The dense set of 13 shots and 109 stations along the 120  km long profile provides a detailed velocity model of the crust.
  The seismic velocity model, obtained by forward and inverse modelling, defines a five-layer crust. A sedimentary layer, 5–8  km thick, is underlain by an upper-crustal layer of variable thickness, with a base generally at a depth of 10–12  km. Two mid-crustal layers are defined, and a lower-crustal layer below 22  km. The Moho lies at a depth of 30–32  km. A low-velocity zone, which coincides with a well-defined gravity low, is observed in the central part of the region and is modelled as a Caledonian granite which intruded upper-crustal basement. The granite may have acted as a buffer to northward-directed Variscan thrusting. The Dingle–Dungarvan Line (DDL) marks a major change in sedimentary and crustal velocity and structure. It lies immediately to the north of the velocity and gravity low, and shows thickness and velocity differences in many of the underlying crustal layers and even in the Moho. This suggests a deep, pre-Variscan control of the structural development of this area. The model is compatible with thin-skinned tectonics, which terminated at the DDL and which incorporated thrusts involving the sedimentary and upper-crustal layers.  相似文献   

4.
While the inversion of electromagnetic data to recover electrical conductivity has received much attention, the inversion of those data to recover magnetic susceptibility has not been fully studied. In this paper we invert frequency-domain electromagnetic (EM) data from a horizontal coplanar system to recover a 1-D distribution of magnetic susceptibility under the assumption that the electrical conductivity is known. The inversion is carried out by dividing the earth into layers of constant susceptibility and minimizing an objective function of the susceptibility subject to fitting the data. An adjoint Green's function solution is used in the calculation of sensitivities, and it is apparent that the sensitivity problem is driven by three sources. One of the sources is the scaled electric field in the layer of interest, and the other two, related to effective magnetic charges, are located at the upper and lower boundaries of the layer. These charges give rise to a frequency-independent term in the sensitivities. Because different frequencies penetrate to different depths in the earth, the EM data contain inherent information about the depth distribution of susceptibility. This contrasts with static field measurements, which can be reproduced by a surface layer of magnetization. We illustrate the effectiveness of the inversion algorithm on synthetic and field data and show also the importance of knowing the background conductivity. In practical circumstances, where there is no a priori information about conductivity distribution, a simultaneous inversion of EM data to recover both electrical conductivity and susceptibility will be required.  相似文献   

5.
The derivation of seismic reflection and transmission coefficients is generally based on the assumption that the medium parameters behave as step functions of depth, at least in a finite region around the interface. However, outliers observed in well logs generally behave quite differently from step functions. In this paper we represent an interface by a self-similar singularity, embedded between two homogeneous half-spaces, and we derive its frequency-dependent normal-incidence reflection and transmission coefficients. For ω  → 0 the expressions for the coefficients reduce to those for a discrete boundary between two homogeneous half-spaces; for ω → ∞ they become frequency-independent. These asymptotic expressions have a relatively simple form and depend on the singularity exponent α .
  The exact as well as the asymptotic expressions are used to evaluate the time-domain reflection and transmission responses of a self-similar interface. Finally, we use a numerical method to model the response of a smoothed version of a self-similar interface (note that the velocity of a smoothed singularity remains finite). It turns out that smoothing has hardly any effect on the response, provided that the smoothing does not affect the scales corresponding to the seismic frequency range.  相似文献   

6.
Knowledge of the age of undated events is not null if a time-order relationship can be found among these events. The knowledge of such a time-ordered sequence can be formalized by using non-informative (uniform) prior probability densities for the ages of undated events and Bayes' theorem to introduce the time-order relationship condition. We show that the conditional probability densities of the ages of events of unknown age are given by various forms of Euler's beta distribution. These distributions yield an estimate of the probability for an undated event to occur in a given age interval.
  We use this method to propose appropriate probabilistic representations of our actual knowledge of the dating of the magnetic polarity reversals during the Cenozoic. These representations take into account the uncertainties arising from irregularities in accretion process and from the quality of a few calibration points. Both types of uncertainties generate large ambiguities in the age of magnetic reversals, which should be taken into consideration when the geomagnetic polarity timescale is used for dating purposes. We propose to use the entropy function to quantify these ambiguities.  相似文献   

7.
We use theoretical and numerical methods to investigate the general pore-fluid flow patterns near geological lenses in hydrodynamic and hydrothermal systems respectively. Analytical solutions have been rigorously derived for the pore-fluid velocity, stream function and excess pore-fluid pressure near a circular lens in a hydrodynamic system. These analytical solutions provide not only a better understanding of the physics behind the problem, but also a valuable benchmark solution for validating any numerical method.
  Since a geological lens is surrounded by a medium of large extent in nature and the finite element method is efficient at modelling only media of finite size, the determination of the size of the computational domain of a finite element model, which is often overlooked by numerical analysts, is very important in order to ensure both the efficiency of the method and the accuracy of the numerical solution obtained. To highlight this issue, we use the derived analytical solutions to deduce a rigorous mathematical formula for designing the computational domain size of a finite element model. The proposed mathematical formula has indicated that, no matter how fine the mesh or how high the order of elements, the desired accuracy of a finite element solution for pore-fluid flow near a geological lens cannot be achieved unless the size of the finite element model is determined appropriately.
  Once the finite element computational model has been appropriately designed and validated in a hydrodynamic system, it is used to examine general pore-fluid flow patterns near geological lenses in hydrothermal systems. Some interesting conclusions on the behaviour of geological lenses in hydrodynamic and hydrothermal systems have been reached through the analytical and numerical analyses carried out in this paper.  相似文献   

8.
This paper presents the application of the electrical resistivity tomography (ERT) method to the investigation of the Tertiary maar structure of Baruth (Germany) known from previous gravimetric surveys. ERT was applied to support the optimum location for a palaeoclimatological drill hole.
  Special modifications of data acquisition, signal processing and inversion are introduced to adapt the method of ERT to the special requirements for the 3-D investigation of structures with horizontal extensions of 1  km or more. More than 5000 dipole–dipole combinations were recorded at three concentric circular electrode arrangements using stand-alone transient data acquisition systems (RefTek).
  We present a fast approximate imaging technique based on the simultaneous iterative reconstruction technique (SIRT). As the complete calculation of the inverse Frechét matrix is avoided, the algorithm is especially suitable for large data and model spaces, where complete inversion is beyond the limits of available computing hardware. The single-step method is applicable to arbitrary irregular electrode layouts. Synthetic tests show that the imaging procedure reconstructs the main features of the subsurface.
  A low-resistivity body could be interpreted as limnic sediments filling the interior of the Tertiary maar crater. Considering the horizontal resistivity gradient, estimates for the lateral and depth extents of the structure were made. An optimum position for a palaeoclimatological borehole was found, and was in good agreement with the gravimetric minimum.  相似文献   

9.
Joint inversion of receiver function and surface wave dispersion observations   总被引：16，自引：0，他引：16  
We implement a method to invert jointly teleseismic P wave receiver functions and surface wave group and phase velocities for a mutually consistent estimate of earth structure. Receiver functions are primarily sensitive to shear wave velocity contrasts and vertical traveltimes, and surface wave dispersion measurements are sensitive to vertical shear wave velocity averages. Their combination may bridge resolution gaps associated with each individual data set. We formulate a linearized shear velocity inversion that is solved using a damped leastsquares scheme that incorporates a priori smoothness constraints for velocities in adjacent layers. The data sets are equalized for the number of data points and physical units in the inversion process. The combination of information produces a relatively simple model with a minimal number of sharp velocity contrasts. We illustrate the approach using noisefree and realistic noise simulations and conclude with an inversion of observations from the Saudi Arabian Shield. Inversion results for station SODA, located in the Arabian Shield, include a crust with a sharp gradient near the surface (shear velocity changing from 1.8 to 3.5 km s1 in 3 km) underlain by a 5kmthick layer with a shear velocity of 3.5 km s1 and a 27kmthick layer with a shear velocity of 3.8 km s1, and an upper mantle with an average shear velocity of 4.7 km s1. The crustmantle transition has a significant gradient, with velocity values varying from 3.8 to 4.7 km s1 between 35 and 40 km depth. Our results are compatible with independent inversions for crustal structure using refraction data.  相似文献   

10.
A high-quality aeromagnetic survey of northern Brittany (line spacing 250  m; flight elevation 150  m) has been used to delineate the Lower Carboniferous St Malo dyke swarm in detail. The dyke swarm fans at its northern and southern ends, and is affected by N60° sinistral transcurrent faults. After restoration of these offsets, the full structure trends in a N–S direction. Small dykes are not imaged, and only one-third of the swarm is evidenced by magnetism. Gravity and magnetic modelling shows that the swarm overlies a single N–S elongated magma chamber.
  The distribution of K max inclinations of anisotropy of magnetic susceptibility (AMS) suggests that the dykes display a fanning magma flow in section. Computed K max inclinations are usually shallower than the measured geological dips, probably because the flow becomes more disturbed as the dyke becomes shallower. We observe that the mean susceptibility values increase when the magma flow is steeper than about 70°.
  A detailed cross-section of St Briac dyke, which is part of the St Malo dyke swarm, reveals that the main carrier of magnetization is magnetite in the centre of the dyke and magnetite + maghaemite on its rims. The emplacement of the St Malo dyke swarm could have been contemporaneous with the N60° shearing which displaced the dyke swarm by about 20  km. The dyke swarm is cut at its southern end by an E–W-orientated fault which probably acted during Late Carboniferous times.  相似文献   

11.
New magnetotelluric data from the Münchberg Gneiss complex in Southern Germany reveal a zone of extremely high electrical conductivity. 1-D modelling of the data is justified in the period range 0.01 to 10  s. At least three layers are required to explain the steepness of the apparent resistivity curves, and the best-fitting models comprise four layers with successively higher conductivities. The layers of highest conductivity at depths between 2.2 and 3.6  km correlate with pronounced bands of high seismic reflectivity (profile DEKORP 85-4N). The Münchberg complex is today widely recognized as a tectonic klippe, consisting of rocks whose metamorphic and stratigraphic order is inverted rather than overturned. The material was transported into its present position by predominantly horizontal tectonic forces along shear zones. We interpret the high conductivity and high reflectivity as remnants of this transport process.  相似文献   

12.
We present a neural network approach to invert surface wave data for a global model of crustal thickness with corresponding uncertainties. We model the a posteriori probability distribution of Moho depth as a mixture of Gaussians and let the various parameters of the mixture model be given by the outputs of a conventional neural network. We show how such a network can be trained on a set of random samples to give a continuous approximation to the inverse relation in a compact and computationally efficient form. The trained networks are applied to real data consisting of fundamental mode Love and Rayleigh phase and group velocity maps. For each inversion, performed on a 2°× 2° grid globally, we obtain the a posteriori probability distribution of Moho depth. From this distribution any desired statistic such as mean and variance can be computed. The obtained results are compared with current knowledge of crustal structure. Generally our results are in good agreement with other crustal models. However in certain regions such as central Africa and the backarc of the Rocky Mountains we observe a thinner crust than the other models propose. We also see evidence for thickening of oceanic crust with increasing age. In applications, characterized by repeated inversion of similar data, the neural network approach proves to be very efficient. In particular, the speed of the individual inversions and the possibility of modelling the whole a posteriori probability distribution of the model parameters make neural networks a promising tool in seismic tomography.  相似文献   

13.
Summary. The problem of determination of the electrical conductivity of the Earth from the geomagnetic induction data is formulated as that of finding the coefficients of the continued fraction expansion of a certain rational fraction representation of the total response (or impedance) of the medium at the surface when this quantity is given for N frequencies. The coefficients of expansion are related to the conductivities of N layers of constant attenuation, where, within each layer the conductivity is assumed to be constant. Thus in this approach the conductivity profile resulting from the inversion of the response function is given as a series of step functions.  相似文献   

14.
In this paper, we will report on the application of Bayesian inference to DC resistivity inversion for 1-D multilayer models. The posterior probability distribution is explored through a Markov process based upon a Gibbs's sampler. The process would lead to unrealistic estimates without additional prior information, which takes the form of a second Markov chain where the transition kernel corresponds to a smoothness constraint. The outcomes are posterior marginal probabilites for each parameter, as well as, if required, joint probabilities for pairs of parameters. We will discuss the main properties of the method in the light of a theoretical example and illustrate its capabilities with some field examples taken from various contexts.  相似文献   

15.
About six separately orientated cores were collected at each of 14 sites distributed throughout the arcuate, west-dipping, 6  km thick, Freetown layered igneous complex. Alternating field and thermal demagnetization both isolate a stable component of remanent magnetism which corresponds to a palaeomagnetic south pole from 13 sites (nine reverse, four normal polarity) at 82.9°S, +32.7°E ( α 95 = 5.6°). This is indistinguishable from that reported in 1971 based on alternating field demagnetization of cores from 10 orientated hand samples.
  The difference between the Freetown pole (age: 193 ± 3  Ma) and other mid-Jurassic poles from West Africa could be due to its greater age. The difference between the whole West African Jurassic pole group and the Karoo pole from southern Africa, however, suggests moderate (∼10°) differential rotation of West Africa relative to the Kaapvaal craton.
  A prevalent magnetic foliation fabric coincides generally with the petrological layering, as might be expected, but a ubiquitous magnetic lineation is predominantly down-dip. This is compatible with a down-dip pyroxene lineation reported to be present in some field outcrops, and interpreted in terms of late-stage deformation during the slow crystallization and cooling of the large igneous body. However, a fold test shows that the igneous layering had already achieved its present attitude before the Complex cooled to ∼570 °C (the maximum blocking temperature of the characteristic remanence).  相似文献   

16.
In the pseudo-Thellier method for relative palaeointensity determinations (Tauxe et al. 1995) the slope of the NRM intensity left after AF demagnetization versus ARM intensity gained at the same peak field is used as a palaeointensity measure. We tested this method on a marine core from the Azores, spanning the last 276  kyr. We compared the pseudo-Thellier palaeointensity record with the conventional record obtained earlier by Lehman et al . (1996 ), who normalized NRM by SIRM. The two records show similar features: intensity lows with deviating palaeomagnetic directions at 40–45  ka and at 180–190  ka. The first interval is associated with the Laschamps excursion, while the 180–190  ka low represents the Iceland Basin excursion (Channell et al. 1997). The pseudo-Thellier method, in combination with a jackknife resampling scheme, provides error estimates on the palaeointensity.
  Spectral analysis of the rock magnetic parameters and the palaeointensity estimates shows orbitally forced periods, particularly 23  kyr for climatic precession. This suggests that palaeointensity is still slightly contaminated by climate. Fuzzy c -means cluster analysis of rock magnetic and geochemical parameters yields a seven-cluster model of predominantly calcareous clusters and detrital clusters. The clusters show a strong correlation with climate, for example samples from detrital clusters predominantly appear during rapid warming. Although both the pseudo-Thellier palaeointensity m a and fuzzy clusters show climatic influences, we have not been able to find an unambiguous connection between the clusters and m a .  相似文献   

17.
Precise time and facies correlations between drilled holes are fundamental for a better understanding of the geological evolution of sedimentary basins. A downhole magnetic measurement device called the geological high-sensitivity magnetic tool (GHMT) has been run within two wells drilled by Gaz de France in the Landes oil-field (southwest France) as part of a gas storage exploration program. The method of interpretation of downhole magnetic measurements yielded a magnetostratigraphy within each well, allowing absolute dating and time correlations between the wells.
  Magnetic susceptibility and natural gamma ray intensity are useful parameters for establishing high-resolution lithological correlations at a basin scale. We present a correlation parameter established from a simultaneous analysis of the susceptibility and the gamma ray logs within each well. The correlation parameter appears to provide a new tool for delineating lithological elements when local lithological changes are too subtle to show clear well-to-well correlations either from susceptibility logs or from gamma ray logs. This new approach is interpreted as a sensitive way to detect relative variations between the detrital and clay content of the penetrated sediment.  相似文献   

18.
Gravity changes are presented from a series of field microgravity surveys conducted at Mt Etna between August 1994 and November 1996, a period including the 1995–1996 explosive summit activity. Data were collected along a microgravity network of 69 stations at a monthly to annual sampling rate, depending on each subarray of the network.
  Results show that seasonal changes in water level within the volcano may induce gravity changes of up to 20  μgal on Etna's southern slope, and indicate that significant magma movement occurred within and below Etna's edifice between 1994 and 1996. In particular, between September 1994 and October 1995, a mass increase of 2 × 1010  kg occurred 2000  m beneath the summit craters. Between October 1995 and July 1996 this mass was lost, while another 2 × 1010  kg was injected at about 1000  m  a.s.l. into the 1989 fracture system. From the gravity data alone, it is not possible to distinguish whether the first shallow intrusion (1994–1995) was then injected laterally into the 1989 fracture, or summit activity was fed by the first shallow intrusion, while new magma entered the 1989 fracture system.
  While magma was being redistributed within the volcanic edifice, measurements along an E–W-trending profile on the southern slope of the volcano detected some 1.5 × 1011  kg of magma accumulating 2–3  km below sea level between October 1995 and November 1996.  相似文献   

19.
Inversion for multiple parameter classes   总被引：2，自引：0，他引：2  
Many geophysical data, such as the frequencies of the free oscillations of the Earth, depend on more than one type of model parameter. For inverse problems depending on multiple parameter classes, an iterative solution procedure is introduced in which each parameter class can be treated in the same way. This approach has considerable advantages where a large number of parameters are employed, but can still be useful for smaller systems.
  The iteration by parameter class commences by solving for the direct dependence on a particular parameter class, and at subsequent iterations the cross-dependences between classes are introduced. The update affects only the right-hand side of the equations, and, because the same sets of equations have to be solved at each iteration, an efficient computational implementation can be made. The largest set of equations that has to be solved at a time corresponds to the number of variables in an individual parameter class rather than the full set of parameters, which confers substantial computational benefits for very large problems.  相似文献   

20.
Anisotropy in multi-offset deep-crustal seismic experiments   总被引：1，自引：0，他引：1  
Modelling of deep-seismic wide-angle data commonly assumes that the Earth is heterogeneous and isotropic. It is important to know the magnitudes of errors that may be introduced by isotropic-based wide-angle models when the Earth is anisotropic. It is equally important to find ways of detecting anisotropy and determining its properties.
  This paper explores the errors introduced by interpreting anisotropic seismic data with isotropic models. Errors in P -wave reflector depths are dependent on the magnitude of the velocity anisotropy and the direction of the fast axis. The interpreted, isotropic, model velocity function is found to correspond closely to the horizontal velocity of the anisotropic medium. An additional observed parameter is the time mismatch , which we define to be the difference between the vertical two-way traveltime to a reflector and the time-converted wide-angle position of the reflector. The magnitude of the time mismatch is typically <1.0  s (when the whole crust is anisotropic) and is found to be closely related to the magnitude and sign of the anisotropic anellipticity. The relationships are extendible to more complicated models, including those with vertical velocity gradients, crustal zonation, and lower symmetry orders.
  A time mismatch may be symptomatic of the presence of anisotropy. We illustrate the observation of a time mismatch for a real multi-offset seismic data set collected north of Scotland and discuss the implications for crustal anisotropy in that region.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号