弹性介质的Lagrange动力学与地震波方程

本文将地球介质看作是弹性介质, 从弹性体的Navier方程出发, 建立均匀弹性介质和非均匀弹性介质的分析动力学方程——Lagrange方程, 利用弹性介质的Lagrange方程导出匀弹性介质和非均匀弹性介质的地震波方程, 为用Lagrange分析动力学研究地球介质中地震波传播规律和解决地震勘探中的有关问题提供基础. 关键词： 地震勘探 弹性介质 Lagrange 方程 地震波方程     

On the theory of multiple scattering of sound waves by spherical particles in liquid and elastic media

Two independent systems of equations are derived for describing the scalar and vector potentials of the sound field in a liquid or elastic medium containing discrete inhomogeneities. One of the systems determines the sound field as the sum of the fields scattered by the particles according to the law of scattering by a single particle with the oscillation amplitudes governed by the properties of the inhomogeneous medium. The other system determines the sound field as the sum of the scattered fields formed in the inhomogeneous medium with the oscillation amplitudes of a particle in a homogeneous medium. Expressions relating the fields that occur in a medium consisting of N particles to the fields in a medium consisting of N − 1 particles are proposed. These expressions may simplify and diversify the methods used for computer simulation of sound fields with the aim to verify the calculations. The results of the study are valid for any particle concentrations under the condition that the scattering by a single particle is determined by its monopole, dipole, and rotary oscillations.     

矩孔光栅的矢量模式理论

本文引入满足均匀矢量亥姆霍兹（Ｈｅｌｍｈｏｌｔｚ）方程的矢量波函数作为基矢，对矩孔光栅的孔内外光场分别进行矢量模式展开和矢量平面波展开，并由耦合条件导出了求解展开系数的方程组，从方程组中求解出相应的振幅系数，可研究光栅的衍射场分布，该方法可研究入入射扬方向和偏振任意时的衍射效率和偏振特性等问题。     

Simulation of the vector wave field of a low-frequency sound source in the ocean: Some important features

An approach to the simulation of low frequency vector wave fields in stratified media (mainly in the ocean) is considered. The approach is characterized by an improved stability with respect to dividing the medium into many layers of arbitrary thickness. The model for the sound field of a point source is based on an integral representation of two-dimensional, cylindrically symmetric vector wave fields in inhomogeneous media, so that the contributions of all types of waves are included automatically. The model medium is subdivided into N horizontally homogeneous layers for which 4(N?1) equations are formulated to satisfy the boundary conditions between adjacent layers. The method of the generalized Schmidt matrix is used to obtain the coefficients of the equations; these coefficients are substituted into the expressions (of the Fourier-Bessel integral type) for the local parameters of the field. The latter are calculated according to the numerical procedure, and the results are used to model the distributions of the acoustic pressure and the horizontal and vertical components of the particle velocity in liquid and elastic media. The instability of the calculation procedure may result in a disagreement between the model and the exact solution. However, the disagreement is shown to occur mainly in models containing excessively thick layers. A way for improving the stability of the numerical model is suggested. The simulation results are compared with the exact analytical solution for the simplest example and with the results obtained according to the commonly used generalized matrix procedure (the benchmark problem). The examples of the practical application of the model for investigating more complex seismoacoustic wave fields in the ocean are presented.     

A full vectorial contrast source inversion scheme for three-dimensional acoustic imaging of both compressibility and density profiles

Imaging the two acoustic medium parameters density and compressibility requires the use of both the acoustic pressure and velocity wave fields, described via integral equations. Imaging is based on solving for the unknown medium parameters using known measured scattered wave fields, and it is difficult to solve this ill-posed inverse problem directly using a conjugate gradient inversion scheme. Here, a contrast source inversion method is used in which the contrast sources, defined via the product of changes in compressibility and density with the pressure and velocity wave fields, respectively, are computed iteratively. After each update of the contrast sources, an update of the medium parameters is obtained. Total variation as multiplicative regularization is used to minimize blurring in the reconstructed contrasts. The method successfully reconstructed three-dimensional contrast profiles based on changes in both density and compressibility, using synthetic data both with and without 50% white noise. The results were compared with imaging based only on the pressure wave field, where speed of sound profiles were solely based on changes in compressibility. It was found that the results improved significantly by using the full vectorial method when changes in speed of sound depended on changes in both compressibility and density.     

