首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The slant-stack migration formula based on the Radon transform is studied with respect to the depth step Δz of wavefield extrapolation. It can be viewed as a generalized trace-interpolation procedure including wave extrapolation with an arbitrary step Δz. For Δz= 0 the formula yields the familiar plane-wave decomposition, while for Δz > 0 it provides a robust tool for migration transformation of spatially undersampled wavefields. Using the stationary phase method, it is shown that the slant-stack migration formula degenerates into the Rayleigh-Sommerfeld integral in the far-field approximation. Consequently, even a narrow slant-stack gather applied before the diffraction stack can significantly improve the representation of noisy data in the wavefield extrapolation process. The theory is applied to synthetic and field data to perform trace interpolation and dip reject filtration. The data examples presented prove that the Radon interpolator works well in the dip range, including waves with mutual stepouts smaller than half the dominant period.  相似文献   

2.
For transient, high frequency, and pulse like excitation of structures in the near field of strong earthquakes, the classical design approach based on relative response spectrum and mode superposition may not be conservative. For such excitations, it is more natural to use wave propagation methods. In this paper (Part I), we review several two-dimensional wave propagation models of buildings and show results for theoretical dispersion curves computed for these models. We also estimate the parameters of these models that would correspond to a seven-story reinforced concrete building in Van Nuys, California. Ambient vibration tests data for this building imply vertical shear wave velocity βz=112 m/s and anisotropy factor βx/βz=0.55 for NS vibrations, and βz=88 m/s and βx/βz=1 for EW vibrations. The velocity of shear waves propagating through the slabs is estimated to be about 2000 m/s. In the companion paper (Part II), we estimate phase velocities of vertically and horizontally propagating waves between seven pairs of recording points in the building using recorded response to four earthquakes.  相似文献   

3.
A migration algorithm appropriate for moderately varying lateral velocity changes is developed as an extension of phase-shift migration by using a variable-length spatial transform. This process significantly reduces the number of lateral wave-numbers necessary to downward continue the data, and it replaces the spatial FFT with a simple recursion relationship. For a given frequency and position x, ten lateral wavenumbers are typically sufficient, and the migration algorithm produces accurate images when the velocity structure V (X, z) changes over a few depth intervals of thickness Δz, with lateral velocity gradients up to 1.4 to 1.0.  相似文献   

4.
The objective of this research was to develop and parameterise a physically justified yet low‐parameter model to quantify observed changes in surface runoff ratios with hillslope length. The approach starts with the assumption that a unit of rainfall‐excess runoff generated at a point is a fraction β of precipitation P (m) which travels some variable distance down a slope before reinfiltrating, depending on the local rainfall, climate, soils, etc. If this random distance travelled Y is represented by a distribution, then a survival function will describe the probability of this unit of runoff travelling further than some distance x (m). The total amount of per unit width runoff Q (m2) flowing across the lower boundary of a slope of length λ (m) may be considered the sum of all the proportions of the units of rainfall excess runoff integrated from the lower boundary x = 0 to the upper boundary x = λ of the slope. Using these assumptions we derive a model Q(λ) = βPμλ/(μ + λ), (μ > 0, 0 ≤ β ≤ 1, λ ≥ 0) that describes the change in surface runoff with slope length, where μ (m) is the mean of the random variable Y. Dividing both sides of this equation by Pλ yields a simple two‐parameter equation for the dimensionless hillslope runoff ratio Qh(λ) = βμ/(μ + λ). The model was parameterised with new rainfall and runoff data collected from three replicates of bounded 2 m wide plots of four different lengths (0.5, 1.0, 2.0 and 4.0 m) for 2 years from a forested SE Australian site, and with 32 slope length–runoff data sets from 12 other published studies undertaken between 1934 and 2010. Using the parameterised model resulted in a Nash and Sutcliffe statistic between observed and predicted runoff ratio (for all data sets combined) of 0.93, compared with –2.1 when the runoff ratio was fixed at the value measured from the shortest plot. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
Results of studies carried out with the help of a three-dimensional seismic model on waves diffracted from edges of varying radius of curvature R and depth h with respect to wave length λ are described. The amplitude decay, travel time, and apparent velocity of the wave diffracted from a sub-surface edge of semi-infinite length are found to depend on the parameters R, h, and distance from the edge on the surface provided the ratio of the parameters to λ are less than some limiting values. The nature of the amplitude decay is independent of R when the depth exceeds 2λ, and independent of h when R exceeds 1.5λ. When these are below the limiting values (h= 2λ and R= 1.5λ), the nature of the decay depends appreciably on R and h. The apparent decay in amplitude on the surface due to geometrical spreading by the diffracting edge is less than that of a cylindrical secondary wave source and decreases with increase in depth of the edge. The nature of the travel time curves of the diffracted waves near the edge depend on R/λ when the depth is within about one λ. Apparent velocity of the wave depends largely on R/λ in the zone of diffraction up to a distance of about one λ from the edge on the surface. Beyond this distance the velocity is almost the same irrespective of R/λ and depend only on h/λ. The width of the zone of diffraction caused by an edge of finite length comparable to λ is more and more narrow as the ratio of the distance of the edge on the surface to its depth increases.  相似文献   

