The reserve growth potential of existing conventional oil reservoirs is huge. This research, through numerical simulation, aims to evaluate pattern size reduction as a strategy for improving waterflood recovery in undersaturated oil reservoirs. A plethora of studies have reported improvements in waterflood recovery resulting from pattern size reduction in heterogeneous reservoirs. The dependence of waterflood recovery on pattern size was attributed to factors such as areal reservoir discontinuity, preferential flooding directions, “wedge-edge” oil recovery, irregular pattern geometry, communication with water-bearing zones, vertical reservoir discontinuity, and project economics (Driscoll, 1974). Though many of these publications relied on decline curve analysis in estimating ultimate oil recovery, simulations completed in this thesis support their findings, specifically for compartmentalized reservoirs, fractured reservoirs, and layered reservoirs. Geostatistically-generated permeability fields were employed in the creation of various types of reservoir models. These models were populated with vertical production and injection wells. Sensitivity analysis was then performed on three development scenarios: 160, 40, and 10 acre five-spots. Based on assigned production and injection constraints, the quantity of oil recovered at simulation termination was used to calculate ultimate recovery efficiency. In homogeneous reservoir models, simulation results suggest that waterflood recovery was independent of pattern size. Similar results were also obtained from models with highly-variable non-zero permeabilities. On the other hand, pattern size reduction was found to enhance oil recovery from reservoir models with a high degree of permeability anisotropy. In such reservoirs, recovery was found to be highly dependent on bottom-hole injection pressures. The higher the injection pressure the larger the quantity of oil bypassed by widely spaced patterns. Likewise, high infill potential exists for reservoir models exhibiting areal discontinuity. In these types of models, the improvement in waterflood recovery resulting from pattern size reduction was directly related to the percentage of imbedded zero-permeability grid blocks. Ultimate oil recovery depended on the percolation of permeable grid blocks between production and injection wells. Increasing well density also enhanced waterflood recovery in vertically discontinuous reservoir models. In such layered reservoirs, the amount oil unswept with large patterns was considerably diminished because of the improved injection profiles associated with tighter patterns. / text