Polarization changes and beam profile deformation of light during acousto-optic interaction in isotropic media

A spatial Fourier transform approach is proposed to investigate the effects of polarization changes and beam profile deformation of light during acousto-optic (AO) interaction in isotropic media. The behaviour of the total scattered optical fields inside the AO cell can be properly described by a vector wave equation of which the permittivity is perturbed by an acoustic wave propagating inside the medium. In the Bragg regime, using a spatial Fourier transform approach, two coupled differential equations can be derived from the wave equation to depict AO interaction in the spatial frequency domain. Analytic solutions, which comprise the effects of changing polarization, beam deformation and propagating diffraction, can be found from the coupled equations. Detailed numerical simulations, including Fourier transforming the incident light profile to calculate the spectra of the scattered light beams and, hence, their profiles in space using the inverse transform, are presented.     

The Mie Scattering at an Inhomogeneous and Arbitrary Shaped Inclusion

We consider the scattering of an electromagnetic wave at an inhomogeneous and arbitrary shaped inclusion describing the inclusion as a non-spherical perturbation. Thus, inhomogenety and arbitrarity in the shape are treated on the same footing. The fields are developed into vector spherical harmonics. The radial parts must be calculated from a coupled set of equations using the methods known from the scattering of a scalar field at a non-spherical potential. Using an operator notation the lengthly expressions can be avoided, which contain the Clebsch-Gordan coefficients.     

A fast linearized pseudo-spherical two orders of scattering model to account for polarization in vertically inhomogeneous scattering-absorbing media

We calculate the reflection matrix for the first two orders of scattering in a vertically inhomogeneous, scattering-absorbing medium. We take full account of polarization and perform a complete linearization (analytic differentiation) of the reflection matrix with respect to both the inherent optical properties of the medium and the surface reflection condition. Further, we compute a scalar-vector correction to the total intensity due to the effect of polarization; this correction is also fully linearized. The solar beam attenuation has been computed for a pseudo-spherical atmosphere.Results from the two orders of scattering (2OS) model have been tested against scalar intensities for an inhomogeneous atmosphere, and against Stokes vector results for a homogeneous atmosphere. We have also performed backscatter simulations of reflected sunlight in the O₂A band for a variety of geometries, and compared our results with those from a full vector multiple scattering code. Our results are exact in the center of strong lines and most inaccurate in the continuum, where polarization is least significant. The s- and p-polarized radiances are always computed very accurately. The effect of gas absorption optical depth, solar zenith angle, viewing geometry, surface albedo and wind speed (in the case of ocean glint) on the intensity, polarization and corresponding weighting functions have been investigated. It is shown that the 2OS model provides fast and reliably accurate polarization corrections to the scalar-model radiance and weighting function fields. The model can be implemented in operational retrieval algorithms as an adjunct radiative transfer code to deal with polarization effects.     

Eigenwaves in an inhomogeneous azimuthally bigyrotropic medium

The problem of natural electromagnetic waves in an azimuthally magnetized bigyrotropic medium which is inhomogeneous in the radial direction is solved analytically. A system of generalized wave equations is obtained in the longitudinal components of the electrical and magnetic fields. Particular solutions of the system are represented in the form of generalized power series. The series exponents are determined and recursion formulas are built up for calculating their coefficients. The passage to particular cases of a homogeneous gyrotropic and inhomogeneous isotropic medium is realized, and the results are compared with known data from the literature.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 11, pp. 46–49, November, 1982.     

Nonstationary two-flux model of radiation transport for optical tomography of scattering media

A nonstationary two-flux model is formulated for the transport of radiation in an inhomogeneous scattering medium and is applied to the situation where such a medium is irradiated by the narrow beam of a pulsed laser. It is shown that when the time distribution of the transmitted photons is measured it is possible simultaneously to reconstruct the two spatial functions (the coefficients of absorption coefficient and of scattering of the radiation by the medium) by means of an inverse Radon transformation and the solution of a system of nonlinear differential equations on the projection lines. An analytic solution is obtained in quadratures for these differential equations. The results constitute a method of solving problems in optical tomography in an inhomogeneous scattering medium Zh. Tekh. Fiz. 67, 61–65 (May 1997)     

