单程波最优化可分近似偏移方法精度分析(英文)     

杜淑媛  陈景波  刘洪 《应用地球物理》2011,8(4):328-336

单程波算子的可分近似形式是关于空间变量和波数变量的函数,这使得快速傅里叶变换的利用成为可能,从而极大提高了计算效率。从函数近似的角度来看,最优可分近似方法与其他几种可分近似方法一样都具有单程波算子可分表达的形式,但最优可分近似方法是其中唯一的函数整体近似的方法,这就造成了该方法在相位误差曲线、脉冲响应、模型偏移结果上表现出与其他可分近似方法不一样的特征:该方法精度较高,且随着阶数的增大,对速度变化的敏感度降低。  相似文献   

松辽盆地榆树林油田下白垩统泉头组扶余油层浅水三角洲沉积特征及其演化     

李洋  朱筱敏  宋英琦  刘芬  李超 《高校地质学报》2013,19(1):23-31

基于Fatti近似的弹性阻抗方程及反演   总被引：9，自引：4，他引：5       下载免费PDF全文

王保丽  印兴耀  张繁昌  李爱山 《地球物理学进展》2008,23(1):192-197

用Connolly的弹性阻抗(EI,elastic impedance)公式进行反演只能直接得到纵、横波速度和密度的信息,然后才可间接计算得到纵横波阻抗等其它的参数数据体,这样便增加了一步误差使数据的准确性降低.本文首先针对该方法的这些缺点和不足,提出了以Zoeppritz方程的Fatti近似为基础一种新的弹性阻抗公式,由该公式可得到比用Fatti近似更准确的反射系数,然后对这个公式进行了标准化以实现不同角度的弹性阻抗间量纲的统一,最后用标准化后的公式进行了反演,从反演得到的不同角度的弹性阻抗数据体中可直接提取得到纵横波阻抗数据体.应用实例表明用这种方法提取得到的纵横波阻抗更加稳定、准确,而且能很好地反映储层信息.这种新的方法是对以Connolly公式为基础的传统方法的改进.  相似文献   

二元函数逼近在斜坡可靠度分析中的应用研究     

范雷  唐辉明  胡斌  倪俊 《岩土力学》2008,29(3):624-628

极限平衡分析方法是斜坡稳定性评价中的常用方法,在长期的工程实践中积累了丰富的经验,但其不能考虑斜坡岩土体中实际存在的不确定性,在应用中具有一定的局限性。可靠度分析方法可有效地考虑斜坡系统内的不确定性和相关性,但因状态函数偏导数的求解比较困难,使可靠度分析方法在实际中应用不便。为解决上述问题,根据二元函数插值逼近原理,在矩形区域上构造拉格朗日不完全双二次多项式逼近状态函数,从而近似地计算状态函数的偏导数,求得状态函数的均值和方差,并利用精度较高的一次二阶矩方法来计算斜坡的可靠指标和破坏概率。据鄂西恩施地区马堡营滑坡实例分析表明,引入二元函数插值逼近的一次二阶矩方法计算结果与剩余推力法及Monte-carlo模拟方法结果一致,其精度可满足工程需求。  相似文献   

Comparison of cell- and vertex-centered discretization methods for flow in a two-dimensional discrete-fracture–matrix system     

H. Hægland  A. Assteerawatt  H.K. Dahle  G.T. Eigestad  R. Helmig 《Advances in water resources》2009

Simulations of flow for a discrete-fracture model in fractured porous rocks have gradually become more practical, as a consequence of increased computer power and improved simulation and characterization techniques. Discrete-fracture models can be formulated in a lower-dimensional framework, where the fractures are modeled in a lower dimension than the matrix, or in an equi-dimensional form, where the fractures and the matrix have the same dimension.  相似文献   

The multiresolution character of collocation   总被引：3，自引：0，他引：3  

