Reservoir heterogeneities play a crucial role in governing reservoir performance and management.Traditionally,detailed and inter-well heterogeneity analyses are commonly performed by mapping seismic facies change in t...Reservoir heterogeneities play a crucial role in governing reservoir performance and management.Traditionally,detailed and inter-well heterogeneity analyses are commonly performed by mapping seismic facies change in the seismic data,which is a time-intensive task.Many researchers have utilized a robust Grey-level co-occurrence matrix(GLCM)-based texture attributes to map reservoir heterogeneity.However,these attributes take seismic data as input and might not be sensitive to lateral lithology variation.To incorporate the lithology information,we have developed an innovative impedance-based texture approach using GLCM workflow by integrating 3D acoustic impedance volume(a rock propertybased attribute)obtained from a deep convolution network-based impedance inversion.Our proposed workflow is anticipated to be more sensitive toward mapping lateral changes than the conventional amplitude-based texture approach,wherein seismic data is used as input.To evaluate the improvement,we applied the proposed workflow to the full-stack 3D seismic data from the Poseidon field,NW-shelf,Australia.This study demonstrates that a better demarcation of reservoir gas sands with improved lateral continuity is achievable with the presented approach compared to the conventional approach.In addition,we assess the implication of multi-stage faulting on facies distribution for effective reservoir characterization.This study also suggests a well-bounded potential reservoir facies distribution along the parallel fault lines.Thus,the proposed approach provides an efficient strategy by integrating the impedance information with texture attributes to improve the inference on reservoir heterogeneity,which can serve as a promising tool for identifying potential reservoir zones for both production benefits and fluid storage.展开更多
Oil saturation is a critical parameter when designing oil field development plans.This study focuses on the change of oil saturation during water flooding.Particularly,a meter-level artificial model is used to conduct...Oil saturation is a critical parameter when designing oil field development plans.This study focuses on the change of oil saturation during water flooding.Particularly,a meter-level artificial model is used to conduct relevant experiments on the basis of similarity principles and taking into account the layer geological characteristics of the reservoir.The displacement experiment’s total recovery rate is 41.35%.The changes in the remaining oil saturation at a millimeter-scale are examined using medical spiral computer tomography principles.In all experimental stages,regions exists where the oil saturation decline is more than 10.0%.The shrinkage percentage is 20.70%in the horizontal well production stage.The oil saturation reduction in other parts is less than 10.0%,and there are regions where the oil saturation increases in the conventional water flooding stage.展开更多
The Lower Cretaceous Manville Group of Upper Mc Murray Formation is one of the main bitumen reservoirs in Athabasca.In this study,the relationship between reservoirs heterogeneity and bitumen geochemical characteristi...The Lower Cretaceous Manville Group of Upper Mc Murray Formation is one of the main bitumen reservoirs in Athabasca.In this study,the relationship between reservoirs heterogeneity and bitumen geochemical characteristics were analyzed through core and microscopic observation,lab analysis,petrophysics and logging data.Based on the sedimentology framework,the formation environment of high-quality oil sand reservoirs and their significance for development were discussed.The results indicate that four types lithofacies were recognized in the Upper Mc Murray Formation based on their depositional characteristics.Each lithofacies reservoirs has unique physical properties,and is subject to varying degrees of degradation,resulting in diversity of bitumen content and geochemical composition.The tidal bar(TB)or tidal channel(TC)facies reservoir have excellent physical properties,which are evaluated as gas or water intervals due to strong degradation.The reservoir of sand bar(SB)facies was evaluated as oil intervals,due to its poor physical properties and weak degradation.The reservoir of mixed flat(MF)facies is composed of sand intercalated with laminated shale,which is evaluated as poor oil intervals due to its poor connectivity.The shale content in oil sand reservoir is very important for the reservoir physical properties and bitumen degradation degree.In the context of regional biodegradation,oil sand reservoirs with good physical properties will suffer from strong degradation,while oil sand reservoirs with relatively poor physical properties are more conducive to the bitumen preservation.展开更多
Currently,limited studies of immiscible water-alternating-CO_(2)(imWACO_(2))flooding focus on the impact of reservoir heterogeneity on reservoir development outcomes.Given this,using the heterogeneous reservoirs in th...Currently,limited studies of immiscible water-alternating-CO_(2)(imWACO_(2))flooding focus on the impact of reservoir heterogeneity on reservoir development outcomes.Given this,using the heterogeneous reservoirs in the Gao 89-1 block as a case study,this study conducted slab core flooding experiments and numerical simulations to assess the impact of reservoir heterogeneity on imWACO_(2)flooding efficiency.It can be concluded that imWACO_(2)flooding can enhance the sweep volume and oil recovery compared to continuous CO_(2)flooding.As the permeability difference increases,the difference in the swept volume between zones/layers with relatively high and low permeability increases.To optimize the exploitation of reservoirs in the Gao 89-1 block,the optimal timing and CO_(2)injection rate for imWACO_(2)flooding are determined at water cut of 40%and 10000 m^(3)/d,respectively.A short injection-production semi-period,combined with multiple cycles of water and CO_(2)injection alternations,is beneficial for enhanced oil recovery from imWACO_(2)flooding.展开更多
The processes of flooding—water flooding, polymer flooding and ternary combination flooding—were simulated respectively on a 2-D positive rhythm profile geological model by using the ASP numerical modeling softw...The processes of flooding—water flooding, polymer flooding and ternary combination flooding—were simulated respectively on a 2-D positive rhythm profile geological model by using the ASP numerical modeling software developed by RIPED (Yuan, et al. 1995). The recovery coefficient, remaining oil saturation, sweep efficiency and displacement efficiency were calculated and correlated layer by layer. The results show that the sweep efficiency and displacement efficiency work different effects on different layers in the severely heterogeneous reservoir. The study shows that the displacement efficiency and sweep efficiency play different roles in different layers for severely heterogeneous reservoirs. The displacement efficiency contributes mainly to the high permeability zones, the sweep efficiency to the low permeability zones, both of which contribute to the middle permeable zones. To improve the sweep efficiency in the low permeability zones is of significance for enhancing the whole recovery of the reservoir. It is an important path for improving the effectiveness of chemical flooding in the severely heterogeneous reservoirs to inject ternary combination slug after profile control.展开更多
Taking reservoir rocks and fluids of the Daqing,Dagang and Changqing oilfields as research objects,the EOR mechanisms and technical approach of polymer flooding were discussed.By comparing the displacement performance...Taking reservoir rocks and fluids of the Daqing,Dagang and Changqing oilfields as research objects,the EOR mechanisms and technical approach of polymer flooding were discussed.By comparing the displacement performances of ordinary polymer,glycerol,polymer in"sheet-net"structure and heterogeneous weak gel at the same viscosity and concentration,the relationship between the viscosity of polymer displacement agents and displacement performance was demonstrated,and the method of improving polymer flooding effect was worked out.The main mechanism of polymer flooding to increase oil recovery is the swept volume expansion of water injection due to polymer retention in porous media.The viscosity of polymer agents has no positive correlation with polymer flooding effect.Although polymer of"sheet-net"structure has strong capacity in increasing viscosity,it has poor compatibility with pore throat structure of reservoir rock,low injectivity and low shear resistance.Heterogeneous weak gel system has higher adsorption and capture capacity in porous media,which is easy to retain in porous media,and can effectively establish seepage resistance in high permeability layers(zones).Compared with polymer solutions with the same viscosity or concentration,it has stronger ability to expand swept volume.Long term injection of polymer flooding agents will inevitably lead to fluid entry profile reversal,and thus worsening of polymer flooding effect.Alternate injection of high retention and low or non-retention displacement agents can further improve the displacement effect of polymer flooding agents.展开更多
During the displacement of water plugging with binary flooding in internally heterogeneous reservoirs,it is essential to understand the distributions of remaining oil as well as the oil displacement mechanisms at diff...During the displacement of water plugging with binary flooding in internally heterogeneous reservoirs,it is essential to understand the distributions of remaining oil as well as the oil displacement mechanisms at different stages.In this study,two types of internally heterogeneous systems,i.e.,vertical and horizontal wells are investigated experimentally through a microscopic approach.The results show that plugging agent types have a greater impact on oil recovery than well types,and foam injection can enhance oil recovery more effectively than gel injection.Additionally,the injection sequence of plugging agents significantly affects oil displacement efficiency.Injecting gel after foam is more beneficial.According to the present results,the main formation mechanisms of remaining oil in each displacement stage are influenced by:capillary force,viscous force,inertial force,shear force,microscopic fingering&channeling.展开更多
In view of the anisotropy,heterogeneity and stress-sensitive permeability in low permeability reservoirs,an analytical well test model was established by introducing the concept of permeability modulus.This model cons...In view of the anisotropy,heterogeneity and stress-sensitive permeability in low permeability reservoirs,an analytical well test model was established by introducing the concept of permeability modulus.This model considered the permeability stress-sensitivity,wellbore storage effect,and the skin effect.The perturbation technique and Laplace transformation were used to solve the mathematical model analytically in Laplace space,and the bottom-hole pressure type curves were plotted and analyzed in real space by using the Stehfest numerical inversion.展开更多
Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information whe...Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information when the reservoir exhibits reservoir heterogeneity and multiphase flow simultaneously.This work proposes a practical-yet-rigorous method to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,and improve the evaluation of reservoir/fracture properties.A new,general,semi-analytical model is proposed that explicitly accounts for multiphase flow,fractalbased reservoir heterogeneity,anomalous diffusion,and pressure-dependent fluid properties.This is achieved by introducing a new Boltzmann-type transformation,the exponent of which includes reservoir heterogeneity and anomalous diffusion.In order to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,the modified Boltzmann variable allows the conversion of three partial differential equations(PDE's)(i.e.,oil,gas and water diffusion equations)into ordinary differential equations(ODE's)that are easily solved using the Runge-Kutta(RK)method.A modified time-power-law plot is also proposed to estimate the reservoir and fracture properties,recognizing that the classical square-root-of-time-plot is no longer valid when various reservoir complexities are exhibited simultaneously.Using the slope of the straight line on the modified time-power-law plot,the linear flow parameter can be estimated with more confidence.Moreover,because of the new Boltzmann-type transformation,reservoir and fracture properties can be derived more efficiently without the need for defining complex pseudo-variable transformations.Using the new semi-analytical model,the effects of multiphase flow,reservoir heterogeneity and anomalous diffusion on rate-decline behavior are evaluated.For the case of approximately constant flowing pressure,multiphase flow impacts initial oil rate,which is a function of oil relative permeability and well flowing pressure.However,multiphase flow has a minor effect on the oil production decline exponent.Reservoir heterogeneity/anomalous diffusion affect both the initial oil production rate and production decline exponent.The production decline exponent constant is a function of reservoir heterogeneity/anomalous diffusion only.The practical significance of this work is the advancement of RTA techniques to allow for more complex reservoir scenarios,leading to more accurate production forecasting and better-informed capital planning.展开更多
This paper demonstrates the use of a commercial simulator as a tool with which to optimize the SAGD (steam-assisted gravity drainage) start-up phase process. The factors affecting the start-up phase are the prime ta...This paper demonstrates the use of a commercial simulator as a tool with which to optimize the SAGD (steam-assisted gravity drainage) start-up phase process. The factors affecting the start-up phase are the prime targets. Among the key investigated factors are wellbore geometry effects, reservoir heterogeneity and circulation phase length. Each of the parameters was investigated via steam chamber development observation along the well pair length and the cross sections in the mid, toe and heel areas. In addition, the cumulative recovery in given time, steam-to-oil ratio and CDOR (calendar day oil rate) production data are used to backup the observations produced in the simulated model. Furthermore, an additional component developed during the research is a statistical modification of a layer cake model with which to create a heterogeneous reservoir to represent real reservoir conditions, based on Monte Carlo's simulation.展开更多
The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a not...The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.展开更多
Terminal fans have formed the sedimentary system of the 2+3 sands of the upper second member, Shahejie formation in the west of the Pucheng Oilfield, Bohai Bay Basin, East China. Based on well logging data and physic...Terminal fans have formed the sedimentary system of the 2+3 sands of the upper second member, Shahejie formation in the west of the Pucheng Oilfield, Bohai Bay Basin, East China. Based on well logging data and physical properties of the reservoir beds, the 2+3 sands were divided into 16 sublayers. The heterogeneity of reservoir beds and distribution of interlayers and seal layers in the 2+3 sands were investigated. The intra-layer heterogeneity and inter-layer heterogeneity primarily belong to the severely heterogeneous classification. The spatial differentiation of sedimentary microfacies resulted in a change of reservoir bed heterogeneity, strong in the middle and southern parts, weak in the northern part. Spatial distribution of interlayers and seal layers is dominated by sedimentary microfacies, and they are thick in north-eastern and middle parts, thin in the south-western part.展开更多
An analysis of drill cores and well logs shows that the main micro-facies of the third member sand bodies of the Qingshankou Formation in Qian'an are subaqueous distributary channel facies, sheet sand facies and suba...An analysis of drill cores and well logs shows that the main micro-facies of the third member sand bodies of the Qingshankou Formation in Qian'an are subaqueous distributary channel facies, sheet sand facies and subaqueous fan facies (olistostrome). Maps showing the distribution of these micro-facies together with inter-channel bay and prodelta mocro-facies are presented for different time-slices (lower, middle and upper parts of the Qingshankou Formation). These maps reveal the instability and change of sediment transport in the Baokang sedimentary system during the depositional period. Sediment transport was from the west in the early stage, from the south in the middle stage and from the northwest in the late stage. Values of thickness, porosity and permeability of the sand bodies in the third member of the Qingshankou Formation show that they have low to medium porosity and low permeability, and are characterized by serious reservoir heterogeneity. The joints between micro-facies and subaqueous fan micro-facies are characterized by the highest heterogeneity, the sheet sand and distal sand bar subfacies come next, and the heterogeneity of the subaqueous distributary channel sand bodies is relatively weak.展开更多
The distribution characteristics of the oil-water contact are the basis for the reservoir exploration and development and reserves evaluation. The reservoir with a tilted oil-water contact has a unique formation mecha...The distribution characteristics of the oil-water contact are the basis for the reservoir exploration and development and reserves evaluation. The reservoir with a tilted oil-water contact has a unique formation mechanism, and the understanding of its distribution and formation mechanism will directly affect the evaluations for the reservoir type, well deployment, selection of well pattern and type, determination of test section, and reserves evaluation. Based on the analysis of reservoir characteristics, petrophysical properties and geological structure in 40 reservoirs worldwide with tilted oil-water contacts, the progress of the research on the formation mechanisms of titled oil-water contacts is summarized in terms of the hydrodynamic conditions, reservoir heterogeneity, neotectonic movement and oil-gas exploitation. According to the formation mechanism of tilted oil-water contacts and the needs of exploration research, different aspects of research methods are summarized and classified, such as the calculation of equipotential surfaces for oil and water in the formation, analysis of formation pressure and analysis of reservoir physical properties and so on. Based upon statistical analysis, it is suggested that the degree of the inclination of the oil-water contact be divided based on the dip of oil-water contact(DipTOWC). The tilted oil-water contact is divided into three categories: large dip(DipTOWC≥55 m/km), medium dip(4 m/km≤DipTOWC55 m/km), and small dip(DipTOWC4 m/km). The classification and evaluation method can be combined with structure amplitude and reservoir property. The formation mechanism of domestic and international reservoirs with tilted oil-water contacts are summarized in this paper, which have important significance in guiding the exploration and development of the oilfield with tilted oil-water contacts, reserves evaluation, and well deployment.展开更多
The characteristics of reservoir heterogeneity of the marine gravity flow tight sandstone from the Miocene Huangliu Formation under abnormally high pressure setting at LD10 area in Yinggehai Basin are studied,and the ...The characteristics of reservoir heterogeneity of the marine gravity flow tight sandstone from the Miocene Huangliu Formation under abnormally high pressure setting at LD10 area in Yinggehai Basin are studied,and the influencing factors on reservoir heterogeneity are discussed,based on modular formation dynamics test,thin sections,XRD analysis of clay minerals,scanning electron microscopy,measurement of pore throat image,porosity and permeability,and high pressure Hg injection,as well as the stimulation of burial thermal history.The aim is to elucidate characteristics of the heterogeneity and the evolution process of heterogeneity of the reservoir,and predict the favorable reservoirs distribution.(1)The heterogeneity of the reservoir is mainly controlled by the cement heterogeneity,pore throat heterogeneity,quality of the reservoir heterogeneity,and the diagenesis under an abnormally high pressure setting.(2)The differences in pore-throat structure caused by diagenetic evolution affected the intergranular material heterogeneity and the pore throat heterogeneity,and finally controlled the heterogeneity of reservoir quality.(3)Compared with the reservoir under normal pressure,abnormally high pressure restrains strength of the compaction and cementation and enhances the dissolution of the reservoir to some extent,and abnormally high pressure thus weakening the heterogeneity of the reservoir to a certain degree.The favorable reservoirs are mainly distributed in the gravity flow sand body under the strong overpressure zone in the middle and lower part of Huangliu Formation.展开更多
This article is aimed to discuss the impact of using two different kinds of surfactant in enhancing oil recovery in heterogeneous reservoirs. With the background of Jidong oilfield, Rui Feng surfactant which could rea...This article is aimed to discuss the impact of using two different kinds of surfactant in enhancing oil recovery in heterogeneous reservoirs. With the background of Jidong oilfield, Rui Feng surfactant which could reach ultra-low interracial tension and combination surfactant RZ-JD80 with strong emulsifying property are chosen to do oil displacement and profile control-oil displacement experiment in homogeneous core and heterogeneous core respectively. The experiment is aimed to study the effect of oil displacement by injecting surfactant individually and the effect after injecting different profile control agent slug before surfactant flooding in heterogeneous cores. The results suggest that injecting Rui Feng surfactant and RZ-JD80 individually could enhance the oil displacement efficiency about 15 percentage points for homogeneous core. For strongly heterogeneous core, it is low efficiency by using either of these two surfactants individually. However, if injected a very little profile control agent slug before surfactant flooding, both of these two kinds of surfactant could enhance the oil recovery by different degree, especially, polymer microsphere plugging^RZ-JD80 flooding composite technology is more adaptable to Gao-63 reservoir. This technology could increase the recovery by 18.52 percentage points aRer surfactant flooding.展开更多
One of the major challenges in reservoir characterization is to estimate the effective porosity and the permeability of the reservoir due to reservoir heterogeneity. Often, the vertical and the horizontal permeabiliti...One of the major challenges in reservoir characterization is to estimate the effective porosity and the permeability of the reservoir due to reservoir heterogeneity. Often, the vertical and the horizontal permeabilities are not considered separately in 3D geo-cellular models and in the reservoir simulations. Conventional reservoir modeling extrapolates all of the small-scale data to full-field scale data without considering the impact of the small-scale geological details, and therefore carries forward inherent errors into the reservoir predictions as a consequence of ignoring the reservoir heterogeneity. Most reservoirs are geologically complex and heterogeneous and that greatly influences reservoir performance. A case study is taken from a CHOPS (cold heavy oil production with sands) field. An innovative method of reservoir heterogeneity estimation has been introduced to illustrate the complex reservoir heterogeneity honouring all of the small-scale geological details in the 3D geological model. This detailed near-wellbore modeling through a synthetic core can provide the realistic quantitative volumetric assumption of the production prediction and improve the EOR (enhanced oil recovery) processes.展开更多
Tight glutenite reservoirs are widely developed in Bohai Bay Basin,East China.They are mostly huge thick and rely on hydraulic fracturing treatment for commercial exploitation.To investigate the propagation behavior o...Tight glutenite reservoirs are widely developed in Bohai Bay Basin,East China.They are mostly huge thick and rely on hydraulic fracturing treatment for commercial exploitation.To investigate the propagation behavior of hydraulic fractures in these glutenite reservoirs,the geological feature of reservoirs in Bohai Bay Basin is studied firstly,including the reservoir vertical distribution feature and the heterogeneous lithology.Then,hydraulic fracturing treatments in block Yan 222 are carried out and the fracturing processes are monitored by the microseismic system.Results show the hydraulic fractures generated in the reservoirs are mostly in X shape.The cause of X-shaped hydraulic fractures in this study is mainly ascribed to(I)the reservoir heterogeneity and(II)the stress shadow effect of two close hydraulic fractures propagating in the same orientation,which is confirmed by the following numerical simulation and related research in detail.This study can provide a reference for the research on the fracturing behavior of the deep thick glutenite reservoirs.展开更多
This study focuses on the heterogeneity of the middle Miocene syn-rift Belayim nullipore(reefal)marine sequences in the Gulf of Suez and its impacts on reservoir quality.The sequences consist of coralline algal reef l...This study focuses on the heterogeneity of the middle Miocene syn-rift Belayim nullipore(reefal)marine sequences in the Gulf of Suez and its impacts on reservoir quality.The sequences consist of coralline algal reef limestones with a highly complex dual-porosity system of primary and secondary porosities of widely varying percentages.To achieve a precise mathematical modeling of these reservoir sequences,a workflow protocol was applied to separate these sequences into a number of hydraulic flow units(HFUs)and reservoir rock types(RRTs).This has been achieved by conducting a conventional core analysis on the nullipore marine sequence.To illustrate the heterogeneity of the nullipore reservoir,the Dykstra-Parsons coefficient(V)has been estimated(V=0.91),indicating an extremely heterogeneous reservoir.A slight to high anisotropy(λ_(k))has been assigned for the studied nullipore sequences.A stratigraphic modified Lorenz plot(SMLP)was applied to define the optimum number of HFUs and barriers/baffles in each of the studied wells.Integrating the permeability-porosity,reservoir quality index-normalized porosity index(RQI-NPI)and the RQI-flow zone indicator(RQIFZI)plots,the discrete rock types(DRT)and the R35 techniques enable the discrimination of the reservoir sequences into 4 RRTs/HFUs.The RRT4 packstone samples are characterized by the best reservoir properties(moderate permeability anisotropy,with a good-to-fair reservoir quality index),whereas the RRT1 mudstone samples have the lowest flow and storage capacities,as well as the tightest reservoir quality.展开更多
Geostatistical data plays a vibrant role for surface-based reservoir modeling through outcrop analogues,which is used to understand three-dimensional(3D)variability of petrophysical properties.The main purpose of this...Geostatistical data plays a vibrant role for surface-based reservoir modeling through outcrop analogues,which is used to understand three-dimensional(3D)variability of petrophysical properties.The main purpose of this study is to improvise the surface-based 3D geo-modeling to demonstrate petrophysical characteristics and heterogeneities of Sandakan reservoirs,NW Borneo.We used point cloud data from Light Detection and Ranging(Li DAR)to build high-resolution virtual outcrop modeling(VOM)onto which we mapped 6 different lithofacies.Porosity and permeability of core plugs were measured to determine the average variance of petrophysical properties for each lithofacies.By utilizing the integration of VOMs analogues and petrophysical properties in Petrel^(TM),we demonstrated the distribution and associations of all lithofacies in pseudo wells that have inherent thin beds heterogeneities in 3D geocellular model.The results concluded that the heterogeneity of thin beds in lithofacies is dependent on porosity and permeability with input dataset.According to the final model,cross-bedding sandstone(CBS),hummocky cross-stratified sandstone(HCSS)and trough cross-bedding sandstone(TCBS)show good reservoir quality due to high porosity ranging from:25.6%to 20.4%and,19.3%-14.5%,and permeability ranging from:74.03 mD to 66.84 mD and,64.86 mD to 21.01 mD.In contrast,massive to weak laminae sandstone(MWLS)and bioturbated sandstone(BS)show fair to poor reservoir quality,caused baffling of surrounding mud sediments in the reservoir lithofacies.Results also revealed that Li DAR based VOM with petrophysical properties can significantly reduce the risk and minimize the cost of reservoir modeling in petroleum industry.展开更多
文摘Reservoir heterogeneities play a crucial role in governing reservoir performance and management.Traditionally,detailed and inter-well heterogeneity analyses are commonly performed by mapping seismic facies change in the seismic data,which is a time-intensive task.Many researchers have utilized a robust Grey-level co-occurrence matrix(GLCM)-based texture attributes to map reservoir heterogeneity.However,these attributes take seismic data as input and might not be sensitive to lateral lithology variation.To incorporate the lithology information,we have developed an innovative impedance-based texture approach using GLCM workflow by integrating 3D acoustic impedance volume(a rock propertybased attribute)obtained from a deep convolution network-based impedance inversion.Our proposed workflow is anticipated to be more sensitive toward mapping lateral changes than the conventional amplitude-based texture approach,wherein seismic data is used as input.To evaluate the improvement,we applied the proposed workflow to the full-stack 3D seismic data from the Poseidon field,NW-shelf,Australia.This study demonstrates that a better demarcation of reservoir gas sands with improved lateral continuity is achievable with the presented approach compared to the conventional approach.In addition,we assess the implication of multi-stage faulting on facies distribution for effective reservoir characterization.This study also suggests a well-bounded potential reservoir facies distribution along the parallel fault lines.Thus,the proposed approach provides an efficient strategy by integrating the impedance information with texture attributes to improve the inference on reservoir heterogeneity,which can serve as a promising tool for identifying potential reservoir zones for both production benefits and fluid storage.
基金supported by the National Science and Technology Major Projects of China for Oil and Gas(Projects Nos.2016ZX05010 and 2016ZX05058).
文摘Oil saturation is a critical parameter when designing oil field development plans.This study focuses on the change of oil saturation during water flooding.Particularly,a meter-level artificial model is used to conduct relevant experiments on the basis of similarity principles and taking into account the layer geological characteristics of the reservoir.The displacement experiment’s total recovery rate is 41.35%.The changes in the remaining oil saturation at a millimeter-scale are examined using medical spiral computer tomography principles.In all experimental stages,regions exists where the oil saturation decline is more than 10.0%.The shrinkage percentage is 20.70%in the horizontal well production stage.The oil saturation reduction in other parts is less than 10.0%,and there are regions where the oil saturation increases in the conventional water flooding stage.
基金sponsored by Major Science and Technology Special Project of CNPC(Grant No.2023ZZ07)。
文摘The Lower Cretaceous Manville Group of Upper Mc Murray Formation is one of the main bitumen reservoirs in Athabasca.In this study,the relationship between reservoirs heterogeneity and bitumen geochemical characteristics were analyzed through core and microscopic observation,lab analysis,petrophysics and logging data.Based on the sedimentology framework,the formation environment of high-quality oil sand reservoirs and their significance for development were discussed.The results indicate that four types lithofacies were recognized in the Upper Mc Murray Formation based on their depositional characteristics.Each lithofacies reservoirs has unique physical properties,and is subject to varying degrees of degradation,resulting in diversity of bitumen content and geochemical composition.The tidal bar(TB)or tidal channel(TC)facies reservoir have excellent physical properties,which are evaluated as gas or water intervals due to strong degradation.The reservoir of sand bar(SB)facies was evaluated as oil intervals,due to its poor physical properties and weak degradation.The reservoir of mixed flat(MF)facies is composed of sand intercalated with laminated shale,which is evaluated as poor oil intervals due to its poor connectivity.The shale content in oil sand reservoir is very important for the reservoir physical properties and bitumen degradation degree.In the context of regional biodegradation,oil sand reservoirs with good physical properties will suffer from strong degradation,while oil sand reservoirs with relatively poor physical properties are more conducive to the bitumen preservation.
基金funded by the National Natural Science Foundation of China(No.51974343)the Open Fund of Hubei Key Laboratory of Drilling and Production Engineering for Oil and Gas(Yangtze University)(No.YQZC202307)+2 种基金the Qingdao Postdoctoral Applied Research Project(No.qdyy20200084)the China Postdoctoral Science Foundation(No.2021M703588)a project of Science&Technology R&D Department of SINOPEC entitled Key Technology for Optimization Design of 3D Development of Shale oil in Dongying Sag(No.P23026).
文摘Currently,limited studies of immiscible water-alternating-CO_(2)(imWACO_(2))flooding focus on the impact of reservoir heterogeneity on reservoir development outcomes.Given this,using the heterogeneous reservoirs in the Gao 89-1 block as a case study,this study conducted slab core flooding experiments and numerical simulations to assess the impact of reservoir heterogeneity on imWACO_(2)flooding efficiency.It can be concluded that imWACO_(2)flooding can enhance the sweep volume and oil recovery compared to continuous CO_(2)flooding.As the permeability difference increases,the difference in the swept volume between zones/layers with relatively high and low permeability increases.To optimize the exploitation of reservoirs in the Gao 89-1 block,the optimal timing and CO_(2)injection rate for imWACO_(2)flooding are determined at water cut of 40%and 10000 m^(3)/d,respectively.A short injection-production semi-period,combined with multiple cycles of water and CO_(2)injection alternations,is beneficial for enhanced oil recovery from imWACO_(2)flooding.
基金This project is supported by the China National Key Basis Research Project (No: G1999022512)
文摘The processes of flooding—water flooding, polymer flooding and ternary combination flooding—were simulated respectively on a 2-D positive rhythm profile geological model by using the ASP numerical modeling software developed by RIPED (Yuan, et al. 1995). The recovery coefficient, remaining oil saturation, sweep efficiency and displacement efficiency were calculated and correlated layer by layer. The results show that the sweep efficiency and displacement efficiency work different effects on different layers in the severely heterogeneous reservoir. The study shows that the displacement efficiency and sweep efficiency play different roles in different layers for severely heterogeneous reservoirs. The displacement efficiency contributes mainly to the high permeability zones, the sweep efficiency to the low permeability zones, both of which contribute to the middle permeable zones. To improve the sweep efficiency in the low permeability zones is of significance for enhancing the whole recovery of the reservoir. It is an important path for improving the effectiveness of chemical flooding in the severely heterogeneous reservoirs to inject ternary combination slug after profile control.
基金Supported by the National Major Special Project of Oil and Gas During the 13th Five-Year Plan Period(NO.2016ZX05058-003-010)General Program of National Natural Science Foundation of China(NO.51574086)Postdoctoral Innovative Talent Support Program of China(NO.BX20190065)。
文摘Taking reservoir rocks and fluids of the Daqing,Dagang and Changqing oilfields as research objects,the EOR mechanisms and technical approach of polymer flooding were discussed.By comparing the displacement performances of ordinary polymer,glycerol,polymer in"sheet-net"structure and heterogeneous weak gel at the same viscosity and concentration,the relationship between the viscosity of polymer displacement agents and displacement performance was demonstrated,and the method of improving polymer flooding effect was worked out.The main mechanism of polymer flooding to increase oil recovery is the swept volume expansion of water injection due to polymer retention in porous media.The viscosity of polymer agents has no positive correlation with polymer flooding effect.Although polymer of"sheet-net"structure has strong capacity in increasing viscosity,it has poor compatibility with pore throat structure of reservoir rock,low injectivity and low shear resistance.Heterogeneous weak gel system has higher adsorption and capture capacity in porous media,which is easy to retain in porous media,and can effectively establish seepage resistance in high permeability layers(zones).Compared with polymer solutions with the same viscosity or concentration,it has stronger ability to expand swept volume.Long term injection of polymer flooding agents will inevitably lead to fluid entry profile reversal,and thus worsening of polymer flooding effect.Alternate injection of high retention and low or non-retention displacement agents can further improve the displacement effect of polymer flooding agents.
文摘During the displacement of water plugging with binary flooding in internally heterogeneous reservoirs,it is essential to understand the distributions of remaining oil as well as the oil displacement mechanisms at different stages.In this study,two types of internally heterogeneous systems,i.e.,vertical and horizontal wells are investigated experimentally through a microscopic approach.The results show that plugging agent types have a greater impact on oil recovery than well types,and foam injection can enhance oil recovery more effectively than gel injection.Additionally,the injection sequence of plugging agents significantly affects oil displacement efficiency.Injecting gel after foam is more beneficial.According to the present results,the main formation mechanisms of remaining oil in each displacement stage are influenced by:capillary force,viscous force,inertial force,shear force,microscopic fingering&channeling.
基金support from the National 973 Program (Grant No. 2006CB705808)PetroChina Young and Middle Aged People Innovation Fund (Grant No. 07E1016)+1 种基金PetroChina Science & Technology Innovation Fund (Grant No. 2008D-5006-02-09)Science & Technology Innovation Fund of Southwest Petroleum University (Grant No. 2007XJZ010)
文摘In view of the anisotropy,heterogeneity and stress-sensitive permeability in low permeability reservoirs,an analytical well test model was established by introducing the concept of permeability modulus.This model considered the permeability stress-sensitivity,wellbore storage effect,and the skin effect.The perturbation technique and Laplace transformation were used to solve the mathematical model analytically in Laplace space,and the bottom-hole pressure type curves were plotted and analyzed in real space by using the Stehfest numerical inversion.
基金The authors would like to acknowledge financial support provided by National Natural Science Foundation of China(No.52074338)We are also grateful to the support of the National Key R&D Program of China(No.2019YFA0708700)+1 种基金National Key Basic Research Program of China(20CX06071A)Bin Yuan would like to thank for the support of Shandong Mountain Tai Scholar Program.Chris Clarkson would like to acknowledge funding support from an NSERC Alliance grant(ALLRP 548652-19)for research related to the topic of this paper.
文摘Rate-transient analysis(RTA)has been widely applied to extract estimates of reservoir/hydraulic fracture properties.However,the majority of RTA techniques can lead to misdiagnosis of reservoir/fracture information when the reservoir exhibits reservoir heterogeneity and multiphase flow simultaneously.This work proposes a practical-yet-rigorous method to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,and improve the evaluation of reservoir/fracture properties.A new,general,semi-analytical model is proposed that explicitly accounts for multiphase flow,fractalbased reservoir heterogeneity,anomalous diffusion,and pressure-dependent fluid properties.This is achieved by introducing a new Boltzmann-type transformation,the exponent of which includes reservoir heterogeneity and anomalous diffusion.In order to decouple the effects of reservoir heterogeneity and multiphase flow during TLF,the modified Boltzmann variable allows the conversion of three partial differential equations(PDE's)(i.e.,oil,gas and water diffusion equations)into ordinary differential equations(ODE's)that are easily solved using the Runge-Kutta(RK)method.A modified time-power-law plot is also proposed to estimate the reservoir and fracture properties,recognizing that the classical square-root-of-time-plot is no longer valid when various reservoir complexities are exhibited simultaneously.Using the slope of the straight line on the modified time-power-law plot,the linear flow parameter can be estimated with more confidence.Moreover,because of the new Boltzmann-type transformation,reservoir and fracture properties can be derived more efficiently without the need for defining complex pseudo-variable transformations.Using the new semi-analytical model,the effects of multiphase flow,reservoir heterogeneity and anomalous diffusion on rate-decline behavior are evaluated.For the case of approximately constant flowing pressure,multiphase flow impacts initial oil rate,which is a function of oil relative permeability and well flowing pressure.However,multiphase flow has a minor effect on the oil production decline exponent.Reservoir heterogeneity/anomalous diffusion affect both the initial oil production rate and production decline exponent.The production decline exponent constant is a function of reservoir heterogeneity/anomalous diffusion only.The practical significance of this work is the advancement of RTA techniques to allow for more complex reservoir scenarios,leading to more accurate production forecasting and better-informed capital planning.
文摘This paper demonstrates the use of a commercial simulator as a tool with which to optimize the SAGD (steam-assisted gravity drainage) start-up phase process. The factors affecting the start-up phase are the prime targets. Among the key investigated factors are wellbore geometry effects, reservoir heterogeneity and circulation phase length. Each of the parameters was investigated via steam chamber development observation along the well pair length and the cross sections in the mid, toe and heel areas. In addition, the cumulative recovery in given time, steam-to-oil ratio and CDOR (calendar day oil rate) production data are used to backup the observations produced in the simulated model. Furthermore, an additional component developed during the research is a statistical modification of a layer cake model with which to create a heterogeneous reservoir to represent real reservoir conditions, based on Monte Carlo's simulation.
基金National Natural Science Foundation of China(51934005,U23B2089)Shaanxi Provincial Natural Science Basic Research Program Project(2024JC-YBQN-0554).
文摘The flow of fluid through the porous matrix of a reservoir rock applies a seepage force to the solid rock matrix.Although the seepage force exerted by fluid flow through the porous matrix of a reservoir rock has a notable influence on rock deformation and failure,its effect on hydraulic fracture(HF)propagation remains ambiguous.Therefore,in this study,we improved a traditional fluid–solid coupling method by incorporating the role of seepage force during the fracturing fluid seepage,using the discrete element method.First,we validated the simulation results of the improved method by comparing them with an analytical solution of the seepage force and published experimental results.Next,we conducted numerical simulations in both homogeneous and heterogeneous sandstone formations to investigate the influence of seepage force on HF propagation.Our results indicate that fluid viscosity has a greater impact on the magnitude and extent of seepage force compared to injection rate,and that lower viscosity and injection rate correspond to shorter hydraulic fracture lengths.Furthermore,seepage force influences the direction of HF propagation,causing HFs to deflect towards the side of the reservoir with weaker cementation and higher permeability.
文摘Terminal fans have formed the sedimentary system of the 2+3 sands of the upper second member, Shahejie formation in the west of the Pucheng Oilfield, Bohai Bay Basin, East China. Based on well logging data and physical properties of the reservoir beds, the 2+3 sands were divided into 16 sublayers. The heterogeneity of reservoir beds and distribution of interlayers and seal layers in the 2+3 sands were investigated. The intra-layer heterogeneity and inter-layer heterogeneity primarily belong to the severely heterogeneous classification. The spatial differentiation of sedimentary microfacies resulted in a change of reservoir bed heterogeneity, strong in the middle and southern parts, weak in the northern part. Spatial distribution of interlayers and seal layers is dominated by sedimentary microfacies, and they are thick in north-eastern and middle parts, thin in the south-western part.
文摘An analysis of drill cores and well logs shows that the main micro-facies of the third member sand bodies of the Qingshankou Formation in Qian'an are subaqueous distributary channel facies, sheet sand facies and subaqueous fan facies (olistostrome). Maps showing the distribution of these micro-facies together with inter-channel bay and prodelta mocro-facies are presented for different time-slices (lower, middle and upper parts of the Qingshankou Formation). These maps reveal the instability and change of sediment transport in the Baokang sedimentary system during the depositional period. Sediment transport was from the west in the early stage, from the south in the middle stage and from the northwest in the late stage. Values of thickness, porosity and permeability of the sand bodies in the third member of the Qingshankou Formation show that they have low to medium porosity and low permeability, and are characterized by serious reservoir heterogeneity. The joints between micro-facies and subaqueous fan micro-facies are characterized by the highest heterogeneity, the sheet sand and distal sand bar subfacies come next, and the heterogeneity of the subaqueous distributary channel sand bodies is relatively weak.
文摘The distribution characteristics of the oil-water contact are the basis for the reservoir exploration and development and reserves evaluation. The reservoir with a tilted oil-water contact has a unique formation mechanism, and the understanding of its distribution and formation mechanism will directly affect the evaluations for the reservoir type, well deployment, selection of well pattern and type, determination of test section, and reserves evaluation. Based on the analysis of reservoir characteristics, petrophysical properties and geological structure in 40 reservoirs worldwide with tilted oil-water contacts, the progress of the research on the formation mechanisms of titled oil-water contacts is summarized in terms of the hydrodynamic conditions, reservoir heterogeneity, neotectonic movement and oil-gas exploitation. According to the formation mechanism of tilted oil-water contacts and the needs of exploration research, different aspects of research methods are summarized and classified, such as the calculation of equipotential surfaces for oil and water in the formation, analysis of formation pressure and analysis of reservoir physical properties and so on. Based upon statistical analysis, it is suggested that the degree of the inclination of the oil-water contact be divided based on the dip of oil-water contact(DipTOWC). The tilted oil-water contact is divided into three categories: large dip(DipTOWC≥55 m/km), medium dip(4 m/km≤DipTOWC55 m/km), and small dip(DipTOWC4 m/km). The classification and evaluation method can be combined with structure amplitude and reservoir property. The formation mechanism of domestic and international reservoirs with tilted oil-water contacts are summarized in this paper, which have important significance in guiding the exploration and development of the oilfield with tilted oil-water contacts, reserves evaluation, and well deployment.
基金Supported by the Research on Exploration and Development Technology and New Exploration Field of High Temperature and Pressure Gas Reservoir in Western South China Sea(CNOOC-KJ135ZDXM38ZJ02ZJ)National Natural Science Foundation of China(41972129)National Science and Technology Key Project(2016ZX05024-005,2016ZX05026-003-005)。
文摘The characteristics of reservoir heterogeneity of the marine gravity flow tight sandstone from the Miocene Huangliu Formation under abnormally high pressure setting at LD10 area in Yinggehai Basin are studied,and the influencing factors on reservoir heterogeneity are discussed,based on modular formation dynamics test,thin sections,XRD analysis of clay minerals,scanning electron microscopy,measurement of pore throat image,porosity and permeability,and high pressure Hg injection,as well as the stimulation of burial thermal history.The aim is to elucidate characteristics of the heterogeneity and the evolution process of heterogeneity of the reservoir,and predict the favorable reservoirs distribution.(1)The heterogeneity of the reservoir is mainly controlled by the cement heterogeneity,pore throat heterogeneity,quality of the reservoir heterogeneity,and the diagenesis under an abnormally high pressure setting.(2)The differences in pore-throat structure caused by diagenetic evolution affected the intergranular material heterogeneity and the pore throat heterogeneity,and finally controlled the heterogeneity of reservoir quality.(3)Compared with the reservoir under normal pressure,abnormally high pressure restrains strength of the compaction and cementation and enhances the dissolution of the reservoir to some extent,and abnormally high pressure thus weakening the heterogeneity of the reservoir to a certain degree.The favorable reservoirs are mainly distributed in the gravity flow sand body under the strong overpressure zone in the middle and lower part of Huangliu Formation.
文摘This article is aimed to discuss the impact of using two different kinds of surfactant in enhancing oil recovery in heterogeneous reservoirs. With the background of Jidong oilfield, Rui Feng surfactant which could reach ultra-low interracial tension and combination surfactant RZ-JD80 with strong emulsifying property are chosen to do oil displacement and profile control-oil displacement experiment in homogeneous core and heterogeneous core respectively. The experiment is aimed to study the effect of oil displacement by injecting surfactant individually and the effect after injecting different profile control agent slug before surfactant flooding in heterogeneous cores. The results suggest that injecting Rui Feng surfactant and RZ-JD80 individually could enhance the oil displacement efficiency about 15 percentage points for homogeneous core. For strongly heterogeneous core, it is low efficiency by using either of these two surfactants individually. However, if injected a very little profile control agent slug before surfactant flooding, both of these two kinds of surfactant could enhance the oil recovery by different degree, especially, polymer microsphere plugging^RZ-JD80 flooding composite technology is more adaptable to Gao-63 reservoir. This technology could increase the recovery by 18.52 percentage points aRer surfactant flooding.
文摘One of the major challenges in reservoir characterization is to estimate the effective porosity and the permeability of the reservoir due to reservoir heterogeneity. Often, the vertical and the horizontal permeabilities are not considered separately in 3D geo-cellular models and in the reservoir simulations. Conventional reservoir modeling extrapolates all of the small-scale data to full-field scale data without considering the impact of the small-scale geological details, and therefore carries forward inherent errors into the reservoir predictions as a consequence of ignoring the reservoir heterogeneity. Most reservoirs are geologically complex and heterogeneous and that greatly influences reservoir performance. A case study is taken from a CHOPS (cold heavy oil production with sands) field. An innovative method of reservoir heterogeneity estimation has been introduced to illustrate the complex reservoir heterogeneity honouring all of the small-scale geological details in the 3D geological model. This detailed near-wellbore modeling through a synthetic core can provide the realistic quantitative volumetric assumption of the production prediction and improve the EOR (enhanced oil recovery) processes.
基金Projects(51879041,51774112,U1810203)supported by the National Natural Science Foundation of ChinaProject(2020M672224)supported by the China Postdoctoral Science FoundationProject(B2020-41)supported by the Doctoral Fund of Henan Polytechnic University,China。
文摘Tight glutenite reservoirs are widely developed in Bohai Bay Basin,East China.They are mostly huge thick and rely on hydraulic fracturing treatment for commercial exploitation.To investigate the propagation behavior of hydraulic fractures in these glutenite reservoirs,the geological feature of reservoirs in Bohai Bay Basin is studied firstly,including the reservoir vertical distribution feature and the heterogeneous lithology.Then,hydraulic fracturing treatments in block Yan 222 are carried out and the fracturing processes are monitored by the microseismic system.Results show the hydraulic fractures generated in the reservoirs are mostly in X shape.The cause of X-shaped hydraulic fractures in this study is mainly ascribed to(I)the reservoir heterogeneity and(II)the stress shadow effect of two close hydraulic fractures propagating in the same orientation,which is confirmed by the following numerical simulation and related research in detail.This study can provide a reference for the research on the fracturing behavior of the deep thick glutenite reservoirs.
基金the Researchers Supporting Project number(RSP-2020/92),King Saud University,Riyadh,Saudi Arabia。
文摘This study focuses on the heterogeneity of the middle Miocene syn-rift Belayim nullipore(reefal)marine sequences in the Gulf of Suez and its impacts on reservoir quality.The sequences consist of coralline algal reef limestones with a highly complex dual-porosity system of primary and secondary porosities of widely varying percentages.To achieve a precise mathematical modeling of these reservoir sequences,a workflow protocol was applied to separate these sequences into a number of hydraulic flow units(HFUs)and reservoir rock types(RRTs).This has been achieved by conducting a conventional core analysis on the nullipore marine sequence.To illustrate the heterogeneity of the nullipore reservoir,the Dykstra-Parsons coefficient(V)has been estimated(V=0.91),indicating an extremely heterogeneous reservoir.A slight to high anisotropy(λ_(k))has been assigned for the studied nullipore sequences.A stratigraphic modified Lorenz plot(SMLP)was applied to define the optimum number of HFUs and barriers/baffles in each of the studied wells.Integrating the permeability-porosity,reservoir quality index-normalized porosity index(RQI-NPI)and the RQI-flow zone indicator(RQIFZI)plots,the discrete rock types(DRT)and the R35 techniques enable the discrimination of the reservoir sequences into 4 RRTs/HFUs.The RRT4 packstone samples are characterized by the best reservoir properties(moderate permeability anisotropy,with a good-to-fair reservoir quality index),whereas the RRT1 mudstone samples have the lowest flow and storage capacities,as well as the tightest reservoir quality.
文摘Geostatistical data plays a vibrant role for surface-based reservoir modeling through outcrop analogues,which is used to understand three-dimensional(3D)variability of petrophysical properties.The main purpose of this study is to improvise the surface-based 3D geo-modeling to demonstrate petrophysical characteristics and heterogeneities of Sandakan reservoirs,NW Borneo.We used point cloud data from Light Detection and Ranging(Li DAR)to build high-resolution virtual outcrop modeling(VOM)onto which we mapped 6 different lithofacies.Porosity and permeability of core plugs were measured to determine the average variance of petrophysical properties for each lithofacies.By utilizing the integration of VOMs analogues and petrophysical properties in Petrel^(TM),we demonstrated the distribution and associations of all lithofacies in pseudo wells that have inherent thin beds heterogeneities in 3D geocellular model.The results concluded that the heterogeneity of thin beds in lithofacies is dependent on porosity and permeability with input dataset.According to the final model,cross-bedding sandstone(CBS),hummocky cross-stratified sandstone(HCSS)and trough cross-bedding sandstone(TCBS)show good reservoir quality due to high porosity ranging from:25.6%to 20.4%and,19.3%-14.5%,and permeability ranging from:74.03 mD to 66.84 mD and,64.86 mD to 21.01 mD.In contrast,massive to weak laminae sandstone(MWLS)and bioturbated sandstone(BS)show fair to poor reservoir quality,caused baffling of surrounding mud sediments in the reservoir lithofacies.Results also revealed that Li DAR based VOM with petrophysical properties can significantly reduce the risk and minimize the cost of reservoir modeling in petroleum industry.