首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
本文给出一种用积分方程法配合有限差分法计算位于导电介质中三维导体的电磁响应的新途径,在二层大地条件下给出了具体计算方法和结果。详细讨论了围岩介质的导电性对导体异常的影响。这对于开展低阻覆盖层地区的电磁法工作具有一定的指导意义。  相似文献   

2.
本文给出一种用积分方程法配合有限差分法计算位于导电介质中三维导体的电磁响应的新途径,在二层大地条件下给出了具体计算方法和结果。详细讨论了围岩介质的导电性对导体异常的影响。这对于开展低阻覆盖层地区的电磁法工作具有一定的指导意义。  相似文献   

3.
The rotating current EM method has been applied to the delineation of two conductive orebodies, Elura near Cobar, NSW, and Thalanga near Charter's Towers, Queensland. The field data were collected in the form of observations of the vertical magnetic field strength ratio and phase difference using a Turam-style receiver with twin vertical coils. By reconstituting this data back to the ring source field and phase, i.e. the observed Hz, phasor, it is possible to present contoured maps of the EM field. Anomaly phasors are obtained by subtracting theoretical phasors from the observed phasors in the complex plane of the Hz phasor. The theoretical phasors for the finite source are based on horizontally layered, half-space earth models, computed at each point of the survey grids, then normalized to selected points of the observed fields. Use is made of the intrinsic circular symmetry of the method in X–Y plots of field versus source-receiver distance to ascertain geoelectric parameters for the earth models. A steel picket fence at Thalanga is modelled by a line source grounded at each end and its Hz, phasor is removed by the same process. A considerable improvement in anomaly delineation is gained over previous Turam-style anomalies and the two survey examples illustrate the limitations of the method in the presence of a conductive overburden (Elura) and its abilities in the absence of a conductive overburden (Thalanga).  相似文献   

4.
We present a closed-form frequency-wave number (ω – k) Green’s function for a layered, elastic half-space under SH wave propagation. It is shown that for every (ω – k) pair, the fundamental solution exhibits two distinctive features: (1) the original layered system can be reduced to a system composed by the uppermost superficial layer over an equivalent half-space; (2) the fundamental solution can be partitioned into three different fundamental solutions, each one carrying out a different physical interpretation, i.e., an equivalent half-space, source image impact, and dispersive wave effect, respectively. Such an interpretation allows the proper use of analytical and numerical integration schemes, and ensures the correct assessment of Cauchy principal value integrals. Our method is based upon a stiffness-matrix scheme, and as a first approach we assume that observation points and the impulsive SH line-source are spatially located within the uppermost superficial layer. We use a discrete wave number boundary element strategy to test the benefits of our fundamental solution. We benchmark our results against reported solutions for an infinitely long circular canyon subjected to oblique incident SH waves within a homogeneous half-space. Our results show an almost exact agreement with previous studies. We further shed light on the impact of horizontal strata by examining the dynamic response of the circular canyon to oblique incident SH waves under different layered half-space configurations and incident angles. Our results show that modifications in the layering structure manifest by larger peak ground responses, and stronger spatial variability due to interactions of the canyon geometry with trapped Love waves in combination with impedance contrast effects.  相似文献   

5.
Conventional processes of extracting magnetotelluric signals from noisy records are reviewed: instrument noises and noises that are generated close to the detectors can be eliminated by the usual auto- and crosscorrelation processes. Identification of coherent noises, such as pulses due to field sources that are not uniform over at least 100 km in oil exploration or 1000 km in crustal studies, is much more tedious. The 5 components Hx, Hy, Hz, Ex, Ey, of the magnetotelluric field have been recorded in many areas in France at different periods of the year, (a) in non-uniform field sources in the vicinity of electric railways and of 50 cycle power lines, and (b) in areas of strong inhomogeneity at depth on the flanks of steep structures and near the sea shore. Means for detecting non-uniformity are reviewed. Measuring the vertical component of magnetic pulses is a good way of estimating field uniformity: if H vertical/H horizontal <10%, the uniform field assumption is valid, and the classical restitution formulas can be used; if H vertical/H horizontal > 10%, uniformity can not be assumed and there is some difficulty in deciding whether non-uniformity is due to the field source or to anisotropy or inhomogeneities at depth. Several ways to solve this difficulty are described. The reliability of calculation of actual resistivity at various depths is examined as a function of the precision of apparent resistivity measurements.  相似文献   

