An Improvement on Modeling of Forced Gravity Drainage in Dual Porosity Simulations Using a New Matrix-Fracture Transfer Function

In fractured oil reservoirs, the gravity drainage mechanism has great potentials to higher oil recovery in comparison with other mechanisms. Recently, the forced gravity drainage assisted by gas injection has also been considered; however, there are few comprehensive studies in the literature. Dual porosity model, the most common approach for simulation of fractured reservoirs, uses transfer function concept to represent the fluid exchange between matrix and its neighborhood fractures. This study compares the results of different available transfer functions with those of fine grid simulations when forced gravity drainage contributes to oil production from a single matrix block. These comparisons can lead to a more sophisticated formulation including the interplay of capillary, gravity and viscous forces. As a result, a new matrix-fracture transfer function can be developed and its results can be tested against the results of fine-grid simulations. Moreover, the reliability of this model for simulation of forced gravity drainage has been demonstrated by performing some sensitivity analysis.     

The Conley Index for Fast-Slow Systems I. One-Dimensional Slow Variable

We develop a qualitative theory for fast-slow systems with a one-dimensional slow variable. Using Conley index theory for singularity perturbed systems, conditions are given which imply that if one can construct heteroclinic connections and periodic orbits in systems with the derivative of the slow variable set to 0, these orbits persist when the derivative of the slow variable is small and nonzero.     

Dual Porosity and Dual Permeability Modeling of Horizontal Well in Naturally Fractured Reservoir

This article is the first investigation on the dual permeability flow issue for horizontal well-production in a naturally fractured dual-porosity reservoir. Based on the inter-porosity flow from matrix system to fracture system and treating the media directly connected with horizontal wellbore as matrix and fracture systems, we established a model of horizontal well-production and then solved the model using some modern mathematical methods, such as Laplace integral transformation, separation of variables, eigenvalue, and eigenfunction. Later in the article, we obtained the standard log–log type curves using numerical simulation and analyzed the transient flow behavior thoroughly, which showed it is dual porosity and dual permeability flow behavior. The numerical simulation results showed that there are obvious differences between dual permeability and single permeability models. The dual permeability flow behavior accelerates energy supplement during production and reduces the classical matrix-fracture (V-shaped) response. We also showed that type curves characteristics are affected by external boundary conditions, the parameter κ, ω _f and λ _mf, etc. The research results show that our model would be a good semi-analytical model supplied to users. Because the single permeability modeling ignores the direct fluid supply from matrix to wellbore, we recommend using the dual permeability modeling to make well testing and rate decline interpretation in real case studies.     

双重孔隙介质非线性流固耦合渗流

本文给出了考虑双孔双涌介质生变形的流固耦合渗流模型。不仅考虑了固结对渗流的影响,同时也考虎了固体变形对渗流参数（孔隙度和渗透率）的影响。这样。渗流就成了双孔双渗介质中非线性流固耦合渗流。在此基础上,本文还推导了双重孔隙介质非线性流固耦合渗流计算。给出了算例并作了对比。结果表明,固体变形引起的介质参数变化对流体渗流早中期过程有重要的影响,对渗流后期影响并不大。这对于石油开采有重要的参考价值。     

Analysis of the Vibrations of Flexible Discrete-Continuous One-Dimensional Systems with Nonpotential Deformation Characteristics

An approach is proposed to solve dynamic problems for discrete-continuous flexible one-dimensional systems with nonpotential deformation characteristics. The approach is based on a generalized stationarity principle. The solution algorithm employs cubic spline functions. A numerical example demonstrates the capabilities of this approach in determining the loads and displacements in such a system under external influence__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 2, pp. 68–74, February 2005.     

Technical Stability of Nonstationary Automatic-Control Systems with Variable Structure

Sufficient conditions of technical stability in measure are established for nonstationary automatic-control systems with variable structure and logic control laws dependent on the mismatch error and its derivatives of finite order for all admissible initial perturbations from a predefined measurable set of initial perturbations. The associated systems of differential equations contain time-dependent coefficients. The logic control laws are described by a variable jump control function of the mismatch coordinate and its derivatives of finite order that is no higher than that of the initial system of equations. Using a signal and its derivatives for control increases the quality of discontinuous control. The relationship between the eigenvalues of the quadratic forms of the corresponding Liapunov functions and the criteria of technical stability is revealed. The general results are applied to a variable-structure system of the third order     

