共查询到20条相似文献，搜索用时 687 毫秒
1.
Summary. The forward solution of the general two-dimensional problem of induction in a model earth comprising a uniformly conducting half-space covered by a thin sheet of variable integrated conductivity is obtained. Unlike some previous treatments of similar problems, the method presented here does not require the field to be separated into its normal and anomalous parts. Both the E - and B -polarization modes of induction are considered and in each case the solution is expressed in terms of the horizontal component of the electric field satisfying, on the surface of the conductor, a singular integral equation whose kernel is a well-known analytic function. A recently published solution of the coast effect is included as a special case. The numerical procedure for solving the integral equations is described and some illustrative calculations are presented. 相似文献
2.
S. Méndez-Delgado E. Gómez-Treviño & M. A. Pérez-Flores 《Geophysical Journal International》1999,137(2):336-352
We present a semi-analytical, unifying approach for modelling the electromagnetic response of 3-D bodies excited by low-frequency electric and magnetic sources. We write the electric and magnetic fields in terms of power series of angular frequency, and show that to obey Maxwell's equations, the fields must be real when the exponent is even, and imaginary when it is odd. This leads to the result that the scattering equations for direct current fields and for fields proportional to frequency can both be explicitly formulated using a single, real dyadic Green's function. Although the underground current flow in each case is due to different physical phenomena, the interaction of the scattering currents is of the same type in both cases. This implies that direct current resistivity, magnetometric resistivity and electric and magnetic measurements at low induction numbers can all be modelled in parallel using basically the same algorithm. We make a systematic derivation of the quantities required and show that for these cases they can all be expressed analytically. The problem is finally formulated as the solution of a system of linear equations. The matrix of the system is real and does not depend on the type of source or receiver. We present modelling results for different arrays and apply the algorithm to the interpretation of field data. We assume the standard dipoledipole resistivity array for the direct current case, and vertical and horizontal magnetic dipoles for induction measurements. In the case of magnetometric resistivity we introduce a moving array composed of an electric dipole and a directional magnetometer. The array has multiple separations for depth discrimination and can operate in two modes. The mode where the predominant current flow runs along the profile is called MMR-TM. This mode is more sensitive to lateral variations in resistivity than its counterpart, MMR-TE, where the mode of conduction is predominantly perpendicular to the profile. 相似文献
3.
T. W. Dawson 《Geophysical Journal International》1983,73(1):83-107
Summary. An analytical solution is obtained for the E-polarization problem of electromagnetic induction in two adjacent half-sheets underlain by a uniform conducting half-space. In this mode the inducing magnetic field is assumed horizontal, uniform and perpendicular to the discontinuity. The same model was previously solved under B-polarization by Dawson & Weaver. The present solution then completes the study of two-dimensional induction in the described model. Further, it extends both the analytic E-polarization solution of Weidelt by the inclusion of an underlying conductor and that of Raval, Weaver & Dawson by the inclusion of arbitrary conductance values for the two surface sheets. The solution may be used as an idealized model of the coast effect and allows detailed study of the field behaviour near the discontinuity. The horizontal magnetic field on each side of the surface layer has a finite jump discontinuity at the interface and the vertical magnetic field exhibits a logarithmic singularity there. If the right-hand conductance (say) becomes infinite, the horizontal magnetic field exhibits an algebraic singularity as the coastline is approached from the right, while the vertical magnetic field does likewise from the left. Calculations are presented for the same two models as discussed in B-polarization by Dawson & Weaver and the results are compared to values obtained from a more general numerical scheme. The electric current distribution inside the conducting half-space is depicted for the second model. 相似文献
4.
C. Sozou 《Geophysical Journal International》1997,130(1):151-156
Laplace's tidal equations for the case of an ocean of constant depth bounded by meridians were considered by two authors at a specific frequency as an eigenvalue problem in the azimuthal wavenumber. A finite spectrum of eigenwavenumbers was found. That eigenvalue problem is re-examined by means of asymptotic techniques and numerical integration of the governing equation of the problem. At low frequencies a formula connecting the frequency and the number of eigensolutions is established. It is shown that at a given frequency the spectrum of eigenwavenumbers is wider than that reported, but (for this type of solution) the meridional boundary conditions are satisfied approximately only for the case of very low frequencies. 相似文献
5.
An analytical approach for the scattering of SH waves by a symmetrical V-shaped canyon: shallow case 总被引：1，自引：0，他引：1
Based on the application of the region-matching technique, an analytical approach is presented for the scattering of plane SH waves from a shallow symmetrical V-shaped canyon, and then a series solution is derived. The analysed region is divided into an enclosed and an open region by introducing a semi-circular auxiliary boundary. In each region, the displacement field can be expressed as infinite sum of appropriate wavefunctions satisfying partial boundary conditions, respectively. The unknown coefficients can be determined by enforcing the continuity conditions in connection with the Graf's addition formula. The frequency- and time-domain responses are both evaluated and displayed for several physical parameters. From graphical results, the effects of the canyon depth on surface ground motion are conspicuous. The proposed series solutions can serve as benchmark for numerical methods, in particular for those at much higher frequencies. 相似文献
6.
D. W. Vasco 《Geophysical Journal International》2000,142(3):970-990
Many geophysical inverse problems derive from governing partial differential equations with unknown coefficients. Alternatively, inverse problems often arise from integral equations associated with a Green's function solution to a governing differential equation. In their discrete form such equations reduce to systems of polynomial equations, known as algebraic equations. Using techniques from computational algebra one can address questions of the existence of solutions to such equations as well as the uniqueness of the solutions. The techniques are enumerative and exhaustive, requiring a finite number of computer operations. For example, calculating a bound to the total number of solutions reduces to computing the dimension of a linear vector space. The solution set itself may be constructed through the solution of an eigenvalue problem. The techniques are applied to a set of synthetic magnetotelluric values generated by conductivity variations within a layer. We find that the estimation of the conductivity and the electric field in the subsurface, based upon single-frequency magnetotelluric field values, is equivalent to a linear inverse problem. The techniques are also illustrated by an application to a magnetotelluric data set gathered at Battle Mountain, Nevada. Surface observations of the electric ( E y ) and magnetic ( H x ) fields are used to construct a model of subsurface electrical structure. Using techniques for algebraic equations it is shown that solutions exist, and that the set of solutions is finite. The total number of solutions is bounded above at 134 217 728. A numerical solution of the algebraic equations generates a conductivity structure in accordance with the current geological model for the area. 相似文献
7.
Kiyoshi Yomogida 《Geophysical Journal International》1985,82(3):511-533
Summary. Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. Using ray-centred coordinates, parabolic equations are obtained for lateral variations while vertical structural variations at a given point are specified by eigenfunctions of normal mode theory as for the laterally homogeneous case. Final results on wavefields close to a ray can be expressed by formulations similar to those for elastic body waves in 2-D laterally heterogeneous media, except that the vertical dependence is described by eigenfunctions of 'local' Love or Rayleigh waves. The transport equation is written in terms of geometrical-ray spreading, group velocity and an energy integral. For the horizontal components there are both principal and additional components to describe the curvature of rays along the surface, as in the case of elastic body waves. The vertical component is decoupled from the horizontal components. With complex parameters the solutions for the dynamic ray tracing system correspond to Gaussian beams: the amplitude distribution is bell-shaped along the direction perpendicular to the ray and the solution is regular everywhere, even at caustics. Most of the characteristics of Gaussian beams for 2-D elastic body waves are also applicable to the surface wave case. At each frequency the solution may be regarded as a set of eigenfunctions propagating over a 2-D surface according to the phase velocity mapping. 相似文献
8.
Summary. The propagation of a pulsed elastic wave in the following geometry is considered. An elastic half-space has a surface layer of a different material and the layer furthermore contains a bounded 3-D inhomogeneity. The exciting source is an explosion, modelled as an isotropic pressure point source with Gaussian behaviour in time.
The time-harmonic problem is solved using the null field approach (the T matrix method), and a frequency integral then gives the time-domain response. The main tools of the null field approach are integral representations containing the free space Green's dyadic, expansions in plane and spherical vector wave functions, and transformations between plane and spherical vector wave functions. It should be noted that the null field approach gives the solution to the full elastodynamic equations with, in principle, an arbitrarily high accuracy. Thus no ray approximations or the like are used. The main numerical limitation is that only low and intermediate frequencies, in the sense that the diameter of the inhomogeneity can only be a few wavelengths, can be considered.
The numerical examples show synthetic seismograms consisting of data from 15 observation points at increasing distances from the source. The normal component of the velocity field is computed and the anomalous field due to the inhomogeneity is sometimes shown separately. The shape of the inhomogeneity, the location and depth of the source, and the material parameters are all varied to illustrate the relative importance of the various parameters. Several specific wave types can be identified in the seismograms: Rayleigh waves, direct and reflected P -waves, and head waves. 相似文献
The time-harmonic problem is solved using the null field approach (the T matrix method), and a frequency integral then gives the time-domain response. The main tools of the null field approach are integral representations containing the free space Green's dyadic, expansions in plane and spherical vector wave functions, and transformations between plane and spherical vector wave functions. It should be noted that the null field approach gives the solution to the full elastodynamic equations with, in principle, an arbitrarily high accuracy. Thus no ray approximations or the like are used. The main numerical limitation is that only low and intermediate frequencies, in the sense that the diameter of the inhomogeneity can only be a few wavelengths, can be considered.
The numerical examples show synthetic seismograms consisting of data from 15 observation points at increasing distances from the source. The normal component of the velocity field is computed and the anomalous field due to the inhomogeneity is sometimes shown separately. The shape of the inhomogeneity, the location and depth of the source, and the material parameters are all varied to illustrate the relative importance of the various parameters. Several specific wave types can be identified in the seismograms: Rayleigh waves, direct and reflected P -waves, and head waves. 相似文献
9.
Summary. Using a single scattering approximation, we derive equations for the scattering attenuation coefficients of P- and S -body waves. We discuss our results in the light of some recent energy renormalization approaches to seismic wave scattering. Practical methods for calculating the scattering attenuation coefficients for various earth models are emphasized. The conversions of P - to S -waves and S- to P -waves are included in the theory. The earth models are assumed to be randomly inhomogeneous, with their properties known only through their average wavenumber power spectra. We approximate the power spectra with piecewise constant functions, each segment of which contributes to the net, frequency-dependent, scattering attenuation coefficient. The smallest and largest wavenumbers of a segment can be plotted along with the wavevectors of the incident and scattered waves on a wavenumber diagram. This diagram gives a geometric interpretation for the frequency behaviour associated with each spectral segment, including a 'transition' peak that is due entirely to the wavenumber limits of the segment. For regions of the earth where the inhomogeneity spectra are concentrated in a band of wavenumbers, it should be possible to observed such a peak in the apparent attenuation of seismic waves. We give both the frequency and distance limits on the accuracy of the theoretical results. 相似文献
10.
Robert I. Odom 《Geophysical Journal International》1986,86(2):425-453
Summary. A coupled mode theory is used to examine surface wave propagation in a laterally inhomogeneous acoustic waveguide. The theory is developed from the equations of motion for the pressure and velocity fields. The presence of lateral inhomogeneities in the form of varying layer thickness causes coupling among the discrete modes of the waveguide and radiation to the continuum. Expressions for the coupling coefficients among all mode types including coupling to the continuum spectrum are derived. The coupling coefficients are proportional to the horizontal derivative of the function describing the interface between layers of constant material properties but varying thickness. The coupled mode equations are solved in approximation for the case of a sinusoidal boundary and a sloping boundary. The results for radiation losses due to interaction with the irregular boundary of the waveguide are presented in analytical form, which clearly show the primary physical effects on the wavefield of the interaction. The far field amplitude of the scattered modes, excited by the interaction of some incident signal with a weak boundary irregularity, is modulated by the spatial Fourier transform of the irregularity. 相似文献
11.
12.
Summary. A layer of constant thickness over a half-space is assumed, and the propagation of head waves is considered for the following two cases: (1) the P -wave velocity varies in the layer in the horizontal direction, and is constant in the half-space: (2) the P -wave velocity varies in the half-space in the horizontal direction, and is constant in the layer. In each case the horizontal velocity gradient is assumed to remain constant. The wave propagation is investigated in the direction of the gradient (direct profile), and opposite to it (reverse profile). Formulae for the travel times and the amplitudes are obtained on the basis of ray-theoretical considerations. Conditions are discussed for the discrimination in a field experiment between the case of a sloping boundary separating the homogeneous media, and the case of an intrinsic horizontal velocity gradient. 相似文献
13.
Summary. A new method for solving problems in three-dimensional electromagnetic induction in which the Earth is represented by a uniformly conducting half-space overlain by a surface layer of variable conductance is presented. Unlike previous treatments of this type of problem the method does not require the fields to be separated into their normal and anomalous parts, nor is it necessary to assume that the anomalous region is surrounded by a uniform structure; the model may approach either an E- or a B -polarization configuration at infinity. The solution is expressed as a vector integral equation in the horizontal electric field at the surface. The kernel of the integral is a Green's tensor which is expressed in terms of elementary functions that are independent of the conductance. The method is applied to an illustrative model representing an island near a bent coastline which extends to infinity in perpendicular directions. 相似文献
14.
Steven N. Ward 《Geophysical Journal International》1981,66(1):23-30
Summary. Propagator matrix solutions to the elastic equations of motion in spherically symmetric, inhomogeneous media with moment tensor sources are recast into a simple and intuitively satisfying form which is applicable to both exact and approximate calculations. The transformed expression benefits from the analogous equations of normal mode excitation, while clearly distinguishing the finer partitions of the displacement field and the more flexible boundary conditions that body wave formulations provide. I believe that this new representation, because of its many advantages, should be favoured as the foundation for elastic wave calculations in a sphere. 相似文献
15.
N. S. Heaps 《Geophysical Journal International》1981,64(1):291-302
Summary. The mathematical basis of a spectral method for the numerical solution of the three-dimensional hydrodynamic equations for tides and surges is described. Vertical eddy viscosity is prescribed in two layers, upper and lower. In each layer, horizontal components of current are expanded through the vertical in terms of a set of eigenfunctions. Coefficients of these expansions are evaluated in the horizontal and through time using a two-dimensional numerical time-stepping procedure. Thence the three-dimensional current structure is determined. 相似文献
16.
A conservative staggered-grid finite difference method is presented for computing the electromagnetic induction response of an arbitrary heterogeneous conducting sphere by external current excitation. This method is appropriate as the forward solution for the problem of determining the electrical conductivity of the Earth's deep interior. This solution in spherical geometry is derived from that originally presented by Mackie et al. (1994 ) for Cartesian geometry. The difference equations that we solve are second order in the magnetic field H , and are derived from the integral form of Maxwell's equations on a staggered grid in spherical coordinates. The resulting matrix system of equations is sparse, symmetric, real everywhere except along the diagonal and ill-conditioned. The system is solved using the minimum residual conjugate gradient method with preconditioning by incomplete Cholesky decomposition of the diagonal sub-blocks of the coefficient matrix. In order to ensure there is zero H divergence in the solution, corrections are made to the H field every few iterations. In order to validate the code, we compare our results against an integral equation solution for an azimuthally symmetric, buried thin spherical shell model ( Kuvshinov & Pankratov 1994 ), and against a quasi-analytic solution for an azimuthally asymmetric configuration of eccentrically nested spheres ( Martinec 1998 ). 相似文献
17.
Elastic scattering from a continuous and laterally unbounded heterogeneous layer has been formulated using the Born approximation. A general solution of the scattered wave equation for the above-stated medium has been given in terms of a Fourier integral over plane waves. Far-field asymptotic expressions for weak elastic scattering by a finite, continuous and inhomogeneous layer have been presented which agree with earlier results. For perturbations of the two elastic parameters and the density having the same form of spatial variation, the spectrum of plane waves scattered from a heterogeneous layer is expressed as a product of an 'elastic scattering factor'and a 'distribution factor'. As in earlier results for small-scale heterogeneity, the scattering pattern depends on various combinations of perturbations of elastic parameters and density. In order to show the general characteristics of the elastic wave scattering, some scattering patterns have been given. 相似文献
18.
When interpreting electromagnetic fields observed at the Earth's surface in a realistic geophysical environment it is often necessary to pay special attention to the effects caused by inhomogeneities of the subsurface sedimentary and/or water layer and by inhomogeneities of the Earth's crust. The inhomogeneities of the Earth's crust are expected to be especially important when the electromagnetic field is generated by a source located in a magma chamber of a volcano. The simulation of such effects can be carried out using generalized thin-sheet models, which were independently introduced by Dmitriev (1969 ) and Ranganayaki & Madden (1980 ). In the first part of the paper, a system of integral equations is derived for the horizontal current that flows in the subsurface inhomogeneous conductive layer and for the vertical current crossing the inhomogeneous resistive layer representing the Earth's mantle. The terms relating to the finite thickness of the laterally inhomogeneous part of the model are retained in the equations. This only marginally complicates the equations, whilst allowing for a significant expansion of the approximation limits.
The system of integral equations is solved using the iterative dissipative method developed by the authors in the period from 1978 to 1988. The method can be applied to the simulation of the electromagnetic field in an arbitrary inhomogeneous medium that dissipates the electromagnetic energy. When considered on a finite numerical grid, the integral equations are reduced to a system of linear equations that possess the same contraction properties as the original equations. As a result, the rate at which the iterative-perturbation sequence converges to the solution remains independent of the numerical grid used for the calculations. In contrast to previous publications on the method, aspects of the algorithm implementation that guarantee its effectiveness and robustness are discussed here. 相似文献
The system of integral equations is solved using the iterative dissipative method developed by the authors in the period from 1978 to 1988. The method can be applied to the simulation of the electromagnetic field in an arbitrary inhomogeneous medium that dissipates the electromagnetic energy. When considered on a finite numerical grid, the integral equations are reduced to a system of linear equations that possess the same contraction properties as the original equations. As a result, the rate at which the iterative-perturbation sequence converges to the solution remains independent of the numerical grid used for the calculations. In contrast to previous publications on the method, aspects of the algorithm implementation that guarantee its effectiveness and robustness are discussed here. 相似文献
19.
Ghislain R. Franssens 《Geophysical Journal International》1983,75(3):669-691
Summary. The calculation of the two-dimensional elasto-dynamic Green's function for a stratified medium is investigated. The solution is represented in the form of an inverse Fourier integral which is to be integrated along a properly chosen path in the complex wavenumber plane. The integrand is computed using a modified propagator matrix method.
This method is based on a mixed formulation using the propagator matrix and the matrix of minors of the propagator matrix (compound matrix). The major advantages of this approach are the elimination of the numerical loss of precision problems associated with the Thomson-Haskell formulation, without losing the attractive tractability and compactness of the propagator matrix method. This modified method is first mathematically derived, and theoretical seismograms are then presented for two examples. 相似文献
This method is based on a mixed formulation using the propagator matrix and the matrix of minors of the propagator matrix (compound matrix). The major advantages of this approach are the elimination of the numerical loss of precision problems associated with the Thomson-Haskell formulation, without losing the attractive tractability and compactness of the propagator matrix method. This modified method is first mathematically derived, and theoretical seismograms are then presented for two examples. 相似文献
20.
A series solution of the plane SH-waves incident on a partially filled semi-circular alluvial valley imbedded in a half-space is presented. Based on the region-matching method, the analysed region is decomposed into two subregions by the interface between two media. The antiplane displacement field of each subregion is expressed in terms of an infinite series of cylindrical wavefunctions with unknown expansion coefficients. After imposing the traction-free condition on the curved valley surface and the matching conditions on the interface with the aid of Graf's addition theorem, the unknown coefficients are obtained. Both the frequency- and time-domain responses are evaluated. In the theoretical derivation of this work, two classical exact series solutions are also included, so the present series solution is more general than those given before. Visible effects of different physical parameters on ground surface motions are illustrated in graphical form. 相似文献