首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 359 毫秒
1.
一种改进的基于非高斯性最大化的预测反褶积算法   总被引：3，自引：1，他引：2  
The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.  相似文献   

2.
Freeze-thaw processes in soils,including changes in frost and thaw fronts(FTFs),are important physical processes.The movement of FTFs affects soil hydrothermal characteristics,as well as energy and water exchanges between the land surface and the atmosphere and hydrothermal processes in the land surface.This paper reduces the issue of soil freezing and thawing to a multiple moving-boundary problem and develops a soil water and heat transfer model which considers the effects of FTF on soil hydrothermal processes.A local adaptive variable-grid method is used to discretize the model.Sensitivity tests based on the hierarchical structure of the Community Land Model(CLM)show that multiple FTFs can be continuously tracked,which overcomes the difficulties of isotherms that cannot simultaneously simulate multiple FTFs in the same soil layer.The local adaptive variable-grid method is stable and offers computational efficiency several times greater than the high-resolution case.The simulated FTF depths,soil temperatures,and soil moisture values fit well with the observed data,which further demonstrates the potential application of this simulation to the land-surface process model.  相似文献   

3.
Based on a detailed analysis of differences between seismic data and well logs, we discuss the problem of matching seismic traces and well logs and present a new matching method based on event search in instantaneous phase which greatly improves seismic resolution. The method is based on flattening events in instantaneous phase to compare the seismic traces to the well log traces with the same phase. We calculate the coefficients using the singular value decomposition method to extrapolate the well logs. As a result, the events in the seismic profile are continuous and match well with well logs. We apply this method to the Mao-2 well in Daqing Oilfield with good results.  相似文献   

4.
The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic(TI)half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf(1985),and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.  相似文献   

5.
In this paper,an approach is developed to optimize the quality of the training samples in the conventional Artificial Neural Network(ANN)by incorporating expert knowledge in the means of constructing expert-rule samples from rules in an expert system,and through training by using these samples,an ANN based on expert-knowledge is further developed.The method is introduced into the field of quantitative identification of potential seismic sources on the basis of the rules in an expert system.Then it is applied to the quantitative identification of the potential seismic sources in Beijing and its adjacent area.The result indicates that the expert rule based on ANN method can well incorporate and represent the expert knowledge in the rules in an expert system,and the quality of the samples and the efficiency of training and the accuracy of the result are optimized.  相似文献   

6.
The theory of elastic wave scattering is a fundamental concept in the study of elastic dynamics and wave motion,and the wave function expansion technique has been widely used in many subjects.To supply the essential tools for solving wave scattering problems induced by an eccentric source or multi-sources as well as multi-scatters,a whole-space transform formula of cylindrical wave functions is presented and its applicability to some simple cases is demonstrated in this study.The transforms of wave functions in cylindrical coordinates can be classifi ed into two basic types: interior transform and exterior transform,and the existing Graf’s addition theorem is only suitable for the former.By performing a new replacement between the two coordinates,the exterior transform formula is fi rst deduced.It is then combined with Graf’s addition theorem to establish a whole-space transform formula.By using the whole-space transform formula,the scattering solutions by the sources outside and inside a cylindrical cavity are constructed as examples of its application.The effectiveness and advantages of the whole-space transform formula is illustrated by comparison with the approximate model based on a large cycle method.The whole-space transform formula presented herein can be used to perform the transform between two different cylindrical coordinates in the whole space.In addition,its concept and principle are universal and can be further extended to establish the coordinate transform formula of wave functions in other coordinate systems.  相似文献   

7.
Because zero-offset VSP (Vertical Seismic Profile) data can only provide the information of rock properties and structure in the area around the Fresnel zone within the well, the scheme of VSP with offset was developed to acquire the reflection information away from the borehole in order to widen the range of VSP survey and to improve the precision of imaging.In this paper, we present a new CDP (Common Depth Point) mapping approach to image the reflecting structure by using offset VSP data. For the processing of offset VSP data, we firstly separated the up-going and down-going wave-fields from VSP data by means of F-K filtering technique, and we can calculate the mapping conditions (position and reflecting traveltime for CDP point) in homogeneous media, and then reconstruct the inner structure of the earth.This method is tested by using the offset VSP data which are used to simulate the case of super-deep borehole by means of finite-difference method. The imaged structure matches the real model very well. The results show that the method present here could accurately image the inner structure of the earth if the deviation of initial velocity model from the true model is less than 10%. Finally, we presented the imaged results for the real offset data by using this method.  相似文献   

