Analytic Wellbore Temperature Model for Transient Gas-Well Testing

2005 ◽  
Vol 8 (03) ◽  
pp. 240-247 ◽  
Author(s):  
A. Rashid Hasan ◽  
C. Shah Kabir ◽  
Dongqing Lin

Summary Questions arise whether bottomhole pressures (BHPs), derived from their wellhead counterpart (WHP), lend themselves to transient analysis. That is because considerable heat exchange may affect the wellbore-density profile, thereby making the WHP translation a nontrivial exercise. In other words, gas density is dependent on both spatial locations in the wellbore and time during transient testing. Fully coupled wellbore/reservoir simulators are available to tackle this situation. However, they are not readily adaptable for their numeric formulations. This paper presents analytic expressions, derived from first principles, for computing time-dependent fluid temperature at any point in the wellbore during both drawdown and buildup testing. The simplicity of the analytic expressions for Tf (z, t) is profound in that one can compute flowing or shut-in BHPs on a spreadsheet. Two tests were considered to verify the new analytic formulae. In one case, measurements were available at both sandface and surface, and partial wellhead information was available in the other case. We explored a parametric study to assess whether a given wellbore/reservoir system will lend itself to wellhead measurements for valid transient analysis. Reservoir flow capacity (kh) turned out to be the most influential variable. Introduction Gas-well testing is sometimes conducted by measuring pressures at the wellhead. Both cost and circumstance (high pressure/high temperature, or HP/HT)often necessitate WHP monitoring or running the risk of having no tests at all. Methods for computing BHP from wellhead pressures for steady flow in gas wells are well established in the literature. For dry-gas wells, the widely used method of Cullender and Smith is most accurate, as confirmed by subsequent studies. For wet gas, either a two-phase model, such as the one offered by Govier and Fogarasi, or the modified Cullender-Smith approach appears satisfactory. However, these methods apply to steady-state gas flow and implicitly presuppose that the wellbore is in thermal equilibrium with the formation. These assumptions may be tested during a transient test. That is because unsteady-state wellbore heat transfer occurs even after the cessation of the wellbore-fluid-storage period. Steady-state fluid flow ordinarily implies the absence of wellbore effects from the viewpoint of transient testing. Consequently, one needs to develop working equations by conserving mass, momentum, and energy in the wellbore to capture physical phenomena. Earlier, we presented a forward model and showed its capability to reproduce BHP, WHP, and wellhead temperature (WHT) given reservoir and wellbore parameters. However, translation of WHP to BHP was not demonstrated clearly. The intent of this work is to present a framework for rigorous computation of BHP from WHP. To achieve this objective, we developed analytic expressions for depth- and time-dependent fluid temperature during both flow and shut-in tests. These temperature relations, in turn, allow computation of gas density and, therefore, pressure at any point in the wellbore.

2018 ◽  
Vol 6 (1) ◽  
pp. SB1-SB9 ◽  
Author(s):  
Marija Macenić ◽  
Tomislav Kurevija

The aim of our research is to use abandoned deep-hydrocarbon reservoirs and dry wells in the Croatian part of the Pannonian Basin as a geothermal energy source. Croatia has been exploring and exploiting hydrocarbon reserves in the Pannonian Basin from the mid-20th century. Therefore, many oil and gas wells are reaching the end of their production phase and many are already abandoned. These wells could be considered for geothermal energy production through the coaxial heat exchanger principle, which is usually used in shallow geothermal energy extraction. Using the abandoned deep well Pčelić located in the Drava subbasin as a case study, we have derived the available energy and fluid temperature changes during 20 years of operation for two cases: one with a constant base heat load throughout the year and the second as a variable heat load depending on outside air temperatures. We determined that the maximum potential heat extraction in a variable system is 1750 MWh per year, with 1.5 MW of peak heating power in winter, depending on the sink temperature, climate, and consumer input data. The maximum theoretical constant heat extraction for possible industrial direct heating could be 400 kW during the entire period of 20 years, with fluid temperature reaching steady state at a favorable 50°C. To define steady-state ratio between extracted heat and consumed energy at the circulating pump, we evaluated seasonal performance factor (SPF) analysis similar to heat pump systems. Lower values of SPF linked to a higher flow rate implies higher energy extraction. Our results show that when using a lower flow, steady-state SPF ratio is as high as 280, and for a higher flow, steady-state SPF ratio drops to only 25.


