首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 156 毫秒
1.
The continuation inverse problem revisited   总被引：1，自引：0，他引：1  
The non-uniqueness of the continuation of a finite collection of harmonic potential field data to a level surface in the source-free region forces its treatment as an inverse problem. A formalism is proposed for the construction of continuation functions which are extremal by various measures. The problem is cast in such a form that the inverse problem solution is the potential function on the lowest horizontal surface above all sources, serving as the boundary function for the Dirichlet problem in the upper half-plane. The desired continuation, at the higher level of interest, must then be in the range of the upward continuation operator acting on this boundary function, rather than being allowed the full freedom of itself being part of a Dirichlet problem boundary function. Extremal solutions minimize non-linear functionals of the continuation function, which are re-expressed as different functionals of the boundary function. A crux of the method is that there is no essential distinction between the upward and downward continuation inverse problems to levels above or below data locations. Casting the optimization as a Lagrange multiplier problem leads to an integral equation for the boundary function, which is readily solved in the Fourier domain for a certain class of functionals. The desired extremal continuation is then given by upward continuation. It is found that for some functionals, application of the Lagrange multiplier theorem requires a further restriction on the set of allowable boundary functions: bandlimitedness is a natural choice for the continuation problem. With this imposition, the theory is developed in detail for semi-norm functionals penalizing departure from a constant potential, in the 2-norm and Sobelev norm senses, and illustrated by application for a small synthetic Deep Tow magnetic field data set.  相似文献   

2.
We use the Direct Solution Method (DSM) together with the modified operators derived by Geller & Takeuchi (1995) and Takeuchi, Geller & Cummins (1996) to compute complete synthetic seismograms and their partial derivatives for laterally heterogeneous models in spherical coordinates. The methods presented in this paper are well suited to conducting waveform inversion for 3-D Earth structure. No assumptions of weak perturbation are necessary, although such approximations greatly improve computational efficiency when their use is appropriate.
An example calculation is presented in which the toroidal wavefield is calculated for an axisymmetric model for which velocity is dependent on depth and latitude but not longitude. The wavefield calculated using the DSM agrees well with wavefronts calculated by tracing rays. To demonstrate that our algorithm is not limited to weak, aspherical perturbations to a spherically symmetric structure, we consider a model for which the latitude-dependent part of the velocity structure is very strong.  相似文献   

3.
Summary. The Lanczos method of separating exponentials is applied to the Fourier transform of seismograms in order to separate the various modes that contribute to the given portion of the seismograms. Phase velocities and amplitudes are obtained as functions of the frequency. When applying the method to artificial seismograms, which are built as an exact superposition of a number of modes, the separation is very accurate. The method was also applied to the surface-wave portion of numerical seismograms for a vertical point force in a layered medium. The phase velocity and amplitude of the fundamental mode are obtained. These functions were taken as the first guess in the Backus—Gilbert generalized inverse procedure and the process converged very rapidly. When a perturbation of the phases and amplitudes is taken as the first guess the process converges to the true model when enough data are available.  相似文献   

4.
We report results from the Seismic Wide-Angle and Broadband Survey carried out over the Mid North Sea High. This paper focuses on integrating the information from a conventional deep multichannel reflection profile and a coincident wide-angle profile obtained by recording the same shots on a set of ocean bottom hydrophones (OBH). To achieve this integration, a new traveltime inversion scheme was developed (reported elsewhere) that was used to invert traveltime information from both the wide-angle OBH records and the reflection profile simultaneously. Results from the inversion were evaluated by producing synthetic seismograms from the final inversion model and comparing them with the observed wide-angle data, and an excellent match was obtained. It was possible to fine-tune velocities in less well-resolved parts of the model by considering the critical distance for the Moho reflection. The seismic velocity model was checked for compatibility with the gravity field, and used to migrate and depth-convert the reflection profile. The unreflective upper crust is characterized by a high velocity gradient, whilst the highly reflective lower crust is associated with a low velocity gradient. At the base of the crust there are several subhorizontal reflectors, a few kilometres apart in depth, and correlatable laterally for several tens of kilometres. These reflectors are interpreted as representing a strike section through northward-dipping reflectors at the base of the crust, identified on orthogonal profiles by Freeman et al. (1988) as being slivers of subducted and imbricated oceanic crust, relics of the mid-Palaeozoic Iapetus Ocean.  相似文献   

5.
  相似文献   

6.
Summary. A general theory of inversion of the reflected travel-time curve is developed. The properties of the travel-time curve and its analytical continuation are examined. The inversion formulae are presented and the structure of the set of velocity—depth functions is studied in detail. The inverse problem or a finite number of reflectors and the problem of separation of the reflected and multiple travel-time curves are both considered.  相似文献   