Phase velocities of plane waves in a pipe with a flow whose velocity varies along the radius

The phase velocities of plane waves in a pipe filled with a moving acoustic medium are studied for different laws of flow velocity variation along the pipe radius. The wave equation is solved by the discretization method, which breaks the entire pipe volume into individual cylinders under the assumption that, within each of the cylinders, the flow velocity of the medium is constant. This approach makes it possible to reduce the solution to the wave problem to solving Helmholtz equations for individual cylinders. Based on boundary conditions satisfied at the boundaries between neighboring cylinders, a homogeneous system of linear algebraic equations is obtained. From this system, with the use of the scattering matrices, a simple dispersion equation is derived for determining the phase velocities of plane waves. The stability of the numerical solution to the dispersion equation with respect to the number of cylinders is investigated. The phase velocities of quasi-homogeneous and inhomogeneous waves in a pipe are numerically calculated and analyzed for different velocities of a moving medium and different laws of flow velocity variation along the radius. It is shown that the variation that occurs in the phase velocity of a homogeneous plane wave in a pipe due to the motion of the medium is identical to the mean flow velocity for different laws of flow velocity variation along the radius. For inhomogeneous plane waves, the phase velocity increment exceeds the mean flow velocity several times and depends on both the law of wave amplitude distribution along the radius and the law of the flow velocity variation along the radius.     

On the determination of the electromagnetic field upon scattering by a small inhomogeneous spherical object

An efficient and fairly simple method of solving the problem of the incidence of a plane electromagnetic wave on an inhomogeneous object with specified spherically symmetric distributions of its electric permittivity and magnetic permeability is presented. The fields inside the object and the integrated scattering and absorption cross sections are found by assuming the object to be small compared to the vacuum wavelength. Since no constraints are imposed on the scales of the fields inside the object, the method is suitable for investigating complex cases, including those associated with the local amplification and absorption of the electromagnetic field in inhomogeneous resonant media.     

Light scattering by dielectric particles with axial symmetry

A new method is suggested for solving the problem of scattering of a plane electromagnetic wave by dielectric particles with the axial symmetry. The method is based on the separation of fields into two parts: the axially symmetric part, which is independent of the angle of rotation, and the part that is not axially symmetric, which vanishes upon averaging over this angle. The scattering problem is solved separately for each of the parts. In the first case, scalar potentials related to the azimuthal components of electromagnetic fields are used, and in the second case, superpositions of the Debye potentials and vertical components of the Hertz vectors are used. The surface integral equations for these potentials are obtained. They are represented as expansions in the spherical wave functions. The infinite systems of linear algebraic equations are obtained for unknown expansion coefficients. Our calculations demonstrated the high efficiency of the new method.     

Self-gravitating wave fields in a Gödel space

The properties of self-gravitating wave fields with integral spin (scalar and vector), compatible with a Gödel type space, are investigated. The simultaneous systems of Einstein's gravitational field equations and the equations corresponding to wave fields in Gödel's metric are solved. For the scalar field, the solutions are obtained for different types of interaction Lagrangians for the gravitational and scalar fields. It is shown that for a massive vector field the relations obtained between the constants lead, within the scope of the strong gravitation theory, to the classical expression for the spin of elementary particles.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 10, pp. 59–63, October, 1981.The authors are grateful to the participants of the theoretical seminar conducted by D. Ivanenko for discussing the results of this work.     

On Inhomogeneous Polarized States near the Phase Transition Point in a Thin Ferroelectric Film

It is shown in terms of the phenomenological Landau theory of phase transitions that a phase transition to an inhomogeneous polar phase preceding in temperature a phase transition to a homogeneous polar state is possible. As a result of solving a boundary eigenvalue problem for the polarization equilibrium equation and electrostatics equations, wave vector k_⊥ characterizing the inhomogeneous phase has been determined and the temperature boundaries of its existence in the dependence on the film thickness and its surface properties have been found.     

