NUMERICAL SIMULATION OF OIL MIGRATION-ACCUMULATION OF MULTILAYER AND ITS APPLICATION

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｏｎｅｏｆｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ.Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃ…     

Parallel arithmetic numerical simulation and application of secondary migration-accumulation of oil resources

From such actual conditions as the effects of characteristics of miltilayer petroleum geology and permeation fluid mechanics, a new numerical model is put forward and coupling splitting-up implicit interactive scheme is formed. For the actual situation of Dongying hollow ( four-layer ) and Tanhai region ( three-layer ) of Shengli Petroleum Field, the numerical simulation test results and the actual conditions are coincident.     

Numerical simulation analysis for migration-accumulation of oil and water

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｅｖｏｌｕｔｉｏｎａｒｙｈｉｓｔｏｒｙｏｆｔｈｅｂａｓｉｎｉｓｃｏｎｄｕｃｔｅｄａｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｅｔｒｏｌｅｕｍｇｅｏｌｏｇｙａｎｄｍｅｃ...     

NUMERICAL METHOD AND APPLICATION FOR THREEDIMENSIONAL NONLINEAR SYSTEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

For the system of multilayer dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources.     

A model for simulation of friction phenomenon between dies and workpiece

Based on the interaction of asperities and upperbound approach a mathematical model for simulation of friction phenomenon between dies and workpiece is proposed. Optimizing the mathematical model with respect to several variables it is found that in addition to adhering, tearing, ploughing, etc., asperities workpiece can move wave-like along the surface layer and under certain circumstances they may disappear. If the asperities wavily move along the surface layer the friction coefficient depends on the geometry of asperities. However, the bonding strength of asperities has no significant influence on friction coefficient. The depth of the plastic deformation layer is related to the geometry of asperities, too. The soundness of the prerequisite of the proposed model and some analytical results were verified by experiments.Department of Materials Engineering, Northwestern Polytechnical UniversityThis is a project of the cooperative program between the Department of Materials Engineering of Northwestern Polytechnical University of the People's Republic of China and the Institute for Metal Forming at Lehigh University of the United States of America.     

Krylov subspace projection method and its application on oil reservoir simulation

ＩｎｔｒｏｄｕｃｔｉｏｎＭａｎｙｐｒｏｂｌｅｍｓｏｆｅｎｇｉｎｅｅｒｉｎｇｉｎｔｅｒｅｓｔ,ｉｎｃｌｕｄｉｎｇｐｒｏｂｌｅｍｓｏｆｏｉｌｒｅｓｅｒｖｏｉｒｓｉｍｕｌａｔｉｏｎ,ｌｅａｄｔｏｔｈｅｔａｓｋｏｆｓｏｌｖｉｎｇａｓｙｓｔｅｍｏｆｌｉｎｅａｒｅｑｕａｔｉｏｎｓＡｘ=ｆ,(1)ｗｈｅｒｅＡｉｓａｎｎ×ｎ,ｒｅａｌ,ｌａｒｇｅｓｐａｒｓｅ,ｎｏｎ＿ＳＰＤｍａｔｒｉｘ.Ｔｈｅｄｉｒｅｃｔｄｅｃｏｍｐｏｓｉｔｉｏｎｍｅｔｈｏｄｉｓｍｏｓｔｅｆｆｉｃｉｅｎｔｓｏｌｖｅｒｗｈｅｎｎｉｓｓｍａｌｌ.Ｈｏｗｅｖｅｒ…     

Careful numerical simulation and analysis of migration-accumulation of Tanhai Region

Numerical simulation of careful parallel arithmetic of oil resources migration-accumulation of Tanhai Region (three-layer) was done. Careful parallel operator splitting-up implicit iterative scheme, parallel arithmetic program, parallel arithmetic information and alternating-direction mesh subdivision were put forward.Parallel arithmetic and analysis of different CPU combinations were done. This numerical simulation test and the actual conditions are basically coincident. The convergence estimation of the model problem has successfully solved the difficult problem in the fields of permeation fluid mechanics, computational mathematics and petroleum geology.     

杆式射流形成的数值模拟研究

利用数值仿真软件LS-DYNA3D,通过数值模拟方法研究了药型罩压跨、射流的形成、延伸和失稳断裂的全过程,数值模拟结果与实验结果吻合较好,得到了一套实用的数值模拟方法及材料模型参数,定性分析了大锥角的药型罩压垮过程,对形成杆式射流的药型罩结构设计具有一定借鉴意义。     