Diffusion of Tracer in Altered Tonalite: Experiments and Simulations with Heterogeneous Distribution of Porosity

Numerical time-domain-diffusion simulations were used for studying the diffusion behavior of tracer molecules in rock matrix with homogeneous and heterogeneous porosity. As the heterogeneous sample in these simulations, a 3D tomographic image of altered tonalite was used, in which the mineral components and the pores resolved by X-ray microtomography were represented by their respective intragranular porosities determined previously by the ¹⁴C-PMMA method. The apparent diffusion coefficient of a tracer in altered tonalite was determined experimentally, and was then used in the simulations. In the altered tonalite analyzed, inclusion of heterogeneity in the porosity increased the diffusion coefficient by 16 %. Altered and pristine feldspar was the main mineral component in the sample (72 %), and it also provided the dominant contribution to tracer diffusion, explaining alone 52 % of the diffusion coefficient. The large pores resolved by microtomography (6 %) and altered and pristine mica (22 %) gave an equal contribution to the diffusion coefficient. The simulation method applied was also validated by comparing the results to both an analytical and a numerical solution to the diffusion equation in a homogenous medium. In addition, the method was compared to discrete-time random-walk simulations in the case of randomly placed overlapping spheres.     

Stability of the Sliding Mode in Nonstationary Automatic-Control Systems with Variable Structure

The paper addresses the stability of the sliding mode in nonstationary automatic-control systems of variable structure. The parameters of dynamic states are assumed to vary arbitrarily and continuously within a given bounded domain. The existence and stability conditions for the sliding mode on the switching hyperplane are established. They define nonstationary sets that contain phase velocity vectors directed away from the sliding hyperplane. The sets do not have common points with the switching hyperplane, except for the origin, for the known extreme values of the discontinuous coefficients of logical controls. This is the essence of the necessary and sufficient conditions for the stable sliding mode of the automatic-control system. The general results are used to analyze the stability conditions for the sliding mode in a nonstationary automatic-control system with a variable structure of the third order     

Three-Dimensional Simulations for Convection Problem in Anisotropic Porous Media with Nonhomogeneous Porosity,Thermal Diffusivity,and Variable Gravity Effects

A convection problem in anisotropic and inhomogeneous porous media has been analyzed. In particular, the effect of variable permeability, thermal diffusivity, and variable gravity with respect to the vertical direction, has been studied. A linear and nonlinear stability analysis of the conduction solution has been performed. The validity of both the linear instability and global nonlinear energy stability thresholds are tested using a three- dimensional simulation. Our results show that the linear threshold accurately predicts on the onset of instability in the basic steady state. However, the required time to arrive at the basic steady state increases significantly as the Rayleigh number tends to the linear threshold.     

电磁对偶有限元及在波导计算中的应用

力学中的Hamilton体系采用对偶变量描述问题。电磁场采用电场和磁场两类变量描述问题。将力学中的Hamilton体系引入到电磁场问题中,电场变量和磁场变量构成对偶变量,把频域电磁场的基本方程导向对偶方程形式,建立电磁场有限元所需的对偶变量变分原理,由此推导出电磁对偶有限元。将电磁对偶有限元应用于电磁波导计算中,可确定电磁波导的传播常数。文中给出了用电磁对偶有限元方法,计算矩形波导不同模式对应的传播常数的数值计算结果。     

多输入多输出振动传递系统的模型及算法

本文在对多输入多输入振动传递系统进行结构分析的基础上，建立了该系统的网络模型和数学规划模型，给出了一种求振动传递总能量的基本最大流算法，并通过分析一简单实例说明此算法的具体应用。     

A (Dual) Network Model for Heat Transfer in Porous Media

Transport in Porous Media - We present a dual network model to simulate coupled single-phase flow and energy transport in porous media including conditions under which local thermal equilibrium...     

Stability of Nonstationary Automatic-Control Systems of Variable Structure in Forced Motion

Conditions of stability in sliding mode are obtained for nonstationary dynamic automatic-control processes of variable structure. The nonstationary parameters of a control process subjected to external perturbations vary in given ranges, whereas the parameters of the switching plane remain constant for the appropriately chosen parameters of discontinuous logic control laws. The corresponding system of differential equations contains coefficients that vary stepwise in time together with the parameters of the controlled process. The general stability criteria for the solutions of variable-structure systems are used to analyze the stability of motion in sliding mode of an automatic variable-structure control system of the second order     

