稠油油藏提高采收率技术数值模拟研究

渤海地区具有十分丰富的稠油资源,其油藏储层胶结疏松、渗透率高、非均质性严重并且边底水活跃,水驱采收率较低.因此,必须采取调剖、堵水和降黏组合技术措施来进行稠油开发,与其它稠油开发措施相比较,稠油乳化降黏技术具有操作简便和经济性好等优点,其中强化分散体系由于其良好的乳化降黏性以及简易的注入工艺受到石油科技工作者的高度重视.为使稠油乳化降黏技术尽快投入矿场试验,针对渤海普通稠油油藏开发实际需求,以油藏工程理论为指导,以数值模拟方法为技术手段,以正交试验原理安排实验,以A油田B1井组为试验平台,对比了不同开发方式对目标井组开发效果的影响.结果表明,单独注入强化分散体系4个月以上增油效果要好于单独调剖,而调剖与注强化分散体系联合作用效果要明显高于单独调剖与单独注强化分散体系之和,达到了"1+1"大于2的功效;高低浓度注入优于单一浓度段塞;水井调剖和油井堵水会明显增加强化分散体系增油有效期和增油量,推荐水井采取"调剖+强化分散体系高低浓度交替注入"组合措施,油井采取"堵水+强化分散体系吞吐"组合措施。     

煤炭资源最佳开采时期研究

煤炭资源是一种不可再生资源,即耗竭性资源.现在的煤炭资源开采策略是一种掠夺式的开采,不仅浪费资源,而且会使我国的煤炭将来很难满足现代化建设的需要.为此,本文采用经济控制论模型试图探讨煤炭资源的最佳开采时期.     

连续铸钢保护渣热传递三维数值模拟

本文讨论了连续铸钢保护渣层热过程的物理 ,数学模型 ,采用三维有限元方法进行了数值计算 ,较成功地模拟了保护渣的热传递过程     

突然扩张方管中三维湍流流动的数值模拟

本文运用SIMPLEC算法计算了突然扩张方管中的三维湍流流动,湍流模型采用k-ε模型。计算结果详细反映了突然扩张方管中三维湍流流场。从本文结果可以看出,由于突然扩张方管几何形状非轴对称,且尺寸有限,边壁对流场的作用是不可忽略的。以往文献中常见的二维突然扩张湍流的数值模拟结果与三维情况有较大差别,在靠近边壁的区域差别很大,因此对于突然扩张方管中湍流流动的数值模拟应用三维模拟。本文计算所得突然扩张截面后主回流区长度与实验结果接近。本文方法可为数值模拟突然扩张方管中湍流流场及各物理参数的分布提供有效工具。     

基于正交设计试验分析法的稠油油藏边水侵入主控因素研究

含强边水的稠油油藏,在蒸汽吞吐开发后期往往出现边水突进现象,致使油井过早水淹,油田高效开发难以维持,因此准确的分析出影响边水侵入的主控因素是进行边水有效治理必须解决的关键问题.结合河南油田BQ67区块的实际油藏参数建立模拟含边水稠油油藏蒸汽吞吐模型,模拟正交设计试验各个方案,通过对结果进行极差、方差分析得出各因素对边水...     

纵向非均质性对海上稠油油田水驱油开发效果的影响研究

为兼顾生产成本与开发效果,海上稠油油田开发大多采用多层合采,层内和层间的差异对水驱开发效果影响显著,针对纵向非均质严重的问题,基于不同渗透率级差、不同平均渗透率,建立了反韵律下9种非均质物理模型,研究层内和层间非均质性对海上稠油油田水驱开发效果的影响.研究结果表明:对于层内非均质模型,渗透率级差对开发效果影响更大,平均渗透率相近时,级差越大,层内非均质性越强,注入水沿高渗透层突进,油层下部驱替效果越差,采出程度越低;层内均质层间非均质时,级差越大采出程度越低,高渗层段注入水形成优势通道抑制低渗透层吸水,高渗透层最高采收率达77.4%,低渗层采收率可低至0.7%,使油层上下部驱替效果相差较大;层内和层间都非均质时,分层开采的开发效果好于多层合采,各层采收率变化幅度与分层开采时差异较大;利用综合影响因子表征层内和层间非均质性的综合影响,其值为0.08时对应采出程度为41.48%,综合影响因子越高,纵向非均质性越强,采出程度越低.     

成本约束下不可再生资源开采微分对策问题研究

