首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
基于单程波偏移算子的地表相关多次波成像   总被引:3,自引:3,他引:0       下载免费PDF全文
在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善.  相似文献   

2.
随着能源和资源勘查开采工作的深入,地形强烈起伏的盆山耦合地区的地震资料处理解释技术正日益成为山地地震勘探面临的重要挑战.逆时偏移方法作为精确的地震偏移成像方法之一,能对地下结构进行高精度成像.逆时偏移的核心是地震波场延拓,由于传统的地震波场延拓技术往往基于水平地表条件,相应的方法在直接处理强地形起伏条件下的地震资料时往往存在一定的精度损失.本文引入一种精度无损的处理起伏边界的模型参数化方法:基于贴体网格的地形"平化"策略发展了与地形有关的地震波波动方程数值模拟方法,采用零延迟归一化互相关成像条件实现了起伏地表条件下的弹性波场逆时偏移成像.对工业界的标准Marmousi模型和盐丘模型进行改造,获得了相应起伏地形条件下的复杂几何模型,开展了起伏地表下的地震偏移成像数值试验.结果表明基于贴体网格"平化"策略的逆时偏移成像方法具有较高的灵活性,可适应不同类型起伏地表采集的地震资料,显示出该方法在地震勘探领域的良好应用前景.  相似文献   

3.
逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求,井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异,本文将相应的变网格算法应用到偏移中。另外,也将Lanczos滤波引入到波场外推计算中,从而可以很好地消除变网格引起的人为反射。试算表明,对比地面、井中和井间的偏移,可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。  相似文献   

4.
为更好地研究复杂构造和速度分布条件下地震资料的精确偏移.在波场模拟方面,本文采用了高阶交错网格有限差分[1]推导了二维各向同性介质声波方程的数值模拟公式.并给出了声波方程完全匹配层吸收边界[2-5](PML).在偏移方面,本文介绍了一种能适应介质速度横向变化的非稳态相移算子(NSPS)[6]及其叠后深度偏移方法.并对简单的四层模型和Marmousi模型进行试算.偏移结果表明,该方法可以有效地处理复杂条件下的各向同性介质声波偏移问题.  相似文献   

5.
王晓  白志明  余丹  纪寿文 《中国地震》2020,36(2):350-358
人工源宽角反射(折射)地震资料具有偏移距较大、信噪比较低等特点,通常用于地震波走时反演重建地壳速度结构。逆时偏移成像方法作为勘探地震学领域获取地下构造形态的有效手段之一,可以有效弥补走时反演方法的不足。本文针对大偏移距宽角反射(折射)地震实验,利用四边形网格谱元法进行波场模拟,结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度获取模型合成地震记录,后采用逆时偏移成像方法将合成地震记录偏移归位,获取地壳几何结构,验证了逆时偏移成像方法在宽角地震资料处理及结果解释中的适用性,为后期实际地震资料的偏移成像提供了理论依据和支持。  相似文献   

6.
复杂构造地震波场分析(英文)   总被引:2,自引:0,他引:2  
在我国的西部地区,地震波场十分复杂,信噪比低。本文采用波动方程正演模拟的方法研究复杂波场的形成原因。在模拟物性差异较大介质中的地震波场时,密度的影响不可忽略,因此,本文用含密度项的声波方程的交错网格有限差分法模拟地震波场并进行分析。设计了一个具有起伏地表、低速覆盖层和高速地层出露的复杂构造,从瞬时波场分布分析了形成复杂波场的原因。低速层对地震波场的影响明显,低速层中产生很强的槽波,低速层顶底形成的多次反射向地下传播又形成了复杂的反射波场。为了验证波场模拟结果的可靠性,对模拟波场用与正演模拟不同的算法进行叠前深度偏移,得到了与已知构造相同的偏移叠加剖面。  相似文献   

7.
多震源地震数据偏移成像方法   总被引:1,自引:0,他引:1       下载免费PDF全文
多震源地震技术是一种高效的地震数据采集方法技术,得到的地震记录是来自多个震源的混合地震数据.本文在多震源波场传播理论和地震波场满足线性叠加原理的基础上,提出了两种多震源地震数据的偏移成像方法.第一种方法是首先对多震源地震数据进行分离,得到各个单震源的地震数据,然后再利用常规的偏移成像方法进行处理;第二种方法是多震源地震数据的直接偏移成像.把本文提出的多震源偏移成像方法应用于数值模拟的多震源地震数据,验证了本文方法的正确性和有效性,直接偏移成像方法较分离后再偏移方法具有更高的计算效率.  相似文献   