7.
A general one-way representation of seismic data can be obtained by substituting a Green's one-way wavefield matrix into a reciprocity theorem of the convolution type for one-way wavefields. From this general one-way representation, several special cases can be derived.
By introducing a Green's one-way wavefield matrix for primaries , a generalized Bremmer series representation is obtained. Terminating this series after the first-order term yields a primary representation of seismic reflection data. According to this representation, primary seismic reflection data are proportional to a reflection operator, 'modified' by primary propagators for downgoing and upgoing waves. For seismic imaging, these propagators need to be inverted. Stable inverse primary propagators can easily be obtained from a one-way reciprocity theorem of the correlation type.
By introducing a Green's one-way wavefield matrix for generalized primaries , an alternative representation is obtained in which multiple scattering is organized quite differently (in comparison with the generalized Bremmer series representation). According to the generalized primary representation, full seismic reflection data are proportional to a reflection operator, 'modified' by generalized primary propagators for downgoing and upgoing waves. Internal multiple scattering is fully included in the generalized primary propagators {either via a series expansion or in a parametrized way). Stable inverse generalized primary propagators can be obtained from the one-way reciprocity theorem of the correlation type. These inverse propagators are the nucleus for seismic imaging techniques that take the angle-dependent dispersion effects due to fine-layering into account.  相似文献   

8.
Upward and downward continuation as inverse problems   总被引：2，自引：0，他引：2  
Summary . The formalism of Backus & Gilbert is applied to the problems of upward and downward continuation of harmonic functions. We first treat downward continuation of a two-dimensional field to a level surface everywhere below the observation locations; the calculation of resolving widths and solution estimates is a straightforward application of Backus—Gilbert theory. The extension to the downward continuation of a three-dimensional field uses a delta criterion giving resolving areas rather than widths. A feature not encountered in conventional Backus—Gilbert problems is the requirement of an additional constraint to guarantee the existence of the resolution integrals. Finally, we consider upward continuation of a two-dimensional field to a level above all observations. We find that solution estimates must be weighted averages of the field not only on this level, but also on a line passing between the observations and sources. Weighting on the lower line may be traded off against resolution on the upper level.  相似文献   

9.
Reciprocity theorems for one-way wavefields   总被引：1，自引：0，他引：1  
Acoustic reciprocity theorems have proved their usefulness in the study of forward and inverse scattering problems. The reciprocity theorems in the literature apply to the two-way (i.e. total) wavefield, and are thus not compatible with one-way wave theory, which is often applied in seismic exploration. By transforming the two-way wave equation into a coupled system of one-way wave equations for downgoing and upgoing waves it appears to be possible to derive 'one-way reciprocity theorems" along the same lines as the usual derivation of the 'two-way reciprocity theorems'. However, for the one-way reciprocity theorems it is not directly obvious that the 'contrast term' vanishes when the medium parameters in the two different states are identical. By introducing a modal expansion of the Helraholtz operator, its square root can be derived, which appears to have a symmetric kernel. This symmetry property appears to be sufficient to let the contrast term vanish in the above-mentioned situation.
The one-way reciprocity theorem of the convolution type is exact, whereas the one-way reciprocity theorem of the correlation type ignores evanescent wave modes. The extension to the elastodynamic situation is not trivial, but it can be shown relatively easily that similar reciprocity theorems apply if the (non-unique) decomposition of the elastodynamic two-way operator is done in such a way that the elastodynamic one-way operators satisfy similar symmetry properties to the acoustic one-way operators.  相似文献   

10.
Summary. In the present paper, we give an approximate analytic solution of the Dirichlet and Neumann problems in the two-dimensional case for a contour of any shape. The formalism is based on a representation of potential functions as a sum of elementary interpolating functions and uses the theory of generalized inverse matrices. Formulae are given in the cases of Cartesian, polar and elliptic co-ordinates. The fact that an analytic expression for continuation is obtained enables one to compute and draw equipotential lines as well as field lines; besides, any further computations that one might want to perform on either the field or the potential can be handled in an analytic way. The formalism can be extended from the Laplace to the Helmholtz equation.
We give examples in the case of magnetostatics, treating in more detail the problem of the distortion of magnetic field lines by an inclusion. We also show how the method allows the computation of conformal mappings in otherwise intricate situations.  相似文献   