研究开采成本等重要因素对不可再生资源开采的影响.假定有N个厂商对不可再生资源进行开采,其面临的市场需求价格弹性为常数,开采成本是开采量的线性函数.利用微分对策理论,分析该寡头市场中各个厂商的均衡开采策略.研究表明,开采成本对是否开采、开采速度、结束时点等有重要影响.同时,本文还与某些经典论文进行了简要比较.     

基于浸入边界-有限元法的流固耦合碰撞数值模拟方法

针对流固耦合碰撞问题,建立了流体中固体与固体碰撞界面解析直接模拟方法,采用清晰界面浸入边界法模拟流体中的动边界问题,避免了传统贴体网格方法在求解流体中存在固体间碰撞问题时网格出现负体积的问题,采用基于罚函数的有限元方法对固体的运动和碰撞进行求解,以分域耦合方式实现流体域和固体域的耦合求解.通过与静止流体中球形颗粒与壁面正碰撞和斜碰撞的实验数据对比,验证了建立的数值模拟方法对流体中固体与固体碰撞数值模拟的正确性,获得了流体域流场在碰撞前后随时间的变化,同时通过该文建立的数值模拟方法也获得了固体域中固体的碰撞力和应力.未来,将把该数值模拟方法应用到流体流动环境中,如固体颗粒对管道的冲蚀、流体诱导海洋立管之间的碰撞、坠物对海底管道的撞击等.     

直接数值模拟／大涡模拟中数值误差影响的研究

通过比较湍流的能谱和总动能,对数值误差(包括混淆误差、离散截断误差)、亚格子模型以及它们之间相互作用对直接数值模拟和大涡模拟的影响进行了系统研究.算例采用了三维各向均匀同性湍流.为了研究复杂几何形状,数值格式采用了谱方法和Pade紧致格式.大涡模型采用了truncated Navior-Stokes(TNS)模型结合Pade离散滤波器.结果表明直接数值模拟中离散误差对结果有很大影响,低阶格式会导致计算发散.而大涡模拟中亚格子模型不仅能表征小尺度对大尺度的影响,而且还缓解了数值误差对计算结果的影响.因而低精度格式也可取得不错的结果.     

三维数值流形方法的理论研究

在二维数值流形方法的基础上,对三维数值流形进行了理论研究．研究了三维覆盖位移函数,进行了三维数值流形的力学分析,给出了三维流形单元的刚度矩阵,详细推导了三维数值流形的Hammer积分及剖分规则,系统地研究三维数值流形的理论体系与数值实现方法．作为数值算例,给出了相应的悬臂梁的计算结果,计算结果表明算法的精度和计算效益较高．     

Temporally regularized direct numerical simulation

Experience with fluid-flow simulation suggests that, in some instances, under-resolved direct numerical simulation (DNS), without a residual-stress model per se but with artificial damping of small scales to account for energy lost in the cascade from resolved to unresolved scales, may be as reliable as simulations based on more complex models of turbulence. One efficient and versatile manner to selectively damp under-resolved spatial scales is by a relaxation regularization, e.g. Stolz and Adams [S. Stolz, N.A. Adams, An approximate deconvolution procedure for large eddy simulation, Phys. Fluids II (1999) 1699-1701]. We consider the analogous approach based on time scales, time filtering and damping of under-resolved temporal features. The paper explores theoretical and practical aspects of temporally damped fluid-flow simulations. We prove existence of solutions to the resulting continuum model. We also establish the effect of the damping of under-resolved temporal features as the energy balance and dissipation and prove that the time fluctuations → 0 in a precise sense. The method is then demonstrated to obtain both steady-state and time-dependent coarse-grid solutions of the Navier-Stokes equations.     

On numerical simulation of tsunami run-up

The improved adaptive grid method is presented for the numerical modelling of propagation of long surface waves in large water areas and their interaction with coasts. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Block‐centered upwind multistep difference method and convergence analysis for numerical simulation of oil reservoir

The three‐dimensional displacement of two‐phase flow in porous media is a preliminary problem of numerical simulation of energy science and mathematics. The mathematical model is formulated by a nonlinear system of partial differential equations to describe incompressible miscible case. The pressure is defined by an elliptic equation, and the concentration is defined by a convection‐dominated diffusion equation. The pressure generates Darcy velocity and controls the dynamic change of concentration. We adopt a conservative block‐centered scheme to approximate the pressure and Darcy velocity, and the accuracy of Darcy velocity is improved one order. We use a block‐centered upwind multistep method to solve the concentration, where the time derivative is approximated by multistep method, and the diffusion term and convection term are treated by a block‐centered scheme and an upwind scheme, respectively. The composite algorithm is effective to solve such a convection‐dominated problem, since numerical oscillation and dispersion are avoided and computational accuracy is improved. Block‐centered method is conservative, and the concentration and the adjoint function are computed simultaneously. This physical nature is important in numerical simulation of seepage fluid. Using the convergence theory and techniques of priori estimates, we derive optimal estimate error. Numerical experiments and data show the support and consistency of theoretical result. The argument in the present paper shows a powerful tool to solve the well‐known model problem.     

