共查询到20条相似文献，搜索用时 78 毫秒
1.
P. HUBRAL 《Geophysical Prospecting》1974,22(4):722-735
The implementation of a stacking filter involves the filtering of each trace with an individual filter and the subsequent summing of all outputs. The actual position of a trace in space as well as certain simultaneous shifts of traces and filter components in time do not influence the process. The resulting output is consequently invariant to various arbitrary coordinate transformations. For a certain useful class of ensembles of non-linear moveout arrival times for signals a particular transformation can be found which transforms a given ensemble into one consisting only of straight lines. It is thus possible to reduce, for instance, the analysis of a stacking filter designed for hyperbola-like moveout curves to the analysis of a velocity filter with linear moveout curves. As the (f—k) transform is a very useful concept to describe a velocity filter, it can consequently be applied to characterize a stacking filter in regard to its performance on input signals with non-linear moveout. 相似文献
2.
R. J. WANG 《Geophysical Prospecting》1977,25(2):342-381
The Wiener prediction filter has been an effective tool for accomplishing dereverberation when the input data are stationary. For non-stationary data, however, the performance of the Wiener filter is often unsatisfactory. This is not surprising since it is derived under the stationarity assumption. Dereverberation of nonstationary seismic data is here accomplished with a difference equation model having time-varying coefficients. These time-varying coefficients are in turn expanded in terms of orthogonal functions. The kernels of these orthogonal functions are then determined according to the adaptive algorithm of Nagumo and Noda. It is demonstrated that the present adaptive predictive deconvolution method, which combines the time-varying difference equation model with the adaptive method of Nagumo and Noda, is a powerful tool for removing both the long- and short-period reverberations. Several examples using both synthetic and field data illustrate the application of adaptive predictive deconvolution. The results of applying the Wiener prediction filter and the adaptive predictive deconvolution on nonstationary data indicate that the adaptive method is much more effective in removing multiples. Furthermore, the criteria for selecting various input parameters are discussed. It has been found that the output trace from the adaptive predictive deconvolution is rather sensitive to some input parameters, and that the prediction distance is by far the most influential parameter. 相似文献
3.
M. R. Thapar 《Pure and Applied Geophysics》1971,92(1):5-18
Summary This paper deals with the problem of finding a Wiener filter when the length of the filter output in not larger than the length of the filter input. Measure for the efficiency of a filter is defined in terms of the relation between the desired filter output and the actual filter output. This measure, called the filter efficiency, is used to find the optimum length of the filter memory function. The relation between the signal-to-noise-ratio (SNR) of the filter input and the SNR of the filter output is discussed. It is shown that there is always some improvement in the SNR through the use of a Wiener filter. 相似文献
4.
拾振器输入信号恢复方法的研究 总被引：1，自引：0，他引：1
为了补偿由于拾振器及其测量系统有限带宽而使测量结果产生的测量误差，本文提出了一种用FIR横向滤波器来恢复拾振器原始输入信号的方法。该方法是将拾振器的输出信号作为FIR横向滤波器的输入，用自适应RLS算法对FIR横向滤波器系数进行辨识而最终获得可以恢复拾振器原始输入信号的反卷积滤波器。通过用几组不同特点的信号对其进行实验检验，结果证明该方法的有效性和较好的泛化能力。 相似文献
5.
用于地震观测的瞬态滤波器仿真研究 总被引：1，自引：0，他引：1
现代地震数据采集器中广泛使用的最小相位滤波器和线性相位滤波器各有不足之处,最小相位滤波引起的“振铃”尾波,线性相位滤波引起的“前缀”波,无法用后续处理的方法消除。本文提出一种专用于地震波形采集的瞬态滤波器设想,该瞬态滤波器既能满足最小相位特性,又能保持较好的线性相位特性。本文用Matlab设计出了可以满足上述条件的瞬态滤波器,并通过仿真比较了正弦接入波通过3种滤波器的输出波形。从仿真结果可以看出,与最小相位和线性相位滤波相比较,该瞬态滤波器的波形失真最小。 相似文献
6.
RALPH A. WIGGINS 《Geophysical Prospecting》1967,15(2):288-296
The design of least-squares optimum filters is based upon minimizing a suitably defined error criterion. The expected value of this error is easily computable after the coefficients of the filter have been determined. When a particular filtering problem is specified, there are several parameters which are specifically not included in the optimization procedure. However, the magnitude of the expected error may be quite sensitive to these parameters. The examination of the relative values of the expected error for variations of these unspecified parameters may lead to a better definition of the filter problem. The parameters which are left unspecified by the general least-square filter definition include: 1. The addition of white noise to the signal autocorrelation to stabilize the filter behavior. 2. The specification of the shape of the desired output of the filter. 3. The specification of the lag between the desired output and the input. Examples are given showing the relationship between these parameters and the value of the expected error. 相似文献
7.
The technique of linear digital filtering developed for the computation of standard curves for conventional resistivity and electromagnetic depth soundings is applied to the determination of filter coefficients for the computation of dipole curves from the resistivity transform function by convolution. In designing the filter function from which the coefficients are derived, a sampling interval shorter than the one used in the earlier work on resistivity sounding is found to be necessary. The performance of the filter sets is tested and found to be highly accurate. The method is also simple and very fast in application. 相似文献
8.
One of the main objectives of seismic digital processing is the improvement of the signal-to-noise ratio in the recorded data. Wiener filters have been successfully applied in this capacity, but alternate filtering devices also merit our attention. Two such systems are the matched filter and the output energy filter. The former is better known to geophysicists as the crosscorrelation filter, and has seen widespread use for the processing of vibratory source data, while the latter is. much less familiar in seismic work. The matched filter is designed such that ideally the presence of a given signal is indicated by a single large deflection in the output. The output energy filter ideally reveals the presence of such a signal by producing a longer burst of energy in the time interval where the signal occurs. The received seismic trace is assumed to be an additive mixture of signal and noise. The shape of the signal must be known in order to design the matched filter, but only the autocorrelation function of this signal need be known to obtain the output energy filter. The derivation of these filters differs according to whether the noise is white or colored. In the former case the noise autocorrelation function consists of only a single spike at lag zero, while in the latter the shape of this noise autocorrelation function is arbitrary. We propose a novel version of the matched filter. Its memory function is given by the minimum-delay wavelet whose autocorrelation function is computed from selected gates of an actual seismic trace. For this reason explicit knowledge of the signal shape is not required for its design; nevertheless, its performance level is not much below that achievable with ordinary matched filters. We call this new filter the “mini-matched” filter. With digital computation in mind, the design criteria are formulated and optimized with time as a discrete variable. We illustrate the techniques with simple numerical examples, and discuss many of the interesting properties that these filters exhibit. 相似文献
9.
Jalal Ferahtia Kamel Baddari Nouredine Djarfour Abdel Kader Kassouri 《Pure and Applied Geophysics》2010,167(11):1389-1404
Seismic noise is a fundamental part of seismic data which cannot be avoided when conducting any seismic survey. It consists
of coherent and random noise. Noise removal or filtering is one of the major concerns in the field of seismic processing.
In this paper, we introduce an image filtering technique based on a detection-estimation algorithm for Gaussian and random
noise removal in seismic data, namely the trilateral filter, based on a statistic called rank-ordered absolute differences.
The non-linear and adaptive behaviour of this filter makes it very robust in the presence of random and coherent noise, in
addition to its computational simplicity and its ability to automatically identify noise in data. We have modified the strategy
of trilateral filtering by adapting the rank-ordered absolute differences formula in order to extract the signal component.
We have successfully used this filter for the removal of surface waves and random spiky noise from synthetic and field data.
Results are very encouraging and show the superiority of this filter compared with other filters, particularly when used recursively. 相似文献
10.
Modeling Full Seismogram Envelopes Using Radiative Transfer Theory with Born Scattering Coefficients
The equation of radiative transfer is used to model the transport of seismic energy in 2-D and 3-D acoustic random media.
Monte-Carlo solutions of this equation using non-isotropic Born scattering coefficients are compared to three analytical solutions:
Markov approximation, radiative transfer theory with isotropic scattering coefficients, and diffusion approximation. Additionally,
we compare to finite differences solutions of the full wave equation in 2-D. We find a good correspondence of radiative transfer
theory to Markov approximation for the case of multiple forward scattering. The comparison to radiative transfer theory with
isotropic scattering coefficients, a model frequently used in data analysis, demonstrates that in the case of forward scattering
the isotropic scattering model is not better than a diffusion approach. To compare radiative transfer theory with non-isotropic
scattering coefficients to finite differences solutions of the full wave equation, the finite source duration and the bandpass
filter process as well as the normalization of absolute amplitudes are explicitely taken into account. We find a good coincidence
of both theories for scattering parameters, which are realistic for usual Earth crust. The theory correctly describes the
unscattered direct wavefront, the envelope broadening caused by multiple forward scattering, as well as the late coda caused
by multiple wide angle scattering. For strong scattering, which can be expected for very heterogeneous media such as strato
volcanoes, the solutions of radiative transfer differ from the more complete solutions of the full wave equation. 相似文献
11.
A partially non-ergodic ground-motion prediction equation is estimated for Europe and the Middle East. Therefore, a hierarchical model is presented that accounts for regional differences. For this purpose, the scaling of ground-motion intensity measures is assumed to be similar, but not identical in different regions. This is achieved by assuming a hierarchical model, where some coefficients are treated as random variables which are sampled from an underlying global distribution. The coefficients are estimated by Bayesian inference. This allows one to estimate the epistemic uncertainty in the coefficients, and consequently in model predictions, in a rigorous way. The model is estimated based on peak ground acceleration data from nine different European/Middle Eastern regions. There are large differences in the amount of earthquakes and records in the different regions. However, due to the hierarchical nature of the model, regions with only few data points borrow strength from other regions with more data. This makes it possible to estimate a separate set of coefficients for all regions. Different regionalized models are compared, for which different coefficients are assumed to be regionally dependent. Results show that regionalizing the coefficients for magnitude and distance scaling leads to better performance of the models. The models for all regions are physically sound, even if only very few earthquakes comprise one region. 相似文献
12.
This study presents an effective method for identifying predictive models and the underlying modal parameters of linear structural systems using only measured output and excitation time histories obtained from dynamic testing. The system under examination is modelled as a first‐order multi‐input multi‐output time‐invariant system, and the structural model is realized using the Eigensystem Realization Algorithm together with the Observer/Kalman filter IDentification algorithm. The identified state‐space model is further refined using a non‐linear optimization technique based on sequential quadratic programming. The numerical examples show that the developed methodology performs very well even in the presence of inadequate instrumentation and measurement noise, and that the methodology is highly capable of creating realistic predictive models of structural systems, as well as estimating their underlying modal parameters. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
13.
We start from the Hankel transform of Stefanescu's integral written in the convolutionintegral form suggested by Ghosh (1971). In this way it is possible to obtain the kernel function by the linear electric filter theory. Ghosh worked out the sets of filter coefficients in frequency domain and showed the very low content of high frequencies of apparent resistivity curves. Vertical soundings in the field measure a series of apparent resistivity values at a constant increment Δx of the logarithm of electrode spacing. Without loss of information we obtain the filter coefficient series by digital convolution of the Bessel function of exponential argument with sine function of the appropriate argument. With a series of forty-one values we obtain the kernel functions from the resistivity curves to an accuracy of better than 0.5%. With the digital method it is possible to calculate easily the filter coefficients for any electrode arrangement and any cut-off frequency. 相似文献
14.
Existing techniques of deconvolution of gravity anomalies are principally based on upward and downward continuation of measured fields. It can be shown that a unique set of linear filters, depending only on geometrical parameters, relates density distribution at a given depth to gravity measured on the surface. A method to compute the filter coefficients is developed. Very accurate reconstitution of theoretical models of intricate shape, prove the validity of the linear relationship. One of these sets of linear filters is applied to a field case of underground quarries. 相似文献
15.
L. MANSINHA 《Geophysical Prospecting》1984,32(6):1155-1166
Forward filters to transform the apparent resistivity function over a layered half-space into the resistivity transform have been derived for a number of sample intervals. The filters have no apparent Gibbs' oscillations and hence require no phase shift. In addition, the end points of the filter were modified to compensate for truncation. The filters were tested on simulated ascending and descending two-layer cases. As expected, “dense” filters with sample spacing of In (10)/6 or smaller performed very well. However, even “sparse” filters with spacing of In (10)/2 and a total of nine coefficients have peak errors of less than 5% for p1:p2 ratios of 10–6 to 106. If a peak error of 5.5% is acceptable, then an even sparser filter with only seven coefficients at a spacing of 3 In (10)/5 may be used. 相似文献
16.
Summary Proof is given of the existence and uniqueness of a discrete, surface, optimum (in the Wiener sense) filter for filtering geophysical fields considering a spectral approach to the construction of the said filter. It is also shown that the filtration coefficients and the filtration errors of the filter, constructed with the help of statistical estimates of the required spectral densities, converge towards their theoretical values. From the mathematical point of view, this paper concludes the building-up of the mathematical model of the discrete Wiener optimum surface filtration, suitable for geophysical fields. 相似文献
17.
基于频域法,利用最新的GOCE卫星重力场模型和卫星测高数据计算了稳态海面动力地形.结合海洋表层漂流浮标的观测结果,对稳态海面动力地形进行了最优空间滤波尺度分析,给出了区域、纬度带和全球稳态海面动力地形的最优空间滤波尺度因子.在此基础上,给出了全球和区域地转流.结果表明:在中高纬度和全球区域,可以分别获得空间尺度优于102km和127km的稳态海面动力地形信息.与海洋表层漂流浮标对比可知,在强流区域,采用稳态海面动力地形得到的地转流速可以解释观测浮标流速的70%;在中高纬度区域,由GOCE重力场得到的地转流略优于对应的GRACE结果;在近赤道区域,由GOCE重力场得到的地转流精度略低于对应的GRACE结果;在北大西洋和阿古拉斯强流区域,由GOCE得到的地转流场明显优于对应的GRACE结果,其精度分别提高了16%和24%. 相似文献
18.
太湖地区农业源污染核算研究进展 总被引：6，自引：4，他引：2
太湖地区水体污染严重,控制太湖污染、改善太湖水质,对发挥太湖的综合服务功能,促进全流域经济可持续发展和社会稳定具有重要的战略意义.充分把握和核算太湖地区农业源污染排放量、入河量及其时空分布是落实太湖地区污染物科学减排从而进一步改善太湖水质的基础.本研究总结和概述了太湖地区农业源污染现有的核算方法、产排污参数、控制途径等,结合太湖地区农业源污染构成、排放与布局的复杂情况,提出综合运用各种核算研究方法,充分考虑各污染来源,确定农业源污染的产污量、排污量和入河量,并通过建立分期、分区、细分类型的农业源污染排放、入河参数,为科学确定农业源污染减排指标以及制定减排措施提供依据. 相似文献
19.
A. H. ELSheikh C. C. Pain F. Fang J. L. M. A. Gomes I. M. Navon 《Stochastic Environmental Research and Risk Assessment (SERRA)》2013,27(4):877-897
A new parameter estimation algorithm based on ensemble Kalman filter (EnKF) is developed. The developed algorithm combined with the proposed problem parametrization offers an efficient parameter estimation method that converges using very small ensembles. The inverse problem is formulated as a sequential data integration problem. Gaussian process regression is used to integrate the prior knowledge (static data). The search space is further parameterized using Karhunen–Loève expansion to build a set of basis functions that spans the search space. Optimal weights of the reduced basis functions are estimated by an iterative regularized EnKF algorithm. The filter is converted to an optimization algorithm by using a pseudo time-stepping technique such that the model output matches the time dependent data. The EnKF Kalman gain matrix is regularized using truncated SVD to filter out noisy correlations. Numerical results show that the proposed algorithm is a promising approach for parameter estimation of subsurface flow models. 相似文献
20.
ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｋｎｏｗｎｔｈａｔｔｈｅｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｉｍｐｅｄａｎｃｅｓｏｆｔｈｅｃｒｕｓｔｍｅｄｉｕｍｍａｙｂｅｃｈａｎｇｅｄｄｕｒｉｎｇｔｈｅｐｒｅｇｎａｎｃｙｏｆａｎｅａｒｔｈｑｕａｋｅ,ａｎｄｍａｙｌｅａｄｔ... 相似文献