首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 359 毫秒
1.
为查明雷琼火山群中九斗洋干玛珥湖的第四纪地层空间展布,以及埋藏古火山形态,为后续研究提供地质背景资料,采用高效、便捷的高密度电阻率法对九斗洋干玛珥湖湖盆区进行勘测,并结合钻孔岩芯的研究进行验证。得出的主要结论为：1）盆地深部地层最高电阻率可达到300 Ω?m以上,盆地中央基岩以下存在高阻地质体,推测为岩浆通道位置。岩浆经过通道溢出后横向展布,覆于沉积地层上,随后在火山口内形成湖泊,堆积第四纪湖相沉积,其火山机构与玛珥湖型火山一致;2）九斗洋干玛珥湖第四纪松散沉积地层的电阻率通常＜60 Ω?m,泥炭层电阻率最低;3）第四纪沉积总厚度为25~50 m,盆地松散沉积的岩浆岩基底总体平整,但尚有波状起伏。研究结果表明,高密度电阻率法结合钻孔验证是研究火山机构形态及火山口湖第四纪沉积空间展布的有效方法。  相似文献   

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

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

4.
Arctic tundra surfaces are dominated by a variety of patterned ground forms. Whereas a large number of studies have described morphology, structure and processes of patterned ground, few have monitored detailed patterned ground dynamics and subsurface environments continuously. We applied electrical resistivity tomography (ERT) to understand near‐surface conditions of two types of patterned ground, ice‐wedge polygons and mudboils in Svalbard, where periglacial processes associated with permafrost are intensively monitored. Automated monitoring shows surface movement characterized by annual cycles of frost heave and thaw settlement, the amounts and rates of which are influenced by the intensity of ice segregation. A time series of ERT shows (1) a distinct resistivity boundary delimiting the active‐layer depth, (2) seasonal variation in resistivity controlled by thermo‐hydrological dynamics and (3) spatial variation in resistivity reflecting desiccation in summer and intensive ice segregation in winter. These results demonstrate ERT as a useful complementary technique for monitoring active‐layer depths and near‐surface hydrological conditions at periglacial patterned ground sites, where automated soil thermal and moisture measurements are limited.  相似文献   

5.
Preferential flow: first results of a full-scale flow model   总被引：1，自引：0，他引：1  
The main goal of a joint project undertaken by the geophysical and hydrological research units of Kiel University is to study preferential flow in a large open-air, full-scale model, looking in particular at near-surface penetration and flow of water through the unsaturated vadose zone. An artificial irrigation device is installed in place of natural rain, and a homogeneous sand body is used instead of natural soil. This provides a reference model for future field experiments. Inside the sand body there are a large number of geophysical and hydrological sensors to measure DC resistivity (using various electrode configurations), water content and water potential (using TDR and tensiometer instruments, respectively). A ground-penetrating radar (GPR) system is installed at the surface, whereas at the bottom several containers and a thin gravel layer are embedded to measure the flow arrival and the discharge of water. Irrigation is varied in intensity, time, area, and salt content (tracer).
  Results of the first six experiments show that the percolation of intruding water can be followed by all techniques and percolation is finally controlled by the discharge measurements. These display some undulations and variations of the water 'front' and agree with the measurements of all other sensors. The redundancy achieved by the use of multiple methods was intended to enable an assessment of the reliability of the techniques used. The true values of electrical resistivity before and after irrigation reflect the distribution of water saturation within the sand body. A numerical 3-D inversion of the apparent resistivity provides information regarding future field experiments, in which it will be possible to install only some of the sensors in order to preserve the natural structure of the soil.  相似文献   

6.
A magnetotelluric study of the Alpine Fault, New Zealand   总被引：1，自引：0，他引：1  
Magnetotelluric soundings have been made at seven locations on a 4  km profile crossing the Alpine Fault in the South Island of New Zealand. The 'distortion' techniques of Groom & Bailey (1989 ) and Lilley (1998a , b ) have been used to derive regional apparent resistivity and phase curves that correspond to electromagnetic induction in orientations parallel and perpendicular to the fault. 2-D inversion of the regional responses reveals that a narrow (<1  km wide) conductive zone is associated with the Alpine Fault. This conductor is most probably related to the heating of deep circulating meteoric water in a region in which enhanced temperatures occur at shallow depth due to the tectonic uplift of the Southern Alps.  相似文献   

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