8.
相对于仅仅提供井中菲涅耳带内地层信息的井口零偏移距VSP(vertical seismic profile)勘探,有偏移距VSP勘探可以得到远离井区域的反射信息,使研究区域拓宽,成像精度提高,且便于对地面勘探实施约束.提出了一种利用有偏移距VSP数据研究地球内部结构的CDP(common depth point )成图方法.其中,利用F-K波场分离技术进行VSP资料的上/下波场分离,并在均匀介质环境下计算成图条件,进而进行地球内部结构的图象重建. 利用本方法对有限差分法正演模拟的超深井VSP地震资料进行了VSP成图处理,并获得较为满意的结果.模拟结果显示,在模型速度与真实速度误差小于10%的情况下,利用本方法仍可以获得比较准确的地球内部结构图象.最后,介绍了本方法处理实际有偏移距VSP观测资料的结果.   相似文献   

9.
介绍一种基于离散粒子理论地震波传播数值模拟的网格剖分计算方法.根据离散粒子理论,将研究区域划分为由一系列相互作用的粒子组成的正六边形网格,这些粒子在它们的接触点处发生相互作用,并用Hooke定律和Newton定律描述.为解决六边形网格带来的网格交错而难以计算以及波场输出问题,将横向网格进行加密,加密处赋予假想的粒子,输出波场时选取偶数行偶数列点或奇数行奇数列点的波场值.均匀介质和层状介质模型的数值模拟结果表明,该网格剖分计算方法能够将离散粒子理论用于模拟弹性波在非均匀各向同性介质中地震波的传播.  相似文献   

10.
起伏地表条件下基于复Pade逼近的叠前深度偏移   总被引:1,自引:0,他引:1  
叠前深度偏移是解决复杂地表和复杂构造地震成像的有效技术,而波场“直接下延”法实现了复杂地表条件下的地震成像.基于上述成果,结合高精度的波场延拓算子,本文提出了一种新的叠前深度偏移方法,这种方法是在波场延拓时,对声波方程中的平方根项进行复Pade逼近,通过推导得到基于复Pade逼近的傅里叶有限差分算子,结合波场“直接下延”法,实现了起伏地表条件下的叠前深度偏移,该算法减少了偏移噪音,从而得到准确、稳定的偏移成像结果.通过理论模型试算和实际资料试处理,验证了该方法的有效性.  相似文献   

11.
Recent theoretical and numerical models for the motion of saltating particles close to stream beds are constituted of three sub-models: a) a set of equations describing the particle “free flight”, b) a sub-model to calculate the post-collision particle velocity, and c) a mathematical representation of the bed roughness. In this paper, a comprehensive three-dimensional (3-D), theoretical/numerical model for bedload motion at large Reynolds numbers is presented. By using geometric considerations and stochastic parameters to characterize collisions with the wall, five new sub-models for representation of bed roughness are, for the first time to the best of our knowledge, proposed and implemented. The emphasis of this paper is on the particle model, for which Basset, Magnus, drag, submerged weight, virtual mass, and lift forces are included. For the range of particle sizes (sands) analyzed herein, it is found that the stream-wise contribution of the Basset force, compared to other forces, may be as large as 60%. Whereas in the wall-normal direction, the Basset force is equally important as the drag force, and it is exceeded only by the submerged weight. It is also found that the best agreement between numerical and experimental results in terms of jump length, jump height, and stream-wise particle velocity is achieved for restitution and friction coefficients of 0.65 and 0.1, respectively. Important conclusions are obtained regarding the lack of realistic prediction with available “roughness” models with small ranges of angles.  相似文献   

12.
Tidal bores may appear in some estuaries when the tides quickly reach a high level. This phenomenon is rare but has a strong impact during its short duration: i.e. the river bed is significantly eroded and sediments are then transported. In this paper, the trajectories of suspended particles induced by this flow are numerically studied. Four undular bores with Froude numbers between 1.1 and 1.2 are studied. Despite similar Froude numbers, various initial flow conditions were selected to produce or not an inversion of the flow direction during the bore passage. The particle trajectories associated with each distinct flow configuration are presented and analyzed. These trajectories, estimated by solving the Maxey-Riley equation, appear to be very different even though the Froude numbers of flows are similar. These observations are important because the Froude number is often used to characterize a tidal bore as it describes well the free surface, however, it cannot describe the sediment transport. Finally, Chen's model of wave-current interactions is adapted to fit the cases studied and is applied to the four bores simulated. The results highlight that this latter model can reproduce the observed trajectories and dissociate their different components. From this model, it is shown that the inertial and Basset history effects can be neglected compared to the gravity and flow entrainment effects due to the viscous drag when one wants to determine the long-term trajectories of suspended particles.  相似文献   

