Porosity and water saturations are the most important petrophysical parameters of hydrocarbon reservoirs that accurate assessment of them in hydrocarbon reservoirs is an effective tool, important and efficient for industry experts, in the context of a comprehensive study of reservoirs and production and management process of reservoir. In this study, using data from five wells of Mansuri oil field, and using the sequential simulation Gaussian method and using Petrel software, the trend of Porosity and water saturation changes in the mentioned field for four zones was simulated. Also the average maps for each zone have been created that results of the simulation parameters in this map showed that highest average porosity is 0.1401 and 0.2756 at least saturation of water is related to zone 1. Finally result of the simulation indicates the Zone 1 is of the best reservoir Zones.
In general, investigation and simulation of petrophysical properties is one of the main aspects of hydrocarbon reservoirs full studies which is performed in most oil companies. Petrophysical properties of reservoir rock such as lithology, porosity, permeability, net pay thickness and percentage of water and oil saturation obtained from well logs interpretations can be modeled by simulation of the reservoir [
Petrophysical modeling is performed by two methods: deterministic and stochastic. SGS method is a stochastic simulation method which implements Geostatistical rules to predict and simulate reservoir parameters in three dimensions of the reservoir. The present study investigates the modeling and the petrophysical parameter of the Asmari formation in Mansuri oil field: South western of Iran (
Modeling software’s present a three-dimensional view of all aspects of the reservoir, such as geology and structure of reservoir, fluid dynamics and well [
One of the most efficient methods to estimate and simulate petrophysical parameters contained in these software’s is SGS method that depends on wells and geophysical data [
Typically, structural modeling of reservoir, is the starting point of the reservoir modeling [
Corner point gridding of the reservoir: the cells in this gridding method are very flexible and can easily show faulting effects. The other advantage of these cells is the ability to create the reservoir flow simulation cells. Generally, these cells are more suitable for geological features and because of this advantage, in this study, these cells were used for gridding of Asmari reservoir. But in the layering of cells proportional method was used in which the number of cell layers was considered constant in the reservoir and their top and bottom surfaces was considered parallel to the top and bottom surfaces of the reservoir layer. The size of each cell was 1 × 120 × 120 meters.
After data up scaling, two fundamental conditions in Geostatistical calculations are: first, the data has a normal distribution and second, the data shows no trend. The spatial relationship between the values of a quantity in taken samples population can be expressed in the form of spatial structure [
Petrophysical data are in the form of logging data in the well trajectory. Before the simulation, the parameters must be averaged in the cell size defined in the reservoir. This gives a value for different parameters for each cell. There are different methods to calculate the average of parameters, numerical averaging method is the best one due to static nature of porosity and water saturation. The formula of numerical averaging method is as follow:
where KA is the average of petrophysical parameter K, Ki is the value of parameter K in the position i and N is the number of samples. In Mansuri oil field using the 5 wells drilled in different locations of the south area of reservoir (
In order to implement stochastic methods, existing trends in the data which might be created due to geological characteristic of the region, must be identified and removed. In order to remove the trends in the data, first, the data must be normalized and then kurtosis of the normal distribution must be increased as much as possible until a range of data which is used accumulate largely around the mean data. Then, no trend remains in the kurtosis of the data. One of other
conditions to implement stochastic method for simulation of petrophysical pa- rameters is to have standard normal distribution of the data. In normal distribution of the data, mean is equal to zero, standard deviation is equal to 1 and their distribution diagram is a Gaussian distribution. Qualitatively, normal distribution of data can be detected using formation statistical tools such as Histograms, skewness and kurtosis of data distribution, based on which the porosity and water saturation data were normalized (
In Geostatistics using the data of a quantity in a known coordination, it is possible to estimate the value of the quantity in a different spot with a known coordinates located within the range of the spatial structure [
where h is variograms value, N(h) is the number of data pairs and xi − xi + h is beginning and end of data pair i.
In general, the variograms in which one axis is variance and the other one is the distance, shows the variability of data versus the distance [
gular and symmetrical arrangement of wells and more homogeneous reservoir, the less variance of the estimation of petrophysical parameters of reservoir rock.
Selection of horizontal directions of the variogram is performed based on the acoustic impedance map which reflects the homogeneity and the heterogeneity of reservoir. So that the major axis of the variogram (X) is in the more homogeneous (low variance) direction and its length is half length of the reservoir and minor axis (Y) is perpendicular to this axis and its length is equal to half the width of the reservoir (
Having defined the mentioned condition, the petrophysical data has been simulated by SGS method which was presented in separate SGS models for each of the fore zones (
Zone | Porosity | Permeability |
---|---|---|
Zone 1 | 0.1401 | 0.2756 |
Zone 2 | 0.0970 | 0.4667 |
Zone 3 | 0.0771 | 0.3915 |
Zone 4 | 0.0642 | 0.6015 |
integrated model of all layers has also been built. Generally, the porosity and water saturation variation maps shows that in porosity values distribution map, its maximum value starts at Southeast of the field, passing through the center it ends at Northwest of the field (
Water saturation model that shows the value of water saturation in oil zones, is mainly of the irreducible water saturation, which is affected by the reservoir capillary pressure, relative permeability and porosity. To build water saturation model in the reservoir co-kriging of reservoir porosity has been used as a secondary parameter. In this model, maximum water saturation belongs to North and Northwest of the field (
To assess the model accuracy, the results should be compared with actual data in wells (validation) (
more accurate will be the developed model which in the Asmari model very good agreement is observed between the model results and data for porosity (
Finally, to identify potential areas of the reservoir, applying cut off limits for two parameters of porosity and water saturation, reservoir and non-reservoir areas model has been obtained.
Cut off limits for porosity and water saturation were considered 5% and 30%, respectively. This means that according to the production characteristics, reservoir areas should be above 5% of porosity and below 30% of water saturation. The model shows that except the southern, central and a part of Southwestern areas, the field has suitable reservoir characteristic which can be taken into account in development plans.
In the present study, the SGS method, has been implemented to produce a 3D
petrophysical model of Asmari reservoir of Mansuri oil field (located in Southwestern of Iran). SGS method is an efficient method to estimate and simulate petrophysical parameters. The petrophysical properties i.e. porosity and water saturation have been used for petrophysical modeling. The data from 5 wells located in Southern section of the field has been used, therefore the results of simulation are limited to Southern area of the studied field. The software, Petrel has been used in this study, which covered whole reservoir simulation stages from the interpretation of seismic data to building static and dynamic models of hydrocarbon reservoirs. Due to high lithology heterogeneity in vertical or well drilling direction in Mansuri oil field, petrophysical parameters change in this direction while because of lithology continuity, variations of these parameters in the horizontal direction are less than those of vertical direction. In data analysis, the more organized well pattern and symmetrical distribution, the less estimation variance for reservoir rock parameters and more accurate spatial structure for the data. The maximum value of porosity starts from south and decreases to the north of the studied section. Also, the maximum amount of water saturation can be found in the north of the section. Zone 1 is the best reservoir section which has the highest value for the porosity and the lowest amount of water saturation. Zone 4 has the lowest reservoir value because of poor porosity and high water saturation. The results of models verification showed that estimation and simulation of the porosity and water saturation using sequential Gaussian Simulation Method (GSM) has been properly performed. Finally, to identify potential areas of the reservoir, 5% and 30% were applied as cut off limits for porosity and water saturation, respectively. The resulting model showed that the southern part possesses better reservoir characteristics comparing to other parts of the studied section.
We would like to express our sincere thanks to the Department of Oil Engineering, Gachsaran Branch, Islamic Azad University, and all our colleagues for critical discussions and support. Thanks also to the anonymous referees for all their gracious and critical points.
Nezhad, H.K. and Tabatabaei, H. (2017) Simulation of Petrophysical Parameters of Asmari Reservoir Using SGS Method in Mansuri Oil Field, Southwest of Iran. Open Journal of Geology, 7, 1188-1199. https://doi.org/10.4236/ojg.2017.78079