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...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.展开更多
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 ...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.展开更多
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 ...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.展开更多
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 pr...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.展开更多
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 res...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.展开更多
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 ga...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.展开更多
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 i...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.展开更多
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 thinni...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.展开更多
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 s...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.展开更多
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 vit...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.展开更多
This paper presents the solutions of three -dimensional unsteady flow in Laplace domain for horizontal well in double porous media with vertical feed, using the method both integral transforms and variable separation....This paper presents the solutions of three -dimensional unsteady flow in Laplace domain for horizontal well in double porous media with vertical feed, using the method both integral transforms and variable separation. The influences of parameters for double porous media ω,λ, attitudes of the horizontal well ZWD. LD and the coefficient of vertical feed a, on the dimensionless pressure behaviors of horizontal well are discussed and calculated by means of the algorithm of Laplace numerical inversion. The diagnosis and specialized graphs are plotted and analysed.展开更多
This paper presents the results of an ongoing investigation into transient pressure pulses using Shan- non entropy. Pressure fluctuations (produced by gas-solid two-phase flow during fluidized dense-phase conveying)...This paper presents the results of an ongoing investigation into transient pressure pulses using Shan- non entropy. Pressure fluctuations (produced by gas-solid two-phase flow during fluidized dense-phase conveying) are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals (Mittal, Mallick, & Wypych, 2014). Two different powders, namely fly ash (median particle diameter 45 μm, particle den- sity 1950 kg/m3. loosely poured bulk density 950 kg/m3) and cement (median particle diameter 15 p,m, particle density 3060 kg/m3, loosely poured bulk density 1070 kg/m3), are conveyed through different pipelines (51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length). The transient nature of pressure fluc- tuations (instead of steady-state behavior) is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities (i.e. 3-5 m/s) and very high velocities (i.e. 11-14 m/s) while low Shannon entropy corresponds to mid-range velocities (i.e. 6-8 m/s).展开更多
Hydraulic fracturing combined with horizontal drilling has been the technology that makes it possible to economically produce natural gas from unconventional shale gas or tight gas reservoirs.Hydraulic fracturing oper...Hydraulic fracturing combined with horizontal drilling has been the technology that makes it possible to economically produce natural gas from unconventional shale gas or tight gas reservoirs.Hydraulic fracturing operations,in particular,multistage fracturing treatments along with horizontal wells in unconventional formations create complex fracture geometries or networks,which are difficult to characterize.The traditional analysis using a single vertical or horizontal fracture concept may be no longer applicable.Knowledge of these created fracture properties,such as their spatial distribution,extension and fracture areas,is essential information to evaluate stimulation results.However,there are currently few effective approaches available for quantifying hydraulic fractures in unconventional reservoirs.This work presents an unconventional gas reservoir simulator and its application to quantify hydraulic fractures in shale gas reservoirs using transient pressure data.The numerical model incorporates most known physical processes for gas production from unconventional reservoirs,including two-phase flow of liquid and gas,Klinkenberg effect,non-Darcy flow,and nonlinear adsorption.In addition,the model is able to handle various types and scales of fractures or heterogeneity using continuum,discrete or hybrid modeling approaches under different well production conditions of varying rate or pressure.Our modeling studies indicate that the most sensitive parameter of hydraulic fractures to early transient gas flow through extremely low permeability rock is actually the fracture-matrix contacting area,generated by fracturing stimulation.Based on this observation,it is possible to use transient pressure testing data to estimate the area of fractures generated from fracturing operations.We will conduct a series of modeling studies and present a methodology using typical transient pressure responses,simulated by the numerical model,to estimate fracture areas created or to quantity hydraulic fractures with traditional well testing technology.The type curves of pressure transients from this study can be used to quantify hydraulic fractures in field application.展开更多
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...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.展开更多
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 unstruct...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.展开更多
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....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.展开更多
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 hi...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.展开更多
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 c...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.展开更多
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 va...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.展开更多
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 ...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 11th Postgraduate Technology Innovation Project of North University of China(No.20141142)
文摘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.
基金The 11th Postgraduate Technological Innovation Project of North University of China(No.20141150)
文摘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.
基金financial support from the special fund of China’s central government for the development of local colleges and universities―the project of national first-level discipline in Oil and Gas Engineering, the National Science Fund for Distinguished Young Scholars of China (Grant No. 51125019)the National Program on Key fundamental Research Project (973 Program, Grant No. 2011CB201005)
文摘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.
基金supported by the Fundamental Research Funds for the Central Universities (No.SW-JTU09CX009)
文摘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.
基金Science Foundation of China University of Petroleum, Beijing (No.YJRC-2011-02)for the financial support during this research
文摘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.
文摘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.
基金the financial support provided by the National Natural Science Foundation of China(No.52104043)。
文摘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.
基金This work is supported by the National Natural Science Foundation of China(No.52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(No.BYESS2023262)Science Foundation of China University of Petroleum,Beijing(No.2462022BJRC004).
文摘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.
基金This work was sponsored by the National Natural Science Foun-dation of China(Grant No.52002265)the China Postdoctoral Science Foundation(Grant No.2022M712930).
文摘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.
基金The authors are grateful to the financial support from China Postdoctoral Science Foundation(2022M712645)Opening Fund of Key Laboratory of Enhanced Oil Recovery(Northeast Petroleum University),Ministry of Education(NEPU-EOR-2021-03).
文摘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.
文摘This paper presents the solutions of three -dimensional unsteady flow in Laplace domain for horizontal well in double porous media with vertical feed, using the method both integral transforms and variable separation. The influences of parameters for double porous media ω,λ, attitudes of the horizontal well ZWD. LD and the coefficient of vertical feed a, on the dimensionless pressure behaviors of horizontal well are discussed and calculated by means of the algorithm of Laplace numerical inversion. The diagnosis and specialized graphs are plotted and analysed.
文摘This paper presents the results of an ongoing investigation into transient pressure pulses using Shan- non entropy. Pressure fluctuations (produced by gas-solid two-phase flow during fluidized dense-phase conveying) are recorded by pressure transducers installed at strategic locations along a pipeline. This work validates previous work on identifying the flow mode from pressure signals (Mittal, Mallick, & Wypych, 2014). Two different powders, namely fly ash (median particle diameter 45 μm, particle den- sity 1950 kg/m3. loosely poured bulk density 950 kg/m3) and cement (median particle diameter 15 p,m, particle density 3060 kg/m3, loosely poured bulk density 1070 kg/m3), are conveyed through different pipelines (51 mm I.D. × 70 m length and 63 mm I.D. × 24 m length). The transient nature of pressure fluc- tuations (instead of steady-state behavior) is considered in investigating flow characteristics. Shannon entropy is found to increase along straight pipe sections for both solids and both pipelines. However, Shannon entropy decreases after a bend. A comparison of Shannon entropy among different ranges of superficial air velocity reveals that high Shannon entropy corresponds to very low velocities (i.e. 3-5 m/s) and very high velocities (i.e. 11-14 m/s) while low Shannon entropy corresponds to mid-range velocities (i.e. 6-8 m/s).
文摘Hydraulic fracturing combined with horizontal drilling has been the technology that makes it possible to economically produce natural gas from unconventional shale gas or tight gas reservoirs.Hydraulic fracturing operations,in particular,multistage fracturing treatments along with horizontal wells in unconventional formations create complex fracture geometries or networks,which are difficult to characterize.The traditional analysis using a single vertical or horizontal fracture concept may be no longer applicable.Knowledge of these created fracture properties,such as their spatial distribution,extension and fracture areas,is essential information to evaluate stimulation results.However,there are currently few effective approaches available for quantifying hydraulic fractures in unconventional reservoirs.This work presents an unconventional gas reservoir simulator and its application to quantify hydraulic fractures in shale gas reservoirs using transient pressure data.The numerical model incorporates most known physical processes for gas production from unconventional reservoirs,including two-phase flow of liquid and gas,Klinkenberg effect,non-Darcy flow,and nonlinear adsorption.In addition,the model is able to handle various types and scales of fractures or heterogeneity using continuum,discrete or hybrid modeling approaches under different well production conditions of varying rate or pressure.Our modeling studies indicate that the most sensitive parameter of hydraulic fractures to early transient gas flow through extremely low permeability rock is actually the fracture-matrix contacting area,generated by fracturing stimulation.Based on this observation,it is possible to use transient pressure testing data to estimate the area of fractures generated from fracturing operations.We will conduct a series of modeling studies and present a methodology using typical transient pressure responses,simulated by the numerical model,to estimate fracture areas created or to quantity hydraulic fractures with traditional well testing technology.The type curves of pressure transients from this study can be used to quantify hydraulic fractures in field application.
文摘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.
基金Project (2016YFB1200602-11) supported by National Key R&D Plan of China
文摘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.
文摘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.
基金the financial support received from the College of Petroleum Engineering and Geosciences at KFUPM through the project SF20006 toward the completion of this work。
文摘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.
文摘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.
文摘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.
文摘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.