11.
Wave propagation is studied in a general anisotropic poroelastic solid. The presence of dissipation due to fluid-viscosity as well as hydraulic anisotropy of pore permeability are also considered. Biot's theory is used to derive a system of modified Christoffel equations for the propagation of plane harmonic waves in porous media. A non-trivial solution of this system is ensured by a determinantal equation. This equation is separated into two different polynomial equations. One is the quartic equation whose roots represent the complex velocities of four attenuating waves in the medium. The other is a eighth-degree polynomial whose roots represent the vertical slowness values for the four waves propagating upward and downward in a finite porous medium. Procedure is explained to associate the numerically obtained roots with the waves propagating in the medium. The slowness surfaces of waves reflected at the boundary of the medium are computed for a realistic numerical model. The behaviours of phase velocity surfaces are analysed with the help of numerical examples.  相似文献   

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

13.
Summary. A new method is presented for the direct inversion of seismic refraction data in dipping planar structure. Three recording geometries, each consisting of two common-shot profiles, are considered: reversed, split, and roll-along profiles. Inversion is achieved via slant stacking the common-shot wavefield to obtain a delay time–slowness (tau– p ) wavefield. The tau– p curves from two shotpoints describing the critical raypath of refracted and post-critically reflected arrivals are automatically picked using coherency measurements and the two curves are jointly used to calculate velocity and dip of isovelocity lines iteratively, thereby obtaining the final two-dimensional velocity model.
This procedure has been successfully applied to synthetic seismograms calculated for a dipping structure and to field data from central California. The results indicate that direct inversion of closely-spaced refraction/wide-aperture reflection data can practically be achieved in laterally inhomogeneous structures.  相似文献   

14.
Velocity estimation remains one of the main problems when imaging the subsurface with seismic reflection data. Traveltime inversion enables us to obtain large-scale structures of the velocity field and the position of seismic reflectors. However, as the media currently under study are becoming more and more complex, we need to know the finer-scale structures. The problem is that below a certain range of velocity heterogeneities, deterministic methods become difficult to use, so we turn to a probabilistic approach. With this in view, we characterize the velocity field as a random field defined by its first and second statistical moments. Usually, a seismic random medium is defined as a homogeneous velocity background perturbed by a small random field that is assumed to be stationary. Thus, we make a link between such a random velocity medium (together with a simple reflector) and seismic reflection traveltimes. Assuming that the traveltimes are ergodic, we use 2-D seismic reflection geometry to study the decrease in the statistical traveltime fluctuations as a function of the offset (the source–receiver distance). Our formulae are based on the Rytov approximation and the parabolic approximation for acoustic waves. The validity and the limits are established for both of these approximations in statistically anisotropic random media. Finally, theoretical inversion procedures are developed for the horizontal correlation structure of the velocity heterogeneities for the simplest case of a horizontal reflector. Synthetic seismograms are then computed (on particular realizations of random media) by simulating scalar wave propagation via finite difference algorithms. There is good agreement between the theoretical and experimental results.  相似文献   

15.
Time-domain electromagnetic migration in the solution of inverse problems   总被引：5，自引：0，他引：5  
Time-domain electromagnetic (TDEM) migration is based on downward extrapolation of the observed field in reverse time. In fact, the migrated EM field is the solution of the boundary-value problem for the adjoint Maxwell's equations. The important question is how this imaging technique can be related to the solution of the geoelectrical inverse problem. In this paper we introduce a new formulation of the inverse problem, based on the minimization of the residual-field energy flow through the surface or profile of observations. We demonstrate that TDEM migration can be interpreted as the first step in the solution of this specially formulated TDEM inverse problem. However, in many practical situations this first step produces a very efficient approximation to the geoelectrical model, which makes electromagnetic migration so attractive for practical applications. We demonstrate the effectiveness of this approach in inverting synthetic and practical TDEM data.  相似文献   

16.
Wavepath traveltime tomography   总被引：1，自引：1，他引：1  
The elastic-wave equation is used to construct sensitivity kernels relating perturbations in elastic parameters to traveltime deviations. Computation of the functions requires a correlation of the forward-propagating seismic wavefield with a backward propagation of the residual wavefield. The computation of the wavefields is accomplished using a finite difference algorithm and is efficiently executed on a CM-2 parallel processor. The source and receiver locations have maximum sensitivity to velocity structure. The sensitivity kernels or wavepaths are well suited for transmission traveltime inversion such as cross-borehole tomography and vertical seismic profiling. Conventional ray tomography and wavepath tomography are applied to a set of P -wave arrival times, from a cross-borehole experiment at Kesterson, California. Because the wavepaths have increased sensitivity near the source and receiver there are differences in resolution of the velocity structure. Both techniques recover the same relative variations in velocity where the coverage is adequate. The wavepath solution is more laterally continuous and the dominant variation is vertical, as is expected for the layered sediments in this region.  相似文献   

