Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

混合层强化混合的数值研究

受 Ｗａｎｇ ＆ Ｆｉｅｄｌｅｒ（１９９７）的实验的启发，采用高阶精度的差分格式，通过数值模拟的方法，研究了二维混合层及限于两平板间的二维混合层（二维受限混合层）入口处加振动对提高混合层混合效率的作用．计算结果表明：对二维混合层，振动的频率越低，在混合层中产生的大尺度涡结构的尺度越大，在频率很低时，涡具有相似性；对限于两平板间的二维混合层，在一定的振动频率下，混合层中产生的涡较大而且破碎得也较好，这将有利于混合．这一结论与 Ｗａｎｇ ＆ Ｆｉｅｄｌｅｒ（１９９７）的实验观测到的结果是一致的．     

一类非线性自治大系统的全局稳定性

本文研究非线性自治大系统 dx_i/dt=sum from i=1 to r(f_(ij)(x_j)) (i=1,…,r) 这里,x_i∈R~(ni),f_(ij)∈C(R~(nj),R~(ni),f_(ij)(o)=0,得到保证其零解为全局一致渐近稳定的充分条件(见定理1)。     

液位检测控制系统设计

本文介绍了单片机和大规模集成电路在工业液位显示、控制中的应用，解决了工业中对液位显示的复杂电路和控制上的单一性，给出了硬件电路和软件设计方法。     

A Limited-Memory Multipoint Symmetric Secant Method for Bound Constrained Optimization

A new algorithm for solving smooth large-scale minimization problems with bound constraints is introduced. The way of dealing with active constraints is similar to the one used in some recently introduced quadratic solvers. A limited-memory multipoint symmetric secant method for approximating the Hessian is presented. Positive-definiteness of the Hessian approximation is not enforced. A combination of trust-region and conjugate-gradient approaches is used to explore useful information. Global convergence is proved for a general model algorithm. Results of numerical experiments are presented.     

A class of methods for solving large,convex quadratic programs subject to box constraints

In this paper we analyze conjugate gradient-type algorithms for solving convex quadratic programs subject only to box constraints (i.e., lower and upper bounds on the variables). Programs of this type, which we denote by BQP, play an important role in many optimization models and algorithms. We propose a new class of finite algorithms based on a nonfinite heuristic for solving a large, sparse BQP. The numerical results suggest that these algorithms are competitive with Dembo and Tulowitzski's (1983) CRGP algorithm in general, and perform better than CRGP for problems that have a low percentage of free variables at optimality, and for problems with only nonnegativity constraints.     

A Subtree-Partitioning Algorithm for Inducing Parallelism in Network Simplex Dual Updates

This paper reports on thedevelopment of a very efficient method for partitioning the networksimplex basis subtree in which dual values must be updated duringa pivot. The partitioning procedure may be concurrently executedby multiple processes. The resulting rapid decomposition of thesubtree allows an arbitrary number of processes to be utilized inthe actual dual update. This approach alleviates a primary limitationof the most efficient parallel network simplex implementationpublished to date. The new code performs at least as well as theprevious implementation on medium-scale problems and reduces averagesolution time by over 34% on large-scale problems.     

A SPARSE SUBSPACE TRUNCATED NEWTON METHOD FOR LARGE-SCALE BOUND CONSTRAINED NONLINEAR OPTIMIZATION

In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices outside of the active set, while the projected gradient method is used to update the active variables. At each iterative level, the search direction consists of three parts, one of which is a subspace truncated Newton direction, the other two are subspace gradient and modified gradient directions. The subspace truncated Newton direction is obtained by solving a sparse system of linear equations. The global convergence and quadratic convergence rate of the algorithm are proved and some numerical tests are given.     

大型注水系统注水站合理排量的确定方法

大型注水管网系统结构复杂 ,涉及的井点、管网节点多 ,且压力、流量互相耦合 ,注水站排量合理分配难度大。以流体力学基本理论为基础 ,利用“压力谷”的性质 ,根据注水管网特性 ,建立了注水管网的基本模型 ,并在满足注水量和管路阻力损失最小的前提下 ,给出了一种确定大型注水系统注水站设计排量和运行排量的方法。实例计算表明 ,利用该方法调节大型注水系统注水站排量 ,可以使注水站和管网同时达到最优运行状态 ,并且能够科学地计算注水站的设计排量。该研究为大型注水系统的合理设计和优化运行提供了一定的理论指导     

The cluster-galaxy cross-correlation and the average infall velocity around cluster in the modified CDM models

The cluster-galaxy cross-correlation and the average infall velocity around cluster for the tilted CDM models, the hot-cold mixed dark matter models (MDM) and theA CDM models are calculated. And it is found that the prediction of the cluster-galaxy cross-correlation of theA CDM models with Ω_A=0.9, Ω_C=0.1,h=1,b = 1 or 2 is much higher than that of observation for the radius larger than 10 h⁻¹ Mpc (H₀ = 100 h.kms⁻¹Mpc⁻¹), and the other models are compatible with the observation on the scale of (3–20) h⁻ Mpc. The results can be used to compare with new observations to restrict the nature of the dark matter and the biasing factor. Project supported by the National Natural Science Foundation of China (Grant No. 19733001).