考虑二次梯度项影响的三孔双渗模型渗流特征

传统的试井模型与物质平衡方程都是不一致的.在非线性偏微分方程中,根据弱可压缩液体的假设,忽略二次梯度项,在试井较长时间内将产生误差.由于二次梯度项的存在,描述多孔介质微可压缩液体流动的压力分布的偏微分方程是非线性的.本文对三重介质渗流特征系统,保留了非线性偏微分方程中的所有项,没有忽略二次梯度项,建立了由基岩系统、裂缝系统及溶洞系统组成的考虑井筒储集和表皮效应对压力影响的三孔双渗模型.采用一阶向前差分,二阶中心差分的方法获得了此类三重介质模型的差分方程,用解非线性方程组的Newton迭代法求得了方程组的数值解.分别讨论了三重介质的参数变化时的无因次压力的变化规律,做出了典型无因次压力曲线图.     

分形煤层气藏数值模拟研究——Douglas-Jones预估-校正法的应用

传统的煤层气动力学模型多建立在欧几里得基础上,难以描述煤层气孔隙结构的复杂性和形状的不规则性。为此,以分形理论为基础,通过引入煤层基质和裂缝的分形维数来刻画煤层气孔隙结构的复杂性和吸附特性,建立了双重分形介质渗流模型,采用Douglas-Jones预估-校正法对非线性方程组进行离散,获得了无限大地层和有限地层定产量生产时拟稳态吸附模型的差分方程,求得数值解。结果表明,Douglas-Jones预估-校正法可以有效解决这类非线性模型的求解问题,获得无限大地层定产量生产时变形双重分形介质模型的数值解;分析各种分形参数下的煤层压力动态,做出了典型压力曲线图。对无限大地层,初期分形维数对压力影响不大,后期分形维数越小,压力越高。对有限地层,初期分形维数的影响明显,且分形维数越大,压力越低。压力随分形指数的减小量呈现先增大后减小的趋势,在末期压力平稳趋向同一值。     

分形油藏不稳定渗流问题的精确解

研究了分形油藏无限大地层和有界地层渗流模型，引入了一类有限广义Ｈａｎｋｅｌ变换，利用这种变换和Ｗｅｂｅｒ变换，在井底定流量和定压生产时，对无限大地层及有界地层（包括封闭和定压地层）六种情况，求得了实空间解析解用双参数（ｄｆ，ｄｓ）来刻画分形油藏的分形特性，分析了分形油藏压力动态特征以及分形参数和边界对压力动态的影响     

具有井筒储集的变形介质双渗模型的压力分析

为更好地研究碳酸盐油藏和低渗油的渗流问题,引入渗透率模数,考虑应力敏感地层中介质的变形,介质的双孔隙度、双渗秀率特征,同时考虑井筒储集的影响,建立新的数学模型。渗透率依赖于孔隙压力变化的流动方程是强非线性的,模型采用Douglas-Jones预估-校正法获得了无限大地层及有界封闭地层的数值解,形成了新的理论图版,并利用这些图版对模型中的有关参数进行了敏感性分析。     

具有井简储集的变形介质双渗模型的压力分析

为更好地研究碳酸盐油藏和低渗油的渗流问题，引入渗透率模数，考虑应力敏感地层中介质的变形，介质的双孔隙度、双渗秀率特征，同时考虑井筒储集的影响，建立新的数学模型。渗透率依赖于孔隙压力变化的流动方程是强非线性的，模型采用Douglas—Jones预估-校正法获得了无限大地层及有界封闭地层的数值解，形成了新的理论图版，并利用这些图版对模型中的有关参数进行了敏感性分析。     

考虑二次梯度项及动边界的双重介质低渗透油藏流动分析

在传统试井模型的非线性偏微分方程中根据弱可压缩流体的假设,忽略了二次梯度项,对于低渗透油藏这种方法是有疑问的.低渗透问题一个显著的特点就是流体的流动边界随着时间不断向外扩展.为了更好地研究双重介质低渗透油藏中流体的流动问题,考虑了二次梯度项及活动边界的影响,同时考虑了低渗透油藏的非达西渗流特征,建立了双重介质低渗透油藏流动模型.采用Douglas-Jones预估-校正差分方法获得了无限大地层定产量生产时模型的数值解,分别讨论了不同参数变化时压力的变化规律及活动边界随时间的传播规律,还分析了考虑和忽略二次梯度项影响时模型数值解之间的差异随时间的变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

变形双重介质广义流动分析

对于碳酸盐油藏和低渗油藏的渗流问题，传统的研究方法都是假设地层渗透率是常数，这假设，对于地层渗透率是压力敏感的情况，对压力的空间变化和瞬时变化将导致较大的误差。本文研究了应力敏感地层中双重介质渗流问题的压力不稳定响应，不仅考虑了储层的双重介质特征，而且考虑了应力敏感地层中介质的变形，建立了应力敏感地层双重介质的数学模型，渗透率依赖于孔隙压力变化的流动方程是强非线性的，采用Douglas-Jones预估－校正法获得了只有裂缝发生形变定产量生产时无限大地层的数值解及定产量生产岩块与裂隙同时发生形变时无限大地层的数值解，并探讨了变形参数和双重介质参数变化时压力的变化规律，给出几种情况下典型压力曲线图版，这些结果可用于实际试井分析。     