The growing demand for energy provides an incentive to pursue unconventional resources. Among these resources, tight gas and shale gas reservoirs have gained significant momentum because recent advances in technology allowed us to produce them at an economical rate. More importantly, they seem likely to contain a significant volume of hydrocarbon. There are, however, many questions concerning hydrocarbon production from these unconventional resources. For instance, in tight gas sandstone, we observe a significant variability in the producibilities of wells in the same field. The heterogeneity is even present in a single well with changes in depth. It is not clear what controls this heterogeneity. In shale gas, the pore connectivity inside the void space is not well explored and hence, a representative pore model is not available. Further, the effects of an adsorbed layer of gas and gas slippage on shale permeability are poorly understood. These effects play a crucial role in assigning a realistic permeability for shale in-situ from a laboratory measurement. In the laboratory, in contrast to in-situ, the core sample lacks the adsorbed layer because the permeability measurements are typically conducted at small pore pressures. Moreover, the gas slippages in laboratory and in-situ conditions are not identical. The present study seeks to investigate these discrepancies. Drainage and imbibition are sensitive to pore connectivity and unconventional gas transport is strongly affected by the connectivity. Hence, there is a strong interest in modeling mercury intrusion capillary pressure (MICP) test because it provides valuable information regarding the pore connectivity. In tight gas sandstone, the main objective of this research is to find a relationship between the estimated ultimate recovery (EUR) and the petrophysical properties measured by drainage/imbibition tests (mercury intrusion, withdrawal, and porous plate) and by resistivity analyses. As a measure of gas likely to be trapped in the matrix during production---and hence a proxy for EUR---we use the ratio of residual mercury saturation after mercury withdrawal (S[subscript gr]) to initial mercury saturation (S[subscript gi]), which is the saturation at the start of withdrawal. Crucially, a multiscale pore-level model is required to explain mercury intrusion capillary pressure measurements in these rocks. The multiscale model comprises a conventional network model and a tree-like pore structure (an acyclic network) that mimic the intergranular (macroporosity) and intragranular (microporosity) void spaces, respectively. Applying the multiscale model to porous plate data, we classify the pore spaces of rocks into macro-dominant, intermediate, and micro-dominant. These classes have progressively less drainage/imbibition hysteresis, which leads to the prediction that significantly more hydrocarbon is recoverable from microporosity than macroporosity. Available field data (production logs) corroborate the higher producibility of the microporosity. The recovery of hydrocarbon from micro-dominant pore structure is superior despite its inferior initial production (IP). Thus, a reservoir or a region in which the fraction of microporosity varies spatially may show only a weak correlation between IP and EUR. In shale gas, we analyze the pore structure of the matrix using mercury intrusion data to provide a more realistic model of pore connectivity. In the present study, we propose two pore models: dead-end pores and Nooks and Crannies. In the first model, the void space consists of many dead-end pores with circular pore throats. The second model supposes that the void space contains pore throats with large aspect ratios that are connected through the rock. We analyze both the scanning electron microscope (SEM) images of the shale and the effect of confining stress on the pore size distribution obtained from the mercury intrusion test to decide which pore model is representative of the in-situ condition. We conclude that the dead-end pores model is more representative. In addition, we study the effects of adsorbed layers of CH₄ and of gas slippage in pore walls on the flow behavior in individual conduits of simple geometry and in networks of such conduits. The network is based on the SEM image and drainage experiment in shale. To represent the effect of adsorbed gas, the effective size of each throat in the network depends on the pressure. The hydraulic conductance of each throat is determined based on the Knudsen number (Kn) criterion. The results indicate that laboratory measurements made with N₂ at ambient temperature and 5-MPa pressure, which is typical for the transient pulse decay method, overestimate the gas permeability in the early life of production by a factor of 4. This ratio increases if the measurement is run at ambient conditions because the low pressure enhances the slippage and reduces the thickness of the adsorbed layer. Moreover, the permeability increases nonlinearly as the in-situ pressure decreases during production. This effect contributes to mitigating the decline in production rates of shale gas wells. Laboratory data available in the literature for methane permeability at pressures below 7 MPa agree with model predictions of the effect of pressure. / text

This study compares four recently developed decline curve methods and the traditional Arps or Fetkovich approach. The four methods which are empirically formulated for shale and tight gas wells are: 1. Power Law Exponential Decline (PLE). 2. Stretched Exponential Decline (SEPD). 3. Duong Method. 4. Logistic Growth Model (LGM). Each method has different tuning parameters and equation forms. The main objective of this work is to determine the best method(s) in terms of Estimated Ultimate Recovery (EUR) accuracy, goodness of fit, and ease of matching. In addition, these methods are compared against each other at different production times in order to understand the effect of production time on forecasts. As a part of validation process, all methods are benchmarked against simulation. This study compares the decline methods to four simulation cases which represent the common shale declines observed in the field. Shale wells, which are completed with horizontal wells and multiple traverse highly-conductive hydraulic fractures, exhibit long transient linear flow. Based on certain models, linear flow is preceded by bilinear flow if natural fractures are present. In addition to this, linear flow is succeeded by Boundary Dominated Flow (BDF) decline when pressure wave reaches boundary. This means four declines are possible, hence four simulation cases are required for comparison. To facilitate automatic data fitting, a non-linear regression program was developed using excel VBA. The program optimizes the Least-Square (LS) objective function to find the best fit. The used optimization algorithm is the Levenberg-Marquardt Algorithm (LMA) and it is used because of its robustness and ease of use. This work shows that all methods forecast different EURs and some fit certain simulation cases better than others. In addition, no method can forecast EUR accurately without reaching BDF. Using this work, engineers can choose the best method to forecast EUR after identifying the simulation case that is most analogous to their field wells. The VBA program and the matching procedure presented here can help engineers automate these methods into their forecasting sheets.