13.
Particle models are often used to simulate the spreading of a pollutant in coastal waters in case of a calamity at sea. Here many different particle tracks starting at the point of release are generated to determine the particle concentration at some critical locations. This Monte Carlo method, however, consumes a large CPU time. Recently, Milstein, et al. [Milstein, GN, Schoenmakers JGM, Spokoiny, V. Transition density estimation for stochastic differential equations via forward–reverse representations. Bernoulli 2004:10(2);281–312] introduced the concept of reverse-time diffusion and derived a reverse particle model from the original forward simulation model. While the original forward model provides results on where the pollutant will go to, the reverse system gives information about where the pollutant came from. The Monte Carlo estimator for the particle concentration can also be based on realizations of this reverse system. In this paper we apply this concept to estimate particle concentrations in coastal waters. The results of the experiments show that the CPU time compared with the classical Monte Carlo method is reduced at least order of magnitude.  相似文献   

14.
《国际泥沙研究》2020,35(1):1-14
A two-dimensional(2 D) computational model for simulation of incipient sediment motion for noncohesive uniform and non-uniform particles on a horizontal bed was developed using the Discrete Element Method(DEM).The model was calibrated and verified using various experimental data reported in the literature and compared with different theories of incipient particle motion.Sensitivity analysis was done and the effects of relevant parameters were determined.In addition to hydrodynamic forces such as drag,shear lift and Magnus force,the particle-particle interaction effects were included in the model.The asymptotic critical mobility number was evaluated for various critical particle Reynolds numbers(R*) in the range of very small and very large R*.The obtained curve is classified into four regions.It was found that in the linear region,the drag force has the principal role on the initiation of motion.Moreover,the critical mobility number is independent of particle diameter.A procedure for estimating the critical shear velocity directly from the information on particle diameter and roughness height was developed.Finally,the mechanism of incipient motion for the different regions was studied and the effect of different forces on the incipient particle motion was obtained.It was found that the maximum effects of lift and Magnus forces were,respectively,less than ten and twenty percent of the total force.The drag force,however,was typically the dominant force accounting for majority of the net hydrodynamic force acting on sediment particles at the onset of incipient motion.  相似文献   

15.
Plume containment using pump-and-treat (PAT) technology continues to be a popular remediation technique for sites with extensive groundwater contamination. As such, optimization of PAT systems, where cost is minimized subject to various remediation constraints, is the focus of an important and growing body of research. While previous pump-and-treat optimization (PATO) studies have used discretized (finite element or finite difference) flow models, the present study examines the use of analytic element method (AEM) flow models. In a series of numerical experiments, two PATO problems adapted from the literature are optimized using a multi-algorithmic optimization software package coupled with an AEM flow model. The experiments apply several different optimization algorithms and explore the use of various pump-and-treat cost and constraint formulations. The results demonstrate that AEM models can be used to optimize the number, locations and pumping rates of wells in a pump-and-treat containment system. Furthermore, the results illustrate that a total outflux constraint placed along the plume boundary can be used to enforce plume containment. Such constraints are shown to be efficient and reliable alternatives to conventional particle tracking and gradient control techniques. Finally, the particle swarm optimization (PSO) technique is identified as an effective algorithm for solving pump-and-treat optimization problems. A parallel version of the PSO algorithm is shown to have linear speedup, suggesting that the algorithm is suitable for application to problems that are computationally demanding and involve large numbers of wells.  相似文献   

16.
Snow water equivalent prediction using Bayesian data assimilation methods   总被引:1,自引:0,他引:1  
Using the U.S. National Weather Service’s SNOW-17 model, this study compares common sequential data assimilation methods, the ensemble Kalman filter (EnKF), the ensemble square root filter (EnSRF), and four variants of the particle filter (PF), to predict seasonal snow water equivalent (SWE) within a small watershed near Lake Tahoe, California. In addition to SWE estimation, the various data assimilation methods are used to estimate five of the most sensitive parameters of SNOW-17 by allowing them to evolve with the dynamical system. Unlike Kalman filters, particle filters do not require Gaussian assumptions for the posterior distribution of the state variables. However, the likelihood function used to scale particle weights is often assumed to be Gaussian. This study evaluates the use of an empirical cumulative distribution function (ECDF) based on the Kaplan–Meier survival probability method to compute particle weights. These weights are then used in different particle filter resampling schemes. Detailed analyses are conducted for synthetic and real data assimilation and an assessment of the procedures is made. The results suggest that the particle filter, especially the empirical likelihood variant, is superior to the ensemble Kalman filter based methods for predicting model states, as well as model parameters.  相似文献   

17.
A 3D Lagrangian model of the saltation of solid spherical particles on the bed of an open channel flow, accounting for turbulence-induced mechanisms, is proposed and employed as the key tool of the study. The differences between conventional 2D models and a proposed 3D saltation model are discussed and the advantages of the 3D model are highlighted. Particularly, the 3D model includes a special procedure allowing generation of 3D flow velocity fields. This procedure is based on the assumption that the spectra of streamwise, vertical and transverse velocity components are known at any distance from the bed. The 3D model was used to identify and quantify effects of turbulence on particle entrainment and saltation. The analysis of particle trajectories focused on their diffusive nature, clarifying: (i) the effect of particle mobility parameter; (ii) the effect of bed topography; and (iii) the effect of turbulence. Specifically, the results of numerical simulations describing the abovementioned effects on the change in time of the variance are presented. In addition, the change in time of the skewness and kurtosis, which are likely to reflect the turbulence influence on the spread of particles, are also shown. Two different diffusion regimes (local and intermediate) for each of the investigated flow conditions are confidently identified.  相似文献   