C. Kotsakis 《Journal of Geodesy》2000,74(3-4):275-290

 An interesting theoretical connection between the statistical (non-stochastic) collocation principle and the multiresolution/wavelet framework of signal approximation is presented. The rapid developments in multiresolution analysis theory over the past few years have provided very useful (theoretical and practical) tools for approximation and spectral studies of irregularly varying signals, thus opening new possibilities for `non-stationary' gravity field modeling. It is demonstrated that the classic multiresolution formalism according to Mallat's pioneering work lies at the very core of some of the general approximation principles traditionally used in physical geodesy problems. In particular, it is shown that the use of a spatio-statistical (non-probabilistic) minimum mean-square-error criterion for optimal linear estimation of deterministic signals, in conjunction with regularly gridded data, always gives rise to a generalized multiresolution analysis in the Hilbert space L <sup>2</sup>(R), under some mild constraints on the spatial covariance function and the power spectrum of the unknown field under consideration. Using the theory and the actual approximation algorithms associated with statistical collocation, a new constructive framework for building generalized multiresolution analyses in L <sup>2</sup>(R) is presented, without the need for the usual dyadic restriction that exists in classic wavelet theory. The multiresolution and `non-stationary' aspects of the statistical collocation approximation procedure are also discussed, and finally some conclusions and recommendations for future work are given. Received: 26 January 1999 / Accepted: 16 August 1999  相似文献   

An enhanced ensemble Kalman filter scheme incorporating model error in sequential coupling between flow and geomechanics     

Edison Caballero  Fernando A. Rochinha  Marcio Borges  Marcio A. Murad 《国际地质力学数值与分析法杂志》2019,43(2):482-500

In this work, we construct a new methodology for enhancing the predictive accuracy of sequential methods for coupling flow and geomechanics while preserving low computational cost. The new computational approach is developed within the framework of the fixed-stress split algorithm procedure in conjunction with data assimilation based on the ensemble Kalman filter (EnKF). In this context, we identify the high-fidelity model with the two-way formulation where additional source term appears in the flow equation containing the time derivative of total mean stress. The iterative scheme is then interlaced with data assimilation steps, which also incorporate the modeling error inherent to the EnKF framework. Such a procedure gives rise to an “enhanced one-way formulation,” exhibiting substantial improvement in accuracy compared with the classical one-way method. The governing equations are discretized by mixed finite elements, and numerical simulation of a 2D slab problem between injection and production wells illustrate the tremendous achievement of the method proposed herein.  相似文献   

Algorithms for Stellar Perturbation Computations on Oort Cloud Comets     

Hans Rickman  Marc Fouchard  Giovanni B. Valsecchi  Christiane Froeschlé 《Earth, Moon, and Planets》2005,97(3-4):411-434

We investigate different approximate methods of computing the perturbations on the orbits of Oort cloud comets caused by passing stars, by checking them against an accurate numerical integration using Everhart’s RA15 code. The scenario under study is the one relevant for long-term simulations of the cloud’s response to a predefined set of stellar passages. Our sample of stellar encounters simulates those experienced by the Solar System currently, but extrapolated over a time of 10¹⁰ years. We measure the errors of perihelion distance perturbations for high-eccentricity orbits introduced by several estimators – including the classical impulse approximation and Dybczyński’s (1994, Celest. Mech. Dynam. Astron. 58, 1330–1338) method – and we study how they depend on the encounter parameters (approach distance and relative velocity). We introduce a sequential variant of Dybczyński’s approach, cutting the encounter into several steps whereby the heliocentric motion of the comet is taken into account. For the scenario at hand this is found to offer an efficient means to obtain accurate results for practically any domain of the parameter space.  相似文献   

On the local multiscale determination of the Earth’s disturbing potential from discrete deflections of the vertical     

T. Fehlinger  W. Freeden  C. Mayer  M. Schreiner 《Computational Geosciences》2008,12(4):473-490

As a first approximation, the Earth is a sphere; as a second approximation, it may be considered an ellipsoid of revolution. The deviations of the actual Earth’s gravity field from the ellipsoidal “normal” field are so small that they can be understood to be linear. The splitting of the Earth’s gravity field into a “normal” and a remaining small “disturbing” field considerably simplifies the problem of its determination. Under the assumption of an ellipsoidal Earth model, high observational accuracy is achievable only if the deviation (deflection of the vertical) of the physical plumb line, to which measurements refer, from the ellipsoidal normal is not ignored. Hence, the determination of the disturbing potential from known deflections of the vertical is a central problem of physical geodesy. In this paper, we propose a new, well-promising method for modelling the disturbing potential locally from the deflections of the vertical. Essential tools are integral formulae on the sphere based on Green’s function with respect to the Beltrami operator. The determination of the disturbing potential from deflections of the vertical is formulated as a multiscale procedure involving scale-dependent regularized versions of the surface gradient of the Green function. The modelling process is based on a multiscale framework by use of locally supported surface curl-free vector wavelets.   相似文献   

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions   总被引：2，自引：0，他引：2  

Chan-Hee Park  Christof Beyer  Sebastian Bauer  Olaf Kolditz 《Environmental Geology》2008,55(8):1755-1766

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.  相似文献