8.
The studies on fluvial reservoir architecture are mainly aimed at outcrop and modern deposition,but rarely at the subsurface reservoir,so there are few effective methods to predict the distribution of subsurface reservoir architectures. In this paper,taking the meandering river reservoir of Guantao formation Gudao Oilfield,Jiyang depression,Baohai Gulf Basin,East China as an example,the archi-tectural modeling method of complex meandering belt reservoir is proposed,that is hierarchy con-straint,pattern fitting and multi-dimensional interaction. Architectures of meandering river reservoir can be divided into three hierarchies: meandering channel sandbody,point bar and lateral accretion body. Different hierarchies of the quantitative architecture pattern are fitted to subsurface well data (including dynamic monitoring data) in different hierarchies through one-dimensional hole,2D profiles and plane and 3D space,which are verified by each other. And then 3D model in different hierarchies is established. At the same time,the quantificational relationship between width of active river and the scale of point bar is set up,and the scale of lateral accretion sand body and shale beddings is con-firmed with horizontal well data. The study not only has significant meaning on the development of geology,but also can improve the oilfield exploitation greatly.  相似文献   

9.
The conventional nonstationary convolutional model assumes that the seismic signal is recorded at normal incidence. Raw shot gathers are far from this assumption because of the effects of offsets. Because of such problems, we propose a novel prestack nonstationary deconvolution approach. We introduce the radial trace （RT） transform to the nonstationary deconvolution, we estimate the nonstationary deconvolution factor with hyperbolic smoothing based on variable-step sampling （VSS） in the RT domain, and we obtain the high-resolution prestack nonstationary deconvolution data. The RT transform maps the shot record from the offset and traveltime coordinates to those of apparent velocity and traveltime. The ray paths of the traces in the RT better satisfy the assumptions of the convolutional model. The proposed method combines the advantages of stationary deconvolution and inverse Q filtering, without prior information for Q. The nonstationary deconvolution in the RT domain is more suitable than that in the space-time （XT） domain for prestack data because it is the generalized extension of normal incidence. Tests with synthetic and real data demonstrate that the proposed method is more effective in compensating for large-offset and deep data.  相似文献   

10.
We present numerical modeling of SH-wave propagation for the recently proposed whole Moon model and try to improve our understanding of lunar seismic wave propagation. We use a hybrid PSM/FDM method on staggered grids to solve the wave equations and implement the calculation on a parallel PC cluster to improve the computing efficiency. Features of global SH-wave propagation are firstly discussed for a 100-km shallow and900-km deep moonquakes, respectively. Effects of frequency range and lateral variation of crust thickness are then investigated with various models. Our synthetic waveforms are finally compared with observed Apollo data to show the features of wave propagation that were produced by our model and those not reproduced by our models. Our numerical modeling show that the low-velocity upper crust plays significant role in the development of reverberating wave trains. Increasing frequency enhances the strength and duration of the reverberations.Surface multiples dominate wavefields for shallow event.Core–mantle reflections can be clearly identified for deep event at low frequency. The layered whole Moon model and the low-velocity upper crust produce the reverberating wave trains following each phases consistent with observation. However, more realistic Moon model should be considered in order to explain the strong and slow decay scattering between various phases shown on observation data.  相似文献   

11.
In acoustic logging-while-drilling （ALWD） finite difference in time domain （FDTD） simulations, large drill collar occupies, most of the fluid-filled borehole and divides the borehole fluid into two thin fluid columns （radius -27 mm）. Fine grids and large computational models are required to model the thin fluid region between the tool and the formation. As a result, small time step and more iterations are needed, which increases the cumulative numerical error. Furthermore, due to high impedance contrast between the drill collar and fluid in the borehole （the difference is 〉30 times）, the stability and efficiency of the perfectly matched layer （PML） scheme is critical to simulate complicated wave modes accurately. In this paper, we compared four different PML implementations in a staggered grid finite difference in time domain （FDTD） in the ALWD simulation, including field-splitting PML （SPML）, multiaxial PML（M- PML）, non-splitting PML （NPML）, and complex frequency-shifted PML （CFS-PML）. The comparison indicated that NPML and CFS-PML can absorb the guided wave reflection from the computational boundaries more efficiently than SPML and M-PML. For large simulation time, SPML, M-PML, and NPML are numerically unstable. However, the stability of M-PML can be improved further to some extent. Based on the analysis, we proposed that the CFS-PML method is used in FDTD to eliminate the numerical instability and to improve the efficiency of absorption in the PML layers for LWD modeling. The optimal values of CFS-PML parameters in the LWD simulation were investigated based on thousands of 3D simulations. For typical LWD cases, the best maximum value of the quadratic damping profile was obtained using one do. The optimal parameter space for the maximum value of the linear frequency-shifted factor （a0） and the scaling factor （β0） depended on the thickness of the PML layer. For typical formations, if the PML thickness is 10 grid points, the global error can be reduced to 〈1% using the optimal PML parameters, and the error will decrease as the PML thickness increases.  相似文献   

