Application of 3D Reservoir Geological Model on Es1 Formation, Block Nv32, Shenvsi Oilfield, China ()
1. Introduction
When the oil has been discovered in a field, many studies will be required to evaluate and comprehend the reservoir heterogeneity, define the extent of the reservoir in the three-dimensional, evaluating the fluid volumetric in the reservoir and identify the suitable method to increase the reservoir fluid recovery. Generally, reservoir characteristics such as lithofacies heterogeneities, porosity, and permeability spatial variation control the reservoir performance, and development plans [1]. The 3D reservoir model represents the reservoir properties, thus building a static reservoir involve to accomplish four tasks which care out by specialists in multiple disciplines, these stages are Structural modeling, Stratigraphic modeling, Lithological modeling, and Petrophysical modeling. Typically, models are not exactly representing the reality, however, merely a computer-aided used to illustrate the property distribution of the reservoir characteristics to assist in the forecast of the reservoir forthcoming result. The reservoir model is an essential tool to identify the appropriate, safest drilling, completion, and recovery selection for the reservoir, and for the field development plane for that reservoir as well as [2]. This research aims to characterize and carry out 3D static modeling of Es1 formation, block Nv32 to make economic management decisions aimed towards the profitability of an oil company by reducing the capital investment by providing encouragement and confidence to the oil investment with the following objectives: Correlate the reservoir across the 22 wells, delineate the hydrocarbon bearing reservoir map out major faults within the field, compute the petrophysical parameters such as porosity, permeability net-to-gross ratio and saturation using the deterministic approach.
2. Study Area
Block Nv32 is located in the west of the Shenvsi oilfield (Dagan oilfield) in the south of Cangzhou City, Hebei province and its structure locate in between Cangxian and Kongdian structures. Due to the complex structure and fault development of the Shenvsi Oilfield, there are many blocks in the oil field. Moreover, block Nv32 mainly contains sandstone and carbonate rocks, sandstone with lithic argillaceous sandstone, feldspathic lithic sandstone, and a small amount of mixed sandstone. The structural area of the sand-oil top structure is 2.0 km2, the structural width is 100 m, the proven oil-bearing area is 1.4 km2, and the petroleum geological reserves are 225 × 104 t. Figure 1 illustrates the location map of block Nv32.
Figure 1. Map location of Es1 formation, block Nv32.
3. Model Workflow
The data used for this research were obtained from block Nv32 in China. These data include structural maps, fault lines, and well logs. Data were checked and organized properly, and then imported into Petrel software to build the 3D reservoir model. The workflow design used for the study and wide range of functional tools in the petrel software includes 3D visualization, well correlation, 3D mapping, and 3D grid design for geology simulation, well log upscaling, petrophysical modeling, data analysis, and volume calculation. Figure 2 illustrates the workflow model of Es1 formation.
Figure 2. Workflow model of Es1 formation.
4. Data Preparation
Data preparation is the foremost stage for constructing a reservoir static model. Geomap3.6 and Petrel software are the chiefly software to construct this model. Based on Petrel’s demand and area characteristics of Es1 formation, block Nv32. Three types of data were required to construct this reservoir static model which can be concluded as: well path data, structure data, and property data. Well path data to describe the pathway of wells across, while structure data is used to create a 3D grid, which includes individual layers data and fault data, on the other hand; Property data is used to generate property model by applying numerical analysis calculation to the data by property definition [3].
5. Well Correlation
Wells correlation has been performed to provide clear depict about the reservoir geological characteristic such as thickness change within Es1 formation units block Nv32, and the change of petrophysical properties such as a change in (Porosity, permeability, oil saturation, and N/G) of different units in Es1 formation. In accordance with data input into Petrel software, the correlation section for twenty-two well was performed. The sequence was subdivided into a lithological unit namely Shale, Sandstone, Limestone according to the sand characteristics of Es1 formation, Nv32 block. The correlating process of Es1 formation wells was made by using different logs such as Spontaneous potential (SP), Caliper (CAL) and Acoustic impedance (AC) logs to perform the well to well correlation. After data were inputted into Petrel software, the correlation section of Es1 formation was performed. Figure 3 illustrates the vertical and horizontal variations in the thickness of Es1 formation and the variations in petrophysical properties as well.
6. Structural Modeling
Structural modeling in Petrel is subdivided into three operations which can be included as follows: Fault modeling, Pillar gridding, and vertical layering. All the mentioned processes were performed sequentially to form a single data model [4]. A structure contour map is an essential tool for interpretation three-dimensional because it provides obvious depictions for the complete three-dimensional form of map horizon. Three-dimensional were constructed a base on the well top for all wells of Es1 formation, block Nv32. Figure 4 illustrates the 2D map of the Es1 formation top horizon.
6.1. Fault Model
The fault model was the primary step for constructing a structure model with petrel workflow tools. Es1 formation mainly consists of fifteen faults. Fault type in bock Nv32 is normal and reverse fault. This fault model has used the fault polygon which Petrel software provides to each type of vertical fault, diagonal fault, curve fault, and various geometry structure fault according to the fault polygon. The initial surface was established by applying fault polygon with various stratification planes, then fault section was corrected by adjusting the fault section, uses breakpoint data obtained by single well correlation, correct the position and shape of fault revision fault and shape, make it convenient and anastomosing. Figure 5 illustrates the fault model of Es1 formation, block Nv32.
Figure 3. Well correlation of Es1 formation, block Nv32.
Figure 4. 2D map of Es1 formation top horizon.
Figure 5. Fault model of Es1 formation, block Nv32.
6.2. The Layering of Geological Model
Identify the thickness and the directions of layers among the horizons of this 3D grid was the final step to construct the structure framework. Moreover, Modern geology required an accurate illustration of layered volumes. Recently, the three-dimensional geological model is the most effective method for constructing the geological model at depth [5]. Each unit in Es1 formation, block Nv32 had been divided into many layers based on the lithological, Petrophysical properties.Es1 × 1 and Es1 × 2 consist of 11 layers, however, Es1 × 3 and Es1 × 4 consist of 10layers which based on the petrophysical properties. Figure 6 illustrates the layering of Es1 formation, block Nv32.
6.3. Facies Modeling
Facies modeling is discrete throughout the model grid. Facies model of Es1 Formation, block Nv32 has constructed the base on the outcome result of the microfacies interpretation of sedimentary environments as obtained from logs. The Es1 formation primarily consists of three types of lithology which are shale, limestone, and sandstone with a proportion of 32.42%, 29.49%, and 38.09% respectively. However, limestone is the predominated reservoir rock in block Nv32. The assigned value method was applied to simulate the sand bodies in Es1 formation. Figures 7-10 illustrate the facies model for each zone of Es1 formation, block Nv32.
Figure 6. Layering of Es1 formation, block Nv32.
6.4. Petrophysical Modeling
Petrophysical property modeling is known as the process of distributing reservoir properties into the geological grid. This is performed individually for each reservoir zone of the geological grid [6].
6.5. Porosity Model
Evaluate oil capacity in the oil reservoir is achieved by determining its porosity. The porosity model of Es1 formation was built based on the outcome result that has been determined from porosity which was obtained from the petrophysical interpretation of block Nv32 wells. In this model, the porosity model was constructed for each unit of Es1 formation, block Nv32. The arithmetic average method was applied to scale up the well logs. Sequential Gaussian Simulation method was used to distribute the porosity in the model. The porosity occurrence distribution is a concern between 0.2% - 36.39% of Es1 formation, block Nv32 with an average porosity value of 17.5%. Table 1 presents the outcome result which indicates that Es1 formation, block Nv32 is a good reservoir. Figure 11 represents the porosity histogram of Es1 formation which indicates that the porosity of Es1 formation is a range between 0.2 - 0.36. Porosity reduces in some parts because of the clay contents, and the heterogonies. Figures 12-15 represent the porosity model of Es1 formation, block Nv32.
Figure 11. Porosity histogram of Es1 formation.
Figure 12. Porosity model zone Es1 × 1, block Nv32.
Figure 13. Porosity model zone Es1 × 2, block Nv32.
Figure 14. Porosity model zone Es1 × 3, block Nv32.
Figure 15. Porosity model zone Es1 × 4, block Nv32.
6.6. Permeability Model
Permeability is an essential parameter in the reservoir rock. Its measure of the ease with which a fluid flows through the connecting pore space of the reservoir rock. The reservoir rock permeability is a vital property to determine the orientation movement and flow rate of the reservoir fluid in the formation. In this model, the permeability model was based on the permeability log generated [7]. The harmonic average method was applied to scale up the well logs. Furthermore, Sequential Gaussian Simulation was used as the statistical method to create the permeability model. The permeability model shows that the permeability of Es1 formation, block Nv32 is concentrated between 0.017 mD to 974.8 mD, having an average permeability of 59.44 mD. The outcome value is indicated that Es1 formation has good permeability, thus Es1 formation is considered as a good reservoir rock. Figures 16-19 represent the permeability model of Es1 formation, block Nv32.
6.7. Oil Saturation Model
The saturation distribution model is known as the percentage between oil, water, and gas in the reservoir pore space. Saturation is an essential tool to identify a high prospective water area [8]. In the absence of core data, a deterministic of oil saturation was achieved from petrophysical interpretation. Sequential Gaussian Simulation was used to construct the oil saturation model of Es1 formation, block Nv32. Figures 20-23 show a 3D perspective view of oil saturation. The saturation model shows that the saturation of Es1 formation is mainly concentrated between 0.00% to 95.9%, having an average saturation of 22.2%. Es1 × 1 zone reveals that the oil saturation is relatively concentrated in northwestern and southeastern. However, oil saturation in Es1 × 2 and Es1 × 3 are concentrated in the west to southwestern. Moreover, Es1 × 4 is highly concentrated in the west to southwestern, also relatively oil saturation is concentrated in the north part.
6.8. Net to Gross Model
Net to Gross ratio is used to eliminate non-productive reservoir rock unit if exist. Moreover, limestone is the predominant reservoir rock of Es1 formation. Thus, Net to gross was set to unity for the volume calculation. The arithmetic average method was applied to scale up the well logs. Moreover, the Sequential Gaussian Simulation method was used to build the Net to Gross model of Es1 formation, block Nv32. Net to Gross model shows that the N/G of Es1 formation is primarily concentrated between 0.01 and 1.00 within an average value of 0.61. Table 2 represents the outcome result of N/G. Figures 24-27 represent the Net to Gross model.
Table 2. Net to gross model result.
7. Reservoir Volumetric
Estimate the quantity of hydrocarbon in the reservoir is known as reservoir volumetric. Reservoir volumetric is an essential part of the geological model because it provides the basis for oilfield exploration and development. Furthermore, reservoir volumetric can examine the reliability of the model [8]. In this model, various data were analyzed and bounded to create a three-dimensional geological model. The Es1 formation, block Nv32 was initially defined by its top formation and base formation, four main horizons, and 42 sub-layers. Petrophysical analysis of all properties and Oil-Water contact was determined. The oil formation factor which was obtained from the laboratory from the PVT analysis was used as 1.21 bbl/STB to calculate hydrocarbon volume initially in place in standard condition (14.7 Psi, 60 F). Also, when the reservoir model of Es1 formation was completed, the structure model and property model were used to calculate oil reserve in term of Stock Tank Oil Initial in Place (STOIIP). The Es1 formation was estimated by applying the equation [9] - [17]. Table 3 shows the volumetric outcome result after modeling has been done.
Table 3. Volumetric calculation results.
Table 3 indicates that oil in zone Es1 × 3 and Es1 × 4 is high and commercial to produce from these zones. Figure 28 and Figure 29 illustrate the oil-bearing zone for Es1 × 3 and Es1 × 4 which reveal a huge area of the oil-bearing zone for the last two zones.
Figure 28. Oil-bearing area for zone Es1 × 3.
Figure 29. Oil-bearing area for zone Es1 × 4.
8. Conclusions
The overall values such as porosity and permeability indicate that Es1 formation, block Nv32 is a good and promising oilfield, because of its good porosity and permeability values.
1) The structural model has been completed for Es1 formation, block Nv32. This model describes the structure trap characteristic, identify the distribution, geometry structure, reservoir attitude, and development of faults and fractures of Es1 formation, block Nv32.
2) Horizons were constructed for Es1 formation, block Nv32, and they are divided into 4 zones namely as Es1 × 1, Es1 × 2, Es1 × 3, and Es1 × 4. Based on the facies, petrophysical properties, layers were built for each zone. Es1 × 1, Es1 × 2 zones were divided into 11 layers. However, Es1 × 3 and Es1 × 4 were divided into10layers, and this was based on the thickness and petrophysical properties.
3) Facies modeling for block Nv32 was built using the assigned value method. This model shows that Es1 formation consists of three types of lithology which are shale, limestone, and sandstone with a proportion of 32.42%, 29.49%, and 38.09% respectively. However, limestone is the predominated rock in this reservoir.
4) The petrophysical model (porosity, permeability, oil saturation, and Net to Gross) of Es1 formation was built based on their values using the Sequential Gaussian Simulation method as a statistical method. Es1 × 3 and Es1 × 4 are the highest petrophysical properties in terms of porosity, permeability.
Acknowledgements
My deepest thanks and gratitude go to my parents, brothers, and professor Guan Zhenliang for their high support and motivation. I would like to convey my special thanks to Dr. Atif Zafar for his sincere advice and continuous support as well as constructive comments which helped me to think critically.