8.
According to the theory of isostasy, the Earth has a tendency to deform its surface in order to reach an equilibrium state. The land-uplift phenomenon in the area of the Fennoscandian Shield is thought to be a process of this kind. The geoid, as an equipotential surface of the Earth's gravity field, contains information on how much the Earth's surface departs from the equilibrium state. In order to study the isostatic process through geoidal undulations, the structural effects of the crust on the geoid have to be investigated.
  The structure of the crust of the Fennoscandian Shield has been extensively explored by means of deep seismic sounding (DSS). The data obtained from DSS are used to construct a 3-D seismic-velocity structure model of the area's crust. The velocity model is converted to a 3-D density model using the empirical relationship that holds between seismic velocities and crustal mass densities. Structural effects are then estimated from the 3-D density model.
  The structural effects computed from the crustal model show that the mass deficiency of the crust in Fennoscandia has caused a geoidal depression twice as deep as that observed from the gravimetric geoid. It proves again that the crust has been isostatically compensated by the upper mantle. In other words, an anomalously high-density upper mantle must exist beneath Fennoscandia.  相似文献   

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

10.
Aquifer-bearing intermontane sedimentary basins of the Sierras Pampeanas in the northwest of Argentina are in general very deep and narrow and contain economically important deposits of Tertiary sediments. This paper presents the results of a study to characterize the sedimentary basin bounded to the west by Sierra de Famatina and to the east by Sierra de Velasco, where an electromagnetic sensing technique, the magnetotelluric (MT) method, was applied. 12 MT sites were deployed along a 30  km E–W transect. Some of the data collected were used to derive a 2-D resistivity model of the basin. The model shows a subsurface trough with a thick (approximately 8  km) sedimentary sequence above it. Anomalous behaviour of the E–W electric-field component ( E y ) was detected in the period range 1–100  s, where the amplitude of this component was below the instrumental noise level. The cause of this anomaly is not known, but it might be due to the presence of an embedded conducting body between 8 and 10  km, which would give rise to N–S current channelling.  相似文献   

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

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

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

14.
b
This paper describes a method in which vertical resistivity sections are generated tomographically from measurements on a linear array of equally spaced electrodes inserted at the ground surface. The array is multiplexed to a resistivity meter which gathers one set of all possible independent apparent resistivity measurements and the geophysical section is then reconstructed by backprojecting these weighted data, along equipotentials, into the subsurface. The technique has been evaluated numerically and in field trials over shallow archaeological structures at Fountains Abbey.  相似文献   

15.
Collocated magnetotelluric (MT) and seismic profiling is emerging as a necessary combined approach for deep and near-surface imaging but the resulting experimental data are typically interpreted separately since no production programs exist for multidimensional joint inversion of MT and seismic data. We present a joint 2-D inversion approach for imaging collocated MT and seismic refraction data with cross-gradient structural constraints. We describe the main features of the algorithm and first apply it to synthetic data generated for a hypothetical complex geological model. For the synthetic data, we find that the scheme leads to models with remarkable structural resemblance and improved estimates of electrical resistivity and seismic velocity. We apply the scheme to near-surface field data to test the consistency of a previously suggested resistivity–velocity interrelationship and its potential use for subsurface lithofacies discrimination or structural classification. The MT-seismic relationship is found to be in excellent accord with that derived previously for DC resistivity and seismic data set at the test site. Our results suggest that joint MT-seismic cross-gradient imaging leads to improved characterization of heterogeneous geological targets at near-surface to mantle depths.  相似文献   

16.
A magnetic study was carried out on lacustrine sediments from the Zoigê basin, Tibetan Plateau, in order to obtain a better understanding of palaeoclimatic changes there. Gyromagnetic remanence (GRM) acquisition is unexpectedly observed during static three-axis alternating field (AF) demagnetization in about 20 per cent of a large number of samples. X-ray diffraction (XRD) analysis on a magnetic extract clearly shows that greigite is the dominant magnetic mineral carrier. Scanning electron microscopy (SEM) reveals that the greigite particles are in the grain size range of 200–300  nm, possibly in the single-domain state. Greigite clumps of about 3  μm size are sealed by silicates. Fitting of XRD peaks yields a crystalline coherence length of about 15  nm, indicating that the particles seen in the SEM are polycrystalline.
  GRM intensities of most samples are of the same order as the NRM, while others show much stronger GRM although their magnetic properties are similar. Variation of the demagnetization sequence confirms that GRM is mainly produced perpendicular to the AF direction. The anisotropy direction can be derived from GRM, but more systematic studies are needed for detailed conclusions. An attempt to correct for GRM failed due to high GRM intensities and because smaller GRM acquisition was also found along the demagnetization axis. Behaviours of acquisition and AF demagnetization of GRM are comparable with those of NRM, ARM, IRM, indicating fine grain sizes of remanence carriers.  相似文献   