12.
Seismic inversion is a highly ill-posed problem, due to many factors such as the limited seismic frequency bandwidth and inappropriate forward modeling. To obtain a unique solution, some smoothing constraints, e.g., the Tikhonov regularization are usually applied. The Tikhonov method can maintain a global smooth solution, but cause a fuzzy structure edge. In this paper we use Huber-Markov random-field edge protection method in the procedure of inverting three parameters, P-velocity, S-velocity and density. The method can avoid blurring the structure edge and resist noise. For the parameter to be inverted, the Huber- Markov random-field constructs a neighborhood system, which further acts as the vertical and lateral constraints. We use a quadratic Huber edge penalty function within the layer to suppress noise and a linear one on the edges to avoid a fuzzy result. The effectiveness of our method is proved by inverting the synthetic data without and with noises. The relationship between the adopted constraints and the inversion results is analyzed as well.  相似文献   

13.
The overexploitation of underground water leads to the continuous drawdown of groundwater levels, change of water quality and dry-up in dynamic water level observation wells. Due to land subsidence, the well pipes uplift and the observation piping systems are damaged. These environmental geology problems can present serious difficulties for the identification of earthquake anomalies by groundwater level observation. Basied on hydrogeological theories and methods, the paper analyzes the relations of the water balance state of aquifers with stress-strain conditions and the water level regime, and then discusses preliminarily the theory and method for identifying well water level anomalies in a groundwater overdraft area. The result shows that we can accurately judge the nature of the anomaly according to the diffusion character of the drawdown funnel in the well area in combination with the aforementioned theory and method and multi-year variation patterns obtained from existing data. The results of the research are helpful for distinguishing the influence of single centralized water pumping from the long-term overdraft of water on the water level, correctly recognizing water level anomalies in the groundwater overdraft area and increasing the level of earthquake analysis and prediction.  相似文献   

14.
The spectral representation method(SRM) is widely used to simulate spatially varying ground motions.This study focuses on the approximation approach to the SRM based on root decomposition,which can improve the efficiency of the simulation.The accuracy of the approximation approach may be affected by three factors: matrix for decomposition,distribution of frequency interpolation nodes and elements for interpolation.The influence of these factors on the accuracy of this approach is examined and the following conclusions are drawn.The SRM based on the root decomposition of the lagged coherency matrix exhibits greater accuracy than the SRM based on the root decomposition of the cross spectral matrix.The equal energy distribution of frequency interpolation nodes proposed in this study is more effective than the counter pith with an equal spacing.Elements for interpolation do not have much of an effect on the accuracy,so interpolation of the elements of the decomposed matrix is recommended because it is less complicated from a computational efficiency perspective.  相似文献   

15.
We propose a domain decomposition method based on the spectral element method(DDM-SEM)for elastic wave computation in frequency domain.It combines the high accuracy of the spectral element method and the high degree of parallelism of a domain decomposition technique,which makes this method suitable for accurate and efficient simulations of large scale problems in elastodynamics.In the DDM-SEM,the original large-scale problem is divided into a number of well designed subdomains.We use the spectral element method independently for each subdomain,and the neighboring subdomains are connected by a frequency-domain version of Riemann transmission condition(RTC)for elastic waves.For the proposed method,we can employ the non-conforming meshes and different interpolation orders in different subdomains to maximize the efficiency.By separating the internal and boundary unknowns of each subdomain,an efficient and naturally parallelizable block LDU direct solver is developed to solve the final system matrix.Numerical experiments verify its accuracy and efficiency,and show that the proposed DDM-SEM can be a promising numerical tool for accurately and effectively solving large and multi-scale problems of elastic waves.It is potentially valuable for the frequency domain seismic inversion where multiple source illuminations are required.  相似文献   