18.
Modeling of suspended sediment particle movement in surface water can be achieved by stochastic particle tracking model approaches.In this paper,different mathematical forms of particle tracking models are introduced to describe particle movement under various flow conditions,i.e.,the stochastic diffusion process,stochastic jump process,and stochastic jump diffusion process.While the stochastic diffusion process can be used to represent the stochastic movement of suspended particles in turbulent flows,the stochastic jump and the stochastic jump diffusion processes can be used to describe suspended particle movement in the occurrences of a sequence of extreme flows.An extreme flow herein is defined as a hydrologic flow event or a hydrodynamic flow phenomenon with a low probability of occurrence and a high impact on its ambient flow environment.In this paper,the suspended sediment particle is assumed to immediately follow the extreme flows in the jump process(i.e.the time lag between the flow particle and the sediment particle in extreme flows is considered negligible).In the proposed particle tracking models,a random term mainly caused by fluid eddy motions is modeled as a Wiener process,while the random occurrences of a sequence of extreme flows can be modeled as a Poisson process.The frequency of occurrence of the extreme flows in the proposed particle tracking model can be explicitly accounted for by the Poisson process when evaluating particle movement.The ensemble mean and variance of particle trajectory can be obtained from the proposed stochastic models via simulations.The ensemble mean and variance of particle velocity are verified with available data.Applicability of the proposed stochastic particle tracking models for sediment transport modeling is also discussed.  相似文献   

19.
In this study, a three-dimensional particle tracking model coupled to a terrain following ocean model is used to investigate the dispersion and the deposition of fish farm particulate matter (uneaten food and fish faeces) on the seabed due to tidal currents. The particle tracking model uses the computed local flow field for advection of the particles and random movement to simulate the turbulent diffusion. Each particle is given a settling velocity which may be drawn from a probability distribution according to settling velocity measurements of faecal and feed pellets. The results show that the maximum concentration of organic waste for fast sinking particles is found under the fish cage and continue monotonically decreasing away from the cage area. The maximum can split into two maximum peaks located at both sides of the centre of the fish cage area in the current direction. This process depends on the sinking time (time needed for a particle to settle at the bottom), the tidal velocity and the fish cage size. If the sinking time is close to a multiple of the tidal period, the maximum concentration point will be under the fish cage irrespective of the tide strength. This is due to the nature of the tidal current first propagating the particles away and then bringing them back when the tide reverses. Increasing the cage size increases the likelihood for a maximum waste accumulation beneath the fish farm, and larger farms usually means larger biomasses which can make the local pollution even more severe. The model is validated by using an analytical model which uses an exact harmonic representation of the tidal current, and the results show an excellent agreement. This study shows that the coupled ocean and particle model can be used in more realistic applications to help estimating the local environmental impact due to fish farms.  相似文献   

20.
X. Qie  S. Soula  S. Chauzy 《Annales Geophysicae》1994,12(12):1218-1228
A numerical model called PICASSO [Production d’Ions Corona Au Sol Sous Orage (French) and Production of Corona Ions at the Ground Beneath Thundercloud (English)], previously designed, is used to describe the evolution of the principal electrical parameters below a thunderstorm, taking into account the major part played by corona ions. In order to improve the model restitution of a real situation, various improvements are performed: an initial vertical distribution of aerosol particles is introduced instead of the previously used uniform concentration; time and space calculation steps are adjusted according to the electric field variation rate; the upper boundary condition is improved; and the coefficients of ion attachment are reconsidered with an exhaustive bibliographic study. The influence of the ion attachment on aerosol particles, on the electric field and charge density aloft, is studied by using three different initial aerosol particle concentrations at ground level and two types of initial vertical distributions: uniform and non-uniform. The comparison between field data and model results leads to adjust the initial aerosol particle concentration over the experimental site at the value of 5000 cm−3 which appears to be highly realistic. The evolutions of the electric field and of the charge density at altitude are greatly influenced by the aerosol concentration. On the contrary, the surface intrinsic field, defined as the electric field that would exist underneath a thundercloud if there were no local charges, is weakly affected when the model is forced by the surface field. A good correlation appears between the success in the triggered lightning attempts and this intrinsic field evaluation. Therefore, when only the surface field is available, the model can be used in a triggered lightning experiment.  相似文献   

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

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