首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 343 毫秒
1.
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.  相似文献   

2.
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.  相似文献   

3.
This paper presents a new derivative-free search method for finding models of acceptable data fit in a multidimensional parameter space. It falls into the same class of method as simulated annealing and genetic algorithms, which are commonly used for global optimization problems. The objective here is to find an ensemble of models that preferentially sample the good data-fitting regions of parameter space, rather than seeking a single optimal model. (A related paper deals with the quantitative appraisal of the ensemble.)
  The new search algorithm makes use of the geometrical constructs known as Voronoi cells to derive the search in parameter space. These are nearest neighbour regions defined under a suitable distance norm. The algorithm is conceptually simple, requires just two 'tuning parameters', and makes use of only the rank of a data fit criterion rather than the numerical value. In this way all difficulties associated with the scaling of a data misfit function are avoided, and any combination of data fit criteria can be used. It is also shown how Voronoi cells can be used to enhance any existing direct search algorithm, by intermittently replacing the forward modelling calculations with nearest neighbour calculations.
  The new direct search algorithm is illustrated with an application to a synthetic problem involving the inversion of receiver functions for crustal seismic structure. This is known to be a non-linear problem, where linearized inversion techniques suffer from a strong dependence on the starting solution. It is shown that the new algorithm produces a sophisticated type of 'self-adaptive' search behaviour, which to our knowledge has not been demonstrated in any previous technique of this kind.  相似文献   

4.
20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.  相似文献   

5.
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.  相似文献   

6.
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.  相似文献   

7.
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 .  相似文献   

8.
A palaeomagnetic investigation has been carried out of rocks from the eastern part of the Voronezh Massif, which constitutes, together with the Ukrainian Shield, the Sarmatian segment in the southern part of the East European Craton. The samples were collected in a quarry close to the town of Pavlovsk (50.4°N, 40.1°E), where a syenitic-granitic body intrudes Archaean units. U–Pb (zircon) dating has yielded an age of 2080  Ma for the intrusion.
  Two characteristic magnetic components, A and B, were isolated by thermal and alternating-field demagnetization. Component A was obtained from granites and quartz syenites (11 samples) and has a mean direction of D = 229°, I = 28°, and a pole position at 12°N, 172°E. This pole is close to a contemporary mean pole (9°N, 187°E) for the Ukrainian Shield, which implies that the Voronezh Massif and the Shield constituted a single entity at 2.06  Ga. These poles differ from contemporaneous poles of the Fennoscandian Shield, indicating that the relative positions of the two shields were different from their present configuration about 2100  Myr ago.
  A component B, isolated only in quartz monzonites (five samples), has a mean direction D = 144°, I = 49°, and a pole position at 4°N, 251°E, which is close to late Sveconorwegian (approximately 900  Ma) poles for Baltica. This suggests that the East European Craton was consolidated some time between 2080 and 900  Ma. Comparison with other palaeomagnetic data permit us to narrow this time span to 1770–1340  Ma.  相似文献   

9.
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.  相似文献   

10.
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.  相似文献   

11.
Estimates of past climate derived from borehole temperatures are assuming a greater importance in context of the millennial temperature variation debate. However, recovery of these signals is usually performed with regularization which can potentially lead to underestimation of past variation when noise is present. In this work Bayesian inference is applied to this problem with no explicit regularization. To achieve this Reversible Jump Markov chain Monte Carlo is employed, and this allows models of varying complexity (i.e. variable dimensions) to be sampled so that it is possible to infer the level of ground surface temperature (GST) history resolution appropriate to the data. Using synthetic examples, we show that the inference of the GST signal back to more than 500 yr is robust given boreholes of 500 m depth and moderate noise levels and discuss the associated uncertainties. We compare the prior information we have used with the inferred posterior distribution to show which parts of the GST reconstructions are independent of this prior information. We demonstrate the application of the method to real data using five boreholes from southern England. These are modelled both individually and jointly, and appear to indicate a spatial trend of warming over 500 yr across the south of the country.  相似文献   

12.
The inversion of recent borehole temperatures has proved to be a successful tool to determine ancient ground surface temperature histories. To take into account heterogeneity of thermal properties and their non-linear dependence on temperature itself, a versatile 1-D inversion technique based on a finite-difference approach has been developed. Regularization of the generally ill-posed problem is obtained by an appropriate version of Tikhonov regularization of variable order. In this approach, a regularization parameter has to be determined, representing a trade-off between data fit and model smoothness. We propose to select this parameter by generalized cross-validation. The resulting technique is employed in case studies from the Kola ultradeep drilling site, and another borehole from northeastern Poland. Comparing the results from both sites corroborates the hypothesis that subglacial ground surface temperatures as met in Kola often are much higher than the ones in areas exposed to atmospheric conditions (Poland).  相似文献   