基于渗透率连续变化的低渗透多孔介质非线性渗流模型

基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合．验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明：非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。     

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

多重介质中弹性渗流的数值模拟

本文用数值模拟系统地研究了多重介质弹性径向渗流问题。为了反映井底附近压力分布的特点,按等比级数划分网格。用跳蛙式算法解3(n 1)个隐式差分方程组。计算结果是可信的,因为:压力的面积平均值与相应的解析解符合得很好,最大误差不超过3％;当井工作时间较短时,多重介质井底压力计算值与均质无限地层解析解符合得相当好;当时间较长时,多重介质计算压力进入拟定常状态。用CJ-709机计算绘制了29幅各种情况的压力分布和流量变化图。计算结果还给出“在定产的条件下,有限封闭圆形地层中各点的压力,经过相同的时间同时到达拟定常状态”这一有意义的结论。     

Exact solutions for nonlinear transient flow model including a quadratic gradient term

The models of the nonlinear radial flow for the infinite and finite reservoirs including a quadratic gradient term were presented. The exact solution was given in real space for flow equation including quadratic gradiet term for both constant-rate and constant pressure production cases in an infinite system by using generalized Weber transform.Analytical solutions for flow equation including quadratic gradient term were also obtained by using the Hankel transform for a finite circular reservoir case. Both closed and constant pressure outer boundary conditions are considered. Moreover, both constant rate and constant pressure inner boundary conditions are considered. The difference between the nonlinear pressure solution and linear pressure solution is analyzed. The difference may be reached about 8% in the long time. The effect of the quadratic gradient term in the large time well test is considered.     

拱桥悬臂施工过程中面内特征值的传递矩阵法

采用传递矩阵法对拱桥悬臂施工过程中的面内特征值问题进行求解,建立了该类桥型施工过程中面内竖弯刚度的评估方法.首先,将索和拱分别视为无垂度的张紧弦和欧拉伯努利梁,基于传递矩阵法基本理论推导了系统的总传递矩阵,考虑拱和索的边界条件以及索拱节点的位移连续性条件得到系统的特征值方程,进而计算出系统的频率和模态.同时,采用有限元...     

Three-dimensional analysis of the flow in a curved hydraulic turbine draft tube

The three-dimensional turbulent flow in a curved hydraulic turbine draft tube is studied numerically. The analysis is based on the steady Reynolds-averaged Navier–Stokes equations closed with the κ-ε model. The governing equations are discretized by a conservative finite volume formulation on a non-orthogonal body-fitted co-ordinate system. Two grid systems, one with 34 × 16 × 12 nodes and another with 50 × 30 × 22 nodes, have been used and the results from them are compared. In terms of computing effort, the number of iterations needed to yield the same degree of convergence is found to be proportional to the square root of the total number of nodes employed, which is consistent with an earlier study made for two-dimensional flows using the same algorithm. Calculations have been performed over a wide range of inlet swirl, using both the hybrid and second-order upwind schemes on coarse and fine grids. The addition of inlet swirl is found to eliminate the stalling characteristics in the downstream region and modify the behaviour of the flow markedly in the elbow region, thereby affecting the overall pressure recovery noticeably. The recovery factor increases up to a swirl ratio of about 0˙75, and then drops off. Although the general trends obtained with both finite difference operators are in agreement, the quantitative values as well as some of the fine flow structures can differ. Many of the detailed features observed on the fine grid system are smeared out on the coarse grid system, pointing out the necessity of both a good finite difference operator and a good grid distribution for an accurate result.     

Work optimization in continuous and discrete systems with complex fluids

We extend to the realm of complex fluids and finite durations the classical problem of extremum work delivered from (or consumed by) the nonequilibrium system composed of a complex fluid, a perfect thermal machine and the environment or an infinite reservoir. The fluid constitutes a valuable resource of a finite flow or amount “a finite resource”, and work production (consumption) takes place sequentially, in stages of “endoreversible” thermal machines. At each stage, heat and mass transfer takes place in boundary layers which play the role of resistances in the system model. For the fluid at flow, total specific work is extremized at constraints which take into account dynamics of heat and mass transport and rate of work generation. Finite rate limits are obtained for the work production and consumption, which provide stronger bounds that those predicted by classical thermodynamics. Optimal work functions, which incorporate an inevitable minimum of the entropy production are found as functions of end states, duration and (in discrete processes) number of stages. Formal analogies between the entropy production expressions for work-assisted and conventional mass transfer operations help formulate optimization models.     

Unsteady flow of a second-order fluid between concentric cylinders

The unsteady motion of an incompressible second-order fluid contained between two finite coaxial cylinders is examined when the outer cylinder is held fixed while the inner cylinder is constrained to execute an arbitrary angular velocity. A solution is obtained in closed form with the aid of transforms and an expression is obtained for the couple experienced. The particular case of a periodic angular velocity is then examined and some numerical work done. There is a marked difference between the results obtained and their classical counterparts.     

