Estimation of Solubility and Diffusivity of Gases in Heavy Oils by Use of Late-Time Pressure-Decay Data: An Analytical/Graphical Approach

SPE Journal ◽  
2015 ◽  
Vol 20 (04) ◽  
pp. 717-728 ◽  
Author(s):  
Francisco J. Pacheco-Roman ◽  
S. Hossein Hejazi

Summary Solubility and diffusivity of gases in heavy oils, quantified by Henry's constant (Hij) and the diffusion coefficient (D), respectively, are essential properties for the design of recovery processes that require the injection of gas or vapor solvents into the reservoir. Data, obtained from various experimental procedures such as the pressure-decay technique (PDT), are used to estimate these two parameters. The PDT uses a pressure/volume/temperature (PVT) cell where the gas-phase pressure declines as gas diffuses into heavy oil following an early- and a late-time regime. Current approaches to analyze data from the conventional PDT are either graphical techniques based on early-time data or full numerical simulation. Early-time data, when the diffusing gas has not reached the bottom of the PVT cell, do not provide enough information to simultaneously estimate both the diffusion coefficient and Henry's constant. Hence, existing graphical procedures are limited to diffusion-coefficient estimation. In this paper, we propose a novel graphical technique to estimate the diffusion coefficient and Henry's constant by use of the late-time data from pressure-decay experiments. Our method is derived from the modeling of gas-phase pressure decay by use of Fick's second law and gas-phase mass-balance equations. We use the integral method to provide an approximate analytical solution to the set of equations. In addition, by use of the resultant solution, we develop a simple graphical method to directly estimate both the diffusion coefficient and Henry's constant. The estimated parameters through the proposed technique for methane/bitumen and carbon dioxide/bitumen experiments are in close agreement with those reported in the literature.

1999 ◽  
Vol 2 (03) ◽  
pp. 271-280 ◽  
Author(s):  
Ekrem Kasap ◽  
Kun Huang ◽  
Than Shwe ◽  
Dan Georgi

Summary The formation-rate-analysis (FRASM) technique is introduced. The technique is based on the calculated formation rate by correcting the piston rate with fluid compressibility. A geometric factor is used to account for irregular flow geometry caused by probe drawdown. The technique focuses on the flow from formation, is applicable to both drawdown and buildup data simultaneously, does not require long buildup periods, and can be implemented with a multilinear regression, from which near-wellbore permeability, p * and formation fluid compressibility are readily determined. The field data applications indicate that FRA is much less amenable to data quality because it utilizes the entire data set. Introduction A wireline formation test (WFT) is initiated when a probe from the tool is set against the formation. A measured volume of fluid is then withdrawn from the formation through the probe. The test continues with a buildup period until pressure in the tool reaches formation pressure. WFTs provide formation fluid samples and produce high-precision vertical pressure profiles, which, in turn, can be used to identify formation fluid types and locate fluid contacts. Wireline formation testing is much faster compared with the regular pressure transient testing. Total drawdown time for a formation test is just a few seconds and buildup times vary from less than a second (for permeability of hundreds of millidarcy) to half a minute (for permeability of less than 0.1 md), depending on system volume, drawdown rate, and formation permeability. Because WFT tested volume can be small (a few cubic centimeters), the details of reservoir heterogeneity on a fine scale are given with better spatial resolution than is possible with conventional pressure transient tests. Furthermore, WFTs may be preferable to laboratory core permeability measurements since WFTs are conducted at in-situ reservoir stress and temperature. Various conventional analysis techniques are used in the industry. Spherical-flow analysis utilizes early-time buildup data and usually gives permeability that is within an order of magnitude of the true permeability. For p* determination, cylindrical-flow analysis is preferred because it focuses on late-time buildup data. However, both the cylindrical- and spherical-flow analyses have their drawbacks. Early-time data in spherical-flow analysis results in erroneous p* estimation. Late-time data are obtained after long testing times, especially in low-permeability formations; however, long testing periods are not desirable because of potential tool "sticking" problems. Even after extended testing times, the cylindrical-flow period may not occur or may not be detectable on WFTs. When it does occur, permeability estimates derived from the cylindrical-flow period may be incorrect and their validity is difficult to judge. New concepts and analysis techniques, combined with 3-D numerical studies, have recently been reported in the literature.1–7 Three-dimensional numerical simulation studies1–6 have contributed to the diagnosis of WFT-related problems and the improved analysis of WFT data. The experimental studies7 showed that the geometric factor concept is valid for unsteady state probe pressure tests. This study presents the FRA technique8 that can be applied to the entire WFT where a plot for both drawdown and buildup periods renders straight lines with identical slopes. Numerical simulation studies were used to generate data to test both the conventional and the FRA techniques. The numerical simulation data are ideally suited for such studies because the correct answer is known (e.g., the input data). The new technique and the conventional analysis techniques are also applied to the field data and the results are compared. We first review the theory of conventional analysis techniques, then present the FRA technique for combined drawdown and buildup data. A discussion of the numerical results and the field data applications are followed by the conclusions. Analysis Techniques It has been industry practice to use three conventional techniques, i.e., pseudo-steady-state drawdown (PSSDD), spherical and cylindrical-flow analyses, to calculate permeability and p* Conventional Techniques Pseudo-Steady-State Drawdown (PSSDD). When drawdown data are analyzed, it is assumed that late in the drawdown period the pressure drop stabilizes and the system approaches to a pseudo-steady state when the formation flow rate is equal to the drawdown rate. PSSDD permeability is calculated from Darcy's equation with the stabilized (maximum) pressure drop and the flowrate resulting from the piston withdrawal:9–11 $$k {d}=1754.5\left({q\mu \over r {i}\Delta p {{\rm max}}}\right),\eqno ({\rm 1})$$where kd=PSSDD permeability, md. The other parameters are given in Nomenclature.