16.
Weak Seismic Signal Extraction Based on the Curvelet Transform   总被引：1，自引：1，他引：0  
Seismic signal denoising is a key step in seismic data processing. Airgun signals are easy to be interfered with by noise when it travels a long distance due to the weak energy of active source signal of the airgun. Aiming to solve this problem, and considering that the conventional Curvelet transform threshold processing method does not use the seismic spectrum information, we independently process the Curvelet scale layer corresponding to valid data based on the characteristics of the Curvelet transform of multi-scale, multi-direction and capable of expressing the sparse seismic signals in order to fully excavate the information features. Combined with the Curvelet adaptive threshold denoising the algorithm, we apply the Curvelet transform to denoising seismic signals while retaining the weak information in the signal as much as possible. The simulation experiments show that the improved threshold denoising method based on Curvelet transform is superior to the frequency domain filtering, wavelet denoising and traditional Curvelet denoising method in detailed information extraction and signal denoising of low SNR signals. The calculation accuracy of the relative wave velocity variation of underground medium is improved.  相似文献   

17.
Based on the Gauss linear frequency-modulated wavelet transform, a new characteristic index is presented, namely time-frequency energy attenuation factor which can reflect the difference features of waveform in earthquake focus mechanism, wave traveling path and its attenuation characteristics in focal area or near field. In order to test its validity, we select the natural earthquakes and explosion or collapse events whose focus mechanisms vary obviously, and some natural earthquakes located at the same site or in a very small area. The study indicates that the time-frequency energy attenuation factors of the natural earthquakes are obviously different with that of explosion or collapse events, and the change of the time-frequency energy attenuation factors is relatively stable for the earthquakes under the normal seismicity background. Using the above-mentioned method, it is expected to offer a useful criterion for strong earthquake prediction by continuous earthquake observation.  相似文献   

18.
The full operator method (FOM) has been proposed to overcome some of the shortcomings of the commonly used operator splitting method (OSM). In particular, the FOM is improved by increasing the accuracy of both the predictor and corrector using the estimated tangent stiffness of the tested structure. The numerical characteristics of the FOM, including stability and accuracy, are investigated in this study. It is shown that FOM is conditionally stable. The stability and accuracy characteristics are dependent on the accuracy of the estimated tangent stiffness and the parameters associated with the acceleration variation in the time-stepping integration method. Mass-spring systems with different types of nonlinearity, including hardening, stiffening, and softening behavior, are used to evaluate the performance of the FOM. It is found that the FOM can capture these types of nonlinearity with satisfactory accuracy. Using a prototype 12-story composite coupled wall system, the influences of the strong nonlinearity of the system as well as the displacement control errors from hydraulic actuators on the performance of the FOM are explored. The results show that the FOM is capable of generating reasonably accurate results despite the presence of strong structural nonlinearity and displacement control errors.  相似文献   

19.
The method of fuzzy mathematics for simultaneous assessment of time and intensity of earthquake hazards has been studied.This method is based on fundamental statistical indices of regional seismicity.Applying the retrieval method of fuzzy information,we can classify the time and intensity into several intervals and classes of seismic activity,then the possible time interval of large earthquakes with magnitude of M≥Ms can be estimated in a given region.Based on the preceding idea,an FRPP program is constructed.For the automatic data processing when this method is used,it is very important to design the statistical process of each index decomposition so that the program could be fit to a different sample discussed.There are some functions in the FRPP program.The man-made impact on results is reduced to the minimum as far as possible.Computation time is saved.There is a menu on which time interval,index,intensity class,and output data all can be selected.The catalog input that can be displayed on the scre  相似文献   

20.
Foreland basins on the western edge of Ordos are characterized by typical loess geomorphologic features. Its surface is severely cut, and its underground geological structures are complicated.Making use of the advantage of CEMP in exploration in frontland as well as successful multi-year exploration experiences in West China, we first employed the CEMP method to conduct exploration on the western edge of the Ordos basin and its surrounding areas. The practices proved that rich and reliable geological results have been achieved with the CEMP method in complicated areas covered with loess. This can provide important evidence for researching the contact relationship among basins and mountains, structural feature of target layers and predict favourable areas.  相似文献   

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

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