Author(s):  
P Yang ◽  
J Cui ◽  
J M Jin ◽  
D Dowson

Time-dependent thermal and non-Newtonian elastohydrodynamic lubrication of an elliptical point contact subjected to a normal harmonic vibration was studied numerically in this work. The contact was idealized as between an infinite plane and a spherical roller. The normal vibration of the roller was described by specifying the centre of the spherical roller to the infinite plane (without deformation) as a cyclic function of time. The shear-thinning rheological property of the lubricant was described by the Eyring model. The time-dependent numerical solution was achieved instant after instant in each period of a vibration. The periodic errors were checked at the end of each vibration cycle until the responses of variables such as pressure, film thickness, and temperature were all periodic functions with the same frequency of the roller's vibration. At each instant, the pressure field was solved with a multi-grid method, the surface deflection produced by pressure was determined with a multi-level multi-integration technique, the non-Newtonian flow of the lubricant was considered by using the equivalent viscosity calculated according to the shear-strain rate along the entrainment direction only, and the temperature field was evaluated with a finite-difference scheme through a column-by-column relaxation process. The computing time for a cyclic solution was 12–15 h on a personal computer with a 3.0 GHz central processing unit. The effects of both the amplitude and the frequency of the vibration were investigated. It was shown that the time-dependent solution is significantly different from the steady-state solution, especially when the amplitude of vibration is large and the frequency of vibration is high. Corresponding to a typical thermal and non-Newtonian case, numerical solutions were also obtained under isothermal and Newtonian, isothermal and non-Newtonian, and thermal and Newtonian conditions. Comparisons between these solutions indicate that, under time-dependent conditions, the effects of thermal and non-Newtonian flow are similar to those under steady-state conditions.


2009 ◽  
Vol 1158 ◽  
Author(s):  
Bang-Hung Tsao ◽  
Katie Sondergelt ◽  
Jacob Lawson ◽  
James Scofield ◽  
Levi Elston

AbstractA three dimensional thermal-fluid and stress model of a single chip SiC power sub-module was generated using ANSYS in order to determine the maximum temperature and deformation under various conditions. The effects of heat flux, working fluid temperature and differential pressure on temperature and thermal stress contours were of particular concern. Steady state analysis with water as the working fluid, a simulated heat flux of 11.12×104 W/m2, an interface coupling film coefficient of either 30 or 200 W/m2-K between the cooling plate and fluid, and ambient film coefficients from 6 W/m2-K to 300 W/m2-K, predicts maximum device junction temperatures between 374 and 316 K, and corresponding deformations from .0351% to .0293%. Under the same boundary and loading conditions, but with air as the working fluid, the deformations reached .0405% to .0296%, with temperatures between 427 and 316 K. Transient analysis also showed junction temperatures in the predicted range and determined the time to reach steady state to be between 150 and 2500 seconds depending on the boundary conditions. Experiments were conducted in order to validate ANSYS results.


Author(s):  
Veikko Taivassalo ◽  
Sirpa Kallio ◽  
Juho Peltola