6.
Acoustic plane wave scattering at a vertical fault structure represents the simplest two-dimensional model of geophysical exploration that can be investigated by analytical techniques. The exact and complete solution, in the time domain, for the scattering of the pressure field of an acoustic plane wave normally incident on a vertical fault structure is determined adapting previous results given for the frequency domain. The wave form of the pressure field of the incident plane wave is expressed by a causal time function that decays exponentially with time at every point above the fault (z<0). The zero-order term of the scattered pressure field has been computed above the fault. This zero-order term consists of an inverse Fourier transform which reduces to a closed expression forx=0, and contains an integral of a Hankel function forx#0. The high frequency part of the inverse Fourier transform forx#0 is computed employing asymptotic expressions for the Hankel function. The integral of the asymptotic expression of the Hankel function reduces to: (i) a Fresnel integral which contains a plane wave term for |x||z|; and (ii) a stationary point plane wave term plus an upper limit term for |x|=O(|z|). For the latter case the plane wave term cancels, leaving a cylindrical wave emanated from the edge of the fault. The wave front is well defined in shape, in phase and in amplitude. The amplitude of the scattered field is discontinuous atx=0, presents a jump and is well defined for |x| small and is rather smooth for |x| large.  相似文献   

7.
One of the most important steps in the conventional processing of reflection seismic data is common midpoint (CMP) stacking. However, this step has considerable deficiencies. For instance the reflection or diffraction time curves used for normal moveout corrections must be hyperbolae. Furthermore, undesirable frequency changes by stretching are produced on account of the dependence of the normal moveout corrections on reflection times. Still other drawbacks of conventional CMP stacking could be listed.One possibility to avoid these disadvantages is to replace conventional CMP stacking by a process of migration to be discussed in this paper. For this purpose the Sherwood-Loewenthal model of the exploding reflector has to be extended to an exploding point model with symmetry to the lineP EX M whereP EX is the exploding point, alias common reflection point, andM the common midpoint of receiver and source pairs.Kirchhoff summation is that kind of migration which is practically identical with conventional CMP stacking with the exception that Kirchhoff summation provides more than one resulting trace.In this paper reverse time migration (RTM) was adopted as a tool to replace conventional CMP stacking. This method has the merit that it uses the full wave equation and that a direct depth migration is obtained, the velocityv can be any function of the local coordinatesx, y, z. Since the quality of the reverse time migration is highly dependent on the correct choice of interval velocities such interval velocities can be determined stepwise from layer to layer, and there is no need to compute interval velocities from normal moveout velocities by sophisticated mathematics or time consuming modelling. It will be shown that curve velocity interfaces do not impair the correct determination of interval velocities and that more precise velocity values are obtained by avoiding or restricting muting due to non-hyperbolic normal moveout curves.Finally it is discussed how in the case of complicated structures the reverse time migration of CMP gathers can be modified in such a manner that the combination of all reverse time migrated CMP gathers yields a correct depth migrated section. This presupposes, however, a preliminary data processing and interpretation.  相似文献   

8.
  相似文献   

9.
It has generally been assumed that diffusive sediment transport on soil‐mantled hillslopes is linearly dependent on hillslope gradient. Fieldwork was done near Santa Barbara, California, to develop a sediment transport equation for bioturbation by the pocket gopher (Thomomys bottae) and to determine whether it supports linear diffusion. The route taken by the sediment is divided into two parts, a subsurface path followed by a surface path. The first is the transport of soil through the burrow to the burrow opening. The second is the discharge of sediment from the burrow opening onto the hillslope surface. The total volumetric sediment flux, as a function of hillslope gradient, is found to be: qs (cm3 cm−1 a−1) = 176(dz/dx)3 − 189(dz/dx)2 + 68(dz/dx) + 34(dz/dx)0·4. This result does not support the use of linear diffusion for hillslopes where gopher bioturbation is the dominant mode of sediment transport. A one‐dimensional hillslope evolution program was used to evolve hillslope profiles according to non‐linear and linear diffusion and to compare them to a typical hillslope. The non‐linear case more closely resembles the actual profile with a convex cap at the divide leading into a straight midslope section. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