Analog circuit design and optimal synchronization of a modified Rayleigh system

This paper addresses the problem of optimization of the synchronization of a chaotic modified Rayleigh system. We first introduce a four-dimensional autonomous chaotic system which is obtained by the modification of a two-dimensional Rayleigh system. Some basic dynamical properties and behaviors of this system are investigated. An appropriate electronic circuit (analog simulator) is proposed for the investigation of the dynamical behavior of the proposed system. Correspondences are established between the coefficients of the system model and the components of the electronic circuit. Furthermore, we propose an optimal robust adaptive feedback which accomplishes the synchronization of two modified Rayleigh systems using the controllability functions method. The advantage of the proposed scheme is that it takes into account the energy wasted by feedback coupling and the closed loop performance on synchronization. Also, a finite horizon is explicitly computed such that the chaos synchronization is achieved at an established time. Numerical simulations are presented to verify the effectiveness of the proposed synchronization strategy. Pspice analog circuit implementation of the complete master–slave controller system is also presented to show the feasibility of the proposed scheme.     

Chebyshev super spectral viscosity solution of a two‐dimensional fluidized‐bed model

The numerical solution of a model describing a two‐dimensional fluidized bed by a Chebyshev super spectral viscosity (SSV) method is considered. The model is in the form of a hyperbolic system of conservation laws with a source term, coupled with an elliptic equation for determining a stream function. The coupled elliptic equation is solved by a finite‐difference method. The mixed SSV/finite‐difference method produces physically shaped bubbles, on a very coarse grid. Fine scale details, which were not present in previous finite‐difference solutions, are present in the solution. Copyright © 2003 John Wiley & Sons, Ltd.     

A layerwise boundary integral equation model for layers and layered media

A hybrid method is presented for the analysis of layers, plates, and multilayered systems consisting of isotropic and linear elastic materials. The problem is formulated for the general case of a multilayered system using a total potential energy formulation. The layerwise laminate theory of Reddy is employed to develop a layerwise, two-dimensional, displacement-based, hybrid boundary element model that assumes piecewise continuous distribution of the displacement components through the system's thickness. A one-dimensional finite element model is used for the analysis of the multilayered system through its thickness, and integral Fourier transforms are used to obtain the exact solution for the in-plane problem. Explicit expressions are obtained for the fundamental solution of a typical infinite layer (element) assuming linear displacement distribution through its thickness. This fundamental solution is given in a closed form in the cartesian space, and it can be applied in the two-dimensional boundary integral equation model to analyze layered structures with finite dimensions. The proposed method provides a simple, efficient, and versatile model for a three-dimensional analysis of thick plates or multilayered systems.     

Depth‐integrated nonhydrostatic free‐surface flow modeling using weighted‐averaged equations

In this study, a depth‐integrated nonhydrostatic flow model is developed using the method of weighted residuals. Using a unit weighting function, depth‐integrated Reynolds‐averaged Navier‐Stokes equations are obtained. Prescribing polynomial variations for the field variables in the vertical direction, a set of perturbation parameters remains undetermined. The model is closed generating a set of weighted‐averaged equations using a suitable weighting function. The resulting depth‐integrated nonhydrostatic model is solved with a semi‐implicit finite‐volume finite‐difference scheme. The explicit part of the model is a Godunov‐type finite‐volume scheme that uses the Harten‐Lax‐van Leer‐contact wave approximate Riemann solver to determine the nonhydrostatic depth‐averaged velocity field. The implicit part of the model is solved using a Newton‐Raphson algorithm to incorporate the effects of the pressure field in the solution. The model is applied with good results to a set of problems of coastal and river engineering, including steady flow over fixed bedforms, solitary wave propagation, solitary wave run‐up, linear frequency dispersion, propagation of sinusoidal waves over a submerged bar, and dam‐break flood waves.     

A finite difference method for elastic wave scattering by a planar crack with contacting faces

This paper presents a finite difference time-domain technique for 2D problems of elastic wave scattering by cracks with interacting faces. The proposed technique introduces cracks into the finite difference model using a set of split computational nodes. The split-node pair is bound together when the crack is closed while the nodes move freely when open, thereby a unilateral contact condition is considered. The development of the open/close status is determined by solving the equation of motion so as to yield a non-negative crack opening displacement. To check validity of the proposed scheme, 1D and 2D scattering problems for which exact solutions are known are solved numerically. The 1D problem demonstrates accuracy and stability of the scheme in the presence of the crack-face interaction. The 2D problem, in which the crack-face interaction is not considered, shows that the proposed scheme can properly reproduce the stress singularity at the tip of the crack. Finally, scattered fields from cracks with interacting faces are investigated assuming a stick and a frictionless contact conditions. In particular, the directivity and higher-harmonics are investigated in conjunction with the pre-stress since those are the basic information required for a successful ultrasonic testing of closed cracks.