2022 ◽  
Author(s):  
Mark Mcclure ◽  
Garrett Fowler ◽  
Matteo Picone

Abstract In URTeC-123-2019, a group of operators and service companies presented a step-by-step procedure for interpretation of diagnostic fracture injection tests (DFITs). The procedure has now been applied on a wide variety of data across North and South America. This paper statistically summarizes results from 62 of these DFITs, contributed by ten operators spanning nine different shale plays. URTeC-123-2019 made several novel claims, which are tested and validated in this paper. We find that: (1) a ‘compliance method’ closure signature is apparent in the significant majority of DFITs; (2) in horizontal wells, early time pressure drop due to near-wellbore/midfield tortuosity is substantial and varies greatly, from 500 to 6000+ psi; (3) in vertical wells, early-time pressure drop is far weaker; this supports the interpretation that early- time pressure drop in horizontal wells is caused by near-wellbore/midfield tortuosity from transverse fracture propagation; (4) the (not recommended) tangent method of estimating closure yields Shmin estimates that are 100-1000+ psi lower than the estimate from the (recommended) compliance method; the implied net pressure values are 2.5x higher on average and up to 5-6x higher; (5) as predicted by theory, the difference between the tangent and compliance stress and net pressure estimates increases in formations with greater difference between Shmin and pore pressure; (6) the h-function and G-function methods allow permeability to be estimated from truncated data that never reaches late-time impulse flow; comparison shows that they give results that are close to the permeability estimates from impulse linear flow; (7) false radial flow signatures occur in the significant majority of gas shale DFITs, and are rare in oil shale DFITs; (8) if false radial signatures are used to estimate permeability, they tend to overestimate permeability, often by 100x or more; (9) the holistic-method permeability correlation overestimates permeability by 10-1000x; (10) in tests that do not reach late-time impulse transients, it is reasonable to make an approximate pore pressure estimate by extrapolating the pressure from the peak in t*dP/dt using a scaling of t^(-1/2) in oil shales and t^(3/4) in gas shales. The findings have direct practical implications for operators. Accurate permeability estimates are needed for calculating effective fracture length and for optimizing well spacing and frac design. Accurate stress estimation is fundamental to hydraulic fracture design and other geomechanics applications.


2011 ◽  
Vol 305 (2) ◽  
pp. 132-144 ◽  
Author(s):  
Seyyed M. Ghaderi ◽  
S. Hamed Tabatabaie ◽  
Hassan Hassanzadeh ◽  
Mehran Pooladi-Darvish

1976 ◽  
Vol 16 (1) ◽  
pp. 99 ◽  
Author(s):  
A.K. Khurana