AbstractCFD simulations of single-phase flows are regularly performed as steady-state utilizing closure models of varying complexity. On the contrary, dense gas-solid flows are usually computed as time dependent. These simulations commonly require a small time step and a fine mesh resulting in costly and time-consuming computations. In case of large industrial circulating fluidized beds (CFB), the steady-state CFD modeling would be an attractive alternative for the transient simulations, if reliable closure models for the time-averaged transport equations were available. The multiphase closure models developed for time-dependent CFB computations are not as such applicable to the steady-state approach. For instance, the fraction of the momentum transfer expressed by the velocities is significantly smaller in the steady-state models than in the transient ones. Therefore, the steady-state simulations rely more on the closure relations and especially on the models for inter-phase momentum transfer and for the Reynolds stress terms.Several attempts to develop closure models for coarsemesh and steady-state simulations have been presented in the literature. In this paper, a novel steady-state simulation approach for a CFB process and a corresponding CFD model are introduced. A successful steady-state simulation for a test case is presented. Compared to the time-dependent simulations, the computing time is reduced by a factor of an order of 1000.


RSC Advances ◽  
2021 ◽  
Author(s):  
Guanzhao Wen ◽  
Xianshao Zou ◽  
Rong Hu ◽  
Jun Peng ◽  
Zhifeng Chen ◽  
...  

Ground- and excited-states properties of N2200 have been studied by steady-state and time-resolved spectroscopies as well as time-dependent density functional theory calculations.


1988 ◽  
Vol 8 (5) ◽  
pp. 1957-1969 ◽  
Author(s):  
R A Shapiro ◽  
D Herrick ◽  
R E Manrow ◽  
D Blinder ◽  
A Jacobson

As an approach to understanding the structures and mechanisms which determine mRNA decay rates, we have cloned and begun to characterize cDNAs which encode mRNAs representative of the stability extremes in the poly(A)+ RNA population of Dictyostelium discoideum amoebae. The cDNA clones were identified in a screening procedure which was based on the occurrence of poly(A) shortening during mRNA aging. mRNA half-lives were determined by hybridization of poly(A)+ RNA, isolated from cells labeled in a 32PO4 pulse-chase, to dots of excess cloned DNA. Individual mRNAs decayed with unique first-order decay rates ranging from 0.9 to 9.6 h, indicating that the complex decay kinetics of total poly(A)+ RNA in D. discoideum amoebae reflect the sum of the decay rates of individual mRNAs. Using specific probes derived from these cDNA clones, we have compared the sizes, extents of ribosome loading, and poly(A) tail lengths of stable, moderately stable, and unstable mRNAs. We found (i) no correlation between mRNA size and decay rate; (ii) no significant difference in the number of ribosomes per unit length of stable versus unstable mRNAs, and (iii) a general inverse relationship between mRNA decay rates and poly(A) tail lengths. Collectively, these observations indicate that mRNA decay in D. discoideum amoebae cannot be explained in terms of random nucleolytic events. The possibility that specific 3'-structural determinants can confer mRNA instability is suggested by a comparison of the labeling and turnover kinetics of different actin mRNAs. A correlation was observed between the steady-state percentage of a given mRNA found in polysomes and its degree of instability; i.e., unstable mRNAs were more efficiently recruited into polysomes than stable mRNAs. Since stable mRNAs are, on average, "older" than unstable mRNAs, this correlation may reflect a translational role for mRNA modifications that change in a time-dependent manner. Our previous studies have demonstrated both a time-dependent shortening and a possible translational role for the 3' poly(A) tracts of mRNA. We suggest, therefore, that the observed differences in the translational efficiency of stable and unstable mRNAs may, in part, be attributable to differences in steady-state poly(A) tail lengths.


1985 ◽  
Vol 50 ◽  
Author(s):  
J-E. Andersson ◽  
O. Persson

AbstractThe results from a large number of single-hole packer tests in crystalline rock from three test sites in Sweden have been analysed statistically. Average hydraulic conductivity values for 25 m long test intervals along boreholes with a maximal length of about 700 m are used in this study. A comparison between steady state and transient analysis of the same test data has been performed.The mean vaule of the hydraulic conductivity determined from steady state analysis was found to be about two to three times higher compared to transient analysis. However, in some cases the steady state analysis resulted in 10 to 20 times higher values compared to the transient analysis. Such divergence between the two analysis methods may be caused by deviations from the assumed flow pattern, borehole skin effects and influence of hydraulic boundaries.