17.
The discrimination between electrolytic and electronic conductors is highly relevant to geological modelling as it allows conclusions to be drawn about the formation and mineral composition of rocks. The induced polarization (IP) method, which compares the electric current injected into the ground with the corresponding earth potential differences can be used for this purpose.
  This paper describes a new method based on the theory that non-linear electrochemical processes on the surface of electronic conductors are responsible for non-linear IP (NLIP) phenomena. This results in multiples of the fundamental frequency being observed in the telluric voltage spectra when a monochromatic current signal is fed into the ground. The non-linearity of the current–voltage characteristic is most effectively described by a spectral method.
  A laboratory experiment was carried out, using an electrolytic trough with a small graphite cylinder serving as an electronic conductor, which clearly demonstrated the validity of the method. A field experiment was undertaken at a borehole of approximately 450  m depth, located in the transition zone of the Tepla-Barrandium and Moldanubicum in East Bavaria. A sinusoidal current was injected into the ground using a logging tool at depths varying between 150 and 450  m. The corresponding potential differences were simultaneously observed along a profile on the surface. Field and laboratory results show a striking similarity. It can be concluded that an extensive electronic conductor—probably graphite—is steeply dipping southwards meeting the borehole at approximately 310  m depth.  相似文献   

18.
Magnetotelluric data from the backarc of the Central Andes in NW Argentinawere re-examined by employing impedance tensor decomposition and 2-D inversion and modelling techniques. The data in the period range of 50–15 000 s were collected on a profile of 220 km length reaching from the Eastern Cordillera across the Santa Barbara System to the Andean foreland of the Argentinean Chaco.
After a dimensionality analysis, data from most sites were treated as regional 2-D. The exception was the eastern section of the profile, where the magnetotelluric transfer functions for periods ≤ 1000 s reflect a 3-D earth. Application of two tensor decomposition schemes yielded a regional strike direction of N–S, which is the azimuth of the Central Andean mountain chains. Several 2-D models were obtained by pseudo- and full 2-D Occam inversion schemes. Special emphasis was placed on the inversion of phase data to reduce the influence of static shifts in the apparent resistivity data. The smooth inversion models all show a good conductor at depth. A final model was then calculated using a finite element forward algorithm.
The most prominent feature of the resulting model is a conductor which rises from depths of 180 km below the Chaco region to 80 km beneath the Santa Barbara System and the Eastern Cordillera. Its interpretation as a rise of the electrical asthenosphere is supported by seismic attenuation studies. Magnetotelluric results, surface heat-flow distribution in the area, and the electrical properties of crustal and mantle rocks suggest that the upper mantle is predominantly ductile beneath the Eastern Cordillera and the western Santa Barbara System. This generally agrees with anelastic seismic attenuation models of the area and is useful in discriminating between models of Q quality factor distribution.  相似文献   

19.
From ACH tomographic models to absolute velocity models   总被引：2，自引：0，他引：2  
The ACH method, a widely used tomographic inverse method, is characterized by the use of relative residuals in order to avoid possible biases coming from outside the target volume. The ACH method thus does not really retrieve the 3-D structure of the target volume, but instead leads to velocity contrasts relative to the layer average of the velocity, this average value remaining unknown ( Aki et al. 1977 ). Two artefacts derive from this particularity: (1) velocity contrasts are known only in the horizontal direction and it is not possible, in a strict mathematical sense, to estimate the contrasts in the vertical direction with ACH alone; (2) negative anomalies are often interpreted as low velocities, whereas negative anomalies may correspond to high velocities if the average value of the corresponding layer is sufficiently high. The converse is true of positive anomalies. We show with synthetic data how these artefacts can affect the interpretation of tomographic images. We propose to correct the artefacts by reintroducing the 1-D regional average model, and show in synthetic experiments how effective this correction can be.
  The application of this procedure to data recorded in the Kunlun region shows that the retrieval of the absolute values of the 3-D velocity model is helpful for interpreting the tomographic images and better defining which features are anomalous.  相似文献   

20.
  相似文献   

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

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