Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Study of the pressure transient behavior of directional wells considering the effect of non-uniform flux distribution

During the production,the fluid in the vicinity of the directional well enters the wellbore with different rates,leading to non-uniform flux distribution along the directional well.However,in all existing studies,it is oversimplified to a uniform flux distribution,which can result in inaccurate results for field applications.Therefore,this paper proposes a semi-analytical model of a directional well based on the assumption of non-uniform flux distribution.Specifically,the direction well is discretized into a carefully chosen series of linear sources,such that the complex well trajectory can be captured and the nonuniform flux distribution along the wellbore can be considered to model the three-dimensional flow behavior.By using the finite difference method,we can obtain the numerical solutions of the transient flow within the wellbore.With the aid of Green's function method,we can obtain the analytical solutions of the transient flow from the matrix to the wellbore.The complete flow behavior of a directional well is perfectly represented by coupling the above two types of transient flow.Subsequently,on the basis of the proposed model,we conduct a comprehensive analysis of the pressure transient behavior of a directional well.The computation results show that the flux variation along the direction well has a significant effect on pressure responses.In addition,the directional well in an infinite reservoir may exhibit the following flow regimes:wellbore afterflow,transition flow,inclined radial flow,elliptical flow,horizontal linear flow,and horizontal radial flow.The horizontal linear flow can be observed only if the formation thickness is much smaller than the well length.Furthermore,a dip region that appears on the pressure derivative curve indicates the three-dimensional flow behavior near the wellbore.

Influence of enlarged section parameters on pressure transients of high-speed train passing through a tunnel

The influence of enlarged section parameters on pressure transients of high-speed train passing through a tunnel was investigated by numerical simulation.The calculation results obtained by the structured and unstructured grid and the experimental results of smooth wall tunnel were verified.Numerical simulation studies were conducted on three tunnel enlarged section parameters,the enlarged section distribution along circumferential direction,the enlarged section area and the enlarged section distribution along tunnel length direction.The calculation results show that the influence of the different enlarged section distributions along tunnel circumferential direction on pressure transients in the tunnel is basically consistent.There is an optimal enlarged section area for the minimum value of the pressure variation amplitude and the average pressure variation in the tunnel.The law of the pressure variation amplitude and the average pressure variation of the enlarged section distribution along tunnel length direction are obtained.

A review of pressure transient analysis in reservoirs with natural fractures,vugs and/or caves

A review of the pressure transient analysis of flow in reservoirs having natural fractures,vugs and/or caves is presented to provide an insight into how much knowledge has been acquired about this phenomenon and to highlight the gaps still open for further research.A comparison-based approach is adopted which involved the review of works by several authors and identifying the limiting assumptions,model restrictions and applicability.Pressure transient analysis provides information to aid the identification of important features of reservoirs.It also provides an explanation to complex reservoir pressuredependent variations which have led to improved understanding and optimization of the reservoir dynamics.Pressure transient analysis techniques,however,have limitations as not all its models find application in naturally fractured and vuggy reservoirs as the flow dynamics differ considerably.Pollard’s model presented in 1953 provided the foundation for existing pressure transient analysis in these types of reservoirs,and since then,several authors have modified this basic model and come up with more accurate models to characterize the dynamic pressure behavior in reservoirs with natural fractures,vugs and/or caves,with most having inherent limitations.This paper summarizes what has been done,what knowledge is considered established and the gaps left to be researched on.

Pressure Transient Analysis of Arbitrarily Shaped Fractured Reservoirs

Reservoir boundary shape has a great influence on the transient pressure response of oil wells located in arbitrarily shaped reservoirs. Conventional analytical methods can only be used to calculate transient pressure response in regularly shaped reservoirs. Under the assumption that permeability varies exponentially with pressure drop, a mathematical model for well test interpretation of arbitrarily shaped deformable reservoirs was established. By using the regular perturbation method and the boundary element method, the model could be solved. The pressure behavior of wells with wellbore storage and skin effects was obtained by using the Duhamel principle. The type curves were plotted and analyzed by considering the effects of permeability modulus, arbitrary shape and impermeable region.

Simultaneous inversion of permeability,skin and boundary from pressure transient test data in three-dimensional single well reservoir model