Technical Stability Conditions for Systems of Variable Structure with Derivatives of the Control Function

We derive sufficient conditions for the technical stability of forced states of automatic control systems with a variable structure depending on the derivatives of the control function. For arbitrary admissible initial perturbations from a measurable set of initial states of a control system with filtering, the technical stability conditions do not depend on the sliding conditions in the given domain of the system parameters. We find how the properties of the eigenvalues of the quadratic forms of the Lyapunov functions are related to the technical stability conditions for automatic control systems of variable structure with filtering. The technical stability conditions are analyzed for automatic control systems with third-order filtering where the forming signal creates eight structures     

Technical Stability of Forced Motion in Nonstationary Automatic-Control Systems of Variable Structure

Sufficient conditions of technical stability are established for nonstationary automatic-control systems of variable structure. The case is considered where for appropriate controls and for all admissible initial perturbations from a predefined measurable set of initial states, the nonstationary parameters of a process under external actions vary within given ranges. The control process depends on the mismatch signal and its derivatives up to an order that is less by one than the order of the equations of motion. The associated system of differential equations contains discontinuous time-dependent coefficients. A relationship is established between the eigenvalues of the quadratic forms of the given Lyapunov functions and the technical-stability criteria. It is shown that the technical-stability conditions do not depend on the sliding-mode conditions for variable-structure control processes     

Soot Particle Size Distribution Functions in a Turbulent Non-Premixed Ethylene-Nitrogen Flame

A scanning mobility particle sizer with a nano differential mobility analyzer was used to measure nanoparticle size distribution functions in a turbulent non-premixed flame. The burner utilizes a premixed pilot flame which anchors a C₂H₄/N₂ (35/65) central jet with Re _D = 20,000. Nanoparticles in the flame were sampled through a N₂-filled tube with a 500- μm orifice. Previous studies have shown that insufficient dilution of the nanoparticles can lead to coagulation in the sampling line and skewed particle size distribution functions. A system of mass flow controllers and valves were used to vary the dilution ratio. Single-stage and two-stage dilution systems were investigated. A parametric study on the effect of the dilution ratio on the observed particle size distribution function indicates that particle coagulation in the sampling line can be eliminated using a two-stage dilution process. Carbonaceous nanoparticle (soot) concentration particle size distribution functions along the flame centerline at multiple heights in the flame are presented. The resulting distributions reveal a pattern of increasing mean particle diameters as the distance from the nozzle along the centerline increases.     

变截面梁动力响应的频域传递矩阵计算方法

为提高变截面梁振动分析的计算效率,提出了基于频域传递矩阵法的动力计算算法．首先,选择线速度、角速度、弯矩和剪力作为求解变量,通过Laplace变换将变截面梁的动力响应时域偏微分方程转换为频域常微分方程;然后,通过求解频域方程并结合协调和边界条件建立变截面梁的频域传递矩阵;通过构造傅里叶级数展开形式的时域响应函数,对变截面梁传递矩阵方法求解的频响函数进行Laplace逆变换,建立了变截面梁的固有特性计算和时域瞬态响应计算方法,最后,借助数值仿真软件,开发了变截面梁动力响应分析的计算程序．完成对算例的仿真计算和分析,并与有限元计算结果进行对比,数值结果验证了该方法的正确性和有效性．     

Oil Displacement for One-Dimensional Three-Phase Flow with Phase Change in Fractured Media

In this article, the numerical simulations for one-dimensional three-phase flows in fractured porous media are implemented. The simulation results show that oil displacement in matrix is dominated by oil–water capillary pressure only under certain conditions. When conditions are changed to decrease the amount of water entering into the fractured media from the boundary of the flow field, water in fracture may be vaporized to superheated steam. In these cases, the appearance of superheated steam in fracture rather than in matrix will decrease the fracture pressure and generate the pressure difference between matrix and fracture, which results in oil flowing from matrix to fracture. Assuming that oil is wetting to steam, the matrix steam–oil capillary pressure will decrease the matrix oil-phase pressure as the matrix steam saturation increases. After the steam–oil capillary pressure finally exceeds the pressure difference due to the appearance of superheated steam in fracture, the oil displacement in matrix will stop. It is also shown that variations of the water relative permeability curve in matrix do not result in different mechanisms for oil displacement in matrix. The simulation results suggest that the amount of liquid water supply from the boundary of flow field fundamentally influence the mechanisms for oil displacement in matrix.