Thesis (Ph. D.)--University of Wyoming, 2006. / Title from PDF title page (viewed on Dec. 21, 2007). Includes bibliographical references (p. 199-216).

<p>Høye poretrykk utgjør en risiko for boreoperasjoner på Norsk Sokkel og internasjonalt. Denne risikoen kan reduseres dersom man har kjennskap til poretrykksforholdene før boring. Poretrykk er også en viktig parameter for felt i produksjon, og kunnskap om hvordan dette utvikler seg over tid vil kunne ha stor betydning for økt oljeutvinning. Seismiske data inneholder informasjon om poretrykket og kan derfor bidra til økt kunnskap på dette området. Avhandlingen tar for seg hvordan hastighets- og amplitudeinformasjon fra seismiske data kan brukes for å estimere poretrykk.</p> / <p>Abnormally high pore pressures in the subsurface pose a hazard to drilling operations worldwide. The problem is not unusual on the Norwegian Continental Shelf. Knowledge of the pore pressure prior to drilling may reduce the risk related to drilling in high pressure zones. Pore pressure is also a vital paramter for producinig fields, and knowledge of how the pressure develops over time can be important for increased oil recovery. Seismic data contain information on the pore pressure and may contribute to increased understanding of subsurface pressure conditions. The thesis deals with methods for estimation of pressure from seismic velocity and amplitude data.</p>

Ranking of multiple realizations is an important step when the processing time for a realization is large. This is the case in reservoir flow simulation and in other areas of geology, environmental and even medical applications. Significant uncertainty exists in all reservoirs especially at unsampled locations where the geological heterogeneity and connectivity are impossible to exactly predict between wells. Geostatistical techniques are used to construct models of static properties such as lithofacies, porosity, permeability and residual fluid saturations and provide multiple equally probable realizations of these properties. The number of realizations that is required for modeling the uncertainty may be large; usually 100 realizations are considered enough to quantify uncertainty. However, this number of realizations is still too high for processing by a flow simulator. This thesis aims at developing a robust and reliable ranking methodology to rank the realizations using a static ranking measure. The outcome is the identification of the high, low, and intermediate ranking realizations for further detailed simulations. The methodology was developed for the steam assisted gravity drainage (SAGD) reservoir application. This thesis will consider the cumulative oil produced (COPrate) and cumulative steam-oil-ratio (CSOR) as the ranking parameters in the flow simulations, hereafter called performance parameter. Connected hydrocarbon volume (CHV) was the parameter that was used in the ranking methodology as the static ranking measure. High calibration between the performance parameters and the CHV would indicate the success of the proposed ranking methodology. The ranking methodology was validated against the results of the flow simulations. The results indicate a mediocre correlation between the SAGD performance parameters and CHV. The ranking methodology was modified by incorporating the average reservoir permeability. Significant improvement in the correlation between the static ranking measure and the SAGD performance parameters resulted. / Petroleum Engineering

Thesis (M.Eng.)--Memorial University of Newfoundland, 2005. / Bibliography: leaves 130-135.

he feasibility of sequestering supercritical CO2 in depleted gas reservoirs. The experimental runs involved the following steps. First, the 1 ft long by 1 in. diameter carbonate core is inserted into a viton Hassler sleeve and placed inside an aluminum coreholder that is then evacuated. Second, with or without connate water, the carbonate core is saturated with methane. Third, supercritical CO2 is injected into the core with 300 psi overburden pressure. From the volume and composition of the produced gas measured by a wet test meter and a gas chromatograph, the recovery of methane at CO2 breakthrough is determined. The core is scanned three times during an experimental run to determine core porosity and fluid saturation profile: at start of the run, at CO2 breakthrough, and at the end of the run. Runs were made with various temperatures, 20°C (68°F) to 80°C (176°F), while the cell pressure is varied, from 500 psig (3.55 MPa) to 3000 psig (20.79 MPa) for each temperature. An analytical study of the experimental results has been also conducted to determine the dispersion coefficient of CO2 using the convection-dispersion equation. The dispersion coefficient of CO2 in methane is found to be relatively low, 0.01-0.3 cm2/min.. Based on experimental and analytical results, a 3D simulation model of one eighth of a 5-spot pattern was constructed to evaluate injection of supercritical CO2 under typical field conditions. The depleted gas reservoir is repressurized by CO2 injection from 500 psi to its initial pressure 3,045 psi. Simulation results for 400 bbl/d CO2 injection may be summarized as follows. First, a large amount of CO2 is sequestered: (i) about 1.2 million tons in 29 years (0 % initial water saturation) to 0.78 million tons in 19 years (35 % initial water saturation) for 40-acre pattern, (ii) about 4.8 million tons in 112 years (0 % initial water saturation) to 3.1 million tons in 73 years (35 % initial water saturation) for 80-acre pattern. Second, a significant amount of natural gas is also produced: (i) about 1.2 BSCF or 74 % remaining GIP (0 % initial water saturation) to 0.78 BSCF or 66 % remaining GIP (35 % initial water saturation) for 40-acre pattern, (ii) about 4.5 BSCF or 64 % remaining GIP (0 % initial water saturation) to 2.97 BSCF or 62 % remaining GIP (35 % initial water saturation) for 80-acre pattern. This produced gas revenue could help defray the cost of CO2 sequestration. In short, CO2 sequestration in depleted gas reservoirs appears to be a win-win technology.

