MODEL AND METHOD OF WELL TEST ANALYSIS FOR WELLS WITH VERTICAL FRACTURE

NomenclatureB———volumecoefficientofformationfluid,dimensionless;C———wellborestoragecoefficient,m3 MPa;Cft ———comprehensivecompressioncoefficientoffluidinfracture ,1 MPa;Ct ———comprehensivecompressioncoefficientofformationfluid,1 MPa ;h———heightofformation ,m ;K———formationpermeability ,μm2 ;Kf ———fracturepermeability ,μm2 ;KfWf ———fractureconductivity ,μm2 ·cm ;Lf ———halfoffracturelength,m ;N———thetotalnumberofselectedtimepoint,dimensionless;P———for…     

A triangular quasi-conforming finite element for transient dynamic analysis

A type of 3 node triangular element is constructed by the Quasi-conforming method, which may be used to solve the equation of a type of inverse problem of wave propagation after Laplace transformation Δu−A ² u=0. The strains in the element are approximated by an exponential function and the string-net function between neighbouring elements is approximated by one dimensional general solution of the equation. Furthermore the strain field satisfies the equation, and therefore in the derivation of the element formulation, no shape function is needed. In this sense, it is a kind of hybrid element. Compared with the ordinary linear triangular element, the new one features higher precision with coarse meshes. Some numerical tests are presented. The project is supported by the National Natural Science Foundation of China.     

A numerical method for solving equilibrium problems of masonry-like solids

This paper proposes a numerical method for the solution of equilibrium problems of solids which do not support tension. Some boundary-value problems are solved numerically and the solution obtained is compared to the exact one.

---

Sommario In questo lavoro viene proposto un metodo numerico per la soluzione di problemi di equilibrio di solidi non resistenti a trazione. Vengono successivamente risolti numericamente alcuni problemi di equilibrio e la soluzione ottenuta è confrontata con quella esatta.

     

Effective numerical methods for elasto-plastic contact problems with friction

Several effective numerical methods for solving the elasto-plastic contact problems with friction are presented. First, a direct substitution method is employed to impose the contact constraint conditions on condensed finite element equations, thus resulting in a reduction by half in the dimension of final governing equations. Second, an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation, which distinguishes two kinds of nonlinearities, and makes the solution unique. In addition, Positive-Negative Sequence Modification Method is used to condense the finite element equations of each substructure and an analytical integration is introduced to determine the elasto-plastic status after each time step or each iteration, hence the computational efficiency is enhanced to a great extent. Finally, several test and practical examples are presented showing the validity and versatility of these methods and algorithms. The Project Supported by National Natural Science Foundation of China.     

气藏压裂水平井非稳态流动分析

应用Green函数和Newman积原理,建立了压裂水平井与气藏耦合的非稳态流动模型,并给出求解方法。所建立的模型可应用多种约束条件,考虑了井筒压降和加速度的影响,并适用于各向异性气藏。实例计算表明：水平井在早期非稳态阶段的流量比拟稳态阶段大得多。早期非稳态流动阶段,各条裂缝之间未发生干扰,流量的大小与裂缝条数成正比例关...     

Dynamic production prediction and parameter identification for gas well with vertical fracture

ＮｏｍｅｎｃｌａｔｕｒｅＢ—ｖｏｌｕｍｅｃｏｅｆｆｉｃｉｅｎｔｏｆｆｏｒｍａｔｉｏｎｆｌｕｉｄ ,ｄｉｍｅｎｓｉｏｎｌｅｓｓ ;Ｃｆｔ—ｃｏｍｐｒｅｈｅｎｓｉｖｅｃｏｍｐｒｅｓｓｉｎｇｃｏｅｆｆｉｃｉｅｎｔｏｆｆｌｕｉｄｉｎｆｒａｃｔｕｒｅ ,(ＭＰａ) - 1;Ｃｐ—ｗｅｉｇｈｔｏｆｓｉｎｇｌｅｌａｙｅｒｐｒｏｐｐａｎｔｐｅｒｕｎｉｔｆｒａｃｔｕｒｅａｒｅａ ,ｋｇ/ｍ2 ;Ｃｔ —ｃｏｍｐｒｅｈｅｎｓｉｖｅｃｏｍｐｒｅｓｓｉｎｇｃｏｅｆｆｉｃｉｅｎｔｏｆｆｏｒｍａｔｉｏｎｆｌｕｉｄ ,(ＭＰａ) - 1;ｄＢＨＮ —Ｂｒ…     

