共查询到20条相似文献，搜索用时 46 毫秒
1.
Conventional two‐way splitting Fourier finite‐difference migration for 3D complex media yields azimuthal anisotropy where an additional phase correction is needed with much increase of computational cost. We incorporate the alternating‐direction‐implicit plus interpolation scheme into the conventional Fourier finite‐difference method to reduce azimuthal anisotropy. This scheme retains the high‐order remnants ignored by the two‐way splitting in the form of a wavefield interpolation in the wavenumber domain. The wavefield interpolation for each step of downward extrapolation is implemented between the wavefields before and after the conventional Fourier finite‐difference extrapolation. As the Fourier finite‐difference migration is implemented in the space and wavenumber dual space, the Fourier transforms between space and wavenumber domain that were needed for the alternating‐direction‐implicit plus interpolation in frequency domain (FD) migration are saved in Fourier finite‐difference migration. Since the azimuth anisotropy in Fourier finite‐difference is much less than that in FD, the application of the alternating‐direction‐implicit plus interpolation scheme in Fourier finite‐difference migration is superior to that in FD migration in handling complex media with large velocity contrasts and steep dips. Impulse responses show that the presented method reduces the azimuthal anisotropy at almost no extra cost. 相似文献
2.
K.R. Rushton 《Journal of Hydrology》1973,18(1):1-19
When digital computers are used in the analysis of time variant groundwater flow in an aquifer, the time is divided into discrete intervals. These discrete time steps can cause oscillations or other errors in the solutions particularly when wells start and stop pumping. A large number of solutions are obtained with a wide range of time steps using three implicit methods, the Crank-Nicholson, the alternating direction implicit and the backward difference methods, and, by comparing the results with theoretical values, recommendations are made to aid the choice of the optimum time interval. 相似文献
3.
K.R. Rushton 《Journal of Hydrology》1974,21(2):153-172
The alternating direction implicit method, when used for aquifer analysis, is liable to lead to erroneous results. By a careful examination of the effect of various factors, such as external boundaries, variable mesh intervals, internal boundaries with fixed heads, re-entrant boundaries and leakage, the occasions when these errors are serious have been identified. Techniques are described which can avoid these errors and lead to reliable numerical solutions. 相似文献
4.
生产油井井下温度场数值模拟分析 总被引：1，自引：0，他引：1
An improved numerical simulation method is presented to calculate the downhole temperature distribution for multiple pay zones in producing oil wells. Based on hydrodynamics and heat transfer theory, a 2-D temperature field model in cylindrical coordinates is developed. In the model, we considered general heat conduction as well as the heat convection due to fluid flow from porous formation to the borehole. We also take into account the fluid velocity variation in the wellbore due to multiple pay zones. We present coupled boundary conditions at the interfaces between the wellbore and adjacent formation, the wellbore and pay zone, and the pay zone and adjacent formation. Finally, an alternating direction implicit difference method （ADI） is used to solve the temperature model for the downhole temperature distribution. The comparison of modeled temperature curve with actual temperature log indicates that simulation result is in general quite similar to the actual temperature log. We found that the total production rate, production time, porosity, thickness of pay zones, and geothermal gradient, all have effects on the downhole temperature distribution. 相似文献
5.
E. T. Vogler C. V. Chrysikopoulos 《Stochastic Environmental Research and Risk Assessment (SERRA)》2001,15(1):33-46
A two-dimensional numerical transport model is developed to determine the effect of aquifer anisotropy and heterogeneity
on mass transfer from a dense nonaqueous phase liquid (DNAPL) pool. The appropriate steady state groundwater flow equation
is solved implicitly whereas the equation describing the transport of a sorbing contaminant in a confined aquifer is solved
by the alternating direction implicit method. Statistical anisotropy in the aquifer is introduced by two-dimensional, random
log-normal hydraulic conductivity field realizations with different directional correlation lengths. Model simulations indicate
that DNAPL pool dissolution is enhanced by increasing the mean log-transformed hydraulic conductivity, groundwater flow velocity,
and/or anisotropy ratio. The variance of the log-transformed hydraulic conductivity distribution is shown to be inversely
proportional to the average mass transfer coefficient. 相似文献
6.
本文提出了一种计算注聚合物井井下温度分布的方法,该方法视井筒内的聚合物溶液为非牛顿幂律流体,考虑注入流体在井筒中同一截面上的速度变化,根据聚合物溶液在多孔介质中的流变性,依据广义达西定律求取渗流速度,基于能量守恒方程建立柱坐标系下注聚合物井井筒内流体、注入层及围岩的二维温度场模型.通过合理的边界条件,将三部分模型耦合起来,采用交替方向半隐式有限差分法求解建立的井下温度场模型.考查了注入量、注入时间、幂律指数和稠度系数、注入液温度等因素对井下温度场分布的影响,结果表明,当注入的聚合物溶液温度低于注入层的原始温度时,随注入量和注入时间的增大以及粘度的减小,注入层的温度降低;注入聚合物溶液温度与注入层原始温度差越大,注入层处的温度剖面异常越明显.本文数值计算结果可用于指导注聚合物井的井温测井应用. 相似文献
7.
《Advances in water resources》2005,28(8):793-806
In this paper, the numerical errors associated with the finite difference solutions of two-dimensional advection–dispersion equation with linear sorption are obtained from a Taylor analysis and are removed from numerical solution. The error expressions are based on a general form of the corresponding difference equation. The variation of these numerical truncation errors is presented as a function of Peclet and Courant numbers in X and Y direction, a Sink/Source dimensionless number and new form of Peclet and Courant numbers in X–Y plane. It is shown that the Crank–Nicolson method is the most accurate scheme based on the truncation error analysis. The effects of these truncation errors on the numerical solution of a two-dimensional advection–dispersion equation with a first-order reaction or degradation are demonstrated by comparison with an analytical solution for predicting contaminant plume distribution in uniform flow field. Considering computational efficiency, an alternating direction implicit method is used for the numerical solution of governing equation. The results show that removing these errors improves numerical result and reduces differences between numerical and analytical solution. 相似文献
8.
The turbulent advection-diffusion mathematical model in three-dimensional space is solved by a mixed finite element finite difference method. Linear finite elements in the vertical direction and central finite differences in the horizontal directions are used coupled with the Galerkin error minimization procedure. The integration in time is performed in fractional steps (one explicit one implicit) by splitting the differential operator. The method is illustrated by application to the three-dimensional movement of suspended sediment. Its accuracy is checked by comparison to analytical solutions and its efficiency is gauged relative to finite elements and implicit finite difference solutions for two-dimensional suspended sediment transport over a dredged channel. 相似文献
9.
Newm ark-更新精细积分法是动力方程求解的隐式的时域逐步积分法,其稳定性条件非常容易满足。与隐式方法相比较,显式积分方法不需要求解耦联的方程组,可以有效地减少内存占用和机时耗费。因此,根据显式积分方法的特点和优点,基于Newm ark-更新精细积分法的基本思想,提出其显式积分格式。对显式积分方法的精度与稳定性进行了初步的分析,指出该显式积分方法具有极好的稳定性,其精度比隐式积分方法的精度稍低。随着时间步长的增加,其精度优于传统的方法。 相似文献
10.
This paper presents the results of the II Tidal Flow Forum experiment of the English Channel and the southern North Sea. The model applies the FADI (falsified alternating direction implicit) scheme and uses the data base of the Tidal Flow Forum. Discrepancies between the model results and the distributed field data of 11 tide stations and 8 tidal current measurement points are shown graphically and are quantified by calculating the RMS (root mean square) errors and the standard deviations. Sensitivity tests have been carried out by changing some parameters (frictions, Coriolis, …). The results which best fit the reference data were obtained by using the Manning's friction law. By doing so, the model can more appropriately adapt the complex bathymetry. The improvements are shown graphically. 相似文献
11.
Comparison between TVD-MacCormack and ADI-type solvers of the shallow water equations 总被引：6，自引：1，他引：6
A total variation diminishing (TVD) modification of the MacCormack scheme is developed for simulating shallow water dynamics on a uniform Cartesian grid. Results obtained using conventional and deviatoric forms of the conservative non-linear shallow water equations (SWEs) are compared for cases where the bed has a varying topography. The comparisons demonstrate that the deviatoric form of the SWEs gives more accurate results than the conventional form, in the absence of numerical balancing of the flux-gradient and source terms. A further comparison is undertaken between the TVD-MacCormack model and an alternating direction implicit (ADI) model for cases involving steep-fronted shallow flows. It is demonstrated that the ADI model is unable to predict trans-critical flows correctly, and artificial viscosity has to be introduced to remove spurious oscillations. The TVD-MacCormack model reproduces all flow regimes accurately. Finally, the TVD-MacCormack model is used to predict a laboratory-scale dyke break undertaken at Delft University of Technology. The predictions agree closely with the experimental data, and are in excellent agreement with results from an alternative Godunov-type model. 相似文献
12.
James Rickett 《Geophysical Prospecting》2013,61(4):874-881
This paper compares three alternative algorithms for simultaneously estimating a source wavelet at the same time as an earth model in full‐waveform inversion: (i) simultaneous descent, (ii) alternating descent and (iii) descent with the variable projection method. The latter is a technique for solving separable least‐squares problems that is well‐known in the applied mathematics literature. When applied to full‐waveform inversion, it involves making the source wavelet an implicit function of the earth model via a least‐squares filter‐estimation process. Since the source wavelet becomes purely a function of medium parameters, it no longer needs to be treated as a separate unknown in the inversion. Essentially, the predicted data are projected onto the measured data in a least‐squares sense at every function evaluation, making use of the fact that the filter estimation problem is trivial when compared to the full‐waveform inversion problem. Numerical tests on a simple 1D model indicate that the variable projection method gives the best result; actually producing results in quality that are very similar to control experiments with a known, correct wavelet. 相似文献
13.
在实际工程结构动力反应分析中,往往由于结构型式十分复杂,常用的两种直接积分方法,即显式积分方法和隐式积分方法,在使用中都存在着一定的局限性,如何将这两种积分方法合理有效地结合起来,是一个十分有意义的研究课题。针对实际工程问题中整体结构计算时间步长的选择往往受局部区域的材料特性、尺寸大小等因素影响的这一现象,提出了一种对结构局部区域进行隐式积分、对其余区域进行显式积分的显隐式积分方法,这种积分格式相对于显式积分格式而言,能显著提高整体结构的计算速度。最后采用两个数值计算实例对这一方法进行验证。 相似文献
14.
We present a vertex-centered finite volume method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media. Fractures are discretely modeled as lower dimensional elements. The method works on unstructured, locally refined grids and on parallel computers with distributed memory. An implicit time discretization is employed and the nonlinear systems of equations are solved with a parallel Newton-multigrid method. Results from two-dimensional and three-dimensional simulations are presented. 相似文献
15.
应力调制与预报研究 总被引：7，自引：0，他引：7
在地球动力作用和地壳运动的过程中构造块体呈现出整体运动的特征,表现在大陆地震的动力学特征具有明显的分区性。因而,大陆每个块体的应力场特征方向与强地震发生的主要构造力方向一致性较好。在它的作用下,块体边界或某些部位出现应力不均匀并有应变能的积累。孕震区应变能由线性积累向临界的非线性转化迹象会表现在弱震活动与主构造方向上固体潮交变应力场的呼应上。 交变应力场可加强孕震区应力集中部位的介质疲劳强度,从而在这些部位触发一些弱震,而且在附加应力场分别为正、负构造应力场方向时发生的弱震频次比例失调,这一过程与构造应力场的呼应显示了强震前弱震震源机制一致性较好的前兆特征,这种情形下的地震活动反映了一种应力调制作用,分析了1970年来的华北资料和1980年来的云南省资料,用直接图像分布、数值和等值线等3种方式讨论,表明它对强震有1~3年的中长期预报意义,R值评估均大于0.5。 相似文献
16.
IINTRODUCTIONEstUariesareprominentcoastalfeatUres.Estuariesareofgreateconomicssignificancetomankind.Attheseareas,manyharborsandwaterchannelshavetobebuiltforeconomicpurposes.ThedesignandconstrUctionofcoastalstrUctUresinestUariesrequireknowledgeofhydrodynamicsaswellassedimenttransportinsuchregions.ThenatUreofestuariesiscontrolledbyvariouscoastalhydrodynamicprocesses.Undertheactionofhydrodynamics,sedimentdepositionsorerosionswilloccurinestuariesornearcoastalstrUCtures.Tomaintainnavigati… 相似文献
17.
This paper presents a new method, called the equivalent force control method, for solving the nonlinear equations of motion in a real‐time substructure test using an implicit time integration algorithm. The method replaces the numerical iteration in implicit integration with a force‐feedback control loop, while displacement control is retained to control the motion of an actuator. The method is formulated in such a way that it represents a unified approach that also encompasses the effective force test method. The accuracy and effectiveness of the method have been demonstrated with numerical simulations of real‐time substructure tests with physical substructures represented by spring and damper elements, respectively. The method has also been validated with actual tests in which a Magnetorheological damper was used as the physical substructure. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
18.
William A. Robertson 《Pure and Applied Geophysics》1964,59(1):93-99
Summary Directions of natural remanent magnetization aftermagnetic cleaning of specimens from monzonite porphyry at Milton are given. Thermal and alternating magnetic field stability tests indicate that the resultant direction is that of the earth's magnetic field at the time of cooling. Comprison of the pole position calculated from this result with other pole positions from rocks of known age confirms the probable Permian age of the intrusion and the wide divergence of Australian from European and North American pole positions of this age. 相似文献
19.
在计算机断层成像（CT）中,伪影会降低重建图像的质量。针对该问题,本文提出一种去除CT图像中同心椭圆伪影的方法。该方法基于方向全变分（DTV）的思想,将椭圆伪影的去除问题建模为能量最小化问题,并通过对椭圆伪影边缘特征的分析,建立适应于椭圆伪影的变分模型。由于提出的模型为具有可分裂结构的非光滑凸优化问题,因此本文使用交替方向乘子法（ADMM）。最后,通过仿真实验验证该模型对于去除椭圆伪影的有效性。 相似文献
20.
Multiphase flow in porous media is described by coupled nonlinear mass conservation laws. For immiscible Darcy flow of multiple fluid phases, whereby capillary effects are negligible, the transport equations in the presence of viscous and buoyancy forces are highly nonlinear and hyperbolic. Numerical simulation of multiphase flow processes in heterogeneous formations requires the development of discretization and solution schemes that are able to handle the complex nonlinear dynamics, especially of the saturation evolution, in a reliable and computationally efficient manner. In reservoir simulation practice, single-point upwinding of the flux across an interface between two control volumes (cells) is performed for each fluid phase, whereby the upstream direction is based on the gradient of the phase-potential (pressure plus gravity head). This upwinding scheme, which we refer to as Phase-Potential Upwinding (PPU), is combined with implicit (backward-Euler) time discretization to obtain a Fully Implicit Method (FIM). Even though FIM suffers from numerical dispersion effects, it is widely used in practice. This is because of its unconditional stability and because it yields conservative, monotone numerical solutions. However, FIM is not unconditionally convergent. The convergence difficulties are particularly pronounced when the different immiscible fluid phases switch between co-current and counter-current states as a function of time, or (Newton) iteration. Whether the multiphase flow across an interface (between two control-volumes) is co-current, or counter-current, depends on the local balance between the viscous and buoyancy forces, and how the balance evolves in time. The sensitivity of PPU to small changes in the (local) pressure distribution exacerbates the problem. The common strategy to deal with these difficulties is to cut the timestep and try again. Here, we propose a Hybrid-Upwinding (HU) scheme for the phase fluxes, then HU is combined with implicit time discretization to yield a fully implicit method. In the HU scheme, the phase flux is divided into two parts based on the driving force. The viscous-driven and buoyancy-driven phase fluxes are upwinded differently. Specifically, the viscous flux, which is always co-current, is upwinded based on the direction of the total-velocity. The buoyancy-driven flux across an interface is always counter-current and is upwinded such that the heavier fluid goes downward and the lighter fluid goes upward. We analyze the properties of the Implicit Hybrid Upwinding (IHU) scheme. It is shown that IHU is locally conservative and produces monotone, physically-consistent numerical solutions. The IHU solutions show numerical diffusion levels that are slightly higher than those for standard FIM (i.e., implicit PPU). The primary advantage of the IHU scheme is that the numerical overall-flux of a fluid phase remains continuous and differentiable as the flow regime changes between co-current and counter-current conditions. This is in contrast to the standard phase-potential upwinding scheme, in which the overall fractional-flow (flux) function is non-differentiable across the boundary between co-current and counter-current flows. 相似文献