Bottom-hole pressure tests conducted in the Kingfish oil reservoir (located in Gippsland Basin - Offshore Victoria) during 1974 and 1975 using a high sensitivity surface recording electronic bottom-hole pressure gauge indicated the presence of sinusoidal pressure oscillations in the reservor. The oscillations are of the order of 0.1 psi in amplitude and their frequency suggests that they are in some way related to tidal phenomena.Although the oscillations do not affect production, they do influence interpretation of pressure build-up and pulse tests. Interpretations of both late time pressure build-up behaviour and pulse tests of small response magnitude and long time lags are considered to be particularly susceptible to errors due to these oscillations if they are not recognized and corrected for. Interpretations of early time pressure build-up data and pulse tests of definite response and relatively short time lags are not regarded as being significantly affected.The physical mechanism causing these pressure oscillations in the reservoirs is not known. However, one of the various possible hypotheses is that the Latrobe Formation sands could be outcropping on the ocean floor at abyssal depths southeast of Kingfish and that the pressure transients generated by changes in the hydrostatic head due to surfate tides are transmitted hydraulically to the reservoir. If this hypothesis is proved to be valid it could influence pressure performance predictions of Gippsland Basin reservoirs.


SPE Journal ◽  
2013 ◽  
Vol 18 (05) ◽  
pp. 982-992 ◽  
Author(s):  
Michael Prats ◽  
R.. Raghavan

Summary The instantaneous source solutions of Prats and Raghavan (2012) and the method of images are used to develop analytic expressions for the pressure distribution in a three-region composite reservoir of finite thickness produced by a finite-length horizontal well that is oriented perpindicular to the interfaces. The composite reservoir is assumed to be infinite in its lateral extents; the outer regions represent the reservoir, and the central region represents a thin natural fracture of relatively high permeability. In most of the cases considered, the well is completed in all three zones. The computational scheme is shown to be both viable and robust. The Shanks (1955) transformation is used to accelerate convergence. Pressure traces are logarithmic in time at early and late times for any well configuration examined here. Early-time pressure characteristics are similar to those discussed in Prats and Raghavan (2012). The duration of the early semilogarithmic responses is mediated not only by the presence of the higher-permeability natural fracture, as before, but also now by the interaction of the upper and lower boundaries with the well. Late-time semilogarithmic responses, however, are distinctly different. Their slope is inversely proportional to the product of the formation thickness and the arithmetic average permeability of the two regions that sandwich the fracture. This result holds even when the well does not cross the natural fracture. We expect this conclusion to apply to a composite system consisting of more than three regions. Observations concerning the late-time slope represent the central finding of this study. A relationship is given for the late-time performance of any horizontal well in terms of that of a vertical well with a constant pseudoskin. Pseudoskin equivalents are reported for all cases discussed.


2019 ◽  
Vol 24 (4) ◽  
pp. 653-663 ◽  
Author(s):  
Xin Wu ◽  
Guangyou Fang ◽  
Guoqiang Xue ◽  
Lihua Liu ◽  
Leisong Liu ◽  
...  

Over the past decade, helicopter-borne transient electromagnetic (HTEM) systems have been rapidly developed. A new HTEM prototype (referred to as a CAS-HTEM) has been developed by the Chinese Academy of Sciences. In terms of hardware, the CAS-HTEM system uses an inflatable structure to carry the transmitting loop, which significantly reduces the weight of the system and makes it easier to transport. A dual gain receiver was innovated to extend the dynamic range of the system. In addition, an observation circuit for transmitting voltage waveform is introduced, so that the derivative waveform of transmitting current with higher SNR could be calculated. In terms of data processing, more reliable early time data could be obtained by band-limited effect removal; a weighted stacking algorithm is introduced to reduce the narrow band noise more effectively and increase the sensitivity of data to the anomaly location; a method based on τ-domain transform is used for late time signal processing. The results of the field test which was carried out in Inner Mongolia were found to be consistent with the drill data, which effectively verified the performance of this HTEM prototype.


SPE Journal ◽  
2016 ◽  
Vol 21 (01) ◽  
pp. 131-143 ◽  
Author(s):  
Huazhou Li ◽  
Daoyong Yang