Spontaneous imbibition plays a very important role in the displacement mechanism of non-wetting fluid in naturally fractured reservoirs. We developed a new 2D two-phase finite element numerical model, as available commercial simulators cannot be used to model small-scale experiments with different boundary conditions as well as complex boundary conditions such as fractures and vugs. Starting with the basic equation of fluid flow, we derived the non-linear diffusion saturation equation. This equation cannot be put in weighted-integral weak variational form and hence Rayleigh-Ritz finite element method (FEM) cannot be applied. Traditionally, the way around it is to use higher order interpolation functions and use Galerkin FEM or reduce the differentiability requirement and use Mixed FEM formulation. Other FEM methods can also be used, but iterative nature of those methods makes them unsuitable for solving large-scale field problems. But if we truncate the non-linear terms and decouple the dependent variables, from the spatial as well as the temporal domains of the primary variable to solve them analytically, the non-linear FEM problem reduces to a simple weighted integral form, which can be put into its corresponding weak form. The advantage of using Rayleigh-Ritz method is that it has immediate effect on the computation time required to solve a particular problem apart from incorporating complex boundary conditions. We compared our numerical models with the analytical solution of this diffusion equation. We validated the FDM numerical model using X-Ray Tomography (CT) experimental data from the single-phase spontaneous imbibition experiment, where two simultaneously varying parameters of weight gain and CT water saturation were used and then went ahead and compared the results of FEM model to that of FDM model. A two-phase field size example was taken and results from a commercial simulator were compared to the FEM model to bring out the limitations of this approach.

Fluid flow mechanisms in a large naturally fractured heterogeneous carbonate reservoir were investigated in this manuscript. A very thin layer with high permeability that produces the majority of production from specific wells and is deemed the Super-K Zone was investigated. It is known that these zones are connected to naturally occurring fractures. Fluid flow in naturally fractured reservoirs is a very difficult mechanism to understand. To accomplish this mission, the Super-K Zone and fractures were treated as two systems. Reservoir management practices and decisions should be very carefully reviewed and executed in this dual continuum reservoir based on the results of this work. Studying this dual media flow behavior is vital for better future completion strategies and for enhanced reservoir management decisions. The reservoir geology, Super-K identification and natural fractures literature were reviewed. To understand how fluid flows in such a dual continuum reservoir, a dual permeability simulation model has been studied. Some geological and production iv data were used; however, due to unavailability of some critical values of the natural fractures, the model was assumed hypothetical. A reasonable history match was achieved and was set as a basis of the reservoir model. Several sensitivity studies were run to understand fluid flow behavior and prediction runs were executed to help make completion recommendations for future wells based on the results obtained. Conclusions and recommended completions were highlighted at the end of this research. It was realized that the natural fractures are the main source of premature water breakthrough, and the Super-K acts as a secondary cause of water channeling to the wellbore.