6.
The inductive response of a conducting horizontal cylinder embedded in a uniform earth is studied using numerical results obtained for an analytical solution for the problem of a conducting cylinder buried in a homogeneous earth for the case of a uniform inducing field. A check of the validity of the numerical results is made by a comparison with analogue model measurements for a number of cases. Numerical results for a range of cylinder radii (a = 1–10 km), depths of burial (d= 0–4 km), conductivity contrasts (σ2= 10?2-10 Sm?1), and source frequencies (f= 10?1-10?4 Hz) of interest in the interpretation of magnetotelluric field measurements are presented. The results indicate that for a uniform inducing field the conductivity and depth of burial of a horizontal cylindrical inhomogeneity are best determined through a measurement of the amplitudes Hy, Hz and Ex and the phases φy and Ψx.  相似文献   

7.
For the computation of the vertical component Hz of the magnetic field of a horizontal A.C. dipole lying on the earth's surface, a recurrence formula is presented for a horizontally stratified half space, to obtain the (n+ 1)-layer case from the w-layer case. By means of several computed diagrams for the two-layer case, Hz can be determined for different ratios of conductivity of the subsoil and that of the overburden. Thereby the distance from the dipole as well as the layer thickness h are expressed in terms of the wave length A of a plain wave in the overburden. Assuming a sufficiently large conductivity difference, the results show that evidence about the subsurface conditions can be obtained if the distance between the measuring coil and the dipole is of the order of A/3, and if the thickness h of the layer varies within the range A/100 < h < A/6. As an example for the 3-layer case, a nonconducting intermediate layer is assumed.  相似文献   

8.
Abstract

An investigation is made of steady thermal convection of a Boussinesq fluid confined in a vertically-mounted rotating cylinder. The top and bottom endwall disks are thermal conductors at temperatures Tt and Tb with δT = Tt ? Tb >0. The vertical sidewall has a finite thermal conductance. A Newtonian heat flux condition is adopted at the sidewall. The Rayleigh number of the fluid system is large to render a boundary layer-type flow. Finite-difference numerical solutions to the full Navier-Stokes equations are obtained. The vertical motions within the buoyancy layer along the sidewall induce weak meridional flows in the interior. Because of the Coriolis acceleration, the meridional flows give rise to azimuthal flows relative to the rotating container. Strong vertical gradients of azimuthal flows exist in the regions near the endwalls. As the stratification effect increases, concentration of flow gradients in thin endwall boundary layers becomes more pronounced. The azimuthal flow field exhibits considerable horizontal gradients. The temperature field develops horizontal variations superposed on the dominant vertical distribution. As either the sidewall thermal conductance or the stratification effect decreases, the temperature distribution tends to the profile varying linearly with height. Comparisons of the sizes of the dynamic effects demonstrate that, in the bulk of flow field, the vertical shear of azimuthal velocity is supported by the horizontal temperature gradient, resulting in a thermal-wind relation.  相似文献   

9.
A two-dimensional analysis is applied to examine the effect that a sloping bedrock half-space has on the amplification of an anti-plane shear wave. The direct boundary integral equation method is used for the two-dimensional analysis. The particular soil–rock configuration investigated includes a homogeneous soil layer underlain by a sloping rock half-space. The rock half-space dips for a horizontal distance L and then becomes horizontal so that the overlying soil layer has a thickness H that remains constant from this point to infinity. The materials in the soil–rock configuration are considered viscoelastic except in the rock half-space below soil layer thickness H, which is considered elastic. This limitation in damping is due to the correction used for the truncation of the half-space boundary. Four cases are used to study the relationship between rock slope and surface displacement, vertical, 1:2, 1:4, 1:8. Surface displacements are determined for each of these cases for half-space incidence angles of 90, 75, and 60°. To allow for applicability to a wide range of problems, results are determined as a function of dimensionless parameters. In addition, solutions from a one-dimensional analysis are compared with the results of the two-dimensional analysis to establish limits outside of which a one-dimensional analysis suffices.  相似文献   