Numerical Simulation of a Turbulent Thermal

Ascent of a large-scale thermal in a standard atmosphere is calculated with the use of the Reynolds equations and the k model of turbulence, which takes into account temperature inhomogeneity and vorticity of the flow, and the Euler equations. Results of numerical calculations of a flow examined experimentally are presented. Gas-dynamic and turbulent flow parameters obtained in calculations and experiments are compared.     

三维爆炸与冲击问题的可视化研究

结合爆炸与冲击物理现象的特点,对三维爆炸场的可视化进行了研究,在三维多流体网格法基础上自行研制了三维爆炸场的可视化工具软件ViSC3D。利用ViSC3D,对两种复杂边界的爆炸场进行了可视化,详细地展现了爆炸场的概貌和局部细节,表明该软件效果较好。     

弹塑性变形油藏中多相渗流的数值模拟

基于流固耦合力学理论,建立了弹性变形油藏中多相渗流的数学模型,该模型考虑了渗流与变形的耦合作用,以及注采交变载荷作下油藏多孔介质的弹性变形特征,给出了耦合数值模拟方法和算例。     

Simulation and Modelling of Turbulent Trailing-Edge Flow

Computations of turbulent trailing-edge flow have been carried out at a Reynolds number of 1000 (based on the free-stream quantities and the trailing-edge thickness) using an unsteady 3D Reynolds-Averaged Navier–Stokes (URANS) code, in which two-equation (k–ε) turbulence models with various low-Re near wall treatments were implemented. Results from a direct numerical simulation (DNS) of the same flow are available for comparison and assessment of the turbulence models used in the URANS code. Two-dimensional URANS calculations are carried out with turbulence mean properties from the DNS used at the inlet; the inflow boundary-layer thickness is 6.42 times the trailing-edge thickness, close to typical turbine blade flow applications. Many of the key flow features observed in DNS are also predicted by the modelling; the flow oscillates in a similar way to that found in bluff-body flow with a von Kármán vortex street produced downstream. The recirculation bubble predicted by unsteady RANS has a similar shape to DNS, but with a length only half that of the DNS. It is found that the unsteadiness plays an important role in the near wake, comparable to the modelled turbulence, but that far downstream the modelled turbulence dominates. A spectral analysis applied to the force coefficient in the wall normal direction shows that a Strouhal number based on the trailing-edge thickness is 0.23, approximately twice that observed in DNS. To assess the modelling approximations, an a priori analysis has been applied using DNS data for the key individual terms in the turbulence model equations. A possible refinement to account for pressure transport is discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

爆炸荷载下岩石破坏的数值流形方法模拟

为了更好地利用数值流形方法对动力学问题进行分析,在对原数值流形方法中的动力学问题求解思想进行分析的基础上,采用动力有限元方法中的Newmark法对该算法进行了改进。改进后的数值流形方法与原来相比具有三个明显的优势：（1）当选择合适的参数后,该方法能够保证解的无条件收敛;（2）可以采用比原算法大得多的时间步长;（3）充分考虑了动力学问题中的阻尼效应。最后通过一个算例说明了改进后的数值流形方法能够很好地模拟岩石在冲击载荷作用下破坏的全过程,克服了有限元法不能模拟岩石破坏后块体运动情况的不足。     

方柱绕流的数值模拟

采用有限差分法,对雷诺数为2.2×10~4的方柱绕流进行了大涡模拟(简称LES)。运用时间分裂控制(Split-Operator)法,将N-S方程分为对流步、扩散步和传播步。对Smagorinsky假设在近壁区的发散问题用两层模型进行处理。对流项用迎风—中心差分格式模拟,压力方程用SOR法迭代求解。计算得到的沿对称线的时均顺流向速度与文献上的实验结果进行了比较,结果吻合较好,同时还对绕方柱流的流场结构进行了分析研究。     

Simulation of multiphase flows in porous media