Summary A novel methodology was developed to determine the molecular-diffusion coefficient for each component of the solvent/CO2 mixture in heavy oil under reservoir conditions on the basis of the pressure-decay theory. Experimentally, molecular-diffusion tests for the solvent/CO2/heavy-oil systems (i.e., pure-CO2/heavy-oil system, C3H8/CO2/heavy-oil system, and n-C4H10/CO2/heavy-oil system) are performed with a DBR pressure/volume/temperature system at constant temperature and decayed pressure. Theoretically, the Peng-Robinson equation of state combined with a 1D diffusion model is developed to describe the diffusion process of solvent/CO2 mixture in heavy oil. The composition analysis in the beginning and the end of pressure-decay experiments for the solvent/CO2/heavy-oil system indicate that the gas-phase solvent fraction decreases as diffusion proceeds, whereas the gas-phase CO2 fraction increases during the tests. One can determine the individual molecular-diffusion coefficient for each component in the mixture by minimizing the discrepancy between the measured composition change and the calculated composition change with the diffusion model. The newly developed methodology is successfully validated with the diffusion tests on the two solvent/CO2 mixtures: C3H8/CO2/heavy-oil system and n-C4H10/CO2/heavy-oil system. As for the solvent/CO2 mixtures tested, the molecular-diffusion coefficient of solvent in heavy oil is found to be significantly larger than that of CO2 in heavy oil. At similar test conditions, the C3H8/CO2/heavy-oil system ends up with a swelling factor of 1.058 after 168 hours of diffusion test, in comparison to 1.031 for the CO2/heavy-oil system.


Geofluids ◽  
2019 ◽  
Vol 2019 ◽  
pp. 1-9
Author(s):  
Yu Zhao ◽  
Chaolin Wang ◽  
Yongfa Zhang ◽  
Qiang Liu

The pressure-pulse decay is a preferred technique for determining permeability of unconventional gas reservoir rocks. The pressure-pulse decay often shows quite different characteristics during the early time and the later time. Most approaches for estimating the permeability proposed in the literature are required to use the later-time pressure-pulse decay measurements. However, the later-time data are often selected subjectively, lacking a universal criterion. In this paper, a method of differentiating the early-time and late-time behavior for pressure-pulse decay test is proposed. The analytical results show that the critical time (dimensionless time) of early-/late-time decay characteristics mainly depends on the volume ratios, and it increases first and then decreases with the volume ratios. The critical time for cases with same chamber sizes is much less than that for cases with unequal chamber sizes. Applicability of the proposed methods is examined using a numerical simulator, TOUGH+REALGASBRINE. The numerical results show that the pressure gradient along the sample varies nonlinearly at the early time and becomes a constant at the late time. Then, the proposed method is applied to real data for permeability estimation. It is found that the early-time behavior is negligible as the volume ratio takes on small values. As the volume ratios increase, the deviation becomes significant and considerable permeability errors will be produced if these early-time data are used.


2019 ◽  
Vol 488 (3) ◽  
pp. 3783-3793 ◽  
Author(s):  
Kornpob Bhirombhakdi ◽  
Ryan Chornock ◽  
Adam A Miller ◽  
Alexei V Filippenko ◽  
S Bradley Cenko ◽  
...  

ABSTRACT SN 2008es is one of the rare cases of a Type II superluminous supernova (SLSN), showing no narrow features in its early-time spectra, and therefore, its powering mechanism is under debate between circumstellar interaction (CSI) and magnetar spin-down. Late-time data are required for better constraints. We present optical and near-infrared (NIR) photometry obtained from Gemini, Keck, and Palomar Observatories from 192 to 554 d after explosion. Only broad H α emission is detected in a Gemini spectrum at 288 d. The line profile exhibits red-wing attenuation relative to the early-time spectrum. In addition to the cooling SN photosphere, an NIR excess with blackbody temperature ∼1500 K and radius ∼1016 cm is observed. This evidence supports dust condensation in the cool dense shell being responsible for the spectral evolution and NIR excess. We favour CSI, with ∼2–3 M⊙ of circumstellar material (CSM) and ∼10–20 M⊙ of ejecta, as the powering mechanism, which still dominates at our late-time epochs. Both models of uniform density and steady wind fit the data equally well, with an effective CSM radius ∼1015 cm, supporting the efficient conversion of shock energy to radiation by CSI. A low amount (≲0.4 M⊙) of 56Ni is possible but cannot be verified yet, since the light curve is dominated by CSI. The magnetar spin-down powering mechanism cannot be ruled out, but is less favoured because it overpredicts the late-time fluxes and may be inconsistent with the presence of dust.


Sign in / Sign up

Export Citation Format

Share Document