电离层层析成像及数值模拟

利用正则化反演法,对电离层层析成像技术进行了数值模拟.模拟结果表明,正则化反演法能够准确反演出电离层赤道异常的电子密度分布.对正则化反演误差的原因进行了分析,定义了代表重建结果精度的参数并研究其随信号数量的变化,结果显示造成反演误差的主要原因是信息量不足,在提供充足的数据条件下,正则化反演法能够得到很高精度的重建结果.数值模拟实验进一步证明了电离层层析技术在反演电离层电子密度分布中有很大的应用价值.     

Direct numerical simulation of compressible isotropic turbulence

Direct numerical simulation (DNS) of decaying compressible isotropic turbulence at turbulence Mach numbers of M_t = 0.2-0.7 and Taylor Reynolds numbers of 72 and 153 is performed by using the 7th order upwind-biased difference and 8th order center difference schemes. Results show that proper upwind-biased difference schemes can release the limit of“start-up” problem to Mach numbers. Compressibility effects on the statistics of turbulent flow as well as the mechanics of shocklets in compressible turbulence are also studied, and the conclusion is drawn that high Mach number leads to more dissipation. Scaling laws in compressible turbulence are also analyzed. Evidence is obtained that scaling laws and extended self similarity (ESS) hold in the compressible turbulent flow in spite of the presence of shocklets, and compressibility has little effect on scaling exponents.     

Turbulent diffusion: Direct numerical simulation

Direct numerical simulations (DNS) of a vibrating grid turbulence at two different Reynolds numbers (based on an amplitude and frequency of the grid) Re = 500, 1000 are presented. The evolution of a turbulent/non-turbulent interface was detected and compared to the theoretical results obtained in the paper [1] using Lie group analysis of the governing equations. Good agreement was found between DNS and theory. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Nonlinear numerical simulation of ferroelectric-ferromagnetic multifunctional composites

The theoretical background of nonlinear constitutive multifield behavior is presented. Nonlinear material models describing the ferroelectric or ferromagnetic behaviors are presented. Both physically and phenomenologically motivated constitutive models have been developed for the numerical calculation of the nonlinear magnetostrictive and ferroelectric behaviors. On this basis, the polarization in the ferroelectric and magnetization in the ferromagnetic respectively magnetostrictive phases are simulated and the resulting effects analyzed. The developed tools enable the prediction of the electromagnetomechanical properties of smart multiferroic composites and supply useful means for their optimization. Goals are to improve the efficiency of ME coupling and to reduce damage associated with the poling process. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Theoretical study and numerical simulation of textiles

We propose a new approach for developing continuum models fit to describe the mechanical behavior of textiles. We develop a physically motivated model, based on the properties of the yarns, which can predict and simulate the textile behavior. The approach relies on the selection of a suitable topological model for the patch of the textile, coupled with constitutive models for the yarn behavior. The textile structural configuration is related to the deformation through an energy functional, which depends on both the macroscopic deformation and the distribution of internal nodes. We determine the equilibrium positions of these latter, constrained to an assigned macroscopic deformation. As a result, we derive a macroscopic strain energy function, which reflects the possibly nonlinear character of the yarns as well as the anisotropy induced by the microscopic topological pattern. By means of both analytical estimates and numerical experiments, we show that our model is well suited for both academic test cases and real industrial textiles, with particular emphasis on the tricot textile.     

Mathematical simulation of oil reservoir properties

The study and computational representation of porous media properties are very important for many industries where problems of fluid flow, percolation phenomena and liquid movement and stagnation are involved, for example, in building constructions, ore processing, chemical industries, mining, corrosion sciences, etc. Nevertheless, these kinds of processes present a noneasy behavior to be predicted and mathematical models must include statistical analysis, fractal and/or stochastic procedures to do it. This work shows the characterization of sandstone berea core samples which can be found as a porous media (PM) in natural oil reservoirs, rock formations, etc. and the development of a mathematical algorithm for simulating the anisotropic characteristics of a PM based on a stochastic distribution of some of their most important properties like porosity, permeability, pressure and saturation. Finally a stochastic process is used again to simulated the topography of an oil reservoir.