The ability to numerically simulate single phase and multiphase flow of fluids in porous media is extremely important in developing an understanding of the complex phenomena governing the flow. The flow is complicated by the presence of heterogeneities in the reservoir and by phenomena such as diffusion, dispersion, and viscous fingering. These effects must be modeled by terms in coupled systems of nonlinear partial differential equations which form the basis of the simulator. The simulator must be able to handle both single and multiphase flows and the transition regimes between the two. A discussion of some of the aspects of modeling dispersion and viscous fingering is presented along with directions for future work.The partial differential equation models are convection-dominated and contain important local effects. An operator-splitting technique is used to address these different effects accurately. Convection is treated by time stepping along the characteristics of the associated pure convection problem, and diffusion is modeled via a Galerkin method for single phase flow and a Petrov-Galerkin technique for multiphase regimes. ELLAM (Eulerian-Lagrangian Localized Adjoint Methods) are discussed to effectively treat the advection-dominated processes. Accurate approximations of the fluid velocities needed in the Eulerian-Lagrangian time-stepping procedure are obtained by mixed finite element methods. Adaptive local grid refinement techniques are then indicated to resolve important local phenomena around wells and large heterogeneities or to resolve the moving internal boundary layers which often govern the mass transfer between phases.     

液固两相流冲洗油管道的冲蚀磨损特性数值模拟及分析

以液固两相流冲洗油管道为研究对象,采用Realizableκ-ε湍流模型、随机轨道模型,结合液固两相流冲蚀磨损试验,建立修正的冲蚀磨损数理模型,数值预测典型工况下冲洗油管道内速度、压力、冲蚀磨损率等流动参数分布情况,分析了冲蚀磨损的形成机制.研究结果表明:受曲率半径影响,冲洗油管道冲蚀磨损速率随曲率半径的增加而减小;由于颗粒惯性及管内二次流影响,弯头中间区域外侧壁面和出口直管段内侧面磨损较为严重,三通管件的最大冲蚀磨损率位于孔口处,数值预测结果与失效解剖案例吻合.本文建立的冲蚀磨损定量预测方法,适用于压力管道的风险评定及寿命预测.     

Simulation and transport phenomena of a ternary two-phase flow

A chemical flood model for a three-component (petroleum, water, injected chemical) two-phase (aqueous, oleic) system is presented. It is ruled by a system of nonlinear partial differential equations: the continuity equation for the transport of each of its components and Darcy's equation for the two-phase flow. The transport mechanisms considered are ultralow interfacial tension, capillary pressure, dispersion, adsorption, and partition of the components between the fluid phases (including solubilization and swelling).The mathematical model is numerically solved in the one-dimensional case by finite differences using an explicit and direct iterative procedure for the discretization of the conservation equations. Numerical results are compared with Yortsos and Fokas' exact solution for the linear waterflood case including capillary pressure effects and with Larson's model for surfactant flooding. The effects of the above-mentioned transport mechanisms on concentration profiles and on oil recovery are also analyzed.     

Effect of flows between formations and capillary forces on the process of oil displacement in a stratified inhomogeneous bed

The results of a numerical investigation of the process of oil displacement in a stratified inhomogeneous formation on the basis of the two-phase flow model with account for capillary forces are presented. It is shown that in many cases the vertical inhomogeneity of oil reservoirs may not be a cause of nonuniform displacement and the non-recovery of large oil reserves by the time of water breakthrough to the extraction surface. The action of the capillary forces is an additional factor leading to equalization of the water propagation front in the inhomogeneous formation, water breakthrough delay, and intensification of the mass transfer between the layers with different permeabilities. Analysis of the contribution of the interlayer flows to the water flooding of low-permeability formation intervals calls into question the practicability of blocking high-permeability inclusions in the neighborhood of pumping wells.     

Simulation of micellar-polymer flooding of a layered oil reservoir of nonuniform thickness

The dynamics of oil displacement from a layered reservoir of nonuniform thickness consisting of two hydrodynamically connected layers of different absolute permeability is studied. Results of numerical calculations are given. The influence of the main determining factors on the oil displacement dynamics is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 6, pp. 115–123, November–December, 2008.     

超高速撞击数值仿真结果分析

为了分析超高速撞击过程的宏观现象和内在机理,对9.53 mm铝球以6.64 km/s的速度撞击2.2 mm厚的铝靶的SPH仿真结果进行了量化分析。结果表明：采用SPH方法以及Steinberg弹塑性模型和Mie-Grneisen状态方程,可获得与试验相符的仿真结果;球形破片开坑或穿孔直径遵循初始快速增加、然后缓慢增加,直至稳定的变化规律;破片/靶板界面的最大撞击压力比材料强度大两个量级以上;靶板阻抗力最大值发生在破片最大直径侵入靶板时刻;碎片云的运动过程具有自相似演化特征,其运动范围不会超出碎片云的包络圆锥范围。