10.
An axisymmetric model of convection in a rotating cylinder in an external uniform magnetic field has been considered. In the considered model, the meridional circulation is created by a nonuniform rotation of the lower boundary relative to the other boundaries. In the considered model, the time of formation of the stationary regime in the magnetic field considerably increases if the vertical density (compressibility) inhomogeneity is taken into account for Ekman numbers of E = E M = 3 × 10−3. This example shows that the compressibility of a medium should be taken into account in the convection and dynamics of the magnetic field when the magnetohydrodynamics of the Earth is analyzed.  相似文献   

11.
The direct boundary integral equation technique is used to study in-plane surface amplification of in-plane seismic body waves for the case of an inhomogeneity in a bedrock half-space. In the studied soil configuration, a soil layer rests on a rock half-space which includes a rock inclusion. The rock inclusion considered is a semi-infinite horizontal rock layer in which its upper boundary borders the soil layer. Materials in the soil–rock configuration are considered viscoelastic except for the section of the rock half-space below the level of the rock inclusion which is considered elastic. A parametric study is performed to determine controlling factors for surface displacement due to in-plane body waves. The study investigates varying the stiffness and the thickness of the rock inclusion for a range of frequencies and wave incidence angles. Anti-plane waves for this type of soil-rock configuration have been addressed in a previous article by Heymsfield (Earthquake Engng. Struct. Dyn. 28 : 841–855 (1999)). Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

12.
The mean tangential stresses at a corrugated interface between a solid, electrically insulating mantle and a liquid core of magnetic diffusivity λ are calculated for uniform rotation of both mantle and core at an angular velocity Ω in the presence of a corotating magnetic field B. The core and mantle are assumed to extend indefinitely in the horizontal plane. The interface has the form z = η(x, y), where z is the upward vertical distance and x, y are the zonal and latitudinal distances respectively. The function η(x, y) has a planetary horizontal length scale (i.e. of the order of the radius of the Earth) and small amplitude and vertical gradient. The liquid core flows with uniform mean zonal velocity U0 relative to the mantle. Ω and B possess vertical and horizontal components.The vertical (poloidal) component Bp is uniform and has a value of 5 G while the horizontal (toroidal) field BT = Bpαz, where α is a constant. When |α| ? 1, the mean horizontal stresses are found to have the same order of magnitude (10?2 N m?2) as those inferred from variations in the decade fluctuations in the length of the day, although the exact numerical values depend on the orientation of Ω as well as on the wavenumbers in the zonal and latitudinal directions.The influence of the steepness (as measured by α) of the toroidal field on the stresses is investigated to examine whether the constraint that the mean horizontal stresses at the core-mantle interface be of the order of 10?2 N m?2 might provide a selection mechanism for the behaviour of the toroidal field in the upper reaches of the outer core of the Earth. The results indicate that the restriction imposed on α is related to the value assigned to the toroidal field deep into the core. For example, if |α| ? 1 then the tangential stresses are of the right order of magnitude only if the toroidal field is comparable with the poloidal field deep in the core.  相似文献   

13.
A geochronological and structural study of dykes is used in order to define the stress field in Sardinia. The Plio-quaternary volcanism of Sardinia (between 5 and <0.2 My old) is mostly fissural and contemporaneous with an intraplate regime which reactivated normal faults of various orientations. Dykes in the Logudoro, Montiferro and Monte Arci areas (west-central and northern Sardinia) are 3.4 to 1.8 My old. Their orientations vary around a north-south mean direction which is the most frequent in the Tyrrhenian Basin. The reactivation of normal faults of various orientations may be linked to a stress regime characterised by a maximum vertical principal stressσ 1 and horizontal principal stressesσ 2 andσ 3 having close values and which can easily interchange. In these conditions, the magmatie liquid can intrude into preexisting fractures of different orientations, especially into the fractures that have an orientation close to the planeσ 1 σ 2 (with N-S orientation in Sardinia). This plane is concordant with the convergent motion of Europe and Africa plates. The influence of the pre-existing fractures on dyke orientations seems more important in this area than in southern France where there is a different tectonic regime (Féralt) andCampredon, 1983).  相似文献   

14.
Summary A numerical model is used to simulate a vertical magnetic dipole (VMD), with time harmonic excitation, buried under a segmented overburden. Results are presented in terms of readily measured field quantities. Of special interest are the ellipticities of the magnetic field in the vertical and horizontal planes. These turn out to be very sensitive to the lateral variation of the conductivity structure. Other parameters, such as the vertical magnetic field on the surface, are much less influenced. It is suggested that the combination of such measurements should lead to an optimum system for source location in application to mine rescue operations.  相似文献   