13.
在离子膜法制碱过程中,电解盐水中的碘会以碘酸盐和高碘酸盐的形式在离子膜内聚积,造成膜堵塞、膜寿命缩短、电流效率降低、电压升高等问题。现有的盐水除碘技术有氧化吸附法、氧化吹除法和还原吹除法,这些方法均需将溶液中的离子碘转化为单质碘。碘离子吸附剂可直接吸附去除溶液中的碘离子,有望用于盐水除碘,既可避免投入氧化剂,又缩短了除碘工艺流程,其盐水除碘性能有待进一步研究。  相似文献   

14.
We use Monte Carlo Markov chains to solve the Bayesian MT inverse problem in layered situations. The domain under study is divided into homogeneous layers, and the model parameters are the conductivity of each layer. We use an a priori distribution of the parameters which favours smooth models. For each layer, the a priori and a posteriori distributions are digitized over a limited set of conductivity values.
  The Markov chain relies on updating the model parameters during successive scanning of the domain under study. For each step of the scanning, the conductivity is updated in one layer given the actual value of the conductivity in the other layers. Thus we designed an ergodic Markov chain, the invariant distribution of which is the a posteriori distribution of the parameters, provided the forward problem is completely solved at each step.
  We have estimated the a posteriori marginal probability distributions from the simulated successive values of the Markov chain. In addition, we give examples of complex magnetotelluric impedance inversion in tabular situations, for both synthetic models and field situations, and discuss the influence of the smoothing parameter.  相似文献   

15.
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.  相似文献   

16.
Taiwan Chelungpu-fault Drilling Project (TCDP) was initiated to understand the physical mechanisms involved in the large displacements of the 1999 Taiwan Chi-Chi earthquake. Continuous measurements of cores (including laboratory work) and a suite of geophysical downhole logs, including P - and S -wave sonic velocity, gamma ray, electrical resistivity, density, temperature, electrical borehole images and dipole-shear sonic imager, were acquired in Hole-A over the depth of 500–2003 m. Integrated studies of cores and logs facilitate qualitative and quantitative comparison of subsurface structures and physical properties of rocks. A total of 10 subunits were divided on the basis of geophysical characteristics. Generally, formation velocity and temperature increase with depth as a result of the overburden and thermal gradient, respectively. Gamma ray, resistivity, formation density, shear velocity anisotropy and density-derived porosity are primarily dependent on the lithology. Zones with changes of percentage of shear wave anisotropy and the fast shear polarization azimuth deduced from Dipole Shear-Imager (DSI) are associated with the appearance of fractures, steep bedding and shear zones. The fast shear wave azimuth is in good agreement with overall dip of the bedding (approximately 30° towards SE) and maximum horizontal compressional direction, particularly in the Kueichulin Formation showing strong shear wave velocity anisotropy. Bedding-parallel fractures are prevalent within cores, whereas minor sets of high-angle, NNW–SSE trending with N- and S-dipping fractures are sporadically distributed. The fault zone at depth 1111 m (FZA1111) is the Chi-Chi earthquake slip zone and could be a fluid conduit after the earthquake. The drastic change in fast shear wave polarization direction across the underlying, non-active Sanyi thrust at depth 1710 m reflects changes in stratigraphy, physical properties and structural geometry.  相似文献   

17.
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).  相似文献   

18.
More than 1000 seismic events in northern Europe at distances of up to 400  km from the detecting network are located using an optimization method in which the global minimum of the traveltime function residuals is searched for using an Interval Arithmetic (IA) method. Epicentres are determined using P waves detected by the Finnish national seismic network: up to 15 stations were used in the analysis. The IA results coincide with locations provided by the University of Helsinki bulletins with a median location bias of 7.6  km.
  A second data set of 59 explosions in the Siilinjärvi mine in central Finland was examined in detail, because the locations of the explosions were known exactly. In this case, the median difference of IA locations was 3.8  km from the average location of mine explosions, while all 59 events were located within 9  km of the 'true' epicentres. The corresponding median error of the University of Helsinki locations was smaller (3.2  km), but some Helsinki locations were well over 10  km from the mine. The convergence towards the global optimum using interval arithmetic was fast when compared with the conventional least-squares approaches for epicentre determinations.  相似文献   

19.
The Loma Prieta earthquake (magnitude 7.0), which occurred in October 1989 in central California, was preceded by a period during which the mean magnitude of background seismicity in a small region near the eventual epicentre was abnormally low. This period may have begun as early as 1979, and it continued until mid-1988, after which the mean magnitude increased to a higher than normal value until the main earthquake. These changes were observed in the seismicity of an area 40  km in radius, centred on the Loma Prieta epicentre, and are consistent with the predictions of fracture mechanics studies. The 1988 change correlates with a reported change in long-term strain.
  A procedure has been developed for resolving such temporal changes in seismicity using CUSUM statistics. It demonstrates that the anomaly was highly significant, on the basis of analyses of two independent catalogues. There was also a significant anomaly before the 1994 Northridge earthquake.
  The hypothesis that large earthquakes are preceded by periods in which the mean magnitude of background activity is abnormally low, in the immediate vicinity of the eventual epicentre, is a tantalizing one. The analysis tool examined here may be useful for resolving such changes. Care needs to be taken, however, in routine surveillance of earthquake populations that contain large aftershock sequences.  相似文献   

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号