Time-Domain Analysis of Millimeter Wave Circulation Around a Magnetised Ferrite Sphere

An efficient numerical method has been devised for the study of wave circulated by a magnetised ferrite sphere. It is a finite-difference time-domain formulation that incorporates Mur's absorbing boundary conditions and a perfectly matched layer. The electromagnetic fields inside the ferrite body are calculated using special updating equations derived from the equation of motion of the magnetization vector and Maxwell's curl equations in consistency. The electromagnetic fields inside ferrite and the power-density distribution on the ferrite's surface normal to the bias external magnetic field are obtained in a wide frequency band with a single time domain run. It is observed that an incident plane wave would circulate around the magnetised ferrite body in an open space as if the ferrite were placed inside a waveguide / microstrip junction circulators.     

The Kirchhoff gauge

We discuss the Kirchhoff gauge in classical electrodynamics. In this gauge, the scalar potential satisfies an elliptical equation and the vector potential satisfies a wave equation with a nonlocal source. We find the solutions of both equations and show that, despite of the unphysical character of the scalar potential, the electric and magnetic fields obtained from the scalar and vector potentials are given by their well-known retarded expressions. We note that the Kirchhoff gauge pertains to the class of gauges known as the velocity gauge.     

Sound propagation and scattering in media with random inhomogeneities of sound speed, density and medium velocity

This review presents both classical and new results of the theory of sound propagation in media with random inhomogeneities of sound speed, density and medium velocity (mainly in the atmosphere and ocean). An equation for a sound wave in a moving inhomogeneous medium is presented, which has a wider range of applicability than those used before. Starting from this equation, the statistical characteristics of the sound field in a moving random medium are calculated using Born-approximation, ray, Rytov and parabolic-equation methods, and the theory of multiple scattering. The results obtained show, in particular, that certain equations previously widely used in the theory of sound propagation in moving random media must now be revised. The theory presented can be used not only to calculate the statistical characteristics of sound waves in the turbulent atmosphere or ocean but also to solve inverse problems and develop new remote-sensing methods. A number of practical problems of sound propagation in moving random media are listed and the further development of this field of acoustics is considered.     

Exact Solutions of Covariant Wave Equations with a Multipole Source Term on Curved Spacetimes

A formalism is presented for calculating exactsolutions of covariant inhomogeneous scalar and tensorwave equations whose source terms are arbitrary ordermultipoles on a curved background spacetime. The developed formalism is based on the theory ofthe higher-order fundamental solutions for wave equationwhich are the distributions that satisfy theinhomogeneous wave equation with the corresponding order covariant derivatives of the Dirac deltafunction on the right-hand side. Like the classicalGreen's function for a scalar wave equation, thehigher-order fundamental solutions contain a direct termwhich has support on the light cone as well as a tailterm which has support inside the light cone. Knowinghow to compute the fundamental solutions of arbitraryorder, one can find exact multipole solutions of wave equations on curved spacetimes. Wepresent complete recurrent algorithms for calculatingthe arbitrary-order fundamental solutions and the exactmultipole solutions in a form convenient for practical computations. As an example we apply thealgorithm to a massless scalar wave field on aparticular Robertson-Walker spacetime.     

On multidimensional ultrasonic scattering in an inhomogeneous elastic background.

This work is concerned with the modeling of elastic wave scattering by solid or fluid-filled objects embedded in an inhomogeneous elastic background. The medium is probed by a monochromatic force and the scattered field is computed (forward problem) or observed (inverse problem) at some known receiver locations. Based on vector integral equations for elastic scattering, a general framework is developed, independent of both the problem geometry and the transmitter-receiver characteristics. This framework encompasses both forward and inverse modeling. In the forward model, a Born approximation for an inhomogeneous background is applied to obtain a closed form expression for the scattered field. In the inverse model, this approximation is also invoked to linearize for the multiparameter characteristic of the object. Finally, an iterative inversion scheme alternating forward and inverse modeling is proposed to improve the resolution and accuracy of the reconstruction algorithm.