15.
The relation in which the vertical and horizontal gradients of potential field data measured along a profile across a two‐dimensional source are a Hilbert transform pair is re‐established using complex domain mathematics. In addition, a relation between the measured field and its vertical gradient in terms of a closed‐form formula is also established. The formula is based on hypersingular or Hadamard's finite‐part integral. To estimate the vertical gradient directly from the field data, Linz's algorithm of computing Hadamard's finite‐part integral is implemented. Numerical experiments are conducted on synthetically generated total magnetic intensity data with a mild level of noise contamination. A model of a magnetically polarised vertical thin sheet buried at a finite depth within a non‐magnetic half‐space was considered in generating the synthetic response. The results from numerical experiments on the mildly noise‐contaminated synthetic response are compared with those from using classical Fourier and robust regularised Hilbert transform‐based techniques.  相似文献   

16.
The amplitude of the horizontal magnetic field in the ground between two parallel wires, both carrying an alternating current in the same direction, is likely to have a saddle point if the separation between the wires is small and the frequency is low. The amplitude has a maximum in the vertical direction and a minimum in the horizontal. Rectangular geological structures in the ground which are centered between the wires have a varying effect on the magnetic fields at the surface. In general, the vertical magnetic field “crosses over” at the center of the structure. A shallow and flat lying conductor displays a broad flat type of profile when the horizontal magnetic field between the wires is measured. Changing the structure to a narrower but more conducting one at depth will provide a more pointed but still broad profile. The phase of the horizontal field is also increased. When the structure is a thin vertical dyke, the amplitude of the horizontal magnetic field anomaly due to the dyke rapidly decreases as the depth of the dyke is increased. The phase of the horizontal field is less sensitive to changes in depth of the dyke but is more sensitive to the conductivity ratio of the dyke and the half-space. The amplitude of the vertical magnetic field anomaly due to the dyke is only slightly influenced by conductivity contrast or the depth of the dyke. The phase of the vertical magnetic field, however, is strongly influenced by the conductivity contrast, particularly if the conductivity frequency product is greater than hundred. In essence, the field behaves like that of the conventional vertical loop source, but the fields are uniform over much larger areas. This suggests the possibility of using dip angle measurements for rapid reconnaissance.  相似文献   

17.
FollowingDmitriev (1960) a rigorous theoretical solution for the problem of scattering by a perfectly conducting inclined half-plane buried in a uniform conductive half-space has been obtained for plane wave excitation. The resultant integral equation for the Laplace transform of scattering current in the half-plane is solved numerically by the method of successive approximation. The scattered fields at the surface of the half-space are found by integrating the half-space Green's function over the transform of the scattering current.The effects of depth of burial and inclination, of the half-plane on the scattered fields are studied in detail. An increase in the depth of burial leads to attenuation of the fields. Inclination introduces asymmetry in the field profiles beside affecting its magnitude. Depth of exploration is greater for quadrature component. An interpretation scheme based on a phasor diagram is presented for the VLF-EM method of exploration for rich vein deposits in a conductive terrain.List of symbols x, y, z Space co-ordinates - Half-space conductivity - 0 Free-space permeability - Excitation frequency (angular) - T Time - h Depth of the half-plane - a Inclination of the half-plane - E x x-Directed total electric field - E x p x-Directed primary electric field - E xo p x-Directed primary electric field atz=0 directly over the half-plane - H y y-Component of total magnetic field - H y p y-Component of primary magnetic field - H y0 p y-Component of primary magnetic field atz=0 directly over the half-plane - H z z-Component of total magnetic field - H z p z-Component of primary magnetic field - J x Surface density ofx-directed scattering current - G Green's function - k 0,K Wave numbers - u,u 0,u 1,u 2 Functions - Space co-ordinate - s Variable in transform domain - Variable of integration - Normalized scattering current - Laplace transform of - N Normalized - , 0, 1, 2 Functions - t Variable of integration - Skin depth - H Total magnetic field - H p Primary magnetic field - H 0 p Primary magnetic field atz=0 directly over the half-plane - M,Q,R,S,U,V Functions - N 1,N 2 Functions  相似文献   

