共查询到20条相似文献，搜索用时 390 毫秒
1.
R.-G. FERBER 《Geophysical Prospecting》1991,39(6):737-755
The calculation of dip moveout involves spreading the amplitudes of each input trace along the source-receiver axis followed by stacking the results into a 3D zero-offset data cube. The offset-traveltime (x–t) domain integral implementation of the DMO operator is very efficient in terms of computation time but suffers from operator aliasing. The log-stretch approach, using a logarithmic transformation of the time axis to force the DMO operator to be time invariant, can avoid operator aliasing by direct implementation in the frequency-wavenumber (f–k) domain. An alternative technique for log-stretch DMO corrections using the anti-aliasing filters of the f–k approach in the x-log t domain will be presented. Conventionally, the 2D filter representing the DMO operator is designed and applied in the f–k domain. The new technique uses a 2D convolution filter acting in single input/multiple output trace mode. Each single input trace is passed through several 1D filters to create the overall DMO response of that trace. The resulting traces can be stacked directly in the 3D data cube. The single trace filters are the result of a filter design technique reducing the 2D problem to several ID problems. These filters can be decomposed into a pure time-delay and a low-pass filter, representing the kinematic and dynamic behaviour of the DMO operator. The low-pass filters avoid any incidental operator aliasing. Different types of low-pass filters can be used to achieve different amplitude-versus-offset characteristics of the DMO operator. 相似文献
2.
An accurate analytical expression for shot-gather dip-moveout (DMO) in the timespace log-stretch domain has until now not been published. We present a simpler, alternative derivation of the exact DMO relationships of Black et al. which correctly take account of the repositioning of the midpoint. A new computationally efficient frequency-wavenumber (F-K) DMO operator for shot profiles is then derived, based on these DMO relationships in the time-space log-stretch domain. The newly derived DMO operator is, unlike most other log-stretch DMO operators) accurate for the full range of reflector dips. Along with other schemes which are performed in the log-stretch domain, it offers considerable time savings over conventional DMO processing. We have compared numerically the impulse response of the new operator with those of a number of other shot-gather DMO operators, and found it to be superior and well match to the theoretical elliptical DMO response. 相似文献
3.
Selection of the wavenumbers k using an optimization method for the inverse Fourier transform in 2.5D electrical modelling 总被引：1，自引：0，他引：1
An optimization method is used to select the wavenumbers k for the inverse Fourier transform in 2.5D electrical modelling. The model tests show that with the wavenumbers k selected in this way the inverse Fourier transform performs with satisfactory accuracy. 相似文献
4.
A short convolutional differentiator (CD) for computing second spatial derivatives in the acoustic wave equation is presented. This differentiator is obtained by tapering the inverse Fourier transform of the band-limited Fourier spectrum of the second-derivative operator. This new filter has been applied to seismogram computations for inhomogeneous media and results are compared with the conventional high-order finite-difference (FD) and Fourier schemes. The operator can be progressively shortened at the model edges to reduce boundary artefacts. The CD method is superior to the conventional FD operator and comparable with the Fourier method in accuracy but faster to run. A strategy to reduce computation time by 20%, which exploits the localized nature of the operator, is given. The method is illustrated using simple 2D models. 相似文献
5.
The factors controlling the chemistry of 69 low-order streams in the Blue Ridge and Valley and Ridge physiographic provinces of Virginia and Maryland were studied over a 13-month period. Principal component analysis was used to examine regional patterns in stream chemistry and to examine the degree to which the chemistry of low-order streams is controlled by the bedrock upon which they flow. Streams clustered into regionally isolated groups, strongly related to bedrock type, with SO2?4 and HCO?3 the chemical variables of most importance. Sulphate concentrations appear to be strongly controlled by climate and hydrology, and sorption in the soils within the watershed. Much of the atmospherically derived SO2?4 accumulates in watersheds during the growing season and is later flushed out. Weathering reactions were found to be particularly important in the production of HCO?3, accounting for 91 per cent on an annual basis, and export of divalent cations from these watersheds, accounting for 48–50 per cent on an annual basis. About half of non-anthropogenic Na+ was derived from weathering of silicates, whereas nearly all K+ was identified with leaching by SO2?4. Water chemistry was strongly related to the rock type in the watershed and the weatherability of the component minerals. Rock type is not a randomly distributed function; instead, it is controlled by geologic factors that result in clusters of similar rock types in a given region. When planning large synoptic studies, it is extremely important to consider that a sampling scheme based on random sampling of a non-randomly distributed function May, not provide the most accurate representation of the variables of interest. Instead, a hierarchical sampling scheme May, be indicated. Our results also suggest that, although one sample in time May, be sufficient to characterize the primary geochemical factors controlling stream chemistry throughout the year, it May, not be sufficient to detect subtle, flow-related alterations in chemistry. 相似文献
6.
Ilya Tsvankin 《Geophysical Prospecting》1997,45(3):479-512
The transversely isotropic (TI) model with a tilted axis of symmetry may be typical, for instance, for sediments near the flanks of salt domes. This work is devoted to an analysis of reflection moveout from horizontal and dipping reflectors in the symmetry plane of TI media that contains the symmetry axis. While for vertical and horizontal transverse isotropy zero-offset reflections exist for the full range of dips up to 90°, this is no longer the case for intermediate axis orientations. For typical homogeneous models with a symmetry axis tilted towards the reflector, wavefront distortions make it impossible to generate specular zero-offset reflected rays from steep interfaces. The ‘missing’ dipping planes can be imaged only in vertically inhomogeneous media by using turning waves. These unusual phenomena may have serious implications in salt imaging. In non-elliptical TI media, the tilt of the symmetry axis may have a drastic influence on normal-moveout (NMO) velocity from horizontal reflectors, as well as on the dependence of NMO velocity on the ray parameter p (the ‘dip-moveout (DMO) signature’). The DMO signature retains the same character as for vertical transverse isotropy only for near-vertical and near-horizontal orientation of the symmetry axis. The behaviour of NMO velocity rapidly changes if the symmetry axis is tilted away from the vertical, with a tilt of ±20° being almost sufficient to eliminate the influence of the anisotropy on the DMO signature. For larger tilt angles and typical positive values of the difference between the anisotropic parameters ε and δ, the NMO velocity increases with p more slowly than in homogeneous isotropic media; a dependence usually caused by a vertical velocity gradient. Dip-moveout processing for a wide range of tilt angles requires application of anisotropic DMO algorithms. The strong influence of the tilt angle on P-wave moveout can be used to constrain the tilt using P-wave NMO velocity in the plane that includes the symmetry axis. However, if the azimuth of the axis is unknown, the inversion for the axis orientation cannot be performed without a 3D analysis of reflection traveltimes on lines with different azimuthal directions. 相似文献
7.
Inspired by the linear filter method introduced by D. P. Ghosh in 1970 we have developed a general theory for numerical evaluation of integrals of the Hankel type: Replacing the usual sine interpolating function by sinsh (x) =a· sin (ρx)/sinh (aρx), where the smoothness parameter a is chosen to be “small”, we obtain explicit series expansions for the sinsh-response or filter function H*. If the input function f(λ exp (iω)) is known to be analytic in the region o < λ < ∞, |ω|≤ω0 of the complex plane, we can show that the absolute error on the output function is less than (K(ω0)/r) · exp (?ρω0/Δ), Δ being the logarthmic sampling distance. Due to the explicit expansions of H* the tails of the infinite summation ((m?n)Δ) can be handled analytically. Since the only restriction on the order is ν > ? 1, the Fourier transform is a special case of the theory, ν=± 1/2 giving the sine- and cosine transform, respectively. In theoretical model calculations the present method is considerably more efficient than the Fast Fourier Transform (FFT). 相似文献
8.
A new single-station model (SSM) for monthly median values of the ionospheric parameters foF2 and M(3000)F2 has been developed. Fourier analysis provides a tool for decomposing the time-varying ionospheric parameters. The 12–month smoothed sunspot number R
12
was used as an external solar characteristic because of its availability and predictability. However, for the first time, the solar activity is described not only by R
12
, but also by the linear coefficient K
R
representing the tendency of the change of solar activity. A general non-linear approximation of the influence of the solar-cycle characteristics R
12
and K
R
and ionospheric parameters foF2 and M(3000)F2 was accepted. The new SSM is applied to several European stations and its statistical evaluation shows better results than the other two SSMs used in the paper. The approach described in the paper does not contradict the use of different synthetic ionospheric indices (as the T-index, MF2–index); the basic aim is to show only that using one additional new characteristic of the solar-cycle variations, such as K
R
, improves the monthly median model. 相似文献
9.
We present the first demonstration of hydraulic tomography (HT) to estimate the three-dimensional (3D) hydraulic conductivity (K) distribution of a fractured aquifer at high-resolution field scale (HRFS), including the fracture network and connectivity through it. We invert drawdown data collected from packer-isolated borehole intervals during 42 pumping tests in a wellfield at the former Naval Air Warfare Center, West Trenton, New Jersey, in the Newark Basin. Five additional tests were reserved for a quality check of HT results. We used an equivalent porous medium forward model and geostatistical inversion to estimate 3D K at high resolution (K blocks <1 m3), using no strict assumptions about K variability or fracture statistics. The resulting 3D K estimate ranges from approximately 0.1 (highest-K fractures) to approximately 10−13 m/s (unfractured mudstone). Important estimated features include: (1) a highly fractured zone (HFZ) consisting of a sequence of high-K bedding-plane fractures; (2) a low-K zone that disrupts the HFZ; (3) several secondary fractures of limited extent; and (4) regions of very low-K rock matrix. The 3D K estimate explains complex drawdown behavior observed in the field. Drawdown tracing and particle tracking simulations reveal a 3D fracture network within the estimated K distribution, and connectivity routes through the network. Model fit is best in the shallower part of the wellfield, with high density of observations and tests. The capabilities of HT demonstrated for 3D fractured aquifer characterization at HRFS may support improved in situ remediation for contaminant source zones, and applications in mining, repository assessment, or geotechnical engineering. 相似文献
10.
We introduce the signal dependent time–frequency distribution, which is a time–frequency distribution that allows the user to optimize the tradeoff between joint time–frequency resolution and suppression of transform artefacts. The signal‐dependent time–frequency distribution, as well as the short‐time Fourier transform, Stockwell transform, and the Fourier transform are analysed for their ability to estimate the spectrum of a known wavelet used in a tuning wedge model. Next, the signal‐dependent time–frequency distribution, and fixed‐ and variable‐window transforms are used to estimate spectra from a zero‐offset synthetic seismogram. Attenuation is estimated from the associated spectral ratio curves, and the accuracy of the results is compared. The synthetic consisted of six pairs of strong reflections, based on real well‐log data, with a modeled intrinsic attenuation value of 1000/Q = 20. The signal‐dependent time–frequency distribution was the only time–frequency transform found to produce spectra that estimated consistent attenuation values, with an average of 1000/Q = 26±2; results from the fixed‐ and variable‐window transforms were 24±17 and 39±10, respectively. Finally, all three time–frequency transforms were used in a pre‐stack attenuation estimation method (the pre‐stack Q inversion algorithm) applied to a gather from a North Sea seismic dataset, to estimate attenuation between nine different strong reflections. In this case, the signal‐dependent time‐frequency distribution produced spectra more consistent with the constant‐Q model of attenuation assumed in the pre‐stack attenuation estimation algorithm: the average L1 residuals of the spectral ratio surfaces from the theoretical constant‐Q expectation for the signal‐dependent time‐frequency distribution, short‐time Fourier transform, and Stockwell transform were 0.12, 0.21, and 0.33, respectively. Based on the results shown, the signal‐dependent time‐frequency distribution is a time–frequency distribution that can provide more accurate and precise estimations of the amplitude spectrum of a reflection, due to a higher attainable time–frequency resolution. 相似文献
11.
The τ-p transform is an invertible transformation of seismic shot records expressed as a function of time and offset into the τ (intercept time) and p (ray parameter) domain. The τ-p transform is derived from the solution of the wave equation for a point source in a three-dimensional, vertically non-homogeneous medium and therefore is a true amplitude process for the assumed model. The main advantage of this transformation is to present a point source shot record as a series of plane wave experiments. The asymptotic expansion of this transformation is found to be useful in reflection seismic data processing. The τ-p and frequency-wavenumber (or f-k) processes are closely related. Indeed, the τ-p process embodies the frequency-wavenumber transformation, so the use of this technique suffers the same limitations as the f-k technique. In particular, the wavefield must be sampled with sufficient spatial density to avoid wavenumber aliasing. The computation of this transform and its inverse transform consists of a two-dimensional Fast Fourier Transform followed by an interpolation, then by an inverse-time Fast Fourier Transform. This technique is extended from a vertically inhomogeneous three-dimensional medium to a vertically and laterally inhomogeneous three-dimensional medium. The τ-p transform may create artifacts (truncation and aliasing effects) which can be reduced by a finer spatial density of geophone groups by a balancing of the seismic data and by a tapering of the extremities of the seismic data. The τ-p domain is used as a temporary domain where the attack of coherent noise is well addressed; this technique can be viewed as ‘time-variant f-k filtering’. In addition, the process of deconvolution and multiple suppression in the τ-p domain is at least as well addressed as in the time-offset domain. 相似文献
12.
A boundary element formulation having discontinuous curved quadratic elements is presented for 2D elastodynamics. The first fundamental solution for static case is subtracted from and added to the first fundamental solution for dynamic case. As both kernels have the same order of singularity, the integral involving the regular expression arising from the subtraction can be calculated. matrix is calculated by employing the well-known rigid-body motion technique. The formulation is performed in Fourier transform space. Based on the formulation presented in this study, a general purpose computer program is developed for elastic or visco-elastic 2D elastodynamic problems. The program performs the analysis in Fourier transform space and can also be used for static analysis by assigning a very small value close to zero for the frequency. The results of some elastodynamic and dynamic soil–structure interaction problems obtained using the present study are compared with those in the literature. 相似文献
13.
Abstract The accuracy of six combined methods formed by three commonly-used soil hydraulic functions and two methods to determine soil hydraulic parameters based on a soil hydraulic parameter look-up table and soil pedotransfer functions was examined for simulating soil moisture. A novel data analysis and modelling approach was used that eliminated the effects of evapotranspiration so that specific sources of error among the six combined methods could be identified and quantified. By comparing simulated and observed soil moisture at six sites of the USDA Soil Climate Analysis Network, we identified the optimal soil hydraulic functions and parameters for predicting soil moisture. Through sensitivity tests, we also showed that adjusting only the soil saturated hydraulic conductivity, Ks , is insufficient for representing important effects of macropores on soil hydraulic conductivity. Our analysis illustrates that, in general, soil hydraulic conductivity is less sensitive to Ks than to the soil pore-size distribution parameter. Editor D. Koutsoyiannis; Associate editor D. Hughes Citation Pan, F., McKane, R.B. and Stieglitz, M., 2012. Identification of optimal soil hydraulic functions and parameters for predicting soil moisture. Hydrological Sciences Journal, 57 (4), 723–737. 相似文献
14.
Palanivel Sathishkumar Kannan Balan Thayumanavan Palvannan Seralathan Kamala‐Kannan Byung‐Taek Oh Susana Rodríguez‐Couto 《洁净——土壤、空气、水》2013,41(7):665-672
Laccase from the white‐rot fungus Pleurotus florida, produced under solid‐state fermentation conditions, was used for the decolorization of reactive dye Remazol Brilliant Blue R (RBBR). RBBR was decolorized up to 46% by P. florida laccase alone in 10 min. In the presence of N‐hydroxybenzotriazole (HBT), the rate of decolorization was enhanced 1.56‐fold. Central composite design of response surface methodology with four variables namely, dye, enzyme, redox mediator concentrations, and time at five levels was applied to optimize the RBBR decolorization. The predicted optimum level of variables for maximum RBBR decolorization (87%) was found to be 52.90 mg L?1 (RBBR), 1.87 U mL?1 (laccase), 0.85 mM (HBT), and 7.17 min (time), respectively. The validation results showed that the experimental value of RBBR decolorization (82%) was close to the predicted one. The disappearance of C–N and C–X groups, and a small shift in N–H groups in Fourier‐transform infra red (FTIR) spectroscopy confirms the degradation of RBBR chromophore by laccase enzyme. The phytotoxicity of RBBR was considerably reduced after the treatment with laccase. RBBR decolorization kinetics; Km and Vmax were calculated to be 145.82 mg L?1 and 24.86 mg L?1 min, respectively. 相似文献
15.
As is frequently cited, dispersivity increases with solute travel distance in the subsurface. This behaviour has been attributed to the inherent spatial variation of the pore water velocity in geological porous media. Analytically solving the advection–dispersion equation with distance-dependent dispersivity is extremely difficult because the governing equation coefficients are dependent upon the distance variable. This study presents an analytical technique to solve a two-dimensional (2D) advection–dispersion equation with linear distance-dependent longitudinal and transverse dispersivities for describing solute transport in a uniform flow field. The analytical approach is developed by applying the extended power series method coupled with the Laplace and finite Fourier cosine transforms. The developed solution is then compared to the corresponding numerical solution to assess its accuracy and robustness. The results demonstrate that the breakthrough curves at different spatial locations obtained from the power series solution show good agreement with those obtained from the numerical solution. However, owing to the limited numerical operation for large values of the power series functions, the developed analytical solution can only be numerically evaluated when the values of longitudinal dispersivity/distance ratio eL exceed 0·075. Moreover, breakthrough curves obtained from the distance-dependent solution are compared with those from the constant dispersivity solution to investigate the relationship between the transport parameters. Our numerical experiments demonstrate that a previously derived relationship is invalid for large eL values. The analytical power series solution derived in this study is efficient and can be a useful tool for future studies in the field of 2D and distance-dependent dispersive transport. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
16.
Seislet变换是一种小波类数学变换方法,主要根据不同小波级数上地震同相轴的局部倾角的不同来分析数据.一般意义上,测线方向上的离散小波变换（DWT）是一种特殊的零局部地震倾角的seislet变换.早期的工作基于低阶版本的离散小波变换来构建seislet变换,在本文中,通过使用Cohen-Daubechies-Feauveau (CDF) 9/7双正交小波变换（常用于JPEG2000压缩标准）作为框架,扩展高阶seislet变换方法.通过分析理论模型和实际数据的处理结果,并对比傅里叶变换、离散小波变换和低阶seislet变换,高阶seislet变换可以为地震数据提供更好的压缩比.因此更加适用于地震数据去噪处理. 相似文献
17.
《Advances in water resources》2007,30(6-7):1571-1582
It has been widely observed in field experiments that the apparent rate of bacterial attachment, particularly as parameterized by the collision efficiency in filtration-based models, decreases with transport distance (i.e., exhibits scale-dependency). This effect has previously been attributed to microbial heterogeneity; that is, variability in cell–surface properties within a single monoclonal population. We demonstrate that this effect could also be interpreted as a field-scale manifestation of local-scale correlation between physical heterogeneity (hydraulic conductivity variability) and reaction heterogeneity (attachment rate coefficient variability). A field-scale model of bacterial transport developed for the South Oyster field research site located near Oyster, Virginia, and observations from field experiments performed at that site, are used as the basis for this study. Three-dimensional Monte Carlo simulations of bacterial transport were performed under four alternative scenarios: (1) homogeneous hydraulic conductivity (K) and attachment rate coefficient (Kf); (2) heterogeneous K, homogeneous Kf; (3) heterogeneous K and Kf with local correlation based on empirical and theoretical relationships; and (4) heterogeneous K and Kf without local correlation. The results of the 3D simulations were analyzed using 1D model approximations following conventional methods of field data analysis. An apparent decrease with transport distance of effective collision efficiency was observed only in the case where the local properties were both heterogeneous and correlated. This effect was observed despite the fact that the local collision efficiency was specified as a constant in the 3D model, and can therefore be interpreted as a scale effect associated with the local correlated heterogeneity as manifested at the field scale. 相似文献
18.
Dip‐moveout (DMO) correction is often applied to common‐offset sections of seismic data using a homogeneous isotropic medium assumption, which results in a fast execution. Velocity‐residual DMO is developed to correct for the medium‐treatment limitation of the fast DMO. For reasonable‐sized velocity perturbations, the residual DMO operator is small, and thus is an efficient means of applying a conventional Kirchhoff approach. However, the shape of the residual DMO operator is complicated and may form caustics. We use the Fourier domain for the operator development part of the residual DMO, while performing the convolution with common‐offset data in the space–time domain. Since the application is based on an integral (Kirchhoff) method, this residual DMO preserves all the flexibility features of an integral DMO. An application to synthetic and real data demonstrates effectiveness of the velocity‐residual DMO in data processing and velocity analysis. 相似文献
19.
《Journal of Atmospheric and Solar》2000,62(2):133-140
In 1995–1996, observations were carried out at Norilsk (geomagnetic latitude and longitude 64.2°N and 160.4°E) to determine dynamic parameters of irregularities in the high-latitude ionosphere. The short-baseline spaced-receiver method that has been implemented at the ionospheric facility of the Norilsk Integrated Magnetic–Ionospheric Station, provides a means of simultaneously measuring parameters of small-scale irregularities (spatial scale of 3–5 km) by the Similar-Fading Method (SFM), as well as of medium-scale irregularities (time scale of 10–30 min, spatial scale of hundreds of kilometres) by the Statistical Angle-of-arrival and Doppler Method (SADM). About 20 h of the observational data for the F2-layer under quiet geomagnetic conditions (Kp < 3), 20 h under disturbed conditions (Kp ≥ 3) and about 15 h for the sporadic E-layer (Kp ≈ 3) were processed. It has been found that the propagation directions and velocities of different-scale irregularities do not coincide. Small-scale irregularities of the F2-layer travel predominantly eastward or westward. The velocity of the F2-layer irregularities is about 100 m/s, and under disturbed conditions it is up to 200–250 m/s. Small-scale irregularities of the sporadic E-layer travel mostly in the northward direction. It is confirmed that the Es-layer is characterised by high velocities of the irregularities (as high as 1000 m/s). Medium-scale irregularities with periods in the range of 10–30 min travel mostly in a southward direction with velocities of 20–40 m/s. 相似文献
20.
Mark Meadows Don Adams Rich Wright † Ali Tura Steve Cole David Lumley 《Geophysical Prospecting》2005,53(2):205-213
Rock physics analysis plays a vital role in time‐lapse seismic interpretation because it provides the link between changes in rock and fluid properties and the resulting seismic data response. In this case study of the Schiehallion Field, we discuss a number of issues that commonly arise in rock physics analyses for time‐lapse studies. We show that:
- 1 Logarithmic fits of dry bulk (Kdry) and shear (Gdry) moduli vs. effective pressure (Peff) are superior to polynomial fits.
- 2 2D surface fits of Kdry and Gdry over porosity (φ) and effective pressure using all the core data simultaneously are more useful and accurate than separate 1D fits over φ and Peff for each individual core.
- 3 One average set (facies) of Kdry(φ, Peff) and Gdry(φ, Peff) can be chosen to represent adequately the entire Schiehallion reservoir.
- 4 Saturated velocities and densities modelled by fluid substitution of Kdry(φ, Peff), Gdry(φ, Peff) and the dry bulk density ρdry(φ) compare favourably with well‐log velocities and densities.
- 5 P‐ and S‐wave impedance values resulting from fluid substitution of Kdry(φ, Peff), Gdry(φ, Peff) and ρdry(φ) show that the largest impedance changes occur for high porosities and low effective pressures.
- 6 Uncertainties in Kdry(φ, Peff) and Gdry(φ, Peff) derived for individual cores can be used to generate error surfaces for these moduli that represent bounds for quantifying uncertainties in seismic modelling or pressure–saturation inversion.