This study presents a novel approach for simultaneous inversion of the key reservoir parameters like horizontal permeability,vertical permeability,skin,and boundary distances for spatial distribution across the grid cells in a 3D single well reservoir model(SWRM).These parameters are first estimated from the standard pressure transient analysis of well test pressure and rate data,which also act as a priori for the inverse problem.A field-worthy layer cake geological model is prepared based on the prior information obtained from pressure transient analysis,followed by a sequential flow simulation of field well test operation.The simulation results provide the model pressure versus rate data as the synthetic data for this study.A cost function is defined incorporating the well test pressure data and model pressure data,which would determine the convergence.The inversion process is to optimize the spatial distribution of reservoir parameters to minimize the difference between the measured pressure transient data and the modelled one,which is obtained from the multiphase fluid flow simulator that solves the implicit black-oil fluid-flow diffusivity equations at every step.A Gauss-Newton(GN)inversion scheme is used for the inversion.The reliability of inversion results depends on the accuracy of priori reservoir parameters fed to the solver,which can be refined if required through uncertainty parameter optimization(UPO).This approach helps to obtain a faster and reliable update of reservoir parameters in a layer cake homogeneous geomodel,hereby introducing the required heterogeneity.This increases the confidence and reliability of a geomodel,which is further used for various production prediction strategies.

A Transient-Pressure-Based Numerical Approach for Interlayer Identification in Sand Reservoirs

Almost all sandstone reservoirs contain interlayers. The identification and characterization of these interlayers iscritical for minimizing the uncertainty associated with oilfield development and improving oil and gas recovery.Identifying interlayers outside wells using identification methods based on logging data and machine learning isdifficult and seismic-based identification techniques are expensive. Herein, a numerical model based on seepageand well-testing theories is introduced to identify interlayers using transient pressure data. The proposed modelrelies on the open-source MATLAB Reservoir Simulation Toolbox. The effects of the interlayer thickness, position,and width on the pressure response are thoroughly investigated. A procedure for inverting interlayer parametersin the reservoir using the bottom-hole pressure is also proposed. This method uses only transient pressuredata during well testing and can effectively identify the interlayer distribution near the wellbore at an extremelylow cost. The reliability of the model is verified using effective oilfield examples.

Well interference evaluation considering complex fracture networks through pressure and rate transient analysis in unconventional reservoirs

Severe well interference through complex fracture networks(CFNs)can be observed among multi-well pads in low permeability reservoirs.The well interference analysis between multi-fractured horizontal wells(MFHWs)is vitally important for reservoir effective development.Well interference has been historically investigated by pressure transient analysis,while it has shown that rate transient analysis has great potential in well interference diagnosis.However,the impact of complex fracture networks(CFNs)on rate transient behavior of parent well and child well in unconventional reservoirs is still not clear.To further investigate,this paper develops an integrated approach combining pressure and rate transient analysis for well interference diagnosis considering CFNs.To perform multi-well simulation considering CFNs,non-intrusive embedded discrete fracture model approach was applied for coupling fracture with reservoir models.The impact of CFN including natural fractures and frac-hits on pressure and rate transient behavior in multi-well system was investigated.On a logelog plot,interference flow and compound linear flow are two new flow regimes caused by nearby producers.When both NFs and frac-hits are present in the reservoir,frac-hits have a greater impact on well#1 which contains frac-hits,and NFs have greater impact on well#3 which does not have frac-hits.For all well producing circumstances,it might be challenging to see divergence during pseudosteady state flow brought on by frac-hits on the logelog plot.Besides,when NFs occur,reservoir depletion becomes noticeable in comparison to frac-hits in pressure distribution.Application of this integrated approach demonstrates that it works well to characterize the well interference among different multi-fractured horizontal wells in a well pad.Better reservoir evaluation can be acquired based on the new features observed in the novel model,demonstrating the practicability of the proposed approach.The findings of this study can help for better evaluating well interference degree in multi-well systems combing PTA and RTA,which can reduce the uncertainty and improve the accuracy of the well interference analysis based on both field pressure and rate data.

PRESSURE TRANSIENT ANALYSIS OF HETEROGENOUS RESERVOIRS WITH IMPERMEABILITY BARRIER USING PERTURBATION BOUNDARY ELEMENT METHOD

The transient flow mathematical model of arbitrary shaped heterogeneous reservoirs with impermeability barrier is proposed in this paper. In order to establish this model, the perturbation method is employed and the solution of model is expanded into a series in powers of perturbation parameter. By using the Boundary Element Method (BEM) and Duhamel principle, wellbore pressure with effects of skins and wellbore storage is obtained. The type curves are plotted and analyzed considering effects of heterogeneity, arbitrary shape and impermeable barriers. Finally, the results obtained by perturbation boundary element method is compared with the analytical solution and is available for the transient pressure analysis of arbitrary shaped reservoirs.

Wireless distributed test system based on transient pressure signal detection and recognition