Recent efforts to optimize oil recovery from Green River reservoirs, Uinta Basin, have stimulated the need for better understanding of the reservoir connectivity at the scale of the operational unit. This study focuses on Green River reservoirs in the Wells Draw study area where oil production response to implemented waterflood is poor and a better understanding of the reservoir connectivity is required to enhance future secondary oil recovery. Correlating the sand bodies between well locations in the area remains difficult at 40-acre well spacing. Thus, interwell connectivity of the reservoirs is uncertain. Understanding the reservoir connectivity in the Wells Draw study area requires integration of all static and dynamic data for generation of probabilistic models of the reservoir at the interwell locations. The objective of this study is two-fold. The first objective was to determine reservoir connectivity at the interwell scale in the Wells Draw study area. To achieve this goal, I used well log and perforation data in the Wells Draw study area to produce probabilistic models of net-porosity for four producing intervals: (1) Castle Peak, (2) Lower Douglas Creek, (3) Upper Douglas Creek, and (4) Garden Gulch. The second objective was to find readily applicable methods for determining interwell connectivity. To achieve this goal, I used sandstone net thickness and perforation data to evaluate interwell connectivity in the Wells Draw study area. This evaluation was done to: (1) assess and visualize connectivity, (2) provide an assessment of connectivity for validating / calibrating percolation and capacitance based methods, and (3) determine flow barriers for simulation. The probabilistic models encompass the four producing intervals with a gross thickness of 1,900 ft and enable simulation assessments of different development strategies for optimization of oil recovery in the Wells Draw study area. The method developed for determining interwell connectivity in Wells Draw study area is reliable and suited to the four producing intervals. Also, this study shows that the percolation based method is reliable for determining interwell connectivity in the four producing intervals.

Høye poretrykk utgjør en risiko for boreoperasjoner på Norsk Sokkel og internasjonalt. Denne risikoen kan reduseres dersom man har kjennskap til poretrykksforholdene før boring. Poretrykk er også en viktig parameter for felt i produksjon, og kunnskap om hvordan dette utvikler seg over tid vil kunne ha stor betydning for økt oljeutvinning. Seismiske data inneholder informasjon om poretrykket og kan derfor bidra til økt kunnskap på dette området. Avhandlingen tar for seg hvordan hastighets- og amplitudeinformasjon fra seismiske data kan brukes for å estimere poretrykk. / Abnormally high pore pressures in the subsurface pose a hazard to drilling operations worldwide. The problem is not unusual on the Norwegian Continental Shelf. Knowledge of the pore pressure prior to drilling may reduce the risk related to drilling in high pressure zones. Pore pressure is also a vital paramter for producinig fields, and knowledge of how the pressure develops over time can be important for increased oil recovery. Seismic data contain information on the pore pressure and may contribute to increased understanding of subsurface pressure conditions. The thesis deals with methods for estimation of pressure from seismic velocity and amplitude data.

In the life of a well, the cement sheath not only provides zonal isolation but also supports casing and increases casing-collapse resistance. Due to the high-pressure, high-temperature (HPHT) conditions, the cement sheath plays an important role in maintaining wellbore integrity. During the production process in HPHT wells, the pressure differential inside the casing and the surrounding formation is larger than the conventional wells. The stress induced by fluid withdrawal in highly compact reservoirs can cause the cement and the casing failure in these wells. These present a greater challenge to the wellbore integrity than the conventional wells. To have reliable data, extensive experimental work on Class G cement was carried out to measure the principal parameters for mechanical structural calculations. The experiment was also set up to simulate conditions under which cement low-cycle fatigue failure could occur. Zero-based cyclic pressure was applied to the casing in the cement low-cycle fatigue test. Three types of cement (72-lbm/ft3, 101-lbm/ft3 and 118-lbm/ft3) were cured and tested at 300ºF to study the cement mechanical properties under high-temperature conditions over the long term. The tests included a 1-year mechanical properties measurement such as compressive strength development; i.e., Young’s modulus and Poisson’s ratio. Finite element methods (FEM) were used to study the casing buckling deformation characteristics of reservoir compaction in some south Texas wells. The 2D and 3D FEM models were built to study the effects of mechanical properties and reservoir compaction on HPHT well integrity. As the confining pressure increases, the cement shows more plasticity and can withstand more pressure cycles. The cement with a higher Poisson’s ratio and lower Young’s modulus showed better low-cycle fatigue behavior. Casing collapse resistance is very sensitive to void location, cement Poisson’s ratio, cement Young’s modulus, and pore pressure. Casing eccentricity and voids shape have minor effect on the casing-collapse resistance. Casing shear failure, tension failure, and buckling failure are the most likely failure modes in reservoir compaction. For different casing wall thickness, the critical buckling strain is almost identical. This study presents a better understanding of casing failure and cement failure in HPHT wells. The results of the study will help improve cement and casing design to maintain wellbore integrity that can in turn be expected to extend throughout the life of the well.