18.
On 10 April 2007, three moderate earthquakes with Mw = 4.9–5.1 occurred in the vicinity of Trichonis Lake (W. Greece). A local network composed of 12 three-component digital seismographs was installed in the epicentral area and recorded more than 1600 events. The double-difference algorithm HYPODD, incorporating both catalog and waveform cross-correlation differential travel-time data, was applied for the successful relocation of 1490 earthquakes. The latter led to the distinction of a main NW-SE trending and NE-dipping zone, as well as of three neighboring faults; a conjugate NW-SE striking and SW-dipping marginal fault mapped along the northeastern flanks of the lake; a E-W trending and south-dipping low-angle normal fault, possibly related to the major Agrinio Fault Zone (AFZ), parallel to the northern bank of the lake; a NE-SW striking and NW-dipping normal fault, likely related to a segment of the active Evinos fault, located south of the lake. Calculation of the Coulomb stress induced by the combination of the 1975 Mw = 6.0 event and the three largest events of 10 April 2007 on the inferred structures, reveals that most of the seismicity lies within the “stress-loaded” region, except for the westernmost activity, which probably belongs to the deep part of the AFZ. A total of 178 reliable focal mechanisms were determined by regional and local body-wave modeling (5 largest events) and P-wave first motion polarity data. The types of the obtained focal mechanisms are predominantly normal and strike-slip, however, numerous earthquakes were found to exhibit reverse faulting. Inversion of focal mechanism data showed that the prevailing principal horizontal component σ3 is quite homogeneous throughout the activated area with a roughly NW-SE trend, parallel to the strike of the Hellenides. On the contrary, the compressional field σ1 appears in two patterns: NE-SW trending onshore and NW-SE trending beneath the lake. This apparent rotation of σ1 by 90° reveals a complex system enclosed by the suggested NW-SE trending antithetic faults in depths between 7 and 9 km. The calculated stress ratios beneath the lake imply that vertical forces are close to the overburden pressure. The overall inferred stress pattern is rather linked to topographic variations, locally imposing increase or decrease of the vertical forces. The presence of the water in the lake possibly plays an additional important role, penetrating through the bedrock, reducing the friction coefficient, while the pore pressure and, consequently, the effective stress increase. Thus, shearing along mature fractures is enhanced, likely yielding the observed diversity.  相似文献   

19.
A half-space containing horizontally multilayered regions of different transversely isotropic elastic materials as well as a homogeneous half-space as the lowest layer is considered such that the axes of material symmetries of different layers and the lowest half-space to be as depth-wise. A rigid circular disc rested on the free surface of the whole half-space is considered to be under a forced either vertical or horizontal vibration of constant amplitudes. Because of the involved integral transforms, the mixed boundary value problems due to mixed condition at the surface of the half-space are changed to some dual integral equations, which are reduced to Fredholm integral equations of second kind. With the help of contour integration, the governing Fredholm integral equations are numerically solved. Some numerical evaluations are given for different combinations of transversely isotropic layers to show the effect of degree of anisotropy of different layers on the response of the inhomogeneous half-space.  相似文献   

20.
In this paper, we describe a non‐linear constrained inversion technique for 2D interpretation of high resolution magnetic field data along flight lines using a simple dike model. We first estimate the strike direction of a quasi 2D structure based on the eigenvector corresponding to the minimum eigenvalue of the pseudogravity gradient tensor derived from gridded, low‐pass filtered magnetic field anomalies, assuming that the magnetization direction is known. Then the measured magnetic field can be transformed into the strike coordinate system and all magnetic dike parameters – horizontal position, depth to the top, dip angle, width and susceptibility contrast – can be estimated by non‐linear least squares inversion of the high resolution magnetic field data along the flight lines. We use the Levenberg‐Marquardt algorithm together with the trust‐region‐reflective method enabling users to define inequality constraints on model parameters such that the estimated parameters are always in a trust region. Assuming that the maximum of the calculated gzz (vertical gradient of the pseudogravity field) is approximately located above the causative body, data points enclosed by a window, along the profile, centred at the maximum of gzz are used in the inversion scheme for estimating the dike parameters. The size of the window is increased until it exceeds a predefined limit. Then the solution corresponding to the minimum data fit error is chosen as the most reliable one. Using synthetic data we study the effect of random noise and interfering sources on the estimated models and we apply our method to a new aeromagnetic data set from the Särna area, west central Sweden including constraints from laboratory measurements on rock samples from the area.  相似文献   

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

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