During the test on transient pressure signal in explosion field,false trigger caused by field interference can lead to test failure.To improve the stability of test system,a signal detection and recognition technology is proposed for transient pressure test system.In the process of signal acquisition,firstly,electrical levels are monitored in real time to find effective abrupt changes and mark them;then the effective data segments are detecdted totected;thus the effective signals can be acquired in turn finally.The experimental results show that the shock wave signal can be collected effectively and the reliability of the test system can be improved after removal of interferences.

Design and implementation of configurable parameters test system for transient pressure measurement

For the shock wave overpressure signal measurement of explosion field, the type and equivalent of the ammunition as well as the distance from test nodes to ground zero have a great influence on the characteristics of shock wave signal. To ensure that shock wave signal with different characteristic parameters can be collected completely in the test of shooting range, a tran- sient pressure signal storage system with configurable parameters is presented. In order to meet the requirements of different test nodes, the system can set parameters of all test nodes by wireless communication in explosion field. The feasibility and re- liability of the system are fully verified through transient pressure measurement in explosion field.

Analysis on setting airshaft at mid-tunnel to reduce transient pressure variation

Setting airshafl is one of the methods efficient to control the aerodynamic effects on a train traveling through a tunnel. By using numeral simulation, analysis is carried out on the effect of airshafl on transient pressures generated in cabin. After setting airshaft, the magnitude of pressure fluctuation in cabin is reduced nearly 40%. By analyzing the role of airshaft in alleviating pressure fluctuation, a formula to determine the optimal airshaft position is deduced.

An efficient algorithm to compute transient pressure responses of slanted wells with arbitrary inclination in reservoirs

Compared with vertical and horizontal wells, the solution and computation of transient pressure responses of slanted wells are more complex. Vertical and horizontal wells are both simplified cases of slanted wells at particular inclination, so the model for slanted wells is more general and more complex than other models for vertical and horizontal wells. Many authors have studied unsteady-state flow of fluids in slanted wells and various solutions have been proposed. However, until now, few of the published results pertain to the computational efficiency. Whether in the time domain or in the Laplace domain, the computation of integration of complex functions is necessary in obtaining pressure responses of slanted wells, while the computation of the integration is complex and time-consuming. To obtain a perfect type curve the computation time is unacceptable even with an aid of high-speed computers. The purpose of this paper is to present an efficient algorithm to compute transient pressure distributions caused by slanted wells in reservoirs. Based on rigorous derivation, the transient pressure solution for slanted wells of any inclination angle is presented. Assuming an infinite-conductivity wellbore, the location of the equivalent-pressure point is determined. More importantly, according to the characteristics of the integrand in a transient pressure solution for slanted wells, the whole integral interval is partitioned into several small integral intervals, and then the method of variable substitution and the variable step-size piecewise numerical integration are employed. The amount of computation is significantly reduced and the computational efficiency is greatly improved. The algorithm proposed in this paper thoroughly solved the difficulty in the efficient and high-speed computation of transient pressure distribution of slanted wells with any inclination angle.

Influence of transient intraocular pressure elevation during laser in situ keratomileusis on rabbit retina thickness

AIMTo utilize tissue micro measurement to study the effect of transient high intraocular pressure (IOP) induced by different durations of suction during laser in situ keratomileusis (LASIK) on rabbit retina thickness.METHODSSixty healthy New Zealand white rabbits were randomly divided into a control group, and 3 negative-pressure suction groups (20s group, 45s group, and 3min group) and each group was comprised of 15 rabbits (30 eyes); the latter 3 groups were the transient high IOP models. The retinal tissue around the papilledema was separated. Hematoxylin and eosin (HE) staining was carried out to generate slices for light microscopy. The changes in the retina thickness values of each layer were measured for all animals in each group at different postoperative recovery periods and compared with the values recorded for the animals in the control group. The thickness of the retinal tissue showed a normal distribution. The ANOVA was performed by using SPSS13.0 statistic software.RESULTSIn the comparison between the 20s and 45s negative-pressure suction groups and the control group, no significant differences were observed, except at 14d. Significant difference was observed between the 3min negative-pressure suction group and the control group, and the retina thickness value of each layer reached a peak at 14d after repair.CONCLUSIONConventional negative suction during LASIK may not lead to significant changes in retinal tissue thickness; however, if the suction duration is increased to 3min, it will cause significant changes in retinal tissue thickness.

Processing and analysis of transient pressure measurements from permanent down-hole gauges

