Browsing by Subject "Reservoir"
Now showing 1 - 12 of 12
Results Per Page
Sort Options
Item A column based variance analysis approach to static reservoir model upgridding(Texas A&M University, 2008-10-10) Talbert, Matthew BrandonThe development of coarsened reservoir simulation models from high resolution geologic models is a critical step in a simulation study. The optimal coarsening sequence becomes particularly challenging in a fluvial channel environment where the channel sinuosity and orientation can result in pay/non-pay juxtaposition in many regions of the geologic model. The optimal coarsening sequence is also challenging in tight gas sandstones where sharp changes between sandstone and shale beds are predominant and maintaining the pay/non-pay distinction is difficult. Under such conditions, a uniform coarsening will result in mixing of pay and non-pay zones and will likely result in geologically unrealistic simulation models which create erroneous performance predictions. In particular, the upgridding algorithm must keep pay and non-pay zones distinct through a non-uniform coarsening of the geologic model. We present a coarsening algorithm to determine an optimal reservoir simulation grid by grouping fine scale geologic model cells into effective simulation cells. Our algorithm groups the layers in such a way that the heterogeneity measure of an appropriately defined static property is minimized within the layers and maximized between the layers. The optimal number of layers is then selected based on an analysis resulting in a minimum loss of heterogeneity. We demonstrate the validity of the optimal gridding by applying our method to a history matched waterflood in a structurally complex and faulted offshore turbiditic oil reservoir. The field is located in a prolific hydrocarbon basin offshore South America. More than 10 years of production data from up to 8 producing wells are available for history matching. We demonstrate that any coarsening beyond the degree indicated by our analysis overly homogenizes the properties on the simulation grid and alters the reservoir response. An application to a tight gas sandstone developed by Schlumberger DCS is also used in our verification of our algorithm. The specific details of the tight gas reservoir are confidential to Schlumberger's client. Through the use of a reservoir section we demonstrate the effectiveness of our algorithm by visually comparing the reservoir properties to a Schlumberger fine scale model.Item A composite tracer analysis approach to reservoir characterization(Texas A&M University, 2005-11-01) Oyerinde, Adedayo StephenIn the quest for production optimization from established resources, there is a continual interest in secondary and tertiary recovery methods. The success of these enhanced recovery methods, however, rely to a large extent on a sound understanding of fluid dynamics and migration paths in the reservoir. To this end, several approaches to reservoir characterization have been put to test with varying degrees of success. The unique ability of tracers to provide direct information on preferential fluid flow paths in the reservoir, and the sensitivity of partitioning tracers to fluid saturation distribution has highlighted the prospects of a detailed reservoir characterization through interwell tracer tests. In a broad sense, analysis of interwell tracer tests fall into two categories, analytical and inverse modeling. While most of the analytical methods are laden with limiting assumptions, the method of moments boasts rigorous formulation and accurate estimates of swept volume and average saturation of bypassed oil. The inverse modeling infers permeability and saturation distribution by matching the tracer response. An extremely effective approach to the inverse modeling methods computes sensitivities based on streamlines. The accurate modeling of tracer flow requires accounting for complex phenomena such as transverse dispersion. Also, it is sometimes desired to model pertinent tracer components through compositional simulation. This necessitates the inclusion of a physical dispersion tensor and, hence, the well established finite difference formulation. In this work, we have coupled the finite difference and streamline simulation techniques for the inversion-based reservoir characterization to take advantage of the robustness of the finite difference formulation and computational efficiency of streamline simulation. We have also extended the formalism of the inversion technique for fluid distribution estimation to scenarios with mobile oil saturations and have attempted integrating the analytical and inverse-modeling techniques to facilitate detailed reservoir characterization. We have demonstrated the feasibility of our approach on both synthetic and field cases.Item Adequate description of heavy oil viscosities and a method to assess optimal steam cyclic periods for thermal reservoir simulation(Texas A&M University, 2006-08-16) Mago, Alonso LuisA global steady increase of energy consumption coupled with the decline of conventional oil resources points to a more aggressive exploitation of heavy oil. Heavy oil is a major source of energy in this century with a worldwide base reserve exceeding 2.5 trillion barrels. Management decisions and production strategies from thermal oil recovery processes are frequently based on reservoir simulation. A proper description of the physical properties, particularly oil viscosity, is essential in performing reliable modeling studies of fluid flow in the reservoir. We simulated cyclic steam injections on the highly viscous Hamaca oil, with a viscosity of over 10,000 cp at ambient temperature, and the production was drastically impacted by up to an order of magnitude when using improper mixing rules to describe the oil viscosity. This thesis demonstrates the importance of these mixing rules and alerts reservoir engineers to the significance of using different options simulators have built in their platforms to describe the viscosity of heavy oils. Log linear and power mixing rules do not provide enough flexibility to describe the viscosity of extra heavy oil with temperature. A recently implemented mixing rule in a commercial simulator has been studied providing satisfactory results. However, the methodology requires substantial interventions, and cannot be automatically updated. We provide guidelines to improve it and suggest more flexible mixing rules that could easily be implemented in commercial simulators. We also provide a methodology to determine the adequate time for each one of the periods in cyclic steam injection: injection, soaking and production. There is a lot of speculation in this matter and one of the objectives of this thesis is to better understand and provide guidelines to optimize oil production using proper lengths in each one of these periods. We have found that the production and injection periods should be similar in time length. Nevertheless, the production period should not be less than the injection period. On the other hand, the soaking period should be as short as possible because it is unproductive time in terms of field oil production for the well and therefore it translates into a negative cash flow for a company.Item Construction of a diagenetic history and identification with quality ranking of reservoir flow units: Grayson field, Columbia County, Arkansas(Texas A&M University, 2007-04-25) Poole, Kathleen ReneeThe purpose of this study was to describe depositional and diagenetic characteristics of the (Jurassic) Smackover formation and subsequently identify and rank the quality of flow units within Grayson field, Columbia County, Arkansas. The field has production from the Smackover, a reservoir which consists mainly of highly altered peloidal grainstones. This was a four part study including a lithological analysis of ten cores, a petrographical study of 97 thin sections, a petrophysical study of reservoir properties from core analyses and borehole logs, and predictive mapping of quality ranked flow units across the field. Examination of the cores and thin sections revealed H1a as the main pore type in Grayson field, which was a hybrid of both depositional and diagenetic processes with dominantly interparticle pores. The lowest ranked reservoir quality corresponded to intraparticle and intercrystalline pore types, which occurred mainly in the wackestone/mudstone and packstone/wackestone facies. The highest ranked reservoir quality corresponded to the H1a pore type which occurred mainly in the grainestone/packstone facies 1 and 2. The reservoir quality maps identified the spatial distribution of the facies within the field, which could be used to locate zones for possible in-fill drilling. These results should aid in the economical development of Grayson field and other similar fields.Item Detection of production-induced time-lapse signatures by geophysical (seismic and CSEM) measurements(2011-05) Shahin, Alireza; Tatham, R. H. (Robert H.), 1943-; Stoffa, Paul L., 1948-; Sen, Mrinal; Goff, John; Smith, BrackinWhile geophysical reservoir characterization has been an area of research for the last three decades, geophysical reservoir monitoring, time-lapse studies, have recently become an important geophysical application. Generally speaking, the main target is to detect, estimate, and discriminate the changes in subsurface rock properties due to production. This research develops various sensitivity and feasibility analyses to investigate the effects of production-induced time-lapse changes on geophysical measurements including seismic and controlled-source electromagnetic (CSEM) data. For doing so, a realistic reservoir model is numerically simulated based on a prograding near-shore sandstone reservoir. To account for the spatial distribution of petrophysical properties, an effective porosity model is first simulated by Gaussian geostatistics. Dispersed clay and dual water models are then efficiently combined with other well-known theoretical and experimental petrophysical correlations to consistently simulate reservoir model parameters. Next, the constructed reservoir model is subjected to numerical simulation of multi-phase fluid flow to replicate a waterflooding scenario of a black oil reservoir and to predict the spatial distributions of fluid pressure and saturation. A modified Archie’s equation for shaly sandstones is utilized to simulate rock resistivity. Finally, a geologically consistent stress-sensitive rock physics model, combined with the modified Gassmann theory for shaly sandstones, is utilized to simulate seismic elastic parameters. As a result, the comprehensive petro-electro-elastic model developed in this dissertation can be efficiently utilized in sensitivity and feasibility analyses of seismic/CSEM data with respect to petrophysical properties and, ultimately, applied to reservoir characterization and monitoring research. Using the resistivity models, a base and two monitor time-lapse CSEM surveys are simulated via accurate numerical algorithms. 2.5D CSEM modeling demonstrates that a detectable time-lapse signal after 5 years and a strong time-lapse signal after 10 years of waterflooding are attainable with the careful application of currently available CSEM technology. To simulate seismic waves, I employ different seismic modeling algorithms, one-dimensional (1D) acoustic and elastic ray tracing, 1D full elastic reflectivity, 2D split-step Fourier plane-wave (SFPW), and 2D stagger grid explicit finite difference (FD). My analyses demonstrate that acoustic modeling of an elastic medium is a good approximation up to ray parameter (p) equal to 0.2 sec/km. However, at p=0.3 sec/km, differences between elastic and acoustic wave propagation is the more dominant effect compared to internal multiples. Here, converted waves are also generated with significant amplitudes compared to primaries and internal multiples. I also show that time-lapse modeling of the reservoir using SFPW approach is very fast compared to FD, 100 times faster for my case here. It is capable of handling higher frequencies than FD. It provides an accurate image of the waterflooding process comparable to FD. Consequently, it is a powerful alternative for time-lapse seismic modeling. I conclude that both seismic and CSEM data have adequate but different sensitivities to changes in reservoir properties and therefore have the potential to quantitatively map production-induced time-lapse changes.Item Effect of flue gas impurities on the process of injection and storage of carbon dioxide in depleted gas reservoirs(Texas A&M University, 2005-11-01) Nogueira de Mago, Marjorie CarolinaPrevious experiments - injecting pure CO2 into carbonate cores - showed that the process is a win-win technology, sequestrating CO2 while recovering a significant amount of hitherto unrecoverable natural gas that could help defray the cost of CO2 sequestration. In this thesis, I report my findings on the effect of flue gas ??impurities?? on the displacement of natural gas during CO2 sequestration, and results on unconfined compressive strength (UCS) tests to carbonate samples. In displacement experiments, corefloods were conducted at 1,500 psig and 70??C, in which flue gas was injected into an Austin chalk core containing initially methane. Two types of flue gases were injected: dehydrated flue gas with 13.574 mole% CO2 (Gas A), and treated flue gas (N2, O2 and water removed) with 99.433 mole% CO2 (Gas B). The main results of this study are as follows. First, the dispersion coefficient increases with concentration of ??impurities??. Gas A exhibits the largest dispersion coefficients, 0.18-0.25 cm2/min, compared to 0.13-0.15 cm2/min for Gas B, and 0.15 cm2/min for pure CO2. Second, recovery of methane at breakthrough is relatively high, ranging from 86% OGIP for pure CO2, 74-90% OGIP for Gas B, and 79-81% for Gas A. Lastly, injection of Gas A would sequester the least amount of CO2 as it contains about 80 mole% nitrogen. From the view point of sequestration, Gas A would be least desirable while Gas B appears to be the most desirable as separation cost would probably be cheaper than that for pure CO2 with similar gas recovery. For UCS tests, corefloods were conducted at 1,700 psig and 65??C in such a way that the cell throughput of CO2 simulates near-wellbore throughput. This was achieved through increasing the injection rate and time of injection. Corefloods were followed by porosity measurement and UCS tests. Main results are presented as follows. First, the UCS of the rock was reduced by approximately 30% of its original value as a result of the dissolution process. Second, porosity profiles of rock samples increased up to 2.5% after corefloods. UCS test results indicate that CO2 injection will cause weakening of near-wellbore formation rock.Item Eocene tectonic controls on reservoir distribution in VLE 196, Block V, Lamar Field, Maracaibo Basin, Venezuela(Texas A&M University, 2006-10-30) Choi, ByeonggooIntegrated interpretation of three-dimensional seismic and well-logging data reveals a prominent ??????pop-up?????? structure associated with the VLE 400 fault on the regional unconformity between the Eocene and Miocene in the VLE 196 field, Maracaibo basin, Venezuela. The VLE 400 fault family, an eastern splay of the left-lateral Icotea fault in the basin, played an important role in hydrocarbon migration and accumulation in the field. Hydrocarbons accumulated to the east of the fault but not to the west. The ??????pop-up?????? structure on the Eocene unconformity has a four -way dip closure, straddling the fault and extending to the west of the fault. Structures of the Misoa Formation, which is the main reservoir developed below the unconformity in the basin, differ from the structure of the unconformity. The structure of the Misoa Formation shows a tilted uplift of the eastern block of the fault dipping toward the east caused by thrust tectonic movements. Thrust movement and following strike-slip movements provided additional accommodation space to the west of the fault and generated expanded thickness of Eocene sediments compared to the area east of the fault. The thickness of the Misoa Formation east and west of the fault shows no significant changes. Expanded sediments overlie the Misoa Formation in the western block in lateral contact with Misoa sediments eastern block act as a lateral seal. Ductile movement of the Guasare Formation shale contributed to the lateral sealing of the fault against the reservoir rocks in the eastern block. Mobilization of the Guasare Formation modified the structure of overlying formations including the anticline of the Eocene unconformity. The growth strata provide useful information of reactivation of e xisting faults, especially subtle movements which are not recognized by conventional seismic interpretation. Growth strata isochrons shows subtle reactivation of the VLE 400 fault family during Miocene time.Item An experimental and simulation study of the effect of geochemical reactions on chemical flooding(2010-12) Chandrasekar, Vikram, 1984-; Delshad, Mojdeh; Pope, Gary A.The overall objective of this research was to gain an insight into the challenges encountered during chemical flooding under high hardness conditions. Different aspects of this problem were studied using a combination of laboratory experiments and simulation studies. Chemical Flooding is an important Enhanced Oil Recovery process. One of the major components of the operational expenses of any chemical flooding project, especially Alkali Surfactant Polymer (ASP) flooding is the cost of softening the injection brine to prevent the precipitation of the carbonates of the calcium and magnesium ions which are invariably present in the formation brine. Novel hardness tolerant alkalis like sodium metaborate have been shown to perform well with brines of high salinity and hardness, thereby eliminating the need to soften the injection brine. The first part of this research was aimed at designing an optimal chemical flooding formulation for a reservoir having hard formation brine. Sodium metaborate was used as the alkali in the formulation with the hard brine. Under the experimental conditions, sodium metaborate was found to be inadequate in preventing precipitation in the ASP slug. Factors affecting the ability of sodium metaborate to sequester divalent ions, including its potential limitations under the experimental conditions were studied. The second part of this research studied the factors affecting the ability of novel alkali and chelating agents like sodium metaborate and tetrasodium EDTA to sequester divalent ions. Recent studies have shown that both these chemicals showed good performance in sequestering divalent ions under high hardness conditions. A study of the geochemical species in solution under different conditions was done using the computer program PHREEQC. Sensitivity studies about the effect of the presence of different solution species on the performance of these alkalis were done. The third part of this research focused on field scale mechanistic simulation studies of geochemical scaling during ASP flooding. This is one of the major challenges faced by the oil and gas industry and has been found to occur when sodium carbonate is used as the alkali and the formation brine present in situ has a sufficiently high hardness content. The multicomponent and multiphase compositional chemical flooding simulator, UTCHEM was used to determine the quantity and composition of the scales formed in the reservoir as well as the injection and production wells. Reactions occurring between the injected fluids, in situ fluids and the reservoir rocks were taken into consideration for this study. Sensitivity studies of the effect of key reservoir and process parameters like the physical dispersion and the alkali concentration on the extent of scaling were also done as a part of this study.Item From outcrop to functional reservoir model : using outcrop data to model the tidally dominated esdolomada sandstone, NE Spain(2012-05) Pinkston, Daniel Patrick; Steel, R. J.; Wood, Lesli J.; Olariu, CornelThe Esdolomada Sandstone member 2 crops out in the Tremp-Graus Basin of north-central Spain and forms the uppermost part of the Eocene Roda Formation. The second Sandstone unit within the Esdolomada member (ESD2) consists of bioturbated and shell-rich, very-fine sandstones as well as stacked sets of fine- to coarse-grained cross-stratified sandstones. The overall upward trend in the member is commonly upward thickening and coarsening of beds into and through the cross-stratified interval, though at some few locations there is no obvious trend or even upward thinning of beds. The internal architecture of the member is one in which groups of beds lie between master surfaces that dip highly obliquely to the migration direction of the individual cross strata. The ESD2 is interpreted to be a shelf tidal sand bar within the overall transgressive Esdolomada Sandstone member. It is likely that these bars migrated in a coast parallel fashion, as suggested by the cross-bed orientations, but also accreted laterally away from the coast along the seaward-dipping master surfaces. LIDAR (light detection and ranging) data collection for the Esdolomada member was attempted along the Isábena River near the village of Roda de Isábena, with a total lateral coverage of approximately 3 kilometers. Detailed outcrop measurements were made in accessible areas along the same transect. Outcrop analogs are the best source of data to understand reservoir heterogeneities and to build reservoir analogs for fluid flow simulations. Sand-rich, offshore tidal sandbodies are usually surrounded by marine mudstones, and are recognized from their very orderly stacking of cross-stratified sets (more orderly than in fluvial settings) , their complex internal architecture of master surfaces dipping obliquely to the direction of migration of the contained cross strata and their significant sandstone/mudstone heterogeneities. Tidal bar systems such as the ESD2 are appealing hydrocarbon prospects for several reasons. Primarily, they are relatively coarse grained, have a high degree of lateral continuity, and are relatively clean sands. In places where sand beds are stacked, they create enough thickness to offer good vertical permeability; however, mud-draped cross-beds can create heterogeneities in this type of system that buffer fluid flow. Due to a fairly unsuccessful attempt to obtain LIDAR coverage of the ESD2, in order to build an analog reservoir model, surfaces were instead based on measured sections and outcrop photomosaics. Using Schlumberger’s Petrel software, facies logs were created from measured section data, and then interpolated to make a facies and porosity model.Item Integrated reservoir study of the Monument Northwest field: a waterflood performance evaluation(2009-05-15) Nduonyi, Moses AsuquoAn integrated full-field study was conducted on the Monument Northwest field located in Kansas. The purpose of this study was to determine the feasibility and profitability of a waterflood using numerical simulation. Outlined in this thesis is a methodology for a deterministic approach. The data history of the wells in the field beginning from spud date were gathered and analyzed into information necessary for building an upscaled reservoir model of the field. Means of increasing production and recovery from the field via a waterflood was implemented. Usually, at the time of such a redevelopment plan or scheme to improve field performance, a tangible amount of information about the reservoir is already known. Therefore it is very useful incorporating knowledge about the field in predicting future behavior of the field under certain conditions. The need for an integrated reservoir study cannot be over-emphasized. Information known about the reservoir from different segments of the field exploration and production are coupled and harnessed into developing a representative 3D reservoir model of the field. An integrated approach is used in developing a 3D reservoir model of the Monument Northwest field and a waterflood is evaluated and analyzed by a simulation of the reservoir model. From the results of the reservoir simulation it was concluded that the waterflood project for the Monument Northwest field is a viable and economic project.Item A life cycle optimization approach to hydrocarbon recovery(2010-12) Parra Sanchez, Cristina, 1977-; Lake, Larry W.; Bickel, James E.The objective of reservoir management is to maximize a key performance indicator (net present value in this study) at a minimum cost. A typical approach includes engineering analysis, followed by the economic value of the technical study. In general, operators are inclined to spend more effort on the engineering side to the detriment of the economic area, leading to unbalanced and occasionally suboptimal results. Moreover, most of the optimization methods used for production scheduling focus on a given recovery phase, or medium-term strategy, as opposed to an integrated solution that allocates resources from discovery to field abandonment. This thesis addresses the optimization of a reservoir under both technical and economic constraints. In particular, the method presented introduces a life cycle maximization approach to establish the best exploitation strategy throughout the life of the project. Deterministic studies are combined with stochastic modeling and risk analysis to assess decision making under uncertainty. To demonstrate the validity of the model, this document offers two case studies and the optimal times associated with each recovery phase. In contrast with traditional depletion strategies, where the optimization is done myopically by maximizing the net present value at each recovery phase, our results suggest that time is dramatically reduced when the net present value is optimized globally by maximizing the NPV for the life of the project. Furthermore, the sensitivity analysis proves that the original oil in place and non-engineering parameters such as the price of oil are the most influential variables. The case studies clearly show the greater economic efficiency of this life cycle approach, confirming the potential of this optimization technique for practical reservoir management.Item Simulation of fluid flow mechanisms in high permeability zones (Super-K) in a giant naturally fractured carbonate reservoir(2009-05-15) Abu-Hassoun, Amer H.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.