17.
Kinematics and dynamics of the southeastern margin of the Tibetan Plateau   总被引：4，自引：0，他引：4  
On the southeastern margin of the Tibetan Plateau lies a large region which seismicity and GPS data show to be actively deforming. This paper describes the active faulting in the region, and how it relates to the velocity field observed with GPS. In places the velocity field is accommodated by rotations about vertical axes, and most or all of the strain at the surface in the region appears to be released seismically. GPS velocities are then compared to velocities calculated using a model for deformation driven by gravitational driving forces. Using rheologies estimated from experimentally derived mineral flow laws, the model provides velocities that are in good agreement with observed GPS velocities. It is not possible to uniquely determine the rheology or flow velocity at depth, and there are two forms of model solution which match the observed horizontal surface velocities. In one of these, vertical planes deform by pure shear, and in the other vertical gradients of horizontal velocity are present within the crust. Two distinct regions of normal-faulting earthquakes are present in the region, and have mechanisms which are most easily explained by gravity-driven deformation.  相似文献   

18.
Wind tunnel modeling and measurements of the flux of wind-blown sand   总被引：1，自引：0，他引：1  
This paper presents a way to empirically fit experimental data for the horizontal flux of various sizes of wind-blown dry sand using data from wind tunnel experiments. We measured vertical wind profiles to derive threshold shear velocity and estimated shear velocity and the flux of sand mass as a function of the height for nine different grain sizes. We propose a fitting model based on the experimental data and a least-squares method and derive an explicit form of sand flux as a function of height and shear velocity for these grain sizes. We also obtained an explicit form of the empirical equation for the measurement of sand transport per unit width and unit time by integrating the empirical equation as a function of height. Finally, we compared the effectiveness of Bagnold's equation, Kawamura's expression and Lettau and Lettau's equation, for predicting sand transport with the results of our empirical equation. The results show that the transport predicted by all of the equations were always lower than the measured results from the empirical equation for all grain sizes and shear velocities. However, the empirical equation matched Bagnold's equation, Kawamura's equation, and Lettau and Lettau's equation if the coefficients in these equations were adjusted instead of using their original coefficients. The empirical equation for sand transport in the present study contradicts previous conclusions generated by Bagnold's equation, which predict that for a given wind drag, the transport of a coarse sand is greater than that of a fine sand.  相似文献   

19.
Summary. The Hatton Bank passive continental margin exhibits thick seaward dipping reflector sequences which consist of basalts extruded during rifting between Greenland and Rockall Plateau. Multichannel seismic reflection profiling across the margin reveals three reflector wedges with a maximum thickness near 7 km, extending from beneath the upper continental slope to the deep ocean basin. We present results of the velocity structure within the dipping reflector sequences at eight locations across the margin, interpreted by synthetic seismogram modelling a set of multichannel expanding spread profiles parallel to the margin. At the top of some reflector sequences, we observe a series of 100 m thick high- and low-velocity zones, which are interpreted as basalt flows alternating with sediments or weathered and rubble layers. At the profile locations, the base of the dipping reflectors correlates with P -wave velocities near 6.5 km s−1. However, elsewhere the reflectors appear to extend significantly deeper than the inferred 6.5 km s−1 velocity contour, indicating that the velocity structure may not be controlled solely by lithological boundaries but also by metamorphic effects. Shear-waves were observed on two lines, permitting the calculation of Poisson's ratio. The decrease in Poisson's ratio from 0.28 to near 0.25 in the upper 5 km of crust may also indicate the effect of metamorphism on seismic properties, or alternatively may be explained by crack closure under load.  相似文献   

20.
Summary. The crustal structure beneath the exposed terranes of southern Alaska has been explored using coincident seismic refraction and reflection profiling. A wide-angle reflector at 8–9 km depth, at the base of an inferred low-velocity zone, underlies the Peninsular and Chugach terranes, appears to truncate their boundary, and may represent a horizontal decollement beneath the terranes. The crust beneath the Chugach terrane is characterized by a series of north-dipping paired layers having low and high velocities that may represent subducted slices of oceanic crust and mantle. This layered series may continue northward under the Peninsular terrane. Earthquake locations in the Wrangell Benioff zone indicate that at least the upper two low-high velocity layer pairs are tectonically inactive and that they appear to have been accreted to the base of the continental crust. The refraction data suggest that the Contact fault between two similar terranes, the Chugach and Prince William terranes, is a deeply penetrating feature that separates lower crust (deeper than 10 km) with paired dipping reflectors, from crust without such reflectors.  相似文献   

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

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