With permanent down-hole gauges （PDGs） widely installed in oilfields around the world in recent years, a continuous stream of transient pressure data in real time is now available, which motivates a new round of research interests in further developing pressure transient processing and analysis techniques. Transient pressure measurements from PDG are characterized by long term and high volume data. These data are recorded under unconstrained circumstances, so effects due to noise, rate fluctuation and interference from other wells cannot be avoided. These effects make the measured pressure trends decline or rise and then obscure or distort the actual flow behavior, which makes subsequent analysis difficult. In this paper, the problems encountered in analysis of PDG transient pressure are investigated. A newly developed workflow for processing and analyzing PDG transient pressure data is proposed. Numerical well testing synthetic studies are performed to demonstrate these procedures. The results prove that this new technique works well and the potential for practical application looks very promising.

PRETTA：A COMPUTER PROGRAM FOR PWR PRESSURIZER’S TRANSIENT THERMODYNAMICS

A computer program PRETTA “Pressurizer Transient Thermodynamics Analysis” was developed for the prediction of pressurizer under transient conditions. It is based on the solution of the conservation laws of heat and mass applied to the three separate and non equilibrium thermodynamic regions. In the program all of the important thermal-hydraulics phenomena occurring in the pressurizer: stratification of the hot water and incoming cold water, bulk flashing and condensation, wall condensation, and interfacial heat and mass transfer have been considered. The bubble rising and rain-out models are developed to describe bulk flashing and condensation, respectively. To obtain the wall condensation rate, a one-dimensional heat conduction equation is solved by the pivoting method. The presented computer program will predict the pressure-time behavior of a PWR pressurizer during a variety of transients. The results obtained from the proposed mathematical model are in good agreement with available data on the CHASHMA nuclear power plant's pressurizer performance.

THE TRANSIENT TWO-DIMENSIONAL FLOW THROUGH DOUBLE POROUS MEDIA

This paper presents the analytical solutions in Laplace domain for two-dimensional nonsteady flow of slightly compressible liquid in porous media with double porosity by using the methods of integral transforms and variables separation. The effects of the ratio of storativities to , interporosity flow parameter on the pressure behaviors for a vertically fractured well with infinite conductivity are investigated by using the method of numerical inversion. The new log-log diagnosis graph of the pressures is given and analysed.

An approximate analytical solution for transient gas flows in a vertically fractured well of finite fracture conductivity

An analytical solution in physical variable space is presented for transient gas flows during constant-rate production from a vertically-fractured well in an infinite homogeneous reservoir with finite fracture conductivity.The solution is based on the short-time asymptotic solution and a new approximate transient elliptical flow solution,which covers transient flows from the bilinear flow regime to the pseudo-radial flow regime.The solution covers the well-known asymptotic solutions in both short-and long-time limits of bilinear and pseudo-radial flows.The analytical model provides a practical and reliable engineering tool to evaluate the fractured reservoir properties,which can be programmed using a spreadsheet.

Research on numerical simulation of transient pressure for high-speed train passing through the most unfavourable length tunnel

The length of high-speed railway tunnel is an important factor affecting transient pressure of high-speed train.When the tunnel length is the most unfavourable,the transient pressure changes in the tunnel and on the surface of the train are the most severe,which may affect the safe operation of the train or damage the structure in the tunnel.Based on the three-dimensional,compressible,unsteady N-S equation and finite volume method,this paper uses the CFD numerical simulation method to study the change and amplitude distribution of the transient pressure on the train surface and the tunnel when a high-speed train passes through the most unfavourable length tunnel.A fast calculation method is proposed to save the cost of calculation;it has great applicability of pressure amplitude.The results show that the pressure distribution in the tunnel and on the surface of the train is affected by the train speed,the length of the train and the position of the measuring point.The minimum negative peak value in the tunnel appears at the position where the superposition phenomenon is the most severe,and the position will change with the speed of the train.There are two negative peak waveforms of the train surface pressure,and the first waveformn is greatly affected by the train speed.It improves a reference for studying the strength requirement of the most unfavourable length tunnels and trains,and ensures the safe operation of trains in tunnels of different lengths.

Water hammer prediction and control:the Green's function method

By Green＇s function method we show that the water hammer （WH） can be analytically predicted for both laminar and turbulent flows （for the latter, with an eddy vis- cosity depending solely on the space coordinates）, and thus its hazardous effect can be rationally controlled and mini- mized. To this end, we generalize a laminar water hammer equation of Wang et al. （J. Hydrodynamics, B2, 51, 1995） to include arbitrary initial condition and variable viscosity, and obtain its solution by Green＇s function method. The pre- dicted characteristic WH behaviors by the solutions are in excellent agreement with both direct numerical simulation of the original governing equations and, by adjusting the eddy viscosity coefficient, experimentally measured turbulent flow data. Optimal WH control principle is thereby constructed and demonstrated.