10.
A one-dimensional non-linear diffusion wave equation is derived from the Saint Venant equations with neglect of the inertia terms. This non-linear equation has no general analytical solution. Numerical schemes are therefore employed to discretize the space and time axes and convert the differential equation to difference form. In this study, the mixing cell method is used to convert the diffusion wave equation to difference form, in which the difference term can be eliminated by selecting an optimal space step size Δx when time step size Δt is given. When the time step size Δt→0, the space step size Δx=Q/(2S0BC]k) where Q is discharge, S0 is bed slope, B is channel width and Ck is kinematic wave celerity, which is the same as the characteristic length proposed by Kalinin and Milyukov. The results of application to two cases show that the mixing cell and linear channel flow routing methods produce hydrographs that are in agreement with the observed flood hydrographs. © 1997 John Wiley & Sons, Ltd.  相似文献   

11.
The gravitational attraction produced by a geological body of irregular shape can be easily determined by dividing it into cubes of small size. The exact expression of the effect of a cube is very complicated, but it can be calculated by using an electronic computer. 4851 values of the gravitational attraction were determined for different positions of a cube with the side of 2l and the center in M(x0, y0, z0), for x0∈[0;20], y0∈[0;20] and z0∈[0;10]. Using these values, templates in double logarithmic representation were drawn, with x0 and z0 as parameters and y0 as abscissa, and with x0 and y0 as parameters and z0 in abscissa; this double set of templates permits a good interpolation for all cubes in the considered domain of M. The use of templates was tested to approximate the effect produced by a theoretical model of spherical shape and in a real case of a three-dimensional salt body of known shape based on a large number of boreholes. In both cases very good results were obtained.  相似文献   

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

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

14.
Summary The best fit curves for upper air mean dry-bulb and dew-point temperatures over Gauhati airport (26°05N, 91°43E, 49 metres a.m.s.l.), for the month of April, have been calculated with the equation,x=A+By+Cy 2,y being the log value in mb of the isobaric level under consideration andx, the mean dry-bulb or dew-point temperature as the case may be, at the isobaric level under consideration. The values of constantsA, B andC for morning dry-bulb and dew-point curves come to be –29.54559, –93.65766 and +37.35048 and –118.84791, –31.15503 and +25.63585 respectively and values of these constants for evening curves come to be –35.86214, –94.15694 and +38.61870 and –127.55970, –29.97192 and +26.36538 respectively. These best fit curves help in finding out mean desired temperatures at any isobaric level in forecasting of thunderstorms and hailstorms, at a station, by dry-bulb and dew-point temperature anomaly technique propsed earlier by the authors.  相似文献   

15.
Modeling and migration couple is one of the most important steps in seismic data processing and interpretation. Absorbing boundary conditions used in the modeling were studied with the wave-equation by different authors. In this study, reflection coefficient analyses of recent solutions are compared to each other for the different incident angles of seismic waves to the modeling boundaries. According to the reflection coefficients correlation, the A3 condition is the most appropriate solution which greatly reduces artificial reflections from the boundaries. However, multi-transmitting Formula is better for relatively high angles between 32–58° with the usage of a special parameter. On the contrary, this formula is not an appropriate condition for angles lower than 32°, although it allows setting the boundary at any preferred angle. Considering that most of the boundaries are set in low angles, A3 solution is still most preferential condition. In this study, it is also aimed to find out the optimum grid intervals for minimizing the ill-posedness arose from the combination of the 45° finite difference migration equation and the B2 absorbing boundary condition for migration. Appropriate values are determined as ωΔx = 0.2 and ωΔz = 0.4 or neighbouring coarser values. It is also concluded that finer mesh spacing can increase the ill-posedness, in contrast to existence of some fine grid size values providing well-posedness. In addition, ill-posedness is obviously standard after ωΔx = 0.6 for all values of ωΔz.  相似文献   

16.
A mathematical expression for potential of a direct current point source in an inhomoge-neous anisotropic earth is derived. The coefficient of anisotropy is given by f= (σr/σz), where σr and σz are the conductivities parallel and perpendicular to the bedding plane. It is assumed that σz varies with depth, whereas σr varies transversely. This potential may be useful in interpretation of geoelectrical data in specified geological situations. Master curves for Wenner and Schlumberger configurations are presented  相似文献   