Finite element analysis for the unsteady nearshore circulation due to wave-current interaction (I)—numerical model

In this paper, a numerical model for predicting the unsteady nearshore circulation due to wave-current interaction was proposed. In addition to the traditional continuity, momentum and energy equations, the dispersion and refraction relations were included in the governing equations. Moreover, the effects of lateral shears, wind, radiation and bottom stresses were analysed in the governing equations. Therefore, we expect that this model may more completely and exactly reflect the law of wave-current interaction. In part (II) we will adopt the selective lumping two-step explicit finite element method to solve the model, and some examples will be presented.     

Flow field modelling near a well with a conductive fracture

In the technology of oil recovery the oil production rate can be increased by generation of a vertical sand-filld conductive fracture on the wall of the well. Oil diffuses through the conductive fracture to the well. In this paper the seepage flow and isothermal deformation fields in both the formation and fracture and the oil production rate at the well are studied by modelling the formation as an infinite poroelastic medium saturated with a one-phase compressible fluid. The fracture is treated as a one-dimensional poroelastic medium. Darcy flows are considered in both the formation and fracture. The plane strain condition is imposed. Our solution is obtained numerically by a finite element method based on a variational principle. The accuracy of the analysis is studied by comparison of the numerical solutions of some problems with their analytical solutions. Since we are dealing with the transient flow problem of an infinite region, an extrapolation technique is employed to find the finite element solution. The production rate of a well with the conductive fracture is compared with that of a well without the conductive fracture.     

Numerical analysis of complex flow of fibre suspensions in a viscoelastic fluid

Research efforts made so far to simulate fibre suspension flows are limited to fibre suspensions in Newtonian fluids. Though short fibre composites are mostly made of polymers, lack of suitable constitutive equations for fibre suspensions in viscoelastic fluids frustrates attempts to simulate flows of these suspensions. A preliminary work done by the author led to a constitutive equation for semiconcentrated fibre suspensions in the Oldroyd-B fluid. This paper describes the mathematical formulation of the flow problems for such a suspension and the numerical procedure to solve them. Some numerical results of flow past a sphere in a tube are also presented. This paper is supported by the University of Melbourne, the National Natural Science Foundation of China and Zhejiang Province.     

A three-dimensional transient two-phase flow analysis with a density pertubation finite volume method

In this paper, we present the analysis of a full three dimensional transient two-phase flow with strong disequilibria. This analysis uses a finite volume scheme with a recently developed density perturbation method. The application test case is a three dimensional extension of a one-dimensional experiment on blowdown pipe. The results are stable with mesh refinements and show realistic three-dimensional effects.     

A finite element model for numerical simulation of thermo-mechanical frictional contact problems

Two kinds of variational principles for numerical simulation of heat transfer and contact analysis are respectively presented. A finite element model for numerical simulation of the thermal contact problems is developed with a pressure dependent heat transfer constitutive model across the contact surface. The numerical algorithm for the finite element analysis of the thermomechanical contact problems is thus developed. Numerical examples are computed and the results demonstrate the validity of the model and algorithm developed. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (50178016, 10225212) and the Foundation for University Key Teacher by the Ministry of Education of China     

A comparison of spatial grids for numerical modelling of flows in near-coastal seas

Three-dimensional algorithms for the numerical computation of flows caused by tides or meteorological forcing are developed for four of Arakawa's spatial grid types using a spectral method in the vertical dimension. Three of the grids, in which the velocity components are computed at the same grid points, offer potential advantages over the commonly used C-grid. The computed results from the four grids are compared for three test problems based on the linearized hydrodynamical equations. It is concluded that the B-grid provides a viable alternative to the C-grid, with significant advantages when a spectral method is used.     

柱塞泵油膜压力分析的高效有限单元法

The paper concerns numerical analysis of pressure distribution of an oil film on the valve plate in the variable height gap of an axial piston pump. The analysis employs the finite element method. For determination of oil pressure variations in the gap, the Reynolds equation, commonly applied in the theory of lubrication, is applied. The equation is solved numerically with the use of self-developed program based on the finite element method. In order to obtain high accuracy of the results, an adaptive mesh refinement based on residual estimations of solution errors is applied. The calculation results are represented as dependent on the geometric and working parameters of the pump.     