2021 ◽  
Author(s):  
Mauricio Espinosa ◽  
Jairo Leal ◽  
Ron Zbitowsky ◽  
Eduardo Pacheco

Abstract This paper highlights the first successful application of a field deployment of a high-temperature (HT) downhole shut-in tool (DHSIT) in multistage fracturing completions (MSF) producing retrograde gas condensate and from sour carbonate reservoirs. Many gas operators and service providers have made various attempts in the past to evaluate the long-term benefit of MSF completions while deploying DHSIT devices but have achieved only limited success (Ref. 1 and 2). During such deployments, many challenges and difficulties were faced in the attempt to deploy and retrieve those tools as well as to complete sound data interpretation to successfully identify both reservoir, stimulation, and downhole productivity parameters, and especially when having a combination of both heterogeneous rocks having retrograde gas pressure-volume-temperature (PVT) complexities. Therefore, a robust design of a DHSIT was needed to accurately shut-in the well, hold differential pressure, capture downhole pressure transient data, and thereby identify acid fracture design/conductivity, evaluate total KH, reduce wellbore storage effects, properly evaluate transient pressure effects, and then obtain a better understanding of frac geometry, reservoir parameters, and geologic uncertainties. Several aspects were taken into consideration for overcoming those challenges when preparing the DHSIT tool design including but not limited to proper metallurgy selection, enough gas flow area, impact on well drawdown, tool differential pressure, proper elastomer selection, shut-in time programming, internal completion diameter, and battery operation life and temperature. This paper is based on the first successful deployment and retrieval of the DHSIT in a 4-½" MSF sour carbonate gas well. The trial proved that all design considerations were important and took into consideration all well parameters. This project confirmed that DHSIT devices can successfully withstand the challenges of operating in sour carbonate MSF gas wells as well as minimize operational risk. This successful trial demonstrates the value of utilizing the DHSIT, and confirms more tangible values for wellbore conductivity post stimulation. All this was achieved by the proper metallurgy selection, maximizing gas flow area, minimizing the impact on well drawdown, and reducing well shut-in time and deferred gas production. Proper battery selection and elastomer design also enabled the tool to be operated at temperatures as high as 350 °F. The case study includes the detailed analysis of deployment and retrieval lessons learned, and includes equalization procedures, which added to the complexity of the operation. The paper captures all engineering concepts, tool design, setting packer mechanism, deployment procedures, and tool equalization and retrieval along with data evaluation and interpretation. In addition to lessons learned based on the field trial, various recommendations will be presented to minimize operational risk, optimize shut-in time and maximize data quality and interpretation. Utilizing the lessons learned and the developed procedures presented in this paper will allow for the expansion of this technology to different gas well types and formations as well as standardize use to proper evaluate the value of future MSF completions and stimulation designs.


Author(s):  
R.A. Gasumov ◽  
◽  
E.R. Gasumov ◽  

The article discusses the modes of movement of gas-liquid flows in relation to the operating conditions of waterlogged gas wells at a late stage of field development. Algorithms have been developed for calculating gas well operation modes based on experimental work under conditions that reproduce the actual operating conditions of flooded wells of Cenomanian gas deposits. The concept of calculating the technological mode of operation of gas wells with a single-row elevator according to the critical velocity of the upward flow is considered based on the study of the equilibrium conditions of two oppositely directed forces: the gravity of water drops directed downward and the lifting force moving water drops with a gas flow directed upward. A calculation was made according to the method of the averaged physical parameters of formation water and natural gas in the conditions of flooded Cenomanian gas wells in Western Siberia. The results of a study of the dependence of the critical flow rate of Cenomanian wells on bottomhole pressure and diameter of elevator pipes are presented.


Sign in / Sign up

Export Citation Format

Share Document