17.
The arrival-time curve of a reflection from a horizontal interface, beneath a homogeneous isotropic layer, is a hyperbola in the x - t-domain. If the subsurface is one-dimensionally inhomogeneous (horizontally layered), or if some or all of the layers are transversely isotropic with vertical axis of symmetry, the statement is no longer strictly true, though the arrival-time curves are still hyperbola-like. In the case of transverse isotropy, however, classical interpretation of these curves fails. Interval velocities calculated from t2 - x2-curves do not always approximate vertical velocities and therefore cannot be used to calculate depths of reflectors. To study the relationship between velocities calculated from t2 - x2-curves and the true velocities of a transversely isotropic layer, we approximate t2 - x2-curves over a vertically inhomogeneous transversely isotropic medium by a three-term Taylor series and calculate expressions for these terms as a function of the elastic parameters. It is shown that both inhomogeneity and transverse isotropy affect slope and curvature of t2 - x2-curves. For P-waves the effect of transverse isotropy is that the t2 - x2-curves are convex upwards; for SV-waves the curves are convex downwards. For SH-waves transverse isotropy has no effect on curvature.  相似文献   

18.
Summary This paper considers an incompressible fluid flowing through a straight, circular tube whose walls are uniformly porous. The flow is steady and one dimensional. The loss of fluid through the wall is proportional to the mean static pressure in the tube. Several formulations of the wall shear stress are considered; these formulations were motivated by the results from Hamel's radial flow problem, boundary layer flows/and boundary layer suction profiles. For each of these formulations exact solutions for the mean axial velocity and the mean static pressure of the fluid are obtained. Sample results are plotted on graphs. For the constant wall shear stress problem, the theoretical solutions compare favorably with some experimental results.Notations A, B, D, E constant parameters - a, b constant parameters - Ai(z), Bi(z) Airy functions - Ai, Bi derivatives of Airy functions - k constant of proportionality betweenV andp - L length of pores - p,p mean static pressure - p 0 static pressure outside the tube - p 0 value ofp atx=0 - Q constant exponent - R inside radius of the tube - T wall shear stress - T 0 shear parameter - t wall thickness - U free stream velocity - ,u mean axial velocity - u 0 value ofu atx=0 - V,V mean seepage velocity through the wall - v 0 mean seepage velocity - x,x axial distance along the tube - z transformed axial distance - z 0 value ofz atx=0 - mean outflow angle through the wall - cos - density of the fluid - wall shear stress - dynamic viscosity of the fluid - over-bar dimensional terms - no bar nondimensional terms The National Center for Atmospheric Research is sponsored by the National Science Foundation  相似文献   

19.
In the western coal-mining area of Ruhrkohle AG, reflection seismic prospecting for the Carboniferous coal measures is severely impaired by structures with halokinetic features. These structures make the interface between Mesozoic and Paleozoic layers, i.e., the top of Zechstein in general, very rugged. Unfortunately the velocity contrast at this interface is very high in that area, the ratio of velocities being 1.5 to 2.0. Therefore, migration and stacking become a problem. Three types of migration are presented:
  • 1 (f, x)-time-migration with vertical time-to-depth conversion as a second step.
  • 2 Kirchhoff migration down to a level determined approximately by the highest points of the top of Paleozoics, i.e., 0.35 s, and Kirchhoff-downward continuation for all times exceeding 0.35 s. Intermediate static corrections for these latter times with subsequent (f, k)-time-migration and final vertical time-to-depth conversion.
  • 3 Direct depth migration in the (f, x)-domain using three interval velocities.
In all cases an intermediate picking of the velocity interfaces is necessary. In case 2 this occurs at an earlier stage of the process than in case 1, and in case 3 at a still earlier stage. The results of the second and third migration procedures are superior to those of the first. Possibilities for misinterpretation of faults are reduced considerably when the second or third migration procedure is applied.  相似文献   

20.
We show that it is possible to estimate the background velocity for prestack depth migration in 2D laterally varying media using a non-linear optimization technique called very fast simulated annealing (VFSA). We use cubic splines in the velocity model parametrization and make use of either successive pairs of shot gathers or several constant-offset sections as input data for the inversion. A Kirchhoff summation scheme based on first-arrival traveltimes is used to migrate/model the input data during the velocity analysis. We evaluate and compare two different measures of error. The first is defined in the recorded data or (x,t) domain and is based on a reflection-tomography criterion. The second is defined in the migrated data or (x,z) domain and is based on a migration-misfit criterion. Depth relaxation is used to improve the convergence and quality of the velocity analysis while simultaneously reducing the computational cost. Further, we show that by coarse sampling in the offset domain the method is still robust. Our non-linear optimization approach to migration velocity analysis is evaluated for both synthetic and real seismic data. For the velocity-analysis method based on the reflection-tomography criterion, traveltimes do not have to be picked. Similarly, the migration-misfit criterion does not require that depth images be manually compared. Interpreter intervention is required only to restrict the search space used in the velocity-analysis problem. Extension of the proposed schemes to 3D models is straightforward but practical only for the fastest available computers.  相似文献   

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

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