A Stokes model with cavitation for the numerical simulation of hydrodynamic lubrication

We present a cavitation model based on the Stokes equation and formulate adaptive finite element methods for its numerical solution. A posteriori error estimates and adaptive algorithms are derived, and numerical examples illustrating the theory are supplied, in particular with comparison to the simplified Reynolds model of lubrication. Copyright © 2010 John Wiley & Sons, Ltd.     

基于离散裂缝的多段压裂水平井数值试井模型及应用

水平井压裂技术已经成为开发低渗透油气藏、页岩气藏和致密气场等非常规油气藏的关键技术。基于离散裂缝模型,对裂缝进行简化,建立了二维多段压裂水平井有限导流数值试井模型,利用有限元方法求解模型,获得多段压裂水平井试井理论曲线和压力场特征。分析表明：多段压裂水平井的试井理论曲线一共分为七个阶段：井筒储存段、裂缝线性流段、裂缝-地层双线性流段、裂缝干扰段、地层线性流段、系统径向流段和边界作用段,其中裂缝-地层双线性流段和裂缝干扰是其典型特征。分析了裂缝数量、裂缝间距、裂缝不对称、裂缝不等长和裂缝部分缺失等因素对试井理论曲线的影响,结果表明：裂缝数量和裂缝间距对试井理论曲线的影响最大。较多的裂缝、较大裂缝间距、对称的裂缝和等长的裂缝有利于降低压裂水平井井底的流动阻力,提高产能。将建立的数值试井模型应用于四川盆地一口多段压裂水平井的压力恢复测试的数值试井解释,结果表明：本文建立的模型可以较好的拟合压力恢复测试数据,可以获得裂缝的导流能力和裂缝长度,为压裂效果评价和压裂设计提供指导。     

A mathematical model and numerical solution technique for a novel adjustable hydrodynamic bearing

An analysis model for a novel adjustable hydrodynamic fluid film bearing is described. The principles of hydrodynamic lubrication are outlined together with an expanded version of the governing pressure field equation as related to the novel bearing. Finite difference approximations are given for the pressure field equation and a temperature model, both related to the fluid film thickness. Relationships of viscosity with temperature and pressure are included. A finite element model and an iterative computational process are described, whereby full simultaneously converged field solutions for fluid film thickness, temperature, viscosity and pressure were obtained, together with oil film forces. The model and solution process were developed to apply to a variety of hydrodynamic bearings and an outline is given of its extensive use in the design and simulation of one version of the novel bearing. Observations are given on the operation, success rates and verifications of the computational process. Copyright © 1999 John Wiley & Sons, Ltd.     

A new method for stress analysis of a cyclically symmetric structure

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

Simulation of size effects by a phase field model for fracture

In phase field fracture models the value of the order parameter distinguishes between broken and undamaged material. At crack faces the order parameter interpolates smoothly between these two states of the material, which can be regarded as phases. The crack evolution follows implicitly from the time integration of an evolution equation of the order parameter, which is coupled to the mechanical field equations. Among other phenomena phase field fracture models are able to reproduce crack nucleation in initially sound materials. For a 1D setting it has been shown that crack nucleation is triggered by the loss of stability of the unfractured, spatially homogeneous solution, and that the stability point depends on the size of the considered structure. This work numerically investigates to which extend size effects are reproduced by the 2D phase field model. Exemplarily, a finite element study of the hole size effect is performed and the simulation results are compared to experimental data.     

A study on numerical methods for transient rotating flow induced by starting blades

Based on the finite volume method, three methods for rotational region treatment were presented and validated by simulating two-dimensional accelerating rotational flows. Separate transient incompressible flows induced by cross-shaped blades during starting process were simulated using the dynamic mesh, sliding mesh and dynamic reference frame methods. The computing performance and stability of the three methods were evaluated by comparing numerical results, and the transient characteristics of the accelerating rotational flow were analysed numerically. Results showed that the starting process affected the flow structure and transient characteristics of the accelerating rotational flows. The sliding mesh method showed higher computational efficiency and accuracy compared with other methods, and could easily be extended to solve three-dimensional transient flows in hydraulic machineries under transient operations, such as start-up and shutdown.